李 盼 李予國(guó)*②(中國(guó)海洋大學(xué)海洋地球科學(xué)學(xué)院，山東青島 266100； ②海洋海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，山東青島 266100)

1 引言

中國(guó)廣闊的海域蘊(yùn)藏著大量的石油和天然氣等資源，地震勘探是探測(cè)海洋油氣資源的主要地球物理方法，利用地震數(shù)據(jù)可以獲得精確的地質(zhì)構(gòu)造信息并預(yù)測(cè)儲(chǔ)層孔隙度，卻很難判斷構(gòu)造中流體的性質(zhì)[1，2]，而海洋可控源電磁法(Marine Controlled Source Electromagnetic，MCSEM)是探測(cè)海底油氣、天然氣水合物以及礦產(chǎn)資源的一種海洋地球物理勘探新方法[3，4]。它通過(guò)探測(cè)含油氣儲(chǔ)層與周圍沉積層之間的電阻率差異引起的電磁異常，研究海底地層的電阻率分布，根據(jù)電阻率與儲(chǔ)層含油氣飽和度的關(guān)系，可推斷其含油氣性。利用該方法可以識(shí)別高阻油氣藏，提高鉆井成功率，降低勘探成本。因此，MCSEM法已成為各大石油勘探公司鉆前儲(chǔ)層評(píng)價(jià)的重要手段之一[5]。

常規(guī)海洋電磁勘探方法在深水海域用水平電偶極源向海底發(fā)射低頻(0.01～10Hz)電磁信號(hào)，由布設(shè)在海底的電磁接收儀器記錄來(lái)自海底地層的電場(chǎng)和磁場(chǎng)信號(hào)[6]。該方法對(duì)于深水區(qū)域海底高阻油氣藏具有很高的分辨率，主要是因?yàn)榫藓竦暮Ｋ畬涌梢杂行諝獠ǖ母蓴_[7]。因此，早期的海洋電磁觀測(cè)都集中在深水區(qū)域(海水深度大于1000m)，并取得了良好的效果[8]。然而，在淺水區(qū)(海水深度小于300m)的應(yīng)用遇到了一定困難。在淺水區(qū)使用常規(guī)的水平電偶極作為發(fā)射源進(jìn)行MCSEM勘探時(shí)，接收到的信號(hào)中空氣波占主導(dǎo)地位，它掩蓋了從海底介質(zhì)中返回的有效信號(hào)，特別是來(lái)自高阻層中的導(dǎo)波信號(hào)，這嚴(yán)重影響了反演解釋結(jié)果[9]。因此，在淺水區(qū)MCSEM資料處理中抑制空氣波的影響，對(duì)提高儲(chǔ)層識(shí)別能力具有重要意義。

合成孔徑技術(shù)在雷達(dá)和聲吶領(lǐng)域已經(jīng)得到廣泛應(yīng)用，其核心思想是通過(guò)不同場(chǎng)源的相互干涉構(gòu)建一個(gè)加長(zhǎng)的合成源，并使其能量向著特定方向加強(qiáng)和傳播。Fan等[10,11]已經(jīng)證明散射場(chǎng)可以通過(guò)適當(dāng)?shù)南辔恍Ｕ驼穹Ｕ晦D(zhuǎn)向和加強(qiáng)。汪軒等[12]利用合成孔徑與徑向消除方法較好地壓制了淺?？煽卦措姶趴碧街锌諝獠ǖ挠绊?。Knaak等[13，14]提出了2D合成孔徑源的概念，并將其應(yīng)用于三維海洋可控源電磁數(shù)據(jù)處理。Yoon等[15]和Zhdanov等[16]改進(jìn)了Fan等的選取校正系數(shù)的方法，通過(guò)參數(shù)優(yōu)化提高合成孔徑源對(duì)電磁異常識(shí)別效果。

目前常用于MCSEM資料解釋的電磁響應(yīng)異常有歸一化振幅、有效異常和可探測(cè)度等。歸一化振幅是最常用的異常響應(yīng)之一[17,18]，但計(jì)算歸一化振幅時(shí)，沒(méi)有考慮觀測(cè)誤差和噪聲干擾等因素的影響[19]。鑒于此，趙慧等[20]把觀測(cè)誤差和噪聲等干擾因素考慮在內(nèi)，提出了有效異常計(jì)算方法，用于確定探測(cè)海底天然氣水合物儲(chǔ)層的MCSEM勘探參數(shù)。在深水區(qū)，基于有效異常進(jìn)行MCSEM響應(yīng)分析非常有效，但在淺水區(qū)效果欠佳。有效異常最明顯的弱點(diǎn)是沒(méi)有包含相位信息，劉婷婷等[21]在考慮相位的基礎(chǔ)上定義了可探測(cè)度，并通過(guò)數(shù)值模擬證明其在淺水區(qū)可有效識(shí)別高阻薄層。

海洋電磁場(chǎng)數(shù)值模擬中，各向同性的地電模型可以很好地近似地下電性結(jié)構(gòu)[22]。因此，本文首先利用MCSEM一維正演程序，計(jì)算淺水區(qū)不同水深模型的電磁場(chǎng)響應(yīng)，比較了使用不同相位和振幅校正時(shí)合成源響應(yīng)的歸一化振幅、有效異常和可探測(cè)度對(duì)高阻層的識(shí)別能力；然后，將可探測(cè)度作為相位校正和振幅校正參數(shù)選擇的基準(zhǔn)，分別利用MCSEM二維和三維正演程序[23-25]計(jì)算二維和三維模型采用相對(duì)最佳組合校正參數(shù)時(shí)的合成源可探測(cè)度。結(jié)果表明利用可探測(cè)度選擇校正參數(shù)能夠有效地提高M(jìn)CSEM合成源干涉法對(duì)高阻體的識(shí)別能力。

2 MCSEM合成源干涉法

在MCSEM勘探中，利用不同發(fā)射源電磁信號(hào)的干涉構(gòu)成一個(gè)加長(zhǎng)的合成源，其表達(dá)式[26]為

(1)

式中： 某一位置r處角頻率為ω的合成源信號(hào)SA(r,ω)是由r1到rN處發(fā)射源的信號(hào)疊加而成；sn(r,rn,ω)表示第n個(gè)源在rn處產(chǎn)生的電磁場(chǎng)信號(hào)，并且每個(gè)源產(chǎn)生的電磁場(chǎng)都經(jīng)過(guò)振幅修正an和相位校正φn。

令an=exp(-c2αΔx)、φn=-c1αΔx，則合成源信號(hào)可以進(jìn)一步表示為[26]

sn(x,xn,ω)

(2)

圖1 MCSEM電磁信號(hào)傳播路徑示意圖

根據(jù)MCSEM信號(hào)主要傳播路徑，可以解釋合成孔徑方法突出高阻異常和抑制空氣波的原因。由圖1可知，電磁信號(hào)從發(fā)射源(Tx)到達(dá)接收站(Rx)主要有三種途徑(分別見圖1中的①、②、③)： (a)發(fā)射源電磁信號(hào)垂直穿過(guò)海水層到達(dá)?？战缑?，然后沿海平面滑行一段距離，再折射垂直向下穿過(guò)海水層到達(dá)接收站，這部分被接收到的電磁信號(hào)稱作空氣波； (b)電磁信號(hào)從發(fā)射源穿過(guò)海水層直接到達(dá)接收站； (c)電磁信號(hào)經(jīng)過(guò)海水層和海底地層后返回至接收站。其中，從海底地層傳播回來(lái)并被接收站接收到的電磁信號(hào)才是有用信號(hào)。可用于研究海底地質(zhì)構(gòu)造及圈閉含油氣情況。

當(dāng)對(duì)合成源的單源電磁信號(hào)進(jìn)行相位校正時(shí)電磁場(chǎng)將發(fā)生轉(zhuǎn)向，電磁信號(hào)傳播路徑會(huì)發(fā)生改變，因此通過(guò)選擇合適的相位校正系數(shù)可以使儲(chǔ)層響應(yīng)得到加強(qiáng)、空氣波得到抑制。其主要原因有兩個(gè)：首先，電磁場(chǎng)能量大部分是垂直于電偶極源向外傳播的，因此當(dāng)發(fā)射源轉(zhuǎn)向高阻體時(shí)，相比發(fā)射源水平時(shí)有更多的能量集中到高阻體，于是接收站接收到的有效信息則相對(duì)增加；其次，當(dāng)電磁場(chǎng)發(fā)生轉(zhuǎn)向時(shí)，豎直向上傳播的能量減少，在海面與空氣的分界面上轉(zhuǎn)換為空氣波的能量相應(yīng)減少，因此空氣波得到抑制。同時(shí)，通過(guò)適當(dāng)?shù)恼穹Ｕ梢哉{(diào)諧不同發(fā)射源電磁信號(hào)之間的相互干涉，使干涉增強(qiáng)，可以進(jìn)一步抑制空氣波，增強(qiáng)目標(biāo)體的異常響應(yīng)。

3 一維模型合成源電磁響應(yīng)異常分析

3.1 電磁響應(yīng)異常

目前用于MCSEM資料解釋的電磁響應(yīng)異常有歸一化振幅、有效異常和可探測(cè)度等。歸一化振幅是指含高阻層模型與背景模型的電磁場(chǎng)振幅之比

(3)

式中：|Eres|指含高阻層模型的電磁場(chǎng)振幅； |Eback|指背景模型的電磁場(chǎng)振幅。對(duì)于深水模型，若R大于1，則表明海底下方可能存在高阻層。但當(dāng)海水深度變淺時(shí)，由于空氣波影響增強(qiáng)，常規(guī)的歸一化振幅對(duì)地下高阻層的識(shí)別能力變?nèi)?。另外，在?jì)算歸一化振幅時(shí)，由于沒(méi)有考慮觀測(cè)誤差和海底電磁數(shù)據(jù)記錄儀及其測(cè)量系統(tǒng)誤差等因素的影響，因此依據(jù)歸一化振幅解釋MCSEM資料時(shí)可能會(huì)產(chǎn)生較大偏差。

裴建新等[19]綜合考慮了海底電磁采集站方位觀測(cè)誤差、數(shù)據(jù)記錄儀和測(cè)量系統(tǒng)誤差的影響，定義了有效異常

(4)

其中噪聲表達(dá)式為[19]

(5)

式中：Enoise為電磁場(chǎng)的噪聲；Erel為電磁場(chǎng)振幅的相對(duì)誤差，通常假定其為背景場(chǎng)振幅值的5%；Erol為由海底電磁采集站電場(chǎng)測(cè)量臂或磁場(chǎng)傳感器方向的不確定性引起的誤差；Eabs為絕對(duì)誤差，即海底電磁站數(shù)據(jù)記錄儀的本底噪聲，通常設(shè)為1.0×10-16V·A-1·m-2。

由于電磁場(chǎng)為復(fù)數(shù)場(chǎng)，它既包含振幅信息也包含相位信息，而歸一化振幅和有效異常都沒(méi)有考慮相位信息。因此，在計(jì)算有效異常的基礎(chǔ)上，為了進(jìn)一步提高對(duì)高阻油氣層的識(shí)別能力，劉婷婷等[21]將相位信息考慮在內(nèi)，對(duì)有效異常公式進(jìn)行了修正，定義了可探測(cè)度

(6)

D反映了對(duì)高阻油氣層的可探測(cè)程度。D值越大，表明高阻油氣層越容易被探測(cè)到。劉婷婷等[21]計(jì)算了不同水深情況下一維水合物模型的MCSEM歸一化振幅、有效異常和可探測(cè)度。結(jié)果表明，在深水環(huán)境下利用歸一化振幅、有效異常和可探測(cè)度均可識(shí)別水合物高阻薄層；但在淺水環(huán)境下，歸一化振幅和有效異常受空氣波影響畸變嚴(yán)重，不能很好地識(shí)別高阻薄層，而空氣波對(duì)可探測(cè)度影響較小，在淺水區(qū)利用可探測(cè)度有利于識(shí)別水合物高阻薄層。

3.2 一維模型合成源響應(yīng)分析

設(shè)計(jì)如圖2所示的海洋一維層狀模型，海水層厚hsw是可變的，分別設(shè)為200、100、50m。海底之下第二層電阻率若等于100Ω·m，則代表高阻模型，若電阻率為1Ω·m，則代表背景模型?；诪榫鶆虬肟臻g。

圖2 變海水層厚一維層狀模型

采用Inline模式(軸向)觀測(cè)系統(tǒng)，51個(gè)水平電偶極源(Tx)均位于海底上方20m處，間距為100m，分布范圍為-2.5～2.5km，發(fā)射頻率為0.25Hz； 201個(gè)接收站(Rx)位于海底，間距為200m，分布范圍為-20～20km。分別計(jì)算海水層厚為50、100和200m三種情況下合成源響應(yīng)的歸一化振幅、有效異常和可探測(cè)度，并對(duì)這些參數(shù)進(jìn)行對(duì)比分析。

當(dāng)對(duì)合成源進(jìn)行相位校正(c2=0)時(shí)，歸一化振幅最大值(Rmax)、有效異常最大值(Smax)、可探測(cè)度最大值(Dmax)隨相位校正系數(shù)c1的變化曲線分別如圖3上所示。由圖可見： ①對(duì)于同一相位校正系數(shù)，三個(gè)電磁響應(yīng)異常最大值都隨著海水層厚減小而降低； ②三個(gè)參數(shù)的最大值都呈現(xiàn)減小—增大—減小的規(guī)律； ③當(dāng)相位校正系數(shù)c1約為0.7時(shí)，校正效果達(dá)到最佳，當(dāng)?shù)扔?.3時(shí)，效果最差；④通過(guò)選擇合適的c1值，高阻層的異常響應(yīng)可以提高10%～20%。當(dāng)對(duì)合成源進(jìn)行振幅校正(即c1=0)時(shí)，Rmax、Smax、Dmax隨c2的變化曲線如圖3下所示。由圖可見： ①對(duì)于同一振幅校正系數(shù)，三個(gè)電磁響應(yīng)異常最大值都隨著海水層厚減小而降低； ②三個(gè)電磁響應(yīng)異常最大值曲線都呈現(xiàn)減小—增大—減小的趨勢(shì)； ③當(dāng)c2約為0.5時(shí)，校正效果最差，而當(dāng)c2=0時(shí)，校正效果最佳。

為了對(duì)比上述三個(gè)參數(shù)在同一條件下的特征差異，對(duì)圖3重新組合，得到圖4。由圖4可見，在三種不同海水層厚下，無(wú)論是進(jìn)行合成源相位校正還是振幅校正，可探測(cè)度參數(shù)的效果最佳，有效異常次之，歸一化振幅最差。

圖3 海水層厚分別為50、100和200m三種情況下合成源產(chǎn)生的Rmax(a)、Smax(b)、Dmax(c)隨c1(上)和c2(下)變化曲線

圖4 海水層厚分別為50(a)、100(b)和200m(c)三種情況下合成源產(chǎn)生的Rmax、Smax、Dmax隨c1(上)和c2(下)變化曲線

圖5 水層厚度分別為50m(上)、100m(中)、200m(下)時(shí)不同校正參數(shù)(c1、c2)組合時(shí)合成源電場(chǎng)Ey分量歸一化振幅最大值Rmax(a)、有效異常最大值Smax(b)、可探測(cè)度最大值Dmax(c)分布圖

綜上可知，在對(duì)不同海水層厚模型合成源進(jìn)行相位校正或振幅校正時(shí)，D相比于R和S對(duì)高阻層具有更高的識(shí)別能力。

下面分析對(duì)合成源進(jìn)行相位與振幅綜合校正時(shí)上述三個(gè)參數(shù)的特征。

圖5為三種海水層厚情形下不同相位校正與振幅校正參數(shù)組合時(shí)得到的合成源電場(chǎng)Ey分量上述三個(gè)參數(shù)最大值的分布圖，激發(fā)頻率為0.25Hz。由圖可知，合成源歸一化振幅最大值Rmax相對(duì)較高的區(qū)域分布范圍較小且區(qū)域不連續(xù)；合成源有效異常最大值Smax和可探測(cè)度最大值Dmax的高值分布范圍較大且相對(duì)連續(xù)，而且隨著海水深度加大，其分布范圍也隨之增大?？梢?，采用振幅校正系數(shù)c1和相位校正系數(shù)c2綜合校正時(shí)，Dmax分布范圍最大、連續(xù)性最好，Smax次之，Rmax最差。

4 二維模型合成源電磁響應(yīng)異常

設(shè)計(jì)的二維模型如圖6所示，海水電阻率為0.3Ω·m，海水層厚為50m，二維體埋深為1000m，厚度為100m，其在y方向的延伸范圍為-2～2km，即寬度為4km。若二維體電阻率為100Ω·m，代表高阻體模型， 若電阻率為1Ω·m，則代表背景模型。

圖6 二維海洋地電模型

采用Inline采集模式，假定51個(gè)水平電偶極源(Tx)均位于海底上方20m處，分布范圍為-9～-4km，間距為100m，發(fā)射頻率為0.25Hz，126個(gè)接收站(Rx)位于海底，分布范圍為-10～15km，間距為200m。分別分析單源、簡(jiǎn)單合成源、相位校正合成源以及相位與振幅組合校正合成源的電磁響應(yīng)異常。

圖7為前述四種場(chǎng)源條件下異常體模型與背景模型的電場(chǎng)Ey分量振幅曲線。由圖7a可見，單個(gè)電偶源條件下，高阻模型與背景模型的電場(chǎng)振幅曲線基本重合，這主要是因?yàn)楹Ｋ畬虞^小，空氣波占據(jù)電磁信號(hào)的主導(dǎo)地位，掩蓋了高阻體異常響應(yīng)。圖7b為51個(gè)簡(jiǎn)單合成源(即c1=0，c2=0)產(chǎn)生的水平電場(chǎng)振幅曲線，圖7c為經(jīng)過(guò)相位校正后(即c2=0)的合成源歸一化振幅曲線，這里參照前述一維模型分析結(jié)果，選擇c1=0.7。由圖7b和圖7c可知，采用簡(jiǎn)單合成源和相位校正合成源時(shí)，高阻模型和背景模型的電場(chǎng)歸一化振幅相比單源情形有所增強(qiáng)，但是這兩種模型的歸一化振幅曲線基本重合在一起，并未見明顯的分離；圖7d為相位和振幅組合校正后合成源的電場(chǎng)歸一化振幅曲線，可見高阻體模型與背景模型的電場(chǎng)歸一化振幅曲線發(fā)生了分離。

圖7 四種場(chǎng)源條件下高阻模型和背景模型的電場(chǎng)Ey分量振幅曲線(a)單源； (b)簡(jiǎn)單合成源； (c)相位校正合成源；(d)相位與振幅綜合校正合成源

為了進(jìn)一步分析合成源的效果，基于圖7繪制了四種場(chǎng)源(單源、簡(jiǎn)單合成源、相位校正合成源及相位與振幅綜合校正合成源)的歸一化振幅、有效異常和可探測(cè)度曲線(圖8)。

圖8 二維模型單源、簡(jiǎn)單合成源、相位校正合成源(c1=0.7)、相位與振幅綜合校正合成源(c1=0.7，c2=0.2)的歸一化振幅R(a)、有效異常S(b)和可探測(cè)度D(c)曲線

由圖8a可知：?jiǎn)卧础⒑?jiǎn)單合成源和相位校正合成源(c1=0.7)三種場(chǎng)源的歸一化振幅曲線基本重合，因此這三種源對(duì)高阻體的識(shí)別能力較弱；相位與振幅組合校正合成源的歸一化曲線存在一個(gè)明顯的負(fù)異常，當(dāng)收發(fā)距大于-1km時(shí)，歸一化振幅隨著收發(fā)距的增大而減小，當(dāng)收發(fā)距為1.2km時(shí)，歸一化振幅值達(dá)到極小值0.4，然后歸一化值開始隨著收發(fā)距的繼續(xù)增大而增大，當(dāng)收發(fā)距為5km時(shí)歸一化振幅基本趨于穩(wěn)定。由于歸一化振幅值小于1，基本無(wú)法識(shí)別高阻體，因此在淺海區(qū)利用歸一化振幅識(shí)別高阻體具有局限性。

由圖8b可知，單源有效異常曲線隨著收發(fā)距的增大未發(fā)生較明顯變化，其最大值為0.31。然而，簡(jiǎn)單合成源、相位校正合成源及相位與振幅組合校正合成源三種場(chǎng)源的有效異常曲線均存在一個(gè)明顯的正異常，其最大值分別為2.45、4.20和11.4。由此可知，在淺水區(qū)，利用有效異?？勺R(shí)別高阻體，且經(jīng)過(guò)相位與振幅組合校正后的合成源有效異常對(duì)高阻體的識(shí)別能力最強(qiáng)，相位校正合成源的效果次之，單源效果最差。

由圖8c可知，簡(jiǎn)單合成源、相位校正合成源及相位與振幅組合校正合成源的可探測(cè)度曲線均存在一個(gè)正異常，可探測(cè)度值分別為2.5、5.56和12.08。由此可見，在淺海區(qū)經(jīng)過(guò)相位與振幅組合校正后的合成源的可探測(cè)度參數(shù)對(duì)高阻體識(shí)別能力最強(qiáng)，相位校正源效果次之，單源效果最不理想。

綜上可知，無(wú)論單源、簡(jiǎn)單合成源、相位校正合成源還是相位與振幅組合校正合成源，其可探測(cè)度對(duì)于識(shí)別高阻體具有最好的效果，有效異常效果次之，歸一化振幅效果最差。

5 三維模型合成源的可探測(cè)度

采用如圖9所示三維模型。若三維體電阻率為100Ω·m，則代表高阻體模型；若電阻率為1Ω·m，則為背景模型。采用Inline采集模式，51個(gè)水平電偶極源(Tx)均沿y方向置于海底上方20m處，范圍為-9～-4km，間距為100m，發(fā)射頻率為0.25Hz。126個(gè)接收站(Rx)也沿y方向布置于海底-10～15km的范圍內(nèi)，間距為200m。

圖9 三維海洋模型

分別計(jì)算單源、簡(jiǎn)單合成源及相位與振幅綜合校正合成源的三維模型的電場(chǎng)振幅(圖10)。圖10中兩條垂直虛線表示高阻體在y方向的延伸范圍。由圖10d可見，經(jīng)相位與振幅綜合校正后的合成源可探測(cè)度得到明顯提高，相比單源情形可探測(cè)度異常值提高了約20倍，且異常覆蓋的范圍大部分位于三維體正上方，為推斷異常體橫向分布范圍提供了依據(jù)?？偟膩?lái)說(shuō)，經(jīng)過(guò)綜合校正后的合成源的可探測(cè)度利于識(shí)別三維異常體。

圖10 三維模型單源(a)、簡(jiǎn)單合成源(b)及相位與振幅綜合校正合成源(c)產(chǎn)生的振幅曲線以及三種源場(chǎng)源的可探測(cè)度曲線(d)

6 實(shí)測(cè)數(shù)據(jù)合成源電磁響應(yīng)

將上述合成源方法應(yīng)用于某測(cè)區(qū)實(shí)測(cè)MCSEM數(shù)據(jù)處理。測(cè)區(qū)海底深度約為770m，高阻儲(chǔ)層埋深約為1500m。測(cè)區(qū)的背景模型根據(jù)地質(zhì)和地球物理資料構(gòu)建。

圖11為實(shí)測(cè)MCSEM數(shù)據(jù)電磁響應(yīng)剖面。由圖11a可知，當(dāng)收發(fā)距較小時(shí)，背景模型的理論電場(chǎng)振幅曲線與實(shí)測(cè)曲線基本重合；隨著收發(fā)距增大，二者逐漸分離，高阻體的影響逐漸凸顯；但當(dāng)收發(fā)距大于5km后，實(shí)測(cè)曲線開始出現(xiàn)振蕩，且隨著收發(fā)距增大，振蕩現(xiàn)象更嚴(yán)重，這是數(shù)據(jù)質(zhì)量變差所致。

構(gòu)建一個(gè)4.1km長(zhǎng)的合成源，相位與振幅綜合校正后的合成源的振幅曲線如圖11b所示。由圖11b可知，當(dāng)收發(fā)距大于5km后，隨著收發(fā)距的增大，實(shí)測(cè)數(shù)據(jù)與背景場(chǎng)的電場(chǎng)振幅曲線逐漸發(fā)生分離，分離程度先增大后減小，即電磁異常清晰可見；同時(shí)，相比大收發(fā)距時(shí)實(shí)測(cè)單源電場(chǎng)振幅曲線的劇烈振蕩，實(shí)測(cè)數(shù)據(jù)的合成源振幅曲線更光滑。

圖11c為可探測(cè)度曲線，由圖可知，在5～10km范圍內(nèi)，單源的可探測(cè)度曲線振蕩較劇烈，而合成源的可探測(cè)度曲線在該區(qū)間平緩變化，即先緩慢增大再逐漸減小，充分說(shuō)明了合成源具有較好的效果。

圖11 實(shí)測(cè)MCSEM數(shù)據(jù)電磁響應(yīng)(a)單源； (b)相位與振幅綜合校正后合成源； (c)可探測(cè)度曲線

7 結(jié)論

本文通過(guò)計(jì)算海洋一維、二維模型正演電磁響應(yīng)，探討了歸一化振幅、有效異常和可探測(cè)度在淺水區(qū)對(duì)高阻層的識(shí)別能力，結(jié)果表明可探測(cè)度參數(shù)的識(shí)別能力最強(qiáng)、有效異常次之、歸一化振幅最差。在此基礎(chǔ)上，建立海洋三維模型，計(jì)算不同源情形下的可探測(cè)度，結(jié)果表明，相位與振幅綜合校正合成源的可探測(cè)度對(duì)異常指示效果最好，僅相位校正合成源的效果次之，單源識(shí)別能力最差。將該方法應(yīng)用于實(shí)際數(shù)據(jù)處理，經(jīng)過(guò)合成源處理后可探測(cè)度能在一定程度上指示高阻異常體的存在。

參考文獻(xiàn)

[1] 嚴(yán)波，李予國(guó)，韓波等.任意方位電偶源的MCSEM電磁場(chǎng)三維正演.石油地球物理勘探，2017，52(4)：859-868.

Yan Bo，Li Yuguo，Han Bo et al.3D forward modeling of MCSEM electromagnetic field of arbitrary azimuth dipole source.OGP，2017，52(4)：859-868.

[2] 徐凱軍，杜潤(rùn)林，劉展.海洋可控源電磁與地震一維聯(lián)合儲(chǔ)層參數(shù)反演.石油地球物理勘探，2016，51(1)：197-203.

Xu Kaijun，Du Runlin，Liu Zhan.Inversion of 1D re-servoir parameters of MCSEM and seismic.OGP，2016，51(1)：197-203.

[3] 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.First Break，2002，20(2)：144-152.

[4] Constable S.Ten years of marine CSEM for hydrocarbon exploration.Geophysics，2010，75(5)：75A67-75A81.

[5] Constable S，Srnka L J.An introduction to marine controlled-source electromagnetic methods for hydrocarbon exploration.Geophysics，2007，72(2)：WA3-WA12.

[6] 何展翔，孫衛(wèi)斌，孔繁恕等.海洋電磁法.石油地球物理勘探，2006，41(4)：451-457.

He Zhanxiang，Sun Weibin，Kong Fanshu et al.Marine electromagnetic approach.OGP，2006，41(4)：451-457.

[7] 周建美，李貅，戚志鵬.淺水域各向異性地層中的瞬變電磁響應(yīng)分析.石油地球物理勘探，2016，51(4)：821-830.

Zhou Jianmei，Li Xiu，Qi Zhipeng.Transient electromagnetic response analysis of anisotropic formation in shallow water.OGP，2016，51(4)：821-830.

[8] Andréis D，MacGregor L.Controlled-source electro-magnetic sounding in shallow water：Principles and applications.Geophysics，2008，73(1)：F21-F32.

[9] 殷長(zhǎng)春，劉云鶴，翁愛(ài)華等.海洋可控源電磁法空氣波研究現(xiàn)狀及展望.吉林大學(xué)學(xué)報(bào)(地球科學(xué)版)，2012，42(5)：1506-1520.

Yin Changchun，Liu Yunhe，Weng Aihua et al.Current status and prospect of marine controlled source electromagnetic airwave.Journal of Jilin University (Earth Science Edition)，2012，42(5)：1506-1520.

[10] Fan Y，Snieder R，Slob E et al.Steering and focusing diffusive fields using synthetic aperture.Epl，2011，95(3)：34006.

[11] Fan Y，Snieder R，Slob E et al.Synthetic aperture controlled source electromagnetics.Geophysical Research Letters，2010，37(13)：69-73.

[12] 汪軒，沈金松.綜合合成孔徑和徑向消除壓制淺?？煽卦措姶趴碧街锌諝獠ㄓ绊?中國(guó)地球科學(xué)聯(lián)合學(xué)術(shù)年會(huì)，2014.

Wang Xuan，Shen Jinsong.Apply synthetic aperture and radial elimination method to depress the airwave of controlled source electromagnetic exploration in shallow sea.Joint Academic Annual Meeting of Chinese Geoscience，2014.

[13] Knaak A，Snieder R，F(xiàn)an Y et al.3D synthetic aperture and steering for controlled-source electromagne-tics.The Leading Edge，2013，32(8)：972-978.

[14] Knaak A.3D synthetic aperture for controlled-source electromagnetics.Geophysics，2015，80(6)：E309-E316.

[15] Yoon D，Zhdanov M S.An optimal synthetic aperture method for the creation of directional sensitivity and removal of the airwave effect in MCSEM data.SEG Technical Program Expanded Abstracts，2014，33：685-690.

[16] Zhdanov M S，Yoon D，Mattsson J.Rapid imaging of towed streamer EM data using the optimal synthetic aperture method.IEEE Geoscience & Remote Sensing Letters，2017，99：1-5.

[17] Constable S，Weiss C J.Mapping thin resistors and hydrocarbons with marine EM methods：Insights from 1D modeling.Geophysics，2006，71(2)：G43-G51.

[18] Mittet R.Normalized amplitude ratios for frequency-domain CSEM in very shallow water.Aesthetic Plastic Surgery，2014，26(11)：47-54.

[19] 裴建新，王啟，張秀麗.海洋CSEM探測(cè)海底天然氣水合物的有效異常研究.石油地球物理勘探，2015，50(1)：177-183.

Pei Jianxin，Wang Qi，Zhang Xiuli.Effective anomaly study of gas hydrate detection using marine CSEM method.OGP，2015，50(1)：177-183.

[20] 趙慧，李予國(guó)，劉婷婷.利用可探測(cè)度探討海洋可控源電磁響應(yīng)的等值性.地球物理學(xué)進(jìn)展，2016，31(2)：928-936.

Zhao Hui，Li Yuguo，Liu Tingting.Study on the equivalence in marine controlled source electromagnetic method using detectivity.Progress in Geophysics，2016，31(2)：928-936.

[21] 劉婷婷，李予國(guó).海洋可控源電磁法對(duì)天然氣水合物高阻薄層的可探測(cè)度.海洋地質(zhì)前沿，2015，31(6)：17-22.

Liu Tingting，Li Yuguo.Detectivity of gas hydrate high-resistivity layer with marine CSEM method.Marine Geology Frontiers，2015，31(6)：17-22.

[22] 劉穎，李予國(guó).層狀各向異性介質(zhì)中任意取向電偶源的海洋電磁響應(yīng).石油地球物理勘探，2015，50(4)：755-765.

Liu Yin，Li Yuguo.Marine electromagnetic response of an arbitrary orientation galvanic source in layered anisotropic medium.OGP，2015，50(4)：755-765.

[23] Li Y G，Li G.Electromagnetic field expressions in the wavenumber domain from both the horizontal and vertical electric dipoles.Journal of Geophysics and Engineering，2016，13(4)：505-515.

[24] 李剛，李予國(guó)，韓波等.基于改進(jìn)的接收點(diǎn)插值算法的頻率域海洋可控源電磁法2.5維正演.地球物理學(xué)報(bào)，2017，60(12)：4887-4900.

Li Gang，Li Yuguo，Han Bo et al.2.5D marine CSEM modeling in the frequency-domain based on an improved interpolation scheme at receiver positions.Chinese Journal of Geophysics，2017，60(12)：4887-4900.

[25] Li G，Li Y G，Han B.Accurate interpolation at recei-ver positions：a novel method for frequency-domain marine CSEM finite-difference modelling.Pure and Applied Geophysics，2017，174(5)：2143-2160.

[26] Fan Y，Snieder R，Slob E et al.Increasing the sensitivity of controlled-source electromagnetics with synthetic aperture.Geophysics，2012，77(2)：E135-E145.