大地测量与地球动力学

龙门山断裂带映秀段地震破裂全周期形变数值模拟

吴若君 ¹^,² , 董培育 ¹^,²^,* , 赵斌 ¹^,²

大地测量与地球动力学 | 特色栏目: 地震大地测量学 2026,46(8): 949-956

收起

大地测量与地球动力学 |特色栏目: 地震大地测量学 2026 , 46 (8) : 949 -956

龙门山断裂带映秀段地震破裂全周期形变数值模拟

全屏

吴若君¹^,², 董培育¹^,²^,*, 赵斌¹^,²

作者信息

1 中国地震局地震研究所, 武汉, 430071

2 中国地震局地震大地测量重点实验室, 武汉, 430071

通讯作者:

董培育, 博士, 副研究员, 主要研究方向为地震大地测量学与地球动力学数值模拟, E-mail: dongpeiyu97@163.com

作者简介:

吴若君，硕士生，主要研究方向为大地测量学与地球动力学数值模拟，E-mail: rjwu0719@163.com

Numerical Simulation of Full-Cycle Earthquake Rupture Deformation in the Yingxiu Section of the Longmenshan Fault Zone

Ruojun WU¹^,², Peiyu DONG¹^,²^,*, Bin ZHAO¹^,²

Affiliations

1 Institute of Seismology, CEA, Wuhan 430071, China

2 Key Laboratory of Earthquake Geodesy, CEA, Wuhan 430071, China

出版时间: 2026-08-15 doi: 10.14075/j.jgg.2025.12.407

文章导航

摘要

收起

以龙门山断裂带映秀段为研究区，构建二维纵剖面黏弹性有限元模型。模型充分考虑研究区域介质的纵向分层结构与横向非均质性，基于滑移弱化摩擦准则，对逆冲型地震的周期性自发破裂全过程进行数值模拟，完整覆盖震间、同震及震后应力调整3个关键阶段。结果表明，在模拟震间速度场与汶川地震震前GPS观测数据达到较好拟合效果的前提下，静摩擦系数与动摩擦系数的差值对同震地表位移量具有显著控制作用，其中拟合效果最优的静摩擦系数与动摩擦系数分别为0.35和0.23；模拟得到的单次地震断层面同震破裂最大滑动量约为7 m。由于本研究震后形变模拟仅考虑中下地壳黏弹性松弛效应，未纳入震后余滑效应的影响，导致近场震后形变模拟值低于实际观测值，而远场模拟结果与观测数据拟合度较高。在本研究设定的计算参数与假设条件下，模拟得到的汶川地震复发周期约为700~770 a，与基于龙门山断裂带映秀段深部滑脱层闭锁模型的估计结果基本一致。

关键词

滑移弱化摩擦准则 / 地震周期变形 / 龙门山断裂带 / 汶川地震 / 数值模拟

Abstract

收起

Taking the Yingxiu section of the Longmenshan fault zone as the study area, a 2D vertical-profile viscoelastic finite element model is constructed. The model fully considers the vertical stratification and lateral heterogeneity of the medium in the study area. Based on the slip-weakening friction law, numerical simulations are carried out for the entire periodic spontaneous rupture process of thrust earthquakes, covering three key stages: interseismic, coseismic, and postseismic stress adjustments.The results show that, on the premise that the simulated interseismic velocity field is in good agreement with the pre-earthquake GPS observations of the Wenchuan earthquake, the difference between the static and dynamic friction coefficients exerts a significant control on the coseismic surface displacement. The optimally fitted static and dynamic friction coefficients are 0.35 and 0.23, respectively. The simulated maximum coseismic slip on the fault plane during a single earthquake is approximately 7 m. Since the postseismic deformation simulation in this study only considers the viscoelastic relaxation effect of the middle and lower crust without incorporating postseismic afterslip, the simulated near-field postseismic deformation is lower than the observed values, while the far-field simulation results agree well with the observations. Under the computational parameters and assumptions adopted in this study, the recurrence interval of the Wenchuan earthquake is simulated to be approximately 700 to 770 years, which is generally consistent with estimates based on the deep detachment layer locking model for the Yingxiu section of the Longmenshan fault zone.

Key words

slip-weakening friction law / seismic cycle deformation / Longmenshan fault zone / Wenchuan earthquake / numerical simulation

引用本文

吴若君, 董培育, 赵斌. 龙门山断裂带映秀段地震破裂全周期形变数值模拟. 大地测量与地球动力学, 2026 , 46 (8) : 949 -956 . DOI: 10.14075/j.jgg.2025.12.407

Ruojun WU, Peiyu DONG, Bin ZHAO. Numerical Simulation of Full-Cycle Earthquake Rupture Deformation in the Yingxiu Section of the Longmenshan Fault Zone[J]. Journal of Geodesy and Geodynamics, 2026 , 46 (8) : 949 -956 . DOI: 10.14075/j.jgg.2025.12.407

正文

收起

地表形变是深部孕震过程及地震活动在地表的直接响应，形变监测能够直观反映地震不同时期的形变特征。然而，大地震的复发周期通常处于百年至千年尺度，而利用现代大地测量学手段(GPS、InSAR)开展的地表形变监测仅持续约30 a，目前尚未覆盖一个完整的地震周期。因此，开展地震周期性循环演化的数值模拟研究，成为了解地震全周期变形特征及地震发生规律的重要手段^[1]。Barbot等^[2]基于速率-状态相关摩擦准则，成功再现了圣安德烈斯断层Parkfield段类似1966年和2004年的M_W 6.0走滑型地震序列，模型预测的同震、震后变形与2004年地震的实际观测值高度吻合。随着地震周期动力学数值模拟研究的不断深入，研究重点逐步扩展至具有复杂构造特征的俯冲带系统，前人从断层本构参数、三维几何形态到热-力学耦合等多个角度细化物理模型，通过构建更贴近实际地质条件的地质结构与物理过程，显著提升了对地震全周期变形过程的解析能力^[3-5]。

2008年汶川地震是有地震历史记载以来，首次发生在大陆内部的高角度逆冲型强震。该地震在龙门山断裂带西南段以逆冲破裂为主，东北段则兼具右旋走滑分量，地表破裂总长度超过200 km^[6]。震前GPS形变监测结果表明，横跨龙门山断裂带的百千米尺度地壳缩短速率约为4 mm/a，表现出较低的地表变形速率特征^[7-10]，远低于绝大多数孕育8级以上逆冲型地震的大型俯冲带^[6]。然而，汶川地震同震破裂产生的最大滑动量却高达7~10 m^[11-12]。目前，不同学者基于不同数据资料和研究方法，对汶川地震复发周期的推测存在显著差异^[13-18](范围为600 a至万年)，且对控制该区域地震周期变形的动力学因素仍不明确。2008年汶川地震留下了丰富的大地测量观测资料，这些资料为揭示该区域的形变特征与动力学机制提供了重要的地表形变约束，对深化认识龙门山断裂带孕震与破裂行为具有重要意义。

本研究基于滑移弱化摩擦破裂准则，在龙门山断裂带映秀段开展逆冲型地震周期变形模拟研究，在明确该区域地震全周期变形特征的同时，重点探讨震间、同震和震后变形的主要控制因素，以期为理解陆内逆冲型地震的物理过程提供参考。

1 模拟计算

收起

1.1 模型建立及参数设置

2008年汶川地震在龙门山断裂带映秀段呈现以逆冲分量为主的破裂形式，本文在震中附近垂直于断层方向建立一个二维纵剖面黏弹性有限元模型(图 1)，模型长度为500 km，深度为100 km，共划分2 365个三角形单元。龙门山断裂带西北侧为川西高原，东南侧为四川盆地。参考Wang等^[12]反演的汶川地震同震破裂模型中的断层几何形态，将龙门山断裂带设置为铲形，15 km以上倾角为55°，15~22 km倾角为7°。上地壳采用理想弹性材料，根据Crust1.0地壳结构模型，设置川西高原一侧厚度为30 km，四川盆地一侧厚度为20 km。中下地壳和岩石圈上地幔为Maxwell体，黏滞系数的设置参考汶川地震震后形变约束下的龙门山地区流变结构^[19-20]，川西高原一侧的中下地壳和上地幔黏滞系数分别设为8.5×10¹⁸Pas和1.0×10¹⁹Pas，四川盆地一侧的中下地壳和上地幔分别设为5.0×10²⁰Pas和1.0×10²¹Pas。

假设研究区域地壳缩短速率稳定且均匀，模型左侧边界设置为均匀的挤压速率边界条件，且不随深度变化，具体数值将通过网格搜索确定；模型右侧边界约束法向(即水平方向)位移为0，切向(竖直方向)自由；模型底边界约束垂向位移为0，水平方向自由；模型上边界设置为自由表面。

1.2 摩擦本构关系及参数

本研究采用滑移弱化摩擦破裂准则^[21]：

(1)

$\tau=\left\{\begin{array}{l}{\left[\mu_s-\left(\mu_s-\mu_d\right) \frac{u}{d_0}\right] \sigma_n, u<d_0} \\\mu_d \sigma_n, u \geqslant d_0\end{array}\right.$

式中，μ_s、μ_d分别为最大静摩擦系数和滑动摩擦系数；u为滑移量(即断层两盘的相对错动量)；d₀为滑移弱化距离，代表静摩擦系数下降到动摩擦系数所产生的滑动量；σ_n为断层面上的有效正应力。震间期，断层处于闭锁状态，应力/应变持续积累，剪应力低于最大摩擦力μ_sσ_n，断层面不发生滑动。当应力超过最大摩擦力阈值，断层开始出现微弱滑移，摩擦力随着位错的增大而线性减小，一旦滑移量u超过滑移弱化距离d₀，静摩擦系数快速弱化为动摩擦系数，导致应力快速释放，断层表现出黏滑失稳破裂。

参考Zhang等^[22]模拟汶川地震同震破裂传播过程使用的断层摩擦系数和初始应力设置，并经过多次测算，最终得到与同震地表水平向位移观测值拟合最优的断层面上摩擦参数的分布情况(图 2)，其中在断层15 km深度以上，静摩擦系数和动摩擦系数分别为0.35和0.23。

2 模拟结果与分析

收起

计算使用有限元开源程序Pylith^[23]，模拟研究区域在长期均匀推挤作用下，由滑移弱化摩擦准则控制的逆冲型地震周期性自发破裂全过程，包括震间、同震和震后长期调整3个阶段，其中不考虑同震破裂传播阶段，仅给出同震破裂滑动量。

以汶川地震震前地表水平向GPS速度观测值为约束，采用网格搜索方法，确定边界挤压速率的最优值。考虑到百千米范围内区域汇聚速率约为4 mm/a^[7-10]，而本研究模型横跨龙门山断裂带500 km，左边界距离龙门山断裂带约300 km，因此设置边界挤压速率搜索下限为4.0 mm/a；根据汶川地震震前GPS速度场^[7-10](图 3)，设置边界挤压速率上限为8.0 mm/a，搜索步长为0.2 mm/a。对于每种边界条件下的模型均计算3万年，并提取震前10 a地表水平速度平均值作为震间速度，依据式(2)计算与GPS观测值之间的拟合残差RMS值。网格搜索确定的最佳边界挤压速率为6.9 mm/a, 在此条件下拟合残差最小, 为1.3 mm/a。

(2)

$\mathrm{RMS}=\sqrt{\frac{1}{n} \sum\limits_{i=1}^n\left(y_i-\hat{y}_i\right)^2}$

式中，y_i为第i个观测值；

$\hat{y}_i$

为第i个模拟值；n为样本数量。

该边界挤压速率结果与Thompson等^[18]考虑深部滑脱层估计的地壳缩短速率(5.7±1.5 mm/a)基本一致，但高于GPS观测揭示的汶川地震震前跨龙门山断裂带500 km东西两侧的相对挤压速率(3.5~6.0 mm/a)^{[8, 24]}。一方面，Thompson等^[18]认为，龙门山断裂带深部滑脱层闭锁导致速度梯度向高原内部偏移，即靠近龙门山断裂带区域的速度梯度小，显示出较小的缩短速率。另一方面，杜方等^[24]研究认为，龙门山断裂带西北侧的龙日坝断裂吸收了青藏高原东向挤出中约2 mm/a的地壳缩短作用。尽管本研究模型未包含龙日坝断裂的影响，实际上搜索出的最佳边界挤压速率中可能已经包含龙日坝断裂吸收区域水平缩短作用这个因素，即实际区域相对挤压速率(3.5~6.0 mm/a)加上龙日坝断裂吸收的水平缩短作用(约2mm/a)^[24]，约等于本文搜索出的边界挤压条件。据此，认为搜索到的最佳边界挤压速率是合理的。

在最佳边界挤压速率条件下，经过1.7万余年的计算，模型运行进入稳定状态，共发生15次地震，每个地震周期内的变形特征基本一致，地震复发周期为700~770 a。因此，选取稳定后第8个地震周期的模拟结果，分析龙门山断裂带一个地震周期内不同阶段(震间长期构造加载、同震瞬态变化、震后长期黏弹性松弛调整)的地表变形特征。

2.1 震间速度场

提取模拟震前10 a的平均地表变形速率作为震间速度模拟值，与GPS实际观测值(图 3(a))进行对比。由于模拟的速度场垂直于断层方向，需要将观测值投影到相同方向上进行对比(图 3(b)，图中黑色虚线代表龙门山断裂带所在位置)。如前所述，模拟值与观测值的拟合残差RMS值约为1.3 mm/a，表明模型与计算结果可靠。但在龙门山断裂带附近的H049点位存在很大差异，实际的龙门山断裂带由汶川-茂县断层、映秀-北川断层以及灌县-安县断层3条断层组成，地表覆盖范围约30~50 km，而本文模型将龙门山断裂带简化为一条断层，可能导致近场站点的拟合残差相对较大。

2.2 同震变形

图 4(a)为2008年汶川地震同震水平位移GPS观测值，红色沙滩球为2008年汶川地震震源机制解，蓝色和灰色箭头表示同震位移场^[12]，蓝色箭头表示本研究所选取的观测点信息，灰色箭头信息未被利用。图 4(a)右上角子图为红色方框内汶川地震近场区域观测值，其中蓝色箭头为靠近剖线位置处的站点观测值。基于上述模型和参数条件计算得到的同震模拟值与观测值如图 4(b)所示，其中蓝色圆点为观测值，红色曲线和红色方块分别表示剖线上及观测点位上的同震模拟值，计算得到二者的误差RMS值约为0.151 m。其中，四川盆地一侧的误差RMS值为0.110 m，表现出较高的一致性。川西高原一侧远离龙门山断裂带约100 km的区域，误差RMS值为0.091 m，近场区域的误差RMS值为0.204 m, 尤其是近场60 km以内的区域模拟误差最大。

断层同震滑动分布模拟结果如图 5所示，在近地表位置滑动量最大，约为7 m；15 km深度以上(倾角55°)滑动量为6~7 m；15 km深度以下滑动量随深度衰减。前人开展的汶川地震动态破裂模拟结果显示，汶川地震逆冲型破裂主要集中在龙门山断裂带西南段震中附近，倾滑占比超过60%，同震逆冲滑动量最大值在5~8 m之间，且断层浅部(约13~18 km深度以上)有较大的滑动量，深部滑动量快速衰减^{[17, 22]}。基于大地测量观测资料反演的断层逆冲滑动分布与动态破裂模拟结果类似，其中Feng等^[25]使用GPS与InSAR联合反演得到的逆冲滑动量最大值约为6.7 m，滑动量随深度增加近似均匀减小；本文模拟值在深浅部有不一致的变化特征；Wang等^[12]使用GPS观测数据反演的断层滑动分布更为复杂，最大逆冲滑动量约为10 m(对应7.2 km深度)。主要原因在于断层的几何形态设定不同，导致计算结果存在差异。此外，上盘距离断层40 km左右的站点，如H045、W109等，模拟结果与观测数据拟合度相对较差，可能是该点本身的地表变形存在右旋分量。

2.3 震后变形

图 6(a)为汶川地震震后6年GPS观测的地表水平向累积位移，模拟值与观测值对比如图 6(b)(不同颜色曲线代表不同时间的累积变形)所示。震后形变主要集中在松潘-甘孜块体一侧，四川盆地一侧震后变形极小，这是由两侧中下地壳流变性质的巨大差异所致。其中近断层约30 km范围内的模拟值与观测值有较大误差，远场模拟值与观测值的拟合度则较高。主要原因在于余滑效应是震后变形的主要机制之一，且主要集中在近场区域和震后短期内，而本研究采用的方法未包含余滑的贡献。汶川地震震后形变的数值模拟结果^{[19, 26-27]}显示，余滑效应在震后1 a造成最大约50 mm的地表变形，占该阶段总变形量的60%，随后逐渐减小，且余滑的贡献主要在近场60 km范围内；在后续时间，黏弹性松弛效应逐渐占据主导地位，且控制远场变形。本研究结果显示, 在最靠近断层的CAOP和WENC两个站点，震后6 a的累积位移模拟值与观测值的误差最大，差异分别约为80 mm和20 mm，约等于余滑效应在这2个站点的贡献值。其中H037、H046和MALO等中远场站点的余滑效应较小，主要受黏弹性松弛效应控制，因此本研究中这3个站点的模拟值与观测值吻合较好。对于H028站点，尽管震后余滑贡献几乎为0，黏弹性松弛占主导，但是受模型计算边界效应的影响，导致误差较大。

3 讨论

收起

3.1 摩擦系数和初始应力的影响

根据滑移弱化摩擦准则，震间断层闭锁阶段，断层上的剪应力在构造加载条件下持续积累，在地震发生前达到τ=μ_sσ_n；地震发生时，断层错动量超过滑移弱化距离d₀，断层即开始发生快速破裂，摩擦系数由静摩擦系数μ_s降低为动摩擦系数μ_d，剪应力τ=μ_dσ_n，由此可得到同震释放的剪应力Δτ=(μ_s－μ_d)σ_n，即为该次地震的应力降。在相同的模型和边界条件下，断层面上正应力的积累速率相同，应力降只由断层的静摩擦系数与动摩擦系数的差值决定，该差值越大，单次地震产生的应力降则越大，即地震震级越大，同震地表位移与断层滑动量也越大，同时地震复发周期也会更长。

此外，断层倾角的大小也会影响正应力和剪应力的分配。倾角越大，投影到断层面上的正应力越大，剪应力越小。因此，高倾角断层发生逆冲型破裂滑动需要的最大主应力和最小主应力之间的差值比相同条件下的走滑和拉张断层更大，且摩擦系数需要更小，才能产生破裂滑动^[28]。由于汶川地震孕震环境特殊，且破裂包含逆冲和右旋走滑分量，本研究中仅考虑逆冲破裂，为简化计算，直接设置断层面上的初始正应力不随深度变化，浅层倾角较大设置较高剪应力，深部剪应力较低。在该条件下，模型能够快速实现第一次破裂。

3.2 地震复发周期的讨论

本文估计的地震破裂复发周期为700~770 a。尽管低于目前多数研究的结论(2 000 a~10 000 a)，但与考虑龙门山断裂带映秀段深部滑脱层影响的估计结果(约600 a)大致相当^[18]。

对地震复发周期的研究，一类是基于弹性回跳理论，利用同震滑动量与断层长期滑动速率^{[14, 18]}，或者同震应力释放量与应力积累速率的关系^[29]，范围为600~6 000 a。假设断层长期处于闭锁状态，发生突然错动，释放能量。这个过程中假设断层的滑动速率在时间尺度上是完全均匀的，以及破裂是一次性将所有的能量全部释放。但是实际上断层滑动速率存在时空上的不均匀性，震级和同震位移之间也并非线性关系，每次地震可能破裂的区域不同，造成震级和滑动量也存在差异。

另一类研究是基于动力学数值模拟实验^{[15, 17]}，估计的地震复发周期为3 000~5 000 a。本文结果与之有较大差异，主要原因是边界挤压速率和摩擦系数的取值不同。朱守彪等^[15]采用速率相关摩擦准则，以古地震调查给出的地震复发周期2 000~6 000 a为约束，当断层面上的静摩擦系数取0.6时，模拟的地震复发周期平均值为3 000 a。Zhang等^[17]建立三维非平面断层模型，基于状态-速率摩擦准则，探讨其对地震破裂过程和长期地震序列的控制作用，设置构造加载速率随断层走向变化(0.8~2.2 mm/a)，认为地震复发间隔主要受断层几何和构造加载速率调控，南段破裂周期约为4 300 a，全段破裂周期约为13 000 a。但是这些研究成果缺乏同震、震间和震后等观测数据的约束。

首先，边界挤压速率是驱动断层应力积累的动力源，直接影响地震复发周期，加载速率越快，达到破裂阈值的时间越短。本研究以震间GPS观测数据为约束，在确定的模型参数条件下，搜索出最佳边界挤压速率为6.9 mm/a，与Thompson等^[18]的结果一致，因此得到的地震复发周期也接近。而Zhang等^[17]设置的加载速率仅为0.8~2.2 mm/a，必然会导致复发周期更长。其次，本研究中采用滑移弱化摩擦准则，主要由静摩擦系数和动摩擦系数控制断层破裂滑动。静摩擦系数越高，表示断层的闭锁程度越强，应力积累越多，可能复发周期越长；动摩擦系数会影响破裂传播速度和应力释放效率。二者的差值越大，同震应力降越高，即同震滑动量越大。本文给出的静摩擦系数略低，参考Zhang等^[22]模拟同震破裂传播使用的静摩擦系数，设置为0.35，可能导致复发间隔较长。再者，其他未考虑的影响因素，如区域历史地震对龙门山断裂带的应力扰动效应等，也可能影响地震复发周期。已有研究认为，鲜水河断裂带上300余年M>6.5的历史地震，对龙门山断裂带中南段具有显著的应力卸载效应^[30]。但是本研究模型未考虑这个因素，可能使得模拟的地震复发周期偏短。最后，本研究建立的二维纵剖面模型只能模拟龙门山断裂带映秀段纯逆冲型地震的破裂，而2008年汶川地震存在一定的右旋走滑分量，且随破裂向东北延伸逐渐变为以右旋走滑为主。该二维模型是通过边界挤压实现断层周期性自发破裂，无法考虑走滑向的影响，给出的是龙门山断裂带映秀段逆冲型地震破裂的复发周期和变形，无法完整描绘汶川地震的真实破裂，且忽略了北川附近另一个主滑动区域的应力释放，因此与前人基于古地震调查或长期滑动亏损得到的地震复发周期有一定差别。

总之，尽管二维模型的计算更高效，也得到部分有意义的结果，但是仍存在较多的缺陷。由于汶川地震破裂因素复杂，需要综合多种因素进行互相约束和校正，以及可能需要采用更精准的三维计算才能获取更高精度的结果。

4 结语

收起

本研究建立一个垂直于龙门山断裂带西南段的二维纵剖面黏弹性模型，基于滑移弱化摩擦准则，模拟逆冲型地震自发破裂全周期过程，探讨影响地震复发周期和地震周期不同阶段地表形变特征的主要动力学因素，研究结论如下：

1) 震间变形及地震复发周期主要受边界挤压速率影响，以震前GPS观测的地表水平速度为约束，网格搜索得到模型最佳边界挤压速率为6.9 mm/a，与考虑深部滑脱层影响而估计的结果大致相当。值得注意的是，本研究二维模型西边界距离龙门山断裂带300 km，跨过龙日坝断裂，该边界挤压速率指的是综合多种因素的模型西边界相对四川盆地的挤压速率。

2) 静摩擦系数和动摩擦系数的差值控制着断层滑动量及地表同震变形，最优静、动摩擦系数分别为0.35和0.23，模拟得到断层同震滑动最大值约为7 m。在震间、同震和震后与观测值对比一致性较好的情况下，地震复发周期为700~770 a。

3) 铲形断层相比平直断层，破裂在地表产生的同震水平向位移更符合实际GPS观测结果，且断层倾角越大意味着断层上盘同震垂向位移越大。此外，断层倾角也会影响断层上的正应力和剪应力的积累速率，从而影响地震复发周期。

基金

收起

国家自然科学基金(42174113)
中国地震局地震研究所和应急管理部国家自然灾害防治研究院基本科研业务费(IS202426355)
应急管理部国家自然灾害防治研究院基本科研业务专项(ZDJ2024-30)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Jiang

J L

, Erickson

B A

, Lambert

V R

, et al. Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip[J]. Journal of Geophysical Research: Solid Earth, 2022, 127 (3)

Barbot

, Lapusta

, Avouac

J P

. Under the Hood of the Earthquake Machine: Toward Predictive Modeling of the Seismic Cycle[J]. Science, 2012, 336 (6082): 707- 710

Barbot

. Asthenosphere Flow Modulated by Megathrust Earthquake Cycles[J]. Geophysical Research Letters, 2018, 45 (12): 6018- 6031

H Y

, Liu

Y J

, Yang

H F

, et al. Modeling Earthquake Sequences along the Manila Subduction Zone: Effects of Three-Dimensional Fault Geometry[J]. Tectonophysics, 2018, 733, 73- 84

Dal Zilio

, Lapusta

, Avouac

J P

, et al. Subduction Earthquake Sequences in a Non-Linear Visco-Elasto-Plastic Megathrust[J]. Geophysical Journal International, 2022, 229 (2): 1098- 1121

张培震, 闻学泽, 徐锡伟, 等. 2008年汶川8.0级特大地震孕育和发生的多单元组合模式[J]. 科学通报, 2009, 54 (7): 944- 953

Zhang

Peizhen

, Wen

Xueze

, Xu

Xiwei

, et al. Tectonic Model of the Great Wenchuan Earthquake of May 12, 2008, Sichuan, China[J]. Chinese Science Bulletin, 2009, 54 (7): 944- 953

Wang

, Zhang

P Z

, Freymueller

J T

, et al. Present-Day Crustal Deformation in China Constrained by Global Posi tioning System Measurements[J]. Science, 2001, 294 (5542): 574- 577

Zhang

P Z

, Shen

Z K

, Wang

, et al. Continuous Deformation of the Tibetan Plateau from Global Positioning System Data[J]. Geology, 2004, 32 (9): 809- 812

Gan

W J

, Zhang

P Z

, Shen

Z K

, et al. Present-Day Crustal Motion within the Tibetan Plateau Inferred from GPS Measurements[J]. Journal of Geophysical Research: Solid Earth, 2007, 112 (B8)

Wang

, Shen

Z K

. Present-Day Crustal Deformation of Continental China Derived from GPS and Its Tectonic Implications[J]. Journal of Geophysical Research: Solid Earth, 2020, 125 (2)

Shen

Z K

, Sun

J B

, Zhang

P Z

, et al. Slip Maxima at Fault Junctions and Rupturing of Barriers during the 2008 Wenchuan Earthquake[J]. Nature Geoscience, 2009, 2 (10): 718- 724

Wang

, Qiao

X J

, Lan

Q G

, et al. Rupture of Deep Faults in the 2008 Wenchuan Earthquake and Uplift of the Longmen Shan[J]. Nature Geoscience, 2011, 4 (9): 634- 640

Densmore

A L

, Ellis

M A

, Li

, et al. Active Tectonics of the Beichuan and Pengguan Faults at the Eastern Margin of the Tibetan Plateau[J]. Tectonics, 2007, 26 (4)

张培震, 徐锡伟, 闻学泽, 等. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因[J]. 地球物理学报, 2008, 51 (4): 1066- 1073

Zhang

Peizhen

, Xu

Xiwei

, Wen

Xueze

, et al. Slip Rates and Recurrence Intervals of the Longmen Shan Active Fault Zone and Tectonic Implications for the Mechanism of the May 12 Wenchuan Earthquake, 2008, Sichuan, China[J]. Chinese Journal of Geophysics, 2008, 51 (4): 1066- 1073

朱守彪, 张培震. 2008年汶川M_S8.0地震发生过程的动力学机制研究[J]. 地球物理学报, 2009, 52 (2): 418- 427

Zhu

Shoubiao

, Zhang

Peizhen

. A Study on the Dynamical Mechanisms of the Wenchuan M_S8.0 Earthquake, 2008[J]. Chinese Journal of Geophysics, 2009, 52 (2): 418- 427

赵静, 江在森, 武艳强, 等. 汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究[J]. 地球物理学报, 2012, 55 (9): 2963- 2972

Zhao

Jing

, Jiang

Zaisen

, Wu

Yanqiang

, et al. Study on Fault Locking and Fault Slip Deficit of the Longmenshan Fault Zone before the Wenchuan Earthquake[J]. Chinese Journal of Geophysics, 2012, 55 (9): 2963- 2972

Zhang

, Liu

Y J

, Li

, et al. Geometric Control on Seismic Rupture and Earthquake Sequence along the Yingxiu-Beichuan Fault with Implications for the 2008 Wenchuan Earthquake[J]. Journal of Geophysical Research: Solid Earth, 2022, 127 (12)

Thompson

T B

, Plesch

, Shaw

J H

, et al. Rapid Slip-Deficit Rates at the Eastern Margin of the Tibetan Plateau Prior to the 2008 M_W7.9 Wenchuan Earthquake[J]. Geophysical Research Letters, 2015, 42 (6): 1677- 1684

Wang

, Shen

Z K

, Wang

Y Z

, et al. Postseismic Deformation of the 2008 Wenchuan Earthquake Illuminates Lithospheric Rheological Structure and Dynamics of Eastern Tibet[J]. Journal of Geophysical Research: Solid Earth, 2021, 126 (9)

Zhu

Y G

, Diao

F Q

, Wang

R J

, et al. Crustal Shortening and Rheological Behavior across the Longmen Shan Fault, Eastern Margin of the Tibetan Plateau[J]. Geophysical Research Letters, 2022, 49 (11)

Andrews

D J

. Rupture Velocity of Plane Strain Shear Cracks[J]. Journal of Geophysical Research, 1976, 81 (32): 5679- 5687

Zhang

Z G

, Zhang

W Q

, Chen

X F

. Dynamic Rupture Simulations of the 2008 M_W7.9 Wenchuan Earthquake by the Curved Grid Finite-Difference Method[J]. Journal of Geophysical Research: Solid Earth, 2019, 124 (10): 10565- 10582

Aagaard

B T

, Knepley

M G

, Williams

C A

. A Domain Decomposition Approach to Implementing Fault Slip in Finite-Element Models of Quasi-Static and Dynamic Crustal Deformation[J]. Journal of Geophysical Research: Solid Earth, 2013, 118 (6): 3059- 3079

杜方, 闻学泽, 张培震, 等. 2008年汶川8.0级地震前横跨龙门山断裂带的震间形变[J]. 地球物理学报, 2009, 52 (11): 2729- 2738

Fang

, Wen

Xueze

, Zhang

Peizhen

, et al. Interseismic Deformation across the Longmenshan Fault Zone before the 2008 M8.0 Wenchuan Earthquake[J]. Chinese Journal of Geophysics, 2009, 52 (11): 2729- 2738

Feng

G C

, Hetland

E A

, Ding

X L

, et al. Coseismic Fault Slip of the 2008 M_W7.9 Wenchuan Earthquake Estimated from InSAR and GPS Measurements[J]. Geophysical Research Letters, 2010, 37 (1)

Diao

F Q

, Wang

R J

, Wang

Y B

, et al. Fault Behavior and Lower Crustal Rheology Inferred from the First Seven Years of Postseismic GPS Data after the 2008 Wenchuan Earthquake[J]. Earth and Planetary Science Letters, 2018, 495, 202- 212

Zhao

D Z

, Qu

C Y

, Shan

X J

, et al. Post-Seismic Deformation of the 2008 Wenchuan Earthquake Reveals a Misaligned Rheological Boundary in the Lower Crust along the Eastern Tibetan Plateau Margin[J]. Geophysical Journal International, 2023, 235 (2): 1353- 1372

周永胜, 何昌荣. 汶川地震区的流变结构与发震高角度逆断层滑动的力学条件[J]. 地球物理学报, 2009, 52 (2): 474- 484

Zhou

Yongsheng

, He

Changrong

. The Rheological Structures of Crust and Mechanics of High-Angle Reverse Fault Slip for Wenchuan M_S8.0 Earthquake[J]. Chinese Journal of Geophysics, 2009, 52 (2): 474- 484

Liu

, Dong

P Y

, Shi

Y L

. Recurrence Interval of the 2008 M_W7.9 Wenchuan Earthquake Inferred from Geodynamic Modelling Stress Buildup and Release[J]. Journal of Geodynamics, 2017, 110, 1- 11

Liu

, Dong

P Y

, Zhu

B J

, et al. Stress Shadow on the Southwest Portion of the Longmen Shan Fault Impacted the 2008 Wenchuan Earthquake Rupture[J]. Journal of Geophysical Research: Solid Earth, 2018, 123 (11): 9963- 9981

2026年第46卷第8期

引用本文

BibTeX

文章信息

doi: 10.14075/j.jgg.2025.12.407

接收时间：2025-12-02
首发时间：2026-10-30
出版时间：2026-08-15

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-12-02

基金

国家自然科学基金(42174113)

中国地震局地震研究所和应急管理部国家自然灾害防治研究院基本科研业务费(IS202426355)

应急管理部国家自然灾害防治研究院基本科研业务专项(ZDJ2024-30)

作者信息

1 中国地震局地震研究所, 武汉, 430071

2 中国地震局地震大地测量重点实验室, 武汉, 430071

通讯作者:

董培育, 博士, 副研究员, 主要研究方向为地震大地测量学与地球动力学数值模拟, E-mail: dongpeiyu97@163.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/ddcl/CN/10.14075/j.jgg.2025.12.407

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT