歐亞大陸東緣存在一個巨型深斷裂系統——基于衛星重力的新發現

黃宗理, 王　典, 嚴加永, 張　懷

1)中國地質科學院礦產資源研究所, 國土資源部成礦作用與資源評價重點實驗室, 北京 100037; 2)國土資源部咨詢研究中心, 北京 100035; 3)中國地質大學(北京), 地下信息探測技術與儀器教育部重點實驗室, 北京 100083; 4)中國科學院計算地球動力學重點實驗室, 北京 100049

?

歐亞大陸東緣存在一個巨型深斷裂系統——基于衛星重力的新發現

黃宗理1, 2), 王典3), 嚴加永1)*, 張懷4)

1)中國地質科學院礦產資源研究所, 國土資源部成礦作用與資源評價重點實驗室, 北京 100037; 2)國土資源部咨詢研究中心, 北京 100035; 3)中國地質大學(北京), 地下信息探測技術與儀器教育部重點實驗室, 北京 100083; 4)中國科學院計算地球動力學重點實驗室, 北京 100049

摘要:通過對衛星重力數據的精細處理, 發現歐亞大陸東緣存在一個規則分布的重力異常組合, 這個異常組合與中國大陸上的已知斷裂嚴格對應。SinoProbe探測計劃實施的反射地震探測出這些斷裂下方存在Moho破裂, 幔源物質上涌, 剩余重力異常, 莫霍面破裂, 高密度的幔源物質上涌和深斷裂之間存在密切聯系。依照這個聯系規律, 對已有斷裂的延伸補充, 顯示出歐亞大陸東緣存在一個巨型深斷裂系統, 主斷裂帶南起中國廣州向北延伸直至鄂霍次克海, 綿延逾3000 km。沿主斷裂帶東側近似等間距的平行分布著9條北東向斷裂, 北東向斷裂向東入海至大陸邊緣。采用數值方法模擬了斷裂系統的形成過程, 在菲律賓板塊和太平洋板塊近NW方向構造力的作用下, 歐亞大陸東緣產生NNE向走滑斷層, 其東側生成彼此平行, 間距大體相等的派生斷層。斷層形成過程中, Moho破裂, 大量地幔物質上涌, 形成中國東部中生代以來的巨量火山巖漿活動, 控制了內生金屬礦床的分布。這個斷裂系統近代還在活動, 導致了郯城8.5級大地震和長白山的火山噴發。

關鍵詞:衛星重力異常; 斷裂系統; 歐亞大陸

本文由國家自然科學基金項目(編號: 41574133; 41104061)和國家深部探測技術與實驗研究專項(編號: SinProbe-08)聯合資助。

由于重力測量可以提供大面積的宏觀圖像和深部地質體的信息, 利用重力場研究區域地質構造的成果越來越多, 如劉光鼎等(1997)根據重力場特征提出中國大陸具有“三橫、兩豎、兩個三角”的宏觀格架; 滕吉文等(1982)給出了“中國地殼構造的基本輪廓”; 馬宗晉等(2006)使用重力水平梯度分析了中國大陸的地質構造; 侯遵澤和楊文采(1997)使用重力場的小波變換研究中國的地殼構造; 涂廣紅等(2006)使用剩余重磁異常研究東北地區的地質構造; 路曉翠等(2012)用小波變換重力資料研究郯廬斷裂等。在對重力場的研究中, 一般認為布格重力異常是由Moho面起伏和地殼內部局部密度不均勻兩種因素引起的。不同尺度的重力異常對應不同規模的地質體, 區域重力異常主要反映了大尺度Moho面的起伏變化, 而地殼內局部密度不均勻體則形成局部重力異常。

我們在研究歐亞大陸東部地區的衛星重力場時,對布格重力異常進行了區域場和局部重力異常的分離。發現在一定尺度的空間域濾波或一定頻帶寬度的頻率域濾波后, 得出的剩余重力異常在歐亞大陸東部邊緣地區出現了一個彼此平行、間距相等的規則組合。進一步的研究發現, 這些剩余重力異常在中國大陸上的部分與已知的地質斷裂完全重合。

1　數據處理

研究區地理位置見圖1a, 經緯度范圍為105°—145°E, 20°—55°N。使用的衛星重力數據和地形數據源自美國全球衛星重力數據庫, 數據編制日期是2011年6月22日, 數據網度1′×1′。經過對衛星數據的去噪濾波, 中間層改正和地形改正后, 得到研究區布格重力異常如圖1b。

本文研究布格重力異常中的局部異常, 重點是Moho附近小尺度變化引起的局部重力異常。為消除區域場的影響, 分別采用空間域的滑動窗口濾波和頻率域的高通濾波兩種算法, 空間域中逐次選擇濾波窗口邊長為14′、16′、18′、20′、22′、24′、26′,頻率域中逐次選擇濾波波長分別是100 km、200 km、300 km、400 km進行試驗。最后確定采用空間域算法, 濾波窗口大小確定為20′×20′(約相當于30 km×30 km)。計算得到的區域重力異常如圖1c和剩余重力異常如圖1d。比較圖1b、c、d三張重力圖可見, 剩余重力異常圖1d已經有效地消除了區域場的影響。

為檢驗衛星重力異常的可靠性, 任選四條剖面將從衛星數據計算出的剩余重力異常與使用地面測量重力數據計算出的剩余異常進行對比。從L1線、L3線、L4線和L5線衛星重力求取的剩余布格重力異常和地面重力求取的剩余異常對比結果(圖2),可以看出這兩種方法得到的剩余重力異常吻合較好,說明使用衛星重力數據計算剩余重力異常是可行的,計算出的結果具有較高的可靠性。

2　剩余重力異常分布特征及其與地質斷裂的對應關系

剩余重力異常在中國東部出現兩條明顯的北北東向的條帶狀異常(圖3), 一條沿大興安嶺、太行山、武陵山一線, 異常強度高規模大, 近南北向的條帶異常在華北地塊北緣, 以及揚子與華北地塊交界處被東西向異常所截斷, 這一個剩余異常帶對應的是中國東部著名的重力梯度帶, 馬宗晉等(2006)認為這個重力梯度帶反映了地殼的陡變, 梯度帶兩側地殼厚度變化達10 km。

第二個條帶狀異常南起廣州, 經南昌, 南京,穿過長江, 經山東半島, 過渤海灣, 經沈陽, 長春,哈爾濱, 從佳木斯以北出境, 沿黑龍江下游進入鄂霍次克海, 綿延逾3000 km。這個近南北向異常在中國東北, 走向轉為北東。異常東側(含異常北段)共有9條走向北東的正異常, 這九條異常彼此平行,兩條異常之間的距離大體相等, 異常向東伸向大陸邊緣, 9條異常構成一個奇妙的規則組合。平行①號近南北向主體異常做一條縱剖面AB, 在剖面上剩余重力異常的規則分布規律表現得更加清楚(圖4)。

將剩余重力異常圖與中國主要斷裂圖(程裕琪, 1994)疊合在一起, 由圖3可以看出:

①號重力異常南段與贛江斷裂, 中段與郯城廬江斷裂, 北段與依蘭—舒蘭斷裂相對應;

圖1　研究區地理位置圖(a)、研究區布格重力異常圖(b)、研究區區域重力異常圖(c)和研究區剩余重力異常圖(d)Fig. 1　Geographic location of the study area (a), Bouguer gravity anomaly (b), regional gravity anomaly (c) and residual gravity anomaly (d)

圖2　衛星重力數據與地面重力數據對比圖Fig. 2　Comparison between satellite gravity and ground gravity

②號重力異常與廣州北東向斷裂對應;

③號重力異常與邵武—河源斷裂相對應;

④號重力異常與紹興—萍鄉, 景德鎮—宜豐斷裂帶相對應;

⑥號重力異常與揚子陸塊北緣斷裂帶相對應;

⑦號重力異常與商南—榮成板塊縫合帶相對應;

⑧號重力異常與鴨綠江斷裂相對應;

⑨號重力異常對應敦化—密山斷裂。

圖3　剩余重力異常與斷裂疊合圖Fig. 3　Overlap map of residual gravity anomaly and faults

圖4　AB剖面的剩余重力異常(紅圈表示剩余重力異常編號)Fig. 4　Residual gravity anomaly along AB profile (red circle indicates serial number of residual gravity)

9條重力異常中, 有8條在中國大陸境內都和已知斷裂吻合得較好, 只有⑤號異常目前沒有已知斷裂與之對應。⑤號異常從南京沿長江入海。已有人根據該區各種地球物理信息, 指出長江下游至泰州為一隱伏深斷裂帶(陳毓川, 2007)(暫稱為長江下游隱伏斷裂)。石油地質工作也已在⑤號異常海域部分發現了一組斷裂(馬文璞, 1992)。

3　剩余重力異常的地球物理解釋

圖5　穿過依蘭—舒蘭斷裂的剩余重力剖面(a)與反射地震圖像(b)Fig. 5　Residual gravity profile (a) and seismic reflection image (b) across the Yilan–Shulan fault

圖6　穿過敦化—密山斷裂的剩余重力剖面(a)與反射地震圖像(b)Fig. 6　Residual gravity profile (a) and seismic reflection image (b) across the Dunhua–Mishan fault

引起重力異常的根本原因是地下存在著高密度物體, 9條重力異常與已知斷裂的嚴格對應說明在這些斷裂下方有高密度物體沿著斷裂分布。問題在于這種綿延數千千米的條帶狀高密度物體是什么？SinoProbe-02專項最近完成的東北和華南兩條深地震反射剖面測量, 獲得了包含Moho面的地殼的精細圖像, 給出了對重力異常的解釋結果。

圖5是切過依蘭—舒蘭斷裂的反射地震圖像與剩余異常剖面L1, 在雙程走時12秒附近連續的Moho反射波界面十分清晰, 但在依蘭—舒蘭斷裂下方, 原來連續的Moho面出現破裂, Moho在破裂處隱隱形成一個直徑約3秒(9 km)的橢圓, 向上隆起約4 km, 其上方出現一系列復雜波形一直貫穿到地表, 構成一個垂直的花狀構造, 從Moho直到地面依蘭—舒蘭斷裂(熊小松等, 2011)。地震圖像清晰地描繪出依蘭—舒蘭斷裂切穿了地殼, 斷裂下方的Moho破裂上隆, 密度較大的幔源物質由破裂處上涌進入地殼, 在地殼內形成幔源和殼?；旌系膸r漿巖巖體, 隆起的Moho和高密度的基性巖漿巖體產生了剩余重力異常。理論計算表明, Moho隆起1 km可產生約3 mGal重力異常, 隆起的Moho和高密度巖漿巖體沿斷裂走向分布, 形成與已知斷裂嚴格對應的條帶狀剩余重力異常。地震圖像中的花冠狀斷裂則是走滑斷層的典型表現。

圖7　穿過贛江斷裂的剩余重力剖面(a)與反射地震圖像(b)Fig. 7　Residual gravity profile (a) and seismic reflection image (b) across the Ganjiang fault

圖6是切過敦化—密山斷裂的地震圖像與剩余異常剖面L2, 圖7是切過贛江斷裂的地震圖像與剩余異常剖面L3。地震圖像的特點和依蘭—舒蘭斷裂完全相同, 說明敦化—密山斷裂和贛江斷裂也是切穿地殼的深斷裂, 深斷裂下的Moho破裂, 幔源物質上涌形成了剩余重力異常, 典型的花狀構造表示這兩條斷層也是走滑斷層(熊小松等, 2009)。

①號異常中段目前還沒有深達Moho面的反射地震資料, 史大年等(2012)對切過⑤號異常做的人工地震寬角反射/折射測量, 橫穿過郯廬斷裂和長江下游隱伏斷裂, 地震測量結果在斷裂下方Moho有約4 km的隆起。朱日祥等(2010)在渤海灣附近做過多條人工地震寬角反射/折射測量, 有兩條剖面分別切過依蘭—舒蘭斷裂南端和郯廬斷裂山東段,在這些斷裂下方, 地震也測出Moho有明顯的隆起和幔源物質上涌。雖然人工地震寬角反射/折射不能提供如同反射地震那樣清晰的地殼結構圖像, 但斷裂下方存在Moho破裂、局部隆起以及幔源物質上涌是清楚的。

地震測量揭示出剩余重力異常和斷裂之間的內在聯系: 斷裂切穿Moho, Moho破裂隆起, 幔源物質上涌, 形成剩余重力異常, 剩余重力異常和深斷裂是密切相關的統一體, 故可以根據剩余重力異常確定深斷裂。目前已知的斷裂都是根據地面地質工作確定的, 只能在陸地上劃定, 也很難了解斷裂的深度。而衛星可以在陸域和海域以相同網度和精度進行重力測量, 由此計算出的剩余重力異常能夠追蹤斷裂的延伸, 判斷斷裂的深度, 給出陸域和海域里斷裂的全貌, 彌補地面地質工作的不足。

4　歐亞大陸東緣的深斷裂體系

依照剩余重力異常與斷裂的對應關系, 可將在中國大陸上地質已經劃定的斷裂沿剩余重力異常向中國境外和海域自然延伸, 把已知的斷裂補充完整。

①號異常北段對應的依蘭—舒蘭斷裂帶, 從剩余重力異?？磻虮毖由斐鼍? 沿黑龍江下游進入鄂霍次克海。①號異常從南端的廣州到北端進入鄂霍次克海, 整體是連續的, 說明①號異常所對應的贛江斷裂帶, 郯廬斷裂帶, 依蘭—舒蘭斷裂帶, 從南到北整體上是連續的。這個連續的斷裂帶在衛星地形圖上也可清晰識別。

②、③、④號三條重力異常對應的廣州北東向斷裂、邵武—河源斷裂、紹興—萍鄉和景德鎮—宜豐斷裂均位于中國華夏地塊之內, 對這一地區的構造存在較多的不同認識。從重力異?？? 以①號異常南段的贛江斷裂為界, ②、③、④三條異常與①號異常南段的贛江斷裂一起構成了一個與西側不同的系統。前人也已發現贛江斷裂東西兩側晚中生代火山巖的分布是不同的(毛建仁, 2013)。該區還分布著一系列平行于海岸線的斷裂, 從重力異?？此鼈兒廷凇ⅱ?、④號異常對應的深斷裂也不是一個系統,而且沒有切穿地殼。

④號異常對應的紹興—萍鄉, 景德鎮—宜春斷裂帶, ⑤號異常對應的長江下游隱伏斷裂帶, ⑥號異常對應的揚子陸塊北緣斷裂帶, ⑦號異常對應的商南—榮成板塊縫合帶應當向東延伸進入大海, 直至大陸邊緣。

⑧號異常對應的鴨綠江斷裂的主體應在沿遼東半島海岸線的海里, 斷裂的東段沿中朝邊界并穿過長白山火山地區。

⑨號異常與敦化—密山斷裂完全對應。

綜合以上地質已經劃定的斷裂和根據剩余重力異常對已知斷裂的追蹤補充, 可以看出歐亞大陸東部邊緣有一條南起廣州北至鄂霍次克海的深斷裂帶, 斷裂帶的東側, 沿北東方向自南向北依次排列著廣州北東向斷裂、邵武—河源斷裂、紹興—萍鄉和景德鎮—宜豐斷裂、長江下游隱伏斷裂、揚子陸塊北緣斷裂, 商南—榮成板塊縫合帶、鴨綠江斷裂、敦化—密山斷裂、依蘭—舒蘭斷裂。這些斷裂彼此平行, 間距大體相等(300 km), 向東延伸至歐亞大陸邊緣。北東向斷裂與郯廬斷裂和贛江斷裂一起構成一個巨型深斷裂系統。

圖8　數值模擬模型參數圖Fig. 8　Parameters of numerical simulation model

圖9　斷層形成過程的數值模擬示意圖Fig. 9　Numerical simulation of fault system formation process

圖10　G點位置圖Fig. 10　10 G point location map

圖11　剩余重力異常與鉛同位素圖Fig. 11　Residual gravity anomaly and lead isotope

5　巨型深斷裂帶系統的數值模擬

將歐亞大陸作為一個整體, 西部穩定不變, 東部邊緣考慮為自由邊界, 根據Critical Coulomb Wedge theory, 采用經典的剖面上的棱柱體模型。模型參數以及示意圖如圖8。底部為滑脫層, 推力施加在左側的塊體邊緣, 不考慮塊體以及滑脫面上的內聚力, 取地殼的平均密度為ρ=2900 kg/m3。根據實際地質體, 模型的外部邊界的控制參數如下:

圖12　剩余重力異常與成礦區帶圖Fig. 12　Residual gravity anomaly and metallogenic belt

β=0.2°,Li=400km,αi=10°,L=3700km,α=0.01°, He=111km,Le=3000km,αe=0.01°

取φB=30°,φD=3°,φRf=15°, 模擬結果是: 在南北向力的作用下, 先生成了廣州北東向的斷裂,隨著繼續向北推擠, 相繼形成了邵武—河源斷裂、紹興—萍鄉和景德鎮—宜豐斷裂、推測的長江下游隱伏斷裂、揚子陸塊北緣斷裂, 商南—榮成縫合帶的再次激活, 鴨綠江斷裂, 隨后斷層在后盤的近端再次激活, 然后在前緣再次生成敦化—密山斷裂以及最后一個依蘭—舒蘭斷裂的北段, 斷層形成過程見圖9。

數值模擬還給出了斷裂的間距和每條斷裂平移的距離。圖10中的G點表示南北向斷裂與北東向斷裂交匯點的位置, 由此可見這9條近北東方向斷層以等間距排列, 間距約為300 km, 各條斷層自南向北逐漸位移, 華南移動的較少, 到商南—榮成板塊縫合帶時移動約550 km, 到依蘭—舒蘭斷裂北端時移動約900 km。

6　巨型深斷裂系統上的地球化學現象

張理剛等(1993)根據巖石中鉛的同位素數據的差異, 將中國東部大陸劃分為不同的地塊, 他稱其為不同的地球化學省(圖11), 對比地球化學圖和剩余重力異常圖可見, 根據地球化學數據劃出的地塊分界和根據地球物理資料劃分的斷裂位置十分一致,以①號南北向重力異常為界, 東西兩側的地球化學分界在走向, 形態上完全不同, 說明東西兩側的地球化學差異很大。在東側, ②、④、⑥、⑦、⑧號剩余重力異常的位置和地球化學省的邊界幾近完全重合。地球化學和地球物理從不同的角度說明了同樣一個地質事實: 中國東部存在一個巨型深斷裂系統, 組成斷裂系統的各條斷裂之間有著深刻的內在聯系, 表現在物理上則是Moho面的局部隆起和規則重力異常分布, 表現在化學上則是由于深部物質上涌造成的化學元素的分區差異, 以及化學分區邊界與物理劃分邊界的一致性。

7　巨型深斷裂系統對金屬礦床分布的控制作用

黃宗理和嚴加永(2011)根據礦床密集分布地區的形狀將中國金屬礦成礦帶分為線狀成礦帶和面狀成礦區兩類(圖12)。線狀成礦帶長度可達上千km,寬度只有50 km至100 km, 面狀成礦區則比線狀成礦帶寬得多。線狀成礦帶多處于板塊邊界或縫合帶上, 地質上表現為深大斷裂, 有條帶狀剩余重力異常存在。面狀成礦區主要分布在在中國東部的東北,山東, 華南等地。為什么會產生面狀成礦區, 當時不太清楚。從本文分析來看, 其實中國東部的幾個面狀成礦區本質上就是一個——由九條斷裂組成的中國東部巨型深斷裂系統。就每一條斷裂而言, 其地質地球物理背景和線狀成礦帶是相同的, 成礦問題的實質都是斷裂切穿地殼, 深部幔源物質上涌,提供了成礦的物質源和能量源, 物理上表現為Moho隆起和重磁異常, 地質上表現為斷裂和巖體控礦。因為各條斷裂相距不遠, 幾條斷裂的綜合作用, 就表現為面狀成礦區了。①號北段和⑨號異常組成東北成礦區, ①號中段和⑦號異常組成山東成礦區, ②到⑤號異常組成長江中下游和華南成礦區。

目前已經在⑦號異常帶發現了山東巨型金礦帶, ⑤號異常帶上發現了著名的長江中下游成礦帶,④號異常帶上發現了德興等巨型銅礦, 在③號異常帶發現了紫金銅金礦集區。最近在依蘭—舒蘭斷裂和敦密斷裂控制的東北小興安嶺地區, 發現了一系列大型鉬礦和鎢礦, 在油氣勘查中也有很重要的新發現。從找礦角度看, 這9條異常都應當重視, 尤其是在異常的交會部位更值得關注。

常印佛和鄧晉福(2013)提出, 從大別到臺灣存在3個隆起區和3個凹陷區, 依次為大別隆起, 下揚子凹陷, 江南隆起, 浙贛凹陷, 武夷隆起, 永梅凹陷, 三條凹陷區中火成巖以幔源物質為主, 形成了一套親幔元素組合, 控制了區內大型鐵銅金礦集區。對照本文可見, 下揚子凹陷對應⑤號剩余重力異常異常, 浙贛凹陷對應④號重力異常, 永梅凹陷大體與②號③號重力異常的東段相對應。三隆三凹的本質其實是巨型深斷裂系統在江南部分的反映。

8　討論

根據數值模擬可以給出斷裂系統的形成模型:已經拼合為整體的歐亞大陸, 在中生代時受到菲律賓板塊和太平洋板塊近NNE方向構造力的作用,其東部邊緣地區首先沿薄弱地區斷開, 贛江斷裂被重新激活, 斷裂東部地塊向北滑移, 進而形成郯廬斷裂和依蘭—舒蘭斷裂, 構成一條縱貫南北的走滑斷層帶。處于自由狀態的斷層帶東盤, 在由南向北滑動中遞次產生一系列派生斷層, 派生斷層彼此平行, 間距大約等于300 km。在東盤整體向北移動的過程中, 發生過多期次的擠壓和伸展的交替過程,造成大量地幔物質沿斷裂帶上涌, 形成中國東部中生代以來的巨量火山巖漿活動。

目前一般認為郯廬斷裂向北穿過渤海進入了東北, 但是否向南穿過長江則看法不一。本文依據地球物理探測的衛星重力和深部地震資料, 認為郯廬斷裂是九條斷裂組成的巨型深斷裂系統的一部分,從整體分析長江以南的贛江等斷裂也屬于這個系統。

數值模擬顯示(圖10)在商南—榮成板塊縫合帶處, 地塊向北平移約550 km, 在鴨綠江斷裂帶處移動約650 km, 敦密斷裂處向北平移了約900 km。徐嘉瑋等(1980)通過對郯廬斷裂東西兩側地塊的時代對比, 得出的商南—榮成縫合帶、鴨綠江斷裂帶的地塊平移的距離分別也是500 km和700 km左右,說明數值模擬結果是正確的。敦密斷裂的平移距離目前還沒有得到觀測數據證實。根據數值模擬結果,地塊在東北長距離的移動, 將會產生如逆沖、地塊疊合等一系列復雜的地質現象, 這對最近在小興安嶺地區火山巖下發現油氣顯示, 對敦密斷裂和依蘭舒蘭斷裂之間殘存的老地塊的來源等問題, 也許都可以提出新的思路。

按照數值模擬的結果, 9條斷裂由南向北形成時代依次變新。一般認為大的斷裂都是長期、多次構造活動的累積結果。目前的地質研究認為贛江斷裂形式于古生代, 活躍于侏羅紀到白堊紀。郯廬斷裂開始形成于三疊紀244~209 Ma, 依蘭—舒蘭斷裂和敦化—密山斷裂形成于中侏羅世170~180 Ma。和數值模擬的結果總體一致。當然整個斷裂系統的形成演化過程及其動力學問題還需要更多地質工作的驗證。

9　結論

(1)衛星重力測量提供了包含陸域和海域統一網度和精度的重力異常, 對衛星布格重力異常進行分離后, 發現隱藏在區域異常背后的剩余重力異常呈現出規則分布的異常組合現象: 一條南北向異常貫穿整個中國大陸東部, 在南北向異常東側, 派生出9條北東向異常, 向東伸向大海, 這9條異常彼此平行, 間隔相等, 與南北向異常共同構成一個統一的系統。

(2)高精度反射地震測量結果表明, 剩余重力異常是深斷裂下方Moho破裂、幔源物質上涌造成的,深斷裂、Moho隆起和剩余重力異常三者密切相關,可根據重力異常對原來只在陸地上劃定的斷裂向海域和境外追蹤補充, 補充完整后的斷裂組合顯示中國東部存在一個巨型深斷裂系統: 主斷裂近南北向,南起廣州, 跨長江, 過渤海, 北出中國國境, 沿黑龍江下游入鄂霍次克海。斷裂北段和東側共有9條北東向斷裂, 北東向斷裂彼此平行, 間隔大致相等,向東延伸至大陸邊緣, 斷裂都是切穿Moho面的深斷裂。

(3)采用計算機數值模擬了巨型深斷裂系統的形成過程, 模擬表明歐亞板塊東緣, 在中生代時受到水平方向力的作用, 沿近南北向產生走滑斷裂,斷裂東盤自南向北滑動, 主斷裂的東側生成一系列派生斷層, 派生斷層基本平行, 間距大體等于300 km。每條斷裂都切穿了地殼, 地幔物質上涌,形成整個中國東部的巨型火山巖漿活動。

(4)巨型深斷裂系統控制著中國東部金屬礦產的分布, 構成東北、山東、華南三個面狀成礦區。9條異常上已經找到巨量金屬礦產, 深入分析9條異常及相互關系可以發現, 這里仍然存在著巨大的找礦潛力。

(5)巨型深斷裂系統在近代還有活動, 發生過1668年的郯城8.5級大地震、海城1975年7.3級地震、長白山的火山噴發等。巨型深斷裂系統地處太平洋板塊的西岸, 與太平洋東岸的圣安德烈斯斷層遙相呼應, 但遠比其宏大。由于中國人口最密集,經濟最發達的地區都處在這個巨型斷裂系統之上,需予以高度關注。

致謝: 感謝中國地質科學院地質研究所高銳研究員、熊小松副研究員, 中國地質科學院礦產資源研究所史大年研究員提供的地震測量資料, 感謝中國地質科學院董樹文研究員, 中國地質大學(北京)姚長利教授對本文寫作的支持和幫助, 感謝中國地質科學院礦產資源研究所汪杰碩士協助繪制了大部分插圖, 感謝兩位匿名審稿人提出的建設性意見。

Acknowledgements:

This study was supported by National Natural Science Foundation of China (Nos. 41574133 and 41104061) and the Science and Technology Project (No. SinProbe-08).

參考文獻:

陳毓川. 2007. 中國成礦體系與區域成礦評價[M]. 北京: 地質出版社: 32-67.

程裕琪. 1994. 中國區域地質概論[M]. 北京: 地質出版社: 75-83.

常印佛, 鄧晉福. 2013. 大別-臺灣走廊域成礦區帶形成的四維結構[M]. 北京: 地質出版社: 85-101.

黃宗理, 嚴加永. 2011. 中國剩余重力異常與金屬礦分布關系研究[J]. 地球學報, 32(6): 652-658.

侯遵澤, 楊文采. 1997. 中國重力異常的小波變換與多尺度分析[J]. 地球物理學報, 40(1): 85-96.

劉光鼎, 郝天姚, 劉伊克. 1997. 中國大地構造宏觀格架及其與礦產資源的關系[J]. 科學通報, 42(2): 113-118.

路曉翠, 張景發, 朱魯, 陳丁, 張鵬, 李麗梅. 2012. 利用小波多尺度分解研究郯廬斷裂帶蘇魯段構造[J]. 地球物理學進展, 27(1): 58-67.

馬文璞. 1992. 區域構造解析: 方法理論和中國板塊構造[M].北京: 地質出版社: 202-289.

馬宗晉, 高祥林, 宋正范. 2006. 中國布格重力異常水平梯度圖的判讀和構造解釋[J]. 地球物理學報, 49(1): 106-114.

毛建仁. 2013. 中國東南部及鄰區中新生代巖漿作用與成礦[M].北京: 科學出版社: 78-103.

史大年, 呂慶田, 徐文藝, 嚴加永, 趙金花, 董樹文, 常印佛. 2012. 長江中下游成礦帶及鄰區地殼結構-MASH成礦過程的P波接收函數成像證據?[J]. 地質學報, 286(3): 389-399.

滕吉文, 王謙身, 魏斯禹. 1982. 中國地殼構造的基本輪廓[J].石油物探, 2: 14-31.

涂廣紅, 江為為, 朱東英, 周立宏, 肖敦清, 高嘉瑞, 袁淑琴. 2006. 中國東北地區剩余重磁異常特征與地質構造及成礦帶的關系[J]. 地球物理學進展, 21(3): 746-755.

熊小松, 高銳, 李秋生, 盧占武, 王海燕, 李文輝, 管燁. 2009.深地震探測揭示的華南地區莫霍面深度[J]. 地球學報, 30(6): 774-786.

熊小松, 高銳, 張興洲, 李秋生, 侯賀晟. 2011. 深地震探測揭示的華北及東北地區莫霍面深度[J]. 地球學報, 32(1): 46-56.

張理剛, 王可法, 陳振勝, 劉敬秀, 李之彤. 1993. 中國東部中生代花崗巖長石鉛同位素組成與鉛同位素省劃分[J]. 科學通報, 03: 254-257.

朱日祥, 徐義剛, 朱光, 張宏福, 夏群科, 鄭天愉. 2010. 華北克拉通破壞[J]. 中國科學: 地球科學版, 42(8): 1135-1159.

References:

CHEN Yu-chuan. 2007. Chinese metallogenic system and regional metallogenic evaluation[M]. Beijing: Geological Publishing House: 32-67(in Chinese).

CHENG Yu-qi. 1994. Introduction to Chinese regional geology[M]. Beijing: Geological Publishing House: 75-83(in Chinese).

CHANG Yin-fo, DENG Jing-fu. 2013. The four-dimensional structure of Dabie Taiwan corridor zone mineralization zone[M]. Beijing: Geological Publishing House: 85-101(in Chinese).

HUANG Zong-li, YAN Jia-yong. 2011. The Relationship between the Residual Gravity Anomaly and the Distribution of Metallic Deposits in China[J]. Acta Geoscientica Sinica, 32(6): 652-658(in Chinese with English abstract).

HOU Zun-ze, YANG Wen-cai. 1997. Wavelet transform and multi-scale analysis on gravity anomalies of China[J]. Acta Geophysica Sinica, 40(1): 85-96(in Chinese with English abstract).

LIU Guang-ding, HAO Tian-yao, LIU Yi-ke. 1997. Tectonic Chinese macro framework and its relationship with the mineral resources[J]. Chinese Science Bulletin, 42(2): 113-118(in Chinese).

LU Xiao-cui, ZHANG Jing-fa, ZHU Lu, CHEN Ding, ZHANG Peng, LI Li-mei. 2012. Study on the structure in the Sulu segment of Tan-Lu fault zone by wavelet multi-scale decomposition[J]. Progress in Geophysics, 27(1): 58-67(in Chinese with English abstract).

MA Wen-pu. 1992. Regional structural analysis: Theory and Chinese plate tectonics[M]. Beijing: Geological Publishing House: 202-289(in Chinese).

MA Zong-jin, GAO Xiang-lin, SONG Zheng-fan. 2006. Analysis and tectonic interpretation to the horizonta-l gradient map calculated from Bouguer gravity data in the China mainland[J]. Chinese Journal of Geophysics, 49(1): 106-114(in Chinese with English abstract).

MAO Jian-ren. 2013. Cenozoic magmatism and mineralization China southeast and adjacent areas[M]. Beijing: Science Press: 78-103(in Chinese).

SHI Dia-nian, Lü Qing-tian, XU Wen-yi, YAN Jia-yong, ZHAO Jin-hua, DONG Shu-wen, CHANG Yin-fo. 2012. Crustal structures beneath the Mid-lower Yangtze metallogenic belt and its adjacent regions in east China - evidences from P-wave receiver function iaging for a MASH metallization process?[J]. Acta Geologica Sinica, 86(3): 389-399(in Chinese with English abstract).

TENG Ji-wen, WANG Qian-shen, WEI Si-yu. 1982. A general picture of the crust of China[J]. Geophysical Prospecting for Petroleum, 2: 14-31(in Chinese with English abstract).

TU Guang-hong, JIANG Wei-wei, ZHU Dong-ying, ZHOU Li -hong, XIAO Dun-qing, GAO Jia-rui, YUAN Shu-qin. 2006. The relationships between the characteristics of Northeast China residual gravity and magnetic anomalies and geological tectonic & metallogenic belt[J]. Progress in Geophysics, 21(3): 746-755(in Chinese with English abstract).

XIONG Xiao-song, GAO Rui, LI Qiu-sheng, LU Zhan-wu, WANG Hai-yan, LI Wen-hui, GUAN Ye. 2009. The Moho Depth of South China Revealed by Seismic Probing[J]. Acta Geoscientica Sinica, 30(6): 774-786(in Chinese with English abstract).

XIONG Xiao-song, GAO Rui, ZHANG Xing-zhou, LI Qiu-sheng, HOU He-sheng. 2011. The Moho Depth of North China and Northeast China Revealed by Seismic Detection[J]. Acta Geoscientica Sinica, 32(1): 46-56(in Chinese with English abstract).

ZHANG Li-gang, WANG Ke-fa, CEHN Zheng-sheng, LIU Jing-xiu, LI Zhi-tong. 1993. Mesozoic granite feldspar Pb isotopic composition of lead isotope provinces in east Chinese[J]. Chinese Science Bulletin, 03: 254-257(in Chinese).

ZHU Ri-xiang, XU Yi-gang, ZHU Guang, ZHANG Hong-fu, XIA Qun-ke, ZHENG Tian-yu. 2012. Destruction of the North China Craton[J]. Sci. China Earth Sci., 42(8): 1135-1159(in Chinese).

A Huge Deep Fault System at the East Edge of Eurasia: The New Tectonic Interpretation Based on Satellite Gravity

HUANG Zong-li1, 2), WANG Dian3), YAN Jia-yong1)*, ZHANG Huai4)
1) MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037; 2) Consulting & Research Center Ministry of Land and Resources, Beijing 100035; 3) Key Laboratory of Geo-detection, Ministry of Education, China University of Geosciences(Beijing), Beijing 100083; 4) Key Laboratory of Computational Geodynamics, Chinese Academy of Sciences, Beijing 100049

Abstract:Through the subtle processing of satellite gravity data, a regular distribution of gravity anomaly combination was found at the east edge of Eurasia. In order to study geological significance of the combination of these anomalies, the authors adopted gravity inversion based seismic constraints, and the result shows that there exist positive relationships between the discovered gravity anomaly, the uplift of the Moho and the deep fault. According to the gravity anomalies and geological information, the authors found a huge deep fault system at the east edge of Eurasia. The main fault stretches more than 3000 km from Guangzhou northward to the sea of Okhotsk. On the eastern side of the main fault zone, there are 9 parallel faults distributed approximatelybook=26,ebook=29equidistantly; these faults trend NE and extend to the continental margin into the sea. The authors used numerical simulation method to verify the formation of this fault system: under the horizontal action of force, SN-trending strike slip faults were generated in the east of Eurasian Continent, and then a series of transform fault deep into the lithosphere came into being. In the formation process of the fault system, a large quantity of mantle matter surged upward, resulted in the uplift of Moho, and then triggered massive volcanic- magmatic activities after Mesozoic in eastern China and controlled the distribution of endogenic metallic deposits. The fault system remains active at present, which led to the magnitude 8.5 earthquake in Tancheng and volcanic eruption in the Changbai Mountains.

Key words:satellite gravity anomaly; fault system; Eurasia

*通訊作者:嚴加永, 男, 1977年生。副研究員。主要從事深部礦產資源探測。

作者簡介:第一 黃宗理, 男, 1947年生。高級工程師?，F從事地球物理與成礦研究。

通訊地址:100037, 北京市西城區百萬莊大街26號。E-mail: zlhuang5653@hotmail.com。 100037, 北京市西城區百萬莊大街26號。E-mail: yanjy@163.com。

收稿日期:2015-05-11; 改回日期: 2015-09-03。責任編輯: 閆立娟。

中圖分類號:P631.123; P542.3

文獻標志碼:A

doi:10.3975/cagsb.2016.01.03