煤矿深部开采动力灾害区域应力效应监测与前兆特征研究.pdf
《煤矿深部开采动力灾害区域应力效应监测与前兆特征研究.pdf》由会员分享,可在线阅读,更多相关《煤矿深部开采动力灾害区域应力效应监测与前兆特征研究.pdf(10页珍藏版)》请在咨信网上搜索。
1、动力灾害的孕育及诱发与区域应力场的力学行为密切相关,煤矿开采扰动作用下区域局部应力场不断调整变化。为了研究区域应力场变化特征对动力灾害的前兆响应,采用变分模态分解和希尔伯特变换的方法对鲍店煤矿采集数据进行时频分析,识别和提取钻孔应变数据中反映动力异常变化的固有模基金项目:山东省科技重大创新项目(2019SDZY02,2019SDYZ05)This research is financially supported by the Major Scientific and Technological Innovation Projects of Shandong Province(Grants 2
2、019SDZY02and2019SDYZ05).第一作者:吴浩源(1988),男,博士,主要从事矿山深部开采动力灾害等方面的研究工作。E-mail:通讯作者:孙源源(1987),男,工程师,主要从事矿山建设工作。E-mail:收稿日期:20230220;修回日期:20230505;责任编辑:王婧第29卷第3期地质力学学报Vol.29No.32023年6月JOURNALOFGEOMECHANICSJun.2023态分量。研究结果表明,钻孔应变观测数据能够有效记录开采扰动相关区域岩体内部的微小变化,解析原始数据滤去干扰信号后获得表征动力异常的特定频段,分解变换后固有模态分量异常波动特征早于动压事件
3、前 27 天出现,表现出形变稳定变化阶段向快速变化阶段“突跳”、形变快速变化阶段“震荡”以及失稳阶段“跌落”的显著特征。依据地震形变三阶段理论,判别动力灾害发生前形变异常开始时刻和“突跳”“震荡”“跌落”作为煤矿动力灾害发生的前兆判据。基于钻孔应变仪对区域应力场的精确观测,构建了应用于煤矿动力灾害超前预警的判据方法,可为煤矿相似开采环境下工作面的安全高效回采提供借鉴和参考。关键词:动力灾害;变分模态分解;钻孔应变;希尔伯特变换;前兆特征中图分类号:TD32;P315.72+7文献标识码:A文章编号:10066616(2023)03035510DOI:10.12090/j.issn.1006-6
4、616.202329050引言近年来,随着煤矿开采深度不断加深,动力灾害的发生日趋频繁。其中,冲击地压灾害的破坏力巨大,造成的后果严重影响了煤矿安全开采。从工作面的尺度出发,煤层在开采过程中常伴随上覆坚硬岩层的突然破裂,引发岩体强烈的冲击破坏,表征为一定范围内的破坏现象(纪洪广等,2003)。近年来,专家学者着眼于工作面尺度之外更大的范围去研究煤矿相关动力灾害的发生机制及预警,将研究范围扩展至区域应力场的范畴(李铁等,2003,刘金海等,2014;张月征等,2016a)。甘肃华亭煤矿于2014 年发生的一次冲击地压灾害,其释放的能量达到了 108J 量级,造成1500m 巷道底鼓 1.5m,与
5、工作面尺度内的顶板结构失稳相比,该灾害表现出破坏力强、波及范围广等特点,其根本原因在于周边区域内煤岩体为冲击地压灾害的发生提供了必要的力能条件(张月征等,2016b)。因此,在工作面开采期间,应从更广的区域范围入手,寻求不同尺度的动力异常响应特征,以期探索深部煤矿开采过程中动力灾害的发生及预警机制。天然地震和矿震都是地层内部的动力学现象,可将天然地震超前预报的观测手段应用到矿震或冲击地压灾害的监测中,其中地壳形变观测是地震地面监测体系中重要的前兆观测项目(邱泽华和石耀霖,2004),监测台站观测到地应力拉伸脉冲异常与地壳地震活动显著相关(邱泽华,2017)。采用S 变换和高通滤波从钻孔应变数据
6、中提取短周期信号,能够发现汶川地震前钻孔应变观测出现异常(刘琦和张晶,2011)。借助钻孔应变观测仪器实现对地壳变形的精确记录,可及时捕捉地壳内部如固体 潮 汐 量 级 的 微 小 变 化(苏 恺 之,2003;邱 泽 华,2010;张聪聪,2015)。钻孔应变仪因其频带宽、精度高,适用于监测煤矿深部开采动力灾害异常响应特征。为了揭示煤层开采过程中区域应力场的演变规律,需要考虑开采扰动与区域应力场协同作用对矿震活动性的相互影响,往往需要捕捉长短周期信号(云龙等,2021)。相比监测煤矿动力灾害的其他设备,钻孔应变仪因其监测频带宽的特性(池顺良,2013),所采集的数据含有高频短周期的噪声和低频
7、长周期的应变固体潮频段以及能够反映矿震活动规律的特定频段,为探究工作面采动环境中矿震孕育与诱发机制,需滤去噪声和应变固体潮频段。针对时域内难以有效区分煤矿动力灾害异常响应信息的问题,文章以鲍店煤矿七采区东翼不规则区域(73上11 工作面)开采期间观测的钻孔应变数据为基础,采用变分模态分解(VMD)方法,获得不同频率的固有模态函数(IMF)分量。通过构建带通滤波通道,滤去高斯噪声信号和固体潮分量,重构反映矿震活动的 IMF 分量,利用希尔伯特(Hilbert)变换将时域内看似杂乱无章却能反映矿压活动的应变信息转换成频域信息,求解该分量的瞬时能量和以天为单位的能量累计变化,分析信号异常处矿震的发生
8、与区域应力场之间的相互作用效应,实现对矿震或冲击地压灾害前兆信息的识别和解析,以期对相似地质条件下煤矿动力灾害的预测预报提供一定指导意义。1工程概况及资料选取鲍店煤矿 73上11 工作面位于井田东北部边界,东部与东滩煤矿四采区(43上09 采空区)、十四采区(14316 采空区和 14302 采空区)相接,北部与兴隆庄356地质力学学报https:/2023煤矿四采区相接,相邻矿井的 3 个采区回采工作已完成,西部为鲍店煤矿 73上02 采空区和 73上04 采空区,南部为十采区、七采区(73上09 采空区)南翼采区,相邻采区已基本采完(图 1)。因此,工作面之间相互关系呈现孤岛工作面开采。其
9、中,七采区东翼南侧已回采,待回采的 73上11 工作面处于不规则三角形区域(图 1)。73上11 工作面主采 3#煤层,煤层采深超过 600m,且距煤层上方 100m 处有厚硬砂岩层,其破断后释放大量弹性应变能,该厚硬岩层的失稳对工作面安全开采影响显著。根据地质勘探结果,已探明鲍店煤矿七采区厚硬砂岩最小厚度为41.75m,最大厚度超过 373.55m。由采区东北部边界向东西翼分界线逐渐增加,73上11 工作面上覆厚硬砂岩层厚度大于 300m,部分区域超过 350m,工作面周边采空区覆岩运动以及应力场调整势必增加 73上11 工作面应力分布的复杂性。综上所述,73上11 工作面开采期间主要受孤岛
10、工作面导致的应力集中和上覆厚硬岩层稳定性的影响,故 73上11 工作面在掘进和回采过程中动力异常响应较为突出,易引发矿震或冲击地压相关动力灾害。14316采空区14302采空区43上09 采空区73上11 工作面73上02 采空区73上04 采空区73上09采空区图1鲍店煤矿 73上11 工作面邻近关系示意图(海拔546.8m)Fig.1Neighboring relationship diagram of the working faceposition(Thealtitudeis546.8m)为了揭示矿震发生前的异常响应特征,选取能量为 105J 的矿震事件作为研究对象,截取 5 月份和1
11、0 月份每次矿震发生前 20 天和震后 10 天作为一个完整的时间序列,基于 VMD 识别和提取钻孔应变数据。2基于 VMD 识别和提取钻孔应变数据工作面开采期间,钻孔应变仪因其具有高分辨率、高精度等特性,采集的数据会受到所处环境中干扰因素的影响,如高斯噪声的即时干扰和应变固体潮的周期性干扰等(文勇等,2012)。消除噪声常见的滤波方法主要有小波变换、傅里叶变换或短时变换(穆峰等,2015),上述方法在处理平稳信号时所表现的效果是满足滤波需求的,但是在处理非平稳信号时误差较大,因非平稳信号的特性是其信号自身频率随时间变化而变化,傅里叶变换在处理相关问题时只能给出总体效果,并不能够反应出某一时
12、刻 信 号 的 局 部 信 息。在 此 基 础 上,Huangetal.(1998)提出经验模态分解(EMD)来处理非平稳、非线性信号,其特点是将原始信号按照频率从高到低的顺序分解为若干个固有模态函数分量,分解后的各模态函数具备其自身固有的属性;然后,再根据希尔伯特变换将信号从时域角度转换至频域角度,这在面对时域内非线性、非平稳信号分析效果差 的 环 境 下,提 供 了 有 力 的 技 术 支 撑。但 随 着EMD 分解方法在信号处理方面的广泛应用,其缺点也暴露出来,如其端点效应和模态混叠严重,不具备严格的数学模型。基于此问题,在 EMD 基础上又发展出很多改进算法,如集合经验模态分解(EEM
13、D;HuangandWu,2008;陈可等,2010)、补充集合经验模态分解(CEEMD;Yehetal.,2010)和改进集合经验模 态 分 解(MEEMD;郑 近 德 等,2013)等。EEMD是通过向原始信号中加入高斯白噪声来抑制经验模态分解过程中存在的模态混叠现象,CEEMD 是通过向目标信号添加成对的符号相反的白噪声,来确保信号分解具有真实的物理意义,MEEMD 结合CEEMD 与排列熵(PE)算法在抑制模态混叠方面取得理想的结果(邢昀和荣剑,2018)。上述相关改进算法虽然在一定程度上解决了模态混叠和端点效应问题(康佳星等,2016),但本质上仍是时域范畴内 的 递 归 求 解。随
14、 后,Dragomiretskiy and Zosso(2014)提出的一种自适应信号处理方法(VMD),通过迭代不断更新各模态函数及中心频率,得到若干具有一定带宽的模态函数,获取变分模态的最优解。VMD 相较于 EMD 是基于频域的非递归求解算法,弥补了 EMD 的缺点。井下工作环境、地质环境复杂多变(谭元隆等,2021),钻孔应变仪采集的数据记录了不同来源、不同特性的混合时频信号,仅仅从应变数据波形特征第3期吴浩源,等:煤矿深部开采动力灾害区域应力效应监测与前兆特征研究357的时变角度难以精确分析局部细节信息。但诸如噪声、应变固体潮等某个单频信号蕴含着固有的优势频率,采用 VMD 的方法,
15、能够将具有不同来源、不同频率的波动信号分离出来,进而针对不同的研究问题,选择敏感频率段的数据作为研究对象。文章采用四分量钻孔应变仪分析采动环境中应变的相对变化规律,四分量钻孔应变仪通过记录元件电容值实现对应变的观测。为界定观测数据质量是否满足研究问题的要求,需要依靠质量评价体系评定数据的可靠性,而钻孔应变数据质量评价体系则是自洽方程(式)。四分量钻孔应变仪采集的数据经过应变换算后得到面应变 Sa(邱泽华等,2009),依靠自洽方程既可以评判数据自身的可靠性,也可以用来识别四分量钻孔应变仪的工作状态。若数据段内整体变化趋势一致,则满足自洽方程:S1+S3=S2+S4S1S2S3S4式中:元件 1
16、 的观测数值;元件 2 的观测数值;元件3 的观测数值;元件4 的观测数值。动力灾害的孕育及其诱发与区域应力场的力学行为密切相关,这里的动力灾害主要涉及冲击地压、矿震以及工作面来压期间发生的顶板异常动力显现事件。针对能级跃迁较大的煤矿动力灾害,钻孔应变曲线走势表现为显著的“阶跃”(张月征等,2016b),即数据曲线在动力灾害发生时刻产生台阶式上升或下降。此外,还存在一些能级跃迁较小的动压事件。5 月 4 日5 月 31 日四分量钻孔应变仪分量 1 和分量 3 根据自洽方程拟合后的观测数据(图 2a),9 月 22 日10 月 13 日分量 2 和分量 4 根据自洽方程拟合后的观测数据(图 2b
17、),均满足自洽方程,图中仅绘制拟合后的观测曲线。观察两个不同时间序列上钻孔应变数据的幅值变化,可知某一时刻发生动力异常事件的典型状态失稳现象在应变数据波动特征上并不能够有效直观地显现,应力传递路径受钻孔应变仪安装位置、震源处煤岩介质力学属性等影响,动力异常事件发生前后微小的响应信号大都淹没在固体潮背景下,故有必要根据时频变换识别或提取区域应力场中细微变化与动力灾害的前兆响应特征。1.951.901.851.801.751.701.65应变/1071.931.921.911.901.89应变/10700.510411041.510421042.510431043.5104410400.51041
18、1041.510421042.51043104时间/min时间/min(a)5 月 4 日5 月 31 日应变观测曲线(b)9 月 22 日10 月 13 日应变观测曲线设备断电造成异常5 月 21 日事件9 月 24 日事件10 月 4 日事件10 月 11 日事件5 月 26 日事件5 月 29 日事件图2时频解析前钻孔应变仪采集原始数据曲线Fig.2Rawdatacurvesbytheboreholestrainmeterbeforetime-frequencyanalysis钻孔应变仪采集的是多频段混合后的数据(欧阳祖熙等,2009),经 VMD 分解后得到若干个固有模态函数(IMF)
19、。定义有意义的瞬时频率其必要条件是函数关于零均值局部对称且过零点的数目与极点的数目一致。据此,Huangetal.(1998)定义了IMF 应满足以下 2 个条件:在整个数据区间内,极值点的数目与过零点的数目相等或至多相差 1 个;在任意一点处,由局部极大值点定义的包络以及由局部极小值点定义的包络的均值为零。此外,为了消除噪声和固体潮等因素对要研究数据段的影响,在 VMD 分解的基础上,采用希尔伯特(Hilbert)变换依次分析 IMF 分量,将时域信息转换至频域信息范畴。按照上述理论对 5 月 1 日5 月 31 日的数据段进行分解,获得各频段 IMF 曲线(图 3)。如表 1 所示,IMF
20、1IMF5 周期皆小于 1h,推测可能是设备自身或周围环境中产生的高斯噪声(Barmanetal.,2016),故剔除 IMF1IMF5 分量,以排除噪声对原始数据的干扰。IMF9 和 IMF10 的频率分别为 2.24105Hz 与 1.16105Hz,转换至周期后分别为 12.4h 和 23.9h,时间维度上与固体潮严格匹配,确定为固体潮的半日波和日波分量。此外,与 IMF9 和 IMF10 邻近的 IMF11 和 IMF12 为固体潮谐波分量。为探究噪声频段对所研频段数据造成的干扰,综上所述,将混杂的各频段数据解析后重构,对比分析降噪前与降噪后的频段曲线(图 4)。降噪信号358地质力学
- 配套讲稿:
如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。