張博輝，賀衛(wèi)中，呂 新

(1.陜西省煤炭科學(xué)研究所，陜西 西安710001;2.陜西省地質(zhì)環(huán)境監(jiān)測(cè)總站，陜西 西安710054;3.長(zhǎng)安大學(xué) 環(huán)境科學(xué)與工程學(xué)院，陜西西安710054)

陜北是我國(guó)西部最重要的煤炭生產(chǎn)基地，采煤與生態(tài)環(huán)境保護(hù)矛盾突出，安全、科學(xué)開采不僅關(guān)系到國(guó)家能源供給，也關(guān)系到區(qū)域水資源和生態(tài)安全[1－3]。由于地下煤炭資源的大規(guī)模開采，破壞了煤層上覆巖層的穩(wěn)定，導(dǎo)致上腹巖層的破碎、垮落，進(jìn)而造成地表的下沉、塌陷和開裂，導(dǎo)致泉井干涸，河流斷流[4]，所以煤炭開采造成的地表移動(dòng)變形破壞是煤礦安全生產(chǎn)面臨的嚴(yán)峻問(wèn)題。如何采用合適的方法對(duì)煤炭開采引起的地表移動(dòng)進(jìn)行準(zhǔn)確預(yù)測(cè)，從而采取行之有效有的防治措施是目前亟待解決的主要問(wèn)題。

多位專家學(xué)者就采煤沉陷預(yù)測(cè)進(jìn)行了研究，1965年，劉寶琛和廖國(guó)華在影響函數(shù)法的基礎(chǔ)上建立了非常實(shí)用的概率積分法[5];1997年鄒友峰對(duì)地表下沉系數(shù)計(jì)算方法進(jìn)行了研究[6];2000年戴華陽(yáng)、王金莊等基于傾角的變化，建立了的矢量預(yù)計(jì)法[7，8];我國(guó)學(xué)著王泳嘉和杜慶華等對(duì)邊界元法進(jìn)行了大量研究，將邊界元法引入巖層移動(dòng)研究，對(duì)均質(zhì)模型進(jìn)行了探討[9]。

杭來(lái)灣井田位于陜北風(fēng)沙灘地區(qū)，采用長(zhǎng)臂式全部垮落法開采方式，高強(qiáng)度開采在地表產(chǎn)生了塌陷和地裂縫等地質(zhì)災(zāi)害[10]?，其地表沉陷的分布規(guī)律和沉陷機(jī)理在陜北風(fēng)沙灘地區(qū)具有典型性。本論文以杭來(lái)灣井田為例，以采煤沉陷基本理論為指導(dǎo)，將野外調(diào)查與室內(nèi)綜合分析相結(jié)合，運(yùn)用概率積分法對(duì)該井田的30101盤區(qū)采煤引起的地表移動(dòng)變形進(jìn)行預(yù)測(cè)，以此為該井田其他盤區(qū)以及類似井田的開采沉陷造成的損害提供科學(xué)依據(jù)[11，12]。

1 研究區(qū)概況

杭來(lái)灣井田屬于榆神礦區(qū)的一期規(guī)劃區(qū)，位于規(guī)劃區(qū)的西南部，是榆神礦區(qū)一期規(guī)劃區(qū)建設(shè)的重點(diǎn)煤礦，井田東西長(zhǎng)約 12.0 km、南北寬約 7.7 km，面積 92 km2。

井田位于鄂爾多斯盆地，區(qū)內(nèi)地質(zhì)構(gòu)造簡(jiǎn)單，無(wú)較大褶皺和斷裂，無(wú)明顯巖漿活動(dòng)痕跡，局部發(fā)育有寬緩的波狀起伏。井田內(nèi)的地層由老至新依次為:三疊系上統(tǒng)瓦窯堡組(T3w)、侏羅系下統(tǒng)富縣組(J1f)、侏羅系中統(tǒng)延安組(J2y)、侏羅系中統(tǒng)直羅組(J2z)、第四系下更新統(tǒng)午城組(Q1w)、中更新統(tǒng)離石組(Q2l)、上更新統(tǒng)薩拉烏蘇組(Q31s)、全新統(tǒng)沖洪積層(Q41al+pl、Q42al)、風(fēng)積層(Q42eol)。

礦井首采區(qū)301盤區(qū)的30101綜采工作面，位于杭來(lái)灣井田的東北部。設(shè)計(jì)礦井初期開采3號(hào)煤層，該煤層最大厚度 10.66 m，最小厚度 8.27 m，平均 9.13 m，煤層傾角 0.5°，為近水平煤層，構(gòu)造簡(jiǎn)單。開采方式為傾斜長(zhǎng)臂式全部垮落法開采。

2 采煤移動(dòng)變形分布規(guī)律

2.1 地表移動(dòng)變形規(guī)律

根據(jù)野外調(diào)查，30101工作面開采后，在其上方形成了一個(gè)大范圍移動(dòng)盆地，面積約為2.2 km2。陜西有色榆林煤業(yè)有限公司的地表移動(dòng)變形監(jiān)測(cè)點(diǎn)監(jiān)測(cè)到地表下沉最大值為7.25 m，傾斜變形最大值為69.1 mm/m，水平變形最大值為31.5 mm/m，主要影響半徑為91～105 m，地表最大下沉速度為 773.3 mm/d，地表移動(dòng)持續(xù)時(shí)間為 547.5 d。

在采空區(qū)邊界上方地表附近形成地裂縫，其發(fā)育程度與距開采邊界的距離以及開采深度有關(guān)。對(duì)本次調(diào)查裂縫的研究區(qū)域選取在一個(gè)主斷面上，其裂縫發(fā)育明顯，并且該區(qū)域的地裂縫在長(zhǎng)時(shí)間停采后已經(jīng)逐漸趨于穩(wěn)定。調(diào)查測(cè)量區(qū)地表上的裂縫共發(fā)育21條，裂縫走向約135°，與煤炭開采邊界走向132°大致平行;裂縫寬度不等;區(qū)域首尾裂縫寬度小(基本小于1 cm)，裂縫兩側(cè)基本無(wú)高差，區(qū)域中間的裂縫寬度較大且兩側(cè)存在高差，采空區(qū)中心的一側(cè)的裂縫低于另一側(cè)的裂縫。根據(jù)地裂縫的調(diào)查結(jié)果及該區(qū)地層資料，經(jīng)過(guò)研究計(jì)算求得產(chǎn)生地裂縫的裂縫角為89°。

2.2 采煤沉陷形成原因

當(dāng)煤層采空以后，采空區(qū)上部的覆巖及采空形成的煤柱邊邦，均形成自由面。原來(lái)的應(yīng)力平衡被破壞，在上覆巖土層重力作用下，覆巖承受的壓力隨著采空區(qū)范圍的擴(kuò)大而增加，當(dāng)這種壓力超過(guò)煤層頂板巖石的承載力以后，頂板巖石破裂塌落形成冒落帶，進(jìn)而在其上部依次形成裂隙帶和彎曲變形帶。在彎曲帶影響下地面上造成塌陷，形成塌陷盆地或地裂縫。

影響地表沉陷的因素主要有兩個(gè):地質(zhì)條件和開采條件。地質(zhì)條件主要包括上覆巖層的巖性、巖層層位、表土層性質(zhì)和煤層傾角等;開采條件主要包括采深采厚比、開采方法以及工作面的推進(jìn)速度等。結(jié)合30101工作面的地質(zhì)條件以及開采方式對(duì)采煤沉陷進(jìn)行分析(見圖1)。

圖1 杭來(lái)灣井田采煤沉陷示意圖

3 移動(dòng)變形預(yù)測(cè)

3.1 相關(guān)參數(shù)計(jì)算

概率積分法是我國(guó)的劉寶琛和廖國(guó)華等多位學(xué)者在影響函數(shù)法和隨機(jī)介質(zhì)理論的基礎(chǔ)上提出的一種比較實(shí)用的開采沉陷預(yù)計(jì)方法[3]。其基本思想是將巖體看作為隨機(jī)介質(zhì)，把地表的移動(dòng)變形看作隨機(jī)事件且服從統(tǒng)計(jì)規(guī)律，把整個(gè)區(qū)域分解成無(wú)限個(gè)微小單元，根據(jù)疊加原理看作開采對(duì)地表的影響等于采區(qū)內(nèi)開采各個(gè)微小單元對(duì)地表的影響之和來(lái)計(jì)算整個(gè)開采引起的地表及巖層變形。利用概率積分法對(duì)地表沉陷進(jìn)行預(yù)計(jì)時(shí)主要用到的參數(shù)有:下沉系數(shù)、水平移動(dòng)系數(shù)、主要影響角正切、開采傳播角及拐點(diǎn)偏移距。李培現(xiàn)等通過(guò)綜合分析國(guó)內(nèi)礦區(qū)的實(shí)測(cè)數(shù)據(jù)，并進(jìn)行回歸計(jì)算，建立了從礦區(qū)地質(zhì)條件確定參數(shù)的統(tǒng)計(jì)方法 ，通過(guò)對(duì)比，表明用該方法計(jì)算出的參數(shù)與實(shí)測(cè)參數(shù)誤差較小。為缺少地表移動(dòng)實(shí)測(cè)資料的情況下獲取概率積分法的參數(shù)提供了技術(shù)保障。

表1 各個(gè)參數(shù)的計(jì)算結(jié)果

3.2 移動(dòng)變形預(yù)測(cè)

根據(jù)該礦區(qū)的地質(zhì)結(jié)構(gòu)、巖性和采空區(qū)平面分布圖，對(duì)30101綜采工作面可能出現(xiàn)的沉陷進(jìn)行了預(yù)計(jì)，應(yīng)用概率積分法對(duì)30101工作面(300 m×4 252 m)進(jìn)行地表移動(dòng)預(yù)測(cè)，并使用MATLAB軟件繪制移動(dòng)變形的預(yù)測(cè)結(jié)果，結(jié)果見圖2和圖3。

由圖2、3，得出的理論下沉值是關(guān)于采空區(qū)中心軸對(duì)稱的橢圓形，沉陷最大值在采空區(qū)中心，向兩側(cè)逐漸減小，在采空區(qū)傾向邊界外100 m處塌陷值小于1 cm，視為沉陷邊界，經(jīng)計(jì)算分析，該工作面的理論邊界角為68.9°。下沉等值線以采空區(qū)中心為對(duì)稱原點(diǎn)呈橢圓分布，并且橢圓的長(zhǎng)軸位于開采工作面中心;由圖4，看出下沉曲線近似一條拋物線，在該工作面地表塌陷的理論最大值為7.17 m，理論值非常接近陜西有色榆林煤業(yè)有限公司監(jiān)測(cè)的實(shí)測(cè)最大下沉值7.25 m;理論主要影響半徑為104.42 m，該理論值在實(shí)測(cè)的91～105 m區(qū)間之內(nèi)。

圖2 開采沉陷下沉值等值線圖(m)

圖3 開采沉陷下沉值三維斜視圖

圖4 走向主斷面的下沉曲線

通過(guò)分析下沉曲線，最大下沉值位于平行于開采傾向的采空區(qū)中心區(qū)域(EOF段);在采空區(qū)傾向邊界以外的100 m處(A、B)塌陷值已小于10 mm，可看作塌陷邊界。所以下沉曲線的分布規(guī)律為:在采空區(qū)中央正上方的平地區(qū)域處地表下沉值最大，從平底盆地兩側(cè)向沉陷區(qū)邊緣逐漸減小，在移動(dòng)盆地邊界A、B處下沉值十分接近0，而且下沉曲線以采空區(qū)中心線對(duì)稱;地表移動(dòng)盆地的最大下沉值已達(dá)到該地質(zhì)采礦條件下的最大值;拐點(diǎn)處(C、D)的下沉值為最大下沉值的一半。

從圖中分析地表裂縫的分布規(guī)律，在近似平底區(qū)域EOF內(nèi)的水平變形、傾斜、曲率為0(或接近0)，此時(shí)該區(qū)域內(nèi)沒(méi)有水平方向的力，所以此區(qū)域內(nèi)的煤礦開采基本不會(huì)發(fā)生引發(fā)地裂縫;在EC和FD之間的區(qū)域該區(qū)域主要是不均勻下沉，地表的曲率為負(fù)而且水平變形也是負(fù)的，所以該區(qū)域的地表主要受壓縮而不會(huì)形成拉裂縫，裂縫主要為不均勻沉陷產(chǎn)生的裂縫;而在AC、BD段地表曲率和水平變形均為正，所以該區(qū)域的地表受拉力作用產(chǎn)生拉裂縫。

通過(guò)概率積分法計(jì)算的地表移動(dòng)變形值見表2，傾斜變形理論最大值為68.9 mm/m，實(shí)測(cè)的最大值為69.1 mm/m;水平變形理論最大值為31.38 mm/m，實(shí)測(cè)最大值為31.5 mm/m。該結(jié)果再一次證明了概率積分法在該井田的適用性與準(zhǔn)確性。

表2 走向主斷面內(nèi)各點(diǎn)的移動(dòng)變形計(jì)算值

我國(guó)一般采用針對(duì)一般磚木結(jié)構(gòu)建筑物的標(biāo)準(zhǔn):i=3 mm/m、ε=2 mm/m、k=0.2 mm/m2，是以這三個(gè)標(biāo)準(zhǔn)的最外一個(gè)值來(lái)確定危險(xiǎn)邊界，根據(jù)表2的計(jì)算結(jié)果，可知在ε=2 mm/m處對(duì)應(yīng)的邊界為最外邊界，此時(shí)距采空區(qū)中心的距離為243 m，即取距采空區(qū)邊界93 m處為危險(xiǎn)邊界，移動(dòng)角為70.4°，同理，確定距開采邊界93 m以外為對(duì)生態(tài)環(huán)境的輕微影響區(qū)，該區(qū)沉陷盆地平緩、在地表沉陷區(qū)范圍不會(huì)出現(xiàn)拉伸裂縫，對(duì)耕地影響很小，對(duì)生態(tài)環(huán)境影響輕微，可自然恢復(fù)。裂縫角對(duì)應(yīng)的邊界是裂縫邊界，移動(dòng)角對(duì)應(yīng)的邊界是移動(dòng)邊界，邊界角對(duì)應(yīng)的地表區(qū)域是移動(dòng)盆地的沉陷邊界。而且沉陷由于地質(zhì)采礦條件簡(jiǎn)單，以及根據(jù)移動(dòng)盆地的對(duì)稱性和實(shí)測(cè)地裂縫的分布，推測(cè)出采區(qū)內(nèi)地裂縫的分布范圍。經(jīng)過(guò)預(yù)測(cè)計(jì)算，移動(dòng)盆地的沉陷范圍面積約為2.193 km2。

(3)通過(guò)野外實(shí)測(cè)數(shù)據(jù)證明運(yùn)用概率積分法對(duì)杭來(lái)灣井田30101綜采工作面地表移動(dòng)變形的預(yù)測(cè)結(jié)果精度高，對(duì)地表移動(dòng)變形的評(píng)價(jià)迅速準(zhǔn)確，為杭來(lái)灣井田及其他類似井田的地表移動(dòng)變形預(yù)測(cè)提供了科學(xué)依據(jù)。

4 結(jié)語(yǔ)

以采煤沉陷基本理論為指導(dǎo)，將野外調(diào)查與室內(nèi)綜合分析相結(jié)合，探討采煤沉陷分布規(guī)律。以杭來(lái)灣井田30101工作面為例，運(yùn)用概率積分法對(duì)采煤引起的地表移動(dòng)變形進(jìn)行了預(yù)測(cè)。

(1)通過(guò)野外實(shí)測(cè)以及概率積分法的預(yù)測(cè)分析結(jié)果，得出煤炭開采引起地表變形的理論邊界角為68.9°，移動(dòng)角為70.4°，裂縫角為 89°，移動(dòng)盆地的沉陷范圍面積是 2.193 km2。

(2)通過(guò)概率積分法計(jì)算的地表移動(dòng)變形值，傾斜變形理論最大值為68.9 mm/m，實(shí)測(cè)的最大值為69.1 mm/m;水平變形理論最大值為31.38 mm/m，實(shí)測(cè)最大值為31.5 mm/m。該結(jié)果證明概率積分法在該井田的適用性與準(zhǔn)確性。

[1]王雙明，黃慶享，范立民，等.生態(tài)脆弱礦區(qū)含(隔)水層特征及保水開采分區(qū)研究[J].煤炭學(xué)報(bào).2010，35(1):7－14.

[2]范立民.論保水采煤?jiǎn)栴}[J].煤田地質(zhì)與勘探.2005，33(5):50－53.

[3]范立民.陜北地區(qū)采煤造成的地下水滲漏及其防治對(duì)策分析[J].礦業(yè)安全與環(huán)保.2007，34(5):62 －64.

[4]馬雄德，王文科，范立民，等.生態(tài)脆弱礦區(qū)采煤對(duì)泉的影響[J].中國(guó)煤炭地質(zhì).2010，22(1):32 －36.

[5]劉寶琛，廖國(guó)華.煤礦地表移動(dòng)的基本規(guī)律[M].北京:中國(guó)工業(yè)出版社.1965.

[6]鄒友峰.地表下沉系數(shù)計(jì)算方法研究[J].巖土工程學(xué)報(bào).1997，19(3):109－112.

[7]戴華陽(yáng).基于傾角變化的開采沉陷模型及其GIS可視化應(yīng)用研究[J].巖石力學(xué)與工程學(xué)報(bào).2002，21(1):148－148.

[8]戴華陽(yáng)，王金莊，蔡美峰.巖層與地表移動(dòng)的矢量預(yù)計(jì)法[J].煤炭學(xué)報(bào).2002，27(5):473 －478.

[9]王泳嘉.離散元法及其在巖石力學(xué)中的應(yīng)用[J].金屬礦山.1986，(8):13 －17.

[10]范立民.榆神府區(qū)煤炭開采強(qiáng)度與地質(zhì)災(zāi)害研究[J].中國(guó)煤炭.2014，40(5):52 －55.

[11]王雙明，范立民，楊宏科.陜北煤炭資源可持續(xù)發(fā)展之開發(fā)思路[J].中國(guó)煤田地質(zhì).2003，15(5):6 ～8，11.

[12]朱占榮，陳利平，郝延錦.特厚砂土層條件下的開采沉陷規(guī)律研究[J].中國(guó)煤炭地質(zhì).2013，25(8);42－45.