高应力区微震监测信号特征分析.pdf
第 2 4 卷 第 2期 2 0 0 7年 O 6月 采矿与安全工程学报 J o u r n a l o f Mi n i n g& S a f e t y En g i n e e r i n g Vo 1 . 24 NO . 2 J u n e 2 0 0 7 文章 编号 1 6 7 3 3 3 6 3 2 0 O 7 O 2 0 1 4 6 0 4 高应力区微震监测信号特征分析 曹安业 ,窦林 名 , 秦 玉红 ,李志华 ,巩 思 园,王 玉刚 中国矿业大学 煤炭 资源 与安全 开采 国家重点实验室 ,江苏 徐州 2 2 1 0 0 8 摘 要 如 何识 别不 同微震 信 号 , 尤其是 快 速识 别 冲 击 危 险信 号 , 以便 及 时 采取 解 危措 施 , 是 微 震 监 测 中亟待 解 决的 问题 . 以三 河尖矿 9 2 0 2采 面过 上覆 煤柱 高应 力 区为 背景 , 通过 微 震信 号 的 时一 频分析技术 , 总结提炼了不同微震信号的重要波形特征. 研究表明 不同微震信号在频率特征、 信 号持 续时间、 释放能量等多方面都存在 区别, 冲击矿压 波形有震幅大、 衰减快、 尾 波较发育、 频带 分布较窄等重要特征. 通过对微震信号的特征分析 , 有助于深入揭示微震监测的震源机理 , 进一 步推 动 冲击矿 压 的预 测预报 工作 . 关键词 冲击矿压 ; 微震;预测预报 ; 动力灾害;幅频特征;时一 频分析 中 图分类 号 T D 3 2 4 文献 标识 码 A Ch a r a c t e r i s t i c o f M i c r o s e i s mi c M o n i t o r i n g S i g n a l i n Hi g h S t r e s s e d Z o n e CAO An y e ,DOU Li n mi n g,QI N Yu h o n g ,LI Z h i h u a ,GONG S i y u a n,W ANG Yu g a n g S t a t e Ke y L a b o r a t o r y o f Co a l Re s o u r c e a n d M i n e S a f e t y ,Ch i n a Un i v e r s i t y o f Mi n i n g Te c h n o l o g y ,Xu z h o u,J i a n g s u 2 2 1 0 0 8 ,Ch i n a Ab s t r a c t I n a n hi g h s t r e s s e d a r e a,Ho w t O i de nt i f y d i f f e r e n t m i c r o s e i s m i c moni t o r i n g s i gn a l s, e s pe c i a l l y ho w t O r e c o g ni z e t ho s e d a n ge r ous s i g na l s f l e e t l y a n d t a ke me a s ur e s i n t i m e i s a n ur g e nt p r ob l e m t o b e s ol v e d. Ba s e d o n t he i ns t a n c e o f No. 9 2 02 wo r k i ng f a c e a c r o s s t he up pe r r e ma i n e d c o a l p i l l a r i n S a n h e j i a n Mi n e ,t h e c h a r a c t e r o f s i g n a l s c a u s e d b y d i f f e r e n t h y p o c e n t e r s ha d b e e n s umm a r i z e d a nd a bs t r a c t e d b y us i n g t he t i me f r e q ue nc y a n a l ys i s me t h od . The r e s u l t s s h o w t h a t t h e s i g n a l s a r e d i f f e r e n t i n f r e q u e n c y,d u r a t i o n,r e l e a s e d e n e r g y,a n d S O o n a n d t h e wa v e f or m o f r o c k b ur s t i s c h a r a c t e r i z e d by b i g s e i s m i c a m p l i t ud e,f a s t t a p e r o f f ,de v e l o pe d t a i l - wa v e a nd na r r o wl y d i s t r i b ut e d b a nd wi d t h. Thi s r e s ul t i s u s e f uI i n pr o mot i n g t he f or e c as t o f r oc k bur s t . Ke y wo r ds r o c k b ur s t;mi c r os e i s m ;f or e c a s t ;d yn a mi c c a t a s t r op he;a mpl i t u de f r e qu e n c y c ha r ~ a c t e r ;t i m e - f r e qu e n c y a na l y s i s 随着我国矿井开采规模和深度的逐年增加 , 煤 矿中发生震动与冲击矿压越来越频繁. 目前 , 对于 冲击矿压的研究主要集中在现场 的预测预报与治 理. 由于冲击矿压 的随机性 、 突发性以及破坏形式 的多样性 , 使得预测预报工作变得非 常复杂 , 预测 的准确性受到很大的影响 ≈ ] . 因此 , 寻求科学可靠 的冲击矿压预测 预报方法 , 提出针对性 的防范措 施 , 是我 国煤炭工业急需解决的课题. 微震监测技 术是利用岩体受力变形和破坏后本身发射 出地震 波来进行监测工程岩体稳定性的技术方法. 该技术 是岩体动力灾害, 特别是冲击矿压预测预报最有效 和最有潜力的监测方法之一. 收稿 日期2 0 0 6一 O 91 O 基金项目国家 自然科学基金项 目 5 0 4 9 0 2 7 3 , 5 0 4 7 4 0 6 8 , 5 0 0 7 4 0 3 0 ; 国家重点基础研究 发展规划 9 7 3 项 目 2 0 0 5 CB 2 2 1 5 0 4 ; 教育部博士点基金项 目 2 0 0 3 0 2 9 0 0 1 7 作者简介 曹安业 1 9 8 2 一 , 男 , 江苏省盐城市人 , 博士, 从事矿 山压力 、 冲击矿压 、 采矿地球物理等方面的研究. E- ma Uc a o a n y e 1 6 3 . c orn Te l 1 3 7 7 5 9 8 5 9 9 5 维普资讯 第 2 期 曹安业等 i高应力区微震监测信号特征分析 1 微震监测与信号处理基本原理 1 . 1 微震 监 测冲 击危 险原 理 煤岩体在开采影响下发生破裂、 滑动过程 中, 会产生一定频率地震波向周 围传播 , 释放聚集的能 量 , 而波的振幅和频率又取决于煤岩体 的强度 、 应 力状态、 破裂尺寸和变形速度等_ 3 ] . 因此 , 每一个微 震 声发射 信号都包含有岩体内部状态变化的丰 富信息 ] . 而冲击矿压的孕育过程实际上就是煤岩 体结构内部损伤裂纹或裂隙发生与扩展 的过程; 发 生冲击则是损伤裂纹或裂隙的失稳扩展与断裂的 过程_ 5 ] 。 通过在空间上不 同方位设置微震拾震器 , 对微震事件进行实时监测 , 进行震源定位 , 并对接 收到的微震信号进行处理、 分析, 可确定各种震动 参数 能量, 频次, 半径 , 力矩 , 压降等 , 以此为基础 判断推理煤岩体应力状态及破坏情况 , 进而评价和 预测冲击矿压等矿井动力灾害危险 ] . 1 . 2 微震信号的时一 频分析技术 谱分析已成为微震研究的一种标准方法; 采用 时一 频分析技术分析微震信号 的功率谱和 幅频特 性, 以便从谱特性进行微震信号 的辨识 , 从 而为矿 井冲击矿压预测预报等动力灾 害提供一条新的线 索 ] . 时间域 内的地震波形模拟分析需要相 当复杂 的技术 , 尚不能做到常规应用. 地震记录或 时间序 列经快速傅里叶变换 F F T 为频率域, 便可得到所 需的振幅谱和相位谱 , 就可求得大部分震 源参数 部分参数在时间域 内 , 亦可确定信号 的频 率特 征, 掌握信号的构成及性质. 把原来在 时域 内以时间 t为变量 的函数 B“ 变换为频域 内以频率 _厂为变量 的函数 B 厂 , 也 就是将原来 的函数分解为一系列振幅不 同的频率 变化的正弦函数 , 得出频域内振幅随频率变化 的 函数 8 厂 . T r B 厂 一 l B H £ e - 2p f t d t , 1 { 式中t 为 B £ 在时域 内延伸的区间. 由于振幅的平方正 比于功率 , 定义单位时间内 功率谱密度 .S 厂 为 S 厂 一 B 厂 。 / T . 2 信号的分析过程通过 ma t l a b完成. 2 不 同类型微震信 号的分析与识别 2 . 1 三 河尖 矿微震 监 测 系统布 置 徐州三河 尖矿是座 冲击灾害严 重 的矿井 , 自 1 9 9 1 年 5月首次在 7 1 1 0面材料 道发生 冲击矿压 以来, 至今已发生破坏性冲击矿压 3 0 余次. 目前, 该矿 9 2 0 2面采深 已超过 8 5 0 m, 顶板为厚层坚硬 砂岩 , 且上覆有 7煤 7 2 0 2面残 留的大煤柱区, 受煤 柱应力传递影响 , 9 2 0 2面过煤柱 区域 时应力 集 中 程度严重 , 冲击危 险性大. 为更好预测预报 冲击矿 压 , 三河尖矿建立了微震监测系统 , 对 9 2 0 2面进行 实时监测 , 并较好地指导 了 9 2 0 2面的冲击矿压防 治实践, 实现了该面的安全顺利 回采。 三河尖矿微震监测 系统 由中心站、 路新庄站 、 南风井站和煤气站等 4个观测站组成 , 其 中中心站 为三分量探头, 其余均为单分量 垂直方向 . 4 个 测站的埋深分别为 0 9 1 , 1 9 3 , 4 0 0 , 1 9 2 m, 各观测站 立体位置关系如图 1 。 各测站均以数字方式传输和 记录信号, 传输速率均 为 3 8 . 4 k b i t / s . 中心站按每 分量 l 0 0 0 s p s 的采样率接收测点数据。 监测网通 过 GP S定位实现微震定位 , 定位精度为数十米. 图 1 三河尖矿微震观测站 布置 m Fi g . 1 Th e a r r a n g e me n t o f o b s e r v a t i o n s t a t i o n o f mi c r o s e i s mi c s y s t e m i n S a n h e j i a n Mi n e 2 . 2 典型微震信号时一 频分析 由于震源 不 同, 各 种微震信 号表现 出不 同特 征. 9 2 0 2面开采过程 中, 微震监测系统接收到的微 震信号有列车通 过信号、 卸压爆破信号、 顶板活动 信号、 冲击矿压信号等 , 这些信号在频率特征 、 信号 持续时间 、 释放能量等多方面都各有特点. 本文 以 其中两个典型的顶板爆破 和冲击矿压信号为例进 行时一 频特征分析. 2 . 2 。 1 冲击矿压信号分析 坚硬顶板在高应力作用下 , 发生变形破坏和破 断, 释放大量弹性能 , 容易造成顶板型冲击矿压. 冲 击产生震动信号的最大不同处在 于其振 幅相 当大 震动能量大 . 在三河尖矿掘进期间 , 曾发生冲击 矿压 , 图 2 , 3即为一次 冲击 时部 分 台站 的记 录信 号. 该信号 接收时 间为 2 0 0 3年 2月 7日 7时 1 2 分 , 冲击发生地点在掘进 9 2 0 2面材料道. 早在 2月 4日微震监测系统就监测到该 区域震动频率数据 呈现异常 , 震动频次由之前 的 9 0次/ d增加为 l l 8 次/d , 并在 2月 5日继续上升达到 1 2 6次/ d后呈 下降趋势. 此外 , 钻屑法检测动力现象明显, 最大煤 维普资讯 1 4 8 采矿与安全工程学报 第 2 4卷 粉量 由 2月 4日的 0 . 9 5 k g / m 逐渐 上 升 至 2月 6 日的 1 . 5 k g / m. 蠢 馨 图 2 典型 冲击矿压微震信号振 幅一 时间曲线 Fi g. 2 The a mpl i t ude - t i me c ur v e s of t he t y pi c a l s i g na l of r oc k bur s t 1 . 4 0 1 . 0 0 0 . 6 0 0 . 2 0 0 .2 0 0 . 1 6 0 . 1 2 0 . 0 8 0 . 0 4 b 路新庄 lL . . . . 0 l 0 0 2 0 0 3 0 0 4 0 0 5 0 0 频率, H z 图 3 典 型冲击矿压微震信号幅频 曲线 Fi g. 3 The a mpl i t u de - f r e qu e nc y c ur v e s of t he t yp i c a l s i gn a l o f r o c k bur s t 信号产生时未进行爆破作业 , 信号持续时间大 约 1 1 0 0 ms , 衰减较快, 尾波较发育, 振幅相当大 , 煤气站观测显示震幅为 2 8 0 0 0 mV, 其余 3个观测 站显示震幅均为 1 0 0 0 0 mV, 是所有观测到的震动 信号中幅值最大的 , 表明冲击使得煤岩体内聚集的 弹性能大量释放. 信号波成分单一 , 频带分布较宽 , 从 0 ~1 0 0 Hz 均有分布. 2 . 2 . 2 顶板爆破信号分析 顶板爆破 强制放顶 以人为地破坏采空 区顶 板为直接 目的, 破坏顶板 的完整性 , 释放因压力聚 集 而产生 的 大量 弹性 能 , 减 少顶 板 因采 动影 响 和时 间效应而产生断裂时对煤层和支架的冲击震动. 图 4 , 5为 2 0 0 3年 6月 1 6日强制放顶产生震 动信号. 该信号强度较大 , 路新庄 站和中心站观测 到震 动幅值 均 达到 5 0 0 0 mV. 震动 持续 时 间短 , 大 约 1 0 0 0 ms , 衰减较快 , 尾波不发育. 震动频率主要 分 布在 0 ~1 5 0 Hz 之 间. 蠢 馨 0 4 0 0 8 0 0 l 2 0 0 1 6 0 0 2 0 0 0 时间/ ms 图 4 典 型顶板爆破微震信号振幅一 时间 曲线 Fi g .4 Th e ampl i t ude - t i me c u r v e s o f t he t yp i c a l s i gn a l o f f or c e d c ol l a p s e we i g ht i n g 奏 馨 0 l 0 0 2 o o 3 0 0 4 O O 5 0 0 频率/ I { z 图 5 典 型顶 板爆 破微 震信号幅频 曲线 Fi g. 5 The a mpl i t ude - f r e q ue nc y c ur v e s o f t h e t ypi c a l s i gna l o f f or c e d c ol l a ps e we i ght i ng 0 们 加 0 们 ∞ 0加 ∞ 锄 0加 仰 如 加 m 0 m 加 如m 0 m J 雠馨 如 0 如 如 m 0 m 0 m 维普资讯 第 2 期 曹安业等 高应力 区微震 监测信号特征分析 1 4 9 由图 4振幅一 时间曲线 可看 出, 该信号 可明显 分为两部分 , 一个较 小信号后 面紧跟一个 较大信 号; 从幅频 曲线也可看出, 频率出现两次高峰 明显 区别于冲击信号 , 表 明信号 由两种不同特性 的波 构成. 分析表明, 信号前部分为炸药爆炸 的震动信 号, 后部分为爆炸引起顶板破坏产生 的震动. 炸药 爆炸能和顶板的震动能叠加在一起 , 释放大量弹性 能, 产生较大震动. 2 . 3 不 同微震 信号 特征 比较 通过傅里叶变换, 分析各种微震信号的幅谱特 性 , 总结 出微地震 波形 的一些重要特征 , 进 而可较 快识别出接收的信号类型, 尤其是快速识别冲击危 险信号便于及时采取有效解危措施 , 并为研究微震 活动规律提供重要依据. 三河尖矿监测到的几种震 动信号波形有以下特点 1 各种信号持续 时间都较短 , 在数百 到千余 ms 之间 , 其 中顶板 来压信号持 续时间最短 , 只有 1 7 0 ms . 列车通过信 号为干扰信号 , 其信号持续时 间长达 1 mi n左右 , 可明显区分. 2 震动信 号大都衰减较快 , 只有 冲击余震信 号衰减较慢. 3 冲击 和顶板来压信号 的震动能量最大 , 信 号振 幅数量级为 1 O , 这是因冲击和来压时岩体内 变形破坏面积大 , 释放大量弹性能. 4 震动信号频率一般分布在 5 0 Hz以下 , 来压 信号频率分布在 1 0 0 Hz以下 , 强制放顶信号频率 在 0 2 0 0 Hz 之间都有分布, 冲击余震信号 由于衰 减较慢 , 高频成分较多 , 频率在 O ~5 0 0 Hz之间都 有分 布. 5 卸压爆破、 强制放顶信号 由于是膨胀源, 其 P波初始相位均向上 ; 冲击矿压信号 P波初动相位 也均 向上 , 表明该矿冲击震源主要产生向外的压缩 波, 造成巷道的张性拉伸破坏. 6 大部分微震信号波成分较简单 , 只有煤体 卸压爆破和顶板爆破的波成分较复杂 , 表 明爆破起 了诱发震动的作用 , 体现了卸压效果. 7 震动波传播路径不 同决定其复杂性质 , 由 于工作面震 动波在煤岩层 、 采空 区、 铁轨 中传播, 波 的相互干扰较大 , 给准确定位带来一定难度. 8 所有 的震动波形 的震源并非集 中于一点 , 而是一种呈立体状 的体震源. ’ 各震动信号的特征对 比见表 1 . 表 1 不 同微震信 号特征比较 Ta b l e 1 Th e c ompa r i s o n b e t we e n f e a t ur e s o f di f f e r e n t mi c r o s e i s mi c s i g n al s 3 结 论 1 冲击矿压微震监测技术具有区域性、 动态、 实时监测 的优点 , 实践证 明, 该技术用于预测和评 价矿区冲击矿压危 险性完全可行 , 具有广泛的发展 前 景. 2 由于震源不 同, 各种微震信号表现出不同 特征. 通过对三河尖矿 9 2 0 2面不同微震信号的时 一 频特征分析 , 总结出不同微震信号的各 自重要特 征 , 为现场快速识别冲击危险信号与及时采取解危 措施提供 了重要依据. 3 如何进一步完善不 同微震信号的特征 比较 指标 , 建立真正普遍化 、 实用可靠 的识别技术 , 并对 微震复杂信号进行识别 及 除噪, 使信 号识 别更快 速、 准确、 有效 , 仍需大量的统计和分析研究. 参考文献 [ 1 ] 窦林名, 何学秋.冲击矿压防治理论与技术[ M] .徐 州 中国矿 业大学出版社 , 2 0 0 1 . E 2 3 窦林名 , 赵 从 国 , 杨思 光 , 等. 煤 矿开 采 冲击 矿压灾 害 防治[ M3 .徐州 中国矿业 大学 出版社 , 2 0 0 6 . [ 3 1 郭然 , 潘 长 良. 有岩爆倾 向硬岩矿 床采矿 理论 与技 术E M3 .北 京 冶金工业 出版社 , 2 0 0 3 . E 4 3 林盛 , 刘 业新 , 李 衍达.基 于时 一频 分解技 术 的全 波列声波测井信号处理[ J 3 .清华大学学报 自然科 学版 , 1 9 9 7 , 3 7 3 6 3 6 6 . LI N S h e n g,LI U Ye - x i n,L I Ya n - d a . F u l l wa v e [ o r m s o n i c d a t a p r o c e s s i n g u s i n g t i me - f r e q u e n c y d e c o mp o s i t i o n [ J ] .J o u r n a l o f Ts i n g h u a Un i v e r s i t y S c i e n c e a n d Te c h n o l o g y ,1 9 9 7, 3 7 3 6 3 6 6 . 下转 第 1 5 4页 维普资讯 1 5 4 采 矿与安全工程学报 第 2 4卷 分析与控制对策[ J ] .能源技术与管理 , 2 0 0 6 2 4 - 5 . L I Gu i - c h e n,ZHANG No n g,KAN J i a - g u a n g,e t a 1 . Th e me c h a n i s m a n a l y s e s a n d c o n t r o l c o u n t e r me a s u r e a b o u t t w o c o a l l a n e w a y s r o o f f a l l a c c i d e n t s [ J ] .E n e r g y Te c h n o l o g y a n d Ma n a g e me n t ,2 0 0 6 2 4 - 5 . [ 4 ] 付跃升.简易锚杆桁 架在煤 巷支护 中 的应 用[ J ] .建 井技术 , 1 9 9 9 , 2 0 2 2 3 2 5 . F U Yu e - s h e n g . Ap p l i c a t i o n o f t h e s i m p l e b o l t t r u s s i n c o a l r o a d wa y s u p p o r t [ J ] . Mi n e C o n s t r u c t i o n Te c h n o l o g y ,1 9 9 9 , 2 0 2 2 3 2 5 . [ 5 ] 张益东.桁架锚 杆与 普通锚 杆对 顶板 的不 同支护作 用 口] _矿 山 压 力 与 顶 板 管 理 , 1 9 9 9 , 1 6 3 / 4 1 5 9 1 61 . ZHANG Yi do ng. The di f f e r e nt a c t i on s t o t he r oo f s u p p o r t o f t r u s s e s a n d c o mmo n b o l t s [ J ] .Gr o u n d P r e s s u r e a n d S t r a t a Co n t r o l ,1 9 9 9,1 6 3 / 4 1 5 9 1 61 . [ 6 ] 张农 , 高 明仕 , 许兴亮.煤巷预 拉力支护 体系及 其 工 程应 用口] .矿 山压 力与 顶板管理 , 2 0 0 2 , 1 9 4 1 4. ZHANG Non g, GAO M i n g - s h i , XU Xi ng - l i an g. Pr e t e nt i on e d s u pp or t i n g s y s t e m o f r oa dwa y a nd i t s e n g i n e e r i n g a p p l i c a t i o n[ J ] . Gr o u n d P r e s s u r e a n d S t r a t a Co n t r o l ,2 0 0 2 , 1 9 4 1 4 . [ 7 ] 张农 , 李 学华 , 高 明仕.迎采 动工作 面沿 空掘 巷预 拉 力支 护 及 工 程 应用 [ J ] .岩石 力学 与 工 程 学 报 , 2 0 0 4, 2 3 1 2 2 1 O O 一 2 1 O 5 . [ 8 ] [ 9 ] ZHANG No n g, LI Xu e - h u a , GAO M i n g s h i . P r e t e n s i o n e d s u p p o r t o f r o a d wa y d r i v e n a l o n g n e x t g o b a n d h e a d i n g a d j a c e n t a d v a n c i n g c o a l r a c e a n d i t s a p p l i c a t i o n [ J ] .C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d En g i n e e r i n g , 2 0 0 4, 2 3 1 2 2 1 0 0 2 1 0 5 . 张农 , 袁 亮 .离层 破碎 型煤 巷 顶板 的控制 原 理 [ J ] .采 矿与安全工程学报 , 2 0 0 6 , 2 3 1 3 4 3 8 . ZHANG Non g, YUAN Li a ng . Co nt r ol p r i nc i pl e o f s e p a r a t i n g a n d b r o k e n r o o f r o c k s t r a t a i n r o a d wa y [ J ] .J o u r n a l o f Mi n i n g S a f e t y E n g i n e e r i n g ,2 0 0 6 , 2 3 1 3 4 3 8 . 宁建 国, 蒋金 泉 , 董 建军 , 等. 回采 巷道顶板 桁架锚 杆 支护角 度 的优 化 研究 [ J ] .矿 山压 力 与顶 板 管 理 , 2 0 0 3, 2 0 1 2 0 2 1 . NI NG J i a n - g u o 。J I ANG J i n - q u a n,D0NG J i a n - j u n , e t a 1 . Th e o p t i mu m s t u d y o f t h e a n gl e o f t h e t r u s s b o l t s u p p o r t i n g t h e r o o f i n t h e a c t u a l mi n i n g r o a d wa y [ J ] .G r o u n d P r e s s u r e a n d S t r a t a C o n t r o l ,2 0 0 3 , 2 0 I 2 0 2 I . [ 1 O ] 李 桂 臣 , 王 凯 , 张农 , 等. 锚带 网索 支 护在 大倾 角、 松软 、 突出煤 巷 中 的应 用 [ J ] .华北 科 技学 院学 报 , 2 0 0 5 , 2 4 1 O 一 1 2 . LI Gu i c h e n ,W ANG Ka i ,Z HANG No n g,e t a 1 . Ap p l i c a t i o n o f a n c h o r b e l t n e t a n d c a b l e s u p p o r t i n l a r ge s l op e a n gl e s o f t an d g a s ou t bur s t c o a l r o a d d r i v i n g [ J ] .J o u r n a l o f N o r t h C h i n a I n s t i t u t e o f S c i e n c e a n d Te c h n o l o g y ,2 0 0 5, 2 4 1 0 1 2 . 上接 第 1 4 9页 [ 5 ] 李玉 , 黄梅 , 廖国华 , 等.冲击地压 发生前微震 活 [ 6 ] 动时空变化的分析特征 [ J ] .北京科技大学学报 , 1 9 9 5 , 1 7 1 1 0 1 3 . LI Yu,HU ANG M e i ,LI AO GUO hua,e t a1 .Fr a c t a l c ha r a c t e r of t he t i me - Spa c e c ha n ge o f mi c r os ei s mi c i t y b e f o r e c o a l o u t b u r s t s o r r o c k b u r s t s [ J ] .J o u r n a l o f U n i v e r s i t y o f S c i e n c e a n d Te c h n o l o g y Be i j i n g,1 9 9 5 , 1 7 1 1 0 - 1 3 . 陆菜平 , 窦林名 , 吴兴荣 , 等.岩体微震 监测 的频 谱分 析与信号识别[ J ] .岩土 工程学 报 , 2 0 0 5 , 2 7 7 7 7 2 一 [ 7 ] [ 8 ] n 5. L U Ca i - p i n g,D0U Li n - mi n g,W U Xi n g - r o n g,e t a 1 . Fr e que nc y s pe c t r um an a l y s i s o n mi c r o s e i s mi c moni t o r i n g a n d s i g n a l d i f f e r e n t i a t i o n o f r o c k ma t e r i a l [ J ] . Ch i n e s e J o u r n a l o f Ge o t e c h n i c a l En g i n e e r i n g,2 0 0 5, 2 7 7 7 7 2 7 7 5 . S L A w0MI R J G, AN DR Z E J K.矿 山地 震 学 引论 [ M] .北 京 地震 出版社 , 1 9 9 8 . 秦玉红.矿井 开采 引发 矿震规 律及其 应用研 究[ D ] . 徐州 中国矿业大 学能源与安全工程学院 , 2 0 0 5 . 维普资讯