二維小波多尺度分析技術在銅山口深部找礦中的應用

梁學堂，張木林

(湖北省地球物理勘察技術研究院，湖北武漢 430056)

二維小波多尺度分析技術在銅山口深部找礦中的應用

梁學堂，張木林

(湖北省地球物理勘察技術研究院，湖北武漢 430056)

小波變換是重力異常分解的有效工具，利用小波多尺度分析方法，可以將重力異常分解到不同尺度空間中。不同尺度的重力異常反映了不同地質體的規模和埋深。作為一種新的有效的位場分離途徑，小波多尺度分析方法為重力資料解釋和研究地殼結構提供了新的思路，在目前廣泛開展的深部找礦工作中將發揮重要的作用。

重力勘探;布格重力異常;小波多尺度分析;對數功率譜;深部找礦

0 引言

重力勘探作為一種經典的物探找礦手段，其良好的抗干擾能力和深穿透特點，在目前廣泛開展的深部找礦中，越來越受到重視。特別是近年來應用高精度自動讀數重力儀(CG-5)，結合高精度GPS測地技術，在礦調項目中重力勘探觀測精度從≥±0.100×10－5m/s2提高到≤ ±0.025 ×10－5m/s2，精度提高了約 4 倍;布格異常總精度也從原來的≥±0.250×10－5m/s2提高到≤±0.100 ×10－5m/s2。因此，隨著重力異常精度的提高，為我們應用重力勘探技術發現更深部的礦床提供了技術上的支撐。

重力勘探同其它大多數物探方法一樣，也有他本身固有的缺點。一般來說，該方法對橫向疊加的異常信號分辯率高，對縱向疊加的異常信號，因體積效應和地質體埋深的增加，往往被淺部的信息所覆蓋而“湮滅”。比如應用重力勘探尋找盆地構造、火山機構、大型巖體及圈定線性構造等效果特別好;但對深部目標地質體(隱伏小巖體或礦致異常等)產生的弱信號異常，因疊加在淺部或較大地質體產生的大幅值、大范圍異常之上，反應則不甚明顯。另外，在沒有深部物性資料的情況之下，引起局部異常的地質體的具體埋深也不易被定量計算。因此，能否從“強大”的背景場中發現來自地下深部弱小的、有用的異常信息，能否在沒有深部物性資料的情況下對目標地質體的埋深實現定量或半定量的計算，將直接影響到應用重力勘探技術找礦的效果［1］。

提取地下深部目標地質體產生的弱信號異常信息，是研究重力勘探技術的重要課題，同時也是目前開展深部找礦工作的要求。從理論上來講，重力場是位場，不論是在空間域或頻率域內，來自不同尺度空間的疊加異常信息應該是可分離、可轉換的，這樣就為應用數學手段分離礦致或有用的弱信號異常信息成為可能。

說是“雨幕”一點也不錯，那陣雨就像電影散場時拉起來的厚重黑幕，整齊地拉成一列，雨水則踏著軍人的正步，齊聲踩過田原，還呼喊著雄壯威武的口令。

1 二維小波多尺度分析技術簡介

小波多尺度分析又稱多分辨分析，它把一個信號分解為逼近部分和細節部分，表示為是逼近部分，Di細節部分，圖1為三層多尺度分析結構圖，其中，S是信號，A1、A2、A3是逼近部分，D1、D2、D3是細節部分。

圖1 三層多尺度分析結構圖Fig.1 Analysis structure graph of multiple scale

根據小波多分辨分析的原理以及小波細節的微分特征，可實現對位場的多尺度分解。地表局部地質體密度不均勻或儀器觀測誤差、地形校正誤差產生的重力異常，是高頻成分;深部的或大的地質體產生的異常是相對低頻成分;深度越大，頻率越低。不同階次的小波細節，對應不同尺度空間的重力異常。一階細節反映淺部的地質體或干擾成分;二階細節反映的是相對更深一些的地質體分布特征;以此類推，階次越高，反映的地質體埋深越大。

2 應用小波對數功率譜計算不同尺度空間地質體的埋深

分解到幾階要根據異常的特征和地質情況來決定，解釋時要賦于小波逼近部分和各階的細節明確的地質意義。小波變換是重力異常分解的有效工具，利用小波多尺度分析方法，可以將重力異常分解到不同尺度空間中，不同尺度的重力異常反映了不同地質體的規模和埋深。作為一種新而有效的位場分離途徑，小波多尺度分析方法為重力資料解釋和研究地殼提供了新的思路，在國內外得到了廣泛的應用。

圖2 對數功率譜曲線圖Fig.2 Curve of logarithmic power spectrum

眾所周知，銅山口地區找礦工作程度很高，淺部小巖體的分布情況已基本查明，但深部隱伏巖體所產生的剩余重力異常因其重力場隨埋深增加而衰減，地表所觀測到的重力異常值一般不會超過 0.500×10－5m/s2。如此弱小的重力異常值，疊加在數十毫伽(×10－5m/s2)的背景場之上，對我們研究目標地質體來說是比較困難的。因此，必需對布格重力異常采取科學、有效的手段進行分離，才能提取我們研究中酸性小隱伏巖體產狀特征所需要的有用信息。

根據小波細節值，由相關軟件計算對數功率譜數據，人工根據對數功率譜數據，在坐標紙上算出深源地質體的似深度H，淺源地質體的似深度h以及縱軸截距之比B/b，計算公式為：

式中：lnE(ω2)、lnE(ω1)分別是兩個不同頻段 ω2、ω1的對數功率譜值;ω是角頻率，ω=nΔω;Δω為基頻，把測區的長邊作為基波波長，如測點方向的邊長一些，則把點數×點距作為基波波長，其倒數即為基頻，反之則用測線方向作為基波波長，Δω=2π/(總點數 ×點距)，淺源地質體似深度的計算與上式一樣。

在縱軸上讀取深源場與淺源場縱軸截距lnB2、lnb2，由下式算出B/b：

圖4～圖7是對銅山口重力異常所進行細節分析及對數功率譜計算結果。因測區范圍只有幾十公里，只對重力異常進行了四階分解，研究深度也只局限在地下1 500 m深度左右。從各階小波分析結果可以看出，一階細節所反映的是地表淺部幾十米之內的信息，重力異常較為零碎，它主要包含了淺部地質的構造信息，同時也包含有地形校正、重力測量及GPS測量等各種方法技術所產生的誤差，其地質意義不大;二階細節反映的是地下200～250 m深度左右的地質構造信息，與剩余重力異常、垂向二階導數異常的特征相似，重力低異常主要反映的是酸性巖漿巖體的位置及空間形態;重力高異常反映的是沉積巖地層的構造信息;三階細節反映的是地下650～700 m深度左右的地質構造信息，淺部的、小型的巖漿巖體的影響已經消除，突出了大型的、深部的巖體的構造特征，特別是三大巖體的邊界及區內西北部沿弧形侵入的花崗閃長斑巖的空間分布特征非常醒目;四階細節反映的是地下1 500～1 600 m深度左右的地質構造信息，從四階細節圖中可以看出，區內大致可分為巖漿巖和沉積巖二個區域，姜橋巖體和殷祖巖體的根部在此深度已連成一體，而測區西北部的花崗閃長斑巖巖體已無反映。據此可以判定，花崗閃長斑巖巖體為“無根”，為淺層侵入體，其侵入深度不會超過1 500 m。

The picture of the discharge was captured using a PCO Dimax HD high-speed camera with an AF NIKKOR 60 mm f/2.8D lens.The frame rate of the camera was set to 33.3 kHz.The exposure time was set to 1.4 ls.

3 二維小波多尺度分析技術在銅山口深部找礦中的應用

銅山口地區，屬于大冶礦田范圍，為長江中下游重要的成礦區帶。該地區礦床主要類型有矽卡巖型、斑巖—矽卡巖復合型銅鉬礦等。成礦地質模式主要與區內中酸性巖體、特別是與燕山期酸性花崗閃長斑巖的侵入密切相關。花崗閃長斑巖的密度值較低，一般為2.50 g/cm3左右，與成礦圍巖二疊系、三疊系地層(地層平均密度約2.70 g/cm3左右)有一定的密度差，這為在本地區開展大比例尺的重力調查工作提供了地球物理前提。因此，圍繞銅山口地區所開展的礦產調查工作中，重力測量的主要目的任務是發現與成礦密切相關的小型花崗閃長斑巖的隱伏位置，并查明其產狀特征［2］。

圖2是一個運用對數功率譜計算場源似深度的一個實例。從圖中可以看出，對應低頻段部分，對數功率譜曲線的斜率越大(指絕對值)，曲線在縱坐標軸上的截距越大，對應的場源似深度越大;中高頻段對數功率譜的曲線特征也是同樣的規律，對應高頻段，曲線發散，主要反映了地表地質體局部不均勻或儀器觀測誤差等造成的干擾成分。具體計算方法及公式如下。

3.1 銅山口地區布格重力異常特征

圖8是根據小波細節分析結果，綜合測區垂向二階導數異常圖得出的地表下500 m深度左右的巖漿巖分布圖。圖中蘭色部分為較低密度的中酸性巖漿巖分布區域，紅色部分為灰巖、大理巖、白云巖等密度相對較大的巖性分布區域。

在育苗之前要科學配制營養土，保證營養土中氮、磷、鉀元素供給充足，為幼苗根系生長奠定堅實基礎。幼苗出生之后，由于根系水分和養分吸收能力較差，可以進行葉面施肥。苗床灌溉少而勤，保持好大棚內光照。

圖3 銅山口測區布格重力異常圖Fig.3 Bouguer gravity anomaly map of Tongshankou survey area

3.2 銅山口地區重力異常小波各階細節特征

資源描述框架 RDF（resource description framework）是由萬維網協會W3C提出的一個語義框架[1]，被廣泛應用在描述語義網[2]中的各類海量數據，可以用三元組（主語、謂語、賓語）的形式來描述語義網上的任何數據。

3.3 銅山口地區隱伏中酸性巖體分布及深部構造特征

從圖3看出，布格異常值總體上呈西高東低特征。在測區銅山鎮周圍，布格異常等值線近北東向弧形排列，弧形中心在測區之外，大概與靈鄉巖體相對應;在測區中部，布格重力異常等值線由北東向弧形排列轉向近南北向的“∫∫∫”形狀;在測區東半部，布格重力異常呈半圓弧排列，圓弧中心位于姜橋鎮南部。另外不難看出，局部重力異常疊加在背景場之上，只是使布格異常等值線扭曲變形，具體細節并不清晰。因此，單從布格異常等值線圖來研究中酸性小巖體的分布特征是困難的。

中酸性巖漿巖(特別是花崗閃長斑巖)是銅山口地區重要的成礦物質來源，與該地區矽卡巖型及斑巖型銅鉬礦床的形成密切相關。查清區內小型中酸性巖體的隱伏位置、深度及其它產狀特征，可以為在本地區開展深部找礦工作提供靶區。從圖8所反映出的巖體分布的信息，一部分目前已得到驗證，與以往所取得的地質認識也是相同的，另一部分尚未被認識，為新的發現。比如：區內已知姜橋巖體與殷祖巖體的分布位置與邊界，在圖8上得到了充分顯現;在測區西北地區銅山口一帶，除已知小型巖體外，本次發現還存在數量可觀的小型隱伏花崗閃長斑巖巖體，這對以后在本地區開展深部找礦工作起著重要的指導意義。

圖4 重力異常一階細節及功率譜曲線Fig.4 Power spectrum curve of first detail of gravity anomaly

圖5 重力異常二階細節及功率譜曲線Fig.5 Power spectrum curve of second detail of gravity anomaly

圖6 重力異常三階細節及功率譜曲線Fig.6 Power spectrum curve of third detail of gravity anomaly

圖7 重力異常四階細節及功率譜曲線Fig.7 Power spectrum curve of fourth detail of gravity anomaly

圖8 測區巖漿巖分布圖Fig.8 Distribution map of magmatic rock in survey area

圖9是根據小波對角線方向細節進行斷裂分析(奇性檢測)得出的深部斷裂構造影像圖。從圖中可以看出，測區深部(隱伏)大的斷裂構造有2條，一條經過銅山口礦區東側近南北向展布;另一條經過胡家灣—大角山近北東東向延伸。

圖9 測區深部斷裂構造圖Fig.9 Structural map of deep fault in survey area

本次所發現的兩條斷裂，為以前地質未曾認知的隱伏深部大斷裂。特別是在銅山口礦區東側近南北向斷裂的發現，可能會對該地區礦床形成的機制產生新的認識與看法。該斷裂結合整個鄂東南地區區域重力異常圖來看，為一條較大型重力梯級帶，北起黃石，經大冶市、銅山口鎮向南延伸至黃沙鋪鎮，測區內展布只是其中一小段。該斷裂以西，為金山店鎮、太和鎮等地形成的片狀重力高值區，并分布有陳子山、大廣山、銅山口等多處礦區;斷裂帶以東，為以陽新巖體、白沙鋪巖體、殷祖巖體形成的重力低值區。由此可見，該斷裂不但規模宏大，對鄂東南地區巖漿巖的分布有控制意義，是重要的巖漿通道，同時也具有控礦特征，是重要的含礦帶。

如果命令結果如上文所示，表明已經安裝好了BIND套件。否則需要執行安裝命令安裝bind和bindutils兩個包，然后將DNS服務設置為在相關運行級別下開啟：

4 結論

應用小波細節分析技術，提取了銅山口地區深部中酸性巖體產生的弱信號異常信息。采用“剝皮”方式，利用不同尺度空間的重力異常特征，刻畫了區內深部巖體和構造在不同尺度空間的形態，取得了豐富的深部找礦信息。目前部分成果已得到了實際檢驗。由此來看，小波變換作為重力異常分解的工具之一，特別是在提取深部目標地質體(或礦致異常)所產生的弱信號異常信息方面，確實是行之有效的。作為一種新而有效的位場分離途徑，小波多尺度分析方法為重力資料解釋和研究地殼結構提供了新的思路，在目前廣泛開展的第二空間找礦工作中必將發揮獨特的、重要的作用。

［1］曾華霖.重力場與重力勘探［M］.北京：地質出版社，2005.

［2］湖北地質礦產局.湖北省區域地質志［M］.北京：地質出版社，1990.

Application of 2D Wavelet Multi-scale Analysis on Depth Mineral Deposit Exploration at Tongshankou

LIANG Xuetang，ZHANG Mulin
(Hubei Institute of Geophysical Exploration Technology，Wuhan，Hubei430056)

Wavelet transformation is an important tool for decomposing the gravity anomaly.With it，gravity anomaly can be transformed into different scale anomalies.And the gravity anomalies with different scales correspond to the different ore bodies with different scales and buries.As a new and effective means for potential field decomposition，the wavelet multi-scale analysis method provides a new way for gravity data interpretation and the earth’s crust study，and plays important roles in the current deep exploration.

gravity exploration;Bouguer gravity anomaly;wavelet multi-scale analysis;logarithmic power spectrum;deep exploration

P631.1

B

1671－1211(2012)03－0263－05

2011－05－24;改回日期：2011－07－28

梁學堂 (1964－)，男，教授級高級工程師，應用地球物理專業，從事區域重力調查、地球物理及數據處理與研究工作。E－mail：lxt19641025@163.com

潘 瀟)