許海濤,李永軍,康慶濤
(1.河北省礦井災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室,河北 三河065201;2.華北科技學(xué)院安全工程學(xué)院,河北三河065201;3.中國(guó)礦業(yè)大學(xué) (北京)煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100083)
華北地區(qū)煤田水文地質(zhì)條件十分復(fù)雜,裂隙以斜交和垂直地層等形式賦存于各類(lèi)地層中,分布面積大,裂隙及溶蝕裂隙發(fā)育,含水性強(qiáng),礦井涌水量大、突水頻率高,煤礦安全生產(chǎn)形勢(shì)嚴(yán)峻。對(duì)鄭州、晉城、徐州、永城及兩淮地區(qū)50多個(gè)礦井的調(diào)研及井下觀測(cè),裂隙是礦井突水主要的導(dǎo)水通道,且多數(shù)淹井、淹工作面的重大突水事件與底板裂隙發(fā)育有關(guān)。煤層底板突水是煤層下伏承壓水沖破底板隔水層的阻隔,以突水、緩發(fā)或滯發(fā)的形式進(jìn)入工作面,造成礦井涌水量增加或淹井的自然災(zāi)害,為礦井生產(chǎn)的五大災(zāi)害之一[1-2]。幾十年來(lái),國(guó)內(nèi)外許多學(xué)者對(duì)礦井突水機(jī)理進(jìn)行了一些有益的探索,取得了大量的研究成果。“突水系數(shù)”[3]、“強(qiáng)滲通道”[4]、“水巖應(yīng)力”[5]、“零位破壞與原位張裂”、“關(guān)鍵層”及“下三帶”等理論學(xué)說(shuō),都從各個(gè)方面揭示了突水發(fā)生的機(jī)理和預(yù)測(cè)方法,對(duì)于礦井安全生產(chǎn)起到了積極的指導(dǎo)作用[6-7]。
隨著我國(guó)一些重點(diǎn)礦區(qū)的奧陶紀(jì)灰?guī)r巖溶水突出問(wèn)題日益受到重視,對(duì)煤層底板承壓水的突水機(jī)理已展開(kāi)了廣泛的研究,礦山壓力是影響底板承壓水突出的重要因素,是底板破壞的力源;了解煤層底板的應(yīng)力分布是研究煤層底板破壞深度的重要前提,有利于加深對(duì)突水機(jī)理的認(rèn)識(shí)。本文擬采用FLAC模擬軟件求解采空區(qū)的圍巖應(yīng)力-應(yīng)變關(guān)系,分析不同階段采空區(qū)底板應(yīng)力變化過(guò)程,并建立應(yīng)力、應(yīng)變隨采場(chǎng)深度變化的本構(gòu)關(guān)系,以解釋其對(duì)底板突水的影響。
工作面開(kāi)挖以后的瞬間,在采場(chǎng)周邊的初始地應(yīng)力被釋放,應(yīng)力由三向狀態(tài)變?yōu)閮上蚧騿蜗驊?yīng)力狀態(tài),使采場(chǎng)周邊附近的圍巖的應(yīng)力重新分布,即底板受力隨著距采場(chǎng)距離而逐漸加大到到原巖應(yīng)力狀態(tài)。由于圍巖應(yīng)力的重新分布引起底板出現(xiàn)變形,在一定條件下,底板的變形發(fā)展導(dǎo)致圍巖破壞、失穩(wěn)。
開(kāi)挖前,底板巖層處于原巖應(yīng)力狀態(tài),如圖1(A)、圖2所示,且視作均布載荷,簡(jiǎn)化底板受力模型。開(kāi)挖后瞬間底板巖層受力如圖1(B)、圖3所示。工作面開(kāi)挖以后的瞬間,采空區(qū)底板巖層上部支承力被去掉,但深部底板巖層的彈塑性恢復(fù)應(yīng)力與原巖應(yīng)力垂直分量相等且均勻分布,底板有向虛線位置發(fā)生位移的趨勢(shì)。由于采空區(qū)底板巖層下部支承力被去掉,而且底板有向虛線位置發(fā)生位移的趨勢(shì),致使底板應(yīng)力平衡被打破,導(dǎo)致在采場(chǎng)周邊的初始地應(yīng)力被釋放,圍巖受力隨著距采空區(qū)底板距離加大而逐漸加大,然后逐漸減小到原巖應(yīng)力狀態(tài)。隨著圍巖應(yīng)力調(diào)整,圍巖的應(yīng)力重新分布成為圖1(CD)、圖4所示應(yīng)力狀態(tài)。
圖1 采場(chǎng)垂直應(yīng)力變化圖
通過(guò)數(shù)值模擬計(jì)算結(jié)果可得底板不同深度測(cè)點(diǎn)的垂直應(yīng)力變化曲線如圖5所示,底板垂直應(yīng)力與底板深度關(guān)系曲線如圖6所示。當(dāng)采場(chǎng)應(yīng)力重新分布基本平衡時(shí),底板處的垂直應(yīng)力為0MPa,隨著深度的增加,垂直應(yīng)力逐步增加,根據(jù)以上數(shù)據(jù)統(tǒng)計(jì),運(yùn)用SPSS數(shù)學(xué)軟件進(jìn)行回歸分析,可得到式(1)。
圖2 未開(kāi)挖時(shí)底板受力分析圖
圖5 底板不同位置處垂直應(yīng)力變化曲線
圖6 底板垂直應(yīng)力與底板深度關(guān)系曲線
式中,Y為垂直應(yīng)力與原巖應(yīng)力的比值;X為采場(chǎng)底板深度,m。
此時(shí)垂直應(yīng)力
式中P0為原巖應(yīng)力,MPa。
當(dāng)y=1時(shí),x的值代表工作面回采對(duì)底板的影響極限深度;當(dāng)y>1時(shí),底板垂直應(yīng)力可由原巖應(yīng)力P0求得。由公式(1)可知:當(dāng)采空區(qū)下方隔水層深度一定時(shí)、厚度一定時(shí),隔水層由于采動(dòng)而引起的力學(xué)變化,如當(dāng)隔水層為底板下方20~50m,則y20=0.06,y50=0.38,即隔水層上部垂直應(yīng)力為:p1=y(tǒng)20×p0,下部垂直應(yīng)力p2=y(tǒng)50×p0。
采空區(qū)底板垂直位移圖見(jiàn)圖7。通過(guò)數(shù)值模擬計(jì)算結(jié)果可得,底板不同深度測(cè)點(diǎn)的垂直位移變化曲線如圖8所示,底板深度與位移關(guān)系曲線如圖9所示。當(dāng)采場(chǎng)應(yīng)力重新分布基本平衡時(shí),底板處的垂直變形最大為104.6mm,隨著深度的增加,垂直位移逐步減??;根據(jù)以上數(shù)據(jù)統(tǒng)計(jì)進(jìn)行回歸分析,其數(shù)學(xué)表達(dá)式見(jiàn)式(2)。
此時(shí)垂直方向的應(yīng)變
圖7 采空區(qū)底板垂直位移圖
圖8 底板下方不同深度垂直位移變化曲線
圖9 底板深度與位移關(guān)系曲線圖
式中:y為底板深度為x時(shí)的垂直位移,m;x為采場(chǎng)底板深度,m;ε為底板深度為x時(shí)的應(yīng)變;[ε]為底板巖體的屈服應(yīng)變。
當(dāng)ε>[ε]時(shí),底板將產(chǎn)生屈服破壞;當(dāng)ε<[ε]時(shí),底板未破壞。
式(2)為本文中數(shù)值模擬結(jié)果的統(tǒng)計(jì)分析,尚未考慮其他因素的影響。
通過(guò)以上分析可知以下內(nèi)容。
1)工作面開(kāi)挖以后的瞬間,在采場(chǎng)周邊的初始地應(yīng)力被釋放,應(yīng)力由三向狀態(tài)變?yōu)閮上蚧騿蜗驊?yīng)力狀態(tài),使采場(chǎng)周邊附近的圍巖的應(yīng)力重新分布,在一定條件下,底板的變形發(fā)展導(dǎo)致圍巖破壞、失穩(wěn)。
2)建立了深度與應(yīng)力的本構(gòu)關(guān)系,當(dāng)采場(chǎng)應(yīng)力重新分布基本平衡時(shí),底板處的垂直應(yīng)力接近于0MPa,隨著深度的增加,垂直應(yīng)力逐步增加。
3)建立了底板深度與應(yīng)變的本構(gòu)關(guān)系,當(dāng)采場(chǎng)應(yīng)力重新分布基本平衡時(shí),底板處的垂向變形最大,隨著深度的增加,垂向位移逐步減小。
4)本公式只是對(duì)本計(jì)算模型數(shù)據(jù)的回歸分析,尚未考慮時(shí)間因素,尚需進(jìn)一步深入研究。
[1]楊天鴻,唐春安,劉紅元,等.承壓水底板突水失穩(wěn)過(guò)程的數(shù)值模型初探[J].地質(zhì)力學(xué)學(xué)報(bào),2003,9(3):281-288.
[2]魏久傳,李白英.承壓水上采煤安全性評(píng)價(jià)[J].煤田地質(zhì)與勘探,2000,28(4):57-59.
[3]沈光寒,李白英,吳戈.礦井特殊開(kāi)采的理論與實(shí)踐[M].第1版.北京:煤炭工業(yè)出版社,1992:285-290.
[4]許學(xué)漢,王杰.煤礦突水預(yù)報(bào)研究[M].第1版.北京:地質(zhì)出版社,1991:228-235.
[5]吳基文,樊成.煤層底板巖體阻水能力原位測(cè)試研究[J].巖土工程學(xué)報(bào),2003,25(1):67-70.
[6]王作宇,劉鴻泉.承壓水上采煤[M].北京:煤礦工業(yè)出版社,1992:102-113.
[7]荊自剛.峰峰二礦開(kāi)采活動(dòng)與底板突水關(guān)系研究[J].煤炭學(xué)報(bào),1984,9(2):20-23.