基于微震监测技术的深井开采地压活动规律研究_图文

1、第26卷 第4期岩石力学与工程学报 V ol.26 No.42007年4月 Chinese Journal of Rock Mechanics and Engineering April ，2007基于微震监测技术的深井开采地压活动规律研究杨承祥12，罗周全1，唐礼忠1，(1. 中南大学 资源与安全工程学院，湖南 长沙 410083；2. 安徽铜都铜业股份有限公司 冬瓜山铜矿，安徽 铜陵 244000摘要：简要介绍微震监测技术基本原理及冬瓜山微震监测系统的组成和结构，运用微震监测系统对冬瓜山矿床首采区段不同时段的地压活动集中区进行圈定，采用量化地震学原理开展岩层应力变形强度分布研究，揭示开采诱

2、发地压活动的时空变化规律，初步评价盘区隔离矿柱和回采采场的地压活动状况及稳定性。研究表明，地压活动相对集中区的变化与井下采掘活动紧密相关，各地压活动相对集中区与采掘工程施工位置相对应，在时空上随采掘活动的改变而发生变化。采场回采区内的地压活动相对较稳定，掘进活动区内的地压活动分布较分散且随采掘活动的结束很快减小或消失，目前首采区段岩层是稳定的。实践证明，微震监测技术为深井矿床开采地压活动规律研究提供一种有效技术手段，研究成果为指导矿山安全生产奠定重要的基础。 关键词：采矿工程；微震监测技术；深井开采；金属矿山；地压活动中图分类号：TD 32 文献标识码：A 文章编号：10006915(2007

3、04081807STUDY ON RULE OF GEOSTATIC ACTIVITY BASED ON MICROSEISMICMONITORING TECHNIQUE IN DEEP MININGYANG Chengxiang12，LUO Zhouquan1，TANG Lizhong1，(1. School of Resources and Safety Engineering，Central South University，Changsha ，Hunan 410083，China ；2. Dongguashan Copper Mine，Anhui Tongdu Copper Co.，L

4、td. ，Tongling ，Anhui 244000，China Abstract ：By introducing the basic principle of the microseismic monitoring technique and the components and structure of the ISS microseismic monitoring system used in Dongguashan Mine，the relative concentrating area of the geostatic activities is compartmentalized

5、 in the first mining area of Dongguashan ore deposit within different periods. The system and the principle of the seismology are adopted to study the intensity distribution of the stress and displacement to discover the temporal and spatial change rules of the geostatic activity produced by the min

6、ing and to evaluate the geostatic activity situation and the stability of the insulating pillar and main mining area. The study indicates that the change of the concentrating area of geostatic activity is dependent on the underground mining activities. Every concentrating area is corresponding to th

7、e position of each mining project. The geostatic activity changes in time and space with the change of mining project position accordingly. The geostatic activity within the mining area of stope is relatively stable. The distribution of geostatic activity in excavating area is separate and reduces q

8、uickly or even disappears with the completion of excavation. The rock mass within the first mining area of Dongguashan ore deposit is stable now. The practice proves that the microseismic monitoring technique is an effective method to study the rule of geostatic activity in deep mining. The result s

9、ettles an important base for Dongguashan Mine to produce safely.Key words：mining engineering；microseismic monitoring technique；deep mining；metal mine；geostatic activity收稿日期：20060519；修回日期：20061115基金项目：国家“十五”科技攻关计划项目(2004BA615A04 ；国家自然科学基金重大项目(50490274作者简介：杨承祥(1965 ，男，博士，1987年毕业于中南工业大学采矿工程专业，现为高级工程师，主

10、要从事矿床深井开采与灾害控制方面的研究工作。E-mail ：tlycq第26卷 第4期 杨承祥，等. 基于微震监测技术的深井开采地压活动规律研究 819 1 引 言近十几年来，随着数字技术、计算机技术、地球物理学和量化地震学的快速发展，微震监测技术在国外特别是南非、加拿大等国的深井矿山中得到了广泛应用1。我国也有少数矿山开展矿山地压活动的微震监测，如门头沟煤矿于20世纪80年代就采用波兰SYLOK 微震监测系统进行地压活动的监测，兴隆庄煤矿于21世纪初采用澳大利亚研制的地震监测系统、凡口铅锌矿于2004年引进加拿大ESG 微震监测系统26进行地压活动的监测。当前，我国矿山进入深井开采已呈必然趋

11、势，为实现安全生产，开展基于微震监测技术的深井矿山岩体破坏和地压活动规律研究显得迫切而重要。冬瓜山矿床赋存于6901 007 m，主井深度为1 120 m(+951 025 m，具有原岩应力高(38 MPa 、岩温高(39 、存在岩爆倾向等复杂开采技术条件，是国内目前开采深度达到1 000 m的大型深埋铜矿山之一，属典型的深井开采矿山7。为保证矿山高效安全开采，“十五”期间，“深井岩爆与地压监测及控制技术研究”课题组引进南非ISS 微震监测系统，建立了以微震监测为主的岩爆与地压综合监测系统，开展对岩体应力活动的动态实时监测，采集岩体活动的振动波形数据，进行地震学处理、分析和可视化研究，初步揭示

12、了矿床因开采引起的岩体应力活动和变形规律。2 微震监测基本原理大量研究资料表明，岩体在破坏之前，必然持续一段时间以声波的形式释放积蓄的能量。矿床开采活动在岩体中引起弹性变形和非弹性变形，在岩体中积蓄的弹性势能在非弹性变形过程中以震动波的形式被逐步或突然释放出去，这种能量释放的强度，随着结构临近失稳而变化，导致岩体内部发生微震(microseismos，MS 事件。这种通过分析微震事件产生的信号(位置、震级等参数 特征，推断开挖过程中的岩体状态和矿岩的力学行为，估测矿岩是否发生破坏，以实现防止、控制或预测潜在的不稳定岩体，从而避免危险事故发生的技术，称为微震监测技术8。微震信号的特征取决于震源性

13、质、所经岩体性质及监测点到震源的距离等。描述微震事件采用定量地震学方法。在震源周围布置一定数量的传感器，组成传感器网络，当岩体内出现微震时，传感器即可将由岩体内部破坏处所产生的应力波信号拾取，通过多点同步数据采集测定各传感器接收到该信号的时刻，连同各传感器坐标及所测波速代入非线性方程：LOC j (x =t j t 0T j (h (1式中：t j 为在测点j 所观测到事件的P 波或S 波的到时时间，t 0为该事件的未知震源时间，T j (h 为P 波或S 波传播到第j 测点的未知传播时间，LOC j (x 为观测到时时间t j 与计算到时时间t 0+T j (h 之差。通过求解5个以上方程构

14、成方程组的唯一解，可确定震源的时空参数，达到定位之目的。矿床开采活动产生应力集中。要监测岩体对开采的反应，就必须连续不断地在时间和空间上量化描述应力和应变状态流动变化的参数。地震事件分析是对出现于一定体积V 和一定时间t 内的许多地震事件数量和强度的分析。常采用4个独立的地震活动参数(E ，M ， 的对数形式对地震活动进行量化和定义。其中，为事件间的平均时间，为相临事件之间的平均距离(包括震源大小 ，M 为地震矩的和，E 为发射能量的和。由这4个基本量可以推得十几个地震参数，如地震应变s 、地震应力s 、视体积V A 、地震Schmidt 数Sc s 等。这些参数描述了同震变形及其相关的应变率

15、、应力和流变特性的统计特性。与本文相关的主要参数和公式为lg E (ergs = 2.4m + 5.8 (2m =23lg M 6. 1 (3A =E /(A (4V A =M /(c 3A =M 2/(c 3GE (5式中：m 为震级；E 为发射的地震能(J；A 为视应力(Pa； 为震源处的平均位移(m；A 为震源面积(m2 ；c 3为尺寸因素，且c 32；G 为刚度；M 为地震矩(N·m 。3 微震监测系统组成和结构冬瓜山微震监测系统由硬件和软件组成。硬件系统包括地震传感器(sensor、数据采集单元地震仪(quake seismometer，QS 、井下通信控制中心(seism

16、ic·820· 岩石力学与工程学报 2007年controller 、地面监测控制中心、终端用户计算机及通讯电源电缆。软件系统包括控制和管理微震监测系统的运行控制软件(RTS；对采集的波形进行地震波波形分析、处理和参数计算，提供地震学分析平台的地震学分析软件(JMTS；以及在三维窗口中显示对采集地震数据分析的各类图像，提供多种参数的时间序列曲线和图表，满足不同空间和时间范围地震活动研究需要的微震事件可视化解释分析软件采用尾砂胶结充填，矿柱采场用尾砂充填。首采区 段划分4个盘区，盘区和采场分布见图2。2004年开始回采，至2006年9月，522#，526#，528#采场已回采

17、结束，但尚未充填。虽然采动范围大，但采空区体积较小且相对孤立。回采区域主要集中于1#2#盘区南端，掘进工作主要是矿体上部的凿岩硐室施工和矿体底部的巷道掘进和拉底工程施工。 4.2 地压活动时空变化及相对集中区圈定由于矿山在不同空间区域和时间段进行不同工序的开采活动，为揭示地压活动的时空变化规律，有必要进行地压活动的时空特征分析1015(JDI6。微震监测系统组成与结构如图1所示。该系统由16个传感器、4个QS 、1个地震仪转发器(QS_Rep、1个地下控制器、1个地表监测控制中心及与之相连的通讯电缆组成9。每个QS 连接4个传感器，其中3个为一维传感器，1个为三维传感器。各传感器采集地震模拟信

18、号通过QS 转换为数字信号传输到井下通信控制中心，再通过光缆传输到地表监测控制中心进行处理和分析。该系统的有效监测范围为600 m×400 m×220 m，矿床首采区域震源定位误差10 m，系统灵敏度为2.0里氏震级。结合矿山实际开采情况，按月划分时间段进行地压活动分析。这里仅以2005年9月和2006年1月的地压活动分布及其集中区圈定来说明研究过程和结果。图3(a，(b分别是2005年9月微震事件空间三维分布透视图及其水平投影图。采用系统中JDI 软件在图3中对首采区域地震事件进行地压活动集中区圈定，图3(c给出首采区地震事件及其空间范围，图3(d为地压活动相对集中区在水

19、平面上的投影范围(线圈所围区域 。从图中可清晰看出地压活动相对比较集中的区域，集中区较小且比较分散。图中球体表示地震事件，球的位置为地震震源位置，球的大小为地震矩对数，球的颜色由深到浅分别表示地震事件发生的先后时间。一个地震事件也可由像素点、圆圈、砂漏等符号表示，其大小可由不同的地震参数如震级、地震矩、震源半径、视应力、释放的能量等多种参数表示，同样，球体的4 地压活动规律研究4.1 首采区段采掘活动冬瓜山矿床首采区段用阶段空场嗣后充填采矿法开采。矿体缓倾斜，沿矿体走向划分盘区，盘区之间暂留隔离矿柱，其宽度为18 m。沿盘区长度方向划分采场，矿房采场长82 m，矿柱采场长78 m，采场宽均为1

20、8 m。先采矿房，后采矿柱，矿房采场监控中心计算机客户端计算机 (信息中心局域网微震监测系统监控软件SuSE 计算机操作系统地震学分析软件 可视化解释分析软件地表 地下光缆QS1670 m水平730 m水平QS2875 m水平QS4井下通信控制中心单对铜绞线875 m水平单对铜绞线730 m水平单对铜绞线QS3875 m水平QS Repeater单对铜绞线单对铜绞线图1 微震监测系统组成与结构Fig.1 Components and structure of microseismic monitoring system第26卷 第4期 杨承祥，等. 基于微震监测技术的深井开采地压活动规律研究

21、821 注：图2中4个角上的数字是矿区地理坐标；单位：m图2 首采区段盘区和采场分布Fig.2 Layout of panels and stopes in the first mining area (a 微震事件空间分布透视图 (b 地压集中区水平投影图(c 微震事件空间分布范围 (d 地压集中区水平投影范围52线图4 不同时段地压活动相对集中区变化 Fig.4 Change of relatively concentrating area of geostaticactivity in different periods活动相对集中区发生了明显的变化。对比井下采掘活动记录可知，地压活动相

22、对集中区的变化与井下采掘活动紧密相关，各地压活动相对集中区与采掘工程位置相对应。首采区各个地压活动相对集中区之间相对独立，表示地压活动较弱，采掘活动引起岩层活动之间的相互影响小。 4.3 岩层应力及变形强度分布采用地震学参数平面震源上的平均位移( 和视应力(A ，在水平面上分别绘制和lg A 的等值线图，研究其分布特征，分析首采区段岩层应力和变形及其强度分布。这里仍以2005年9月和2006年1月2个时段的地压活动研究为例加以讨论。图3 2005年9月地压活动相对集中区Fig.3 Relatively concentrating area of geostatic activitiesin S

23、ept. 2005颜色亦可用于表达其他地震参数。采用同样的方法可获得2006年1月的地压活动相对集中区的空间分布及其在水平面上的投影范围。把2005年9月和2006年1月的地压活动相对集中区重叠后如图4所示，其中2005年9月的地压活动相对集中区范围线用实线表示。显然，地压·822· 岩石力学与工程学报 2007年图5是2005年9月首采区段的和lg A 等值线图。从图5(a可以看出，井下采掘引起的位移分布靠近52与53线的回采区、盘区北部和东部区域。最大位移集中区域位于528#采场的53线顶板岩体，位移约为0.04 m，其他位移区的位移都小于0.01 m，发生在首采区域采

24、掘范围附近。这些变形与这期间采场顶板凿岩和采场底部结构施工以及回采爆破有关。图5(b表明，较大的视应力主要分布于采场开采活动区，其量值大约为5.8 Pa。值得注意的是，在56线8#10#采场隔离矿柱附近引起了较大的视应力集中区，量值约为8.3 Pa。经现场调查发现，该处掘进岩石为有岩爆倾向的石英闪长岩，地质钻机在此附近进行探矿活动。结合位移分布可以看出，应力释放除在很小的范围内引起较大位移外，在绝大部分范围内并没有产生明显的位移，说明其量值并不高。56线附近局部出现的视应力集中区与该处的位移发生区基本吻合，说明有一定的位移，但位移并不大。 (a (单位：m (b lgA (单位：Pa图5 20

25、05年9月岩体应力及平均位移等值线图 Fig.5 Isograms of and lgA in the first mining area in Sept.2005图6是2006年1月首采区段的和lg A 等值线图。对比图6(a与图5(a可以看出，该时间段内首采区段位移有两个特点：一是发生位移的区域缩小了，说明位移区域发生了明显变化；二是回采范围发生较大位移区内位移增大到0.046 m。比较图6(b与图5(b可以看出，仍然有两个相对较大的应力集中区，其中量值较大的集中区位于回采范围内，但该集中区的范围相对较小，视应力有所下降，量值分别约为5.2和5.9 Pa。回采区域应力主要集中在采场位置，最

26、大应力集中区在有岩爆倾向的石英闪长岩中。上述分析表明，通过地震位移和视应力参数可了解首采区段岩层应力和变形大小的变化及其分布，利用其等值线图可以明确划分岩层变形和应力 (a (单位：m (b lgA (单位：Pa图6 2006年1月岩体应力及平均位移等值线图 Fig.6 Isograms of and lgA in the first mining area in Jan. 2006分布集中区域。首采区段目前应力和位移活动区主要分布于回采和掘进作业区域，相对而言，位移较大的区域主要集中于回采区，掘进区位移作业一旦结束位移就停止。此外，应力和变形是波动的，并不是随着开采的进行一定会增大，也存在减

27、小的变化。 4.4 盘区隔离矿柱应力及变形研究监测盘区矿柱的应力活动和稳定状态对矿区安全生产具有重要意义。目前，采场回采活动主要分布于54线隔离矿柱两侧，因此有必要对54线隔离矿柱围岩中的应力和变形及其分布进行研究。绘制不同时段内隔离矿柱岩体的位移、视应力对数lg A 、累积视体积(cumulative apparentvolume ，CA V 半径等值线图，分析不同时间段内隔离矿柱岩体的应力和变形分布。图7为2006年6月30日前90 d和前300 d隔离矿柱中的位移与lg A 等值线图、位移与CA V 半径等值线图。从前300 d的累积图形来看，位移总量较小，最大值为0.015 m。相对较

28、大的累积位移主要分布在526#527#采场和545#546#采场区域，且靠近采场的顶板岩层，这是由该区域采场回采比较集中所致；另外，在采场的西北端开挖及其变形引起的位移也较大。最大视应力位于隔离矿柱的WN 端，量值为5.6 Pa。结合位移分布特征和累积视体积分布特征，说明该区域岩体位移由该区域井巷围岩应力集中形成。从长期来看，累积视体积约等于被开挖出的岩体体积，显然，最大累积视体积集中于回采区。从图7可以看出，不同时间段内位移、视应力和累积视体积的空间和强度分布发生了明显变化，到后期(前90 d，隔离矿柱WN 端应力、视位移及累积视体积都减小或消失了，进一步证明隔离矿柱WN 端的地压活动及岩层

29、应力和变形是由井巷掘进所引起。第 26 卷 第 4 期 杨承祥，等. 基于微震监测技术的深井开采地压活动规律研究 823 (a 前 300 d 位移和视应力等值线图 (b 前 300 d 位移和 CAV 半径等值线图 (c 前 90 d 位移和视应力对数等值线图 (d 前 90 d 位移和 CAV 半径等值线图 Fig.7 图 7 隔离矿柱应力及变形分布(位移单位：m；应力单位：Pa Distributions of stress and displacement of insulating pillar(unit of displacement：m；unit of stress：Pa 4.5

30、 采场围岩地压活动研究 用分析隔离矿柱的方法对 522#试验采场进行 地压活动研究。图 8 是系统运行(2005 年 9 月 1 日 起30，90 和 300 d 后，522#采场纵剖面上的围岩 应力和位移分布图。从图中可以看出，系统运行后 30 d 内产生位移的范围较小，且集中于采场的上下 端；随着开采的进行，产生位移的范围从该位置沿 底部结构水平逐渐扩大， 但是累积位移量变化不大， 量值小于 0.003 m。围岩应力变化相对比较明显，与 采场开采过程密切相关，随着开挖体积不断增大， 视应力集中区范围扩大，其量值增大，但最大视应 力值并不随其增大，相反还减小，30 d 为 8.5 Pa， 9

31、0 d 为 8.0 Pa，300 d 为 6.0 Pa。从位移和视应力变 化范围可看出，开挖过程中岩体的状态，产生位移 和应力的范围在扩大，对周围岩体影响在增强。位 移和视应力量值较小，说明采场的地压活动较弱， 采场的岩层是稳定的。 (b 系统运行后 90 d (c 系统运行后 300 d 图8 (a 系统运行后 30 d 522#采场纵剖面上的围岩应力和位移分布(位移单 位：m；应力单位：Pa Fig.8 Distributions of stress and displacement of rockmass on longitudinal section of stope #522(uni

32、t of displacement：m；unit of stress：Pa ·824· 岩石力学与工程学报 2007 年 Chinese Journal of Rock Mechanics and Engineering ， 2002 ， 5 结论与建议 6 21(Supp.2：2 6092 612.(in Chinese 赵向东，王育平，陈 波，等. 微地震研究及在深部采动围岩监测 中的应用J. 合肥工业大学学报(自然科学版，2003，26(3：363 367.(ZHAO Xiangdong，WANG Yuping，CHEN Bo，et al. On the micros

33、eismic technology and its application to monitoring of rock stability in deep mine caving processesJ. Journal of Hefei University of Technology(Natural Science， 2003， 26(3： 363367.(in Chinese 7 杨承祥，罗周全. 深井高应力高温缓倾斜矿床开采技术研究J. 金 属矿山，2005，(4：710.(YANG Chengxiang，LUO Zhouquan. Research on mining technolo

34、gy for high stress high stress and high temperature gently in inclined deep depositJ. Metal Mine， 2005， (4： 710.(in Chinese 8 MENDECKI A J. Seismic monitoring in minesM. London，UK： Chapman and Hall Press，1997. 9 唐礼忠，杨承祥，潘长良. 大规模深井开采微震监测系统站网布置 优化J. 岩石力学与工程学报，2006，25(10：2 0362 042.(TANG Lizhong，YANG C

35、hengxiang，PAN Changliang. Optimization of microseismic monitoring network for large-scale deep well miningJ. Chinese Journal of Rock Mechanics and Engineering，2006，25(10： 2 0362 042.(in Chinese 10 EBRAHIM-TROLLOPE R，JOOSTE Y. Seismic hazard quantifycationC/ POTVIN Y，HUDYMA M ed. Controlling Seismic

36、Risk， Proc. of the 6th International Symposium on Rockburst and Seismicity in Mines. Nedlands：Australian Centre for Geomechanics Press，2005： 157164. (1 冬瓜山数字化微震监测系统应用研究初步 揭示了现有开采条件下首采区段地压活动的规律。 研究表明，首采区段应力和位移活动区主要分布在 回采和掘进作业区域。 (2 目前的地压活动主要由采掘活动引起，各 地压活动集中区与采掘工程相对应，在时空上随采 掘活动的改变而发生变化；累积位移量和应力都较 小，说明

37、地压活动较弱，总体来看，岩层是稳定的。 (3 基于微震监测技术矿床开采诱发的地压活 动规律研究有待于进一步深入。需要不断开展微震 事件的震源空间分布、震级、事件发生频率等与地 质构造及开挖过程之间的关系研究，地压活动和深 井围岩应力分布状态的特征及变化形式研究，预测 岩爆和地压危险区和危险程度等研究。实践证明， 微震监测技术为深井矿床开采地压活动规律研究提 供了一种有效技术手段。可以预见，随着研究工作 的深入，该方法可为采取合理有效的措施控制矿山 工程岩体破坏与岩爆地压灾害提供依据，研究成果 可为指导矿山安全高效生产奠定基础。 参考文献(References： 1 ORTLEPP W D. R

38、aSiM comes of agea review of the contribution to the understanding and control of mine rockburstsC/ POTVIN Y， HUDYMA M ed. Controlling Seismic Risk ， Proc. of the 6th International Symposium on Rockburst and Seismicity in Mine. Nedlands：Australian Centre for Geomechanics Press，2005：320. 2 鲁振华，张连成. 门

39、头沟矿微震的近场监测效能评估J. 地震， 1989，(5：3239.(LU Zhenhua，ZHANG Liancheng. Evaluation of efficiency of near field microseismic monitoring in Mentougou MineJ. Earthquake，1989，(5：3239.(in Chinese 3 姜福兴， XUN Luo. 微震监测技术在矿井岩层破裂监测中的应用J. 岩土工程学报，2002，24(2：147149.(JIANG Fuxing，XUN Luo， Application of microseismic monit

40、oring technology of strata fracturing to underground coal mineJ. Chinese Journal of Geotechnical Engineering，2002，24(2：147149.(in Chinese 4 李庶林，尹贤刚，郑文达，等. 凡口铅锌矿多通道微震监测系统及其 应用研究J. 岩石力学与工程学报，2005，24(12：2 0482 053.(LI Shulin ， YIN Xiangang ， ZHENG Wenda ， et al. Research on multi-channel microseismic monitoring system and its application to Fankou Lead-zinc MineJ. Chinese Journal of Rock Mechanics and Engineering，2005，24(12：2 0482 053.(in Chinese 5 赵向东，陈 波，姜福兴. 微地震工程应用研究J. 岩石力学与工 程学报， 2002， 21(增 2： 6092 612.(ZHAO Xian