煤矸石井下处理与冒落区流态化充填技术

李亮 黄庆享 吴杰 左小 张伟龙 张杰 胡俊峰 杨玉玉

摘 要：為解决韩城矿区煤矸石排放引起的制约生产、生态环境恶化的矛盾，在总结煤矸石井下处理与利用方法的基础上，系统阐述了冒落区流态化充填技术原理和优势，采用理论计算和瞬变电磁法探查方法，揭示了试验工作面冒落区残余空间分布规律，确定了基于“高、低位协同共充”和“空洞-空隙-孔隙”多类型残余空间精准高效利用的冒落区流态化充填方式，并在桑树坪二号井实施了工业性试验，结果表明：试验累计充填矸石浆体4 949 m³，处理矸石6 414 t，验证了冒落区流态化充填矸石的可行性，发现了矸石浆体在冒落区扩散过程中存在显著的“通道效应”，采用“高、低位协同共充”技术可显著降低或解除“通道效应”和高效利用冒落区残余空间，但该技术增大了钻孔工程量，增加了充填成本，因此在工程尺度上合理确定高、低位充填钻孔参数提出要求，进而确保达到较好的充填效果。研究成果为长壁开采冒落区矸石流态化充填技术提供理论支撑。

关键词：煤矸石；冒落区；残余空间；流态化充填；浆体；扩散规律

中图分类号：TD 823

文献标志码：A

文章编号：1672-9315（2022）05-0865-09

DOI：10.13800/j.cnki.xakjdxxb.2022.0504开放科学（资源服务）标识码（OSID）：

Underground treatment of coal gangue and fluidization filling technology in caving area

LI Liang^1，2，3，HUNAG Qingxiang¹，WU Jie⁴，ZUO Xiao^2，3，ZHANG Weilong^1，2，3，ZHANG Jie⁴，HU Junfeng⁴，YANG Yuyu^2，3

（1.College of Energy Science and Engineering，Xian University of Science and Technology，Xian 710054，China；2.Shaanxi Coal and Chemical Industry Technology Research Institute Co.，Ltd.，Xian 710100，China；3.The National Joint Engineering Research Center of the Green，Safe and Efficient Coal Mining，Xian 710065，China；4.Shaanxi Shanmei Hancheng Mining Co.，Ltd.，Hancheng 715400，China）

Abstract：In order to solve the contradiction of restricting production and worsening ecological environment caused by the discharge of coal gangue in Hancheng mining area，

the principle and advantages of fluidized filling technology in caving area are systematically expounded by summarizing  the underground treatment and utilization methods of coal gangue.The theoretical  calculation and transient electromagnetic method are applied to reveal the residual space distribution law of the caving area in the test working face，with the fluidization filling formula of the caving area  determined，which is based on the “high and low collaborative filling” and the multi type residual space accurate efficient use of “cavity-void-pore”.An industrial test was carried out in No.2 well of Sang Shuping.The results show that the test accumulated 4 949 m³of gangue slurry and 6 414 t of gangue，which verifies the feasibility of fluidized gangue filling in the caving area.It is found that there is a significant “channel effect” in the diffusion process of gangue slurry in the caving area，and the “high and low collaborative filling” technology can significantly reduce or eliminate the “channel effect” and help to make efficient use of the residual space in the caving area，but it increases the drilling quantities and the filling cost.Therefore，it puts forward requirements for reasonably determining the drilling parameters of high and low level filling in the engineering scale，so as to ensure a better filling effect.The research results provide theoretical support for the fluidization filling technology of gangue in the caving area of longwall mining.

Key words：coal gangue；caving area；residual space；fluidization filling；slurry；diffusion law

0 引 言近年来，随着大型现代化矿井建设的推进，煤矿排矸呈集中化、高产化和规模化的发展趋势，开采产生的矸石堆放至地面引发的环境污染问题愈发突出^[1]。据不完全统计，目前中国矸石累计堆放量超过60亿t，形成矸石山1 500～1 700座，占地20余万亩，且以约5亿～8亿t/a的排放量逐年增加^[2]。韓城矿区的煤矸石年排放量已突破100万t，绝大多数煤矸石仍以地面堆积的方式进行处理，不但制约生产，而且煤矸石地表堆积对矿区生存环境和条件带来的威胁与危害愈加严重。在进入“十四五”时期，能源供需格局深刻变化，《煤炭工业“十四五”高质量发展指导意见》等文件中均指出提高煤矸石综合利用水平，推动煤炭绿色低碳发展。不设地面排矸场、减少或杜绝矸石升井是发展低碳经济的必然趋势。众多学者和煤矿企业在源头控制矸石排放、保护矿区生态环境方面亦做了大量的工作，从矸石处理与利用的众多实践中得出，现有以煤矸石发电、铺路、制建筑材料等为主的煤矸石地面综合处理技术已经不能满足矿区日益增大的煤矸石处理量的要求，仅能作为煤矸石规模化处理的补充技术^[2-5]。煤矸石井下充填处理技术具有处理效率高、产矸点就是处理点、处理集中等优点，以综合机械化固体充填技术和胶结充填技术为核心的处理技术应用较为广泛，但受充填效果要求、装备工艺及自动化程度等方面的限制，也存在投资大、成本高、工艺复杂、处理能力小等缺点，主要用于“三下”采煤^[6-11]。神东矿区已在15个矿井建成地下水库35座，总设计储水量达2 499.5万m³，储水系数一般为 0.15～0.25，这对地下采动空间进行了成功的探索和实践^[12]。冒落区矸石流态化充填技术是地下采动空间储废新技术的一种，是将矸石破碎后简单胶结或直接加水制浆泵送充填至井下已冒落采空区的新型充填技术，既满足矿井矸石生态处理的要求，又对采空区灾害治理有一定促进作用，充填工艺简单、投资小、“采—充”工艺独立，在低成本规模化处理矸石方向的应用前景已初步显现。钱鸣高院士等通过开采后岩层移动实测，得到了岩层内部位移曲线，并揭示了采场上覆岩层下沉规律和采空区覆岩破断的“O型圈”形态，为精准利用采空区空间奠定了基础^[13]。李兴尚利用分形理论研究条带开采冒落矸石堆体形态、空隙分布特征，并对某煤矿冒落区注浆充填量进行了预计^[14]。张宏贞通过建立覆岩空间结构类型的残留空洞、空隙预测模型，将采空区覆岩结构类型分为“匚”、“工”、“伞”型3种基本类型^[15]。胡胜勇等研究了采空区应力分布、岩石碎胀特性和垮落带空隙率的关系^[16]。黄炳香等通过相似模拟试验研究了冒落带和裂隙带的孔隙率变化规律^[17]。王建学等研究了冒落矸石空隙注浆胶结充填减沉技术的可行性^[18]。张天军等研究了矸石膏体充填材料粒径对其力学性能的影响，根据Talbol理论进行连续级配设计，配制了矸石膏体浆料^[19]。古文哲等提出了煤矿固体废弃物浆体处理的技术构思，在龙王沟煤矿和锦界电厂进行了应用^[20]。上述研究，为矸石冒落区流态化充填技术的应用提供了良好借鉴，但是，对地下采动空间的产生与发展规律的研究及利用尚无系统性成果。基于砌体梁理论和采空区上覆岩层下沉规律，提出利用工作面冒落采空区残余空间处理矿井矸石，研究冒落区残余空间分布规律和冒落区矸石流态化充填技术，并在桑树坪二号井开展了工业性试验。

1 煤矸石井下处理与利用

1.1 煤矸石井下处理与利用方法西部矿区煤炭开采普遍具有采高大、单产高等高强度开采特点，造成采动损伤更加剧烈，加之西部矿区生态脆弱，面临高强度开采和生态损伤剧烈与修复难的突出矛盾^[21-23]。韩城矿区位于渭北煤田东北端，属典型的渭北黄土高原低山丘陵地貌特征，矸石地面堆积处理同样产生以上问题，煤矸石的综合利用及生态治理已是一个刻不容缓、亟待解决的问题，因此，韩城矿区煤矸石处理与利用应遵循减量化排放、源头治理、“采—充”平衡、经济高效和环境效益最佳的原则。基于以上原则，矿区矸石井下处理与利用方法如图1所示。

1.2 冒落区矸石流态化充填技术原理在煤层开采后，覆岩由下至上分别为冒落带、裂隙带、弯曲下沉带。基于采动损伤空间守恒原理，提出采动损伤空间耗散分区模型^[24]如下。

V_M=V₀+V₁+V₂+V₃    （1）

式中 V₀为地表下沉的总空间；V₁为冒落带采动损伤耗散空间；V₂为裂隙带内采动损伤耗散空间；V₃为弯曲下沉带内损伤耗散空间。根据空间守恒原理，采动损伤总空间V_M等于煤层采出空间。冒落带由冒落矸石组成，即使工作面充分采动后，内部仍然残余较大的空间，加之破碎煤岩体自身的碎胀性，使得冒落带为采动后空间最大的区域，这些空间是随着工作面推采动态变化的，且以空洞、空隙、孔隙等多种类型存在，称之为冒落区残余空间。在采空区边界附近由于悬臂岩梁的存在会产生大的空洞区域，深部则主要以空隙、孔隙区域分布，并且距离开采边界越远，空洞、空隙、孔隙量则越小^[14-15]。煤层开采后产生大量的冒落区残余空间可作全矿井矸石处理的天然空间，因此，可通过邻近巷道在冒落带范围内协同布置高位、低位充填孔进行矸石流态化充填，实现冒落区“空洞-空隙-孔隙”多类型残余空间精准高效利用，即“高、低位协同共充技术”。高位孔终孔高度一般以冒落带上限为准，水平孔终孔高度与采高一致，具体可根据地质采矿条件调整。技术原理如图2所示。

冒落区流态化充填系统布置地面或井下矸石仓、破碎站、制浆站、充填泵站及管路，如图3所示。制备好的浆体通过泵站和管路输送至充填作业点，利用高位钻孔和低位钻孔协同共充至冒落区。

2 冒落区残余空间分布规律

2.1 试验工作面概况桑树坪二号井3305工作面采用综合机械化放顶煤开采工艺，煤厚6 m，整体倾角为1°，工作面沿煤层倾向布置，切眼宽度165 m，可采长度870 m，可采面积143 550 m²。顶底板概况见表1。

2.2 冒落区残余空间理论计算冒落区内部的残余空间通常是不规则的、是无序分布的。根据冒落岩块分维数与孔隙率的关系，得到孔隙率与分维数的关系式为^[14]

式中 P为孔隙率，%；V为冒落区岩块总体积，m³；M为冒落矸石岩块质量，kg；ρ₀为冒落岩块的密度，kg/m³；r_min为岩层冒落堆体中破碎最严重、数目最多，尺寸最小的岩块的直径，m；ρ为冒落区松散堆体密度，kg/m³；D为分维数；H为垮落带高度，m；h为开采煤层高度，m；b为开采煤层倾向长度，m；L为开采煤层走向长度，m。根据试验工作面地层条件，结合冒落区矸石取样观察，工作面伪顶及直接顶较易垮落，冒落区矸石中泥岩、砂质泥岩、细砂岩的破碎较严重，如图4所示。根据工作面冒落区残余空间“空洞—空隙—孔隙”分布规律，将采空区按顶板岩层冒落量及矸石密积程度划分为Ⅰ，Ⅱ，Ⅲ，Ⅳ区。根据Menger海绵分形模型，計算得出冒落区矸石的统计分维值D为

2.765，则由式（2）计算可得Ⅰ，Ⅱ，Ⅲ，Ⅳ区对应孔隙率分别为：32.6%，20.27%，12.33%，8.60%。

采空区孔隙体积主要分布在冒落区，对孔隙率由冒落带在垂直方向上积分得空间滞留高度^[25]。

式中S为工作面面积，m²；h_v为滞留高度，m；V_c为冒落区残余空间体积，m³。将计算得到的冒落区孔隙率代入式（7）中可得滞留高度在冒落区不同位置分别为：2.47，1.54，0.94，0.65 m。理论计算整个试验工作面冒落区残余空间体积约为200 970 m³，占整个采空区体积比为23.3%，以冒落区走向方向单侧长度1 m内的料浆充填量用延米充填量表征，则冒落区单侧延米充填量为97.1 m³。

2.3 冒落区残余空间分布规律探查采用瞬变电磁法对试验工作面邻近的3307工作面冒落区进行残余空间地表电法探查，该工作面已回采完毕，地质采矿条件与3305试验工作面基本一致。依据实际地形和实地踏勘情况，在3307工作面设计瞬变电磁主测线21条，在主测线附近布置旁测线23条，瞬变电磁测线合计44条，瞬变电磁物理点合计1 806个，探测直径按D=H/10圈定，实际探测总面积40 484.5 m²，测点布置及探测范围如图5所示。

探测得出所有瞬变电磁测线（TMP-Z1～ TMP-Z21，TMP-P1～ TMP-P23）的视电阻率断面成果图，图6为测线TMP-Z14，TMP-Z19瞬变电磁视电阻率断面成果。

根据采空区顶板“O-X”破断规律，选取工作面南半区为代表研究区，提取所有测点250 m高程位置的视电阻率值，并根据总探测范围，绘制工作面250 m深度位置视电阻率切片成果如图7所示，3307工作面地质解译如图8所示。根据工作面地质解译平面图，对研究区内已推断异常区规模进行统计如下：推断空区（含水）异常区域共4处，总面积约为3 785.0 m²，占比约15.81%，推断空区（不含水）异常区域共5处，总面积约为1 847.2 m²，占比约7.72%。针对瞬变电磁法探测原理、探测设备以及成果解译过程中对第四系覆盖层、基岩层、煤层、含水空区、不含水空区的异常识别特征文中不再赘述。

根据上述异常区分布规律可知，冒落区残余空间多集中在边界附近，形态符合“O”形圈分布规律，实际探测到的冒落区残余空间面积占整个研究区的23.5%，探测结果为精准利用冒落区残余空间提供了指导。

3 现场充填试验

3.1 工艺流程本次试验采用在相邻巷道跨煤柱施工充填钻孔对试验工作面冒落区进行注浆充填。为简化试验流程，试验所用矸石浆体均在商混站配制后罐车输送至地面泵站，经下料斗二次搅拌后通过工业充填泵泵送至井下充填钻孔充填冒落区。试验采用KOS25100HP工业充填泵，最高泵送压力14 MPa，最大输送量150 m³/h，充填主管路为DN245×22Q345B（16 Mn），总长1 840 m，总高程差约175 m。现场充填试验工艺流程如图9所示。

3.2 充填材料选择为了提高冒落区残余空间的利用率，需将矸石破碎到一定粒径级配作为骨料，并按比例配制

成具有流动性好、扩散能力强、泌水率小的浆体。

试验矸石的破碎粒径＜3 mm，粒径占比如图10所示。

充填试验选用H-2，S-1和Y-1共3种配比浆体，所选用浆体矸石占比大、泌水率小、成本低，但因矸石参量不同，不同配比浆体的流动性和扩散能力有差异。为避免堵管，试验添加一定量的添加剂，同时配置灰浆用于充填前后润管，浆体配比参数及特点见表2。

3.3 充填钻孔设计依据矿井前期研究成果，3#煤层覆岩冒落带高度为21～26 m。基于试验工作面冒落带高度发育规律，设计高位孔如下：高位孔布置5#，6#，7#共3个，间距6 m，仰角分别为31°，45°，40°，实际施工孔深分别是18，20，40 m，终孔高度分别为15.3，20.1，31.7 m，终孔位置与煤柱水平距离为5.1，2.6，19.1 m。为充分考察浆体扩散规律，低位水平充填孔布置1#，2#，3#，12#，13#共5个，间距12，15 m，用12#，13#，4#水平孔兼作应急废料孔，充填钻孔布置如图11所示。

3.4 试验结果分析试验持续时间15 d，合计充填浆体4 949 m³，处理矸石约6 414 t。其中，高位孔充填2 421 m³，低位孔充填浆体总体积为

2 545 m³，各充填孔试验充填量见表3。

通过试验数据分析，得出了以下结论。

1）2#水平孔充填全过程孔口无压力，平均流量为120 m³/h，回撤通道有少量浆体溢出，

扩散距离为42 m，表现出无压阶段持续时间长、流量大的特点，主要以走向方向无压淌流式扩散为主。

2）5#高位孔首次充填1 063 m³，充填起始阶段压力有小微波动，二次充填浆体累计2 183 m³时孔口压力仍未出现陡增现象，充填期间13#平孔有少量水析出，浆体的实际扩散距离约为47 m，相比而言，高位孔充填覆盖冒落区垂向空间大，更利于浆体水平方向扩散以及冒落区多类型残余空间的精准高效利用。

3）12#，13#低位孔受邻近2#，5#孔充填的影响，充填量偏低，浆体沿倾向向深部扩散的动力不足，深部残余空间利用较低，结合2#，5#孔充填过程处于无压状态的特点，试验发现浆体扩散存在“通道效应”，即浆体在扩散过程中，普遍会选择阻力最小的“空洞”区域扩散。通过高位孔、低位孔协同可显著降低或解除“通道效应”和高效利用冒落区残余空间，但增大了钻孔工程量，因此对工程尺度上合理确定高、低位充填鉆孔参数提出要求。

4）通过冒落区高、低位充填孔均可以安全高效充填矸石浆体，试验得出本次试验工作面冒落区单侧的延米充填量为31.9～58.9 m³。由于试验区为整个冒落区的一部分，相当于无封闭边界约束，无压充填浆体呈椎体形态扩展，在采用单孔充填时，存在孔口正下方浆体堆积厚度较大、两侧浆体则较薄的特点，因此，本次试验冒落区残余空间未完全利用，实际延米充填量低于理论延米充填量。

4 结 论

1）指出韩城矿区煤矸石处理与利用应遵循减量化排放、源头治理、“采-充”平衡、经济高效和环境效益最佳的原则，基于冒落区残余空间分布规律揭示冒落区矸石流态化充填技术原理。

2）根据冒落岩块孔隙率与分维数关系，结合冒落区残余空间的分区特性，得出试验工作面理论预计残余空间占采空区体积比约为23.3%。同时采用瞬变电磁法对试验工作面冒落区残余空间分布规律进行了地表探查研究，推断研究区内空区（含水）异常区域4处，推断空区（不含水）异常区域5处，实际探测到的冒落区残余空间面积占整个研究区的23.5%，探测结果为精准利用冒落区残余空间提供了指导。

3）工业性试验验证冒落区高、低位充填孔均可以安全高效充填矸石浆体，总结了高、低位充填孔充填的工艺特点，得出试验工作面冒落区单侧的延米充填量为31.9～58.9 m³。研究成果为韩城矿区大规模低成本处理矸石奠定了空间容量设计依据。

参考文献（References）：

[1]曹健，黄庆享.煤炭转型期的环境保护和绿色发展[J].西安科技大学学报，2017，37（6）：779-783.CAO Jian，HUANG Qingxiang.Environmental protection and green development in coal transformation period[J].Journal of Xian University of Science and Technology，2017，37（6）：779-783.

[2]周楠，姚依南，宋卫剑，等.煤矿矸石处理技术现状与展望[J].采矿与安全工程学报，2020，37（1）：136-146.ZHOU Nan，YAO Yinan，SONG Weijian，et al.Present situation and prospect of coal gangue treatment technology[J].Journal of Mining & Safety Engineering，2020，37（1）：136-146.

[3]狄升贯.高速公路煤矸石路基路用性能研究[D].西安：长安大学，2008.DI Shengguan.Study on the road performance of coal gangue subgrade of highway[D].Xian：Changan University，2008.

[4]邱继生，杨占鲁，关虓，等.煤矸石陶粒混凝土微观孔结构特征及抗压强度[J].西安科技大学学报，2020，40（1）：110-117.QIU Jisheng，YANG Zhanlu，GUAN Xiao，et al.Microscopic pore structure and compressive strength of coal gangue ceramsite concrete[J].Xian University of Science and Technology，2020，40（1）：110-117.

[5]宋振骐，崔增娣，夏洪春，等.无煤柱矸石充填绿色安全高效开采模式及其工程理论基础研究[J].煤炭学报，2010，35（5）：705-710.SONG Zhenqi，CUI Zengdi，XIA Hongchun，et al.The fundemental theoretial and engineering research on the green safe no coal pillar mining model by mainly using coal gangue backfill[J].Journal of China Coal Society，2010，35（5）：705-710.

[6]孙希奎.矿山绿色充填开采发展现状及展望[J].煤炭科学技术，2020，48（9）：48-55.SUN Xikui.Present situation and prospect of green backfill mining in mines[J].Coal Science and Techno-logy，2020，48（9）：48-55.

[7]缪协兴，张吉雄，郭广礼.综合机械化固体充填采煤方法与技术研究[J].煤炭学报，2010，35（1）：1-6.MIAO Xiexing，ZHANG Jixiong，GUO Guangli.Study on waste-filling method and technology in fully-mechanized coal mining[J].Journal of China Coal Society，2010，35（1）：1-6.

[8]张吉雄，巨峰，李猛，等.煤矿矸石井下分选协同原位充填开采方法[J].煤炭学报，2020，45（1）：131-140.ZHANG Jixiong，JU Feng，LI Meng，et al.Method of coal gangue separation and coordinated in-situ backfill mi-ning[J].Journal of China Coal Society，2020，45（1）：131-140.

[9]周华强，侯朝炯，孙希奎，等.固体废物膏体充填不迁村采煤[J].中国矿业大学学报，2004，33（2）：154-158.ZHOU Huaqiang，HOU Chaojiong，SUN Xikui，et al.Solid waste paste filling for none-village-relocation coal mining[J].Journal of China University of Mining & Technology，2004，33（2）：154-158.

[10]赵才智.煤矿新型膏体充填材料性能及其应用研究[D].徐州：中国矿业大学，2008.ZHAO Caizhi.Study on coal mine new paste filling material properties and its application[D].Xuzhou：China University of Mining and Technology，2008.[11]刘恒凤，陈登红，陈鹏.高水复合矸石充填技术研究[J].煤炭技术，2016，35（1）：94-96.LIU Hengfeng，CHEN Denghong，CHEN Peng.Study of high-water material composite waste-filling method and technology[J].Coal Technology，2016，35（1）：94-96.

[12]陈苏社.神东矿区井下采空区水库水资源循环利用关键技术研究[D].西安：西安科技大学，2016.CHEN Sushe.Research on the key technology of water resources recycling utilization in the underground goaf reservoir in shendong mining area[D].Xian：Xian University of Science and technology，2016.

[13]钱鸣高，许家林.覆岩采动裂隙分布的“O”形圈特征研究[J].煤炭学报，1998，23（5）：466-469.QIAN Minggao，XU Jialin.Study on “O” ring characteristics of mining induced fracture distribution in overburden[J].Journal of China Coal Society，1998，23（5）：466-469.

[14]李兴尚.建筑物下条带开采冒落区注浆充填减沉技术的理论研究[D].徐州：中国矿业大学，2008.LI Xingshang.Study on mechanism of the grouting backfill in caving area with strip mining under buildings[D].Xuzhou：China University of Mining and Technology，2008.

[15]张宏贞.老采空区残留空洞空隙分布规律研究与应用[D].徐州：中国矿业大学，2013.ZHANG Hongzhen.Study and application on the distribution law of old goaf residual cavity and void[D].Xuzhou：China University of Mining and Technology，2013.

[16]胡胜勇，韩丹丹，冯国瑞，等.应力传递对垮落岩体空隙率空间分布的影響[J].采矿与安全工程学报，2020，37（4）：723-730.HU Shengyong，HAN Dandan，FENG Guorui，et al.The effect of stress transfer on the spatial distribution of void ratios in crushed rock mass[J].Journal of Mining & Safety Engineering，2020，37（4）：723-730.

[17]黄炳香，刘长友，许家林.采动覆岩破断裂隙的贯通度研究[J].中国矿业大学学报，2010，39（1）：45-49.HUANG Bingxiang，LIU Changyou，XU Jialin.Research on through-degree of overlying strata fracture fissure induced by mining[J].Journal of China University of Mining and Technology，2010，39（1）：45-49.[18]王建學，刘天泉.冒落矸石空隙注浆胶结充填减沉技术的可行性研究[J].煤矿开采，2001（1）：44-46.WANG Jianxue，LIU Tianquan.Feasibility study on subsidence reducing technology by filling with grouting cement in space of caving guange[J].Coal Mining Technology，2001（1）：44-46.

[19]张天军，王乾，羽玥，等.矸石膏体粗集料级配优化试验研究[J].西安科技大学学报，2017，37（4）：492-499.ZHANG Tianjun，WANG Qian，YU Yue，et al.Experimental investigation on coaras materal gradation optimization of gangue paste[J].Journal of Xian University of Science and Technology，2017，37（4）：492-499.

[20]古文哲，朱磊，刘治成，等.煤矿固体废弃物流态化浆体充填技术[J].煤炭科学技术，2021，49（3）：83-91.GU Wenzhe，ZHU Lei，LIU Zhicheng et al.Fluidization slurry backfilling technology of coal mine solid waste[J].Coal Science and Technology，2021，49（3）：83-91.

[21]黄庆享，赵萌烨，黄克军.浅埋煤层群开采顶板双关键层结构及支护阻力研究[J].中国矿业大学学报，2019，48（1）：71-77.HUANG Qingxiang，ZHAO Mengye，HUANG Kejun.Study of roof double key strata structure and support resistance of shallow coal seams group mining[J].Journal of China University of Mining & Technology，2019，48（1）：71-77.

[22]范立民，张晓团，向茂西，等.浅埋煤层高强度开采区地裂缝发育特征——以陕西榆神府矿区为例[J].煤炭学报，2015，40（6）：1442-1447.FAN Limin，ZHANG Xiaotuan，XIANG Maoxi，et al.Characteristics of ground fissure development in high intensity mining area of shallow seam in Yushenfu coal field[J].Journal of China Coal Society，2015，40（6）：1442-1447.

[23]卞正富，雷少刚，刘辉，等.风积沙区超大工作面开采生态环境破坏过程与恢复对策[J].采矿与安全工程学报，2016，33（2）：305-310.BIAN Zhengfu，LEI Shaogang，LIU Hui，et al.The process and countermeasures for ecological damage and restoration in coal mining area with super-size mining face at aeolian sandy site[J].Journal of Mining & Safety Engineering，2016，33（2）：305-310.

[24]李全生，张村.基于采动空间守恒的西部矿区高强度开采损伤传导模型及应用[J].采矿与安全工程学报，2021，38（1）：1-8.LI Quansheng，ZHANG Cun.Damage conduction model of high intensity mining in western mining area based on conservation of mining space and its application[J].Journal of Mining & Safety Engineering，2021，38（1）：1-8.

[25]王玉涛.采空区多孔介质空隙率与渗透特性三维空间动态分布模型[J].中国安全生产科学技术，2020，16（10）：40-46.WANG Yutao.Three-dimensional spatial dynamic distribution model on porosity and permeability characteristics of porous media in goaf[J].Journal of Safety Science and Technology，2020，16（10）：40-46.