白音華三號礦三維地質模型構建及凸臺邊坡穩定性分析

侯 偉

（國家電投集團內蒙古白音華煤電有限公司 露天礦，內蒙古 錫林郭勒 026209）

礦區三維地質建模不僅僅是二維地質建模向三維的延伸，三維地質建模技術突破了傳統地質構造二維表示的局限性，它能更直觀地描述露天礦山巖層及地質體的空間和結構位置關系，是實現露天礦“透明地質”的基礎。露天煤礦巖層多為沉積巖層，邊坡三維模型的構建對于礦山生產規劃及露天礦邊坡數值模擬等工作具有重要的意義。隨著人們對露天礦區地質精細程度的增加，研究人員也逐步意識到三維精細化地質模型在礦區生產的必要性，提出了一系列三維數值建模的措施。李春生等[1]以ArcGIS 軟件作為開發平臺，以礦區鉆孔數據庫為基礎，采用基于層面的約束不規則三角網（TIN）建立了露天礦地質模型；李佳君[2]基于三維礦業軟件DIMINE三維礦業軟件建立了白云鄂博礦區的地表DTM 模型和鉆孔數據庫模型，進一步對數據看進行品位整合后通過地質解譯建立礦體三維模型。三維地質模型的誤差分析一直是技術人員關注的重點，模型的準確度直接決定了地質模型用于指導生產的適用性。針對三維地質模型可能存在的誤差，劉光偉等[3]分析了三維建模過程中誤差來源，提出采用C－F 方法對三維地質模型各階段不確定性進行統一度量，以實現三維地質模型精度有效評價；朱良峰等[4]提出了三維地質結構模型精度評估、誤差檢測、動態修正的總體研究框架，為三維地質結構模型度與誤差修正的理論體和方法體系奠定了基礎。

露天礦邊坡穩定性研究經歷了由經驗分析、定性分析、單一評價理論到理論分析、定量分析、多種方法綜合評價的轉變。其中定性分析方法主要有工程地質分析法和圖解法[5]，定量分析方法主要有極限平衡法[6-8]和數值分析方法[9-12]等。其中數值分析方法已成為邊坡穩定性分析的一種重要的方法。隨著計算機技術的日益發展和成熟，有限單元法、離散元法、有限差分法、數值流形法等常用數值計算方法也在邊坡穩定性分析中得到了充分的發展，數值分析方法可以從較大范圍考慮邊坡介質的復雜性，全面地分析分析邊坡的應力應變狀況，有助于對邊坡變形和破壞機理的認識。

為此，以白音華三號露天礦已有的地質鉆孔信息為基礎，基于3Dmine 軟件和Rhino 軟件，結合無人機航拍圖，建立了能夠反映露天礦區巖土體空間分布的三維地質模型；進一步在白音華三號礦新布置了2 個鉆孔，經過鉆孔巖心圖與三維模型中同一位置巖性的分布對比，驗證了模型具有較高的準確性；在對地質信息充分了解后，建立了白音華三號礦凸臺邊坡三維數值模擬，基于三維數值模擬得到了邊坡主要的變形破壞區域。

1 工程概況

白音華三號礦隸屬于國家電投集團內蒙古能源有限公司，礦山位于內蒙古自治區錫林郭勒盟西烏珠穆沁旗白音華鎮管轄范圍內，白音華三號礦面積約35.1 km2，年產量1 500 萬t。由于礦權原因，在礦區的東北部形成了1 個三面解除約束的凸臺。白音華露天礦共有3 個煤組，其中凸臺下部壓有大量的優質煤炭，但是凸臺的北面出現了部分區域出現了明顯的變形，邊坡底部產生了0.5 m 的底鼓，凸臺的穩定性制約著煤炭的高效開采。

礦區所在的白音華煤田為大興安嶺南段西側的山間斷陷盆地，北東—北北東向長條狀展布。地塹式，盆地規模600 km2，長寬比近5：1，最大沉積厚度近1 400 m，含煤性好，富煤帶展布方向與盆地長軸方向一致。白音華三號露天礦位于煤田的中部，露天礦巖（土）體地層較為簡單，白音華三號礦典型巖土層分布如圖1。

圖1 白音華三號礦典型地質圖

巖土自上而下分布為第四系地層、第三系地層、泥巖層、煤層、泥砂巖互層等主要巖層。第四系地層主要是由腐殖土、砂土、細砂、粗砂、粉砂、亞砂土、砂礫和少量亞黏土、黃土等組成；第三系主要是由黏土或砂質黏土構成，第三系底部含有灰白－灰黃色砂礫層，主采煤層由2－1＃中煤層、2－1＃下煤層、3－1＃煤層、3－2＃煤層、3－3＃煤層5 層煤組成。地層厚度帶的總體方向為北東向，中部薄，西南、東北兩端厚，最大厚度帶偏向于井田西北一側，中部的厚度變薄帶為一個同沉積隆起。

2 白音華三號礦三維地質模型

利用無人機對整個礦區進行掃描得到礦區的高程線，基于高程線生成白音華三號礦DTM 模型，隨后在Rhino 中建立與表面大小相同的實體得到礦坑表面模型，在表面模型的基礎上建立礦區地質模型。

礦區地質模型主要是根據礦區已有的巖層結構剖面圖進行構建。Rhino 軟件建模的核心是NURBS曲面技術，NURBS 是指能夠用數學的方式精確的描述所有造型，即從簡單的2D 線到復雜的3D 有機自由曲面與實體。將巖性分界面與礦區表面模型進行整合，建立礦山巖層實體模型如。白音華露天礦三維地質體模型如圖2。

圖2 白音華露天礦三維地質體模型

對三維模型進行貼圖能夠更直觀反映礦坑真實情況，對于分析邊坡雷達點變形、裂隙分布等地質信息有重要意義。貼圖過程是在Rhino 軟件中完成，根據無人機航拍的高清礦區圖片，對所建立的白礦模型進行巖性賦予，選擇保存的無人機圖片并調整貼圖軸大小實現貼圖。白音華三號礦地表－地下精細化三維地質模型如圖3。

圖3 白音華三號礦地表－地下精細化三維地質模型

3 白音華三號礦三維地質模型準確性驗證

為了驗證模型的準確性，在白音華三號礦凸臺邊坡的西北部設計施工了2 個鉆孔，2 個鉆孔標高及鉆孔信息見表1。實際施工中根據巖石性質選擇鉆進方法，鉆頭規格選用?113 mm，巖心管選用?89 mm，回次進尺長度均未超過巖心管有效容納長度，鉆進松散、破碎的巖層時，根據實際情況采用低轉速、小泵量，保證取心質量。

表1 新增鉆孔信息表

基于新增的2 個鉆孔巖心巖性分布圖對所建模型準確度進行檢驗。首先根據鉆孔巖心圖標注出不同深度的巖性，然后將模型中對應深度的巖性與鉆孔中編錄的巖性進行對比。白音華三號礦模型準確性檢驗圖如圖4。圖4 中，四邊形的顏分別代表了鉆孔揭露巖性和模型表征的巖性，四變形的高度代表了不同巖性（砂巖、泥巖、煤層）。

圖4 白音華三號礦模型準確性檢驗圖

由對比可以看出：結果顯示模型能相對準確地反應巖土體的賦存狀態，尤其是在含煤的深度誤差較小，這對于煤層的開采具有重要的作用，通過對比，模型準確度較高。

4 白音華三號礦凸臺邊坡穩定性分析

主要通過FLAC3D軟件實現凸臺邊坡穩定性分析。其求解過程采用顯示有限差分格式，人機交互界面易于操作。針對應力場的控制微分方程進行求解，既不需要像極限平衡法一樣通過迭代求解，也不需要通過剛度矩陣，只需借助建立相應的本構關系，直接由應變求解應力。

建立的凸臺邊坡三維數值模型如圖5。

圖5 白音華三號礦凸臺邊坡三維計算模型

模型總高度為540 m，模型計算忽略巖體內部的構造應力僅在純重力作用下求解平衡，設置加載方向為豎直向下的重力加速為9.8 m/s2，對研究區域模型底面邊界固定x、y、z 位移，頂面設定為自由面，左側和右側邊界水平位移限制x 方向的水平位移；前、后兩邊界限制y 方向的水平位移。模型的z 方向豎直向上；x 軸正向指定為大地坐標系的“東”方向，即為模型的右側；正y 方向對應于大地坐標系的“北”方向，即指向邊坡的背面。

首先基于彈性本構模型求解凸臺邊坡的初始地應力，獲得初始地應力后基于摩爾－庫倫彈塑性本構模型求解邊坡的主要破壞區域。當模擬中的最大不平衡力與典型內力比值下限達到10－5時，認為系統達到平衡。

初始平衡狀態最大主應力如圖6，初始平衡狀態最小主應力云圖如圖7，初始平衡狀態最大剪應力云圖如圖8。

圖6 初始平衡狀態最大主應力云圖

圖7 初始平衡狀態最小主應力云圖

圖8 初始平衡狀態最大剪應力云圖

由圖6 可知：初始平衡狀態最小主應力隨深度增加而增加，等值線方向近乎與邊坡平行，在邊坡內部最大主應力方向與重力方向一致，在臨近坡表位移變為與坡面平行。

由圖7 可知：初始平衡狀態最小主應力隨埋深增加逐漸增大，同一水平靠近邊坡的位置應力小，遠離邊坡的位置應力大，最小主應力的方向在坡體內部一般為水平方向，而在臨近坡表調整為與坡面垂直。

由圖8 可知：坡面頂部剪應力較小，坡腳處出現了剪應力集中現象。

邊坡塑性區分布可以反映邊坡的破壞，凸臺邊坡塑性變形區域分布如圖9。

圖9 凸臺邊坡塑性變形云圖

由圖9 中可以看出邊坡局部產生了塑性變形，塑性變形主要集中在凸臺邊坡北面及西面邊坡坡底處，凸臺南面邊坡坡底的高程較北面坡底較高，凸臺南面邊坡相對于北面邊坡高度較小，因此在凸臺南面邊坡底部產生的破壞明顯少于北面邊坡底部；凸臺邊坡北面與西面的轉角處邊坡產生破壞較為嚴重，原因是凸臺北面與凸臺相向的邊坡為開采邊坡，開采邊坡的存在使得在凸臺北面與西面邊坡底部轉角處形成了應力集中，在該區域由于應力集中，底部巖體極易產生破壞，同時現場再該區域觀察到明顯的底鼓現象，也驗證了應力集中現象的存在；此外在凸臺邊坡頂部有部分塑性區產生，這些區域與現場變形較大的區域基本吻合，驗證了模擬的準確性，同時為邊坡監測及生產提供了參考。

4 結語

基于白音華三號礦的地下鉆孔信息及無人機航拍圖信息，利用3Dmine 軟件和Rhino 軟件建立了白音華三號露天煤礦三維地質模型。通過新增的鉆孔巖心數據驗證了模型的準確性，通過與新增鉆孔信息的對比認為模型能夠相對較真實地反映礦山的巖土體空間分布特征?；诎滓羧A三號礦的邊坡三維數值模型和Flac3D軟件研究得出邊坡破壞主要集中在邊坡坡底、凸臺邊坡坡頂和凸臺邊坡北側轉角處邊坡，模擬結果與現場變形區相吻合。