沈金松, 汪 轩，魏帅帅，李 曼

1.中国石油大学(北京)地球物理与信息工程学院，北京 102249 2.油气资源与探测国家重点实验室，北京 102249 3.中国石油集团公司物探重点实验室，北京 102249



浅水域可控源海洋电磁测量中空气波的渐近表达与波场分离方法

沈金松1,2,3, 汪 轩1，魏帅帅1，李 曼1

1.中国石油大学(北京)地球物理与信息工程学院，北京 102249 2.油气资源与探测国家重点实验室，北京 102249 3.中国石油集团公司物探重点实验室，北京 102249

可控源海洋电磁(MCSEM)勘探中空气波对海底电磁响应的影响已为业界所重视，它是由水平电偶极子源发射的电磁信号沿空气-海水界面传播与来自海底地层的有效信号相互作用产生的。在浅水域勘探时，空气波淹没来自地层的有效信号，使浅水域MCSEM实现油气层识别产生困难。基于层状介质模型，采用电磁场的模式分解理论导出了半空间电阻率模型的空气波表达式，利用该式将空气波在海水层的传播近似用海水-空气界面与海底地层之间的多次鸣震表达。结合源和接收器两边电磁场的传播特征，导出了有限水深时空气波近似关系，用于近似模拟空气波响应。为了对比不同方法压制空气波的效果，基于上行波场和下行波场的分解方法，获得了含油气高阻储层上上行波的异常幅度增大数倍的结果，显示了波场分解方法压制空气波的良好效果。最后，利用不含油气层的背景模型和含油气储层模型电磁响应的数值模拟结果，比较空气波渐近表达和波场分离2种空气波压制方法可知：对于水平层状模型后者效果更好；前者可适用于崎岖海底地层的数据处理，后者只适用于水平海底地层。

浅水域；海洋可控源电磁勘探；鸣震近似；波场分离；空气波压制

0 引言

可控源海洋电磁(MCSEM)勘探实现油气储层识别和描述的主要原理在于电磁能量在海底地层中传播时遇到不同电阻率层表现出衰减和相位的差异[1-4]。油气储层电阻率比围岩高几到几十倍，电磁能量在储层中衰减较慢，因此，储层中的电磁场幅度比围岩中大[5-6]。然而，在浅水域实施MCSEM勘探时，发射源在空气-海水界面激发的空气波在大偏移距的海底接收器响应信号中占主导，使来自海底地层的有效信号淹没于强空气波信号中，难以区分高阻层地层引起的异常[7-8]。相应地，在深水域，由于空气波在海水层中的双程衰减，对接收信号的影响可以忽略。因此，MCSEM方法在海洋油气勘探中的成功应用目前主要限于深水域(水深>300 m)[9-10]。

图1给出了海洋电磁信号的基本传播路径和模式，其中的“空气波”是受弹性波传播的启发，最初是由Chave和Cox[11]命名的，主要指由源到海水面垂直扩散的信号产生、在空气层以波的形式传播、通过海水层扩散返回海底的信号，它在空气层传播时几乎没有衰减。与深水域相比，在浅水域时接收器记录的空气波要强得多，这是浅水域海洋电磁勘探的困难所在。

图1 水平电偶极子源发射低频电磁波在海水表面、水中和海底地层传播示意图Fig.1 Schematic diagram of propagation paths of low-frequency electromagnetic wave emitted from a horizontal electric dipole to the sea surface, seawater and undersea stratum

关于“空气波”性质的研究，早在50年前Wait[12]就曾论述过，后来Baos[13]、Bannister[14]和King等[15]也已开展过研究，只是他们没有用“空气波”的名称而已，他们将之称为“侧面波”并给出了半空间水层之上空气层中的空气波表达式。近年来，Constable等[1]、MacGregor等[4]、Løseth等[16]基于层状介质的TE(横电)和TM(横磁)模式分解理论，分析了空气波产生的机制和影响大小。Um等[5]用电流密度剖面分析了不同海水深度上空气波能量占主导的临界偏移距位置，得到了垂直电场分量不受空气波影响的重要认识。Ziolkowski等[17]提出了时间域的信号分离等多种压制空气波方法。Amundsen等[8]基于“空气波”的波动性质，将接收的电磁信号分离为上行和下行的电磁波场，通过去除下行波场实现空气波的压制。Lu等[18]提出了应用数值模拟的方法估计空气波信号，再从接收信号中减去的方法，减弱空气波的影响。国内，殷长春等[19]、沈金松等[20]和刘长胜等[21]基于一维层状介质的电磁场响应，分析了空气波随水深变化对接收器响应的影响和相应的压制方法。

本研究基于电磁场的模式分解理论导出海底接收信号中空气波与空气层、海水层、海底地层的相互作用关系，用空气波产生的物理机制解释空气波与海底地层信号的复杂耦合特征。在Nordskag等[22]研究的基础上，利用有限水深模型，建立空气波在海底地层和海水界面之间交混回响衰减模式，导出有限水深时空气波分量的近似表达式，构建模拟电磁场相减压制空气波的方法。同时，还依据Amundsen等[8]的波场分离方法，推导海洋电磁测量数据的上行和下行波场分离关系，讨论对剖面数据实现空气波压制的效果。通过层状背景模型和高阻含油气储层模型的模拟数据比较2种空气波压制方法的效果，分析各自的优缺点。

1 电磁场基本方程

由频率域Maxwell方程，低频水平电偶极子(HED)的电磁场可表示为[23]

,

在层状均匀介质中，方程(1)的电场和磁场可分解为横电和横磁2种模式及相应的矢量位函数表示[20,23]：

ETE=-(

其中：ETE、ETM和HTE、HTM分别表示横电(TE)和横磁(TM)模式的电场和磁场强度；F和A为2种模式对应的位函数。F和A满足如下非均匀Helmholtz方程[23]：

对于图2所示的无限水深单界面的情况，若以海水-空气界面为z=0，向下为正，利用电磁场连续性条件，即可求得海水层中HED源激发的空气和海水层中的电磁位[20]：

图2 无限水深MCSEM分析介质模型Fig.2 Model of infinite water depth for MCSEM analysis

2 海水半空间模型中空气波的渐近表达式

为更好地分析源信号与空气界面、海水层的相互作用，以半空间海水层模型(图2)的电场径向分量为例对空气波的表达式进行分析。根据式(6)——式(9)，TM和TE两种模式的径向电场分量可表示如下[20]：

式中：k0≈0m-1；k1=(iωμ0σ1)1/2；ω是圆频率；θ是源的方位角；zs和zr分别为源和接收器的深度。

3 有限水深模型中空气波的近似表达与压制

3.1 有限水深模型中空气波的多次鸣震近似

基于远场近似关系式(13)知道，在半空间海水层的情况下，空气波主要与4个因素相关，即海水层的电导率σ1、源与接收器的深度zs和zr，以及源与接收器之间的偏移距ρ。对于有限水深，Bannister[14]对侧面电磁波的研究表明，水平电偶极子源的TE模式分量基本由源垂直向上传播，在海水与空气界面感应产生空气波，在海底地层和海水界面之间形成一系列快速衰减的交混回响信号。同样，由源激发的垂直向下传播的TE模式信号到达海底也形成向上反射的交混回响衰减信号。类似地，在接收器一侧，最初垂直向下的TE模式信号在海底地层和海水面之间也发生交混回响。然而，由于海水层电导率较大，其趋肤深度较小，在海水层中对空气波分量有较大贡献的交混回响次数不会太大。因此，对于有限水深的空气波传播，笔者从发射器一侧、海水空气界面和接收器一侧3个路径上电磁信号的传播和交混回响来分析。

图3给出了源激发的空气波在海水和海底界面之间反射的示意图。海水的电导率为σ1，深度为zb处海底地层的电导率为σ2。图3a给出由zs源处激发向上扩散的信号到达z=0处的海水面。由式(13)第一个指数项，发射器一侧的影响可以近似地表示为[24]

式中，复波数k1决定了信号的衰减。

a.来自源的初始上行波场；b, c.由初始上行波场感应的多次交混回响；d.初始下行波场和海底反射在海水面感应的空气波；e,f.由下行波场诱发的多次交混回响。图3 源侧反射和交混回响信号Fig.3 Signals of the reflection and reverberation near the source side

图3d为源产生的信号向下扩散再在海底反射向上传播，并在海水面诱发空气波。参考式(14)在z=0处的空气波可用下式近似：

图3b和图3c给出了图3a中描述的向上传播的源信号在空气波产生前所发生的一次或二次鸣震，对应的近似关系分别为[24]：

同样，图3e和图3f给出了图3d中最初向下传播的源信号在产生空气波前发生了一次或二次鸣震，对应的关系分别为[24]：

exp(-ik1zs)·R·

exp(-ik1zs)·R·

空气波交混回响过程的无限持续所形成的总信号为[7,24]

其中，

S={1+R·exp[2ik1(zb-zs)]}[1+R·

exp(2ik1zb)+R2·exp(4ik1zb)+…+

Rn·exp(2nik1zb)+…]={1+R·

exp[2ik1(zb-zs)]}/{1-R·exp(2ik1zb)} 。

S是一个滤波器，它表示源一侧产生的所有空气波分量。两个花括号中的项依次代表来自源的初始信号和海底反射信号。

类似地，图4中给出了空气波分量由海底反射和鸣震在接收器一侧的改造作用。在图4a中，空气波分量到达zr处的接收器；然后，它在海底反射又回到接收器。这一过程可以近似地描述为[24]

图4b和图4c分别考虑海水层中的一次和二次鸣震，其数学表达为[24]

exp(-ik1zr){1+R·

exp(-ik1zr){1+R·

最后，考虑无限多次鸣震的衰减，接收器一侧的反射和鸣震级数可表示为[24]

对于有限水深模型的空气波，除了式(13)中的项外，还包括了源一侧和接收器一侧的反射和鸣震：

方程(25)假设了HED源激发的TE模式场最初是在海水层中垂直向上和向下传播的，且也在海水层中激发一系列垂直传播并以指数形式衰减的反射和鸣震级数，每个TE模式分量在海水表面诱发幅度逐渐衰减的空气波分量。在接收器一侧，也发生与源一侧类似的反射和鸣震。

a.接收器记录的来自空气波的下行波场和海底反射；b,c.接收器一侧经过多次鸣震改造的空气波下行波场和海底反射。图4 接收器一侧反射和交混回响信号Fig.4 Signals of the reflection and reverberation near the receiver side

3.2 崎岖海底和海水电导率横向变化模型中空气波的近似

设源一侧的水深为zbs，海底的反射系数为Rs：

式中：σ1s为源位置海水的电导率；σ2s为zbs处海底地层的电导率。将Rs代入式(20)中的S[7,22,24]：

同样，在接收器一侧，若海水深为zbr及海底反射系数为Rr：

式中：σ1r为接收器位置海水的电导率；σ2r为zbr处海底地层的电导率。用zbr和Rr代替方程(24)中的zb和R，可以得到新的R[24]：

利用式(26)和式(27)的修正关系，可在源和接收器附近为局部平坦崎岖海底(图5)实现空气波压制。它在除海底反射外局部崎岖区域不发生较大散射时，可以较好地模拟崎岖海底的空气波。需要注意的是，前文分析的鸣震近似方法主要考虑了发射和接收机附近电性分布变化的影响，同时在式(26)和式(27)中也考虑了源和接收器处的海水深度变化，但没有考虑发射接收间距较大情况下电性变化和海底起伏较大时对空气波的影响，这也是本方法有待改进的致命缺陷之一。

图5 源和接收器在深度存在较大差异的崎岖海底Fig.5 Rugged seabed model with different depth levels of source and receiver positions

3.3 鸣震近似关系空气波压制效果

为分析前述空气波近似关系(式(27))的空气波压制效果，利用有限水深模型(图6a)计算的近似空气波的可探测性与无限水深模型(图6b)进行对比，分析交混回响近似。

a.100 m水深五层和三层模型；b.15 000 m水深五层和三层模型。图6 空气波压制效果模型示意图Fig.6 Models for indication of the air wave mitigation effect

图6a中左右两模型的差异为储层的有无：左图是五层模型，存在一层电导率为0.02 S/m的储层；右图是三层模型，无储层。图6b中两模型与图6a的差异为海水层的厚度，图6b中两模型海水深度都为15 000 m，用于近似模拟无限水深模型，计算无空气波时的电场幅度。

a.水深100 m；b.五层模型。图7 鸣震近似压制空气波结果对比Fig.7 Comparison of the air wave mitigation results by using reverberation approximate expression

图7a给出了100 m水深时五层模型和三层模型的原始电场径向分量曲线与鸣震近似关系压制空气波后的电场径向分量曲线。由图7a看到：2种模型测量的原始数据(未压制)只在4 000～6 000 m的偏移距上出现微小的幅度差，显示了油气储层的存在；用鸣震近似关系压制后，电场幅度在3 000～5 000 m的偏移距上得到微小提升。

图7b给出了五层模型在15 000、100 m水深时的原始电场响应与鸣震近似关系压制空气波后100 m水深的标准化电场幅度曲线(分别用各自的三层模型标准化)。从图7b中看出：由于空气波的影响，水深100 m时，原始电场曲线(未压制)难以显示油气储层，而经过鸣震近似关系压制空气波后，在3 000～5 400 m的偏移距上显示了小幅度提升，但在5 400～10 000 m又呈现小幅度下降，该空气波压制方法效果仍有待改善；将水深100 m鸣震近似后的标准化幅度曲线与水深15 000 m的原始电场响应(未压制)的标准化幅度曲线对比知道，前文实现的鸣震近似空气波压制方法仍然难以较好地消除空气波的影响，这主要是由于前文的鸣震近似方法未考虑占比较大的非垂直传播散射场的影响。

4 上行波与下行波分离及空气波压制方法

4.1 散射电磁波的上行与下行波场分离

前文的鸣震近似压制空气波的方法主要将空气波的传播近似看成海水层中垂直传播和空气层中的测量传播的波场，这是将散射场近似成波场的考虑方法。类似地，将散射场近似看成向上和向下衰减的波场，Løseth等[25]和Andréis等[26]提出了将接收器记录的散射场分解成上行和下行衰减场，进而压制空气波的方法。下面以他们[25-26]的波场分解理论为基础，分析上行和下行散射场分解压制空气波的方法。

假设接收器记录的电磁响应是在某个电导率、磁导率和介电常数为已知常数的平面上获得的，设x=(x1,x2,x3)=(x,y,z)，记深度轴向下为正，与式(1)相似，各向同性无源介质中电磁场方程(1)退化为[27-28]

×E=-iωμH ,

其中：

Eh=(E1,E2)T，Hh=(-H2,H1)T。

为了导出电场和磁场之间的关系以及它们与上行和下行分量之间的关系，需要找出系统矩阵A的特征函数，即如下方程的解：

其中：λN是特征值；而lN是对应λN的特征向量。4个特征值之间有如下关系：

它们可以组合成4×4的对角阵：

Λ=diag[λ1,λ2,λ3,λ4]=

类似地，将特征向量作为列向量也可以形成4×4的矩阵：

那么方程(31)的特征解可表示为

在介质空间的均匀区域，系统矩阵A是常数，从而特征矢量矩阵L与x3无关。对于电磁场矢量b的一阶微分方程(29)可以转化成下面的一阶微分方程：

引入4×1的列向量w=L-1b，方程(36)对每个w元素变为对上行波场和下行波场解耦的微分方程，分别代表2个上行波场和2个下行波场的微分方程。因此,w中的元素代表了上行波和下行波场。对特征向量作如下标准化可以得到

其中，L-1的作用是将电磁场矢量分解成上行波和下行波分量，因此，L-1称为分解矩阵。

事实上，分解只是物理变量的简单线性变换，它是一个稳定的计算过程。由b=Lw知道，L的作用是将上行和下行衰减场分量组合成电磁场矢量，L称为组合矩阵。在波动理论中，λN代表相位慢度，lN表示极化矢量[29]，将电磁场矢量分解成上行和下行衰减场的具体实现过程见文献[8]。

4.2 电磁波场分离压制空气波的效果

a.水深100 m；b.五层模型。图8 上行波和下行波分离压制空气波结果对比Fig.8 Comparison of the air wave mitigation results by using decomposition of up going and down going wave fields

为了考察电磁波场分离压制空气波的效果并与前文的鸣震近似压制空气波效果对比，这里仍然沿用图6a所示的层状介质模型。偶极子源发射0.25 Hz的正弦波。

2种模型下，利用一维频域模拟算法计算得到的E1幅度曲线如图8a所示。从图8a看到:2种模型测量的原始数据(未压制)只在3 000～6 500 m的偏移距上出现微小的幅度差，显示了油气储层存在；上、下行波场分离与空气波压制后，电场幅度在空气波较强的大偏移距上得到一定提升。

图8b给出了五层模型在15 000、100 m水深的原始电场响应(未压制)与水深100 m时上下行波分离压制后的标准化电场幅度曲线。从图8b中看出：水深100 m时，由于空气波的影响，原始电场曲线(未压制)难以显示油气储层，而经过下行波分离和压制后，上行波场幅度有所提升;在3 000～8 500 m的偏移距上显示了小幅度异常，显示了油气储层的存在。然而，将水深100 m时上下行波分离压制后的标准化幅度曲线与水深15 000 m时原始电场(未压制)的标准化幅度曲线对比可知，前文描述的上行、下行波场分离的空气波压制方法仍然难以完全压制空气波的影响。这主要是由于前文的分离方法难以将与地下地层信号相互耦合的空气波影响消除，而且这些影响在空气波响应中占较大比重。

5 结论和认识

基于Maxwell方程和模式分解理论导出了无限水深时空气波的完整表达式，据此分析了空气波的作用机理和近似表达式。由无限水深空气波的完整表达式知道，在MCSEM测量中空气波分量主要来自于HED源场的TE模式，它几乎从源垂直向上及向下传播；在海水面，上行波信号产生空气波，在海水面和海底之间激发一系列快速衰减的鸣震信号，其中每个信号在海水面又激发额外的空气波信号。从激发源一侧和接收器一侧空气波反射和鸣震的近似数学表达出发，建立了层状介质和崎岖海底的空气波近似耦合关系，给出了可用于崎岖海底空气波近似分析的方法。数值模拟结果显示，鸣震近似空气波模拟方法，只是模拟了垂直入射和反射的空气波部分，它们只占空气波中的极小比例，因此，空气波压制效果不甚理想，若期望进一步改善效果，只有考虑所有非垂直入射和与地下地层信号耦合的多次空气波才可能实现。另外，基于电磁场上行波和下行波分解的方法，虽然从理论上可以较好地压制空气波，但数值模拟结果显示，该方法对空气波压制及提升有效信号的比例仍然有限。这主要是由于上行波场中仍然包含很大部分空气波与地层信号相互干涉的多次波场。对比文中2种压制空气波方法的应用效果发现，对于实际测量的情况，研究如何分离及消除空气波与地层信号相互干涉的多次波场，而尽可能不损伤地层有效信息的压制空气波的方法，才是解决浅水域MCSEM勘探的根本出路。

[1] Constable S, Weiss C J. Mapping Thin Resistors and Hydrocarbons with Marine EM Methods: Insights from 1D Modeling[J]. Geophysics, 2006, 71(2): 43-51.

[2] Eidesmo T, Ellingsrud S, MacGregor L M, et al. Sea Bed Logging (SBL): A New Method for Remote and Direct Identification of Hydrocarbon Filled Layers in Deepwater Areas[J]. First Break, 2002, 20(3):144-152.

[3] Ellingsrud S, Eidesmo T, Johansen S, et al. Remote Sensing of Hydrocarbon Layers by Sea Bed Logging-SBL: Results from a Cruise Offshore Angola[J]. The Leading Edge, 2002, 21(10):972-982.

[4] MacGregor L M, Thomkins M, Andreis D. Electromagnetic Surveying for Hydrocarbon Reservoirs: U K,20050077902[P]. 2005-04-14.

[5] Um E S, Alumbaugh D L. On the Physics of the Marine Controlled-Source Electromagnetic Method[J]. Geophysics, 2007, 72(2): 13-26.

[6] Chave A D, Constable S C, Edwards R N. Electrical Exploration Methods for the Seafloor[C]//Nabighian M N.Electromagnetic Methods in Applied Geophy-sics. Tulsa: Society of Exploration Geophysicists, 1991:931-966.

[7] Chen J, Alumbaugh D L. Three Methods for Mitigating Airwaves in Shallow Water Marine Controlled-Source Electromagnetic Data[J]. Geophysics, 2011, 76(2): 89-99.

[8] Amundsen L, Løseth L, Mittet R, et al. Decomposition of Electromagnetic Fields into Up Going and Down Going Components[J].Geophysics,2006, 71(5):211-223.

[9] Constable S. Ten Years of Marine CSEM for Hydrocarbon Exploration[J]. Geophysics, 2010, 75(5):7567-7581.

[10] Maaø F A, Nguyen A K. Enhanced Subsurface Response for Marine CSEM Surveying[J]. Geophysics, 2010,75(3):7-10.

[11] Chave A D, Cox C S. Controlled Electromagnetic Sources for Measuring Electrical Conductivity Beneath the Oceans: 1: Forward Problem and Model Study[J]. Journal of Geophysical Research, 1982, 87(B7): 5327-5338.

[12] Wait J R. The Electromagnetic Fields of a Horizontal Dipole in the Presence of a Conducting Half-Space[J]. Canadian Journal of Physics, 1961, 39(7):1017-1027.

[14] Bannister P R. New Simplified Formulas for ELF Sub-surface-to-Subsurface Propagation[J]. IEEE Journal of Oceanic Engineering, 1984, 9(3):154-163.

[15] King R W P, Owens M, Wu T T. Lateral Electromagnetic Waves: Theory and Applications to Communications, Geophysical Exploration, and Remote Sensing[M]. Newyork: Springer-Verlag Inc,1992.

[16] Løseth L O, Amundsen L, Jenssen A J K. A Solution to the Airwave-Removal Problem in Shallow-Water Marine EM[J]. Geophysics, 2010, 75(5):37-42.

[17] Ziolkowski A, Wright D. Removal of the Airwave in Shallow Marine Transient EM Data[C]// 77th Annual International Meeting. San Antonio: Society of Exploration Geophysicists, 2007:534-538.

[18] Lu X, Srnka L J, Carazzone J J. Method for Removing Air Wave Effect from Offshore Frequency Domain Controlled-Source Electromagnetic Data: U S, O2005/010 560A1[P].2005-10-20.

[19] 殷长春，刘云鹤，翁爱华，等, 海洋可控源电磁法空气波研究现状及展望[J]. 吉林大学学报：地球科学版， 2012，42(5)，1506-1520. Yin Changchun, Liu Yunhe, Weng Aihua, et al. Research on Marine Controlled Source Electromagnetic Method Airwave[J]. Journal of Jilin University: Earth Science Edition, 2012, 42(5):1506-1520.

[20] 沈金松, 詹林森, 王鹏飞, 等.可控源海洋电磁勘探中空气波影响的理论分析和数值模拟[J]. 地球物理学报, 2012, 55(7):2473-2488. Shen Jinsong, Zhan Linsen, Wang Pengfei, et al. Theoretic Analysis and the Numerical Simulation of the Effects of the Air Wave Interaction on the Marine Controlled Source Electromagnetic Exploration[J]. Chinese Journal of Geophysics, 2012, 55(7): 2473-2488.

[21] 刘长胜, 林君. 海底表面磁源瞬变响应建模及海水影响分析[J]. 地球物理学报, 2006, 49(6): 1891-1898. Liu Changsheng, Lin Jun. Transient Electromagnetic Response Modeling of Magnetic Source on Seafloor and the Analysis of Seawater Effect[J]. Chinese Journal of Geophysics, 2006, 49(6):1891-1898.

[22] Nordskag J I, Amundsen L, Ursin B. Decomposition of Marine Electromagnetic Fields into TE and TM Modes for Enhanced Interpretation[C]//EGM International Workshop. [S. l.]: EGM, 2007.

[23] Ward S H, Hohmann G W. Electromagnetic Theory for Geophysical Applications[C]//Ward H S, Hohmann W G. Electromagnetic Methods in Applied Geophysics. [S.l.]:SEG, 1988:131-311.

[24] Weiss C J. The Fallacy of the Shallow-Water Problem in Marine CSEM Exploration[J]. Geophysics, 2007, 72(6): 93-97.

[25] Løseth L O, Pedersen H M, Schaug-Pettersen T, et al. A Scaled Experiment for the Verification of the Seabed Logging Method[J]. Journal of Applied Geophysics, 2007, 64(3/4):47-55.

[26] Andréis D, MacGregor L. Controlled-Source Electromagnetic Sounding in Shallow Water: Principles and Applications[J]. Geophysics, 2008, 73(1):21-32.

[27] Ziolkowski A. Developments in the Transient Electromagnetic Method[J]. First Break, 2007, 25(6): 99-106.

[28] Mittet R. Normalized Amplitude Ratios for Frequency-Domain CSEM in Very Shallow Water[J]. First Break, 2008, 26(11):47-54.

[29] Wapenaar K, Slob E, Snieder R. Seismic and Electromagnetic Controlled-Source Interferometry in Dissipative Media[J]. Geophysical Prospecting, 2008, 56(3): 419-434.

Methods of Air Wave Approximation and Decomposition for Marine Controlled Source Electromagnetic Measurements in Shallower Water Domain

Shen Jinsong1,2,3, Wang Xuan1, Wei Shuaishuai1, Li Man1

1.FacultyofGeophysicsandInformatics,ChinaPetroleumUniversity(Beijing)，Beijing102249，China2.StateKeyLaboratoryofPetroleumResourceandProspecting,Beijing102249,China
3.CNPCKeyLabofGeophysicalExploration,Beijing102249，China

In marine controlled-source electromagnetic (MCSEM) exploration, the influence on the EM response of the airwave has been paid great attention by the industry and the air wave is one kind of EM signal from the horizontal electric dipole (HED) source that is reflected from or travels along the air-sea interface and return back to receivers. In shallow water domains, the airwave dominates over the effective signal from seabed sediments, and makes it very difficult to identify the hydrocarbon reservoirs. On the basis of the layer model, we have derived airwave expression of half space resistivity model by using the mode decomposition theory of the EM field and the air wave in the sea water layer is approximately expressed as a series of wave reverberations that travel down and backward between the sea surface and the seabed. According to the propagations of EM field near sources and receivers, we have obtained airwave expression of a finite water layer. In order to compare the effectiveness of the various air wave mitigation methods, we have adopted the decomposition method of EM fields into upgoing and downgoing components, and have acuired larger anomaly of several times than that of the raw EM field above the resistive formation, and indicated the favorable effect of airwave attenuation by using of wave separation. Finally, simulation results of electromagnetic responses from different background and reservoir models have been used to check the effectiveness and adaptability of the two air wave mitigation methods of wave reverberations approximation and wave field decomposition, numerical results show that for model of horizontal layers, the latter is more effective, however, the former is applicable to rugged sea floor formation and the latter only fits for horizontal layers.

shallower water domain; marine controlled-source electromagnetic exploration; wave reverberations approximation; air wave separation; air wave mitigation

10.13278/j.cnki.jjuese.201501302.

2014-04-06

国家自然科学基金项目(41374141)；国家“973”计划项目(2013CB228605)

沈金松(1964——)，男，研究员，博士，主要从事地球物理测井及电磁探测理论方法和应用研究，E-mail:shenjinsong@cup.edu.cn。

10.13278/j.cnki.jjuese.201501302

P631.3

A

沈金松，汪轩，魏帅帅，等. 浅水域可控源海洋电磁测量中空气波的渐近表达与波场分离方法.吉林大学学报:地球科学版,2015,45(1):282-292.

Shen Jinsong, Wang Xuan, Wei Shuaishuai, et al. Methods of Air Wave Approximation and Decomposition for Marine Controlled Source Electromagnetic Measurements in Shallower Water Domain.Journal of Jilin University:Earth Science Edition,2015,45(1):282-292.doi:10.13278/j.cnki.jjuese.201501302.