应力增量触发断层岩体能量释放模拟与地震成因探讨——以龙门山断

地    震    学    报

ACTA  SEISMOLOGICA  SINICA

Vol. 41，No. 4

Jul.，2019

师皓宇，马念杰，石建军，李楠，马骥. 2019. 应力增量触发断层岩体能量释放模拟与地震成因探讨—以龙门山断裂带为例. 地震学报，41（4）：502−511. doi：10.11939/jass.20180151.

Shi H Y，Ma N J，Shi J J，Li N，Ma J. 2019. Simulation on energy release of fault rock mass triggered by stress incrementand discussion on seismogenesis：Taking Longmenshan fault zone as an example. Acta Seismologica Sinica，41（4）：502−511. doi：10.11939/jass.20180151.

应力增量触发断层岩体能量释放模拟

与地震成因探讨

—以龙门山断裂带为例

*

师皓宇1，2）    马念杰2）    石建军1）    李    楠1），     马    骥2）

1） 中国北京 101601 华北科技学院安全工程学院

2） 中国北京 100083 中国矿业大学（北京）资源与安全工程学院

摘要    本文以龙门山断裂带为背景，基于岩体应变能基本理论，使用FLAC软件模拟地震能量源和能量释放形式，计算结果显示：在0.01 MPa水平应力增量作用下，龙门山断裂带及附近区域可释放的应变能约为3.24×1013 J；使得断层面之间发生滑移，克服断层面滑动摩擦所需消耗的能量约为2.10×1013 J；岩体在重力方向上产生位移，克服重力做功所消耗的能量约为1.14×1013 J。由此可推断：在一定区域内，应力触发释放能量值与克服断层面滑动摩擦和克服重力做功所消耗的能量之和大致相当；应变能可能会在某一区域范围内集中释放，形成地震效应。本次应力增量触发断层周围岩体能量释放事件中，在映秀—北川断裂与灌县—安县断裂之间的局部区域集中释放的能量为7.67×1012 J，相当于一次MS5.39地震发生所释放的能量。关键词 应力增量    应变能    龙门山断裂带    数值模拟

doi：10.11939/jass.20180151        中图分类号：P315.72＋7        文献标识码：A        

Simulation on energy release of fault rock mass triggered bystress increment and discussion on seismogenesis：Taking

Longmenshan fault zone as an example

Shi Haoyu1，2）    Ma Nianjie2）    Shi Jianjun1）    Li Nan1），     Ma Ji2）

1） North China Institute of Science and Technology，Beijing 101601，China

2） College of Resources and Safety Engineering，China University of Mining and Technology

（Beijing），Beijing 100083，China

Abstract：Taking Longmenshan fault zone as the research background，basing on theory ofstrain energy of rock mass，this paper simulated energy source and energy release form of asimulated earthquake by using the software FLAC. Calculation results show that the total strainenergy that can be released from the Longmenshan fault zone and its vicinity is approximately

* 基金项目　国家自然科学基金（51674119）和高校基本科研业务费（3142018022）联合资助. 收稿日期　2018−12−27收到初稿，2019−03−06决定采用修改稿. 通信作者　e-mail：linan78@ncist.edu.cn

4 期

师皓宇等：应力增量触发断层岩体能量释放模拟与地震成因探讨

—以龙门山断裂带为例

503

3.24×1013 J under the action of horizontal stress increment 0.01 MPa. When stress triggeringcauses slip between fault planes，energy consumption for overcoming sliding friction on the sur-face is about 2.10×1013 J on fault plane. When the rock mass moves in the direction of gravity，energy consumption for overcoming gravity is 1.14×1013 J. Therefore it is deduced that，in acertain region，the energy released by stress triggering is approximately equal to the sum of con-sumption energy for overcoming the friction on fault surface and that for overcoming thegravity of regional rock mass. Strain energy may be released in a certain range，resulting inseismic effect. In the event of energy release of rock mass around fault triggered by stress incre-ment，as for the the local region between the Yingxiu-Beichuan fault and Guanxian-Anxianfault in this paper，energy released centrally is about 7.67×1012 J，equivalent to that for occur-rence of an earthquake with MS5.39.

Key words：stress increment；strain energy；Longmenshan fault zone；numerical simulation

引言

能量转化是材料物理过程的本质特征，岩石的破坏从根本上而言是能量变化过程中的一种失稳现象，是达到强度极限时岩石内部的弹性应变能释放的结果。处于不同应力状态的岩石对应于不同的能量状态，从弹性阶段到塑性阶段再到破坏，始终与外界进行着能量交换（温韬等，2016）。能量积聚过程实质上是应力变化过程，岩体能量释放实质上是地应力达到岩体承载极限，在一定触发条件下应力变化所致，且岩石变形或破坏过程中的能量耗散、释放与岩石破坏具有一定的内在联系（谢和平等，2005a，b；朱维申等，2001；许等，2011）。

我国地震频发，大部分M7.0以上的地震均发生在断裂带的边界附近（张国民等，2003；张培震等，2003），其中部分地震是受某些因素触发而引起地壳中的先存断层错动，或在某些特定区域内形成新生断层，从而释放能量，触发地震（Scholz，2002；徐锡伟等，2003；师皓宇，马念杰，2018；师皓宇等，2018）。断层间的相互运动、温度变化等必然造成岩体应力的变化（Meade，2007），进而对断层附近区域的地震起到触发作用（Velasco et al，2008；Kilb et al，2000），例如1992年美国加州兰德斯 （Landers） MW7.3地震、2011年日本东北（Tohoku-Oki）MW9.0 地震以及2017年我国九寨沟MW6.5地震均有较明显的同震应力触发作用（许才军等，2018）。1973—1976年在巴颜喀拉地块东边界的虎牙断裂带上发生的4次强震均存在显著的应力触发效应，这些地震改变了龙门山断裂带中南段及其周围区域的应力状态，对其后的2008年汶川MS8.0地震和2013年芦山MS7.0地震的发生具有较大影响（屈勇，朱航，2017）。

自二十世纪初，研究人员就地震的能量来源相继提出了诸多假说：Reid （1910）的弹性回跳假说认为地震能量是断层两侧岩体因地壳变形而产生和储存的弹性变形能；Bridgman（1945）的相变说认为地震能量来源于地下物质在临界温度和压力作用下，使得周围岩体应力状态发生改变而激发的地震波。然而地震、断层、能量之间尚未建立起符合逻辑的数学关系。数值模拟作为一种理想状态下的计算方法，虽然有一定的局限性，但能够再现某些特定状态的现象和结果（赵由佳等，2018a）。陶玮等（2011）以区域应力积累为基础，模拟了龙门山断裂带能量释放、断层位错等过程，本文拟采用FLAC软件对龙门山断裂带进行模拟计算，研究应力增量事件对深部岩体应变能量的巨大影响，以期获得断层岩体释放的应变能与断层面滑动摩擦和克服重力做功之间的关系，试图阐明应力增量对地震的触发机制。 

504

地        震        学        报41 卷

1 岩体应变能计算方法

1.1 岩体应变能释放机理

岩体空间任意一点（x，y，z）受到三向主应力σ1，σ2，σ3作用，如图1a所示；当对其施加一定的应力增量时，整个岩体的空间应力必然重新分布，其主应力状态如图1b所示；当应力重新达到平衡时，其主应力状态如图1c所示。该点在应力变化前后的应变能密度ua（x，y，z）和ub（x，y，z）（徐芝纶，1980）分别为

ua(x，y，z)＝ub(x，y，z)＝

122[σa1＋σ2a2＋σa3－2µ(σa1σa2＋σa1σa3＋σa2σa3)]，2E

122[σb1＋σ2b2＋σb3－2µ(σb1σb2＋σb1σb3＋σb2σb3)]，2E

（1）（2）

式中：σa1，σa2，σa3和σb1，σb2，σb3分别为应力变化前、后的最大、中间、最小主应力，E为岩体的弹性模量，μ为泊松比。则该位置处的能量密度差值ue（x，y，z）为

ue(x，y，z)＝ua(x，y，z)－ub(x，y，z).

（3）

当ue（x，y，z）＞0时，该点能量密度减小，表示该点释放能量；当ue(x，y，z)＜0时，该点能量密度增大，表示该点积聚能量。由式（1）可知，影响应变能密度的参数有弹性模量、泊松比、三向主应力。一定区域内所释放或积聚的能量值为

󰀂W＝ue(x，y，z)dV.

Ω

（4）

2

σ3

σ2

σa

3

σa

σb

3

σb

2

σ1

σ2

σ1

σa

1

σa

σa

1

σb

1

σb

σb

1

22

（a）σ3

（b）

σa

3

（c）σb

3

图 1    三向受力状态示意图

（a）原始应力状态；（b）增加Δσx后的应力状态;（c）平衡后的应力状态

Fig. 1    Schematic diagram of three-dimensional stress state

（a） Original stress state；（b） The stress state after adding Δσx；（c） The stress state after being balanced

由式（4）可知，释放的能量值与体积正相关。根据弹性力学可知，式（1）和（2）中的应变能密度u的单位为Pa，则式（4）中能量W的单位为J。由于实际的地壳中存在构造各异的复杂岩体，任意两点的能量密度都不尽相同，因此难以采用式（3）和（4）计算特定区域的能量释放或积聚值。而数值模拟则可将一定区域的地壳岩体划分为若干单元，当计算单元数量无限多、单元体积无限小时，计算结果即逼近于真值，从而获取相对真实的结果。1.2 龙门山断裂带数值模型

龙门山地处青藏高原东缘与四川盆地的交接部位（付碧宏等，2008；颜照坤等，2014），如图2所示。本文数值模拟中，选取垂直于龙门山断裂带走向的AA′剖面，该剖面的走向长 度为160 km，包含马尔康至四川盆地的部分区域，横跨龙门山的汶川—茂县断裂（F1）、映秀—北川断裂（F2）和灌县—安县断裂（F3）。模拟所用模型高40 km，厚1 km，汶川地震的震源

4 期

师皓宇等：应力增量触发断层岩体能量释放模拟与地震成因探讨

—以龙门山断裂带为例

103°E

F5

F2

104

F4

105°

广元

505

106°

大致位于深14 km的F2断层附近，地震频发区域基本包含于本文的选取范围。数值计算建模如图3所示，模型的上部为自由边界，模型的下部简化为z向位移边界；模型右侧为x向位移边界。龙门山断裂带呈铲型，近地表断层倾角由近36°逐渐增大，在西南端接近83°，并在地表附近出现滑动峰（Wan et al，2017）。对整个模型范围内的岩体施加重力应力场，使模型具备初始应力环境。模型左侧施加水平移动速度，且该速度随深度的增加而增大，上部加载的正向水平速度V1为2 mm/ a，下部加载的正向水平速度V2为4 mm/a （朱守彪，

地形的起伏简化为折线（邓起东等，1994）；

门 山A汶川

带茂县

32°N

江油绵阳

断F3A′F1

汶川MS8.0龙 裂都江堰

成都

F6

F1: 汶川—茂县断裂F2: 映秀—北川断裂

F3: 灌县—安县断裂31°F4: 青川断裂F5: 岷江断裂F6: 龙泉山断裂

地震断裂活动断层

图 2    龙门山及附近区域断层分布图Fig. 2    Distribution map of faults in

Longmenshan and its vicinity

张培震，2009；赖锡安等，2000；王连捷等，2009），位移量按加载步逐步施加。模型的物理力学参数详见表1 （沈明荣，陈建峰，2006；熊魂等，2015）。

V1=2 mm/a

40模型高度/km3020100

40

汶川－茂县断裂 F1映80

秀－北川断裂 F2县－县安断 F3裂灌V2=4 mm/a

120160

水平距离/km

图 3    数值计算模型图

Fig. 3    The numerical calculation model

表 1    模型岩体物理力学参数

Table 1    Physico-mechanical parameters for rock mass of the numerical model

弹性模量/GPa地表40

底部106

抗拉强度/MPa12

内聚力/MPa16

摩擦角/°35

泊松比0.286

密度/（103 kg·m−3）

2 650

重力加速度/（m·s−2）9.8

法向刚度/GPa

1

断层面

切向刚度/GPa摩擦角/°

0.5

10

虽然板块运动的起因与动力不得而知，但板块运动是客观存在的。持续的板块运动必然会使板块内部产生较大的应力变化，形成偏高水平的应力场，促使断裂带的形成和断层的持续滑移(师皓宇，马念杰，2018)。库仑破裂应力变化的典型值一般为0.1—1 MPa （Stein，1997；Harris，1998），宁夏海原MS8.5大地震引起甘肃古浪地震的断层面在滑动方向上产生0.01 MPa的静态库仑应力变化，并触发了甘肃古浪MS8.0大地震(傅征祥，刘桂萍，1999；万永革等，2002)，因此0.01 MPa触发应力在地震研究中是一个重要数值。本文将计算区域板块运动中的某一时刻（30万年），在整个模型上施加0.01 MPa水平应力增量的情况下，模拟一次应力增量事件对龙门山断裂带周围岩体能量释放、岩体形变和断层滑移的影响。

506

地        震        学        报41 卷

1.3 数据处理方法

模拟计算结果的波动性必然带来一定的误差，计算模型需要较高的计算精度。在施加水平应力增量前后，设定不平衡率为10−9，计算结果则逼近其真值。计算思路如下：

1） 建立数值计算模型，使计算模型的初始状态达到稳态；

2） 在水平方向给予一定的应力增量∆σx，使计算模型不平衡率达到10−9；3） 导出计算模型中的主应力、弹性模量、泊松比等数据。

对于计算模型的单元体而言，其应力分布均匀，根据式（4），则第i个单元体的能量释放值Wei为

Wei＝uei·Vi，

体个数为n时，则该区域释放的应变能为

（5）

式中，uei为第i个单元体的能量释放密度值，Vi为第i个单元体的体积。一定区域内的单元

n∑i＝1

We＝Wei.

（6）

2 岩体应变能释放密度分布

从计算模型中导出所有单元体的主应力值、体积、弹性模量、泊松比等参数，代入式（1）和（2），即可得到每一单元体触发前后的应变能、能量密度及能量释放密度。深部岩体的应变能巨大，当对其施加一定的应力增量时，其内部应力重新分布，会导致局部应变能减小，减小的应变能必将被释放或转移，结果显示触发前总应变能可达2.61×1019 J。将模型中所有单元体的能量释放密度值通过Surfer处理，得到如图4所示的能量释放密度分布图。由图4可见，当触发应力为0.01 MPa时，能量释放区域主要集中在F2与F3之间的某一区域，将范

3能量释放·密m-)度/(J模型高度/km围x∈（100，130）且y∈（14，24）所属范围命名为区域A。计算模型因水平应力增加而导致计

40302010020406080100水平距k离m/120140能量释放区 A10048－4F1F2F3－56－108160－160 

图 4    能量释放密度图（断裂名称同图3）

Fig. 4    Distribution map of energy release density （The faults are the same as those in Fig.3）

13951－3－7最大主应力变化值/kPa40模型高度/km302010020406080水平距离/km10012014016017 

图 5    应力触发前后最大主应力变化分布图

Fig. 5    Distribution map of maximum principal stress variation before and after stress triggering

4 期

师皓宇等：应力增量触发断层岩体能量释放模拟与地震成因探讨

—以龙门山断裂带为例

507

算模型应变能总量增加，同时释放出的巨量应变能约为 3.24×1013 J。图4中区域A为能量释放密度值较高区域，该区域释放的能量约为7.67×1012 J，约占释放总能量的23.6%。

应力触发前后的最大主应力变化如图5所示，对比图4与图5可知，应变能释放区与主应力变化区的位置基本对应。假如把此次能量释放事件作为一次地震，按照里克特级数的震级与能量计算公式MS＝(lgW－4.8)/1.5计算，其释放的能量相当于MS5.39地震发生所释放的能量，而区域A则可能为震源位置。因此微小应力触发一次地震事件理论上是可能的。

从应力增加至应力平衡，即不平衡率重期T，区域A释放的能量为7.67×10 J，将整个周期平均分为10段，计算每一段释放的能量总值，结果如图6所示，可见：第一段内，释放能量为5.61×1012 J，约占总能量的73.2%；第二段内释放能量累计值为6.75×1012 J，约占总能量的88.1%。因此，在整个能量释放周期内，能量释放主要集中在前两段，随后能量释放速度逐步放缓。

12

876能量/(1012 J)5432100

2

4

6

能量释放周期分段序号

释放能量累加值新达到10−9时，作为一个完整的能量释放周

不同周期段能量释放值810

图 6    能量释放全周期变化曲线图

Fig. 6    Variation of energy release curve in a full cycle

3 断层面摩擦做功消耗能量

断层间的滑移是地震发生所造成的普遍现象(赵由佳等，2018b)，大部分地震滑动发生在两个陡倾断层面的浅层地壳中(Wang et al，2011)。断层滑移量和滑移面积是计算地震矩的重要指标，对于龙门山断裂带南段发生的芦山地震，其震源处最大的滑移量为1.5 m （刘成利等，2013），而汶川地震的同震垂直位移达12 m (王卫民等，2008)。本模型的计算结果同样存在断层之间的滑移现象，地震发生后，灌县—安县断裂（F3）的上下盘之间的滑移量明显高于其它两条断层。提取模型所有节点的位移数据可知，岩体变形主要集中在映秀—北川断裂（F2）与灌县—安县断裂（F3）之间，能量释放区域的上方，因此地震发生后，断层上盘区域往往是重灾区，这在许多逆冲型地震中均有体现。断层发生滑移必然耗散能量，单位面积上发生滑动时，克服摩擦做功的计算相对简单，但由于整个计算模型不同位置的位移量各不相同，因此后期数据的处理较为困难。

本文计算模型中，汶川—茂县断裂（F1）的断层面有194个单元面，映秀—北川断裂（F2）的断层面有214个单元面，灌县—安县断裂（F3）的断层面有214个单元面，第i个断层单元面上克服摩擦做功的计算公式为

Wfi＝Ni·di＝σNi·Si·tanφi·di

元的摩擦角。则整个断层面上克服摩擦所做的功为

n∑Wf＝Wfi.

i＝1

（7）

式中：Wfi为克服摩擦所做的功，σNi为法向力，di为位移量，Si为接触单元面积，φi为接触单

（8）

通过Fish语言导出计算模型中接触面的正应力、位移、面积，并将其代入式（7）和式（8），

508

地        震        学        报41 卷

经计算可得，本次应力触发事件所导致的3条断层克服摩擦所做的功总值为2.10×1013 J，其中克服F1，F2和F3断层摩擦所做的功分别为1.65×1012 J，5.65×1012 J和1.37×1013 J。

4 克服重力做功耗损能量

因应力重新分布而导致地壳岩体发生形变，处于重力场中的岩体一定会在竖直方向上产生正向或负向位移，因而必定在重力作用下做功或克服重力做功。导出计算模型中节点的位移值，并将与节点对应的位移值导入Surfer软件处理得到位移，如图7所示，可见：在F2与F3断层之间的部分岩体具有明显的垂直向上的位移。但FLAC计算软件所导出的单元块体积是基于单元体（zone），而导出的位移是基于节点（grid point），因此两种数据不能完全对应，需分区计算，即沿水平方向分为8段，沿竖直方向分为4段，分别导出所有单元块体的体积和模型节点的位移平均值，对计算模型进行分区域计算，得到克服重力所做的功的总值。对于单元体而言，克服重力所做的功Wgi为

Wgi＝mi·g·di＝Vi·ρ·g·di，

式中mi为单元体重量，g为重力加速度，di为单元体位移量，Vi为单元体体积。

（9）

克服重力做功的能量密度分布如图8所示。克服重力做功主要集中在F2与F3断层之1.2×1013 J/m3。

403020100

20406080100水平距k离m/

120

140

F1

F2

F3

垂直位移量/mm间，即岩体克服重力做功的密度与岩体竖向位移有关，克服重力做功的最大能量密度为

模型高度/km4.03.02.01.00

160

－1.0

图 7    岩体垂直位移分布图

Fig. 7    Distribution of vertical displacement of rock mass

40模型高度/km3020100

20

40

60

80

水平距离/km

100

120

140

160

F1

F2

F3

12.09.57.55.53.51.5－0.5－2.5

克服重力做功/(104 J·m-3) 

图 8    克服重力做功能量密度图

Fig. 8    Energy density map for overcoming gravity

5 讨论与结论

因应力变化而促使岩体释放的能量必然以岩体动能、重力势能和克服断层面摩擦做功等形式向外释放（Dahlen，1977），或破坏岩石之间的晶格结构产生热量，或引起岩体的震动并以波的形式传播出去(Gudmundsson，2014)。从本模型的计算结果来看，模型释放的能量We为3.24×1013 J，克服断层面摩擦做功消耗的能量Wf为2.10×1013 J，因岩体形变而克服重

4 期

师皓宇等：应力增量触发断层岩体能量释放模拟与地震成因探讨

—以龙门山断裂带为例

509

力做功耗损的能量Wg为1.14×1013 J，即We≈Wf＋Wg，这表明微小的水平应力增量将导致断裂带岩体释放巨量的应变能，这些能量释放主要通过两种形式，一是用于克服摩擦做功，使断层面发生滑移；二是克服重力做功，使地壳岩体发生形变。

将整个周期的能量释放全过程划分为10段，从模拟结果来看，第一段内能量释放值约占总能量的73.2%，表明能量更容易在前期很短的时间内集中释放，具备地震发生的基本条件，符合主余震型地震的特征。以本次事件计算为例，0.01 MPa水平应力可触发MS5.39地震。因本文所模拟的断层走向长度仅1 km，而汶川地震断裂长达300 km左右，释放的能量值可达2.28×1015 J，相当于MS7.04地震。地震能量的释放随着触发应力的增大而增大，经测算，当触发应力达到1 MPa时，释放的能量值则可造成大概MS8.38地震的发生。

本文研究结果能够说明能量释放的空间位置和能量量级，且该能量释放具有完整周期，但由于模拟计算时间属于相对时间，尚不能与实际时间准确对应，因此不能定性为地震能。但模型应变能的减小值与断层面滑动摩擦做功和克服重力做功之和相吻合，这一结果证明因应力变化而造成地壳岩体释放的应变能可能是断层滑移、地壳形变甚至是岩体弹塑性状态改变、地震等现象发生的重要能量源，论证了微小应力触发地震的可能性和科学性，解释了断层附近的地震能量来源与释放形式。

参 考 文 献

邓起东，陈社发，赵小麟. 1994. 龙门山及其邻区的构造和地震活动及动力学[J]. 地震地质，16（4）：3–403.

Deng Q D，Chen S F，Zhao X L. 1994. Tectonics，seismicity and dynamics of Longmenshan mountains and its adjacent

regions[J]. Seismology and Geology，16（4）：3–403 （in Chinese）.

付碧宏，时丕龙，张之武. 2008. 四川汶川MS8.0大地震地表破裂带的遥感影像解析[J]. 地质学报，82（12）：1679–1687.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[J]. Acta Geologica Sinica，82（12）：1679–1687 （in Chinese）.

傅征样，刘桂萍. 1999. 海原大地震可能触发古浪大地震的力学机制[C]//中国地震学会成立20周年纪念文集. 北京：中国

地震学会：234−243.

Fu Z X，Liu G P. 1999. The mechanism of great Gulang earthquake triggered probably by the great Haiyuan earthquake[C]//Pro-ceedings of the 20th Anniversary Collection of the Chinese Seismological Society. Beijing：Seismological Press：234−243 （inChinese）.

赖锡安，许菊生，卓力格图，刘经南，施闯，姜卫平. 2000. 中国主要构造块体现今运动的基本特征[J]. 中国地震，

16（3）：213–222.

Lai X A，Xu J S，Zhuoli G T，Liu J N，Shi C，Jiang W P. 2000. The fundamental characteristics of principal tectonic blocks

present movement in Chinese mainland[J]. Earthquake Research in China，16（3）：213–222 （in Chinese）.

刘成利，郑勇，葛粲，熊熊，许厚泽. 2013. 2013年芦山7.0级地震的动态破裂过程[J]. 中国科学：地球科学，43（6）：

1020–1026.

Liu C L，Zheng Y，Ge C，Xiong X，Xu H Z. 2013. Rupture process of the MS7.0 Lushan earthquake，2013[J]. Science China

Earth Sciences，56（7）：1187–1192.

屈勇，朱航. 2017. 巴颜喀拉块体东—南边界强震序列库仑应力触发过程[J]. 地震研究，40（2）：216–225.

Qu Y，Zhu H. 2017. Coulomb stress triggering process of major earthquake sequence in the eastern and southern boundaries of the

Bayan Har block[J]. Journal of Seismological Research，40（2）：216–225 （in Chinese）.沈明荣，陈建峰. 2006. 岩体力学[M]. 上海：同济大学出版社：52−54.

Shen M R，Chen J F. 2006. Rock Mass Mechanics[M]. Shanghai：Tongji University Press：52−54 （in Chinese）.

师皓宇，马念杰. 2018. 龙门山断裂带及附近区域地貌形成与地应力演化机制研究[J]. 地震学报，40（3）：332–340.Shi H Y，Ma N J. 2018. Geomorphic formation and crustal stress evolution mechanism in the Longmenshan fault zone and its adja-

510

地        震        学        报41 卷

cent regions[J]. Acta Seismologica Sinica，40（3）：332–340 （in Chinese）.

师皓宇，马念杰，马骥. 2018. 龙门山断裂带形成过程及其地应力状态模拟[J]. 地球物理学报，61（5）：1817–1823.Shi H Y，Ma N J，Ma J. 2018. Numerical simulation for the formation process of the Longmenshan fault zone and its crustal

stress state[J]. Chinese Journal of Geophysics，61（5）：1817–1823 （in Chinese）.

陶玮，胡才博，万永革，沈正康，王康. 2011. 铲形逆冲断层地震破裂动力学模型及其在汶川地震研究中的启示[J]. 地球物

理学报，54（5）：1260–1269.

Tao W，Hu C B，Wan Y G，Shen Z K，Wang K. 2011. Dynamic modeling of thrust earthquake on listric fault and its inference to

study of Wenchuan earthquake[J]. Chinese Journal of Geophysics，54（5）：1260–1269 （in Chinese）.万永革，吴忠良，周公威，黄静，秦立新. 2002. 地震应力触发研究[J]. 地震学报，24（5）：533–551.

Wan Y G，Wu Z L，Zhou G W，Huang J，Qin L X. 2002. Research on seismic stress triggering[J]. Acta Seismologica Sinica，

24（5）：533–551 （in Chinese）.

王连捷，崔军文，周春景，孙东生，王薇，唐，钱华山. 2009. 汶川5·12地震发震机理的数值模拟[J]. 地质力学学报，

15（2）：105–113.

Wang L J，Cui J W，Zhou C J，Sun D S，Wang W，Tang Z M，Qian H S. 2009. Numerical modeling for Wenchuan earthquake

mechanism[J]. Journal of Geomechanics，15（2）：105–113 （in Chinese）.

王卫民，赵连锋，李娟，姚振兴. 2008. 四川汶川8.0级地震震源过程[J]. 地球物理学报，51（5）：1403–1410.

Wang W M，Zhao L F，Li J，Yao Z X. 2008. Rupture process of the MS8.0 Wenchuan earthquake of Sichuan，China[J]. Chinese

Journal of Geophysics，51（5）：1403–1410 （in Chinese）.

温韬，唐辉明，刘佑荣，王康，杨呈刚. 2016. 不同围压下板岩三轴压缩过程能量及损伤分析[J]. 煤田地质与勘探，

44（3）：80–86.

Wen T，Tang H M，Liu Y R，Wang K，Yang C G. 2016. Energy and damage analysis of slate during triaxial compression under

different confining pressures[J]. Coal Geology & Exploration，44（3）：80–86 （in Chinese）.

谢和平，鞠杨，黎立云. 2005a. 基于能量耗散与释放原理的岩石强度与整体破坏准则[J]. 岩石力学与工程学报，24（17）：

3003–3010.

Xie H P，Ju Y，Li L Y. 2005a. Criteria for strength and structural failure of rocks based on energy dissipation and energy release

principles[J]. Chinese Journal of Rock Mechanics and Engineering，24（17）：3003–3010 （in Chinese）.

谢和平，彭瑞东，鞠杨，周宏伟. 2005b. 岩石破坏的能量分析初探[J]. 岩石力学与工程学报，24（15）：2603–2608.Xie H P，Peng R D，Ju Y，Zhou H W. 2005b. On energy analysis of rock failure[J]. Chinese Journal of Rock Mechanics and En-gineering，24（15）：2603–2608 （in Chinese）.

熊魂，付小敏，王从颜，宾婷婷，沈忠，黄兴建. 2015. 砂岩在不同围压条件下变形特征的试验研究[J]. 中国测试，43（3）：

113–116.

Xiong H，Fu X M，Wang C Y，Bin T T，Shen Z，Huang X J. 2015. Experimental study of sandstone under different confining

pressure deformation characteristics[J]. China Measurement & Testing，43（3）：113–116 （in Chinese）.

许才军，汪建军，熊维. 2018. 地震应力触发回顾与展望[J]. 武汉大学学报（信息科学版），43（12）：2085–2092.

Xu C J，Wang J J，Xiong W. 2018. Retrospection and perspective for earthquake stress triggering[J]. Geomatics and Information

Science of Wuhan University，43（12）：2085–2092 （in Chinese）.

许，牛双建，靖洪文，杨圣奇，王文龙. 2011. 砂岩加卸载条件下能耗特征试验研究[J]. 岩土力学，32（12）：

3611–3617.

Xu G A，Niu S J，Jing H W，Yang S Q，Wang W L. 2011. Experimental study of energy features of sandstone under loading and

unloading[J]. Rock and Soil Mechanics，32（12）：3611–3617 （in Chinese）.

徐锡伟，闻学泽，郑荣章，马文涛，宋方敏，于贵华. 2003. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科

学：D辑，33（增刊1）：151–162.

Xu X W，Wen X Z，Zheng R Z，Ma W T，Song F M，Yu G H. 2003. Pattern of latest tectonic motion and its dynamics for

active blocks in Sichuan-Yunnan region，China[J]. Science in China：Series D，46（S2）：210–226.徐芝纶. 1980. 弹性力学简明教程[M]. 北京：人民教育出版社：93−95.

Xu Z L. 1980. Brief Course of Elastic Mechanics Theory[M]. Beijing：People’s Education Press：93−95 （in Chinese）.颜照坤，李勇，赵国华，周荣军，李敬波，张威，郑立龙，李奋生，闫亮. 2014. 从龙门山地质地貌分段性探讨芦山地震与

4 期

师皓宇等：应力增量触发断层岩体能量释放模拟与地震成因探讨

—以龙门山断裂带为例

511

汶川地震的关系[J]. 自然杂志，36（1）：51–58.

Yan Z K，Li Y，Zhao G H，Zhou R J，Li J B，Zhang W，Zheng L L，Li F S，Yan L. 2014. The relationship between Lushan

earthquake and Wenchuan earthquake by segmentation of geology and geomorphology of Longmen Shan[J]. Chinese Journalof Nature，36（1）：51–58 （in Chinese）.

张国民，田勤俭，王辉. 2003. 可可西里—东昆仑活动构造带强震活动研究[J]. 地学前缘，10（1）：39–46.

Zhang G M，Tian Q J，Wang H. 2003. Strong earthquake activities in Kekexili-East Kunlun mountains active fault zone，Northw-est China[J]. Earth Science Frontiers，10（1）：39–46 （in Chinese）.

张培震，邓起东，张国民，马瑾，甘卫军，闵伟，毛凤英，王琪. 2003. 中国的强震活动与活动地块[J]. 中国科学：

D辑，33（增刊1）：12–20.

Zhang P Z，Deng Q D，Zhang G M，Ma J，Gan W J，Min W，Mao F Y，Wang Q. 2003. Active tectonic blocks and strong earth-quakes in the continent of China[J]. Science in China：Series D，46（S2）：13–24.

赵由佳，张国宏，单新建，尹昊，屈春燕. 2018a. 考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟[J]. 地球物

理学报，61（5）：1853–1862.

Zhao Y J，Zhang G H，Shan X J，Yin H，Qu C Y. 2018a. Numerical simulation of the strong ground motion of the 2008 Wen-chuan earthquake incorporated with topography and barrier rupture model[J]. Chinese Journal of Geophysics，61（5）：1853–1862 （in Chinese）.

赵由佳，张国宏，张迎峰，单新建，屈春燕. 2018b. 基于连续-离散单元法的汶川地震动力学二维自发破裂全周期模拟研

究[J]. 地震地质，40（1）：12–26.

Zhao Y J，Zhang G H，Zhang Y F，Shan X J，Qu C Y. 2018b. Two-dimensional whole cycle simulation of spontaneous rupture

of the 2008 Wenchuan earthquake using the continuous-discrete element method[J]. Seismology and Geology，40（1）：12–26（in Chinese）.

朱守彪，张培震. 2009. 2008年汶川MS8.0地震发生过程的动力学机制研究[J]. 地球物理学报，52（2）：418–427.

Zhu S B，Zhang P Z. 2009. A study on the dynamical mechanisms of the Wenchuan MS8.0 earthquake，2008[J]. Chinese Journal

of Geophysics，52（2）：418–427 （in Chinese）.

朱维申，李术才，程峰. 2001. 能量耗散模型在大型地下洞群施工顺序优化分析中的应用[J]. 岩土工程学报，23（3）：

333–336.

Zhu W S，Li S C，Cheng F. 2001. Application of energy dissipation model to optimization of construction order for large under-ground caverns[J]. Chinese Journal of Geotechnical Engineering，23（3）：333–336 （in Chinese）.Bridgman P W. 1945. Polymorphic transitions and geological phenomena[J]. Am J Sci，243（2）：90–97.Dahlen F A. 1977. The balance of energy in earthquake faulting[J]. Geophys J R astr Soc，48（2）：239–261.Gudmundsson A. 2014. Elastic energy release in great earthquakes and eruptions[J]. Front Earth Sci，2：10.

Harris R A. 1998. Introduction to special section：Stress triggers，stress shadows，and implications for seismic hazard[J]. J Geo-phys Res，103（B10）：24347–24358.

Kilb D，Gomberg J，Bodin P. 2000. Triggering of earthquake aftershocks by dynamic stresses[J]. Nature，408（6812）：570–574.Meade B J. 2007. Present-day kinematics at the India-Asia collision zone[J]. Geology，35（1）：81–84.

Reid H F. 1910. Mechanics of the Earthquake，the California Earthquake of April 18，1906[R]. Washington DC：Carnegie Insti-tution of Washington：192.

Scholz C H. 2002. The Mechanics of Earthquakes and Faulting[M]. 2nd ed. Cambridge：Cambridge University Press：504.Stein R S，Barka A A，Dieterich J H. 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress trig-gering[J]. Geophys J Int，128（3）：594–604.

Velasco A A，Hernandez S，Parsons T，Pankow K. 2008. Global ubiquity of dynamic earthquake triggering[J]. Nat Geosci，

1（6）：375–379.

Wan Y G，Shen Z K，Burgmann R，Sun J B，Wang M. 2017. Fault geometry and slip distribution of the 2008 MW7.9 Wen-chuan，China earthquake，inferred from GPS and InSAR measurements[J]. Geophys J Int，208（2）：748–766.

Wang Q，Qiao X J，Lan Q G，Freymueller J，Yang S M，Xu C J，Yang Y L，You X Z，Tan K，Chen G. 2011. Rupture of deep

faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan[J]. Nat Geosci，4（9）：634–0.