煤矿深部开采动力灾害区域应力效应监测与前兆特征研究.pdf

动力灾害的孕育及诱发与区域应力场的力学行为密切相关，煤矿开采扰动作用下区域局部应力场不断调整变化。为了研究区域应力场变化特征对动力灾害的前兆响应，采用变分模态分解和希尔伯特变换的方法对鲍店煤矿采集数据进行时频分析，识别和提取钻孔应变数据中反映动力异常变化的固有模基金项目：山东省科技重大创新项目（2019SDZY02，2019SDYZ05）This research is financially supported by the Major Scientific and Technological Innovation Projects of Shandong Province(Grants 2019SDZY02and2019SDYZ05).第一作者：吴浩源（1988），男，博士，主要从事矿山深部开采动力灾害等方面的研究工作。E-mail:通讯作者：孙源源（1987），男，工程师，主要从事矿山建设工作。E-mail:收稿日期：20230220；修回日期：20230505；责任编辑：王婧第29卷第3期地质力学学报Vol.29No.32023年6月JOURNALOFGEOMECHANICSJun.2023态分量。研究结果表明，钻孔应变观测数据能够有效记录开采扰动相关区域岩体内部的微小变化，解析原始数据滤去干扰信号后获得表征动力异常的特定频段，分解变换后固有模态分量异常波动特征早于动压事件前 27 天出现，表现出形变稳定变化阶段向快速变化阶段“突跳”、形变快速变化阶段“震荡”以及失稳阶段“跌落”的显著特征。依据地震形变三阶段理论，判别动力灾害发生前形变异常开始时刻和“突跳”“震荡”“跌落”作为煤矿动力灾害发生的前兆判据。基于钻孔应变仪对区域应力场的精确观测，构建了应用于煤矿动力灾害超前预警的判据方法，可为煤矿相似开采环境下工作面的安全高效回采提供借鉴和参考。关键词：动力灾害；变分模态分解；钻孔应变；希尔伯特变换；前兆特征中图分类号：TD32；P315.72+7文献标识码：A文章编号：10066616（2023）03035510DOI：10.12090/j.issn.1006-6616.202329050引言近年来，随着煤矿开采深度不断加深，动力灾害的发生日趋频繁。其中，冲击地压灾害的破坏力巨大，造成的后果严重影响了煤矿安全开采。从工作面的尺度出发，煤层在开采过程中常伴随上覆坚硬岩层的突然破裂，引发岩体强烈的冲击破坏，表征为一定范围内的破坏现象（纪洪广等，2003）。近年来，专家学者着眼于工作面尺度之外更大的范围去研究煤矿相关动力灾害的发生机制及预警，将研究范围扩展至区域应力场的范畴（李铁等，2003，刘金海等，2014；张月征等，2016a）。甘肃华亭煤矿于2014 年发生的一次冲击地压灾害，其释放的能量达到了 108J 量级，造成1500m 巷道底鼓 1.5m，与工作面尺度内的顶板结构失稳相比，该灾害表现出破坏力强、波及范围广等特点，其根本原因在于周边区域内煤岩体为冲击地压灾害的发生提供了必要的力能条件（张月征等，2016b）。因此，在工作面开采期间，应从更广的区域范围入手，寻求不同尺度的动力异常响应特征，以期探索深部煤矿开采过程中动力灾害的发生及预警机制。天然地震和矿震都是地层内部的动力学现象，可将天然地震超前预报的观测手段应用到矿震或冲击地压灾害的监测中，其中地壳形变观测是地震地面监测体系中重要的前兆观测项目（邱泽华和石耀霖，2004），监测台站观测到地应力拉伸脉冲异常与地壳地震活动显著相关（邱泽华，2017）。采用S 变换和高通滤波从钻孔应变数据中提取短周期信号，能够发现汶川地震前钻孔应变观测出现异常（刘琦和张晶，2011）。借助钻孔应变观测仪器实现对地壳变形的精确记录，可及时捕捉地壳内部如固体 潮 汐 量 级 的 微 小 变 化（苏 恺 之，2003；邱 泽 华，2010；张聪聪，2015）。钻孔应变仪因其频带宽、精度高，适用于监测煤矿深部开采动力灾害异常响应特征。为了揭示煤层开采过程中区域应力场的演变规律，需要考虑开采扰动与区域应力场协同作用对矿震活动性的相互影响，往往需要捕捉长短周期信号（云龙等，2021）。相比监测煤矿动力灾害的其他设备，钻孔应变仪因其监测频带宽的特性（池顺良，2013），所采集的数据含有高频短周期的噪声和低频长周期的应变固体潮频段以及能够反映矿震活动规律的特定频段，为探究工作面采动环境中矿震孕育与诱发机制，需滤去噪声和应变固体潮频段。针对时域内难以有效区分煤矿动力灾害异常响应信息的问题，文章以鲍店煤矿七采区东翼不规则区域（73上11 工作面）开采期间观测的钻孔应变数据为基础，采用变分模态分解（VMD）方法，获得不同频率的固有模态函数（IMF）分量。通过构建带通滤波通道，滤去高斯噪声信号和固体潮分量，重构反映矿震活动的 IMF 分量，利用希尔伯特（Hilbert）变换将时域内看似杂乱无章却能反映矿压活动的应变信息转换成频域信息，求解该分量的瞬时能量和以天为单位的能量累计变化，分析信号异常处矿震的发生与区域应力场之间的相互作用效应，实现对矿震或冲击地压灾害前兆信息的识别和解析，以期对相似地质条件下煤矿动力灾害的预测预报提供一定指导意义。1工程概况及资料选取鲍店煤矿 73上11 工作面位于井田东北部边界，东部与东滩煤矿四采区（43上09 采空区）、十四采区（14316 采空区和 14302 采空区）相接，北部与兴隆庄356地质力学学报https:/2023煤矿四采区相接，相邻矿井的 3 个采区回采工作已完成，西部为鲍店煤矿 73上02 采空区和 73上04 采空区，南部为十采区、七采区（73上09 采空区）南翼采区，相邻采区已基本采完（图 1）。因此，工作面之间相互关系呈现孤岛工作面开采。其中，七采区东翼南侧已回采，待回采的 73上11 工作面处于不规则三角形区域（图 1）。73上11 工作面主采 3#煤层，煤层采深超过 600m，且距煤层上方 100m 处有厚硬砂岩层，其破断后释放大量弹性应变能，该厚硬岩层的失稳对工作面安全开采影响显著。根据地质勘探结果，已探明鲍店煤矿七采区厚硬砂岩最小厚度为41.75m，最大厚度超过 373.55m。由采区东北部边界向东西翼分界线逐渐增加，73上11 工作面上覆厚硬砂岩层厚度大于 300m，部分区域超过 350m，工作面周边采空区覆岩运动以及应力场调整势必增加 73上11 工作面应力分布的复杂性。综上所述，73上11 工作面开采期间主要受孤岛工作面导致的应力集中和上覆厚硬岩层稳定性的影响，故 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 月份和10 月份每次矿震发生前 20 天和震后 10 天作为一个完整的时间序列，基于 VMD 识别和提取钻孔应变数据。2基于 VMD 识别和提取钻孔应变数据工作面开采期间，钻孔应变仪因其具有高分辨率、高精度等特性，采集的数据会受到所处环境中干扰因素的影响，如高斯噪声的即时干扰和应变固体潮的周期性干扰等（文勇等，2012）。消除噪声常见的滤波方法主要有小波变换、傅里叶变换或短时变换（穆峰等，2015），上述方法在处理平稳信号时所表现的效果是满足滤波需求的，但是在处理非平稳信号时误差较大，因非平稳信号的特性是其信号自身频率随时间变化而变化，傅里叶变换在处理相关问题时只能给出总体效果，并不能够反应出某一时 刻 信 号 的 局 部 信 息。在 此 基 础 上，Huangetal.（1998）提出经验模态分解（EMD）来处理非平稳、非线性信号，其特点是将原始信号按照频率从高到低的顺序分解为若干个固有模态函数分量，分解后的各模态函数具备其自身固有的属性；然后，再根据希尔伯特变换将信号从时域角度转换至频域角度，这在面对时域内非线性、非平稳信号分析效果差 的 环 境 下，提 供 了 有 力 的 技 术 支 撑。但 随 着EMD 分解方法在信号处理方面的广泛应用，其缺点也暴露出来，如其端点效应和模态混叠严重，不具备严格的数学模型。基于此问题，在 EMD 基础上又发展出很多改进算法，如集合经验模态分解（EEMD；HuangandWu，2008；陈可等，2010）、补充集合经验模态分解（CEEMD；Yehetal.，2010）和改进集合经验模 态 分 解（MEEMD；郑 近 德 等，2013）等。EEMD是通过向原始信号中加入高斯白噪声来抑制经验模态分解过程中存在的模态混叠现象，CEEMD 是通过向目标信号添加成对的符号相反的白噪声，来确保信号分解具有真实的物理意义，MEEMD 结合CEEMD 与排列熵（PE）算法在抑制模态混叠方面取得理想的结果（邢昀和荣剑，2018）。上述相关改进算法虽然在一定程度上解决了模态混叠和端点效应问题（康佳星等，2016），但本质上仍是时域范畴内 的 递 归 求 解。随 后，Dragomiretskiy and Zosso（2014）提出的一种自适应信号处理方法（VMD），通过迭代不断更新各模态函数及中心频率，得到若干具有一定带宽的模态函数，获取变分模态的最优解。VMD 相较于 EMD 是基于频域的非递归求解算法，弥补了 EMD 的缺点。井下工作环境、地质环境复杂多变（谭元隆等，2021），钻孔应变仪采集的数据记录了不同来源、不同特性的混合时频信号，仅仅从应变数据波形特征第3期吴浩源，等：煤矿深部开采动力灾害区域应力效应监测与前兆特征研究357的时变角度难以精确分析局部细节信息。但诸如噪声、应变固体潮等某个单频信号蕴含着固有的优势频率，采用 VMD 的方法，能够将具有不同来源、不同频率的波动信号分离出来，进而针对不同的研究问题，选择敏感频率段的数据作为研究对象。文章采用四分量钻孔应变仪分析采动环境中应变的相对变化规律，四分量钻孔应变仪通过记录元件电容值实现对应变的观测。为界定观测数据质量是否满足研究问题的要求，需要依靠质量评价体系评定数据的可靠性，而钻孔应变数据质量评价体系则是自洽方程（式）。四分量钻孔应变仪采集的数据经过应变换算后得到面应变 Sa（邱泽华等，2009），依靠自洽方程既可以评判数据自身的可靠性，也可以用来识别四分量钻孔应变仪的工作状态。若数据段内整体变化趋势一致，则满足自洽方程：S1+S3=S2+S4S1S2S3S4式中：元件 1 的观测数值；元件 2 的观测数值；元件3 的观测数值；元件4 的观测数值。动力灾害的孕育及其诱发与区域应力场的力学行为密切相关，这里的动力灾害主要涉及冲击地压、矿震以及工作面来压期间发生的顶板异常动力显现事件。针对能级跃迁较大的煤矿动力灾害，钻孔应变曲线走势表现为显著的“阶跃”（张月征等，2016b），即数据曲线在动力灾害发生时刻产生台阶式上升或下降。此外，还存在一些能级跃迁较小的动压事件。5 月 4 日5 月 31 日四分量钻孔应变仪分量 1 和分量 3 根据自洽方程拟合后的观测数据（图 2a），9 月 22 日10 月 13 日分量 2 和分量 4 根据自洽方程拟合后的观测数据（图 2b），均满足自洽方程，图中仅绘制拟合后的观测曲线。观察两个不同时间序列上钻孔应变数据的幅值变化，可知某一时刻发生动力异常事件的典型状态失稳现象在应变数据波动特征上并不能够有效直观地显现，应力传递路径受钻孔应变仪安装位置、震源处煤岩介质力学属性等影响，动力异常事件发生前后微小的响应信号大都淹没在固体潮背景下，故有必要根据时频变换识别或提取区域应力场中细微变化与动力灾害的前兆响应特征。1.951.901.851.801.751.701.65应变/1071.931.921.911.901.89应变/10700.510411041.510421042.510431043.5104410400.510411041.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）。定义有意义的瞬时频率其必要条件是函数关于零均值局部对称且过零点的数目与极点的数目一致。据此，Huangetal.（1998）定义了IMF 应满足以下 2 个条件：在整个数据区间内，极值点的数目与过零点的数目相等或至多相差 1 个；在任意一点处，由局部极大值点定义的包络以及由局部极小值点定义的包络的均值为零。此外，为了消除噪声和固体潮等因素对要研究数据段的影响，在 VMD 分解的基础上，采用希尔伯特（Hilbert）变换依次分析 IMF 分量，将时域信息转换至频域信息范畴。按照上述理论对 5 月 1 日5 月 31 日的数据段进行分解，获得各频段 IMF 曲线（图 3）。如表 1 所示，IMF1IMF5 周期皆小于 1h，推测可能是设备自身或周围环境中产生的高斯噪声（Barmanetal.,2016），故剔除 IMF1IMF5 分量，以排除噪声对原始数据的干扰。IMF9 和 IMF10 的频率分别为 2.24105Hz 与 1.16105Hz，转换至周期后分别为 12.4h 和 23.9h，时间维度上与固体潮严格匹配，确定为固体潮的半日波和日波分量。此外，与 IMF9 和 IMF10 邻近的 IMF11 和 IMF12 为固体潮谐波分量。为探究噪声频段对所研频段数据造成的干扰，综上所述，将混杂的各频段数据解析后重构，对比分析降噪前与降噪后的频段曲线（图 4）。降噪信号358地质力学学报https:/2023段与原始信号背景趋势走向一致，对钻孔应变数据的局部异常突跳点捕捉效果良好，降噪信号能够从局部微小信号上反映钻孔应变的波动特征。利用 VMD 对含噪信号按照频率从高到低依次分解，得到若干个固有模态分量。通过构建带通滤波器，滤去高频噪声项和低频固体潮背景趋势项后，将其余 IMF 分量重构。通过上述方法，分析用于监测矿山动力灾害和区域应力场变化的钻孔应变数据。依据变分模态分解原理，既能排除原始信号中的高频噪声，又能剔除低频固体潮背景数据。因此，该方法充分保留了特定信号段非线性、非平稳的固有属性，经频谱分析得到 IMF1IMF12 频率与周期的优势分布，对比每个 IMF 的固有属性和信号时频差异性，从而实现钻孔应变数据的降噪处理。3基于希尔伯特（Hilbert）变换的区域应力异常提取方法冲击地压和矿震都属于煤矿常见的动力灾害，202IMF 100.510411041.510421042.510431043.51044104时间/min202IMF 200.510411041.510421042.510431043.51044104时间/min202IMF 300.510411041.510421042.510431043.51044104时间/min202IMF 500.510411041.510421042.510431043.51044104时间/min（a）采集环境中产生的高斯噪声频段(IMF1IMF5)1012IMF 400.510411041.510421042.510431043.51044104时间/min42024IMF 600.510411041.510421042.510431043.51044104时间/min101IMF 800.510411041.510421042.510431043.51044104时间/min（b）与矿震相关频段(IMF6IMF8)24024IMF 700.510411041.510421042.510431043.51044104时间/min第3期吴浩源，等：煤矿深部开采动力灾害区域应力效应监测与前兆特征研究359两者既有相似之处，也有其内在的差异。通常由于积聚的弹性应变能瞬间释放致使周围煤岩体的破坏，从而引起冲击地压和矿震。前者多发生在巷道或工作面以及开采扰动近源区，并伴随有煤岩体抛出的动力现象，而后者是开采扰动与区域应力场共同作用诱发的地震，多发生在工作面开采扰动远源区，其在表征动力破坏程度上弱于冲击地压，且不发生煤岩体的抛出现象。此外，从二者的发生频次上分析，冲击地压灾害数要远小于矿震，但当满足一定的临界条件时，矿震会诱发冲击地压灾害的发生。因此，区域应力场监测是实现矿震预测、预警的基础和前提，且有利于揭示冲击地压的致灾机理（姜耀东等，2014）。为了得到区域应力场的实时变化，将经过 VMD分解后得到的精细信号段进一步分析和处理。常用的方法有傅里叶变换和希尔伯特变换。其中，傅里叶变换的优势表现为对平稳或线性信号的处理，若使用该方法处理非平稳、非线性的信号，则需进行特定的变形后方可使用，如：小波变换和短时傅里叶变换。但经过小波变换和短时傅里叶变换解析的信号会出现负频率，负频率在分析相关问题时不具有任何物理意义，而希尔伯特变换在解析过程表1钻孔应变数据 VMD 分量频率周期表Table1FrequencyperiodictableofVMDcomponentofboreholestraindataIMF123456频率/Hz3.981032.621031.191036.641043.171042.41104周期/h0.070.100.230.410.891.15IMF789101112频率/Hz1.531046.371052.241051.161051.511067.49107周期/h1.84.312.423.91853701201IMF 900.510411041.510421042.510431043.51044104时间/min505IMF 1000.510411041.510421042.510431043.51044104时间/min051510IMF 1200.510411041.510421042.510431043.51044104时间/min505IMF 1100.510411041.510421042.510431043.51044104时间/min（c）固体潮及相关谐波频段(IMF9IMF12)图35 月 1 日5 月 31 日钻孔应变数据 VMD 分解曲线图Fig.3VMDdecompositioncurves(IMF1IMF13)ofboreholestraindatafromMay1sttoMay31st7.497.487.477.467.457.447.437.427.417.407.39应变/1070200400600800100012001400时间/min含噪信号降噪信号图4降噪前后对比效果图Fig.4Comparisondiagrambeforeandafternoisereduction360地质力学学报https:/2023/2中虚数部分较实数部分滞后，借用此特性可消除负频率的影响。因此，希尔伯特变换在解析非线性、非平稳信号时，有其自身特有的优势。根据矿震事件发生的时间，以分钟为单位提取钻孔应变数据，经希尔伯特变换后得到矿震分量的瞬时幅值、瞬时能量和每日累计能量（图 5，图 6）。由 5 月 1 日 5 月 31 日 和 9 月 22 日 10 月 13 日的钻孔应变瞬时幅值时程变化曲线可知（图5a，图6a），在时间序列内统计矿震发生频次共计 6 次，分别为5 月 21 日、5 月 26 日 和 5 月 29 日 及 9 月 24 日、10 月 4 日和 10 月 11 日，从瞬时幅值和累计能量两个角度均反应了相同的变化规律。其中，日累计能量的变化趋势说明区域应力场内能量在不断的积聚、释放和耗散（图 5b，图 6b）。1.61.81.41.21.00.80.60.40.20日累计能量/1010J5 月 4 日5 月 5 日5 月 6 日5 月 7 日5 月 8 日5 月 9 日5 月 10 日5 月 11 日5 月 12 日5 月 13 日5 月 14 日5 月 15 日5 月 16 日5 月 17 日5 月 18 日5 月 19 日5 月 20 日5 月 21 日5 月 22 日5 月 23 日5 月 24 日5 月 25 日5 月 26 日5 月 27 日5 月 28 日5 月 29 日5 月 31 日5 月 30 日日期4.03.53.02.52.01.51.00.5应变/1040.510411041.510421042.510431043.51044104时间/min5 月 21 日事件5 月 26 日事件5 月 29 日事件（a）“突跳”“震荡”“跌落”前兆响应模式（b）日累计能量分布图图55 月 1 日5 月 31 日重构分量动力异常响应幅值和能量时程变化Fig.5DynamicanomalyresponseamplitudeandenergytimehistorychangeofreconstructedcomponentsfromMay4thtoMay31st2.01.81.61.41.21.00.80.60.40.2应变/10400.510411041.510421042.51043104时间/min0.350.400.300.250.200.150.100.050日累积能量/1010J9 月 22 日9 月 28 日10 月 3 日10 月 8 日10 月 12 日9 月 26 日10 月 1 日10 月 6 日10 月 11 日9 月 24 日9 月 29 日10 月 4 日10 月 9 日9 月 23 日9 月 27 日10 月 2 日10 月 7 日9 月 25 日9 月 30 日10 月 5 日10 月 10 日10 月 13 日日期9 月 24 日事件10 月 4 日事件10 月 11 日事件（a）“突跳”“震荡”“跌落”前兆响应模式（b）日累计能量分布图图69 月 22 日10 月 13 日重构分量动力异常响应幅值和能量时程变化Fig.6DynamicanomalyresponseamplitudeandenergytimehistorychangeofreconstructedcomponentsfromSeptember22ndtoOctober13rd 3.1从幅值曲线波动形态上分析矿震与地震都是地下岩层运动后产生形变进而导致弹性应变能的积累并释放（蒋骏等，2000；彭华等，2008，2011；王凯英，2018），两者在发生机制上具 有 相 似 的 特 征。Mescherikov（1968）、Yoichiro（1974）通过对大量地震资料的归纳总结发现，在地震的孕育和诱发过程中存在、和 三个阶段，运用该理论在后续的地震跟踪预测中取得了显著效果。为解释地震形变三阶段的物理实质，冯德益（1983）认为 阶段以岩层的微小破裂为主，持续扩展乃至相互贯通，导致最终发生非震滑动或岩层的蠕变；阶段则是因断层活化而引起的突变异常，并最终进入 阶段，表征为大量积蓄弹性应变能的释放。如图 7 所示，地震形变三阶段的曲线走势与岩石破坏的应力应变曲线相似，运用该理论分析煤矿矿震的活动规律，以期探究矿震发生的前兆特征。以 5 月 21 日矿震事件为例，应变幅值“有升有降”走势契合地震发生的前期、中期 和后期 三个形变阶段。阶段：应变稳定变化，该阶段虽然也会造成煤岩破裂，但区域应力场总体稳定，能量积聚速率与能量释放速率趋于一致，不足以引发动力灾害。阶段：应变急剧变化，表现形式为“震荡”+“突跳”，这是相关区域岩层破裂产生非稳态应力调整现象的表征，能 量 积 聚 与 能 量 释 放 处 于 临 界 状 态，如 5 月21 日矿震前期先发生“震荡”后发生“突跳”，而5 月 26 日矿震 阶段应力异常表征为“震荡”、“震第3期吴浩源，等：煤矿深部开采动力灾害区域应力效应监测与前兆特征研究361荡”+“突跳”2 种应变异常表现形式受矿震震源处和“力链”范围内煤岩介质属性、应力传递路径、开采扰动影响范围，以及震源、工作面和钻孔应变安装设备三者相对位置的影响，区分 2 种异常形式有利于了解区域应力场的赋存状态，并为未来的煤矿开采工作提供技术指导。阶段：大量积聚弹性应变能瞬间释放，造成煤岩体破裂并引发矿震。应变时间O图7地震形变三阶段物理模型（宋治平等，2004）Fig.7Three-stagephysicalmodelofseismicdeformation(Songetal.,2004)3.2从形变异常开始时间分析钻孔应变数据经 VMD 识别并提取能够反映矿震的 IMF 分量。该 IMF 分量经希尔伯特变换后所得时间幅值变化关系，亦是矿震震源、工作面开采位置和钻孔应变仪安装位置三部分组成的系统所形成的一定空间范围内区域应力场变化的表征。不同于天然地震发生前 2 年应变出现异常的规律（宋先月等，2003)，解析出矿震形变异常早于其发生27 天（表 2，表 3），这是天然地震与煤矿矿震的孕震环境不同所造成的。工作面开采扰动促使区域应力场不断的调整，对比观察开采扰动作用下的区域应力场幅值随时间的变化，形变的三个阶段中稳定变化阶段（）的应变格值为 0.51010，形变快速变化阶段（）的应变格值为稳定变化阶段的 3 倍以上，呈现为“突跳”的变化特征。动力灾害发生前，区域应力场内应力（应变）显著升高，积蓄大量的弹性应变能，当相关区域的可释放弹性应变能超过岩石表面应变能时，即可沿破裂面瞬间释放。基于鲍店煤矿典型动力灾害矿震事件，通过分析钻孔应变数据对区域应力场异常的响应，认为应变幅值由形变稳定变化阶段向快速变化阶段的幅值阶跃为形变异常开始时间。随后，区域应力场在开采扰动作用下不断调整直至进入震前失稳阶段，能量释放呈现出动力灾害形式。因此，对应变数据进行 VMD 和 Hilbert 变换后，可以识别并提取能够反映煤矿动力异常响应的 IMF 分量，并在地震形变三阶段的理论基础上，统计分析其幅值和累积能量的时程变化规律，进而尝试建立煤矿动力灾害的预测预报模式。4存在的不足与展望（1）基于鲍店煤矿所在区域地应力监测数据，采用 Hilbert 变换方法进行了解析，获得了煤矿动力灾害发生前明显的前兆异常响应特征，但不能代表不同构造环境下的应变前兆特征，未来还需进一步丰富不同构造环境下的应变监测数据采集与分析工作，以期归纳总结出煤矿动力灾害的预测预报模式。（2）在对分解出的若干个 IMF 分量进行异常信号提取和识别工作中，由于煤矿井下环境复杂多变，设备采集到的原始数据受到严重干扰，以上研究中所采用的方法仅是针对具体问题的定向研究，在一定程度上删除了干扰信号的影响，但是使用的算法有其自身的局限性，存在部分数据段无法消除干扰的问题，后期还需进一步优化算法，以期提高数据分析的准确性。5结论采用 VMD 和 Hilbert 变换的方法，分析了鲍店煤矿开采期区域应力场变化的钻孔应变监测数据，结论如下：（1）将 VMD 方法应用到钻孔应变数据处理中，能够将混频信号分解为若干个具有特定频段的固有模态分量，为后续动力异常响应特征的识别和提表25 月份形变异常前兆模式特征表Table2CharacteristicsofdeformationanomalyprecursorymodeinMay发生时间形变异常开始时间应力调整持续时间表现形式5月21日5月11日10d“震荡”+“突跳”5月26日5月24日2d“震荡”5月29日5月27日2d“突跳”+“震荡”表310 月份形变异常前兆模式特征表Table3CharacteristicsofdeformationanomalyprecursorymodeinOctober发生时间形变异常开始时间应力调整持续时间表现形式9月24日9月22日2d“突跳”10月4日9月27日7d“突跳”10月11日10月6日5d“突跳”+“震荡”362地质力学学报https:/2023取提供无模态混叠的数据源。（2）通过构建带通滤波器，消除高斯噪声和应变固体潮干扰项，重构信号获得能够反映区域应力场对矿震、冲击地压等动力灾害响应的 IMF 分量，并采用 Hilbert 变换方法处理特定频段 IMF 分量。（3）依据地震形变三阶段理论，定义形变异常开始时间为区域应力场进入异常状态时刻，并经过形变快速变化阶段的“震荡”和“突跳”后，释放出积聚的弹性应变能，满足岩石强度失效和变形破坏准则时，即能引发动力灾害。（4）基于钻孔应变仪对区域应力场的精确观测，将形变异常开始时间作为动力灾害事故发生前的响应时刻，为未来煤矿开采工程的安全性和稳定性评价提供了一种切实可行的前兆信息判别模式。ReferencesBARMANC,GHOSED,SINHAB,etal.，2016.Detectionofearthquakein-ducedradonprecursorsbyHilbertHuangTransformJ.JournalofAp-pliedGeophysics，133：123-131.CHENK,LIY,CHENL，2010.Ensembleempiricalmodedecompositionforpower quality detection applicationsJ.Computer Simulation，27（3）：263-266.（inChinesewithEnglishabstract）CHISL，2013.StrainanomaliesbeforeWench