采动覆岩微震分区演化特征的数值模拟研究.pdf
《采动覆岩微震分区演化特征的数值模拟研究.pdf》由会员分享,可在线阅读,更多相关《采动覆岩微震分区演化特征的数值模拟研究.pdf(11页珍藏版)》请在咨信网上搜索。
1、采动覆岩微震分区演化特征的数值模拟研究袁国涛1,张明伟1,王杰2,卫俊1,杨坤1(1.中国矿业大学深部岩土力学与地下工程国家重点实验室,江苏徐州221000;2.中煤科工集团西安研究院(集团)有限公司,陕西西安710065)摘要:为了探究采动覆岩的破裂分区特征与微震演化规律,以某矿井工程地质条件为背景,根据矩张量与颗粒流理论构建了微震模拟方法,实现了对回采工作面覆岩微震演化特征的模拟研究。研究表明:采动覆岩的微震事件矩震级集中在2.71,微震破裂强度符合常见的正态分布性质,微震发生频次与微震矩震级的关系也满足经典的矩震级频度关系式;微震事件在采场横向上表现出明显的区域分布特征,采空区左右煤壁中
2、的微震事件先于采空区顶板发生,微震事件围绕采掘空间向两侧发展,且存在一定的超前性与滞后性;当工作面开采一定距离后,采空区上覆岩层的垂向应力表现出“双峰值”特征,微震事件迅速向上延伸,并在上覆软弱岩层内集聚、发育,使得微震事件在竖向上同样表现出明显的分层分布特征;根据采场的应力、位移以及微震分布特征,可以将上覆岩层划分为“横四区”与“竖三带”,各个破裂区域的微震事件破裂类型占比各不相同,煤壁支撑区与离层区分别为剪切型与拉伸型微震事件的主要占比区域。此外,对研究区域建立了微震监测系统,现场监测与数值模拟结果基本吻合,表明构建的微震模拟方法适用性较好,可以从细观层面再现采动覆岩的微震演化过程,研究结
3、果为揭示采动裂隙演化机理提供一定的理论基础。关键词:采动覆岩;微震;矩张量;颗粒流;破裂类型中图分类号:TD323文献标志码:A文章编号:02532336(2023)08003611Numerical simulation study on sub-regional evolution of microseismiccharacteristics of mining overburden rockYUANGuotao1,ZHANGMingwei1,WANGJie2,WEIJun1,YANGKun1(1.State Key Laboratory of Deep Geotechnical Mech
4、anics and Underground Engineering,China University of Mining&Technology,Xuzhou 221000,China;2.Xian Research Institute Co.,Ltd,China Coal Technology and Engineering Group Corp.,Xian 710065,China)Abstract:Inordertoexplorethefracturezonecharacteristicsandmicroseismicevolutionwithinoverburdenrockabovemi
5、ningcoalseam,basedonthegeologicalengineeringconditionsofamine,amicroseismicsimulationmethodwasconstructedaccordingtothemomenttensorandparticleflowtheory,andthemicroseismicevolutioncharacteristicsofoverburdenrockabovepanelsweresimulated.Theres-ultsshowthatthemagnitudeofmicroseismiceventsinmining-dist
6、urbedoverburdenisbetween2.7and1,andthemicroseismicfrac-turestrengthconformstotheGaussiandistribution.Therelationshipbetweenmicroseismicfrequencyandmomentmagnitudealsosatis-fiestheclassicalmomentmagnitude-frequencyrelationship.Themicroseismiceventsshowprominentregionaldistributioncharacteristicsinthe
7、horizontaldirectionofthepanel.Themicroseismiceventsintheleftandrightcoalwallsofthegoafoccurearlierthantheroofabovegoaf,andthemicroseismiceventsdeveloptobothsidesaroundtheminingspace,andthereisacertainadvanceandlag.Whenthepanelisminedforacertaindistance,theverticalstressoftheoverburdenabovethegoafsho
8、wsadoublepeakcharacteristic,andthemicroseismiceventsextendupwardrapidlyandgatheranddevelopintheweakoverburdensothatthemicroseismiceventsalsoshowprominentlayereddistributioncharacteristicsintheverticaldirection.Accordingtothestress,displacement,andmicroseismicdistributioncharacteristics收稿日期:20220509责
9、任编辑:常琛DOI:10.13199/ki.cst.2022-0693基金项目:国家自然科学基金资助项目(52074260);中央高校基本科研业务费专项资金资助项目(2019QNA26)作者简介:袁国涛(1993),男,山东菏泽人,博士研究生。E-mail:通讯作者:张明伟(1984),男,山东潍坊人,副研究员,博士。E-mail:第51卷第8期煤炭科学技术Vol.51No.82023年8月CoalScienceandTechnologyAug.2023袁国涛,张明伟,王杰,等.采动覆岩微震分区演化特征的数值模拟研究J.煤炭科学技术,2023,51(8):3646.YUANGuotao,ZHA
10、NGMingwei,WANGJie,et al.Numericalsimulationstudyonsub-regionalevolutionofmi-croseismiccharacteristicsofminingoverburdenrockJ.CoalScienceandTechnology,2023,51(8):3646.36withinoverburden,theoverburdencanbedividedintofourhorizontalzonesandthreeverticalzones.Theproportionofmicroseismiceventsandfracturet
11、ypesineachfractureareaisdifferent.Thecoalwallsupportedandbedseparationareasarethelargeproportionareasofshearandtensilemicroseismicevents,respectively.Inaddition,amicroseismicmonitoringsystemwasestablishedforthestudyarea.Thefieldmonitoringandnumericalsimulationresultsareconsistent,indicatingthattheco
12、nstructedmicroseismicsimulationmethodhasagoodapplicabilityandcanreproducethemicroseismicevolutionprocessofmining-disturbedoverburdenfromthemesoscopiclevel.Theresearchresultsprovideatheoreticalbasisforrevealingtheevolutionmechanismofmining-inducedfractures.Key words:miningoverburdenrock;microseism;mo
13、menttensor;particleflow;fracturetype0引言采动覆岩裂隙延展规律研究作为工程界十分关注的课题,与采掘过程中发生的各类工程灾害、环境灾害等问题密切相关12。因此,有效掌握开采煤层引起的覆岩破坏特征及裂隙演化规律,对煤矿的安全生产和绿色开采具有十分重要的意义。微震监测技术起源于地球物理学领域3,因具有大区域性、高连续性等特点,近几年被广泛应用于岩体破裂过程的实时监测中,并取得了诸多有益的研究成果。CHENG 等4基于矿山地质资料与微震监测技术,研究了采动覆岩横向和垂向的微震事件分布规律,并在此基础上构建了覆岩在横向和垂向的破裂分区方法。侯恩科等5利用井地联合微震监
14、测技术,从微震事件的能量和密度角度探讨深埋煤层开采顶板覆岩破坏过程及导水裂隙带发育高度。孟祥军等6运用微震监测技术对工作面推进过程中微震事件进行实时记录,并分析了采动覆岩的裂隙发育形态特征和演化趋势,为高位钻场瓦斯抽采参数设计与优化提供一定的依据。综上可知,前人关于采动覆岩裂隙已进行了大量微震监测现场实践工作,并取得了丰富的研究成果。但是对于采动覆岩微震演化特征背后的细观机制研究较少,加上微震监测技术在实际工况应用中,存在定位误差偏大7、源能计算方法不固定8等特点,因此急需寻求另一种研究方法对此加以完善与补充,从细观层面上揭示采动覆岩微震演化特征背后的力学机制。颗粒流离散元计算方法可以从细观层
15、面再现岩体破裂过程,已被广泛应用于各类围岩力学特性与工程的模拟910,此外,有学者将该方法与矩张量理论结合,成功实现了室内尺度的岩体破裂过程的声发射模拟1113。鉴于此,基于颗粒流离散元计算方法建立了采动覆岩的数值模型,并在此基础上采用矩张量理论的微震计算方法,对不同破裂分区条件下的微震事件演化规律进行分析,深入探讨采动覆岩的细观破裂过程,为从微震监测角度揭示采动裂隙演化机理提供一定的理论基础。1采动覆岩的破裂分区模型l1123l2大量的研究表明,煤层开采后,上覆岩体的应力平衡被打破,应力重新分布。图 1 为采场上覆岩层的破裂分区模型考虑到煤壁支撑区的影响,在水平方向可以将顶板岩层划分为 4
16、个区,从左到右分别为原岩区(A 区)、煤壁支撑区(B 区)、离层区(C 区)和再压实区(D 区)4,14。4 个区域的应力状态不同,其岩体运移规律以及破裂特征也有所不同。原岩区位于采空区前方,超前距离较远,处于应力平衡状态,在回采期间一般不会产生新的裂隙。煤壁支撑区紧靠回采工作面前方,受超前应力的影响,压力峰值出现在工作面前方一定距离处,工作面到支撑压力峰值位置的煤岩体处于极限平衡状态,主要由剪应力控制,煤壁支撑区范围、倾角与顶板岩层刚度及所受超前应力等因素有关。离层区主要由拉应力主导,该区域岩层在重力应力作用下发生弯曲变形,由于煤岩体力学性质的差异,所表现的弯曲变形程度不同,从而导致层间裂隙
17、的形成15。此外,在拉应力作用下岩层会产生断裂破坏,并不断向上发展,随着原生裂隙与次生裂隙的不断延展,岩层之间也产生了复杂的剪应力,进而加速了裂隙和裂缝的发展,从而形成一定角度的裂隙并向上扩展。离层区裂隙的发展角、以及范围 都与上覆岩体的刚度有关。再压实区为应力再平衡区,采空区顶板岩体经过垮落、裂隙延展等过程,在自重应力的作用下逐渐达到区域平衡状态,由于该部分区域是通过前两个区域过渡而来,因而此区域岩体破坏最为严重,整体弯曲下沉量最大。A 区B 区123l1l2C 区D 区IIIIII图1采场上覆岩层的破裂分区模型Fig.1Fracturezoningmodelofstrataonmining
18、site袁国涛等:采动覆岩微震分区演化特征的数值模拟研究2023年第8期37在垂直方向上,根据上覆岩体破断类型和裂隙发育情况,从下到上依次分为垮落带 I、裂缝带 II 和弯曲下沉带 III14,如图 1 所示。垮落带位于煤层上方,由杂乱无序的破碎煤岩体组成,纵向破断裂隙与横向离层裂隙在此区域纵横交错分布。裂缝带位于垮落带上方,其形成过程伴随着岩体的破断运移以及应力重分布,该区域存在 2 种类型的裂隙,第 1 类为采空区上覆岩层周期性破断下沉运动产生的破断裂隙,主要由拉应力控制,裂隙间的水平间距与岩层的破断步距近似相等;第 2 类为煤壁支撑区产生的岩层压剪裂隙,主要由剪应力控制,这类裂隙的分布相
19、对杂乱无序,裂缝尺度通常小于前者,分布密度却要高于前者1617。弯曲下沉带位于裂缝带上方的部分岩层与松散层区域,该区域岩体不发生破坏,只产生一定程度的弯曲变形。此外,有学者认为,关键层的存在对上覆岩层的变形和断裂过程起到决定性作用18,在关键层不发生断裂的前提下,关键层可以有效防止下方岩体破坏与运移传递到关键层之上的地层,最终导致了“竖三带”区域边界通常与关键地层一致。综上所述,根据采场上覆岩层的应力、移动破断以及裂隙分布特征,可以简单的将上覆岩层划分为“横四区”与“竖三带”,通过深入研究各个分区的采动裂隙形成过程及类型,为采动覆岩破裂分区治理奠定坚实的基础。2数值模型建立与计算2.1工程概况
20、f以兖州矿区某矿 325 工作面为研究背景。325工作面埋深约 270m,走向长度约 1100m,倾向宽度约 100m,煤层平均厚度 8.2m,煤层结构简单,赋存稳定,采用长壁式综采放顶煤开采。该工作面为单斜构造,煤层倾角平均约 4,总体呈现南高北低的趋势。煤层直接顶为粉砂岩,坚固性系数 为 3.06.0,属于中硬覆岩,工作面钻孔柱状图如图 2 所示。2.2PFC 与微震离散元模拟实现PFC 属于 Itasca 公司开发的力学分析软件,随着几十年的不断发展完善,已成为目前岩土、采矿领域相关问题数值计算最为有效的研究工具1011。该软件主要是通过牛顿第二定律研究颗粒(簇)以及接触之间的力与位移的
21、关系,从细观的角度来分析材料的宏观力学特性,进而实现煤岩材料裂纹萌生与扩展过程的模拟。然而,由于 PFC 中的颗粒尺寸几乎相同,因此颗粒之间的黏结键断裂所释放的能量几乎相同,这与真实的微震规律有所不同。为了更加可靠的描述微震破裂强度分布规律,一般将发生在特定时间和空间上相近的黏结断裂视为同一微震事件,当微震事件由多条微裂纹组成时,则将裂纹群的几何中心作为微震事件的质心11。在岩石力学以及采矿领域,矩张量理论已成为研究岩体内部破裂震源机制及时空演化规律的主流工具1920。一般来说,任何震源的应力状态都可以通过矩张量进行描述,在地震学中,矩张量可以通过记录震源释放出的动力波转化获得20;在 PFC
22、 软件中,微震事件所属颗粒的位置、力和位移都是已知的,那么根据黏结破坏时周围颗粒接触力的变化很容易计算出微震事件的矩张量11,13。Mij当黏结键在外部负载下发生断裂时,每个黏结键断裂都代表一个微裂纹。微裂纹两端原先接触的颗粒为源颗粒,随着微裂纹的产生,裂纹周围的源颗粒将发生移动,进而引起接触力发生变化。将源颗粒上所有接触的接触力变化量乘以对应的力臂,并进行求和运算,得到表征微震事件的矩张量分量11,13,其表达式为Mij=SFij(1)FiijjS式中,为接触力变化量的第 个分量;为接触点与微裂纹中心距离的第 个分量;为破裂面,图 3 为PFC 中微震事件模拟的实现过程示意。岩性柱状层厚/m
23、底部距煤层高/m第四系松散层183.1985.84粗砂岩1.5184.33黏土岩12.0172.32中砂岩4.3268.00黏土岩2.1365.87中细砂岩互层17.4748.40黏土岩2.3746.03粉砂岩11.6934.34黏土岩2.4531.89粉砂岩4.0327.862 号煤0.6527.21黏土岩1.8525.36中细砂岩互层14.3810.98粉砂岩10.9803 号煤8.14泥岩1.191.19粉细砂岩互层12.0813.27图2325 工作面钻孔柱状Fig.2DrillinghistogramnearNo.325workingface2023年第8期煤炭科学技术第51卷38微
24、震事件所属源颗粒组成微震事件的裂纹群微震事件的矩张量微震事件所属源颗粒的接触力图3PFC 中微震事件模拟的实现过程Fig.3RealizationprocessofmicroseismiceventsimulationinPFCM0微震现场试验和监测结果表明,微震事件释放的能量随时间不断演化。那么对应的矩张量也可以表示为时间的函数2021。在 PFC 模拟方法中,可以通过计算每一时步下的矩张量,获得具有全程时间相关性的矩张量。然而,在实际的运算中,若存储所有微震事件的全程时间相关性的矩张量,将耗费大量的内存及存储资源。为了提高计算效率,采用具有最大标量力矩值时的矩张量作为每一微震事件的矩张量并
25、存储。根据矩张量矩阵,标量力矩的表达式11为:M0=(3j=1m2j2)1/2(2)mjjM式中:为矩张量矩阵的第 个特征值。根据微震事件矩张量的最大标量力矩,可由式(3)计算微震事件的矩震级:M=23lgM06(3)M矩震级反映了微震事件的破裂强度。值得说明的是,上述矩震级求解公式借鉴了地震学中矩震级标度的定义方法,它是一个均匀震级模型,适合于震级尺度范围很宽的统计,目前被广泛应用于地震以及微震事件的破裂强度表征11,13,21。此外,为了量化矩张量的破裂类型,需对微震矩张量做一步分解工作,进而建立矩张量参量和岩体破裂类型之间的关系。本文采用 FEIGNIER21提出的破裂类型判据,该判据对
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 采动覆岩微震 分区 演化 特征 数值 模拟 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。