蘇二換，陳 濤，鈕景付，劉 宇

（1.神華準(zhǔn)能集團(tuán)有限責(zé)任公司 哈爾烏素露天煤礦，內(nèi)蒙古 鄂爾多斯 017000；2.中國(guó)礦業(yè)大學(xué)，江蘇 徐州 221116）

進(jìn)入21世紀(jì)以來(lái)，隨著國(guó)家經(jīng)濟(jì)高速發(fā)展和對(duì)煤炭需求量的不斷增加，我國(guó)新疆、內(nèi)蒙古等地區(qū)的露天煤礦得到了快速發(fā)展。由于工程地質(zhì)條件的限制，煤炭開(kāi)采不斷向深部開(kāi)拓。在井工開(kāi)采煤炭資源的過(guò)程中，由于受充填開(kāi)采等技術(shù)經(jīng)濟(jì)成本的限制，一些煤礦在開(kāi)采后留下大量采空區(qū)[1-2]。特別在20世紀(jì)80年代，一些個(gè)人或集體為了節(jié)省開(kāi)采成本，只采富礦、采后不回填的開(kāi)采現(xiàn)象屢見(jiàn)不鮮，這樣在許多礦區(qū)附近留下了大量的不明采空區(qū)。在露天礦區(qū)附近的不明采空區(qū)受爆破震動(dòng)、大型機(jī)械碾壓擾動(dòng)及端幫排土場(chǎng)連續(xù)堆載等影響，造成露天礦采空區(qū)坍塌、地表沉陷及開(kāi)裂等現(xiàn)象發(fā)生[3-5]。由于地下采空區(qū)引發(fā)的露天礦道路塌陷、邊坡失穩(wěn)已經(jīng)嚴(yán)重威脅到露天開(kāi)采的生產(chǎn)安全[6-7]。因此，研究不均勻上覆荷載對(duì)露天礦采空區(qū)沉陷的影響規(guī)律，具有重要的理論意義和工程實(shí)用價(jià)值。

為探究此類(lèi)采空區(qū)沉陷的影響規(guī)律，許多學(xué)者運(yùn)用不同的方法對(duì)其開(kāi)展了研究[8-19]，并取得許多成果，這些研究成果為后續(xù)的露天礦境界內(nèi)的采空區(qū)沉陷研究奠定了理論基礎(chǔ)；但是露天礦區(qū)內(nèi)采空區(qū)的不連續(xù)性及形狀差異較大，對(duì)于不均勻上覆荷載對(duì)露天礦采空區(qū)（井工開(kāi)采余留下的廢棄巷道）沉陷的影響規(guī)律有待深入研究。為此，基于開(kāi)采沉陷理論與技術(shù)，建立采空區(qū)沉陷模型，通過(guò)建立數(shù)值分析了排土場(chǎng)堆載過(guò)程中產(chǎn)生的不均勻載荷對(duì)采空區(qū)、不同外排土場(chǎng)形狀對(duì)采空區(qū)的地表沉陷過(guò)程中的應(yīng)力分布形態(tài)及變形規(guī)律，為實(shí)現(xiàn)露天礦安全高效開(kāi)采提供保障。

1 采空區(qū)沉陷機(jī)理及模型分析

當(dāng)煤層開(kāi)采后，采空區(qū)周?chē)脑瓗r應(yīng)力受到擾動(dòng)，引起應(yīng)力重新分布，采空區(qū)頂板在上覆巖層的壓力及不均勻載荷下發(fā)生顯著的下沉、彎曲和離層，并由下向上發(fā)展至地表引起地表的移動(dòng)，從而造成地表沉陷現(xiàn)象[20]。在采空區(qū)巖層的下沉過(guò)程中，在豎直方向上，由于巖層自重和上覆巖層的壓力作用，破裂巖石逐漸被壓實(shí)，完整巖石被重新壓縮，地表的活動(dòng)逐漸趨于穩(wěn)定；同時(shí)，在工作面推進(jìn)方向上，采空區(qū)上覆巖層周期性重復(fù)下沉彎曲、破裂和移動(dòng)過(guò)程，造成地表塌陷區(qū)域不斷變大。

根據(jù)開(kāi)采沉陷巖體內(nèi)部移動(dòng)變形規(guī)律，建立露天礦排土場(chǎng)堆載下采空區(qū)上覆巖體的結(jié)構(gòu)模型，露天礦排土場(chǎng)下覆采空區(qū)巖層移動(dòng)模型如圖1。該模型反映了露天礦排土場(chǎng)地表沉陷穩(wěn)定后的覆巖結(jié)構(gòu)特征。模型最下部頂板破斷后的巖塊呈不規(guī)則垮落，巖層喪失了結(jié)構(gòu)連續(xù)性，稱(chēng)為垮落帶?？迓鋷е?，巖層雖然斷裂，但巖塊間仍然排列有序，保持自身“砌體梁”結(jié)構(gòu)，稱(chēng)為裂縫帶。在裂縫帶之上，各巖層保持原有力學(xué)性質(zhì)和結(jié)構(gòu)連續(xù)性，以整體形式彎曲下沉，稱(chēng)為彎曲帶，該帶的上下巖層在豎直方向上的下沉差值很小。最上部的是露天礦剝離物堆積而成的排土場(chǎng)，排土場(chǎng)將會(huì)對(duì)下覆巖層造成不均勻載荷，致使原本下沉穩(wěn)定的采空區(qū)發(fā)生二次沉陷，直接影響露天礦排土場(chǎng)邊坡的穩(wěn)定性。

圖1 露天礦排土場(chǎng)下覆采空區(qū)巖層移動(dòng)模型Fig.1 Rock strata movement model of overlying goaf in open-pit mine dump

露天礦的采空區(qū)種類(lèi)眾多，包括普通采空區(qū)、采空巷道、采空積水區(qū)和采空自燃區(qū)。對(duì)不同幾何尺寸的采空區(qū)巷道，建立采空區(qū)斷面力學(xué)模型，并將排土場(chǎng)堆積載荷附加到采空區(qū)上覆巖層表面。露天礦由于先前小煤窯的私挖亂采造成采空區(qū)幾何結(jié)構(gòu)多為“平頂”采空區(qū)，其力學(xué)模型屬于組合梁結(jié)構(gòu)，對(duì)于這種力學(xué)模型，其斷裂破壞過(guò)程中的力學(xué)機(jī)理是：層狀梁結(jié)構(gòu)破壞取決于懸空巖體的抗拉強(qiáng)度、搭接位置巖體的抗剪強(qiáng)度。

根據(jù)平頂采空區(qū)力學(xué)模型以及不均勻上覆荷載，分別進(jìn)行不同上覆載荷幾何狀態(tài)的采空區(qū)沉陷數(shù)值模擬過(guò)程中的應(yīng)力及變形位移計(jì)算，分析不同條件下采空區(qū)地表沉陷規(guī)律。采用FLAC3D軟件建立的平頂采空區(qū)力學(xué)模型如圖2，模型中采空區(qū)的寬度為30 m，高為6 m（根據(jù)某露天煤礦采空區(qū)勘探結(jié)果的平均尺寸[21]），將試驗(yàn)所得的力學(xué)參數(shù)賦入模型中，進(jìn)行不均勻上覆荷載下的采空區(qū)安全性模擬，分析其地表沉陷規(guī)律。為了防止地層產(chǎn)狀影響本研究對(duì)于不均勻上覆載荷對(duì)采空區(qū)沉陷普遍規(guī)律的認(rèn)識(shí)，特對(duì)地層產(chǎn)狀和數(shù)量進(jìn)行簡(jiǎn)化處理。煤巖物理力學(xué)參數(shù)匯總表見(jiàn)表1。

圖2 平頂采空區(qū)力學(xué)模型Fig.2 Mechanical model of flat-topped goaf

表1 煤巖物理力學(xué)參數(shù)匯總表Table 1 Summary of physical and mechanical parameters of coal and rock

2 堆載不平衡應(yīng)力對(duì)采空區(qū)地表沉陷的影響

在排土場(chǎng)堆載初期，由于堆載高度較小，排土場(chǎng)產(chǎn)生的不均勻載荷較低，傳遞到采空區(qū)上的應(yīng)力較低，對(duì)采空區(qū)的穩(wěn)定性影響較小。隨著排土場(chǎng)堆載高度增加，堆載體類(lèi)似一個(gè)臺(tái)體形狀，由堆載體自重應(yīng)力引起的豎向應(yīng)力及側(cè)向應(yīng)力不斷增大，當(dāng)傳遞到外排土場(chǎng)基底下方的采空區(qū)時(shí)，造成采空區(qū)所受擠壓應(yīng)力增大，導(dǎo)致采空區(qū)塌陷，直接影響礦山安全生產(chǎn)。

2.1 排土場(chǎng)堆載高度對(duì)采空區(qū)地表沉陷的影響

采用Rhino軟件在平頂采空區(qū)模型上建立堆載高度為0、20、40、60、80 m的排土場(chǎng)模型并進(jìn)行網(wǎng)格劃分，再導(dǎo)入FLAC3D數(shù)值模擬軟件進(jìn)行地表沉降位移、應(yīng)力模擬分析，分別在采空區(qū)正上方距離地表0、20、40、60、80、100、120、140、160 m處 布 置 監(jiān) 測(cè)點(diǎn)，得到不同位置的位移、應(yīng)力曲線，對(duì)比得到排土場(chǎng)堆載不同高度下豎直方向上觀測(cè)線下沉規(guī)律。排土場(chǎng)堆載不同高度下豎直方向上的地表沉降量曲線如圖3，排土場(chǎng)堆載不同高度下豎直方向上的應(yīng)力變化曲線如圖4。

圖3 排土場(chǎng)堆載不同高度下豎直方向上的地表沉降量曲線Fig.3 Surface settlement curves in the vertical direction under different heights of dump loading

圖4 排土場(chǎng)堆載不同高度下豎直方向上的應(yīng)力變化曲線Fig.4 Stress change curves in the vertical direction under different heights of dump loading

從地表沉降位移曲線可以看出，當(dāng)堆載高度從0 m增加到80 m的過(guò)程中，地表面的沉降量由24.66 mm增大到42.08 mm，受到排土場(chǎng)自然沉降影響較大，其中堆載高度80 m時(shí)較為明顯。排土場(chǎng)堆載高度40～80 m時(shí)，觀測(cè)點(diǎn)距離地表面40～80 m的監(jiān)測(cè)數(shù)據(jù)較為接近，是由于砂巖層在下沉彎曲變形過(guò)程中保持整體結(jié)構(gòu)連續(xù)性，起到關(guān)鍵層作用。隨著排土場(chǎng)堆載高度的加大，采空區(qū)上覆不均勻載荷隨之增大，地表沉降量逐漸增大，其中采空區(qū)頂板的下沉累計(jì)量由未堆載排棄物的38.15 mm增加到堆載高度為80 m的81.30 mm。

從豎直方向應(yīng)力變化曲線可以看出，地表監(jiān)測(cè)點(diǎn)位0 m時(shí)的應(yīng)力值和堆載高度沒(méi)有關(guān)系，這是地表面巖土層受采空區(qū)變形的影響較小，其應(yīng)力值主要受自重應(yīng)力影響。在采空區(qū)頂板處的應(yīng)力值為正值，由于砂巖層在下沉彎曲過(guò)程中起到一個(gè)“梁模型”的作用，致使“梁模型”的下部受到拉應(yīng)力，且拉應(yīng)力值隨堆載高度增加逐漸增大。豎直方向應(yīng)力最大值并不是出現(xiàn)在采空區(qū)頂板表面，而受堆載高度影響。堆載高度由0 m增加到80 m的過(guò)程中，出現(xiàn)最大應(yīng)力值的觀測(cè)點(diǎn)分別為：距離地表40、60、70、80、100 m，在巖層中表現(xiàn)為相對(duì)地表往巖層上部移動(dòng)，這說(shuō)明上覆不均勻載荷加大，巖層受到破壞的范圍變大。隨著堆載高度增加，豎直方向應(yīng)力呈現(xiàn)壓應(yīng)力隨埋深增加到最大值后遞減為0后變?yōu)槔瓚?yīng)力，所以在露天礦排土場(chǎng)下覆采空區(qū)巖層移動(dòng)模型中必然曾在1個(gè)應(yīng)力為0的面。通過(guò)分析排土場(chǎng)連續(xù)堆載下采空區(qū)應(yīng)力分布特征，可以看出隨堆載高增加，在采空區(qū)的周?chē)艿綌D壓呈現(xiàn)拉應(yīng)力且影響范圍越大，應(yīng)力集中越明顯。如排土場(chǎng)連續(xù)堆載下采空區(qū)應(yīng)力分布特征如圖5。

圖5 排土場(chǎng)連續(xù)堆載下采空區(qū)應(yīng)力分布特征Fig.5 Stress distribution characteristics of goaf under continuous loading of dump

2.2 排土場(chǎng)下的不同位置采空區(qū)的沉陷規(guī)律

建立排土場(chǎng)堆載高度為60 m，采空區(qū)位置1到位置5分別位于排土場(chǎng)正下方左側(cè)300、225、150、75、0 m的模型，排土場(chǎng)下的不同位置采空區(qū)模型如圖6。

圖6 排土場(chǎng)下的不同位置采空區(qū)模型Fig.6 Goaf models at different locations under the dump

分別進(jìn)行FLAC3D數(shù)值模擬，并在采空區(qū)正上方距離地表0、20、40、60、80、100、120、140 m處布置監(jiān)測(cè)點(diǎn)，得到不同位置的采空區(qū)變形位移、應(yīng)力曲線，對(duì)比得到排土場(chǎng)下的不同位置采空區(qū)豎直方向上觀測(cè)線下沉規(guī)律。排土場(chǎng)下不同位置的采空區(qū)豎直方向上的地表沉降量曲線如圖7，排土場(chǎng)下不同位置的采空區(qū)豎直方向上的應(yīng)力變化曲線如圖8。

圖7 排土場(chǎng)下不同位置的豎直方向上的地表沉降量曲線Fig.7 Surface settlement curves in vertical direction of goaf at different positions under waste dump

圖8 排土場(chǎng)下不同位置的豎直方向上的應(yīng)力變化曲線Fig.8 Stress change curves in the vertical direction of goaf at different positions under waste dump

從圖7的地表沉降位移曲線可以看出，采空區(qū)位置從位置1到位置5過(guò)程中，地表沉降量由24.56 mm增大到34.10 mm，采空區(qū)位置1與未堆載排棄物下的地表沉降量幾乎差不多，采空區(qū)位置2地表沉降量達(dá)到28.00 mm，說(shuō)明采空區(qū)水平距離排土場(chǎng)300 m以上，排土場(chǎng)載荷對(duì)采空區(qū)變形影響較小。當(dāng)采空區(qū)距離排土場(chǎng)越接近，采空區(qū)正上方的下沉量越大，且靠近采空區(qū)的下部巖層的下沉量也越大。隨著采空區(qū)位置水平距離排土場(chǎng)越近，采空區(qū)上覆不均勻載荷隨之增大，地表沉降量逐漸增大，其中采空區(qū)頂板的下沉累計(jì)量由位置1的38.25 mm增加到堆載高度為60 m的70.80 mm。從圖8豎直方向應(yīng)力變化曲線可以看出，采空區(qū)位置從位置1到位置5過(guò)程中，排土場(chǎng)載荷產(chǎn)生的側(cè)向應(yīng)力引起采空區(qū)上方的應(yīng)力最大值增大，其中位置1到位置3過(guò)程中，采空區(qū)最頂部并未堆載排棄物，此時(shí)3個(gè)位置的垂直應(yīng)力出現(xiàn)在采空區(qū)上方巖層的中部位置，且地表沉降量約為采空區(qū)頂板累計(jì)沉降量的2倍，這與黃樂(lè)亭[22]得出的結(jié)論一致。

3 不同排土場(chǎng)形狀對(duì)采空區(qū)地表沉陷規(guī)律的影響

剝離物在界外沿幫排棄，隨著堆載工程量的不斷增大，外排土場(chǎng)高度逐漸增加，形成高陡外排土場(chǎng)，常見(jiàn)的外排土場(chǎng)有方形基礎(chǔ)和圓形基礎(chǔ)2種形式[23]。根據(jù)這2種基礎(chǔ)分別建立方形和圓形排土場(chǎng)堆載下的采空區(qū)三維模型，2個(gè)模型排土場(chǎng)的高度為50 m，排棄物的體積一樣（即保證2個(gè)排土場(chǎng)的自重一樣大），分別進(jìn)行FLAC3D數(shù)值模擬，并在采空區(qū)正上方距離地表0、20、40、60、80、100、130 m處布置監(jiān)測(cè)點(diǎn)，得到不同位置的位移、應(yīng)力曲線，對(duì)比得到方形和圓形排土場(chǎng)堆載下的采空區(qū)豎直方向上觀測(cè)線下沉規(guī)律。不同排土場(chǎng)形狀下豎直方向上的地表沉降量曲線如圖9，不同排土場(chǎng)形狀下豎直方向上的應(yīng)力變化曲線如圖10。

圖9 不同排土場(chǎng)形狀下豎直方向上的地表沉降量曲線Fig.9 Surface settlement curves in the vertical direction under different dump shapes

圖10 不同排土場(chǎng)形狀下豎直方向上的應(yīng)力變化曲線Fig.10 Stress change curves in the vertical direction under different dump shapes

從地表沉降位移曲線可以看出，在采空區(qū)上覆載荷一樣大小，但方形排土場(chǎng)和圓形排土場(chǎng)的幾何結(jié)構(gòu)差異產(chǎn)生的不均勻應(yīng)力不一致，導(dǎo)致方形排土場(chǎng)的地表下沉位移量比圓形排土場(chǎng)的大。2種排土場(chǎng)結(jié)構(gòu)影響的地表面下沉為幾乎一樣大小且整體的下沉速度較為均勻。從豎直方向應(yīng)力變化曲線可以看出，2種排土場(chǎng)結(jié)構(gòu)對(duì)采空區(qū)上覆巖層的垂直應(yīng)力影響較小且兩者幾乎接近，也呈現(xiàn)出豎直方向應(yīng)力的壓應(yīng)力隨埋深增加到最大值后遞減為0后變?yōu)槔瓚?yīng)力。

為了更好揭示排土場(chǎng)結(jié)構(gòu)對(duì)采空區(qū)變形的影響，在采空區(qū)長(zhǎng)度上頂板處每間隔50 m布置1個(gè)位移和應(yīng)力監(jiān)測(cè)點(diǎn)，根據(jù)模擬分析得到監(jiān)測(cè)數(shù)據(jù)曲線，不同排土場(chǎng)形狀下采空區(qū)長(zhǎng)度上頂板處下沉曲線如圖11，不同排土場(chǎng)形狀下采空區(qū)長(zhǎng)度上頂板處應(yīng)力變化曲線如圖12。

圖11 不同排土場(chǎng)形狀下采空區(qū)長(zhǎng)度上頂板處下沉曲線Fig.11 Subsidence curves at upper roof of goaf length under different dump shapes

圖12 不同排土場(chǎng)形狀下采空區(qū)長(zhǎng)度上頂板處應(yīng)力變化曲線Fig.12 Stress change curves at upper roof of goaf length under different dump shapes

由圖11可以看出，2種排土場(chǎng)結(jié)構(gòu)下的頂板下沉量都在采空區(qū)長(zhǎng)度上頂板的中間位置達(dá)到最大，且下沉位移曲線呈V型對(duì)稱(chēng)。方形排土場(chǎng)對(duì)采空區(qū)頂板的下沉位移量比圓形排土場(chǎng)的大，且相同埋深位置的2種結(jié)構(gòu)下的位移差值逐漸增加，在采空區(qū)長(zhǎng)度上頂板的中間位置達(dá)到最大。由圖12可以看出，2種排土場(chǎng)結(jié)構(gòu)下的頂板應(yīng)力都為拉應(yīng)力，且應(yīng)力曲線呈M型對(duì)稱(chēng)。在采空區(qū)長(zhǎng)度上頂板中間位置，圓形排土場(chǎng)造成的采空區(qū)頂板應(yīng)力比方形排土場(chǎng)的小。

分析上述地表沉降位移曲線和應(yīng)力曲線圖可知，在露天礦剝離量相同的情況下，選擇圓形排土場(chǎng)對(duì)下覆巖層中采空區(qū)變形影響更小，更有利于排土場(chǎng)邊坡的穩(wěn)定性。

4 結(jié) 語(yǔ)

1）隨著排土場(chǎng)連續(xù)堆載高度加大，采空區(qū)上覆不均勻載荷隨之增大，地表沉降量逐漸增大，豎直方向應(yīng)力呈現(xiàn)壓應(yīng)力隨埋深增加到最大值后遞減為0后變?yōu)槔瓚?yīng)力，在露天礦排土場(chǎng)下伏采空區(qū)的上覆巖層中存在1個(gè)應(yīng)力為零的面。

2）隨著采空區(qū)位置水平距離排土場(chǎng)越近，地表沉降量逐漸增大，采空區(qū)水平距離排土場(chǎng)300 m以上，排土場(chǎng)載荷對(duì)采空區(qū)變形影響較小，采空區(qū)最頂部并未堆載排棄物，垂直應(yīng)力出現(xiàn)在采空區(qū)上方巖層的中部位置，且地表沉降量約為采空區(qū)頂板累計(jì)沉降量的2倍。

3）方形排土場(chǎng)的地表下沉位移量比圓形排土場(chǎng)的大，2種排土場(chǎng)結(jié)構(gòu)下的頂板下沉量都在采空區(qū)長(zhǎng)度上頂板的中間位置達(dá)到最大，且下沉位移曲線呈V型對(duì)稱(chēng)，應(yīng)力位移曲線呈M型。圓形排土場(chǎng)對(duì)下覆巖層中采空區(qū)變形影響更小，更有利于露天礦排土場(chǎng)的穩(wěn)定性。