基于遍布節理模型的豎井工程破壞機理研究

翟淑花，高 謙，丁桂伶

（1.北京市地質研究所，北京 100120；2.北京科技大學金屬礦山高效開采與安全教育部重點實驗室，北京 100083）

基于遍布節理模型的豎井工程破壞機理研究

翟淑花1，高 謙2，丁桂伶1

（1.北京市地質研究所，北京 100120；2.北京科技大學金屬礦山高效開采與安全教育部重點實驗室，北京 100083）

為揭示金川二礦區14行風井的破壞機理，本文首先根據井筒工程地質條件，建立了集1#礦體和14行風井在內的三維數值模型；然后，采用遍布節理模型分析了采動對14行風井穩定性的影響，從井壁位移和應力方面探索了金川二礦區14行風井的破壞位置和破壞時間。結果表明，導致金川二礦區14行風井破壞的直接原因為采動引起豎井圍巖沿弱面的滑移，除此之外，復雜的工程地質條件和低強度的支護形式也是造成14行風井破壞的原因，為14行風井返修加固提供依據。

遍布節理模型；巖體移動；破壞機理；安全系數

0 引言

豎井是礦區建設中的咽喉工程，其穩定性不僅關乎工人的生命安全，還與礦山投資成敗息息相關，現有調查資料的統計結果表明，迄今為止我國已有十多個大中型礦山的豎井井筒受開采影響而產生變形和破壞，帶來了巨大的經濟損失。為此，許多科研工作者借助于理論分析、數值模擬、相似試驗模擬、原位測試以及統計理論等分析方法［1～16］對井筒的破壞機理進行了研究和探索，并取得了卓有成效的成績，但這些成果大多以煤礦厚沖積層中的豎井為研究對象，且提出了得到廣泛認可的厚沖積層沉降引起井壁豎直附加應力增加的研究理論。但針對金屬礦山礦井破壞機理的研究甚少，特別是采動影響下井筒變形破壞機理。

金川是我國最大的鎳生成基地和鉑族金屬提煉中心，是集采、選、冶和金屬加工于一體的特大型有色金屬聯合企業。金川礦山地質條件復雜，構造發育，礦巖破碎，巖體穩定性極差，在全國范圍內，屬于極難開采的礦山之一。金川二礦區采用的是充填法采礦，根據中國科學院地質與地球物理研究所為期5年半的GPS監測結果顯示［17］，金川二礦區地表巖體移動非常嚴重，地表沉降最大值已經達到889mm，并表現出極大的不均勻性和不對稱性，這給二礦區內的豎井工程的穩定性帶來了極大的威脅。2005年3月9日至22日金川礦區14行風井突然發生冒落事故，致使650m井筒被充填，距地表僅存170m井筒，造成工程不能使用，給金川二礦區帶來巨大的經濟損失。因此，開展14行風井的破壞機理研究既具有理論意義也具有經濟價值。

1 14行風井工程地質條件及破壞形式

14行風井位于14行勘探線的下盤，屬二礦區礦山改擴建工程，是二礦區專用回風井。承擔著二礦區1150m、1000m兩個主運輸水平，14個生產盤區的回風任務。14行風井工程穿過的工程圍巖有超基性巖、花崗巖、大理巖、混合巖等，局部巖體裂隙及小斷裂構造發育，在井筒底部1000m中段、距井筒下底約90m處有一斷層F16出露，該段井筒工程地質條件較差，詳見圖1。

二礦區14行風井破壞后，通過對上部尚未填實部分井筒進行調查，發現井筒破壞呈現以下幾種主要特征：

（1）井壁錯動開裂

二礦區14行風井井壁除了有環狀裂縫外，還具有豎向裂縫和X形交叉裂縫。裂縫主要出現在井壁西北及東南幫，以豎向裂縫為主。

圖1 14行風井井筒工程地質剖面示意圖

（2）片冒

片冒也是14行風井井壁圍巖變形破壞中常見的一種破壞形式，有的片幫比較輕微，僅在井壁出現貫通性的縱向裂縫，并有輕微的掉塊；有的則比較嚴重，混凝土支護破裂并出現離層，與圍巖之間形成空腔，根據混凝土破壞特征來看，它們的破壞機制為剪切破裂或拉張破裂。

（3）冒落物沖擊井壁破壞

在井下大約170m處，在南東40°左右位置井壁發現了一個窟窿，可能為對面井壁冒落物掉下所砸壞，而且該處堆積的冒落物東南高，距井口約171m，西北低，距井口約173m，因此也可以初步判定片冒主要發生在井壁的西北幫。井壁具體破壞情況如圖2所示。

圖2 二礦區14行風井井壁變形破壞照片

2 14行風井變形破壞仿真分析

由于巖體結構的復雜性，要建立完全反映巖體結構特征的模型是不現實的，因此對具體工程而言必須進行適當簡化，但同時必須認識到巖體強度由結構面強度控制。鑒于14行風井的實際地質情況的復雜性和三維建模的難度，本文論文采用FLAC3D自帶的遍布節理模型來模擬弱面對14行風井穩定性的影響，遍布節理模型實際上是莫爾-庫侖模型(Mohr-Cou lom b m odel)的擴展，即在莫爾-庫侖體中增加節理面，此節理面也服從莫爾-庫侖屈服準則。該模型同時考慮巖體和節理的物理力學屬性，破壞可能首先出現在巖體中或沿節理面，或二者同時破壞，其主要取決于巖體應力狀態、節理產狀、巖體及節理力學性質等。

2 .1 數值模型的建立

（1）遍布節理模型

遍布節理模型中巖體適用于莫爾-庫倫屈服準則，與節理面相關的屈服準則及其數值模擬計算的實現如下：

節理面方位和應力狀態可由局部坐標表示，廣義坐標和局部坐標下應力關系可由式（1）表示：

式中， θ為節理傾角

相應地，局部坐標下彈性應力和應變增量可由式（2）表示：

式中，K為體積模量；G為剪切模量

根據莫爾-庫倫準則，局部坐標下屈服包絡線AB可表示為 fs= 0；拉伸破壞包絡線BC可表示為 ft= 0；且有以下函數式存在：

式中，jc，jφ以及jtσ 分別表示節理面的凝聚力、內摩擦角及抗拉強度。當巖體應力處于穩定區域時，巖體呈現彈性狀態，不需要進行塑性修正，而進入屈服區域時，根據關聯流動法則，需進行修正。局部坐標與廣義坐標轉換后剪切破壞修正后的應力增量關系可表示為式（4）：

拉伸破壞(BC 段)修正后的應力增量關系可表示為式（5）：

（2）計算模型的建立

金川1#礦體是目前二礦區的主采礦體，礦體長1.6km，大寬度為200m，平均寬98m，延伸千余米。為研究采動對14行風井的影響，建立了包括1#礦體和14行風井在內的真三維模型。三維數值模型見圖3，模型尺寸為3600m×4000m×1700m，根據研究需要對14行附近圍巖網格進行加密，且采用映射型網格，單元總數140085，節點總數31388。模型x、y、z軸坐標范圍分別為（-1200，2400），（5000，9000），（0，1700）。計算中，采用實體單元模擬混凝土支護，混凝土強度等級為C30，巖體參數采用GPS監測數據反分析所得到的等效巖體參數，本構關系取為理想彈塑性，巖體及礦體采用莫爾-庫侖模型，鑒于14行風井所穿巖層節理分布較為密布，因此，井筒周圍一定范圍內的巖體采用遍布節理模型(ubiquitous-join tm odel)進行模擬（表1）。C30混凝土襯砌采用彈性模型。

表1 遍布節理模型參數

圖3 金川二礦區1號礦體及14行風井三維透視圖

（3）計算方法

計算分3步進行：首先計算在自重體積條件下使模型達到平衡狀態，獲得該條件下的初始應力場；在此基礎上將位移清零，按照二礦區實際采礦計劃（表2）進行數值模擬，根據開挖礦體前后應力場及位移場的變化，來分析地下開采對豎井變形、破壞的影響。

2.2 模擬結果分析

（1）變形分析

圖4為井筒沉降變形隨回采過程的變化曲線，揭示了井筒沉降變化規律。由于井深197m、標高為1520m附近的圍巖節理裂隙較為發育，因此，此段井壁沉降變形曲線在此發生突變，在靠近采空區標高1350m有個很緩的轉折。其次，在數值上，單中段回采時，不同回采步（1、2、3、4步）各回采步沉降差不是很大，雙中段開采（5、6步）較單中段開采時沉降數值變化較大，特別是第6步回采時，沿井壁軸線沉降曲線出現了大的變化。因此，可以初步判斷標高1520m處出現的沉降突變點很可能就是井壁最先破壞的位置，這與現場的井筒僅剩余170m情況相吻合。

（2）基于單元最小安全系數的井筒局部失穩分析

巖體強度破壞理論是判斷巖體失穩的重要判據，即在低約束壓力條件下，當巖體內某斜截面的剪應力超過破壞理論規定的滑動界限范圍時，巖體就發生剪切極限破壞。對于巖土介質，Moh r-Cou lom b服條件是工程界應用最為廣泛的屈服條件之一，其主應力表示形式以及安全系數sF可以表示為：

表2 二礦區開采順序表

其中，1σ為最大主應力，3σ為最小主應力，c為內聚力， φ為內摩擦角。

如果巖體中存在節理面，則對應于節理面破壞，其主應力表示形式以及安全系數sF為：

圖4 井壁沉降隨礦體開采過程的變化曲線

式中，1σ為最大主應力，3σ為最小主應力，jC是節理面粘結強度，jφ是節理面摩擦角， β為節理面傾角。

當 Fs＞1表示未破壞（屈服面內部）； Fs＜1表示已破壞（屈服面外部）； Fs=1表示處于臨界狀態（屈服面上部）。參照文獻［18］中的FLAC3D安全系數的FISH語言編寫方法，若巖體中任一點的應力狀態滿足屈服條件時，巖體在該點發生破壞。在進行強度破壞判斷時，采用抗拉強度準則和抗剪強度準則，拉應力當時，即認為巖體發生拉斷破壞；剪應力區：當單元的局部安全系數sF＜1時，即判斷該單元發生剪切破壞。采用該方法給出第6步回采時14行風井沿井筒標高（1000～1700m）最小安全系數變化曲線，見圖5。

圖5 井筒安全系數隨采礦步驟變化圖

由圖5可知，當回采至第6步時，井筒（標高1520～1480m）局部安全系數為0.984小于1，時間大約為2005年，這與風井破壞的時間相吻合，同時，在靠近采空區段，安全系數卻呈現逐漸增大的趨勢，這是由于開采引起卸荷作用對井筒起了積極作用，這也是14行風井的破壞沒有發生在1150m段而是發生在距離井口178m的原因。由此可見，14風井井筒的破壞完全是由于1520～1480m處的節理受采動影響發生了剪切破壞造成的，這一點與現場觀測情況相吻合。

3 結論

（1）建立了集14行風井和大型礦體為一體的三維幾何模型，采用遍布節理力學模型研究了井周密集節理對井筒穩定性的影響，揭示了采動影響下14行風井的破壞機理，結果表明，隨開采的不斷深入，14行風井井壁變形在1520～1480m段位移發生突變，當開采水平為1190～1170m和1130～1110m時，其局部安全系數降為0.984，井筒破壞，這與井筒實際的破壞位置和破壞時間相吻合，說明本次模擬的可靠性。

（2）通過分析可知，因14行風井局部節理密集，在采動的不斷擾動下，14行風井圍巖沿弱面（斷裂面、軟弱層面）發生了剪切破壞，是造成14行風井失穩垮塌的直接原因。除此之外，復雜的地質環境與條件，低強度的支護型式也是造成14行風井失穩垮塌的原因。

［1］汴正修．黃淮地區井筒破裂的地質因素［J］．煤炭科學技術，1993，(2)：56～60．

［2］馮開源．厚表土凍結施工井壁的裂露［J］．煤炭科學技術，1991，(4)：11～13．

［3］顧盂寒．立井井壁破壞機理與治理的探討［J］．煤炭科學技術，1991，(6)：44～47．

［4］楊 維，崔廣心．特殊地層條件下豎井井壁破壞機理及防治技術之一［J］．中國礦業大學學報，1996，25(4)：1～5．

［5］虞咸祥．大屯礦區李靖井壁斷裂原因分析［J］．煤炭科學技術，1991，(06)：41～44．

［6］崔廣心．特殊地層條件下豎井井壁破壞機理［J］．建井技術，1998，19(2)：29～32．

［7］崔廣心．特殊地層條件下豎井井壁破壞機理及防治技術［J］．建井技術，1998，19(1)：28～32．

［8］崔廣心，程錫祿．徐淮地區井壁破壞原因的初步研究［J］．煤炭科學技術，1991，(08)：46～50．

［9］呂恒林，崔廣心．深厚表土層中井壁結構破裂的力學機理［J］．中國礦業大學報，1999，28(6)：539～543．

［10］經來盛．表土沉降對井壁破裂的影響及防破裂措施的研究［J］．煤炭學報，2001，26(1)：49～53．

［11］黃志平．立井井壁破裂原因解釋［J］．煤炭科學技術，1991，(09)：50～53．

［12］張忠華．立井表土段井壁破裂原因分析［J］．煤炭科學技術，1992，(02)：43～44．

［13］陸孝軍．立井井壁破裂的原因分析［J］．煤炭科學技術，1994，22(12)：41～44．

［14］隧惠權，范學理．立井采動變形損壞特征及監測技術研究［J］．礦山測量，2002，(2)：35～37．

［15］李文秀，趙勝濤，梁旭黎等．魯中礦區地下開采對豎井井塔的影響分析［J］．巖石力學與工程學報，2006，25(1)：774～778．

［16］劉環宇，李 曉，王思敬等．厚沖積層立井井筒非采動破裂過程的三維耦合數值模擬［J］．工程地質學報，2004，12(S)：529～532．

［17］中國科學院地質與地球物理研究所，金川礦區地表巖移GPS監測、巖體移動變形規律與采動影響［R］，研究成果報告，2007．

［18］李樹忱，李術才，徐幫樹．隧道圍巖穩定分析的最小安全系數法［J］，巖土力學，28（03）：549～554．

Research on the Failing Mechanics of Shaft Engineering Based on the Ubiquitous
Joint M odel

ZHAI Shuhua1，GAO Qian2，DING Guiling1

(1. Beijing institution of Geology, Beijing 100120；2.State Key Laboratory of High-effi cient M ining and Safety of Metal M ines, Beijing University of Science and Technology, Beijing 100083)

In order to show the failing mechanics of No.14 ventilation shaft in Jinchuan m ine district, fi rstly, according to the geological conditions of shaft, three-dimension model was established to integrate No. 1 ore body and No.14 ventilation shaft. Then, the effect of m ining on the stability of No.14 ventilation shaft was evaluated numerically based on the ubiquitous joint model, the displacement and stress distribution of shaft were obtained to predict the failing location and time, which shows that the direct reason of failure of No.14 ventilation shaft is the sliding of surrounding rock mass induced by m ining. Besides, complex geological conditions and low-strength supporting type are also the factors causing the damage of shaft, which offer some basis for the repairing of shaft.

Ubiquitous joint model; Rock mass movement; Failing mechanics; Safety factor

TD265.3+3

A

1007-1903（2015）01-0061-06

翟淑花（1979- ），女，博士，高級工程師，主要從事巖土工程可靠性研究。