[go: up one dir, main page]

CN109061723A - 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统 - Google Patents

一种隧洞岩爆孕育过程的微震源高精度定位方法及系统 Download PDF

Info

Publication number
CN109061723A
CN109061723A CN201810481265.6A CN201810481265A CN109061723A CN 109061723 A CN109061723 A CN 109061723A CN 201810481265 A CN201810481265 A CN 201810481265A CN 109061723 A CN109061723 A CN 109061723A
Authority
CN
China
Prior art keywords
rock burst
microseismic
microquake sources
event
preparation process
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810481265.6A
Other languages
English (en)
Other versions
CN109061723B (zh
Inventor
丰光亮
冯夏庭
陈炳瑞
肖亚勋
李桐
朱新豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Institute of Rock and Soil Mechanics of CAS
Original Assignee
Wuhan Institute of Rock and Soil Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Wuhan Institute of Rock and Soil Mechanics of CAS
Priority to CN201810481265.6A priority Critical patent/CN109061723B/zh
Publication of CN109061723A publication Critical patent/CN109061723A/zh
Application granted granted Critical
Publication of CN109061723B publication Critical patent/CN109061723B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288—Event detection in seismic signals, e.g. microseismics

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Emergency Management (AREA)
  • Business, Economics & Management (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种隧洞岩爆孕育过程的微震源高精度定位方法及系统，通过可实测空间位置的岩爆事件揭示微震源定位方程中三维速度之间的关系，基于该三维速度关系减少隧洞岩爆孕育过程中微震源定位方程的未知数，从而实现基于三维速度模型的隧洞岩爆孕育过程的微震源高精度定位。该方法避免了简化速度参数带来的定位误差，且基于岩爆事件所获取的三维速度关系能够准确恰当用于岩爆孕育过程的微震源定位，从而有效提高了岩爆孕育过程微震源定位的精度。本发明可用于隧洞工程岩爆监测分析。

Description

一种隧洞岩爆孕育过程的微震源高精度定位方法及系统
技术领域
本发明涉及岩土工程领域，更具体涉及一种隧洞岩爆孕育过程的微震 源高精度定位方法及系统。
背景技术
微震监测技术利用微震传感器感知与微震采集设备采集岩石破裂产生 的弹性波信号，经过分析，获得岩石破裂的空间位置、能量、发生时间等 信息，进而科学的解释可能发生的岩石破坏现象。因其具有监测范围广、 自动化、信息化、智能化等特点，且能够全方位实时采集隧洞岩爆孕育过 程的微破裂信息，已被广泛用于隧洞岩爆的监测分析、预警与控制。
微震源定位是微震监测的基础，是基于微震监测技术揭示潜在岩爆风 险位置以及进行岩爆孕育过程分析的重要依据。但在隧洞工程中，因微震 源定位方程中未知数过多且参数相互关联，易造成微震源定位求解不稳定 甚至难以求解。为提高微震源定位求解的稳定性，常对微震源定位过程中 的波速参数进行简化处理，假设微震源到各个微震传感器波速一致，从而 减少未知数，使得定位求解相对稳定。但简化后的波速与实际岩体波速将存在差异，因此会降低微震源定位精度，进而影响对岩爆孕育过程的准确 分析。
发明内容
本发明的目的就是为了解决上述隧洞岩爆孕育过程的微震源定位不稳 定和精度不高的问题，提供一种隧洞岩爆孕育过程的微震源高精度定位方 法及系统。
为了实现上述目的，本发明所设计的隧洞岩爆孕育过程的微震源高精 度定位方法，其特殊之处在于：
通过获取岩爆事件参数，建立基于岩爆事件的三维速度模型微震源定 位方程：
其中：R(a,b,c)为岩爆事件空间坐标，(xi,yi,zi)为第i个微震传感器的 空间坐标，第i个微震传感器监测到的P波和S波到时，第 i个微震传感器的P波和S波三维速度，tR为岩爆事件发生时间；
对方程(1)进行求解，消除未知数tR，获得基于岩爆事件的三维速度 之间的关系
其中，u和v分别为第u个和v个微震传感器编号；
利用上述基于岩爆事件的三维速度之间的关系(2)对该次岩爆孕育过程 的所有微震事件逐一进行定位。首先对该次岩爆孕育过程中发生的第1个 微震事件进行分析，拾取第i个微震传感器对第1个微震事件监测到的P波 和S波到时该微震事件到第i个微震传感器的P波和S波三维速 度与岩爆事件一样，为记该微震事件的空间位置为(a1,b1,c1)， 发生时间为t1。
将第1个微震事件的上述信息代入基于三维速度模型的微震源定位方 程，可得
通过基于岩爆事件的三维速度之间的关系(2)减少第1个微震事件的微 震源定位方程(3)中的未知数，使方程数不小于未知数，进而可对定位方程(3)进行求解，获取极小值，从而精确定位第1个微震事件的空间位置 (a1,b1,c1)。
进一步，按照上述岩爆孕育过程中第1个微震事件的定位方法，对第 2个直至最后一个微震事件逐一进行定位，从而精确获取岩爆孕育过程的所 有微震源空间位置。
本发明还设计了一种隧洞岩爆孕育过程的微震源高精度定位系统，其 特征在于：包括微震监测系统、通讯模块、处理器模块、显示模块、人机 交互模块；
所述微震监测系统捕捉岩爆孕育过程的微震监测数据；
所述通讯模块将微震监测系统所测微震数据传输给处理器模块，
所述处理器模块利用如上所述的隧洞岩爆孕育过程的微震源高精度定 位方法对微震源定位置进行计算，
所述人机交互模块用于输入已知参数，
所述显示模块用于显示微震源位置数据。
进一步度，该系统的处理器模块、显示模块和人机交互模块为工控机， 可通过有线或者无线的方式进行数据的传输。
本发明的优点在于：
1、利用岩爆事件揭示定位方程中的三维速度关系，解决了三维速度难 以获取的问题。该速度关系有效减少了隧洞岩爆孕育过程中微震源定位方 程的未知数，从而使基于三维速度模型的隧洞岩爆孕育过程微震源的定位 稳定可行。
2、三维速度模型更符合隧洞岩体非均质性结构条件，避免了简化隧洞 岩爆孕育过程的微震源定位波速参数带来的误差，从而大幅度提高岩爆孕 育过程的微震源定位精度。
3、岩爆事件与岩爆孕育过程的微震源到微震传感器的震动路径基本一 致，保证了基于岩爆事件所获取的三维速度关系能够准确恰当用于岩爆孕 育过程的微震源定位，从而进一步有效提高隧洞岩爆孕育过程的微震源定 位精度。
具体实施方式
下面结合具体实施例对本发明作进一步的详细描述：
某隧洞发生一起中等岩爆后，经测量其空间坐标为R(127.61,7.6, -33.04)。共有8个微震传感器监测到该次岩爆发生时的震动信号，8个微震 传感器的空间坐标依次分别为(32,9.9,-37.2)、(34,0.3,-28.3)、(30,-9.7, -37.8)、(33,-0.4,-46.3)、(75,10.2,-37.8)、(72,0.3,-29.7)、(76,-9.6,-37.5)和(73, -0.5,-47.1)，8个微震传感器监测到的P波和S波到时依次分别为 (50.01679440,50.02900851)、(50.01656291,50.02860867)、(50.01823320, 50.03149370)、(50.01683596,50.02908029)、(50.00881481,50.01511110)、 (50.00954002,50.01635432)、(50.01055402,50.01809260)和(50.00951394, 50.01630961)。记岩爆位置到第i个微震传感器的P波和S波三维速度为岩爆事件发生时间为tR。
将以上信息代入基于三维速度模型的微震源定位方程后，可得
对方程(1)进行求解，消除未知数tR后，获得基于岩爆事件的三维速度 之间的关系
进一步，利用上述基于岩爆事件的三维速度之间的关系(2)对该次岩爆 孕育过程的所有微震事件逐一进行定位。首先对该次岩爆孕育过程中发生 的第1个微震事件进行分析，拾取8个微震传感器对第1个微震事件监测 到的P波和S波到时，分别为(0.01702404,0.02940516)、(0.01682144, 0.02905521)、(0.01852366,0.03199542)、(0.01708883,0.02951708)、 (0.00902658,0.01547414)、(0.00980215,0.01680368)、(0.01089018,0.01866889)和(0.00976144,0.01673390)，该微震事件到第i个微震传感器的 P波和S波三维速度与岩爆事件一样，为记该微震事件的位置为 (a1,b1,c1)，发生时间为t1。
将第1个微震事件的上述信息代入基于三维速度模型的微震源定位方 程后，可得
通过基于岩爆事件的三维速度之间的关系(2)减少第1个微震事件的微 震源定位方程(3)的未知数，使方程数不小于未知数，进而可对定位方程(3) 进行求解，从而精确定位出第1个微震事件的位置(128.95,9.10,-33.17)。
进一步，按照上述岩爆孕育过程中第1个微震事件的定位方法，对第 2个直至该次岩爆孕育过程的最后一个微震事件逐一进行定位，从而精确获 取岩爆孕育过程的所有微震源空间位置。
本发明所设计的隧洞岩爆孕育过程的微震源高精度定位系统，包括微 震监测系统、通讯模块、处理器模块、显示模块、人机交互模块；
所述微震监测系统捕捉岩爆孕育过程的微震监测数据，
所述通讯模块将微震监测系统所测微震数据传输给处理器模块，
所述处理器模块利用权利要求1所述的方法对微震源定位置进行计 算，
所述人机交互模块用于输入已知参数，
所述显示模块用于显示微震源位置数据。
该系统的处理器模块、显示模块和人机交互模块为工控机，可通过有 线或者无线的方式进行数据的传输。使用者将岩爆坐标和微震传感器坐标 通过人机交互模块输入本系统，其他数据可通过微震监测系统和通讯模块 获取，系统获取数据后，利用上述隧洞岩爆孕育过程的微震源高精度定位 方法计算得出微震源的空间位置，并将结果呈现在显示模块上。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明 所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或 补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权 利要求书所定义的范围。

Claims (3)

1.一种隧洞岩爆孕育过程的微震源高精度定位方法，其特征在于：
通过获取岩爆事件参数，建立基于三维速度模型的微震源定位方程：
其中：R(a,b,c)为岩爆事件空间坐标，(xi,yi,zi)为第i个微震传感器的空间坐标，第i个微震传感器监测到的P波和S波到时，第i个微震传感器的P波和S波三维速度，tR为岩爆事件发生时间；
对方程(1)进行求解，消除未知数tR，获得基于岩爆事件的三维速度之间的关系
其中，u和v分别为第u个和v个微震传感器编号；
对待定位的微震源进行所述的三维速度模型的微震源定位方程的建立，并利用所述的岩爆事件的三维速度之间的关系对微震源进行定位计算，获取极小值，从而获得微震事件的空间位置。
2.一种隧洞岩爆孕育过程的微震源高精度定位系统，其特征在于：包括微震监测系统、通讯模块、处理器模块、显示模块、人机交互模块；
所述微震监测系统通过通讯模块将所获数据传输给处理器模块，
所述处理器模块利用权利要求1所述的方法对微震源定位置进行计算，
所述人机交互模块用于输入已知参数，
所述显示模块用于显示微震源位置数据。
3.根据权利要求2所述的隧洞岩爆孕育过程的微震源高精度定位系统，其特征在于：该系统的处理器模块、显示模块和人机交互模块为工控机。
CN201810481265.6A 2018-05-18 2018-05-18 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统 Active CN109061723B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810481265.6A CN109061723B (zh) 2018-05-18 2018-05-18 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810481265.6A CN109061723B (zh) 2018-05-18 2018-05-18 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统

Publications (2)

Publication Number Publication Date
CN109061723A true CN109061723A (zh) 2018-12-21
CN109061723B CN109061723B (zh) 2020-07-10

Family

ID=64820147

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810481265.6A Active CN109061723B (zh) 2018-05-18 2018-05-18 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统

Country Status (1)

Country Link
CN (1) CN109061723B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110018165A (zh) * 2019-05-08 2019-07-16 中国科学院武汉岩土力学研究所 隧道间歇型岩爆孕育演化过程的监测方法
CN110456413A (zh) * 2019-07-15 2019-11-15 中国科学院武汉岩土力学研究所 预判岩爆类型的方法
CN111897003A (zh) * 2020-08-26 2020-11-06 中国科学院武汉岩土力学研究所 一种考虑传感器阵列影响的微震源定位方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008124759A1 (en) * 2007-04-09 2008-10-16 Baker Hughes Incorporated Location of a microseismic event using the combination of a vsp array and a surface array
CN101770038A (zh) * 2010-01-22 2010-07-07 中国科学院武汉岩土力学研究所 矿山微震源智能定位方法
CN102129063A (zh) * 2010-12-23 2011-07-20 中南大学 一种微震源或声发射源的定位方法
CN106154334A (zh) * 2015-04-13 2016-11-23 中石化石油工程地球物理有限公司胜利分公司 基于网格搜索的井下微地震事件实时反演定位方法
CN107290722A (zh) * 2017-06-29 2017-10-24 中国石油大学(北京) 微震源的定位方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008124759A1 (en) * 2007-04-09 2008-10-16 Baker Hughes Incorporated Location of a microseismic event using the combination of a vsp array and a surface array
CN101770038A (zh) * 2010-01-22 2010-07-07 中国科学院武汉岩土力学研究所 矿山微震源智能定位方法
CN102129063A (zh) * 2010-12-23 2011-07-20 中南大学 一种微震源或声发射源的定位方法
CN106154334A (zh) * 2015-04-13 2016-11-23 中石化石油工程地球物理有限公司胜利分公司 基于网格搜索的井下微地震事件实时反演定位方法
CN107290722A (zh) * 2017-06-29 2017-10-24 中国石油大学(北京) 微震源的定位方法和装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110018165A (zh) * 2019-05-08 2019-07-16 中国科学院武汉岩土力学研究所 隧道间歇型岩爆孕育演化过程的监测方法
CN110018165B (zh) * 2019-05-08 2021-01-01 中国科学院武汉岩土力学研究所 隧道间歇型岩爆孕育演化过程的监测方法
CN110456413A (zh) * 2019-07-15 2019-11-15 中国科学院武汉岩土力学研究所 预判岩爆类型的方法
CN111897003A (zh) * 2020-08-26 2020-11-06 中国科学院武汉岩土力学研究所 一种考虑传感器阵列影响的微震源定位方法
CN111897003B (zh) * 2020-08-26 2021-09-14 中国科学院武汉岩土力学研究所 一种考虑传感器阵列影响的微震源定位方法

Also Published As

Publication number Publication date
CN109061723B (zh) 2020-07-10

Similar Documents

Publication Publication Date Title
CN110095738B (zh) 磁梯度仪参数设计方法
US10508970B2 (en) System for precision measurement of structure and method therefor
US9836993B2 (en) Realistic training scenario simulations and simulation techniques
CN109974749B (zh) 一种三轴转台综合指向误差的评定方法
CN108225370B (zh) 一种运动姿态传感器的数据融合与解算方法
CN102129063A (zh) 一种微震源或声发射源的定位方法
CN100554884C (zh) 挠性陀螺仪最优八位置标定方法
CN109061723A (zh) 一种隧洞岩爆孕育过程的微震源高精度定位方法及系统
CN102539107A (zh) 一种实现风洞试验信号精确同步的方法
CN103808349B (zh) 矢量传感器的误差校正方法和装置
CN106959105B (zh) 罗盘校准方法和装置
CN107861164B (zh) 三维探地雷达的数据处理的方法及装置、三维探地雷达
CN109239654A (zh) 基于神经网络的时差定位结果纠偏方法
US20230324574A1 (en) Stereoscopic identification method and apparatus for disturbance stress evolution process of underground cave surrounding rock
CN107861096A (zh) 基于声音信号到达时间差的最小二乘测向方法
CN108768605A (zh) 管道内检测漏磁和惯导子系统数据的在线同步方法
CN110096779B (zh) 一种伺服机构动特性分析方法
CN104964689A (zh) 基于声源定位、ins的爬壁机器人导航定位系统及方法
CN204788411U (zh) 基于声源定位、ins的爬壁机器人导航定位系统
CN102654515B (zh) 三轴加速度传感器的z敏感轴的校准方法
CN106250649A (zh) 针对星载高精度载荷安装结构在轨变形指向精度预示方法
CN108562269A (zh) 一种相对高度测量方法及装置
CN107084712A (zh) 数据处理方法和装置以及罗盘校准方法和装置
CN104635246A (zh) 卫星导航信号的动态范围检测系统及检测方法
Liu et al. A large scale 3D positioning method based on a network of rotating laser automatic theodolites

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant