CN103984017B - 一种微震震源自动定位方法 - Google Patents
一种微震震源自动定位方法 Download PDFInfo
- Publication number
- CN103984017B CN103984017B CN201410236187.5A CN201410236187A CN103984017B CN 103984017 B CN103984017 B CN 103984017B CN 201410236187 A CN201410236187 A CN 201410236187A CN 103984017 B CN103984017 B CN 103984017B
- Authority
- CN
- China
- Prior art keywords
- deviation
- effective
- positioning method
- linear function
- passage
- 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.)
- Active
Links
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种微震震源自动定位方法。包括:根据获取的条有效通道的信息构建出个四元一次方程组;为有效通道数,;对所述个四元一次方程组中每一个四元一次方程组求解,得到组解;把所述组解分别带入震源定位目标方程中,得到个对应每条有效通道的偏差值;判断所述个偏差值是否均满足误差范围,如果是,输出所述震源定位目标方程中的偏差目标函数值,否则,改变个偏差值中不满足误差范围的偏差值对应的有效通道信息,然后执行根据所述个有效通道的信息构建个四元一次方程组的步骤。本发明的定位方法通过计算机来实现,省时省力。
Description
技术领域
本发明涉及微震检测领域,尤其涉及一种微震震源自动定位方法。
背景技术
在煤矿开采活动中有裂缝等岩层破坏事件时,岩体的瞬间破裂会激发弹性波。这些弹性波携带着震源的信息,依赖岩体弹性介质向四周传播。矿山微震监测就是利用震动传感器在远处测量这些弹性波信号,并根据所监测的微震信号特征来确定破裂的发生时间、空间位置、尺度、强度及性质,做出有效合理的预警,具有远程、三维和实时的特点,因此在矿山、岩土、水利水电等工程中具有广阔的应用前景,日益受到人们的重视。而矿山微震检测技术的核心是真元定位,能否进行快速、准确的确定为是微震检测系统能否发挥作用的关键。
传统的震源定位计算会将震源定位的目标方程转换成对四元一次方程求解的问题,理论上当给出已知的检波器坐标和标尺位置,就能求出方程组的解,实际上,虽然构造的多个方程组能求出多个解,但是未必有满足目标函数的解,这就要求检波器的坐标和我们对标尺位置的标注尽量准确,这样,当标出位置后,发现某一条通道目标函数较大时,就要重新对该条通道进行标注。对于多条通道,每一条通道都要标出准确位置,既费时又费力。
发明内容
为解决上述技术问题,本发明的目的是提供一种微震震源自动定位方法,通过计算机来实现定位计算,达到省时省力的目的。
本发明的一种微震震源自动定位方法,包括:
根据获取的条有效通道的信息构建出个四元一次方程组;为有效通道数,;
对所述个四元一次方程组中每一个四元一次方程组求解,得到组解;
把所述组解分别带入震源定位目标方程中,得到个对应每条有效通道的偏差值;
判断所述个偏差值是否均满足误差范围,如果是,输出所述震源定位目标方程中的偏差目标函数值,否则,改变所述个偏差值中不满足误差范围的偏差值对应的有效通道信息,然后执行根据所述个有效通道的信息构建个四元一次方程组的步骤。
其中,有效通道信息包括:标尺相对位置、标尺相对时间和检波器坐标。
进一步地,得到个对应每条有效通道的偏差值之后,输出所述震源定位目标方程中的偏差目标函数值之前,还包括:
根据每一条有效通道偏差值计算个对应每条有效通道的均方值;
判断所述个均方值中最小均方值是否小于预设阈值,如果是,则执行判断所述个偏差值是否均满足误差范围的步骤,否则,所述震源定位目标方程无解。
具体地,根据获取的条有效通道的信息构建出个四元一次方程组,包括:
根据获取的条有效通道的信息构建出个定位计算公式;
分别对所述个定位计算公式两边平方,得到个平方后的定位计算公式;
将所述个平方后的定位计算公式两两相减,得到个四元一次方程;
构建出个四元一次方程组。
其中定位计算公式为,为震源坐标,为第个测点的坐标,为震源震动时间,为P波到达第个测点的时间,为P波的传播速度。
其中,震源定位目标方程为
,其中,为偏差目标函数,为有效通道数,为震源坐标,为震源震动时间,为P波到达第个测点的时间,为第个测点的坐标,为P波的传播速度。
具体地,改变个偏差值中不满足误差范围的偏差值对应的有效通道信息,包括:
判断个偏差值中不满足误差范围的偏差值是否大于0,如果是,则减小不满足误差范围的偏差值对应的有效通道标尺位置,否则,增大不满足误差范围的偏差值对应的有效通道标尺位置。
具体地,对个四元一次方程组中每一个四元一次方程组采用全选主元高斯-约当消去法求解。
借由上述方案,本发明至少具有以下优点:
本发明实施例只须动标注有效通道信息一次,然后就可以通过计算机代替手工方法寻找最优解,在某条有效通道信息偏差值不满足误差范围时,计算机也自动修改有效通道信息,重新寻找最优解,即偏差目标函数值也就是震源坐标实际和理论的差值,相对于现有技术的手工寻找最优解省时省力。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
图1是本发明实施例提供的一种微震震源自动定位方法示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
参见图1,一种微震震源自动定位方法,包括如下步骤:
1、根据获取的条有效通道信息构建出个四元一次方程组;为有效通道数,;
有效通道信息可以是通过自动标注获取的信息,每条有效通道信息包括:标尺相对位置、标尺相对时间和检波器坐标。
标尺是指发生震动事件的那一刻所在的时间段(本发明实施例中定位5秒钟)中,标尺标注的位置就是这5秒中微震动起震的那一刻。假如事件发生在2014.1.2.15:20:30,该时间段即为15:20:28~15:20:32这5秒数据,标尺就应该标在微震动起震位置15:20:30这一刻。
获取上述条有效通道信息后,构建个定位计算公式:
具体地,定位基本公式为:
(1)
其中,为震源坐标,为第个测点的坐标,为震源震动时间,为P波到达第个测点的时间,为P波均匀的传播速度。有条有效通道,则有个定位计算公式(1)。
分别对上述个定位计算公式(1)两边平方,得到个平方后的定位计算公式;
定位基本公式为公式(1)时,公式(1)两边平方后则变为
(2)
有个定位计算公式(1),则对应有个公式(2)。
将上述个公式(2)两两相减,得到个四元一次方程的标准形式,于是,可以构建出(表示平面定位)个四元一次方程组。其中,C是求组合个数。
本发明实施例以为例进行说明,根据4条有效通道信息构建4个定位计算公式(1),则对应得到4个公式(2)。将4个公式(2)两两相减,得到个四元一次方程,因此,构建出个四元一次方程组。
2、对上述构建的个方程组中每一个方程组进行求解,得到组解。
具体地,可以采用全选主元高斯-约当消去法、复系数方程组的全选主元高斯消去法、复系数方程组的全选主元高斯-约当消去法、对角线方程组的追赶法、大型稀疏方程组的全选主元高斯-约去消去法、求解托伯利兹方程组的列文逊方法等。本发明实施例以采用全选主元高斯-约当消去法进行求解。
3、把步骤2中得到的组解分别带入震源定位目标方程中,得到个对应每条有效通道的偏差值;
4、判断上述个偏差值是否均满足误差范围,如果是,则执行步骤5;否则,执行步骤6;
5、输出震源定位目标方程中的偏差目标函数值,偏差目标函数值也就是震源坐标实际和理论的差值。
6、改变个偏差值中不满足误差范围的偏差值对应的有效通道信息,然重新执行步骤1。
其中,震源定位目标方程为
(3)
公式(3)的约束条件为。
公式(3)中为偏差目标函数,为有效通道数,为震源坐标,为震源震动时间,为P波到达第个测点的时间,为第个测点的坐标,为P波的传播速度。
根据公式(3)可得到每一条有效通道的偏差值,即。
进一步地,在步骤3之后,步骤4之前,还包括:
根据每一条有效通道偏差值计算个对应每条有效通道的均方值;
判断上述个均方值中最小均方值是否小于预设阈值,如果是,则执行步骤4,否则,震源定位目标方程无解。
若每一条有效通道的偏差值分别为,则每条有效通道的均方值为。
通过判断上述个均方值中最小均方值是否小于预设阈值,可以确定公式(3)是否存在最优解,如果是,则公式(3)存在最优解;如果否,则公式(3)不存在满足条件的解,直接结束算法,这样可节省计算时间。
进一步地,步骤4中改变个偏差值中不满足误差范围的偏差值对应的有效通道信息包括:
判断个偏差值中不满足误差范围的偏差值是否大于0,如果是,则减小不满足误差范围的偏差值对应的有效通道标尺位置,否则,增大该不满足误差范围的偏差值对应的有效通道标尺位置。
本发明实施例采用该震源自动定位方法进行了10次测试,并计算每一次算法的循环次数和耗时,结果如下表所示:
从上表可以看出,自动求解算法的循环耗时一般都在1S之内,而且保证误差值均在10ms以内,达到求最优解的要求。
本发明实施例只须动标注有效通道信息一次,然后就可以通过计算机代替手工方法寻找最优解,在某条有效通道信息偏差值不满足误差范围时,计算机自动修改不满足误差范围的偏差值对应的有效通道信息,重新寻找最优解,即偏差目标函数值也就是震源坐标实际和理论的差值,相对于现有技术的手工寻找最优解省时省力。
以上所述仅是本发明的优选实施方式,并不用于限制本发明,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。
Claims (7)
1.一种微震震源自动定位方法,其特征在于:
根据获取的n条有效通道的信息构建出k个四元一次方程组;n为有效通道数,n≥4;
对所述k个四元一次方程组中每一个四元一次方程组求解,得到l组解,且l满足条件l≤k;
把所述l组解分别带入震源定位目标方程中,得到n个对应每条有效通道的偏差值,且l满足条件l≤k;
判断所述个偏差值是否均满足误差范围,如果是,输出所述震源定位目标方程中的偏差目标函数F(x0,y0,z0)值,否则,改变所述n个偏差值中不满足误差范围的偏差值对应的有效通道信息,然后执行根据所述n个有效通道的信息构建k个四元一次方程组的步骤;
得到n个对应每条有效通道的偏差值之后,输出所述震源定位目标方程中的偏差目标函数值之前,还包括:
根据每一条有效通道偏差值计算n个对应每条有效通道的均方值;
其中,每一条有效通道的所述偏差值分别为F(1),F(2)…F(n),则每条有效通道的均方值为[F(1)2+F(2)2+...F(n)2]/n;
判断所述n个均方值中最小均方值是否小于预设阈值,如果是,则执行判断所述n个偏差值是否均满足误差范围的步骤,否则,所述震源定位目标方程无解。
2.根据权利要求1 所述的微震震源自动定位方法,其特征在于,所述有效通道信息包括:标尺相对位置、标尺相对时间和检波器坐标。
3.根据权利要求1 所述的微震震源自动定位方法,其特征在于,根据获取的n 条有效通道的信息构建出k 个四元一次方程组,包括:
根据获取的n 条有效通道的信息构建出n 个定位计算公式;
分别对所述n 个定位计算公式两边平方,得到n 个平方后的定位计算公式;
将所述n 个平方后的定位计算公式两两相减,得到m=个四元一次方程;构建出k=个四元一次方程组。
4.根据权利要求3所述的微震震源自动定位方法,其特征在于,所述定位计算公式为=v(ti-t0),其中(x0,y0,z0)为震源坐标,(xi,yi,zi)为第i个测点的坐标,t0为震源震动时间,ti为P 波到达第i个测点的时间,v为P 波的传播速度。
5.根据权利要求1 所述的微震震源自动定位方法,其特征在于,所述震源定位目标方程为F(x0,y0,z0)= P,其中,F(x0,y0,z0)为偏差目标函数,n为有效通道数,(x0,y0,z0)为震源坐标,t0为震源震动时间,ti为P 波到达第i个测点的时间,(xi,yi,zi)为第i个测点的坐标,v为P 波的传播速度,P表示此偏差目标函数为P波的目标函数。
6.根据权利要求1 所述的微震震源自动定位方法,其特征在于,所述改变n 个偏差值中不满足误差范围的偏差值对应的有效通道信息,包括:
判断所述n 个偏差值中不满足误差范围的偏差值是否大于0,如果是,则减小不满足误差范围的偏差值对应的有效通道标尺位置,否则,增大不满足误差范围的偏差值对应的有效通道标尺位置。
7.根据权利要求1 所述的微震震源自动定位方法,其特征在于,对所述k个四元一次方程组中每一个四元一次方程组采用全选主元高斯-约当消去法求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410236187.5A CN103984017B (zh) | 2014-05-30 | 2014-05-30 | 一种微震震源自动定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410236187.5A CN103984017B (zh) | 2014-05-30 | 2014-05-30 | 一种微震震源自动定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103984017A CN103984017A (zh) | 2014-08-13 |
CN103984017B true CN103984017B (zh) | 2017-01-18 |
Family
ID=51276059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410236187.5A Active CN103984017B (zh) | 2014-05-30 | 2014-05-30 | 一种微震震源自动定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103984017B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403914A (zh) * | 2014-08-26 | 2016-03-16 | 长沙矿山研究院有限责任公司 | 一种微震震源联合定位方法 |
CN105785436B (zh) * | 2016-03-17 | 2018-08-14 | 北京矿冶研究总院 | 一种矿用微震监测方法 |
CN107765340B (zh) * | 2017-10-18 | 2019-08-27 | 北京科技大学 | 基于电磁辐射和微震联合监测的井下煤岩破裂源定位方法 |
CN108831081A (zh) * | 2018-07-26 | 2018-11-16 | 广东电网有限责任公司 | 一种基于源定位技术电缆防外力破坏预警方法 |
CN113176606B (zh) * | 2021-06-02 | 2023-09-26 | 中国恩菲工程技术有限公司 | 微震震源定位方法、系统、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129063B (zh) * | 2010-12-23 | 2012-10-10 | 中南大学 | 一种微震源或声发射源的定位方法 |
CN102778668A (zh) * | 2012-07-23 | 2012-11-14 | 中煤科工集团西安研究院 | 矿山被动震源快速精确定位方法 |
-
2014
- 2014-05-30 CN CN201410236187.5A patent/CN103984017B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103984017A (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103984017B (zh) | 一种微震震源自动定位方法 | |
CN105740203B (zh) | 多传感器无源协同测向定位方法 | |
CN103759732B (zh) | 一种角度信息辅助的集中式多传感器多假设跟踪方法 | |
CN105842735B (zh) | 具有复杂速度分布的区域岩体微震震源定位方法 | |
CN103697999B (zh) | 一种高应力硬岩tbm施工隧道微震波速实时获取方法 | |
CN105223576A (zh) | 一种基于单矢量潜标的线谱信号目标自动检测方法 | |
CN105589100A (zh) | 一种微地震震源位置和速度模型同时反演方法 | |
CN102495425A (zh) | 一种基于能量的微地震震源自动定位方法 | |
CN110174643A (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
CN103389039A (zh) | 一种基于fpga和fpaa的高速高精度光幕检测装置 | |
CN104684081A (zh) | 基于距离聚类遴选锚节点的无线传感器网络节点定位算法 | |
CN103410516A (zh) | 小直径隧洞全断面开挖围岩变形预警工法 | |
CN104749626B (zh) | 微地震事件点的定位方法 | |
CN104776827B (zh) | Gps高程异常数据的粗差探测方法 | |
CN103033716B (zh) | 一种电网综合负荷模型中各负荷成分所占比例的计算方法 | |
CN105954741A (zh) | 一种基于多假设拟蒙特卡罗的多目标无源协同定位方法 | |
US20230184983A1 (en) | Vector-resistivity-based real-time advanced detection method for water-bearing hazard body | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN103292773A (zh) | 一种基于最小区域的对称度误差评定方法 | |
CN103742156A (zh) | 深埋硬岩隧道贯通前相向改单向掘进时机与方式确定方法 | |
CN103926569B (zh) | 三坐标雷达网基于交叉定位点关联的集中式与分布式压制干扰鉴别方法 | |
CN110514963A (zh) | 一种改进双端行波故障测距方法 | |
CN104180822B (zh) | 一种变形监测基准点稳定性检验方法 | |
CN105043390A (zh) | 基于泛克里金法的重力场插值方法 | |
CN103177192B (zh) | 五孔探针角度特性曲线数据处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |