羅鳴， 裴建新， 葉益信

1 中國海洋大學(xué)海洋地球科學(xué)學(xué)院, 青島 266100 2 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室, 青島 266100 3 東華理工大學(xué)地球物理與測(cè)控技術(shù)學(xué)院, 南昌 330013

0 引言

海洋可控源電磁法(Controlled-Source ElectroMagnetic method, CSEM)是一種有效探測(cè)海底油氣及礦產(chǎn)資源的海洋地球物理勘探方法(Constable，2010；柳建新等，2021).雖然海洋CSEM方法已被廣泛地應(yīng)用于洋中脊和大陸邊緣地質(zhì)構(gòu)造研究以及海底油氣資源和天然氣水合物探測(cè)中，但該方法對(duì)各向異性地層的分辨能力依然不明確(Abubakar et al.，2010；Wirianto et al.，2010；周建美等，2017; Li et al.,2020).電磁場(chǎng)關(guān)于地下介質(zhì)電阻率靈敏度的計(jì)算方法是求取電磁場(chǎng)分量對(duì)介質(zhì)電阻率的導(dǎo)數(shù)，靈敏度分析是一種了解地電模型電磁場(chǎng)對(duì)探測(cè)目標(biāo)分辨能力的有效方法，傳統(tǒng)的地球物理反演方法也需要精確計(jì)算電磁場(chǎng)關(guān)于地下介質(zhì)電阻率的靈敏度.層狀介質(zhì)的正反演是最基本且應(yīng)用廣泛的數(shù)值模擬方法，同時(shí)能為二維和三維數(shù)值模擬提供參考.目前，電磁學(xué)家們已開發(fā)出了層狀介質(zhì)海洋可控源電磁響應(yīng)解析解的計(jì)算方法(Chave，1983；Everett and Constable，1999；Tompkins，2005；L?seth and Ursin,2007；Hunziker et al.，2015).由此，研究層狀介質(zhì)靈敏度的解析計(jì)算方法是有必要且可行的.

海底地層常常呈現(xiàn)電阻率各向異性(Li et al., 2013；殷長春等，2014；羅鳴和李予國，2015；蔡紅柱等，2015；Yang and Qin，2020)，垂直各向異性(VTI)是自然界中廣泛存在的各向異性情況，VTI反演是目前主要的電阻率各向異性反演解釋方法之一(Constable et al.,1987; Zhdanov et al., 2014；彭榮華等，2019).然而，在已有的地球物理文獻(xiàn)中，VTI反演存在高阻薄層橫向電阻率重構(gòu)效果不佳的問題(Ramananjaona，2011；羅鳴等，2016；趙寧等，2017).對(duì)于該問題，鮮有文獻(xiàn)從靈敏度分析的角度討論可控源電磁對(duì)電阻率各向異性薄層的分辨能力(Constable and Weiss，2006；Christensen and Dodds，2007；Medina et al.，2010).Brown等(2012)通過分析層狀VTI模型電磁響應(yīng)的有效異常討論了可控源電磁法對(duì)高阻薄層的分辨能力，然而，有效異常是關(guān)于各向同性模型電磁響應(yīng)與各向異性模型電磁響應(yīng)差異的參數(shù)，其本質(zhì)為正演響應(yīng)的特征分析，而非靈敏度分析.

目前，關(guān)于層狀介質(zhì)的靈敏度解析計(jì)算的文獻(xiàn)并不多，Chave(1984)、Flosadóttir和Constable(1996)、Key(2009)推導(dǎo)了一維電阻率各向同性介質(zhì)的靈敏度解析計(jì)算公式.Streich和Becken(2011)推導(dǎo)了垂直各向異性介質(zhì)柱坐標(biāo)系下x方向和z方向電偶極源的靈敏度計(jì)算公式并分析了陸地模型的靈敏度特征，但文獻(xiàn)中未推導(dǎo)y方向電偶極源的電磁場(chǎng)表達(dá)式和靈敏度解析公式.在模擬和解釋海洋CSEM資料時(shí)，地球物理數(shù)值模擬常在笛卡爾直角坐標(biāo)系下進(jìn)行，且通常假定發(fā)射源為理想的水平電偶極源.然而，在實(shí)際的海洋可控源電磁勘探作業(yè)中，由于洋流、波浪、潮汐等海水運(yùn)動(dòng)的影響，發(fā)射源可能會(huì)發(fā)生旋轉(zhuǎn)和傾斜等.復(fù)雜姿態(tài)的電偶極源可分解為三個(gè)相互正交方向的電偶極源，通過計(jì)算并矢量疊加三個(gè)方向發(fā)射源分量的電磁場(chǎng)，即可獲得復(fù)雜姿態(tài)下電偶極源產(chǎn)生的電磁場(chǎng)(羅鳴等，2017).因此，僅利用單一水平方向和垂直方向電偶源的電磁場(chǎng)和靈敏度計(jì)算公式并不能滿足復(fù)雜數(shù)值模擬的要求.

本文基于L?seth和Ursin(2007)提出的柱坐標(biāo)系下電阻率垂直各向異性層狀介質(zhì)單一水平方向電偶極源和垂直電偶極源電磁場(chǎng)的計(jì)算方法，推導(dǎo)了笛卡爾直角坐標(biāo)系下三個(gè)正交方向電偶極源電磁場(chǎng)，并在此基礎(chǔ)上，詳細(xì)導(dǎo)出了電磁場(chǎng)分量關(guān)于各向異性電導(dǎo)率的靈敏度解析表達(dá)式.該方法適用于任意位置激發(fā)和任意位置接收的觀測(cè)系統(tǒng).通過模擬電阻率各向異性海洋地電模型的電場(chǎng)靈敏度，分析不同方向電偶源電場(chǎng)對(duì)海底地層的橫向電阻率和垂向電阻率的靈敏度特征，比較不同方向電偶源對(duì)地層各向異性電阻率的分辨能力；基于模擬結(jié)果，進(jìn)一步分析電場(chǎng)關(guān)于海底地層各向異性率的靈敏度特征，探討可行的電阻率各向異性反演策略.

1 層狀VTI介質(zhì)海洋可控源電磁場(chǎng)表達(dá)式

(1)

其中，σh為水平方向的電導(dǎo)率，σv為垂直方向的電導(dǎo)率.在VTI介質(zhì)中，電磁場(chǎng)能夠分解為TE極化場(chǎng)和TM極化場(chǎng)，由此可得在第j層界面處的TE模式和TM模式反射系數(shù)r和透射系數(shù)t(圖1)，其表達(dá)式如下:

圖1 層狀VTI介質(zhì)模型示意圖Fig.1 Schematic diagram of layered VTI medium model

(2)

(3)

(4)

(5)

其中，μ0為真空中的磁導(dǎo)率，pρ=λ/ω，λ為Hankel變換濾波系數(shù).

圖2 反射系數(shù)R和透射系數(shù)T遞推示意圖Fig.2 Diagram of recursion of reflection coefficient R and transmission coefficient T

(6)

(7)

復(fù)雜姿態(tài)的電偶極源可分解為三個(gè)相互正交方向的電偶極源，通過計(jì)算并矢量疊加三個(gè)方向發(fā)射源分量的電磁場(chǎng)，即可獲得復(fù)雜姿態(tài)下電偶極源產(chǎn)生的電磁場(chǎng).假設(shè)發(fā)射源與接收點(diǎn)所在層具有電阻率各向同性或垂直各向異性(VTI)，依據(jù)電磁波傳播理論遞推求得傳播矩陣，進(jìn)一步可推導(dǎo)笛卡爾直角坐標(biāo)系下x、y和z方向電偶極源任意位置激發(fā)，任意接收位置的電磁場(chǎng)表達(dá)式(L?seth and Ursin,2007；羅鳴和李予國，2015).

1.1 x方向水平電偶極源(xHED, Horizontal Electric Dipole)電磁場(chǎng)表達(dá)式

假設(shè)一個(gè)x方向發(fā)射源位于(xs,ys,zs)，發(fā)射源長度為lx，發(fā)射電流為I安培，該發(fā)射源在層狀VTI介質(zhì)接收站(x,y,z)處產(chǎn)生的電磁場(chǎng)為

(8)

(9)

(10)

(11)

(12)

(13)

(14)

其中，當(dāng)X=A,a1=1,a2=1；當(dāng)X=D,a1=-1,a2=1.

②當(dāng)接收站在源的下方時(shí)，即z≥zs(圖2)，有

(15)

1.2 y方向水平電偶極源(yHED)電磁場(chǎng)表達(dá)式

假設(shè)一個(gè)y方向發(fā)射源位于(xs,ys,zs)，發(fā)射源長度為ly，發(fā)射電流為I安培，該發(fā)射源在層狀VTI介質(zhì)接收站(x,y,z)處產(chǎn)生的電磁場(chǎng)為

(16)

(17)

(18)

(19)

(20)

(21)

1.3 垂直電偶極源(VED, Vertical Electric Dipole)電磁場(chǎng)表達(dá)式

假設(shè)一個(gè)垂直方向發(fā)射源位于(xs,ys,zs)，發(fā)射源長度為lz，發(fā)射電流為I安培，該發(fā)射源在層狀VTI介質(zhì)接收站(x,y,z)處產(chǎn)生的電磁場(chǎng)為

(22)

(23)

(24)

(25)

×λ3J0(λr)dλ,

(26)

(2)針對(duì)實(shí)驗(yàn)Ⅰ,有同學(xué)提出“Fe3+產(chǎn)生的原因可能是Cl-在陽極放電,生成的Cl2將Fe2+氧化”,寫出能支持這種假設(shè)的理由和相關(guān)反應(yīng)的化學(xué)方程式:___。

在模擬VTI介質(zhì)的空間域電磁場(chǎng)響應(yīng)時(shí)，可通過數(shù)字濾波漢克爾變換實(shí)現(xiàn)由波數(shù)域到空間域的轉(zhuǎn)換(Anderson,1982).

2 偏導(dǎo)數(shù)的表達(dá)式

2.1 特征值的偏導(dǎo)數(shù)

(27)

2.2 反射系數(shù)和透射系數(shù)的偏導(dǎo)數(shù)

①m層界面的反射、透射系數(shù)對(duì)σ{h,v}m的導(dǎo)數(shù)有

(28)

(29)

和

(30)

(31)

同理地，②m-1層界面的反射、透射系數(shù)對(duì)σ{h,v}m的導(dǎo)數(shù)有

(32)

(33)

和

(34)

(35)

2.3 遞歸反射響應(yīng)和透射響應(yīng)的偏導(dǎo)數(shù)

表1 方向向下遞歸響應(yīng)的遞推范圍Table 1 Recursion parameters of downward recursive responses

(36)

(37)

(38)

其中，Em=eiωpzmdm.

(39)

(40)

(41)

表2 方向向上遞歸響應(yīng)的遞推范圍Table 2 Recursion parameters of upward recursive responses

(42)

(43)

(44)

其中，Em=eiω pzmdm

(45)

(46)

(47)

2.4 系數(shù)矩陣(RA,RD,RB,RC)的偏導(dǎo)數(shù)

系數(shù)矩陣R在不同發(fā)射源情況下取值不同：當(dāng)發(fā)射源為HED時(shí)，電場(chǎng)和磁場(chǎng)分量的R分別取RA和RD;當(dāng)發(fā)射源為VED時(shí)，電場(chǎng)和磁場(chǎng)分量的R分別取RB和RC.

①當(dāng)接收站在發(fā)射源的上方時(shí)，即z

(48)

分別對(duì)系數(shù)矩陣(RA,RD,RB,RC)進(jìn)行討論：當(dāng)X=A,a1=1,a2=1；當(dāng)X=D,a1=-1,a2=1；當(dāng)X=B,a1=1,a2=-1；當(dāng)X=C,a1=-1,a2=-1.其偏導(dǎo)數(shù)為

(49)

其中，

②當(dāng)接收站在發(fā)射源的下方時(shí)，即z>zs,

(50)

分別對(duì)系數(shù)矩陣(RA,RD,RB,RC)進(jìn)行討論：當(dāng)X=A,b1=1,b2=1；當(dāng)X=D,b1=-1,b2=1；當(dāng)X=B,b1=1,b2=-1；當(dāng)X=C,b1=-1,b2=-1.其偏導(dǎo)數(shù)為

(51)

在計(jì)算偏導(dǎo)數(shù)時(shí)，如果遞推地層為包含多界面的復(fù)合層，反射和透射響應(yīng)采用遞歸算式；如果是單一界面，則采用單一界面的反射和透射系數(shù)計(jì)算公式.

2.5 電磁場(chǎng)的偏導(dǎo)數(shù)

x方向水平電偶極源：

(52)

(53)

(54)

(54)

(56)

(57)

y方向水平電偶極源：

(58)

(59)

(60)

(61)

(62)

(63)

垂直電偶極：

(64)

(65)

(66)

(67)

(68)

3 算例

在海洋可控源電磁反演中，通常選擇電磁場(chǎng)關(guān)于地層電阻率對(duì)數(shù)的導(dǎo)數(shù)作為雅各比矩陣，因此，本文將靈敏度J定義為式(69)所示形式，并在后續(xù)的靈敏度驗(yàn)算和分析中均使用該形式計(jì)算的靈敏度進(jìn)行討論

(69)

其中，F(xiàn)為電磁場(chǎng)響應(yīng)，ρh和ρv分別為橫向和垂向電阻率.

3.1 靈敏度驗(yàn)算

為驗(yàn)證本文所提出偏導(dǎo)數(shù)算法的正確性，我們以圖3所示一維電阻率各向同性模型為例，將本文算法計(jì)算的結(jié)果與各向同性理論計(jì)算的結(jié)果進(jìn)行對(duì)比.空氣波是海洋可控源電磁響應(yīng)的主要干擾場(chǎng)之一，為了針對(duì)性地分析電磁場(chǎng)關(guān)于不同海底介質(zhì)各向異性電阻率(ρh和ρv)的靈敏度，我們將海水層設(shè)為10 km以消除空氣波的影響.假設(shè)深度為10 km的海水層電阻率為0.3 Ωm，埋深為1 km、厚度為300 m高阻薄層的電阻率為100 Ωm，覆蓋層和基巖的電阻率均為1 Ωm.假設(shè)一個(gè)電偶極源布設(shè)于海底正上方50 m處，坐標(biāo)為(0,0,9950)，76個(gè)接收站等間距地布放于y方向測(cè)線0～15 km范圍的海底(測(cè)線位置x=0).發(fā)射頻率為0.25 Hz.

圖3 一維海洋電阻率各向同性地電模型Fig.3 1D canonical isotropic conductivity model

(70)

圖4顯示了利用本文各向異性算法計(jì)算的三個(gè)正交方向電偶極源的電場(chǎng)靈敏度與各向同性算法計(jì)算結(jié)果的對(duì)比結(jié)果.圖中，藍(lán)線、紅線和綠線分別對(duì)應(yīng)電場(chǎng)關(guān)于覆蓋層、高阻層和基巖的靈敏度曲線，‘+’線為各向同性算法計(jì)算的結(jié)果(Li and Li，2017).由圖可見，二者相對(duì)誤差均非常小，由此說明，本文提出的垂直各向異性介質(zhì)靈敏度算法是正確的.

圖4 Ex(第一行)、Ey(第二行)和Ez(第三行)關(guān)于海底地層的橫向電阻率(第一列)、垂向電阻率(第二列)和各向同性電阻率(第三列)的靈敏度曲線；實(shí)線為本文算法的計(jì)算結(jié)果，‘+’線為各向同性算法的計(jì)算結(jié)果(Li and Li，2017)，第四列為二者的相對(duì)誤差.Fig.4 Sensitivities of Ex (the 1st row), Ey (the 2nd row) and Ez (the 3rd row) to the horizontal resistivities (the 1st column), vertical resistivities (the 2nd column) and the sum of them (the 3rd column) for the sediment (blue curves), resistive reservoir (red curves) and basement (green curves) layers. The 3rd column shows the comparison between the sensitivities computed from this paper (solid lines) and the isotropic algorithm (‘+’ lines) (Li and Li, 2017) and the 4th column shows the relative error between them.

3.2 靈敏度分析

目前，海洋CSEM方法主要應(yīng)用于探測(cè)數(shù)公里深度的油氣和礦產(chǎn)資源(Constable，2010)，富含油氣的地層通常呈現(xiàn)出高于圍巖的電阻率值，該類模型可近似為高阻薄層鑲嵌于低阻圍巖的地電模型(簡(jiǎn)稱為“高阻薄層模型”)；而富含硫化物等礦體的地層通常呈現(xiàn)出低于圍巖的電阻率值(Peng et al., 2020)，此類模型可近似為高導(dǎo)薄層鑲嵌于高阻圍巖的地電模型(簡(jiǎn)稱為“高導(dǎo)薄層模型”).因此，我們分別建立高阻薄層模型和高導(dǎo)薄層模型，并設(shè)定空氣和海水層為電阻率各向同性，海底介質(zhì)呈現(xiàn)電阻率各向異性.基于該兩個(gè)模型，分析不同方向電偶源和復(fù)雜姿態(tài)電偶源電場(chǎng)分量關(guān)于海底地層的橫向電阻率和垂向電阻率靈敏度特征.假設(shè)一個(gè)電偶源布設(shè)于海底正上方50 m處，151個(gè)接收站等間距地布放于y方向測(cè)線-15至15 km范圍的海底(測(cè)線位置x=0).發(fā)射頻率為0.25 Hz.

3.2.1 單一方向發(fā)射源的靈敏度

(1)y方向水平電偶源(yHED)

本小節(jié)分析yHED的靈敏度分布特征，此時(shí)發(fā)射源方向與測(cè)線方向相同,即為軸向觀測(cè)模式.為了更充分地對(duì)比和分析薄層對(duì)圍巖各向異性靈敏度的影響，我們首先分析海底介質(zhì)為高阻各向異性半空間和高導(dǎo)各向異性半空間的靈敏度特征.假設(shè)高導(dǎo)各向異性半空間模型的下半空間橫向和垂向電阻率分別為ρh=1 Ωm和ρv=4 Ωm，高阻各向異性半空間模型的下半空間橫向和垂向電阻率分別為ρh=25 Ωm和ρv=100 Ωm.圖5為高導(dǎo)各向異性半空間模型(左)和高阻各向異性半空間模型(右)水平電場(chǎng)Ey關(guān)于橫向電阻率ρh靈敏度(a)、垂向電阻率ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.其中，黑線為電場(chǎng)log10|Ey|的等值線，圖5c中白線為水平電場(chǎng)Ey關(guān)于橫向電阻率ρh的靈敏度(用?ρh表示)與水平電場(chǎng)Ey關(guān)于垂向電阻率ρv的靈敏度(用?ρv表示)相等的等值線(即?ρh/?ρv=1).由圖5a和b可見，?ρh和?ρv均在發(fā)射源處最大，并隨著收發(fā)距和深度的增大而不斷衰減；高阻各向異性半空間模型的?ρh和?ρv具有與高導(dǎo)各向異性半空間模型相似的變化趨勢(shì)；在同一半空間模型中，?ρh在水平方向的衰減速度較于?ρv更快，然而在深度方向上?ρh比?ρv更慢.由圖5c可見，高阻各向異性半空間模型海底介質(zhì)的?ρh/?ρv變化趨勢(shì)與高導(dǎo)半空間模型類似，表現(xiàn)為發(fā)射源下方?ρh/?ρv大于1，發(fā)射源兩側(cè)小于1，并隨著收發(fā)距增大而不斷衰減；在海底1 km的深度，高導(dǎo)各向異性半空間模型(圖5c1)海底介質(zhì)的?ρh/?ρv在2 km收發(fā)距處出現(xiàn)極大值(約為7)，在15 km收發(fā)距處出現(xiàn)極小值(約為0.01)，而高阻各向異性半空間模型(圖5c2)海底介質(zhì)的?ρh/?ρv在3 km收發(fā)距處出現(xiàn)極大值(約為10)，在15 km收發(fā)距處出現(xiàn)極小值(約為0.1).

圖5 高導(dǎo)半空間模型(左)和高阻半空間模型(右)yHED水平電場(chǎng)Ey關(guān)于橫向電阻率ρh靈敏度(a)、垂向電阻率ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.需注意，(a)和(b)具有相同的色棒范圍Fig.5 The distribution of sensitivities of horizontal electric field Ey with respect to ρh (a) and ρv sensitivity (b) and their ratios (c) in yoz plane for the high conductivity half-space model (left) and the high resistivity half-space model (right). Note that (a) and (b) share the same color bar

為了說明海底分別存在高阻和高導(dǎo)薄層時(shí)電磁場(chǎng)關(guān)于地下介質(zhì)的橫向電阻率和垂向電阻率的靈敏度特征，我們分別計(jì)算了高阻薄層模型和高導(dǎo)薄層模型的靈敏度.所選擇的高阻薄層模型在圖3所示模型基礎(chǔ)上，設(shè)定圍巖的橫向和垂向電阻率分別為ρh=1 Ωm和ρv=4 Ωm，薄層的橫向和垂向電阻率分別為ρh=25 Ωm和ρv=100 Ωm；而高導(dǎo)薄層模型則設(shè)定圍巖的橫向和垂向電阻率分別為ρh=25 Ωm和ρv=100 Ωm，薄層的橫向和垂向電阻率分別為ρh=1 Ωm和ρv=4 Ωm，其他參數(shù)不變.圖6為高阻薄層模型(左)和高導(dǎo)薄層模型(右)水平電場(chǎng)Ey關(guān)于橫向電阻率ρh靈敏度(a)、垂向電阻率ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.其中，黑線為電場(chǎng)log10|Ey|的等值線，圖6c中白線為?ρh/?ρv=1的等值線.對(duì)比圖5和圖6可見，當(dāng)?shù)叵陆橘|(zhì)存在具有電性差異的薄層時(shí)，?ρh和?ρv分布都發(fā)生了改變：1)薄層模型和半空間模型的?ρh分布總體差異不大，主要是層狀模型覆蓋層靈敏度和薄層所在深度靈敏度的輕微改變；2)高阻薄層模型覆蓋層?ρv與半空間模型對(duì)應(yīng)位置的?ρv(圖5a1)接近，而在薄層所在深度呈現(xiàn)出極大值，且薄層以下的圍巖?ρv略低于半空間模型(圖5b1)對(duì)應(yīng)位置的?ρv；3)高導(dǎo)薄層模型覆蓋層?ρv顯著高于半空間模型(圖5a2)對(duì)應(yīng)位置的?ρv，且在薄層所在深度呈現(xiàn)出極小值，薄層以下的圍巖?ρv略低于半空間模型(圖5b2)對(duì)應(yīng)位置的?ρv；4)高阻薄層模型的?ρh/?ρv在薄層所在深度(海底以下1 km)出現(xiàn)了10-4量級(jí)的極小異常值,即高阻薄層的?ρh遠(yuǎn)低于?ρv，而高導(dǎo)薄層模型的?ρh/?ρv在薄層所在深度(海底以下1 km)出現(xiàn)了103量級(jí)的極大異常值,即高阻薄層的?ρh遠(yuǎn)高于?ρv.

圖6 高阻薄層模型(左)和高導(dǎo)薄層模型(右)yHED水平電場(chǎng)Ey關(guān)于ρh靈敏度(a)、ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.需注意，(a)和(b)具有相同的色棒范圍Fig.6 The distribution of sensitivities of Ey with respect to ρh (a) and ρv (b) and their ratios (c) in yoz plane for resistive thin layer model (left) and conductive thin model (right). Note that (a) and (b) share the same color bar

綜合以上算例結(jié)果可知，當(dāng)海底僅存在一種介質(zhì)時(shí)，無論是高導(dǎo)介質(zhì)還是高阻介質(zhì)，電場(chǎng)關(guān)于海底介質(zhì)的橫向電阻率靈敏度?ρh和垂向電阻率靈敏度?ρv差異不大，且?ρh/?ρv沒有突變.當(dāng)?shù)叵陆橘|(zhì)存在具有電性差異的薄層時(shí)，薄層將影響圍巖電阻率靈敏度分布，且在薄層所在深度出現(xiàn)?ρh/?ρv突變，具體表現(xiàn)為：高阻薄層對(duì)低阻圍巖的靈敏度分布影響較小，高阻薄層?ρh/?ρv呈現(xiàn)出低值異常；高導(dǎo)薄層對(duì)高阻覆蓋層?ρv影響劇烈，導(dǎo)致高阻覆蓋層?ρv明顯高于半空間?ρv，且高導(dǎo)薄層?ρh/?ρv呈現(xiàn)出高值異常.

(2)x方向水平電偶源(xHED)

當(dāng)發(fā)射源取向?yàn)閤方向，此時(shí)發(fā)射源方向與測(cè)線方向垂直,即為赤道觀測(cè)模式.圖7為高阻薄層模型(左)和高導(dǎo)薄層模型(右)水平電場(chǎng)Ex關(guān)于橫向電阻率ρh靈敏度(a)、垂向電阻率ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.其中，黑線為電場(chǎng)log10|Ex|的等值線，圖7c中白線為?ρh/?ρv=1的等值線.對(duì)比圖6和圖7可見：1)赤道觀測(cè)模式的電場(chǎng)分量較于軸向觀測(cè)模式電場(chǎng)分量衰減更慢(對(duì)比電場(chǎng)等值線)；2)赤道觀測(cè)模式的?ρh、?ρv和?ρh/?ρv的變化趨勢(shì)和幅值與軸向觀測(cè)模式都非常相近，即赤道觀測(cè)模式和軸向觀測(cè)模式對(duì)海底地層的各向異性電阻率分辨能力相近.

圖7 高阻薄層模型(左)和高導(dǎo)薄層模型(右)xHED水平電場(chǎng)Ex關(guān)于ρh靈敏度(a)、ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.需注意，(a)和(b)具有相同的色棒范圍Fig.7 The distribution of sensitivities of Ey with respect to ρh (a) and ρv (b) and their ratios (c) in yoz plane for resistive thin layer model (left) and conductive thin model (right). Note that (a) and (b) share the same color bar

(3)垂直方向水平電偶源(VED)

圖8為高阻薄層模型(左)和高導(dǎo)薄層模型(右)垂直電場(chǎng)Ez關(guān)于橫向電阻率ρh靈敏度(a)、垂向電阻率ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.其中，黑線為電場(chǎng)log10|Ez|的等值線，圖8c中白線為?ρh/?ρv=1的等值線.對(duì)比圖6、圖7和圖8可見：1)垂直方向電偶源電場(chǎng)分量Ez在水平方向衰減速度均快于赤道觀測(cè)模式和軸向觀測(cè)模式電場(chǎng)分量(對(duì)比電場(chǎng)等值線)；2)垂直方向電偶源的?ρh、?ρv和?ρh/?ρv的變化趨勢(shì)與軸向和赤道觀測(cè)模式相近，但靈敏度隨著收發(fā)距的增大而衰減更快.需注意的是，對(duì)于海底各向異性地層，垂直方向電偶源?ρh/?ρv較于軸向和赤道觀測(cè)模式更加接近于1，即垂直電偶源電場(chǎng)關(guān)于海底地層的橫向電阻率和垂向電阻率靈敏度的差異較于水平電偶源更小，由此說明，垂直電偶源對(duì)海底地層各向異性電阻率的分辨能力高于水平電偶源.

圖8 高阻薄層模型(左)和高導(dǎo)薄層模型(右)VED垂直電場(chǎng)Ez關(guān)于ρh靈敏度(a)、ρv靈敏度(b)，以及二者的比值(c)在yoz平面的分布圖.需注意，(a)和(b)具有相同的色棒范圍Fig.8 The sensitivities distribution of Ey with respect to ρh (a) and ρv (b) and their ratios (c) in yoz plane for resistive thin layer model (left column) and conductive thin model (right column). Note that (a) and (b) share the same color bar

3.2.2 復(fù)雜姿態(tài)電偶極源的靈敏度

復(fù)雜姿態(tài)的電偶極源可通過計(jì)算并矢量疊加三個(gè)正交方向發(fā)射源分量的電磁場(chǎng)以獲得總電磁場(chǎng)，其對(duì)應(yīng)的靈敏度也可通過三個(gè)正交方向的電偶源所激發(fā)電磁場(chǎng)的靈敏度疊加得到：

J=JxHED+JyHED+JVED,

式中，J表示復(fù)雜姿態(tài)電偶極源的電磁場(chǎng)靈敏度，JVED表示垂直電偶極源(VED)的電磁場(chǎng)靈敏度，JxHED和JyHED分別表示沿x軸和y軸取向水平電偶極源的電磁場(chǎng)靈敏度.

假設(shè)一個(gè)y方向水平偶極源受到海流等綜合因素的影響，其姿態(tài)出現(xiàn)了旋轉(zhuǎn)和傾斜，假設(shè)其方位旋轉(zhuǎn)了30°(z軸不動(dòng)，xoy平面順時(shí)針旋轉(zhuǎn))，并發(fā)生30°傾斜(x軸不動(dòng)，yoz平面順時(shí)針旋轉(zhuǎn)).圖9為高阻薄層模型(左)和高導(dǎo)薄層模型(右)復(fù)雜姿態(tài)電偶極源Ex分量(a)、Ey分量(b)和Ez分量(c)關(guān)于ρh靈敏度與ρv靈敏度的比值在yoz平面的分布圖，其中黑線為電場(chǎng)等值線，白線為?ρh/?ρv=1的等值線.由圖9可見，由于電偶源姿態(tài)變化，Ey和Ez電場(chǎng)等值線和靈敏度均呈現(xiàn)非對(duì)稱分布；對(duì)比圖9a、b和c發(fā)現(xiàn)，Ex分量的?ρh/?ρv與Ey分量的?ρh/?ρv幅值非常接近，而Ez分量的?ρh/?ρv較于Ex和Ey分量都更加接近于1.從本算例結(jié)果發(fā)現(xiàn)，當(dāng)電偶源發(fā)生姿態(tài)變化時(shí)，不同方向電場(chǎng)分量具備相應(yīng)方向電偶源電場(chǎng)的特征；較于水平電場(chǎng)分量，垂直電場(chǎng)分量關(guān)于海底地層的橫向和垂向電阻率靈敏度差異最小，這種現(xiàn)象在薄層中更為明顯.

圖9 高阻薄層模型(左)和高導(dǎo)薄層模型(右)復(fù)雜姿態(tài)電偶極源Ex分量(a)、Ey分量(b)和Ez分量(c)關(guān)于ρh靈敏度與ρv靈敏度的比值在yoz平面的分布圖.需注意，a和b具有相同的色棒范圍Fig.9 The distribution of ratios of ? ρh and ? ρv for Ex (a), Ey (b) and Ez (c) in yoz plane for resistive thin layer model (left) and conductive thin model (right). Note that all subgraph share the same color bar

3.2.3 電磁場(chǎng)關(guān)于各向異性率λ的靈敏度

從以上算例發(fā)現(xiàn)，薄層的橫向電阻率靈敏度(?ρh)和垂向電阻率靈敏度(?ρv)差異明顯高于圍巖的?ρh和?ρv差異，且表現(xiàn)為高阻薄層?ρh遠(yuǎn)小于?ρv，而高導(dǎo)薄層?ρh遠(yuǎn)大于?ρv.?ρh和?ρv的巨大差異可能導(dǎo)致反演很難重構(gòu)出靈敏度較低的電阻率值.假設(shè)反演能夠較為準(zhǔn)確地重構(gòu)海底地層的橫向電阻率ρh或垂向電阻率ρv，此時(shí)若能夠恢復(fù)較為準(zhǔn)確的各向異性率λ，則可通過二者的關(guān)系λ重構(gòu)出靈敏度較低的電阻率值.因此，我們進(jìn)一步分析電場(chǎng)關(guān)于各向異性率λ的靈敏度特征.

電磁場(chǎng)關(guān)于各向異性率λ的靈敏度存在以下關(guān)系：

(71)

采用上圖9所示算例相同的電偶源，即yHED方位旋轉(zhuǎn)30°，并傾斜30°.圖10為高阻薄層模型(左)和高導(dǎo)薄層模型(右)復(fù)雜姿態(tài)電偶極源Ex分量(a)、Ey分量(b)和Ez分量(c)關(guān)于各向異性率λ靈敏度在yoz平面的分布圖,其中黑線為電場(chǎng)等值線.對(duì)比圖10a、b和c發(fā)現(xiàn)，高阻薄層模型三分量電場(chǎng)關(guān)于各向異性率λ的靈敏度特征與ρv的靈敏度特征(圖6c1, 圖7c1, 圖8c1)非常接近，圍巖的?Ey/?λ略高于?Ex/?λ和?Ez/?λ，需注意的是，高阻薄層的三分量電場(chǎng)關(guān)于λ的靈敏度明顯高于ρh的靈敏度；高導(dǎo)薄層模型圍巖的三分量電場(chǎng)關(guān)于各向異性率λ的靈敏度特征與ρh的靈敏度特征(圖6c2, 圖7c2, 圖8c2)接近，但高導(dǎo)薄層的三分量電場(chǎng)關(guān)于λ的靈敏度明顯高于ρv的靈敏度.從本算例結(jié)果發(fā)現(xiàn)，電場(chǎng)分量關(guān)于λ的靈敏度綜合了橫向電阻率靈敏度?ρh和垂向電阻率靈敏度?ρv的特征；高阻薄層電場(chǎng)分量關(guān)于λ的靈敏度高于ρh的靈敏度，且高導(dǎo)薄層電場(chǎng)分量關(guān)于λ的靈敏度高于ρv的靈敏度，由此可知，反演中不直接重構(gòu)低靈敏度的各向異性電阻率值，而通過反演各向異性率λ間接恢復(fù)低靈敏度的各向異性電阻率值將是一個(gè)可行的反演策略.

圖10 高阻薄層模型(左)和高導(dǎo)薄層模型(右)復(fù)雜姿態(tài)電偶極源Ex分量(a)、Ey分量(b)和Ez分量(c)關(guān)于各向異性率λ靈敏度在yoz平面的分布圖.需注意，a和b具有相同的色棒范圍Fig.10 The distribution of ratios of? ρh and ? ρv for Ex (a), Ey (b) and Ez (c) in yoz plane for resistive thin layer model (left column) and conductive thin model (right column). Note that all subgraph share the same color bar

4 結(jié)論

本文推導(dǎo)了笛卡爾直角坐標(biāo)系下，電阻率垂直各向異性介質(zhì)中三個(gè)正交方向電偶極源電磁場(chǎng)關(guān)于各向異性電導(dǎo)率的靈敏度表達(dá)式.模擬了海底介質(zhì)為高阻半空間模型和高導(dǎo)半空間模型，以及高導(dǎo)圍巖中存在高阻薄層模型和高阻圍巖中存在高導(dǎo)薄層模型的靈敏度.計(jì)算結(jié)果表明，當(dāng)海底各向異性半空間存在具有電性差異的薄層時(shí)，薄層將顯著影響地下介質(zhì)各向異性電阻率的靈敏度分布；高阻薄層所在深度橫向電阻率靈敏度和垂向電阻率靈敏度比值呈現(xiàn)低值異常；高導(dǎo)薄層引起高阻覆蓋層靈敏度明顯高于高阻半空間靈敏度，且高導(dǎo)薄層所在深度橫向電阻率靈敏度和垂向電阻率靈敏度比值呈現(xiàn)出高值異常；垂直電偶源對(duì)海底地層各向異性電阻率的分辨能力高于水平電偶源.對(duì)比各向異性率和各向異性電阻率靈敏度可見，電場(chǎng)分量關(guān)于各向異性率的靈敏度綜合了橫向電阻率靈敏度和垂向電阻率靈敏度的特征，通過反演各向異性率間接恢復(fù)低靈敏度的各向異性電阻率值是一個(gè)可行的反演策略.

致謝感謝匿名審稿人提出的寶貴意見.