礦區(qū)金屬礦勘探AMT法運用

1、礦區(qū)金屬礦勘探amt法運用音頻大地電磁測深(amt, audiomagnetotelluric)法主要用于金屬礦、地 熱田、工程和埋藏不深的油氣田勘探。20世紀60年代初，strangway等 提出le音頻人地電磁法。在20世紀70年代，strangway等在應用音頻大 地電磁測深法尋找金屬硫化礦床方面做了大量t作，取得了有意義的成果 l_2jo音頻大地電磁測深法是通過觀測由天然平面電磁波信號以確定地下 的電阻率值的方法，具有対低阻體有較好的分辨率的優(yōu)點，高效低成本。 由于頻率較高，對淺部的分辨率較高，更適于資源勘探，是重要的地球物 理探勘手段。缺點是場源不可控制并且信號微弱，易受自然環(huán)境的影

2、響， 尤其是在礦山、城區(qū)附近很難開展工作。下面，筆者利用軟件處理解釋某 礦區(qū)的amt數(shù)據(jù),得到tm極化模式(橫向極化模式)下的二維視電阻率斷 面圖，根據(jù)對比地質剖面和鉆孔資料，結合超基性巖體的礦化規(guī)律，論述 了 amt法在某礦區(qū)金屬礦勘查屮的應用效果。1礦區(qū)地質與地球物理特征11地質特征礦區(qū)地處太平洋板塊向亞洲板塊俯沖，消減形成一系列北東一北北東向深 斷裂以及東西向和北西向3組斷裂構造并存的格子狀斷裂系統(tǒng)。這些格子 狀排列的斷隆、斷陷帶和格子狀斷裂構造系統(tǒng)是該區(qū)主要的構造形式，它 們不僅控制了中生代斷陷盆地，也控制了區(qū)內礦產(chǎn)的分布格局。該礦區(qū)地 下巖體為中元古代超基性巖，山斜方輝右和貴橄欖石組

3、成，巖石強烈鎂鐵 閃石化，與其伴生較多金屬礦物沿其礦物節(jié)理或空隙分布，屬鎂質超基性 巖。地表見超基性巖，呈北東向透鏡狀產(chǎn)出，總體向北西陡傾斜，自北向 南西側伏。巖石主要礦物成分為橄欖石、普通輝石和片狀金云母。橄欖巖 面呈灰黑色，塊狀構造，較強蛇紋石化。1. 2地球物理特征該礦區(qū)超基性巖體具低阻高極化異常特征。在礦區(qū)內共釆集了 3個鉆孔(井 一 1、井一 2、井一 3)巖芯、探槽及地表露頭巖石標本，根據(jù)采集的巖芯 和巖石標本發(fā)現(xiàn)該區(qū)的蝕變大理巖、橄欖巖的密度特別高，平均值達3. oog /cm。以上；角巖、輝綠巖的密度較高，平均密度值達2. 90g/cmo以上; 此類巖石若冇一定的規(guī)模，則能引起

4、明顯的局部重力高異常。在磁性方面, 輝綠巖的磁性最強，磁化率值達2909x4nx10 一(si)；橄欖巖、花崗巖 的磁性較強，磁化率平均值分別為1599x4nx10 一 s(si)和1162x4nx lo-s(si)；若此類巖石埋藏較淺，在地面則能引起峰值不同的局部高磁 異常。2amt原理及數(shù)據(jù)處理2. 1amt方法原理釆用標量的方式，即在測點同時觀測相互垂直的電場(e)和磁場(e)分量， 并按：/o計'算視電阻率，繪出視電阻率曲線，并對它進行解釋口。 式中，p 口是視電阻率；是電磁波的頻率；flo是真空中的磁導率。2. 2反演算法doccam反演occam反演實際是求一個多層地球模型

5、的最光滑解，即在一 定的擬合誤差標準下使模型的粗糙度最小。occam反演以其穩(wěn)定的收斂性 在mt(magnetotelluric,大地電磁測深法)中得到廣泛應用，其穩(wěn)定收斂 性主要來自2個方面：對數(shù)等間隔的垂直剖分及求最光滑模型口。2)非線性共輒梯度法該算法是為了計算mt的2d反演問題的正則化解，使 h標函數(shù)最小以能使錯位數(shù)據(jù)和電阻率的空間二次導數(shù)規(guī)則化4,即校 正錯位數(shù)據(jù)和模型的光滑度。二維非線性共軌梯度法反演廣泛應用于mt 的2d反演中。2. 3數(shù)據(jù)處理1) 數(shù)據(jù)的編輯野外采集工作，由于干擾和觀測誤差的存在，某些頻段的頻 點數(shù)據(jù)有時會出現(xiàn)非正常的跳躍和飛點。根據(jù)相鄰測點的曲線特征，對每 一

6、個測點進行編輯，對測點曲線進行圓滑或丟棄誤差太大的頻點，保留高 質量的頻點數(shù)據(jù)。2) 1d反演根據(jù)實測的視電阻率喘線分層，并賦予每層電阻率值和厚度值， 以此作為初始模型進行occam反演，并對初始模型不斷修改再反演，直到 計算的視電阻率曲線和實測的視電阻率曲線擬合差達到所要求的精度。1d 反演形成2d擬斷面，并能得到視電阻率剖斷面。2)2d反演以ldoccam反演后得到的視電阻率剖斷面作為2d反演的初始模 型剖分網(wǎng)格，并設置反演參數(shù)進行2d非線性共軌梯度反演，最終得出2d 反演視電阻率斷面圖。3成果資料解釋在數(shù)據(jù)處理進行2d反演時，需要設定一些參數(shù)，以約束反演模型和迭代, 比如平滑因子r、門檻

7、誤差、加權函數(shù)a和等。在tm極化模式下取平滑因 子r 1,門檻誤差tmrho=10和tmphase 一 5,加權函數(shù)a 0.3和fl=0, 迭代100次后，均方根誤差rms-2. 123o最后得出2d反演視電阻率斷面圖和2d模型的實測的與計算的視電阻率擬斷面擬合比較圖，如圖1和 圖2所示。圖2中，上半部為實測的2d視電阻率擬斷面圖，下半部為計 算的2d視電阻率擬斷面圖。由圖2可看出，計算的視電阻率值和相應深 度絕大部分與實測視電阻率值和相應深度相符合，擬合效果是令人滿意 的。鉆孔井一 1 井一 2和井一 3的位置分別處丁測點l0 1300> 10一 1150和101050附近，開孔角皆為

8、80。地形最高點在測點lo 一 1350 附近，地質剖面中的鉆孔位置和角度如圖3所示。由圖1中可看出，反演 不足之處是細節(jié)的反映不強，尤其是近地表細節(jié)反映不詳細。測點l0-950 到l0-1400卜大概350m處出現(xiàn)高阻異常，視電阻率等值線梯度變化比較 人。根據(jù)地質及鉆孑l資料，井一 1孔遇角巖，人理巖和角巖夾層，部分 橄欖巖，鉆孔底部遇花崗巖；井一 2孔大部分為橄欖巖。井一 3孔大部分 為輝石橄欖巖和輝綠巖，少量花崗斑巖、角巖、石英正長斑巖和黑云母花 崗巖等。在視電阻率深度斷面圖中，近地表巖層嚴重風化破碎含水，在電 性上都表現(xiàn)為低阻；井一1孔底部遇花崗巖的部位為高阻區(qū)，在圖1屮顯 示為高程450m左右，此區(qū)域以下為高阻體。井一 1孔至井一 3孔之間大 部分為輝石橄欖巖，表2現(xiàn)為低阻，在圖1中的顯示為測點l0-1300至 測點lo 一 1050之間高程400m以上為低阻，與地質剖面圖對應的很好。根據(jù)地質資料，超基性巖體呈北東方向展布，總體向北西陡傾斜，自北東 向南西側伏，超基性巖主要集屮分布在剖面的屮部。圖3和圖1比