CN108535778B - 一种地震波初至检测方法及系统 - Google Patents
一种地震波初至检测方法及系统 Download PDFInfo
- Publication number
- CN108535778B CN108535778B CN201810768697.5A CN201810768697A CN108535778B CN 108535778 B CN108535778 B CN 108535778B CN 201810768697 A CN201810768697 A CN 201810768697A CN 108535778 B CN108535778 B CN 108535778B
- Authority
- CN
- China
- Prior art keywords
- data
- time window
- earthquake record
- record sequence
- seismic
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 39
- 238000000546 chi-square test Methods 0.000 claims abstract description 12
- 239000013589 supplement Substances 0.000 claims description 17
- 230000001502 supplementing effect Effects 0.000 claims description 5
- 235000013399 edible fruits Nutrition 0.000 claims description 4
- 238000000034 method Methods 0.000 abstract description 19
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/41—Arrival times, e.g. of P or S wave or first break
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种地震波初至检测方法及系统。该方法包括:获取检波器采集的待检测的地震记录序列;建立第一时间窗和紧随第一时间窗的第二时间窗,第一时间窗和第二时间窗的容量均为奇数;对地震记录序列的左侧和右侧进行数据填充;从第一时间窗的中心对准地震记录序列的第一个数据滑动第一时间窗和第二时间窗,直至第一时间窗的中心对准地震记录序列的最后一个数据,每滑动一个数据的步长,对第一时间窗内的数据与第二时间窗内的数据进行卡方检验,从而确定地震记录序列中的各个数据处是否有地震波动到达。本发明的方法及系统,能够有效避免噪声对检测精度的影响,提高检测准确度。
Description
技术领域
本发明涉及地震波初至检测技术领域,特别是涉及一种地震波初至检测方法及系统。
背景技术
地震波初至在矿物资源开发、地震定位以及地震层析成像中都起着十分关键的作用。在石油勘探技术中,地震波信号可以用来了解地质构造,例如石油勘探就需要利用人工地震信号,其中准确检测到地震波到达检波点的初至时刻尤为重要。
从根本上说,地震信号的初至波既包含地震信号又包含噪声信号,是两者的叠加,初至拾取的精确性是地震资料处理中基础而又重要的环节。传统的初至拾取主要依靠手工,但人工拾取耗时且精确度低。为了提高地震波初至的提取效率与提取速度,现在已经出现了自动拾取的方法,例如能量比法、最大振幅法、分形维法等。
初至自动拾取方法中最常用的是基于能量比的方法。根据能量比达到最大值时的时间点就可以估算出地震波初至到来的时刻。但值得注意的是,当信噪比相对较低时这种基于能量比的方法就十分不理想。与能量比法原理相类似,振幅特性也常被用来检测初至,利用这种方法,根据地震波初至前后的振幅比的变化来拾取初至。一般来说,当信噪比足够高时这些方法是有效的。但我国石油勘探环境并不都是理想的,地表复杂的地区噪声干扰非常严重,这导致了现有方法对低信噪比情况下的初至拾取并不理想,精度并不高。
发明内容
本发明的目的是提供一种地震波初至检测方法及系统,有效避免噪声对检测精度的影响,提高检测准确度。
为实现上述目的,本发明提供了如下方案:
一种地震波初至检测方法,包括:
获取检波器采集的待检测的地震记录序列;
建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数;所述第二时间窗紧随所述第一时间窗;
对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据;
从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达;
记录所有有地震波动到达的数据在所述地震记录序列中的序号,并计算每个有地震波动到达的数据对应的初至。
可选的,所述对所述地震记录序列的左侧和右侧进行数据填充,具体包括:
将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;
将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
可选的,所述每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达,具体包括:
求取所述第一时间窗内所有数据的最大值和最小值,得到取值区间;
将所述取值区间划分为K段,得到K段小区间;
计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k);其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;k为小区间序号,k=1,2,…,K;
计算数量差值比率:
根据参数K和置信度α,通过查询卡方表获取卡方值χ2(α,K-1);
判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;
若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,则确定在所述地震记录序列的第n个数据处没有地震波动到达。
本发明还公开一种地震波初至检测系统,包括:
序列获取模块,用于获取检波器采集的地震记录序列;
时间窗建立模块,用于建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数;所述第二时间窗紧随所述第一时间窗;
数据填充模块,用于对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据;
地震波初至检测模块,用于从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达;
记录与初至计算模块,用于记录所有有地震波动到达的数据在所述地震记录序列中的序号,并计算每个有地震波动到达的数据对应的初至。
可选的,所述数据填充模块,具体包括:
序列左填充单元,用于将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;
序列右填充单元,用于将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
可选的,所述地震波初至检测模块,具体包括:
区间计算单元,用于求取所述第一时间窗内所有数据的最大值和最小值,得到取值区间;
区间划分单元,用于将所述取值区间划分为K段,得到K段小区间;
区间数据量计算单元,用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k);其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;k为小区间序号,k=1,2,…,K;
比率计算单元,用于计算数量差值比率:
卡方值查询单元,用于根据参数K和置信度α,通过查询卡方表获取卡方值χ2(α,K-1);
判断单元,用于判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;
地震初至确定单元,用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,则确定在所述地震记录序列的第n个数据处没有地震波动到达。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明的地震波初至检测方法及系统,利用卡方检验的方法识别地震波初至,关注地震记录数据落入某个区间的数量,能够避免依赖于地震数据本身的大小,从而有效避免噪声的影响,提高检测准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明地震波初至检测方法实施例的方法流程图;
图2为本发明地震波初至检测系统实施例的系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种地震波初至检测方法及系统,在低信噪比的情况下有效避免噪声对检测精度的影响,提高检测准确度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明地震波初至检测方法实施例的方法流程图。
本发明将时间窗口内采样值落入某一区间的次数作为指标,统计次数差值之和,并将此差值之和与卡方值相比较,将比较结果作为地震初至是否到达的标志。
参见图1,该地震波初至检测方法,包括:
步骤101:获取检波器采集的待检测的地震记录序列。
待检测的地震记录序列为p0,p1,…,pN-1,其中N为待检测的地震记录序列的长度,即待检测的地震记录序列所包含的地震数据的数量。
步骤102:建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数,且L满足31≤L≤41;所述第二时间窗紧随所述第一时间窗。
步骤103:对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据。
其中,对所述地震记录序列的左侧和右侧进行数据填充,具体包括:将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
步骤104:从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达。
所述第一时间窗内的数据序列为其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;pn为位于所述第一时间窗中心的数据;P表示第一时间窗的数据序列,p表示地震记录序列内的各个数据;由于第二时间窗紧随第一时间窗,则所述第二时间窗内的数据序列为
其中,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达,具体包括:
求取所述第一时间窗Pn内所有数据的最大值和最小值,得到取值区间;其中最小值记为p(n) min=min[Pn],最大值记为p(n) max=max[Pn],则取值区间为[p(n) min,p(n) max]。
将所述取值区间划分为K段,得到K段小区间。其中K的取值一般为10;各个小区间的序号为k,划分之后的每个小区间可以表示为:
k=1,2,…,K
n=1,2,…,N
计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k)。
计算数量差值比率:
根据参数K和置信度α(α一般取为95%),通过查询卡方表获取卡方值χ2(α,K-1);
判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;
若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,即γ(n)≥χ2(α,K-1),则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,即γ(n)<χ2(α,K-1),则确定在所述地震记录序列的第n个数据处没有地震波动到达。
步骤105:记录所有有地震波动到达的数据在所述地震记录序列中的序号nFB,并计算每个有地震波动到达的数据对应的初至。各个有地震波动到达的数据所对应的地震波初至为t=(nFB-1)*T。其中T为检波器采样时间间隔,单位为秒。
地震波初至检测所面临的最大问题是噪声对检测精度的不良影响。噪声超过一定程度,会造成检测精度极速下降。而本发明则考察地震记录数据落入区间的次数,这种统计方法可以在一定程度上避免噪声的影响。
因此,本发明是对相邻的两个数据序列(P和Q)的分布进行统计刻画,如果两者差值之和小于卡方值,则认为这两个数据序列(P和Q)具有相同的分布,没有地震波动到达;如果大于卡方值,则认为这两个数据序列(P和Q)的分布不同,意味着有地震波动到达。利用分布是否相同,可以在一定程度上避免噪声的影响,提高检测精度,并且计算方法十分简单,原理非常明确。
图2为本发明地震波初至检测系统实施例的系统结构图。
参见图2,该地震波初至检测系统,包括:
序列获取模块201,用于获取检波器采集的待检测的地震记录序列。
时间窗建立模块202,用于建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数;所述第二时间窗紧随所述第一时间窗。
数据填充模块203,用于对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据。
所述数据填充模块203,具体包括:序列左填充单元和序列右填充单元。其中序列左填充单元,用于将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;序列右填充单元,用于将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
地震波初至检测模块204,用于从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达。
所述地震波初至检测模块204,具体包括:
区间计算单元,用于求取所述第一时间窗内所有数据的最大值和最小值,得到取值区间;区间划分单元,用于将所述取值区间划分为K段,得到K段小区间;区间数据量计算单元,用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k);其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;k为小区间序号,k=1,2,…,K;比率计算单元,用于计算数量差值比率:卡方值查询单元,用于根据参数K和置信度α,通过查询卡方表获取卡方值χ2(α,K-1);判断单元,用于判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;地震初至确定单元,用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,则确定在所述地震记录序列的第n个数据处没有地震波动到达。
记录与初至计算模块205,用于记录所有有地震波动到达的数据在所述地震记录序列中的序号,并计算每个有地震波动到达的数据对应的初至。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明的地震波初至检测方法及系统,利用卡方检验的方法识别地震波初至,关注地震记录数据落入某个区间的数量,能够避免依赖于地震数据本身的大小,从而有效避免噪声的影响,提高提高检测准确度。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种地震波初至检测方法,其特征在于,包括:
获取检波器采集的待检测的地震记录序列;
建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数;所述第二时间窗紧随所述第一时间窗;
对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据;
从所述第一时间窗的中心对准填充前的所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准填充前的所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达;
记录所有有地震波动到达的数据在所述地震记录序列中的序号,并计算每个有地震波动到达的数据对应的初至;
所述对所述地震记录序列的左侧和右侧进行数据填充,具体包括:
将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;
将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
2.根据权利要求1所述的一种地震波初至检测方法,其特征在于,所述每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达,具体包括:
求取所述第一时间窗内所有数据的最大值和最小值,得到取值区间;
将所述取值区间划分为K段,得到K段小区间;
计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k);其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;k为小区间序号,k=1,2,…,K;
计算数量差值比率:
根据参数K和置信度α,通过查询卡方表获取卡方值χ2(α,K-1);
判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;
若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,则确定在所述地震记录序列的第n个数据处没有地震波动到达。
3.一种地震波初至检测系统,其特征在于,包括:
序列获取模块,用于获取检波器采集的待检测的地震记录序列;
时间窗建立模块,用于建立第一时间窗和第二时间窗,所述第一时间窗和所述第二时间窗的容量均为L个数据,其中L为奇数;所述第二时间窗紧随所述第一时间窗;
数据填充模块,用于对所述地震记录序列的左侧和右侧进行数据填充,使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时,所述第一时间窗和所述第二时间窗内均填充满L个数据;
地震波初至检测模块,用于从所述第一时间窗的中心对准填充前的所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗,直至所述第一时间窗的中心对准填充前的所述地震记录序列的最后一个数据,每滑动一个数据的步长,对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验,从而确定所述地震记录序列中的各个数据处是否有地震波动到达;
记录与初至计算模块,用于记录所有有地震波动到达的数据在所述地震记录序列中的序号,并计算每个有地震波动到达的数据对应的初至;
所述数据填充模块,具体包括:
序列左填充单元,用于将所述地震记录序列的左侧补充至少个数据;左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同;
序列右填充单元,用于将所述地震记录序列的右侧补充至少个数据;右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
4.根据权利要求3所述的一种地震波初至检测系统,其特征在于,所述地震波初至检测模块,具体包括:
区间计算单元,用于求取所述第一时间窗内所有数据的最大值和最小值,得到取值区间;
区间划分单元,用于将所述取值区间划分为K段,得到K段小区间;
区间数据量计算单元,用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n) P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n) Q(k);其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号;k为小区间序号,k=1,2,…,K;
比率计算单元,用于计算数量差值比率:
卡方值查询单元,用于根据参数K和置信度α,通过查询卡方表获取卡方值χ2(α,K-1);
判断单元,用于判断所述数量差值比率是否大于或等于所述卡方值,得到第一判断结果;
地震初至确定单元,用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值,则确定在所述地震记录序列的第n个数据处有地震波动到达;若所述第一判断结果表示所述数量差值比率小于所述卡方值,则确定在所述地震记录序列的第n个数据处没有地震波动到达。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810768697.5A CN108535778B (zh) | 2018-07-13 | 2018-07-13 | 一种地震波初至检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810768697.5A CN108535778B (zh) | 2018-07-13 | 2018-07-13 | 一种地震波初至检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108535778A CN108535778A (zh) | 2018-09-14 |
CN108535778B true CN108535778B (zh) | 2019-07-05 |
Family
ID=63488277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810768697.5A Expired - Fee Related CN108535778B (zh) | 2018-07-13 | 2018-07-13 | 一种地震波初至检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108535778B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110361779B (zh) * | 2019-07-14 | 2021-10-01 | 广东石油化工学院 | 一种基于卡方分布的微震事件检测方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104076392A (zh) * | 2014-05-28 | 2014-10-01 | 中国矿业大学(北京) | 基于网格搜索和牛顿迭代的微震震源定位联合反演方法 |
CN107526713A (zh) * | 2017-07-04 | 2017-12-29 | 北京航天易联科技发展有限公司 | 一种被动式太赫兹人体安检仪成像时间的确定方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6686738B2 (en) * | 2001-04-17 | 2004-02-03 | Baker Hughes Incorporated | Method for determining decay characteristics of multi-component downhole decay data |
-
2018
- 2018-07-13 CN CN201810768697.5A patent/CN108535778B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104076392A (zh) * | 2014-05-28 | 2014-10-01 | 中国矿业大学(北京) | 基于网格搜索和牛顿迭代的微震震源定位联合反演方法 |
CN107526713A (zh) * | 2017-07-04 | 2017-12-29 | 北京航天易联科技发展有限公司 | 一种被动式太赫兹人体安检仪成像时间的确定方法和装置 |
Non-Patent Citations (4)
Title |
---|
A method for microseismic event detection and P-phase picking;G-Akis Tselentis;《SEG San Antonio 2011 Annual Meeting》;20111231;第1638-1642页 |
A multi-window algorithm for real-time automatic detection and picking of P-phases of microseismic events;Zuolin Chen and Robert R. Stewart;《CREWES Research Report》;20061231;第1-9页 |
Strategy for automated analysis of passive microseismic data based on S-transform, Otsu’s thresholding, and higher order statistics;G-Akis Tselentis,et al.;《GEOPHYSICS》;20121231;第43-54页 |
基于STA/LTA 方法的微地震事件自动识别技术;段建华,等;《煤田地质与勘探》;20150228;第76-81页 |
Also Published As
Publication number | Publication date |
---|---|
CN108535778A (zh) | 2018-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108535773B (zh) | 一种微地震事件检测方法及系统 | |
RU2694621C1 (ru) | Способ и устройство для обработки сейсмических данных | |
CN105425292B (zh) | 一种油气预测方法及装置 | |
Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
WO2017024523A1 (zh) | 一种射线弹性参数的反演方法 | |
CN105388518A (zh) | 一种质心频率与频谱比联合的井中地震品质因子反演方法 | |
CN108535777B (zh) | 一种地震波初至检测方法及系统 | |
CN112147684B (zh) | 表征同沉积断层活动强度的方法和装置 | |
CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
CN109884709A (zh) | 一种基于面波旅行时层析的转换波静校正方法 | |
CN108897038B (zh) | 一种微地震事件检测方法及系统 | |
CN105549083B (zh) | 一种确定地震观测系统重复性整体量度的方法和装置 | |
CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
CN108254788A (zh) | 一种地震波初至拾取方法及系统 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN108535778B (zh) | 一种地震波初至检测方法及系统 | |
CN107576985B (zh) | 一种地震反演的方法和装置 | |
CN106353807B (zh) | 裂缝识别方法和装置 | |
Rodríguez‐Pradilla et al. | Automated microseismic processing and integrated interpretation of induced seismicity during a multistage hydraulic‐fracturing stimulation, Alberta, Canada | |
Bernardino et al. | A multicomponent Isabella anomaly: Resolving the physical state of the Sierra Nevada upper mantle from Vp/Vs anisotropy tomography | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
CN102778691B (zh) | 一种计算检波器组内静校正时差的方法 | |
Hildyard et al. | T pd, a damped predominant period function with improvements for magnitude estimation | |
CN103498658B (zh) | 一种伽马测井数据处理方法 | |
Finney | The reality of GSSPs |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190705 Termination date: 20200713 |