准格尔煤田龙王沟煤矿特厚煤层综放开采地表沉陷规律

张海君 吴奕枢 王飞 谢昆仑

摘 要：在特厚煤层综放开采过程中，工作面各上覆岩层受其影响相继发生变形、垮落，使地表出现不均匀沉降，产生塌陷、裂缝等，严重威胁矿井及其周边地区人员和生产设施安全。因此，地表沉陷变化规律的研究对沉陷区范围预测和治理均有重要意义。文中以准格尔煤田龙王沟煤矿61601特厚煤层综放工作面为研究对象，建立基于概率积分法的地表沉陷预计模型，并对计算结果加以分析，结果显示特厚煤层综放开采地表下沉情况基本符合概率积分法预计模型，但其发生沉陷的區域范围更小，沉陷区域边界坡度更大，特厚煤层综放开采地表沉陷情况相较一般开采呈现出下沉速度快、下沉幅度剧烈、变形分别集中的特征。同时运用FLAC^3D数值模拟技术对地表沉陷形成机制和动态变化过程进行补充，FLAC^3D数值模拟能更为全面地体现地表沉陷规律。此外，文中还进一步结合FLAC^3D数值模拟和实测沉陷结果分析了不同推进长度下走向和倾向主断面地表下沉曲线、走向地面测点下沉速度变化情况、开采区域总体位移情况等，归纳出特厚煤层综放开采地表沉陷规律及特征，为合理规划煤矿开采作业方案提供理论依据。关键词：特厚煤层；地表沉陷；FLAC^3D；概率积分法中图分类号：TD 822

文献标志码：A

文章编号：1672-9315（2022）05-0874-10

DOI：10.13800/j.cnki.xakjdxxb.2022.0505开放科学（资源服务）标识码（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 引 言井下煤炭资源采出后，煤层顶板悬露，顶板原岩应力遭到扰动，顶板岩石在重力作用下将出现移动变形。随着工作面的继续推进，顶板出现破坏、冒落，导致地表沉陷现象的发生，将会对矿区及矿区周边的建筑物、生产设施和人员的安全造成威胁。随着中国煤炭资源的持续勘探和开发，煤炭开采所面对的矿区环境更加复杂，如何有效预测和防范开采造成的地表沉降危害，已成为矿区整体规划设计的一个重点。井下开采所致的工作面顶板岩移及地表沉陷现象是一个涉及复杂岩体力学的变化过程。以往国内外学者提出了许多地表沉陷模型和理论^[1-4]，典型的如基于观测数据的正态理论、二等分线理论、自然斜面理论和开采沉陷“拱理论”^[5-6]、工作面上覆岩层“三带”分布理论^[7]等。随着科学技术水平的不断提升，开采沉陷理论得到进一步发展。Salamon将影响函数预计法与连续介质理论相结合，开创了用于探究开采沉陷规律的边界元法（BEM）^[8]；DONNELLY等基于SWIFT技术对矿山地表沉陷问题进行预测与分析^[9]，BARYAKH等基于线性粘弹性理论，利用蠕变函数开发了一种动态地表沉陷预计程序^[10]；崔希民等应用流变模型对地表沉降规律进行分析^[11]；朱广轶等基于概率积分法提出了一种动态地表沉陷时间函数，并通过实例分析证明其具有比Knote函数更优的精准度^[12]；余学义等根据观测数据研究单一工作面开采和多工作面开采后地表的沉陷规律，并分别评估了2种开采场景下地表的损害情况^[13]；刘玉成对Knote时间函数加以改良，建立了一种更加精准的地表沉陷预计模型^[14]；石晓宇等探讨了2种灰色模型在地表沉陷监测与预计方面的应用，并指出二者适用的开采情况^[15]；黄明江等指出地表下沉曲线符合分形增长规律，并借助分形插值法实现了有限观测数据下对地表沉陷连续变化情况的拟合^[16]。此外，一些学者还采用GIS，InSAR，ArcObjects，RTK，无人机激光雷达等新型观测技术和地理信息系统对矿山地表沉陷情况展开监测与分析^[17-20]。中国许多大型矿区属于厚煤层，煤炭储量相当丰富。近年来，综放开采技术被运用于厚煤层的开采过程中，因其具有产量大、效率高、成本低等特点，而备受业内认可。关于厚煤层开采引起的覆岩和地表移动规律等方面目前已有不少研究成果，如郭文兵等对“三软”特厚煤层地表沉陷情况进行观测，并基于概率积分法运用Matlab拟合了地表下沉曲线^[21]；胡青峰等建立了塔山煤矿特厚煤层倾向主断面相似模拟实验，并研究了重复采动下覆岩变形和地表沉陷规律^[22]；高超等针对东坡煤矿浅埋深、特厚煤层综采等特点，结合地表实测数据，对地表动态移动变形特征进行分析^[23]；赵兵朝等通过对郭家河煤矿特厚煤层综放开采实地观测和分析求取了地表移动变形有关参数，并综合运用反演模拟和边坡稳定性分析对所得参数进行修正^[24]，但所研究的对象较少涉及20 m以上的大采放比特厚煤层综放开采，故仍需对此类作业引起的覆岩与地表移动规律进行系统地研究。

文中将根据准格尔煤田龙王沟煤矿特厚煤层一次采全厚这一开采特征，建立基于概率积分法的地表沉陷预计模型，运用该模型分析煤层开采引起的地表移动规律与变形特征，同时讨论该模型存在的局限性，并采用FLAC^3D数值模拟方法加以补充，从而获得较为准确、全面的特厚煤层综放开采的地表沉陷规律及其特征，为合理规划煤矿开采作业方案，保护矿区人员生产生活的安全提供指导方向和理论支持。

1 龙王沟煤矿61601工作面概述龙王沟井田位于内蒙古自治区鄂尔多斯市准格尔煤田中北部，面积41.35 km²，属典型黄土高原地貌。依據地表出露及钻孔揭露情况，该井田地层自下而上分别为：奥陶系中统马家沟组（O₂m），石炭系上统太原组（C₂t），二叠系下统山西组（P₁s）、下石盒子组（P₁x），二叠系上统上石盒子组（P₂s）、石千峰组（P₂sh），三叠系下统刘家沟组（T₁l），第三系上新统（N₂）和第四系（Q）等。煤田地层结构形态为一近似于南北走向、呈波状起伏、向西侧倾斜的单斜构造，地层倾角较小，低于10°。61601工作面布置于龙王沟煤田太原组上部6煤处，为首采面，平均采深399.9 m，工作面倾向长度254.6 m，走向可采长度624 m，煤均厚23.15 m，煤层稳定。工作面设计采高5.1 m，顶煤厚18.05 m，采放比1∶3.54。煤层倾角0～2°，接近水平煤层。煤层顶底板各岩层如图1所示。61601工作面开采选用后退式走向长壁采煤方法，综采放顶煤一次采全高开采工艺，完全垮落法管理顶板。为获得特厚煤层开采条件下地表沉陷的实测数据，分析地表沉陷规律，根据《煤矿测量规程》在61601工作面外设置地表沉陷观测站，根据《煤矿测量规程》规定，当工作面走向长度大于0.9H₀（H₀为煤层埋深，399.9 m）时，可只设半条走向观测线，故本工作面设置半条走向观测线。观测点按便于观测塌陷区的下沉和能够准确反映工作面下沉情况进行布置，并以观测站为起点布置观测点，首个观测点距开切眼上方地表水平距离288 m。之后每隔25 m设置一个观测点，直至工作面走向推进长度中点之后25 m，以确保观测数

据能够监测到下沉量的极值点。观测线全长625 m，

各点标记为T1，T2，T3，…，T26。各测点的布置如图2所示，图中矩形部分为61601工作面示意图，三角形中心为观测点位置。走向测线下沉测量工作起始于2019年1月17日，测量周期1年5个月，各测点平均每月测量一次，直至测量完毕。

2 基于概率积分法的地表沉陷预计模型概率积分法是基于随机介质理论发展而来的岩体移动与变形计算模型，经典的概率积分法表述的地表任意点下沉位移计算公式见式（1）。

公式（1）表示地表任意点（x，y）在开采区域为（x₁-x₀）*（y₁-y₀）的煤层开采后产生的下沉量。W_max为地表最大下沉量；x₀，x₁，y₀，y₁分别是煤层沿倾向开采和走向开采的起止点；r为主要影响半径，可由公式（2）得到。

r=Htanβ（2）

式中 H为采深；tanβ为主要影响角正切值。文中使用Python Integrate模块下的Dblquad函数编写基于概率积分法的地表任意点移动与变形模拟程序，设地表为一倾向长945 m，走向长1 392 m的水平平面，并取左下角为坐标原点（0，0）。平面采用矩形网格进行划分，依照61601工作面实测数据确定开采范围和岩移参数，煤层厚23.15 m，下沉系数0.53，主要影响角正切值2.9。地表下倾向开采区域位于345～600 m之间，走向开采区域位于384～1 008 m之间。

模拟结果将保存为记录有地表各节点坐标及下沉量的数据文件。将该文件导入Surfer中，绘制地表开采下沉量等值线图，如图3所示，同时比较计算值与实测值，绘制相应的地表下沉曲线图，如图4所示（部分观测点受采动影响发生破坏，这些点未在图中绘出）。

比照实际观测成果可知，61601工作面地表下沉变形趋势总体上与概率积分法沉陷预计模型相吻合。从图3可看出，采动后地表出现了一个以工作面所在位置上方为中心并向外扩大的类椭圆形沉降区域。地表下沉值较大，达到-13.32 m，等值线在-2.5～-10 m之间密集分布，显示此区间内地表高程存在较大幅度的下移变化，形成了坡度较陡的地表沉降区域边界，下沉量大于-10 m的等值线则较为稀疏，表明地表形成了上宽下窄的“漏斗状”下沉盆地。从模型的岩移参数上看，相比普通煤层的开采^[25-26]，特厚煤层综放开采的地表下沉系数较大，边界角值比较小，主要影响角正切值显著偏大，反映出大采高、快速推进的综采放顶煤开采地表下沉速度快，岩体动态变形幅度大、运动剧烈的特征。此外，对比发现实测值具有更为陡峭的地表下沉边界，更为集中的地表下沉区域，实测地表发生沉陷的范围更小，下沉曲线在盆地边界收敛更快，充分地体现出特厚煤层综放开采地表沉陷相比一般地表沉陷规律具有下沉幅度剧烈、变形集中等特点。概率积分法具有公式简单、计算高效快捷的特点，其只需借助较少的岩移参数，就能计算地表沉陷情况，为研究地质条件简单、岩层达到稳定状态的工作面沉陷状况提供了快速预计和分析的方法。但是概率积分法是一种基于随机介质理论的理想化模型，其不能充分考虑开采区各岩层的空间构型、力学特征的差异对地表移动

变形产生的影响；同时作为一种静态地表沉陷预计方法，概率积分法难以直接建立地表移动变化趋势与进尺长度、时间等因素之间的关系等。为提高地表沉陷预计模型的精准性，进一步解释实地观测和概率积分法模型预计结果的形成机理，从动态的角度探究开采沉陷随时间的变化规律，这里讨论采用FLAC^3D数值模拟方法对地表沉陷预计模型进行补充。

3 基于FLAC^3D的61601工作面地表沉陷规律模拟

3.1 FLAC^3D数值模拟模型的建立根据61601工作面相关地质资料建立几何模型，如图5所示。该模型沿X轴方向全长1 195.4 m，對应开采工作面的倾向方向；沿Y轴方向全长2 000 m，对应开采工作面的走向方向。模型依据各地层实地勘测到的岩土组成情况分别建立了28个地层，并根据地质结构添加部分节理接触面以提高精度，总深度417.8 m。其中开采工作面位于第24岩层（六煤）的下部，开采区域所在位置为：沿X轴方向0～254.6 m，沿Y轴方向0～624 m。模型选用矩形网格进行划分，根据矿区岩层倾角很小这一特征将模型岩层形态设置为水平岩层。

以往的研究表明，直接将岩土力学试验获得的力学参数应用于数值模拟当中往往容易发生错误，通常需取实际力学参数的1/3～1/5作为数值模拟的模拟值，并进行反演分析多次计算方可取得较为准确的模型力学参数^[27-28]。根据模拟比例经验值运用反演计算确定岩层模型最优力学参数，并导入模型中，主要岩层模型的力学参数见表1。

本次模拟基于摩尔-库伦弹塑性模型进行预测计算，分别记录工作面推进169 m（开挖第37天）、244 m（开挖第77天）、342 m（开挖第181天）、496 m（开挖第310天）和624 m（开挖第369天）时地表的走向下沉情况。开挖模型如图6所示，其中2个三角形的位置对应地表下沉观测点的起止范围（T1～T26），模型下部的矩形代表61601工作面投影，模型开挖进度与61601工作面实际开挖进度保持一致。

3.2 结果及分析选取工作面进尺169，244，342，496和624 m时的走向地表Z向位移拟合数据，与地表观测站实测下沉量进行比照（表2），绘制下沉曲线瀑布图，如图7所示（部分观测点受采动影响发生破坏掩埋无法记录，这些点未在图中绘出）。

从图7（e）数据可看出，开采完毕后FLAC^3D地表沉陷预计模型的最大下沉量为-13.39m，与实测值（-13.32 m）相差小于5%，可见FLAC^3D数值模拟在最大下沉量方面显示出了很好的拟合度。实测的地表最大下沉点位于距观测线起点525 m的T22号观测点附近，而数值模拟的最大下沉点位于地表平面的中心点（T25点），这可能是因为自然条件下矿区岩体仍存在一定倾角以及岩层间存在断层等不连续结构所致。此外，与模拟值相比，开采初期（图7（a）和（b））实测下沉曲线在拐点后的切线斜率迅速增大，形成更陡峭的下沉曲线，但随着工作面的推进，模拟下沉曲线与实测数据间的差距持续缩小（图7（c）和（d）），至停采线位置时（图7（e）），二者已基本趋于一致。总体而言，FLAC^3D对特厚煤层综放开采地表沉陷的数模模拟效果相比经典概率积分法预计模型更为精确，可较为全面地拟合开采过程中地表下沉趋势和最大下沉情况。

从图7可看出，随着工作面的推进，地表岩层也随之发生移动变形，并逐渐形成一个规模远大于采空区的下沉区域。分析地表走向发生明显沉降（＞1 m）区域的各测点实测下沉速度随工作面推进的变化情况（图8）可知，观测线一侧的地表在工作面推进169～244 m和496 m处形成2个沉陷速度峰值区域，其中T14和T15测点的最大下沉速度在推进距离169 m附近取得最大值，T15～T18测点最大下沉速度在推进距离244 m附近取得最大值，并在T17点处出现地表最大下沉速度12.72 cm/d，以上各点在244 m后下沉速度开始迅速降低，表明开采推进169～244 m时地表及工作面各上覆岩层处于剧烈的运动当中。采空区中部区域的地表各测点下沉速度则在推进距离496 m左右达到最大值，形成第2个较为明显的峰值区域，但低于首次出现峰值区域的下沉速度。地表走向下沉区域各点的沉降速度均呈现先上升后下降的周期性变化，且总体上距开切眼越近的点，其达到最大沉降速度的时间越早，开采中地表沉降速度最大值所在位置均滞后进尺线一定距离，平均最大下沉速度滯后距为132 m。当工作面推进到169 m处时，地表下沉较为平缓，随后地表开始急剧下沉，至推进长度244 m时最大下沉值达到-7.745 m，地表移动变形异常活跃，此后61601工作面长度中点（推进312 m）之前的下沉边界区各点下沉速度开始放缓并逐渐趋于稳定，中点之后的地表各点下沉速度开始加快，但此后也不断减慢，到开采作业结束前地表各点下沉速度均下降至趋近于零，地表总体移动活跃时间较短。造成这一特征的原因是由于快速推进，使覆岩由下而上传递变形的时间缩短，推进速度越快，上覆岩层越接近整体连续变形，动态变形过程相对缩短，压实效果愈好，从而使移动期缩短。此外，通过现场实测发现，工作面切眼位置两侧下沉盆地边缘出现动态裂缝，裂缝方向与采空区边界方向基本一致，宽度一般200～300 mm，且以工作面切眼位置最为发育，这表明下沉盆地边缘区域处于61601工作面裂缝带中。工作面走向边界角、移动角、裂隙角、超前影响角见表3。综合分析可知，特厚煤层综放开采下，地表沉陷呈现出盆地边界陡峭、变形分布集中，地表下沉速度大且变化幅度大等特征，这是由于工作面快速推进，煤层上覆岩层各层次的下沉速度普遍快速增加，相对悬露的时间减少，使得岩层发生形变、垮落并进而影响地表结构的范围趋于集中。另外，特厚煤层综放开采的采高大，岩体破碎、断裂区域扩大，形成的冒落带、裂缝带相应增高，弯曲带相对减小，也使地表移动变形分布相对集中。

运用FLAC^3D模拟倾向主断面处地表沉陷随推进长度的变化情况，并记录不同推进长度下地表最大下沉值（表4）。观察可发现，推进长度≤342 m时，地表下沉较不明显，下沉值均低于1 m。当工作面推进到395 m时，地表最大下沉值首次＞1 m，开始进入下沉活跃阶段。绘制工作面推进长度在244，342，395，496和624 m处工作面及其周围（420～795 m）地表下沉曲线图（图9），可知开采结束后倾向主断面地表沉陷同样具有下沉速度大、下沉幅度剧烈、变形分布集中等特厚煤层综放开采特征，且相比走向主断面在距采空区中点上方地表的任意观测点长度之间的倾斜程度更大，形成的下沉盆地边缘也更陡峭。开采初期，采空区范围较小且未到达倾向主断面正上方，地表受影响发生形变的范围及程度相对较小，倾向下沉现象直到工作面推进过半以后才较为明显，且随着工作面的不断推进并最终到达停采线，下沉曲线的整体对称性出现较为显著的下降。倾向主断面的下沉变化持续时间比走向主断面更短，这是因为在工作面推采过程中，地表倾向主断面上各点受采动影响较走向主断面更晚。

采用FLAC^3D和Tecplot 360绘制开采区域Z向位移云图和走向主断面Z向位移云图，如图10，11所示。与图7（a）和图7（b）对比可知，工作面推进196 m之前，煤层处于初始采动状态，走向下沉区域呈现上部平缓，下部深陷的“V”字型断面，且下沉点集中分布于最大下沉点两侧100 m的范围内，地表未形成明显的下沉盆地。此后，采空区范围不断扩展，工作面上覆岩层受应力影响发生变形断裂，并向上传递至地表，地表继而快速进入沉降活跃期，活跃期内的移动与变形剧烈且集中，走向地表下沉区域进一步扩大并形成底部开阔的“U”字型断面（图10，图11），下沉盆地逐渐发育。随着工作面上覆岩层逐渐垮落并充填采空区，地表移动进入衰减阶段，开采至临界开采尺寸时，地表达到充分采动状态，此时开切眼附近的下沉量已基本保持稳定，地表下沉盆地近似椭圆形，底部区域保持平坦并逐渐向两侧扩展。

4 结 论

1）特厚煤层综放开采地表下沉情况基本符合概率积分法预计模型，但其发生沉陷的区域范围更小，沉陷区域边界坡度更大；FLAC^3D数值模拟在下沉趋势和下沉量方面与实测值吻合度更高，且能从地表下沉速度、下沉盆地发育状况等动态变化量出发分析地表沉陷规律，对概率积分法地表沉陷预计模型具有较好地补充作用。

2）特厚煤层综放开采具有开采厚度大、推进速度快的特点，其地表沉陷情况相较一般开采呈现出下沉速度快、下沉幅度剧烈、变形分别集中的特征，形成的地表下沉盆地边界陡峭、中部区域扩大，地表最大下沉值达到-13.32 m，下沉系数0.53，主要影响角正切值2.9，主要影响半径137.90 m。

3）地表走向主断面各观测点的沉降速度均先后出现先上升后下降的变化，并在观测线一侧取得2个较为明显的峰值区域，最大下沉速度12.72 cm/d，最大下沉速度滞后距平均132 m。倾向主断面在工作面中点上方地表的任意观测点长度之间的倾斜程度比走向主断面更大，下沉盆地边缘更陡峭，下沉变化持续时间也更短。地表移动变形动态特征方面，工作面推进196 m之前走向下沉区域呈现“V”字型断面，且随着开采逐渐扩大，整体形成底部开阔的“U”字型断面并进一步发育为盆地，达到充分采动状态时，地表下沉盆地近似椭圆形，底部区域平坦，并逐渐向两侧扩展。

参考文献（References）：

[1]中国矿业学院.煤矿岩层与地表移动[M].北京：煤炭工业出版社，1981.

[2]刘宝琛，廖国华.煤矿地表移动的基本规律[M].北京：中国工业出版社，1965.

[3]邹友峰，邓喀中，马伟民.矿山开采沉陷工程[M].徐州：中国矿业大学出版社，2003.

[4]滕永海，唐志新，郑志刚.综采放顶煤地表沉陷规律研究及应用[M].北京：煤炭工业出版社，2009.

[5]王树元，赵德勤.矿山开采破坏学[M].北京：冶金工业出版社，1993.

[6]贺跃光.工程开挖引起的地表移动与变形模型及监测技术研究[D].长沙：中南大学，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]崔希民，鄧喀中.煤矿开采沉陷预计理论与方法研究评述[J].煤炭科学技术，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]朱广轶，朱乐君，郭影.地表沉陷动态时间函数研究[J].西安科技大学学报，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]余学义，王昭舜，杨云.大采深综放开采地表移动变形规律[J].西安科技大学学报，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]刘玉成.开采沉陷的动态过程及基于关键层理论的沉陷模型[D].重庆：重庆大学，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种灰色模型在矿区地表沉陷预计中的适用性[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]黄明江，汤伏全.开采引起的地表动态下沉分形增长规律[J].西安科技大学学报，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的开采沉陷预计参数反演模型研究[J].武汉大学学报（信息科学版），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的地表沉陷信息可视化技术研究与实现[J].中国矿业，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].科学技术与工程，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]郭文兵，李圣军.“三软”特厚煤层开采地表沉陷规律研究[J].河南理工大学学报（自然科学版），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]胡青峰，崔希民，劉文锴，等.特厚煤层重复开采覆岩与地表移动变形规律研究[J].采矿与岩层控制工程学报，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]高超，徐乃忠，刘贵，等.特厚煤层综放开采地表沉陷规律实测研究[J].煤炭科学技术，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]赵兵朝，刘阳，贺卫中，等.郭家河煤矿综放开采地表岩移参数分析[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]谭志祥，王宗胜，李运江，等.高强度综放开采地表沉陷规律实测研究[J].采矿与安全工程报，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]郭广礼，查剑峰.矿山开采沉陷学[M].徐州：中国矿业大学出版社，2020.

[27]刘俊.急倾斜特厚煤层围岩应力演化规律数值模拟研究[J].煤炭技术，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]王哲，涂敏，李利青.小回沟煤矿基于正交试验设计的岩体力学参数反演[J].采矿与安全工程学报，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.