固液耦合模式下采動誘發(fā)斷層界面應(yīng)力變化規(guī)律的模擬研究

張培森 ,楊 健 ,王明輝 ,王 浩

(1.山東科技大學(xué)資源與環(huán)境工程學(xué)院，山東 青島 266590；2.礦山災(zāi)害預(yù)防控制-省部共建國家重點實驗室培育基地，山東 青島 266590；3.山西天地王坡煤業(yè)有限公司，山西 晉城 048021)

由于煤炭為一次性能源，隨著煤炭開采深度及強度的不斷加大，淺部的資源已頻臨枯竭，煤炭開采以15～30m/a的速度向深部發(fā)展，當(dāng)前采深大于600m的已占30%，超過千米的深井達(dá)39處，預(yù)計到2030年我國很多煤礦的開采深度將達(dá)到1000～1500m。深部煤層的賦存條件也越趨復(fù)雜，開采難度大大提高，煤礦事故，特別是斷層突水事故是當(dāng)前深部開采所面對的主要難題之一[1-2]。

采動誘發(fā)斷層活化是采礦學(xué)科中典型的非線性力學(xué)問題，是指含煤地層中的斷層在采動影響下產(chǎn)生滑移的現(xiàn)象[3]。斷層活化使得斷層帶以及兩盤靠近斷層的部位產(chǎn)生新的裂隙，原有的裂隙進一步擴大，最后形成突水通道，為工作面突水提供了重要的條件[3]。長期以來，國內(nèi)不少學(xué)者對斷層突水進行了大量的研究[1-10]，也取得了很多有價值的成果，并進行了大量有關(guān)斷層突水的相似材料模擬試驗、理論分析及數(shù)值模擬研究，而有關(guān)固液耦合模式下采動誘發(fā)斷層界面應(yīng)力變化對斷層活化突水的影響研究卻不多。本文以五溝煤礦F16斷層地質(zhì)條件為背景，采用大型巖土類數(shù)值模擬軟件FLAC3D的固液耦合模式對采動誘發(fā)斷層界面應(yīng)力變化規(guī)律進行數(shù)值模擬研究，為斷層安全防水煤巖柱的留設(shè)及斷層突水防治提供理論參考。

1 模型建立及計算方案

1.1 模型構(gòu)建及計算參數(shù)選取

五溝煤礦F16斷層落差為70m，開采煤層為10煤，平均厚度為4m。計算模型選取工作面的推進方向為X軸，工作面的傾斜方向為Y軸，模型的豎直方向為Z軸?？紤]到邊界效應(yīng)和計算速度因素，最終確定計算模型大小為600m×10m×200m。根據(jù)研究目的，對斷層、煤層及煤層直接頂和直接底網(wǎng)格進行加密處理，具體網(wǎng)格劃分如圖1所示。模型四周邊界均施加水平位移約束，底邊界均施加水平位移及垂直位移約束，上部邊界為自由面，老頂上方載荷與上覆巖層的重力(∑γh)有關(guān)。

經(jīng)過大量巖石力學(xué)試驗以及工程類比可知，Mohr-Coulomb準(zhǔn)則能很好的反映出巖石破壞的特性，故本次計算選取MOHR模型及FL_ISO模型進行分析。

計算采用莫爾-庫侖(Mohr-Coulomb)屈服準(zhǔn)則判斷巖體的破壞。

式中：σ1、σ3分別是最大和最小主應(yīng)力;φ、c分別是摩擦角和粘結(jié)力。當(dāng)fs>0時，材料將發(fā)生剪切破壞。在通常應(yīng)力狀態(tài)下，巖體的抗拉強度很低，因此，可根據(jù)抗拉強度準(zhǔn)則(σ3≥σt)判斷巖體是否產(chǎn)生拉斷破壞。

根據(jù)現(xiàn)場的地質(zhì)調(diào)查和相關(guān)試驗研究所提供的巖石力學(xué)試驗數(shù)據(jù)，在考慮巖石尺度效應(yīng)的基礎(chǔ)上，模擬計算采用的主要巖層的物理力學(xué)參數(shù)見表1。

圖1 模型網(wǎng)格劃分

表1 主要巖層的物理力學(xué)參數(shù)

1.2 計算方案

針對五溝煤礦F16斷層的地質(zhì)條件，本文擬采用三維有限差分?jǐn)?shù)值軟件FLAC3D的固液耦合模式對采動誘發(fā)斷層界面應(yīng)力變化規(guī)律進行數(shù)值模擬研究。斷層界面應(yīng)力變化規(guī)律與多種因素有關(guān)，在分析現(xiàn)有研究資料的基礎(chǔ)上，本文重點分析斷層傾角、斷層破碎帶寬度及承壓水水壓不同時在斷層上盤煤層開采過程中斷層界面應(yīng)力變化規(guī)律。

本次計算方案共分析11個計算工況，分為三類：第一類為斷層破碎帶寬度為4m，含水層水壓為3MPa保持不變的前提下，斷層傾角分別為30°、40°、50°、60°、70°；第兩類為斷層傾角為70°，含水層水壓為3MPa保持不變的前提下，斷層破碎帶寬度分別為4m、6m、8m；第三類為斷層破碎帶寬度為6m，斷層傾角為70°保持不變的前提下，含水層水壓分別為2MPa、3MPa、4MPa。

2 計算結(jié)果分析

開采煤層位于斷層上盤，考慮到邊界效應(yīng)，工作面從x=100m處開挖推進。從左向右，向斷層推進，直到距離斷層50m時回采工作結(jié)束。在斷層上下盤同一水平布置布設(shè)一對監(jiān)測點，通過分析這一對監(jiān)測點應(yīng)力的變化來研究斷層界面的應(yīng)力變化規(guī)律。

2.1 斷層傾角不同對斷層界面應(yīng)力場變化規(guī)律的影響

圖2～6分別描述了斷層傾角不同時，在煤層開采過程中，斷層界面上同一標(biāo)高處不同監(jiān)測點上的正應(yīng)力及剪應(yīng)力隨計算時步增加而變化的情況。

圖2 斷層傾角為30°時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖3 斷層傾角為40°時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖4 斷層傾角為50°時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖5 斷層傾角為60°時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖6 斷層傾角為70°時煤層開采過程中斷層界面上的應(yīng)力變化曲線

從圖2～6中可以看出，在距離斷煤交線50m時，斷層的上盤界面所受的正應(yīng)力大于下盤斷層界面所受的正應(yīng)力，斷層上盤界面所受的剪應(yīng)力大于斷層下盤界面所受的剪應(yīng)力；隨著斷層傾角的增加，在工作面推進至斷層50m時，斷層上、下盤界面所受的正應(yīng)力逐漸減小，而斷層上、下盤所受的剪應(yīng)力一開始隨著斷層傾角的增加而增加，在斷層傾角超過50°之后又隨著斷層傾角的增加而減小。

2.2 斷層破碎帶寬度不同對斷層界面應(yīng)力場變化規(guī)律的影響

圖7～9分別描述了斷層破碎帶寬度不同時，在煤層開采過程中，斷層界面上同一標(biāo)高處不同監(jiān)測點上的正應(yīng)力及剪應(yīng)力隨計算時步增加而變化的情況。

圖7 斷層破碎帶寬度為4m時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖8 斷層破碎帶寬度為6m時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖9 斷層破碎帶寬度為8m時煤層開采過程中斷層界面上的應(yīng)力變化曲線

從圖7～9中可以看出，在工作面距離斷層50m時，斷層上盤界面所受的正應(yīng)力大于斷層下盤界面所受的正應(yīng)力，斷層上盤所受的剪應(yīng)力大于斷層下盤所受的剪應(yīng)力；隨著斷層破碎帶寬度的增大，在斷層破碎帶寬度從4m增加到8m時，斷層上、下盤界面所受的正應(yīng)力越來越大，而斷層上、下界面所受的剪應(yīng)力逐漸減小。

2.3 含水層水壓不同對斷層界面應(yīng)力場變化規(guī)律的影響

圖10～12分別描述了含水層水壓不同時，在煤層開采過程中，斷層界面上同一標(biāo)高處不同監(jiān)測點上的正應(yīng)力及剪應(yīng)力隨計算時步增加而變化的情況。

圖10 含水層水壓為2MPa時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖11 含水層水壓為3MPa時煤層開采過程中斷層界面上的應(yīng)力變化曲線

圖12 含水層水壓為4MPa時煤層開采過程中斷層界面上的應(yīng)力變化曲線

從圖10～12中可以看出，在距離斷煤交線50m時，斷層的上盤界面所受的正應(yīng)力大于斷層下盤界面所受的正應(yīng)力，斷層上盤界面所受的剪應(yīng)力大于斷層下盤界面所受的剪應(yīng)力；隨著含水層水壓的增加，在工作面推進至斷層50m時，斷層上、下盤界面所受的正應(yīng)力逐漸增大，斷層上、下盤所受的剪應(yīng)力也隨之增加。

3 結(jié)論

1)在相同工況條件下，斷層的上盤斷層界面所受的正應(yīng)力大于下盤斷層界面所受的正應(yīng)力，斷層上盤界面所受的剪應(yīng)力大于斷層下盤界面所受的剪應(yīng)力；斷層傾角越大，在工作面推進相同距離時，斷層上、下盤界面所受的正應(yīng)力越趨減小，而斷層上、下盤所受的剪應(yīng)力一開始隨著斷層傾角的增加而增加，在斷層傾角超過50°之后又隨著斷層傾角的增加而減小。

2)斷層破碎帶寬度越大，斷層上、下盤界面所受的正應(yīng)力越大，而斷層上、下界面所受的剪應(yīng)力逐漸減小。

3)隨著含水層水壓的增加，在工作面推進至斷層50m時，斷層上、下盤界面所受的正應(yīng)力逐漸增大，斷層上、下盤所受的剪應(yīng)力也隨之增加；水壓越大，斷層突水危險性就越大。

[1] 黃存捍.采動斷層突水機理研究[D].長沙:中南大學(xué), 2010.

[2] 何滿潮, 謝和平，彭蘇萍，等.深部開采巖體力學(xué)研究[J].巖石力學(xué)與工程學(xué)報, 2005, 24(16):2803-2813.

[3] 于廣明, 謝和平，楊倫，等.采動斷層活化分形界面效應(yīng)的數(shù)值模擬研究[J].煤炭學(xué)報, 1998, 23(4):396-400.

[4] 卜萬奎.采場底板斷層活化及突水力學(xué)機理研究[D].徐州:中國礦業(yè)大學(xué), 2009.

[5] 黎良杰,錢鳴高,李樹剛.斷層突水機理分析[J].煤炭學(xué)報,1996,21(2):119-123.

[6] 卜萬奎,茅獻彪.斷層傾角對斷層活化及底板突水的影響研究[J].巖石力學(xué)與工程學(xué)報,2009,28(2):386-394.

[7] Ran J Q.Shear slig ing failure of the jointed roof in laminated rock mass[J].Rock Meth Rock Engng,1994,27(4):235-251.

[8] 李凱,茅獻彪,陳龍，等.采動對承壓底板斷層活化及突水危險性的影響分析[J].力學(xué)季刊, 2011,32(2):261-268.

[9] 張金才, 張玉卓, 劉天泉.巖體滲流與煤層底板突水[M].北京:地質(zhì)出版社, 1997.

[10] 盧興利, 尤春安, 孫鋒，等.斷層保護煤柱合理留設(shè)的數(shù)值模擬分析[J].巖土力學(xué),2006,27(zk):239-242.