转换波共检波点叠加静校正方法
技术领域
本发明涉及石油天然气地震勘探领域,更具体地讲,涉及一种能够精确地获取转换波的检波点静校正量的基于迭代的最优转换波共检波点叠加静校正方法。
背景技术
近年来,多波勘探得到了越来越多的重视。特别是多波勘探在岩性分析、流体识别、裂缝检测、气云区成像等方面显示出其独特优势,在实际工业生产在得到了越来越多的重视。转换波处理与传统的纵波处理有很多不同,其关键步骤通常是静校正量的求取。转换波静校正量是由炮点的纵波静校正量和检波点的横波静校正量组成,转换波炮点静校正可以用传统方法求得,而检波点横波静校正则需要专门的方法进行求取,这是因为横波与纵波在近地表传播速度差异极大,相同的位置上横波静校正量通常比纵波静校正量大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、磁盘或磁带等)上被存储。将理解:存储装置和存储介质是适于存储包括指令的程序的机器可读存储器的实施例,其中,当所述指令被执行时,所述指令实施本发明的实施例。因此,本发明的实施例提供一种包括代码的程序和存储所述程序的机器可读存储器,其中,所述代码用于实施如上所述的根据本发明的示例性实施例的转换波共检波点叠加静校正方法。进一步地说,所述程序可经由任何介质(诸如通过有线或无线连接携带的通信信号)被电子地传输,并且本发明的实施例适当地包含所述程序。
虽然已经显示和描述了一些实施例,但是本领域技术人员应该理解,在不脱离本发明的原理和精神的情况下,可以对这些实施例进行修改,本发明的范围由权利要求及其等同物限定。