移动阅读李树刚1，魏宗勇1，林海飞1，赵鹏翔1，肖 鹏1，郝昱宇2
(1.西安科技大学 安全科学与工程学院，陕西 西安 710054; 2.西安科技大学 工程训练中心，陕西 西安 710054)
摘 要:为进一步解决煤与瓦斯共采模型实验研究手段不足的问题，自主研制了一套煤与瓦斯共采三维大尺度物理模拟实验系统。该系统采用模块化设计，高度集成机、电、液、气于一体，主要由大尺度箱体(3.0 m×2.5 m×1.8 m)与基座、自动液压开采、柔性加载、自动通风、瓦斯抽采、瓦斯注入以及综合数据采集与控制等7个子单元构成。按几何相似比1∶100计，加载单元可模拟最大采深2 105 m，开采单元可模拟采高0～12 m以及推进距离200 m;通风单元可模拟U型、U+L型、Y型等多种通风方式以及实现不同风量通风;抽采单元可模拟高位巷、高位钻孔、地面抽采等多种立体化抽采方式;瓦斯注入单元采用独立注入方式，实现不同瓦斯涌出量、不同位置的瓦斯涌出;综合数据采集与控制单元实现覆岩裂隙、矿山压力、瓦斯运移、瓦斯抽采等表征参数的采集以及对整个实验系统进行自动控制。该实验系统可进行工作面煤层开采、通风、瓦斯涌出与抽采等功能的模拟，实现煤层开采过程中覆岩裂隙演化、矿山压力分布、卸压瓦斯运移、瓦斯抽采等科学问题的一体同步研究。运用该系统对山西某矿302工作面开采过程进行模拟实验，得到了该矿条件下基本顶初次来压步距45 m，周期来压步距20 m，覆岩破坏在空间上呈椭圆抛物形态等覆岩破断与裂隙演化规律;工作面推进过程中应力峰值不断前移，应力集中系数2.11～2.63，超前工作面距离6～11 m等动态应力变化规律;在卸压瓦斯储集与分布规律方面，得到采空区后部76～120 m瓦斯浓度增加较快，120 m之后趋于稳定，采空区上部5～60 m裂隙带中瓦斯浓度逐渐增加，裂隙带最上层瓦斯浓度达到65%～68%。实验结果表明，该系统能够较好进行工作面煤与瓦斯共采全过程的模型实验研究。
关键词:大尺度;煤与瓦斯共采;裂隙演化;瓦斯运移;三维模拟实验
中图分类号:TD712
文献标志码:A
文章编号:0253-9993(2019)01-0236-10
移动阅读李树刚，魏宗勇，林海飞，等.煤与瓦斯共采三维大尺度物理模拟实验系统的研制与应用[J].煤炭学报,2019,44(1):236-245.doi:10.13225/j.cnki.jccs.2018.1635
LI Shugang,WEI Zongyong,LIN Haifei,et al.Research and development of 3D large-scale physical simulation experimental system for coal and gas co-extraction and its application[J].Journal of China Coal Society,2019,44(1):236-245.doi:10.13225/j.cnki.jccs.2018.1635
收稿日期:2018
1212
修回日期:20180107
责任编辑:常明然
基金项目:国家自然科学基金科学仪器专项资助项目(51327007)；国家自然科学基金重点资助项目(51734007)；国家自然科学基金青年基金资助项目(51704228)
作者简介:李树刚(1963—)，男，甘肃会宁人，教授，博士生导师。Tel:029-85583143，E-mail:list63@163.com
通讯作者:魏宗勇(1986—)，男，陕西汉中人，工程师，博士研究生。E-mail:342583243@qq.com
LI Shugang1,WEI Zongyong1,LIN Haifei1,ZHAO Pengxiang1,XIAO Peng1,HAO Yuyu2
(1.College of Safety Science and Engineering,Xi’an University of Science and Technology,Xi’an 710054,China; 2.Engineering Training Centre,Xi’an University of Science and Technology,Xi’an 710054,China)
Abstract:In order to overcome the shortcomings of experimental research methods for coal and gas co-extraction model,a set of 3D large-scale physical simulation experimental system for coal and gas co-extraction was developed.The system adopts modular design and highly integrates the machine,electricity,liquid and gas functions in one.It is mainly composed of seven sub-units including large-scale box body (3.0 m×2.5 m×1.8 m) and base,automatic hydraulic mining,flexible loading,automatic ventilation,gas extraction,gas injection,and integrated data acquisition and control.Following the geometric similarity ratio of 1∶100,the loading unit can simulate a maximum mining depth of 2105 m,the mining unit can simulate a mining height of 0-12 m and a face retreat distance of 200 m.The ventilation unit can simulate various ventilation modes such as U-type,U+L-type,Y-type with different flow quantities.The gas extraction unit can simulate gas drainage using different kinds of methods and a combination of them,such as a high-level gas drainage roadway,cross measure borehole and surface borehole.The gas injection unit adopts independent injection mode to realize gas emission at different positions and different gas emission rates.Comprehensive data acquisition and control unit is able to obtain the characterization parameters such as overburden fissures,strata pressure,gas migration,gas extraction and automatic control of the whole experimental system.Therefore,the experimental system can simulate coal seam mining process under different ventilation,gas emission and extraction conditions,making it possible to simultaneously study the overburden fracture evolution,mine pressure distribution,pressure relieved gas migration,gas extraction and other scientific issues during coal mining operation.The system was used to simulate the mining process at the 302 working face of a mine in Shanxi province,China.The laws of overburden breakage and fracture evolution were obtained,such as,the first weighting of the main roof was 45 m,the periodical weighting interval is 20 m,and the overburden failure shows a 3D elliptical parabolic pattern.The dynamic variation of stress was obtained,indicating the peak stress moves forward continuously with the stress concentration coefficient variation between 2.11 and 2.63,and the peak stress occurrence at about 6-11 m ahead of the working face.In terms of the storage and distribution of pressure relieved gas,it was found that the gas concentration increased rapidly at 76-120 m behind face in the goaf and became stable beyond 120 m.The gas concentration increased gradually from the lower section to the higher section of the fractured zone which was about 5-60 m above the extraction level,reaching about 65%-68% at the top of fractured zone.The experimental results show that the system can well carry out the model experimental study of the whole process of coal and gas co-extraction in working face.
Key words:large scale;coal and gas co-extraction;fracture evolution;gas migration;3D simulation exper-iment
煤炭是我国的主体能源，根据国家数据统计在一次能源生产和消费中比重分别占69.6%和62%[1-2]。随着煤层开采强度的增大、开采深度的增加，瓦斯含量及瓦斯涌出量逐渐增大，瓦斯灾害也日趋严重[3-4]。由于我国大部分矿区煤层透气性较差，预抽效果不甚理想，采后卸压抽采为实现煤与瓦斯共采的一种有效手段，其关键科学问题是采动覆岩裂隙动态演化规律及卸压瓦斯运移规律。为研究该问题，物理模拟实验、数值模拟以及现场工业性实验等手段得以应用，相较于后两者，物理模拟实验具有可靠性高、直观性强与重复性好等特点，已成为研究煤与瓦斯共采科学问题的主要方法[5]。
早期，物理模拟实验主要以平面应力模型或平面应变模型为主，无法真正体现三维空间状态下的规律[6-7]。近年来，为了解决二维模型存在的问题，研制了一批仿真程度更高的三维物理模拟实验平台。李仲奎等[8-9]发明了三维多主应力面加载实验系统，成功完成了复杂洞室群模型实验。胡耀青等[10-11]建立的大型三维固流耦合模拟实验系统，研究了覆岩位移和应力变化规律，为矿井煤层开采覆岩活动及地表沉降的研宄提供了基础指导。伍永平等[12-14]依托自主研发的多功能立体支架实验装置，探究了急斜煤层开釆煤岩多场耦合失稳机理、围岩变形破坏及应力分布特征。李树刚等[15-17]主持研制了煤与瓦斯共采小型三维物理模拟实验台，探究了煤层开采过程中采动应力分布和覆岩破坏规律。尹光志等[18-20]应用三维物理模拟的方法对煤与瓦斯突出进行了深入的研究。已有研究取得了丰富成果，初步实现了实际工程原型的模型试验，但以往研究主要集中于地应力对围岩巷道影响，或侧重于高应力、高瓦斯压力下煤与瓦斯突出规律，对于煤与瓦斯共采物理模拟涉及较少;另外，覆岩裂隙演化、矿山压力以及卸压瓦斯储运与抽采一体同步研究的大型三维物理模拟实验平台功能实现仍需完善。
为进一步解决煤与瓦斯共采模型实验研究的不足，笔者在充分借鉴前人经验的基础上，自主研制了一套煤与瓦斯共采三维大尺度物理模拟实验系统。该系统可相对真实的再现工作面煤层开采、通风、瓦斯涌出及抽采环境，将覆岩裂隙演化、矿山压力分布以及卸压瓦斯储运与抽采等进行同时同步研究，从而为实现煤与瓦斯共采提供一定的理论依据。
煤与瓦斯共采三维大尺度物理模拟实验系统能够进行煤层开采条件下覆岩裂隙演化、矿山压力分布、卸压瓦斯储运及抽采的同步化物理模拟实验研究。实验平台采用模块化设计方法，主要由7个单元构成，分别是大尺度箱体与基座、自动液压开采单元、柔性加载单元、自动通风单元、瓦斯注入单元、瓦斯抽采单元、综合数据采集与控制单元，可实现应力场、瓦斯场、裂隙场等多场表征参数的动态采集及控制(图1)。
图1 实验系统总设计构架
Fig.1 General design framework of experimental system
1.2.1 大尺度实验箱体与基座
根据目前煤矿综合机械化采煤工作面的面长与推进距离，以及模拟实验的要求，确定实验模型尺寸为3.0 m×2.5 m×1.8 m，如图2所示。箱体采用高强度加强筋焊接钢板拼接而成，模型正面与背面分别设计螺纹孔，用于测试覆岩裂隙中瓦斯压力以及瓦斯浓度;模型采空区侧可接入多路传感器线，用于测试煤层顶底板应力、瓦斯浓度;大巷侧设计有通风孔、瓦斯抽采孔，可模拟多种通风方式以及模拟高位巷、高位钻孔抽采。大尺度实验箱体的基座采用钢筋混凝土结构，承载整个实验系统。
图2 大尺度实验箱体
Fig.2 Large-scale experimental box
1.2.2 自动液压开采单元
自动液压开采单元用于模拟封闭状态下煤层开采全过程。该单元主要由特制条钢、高承载液压缸、电磁换向阀、位移传感器、PLC电气控制组成闭环控制系统，实现封闭空间中煤层开采及煤层开采速度控制的难题，也可模拟不同煤层厚度。实现原理是每2个液压缸支撑1个特制条钢，构成刚性联结，液压缸缸体以法兰方式联结于箱体底部;实验开始前将所有液压缸升至预定高度形成未开采煤层，然后铺设模型，煤层的开采过程就是逐个降下特制条钢。以1∶100几何相似比，可实现工作面长200 m、采高0～12 m、推进距离200 m的煤层开采过程真实模拟(图3)。
图3 自动液压开采单元
Fig.3 Automatic hydraulic mining unit
1.2.3 柔性加载单元
随模型实验煤层开采的不断推进，物理模型顶面会产生不均匀的弯曲下沉，为实现煤层开采后模型顶部的均布加载，研制了液压柔性加载单元。柔性加载单元由大型液压站、高强度刚性龙门、特制压头、弹簧组以及特制橡胶构成，采用伺服液压控制模式精确控制加载力(图4)。加载单元整体安装在直线导轨上，并由减速机-齿轮-链条拖动。当铺设模型时将加载单元移动到箱体外侧，利于模型装填;加载时，将加载单元移动到模型上方，采用96组高强螺栓与基座固定。柔性加载单元最大加载力为3×106 N，按几何相似比1∶100以及新型相似材料容重[16]计，最大模拟采深2 105 m。
图4 柔性加载单元
Fig.4 Flexible loading unit
1.2.4 自动通风单元
自动通风单元主要模拟矿井通风，实现对采煤工作面的精准配风，可模拟U型、U+L型、Y型等多种通风方式。该单元主要由通风巷道、常闭电磁阀、HG-370型风机、各种传感器(流量、浓度、负压传感器)、通风管路以及控制系统等组成，风量采用调节器控制为比例、积分、微分的PID控制方式，实现通风量的精确控制。
U型通风设计采用自成巷道原理，实验模型装填前，将通风特制条钢液压缸升至进风、回风巷道的比例高度;模型装填完成、晾干以后，降下进风巷道、回风巷道条钢，形成U型通风系统，原理如图5所示。
图5 U型通风设计原理
Fig.5 Design principle of U-type ventilation
U+L型通风设计原理:以U型通风巷道为基础，采用有机玻璃制成L型巷道模型，巷道断面尺寸为50 mm×40 mm(宽×高)。U型与L型之间联络巷用电磁阀控制开启、封闭，联络巷之间距离设计为35 cm，共设计6条联络巷(图6)。Y型通风主要采用两进一回通风方式，沿空埋设管路。
图6 U+L型通风设计原理
Fig.6 Design principle of U+L-type ventilation
1.2.5 瓦斯注入单元
瓦斯注入单元用于模拟采区瓦斯涌出规律，可实现采空区、工作面等不同位置的瓦斯涌出，也可实现不同瓦斯涌出量的模拟。该单元为实验模型提供瓦斯分源涌出，主要由瓦斯钢瓶、瓦斯注入控制柜、管路等组成，其中瓦斯注入控制柜内含减压阀、储气罐、气体压力表与出口压力控制阀、40路气体质量流量控制器、控制板卡、显示器等，高度集成。瓦斯气体经质量流量控制器控制流量，由实验箱体底部注入特制开采条钢，开采条钢上密集布置筛孔，以利于均匀散溢瓦斯(图7)。
图7 瓦斯注入单元
Fig.7 Gas injection unit
1.2.6 瓦斯抽采单元
瓦斯抽采单元主要模拟高位巷、高位钻孔、地面钻孔以及上隅角瓦斯插管或埋管抽采，研究最佳抽采位置及不同抽采方式的瓦斯抽采效率。该系统主要由抽采管路、调节阀、调速抽采泵以及瓦斯浓度传感器、流量传感器、负压传感器以及采集与控制部分构成。瓦斯抽采流量、负压、浓度数据可实时采集及控制(图8)。
图8 瓦斯抽采单元
Fig.8 Gas extraction unit
1.2.7 综合数据采集与控制单元
采动覆岩破断及裂隙演化的表征参数采集主要采用高频微震监测系统采集震动信号与高清窥视系统观测裂隙发育。
高频微震监测系统主要由高频微震传感器、前置放大器、信号采集系统和主机分析系统组成(图9)。系统能够实现信号特征参数和波形的实时同步采集和分析，并可对信号进行三维空间定位。
高清窥视系统主要由钻孔装置、清孔装置和高清工业内窥镜组成。高清窥视系统可对相似材料模型进行任意高度、角度和方位的钻孔，并可对模型内部裂隙的发育情况进行清晰成像和记录，得到可靠的裂隙参数数据。
图9 高频微震监测系统构成
Fig.9 Composition of high frequency microseismic monitoring system
矿山压力分布、卸压瓦斯运移、瓦斯抽采以及通风表征数据采集，主要使用压力传感器、瓦斯浓度传感器、瓦斯压力传感器、流量传感器等多种传感器;采用以SQL数据库为基础，TCP/IP网络通讯为桥梁，C++build为方法的设计模式(图10)，编制了综合实验数据采集软件。其中SQL数据库作为系统的数据存储单元，具有专业的数据操作命令集，是一种功能齐全的数据库语言。通过该数据库可快速实现数据查询、存储、显示功能。TCP/IP网络通讯构架可以实现大型数据的快速传输，并为系统上网提供了渠道。通过该设计模式完成了煤与瓦斯共采三维实验监测软件的实时显示、数据处理、在线存储、网络互联的实验需求。
图10 煤与瓦斯共采数据采集构架
Fig.10 Data acquisition framework for coal and gas co-extraction
煤与瓦斯共采三维大尺度物理模拟实验系统由各功能单元有机结合(图11)，高度集成、协调运行，能够达到工作面的煤层开采、通风、瓦斯涌出及瓦斯抽采功能，依托该实验系统根据工作面实际条件设计实验，实现将覆岩移动破断、上覆岩层裂隙演化、矿山压力分布、卸压瓦斯运移、瓦斯抽采等规律一体化同步研究。
在模拟实验过程中，利用高频微震监测系统采集震动信号、高清窥视系统探测裂隙发育相结合，得到煤层开采覆岩裂隙演化特征、空隙率分布特征，分析瓦斯最优运移通道。利用在模型关键位置立体化埋深的压力传感器，研究煤层开采过程中矿山压力分布规律及动态变化特征。对模型注入瓦斯并分别进行U，U+L，Y型等方式通风，结合煤层开采时空效应研究覆岩裂隙中瓦斯的浓度、压力、流量分布，揭示卸压瓦斯运移机理;并对模型进行多种方式瓦斯抽采，确定最佳的抽采位置以及不同瓦斯抽采方式的瓦斯抽采效率，为煤与瓦斯安全共采提借鉴。
图11 煤与瓦斯共采三维大尺度物理模拟实验系统
Fig.11 3D large-scale physical simulation experimental system for coal and gas co-extraction
实验选取山西某矿302工作面为原型，工作面走向布置，面长198 m，走向长度2 025 m，地表为丘陵山地。地面标高+1 483～+1 663 m，井下标高+1 152～+1 224 m。302工作面开采煤层为太原组15号，煤厚4.35～5.24 m，煤层倾角3°～14°，平均8.5°。采用综合机械化一次采全高采煤法，全部垮落法管理顶板。
根据几何相似、运动相似、时间相似，结合现场实际情况及实验条件，确定几何相似比1∶100，风速相似比1∶1，其他相似比见表1。
表 1 相似比例
Table 1 Similarity ratio
相似材料模型的制作过程为:将沙子筛分成满足实验要求的粒径，其次秤取预定量的沙子、水泥、石膏和淀粉，按照配比加入适量的水进行充分混合。然后通过相似材料输送装置将材料送入实验箱体中，将材料铺平后，采用自动夯实电锤将材料夯实至预定容重和厚度即可完成一层的铺装。岩层之间均匀铺撒云母片作为分层，如此循环往复可完成模型的制作。在制作模型的过程中，将高频微震传感器、应力传感器、瓦斯浓度传感器、瓦斯压力测试管路、瓦斯抽采管路等按预设位置布置在模型中。
在工作面推进过程中，高频微震监测系统对模型破裂产生的微震信号进行实时采集分析，并对较大的微震信号进行时空定位。
工作面推进至45 m时，采空区覆岩裂隙大量增多，破坏分布范围进一步增大，在覆岩13 m处产生较大破坏，覆岩发生大面积破断现象，形成较多纵向穿层裂隙，破坏一直延伸至28 m高度处，产生若干离层裂隙，覆岩空间破坏分布形态呈现出近似椭圆抛物状。由此时采空区覆岩活动现象可判定为工作面初次来压。
如图12所示工作面推进至64 m时，采空区覆岩在23 m高度处增加较多破坏，覆岩发生较大破断现象，纵向裂隙发育，最大破坏高度增加至33 m，覆岩离层裂隙向上延伸。此阶段破坏范围及强度较大，裂隙在空间上呈椭圆抛物的分布更加明显。此时采空区覆岩活动范围较大，且活动较强烈，为工作面第1次周期来压。
图12 工作面推进64 m微震事件分布
Fig.12 Distribution of microseismic events of working face advancing 64 m
之后工作面每开采20 m左右，采空区覆岩发生周期性破断现象，沿走向破坏范围不断增大，最后垮落高度稳定在60 m，在空间上集中破坏分布基本呈椭球体形态。
采煤工作面在推进过程中，采空区覆岩自下而上依次发生垮落、离层、弯曲下沉等过程，使采场一定范围内应力发生变化，破坏了原有的应力平衡状态，相应的在煤壁前方一定范围内的顶板出现垂直应力集中区和卸压区，在卸压区内的顶板岩层之间及层内产生具有一定规律的裂隙。以应力集中系数(K，煤层现在应力与原始应力之间的比值)来定量描述采动过程中覆岩应力的动态变化。根据应力传感器采集数据，绘制了工作面推进124～200 m范围煤层底板沿走向应力分布规律，如图13～16所示。
图13 工作面推进124～150 m时垂直应力
Fig.13 Vertical stress of working face advancing 124-150 m
图14 推进154～180 m时垂直应力
Fig.14 Vertical stress of working face advancing 154-180 m
图15 推进184～200 m时垂直应力
Fig.15 Vertical stress of working face advancing 180-200 m
图16 工作面推进至200 m应力分布
Fig.16 Vertical stress of working face advancing to 200 m
由图13～16可知，随着工作面推进，垂直应力是动态变化的，受采动影响煤壁前方形成了随工作面推进而不断前移的应力集中现象。在工作面推进过程中应力峰值超前工作面距离6～11 m，应力峰值集中系数为2.11～2.63;超前支撑压力影响范围为47～68 m;未受采动影响区位于工作面前70 m以远，此区受采动影响较小。
当工作面推进至200 m，按照相似比例，对模型进行通风、瓦斯注入，通风方式为U型、风量为0.3 m3/min、瓦斯涌出量0.001 8 m3/min。图17是距离煤层底板5 m处，沿工作面方向距离进风巷道0，40，120，160 m处，沿采空区深部方向2，28，76，124,172，198 m瓦斯体积分数分布情况。
图17 采空区深度方向瓦斯体积分数分布
Fig.17 Gas concentration distribution in depth direction of goaf
由图17可知，距工作面2，28，76 m处瓦斯体积分数依次递增，且受到风流影响增加趋势较为平缓;从76～120 m瓦斯体积分数急剧增加，曲线斜率明显高于前部;在深部采空区，体积分数在32%～35%，瓦斯体积分数较为稳定，这是因为在深部采空区裂隙密度较浅部最大，不会受到风流的影响，且瓦斯解吸速率达到平衡。
图18是U型通风方式下，通风量0.3 m3/min、瓦斯涌出量0.001 8 m3/min时，覆岩裂隙高度方向上瓦斯体积分数的分布情况。瓦斯体积分数传感器按层位布置，分别距离煤层底板5，24，60 m。
图18 覆岩裂隙不同高度瓦斯体积分数分布
Fig.18 Distribution of gas concentration at different heights in overburden cracks
由图18可知，覆岩裂隙中距煤层底板5，24,60 m瓦斯体积分数依次增高，在裂隙带顶端的瓦斯体积分数达到最大值。在60 m高度上进风侧与回风侧、沿着采空区深度方向瓦斯体积分数基本一致，瓦斯体积分数达到65%～68%。这是由于瓦斯在采动裂隙中升浮-扩散、渗流，最终富集在裂隙带顶端。为验证物理模拟的准确性，采用FLUENT数值模拟软件建立模型，得到采空区瓦斯分布规律，模结果如图19所示。从结果来看，物理模拟与数值模拟的采空区中瓦斯体积分数分布基本一致。
图19 U型通风采空区瓦斯分布
Fig.19 Gas distribution in U-type ventilated goaf
(1)煤与瓦斯共采三维物理模拟实验系统采用模块化设计，主要由大尺度箱体与基座、自动液压开采单元、柔性加载单元、自动通风单元、瓦斯注入单元、瓦斯抽采单元以及综合数据采集与控制单元等7部分组成。
(2)实验系统可实现覆岩裂隙演化规律、矿山压力分布规律、卸压瓦斯储运规律、瓦斯抽采规律等一体同步模拟研究。
(3)运用该系统进行验证性实验，得到山西某矿302工作面基本顶初次来压步距45 m，周期来压步距20 m，覆岩破坏在空间上呈椭圆抛物形态;工作面推进过程中应力峰值超前工作面距离6～11 m，集中应力系数2.11～2.63;采空区后部76～120 m瓦斯体积分数增加较快，120 m之后趋于稳定，采空区上部5～60 m裂隙带中瓦斯体积分数逐渐增加，裂隙带最上层瓦斯体积分数达到65%～68%，所得结果采用数值模拟软件进行验证，其结果基本一致。
参考文献
[1] 谢和平，鞠杨，高明忠，等.煤炭深部原位流态化开采的理论与技术体系[J].煤炭学报，2018，43(5):1210-1219.
XIE Heiping，JU Yang，GAO Mingzhong，et al.Theories and technologies for in-situ fluidized mining of deep underground coal resources[J].Journal of China Coal Society，2018，43(5):1210-1219.
[2] 袁亮.我国深部煤与瓦斯共采战略思考[J].煤炭学报，2016，41(1):1-6.
YUAN Liang.Strategic thinking of simultaneous exploitation of coal and gas in deep mining[J].Journal of China Coal Society，2016，41(1):1-6.
[3] 袁亮.煤炭精准开采科学构想[J].煤炭学报，2017，42(1):1-7.
YUAN Liang.Scientific conception of precision coal mining[J].Journal of China Coal Society，2017，42(1):1-7.
[4] 祁和刚.深部高应力巷道综合卸压技术研究与实践[J].采矿与安全工程学报，2016，33(6):1023-1029.
QI Hegang.Research and practice on integrated pressure releasing technology in deep coal mine rock roadway under high stress[J].Journal of Mining & Safety Engineering，2016，33(6):1023-1029.
[5] 杨仁树，张宇菲，王梓旭，等.新型地质力学模型实验系统的研制与应用[J].煤炭学报，2018，43(2):398-404.
YANG Renshu，ZHANG Yufe，WANG Zixu，et al.A newly-built geomechanical model test system and its application[J].Journal of China Coal Society，2018，43(2):398-404.
[6] 王炯，朱道勇，何满朝，等.切顶卸压自动成巷岩层运动规律物理模拟实验[J].岩石力学与工程学报，2018，371(11):2537-2547.
WANG Jiong，ZHU Daoyong，HE Manchao，et al.Physical simulation experiment on the movement of rock strata upon automatic roadway forming by roof cutting and pressure releasing[J].Chinese Journal of Rock Mechanics and Engineering，2018，371(11):2537-2547.
[7] 张强勇，李术才，郭小红，等.铁晶砂胶结新型岩土相似材料的研制及其应用[J].岩土力学，2008，29(8):2126-2130.
ZHANG Qiangyong，LI Shucai，GUO Xiaohong，et al.Research and development of new typed cementitious geotechnical similar material for iron crystal sand and its application[J].Rock and Soil Mechanics，2008，29(8):2126-2130.
[8] 李仲奎，卢达溶，中山元，等.三维模型实验新技术及其在大型地下洞群研究中的应用[J].岩石力学与工程学报，2003，22(9):1430-1436.
LI Zhongkui，LU Darong，NAKAYAMA H，et al.New technology and application of 3D model test for large scaled underground group caverns[J].Chinese Journal of Rock Mechanics and Engineering，2003，22(9):1430-1436.
[9] 李仲奎，王爱民，王克忠,等.地下工程三维地质力学模型制作逆向控制技术研究与应用[J].岩石力学与工程学报，2009，28(9):1730-1734.
LI Zhongkui，WANG Aimin，WANG Kezhong，et al.Study and application of inverse controlling technology for 3D geomechanical model construction of underground engineering[J].Chinese Journal of Rock Mechanics and Engineering，2009，28(9):1730-1734.
[10] 胡耀青，赵阳升，杨栋，等.带压开采顶板破坏规律的三维相似模拟研宄[J].岩石力学与工程学报，2003，22(8):1239-1243.
HU Yaoqing，ZHAO Yangsheng，YANG Dong,et al.Study of 3D simulation on breakage for coal floor in mining above aquifer[J].Chinese Journal of Rock Mechanics and Engineering，2003，22(8):1239-1243.
[11] 胡耀青，赵阳升，杨栋.三维固流藕合相似模拟理论与方法[J].辽宁工程技术大学学报，2007，26(2):204-206.
HU Yaoqing，ZHAO Yangsheng，YANG Dong.Simulation theory & method of 3D solid-liquid coupling[J].Journal of Liaoning Technical University，2017，26(2):204-206.
[12] 伍永平，来兴平，曹建涛，等.多场耦合下急斜煤层开采三维物理模[J].西安科技大学学报，2009，29(6):647-653
WU Yongping，LAI Xingping，CAO Jiantao,et al.3D physical model of steep coal seam mining under multi-field coupling[J].Journal of Xi’an University of Science and Technology，2009，29(6):647-653.
[13] 来兴平，伍永平，曹建涛，等.复杂环境下围岩变形大型三维模拟实验[J].煤炭学报，2010,35(1):31-36.
LAI Xingping，WU Yongping，CAO Jiantao,et al.Experiment on rock-mass deformation of large scale 3D-simulation in complex environment[J].Journal of China Coal Society，2010，35(1):31-36.
[14] 单鹏飞.矿山高陡边坡稳定性三维物理模拟实验研究[D].西安:西安科技大学2013.
SHAN Pengfei.Research on stabilization of high and steep slopes in open-pit mines based on three-dimensional physical simulation experiments[D].Xi’an:Xi’an University of Science and Technology,2013.
[15] 魏宗勇，李莉，李树刚，等.覆岩采动裂隙中瓦斯运移三维实验台的研制与应用[J].煤矿安全，2015,46(7):5-8.
WEI Zongyong，LI Li，LI Shugang,et al.Research and development of three dimensional test bench on gas migration in the mining fissure of overburden rock and its application[J].Safety in Coal Mines，2015，46(7):5-8.
[16] 李树刚，赵鹏翔，林海飞，等.煤岩瓦斯“固-气”耦合理模拟相似材料特性实验研究[J].煤炭学报，2015，40(1):80-86.
LI Shugang，ZHAO Pengxiang，LIN Haifei，et al.Study on character of physical simulation similar material of coal-rock and gas solid-gas coupling[J].Journal of China Coal Society，2015，40(1):80-86.
[17] 李树刚，林海飞，赵鹏翔，等.采动裂隙椭抛带动态演化及煤与甲烷共采[J].煤炭学报，2014，39(8):1455-1462.
LI Shugang，LIN Haifei，ZHAO Pengxiang，et al.Dynamic evolution of mining fissure elliptic paraboloid zone and extraction coal and gas[J].Journal of China Coal Society，2014，39(8):1455-1462.
[18] 李术才，李清川，王汉鹏，等.大型真三维煤与瓦斯突出定量物理模拟试验系统研发[J].煤炭学报，2018，43(S1):121-129.
LI Shucai，LI Qingchuan，WANG Hanpeng，et al.A large-scale three-dimensional coal and gas outburst quantitative physical modeling system[J].Journal of China Coal Society，2018，43(S1):121-129.
[19] 刘东，许江，尹光志，等.多场耦合煤矿动力灾害大型模拟试验系统研制与应用[J].岩石力学与工程学报，2013，32(5):967-975.
LIU Dong，XU Jiang，YIN Guangzhi，et al.Development and application of multifield coupling testing system for dynamic disaster in coal mine[J].Chinese Journal of Rock Mechanics and Engineering，2013，32(5):967-975.
[20] 袁瑞甫，李怀珍.含瓦斯煤动态破坏模拟实验设备的研制与应用[J].煤炭学报，2013，38(1):118-123.
YUAN Ruifu，LI Huaizhen.Development and application of simulation test apparatus for gassy coal dynamic failure[J].Journal of China Coal Society，2013，38(1):118-123.