肖鵬飛， 齊文杰*， 高永才， 徐鳳嬌， 李弘， 俞建寶

(1.中國石油化工股份有限公司石油物探技術(shù)研究院， 南京 211103； 2.成都理工大學(xué)地球物理學(xué)院， 成都 610059)

可控源音頻大地電磁測深法(controlled-source audio magnetotelluric method，CSAMT)，是在大地電磁測深法(magnetotelluric sounding，MT)和音頻大地電磁測深法(audio magnetotelluric method，AMT)的基礎(chǔ)上發(fā)展起來的一種人工源頻率域測深方法。近年來，隨著地質(zhì)找礦向縱深發(fā)展，具有探測深度大、勘探能力強(qiáng)、分辨率高等特點(diǎn)的可控源音頻大地電磁法，在油氣勘查、金屬礦普查、水文及環(huán)境地質(zhì)等多個領(lǐng)域顯示出了廣闊的應(yīng)用前景[1]。CSAMT法觀測電磁場的頻率、場強(qiáng)和方向可由人工控制，克服了因天然場相對較弱、人文干擾較大時觀測的困難，其極化方向明顯，信噪比高，易于觀測[2]。

但同時也帶來一些負(fù)面影響，如與場源相關(guān)的場源效應(yīng)。CSAMT的場源效應(yīng)主要表現(xiàn)為：由于靠近場源而產(chǎn)生的非平面波效應(yīng)；由于場源下面的地質(zhì)情況或者場源與測深點(diǎn)之間的地質(zhì)情況而產(chǎn)生的場源附加效應(yīng)；由于場源與測深點(diǎn)之間的地質(zhì)異常體被投射而產(chǎn)生的陰影效應(yīng)等[3]。其中，最為突出的是近場效應(yīng)問題，直接制約了CSAMT法的有效探測深度，嚴(yán)重影響了CSAMT法的勘探能力[4]。底青云等[5]認(rèn)為CSAMT資料處理的質(zhì)量依賴于近場校正的效果。

中外學(xué)者針對CSAMT方法的近場校正技術(shù)進(jìn)行了相關(guān)研究。通常的近場校正方法有如下幾種。

一是三角形校正法。其原理是利用波區(qū)和近區(qū)公式計(jì)算的視電阻率曲線在中間區(qū)形成三角形，其形態(tài)與大地電阻率無關(guān)，只與收發(fā)距有關(guān)，這種方法屬于經(jīng)驗(yàn)性的，對于簡單的地質(zhì)情況可以得到較好的效果，但對于復(fù)雜地質(zhì)情況，其校正效果還難以評價[4]。

二是K值校正法。該法利用均勻大地表面電磁場的精確表達(dá)式和波區(qū)漸進(jìn)表達(dá)式的比值，計(jì)算不同收發(fā)距時的修正系數(shù)，然后對卡尼亞電阻率進(jìn)行校正，這種方法的校正效果類似三角形校正法，對層狀介質(zhì)的理論計(jì)算也不夠完善，且增加了磁場垂直分量的觀測[6]。

三是山下改正法。該法也是目前CSAMT非遠(yuǎn)場改正中常用的方法之一，雖然改正后的曲線較未改正曲線有接近MT曲線的趨勢，但應(yīng)指出這種方法在過渡區(qū)改正不足，而在近場區(qū)又改正過頭，導(dǎo)致測深曲線發(fā)生畸變，分析認(rèn)為這是由于將均勻大地CSAMT非遠(yuǎn)區(qū)場的改正方法應(yīng)用于層狀大地非遠(yuǎn)區(qū)場的改正，由于理論前提不同而導(dǎo)致改正效果欠佳[7]。

四是遠(yuǎn)場因子-近場因子(Kf-Kn)系數(shù)校正法。該法是由加拿大鳳凰公司推出的，由于該法本身的近似性及其隨觀測條件而變化，使其校正精度也不高，在地層分布較復(fù)雜時，校正結(jié)果與MT曲線偏差較大，分析認(rèn)為這也是由于將均勻大地的改正方法應(yīng)用于層狀大地所致[8]。

五是全頻域視電阻率法。當(dāng)CSAMT測量不滿足波區(qū)條件時，運(yùn)用電場和磁場分量與收發(fā)距、方位角和電極距之間的關(guān)系，對Cagniard電阻率的計(jì)算公式進(jìn)行了修正，近年來這一校正方法得到了較為廣泛的關(guān)注和應(yīng)用[9]。

張凱等[2]國內(nèi)學(xué)者認(rèn)為，所有近場校正方法都不是最有效和最完善的，處理后的近場成分還是無法消除，仍然影響著對深部信息的提取和認(rèn)識。在這方面，楊毅等[10]于2009年利用積分方程法研究了CSAMT中存在的源效應(yīng)問題，認(rèn)為源效應(yīng)問題非常復(fù)雜，不存在簡單的校正方法，實(shí)際工作應(yīng)盡量在遠(yuǎn)區(qū)。

鑒于近場校正問題的復(fù)雜性和困難程度，基于一維正演模擬，借鑒Kf-Kn系數(shù)校正法的思路，現(xiàn)從如下兩方面入手，一方面，提出視電阻率數(shù)據(jù)近場校正與反演處理協(xié)同聯(lián)進(jìn)的資料處理策略——即反演遞歸校正法，以往所有的近場校正方法都是將校正與反演分割開來單獨(dú)進(jìn)行的，考慮到近場影響難以一步校正到位的現(xiàn)實(shí)問題，采用分步收斂校正策略，即校正結(jié)果與反演地電模型以高度耦合的方式，在同步收斂中最終實(shí)現(xiàn)比較理想的校正和反演處理效果；另一方面，考慮到在很多工區(qū)存在一定量的測井、地震、地質(zhì)等已知信息，為降低CSAMT實(shí)測資料處理結(jié)果的非唯一性，又提出基于構(gòu)建研究區(qū)層狀背景地電模型的近場校正方案。最后，以一條地?zé)峥碧綄?shí)測剖面的數(shù)據(jù)處理為例，通過以上兩種校正方法的聯(lián)合運(yùn)用，驗(yàn)證方法的可行性和有效性。旨在進(jìn)一步深化對近場校正方法的研究，為推動CSAMT實(shí)測資料處理技術(shù)的進(jìn)步做出探討。

1 基于均勻半空間正演模擬的查表校正法

1.1 方法原理

對CSAMT勘探，理論上只要使收發(fā)距r足夠大，就可以使觀測點(diǎn)處于遠(yuǎn)區(qū)或波區(qū)，但是在實(shí)際工作中，由于發(fā)射功率的限制，加之電磁場隨r的增大以高次冪而衰減，所以當(dāng)r足夠大以后，往往使電磁場的信號變得微弱，不能保證觀測精度。因此，在r有限大的情況下，在采用MT常規(guī)方法處理實(shí)測資料之前，需要做近場和過渡場校正處理。

均勻半空間地面水平電偶極子源的電磁場[8]的坐標(biāo)方位如圖1所示，設(shè)電偶極子AB與X軸方向一致，地面上任一點(diǎn)M與AB中點(diǎn)O之距離為r，OM連線與AB線(即X軸)之間的夾角為θ，均勻半空間電阻率為ρ1，在多層模型中，ρ1代表第一層電阻率。

圖1 水平偶極及坐標(biāo)示意圖Fig.1 Schematic diagram of horizontal dipoles and coordinates

在遠(yuǎn)區(qū)有

(1)

(2)

在近區(qū)有

(3)

(4)

式中：i為復(fù)數(shù)；Hy為磁場分量；Ex為電場分量；I為供電電流強(qiáng)度，A；AB為供電偶極長度，m；r為場源到接收點(diǎn)之間的距離，km。

在近區(qū)情況下，由式(3)有|Hy|∝1/r2。

導(dǎo)入因子Kf/5，令

(5)

導(dǎo)入近場因子Kn，令

(6)

通過計(jì)算遠(yuǎn)場因子Kf和近場因子Kn與fr|Hy/Ex|(發(fā)射頻率、收發(fā)距與阻抗倒數(shù)模三者的乘積，下文簡稱綜合因數(shù))的關(guān)系，如圖2所示，其關(guān)系竟出乎意料的簡單，Kn在近場時為一常數(shù)0.63，Kf在遠(yuǎn)場為一常數(shù)1.0。從實(shí)際運(yùn)用的f、r和測定的|Hy/Ex|值可以把場分為遠(yuǎn)場區(qū)、近場區(qū)和過渡場區(qū)。

在遠(yuǎn)場區(qū)，因?yàn)镵f=1，表明式(5)與Cagniard方程相同，而在近場區(qū)和過渡場區(qū)則必須用Kf和Kn系數(shù)進(jìn)行校正。

圖2 Kf和Kn與fr|Hy/Ex|關(guān)系曲線Fig.2 Kf and Kn and fr|Hy/Ex| relation curve

1.2 校正步驟

CSAMT一維正演模擬采用的是斯克里普斯海洋研究所Kerry Key公布的Dipole1D程序[11-13]。所述近場校正方法的研究是以赤道偶極(θ=90°)情形為例展開的(其他θ角度值觀測數(shù)據(jù)的校正方法是類似的)，且僅討論應(yīng)用遠(yuǎn)場因子Kf對實(shí)測視電阻率進(jìn)行校正的問題。基于均勻半空間正演模擬的查表校正法之?dāng)?shù)據(jù)處理過程如下。

步驟1計(jì)算遠(yuǎn)場因子Kf和綜合因數(shù)fr|Hy/Ex|的對應(yīng)關(guān)系數(shù)據(jù)表。

首先根據(jù)生產(chǎn)施工(或正演模擬)所使用的發(fā)射—觀測系統(tǒng)，對模型文件設(shè)置如下幾何參數(shù)：包括發(fā)射偶極長度、電偶極子計(jì)算的高斯求積積分點(diǎn)數(shù)、發(fā)射偶極XYZ坐標(biāo)及其水平角和傾角、測線端點(diǎn)XYZ坐標(biāo)等。測線端點(diǎn)的布置不僅要覆蓋實(shí)際測線觀測范圍，還要適當(dāng)超出一定距離，以盡量擴(kuò)大綜合因數(shù)分布區(qū)間，方便后續(xù)對Kf校正值做查表插值計(jì)算(若實(shí)際綜合因數(shù)大小超出了數(shù)據(jù)表中的區(qū)間范圍，則要采用外推方法計(jì)算Kf值，這將會導(dǎo)致校正誤差的增大)。

然后設(shè)置物性參數(shù)。均勻半空間的物性參數(shù)僅有一個，即電阻率值ρ1，但在計(jì)算Kf和綜合因數(shù)的對應(yīng)關(guān)系數(shù)據(jù)表時，同樣是出于盡量擴(kuò)大綜合因數(shù)區(qū)間范圍的目的，研究10、100、1 000、10 000 Ω·m由低到高幾種有代表性的均勻半空間情形。

當(dāng)模型參數(shù)文件設(shè)置完成后，執(zhí)行Dipole1D程序，對構(gòu)建好的一系列均勻半空間地電模型，分別進(jìn)行正演模擬(發(fā)射頻率采用實(shí)測數(shù)據(jù)文件中的頻率表)，最后將這幾種模型計(jì)算得到的Kf校正值，合并成一個總表，并按綜合因數(shù)由小到大的順序排序后備用，Kf具體曲線形態(tài)參見圖2中的藍(lán)色曲線。

步驟2對實(shí)測視電阻率數(shù)據(jù)進(jìn)行場源影響校正處理。

首先，從實(shí)際運(yùn)用的發(fā)射頻率f、收發(fā)距r和測定的Ex和Hy電磁場數(shù)據(jù)，計(jì)算得到各測點(diǎn)各頻率的實(shí)際綜合因數(shù)值；然后，通過查表經(jīng)內(nèi)插外推確定Kf值；最后，將各實(shí)測視電阻率值與求得的Kf相乘，便完成了對實(shí)測數(shù)據(jù)的近場校正處理。

CSAMT法的野外觀測資料通過近場和過渡場校正以后，就得到了相當(dāng)于遠(yuǎn)場的視電阻率結(jié)果，于是就可以運(yùn)用目前已完善的 MT 或 AMT的所有資料處理及解釋方法，開展后續(xù)的研究工作了。

1.3 算例分析

為驗(yàn)證上述查表法對CSAMT視電阻率數(shù)據(jù)的場源影響校正效果，設(shè)計(jì)了3個地電模型，各模型正演模擬數(shù)據(jù)及其校正結(jié)果如下。

1.3.1 模型一：Q型三層地電模型

模型的幾何參數(shù)和物性參數(shù)分別為h1=300 m(第一層厚度)，ρ1=10 000 Ω·m(第一層電阻率)，h2=300 m(第二層厚度)，ρ2=1 000 Ω·m(第二層電阻率)，ρ3=100 Ω·m(第三層電阻率)，發(fā)射偶極長度為4 000 m，分別對收發(fā)距r為4 km和8 km兩種情形進(jìn)行了研究，模擬及校正結(jié)果曲線如圖3所示。

一維正演模擬結(jié)果，參見圖3(a)、圖3(b)中的紅色曲線，通過與MT參考曲線對比可知，當(dāng)r=4 km時，從200 Hz開始，ρs曲線便進(jìn)入了過渡場區(qū)，并在最低頻率0.125 Hz的近場區(qū)，曲線上升至2 000 Ω·m；當(dāng)r=8 km時，ρs曲線則是從12.5 Hz才逐漸進(jìn)入過渡場區(qū)的，而在最低頻點(diǎn)曲線也只上升到了600 Ω·m左右。可見，在其他條件不變的情況下，適當(dāng)增大收發(fā)距，便可明顯降低進(jìn)入過渡場區(qū)的頻率閾值及近場影響程度，有利于觀測到更大頻率區(qū)間的遠(yuǎn)場數(shù)據(jù)。

查表法Kf校正結(jié)果，參見圖3(a)、圖3(b)中的藍(lán)色曲線，通過與MT參考曲線對比可知，兩種收發(fā)距情況下，都取得了非常理想的校正效果，而r=8 000 m時，校正曲線與MT曲線重合度更高。分析認(rèn)為，視電阻率數(shù)據(jù)校正效果上佳的原因，主要是由Kf與綜合因數(shù)關(guān)系表的適用性決定的。由本算例，可以把收發(fā)距越大校正效果越佳這一特征，作為評價過渡場和近場校正過程中，所使用上述關(guān)系表質(zhì)量優(yōu)劣的一種判斷標(biāo)準(zhǔn)，即滿足這一特征說明Kf與綜合因數(shù)關(guān)系表質(zhì)優(yōu)、適用性好，否則說明質(zhì)劣、適用性差。

圖3 Q型三層地電模型模擬數(shù)據(jù)及場源影響校正結(jié)果Fig.3 Forward modeling data of Q type three-layer geoelectric model and the field source influence correction results

1.3.2 模型二：H型三層地電模型

模型的幾何參數(shù)和物性參數(shù)分別為第四系(Q)和第三系(R)地層厚度h1=300 m，電阻率ρ1=300 Ω·m；白堊系(K)地層厚度h2=300 m，電阻率ρ2=100 Ω·m；侏羅系(J)地層電阻率ρ3=500 Ω·m，發(fā)射偶極長度為4 000 m，分別對收發(fā)距r為4 km和8 km兩種情形進(jìn)行了研究，模擬及校正結(jié)果曲線如圖4所示。

一維正演模擬結(jié)果，參見圖4(a)、圖4(b)中的紅色曲線，通過與MT參考曲線對比可知，當(dāng)r=4 km時，從100 Hz開始，ρs曲線便進(jìn)入了過渡場區(qū)，然后隨頻率降低而高角度上升，并在最低頻率0.125 Hz，視電阻率值高達(dá)31 500 Ω·m以上；當(dāng)r=8 km時，ρs曲線則從55 Hz逐漸進(jìn)入過渡場區(qū)，在最低頻點(diǎn)曲線上升為12 500 Ω·m左右。這與模型一中所分析收發(fā)距對觀測數(shù)據(jù)的影響特征是類似的。需要說明的是，在H型地電模型的模擬曲線中，出現(xiàn)了“負(fù)尖峰”現(xiàn)象，這種“負(fù)尖峰”的出現(xiàn)主要是由于電場在過渡區(qū)的變化特點(diǎn)所致。當(dāng)電場由遠(yuǎn)場進(jìn)入過渡區(qū)時，平面波向非平面波轉(zhuǎn)化，使得電場變小、曲線下降。當(dāng)繼續(xù)降低頻率時，由于基底高阻地層的影響，使得電場曲線形成一個極小點(diǎn)后逐漸增大。這樣相應(yīng)地就在CSAMT視電阻率曲線上產(chǎn)生了頻率范圍很窄的“負(fù)尖峰”[7]。

查表法Kf校正結(jié)果，參見圖4(a)、圖4(b)中的藍(lán)色曲線，通過與MT參考曲線對比可知，兩種收發(fā)距情況下，都取得了比較理想的校正效果，而r=8 km時，校正曲線的尾支與MT曲線重合度更高(缺點(diǎn)是對“負(fù)尖峰”畸變數(shù)據(jù)沒有起到校正作用)。根據(jù)上述對模型一校正效果的分析，可以認(rèn)為算例在過渡場和近場校正過程中，所采用的Kf數(shù)據(jù)表質(zhì)優(yōu)、適用性好，同樣體現(xiàn)了收發(fā)距越大校正效果越佳這一顯著特征。

圖4 H型三層地電模型模擬數(shù)據(jù)及場源影響校正結(jié)果Fig.4 Forward modeling data of H type three-layer geoelectric model and the field source influence correction results

1.3.3 模型三：電性呈“高-低-高-低”分布的四層地電模型

模型的幾何參數(shù)和物性參數(shù)分別為h1=200 m，ρ1=100 Ω·m，h2=200 m，ρ2=50 Ω·m，h3=900 m，ρ3=10 000 Ω·m，薊縣系(Jx)地層電阻率ρ4=100 Ω·m，發(fā)射偶極長度為4 000 m，分別對收發(fā)距r為4 km和8 km兩種情形進(jìn)行了研究，模擬及校正結(jié)果曲線如圖5所示。

一維正演模擬結(jié)果，參見圖5(a)、圖5(b)中的紅色曲線，通過與MT參考曲線對比可知，當(dāng)r=4 km時，從300 Hz開始ρs曲線便進(jìn)入了過渡場區(qū)，跨過“負(fù)尖峰”后，隨頻率降低而高角度上升，并在最低頻率0.125 Hz，視電阻率值高達(dá)63 000 Ω·m以上；當(dāng)r=8 km時，ρs曲線則是從250 Hz逐漸進(jìn)入過渡場區(qū)的，同樣在跨過頻帶稍寬的“負(fù)尖峰”后，隨頻率降低以稍低的斜率增大，在最低頻點(diǎn)曲線上升為31 500 Ω·m左右。這與模型二中所分析收發(fā)距對觀測數(shù)據(jù)影響特征是類似的。

查表法Kf校正結(jié)果，參見圖5(a)、圖5(b)中的藍(lán)色曲線，通過與MT參考曲線對比可知，兩種收發(fā)距情況下，校正效果均較差，特別是低頻段偏差較大，r=4 km時，在最低頻點(diǎn)偏差約0.7個對數(shù)單位；而r=8 km時，這一偏差則達(dá)到了0.85，同樣對“負(fù)尖峰”畸變數(shù)據(jù)沒有起到校正作用。據(jù)此，可以認(rèn)為算例在過渡場和近場校正過程中，所采用的Kf數(shù)據(jù)表質(zhì)良差、適用性欠佳，出現(xiàn)了隨著收發(fā)距增大而校正偏差亦隨之增大的異?，F(xiàn)象。為解決這一問題，又提出了如下兩種校正方法：反演遞歸校正法和設(shè)計(jì)層狀背景模型校正法，希望能在一定程度上改進(jìn)Kf數(shù)據(jù)表的質(zhì)量，從而達(dá)到提高場源影響校正效果的目的。

圖5 四層地電模型模擬數(shù)據(jù)及場源影響校正結(jié)果Fig.5 Forward modeling data of four-layer geoelectric model and the field source influence correction results

2 反演遞歸校正法

基于均勻半空間正演模擬的查表校正法所列算例表明，簡單的三層地電模型取得了比較理想的校正效果，而對稍為復(fù)雜的四層模型，則校正效果欠佳，分析認(rèn)為出現(xiàn)這種實(shí)驗(yàn)結(jié)果的原因，主要是由Kf數(shù)據(jù)表的質(zhì)量，亦即其對所論地電模型的適用性所造成的。為了提高對實(shí)測數(shù)據(jù)的校正效果，重點(diǎn)應(yīng)在改進(jìn)Kf數(shù)據(jù)表的質(zhì)量，基于這種認(rèn)識，提出并開展了反演遞歸校正法的研究。

2.1 校正步驟

所謂反演遞歸校正，就是以均勻半空間Kf視電阻率校正結(jié)果為基礎(chǔ)，通過一維OCCAM反演[14-15]，用層狀介質(zhì)取代均勻半空間來重新計(jì)算Kf數(shù)據(jù)表，從而對實(shí)測數(shù)據(jù)開展進(jìn)一步場源影響校正的方法，具體數(shù)據(jù)處理過程如下。

步驟1對均勻半空間Kf校正結(jié)果做一維MT反演，需要說明的是，對于H型曲線，由于校正后的視電阻率仍然存在“負(fù)尖峰”，要對這種畸變數(shù)據(jù)進(jìn)行編輯預(yù)處理后再做MT反演。

步驟2采用一維反演結(jié)果，建立CSAMT層狀地電模型，即形成Dipole1D程序所需要的模型參數(shù)控制文件，然后進(jìn)行正演計(jì)算，得到模擬的電磁場各分量場值和視電阻率等數(shù)據(jù)。

步驟3采用一維迭代反演的擬合數(shù)據(jù)，作為地電模型響應(yīng)理論數(shù)據(jù)，用于標(biāo)定計(jì)算Kf值，從而得到新的校正系數(shù)數(shù)據(jù)表。

步驟4使用該Kf數(shù)據(jù)表，對觀測數(shù)據(jù)再次進(jìn)行近場和過渡場校正。

步驟5如果需要，可以繼續(xù)對校正結(jié)果進(jìn)行一維MT反演，然后重復(fù)步驟2～步驟5，直至獲得滿意的校正結(jié)果或者校正結(jié)果已基本不再變化為止。

2.2 算例分析

對于模型二，按上述介紹的數(shù)據(jù)處理步驟，對模型二的模擬數(shù)據(jù)做了進(jìn)一步的近場校正處理，如圖6所示?？梢钥闯?，新的校正結(jié)果較均勻半空間Kf初步校正結(jié)果效果更佳，表現(xiàn)在中頻段和低頻段與MT參考曲線偏差更小，特別是在曲線尾支，反演遞歸校正結(jié)果較均勻半空間Kf校正結(jié)果，在降低偏差方面取得了明顯的效果，缺點(diǎn)仍然是對“負(fù)尖峰”畸變數(shù)據(jù)沒有起到校正作用。

由此可初步認(rèn)為，這種場源影響校正與MT反演緊密結(jié)合的CSAMT資料處理策略是可行且有效的。當(dāng)然，這種有效性是建立在初步校正時，Kf數(shù)據(jù)表的基本適用前提下的。

對于模型三，同樣的，也對模型三的模擬數(shù)據(jù)做了進(jìn)一步的校正處理，如圖7所示?？梢钥闯觯碌男ＵY(jié)果較均勻半空間Kf初步校正結(jié)果有所改進(jìn)，主要表現(xiàn)在中頻段與MT參考曲線偏差有所減小，但這種減小幅度有限，總體上改進(jìn)效果不明顯，而且對“負(fù)尖峰”畸變數(shù)據(jù)沒有起到校正作用。本算例沒有取得同模型二相近的改進(jìn)效果，分析認(rèn)為這主要是由初步校正時，Kf數(shù)據(jù)表的質(zhì)量差、不適用造成的。

兩個算例的實(shí)驗(yàn)結(jié)果表明，反演遞歸校正法是建立在均勻半空間Kf校正有效性基礎(chǔ)之上的，即初步校正時所采用的Kf數(shù)據(jù)表適用性若較好，則反演遞歸校正效果將會進(jìn)一步提高；反之，若初步校正時所采用的Kf數(shù)據(jù)表質(zhì)量較差，則會對后續(xù)校正效果，產(chǎn)生比較嚴(yán)重的制約作用，這也是本方法在處理實(shí)測資料時應(yīng)予以注意的一個局限性。

圖6 模型二校正結(jié)果對比Fig.6 Comparison of correction results in Model 2

圖7 模型三校正結(jié)果對比Fig.7 Comparison of correction results in Model 3

3 設(shè)計(jì)層狀背景模型校正法

反演遞歸校正法所列算例表明，與均勻半空間Kf校正結(jié)果相比，效果有所改進(jìn)，特別是當(dāng)初步校正所采用Kf數(shù)據(jù)表質(zhì)量較好時，這種改進(jìn)效果還是較為顯著的；但情況相反時，則改進(jìn)效果不明顯。為了進(jìn)一步提高對實(shí)測數(shù)據(jù)的處理能力，需要提高校正所用Kf數(shù)據(jù)表的質(zhì)量，提出了設(shè)計(jì)層狀背景模型的校正方法。

3.1 校正步驟

所謂設(shè)計(jì)層狀背景模型校正法，就是根據(jù)研究區(qū)以及鄰區(qū)所掌握的各種已知資料，通過整理、分析等建立起區(qū)內(nèi)的層狀背景地電模型，以計(jì)算出適用于本區(qū)的、質(zhì)量較高的Kf數(shù)據(jù)表，從而促進(jìn)對CSAMT實(shí)測數(shù)據(jù)場源影響校正效果的提高，具體步驟如下。

步驟1根據(jù)研究區(qū)和鄰區(qū)已經(jīng)掌握的地質(zhì)、物性、測井、地震等各種資料，對研究區(qū)探測深度范圍內(nèi)的主要地層開展各地層代表巖性物性資料的整理、分析工作，當(dāng)所收集資料不充足時，要補(bǔ)充一定量的野外樣品采集和實(shí)驗(yàn)室物性測試數(shù)據(jù)，以結(jié)合區(qū)內(nèi)各主要地層所統(tǒng)計(jì)的平均厚度等成果，完成區(qū)內(nèi)層狀背景地電模型的構(gòu)建工作。

步驟2采用上述構(gòu)建完成的層狀背景地電模型，形成Dipole1D程序所需要的模型參數(shù)控制文件，然后進(jìn)行正演計(jì)算，得到模擬的電磁場各分量場值和視電阻率等數(shù)據(jù)。

步驟3采用上述層狀背景地電模型的MT一維正演視電阻率數(shù)據(jù)，來標(biāo)定計(jì)算Kf值，形成適用于研究區(qū)的校正系數(shù)數(shù)據(jù)表。

步驟4使用該Kf數(shù)據(jù)表，對區(qū)內(nèi)實(shí)測的CSAMT視電阻率數(shù)據(jù)進(jìn)行校正處理。

3.2 算例分析

為近似驗(yàn)證，通過設(shè)計(jì)層狀背景模型的方式，對觀測數(shù)據(jù)的校正效果，對模型三的地電參數(shù)做了一些調(diào)整，將調(diào)整后的地電參數(shù)作為設(shè)計(jì)層狀背景模型的參數(shù)(表1)，以從理論角度初步研究這種校正方法的有效性和可行性。

表1 模型三真實(shí)地電參數(shù)與設(shè)計(jì)層狀模型參數(shù)對比表Table 1 Real geoelectric parameters of Model 3 and design layered background model parameters

由表1中的設(shè)計(jì)地電參數(shù)，標(biāo)定計(jì)算并形成Kf數(shù)據(jù)表后，對模型三的模擬數(shù)據(jù)，校正處理的結(jié)果如圖8所示?？梢钥闯?，設(shè)計(jì)背景模型校正結(jié)果與MT模擬數(shù)據(jù)，在曲線形態(tài)上有很好的可比性，特別是在高頻段和中頻段兩者重合性較好，僅尾支偏差較大(這與設(shè)計(jì)模型給定的ρ4偏大有很大關(guān)系)，而且對反演遞歸校正結(jié)果基本沒有改正的“負(fù)尖峰”畸變數(shù)據(jù)，也給出了很好的校正結(jié)果。

值得說明的是，算例僅是通過對真實(shí)地電參數(shù)做了有限幅度的調(diào)整，來作為設(shè)計(jì)模型開展近似研究而得到的校正結(jié)果。在實(shí)際應(yīng)用中，具體校正效果如何，將完全取決于所構(gòu)建層狀背景模型對研究區(qū)真實(shí)地電分布的近似程度，毫無疑問近似程度越高則近場校正效果越好，反之亦然。

圖8 模型三校正結(jié)果對比Fig.8 Comparison of correction results in Model 3

4 實(shí)例效果分析

以河北省廊坊市地?zé)豳Y源CSAMT勘探實(shí)測數(shù)據(jù)為例介紹場源影響校正的效果。CSAMT法選用的儀器為美國Zonge公司生產(chǎn)的GDP-32系統(tǒng)，發(fā)射機(jī)為Zonge公司的GGT-10大功率發(fā)射機(jī)，選用頻段為0.125～8 192 Hz。所處理L03測線數(shù)據(jù)，施工參數(shù)為：發(fā)射電偶極AB長3 km，收發(fā)距r為7 km，測線長4.2 km，點(diǎn)距為50 m。工作目的是了解勘探區(qū)段內(nèi)斷裂構(gòu)造的分布狀況及其位置，同時對地層進(jìn)行劃分。工區(qū)為第四系地層所覆蓋，地形平坦，交通便利，人文干擾相對較小。實(shí)測數(shù)據(jù)經(jīng)編輯平滑、靜態(tài)校正、空間濾波等預(yù)處理后，視電阻率等值線擬斷面圖如圖9所示。

可以看出，經(jīng)各項(xiàng)預(yù)處理后的CSAMT數(shù)據(jù)，當(dāng)頻率低于10 Hz時，卡尼亞視電阻率等值線分布密集，且量值隨頻率降低呈高梯度迅速增大趨勢，這便是典型的場源影響特征。對此，采用如下步驟，對場源影響予以校正。

步驟1根據(jù)所收集到的鄰區(qū)地?zé)峋⑶叭宋镄越y(tǒng)計(jì)研究成果等資料，經(jīng)綜合分析，構(gòu)建工作區(qū)層狀背景模型如下：第四系(Q)和第三系(R)地層厚度h1=500 m，電阻率ρ1=30 Ω·m；白堊系(K)地層厚度h2=700 m，電阻率ρ2=50 Ω·m；侏羅系(J)地層厚度h3=900 m，電阻率ρ3=100 Ω·m；薊縣系(Jx)地層電阻率ρ4=200 Ω·m。其中，K、J兩地層為熱儲層Jx的蓋層，Jx為本工區(qū)的熱儲層。

步驟2采用該層狀背景地電模型，先后通過Dipole1D程序和MT一維正演的模擬計(jì)算，形成適用于工作區(qū)的初始校正系數(shù)Kf數(shù)據(jù)表，并對圖9所示的CSAMT視電阻率數(shù)據(jù)逐測點(diǎn)進(jìn)行場源影響初步校正處理。

步驟3對校正結(jié)果做MT一維OCCAM單次迭代反演，由反演結(jié)果建立CSAMT層狀地電模型，然后使用Dipole1D程序正演模擬結(jié)果和迭代反演的擬合數(shù)據(jù)，來逐測點(diǎn)標(biāo)定計(jì)算更新后的Kf數(shù)據(jù)表，并采用該質(zhì)量優(yōu)化后的校正系數(shù)表對實(shí)測數(shù)據(jù)逐測點(diǎn)進(jìn)行場源影響遞歸校正處理。

步驟4重復(fù)步驟3，直至各測點(diǎn)反演迭代滿足如下條件則終止迭代：①均方根擬合差小于0.1；②相鄰兩次迭代擬合差相對變化小于1.0%；③迭代次數(shù)達(dá)到最大值10。

經(jīng)上述遞歸校正和迭代反演處理，最終的卡尼亞視電阻率場源影響校正結(jié)果和電阻率反演斷面分別如圖10和圖11所示。

由圖10可知，經(jīng)遞歸校正后，在擬斷面中頻率低于10 Hz時，卡尼亞視電阻率等值線分布密度明顯降低，且量值隨頻率降低的增大梯度趨于平緩——即場源對CSAMT數(shù)據(jù)的近場影響特征已基本消失，亦即比較理想的將CSAMT數(shù)據(jù)校正到了平面波場數(shù)據(jù)，這種數(shù)據(jù)就可以采用MT成熟的資料反演技術(shù)來處理了。

圖9 CSAMT視電阻率擬斷面圖Fig.9 CSAMT apparent resistivity pseudo-section

圖10 場源影響校正結(jié)果視電阻率擬斷面圖Fig.10 Field source influence correction result apparent resistivity pseudo section diagram

圖11 地電—地質(zhì)綜合解釋斷面圖Fig.11 Geoelectric-geological comprehensive interpretation section

由圖11可知，從電阻率反演斷面成果圖來看，共顯示出4個大的電性層，結(jié)合鄰區(qū)鉆孔資料，推測了各電性層的地質(zhì)屬性。第一層：電阻率介于30～130 Ω·m，其深度為0～700 m，推測為第四系(Q)地層和第三系地層(R)的反映，該層電阻率變化較大，與淺表電性不均勻有關(guān)。第二層：電阻率介于50～120 Ω·m，其深度為450～1 700 m，推測為白堊系(K)的反映。第三層：電阻率介于100～140 Ω·m，其深度為830～2 700 m，推測為侏羅系(J)的反映。第四層：電阻率大于130 Ω·m，推測為薊縣系地層(Jx)的反映，頂板埋深1 050～2 700 m。同時，根據(jù)反演結(jié)果電阻率梯級帶的分布，對斷層構(gòu)造進(jìn)行了推測，在圖11中分別標(biāo)示為F1、F2和F3，3條斷層將CSAMT地電斷面分成4個斷塊。對熱儲層薊縣系來講，地層頂板埋深起伏較大，總的趨勢是小號點(diǎn)方向高、大號點(diǎn)方向低。

通過以上對CSAMT視電阻率數(shù)據(jù)的場源影響遞歸校正和迭代反演處理，基本查明了勘探區(qū)各地層的賦存情況和斷裂構(gòu)造的空間展布特征，為區(qū)內(nèi)后續(xù)地?zé)豳Y源的勘探提供了可靠的物探依據(jù)。

5 結(jié)論與認(rèn)識

(1)CSAMT勘探，為了保障觀測精度，人工場源與接收點(diǎn)的距離不可能達(dá)到無窮遠(yuǎn)，這樣實(shí)測數(shù)據(jù)難免受到過渡場與近區(qū)場的影響，為簡化、且有效地處理解釋CSAMT資料，非遠(yuǎn)場視電阻率數(shù)據(jù)的校正方法研究，目前而言仍然是解決實(shí)際生產(chǎn)問題所不可忽視、且亟待提高的一項(xiàng)重要實(shí)用技術(shù)。

(2)實(shí)現(xiàn)的基于均勻半空間正演模擬的Kf查表校正法，對于簡單的地電模型模擬數(shù)據(jù)，取得了比較理想的校正效果，分析認(rèn)為這主要是由Kf數(shù)據(jù)表對所論模型的適用性決定的，即當(dāng)所用Kf數(shù)據(jù)表適用性好、質(zhì)量優(yōu)良時，就會取得好的校正效果，反之亦然。

(3)反演遞歸校正法是在均勻半空間Kf初步校正基礎(chǔ)上，希望以層狀大地介質(zhì)代替均勻半空間，為進(jìn)一步提高校正效果而提出來的，理論算例表明，這種策略是可行且有效的，特別是當(dāng)初步校正所用Kf數(shù)據(jù)表質(zhì)量較好時，校正效果的改進(jìn)還是比較顯著的，當(dāng)然若情況相反，則改進(jìn)不明顯。

(4)設(shè)計(jì)層狀背景模型校正法是希望以層狀大地介質(zhì)代替均勻半空間，為進(jìn)一步提高Kf數(shù)據(jù)表的適用性而提出的另外一種場源影響校正策略，這種校正方法是建立在對工區(qū)所掌握的各種地質(zhì)、物性等資料綜合研究成果基礎(chǔ)之上的，即所掌握的已知資料越充分、綜合研究程度越高，所構(gòu)建的研究區(qū)層狀背景模型對真實(shí)的大地電性分布近似性越好，則對實(shí)測視電阻率數(shù)據(jù)的校正效果就會越佳，反之亦然。

(5)對CSAMT實(shí)測資料場源影響校正方法的研究是在一維正演模擬的基礎(chǔ)上開展的，場源影響校正問題并不是一個孤立的問題，而應(yīng)把它放在CSAMT勘探理論與實(shí)踐這個整體之中來加以系統(tǒng)考慮，譬如將研究區(qū)層狀背景模型構(gòu)建以及MT一維反演等環(huán)節(jié)與場源影響校正問題結(jié)合起來加以系統(tǒng)研究的思路，就是希望將場源影響校正問題以更高的鍥合度，通過系統(tǒng)整體的思路，來促進(jìn)這一實(shí)用技術(shù)研究水平及解決實(shí)際生產(chǎn)問題能力的協(xié)同提高。

(6)通過CSAMT地?zé)峥碧綄?shí)例，簡要介紹了設(shè)計(jì)層狀背景模型校正法和反演遞歸校正法在實(shí)測數(shù)據(jù)處理過程中的聯(lián)合運(yùn)用，并初步驗(yàn)證了場源影響校正與迭代反演協(xié)同處理策略的可行性和有效性。需要說明的是，所提出的場源影響校正方案僅適用于因觀測靠近場源而產(chǎn)生的非平面波效應(yīng)。然而，可控源音頻大地電磁測深法受場源影響因素較多，比如由場源下面的地質(zhì)情況或者場源與測深點(diǎn)之間的地質(zhì)情況而產(chǎn)生的場源附加效應(yīng)、由場源與測深點(diǎn)之間的地質(zhì)異常體被投射而產(chǎn)生的陰影效應(yīng)等。因此，所提出的場源影響校正方法的實(shí)用性與可行性仍需進(jìn)一步的探討和實(shí)驗(yàn)研究。