CN105022091B - 一种无预测速的远场震源快速定位方法 - Google Patents
一种无预测速的远场震源快速定位方法 Download PDFInfo
- Publication number
- CN105022091B CN105022091B CN201510484618.4A CN201510484618A CN105022091B CN 105022091 B CN105022091 B CN 105022091B CN 201510484618 A CN201510484618 A CN 201510484618A CN 105022091 B CN105022091 B CN 105022091B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- msqrt
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012360 testing method Methods 0.000 title claims abstract description 9
- 230000000694 effects Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000035939 shock Effects 0.000 claims description 4
- 230000002068 genetic effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000004807 localization Effects 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 description 6
- 239000003245 coal Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明的目的是提供一种无预测速的远场震源快速定位方法,本发明在经典TDOA定位方法的基础上,提出一种新的数学模型,在波速未知的情况下,使用降维的方式,用双曲线的渐近线构建寻优目标函数,得到最优解可以快速稳定的收敛,从而快速得到震源的坐标位置。
Description
技术领域
本发明涉及工业监测定位领域,尤其是一种适用于均匀介质环境下远场微震震源快速定位方法,适用于煤矿、石油井下等以震动波、声波、电磁波作传输媒介的工业监测定位系统。
背景技术
微震监测技术是煤矿灾害事故预测的重要手段。研究发现震动活动的频次、能量分布与井下冲击地压、煤与瓦斯突出等矿井动力灾害关系密切。煤矿微震监测是将微震监测技术应用在煤矿井下,记录震动活动频率、能量、波形等信息,利用震波记录初至时间信息定位煤矿井下微震事件的位置、反演煤岩体的纵波速度、分析应力的分布特征等,为预测冲击危险区域提供支持。
微地震监测的基本做法是通过在井中或地面布置检波器排列接收生产活动所产生或诱导的微小地震事件,并通过对这些事件的反演求取微地震震源位置等参数,然后,应用这些参数对生产活动进行监控或指导。长期以来,对微震震源定位方法、提高对微震源定位的准确性和精度的研究,一直是微震监测技术研究的重要内容。震源定位方法很多,主要包括几何方法、物理方法与数学方法等。目前国内外使用的震动源定位方法,大多是以预先测定平均速度或给出平均速度模型为前提,故在现场进行定位时,平均速度测量的准确性直接影响着定位精度。
目前最为常用无需预先测速的震动源定位的以二维函数做寻优目标的到达时间差定位方法,但是这个方法的求解速度和效率都不够,研究发现震源位于远场时,二维目标函数解值在收敛域附近会发生持续震荡,严重影响了求解的稳定性和速度。
因此,有必要研究出一种无预测速的远场震源快速定位方法,从而解决现有技术当中的上述问题。
发明内容
本发明的目的是提供一种无预测速的远场震源快速定位方法,从而能够快速确定震源的位置。
本发明在经典TDOA定位方法的基础上,提出一种新的数学模型,在波速未知的情况下,使用降维的方式,用双曲线的渐近线构建寻优目标函数,得到最优解可以快速稳定的收敛,从而快速得到震源的坐标位置。
本发明请求保护一种无预测速的远场震源快速定位方法,该方法包括如下步骤:
S101,获取N个检波器的位置坐标(x1,y1)、(x2,y2)、(x3,y3)、...(xn,yn),该N个检波器(N≥4)成等距线型排列在待检测定位的环境中,如图2所示;
S102,在待测区域发生微震活动时,各检波器将接收到震动波信号,连接各检波器的采集处理模块将震动信号集中处理,得到各相邻检波器震波的到达时间差,分别记为t12,t23…,t(N-1)N;
S103,将每3个检波器定义为一组检波器,如(x1,y1)、(x2,y2)、(x3,y3)为第一组检波器,对应到时差为t12、t23,得待测位置横坐标与波速一组关系式x=f(v);
认为(x4,y4)、(x5,y5)、(x6,y6)为第二组检波器,对应到时差为t45、t56,根据(1)式得另一组横坐标与波速关系式x′=f(v);
对每两组检波器求分别求|x-x′|的值;
使用遗传算法寻优策略,设定目标函数为所有|x-x′|的和,即为
找到使目标函数为0的v即可认为求解过程完成。
S104,将求得的v和各已知参数代入式(1)和渐近线公式(3)即可求得目标坐标。
附图说明
图1为本发明方法的流程示意图;
图2为实施例中检波器位置的示意图。
具体实施方式
下面通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
以下通过图1对本发明的方法作进一步的介绍。本发明提供一种无预测速的远场震源快速定位方法,该方法包括如下步骤:
S101,获取N个检波器的位置坐标(x1,y1)、(x2,y2)、(x3,y3)、...(xn,yn),该N个检波器(N≥4)成等距线型排列在待检测定位的环境中,如图2所示;
S102,在待测区域发生微震活动时,各检波器将接收到震动波信号,连接各检波器的采集处理模块将震动信号集中处理,得到各相邻检波器震波的到达时间差,分别记为t12,t23…,t(N-1)N;
S103,将每3个检波器定义为一组检波器,如(x1,y1)、(x2,y2)、(x3,y3)为第一组检波器,对应到时差为t12、t23,得待测位置横坐标与波速一组关系式x=f(v);
认为(x4,y4)、(x5,y5)、(x6,y6)为第二组检波器,对应到时差为t45、t56,根据(1)式得另一组横坐标与波速关系式x′=f(v);
对每两组检波器求分别求|x-x′|的值;
使用遗传算法寻优策略,设定目标函数为所有|x-x′|的和,即为
找到使目标函数为0的波速v即可认为求解过程完成。
S104,纵坐标y的计算公式为
将求得的v和各已知参数代入公式(1)和公式(3)即可求得目标坐标(x,y)。
实施例一:
假设检波器坐标为(970,0)、(980,0)、(990,0)、(1000,0),震源坐标为(600,1000),波速为3000m/s,通过数值计算可得,各检波器之间的到时差t21=1.170400ms、t32=1.197628ms、t43=1.224583ms。以本例对二维定位问题进行详细说明,实际定位时,已知量为四个检波器的坐标及传感器触发记录的时刻,微震源的位置未知,在这里之所以给出,目的是通过本发明提出的方法进行验证。具体实施步骤如下:
1、列出四个检波器的坐标值:(970,0)、(980,0)、(990,0)、(1000,0);
2、列出经采集处理后得到的各检波器间的到到时差t21=1.170400ms、t32=1.197628ms、t43=1.224583ms;
3、以标准GA作为寻优方法,以式(2)|x-x′|作为目标函数进行计算。算法运行于MATLAB环境,使用英国谢菲尔德大学开发的GATBX工具箱编程实现,取v∈(1,8165);
通过分析计算的结果,可以发现本发明的方法收敛速度较快且稳定,一般在迭代20次左右即收敛,并且计算最优的v非常稳定,其值恒定为3002.4,已经较接近真实解。
4、将得到稳定的v值代入式(1)、(3),即可得到待测震源位置坐标,计算得到坐标为(600.1,998.8),近似于真实值。
应当理解的是,本发明的应用不限于上述的举例,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (1)
1.一种无预测速的远场震源快速定位方法,其特征在于,该方法包括如下步骤:
S101,获取N个检波器的位置坐标(x1,y1)、(x2,y2)、(x3,y3)、...(xn,yn),该N个检波器成等距线型排列在待检测定位的环境中,其中N≥4;
S102,在待测区域发生微震活动时,各检波器将接收到震动波信号,连接各检波器的采集处理模块将震动信号集中处理,得到各相邻检波器震波的到达时间差,分别记为t12,t23…,t(N-1)N;
S103,将每3个检波器定义为一组检波器,如(x1,y1)、(x2,y2)、(x3,y3)为第一组检波器,对应到时差为t12、t23,得待测位置横坐标与波速一组关系式x=f(v);
<mrow>
<mi>x</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>23</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>12</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>23</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>12</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
定义(x4,y4)、(x5,y5)、(x6,y6)为第二组检波器,对应到时差为t45、t56,根据公式(1)得另一组横坐标与波速关系式x′=f(v);
<mrow>
<msup>
<mi>x</mi>
<mo>,</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>56</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>45</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>56</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mn>45</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
对每两组检波器求分别求|x-x′|的值;
使用遗传算法寻优策略,设定目标函数为所有|x-x′|的和,即为
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mfrac>
<mi>N</mi>
<mn>6</mn>
</mfrac>
</mrow>
</munderover>
<msub>
<mrow>
<mo>|</mo>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mo>|</mo>
</mrow>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mfrac>
<mi>N</mi>
<mn>6</mn>
</mfrac>
</mrow>
</munderover>
<mo>|</mo>
<mfrac>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mrow>
<mn>5</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>4</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>6</mn>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
找到使目标函数为0的波速v即可认为求解过程完成;
S104,纵坐标y的计算公式为
将求得的v和各已知参数代入公式(1)和公式(3)即可求得震源的坐标(x,y)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510484618.4A CN105022091B (zh) | 2015-08-07 | 2015-08-07 | 一种无预测速的远场震源快速定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510484618.4A CN105022091B (zh) | 2015-08-07 | 2015-08-07 | 一种无预测速的远场震源快速定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105022091A CN105022091A (zh) | 2015-11-04 |
CN105022091B true CN105022091B (zh) | 2017-08-25 |
Family
ID=54412190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510484618.4A Expired - Fee Related CN105022091B (zh) | 2015-08-07 | 2015-08-07 | 一种无预测速的远场震源快速定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105022091B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105759311B (zh) * | 2016-01-25 | 2017-11-14 | 西南交通大学 | 一种近实时地震震源位置定位方法 |
CN106990435B (zh) * | 2017-06-07 | 2018-11-23 | 中煤科工集团西安研究院有限公司 | 一种减弱依赖初至拾取精度的微震定位方法及装置 |
CN108957404A (zh) * | 2018-05-18 | 2018-12-07 | 中国科学技术大学 | 一种地震定位方法、装置及系统 |
CN110687593B (zh) * | 2019-10-12 | 2021-08-17 | 中国矿业大学 | 二维小波域矿震监测数据反演方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129063B (zh) * | 2010-12-23 | 2012-10-10 | 中南大学 | 一种微震源或声发射源的定位方法 |
CN102262220B (zh) * | 2011-04-28 | 2013-07-17 | 中南大学 | 一种基于非线性拟合的微震源或声发射源的定位方法 |
CN102435980B (zh) * | 2011-09-15 | 2013-05-08 | 中南大学 | 一种基于解析求解的声发射源或微震源定位方法 |
CN103412283A (zh) * | 2013-08-27 | 2013-11-27 | 中国地质调查局水文地质环境地质调查中心 | 地面塌陷微地震三维定位监测方法和系统 |
CN103605151B (zh) * | 2013-11-20 | 2016-03-02 | 中北大学 | 基于相位测量的分布式群波浅层微震定位方法 |
-
2015
- 2015-08-07 CN CN201510484618.4A patent/CN105022091B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105022091A (zh) | 2015-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105022031B (zh) | 一种区域岩体微震震源的分层速度定位方法 | |
Dong et al. | Velocity-free MS/AE source location method for three-dimensional hole-containing structures | |
CN102129063B (zh) | 一种微震源或声发射源的定位方法 | |
CN104656141B (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN105589100B (zh) | 一种微地震震源位置和速度模型同时反演方法 | |
CN105842735B (zh) | 具有复杂速度分布的区域岩体微震震源定位方法 | |
CN105022091B (zh) | 一种无预测速的远场震源快速定位方法 | |
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
CN109738940B (zh) | 一种存在空区条件下的声发射/微震事件定位方法 | |
CN109991658B (zh) | 一种基于“震源-台站”速度模型的微地震事件定位方法 | |
CN104459763A (zh) | 紧支集小波探测地下空洞位置的方法及系统 | |
CN103323530A (zh) | 利用爆破作业震动波检测巷道群围岩稳定性的方法 | |
CN103852785B (zh) | 地层各向异性的评价方法 | |
CN106249297A (zh) | 基于信号估计的水力压裂微地震震源定位方法及系统 | |
Huang et al. | Source identification of microseismic events in underground mines with interferometric imaging and cross wavelet transform | |
CN103076397A (zh) | 土木工程墙板结构损伤检测的微动测试方法与系统 | |
WO2016032353A1 (ru) | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам | |
CN104749630A (zh) | 构建微地震监测速度模型的方法 | |
Ortiz et al. | Using infrasound waves to monitor tropospheric weather and crater morphology changes at Volcán Tungurahua, Ecuador | |
CN102298158B (zh) | 断层探查装置以及断层探查方法 | |
CN106324670A (zh) | 一种微地震监测系统中的震源定位的方法 | |
CN103364834A (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
CN104570078A (zh) | 一种基于频率域倾角的相似性横向变化率的溶洞检测方法 | |
CN104536046A (zh) | 基于地震记录的震源激发信号一致性评价方法 | |
CN113217109B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170825 Termination date: 20210807 |
|
CF01 | Termination of patent right due to non-payment of annual fee |