任仰輝

（煤炭科學(xué)研究總院有限公司，北京 100013）

煤炭作為我國主體能源，在我國能源構(gòu)成中長期處于基礎(chǔ)地位。長期以來煤炭資源的高強(qiáng)度開采使得地層淺部容易開采的煤層已所剩無幾，不得不轉(zhuǎn)而開采深部煤炭，而深部開采往往受到底板承壓水的威脅。深部開采初期，由于對承壓水的賦存特點(diǎn)及對煤炭開采的影響研究較少，缺乏相應(yīng)技術(shù)支撐，造成承壓水上開采工作面底板突水事故時(shí)有發(fā)生?！笆濉逼陂g發(fā)生重大水害事故17 起（死亡229 人），占同期全國煤礦重特大事故23.6%。為實(shí)現(xiàn)承壓水上煤炭資源的安全回收，國內(nèi)學(xué)者進(jìn)行了大量研究：文獻(xiàn)[1]提出了承壓水上采煤的框架結(jié)構(gòu)及該學(xué)科的主要研究內(nèi)容和方法；文獻(xiàn)[2]系統(tǒng)梳理了我國煤礦水害防治技術(shù)新進(jìn)展、新方法，并對未來水害防治技術(shù)體系發(fā)展進(jìn)行了展望；文獻(xiàn)[3-6]利用相似材料模擬帶壓開采煤層底板破壞與裂隙發(fā)育規(guī)律；文獻(xiàn)[7]以微震監(jiān)測為基礎(chǔ)，以微震能量密度及連通度表征底板巖層采動(dòng)裂隙的導(dǎo)通性，為推動(dòng)應(yīng)用微震監(jiān)測技術(shù)防治底板突水提供了新路徑；文獻(xiàn)[8-9]建立了底板突水危險(xiǎn)性評價(jià)方法，并提出了底板奧灰水防治關(guān)鍵技術(shù)；文獻(xiàn)[10]構(gòu)建了周期來壓時(shí)工作面底板力學(xué)模型，并結(jié)合鉆孔注水和直流電法等方法實(shí)測了承壓水上開采底板破壞深度；文獻(xiàn)[10-13]采用數(shù)值模擬方法研究了底板破壞和突水通道形成過程，并建立力學(xué)模型分析了突水危險(xiǎn)性，研究了隱伏小斷層對突水的影響；文獻(xiàn)[14-17]系統(tǒng)研究了承壓水上膏體（似膏體）充填開采技術(shù)體系、底板采動(dòng)效應(yīng)等；文獻(xiàn)[18]通過建立力學(xué)分析模型，分析了充填開采采空區(qū)垮落帶、斷裂帶巖體滲透性特征，研究了隔水關(guān)鍵層對水害防治的作用；文獻(xiàn)[19-20]研究了固體充填材料影響巖層控制效果的機(jī)理，建立了固體充填開采充填材料選擇設(shè)計(jì)方法，設(shè)計(jì)了煤礦井下煤矸分選與充填開采一體化系統(tǒng)。以上研究成果為受承壓水威脅采場的底板控制提供了重要參考，但深部開采面臨的高應(yīng)力、高承壓水等特殊條件，仍然對深部煤炭資源的安全高效經(jīng)濟(jì)開采帶來嚴(yán)峻挑戰(zhàn)。同時(shí)，充填開采在減少煤層底板變形方面效果明顯，且隨著技術(shù)進(jìn)步已得到越來越多的應(yīng)用。因此，進(jìn)一步研究深部承壓水上煤層采用充填開采時(shí)的底板破壞效應(yīng)具有重要現(xiàn)實(shí)意義和應(yīng)用價(jià)值。為此，采用理論計(jì)算、數(shù)值分析和現(xiàn)場觀測相結(jié)合的方法對工作面承壓水上充填開采底板破壞特征進(jìn)行了綜合研究，研究結(jié)果有助于摸清充填開采對底板水害的控制作用，為保障承壓水上的安全開采提供參考依據(jù)。

1 工作面概況

某礦A 工作面地面標(biāo)高+74～+75 m，工作面標(biāo)高-530～-600 m，主采七號煤，煤層產(chǎn)狀大致為走向北西，傾向北東，傾角5°～11°，平均7°，七號煤層賦存較為穩(wěn)定，煤厚平均1.5 m。工作面直接頂為泥巖，黏結(jié)性較好，可隨采隨垮落，基本頂為細(xì)砂巖，厚6 m，強(qiáng)度中等，屬易管理頂板。7號煤層下距底板四灰含水層22 m 左右，四灰含水層水量不大。下距奧陶系石灰?guī)r含水層40 m 左右，奧灰含水層水量豐富。在四灰和奧灰含水層之間，賦存有厚度約17 m 的以泥巖為主的巖層，隔水性較好。

2 充填開采底板破壞深度分析

充填開采是指煤層采出后，利用回采和掘進(jìn)矸石、膏體充填料、高水材料或高濃度膠結(jié)材料等充填體充填采空區(qū)，使充填體和周邊煤柱共同支撐上覆基巖載荷。由于充填工藝、充填材料和工作面管理等問題，充填體與頂板巖層之間一般存在一定的空頂高度?？枕敻叨鹊拇嬖谑钩涮铙w不能及時(shí)發(fā)揮承載能力，因此工作面推過之后的一定時(shí)期內(nèi)，上方基巖載荷仍由周邊煤柱和工作面前方煤壁支撐，導(dǎo)致支承壓力集中現(xiàn)象的產(chǎn)生。

對于采空區(qū)底板中的任意一點(diǎn)M，其受力一般可認(rèn)為由4 個(gè)部分構(gòu)成，包括M之上底板巖體的自重、充填體之上的頂板載荷、充填體的重力作用和工作面煤壁前面一定距離頂板巖層的支承壓力作用。

根據(jù)圣維南原理，作用在物體的一小部分邊界上的面力，變換為分布不同但作用等效的面力，那么近處的應(yīng)力分布將有顯著的改變，但遠(yuǎn)處所受的影響可以忽略不計(jì)[21]，因此對M點(diǎn)應(yīng)力狀態(tài)造成影響的非線性分布的支承壓力可簡化為分段載荷，充填開采底板應(yīng)力簡化模型如圖1。

圖1 充填開采底板應(yīng)力簡化模型Fig.1 Simplified model of floor stress in back filling mining

圖1 中，q為煤壁前方支承壓力的等效應(yīng)力，其值等于ρgH（n+1）/2；ρ為煤巖體密度；H為煤層埋藏深度，n為支承壓力集中系數(shù)，一般取2～3。根據(jù)相關(guān)力學(xué)原理可知q的作用寬度L=2xa。由土力學(xué)知識可知，M點(diǎn)應(yīng)力狀態(tài)是支承壓力同頂板載荷、充填體重力、底板載荷共同決定的，對各因素對M點(diǎn)應(yīng)力狀態(tài)的影響進(jìn)行疊加可得底板巖層中任意一點(diǎn)M的主應(yīng)力為：

式中：σ1為最大主應(yīng)力；σ3為最小主應(yīng)力；ρc為充填體密度；g為重力加速度；Hc為充填體穩(wěn)定后高度；α為最大主應(yīng)力與x軸的夾角；z為M點(diǎn)到煤層底版的深度。

代入Mohr-Coulomb 破壞準(zhǔn)則，即σ1-Kσ3=Rc，整理后可得：

式中：Rc為巖體單軸抗壓強(qiáng)度；K=；φ0為底板巖體內(nèi)摩擦角。

對式（2）進(jìn)行求導(dǎo)，即可得到底板破壞深度達(dá)到極大值時(shí)的α角：

由式（3）可得當(dāng)角度α滿足式（4）時(shí)，底板破壞深度達(dá)到極大值。

280 Establishment and improvement of abdominal heterotopic heart transplantation model in mice

將式（4）及q=ρgH（n+1）/2 代入式（2）可求得充填開采工作面底板巖層最大破壞深度Hm為：

由式（5）可以看出，采用充填開采工藝時(shí)，煤層底板最大破壞深度與煤層埋深、煤壁前方支承壓力成正比，與底板巖體強(qiáng)度、充填體密度及充填高度負(fù)相關(guān)，將A 工作面及煤巖體物理力學(xué)參數(shù)代入式（5）可得該工作面底板的最大破壞深度為20 m，由于在承壓水導(dǎo)升帶與采動(dòng)影響破壞帶之間存在厚17 m 的隔水關(guān)鍵層，由理論計(jì)算可以看出，采用充填采煤法控制底板變形時(shí)，A 工作面沒有突水危險(xiǎn)，可以實(shí)現(xiàn)七號煤層的安全回采。

3 承壓水上充填開采的數(shù)值模擬

根據(jù)A 工作面頂?shù)装鍘r層綜合柱狀圖建立的數(shù)值模型如圖2，模型高157 m，走向長400 m，模型四周和底部固定位移約束，頂部施加8 MPa的應(yīng)力邊界條件，底板奧灰?guī)r中施加3.5 MPa 的孔隙水壓，煤巖體破壞準(zhǔn)則選取庫倫-摩爾準(zhǔn)則，采用充填采煤工藝，充填材料選取應(yīng)變硬化模型，并假設(shè)充填工作與工作面開采并行。煤巖及充填體物理力學(xué)參數(shù)見表1。

表1 煤巖層物理力學(xué)參數(shù)表Table 1 Physical and mechanical parameters of coal and rock

圖2 數(shù)值模型Fig.2 Numerical calculation model

3.1 推進(jìn)不同距離頂?shù)装逵行?yīng)力分布

充填開采情況下，工作面推進(jìn)至不同距離時(shí)，工作面頂?shù)装逵行?yīng)力分布如圖3。

圖3 推進(jìn)至不同距離時(shí)圍巖有效應(yīng)力分布Fig.3 Effective stress distribution of surrounding rock at different distances

由于充填體的承載作用，煤壁前方支承壓力集中系數(shù)很小，工作面推進(jìn)50 m 時(shí)約為1.22，推進(jìn)至150 m 時(shí)應(yīng)力集中系數(shù)約為1.51，增幅較小。當(dāng)推進(jìn)距離不大時(shí)，由于頂板下沉量較小，頂板同充填體之間的接觸耦合作用不強(qiáng)，充填體發(fā)揮的承載能力有限，如圖3（a），此時(shí)同垮落法管理頂板類似，工作面頂?shù)装宄霈F(xiàn)較大范圍的卸壓區(qū)；而當(dāng)工作面推進(jìn)至150 m 時(shí)，如圖3（b），由于頂板已多次垮落，覆巖充分下沉，充填體的承載能力得以充分發(fā)揮，此時(shí)采空區(qū)中部垂直應(yīng)力已恢復(fù)至原巖應(yīng)力水平。

另外，由于承壓水的存在，距離七號煤35 m處的奧灰?guī)r承受的有效應(yīng)力水平很低，特別是開挖卸荷后，石灰?guī)r中出現(xiàn)較大范圍的低應(yīng)力區(qū)，應(yīng)力狀態(tài)的改變致使奧灰?guī)r容易發(fā)生破壞，可能使巖體中的裂隙不斷發(fā)育并相互貫通，并形成承壓水導(dǎo)升帶。

3.2 推進(jìn)不同距離時(shí)頂?shù)装宕怪蔽灰品植?/h3>

圖4 推進(jìn)至不同距離時(shí)圍巖位移場Fig.4 Displacement field of surrounding rock at different distances

工作面推出開切眼50 m 時(shí)，頂板最大下沉量320 mm，頂板同充填體沒有充分接觸，因此充填體對底板的控制效果較差，此時(shí)，最大底鼓量為180 mm；工作面推進(jìn)至150 m 時(shí)，頂板最大下沉量達(dá)到720 mm，頂板同充填體之間的空頂高度已不存在，基巖同充填體密實(shí)接觸，充填體被壓實(shí)，對底板的控制效果轉(zhuǎn)好，此時(shí)底板最大鼓起高度為120 mm，說明在充填體的作用下，已采空間頂?shù)装逦灰剖艿郊s束，卸壓范圍減小。

3.3 推進(jìn)至不同距離底板破壞深度

工作面推出開切眼不同距離時(shí)，采場圍巖塑性破壞區(qū)分布如圖5。

圖5 推進(jìn)至不同距離時(shí)圍巖塑性區(qū)分布Fig.5 Distribution of surrounding rock plastic zone at different distances

由于充填作用，底板破壞深度較小，且隨著工作面的推進(jìn)，底板破壞深度并沒有沿垂直方向發(fā)展，推進(jìn)至50、150 m 時(shí)底板破壞深度大致相同，均為18 m 左右。在承壓水導(dǎo)升帶與采動(dòng)影響破壞帶之間存在厚17 m 的隔水關(guān)鍵層，說明在使用充填采煤工藝的條件下，A 工作面沒有突水危險(xiǎn)。

3.4 非充填開采情況下數(shù)值模擬結(jié)果對比

為進(jìn)一步對比充填開采的效果，進(jìn)行了非充填開采情況下的模擬計(jì)算。除未進(jìn)行工作面的充填開采外，其他條件與充填開采計(jì)算模型保持一致。非充填開采情況下工作面推進(jìn)到不同距離時(shí)的圍巖塑性區(qū)分布如圖6。

圖6 非充填開采工作面推進(jìn)至不同距離時(shí)圍巖塑性區(qū)分布Fig.6 Distribution of surrounding rock plastic zone when the non backfilling face advances to different distances

工作面推進(jìn)50 m 時(shí)，底板破壞深度在20 m左右，超出采用充填開采時(shí)底板破壞深度約2 m，此時(shí)采動(dòng)影響尚未波及底板奧灰含水層；當(dāng)工作面推進(jìn)150 m 后，隨著底板巖層塑性區(qū)的擴(kuò)展，受采動(dòng)影響的底板巖層破壞范圍已接近40 m，超出采用充填開采時(shí)底板破壞深度約12 m。

數(shù)值模擬結(jié)果表明，非充填開采情況下，工作面推進(jìn)到150 m 時(shí)，與充填開采相比，底板破壞深度明顯增加，已與奧灰含水層導(dǎo)通，極有可能誘發(fā)奧灰突水事故。也就是說，采用充填開采，可以有效減小煤層開采對底板的擾動(dòng)，減小底板破壞深度，有效防治底板突水事故。

4 底板采動(dòng)效應(yīng)的現(xiàn)場實(shí)測

A 工作面采用充填采煤工藝后，為得到在采動(dòng)和承壓水的雙重影響下，底板巖層的實(shí)際變形破壞特征，在距離A 工作面底板約29 m 處的八號煤泄水巷道中超前工作面安裝動(dòng)態(tài)儀，共設(shè)置12個(gè)測站，通過記錄動(dòng)態(tài)儀數(shù)據(jù)，觀測泄水巷道的頂?shù)装逡平考跋锏赖钠茐那闆r。測點(diǎn)2 至測點(diǎn)12 等11 個(gè)測點(diǎn)的頂?shù)装逡平孔兓€（測點(diǎn)1由于出現(xiàn)故障，未獲得有效數(shù)據(jù)），測點(diǎn)頂?shù)装逡平侩S工作面推進(jìn)距離的變化曲線如圖7。

圖7 測點(diǎn)頂?shù)装逡平侩S工作面推進(jìn)距離的變化曲線Fig.7 Change curves of roof and floor movement with advancing distance of working face

由圖7 可以看出：在A 工作面推進(jìn)距離小于90 m 之前，各測點(diǎn)巷道的頂?shù)装逡平慷己苄。静皇芷溟_采活動(dòng)的影響；當(dāng)A 工作面推進(jìn)距離達(dá)到90 m 后，由于6#測點(diǎn)位于工作面正下方，且存在斷層，其頂?shù)装逡平繑?shù)值變化突然增大，受采動(dòng)影響較為劇烈；在A 工作面推進(jìn)距離達(dá)到100 m 時(shí)，位于A 工作面回風(fēng)巷正下方的4#測點(diǎn)數(shù)值也發(fā)生顯著變化，但頂?shù)装逡平繛?5 mm后變化不大。

以上代表性測點(diǎn)的觀測數(shù)據(jù)表明：在A 工作面開采初期，其開采活動(dòng)對八號煤泄水巷道基本沒有影響，而當(dāng)A 工作面推進(jìn)90 m 后，個(gè)別測點(diǎn)由于A 工作面采空區(qū)正下方或處于A 工作面回風(fēng)巷正下方，受到較大影響；但是在A 工作面推進(jìn)過程中，泄水巷道并沒有出現(xiàn)很大變形，一直可以利用，說明在充填開采條件下，七號煤的采動(dòng)影響并不會(huì)造成底板突水，可以實(shí)現(xiàn)安全回采。

5 結(jié) 語

1）結(jié)合巖石力學(xué)、土力學(xué)及彈性力學(xué)基本原理，建立了承壓水上充填開采底板破壞深度分析的力學(xué)模型，并推導(dǎo)出了底板最大破壞深度的求解公式，由該公式可計(jì)算得出A 工作面充填開采時(shí)底板最大破壞深度為20 m，發(fā)生在工作面后方，工作面沒有突水危險(xiǎn)，可實(shí)現(xiàn)安全回采。

2）通過數(shù)值模擬得到了充填體對采場圍巖應(yīng)力場、位移場及塑形破壞區(qū)域的控制作用。在充填體作用下，采場圍巖的采動(dòng)卸壓范圍、頂?shù)装逡平棵黠@減小，圍巖應(yīng)力集中程度降低，底板破壞深度并未隨著推進(jìn)距離的增加而持續(xù)縱向發(fā)展，底板破壞深度約18 m，同理論計(jì)算結(jié)果基本吻合。

3）A 工作面推進(jìn)過程中，對底板中八號煤層泄水巷道的變形破壞特征進(jìn)行了現(xiàn)場實(shí)測，在七號煤采動(dòng)影響下，泄水巷道沒有因發(fā)生較大變形而廢棄，只有局部出現(xiàn)片幫冒頂現(xiàn)象，說明充填開采條件下，A 工作面的底板采動(dòng)效應(yīng)得到很大程度的緩和。