2. "陆相页岩油气成藏及高效开发"教育部重点实验室, 大庆 163318
2. Key Laboratory of Continental Shale Hydrocarbon Accumulation and Efficient Development, Northeast Petroleum University, Daqing 163318, China
对于实际地震资料来说,多次反射的存在可能会使偏移成像结果产生强的假象,对后续地震资料的速度分析和解释造成一系列问题(黄饶等,2009;郭梦秋等,2019).多次波是指向下反射一次或多次的地震波,包括自由表面多次波和层间多次波.表面多次波是在自由表面上产生下行反射的地震波,在自由表面以下向下反射的多次波称为层间多次波.目前,多次波压制方法可分为滤波法和基于波动方程的预测相减法两大类.滤波法利用多次波的周期性和可分性来区分一次波和多次波,通常要求多次波和一次波之间有良好的可区分性.基于波动方程的预测相减法是一种利用波动方程预测多次波的数据驱动方法,将预测的多次波通过自适应匹配滤波从地震记录中减去,达到压制多次波的目的.由于层间多次波与一次有效波在传播时间和速度上的差异较小等因素,层间多次波的识别和压制难度较大,层间多次波衰减是当今国内外地球物理学专家和地震资料处理工作者所面临的难题.
层间多次波的压制方法不像自由表面多次波的压制方法(Verschuur et al., 1992;王维红等,2007;石颖等,2013)发展成熟,工业上还未形成较为系统的层间多次波压制技术.目前,对于层间多次波的压制研究主要为预测和减去层间多次波的方法,如Berkhout和Verschuur(1997)提出的层间多次波压制法,Weglein等(1997)提出的逆散射级数法(Inverse Scattering Series method,简称ISS),以及Jakubowicz(1998)提出的地表数据驱动法等.对于Berkhout和Verschuur提出的层间多次波压制方法是与自由表面相关多次波衰减(Surface Related Multiple Elimination,简称SRME)算法的扩展,将波场外推到产生层间多次波的界面,并采用SRME方法对多次波进行衰减,该算法需要了解地表和多次波产生界面之间的平滑速度模型.Berkhout和Verschuur(2005)在SRME方法的基础上,通过利用共聚焦点(Common Focus Point,简称CFP)道集使得SRME方法应用于地下散射点,实现层间多次波的预测.然而,该方法仍需要识别输入数据中发生层间多次反射的边界或层位的模型信息,需要复杂的基准面重建.对于ISS,利用自由表面、点散射以及地震波在参考介质中的传播预测多次波,Löer等(2016)基于该理论提出层间多次波可以用近似振幅进行预测.该方法虽然不需要任何的先验背景信息,但面临计算量大的问题.金德刚等(2008)推导了1.5D时空域ISS层间多次波预测算法,提高了计算效率,但不适用于三维区域.Jakubowicz(1998)利用逆时数据和一次反射地震数据预测相应的层间多次波,提出了一种显式的包括两道卷积和单道相关的层间多次波预测方法.目前的层间多次波预测方法大多基于这一思想.然而,这些方法需要精确地拾取地下层位,并且依赖于自适应减法来压制地震数据的层间多次波.吴静等(2013)将该方法扩展到多个界面产生的层间多次波衰减.但是该方法仍不能较为准确的预测出层间多次波的振幅.总体来说,反射界面选择和自适应减法是不可避免的.如果速度模型不准确,界面选择会受到影响.自适应减法通常是在最小二乘意义下执行的,最小能量准则的限制可能会对一次反射波造成损伤.
当前,国际上新发展起来一种针对层间多次波的衰减方法——基于Marchenko理论的层间多次波压制策略.该方法也是一种基于波动方程的压制技术,但不需要识别层间多次波产生界面.Broggini和Snieder(2012)将Marchenko理论引入地球物理领域.该理论的核心是一组方程耦合的Marchenko方程,这些方程可以迭代求解所谓的聚焦函数.聚焦函数可在地下介质中聚焦于任何期望点,并且充当产生或记录格林函数的虚拟点源或虚拟接收器.Broggini等(2012)和Wapenaar等(2013)将其发展到二维和三维.Meles等(2015)首次将Marchenko理论自聚焦技术用于压制层间多次波.其方法是通过虚源点与地表之间的直达波记录和地表反射响应得到上下行格林函数.然而,Marchenko理论自聚焦压制层间多次波技术仍然需要利用全局或局部匹配滤波器来从测量数据中减去预测的层间多次反射波(匡伟康等,2018),最小能量准则可能会损伤或破坏一次反射波.Meles等(2016)提出了一种基于Marchenko理论基准面重定和卷积干涉的近似一次反射获得方案,该方法仍不能避免估计下行聚焦函数的直接到达值.初始下行聚焦函数的估计需要首先建立一个平滑的速度模型,然后再应用该方法进行预测和减去.因此,现有的方法局限于自适应匹配滤波器以及层间多次波压制过程中的模型信息.该理论近些年被扩展到地表相关多次波(Singh et al., 2015, 2016;Slob and Wapenaar, 2017;Zhang and Slob, 2019)、弹性介质(Wapenaar,2014;Da Costa Filho et al., 2014;Urruticoechea and Wapenaar, 2017),以及其他许多方面.Van Der Neut和Wapenaar(2016)以及Zhang和Slob(2019)提出通过将特定深度上所有聚焦点的聚焦函数投影回采集表面上的接收器位置,消除了对初始下行聚焦函数的估计,并且不需要进行自适应减法,是一种完全数据驱动的层间多次波压制技术.本文依据该方法提出了无需先验信息和自适应相减的层间多次波压制策略,只利用地表获取的地震记录直接得到一次反射.带有高速薄透镜状地质体的数值模型测试证明了本方法的适用性.
1 理论方法依据单边格林函数表达式和Marchenko方程,在数据域中将聚焦波场从任意深度水平投影到地表,也就是说,来自聚焦深度的一次反射作为投影聚焦波场中最后一个事件出现.该方法仍是一种迭代法压制层间多次波的方案,使得利用地表获取的地震记录直接得到一次反射.
通常用t表示时间,z表示深度,(x, y)表示水平坐标.∂D0表示z0=0时的自由表面.zi表示地表以下任意深度,用∂Di表示该平面.将地表记录到的脉冲反射响应表示为R(x′0, x0, t),其中x0表示震源位置,x′0表示接收器位置,炮点、检波点均位于自由表面∂D0.格林函数G(xi, x0, t)是由在x0处激发的脉冲源定义的,并且接收器位于焦点xi处.格林函数是在与地表激发地表接收得到的单边反射响应相同的物理介质中定义的.Wapenaar等(2014)将z0<z<zi定义为截断介质.在截断介质中,介质的性质等价于物理介质的性质,在截断的介质外是无反射的理想介质.聚焦函数f1(x0, xi, t)是截断介质中波动方程的解,在t=0时聚焦于焦点xi处.聚焦函数是在截断介质和理想介质中定义的.聚焦函数和格林函数可以分为上行和下行两部分.
依据三维形式的通量归一化波场的单向互易定理,并将其用于深度z0和zi.当采集面z0以上的介质无反射时,格林函数可以表示为(Slob et al., 2014;Wapenaar et al., 2014):
(1) |
(2) |
其中上标“+”和“-”分别表示下行和上行波场.把下行聚焦函数和格林函数写成一个直接到达部分和一个尾波加和的形式:
(3) |
(4) |
其中f1d+和Gd+表示直接到达部分,f1m+和Gm+表示之后的尾波部分.初始下行聚焦函数可以理解为初始直达下行格林函数的逆函数,如Wapenaar等(2014)所述:
(5) |
式中,δ(xH)表示空间域的一个空间带限二维delta函数,δ(t)是一个时间域的时间delta函数.
根据式(5)可知,格林函数下行波场直接到达部分可以通过卷积积分将聚焦函数初始值折叠为一个delta函数,其中,Gd+和f1d+互为逆函数.根据Van Der Neut和Wapenaar(2016),将这个卷积积分算子应用到式(1)、(2)中,并根据反射层的反射是通过是格林函数上行波场中第一个事件与初始直达格林函数下行部分进行卷积得到为依据,可以得到:
(6) |
(7) |
其中,
(8) |
(9) |
式中t2表示从自由表面x′0点到聚焦平面深度zi,再回到自由表面x″0的最小双程走时.对于v-和vm+,式(8)、(9)可以只使用单边反射响应来求解,f1d+可以变为delta函数.由式(8)、(9)得到的v-和vm+可用来计算由式(6)得到的U-:
(10) |
将式(8)代入式(9),可以得到vm+表达式,对其进行简化,用R表示地表接收地震记录与波场的卷积积分算子,R*表示相关积分算子,最终得到vm+表达式:
(11) |
根据(Van Der Neut and Wapenaar, 2016)将式(11)展开为诺依曼级数,得到vm+表达式:
(12) |
将式(12)代入式(10),得到U-的最终表达式:
(13) |
如图 1以一个三层水平层状模型为例展示了U-的物理意义,图中以绿色五角星作为聚焦点,红色爆炸形标志代表震源.当聚焦点在实际的反射层上时,式(1)中的格林函数为聚焦点上的脉冲源在地表所产生的上行波场,如图 1a所示.将该反射定义为G-的第一个事件,并通过将格林函数与Gd+进行卷积,将∂Di处的所有焦点重新计算到地表一个点上,得到U-,如图 1b所示.并对所有时刻计算并记录U-,即得到只包含一次反射的地震记录.
该方法不需要任何模型信息也不必进行自适应减法.求得的数据集一次有效波保存的更加完整,层间多次波压制效果也更好.
2 模型测试 2.1 水平层状模型试算为验证本文方法的有效性,首先以一个简单的8层水平层状二维模型为例,进行了基于Marchenko理论一步法直接压制层间多次波方法的测试,模型如图 2所示,模型大小为6000 m×2000 m,图 2a为模型的速度信息,图 2b为模型的密度信息.
数值例子利用基于声波方程的高阶有限差分正演模拟,以主频25 Hz雷克子波激发,炮点位置与检波点位置重合,共501炮激发,每炮501道接收,炮间距和道间距均为5 m,采样点数1024、时间采样间隔4 ms,网格大小2.5 m×2.5 m,在模型顶部应用吸收边界条件.图 3a为地表激发地表接收的包含有效波和层间多次波的单炮地震记录,并且去除了直达波.
采用Marchenko理论一步法直接压制层间多次波,得到的单炮地震记录如图 3b所示.当将此结果与图 3a中层间多次波压制前的地表接收地震记录进行比较时,特别是在图 3箭头所指示的位置,可以很明显的看到,压制后的单炮地震记录在不损伤一次波的基础上层间多次波得到了有效衰减.
2.2 带高速透镜状地质体模型试算为进一步验证本文方法的适用性,接下来对夹薄高速透镜体的二维层状模型进行测试.速度模型如图 4所示,图 4a为速度模型,图 4b为密度模型,模型大小为6000 m×2000 m,共有8个反射边界.速度模型浅层深层均为水平层状,中间层中存在一个高速的薄透镜状地质体.
图 5a为利用有限差分法以主频25 Hz雷克子波正演得到的地表激发地表接收的单炮地震记录.采样记录时间4.092 s,采样间隔4 ms,网格大小2.5 m×2.5 m,道间距和炮间距均为5 m,总炮数为501,每炮接收道数501.其中,正演得到的原始地震记录中不包含自由表面多次波和直达波.
图 5b显示了对带高速透镜状地质体模型的2D合成数据应用基于Marchenko理论自聚焦方法得到的层间多次波预测相减后的单炮地震记录,采用基于Marchenko理论一步法直接压制层间多次波后的单炮地震记录如图 5c所示,在将两种方法的压制结果进行比较时,可以清楚地看到,基于Marchenko理论一步法压制层间多次波得到的效果更好,层间多次波被很好的消除,该方法可以在不影响一次反射波的前提下将层间多次波从地震记录中衰减干净,压制后的单炮地震记录更清晰,如图 5黄色箭头处所示.观察图 5b的压制结果,可以明显看出,尽管Marchenko理论自聚焦法也可压制层间多次波,但与图 5c所示结果相比,压制效果相对较差,当层间多次波和一次反射相互重叠时,匹配相减滤波器无法区分它们,因此它们都会被滤波器去除.如图 5所示,在应用预测相减法压制层间多次波后(图 5b),红色箭头所指示的有效波的能量受到损失.分别抽取层间多次波压制前和两种方法压制层间多次波后的零偏移距单道地震记录,如图 6所示,层间多次波压制前的地震记录道用浅蓝色实线表示,基于Marchenko理论一步法压制层间多次波后的地震道用橙黄色粗实线表示,基于Marchenko理论自聚焦方法压制层间多次波后的地震道用深蓝色虚线表示,从两种方法压制层间多次波前后单道地震记录对比图中可以更加直观的看出基于Marchenko理论一步法压制层间多次波的效果要优于自聚焦预测相减法压制层间多次波的效果.图 7所示的局部放大图可以更清晰的显示出黄色箭头指示的层间多次波压制区.图 7结果表明,在应用基于Marchenko理论的一步法压制层间多次波的结果中层间多次波被消除,而应用基于Marchenko理论自聚焦方法得到的层间多次波压制结果中存在明显的多次波能量残余.并且,图 7中红色箭头所指示的有效波,在应用Marchenko理论自聚焦压制层间多次波方法后能量有所衰减,而基于Marchenko理论一步法压制层间多次波使得一次反射在时间、振幅、相位等方面都得到了很好的保留和恢复,更接近原始模型中的地震记录道.层间多次波压制前后的共偏移距剖面如图 8a、b所示.图 8a中的剖面包含多处多次波噪声.然而,图 8b中得到的剖面几乎没有层间多次反射.
基于Marchenko理论自聚焦方法压制层间多次波要求必须提供一个下行聚焦函数的初始估计(先验信息)作为迭代的初值,这在实际地震资料处理过程中带来了诸多不便.为有效解决该问题,发展了完全数据驱动的基于Marchenko理论的一步法层间多次波压制技术,通过将Marchenko方程的两侧与初始格林函数卷积,将其投影到地表,可以只利用地表获取的地震记录,在双程走时域中迭代计算一次反射波,以消除所有阶次的层间多次波.
该技术将传统的虚源点投影到地表,无需虚源点的人工选择,也就是无需初始模型参数信息,另外也无需传统方法的多次波压制的匹配滤波处理步骤,可自动实现层间多次波压制.
二维数值模型试算表明,该方法在没有任何模型信息的情况下,可有效地去除地面地震资料中的层间多次波.和传统方法对比表明,基于Marchenko理论一步法直接压制层间多次波技术的计算精度更高,可在不损伤有效波的前提下有效消除层间多次波,具有很好的保持幅性.具有很好的发展前景,可望有效应用于实际地震资料的多次波压制处理.
Berkhout A J, Verschuur D J. 1997. Estimation of multiple scattering by iterative inversion, Part I: Theoretical considerations. Geophysics, 62(5): 1586-1595. DOI:10.1190/1.1444261 |
Berkhout A J, Verschuur D J. 2005. Removal of internal multiples with the common focus-point (CFP) approach. Part 1: Explanation of the theory. Geophysics, 70(3): V45-V60. |
Broggini F, Snieder R. 2012. Connection of scattering principles: A visual and mathematical tour. European Journal of Physics, 33(3): 593-613. DOI:10.1088/0143-0807/33/3/593 |
Broggini F, Snieder R, Wapenaar K. 2012. Focusing the wavefield inside an unknown 1D medium: Beyond seismic interferometry. Geophysics, 77(5): A25-A28. DOI:10.1190/geo2012-0060.1 |
Da Costa Filho C A, Ravasi M, Curtis A, et al. 2014. Elastodynamic Green's function retrieval through single-sided Marchenko inverse scattering. Physical Review E, 90(6): 063201. DOI:10.1103/PHYSREVE.90.063201 |
Guo M Q, Ba J, Max M, et al. 2019. Q compensation and reverse time migration imaging for structural complex area beneath shallow gas of marine gas field: a case study at Santos Basin, Brazil. Progress in Geophysics (in Chinese), 34(2): 717-723. DOI:10.6038/pg2019BB0421 |
Huang R, Chen X H, Li J Y. 2009. Simulation and effect analysis of interbed multiples. Progress in Geophysics (in Chinese), 24(3): 981-987. DOI:10.3969/j.issn.1004-2903.2009.03.022 |
Jakubowicz H. 1998. Wave equation prediction and removal of interbed multiples. //68th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
Jin D G, Chang X, Liu Y K. 2008. Algorithm improvement and strategy of internal multiples prediction based on inverse scattering series method. Chinese Journal of Geophysics (in Chinese), 51(4): 1209-1217. |
Kuang W K, Hu T Y, An S P, et al. 2018. Internal multiples attenuation using Marchenko autofocusing. Progress in Geophysics (in Chinese), 33(6): 2347-2357. DOI:10.6038/pg2018BB0474 |
Löer K, Curtis A, Meles G A. 2016. Relating source-receiver interferometry to an inverse-scattering series to derive a new method to estimate internal multiples. Geophysics, 81(3): Q27-Q40. DOI:10.1190/geo2015-0330.1 |
Meles G A, Löer K, Ravasl M, et al. 2015. Internal multiple prediction and removal using Marchenko autofocusing and seismic interferometry. Geophysics, 80(1): A7-A11. DOI:10.1190/geo2014-0408.1 |
Meles G A, Wapenaar K, Curtis A. 2016. Reconstructing the primary reflections in seismic data by Marchenko redatuming and convolutional interferometry. Geophysics, 81(2): Q15-Q26. DOI:10.1190/geo2015-0377.1 |
Shi Y, Wang W H, Li Y, et al. 2013. 3D surface-related multiple prediction approach investigation based on wave equation. Chinese Journal of Geophysics (in Chinese), 56(6): 2023-2032. DOI:10.6038/cjg20130623 |
Singh S, Snieder R, Behura J, et al. 2015. Marchenko imaging: Imaging with primaries, internal multiples, and free-surface multiples. Geophysics, 80(5): S165-S174. DOI:10.1190/geo2014-0494.1 |
Singh S, Snieder R, Van Der Neut J, et al. 2016. Accounting for free-surface multiples in Marchenko imaging. Geophysics, 82(1): R19-R30. |
Slob E, Wapenaar K, Broggini F, et al. 2014. Seismic reflector imaging using internal multiples with Marchenko-type equations. Geophysics, 79(2): S63-S76. DOI:10.1190/geo2013-0095.1 |
Slob E C, Wapenaar K. 2017. Theory for Marchenko imaging of marine seismic data with free surface multiple elimination. //87th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
Urruticoechea R, Wapenaar C P A. 2017. Elastodynamic single-sided homogeneous Green's function representation: theory and examples. //87th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts.
|
Van Der Neut J, Wapenaar K. 2016. Adaptive overburden elimination with the multidimensional Marchenko equation. Geophysics, 81(5): T265-T284. DOI:10.1190/geo2016-0024.1 |
Verschuur D J, Berkhout A J, Wapenaar C P A. 1992. Adaptive surface-related multiple elimination. Geophysics, 57(9): 1166-1177. DOI:10.1190/1.1443330 |
Wang W H, Cui B W, Liu H. 2007. Research progress in surface-related multiple attenuation. Progress in Geophysics (in Chinese), 22(1): 156-164. |
Wapenaar K. 2014. Single-sided Marchenko focusing of compressional and shear waves. Physical Review E, 90(6): 063202. DOI:10.1103/PHYSREVE.90.063202 |
Wapenaar K, Broggini F, Slob E, et al. 2013. Three-dimensional single-sided Marchenko inverse scattering, data-driven focusing, Green's function retrieval, and their mutual relations. Physical Review Letters, 110(8): 084301. DOI:10.1103/PhysRevLett.110.084301 |
Wapenaar K, Thorbecke J, Van Der Neut J, et al. 2014. Marchenko imaging. Geophysics, 79(3): WA39-WA57. DOI:10.1190/geo2013-0302.1 |
Weglein A B, Gasparotto F A, Carvalho P M, et al. 1997. An inverse scattering series method for attenuating multiples in seismic reflection data. Geophysics, 62(6): 1975-1989. DOI:10.1190/1.1444298 |
Wu J, Wu Z Q, Hu T Y, et al. 2013. Seismic internal multiple attenuation based on constructing virtual events. Chinese Journal of Geophysics (in Chinese), 56(3): 985-994. DOI:10.6038/cjg20130326 |
Zhang L L, Slob E. 2019. Free-surface and internal multiple elimination in one step without adaptive subtraction. Geophysics, 84(1): A7-A11. DOI:10.1190/geo2018-0548.1 |
郭梦秋, 巴晶, MaxMao, 等. 2019. 海上气田浅层气下覆复杂构造地区Q值补偿及逆时偏移成像——以巴西桑托斯盆地应用实例. 地球物理学进展, 34(2): 717-723. DOI:10.6038/pg2019BB0421 |
黄饶, 陈小宏, 李景叶. 2009. 层间多次波模拟与影响分析. 地球物理学进展, 24(3): 981-987. DOI:10.3969/j.issn.1004-2903.2009.03.022 |
金德刚, 常旭, 刘伊克. 2008. 逆散射级数法预测层间多次波的算法改进及其策略. 地球物理学报, 51(4): 1209-1217. DOI:10.3321/j.issn:0001-5733.2008.04.032 |
匡伟康, 胡天跃, 安圣培, 等. 2018. 基于Marchenko自聚焦的层间多次波压制. 地球物理学进展, 33(6): 2347-2357. DOI:10.6038/pg2018BB0474 |
石颖, 王维红, 李莹, 等. 2013. 基于波动方程三维表面多次波预测方法研究. 地球物理学报, 56(6): 2023-2032. DOI:10.6038/cjg20130623 |
王维红, 崔宝文, 刘洪. 2007. 表面多次波衰减的研究现状与进展. 地球物理学进展, 22(1): 156-164. DOI:10.3969/j.issn.1004-2903.2007.01.022 |
吴静, 吴志强, 胡天跃, 等. 2013. 基于构建虚同相轴压制地震层间多次波. 地球物理学报, 56(3): 985-994. DOI:10.6038/cjg20130326 |