陜北典型地質(zhì)條件下煤炭開采方法的數(shù)值模擬研究

張 澤

(陜西匯森煤業(yè)開發(fā)有限責任公司，陜西 西安 710065)

0 引言

煤炭作為我國工業(yè)發(fā)展的支柱，在過去的幾十年中，產(chǎn)量是其行業(yè)發(fā)展的主要衡量指標。然而，以產(chǎn)量定輸贏的粗放型開采模式造成了眾多環(huán)境及安全問題。當前，隨著社會的進步，如何實現(xiàn)煤炭的綠色高效開采對于我國國民經(jīng)濟的發(fā)展至關(guān)重要[1-2]。走向長壁開采作為近年來我國各大礦井最常應用的開采方法，能夠最大化地提升煤炭產(chǎn)量是其主要優(yōu)勢，但是由于其開采方法近似是將地層中的煤炭進行整層剝離，由此導致的地表沉陷、地下水斷流及地表植被匱竭業(yè)已成為綠色開采的桎梏所在[3-4]。

在國家既需要煤炭又需要煤炭開采不破壞環(huán)境的大前提下，采用對環(huán)境破壞程度較小的條帶式采煤方法成為一種可行方式，而條帶采煤方法中煤柱的穩(wěn)定性是學者們研究的重點。劉貴[5-7]運用彈塑性理論優(yōu)化了A.H Wilson[8]計算公式，定量化地提出了條帶開采的煤柱尺寸表達式并在現(xiàn)場中得到了檢驗。郭惟嘉等[9-10]運用物理相似模擬及現(xiàn)場實測的方式，研究了深井寬條帶開采中的地表沉陷及煤柱穩(wěn)定性問題。朱衛(wèi)兵等[11-13]通過運用離散元數(shù)值模擬手段，分析了村莊下煤礦的條帶開采方案，得出了條帶采寬與垂直支承應力間的線性關(guān)系，從而為礦井的高效生產(chǎn)提供了技術(shù)支撐。

針對我國西部地區(qū)水資源匱乏的現(xiàn)狀，基于厚松散表土層及淺埋煤層這一陜北典型地質(zhì)條件為基礎(chǔ)，運用有限元數(shù)值計算軟件FLAC3D，針對涼水井煤礦4-2煤層，對比分析條帶開采與長壁開采覆巖結(jié)構(gòu)的演化規(guī)律，并重點分析采用不同采、留寬參數(shù)條件下條帶煤柱的穩(wěn)定性，從而為類似礦井的綠色回采給予理論支撐。

1 礦井及工作面概況

涼水井煤礦位于陜西神木市境內(nèi)，其礦區(qū)地表為一巨厚松散層，平均厚度40 m。4-2煤層為礦井主采煤層，平均埋深110 m，均厚3.0 m，傾角0°～3°，為淺埋煤層。基本頂主要為細粒砂巖及粉砂巖，巖性以粉砂巖—細粒砂巖為主，厚度0.55～26.72 m，平均7 m，偶見泥巖偽頂，屬于穩(wěn)定頂板。4-2煤層直接底巖性以粉砂巖為主，厚度0.51～24.03 m，平均12 m。

模擬研究42113工作面。工作面長度310 m，走向長度4 380 m。由于地表松散層較厚且煤層埋深較淺，在采用走向長壁采煤法回采時，地表沉陷較為明顯并偶爾伴隨有潰沙現(xiàn)象，回采形成的覆巖裂縫已經(jīng)貫穿至地表，當有地面降水時，由于此裂縫帶的存在，造成水土流失嚴重，雨水還會導入井下，造成水害。

2 模型構(gòu)建

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

數(shù)值計算主要通過對比同一地質(zhì)條件下走向長壁與條帶開采的區(qū)別來展開研究。模擬采用FLAC3D5.0有限元數(shù)值軟件。由于沿著回采方向，無論條帶開采還是走向長壁開采均呈現(xiàn)周期性來壓特征，故所構(gòu)建模型可看作沿著工作面回采方向的剖面，如圖1(a)所示A-A剖面。實驗設(shè)計模型長度300 m(沿工作面推進方向)，高度140 m，模擬4-2煤層埋深110 m，如圖1(b)所示。根據(jù)相關(guān)鉆孔信息及巖石力學試驗，對各煤巖層進行優(yōu)化賦參，煤巖層的物理力學參數(shù)見表1。

a-模型走向剖面位置示意圖；b-數(shù)值模型圖1 42113工作面地質(zhì)條件下模型構(gòu)建示意圖

2.2 約束條件

針對現(xiàn)場地質(zhì)情況，為保證模擬的合理有效，對模型兩側(cè)邊界進行水平位移約束。在此地質(zhì)條件下，運用走向長壁開采及條帶開采2種不同模式進行開挖。

走向長壁法開采：模型開挖步距為4 m，開挖一步運算一步，運算平衡標準為最大不平衡力小于10-3N，共連續(xù)開挖60步，開挖總步距240 m，最終開挖運算完成如圖3(a)所示。

表1 各煤巖層物理力學參數(shù)

條帶開采：如圖2(b)所示，采20 m留10 m，開采20 m運算一次至平衡，然后隔10 m再開采20 m，以此為一個周期共運算8次，開采加留設(shè)總長度240 m；如圖2(c)所示，采12 m留8 m，開采12 m運算一次至平衡，然后隔8 m再開采12 m，以此為一個周期共運算12次，開采加留設(shè)總長度240 m；如圖2(d)所示，采10 m留10 m，開采10 m運算一次至平衡，然后隔10m再開采10 m，以此為一個周期共運算12次，開采加留設(shè)總長度240 m。

a-走向長壁開采；b-采20 m留10 m；c-采12 m留8 m；d-采10 m留10 m圖2 不同開采條件下模型開挖運算至平衡圖

3 模擬結(jié)果分析

3.1 模型建立的合理準確性

圖3為未開采時的煤巖層賦參至平衡垂直應力云圖。由圖3可知，在數(shù)值模型建立以后，各煤巖層受力狀態(tài)均勻，垂直應力等值線呈近似直線狀態(tài)分布，表明模型中沒有未耦合點導致的應力集中現(xiàn)象出現(xiàn)，巖層應力逐層傳遞，從而可知模型建立的科學合理，人為因素的建模誤差極小。

圖3 未開采條件下煤巖層賦參后運算平衡垂直應力云圖

3.2 垂直應力分布規(guī)律

應力場的變化是覆巖結(jié)構(gòu)變化以及裂隙帶產(chǎn)生的根源，因此，要分析走向長壁與條帶開采的本質(zhì)區(qū)別，必須掌握二者應力場的分布特征及其同異點。

圖4為不同開采條件下煤巖垂直應力分布云圖。表明工作面無論以何種方式回采，由于其開采擾動，必然產(chǎn)生應力集中現(xiàn)象。2種開采模式中，均為采空區(qū)頂、底板中部區(qū)域所受應力最小，兩支承點所受應力最大。

走向長壁開采中的應力分布呈現(xiàn)出駝峰狀的特征，兩支承點最大受力區(qū)域達6 m寬，最大受力區(qū)域的垂向應力在21.2～20.0 MPa之間。除此之外，將圖4(a)與圖3對比，采空區(qū)上方覆巖的受力也不再有均勻的等值線，表明開采對上覆巖層的受力狀態(tài)影響明顯。而條帶開采中，采空區(qū)的應力分布狀態(tài)呈現(xiàn)出周期性的拱形分布，類似于橋下的孔洞，這與條帶開采后的圍巖結(jié)構(gòu)及其覆巖支撐密切相關(guān)。由于覆巖有多個窄條帶煤柱進行支承，因此，雖然各支承點的應力分布明顯大于未開采區(qū)域，但與長壁開采要小很多，其最大垂向應力在6.48～8.5 MPa之間，其支承應力僅為長壁開采的32.4%～40.1%。圖4(b)～(d)與圖3對比表明，位于采空區(qū)上方50 m及以上位置，依然有明顯的近似直線的應力等值線存在，表明松散層及最上部的基巖僅僅受到極小的開采擾動的影響。

將條帶開采中的單個煤柱取出，對比分析3種條帶開采參數(shù)下煤柱的受力特征如圖5所示。采20 m留10 m的煤柱內(nèi)部受力明顯呈現(xiàn)出駝峰狀特征，最大應力達8.5 MPa，表明煤柱兩側(cè)受力相對較大。采12 m留8 m與采10 m留10 m煤柱內(nèi)部受力狀態(tài)類似于梯形。采12 m留8 m煤柱內(nèi)部應力最小，僅為4 MPa。因此，從受力角度分析，采12 m留8 m=采10 m留10 m>采20 m留10 m。

3.3 覆巖變形及地表下沉規(guī)律

位移場的變化是覆巖運移以及地表下沉最直接的體現(xiàn)，地表下沉位移量越小，則地面產(chǎn)生沉陷進而導致水土流失的可能性也越小。因此，分析不同開采方法下的位移場變化特征對于保障綠色開采極為關(guān)鍵。

a-走向長壁開采；b-采20 m留10 m；c-采12 m留8 m；d-采10 m留10 m圖4 不同開采條件下垂直應力分布云圖

圖6為不同開采條件下位移分布云圖。由圖6可知，走向長壁開采平衡后，其位移等值線圖呈現(xiàn)出明顯的急傾斜特征，其中采空區(qū)中部甚至已經(jīng)接近垂直，表明厚表土松散層及淺埋煤層條件下，長壁開采后的覆巖頂板會產(chǎn)生周期性的切落，表現(xiàn)出臺階下沉現(xiàn)象，因此才會出現(xiàn)采空區(qū)上覆直接頂?shù)奈灰婆c地表位移一致的現(xiàn)象。在此基礎(chǔ)上結(jié)合圖4(a)可知，其支撐區(qū)域應力極大正是由于這種臺階下沉式的切落導致的，這是由于采空區(qū)覆巖所形成的短砌體梁鉸接結(jié)構(gòu)的穩(wěn)定性較小導致覆巖整體下沉。因此，雖說為淺埋煤層，在此地質(zhì)條件下運用走向長壁法開采也極有可能出現(xiàn)明顯的動壓現(xiàn)象。由圖7也可以看出，走向長壁的地表位移量遠遠大于條帶開采的位移量，以二者最大位移量的比值來看，走向長壁的地表下沉是條帶開采最大下沉量的6.428倍。

條帶開采中，采空區(qū)覆巖的位移分布狀態(tài)越往淺部越呈現(xiàn)出周期性近似直線的弧線形分布(圖6(b)～(d))，其中，采20 m留10 m的最大下沉為28 mm，采12 m留8 m最大下沉為25 mm，采10 m留10 m最大下沉為23 mm(如圖7)。尤其是采10 m留10 m的方式，地表下沉曲線已經(jīng)極為類似于一條直線，這表明條帶開采能夠極為有效地遏制地表沉陷，而且其覆巖依舊存在極為穩(wěn)定的結(jié)構(gòu)，松散層未明顯受到開采擾動的影響。

3.4 綜合分析

在涼水井煤礦現(xiàn)采地質(zhì)條件基礎(chǔ)上，綜合覆巖變形及應力變化特征，條帶開采能夠有效地防治巖層移動導致的地表沉陷現(xiàn)象，從而防止水土流失，實現(xiàn)綠色開采。而就采12 m留8 m、采10 m留10 m、采20 m留10 m這3種開采方式而言，采20 m留10 m與采12 m留8 m的經(jīng)濟效益最好，采12 m留8 m的煤柱穩(wěn)定性最好，采10 m留10 m與采12 m留8 m的地表下沉程度最輕，故綜合考慮，推薦采12 m留8 m的條帶開采方式。

a-走向長壁開采；b-采20 m留10 m；c-采12 m留8 m；d-采10 m留10 m圖6 不同開采條件下垂直應力分布云圖

圖7 不同開采方法下地表位移分布云圖

4 結(jié)論

(1)在厚表土松散層及淺埋煤層條件進行長壁開采，采空區(qū)上覆直接頂?shù)奈灰婆c地表位移極為相近，原因是采空區(qū)覆巖所形成的短砌體梁鉸接結(jié)構(gòu)的穩(wěn)定性較差導致的覆巖周期性的切落，從而極易發(fā)生明顯的地表沉陷及工作面動壓顯現(xiàn)。

(2)條帶開采能有效遏制地表下沉，采空區(qū)上覆巖層在條帶開采條件下支撐點較多，覆巖依然能形成穩(wěn)定結(jié)構(gòu)，其開采擾動對表土基巖層影響極小。

(3)通過對采12 m留8 m、采10 m留10 m、采20 m留10 m 3種不同條帶開采參數(shù)方式下的經(jīng)濟、安全及環(huán)保因素綜合考慮，涼水井煤礦若采用條帶開采方法，運用采12 m留8 m的條帶開采參數(shù)最為合適，從而為涼水井煤礦及類似礦井開展綠色可持續(xù)開采提供一定的理論支撐。