黃宗理, 王　典, 嚴(yán)加永, 張　懷

1)中國(guó)地質(zhì)科學(xué)院礦產(chǎn)資源研究所, 國(guó)土資源部成礦作用與資源評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室, 北京 100037; 2)國(guó)土資源部咨詢研究中心, 北京 100035; 3)中國(guó)地質(zhì)大學(xué)(北京), 地下信息探測(cè)技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100083; 4)中國(guó)科學(xué)院計(jì)算地球動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京 100049

?

歐亞大陸東緣存在一個(gè)巨型深斷裂系統(tǒng)——基于衛(wèi)星重力的新發(fā)現(xiàn)

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

1)中國(guó)地質(zhì)科學(xué)院礦產(chǎn)資源研究所, 國(guó)土資源部成礦作用與資源評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室, 北京 100037; 2)國(guó)土資源部咨詢研究中心, 北京 100035; 3)中國(guó)地質(zhì)大學(xué)(北京), 地下信息探測(cè)技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100083; 4)中國(guó)科學(xué)院計(jì)算地球動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京 100049

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

關(guān)鍵詞:衛(wèi)星重力異常; 斷裂系統(tǒng); 歐亞大陸

本文由國(guó)家自然科學(xué)基金項(xiàng)目(編號(hào): 41574133; 41104061)和國(guó)家深部探測(cè)技術(shù)與實(shí)驗(yàn)研究專項(xiàng)(編號(hào): SinProbe-08)聯(lián)合資助。

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

我們?cè)谘芯繗W亞大陸東部地區(qū)的衛(wèi)星重力場(chǎng)時(shí),對(duì)布格重力異常進(jìn)行了區(qū)域場(chǎng)和局部重力異常的分離。發(fā)現(xiàn)在一定尺度的空間域?yàn)V波或一定頻帶寬度的頻率域?yàn)V波后, 得出的剩余重力異常在歐亞大陸東部邊緣地區(qū)出現(xiàn)了一個(gè)彼此平行、間距相等的規(guī)則組合。進(jìn)一步的研究發(fā)現(xiàn), 這些剩余重力異常在中國(guó)大陸上的部分與已知的地質(zhì)斷裂完全重合。

1　數(shù)據(jù)處理

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

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

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

2　剩余重力異常分布特征及其與地質(zhì)斷裂的對(duì)應(yīng)關(guān)系

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

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

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

①號(hào)重力異常南段與贛江斷裂, 中段與郯城廬江斷裂, 北段與依蘭—舒蘭斷裂相對(duì)應(yīng);

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

圖2　衛(wèi)星重力數(shù)據(jù)與地面重力數(shù)據(jù)對(duì)比圖Fig. 2　Comparison between satellite gravity and ground gravity

②號(hào)重力異常與廣州北東向斷裂對(duì)應(yīng);

③號(hào)重力異常與邵武—河源斷裂相對(duì)應(yīng);

④號(hào)重力異常與紹興—萍鄉(xiāng), 景德鎮(zhèn)—宜豐斷裂帶相對(duì)應(yīng);

⑥號(hào)重力異常與揚(yáng)子陸塊北緣斷裂帶相對(duì)應(yīng);

⑦號(hào)重力異常與商南—榮成板塊縫合帶相對(duì)應(yīng);

⑧號(hào)重力異常與鴨綠江斷裂相對(duì)應(yīng);

⑨號(hào)重力異常對(duì)應(yīng)敦化—密山斷裂。

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

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

9條重力異常中, 有8條在中國(guó)大陸境內(nèi)都和已知斷裂吻合得較好, 只有⑤號(hào)異常目前沒有已知斷裂與之對(duì)應(yīng)。⑤號(hào)異常從南京沿長(zhǎng)江入海。已有人根據(jù)該區(qū)各種地球物理信息, 指出長(zhǎng)江下游至泰州為一隱伏深斷裂帶(陳毓川, 2007)(暫稱為長(zhǎng)江下游隱伏斷裂)。石油地質(zhì)工作也已在⑤號(hào)異常海域部分發(fā)現(xiàn)了一組斷裂(馬文璞, 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條重力異常與已知斷裂的嚴(yán)格對(duì)應(yīng)說(shuō)明在這些斷裂下方有高密度物體沿著斷裂分布。問題在于這種綿延數(shù)千千米的條帶狀高密度物體是什么？SinoProbe-02專項(xiàng)最近完成的東北和華南兩條深地震反射剖面測(cè)量, 獲得了包含Moho面的地殼的精細(xì)圖像, 給出了對(duì)重力異常的解釋結(jié)果。

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

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

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

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

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

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

依照剩余重力異常與斷裂的對(duì)應(yīng)關(guān)系, 可將在中國(guó)大陸上地質(zhì)已經(jīng)劃定的斷裂沿剩余重力異常向中國(guó)境外和海域自然延伸, 把已知的斷裂補(bǔ)充完整。

①號(hào)異常北段對(duì)應(yīng)的依蘭—舒蘭斷裂帶, 從剩余重力異?？磻?yīng)向北延伸出境, 沿黑龍江下游進(jìn)入鄂霍次克海。①號(hào)異常從南端的廣州到北端進(jìn)入鄂霍次克海, 整體是連續(xù)的, 說(shuō)明①號(hào)異常所對(duì)應(yīng)的贛江斷裂帶, 郯廬斷裂帶, 依蘭—舒蘭斷裂帶, 從南到北整體上是連續(xù)的。這個(gè)連續(xù)的斷裂帶在衛(wèi)星地形圖上也可清晰識(shí)別。

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

④號(hào)異常對(duì)應(yīng)的紹興—萍鄉(xiāng), 景德鎮(zhèn)—宜春斷裂帶, ⑤號(hào)異常對(duì)應(yīng)的長(zhǎng)江下游隱伏斷裂帶, ⑥號(hào)異常對(duì)應(yīng)的揚(yáng)子陸塊北緣斷裂帶, ⑦號(hào)異常對(duì)應(yīng)的商南—榮成板塊縫合帶應(yīng)當(dāng)向東延伸進(jìn)入大海, 直至大陸邊緣。

⑧號(hào)異常對(duì)應(yīng)的鴨綠江斷裂的主體應(yīng)在沿遼東半島海岸線的海里, 斷裂的東段沿中朝邊界并穿過長(zhǎng)白山火山地區(qū)。

⑨號(hào)異常與敦化—密山斷裂完全對(duì)應(yīng)。

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

圖8　數(shù)值模擬模型參數(shù)圖Fig. 8　Parameters of numerical simulation model

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

圖10　G點(diǎn)位置圖Fig. 10　10 G point location map

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

5　巨型深斷裂帶系統(tǒng)的數(shù)值模擬

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

圖12　剩余重力異常與成礦區(qū)帶圖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°, 模擬結(jié)果是: 在南北向力的作用下, 先生成了廣州北東向的斷裂,隨著繼續(xù)向北推擠, 相繼形成了邵武—河源斷裂、紹興—萍鄉(xiāng)和景德鎮(zhèn)—宜豐斷裂、推測(cè)的長(zhǎng)江下游隱伏斷裂、揚(yáng)子陸塊北緣斷裂, 商南—榮成縫合帶的再次激活, 鴨綠江斷裂, 隨后斷層在后盤的近端再次激活, 然后在前緣再次生成敦化—密山斷裂以及最后一個(gè)依蘭—舒蘭斷裂的北段, 斷層形成過程見圖9。

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

6　巨型深斷裂系統(tǒng)上的地球化學(xué)現(xiàn)象

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

7　巨型深斷裂系統(tǒng)對(duì)金屬礦床分布的控制作用

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

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

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

8　討論

根據(jù)數(shù)值模擬可以給出斷裂系統(tǒng)的形成模型:已經(jīng)拼合為整體的歐亞大陸, 在中生代時(shí)受到菲律賓板塊和太平洋板塊近NNE方向構(gòu)造力的作用,其東部邊緣地區(qū)首先沿薄弱地區(qū)斷開, 贛江斷裂被重新激活, 斷裂東部地塊向北滑移, 進(jìn)而形成郯廬斷裂和依蘭—舒蘭斷裂, 構(gòu)成一條縱貫?zāi)媳钡淖呋瑪鄬訋АＬ幱谧杂蔂顟B(tài)的斷層帶東盤, 在由南向北滑動(dòng)中遞次產(chǎn)生一系列派生斷層, 派生斷層彼此平行, 間距大約等于300 km。在東盤整體向北移動(dòng)的過程中, 發(fā)生過多期次的擠壓和伸展的交替過程,造成大量地幔物質(zhì)沿?cái)嗔褞嫌? 形成中國(guó)東部中生代以來(lái)的巨量火山巖漿活動(dòng)。

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

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

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

9　結(jié)論

(1)衛(wèi)星重力測(cè)量提供了包含陸域和海域統(tǒng)一網(wǎng)度和精度的重力異常, 對(duì)衛(wèi)星布格重力異常進(jìn)行分離后, 發(fā)現(xiàn)隱藏在區(qū)域異常背后的剩余重力異常呈現(xiàn)出規(guī)則分布的異常組合現(xiàn)象: 一條南北向異常貫穿整個(gè)中國(guó)大陸東部, 在南北向異常東側(cè), 派生出9條北東向異常, 向東伸向大海, 這9條異常彼此平行, 間隔相等, 與南北向異常共同構(gòu)成一個(gè)統(tǒng)一的系統(tǒng)。

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

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

(4)巨型深斷裂系統(tǒng)控制著中國(guó)東部金屬礦產(chǎn)的分布, 構(gòu)成東北、山東、華南三個(gè)面狀成礦區(qū)。9條異常上已經(jīng)找到巨量金屬礦產(chǎn), 深入分析9條異常及相互關(guān)系可以發(fā)現(xiàn), 這里仍然存在著巨大的找礦潛力。

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

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

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).

參考文獻(xiàn):

陳毓川. 2007. 中國(guó)成礦體系與區(qū)域成礦評(píng)價(jià)[M]. 北京: 地質(zhì)出版社: 32-67.

程裕琪. 1994. 中國(guó)區(qū)域地質(zhì)概論[M]. 北京: 地質(zhì)出版社: 75-83.

常印佛, 鄧晉福. 2013. 大別-臺(tái)灣走廊域成礦區(qū)帶形成的四維結(jié)構(gòu)[M]. 北京: 地質(zhì)出版社: 85-101.

黃宗理, 嚴(yán)加永. 2011. 中國(guó)剩余重力異常與金屬礦分布關(guān)系研究[J]. 地球?qū)W報(bào), 32(6): 652-658.

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

劉光鼎, 郝天姚, 劉伊克. 1997. 中國(guó)大地構(gòu)造宏觀格架及其與礦產(chǎn)資源的關(guān)系[J]. 科學(xué)通報(bào), 42(2): 113-118.

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

馬文璞. 1992. 區(qū)域構(gòu)造解析: 方法理論和中國(guó)板塊構(gòu)造[M].北京: 地質(zhì)出版社: 202-289.

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

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

史大年, 呂慶田, 徐文藝, 嚴(yán)加永, 趙金花, 董樹文, 常印佛. 2012. 長(zhǎng)江中下游成礦帶及鄰區(qū)地殼結(jié)構(gòu)-MASH成礦過程的P波接收函數(shù)成像證據(jù)?[J]. 地質(zhì)學(xué)報(bào), 286(3): 389-399.

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

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

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

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

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

朱日祥, 徐義剛, 朱光, 張宏福, 夏群科, 鄭天愉. 2010. 華北克拉通破壞[J]. 中國(guó)科學(xué): 地球科學(xué)版, 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

*通訊作者:嚴(yán)加永, 男, 1977年生。副研究員。主要從事深部礦產(chǎn)資源探測(cè)。

作者簡(jiǎn)介:第一 黃宗理, 男, 1947年生。高級(jí)工程師?，F(xiàn)從事地球物理與成礦研究。

通訊地址:100037, 北京市西城區(qū)百萬(wàn)莊大街26號(hào)。E-mail: zlhuang5653@hotmail.com。 100037, 北京市西城區(qū)百萬(wàn)莊大街26號(hào)。E-mail: yanjy@163.com。

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

中圖分類號(hào):P631.123; P542.3

文獻(xiàn)標(biāo)志碼:A

doi:10.3975/cagsb.2016.01.03