基于大比例尺物探技術的類比求同法對隱伏礦床(體)的預測與評價研究

周 寧， 雷 科

(四川省地質礦產勘查開發局 物探隊，成都 610072)

隱伏礦床(體)的預測在成礦預測中難度相對較大，一般通過在礦區范圍內開展大比例尺成礦預測研究來實現。在相似的成礦地質環境具有相似的成礦地質產物的前提下，通過對已知礦床成礦地質條件、成礦規律及成礦標志認識的基礎上，類比分析，進行深部以及外圍隱伏礦床(體)的預測與評價。

就目前而言，成因完全相同的礦產幾乎沒有，即使礦化類型相同的兩個礦床，其特征也存在差異。因此，單純依靠已知礦床綜合模型對未知區進行找礦預測，其預測結果會有所偏差。另外，已知礦區的勘查、研究程度遠遠高于未知區，工作程度的不均衡性，同時已知區與未知區由于物性密度、磁性以及電性參數的不均一性，都會使兩個區之間資料難以逐一對比。已知礦產模型的認識(包括成礦模式和成礦識別標志)和建立是預測與評價的基礎，其又不完全等同于找礦的預測模型，兩者之間聯系和差異共存。針對隱伏礦床(體)的預測與評價，就預測模型的建立而言，其內容既與已知礦床模型相聯系(即能客觀反映礦床(體)的客觀存在)，又能通過其他某一項資料進行比較分析，提取類似標志和因素。

因此作者針對工程實際，在某地區一大型火山巖成因塊狀硫化物礦床上，通過已知礦產模型的認識分析，結合已知成礦地質環境，利用大比例尺重力、磁法以及電法等綜合物探方法對礦區外圍進行了類比求同分析，對外圍隱伏礦床(體)進行了預測與遠景評價。

1 物探方法簡介

1.1 重力勘探

由于地下巖(礦)石密度分布不均勻而引起的重力變化產生重力異常，它是地質體的剩余質量所產生的引力在重力方向的分量。重力勘探[2]就是以地球的重力場作為被探測物體的引力場，利用重力儀器觀測被探測物體受地球的吸引所產生的重力異常來達到探測目的。

通過高精度重力儀器獲得的觀測數據，經過混合零點位移和固體潮改正以后得到各測點相對于總基點的相對重力值，其包含了因地下密度不均勻的地質體所產生的異常、因各測點周圍地形不同、所處位置的緯度不同等因素的影響，為了單純獲得并比較各測點處重力異常及其大小，必須將各測點的相對重力值按照同一個標準進行相應改正，從而獲得單一地質因素，即地下地質體密度分布不均勻所引起的重力異常值。

重力異常改正包括以下幾個方面：

1)地形改正。測點周圍的地形起伏將會對測點處的重力值產生吸引從而減小其值，故就大比例尺重力勘探而言，對測點需進行近區(0 m～20 m)、中區(20 m～500 m)、遠區(500 m～2 000 m)分別進行改正計算。

2)布格改正。它包括中間層改正和高度改正。由于測點還受到中間層及不同高度的影響，因此在對重力值進行地形改正以后還需要進行布格改正。

3)緯度改正。經過地形、布格改正以后，由于所處緯度不同還存在其正常重力值的不同，這一影響必須去掉。

對觀測值進行地形、布格以及緯度改正以后獲得布格重力異常值，即：

ΔgB=g-g0+δgB+δgT

其中g為測點重力值；g0為正常重力值；δgB為布格改正值；δgT為地形改正值。

1.2 磁法勘探

由于地下巖(礦)石磁性分布不均勻從而引起磁性變化所產生磁異常。磁法勘探以地球的磁場作為被探測物體的磁化場，通過磁力儀其觀測被探測物體被地磁場磁化后所產生的磁異常來達到探測目的。

大比例尺磁測工作中觀測的是磁場總場，因而為了獲得單一因素引起的磁異常，需要對觀測值進行日變改正、磁基點改正、緯度改正和高度改正后才能獲得磁異?！鱐。

1.3 電法勘探

本次研究主要是通過V8多功能電法儀，采用可控源音頻大地電磁測深(CSAMT)方法[3]?？煽卦匆纛l大地電磁測深法CSAMT方法工作原理是基于地下巖(礦)石電性差異性，利用接地偶極按一定的頻率序列(最大范圍2-13Hz～213Hz)向地下發送一次脈沖電流，觀測端(測深點)位于距場源較遠地段(依觀測裝置、目標勘查深度而定)，通過觀測不同發射頻率下電磁場的正交電磁分量及其相位差，計算出不同頻率下的視電率。由于不同頻率的激勵場具有不同的的趨膚深度，因而觀測結果可以反映測點下電阻率隨深度的變化特征。通過對各測深點數據進行匯總、處理及反演計算，則可以得到整個測區內電阻率的空間分布狀態，為進一步地地質解釋提供詳實可靠的深部資料。

CSAMT方法是基于電磁波傳播理論和麥克斯韋方程組[4]導出的水平電偶源在地面上的電場及磁場公式：

式中I為供電電流強度；AB為供電偶極長度；r為場源到接收點之間的距離。

將上面所列式子中的電場Ex與磁場Hy相比，并經過一些簡單運算，就可獲得地下的視電阻率ρs公式：

式中f代表頻率。由視電阻率公式可見，只要在地面上能觀測到兩個正交的水平電磁場(Ex，Hy)就可獲得地下的視電阻率ρs，也稱卡尼亞電阻率。

工作中可以通過人工調整二次場觀測頻率進而采集各觀測點不同頻率下不同方位的電、磁場振幅及相位數據，通過各種復雜的數據處理、反演手段，最終反映出地下電阻率三維分布特征，從而達到測深的目的。

2 礦區成礦模型

成礦作用的發生必然與其周圍的巖石在地質、地球物理、地球化學等方面具有一定的聯系，并以一定的特征表現出來?？偨Y前人資料，對于本次研究區成礦模型[5]，其成礦條件、控礦因素、成礦標志以及地球物理成礦信息標志如下：

2.1 地質特征

1)地層。為中泥盆統中酸性火山碎屑巖層，火山-沉積巖夾碳酸鹽巖建造，是礦區主要的賦礦層位。

2)火山巖。含礦層位為遠火山巖相。巖性主要為火山角礫巖、晶屑凝灰巖、凝灰巖夾集塊巖、角礫凝灰巖、沉凝灰巖，附近有中基性火山巖侵入。

3)構造。處于北北西向、北西向兩組壓性或壓扭性斷裂帶內。

2.2 成礦標志

1)礦化標志。地表有褐鐵礦化、孔雀石化。

2)蝕變標志。其蝕變類型有硅化、綠泥石化、絹云母化、黃鐵礦化(褐鐵礦化)、綠簾石化。

2.3 地球物理標志

1)重力。在礦化蝕變帶上一側有明顯的剩余重力異常，△Gmax=0.3×10-5m/s2。

2)磁法。在重力異常一側及礦體對應位置有明顯的長軸條帶狀的高磁異常，△Tmax≥1 000 nT。

3)電法。在重力異常一側及礦體對應位置表現為明顯的中低電阻率異常。

3 綜合數據處理

3.1 綜合數據處理流程

根據研究區實際地質環境，結合預測和遠景評價的目的任務，有選擇性的對物探數據進行了各種數據整理與處理，其流程如圖1所示。

圖1 數據處理流程圖Fig．1 Data processing flow chart

4 預測區地質-地球物理模型建立

為了完成本次研究，在某地區通過1:10 000重力、磁法掃面工作，對預測區進行了基礎性物探調查，通過符合本區地質環境的數據處理[6]過后，針對提取的有意義異常地段開展了1∶10 000可控源音頻大地電磁測深(CSAMT)工作。根據所取得工作成果，現總結其地質地球物理模型如下：

4.1 地質特征

1)地層。為中泥盆統中酸性火山碎屑巖層，火山-沉積巖夾碳酸鹽巖建造。

2)火山巖。具火山口相和石英角斑質凝灰巖、晶屑凝灰巖，沉凝灰巖?；曰鹕綆r與中酸性火山巖廣泛發育。

3)構造。受近南北向或北西向中—強變形變質構造帶控制。

4.2 成礦標志

1)礦化標志。地表有褐鐵礦化、見鐵帽現象。

2)蝕變標志。圍巖見黃鐵礦化、絹云母化、硅化(次生石英巖化)、綠泥石化蝕變組合。

4.3 地球物理標志

1)重力。在蝕變帶對應一側有明顯的剩余重力高異常，△Gmax≥0.2×10-5m/s2；

2)磁法。在蝕變帶對應位置有明顯的長軸條帶狀的高磁異常，△Tmax=4 500 nT。

3)電法。在蝕變帶對應一側位置表現為明顯的中低電阻率異常。

5 類比分析應用研究

據研究區所取得的成果和建立的地質地球物理模型，與已知區進行了類比分析，其情況如表1及圖2、圖3所示。

表1 類比分析表

圖2 已知礦床(體)特征Fig．2 Known deposits (body) features(a)剩余重力異常；(b)磁異常；(c)礦體位置及地質

經過已知礦床(體)地質-地球物理模型的對比分析可知，研究區所處的地質環境與已知區近于相同，具有相同的礦化和蝕變特征，同樣表現出了高重、高磁、中低電阻率的地球物理特征。因此鑒于預測區與已知礦床(體)的相似特點，推斷預測在異常區范圍內存在一隱伏礦床(體)，成礦可能性很大,后經鉆孔驗證在該異常區內多個鉆孔約在210 m處見礦。

6 結論與建議

類比求同分析方法在對隱伏礦床(體)的預測與評價中具有重要的作用，對隱伏礦床(體)在成礦預測中難度較大的實際情況，根據預測區地質環境和物性差異，通過大比例尺的物探方法技術，進一步剖析預測區與已知區的地質-地球物理特征的相似性，建立兩者之間的類比關系，從而對隱伏礦床(體)進行推斷預測將是以后成礦預測的又一重要指導。

本次研究不僅涉及地質資料的深度解析，同時涵蓋了物探技術[7-8]范疇，因此為了更好地提出預測模型，更加準確地完成類比分析，就需要對地質和物探方法的多解性進行經驗試驗，針對實際情況提取出單一地質因素引起的地質-地球物理模型場，以期對隱伏礦床(體)作出更加合理的預測與評價。

參考文獻：

[1] 李才明，李軍.重磁勘探原理與方法[M]. 北京：科學出版社，2013.

[2] 肖宏躍，雷宛. 地電學教程[M]. 北京：地質出版社，2008.

[3] 王全明，馮京，莊道澤，等. 銅鋅礦床特征與隱伏礦類比求同法預測[M]. 北京：地質出版社，2000.

[4] 李金銘.地電場與電法勘探[M]. 北京：地質出版社，2005.

[5] 王家華. 克利金地址繪圖技術/計算機的模型和算法[M]. 北京：石油工業出版社，1999.

圖3 已知礦床(體)特征Fig．3 Geophysical forecast - geological features(a)剩余重力異常；(b)磁異常；(c)橫穿異常剖面視電阻率異常；(d)異常區地質圖

[6] 姚鳳娘，鄭明華. 礦床學基礎教程[M]. 北京：地質出版社，1980.

[7] 曾懷恩，黃聲享. 基于Kriging方法的空間數據插值研究[J]. 測繪工程，2007，16(5)：5-8，13.

[8] 王金玲，張東明. 空間數據插值算法比較分析[J]. 礦山測量，2010(2)：55-57.