陳清運，楊從兵，王水平，張電吉

(1.武漢工程大學環(huán)境與城市建設學院,湖北 武漢 430074;2.磷資源開發(fā)利用教育部工程研究中心,湖北,武漢 430074;3.武鋼集團礦業(yè)有限責任公司程潮鐵礦,湖北 鄂州 436000;4.武鋼集團礦業(yè)有限責任公司金山店鐵礦,湖北 大冶 435116)

0 引 言

巖體是具有一定結構、賦存于一定地質環(huán)境中的地質體.巖體特性除受巖石力學性質、結構面力學性質和地質環(huán)境等因素影響外，還存在空間尺寸效應.巖體力學參數是進行巖土工程設計和穩(wěn)定性分析的基礎和關鍵，但是巖體力學參數不像巖石力學參數那樣容易獲取，確定受開采擾動的巖層力學參數就更難了.目前經常采用的方法有經驗折減系數法、大型原位試驗法、位移反分析法和巖體分類法等.

經驗折減系數法是在對巖體節(jié)理性狀定量描述的基礎上，通過弱化巖石強度求取巖體強度，雖然方法簡單，指標也易于獲得，但是精度常常達不到巖土工程設計和穩(wěn)定法分析的要求.大型巖體原位試驗，包括原位承壓板變形試驗和剪切試驗，用來研究巖體的變形特性、尺寸效應和軟弱結構面力學性質等，但是試驗費用高，次數有限.位移反分析法以工程開挖過程中圍巖變形觀測數據為基礎，通過反演分析獲得巖體力學參數和環(huán)境力學參數，因綜合了工程巖體各方面的信息，具有重要的參考價值.巖體分類法估計巖體力學參數是目前發(fā)展最快應用也最廣的方法，Bieniawski(1976，1989)提出了地質巖體質量指標(RMR)，并依此進行地質巖體分類，對不同的類別分別賦予不同的巖體力學參數值.Hoek和Brown(1980)，基于大量的試驗數據，提出了巖體經驗強度準則，建立了巖體的材料常數m、s、a與巖體質量指標之間的聯系，提出了一套估計巖體力學參數的方法.隨后，Hoek和Brown(1994、2004)提出了以地質強度指標(GSI)為基礎的廣義經驗強度準則，提出了一套巖體分類和力學參數估計的方法[1-2].

金山店鐵礦東區(qū)地表建筑物與構筑物較多，地下開采對地表影響范圍與程度是礦山關注的重要問題，正確確定巖層力學參數是合理評價地下開采巖層穩(wěn)定性的關鍵.

1 工程地質背景

金山店鐵礦張福山礦床以25勘探線為界將礦區(qū)分為東西兩大采區(qū)，在西采區(qū)又以16線為界，分為西Ⅰ采區(qū)和西Ⅱ采區(qū).西Ⅰ采區(qū)的工程地質特征簡述如下：

礦區(qū)地層主要由三迭系中—上統(tǒng)(T2-3)及侏羅系中—下統(tǒng)(J1-2W)地層組成.其中分布最廣的地層是三迭系中—上統(tǒng)薄圻群(T2-3pq)，為一套砂頁巖組合地層，分布于礦體上盤.對西Ⅰ采區(qū)影響較大的地質構造有F103和F104斷層.F103斷層走向北西，傾向南西，切割礦化帶及礦體，為一壓扭性斷層，斷層較緊閉，局部充填糜棱巖.F104斷層分布于礦體上盤與蒲圻群地層的接觸界面上，斷層走向與礦體走向近一致，傾向南，傾角65～85°，上陡下緩，具有壓扭性特征，斷裂帶中發(fā)育有碎裂巖和糜棱巖.巖漿巖在區(qū)內以石英閃長巖分布最廣，組成礦體下盤圍巖.

西Ⅰ采區(qū)礦體特征為：Ⅰ號礦體是西Ⅰ采區(qū)的主要礦體，東西長350 m，走向77～126°，傾向197～ 210°，傾角63～ 85°，礦體平均厚24.90 m.礦體出露地表，呈似層狀、脈狀，分枝復合現象明顯.Ⅲ號礦體與I礦體基本平行，走向長度為120 m，平均厚度15 m，兩礦體平均間距為35 m，Ⅲ號礦體位于I號礦體的上盤.

此前對19～23-1勘探線之間出露地表的礦體進行了露天開采，最終開采境界為+40 m水平，并逐年對露天坑進行了回填，礦山目前采用無底柱分段崩落法進行回采.

2 上覆巖層力學參數研究

2.1 巖體力學參數

a. 基于巖體質量指標(RMR)的巖體力學參數估計.Bieniawski地質力學分類法采用了5個分類參數：完整巖石材料的強度(或巖石點荷載強度指標)、巖石質量指標(RQD)、節(jié)理間距、節(jié)理狀態(tài)和地下水條件.因每個參數對巖體性狀影響不同，賦予了不同的權值，參數指標之和為巖體質量指標，根據巖體質量指標將巖體分五級.Serafim和Pereira (1983)，基于Hoek-Brown巖體經驗強度準則，根據巖體質量指標(RMR)，提出了估計巖體強度的計算方法，見表1.

表1 RMR法估計巖體力學參數計算體系

根據表1，編制相應程序進行計算，結果見表4.

b. 基于地質強度指標(GSI)的巖體力學參數估計.1995年，Hoek和Brown提出了地質強度指標GSI，巖體的地質強度指標與巖體結構、巖塊的嵌鎖狀態(tài)和巖體中不連續(xù)面的性狀有關.由于該法在巖體結構等級和不連續(xù)面性狀等級的劃分上過于粗糙，Sonmez和Ulusay對這兩個指標進行了細化和定量化，同時考慮不同的開挖方式對巖體地質強度指標GSI的影響，提出了巖體地質強度指標GSI的修正方法.

隨后，Hoek提出了適用范圍更廣的Hoek-Brown巖體經驗強度準則，并運用Balmer(1952)等以Mohr-Coulomb準則表示巖體力學極限平衡數學模型公式，提出了估計巖體強度 公式，見表2.根據表2，編制相應程序進行計算，結果見表4.

表2 GSI法估計巖體力學參數計算體系

c. 原位變形強度測試[3].現場承壓板變形試驗是為了確定角巖和石英閃長巖變形模量，這兩種巖石組成金山店鐵礦上下盤主要圍巖，其變形模量是影響巖層移動的主要因素.原位變形試驗選在-130 m水平的7#穿脈中進行的，上盤角頁巖中選一點，下盤石英閃長巖中選一點.試洞準備：上盤在距礦體40 m處的角巖中，下盤在距礦體50 m處的石英閃長巖中，于巷道壁上開挖出寬×高×深=1.5 m×1.4 m×2.0 m的試洞.儀器設備：加壓系統(tǒng)有液壓千斤頂(150T)、油壓泵、管路、壓力表；測量系統(tǒng)有測表支架、千分表、百分表、磁鐵表架、測量標點、基準工字鋼梁；傳力系統(tǒng)有剛性承壓板(Φ50 cm)，墊板、傳力柱、楔形塊.

依據所測得的應力—應變曲線及計算公式(1)，即可求得巖體的變形(彈性)模量，試驗結果見表3.

(1)

式(1)中：Em—巖體變形(彈性)模量(Pa)；

w0—巖體變形(m)；

p—承壓板上單位面積壓力(Pa)；

d—承壓板直徑(m)；

μ—泊松比.

表3 現場承壓板變形試驗結果

表4 巖體力學參數估計結果

由表4可知，采用地質力學分類法和經驗強度準則法進行巖體力學參數估計，其值有一定的誤差，取算術平均數，結果見表5.

表5 礦區(qū)巖體力學參數研究結果

2.2 結構面力學參數研究

結構面力學參數的研究是巖體力學參數研究的一個重要組成部分.經對張福山礦區(qū)幾種主要巖石的結構面進行了力學參數的測定，其結果如表6所示[4-5].

表6 巖體結構面力學參數

3 巖體力學參數工程驗證

3.1 計算模型的建立

計算范圍的確定：西Ⅰ采區(qū)開采18-1～23-1線之間的礦體，長200 m，包括Ⅰ、Ⅲ號兩礦體.主要研究西Ⅰ采區(qū)-141.5 m水平以上開采時巖層移動變形，但為今后研究方便，本次建立數值計算模型時從-500 m開始建立模型.根據類似工程條件下實測的地表邊界角：垂直礦體走向為57°、沿著礦體走向為80°.以該值為初值，圈定模型的計算范圍，并適當放大，以計算過程中應力重分布不影響到模型的邊界為原則，最后確定數值計算模型的范圍為：

X：3 334 600～3 335 600、Y：38 577 300～ 38 578 400、Z：-500～60，即長×寬×高=1 100 m×1 000 m×560 m.

工程地質和邊界條件：影響西Ⅰ采區(qū)的斷層有F103和F104，它們與礦體產狀基本一致或小角度斜穿礦體，對上下盤圍巖的影響比較小，因此在建模過程中可以不考慮斷層的影響；為了模擬露天坑中散體的力學性態(tài)，在建立三維數值計算模型時，上下盤以65°、端部以80°圈定露天坑散體，將之作為彈性體看待；破壞準則采用摩爾-庫侖模型(Mohr-Coulomb)；礦床在開采前要進行地下水的疏干，計算中可以不考慮地下水的影響；初始水平原巖應力按1.5倍自重應力賦初值.由于礦體不規(guī)則，用FLAC3D建立數值計算模型比較困難，因此先用ANSYS建立數值模型，然后導入FLAC中建立用于FLAC3D計算的模型，采用四面體單元對礦體和結構面進行擬合.在建立模型過程中主要依據礦區(qū)地形圖、剖面圖、各分段和階段的平面圖.建立的模型如圖1～2所示，單元數為88 625個，節(jié)點數為15 576個.

圖1 西Ⅰ采區(qū)礦體模型圖

圖2 西Ⅰ采區(qū)數值計算模型網格圖

3.2 計算結果及分析

為了反映礦區(qū)整體地表變形情況，在三條測線上都選擇了測點，同時兼顧到上下盤；另外，測線上位移變化較小的測點盡量不選，主要選擇位移變化比較活躍的測點；同時，不選位移變化奇異的測點.綜合分析，最后選擇A4、A5、B3、B12、C2等5點作為測點，主要是水平位移.

當西Ⅰ采區(qū)開采至-141.5 m水平時，在地表測點上，垂直礦體走向方向的水平位移計算值與實測值結果見表7，并將表7的結果繪成圖3.地表移動變形規(guī)律見圖4所示.

表7 計算值與實測值對比表

圖3 地表垂直礦體水平位移計算值與實測值對比圖

圖4 地表垂直礦體水平位移等值線圖

由圖3可見，兩條曲線吻合比較好，變化規(guī)律一致，測點實測值與計算值相對誤差在12％以下，對于采礦工程而言，這個精度是可以接受的.圖4反映的地表變形范圍與實際地表變形范圍一致.因此，巖體力學參數值可以作為東區(qū)地下開采地表變形的計算基礎數據.本結果為東區(qū)的地表變形分析提供重要的基礎數據，保證東區(qū)分析的正確性.

參考文獻：

[1]布雷迪 B H G, 布朗 E T.地下采礦巖石力學[M].馮樹仁，譯.北京：煤炭工業(yè)出版社,1990:318-321.

[2]Brady B H G,Brown E T. Rock Mechanics for underground mining[M].3 edition. LONDON： KLUWER ACADEMIC PUBLISHERS,2004:135-139.

[3]中國科學院武漢巖土力學研究所.金山店鐵礦張福山礦床地下采礦陷落角、錯動角研究[R].武漢：中國科學院武漢巖土力學研究所，1998.

[4]陳清運,何玉早.地下開采地表變形數值模擬研究[J].金屬礦山,2004(6):19-21.

[5]陳清運,徐臘明.數值模擬法在邊坡穩(wěn)定性分析中的應用[J].地下空間,1999,19(5):772-778.