文章快速检索     高级检索
  大地测量与地球动力学  2019, Vol. 39 Issue (5): 458-463  DOI: 10.14075/j.jgg.2019.05.004

引用本文  

侯强, 丁小军, 赵宏, 等. 青藏高原东南缘构造应力场重构[J]. 大地测量与地球动力学, 2019, 39(5): 458-463.
HOU Qiang, DING Xiaojun, ZHAO Hong, et al. Tectonic Stress Field Reconstruction in Southeastern Tibetan Plateau[J]. Journal of Geodesy and Geodynamics, 2019, 39(5): 458-463.

项目来源

中央高校基本科研业务费专项基金(CUGL120234)。

Foundation support

Special Fund for Basic Scientific Research of Central Universities, No.CUGL120234.

第一作者简介

侯强,副教授,主要从事地球动力学与统计地震学研究,E-mail:houq2008@cug.edu.cn

About the first author

HOU Qiang, associate professor, majors in geodynamics and statistical seismology, E-mail:houq2008@cug.edu.cn.

文章历史

收稿日期:2018-05-28
青藏高原东南缘构造应力场重构
侯强1     丁小军1     赵宏2     欧明霖1     张博康1     
1. 中国地质大学(武汉)机械与电子信息学院,武汉市鲁磨路388号,430074;
2. 云南大学资源环境与地球科学学院,昆明市翠湖北路2号650091
摘要:采用有限元方法,针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件,分别进行优化分析后处理,将对应台站形变模拟值与GPS实测值进行误差分析,最终选取最佳方案作为古构造应力场。结果表明,青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用,特别是印度板块NNE向强烈碰撞作用,成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源,控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源:西北部青藏高原侧向挤压造成的WE向应力约105 MPa;西南部直接来自于印度板块的NE-WS向应力约70 MPa;南部NS向作用力33 MPa;东南部扬子块体侧向NW-SE阻挡力56 MPa;北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘,形成现今复杂应力场。
关键词青藏高原东南缘古应力场有限元数值模拟

20世纪70年代Tapponnier等[1]和Molnar等[2]提出的青藏高原侧向挤出滑移运动图像及力学机制,对认识亚欧地壳运动及我国西南地区地壳块体的相互运动很有意义。鄢家全等[3]、邓起东等[4]阐述了中国构造应力场的总体特征。罗灼礼[5]分析了构造应力场和地震活动的关系。阚荣举等[6]讨论了青藏高原东南缘现代构造应力场分区及构造活动特征。许忠淮等[7]推断青藏高原东南缘平均主应力轴的方向。钟继茂等[8]根据川滇多个地震区的256次中强地震震源机制解,给出各区应力张量的定量分析结果。

地球物理探测揭示的强震孕育的发震结构和物性条件需要与介质的力学性质相联系[9]。现代构造应力研究是目前地震学和地球动力学研究的热点,但视野过于局限,只有全面了解区域应力的构造演化过程,才能充分了解地震孕育机制,预判应力发展趋势并分析未来地震危险性。因此,古应力场的研究具有重大的理论和实践意义,也是本文探讨的主要内容。

古构造应力场定量化研究方法很多,但受适用条件等因素限制,效果均不佳。本文采用有限元方法反演应力场。地质体的变形具有变形机制复杂、变形经历时间长和变形大等特点,到目前还没有建立一个能精确描述其复杂变形过程的本构方程,但在一定限制条件(如地壳浅层的短期变形)下, 地质体大体符合弹性(线性或非线性)变形特征。

本文介绍有限元数值模拟方法、数据、建模过程,详细描述后处理中的优化分析过程,给出9种不同的应力边界条件方案。对每种方案进行优化分析后处理,并将相应台站形变模拟值与GPS实测值进行误差分析,选取最佳方案作为青藏高原东南缘古构造应力场。

1 有限元数值模拟

ABAQUS是具有齐全的单元模式、丰富的材料模型以及精确的分析过程的非线性有限元分析软件,特别适用于岩土工程方面的模型分析。

创建有限元模型一般需要建立各个部件,并赋予物理材料性质,然后将其装配形成模型整体。此流程操作时存在一定弊端:2个部件相结合的边界不易吻合,导致后续未能严丝合缝地装配。为此,ABAQUS提供了草图功能,用户可将整个模型内外部所有边界预先画入草图,由于软件不检查模型的自封闭性,后面建立部件时可直接导入草图。

对青藏高原东南缘进行有限元建模前需要先弄清该区的主要断裂和块体分区。

1.1 主要断裂和块体分区

青藏高原东南缘发育多条大型断裂带,围绕菱形块体东侧的是玉树-鲜水河-小江断裂带,包括甘孜-玉树断裂、鲜水河断裂、安宁河断裂、则木河断裂以及小江断裂,属于一级构造边界带;菱形块体西侧边界由澜沧江断裂和红河断裂带构成;同时,内部发育的小金河-丽江断裂将菱形块体一分为二,周围的大断裂则主要有龙门山断裂带和南汀河断裂带。本文将龙门山断裂带北侧的岷江断裂也考虑进建模过程中,因为从速度场结果来看,速度值在岷江断裂附近明显减小,说明两侧地块的物性存在较大差别,应区分对待。

图 1标注了青藏高原东南缘的主要活动断裂带,这些断裂将整个地区划分为5个一级块体:马尔康块体、川滇菱形块体、扬子块体、滇西南块体和拉萨块体。其中,岷江断裂又将马尔康块体划分为松潘-甘孜块体和秦岭块体2个次级地块,小金河-丽江断裂则将菱形块体分为北侧的滇西北块体和南侧的滇中块体[10]

Ⅰ马尔康块体; Ⅰ1松潘-甘孜块体; Ⅰ2秦岭块体; Ⅱ川滇菱形块体; Ⅱ1滇西北块体; Ⅱ2滇中块体 图 1 青藏高原东南缘主要断裂带及块体分布 Fig. 1 The major fault zones and block partition in southeastern Tibetan plateau

对照图 1,选取分区的主要控制点,要求既能反映模型主要特征,又不可过于简略以致与实际不符。连接各控制点完成草图,之后创建部件时导入并逐个建立块体,同时将垂向拉伸长度设置为20 km,即上地壳平均厚度(模型认为断裂带主要分布于上地壳)。模型中除龙门山断裂带以外均设置为直立断层。考虑龙门山断裂带浅陡深缓特征,上地壳浅部倾角约60°~80°,深部则放缓至30°~40°,故设置龙门山断裂倾角为70°[11]

1.2 介质参数

建模时选用线弹性模型,主要涉及的物理量为杨氏模量和泊松比。设置断裂带则需要对每个块体进行介质分区,沿断裂带在块体内部作宽度10 km的平行线,将这10 km区域作为弱化带处理,其杨氏模量为所在块体杨氏模量的1/3。各块体介质参数如表 1所示[9, 12-13],其中滇中块体杨氏模量较其他块体偏大,原因可能为滇中块体与扬子块体本属同一块体。

表 1 块体介质参数 Tab. 1 The medium parameters of blocks
1.3 网格划分

ABAQUS划分网格有多种策略,首先需要在部件上建立种子来控制网格节点的密度。本文选择全局种子尺寸为30 km,断裂带附近由于精度要求高,故将断裂附近边界种子尺寸适当缩小。选择自由网格划分技术,通过四面体单元来划分网格,共建立61 833个网格、31 822个节点。

1.4 分析步及相互作用

一般来说可建立3个分析步:首先是系统默认的初始步,在该步中可添加模型的初始应力场,以及一直存在的载荷或约束;通常会在第2个分析步建立接触,使模型能够在接触稳定后过渡至后面的步骤,否则可能会出现一些不稳定或是不收敛的解;最后加载边界条件,此时将初始时间增量定为总时间的1/10,从而使载荷平滑地加载到模型,有效减小计算不收敛发生的概率。

图 2 青藏高原东南缘上地壳三维有限元模型 Fig. 2 Three dimensional finite element model of upper crust in southeastern Tibetan plateau

图 3 模型应力边界条件 Fig. 3 Stress boundary conditions of the model

定义相互作用时首先选定作用属性为“接触”,接下来添加切向行为,即沿接触面的行为,选用罚函数摩擦模型,并将摩擦系数设置为0.6[14]。选择模型的主从接触面共9对,勾选“仅调整至消除过盈量”,从而使得主从接触面在分析开始时能够精确接触。本文以4 Ma时瞬时平衡态的位移场和应力场图像作为计算结果进行分析。

1.5 应力边界条件

中国大陆现代构造应力场主要起源于中国大陆周围板块间的相互作用。印度板块NNE向的强烈碰撞作用是中国大陆构造应力场最主要的动力源,它控制了中国大陆构造应力场和构造运动的总体特征。对青藏高原东南缘来说,西北部受到青藏高原隆起及其向东物质挤出的作用;西南部受印度板块的挤压作用;东部受到华南块体的阻挡作用。其中,稳定的四川盆地的阻挡作用导致物质东流的方向发生偏转,从向东转为向东南。青藏高原东南缘的应力场大致存在3个分区:27°N以北的区域应力场方向为近EW向;27°N以南、100°E以西的区域应力场方向为NE-SW向;27°N以南、100°E以东的区域应力场方向为NW-SE向。动力边界条件的具体数据参考前人研究。陈化然等[15]根据震源机制解给出具体数据为:F1=105 MPa,F2=70 MPa,F3=33 MPa,F4=66 MPa,F5=100 MPa。本文据此对5个参数进行调节,确定最佳边界条件作为青藏高原东南缘的古构造应力。模型上表面为自由表面,底部垂直方向固定,水平方向自由。

2 模拟结果后处理

将场输出请求(field output request)设置为我们关心的物理量,创建相应模型的Job并提交给软件进行计算,即可得到模拟结果。本次模拟运行内存为8 G,计算时间约10 min。运行完后,使用Visualization功能即可将计算结果进行可视化处理。

由于本文需要从形变模拟结果中提取出位于研究区域内81个GPS测站所对应的位移模拟值进行匹配,故选用Python语言进行场输出数据的后处理。ABAQUS有限元程序通过集成Python脚本向二次开发用户提供丰富的库函数,可直接操纵ABAQUS内核,实现建模、划分网格、指定材料属性、提交作业、后处理分析等功能。

后处理的第一步是建立模型结果文件odb,然后从odb对象中读取step的最后一步,并从中得到最后一帧(frames[-1])的fieldOutputs场输出结果数据。该数据中包含创建场输出时勾选的所有变量结果,这里选择“U”,即形变结果。最后循环读取每个实例中每个节点的东向形变U1和北向形变U2。需要注意的是,ABAQUS是按照块体对节点单独编号的,而不是整体编号,所以需要清楚每个块体有多少个节点,GPS测站所在位置位于哪个实例块体的哪个节点编号上,这样才能准确读取对应台站的模拟值并进行后续匹配。本文模型中,松潘-甘孜块体共10 732个节点,其中有12个对应GPS台站;秦岭块体共5 724个节点,其中有3个对应GPS台站;滇西北块体共4 624个节点,其中有4个对应GPS台站;滇中块体共6 363个节点,其中有11个对应GPS台站;扬子块体共17 437个节点,其中有22个对应GPS台站;滇西南块体共6 789个节点,其中有9个对应GPS台站;拉萨块体共10 164个节点,其中有10个对应GPS台站。利用经纬度换算确定每个GPS台站对应的InstanceName和NodeID,然后使用Python查找这2个条件下对应点的位移值U1U2,保存在文件中,方便后续作误差分析。

3 边界条件误差分析

前面提到,本文对5个参数F1~F5进行调节,确定最佳边界条件作为青藏高原东南缘的古构造应力。调节的原则是增大或减小10 MPa,由此得到以下9种应力边界条件选择方案。对于每一种方案,采用上面的方法提取出对应台站的模拟值abaqus_U1i和abaqus_U2i(i=1, 2, …70),对应GPS测站的真实值为gps_U1i和gps_U2i(i=1, 2, …70),由此可以计算形变大小模拟值abaqus_Ui和实测值gps_Ui,形变方向模拟值abaqus_Uazi和实测值gps_Uazi,利用最小二乘原理计算二者之间的误差,公式如下:

$ \delta U = \\ \sqrt {{{\left( {{\rm{abaqus}}\_{U_1}-{\rm{gps}}\_{U_2}} \right)}^2} + {{\left( {{\rm{abaqus}}\_{U_2}-{\rm{gps}}\_{U_2}} \right)}^2} + \cdots }, \\ \;\left( {i = 1, 2, \cdots, 70} \right) $ (1)
$ \delta {U_{{\rm{az}}}} =\\ \sqrt {{{\left( {{\rm{abaqus}}\_{U_{{\rm{a}}{{\rm{z}}_1}}}-{\rm{gps}}\_{U_{{\rm{a}}{{\rm{z}}_2}}}} \right)}^2} + {{\left( {{\rm{abaqus}}\_{U_{{\rm{a}}{{\rm{z}}_2}}}-{\rm{gps}}\_{U_{{\rm{a}}{{\rm{z}}_2}}}} \right)}^2} + \cdots },\\ \;\left( {i = 1, 2, \cdots, 70} \right) $ (2)
4 讨论

表 2可知,方案4给出的应力边界条件得到的模拟结果与GPS实测值最为匹配,形变大小误差为1.432 8 MPa,形变方向误差为16.372 1°。图 4展示了该边界条件下模拟得到的位移场以及与GPS台站实测值的对比图,图 4(a)中箭头方向代表位移矢量方向,长短以及颜色均代表位移大小;图 4(b)为各基准站模拟结果与GPS实测值的对比,黑色箭头为GPS实测值,红色为模拟值。对比二者可以看出,方案4的数值模拟结果较好地反映了青藏高原东南缘西北部青藏高原推挤、东侧扬子块体阻挡的格局,同时位移方向在菱形块体及拉萨块体内部绕东构造结顺时针旋转的特性。图 4(b)中大部分站点在大小和方向上吻合较好,仅北侧松潘-甘孜块体内部以及南侧滇西南块体内部少数存在一定差异,其原因可能是这些地块内部存在某些特殊构造[14]

表 2 应力边界条件方案及误差 Tab. 2 The medium parameters of blocks

图 4 形变场模拟结果 Fig. 4 The simulated result of displacement field

ABAQUS中绝对主应力最大值的正值代表主张应力,负值代表主压应力。观察青藏高原东南缘最大主应力方向,同样在“Y”字形结构下方呈现高压区;拉萨块体东侧与小金河-丽江断裂带西段交界区域存在一定东西向压力,通过澜沧江断裂带在北侧形成近南北向张力。滇西南块体整体呈现拉张趋势,与拉萨块体交界处的南汀河断裂带西段北侧承受着较大的压力,在通过断裂带挤压过程中转化为张力;红河断裂带中段存在来自滇中块体的南向压力,虽然数值较小,但仍然不容忽视。图 5(b)给出了利用时空应力场反演(spacial and temporal stress inversion, SATSI)程序以及青藏高原东南缘1 131个震源机制解得到的该区现今应力图像。可以看出,青藏高原东南缘西北部受到来自青藏高原NEE向挤压作用;东部则有扬子块体阻挡产生的NWW、NW向应力,该阻挡作用导致来自高原的E向应力在菱形块体及其西侧发生顺时针旋转,菱形块体内部呈现北部SE向、南部SSE向应力分布的情况;另一方面,西南部由于印度板块持续推挤产生NE向应力。

图 5 应力场模拟结果 Fig. 5 The simulated result of stress field
5 结语

分析认为,青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用,特别是印度板块NNE向强烈碰撞作用,成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源,控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括以下几个力源:西北部青藏高原侧向挤压造成的WE向应力约105 MPa;西南部直接来自于印度板块的NE-WS向应力约70 MPa;南部NS向作用力33 MPa;东南部扬子块体侧向NW-SE阻挡力56 MPa;北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘,经过数百万年的不断演化和构造运动,最终形成现今的复杂应力场。

本文关于应力边界条件的初始设定来自其他文献,仅分为5个区域,未来可参考更多研究成果同时加入程序计算以提高精确度。同时,调整应力作用时,仅考虑了古应力的大小,未调整方向,未来希望在二次开发过程中不断调整其大小和方向,实现结果的自动匹配和优化分析,从而得出更为精确、可靠的古应力场。

参考文献
[1]
Tapponnier P, Molnar P. Slip-Line Field Theory and Large-Scale Continental Tectonics[J]. Nature, 1976, 264(5584): 319-324 (0)
[2]
Molnar P, Tapponnier P. Cenozoic Tectonics of Asia: Effects of a Continental Collision: Features of Recent Continental Tectonics in Asia Can Be Interpreted as Results of the India-Eurasia Collision[J]. Science, 1975, 189(4201): 419-426 (0)
[3]
鄢家全, 时振梁, 汪素云, 等. 中国及邻区现代构造应力场的区域特征[J]. 地震学报, 1979, 1(1): 9-24 (Yan Jiaquan, Shi Zhenliang, Wang Suyun, et al. Some Features of the Recent Tectonic Stress Field of China and Environs[J]. Acta Seismologica Sinica, 1979, 1(1): 9-24) (0)
[4]
邓起东, 张裕明, 许桂林, 等. 中国构造应力场特征及其与板块运动的关系[J]. 地震地质, 1979, 1(1): 11-22 (Deng Qidong, Zhang Yuming, Xu Guilin, et al. On the Tectonic Stress Field in China and Its Relation to Plate Movement[J]. Seismology and Geology, 1979, 1(1): 11-22) (0)
[5]
罗灼礼. 印度洋板块向北顶撞与我国西部及其邻区现代构造应力场和地震活动的关系[J]. 地震工程学报, 1979, 1(4): 11-21 (Luo Zhuoli. The Relationship between the Northward Collision of the Indian Ocean Plate and the Modern Tectonic Stress Field and Seismic Activity in the Western Part of China and Its Adjacent Areas[J]. Journal of Earthquake Engineering, 1979, 1(4): 11-21) (0)
[6]
阚荣举, 张四昌, 晏凤桐, 等. 我国西南地区现代构造应力场与现代构造活动特征的探讨[J]. 地球物理学报, 1977, 20(2): 96-109 (Kan Rongju, Zhang Sichang, Yan Fengtong, et al. Present Tectonic Stress Field and Its Relation to the Characteristics of Recent Tectonic Activity in Southwest China[J]. Chinese Journal of Geophysics, 1977, 20(2): 96-109) (0)
[7]
许忠淮, 汪素云, 黄雨蕊, 等. 由多个小震推断的青甘和川滇地区地壳应力场的方向特征[J]. 地球物理学报, 1987, 30(5): 476-486 (Xu Zhonghuai, Wang Suyun, Huang Yurui, et al. Directional Characteristics of the Crustal Stress Field in the Southeastern Margin of Qinggan and Qinghai-Tibet Plateau Inferred from Multiple Small Earthquakes[J]. Chinese Journal of Geophysics, 1987, 30(5): 476-486 DOI:10.3321/j.issn:0001-5733.1987.05.005) (0)
[8]
钟继茂, 程万正. 由多个地震震源机制解求川滇地区平均应力场方向[J]. 地震学报, 2006, 28(4): 337-346 (Zhong Jimao, Cheng Wanzheng. Determination of Directions of the Mean Stress Field in Sichuan-Yunnan Region from a Number of Focal Mechanism Solutions[J]. Acta Seismologica Sinica, 2006, 28(4): 337-346 DOI:10.3321/j.issn:0253-3782.2006.04.001) (0)
[9]
马宏生.川滇地区强震孕育的深部动力环境研究[D].北京: 中国地震局地球物理研究所, 2007 (Ma Hongsheng. Study on the Deep Dynamic Environment of Strong Earthquakes of Sichuan-Yunnan Region[D]. Beijing: Institute of Geophysics, CEA, 2007) http://cdmd.cnki.com.cn/Article/CDMD-85401-2007205365.htm (0)
[10]
徐锡伟, 闻学泽, 郑荣章, 等. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科学D辑:地球科学, 2003, 33(增1): 151-162 (Xu Xiwei, Wen Xueze, Zheng Rongzhang, et al. The Latest Tectonic Change Pattern and Its Dynamic Source of Active Blocks in the Sichuan-Yunnan Region[J]. Science in China Series D:Earth Sciences, 2003, 33(S1): 151-162) (0)
[11]
张竹琪, 张培震, 王庆良. 龙门山高倾角逆断层结构与孕震机制[J]. 地球物理学报, 2010, 53(9): 2068-2082 (Zhang Zhuqi, Zhang Peizhen, Wang Qingliang. The Structure and Seismogenic Mechanism of Longmenshan High Dip-Angle Reverse Fault[J]. Chinese Journal of Geophysics, 2010, 53(9): 2068-2082 DOI:10.3969/j.issn.0001-5733.2010.09.007) (0)
[12]
杨兴悦, 陈连旺, 杨立明, 等. 巴颜喀拉块体强震动力学过程数值模拟[J]. 地震学报, 2013, 35(3): 304-314 (Yang Xingyue, Chen Lianwang, Yang Liming, et al. Numerical Simulation on Strong Earthquake Dynamic Process of Bayan Har Block[J]. Acta Seismologica Sinica, 2013, 35(3): 304-314 DOI:10.3969/j.issn.0253-3782.2013.03.003) (0)
[13]
李红.构造应力场、活动断裂及区域地震活动性的数值模拟研究[D].北京: 中国地震局地壳应力研究所, 2008 (Li Hong. Numerical Simulation on Tectonic Stress Field, Active Fault and Regional Seismicity[D]. Beijing: Institute of Crustal Dynamics, CEA, 2008) (0)
[14]
李玉江, 陈连旺, 陆远忠, 等. 汶川地震的发生对周围断层稳定性影响的数值模拟[J]. 地球科学——中国地质大学学报, 2013, 38(2): 398-410 (Li Yujiang, Chen Lianwang, Lu Yuanzhong, et al. Numerical Simulation of the Influences of Wenchuan Earthquake on the Stability of Faults in the neighborhood[J]. Earth Science—Journal of China University of Geosciences, 2013, 38(2): 398-410) (0)
[15]
陈化然, 陈连旺, 马宏生, 等. 川滇地区应力场演化与强震间相互作用的三维有限元模拟[J]. 地震学报, 2004, 26(6): 567-575 (Chen Huaran, Chen Lianwang, Ma Hongsheng, et al. 3-D Finite Element Modeling for Evolution of Stress Field and Interaction among Strong Earthquakes in Sichuan-Yunnan Region[J]. Acta Seismologica Sinica, 2004, 26(6): 567-575 DOI:10.3321/j.issn:0253-3782.2004.06.001) (0)
Tectonic Stress Field Reconstruction in Southeastern Tibetan Plateau
HOU Qiang1     DING Xiaojun1     ZHAO Hong2     OU Minglin1     ZHANG Bokang1     
1. Faculty of Mechanical and Electronic Information, China University of Geosciences, 388 Lumo Road, Wuhan 430074, China;
2. School of Resource Environment and Earth Science, Yunnan University, 2 North-Cuihu Road, Kunming, 650091, China
Abstract: Using finite element method on the southeast margin of the Tibetan plateau, this paper establishes a more detailed and accurate three dimensional finite element model of elastic, selects the stress of the nine different boundary conditions, applies optimization analysis for each project in post-processing, analyzes the corresponding station deformation GPS simulated values and measured error, and finally chooses the best solution as palaeotectonic stress field. Results show that around four million years ago the ancient stress field of the southeast margin of Tibetan plateau mainly originated from intense collision of the Indian plate NNE with the Chinese mainland plate. The mainland of China, especially the Tibetan plateau in southwest region, is the main power source of the tectonic stress field, and controls how the various blocks interact and their motion patterns. The paleostress field on the southeastern edge of the Tibet plateau mainly includes the following sources: The NE-WS directional stress in the southwest is about 70 MPa directly from the Indian plate. Southern NS force 33 MPa; lateral NW-SE blocking force of southeast Yangtze block is 56 MPa. At the same time, the strong EW blocking force of the Yangtze block in the eastern north is 90 MPa, and these forces act together on the southeastern edge of the Tibetan plateau, forming a complex stress field.
Key words: southeastern Tibetan plateau; paleotectonic stress field; finite element method