南黄海崂山隆起区水深较浅，新近系和古近系界面为一强阻抗界面(T₂)，多次波类型较多，主要为与海底相关的短周期多次波、海底相关的长周期多次波和层间多次波(图 2)。选取穿越崂山隆起南北向地震测线A进行分析(图 1)，从最小偏移距的共炮检距剖面上看(图 3)，多次波干扰比较严重，在新生界以下无法识别有效反射信号，海底相关多次波与强阻抗界面(新近系和古近系界面)产生的多次波混叠在一起，多次波类型以短周期为主。

图  2  南黄海崂山隆起多次波示意图

Figure 2.  Sketch showing various multiples on the Laoshan uplift of south Yellow Sea

图  3  最小偏移距的共炮检距剖面

Figure 3.  Common-offset profile of minimum offset

南黄海崂山隆起区水深较浅，同时存在强阻抗界面，使得地震资料多次波类型丰富，为了更好地压制多次波，首先通过自相关拾取精确的海底反射时间(海底多次波周期)，然后采用DWD+SRME组合方法，压制海底相关的多次波；对于层间多次波，采用高精度Radon多次波衰减方法来进行压制。

该地震测线水深范围为26~56 m，由于水深较浅，海底反射波与直达波旅行时差较小，不能完全区分开，在共炮检距剖面上(图 4)只有左侧能够识别直达波和海底反射波，但在右侧直达波与海底反射波重叠在一起，无法拾取准确的海底反射时间，不能反映海底多次波的真实周期。

图  4  共炮检距剖面

Figure 4.  Common-offset profile

众所周知自相关谱上可以直观地显示多次波的周期变化，因此，在共炮检距剖面拾取海底反射时间的基础上，选取合适的时窗，进行自相关分析，在自相关谱上，在已拾取的海底反射时间约束下，提取地震最大振幅属性来自动拾取海底多次波周期。由图 5可以看出，仅根据共炮检距剖面上拾取的时间与真实的海底多次波周期还存在一定差异，而通过自相关谱上最大振幅属性提取获得的时间与自相关谱上多次波同相轴非常吻合，表明提取的海底多次波周期非常准确，为下一步海底多次波的压制奠定了基础。

图  5  自相关谱上自动提取海底多次波周期

Figure 5.  Submarine multiple cycles automatically extracted from autocorrelation spectra

DWD(Deterministic water-layer demultiple)是浅水区压制多次波的有效方法，该方法与传统反褶积压制多次波的方法相比，避免了反褶积对资料造成的“伤害”；DWD与SRME方法相比，对近道数据的要求更低，更易于实现，并且更适合浅水海洋环境地震资料的处理^[8-10]。DWD方法首先通过自相关分析获得海底多次波周期时间，然后变换到Taup域计算海底多次波周期，将数据静校正到一阶多次波，乘以测算的反射系数，采用自适应相减，再经过Taup反变换，得到海底多次波模型，最后通过自适应相减，来压制海底相关的短周期多次波(图 6)。

图  6  DWD方法工作原理

Figure 6.  The principle of DWD method

DWD方法对海底相关的短周期多次波可以较好地压制，但还存在自由表面相关的长周期多次波，该类多次波可以用SRME方法来进行压制，在使用SRME方法前，需要对海底反射进行切除^[11-13]，避免SRME对短周期多次波再次预测。

SRME方法基于波动方程的多次波衰减方法^[14]。在均匀介质模型下，把地震波在地下介质的脉冲响应定义为x₀(t)，它包含了地下所有有效波和层间多次波，那么一阶表面多次波可表示为：

同理，二阶表面多次波表示为：

因此，所有表面相关多次波可以由整个响应与一次反射褶积得到，这种通过地震数据自身褶积，一次反射波变成多次波，低级多次波变成高级多次波，这样就可以同时预测出自海底表面向下反射的一次或多次的所有与表面有关的多次波。SRME方法适用于压制长周期的海底相关多次波，最大的优势是预测过程不需要知道地下介质信息^{[14, 15]}，缺点是要求待处理的地震数据中包含从零炮检距到最大炮检距的全波场地震数据^[14]。

Radon变换是基于一次波与多次波速度时差的滤波方法^{[16, 17]}。对CMP道集进行动校后，一次波同相轴被拉平，多次波同相轴呈向下弯曲状，把动校后的CMP道集进行Radon正变换后, 一次波能量分布在零P道附近，而多次波能量呈弯曲条带状分布在远离零P道的地方^[16]，这样就可以在τ-p域中把多次波从地震资料中分离出来。该方法分辨能力强、振幅和相位的保真性好，其缺点是可能损伤浅层的高频成分，且对近道的多次波压制相对不足^[14]。

抛物线Radon正变换公式：

反变换公式：

式中：u(t, x)是在时空域数据，v(τ, q)是抛物线Radon域数据，q是以零偏移距为中心的抛物线的曲率；τ是抛物线在时间轴的截距。

南黄海崂山隆起区水深较浅，海底相关的多次波发育，首先采用确定性水层多次波压制方法(DWD)，对海底相关的短周期多次波进行压制，通过DWD方法压制多次波后，与海底相关的短周期多次波得到了较好的压制(图 7(b))，特别是强阻抗界面(新近系和古近系界面)的海底相关短周期多次波压制效果明显，从图 7(b)蓝色虚线框中可以看出，海底相关多阶多次波得到了衰减，这种海底相关的短周期多次波压制，改善了地震剖面的波组特征。

图  7  DWD压制多次波效果分析

Figure 7.  Effectiveness analysis of multiple waves suppressed by DWD method

在DWD方法去除短周期海底相关多次波的基础上，采用SRME技术对自由表面长周期的多次波进行压制，图 8(a)为DWD+SRME组合去除多次波后的叠加剖面，从剖面上可以看出海底相关的多次波得到了很好的压制，波组特征较好，与常规的SRME方法相比(图 8(b))，剖面上1~1.5 s处，DWD+SRME组合去多次波法对强阻抗界面产生的多次波压制较好；而采用Taup域预测反褶积去除多次波(图 8(c))，剖面上深部还存在比较明显的多次波。

图  8  DWD+SRME组合去多次波对比分析

Figure 8.  Comparison analysis after multiples wave depression with combined DWD and SRME methods

分别对以上3种去多次波后的叠加剖面做自相关分析，从自相关谱上(图 9)来看，DWD+SRME组合法压制海底相关多次波取得了良好的效果，图 9(b)中主要能量集中于零时间处，两侧的同相轴较少；图 9(c)中，零时间的两侧还存在短周期的同相轴，说明SRME方法没有完全去除短周期的海底相关多次波；Taup域预测反褶积方法后，自相关谱上存在间隔时间较大的同相轴，表明长周期多次波衰减不彻底。

图  9  DWD+SRME组合去多次波后自相关谱对比分析

Figure 9.  Comparison analysis of Spectrum autocorrelation after DWD and SRME methods to suppress multiples wave

通过对DWD+SRME组合法、SRME方法和Taup域预测反褶积去除多次波的效果对比，表明在浅水区，采用DWD+SRME组合法压制海底相关多次波效果较好，DWD方法主要压制海底相关的短周期多次波，SRME方法用来衰减自由表面相关的长周期多次波。

在压制海底相关多次波的基础上，对于剩余的层间多次波采用Radon变换来进行衰减，由于地震资料信噪比较低，在速度谱上(图 10(a))，1.2s以下无法识别有效的能量团，存在严重的多次波干扰，有效信息基本被多次波湮没，很难从速度谱上获取可靠的速度，经过DWD+SRME组合去除多次波后(图 10(b))，速度谱上多次波能量得到一定衰减，1s以下的有效波能量团逐步显现出来，采用Radon变换去除多次后(图 10(c))，有效波的能量信息得到很大的提升，多次波信息基本被压制干净。

图  10  Radon变换去除多次波速度谱分析

Figure 10.  Velocity spectrum analysis after multiples wave depression with Radon transform

图 11为Radon变换衰减多次波后，经过叠前时间偏移的叠加剖面，剖面上多次波基本被压制，在2.5s附近有效反射同相轴得到了体现，可以看出在南黄海低信噪比资料中，多次波压制非常关键，采用DWD+SRME组合压制多次波和Radon变换衰减多次波方法，在该研究区是比较适用的。

图  11  Radon变换去多次波后叠前时间偏移剖面

Figure 11.  Prestack time migration profile after multiples wave depression with Radon tranform

南黄海盆地崂山隆起地震资料深部中—古生界地震成像较弱，波组不连续—较连续，原因主要是经历了多期构造运动改造。在晚印支—早燕山时期华南板块继续向华北板块之下俯冲，陆-陆碰撞造成秦岭-大别-胶南造山带隆升，产生前陆盆地带，造山带南侧发生了扬子板块北缘前陆盆地带^[18]，该期的构造运动使得崂山隆起基本定型，同时断裂的逆冲活动加剧，造成部分海相中—古生界的缺失^[19]；晚燕山至早喜马拉雅期扬子地区为张性构造环境，由于太平洋板块对欧亚板块俯冲产生的弧后拉张效应与燕山期挤压冲断后的应力松弛有关^[20]，崂山隆起的北缘早期的逆冲断裂呈现负反转，形成早白垩世—古近纪的断陷；渐新世—中新世时期，印度板块和欧亚板块最终碰撞拼接，以及菲律宾板块的吕宋弧与欧亚板块斜向碰撞，南方遭受强烈挤压^[21]，古近系的回返，伸展构造被掀斜或褶皱、剥蚀，形成古近系和新近系之间的角度不整合，断裂活动减弱，崂山隆起形成新近纪以来的区域性覆盖沉积。

崂山隆起为刚性的结晶基底，为稳定地块^[22]，基底之上具有上、中、下三元结构^[23](图 12，地质剖面测线位置图见图 1)，上部为第四系—新近系构造层；中部为古近系—泥盆系构造层，位于下志留统高家边组滑脱面之上，构造相对复杂，以滑脱、冲断结构为主，褶皱及逆冲叠瓦断裂发育；下部为志留系—震旦系构造层，位于下志留高家边组滑脱面之下，内幕反射层比较简单，断裂不发育，地层较平缓，构造变形相对较弱^[24]。崂山隆起在加里东和海西期造成了区内幅度不大的差异沉降，其中在加里东期开始初步隆升，海西期出现凹凸雏形，局部还形成了逆冲作用不强的逆断裂，但总体上这两期的构造运动对崂山隆起的改造作用不强。

崂山隆起总体是一个比较平缓的箱状背斜构造，由于逆冲推覆作用，隆起上发育多个逆冲背斜构造(图 12)，EW向断裂是崂山隆起上的主体断裂，断裂规模较大，控制了崂山隆起的边界和主要构造轮廓，早期NE向和NW向断裂为一对共轭断裂，后期被卷入盖层滑脱。隆起上海相中、古生代地层发育齐全，且厚度较大，古生代发育一套陆棚—台缘斜坡相沉积，在印支—早燕山时期大规模隆升、剥蚀，局部缺失上古生界、中生界大部分地层。

图  12  南黄海崂山隆起地震解释剖面

Figure 12.  Seismic interpretation profiles of Laoshan Uplift of the South Yellow Sea Basin

(1) 南黄海地震资料多次波类型丰富，除了海底相关的多次波外，还存在强阻抗界面(新近系和古近系界面)产生的层间多次波，另外信噪比较低，深部有效的反射信号弱，给多次波压制、目的层成像工作带来了极大困难。

(2) 采用DWD+SRME组合法应用于南黄海地震资料，通过与SRME方法、Taup域预测反褶积方法去除多次波的效果对比分析，表明DWD+SRME组合法适用于南黄海地震资料，该方法有效地去除了海底相关多次波，并且没有损害有效的一次反射信号。采用Radon变换较好地衰减了层间多次波，在速度谱上，有效信号能量团得到了加强，采用叠前时间偏移技术，使深部地层有效反射信号得到较好成像。

(3) 通过地震资料构造分析崂山隆起是一个比较平缓的箱状背斜构造，背斜上发育多个逆冲断裂，在加里东期初具雏形，印支期隆起形成，整体抬升，形成比较平缓的格局；晚印支—早燕山隆起大幅度抬升，格局基本定型，并形成滑脱构造；晚燕山—喜山早期，北缘早期逆冲断裂负反转活动，形成晚白垩—古近系断裂；新近纪开始断裂活动减弱，崂山隆起形成新近纪以来的区域性覆盖沉积。