中国大陆现今地应力场黏弹性球壳数值模拟综合研究

范桃园，龙长兴，杨振宇，陈群策，吴中海，邵兆刚，仝亚博

国土资源部新构造运动与地质灾害重点实验室，中国地质科学院地质力学研究所，北京 100081

中国大陆现今地应力场黏弹性球壳数值模拟综合研究

范桃园，龙长兴，杨振宇，陈群策，吴中海，邵兆刚，仝亚博

国土资源部新构造运动与地质灾害重点实验室，中国地质科学院地质力学研究所，北京 100081

中国大陆现今实测地应力场的状态与板块构造环境、活动断裂带分布、地形地貌以及地壳结构呈现一定相关性.在中国大陆西缘，印度洋板块与欧亚板块陆发生陆碰撞，在中国大陆东缘，菲律宾海板块、太平洋板块俯冲到欧亚板块之下.中国大陆内部被大型活动断裂带分割为多个块体，各个块体的地壳结构和厚度呈不均匀分布，地形地貌起伏具有很大的差异.笔者以中国大陆块体模型为基础，把板块构造作用和重力势作为主要影响地应力状态的两个主要要素，在现今活动构造、GPS和实测地应力等成果的约束下，利用线性黏弹体球壳有限元模拟分析了中国大陆现今地应力场的分布特征和控制因素.结果表明：（1）构造应力场总体上呈现出西部挤压，东部拉张的特征，印度板块与欧亚板块的持续碰撞形成了青藏高原及其周缘的挤压性质的构造应力场，而东部菲律宾板块与太平洋板块的俯冲形成了黄海、东海和环渤海区域的拉张性质的构造应力场，中间为拉张环境和挤压环境的过渡，最大主应力的方向受到板块构造环境和活动构造分布的控制；（2）重力的影响主要体现在地形梯度大和地壳厚度结构变化大的地壳浅部区域，在藏南、滇西北局部地区的地壳浅部由于受到重力势控制，呈现为张性应力场，在塔里木地区由于重力势引起的应力场与构造应力场同为挤压性质，因此该区的挤压强度得以增加；（3）中国大陆浅部地应力场的状态主要受到区域板块构造环境、块体边界活动构造带的展布和地形的控制，总体上以南北构造带为界，西部以较强的压性构造环境为主，东部为较弱的压性构造环境，藏南和滇西北局部地区存在有张性构造环境；构造应力对地应力的贡献比重随着深度增加而增加；（4）采用黏弹性模型的构造应力场模拟结果比完全弹性模型的模拟结果能够更好地与实测地应力场相吻合，利用完全弹性模型分析由地震等诱发的地应力瞬时变化是有效的；（5）青藏高原东南缘最大主应力方向发生了较大的偏转，其主要控制因素有：印度板块持续的碰撞、中下地壳对上地壳拖曳以及印度板块通过实皆断裂对欧亚板块的剪切拉伸作用.中国大陆现今地应力场是整个地壳岩石黏弹特性长期演化和断裂活动的结果，是地应力场动态演化过程中在现今时间点上的状态，受到板块构造环境、大陆内部活动断裂分布、地形地貌和地壳结构等因素不同程度的控制，模拟结果为中国大陆地应力场提供了一个定量的参考模型.

中国大陆，地应力场，活动构造，黏弹性，有限元模拟

1 引 言

地应力是地震活动、造山事件、沉积盆地形成以及地下流体运动等地质过程的直接控制因素，也是采矿、水利、桥梁等工程建设中必须考虑的地质条件.中国大陆属于欧亚板块的一部分，被夹持在印度板块、菲律宾海板块、太平洋板块和西伯利亚—蒙古块体之间，受印度板块与欧亚板块的碰撞和菲律宾板块的俯冲的影响，其内部构造运动极为活跃，板块构造环境所控制的地应力场是活动构造的直接动力因素.通过活动断裂、原地应力测量和地震活动的震源机制等研究能够反映局部的地应力场状态，并在一定程度上反映区域地应力场的分布趋势特征［1］，尽管中国大陆地应力场在区域上受控于周缘的板块构造环境，而在局部地区则显示了活动断裂构造和块体构造单元地壳结构的影响，地应力场在不同的空间尺度具有不完全一致的控制因素.活动构造与地应力场的关系相互影响、密不可分，一方面，活动构造的活动性体现了局部应力场的特征；另一方面，活动构造的存在会造成局部地应力场的调整，从而影响地应力场在局部地区的分布状态［2－5］.地形起伏和地壳介质密度分布不均匀引起的重力势所产生的地应力场也是地应力分析中重要的因素之一，并引起许多学者的关注和研究［6－8］.中国大陆在不同区域地壳厚度分布不均匀，内部结构也存在很大的差异［9－16］，青藏高原的地壳厚达70多公里，而东部华北地区只有30km左右；地形地貌起伏很大，青藏高原平均高程达4km.地应力场在不同区域、不同深度其主要控制因素也有所区别.众多学者［17－19］利用数值模拟方法从总体上分析了中国大陆构造应力场的特征和成因，模拟过程中主要有以下几点考虑不够：（1）模拟过程中都采用平面直角坐标系，而对于中国大陆整体上经纬度跨度都比较大，选用球壳模型更能够接近实际；（2）在模拟过程中，边界条件和模拟结果约束或者只选用GPS，或者只选用实际地应力场观测数据，没有把二者给予统一考虑；（3）重力作为影响地应力场分布的重要因素之一，缺乏对其定量的分析.完整地理解现今地应力场的状态，不仅需要现今的构造环境和岩石圈展布特征，而且需要了解岩石圈在最近地质历史时期内的运动状态和动力学演化过程，特别是活动构造所反映的地应力场状态和演化.本文以中国大陆主要活动断裂所控制的块体模型为基础，采用黏弹性壳体有限单元，综合考虑GPS、地应力实测结果以及重力等影响，分析中国大陆地应力场的特征和控制因素.

2 基本方法

从岩石力学上来说，尽管整个地壳的岩石主要表现为弹性，但是目前在中地壳普遍存在一个岩石强度相对弱的区域.地壳上部的构造块体单元可以看作是完全弹性介质，活动断裂作为块体边界具有与构造块体本身完全不同的性质，活动断裂运动除了伴随地震发生瞬间明显的断层滑动外，更多的运动则表现为大量的无震滑动和长期蠕动等形式［20］，活动断裂带可以视为调整块体弹性应力的软弱带.

现今中国大陆地应力场是在周缘板块构造环境条件下，地应力场长期演化过程中的一个瞬间状态.地壳岩石在大时间尺度上通常可用标准线性体来近似，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

中国大陆块体通过板块边界受到相邻板块的作用力，这些作用力一方面通过完全弹性变形作用于中国大陆内部，另一方面，相邻板块持续的作用通过岩石的黏滞性对中国大陆内部作用，这种作用不同于完全弹性，在不同的时间不同块体具有不同的状态.由于受到断层库仑摩擦破裂强度的限制，块体内能够维持的地应力通常符合Byerlee准则［25］.当板块边界作用于大陆地壳内部时，板块边界的作用力不会超过Byerlee准则所确定最大应力，所以边界作用力在一定的地质历史时间内，当板块构造环境没有发生变化时，应为一个持续恒定的应力，这种边界作用力除了瞬时的弹性应力作用于块体内部，还会由于不同块体具有不同岩石黏滞性，导致在不同的时段具有不同的应力状态，并最终达到平衡，大陆地壳整体都处于一个临界的断裂平衡的状态，全球广泛存在有水库蓄水和流体注入引起的地震活动、一个地震的发生会触发另一个地震的发生，原地深孔地应力测量结果符合Byerlee准则，这些现象恰好说明了目前大陆内部，特别是活动构造断裂附近处于临界断裂平衡状态［26－27］.

地应力场通常为构造应力场与重力势引起的应力场的叠加，重力势引起的应力场主要是由地形起伏和岩石层密度不均匀造成的，当地势平缓且地下岩石分布均匀时，重力对地应力场影响可以忽略不计，构造应力场通常由于受到岩石强度的限制，在不同的构造环境具有不同的变化规律.

3 地质模型

3.1 中国大陆活动构造块体模型

中国大陆是由不同活动块体组成，块体之间以大型的活动断裂带为分界，板内块体具有不同活动程度的次级断裂［28－29］.中国大陆及其邻区的活动地块主要有拉萨、羌塘、巴颜喀拉、柴达木、祁连、川滇、塔里木、天山、准噶尔、华北等（图2）.随着现代GPS观测资料的不断积累，初步构建了中国地壳活动的运动模型［30－33］，并在此基础上分析了其应变场的特征［34－35］，中国大陆的现今构造变形主要以刚性地块运动为主，如塔里木、鄂尔多斯、华南等地块；只有少数块体内部存在比较明显的连续变形，如青藏高原和天山.块体分界的大型活动断裂区是变形的主要区域，边界断裂区的应变率远大于块体内部的应变率，通常边界断裂带都具有数十公里的宽度，发育了中国大陆的主要地震活动［29，36］.周缘的印度板块、菲律宾板块和太平洋板块的运动是中国大陆板块内部构造活动的主要动力来源.印度板块与欧亚板块为陆陆碰撞，形成挤压环境；菲律宾板块和太平洋板块与欧亚板块为洋陆碰撞，太平洋洋壳和菲律宾的洋壳向欧亚板块之下俯冲，形成了日本海和冲绳海槽，形成拉张的边界环境.

3.2 中国大陆地壳结构

中国及邻区地壳厚度分布的特点为从东到西逐渐增厚，在构造块体内部地壳厚度通常变化较缓，如蒙古高原、华北、华中、华南、青藏高原内部等地区（图2）.西部青藏高原的地壳厚度平均为60～70km，局部达到74km，而东部地区的厚度平均为30～36km，沿海地区仅有28～30km；地壳局部厚度变化与构造环境密切相关，大型地块的边界带处的地壳厚度梯度也比较大，南北构造带以东地壳厚度为40km，而在南北构造带以西达到60km，在很短的距离上地壳厚度发生了很大的变化.中国大陆的地壳在沉降带或盆地发育地方的地壳较周围薄一些，如塔里木、准格尔、柴达木、鄂尔多斯等盆地，而抬升区或造山带处的地壳较周围厚一些，如青藏高原、天山、太行山等地区［37－39］.稳定的克拉通塔里木、华南、鄂尔多斯地块处，其地壳内低速层不发育，而其它地区都普遍有地壳低速层发育，特别在青藏高原的中地壳低速层整体上普遍存在［9］.而细致地壳的流变结构成果［40］为本文进行地应力场的数值分析提供了直接的参数依据.

4 有限元模型的建立

4.1 有限元网格模型

在模拟过程中，以中国大陆块体模型为基础，采用三维球壳黏弹线性体单元，对块体模型进行无结构网格剖分（图3），模型的地形分布参照GTOPO30数字高程模型，按照区域板块构造环境，设定边界条件和加载.块体边界的大规模活动断裂带的网格单元划分相对细密，从而保证边界断裂带的数值计算的稳定性.各个块体和断裂带分层岩石物性参数见表1，各个块体分为上地壳、中地壳和下地壳三层，表中的杨氏模量和泊松比为瞬态泊松比，剪切模量系数是指岩石剪切模量在经过足够长的时间后，其长期剪切松弛模量与瞬时剪切模量的近似比值，长时间的泊松比近似为0.5，岩石密度上地壳取为2750kg／m3，中下地壳为3300kg／m3.

图2 中国大陆构造块体分区及板块构造环境Fig.2 China mainland tectonic blocks and its surrounding plate tectonics

图3 有限元网格模型及边界条件Fig.3 FEM mesh model and boundary condition

4.2 边界条件和加载

模型北边界中段为南北向和东西向均为约束边界；北边界东西两端和南边界东段均为南北向约束，东西向自由；西边界北段和东边界南端均为东西向约束，南北向自由；其他边界均按照构造应力环境加载应力.

根据Anderson断层理论可以确定不同性质的断层所在区域的构造应力特征［27，41］.在正断层、走滑断层和逆冲断层活动的不同区域，其应力特征分别为：

其中σ1为最大主应力，σ3为最小主应力，SHmax为水平最大主应力，Shmin为水平最小主应力，SV为垂向应力，Pf为孔隙流体压力，μ为断层临界摩擦系数，根据Byerlee实验数据［25］和深部地应力测量结果［27］，其取值范围在0.6～1.0之间.由此可得出正断层区域，垂直断层的方向为最小主应力方向，也是水平最小主应力，其大小约为

表1 模型岩石物性参数Table 1 Physical property parameters of rocks in model

在逆冲断层区域，垂直断层的方向为最大主应力方向，也是水平最大主应力，其大小约为

通常岩石流体密度比ρrock／ρliquid约为3.0，

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

由于板块边界的最大主应力方向与板块边界近似垂直，所以按照板块构造环境对边界加载最大主应力或最小主应力的应力载荷，除去重力作用的影响，边界条件见表2.

表2 边界加载应力表Table 2 The load on the boundary

5 模拟结果分析

图4为在边界加载压力约10万年后的近地表的速度分布，在这一时间点上，模拟速度整体上能够很好地与现今实测的GPS速度场吻合，因此，现今构造应力场的分布特征可参考该时点的应力状态模拟结果.图5为在只加载边界应力，而不考虑重力作用，约10万年后的距地表1000m深度的构造应力场模拟结果，总体上构造应力场在西部为挤压，东部为拉张，构造应力场最大主应力的方向受到板块构造环境和活动构造分布的控制，印度板块与欧亚板块的持续碰撞形成了青藏高原及其周缘的挤压性质的构造应力场，而东部菲律宾板块与太平洋板块的俯冲形成了黄海、东海和环渤海区域的拉张性质构造应力场，中间为拉张环境和挤压环境的过渡区域.

图4 模拟速度与现今实测GPS速度对比图图中白色的箭头为实测的GPS速度场，参照系为稳定欧亚大陆板块（据王琪等［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以上；第二级阶梯位于青藏高原与大兴安岭─太行山─巫山─雪峰山之间，包括内蒙古高原、黄土高原、云贵高原和塔里木盆地、准噶尔盆地、四川盆地等地区，海拔一般为1000～2000m，惟四川盆地较低，海拔在500m以下；第三级最低阶梯，在大兴安岭─太行山─巫山─雪峰山以东，自北而南，有海拔200m以下的东北平原、华北平原和长江中下游平原.重力势引起的应力场与构造应力场的结果大致相反，在青藏高原的藏东南、四川西部和云南西北地区为拉张性质的应力场，而在塔里木盆地和东部的平原区为挤压性质的应力场，应力场的方向受地势展布特征控制.图7为1000m深度位置的中国大陆地应力场，是由构造应力的模拟结果与重力作用模拟结果进行张量求和合成，可以作为浅部地应力场的定量参考模型，地壳深部的地应力场应与构造应力场相近.模拟地应力场与实测的地应力分布特征［43］（图8）也比较一致，水平最大主应力的方向主要受到板块构造环境和活动断裂带的控制，总体上以南北构造带为界，以西以较强的压性构造环境为主，以东为较弱的压性地应力环境，由于受到重力的影响，藏南和滇西北局部地区的浅部存在有张性地应力环境，而在塔里木区域，由于重力势引起的应力场与构造应力场都为挤压性质，因此挤压强度得以增加.图9为初期瞬时的模拟结果，可作为完全弹性结果，完全弹性模型的构造应力场模拟结果与实测的地应力场存在一定的偏差，如在龙门山断裂带的北段地区，水平最大主应力的方向总体为北西—南东向，而完全弹性模型的结果为北东—南西向.青藏高原南缘在重力作用下形成了张性应力环境，但是由于重力作用影响有限，在深部的张性环境应与深部物质上涌等局部构造活动有更为密切的关系.

图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 中国现代构造应力场图（谢富仁等［43］，2003）Fig.8 The measured crust stress of China（reference the Crust Stress Databese of Institute of Crustal Dynamics［43］）

中国大陆地应力场和现今活动构造的性质总体上主要受到印度洋板块的碰撞挤压和太平洋板块与菲律宾板块俯冲的影响.由于受到华南板块、塔里木克拉通以及鄂尔多斯块体的阻挡，使得印度洋板块碰撞作用集中在青藏高原及其周边，青藏高原主体多为活动断裂的走滑，块体内部的连续变形较主要断裂活动要小很多.华北块体、山西地堑除了受到印度洋板块碰撞的影响之外，菲律宾和太平洋俯冲板块所形成的弧后张性环境的影响不容忽视.

图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，所以岩石块体内部处于一个稳定的无变形状态.一方面，上地壳的变形受控于自身的构造块体结构，在板块构造环境控制下，如果不考虑其底部中下地壳的作用，在仅受板块构造环境的作用，其应力状态和变形特征应接近于完全弹性模型的数值分析结果；另一方面，当柔性较强的中下地壳受到与上地壳同样的板块构造环境的作用时，中下地壳会表现出更强的运动性，从而对上地壳的变形和应力分布产生影响.例如在青藏高原，由于中下地壳向两侧挤出，导致上地壳被中下地壳拖曳，使得上地壳在活动断裂的控制下随中下地壳运动，现今中下地壳对上地壳的变形模式和应力分布特征的影响作用不容忽视.

青藏高原东南部主应力方向出现转折，主要有以下几个因素控制：（1）由于印度板块持续的碰撞，使得青藏高原整体受到塔里木克拉通、华南板块和鄂尔多斯块体的阻挡，使得青藏高原随着活动断裂的滑动向两侧挤出［44］，在挤出过程中，刚性的上地壳块体构造的运动受边界断裂的控制；（2）由于中下地壳较上地壳有更强的柔性，所以挤出时上地壳的运动受到中下地壳拖曳的影响；（3）印度板块向北运动过程中，通过缅甸的实皆断裂（Sagaing Fault）剪切拉伸作用于欧亚板块，使得上地壳中青藏高原东南部应力和地表运动方向都发生了转向.

6 结论与讨论

本文基于线性黏弹体模型对现今中国大陆的构造应力场进行了数值模拟分析，其前提假设是中国大陆周边的板块构造环境相对稳定，且维持了足够长的时间，从而能够保证模拟分析过程中边界条件的一致性；大型活动断裂所决定的块体结构保持稳定，从而保证在模拟分析过程中地质模型的物性参数保持一致.该结果主要考虑了板块构造作用和区域性的地貌和地壳结构所产生的重力的影响，在进一步的研究中需要考虑底部地幔上涌和岩石圈热演化等影响因素，深入讨论地应力场控制因素.

整体上，中国大陆构造应力场的分布特征与活动构造的展布密切相关，是在板块构造环境控制下，整个地壳岩石黏弹特性长期演化和活动断裂无震蠕动的结果.而伴随地震发生的断层滑动引起的地应力场的瞬时或短期调整，可参考完全弹性分析的结果.陈群策等（私人通信）最近在广元附近进行的汶川地震后的地应力测量，所测的地应力水平最大主应力方向为北东—南西向，而该地区通常的区域应力场的分布特征为北西—南东向.应力场作为动态演化的过程，特别是断层附近的地应力场，由于受到断层蠕动和滑动的影响，地应力的变化会更为突出.郭启良［45］和廖椿庭［46］分别在汶川地震前后和昆仑山地震前后的地应力测量所反映的地应力的变化，也同样说明，活动断裂附近的地应力场受到局部构造的影响，从而体现了不稳定性.随着地应力监测网络的建立，可以有效地利用地应力动态的信息，更好地认识区域构造应力与局部地质构造和岩石物性之间的关系.

致 谢 感谢匿名评审人提出的宝贵意见！

（References）

［1］ 谢富仁，崔效锋，赵建涛等.中国大陆及邻区现代构造应力场分区.地球物理学报，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］ 刘峡，马瑾，傅容珊等.华北地区现今地壳运动动力学初步研究.地球物理学报，2010，53（6）：1418－1427.Liu X，Ma J，Fu 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］ 邓阳凡，李守林，范蔚茗等.深地震测深揭示的华南地区地壳结构及其动力学意义.地球物理学报，2011，54（10）：2560－2574，doi：10.3969／j.issn.0001－5733.2011.10.013.Deng Y F，Li S L，Fan 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］ 葛粲，郑勇，熊熊.华北地区地壳厚度与泊松比研究.地球物理学报，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］ 滕吉文，张洪双，孙若昧等.青藏高原腹地东西分区和界带的地球物理场特征与动力学响应.地球物理学报，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］ 刘琼林，王椿镛，姚志祥等.华北克拉通中西部地区地壳厚度与波速比研究.地球物理学报，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］ 李飞，周仕勇，苏有锦等.川滇及邻区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］ 于湘伟，陈运泰，张怀.京津唐地区地壳三维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］ 汪素云，陈培善.中国及邻区现代构造应力场的数值模拟.地球物理学报，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］ 朱守彪，石耀霖.中国大陆及邻区构造应力场成因的研究.中国科学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］ 石耀霖.关于应力触发和应力影概念在地震预报中应用的一些思考.地震，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］ 程佳，刘杰，甘卫军等.1997年以来巴颜喀拉块体周缘强震之间的黏弹性触发研究.地球物理学报，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年汶川地震近断裂区域变形场的空间分布.地球物理学报，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，Freymueller J T，et al.Present－day crustal deformation in China constrained by global positioning system measurements.Science，2001，294（5542）：574－577.

［32］ 张培震，甘卫军，沈正康等.中国大陆现今构造作用的地块运动和连续变形耦合模型.地质学报，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观测资料划分现今地壳活动块体的方法.大地测量与地球动力学，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，Freymueller 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，Fan J C.Crustal structure of mainland China from deep seismic sounding data.Tectonophysics，2006，420（1－2）：239－252.

［38］ 蔡学林，朱介寿，曹家敏等.中国大陆及邻区岩石圈地壳三维结构与动力学型式.中国地质，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］ 王琪，张培震，马宗晋.中国大陆现今构造变形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］ 谢富仁，陈群策，崔效锋等.中国大陆地壳应力环境研究.北京：地质出版社，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，FE modeling

10.6038／j.issn.0001－5733.2012.04.020

P313

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

中央级科研院所基本科研业务专项资助项目（DZLXJK201101、DZLXJK200801），中国地质调查项目（1212010911049，1212011120163）和国土资源部公益性行业科研专项（SinoProbe－06）共同资助.

范桃园，男，1971年生，2001年毕业于中国科学院研究生院获博士学位，现为中国地质科学院地质力学研究所副研究员，主要从事地球动力学、地壳应力场方面的研究.E－mail：fan3ty＠gmail.com

范桃园，龙长兴，杨振宇等.中国大陆现今地应力场黏弹性球壳数值模拟综合研究.地球物理学报，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.

（本文编辑 胡素芳）