張量CSAMT數(shù)據(jù)處理技術初步研究與示范應用

孟慶奎，林品榮，李 勇，李建華，朱宏偉，李 蕩

1.中國國土資源航空物探遙感中心，北京 10083 2.中國地質科學院地球物理地球化學勘查研究所，河北 廊坊 065000

?

張量CSAMT數(shù)據(jù)處理技術初步研究與示范應用

孟慶奎1，林品榮2，李 勇2，李建華2，朱宏偉1，李 蕩2

1.中國國土資源航空物探遙感中心，北京 10083 2.中國地質科學院地球物理地球化學勘查研究所，河北 廊坊 065000

面對大深度和三維探測，標量CSAMT(controlled source audio-frequency magnetotellurics)已顯力不從心，而張量CSAMT正逐漸成為研究熱點。在前人研究基礎上，筆者就張量CSAMT數(shù)據(jù)處理方面提出Mohr圓分析和傳統(tǒng)坐標旋轉兩種處理技術，并首次創(chuàng)新性地使用我國自主研發(fā)的綜合電法儀開展了張量CSAMT示范應用。結果表明：Mohr圓分析和傳統(tǒng)坐標旋轉法相結合應用效果良好，可定性挖掘地下地質體的縱向和橫向構造特征，主軸方位的視電阻率和阻抗相位與目標體對應良好，并且電阻率反演結果與鉆孔等地質資料綜合分析對建立試驗區(qū)的地質地球物理模型是有益的。

張量CSAMT；坐標旋轉；Mohr圓；示范應用

0 引言

在音頻大地電磁/大地電磁(AMT/MT)和標量可控源音頻大地電磁(標量CSAMT)的理論基礎之上，張量可控源音頻大地電磁(張量CSAMT)法得以逐步發(fā)展。20世紀50年代初，基于天然源平面電磁波垂直入射大地及其趨膚效應原理，前蘇聯(lián)的Tichonov和法國的Cagniard提出了大地電磁法。70年代初，加拿大多倫多大學Strangberg教授和他的研究生Goldstein提出了加載人工源的音頻大地電磁法(CSAMT)[1-4]。自此，CSAMT得到了長足發(fā)展和廣泛應用，并多集中于標量測量方法和處理技術；然而直到目前，國內(nèi)外仍鮮有關于張量CSAMT的研究報道。大深度和三維探測已成大勢所趨[5-12]，張量CSAMT作為一種探測地質信息的有效手段，正逐步引起國內(nèi)外專家學者的重視。其中具有代表性的為Li等[13-16]的“張量可控源大地電磁法”，給出了該方法的定義，詳細說明了可控電偶源阻抗張量和傾子向量的推導過程，并最終得出了結果表達式。Boerner 和Wright[17]將張量CSAMT運用于紐芬蘭中部布切恩斯礦山的勘探中，取得了顯著的效果。2011年，德國Metronix地球物理研究所Bernhard Friedrichs教授[18]率領團隊成功研制了世界上第一套可實用化開展張量CSAMT測量的綜合電磁法儀(GMS-07e)，并詳細介紹了其在探測金屬硫化物礦產(chǎn)中的應用。從此揭開了張量CSAMT發(fā)展新篇章。本文通過Mohr圓分析和傳統(tǒng)坐標旋轉兩種處理方法，對張量CSAMT數(shù)據(jù)處理技術進行了初步研究。

1 張量CSAMT數(shù)據(jù)處理技術

針對張量CSAMT的信號特點，基于現(xiàn)已成熟的MT，本文圍繞以下兩種方法——坐標旋轉法和Mohr圓分析法，對張量CSAMT數(shù)據(jù)進行處理。

1.1 坐標旋轉法

在坐標旋轉過程中，阻抗張量元素的變換關系如下：

(1)

式中：Zij和Zij(θ)(i,j=x,y)分別為坐標旋轉前和旋轉后的阻抗張量元素，阻抗張量元素有其明確的物理意義，Zij表示在場源作用下，i方向產(chǎn)生的電場分量加權線性疊加結果和j方向產(chǎn)生的磁場分量加權線性疊加結果的比值；θ為坐標系順時針旋轉角度。式(1)展開如下：

(2)

其中：

(3)

根據(jù)主軸上阻抗張量元素所具有的特點，借助解析法可得主軸方位為

tg4θ0=

(4)

式中：θ0為主軸方位；*表示相應量的共軛復數(shù)。

據(jù)式(4)求出主軸方位之后，可分別得到主軸上的波阻抗、視電阻率和阻抗相位：

(5)

(6)

(7)

式中：ZTM、ZTE分別為TM和TE模式下的波阻抗；ρxy、ρyx分別為TM和TE模式下的視電阻率；φxy、φyx分別為TM和TE模式下的阻抗相位；T為周期；imag和real分別為取虛部和取實部函數(shù)。

引用MT中兩個常用的判別參數(shù)二維偏離度和橢率對地下異常體進行維度判定。

①二維偏離度

MT中定義的二維偏離度S如下：

(8)

由于式(8)中分母、分子都與坐標方位無關，所以可用其表達地電結構的特征。鑒于在二維介質中S=0，三維介質中S>0;故S值越小，其體現(xiàn)的異常體的二維性越明顯。一般地，當S≤0.5時，可認為其近似是二維的。

②橢率

阻抗張量元素在復平面內(nèi)，對于三維介質的情況，其隨著θ角變化的軌跡為橢圓；對于二維介質，相應的橢圓退化為直線；對于一維介質，相應的直線退化為點。橢圓函數(shù)在復平面上短軸B和長軸A之比稱為橢率β。當β=1時，地下異常體是等軸狀三維構造；當β由1逐漸變?yōu)?時，地下異常體由等軸狀三維構造逐漸演變?yōu)槎S構造。

1.2 Mohr圓分析法

Mohr圓可將MT阻抗張量不變量有機地結合在一起，因為MT阻抗張量不變量具有信息的多樣性和穩(wěn)定性等優(yōu)點，所以采用Mohr圓作為研究MT阻抗張量性質的圖示工具，能夠簡便、直觀地得到二維偏離度和各向異性等重要信息，可為研究地電結構橫向和縱向的變化提供便利[19-21]。

然而，在礦山、城市附近等地區(qū)很難有效地應用Mohr圓分析法，因為MT使用天然場源，信號微弱，易受人為噪聲干擾。而張量CSAMT采用的是人工場源，可有效提高信噪比，引入Mohr圓分析法，可為地下地質構造解釋提供有效依據(jù)。

測量軸的阻抗張量為Z，順時針旋轉θ角后變?yōu)閆′，將式(2)中第一項整理如下：

(9)

令

(10)

則

(11)

式中：

(12)

統(tǒng)一對其他分量進行上述變換，結果列舉如下：

(13)

(14)

半徑為

(15)

圖1 第一類實部Mohr圓圖示Fig.1 Real part Mohr circle of the first kind

當圓心C偏離橫軸時，體現(xiàn)地下構造的三維性，且偏離橫軸愈遠，三維性愈強；當圓心C在橫軸上且Mohr圓半徑不為0時，則地下構造是二維的；當Mohr圓退化為橫軸上的點時，地下構造為一維情形。

依據(jù)式(14)和式(15)繪制的Mohr圓可作為一種有效的圖示分析工具，用來研究阻抗張量不變量的性質。就單個測點而言，據(jù)不同頻率(f)得到的Mohr圓，可分析地下構造的縱向變化特征。就整條剖面而言，按上述方式畫出所有測點的Mohr圓，再對比測點間相同頻率Mohr圓的位置和形態(tài)，可討論地下構造的橫向變化特征。

2 張量CSAMT在新疆某礦區(qū)的應用示范

圖2 試驗區(qū)及鄰區(qū)地質圖Fig.2 Geologic map of experiment and nearby area

試驗區(qū)位于新疆中天山地塊(Ⅱ級)東南緣。區(qū)內(nèi)出露的主要地層為中元古界的中、深變質巖和上古生界的火山熔巖及火山碎屑巖。研究區(qū)中部偏東存在一斷裂F4，近似弧形展布，出露長度大于8 km，整體走向約60°，航片可見明顯線性構造，產(chǎn)狀向北西陡傾。此條斷裂為含礦基性-超基性巖漿提供上升通道，是研究區(qū)內(nèi)最主要的導巖和容礦構造。斷裂北側普遍存在一套片理化的云母斜長片巖和變粒巖系，以及一套糜棱巖化的花崗質巖系，后者中存在若干含銅鎳雜巖體。本次張量CSAMT針對∑19號雜巖體開展采集工作，該雜巖體出露于相對低洼地帶，呈長方形展布，總體走向約80°，出露地表長度約300 m，寬度為10～150 m。雜巖體產(chǎn)狀變化較大，北緣向南陡傾，南緣向北陡傾，傾角約85°～88°。經(jīng)鉆孔驗證雜巖體呈巖墻產(chǎn)出，其北側圍巖為糜棱巖化花崗巖，南側為片麻狀花崗巖。雜巖體淺部具全巖蝕變，深部與圍巖接觸的橄欖巖具弱礦化，部分可達邊界品位。試驗區(qū)及鄰區(qū)地質狀況如圖2所示。

表1為測定的巖(礦)電性參數(shù)統(tǒng)計。區(qū)內(nèi)超基性巖相對電阻率一般為87～1 857 Ω·m，異常段極化率值一般為1.05%～1.41%，最高可達2.74%，呈低阻高極化特征；輝長巖、片巖、花崗巖等均呈高阻低極化特征。

表1 巖(礦)石電性參數(shù)統(tǒng)計

2.1 工作部署及參數(shù)設置

試驗區(qū)張量CSAMT法工作布置見圖3。張量CSAMT測線長度為600 m。為使測區(qū)測點接收到最強的張量CSAMT電磁場信號，供電點A、B和C的位置選擇滿足理論要求。

圖3 試驗區(qū)張量CSAMT工作布置圖Fig.3 Layout of tensor CSAMT in experiment area

圖4 張量CSAMT工作裝置示意圖Fig.4 Tensor CSAMT working device schematic diagram

張量CSAMT供電與接收裝置見圖4。發(fā)射端采用“L”型源，供電電極AB平行于測線，垂直于供電電極BC；接收端采用十字型，M1N1沿測線布置，M2N2垂直于測線布置。各參數(shù)設置見表2。

2.2 阻抗張量不變量及Mohr圓結果分析

用阻抗張量不變量實部繪制的各測點(214--270)Mohr圓如圖5所示。分析圖5可知，測點218、222所有頻點的Mohr圓圓心均較遠地偏離直線Zxxr=0，表明這兩個測點附近的地下電性結構具有較強的三維性； 而其他測點的Mohr圓圓心都落在直線Zxxr=0附近，表明這些測點附近的地下電性結構趨于二維構造。該推斷與下述傳統(tǒng)坐標旋轉法的分析結果相符。 214、 218、 222、 230 等測點的Mohr圓半徑明顯比其他測點的更大，可判斷其各

向異性程度更為強烈。就單個測點而言，Mohr圓半徑總是隨頻率的減小而減小，可推斷測區(qū)內(nèi)各項異性程度于淺部較大、而深部較小，這體現(xiàn)了近地表巖體，特別是雜巖體上部呈全巖蝕變的地質現(xiàn)象。在∑19號雜巖體的中心地段，242測點各頻點的Mohr圓圓心幾近重合，這與它附近測點的Mohr圓形態(tài)存在顯著差別，且238測點和246測點對應頻點的Mohr圓以242測點呈對稱分布，該現(xiàn)象有待進一步研究。

2.3 坐標旋轉法結果分析

主軸方位的視電阻率和阻抗相位擬合斷面如圖6a,b所示。分析圖6a,b可知：238測點至250測點呈低阻高相位特征，對比地質剖面草圖(圖6c)可知，該段恰好為∑19號雜巖體主體地段，特別在242測點附近視電阻率達100 Ω·m，阻抗相位達45°。而據(jù)表1可知，雜巖體中輝長巖呈高值電阻率特性，而這一特性已被完全掩蓋，故可推測雜巖體內(nèi)部存在銅鎳礦的可能性極大。目前，這一推測已由ZK7-1號垂直鉆孔資料所證實。雜巖體南側254測點附近呈較低視電阻率、較高阻抗相位特性，體現(xiàn)了雜巖體與圍巖接觸帶處的礦化蝕變現(xiàn)象；266測點附近也呈較低視電阻率、較高阻抗相位特性，推斷為糜棱巖化花崗巖弱礦化所致；據(jù)視電阻率擬斷面圖推斷，∑19號雜巖體產(chǎn)狀近于直立，這與已知的地質資料相吻合。

通過對214--270各測點二維偏離度曲線的分析，可知218、222、230三個測點的二維偏離度數(shù)值較大，而其他測點的二維偏離度均小于0.5，據(jù)此推測小號測點附近地下介質主要呈三維構造，測線中部和大號測點附近地下介質主要呈二維構造。再對214--270各測點主軸方位曲線分析，可知所有測點的主軸方位均為-10°～10°，據(jù)此推斷地下介質的主軸方位與測線方向大體相同，以上推斷符合并豐富了已知地質資料。

表2 張量CSAMT工作參數(shù)設置

圖5 實部Mohr圓圖示Fig.5 Real part Mohr circle

a.視電阻率擬斷面圖 ；b.阻抗相位擬斷面圖；c.地質剖面草圖。圖6 張量CSAMT成果圖Fig.6 Result of tensor CSAMT

圖7 二維反演電阻率與鉆孔資料對比效果圖Fig.7 Contrast effect of two dimensional inversion resistivity and borehole data

2.4 數(shù)據(jù)反演與綜合解釋

從上述兩種方法的定性分析，可初步了解∑19號雜巖體及附近地質體的電性分布及構造特征，對研究區(qū)的整體把握是有益的。為進一步挖掘有用信息，對數(shù)據(jù)進行了二維反演，并與已知的鉆孔資料對比分析(圖7)。據(jù)鉆孔資料分析得知：該雜巖體為向北西向側伏的基性--超基性巖體，巖體在側伏向上巖相分帶明顯，從上到下依次為輝長巖、輝石巖、橄輝巖和橄欖巖相；礦區(qū)內(nèi)已開采銅鎳礦與該礦特征十分相似，見上下兩個鎳礦體(圖7)，同時鎳礦體有變厚變富的趨勢。據(jù)反演結果可知：雜巖體與圍巖界限明顯，地表附近輝長巖體成薄層狀；242測點處的低阻異常與上部銅鎳礦體對應良好，可圈定礦體具體分布范圍；下部銅鎳礦體位于片麻狀花崗巖與橄欖巖、輝橄巖接觸帶部位，據(jù)此可推測，低阻和高阻過渡帶也是不可忽視的找礦有利地段。

3 結論

1)本文開展的張量CSAMT數(shù)據(jù)處理技術初步研究與應用示范，對采集的數(shù)據(jù)進行預處理，得到測量軸的阻抗張量元素，進而從坐標旋轉和Mohr圓兩個方面展開分析，完成了相關公式的推導并采用Fortran語言編程實現(xiàn)。

2)基于本文提出的張量CSAMT數(shù)據(jù)處理技術，首次使用我國自主研發(fā)的DEM電磁探測系統(tǒng)開展張量CSAMT示范研究，對國產(chǎn)儀器性能的完善和應用推廣起到一定的積極作用。

3)Mohr圓分析法和傳統(tǒng)坐標旋轉法在新疆某礦區(qū)的應用效果顯著：測區(qū)內(nèi)各項異性程度于淺部較大，而深部較小，這體現(xiàn)了近地表巖體、特別是雜巖體上部呈全巖蝕變的地質現(xiàn)象；推測小號測點附近地下介質主要呈三維構造，測線中部和大號測點附近地下介質主要呈二維構造；推斷雜巖體內(nèi)部存在銅鎳礦的可能性極大，該推測已由鉆孔資料證實。

4)對比分析了二維電阻率反演結果與鉆孔等地質資料，對礦區(qū)建立相應的地質地球物理找礦模型是有益的。

5)本文僅對張量CSAMT數(shù)據(jù)處理技術進行了初步研究和應用示范，今后在大數(shù)據(jù)和多參量的定量反演和成果解釋方面將做更多的嘗試和研究工作。

本文實測數(shù)據(jù)由中國地質科學院地球物理地球化學勘查研究所電磁綜合研究室陣列電磁法項目組采集并提供,在此表示誠摯的謝意！

[1] 王家映.我國大地電磁測深研究新進展[J].地球物理學報,1997,40(增刊1):206-216. Wang Jiaying. New Development of Magnetotelluric Sounding in China[J]. Chinese Journal of Geophysics, 1997,40(Sup.1):206-216.

[2] 魏文博.我國大地電磁測深新進展及瞻望[J].地球物理學進展,2002,17(2):245-254. Wei Wenbo. New Advance and Prospect of Magnetotelluric Sounding (MT) in China[J]. Progress in Geophysics, 2002,17(2):245-254.

[3] 嚴家斌.大地電磁信號處理理論及方法研究[D].長沙:中南大學,2003. Yan Jiabin. The Study on Theory and Method of Magnetotelluric Signal Processing[D]. Changsha: Central South University,2003.

[4] 吳璐萍,石昆法.可控源音頻大地電磁法在地下水勘查中的應用研究[J].地球物理學報,1996,39(5):62-69. Wu Luping, Shi Kunfa. Application of CSAMT to the Search for Groundwater[J]. Chinese Journal of Geophysics, 1996,39(5):62-69.

[5] 林品榮,趙子言.分布式被動源電磁法系統(tǒng)及其應用[J]. 地震地質，2001，23(2):138-142. Lin Pinrong, Zhao Ziyan. A Brief Description of Distributive Passive Electromagnetic System and Its Application[J]. Seismology and Geology, 2001, 23(2): 138-142.

[6] 張賽珍,王慶乙,羅延鐘. 中國電法勘探發(fā)展概況[J]. 地球物理學報，1994 (增刊1):408-424. Zhang Saizhen, Wang Qingyi, Luo Yanzhong. An Overview on the Development of the Electrical Prospecting Method in China[J]. Chinese Journal of Geophysics, 1994 (Sup.1): 408-424.

[7] 李勇,吳小平,林品榮. 大地電磁測深三維正演模擬及異常特征研究[J]. 地球物理學進展，2012，27 (6):2452-2463. Li Yong, Wu Xiaoping, Lin Pinrong. A Study on Three-Dimension Forward Modeling and Anomaly Features of Magnetotelluric Sounding[J]. Progress in Geophysics, 2012, 27(6):2452-2463.

[8] 李勇，林品榮, 李桐林, 等. 基于異常復電位2.5維CR有限元數(shù)值模擬[J]. 吉林大學學報：地球科學版, 2011,41(5): 1596-1604. Li Yong, Lin Pinrong, Li Tonglin, et al. Finite Element Method for Solving Anomalous Complex Potential of 2.5-D Complex Resistivity[J]. Journal of Jilin University: Earth Science Edition, 2011, 41(5):1596-1604.

[9] 郝興中, 楊毅恒, 李英平，等. 綜合找礦方法在覆蓋區(qū)的應用：以山東省單縣大劉莊鐵礦勘查為例[J]. 吉林大學學報：地球科學版, 2013,43(2): 641-648. Hao Xingzhong, Yang Yiheng, Li Yingping, et al. Application of Integrated Ore Prognosis Method in Covered Area: A Case Study of Shanxian Iron Deposit in Shandong Province, China[J]. Journal of Jilin University: Earth Science Edition, 2013, 43(2):641-648.

[10] 張文秀, 周逢道, 林君, 等. 分布式電磁探測系統(tǒng)在深部地下水資源勘查中的應用[J]. 吉林大學學報：地球科學版, 2012,42(4): 1207-1213. Zhang Wenxiu, Zhou Fengdao, Lin Jun, et al. Application of Distributed Electromagnetic System in Deep Groundwater Prospecting[J]. Journal of Jilin University: Earth Science Edition, 2012, 42(4):1207-1213.

[11] 于向前,趙義平,王明新，等. 音頻大地電磁法與核磁共振法結合劃分含水層的試驗[J]. 吉林大學學報：地球科學版, 2014,44(1): 350-358. Yu Xiangqian, Zhao Yiping, Wang Mingxin, et al. Combination of Audio Magnetotelluric and Nuclear Magnetic Resonace Used to Aquifer Division[J]. Journal of Jilin University: Earth Science Edition, 2014, 44(1): 350-358.

[12] 湯井田, 周聰, 張林成. CSAMT電場y方向視電阻率的定義及研究[J]. 吉林大學學報：地球科學版, 2011,41(2): 552-558. Tang Jingtian, Zhou Cong, Zhang Lincheng. A New Apparent Resistivity of CSAMT Defined by Electric Fieldy-Direction[J]. Journal of Jilin University: Earth Science Edition, 2011, 41(2): 552-558.

[13] Li X. Controlled-Source Tensor Magnetotellurics and Its Application to the Detection of Azimuthal Anisotropy in the Earth’s Crust[D]. Uppsala: Uppsala University, 1994.

[14] Li X, Pedersen L B. Controlled Source Tensor Magnetotellurics[J]. Geophysics, 1991, 56：1456-1461.

[15] Li X, Pedersen. The Electromagnetic Response of an Azimuthally Anisotropic Half-Space[J]. Geophysics, 1991, 56:1462-1473.

[16] Li X, Pedersen. Controlled-Source Tensor Magnetotelluric Responses of a Layered Earth with Azimuthal Anisotropy[J]. Geophysics J Int, 1992,111:91-103.

[17] Boerner D E, Wright J A. Tensor CSAMT Studies at the Buchans Mine in Central New Foundland[J]. Geophysics, 1993, 38(1):12-17.

[18] Bernhard F. Multi Dipole CSAMT[C]//The 10th China International Geo-Electromagnetic Workshop. Nangchang: Chinese Geophysical Society, 2011:96-99.

[19] 晉光文,孫潔,江釗.大地電磁阻抗張量不變量及其Mohr圓分析[J].地震地質,1995,17(4)：439-445. Jin Guangwen, Sun Jie, Jiang Zhao. Mohr Circles and Magnetotelluric Impedance Tensor Invariant Analysis[J]. Seismology and Geology, 1995,17(4)：439-445.

[20] 晉光文,孫潔,江釗.區(qū)域三維大地電磁阻抗張量畸變效應的Mohr圓分析[J].地震地質,1996,18(4)：382-390. Jin Guangwen, Sun Jie, Jiang Zhao. Mohr Circle Analysis 3D Magnetotelluric Impedance Data Affected by Distortion[J]. Seismology and Geology, 1996,18(4)：382-390.

[21] Lilley F E M. Magnetotelluric Analysis Using Mohr Circle[J]. Geophysics,1993,58(10):1498-1506.

Preliminary Study and Demonstration of Tensor CSAMT Data Processing Technology

Meng Qingkui1, Lin Pinrong2, Li Yong2, Li Jianhua2, Zhu Hongwei1, Li Dang2

1.ChinaAeroGeophysicalSurveyandRemoteSensingCenterforLandandResources,Beijing100083,China2.InstituteofGeophysicalandGeochemicalExploration,ChineseAcademyofGeologicalScience,Langfang065000,Hebei,China

Scalar CASMT is inadequate to the face of great depth and three-dimensional detection; while tensor CSAMT is becoming a research focus. On the basis of previous studies, two kinds of tensor CSAMT processing technologies are proposed; which are Mohr analysis and the traditional coordinate rotation. By innovative use of the integrated electrical instruments developed in China for the first time, we demonstrate the application of tensor CSAMT. The results show that through Mohr analysis and traditional coordinate rotation, the vertical and horizontal structure of underground geological bodies can be qualitatively analyzed; apparent resistivity and impedance phase in principal axis orientation correspond with the targets; and the comprehensive analysis of resistivity inversion and drilling information are profitable for the establishment of geological-geophysical models on the test area.

tensor CSAMT; coordinate rotation; Mohr circle; demonstration application

10.13278/j.cnki.jjuese.201506302.

2015-01-23

國家重大科學儀器設備開發(fā)專項(2011YQ050060)；中國地質調(diào)查局地質調(diào)查項目(12120115039401)

孟慶奎(1987--)，男，助理工程師，主要從事應用地球物理方法研究和數(shù)據(jù)處理解釋工作，E-mail:qingkui_meng@163.com。

10.13278/j.cnki.jjuese.201506302

P631.3

A

孟慶奎，林品榮，李勇，等. 張量CSAMT數(shù)據(jù)處理技術初步研究與示范應用.吉林大學學報:地球科學版,2015,45(6):1846-1854.

Meng Qingkui, Lin Pinrong, Li Yong, et al. Preliminary Study and Demonstration of Tensor CSAMT Data Processing Technolog.Journal of Jilin University:Earth Science Edition,2015,45(6):1846-1854.doi:10.13278/j.cnki.jjuese.201506302.