基于归一化互相关的三维Marchenko成像方法技术

本申请提出了基于归一化互相关的三维Marchenko成像方法。根据本申请的技术方案以三维速度模型和三维工区地震资料作为输入，以迭代或直接反演的方法对Marchenko方程进行数值求解得到上下行格林函数，显著提高了成像精度。同时，由于本申请采用了选取优势频带及局部孔径地震数据进行格林函数重构并成像，降低了内存及计算量，在保证计算精度的同时，提高了计算效率。了计算效率。了计算效率。

【技术实现步骤摘要】
基于归一化互相关的三维Marchenko成像方法


[0001]本专利技术属于地震资料叠前深度偏移成像领域，更具体地，涉及基于归一化互相关的三维Marchenko成像方法。

技术介绍

[0002]多次波(尤其是层间多次波)是偏移成像难点。常规偏移成像方法以单散射假设为前提，利用一次反射波进行成像，无法直接对多次波进行处理，多次波在成像剖面通常表现为虚假构造或噪音干扰，此类成像方法在成像前需要对地震道集进行多次波压制，然后进行偏移成像处理。Marchenko成像方法是一种数据驱动的偏移成像方法，既可实现一次反射波偏移归位，又可有效处理层间多次波。目前，业界主要集中在二维Marchenko成像方法的研究与应用。由于二维测线仅仅是对地下构造某一方向切片显示，二维地震资料处理的前提是所有的地震波来自于同一个垂面，当地下构造不止分布于同一垂直平面内时，二维地震数据中会出现侧面波信息。常规的二维地震资料处理方法难以有效处理侧面波信息。

技术实现思路

[0003]有鉴于此，本申请提出了一种基于归一化互相关的三维Marchenko成像方法，具有更本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于归一化互相关的三维Marchenko成像方法，其特征在于，所述方法包括：步骤1，以初始格林函数G
d
(x
i
，x0，
‑
t)作为输入，并利用逆时变换得到初始上行聚焦函数其中，x0、x
i
，均表示三维位置坐标矢量，t表示时间；步骤2，将初始上行聚焦函数和校正后的地震数据作为输入，基于聚焦函数和格林函数通过迭代反演得到三维Marchenko聚焦函数；步骤3，选择指定频带和指定偏移距孔径内的地震数据和三维Marchenko聚焦函数，得到所述指定频带和指定偏移距孔径的上行格林函数G
‑
(x
i
，x0，ω)和下行格林函数G
+
(x
i
，x0，ω)，然后做傅里叶变换得到G
‑
(x
i
，x0，t)和G
+
(x
i
，x0，t)；步骤4，基于下式对初始格林函数和上行格林函数进行归一化互相关，得到三维Marchenko成像结果：其中，G
‑
(x
i
，x0，t)表示上行格林函数。2.根据权利要求1所述的方法，其特征在于，所述步骤1具体包括：基于下式得到地表激发点x0到地下聚焦点x
i
之间的投射响应的逆之间的投射响应的逆基于下式得到初始上行聚焦函数基于下式得到初始上行聚焦函数其中，x
′0表示三维位置坐标矢量，R(x0，x
′0，t)表示x
′0处激发、x0处接收得到的地表反射响应，D0表示积分空间，*表示褶积运算。3.根据权利要求2所述的方法，其特征在于，所述步骤2具体包括：基于下式进行迭代反演：演：演：其中，k表示迭代次数，表示第k
‑
1次迭代的上行聚焦函数，表示第k+1次迭代的上行聚焦函数，表示散射尾波。4.根据权利要求3所述的方法，其特征在于，在步骤3中，基于下式得到所述指定频带和指定偏移距孔径的上行格林函数G
‑
(x
i
，x0，ω)和下行格林函数G
+
(x
i
，x0，ω)：其中，ω表示频率，取所述指定频带内的数值；表示三维Marchenko上行聚焦函数；表示三维Marchenko下行聚焦函数。5.一种电子设备，其特征在于，所述电子设备包括：存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1
‑
4中任一项所述的方法。6.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1

【专利技术属性】
技术研发人员：韩冬，许璐，崔月，王鸣谦，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：