胡振琪，王新靜，賀安民

(1.中國礦業(yè)大學(xué)(北京) 土地復(fù)墾與生態(tài)重建研究所,北京 100083；2.神華神東煤炭集團(tuán)有限責(zé)任公司,陜西 神木 719315)

風(fēng)積沙區(qū)采煤沉陷地裂縫分布特征與發(fā)生發(fā)育規(guī)律

胡振琪1，王新靜1，賀安民2

(1.中國礦業(yè)大學(xué)(北京) 土地復(fù)墾與生態(tài)重建研究所,北京 100083；2.神華神東煤炭集團(tuán)有限責(zé)任公司,陜西 神木 719315)

煤炭開發(fā)可能加劇對生態(tài)脆弱區(qū)環(huán)境的損傷，風(fēng)積沙區(qū)地下煤炭開采對地表環(huán)境的影響主要是沉陷地裂縫，且目前的研究缺乏對地裂縫從發(fā)育到湮滅全過程的研究，筆者通過井上下相結(jié)合的空間坐標(biāo)控制體系和自主研發(fā)的動態(tài)地裂縫監(jiān)測方法，對補(bǔ)連塔12406綜采工作面地裂縫進(jìn)行持續(xù)動態(tài)監(jiān)測，提出了邊緣裂縫的分布規(guī)律和動態(tài)裂縫的發(fā)生發(fā)育規(guī)律及其與地質(zhì)采礦條件之間的關(guān)系。研究表明：① 動態(tài)裂縫超前于當(dāng)前工作面向前發(fā)展，裂縫超前距與工作面日進(jìn)尺量呈現(xiàn)明顯的線性正相關(guān)，平均的裂縫超前距為10.359m，超前裂縫角為近似垂直角；② 高強(qiáng)度開采動態(tài)地裂縫的兩側(cè)無明顯落差，裂縫的寬度值在采動過程中有顯著的周期性變化，呈“M”型雙峰波形，且第1個峰值明顯大于第2個峰值，峰值比為1.4～3.8；裂縫發(fā)育周期T包含兩個時長近似相等的“開裂—閉合”過程，約為18d,在此基礎(chǔ)上，結(jié)合地質(zhì)采礦條件，建立了T的通用函數(shù)模型，表明動態(tài)地裂縫具有快速閉合的自修復(fù)特征，不需要人工修復(fù)；③ 邊緣裂縫以“帶狀”、“O”形圈的形態(tài)分布在工作面開采邊界的內(nèi)側(cè)，裂縫帶寬為46～50m，裂縫帶整體向內(nèi)收縮，臨近工作面采動會減輕原有地裂縫的影響，減輕約40%，開采結(jié)束后邊緣裂縫仍然存在，是主要地表環(huán)境損傷區(qū)和重點(diǎn)人工修復(fù)區(qū)域。

開采沉陷；地裂縫；生態(tài)損傷；風(fēng)積沙區(qū)

地裂縫是煤炭開采后地表破壞的形式之一，也是風(fēng)積沙區(qū)采煤沉陷對地表環(huán)境損傷最直觀的表現(xiàn)形式，其發(fā)生以及發(fā)育規(guī)律與煤礦的開發(fā)建設(shè)活動密切相關(guān)，受到社會各界的廣泛關(guān)注[1]。目前，針對該地區(qū)的采煤沉陷對土地生態(tài)環(huán)境的影響規(guī)律及程度，國內(nèi)學(xué)者開展了諸多研究，形成了兩種意見相悖的觀點(diǎn)，相當(dāng)一部分學(xué)者認(rèn)為采煤沉陷對風(fēng)積沙區(qū)土地的生態(tài)環(huán)境造成了明顯的影響[2-6]，且部分指標(biāo)(如土壤含水量)受影響的周期達(dá)2a以上；也有人認(rèn)為采煤沉陷對該地區(qū)土地生態(tài)環(huán)境的影響很小或者沒有影響[7-9]。出現(xiàn)這種差異的原因在于：① 研究的時間節(jié)點(diǎn)不同，風(fēng)積沙區(qū)土地受損有特定的演變過程，例如張口地裂縫，隨采動過程中呈現(xiàn)“開裂、閉合”的現(xiàn)象，其發(fā)育規(guī)律必然會對上述研究結(jié)果產(chǎn)生較大的影響；② 調(diào)查取樣的方法不同，一些研究的調(diào)查樣分布在下沉盆地內(nèi)部，而其他研究則分布在下沉盆地的邊界，有大量的地裂縫分布，不同調(diào)查區(qū)對采煤沉陷的響應(yīng)特征以及程度必然存在差異。因此，開展該地區(qū)采煤地裂縫分布特征及其發(fā)生發(fā)育規(guī)律的研究，是研究煤炭開采對地表環(huán)境影響的前提和必要條件。

對于采動引起的地裂縫的分布規(guī)律與特征，眾多學(xué)者采用現(xiàn)場實(shí)測、理論分析等方面進(jìn)行諸多研究，表明采動裂縫分布形態(tài)隨工作面的推進(jìn)而發(fā)生變化。欒長青等[10]認(rèn)為地裂縫的空間分布與煤炭開采的方位有關(guān)，地裂縫的發(fā)展與開采進(jìn)度也有密切聯(lián)系，李彥軍[11]較為詳細(xì)地敘述了靖遠(yuǎn)礦區(qū)地裂縫的發(fā)育狀態(tài)，并對其進(jìn)行穩(wěn)定性分析，杜善周[12]指出神東某工作面上方動態(tài)裂縫以圓錐狀超前發(fā)展，超前裂縫角呈高角度，李曉等[13]通過現(xiàn)場實(shí)測發(fā)現(xiàn)金川二礦地裂縫兩側(cè)土體的位移具有三維特征。吳侃等[14-15]則利用土體破壞準(zhǔn)則，建立了適宜東部地區(qū)條件下開采地裂縫分布規(guī)律以及半定量的分析方法，部分研究人員利用三維激光掃描[16]、InSAR/D-InSAR[17-18]等手段獲取到采動地裂縫的分布范圍。上述的研究成果多集中在中東部等礦區(qū)，主要針對在采動過后地裂縫的分布形態(tài)，相對而言，對于風(fēng)積沙區(qū)地裂縫的分布特征以及發(fā)生發(fā)育動態(tài)變化規(guī)律的研究相對較少。

根據(jù)開采沉陷理論與方法，采煤地裂縫主要區(qū)分為邊緣裂縫和動態(tài)裂縫。邊緣裂縫一般在開采工作面的外邊緣區(qū)，動態(tài)裂縫位于工作面上方地表，平行于工作面，并隨著工作面的推進(jìn)不斷產(chǎn)生和閉合。但對于風(fēng)積沙區(qū)，特別在高強(qiáng)度開采條件下，動態(tài)裂縫的發(fā)生發(fā)育規(guī)律以及它與開采過程之間的關(guān)系和定量描述、邊緣裂縫的發(fā)生發(fā)育特征以及裂縫深度與寬度的關(guān)系等科學(xué)問題，尚沒有較為系統(tǒng)的報道，需要得到解決。過去的研究大都是工作面終采后進(jìn)行裂縫的調(diào)查，缺乏裂縫發(fā)生發(fā)育全過程的監(jiān)測，因而，對上述問題始終無法很好地回答。

筆者通過建立井上下相結(jié)合的空間坐標(biāo)控制體系，對補(bǔ)連塔12406綜采工作面的地裂縫在采動過程中的演變特征進(jìn)行持續(xù)的動態(tài)監(jiān)測，提出了地裂縫的分布特征與發(fā)生發(fā)育規(guī)律，并在此基礎(chǔ)上，結(jié)合采礦地質(zhì)條件，建立了動態(tài)地裂縫發(fā)育周期的定量函數(shù)模型，為該區(qū)域煤炭開采引起土地生態(tài)環(huán)境損傷研究提供時間位置和時間基準(zhǔn)等技術(shù)參數(shù)。

1 研究區(qū)概況和監(jiān)測方法

1.1 研究區(qū)概況

12406綜采工作面隸屬神東礦區(qū)補(bǔ)連塔煤礦，地處毛烏素沙漠東北邊緣，下山方向有老采空區(qū)分布，采用長臂開采、垮落式管理頂板方式，2011年4月開始回采，2011年12月完成全部回采工作。工作面長度300.5m，走向長度3 592m，日平均開采速度約為12m/d。工作面開采1-2煤，平均煤厚4.81m，煤層埋深190～220m，煤層傾角1°～3°，屬近水平煤層，上覆基巖的厚度160～200m，其中基本頂平均厚度181m，以粉砂巖為主，夾有砂質(zhì)泥巖薄層，直接頂平均厚度5m，水平層理，中間夾薄層細(xì)砂巖，偽頂平均厚度4.8m，近地表上覆松散層厚度8～27 m。由于較硬覆巖的存在，煤炭開采導(dǎo)致地表下沉的最大值約為2500mm，下沉系數(shù)約為0.55，數(shù)值偏小，塌陷程度較輕，地裂縫為主要的地表破壞形式。

為了消除地形起伏對裂縫發(fā)育形態(tài)的附加影響，利用ArcGIS對研究區(qū)1∶5000的地形圖進(jìn)行坡度分析，選取距離開切眼300～500m、地勢相對平坦的區(qū)域作為動態(tài)裂縫的重點(diǎn)監(jiān)測區(qū)，而邊緣裂縫監(jiān)測區(qū)域則為從開切眼一直到距離開切眼500m處，待工作面終采后，完成最后觀測。

1.2 準(zhǔn)備工作

研究組在工作面開采之前，采用剖面法布設(shè)了地表移動觀測線，并將其視為裂縫觀測基準(zhǔn)線。兩條基準(zhǔn)線相互垂直正交，分布于下沉盆地的主斷面上。其中，走向觀測線、傾向觀測線布設(shè)長度分別為900m和1100m，點(diǎn)間距20～30m，共計99個地表移動觀測點(diǎn)，另外在采動影響范圍加設(shè)了8個控制點(diǎn)，觀測線布設(shè)方式如圖1所示。

圖1 裂縫觀測基準(zhǔn)線布設(shè)示意Fig.1 Layout of fissure observation reference line

1.3 地裂縫的監(jiān)測方法

裂縫監(jiān)測采用動態(tài)監(jiān)測法，即動態(tài)裂縫選擇新發(fā)現(xiàn)的裂縫進(jìn)行發(fā)育周期全過程的觀測，邊緣裂縫也是隨工作面推進(jìn)不斷地觀測，直至整個工作面結(jié)束后對邊緣裂縫進(jìn)行系統(tǒng)監(jiān)測。

1.3.1動態(tài)裂縫的監(jiān)測

(1)當(dāng)工作面推進(jìn)至監(jiān)測區(qū)域的臨界位置時，利用全站儀和鄰近的地表移動觀測站，測定工作面前方地表新出現(xiàn)的地裂縫，形成裂縫分布圖件；(2)將其中距離工作面推進(jìn)位置最遠(yuǎn)的裂縫視為最前端裂縫，并作為監(jiān)測對象，在裂縫線上選擇若干取樣點(diǎn)，量取裂縫寬度的初始距離，并在裂縫兩端布設(shè)若干成對的控制棒，用于進(jìn)行裂縫寬度變化監(jiān)測；(3)每隔1～2d量取裂縫兩端控制棒的距離，進(jìn)而求得裂縫在采動過程中的演變特征，直至裂縫完全閉合裂縫結(jié)束量取工作，監(jiān)測裂縫寬度的同時記錄觀測時間和采煤工作面推進(jìn)的位置；(4)重復(fù)上述步驟，在不同開采階段，得到一系列新裂縫發(fā)生位置與開采的相對關(guān)系以及多條裂縫發(fā)育周期的數(shù)據(jù)；(5)地表移動變形監(jiān)測也與之同步進(jìn)行，以獲取裂縫變化與沉陷過程的關(guān)系。

1.3.2邊緣裂縫的監(jiān)測

邊緣裂縫的監(jiān)測重點(diǎn)側(cè)重于分布范圍，記錄包括最外側(cè)、最內(nèi)側(cè)以及主裂縫(最大裂縫)的空間位置、寬度以及落差等相關(guān)屬性，并形成邊緣裂縫分布圖件和屬性表；選取相對穩(wěn)定的邊緣裂縫的若干點(diǎn)進(jìn)行深度監(jiān)測，監(jiān)測點(diǎn)與工作面推進(jìn)位置的距離應(yīng)大于200m，量取監(jiān)測點(diǎn)處裂縫的寬度，同時采用石膏漿進(jìn)行灌注，待膏體固結(jié)后，進(jìn)行開挖，量取裂縫的擴(kuò)展深度。

2 結(jié)果與討論

2.1 動態(tài)裂縫的分布特征與發(fā)生發(fā)育規(guī)律

通過對觀測區(qū)域內(nèi)每天新出現(xiàn)的地裂縫的空間位置進(jìn)行測量，總計觀測到55條動態(tài)裂縫，并對5條最前端裂縫(每次觀測時段內(nèi)出現(xiàn)的最前方的裂縫)的演變特征進(jìn)行了持續(xù)動態(tài)監(jiān)測，直至裂縫完全閉合為此。

2.1.1裂縫產(chǎn)生位置及其與采礦的關(guān)系

新出現(xiàn)的動態(tài)裂縫主要分布在工作面走向中心線附近，且平行于工作面。在每個觀測時段，裂縫以一定間距向前發(fā)生于上方地表，呈“帶狀”形態(tài)分布，中間分布若干斷裂的裂縫，且表征相對較弱。隨著工作面的向前推進(jìn)，前方地表不斷產(chǎn)生新的地裂縫，先前出現(xiàn)的裂縫開始擴(kuò)展，裂縫的寬度和長度不斷增加，并在采動影響減弱后，迅速閉合。

通過井上下相結(jié)合的坐標(biāo)控制系統(tǒng)，結(jié)合工作面生產(chǎn)進(jìn)度表以及對應(yīng)觀測時段的工作面位置，最前端地裂縫產(chǎn)生位置與采礦過程的相關(guān)關(guān)系見表1。

最前端裂縫均超前于工作面推進(jìn)距離，且動態(tài)裂縫的超前距(Sc)隨著工作面日進(jìn)尺量(v，也稱推進(jìn)速度)的增大而增大，相關(guān)性分析顯示，二者呈現(xiàn)明顯的線性正相關(guān)，線性回歸模型為Sc=0.637v+2.715，R2=0.862。以工作面平均日進(jìn)尺量12m計，動態(tài)裂縫的發(fā)生位置平均超前工作面10.359m，則超前裂縫角為：δ=87.035°，呈近似垂直角的形態(tài)分布。這也與杜善周等[12]的研究結(jié)論一致。地裂縫的產(chǎn)生與采動巖體破壞、應(yīng)力分布特征密切相關(guān)，高強(qiáng)度綜放開采使覆巖拉伸富集區(qū)位于煤壁的前后方[19]，采動巖體裂隙場以高角度甚至垂直巖層層面的裂隙為主[20]，由于12406風(fēng)積沙松散層的存在，其抗剪能力很差，動態(tài)地裂縫的發(fā)生位置往往反映了覆巖變形以及裂隙特征。

表1工作面開采進(jìn)度和最前端裂縫位置對應(yīng)關(guān)系
Table1Correspondencerelationshipbetweentheworkingfaceprogressandforefrontdynamicfissurem

觀測時間工作面推進(jìn)位置最前端裂縫位置日進(jìn)尺量裂縫超前距2011-05-19431 000440 20011 0009 2002011-05-20440 500449 5208 5009 0202011-05-22462 000472 13812 00010 1382011-05-23473 000481 37011 0008 3702011-05-24484 300493 17011 3008 8702011-05-25497 900508 44013 60011 5402011-05-26516 000529 24318 10014 243

2.1.2動態(tài)裂縫的發(fā)育周期

傳統(tǒng)觀點(diǎn)一般認(rèn)為，動態(tài)地裂縫隨著工作面推進(jìn)先張開而后逐漸閉合，裂縫寬度一般呈現(xiàn)由小變大、最終閉合的單峰周期，但風(fēng)積沙區(qū)補(bǔ)連塔12406工作面在高強(qiáng)度快速開采條件下，動態(tài)裂縫呈現(xiàn)出雙峰周期，且發(fā)育周期很短，實(shí)測數(shù)據(jù)如圖2所示。

圖2 動態(tài)裂縫發(fā)育周期分布Fig.2 Development cycle of dynamic ground fissures

圖2表明，初始產(chǎn)生時間不同的5條動態(tài)裂縫，各個裂縫寬度值有所不同，存在一定的差異性，裂縫寬度的最大值變化范圍為12.3～31.3 mm，平均值為19.66 mm，但其發(fā)育過程具有明顯的相似性，在采動的過程中，裂縫寬度值均呈現(xiàn)出由小變大、然后迅速變小達(dá)到初次閉合、再次開裂變大、最終完全閉合的趨勢。5條裂縫寬度值均呈現(xiàn)“M”型，包含兩個波峰，形成兩個“開裂—閉合”的過程，且第1個“開裂—閉合”過程中裂縫寬度的峰值均明顯高于第2個過程，峰值比為1.4～3.8，平均值約為2.5.

從時間尺度來看，5條裂縫的從產(chǎn)生到完全閉合的全發(fā)育周期分別為17，19，18，19，18d，平均18.2d，方差為0.84，變異系數(shù)4.61%，離散程度很小，動態(tài)裂縫發(fā)育周期基本趨于一致。從裂縫發(fā)育的兩個“開裂—閉合”的過程來看，兩個過程所經(jīng)歷的時長分別為8，10，9，9，9和9，9，9，10，9d，平均值分別為9d和9.2d，方差分別為0.71和0.45，變異系數(shù)分別為7.89%和4.89%，在檢驗(yàn)水平0.05下，兩個過程的時長差異不顯著(p<0.05)。由此可判定，對于12406工作面，動態(tài)裂縫的發(fā)育周期約為18d，且包含兩個時長近似相等的“開裂—閉合”過程。

從整體上看，12406工作面動態(tài)裂縫的發(fā)育時間和擴(kuò)展程度相對比較小，這應(yīng)該是風(fēng)積沙區(qū)高強(qiáng)度開采導(dǎo)致的，上覆較硬的巖層以及工作面的快速推進(jìn)對斷裂帶的發(fā)育起到了一定抑制作用[21]，且在快速推進(jìn)過程中，斷裂的覆巖極易形成的暫時穩(wěn)定的“力學(xué)平衡結(jié)構(gòu)”，這在一定程度也阻止了地面下沉臺階的形成，相對厚的松散層存在，也弱化了覆巖結(jié)構(gòu)變化對地表影響，縮短了裂縫擴(kuò)展的空間，進(jìn)一步促使裂縫的快速壓實(shí)閉合[22]。

2.1.3裂縫發(fā)育周期與地質(zhì)采礦條件的關(guān)系

上述的動態(tài)裂縫的發(fā)育周期應(yīng)與巖性、采深以及工作面推進(jìn)速度等采礦地質(zhì)條件有關(guān)，研究組在對動態(tài)裂縫發(fā)育周期監(jiān)測的同時，也對裂縫附近的地表移動觀測站進(jìn)行了多期的加密水準(zhǔn)觀測，以此將二者的變化過程相結(jié)合，建立動態(tài)裂縫的較為通用的定量函數(shù)模型，進(jìn)而可推廣至其他相似礦區(qū)。圖3為距離開切眼440m的 B42地表移動點(diǎn)及其附近動態(tài)裂縫1的發(fā)育特征與地表下沉規(guī)律的關(guān)系。

圖3 裂縫發(fā)育過程與地表沉陷規(guī)律的關(guān)系Fig.3 Relationship between development process of fissure and ground subsidence

當(dāng)工作面推進(jìn)至431m時，裂縫首次出現(xiàn)在B42附近，裂縫超前工作面的距離約為9m。在裂縫發(fā)育的初始階段，裂縫寬度隨地表下沉值量以及下沉速度的增大而增大，下沉速度急劇增加至209.667 mm/d時，該處地裂縫呈現(xiàn)最大表征，而當(dāng)該點(diǎn)的地表下沉速度達(dá)到最大值268.5mm/d時(工作面推進(jìn)位置538m)，此時，該處地裂縫則首次閉合，從下沉速度曲線的分段斜率來看，地表下沉最劇烈時，地裂縫經(jīng)歷了首次“開裂—閉合”階段；隨著地表下沉值的繼續(xù)增大，巖層蠕動逐漸變小，該處下沉速度開始迅速遞減，裂縫在先前斷裂處重新斷開，可能是前方地表下沉速度明顯大于該處導(dǎo)致的。后期，當(dāng)?shù)乇硐鲁林第呌谧畲笾?477 mm，下沉速度小于9mm/d時，地表活動開始進(jìn)入衰退期時，裂縫出現(xiàn)再次閉合現(xiàn)象。裂縫從開裂到完全閉合整個周期(17 d)，工作面推進(jìn)長度為209m，約為平均采深。

從上述分析可以看出，裂縫發(fā)育過程與地表下沉速度密切相關(guān)。當(dāng)裂縫出現(xiàn)第1次“開裂—閉合”過程，裂縫處的地表下沉速度增大并趨于最大值，在此過程中，工作面的持續(xù)推進(jìn)量由裂縫的超前距(Sc)和最大下沉速度滯后距L二者之和組成，由于裂縫發(fā)育周期包含兩個時長近似對等“開裂—閉合”過程，故此，裂縫發(fā)育與地質(zhì)采礦條件的關(guān)系可近似抽象為圖4，圖中，Φ為最大下沉速度滯后角；δ為動態(tài)裂縫超前角；H0為煤層埋深(采深)，m；v為平均開采速度，m/d；A，B，C和A′，B′，C′分別代表裂縫開裂、初次閉合和完全閉合的3個階段以及對應(yīng)的工作面推進(jìn)位置。

根據(jù)動態(tài)裂縫的發(fā)育規(guī)律，其發(fā)育周期為

(1)

Sc=H0/tanδ，L=H0/tanΦ

故此，動態(tài)裂縫發(fā)育周期模型為

(2)

圖4 裂縫發(fā)育周期與地質(zhì)采礦條件的關(guān)系Fig.4 Relationship between fissure development cycle and mining geological conditions

2.2 邊緣裂縫的發(fā)生發(fā)育及分布特征

圖5為邊緣裂縫的分布，其中黑色部分為工作面推進(jìn)550m時的裂縫分布情況，調(diào)查結(jié)果顯示，多條邊緣裂縫以“帶狀”形式平行于開采邊界分布。與動態(tài)裂縫不同，邊緣裂縫發(fā)生的位置滯后于工作面推進(jìn)位置，滯后距離約50m，在開采的過程中，距離工作面走向中心線90～100m的范圍最先出現(xiàn)新的邊緣裂縫，并隨著工作面的不斷推進(jìn)，裂縫向前、向外擴(kuò)展，最終在開采邊界形成裂縫帶。受鄰近采空區(qū)的影響，下山方向上裂縫的分布密度明顯小于上方方向，且同一觀測時期，下山方向的邊緣裂縫的發(fā)生位置也滯后于上山方向，滯后的距離為20～31m。工作面終采后，邊緣裂縫以“O”形圈分布于地表，進(jìn)一步體現(xiàn)了覆巖采場裂縫的分布形態(tài)[23]，開切眼與終采線位置處的裂縫近似對稱分布，下山方向的裂縫向工作面內(nèi)部收縮。

圖5 邊緣裂縫分布Fig.5 Distribution of marginal fissures

從現(xiàn)場觀測結(jié)果看，開切眼附近的裂縫開裂程度最大，形成兩條大的主裂縫，裂縫間距2.154m，裂縫中間地面下凹呈“臺階式”形狀，在接近工作面上下山邊界時，兩條主裂縫間距逐漸變小，并與上下山方向上的裂縫貫通，表2為工作面終采后穩(wěn)定邊緣裂縫相關(guān)屬性的統(tǒng)計。

表2穩(wěn)定邊緣裂縫位置及相關(guān)屬性統(tǒng)計
Table2Thepositionsandotherpropertiesofstablemarginalfissures

裂縫區(qū)域裂縫類別距邊界距離/m裂縫最大寬度/cm裂縫落差/cm裂縫帶寬/m裂縫角/(°)最外側(cè)裂縫5 580 450 23開切眼附近主裂縫-8 890～-11 04421 346 5245 8788 40最內(nèi)側(cè)裂縫-39 990 770 12最外側(cè)裂縫8 370 330 17上山方向主裂縫-12 2911 251 4649 5787 60最內(nèi)側(cè)裂縫-36 350 280 15最外側(cè)裂縫-11 740 290 10下山方向主裂縫-29 835 300 7628 5193 36最內(nèi)側(cè)裂縫-50 250 330

注：“-”表示裂縫位于工作面邊界內(nèi)側(cè)。

在沒有鄰近工作面開采的影響下，主裂縫分布在工作面內(nèi)部距離邊界10m左右的范圍內(nèi)，裂縫的寬度以及落差相對較大，除開切眼附近，地表均無明顯的下沉臺階，裂縫的帶寬在46～50m。下山方向受老采空區(qū)的影響，裂縫均位于開采邊界內(nèi)部，其分布范圍約為上山方向的60%，且裂縫的最大寬度、落差等相關(guān)屬性均相對較弱，顯示特征不明顯。傳統(tǒng)意義下，邊緣裂縫主要分布在塌陷盆地的外邊緣區(qū)，而風(fēng)積沙區(qū)高強(qiáng)度開采(寬深比大、回采率高、推進(jìn)速度快)引起的邊緣裂縫則主要分布在工作面邊界的內(nèi)部，整體向內(nèi)收縮，裂縫角呈高角度甚至近似于垂直角，地表移動范圍減小且變形集中，譚志祥等[24]針對中東部高強(qiáng)度綜放開采地表沉陷規(guī)律研究也得到相似的結(jié)論。

2.3 裂縫寬度與深度的關(guān)系

工作面內(nèi)部的動態(tài)裂縫隨著開采工作的進(jìn)行，最終完全呈現(xiàn)閉合狀態(tài)，而邊緣裂縫在無外力擾動下，會持續(xù)存在于地表，故選擇穩(wěn)定的邊緣裂縫帶的若干作為研究對象，進(jìn)行石灰漿灌注(除無法開挖的主裂縫)，共開挖40個注漿點(diǎn)，其中開切點(diǎn)附近15個，上山方向25個，裂縫寬度為13.2～58.8mm，深度為68.4～1120.8mm，如圖6所示。

圖6 邊緣裂縫寬度與深度對應(yīng)信息Fig.6 Correspondent information of width and depth of marginal fissures

以往的研究表明，裂縫的寬度(d)和深度(h)二者經(jīng)驗(yàn)[25]公式為

其中,A為系數(shù)，黃土耕地取8，風(fēng)化基巖取15。圖6說明，補(bǔ)連塔12406工作面裂縫寬度與深度的關(guān)系比較復(fù)雜，并不符合上述經(jīng)驗(yàn)公式，其原因在于裂縫深度的監(jiān)測比較困難，監(jiān)測數(shù)據(jù)誤差大，因此，需要進(jìn)一步深入探討裂縫深度的監(jiān)測方法，從而進(jìn)一步完善裂縫與寬度的關(guān)系。

3 結(jié) 論

(1)動態(tài)地裂縫往往超前于當(dāng)前工作面的開采位置向前發(fā)展，最前端動態(tài)裂縫的超前距離(Sc)與日推進(jìn)速度(v)呈現(xiàn)明顯的線性正相關(guān)關(guān)系：Sc=0.637v+2.715，平均的裂縫超前距為10.359m，超前裂縫角δ為近似垂直角。

(2)發(fā)現(xiàn)動態(tài)地裂縫的發(fā)育為包含兩個時長近似相等的“開裂—閉合”過程，即裂縫的發(fā)育呈 “M”型的雙峰波形，且第1個峰值明顯大于第2個峰值，峰值比為1.4～3.8；裂縫發(fā)育周期約為18d，此時工作面推進(jìn)距離約等同于平均采深，說明高強(qiáng)度開采動態(tài)地裂縫發(fā)育較快，并呈現(xiàn)明顯的自修復(fù)現(xiàn)象。

(3)動態(tài)裂縫發(fā)育過程與地質(zhì)采礦條件密切相關(guān)，地表下沉速度趨于最大值時，該處的動態(tài)地裂縫則首次閉合，裂縫發(fā)育周期T與裂縫超前距(Sc)以及最大下沉速度滯后距L成正比，與日平均開采速度v成反比，并基于此，建立了T與采礦地質(zhì)條件的函數(shù)模型：T=2H0(1/tanδ+ 1/tanΦ)/v。

(4)邊緣裂縫以“帶狀”形式、“O”形圈的形態(tài)分布在工作面的開采邊界，主要分布在工作面內(nèi)部距邊界40m的范圍內(nèi)，裂縫帶寬為46～50m，受臨近工作面采動影響的邊緣裂縫發(fā)育帶寬度明顯減少，帶寬減為28.51m，減少了約40%。說明風(fēng)積沙區(qū)高強(qiáng)度開采導(dǎo)致邊緣裂縫帶整體向工作面內(nèi)部收縮，臨近工作面采動會減輕原有地裂縫的影響，具有一定的自修復(fù)功能。

(5)裂縫寬度和深度不符合以往的經(jīng)驗(yàn)公式，二者之間的關(guān)系相對復(fù)雜，裂縫深度的探測方法還有待進(jìn)一步深入研究。

風(fēng)積沙區(qū)高強(qiáng)度的開采使地裂縫發(fā)生發(fā)育規(guī)律呈現(xiàn)出了全新的特征，動態(tài)地裂縫具有快速閉合的自修復(fù)特征，不需要人工修復(fù)，而邊緣裂縫呈帶狀向工作面內(nèi)部收縮，是地表環(huán)境損傷和人工修復(fù)的重點(diǎn)區(qū)域，本文的研究僅僅是拋磚引玉，期待更多的實(shí)例研究成果和更為深入探討在今后陸續(xù)發(fā)表。

[1] 范立民.煤礦地裂縫研究(三輯)[M].北京:地震出版社,1995.

[2] 趙紅梅.采礦塌陷條件下包氣帶土壤水分布與動態(tài)變化特征研究[D].北京:中國地質(zhì)科學(xué)院,2006.

[3] 臧蔭桐,汪 季,丁國棟,等.采煤沉陷后風(fēng)沙土理化性質(zhì)變化及其評價研究[J].土壤學(xué)報,2010,47(2):262-269. Zang Yintong,Wang Ji,Ding Guodong,et al.Variation of physic-chemical properties of aeolian sandy soil at coal mining subsidence and its evaluation[J].Acta Pedologica Sinica,2010,47(2):262-269.

[4] 王文龍,李占斌,張平倉.神府東勝煤田開發(fā)中誘發(fā)的環(huán)境災(zāi)害問題研究[J].生態(tài)學(xué)雜志,2004,23(1):34-38. Wang Wenlong,Li Zhanbin,Zhang Pingcang.Environmental disaster issues induced by coal exploitation in Shenfu Dongsheng coalfield[J].Chinese Journal of Ecology,2004,23(1):34-38.

[5] 張發(fā)旺,侯新偉,韓占濤,等.煤炭開發(fā)引起水土環(huán)境演化及其調(diào)控技術(shù)[J].地球?qū)W報,2011,22(4):345-350. Zhang Fawang,Hou Xinwei,Han Zhantao,et al.Water-soil environmental evolution induced by coal mining and its regulating and controlling techniques[J].Acta Geoscientica Sinica,2011,22(4):345-350.

[6] 李文平,段中會,華解明,等.陜北榆神府礦區(qū)地質(zhì)環(huán)境現(xiàn)狀及采煤效應(yīng)影響預(yù)測[J].工程地質(zhì)學(xué)報,2000,8(3):324-333. Li Wenping,Duan Zhonghui,Hua Jieming,et al.Evaluation of present geological environment and prediction of its variation caused by mining in Yushenfu mine area of North Shanxi[J].Journal of Engineering Geology,2000,8(3):324-333.

[7] 呂晶潔,胡春元,賀 曉.采煤塌陷對固定沙丘土壤水分動態(tài)的影響研究[J].干旱區(qū)資源與環(huán)境,2005,19(7):152-156. Lü Jingjie,Hu Chunyuan,He Xiao.Study on the influences of coal mining collapse on the soil water dynamics of fixed dune[J].Journal of Arid Land Resources and Environment,2005,19(7):152-156.

[8] 周 瑩,賀 曉,徐 軍,等.半干旱區(qū)采煤沉陷對地表植被組成及多樣性的影響[J].生態(tài)學(xué)報,2009,29(8):4517-4525. Zhou Ying,He Xiao,Xu Jun,et al.Effects of coal mining subsidence on vegetation composition and plant diversity in semi-arid region[J].Acta Ecological Sinica,2009,29(8):4517-4525.

[9] 雷少剛.荒漠礦區(qū)關(guān)鍵環(huán)境要素的監(jiān)測與采動影響規(guī)律研究[J].煤炭學(xué)報,2010,35(9):1587-1588. Lei Shaogang.Monitoring and analyzing the mining impacts on key environmental elements in desert area[J].Journal of China Coal Society,2010,35(9):1587-1588.

[10] 欒長青,唐益群,趙法鎖,等.陜西韓城礦區(qū)地表沉降變形研究[J].自然災(zāi)害學(xué)報,2007,16(3):81-85. Luan Changqing,Tang Yiqun,Zhao Fasuo,et al.Study on surface subsidence and deformation in mineral region in Hancheng City,Shaanxi Province[J].Journal of Natural Disasters,2007,16(3):81-85.

[11] 李彥軍.靖遠(yuǎn)礦區(qū)地面塌陷地裂縫特征與治理研究[J].山西建筑,2008,34(29):144-145. Li Yanjun.The characteristic of ground collapse and ground fissure and treatment research in Jingyuan mining area[J].Shanxi Architecture,2008,34(29):144-145.

[12] 杜善周.神東礦區(qū)大規(guī)模開采的覆巖破壞與地表移動規(guī)律研究[D].北京:中國礦業(yè)大學(xué)(北京),2010.

[13] 李 曉,路世豹,廖秋林,等.充填法開采引起的地裂縫分布特征與現(xiàn)場監(jiān)測分析[J].巖石力學(xué)與工程學(xué)報，2006,25(7):1361-1369. Li Xiao,Lu Shibao,Liao Qiulin,et al.Distribution characteristics and field monitoring of ground fissures caused by backfill mining[J].Chinese Journal of Rock Mechanics and Engineering,2006,25(7):1361-1369.

[14] 吳 侃,胡振琪,常 江,等.開采引起的地表裂縫分布規(guī)律[J].中國礦業(yè)大學(xué)學(xué)報,1997,26(2):56-59. Wu Kan,Hu Zhenqi,Chang Jiang,et al.Distribution law of ground crack induced by coal mining[J].Journal of China University of Mining & Technology,1997,26(2):56-59.

[15] 吳 侃,周 鳴,胡振琪.開采引起的地表裂縫深度和寬度預(yù)計[J].遼寧工程技術(shù)大學(xué)學(xué)報(自然科學(xué)版),1997,16(6):649-652. Wu Kan,Zhou Ming,Hu Zhenqi.The prediction of ground fissure depth and width by mining[J].Journal of Liaoning Technical University(Natural Science Edition),1997,16(6):649-652.

[16] 李 亮,吳 侃,陳冉麗,等.小波分析在開采沉陷區(qū)地表裂縫信息提取的應(yīng)用[J].測繪科學(xué),2010,35(1):165-168. Li Liang,Wu Kan,Chen Xiaoli,et al.Information extraction of surface crack position in mining subsidence area based on wavelet transform[J].Science of Surveying and Mapping,2010,35(1):165-168.

[17] 趙超英,張 勤,丁曉利,等.利用InSAR技術(shù)定位西安活動地裂縫[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2009,34(7):809-814. Zhao Chaoying,Zhang Qin,Ding Xiaoli,et al.Positioning of Xi’an active ground fissures with SAR interferometry[J].Geomatics and Information Science of Wuhan University,2009,34(7):809-814.

[18] Herrera G,lvarez Fernández M I,Tomás R,et al.Forensic analysis of buildings affected by mining subsidence based on differential interferometry[J].Engineering Failure Analysis,2012,24(1):67-76.

[19] 楊 科,謝廣祥.綜放開采采動裂隙分布及其演化特征分析[J].礦業(yè)安全與環(huán)保,2009,36(4):1-4. Yang Ke,Xie Guangxiang.Distribution of mining fissures induced by fully mechanized caving mining and analysis on their evolvement characteristics[J].Ming Safety & Environmental Protection,2009,36(4):1-4.

[20] 張玉軍,李鳳明.高強(qiáng)度綜放開采采動覆巖破壞高度及裂隙發(fā)育演化監(jiān)測分析[J].巖石力學(xué)與工程學(xué)報,2011,30(S1):2294-2301. Zhang Yujun,Li Fengming.Monitoring analysis of fissure development evolution and height of overburden failure of high tension fully-mechanized caving mining[J].Chinese Journal of Rock Mechanics and Engineering,2011,30(S1):2294-2301.

[21] 彭永偉,齊慶新,李宏艷,等.高強(qiáng)度地下開采對巖體斷裂帶高度影響因素的數(shù)值模擬分析[J].煤炭學(xué)報,2009,34(2):145-149. Peng Yongwei,Qi Qingxin,Li Hongyan,et al.Numerical simulation and analysis on influencing factors of fractured zone height of high-strength underground mining[J].Journal of China Coal Society,2009,34(2):145-149.

[22] 范鋼偉,張東升,馬立強(qiáng).神東礦區(qū)淺埋煤層開采覆巖移動與裂隙分布特征[J].中國礦業(yè)大學(xué)學(xué)報,2011,40(2):196-202. Fan Gangwei,Zhang Dongsheng,Ma Liqiang.Overburden movement and fracture distribution induced by longwall mining of shallow coal seam in the Shendong coalfield[J].Journal of China University of Mining & Technology,2011,40(2):196-202.

[23] 錢鳴高,許家林.覆巖采動裂隙分布的“O”形圈特征研究[J].煤炭學(xué)報,1998,23(5):466-469. Qian Minggao,Xu Jialin.Study on the “O Shape” circle distribution characteristics of mining induced fractures in the overlaying strata[ J].Journal of China Coal Society,1998,23(5):466-469.

[24] 譚志祥,王宗勝,李運(yùn)江,等.高強(qiáng)度綜放開采地表沉陷規(guī)律實(shí)測研究[J].采礦與安全工程學(xué)報,2008,25(1):59-63. Tan Zhixiang,Wang Zongsheng,Li Yunjiang,et al.Field research on ground subsidence rules of intensive fully mechanized mining by sublevel caving[J].Journal of Mining & Safety Engineering,2008,25(1):59-63.

[25] 趙勤正,張和生.基于特征的井工開采土地破壞程度特征因子選取[J].能源環(huán)境保護(hù),2003,17(6):54-57. Zhao Qinzheng,Zhang Hesheng.Feature factor of land destroy degree are chosen in mining on base of features in semi-dry mountain areas[J].Energy Environmental Protection,2003,17(6):54-57.

Distributioncharacteristicanddevelopmentrulesofgroundfissuresduetocoalmininginwindyandsandyregion

HU Zhen-qi1,WANG Xin-jing1,HE An-min2

(1.InstituteofLandReclamationandEcologicalRestoration,ChinaUniversityofMiningandTechnology(Beijing),Beijing100083,China；2.ShenhuaShendongCoalGroupCorporationLimited,Shenmu719315,China)

Coal mining might cause serious environmental damage in ecologically fragile region.In sandy and windy region,ground fissure is the most typical characteristic of ecological damages caused by underground coal mining.However,few works and reports are related to the generation and whole development process of ground fissures caused by high tension and rapid caving mining.Through the united space coordinates control system at surface-underground,a comprehensive observation method was invented to monitor ground fissure on the No.12406 working face of Bulianta coalfield,distribution characteristic of marginal fissures,genesis and development rules of dynamic fissures and their relationship with mining geological conditions were revealed.The results show that:① Dynamic fissures always develop ahead of the current working location,the crossover distance has a significant positive linear correlation with size of daily footage.The average fissure crossover distance is 10.359m and advance angle of break is approach to vertical angle.② Under the high tension and rapid caving mining,dynamic ground fissures havn’t obvious heads in vertical direction,those widths had similar cyclical changes during the mining process and present ‘M’ type with two different peaks.The first peak(maximum width)is 1.4-3.8times of the second one.The life-cycle of a dynamic fissure is 18days，including two cracking-closure processes with similar duration.On this basis,a common function model for the life-cycle was established based on mining and geological conditions.It indicates that dynamic fissures have fast ‘self-healing’features and manual repairmen might not be necessary.③ Marginal fissures occurrs in strips and the ‘O’ shape inside of boundary working face,and the width of strip of marginal fissures is about 46-50m.If the adjacent coal face excavates,the width of the strip of marginal fissures in this side will be narrowed 40% due to the double disturbance.The distribution region of marginal fissures will exist long time and need manual restoration.

mining subsidence;ground fissures;ecological damages;windy and sandy regions

10.13225/j.cnki.jccs.2013.1289

國家自然科學(xué)基金委員會-神華集團(tuán)有限公司煤炭聯(lián)合基金資助項(xiàng)目(U1361203);中國礦業(yè)大學(xué)(北京)博士生拔尖創(chuàng)新人才培育基金資助項(xiàng)目(800015Z628)

胡振琪(1963—)，男，安徽五河人，教授，博士生導(dǎo)師。E-mail:huzq@cumtb.edu.cn

TD167

A

0253-9993(2014)01-0011-08

胡振琪，王新靜，賀安民.風(fēng)積沙區(qū)采煤沉陷地裂縫分布特征與發(fā)生發(fā)育規(guī)律[J].煤炭學(xué)報,2014,39(1):11-18.

Hu Zhenqi,Wang Xinjing,He Anming.Distribution characteristic and development rules of ground fissures due to coal mining in windy and sandy region[J].Journal of China Coal Society,2014,39(1):11-18.doi:10.13225/j.cnki.jccs.2013.1289