我国煤层赋存条件差异大、整体透气性差,瓦斯直接抽采效果不佳,因此卸压瓦斯抽采显得十分重要[1-3]。而采动覆岩压实区对卸压瓦斯抽采系统的各项布置参数具有重要的影响,掌握其相关规律有助于提高抽采系统的瓦斯治理效果。目前,围绕采后岩层变化特征及其最终形态,李树刚等[4-6]通过物理相似模拟及数值模拟发现,煤层开采之后采动覆岩裂隙空间形态可以用椭抛带表征。伍永平等[7-9]、来兴平等[10-12]、罗生虎等[13-14]建立了大倾角煤层走向长壁开采岩体结构空间模型,对非对称煤岩的滑落失稳,进而诱发动力学灾害的机理进行了深入分析,并确定了维持大倾角工作面稳定性的区段煤柱合理尺寸。徐玉胜等[15]结合卸荷岩体力学与损伤力学分析研究了采高对采空区顶板卸荷及瓦斯通道演化的影响,得到了不同采高下瓦斯通道的卸荷损伤范围。刘洪永等[16]结合数值仿真对不同推进速度下的优势瓦斯通道的诱导与控制作用进行研究,建立了采动优势瓦斯通道带的时空形态理论模型。基于关键层理论,吴仁伦等[17]就工作面宽度对煤层群开采瓦斯卸压运移“三带”范围的影响进行了深入研究,提出不同的工作面宽度将极大地影响覆岩关键层的移动、破坏形态,进而进一步影响瓦斯卸压运移“三带”范围。张春雷等[18]结合UDEC数值模拟对近距离煤层群开采不同层间岩层结构下围岩裂隙演化规律进行了研究,发现亚关键层的层数、位置会影响覆岩裂隙的发育和分布,但不会影响裂隙的发展高度。7目前从上覆岩层的厚度[19-20]、强度[21-23]以及关键层层位[24]、工作面的采高[25-26]、推进速度[27]、推进距离[28]等方面对采动覆岩裂隙的演化规律及机理的研究较多,但煤层倾角这一对覆岩垮落形态及规律的重要影响因素开展研究较少。以山西和顺某高瓦斯煤矿主采工作面为试验原型,研究在不同煤层倾角条件下仰斜综采工作面覆岩裂隙网络压实区的演化规律,并结合采动覆岩裂隙椭抛带理论,构建考虑煤层倾角的采动覆岩裂隙演化模型,为仰斜综采工作面卸压瓦斯富集区精准探测提供一定的基础。
试验以山西和顺某高瓦斯煤矿主采工作面为原型,开采水平为+1 150 m,平均埋深410 m,煤层平均厚度5.1 m,工作面采用走向长壁综合机械化一次采全高采煤法。直接顶为泥岩,局部含有砂质;基本顶为中砂岩。试验采用平面走向模型,煤岩层物理力学参数及上覆岩层高度见表1、表2。
表1 煤岩层物理力学参数
Table 1 Physical-mechanical parameters of coal and rock
序号岩性容重/(kN·m-3)弹性模量/GPa抗压强度/MPa泊松比1泥岩20.82020.50.1952砂质泥岩26.45648.80.2783中砂岩26.65065.10.2804细砂岩26.24369.00.2605粉砂岩26.05458.50.2536石灰岩26.54791.20.2307铝质泥岩13.04016.00.2508煤14.61413.50.275
表2 上覆岩层模拟高度
Table 2 Simulation height of overlying strata
序号岩性厚度/m序号岩性厚度/m1砂质泥岩2.41613号煤0.52铝质泥岩0.217K3石灰岩3.0315号煤5.118砂质泥岩1.04泥岩0.519细砂岩0.85中砂岩7.120铝质泥岩0.56粉砂岩3.021细砂岩2.27砂质泥岩3.722砂质泥岩4.5814下号煤0.82312号煤1.09粉砂岩1.024砂质泥岩3.410泥岩1.725粉砂岩0.91114号煤0.82611号煤0.312K2石灰岩5.527K4石灰岩4.713粉砂岩2.828砂质泥岩7.814细砂岩6.129细砂岩4.015砂质泥岩2.030砂质泥岩2.7
试验利用自主研发的变化自角物理相似模拟试验台(图1),搭建相似比例为1∶100的物理模型。首先,在试验台底部均匀铺设适宜的应力传感器,然后,按试验具体要求将模型搭建在其上方,待模型自然风干后,再布置相应的位移测点及顶部配重。根据物理相似模拟的特点,模型需满足几何、容重、应力、时间、强度以及泊松比相似的条件,见表3。
表3 模型相似常数
Table 3 Similarity constant of model
参数长宽高容重应力时间强度泊松比相似常数1001.5150101501.0
图1 物理相似模拟试验模型
Fig.1 Physically similar simulation test model
通过调整试验平台的角度(0°、15°、30°),搭建3个只有倾角不同的试验模型。在模拟工作面推进时,两侧先各留宽10 cm煤柱,然后分别按2、3 cm依次循环推进的形式沿煤层走向方向仰斜回采。每次开采完成后,待岩层活动稳定之后拍照记录,并测量各位移测点的下沉量以及岩层的离层量。同时,通过预先埋设的应力传感器采集煤层底板应力。
工作面回采过程中,采空区上覆岩层内的原始应力遭到破坏,当工作面推进到一定距离时,高层位岩层逐渐产生弯曲断裂,加大了低水平岩层的竖向受力。同时下方悬臂梁随工作面的推进不断变长,导致下部岩层无法承受上部岩层的应力作用而产生周期性垮落,如图2所示。
图2 来压分布
Fig.2 Distribution of periodic pressure
不同煤层倾角条件下周期来压显现时的推进距离不同,0°的平均来压步距为12.5 m,分别是15°以及30°的1.08、1.20倍。压实区是采空区上覆岩层受采动影响垮落堆积后,其中瓦斯渗透、运移明显更困难的局部区域,一般在1~2次周期来压后初步形成。煤层倾角越大,周期来压到来所需推进距离越短,因此,压实区的发育周期不断变短,为之后的缓慢充分的发育提供了更为充足的时间。
覆岩离层量即相邻岩层间的裂隙宽度,是定量描述采动过程中覆岩离层的动态变化的重要指标。覆岩离层量分布如图3所示,在受采动影响的采空区边缘,覆岩的离层量较大,而采空区中部因高层位的岩层破断及垮落,挤压低层位的岩层,迫使离层裂隙压缩变窄,甚至闭合,从而离层量较小。根据文献[6]中对压实区的判定准则,即压实区采动覆岩离层量相对于两侧边界较小且在采空区底板会形成应力集中现象,两者结合可判定压实区边界。因此,结合图3,可大致判定中部的压实区和边界的裂隙区所在区域。
图3 覆岩离层量分布
Fig.3 Abscission layer distribution of rock
由于裂隙区的离层量普遍大于压实区,故在二者交界区域离层量的大小会有较明显的变化。随着煤层倾角的增加,中部压实区和两侧裂隙区所对应的离层量也逐渐增加,但受岩层自身硬度、厚度以及层位的影响,最大离层量的增加值在减小。而在压实区中,受采动影响后不同层位的岩层垮落形成的堆砌状态以及所受到的外力大小及方向都存在一定差别,因此,压实区覆岩离层量的数值会有一定的波动。
通过分析模型底板处应力传感器监测到的回采期间的应力数据,如图4所示,发现采空区底板的应力相较于回采前明显变小,且应力集中系数在0.05~0.45波动。由于中部的压实区直接受到上位岩层向下的重力,而两侧的裂隙区得益于垮落时形成的铰接结构的支撑作用,中部的压实区应力集中系数要比两侧裂隙区的大一些,在3个区域的交界处,其应力集中系数会有较为明显的变化,如图4中的局部放大的区域所示。同时,结合覆岩离层量的分布规律可知,不同煤层倾角压实区的宽度扩展5~7 m。
图4 工作面底板应力分布
Fig.4 Stresses distribution of bottom plater
根据文献[6]中对压实区的判定准则,再结合模拟试验中周期来压具体时段、覆岩裂隙发育情况、以及底板应力分布规律,绘制不同煤层倾角的压实区演化形态(图5)。
图5 不同煤层倾角压实区演化形态
Fig.5 Different mining angle of evolution form of compacted area
由图5可知,随着回采工作面的不断推进,压实区也不断的扩展。周期来压依次显现时,压实区明显在沿工作面方向扩展,同时也在向更高层位发展。在压实区向更高层位发展时,覆岩会形成一定的垮落角,且垮落角在逐渐地减小,使得压实区的左右边界呈椭圆抛物线状。不同煤层倾角条件下所表现出的压实区形态有着明显的差异,主要体现在高度、垮落角以及宽度这3个方面。
如图6所示,压实区的高度随着工作面的推进持续增大,同一推进距离下,煤层倾角越大,压实区的高度越高。压实区高度拟合直线斜率初期较小,这是由于压实区在形成初期,受采动影响的覆岩范围较小,压实区整体规模相差不大,即不同煤层倾角下的压实区的高度在采煤初期差值不大。当工作面推进至一定距离时(55~65 m),拟合直线斜率剧增,如图6中阴影区域所示,这是由于受煤层倾角影响,平均来压步距减小,煤层倾角较大的工作面发生周期来压的次数更多,压实区发育更为充分,导致其高度明显增大。当工作面推进距离达到80 m后,斜率又趋于稳定,不同煤层倾角下的压实区的高度其差变化不再明显,受推进距离的影响减小。
图6 不同煤层倾角压实区高度拟合直线斜率变化趋势
Fig.6 Change trend of height slope in compacted zone of different coal seam inclination angle
由于压实区的高度受到工作面推进距离和煤层倾角的共同影响,为了能定量化描述压实区高度与煤层倾角及工作面推进距离之间的关系,将三者进行多元拟合,得到拟合关系:
H=25.812ln(L-10.057)+0.178θ-69.216,R2=0.949
(1)
式中:H为压实区高度,m;L为工作面推进距离,m;θ为煤层倾角,(°);R2为决定系数。
如图7所示,随着煤层倾角的增大,开切眼侧垮落角α逐渐减小,工作面侧垮落角β逐渐增大,但随工作面推进距离的增加,两侧垮落角都在逐渐减小。煤层倾角越大,下滑作用越强,受采动影响的上覆岩层破断及垮落后沿煤层底板向下滑动或滚落填充下部采空区,且煤层倾角越大,填充下部采空区的岩块堆砌的高度越高,导致开切眼侧的上覆岩层垮落空间相对减少,并且受到下部填充岩块的支撑,使得岩层间形成较大的离层却未发生垮落,而以覆岩离层量为判定依据的开切眼侧压实区边界则向工作面侧进一步弯曲倾斜。因此,对于开切眼侧的垮落角α而言,煤层倾角越大,其值越小。
图7 不同煤层倾角压实区垮落角演化规律
Fig.7 Different of mining angle of evolution of collapse angle
煤层倾角越大,下滑作用越强,岩块堆砌的速度越快,导致填充下部采空区的岩块在滑动或滚落的过程中相互摩擦、碰撞及挤压程度加剧,岩块的单个体积更小,单位空间内堆砌的岩块更多,进而进一步增大了工作面侧上覆岩层的垮落空间,难以得到下部垮落岩块的支撑,使其不易形成较大的离层量,从而工作面侧的压实区边界向开切眼侧的弯曲倾斜程度相对减小。因此,对于工作面侧的垮落角β而言,煤层倾角越大,其值越大。
由于压实区的垮落角受到工作面推进距离和煤层倾角的共同影响,为了能定量化描述压实区垮落角与煤层倾角及工作面推进距离之间的关系,将三者进行多元拟合,得到拟合关系:
α=-0.640ln L-0.003θ+3.576,R2=0.981
(2)
β=-0.711ln L+0.008θ+3.878,R2=0.960
(3)
由于压实区的垮落角是逐渐减小的,这使得压实区的边界存在一定的弧度,但同时受到下滑作用的影响,从而左右边界呈不对称椭抛线状。而压实区形成之后,其扩展的部分都是在前一次的基础上进一步形成的,因此根据式(2)和式(3)所得到的垮落角拟合公式,可以建立如图8所示的简化模型。
图8 走向压实区简化模型
Fig.8 Simplified model of compacted area
进一步得到压实区宽度与垮落角及垮落高度之间的关系。由于中部压实区边界呈不对称椭抛线状,近似认为该模型是由2部分组成的,可得到以下关系:
Li=Lαi+Lβi
(4)
Lαi=ΔHα1cot α1+ΔHα2cot α2+…+
ΔHαncot αn=Hα1cot α1+(Hα2-Hα1)×
cot α2+…+(Hαn-Hαn-1)cot αn=
(5)
则:
(6)
式中:Li为第i次形成的压实区宽度,m;Lαi为第i次形成的压实区开切眼侧的宽度,m;Lβi为第i次形成的压实区工作面侧的宽度,m;ΔHαi为第i次压实区扩展时开切眼侧高度的增量,m;Hαi为压实区第i次扩展时开切眼侧的高度,m;Hβi为压实区第i次扩展时工作面侧的高度,m;αi为压实区第i次扩展时开切眼侧形成的垮落角,(°);βi为压实区第i次扩展时工作面侧形成的垮落角,(°)。
根据所得的压实区高度、垮落角与工作面推进距离和煤层倾角的关系,将其代入式(4)中,便可得到压实区宽度与工作面推进距离和煤层倾角的关系。
通过物理相似模拟试验可以发现,采空区中部的压实区受下滑作用影响而整体呈不对称的类椭圆抛物面。根据试验所得数据,将压实区倾角模型旋转θ(θ为煤层倾角)建立了如图9所示的椭抛带压实区演化模型,同时,基于采动裂隙椭抛带理论,结合式(1)—式(6),构建了受煤层倾角及工作面推进距离共同影响的压实区空间形态数学表征方程。
图9 采动裂隙椭抛带压实区空间几何模型
Fig.9 Geometric model of compacted area in elliptic paraboloid zone
(7)
式中:x1、y1、z1为图9中压实区的坐标值;A1、A2分别为进、回风巷处,压实区底部边界与其距离,m;Lk为工作面宽度,m;La为压实区底部边界至开切眼处的距离,m;Kc为压实区范围内岩层破断碎胀系数;H为压实区高度,m。
1)随着煤层倾角的增大,相同区域内的采动覆岩离层量随之增大,结合覆岩离层量与底板应力分布规律,可知不同煤层倾角下压实区的宽度扩展5~7 m。
2)随着煤层倾角的增大,平均来压步距逐渐减小(12.5 m>11.6 m>10.4 m),工作面推进至55~65 m时,煤层倾角较大的工作面发生周期来压次数更多,压实区之间的高差急剧增大。
3)同一推进距离下,煤层倾角越大,压实区工作面侧垮落角越大,开切眼侧垮落角越小,结合压实区高度演化规律,建立受煤层倾角因素影响的走向压实区简化模型,对压实区宽度分段表示。
4)建立了考虑煤层倾角因素的采动裂隙椭抛带压实区数学方程,为进一步研究受煤层倾角因素影响的采动覆岩裂隙发育规律,准确定位仰斜综采工作面卸压瓦斯富集区提供一定的理论依据。
[1] 袁 亮.煤炭精准开采科学构想[J].煤炭学报,2017,42(1):1-7.
YUAN Liang.Scientific conception of precise coal mining [J].Journal of China Coal Society,2017,42(1):1-7.
[2] 郝全明,柴 敏.近距离煤层群下邻近层增透效应及卸压瓦斯治理技术研究[J].煤炭技术,2018,37(5):141-143.
HAO Quanming,CHAI Min.Study on the anti-reflection effect and the technology of pressure relief gas control in the adjacent coal seam group[J].Coal Technology,2018,37(5):141-143.
[3] 闫 浩,张吉雄,鞠 杨,等.上保护层开采下充实率控制裂隙发育规律及瓦斯抽采研究[J].采矿与安全工程学报,2018,35(6):1262-1268.
YAN Hao,ZHANG Jixiong,JU Yang,et al.Study on the development law of the filling rate and the gas drainage in the upper protective layer mining[J].Journal of Mining and Safety Engineering,2018,35(6):1262-1268.
[4] 李树刚,林海飞,赵鹏翔,等.采动裂隙椭抛带动态演化及煤与甲烷共采[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.
[5] 李树刚,徐培耘,赵鹏翔,等.采动裂隙椭抛带时效诱导作用及卸压瓦斯抽采技术[J].煤炭科学技术,2018,46(9):146-152.
LI Shugang,XU Peiyun,ZHAO Pengxiang,et al.Aging induction and pressure relief gas extraction technology of mining-induced fracture ellipsoid strip [J].Coal Science and Technology,2008,46(9):146-152.
[6] 李树刚,徐培耘,赵鹏翔,等.综采工作面覆岩压实区演化采高效应分析及应用[J].煤炭学报,2018,43(S1):112-120.
LI Shugang,XU Peiyun,ZHAO Pengxiang,et al.Analysis and application of evolutionary mining height effect in overburden compacted zone of fully mechanized mining face [J].Journal of China Coal Society,2018,43(S1):112-120.
[7] 伍永平,皇甫靖宇,解盘石,等.基于大范围岩层控制技术的大倾角煤层区段煤柱失稳机理[J].煤炭学报,2018,43(11):3062-3071.
WU Yongping,HUANGFU Jingyu,XIE Panshi,et al.Instability mechanism of coal pillars in large dip angle coal seam based on large-scale rock formation control technology[J].Journal of China Coal Society,2018,43(11):3062-3071.
[8] 伍永平,胡博胜,解盘石,等.大倾角工作面飞矸冲击损害及其控制[J].煤炭学报,2018,43(10):2694-2702.
WU Yongping,HU Bosheng,XIE Panshi,et al.Impact damage and control of the large dip angle of the flying shoal[J].Journal of China Coal Society,2018,43(10):2694-2702.
[9] 伍永平,郎 丁,解盘石,等.大倾角软煤综放工作面倾斜方向顶煤区划破坏特征[J].采矿与安全工程学报,2018,35(3):553-560.
WU Yongping,LANG Ding,XIE Panshi,et al.Destruction characteristics of top coal division in inclined direction of fully mechanized caving face with large dip angle soft coal[J].Journal of Mining and Safety Engineering,2018,35(3):553-560.
[10] 来兴平,刘简宁,崔 峰,等.急斜煤层顶板裂隙扩展诱导能量时-空演变特征[J].西安科技大学学报,2018,38(4):562-568.
LAI Xingping,LIU Jianning,CUI Feng,et al.Time-space evolution characteristics of induced energy propagation in cracked roof of steep coal seam[J].Journal of Xi’an University of Science and Technology,2018,38(4):562-568.
[11] 来兴平,杨毅然,单鹏飞,等.急斜煤层顶板应力叠加效应致灾特征综合分析[J].煤炭学报,2018,43(1):70-78.
LAI Xingping,YANG Yiran,SHAN Pengfei,et al.A comprehensive analysis of the disaster-causing characteristics of the roof stress in the steep coal seam[J].Journal of China Coal Society,2018,43(1):70-78.
[12] 来兴平,杨毅然,王宁波,等.急斜煤岩体动力失稳时空演化特征综合分析[J].岩石力学与工程学报,2018,37(3):583-592.
LAI Xingping,YANG Yiran,WANG Ningbo,et al.Synthesis analysis of spatio-temporal evolution characteristics of dynamic instability of steep coal-rock mass[J].Chinese Journal of Rock Mechanics and Engineering,2018,37(3):583-592.
[13] 罗生虎,伍永平,王红伟,等.大倾角煤层长壁开采底板非对称破坏形态与滑移特征[J].煤炭学报,2018,43(8):2155-2161.
LUO Shenghu,WU Yongping,WANG Hongwei,et al.Asymmetric failure modes and slip characteristics of longwall mining in large dip angle coal seams[J].Journal of China Coal Society,2018,43(8):2155-2161.
[14] 罗生虎,伍永平,刘孔智,等.大倾角大采高综采工作面煤壁非对称受载失稳特征[J].煤炭学报,2018,43(7):1829-1836.
LUO Shenghu,WU Yongping,LIU Kongzhi,et al.Asymmetric load-bearing instability characteristics of coal wall with large dip angle and large mining height[J].Journal of China Coal Society,2018,43(7):1829-1836.
[15] 徐玉胜,李春元,张 勇,等.不同采高下瓦斯通道卸荷损伤演化及抽采验证[J].煤炭学报,2018,43(9):2501-2509.
XU Yusheng,LI Chunyuan,ZHANG Yong,et al.Evolution and extraction verification of unloading damage in gas channels under different mining heights [J].Journal of China Coal Society,2008,43(9):2501-2509.
[16] 刘洪永,程远平,周红星,等.综采长壁工作面推进速度对优势瓦斯通道的诱导与控制作用[J].煤炭学报,2015,40(4):809-815.
LIU Hongyong,CHENG Yuanping,ZHOU Hongxing,et al.Induction and control effect of advancing speed of longwall working face in fully mechanized mining on dominant gas channel [J].Journal of China Coal Society,2015,40(4):809-815.
[17] 吴仁伦,王亚飞,徐东亮,等.工作面面宽对煤层群开采瓦斯卸压运移“三带”范围的影响[J].采矿与安全工程学报,2017,34(1):192-198.
WU Renlun,WANG Yafei,XU Dongliang,et al.Influence of aspect width on the “three zones”of gas pressure relief and migration in coal seam group mining [J].Journal of Mining and Safety Engineering,2017,34(1):192-198.
[18] 张春雷,李占平,康 强.近距煤层群不同层间岩层结构下围岩裂隙演化规律[J].煤矿安全,2018,49(9):91-95.
ZHANG Chunlei,LI Zhanping,KANG Qiang.The evolution law of surrounding rock fissures under different interlayer rock stratum structures[J].Safety in Coal Mine,2018,49(9):91-95.
[19] WANG W,CHENG Y,WANG H,et al.Fracture failure analysis of hard-thick sandstone roof and its controlling effect on gas emission in underground ultra-thick coal extraction[J].Engineering Failure Analysis,2015,54:150-162.
[20] ZHAO Y,WANG S,ZOU Y,et al.Pressure-arching characteristics of fractured strata structure during shallow horizontal coal mining[J].Tehnicki Vjesnik,2018,25(5):1457-1466.
[21] 张培鹏,蒋力帅,刘绪峰,等.高位硬厚岩层采动覆岩结构演化特征及致灾规律[J].采矿与安全工程学报,2017,34(5):852-860.
ZHANG Peipeng,JIANG Lishuai,LIU Xufeng,et al.Structural evolution characteristics and disaster law of mining overburden in high hard and thick strata [J].Journal of Mining and Safety Engineering,2017,34(5):852-860.
[22] NING J,WANG J,JIANG L,et al.Fracture analysis of double-layer hard and thick roof and the controlling effect on strata behavior:A case study[J].Engineering Failure Analysis,2017,81:117-134.
[23] WANG Xinfeng,GAO Mingzhong,LI Longqin.Time-space coupling law of mining stress,overburden movement and fracture field distribution in deep stope [J].Journal of mining and safety engineering,2016,33(4):604-610.
[24] 吴仁伦.关键层对煤层群开采瓦斯卸压运移“三带”范围的影响[J].煤炭学报,2013,38(6):924-929.
WU Renlun.Influence of key strata on the "three zones" of gas pressure relief and migration in coal seam group mining [J].Journal of China Coal Society,2013,38(6):924-929.
[25] 李树刚,李志梁,林海飞,等.采高对采动裂隙演化规律的影响研究[J].矿业安全与环保,2015,42(5):25-28.
LI Shugang,LI Zhiliang,LIN Haifei,et al.Study on the influence of mining height on the evolution law of mining fissures[J].Mining Safety and Environmental Protection,2015,42(5):25-28.
[26] 吴仁伦,王继林,折志龙,等.煤层采高对采动覆岩瓦斯卸压运移“三带”范围的影响[J].采矿与安全工程学报,2017,34(6):1223-1231.
WU Renlun,WANG Jilin,CHAI Zhilong,et al.Influence of coal seam height on “three zones”of mining-induced overburden gas pressure relief movement [J].Journal of Mining and Safety Engineering,2017,34(6):1223-1231.
[27] 赵鹏翔,卓日升,李树刚,等.综采工作面推进速度对瓦斯运移优势通道演化的影响[J].煤炭科学技术,2018,46(7):99-108.
ZHAO Pengxiang,ZHUO Risheng,LI Shugang,et al.The influence of the speed of fully mechanized mining face on the evolution of the dominant channel of gas migration[J].Coal Science and Technology,2018,46(7):99-108.
[28] 王新丰,高明中,李隆钦.深部采场采动应力、覆岩运移以及裂隙场分布的时空耦合规律[J].采矿与安全工程学报,2016,33(4):604-610.
WANG Xinfeng,GAO Mingzhong,LI Longqin.Spatial and temporal coupling of mining stress,overburden movement and fracture field distribution in deep stope[J].Journal of Mining and Safety Engineering,2016,33(4):604-610.