一种地震反射波联合成像方法及系统技术方案

本发明专利技术公开了一种地震反射波联合成像方法及系统，涉及陆上或海洋地震勘探领域。所述地震反射波联合成像方法包括：对获取到的初始地震数据进行处理，得到叠前共成像点地震数据；对叠前共成像点地震数据中的表面相关多次波进行压制处理，得到包含一次反射波的零偏移距地震数据；对叠前共成像点地震数据和零偏移距地震数据进行互相关处理，得到包含表面相关多次波的虚反射地震数据；对叠前共成像点地震数据和虚反射地震数据进行互相关处理，得到数据驱动动较正后的共成像点地震数据；基于数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。本发明专利技术能提高地震反射波的成像质量。本发明专利技术能提高地震反射波的成像质量。本发明专利技术能提高地震反射波的成像质量。

【技术实现步骤摘要】
一种地震反射波联合成像方法及系统


[0001]本专利技术涉及陆上或海洋地震勘探领域，特别是涉及一种地震反射波联合成像方法及系统。

技术介绍

[0002]常规的地震成像方法仅对地震一次反射波进行成像，而地震多次反射波作为来自于介质内部的反射信息常常被当成主要的叠前噪音，严重影响了地震一次反射波的成像质量和可靠性，破坏了目标反射波场的振幅和相位信息，增加了目标层位构造解释和物性反演的不确定性。因此，需要在常规的地震数据处理中需要对地震多次波进行压制。而多次波中又以表面相关多次波能量最强，尤其在常规海洋地震勘探中，由于海水层的存在，地震波会在海水层中来回反射多次，进而形成强能量的表面相关多次波，其严重影响了一次反射波的精准成像。因此，需要对表面相关多次波进行压制。但是表面相关多次波本身亦包含了来自于介质内部的反射信息，尤其是可以增加成像质量的广角反射信息，其可以弥补一次反射波能量的不足，大大增加反射波成像剖面的信噪比和分辨率。目前，常规的基于速度模型和一次反射波的成像方法不仅会产生信号畸变，而且不能有效利用表面相关多次波，地震数据的成像质量有待提高。

技术实现思路

[0003]基于此，本专利技术实施例提供一种地震反射波联合成像方法及系统，采用表面相关多次波与一次波联合成像，以提高地震反射波的成像质量。
[0004]为实现上述目的，本专利技术提供了如下方案：
[0005]一种地震反射波联合成像方法，包括：
[0006]获取初始地震数据；所述初始地震数据为加载有观测系统信息的地震数据；<br/>[0007]对所述初始地震数据进行处理，得到叠前共成像点地震数据；
[0008]对所述叠前共成像点地震数据中的表面相关多次波进行压制处理，得到包含一次反射波的零偏移距地震数据；
[0009]对所述叠前共成像点地震数据和所述零偏移距地震数据进行互相关处理，得到包含表面相关多次波的虚反射地震数据；
[0010]对所述叠前共成像点地震数据和所述虚反射地震数据进行互相关处理，得到数据驱动动较正后的共成像点地震数据；
[0011]基于所述数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。
[0012]可选的，所述对所述初始地震数据进行处理，得到叠前共成像点地震数据，具体包括：
[0013]对所述初始地震数据进行地震鬼波压制处理，得到鬼波压制后的地震数据；
[0014]利用所述鬼波压制后的地震数据，拾取初至波旅行时，并基于所述初至波旅行时、
折射静校正量求取算子和折射静校正量加载算子，得到折射静校正后的地震数据；
[0015]对所述折射静校正后的地震数据进行噪声压制，得到噪音压制后的地震数据；
[0016]对所述噪音压制后的地震数据进行宽频反褶积处理，得到宽频反褶积处理后的地震数据；
[0017]对所述宽频反褶积处理后的地震数据进行地表一致性振幅补偿，得到地表一致性振幅补偿后的地震数据；
[0018]对所述地表一致性振幅补偿后的地震数据进行叠前偏移处理，得到偏移后的地震叠前数据体；
[0019]对所述叠前偏移后的地震叠前数据体进行反动校正处理，得到叠前共成像点地震数据。
[0020]可选的，所述基于所述数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果，具体包括：
[0021]对所述数据驱动动较正后的共成像点地震数据进行加权移动平均滤波处理，得到滤波后的地震数据；所述滤波后的地震数据为包含目标反射波场的地震数据；
[0022]利用所述滤波后的地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。
[0023]可选的，所述零偏移距地震数据的计算公式为：
[0024]U
P
(0；ω)＝SRME(U
G
(x；ω))；
[0025]其中，U
P
(0；ω)表示零偏移距地震数据；SRME()表示表层相关多次波压制算子；U
G
(x；ω)表示叠前共成像点地震数据；x表示偏移距；ω表示角频率。
[0026]可选的，所述虚反射地震数据的计算公式为：
[0027]V
G
(x；ω)＝U
G
(x；ω)U
P*
(0；ω)；
[0028]其中，V
G
(x；ω)表示虚反射地震数据；U
G
(x；ω)表示叠前共成像点地震数据；U
P*
(0；ω)表示U
P
(0；ω)的复共轭；U
P
(0；ω)表示零偏移距地震数据；x表示偏移距；ω表示角频率。
[0029]可选的，所述数据驱动动较正后的共成像点地震数据的计算公式为：
[0030]R
G
(x；ω)＝U
G
(x；ω)V
G*
(x；ω)；
[0031]其中，R
G
(x；ω)表示数据驱动动较正后的共成像点地震数据；U
G
(x；ω)表示叠前共成像点地震数据；V
G*
(x；ω)表示V
G
(x；ω)的复共轭；V
G
(x；ω)表示虚反射地震数据；x表示偏移距；ω表示角频率。
[0032]可选的，所述滤波后的地震数据的计算公式为：
[0033]P(x；ω)＝M(R
G
(x；ω))；
[0034]其中，P(x；ω)表示滤波后的地震数据；M()表示加权移动平均滤波算子；R
G
(x；ω)表示数据驱动动较正后的共成像点地震数据；x表示偏移距；ω表示角频率。
[0035]可选的，所述成像结果为：
[0036][0037]其中，E(ω)表示成像结果；P(x；ω)表示滤波后的地震数据；x表示偏移距；ω表示角频率。
[0038]本专利技术还提供了一种地震反射波联合成像系统，包括：
[0039]数据获取模块，用于获取初始地震数据；所述初始地震数据为加载有观测系统信息的地震数据；
[0040]数据处理模块，用于对所述初始地震数据进行处理，得到叠前共成像点地震数据；
[0041]多次波压制处理模块，用于对所述叠前共成像点地震数据中的表面相关多次波进行压制处理，得到包含一次反射波的零偏移距地震数据；
[0042]第一互相关处理模块，用于对所述叠前共成像点地震数据和所述零偏移距地震数据进行互相关处理，得到包含表面相关多次波的虚反射地震数据；
[0043]第二互相关处理模块，用于对所述叠前共成像点地震数据和所述虚反射地震数据进行互相关处理，得到数据驱动动较正后的共成像点地震数据；
[0044]成像模块，用于基于所述数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。
[0045]可本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种地震反射波联合成像方法，其特征在于，包括：获取初始地震数据；所述初始地震数据为加载有观测系统信息的地震数据；对所述初始地震数据进行处理，得到叠前共成像点地震数据；对所述叠前共成像点地震数据中的表面相关多次波进行压制处理，得到包含一次反射波的零偏移距地震数据；对所述叠前共成像点地震数据和所述零偏移距地震数据进行互相关处理，得到包含表面相关多次波的虚反射地震数据；对所述叠前共成像点地震数据和所述虚反射地震数据进行互相关处理，得到数据驱动动较正后的共成像点地震数据；基于所述数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。2.根据权利要求1所述的一种地震反射波联合成像方法，其特征在于，所述对所述初始地震数据进行处理，得到叠前共成像点地震数据，具体包括：对所述初始地震数据进行地震鬼波压制处理，得到鬼波压制后的地震数据；利用所述鬼波压制后的地震数据，拾取初至波旅行时，并基于所述初至波旅行时、折射静校正量求取算子和折射静校正量加载算子，得到折射静校正后的地震数据；对所述折射静校正后的地震数据进行噪声压制，得到噪音压制后的地震数据；对所述噪音压制后的地震数据进行宽频反褶积处理，得到宽频反褶积处理后的地震数据；对所述宽频反褶积处理后的地震数据进行地表一致性振幅补偿，得到地表一致性振幅补偿后的地震数据；对所述地表一致性振幅补偿后的地震数据进行叠前偏移处理，得到偏移后的地震叠前数据体；对所述叠前偏移后的地震叠前数据体进行反动校正处理，得到叠前共成像点地震数据。3.根据权利要求1所述的一种地震反射波联合成像方法，其特征在于，所述基于所述数据驱动动较正后的共成像点地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果，具体包括：对所述数据驱动动较正后的共成像点地震数据进行加权移动平均滤波处理，得到滤波后的地震数据；所述滤波后的地震数据为包含目标反射波场的地震数据；利用所述滤波后的地震数据，采用水平叠加法得到目标反射波场中共反射点处的成像结果。4.根据权利要求1所述的一种地震反射波联合成像方法，其特征在于，所述零偏移距地震数据的计算公式为：U
P
(0；ω)＝SRME(U
G
(x；ω))；其中，U
P
(0；ω)表示零偏移距地震数据；SRME()表示表层相关多次波压制算子；U
G
(x；ω)表示叠前共成像点地震数据；x表示偏移距；ω表示角频率。5.根据权利要求1所述的一种地震反射波联合成像方法，其特征在于，所述虚反射地震数据的计算公式为：
V
G
(x；ω)＝U
G
(x；ω)U
P*
(0；ω)；其中，V
G
(x；ω)表示虚反射地震数据；U
G
(x；ω)表示叠前共成像点地震数据；U
P*
(0；ω)表示U
P
(0；ω)的复共轭；U
P
(0；ω)表示零偏移距地震数据；x表示偏移距；ω表示角频...

【专利技术属性】
技术研发人员：曹成寅，乔宝平，李子伟，潘自强，黄昱丞，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：