CN112526602A - 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 - Google Patents
一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 Download PDFInfo
- Publication number
- CN112526602A CN112526602A CN202011281107.XA CN202011281107A CN112526602A CN 112526602 A CN112526602 A CN 112526602A CN 202011281107 A CN202011281107 A CN 202011281107A CN 112526602 A CN112526602 A CN 112526602A
- Authority
- CN
- China
- Prior art keywords
- wave
- picking
- signal
- arrival time
- long
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000000694 effects Effects 0.000 title claims abstract description 23
- 238000005070 sampling Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000012935 Averaging Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000007774 longterm Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
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/282—Application of seismic models, synthetic seismograms
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
发明提供一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法。该方法包括提取微震信号波形数据、初步判断P波到时所处区间和拾取P波到时等步骤。本方法借助STA/LTA法初步判断P波到时所处区间以及信号与噪声水平差异导致的自回归模型方差激增效应,极大地降低了拾取错误率,提高了拾取效率,同时也增强了P波拾取的稳定性。此方法具有拾取精度高、速度快、稳定性好、适用性强等特点。
Description
技术领域
本发明涉及信号处理技术领域,特别涉及一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法。
背景技术
随着我国经济建设的快速发展,能源开采及存储工程、隧道工程、核废料深埋工程不断兴建,矿震、岩爆等地压灾害剧增,严重威胁了国民生命财产安全,因此有效预测岩爆等灾害已成为刻不容缓的研究课题。微震监测作为一种有效监测技术,已广泛应用于矿震、岩爆等灾害的预警中。微震监测主要包括波形识别、震相拾取、震源定位、震源机制分析和微震活动预测等,其中震相拾取是震源定位以及震源机制分析中最为关键的一步。因S波常叠加于P波尾部,难以识别,故震相拾取一般选用P波到时拾取,其准确性直接影响到震源定位和震源机制分析的准确性,其处理速度直接影响微震监测的效率。
P波到时拾取最初采用人工拾取,虽保证了拾取的精度,但效率太低。因此,国内外学者相继提出了一系列P波到时自动拾取方法,提高了拾取的效率,但它们都存在各自的局限:STA/LTA算法中长短时间窗口的长度、阈值的选择都直接影响了识别效果,且低信噪比或初动不明显时,识别效果较差;小波变换易陷入Fourier分析谱泄露的局限;高阶统计法对波形清晰程度要求较高;PAI-S/K方法引入微震波形峰度和偏斜度函数,计算简单、速度快,但抗噪不佳;分形维数方法和神经网络法计算耗时长。基于自回归模型(AR模型)的AR-AIC算法方法拾取精度较高,但采用双序列AR模型导致拾取效率不高,同时对于低信噪比、刺突信号拾取效果较不稳定。
因此,急需研究出一种拾取精度高、速度快、稳定性好的自动拾取方法。
发明内容
本发明的目的是提供一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,以解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,包括以下步骤:
1)从微震信号中截取待拾取信号波形数据x(n)。其中,n=1,2,…,N。N为微震信号的采样点个数。N=1000~2000。
2)采用STA/LTA法定位包含P波到时的区间X′。
3)拾取X′的VIC自相关最小值点作为P波到时。
进一步,步骤2)具体包括以下步骤:
2.1)给定一个滑动的长时间窗,在长时间窗内取一个短时间窗。所述长时间窗和短时间窗的终点重合。短时窗的长度为ST。长时窗的长度为LT。用短时窗信号平均值STA和长时窗信号平均值LTA之比反映信号振幅的变化。STA/LTA算法公式为:
S(k)=STA(k)/LTA(k) (1)
式中,i为采样时刻,X(i)为在i时刻的关于微震信号的特征函数值。k为采样的第k个点。X(i)=|xi|。
2.2)拾取STA/LTA法首个阈值触发点k1。STA/LTA阀值为3。
2.3)取K=k1+LT-1作为初步P波作为初步P波初至点。
2.4)利用初步P波初至点左右各推移250采样点,得到P波到时所处区间X′=[X(K-250),X(K+250)]。
进一步,步骤3)中由公式(2)计算x(n)的自相关值VIC(k)。选取使VIC最小的k作为P波到时。
进一步,步骤3)之后,还具有对拾取的P波到时进行评价和后续分析处理的相关步骤。
本发明还提供一种P波到时拾取设备,包括:
消息交换模块,用于接收并截取微震信号。
存储器,用于存储计算机程序。
处理器,用于执行所述存储器中存储的计算机程序,实现如权利要求1~4任一所述的P波到时拾取方法。
数据共享模块,用于对P波到时进行报送。
本发明的技术效果是毋庸置疑的:
A.解决了AR-AIC法拾取精度较高,但计算效率低下且对低信噪比、刺突信号拾取不稳定的技术问题;
B.借助STA/LTA法初步判断P波到时所处区间以及信号与噪声水平差异导致的自回归模型方差激增效应,极大地降低了拾取错误率,提高了拾取效率,同时也增强了P波拾取的稳定性;
C.拾取精度高、速度快、稳定性好且适用性强。
附图说明
图1为方法流程图;
图3为AR-VIC法P波到时拾取示意图;
图4为常规信号AR-AIC拾取图;
图5为刺突信号AR-AIC拾取图;
图6为低信噪比信号AR-AIC拾取图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1,本实施例公开一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,包括以下步骤:
1)从微震信号中截取待拾取信号波形数据x(n)。其中,n=1,2,…,N。N为微震信号的采样点个数。N=1000~2000。
2)采用STA/LTA(short term averaging/long term averaging,长短时窗能量比)法定位包含P波(primary wave)到时的区间X′。
2.1)给定一个滑动的长时间窗,在长时间窗内取一个短时间窗。所述长时间窗和短时间窗的终点重合。短时窗的长度为ST,通常取8-15,长时窗的长度为LT,通常取(5-10)ST。用短时窗信号平均值STA和长时窗信号平均值LTA之比反映信号振幅的变化。STA/LTA算法公式为:
S(k)=STA(k)/LTA(k) (1)
式中,i为采样时刻,X(i)为在i时刻的关于微震信号的特征函数值,k为采样的第k个点,其中X(i)=|xi|。
2.2)拾取STA/LTA法首个阈值触发点k1。STA/LTA阀值为3。
2.3)取K=k1+LT-1作为初步P波初至点。本实施例中ST取8,LT取50。
2.4)利用初步P波初至点左右各推移250采样点,得到P波到时所处区间X′=[X(K-250),X(K+250)]。
3)拾取X′的VIC自相关最小值点作为P波到时。由公式(2)计算x(n)的自相关值VIC(k)。选取使VIC最小的k作为P波到时。
参见图2,波形数据其AR模型方差的对数lnσ随着N的增加呈现出稳定-激增-稳定的变化趋势,那么考虑方差激增效应,并引入一个减小因子nlnN/N,从而构造VIC值,使其第一阶段的稳定态变为下降状态,前期下降快而后期下降趋缓,同时保留第三阶段的稳定态,则第一阶段的VIC的极小值点有理由认为是P波到时点。为定阶AR模型利用最小二乘法估计得到的标准差估计值,a为方差激增系数,目的在于放大方差激增幅度,保证第三阶段VIC最终值始终小于第一阶段的VIC极小值。噪声与信号幅值数量级相差越小,其值越大,但不能过大,否则波动将变得很大。a可以考虑取8/log(b),b为信号与噪声最大幅值之比,此时大于20,只要噪声长度x1不过小(x1<100),信号长度x2不过长(x2>2000),nlnN/N肯定小于20。n为减小因子的放大系数,n取太小会造成后期减小因子减幅过小,使得整个VIC值受方差控制,方差波动较大时,一旦方差的最小值出现在VIC拐点前,此点VIC值可能最小,从而造成识别错误。因此噪声波动较大,n应取大值。N为AR模型的容量,经大量试算,n取整个序列容量时AR-VIC法识别效果最佳。
值得说明的是,AR-AIC法采用[X′(1),X′(k)]和[X′(k+1),X′(n)]两个序列的AR模型的AIC值之和的最小值所对应的k作为P波到时,而本发明的STA/LTA_VIC法仅需[X′(1),X′(k)]序列便可完成P波到时拾取,故本发明方法效率更佳。
图4~6中,(a)微震信号原始波形和STA/LTA_VIC法拾取示意图,虚线为人工和STA/LTA_VIC法拾取P波到时(两者相同);(b)整个微震信号对应的AIC值,括号中数字分别对应可能拾取点和该点AIC值;(c)整个微震信号对应的STA/LTA值,括号中数字分别对应可能拾取点和该点STA/LTA值;(d)P波到时所在区间的信号的AIC值,括号中数字分别对应可能拾取点和该点AIC值。
图4至图6表明:相比于AR-AIC法、STA/LTA法和STA/LTA-AIC,STA/LTA-VIC法拾取常规信号、刺突信号和低信噪比信号的P波到时精度更高。将其结果汇于下表1。实际微震信号拾取效果(表1)表明STA/LTA-VIC拾取效果在这几种方法中最优。表1为三种方法拾取结果。从表1可看出,STA/LTA_VIC对常规信号拾取精度高,且有效解决了低信噪比和刺突对AIC拾取结果稳定的影响。
表1
表1中,微震序号1~3分别对应常规信号、刺突信号和低信噪比信号,Manual表示手动拾取,SL表示STA/LTA法,SL-A表示STA/LTA_AIC法,SL-V表示STA/LTA_VIC法。
表2是STA/LTA法、STA/LTA_AIC法和STA/LTA_VIC法拾取误差表,图中50组微震信号为加入一定信噪比的高斯噪声的声发射信号。由表2知STA/LTA_VIC法拾取精度最高,其次是STA/LTA_AIC法,STA/LTA法最差,且在低信噪比下,STA/LTA_VIC法仍有极佳的拾取精度,此外随着信噪比的降低,三种方法拾取结果的差异逐渐增大。
表2
4)对拾取的P波到时进行评价和后续分析处理。微震时差定位法利用拾取的P波到时进行微震定位。通过识别P波到时从而提取出P波初动振幅,继而进行微震类型的判别。
实施例2:
本实施例公开一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,包括以下步骤:
1)从微震信号中截取待拾取信号波形数据x(n)。其中,n=1,2,…,N。N为微震信号的采样点个数。N=1000~2000。
2)采用STA/LTA(short term averaging/long term averaging,长短时窗能量比)法定位包含P波(primary wave)到时的区间X′。
3)拾取X′的VIC自相关最小值点作为P波到时。
实施例3:
本实施例公开一种P波到时拾取设备包括:
消息交换模块,用于接收并截取微震信号。
存储器,用于存储计算机程序。
处理器,用于执行所述存储器中存储的计算机程序,实现如实施例1或2任一所述的P波到时拾取方法。
数据共享模块,用于对P波到时进行报送。
Claims (5)
1.一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,其特征在于,包括以下步骤:
1)从微震信号中截取待拾取信号波形数据x(n);其中,n=1,2,…,N;N为微震信号的采样点个数;N=1000~2000;
2)采用STA/LTA法定位包含P波到时的区间X′;
3)拾取X′的VIC自相关最小值点作为P波到时。
2.根据权利要求1所述的一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,其特征在于,步骤2)具体包括以下步骤:
2.1)给定一个滑动的长时间窗,在长时间窗内取一个短时间窗;所述长时间窗和短时间窗的终点重合;短时窗的长度为ST;长时窗的长度为LT;用短时窗信号平均值STA和长时窗信号平均值LTA之比反映信号振幅的变化;STA/LTA算法公式为:
S(k)=STA(k)/LTA(k) (1)
式中,i为采样时刻,X(i)为在i时刻的关于微震信号的特征函数值;k为采样的第k个点;X(i)=|xi|;
2.2)拾取STA/LTA法首个阈值触发点k1;STA/LTA阀值为3;
2.3)取K=k1+LT-1作为初步P波作为初步P波初至点;
2.4)利用初步P波初至点左右各推移250采样点,得到P波到时所处区间X′=[X(K-250),X(K+250)]。
4.根据权利要求1所述的一种基于长短时窗和AR模型方差激增效应的P波到时拾取方法,其特征在于:步骤3)之后,还具有对拾取的P波到时进行评价和后续分析处理的相关步骤。
5.一种P波到时拾取设备,其特征在于,包括:
消息交换模块,用于接收并截取微震信号;
存储器,用于存储计算机程序;
处理器,用于执行所述存储器中存储的计算机程序,实现如权利要求1~4任一所述的P波到时拾取方法;
数据共享模块,用于对P波到时进行报送。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011281107.XA CN112526602B (zh) | 2020-11-16 | 2020-11-16 | 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011281107.XA CN112526602B (zh) | 2020-11-16 | 2020-11-16 | 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526602A true CN112526602A (zh) | 2021-03-19 |
CN112526602B CN112526602B (zh) | 2023-10-20 |
Family
ID=74980991
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011281107.XA Active CN112526602B (zh) | 2020-11-16 | 2020-11-16 | 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112526602B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113189644A (zh) * | 2021-04-30 | 2021-07-30 | 哈尔滨工业大学(威海) | 一种微震震源定位方法及系统 |
CN114002733A (zh) * | 2021-10-27 | 2022-02-01 | 武汉科技大学 | 微震波信号初至到时自动拾取方法及微震监测装置 |
CN115146678A (zh) * | 2022-07-04 | 2022-10-04 | 长江水利委员会长江科学院 | 一种爆破振动信号的p波振相初至识别方法和系统 |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5995448A (en) * | 1996-11-20 | 1999-11-30 | Krehbiel; Steven | Method for mapping seismic reflective data |
RU2480740C1 (ru) * | 2011-11-25 | 2013-04-27 | Общество с ограниченной ответственностью "ИНКОТЕС" | Способ измерения задержек распространения ультразвуковых волн в анизотропных средах и устройство для его осуществления |
JP2014169960A (ja) * | 2013-03-05 | 2014-09-18 | Taisei Corp | 地震の主要動到達時刻の予測方法 |
US20150164957A1 (en) * | 2012-05-18 | 2015-06-18 | The University Of Chicago | Conversion of cardiomyocytes into fast conducting cardiomyocytes or slow conducting nodal cells |
CN104914468A (zh) * | 2015-06-09 | 2015-09-16 | 中南大学 | 一种矿山微震信号p波初至时刻联合拾取方法 |
JP2015215221A (ja) * | 2014-05-09 | 2015-12-03 | 大成建設株式会社 | 地震の主要動の到達判定方法および判定システム |
CN105223614A (zh) * | 2015-09-23 | 2016-01-06 | 中南大学 | 一种基于dwt_sta/lta的含噪信号p波初至峰度拾取方法 |
CA2910878A1 (en) * | 2014-07-28 | 2016-01-28 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus and method for selecting one of a first encoding algorithm and a second encoding algorithm using harmonics reduction |
CN105487114A (zh) * | 2015-12-08 | 2016-04-13 | 中南大学 | 一种微震信号p波初至点综合拾取方法 |
CN105527650A (zh) * | 2016-02-17 | 2016-04-27 | 中国科学院武汉岩土力学研究所 | 一种工程尺度下微震信号及p波初至自动识别算法 |
CN106896407A (zh) * | 2017-03-28 | 2017-06-27 | 吉林大学 | 一种基于近似负熵的微地震信号初至拾取方法 |
CN108919353A (zh) * | 2018-07-03 | 2018-11-30 | 华北科技学院 | 一种微震波形初至到时的自动分级拾取与优选方法 |
CN109283576A (zh) * | 2018-11-26 | 2019-01-29 | 辽宁工程技术大学 | 一种以振幅为特征函数的自动拾取p波震相方法 |
CN110018515A (zh) * | 2019-03-13 | 2019-07-16 | 辽宁工程技术大学 | 一种用振幅变化长短时均值比实现地震预警中p波自动拾取的方法 |
KR102003382B1 (ko) * | 2018-03-23 | 2019-07-24 | 동일테크주식회사 | 실시간 지진 판별 및 조기 경보 기능을 구비한 지진동 기록 장치 및 이에 사용되는 지진동 판별 방법 |
CN110646844A (zh) * | 2019-09-30 | 2020-01-03 | 东北大学 | 基于波形包络线的隧道岩石破裂微震s波到时拾取方法 |
CN110907991A (zh) * | 2019-12-11 | 2020-03-24 | 重庆大学 | 基于数据场势值的震源定位方法、系统及可读存储介质 |
CN111175810A (zh) * | 2019-07-05 | 2020-05-19 | 中南大学 | 微震信号到时拾取方法、装置、设备及存储介质 |
CN111736208A (zh) * | 2020-06-24 | 2020-10-02 | 重庆大学 | 变权重联合P波和S波初至数据的微震事件Bayes定位方法、系统及介质 |
-
2020
- 2020-11-16 CN CN202011281107.XA patent/CN112526602B/zh active Active
Patent Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5995448A (en) * | 1996-11-20 | 1999-11-30 | Krehbiel; Steven | Method for mapping seismic reflective data |
RU2480740C1 (ru) * | 2011-11-25 | 2013-04-27 | Общество с ограниченной ответственностью "ИНКОТЕС" | Способ измерения задержек распространения ультразвуковых волн в анизотропных средах и устройство для его осуществления |
US20150164957A1 (en) * | 2012-05-18 | 2015-06-18 | The University Of Chicago | Conversion of cardiomyocytes into fast conducting cardiomyocytes or slow conducting nodal cells |
JP2014169960A (ja) * | 2013-03-05 | 2014-09-18 | Taisei Corp | 地震の主要動到達時刻の予測方法 |
JP2015215221A (ja) * | 2014-05-09 | 2015-12-03 | 大成建設株式会社 | 地震の主要動の到達判定方法および判定システム |
CA2910878A1 (en) * | 2014-07-28 | 2016-01-28 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus and method for selecting one of a first encoding algorithm and a second encoding algorithm using harmonics reduction |
CN104914468A (zh) * | 2015-06-09 | 2015-09-16 | 中南大学 | 一种矿山微震信号p波初至时刻联合拾取方法 |
CN105223614A (zh) * | 2015-09-23 | 2016-01-06 | 中南大学 | 一种基于dwt_sta/lta的含噪信号p波初至峰度拾取方法 |
CN105487114A (zh) * | 2015-12-08 | 2016-04-13 | 中南大学 | 一种微震信号p波初至点综合拾取方法 |
CN105527650A (zh) * | 2016-02-17 | 2016-04-27 | 中国科学院武汉岩土力学研究所 | 一种工程尺度下微震信号及p波初至自动识别算法 |
CN106896407A (zh) * | 2017-03-28 | 2017-06-27 | 吉林大学 | 一种基于近似负熵的微地震信号初至拾取方法 |
KR102003382B1 (ko) * | 2018-03-23 | 2019-07-24 | 동일테크주식회사 | 실시간 지진 판별 및 조기 경보 기능을 구비한 지진동 기록 장치 및 이에 사용되는 지진동 판별 방법 |
CN108919353A (zh) * | 2018-07-03 | 2018-11-30 | 华北科技学院 | 一种微震波形初至到时的自动分级拾取与优选方法 |
CN109283576A (zh) * | 2018-11-26 | 2019-01-29 | 辽宁工程技术大学 | 一种以振幅为特征函数的自动拾取p波震相方法 |
CN110018515A (zh) * | 2019-03-13 | 2019-07-16 | 辽宁工程技术大学 | 一种用振幅变化长短时均值比实现地震预警中p波自动拾取的方法 |
CN111175810A (zh) * | 2019-07-05 | 2020-05-19 | 中南大学 | 微震信号到时拾取方法、装置、设备及存储介质 |
CN110646844A (zh) * | 2019-09-30 | 2020-01-03 | 东北大学 | 基于波形包络线的隧道岩石破裂微震s波到时拾取方法 |
CN110907991A (zh) * | 2019-12-11 | 2020-03-24 | 重庆大学 | 基于数据场势值的震源定位方法、系统及可读存储介质 |
CN111736208A (zh) * | 2020-06-24 | 2020-10-02 | 重庆大学 | 变权重联合P波和S波初至数据的微震事件Bayes定位方法、系统及介质 |
Non-Patent Citations (5)
Title |
---|
WANG GL 等: "Decay constants of P-wave mesons", 《PHYSICS LETTERS B》, vol. 650, no. 1, pages 15 - 21, XP022109062, DOI: 10.1016/j.physletb.2007.05.001 * |
王桂林 等: "基于改进灰狼算法的岩石声发射定位研究", 《工程地质学报》, vol. 30, no. 1, pages 234 - 241 * |
王进铭: "矿山微震信号P波初至时间拾取与震源定位方法研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》, no. 2, pages 021 - 494 * |
董耀;高鹏举;金路;李恒;肖娟;: "微动探测在城市地质勘查中的应用研究", 能源与环保, no. 12, pages 92 - 96 * |
袁明德: "工程地震勘探技术的进展", 地球物理学进展, no. 04, pages 127 - 132 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113189644A (zh) * | 2021-04-30 | 2021-07-30 | 哈尔滨工业大学(威海) | 一种微震震源定位方法及系统 |
CN113189644B (zh) * | 2021-04-30 | 2023-03-07 | 哈尔滨工业大学(威海) | 一种微震震源定位方法及系统 |
CN114002733A (zh) * | 2021-10-27 | 2022-02-01 | 武汉科技大学 | 微震波信号初至到时自动拾取方法及微震监测装置 |
CN114002733B (zh) * | 2021-10-27 | 2024-01-23 | 武汉科技大学 | 微震波信号初至到时自动拾取方法及微震监测装置 |
CN115146678A (zh) * | 2022-07-04 | 2022-10-04 | 长江水利委员会长江科学院 | 一种爆破振动信号的p波振相初至识别方法和系统 |
CN115146678B (zh) * | 2022-07-04 | 2023-05-09 | 长江水利委员会长江科学院 | 一种爆破振动信号的p波振相初至识别方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112526602B (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112526602A (zh) | 一种基于长短时窗和ar模型方差激增效应的p波到时拾取方法 | |
CN104914468B (zh) | 一种矿山微震信号p波初至时刻联合拾取方法 | |
CN107479094B (zh) | 一种实现地震预警的方法 | |
CN110907991B (zh) | 基于数据场势值的震源定位方法、系统及可读存储介质 | |
CN105487114B (zh) | 一种微震信号p波初至点综合拾取方法 | |
CN103646648B (zh) | 一种噪声功率估计方法 | |
CN109581313B (zh) | 一种脉组频率捷变雷达的目标检测方法 | |
CN105223482A (zh) | 局部放电特高频信号波形的小波分解二值去噪方法 | |
CN112364296B (zh) | 一种基于深度学习的p波到时自动拾取方法 | |
CN112332807B (zh) | 一种弱包络信号检测方法和系统 | |
US20170296081A1 (en) | Frame based spike detection module | |
CN114002733A (zh) | 微震波信号初至到时自动拾取方法及微震监测装置 | |
CN110929842A (zh) | 非合作无线电信号突发时间区域精确智能检测方法 | |
CN112649512B (zh) | 一种岩体声发射初至自适应识别方法 | |
CN110275190A (zh) | 一种基于查表法的北斗信号高精度捕获方法 | |
CN111538086B (zh) | 一种提高地震数据初至波质量的初至自动拾取方法 | |
CN115508805B (zh) | 一种机载激光雷达海洋探测回波信号分类方法 | |
CN116776085A (zh) | 基于模糊聚类和赤池信息准则的微震信号到时拾取方法 | |
CN109270536B (zh) | 一种多普勒天气雷达的信号处理器 | |
CN115510898A (zh) | 一种基于卷积神经网络的舰船声尾流检测方法 | |
CN114114400A (zh) | 微地震事件有效信号拾取方法 | |
CN110007342B (zh) | 一种用于低信噪比地震信号的时频域直接拾取初至方法及系统 | |
CN109298447B (zh) | 一种低信噪比声发射信号的到时拾取方法 | |
CN113219526A (zh) | 一种地震波初至自动拾取与筛选方法 | |
CN114325722B (zh) | 基于水下声信标信号多脉冲累积的高增益检测方法和系统 |
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 |