Microsoft Word - 4曹安业h.doc

Size: px

Start display at page:

Download "Microsoft Word - 4曹安业h.doc"

给滕
7 years ago
Views:

1 第 32 卷第 1 期采矿与安全工程学报 Vol.32 No 年 01 月 Journl of Mining & Sfety Engineering Jn 文章编号 : (2015) 孤岛面开采强矿震异常区的被动声波探测技术及应用曹安业 1, 井广成 1, 窦林名 1, 巩思园 1, 赵永亮 2, 张贞良 (1. 中国矿业大学矿业工程学院, 深部煤炭资源开采教育部重点实验室, 煤炭资源与安全开采国家重点实验室, 江苏徐州 ;2. 兖州煤业股份有限公司鲍店煤矿, 山东邹城 ) 摘要针对鲍店煤矿 103 上 05 孤岛面综放开采面临的高应力集中与强矿震威胁, 论文在分析声波波速与加载载荷的关系基础上, 建立了基于矿震被动声波探测的孤岛面开采矿震异常区评价方法, 提出了采用 P 波波速异常与波速梯度异常 2 个指标进行采动应力集中程度的评价, 并给出了判别标准, 最后开展了现场实践研究表明 :P 波波速与应力场具有正相关特性,103 上 05 孤岛工作面临近小采空区侧常出现高波速与高波速梯度异常区, 与后续开采阶段的强矿震异常区基本重合实践证明该方法能够对孤岛面开采强矿震危险性与危险区域做出有效动态评价, 强矿震活动的预警准确率较高矿震被动声波探测技术提升了煤矿强矿震冲击矿压等煤岩动力灾害的微震识别水平关键词孤岛工作面 ; 矿震异常 ; 波速 ; 声波探测 ; 波速梯度中图分类号 TD 324 文献标志码 A DOI /j.cnki.jmse Seismic hzrd ssessment in complex islnd colfce by computed tomogrphy CAO Anye 1,JING Gungcheng 1,DOU Linming 1,GONG Siyun 1, ZHAO Yongling 2,ZHANG Zhenling 2 (1.School of Mines,Key Lbortory of Deep Col Resource Mining,Ministry of Eduction of Chin, Stte Key Lbortory of Col Resources nd Sfe Mining,Chin University of Mining & Technology, Xuzhou,Jingsu ,Chin;2.Bodin Col Mine,Ynzhou Col Mining Compny Limited, Zoucheng,Shndong ,Chin) Abstrct According to the high mining stress nd seismic hzrd risk in excvtion opertion of Islnd Colfce in Bodin col mine, the reltion between P-wve velocity nd stress ws nlyzed, nd the seismic hzrd ssessment method by using pssive computed tomogrphy in Islnd Colfce ws estblished. Menwhile, the indexes of P-wve velocity nomly nd velocity grdient nomly were proposed to evlute the stress concentrtion or seismic hzrd zones, nd the discrimintion stndrds were lso constructed. Finlly, the ssessment method ws pplied in Is- 收稿日期 : 基金项目 : 国家自然科学基金青年科学基金项目 ( ); 江苏高校优势学科建设工程项目 (SZBF B35); 中央高校基本科研业务专项资金项目 (2012BZB17,2010QNA30); 国家十二五科技支撑计划项目 (2012BAK04B06); 国家重点基础研究发展计划 (973) 项目 (2010CB226805) 作者简介 : 曹安业 (1982 ), 男, 江苏省盐城市人, 副教授, 博士生导师, 从事矿山压力冲击矿压采矿地球物理等方面的研究 E-mil:conye@163.com Tel:

2 第 1 期曹安业等 : 孤岛面开采强矿震异常区的被动声波探测技术及应用 21 lnd Colfce. The results show tht P-wve velocity nd stress presents chrcteristics of positive correltion. The high velocity nomlies nd velocity grdient nomlies re lwys locted in the regions djcent to the nrrow gofing side of Islnd Colfce, which re bsiclly similr to the nomly res of seismic hzrds. The tomogrphy detection prctices show tht it is fesible for seismic hzrd ssessment nd nomly region division during islnd colfce mining, nd the hzrd ssessment ccurcy is high. Thus, the pssive computed tomogrphy cn improve the monitoring nd pre-wrning level for seismic hzrds in underground mining. Key words islnd colfce;seismic hzrd;wve velocity;tomogrphy detection;velocity grdient 矿震冲击矿压等动力现象是煤岩变形破裂动态演化导致的突发性动力显现, 对其诱发动态过程的有效监测预警是防止灾害发生的关键, 是迫切需要解决的安全难题之一目前, 冲击矿压的监测方法主要包括微震法声发射电磁辐射应力监测钻屑量等 [1-4] 其中, 钻孔应力法钻屑量等主要监测诱发冲击危险的煤体静载荷, 微震法声发射等主要监测诱发冲击危险的动载荷微震监测法是利用煤岩体受力变形破坏后激发的震动波来进行监测工程岩体稳定性的技术方法相比于其他监测手段, 该方法具有区域动态三维实时监测的特点, 可为评价全矿范围内的强矿震冲击矿压等危险提供依据 [1,5] 目前, 微震监测预警主要是根据震动频次能量等微震活动的时序变化或微震空间演化进行冲击矿压矿震危险的定性评价近年, 声波探测技术作为一种新的地球物理方法被逐渐应用于矿山工程与地质诊断 [6], 特别是将声波探测技术和微震实时监测相结合形成的矿震被动声波探测方法, 是冲击矿压或强矿震危险评价和预测的最新发展方向基于矿山微震监测, 采用矿震被动声波探测方法的监测范围可达 2000 m 2000 m 左右 [7] 该方法已在国外一些高冲击或强矿震矿井开展了初步研究与应用 [7-12] 例如,Lurk [7] 利用被动声波探测技术对波兰 Zbrze Bielszowice 等煤矿矿震异常区进行评价, 发现高能量矿震往往出现在高波速区 [8] Kry Luxbcher 等基于微震监测, 对美国西部某煤矿长壁工作面开采过程中的矿震波速场进行了反演, 发现震动波速与采动应力分布具有较好的吻 [9] 合性 Friedel 等研究了美国某深部金矿开采过程中的波速分布区域, 得出低波速区往往分布于巷道采场等空区附近, 而高波速主要分布于矿柱中 [10] 心等高岩爆危险区域 Scott 等研究了美国华盛顿州 Homestke 矿 m 深部矿柱内的波速分布特征, 得出低波速区主要分布于钻孔巷道及相邻矿柱边缘, 而高波速区主要分布在高应力矿柱内因此, 矿震被动声波探测技术用于探测采动应力分布或评价冲击矿压岩爆或强矿震危险是可行的, 该方法可初步实现冲击矿压或强矿震危险的微震监测由定性评价向定量预警转变煤矿开采过程中时常会出现孤岛工作面开采孤岛工作面由于多侧采空, 尤其是在坚硬巨厚顶板条件下, 覆岩断裂范围和运动剧烈程度较大, 开采诱发强矿震冲击矿压的危险性也更高因此, 孤岛工作面的高冲击或强矿震危险力源主要以覆岩破断的动载扰动为主微震监测应作为孤岛面开采矿震异常区探测或动压危险评价的主要手段本文以兖州鲍店煤矿 103 上 05 孤岛综放面开采为研究对象, 基于矿山微震监测, 建立孤岛工作面开采的被动声波探测模型与预警准则, 动态评价该孤岛工作面开采过程中矿震异常区的分布与演化过程, 为孤岛工作面强矿震冲击矿压的有效预警提供依据 1 工程背景矿震活动对鲍店煤矿的安全生产造成很大威胁, 该矿二采区 2305,2310 等工作面开采过程中曾发生过强矿震摧毁密闭事故, 而最大矿震发生于十采区 103 上 04 面开采过程中, 达到震级 3.7 级 103 上 05 工作面位于十采区南部, 地面平均标高 m, 井下平均标高 -428 m, 基本顶为 m 厚的粗砂岩, 煤层 130 m 上方赋存有厚达 200 m 的巨厚粗砂岩工作面南部与 103 上 06 采空区毗邻, 北部与 04,03,02,01 采空区相邻, 东部为井田边界, 西部与二采区采空区相邻 103 上 05 工作面右侧区域于 2004 年已回采 945 m, 后续左侧回采阶段将形成四面采空的孤岛结构同时, 放眼整个矿区, 十采区周边的二采区四采区也已经基本回采完毕, 因此,103 上 05 工作面又是一个矿井较大范

3 22 采矿与安全工程学报第 32 卷围内的大孤岛结构, 见图上 05 工作面这样一个矿井大孤岛与采区小孤岛双重影响下的特殊孤岛结构, 其覆岩运动和矿压显现特征与普通工作面开采有较大差别, 工作面附近应力集中程度较高 ( 见图 2), 矿震活动的动态扰动可能更为强烈图上 05 工作面布置 Fig.1 Lyout of LW 图上 05 孤岛工作面垂直应力模拟云图 Fig.2 Simultion results of verticl stress round LW 图 3 为 103 上 05 工作面内 99-1 钻孔综合柱状图, 根据关键层理论, 可判断工作面 m 厚的基本顶为低位亚关键层 1, 距煤层 106 m 厚约 m 的粉砂岩层为高位亚关键层 2, 距煤层约 130 m 总厚达 206 m 的巨厚红层为高位主关键层 103 上 05 工作面北侧采空区宽度范围 ( 约 750 m) 远大于其开采深度 ( 约 450 m), 顶板垮落充分, 属充分采动, 南侧采空区宽度 ( 约 200 m) 相比较小, 开采尚未能引起高位主关键层的大范围破断运动因此,103 上 05 工作面属一侧充分采动一侧非充分采动的非对称孤岛工作面, 其回采前上覆岩层结构示意如图 4 所示 103 上 05 工作面回采后,05, 06 工作面采空区覆岩易形成一个较大结构拱, 引起上覆巨厚岩层的整体性破断运动图钻孔综合柱状图 Fig.3 Synthesis column of borehole 99-1 图上 05 工作面回采前覆岩结构示意图 Fig.4 Sketches of overlying strt structure bove LW 孤岛面开采矿震异常区的被动声波探测技术 2.1 煤岩应力状态与 P 波波速的关系研究表明, 矿震声波波速与煤岩类型破裂面角度各向异性应力场孔隙率弹性模量等密切相关因此, 采动应力场分布与动态迁移可通过测量声波波速 ( 主要为 P 波 ) 来反映 [11,13-19] Adms [13] 等在 20 世纪 20 年代就发现被动声波探测是建立在煤岩体弹性模量与波速密切相关的基础之上的 [14] Yle 等研究发现岩石随应力增加孔隙率会降低, [15] 从而导致波速增加 Nur 等研究了 P,S 波波速与施加载荷方向的关系, 得出 P,S 波均随应力增加而增大, 且增加幅度与应力方向及 P 波传播方向相关, 即当波的传播方向与载荷方向平行时波速变化最大此外,Jones,Toksoz [16-17] 等分别研究了孔

4 第 1 期曹安业等 : 孤岛面开采强矿震异常区的被动声波探测技术及应用 23 隙压力孔隙流体等对波速的影响煤矿冲击矿压强矿震等煤岩动力现象的发生都与应力有着密切关系, 矿震被动声波探测就是基于矿山微震监测, 通过获得反演区域内矿震活动的震动波波速分布, 评价矿震异常区或冲击危险性对十采区顶板岩样进行单轴压缩与单轴循环加卸载下的声波波速测试 [18] 图 5 为岩样加载过程中对应的应力与 P 波波速的拟合曲线, 由图可见, 轴向应力与 P 波波速间具有较明显的幂函数关系其中, 单轴加载下 P 波波速在载荷初始阶段有较高变化梯度, 随着载荷不断增加, 波速上升幅度减缓并逐渐趋于稳定 ; 单轴循环加卸载过程中 P 波波速随着应力变化产生相应变化, 且加卸载在高应力变化时波速变化慢, 而在低应力时波速变化快 [18] 波速 /(m s -1 ) 波速 /(m s -1 ) 测量波速试验公式 /MP () 单轴压缩条件下应力与 P 波波速关系 3200 纵波波速应力时间 /s (b) 单轴循环加卸载条件下应力与 P 波波速关系图 5 单轴加载条件下应力与 P 波波速之间的相关性曲线 Fig.5 Reltion between stress nd P-Wve velocity in unixil compression test 2.2 矿震被动声波探测的技术原理利用震动波速研究开采区域煤岩体应力状态等信息可由主动源与矿震 2 种激发源得到由爆破或重锤敲击等主动源激发震动形成的成像技术称为主动声波探测技术 ; 而基于开采活动诱发的矿山震动 ( 未知其初始时刻与位置 ) 形成的成像技术称为被动声波探测技术主动声波探测技术必须人工激发震源和携带仪器下井, 劳动量大, 不利于危险区域的长期探测研究被动声波探测技术因矿震事件的能量要比人工激发源高很多, 研究范围较大, 且由于矿震活动一直伴随开采发生, 可利用矿震震动波进行长期动态探测 1/MP 根据微震监测台站与矿震震源之间的距离 l 和台站接收到的 P 波初至传播时间 t 来反演区域内的波速分布 v(x,y,z) 或慢度 s(x,y,z)=1/v(x,y,z) 假设第 i 个 P 波的传播路径为 l i, 其传播时间为 t i, 具体关系可表示为 [11] : l v = vt = l (1) t dl ti = = s( x, y, z) dl li v x, y, z (2) li ( ) m ( 1,, ) (3) t = d s i = n i ij j j = 1 式中 :t i 为 P 波传播时间,s;l i 为第 i 条 P 波射线路径 ;v(x,y,z) 为 P 波传播速度,m/s;s(x,y,z)=1/v(x,y,z) 为慢度,s/m;d ij 为第 i 条 P 波射线穿过第 j 个网格的长度 ;n 为射线总数 ;m 为网格数量上式 (3) 写成矩阵形式如下 : DS = T (4) 式中 :D=(d ij ) m n 为距离矩阵 ;T=(t i ) m 1 为传播时间向量 ;S=(S j ) m 1 为慢度列向量通过求解该矩阵就可以得到离散慢度分布, 从而实现开采区域的声波速度场反演成像实际应用中要反复求解 (4) 式来得到重建区域的速度场, 通过速度场的异常反演并判断强矿震或冲击矿压危险区域矿震震动波速的反演示意图如图 6 所示图 6 声波波速反演示意图 Fig.6 Sketch mp of seismic wve velocity tomogrphy 由于矿震震源位置未知, 导致计算的震源位置与实际位置相比会有一定偏差, 对波速反演结果会有一定影响当然, 一般矿震都集中在当前开采区域, 尤其是高冲击或强矿震危险区域微震台网往往布设较密, 能够形成较好的射线覆盖范围, 反演结果较可靠矿震被动声波反演时, 需事先根据常值速度模型估计震源位置, 并作为模型输入的初始值 ; 求解过程中, 每改变一次波速模型, 需重新进行震源位置计算, 并修改传播射线起点位置 ; 按上述方法反复迭代, 直到残差小于给定门槛值或达到最大迭代次数本文采用 SIRT 迭代反演算法 [12]

5 24 采矿与安全工程学报第 32 卷 2.3 矿震异常区的被动声波探测模型矿震异常区预警评价的基础是确定煤岩层中应力状态与应力集中程度研究表明, 高应力或应力集中区域, 相对其他区域将出现 P 波波速的正异常, 其异常值由下式计算得到 : vp vp An = (5) vp 式中 :v p 为区域内某点 P 波波速值 ; v p 为模型波速平均值煤层开采形成的采空区导致覆岩载荷加载到相邻支撑区域, 形成一侧应力降低与一侧高应力集中区, 采空区对应低波速区, 应力集中区域则对应高波速区 2 个区域之间是从高应力向低应力过渡区域, 即为高应力差区域矿震异常区不仅发生在高波速区域, 也易发生在波速梯度变化区域 [11,18] 因此,103 上 05 工作面开采前覆岩内波速分布示意可见图 7, 尤其是工作面小采空区侧由于主关键层尚未断裂, 存在高应力 ( 差 ) 与高波速 ( 梯度 ) 区域表 2 vg 异常变化与矿震异常之间的关系 Tble 2 Reltion between vg vlue nd seismic hzrd 矿震异常指标异常对应的动压危险性 vg 异常 /% 0 无 <5 1 弱中等强 >25 3 孤岛面开采矿震异常区的被动声波探测实践 3.1 矿震被动声波探测方案鲍店矿已安装波兰 16 通道 SOS 微震监测系统, 而能够对 103 上 05 工作面进行有效监测的台站约有 9 个, 如图 8 实心圆所示图上 05 工作面附近微震监测台站布置 Fig.8 Lyout of geophone sttions round LW10305 图上 05 工作面覆岩中波速分布示意图 Fig.7 Sketch of velocity distribution in strt structure of LW10305 对波速梯度变化, 可采用 vg 值描述相邻节点间波速的变化程度, 其异常值可表达为 vg vg A = (6) n vg 式中 vg 为模型中构建的 vg 平均值表 1,2 分别为波速正异常波速梯度异常与矿震危险之间的关系 [18,20] 表 1 波速正异常变化与应力集中程度关系表 Tble 1 Reltion between positive nomlies of P-wve speed nd stress concentrtion 矿震异常指标应力集中特征波速正异常 /% 0 无 <5 1 弱 5~15 2 中等 15~25 3 强 >25 矿震被动声波探测选取 103 上 05 工作面于 2012 年 3 月 2012 年 10 月开采期间的震动波形作为反演研究数据所有 P 波到时都由人工进行标定, 由此确定震源分布对形成的射线进行波速统计, 统计得到每个波速区间内的射线条数由于地面探头的存在, 建立层状波速模型进行计算, 并选用射线弯曲算法中的混合算法网格划分为 ,X,Y,Z 方向间距为 30 m 30 m 125 m, 模型波速在 2~6 km/s 范围内梯度分布为探测下一阶段的强矿震异常区, 采用对本时段 ( 以 1 个月为 1 个反演周期 ) 的 P 波波速反演结果进行规律分析, 以达到对下一开采阶段强矿震危险的预警目的 3.2 矿震异常区的被动声波探测实践选取 103 上 05 工作面煤层 -420 m 水平切片的波速值和波速梯度变化系数 vg 云图作为工作面矿震异常区的探测评价 103 上 05 工作面于 2012 年 3 月 1 日至 3 月 31 日初采阶段监测的有效震动 644 个 ( 见图 9), 选择 105 个适合进行波速反演的波形进行声波探测, 射线条数 541, 每个波形的射线个数均大于 4 依据 3 月份微震活动进行的 P 波波速反演结果如图 9b,c

6 第 1 期曹安业等 : 孤岛面开采强矿震异常区的被动声波探测技术及应用上 04 左 103 上 05 左 103 上 06 左能量 /J 0~ ~ ~ ~10 7 属于反演不可信区域因此, 随着工作面推进, 下一阶段矿震异常区仍主要处于工作面前方 05,06 面交汇处, 但异常区明显扩大能量 /J 0~ ~ ~ ~10 7 y 方向 /m y 方向 /m () 矿震活动平面分布 (2012 年 3 月 ) (b) 波速场反演 (c) vg 异常指数反演波速 /(km s -1 ) 图上 05 工作面 3 月份 P 波波速反演结果 Fig.9 P-wve velocity tomogrphy results in LW10305 ccording to seismic events occurred in Mr., 2012 由图 9 可见, 高波速及高 vg 异常区主要分布在 103 上 05 工作面与 103 上 06 采空区交接处, 处于反演的可信区域之内, 即工作面前方小采空区侧处于矿震异常区, 强矿震危险较高图 10 为该工作面 4 月份发生的能量 >10 3 J 矿震分布情况, 可见下一回采阶段出现的高能量矿震主要分布在反演的高波速与高 vg 异常区域, 佐证了 3 月份矿震声波反演结果的准确性图 11 为工作面 4 月份微震监测数据的 P 波波速反演结果由图可见, 高波速及高 vg 异常区主要分布于 103 上 05 工作面与 103 上 06 采空区交接处采区煤柱及大采空区局部区域结合矿震射线平面投影 ( 图 12),P 波射线切割密集区处于 05,06 面交接处, 属于反演可信区域, 其他区域射线分布较少, 异常值 103 上 04 左 103 上 05 左 103 上 06 左图上 05 工作面 4 月较大能量矿震分布 (>10 3 J) Fig.10 Strong tremors with energies lrger thn 10 3 J occurred in LW10305 during Apr., 2012 () 波速场反演 (b) vg 异常指数反演图上 05 工作面 4 月份 P 波波速反演结果 Fig.11 P-wve velocity tomogrphy results in LW10305 ccording to seismic events occurred in Apr., 2012 图上 05 工作面 4 月份矿震射线平面投影 Fig.12 P-wve ry pth of events occurred in Apr., 2012

7 26 采矿与安全工程学报第 32 卷 103 上 05 工作面 5 月份出现的能量 >10 3 J 强矿震主要分布在反演的 05 工作面前方轨道顺槽侧及 06 采空区内, 处于反演的可信区域内 ( 图 13) 声波被动探测能够对 103 上 05 孤岛工作面开采过程中的矿震异常区及动压危险进行有效动态评价图 14 为 103 上 05 工作面 2012 年 9 月末采阶段微震监测数据的 P 波波速成像结果采区煤柱附近, 属于反演可信区域, 其他区域射线分布较少, 属于反演盲区后续开采阶段发生的能量 >10 3 J 强矿震活动也大多处于反演的矿震异常区域内 ( 图 16), 动态评价效果较好图上 05 工作面 9 月份矿震射线平面投影 Fig.15 P-wve ry pth of events occurred in Sep., 2012 图上 05 工作面 5 月较大能量矿震分布 (>10 3 J) Fig.13 Strong tremors with energies lrger thn 10 3 J occurred in LW10305 during My, 2012 () 波速场反演 (b) vg 异常指数反演图上 05 工作面 9 月份 P 波波速反演结果 Fig.14 P-wve velocity tomogrphy results in LW10305 ccording to seismic events occurred in Sep., 2012 由图可见, 随着工作面逐渐临近停采线位置, 高波速及高 vg 异常区主要分布于在 05,06,04 采空区贯通区域采区边界煤柱区域及 01,02 采空区内结合矿震射线平面投影 ( 图 15),P 波射线切割密集区处于 05 工作面与 06,04 采空区交接处及 103 上 04 左 103 上 05 左 103 上 06 左能量 /J 0~ ~ ~ ~10 7 图上 05 工作面 10 月较大能量矿震分布 (>10 3 J) Fig.16 Strong tremors with energies lrger thn 10 3 J occurred in LW10305 during Oct., 结论 1) 对于同一性质的煤岩体,P 波波速与应力场或应力差具有正相关特性 103 上 05 孤岛工作面小采空区侧由于主关键层尚未断裂, 存在高应力 ( 差 ) 与高波速 ( 梯度 ) 区域, 强矿震活动较异常 2) 根据 103 上 05 工作面 P 波波速反演实践, 由波速异常和波速梯度异常系数等组成的被动声波探测技术可对孤高工作面开采的矿震异常区做出有效动态评价, 预警准确率较高 3) 被动声波探测技术的评价效果取决于波速模型的建立迭代算法的选择计算数据的筛选及反演可信区的识别等, 该技术的发展与应用有效提升了强矿震冲击矿压等动力灾害的微震识别水平参考文献 : [1] 姜福兴, 杨淑华, 成云海, 等. 煤矿冲击地压的微地震监测研究 [J]. 地球物理学报,2006,49(5):

8 第 1 期曹安业等 : 孤岛面开采强矿震异常区的被动声波探测技术及应用 27 JIANG Fuxing,YANG Shuhu,CHENG Yunhi,et l.a study on microseismic monitoring of rock burst in col mine[j].chinese Journl Geophysics,2006,49(5): [2] 齐庆新, 李首滨, 王淑坤. 地音监测技术及其在矿压监测中的应用研究 [J]. 煤炭学报,1994,19(3): QI Qingxin,LI Shoubin,WANG Shukun.Appliction of AE technique in monitoring ground pressure[j].journl of Chin Col Society,1994,19(3): [3] 曲效成, 姜福兴, 于正兴, 等. 基于当量钻屑法的冲击地压监测预警技术研究及应用 [J]. 岩石力学与工程学报,2011,30(11): QU Xiocheng,JIANG Fuxing,YU Zhengxing,et l. Rockburst monitoring nd precution technology bsed on equivlent drilling reserch nd its ppliction[j].chinese Journl of Rock Mechnics nd Engineering,2011,30(11): [4] HE Xueqiu,CHEN Wenxue,NIE Bisheng,et l. Electromgnetic emission theory nd its ppliction to dynmic phenomen in col-rock[j].interntionl Journl of Rock Mechnics nd Mining Sciences,2011,48(8): [5] 曹安业, 窦林名, 秦玉红, 等. 高应力区微震监测信号特征分析 [J]. 采矿与安全工程学报,2007,24(2): CAO Anye,DOU Linming,QIN Yuhong,et l.anlysis on chrcter of microseismic monitoring signl in high stress zone[j].journl of Mining & Sfety Engineering, 2007,24(2): [6] ZHAO Yongguo,LI Qin,GUO Hong,et l.seismic ttenution tomogrphy in frequency domin nd its ppliction to engineering[j].science in Chin:Series D, 2000,43(4): [7] LURKA A.Loction of high seismic ctivity zones nd seismic hzrd ssessment in Zbrze Bielszowice col mine using pssive tomogrphy[j].journl of Chin University Mining & Technology,2008,18(2): [8] LUXBACHER K,WESTMAN E,SWANSON P,et l. Three-dimensionl time-lpse velocity tomogrphy of n underground longwll pnel[j].interntionl Journl of Rock Mechnics nd Mining Sciences,2008,45(4): [9] FRIEDEL M J,SCOTT D F,JACKSON M J,et l.3-d tomogrphic imging of nomlous stress conditions in deep US gold mine[j].journl of Applied Geophysics, 1996,36(1):1-17. [10] SCOTT D F,WILLIAMS T J,FRIEDEL M J,et l. Reltive stress conditions in n underground pillr, Homestke Mine,Led,SD[J].Interntionl Journl of Rock Mechnics nd Mining Sciences,1997,34: [11] GIBOWICZ S J,KIJKO A.An introduction to mining seismology[m].sn Diego:Acdemic Press,1994. [12] HOSSEINI N,ORAEE K,SHAHRIAR K,et l.pssive seismic velocity tomogrphy on longwll mining pnel bsed on simultneous itertive reconstructive technique (SIRT)[J].Journl of Centrl South University,2013, 19(8): [13] ADAMS L,WILLIAMSON E.On the compressibility of minerls nd rocks t high pressures[j].journl of the Frnklin Institute,1923,195: [14] YALE D.Recent dvnces in rock physics[j].geophysics,1985,50(12): [15] NUR A,SIMMONS G.Stress-induced velocity nisotropy in rock:n experimentl study[j].journl of Geophysicl Reserch,1969,74(27): [16] JONES L,WANG H.Ultrsonic velocities in Cretceous shles from the Williston bsin[j].geophysics,1981, 46(3): [17] TOKSÖZ M,CHENG C,TIMUR A.Velocities of seismic wves in porous rocks[j].geophysics,1976, 41(4): [18] 巩思园. 矿震震动波波速层析成像原理及其预测煤矿冲击危险应用实践 [D]. 徐州 : 中国矿业大学,2010. [19] 曹安业. 采动煤岩冲击破裂的震动效应及其应用研究 [D]. 徐州 : 中国矿业大学,2009. [20] 窦林名, 何学秋. 采矿地球物理学 [M]. 北京 : 中国科学文化出版社,2002:23-26.

Microsoft Word - 10朱广安h.doc

Microsoft Word - 10朱广安h.doc 第 33 卷第 4 期采矿与安全工程学报 Vol.33 No.4 216 年 7 月 Journal of Mining & Safety Engineering July 216 文章编号 :1673-3363-(216)4-63-6 1 深埋复杂不规则孤岛工作面冲击矿压机制研究朱广安, 窦林名, 刘阳, 苏振国, 李慧 ( 中国矿业大学煤炭资源与安全开采国家重点实验室, 矿业工程学院, 江苏徐州