地球物理学报  2018, Vol. 61 Issue (5): 1817-1823   PDF    
龙门山断裂带形成过程及其地应力状态模拟
师皓宇1,2, 马念杰1, 马骥1     
1. 中国矿业大学(北京)资源与安全工程学院, 北京 100083;
2. 华北科技学院安全工程学院, 北京 101601
摘要:以龙门山附近区域运动特征以及深部岩体力学特性为基本条件,采用FLAC模拟软件计算模拟分析了龙门山断裂形成过程和塑性区分布的几何特征.研究结果认为:汶川—茂县断裂F1、映秀—北川断裂F2和灌县—安县断裂F3是从中部开始发育,塑性区上侧逐渐发展至地表;F1断层首先出现,以后才依次形成F2和F3断层.这些塑性区的发育方向与水平面的夹角大致呈40°左右,并出现了各自的共轭断层,与实际的龙门山断裂大致吻合,说明区域构造运动是龙门山断层产生的主因.论文还揭示了区域构造运动过程中区块内主应力的变化规律,主要分为2个阶段,一是在区块挤压初期,主应力随着区块位移量的增大而增大;二是主应力相对稳定阶段,区块岩体内的主应力保持相对稳定的状态,达到临界破坏岩体的最大最小主应力比值一般介于3.31~4.39.
关键词: 板块运动      龙门山断裂带      形成过程      数值模拟      地应力状态     
Numerical simulation for the formation process of the Longmenshan fault zone and its crustal stress state
SHI HaoYu1,2, MA NianJie1, MA Ji1     
1. College of Resources and Safety Engineering, China University of Mining and Technology(Beijing), Beijing 100083, China;
2. North China Institute of Science and Technology, Beijing 101601, China
Abstract: Based on characteristics of tectonic movement around the Longmenshan Mountains and mechanical properties of deep rock mass, we applied the software FLAC to analyze the formation process of Longmenshan fault zone and the geometrical distribution characteristics of plastic belts. The study suggests that the Wenchuan-Maoxian fault F1, Yingxiu-Beichuan fault F2 and Guanxian-Anxian fault F3 began to shape from the middle, then the upper plastic region gradually developed to the surface. The fault F1 first appeared, successively followed by F2 and F3. The angle between these plastic zones's development directions and the horizontal plane was approximately 40 degrees, with their own conjugate faults which were roughly consistent with the actual Longmenshan fault zone. It indicates that regional tectonic movement was the main cause for the generation of the Longmenshan fault. This study also revealed the variation law of block's principal stress during regional tectonic movement. It included 2 stages. Firstly, in the early stage of block compression, principal stress increased with the growing block's displacement. Secondly, the principal stress remains relatively stable, while the ratio of the maximum to minimum principal stress as the critical failure rock was generally between 3.31~4.39.
Key words: Plate motion    Longmenshan fault zone    Formation process    Numerical simulation    Crustal stress state    
0 引言

板块的绝对或相对运动控制着板块边界的断裂带形成和应力场特征.印度洋板块与欧亚板块的碰撞-挤压作用,形成了喜马拉雅弧形造山带,其东部弧顶沿N—NE方向插入青藏高原的东北缘(Molnar and Tapponnier, 1975England,1986England and Houseman, 1988Harrison et al., 1992付碧宏等,2008),造成了巴颜喀拉块体和龙门山断裂系深、浅部构造强烈活动和变形,而高原地壳、地幔物质以大型走滑断裂为通道边界向E—ES方向运移(滕吉文等,2014),导致青藏高原持续隆升,致使周缘发生激烈的地壳形变和频繁的大地震活动等构造运动,并控制了中国西部乃至其以北的广大地区(Van Der Rob et al., 1997Borwn et al., 2002),如图 1所示.根据火山岩化学全分析资料和GPS监测计算研究现代板块边界的运动速度,印支、燕山与四川期板内平均缩短速度均在5.4~5.6 cm·a-1,喜马拉雅期在中国西部平均缩短速度达5.2 cm·a-1,新构造期仅为2.5~2.8 cm·a-1,喜马拉雅地区表面运动速度都呈现N—NE方向运动(Sugisaki,1976万天丰,1993万天丰和任之鹤,1999叶正仁和王建,2004张培震,2008).

图 1 龙门山活动断裂与地震断层空间分布图(据付碧宏等, 2008) Fig. 1 Spatial distribution map of the active faults and seismic fault in the Longmenshan fault zone (Fu et al., 2008)

板块之间的相互运动,必然造成板块岩体应力的变化,通过“世界应力图”发现全球大部分地区的最大水平主应力方向与板块绝对运动迹线保持较好的一致性,反映出构造应力与板块运动的关系密切(谢富仁等,2003).沿着喜马拉雅地区的印度洋板块和欧亚板块的碰撞运动产生强烈的挤压构造应力(徐纪人等,2008),当岩体应力超过其强度极限时必然进入塑性状态(徐芝纶,1980沈明荣和陈建峰,2006).以上研究表明,地壳板块之间存在相互运动,且运动方向与最大主应力方向具有较好的相关性,岩体的主应力状态决定着岩体的弹塑性状态.但板块塑性区发育的几何特征与主应力之间有何关系?塑性区的发育与断层的形成有何关系?这些问题仍需要我们深入研究.龙门山断裂带是在一定深处收敛成为一条NE方向的深部剪切断裂带(滕吉文等,2014),本文拟以板块运动为研究的基本条件,依据龙门山附近区域的构造运动特性以及深部岩体力学特性,采用FLAC模拟软件模拟龙门山断裂形成过程和塑性区分布的几何特征,研究区域构造运动对于龙门山断裂带形成及地应力状态的影响.

1 龙门山地形与断裂带特征

龙门山地处青藏高原东缘与四川盆地的交接部位, 山-盆构造地貌反差强烈, 是青藏高原周边地形梯度变化最大的地区, 在不到50 km的范围内地形高差达约4000 m(邓起东等,1994; 苏生瑞等,2002付碧宏等,2008颜照坤等,2014),如图 2所示.龙门山断裂带主要由大致平行的3条逆冲、逆冲兼走滑性质的断裂组成,分别为:汶川—茂县断裂F1、映秀—北川断裂F2和灌县—安县断裂F3,该断裂带长500多公里,宽30~50 km,大致呈北东—南西向延伸构造(邓起东等, 1994; Burchfiel et al, 1995; 陈国光等, 2007).

图 2 龙门山断裂地形地貌剖面图(颜照坤等,2014) Fig. 2 Topographic profile of the Longmenshan fault zone (Yan et al., 2014)

青藏高原Moho面的深度为60~70 km, 四川盆地Moho面的深度为40 km左右,高原内部物质在推挤压力与重力的作用下向四周运移, 出现了高原物质东流现象(曾融生等, 1992; Royden et al., 1997; 朱守彪和张培震,2009),具有流变性质的下地壳速度V2大于上地壳速度V1(朱守彪和张培震,2009),如图 3所示; 根据GPS观测结果显示, 龙门山断裂带附近地区位移速度差约为2 mm·a-1(赖锡安等,2004王连捷等,2009).

图 3 龙门山断裂带运动模型示意(朱守彪等,2009) Fig. 3 Motion model of the Longmenshan fault zone (Zhu et al., 2009)
2 龙门山断裂带的数值模拟分析 2.1 模型的几何尺寸与边界条件

本文选用FLAC3D5.0数值模拟软件,它是基于拉格朗日差分法的一种显式有限差分程序(Fast Lagrangian Analysis of Continua),内置了摩尔-库仑准则,适宜于处理大尺度、大变形工程和地质问题(杨铭键等,2012王涛等,2015).

计算机数值模拟分析模型的主要平面与龙门山断裂带走向垂直,几何尺寸如图 4所示.走向长度为160 km,包含马尔康至四川盆地的部分区域;高度为40 km,与龙门山断裂带下部的Moho面大致相同;厚度取1 km.

图 4 龙门山断裂带数值计算模型 Fig. 4 Numerical calculation model of Longmenshan fault

模型的边界条件和加载方式如图 4所示.上部为自由边界,地形的起伏简化为折线;底部大体与Moho面的深度相同,简化为位移边界条件,在x方向可以运动,z方向为固定支座,即w=0;右侧是板块运动相对静止的边界,简化为在z方向可以运动,其他方向为固定铰支座,即u=v=0.在模型左侧边界施加位移前,对整个模型范围内的岩体施加重力应力场,使模型具备初始应力环境.

模型左侧的移动量按加载步逐步施加,即每个加载步施加1a的移动量.模型的左侧边界相对于右侧可以移动,移动速度取两个边界实际移动速度的差值;考虑到地表与Moho面处的移动速度的不同,在左侧边界上的水平速度根据深度的增加成线性增长,模型上部加载的正向水平速度2 mm·a-1,底部加载的正向水平速度4 mm·a-1.

为研究模型不同位置的主应力变化情况,在模型的水平位置中部(80 km)处,地表下1 km、11 km、21 km、31 km位置处设置最大最小主应力追踪测点.

2.2 地壳岩体物理力学性质

岩石破坏选用摩尔-库仑准则,其基本的物理力学参数包括弹性模量、泊松比、内聚力、内摩擦角、剪胀角和抗拉强度.

迄今为止,对地壳深部岩石力学性质的研究不是很多,仅限于几公里范围内(苏生瑞等,2002).一般认为,随着深度的增加,岩体更加致密,岩体的密度和弹性模量等都会随深度的增大而增大(蔡美峰和王双红,1997秦向辉等,2012柳畅等,2014; 熊魂等,2015).本文数值模拟分析模型的岩体弹性模量参考花岗岩取值,为4.0×1010~10.6×1010Pa(沈明荣和陈建峰, 2006),对整个模型按均匀梯度对弹性模量和密度进行赋值,其他参数为定值,物理力学参数如表 1所示.

表 1 模型岩体物理力学参数 Table 1 Physical parameters of rock for the model
3 龙门山断层形成过程

图 5描述了龙门山断裂带3条主要断裂的形成过程.由该图可见,区域构造运动6万年之前,龙门山地表出现塑性破坏区,地壳下方岩体保持完好,处于弹性状态;区域构造运动6.39万年时,在马尔康下方深度约32.7 km、水平位置为29 km处开始有塑性区产生,此处为F1断层发育初始位置;到6.72万年时,其右上方的塑性区发育较为明显,F1断层初露端倪;区域构造运动7万年时,F1断层发育至地表,在莫霍面附近呈现不对称共轭塑性区;到10.07万年时,地表下方22.0 km、水平57 km处出现塑性区,此处为F2断层发育初始位置;到10.25万年时,F2断层由中部向两侧发育;到11万年时,F2断层发育至地表,F2的底部同时发育不对称共轭断层;到14.51万年时,在深度为24.2 km、水平69.9 km位置处出现塑性区,此处为F3断层发育初始位置;到15万年时,塑性区继续发育;到16万年时,F3断层基本发育至地表;到17万年时,其底部同时发育了不对称共轭断层.

图 5 龙门山断裂形成过程模拟图 Fig. 5 Simulation of formation process for the Longmenshan fault zone

因此可认为,F1、F2、F3断层是从中部向上下两侧开始发育,上侧逐渐发展至地表,而且F1、F2和F3断层依次形成断层.这些塑性区的发育方向与水平面的夹角大致呈40°左右,并出现了各自的共轭断层,与实际的龙门山断裂大致吻合(胡幸平等,2012).

图 5还可以看出,经过长期的区域构造运动作用,在区块岩体内形成了纵横交错的塑性区,它们随着区域构造运动而不断发生变化,新的塑性区和断层不断产生.这种断层不断产生的过程也许会引发不同的地震,当区域构造运动活跃时,产生的地震会增多,反之亦然;同时,随着区域构造运动,地表的破坏深度也具有很大差别,破坏深度大的区域或许形成谷地,而破坏深度小的地域成为山峰.

4 断层形成过程地应力状态分析

提取测点1~测点4的最大最小主应力数值,处理为随时间变化的曲线图. 图 6a为测点1的最大最小主应力曲线,区块运动4万年之前, 最大最小主应力基本随时间的增大而增大,即随着区域构造运动的持续,主应力随之增大,符合岩体应力随着应变量的增大而增大的规律.当区域板块运动到一定程度时,测点的主应力数值达到相对稳定的状态,但达到稳定状态的时间不同,如图 6b图 6d所示,测点1~测点4的主应力达到稳定状态的时间大致为4万年、10万年、12万年、17万年.地壳岩体水平挤压会导致水平应力的增大,而地壳岩体变形、塑性破坏等又会导致地应力的释放,因此,地壳岩体内部因区域构造运动而增加的主应力与释放的主应力可能形成相对平衡的状态,该平衡机制类似于自组织临界性,即一个具有持续能量供给且由许多基本单元组成的系统,当单元之间具有相互作用时,会自发地演化到某个临界状态(Bak et al., 1987易桂喜等,2013).

图 6 不同测点最大(σ1)、最小(σ3)主应力随时间变化曲线图 (a)深度为1 km; (b)深度为11 km; (c)深度为21 km; (d)深度为31 km. Fig. 6 Curves of maximum (σ1) and minimum (σ3) principal stresses varying with time at different sites (a) Depth=1 km; (b) Depth=11 km; (c) Depth=21 km; (d) Depth=31 km.

为研究断层形成前、F1断层贯通、F2断层贯通和F3断层贯通等几个关键时间点的主应力变化规律,截取6万年、7万年、11万年和16万年时的4个测点的最大最小主应力数值,计算其比值如表 2所示.处于不同深度的各测点最大最小主应力虽然数值变化较大,但比值变化较小,测点1在6万年~16万年的最大最小主应力比值介于4.14~4.39,平均为4.3.深部测点也大致呈同样的规律,测点2、测点3和测点4的最大最小主应力比值变化幅度同样较小,平均比值分别为3.41、3.43和3.39.

表 2 测点主应力及其比值汇总表 Table 2 Summary table on principal stress and its ratio of measuring point

根据蝶形塑性区理论(Zhao et al., 2013马念杰等,2015),塑性区的发生发展主要是由偏应力所决定,偏应力能控制岩体的破坏,当最大最小主应力比值η达到某一临界条件(K1)后,圆孔(或地质变异体,或最先出现的塑性体)周边的塑性区呈“X”形共轭分布,塑性区半径Rmax增长趋势会发生急剧扩展现象.从本模型计算的结果来看,模型的最大最小主应力比值η在3.31~4.39之间,表明持续的水平运动,可能会致使最大最小主应力的持续变化,当其比值达到临界状态时,塑性区的发展将发生无限扩展,这也是断层从深部发育至地表的原因之所在.

5 结论

(1) 板块水平运动是地壳主应力变化的主要原因,岩体内应力增大至岩体强度极限时,必将产生塑形破坏.在一定应力环境中岩体的塑形区发育具备2个特性:一是方向特性,即断层的发育延伸具有一定的方向性,本计算模型的塑性区发育倾角约为40°左右;二是共轭特性,在大塑性区附近会形成近似垂直相交的不对称共轭塑性区.

(2) F1、F2、F3断层的发育是从中部开始,逐渐发展至地表,F1断层首先出现,F2和F3断层依次出现.模拟形成的断层几何特征与实际的龙门山3条断层大致吻合.

(3) 区块内主应力变化分为2个阶段,一是在区块挤压初期,主应力随着区块内应变量的增大而增大;二是主应力相对稳定阶段,当区块挤压持续进行时,区块岩体内的主应力保持相对稳定的状态,表明区块岩体内部存在应力积聚-应力释放的平衡机制.

(4) 不同深度岩层中的最大最小主应力比值均存在临界状态,浅部较大而深部略小,较高的偏应力环境可能是断层形成的重要原因.

致谢

感谢本文的评审专家提出的中肯建议,使得文章内容更加完善.

参考文献
Bak P, Tang C, Wiesenfeld K. 1987. Self-organized criticality: An explanation of the 1/f noise. Physical Review Letters, 59(4): 381-384. DOI:10.1103/PhysRevLett.59.381
Brown E T, Bendick R, Bourlès L D, et al. 2002. Slip rates of the Karakorum fault, Ladakh, India, determined using cosmic ray exposure dating of debris flows and moraines. Journal of Geophysical Research:Solid Earth, 107(B9): 2192. DOI:10.1029/2000JB000100
Burchfiel B C, Chen Z L, Liu Y P, et al. 1995. Tectonics of the Longmen Shan and adjacent regions, Central China. International Geology Review, 37(8): 661-735. DOI:10.1080/00206819509465424
Cai M F, Wang S H. 1997. Relation between ground stress behavior and properties of surrounding rock. China Mining Magazine, 6(6): 38-41.
Chen G G, Ji F J, Zhou R J, et al. 2007. Primary research of activity segmentation of Longmenshan fault zone since late-Quaternary. Seismology and Geology, 29(3): 657-673.
Deng Q D, Chen S F, Zhao X L. 1994. Tectonics, seismicity and dynamics of Longmenshan mountains and its adjacent regions. Seismology and Geology, 16(4): 389-403.
England P C. 1986. Finite strain calculations of continental deformation: 2. Comparison with the India-Asia collision zone. Journal of Geophysical Research: Solid Earth, 91(B3): 3664-3676. DOI:10.1029/JB091iB03p03664
England P C, Houseman G. 1988. The Mechanics of the Tibetan Plateau. London: Philosophical Transactions of the Royal Society.
Fu B H, Shi P L, Zhang Z W. 2008. Spatial characteristics of the surface rupture produced by the MS8.0 Wenchuan earthquake using high-resolution remote sensing imagery. Acta Geologica Sinica, 82(12): 1679-1687.
Harrison T M, Copeland P, Kidd W S F, et al. 1992. Raising Tibet. Science, 255(5052): 1663-1670. DOI:10.1126/science.255.5052.1663
Hu X P, Cui X F, Ning J Y, et al. 2012. Preliminary study on tectonic deformation models in the Longmenshan region based on focal mechanism solutions of the Wenchuan earthquake sequence. Chinese Journal of Geophysics, 55(8): 2561-2574. DOI:10.6038/j.issn.0001-5733.2012.08.008
Lai X A, Huang L R, Xu J S. 2004. Present-Day Crustal Movement in China Constrained. Beijing: Seismological Press.
Liu C, Shi Y L, Zhu B J, et al. 2014. Crustal rheology control on the mechanism of the earthquake generation at the Longmen Shan fault. Chinese Journal of Geophysics, 57(2): 404-418. DOI:10.6038/cjg20140207
Ma N J, Li J, Zhao Z Q. 2015. Distribution of the deviatoric stress field and plastic zone in circular roadway surrounding rock. Journal of China University of Mining & Technology, 44(2): 206-213.
Molnar P, Tapponnier. 1975. Cenozoic tenozoic tectonics of Asia: Effects of a continental collision. Science, 189(4201): 419-426. DOI:10.1126/science.189.4201.419
Qin X H, Tan C X, Sun J Z, et al. 2012. Experimental study of relation between in-situ crustal stress and rock elastic modulus. Rock and Soil Mechanics, 33(6): 1689-1695.
Royden L H, Burchfiel B C, King R W, et al. 1997. Surface deformation and lower crustal flow in eastern Tibet. Science, 276(5313): 788-790. DOI:10.1126/science.276.5313.788
Sugisaki R. 1976. Chemical characteristics of volcanic rocks: relation to plate movements. Lithos, 9(1): 17-30. DOI:10.1016/0024-4937(76)90053-0
Teng J W, Pi J L, Yang H, et al. 2014. Wenchuan-Yingxiu MS8.0 earthquake seismogenic faults and deep dynamic response. Chinese Journal of Geophysics, 57(2): 392-403. DOI:10.6038/cjg20140206
Van Der VooR, Spakman W, Bijwaard H. 1997. Tethyan subducted slabs under India. Earth and Planetary Science Letters, 171(1): 7-20.
Wan T F, Ren Z H. 1999. Research on the intraplate deformation velocity of China in meso-Cenozoic. Geoscience, 13(1): 83-92.
Wang L J, Cui J W, Zhou C J, et al. 2009. Numerical modeling for Wenchuan earthquake mechanism. Journal of Geomechanics, 12(5): 105-113.
Xie F R, Cui X F, Zhao J T. 2003. Analysis of global tectonic stress field. Earth Science Frontiers, 10(S1): 22-30.
Xiong H, Fu X M, Wang C Y, et al. 2015. Experimental study of sandstone under different confining pressure deformation characteristics. China Measurement & Test, 41(3): 113-116, 120.
Xu J R, Zhao Z X, Yuzo I. 2008. Regional characteristics of crustal stress field and tectonic motions in and around Chinese mainland. Chinese Journal of Geophysics, 51(3): 770-781.
Yan Z K, Li Y, Zhao G H, et al. 2014. The relationship between Lushan earthquake and Wenchuan earthquake by segmentation of geology and geomorphology of Longmen Shan. Chinese Journal of Nature, 36(1): 51-58.
Yang M J, Yu X B, Li J H. 2012. Comparative analysis for slope stability by ANSYS and FLAC. Science Technology and Engineering, 12(24): 6241-6244, 6251.
Ye Z R, Wang J. 2004. Dynamics of present-day crustal movement in the China mainland. Chinese Journal of Geophysics, 47(3): 456-461.
Yi G X, Wen X Z, Xin H, et al. 2013. Stress state and major-earthquake risk on the southern segment of the Longmen Shan fault zone. Chinese Journal of Geophysics, 56(4): 1112-1120. DOI:10.6038/cjg20130407
Zeng R S, Zhu J S, Zhou B, et al. 1992. Seismicity and focal mechanism in Tibetan Plateau and its implications to lithospheric flow. Acta Seismologica Sinica, 14(S1): 523-533.
Zhao Z Q, Ma N J, Jia H S, et al. 2013. Partitioning characteristics of gas channel of coal-rock mass in mining space and gas orientation method. International Journal of Mining Science and Technology, 23(6): 873-877. DOI:10.1016/j.ijmst.2013.10.014
Zhu S B, Zhang P Z. 2009. A study on the dynamical mechanisms of the Wenchuan MS8.0 earthquake, 2008. Chinese Journal of Geophysics, 52(2): 418-527.
蔡美峰, 王双红. 1997. 地应力状态与围岩性质的关系研究. 中国矿业, 6(6): 38–41.
陈国光, 计凤桔, 周荣军, 等. 2007. 龙门山断裂带晚第四纪活动性分段的初步研究. 地震地质, 29(3): 657–673.
邓起东, 陈社发, 赵小麟. 1994. 龙门山及其邻区的构造和地震活动及动力学. 地震地质, 16(4): 389–403.
付碧宏, 时丕龙, 张之武. 2008. 四川汶川MS8.0大地震地表破裂带的遥感影像解析. 地质学报, 82(12): 1679–1687. DOI:10.3321/j.issn:0001-5717.2008.12.005
胡幸平, 崔效锋, 宁杰远, 等. 2012. 基于汶川地震序列震源机制解对龙门山地区构造变形模式的初步探讨. 地球物理学报, 55(8): 2561–2574. DOI:10.6038/j.issn.0001-5733.2012.08.008
赖锡安, 黄立人, 徐菊生. 2004. 中国大陆现今地壳运动. 北京: 地震出版社.
柳畅, 石耀霖, 朱伯靖, 等. 2014. 地壳流变结构控制作用下的龙门山断裂带地震发生机理. 地球物理学报, 57(2): 404–418. DOI:10.6038/cjg20140207
马念杰, 李季, 赵志强. 2015. 圆形巷道围岩偏应力场及塑性区分布规律研究. 中国矿业大学学报, 44(2): 206–213.
秦向辉, 谭成轩, 孙进忠, 等. 2012. 地应力与岩石弹性模量关系试验研究. 岩土力学, 33(6): 1689–1695.
沈明荣, 陈建峰. 2006. 岩体力学. 上海: 同济大学出版社.
苏生瑞, 黄润秋, 王士天. 2002. 断裂构造对地应力场的影响及其工程应用. 北京: 科学出版社.
滕吉文, 皮娇龙, 杨辉, 等. 2014. 汶川—映秀MS8.0地震的发震断裂带和形成的深层动力学响应. 地球物理学报, 57(2): 392–403. DOI:10.6038/cjg20140206
万天丰. 1993. 中国东部中、新生代板内变形构造应力场及其应用. 北京: 地质出版社.
万天丰, 任之鹤. 1999. 中国中、新生代板内变形速度研究. 现代地质, 13(1): 83–92.
王连捷, 崔军文, 周春景, 等. 2009. 汶川5.12地震发震机理的数值模拟. 地质力学学报, 12(5): 105–113.
王涛, 韩煊, 赵先宇, 等. 2015. FlAC3D数值模拟方法及工程应用——深入剖析FLAC3D5.0. 北京: 中国建筑工业出版社.
谢富仁, 崔效锋, 赵建涛. 2003. 全球应力场与构造分析. 地学前缘, 10(S1): 22–30.
熊魂, 付小敏, 王从颜, 等. 2015. 砂岩在不同围压条件下变形特征的试验研究. 中国测试, 41(3): 113–116, 120. DOI:10.11857/j.issn.1674-5124.2015.03.026
徐纪人, 赵志新, 石川有三. 2008. 中国大陆地壳应力场与构造运动区域特征研究. 地球物理学报, 51(3): 770–781.
徐芝纶. 1980. 弹性力学简明教程. 北京: 高等教育出版社.
颜照坤, 李勇, 赵国华, 等. 2014. 从龙门山地质地貌分段性探讨芦山地震与汶川地震的关系. 自然杂志, 36(1): 51–58.
杨铭键, 余贤斌, 黎剑华. 2012. 基于ANSYS与FLAC的边坡稳定性对比分析. 科学技术与工程, 12(24): 6241–6244, 6251. DOI:10.3969/j.issn.1671-1815.2012.24.058
叶正仁, 王建. 2004. 中国大陆现今地壳运动的动力学机制. 地球物理学报, 47(3): 456–461.
易桂喜, 闻学泽, 辛华, 等. 2013. 龙门山断裂带南段应力状态与强震危险性研究. 地球物理学报, 56(4): 1112–1120. DOI:10.6038/cjg20130407
曾融生, 朱介寿, 周兵, 等. 1992. 青藏高原及其东部邻区的三维地震波速度结构与大陆碰撞模型. 地震学报, 14(S1): 523–533.
张培震. 2008. 青藏高原东缘川西地区的现今构造变形、应变分配与深部动力过程. 中国科学D辑:地球科学, 38(9): 1041–1056.
朱守彪, 张培震. 2009. 2008年汶川MS8.0地震发生过程的动力学机制研究. 地球物理学报, 52(2): 418–427.