長白山地區航磁異常居里面計算及其地質意義

蘇曉軼，李福文，屈旭鈞，張天驥

1．長白山天池火山監測站，吉林 安圖 133613；2．吉林省勘查地球物理研究院，長春 130012；3．吉林省地質勘查基金管理中心，長春 130061

0 引言

居里面是指磁性礦物在一定溫度下從有磁性轉變為無磁性的居里(點)等溫度面, 簡稱居里面,是磁性層的下界面,反映了巖石圈的熱狀態而不是確定的地質構造面。申寧華等[1]曾指出居里面埋深對油氣生成、保存有影響,反演居里面深度對大地構造研究和油氣資源評價與勘探具有重要意義。此外,反演居里面對地熱能的開發利用以及火山噴發的監測與預報等也具有重要意義。

前人曾對居里面的反演做過很多研究,劉天佑[2]曾計算松遼盆地的居里面深度，郝天珧等[3]和江為為等[4]曾對賀蘭山—龍門山以東的中國陸海及鄰區進行居里面反演,并利用反演結果劃分洋、陸構造單元。胡旭芝等[5]、吳招才等[6]、高德章[7]、南方舟等[8]用Spector-Grant的功率譜法反演東北、南海、東海及華北地區居里面,并研究了這些地區的地殼熱狀態。申寧華[9]推導了快速計算界面位場異常的正反演公式，其原理及性質與Parker[10]、Oldenburg[11]算法基本一致，該算法具有速度快、功能強等特點。

由于功率譜法原理假設條件與客觀差距較大，實際應用效果不如Parker方法。本文采用改進的Parker-Oldenburg反演迭代算法，計算了長白山地區居里面深度，實際應用效果較好。利用居里面分析該區的熱狀態,并準確預測了研究區內多個地熱有利區。

1 居里面反演的Parker-Oldenburg理論

設地下有物質界面A，形狀如圖1所示，假定其上部物質的磁化強度M為零，下部的磁化強度不為零，還假定M的方向垂直向下。A面下某點Q(ξ,η,ζ)處體積元dV對地表測點P(x,y,0)在垂直磁化情況下所產生的垂直磁異常為ΔZ，設直角坐標系Z軸垂直向下，則推導過程為：

(1)

式(1)中，μ0為真空磁導率，T由式(2)表示：

(2)

假設A面的平均深度為h，在ζ=h處將K用泰勒級數展開得到：

(3)

式(3)中，K(n)(x-ξ,y-η,h)是在ζ=h處K對ζ的n階導數，利用上式將K對ζ從(h+Δh)到h積分，其中Δh是A面相對于h的距離，在ζ=h以上Δh是負值，在其下面Δh為正值，且Δh是(ξ，η)的函數由式(2)、(3)可以得到：

(4)

設M在A以下只是(ξ,η)函數，則A以下全部物質對P點產生的磁異常為：

(5)

(6)

(7)

(8)

(9)

迭代反演的基本原理:根據場值計算界面起伏的初始值，根據該初始界面模型進行正演計算，利用正演計算值與實測值之差修正界面起伏。利用修正后的界面進行正演計算，將正演值與實測值進行比較，找出參差，并對交錯進行修正和逐次迭代，直至界面模型和實測場擬合的均方差小于預設要求，得到最終界面分布。因此，界面模型的選擇和相應正演計算方法的精度直接影響計算結果的準確性。合理地給出界面波動的初值和修正值將直接影響迭代的收斂速度。

本文涉及的單層界面反演利用界面反演中比較常用的Parker-Oldenburg方法，該方法采用頻率域計算公式，速度快、精度高，且使用方便。但反演公式包含由上而下延拓因子，這些因子通過放大高頻干擾成分使反演過程變得不穩定。利用式(8)反演單界面A的起伏Δh(ξ,η)，可以將公式進行移項，寫成：

(10)

若令等式右端的Δh值為第j次的近似值，而等式左端的Δh為第(j+1)次的近似值，那么這就是Parker-Oldenburg提出的迭代反演公式(左端的△h經傅里葉變換后，得到空間域的Δh(ξ,η)值)。但是這樣的迭代，因每次求得是△h值，而非修正量，不易保證逐次逼近收斂。△h<

(11)

(12)

(13)

(14)

(15)

這種迭代所計算的值與實際值相差越來越少，因為每次求Δh(ξ,η)增量的線性逼近，使迭代達到收斂。只需加一適當的正則化因子防止高頻因子的發散作用。實際證明這種迭代反演的方法結果較好。

以上在求Δh(ξ,η)時，須事先已知A界面的大體平均深度h及界面以下的磁化強度M。如果存在雙界面，則還需知道H和ΔH(ξ,η)才能計算出界面異常。

圖1 界面A示意圖Fig.1 Diagram of interface A

2 長白山地區居里面反演結果與分析

2.1 航磁異常特征

長白山地區的航磁異常就異常強度而言，研究區可分為南、北兩個異常區，南部為高磁異常區(+200～+550 nT)，集中出現在南端的漫江—長白縣一帶；北部廣大地區為中-低磁異常區(-200～+200 nT)。根據該區的磁異常形態與構造單元的對應關系，研究區可分為六個磁異常區(圖2)：蒲柴河—露水河鎮以東，松江鎮以西地區，以夾皮溝地塊為主體，低緩的低磁異常，呈北西向展布(Ⅰ區)；撫松—泉陽鎮以北，蒲柴河—露水河鎮以西地區，龍崗地塊為主體，磁異常正負交替，呈北東向展布(Ⅱ區)；撫松—泉陽鎮以南，花山鎮—松江河以北，板石溝地塊為主體，環狀高磁異常為主(Ⅲ區)；花山鎮—松江河以南，松江鎮以西，老嶺地塊為主體，大面積低磁異常為主(IV區)；望天鵝以南，長山縣以北，長白地塊為主體，強磁異常，呈近東西向展布(V區)；松江鎮以東，和龍地塊為主體，磁異常正負交替，呈北西向展布(VI區)。

Ⅰ.夾皮溝地塊；Ⅱ.龍崗地塊；Ⅲ.板石溝地塊；Ⅳ.老嶺地塊；V.長白地塊；VI.和龍地塊.圖2 長白山地區航磁異常圖Fig.2 Aeromagnetic anomaly map of Changbai Mountain area

航磁異常形態特點不但反映了龍崗地塊與和龍等地塊及古生代褶皺斷裂帶之間的構造演化關系，而且反映了該區淺表層的地質巖礦石組成特點。

(1)北部低緩異常區(圖2 Ⅰ區)

該區為低緩異常區，低緩正值場，在大荒溝—老嶺一帶正磁場為太古代鉀長花崗片麻巖引起，在興旺屯—老森嶺一帶正磁場為中元古代花崗閃長巖引起。金山南負場區為侏羅世鉀長花崗巖引起，四方頂子現負場區推測為深部存在酸性巖體。該區兩條深大斷裂，磁異常明顯。

(2)西北部正負交替磁異常區(圖2 Ⅱ區)

該區磁場為正負交替場，區內有多個時代地質出露，場態復雜，異常強度變化較大，強度在-300～+300 nT之間。靖宇附近為正磁異常，范圍大、形狀不規則，強度為+50～+350 nT之間，對應的主要為侏羅系火山巖和新生代玄武巖，其中高值異常為燕山期花崗巖和閃長巖，并有一處大營子小型鉛鋅礦床。萬良附近正磁異常表現為南北兩個低緩異常，區內南部為侏羅系安山巖，北部為中太古代鉀長花崗質片麻巖，異常北側有三處小型鐵礦床。抽水村北為正磁異常，異常近南北方向，南部異常較大，北部異常較小，呈亞葫蘆狀，為低緩異常，強度+350 nT，地表大部出露太古代表殼巖，對尋找“鞍山式”沉積變質鐵礦有利。

(3)西部高磁異常區(圖2 Ⅲ區)

該區磁場以環形正磁場為主，在松江河附近為正磁異常，范圍大，形態較規律，對應的主要是侏羅系火山巖和新生代玄武巖，其中高值異常為燕山期花崗巖和閃長巖。

(4)中部低磁異常區(圖2 Ⅳ區)

該區磁場以低磁為主，但異常強度變化較大，強度在-300～+200 nT之間。該區幾乎被火山巖層覆蓋，但由于大量酸性巖體的侵入，導致大面積低磁異常。低磁異常區為酸性侵入巖區，高磁異常為第三第四紀基性火山巖覆蓋較厚區。在正磁場中從望天鵝—天池一帶出現多處環形異常，區內主要出露為玄武巖，磁場體現了玄武巖特征，出現的環形異常推測為火山構造引起。

(5)南部強磁區(圖2 V區)

該區磁場以正磁場為主，六道溝北局部正磁異常，呈橢圓形、走向北東，強度+400 nT，異常為侏羅系果松組安山巖，晚侏羅世閃長巖引起，六道溝小型矽卡巖型銅鉬礦床位于異常中。長白縣八道溝—十四道溝之間，高磁異常呈東西向展布，主要為裸露的龍崗群古老基地巖系。

(6)東北部負交替磁異常區(圖2 VI區)

該區磁場為正負交替場，但異常強度變化不大，強度在-100～+250 nT之間。該區北部有一條深大斷裂，異常明顯。高磁異常主要為一套閃長質、花崗質深成侵入體。

2.2 長白山地區居里面及其分析

利用本文提及的磁性界面求解方法，計算了長白山地區的居里面深度(圖3)。長白山地區磁性層底面即居里面深度范圍在12.8～15.8 km，靖宇—仁義砬子地區局里面深度較淺，在12.8 km；松江河以及長白山天池西南居里面深度較大，在14～15 km；長白山天池地區居里面較淺，深度在12.8 km±。

由計算結果分析：長白山天池作為一個活動的火山，居里面較淺是合理的。大地電磁測深的結果揭示[12]，該處深為12 km的地方，有低電阻的特征，推測是未固結的巖漿囊。靖宇—花山鎮一帶居里深度也較淺，花山鎮附近有老三隊溫泉，表明深部地殼處于熱狀態，而靖宇地區作為臺區，表層出露的巖石是太古代地層，但其深部可能存在深部巖漿的底侵，造成該區居里面較淺。長白山的東北部和西南部，居里深度較深，該區分布著較厚的元古界和古生界，居里深度較深是合理的。

圖3 長白山地區居里面深度圖Fig.3 Curie depth map of Changbai Mountain area

3 結論

(1)長白山地區由于地溫梯度較高，居里深度總體較淺，居里深度為12.8～15.8 km，較全球的居里深度20 km為淺。

(2)長白山地區的居里深度平面上有變化，最深處在松江鎮右側，最淺處在長白山與靖宇—臨

江—撫松一帶，與當前地熱出露狀態符合。

(3)長白山地區的地熱主要有利地區除在天池地區以外，松江河以西的撫松、仙人橋、老三隊地區都屬居里深度較淺地區，為較好的地熱遠景區。