李鶴鶴，冀宇鑫，宋高峰

(北方工業(yè)大學(xué) 土木工程學(xué)院，北京 100144)

0 引言

充填采煤技術(shù)對覆巖移動和地表沉降具有明顯的控制作用，同時也是綠色開采技術(shù)的重要組成部分[1]。充填采煤技術(shù)是采煤過程結(jié)合充填過程同時進(jìn)行的一項采煤技術(shù)，其主要把采煤過程中產(chǎn)生的廢棄物矸石作為充填材料，對采煤工作面后方采空區(qū)進(jìn)行填充，不僅可以對采煤過程中產(chǎn)生的矸石進(jìn)行重新利用，而且可以有效解決由矸石帶來的占用耕地、污染水質(zhì)、破壞土壤等環(huán)境污染問題[2-3]。

充填開采技術(shù)的廣泛運(yùn)用對生態(tài)環(huán)境保護(hù)具有重要意義，符合安全生產(chǎn)、生態(tài)文明建設(shè)及保護(hù)環(huán)境的發(fā)展理念[4]。學(xué)者們在充填開采技術(shù)方面進(jìn)行了不同角度的研究分析，左建平等[5]通過建立充填開采巖層曲率模型，對充填開采下煤體覆巖連續(xù)變形移動規(guī)律進(jìn)行了分析；許猛堂等[6]根據(jù)超高水材料的基本性能試驗，利用數(shù)學(xué)與專業(yè)理論建立動力學(xué)模型，進(jìn)而對超高水材料長壁工作面充填開采頂板活動規(guī)律進(jìn)行了分析；曹正正等[7]通過運(yùn)用FLAC軟件建立數(shù)值計算模型，對不同的充實率條件下覆巖移動與地表下沉的規(guī)律進(jìn)行了分析；趙相海等[8]通過分析地表沉陷機(jī)理及充填開采覆巖破壞特征，得到影響地表沉陷主要的3個因素并進(jìn)一步進(jìn)行了分析；劉長友等[9]通過結(jié)合現(xiàn)場地質(zhì)條件，對充填開采時上覆巖層的活動規(guī)律和穩(wěn)定性進(jìn)行了分析；賈林剛等[10]對長壁充填開采充填體穩(wěn)定性進(jìn)行了分析；王荀勇[11]對矸石充填開采技術(shù)對地表變形的影響進(jìn)行了分析；王永杰等[12]對分階段充填開采覆巖應(yīng)力演化和地表減沉效果進(jìn)行了分析。

前人研究雖然取得了豐碩的成果，對充填體不同剛度條件下垂直應(yīng)力及巖層移動規(guī)律的研究仍可進(jìn)一步深入。為此，擬采用PHASE 2D有限元軟件，對充填體不同剛度條件下采場應(yīng)力分布及覆巖移動規(guī)律進(jìn)行研究。

1 數(shù)值模型的建立

1.1 工況概括

數(shù)值模擬以某煤礦11030大采高綜采工作面為工程背景，研究了不同充填體剛度條件下長壁工作面前方煤體塑性區(qū)演化規(guī)律、工作面前方支承壓力分布規(guī)律和覆巖變形演化規(guī)律。大采高工作面傾斜長度161 m，走向長度1 600 m，煤層厚度5.39～6.93 m，煤層傾角2.8°～13.5°，煤體內(nèi)生裂隙發(fā)育，煤質(zhì)堅硬。

1.2 模型介紹

運(yùn)用PHASE 2D有限元軟件建立數(shù)值模型，模型平面形狀為長方形，水平方向長400 m，垂直方向高100 m，如圖1所示。模型左、右垂直方向邊界施加水平方向鏈桿約束，下部邊界施加垂直向上鏈桿約束，防止底部發(fā)生向下垂直位移。建立數(shù)值模型中總共包括13 813個單元，27 235個節(jié)點(diǎn)。模型工作面推進(jìn)為200 m，推進(jìn)范圍為水平方向100～300 m。模擬煤層深度為270 m，煤層厚度為3 m，左右兩側(cè)各留100 m煤柱，煤層底部距離模型底端30 m。由于模型中煤層上方有70 m厚的巖層，故在模型頂部施加了5.5 MPa的補(bǔ)償應(yīng)力，以模擬模型上方200 m厚巖層的重量[13]。建立了3個數(shù)值模型，分別研究充填體彈性模量為2 000 MPa、3 000 MPa和4 000 MPa時采場應(yīng)力分布及巖層移動變化規(guī)律。

圖1 PHASE 2D數(shù)值模型Fig.1 PHASE 2D numerical model

模型地層結(jié)構(gòu)從上到下依次為粉砂巖、粗砂巖、細(xì)粉砂巖、泥巖、煤層、泥巖、粉砂巖、石灰?guī)r。模型采用霍克-布朗強(qiáng)度準(zhǔn)則，各個巖層煤巖體參數(shù)見表1。其中，v為泊松比；σci為巖石單軸抗壓強(qiáng)度，MPa；GSI為地質(zhì)強(qiáng)度系數(shù)；mi為巖石常數(shù)；mb為巖石折減系數(shù)；s，a為巖體常數(shù)；Em為彈性模量，MPa。

表1 煤巖霍克-布朗強(qiáng)度參數(shù)Table 1 Hoek-Brown strength parameters of coal

2 數(shù)值模型結(jié)果

2.1 工作面前方煤體塑性區(qū)云圖分析

模型推進(jìn)200 m后，工作面前方煤體塑性區(qū)如圖2所示。工作面煤體塑性區(qū)形狀大體呈現(xiàn)為上部寬、下部窄的特征。當(dāng)彈性模量E為2 000 MPa時，工作面煤體的塑性區(qū)寬度為1.98 m；當(dāng)彈性模量增大到3 000 MPa時，工作面煤體的塑性區(qū)寬度為1.54 m；當(dāng)彈性模量進(jìn)一步增大為4 000 MPa時，工作面煤體的塑性區(qū)寬度仍然為1.54 m左右。因此，隨著彈性模量的增大，工作面前方煤體塑性區(qū)的寬度不斷減小，然而當(dāng)塑性區(qū)寬度降低到一定程度時，增大充填體彈性模型很難進(jìn)一步降低工作面煤體塑性區(qū)寬度。

圖2 不同剛度工作面前方煤體塑性區(qū)規(guī)律Fig.2 Distribution law of plastic zone of coal mass in front of working face with different stiffness

2.2 垂直應(yīng)力演化規(guī)律

推進(jìn)200 m后，工作面前方支承壓力分布規(guī)律如圖3所示。

圖3 工作面前方支承壓力分布規(guī)律Fig.3 Distribution law of support pressure in front of working face

在工作面前方，支承壓力增高系數(shù)先增大后降低，在距離工作面煤壁50～100 m以外，支承壓力逐漸恢復(fù)到原巖應(yīng)力水平。充填體剛度為2 000 MPa、3 000 MPa、4 000 MPa時，對應(yīng)的支承壓力集中系數(shù)峰值分別為1.8、1.71、1.67。隨著彈性模量的增大，工作面前方支承壓力整體表現(xiàn)為逐漸減小，但支承壓力的減小程度逐漸降低。

2.3 工作面支承壓力峰值

圖4給出了工作面推進(jìn)距離與支承壓力峰值增高系數(shù)的關(guān)系。隨著工作面的推進(jìn)，支承壓力峰值表現(xiàn)為先逐漸增大后趨于穩(wěn)定的變化規(guī)律。當(dāng)彈性模量為2 000 MPa時，支承壓力峰值增高系數(shù)最大值為1.8；當(dāng)彈性模量為3 000 MPa時，支承壓力峰值增高系數(shù)最大值為1.71；當(dāng)彈性模量為4 000 MPa時，支承壓力峰值增高系數(shù)最大值為1.67。隨著充填體剛度的增大，支承壓力峰值整體有所減小，且隨著充填體剛度的等額增加，支承壓力峰值整體減小幅度有所降低。因此，當(dāng)充填體具有較大剛度時，采空區(qū)能夠分擔(dān)更多的上覆巖層自重，有利于保障圍巖穩(wěn)定性，但隨著充填體剛度的進(jìn)一步增大，這一趨勢有所減弱。

圖4 充填體不同剛度下支承壓力峰值發(fā)展規(guī)律Fig.4 Development law of support pressure under different stiffness of filling materials

2.4 工作面前后方直接頂下沉量

工作面推進(jìn)200 m，工作面前后方5 m范圍內(nèi)的直接頂下沉量如圖5所示。不難看出，在工作面前方，距離工作面越遠(yuǎn)，直接頂下沉量越來越小；在工作面后方，距離工作面越遠(yuǎn)，直接頂下沉量越來越大。當(dāng)彈性模量為2 000 MPa時，工作面后方直接頂下沉量最大值為20 mm，工作面前方直接頂下沉量逐漸減小為0；當(dāng)彈性模量增大到3 000 MPa、4 000 MPa時，在工作面前后方直接頂下沉量曲線幾乎重合，且工作面后方直接頂下沉量最大值為17 mm，說明在這2種充填體剛度條件下，對直接頂下沉量的影響不再明顯。

圖5 工作面前后方直接頂下沉量Fig.5 The amount of direct roof subsidence in front and rear of the working face

2.5 不同剛度充填體垂直位移變化規(guī)律

工作面推進(jìn)200 m，充填體不同剛度條件下覆巖變形規(guī)律如圖6所示。垂直位移呈現(xiàn)中間大、兩端小的變形規(guī)律，其中最大位移出現(xiàn)在采空區(qū)的中部位置，垂直位移曲線關(guān)于最大垂直位移對稱。當(dāng)彈性模量E為2 000 MPa時，垂直位移最大值為45 mm；當(dāng)彈性模量為3 000 MPa時，垂直位移最大值為40 mm；當(dāng)彈性模量4 000 MPa時，垂直位移最大值為37 mm。隨著充填體剛度的增大，垂直位移整體有所減小，但等額增大充填體剛度值，對垂直位移的影響有所降低。

圖6 充填體不同剛度下垂直位移曲線Fig.6 Vertical displacement curve of filling materials with different stiffness

3 結(jié)論

(1)當(dāng)彈性模量分別為2 000 MPa、3 000 MPa、4 000 MPa時，對應(yīng)工作面煤體的塑性區(qū)寬度分別為1.98 m、1.54 m、1.54 m。因此，隨著彈性模量的增大，工作面前方煤體塑性區(qū)的寬度不斷減小，并且當(dāng)彈性模量達(dá)到3 000 MPa，繼續(xù)增大充填體彈性模量幾乎不再對工作面媒體塑性區(qū)寬度有明顯影響。

(2)當(dāng)彈性模量分別為2 000 MPa、3 000 MPa、4 000 MPa時，對應(yīng)支承壓力峰值增高系數(shù)最大值分別為1.8、1.71、1.67。增大充填體彈性模量，支承壓力峰值整體表現(xiàn)為有所減小，但隨著充填體彈性模量的等額增加，支承壓力峰值整體減小幅度有所降低，說明隨著充填體彈性模量的增大，對支承壓力峰值的影響程度有所降低。

(3)工作面推進(jìn)200 m后，當(dāng)彈性模量E為2 000 MPa時，直接頂下沉量最大值為20 mm，垂直位移最大值為45 mm；當(dāng)彈性模量增大為3 000 MPa、4 000 MPa時，對應(yīng)直接頂下沉量曲線幾乎一樣，直接頂下沉量最大值都為17 mm，垂直位移最大值分別為40 mm、37 mm。因此，增大充填體彈性模量可有效控制直接頂下沉量和覆巖垂直位移，但隨著充填體彈性模量的增大，對直接頂下沉量及覆巖垂直位移的影響程度不斷降低直至不再明顯。