随着开采深度的不断增加，煤岩体在高地应力条件下的破坏与失稳行为日趋复杂化，动态灾害风险显著升高。尽管微震（MS）监测技术已得到广泛应用，但现有方法主要聚焦于事件能量、频次等单一维度指标，难以准确捕捉煤岩体破裂失稳的完整演化过程及其前兆特征。本研究基于平煤八矿己15-21050工作面近五个月的微震监测数据，系统分析了微震信号频谱特征、能量响应与裂隙结构之间的时空耦合关系。研究结果表明：（1）微震事件频谱以低频单峰型为主导，表明破裂过程受大尺度、长历时事件控制，与长壁开采中的渐进式应力释放机制一致；（2）高频复杂频谱虽占比较低，但集中于构造活跃区域，对局部应力异常和结构活化具有高度敏感性；（3）频谱复杂性指数（SCI）在失稳孕育阶段达到峰值，显著早于能量积聚和裂隙连通性增强阶段，显示出对微裂纹活化的强早期预警能力和敏感性；（4）裂隙方向集中度指数（Dc）在强能量释放前显著增大，表明其能够有效捕捉能量释放的几何响应特征，可作为强失稳的有效物理前兆。基于上述发现，研究人员构建了熵权多维风险预测模型，现场应用显示煤岩体失稳风险评分由2023年12月的0.27升至2024年2月的0.76，与2024年3月煤炮事件、顶板变形及高能微震活动的时空分布高度吻合，验证了模型的预测准确性和工程适用性。
## 研究背景、问题与研究意义
随着煤炭开采深度的持续增加，原位应力水平显著升高，煤岩体的力学响应和破坏模式发生明显改变，呈现出与浅部开采截然不同的特征。在高应力和强采动扰动的耦合作用下，煤岩体更易发生失稳，可能诱发冲击地压、煤与瓦斯突出等动态灾害，对矿山安全高效生产构成严重威胁。因此，准确捕捉煤岩体破裂与失稳的演化规律，实现动态灾害的早期预测，仍是深部开采面临的重大挑战。
为阐明深部动态灾害的复杂致灾机理，学术界提出了多种理论，包括刚度理论、能量理论、三准则理论、"三因素"理论、动静载叠加诱冲理论以及冲击地压与煤与瓦斯突出统一失稳理论等。这些基础研究从能量积聚与释放、系统稳定性、应力-强度耦合等不同角度揭示了煤岩体失稳的内在机制。同时，针对煤与瓦斯突出机理，也提出了 gas-dominated、stress-dominated、化学主导、流变突变、粘滑失稳及球壳失稳等多种代表性假说。这些理论进展为动态灾害预测方法的发展奠定了坚实基础。
随着矿山安全保障需求的增长，研究重心逐渐从纯机理分析转向灾害监测与预警。监测方法从传统的钻屑分析、煤岩变形测量等岩石力学手段，发展到微震（MS）、声发射（AE）、电磁辐射（EMR）、视电阻率等动态连续的地球物理技术。微震监测因其实时性强、空间覆盖范围广、抗干扰能力强等优势，已成为揭示煤岩体破裂演化与能量释放机制的关键技术，广泛应用于煤矿、隧道、水电站等地下工程结构的稳定性评估。
然而，现有基于微震监测的煤岩体失稳风险预测方法主要聚焦于微震事件的时间、空间、震级等特征的定性或定量分析，发展出时间集中度、活动性指数、震源集中度、b值、分形维数、地震平静期、能量指标、速度异常、视应力与视体积等一系列衍生指标。尽管机器学习、聚类分析、降维等智能算法的引入极大提升了微震数据的智能化处理水平，但现有指标和方法对微震波形中蕴含的复杂物理机制和频谱特征关注有限。微震波形蕴含丰富信息，短时傅里叶变换（STFT）、小波包变换、希尔伯特-黄变换（HHT）等时频分析技术已被广泛用于提取不同频带的能量分布，但现有研究多集中于优势频率、频带能量等单一时频参数的提取与应用，对微震频谱特征的时空分布探索不足，更缺乏频谱复杂性、能量与裂隙几何之间耦合协同演化的系统性研究。这一局限性制约了非平稳煤岩体破坏过程的动态精准识别，降低了灾害预测模型的时效性和适用性。
针对上述挑战，本研究以河南神马股份有限公司八矿己
15-21050工作面为工程背景，基于近四个月的连续微震监测数据，开展微震波形的频谱特征提取、频谱聚类和空间演化映射，创新性地将频谱特征纳入能量场和裂隙几何的时空动态分析，提出了频谱复杂性指数（SCI），建立了融合频谱、能量与裂隙几何特征的多维指标体系，构建了熵权多指标煤岩体动态失稳预测模型。
## 主要技术方法
本研究采用的关键技术方法包括：（1）基于快速傅里叶变换（FFT）与统计描述符的微震信号频谱特征提取，包括优势频率（fd）、频谱峭度（K(f)）和频谱偏度（S(f)）；（2）滑动能量比、局部峰值跳跃、均方根（RMS）增强比及到达时一致性约束等多准则信号筛选方法；（3）基于K-means算法的微震信号无监督聚类分析，结合轮廓系数确定最优聚类数；（4）基于谱特征调制的局部微震能量场计算方法，引入多因素扰动半径调整函数Rsi替代传统恒定扰动半径；（5）基于时间、频谱、空间、方向、能量等多约束条件的裂隙追踪模型；（6）基于信息熵的频谱复杂性指数（SCI）和基于圆统计学的裂隙方向集中度指数（Dc）；（7）熵权法（EWM）确定多指标权重，构建综合风险预测评分模型。样本数据来源于平煤八矿己
15-21050工作面2023年11月至2024年2月期间的现场微震监测数据。
## 研究结果
### 微震信号频谱特征的时间演化
煤岩体失稳早期频繁发生的小尺度、低能量破裂事件具有破裂历时短、能量释放迅速集中的特点，表现为较高的优势频率；随着破坏向大尺度、长历时能量释放过程发展，优势频率逐渐降低。频谱峭度K(f)量化了局部频带内能量的瞬态集中程度，高K(f)值指示加速破裂或瞬时微裂纹扩展；频谱偏度S(f)反映频谱对称性，正值表示右偏（高频主导）谱，与快速剧烈破裂相关，负值表示低频成分增强，对应渐进或低能量破裂过程。基于上述参数，研究将微震事件划分为单峰型、多峰型、频率拖尾型、高频快速衰减型和低频持续型五种典型频谱类型。
在2253个记录事件中，2231个被成功识别分类。统计分析表明，低频单峰型频谱占主导地位，表明观测期间煤岩体破裂过程主要由长历时、慢破裂事件控制，与采煤工作面推进过程中的渐进式能量释放和应力卸压机制相符。中高频、多峰及复合频谱类型虽占比较小，但指示了局部复杂破裂机制的存在和多样的能量释放行为。
### 微震频谱聚类特征的时空演化
采用K-means无监督聚类方法，基于峰值频率fpeak、频谱峭度K(f)、频谱偏度S(f)、低/中/高频带相对能量比R、优势频率持续时间Tp和频谱形态类型Shape等构建多维特征向量，将微震事件划分为四类：K1型（集中低频型，稳定低能量响应、塑性变形）、K2型（高频衰减型，快速能量释放、脆性破坏）、K3型（频率拖尾型，持续破裂、非均质弱结构）和K4型（复杂多峰型，多源激发、累积扰动）。
构建工作面推进过程中四个阶段的频谱-空间聚类图，揭示不同阶段微震事件的阶段性响应特征：第一阶段（推进约70 m）以K1型低频单峰频谱为主，空间分布集中于工作面前方0-100 m范围内，聚类密度低、空间结构稳定，整体应力场处于平衡状态，但在100-250 m构造弱化带出现K3型"频率拖尾"事件的局部聚集；第二阶段（推进157.6 m）扰动强度显著增加，微震活动更为频繁，K3、K4型高频及复杂频谱事件呈块状群集出现，沿工作面形成多个密集聚类核，系统进入能量积聚和裂隙活跃阶段；第三阶段系统暂时稳定，K1型频谱重新占主导，但210-300 m构造带残留K2型聚类，表明边界应力效应的持续影响；第四阶段频谱响应进一步减弱，K1型事件优势增强，但在110-250 m范围内存在条带状K3、K4型残留聚类，呈现"滞后响应"特征。西部断层带和构造带在所有阶段均持续呈现K2-K4型危险聚类条带，代表系统性应力集中和裂隙传播区域，是微震聚类和动态灾害的主要激发源。
### 微震能量场的时空演化
基于谱特征调制的局部能量场计算方法，构建多阶段能量场并进行分析。初始阶段能量释放整体较低，能量场呈局部、孤立分布；随着推进至157.6 m，能量场显著增强，分布扩展为多个条带状区域，峰值能量超过10
4 J；至2024年1月，能量场进一步增至10
4-10
5 J，主高能区逐渐向工作面中部收敛，呈现"少事件、高能量"趋势；至2月下旬，能量场沿采掘方向形成连续高能带，峰值接近10
5 J，事件密度显著增加，空间分布与煤层厚度异常和断层结构高度耦合。统计表明，大于10
4 J的高能事件占比从第一阶段0%增至第四阶段6%，能量场峰值强度从2023年12月约10
4 J升至2024年2月约10
5 J。
### 煤岩体裂隙的多约束识别与演化特征
基于时间、频谱、空间、方向和能量等多约束条件构建裂隙追踪模型，约束条件包括：空间距离≤100 m、传播方向偏转角≤90°、相邻事件时间间隔≤1 d、相邻事件能量比≤5、分支仅允许在偏转角小于30°且时间间隔不超过1 d时发生。从2023年11月至2024年2月数据中，共识别出67条有效裂隙。
统计分析显示：2023年11月裂隙数量多但平均长度较短（57.7 m），方向集中度Dc=0.254，中位能量较低，空间分散；2023年12月总裂隙长度显著增至约1065 m，Dc升至0.454，裂隙方向 increasingly clustered along NW-SE方向；2024年1月平均裂隙长度升至81.9 m，中位能量约1140 J，Dc降至0.221，总裂隙长度和能量持续增加，几何连通性增强；2024年2月中位能量升至约2350 J，Dc降至0.106，裂隙规模和能量持续增加但方向一致性减弱，转入多向扩展阶段。Dc增大通常先于能量峰值出现，表明裂隙连通性可作为强能量释放的前兆。
### 频谱复杂性、能量场与裂隙演化的耦合特征
构建频谱-能量-裂隙耦合分析图，揭示三者的交互作用。高能区（>10
4 J）与K2-K4型复杂及高频微震事件高度吻合，约72.4%的K2-K4型事件分布于这些高能区域，表明高频复杂频谱信号主要反映应力集中区的快速、局部化能量释放。裂隙连通性相对于频谱演化存在明显的时间滞后：第二阶段复杂频谱事件密集分布、频谱结构复杂性增加时，裂隙数量略减、总裂隙长度缓慢增长；第三阶段复杂频谱事件比例下降约40%，但总裂隙长度较前一阶段增长超过60%，表明微震频谱变化主要反映破裂起始和微裂纹激活，而裂隙连通和能量集中为后续响应，呈现"频谱响应先于结构发展"的异步演化模式。能量场与裂隙系统的空间耦合进一步显示，高能带与各阶段优势裂隙方向（如12月沿NW-SE方向聚集、1月向N-NW方向旋转）广泛对应；第三、四阶段高能带与优势裂隙方向部分重叠，能量峰值达约10
5 J时Dc降至0.1左右，表明裂隙传播更趋无方向性，主应力方向影响减弱，能量呈多向扩散，系统由早期同步耦合转向后期部分解耦。
基于信息熵的频谱复杂性指数SCI，结合频谱聚类熵和K2-K4型高频复杂频谱事件比例计算，定量表征频谱分布的非均匀性和演化复杂性。对比分析表明，2023年11月至12月SCI持续升高，反映频谱结构复杂性显著增强，系统处于微裂纹激活和应力扰动加剧的孕育阶段；2024年1月后SCI快速下降，而能量场峰值和总裂隙长度同步增加，反映系统由局部扰动向宏观失稳转变。SCI峰值显著早于高能区形成和裂隙几何连通阶段，展现出对煤岩体失稳演化阶段的强前瞻性和敏感性。
### 熵权多指标驱动的煤岩体失稳预测
基于频谱-能量-裂隙耦合演化规律，建立包含频谱复杂性、能量释放强度和裂隙演化特征的多维层次化风险预测指标体系，采用熵权法客观确定各指标权重。指标包括：频谱复杂性指数（SCI）、SCI变化率（SCI-VR）、中位能量（MIE）、最大能量（MXE）、能量场峰值（EFP）、总裂隙长度（TFL）、平均裂隙长度（AFL）、裂隙增长率（FGR）和裂隙方向集中度（Dc）。熵权分析表明，能量场峰值（权重0.1636）和最大能量（权重0.1589）是最主要的预测指标，SCI（权重0.0912）及其变化率（权重0.1177）在早期破裂阶段展现出高敏感性。
煤岩体失稳风险等级划分为：I级安全（0-0.3）、II级注意（0.3-0.6）、III级警告（0.6-0.8）、IV级高风险（>0.8）。应用于己
15-21050工作面微震数据，风险评分由2023年12月的0.27快速升至2024年2月的0.77，风险等级由I级升至III级，反映系统稳定性持续劣化、失稳风险显著上升。
### 预测结果验证
2024年3月现场跟踪调查期间，月推进速率由2月150 m降至130 m，高能（>5×10
4 J）微震事件由7次降至4次，整体动态活动看似减弱，但现场记录到三次连续可闻煤炮，地下勘查发现异常区段存在围岩结构与支护系统协同破坏的复合失稳实例，确认该段已进入高风险破坏阶段。微震能量时间演化显示，3月日均能量虽较2月峰值略有下降，但仍显著高于2023年11月至2024年1月水平，多日保持高能状态，与模型预测的风险指数升至III级趋势高度吻合。现场数据证实，该无标签风险预测模型可提前约1个月识别煤岩体失稳趋势，且通过调整滑动时间窗口参数支持周滚动预测。
## 讨论与结论
研究讨论部分系统对比了本研究与以往工作的差异。现有研究主要建立了基于有限数量预警指标的煤岩体动态风险识别模型，如微震事件数量、能量释放、b值、潜在最大震级等统计或经验指标，以及多指标加权与模糊评价方法，但对频谱复杂性、能量场演化和裂隙几何响应的表征仍显不足。近年来虽有研究引入机器学习和深度学习方法进行微震事件识别和动态灾害风险预测，如随机森林（RF）、支持向量机（SVM）、深度卷积神经网络（DCNN）自动分类微震事件类型或空间聚类，无监督聚类方法进行频谱特征降维和震源机制分类，以及长短期记忆（LSTM）和卷积神经网络（CNN）模型进行微震-声发射-电磁辐射多源信号联合分析与预警，但这些研究普遍仍缺乏将频谱复杂性、能量场特征与裂隙结构参数在统一评价体系中进行定量耦合分析。
本研究提出的多维熵权融合模型将SCI及其变化率与能量场峰值、能量统计特征以及总裂隙长度、裂隙增长率、方向集中度等结构指标纳入统一评价体系，客观确定各指标权重。应用结果表明，该模型能够敏感反映煤岩体从相对稳定状态到高风险阶段的演化，并对高风险阶段的起始呈现早期响应。相较于传统仅依赖有限能量或统计特征的方法，该频谱-能量-裂隙多源信息融合评价框架在特征构建的物理完备性、早期失稳识别和滚动预测适应性方面具有优势，为未来工程案例的区域风险分区和动态预警提供了可扩展的方法论基础。该模型适用于高应力、微震活动显著、裂隙发育复杂的采掘区域，在稳定区域微震活动较少时性能可能受限，未来需针对不同地质条件进行定制化调整和验证。
研究结论如下：
（1）己
15-21050工作面开采过程中微震频谱分析表明，煤岩体破裂以低频单峰型事件为主，高频复杂频谱模式随扰动增强而增加。基于谱特征的局部能量建模方法揭示了从孤立能量结构到条带能量结构的演化，高能区与煤层厚度异常和断层影响区空间吻合。裂隙识别结果进一步证明了从多点激活到定向传播再到高能释放的阶段性演化。
（2）频谱复杂性、能量与裂隙几何的时空耦合呈现明显的阶段依赖性关系。频谱复杂性在能量峰值和裂隙连通性之前增大，表现出"频谱先行、结构滞后"的模式。随着扰动加剧，能量峰值与优势裂隙方向吻合，表明强耦合；在高能释放阶段，裂隙传播方向性减弱，能量多向扩散，系统逐渐由耦合转向部分解耦演化。
（3）提出的裂隙方向集中度指数Dc和SCI均能有效进行早期失稳识别。Dc通常在能量峰值前上升，表明其可作为潜在前兆。SCI在能量积聚和裂隙连通之前达到峰值，展现出强敏感性和早期预警潜力，为煤岩体失稳预测提供了可靠的物理依据。
（4）通过融合频谱、能量和裂隙特征的熵权多指标赋值，构建了无标签煤岩体失稳风险预测模型。与以往聚焦单一维度指标的模型不同，该模型能够捕捉煤岩体失稳的完整演化过程。在己
15-21050工作面案例中，风险评分由2023年12月的0.27升至2024年2月的0.76，随后2024年3月的现场监测识别出煤炮、顶板变形和大量高能微震事件，验证了模型的准确性和工程适用性。
