承压水上充填开采底板破坏特征研究.pdf

187移动扫码阅读CoalMine(9):187-193REN Yanghui.Failure characteristics of floor above aquifer using backfilling mining methodJ.Safety in任仰辉.承压板破环特征研究(9):187-193.SafetyinCoalMinesSep.20232023年9 月煤砺发全No.9Vol.54第9 期第5 4 卷DOI:10.13347/ki.mkaq.2023.09.025承压水上充填开采底板破坏特征研究任仰辉（煤炭科学研究总院有限公司，北京1 0 0 0 1 3）摘要：为探索充填开采对底板破坏深度的影响及对突水灾害的控制效果，结合理论计算、数值分析和煤矿现场观测等方法对某矿A工作面承压水上充填开采底板破坏特征进行了综合研究；通过理论计算推导出了充填开采工作面底板最大破坏深度计算公式；通过数值模拟对充填开采和非充填开采时底板破坏情况进行了对比研究，通过现场实测分析了充填开采时底板岩层的实际变形破坏特征。结果表明：充填开采时，采场围岩的采动卸压范围、顶底板移近量明显减小，围岩应力集中程度降低，底板破坏深度并未随着推进距离的增加而持续纵向发展；由于充填体的应力传递和支撑作用，A工作面开采对底板的损坏影响得到减缓，破坏深度未连通下部含水层，可实现安全回采。关键词：承压水；充填开采；数值模拟；煤层底板；破坏特征中图分类号：TD823.83文献标志码：A文章编号：1 0 0 3-4 9 6 X（2 0 2 3）0 9-0 1 8 7-0 7Failure characteristics of floor above aquifer using backfilling mining methodRENYanghui(Chinese Insititue of Coal Science,Beijing 100013,China)Abstract:In order to explore the influence of backfiling mining on the damage depth of floor and the control effect of water inrushdisaster,combined with the methods of theoretical calculation,numerical analysis and field observation in coal mine,the failure char-acteristics of the floor in backflling mining above confined water in A working face of a mine are comprehensively studied.Theoret-ical calculation formula for determining maximum damage depth of the floor in backfilling workface is proposed.Through numeric-al simulation,a comparative study is carried out on the floor failure of backflling mining and non backflling mining,and the actualdeformation and failure characteristics of floor rock during flling mining are analyzed through field measurement.The results showthat surounding rock stress concentration and roof-to-floor convergence is decreasing because of the existence of the filling body;floor damage depth remain stable,no longer development in the vertical direction with advancing of the workface,mining influenceof A workface to the floor is eased for stress transfer function of the flling body,the damage depth of the floor does not reach thelimestone aquifer,and the coal seam can be safely mined.Key words:confined water;backfilling mining;numerical simulation;coal seam floor;failure characteristics煤炭作为我国主体能源，在我国能源构成中长期处于基础地位。长期以来煤炭资源的高强度开采使得地层浅部容易开采的煤层已所剩无几，不得不转而开采深部煤炭，而深部开采往往受到底板承压水的威胁。深部开采初期，由于对承压水的赋存特点及对煤炭开采的影响研究较少，缺收稿日期：2 0 2 2-0 7-2 2责任编辑：朱蕾基金项目：中国工程院咨询研究资助项目（2 0 2 2-XZ-36）作者简介：任仰辉（1 9 8 4 一），男，山东郓城人，高级工程师，硕士，主要从事煤炭发展战略与工程科技研究、煤矿安全、科技管理等方面的工作。E-mail:188Safety inCoal MinesSep.20232023年9 月No.9煤砺发全Vol.54第9 期第5 4 卷乏相应技术支撑，造成承压水上开采工作面底板突水事故时有发生。“十二五”期间发生重大水害事故1 7 起（死亡2 2 9 人），占同期全国煤矿重特大事故2 3.6%。为实现承压水上煤炭资源的安全回收，国内学者进行了大量研究：文献 1 提出了承压水上采煤的框架结构及该学科的主要研究内容和方法；文献 2 系统梳理了我国煤矿水害防治技术新进展、新方法，并对未来水害防治技术体系发展进行了展望；文献 3-6 利用相似材料模拟带压开采煤层底板破坏与裂隙发育规律；文献 7 以微震监测为基础，以微震能量密度及连通度表征底板岩层采动裂隙的导通性，为推动应用微震监测技术防治底板突水提供了新路径；文献 8-9 建立了底板突水危险性评价方法，并提出了底板奥灰水防治关键技术；文献 1 0 构建了周期来压时工作面底板力学模型，并结合钻孔注水和直流电法等方法实测了承压水上开采底板破坏深度；文献10-13采用数值模拟方法研究了底板破坏和突水通道形成过程，并建立力学模型分析了突水危险性，研究了隐伏小断层对突水的影响；文献 1 4-1 7 系统研究了承压水上膏体（似膏体）充填开采技术体系、底板采动效应等；文献 1 8 通过建立力学分析模型，分析了充填开采采空区垮落带、断裂带岩体渗透性特征，研究了隔水关键层对水害防治的作用；文献 1 9-2 0 研究了固体充填材料影响岩层控制效果的机理，建立了固体充填开采充填材料选择设计方法，设计了煤矿井下煤研分选与充填开采一体化系统。以上研究成果为受承压水威胁采场的底板控制提供了重要参考，但深部开采面临的高应力、高承压水等特殊条件，仍然对深部煤炭资源的安全高效经济开采带来严峻挑战。同时，充填开采在减少煤层底板变形方面效果明显，且随着技术进步已得到越来越多的应用。因此，进一步研究深部承压水上煤层采用充填开采时的底板破坏效应具有重要现实意义和应用价值。为此，采用理论计算、数值分析和现场观测相结合的方法对工作面承压水上充填开采底板破坏特征进行了综合研究，研究结果有助于摸清充填开采对底板水害的控制作用，为保障承压水上的安全开采提供参考依据。1工作面概况某矿A工作面地面标高+7 4 +7 5 m，工作面标高-5 3 0 -6 0 0 m，主采七号煤，煤层产状大致为走向北西，倾向北东，倾角5 1 1，平均7，七号煤层赋存较为稳定，煤厚平均1.5 m。工作面直接顶为泥岩，黏结性较好，可随采随垮落，基本顶为细砂岩，厚6 m，强度中等，属易管理顶板。7号煤层下距底板四灰含水层2 2 m左右，四灰含水层水量不大。下距奥陶系石灰岩含水层4 0 m左右，奥灰含水层水量丰富。在四灰和奥灰含水层之间，赋存有厚度约1 7 m的以泥岩为主的岩层，隔水性较好。2充填开采底板破坏深度分析充填开采是指煤层采出后，利用回采和掘进石、膏体充填料、高水材料或高浓度胶结材料等充填体充填采空区，使充填体和周边煤柱共同支撑上覆基岩载荷。由于充填工艺、充填材料和工作面管理等问题，充填体与顶板岩层之间一般存在一定的空顶高度。空顶高度的存在使充填体不能及时发挥承载能力，因此工作面推过之后的一定时期内，上方基岩载荷仍由周边煤柱和工作面前方煤壁支撑，导致支承压力集中现象的产生。对于采空区底板中的任意一点M，其受力一般可认为由4 个部分构成，包括M之上底板岩体的自重、充填体之上的顶板载荷、充填体的重力作用和工作面煤壁前面一定距离顶板岩层的支承压力作用。根据圣维南原理，作用在物体的一小部分边界上的面力，变换为分布不同但作用等效的面力，那么近处的应力分布将有显著的改变，但远处所受的影响可以忽略不计 2 1 ，因此对M点应力状态造成影响的非线性分布的支承压力可简化为分段载荷，充填开采底板应力简化模型如图1。图1 中，q为煤壁前方支承压力的等效应力，其值等于pgH（n+1)/2；p 为煤岩体密度；H为煤层埋藏深度，n为支承压力集中系数，一般取2 3。根据相关力学原理可知q的作用宽度L=2xa。由土力学知识可知，M点应力状态是支承压力同顶板载荷、充填体重力、底板载荷共同决定的，对各因素对M点应力状态的影响进行叠加可得底板岩层中任意一点M的主应力为：(+sina)+pegHe+pgz元(1)9(-sin)+pegH+pgz元式中：为最大主应力；为最小主应力；189Safety in Coal MinesSep.20232023年9 月Vol.54No.9煤砺发全第9 期第5 4 卷npgH充填体PgH充填体pgHXX2NNMML（a）支撑压力示意图（b）支撑压力等效应力示意图图1充填开采底板应力简化模型Fig.1Simplified model of floor stress in back filling miningp为充填体密度；g为重力加速度；H.为充填体稳定后高度；为最大主应力与x轴的夹角；z为M点到煤层底版的深度。代人Mohr-Coulomb破坏准则，即1-K3=R，整理后可得：9K+1R。PZ=sin-H(2)元p8K-1pg(K-1)P式中：R。为岩体单轴抗压强度；K=1+sinpo1-singoo为底板岩体内摩擦角。对式(2)进行求导，即可得到底板破坏深度达到极大值时的角：dzK+1cos-1=0(3)da元pgK-1由式(3)可得当角度满足式(4)时，底板破坏深度达到极大值。K-1=cos51(4)K+1将式(4)及q=pgH(n+1)/2代人式(2)可求得充填开采工作面底板岩层最大破坏深度H为：n+12VKK-1)RPH。HmHCOS2元K-1K+1)pg(K-1)P(5)由式（5)可以看出，采用充填开采工艺时，煤层底板最大破坏深度与煤层埋深、煤壁前方支承压力成正比，与底板岩体强度、充填体密度及充填高度负相关，将A工作面及煤岩体物理力学参数代人式（5）可得该工作面底板的最大破坏深度为20m，由于在承压水导升带与采动影响破坏带之间存在厚1 7 m的隔水关键层，由理论计算可以看出，采用充填采煤法控制底板变形时，A工作面没有突水危险，可以实现七号煤层的安全回采。3承压水上充填开采的数值模拟根据A工作面顶底板岩层综合柱状图建立的数值模型如图2，模型高1 5 7 m，走向长4 0 0 m，模型四周和底部固定位移约束，顶部施加8 MPa的应力边界条件，底板奥灰岩中施加3.5 MPa的孔隙水压，煤岩体破坏准则选取库伦-摩尔准则，采用充填采煤工艺，充填材料选取应变硬化模型，并假设充填工作与工作面开采并行。煤岩及充填体物理力学参数见表1。泥岩大青灰岩粉砂岩ULST开切眼粉砂岩泥岩互层本溪灰岩铝王泥岩奥陶系石灰岩煤400m细砂岩图2数值模型Fig.2Numerical calculationmodel表1 煤岩层物理力学参数表Table1Physical and mechanical parameters ofcoal and rock参数弹性模黏聚力/内摩擦 抗拉强密度/泊松比岩性量/GPaMPa角/)度/MPa(kg m-3)粉砂岩1.050.854.90401.202.660煤0.480.151.00250.031400石灰岩4.492.824.00390.952.650细砂岩0.980.925.00421.102.650泥岩1.780.511.90300.892610灰岩7.853.854.00381.792.760充填体0.030.022.55231.401870190SafetyinCoal MinesSep.20232023年9 月No.9Vol.54煤矿发全第9 期第5 4 卷3.1推进不同距离顶底板有效应力分布充填开采情况下，工作面推进至不同距离时，工作面顶底板有效应力分布如图3。-12.287-12MPa工作面20.9417(a)推进5 0 m-15.173-14MPa工作面21.6027(b)推进 1 5 0 m图3推进至不同距离时围岩有效应力分布Fig.3Effective stress distribution of surrounding rock atdifferent distances由于充填体的承载作用，煤壁前方支承压力集中系数很小，工作面推进5 0 m时约为1.2 2，推进至1 5 0 m时应力集中系数约为1.5 1，增幅较小。当推进距离不大时，由于顶板下沉量较小，顶板同充填体之间的接触耦合作用不强，充填体发挥的承载能力有限，如图3(a），此时同垮落法管理顶板类似，工作面顶底板出现较大范围的卸压区；而当工作面推进至1 5 0 m时，如图3(b），由于顶板已多次垮落，覆岩充分下沉，充填体的承载能力得以充分发挥，此时采空区中部垂直应力已恢复至原岩应力水平。另外，由于承压水的存在，距离七号煤3 5 m处的奥灰岩承受的有效应力水平很低，特别是开挖卸荷后，石灰岩中出现较大范围的低应力区，应力状态的改变致使奥灰岩容易发生破坏，可能使岩体中的裂隙不断发育并相互贯通，并形成承压水导升带。3.2推进不同距离时顶底板垂直位移分布工作面推进至不同距离时，顶底板位移场分布如图4。工作面推出开切眼5 0 m时，顶板最大下沉量320mm，顶板同充填体没有充分接触，因此充填体对底板的控制效果较差，此时，最大底鼓量为180mm；工作面推进至1 5 0 m时，顶板最大下沉量达到7 2 0 mm，顶板同充填体之间的空顶高度已不存在，基岩同充填体密实接触，充填体被压实，320.53300mm工作面150182.37(a)推进5 0 m-719.85-700mm工作面100132.68(b)推进 1 5 0 m图4#推进至不同距离时围岩位移场Fig.44Displacement field of surrounding rock atdifferentdistances对底板的控制效果转好，此时底板最大鼓起高度为1 2 0 mm，说明在充填体的作用下，已采空间顶底板位移受到约束，卸压范围减小。3.3推进至不同距离底板破坏深度工作面推出开切眼不同距离时，采场围岩塑性破坏区分布如图5。Noneshear-n shear-pshear-p工作面承压水头隔水关键层(a)推进 5 0 mNoneshear-n shear-pshear-p工作面承压水头隔水关键层(b)推进1 5 0 m图5#推进至不同距离时围岩塑性区分布Fig.5Distribution of surrounding rock plastic zone atdifferentdistances由于充填作用，底板破坏深度较小，且随着工作面的推进，底板破坏深度并没有沿垂直方向发展，推进至5 0、1 5 0 m时底板破坏深度大致相同，均为1 8 m左右。在承压水导升带与采动影响破坏带之间存在厚1 7 m的隔水关键层，说明在使191SafetyinCoalMines2023年9 月Sep.2023煤砺发全No.9Vol.54第9 期第5 4 卷用充填采煤工艺的条件下，A工作面没有突水危险3.4非充填开采情况下数值模拟结果对比为进一步对比充填开采的效果，进行了非充填开采情况下的模拟计算。除未进行工作面的充填开采外，其他条件与充填开采计算模型保持一致。非充填开采情况下工作面推进到不同距离时的围岩塑性区分布如图6。承压水头(a)推进 5 0 m承压水头(b)推进1 5 0 m图6 非充填开采工作面推进至不同距离时围岩塑性区分布Fig.6Distribution of surrounding rock plastic zonewhenthe non backfilling face advances to different distances工作面推进5 0 m时，底板破坏深度在2 0 m左右，超出采用充填开采时底板破坏深度约2 m，此时采动影响尚未波及底板奥灰含水层；当工作面推进1 5 0 m后，随着底板岩层塑性区的扩展，受采动影响的底板岩层破坏范围已接近4 0 m，超出采用充填开采时底板破坏深度约1 2 m。数值模拟结果表明，非充填开采情况下，工作面推进到1 5 0 m时，与充填开采相比，底板破坏深度明显增加，已与奥灰含水层导通，极有可能诱发奥灰突水事故。也就是说，采用充填开采，可以有效减小煤层开采对底板的扰动，减小底板破坏深度，有效防治底板突水事故。4底板采动效应的现场实测A工作面采用充填采煤工艺后，为得到在采动和承压水的双重影响下，底板岩层的实际变形破坏特征，在距离A工作面底板约2 9 m处的八号煤泄水巷道中超前工作面安装动态仪，共设置1 2个测站，通过记录动态仪数据，观测泄水巷道的顶底板移近量及巷道的破坏情况。测点2 至测点12等1 1 个测点的顶底板移近量变化曲线（测点1由于出现故障，未获得有效数据），测点顶底板移近量随工作面推进距离的变化曲线如图7。140120100807#608#9#4010#2011#12#020406080100120140160180工作面推采距离/m图7测点顶底板移近量随工作面推进距离的变化曲线Fig.7Change curves of roof and floor movement withadvancing distance of working face由图7 可以看出：在A工作面推进距离小于90m之前，各测点巷道的顶底板移近量都很小，基本不受其开采活动的影响；当A工作面推进距离达到9 0 m后，由于6 测点位于工作面正下方，且存在断层，其顶底板移近量数值变化突然增大，受采动影响较为剧烈；在A工作面推进距离达到100m时，位于A工作面回风巷正下方的4#测点数值也发生显著变化，但顶底板移近量为4 5 mm后变化不大。以上代表性测点的观测数据表明：在A工作面开采初期，其开采活动对八号煤泄水巷道基本没有影响，而当A工作面推进9 0 m后，个别测点由于A工作面采空区正下方或处于A工作面回风巷正下方，受到较大影响；但是在A工作面推进过程中，泄水巷道并没有出现很大变形，一直可以利用，说明在充填开采条件下，七号煤的采动影响并不会造成底板突水，可以实现安全回采。5结语1)结合岩石力学、土力学及弹性力学基本原理，建立了承压水上充填开采底板破坏深度分析的力学模型，并推导出了底板最大破坏深度的求解公式，由该公式可计算得出A工作面充填开采时底板最大破坏深度为2 0 m，发生在工作面后方，工作面没有突水危险，可实现安全回采。192Safetyin Coal MinesSep.20232023年9 月煤砺发全No.9Vol.54第9 期第5 4 卷2)通过数值模拟得到了充填体对采场围岩应力场、位移场及塑形破坏区域的控制作用。在充填体作用下，采场围岩的采动卸压范围、顶底板移近量明显减小，围岩应力集中程度降低，底板破坏深度并未随着推进距离的增加而持续纵向发展，底板破坏深度约1 8 m，同理论计算结果基本吻合。3)A工作面推进过程中，对底板中八号煤层泄水巷道的变形破坏特征进行了现场实测，在七号煤采动影响下，泄水巷道没有因发生较大变形而废弃，只有局部出现片帮冒顶现象，说明充填开采条件下，A工作面的底板采动效应得到很大程度的缓和。参考文献(References):1王作宇，刘鸿泉.承压水上采煤学科理论与实践 J.煤炭学报，1 9 9 4,1 9(1)：4 0-4 8.WANG Zuoyu,LIU Hongquan.Theory and practice ofcoal mining discipline on confined waterJ.Journal ofChina Coal Society,1994,19(1):40-48.2靳德武.我国煤矿水害防治技术新进展及其方法论思考 J.煤炭科学技术，2 0 1 7,4 5(5)：1 4 1-1 4 7.JIN Dewu.New development of water disaster preven-tion and control technology in China coal mine and con-sideration on methodologyJ.Coal Science and Tech-nology,2017,45(5):141-147.3张培森，安羽枫，武守鑫，等.带压开采煤层底板破坏突水规律 J.煤矿安全，2 0 1 9,5 0(2)：2 5-2 9.ZHANG Peisen,AN Yufeng,WU Shouxin,et al.Waterinrush laws of floor failure in coal seam with pressureminingJ.Safty in Coal Mines,2019,50(2):25-29.4王厚柱，鞠远江，秦坤坤，等.深部近距离煤层开采底板破坏规律实测对比研究 J.采矿与安全工程学报，2020,37(3):553-561.WANG Houzhu,JU Yuanjiang,QIN Kunkun et al.Acomparative study on floor failure law in deep and short-distance coal seam miningJ.Journal of Mining&Safety Engineering,2020,37(3):553-561.5李杨杨，张士川，孙熙震，等.煤层采动底板突水演变过程可视化试验平台研制与试验研究 J.煤炭学报，2021,46(11):3515-1524.LI Yangyang,ZHANG Shichuan,SUN Xizhen,et al.Development and experimental study on visualizationtest platform for water inrush evolution process of coalseam mining floorJ.Journal of China Coal Society,2021,46(11):3515-1524.6宋文成，梁正召，刘伟韬，等.采场底板破坏特征及稳定性理论分析与试验研究 J.岩土工程学报，2 0 1 9,38(1):2208-2218.SONG Wencheng,LIANG Zhengzhao,LIU Weitao,etal.Theoretical analysis and experimental investigationon failure characteristics and stability of stope floors J.Chinese Journal of Rock Mechanics and Engineering,2019,38(1):2208-2218.7靳德武，段建华，李连崇，等.基于微震的底板采动裂隙扩展及导水通道识别技术研究 J.工程地质学报，2021,29(4):962-971.JIN Dewu,DUAN Jianhua,LI Lianchong,et al.Micro-seismicity based research for mining induced fracturepropagation and water pathway identification techno-logy of floorJ.Journal of Engineering Geology,2021,29(4):962-971.8朱开鹏.基于双指标的底板突水危险性评价方法 J.煤矿安全，2 0 2 0,5 1(8)：2 3 2-2 3 8.ZHU Kaipeng.Evaluation method for floor water-in-rush risk based on double-indexJ.Safety in CoalMines,2020,51(8):232-238.9庞贵艮.保德煤矿底板奥灰水害防治关键技术 J.煤矿安全,2 0 2 0,5 1(1)：7 5-7 9.PANG Guigen.Key technology of prevention and con-trol of water disaster of ordovician limestone in BaodeCoal Mine FloorJ.Safety in Coal Mines,2020,51(1):75-79.10 商薛卫峰，王苏健，黄克军，等。承压水矿区切顶沿空留巷底板破坏理论与应用 J.煤炭学报，2 0 2 0,4 5(S2)：581-588.XUE Weifeng,WANG Sujian,HUANG Kejun,et al.Theoretical analysis and field measurement of floor fail-ure in gob side entry of cutting roof in confined watermine areaJ.Journal of China Coal Society,2020,45(S2):581-588.【1 1 宋文成，梁正召.承压水上开采倾斜底板破坏特征与突水危险性分析 J.岩土力学，2 0 2 0,4 1(2)：6 2 4-6 3 4.SONG Wencheng,LIANG Zhengzhao.Investigation onfailure characteristics and water inrush risk of inclinedfloor mining above confined aquiferJ.Rock and SoilMechanics,2020,41(2):624-634.【1 2 冯宇，林志斌，王文贞.采动影响下底板隐伏小断层突水灾变规律研究 J.煤矿安全，2 0 2 2,5 3（8：2 0 0-2 0 9.FENG Yu,LIN Zhibin,WANG Wenzhen.Research onwater inrush disaster law from coal mine floor with con-cealed small fault under miningJ.Safety in CoalMines,2022,53(8):200-209.1 3 王进尚，姚多喜，黄浩.煤矿隐伏断层递进导升突水的临界判据及物理模拟研究 J.煤炭学报，2 0 1 8,4 3(7)：193Safetyin Coal MinesSep.20232023年9 月No.9Vol.54煤砺安全第9 期第5 4 卷2014-2020.WANG Jinshang,YAO Duoxi,HUANG Hao.Criticalcriterion and physical simulation research on progress-ive ascending water inrush in hidden faults of coalmines J.Journal of China Coal Society,2018,43(7):2014-2020.1 4 廖志恒.承压水上上膏体充填开采底板采动破坏特征J.煤矿安全,2 0 1 8,4 9(4):1 8 5-1 8 8.LIAO Zhiheng.Mining-induced damage characteristicsof floor by paste filling mining on confined water J.Safety in Coal Mines,2018,49(4):185-188.1 5 邱梅，施龙青，于小鸽，等.似膏体充填开采防止高承压底板水突出研究 J.采矿与安全工程学报，2 0 1 6,33(3):409-422.QIU Mei,SHI Longqing,YU Xiaoge,et al.On prevent-ing the high pressure water-inrush from floor with paste-like backfill miningJ.Journal of Mining&Safety En-gineering,2016,33(3):409-422.1 6】轩召军，常庆粮，膏体充填开采底板破坏机理及分区治理技术 J.煤矿安全，2 0 2 2,5 3(8)：2 1 0-2 1 7.XUAN Zhaojun,CHANG Qingliang.Study on floorfailure mechanism of paste backfill mining and divisionarea control technology J.Safety in Coal Mines,2022,53(8):210-217.1 7 孙希奎，常庆粮，施现院，等.膏体充填开采底板破坏规律充填效应及实测研究 J.煤矿安全，2 0 1 8,4 9(1)：76-84.SUN Xikui,CHANG Qingliang,SHI Xianyuan,et al.Study on filling effect and measurement of floor failurein paste backfill mining J.Safety in Coal Mines,2018,49(1):76-84.1 8 孙强，刘恒凤，张吉雄，等.充填开采隔水关键层渗流稳定性力学模型及分析 J.采矿与安全工程学报，2022,39(2):273-281.SUN Qiang,LIU Hengfeng,ZHANG Jixiong,et al.Mechanical model and analysis of seepage stability inkey aquiclude strata during backfill miningJ.Journalof Mining&Safety Engineering,2022,39(2):273-281.1 9 王云搏，张强，孟国豪，等.基于固体充填开采的充填材料选择设计方法.采矿与岩层控制工程学报，2022,4(3):033035.WANG Yunbo,ZHANG Qiang,MENG Guohao,et al.Methods of selecting filling materials for solid fillingmining J.Journal of Mining and Strata Control Engin-eering,2022,4(3):033035.【2 0】赵振宇，张开加，董宇.鑫岩煤矿井下煤研分选与充填开采一体化系统设计 J.煤炭工程，2 0 2 2,5 4(1)：11-17.ZHAO Zhenyu,ZHANG Kaijia,DONG Yu.Integratedsystem design of coal and gangue separation and back-fill mining in Xinyan Coal MineJ.Coal Engineering,2022,54(1):11-17.【2 1 任仰辉.复杂条件下底板采动效应特征及水害综合防治技术研究 D.青岛：山东科技大学，2 0 0 9:1 9-2 0.