中國大陸現(xiàn)今地應力場黏彈性球殼數值模擬綜合研究

范桃園，龍長興，楊振宇，陳群策，吳中海，邵兆剛，仝亞博

國土資源部新構造運動與地質災害重點實驗室，中國地質科學院地質力學研究所，北京 100081

中國大陸現(xiàn)今地應力場黏彈性球殼數值模擬綜合研究

范桃園，龍長興，楊振宇，陳群策，吳中海，邵兆剛，仝亞博

國土資源部新構造運動與地質災害重點實驗室，中國地質科學院地質力學研究所，北京 100081

中國大陸現(xiàn)今實測地應力場的狀態(tài)與板塊構造環(huán)境、活動斷裂帶分布、地形地貌以及地殼結構呈現(xiàn)一定相關性.在中國大陸西緣，印度洋板塊與歐亞板塊陸發(fā)生陸碰撞，在中國大陸東緣，菲律賓海板塊、太平洋板塊俯沖到歐亞板塊之下.中國大陸內部被大型活動斷裂帶分割為多個塊體，各個塊體的地殼結構和厚度呈不均勻分布，地形地貌起伏具有很大的差異.筆者以中國大陸塊體模型為基礎，把板塊構造作用和重力勢作為主要影響地應力狀態(tài)的兩個主要要素，在現(xiàn)今活動構造、GPS和實測地應力等成果的約束下，利用線性黏彈體球殼有限元模擬分析了中國大陸現(xiàn)今地應力場的分布特征和控制因素.結果表明：（1）構造應力場總體上呈現(xiàn)出西部擠壓，東部拉張的特征，印度板塊與歐亞板塊的持續(xù)碰撞形成了青藏高原及其周緣的擠壓性質的構造應力場，而東部菲律賓板塊與太平洋板塊的俯沖形成了黃海、東海和環(huán)渤海區(qū)域的拉張性質的構造應力場，中間為拉張環(huán)境和擠壓環(huán)境的過渡，最大主應力的方向受到板塊構造環(huán)境和活動構造分布的控制；（2）重力的影響主要體現(xiàn)在地形梯度大和地殼厚度結構變化大的地殼淺部區(qū)域，在藏南、滇西北局部地區(qū)的地殼淺部由于受到重力勢控制，呈現(xiàn)為張性應力場，在塔里木地區(qū)由于重力勢引起的應力場與構造應力場同為擠壓性質，因此該區(qū)的擠壓強度得以增加；（3）中國大陸淺部地應力場的狀態(tài)主要受到區(qū)域板塊構造環(huán)境、塊體邊界活動構造帶的展布和地形的控制，總體上以南北構造帶為界，西部以較強的壓性構造環(huán)境為主，東部為較弱的壓性構造環(huán)境，藏南和滇西北局部地區(qū)存在有張性構造環(huán)境；構造應力對地應力的貢獻比重隨著深度增加而增加；（4）采用黏彈性模型的構造應力場模擬結果比完全彈性模型的模擬結果能夠更好地與實測地應力場相吻合，利用完全彈性模型分析由地震等誘發(fā)的地應力瞬時變化是有效的；（5）青藏高原東南緣最大主應力方向發(fā)生了較大的偏轉，其主要控制因素有：印度板塊持續(xù)的碰撞、中下地殼對上地殼拖曳以及印度板塊通過實皆斷裂對歐亞板塊的剪切拉伸作用.中國大陸現(xiàn)今地應力場是整個地殼巖石黏彈特性長期演化和斷裂活動的結果，是地應力場動態(tài)演化過程中在現(xiàn)今時間點上的狀態(tài)，受到板塊構造環(huán)境、大陸內部活動斷裂分布、地形地貌和地殼結構等因素不同程度的控制，模擬結果為中國大陸地應力場提供了一個定量的參考模型.

中國大陸，地應力場，活動構造，黏彈性，有限元模擬

1 引 言

地應力是地震活動、造山事件、沉積盆地形成以及地下流體運動等地質過程的直接控制因素，也是采礦、水利、橋梁等工程建設中必須考慮的地質條件.中國大陸屬于歐亞板塊的一部分，被夾持在印度板塊、菲律賓海板塊、太平洋板塊和西伯利亞—蒙古塊體之間，受印度板塊與歐亞板塊的碰撞和菲律賓板塊的俯沖的影響，其內部構造運動極為活躍，板塊構造環(huán)境所控制的地應力場是活動構造的直接動力因素.通過活動斷裂、原地應力測量和地震活動的震源機制等研究能夠反映局部的地應力場狀態(tài)，并在一定程度上反映區(qū)域地應力場的分布趨勢特征［1］，盡管中國大陸地應力場在區(qū)域上受控于周緣的板塊構造環(huán)境，而在局部地區(qū)則顯示了活動斷裂構造和塊體構造單元地殼結構的影響，地應力場在不同的空間尺度具有不完全一致的控制因素.活動構造與地應力場的關系相互影響、密不可分，一方面，活動構造的活動性體現(xiàn)了局部應力場的特征；另一方面，活動構造的存在會造成局部地應力場的調整，從而影響地應力場在局部地區(qū)的分布狀態(tài)［2－5］.地形起伏和地殼介質密度分布不均勻引起的重力勢所產生的地應力場也是地應力分析中重要的因素之一，并引起許多學者的關注和研究［6－8］.中國大陸在不同區(qū)域地殼厚度分布不均勻，內部結構也存在很大的差異［9－16］，青藏高原的地殼厚達70多公里，而東部華北地區(qū)只有30km左右；地形地貌起伏很大，青藏高原平均高程達4km.地應力場在不同區(qū)域、不同深度其主要控制因素也有所區(qū)別.眾多學者［17－19］利用數值模擬方法從總體上分析了中國大陸構造應力場的特征和成因，模擬過程中主要有以下幾點考慮不夠：（1）模擬過程中都采用平面直角坐標系，而對于中國大陸整體上經緯度跨度都比較大，選用球殼模型更能夠接近實際；（2）在模擬過程中，邊界條件和模擬結果約束或者只選用GPS，或者只選用實際地應力場觀測數據，沒有把二者給予統(tǒng)一考慮；（3）重力作為影響地應力場分布的重要因素之一，缺乏對其定量的分析.完整地理解現(xiàn)今地應力場的狀態(tài)，不僅需要現(xiàn)今的構造環(huán)境和巖石圈展布特征，而且需要了解巖石圈在最近地質歷史時期內的運動狀態(tài)和動力學演化過程，特別是活動構造所反映的地應力場狀態(tài)和演化.本文以中國大陸主要活動斷裂所控制的塊體模型為基礎，采用黏彈性殼體有限單元，綜合考慮GPS、地應力實測結果以及重力等影響，分析中國大陸地應力場的特征和控制因素.

2 基本方法

從巖石力學上來說，盡管整個地殼的巖石主要表現(xiàn)為彈性，但是目前在中地殼普遍存在一個巖石強度相對弱的區(qū)域.地殼上部的構造塊體單元可以看作是完全彈性介質，活動斷裂作為塊體邊界具有與構造塊體本身完全不同的性質，活動斷裂運動除了伴隨地震發(fā)生瞬間明顯的斷層滑動外，更多的運動則表現(xiàn)為大量的無震滑動和長期蠕動等形式［20］，活動斷裂帶可以視為調整塊體彈性應力的軟弱帶.

現(xiàn)今中國大陸地應力場是在周緣板塊構造環(huán)境條件下，地應力場長期演化過程中的一個瞬間狀態(tài).地殼巖石在大時間尺度上通常可用標準線性體來近似，Maxwell體為標準線性體的一個特例（圖1）［21－24］，作為黏彈性的巖石具有應力松弛和應變蠕變的特性.其本構關系在Laplace域可表示為

其中

τ為弛豫時間：

圖1 Maxwell標準線性黏彈體（a）及其蠕變（b）和松弛曲線（c）Fig.1 （a）The Maxwell form of standard linear solid；（b）The creep compliance；（c）The relaxation modulus

中國大陸塊體通過板塊邊界受到相鄰板塊的作用力，這些作用力一方面通過完全彈性變形作用于中國大陸內部，另一方面，相鄰板塊持續(xù)的作用通過巖石的黏滯性對中國大陸內部作用，這種作用不同于完全彈性，在不同的時間不同塊體具有不同的狀態(tài).由于受到斷層庫侖摩擦破裂強度的限制，塊體內能夠維持的地應力通常符合Byerlee準則［25］.當板塊邊界作用于大陸地殼內部時，板塊邊界的作用力不會超過Byerlee準則所確定最大應力，所以邊界作用力在一定的地質歷史時間內，當板塊構造環(huán)境沒有發(fā)生變化時，應為一個持續(xù)恒定的應力，這種邊界作用力除了瞬時的彈性應力作用于塊體內部，還會由于不同塊體具有不同巖石黏滯性，導致在不同的時段具有不同的應力狀態(tài)，并最終達到平衡，大陸地殼整體都處于一個臨界的斷裂平衡的狀態(tài)，全球廣泛存在有水庫蓄水和流體注入引起的地震活動、一個地震的發(fā)生會觸發(fā)另一個地震的發(fā)生，原地深孔地應力測量結果符合Byerlee準則，這些現(xiàn)象恰好說明了目前大陸內部，特別是活動構造斷裂附近處于臨界斷裂平衡狀態(tài)［26－27］.

地應力場通常為構造應力場與重力勢引起的應力場的疊加，重力勢引起的應力場主要是由地形起伏和巖石層密度不均勻造成的，當地勢平緩且地下巖石分布均勻時，重力對地應力場影響可以忽略不計，構造應力場通常由于受到巖石強度的限制，在不同的構造環(huán)境具有不同的變化規(guī)律.

3 地質模型

3.1 中國大陸活動構造塊體模型

中國大陸是由不同活動塊體組成，塊體之間以大型的活動斷裂帶為分界，板內塊體具有不同活動程度的次級斷裂［28－29］.中國大陸及其鄰區(qū)的活動地塊主要有拉薩、羌塘、巴顏喀拉、柴達木、祁連、川滇、塔里木、天山、準噶爾、華北等（圖2）.隨著現(xiàn)代GPS觀測資料的不斷積累，初步構建了中國地殼活動的運動模型［30－33］，并在此基礎上分析了其應變場的特征［34－35］，中國大陸的現(xiàn)今構造變形主要以剛性地塊運動為主，如塔里木、鄂爾多斯、華南等地塊；只有少數塊體內部存在比較明顯的連續(xù)變形，如青藏高原和天山.塊體分界的大型活動斷裂區(qū)是變形的主要區(qū)域，邊界斷裂區(qū)的應變率遠大于塊體內部的應變率，通常邊界斷裂帶都具有數十公里的寬度，發(fā)育了中國大陸的主要地震活動［29，36］.周緣的印度板塊、菲律賓板塊和太平洋板塊的運動是中國大陸板塊內部構造活動的主要動力來源.印度板塊與歐亞板塊為陸陸碰撞，形成擠壓環(huán)境；菲律賓板塊和太平洋板塊與歐亞板塊為洋陸碰撞，太平洋洋殼和菲律賓的洋殼向歐亞板塊之下俯沖，形成了日本海和沖繩海槽，形成拉張的邊界環(huán)境.

3.2 中國大陸地殼結構

中國及鄰區(qū)地殼厚度分布的特點為從東到西逐漸增厚，在構造塊體內部地殼厚度通常變化較緩，如蒙古高原、華北、華中、華南、青藏高原內部等地區(qū)（圖2）.西部青藏高原的地殼厚度平均為60～70km，局部達到74km，而東部地區(qū)的厚度平均為30～36km，沿海地區(qū)僅有28～30km；地殼局部厚度變化與構造環(huán)境密切相關，大型地塊的邊界帶處的地殼厚度梯度也比較大，南北構造帶以東地殼厚度為40km，而在南北構造帶以西達到60km，在很短的距離上地殼厚度發(fā)生了很大的變化.中國大陸的地殼在沉降帶或盆地發(fā)育地方的地殼較周圍薄一些，如塔里木、準格爾、柴達木、鄂爾多斯等盆地，而抬升區(qū)或造山帶處的地殼較周圍厚一些，如青藏高原、天山、太行山等地區(qū)［37－39］.穩(wěn)定的克拉通塔里木、華南、鄂爾多斯地塊處，其地殼內低速層不發(fā)育，而其它地區(qū)都普遍有地殼低速層發(fā)育，特別在青藏高原的中地殼低速層整體上普遍存在［9］.而細致地殼的流變結構成果［40］為本文進行地應力場的數值分析提供了直接的參數依據.

4 有限元模型的建立

4.1 有限元網格模型

在模擬過程中，以中國大陸塊體模型為基礎，采用三維球殼黏彈線性體單元，對塊體模型進行無結構網格剖分（圖3），模型的地形分布參照GTOPO30數字高程模型，按照區(qū)域板塊構造環(huán)境，設定邊界條件和加載.塊體邊界的大規(guī)模活動斷裂帶的網格單元劃分相對細密，從而保證邊界斷裂帶的數值計算的穩(wěn)定性.各個塊體和斷裂帶分層巖石物性參數見表1，各個塊體分為上地殼、中地殼和下地殼三層，表中的楊氏模量和泊松比為瞬態(tài)泊松比，剪切模量系數是指巖石剪切模量在經過足夠長的時間后，其長期剪切松弛模量與瞬時剪切模量的近似比值，長時間的泊松比近似為0.5，巖石密度上地殼取為2750kg／m3，中下地殼為3300kg／m3.

圖2 中國大陸構造塊體分區(qū)及板塊構造環(huán)境Fig.2 China mainland tectonic blocks and its surrounding plate tectonics

圖3 有限元網格模型及邊界條件Fig.3 FEM mesh model and boundary condition

4.2 邊界條件和加載

模型北邊界中段為南北向和東西向均為約束邊界；北邊界東西兩端和南邊界東段均為南北向約束，東西向自由；西邊界北段和東邊界南端均為東西向約束，南北向自由；其他邊界均按照構造應力環(huán)境加載應力.

根據Anderson斷層理論可以確定不同性質的斷層所在區(qū)域的構造應力特征［27，41］.在正斷層、走滑斷層和逆沖斷層活動的不同區(qū)域，其應力特征分別為：

其中σ1為最大主應力，σ3為最小主應力，SHmax為水平最大主應力，Shmin為水平最小主應力，SV為垂向應力，Pf為孔隙流體壓力，μ為斷層臨界摩擦系數，根據Byerlee實驗數據［25］和深部地應力測量結果［27］，其取值范圍在0.6～1.0之間.由此可得出正斷層區(qū)域，垂直斷層的方向為最小主應力方向，也是水平最小主應力，其大小約為

表1 模型巖石物性參數Table 1 Physical property parameters of rocks in model

在逆沖斷層區(qū)域，垂直斷層的方向為最大主應力方向，也是水平最大主應力，其大小約為

通常巖石流體密度比ρrock／ρliquid約為3.0，

當μ＝0.6時，Shmin≥0.55ρrockgh，SHmax≤2.41ρrockgh，當μ＝1.0時，Shmin≥0.45ρrockgh，SHmax≤4.22ρrockgh.

由于板塊邊界的最大主應力方向與板塊邊界近似垂直，所以按照板塊構造環(huán)境對邊界加載最大主應力或最小主應力的應力載荷，除去重力作用的影響，邊界條件見表2.

表2 邊界加載應力表Table 2 The load on the boundary

5 模擬結果分析

圖4為在邊界加載壓力約10萬年后的近地表的速度分布，在這一時間點上，模擬速度整體上能夠很好地與現(xiàn)今實測的GPS速度場吻合，因此，現(xiàn)今構造應力場的分布特征可參考該時點的應力狀態(tài)模擬結果.圖5為在只加載邊界應力，而不考慮重力作用，約10萬年后的距地表1000m深度的構造應力場模擬結果，總體上構造應力場在西部為擠壓，東部為拉張，構造應力場最大主應力的方向受到板塊構造環(huán)境和活動構造分布的控制，印度板塊與歐亞板塊的持續(xù)碰撞形成了青藏高原及其周緣的擠壓性質的構造應力場，而東部菲律賓板塊與太平洋板塊的俯沖形成了黃海、東海和環(huán)渤海區(qū)域的拉張性質構造應力場，中間為拉張環(huán)境和擠壓環(huán)境的過渡區(qū)域.

圖4 模擬速度與現(xiàn)今實測GPS速度對比圖圖中白色的箭頭為實測的GPS速度場，參照系為穩(wěn)定歐亞大陸板塊（據王琪等［42］和中國地殼運動觀測網絡），黑色箭頭為數值模擬結果.Fig.4 Comparison between the modeling velocity and the measured velocity with GPSThe white arrows are GPS measured velocity，the reference frame is stable Europe－Asia continental plate（reference Wang［42］and the China Crust Movement Observation Network），the black arrows are modeling results.

圖5 中國大陸構造應力場水平最大主應力模擬結果分布圖（黏彈性模擬，1000m深度）其中紅色箭頭為水平最大主應力的方向，云圖為水平最大主應力值，h為深度，紅色為張應力，藍色為壓應力Fig.5 Modeling tectonic horizontal maximum principal stress of China mainland with viscoelastic model at 1000 meters depthRed line arrows are the direction of maximum principal stress，cloud image is the value of maximum principal stress，h is depth，red zone is extensive，blue zone is compressive.

圖6 中國大陸重力引起的水平最大主應力模擬結果分布圖（1000m深度）其中紅色箭頭為水平最大主應力的方向，云圖為水平最大主應力值，h為深度，紅色為張應力，藍色為壓應力.Fig.6 Modeling maximum principal stress of China mainland induced by gravity（at depth 1000m）Red line arrows are the direction of maximum principal stress，cloud image is the value of maximum principal stress，h is depth，red zone is extensive，blue zone is compressive.

圖6為中國大陸在1000m深度重力勢引起的偏應力模擬結果，由于中國大陸地勢西高東低，整體上形成三級階梯，第一級最高的階梯為青藏高原，平均海拔在4000m以上；第二級階梯位于青藏高原與大興安嶺─太行山─巫山─雪峰山之間，包括內蒙古高原、黃土高原、云貴高原和塔里木盆地、準噶爾盆地、四川盆地等地區(qū)，海拔一般為1000～2000m，惟四川盆地較低，海拔在500m以下；第三級最低階梯，在大興安嶺─太行山─巫山─雪峰山以東，自北而南，有海拔200m以下的東北平原、華北平原和長江中下游平原.重力勢引起的應力場與構造應力場的結果大致相反，在青藏高原的藏東南、四川西部和云南西北地區(qū)為拉張性質的應力場，而在塔里木盆地和東部的平原區(qū)為擠壓性質的應力場，應力場的方向受地勢展布特征控制.圖7為1000m深度位置的中國大陸地應力場，是由構造應力的模擬結果與重力作用模擬結果進行張量求和合成，可以作為淺部地應力場的定量參考模型，地殼深部的地應力場應與構造應力場相近.模擬地應力場與實測的地應力分布特征［43］（圖8）也比較一致，水平最大主應力的方向主要受到板塊構造環(huán)境和活動斷裂帶的控制，總體上以南北構造帶為界，以西以較強的壓性構造環(huán)境為主，以東為較弱的壓性地應力環(huán)境，由于受到重力的影響，藏南和滇西北局部地區(qū)的淺部存在有張性地應力環(huán)境，而在塔里木區(qū)域，由于重力勢引起的應力場與構造應力場都為擠壓性質，因此擠壓強度得以增加.圖9為初期瞬時的模擬結果，可作為完全彈性結果，完全彈性模型的構造應力場模擬結果與實測的地應力場存在一定的偏差，如在龍門山斷裂帶的北段地區(qū)，水平最大主應力的方向總體為北西—南東向，而完全彈性模型的結果為北東—南西向.青藏高原南緣在重力作用下形成了張性應力環(huán)境，但是由于重力作用影響有限，在深部的張性環(huán)境應與深部物質上涌等局部構造活動有更為密切的關系.

圖7 中國大陸地應力場的水平主應力模擬結果分布圖（1000m深度）（a）水平最大主應力，其中紅色箭頭為水平最大主應力的方向，云圖為水平最大主應力值，h為深度，紅色為張應力，藍色為壓應力；（b）水平最小主應力，其中紅色箭頭為水平最小主應力的方向，云圖為水平最小主應力值，h為深度，紅色為張應力，藍色為壓應力；（c）為剪應力，其中紅色箭頭為水平最大主應力的方向，云圖為剪應力值，h為深度.Fig.7 Modeling crustal horizontal stress in China mainland（at depth 1000m）（a）Horizontal maximum principal stress，Red line arrows are the direction of maximum principal stress，cloud image is the value of maximum principal stress，his depth，red zone is extensive，blue zone is compressive；（b）Horizontal minim principal stress，Red line arrows are the direction of minim principal stress，cloud image is the value of minim principal stress，his depth，red zone is extensive，blue zone is compressive；（c）Shear stress，Red line arrows are the direction of maximum principal stress，cloud image is the value of shear stress，his depth

圖8 中國現(xiàn)代構造應力場圖（謝富仁等［43］，2003）Fig.8 The measured crust stress of China（reference the Crust Stress Databese of Institute of Crustal Dynamics［43］）

中國大陸地應力場和現(xiàn)今活動構造的性質總體上主要受到印度洋板塊的碰撞擠壓和太平洋板塊與菲律賓板塊俯沖的影響.由于受到華南板塊、塔里木克拉通以及鄂爾多斯塊體的阻擋，使得印度洋板塊碰撞作用集中在青藏高原及其周邊，青藏高原主體多為活動斷裂的走滑，塊體內部的連續(xù)變形較主要斷裂活動要小很多.華北塊體、山西地塹除了受到印度洋板塊碰撞的影響之外，菲律賓和太平洋俯沖板塊所形成的弧后張性環(huán)境的影響不容忽視.

圖9 中國大陸構造應力場水平最大主應力模擬結果分布圖（完全彈性模擬）其中紅色箭頭為水平最大主應力的方向，云圖為水平最大主應力值，h為深度，紅色為張應力，藍色為壓應力.Fig.9 Modeling tectonic horizontal maximum principal stress in China mainland with pure elastic modelRed line arrows are the direction of maximum principal stress，cloud image is the value of maximum principal stress，h is depth，red zone is extensive，blue zone is compressive.

巖石蠕變曲線表明，在給定的應力條件下，巖石蠕變逐漸趨于一個常數，這就意味著隨時間的演化，巖石塊體整體可能會運動，但是由于其應變率逐漸趨于0，所以巖石塊體內部處于一個穩(wěn)定的無變形狀態(tài).一方面，上地殼的變形受控于自身的構造塊體結構，在板塊構造環(huán)境控制下，如果不考慮其底部中下地殼的作用，在僅受板塊構造環(huán)境的作用，其應力狀態(tài)和變形特征應接近于完全彈性模型的數值分析結果；另一方面，當柔性較強的中下地殼受到與上地殼同樣的板塊構造環(huán)境的作用時，中下地殼會表現(xiàn)出更強的運動性，從而對上地殼的變形和應力分布產生影響.例如在青藏高原，由于中下地殼向兩側擠出，導致上地殼被中下地殼拖曳，使得上地殼在活動斷裂的控制下隨中下地殼運動，現(xiàn)今中下地殼對上地殼的變形模式和應力分布特征的影響作用不容忽視.

青藏高原東南部主應力方向出現(xiàn)轉折，主要有以下幾個因素控制：（1）由于印度板塊持續(xù)的碰撞，使得青藏高原整體受到塔里木克拉通、華南板塊和鄂爾多斯塊體的阻擋，使得青藏高原隨著活動斷裂的滑動向兩側擠出［44］，在擠出過程中，剛性的上地殼塊體構造的運動受邊界斷裂的控制；（2）由于中下地殼較上地殼有更強的柔性，所以擠出時上地殼的運動受到中下地殼拖曳的影響；（3）印度板塊向北運動過程中，通過緬甸的實皆斷裂（Sagaing Fault）剪切拉伸作用于歐亞板塊，使得上地殼中青藏高原東南部應力和地表運動方向都發(fā)生了轉向.

6 結論與討論

本文基于線性黏彈體模型對現(xiàn)今中國大陸的構造應力場進行了數值模擬分析，其前提假設是中國大陸周邊的板塊構造環(huán)境相對穩(wěn)定，且維持了足夠長的時間，從而能夠保證模擬分析過程中邊界條件的一致性；大型活動斷裂所決定的塊體結構保持穩(wěn)定，從而保證在模擬分析過程中地質模型的物性參數保持一致.該結果主要考慮了板塊構造作用和區(qū)域性的地貌和地殼結構所產生的重力的影響，在進一步的研究中需要考慮底部地幔上涌和巖石圈熱演化等影響因素，深入討論地應力場控制因素.

整體上，中國大陸構造應力場的分布特征與活動構造的展布密切相關，是在板塊構造環(huán)境控制下，整個地殼巖石黏彈特性長期演化和活動斷裂無震蠕動的結果.而伴隨地震發(fā)生的斷層滑動引起的地應力場的瞬時或短期調整，可參考完全彈性分析的結果.陳群策等（私人通信）最近在廣元附近進行的汶川地震后的地應力測量，所測的地應力水平最大主應力方向為北東—南西向，而該地區(qū)通常的區(qū)域應力場的分布特征為北西—南東向.應力場作為動態(tài)演化的過程，特別是斷層附近的地應力場，由于受到斷層蠕動和滑動的影響，地應力的變化會更為突出.郭啟良［45］和廖椿庭［46］分別在汶川地震前后和昆侖山地震前后的地應力測量所反映的地應力的變化，也同樣說明，活動斷裂附近的地應力場受到局部構造的影響，從而體現(xiàn)了不穩(wěn)定性.隨著地應力監(jiān)測網絡的建立，可以有效地利用地應力動態(tài)的信息，更好地認識區(qū)域構造應力與局部地質構造和巖石物性之間的關系.

致 謝 感謝匿名評審人提出的寶貴意見！

（References）

［1］ 謝富仁，崔效鋒，趙建濤等.中國大陸及鄰區(qū)現(xiàn)代構造應力場分區(qū).地球物理學報，2004，47（4）：654－662.Xie F R，Cui X F，Zhao J T，et al.Regional division of the recent tectonic stress field in China and adjacent areas.Chinese J.Geophys.（in Chinese），2004，47（4）：654－662.

［2］ 劉峽，馬瑾，傅容珊等.華北地區(qū)現(xiàn)今地殼運動動力學初步研究.地球物理學報，2010，53（6）：1418－1427.Liu X，Ma J，F(xiàn)u R S，et al.Primary study on the dynamics of the present－day crustal motion in North China Region.Chinese J.Geophys.（in Chinese），2010，53（6）：1418－1427.

［3］ 朱守彪，張培震，石耀霖.華北盆地強震孕育的動力學機制研究.地球物理學報，2010，53（6）：1409－1417.Zhu S B，Zhang P Z，Shi Y L.A study on the mechanisms of strong earthquake occurrence in the North China Basin.Chinese J.Geophys.（in Chinese），2010，53（6）：1409－1417.

［4］ 石耀霖，曹建玲.庫侖應力計算及應用過程中若干問題的討論——以汶川地震為例.地球物理學報，2010，53（1）：102－110.Shi Y L，Cao J L.Some aspects in static stress change calculation——case study on Wenchuan earthquake.Chinese J.Geophys.（in Chinese），2010，53（1）：102－110.

［5］ 酆少英，高銳，龍長興等.銀川地塹地殼擠壓應力場：深地震反射剖面.地球物理學報，2011，54（3）：692－697，doi：10.3969／j.issn.0001－5733.2011.03.008.Feng S Y，Gao R，Long C X，et al.The compressive stress field of Yinchuan garben：Deep seismic reflection profile.Chinese J.Geophys.（in Chinese），2011，54（3）：692－697，doi：10.3969／j.issn.0001－5733.2011.03.008.

［6］ Zoback M L.First－and second－order patterns of stress in the lithosphere：The World Stress Map project.J.Geophys.Res.，1992，97（8）：11703－11728.

［7］ Flesch L M，Haines A J，Holt W E.Dynamics of the India－Eurasia collision zone.J.Geophys.Res.，2001，106（B8）：16435－16460.

［8］ 朱守彪，石耀霖.青藏高原地形擴展力以及下地殼對上地殼的拖曳力的遺傳有限單元法反演.北京大學學報：自然科學版，2005，41（2）：225－234.Zhu S B，Shi Y L.Genetic algorithm finite element inversion of topographic spreading forces and drag forces of lower crust to upper crust in Tibetan Plateau.Acta Scientiarum Naturalum Universitis Pekinesis（in Chinese），2005，41（2）：225－234.

［9］ 金勝，魏文博，汪碩等.青藏高原地殼高導層的成因及動力學意義探討——大地電磁探測提供的證據.地球物理學報，2010，53（10）：2376－2385，doi：10.3969／j.issn.0001－5733.2010.10.011.Jin S，Wei W B，Wang S，et al.Discussion of the formation and dynamic signification of the high conductive layer in Tibetan crust.Chinese J.Geophys.（in Chinese），2010，53（10）：2376－2385，doi：10.3969／j.issn.0001－5733.2010.10.011.

［10］ 鄧陽凡，李守林，范蔚茗等.深地震測深揭示的華南地區(qū)地殼結構及其動力學意義.地球物理學報，2011，54（10）：2560－2574，doi：10.3969／j.issn.0001－5733.2011.10.013.Deng Y F，Li S L，F(xiàn)an W M，et al.Crustal structure beneath South China revealed by deep seismic soundings and its dynamics implications.Chinese J.Geophys.（in Chinese），2011，54（10）：2560－2574，doi：10.3969／j.issn.0001－5733.2011.10.013.

［11］ 葛粲，鄭勇，熊熊.華北地區(qū)地殼厚度與泊松比研究.地球物理學報，2011，54（10）：2538－2548，doi：10.3969／j.issn.0001－5733.2011.10.011.Ge C，Zheng Y，Xiong X.Study of crustal thickness and Poisson ratio of the North China Craton.Chinese J.Geophys.（in Chinese），2011，54（10）：2538－2548，doi：10.3969／j.issn.0001－5733.2011.10.011.

［12］ 滕吉文，張洪雙，孫若昧等.青藏高原腹地東西分區(qū)和界帶的地球物理場特征與動力學響應.地球物理學報，2011，54（10）：2510－2527，doi：10.3969／j.issn.0001－5733.2011.10.009.Teng J W，Zhang H S，Sun R M，et al.Geophysical field characteristics and dynamic response of segmentations in eastwest direction and their boundary zone in central Tibetan plateau.Chinese J.Geophys.（in Chinese），2011，54（10）：2510－2527，doi：10.3969／j.issn.0001－5733.2011.10.009.

［13］ 劉瓊林，王椿鏞，姚志祥等.華北克拉通中西部地區(qū)地殼厚度與波速比研究.地球物理學報，2011，54（9）：2213－2224，doi：10.3969／j.issn.0001－5733.2011.09.003.Liu Q L，Wang C Y，Yao Z X，et al.Study on crustal thickness and velocity ratio in mid－western North China Craton.Chinese J.Geophys.（in Chinese），2011，54（9）：2213－2224，doi：10.3969／j.issn.0001－5733.2011.09.003.

［14］ 陳石，王謙身，祝意青等.青藏高原東緣重力導納模型均衡異常時空特征.地球物理學報，2011，54（1）：22－34，doi：10.3969／j.issn.0001－5733.2011.01.004.Chen S，Wang Q S，Zhu Y Q，et al.Temporal and spatial features of isostasy anomaly using gravitational admittance model at eastern margin of Tibetan Plateau.Chinese J.Geophys.（in Chinese），2011，54（1）：22－34，doi：10.3969／j.issn.0001－5733.2011.01.004.

［15］ 李飛，周仕勇，蘇有錦等.川滇及鄰區(qū)Pn波速度結構和各向異性研究.地球物理學報，2011，54（1）：44－54，doi：10.3969／j.issn.0001－5733.2011.01.006.Li F，Zhou S Y，Su Y J，et al.Study on Pn－wave velocity structure and anisotropy in the Sichuan－Yunnan and its adjacent areas.Chinese J.Geophys.（in Chinese），2011，54（1）：44－54，doi：10.3969／j.issn.0001－5733.2011.01.006.

［16］ 于湘?zhèn)ィ愡\泰，張懷.京津唐地區(qū)地殼三維P波速度結構與地震活動性分析.地球物理學報，2010，53（8）：1817－1828.Yu X W，Chen Y T，Zhang H.Three－dimensional crustal P－wave velocity structure and seismicity analysis in Beijing－Tianjin－Tangshan Region.Chinese Journal Geophysics（in Chinese），2010，53（8）：1817－1828.

［17］ 汪素云，陳培善.中國及鄰區(qū)現(xiàn)代構造應力場的數值模擬.地球物理學報，1980，23（1）：35－45.Wang S Y，Chen P S.A numerical simulation of the present tectonic stress field in China and its vicinity.Chinese J.Geophys.（in Chinese），1980，23（1）：35－45.

［18］ 朱守彪，石耀霖.中國大陸及鄰區(qū)構造應力場成因的研究.中國科學D輯，2006，36（12）：1077－1083.Zhu S B，Shi Y L.Study the controls of the tectonic stress of China mainland.Science in China Series D：Earth Sciences（in Chinese），2006，36（12）：1077－1083.

［19］ Liu M，Yang Y Q，Shen Z K，et al.Active tectonics and intracontinental earthquakes in China：The kinematics and geodynamics.∥Stein S，Mazzotti S eds.Continental Intraplate Earthquakes：Science Hazard and Policy Issues.Geological Society of America Special Paper.2007：299－318，doi：10.1130／2007.2425（19）.

［20］ 石耀霖.關于應力觸發(fā)和應力影概念在地震預報中應用的一些思考.地震，2001，21（3）：1－7.Shi Y L.Stress triggers and stress shadows：How to apply these concepts to earthquake prediction.Earthquake（in Chinese），2001，21（3）：1－7.

［21］ 程佳，劉杰，甘衛(wèi)軍等.1997年以來巴顏喀拉塊體周緣強震之間的黏彈性觸發(fā)研究.地球物理學報，2011，54（8）：1997—2010，doi：10.3969／j.issn.0001－5733.2011.08.007.Cheng J，Liu J，Gan W J，et al.Coulomb stress interaction among strong earthquakes around the Bayan Har block since the Manyi earthquake in 1997.Chinese J.Geophys.（in Chinese），2011，54（8）：1997—2010，doi：10.3969／j.issn.0001－5733.2011.08.007.

［22］ 張晁軍，曹建玲，石耀霖.從震后形變探討青藏高原下地殼黏滯系數.中國科學（D輯），2008，38（10）：1250－1257.Zhang C J，Cao J L，Shi Y L.Study the viscosity of crust using the post－earthquake deformation.Science in China：Earth Sciences（Series D）（in Chinese），2008，38（10）：1250－1257.

［23］ 白玉柱，徐錫偉，徐杰等.2008年汶川地震近斷裂區(qū)域變形場的空間分布.地球物理學報，2011，54（7）：1798－1804，doi：10.3969／j.issn.0001－5733.2011.07.013.Bai Y Z，Xu X W，Xu J，et al.A research on the distribution of deformation fields near the fault of 2008Wenchuan Earthquake.Chinese J.Geophys.（in Chinese），2011，54（7）：1798－1804，doi：10.3969／j.issn.0001－5733.2011.07.013.

［24］ Turcotte D L，Schubert G.Geodynamics（2nd Edition）.New York：Cambridge University Press，2002.

［25］ Byerlee J D.Friction of rocks.Pure and Applied Geophysics，1978，116：615－629.

［26］ Townend J，Zoback M D.How faulting keeps the crust strong.Geology，2000，28（5）：399－402.

［27］ Zoback M D，Townend J.Implications of hydrostatic pore pressures and high crustal strength for the deformation of intraplate lithosphere.Tectonophysics，2001，336（1－4）：19－30.

［28］ 鄧起東，張培震，冉勇康等.中國活動構造基本特征.中國科學（D輯），2002，32（12）：1020－1030.Deng Q D，Zhang P Z，Ran Y K，et al.The basic character of China active tectonics.Science in China：Earth Sciences（Series D）（in Chinese），2002，32（12）：1020－1030.

［29］ 張培震，鄧起東，張國民等.中國大陸的強震活動與活動地塊.中國科學（D輯），2003，33（增刊）：12－20.Zhang P Z，Deng Q D，Zhang G M，et al.Active tectonic blocks and strong earthquakes in the continent of China.ScienceinChina（Series D）（in Chinese），2003，46（Suppl）：13－24.

［30］ 唐方頭，宋鍵，曹忠權等.最新GPS數據揭示的東構造結周邊主要斷裂帶的運動特征.地球物理學報，2010，53（9）：2119－2128.Tang F T，Song J，Cao Z Q，el al.The movement characters of main faults around Eastern Himalayan Syntaxis revealed by the latest GPS data.Chinese J.Geophys.（in Chinese），2010，53（9）：2119－2128.

［31］ Wang Q，Zhang P Z，F(xiàn)reymueller J T，et al.Present－day crustal deformation in China constrained by global positioning system measurements.Science，2001，294（5542）：574－577.

［32］ 張培震，甘衛(wèi)軍，沈正康等.中國大陸現(xiàn)今構造作用的地塊運動和連續(xù)變形耦合模型.地質學報，2005，79（6）：748－756.Zhang P Z，Gan W J，Shen Z K，et al.A coupling model of rigid－block movement and continuous deformation：patterns of the present－day deformation of China′s continent and its vicinity.Acta Geologica Sinica（in Chinese），2005，79（6）：748－756.

［33］ 石耀霖，朱守彪.利用GPS觀測資料劃分現(xiàn)今地殼活動塊體的方法.大地測量與地球動力學，2004，24（2）：1－5.Shi Y L，Zhu S B.Method for division of present active crustal blocks by GPS survey data.Journal of Geodesy and Geodynamics（in Chinese），2004，24（2）：1－5.

［34］ 石耀霖，朱守彪.用GPS位移資料計算應變方法的討論.大地測量與地球動力學，2006，26（1）：1－8.Shi Y L，Zhu S B.Discussion on method of calculating strain with GPS displacement data.Journal of Geodesy and Geodynamics（in Chinese），2006，26（1）：1－8.

［35］ Chen Q Z，F(xiàn)reymueller J T，Wang Q，et al.A deforming block model for the present－day tectonics of Tibet.J.Geophys.Res.，2004，109：B01403，doi：10.1029／2002JB002151.

［36］ 張國民，馬宏生，王輝等.中國大陸活動地塊邊界帶與強震活動.地球物理學報，2005，48（3）：602－610.Zhang G M，Ma H S，Wang H，et al.Boundaries between active－tectonic blocks and strong earthquakes in the China mainland.Chinese J.Geophys.（in Chinese），2005，48（3）：602－610.

［37］ Li S L，Mooney W D，F(xiàn)an J C.Crustal structure of mainland China from deep seismic sounding data.Tectonophysics，2006，420（1－2）：239－252.

［38］ 蔡學林，朱介壽，曹家敏等.中國大陸及鄰區(qū)巖石圈地殼三維結構與動力學型式.中國地質，2007，34（4）：543－557.Cai X L，Zhu J S，Cao J M，et al.3Dstructure and dynamic types of the lithospheric crust in continental China and its adjacent regions.Geology in China（in Chinese），2007，34（4）：543－557.

［39］ 馬杏垣.中國巖石圈動力學地圖集.北京：中國地圖出版社，1989：1－68.Ma X Y.Atlas of the Lithosphere Dynamics of China（in Chinese）.Beijing：Cartographic Press，1989：1－68.

［40］ 石耀霖，曹建玲.中國大陸巖石圈等效粘滯系數的計算和討論.地學前緣，2008，15（3）：82－95.Shi Y L，Cao J L.Effective viscosity of China continental lithosphere.Earth Science Frontiers（in Chinese），2008，15（3）：82－95.

［41］ Zoback M D.Reservoir Geomechanics.New York：Cambridge University Press，2007：123－139.

［42］ 王琪，張培震，馬宗晉.中國大陸現(xiàn)今構造變形GPS觀測數據與速度場.地學前緣，2002，9（2）：415－429.Wang Q，Zhang P Z，Ma Z J.GPS database and velocity field of contemporary tectonic deformation in continental China.Earth Science Frontiers（in Chinese），2002，9（2）：415－429.

［43］ 謝富仁，陳群策，崔效鋒等.中國大陸地殼應力環(huán)境研究.北京：地質出版社，2003.Xie F R，Chen Q C，Cui X F，et al.Crustal Stress in China（in Chinese）.Beijing：Geological Press，2003.

［44］ Cook K L，Royden L H.The role of crustal strength variations in shaping orogenic plateaus，with application to Tibet.J.Geophys.Res.，2008，113，B08407，doi：10.1029／2007JB005457.

［45］ 郭啟良，王成虎，馬洪生等.汶川Ms8.0級大震前后的水壓致裂原地應力測量.地球物理學報，2009，52（5）：1395－1401，doi：10.3969／j.issn.0001－5733.2009.05.029.Guo Q L，Wang C H，Ma H S，et al.In－site hydro－fracture stress measurement before and after the Wenchuan Ms8.0 earthquake of China.Chinese J.Geophys.（in Chinese），2009，52（5）：1395－1401，doi：10.3969／j.issn.0001－5733.2009.05.029.

［46］ Liao C T，Zhang C S，Wu M L.Stress change near the Kunlun fault before and after the Ms8.1Kunlun earthquake.Geophysical Research Letters，2003，30（20）：2027－2030，doi：10.1029／2003GL018106.

Comprehensive modeling on the present crustal stress of China mainland with the viscoelastic spherical shell

FAN Tao－Yuan，LONG Chang－Xing，YANG Zhen－Yu，CHEN Qun－Ce，WU Zhong－Hai，SHAO Zhao－Gang，TONG Ya－Bo
Key Laboratory of Geotectonic Movement ＆Geohazard，Institute of Geomechanics，Chinese Academy of Geological Sciences，Beijing100081，China

The observed crustal stress of China mainland is related to the surrounding platetectonics，active tectonics，topography and the crustal structure.On the west of China mainland the Eurasia plate collides with the India plate，on the east the Philippine plate and Pacific plate subducts under the Eurasia plate.Large active faults divide China mainland into many blocks.The topography，crustal structure and thickness vary with the different block.Taking GPS and measured crust stress as restriction，we model the crustal stress using the viscoelastic spherical shell finite element（FE）method including the gravity factor.The finite element mesh was adapted to the active tectonics model，the topography and crust structure.The results are as follows.（1）The tectonic stress generated by the plate movement varies with different tectonic area and depth.Qingzang area is compressive，the Yellow Sea，East China Sea and the Bohai Sea are extensive，the middle area including North China and South China is transition zone.The direction of maximum principal stress is controlled by the plate movement and the active faults.（2）In the area of large topography gradient and crustal thickness gradient the gravity plays an important role in the distribution of shallow crustal stress.The shallow crustal stress in southern Tibet is extensive and is induced by the gravity，in the Tarim area the gravity influence enforces the tectonic compressive stress so the stress is strengthened.（3）The surround plate movement，interior active faults and topography determine the characters of shallow crustal stress of China mainland.The north－south seismic zone is the division of crustal stress state，west of which is strongly compressive zone，east of which is weakly compressive zone，some areas of southern Tibet and northwestern Yunnan are extensive zone.The tectonic stress contribution to crustal stress becomes more and more important as the depth increases.（4）The result of viscoelastic FE modeling fits better to the measured crustal stress than the elastic FE modeling.The result of elastic FE modeling is helpful to analyze the instantaneous stress change generated by the earthquake.（5）The crustal stress of southeast Qingzang plateau varies greatly with space.The crustal movement and the direction of maximum principal stress rotate around the Eastern Himalayan Syntaxis.The predominant controlling factors include plate collision of India plate and Eurasia plate，the India plate shear dragging through Sagaing fault and the lower crust dragging.The present crustal stress of China mainland is a time slice of viscoelatic stress evolution of crustal blocks and active tectonics.The crustal movement and deformation are related to the viscoelasticity of crust rock and the active tectonic aseismic creep.The active tectonic deformation is larger than the blocks.The surround plate tectonics，inner active tectonics，topography and crustal structure contribute to the crustal stress with different ratio in different area and different depth.Near the active fault the direction of maximum principal stress changes due to the weak fault zone.We provide a quantitative model of crustal stress of the China mainland which can be used for the crustal stress extrapolation.

China mainland，Crustal stress，Active tectonics，Viscoelastic，F(xiàn)E modeling

10.6038／j.issn.0001－5733.2012.04.020

P313

2011－08－17，2012－02－20收修定稿

中央級科研院所基本科研業(yè)務專項資助項目（DZLXJK201101、DZLXJK200801），中國地質調查項目（1212010911049，1212011120163）和國土資源部公益性行業(yè)科研專項（SinoProbe－06）共同資助.

范桃園，男，1971年生，2001年畢業(yè)于中國科學院研究生院獲博士學位，現(xiàn)為中國地質科學院地質力學研究所副研究員，主要從事地球動力學、地殼應力場方面的研究.E－mail：fan3ty＠gmail.com

范桃園，龍長興，楊振宇等.中國大陸現(xiàn)今地應力場黏彈性球殼數值模擬綜合研究.地球物理學報，2012，55（4）：1249－1260，

10.6038／j.issn.0001－5733.2012.04.020.

Fan T Y，Long C X，Yang Z Y，et al.Comprehensive modeling on the present crustal stress of China mainland with the viscoelastic spherical shell.Chinese J.Geophys.（in Chinese），2012，55（4）：1249－1260，doi：10.6038／j.issn.0001－5733.2012.04.020.

（本文編輯 胡素芳）