3.安徽理工大学 地球与环境学院，安徽 淮南232001)摘要：为研究断层对采动应力演化控制作用及发生滑移特征，利用FLAC3D软件对断层带两盘 的接触状态以库仑剪切模型的接触面单元代替，模拟随着工作面沿断层倾向和沿断层走向推进
时断层接触面上法向应力和剪切应力的演化特征、断层滑移危险性演化规律。结果表明：断层面 上法向应力和剪切应力的演化具有明显的时空特征，且断层面上剪切应力对采动影响的敏感性 强于法向应力；当工作面沿断层倾向或断层走向推进时，开采层上方断层先于开采层下方断层
“活化”，开采层上方附近断层在采空区距离断层20 m时最容易出现“活化”，开采层下方断层在 采空区距离断层小于10 m时滑移危险性最高；相比之下，工作面沿断层走向推进时开采扰动对
开采层附近断层产生的影响较工作面沿断层倾向推进时大。关键词：逆断层；法向应力；剪切应力；滑移危险性；断层“活化”；采动应力中图分类号:TD325 文献标志码：A 文章编号：1003-496X(2019) 12-0174-05Numerical Simulation Study on Influence of Faults on Mining StressDOU Zhongsi1, TIAN Nuocheng2, WU Jiwen3(1 .School of Earth Sciences, East China University of Technology, Nonehang 330013, China;2. School of Resource and
Environmentxd Engineering, Hefei University of Technology, Hefei 230009, China\\3.School of Earth and Environment, Anhui
University of Science and Technology, Huainan 232001, China)Abstract: To study the controlling effect of faults on the evolution of mining stress and the characteristics of slip occurrence, an interface element with a Coulomb shear model was used to simulate the interface at the fault by using FLAC3D numerical
simulation software and the processes of face mining along fault trend and fault dip were simulated respectively. The development
of normal stress and shear stress as well as slip risk at the fault was studied. Simulation results showed that: the development of
normal stress and shear stress at the fault has obvious spatial -temporal features and the sensitivity of shear stress on the fault
surface is stronger than the normal stress. Fault above mining layer activates sooner than fault below mining layer when face
advances along the fault trend or fault dip. Fault near the upper side of mining layer activates easiest when the distance of mined-
out area and fault is 20 m and the slip risk of fault near the below of mining layer is highest when the distance of mined-out area and fault is less 10 ni. By contrast with face advancing along the fault trend, mining disturbance makes a larger effect on fault
near mining layer when face advances along the fault dip.Key words: inverse fault; normal stress; shear stress; slip risk; fault activation; mining stress近年来，对断层“活化”诱因研究方法主要有理 力偏转影响而减小。论解析冋擞值模拟Z、相似模拟UW3种方法。赵 毅鑫等何采用库仑破裂应力增量为断层单元活化判
李鹏飞等阿通过FLAC3D数值软件分析断层倾 角、深度、破碎带宽度等对断层“活化”活化，表明： 倾角越小、采深越大、破碎带越宽、落差大的断层更 容易出现断层活化。吴基文等M基于相似物理材料 模拟试验得出断层“活化”是当断层煤柱宽度减小到
据，引入损伤因子评价断层面活化滑移破坏程度， 研究表明：采动诱发主应力偏转使得沿断层面剪应
力增加，而垂直于断层的正应力又受开采卸压和应一定程度时，断层面上、下盘岩体会发生相对错动。
基金项目：江西省数字国土重点实验室开放基金资助项目(DLLJ201903)崔晓荣等冋以龙东煤矿F张断层为原型进行相似模
拟实验,表明:随着工作面的推进，超前支撑应力在
•174-第50卷第12期2019年12月Safety in Coal MinesVol.50 No. 12Dec. 2019煤层底板下30 m处反应明显；当煤柱宽度减小至
30 m时，底板集中应力达到最大值127.64 kPa,应 力集中系数2.8；煤柱宽度20 m时，煤层底板深部
近断层带附近开始明显卸压，同时断层带两侧应力 变化岀现明显差异，显示出断层带对应力传递的影
响，上下盘表现出不同的采动特征，容易造成断层 的“活化”。但是这些研究工作缺乏对开采过程中断 层面上法向应力和剪切应力的变化规律的分析，尤
其是工作面沿逆断层走向开采时断层面上法向应力
和剪切应力演化规律的分析研究。而目前解析法和 相似物理材料模拟法还难以实现对开挖过程中断层
面应力状态的理论计算或仪器监测。1模型建立与计算方案1.1数值模型建立建立的三维工程地质模型如图1,模拟长(％向)
373 m、宽(y向)250 m、高(z向)150 m,共划分为
562 500个单元,592 722个节点。Block Group灰岩泥岩砂岩7?煤粉砂岩Interface locations图1三维数值模型
边界条件:模型前后左右固支，限制水平运动; 模型底面固支，限制垂向运动;模型顶面为自由面,
不予约束。在模型顶面施加等效荷载，模拟上覆岩 土层的自重荷载，在水平方向上施加由自重应力产
生的侧向应力。断层模拟方法:考虑到逆断层断层带宽度较小,
断层带处采用接触面构建的。主要参数为：法向刚 度5.9 GN/m,剪切刚度7.1 GN/m,内摩擦角22。,黏 聚力0.5 MPao模型构建是基于莫尔-库伦判别准则
来判识岩体的破坏。各岩层的物理力学参数见表1。l. 2计算方案当工作面沿断层倾向开采时，以距离模型左侧 边界50 m作为约束煤柱，工作面由远及近向断层
依次开采，直至开采接近断层。工作面宽度设为 150 m,初次来压步距为30 m,周期来压步距为20m, 初次填充步距为30 m,周期填充步距为20 m,工表1岩体物理力学参数岩性体积模 剪切模 内摩擦 黏聚力抗拉强 密度心量/GPa量/GPa角/(。)/MPa度/MPacm'3)泥岩8.804.50281.000.852.485砂岩22.0024.00406.516」22.6817?煤3.331.11251.000.151.400粉砂岩13.0011.00342.802.362.619灰岩17.901&90407」04.202.600作面沿乂向共推进270 m,考察工作面由远及近向 断层推进过程中断层面应力、位移演化规律。当工作面沿断层走向开采时，以距离模型前侧
边界50 m作为约束煤柱，工作面宽度仍设为150
m,为便于与工作面沿断层倾向开采对比，工作面沿
y向共推进150 m,填充步距同上，考察留设不同宽 度煤柱下断层面应力、位移演化规律。监测点布置在断层面y向中心线开采层上方
10.50 m和开采层下方10 m处，分别标记为A、B、C
监测点。2沿断层倾向开采断层采动效应分析2.1断层采动应力分析工作面沿断层倾向开采时应力演化规律如图2。s152s
9 6♦A测点
3测点
亠C测点) 60 90工作面距断层距离/m(a )法向应力0
30 60 90 120 150工作面距断层距离/m(h)剪切应力图2工作面沿断层倾向开采时应力演化规律工作面由远及近向断层推进时断层面上监测点
法向应力演化规律如图2(a)0可以看出三监测点法 向应力演化规律具有明显不同的曲线形态，说明断 层面上法向应力的变化具有明显的时空特性。随着
•175-第50卷第12期2019年12月Safety in Coal MinesVol.50 No. 12Dec. 2019工作面的推进，开采层上方50 m的A监测点法向
应力总体呈先逐渐减小后趋于平缓的趋势，当工作
面推进至距离断层10 m时，法向应力达到最小值4.5 MPa,较初始值降低约67%;开采层上方10 m
的B监测点法向应力总体呈先缓慢增加后快速减 小的趋势，在工作面距断层120-150 m范围时法向
应力约为13.9 MPa,在工作面距离断层30 m时达 到最大值15.83 MPa,应力集中系数为1.14,而后法
向应力快速减小，应力得到释放，工作面接近断层 时法向应力降低至1.9 MPa,且小于A监测点在工
作面接近断层时的法向应力值；开采层下方10 m 的C监测点法向应力在工作面距断层150-50 m范
围时法向应力约为13.65 MPa,在工作面距离断层
10 m时达到最小值12.83 MPa,较初始值降低约
6%,而后增加到最大值13.8 MPa。工作面由远及近向断层推进时断层面上监测点 剪切应力演化规律如图2(b)。A监测点剪切应力呈
先增加后减小再增加的趋势，在工作面距离断层80
m时剪切应力达到最大值3.66 MPa,剪切应力集中
系数为1.22,在工作面距离断层30 m时达到最小
值0.32 MPa,而后逐渐增加至2.46 MPa;B监测点 剪切应力演化规律较为复杂，在工作面距断层40~
150 m范围时剪切应力总体呈逐渐减小的趋势，而 后逐渐增加，在工作面距断层20 m时达到最大值
3.45 MPa,略高于其原岩应力水平，随后剪切应力快 速减小，应力快速释放，在工作面距断层10 m时达 到最小值0.44 MPa;C监测点剪切应力在工作面距
断层70-150 m范围时约为3.21 MPa,表明受开采 活动影响较小，之后剪切应力随着工作面推进而呈 现逐渐增加趋势，在接近断层时达到最大值6.36
MPa,应力集中系数为1.98。对比图2(a)和图2(b)可以看出：法向应力曲线
(a)的变化幅度明显小于剪切应力曲线(b)变化幅
度，以C监测点为例，其法向应力在工作面距离断
层距离70 m时基本无变化，宜至距50 m时开始出 现轻微变化，而其剪切应力在工作面距离断层70 m
时开始出现明显变化。这充分说明随着工作面推
进，距断层一定范围之内，断层带上的剪应力先收
到影响，随后法向应力才受到一定影响。2.2断层滑移程度分析依据摩尔-库伦定律知：接触面的摩擦在很大 程度上由剪切应力与法向应力的比值决定的，因此
将模拟断层面上监测点剪切应力与法向引力的比值 作为衡量断层滑移危险性的指标，此指标值越大，
•176-说明断层滑移危险性越大，反之亦反。随着工作面推进，监测点剪切应力与法向应力
的比值总体演化规律与监测点的剪切应力随工作面 距断层距离的总体演化规律近似，工作面沿断层倾 向开采时剪切应力与法向应力比值如图3O图3工作面沿断层倾向开采时剪切应力与法向应力比值随着工作面推进A监测点变化曲线趋势是先
增加后减小再增加，在工作面距离断层60 m时，比 值首次达到最大值0.37,此时此处的断层最易于首 次出现滑动，出现断层活化的可能性极高；B监测点
剪切应力与法向应力的比值在工作面距离断层20
m时达到极大值0.253,此时此处的断层滑移危险性 最高，随后比值快速减小，在工作面距离断层10 m 时达到最小值0.06,在工作面接近断层时比值快速
增大至0.69,且大于A监测点在工作面接近断层时
的比值;C监测点切应力与法向应力的比值随工作 面的推进呈现逐渐增加的趋势，说明整个过程中剪
切应力增加的趋势占主导，在工作面距离断层10 m
时达到最大值0.46后趋于稳定，表征C监测点在工
作面距离断层10 m时最容易出现滑移，断层活化 的可能性最高。随着工作面沿断层倾向推进，最先出现断层活
化的部位在开采层上方距开采层较远处，之后开采 层上方距开采层较近处活化危险性增大，开采层下
方距开采层较近处在工作面距离断10 m时活化危 险性最高。工作面沿断层倾向开采时竖向滑移量演化规律
如图4。图4中曲线反应了与图3中相同的结论，开 采层上方50 m的A监测点先于开采层上方10 m 的B监测点和开采层下方10 m的C监测点出现滑
动，在工作面接近断层时，监测点C的竖向位移量 明显小于A监测点和B监测点的竖向滑移量。3沿断层走向开采断层采动效应分析3.1断层采动应力分析第50卷第12期2019年12月Safety in Coal MinesVol.50 No. 12Dec. 2019o
图4工作面沿断层倾向开采时竖向滑移量演化规律工作面沿断层走向开采时应力演化规律如图5。81 B5I 1
2 、9-R词 6
■A测点-z ♦B测点 也3 亠C测点o
20
40 60 80 100
120断层煤柱宽度/m(a)法向应力♦A测点 十1＜测点 亠C测点断层煤柱宽度/m(h)剪切应力图5工作面沿断层走向开采时应力演化规律工作面沿断层走向推进时监测点法向应力演化
规律与工作面由远及近向断层推进时法向应力变化
规律十分相似，但也有细微差别，如图5(a)o A监测
点的法向应力在断层煤柱宽度为20 m时达到极小 值4.34 MPa,略小于工作面沿断层倾向开采距离断 层20 m时此监测点的法向应力，随着保留煤柱厚 度的减小，法向应力呈逐渐增加趋势，与工作面距
离断层20 m范围内法向应力趋于平缓的趋势略有 不同；B监测点的法向应力在断层煤柱宽度为30 m 时达到最大值15.79 MPa,应力集中系数为1.11,略 小于工作面沿断层倾向开采距离断层30 m时此监
测点的应力集中系数，而后法向应力快速减小，应
力得到释放，断层煤柱宽度为0 m时法向应力降低 至3.63 MPa,小于A监测点在断层煤柱宽度为0 m 时的法向应力值，大于工作面沿断层倾向开采接近
断层时的法向应力值;C点法向应力在工作面距断
层50-150 m范围内基本相等，基本未受到工作面
回采的影响，在断层煤柱宽度为10 m时达到最小 值13.23 MPa,略大于工作面沿断层倾向开采距离 断层10 m时此监测点的法向应力，而后增加到
13.98 MPa,接近其原岩应力水平。监测点剪切应力演化规律与工作面由远及近向 断层推进时监测点剪切应力演化规律也十分相似，
但也有细微的差别，如图5(b)0 A监测点在断层煤
柱宽度为80 m时剪切应力达到最大值3.49 MPa,
略大于工作面沿断层倾向开采距离断层80 m时的 剪切应力，在断层煤柱宽度为30 m时达到最小值
0.21 MPa,略小于工作面沿断层倾向开采距离断层
30 m时的剪切应力；B监测点在断层煤柱宽度为20
m时剪切应力达到最大值5.08 MPa,远高于其原岩 应力水平，随后剪切应力快速减小,应力快速释放， 在断层煤柱为0 m时达到最小值2.00 MPa,与工作
面沿断层倾向开采距离断层0~20 m范围时演变规 律不同；当断层煤柱宽度大于90 m时，剪切应力基 本无变化，稳定在3.20 MPa,说明此范围内基本不 受工作面开采影响，随着断层煤柱的减小，剪切应
力逐渐增加，在断层煤柱为0 m时达到最大值6.43
MPa,略大于工作面沿断层倾向开采接近断层时的
剪切应力。同样对比断层面上监测点的法向应力和
剪切应力演化规律曲线可以看岀，断层面上监测点
的剪切应力对工作面回采影响的敏感性略强于法向 应力。3.2断层滑移程度分析工作面沿断层走向开采时剪切应力与法向应力
比值如图6O图6工作面沿断层走向开采时剪切应力与法向应力比值顺断层倾向、逆断层倾向2种方式开采下，剪切 应力与法向应力的比值随工作面距断层距离的演化 规律总体近似。A监测点剪切应力与法向应力的比
值随断层煤柱的减小呈现先增加后减小再增加最后•177-第50卷第12期2019年12月Safety in Coal MinesVol.50 No. 12Dec. 2019缓慢减小的趋势，与其剪切应力的演化规律略有差 别，在断层煤柱宽度为60 m时，比值首次达到极大 值0.35，同样此时的断层最易于首次出现滑动;B监 测点剪切应力与法向应力的比值在断层煤柱宽度为
20 m时达到极大值，此时此处的断层滑移危险性最 高，随后减小再快速增加，与工作面沿断层倾向开 采距离断层10 m时不同，在断层煤柱宽度为10 m 时，比值虽达到极小值，但比值仍明显大于其初始 值。C监测点切应力与法向应力的比值随断层煤柱
宽度减小呈现的变化趋势与其剪切应力的变化趋势 近似，比值在断层煤柱宽度为10 m处达到最大值
0.45后趋于稳定，说明C监测点在断层煤柱宽度为
10 m时最容易出现滑移。随着工作面沿断层倾向推进，最先出现断层活 化的部位在开采层上方距开采层较远处，随后开采 层上方距开采层较近处活化危险性增大，开采层下
方距开采层较近处在工作面距离断层10 m时活化 危险性最高。工作面沿断层走向开采时竖向滑移量演化规律
如图7。图7中曲线反应了与图6中相同的结论，开 采层上方50 m的A监测点先于开采层上方10 m
的B监测点和开采层下方10 m的C监测点出现滑
动，在工作面接近断层时,A监测点和B监测点的最 终滑移量十分接近，并远大于C监测点的竖向滑移
量。对比工作面沿断层倾向开采和沿断层走向开采 监测点竖向滑移量可以发现，在采空区与断层距离 相同的条件下，沿断层走向开采时开采层上方10 m
的监测点和开采层下方10 m的C监测点竖向滑移
量略大于沿断层倾向开采时B、C监测点的竖向滑 移量，说明沿断层走向推进时工作面回采对开采层 附近断层产生的影响较工作面沿断层倾向推进大。图7工作面沿断层走向开采时竖向滑移量演化规律4结论1)无论是沿断层倾向开采还是沿断层走向开•178-采,开采层上方断层先于开采层下方断层活化,开采 层上方附近断层在采空区距离断层20 m时最容易 出现活化,开采层下方断层在采空区距离断层10 m
时滑移危险性最高。2) 随着工作面推进，断层带上剪切应力先于法 向应力出现变化。3) 工作面沿断层走向推进对开采层附近断层产
生的影响较工作面沿断层倾向推进时大，为采场附
近巷道布置及断层保护煤柱提供理论依据。参考文献：[1]
鲁海峰，沈丹，姚多喜，等.断层影响下底板采动临界 突水水压解析解[J].采矿与安全工程学报,2014(6)：
887-895.[2] Jiang B, Huang D, Zhu X, et al. Canonical variate
analysis-based contributions for fault identification [ J ].
Journal of Process Control, 2015, 26: 17-25.[3] 林远东，涂敏，刘文震，等.基于梯度塑性理论的断层 活化机理[J].煤炭学报,2012,37(12):2060-2064.[4] 张卫强，孙强，朱术云，等.采动条件下顶板逆断层失 稳-活化机制分析[J].煤矿安全,2013,44(3):179.[5]
姜耀东，王涛，赵毅鑫，等.采动影响下断层活化规律 的数值拟研究[J].中国矿业大学学报,2013,42(1)：
1-5.[6]
赵华安，赵光明，孟祥瑞.断层影响下合理煤柱尺寸留 设数值模拟分析[J ].中国安全生产科学技术,2014,10
( 11 )27-33.[7]
曹兰柱，王珍，王东，等.露天矿含断层逆倾软岩边坡 滑移模式及稳定性研究[J].安全与环境学报,2018,
18(2):457-461.[8] 李凯,茅献彪，陈龙，等.采动对承压底板断层活化及 突水危险性的影响分析[J].力学季刊,2011,32(2)：
261-26&[9] Adams M, Kempka T, Chabab E, et al. An innovative
computationally efficient hydromechanical coupling approach for fault reactivation in geological subsurface uti-
lization[J]. Computers & Geosciences,2017: [s.n.].[10] 张周鑫，刘启蒙，刘丽红.大断层采动“活化”导水特
性研究[J].中国安全生产科学技术,2014,10(7)：
44-4&[11] 徐晓惠，吕进国，刘闯，等.采动影响下逆断层特征参 数对断层活化的作用规律[J].重庆大学学报，
2015,38(3):107-115.[12]
李振华，翟常志，李龙飞.带压开采煤层底板断层活 化突水机理试验研究[J].中南大学学报(自然科学
版),2015,46 (5):1806-1811.(下转第183页)第50卷第12期Vol.50 No. 122019年12月Safety in Coal MinesDec. 2019垩系下统志丹群孔隙、裂隙承压水是工作面的宜接 移动规律相似模拟[J].西安科技大学学报,2013,
充水水源，回采前应提前做好充足防治水准备。33(3):249-253.[11]
田坤云，孙文标,魏二剑•上保护层开采保护范围确 参考文献：定及数值模拟[J].辽宁工程技术大学学报(自然科
[1]
武强.我国矿井水防控与资源化利用的研究进展、问 学版),2013,32(1):7-13.题和展望[J].煤炭学报,2014,39(5):795-805.[12]
栾元重，李静涛,班训海，等.近距煤层开采覆岩导水 [2] 张志龙，高延法，武强，等.浅谈矿井水害立体防治技
裂缝带高度观测研究[J].采矿与安全工程学报，
术体系[J]•煤炭学报,2013,38(3):378-383.2010,27(1):139-142.[3]
董书宁.对中国煤矿水害频发的几个关键科学问题的 [13]
赵高博，郭文兵,杨达明，等.综放开采覆岩破坏模型 探讨[J].煤炭学报,2010,35(1):66-71.及导水裂隙带高度研究[J].中国安全科学学报，
[4]
虎维岳.深部煤炭开采地质安全保障技术现状与研究
2017,27(11):144-149.方向[J].煤炭科学技术,2013,41(8):1-6.[14] 殷伟，张强，韩晓乐，等.混合综采工作面覆岩运移规 [5]
武强，赵苏启，孙文洁,等.中国煤矿水文地质类型划 律及空间结构特征分析[J].煤炭学报,2017,42(2)：
分与特征分析[J].煤炭学报,2013,38(6):901-905.388-396.[6] 张浩,伍永平.慈林山矿煤层群下部采场覆岩运移规 [15] 高喜才，伍永平,曹沛沛，等.大倾角煤层变角度综放
律[J].煤矿安全,2018,49(5):51-54.工作面开采覆岩运移规律[J].采矿与安全工程学
[7] 殷伟，张强，韩晓乐，等.混合综采工作面覆岩运移规 报,2016,33(3):381-386.律及空间结构特征分析[J].煤炭学报,2017,42(2)：
[16]
赵高博，郭文兵,杨达明，等.综放开采覆岩运移及导
388-396.水裂隙带高度研究[J].中国科技论文,2017,12
[8] 吴博文，张宏伟,神文龙，等.麻家梁矿采动双硬顶板 (21):2425-2430.工作面覆岩运移规律[J].煤矿安全,2017,48(10)：
51-54.作者简介：丁 航( 1987—),男，工程师，工学学士，
[9] 樊永山，张胜云.近距离煤层群下行开采下煤层覆岩
2010年毕业于辽宁工程技术大学，现任神东煤炭集团补连
运移规律模拟[J].辽宁工程技术大学学报(自然科学 塔煤矿掘锚一队队长。版),2015,34(8):887-891.[10]
马瑞，来兴平,曹建涛，等.浅埋近距煤层采空区覆岩
(收稿日期：2018-10-11;责任编辑：朱 蕾)(上接第178页)[17] 吴基文，童宏树,童世杰，等.断层带岩体采动效应的 [13 ] Zhang R, Jiang Z, Zhou H, et al. Groundwater out -
相似材料模拟研究[J].岩石力学与工程学报,2007,
bursts from faults above a confined aquifer in the coal 26(S2):4170-4175.mining[j]. Natural Hazards, 2014, 71(3): 1861.[18] 翟晓荣，吴基文,沈书豪，等.断层带边界岩体采动应
[14]
王涛，王矍华,姜耀东，等.开采扰动下断层滑移过程 力特征相似材料模拟研究[J].中国安全生产科学 围岩应力分布及演化规律的试验研究[J].中国矿业 技术,2014(5):56-61.大学学报,2014,43(4):588-592.[15] 赵毅鑫，卢志国，朱广沛，等.考虑主应力偏转的采动
作者简介：窦仲«(1978—),男，安徽六安人，副教授，
诱发断层活化机理研究[J].中国矿业大学学报，
博士，2017年毕业于安徽理工大学，主要从事煤矿工程地 2018,47(1):73-80.质教学和研究工作。[16]
李鹏飞，刘启蒙，陈秀艳.断层活化的主控因素关联 度分析[J].矿业安全与环保,2018,45(5):109-114.(收稿日^：2019-04-19;责任编辑：陈 洋)•183-
因篇幅问题不能全部显示，请点此查看更多更全内容