一种分频叠前时间偏移方法及系统技术方案

本发明专利技术提出了一种分频叠前时间偏移方法及系统，属于地震勘探领域。该方法包括：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。通过本发明专利技术的方法，有效压制了常规地震资料分频后所产生的“信号震荡”现象，且分频地震资料包含更丰富的地质信息。实施例证明了该专利技术的有效性和适用性。

A Frequency Division Prestack Time Migration Method and System

The invention provides a frequency division prestack time migration method and system, which belongs to the field of seismic exploration. The method includes frequency division of prestack CMP gathers and prestack time migration of prestack CMP gathers with different frequency bands. The method of the invention effectively suppresses the \signal oscillation\ phenomenon caused by frequency division of conventional seismic data, and the frequency division seismic data contains more abundant geological information. The embodiment demonstrates the validity and applicability of the patented technology.

【技术实现步骤摘要】
一种分频叠前时间偏移方法及系统
本专利技术属于石油地球物理勘探技术，具体涉及一种分频叠前时间偏移方法及系统。
技术介绍
基于地震资料的分频技术一直是地震资料解释的研究热点。将地震资料变换到不同的频带内，可以挖掘更多的地质信息和特征。常规的分频技术均是基于叠后地震资料开展的分频，地震剖面上常出现“信号震荡”的现象且分频信号的有效频带较窄。与叠后地震资料相比，叠前地震资料包含的有效信息较多，因此，叠前分频技术也会保留较多的有效信息。针对现有技术的缺点和不足，为此本专利专利技术了一种分频叠前时间偏移的方法及系统，用于得到质量更高、信息更丰富的分频地震资料。
技术实现思路
由于常规叠后地震资料的分频技术常出现“信号震荡”现象且分频信号的有效频带较窄，为此，本专利专利技术一种分频叠前时间偏移方法及系统，用于压制分频后产生的“信号震荡”现象，增加分频地震资料的有效信息。本专利技术首先对叠前CMP道集进行分频处理，得到不同频带的叠前CMP道集，结合全频带的叠前时间偏移处理流程，进一步开展分频带的偏移处理，最终得到不同频带的叠后地震数据体。根据本专利技术的一个方面，提供一种分频叠前时间偏移方法，该方法包括：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。进一步地，对叠前CMP道集进行分频处理包括：通过傅里叶变换将叠前CMP道集变换到频率域；设叠前CMP道集为输入数据s(t)，通过傅里叶变换将其变换到频率域，其中，s(w)为s(t)的振幅谱，e-iwt为基函数。进一步地，将变换到频率域的叠前CMP道集与Butterworth滤波器进行相乘，得到不同频带范围内的CMP道集振幅谱，再通过傅里叶反变换则得到不同频带范围内的叠前CMP道集。进一步地，所述Butterworth滤波器公式如下：其中，n为滤波器阶数，wc截止频率，w0中间频率，w频带分布范围。进一步地，所述不同频带范围内的叠前CMP道集为：其中，f(t)为不同频带的叠前CMP道集。根据本专利技术的另一个方面，提供一种分频叠前时间偏移系统，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。在对叠前CMP道集进行分频处理之后，结合已有的Kirchhoff叠前时间偏移处理流程对其进行分频带偏移，最终得到不同频带的叠后数据体。Kirchhoff叠前时间偏移处理流程为现有的方法，在此不再赘述。相对于已有技术，本专利通过一种分频叠前时间偏移方法，有效压制了分频过程中产生的“信号震荡”现象，同时增加了分频信号的有效信息。附图说明通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了分频叠前时间偏移处理流程图。图2示出了本专利技术一个实施例的叠前时间偏移处理流程。图3示出了不同滤波器的频率响应特征曲线比较图。图4示出了不同频带的叠前CMP道集，左为过A井的叠后全频带地震剖面，右上为过A井的全频带叠前CMP道集，右下分别为10-15Hz，30-35Hz，55-60Hz的叠前CMP道集。图5示出了叠前分频与叠后分频偏移剖面对比图。图6示出了不同频带分频偏移剖面对比图，左上为全频带地震剖面，右上为10-15Hz地震剖面，左下为30-35Hz地震剖面，右下为55-66Hz地震剖面。具体实施方式下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。本专利提供了一种频叠前时间偏移方法及系统，属于地震勘探领域。本方法包括：(1)搭建合理的叠前时间偏移处理流程；(2)构建Butterworth滤波器；(3)对叠前CMP道集进行分频处理；(4)对不同频带的叠前CMP道集进行叠前时间偏移处理。通过本专利的方法，有效压制了常规地震资料分频后所产生的“信号震荡”现象，且分频地震资料包含更丰富的地质信息。本公开提出了一种分频叠前时间偏移方法，该方法包括：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。在本方法中，首先对经过预处理后的叠前CMP道集进行分频处理，设叠前CMP道集为输入数据s(t)，通过傅里叶变换将其变换到频率域，其中：s(w)为s(t)的振幅谱，e-iwt为基函数。然后设计Butterworth滤波器(2)，公式如下：其中：n为滤波器阶数，wc截止频率，w0中间频率，w频带分布范围。将叠前CMP道集与Butterworth滤波器进行相乘，得到不同频带范围内的CMP道集振幅谱，通过傅里叶反变换则得到不同频带范围内的叠前CMP道集。其中：f(t)为不同频带的叠前CMP道集。在对叠前CMP道集进行分频处理之后，结合已有的Kirchhoff叠前时间偏移处理流程对其进行分频带偏移，最终得到不同频带的叠后数据体。根据本专利技术的另一个实施方式，提供一种分频叠前时间偏移系统，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。为便于理解本专利技术实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本专利技术，其任何具体细节并非意在以任何方式限制本专利技术。通过新疆某工区的实际资料来说明分频叠前时间偏移的有效性，该方法的技术流程如图1所示，图1示出了分频叠前时间偏移处理流程图。对叠前全频道集进行分频处理获得多个频带的叠前道集，然后进行PSTM(即Kirchhoff叠前时间偏移处理)，最后获得不同频带的叠后数据体。具体地，首先建立合适的偏移处理流程，如图2所示，图2示出了本专利技术一个实施例的叠前时间偏移处理流程。通过叠前时间偏移处理得到质量较好的叠前CMP道集。图2所示是现有的处理流程，第一步是对原始地震资料进行质量分析与评价，以设计合理的处理流程；第二步为层析静校正，消除复杂起伏地表对地震资料造成的影响；第三步为叠前保幅、保真去噪，以消除不同类型的噪音对地震资料成像的影响；第四步为子波一致性反褶积，以提高地震资料的分辨率；第五步为分频剩余静校正，进一步消除复杂起伏地表对地震资料造成的影响；最终得到较为可靠的叠前CMP道集，为分频处理和叠前时间偏移处理做好数据准备。再利用Butterworth滤波器对经过预处理后的叠前CMP道集进行分频处理，与其他滤波器相比，该滤波器具有较光滑的边界和较平坦的带宽(如图3所示)，可有效压制“信号震荡”现象。以过A井的叠前CMP道集为例进行分频处理，得到不同频带的叠前CMP道集(如图4所示)，再结合已有的偏移处理流程，进行分频带的叠前时间偏移处理，最终得到不同频带的叠后地震数据体(如图5-6所示)。由图5可见，由分频叠前时间偏移产生的30-35Hz地震剖面与全频带的地震剖面较相似，且剖面上的“震荡现象”明显弱于某软件分本文档来自技高网...

【技术保护点】
1.一种分频叠前时间偏移方法，其特征在于，该方法包括：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。

【技术特征摘要】
1.一种分频叠前时间偏移方法，其特征在于，该方法包括：对叠前CMP道集进行分频处理；对不同频带的叠前CMP道集进行叠前时间偏移处理。2.根据权利要求1所述的分频叠前时间偏移方法，其特征在于，对叠前CMP道集进行分频处理包括：通过傅里叶变换将叠前CMP道集变换到频率域；设叠前CMP道集为输入数据s(t)，通过傅里叶变换将其变换到频率域，其中，s(w)为s(t)的振幅谱，e-iwt为基函数。3.根据权利要求2所述的分频叠前时间偏移方法，其特征在于，将变换到频率域的叠前CMP道集与Butterworth滤波器进行相乘，得到不同频带范围内的CMP道集振幅谱，再通过傅里叶反变换则得到不同频带范围内的叠前CMP道集。4.根据权利要求3所述的分频叠前时间偏移方法，其特征在于，所述Butterworth滤波器公式如下：其中，n为滤波器阶数，wc截止频率，w0中间频率，w频带分布范围。5.根据权利要求4所述的分频叠前时间偏移方法，其特征在于，所述不同频带范围内的叠前CMP道集为：其中，f(t)为不同频带的叠前CMP道集。6.一种分频叠前时间偏移系统，其特征在于，该系统包...

【专利技术属性】
技术研发人员：许凯，唐金良，王世星，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11