2017, Vol. 32
2. 中国地质科学院, 北京 100037
3. 国土资源部海洋油气资源和环境地质重点实验室, 青岛海洋地质研究所, 青岛 266071
4. 海洋国家实验室海洋矿产资源评价与探测技术功能实验室, 青岛 266071
5. 斯伦贝谢科技服务 (北京) 有限公司, 北京 100015
2. Graduate School of Chinese Academy of Geological Sciences, Beijing 100037, China
3. Key Laboratory of Marine Hydrocarbon Resource and Environmental Geology, Ministry of Land and Resources, Qingdao Institute of Marine Geology, CGS, Qingdao 266071, China
4. Laboratory for Marine Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China
5. Schlumberger Technology Services (Beijing) Ltd, Beijing 100015, China
大陆架是陆地向海洋的自然延伸，由于其丰富的油气资源储量，已成为海洋油气资源勘探的重要区域.我国海域拥有宽阔的大陆架，近海油气资源丰富.通过多年的勘探开发，已在渤海、东海陆架、珠江口、北部湾、莺歌海等沉积盆地相继实现了油气产量 (姜亮，2003；刘光鼎和陈洁，2005；范春花等，2013).由于近海大陆架区水浅，海面和海底这两个强反射系数界面的存在使得自由表面多次波异常发育.浅水自由表面多次波能量强、周期短、阶次高，严重降低地震资料信噪比，影响对有效波的正确识别.若不进行有效压制，会给后期的速度分析、偏移成像等环节带来困难，最终影响地震资料的成像精度.
目前，实际处理中多用基于反馈迭代思想 (肖二莲等，2010；范宝仓和谢涛，2012；金明霞等，2014) 的表面多次波压制方法 (Surface-Related Multiple Elimination，SRME) 或是利用多次波周期性的线性拉东域 (Tau-p域) 反褶积方法 (赵昌垒等，2013) 来进行浅水区自由表面多次波压制.其中，SRME方法数据驱动、不依赖于速度和层位等先验信息的特点使其迅速推广，在海洋深水地震资料去多次中普遍应用.但是该方法应用于浅水资料存在近偏移距海底反射估计困难、多次迭代易引入振幅误差等缺陷，导致其在浅水区应用效果不佳.Tau-p域反褶积方法利用多次波在线性拉东域更好的周期特性对其进行压制，但应用于浅水资料时有损伤有效波的风险 (刘俊等，2016；张治忠等，2016).
本文首先阐述了自由表面多次波的定义并对其进行了分类，对常用的预测反褶积及SRME方法进行了数学推导并指出了这两种方法应用于浅水地震资料时存在的主要技术缺陷.针对这些方法的局限性，使用确定性水层多次波压制方法 (Deterministic Water-layer Demultiple，DWD) 来压制水层 (海面和海底间) 多次波，对该方法的实现过程和参数的选取进行了讨论和说明.该方法利用水体反射周期来预测多次波，能够避免浅水近偏移距数据缺失对SRME技术的影响.同时能够去除与一次波速度差异较小的多次波且不伤害有效波.应用此方法后将海底切除，利用SRME技术来压制残留长周期自由表面多次波.将DWD方法与SRME方法结合，组合压制浅水区的自由表面多次波.来自东海陆架盆地的实际地震资料处理结果表明，该技术组合能够较好地解决浅水区自由表面多次波压制问题.结合高精度拉东变换去除其余长周期多次波后，地震资料的成像质量得以提升.
1 自由表面多次波压制 1.1 自由表面多次波定义及分类自由表面多次波指地下介质反射的地震波到达自由表面 (海洋中的为海水面) 后，至少在自由表面发生一次下行反射，然后在地下经一定的传播路径后重新返回自由表面所形成的地震波 (Dragoset and Jeričević, 1998).也就是说，只要地震波传播过程中在自由表面存在下行反射点即可定义为自由表面多次波.图 2为海洋地震勘探常见的多次波传播路径.根据上述定义，除了 (g) 的下行反射发生在海底与地层界面间，被定义为层间多次波外，其余多次波均至少在自由界面发生一次下行反射，可定义为自由表面多次波.此图同样可见，自由表面多次波类型丰富，特别是在浅水数据中，主要的多次波干扰是由自由表面所产生的 (黄新武等, 2005).如图 3实际浅水地震数据的近道剖面所见，有效反射湮没在自由表面多次波中无法识别.如果能够较好地压制自由表面多次波，也就压制了大部分的浅水多次波能量.因此浅水地震资料中自由表面多次波的压制是多次波压制乃至地震成像的重要环节.
|
图 2 多次波传播路径图 Figure 2 Ray paths of some multiple events |
|
图 3 浅水地震数据近道剖面 Figure 3 Near offset profile of seismic data in shallow water area |
对自由表面多次波进行分类，将地震波传播路径中最后一次下行反射发生在自由界面，再从海底向上反射直接被检波器接收的自由表面多次波称为水层多次波，传播路径如图 2a-d.将传播路径如图 2e-f类型的自由表面多次波称为非水层多次波.
1.2 主要压制方法自由表面多次波压制方法分为滤波法和波动理论方法两大类 (Weglein, 1999; 石颖和邢小林, 2011; 宋家文等, 2014)，如表 1.
|
表 1 自由表面多次波压制方法 Table 1 suppression methods of surface-related multiple |
在滤波类方法中，F-K滤波、抛物和双曲Radon变换滤波利用一次波和多次波速度差异，将动校正后的道集变换域后设计滤波函数来进行多次波的分离和制.CMP叠加利用动校正后一次波被拉平，多次波仍有剩余时差无法同相叠加的特性来压制多次波.预测反褶积利用多次波的周期性，设计预测滤波器来消除地震数据中的周期成分.由于时空域远偏移距数据中的多次波周期性不佳，又将该方法拓展到了Tau-p域，利用水平层状介质中传播角度相同的多次波严格的周期性特征来进行压制.Tau-p域预测反褶积是目前较为常用的滤波类多次波压制方法，其数学实现过程如下：
首先将时空域的道集数据变换到Tau-p域，其离散形式为
|
(1) |
其中τ, p为截距和慢度，n为变换所参与的道数.f(x, t) 为时空域道集，F(τ, p) 为其Tau-p变换，且有
|
(2) |
|
图 1 宽阔的大陆架区域 (红色段为实际测线位置) Figure 1 Wide continental shelf (location of real data is shown as red line) |
对各慢度进行预测反褶积去除周期性多次波.假设多次波周期为α，将记录F (τ, p) 通过预测滤波器aτ得到τ+α的估计记录
|
(3) |
其中ri(i=0, …, n－1) 为自相关记录，rα+i(i=0, …, n－1) 为输入与期望输出互相关，λ为白噪因子.求解上述方程组即可得到预测滤波器，将求得的滤波器与原始数据褶积获得预测多次波，将其从原始数据中减去得到Tau-p域去除多次波后的数据.对其进行Tau-p反变换最终获得时空域衰减多次波后的记录fa(x, t)，即：
|
(4) |
上述推导过程可见该方法具有不依赖速度信息、参数简单、计算速度快等优点，但该方法运用于陆架浅水资料存在以下问题：
① 如果数据波场中包含周期与多次波接近的有效波成分，基于多次波周期性的预测反褶积方法对有效波是有损伤的；
② 如果同时出现如图 2a和图 2b所示的全水层周期多次波和peg-leg多次波，单步长的预测反褶积方法是无法压制彻底的.
在波动理论类方法中，波场延拓法基于模型驱动，需要水深及水速等先验信息.而逆散射级数法和反馈迭代方法则为数据驱动，无需先验信息.反馈迭代去除多次波方法中，荷兰Delft大学的DELPHI小组为SRME (Surface-Related Multiple Elimination) 方法的提出及应用作出了杰出的贡献 (Berkhout and Verschuur, 1997; Verschuur and Berkhout, 1997).SRME为现阶段较为常用的波动理论类多次波压制方法.其核心思想为自由表面多次波可以由全波场数据的一次波褶积得到.褶积一次可得到一阶表面多次波，褶积两次可以得到二阶表面多次波 (史文英等, 2013).主要数学表达式为
|
(5) |
其中M为多次波，D为包含一次波和多次波的实际记录，P为一次波，S为震源子波，⊗为褶积符号，又有：
|
(6) |
将式 (5) 代入式 (6) 可以得到一次波估算为
|
(7) |
假设震源子波不变，自由界面的反射系数为-1，第一次迭代时的一次波用实际记录代替，反子波用最小能量假定统计估算，忽略子波项后实际迭代过程可以表示为
|
(8) |
|
(9) |
|
(10) |
式 (10) 中n为迭代次数，f表示自适应滤波反褶积.获得预测的多次波后，利用基于能量准则的自适应相减方法得到去除多次波的数据记录.该方法已被广泛应用于深水数据多次波去除，但应用于浅水地震资料时仍存在困难：
① 该方法要求输入数据为规则采样的全波场数据，这在浅水条件下很难满足.由于最小偏移距的存在必须先对输入道集进行外插至零偏移距.然而由于水深较浅无法区分直达波和海底反射，导致近偏移距处反射估算困难，插值后的海底反射不够准确.此外，海上二维地震资料采集时羽角的存在同样使得数据难以满足规则采样.
② 推导公式 (10) 可以看出，SRME方法是通过迭代的方式来求取各阶多次波.浅水多次波的多阶性使得计算过程需要多次迭代，容易使预测模型引入振幅误差.进而导致原始记录与多次波模型能量不匹配，在自适应相减时存在伤害有效波的风险.
综上所述，预测反褶积和SRME方法在压制浅水自由表面多次波时都存在较大的局限性.为了更好地压制陆架浅水区的自由表面多次波，提高地震数据成像质量，本文采用确定性水层多次波压制方法和SRME组合的压制思路.确定性水层多次波压制方法通过在Tau-p域采用基于水体模型波场延拓的多次波预测算法实现浅水水层多次波的压制，结合SRME方法，能够较好地解决浅水区自由表面多次波压制问题.
2 确定性水层多次波压制方法 2.1 方法原理及实现过程确定性水层多次波压制，顾名思义，用于压制传播路径为图 2a-d类型的水层多次波.主要原理为：假定地震波在水层 (海面和海底) 间再传播一个反射走时，则海底反射变成了一阶水层多次，一阶水层多次波变成二阶水层多次.因此，通过将原始记录向下延拓一个水层走时，再乘以海底反射系数即可获得多次波模型 (假设地震子波在此过程中稳定).获得多次波模型后通过自适应相减的方法从原始记录中减去从而得到去除水层多次波后数据.类似的思想最早出现于20世纪80年代 (Morley and Claerbout, 1983; Berryhill and Kim, 1986; Wiggins, 1988).由于该方法需要水速及水深信息及一维层状介质的假设限制了其使用范围.然而，在压制陆架浅水区自由表面多次波时该方法又具有优势：
(1) 陆架浅水区海底相对平缓无剧烈水深变化，满足近似层状介质假设.
(2) 船载测深设备能够获得同步水深信息并可将其记于地震数据道头.现有测深仪精度较高加上水浅回声能量强，获得的水深数据真实可信.
(3) 根据海水速度与温度、盐度、深度关系可知，浅水区海水速度稳定可近似为常数 (Han et al., 2012)，可将其视为已知.
(4) 该方法利用水体反射周期作为多次波预测算子，从而避免了浅水近偏移距数据缺失对SRME技术的影响.
(5) 该方法有能力衰减与一次波速度差异很小的自由表面多次波成分，且更好地保护一次波.
确定性水层多次波压制方法流程如图 4所示：首先对输入炮集数据利用高精度Radon变换外插近偏移距道获得全偏移距数据 (Kabir and Verschuur, 1995; 黄新武等, 2003)，插值结果如图 4b(虚线框中为插值道).对数据进行Tau-p正变换，将变换后Tau-p内各偏移距对应的水层走时作为时移静校正量，使记录整体下移一个水层反射时间.可以看到，移动前后同一位置的同相轴极性相反，如图 4c和图 4d，表明数据又进行了一次水层反射.然后将时移后的记录与反射系数相乘，即可获得预测的水层多次波模型1，如图 4e.此时就可以将原始数据与水层多次波模型进行自适应相减获得去水层多次波的数据.同样也可以将相减后的数据重复 (d)、(e) 两步，获得水层多次波预测模型2.根据实际数据去水层多次波的效果，该预测模型过程可迭代多次.将得到的多个预测模型相加获得最终的水层多次波预测模型，如图 4g所示.利用Tau-p反变换将多次波模型转换到t-x域，与插道后的原始数据进行自适应相减获得去除水层多次波后的记录，去除外插近道后输出数据如图 4i.
|
图 4 确定性水层多次去除方法流程 Figure 4 The process of DWD |
水层多次波的预测需要较为精确的海底反射时间.将经过水速动校正的近道剖面进行自相关，在自相关剖面上拾取海底反射同相轴 (图 5)，自相关剖面能够避免系统延迟等因素的影响.在拾取前应对数据进行枪缆深度校正，将炮点和接收点统一到海平面.同时，将拾取的海底反射时间与通过测深仪获得的结果进行对比，以验证拾取的准确性.
|
图 5 自相关剖面确定海底反射时间 Figure 5 Reflection time of sea bottom obtained from autocorrelation profile |
将T-X域的道集数据转换到Tau-p域后，需要计算各偏移距处的海底反射时间，计算公式为
|
(11) |
其中tw为零偏移距处的海底双程走时，p为经Tau-p变换后的斜率值，vw为海水速度，x为参考偏移距.根据公式 (11) 得到的Tau-p域各偏移距处海底反射时间如图 6.
|
图 6 Tau-p域各偏移距海底反射时间 Figure 6 Reflection time of sea bottom with offset in Tau-p domain |
为了获得多次波预测模型，需用向下延拓后的数据乘以海底反射系数.海水的速度和密度相对稳定，海底反射系数与海底下第四纪沉积物的速度和密度参数有关，海底反射系数取值范围一般为0.3~0.6之间.
3 实例分析实际多道地震数据来自东海陆架盆地南部浅水区，测线位置见图 1红色线段，平均水深70 m左右.测线沿近岸自西向东跨越瓯江凹陷、雁荡低凸起到达闽江凹陷，由于水浅且海底底质较硬，大量发育多次波，其中自由表面多次波为主要的多次波成分.在去多次波前进行了涌浪、水鸟、电缆噪声的去除，噪声去除过程中注意不破坏多次波的传播特性.多次波压制时对SRME方法、Tau-p域反褶积方法及DWD和SRME组合方法进行了测试对比.图 7a为典型的原始单炮记录及自相关，自相关窗口取800~2000 ms，可见该数据存在大量短周期多次波.图 7b为只进行SRME后的炮集记录，图 7c为只进行SRME后的炮集记录.可以看到，针对有效反射下跟随的水层多次波，SRME方法去除效果不佳，反褶积方法能够去除大量的短周期自由表面多次波，但无论是从单炮、自相关还是Tau-p域数据都可以看到DWD和SRME组合方法的多次波压制效果更好且一次波特征得到更好的保留.
|
图 7 自由表面多次波去除对比 (a) 原始单炮；(b) 仅做SRME；(c) 仅做Tau-p反褶积；(d) DWD+SRME；(e)-(h) 为单炮对应的Tau-p域记录. Figure 7 Methods comparison of surface-related multiple removal (a) Raw data; (b) SRME only; (c) Deconvolution in Tau-p domain only; (d) DWD+SRME; (e)-(h) Data in Tau-p domain. |
对数据采用DWD方法压制水层多次波后需在道集上切除近海底反射区域，再利用SRME方法来压制传播路径如图 2e所示的长周期自由表面多次波.切除海底反射的目的在于防止SRME重复预测水层多次波.DWD和SRME组合去除自由表面多次波前后的首道剖面如图 8所示.图 8a可见由于水浅造成多阶的自由表面多次波震荡干涉形成大量的异常振幅条带，使用组合压制方法后异常振幅条带基本消除.在首道剖面上的200~1200 ms及2100~3100 ms时窗内进行自相关分析，用以评价浅层和中深层的多次波压制效果.结果分析表明DWD前后两个时窗内的短周期多次波得到明显压制，SRME后浅层时窗差异不大，而中深层时窗内的压制效果进一步提升，这也体现了SRME方法能够针对长周期自由表面多次波进行压制.图 9为自由表面多次波压制前后的叠加剖面局部放大图.去多次波前的叠加剖面存在明显的多轴现象，进行DWD和SRME方法组合压制后，多轴现象消失，有效波特征得以凸显.图 10为组合压制前后的频谱变化，可以看到无论无论是中浅层还是中深层，多次波压制后的频谱 (蓝色线) 陷波得到大幅补偿，这也充分证明了该组合方法对于多次波压制的有效性.
|
图 8 组合压制自由表面多次波前后的首道剖面 (a) 去多次波前；(b) DWD后；(c) DWD+SRME后，点线框为自相关区域，自相关步长200 ms. Figure 8 First trace profile before and after surface-related multiple suppression (a) Before demultiple; (b) After DWD; (c) After DWD+SRME, autocorrelation area showed in dotted line, the autocorrelation step is 200 ms. |
|
图 9 DWD和SRME组合压制多次波前后叠加剖面局部对比 (a) 去多次波前；(b) DWD+SRME后. Figure 9 Partial zoom-in of stacking profile before and after DWD+SRME (a) Raw data; (b) After DWD+SRME. |
|
图 10 DWD和SRME组合压制多次波前 (红色) 后 (蓝色) 频谱 Figure 10 Spectrum before and after DWD+SRME |
|
图 11 速度谱 (a) 原始数据；(b) DWD后；(c) SRME后；(d) RADON变换去多次波后). Figure 11 Velocity spectrum (a) Raw data; (b) After DWD; (c) After RADON transform; (d) E-muitiple. |
自由表面多次波压制后，有效波速度能量团更加聚焦，深部能量团得以显现，有利于提高速度拾取的准确性.在拾取较为准确的速度后，利用基于速度差异的高精度拉东变换去多次波方法去除剩余的多次波能量，就完成了该浅水数据的多次波压制.
有效压制多次波是海洋地震成像的基础和关键.完成多次波压制后通过多轮的叠前时间偏移迭代，最终获得的偏移剖面能够正确地反应东海陆架盆地南部各地质单元的接触关系 (图 12).相比于原始剖面凹陷基底的“多轴”现象，进行多次波压制和偏移成像后凹陷基底反射连续，波阻关系清晰，尤其是偏移剖面揭示出的部分中生代地层有效反射，给该海域中生代油气资源调查提供了可靠的基础资料.
|
图 12 叠加剖面 (a) 原始叠加；(b) 去多次波后；(c) 偏移剖面. Figure 12 Stacking profile (a) Raw data; (b) After multiple removal; (c) PSTM section. |
我国大陆架浅水区分布广，油气资源潜力大，提升该区地震资料成像品质是开展后续研究工作的现实需求.多次波压制是海洋地震资料处理的基础和关键，而陆架浅水区直达波和海底反射无法区分，自由表面多次波能量强、周期短、阶次多的特点使得常用的基于反馈迭代的SRME方法存在多次波预测不准、容易引入异常振幅的缺陷，同时预测反褶积方法存在损伤有效波的危险.
4.2利用基于模型驱动的确定性水层多次波压制方法来压制水层多次波，使用SRME方法在切除海底的数据上压制长周期自由表面多次波.使用该方法组合能够较好地解决浅水区自由表面多次波压制问题.来自东海陆架盆地浅水区的数据处理结果表明了该组合压制方法的有效性.
4.3一维层状介质假设使确定性水层多次波压制方法更适用于海底平缓的陆架区，进一步的研究应拓展其对于崎岖海底的适用性.
致谢 感谢斯伦贝谢北京技术支持中心李蓉在实际数据处理中给予的帮助，感谢“发现2号”船全体成员在多道地震数据采集过程中付出的努力.| [] | Berkhout A J, Verschuur D J. 1997. Estimation of multiple scattering by iterative inversion, Part Ⅰ:Theoretical considerations[J]. Geophysics, 62(5): 1586–1595. DOI:10.1190/1.1444261 |
| [] | Berryhill J R, Kim Y C. 1986. Deep-water peg legs and multiples:Emulation and suppression[J]. Geophysics, 51(12): 2177–2184. DOI:10.1190/1.1442070 |
| [] | Dragoset W H, JeričevićŽ. 1998. Some remarks on surface multiple attenuation[J]. Geophysics, 63(2): 772–789. DOI:10.1190/1.1444377 |
| [] | Fan B C, Xie T. 2012. The application of partial SRME in shallow sea[J]. Journal of Oil and Gas Technology (in Chinese), 34(12): 73–75. |
| [] | Fan C H, Wang Y M, Liu H. 2013. A review of shallow continental shelf sedimentary characteristics[J]. Geological Science and Technology Information (in Chinese), 32(2): 29–34. |
| [] | Han F X, Sun J G, Wang K. 2012. The influence of sea water velocity variation on seismic traveltimes, ray paths and amplitude[J]. Applied Geophysics, 9(3): 319–325. DOI:10.1007/s11770-012-0344-2 |
| [] | Huang X W, Sun C Y, Niu B H, et al. 2005. Surface-related multiple prediction and suppression based on data-consistence:A theoretical study and test[J]. Chinese J. Geophys. (in Chinese), 48(1): 173–180. DOI:10.3321/j.issn:0001-5733.2005.01.023 |
| [] | Huang X W, Wu L, Niu B H. 2003. Reconstruction of seismic traces by parabolic radon transform[J]. Journal of China University of Mining & Technology (in Chinese), 32(5): 534–539. |
| [] | Jiang L. 2003. Exploration status and perspective of petroleum resources in east China sea shelf basin[J]. China Offshore Oil and Gas (Geology) (in Chinese), 17(1): 1–5. |
| [] | Jin M X, Yi S C, Wang Z. 2014. The analysis and application of multiple wave forward modeling and the influencing factors of SRME technology[J]. Geophysical and Geochemical Exploration (in Chinese), 38(3): 566–570. DOI:10.11720/wtyht.2014.3.26 |
| [] | Kabir M N N, Verschuur D J. 1995. Restoration of missing offsets by parabolic Radon transform[J]. Geophysical Prospecting, 43(3): 347–368. DOI:10.1111/gpr.1995.43.issue-3 |
| [] | Liu G D, Chen J. 2005. Potential analysis of petroleum explortation in residual basins of the China sea[J]. Progress in Geophysics (in Chinese), 20(4): 881–888. DOI:10.3969/j.issn.1004-2903.2005.04.001 |
| [] | Liu J, Wu S Y, Gao J Y, et al. 2016. An effectiveness analysis of multiple depression technique in the Shallow water of the central uplift in South Yellow Sea[J]. Geophysical and Geochemical Exploration (in Chinese), 40(3): 568–577. DOI:10.11720/wtyht.2016.3.20 |
| [] | Morley L, Claerbout J. 1983. Predictive deconvolution in shot-receiver space[J]. Geophysics, 48(5): 515–531. DOI:10.1190/1.1441483 |
| [] | Shi W Y, Li L, Yuan Q S, et al. 2013. The application of series SRME method to multiple wave attenuation in southwest Weizhou area[J]. Geophysical and Geochemical Exploration (in Chinese), 37(5): 911–915. DOI:10.11720/j.issn.1000-8918.2013.5.28 |
| [] | Shi Y, Xing X L. 2011. Investigation progress on surface-related multiple suppression:Review and outlook[J]. Progress in Geophysics (in Chinese), 26(6): 2046–2054. DOI:10.3969/j.issn.1004-2903.2011.06.020 |
| [] | Song J W, Verschuur D J, Chen X H. 2014. Research status and progress in multiple elimination[J]. Progress in Geophysics (in Chinese), 29(1): 240–247. DOI:10.6038/pg20140134 |
| [] | Verschuur D J, Berkhout A J. 1997. Estimation of multiple scattering by iterative inversion, Part Ⅱ:Practical aspects and examples[J]. Geophysics, 62(5): 1596–1611. DOI:10.1190/1.1444262 |
| [] | Weglein A B. 1999. Multiple attenuation:An overview of recent advances and the road ahead[J]. The Leading Edge, 18(1): 40–44. DOI:10.1190/1.1438150 |
| [] | Wiggins J W. 1988. Attenuation of complex water-bottom multiples by wave-equation-based prediction and subtraction[J]. Geophysics, 53(12): 1527–1539. DOI:10.1190/1.1442434 |
| [] | Xiao E L, Chen Y, Wan H, et al. 2010. Surface-related multiple elimination on marine seismic data[J]. Progress in Geophysics (in Chinese), 25(3): 1057–1064. DOI:10.3969/j.issn.1004-2903.2010.03.046 |
| [] | Zhang Z Z, Li S F, Fang Z Y, et al. 2016. The combination of SWMA and SRME for surface-related multiple attenuation——case study in the a block of the Zhujiangkou basin[J]. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 38(2): 244–251. DOI:10.3969/j.issn.1001-1749.2016.02.16 |
| [] | Zhao C L, Ye Y M, Yao G S, et al. 2013. Prediction deconvolution in linear radon domain on the application of ocean multiples attenuation[J]. Progress in Geophysics (in Chinese), 28(2): 1026–1032. DOI:10.6038/pg20130256 |
| [] | 范宝仓, 谢涛. 2012. 部分去除表面多次波技术在浅水海域的应用[J]. 石油天然气学报, 34(12): 73–75. DOI:10.3969/j.issn.1000-9752.2012.12.016 |
| [] | 范春花, 王英民, 刘豪. 2013. 浅海陆架沉积特征研究综述[J]. 地质科技情报, 32(2): 29–34. |
| [] | 韩复兴, 孙建国, 王坤. 2012. 海水速度变化对地震波走时、射线路径及振幅的影响[J]. Applied Geophysics, 9(3): 319–325. DOI:10.1007/s11770-012-0344-2 |
| [] | 黄新武, 孙春岩, 牛滨华, 等. 2005. 基于数据一致性预测与压制自由表面多次波——理论研究与试处理[J]. 地球物理学报, 48(1): 173–180. DOI:10.3321/j.issn:0001-5733.2005.01.023 |
| [] | 黄新武, 吴律, 牛滨华. 2003. 基于抛物线拉东变换的地震道重构[J]. 中国矿业大学学报, 32(5): 534–539. |
| [] | 姜亮. 2003. 东海陆架盆地油气资源勘探现状及含油气远景[J]. 中国海上油气 (地质), 17(1): 1–5. |
| [] | 金明霞, 易淑昌, 王征. 2014. 多次波正演与SRME技术影响因素分析及应用[J]. 物探与化探, 38(3): 566–570. DOI:10.11720/wtyht.2014.3.26 |
| [] | 刘光鼎, 陈洁. 2005. 中国海域残留盆地油气勘探潜力分析[J]. 地球物理学进展, 20(4): 881–888. DOI:10.3969/j.issn.1004-2903.2005.04.001 |
| [] | 刘俊, 吴淑玉, 高金耀, 等. 2016. 南黄海中部浅水区多次波衰减技术及其效果分析[J]. 物探与化探, 40(3): 568–577. DOI:10.11720/wtyht.2016.3.20 |
| [] | 史文英, 李列, 袁全社, 等. 2013. 串联SRME在涠西南地区多次波衰减中的应用[J]. 物探与化探, 37(5): 911–915. DOI:10.11720/j.issn.1000-8918.2013.5.28 |
| [] | 石颖, 邢小林. 2011. 表面多次波压制的研究进展:回顾与展望[J]. 地球物理学进展, 26(6): 2046–2054. DOI:10.3969/j.issn.1004-2903.2011.06.020 |
| [] | 宋家文, VerschuurD J, 陈小宏. 2014. 多次波压制的研究现状与进展[J]. 地球物理学进展, 29(1): 240–247. DOI:10.6038/pg20140134 |
| [] | 肖二莲, 陈瑜, 万欢, 等. 2010. SRME多次波衰减方法在海洋地震资料中的应用[J]. 地球物理学进展, 25(3): 1057–1064. DOI:10.3969/j.issn.1004-2903.2010.03.046 |
| [] | 张治忠, 李三福, 方中于, 等. 2016. SWMA与SRME组合衰减自由界面多次波技术——在珠江口盆地A区的应用[J]. 物探化探计算技术, 38(2): 244–251. DOI:10.3969/j.issn.1001-1749.2016.02.16 |
| [] | 赵昌垒, 叶月明, 姚根顺, 等. 2013. 线性拉东域预测反褶积在海洋多次波去除中的应用[J]. 地球物理学进展, 28(2): 1026–1032. DOI:10.6038/pg20130256 |