楊 銳，王 猛，王功祥，鄧 明，景建恩，陳 凱

(1.中國(guó)地質(zhì)大學(xué)(北京) 地球物理與信息技術(shù)學(xué)院，北京 100083;2.中國(guó)地質(zhì)調(diào)查局 廣州海洋地質(zhì)調(diào)查局,廣州 510075 )

0 引言

近十年在國(guó)際海洋可控源電磁法相對(duì)成熟的背景下[1-3]，國(guó)內(nèi)海洋電磁法發(fā)展迅速[4-9]。作為電磁法的一個(gè)應(yīng)用分支，海洋可控源電磁法在海洋油氣藏邊界圈定和天然氣水合物勘探上扮演了越來(lái)越重要的角色。通過(guò)拖曳式或坐底式電磁發(fā)射機(jī)發(fā)射電磁場(chǎng)，海底多分量電磁接收機(jī)接收電磁場(chǎng)響應(yīng)，探測(cè)海底地層電性差異導(dǎo)致的電磁場(chǎng)異常分布特征，反演計(jì)算地下電阻率信息推算地質(zhì)結(jié)構(gòu)[11-13]。然而，電磁法的反演計(jì)算過(guò)程十分復(fù)雜且非常耗時(shí)，一次二維斷面的計(jì)算通常會(huì)花費(fèi)幾天甚至更長(zhǎng)的時(shí)間，因此在海上無(wú)法快速得知地下電性結(jié)構(gòu)信息。

David Myer[14]在一次試驗(yàn)中采用了歸一化幅值直接計(jì)算擬斷面圖進(jìn)行定性解釋；Swidinsky等[15]在研究過(guò)程中也有涉及過(guò)快速定性成像的方法，但這兩種方法對(duì)數(shù)據(jù)預(yù)處理過(guò)程要求較高，前期處理相對(duì)復(fù)雜耗時(shí);劉穎等[10]提出過(guò)一種海洋可控源電磁法視電阻率的定義和計(jì)算方法，該方法以空氣-海水-大地三層模型定義視電阻率計(jì)算公式，提高了定性解釋效率。筆者提出了一種以海水-大地半空間模型定義視電阻率的方法進(jìn)一步簡(jiǎn)化計(jì)算公式，并集成了一套將海試數(shù)據(jù)直接快速成像的算法。

1 快速成像計(jì)算方法

1.1 視電阻率定義與計(jì)算公式

為避開(kāi)反演過(guò)程中網(wǎng)格化后的大量運(yùn)算，直接對(duì)海底結(jié)構(gòu)做簡(jiǎn)化處理，類似于陸地電磁法對(duì)視電阻率的處理方式，將海水和地層看作兩個(gè)均勻各向同性的半空間，從麥克斯韋方程組推導(dǎo)地層電導(dǎo)率與海底電磁場(chǎng)響應(yīng)的關(guān)系式[16]如下：

(1)

其中：Eρ為海底電磁場(chǎng)響應(yīng)；Ids為電偶極矩；σ0為海水電導(dǎo)率；σ為地層電導(dǎo)率；ρ為收發(fā)距；φ為電偶極子與ρ夾角；γ2=iμω(σ+iεω)。通過(guò)公式(1)計(jì)算出來(lái)的等效地層電阻率即為視電阻率。由于電阻率與電磁場(chǎng)響應(yīng)并非線性關(guān)系，不能直接計(jì)算，所以采用迭代法或二分法求解。

海洋可控源電磁法拖曳式發(fā)射或者坐底式發(fā)射方式均是混頻時(shí)間域序列，在數(shù)據(jù)預(yù)處理階段，發(fā)射電流與電磁場(chǎng)值通過(guò)標(biāo)定、FFT、去噪及歸一化等步驟得到單頻時(shí)間域數(shù)據(jù)，再合并航道GPS數(shù)據(jù)推算收發(fā)距。令函數(shù)：

(2)

設(shè)定誤差限，代入?yún)?shù)，輸入初始電阻率范圍，二分迭代計(jì)算得到某一頻率下對(duì)應(yīng)收發(fā)距和點(diǎn)位的視電阻率值。

1.2 視電阻率快速成像算法

海洋可控源電磁法通常會(huì)經(jīng)過(guò)數(shù)據(jù)預(yù)處理得到MVO和PVO曲線，分別反應(yīng)電磁場(chǎng)振幅和相位信息，但這些信息并不能直觀描述電性分布。通過(guò)視電阻率與海底電磁場(chǎng)響應(yīng)的對(duì)應(yīng)關(guān)系，代入?yún)?shù)和電磁場(chǎng)振幅得到視電阻率，其對(duì)應(yīng)坐標(biāo)位置如圖1所示為發(fā)射點(diǎn)位和接收點(diǎn)位的中點(diǎn)，視深度為收發(fā)距的一半。

當(dāng)一條測(cè)線完成，一個(gè)斷面的數(shù)據(jù)就能覆蓋大部分地下區(qū)域，呈現(xiàn)倒三角或者倒梯形的數(shù)據(jù)集，視電阻率根據(jù)其幾何測(cè)深的位置會(huì)反映出電性差異，根據(jù)實(shí)際需求，選擇目標(biāo)區(qū)域網(wǎng)格化插值，作視電阻率擬斷面圖即為快速成像最終結(jié)果。

圖1 視電阻率擬斷面計(jì)算原理圖

圖2 一維高阻模型

圖3 一維模型MVO曲線

2 模型驗(yàn)證

海洋可控源電磁法常用于淺層油氣和天然氣水合物勘探，以彌補(bǔ)海洋地震方法對(duì)非層狀介質(zhì)探測(cè)物性特征判斷不準(zhǔn)的問(wèn)題。為驗(yàn)證視電阻率擬斷面圖快速成像算法的可靠性，以一維和二維的天然氣水合物高阻模型為例模擬計(jì)算。

圖4 一維模型快速成像

圖5 二維模型

圖6 二維模型MVO曲線

圖7 二維模型快速成像

圖2為一維高阻模型示意圖，海底位于z軸原點(diǎn)，電磁接收站位為x軸方向2 500 m～4 000 m，間隔150 m，發(fā)射點(diǎn)位為沿x軸方向-5 000 m～5 000 m，間隔50 m。海水層電阻率為0.3 Ω·m，海底圍巖電阻率為1 Ω·m，高阻薄層電阻率為100 Ω·m，厚度為100 m，位于海底500 m處。圖3是借助Kerry Key[17]的一維Occam正演程序計(jì)算的歸一化場(chǎng)值隨收發(fā)距變化的曲線(MVO曲線)，圖4展示了一維模型的快速成像(頻率為0.5 Hz)。從圖4可以看出，400 m～500 m處有明顯的高阻異常，分層效果很好，說(shuō)明該算法對(duì)一維高阻薄層有較好的識(shí)別作用。

二維地層比一維層狀結(jié)構(gòu)更貼合實(shí)際地質(zhì)情況(圖5)，本次二維模型設(shè)定了兩個(gè)高阻異常體，分別為沿x軸方向的高阻體A(-3 000 m至-2 000 m，電阻率為50 Ω·m，深度為100 m，厚度為50 m)與高阻體B(0 m至3 000 m，電阻率為10 Ω·m，深度為100 m，厚度為50 m)。海底依然設(shè)為z軸原點(diǎn)，接收機(jī)站位為x軸-3.5 km～3.5 km，點(diǎn)距為500 m，發(fā)射機(jī)點(diǎn)位為x軸方向-5 km～5 km，點(diǎn)距為50 m。海水電阻率為0.3 Ω·m，海底圍巖電阻率為1 Ω·m。圖6為當(dāng)前模型的MVO曲線，可以看出在高阻異常體的對(duì)應(yīng)位置曲線有明顯變化，但并不直觀反映電阻率大小，圖7為二維模型視電阻率快速成像結(jié)果(頻率為2 Hz)，如圖7所示，兩個(gè)高阻體均出現(xiàn)了明顯的視電阻率異?，F(xiàn)象，且橫向分辨率高，橫向邊界與模型相近，但是縱向邊界分辨能力比較差，視深度大約在500 m至1 000 m之間，與模型的高阻體實(shí)際深度100 m相差較大。

綜合圖4和圖7兩個(gè)斷面圖的結(jié)果，可以得出結(jié)論，視電阻率對(duì)電性異常有較明顯地反應(yīng)，在成像后能劃分出高阻體的近似區(qū)域。尤其是橫向分辨率高，圈定高阻體橫向邊界較準(zhǔn)確，而視深度則與真實(shí)深度有一定差距，無(wú)法確定與實(shí)際高阻體埋深的關(guān)系，故該算法可以作為一種簡(jiǎn)單快速且直觀的海洋可控源電磁法資料定性解釋方法。

3 海試數(shù)據(jù)快速成像與傳統(tǒng)反演對(duì)比

圖8 工區(qū)與測(cè)線點(diǎn)位分布圖

圖9 拖曳式發(fā)射方式快速成像

圖10 坐底式發(fā)射方式快速成像

圖11 反演結(jié)果對(duì)比

2015年3月，在國(guó)家重點(diǎn)研發(fā)計(jì)劃的支持下，搭載自主研發(fā)的海洋可控源電磁法發(fā)射機(jī)與接收機(jī)進(jìn)行了天然氣水合物重點(diǎn)靶區(qū)的探測(cè)試驗(yàn)。工區(qū)位于瓊東南海域，前期開(kāi)展過(guò)海洋地震勘探作業(yè)，水深在1 300 m～1 400 m之間，四個(gè)接收機(jī)站位沿北西方向布設(shè)，間隔約500 m，測(cè)線長(zhǎng)約20 km，與海洋地震測(cè)線重合，共進(jìn)行了往返兩次拖曳發(fā)射以及坐底式發(fā)射(6個(gè)發(fā)射站位，間距約200 m)，工區(qū)與測(cè)線點(diǎn)位分布情況如圖8所示。發(fā)射信號(hào)覆蓋0.5 Hz、1.5 Hz、8 Hz，拖曳發(fā)射電流約為260 App，坐底式發(fā)射電流約為100 App。

圖9和圖10分別是拖曳式發(fā)射與坐底式發(fā)射電磁場(chǎng)視電阻率擬斷面圖(頻率均為1.5 Hz)，R1、R2、R3和R4為由南至北的四個(gè)接收機(jī)站位，R1為坐標(biāo)原點(diǎn)，可以看出四個(gè)站位間隔約為500 m。在視深度約為200 m～300 m的地方，兩種方式等效視電阻率異常幾乎完全重合，均表現(xiàn)出了相同形態(tài)的高阻異常，且在R2和R3下方的高阻異常較R1和R4略淺。另外值得注意的是，坐底式數(shù)據(jù)異常表現(xiàn)更明顯，可能跟發(fā)射源更靠近海底致發(fā)射信號(hào)信噪比更高有關(guān)。運(yùn)算時(shí)間上，拖曳式發(fā)射方式計(jì)算單個(gè)擬斷面圖平均耗時(shí)87.73 s，而計(jì)算單個(gè)坐底式發(fā)射方式擬斷面圖僅需要16.70 s。

圖11為該測(cè)線上同一位置傳統(tǒng)處理流程的地震層析成像和電磁法Occam反演[9]結(jié)果，前者顯示了較為明顯的上下強(qiáng)反射界面，R1、R2與R3下方約100 m處有強(qiáng)反射，R3與R4之間下方約200 m處也出現(xiàn)了反射界面，后者電磁法反演結(jié)果表明，R2與R3之間下方100 m～300 m出現(xiàn)高阻異常，而兩邊異常相對(duì)較深。兩種反演資料處理結(jié)果與視電阻率異常形態(tài)吻合度高，故海洋可控源電磁法擬斷面圖快速成像在保證計(jì)算時(shí)間大大減少的情況下，能夠?qū)Φ叵码娦援惓Ｓ斜容^明顯且準(zhǔn)確的識(shí)別能力。

4 結(jié)論

通過(guò)對(duì)海洋可控源電磁法視電阻率定義以及擬斷面圖模型和海試數(shù)據(jù)驗(yàn)證，得出以下幾點(diǎn)結(jié)論：

1)2015年海試數(shù)據(jù)處理結(jié)果表明，4個(gè)接收站位地下表現(xiàn)出深淺不一的高阻異常，且視電阻率擬斷面圖與傳統(tǒng)地震和電磁法反演表現(xiàn)為相似的異常形態(tài)。

2)海洋可控源電磁法拖曳式發(fā)射和坐底式發(fā)射均能用于擬斷面圖快速成像算法，兩種方式處理結(jié)果表現(xiàn)為幾乎相同的電性異常，但后者比前者計(jì)算速度更快。

3)該算法能夠很好地識(shí)別地下高阻異常，橫向分辨率較高，運(yùn)算時(shí)間短，能夠快速分析海底電性分布生成視電阻率擬斷面圖，為海洋可控源電磁法定性解釋提供了一種簡(jiǎn)單快速且直觀的算法。