張海君 吳奕樞 王飛 謝昆侖

摘 要：在特厚煤層綜放開采過程中，工作面各上覆巖層受其影響相繼發(fā)生變形、垮落，使地表出現(xiàn)不均勻沉降，產(chǎn)生塌陷、裂縫等，嚴(yán)重威脅礦井及其周邊地區(qū)人員和生產(chǎn)設(shè)施安全。因此，地表沉陷變化規(guī)律的研究對沉陷區(qū)范圍預(yù)測和治理均有重要意義。文中以準(zhǔn)格爾煤田龍王溝煤礦61601特厚煤層綜放工作面為研究對象，建立基于概率積分法的地表沉陷預(yù)計(jì)模型，并對計(jì)算結(jié)果加以分析，結(jié)果顯示特厚煤層綜放開采地表下沉情況基本符合概率積分法預(yù)計(jì)模型，但其發(fā)生沉陷的區(qū)域范圍更小，沉陷區(qū)域邊界坡度更大，特厚煤層綜放開采地表沉陷情況相較一般開采呈現(xiàn)出下沉速度快、下沉幅度劇烈、變形分別集中的特征。同時運(yùn)用FLAC^3D數(shù)值模擬技術(shù)對地表沉陷形成機(jī)制和動態(tài)變化過程進(jìn)行補(bǔ)充，F(xiàn)LAC^3D數(shù)值模擬能更為全面地體現(xiàn)地表沉陷規(guī)律。此外，文中還進(jìn)一步結(jié)合FLAC^3D數(shù)值模擬和實(shí)測沉陷結(jié)果分析了不同推進(jìn)長度下走向和傾向主斷面地表下沉曲線、走向地面測點(diǎn)下沉速度變化情況、開采區(qū)域總體位移情況等，歸納出特厚煤層綜放開采地表沉陷規(guī)律及特征，為合理規(guī)劃煤礦開采作業(yè)方案提供理論依據(jù)。關(guān)鍵詞：特厚煤層；地表沉陷；FLAC^3D；概率積分法中圖分類號：TD 822

文獻(xiàn)標(biāo)志碼：A

文章編號：1672-9315（2022）05-0874-10

DOI：10.13800/j.cnki.xakjdxxb.2022.0505開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Surface subsidence law in fully-mechanized caving mining of extra-thick seam in Longwanggou Coal Mine of Zhungeer coalfield

ZHANG Haijun¹，WU Yishu²，WANG Fei²，XIE Kunlun²

（1.Ordos Guoyuan Mining Development Co.，Ltd.，Ordos 017000，China；2.School of Mechanics and Safety Engineering，Zhengzhou University，Zhengzhou 450001，China）

Abstract：During the fully-mechanized caving mining of extra-thick coal seams，the overlying strata of the working face are successively deformed and collapsed under the influence of mining，resulting in uneven settlement，collapse，cracks and other accidents on the surface，which seriously threatens the safety of personnel and production facilities in the mine and its surrounding areas.Therefore，it is of great significance to study the variation law of surface subsidence for the prediction and control of the subsidence area.In this paper，the surface subsidence prediction model based on the probability integral method was established for the 61601 extra-thick coal seam comprehensive mining face of Longwanggou coal mine in Zhungeer coalfield.The results showed that the surface subsidence of the fully-mechanized caving mining of extra-thick coal seam was basically in accordance with the prediction model of the probability integral method，but the area where subsidence occurs is smaller and the slope of the boundary of the subsidence area is larger.The surface subsidence of fully-mechanized caving mining of extra-thick coal seams presents the characteristics of fast subsidence speed，intense subsidence magnitude and concentrated deformation respectively compared with the general mining. And the FLAC^3Dnumerical simulation technology was applied to supplement the formation mechanism as well as dynamic change process of surface subsidence，which show that FLAC^3Dnumerical simulation can reflect the surface subsidence law more comprehensively.In addition，a detailed analysis has been made has been made of the followings：the surface subsidence curves of the strike and inclination main sections under different propulsion lengths，the changes of the subsidence speed of the strike ground measurement points，and the overall displacement of the mining area by combining the results of FLAC^3Dnumerical simulation and the measured subsidence，with the surface subsidence law summarized and the characteristics revealed of comprehensive mining in extra-thick coal seams，which provide a theoretical basis for the reasonable planning of coal mining operations.

Key words：extra-thick coal seam；surface subsidence；FLAC^3D；probability integral method

0 引 言井下煤炭資源采出后，煤層頂板懸露，頂板原巖應(yīng)力遭到擾動，頂板巖石在重力作用下將出現(xiàn)移動變形。隨著工作面的繼續(xù)推進(jìn)，頂板出現(xiàn)破壞、冒落，導(dǎo)致地表沉陷現(xiàn)象的發(fā)生，將會對礦區(qū)及礦區(qū)周邊的建筑物、生產(chǎn)設(shè)施和人員的安全造成威脅。隨著中國煤炭資源的持續(xù)勘探和開發(fā)，煤炭開采所面對的礦區(qū)環(huán)境更加復(fù)雜，如何有效預(yù)測和防范開采造成的地表沉降危害，已成為礦區(qū)整體規(guī)劃設(shè)計(jì)的一個重點(diǎn)。井下開采所致的工作面頂板巖移及地表沉陷現(xiàn)象是一個涉及復(fù)雜巖體力學(xué)的變化過程。以往國內(nèi)外學(xué)者提出了許多地表沉陷模型和理論^[1-4]，典型的如基于觀測數(shù)據(jù)的正態(tài)理論、二等分線理論、自然斜面理論和開采沉陷“拱理論”^[5-6]、工作面上覆巖層“三帶”分布理論^[7]等。隨著科學(xué)技術(shù)水平的不斷提升，開采沉陷理論得到進(jìn)一步發(fā)展。Salamon將影響函數(shù)預(yù)計(jì)法與連續(xù)介質(zhì)理論相結(jié)合，開創(chuàng)了用于探究開采沉陷規(guī)律的邊界元法（BEM）^[8]；DONNELLY等基于SWIFT技術(shù)對礦山地表沉陷問題進(jìn)行預(yù)測與分析^[9]，BARYAKH等基于線性粘彈性理論，利用蠕變函數(shù)開發(fā)了一種動態(tài)地表沉陷預(yù)計(jì)程序^[10]；崔希民等應(yīng)用流變模型對地表沉降規(guī)律進(jìn)行分析^[11]；朱廣軼等基于概率積分法提出了一種動態(tài)地表沉陷時間函數(shù)，并通過實(shí)例分析證明其具有比Knote函數(shù)更優(yōu)的精準(zhǔn)度^[12]；余學(xué)義等根據(jù)觀測數(shù)據(jù)研究單一工作面開采和多工作面開采后地表的沉陷規(guī)律，并分別評估了2種開采場景下地表的損害情況^[13]；劉玉成對Knote時間函數(shù)加以改良，建立了一種更加精準(zhǔn)的地表沉陷預(yù)計(jì)模型^[14]；石曉宇等探討了2種灰色模型在地表沉陷監(jiān)測與預(yù)計(jì)方面的應(yīng)用，并指出二者適用的開采情況^[15]；黃明江等指出地表下沉曲線符合分形增長規(guī)律，并借助分形插值法實(shí)現(xiàn)了有限觀測數(shù)據(jù)下對地表沉陷連續(xù)變化情況的擬合^[16]。此外，一些學(xué)者還采用GIS，InSAR，ArcObjects，RTK，無人機(jī)激光雷達(dá)等新型觀測技術(shù)和地理信息系統(tǒng)對礦山地表沉陷情況展開監(jiān)測與分析^[17-20]。中國許多大型礦區(qū)屬于厚煤層，煤炭儲量相當(dāng)豐富。近年來，綜放開采技術(shù)被運(yùn)用于厚煤層的開采過程中，因其具有產(chǎn)量大、效率高、成本低等特點(diǎn)，而備受業(yè)內(nèi)認(rèn)可。關(guān)于厚煤層開采引起的覆巖和地表移動規(guī)律等方面目前已有不少研究成果，如郭文兵等對“三軟”特厚煤層地表沉陷情況進(jìn)行觀測，并基于概率積分法運(yùn)用Matlab擬合了地表下沉曲線^[21]；胡青峰等建立了塔山煤礦特厚煤層傾向主斷面相似模擬實(shí)驗(yàn)，并研究了重復(fù)采動下覆巖變形和地表沉陷規(guī)律^[22]；高超等針對東坡煤礦淺埋深、特厚煤層綜采等特點(diǎn)，結(jié)合地表實(shí)測數(shù)據(jù)，對地表動態(tài)移動變形特征進(jìn)行分析^[23]；趙兵朝等通過對郭家河煤礦特厚煤層綜放開采實(shí)地觀測和分析求取了地表移動變形有關(guān)參數(shù)，并綜合運(yùn)用反演模擬和邊坡穩(wěn)定性分析對所得參數(shù)進(jìn)行修正^[24]，但所研究的對象較少涉及20 m以上的大采放比特厚煤層綜放開采，故仍需對此類作業(yè)引起的覆巖與地表移動規(guī)律進(jìn)行系統(tǒng)地研究。

文中將根據(jù)準(zhǔn)格爾煤田龍王溝煤礦特厚煤層一次采全厚這一開采特征，建立基于概率積分法的地表沉陷預(yù)計(jì)模型，運(yùn)用該模型分析煤層開采引起的地表移動規(guī)律與變形特征，同時討論該模型存在的局限性，并采用FLAC^3D數(shù)值模擬方法加以補(bǔ)充，從而獲得較為準(zhǔn)確、全面的特厚煤層綜放開采的地表沉陷規(guī)律及其特征，為合理規(guī)劃煤礦開采作業(yè)方案，保護(hù)礦區(qū)人員生產(chǎn)生活的安全提供指導(dǎo)方向和理論支持。

1 龍王溝煤礦61601工作面概述龍王溝井田位于內(nèi)蒙古自治區(qū)鄂爾多斯市準(zhǔn)格爾煤田中北部，面積41.35 km²，屬典型黃土高原地貌。依據(jù)地表出露及鉆孔揭露情況，該井田地層自下而上分別為：奧陶系中統(tǒng)馬家溝組（O₂m），石炭系上統(tǒng)太原組（C₂t），二疊系下統(tǒng)山西組（P₁s）、下石盒子組（P₁x），二疊系上統(tǒng)上石盒子組（P₂s）、石千峰組（P₂sh），三疊系下統(tǒng)劉家溝組（T₁l），第三系上新統(tǒng)（N₂）和第四系（Q）等。煤田地層結(jié)構(gòu)形態(tài)為一近似于南北走向、呈波狀起伏、向西側(cè)傾斜的單斜構(gòu)造，地層傾角較小，低于10°。61601工作面布置于龍王溝煤田太原組上部6煤處，為首采面，平均采深399.9 m，工作面傾向長度254.6 m，走向可采長度624 m，煤均厚23.15 m，煤層穩(wěn)定。工作面設(shè)計(jì)采高5.1 m，頂煤厚18.05 m，采放比1∶3.54。煤層傾角0～2°，接近水平煤層。煤層頂?shù)装甯鲙r層如圖1所示。61601工作面開采選用后退式走向長壁采煤方法，綜采放頂煤一次采全高開采工藝，完全垮落法管理頂板。為獲得特厚煤層開采條件下地表沉陷的實(shí)測數(shù)據(jù)，分析地表沉陷規(guī)律，根據(jù)《煤礦測量規(guī)程》在61601工作面外設(shè)置地表沉陷觀測站，根據(jù)《煤礦測量規(guī)程》規(guī)定，當(dāng)工作面走向長度大于0.9H₀（H₀為煤層埋深，399.9 m）時，可只設(shè)半條走向觀測線，故本工作面設(shè)置半條走向觀測線。觀測點(diǎn)按便于觀測塌陷區(qū)的下沉和能夠準(zhǔn)確反映工作面下沉情況進(jìn)行布置，并以觀測站為起點(diǎn)布置觀測點(diǎn)，首個觀測點(diǎn)距開切眼上方地表水平距離288 m。之后每隔25 m設(shè)置一個觀測點(diǎn)，直至工作面走向推進(jìn)長度中點(diǎn)之后25 m，以確保觀測數(shù)

據(jù)能夠監(jiān)測到下沉量的極值點(diǎn)。觀測線全長625 m，

各點(diǎn)標(biāo)記為T1，T2，T3，…，T26。各測點(diǎn)的布置如圖2所示，圖中矩形部分為61601工作面示意圖，三角形中心為觀測點(diǎn)位置。走向測線下沉測量工作起始于2019年1月17日，測量周期1年5個月，各測點(diǎn)平均每月測量一次，直至測量完畢。

2 基于概率積分法的地表沉陷預(yù)計(jì)模型概率積分法是基于隨機(jī)介質(zhì)理論發(fā)展而來的巖體移動與變形計(jì)算模型，經(jīng)典的概率積分法表述的地表任意點(diǎn)下沉位移計(jì)算公式見式（1）。

公式（1）表示地表任意點(diǎn)（x，y）在開采區(qū)域?yàn)椋▁₁-x₀）*（y₁-y₀）的煤層開采后產(chǎn)生的下沉量。W_max為地表最大下沉量；x₀，x₁，y₀，y₁分別是煤層沿傾向開采和走向開采的起止點(diǎn)；r為主要影響半徑，可由公式（2）得到。

r=Htanβ（2）

式中 H為采深；tanβ為主要影響角正切值。文中使用Python Integrate模塊下的Dblquad函數(shù)編寫基于概率積分法的地表任意點(diǎn)移動與變形模擬程序，設(shè)地表為一傾向長945 m，走向長1 392 m的水平平面，并取左下角為坐標(biāo)原點(diǎn)（0，0）。平面采用矩形網(wǎng)格進(jìn)行劃分，依照61601工作面實(shí)測數(shù)據(jù)確定開采范圍和巖移參數(shù)，煤層厚23.15 m，下沉系數(shù)0.53，主要影響角正切值2.9。地表下傾向開采區(qū)域位于345～600 m之間，走向開采區(qū)域位于384～1 008 m之間。

模擬結(jié)果將保存為記錄有地表各節(jié)點(diǎn)坐標(biāo)及下沉量的數(shù)據(jù)文件。將該文件導(dǎo)入Surfer中，繪制地表開采下沉量等值線圖，如圖3所示，同時比較計(jì)算值與實(shí)測值，繪制相應(yīng)的地表下沉曲線圖，如圖4所示（部分觀測點(diǎn)受采動影響發(fā)生破壞，這些點(diǎn)未在圖中繪出）。

比照實(shí)際觀測成果可知，61601工作面地表下沉變形趨勢總體上與概率積分法沉陷預(yù)計(jì)模型相吻合。從圖3可看出，采動后地表出現(xiàn)了一個以工作面所在位置上方為中心并向外擴(kuò)大的類橢圓形沉降區(qū)域。地表下沉值較大，達(dá)到-13.32 m，等值線在-2.5～-10 m之間密集分布，顯示此區(qū)間內(nèi)地表高程存在較大幅度的下移變化，形成了坡度較陡的地表沉降區(qū)域邊界，下沉量大于-10 m的等值線則較為稀疏，表明地表形成了上寬下窄的“漏斗狀”下沉盆地。從模型的巖移參數(shù)上看，相比普通煤層的開采^[25-26]，特厚煤層綜放開采的地表下沉系數(shù)較大，邊界角值比較小，主要影響角正切值顯著偏大，反映出大采高、快速推進(jìn)的綜采放頂煤開采地表下沉速度快，巖體動態(tài)變形幅度大、運(yùn)動劇烈的特征。此外，對比發(fā)現(xiàn)實(shí)測值具有更為陡峭的地表下沉邊界，更為集中的地表下沉區(qū)域，實(shí)測地表發(fā)生沉陷的范圍更小，下沉曲線在盆地邊界收斂更快，充分地體現(xiàn)出特厚煤層綜放開采地表沉陷相比一般地表沉陷規(guī)律具有下沉幅度劇烈、變形集中等特點(diǎn)。概率積分法具有公式簡單、計(jì)算高效快捷的特點(diǎn)，其只需借助較少的巖移參數(shù)，就能計(jì)算地表沉陷情況，為研究地質(zhì)條件簡單、巖層達(dá)到穩(wěn)定狀態(tài)的工作面沉陷狀況提供了快速預(yù)計(jì)和分析的方法。但是概率積分法是一種基于隨機(jī)介質(zhì)理論的理想化模型，其不能充分考慮開采區(qū)各巖層的空間構(gòu)型、力學(xué)特征的差異對地表移動

變形產(chǎn)生的影響；同時作為一種靜態(tài)地表沉陷預(yù)計(jì)方法，概率積分法難以直接建立地表移動變化趨勢與進(jìn)尺長度、時間等因素之間的關(guān)系等。為提高地表沉陷預(yù)計(jì)模型的精準(zhǔn)性，進(jìn)一步解釋實(shí)地觀測和概率積分法模型預(yù)計(jì)結(jié)果的形成機(jī)理，從動態(tài)的角度探究開采沉陷隨時間的變化規(guī)律，這里討論采用FLAC^3D數(shù)值模擬方法對地表沉陷預(yù)計(jì)模型進(jìn)行補(bǔ)充。

3 基于FLAC^3D的61601工作面地表沉陷規(guī)律模擬

3.1 FLAC^3D數(shù)值模擬模型的建立根據(jù)61601工作面相關(guān)地質(zhì)資料建立幾何模型，如圖5所示。該模型沿X軸方向全長1 195.4 m，對應(yīng)開采工作面的傾向方向；沿Y軸方向全長2 000 m，對應(yīng)開采工作面的走向方向。模型依據(jù)各地層實(shí)地勘測到的巖土組成情況分別建立了28個地層，并根據(jù)地質(zhì)結(jié)構(gòu)添加部分節(jié)理接觸面以提高精度，總深度417.8 m。其中開采工作面位于第24巖層（六煤）的下部，開采區(qū)域所在位置為：沿X軸方向0～254.6 m，沿Y軸方向0～624 m。模型選用矩形網(wǎng)格進(jìn)行劃分，根據(jù)礦區(qū)巖層傾角很小這一特征將模型巖層形態(tài)設(shè)置為水平巖層。

以往的研究表明，直接將巖土力學(xué)試驗(yàn)獲得的力學(xué)參數(shù)應(yīng)用于數(shù)值模擬當(dāng)中往往容易發(fā)生錯誤，通常需取實(shí)際力學(xué)參數(shù)的1/3～1/5作為數(shù)值模擬的模擬值，并進(jìn)行反演分析多次計(jì)算方可取得較為準(zhǔn)確的模型力學(xué)參數(shù)^[27-28]。根據(jù)模擬比例經(jīng)驗(yàn)值運(yùn)用反演計(jì)算確定巖層模型最優(yōu)力學(xué)參數(shù)，并導(dǎo)入模型中，主要巖層模型的力學(xué)參數(shù)見表1。

本次模擬基于摩爾-庫倫彈塑性模型進(jìn)行預(yù)測計(jì)算，分別記錄工作面推進(jìn)169 m（開挖第37天）、244 m（開挖第77天）、342 m（開挖第181天）、496 m（開挖第310天）和624 m（開挖第369天）時地表的走向下沉情況。開挖模型如圖6所示，其中2個三角形的位置對應(yīng)地表下沉觀測點(diǎn)的起止范圍（T1～T26），模型下部的矩形代表61601工作面投影，模型開挖進(jìn)度與61601工作面實(shí)際開挖進(jìn)度保持一致。

3.2 結(jié)果及分析選取工作面進(jìn)尺169，244，342，496和624 m時的走向地表Z向位移擬合數(shù)據(jù)，與地表觀測站實(shí)測下沉量進(jìn)行比照（表2），繪制下沉曲線瀑布圖，如圖7所示（部分觀測點(diǎn)受采動影響發(fā)生破壞掩埋無法記錄，這些點(diǎn)未在圖中繪出）。

從圖7（e）數(shù)據(jù)可看出，開采完畢后FLAC^3D地表沉陷預(yù)計(jì)模型的最大下沉量為-13.39m，與實(shí)測值（-13.32 m）相差小于5%，可見FLAC^3D數(shù)值模擬在最大下沉量方面顯示出了很好的擬合度。實(shí)測的地表最大下沉點(diǎn)位于距觀測線起點(diǎn)525 m的T22號觀測點(diǎn)附近，而數(shù)值模擬的最大下沉點(diǎn)位于地表平面的中心點(diǎn)（T25點(diǎn)），這可能是因?yàn)樽匀粭l件下礦區(qū)巖體仍存在一定傾角以及巖層間存在斷層等不連續(xù)結(jié)構(gòu)所致。此外，與模擬值相比，開采初期（圖7（a）和（b））實(shí)測下沉曲線在拐點(diǎn)后的切線斜率迅速增大，形成更陡峭的下沉曲線，但隨著工作面的推進(jìn)，模擬下沉曲線與實(shí)測數(shù)據(jù)間的差距持續(xù)縮?。▓D7（c）和（d）），至停采線位置時（圖7（e）），二者已基本趨于一致。總體而言，F(xiàn)LAC^3D對特厚煤層綜放開采地表沉陷的數(shù)模模擬效果相比經(jīng)典概率積分法預(yù)計(jì)模型更為精確，可較為全面地?cái)M合開采過程中地表下沉趨勢和最大下沉情況。

從圖7可看出，隨著工作面的推進(jìn)，地表巖層也隨之發(fā)生移動變形，并逐漸形成一個規(guī)模遠(yuǎn)大于采空區(qū)的下沉區(qū)域。分析地表走向發(fā)生明顯沉降（＞1 m）區(qū)域的各測點(diǎn)實(shí)測下沉速度隨工作面推進(jìn)的變化情況（圖8）可知，觀測線一側(cè)的地表在工作面推進(jìn)169～244 m和496 m處形成2個沉陷速度峰值區(qū)域，其中T14和T15測點(diǎn)的最大下沉速度在推進(jìn)距離169 m附近取得最大值，T15～T18測點(diǎn)最大下沉速度在推進(jìn)距離244 m附近取得最大值，并在T17點(diǎn)處出現(xiàn)地表最大下沉速度12.72 cm/d，以上各點(diǎn)在244 m后下沉速度開始迅速降低，表明開采推進(jìn)169～244 m時地表及工作面各上覆巖層處于劇烈的運(yùn)動當(dāng)中。采空區(qū)中部區(qū)域的地表各測點(diǎn)下沉速度則在推進(jìn)距離496 m左右達(dá)到最大值，形成第2個較為明顯的峰值區(qū)域，但低于首次出現(xiàn)峰值區(qū)域的下沉速度。地表走向下沉區(qū)域各點(diǎn)的沉降速度均呈現(xiàn)先上升后下降的周期性變化，且總體上距開切眼越近的點(diǎn)，其達(dá)到最大沉降速度的時間越早，開采中地表沉降速度最大值所在位置均滯后進(jìn)尺線一定距離，平均最大下沉速度滯后距為132 m。當(dāng)工作面推進(jìn)到169 m處時，地表下沉較為平緩，隨后地表開始急劇下沉，至推進(jìn)長度244 m時最大下沉值達(dá)到-7.745 m，地表移動變形異?；钴S，此后61601工作面長度中點(diǎn)（推進(jìn)312 m）之前的下沉邊界區(qū)各點(diǎn)下沉速度開始放緩并逐漸趨于穩(wěn)定，中點(diǎn)之后的地表各點(diǎn)下沉速度開始加快，但此后也不斷減慢，到開采作業(yè)結(jié)束前地表各點(diǎn)下沉速度均下降至趨近于零，地表總體移動活躍時間較短。造成這一特征的原因是由于快速推進(jìn)，使覆巖由下而上傳遞變形的時間縮短，推進(jìn)速度越快，上覆巖層越接近整體連續(xù)變形，動態(tài)變形過程相對縮短，壓實(shí)效果愈好，從而使移動期縮短。此外，通過現(xiàn)場實(shí)測發(fā)現(xiàn)，工作面切眼位置兩側(cè)下沉盆地邊緣出現(xiàn)動態(tài)裂縫，裂縫方向與采空區(qū)邊界方向基本一致，寬度一般200～300 mm，且以工作面切眼位置最為發(fā)育，這表明下沉盆地邊緣區(qū)域處于61601工作面裂縫帶中。工作面走向邊界角、移動角、裂隙角、超前影響角見表3。綜合分析可知，特厚煤層綜放開采下，地表沉陷呈現(xiàn)出盆地邊界陡峭、變形分布集中，地表下沉速度大且變化幅度大等特征，這是由于工作面快速推進(jìn)，煤層上覆巖層各層次的下沉速度普遍快速增加，相對懸露的時間減少，使得巖層發(fā)生形變、垮落并進(jìn)而影響地表結(jié)構(gòu)的范圍趨于集中。另外，特厚煤層綜放開采的采高大，巖體破碎、斷裂區(qū)域擴(kuò)大，形成的冒落帶、裂縫帶相應(yīng)增高，彎曲帶相對減小，也使地表移動變形分布相對集中。

運(yùn)用FLAC^3D模擬傾向主斷面處地表沉陷隨推進(jìn)長度的變化情況，并記錄不同推進(jìn)長度下地表最大下沉值（表4）。觀察可發(fā)現(xiàn)，推進(jìn)長度≤342 m時，地表下沉較不明顯，下沉值均低于1 m。當(dāng)工作面推進(jìn)到395 m時，地表最大下沉值首次＞1 m，開始進(jìn)入下沉活躍階段。繪制工作面推進(jìn)長度在244，342，395，496和624 m處工作面及其周圍（420～795 m）地表下沉曲線圖（圖9），可知開采結(jié)束后傾向主斷面地表沉陷同樣具有下沉速度大、下沉幅度劇烈、變形分布集中等特厚煤層綜放開采特征，且相比走向主斷面在距采空區(qū)中點(diǎn)上方地表的任意觀測點(diǎn)長度之間的傾斜程度更大，形成的下沉盆地邊緣也更陡峭。開采初期，采空區(qū)范圍較小且未到達(dá)傾向主斷面正上方，地表受影響發(fā)生形變的范圍及程度相對較小，傾向下沉現(xiàn)象直到工作面推進(jìn)過半以后才較為明顯，且隨著工作面的不斷推進(jìn)并最終到達(dá)停采線，下沉曲線的整體對稱性出現(xiàn)較為顯著的下降。傾向主斷面的下沉變化持續(xù)時間比走向主斷面更短，這是因?yàn)樵诠ぷ髅嫱撇蛇^程中，地表傾向主斷面上各點(diǎn)受采動影響較走向主斷面更晚。

采用FLAC^3D和Tecplot 360繪制開采區(qū)域Z向位移云圖和走向主斷面Z向位移云圖，如圖10，11所示。與圖7（a）和圖7（b）對比可知，工作面推進(jìn)196 m之前，煤層處于初始采動狀態(tài)，走向下沉區(qū)域呈現(xiàn)上部平緩，下部深陷的“V”字型斷面，且下沉點(diǎn)集中分布于最大下沉點(diǎn)兩側(cè)100 m的范圍內(nèi)，地表未形成明顯的下沉盆地。此后，采空區(qū)范圍不斷擴(kuò)展，工作面上覆巖層受應(yīng)力影響發(fā)生變形斷裂，并向上傳遞至地表，地表繼而快速進(jìn)入沉降活躍期，活躍期內(nèi)的移動與變形劇烈且集中，走向地表下沉區(qū)域進(jìn)一步擴(kuò)大并形成底部開闊的“U”字型斷面（圖10，圖11），下沉盆地逐漸發(fā)育。隨著工作面上覆巖層逐漸垮落并充填采空區(qū)，地表移動進(jìn)入衰減階段，開采至臨界開采尺寸時，地表達(dá)到充分采動狀態(tài)，此時開切眼附近的下沉量已基本保持穩(wěn)定，地表下沉盆地近似橢圓形，底部區(qū)域保持平坦并逐漸向兩側(cè)擴(kuò)展。

4 結(jié) 論

1）特厚煤層綜放開采地表下沉情況基本符合概率積分法預(yù)計(jì)模型，但其發(fā)生沉陷的區(qū)域范圍更小，沉陷區(qū)域邊界坡度更大；FLAC^3D數(shù)值模擬在下沉趨勢和下沉量方面與實(shí)測值吻合度更高，且能從地表下沉速度、下沉盆地發(fā)育狀況等動態(tài)變化量出發(fā)分析地表沉陷規(guī)律，對概率積分法地表沉陷預(yù)計(jì)模型具有較好地補(bǔ)充作用。

2）特厚煤層綜放開采具有開采厚度大、推進(jìn)速度快的特點(diǎn)，其地表沉陷情況相較一般開采呈現(xiàn)出下沉速度快、下沉幅度劇烈、變形分別集中的特征，形成的地表下沉盆地邊界陡峭、中部區(qū)域擴(kuò)大，地表最大下沉值達(dá)到-13.32 m，下沉系數(shù)0.53，主要影響角正切值2.9，主要影響半徑137.90 m。

3）地表走向主斷面各觀測點(diǎn)的沉降速度均先后出現(xiàn)先上升后下降的變化，并在觀測線一側(cè)取得2個較為明顯的峰值區(qū)域，最大下沉速度12.72 cm/d，最大下沉速度滯后距平均132 m。傾向主斷面在工作面中點(diǎn)上方地表的任意觀測點(diǎn)長度之間的傾斜程度比走向主斷面更大，下沉盆地邊緣更陡峭，下沉變化持續(xù)時間也更短。地表移動變形動態(tài)特征方面，工作面推進(jìn)196 m之前走向下沉區(qū)域呈現(xiàn)“V”字型斷面，且隨著開采逐漸擴(kuò)大，整體形成底部開闊的“U”字型斷面并進(jìn)一步發(fā)育為盆地，達(dá)到充分采動狀態(tài)時，地表下沉盆地近似橢圓形，底部區(qū)域平坦，并逐漸向兩側(cè)擴(kuò)展。

參考文獻(xiàn)（References）：

[1]中國礦業(yè)學(xué)院.煤礦巖層與地表移動[M].北京：煤炭工業(yè)出版社，1981.

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

[3]鄒友峰，鄧喀中，馬偉民.礦山開采沉陷工程[M].徐州：中國礦業(yè)大學(xué)出版社，2003.

[4]滕永海，唐志新，鄭志剛.綜采放頂煤地表沉陷規(guī)律研究及應(yīng)用[M].北京：煤炭工業(yè)出版社，2009.

[5]王樹元，趙德勤.礦山開采破壞學(xué)[M].北京：冶金工業(yè)出版社，1993.

[6]賀躍光.工程開挖引起的地表移動與變形模型及監(jiān)測技術(shù)研究[D].長沙：中南大學(xué)，2003.HE Yueguang.Modeling and monitoring of ground surface movements due to excavation[D].Changsha：Central South University，2003.

[7]LIOSIS N，MARPU P，PAVLOPOULOS K，et al.Ground subsidence monitoring with SAR interferometry techniques in the rural area of Al Wagan，UAE[J].Remote Sensing of Environment，2018，216：276-288.

[8]SALAMON M D G. Elastic analysis of displacements and stresses induced by the mining of seam or reef deposits[J].Journal of the South African Institute of Mining and Metallurgy，1963，64（4）：128-149.

[9]DONNELLY L，CRUZ H，ASMAR I，et al.The monitoring and prediction of mining subsidence in the Amaga，Angelopolis，Venecia and Bolombolo Regions，Antioquia，Colombia[J].Engineering Geology，2001，59（1-2）：103-114.

[10]BARYAKH A，TELEGINA E，SAMODELKINA N，et al.Prediction of the intensive surface subsidences in mining potash series[J].Journal of Mining Science，2005，41（4）：312-319.

[11]崔希民，鄧喀中.煤礦開采沉陷預(yù)計(jì)理論與方法研究評述[J].煤炭科學(xué)技術(shù)，2017，45（1）：160-169.CUI Ximin，DENG Kazhong.Research review of predicting theory and method for coal mining subsidence[J].Coal Science and Technology，2017，45（1）：160-169.

[12]朱廣軼，朱樂君，郭影.地表沉陷動態(tài)時間函數(shù)研究[J].西安科技大學(xué)學(xué)報(bào)，2009，29（3）：329-332.ZHU Guangyi，ZHU Lejun，GUO Ying.Research on dynamic time function of surface subsidence[J].Journal of Xian University of Science and Technology，2009，29（3）：329-332.

[13]余學(xué)義，王昭舜，楊云.大采深綜放開采地表移動變形規(guī)律[J].西安科技大學(xué)學(xué)報(bào)，2019，39（4）：555-563.YU Xueyi，WANG Zhaoshun，YANG Yun.Surface movement and deformation law of fully-mechanized caving mining with large mining depth[J].Journal of Xian University of Science and Technology，2019，39（4）：555-563.

[14]劉玉成.開采沉陷的動態(tài)過程及基于關(guān)鍵層理論的沉陷模型[D].重慶：重慶大學(xué)，2010.LIU Yucheng.Study on the dynamic course of the surface subsidence and the model based on theory of key rock stratum[D].Chongqing：Chongqing University，2010.

[15]石曉宇，張燕海，楊可明，等.2種灰色模型在礦區(qū)地表沉陷預(yù)計(jì)中的適用性[J].工礦自動化，2020，46（5）：28-33，81.SHI Xiaoyu，ZHANG Yanhai，YANG Keming，et al.Applicability of two grey models in earth-surface subsidence prediction of mining area[J].Journal of Mine Automation，2020，46（5）：28-33，81.

[16]黃明江，湯伏全.開采引起的地表動態(tài)下沉分形增長規(guī)律[J].西安科技大學(xué)學(xué)報(bào)，2009，29（2）：131-135.HUANG Mingjiang，TANG Fuquan.Fractal growth law of surface dynamic subsidence caused by mining[J].Journal of Xian University of Science and Technology，2009，29（2）：131-135.

[17]YU X，LV W，JIANG F，et al.Research on the automatic monitoring system for coal mining subsidence[J].Applied Mechanics and Materials，2014，644-650：1355-1360.

[18]王磊，張鮮妮，池深深，等.融合InSAR和GA的開采沉陷預(yù)計(jì)參數(shù)反演模型研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2018，43（11）：1635-1641.WANG Lei，ZHANG Xianni，CHI Shenshen，et al.Parameter inversion model for mining subsidence prediction based on fusion of InSAR and GA[J].Geomatics and Information Science of Wuhan University，2018，43（11）：1635-1641.

[19]趙林林，佟光成，申靈玲，等.基于ArcObjects的地表沉陷信息可視化技術(shù)研究與實(shí)現(xiàn)[J].中國礦業(yè)，2007（11）：33-35.ZHAO Linlin，TONG Guangcheng，SHEN Lingling，et al.Study and realization of mining subsidence information visualization based on ArcObjects[J].China Mining Magazine，2007（11）：33-35.

[20]亓立壯，安士凱，田超，等.基于無人機(jī)激光雷達(dá)技術(shù)的開采沉陷監(jiān)測方法與參數(shù)反演[J].科學(xué)技術(shù)與工程，2022，22（12）：4752-4761.QI Lizhuang，AN Shikai，TIAN Chao，et al.Mining subsidence monitoring method and parameter inversion based on unmanned aerial vehicle lidar technology[J].Science Technology and Engineering，2022，22（12）：4752-4761.

[21]郭文兵，李圣軍.“三軟”特厚煤層開采地表沉陷規(guī)律研究[J].河南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，33（6）：733-737.GUO Wenbing，LI Shengjun.Study pattern of surface subsidence characteristics of “three-soft” and especially thick coal seam mining[J].Journal of Henan Polytechnic University（Natural Science），2014，33（6）：733-737.

[22]胡青峰，崔希民，劉文鍇，等.特厚煤層重復(fù)開采覆巖與地表移動變形規(guī)律研究[J].采礦與巖層控制工程學(xué)報(bào)，2020，2（2）：31-39.HU Qingfeng，CUI Ximin，LIU Wenkai，et al.Law of overburden and surface movement and deformation due to mining super thick coal seam[J].Journal of Mining and Strata Control Engineering，2020，2（2）：31-39.

[23]高超，徐乃忠，劉貴，等.特厚煤層綜放開采地表沉陷規(guī)律實(shí)測研究[J].煤炭科學(xué)技術(shù)，2014，42（12）：106-109.GAO Chao，XU Naizhong，LIU Gui，et al.Research on actual measurement of surface subsidence law for fully-mechanized top coal caving in extra thick seam[J].Coal Science and Technology，2014，42（12）：106-109.

[24]趙兵朝，劉陽，賀衛(wèi)中，等.郭家河煤礦綜放開采地表巖移參數(shù)分析[J].煤礦安全，2018，49（8）：254-257.ZHAO Bingchao，LIU Yang，HE Weizhong，et al.Parameters analysis of surface rock movement in Guojiahe mine fully mechanized caving mining[J].Safety in Coal Mines，2018，49（8）：254-257.

[25]譚志祥，王宗勝，李運(yùn)江，等.高強(qiáng)度綜放開采地表沉陷規(guī)律實(shí)測研究[J].采礦與安全工程報(bào)，2008，25（1）：59-62.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-62.

[26]郭廣禮，查劍峰.礦山開采沉陷學(xué)[M].徐州：中國礦業(yè)大學(xué)出版社，2020.

[27]劉俊.急傾斜特厚煤層圍巖應(yīng)力演化規(guī)律數(shù)值模擬研究[J].煤炭技術(shù)，2021，40（8）：76-80.LIU Jun.Numerical simulation research on stress evolution law of surrounding rock in steeply inclined and extra-thick coal seam[J].Coal Technology，2021，40（8）：76-80.

[28]王哲，涂敏，李利青.小回溝煤礦基于正交試驗(yàn)設(shè)計(jì)的巖體力學(xué)參數(shù)反演[J].采礦與安全工程學(xué)報(bào)，2018，35（2）：397-401.WANG Zhe，TU Min，LI Liqing.Back analysis of rock mechanical parameters in Xiaohuigou coal mine based on orthogonal experimental design[J].Journal of Mining and Safety Engineering，2018，35（2）：397-401.