CN109143367B - 一种初至波自动拾取方法 - Google Patents
一种初至波自动拾取方法 Download PDFInfo
- Publication number
- CN109143367B CN109143367B CN201710506333.5A CN201710506333A CN109143367B CN 109143367 B CN109143367 B CN 109143367B CN 201710506333 A CN201710506333 A CN 201710506333A CN 109143367 B CN109143367 B CN 109143367B
- Authority
- CN
- China
- Prior art keywords
- seismic
- seed
- point
- area
- arrival
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000000203 mixture Substances 0.000 claims description 2
- 238000012937 correction Methods 0.000 description 7
- 230000003068 static effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
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
公开了一种初至波自动拾取方法。该方法包括:1)在地震记录的第一区域中确定一个种子点作为生长的起始点;2)对地震道进行顺序扫描,将所述种子点周围邻域中与其具有相同或者相似性质的点合并到第一区域中;3)寻找第一区域的边界,将其确定为初至时间点。本发明基于地震道的局部相似性形成了一种新的地震属性道,并将新的地震属性道与地震道融合形成融合道,然后将融合道作为输入,用区域生长的方法来确定地震波初至前的区域,从而确定出初至,该方法具有抗噪能力强,拾取精度高的特点。
Description
技术领域
本发明涉及油气及煤层气地震勘探与开发领域,更具体地,涉及一种初至波自动拾取方法。
背景技术
静校正是提高叠加剖面信噪比和垂向分辨率的一项关键技术。常规折射波静校正处理是建立在水平均匀层状介质假设前提下的。当勘探区域的地表起伏大,近地表厚度、速度变化复杂时,折射静校正就无法适用了。层析静校正能反演任意介质的近地表速度模型,是目前解决复杂地表静校正最理想的静校正方法,也是最常用的方法。而初至拾取又是层析静校正处理的第一步,它的准确与否会直接影响到层析反演结果。最初的地震波初至拾取是以人机交互的方式人工拾取的。因此初至拾取就成为一项费时而繁琐的工作。近年来,出现了许多自动初至拾取方法,例如相关法、能量比法、最大振幅法、分形维法以及神经网络法等。其中,能量比法以其简单实用而得到了广泛应用。但是在复杂近地表条件下,由于初至波能量特征、波形特征、相位特征变化大,又常常受噪声影响,初至波自动拾取方法的精度和效率仍然无法满足要求。因此有必要提出一种能够提高初至波自动拾取的准确性、稳定性和抗干扰能力的初至波拾取方法。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提出了一种初至波自动拾取方法,其能够克服现有评价方法的缺点和不足。
根据本发明的初至波自动拾取方法,可以包括以下步骤:
1)在地震记录的第一区域中确定一个种子点作为生长的起始点;
2)对地震道进行顺序扫描,将所述种子点周围邻域中与其具有相同或者相似性质的点合并到第一区域中;
3)寻找第一区域的边界,将其确定为初至时间点。
优选地,在步骤1)中,指定全局最小的初至时间之前的时间点为第一个区域的种子点。
优选地,在步骤1)中,通过以下公式确定第一个区域的种子点:
其中tseed,j为第j道种子点时间,xj为第j道炮检距,v为勘探区域近地表速度值。
优选地,在步骤2)中,通过以下公式判断种子点周围邻域的点是否与其具有相同或者相似性质:
其中i为种子点周围邻域的点,yi为i点的区域值,si为i点的地震数据,sth为设定的阈值,区域值为1表示该区域为第一区域,区域值为-1表示该区域为第二区域。
优选地,si为地震属性融合道的地震数据,其通过以下公式定义:
si=s1i·w1+s2i·w2+s3i·w3 (3)
其中,s1i为地震道的数据,s2i为一个地震属性道的数据,s3i为另一种地震属性道的数据,w1、w2、w3为加权值。
优选地,w1=1,w2=1,w3=1。
优选地,s2i为地震道的局部标准方差,s3i为地震道的局部互相关系数。
优选地,阈值sth通过以下公式定义:
sth=f1·s1seed+f2·s2seed+f3·s3seed (4)
其中,s1seed、s2seed、s3seed分别为所述种子点处s1i,s2i和s3i的值,f1,f2, f3为加权值。
优选地,f1取1~2,f2=10,f3=10。
优选地,地震道的局部标准方差Std(t)通过以下公式计算:
本发明提出一种地震波初至自动拾取方法,其基于地震道的局部相似性形成了一种新的地震属性道,并将新的地震属性道与地震道融合形成融合道,然后将融合道作为输入,用区域生长的方法来确定地震波初至前的区域,从而确定出初至,该方法具有抗噪能力强,拾取精度高的特点。
本发明的方法具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的,或者将在并入本文中的附图和随后的具体实施例中进行详细陈述,这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1为根据本发明的示例性实施方案的初至波自动拾取方法的流程图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
以下参考图1详细描述根据本发明示例性实施例的初至波自动拾取方法。
该方法主要包括:
步骤1:在地震记录的第一区域中确定一个种子点作为生长的起始点。
第一步,对于每个需要分割的区域找一个种子点作为生长的起始点。对于地震记录而言,可以事先浏览地震记录,然后指定全局最小的初至时间之前的时间点为第一个区域的种子点,也可以通过下式粗略的但是自适应地给出每一道第一个区域的种子点:
其中tseed,j为第j道种子点时间,xj为第j道炮检距,v为勘探区域近地表速度的一个粗略估计值。
由于地震记录是分为两个区域,因此只需要给定一个区域的种子点,然后分割出第一个区域,那么剩下的就是第二个区域。所以不用确定第二个区域的种子点。如果将第一个区域的值都赋值为1,而第二个区域的值都赋值为-1,那么种子点处的值就是1,而其余点处的值都是-1。
步骤2:对地震道进行顺序扫描,将所述种子点周围邻域中与其具有相同或者相似性质的点合并到第一区域中。
将种子点周围邻域中与种子点具有相同或者相似性质的点合并到种子点所在的区域中,即把这些点处的区域值也设置为1。
判断种子点周围邻域的点是否与种子点具有相同或者相似性质,是需要一个准则。本发明中,是给定一个阈值,即通过下式来划分区域:
其中i为种子点周围邻域的点,yi为i点的区域值,si为i点的地震数据,sth为设定的阈值,区域值为1表示该区域为第一区域,区域值为-1表示该区域为第二区域。
对于信噪比比较高的地震记录,将阈值设置为种子点处地震记录的1~2倍即可很好的划分出两个区域。但是对于信噪比低的地震记录,这种方法还有可能失败。为此,本发明进一步提出一种地震属性融合道,即几种地震属性道进行融合,融合方式可以采用如下形式:
si=s1i·w1+s2i·w2+s3i·w3 (3)
其中,s1i为地震道的数据,s2i为一个地震属性道的数据,s3i为另一种地震属性道的数据,w1、w2、w3为加权值。w1,w2,w3为加权值,这些权值可以都设置为1。而s2i和s3i可以是局部属性道,例如地震道的局部标准方差和地震道局部互相关系数。
当(2)式中的输入不是地震道,而是(3)式中的融合道时,(2)式中的阈值也要做相应的改变。取s1i为地震道,s2i和s3i分别是地震道的局部标准方差和地震道局部互相关系数的最大值,则此时阈值的设置如下:
sth=f1·s1seed+f2·s2seed+f3·s3seed (4)
其中,s1seed、s2seed、s3seed分别为所述种子点处s1i,s2i和s3i的值,f1,f2, f3为加权值。f1取1~2,而f2和f3则可以取10
地震道的局部标准方差可以按下式计算:
其中Nw为一个滑动窗口的长度,si为i点的地震数据,sw为窗口内地震数据的平均值,t为时间。
地震道局部互相关系数的最大值是两个长度相同,但是起始时间相隔一定距离的两个窗口内的地震道的互相关系时的最大值,通常起始时间的间隔取1 个采样点。
步骤3:寻找第一区域的边界,将其确定为初至时间点。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (8)
1.一种初至波自动拾取方法,其特征在于,所述方法包括:
1)在地震记录的第一区域中确定一个种子点作为生长的起始点;
2)对地震道进行顺序扫描,将所述种子点周围邻域中与其具有相同或者相似性质的点合并到第一区域中;
3)寻找第一区域的边界,将其确定为初至时间点;
在步骤2)中,通过以下公式判断种子点周围邻域的点是否与其具有相同或者相似性质:
其中i为种子点周围邻域的点,yi为i点的区域值,si为i点的地震数据,sth为设定的阈值,区域值为1表示该区域为第一区域,区域值为-1表示该区域为第二区域;
si为地震属性融合道的地震数据,其通过以下公式定义:
si=s1i·w1+s2i·w2+s3i·w3 (3)
其中,s1i为地震道的数据,s2i为一个地震属性道的数据,s3i为另一种地震属性道的数据,w1、w2、w3为加权值。
2.根据权利要求1所述的初至波自动拾取方法,其特征在于,在步骤1)中,指定全局最小的初至时间之前的时间点为第一个区域的种子点。
4.根据权利要求1所述的初至波自动拾取方法,其特征在于,w1=1,w2=1,w3=1。
5.根据权利要求1或4所述的初至波自动拾取方法,其特征在于,s2i为地震道的局部标准方差,s3i为地震道的局部互相关系数。
6.根据权利要求5所述的初至波自动拾取方法,其特征在于,阈值sth通过以下公式定义:
sth=f1·s1seed+f2·s2seed+f3·s3seed (4)
其中,s1seed、s2seed、s3seed分别为所述种子点处s1i,s2i和s3i的值,f1,f2,f3为加权值。
7.根据权利要求6所述的初至波自动拾取方法,其特征在于,f1取1~2,f2=10,f3=10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710506333.5A CN109143367B (zh) | 2017-06-28 | 2017-06-28 | 一种初至波自动拾取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710506333.5A CN109143367B (zh) | 2017-06-28 | 2017-06-28 | 一种初至波自动拾取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143367A CN109143367A (zh) | 2019-01-04 |
CN109143367B true CN109143367B (zh) | 2020-06-19 |
Family
ID=64805551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710506333.5A Active CN109143367B (zh) | 2017-06-28 | 2017-06-28 | 一种初至波自动拾取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143367B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111694053A (zh) * | 2019-03-14 | 2020-09-22 | 中国石油天然气股份有限公司 | 初至拾取方法及装置 |
CN110687592B (zh) * | 2019-09-25 | 2021-03-16 | 中石化石油工程技术服务有限公司 | 一种微地震震相识别初至拾取方法、设备和存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7391675B2 (en) * | 2004-09-17 | 2008-06-24 | Schlumberger Technology Corporation | Microseismic event detection and location by continuous map migration |
CN102313901B (zh) * | 2010-06-29 | 2014-01-15 | 中国石油天然气集团公司 | 一种初至波迭代拾取的方法 |
CN104570087B (zh) * | 2013-10-29 | 2017-05-10 | 中国石油天然气集团公司 | 利用瞬时属性进行地震数据初至及同相轴拾取的方法 |
CN106154332B (zh) * | 2015-05-13 | 2018-08-31 | 中国石油化工股份有限公司 | 一种井中微地震纵横波事件初至识别方法 |
CN106054251B (zh) * | 2016-06-20 | 2018-08-10 | 中国石油天然气集团公司 | 一种初至波拾取方法及装置 |
-
2017
- 2017-06-28 CN CN201710506333.5A patent/CN109143367B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109143367A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mucciarelli et al. | The stability of the horizontal-to-vertical spectral ratio of triggered noise and earthquake recordings | |
CN109884707B (zh) | 近地表分层时深曲线静校正方法 | |
CN109143367B (zh) | 一种初至波自动拾取方法 | |
CN104793247A (zh) | 近地表速度建模方法及建模装置 | |
CN104570076A (zh) | 一种基于二分法的地震波初至自动拾取方法 | |
CN105116448A (zh) | 一种转换波方位各向异性校正方法及装置 | |
CN104422959A (zh) | 一种检测储层边界的曲率属性方法 | |
CN106291681A (zh) | 地震层位自动追踪方法及装置 | |
CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
CN105093329A (zh) | 地震数据的低频信息的自动补偿方法及装置 | |
CN104570122A (zh) | 基于反射波的地震数据静校正方法及装置 | |
Yan et al. | Seismic imaging and velocity analysis for an Alberta Foothills seismic survey | |
CN116027419B (zh) | 利用地震波形特征识别与判定含气砂岩储层的方法 | |
CN105093327A (zh) | 多分量地震数据的矢量均值滤波方法 | |
CN109143364B (zh) | 初至时间拾取方法 | |
Zhang et al. | First break of the seismic signals in oil exploration based on information theory | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN109061737B (zh) | 一种基于合成地震记录的储层预测方法及装置 | |
CN112147699B (zh) | 结合过井剖面的vsp深度域层速度模型构建方法及系统 | |
Gong et al. | Automatic time picking of microseismic data based on shearlet-AIC algorithm | |
Hawkins et al. | Addressing anisotropy in 3-D prestack depth migration: A case study from the Southern North Sea | |
CN103837890B (zh) | 获取地震数据的方法及设备 | |
Qin et al. | First-arrival automatic picking based on improved energy ratio method and outlier detection theory | |
CN114296137B (zh) | 一种地震q场建立方法及系统 | |
CN110857997A (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 |