频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

王冲， 顾汉明*， 许自强， 方中于， 唐永杰， 张健男， 胡发动

1 中国地质大学地球物理与空间信息学院， 地球内部多尺度成像湖北省重点实验室， 武汉 430074 2 中海油能源发展工程技术物探技术研究所， 广东湛江 524057



频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用

王冲1， 顾汉明1*， 许自强2， 方中于2， 唐永杰1， 张健男2， 胡发动2

1 中国地质大学地球物理与空间信息学院， 地球内部多尺度成像湖北省重点实验室， 武汉 430074 2 中海油能源发展工程技术物探技术研究所， 广东湛江 524057

海上倾斜缆采集技术具有多样的陷波特征，通过去鬼波处理可获得宽频数据.针对海水面波浪起伏及缆深误差引起的鬼波延迟时间估计误差以及崎岖海底和目的层深度变化使得鬼波和一次反射波的振幅差异系数随偏移距的变化而难以给定一个固定值的问题，本文推导出频率慢度域中鬼波滤波算子以及自适应迭代反演求解上行波算法，该鬼波滤波算子与不同水平慢度对应的鬼波和一次反射波的振幅差异系数以及鬼波延迟时间有关.并基于计算出的理论下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代最优计算该振幅差异系数和鬼波延迟时间.合成的及某海上采集的倾斜缆数据去鬼波处理结果表明，频率慢度域自适应迭代反演算法能较好地去除海上变深度缆鬼波，能达到拓宽地震记录频带目的.

倾斜缆采集； 频率慢度域； 自适应迭代反演； 鬼波压制； 宽频处理

1 引言

拓宽地震勘探信号的频带，不仅能够增强地震成像效果，提高分辨率，而且能够提高波阻抗反演和全波形反演的精度(ten Kroode et al.，2013；Soubaras et al.，2013).而海上拖缆地震采集的记录受海水面的虚反射(也称鬼波)影响，存在陷波特性，使得地震记录的频带变窄，不仅降低了地震剖面的分辨率，而且缺失了地震反演所需要的低频分量(Kragh et al.，2010；李套山等，1997).因此，近几年来，国内外开展了海水宽频带地震采集新方法试验(吴志强等，2013；吴志强，2014; 余本善和孙乃达，2015)，旨在充分利用地震采集数据的特点和有针对性的数据处理手段压制鬼波.这些试验包括海上双传感器地震拖缆采集(Carlson et al.，2007)、上下拖缆采集技术(Hill et al.，2006；刘春成等，2013；赵仁永等，2011)，以及变深度拖缆(VDS)或倾斜缆采集技术(Soubaras and Dowle，2010；Sablon et al.，2011)，其中倾斜缆采集技术因其充分利用同一条拖缆中陷波频率的多样性，通过在F-XY域、T-X域及τ-p域运用专门的宽频处理算法有效压制鬼波(Soubaras，2010；Soubaras，2012；Lin et al.，2011；Wang and Peng, 2012; Wang et al.， 2013；Song et al.，2015)，能较好地拓宽地震记录频带，因此，该技术在国内得到了研究和推广应用，并取得了较好的应用效果(谢玉洪等，2012；张振波和许自强，2014；许自强等，2015；王冲等，2016).现有的技术在推导鬼波滤波算子时仅假定海水面反射系数为-1，没有考虑到崎岖海底和目的层深度使得鬼波和一次反射波的振幅差异系数不再是一个恒定值而是随偏移距的变化而变化，进而影响了压制鬼波的效果，尤其在远偏移距时此问题更严重，为此，本文推导出频率慢度域中鬼波滤波算子以及最小二乘反演算法求解上行波算法，该鬼波滤波算子与不同水平慢度对应的鬼波延迟时间和海水面反射系数有关.并针对实际缆深误差引起的鬼波延迟时间估计误差以及鬼波与一次反射波振幅差异随偏移距和目的层深度的变化而难以给定一个固定值的问题，基于计算出的下行波与实际下行波之间的平方误差最小理论实现自适应反演迭代计算最优鬼波与一次反射波振幅差异系数和鬼波延迟时间.合成及某海上采集的倾斜缆数据处理结果表明，该方法能较好地压制变深度缆由海水面产生的鬼波，达到了拓宽地震记录频带的目的.

2 频率慢度域鬼波滤波算子

如图1所示，粗实线显示的倾斜拖缆在A点处的沉放深度为zi，偏移距为xi，i=1,…N，N为单炮记录总道数.不妨假设某一出射角为θj的上行平面波到达A点处的波前面为AB，该波前面的截距时间为τ，波场为uU(pj,τ)，pj为水平方向(沿着x方向)的慢度，且有

(1)

图1 斜缆处接收的上行波和下行波几何关系示意图Fig.1 The raypath geometry map of the up-going wave and down-going wave received at VDS

该平面波在到达海水面C点处，海水面反射下行达到A点处的下行波为uD(pj,τ)，在A点处接收到的总波场uS(pj,τ)是上行波场和下行波场之和，也即

uS(pj,τ)=uU(pj,τ)+uD(pj,τ).

(2)

由于下行波场uD(pj,τ)是上行波场在海水面反射向海底方向传播的下行波场，对于复杂构造介质尤其崎岖海底，地下同一反射界面上不同反射点反射的振幅是不同的，而在同一接收点上接收到的上行波是深层反射界面上某点处的一次反射波，同一接收点上接收的鬼波是深层反射界面上另外一点处的一次反射波上行至海水面并经海水面反射的下行波，而且这两个上行波在海底出射点也不同，导致在同一接收点上接收到的一次反射波和海水面反射的鬼波除了符号相反外，还存在振幅差异，该差异与构造复杂程度和目的层深度有关，也即鬼波和一次反射波的振幅差异随偏移距的变化而变化，因此，不妨假定鬼波与一次反射波振幅差异用系数Rj=R(pj)表示，其值大小与出射角有关，根据平面波传播理论以及图1中的几何关系，可以得到

uD(pj,τ)=RjuU(pj,τ-tBC-tAC),

(3)

(4)

令垂直方向(沿着z方向)的慢度为qj，则有

(5)

式(5)代入式(4)，则式(3)可以写成

uD(pj,τ)=R(pj)uU(pj,τ-τj),

(6)

这里，

(7)

式(6)代入式(2)得

uS(pj,τ)=uU(pj,τ)+RjuU(pj,τ-τj),

(8)

则对于M个水平慢度pj(j=1,…M)在接收点A处的总波场d(xi,t)可以表示为

(9)

基于时间域离散线性τ-p正反变换可知，式(9)中的uU(pj,τ)就是A点处的上行波场u(x,t)的τ -p正变换uU(p,τ)的第pj分量，且有τ=t-pjxi.

对式(9)做关于时间t的傅里叶变换，得

(10)

这里，G(xi,pj,ω)为鬼波滤波算子在频率慢度域表示形式，即

G(xi,pj,ω)=[1+Rje-jω τj]e-jωpjxi.

(11)

3 自适应最小二乘反演迭代去鬼波算法和步骤

对式(3)做关于时间t的傅里叶变换，得

UD(pj,ω)=RjU(pj,ω)e-jω τj,

(12)

对于倾斜缆N道而言，对每一个频率分量ω，式(10)可表示为矩阵形式：

D(x,ω)=G(x,p,ω)U(p,ω),

(13)

这里，D(x,ω)为N个元素组成列向量，即

D(x,ω)=[D(x1,ω),D(x2,ω),…,D(xN,ω)]T,

(14)

U(p,ω)为M个元素组成列向量，即

U(p,ω)=[U(p1,ω),U(p2,ω),…,U(pM,ω)]T,

(15)

G(x,p,ω)是由N×M个元素组成的矩阵，其元素可表示为

Gi,j(x,p,ω)=G(xi,pj,ω)=[1+Rje-jω τj]e-jωpjxi.

(16)

在利用式(7)计算延迟时间以及给定海水面反射系数前提下，就可采用最小二乘反演方法求解线性方程组(13)，可得(Yilmaz，1989；Song et al.，2015)

U(p,ω)= [G(x,p,ω)TG(x,p,ω)+εI]-1

×G(x,p,ω)TD(x,ω),

(17)

式中，ε为阻尼系数，I为单位矩阵.

由于拖缆沉放深度通常有误差，导致由式(7)计算的鬼波延迟时间τj有误差；其次，鬼波与一次反射波的振幅差异系数与平面波出射角有关，难以用一个恒定的值表示，也即式(16)难以准确计算，为此，本文采用自适应迭代计算最优振幅差异系数Rj和鬼波延迟时间τj，然后采用式(17)求解上行波，也即压制鬼波后的波场.

一方面，从式(12)可以看出，对某个慢度pj传播的上行波(一次反射波)施加一个算子Rje-jω τj便得到了该慢度传播的下行波(鬼波)；另一方面，由于倾斜缆接收到的总波场是上行波场和下行波场之和，因而如果上行波已知的话，就可以从下式计算下行波场，即

(18)

设第k次迭代时，对于某个水平慢度pj，可以基于式(17)计算出所有频率分量ω对应的上行波U(k)(pj,ω)，则下行波可通过下式求得，即

(19)

(20)

为了最小二乘求解式(20)，令

(21)

化简得

(22)

频率域τ -p变换自适应最小二乘迭代反演压制鬼波方法的具体步骤可概括如下：

(2) 对于所有的频率分量ω、所有的道和水平慢度pj，先假定振幅差异系数Rj为-1，并用式(7)计算鬼波延迟时间τj，利用式(11)计算得到G(xi,pj,ω)，进而可以得到矩阵G(x,p,ω)；

(3) 对于所有的频率分量ω和水平慢度pj，利用式(16)求解得到初始的上行波场向量U(p,ω)；

(5) 对于所有的频率分量ω和水平慢度pj，利用式(11)计算得到矩阵G(x,p,ω)，进而基于式(17)求解得到最终的上行波场向量U(p,ω)；

4 合成倾斜缆炮集记录试算

为了验证变深度缆最小二乘反演迭代去鬼波技术对由海水自由反射界面所引起的鬼波的压制效果，对一个复杂盐丘模型(图2a)的合成变深度缆记录进行去鬼波试算，正演模拟时，采用主频为30Hz的雷克子波，记录道数为300道，道间距为10m，采样间隔为2ms，记录长度为3.0s，偏移距为0m，单边放炮.斜缆采集是不同接收点处的拖缆沉放深度不同，斜缆的形态有多种，通常采用抛物线形态(AmundsenandLandrø, 2014)，其数学表达式如下：

图2 复杂盐丘速度模型(a)及拖缆沉放深度随偏移距变化(b)示意图Fig.2 Velocity model for complex salt (a) and the depth for streamer sinking versus offset (b)

(23)

这里可以看出，抛物线形态斜缆由拖缆沉放最小深度或起始深度zmin、最大沉放深度zmax以及拖缆倾斜段水平距离xc等三个参数控制，深度z是沿着Z轴测量的，炮检距x是沿着X轴测量的，第一式提供斜缆沉放深度数据，而第二式提供水平缆沉放深度.起始深度zmin通常为几米到十多米，xc也即拖缆倾斜段起始点到终止点之间的水平长度，其值从几百米到几千米.如图2b所示的拖缆三个参数分别为6m、50m和1500m.

图3a为基于有限差分法波动方程正演模拟的变深度缆单炮记录，炮点位置在400m，明显看出，在第一道时间在1.1s处有一组海底界面上的反射波同相轴，在1.5s、1.8s、2.1s、2.5s等处分别有一组地层界面上的反射波同相轴，可以看出每组反射波同相轴具有两个强相位特征，第一个强相位对应一次反射波，而续至的第二个强相位对应海水面的鬼波，而且随着炮检距的增大，一次波和鬼波逐渐明显分开又逐渐靠近，也可以从图4a和4c显示的局部放大结果看出这个特征.图3d为τ-p域炮集记录，也可以看出这些特征，随着慢度的增大，一次波和鬼波逐渐明显分开又逐渐靠近.

图6为原始和去鬼波后的t-x域单炮记录的频率波数谱(FK谱)，明显看出，原始炮集记录的FK谱(图6a)存在明显的陷波特征(箭头所示)，而去鬼波后的炮集记录的FK谱(图6b)上的陷波得到了消除，能量得到补偿，频谱连续，频带得到了拓宽，更进一步验证了自适应最小二乘反演迭代去鬼波算法可以有效压制大偏移距接收缆鬼波，弥补了由于鬼波存在引起的频谱缺失.

图3 去鬼波前后的合成单炮记录(a) 原始t-x域单炮记录； (b) 去鬼波后的t-x域单炮记录； (c) 残差或鬼波的t-x域单炮记录；(d) 原始τ -p域单炮记录； (e) 去鬼波后的τ -p域单炮记录； (f) 残差或鬼波的τ -p域单炮记录.Fig.3 Synthetic shot gather with a VDS acquisition for (a) before deghosted gather in t-x domain, (b) after deghosted gather in t-x domain, (c) residuals or ghost gather in t-x domain, (d) before deghosted gather in τ -p domain, (e) after deghosted gather in τ -p domain, (f) residuals or ghost gather in τ-p domain

图4 与图3对应的合成单炮原始记录t-x域局部放大显示(a) 原始单炮记录(近偏移距)； (b) 去鬼波后的单炮记录(近偏移距)；(c) 原始单炮记录(远偏移距)； (d) 去鬼波后的单炮记录(远偏移距).Fig.4 Synthetic shot gather enlarged display in Fig.3 for (a) before deghosted gather with near source, (b) after deghosted gather with near source, (c) before deghosted gather with far source, (d) after deghosted gather with far source

图5 去鬼波前后的合成单炮记录的频谱(a) 原始炮集记录各道振幅谱； (b) 去鬼波后的各道振幅谱； (c) 原始炮集记录各慢度对应的振幅谱； (d) 去鬼波后的各慢度对应的振幅谱.Fig.5 The amplitude spectra for (a) each trace in the gather before deghosted, (b) each trace in the gather after deghosted, (c) each slowness in the gather before deghosted, (d) each slowness in the gather after deghosted

图6 去鬼波前(a)后(b)的合成炮集记录FK谱Fig.6 F-K spectra for synthetic shot gather for (a) before deghosted and (b) after deghosted

图7 实际斜缆采集的单炮原始记录、去鬼波后的记录以及残差对比(a) 原始t-x域单炮记录； (b) 去鬼波后的t-x域单炮记录； (c) 残差或鬼波的t-x域单炮记录；(d) 原始τ -p域单炮记录； (e) 去鬼波后的τ -p域单炮记录； (f) 残差或鬼波的τ -p域单炮记录.Fig.7 Actual shot gather with a VDS acquisition for (a) before deghosted gather in t-x domain, (b) after deghosted gather in t-x domain, (c) residuals or ghost gather in t-x domain, (d) before deghosted gather in τ -p domain, (e) after deghosted gather in τ-p domain, (f) residuals or ghost gather in τ-p domain

图8 去鬼波前后的实际斜缆采集的单炮记录的频谱(a) 原始炮集记录各道振幅谱； (b) 去鬼波后的各道振幅谱； (c) 原始炮集记录各慢度对应的振幅谱； (d) 去鬼波后的各慢度对应的振幅谱.Fig.8 The amplitude spectra of actual shot gather with a VDS acquisition for (a) each trace in the gather before deghosted, (b) each trace in the gather after deghosted, (c) each slowness in the gather before deghosted,(d) each slowness in the gather after deghosted

图9 实际斜缆采集的去鬼波前后的炮集记录FK谱(a) 原始炮集波场FK谱； (b) 去鬼波后的波场FK谱； (c) 去鬼波前的海底反射FK谱； (d) 去鬼波后的海底反射FK谱.Fig.9 F-K spectra for actual shot gather with a VDS acquisition for (a) before deghosted and (b) after deghosted

图10 原始(a)与去鬼波(b)共偏移距剖面Fig.10 Common offset gather before deghosted (a), and after deghosted gather (b)

图11 去鬼波前后的共偏移距剖面频谱对比Fig.11 Amplitude spectra for common offset gather

5 实际倾斜缆采集数据试算

实际海上斜缆采集观测系统参数道间距为12.5m，炮间距为25m，接收长度6.25s，采样率为2ms，震源深度为5m，斜缆三个参数分别为5m、50m和3000m.由于海况影响，实际采集时电缆沉放深度与设计的理论深度之间存在一定的误差，实际采集时海水面也是波浪起伏，导致鬼波实际到达接收点的时间与理论时间存在一个差值，鬼波与一次波的振幅也不同.

图7a和图7d分别显示了涌浪干扰滤波后的t-x域、τ -p域倾斜缆单炮记录，明显可以看出每组反射波同相轴具有两个强相位特征，第一个强相位对应一次反射波，而续至的第二个强相位对应海水面的鬼波，而且随着炮检距的增大，一次波和鬼波逐渐明显分开又逐渐靠近.

基于自适应最小二乘反演迭代去鬼波算法，对原始单炮记录进行去鬼波处理，图8显示了去鬼波前后的不同域记录的振幅谱，与倾斜缆合成记录振幅谱特征类似，对比图中箭头位置处的振幅谱，可以看出去鬼波前的记录由于鬼波存在产生陷波特征，而在去鬼波后的记录频谱上，陷波能量得到较好的补偿，说明了自适应最小二乘反演迭代去鬼波算法可以较好地消除倾斜缆处对应的海水面的鬼波.

图7b、7e分别为实际斜缆采集的去鬼波后的t-x域和τ-p域倾斜缆单炮记录，图7c、7f分别为t-x域和τ-p域残差或所去除的鬼波数据，对比明显可以看出，去鬼波后的单炮记录上的反射波主要为一个正相位和一个负相位特征，远近偏移距道上的鬼波都能得到较好的压制，进一步说明了自适应最小二乘反演迭代去鬼波算法可以较好地消除倾斜缆处对应的海水面的鬼波.

图9为原始和去鬼波后的t-x域单炮记录的频率波数谱(FK谱)，明显看出，整炮原始记录的FK谱(图9a)存在明显的陷波特征(椭圆圈所示)，而去鬼波后的整炮记录的FK谱(图9b)上的陷波得到了消除，能量得到补偿，频谱连续，频带得到了拓宽.从图9c和图9d显示的炮集记录中对应海底反射的FK谱更明显看出这些特征，进一步验证了自适应最小二乘反演迭代去鬼波算法可以有效压制大偏移距接收缆鬼波，弥补了由于鬼波存在引起的频谱缺失.

图10显示出了偏移距为2800m的原始及去鬼波后的共偏移距剖面，可以明显看出，原始共偏移距剖面上反射波出现了两条同相轴，使得分辨率降低，尤其对深层反射同相轴影响更为严重；而去鬼波结果中压制了其中的一条同相轴，分辨率得到了提高.这也可以从图11显示出的整条地震记录剖面的一维平均叠加频谱分析结果看出，在原始记录频谱图(细实线)上有多个陷波频率，而在去鬼波压制后的频谱图(粗实线)上，陷波消去了，说明本文提出的去鬼波算法能充分利用变深度缆陷波频率发散特征，有效弥补陷波，扩宽频带.

还可以从压制鬼波前后的共中心点道集记录及其自相关谱(图12)分析去鬼波的效果，首先采用波场延拓法对变深度缆采集的原始和去鬼波炮集数据进行校正到水平缆，然后进行抽道集和动校正，得到动校正后的共中心点道集记录.对比图12a和图12b，从动校正后的原始共中心点道集记录，可以明显看出存在于一次反射波后的鬼波，且存在明显的剩余时差，而动校正后的去鬼波后的道集记录上，只有一次反射波同相轴，鬼波得到较好的压制.图12c和图12d分别对应去鬼波前后的CMP道集上黑色粗框窗口内各道的自相关图.对比可以看出，在图12c上，由于鬼波的存在，使得每一条自相关曲线上都呈现出多个极值，而在图12d上，由于鬼波得到压制，使得每一条自相关曲线上只有一个极值，这些图件进一步说明了频率慢度域自适应最小二乘反演迭代去鬼波算法在压制变深度缆鬼波的应用效果.

图13显示了去震源和接收缆鬼波前后的偏移叠加剖面，对比可以看出，原始偏移叠加剖面上海底及目的层由正负峰组成的多条同相轴在去鬼波后的最终处理剖面上变成了单峰组成的同相轴，海水面引起的鬼波得到明显压制，地震记录的分辨率得到明显提高，达到了宽频处理的目的.

6 结论

通过模型变深度缆模拟记录及实际采集的变深度缆去鬼波处理结果分析，可以得出如下结论：

(1) 频率慢度域导出的鬼波算子中同一水平慢度对应不同深度反射层的鬼波延迟时间以及鬼波与一次波的振幅差异系数相同，克服了在时空域中鬼波算子中在同一炮检距下难以给定一个固定的鬼波延迟时间以及鬼波振幅差异系数的问题.

(2) 在频率慢度域可以通过理论下行波与实际下行波之间的平方误差最小实现自适应反演迭代计算得到鬼波最优延迟时间以及振幅差异系数，进而通过频率慢度域阻尼最小二乘反演算法计算得到压制鬼波后的上行波，合成记录验证了该算法的可靠性.

(3) 将频率慢度域自适应最小二乘反演迭代算法应用于某海上变深度缆去鬼波技术后的炮集记录，其振幅谱频带明显比原始记录的振幅谱频带宽，在最终处理的偏移叠加剖面上，中高频谱振幅能量得到提高，使得频带变得更宽，中深层反射的分辨率得到明显提高，信噪比也得到提高，断点更清晰.

Carlson D, Long A, Söllner W, et al. 2007. Increased resolution and penetration from a towed dual-sensor streamer.FirstBreak, 25(12): 71-77.

Hill D, Combee L, Bacon J. 2006. Over/Under acquisition and data processing: the next quantum leap in seismic technology?.FirstBreak, 24(6): 81-95.

Kragh E, Muyzert E, Curtis T, et al. 2010. Efficient broadband marine acquisition and processing for improved resolution and deep imaging.TheLeadingEdge, 29(4): 464-469.

Amundsen L, Landrø M. 2014. Broadband seismic technology and beyond, Part VII: CGG′s BroadSeis—a change of thinking.GeoExPro, 11(1): 50-52,54.

Li T S, Gao H Q, Liu Z X, et al. 1997. The acquisition of ghost information and the study of its generation and frequency response.GeophysicalProspectingforPetroleum(in Chinese), 36(4):38-44.Lin D, Sablon R, Gao Y, et al. 2011. Optimizing the processing flow for variable-depth streamer data.FirstBreak, 29(9): 89-95. Liu C C, Liu Z B, Gu H M. 2013. The determination of optimal sinking depths of over/under streamers in offshore survey by merge operator.GeophysicalProspectingforPetroleum(in Chinese), 52(6): 623-629.

Sablon R, Russier D, Zurita O, et al. 2011. Multiple attenuation for variable-depth streamer data: From deep to shallow water. ∥81st Annual International Meeting, Society of Exploration Geophysicists.ExpandedAbstracts, 3505-3509.

Song J G, Gong Y L, Li S. 2015. High-resolution frequency-domain Radon transform and variable-depth streamer data deghosting.AppliedGeophysics, 12(4): 564-572.

Soubaras R, Dowle R. 2010. Variable-depth streamer—a broadband marine solution.FirstBreak, 28(12): 89-96.

Soubaras R. 2010. Deghosting by joint deconvolution of a migration and a mirror migration. ∥80th Annual International Meeting, Society of Exploration Geophysicists. Expanded Abstracts, 3406-3410.

Soubaras R. 2012. Pre-stack deghosting for variable-depth streamer data. ∥82nd Annual International Meeting, Society of Exploration Geophysicists. Expanded Abstracts, 1-5. Soubaras R and Lafet Y. 2013. Variable-depth streamer acquisition: Broadband data for imaging and inversion.Geophysics,78(2), WA27-WA39.

ten Kroode F, Bergler S, Corsten C, et al. 2013. Broadband seismic data—the importance of low frequencies.Geophysics, 78(2): WA3-WA14.

Wang C, Gu H M, Xu Z Q, et al. 2016. The application of least-squares inversion iteration algorithm to deghost for marine variable-depth streamer data.ChineseJ.Geophys. (in Chinese), 59(5): 1790-1803, doi: 10.6038/cjg20160522.Wang P, Peng C. 2012. Premigration deghosting for marine streamer data using a bootstrap approach. ∥82nd Annual International Meeting, Society of Exploration Geophysicists. Expanded Abstracts, 1-5.

Wang P, Ray S, Peng C, et al. 2013. Premigration deghosting for marine streamer data using a bootstrap approach in tau-p domain. ∥83rd Annual International Meeting, Society of Exploration Geophysicists. Expanded Abstracts, 4221-4225.

Wu Z Q, Yan G J, Tong S Y, et al. 2013. New advances in marine

图12 去鬼波前(a)、后(b)的共中心点道集动校正后的记录及其自相关谱(c,d)(a) 原始共中心点道集动校正后的记录； (b) 去鬼波后的共中心点道集动校正后的记录； (c) 对应图(a)框图内的记录自相关谱； (d) 对应图(b)框图内的记录自相关谱.Fig.12 NMO corrected CMP gather for (a) before deghosted, (b) after deghosted, autocorrelation spectra for (c) before deghostred and (d) after deghosted

seismic data acquisition technology and enlightenment to domestic offshore petroleum seismic exploration.ProgressinGeophys. (in Chinese), 28(6): 3056-3065, doi: 10.6038/pg20130629. Wu Z Q. 2014. New advances in marine broadband seismic exploration.OilGeophysicalProspecting(in Chinese), 49(3): 421-430.Xie Yuhong, Li Lie, Yuan Quanshe. 2012. Broadband marine seismic exploration in Qiongdongnan Basin deepwater areas.OilGeophysicalProspecting(in Chinese), 4(3): 431-435.

Xu Z Q, Fang Z Y, Gu H M, et al. 2015. The application of optimal deghosting algorithm on marine variable-depth streamer data.GeophysicalProspectingforPetroleum(in Chinese), 54(4): 404-413.

Yu B S, Sun N D. 2015. Latest development of marine braodband seismic acquisition technology.OilForum(in Chinese), 34(1): 41-45.Yilmaz O. 1989. Velocity-stack processing.GeophysicalProspecting, 37(4):357-382.Zhang Z B, Xu Z Q. 2014. Discussion on seismic data processing techniques series in deepwater area: A case study of Liwan 3-1 gas field.ProgressinGeophys. (in Chinese), 29(5): 2320-2325, doi: 10.6038/pg20140549.

Zhao R Y, Zhang Z B, Xuan Y H. 2011. Application of over/under streamer and over/under source seismic acquisition in the Pearl River Mouth Basin.OilGeophysicalProspecting, 46(4): 517-521.

附中文参考文献

李套山,高洪强,刘振夏等. 虚反射信息的采集及其形成机制、频率响应的理论探讨. 石油物探,1997,36(4):38-44+27.

刘春成, 刘志斌, 顾汉明. 2013. 利用上/下缆合并算子确定海上上/下缆采集的最优沉放深度组合. 石油物探, 52(6): 623-629.

王冲, 顾汉明, 许自强等. 2016. 最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用. 地球物理学报, 59(5): 1790-1803, doi: 10.6038/cjg20160522.

吴志强, 闫桂京, 童思友等. 2013. 海洋地震采集技术新进展及对我国海洋油气地震勘探的启示. 地球物理学进展, 28(6): 3056-3065, doi: 10.6038/pg20130629.

吴志强. 2014. 海洋宽频带地震勘探技术新进展. 石油地球物理勘探, 49(3): 421-430.

谢玉洪, 李列, 袁全社. 2012. 海上宽频地震勘探技术在琼东南盆地深水区的应用. 石油地球物理勘探, 47(3): 431-435.

许自强, 方中于, 顾汉明等. 2015. 海上变深度缆数据最优化压制鬼波方法及其应用. 石油物探, 54(4): 404-413.

余本善, 孙乃达. 2015. 海上宽频地震采集技术新进展. 石油科技论坛, 34(1): 41-45.

图13 去鬼波前后的偏移叠加剖面对比(a) 原始偏移叠加剖面； (b) 去鬼波后的最终处理的偏移叠加剖面.Fig.13 Migration stack for (a) before deghosted and (b) after deghosted and final processed

张振波, 许自强. 2014. 深水区地震资料处理技术系列探讨——以荔湾3-1气田为例. 地球物理学进展, 29(5): 2320-2325, doi: 10.6038/pg20140549.

赵仁永, 张振波, 轩义华. 2011. 上下源、上下缆地震采集技术在珠江口的应用. 石油地球物理勘探, 46(4): 517-521.

(本文编辑 胡素芳)

The application of adaptive iteration inversion algorithm to deghost for marine variable-depth streamer data in frequency-slowness domain

WANG Chong1, GU Han-Ming1*， XU Zi-Qiang2, FANG Zhong-Yu2, TANG Yong-Jie1, ZHANG Jian-Nan2, HU Fa-Dong2

1HubeiSubsurfaceMulti-scaleImagingKeyLaboratory,InstituteofGeophysicsandGeomatics，ChinaUniversityofGeosciences，Wuhan430074，China2CNOOCEnergyTechnology&Services-Development&ProspectingGeophysicalCo.,GuangdongZhanjiang524057,China

The variable-depth streamer (VDS) acquisition technology in marine has receiver notch diversity. The broadband data can be obtained by special deghosting method. The choppy sea sarface and error of streamer depth will cause the error of the ghost-delay time with respect to primary reflection. The amplitude difference coefficient between ghost and primary reflection varies with offset due to rugged seabed and target depth variation, so the coefficient is difficult to be given a fixed value. The ghost filtering operator and adaptive iteration inversion algorithm to solve up-going wave in frequency-slowness domain are deduced. The ghost filtering operator is related to the ghost-delay time and the amplitude difference coefficient between ghost and primary reflection at each horizontal slowness. Then, the amplitude difference coefficient and ghost-delay time are optimally inverted by the adaptive iteration inversion algorithm through a least-squares fit of the theoretical down-going wave to actual down-going wave. The results of deghosting for synthetic and marine VDS seismic data show that the method can availably deghost for marine VDS data and achieve the purpose of broadband seismic data.

Variable-depth streamer acquisition； Frequency-slowness domain； Adaptive iteration inversion； Deghosting； Broadband processing

10.6038/cjg20161227.

国家重大科技专项子专题“南海深水区油气勘探地球物理关键技术”(2016ZX05026-001)资助.

王冲，1992年生，中国地质大学(武汉)博士在读，主要从事地震勘探方法研究.E-mail：wangchongluck@163.com

*通讯作者 顾汉明，1963年生，男，教授，博士生导师，主要从事油气地震勘探与开发研究.E-mail：hmgu@cug.edu.cn

10.6038/cjg20161227

P631

2016-03-08，2016-07-19收修定稿

王冲， 顾汉明， 许自强等. 2016. 频率慢度域自适应迭代反演算法压制海上倾斜缆鬼波方法及其应用. 地球物理学报，59(12):4677-4689,

Wang C, Gu H M， Xu Z Q, et al. 2016. The application of adaptive iteration inversion algorithm to deghost for marine variable-depth streamer data in frequency-slowness domain.ChineseJ.Geophys. (in Chinese),59(12):4677-4689,doi:10.6038/cjg20161227.