中國東部海域居里面深度特征

(1.山東科技大學 地球科學與工程學院，山東 青島 266590； 2.中國冶金地質總局 青島地質勘查院，山東 青島 266109；3.自然資源部第二海洋研究所 海底科學重點實驗室，浙江 杭州 310012)

作為大家熟知的地熱學界面，居里面能很好地反映出地殼熱流及磁性層的分布狀態，也在一定程度上反映出塊體的構造是否穩定，在礦產資源、地熱、大地構造及火山活動等研究中具有重要意義。已有眾多學者應用不同方法計算過我國東部及鄰近海域的居里面。早在1982年郝書儉等[1]就利用譜分析方法計算了渤海水域的居里面，認為該區域的居里面深度范圍為16～26 km。后來江為為等[2]應用中國東部磁異常數據反演計算了中國東部及其鄰域居里面深度的起伏特征，認為中國東海海域居里界面深度在10～17 km之間變化。而韓波等[3]基于功率譜法得到的東海區域居里面深度為15～29 km。此外，李春峰等[4]計算的東部海域的居里面深度為19.6～48.9 km，平均深度約31.7 km?？梢姴煌瑢W者的計算結果存在差異。2016年，熊盛青等[5]根據最新編制的1/100 萬航陸域磁異常數據，編制了中國陸域居里面深度圖，并分析了居里面與地溫梯度和大地熱流的關系。鑒于此，擬利用中國東部海域，包括渤海、黃海和東海的高精度船測磁異常數據，采用功率譜方法計算居里面，并繪制整個中國東部海域居里面深度圖，分析各海域居里面與主要斷裂、海底熱流之間的關系。

1 數據與方法

1.1 數據來源

海域磁異常數據為船載磁測資料，為了保證計算的完整性及其效果，對周邊陸地區域補充了第一版世界數字化磁異常圖(World Digital Magnetic Anomaly Map，WDMAM)[6]，數據高度陸區為大地水準面上方的5 km，分辨率為3弧分×3弧分的網格數據。采用迭代向下延拓的方法對數據進行融合，將陸域的磁力數據向下延拓到海平面，通過網格化使之與船載數據的網格間隔一致。研究范圍為114°E～129°E，22°N～41°N，采用墨卡托投影，WGS84橢球，中央經線117°E，標準緯線28°N。研究區的磁異常如圖1所示。

居里面作為一個重要的熱力學界面，能很好地指示上地幔和地殼的熱狀態。為了研究居里面和熱流之間的相關性，主要參考中國陸地區域大地熱流數據匯編的1～4版[7-10]、楊樹春等[11]南黃海的八個熱流數據以及全球熱流數據庫[12]中的部分熱流數據，收集了中國東部及鄰近海域900個熱流測量數據。

1.2 小波多尺度分解

功率譜法反演居里面深度需要計算磁性層頂面深度，而淺部干擾和數據中的噪聲對居里點深度的計算結果影響又較大，但可以采用圓滑濾波去掉磁異常高頻成分或采用向上延拓的方法來減小影響。眾所周知，小波多尺度分解能夠將重力、磁力數據精細地分解到不同尺度，用來反映不同尺度和深度的數據特點，故常被用于區域場的分解和分析。另外，由于研究區整體緯度跨度較大，磁傾角在29°～60°之間變化，采用固定磁傾角化極很難真實地反映磁異常變化，因此將研究區從北向南進行了分帶，并采用連續化極使化極窗口大小為200 km×200 km。綜上，對研究區化極磁異常做了6階的小波分解，得到各階的小波逼近和細節，通過頻譜分析，得出各自的似深度(圖2)。1階細節未出現低頻段的“直線”特征，說明1階細節主要反映出高頻隨機干擾；2階細節至6階細節的似深度分別為7.80、15.29、45.47、90.37、159.12 km；而1～6階磁異常逼近的似深度分別為12.34、14.78、44.47、102.66、140.03、155.57 km，可以看到3階細節和2階逼近似深度很接近，4階細節和3階逼近，5階細節和4階逼近以及6階細節和5、6階逼近的場源似深度接近，說明二階分解之后，小波逐漸將深部的信息反映出來。其中2階細節反映出淺部的不均勻磁性體，3階細節反映了磁性基底的變化。據此選取小波二階逼近結果作為反演居里面的資料。

圖1 中國東部近海磁異常圖

圖2 頻譜分析

1.3 居里面深度計算

采用針對組合磁異常的統計功率譜分析方法[13]，Zt為棱柱磁性體的頂層深度，Z0為中心深度，根據公式Zb=2Z0-Zt換算出磁性層底層深度值。居里面計算還需考慮滑動窗口，窗口太小不能獲得足夠的頻譜信息，窗口過大則降低了空間域的分辨率。Okubo等[14]認為窗口的尺度應該在磁源質心深度的12至13倍以上，才能獲得較為理想的居里面深度值。一般在計算時將窗口大小固定，以實現快速計算。本研究認為研究區磁性層平均中心深度大約為12 km，故采用子窗口大小為160 km×160 km。每個子窗口中，低頻段采用0.003～0.051頻段擬合Z0，高頻段則選擇0.057～0.076頻段計算Zt，計算結果見圖3。

圖3 中國東部及鄰近海域居里面深度及熱流點分布圖

2 結果分析與討論

2.1 居里面深度分析

從圖3可以看到，渤海區域對比華北陸區的居里面深度整體呈上隆狀態，淺于28 km。但在渤海盆地內居里面與基底起伏具有鏡像對稱關系，即坳陷區居里面較深。遼東灣坳陷為NE向的下坳區，深度為28 km上下；渤中坳陷居里面下凹不明顯。居里面沿F1斷裂呈NW向展布，深度為18～22 km。F2斷裂從濰坊進入萊州灣和渤海灣，控制渤海灣東部邊界，其在渤海區域是一條左行走滑斷裂，斷裂帶東側磁異常表現為寬緩負異常，西側整體表現為高磁異常區，幅值在100 nT左右，斷裂帶兩側為大面積的居里面隆起區。華北克拉通東部晚中生代以來發生克拉通破壞，巖石圈減薄已得到廣泛證實。Chen等[15]計算出的渤海海域、魯西隆起及F2斷裂帶附近的巖石圈厚度為60～80 km；Pn波速度偏低[16]顯示出巖石圈的拉伸減薄以及地幔物質的向上侵入。F2斷裂帶渤海段附近的地殼廣泛發育低速層[17]，從解析信號圖上看，斷裂呈串珠狀比較孤立的巖體分布特點[18]。作為F1斷裂與F2斷裂交匯的區域，渤中坳陷的重、磁異常均較高[19]。推測中新生代克拉通東部巖石圈減薄時期交匯處可能作為軟流圈熱物質的上涌通道，巖漿沿斷裂交匯處上涌改造基底，增強了基底的密度和磁化強度并抬高了該區域的居里面，因此渤中坳陷居里面下凹并不明顯。

北黃海居里面起伏狀態整體平緩，深度較淺，在17～25 km之間變化。北黃海海域的自由空間重力異常跟布格重力異常相似，表現為較強的正重力異常[20]，反映了地幔上隆的重力效應。磁異常卻為大面積的低緩負異常，在低值背景場上出現零碎的NW向正異常，可能是受到F5斷裂帶的影響。與渤海類似，屬于中國東部裂谷拉伸作用影響的巖石圈減薄區。楊艷秋等[21]通過分析重、磁資料認為北黃海盆地之下隱伏有埋深較淺的高密度、弱磁性結晶基底。這些反映了北黃海較弱的構造和巖漿活動。

南黃海北部盆地居里面起伏整體變化不大(圖3)，范圍為25～29 km。蘇北及南黃海南部盆地居里面深度較大，范圍為25～32 km。介于黃海北部盆地和南部盆地之間的中部隆起帶表現為大塊狀正異常，在中間和東端被負異常切割而不完整，磁異?！胺侄尾贿B片”[22]，居里面在西部為下凹區，最深達33 km，走向EW向，東部相對上隆，深度為28 km，走向不明顯。關于南黃海中部隆起的大面積團塊狀正磁異常，有學者認為是古陸塊的體現，有學者則認為該區存在上地幔隆起，引起深部熱活動加劇，地幔物質上涌侵入下地殼，增加了下地殼基性物質成分[23]。吳其反等[24]通過三維反演發現20～25 km深度存在一規模巨大的NW向展布的強磁性體，推測為華北變質巖系的殘留體。吳招才等[25]計算了該區的磁源重力異常，認為該區可與華北地塊古陸塊的磁源重力異常規模相比。該區莫霍面整體抬升區，深度為24～26 km；居里面則是西部較深而東部較淺，東部的地震剖面[18,22]也觀測到深部巖體侵入。因此推測南黃海中部隆起團塊狀正磁異常是古陸塊和巖漿侵入的綜合反映。

夾在F3斷裂與F4斷裂之間的千里巖隆起帶居里面較淺(圖3)，深度為19～24 km。千里巖隆起存在高速推覆體，重力異常值大于20 mGal，屬于高密度變質巖特征。千里巖隆起帶居里面淺和密度高的特征可能與侏羅紀以來太平洋板塊俯沖形成晚期較強的張裂或深部巖漿與熱流活動有關[4]。邱寧[26]則認為在揚子板塊與華北板塊發生碰撞時，揚子板塊深俯沖過程中拉動相鄰的華北板塊的一部分物質向下俯沖，引起上覆塊體物質遭到了俯沖剝蝕，形成了不對稱的雙向俯沖，板塊俯沖受阻后發生滯留、下沉引起地幔擾動，而在老的構造活動造成的地殼減薄帶地幔流上涌，使得該區的居里面上隆。

F6斷裂帶磁異常整體呈NE走向的正值分布(圖3)，內部正、負變化劇烈，正負異常圈閉零散、細小，方向無規律。居里面沿F6斷裂帶上隆，深度為16～22 km，代表一個重要的塊體邊界，但居里面的上隆趨勢在濟州島北側發生改變，變為下凹。

浙閩隆起帶的居里面呈NNE向上隆(圖3)，深度為16～24 km，與中國東南沿海晚中生代以及新生代近NNE向火山巖帶大體對應；隆起帶的重、磁場也表現為NNE向的細節形態[19]，顯示出其受到了燕山晚期和喜山早期的NNE向的斷裂切割和巖漿侵入作用。東海區域磁異常大致在北緯29°的南、北兩側分布不同，南側東西分帶性明顯，北側的低磁異常分帶性則不明顯。其居里面呈較為明顯的“東西分帶”特征，基本與東海的構造劃分相一致。相對于浙閩隆起帶，東海陸架盆地的甌江凹陷、武夷低凸起西南部、西湖凹陷東南部、基隆凹陷及福江凹陷南部的居里面深度均呈下凹狀態，深度為26～29 km。陸架盆地由西向東形成的年代逐漸變新, 西部坳陷主體形成于白堊紀，東部坳陷帶主體則形成于始新世,漸新世以來成為陸架盆地的沉積中心，具有較深的巖漿巖分布。釣魚島隆起帶整體表現為高磁異常，其屬于中新世古東海陸架邊緣地質體分離出的一部分[27]。隆起帶的中部和南部居里面表現為28 km左右的下凹帶，體現了燕山期和喜山期強烈的巖漿作用改造使得該處地殼增厚和基底隆起。

因此，研究區居里面深度變化特征較好地反映了該區域的深部構造格局，在構造環境穩定、地殼厚度正常的區域，居里面主要顯示為塊狀坳陷且起伏變化平緩，這可能與相對較弱的中、新生代構造活動有關；而在構造活動區邊緣或者中、新生代造山系，居里面主要表現為淺的區域性隆起。

2.2 居里面與熱流的關系

一般來說，高熱流區域居里面深度小，低熱流區域居里面深度大。然而在構造復雜的區域以及古老的地區中，可能體現不出這種相關性，這與許多因素相關，比如當地的地溫梯度、巖石生熱率、地幔熱流、斷層活動和巖性特征等均會擾亂傳導關系。中國東部海域(不包括沖繩海槽及以東區域)構造比較穩定，故這種地表熱流與深層地熱之間的簡單傳導關系可能不明顯。

京津唐地區的熱流值變化較大，為0～120 mW/m2，居里面深度與熱流的分布基本一致，可能是由于該地區西部和北部為大面積的山區，居里面較深，而其他地區屬于平原，居里面深度下凹和上隆呈塊狀分布，起伏變化較大。陸區合肥盆地居里面較深，大約30 km，熱流值大部分在70 mW/m2以內；而其東部的蘇北盆地的熱流值范圍在54～83 mW/m2之間，平均68 mW/m2，和南黃海南部盆地的平均熱流值相同。這是因為南黃海南部盆地與蘇北盆地均位于揚子準地臺基底之上,無論是形成時代和還是演化過程均具有相似性，現在正在經歷熱冷卻過程[11]。隨著熱流值的增高，可明顯看出熱流分布區域逐步向東部海域轉移，居里面深度也逐步變淺。

渤海灣的熱流值與居里面深度有很好的對應關系，其大地熱流最高85 mW/m2，最低50.3 mW/m2；渤中凹陷熱流值略微高于60 mW/m2，與大洋平均熱流值相當，渤海地區熱流值大于70 mW/m2大體沿NE向展布，從遼東灣到渤南凸起以及墾東凸起高熱流分布呈現出與盆地居里面起伏形態大體對應的特征，反映該盆地中新生代曾經歷的裂谷作用及巖漿活動。

東海陸架盆地熱流值55～88 mW/m2，平均熱流值70.5 mW/m2，屬于正常熱流值區域，熱流點分布在甌江凹陷、漁礁凸起東南部以及西湖凹陷，位于居里面深度下凹區域。120 mW/m2以上的高熱值主要分布在沖繩海槽區域，而中、低熱流值在研究區內均有分布。沖繩海槽是高熱流值集中地區,也是低熱流值廣泛分布的地區，表現為低熱流值和高熱流值在同一區域出現,可能是在垂直于海溝軸的方向上熱流總和趨于平衡[28]，即在熱液噴出口周圍的小范圍的斷層作用使冷水下滲,形成低熱流循環，該區域居里面表現為隆起。

圖4為研究區的居里點埋深與熱流值分布散點圖，為了繪圖方便，將熱流值大于350 mW/m2的數據取350 mW/m2，大于300 K·km-1的數據取300 K·km-1。通過計算研究區居里面深度與熱流、地溫梯度的自相關系數(Pearson相關系數)后發現，結果分別為-0.389 6與-0.516 0，說明該區域居里面深度與熱流值之間為弱負相關，與地溫梯度為中等程度負相關，但相互間沒有很好的線性關系。隨著居里面深度逐漸增加，地溫梯度和熱流值均降低；居里面深度大于29 km時，東部海域的熱流值明顯降低，小于100 mW/m2；而在居里面隆起的區域(<28 km)，東部海域熱流值變化范圍大。地溫梯度也存在類似特征，居里面深度大于29 km時，地溫梯度大部分小于50 K·km-1。

圖4 居里深度與熱流值及地溫梯度分布散點圖

3 結論

選取中國東部及鄰近海域磁力數據的二階小波逼近結果反演居里面，對居里面及熱流特征進行了分析，得到結論如下：

1) 中國東部海域居里面深度為16.8～36.6 km。在構造環境穩定、地殼厚度正常的區域，如蘇北—南黃海盆地、東海陸架盆地，居里面深度在28 km以上，起伏變化平緩，可能與相對較弱的中、新生代構造活動有關；在構造活動區邊緣或者中、新生代造山系，如渤海盆地、千里巖隆起區及沖繩海槽，居里面深度則以上隆為特征。

2) 居里面深度與熱流、地溫梯度的自相關系數的計算結果分別為-0.389 6與-0.516 0，顯示出居里點深度與熱流值為弱的負相關性，與地溫梯度為中等程度負相關性，均未呈現很好的線性關系。

3) 南黃海中部隆起的居里面深度西部下凹而東部相對上隆，綜合前人工作及重磁震資料，認為可能是古陸核與巖漿侵入的綜合反映。