CN103487834A - 转换波共检波点叠加静校正方法 - Google Patents

转换波共检波点叠加静校正方法 Download PDF

Info

Publication number
CN103487834A
CN103487834A CN201310408061.7A CN201310408061A CN103487834A CN 103487834 A CN103487834 A CN 103487834A CN 201310408061 A CN201310408061 A CN 201310408061A CN 103487834 A CN103487834 A CN 103487834A
Authority
CN
China
Prior art keywords
wave
section
layer position
static correction
correction amount
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
Application number
CN201310408061.7A
Other languages
English (en)
Other versions
CN103487834B (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
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201310408061.7A priority Critical patent/CN103487834B/zh
Publication of CN103487834A publication Critical patent/CN103487834A/zh
Application granted granted Critical
Publication of CN103487834B publication Critical patent/CN103487834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种转换波共检波点叠加静校正方法,包括以下步骤:(a)在纵波和转换波共检波点叠加剖面上分段拾取满足预定条件的层位;(b)执行自动层拉平,并在设置的时窗内自动精细拾取相应的层位段;(c)计算拾取的层位段的速度比值;(d)以所述速度比值将纵波剖面上拾取的层位段投影到转换波剖面,并将转换波剖面上拾取的层位段与投影到转换波剖面上的层位段之间的时间差确定为检波点静校正量;(e)将确定的检波点静校正量应用到转换波剖面,以便对检波点进行静校正。

Description

转换波共检波点叠加静校正方法
技术领域
本发明涉及石油天然气地震勘探领域,更具体地讲,涉及一种能够精确地获取转换波的检波点静校正量的基于迭代的最优转换波共检波点叠加静校正方法。
背景技术
近年来,多波勘探得到了越来越多的重视。特别是多波勘探在岩性分析、流体识别、裂缝检测、气云区成像等方面显示出其独特优势,在实际工业生产在得到了越来越多的重视。转换波处理与传统的纵波处理有很多不同,其关键步骤通常是静校正量的求取。转换波静校正量是由炮点的纵波静校正量和检波点的横波静校正量组成,转换波炮点静校正可以用传统方法求得,而检波点横波静校正则需要专门的方法进行求取,这是因为横波与纵波在近地表传播速度差异极大,相同的位置上横波静校正量通常比纵波静校正量大2~10倍。与纵波相比,横波实际上不受近地表潜水面的影响,因此纵波和横波静校正量不存在简单的比例关系,使用纵波静校正量乘以某一系数来对横波静校正量进行近似通常是不能用于精细成像的。
检波点横波静校正方法可以分为两大类:一类是基于建立速度模型再求取静校正量的方法,另一类是基于数据直接得到静校正量的方法。模型法主要依赖于近地表速度模型。基于地震数据本身的静校正方法,即,共检波点叠加相干法,虽然在求取静校正量的方法上得到了最好的处理效果,但需要很多的手工操作,不利于规模应用,而且仅适用于构造较为简单的情况。
发明内容
因此,本发明的目的在于提供一种基于迭代的最优转换波共检波点叠加静校正方法,从而为转换波处理提供更精细和更高效的静校正处理方法。
根据本发明的一方面,提供一种转换波共检波点叠加静校正方法,包括以下步骤:(a)在纵波和转换波共检波点叠加剖面上分段拾取满足预定条件的层位;(b)执行自动层拉平,并在设置的时窗内自动精细拾取相应的层位段;(c)计算拾取的层位段的速度比值;(d)以所述速度比值将纵波剖面上拾取的层位段投影到转换波剖面,并将转换波剖面上拾取的层位段与投影到转换波剖面上的层位段之间的时间差确定为检波点静校正量;(e)将确定的检波点静校正量应用到转换波剖面,以便对检波点进行静校正。
优选地,所述转换波共检波点叠加静校正方法还包括以下步骤:在步骤(a)之前执行纵波和转换波共检波点叠加。
优选地,步骤(e)包括:(e1)将确定的检波点静校正量从低频至高频划分为多个频率成分;(e2)按照从低频至高频的顺序将划分的检波点静校正量的多个频率成分依次应用到转换波剖面,其中,在步骤(e2)中,每当应用划分的检波点静校正量的一个频率成分之后,判断检波点静校正量是否稳定,如果检波点静校正量不稳定,则应用划分的检波点静校正量的下一个频率成分,如果检波点静校正量稳定,则终止所述转换波共检波点叠加静校正方法,其中,如果在应用划分的检波点静校正量的最高频率成分之后,检波点静校正量仍然不稳定,则返回步骤(a),以进行下一轮的检波点静校正。
优选地,如果检波点静校正量的当前频率成分与检波点静校正量的前一个频率成分之间的差值大于预定阈值,则确定检波点静校正量不稳定。
优选地,所述预定条件表示通过拾取层位段获得的二维曲线段具有良好的连续性。
优选地,在步骤(c)中,根据如下等式来计算层位段的速度比值:
r(CRP)=2T0ps/T0pp-1,
其中,T0ps表示与转换波剖面上拾取的层位段对应的时间,T0pp表示与纵波剖面上拾取的层位段对应的时间,r(CRP)表示拾取的层位段的速度比值。
优选地,对层位段的速度比值r(CRP)进行空间上的平滑处理,来消除空间非构造因素的干扰。
附图说明
通过下面结合附图对实施例进行的描述,本发明的这些和/或其他方面和优点将会变得清楚和更易于理解,其中:
图1是示出根据本发明的示例性实施例的转换波共检波点叠加静校正方法的流程图。
具体实施方式
现在对本发明实施例进行详细的描述,其示例表示在附图中,其中,相同的标号始终表示相同部件。下面通过参照附图对实施例进行描述以解释本发明。
图1是示出根据本发明的示例性实施例的转换波共检波点叠加静校正方法的流程图。
参照图1,在步骤S101中,执行纵波和转换波共检波点叠加。这里,可根据现有技术容易地执行纵波和转换波共检波点叠加。
在步骤S102中,在纵波和转换波共检波点叠加剖面上分段拾取满足预定条件的层位,即,针对共检波点,同时在纵波剖面和转换波剖面上分段拾取相应的层位。这里,对于构造简单的地震剖面,一个层位就是地下某一时期沉积而形成的反射界面,在视觉上是从左到右的一个波峰或波谷连续拾取而获得的曲线,而由于沉积的复杂性、后期构造的破坏以及转换波采集技术的原因,通常一些层位连续性不佳。因此,本发明提出“分段”的概念,也就是拾取不同层位(时间)在空间上连续性较好的段。在这种情况下,层位段就是一个时间-空间的二维曲线段。所述预定条件表示通过拾取层位段获得的二维曲线段具有良好的连续性。因此,在步骤S102中,可以在纵波和转换波共检波点叠加剖面上拾取具有良好连续性的层位段。
此外,由于一个剖面在同一空间区域从浅到深会有很多连续的反射同相轴,因此每一个反射同相轴对应一个层位。在这种情况下,拾取的层位段可以不是同一层位。
接下来,在步骤S103中,执行自动层拉平,并在设置的时窗内自动精细拾取相应的层位段。由于在步骤S102中拾取的层位段一般都是扭曲的,特别是对于转换波来说,上下跳动比较严重。因此,可通过根据现有技术执行自动层拉平,通过整体上下移动拾取的地震道,使得拾取的层位点都在同一水平面上。通过层拉平以后,转换波本身上下跳动比较严重的情况可以得到改善。因此,在此基础上给定一个较小的时窗范围进行精细拾取,可以获得更好的成像效果,以便更好地进行后续处理。
其后,在步骤S104中,计算拾取的层位段的速度比值。具体地讲,可根据如下等式(1)来计算层位段的速度比值:
r(CRP)=2T0ps/T0pp-1             (1),
其中,T0ps表示与转换波剖面上拾取的层位段对应的时间,T0pp表示与纵波剖面上拾取的层位段对应的时间,r(CRP)表示空间-时间位置的速度比值,即,拾取的层位段的速度比值。这里,r(CRP)需要通过空间上的平滑处理,来消除空间非构造因素(主要是静校正)的干扰,以便后面计算的静校正更准确。
接着,在步骤S105中,以所述速度比值将纵波剖面上拾取的层位段投影到转换波剖面,并将转换波剖面上拾取的层位段与投影到转换波剖面上的层位段之间的时间差确定为检波点静校正量。
然后,在步骤S106中,将确定的检波点静校正量应用到转换波剖面,以便对检波点进行静校正。具体地讲,在步骤S106中,首先将确定的检波点静校正量从低频至高频划分为多个频率成分,然后按照从低频至高频的顺序将划分的检波点静校正量的多个频率成分依次应用到转换波剖面,以便对检波点进行静校正。这里,可通过横向平滑方法来提取检波点静校正量的低频成分,并通过从原始值减去经横向平滑得到的低频值来提取检波点静校正量的中高频成分。更具体地,每当应用划分的检波点静校正量的一个频率成分之后,需要判断检波点静校正量是否稳定。如果检波点静校正量不稳定,则应用划分的检波点静校正量的下一个频率成分,如果检波点静校正量稳定,则结束所述转换波共检波点叠加静校正方法。此外,如果在应用划分的检波点静校正量的最高频率成分之后,检波点静校正量仍然不稳定,则返回步骤S102,以进行下一轮的检波点静校正,从而进行更多的迭代处理。这里,如果检波点静校正量的当前频率成分与检波点静校正量的前一个频率成分之间的差值小于预定阈值,则可确定检波点静校正量稳定,如果所述差值大于预定阈值,则可确定检波点静校正量不稳定。
根据本发明的示例性实施例的转换波共检波点叠加静校正方法能够应用于二维和三维转换波高精度检波点静校正处理,对油气、矿产资源勘探有重要应用价值。此外,该方法特别适用于复杂山地的构造转换波静校正处理,具有计算步骤简单、计算效率高和处理效果好等特点,在目前西部的复杂山地地震资料处理中具有广阔的应用前景。
将理解:可按照硬件、软件或硬件和软件的组合的形式来实现如上所述的本发明的示例性实施例。可在非暂时性计算机可读存储介质中存储任何所述软件。非暂时性计算机可读存储介质存储一个或多个程序(软件模块),所述一个或多个程序包括指令,其中,当由在电子装置中的一个或多个处理器执行所述指令时,所述指令促使电子装置执行本发明的方法。任何所述软件可按照易失性或非易失性存储器(诸如,不论是否是可擦除的或可重写的存储装置,比如只读存储器(ROM))的形式被存储,或者按照存储器(诸如,随机存取存储器(RAM)、存储芯片、装置或集成电路)的形式被存储,或者在光学或磁可读介质(诸如,CD、DVD、磁盘或磁带等)上被存储。将理解:存储装置和存储介质是适于存储包括指令的程序的机器可读存储器的实施例,其中,当所述指令被执行时,所述指令实施本发明的实施例。因此,本发明的实施例提供一种包括代码的程序和存储所述程序的机器可读存储器,其中,所述代码用于实施如上所述的根据本发明的示例性实施例的转换波共检波点叠加静校正方法。进一步地说,所述程序可经由任何介质(诸如通过有线或无线连接携带的通信信号)被电子地传输,并且本发明的实施例适当地包含所述程序。
虽然已经显示和描述了一些实施例,但是本领域技术人员应该理解,在不脱离本发明的原理和精神的情况下,可以对这些实施例进行修改,本发明的范围由权利要求及其等同物限定。

Claims (7)

1.一种转换波共检波点叠加静校正方法,包括以下步骤:
(a)在纵波和转换波共检波点叠加剖面上分段拾取满足预定条件的层位;
(b)执行自动层拉平,并在设置的时窗内自动精细拾取相应的层位段;
(c)计算拾取的层位段的速度比值;
(d)以所述速度比值将纵波剖面上拾取的层位段投影到转换波剖面,并将转换波剖面上拾取的层位段与投影到转换波剖面上的层位段之间的时间差确定为检波点静校正量;
(e)将确定的检波点静校正量应用到转换波剖面,以便对检波点进行静校正。
2.根据权利要求1所述的转换波共检波点叠加静校正方法,还包括以下步骤:在步骤(a)之前执行纵波和转换波共检波点叠加。
3.根据权利要求1所述的转换波共检波点叠加静校正方法,其中,步骤(e)包括:
(e1)将确定的检波点静校正量从低频至高频划分为多个频率成分;
(e2)按照从低频至高频的顺序将划分的检波点静校正量的多个频率成分依次应用到转换波剖面,
其中,在步骤(e2)中,每当应用划分的检波点静校正量的一个频率成分之后,判断检波点静校正量是否稳定,如果检波点静校正量不稳定,则应用划分的检波点静校正量的下一个频率成分,如果检波点静校正量稳定,则终止所述转换波共检波点叠加静校正方法,
其中,如果在应用划分的检波点静校正量的最高频率成分之后,检波点静校正量仍然不稳定,则返回步骤(a),以进行下一轮的检波点静校正。
4.根据权利要求3所述的转换波共检波点叠加静校正方法,其中,如果检波点静校正量的当前频率成分与检波点静校正量的前一个频率成分之间的差值大于预定阈值,则确定检波点静校正量不稳定。
5.根据权利要求1所述的转换波共检波点叠加静校正方法,其中,所述预定条件表示通过拾取层位段获得的二维曲线段具有良好的连续性。
6.根据权利要求1所述的转换波共检波点叠加静校正方法,其中,在步骤(c)中,根据如下等式来计算层位段的速度比值:
r(CRP)=2T0ps/T0pp-1,
其中,T0ps表示与转换波剖面上拾取的层位段对应的时间,T0pp表示与纵波剖面上拾取的层位段对应的时间,r(CRP)表示拾取的层位段的速度比值。
7.根据权利要求6所述的转换波共检波点叠加静校正方法,其中,对层位段的速度比值r(CRP)进行空间上的平滑处理,来消除空间非构造因素的干扰。
CN201310408061.7A 2013-09-09 2013-09-09 转换波共检波点叠加静校正方法 Active CN103487834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310408061.7A CN103487834B (zh) 2013-09-09 2013-09-09 转换波共检波点叠加静校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310408061.7A CN103487834B (zh) 2013-09-09 2013-09-09 转换波共检波点叠加静校正方法

Publications (2)

Publication Number Publication Date
CN103487834A true CN103487834A (zh) 2014-01-01
CN103487834B CN103487834B (zh) 2016-07-13

Family

ID=49828203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310408061.7A Active CN103487834B (zh) 2013-09-09 2013-09-09 转换波共检波点叠加静校正方法

Country Status (1)

Country Link
CN (1) CN103487834B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104199103A (zh) * 2014-08-15 2014-12-10 中国石油天然气集团公司 转换波静校正方法及装置
CN104635270A (zh) * 2015-02-16 2015-05-20 中国地质大学(北京) 基于合成记录约束的转换波静校正方法
CN106443791A (zh) * 2016-11-11 2017-02-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 求取倾斜地层或各向异性地层横波剩余静校正量的方法
CN111751876A (zh) * 2020-07-07 2020-10-09 中油奥博(成都)科技有限公司 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置
CN111965701A (zh) * 2019-05-20 2020-11-20 中国石油天然气集团有限公司 近地表结构反演方法及系统
CN112946745A (zh) * 2019-12-11 2021-06-11 中国石油天然气集团有限公司 一种基于引力拟合的转换波剩余静校正方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5157638A (en) * 1992-01-13 1992-10-20 Conoco Inc. Method of deriving statics corrections from common reflection point gathers
CN1773310A (zh) * 2004-11-09 2006-05-17 中国石油天然气集团公司 转换波延迟时静校正方法
CN101419292A (zh) * 2007-10-25 2009-04-29 中国石油天然气集团公司 采用纵波源多分量地震数据生成横波地震剖面的方法
CN102692645A (zh) * 2012-06-01 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用纵波、转换波数据联合反演储层纵横波速度比的方法
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法
CN102798888A (zh) * 2011-05-27 2012-11-28 中国石油天然气集团公司 一种利用非零井源距数据计算纵横波速度比的方法
CN102854526A (zh) * 2011-06-30 2013-01-02 中国石油化工股份有限公司 一种多分量地震资料的处理方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5157638A (en) * 1992-01-13 1992-10-20 Conoco Inc. Method of deriving statics corrections from common reflection point gathers
CN1773310A (zh) * 2004-11-09 2006-05-17 中国石油天然气集团公司 转换波延迟时静校正方法
CN101419292A (zh) * 2007-10-25 2009-04-29 中国石油天然气集团公司 采用纵波源多分量地震数据生成横波地震剖面的方法
CN102798888A (zh) * 2011-05-27 2012-11-28 中国石油天然气集团公司 一种利用非零井源距数据计算纵横波速度比的方法
CN102854526A (zh) * 2011-06-30 2013-01-02 中国石油化工股份有限公司 一种多分量地震资料的处理方法
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法
CN102692645A (zh) * 2012-06-01 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用纵波、转换波数据联合反演储层纵横波速度比的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
彭苏萍: "用共检波点时窗道集做转换波静校正", 《中国地球物理学会年刊2002-中国地球物理学会第十八届年会论文集》 *
童庆: "转换波CRP叠加剖面构造静校正方法研究", 《成都理工大学硕士学位论文》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104199103A (zh) * 2014-08-15 2014-12-10 中国石油天然气集团公司 转换波静校正方法及装置
CN104199103B (zh) * 2014-08-15 2017-01-25 中国石油天然气集团公司 转换波静校正方法及装置
CN104635270A (zh) * 2015-02-16 2015-05-20 中国地质大学(北京) 基于合成记录约束的转换波静校正方法
CN104635270B (zh) * 2015-02-16 2017-03-08 中国地质大学(北京) 基于合成记录约束的转换波静校正方法
CN106443791A (zh) * 2016-11-11 2017-02-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 求取倾斜地层或各向异性地层横波剩余静校正量的方法
CN106443791B (zh) * 2016-11-11 2018-07-10 中国石油集团东方地球物理勘探有限责任公司 求取倾斜地层或各向异性地层横波剩余静校正量的方法
CN111965701A (zh) * 2019-05-20 2020-11-20 中国石油天然气集团有限公司 近地表结构反演方法及系统
CN111965701B (zh) * 2019-05-20 2023-09-26 中国石油天然气集团有限公司 近地表结构反演方法及系统
CN112946745A (zh) * 2019-12-11 2021-06-11 中国石油天然气集团有限公司 一种基于引力拟合的转换波剩余静校正方法及系统
CN111751876A (zh) * 2020-07-07 2020-10-09 中油奥博(成都)科技有限公司 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置
CN111751876B (zh) * 2020-07-07 2022-05-20 中油奥博(成都)科技有限公司 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置

Also Published As

Publication number Publication date
CN103487834B (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
CN103487834A (zh) 转换波共检波点叠加静校正方法
RU2760102C2 (ru) Способ и система автоматического распознавания центра залежи в карстовой пещере
CN107884822A (zh) 一种提高矿用微震震源定位精度的方法
CN105093319B (zh) 基于三维地震数据的地面微地震静校正方法
CN107783186B (zh) 一种近地表速度建模方法及装置
CN104181599A (zh) 一种基于近地表层的折射波静校正处理方法以及系统
RU2631407C1 (ru) Способ и устройство для обработки сейсмических сигналов
CN102692648A (zh) 一种基于大炮初至的折射波剩余静校正方法
CN102879819A (zh) 保持地震波场运动学特征的地震数据处理方法及装置
CN104375175A (zh) 倾角传播法层位自动追踪方法
CN107831542A (zh) Ddw高精度深度域井震匹配方法
CN102879820A (zh) 基于三角网格的三维表层模型构建方法
CN103941286A (zh) 一种基于平面导航的快速三维层位解释方法
CN104199103A (zh) 转换波静校正方法及装置
CN105301636A (zh) 速度模型的建立方法和装置
CN104330827A (zh) 一种地表模型静校正的处理方法
CN104155694A (zh) 一种反射转换横波共检波点叠加剖面的剩余静校正方法
CN106483537A (zh) 一种卫星定位轨迹的优化方法
CN104182651B (zh) 用于三分量检波器接收的微地震事件方位角自动质控方法
CN105301638B (zh) 一种提取风化层底界面的方法和装置
CN103576187B (zh) 一种利用反射成像时差确定多波静校正量的方法
CN109725354B (zh) 各向异性速度建模方法及系统
CN104459789A (zh) 一种地震波初至拾取的方法
CN105160156A (zh) 一种融合数据预处理的垂测电离图反演方法
CN105184039A (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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180206

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200916

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

Co-patentee after: BGP Inc., China National Petroleum Corp.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

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