开滦矿区深部开采中巷道围岩稳定性研究.pdf
第 2 4卷第 1 1 期 2 0 0 5年 6月 岩石力学 与工程 学报 C h i n e s e J o u r n a l o fR o c k Me c h a n i c s a n d E n g i n e e r i n....g Vo 1 . 2 4 No . 1 1 J u n e , 2 o 0 5 开滦矿 区深部开采中巷道围岩稳定性研究 姜耀 东 ,刘文岗 ,赵毅鑫 ,殷作如 。 , 1 .中国矿业大学 北京校区力学与建筑工程学院 ,北京 李建民 2 ,邓智毅 2 ,洪溢清 1 0 0 0 8 3 2 .开滦矿业集团,河北 唐 I J 1 0 6 3 0 0 0 摘要 在对开滦矿区赵各庄矿、唐 山矿深部开采过程中巷道 围岩变形、破坏特征和矿井动力显现观测的基础上 , 结合数值模拟研究发现 1 开滦矿区在千米深开采时水平构造应力远大于 自重应力,原岩应力值远大于巷道围岩 强度值,巷道将不可避免地要遭受破坏并表现 为大变形强流变和严重底鼓; 2 随着开采深度增大 ,软弱煤岩层发 生冲击地 压的危险性随之增大 ,通常 以煤层的整体突出为特征 , 其 发生机理 和高应力状态下硬岩 煤 中的岩爆 有区别; 3 原岩应力状态对巷道围岩的变形破坏有着重要影响,而巷道形状对巷道围岩的破坏形式影响不大 。 关键词岩石力学;深部开采;原岩应力;冲击地压;巷道大变形;巷道底鼓 中图分类号T D 3 2 2 . 4 文献标识码 A 文章编号1 0 0 0 6 9 1 5 2 0 0 5 l 1 1 8 5 7 0 6 S TUDY oN S URRoUNDI NG RoCK S TABI LI TY oF DEEP M I NI NG I N KAI LUAN M I NI NG GRoUP J I ANG Ya o d o n g l , LI U W e n . g a n g , ZHAO Yi x i n , YI N Z u o r u ’ 2 , LI J i a n mi n 、 2 , DENG Zh i y i . f 1 .HONG Yi . q i n g t 2 r 1 . S c h o o l o f Me c h a n i c s ,A r c h i t e c t u r e a n d Ci v i l E n g i n e e r i n g ,C h i n a U n i v e r s i t y o fMi n i n g a n d T e c h n o l o g y ,B e ij i n g 1 0 0 0 8 3 ,Ch i n a 2 . K a i l u a n Mi n i n g Gr o u p C o . ,T a n g s h a n 0 6 3 0 0 0 ,C h i n a J Abs t r a c t Ba s e d o n t h e o bs e r v a t i o n a n d n u me ric a l s i mu l a t i o n o f r o a d wa y f a i l u r e i n d e e p mi n i n g , i t i S f o u n d t h a t t h e h o r i z o n t a l t e c t o n i c s t r e s s i S mu c h g r e a t e r t h a n t h e o rig i n a l v e r t i c a l r o c k s t r e s s a n d t h e i n s i t u s t r e s s i S b e y o n d the s e n g th o f r o c k s u r r o u n d i n g r o a d wa y i n Ka i l u a n mi n i n g a r e a . Th u s , t h e l arg e d e f o r ma t i o n o f r o a d wa y a n d c a t a s t r o p h i c fl o o r h e a v e are d e e me d t o b e i n e v i t a b l e . W i t h t h e i n c r e a s i n g o f mi n i n g d e p t h, c o a l b u mp s b e c o me mo r e s e v e r e an d v i o l e n t . T r a n s l a t o r y c o a l b u mp s are t h e ma i n d y n a mi c h a z ard s i n d e e p mi n i n g o f Ka i l u a n mi n i n g g r o u p, wh i c h h a v e d i f f e r e n t me c h a n i s m c o mp are d t o the r o c k b u r s t s h a p p e n e d i n h ard r o c k u n d e r l arg e s c a l e i n . s i t u s t r e s s . F u r t h e rm o r e , i t i s n o t i c e d t h a t the fra c t u r e z o n e o f s u rro u n d i n g r o a d wa y i s ma i n l y i n flu e n c e d b y t h e i n . s i t u s t r e s s . Ho we v e r , the r o a d wa y s h a p e h a s l i t t l e e f f e c t s o n t h e fra c t u r e z o n e o f s u rro u n d i n g r o a d wa y . Ke y wo r d s r o c k me c h a n i c s d e e p mi n i n g; g r o u n d s t r e s s ; c o a l b u mp s ;l arg e d e f o rm a t i o n o f r o a d wa y ;flo o r h e a v e 1 引 言 我国的煤炭资源埋深在 1 0 0 0 m 以下的为 2 .9 5 万亿 t ,占煤炭资源总量的 5 3 %。目前煤矿开采深 度以每年 8 ~1 2 m 的速度增加,可以预计在未来 2 0 a我 国很 多煤矿 将进入到 1 0 0 0 1 5 0 0 m 的深 度。如开滦矿区唐山矿 目前的开采深度为一 9 6 0 m, 赵各庄矿 开采水平 已达 一1 1 0 0 ~ 一1 1 5 0 m。国 内外 和开滦矿 区的经验都表明,随着采深的增加 ,地质 收稿 日期2 0 0 4 0 2一O 2 修 回日期 2 0 0 40 41 8 基金项 目国家 自然科 学基金重大项 目 5 049 0 2 7 2 国家重点基础研 究发展 规 J 9 7 3 项 目 2 0 0 2 C B 4 1 2 7 0 5 高等学校博士学科点专项科研基金 2 0 0 3 0 2 9 0 0 0 1 教育部 重大科技资助项 目基金 1 040 5 2 7 作者简介姜耀东 1 9 5 8 一 ,男, 博士,1 9 8 2 年毕业于中国矿业大学,现任教授、 博士生导师, 主要从事岩石力学与采矿工程方面的教学与研究工作。 E ma i l j i a n g y d c u mt b . e d u .c r l 。 维普资讯 l 8 5 8 岩石力学与工程学报 2 0 0 5年 环境更加复杂,地应力增大、涌水量加大、地温升 高,将导致突发性工程灾害和重大恶性事故增加、 作业环境恶化和生产成本急剧增加等一系列问题, 为深部资源开采提 出了严重的挑战u J 。 深部煤炭资源开采的关键技术难题之一是巷道 围岩 控制 问题 ,主要表现 为 2个方面 冲 击地压和 巷道围岩的大变形和强流变性。据调查,开滦唐山 矿在浅部开采时没有发生过冲击地压事故,但在采 深超过 5 3 0 m后冲击地压灾害严重,且随着开采深 度的增加,冲击地压的频度和强度都在增大。开滦 赵各庄矿十水平巷道 采深一 8 7 6 m 返修率为 4 9 .9 %, 十一水平巷道返修率为 7 0 %,进入十二水平,绝大 多数需翻修。 本文在对开滦矿区2个开采深度均超过 1 0 0 0 m 的赵各庄矿和唐山矿巷道围岩矿压显现规律调研 的 基础上,对煤矿深部开采过程中巷道围岩变形、破 坏特征和冲击地压动力显现特征进行了总结,同时 通过计算分析了地应力状态对巷道围岩破坏的影响 关 系。 2 深部开采中巷道围岩稳定性分析 2 . 1 原岩应力特征 原岩应力是引起采矿工程围岩、支护变形和破 坏、产生矿井动力现象的根本作用力。 赵各庄矿采用应力解除法, 测得十一水平 距地 表深度 9 5 2 m 最大主应力为水平应力,其大小为 8 5 .4 MP a ,方 向为倾角 一 0 . 8 7 。 ,方位角 2 7 3 .6 。 , 基本为水平方向,在平面上为近东西方向。 自重应 力 2 8 . 5 MP a ,水平应力是 自重应力的 2 . 9 7倍。 唐 山矿采用孔壁应变法测的十三水平 一 7 5 0 m 最大主应力‘ , 1 5 6 .7 MP a ,其方向为 N E E 7 6 .9 。 ,白 重应力 1 4 .7 MP a ,构造应力是 自重应力的 3 . 8 8 倍 。 2 个矿 的原 岩应力测量 结果表 明,水平 应力远 大于 自重应力,说明地质构造应力在原岩应力场中 占主导地位。开滦矿区的荆各庄矿、钱家营矿和范 各庄矿 的原岩应 力测量结果也相 似 。这主要 是 由于 开滦矿区处于东西 向构造带与 N E,NN E 向构造带 的交汇地区,全区地貌的形成受不同方向的断裂活 动影响, 从而形成了开滦矿区存在 N E E S WW 向的 近水平 的挤 压应力场 。 煤矿开采中所遇到的岩层大多为沉积岩,其强 度 通常较低 ,以赵 各庄 矿十二 煤层为例 ,其 顶底 板 岩层 的单 轴抗压强度 如表 1所示 。 表 1 十二煤层顶底板岩层的抗压 强度 T a b l e 1 Co mp r e s s i v e s t ren g t h s o fr o o fa ndfloo rr o c k i n c 0 a i s e a m No . 1 2 考虑到岩石强度和岩体强度的差异,十二煤层 顶底板岩层实际岩体的抗压强度要比表 1 中所列岩 石的强度低得多。由于原岩应力值本身就大于岩体 的强度,加之采动应力和应力集中的影响,因此, 在其中掘进的巷道围岩将遭遇数倍于其强度 的应力 状态,结果是巷道围岩不可避免地要遭受破坏。 由此可见,在深部煤矿开采中,由于高地应力 和采动应力的双重影响,同时部分巷道不得不在破 碎岩 煤 体 中掘进 如沿 空掘 巷 ,巷道 围岩 的破坏将 是非 常普遍 尤 其是 矿 山巷 道 不可 能用 高成 本 的强 力支护 。 但围岩的破坏并不意味着巷道 的失稳,在 某些情况下却是利用围岩的局部破坏 应力转移 达 到巷道稳定 的 目的 。只 要支护得 当 ,围岩 己破坏 的 巷道仍能保持稳定性I3 】 ,这也是巷道工程与地面 建筑的根本区别之一。因此深部煤矿开采中巷道围 岩控制主要是如何研究 已遭遇变形破坏围岩体的稳 定性 问题 。 2 . 2 巷道围岩的变形破坏特征 1 赵 各庄矿 3 1 3 9上 山巷道变形观 测 3 1 3 9工作面位于十三水平西翼 1 ~3石门,属 十 二 煤 层 。煤 层 沉 积 较 稳 定 ,煤 层 厚 度9 . 1 2 ~ l 4 . 2 7 m, 倾角 1 7 。 ~3 8 。 。老顶深灰色粉砂岩,平均 厚度 3 .4 1 m。直接顶黑色腐泥质粘土岩,平均厚度 5 . 2 0 m。老底灰色细砂岩,平均厚度 2 . 3 6 m,硅质 胶结,成分以石英、长石为主,局部夹粉砂岩,水 平层 理 ,层面有 植物碎屑化 石 。 赵各庄矿十三水平 一1 1 0 0 m 以上采区巷道基 本为 1 8 2 9 U型钢支护, 巷道返修量大, 支护成本 高 。目前 3 1 3 9上山巷道采用锚网和锚索联合支护, 支 护设计如表 2所示 。 对 3 1 3 9上 山巷道锚杆 网和 锚索联合支 护进行 8 个 月巷道变形观 测见表 3 。结果表 明巷道 垂直变形 维普资讯 第 2 4卷第 1 1期 姜耀东等.开滦矿区深部开采中巷道 围岩稳定性研究 表 2 3 1 3 9 上 山支护参数 Ta b l e 2 31 3 9 r o a d wa y s u pp o r t p a r a m e t e r s 表 3 3 1 3 9上山巷道变形观测 T a b l e 3 Obs e r v a t i o n d a t a o f c o n v e r g e nc e de f o r ma t i o n i n 31 3 9 r o a dwa y 注 顶底板移近 量和底鼓量 的单位均为 IT t I n 。 远大于水平变形 ,约为水平变形量 的 2 . 3 ~3 . 8倍 , 在巷道垂直变形中底鼓量约占 7 5 %~9 0 %。 2 唐 山矿 巷道变形观测 唐 山矿 目前 已开拓至十 四水平 一9 6 0 m ,有 3 个主要可采煤层 ,分别为五煤层和八、九合区煤层 及十二煤层L5 J 。表 4为在上述煤层中锚杆支护巷道 围岩变形数据观测统计,表明巷道围岩具有大变形 强流变性, 典型的巷道围岩变形特征有 3种情况 如 图 1 所示 图 1 a 为五煤层中的巷道变形,主要表 现为水平变形大 两帮相对移近 量大 ; 图 1 b 为八、 九合区煤层中巷道的变形,主要表现为严重底鼓; 图 1 c 为十二煤层中巷道的变形,主要表现为严重 冒顶 。 如前所述,在开滦矿区深部开采条件下,巷道 围岩将不可避免地遭受破坏。研究表明,煤矿深井、 软岩、动压巷道围岩的主要破坏形式和变形机制为 挤压流动变形 如图 2所示 ,其特点是巷道的围岩 为 已经遭受过变形破 坏 的软弱 碎裂岩体 ,在深部高 应力和采动应力作用下,这些软弱碎裂围岩变形破 坏过程 中的体积碎胀 流动导致 巷道发 生大位移 ,其 表 4 各煤层巷道收敛变形观测数据 T a bl e 4 Ob s e rv a tio n d a ta o f ro a dwa y c o n v e r g e n c e d e f o r i l l a t i o n i n di f f e r e nt c o a l s e a ms 圆 3 8 00 。 煤层中巷道 底鼓 变形大 单位IT l I n c 十二煤层中 巷道垂直变形大 图 1 巷道变形破坏特征 Fi g . 1 De f o r ma t i o n a n d r u p t ur e c h a r a c t e r i s ti c s o f r o a d wa y s a b 图 2 煤矿深井巷道围岩的挤压流动变形 F i g. 2 S q u e e z i n g a n d flo wi n g d e f o rm a t i o n o f s u r r o u n d i n g r o c k i n d e e p mi n i n g 特征为大变形强流变,巷道底鼓严重 J 。 2 . 3 冲击地压显现特征 冲击地压、煤与瓦斯突出和矿震统称为矿井动 力现象,是煤矿的重大突发灾害之一。随着矿井采 深加大,开滦矿区冲击地压灾害程度逐渐严重。赵 各庄矿深 部地 压动力现象 统计如表 5 【 7 ] 。 根据现 场调查 ,赵 各庄矿深 部冲 击地压动力 显 现有如下特 征 维普资讯 1 8 6 0 岩石力学与工程学报 2 0 0 5笠 表 5 开滦赵各庄矿冲击地压现象统计表 Ta b l e 5 Dy na m i c p he n o m e na s t a t i s t i c s i n Zha o g e z h ua n g c o a l m i ne 叁 地 点 标 高 ,m 强 度 显 现 情 况 2 3 6 lI 1 一9 0 6 . 6 2 一煤炮 线 , 在下面中 一 . 次巨响煤炮 ’ ; 3 盎 一 - 0 0 2 .3 ,次6 凛 、 底 二 中 巷 及 下 出 口 一 1 . 炮 , 次 较 大 煤 一 ’ ⋯ 位置 炮 一 ’ ⋯⋯ ⋯ 2 东 掘 大 、小煤炮 多次 巷道底鼓 ,瓦斯大 s ‘兰 一 - ⋯ 。 备 道 2 3 3 、 进半 ⋯. j ’ 有 震 感 s 一 z 3 次 斜轴部 1 动 力现象 的发 生具有 区域性 。赵 各庄矿 发 生 出现 煤炮 巨响 等地压动力现象的区域, 主要分 布在十 一水平 、十 二水平井 口西 翼开 阔向斜 的轴 部 区域 。如 2 1 3 7 ,2 2 3 7东工作面地压动力现象处在开 阔向斜轴 部 ,2 3 3 7西工作面 地压动力现 象处在次级 向斜 区域 ; 2 赵各庄矿发生煤炮等地压动力现象的临界 深度为一 8 7 6 m, 在临界深度以上没有发生过冲击地 压 3 冲击地压伴有瓦斯异常涌出现象,说 明煤 体在一定范 围内破 坏,导致 瓦斯 涌 出; 4 冲击地 压发生前 ,巷 道底鼓严重 。 唐 山矿 的调 查 也有 类 似 的结 论 1 在 采 深 一 5 3 0 m 以上 未 发现 冲击 地 压现 象 ; 2 在一5 3 O ~ 一 6 0 0 m时发生在 “ 孤岛” 、 “ 半孤岛”煤柱开采的工 作面,共 发生冲击地压现 象 4 7次; 3 当采深达到 一 6 3 0 m 后, 在 正常开采顺序 的工作面 中也发生冲击 地压现象 ,有 4 8 %的冲 击地压发 生在顺槽超前巷道 内; 4 特厚煤层第一分层工作面开采时,冲击地 压严重 ,在 第二 、三 、四分层开采 时则没有发生过 冲击地压 。 开滦矿区的冲击地压显现特征分析至少说明 2 个重要 问题 。首先煤矿深部开采 中即使是软岩或 “ 三软 ”矿 区也有可 能发生冲击地 压,这种冲击 地 压通常 以煤层 的整体突 出为特征 H_ L i p p m a n n , 1 9 9 0 , 在其发生机理上和通常所指的高应力状态下 硬 岩 中的岩爆有 很大 区别 j 。其次 ,高应 力状态和 煤层与顶底板岩层结构之间刚度的差异是导致冲击 地 压 的根本 原因 。特别 指 出,高应 力状 态在冲击地 压 发生过程 中起 着关键主 导作 用 ,否 则就不能解释 为什么在临界深度以下、构造应力区域 内,在 “ 孤 岛” 、“ 半孤岛”煤柱;在顺槽超前巷道 内,第一分 层工作面 开采时冲击地压严 重这一现象 。 3 高地应 力状态对巷道 变形破坏 的 影响研究 如 前所述 ,煤矿 开采 中所遇 到的岩层大 多为强 度较低 的沉积岩 ,进 入深部 开采 后巷道 围岩将遭 遇 数倍于其强度的应力状态,结果是巷道围岩不可避 免地要 遭受破坏 。为 了研究不 同的高应力状 态对 巷 道变形破坏的影响,本文采用 F L A C 2 D 进行数值模 拟 J ,得 到了一些非常有趣 的结论 。 3 . 1 模型的建立 以开滦矿区深部开采为模拟对象,研究不同的 高应 力状态对不 同形状巷道 的变 形破坏影响 。 1 应力状态 静水应 力状态 ,水平 应力大 于垂 直应力 O “h2 a ,垂直应 力大于 水平 应力 O “ v 2 O “ h 。 2 巷道形状 圆形、半圆拱 、梯形和矩 形 。 为了使问题简单化,在模型中假设顶底板岩性 相 同。 3 . 2 模拟结果分析 1 高应力状态对 围岩破坏 的影响 以圆形巷道为例,图 3为静水应力状态下围岩 中的塑性 区形状 ,破坏 区 塑性 圈 为圆形。 图 3 静水压力时主应力及塑性区分布图 F i g . 3 Di s t r i b u t i o n d i a g r a m o f p r i n c i p a l s t r e s s a n d p l a s t i c r e g i o n a t h y d r o s t a t i c s tre s s f i e l d 图 4 a 为水平 应力大 于垂直应 力 O “ h2 下 围岩 中的塑性 区形状 ,其 中浅的阴影部分呈 垂直蝴 蝶形 ,浅的 阴影部分代 表在 该区域 的应力 曾达到塑 性状态但经过应力调整后又回到了弹性状态;深的 维普资讯 维普资讯 1 8 6 2 岩石力学与工程学报 2 0 0 5 盆 状 关 系不 大 ,对 深 井 巷道 围岩 控制 有 积 极 的指 导 意 义 。 4 结语 1 原岩应 力测量 结 果表 明 ,开滦 矿 区在千米 开采 时水平构造 应力远大 于 白重应 力,原岩应 力值 远大于巷道围岩强度值,巷道将不可避免地要遭受 破坏,表现为大变形强流变和严重底鼓; 2 随着开采深度 的增加,软弱煤岩层发生冲 击地压的危险性增大,通常以煤层的整体突出为特 征,其发生机理和高应力状态下硬岩 煤 中的岩爆 有区别 ; 3 在深 井开采 中,原岩 应力状 态对 巷道 围岩 的变形破坏有着重要影响,而巷道形状对巷道围岩 的破坏形式影 响不大 。 当水平应力大 于垂直应力 时, 巷道的破坏主要发生在顶底板方 向;当垂直应力大 于水平应力时,巷道的破坏主要发生在两帮方向, 原岩应力为静水压力时则巷道破坏表现为四周来 压 。 参考文献 R e f e r e n c e s [ 2 】 [ 3 】 钱鸣高,刘听成.矿 山压力及 其控 制 修订本 [ M】 .北京煤炭工业 出版社, 1 9 9 6 . Q i a n Mi n g g a o , L i u T i n g c h e n g . Gr o u n d P r e s s u r e a n d I t s Co n tr o l Re v i s e d [ M1 . Be i j i n g C h i n a C o a l I n d u s t r y P u b l i s h i n g H o u s e . 1 9 9 6 . i n Ch i n e s e 钱鸣高,缪协兴.深部采动岩体中的关键 层活动对冲击地压和突水 等灾害的影响[ AI .见 香 山科学会议第 1 7 5次学术大会[ CI . 【 S . 1 . 】 [ S . n . 】 ,2 0 0 1 .1 01 3 . Qi a n Mi n g g a o ,Mi a o X i e x i n g . E ff e c t s o f t h e ke y s tra m mov e me n t i n the d e e p mi n i n g i n du c e d r o c k ma s s o n the c o a l b u mp s and wa t e r b u r s t [ A] . I n T h e 1 7 5 t h Xi an g s h an S y mp o s i u m o n Mi nin gS c i e n c e an d T ech n o l o g y [ C] . 【 S . 1 . 】 [ S . n . 】 ,2 0 0 1 . 1 01 3 . i n C h i n e s e 何满潮 ,景海河 ,孙晓 明,等.软岩工程力学[ MI .北京 科学 出 版社 , 2 0 0 2 . He Manc h a o , J i n g Ha i h e , S u n X i a o mi n g , e t a1. Mech a n i c s o f S o f t Rock E n g i n e e ri n g [ M] .B e i j i n gS c i e n c e Pre s s ,2 0 0 2 . i n C h i n e s e [ 4 】 刘高,聂德新,韩文峰.高应力软岩巷道围岩变形破坏研 究[ J 】 岩石力学与工程学报, 2 0 0 0 ,1 9 6 7 2 6 7 3 0 . L i u Ga o , Ni e De x i n Ha n W e n f e n g . De f o r m a t i o n a n d f a i l u r e o f s u r r o u n d i n g r oc k s o f r o a d wa y i n h i g h s tr e s s e d s o ft r o c k [ J ] .Ch i n e s e J o u r n a l o f R oc k Me c h ani c s an d E n g i n eeri n g ,2 0 0 0 , 1 9 6 7 2 6 7 3 0 . i n Ch i n e s e [ 5 】 贾德毅 ,杨中东,冯玉,等.唐 I 【 J 矿地质报 告[ RI .开滦开滦集 团唐 山分公司 ,1 9 9 9 . J i a D e y i ,Y a n g Z h o n g d o n g ,F e n g Y u ,e t a1. Ge o l o g i c r e po r t o f T ang s h an c o a l mi n e [ R 1 . K a i l u a n T ang s h an B r anc h o f K a i l u an Mi n i n g Gr o u p Co, 1 9 9 9 . i n C h i n e s e [ 6 】 姜耀东 ,陆士 良.巷道底鼓机理分析[ J 1 . 煤炭学报 ,1 9 9 4,1 9 4 3 4 33 5 1 . J i angY a o d o n g , L u S h i l i a n g . I n v e s t i g a t i o n o f me c h a n i s m o f f l o o r h e a v e o f r o a d w a y [ J ] . J o u rna l o f C h i n a C o a l S o c i e ty, 1 9 9 4 , 1 9 4 3 4 33 5 1 . i n C h i n e s e [ 7 】 周凤增 ,刘占春,宋恩春,等.开滦赵各庄矿深部地质特 征与灾害 防 止 对 策 [ R】 .开滦 集 团 赵 各 庄 矿 , 1 9 9 6 . Z h o u F e n g z e n g,L i u Zh a n c h u n, S o n g En c h u n, e t a1.I n v e s ti g a t i o n o f t h e c h a r a c t e r s o f g e o l o g y an d c o n t r o l o f h a z ard s i n d e e p mi n i n g [ R] . Ka i l u an Mi n i n g Gr o u p Co . , 1 9 9 6 . i n Ch i n e s e 【 8 】 齐 庆新 ,陈 尚本 ,王怀 新 ,等 .冲击 地 压、 岩爆 、矿 震的 关系 及其数值模拟研究[ J 】 . 岩石力学与工程学报. 2 0 0 3 ,2 2 1 1 1 8 5 2 1 8 5 8 . Qi Qi n g x i n , C h e n S h a n g b e n , Wa n g Hu a i x i n , e t a1. S t u d y o n the r e l a t i on s a mo n g c o a l b u mp, r oc k b u r s t an d mi n i n g tre mo r wi th n u me ri c a l s i mu l a t i o n [ J ] .C h i n e s e J o u rna l o f R oc k M e c h ani c s and E n g i n e e ri n g . 2 0 0 3 ,2 2 1 1 1 8 5 21 8 5 8 . i n Chin e s e [ 9 】Co g a n J , P a n t D, J o s h i G M .F LA C a p p l i c a t i o n s f o r t h e K ali Ga n d a k i ’ A ’ h y d r o e l e c t r i c p r o j e c t , Ne p a l [ A] . I n F L AC and N u me ri c a l M o d e l i n g i n Ge o me c h a n i c s [ C] . R o t t e r d a mA. A. B a l k e ma ,1 9 9 9 . 7178 【 l 0 】 多尔恰 尼诺夫 I4 A.构造 应力与井巷 工程稳 定性 [ MI .赵 悖义译. 北京煤炭工 业出版 社,1 9 8 4 . Do r h n o v H A. T e c t o n i c S t r e s s and S t a b i l i t y o f R o a d w a y s [ M] . T r an s l a t e d b y Z h a o C h u n y i . Be i ,ii n g Ch i n a Co al I n d u s t r y Pu b l i s h i n g Ho u s e , 1 9 8 4 . i n Ch i n e s e 维普资讯