采空区积水条件下工作面推进速度对煤岩层稳定性影响研究.pdf

第 卷第期中 国 煤 炭V o l N o 年 月C H I N AC O A LA u g u s t 煤矿安全移动扫码阅读引用格式:金珠鹏,徐亚明采空区积水条件下工作面推进速度对煤岩层稳定性影响研究 J中国煤炭,()D O I:/j c n k i c c m J I NZ h u p e n g,XUY a m i n g R e s e a r c ho n t h e i n f l u e n c eo fw o r k i n g f a c e a d v a n c i n gs p e e do n t h e s t a b i l i t yo f c o a la n dr o c ks t r a t au n d e r t h e c o n d i t i o no fw a t e r a c c u m u l a t i o n i ng o a fJC h i n aC o a l,():D O I:/j c n k i c c m 基金项目:黑龙江省“揭榜挂帅”科技攻关项目(Z X J A ,Z X J A )作者简介:金珠鹏(),男,辽宁本溪人,副教授,博士研究生,主要从事矿山压力与岩层控制方面的研究.E m a i l:q q c o m通讯作者:徐亚明(),男,山东枣庄人,硕士研究生,主要从事矿山压力与岩层控制方面的研究.E m a i l:q q c o m采空区积水条件下工作面推进速度对煤岩层稳定性影响研究金珠鹏,徐亚明(黑龙江科技大学矿业工程学院,黑龙江省哈尔滨市,)摘要针对龙湖煤矿 号煤层受采空区积水影响下的安全开采问题,基于流固耦合及弹塑性力学理论,采用室内力学实验、数值模拟等方法,分析在采空区积水条件下,工作面推进速度对煤岩层稳定性的影响规律.结果表明:随着工作面的推进,工作面处应力集中现象越发显著,且推进速度越快应力峰值越大;煤岩层的位移沉降规律经历了积水影响异常区、工作面沉降影响区以及上下煤层沉降贯通区个阶段,工作面推进速度越快煤层开采区域位移沉降量就越大;水压力随着工作面推进速度的增加而略有增大,但不同的推进速度对煤层开采水力渗流影响不大;随着工作面与积水区距离逐渐缩短,最终诱发上部积水异常范围产生塑性区并与下部煤层贯通;工作面推进速度越快开采扰动产生的煤岩层损伤范围就越大,塑性应变峰值同样增大.关键词采空区积水;数值模拟;工作面推进速度;煤岩层稳定性;应力集中;位移沉降中图分类号T D ;T D 文献标志码AR e s e a r c ho nt h e i n f l u e n c eo fw o r k i n gf a c ea d v a n c i n gs p e e do nt h e s t a b i l i t yo f c o a la n dr o c ks t r a t au n d e r t h ec o n d i t i o no fw a t e ra c c u m u l a t i o n i ng o a fJ I NZ h u p e n g,XUY a m i n g(S c h o o l o fM i n i n gE n g i n e e r i n g,H e i l o n g j i a n gU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,H a r b i n,H e i l o n g j i a n g ,C h i n a)A b s t r a c t F o r t h e s a f em i n i n go fN o c o a l s e a ma f f e c t e db yw a t e r a c c u m u l a t i o n i ng o a f i nL o n g h uC o a lM i n e,t h e i n f l u e n c el a wo f t h e a d v a n c i n gs p e e do f t h ew o r k i n g f a c eo n t h e s t a b i l i t yo f c o a l a n dr o c ks t r a t au n d e r t h e i n f l u e n c eo fw a t e r a c c u m u l a t i o n i ng o a fw a sa n a l y z e d b a s e do nf l u i d s o l i dc o u p l i n ga n de l a s t i c p l a s t i c m e c h a n i c st h e o r i e s,i n d o o r m e c h a n i c a lt e s t s,n u m e r i c a ls i m u l a t i o na n do t h e rm e t h o d s T h er e s u l t ss h o wt h a t a s t h ew o r k i n gf a c ea d v a n c e s,t h es t r e s sc o n c e n t r a t i o na t t h ew o r k i n gf a c eb e c o m e sm o r e a n dm o r e s i g n i f i c a n t,a n d t h e s t r e s sp e a k i s i n c r e a s i n gw i t h t h e a d v a n c i n gs p e e d;t h ed i s p l a c e m e n t s e t t l e m e n t l a wo fc o a l a n dr o c ks t r a t ag o e st h r o u g ht h r e es t a g e s:t h ea n o m a l ya r e aa f f e c t e db yw a t e ra c c u m u l a t i o n,t h ei n f l u e n c e da r e ab yt h es e t t l e m e n to f t h ew o r k i n gf a c ea n dt h es e t t l e m e n tp e n e t r a t i o na r e ao ft h eu p p e ra n dl o w e rc o a ls e a m s,a n dt h ed i s p l a c e m e n ts e t t l e m e n to f t h ec o a ls e a m m i n i n ga r e ai si n c r e a s i n gw i t ht h ea d v a n c i n gs p e e d;t h ew a t e rp r e s s u r ei n c r e a s e ss l i g h t l yw i t ht h ea d v a n c i n gs p e e do f t h ew o r k i n g f a c e,w h i l ed i f f e r e n t a d v a n c i n gs p e e dh a s l i t t l e i n f l u e n c e o nh y d r a u l i c s e e p a g e i nc o a l s e a m m i n i n g;a s t h ed i s t a n c eb e t w e e nt h ew o r k i n gf a c ea n dt h ew a t e ra c c u m u l a t i o na r e ag r a d u a l l ys h o r t e n s,t h eu p p e rw a t e ra c c u m u l a t i o na b n o r m a l r a n g e i s f i n a l l y i n d u c e d t op r o d u c e ap l a s t i c z o n e a n dp e n e t r a t ew i t h t h e l o w e r c o a l s e a m;t h e c o a l a n d r o c ks t r a t ad a m a g er a n g eg e n e r a t e db ym i n i n gd i s t u r b a n c ea n dt h ep e a kp l a s t i cs t r a i nb o t h i n c r e a s ew i t ht h ea d v a n c i n gs p e e do f t h ew o r k i n gf a c e K e y w o r d s w a t e ra c c u m u l a t i o ni ng o a f;n u m e r i c a ls i m u l a t i o n;a d v a n c i n gs p e e do fw o r k i n gf a c e;s t a b i l i t yo fc o a la n dr o c ks t r a t a;s t r e s sc o n c e n t r a t i o n;d i s p l a c e m e n t s e t t l e m e n t引言矿井水害已成为煤矿事故中仅次于瓦斯的第二大灾害,而采空区突水占煤矿水害事故的比重非常大.采空区突水具有时间短、水量大、破坏性强等特点.我国煤矿开采历史悠久,由于早期开采方法与技术落后,一些小煤矿为了开采便利,采用巷道式、房柱式采煤方法对下煤层随意进行采掘,产生了大量错综复杂的采空区.采空区形成后,岩层的自然应力状态和平衡状态受到破坏,岩层内部的应力需要重新分布以达到新的平衡状态,这导致采空区上覆岩层变形、垮落、断裂、离层和移动,形成裂缝带.严重时,岩层的裂隙可能贯通含水层,导致水沙溃入工作面,损坏工作面和井下设备,造成人员伤亡.因此,采空区积水是威胁矿井安全生产的因素之一,研究采空区积水条件下工作面推进速度对煤岩层稳定性的影响规律,对矿井的安全生产具有重要意义.笔者以黑龙江龙煤七台河矿业有限责任公司龙湖煤矿(以下简称“龙湖煤矿”)积水勘测项目为工程背景,通过瞬变电磁法对该矿区北 东一采区及南二采区进行探测,分析采空区积水异常区分布情况.煤层概况选取龙湖煤矿、号煤层作为研究对象,通过理论分析、室内力学实验、数值模拟等研究方法,分析了在采空区积水条件下,工作面推进速度对煤岩层稳定性的影响规律.号煤层为已开挖煤层,煤层平均厚度 m,存在大面积采空区积水;号煤层正在施工开采,位于 号煤层下方,煤层平均厚度 m.龙湖煤矿煤种为焦煤,煤层结构简单,煤层厚度变化不大,较稳定.号煤层顶板为粉砂岩、中砂岩,底板为粉砂岩,质地较坚硬,裂隙不发育,对工作面无影响.研究区域地面平均标高 m,号煤层工作面标高 m,可采走向长度 m,倾斜长度 m,倾角.研究区段煤岩层柱状如图所示.图研究区段煤岩层柱状数值模型的建立 煤岩体力学参数测定结合图煤层柱状图,对相应的代表区域进行现场取样.在实验室制作岩石试样后通过室内力学实验,最终获取各煤岩层物理力学参数,见表.煤岩层流固耦合及弹塑性力学理论 渗流 应力耦合理论()基本理论.流固耦合理论在地下工程中是研究流体与岩体相互作用的一种重要理论依据.岩体固体在渗流场的作用下会产生相应的位移及变金珠鹏等:采空区积水条件下工作面推进速度对煤岩层稳定性影响研究 年第期形,同时位移与变形也会影响渗流场.表现场取样物理力学参数实验结果序号煤岩层密度/(gc m)弹性模量/G P a泊松比粘聚力/k P a内摩擦角/()粉砂岩 中砂岩 号煤层 中砂岩 粉砂岩 号煤层 粉砂岩 当应力为拉应力且孔隙压力达到饱和状态时,岩体有效应力是正值.其表达式为:i ji j i jx pw()paKVKS()式中:i j 有效应力;i j 总应力;x 与饱和度以及表面张力相关参数;pw 孔隙水压力;pa 孔隙气压力;i j K r o n e c k e r符号;B i o t系数,一般取;KV 岩石体积的压缩模量;KS 固体颗粒的压缩模量.()平衡方程.在地层结构中,岩土体应力平衡方程可以通过虚功原理描述,也就是说,在特定的时间点,岩土体受到体力和面力的影响,会导致它们之间的虚功保持不变,即:V TddVV uTdfdVS uTdtdS()式中:t 面力;f 体力;V 岩石体积;S 固体颗粒;T 特定时间点;应变;虚应变;u 虚位移.弹塑性力学理论在地下工程中,煤岩层的承载受力通常分为弹性阶段、塑性阶段和损伤破坏阶段,其应力应变演化规律具有明显的非线性特征.当岩石受地层活动影响所受到的应力值低于极限状态时,应力 应变曲线服从胡克定律,为弹性变化阶段;随着应力的持续增大,超过岩石的弹性极限时,它的应力 应变变化不再按照原来的胡克模式发生,这时进入塑性变化阶段.在材料开始屈服以后,煤岩层因受到周边未屈服物质的阻挡而保持稳定,不会立刻崩溃.根据流体力学原理,在这个阶段,材料的应力 应变关系可以用应变增量来描述,其中,弹性应变增量反映了煤岩石材料的可恢复变形能力,而塑性应变增量反映了岩石的损伤和破坏程度.通常,理想的弹性塑料的初始屈服点与破坏点的应力相同,但应变却有所差异.软硬岩土在初始屈服后,会经历塑性阶段,最终形成破坏面.典型岩土应力 应变演化规律如图所示.图煤岩体典型应力 应变曲线 模型建立数值计算参数参照龙湖煤矿 号、号煤层.模型尺寸选定沿煤层走向长度 m,沿工作面倾斜长度 m,模型高度 m,即整体模型尺寸为 m m m(长宽高).根据积水勘测项目探测结果设定 号煤层为已开挖采空区积水煤层,煤层厚度 m,煤层倾角,采空区走向长度 m,中部 m为积水异常区;未开挖 号煤层厚度 m,煤层倾角,工作面倾斜长度 m,沿煤层走向分步开挖,最终挖掘长度 m.为模拟上覆岩层重量,将其重量转化为 M P a等效均布荷载作用于模型顶部.各岩层属性及深度参照 年第期 中 国 煤 炭 第 卷图的煤层柱状,将煤层及各类型岩层分区,得到数值模型网格划分如图所示.图数值模拟网格划分以工作面推进速度m/d工况模型为例,考虑到 采 空 区 积 水 效 应 的 影 响,网 格 单 元 设 定 为A B AQU S用于计算地下孔隙流体应力的八节点六面体单元(C D P),模型共划分单元 个,节点数量 个.对于模型位移边界条件,限制左右两侧X方向的水平位移,模型前后Z方向的位移,以及模型底部Y方向的水平位移.此外,由于上部煤层存在 m积水异常区,通过在该区域施加 MP a的孔隙水压力荷载来实现积水异常的不利荷载影响.该模型采用弹塑性本构模型,结构服从莫尔库仑(M o h r C o u l o m b)准则.模拟方案及数值模拟结果在上部煤层采空区积水影响下,为探究下部煤层开采时工作面不同推进速度对层间煤柱的扰动影响,根据不同的推进速度逐次开挖,模拟工作面推进 m时煤柱破坏变形规律特征.分别选取工作面推进速度、m/d的种工况展开分析,在数值模拟的静力学分析中,为保证不同的工作面推进速度应用于数值模拟不同工况之中,故模拟将工作面推进速度等效成开挖步距来表示.为减小边界条件效应,模型沿开挖走向两端各留设 m煤柱,沿倾向各留设 m煤柱.分析工作面以、m/d的速度推进工况下,其位移分布、水压力状态、应力特征以及煤层周边破坏演化规律,为便于后续对比分析,每组工况分别选取工作面推进距离、m条件时的数据云图展开讨论与分析.不同工况下煤岩层应力变化为了研 究 不 同 工 况 开 挖 后 煤 岩 层 的 应 力 变化,模拟计算不同的开采工况,鉴于篇幅限制,以工作面推进m/d工况为例进行分析.工作面推进速度m/d工况下不同推进距离条件下煤岩层应力分布如图所示.由图可以看出,在采空区影响及承压水压共同作用下,煤层开采后煤岩体的原始应力场遭受破坏,工作面附近应力集中,并随工作面逐步推进,工作面应力集中现象越显著,最终工作面推进 m时,应力峰 值 达 到 最 大.推 进 速 度 越 快,应 力 峰 值 也越大.不同工况下煤岩层位移变化工作面推进m/d工况下不同推进距离条件下煤岩层位移分布如图所示.由图可以看出,在上覆采空区积水影响下,随着煤层开采,工作面推进过程中煤岩层的位移沉降规律经历了积水影响异常区、工作面沉降影响区以及上下煤层沉降贯通区个阶段过程.推进速度越快,煤层位移沉降量就越大.金珠鹏等:采空区积水条件下工作面推进速度对煤岩层稳定性影响研究 年第期图工作面推进m/d时的煤岩层应力分布情况图工作面推进m/d时的煤岩层位移分布情况 不同工况下煤岩层水压力变化工作面推进m/d工况下不同推进距离条件下煤岩层水压力分布如图所示.由图可以看出,随着工作面的推进,由上部煤层积水区产生的水流经过岩层裂隙涌入下部开采区,增加下部煤层突涌水的风险.对于工作面不同推进速度工况条件下,水压力随推进速度的增加而略有增大.不同工况下煤岩层岩层塑性应变工作面推进m/d工况时不同推进距离条件下煤岩层塑性应变分布如图所示.由图可以看出,塑性应变是反映地层由弹性阶段过渡至塑性损伤破坏的一个重要指标.在工作面推进初期,推进距离小于 m,由于煤层工作面远离上覆煤层积水异常区,因此仅工作面附近煤岩层产生扰动,存在塑性变形.当下覆煤层继续开采逐步靠近积水区时(推进 m),受积水异常区及下部采空区开采扰动的双重影响,其上覆煤层同样出现了塑性区.最终,工作面推进至 m后,上下煤层塑性区形成贯通,最大塑 年第期 中 国 煤 炭 第 卷性应变值达到了 .图工作面推进m/d时的煤岩层水压力分布情况图工作面推进m/d时的煤岩层塑性应变分布情况金珠鹏等:采空区积水条件下工作面推进速度对煤岩层稳定性影响研究 年第期 不同工况下层间岩层稳定性分析为分析整个煤层开采过程对煤岩层结构稳定性的影响,故选取距离初始工作面 m时煤层顶部特征点作为典型监测点,分析工作面不同的推进速度对煤层开采过程的应力、位移、水压力以及塑性应变的演化规律.测点位置如图所示.图测点位置开采过程监测点应力变化曲线如图所示.由图可以看出,在煤层开采过程中,测点应力在种工况条件下均呈现出先上升后下降的规律,当推进距离小于 m时,工作面与监测点的距离逐步靠近,受到的煤岩层扰动产生的应力不断积聚,当工作面达到监测点所在截面时(开采 m),应力集中现象最为显著.工作面继续向前推进后,由于监测点两侧均为邻空面,产生应力释放,应力值急剧减小,最后趋于稳定.同时,可以看出不同工作面推进距离工况下,应力峰值大小也存在差异,呈现出推进速度越快,应力峰值越大的规律.当工作面推进距离 m时,工作面推进速度为、m/d时应力峰值分别为 、MP a,差异性明显.开采过程监测点位移变化曲线如图 所示.由图 可以看出,在煤层开采过程中位移变化规律划分为稳定变化阶段以及骤升阶段.在工作面未达到监测点之前(推进距离小于 m),此时该测点受开采扰动的影响较小,位移增量小于 c m,当开采距离大于 m后,监测点下部采空,受上部煤层顶板、水压力以及地层荷载的共同影响,位移沉降急剧上升,直至开采结束.在煤层开采时,工作面推进速度越快,煤层测点位移沉降量就越大.当煤层开采完成后(推进距离 m),测点在推进速度、m/d时的位移沉降量分别为 、c m.可以看出工作面推进速度为m/d时位移远大于其他工况,因此煤层开采应注意工作面推进速度以确保施工安全.图开采过程监测点应力变化曲线图 开采过程监测点位移变化曲线开采过程监测点水压力变化曲线如图 所示.由图 可以看出,种工况下水压力变化规律均呈现出稳定变化阶段及骤升阶段种变化阶段,且分界点均位于工作面推进至测点所在截面(推进距离 m),此时,由于煤层采空,导致上部积水异常区涌水沿着此煤层顶板裂隙涌入采空区,导致水压骤升.最终水压力随工作面推进速度的增加而略有增大,当推进距离达到 m时,工作面推进速度为、m/d工况时水压力分别为 、k P a.开采过程监测点塑性应变变化曲线如图 所示.由图 可以看出,在工作面未推进至监测面时(推进距离 m),开采扰动并未波及监测点,因此基本未产生塑性应变.当工作面推进至此处 年第期 中 国 煤 炭 第 卷时,塑性应变急剧上升,发生塑性应变突变,之后随着煤层开采又平稳增加.此外,工作面推进速度越快,应变突变增量就越大.最终开采完毕后,该测点在工作面推进速度为、m/d时塑性应变分别为 、,可见不同的工作面推进速度对煤岩层塑性应变的影响程度显著.图 开采过程监测点水压力变化曲线图 开采过程监测点塑性应变变化曲线结论通过有限元软件A B AQU S模拟在上方有采空区积水的条件下,下方煤层开采时工作面分别以、m/d的推进速度对煤岩层结构安全稳定性的影响.分析得出在采空区积水及煤层不同开采速度共同作用影响下,煤岩层的应力、位移、水压力以及塑性应变规律.()煤层开采后煤岩体的原始应力场遭受破坏,工作面附近产生应力集中,随着工作面的推进,工作面处应力集中现象越发显著,且推进速度越快,应力峰值也越大.()工作面推进过程中煤岩层的位移沉降规律经历了积水影响异常区、工作面沉降影响区以及上下煤层沉降贯通区个阶段过程,工作面推进速度越快,煤层开采区域位移沉降量就越大.()水压力随工作面推进速度的增加而略有增大,但不同推进速度对煤层开采的水力渗流影响不大.()对于煤岩层损伤方面,随工作面与积水区距离逐渐缩短,最终诱发上部积水异常范围产生塑性区并与下部煤层贯通.工作面推进速度越快,开采扰动产生的煤岩层损伤范围就越大,塑性应变峰值同样增大.参考文献:景国勋,秦瑞琪 年我国煤矿水害事故相关因素特征分析 J安全与环境学报,():张楠慈林山矿上覆采空区积水煤层开采突水危险性评 价 J煤 炭 与 化 工,():,CHE NZ i h a n,Z OU Q u a n l e,G AO S h i k a n g,e ta l L o n g t e r m e v o l u t i o no fo v e r l y i n gr o c kf r a c t u r e si nm i n e d o u t a r e a s a n d i t s e f f e c t o ng a s f l o wc o n d u c t i v i t yJF u e l,:C U IC h u a n b o,J I A O Z h i p e n g,Z HO U Y u y i n g,e t a l E x p e r i m e n t a l s t u d yo f t h r e e d i m e n s i o n a l t i m ed o m a i no f w a t e r e v a p o r a t i o n f r o m g o a fJ A p p l i e dS c i e n c e s,():L I U W e i,Z HAN G F e n g j i e,G AO T i e g a n g,e ta l E f f i c i e n tp r e v e n t i o no fc o a ls p o n t a n e o u sc o m b u s t i o nu s i n gc o o l i n gn i t r o g e n i n j e c t i o n i na l o n g w a l l g o b:A na p p l i c a t i o nc a s eJE n e r g y,:冯国瑞,张玉江,戚庭野,等中国遗煤开采现状及 研 究 进 展 J煤 炭 学 报,():王珩珩关于未来煤炭开采技术发展方向的探讨J当代化工研究,():刘福胜基于流固耦合及复变函数分析的隧道渗流问题研究 D广州:华南理工大学,郑颖人岩土塑性力学原理:广义塑性力学 M北京:中国建筑工业出版社,(责任编辑张艳华)金珠鹏等:采空区积水条件下工作面推进速度对煤岩层稳定性影响研究 年第期