MN長度誤差及方向偏差對EH4視電阻率的影響

羅安華，李文堯

MN長度誤差及方向偏差對EH4視電阻率的影響

羅安華，李文堯

（昆明理工大學 國土資源工程學院，昆明 650093）

在EH4工作中，MN誤差會導致視電阻率誤差。有學者曾推導出MN長度誤差引起的視電阻率相對誤差公式，但公式的推導過程與實際工作有一定的偏差。MN方向偏差、MN長度誤差及方向偏差同時存在對EH4視電阻率的影響，卻未見報道。因此，這里對MN長度誤差、MN方向偏差、MN長度誤差及方向偏差同時存在引起的EH4視電阻率相對誤差進行理論推導和數值計算。研究表明，EH4視電阻率誤差隨著MN誤差的增大而增大；當存在MN長度誤差時，MN長度誤差為負的視電阻率相對誤差比MN長度誤差為正時大。

EH4；高頻大地電磁法；視電阻率誤差；MN長度誤差；MN方向偏差

0 前言

EH4是以巖（礦）石間的電性差異為基礎，通過觀測時間域電磁場信號，經傅里葉變換，轉換為頻率域內電磁場信號，來查明地質構造或有用礦產的地球物理勘探方法。由于EH4①具有能同時使用人工電磁場和天然電磁場兩種場源；②具有有源電磁法的穩定性及無源電磁法的節能和輕便；③同時接收和分析X、Y兩個方向的電場和磁場，反演X－Y電導率張量剖面，對判斷二維構造特別有利；④觀測時間短，能在短時間內完成較大深度的測量，實時數據處理和顯示，資料解釋簡捷，圖像直觀等優點。因此EH4被廣泛應用于地殼結構研究［1～3］，地下水調查［4］，環境調查［5］，金屬礦勘探［6、7］，地熱勘探［8、9］以及工程研究［10、11］等方面。

李汝傳［12］對MN長度誤差產生的卡尼亞視電阻率誤差進行理論推導和數值計算研究，研究結果表明，卡尼亞視電阻率相對誤差為MN長度誤差的二倍。周劍［6］、曾玉嬌［13］對MN長度誤差產生的卡尼亞視電阻率誤差進行總結，并計算出多組由MN長度誤差引起的卡尼亞視電阻率相對誤差。劉春明等［14］對EH4野外工作方法進行研究，討論影響EH4勘探效果的各類因素，但并未推導MN長度誤差引起的EH4視電阻率相對誤差的公式。

由于MN長度誤差引起的視電阻率相對誤差公式的推導過程與實際工作有一定的偏差，MN方向偏差、MN長度誤差及方向偏差同時存在對EH4視電阻率影響的研究未見報道。因此，作者在本文對MN長度誤差、MN方向偏差、MN長度誤差及方向偏差同時存在所引起的EH4視電阻率誤差進行了研究。

1 EH4工作的基本原理［14～18］

EH4通過測量MN方向的電場E和與之垂直的磁場H計算卡尼亞視電阻率。大地電磁測深理論以Maxwell方程組為核心，將天然電場和與之垂直磁場分量的比值定義為地面波阻抗：

當交變電磁場以波的形式在地下介質中傳播時，由于電磁感應的作用，地面電磁場的觀測值包含有地下介質電阻率分布的信息，卡尼亞視電阻率計算公式為：

由式（2）計算所得到的電阻率為視電阻率，f為電磁波頻率。

由于電磁場的集膚效應，不同周期的電磁場信號具有不同的趨膚深度或穿透深度，電磁波的趨膚深度表達式為：

式中 ω為角頻率；μ為磁導率；σ為電導率。

由于趨膚深度并不代表有效探測深度，Bostic定義的有效探測深度經驗計算公式為：

電磁波的趨膚深度取決于大地電阻率和電磁波頻率，在大地電阻率結構穩定時，通過改變電磁波頻率可以對大地進行連續的垂直測深。

2 MN長度誤差對EH4視電阻率的影響

在野外測量工作中，實際測量的MN長度與設計的MN長度難免存在誤差。

［12］，MN長度誤差產生的卡尼亞視電阻率相對誤差公式為：

在實際工作中，通過測量MN間的電位差V，并用MN對V歸一換算成電場值：

將式（6）代入式（5）得：

式（5）、式（7）中Vc為電極距MN無誤差時測量的電位差；Vr為電極距MN存在長度誤差ΔMN時測量的電位差；ΔV＝Vr－Vc；MN為電極距長度；ΔMN為電極距長度誤差。

在實際工作中，無論MN是否存在長度誤差，測量和用于計算視電阻率的電位差只有一個，視電阻率誤差是由MN長度誤差引起，而非電位差誤差引起。公式（7）用電位差誤差推導所得，推導過程與實際工作有一定的偏差。因此，應重新推導視電阻率相對誤差公式。

設實際測量電極距長度為MN，電極距長度誤差為ΔMN，則設計電極距長度為MN＋ΔMN，實際測量MN間的電位差為V，將式（6）代入式（2）后可得MN間的實際視電阻率表達式（8）。

由于未發現電極距長度誤差ΔMN，因此，誤將實測電極距長度MN替換為設計電極距長度MN＋ΔMN，將式（6）代入式（2）可得MN間含電極距長度誤差的視電阻率表達式（9）。

視電阻率相對誤差表達式為：

將式（8）和式（9）代入式（10），得出視電阻率相對誤差表達式為：

當MN誤差相同時，由式（12）和式（7）計算視電阻率相對誤差見下頁表1。

計算結果表明：①視電阻率誤差隨著MN長度誤差的增大而增大；②式（12）計算結果顯示，MN長度誤差為負的視電阻率相對誤差較MN長度誤差為正時大；③式（12）與式（7）計算的結果有差別。

3 MN方向偏差對EH4視電阻率的影響

MN方向偏差示意圖見下頁圖1。

在野外工作的過程中，并不能保證MN的測量方位與設計測線方位（Ex代表測線方位）完全一致，設MN長度測量準確，而MN測量方位與設計測線方位存在一夾角θ時，MN的方位變為M1N1，所測量的M1N1之間的電位差為MN0間的電位差，即M1N1間的電位差比MN間的電位差小。根據式（2）可知，MN方向偏差將對EH4視電阻率產生影響。

在實際工作中，設MN長度測量準確，MN與測線方向有一個夾角θ，測量M1N1間電位差為V，將式（6）代入式（2），可得實際視電阻率表達式：

由于未發現MN存在方向偏角θ。因此，誤將實際電極距長度MN cosθ替換為設計電極距長度MN，將式（6）代入式（2），可得含MN方向偏差的視電阻率表達式：

將式（13）和式（14）代入式（10），得視電阻率相對誤差表達式為：

當MN方向偏角為θ時，EH4視電阻率誤差計算結果見表2。

計算結果表明：當MN長度測量準確時，MN方向偏差θ越大，所產生的視電阻率誤差越大。

4 MN長度誤差及方向偏差對EH4視電阻率的影響

在實際工作中，MN長度誤差和MN方向偏差是同時存在的。當MN同時存在長度誤差ΔMN和方向偏差θ時，由式（11）和式（15）可得視電阻率相對誤差表達式：

當MN誤差為ΔMN，方向偏差為θ時，視電阻率相對誤差計算結果見表3。

計算結果表明：①EH4視電阻率誤差隨著MN長度誤差的增大而增大；②當MN長度誤差為正時，視電阻率相對誤差隨著MN方向偏差的增大而增大；當MN長度誤差為負時，視電阻率相對誤差隨著MN方向偏差的增大而減小；③MN長度誤差為負的視電阻率相對誤差比MN長度誤差為正時大。

5 實測數據反應計算與公式（12）理論計算結果對比

結合實測數據驗證公式（12），本次EH4工作的MN長度為80 m，野外數據采集所用儀器是美國Geometrics公司和EMI公司聯合研制的雙源型電磁系統，儀器型號為Stratagem，采用低頻探頭測量，反演軟件為IMAGEM程序。

在DOS系統中編輯“＠L”文件中的MN長度值后，用IMAGEM程序計算編輯數據后的視電阻率值。

實測數據反演計算和式（12）理論計算視電阻率相對誤差結果見表4。

計算結果表明：實測數據反演計算與公式（11）計算視電阻率相對誤差一致。

6 結論

（1）作者在本文推導了MN長度誤差產生的視電阻率相對誤差公式（12）、MN方向偏差產生的視電阻率相對誤差公式（15），以及MN長度誤差和方向偏差產生的視電阻率相對誤差公式（16）。并結合實測數據，驗證了作者在本文中推導的公式（12）是正確的。

（2）由于EH4視電阻率誤差隨著MN長度誤差的增大而增大，所以在工作中，應盡量減小MN長度誤差。

（3）由于視電阻率誤差隨著MN方向偏差的增大而增大，而MN方向偏差引起的EH4視電阻率誤差比MN長度誤差引起的視電阻率誤差小，所以在工作中，應減小MN方向偏差。

（4）MN長度誤差及方向偏差同時存在對EH4視電阻率誤差的影響規律為：①EH4視電阻率誤差隨著MN長度誤差的增大而增大；②當MN長度誤差為正時，視電阻率相對誤差隨著MN方向偏差的增大而增大，當MN長度誤差為負時，視電阻率相對誤差隨著MN方向偏差的增大而減小；③MN長度誤差為負的視電阻率相對誤差比MN長度誤差為正時大。

（5）減小MN誤差，可以減小EH4視電阻率誤差，也可以減少異常推斷解釋的多解性，使工作成果更可靠。

致謝：

感謝昆明鑫地地質勘探有限公司提供EH4野外測量數據；感謝楊瑨卿在野外數據采集過程中提供的幫助。

參考文獻：

［1］JORDING A，FERRARI L，ARZATE J．Crustal variations and terrane boundaries in southern Mexico as imaged by magnetotelluric transfer functions［J］．TECTONOPHYSICS，2000（327）：1．

［2］張作倫，曾慶棟，于昌明，等．利用EH4、GDP－32Ⅱ預測內蒙新民煤田煤系地層空間特征［J］．地球物理學進展，2009，1（24）：303．

［3］李艷麗，林春，于建國，等．EH4電磁成像系統在杭州灣地區晚第四紀地層中的應用［J］．地質評論，2007，3（53）：413．

［4］王玉和，王厚臣，程久龍，等．EH4電磁方法在斷層定位及導水性探測中的應用［J］．山東科技大學學報：自然科學版，2010，4（29）：35．

［5］周劍，白宜誠．EH4電導率成像系統在水電隧洞工程地質災害預報中的應用［J］．礦產與地質，2006，4－5（20）：560．

［6］周劍．高頻大地電磁測深法在金屬礦產探測中的應用研究［D］．中南大學碩士學位論文．2006．

［7］王沖，董平，孫斌，等．EH4電磁測深在江西城門山礦區深部及外圍找礦中的應用［J］．地質與勘探，2009，6（45）：742．

［8］李朝林，裴忠，金秀芹，等．大地電磁測深應用于地熱勘探［J］．中國煤炭地質，2008，9（20）：68．

［9］謝可明，宋喜林．高頻電磁測深法在地熱勘探中的初步應用［J］．物探裝備，2004，2（14）：118．

［10］趙智能，劉建，邢保國．EH4電磁成像系統在高速公路隧道工程勘察中的應用研究［J］．工程地球物理學報，2010，3（7）：323．

［11］祝順義，陳慶凱．EH4電磁成像系統在水庫大壩勘查滲漏的應用［J］．有色礦冶，2004，2（20）：10．

［12］李汝傳．MN長度誤差及接地電阻對卡尼亞電阻率的影響［J］．地質與勘探，2004，增刊（40）：188．

［13］曾玉嬌．EH4在云南金牛廠鉛鋅礦的應用研究［D］．昆明理工大學碩士學位論文．2009．

［14］劉春明，柳建新，李夕兵．EH4野外工作方法的研究與應用［A］．中國地球物理學會．第8屆中國國際地球電磁學討論會論文集［C］．湖北：中國地球物理學會地球電磁學專業委員會：2007．

［15］GAMBLE TM，GOUBAU WM，CLARKE J．Error analysis for remote reference agnetotelluric［J］．Geophysics，1979，44（5）：959．

［16］SIMS WE，BOSTICK FX，SMITH HW．The estimation of magnetotelluric impedance tensor elements from measured data［J］．Geophysics，1971，36（5）：538．

［17］KAO DW，RANKIN D．Enhancement of signal－to－noise ratio in magnetotelluric data［J］．Geophysics，1977，42（1）：103．

［18］石英俊，劉國棟，吳廣耀，等．大地電磁測深法教程［M］．北京：地震出版社，1985．

P 631．3＋25

A

10．3969／j．issn．1001－1749．2012．04．05

2011－10－21 改回日期：2012－03－22

1001—1749（2012）04—0397—05

羅安華（1984－），男，碩士，昆明理工大學地球探測與信息技術專業。