淺水域可控源海洋電磁測量中空氣波的漸近表達與波場分離方法

沈金松, 汪 軒，魏帥帥，李 曼

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.