孫玉軍，吳中海，賈鳳琴

(1.中國地質科學院地質力學研究所，北京 100081；2.海南省地質環(huán)境監(jiān)測總站，?？?570216)

?

長江經(jīng)濟帶地區(qū)巖石圈熱流變結構和深部動力學

孫玉軍1，吳中海1，賈鳳琴2

(1.中國地質科學院地質力學研究所，北京 100081；2.海南省地質環(huán)境監(jiān)測總站，海口 570216)

通過Crust 2.0模型構建有限元三維數(shù)值模型，以地表觀測溫度、深部反演溫度和地表熱流作為約束，計算了長江經(jīng)濟帶地區(qū)巖石圈溫度結構；在溫度結構的基礎上，通過GPS觀測數(shù)據(jù)得到的地表應變率和選取代表性巖石物性，計算了長江經(jīng)濟帶巖石圈流變強度和等效粘滯性系數(shù)。結果表明：長江經(jīng)濟帶地區(qū)巖石圈表現(xiàn)出明顯的橫向不均勻性，其中四川盆地表現(xiàn)為低溫、高強度和高粘滯性的特征，相同深度，四川盆地核心區(qū)比周圍地塊的溫度低100～300 ℃，強度和粘滯性分別比周緣高1～2個數(shù)量級；在溫度、強度和粘滯性的過渡地帶構造活動性比較強烈。結合深部地球物理觀測綜合分析認為，長江經(jīng)濟帶東、西部分別受到太平洋板塊深俯沖和印度-歐亞大陸板塊碰撞的影響，其深部影響范圍分別達到四川盆地的東、西邊界，深部動力學過程可能導致了巖石圈的橫向不均勻性。

長江經(jīng)濟帶；巖石圈；熱流變結構

1 研究背景

長江經(jīng)濟帶地區(qū)主要是指長江干流及其支流流經(jīng)的區(qū)域，依托黃金水道建設長江經(jīng)濟帶是我國未來經(jīng)濟社會發(fā)展的重大戰(zhàn)略目標之一。從構造地質角度看，長江經(jīng)濟帶與華南大陸的大部分相重合，包含了揚子地塊的主體以及華夏地塊和華北地塊的部分地區(qū)，有著復雜的地質構造演化歷史[1～2](見圖1)。中新生代以來，該區(qū)處于青藏高原形成及擴展、太平洋板塊西向俯沖和印度-澳大利亞板塊北向差異運動的共同作用下，因此構造變形強烈，動力學演化過程復雜[1,3]。對長江經(jīng)濟帶地區(qū)的巖石圈熱流變結構和深部地球物理特征進行研究不僅有助于理解該區(qū)域的動力學演化機制，同時可以對新構造和活動構造調查提供深部約束。

箭頭代表相對于穩(wěn)定的歐亞大陸板塊，該區(qū)域及周緣的GPS觀測速度值[4～5]，紅色實線代表活動斷裂[6]圖1 長江經(jīng)濟帶地區(qū)構造背景圖(據(jù)文獻[4～6]修編)Fig.1 Tectonic setting map around the Yangtze River Economic Belt

為獲取該區(qū)域的巖石圈熱流變結構，本研究基于Crust 2.0模型構建了合理的三維數(shù)值模型，以地表觀測溫度、深部反演溫度和地表熱流作為約束計算長江經(jīng)濟帶巖石圈的深部溫度結構，并在此基礎上，通過GPS觀測數(shù)據(jù)得到的地表應變率，選取不同深度代表性巖石物性，計算了長江經(jīng)濟帶地區(qū)巖石圈流變強度和等效粘滯性系數(shù)?；谶@些結果結合深部地球物理研究結果對長江經(jīng)濟帶地區(qū)的巖石圈特征及深部動力學過程進行了綜合分析。

2 計算方法

2.1 溫度場

根據(jù)Crust 2.0模型和地表地形數(shù)據(jù)，在球坐標系下建立了中國大陸巖石圈三維有限元模型。通過地震波波速獲取的巖石圈上地幔溫度為計算巖石圈熱結構提供了強有力的深部約束[7～8]。采用此方法，An and Shi[9～10]計算了中國大陸及鄰區(qū)巖石圈上地幔溫度。本文選取他們在100 km深度的溫度結果作為計算巖石圈熱結構的深部約束，同時將中國大陸地區(qū)195個地表氣象臺站提供的地表地溫數(shù)據(jù)作為地表的溫度約束。根據(jù)三維穩(wěn)態(tài)熱傳導方程，考慮到熱參數(shù)的各向不均勻性，計算了中國大陸及鄰區(qū)巖石圈三維熱結構。該計算結果與地表熱流觀測數(shù)據(jù)進行了對比，其中89%的數(shù)據(jù)誤差控制在20%以內(熱流實際觀測的誤差為20%[11])，本文中計算巖石圈蠕變強度中需要的溫度采用此結果，詳見文獻[12]。

2.2 巖石圈強度和等效粘滯性系數(shù)

巖石圈強度和等效粘滯性系數(shù)是巖石圈流變結構的2個基本參數(shù)。一般來說，巖石的失效形式有摩擦滑動、脆性破裂和蠕變，而這3種形式也被用來定義巖石圈的強度[13]。

①摩擦滑動：假定2σ2=σ1-σ3，孔隙壓力為0，本文采用走滑斷層的強度表達式來計算巖石圈的摩擦滑動強度[14～16]。

(1)

式中：τf為摩擦滑動強度；z為深度，km；σ1、σ2和σ3和分別為最大、中間和最小主應力，MPa。

②脆性破裂：如果所考慮的巖石預先并不存在破裂面或者在破裂面之間完整巖石的尺寸比較大，則巖石的整體強度需要用脆性破裂來約束[13]。這種強度機制與巖石樣品的尺度、所處的溫度、壓力和應變率密切相關。在溫度為800 ℃以下可以表示為如下形式[15～18]：

(2)

③蠕變：當巖石圈溫度和壓力足夠高時，蠕變往往占主導地位。蠕變強度與溫度、應變率和巖石屬性密切相關。蠕變強度可以采用下式計算[19～21]：

(3)

在一定溫度和壓力下，以上3種強度機制中的一種將會占主導地位。本文在計算巖石圈強度的過程中，3種強度機制都進行了計算，在一定溫度和壓力條件下，選取其中最小的一個作為巖石圈的強度，即：

(4)

得到巖石圈的強度后，巖石圈的等效粘滯性系數(shù)可以表達為如下形式[22]：

(5)

在溫度計算[12]的基礎上，通過GPS觀測數(shù)據(jù)得到的地表應變率[23]，選取不同深度代表性巖石物性[24]，計算得到了長江經(jīng)濟帶巖石圈流變強度和等效粘滯性系數(shù)。

3 長江經(jīng)濟帶地區(qū)巖石圈熱流變結構

3.1 長江經(jīng)濟帶地區(qū)巖石圈熱結構特征

不同深度的溫度結果(見圖2)顯示，長江經(jīng)濟帶地區(qū)的四川盆地在相同深度的溫度較低，與之相鄰的松潘—甘孜地塊、印支地塊、華夏地塊和中朝地塊在相同深度溫度均較高。從具體溫度數(shù)值來看，在相同深度，四川盆地核心區(qū)比周圍地塊的溫度低100～300 ℃。溫度場結果顯示，四川盆地巖石圈呈現(xiàn)冷而厚的特征，這使其構造變形弱，相對穩(wěn)定；而在四川盆地與其他地塊的過渡地帶，構造變形均較強烈，這與地表活動構造分布基本一致。從垂向剖面(見圖3a)看，四川盆地的巖石圈下部(大于80 km)溫度整體較低，而處于其兩側的青藏高原和華夏地塊巖石圈下部溫度較高。長江經(jīng)濟帶下游的華夏地塊和中朝地塊的深部高溫區(qū)可能與太平洋板塊俯沖引起的巖石圈減薄因素有關，而上游的印支地塊可能與印度板塊的俯沖有關，松潘—甘孜地塊的高溫區(qū)可能與青藏高原物質東向運移有關。因此，該結果顯示揚子地塊的東、西兩部分溫度分布特征差異較大。

圖2 長江經(jīng)濟帶地區(qū)巖石圈不同深度溫度場Fig.2 Temperature distribution at different depths of the lithosphere in the Yangtze River Economic Belt

圖3 長江經(jīng)濟帶地區(qū)各個剖面巖石圈溫度、強度和等效粘滯性系數(shù)分布(剖面位置見圖1)Fig.3 Distributions of temperature, strength and viscosity of the lithosphere at different profiles in the Yangtze River Economic Belt

3.2 長江經(jīng)濟帶地區(qū)巖石圈流變結構特征

巖石圈的強度顯示出明顯的橫向不均勻性(見圖4)。在淺部(深度小于60 km)，青藏高原周緣的巖石圈強度較低(小于10 MPa)，揚子地塊、華夏地塊和中朝地塊巖石圈強度都較高(大于10 MPa)，而四川盆地及周緣地區(qū)的巖石圈強度最高，基本處于100 MPa左右。在深部(深度大于60 km)，只有四川盆地及周緣地區(qū)巖石圈強度高于1 MPa，其他地區(qū)均低于該值。不同地塊之間存在比較明顯的強度過渡帶，特別是在40 km深度處青藏高原和四川盆地之間的過渡帶比較明顯，強度差異大(見圖4b)，同時該過渡帶與重力梯度帶、地形階梯帶、地震波速變化帶和構造活動性的過渡帶也相對一致。垂向上看(見圖3b)，青藏高原的巖石圈強度整體較低，即相比周緣地塊，其地殼和巖石圈地幔強度均較低，而且在60 km以下基本都低于1 MPa；在四川盆地及周緣100 km以上巖石圈強度基本都高于1 MPa，地殼和地幔巖石圈頂部強度最高，在10～100 MPa；其他區(qū)域地殼和地幔巖石圈頂部強度也較高，但巖石圈地幔深部(大于80 km)，巖石圈強度明顯降低至0.1 MPa左右。這與這些地區(qū)的溫度分布特征比較相似。計算得到的強度剖面也顯示了明顯的分層性(見圖3b)，巖石圈強度呈現(xiàn)出脆-韌性相間分布的特征，下地殼均呈現(xiàn)出韌性特征，而且弱強度的下地殼廣泛存在于研究區(qū)，下地殼的流變強度一般比上地殼和巖石圈地幔上部低1～2個數(shù)量級。

圖4 長江經(jīng)濟帶地區(qū)巖石圈不同深度的強度分布Fig.4 Lithospheric strength distribution at different depths in the Yangtze River Economic Belt

巖石圈等效粘滯性分布特征與巖石圈強度的分布特征比較一致。橫向上看(見圖5)，在淺部(深度小于60 km)，青藏高原的等效粘滯性系數(shù)(一般低于1023Pa·s)比揚子地塊、華夏地塊和中朝地塊(一般高于1023Pa·s)要低；在深部(深度大于60 km)，只有四川盆地及周緣巖石圈強度高于1022Pa·s，其他地區(qū)均低于該值。不同地塊之間粘滯性系數(shù)也存在比較明顯的過渡帶。垂向上看(見圖3c)，上地殼、中地殼和巖石圈地幔上部的等效粘滯性系數(shù)較高，下地殼的等效粘滯性系數(shù)較低，厚而“弱”的下地殼在青藏高原的特征最為明顯；從量級上看(見圖3c)，下地殼的等效粘滯性系數(shù)一般比上地殼和巖石圈地幔的要低1～2個數(shù)量級。青藏高原的上地殼和中地殼等效粘滯性系數(shù)在1021～1022量級，而揚子地塊、華夏地塊和中朝地塊的上地殼和中地殼的等效粘滯性系數(shù)在1023～1024量級；青藏高原的下地殼等效粘滯性系數(shù)只有1019～1021Pa·s，而揚子地塊、華夏地塊和中朝地塊的下地殼等效粘滯性系數(shù)在1021～1022量級。

3.3 長江經(jīng)濟帶地區(qū)深部地球物理特征

從長江經(jīng)濟帶地區(qū)的深部地震層析成像結果(見圖6)來看，所獲取的熱流變結構與其有較好的一致性。高精度P波層析成像結果[25～26]顯示，揚子地塊的東、西部有較大差異。在揚子地塊西部，主要以四川盆地的低速異常為主，巖石圈厚而穩(wěn)定，受其西部印度板塊俯沖的影響，四川盆地下部的軟流圈表現(xiàn)為低速異常，但是在400 km深度之上整體為低速異常，與本文所獲取的溫度場、巖石圈強度和等效粘滯性比較一致，因此，印度板塊俯沖對揚子地塊的影響可能只限于四川盆地及以西地區(qū)。在揚子地塊中東部，除了揚子地塊東部100 km深度以上部分表現(xiàn)為高速異常外，其他400 km深度以上區(qū)域均表現(xiàn)為明顯的低速異常，這可能是由于太平洋板塊深俯沖引起的，地幔深部層析成像結果顯示深俯沖的太平洋板塊在410～660 km深度范圍內可能一直延伸至四川盆地以東區(qū)域，使得其上覆地幔加熱擾動而呈現(xiàn)低速異常，該結果顯示太平洋板塊的深俯沖可能向西影響到四川盆地的東邊界。從本文所獲取的巖石圈熱結構和流變特征來看，與深部地震層析成像結果比較一致。

圖6 長江經(jīng)濟帶地區(qū)及周緣深部地震層析成像結果[25](剖面沿著北緯30°)Fig.6 Deep seismic tomography results in the Yangtze River Economic Belt and adjacent area

4 討論和結論

本文通過構建三維數(shù)值模型，結合地表溫度觀測、熱流觀測和深部巖石圈地幔溫度約束，獲取了長江經(jīng)濟帶地區(qū)巖石圈溫度場結果；在此基礎上根據(jù)巖石圈分層結構和實驗室獲取的典型代表性巖石的流變參數(shù)，依據(jù)合理的流變強度機制約束，計算分析了長江經(jīng)濟帶地區(qū)巖石圈強度和等效粘滯性系數(shù)。所獲取的巖石圈溫度、強度和粘滯性系數(shù)雖然存在一定的不確定性[12,24]，但與深部地震層析成像獲取的結果具有較好的一致性。

相對于青藏高原，揚子地塊和華夏地塊整體表現(xiàn)為低溫、高強度和高粘滯性的特征。但揚子地塊東部和西部在溫度、巖石圈強度和等效粘滯性系數(shù)等特征上存在較大差別，四川盆地及其周緣地殼和地幔巖石圈整體表現(xiàn)為低溫、高強度和高粘滯性的特征，構造穩(wěn)定變形弱；而揚子地塊中、東部地區(qū)表現(xiàn)為高溫、低強度和低粘滯性的特征。結合深部地震波層析成像結果初步認為，印度-歐亞板塊的碰撞動力作用可能僅僅影響到四川盆地以西地區(qū)，而太平洋板塊的深俯沖影響可能一直延伸至四川盆地的東邊界，長江經(jīng)濟帶地區(qū)東、西方向受到這兩大動力作用影響，在溫度、巖石圈強度和粘滯性方面表現(xiàn)出較大差異性。

[1] 舒良樹．華南構造演化的基本特征[J]．地質通報，2012，31(7)：1035～1053.

SHU Liang-shu. An analysis of principal features of tectonic evolution in South China Block[J]. Geological Bulletin of China, 2012, 31(7): 1035～1053.

[2] 張國偉，郭安林，王岳軍，等．中國華南大陸構造與問題[J]．中國科學：地球科學，2013，43(10)：1553～1582.

ZHANG Guo-wei, GUO An-lin, WANG Yue-jun, et al. Tectonics of South China continent and its implications[J]. Science China: Earth Sciences, 2013, 43(10): 1553～1582.

[3] 張岳橋，徐先兵，賈東，等．華南早中生代從印支期碰撞構造體系向燕山期俯沖構造體系轉換的形變記錄[J]．地學前緣，2009，16(1)：234～247.

ZHANG Yue-qiao, XU Xian-bing, JIA Dong, et al. Deformation record of the change from Indosinian collision-related tectonic system to Yanshanian subduction-related tectonic systemin South China during the Early Mesozoic[J]. Earth Sciences Frontiers, 2009, 16(1): 234～247.

[4] Gan W, Zhang P, Shen Z-K, et al. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J]. Journal of Geophysical Research, 2007, 112: B08416.

[5] Kreemer C, Blewitt G, Klein E C. A geodetic plate motion and Global Strain Rate Model[J]. Geochemistry, Geophysics, Geosystems, 2014, 15(10): 3849～3889.

[6] Deng Q, Zhang P, Ran Y, et al. Basic characteristics of active tectonics of China[J]. Science in China: Earth Sciences, 2003, 46(4): 356～372.

[7] Goes S, Govers R, Vacher P. Shallow mantle temperatures under Europe from P and S wave tomography[J]. Journal of Geophysical Research, 2000, 105(B5): 11153～11169.

[8] Goes S, Lee S v d. Thermal structure of the North American uppermost mantle inferred from seismic tomography[J]. Journal of Geophysical Research, 2002, 107(B3): 2050.

[9] An M, Shi Y. Lithospheric thickness of the Chinese continent[J]. Physics of the Earth and Planetary Interiors, 2006, 159(3/4): 257～266.

[10] An M, Shi Y. Three-dimensional thermal structure of the Chinese continental crust and upper mantle[J]. Science in China: Series D, 2007, 50(10): 1441～1451.

[11] Hu S, He L, Wang J. Heat flow in the continental area of China: a new data set[J]. Earth and Planetary Science Letters, 2000, 179: 407～419.

[12] Sun Y, Dong S, Zhang H, et al. 3D thermal structure of the continental lithosphere beneath China and adjacent regions[J]. Journal of Asian Earth Sciences, 2013, 62: 697～704.

[13] Kohlstedt D L, Evans B, Mackwell S J. Strength of the lithosphere: Constraints imposed by laboratory experiments[J]. Journal of Geophysical Research, 1995, 100(B9): 17587～17602.

[14] 臧紹先，李昶，魏榮強．巖石圈流變機制的確定及影響巖石圈流變強度的因素[J]．地球物理學進展，2002，17(1)：50～60.

ZANG Shao-xian, LI Chang, WEI Rong-qiang. The determination of the deformational mechanism of the lithosphere and the factors influenced the lithospheric rheological strength[J]. Progress in Geophysics, 2002, 17(1): 50～60.

[15] Zang S X, Wei R Q, Liu Y G. Three-dimensional rheological structure of the lithosphere in the Ordos block and its adjacent area[J]. Geophysical Journal International, 2005, 163(1): 339～356.

[16] Zang S X, Wei R Q, Ning J Y. Effect of brittle fracture on the rheological structure of the lithosphere and its application in the Ordos[J]. Tectonophysics, 2007, 429: 267～285.

[17] 魏榮強，臧紹先．巖石破裂強度的溫度和應變率效應及其對巖石圈流變結構的影響[J]．地球物理學報，2006，49(6)：1730～1737.

WEI Rong-qiang, ZANG Shao-xian. Effects of temperature and strain rate on the fracture strength of rock and their influences on the rheological structure of the lithosphere[J]. Chinese Journal of Geophysics, 2006, 49(6): 1730～1737.

[18] Shimada M. Lithosphere strength inferred from fracture strength of rocks at high confining pressures and temperatures[J]. Tectonophysics, 1993, 217: 55～64.

[19] Kirby S H. Rheology of the lithosphere[J]. Reviews of Geophysics and Space Physics, 1983, 21(6): 1458～1487.

[20] Ranalli G, Murphy D C. Rheological stratification of the lithosphere[J]. Tectonophysics, 1987, 132: 281～295.

[21] Weertman J. The creep strength of the Earth′s mantle[J]. Reviews of Geophysics and Space Physics, 1970, 8(1): 145～168.

[22] Kirby S H, Kronenberg A K. Rheology of the lithosphere: Selected topics[J]. Review of Geophysics, 1987, 25(6): 1219～1244.

[23] Zhu S, Shi Y. Estimation of GPS strain rate and its error analysis in the Chinese continent[J]. Journal of Asian Earth Sciences, 2010, 40(1): 351～362.

[24] 孫玉軍，董樹文，范桃園，等．中國大陸及鄰區(qū)巖石圈三維流變結構[J]．地球物理學報，2013，56(9)：2936～2946.

SUN Yu-jun, DONG Shu-wen, FAN Tao-yuan, et al. 3D rheological structure of the continental lithosphere beneath China and adjacent regions[J]. Chinese Journal of Geophysics, 2013, 56(5): 546～558.

[25] Huang J, Zhao D. High-resolution mantle tomography of China and surrounding regions[J]. Journal of Geophysical Research, 2006, 111(B9): B09305.

[26] Wei W, Xu J, Zhao D, et al. East Asia mantle tomography: New insight into plate subduction and intraplate volcanism[J]. Journal of Asian Earth Sciences, 2012, 60: 88～103.

LITHOSPHERIC THERMAL-RHEOLOGICAL STRUCTURE AND DEEP GEODYNAMICS IN THE YANGTZE RIVER ECONOMIC BELT

SUN Yu-jun1, WU Zhong-hai1, JIA Feng-qin2

(1.Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China;2.HainanProvinceGeologicalEnvironmentalMonitoringStation,Haikou570216,China)

Constructing the finite three-dimensional numerical model through Crust 2.0, under constraints of the surface temperature, deep inversion temperature and surface heat flow, we have calculated the thermal structure of the continental lithosphere beneath Yangtze River Economic Belt and adjacent regions. On the basis of 3D thermal structure, deriving surface strain rates from GPS observation, and selecting representative rock physical properties, we have calculated the strength and effective viscosity of the continental lithosphere. The results show that there is obvious lateral heterogeneity in the Yangtze River Economic Belt. The Sichuan basin has low temperature, high strength and high viscosity. In the same depth, temperature beneath the Sichuan basin is lower than ambient region for 100～300 ℃. Strength and viscosity is lower for 1～2 orders of magnitude. The transition zone of temperature, strength and viscosity has strong active tectonics. Combining with deep geophysical observation, the study proposed that the Pacific Plate subduction and India-Eurasia Plate collision played important roles in the deep dynamic process respectively in the eastern and western parts of the Yangtze River Economic Belt. The effective sphere reach the eastern and western boundaries of Sichuan basin, deep dynamic process may lead to the horizontal heterpgeneity of the lithosphere.

the Yangtze River Economic Belt; lithosphere; thermal-rheological structure

1006-6616(2016)03-0421-09

2016-03-06

中國地質調查局地質調查項目“長江經(jīng)濟帶活動構造與區(qū)域地殼穩(wěn)定性調查”(DD20160268)

孫玉軍(1983-)，男，博士，副研究員，主要從事地球動力學方面的研究。E-mail：sunyujunabc@163.com

P542+.5

A