CN108387936B - 一种可控震源资料初至高效拾取方法及层析静校正方法 - Google Patents

一种可控震源资料初至高效拾取方法及层析静校正方法 Download PDF

Info

Publication number
CN108387936B
CN108387936B CN201810227396.1A CN201810227396A CN108387936B CN 108387936 B CN108387936 B CN 108387936B CN 201810227396 A CN201810227396 A CN 201810227396A CN 108387936 B CN108387936 B CN 108387936B
Authority
CN
China
Prior art keywords
arrival
picking
data
vibroseis
vibroseis data
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
Application number
CN201810227396.1A
Other languages
English (en)
Other versions
CN108387936A (zh
Inventor
杜佳骏
夏凌
李若禹
彭路
候莉敏
冯慧
石冷睿
庄重
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BGP Inc filed Critical BGP Inc
Priority to CN201810227396.1A priority Critical patent/CN108387936B/zh
Publication of CN108387936A publication Critical patent/CN108387936A/zh
Application granted granted Critical
Publication of CN108387936B publication Critical patent/CN108387936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/52—Move-out correction
    • 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)
  • Apparatus For Radiation Diagnosis (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明提供了一种可控震源资料初至高效拾取方法。所述方法包括：自动拾取可控震源资料的初至视速度，获得初至视速度变化趋势；根据初至视速度变化趋势确定控制点，人工拾取控制点的初至视速度，根据所述控制点的初至视速度，通过插值获得初至视速度数据体；根据初至视速度数据体对可控震源资料进行动校正，对经过动校正之后的可控震源资料完成初至自动拾取。根据本发明方法能够进行高效、高质量的初至自动拾取；能够获得高质量的层析静校正叠加剖面；能够很好地解决低频静校正的问题；能够在表层静校正处理中起着至关重要的作用。

Description

一种可控震源资料初至高效拾取方法及层析静校正方法
技术领域
本发明属于地震资料处理技术领域，更具体地讲，涉及一种利用空变和时变初至视速度拾取方法、智能相关能量比法和多核并行运算方式的可控震源资料初至高效拾取方法及层析静校正方法。
背景技术
自动检测初至时间的方法就是寻找地震道振幅能量的突然增大，可以通过计算地震道每一个采样点前后能量比率的最大值来实现，与此同时，从其它的地震道属性，例如，振幅、标准偏差、瞬时频率也可以计算得到类似的比率。然而，当地震资料的初至较弱或者信噪比较低时这些单道技术就会失效。因此，必须从典型地震记录中获得的多余信息来获得准确的初至时间，也就是说，拾取初至时间会利用单道和相邻道的信息。
炮相关自动拾取法综合利用了相关技术和地震属性比来拾取折射初至时间，拾取时间由炮相关叠加和道叠加得到。通过叠加可以提高拾取的信噪比，提高拾取的质量和一致性，然而，基本的假设是炮相关叠加和道叠加比单独拾取更简单和可靠。
可控震源资料初至的特点主要表现在，近偏移距初至杂乱，远、中、近偏移距初至视速度变化剧烈，中、远偏移距初至能量弱，这些因素导致初至难于识别。常规初至拾取方法采用的是常速线性动校时窗，不能适应可控震源资料远、中、近偏移距初至视速度变化剧烈的特点，对可控震源资料初至的自动拾取效果不理想。
发明内容
针对现有技术中存在的问题，本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的之一在于提供一种利用空变和时变初至视速度拾取方法、智能相关能量比法和多核并行运算方式的可控震源资料初至高效拾取方法。
例如，本发明的另一目的在于提供一种基于时变和空变线性动校速度初至自动拾取的层析静校正方法。
为了实现上述目的，本发明一方面提供了一种可控震源资料初至高效拾取方法。所述方法包括：
自动拾取可控震源资料的初至视速度，获得初至视速度变化趋势。
根据初至视速度变化趋势确定控制点，人工拾取控制点的初至视速度，根据所述控制点的初至视速度，通过插值获得初至视速度数据体。
根据初至视速度数据体对可控震源资料进行动校正，对经过动校正之后的可控震源资料完成初至自动拾取。
为了实现上述目的，本发明另一方面提供了一种基于时变和空变线性动校速度初至自动拾取的层析静校正方法。所述方法包括：采用可控震源资料初至高效拾取方法，完成可控震源资料的初至拾取；计算层析静校正量，获得层析静校正叠加剖面。
在本发明的一个示例性实施例中，所述初至视速度的拾取，二维可控震源资料可以沿检波线左半支和检波线右半支拾取，三维可控震源资料可以沿4个象限方向分别拾取。
在本发明的一个示例性实施例中，所述插值可以采用三角函数线性插值方法。
在本发明的一个示例性实施例中，所述控制点可以根据初至视速度变化趋势确定，初至视速度变化增大，控制点增多，初至视速度变化减小，控制点减少。
在本发明的一个示例性实施例中，所述初至自动拾取可以采用智能相关能量比法。
在本发明的一个示例性实施例中，所述初至高效拾取方法可以采用多核并行运算方式。
与现有技术相比，本发明的有益效果包括：能够进行高效、高质量的初至自动拾取；能够获得高质量的层析静校正叠加剖面；能够很好地解决低频静校正的问题；能够在表层静校正处理中起着至关重要的作用。
附图说明
图1示出了根据本发明的方法的一个示例性实施例的流程图。
图2示出了根据本发明的方法的人工拾取初至视速度示意图。
图3示出了根据本发明的方法的控制点选取示意图。
图4示出了根据常规初至自动拾取方法的单炮初至拾取示意图。
图5示出了根据本发明的方法的单炮初至拾取示意图。
图6示出了根据常规初至自动拾取方法计算得到的层析静校正叠加剖面示意图。
图7示出了根据本发明的层析静校正方法计算得到的层析静校正叠加剖面示意图。
具体实施方式
在下文中，将结合示例性实施例和附图来详细说明本发明的可控震源资料初至高效拾取方法及层析静校正方法。
图1示出了根据本发明的示例性实施例的可控震源资料初至高效拾取方法及层析静校正方法的流程图。如图1所示，在一个示例性实施例中，本发明的可控震源资料初至高效拾取方法及层析静校正方法可以通过如下步骤实现：
步骤(1)，首先，选取工区内一个合适的震源点来人工拾取初至视速度，其中，二维可控震源资料沿检波线左半支和检波线右半支两个方向人工拾取初至视速度。三维可控震源资料将360°方位角初至分成四个象限，沿4个象限方向分别时变拾取不同方向的初至视速度。如图2所示为三维可控震源资料的初至视速度人工拾取示意图，其中，横轴代表线道号，纵轴代表时间，图中的第1部分为沿第一象限方向拾取，如图黑色箭头所指的黑线；图中的第2部分为沿第三象限方向拾取，如图中黑色箭头所指的黑线；图中的第3部分为沿第二象限方向拾取，如图中黑色箭头所指的黑线；图中的第4部分为沿第四象限方向拾取，如图中黑色箭头所指的黑线。完成该震源点的初至拾取之后，通过自动拾取的方式完成全区可控震源资料的初至视速度拾取，获得全区的I类初至视速度数据体。最后，通过对I类初至视速度数据体在垂向方向进行加权计算，得到可控震源资料的初至视速度变化趋势，可以较好地反映全区的初至视速度变化。需要说明的是，初至视速度变化趋势反映视速度的变化剧烈程度，可以用于控制点的选择(在本发明中控制点为根据初至视速度变化趋势选取的用于人工拾取视速度的一个或多个震源点)。还需要注意的是，控制点的选取在视速度变化大的区域可以多选控制点，视速度稳定的区域控制点可以选取少一些。控制点选取的原则是人工拾取的初至视速度要能控制全区初至视速度变化趋势，如果全区的初至视速度比较稳定，可以只选取一个控制点。
步骤(2)，根据步骤(1)获得的初至视速度变化趋势选取合适的控制点，如图3所示，图中黑色八角星代表选取的控制点，其中，横轴代表X坐标，纵轴代表Y坐标。图中的不同颜色分布代表了初至视速度变化趋势，颜色较统一地方说明初至视速度较稳定，颜色分布差异较大说明工区不同位置初至视速度差异较大。在选取控制点的时候可以在颜色差异大的区域适当增加速度控制点密度。对选取的所有控制点进行人工拾取初至视速度，如图2所示，根据人工拾取的控制点的初至视速度，可以通过三角函数线性插值的方式获得II类初至视速度数据体，该II类初至视速度数据体可以更好地反映全区的视速度变化趋势。从步骤(1)到步骤(2)可以理解为是一个优化的过程，可以不断循环该过程以达到初至视速度拾取的最优化，为后续计算提供更加准确的数据。
步骤(3)，根据步骤(2)获得的II类初至视速度数据体对可控震源资料进行动校正，对经过动校正之后的可控震源资料可以通过智能相关能量比法完成初至自动拾取，如图4、图5所示，其中，横轴代表X坐标，纵轴代表时间。图4为采用常规方法完成的初至自动拾取，图5为采用本发明的方法完成的初至自动拾取，通过对比可以看出(如图中黑色箭头所指)，采用本发明的方法进行初至自动拾取的精度更高，例如，同相轴更加连续等。
步骤(4)，可以对采用本发明的方法完成初至自动拾取之后的可控震源资料进行层析静校正，获得时间域层析静校正叠加剖面，如图6、图7所示，其中，横轴代表横纵测线号，纵轴代表时间。图6为采用常规初至自动拾取方法之后计算得到的层析静校正叠加剖面，图7为采用本发明的方法完成初至自动拾取之后计算得到的层析静校正叠加剖面。通过对比可以看出(如图中的黑色椭圆圈内)，采用本发明的方法完成初至自动拾取之后计算得到的层析静校正量的应用效果更好，低频静校正问题可以得到很好地解决，例如，同相轴更加连续，能量更强等。
需要说明的是，在步骤(1)至步骤(4)的计算过程中，可以采用多核并行的方式，通过多核并行的计算方式可以提高计算效率。
综上所述，采用本发明的方法，能够进行高效、高质量的初至自动拾取；能够获得高质量的层析静校正叠加剖面；能够很好地解决低频静校正的问题；能够在表层静校正处理中起着至关重要的作用。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。

Claims (7)

1.一种可控震源资料初至高效拾取方法，其特征在于，所述方法包括以下步骤：
选取震源点并人工拾取所述选取震源点的初至视速度，自动拾取可控震源资料的初至视速度，获得初至视速度变化趋势；
根据初至视速度变化趋势确定控制点，人工拾取控制点的初至视速度，根据所述控制点的初至视速度，通过插值获得初至视速度数据体；
根据初至视速度数据体对可控震源资料进行动校正，对经过动校正之后的可控震源资料完成初至自动拾取。
2.根据权利要求1所述的可控震源资料初至高效拾取方法，其特征在于，所述初至视速度的拾取，二维可控震源资料沿检波线左半支和检波线右半支拾取，三维可控震源资料沿4个象限方向分别拾取。
3.根据权利要求1所述的可控震源资料初至高效拾取方法，其特征在于，所述插值采用三角函数线性插值方法。
4.根据权利要求1所述的可控震源资料初至高效拾取方法，其特征在于，所述控制点根据初至视速度变化趋势确定，初至视速度变化增大，控制点增多，初至视速度变化减小，控制点减少。
5.根据权利要求1所述的可控震源资料初至高效拾取方法，其特征在于，所述对经过动校正之后的可控震源资料完成初至自动拾取步骤中，所述初至自动拾取采用智能相关能量比法。
6.根据权利要求1所述的可控震源资料初至高效拾取方法，其特征在于，所述初至高效拾取方法采用多核并行的方式进行运算。
7.一种层析静校正方法，其特征在于，所述方法还包括以下步骤：
采用如权利要求1至6中任意一项所述的可控震源资料初至高效拾取方法，完成可控震源资料的初至拾取；
计算层析静校正量，获得层析静校正叠加剖面。
CN201810227396.1A 2018-03-20 2018-03-20 一种可控震源资料初至高效拾取方法及层析静校正方法 Active CN108387936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810227396.1A CN108387936B (zh) 2018-03-20 2018-03-20 一种可控震源资料初至高效拾取方法及层析静校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810227396.1A CN108387936B (zh) 2018-03-20 2018-03-20 一种可控震源资料初至高效拾取方法及层析静校正方法

Publications (2)

Publication Number Publication Date
CN108387936A CN108387936A (zh) 2018-08-10
CN108387936B true CN108387936B (zh) 2020-01-10

Family

ID=63067693

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810227396.1A Active CN108387936B (zh) 2018-03-20 2018-03-20 一种可控震源资料初至高效拾取方法及层析静校正方法

Country Status (1)

Country Link
CN (1) CN108387936B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109685127A (zh) * 2018-12-17 2019-04-26 郑州云海信息技术有限公司 一种并行深度学习初至拾取的方法和系统
CN112099090B (zh) * 2019-06-17 2021-07-30 中国海洋大学 地震资料视速度域非一致性长波长静校正方法
CN111060966B (zh) * 2019-12-31 2022-03-18 恒泰艾普（北京）能源科技研究院有限公司 一种用于静校正技术的初至拾取方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN103837891A (zh) * 2014-03-26 2014-06-04 宋维琪 微地震初至的高精度拾取方法
CN106054251A (zh) * 2016-06-20 2016-10-26 中国石油天然气集团公司 一种初至波拾取方法及装置
CN106814391A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 基于菲涅尔体层析反演的地面微地震事件定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN103837891A (zh) * 2014-03-26 2014-06-04 宋维琪 微地震初至的高精度拾取方法
CN106814391A (zh) * 2015-11-27 2017-06-09 中国石油化工股份有限公司 基于菲涅尔体层析反演的地面微地震事件定位方法
CN106054251A (zh) * 2016-06-20 2016-10-26 中国石油天然气集团公司 一种初至波拾取方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《VTI介质平移初至走时动校正双正参数反演》;张博等;《石油物探》;20120331;第51卷(第2期);全文 *
《可控震源地震勘探新方法研究与应用》;倪宇东;《中国博士学位论文全文数据库基础科学辑》;20130515;全文 *
《基于FCM的微地震初至自动拾取算法研究》;朱丹;《中国优秀硕士学位论文全文数据库工程可以1辑》;20170915;全文 *

Also Published As

Publication number Publication date
CN108387936A (zh) 2018-08-10

Similar Documents

Publication Publication Date Title
CN108387936B (zh) 一种可控震源资料初至高效拾取方法及层析静校正方法
CN101609160B (zh) 一种抗背景噪音强的地震数据初至波自动拾取方法
CN110956078B (zh) 一种电力线检测方法及装置
CN108830899B (zh) 一种基于几何特征的快速自适应角点结合的圆心检测方法
Trad et al. Fast and robust deblending using apex shifted Radon transform
CN114418953B (zh) 一种基于时间历程累计图像的声呐目标检测方法和系统
CN113093280B (zh) 基于相干函数控制的虚反射走时电缆等浮校正方法
CN104865570B (zh) 一种快速的动态规划检测前跟踪方法
CN104570076A (zh) 一种基于二分法的地震波初至自动拾取方法
CN104155694B (zh) 一种反射转换横波共检波点叠加剖面的剩余静校正方法
CN111538082B (zh) 一种地震波时频域初至自动拾取方法
CN109696705A (zh) 深海obs上下行波场成像数据校正方法及装置
CN104764407A (zh) 一种电缆护套厚度的精细测量方法
CN104932011A (zh) 一种拾取地震波初至的方法及装置
CN106842325B (zh) 利用检波点静校正迭代法提高初至波拾取效率及准确率的方法
CN109782341B (zh) 一种基于信息散度的地震数据属性提取方法
CN111624653A (zh) 地震速度自动拾取方法及装置
CN111596359B (zh) 一种倾斜地层宽线数据校正成像方法及装置
CN106569278B (zh) 一种多道相似相干速度谱计算方法
CN111781644B (zh) 浅中地层地震数据的线性干扰衰减方法及装置
CN106646614A (zh) 面向时差校正的叠合波分离方法
CN105116446B (zh) 一种用反射层位约束进行速度分析的方法
CN109884701A (zh) 地质体散射角度导向深度成像方法
CN111856557A (zh) 深度域合成地震记录制作方法及装置
CN111562616A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201116

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.