深海拖曳式瞬變電磁的響應規律

周勝 ，席振銖，宋剛，李波，左立標，侯海濤

(1. 中南大學 地球科學與信息物理學院，湖南 長沙，410083；2. 長沙五維地科勘查技術有限責任公司，湖南 長沙，410205；3. 北京先驅高科技開發公司，北京，100081)

海底熱液硫化物礦是繼發現錳結核后的又一種海底礦產資源，賦存于水深1.2～3.7 km的熱液活動區，富含Cu，Zn，Pb，Ag，Au，Co，Cd和Mn等金屬，主要分布在環太平洋構造帶中的西太平洋，以及西南太平洋的弧后擴張盆地、東太平洋海隆和大西洋洋中脊區域，堆積在長度為數百米、寬度為近百米的范圍內，品位高、埋深淺，資源開發前景廣闊，已成為各國海洋資源調查的重點之一[1]。針對深海堆積型熱液硫化礦床的特定環境，實施“快速、有效、便于實施”的深海拖曳式瞬變電磁探測系統十分必要。深海拖曳式瞬變電磁是借鑒陸地瞬變電磁法(Transient electromagnetic methods)，具有快速、有效和便于實施的特點，其整個系統主要包括工作船、水下儀器艙拖體單元和水下發射及接收線圈拖體單元 3部分。Edwards等[2?5]利用電偶源瞬變電磁偶極偶極方法進行了海底電導率填圖，開創了海洋瞬變電磁的先河；劉長勝等[6]計算了海水為均勻半空間的瞬變電磁響應；李慧[7]計算了淺海瞬變電磁垂直偶極子和中心回線裝置的響應曲線，并與陸地情況對比分析了海水對感應電動勢的影響；Li等[8]對淺海瞬變電磁法進行了一維數值模擬，研究了天波對時間域可控源電磁(CSEM)響應的影響。在此，本文作者建立深海熱液金屬硫化物礦床的地電模型，運用瞬變電磁全空間理論計算方法，模擬計算深海瞬變電磁響應，總結其規律。

1 深海堆積型熱液硫化礦的電性模型

在深海環境中，海底堆積型熱液硫化礦被海水包圍，其主要圍巖為海底沉積物、碎玄武巖和完整玄武巖。深海近海底水溫基本保持在 2～4 ℃。根據Accerboni等[9]提出的海水電導半經驗估算公式可以得到深海近海底海水的電導率變化范圍為 3.2～3.4 S/m。海底沉積物、碎玄武巖和完整玄武巖的電導率主要取決于孔隙度、滲透率、溫度、年代和各向異性等因素，其變化范圍分別為 0.750～2.00，0.100～0.500和 0.030～0.001 S/m(表 1)。海底熱液硫化礦電導率變化范圍為0.800～10.000 S/m，這取決于礦物類型、金屬含量、溫度和孔隙密度等因素。但在通常情況下，堆積型熱液硫化礦的電導率大約是海水的3倍，是海底沉積物和基巖的5倍多，具備基于電導率為基礎的電法勘探的物性條件，可以運用瞬變電磁法等方法方法有效探測深海堆積型熱液硫化礦。海水的浸泡作用使得近海底電性比裸露在空氣介質中近地表電性均勻得多。根據海水電導率和深海堆積型熱液硫化礦及其圍巖電性特征，以及深海堆積型熱液型硫化礦的分布特征，可以建立如圖1所示的近似層狀地電模型。

表1 海底巖礦石電導率Table 1 Conductivities of rocks on seafloor

圖1 深海堆積型熱液硫化礦電性模型Fig.1 Electrical model of deep-sea hydrothermal sulfide ore deposit

2 深海瞬變電磁計算方法

深海環境中瞬變電磁傳播是一個全空間過程，見圖2。當水平發送回線處于深海環境時，電磁波既向海底地層中傳播，也向上在海水中傳播，由于海水不同于空氣，其電阻率低，感應產生二次場很強，因此，海水感應的二次場的影響是不可忽略的。Krivochieva等[10?11]討論了層狀介質中瞬變電磁全空間計算方法，借鑒其計算方法，分析深海環境下瞬變電磁場的傳播過程(見圖2)，建立深海堆積型熱液硫化礦層狀的電性模型(見圖1)，模擬計算拖曳式深海瞬變電磁響應。需說明的是：深海拖曳式瞬變電磁發送回線的面積是有限的，相對于在深海水平方向上有100 m延伸的硫化物礦床，發送回線的面積小得多，可以將深海硫化物礦床簡化為層狀模型進行模擬計算。

如模型4，當發射源位于海水中的第i分層時，矢量位A*在每一層介質中都滿足亥姆霍茲方程[12?16]，并且A*只有Z分量。

圖2 全空間瞬變場的傳播示意圖Fig.2 Diagram of whole space transient spread

圖3 瞬變電磁全空間層模型Fig.3 Whole space layer model of transient electromagnetic

則：

其中：j為層數；PM=ISN為磁偶極距；I為發射電流；S為發射線圈面積；N為發射線圈匝數；Cj和Dj為待求系數，當z→?∞時，D1=0；當z→∞時，Cn=0，根據邊界條件：

建立方程組，可以求得 Zj=?i ω μ/uj，j=1，2，3，…，n；Z(0)=Z0；Z(5)=Z5；

由快速漢克爾變換[14]計算頻率域電磁場Hz，用余弦變換多項式法計算

3 深海拖曳式瞬變電磁響應規律

建立一維深海地電模型，如圖1所示。取導電率σ0=0.000 01 S/m，深度 h0→?∞；σ1=3 S/m，h1=2 km；σ2=3 S/m；σ4=0.1 S/m，h4=500 m；σ5=0.002 S/m，h5→∞，I=10 A，發送回線半徑為5 m，回線匝數50匝，采用重疊回線裝置。當發射源位于海底即 h2=h3=0 m，不存在低阻熱液硫化礦層時，得到深海瞬變電磁響應曲線和不考慮海水影響的瞬變電磁響應如圖4所示。從圖4可以看出：海水的作用對瞬變電磁響應曲線形態影響較小，曲線形態與不考慮海水影響的瞬變電磁的形態基本相似，但海水的存在導致 TEM 響應衰減較慢，響應幅值增大，尤其是在晚延時階段，深海瞬變電磁相比無海水情況的瞬變電磁響應幅值增加了近 2個數量級。當發射源位于海底且存在硫化礦層時，取h2=0 m，h3=20 m，得到如圖5所示的瞬變電磁響應曲線。從圖5可見：深海瞬變電磁響應與不考慮海水影響的瞬變電磁響應曲線的變化趨勢一致，但是，深海瞬變電磁早延時和晚延時階段響應幅值顯著增強，為不考慮海水情況的2倍多，海洋與不考慮海水影響的瞬變電磁響應在1～100 ms時窗內，異常基本一致。

圖4 無礦層瞬變電磁響應對比曲線Fig.4 Contrast response curves of transient electromagnetic sounding

圖5 含礦層時瞬變電磁響應對比曲線Fig.5 Contrast transient electromagnetic response curves of hydrothermal sulfide ores

4 拖曳高度對瞬變電磁響應的影響規律

如圖1所示，取電導率 σ0=0.000 01 S/m，深度h0→?∞；σ1=3 S/m，h1=2 000 m；σ2=3 S/m；h2=H(H為拖曳高度)；σ3=10 S/m，h3=20 m，σ4=0.1 S/m，h4=500 m；σ5=0.002 S/m，h5→∞，發送回線半徑為5 m，回線匝數50匝，采用重疊回線裝置。當發送回線位于海底即H=0 m時，硫化礦層的異常響應非常明顯(見圖6(a))，從早延時到大約100 ms的觀測窗口都可以觀測到低阻礦層的深海瞬變電磁響應，比深海背景場響應(即不存在低阻熱液硫化礦層是的瞬變電磁響應)高2～3倍；隨著拖曳高度的增加(見圖6(b)，(c)和(d))，深海低阻礦層的瞬變電磁異常響應向后延時，并迅衰減速接深海背景場。拖曳高度在0～100 m范圍內變化，對29.75 ms時間道進行抽道，如圖7所示。從圖7可見：隨著拖曳高度的升高，深海低阻層瞬變電磁響應與深海背景場迅速接近。理論上，隨著拖曳高度的增加，通過延時觀測可以達到探測海底硫化礦的目的，但是，海水的電阻率極低，電磁波的衰減嚴重，當拖曳高度超過 50 m時，異常響應幅值只有約為背景場值的7%，已經相當微弱。為了突出海底淺部硫化礦層的異常特征，拖曳高度應盡量接近海底，不宜超過50 m。

圖6 不同拖曳高度下深海瞬變電磁響應對比曲線Fig.6 Contrast of transient electromagnetic responses to different towed heights

圖7 不同拖曳高度下29.75 ms含礦層瞬變電磁響應Fig.7 Transient electromagnetic response of different towed heights at 29.75 ms including ore layer

5 深海拖曳式瞬變電磁響應的時窗范圍

在圖1中，設參數與前面的相同，拖曳高度H是可變的，當拖曳高度H為10 m和50 m時含礦層與海洋背景瞬變電磁響應對比曲線如圖8所示。通過分析圖8發現：海水的影響使得瞬變電磁響應在小于關斷延時1 ms時，難以與海洋背景場相區分；從關斷5 ms左右開始，逐步顯示熱液硫化礦層的異常，據此可以清楚地分辨硫化礦層異常響應；關斷100 ms以后，瞬變電磁響應又逐漸趨于海洋背景場，難以分辨熱液硫化礦層的異常響應。可見：針對深海熱液硫化礦，深海拖曳式瞬變電磁測深最佳觀測窗口為 1～100 ms時窗。

圖8 不同拖曳高度下深海瞬變電磁響應曲線Fig.8 Transient electromagnetic responses to different towed heights

6 結論

(1) 深海與無海水瞬變電磁相比，響應衰減曲線的形態相似，但前者明顯衰減緩慢且響應幅值明顯增強。

(2) 拖曳高度對海洋瞬變電磁探測的影響不可忽略。隨著拖曳高度的升高，低阻層的異常響應向后延遲，并迅速衰減，最后幾乎淹沒在海洋背景場中。當拖曳高度超過50 m時，異常響應極其微弱，在實際應用中將難以觀測。為了探明深海熱液硫化礦，深海拖曳式瞬變的拖曳高度應低于50 m。

(3) 海水的存在使得深海瞬變電磁衰減緩慢，在拖曳高度不高于50 m時，為了突出深海熱液硫化礦的瞬變電磁響應，拖曳式瞬變電磁的觀測時窗應滿足涵蓋關斷后延時1～100 ms的時間窗口。

致謝：承蒙牛之璉教授和鮑光淑教授對本文提出了寶貴的修改意見，在此致以誠摯的謝意！

[1]何繼善, 鮑力知. 海洋電磁法研究的現狀和進展[J]. 地球物理學進展, 1999, 14(1): 7?39.HE Ji-shan, BAO Li-zhi. The situation and progress of marine electromagnetic method research[J]. Progress in Geophysics,1999, 14(1): 7?39.

[2]Edwards R N, Chave A D. Atransient electric dipole-dipole method for mapping the conductivity of the sea floor[J].Geophysics, 1986, 51(4): 984?987.

[3]Cheosman S J, Edwards R N, Chave A D. On the theory of sea-floor conductivity mapping using transient electromagnetic systems[J]. Geophysics, 1987, 52(2): 204?217.

[4]Edwards R N. Two-dimensional modeling of a towed in-line electric dipole-dipole sea-floor electromagnetic system: The optimum time delay or frequency for target resolution[J].Geophysics, 1988, 53(6): 846?853.

[5]楊建文, Edwards R N. 用于海底電導率填圖的可控源時域電磁法[J]. 中國有色金屬學報, 1998, 8(4): 705?713.YANG Jian-wen, Edwards R N. Controlled source time-domain electromagnetic methods for sea-floor conductivity mapping[J].The Chinese Journal of Nonferrous Metals, 1998, 8(4): 705?713.

[6]劉長勝, 林君. 海底表面磁源瞬變響應建模及海水影響分析[J]. 地球物理學報, 2006, 49(6): 1891?1898.LIU Chang-sheng, LIN Jun. Transient electromagnetic response modeling of magnetic source on seafloor and the analysis of seawater effect[J]. Chinese J Geophys, 2006, 49(6): 1891?1898.

[7]李慧. 海洋瞬變響應理論計算及淺海底瞬變電磁探測技術研究[D]. 長春: 吉林大學儀器科學與電氣工程學院, 2006:27?48.LI Hui. Theoretical calculation of marine transient electromagnetic response and research of shallow sea-floor detection technology[D]. Changchun: Jilin University. College of Instrumentation and Electrical Engineering, 2006: 1?140.

[8]Li Y G, Constable S. Transient electromagnetic in shallow water:insights from 1D modeling[J]. Chinese J Geophys, 2010,53(3): 737?742.

[9]Accerboni E, Mosetti F. A physical relationship among salinity,temperature and electrical conductivity of sea water[J]. Boll Geofis Teor Appl, 1967, 34(9): 87?96.

[10]Krivochieva S, Chouteau M. Whole space modeling of a layered earth in time-domain electromagnetic measurements[J]. Journal of Applied Geophysics, 2002, 50(4): 375?391.

[11]HU Bo, YUE Jian-hua, YANG Hai-yan. Algorithm study of transient response of vertical magnetic bipolar source in whole space plane layered medium[J]. Mining Science and Technology 2009, 19(1): 120?123.

[12]牛之璉. 時間域電磁法原理[M]. 長沙: 中南大學出版社,2007: 37?66.NIU Zhi-lian. The theory of time domain electromagnetic method[M]. Changsha: Central South University Press, 2007:37?66.

[13]方文藻, 李予國, 李貅. 瞬變電磁測深原理[M]. 西安: 西北工業大學出版社, 1993: 9?50.FANG Wen-zao, LI Yu-guo, LI Xiu. Principles of transient electromagnetic sounding[M]. Xi’an: Northwestern Polytechnical University Press, 1993: 9?50.

[14]考夫曼 A A, 凱勒 G V. 頻率域與時間域電磁測深[M]. 王建謀, 譯. 北京:地質出版社, 1987: 216?279.Kaufman A A, Keller G V. Frequency and transient soundings[M]. WANG Jian-mou, transl. Beijing: Geology Press,1987: 216?279.

[15]米薩克 N·納比吉安. 勘查地球物理電磁法: 第一卷 理論[M].趙經祥, 譯. 北京: 地質出版社, 1992: 195?211.Nabighian M N. Electromagnetic methods in applied geophysics:Volume 1, Theory[M]. ZHAO Jing-xiang, transl. Beijing:Geology Press, 1992: 195?211.

[16]Guptasama D, Singh B. New digital linear filters for Hankel J0 and J1 transforms[J]. Geophysical Prospecting, 1997, 45(4):745?762.