弯折断层系统上超剪切破裂传播的产生条件

断层的破裂速度是描述地震震源过程的重要物理量.如果震源破裂的传播速度超过剪切波速,将会对地震波场产生影响,造成更大的破坏性.超剪切破裂的产生受多种因素影响,断层的几何形状是因素之一.本文针对弯折断层的情况,采用三维空间非结构化网格的边界积分方法计算参数空间中的破裂相图,从中分析超剪切破裂的产生条件.以15°、25°和40°为例,得到了不同断层弯折角度的破裂相图.在本文的初始应力设置下,通过对不同的无量纲化临界滑动弱化位移Dc和初始剪应力Te参数组合的结果进行交叉对比发现,对于弯折面处于压缩区的断层模型,不可持续传播的自发停止破裂的发生条件与弯折角无关.而对于可持续传播破裂,其在平面断层的传播速度也不受弯折角影响;在弯折部分,随弯折角度增大,破裂传播速度越小,正应力越大,破裂强度越大,破裂越难以越过弯折交界线继续传播（如40°）.对比三个不同弯折角的相图,弯折角越小,越容易发生超剪切破裂,即发生超剪切的参数空间越大.同时,随着初始剪应力的增大,超剪切不仅可以发生在弯折面上,甚至在平面部分就可以发生.总体而言,Dc较小、Te较大时,破裂传播速度更大,更容易形成超剪切破裂.另外,因克服弯折交界处的正应力而产生的错位延迟效应也与弯折角度正相关.     

角度对于Y型分叉断层自发破裂传播过程的影响

断层的几何形态是地震破裂传播过程的控制因素之一,从而影响着地震的危险性.Y型分叉断层是断层的多种复杂几何形态当中常见的一种,研究断层的分叉特征对震源破裂传播的影响,对于深入认识复杂几何形态断层的动力学特征具有重要意义.本文利用边界积分方程方法,模拟了612个分叉面,通过改变分叉断层面的分叉角度来分析断层分叉对传播的影响,并定量分析了分叉面之间以及主断层对分叉面的应力作用机制.模拟结果表明：分叉断层的一个断层分叉的破裂,既受到主断层作用的影响,也受到该断层另一分叉的作用的影响,是两者共同作用的结果.其中,主断层的作用基本上只与该断层分叉与主断层延长面的夹角有关,而与该断层另一分叉关系不大;该断层另一分叉的作用主要与两个断层分叉的夹角有关,但同时也要考虑其破裂状况.对于破裂以超剪切速度到达断层分叉处的情况,主断层对于大角度和小角度分叉的破裂促进作用较强,而对中等角度分叉的破裂促进作用较弱;该断层另一分叉对破裂的作用随着两个断层分叉之间夹角的增加,由强烈的抑制转为促进.

     

局部山体地形对断层动力学破裂过程的影响研究

基于曲线网格有限差分方法研究了垂直走滑断层在不同山体地形情况下的动力学破裂传播,分析并讨论了局部山体地形对断层破裂过程及相应地面地震动的影响,得到了各模型断层面的动力学破裂过程及相应的地表峰值速度特征。研究结果表明,山体地形尺度（山体高度及底部延展距离等）对断层动力学破裂过程影响较大,进而影响到相应的地面地震动分布。当山体地形处于自由地表上亚剪切向超剪切转换的位置附近时,山体地形会阻碍断层面上自由地表超剪切的产生。一般而言,对于具有一定埋深的断层,当山体地形底部延展距离一定时,山体高度越高,其对自由地表超剪切的阻碍程度越大;当山体高度一定时,地形底部延展距离越大,越会阻碍自由地表超剪切的产生,这种破裂过程的变化会导致相应地面地震动呈现不同特征的分布。此外,还探讨了断层破裂过程及相应地震动对成核区外初始剪切应力变化的响应,结果显示,当初始剪切应力较高时,高应力降引起的超剪切破裂会对断层破裂及相应的地震动分布起主导作用。      

青海玉树地震(M_s=7.1)产生超剪切破裂过程的动力学机制研究

由于2010年玉树地震(Ms=7.1)产生了超剪切地震破裂,所以地震灾害特别严重.国内外地球科学家对该地震产生超剪切破裂过程的物理机制一直非常关注,但至今没有给出满意的解答.为此,文中根据玉树地震发震断层的实际几何构建有限单元数值模型,模型中的断层由2个断层段构成,它们之间有约10°的夹角,形成断层拐折.模拟结果表明,玉树地震的破裂由2个子事件组成;当破裂在震源所在的断层上成核后,先在第一个断层段上传播,其速度为亚剪切波速度;当破裂一旦越过断层拐折,在第二个断层段上传播时,破裂速度就立即转变为超剪切波速度.计算结果显示,当断层发生超剪切破裂时,断层上的位错幅度、破裂产生的地震波速度及加速度都会显著增大,从而造成地震灾害大大增加,这很可能是玉树地震的震害特别严重的重要原因.从模拟实验中还看到,若是模型中的断层没有发生拐折,在模型的其他参数都保持不变的情况下,破裂速度不会发生变化.但是,若初始应力场的方位与断层之间的夹角发生变化,这时断裂系统中尽管存在断层拐折,也不是一定能产生超剪切破裂.只有当初始应力方位与断层之间的夹角以及断层走向变化的偏角二者之间的关系恰到好处时,断层拐折才有可能促使断层破裂由亚剪切转化为超剪切破裂.所以,玉树地震之所以能产生超剪切地震破裂,恰恰是发震断层几何与初始应力场方位之间的关系达到某种"最佳状态"的结果.这也可能是天然地震中超剪切破裂事件稀少的原因之一.因此,研究超剪切地震破裂过程的动力学机制,对于深入研究地震震源过程、地震灾害评估等有着非常重要的科学意义.     

2008年汶川大地震中北川地区极重震害的物理机制研究

2008年汶川特大地震中北川地区受灾特别严重,而该地区远离汶川主震震中超过100多千米.尽管汶川地震发生已近10年,但这种不正常现象一直困扰着地学工作者,至今没有给出合理的解答.为此,本文利用有限单元方法模拟汶川地震的主要发震断裂带（即：映秀—北川断层）的自发破裂动力学过程.模型中,几乎平直的映秀—北川断裂带在高川地区发生拐折,形成高川右弯.模拟结果显示,破裂沿着映秀—北川断裂带向东北方向前行,当破裂前端到达高川右弯时,破裂形态发生剧烈变化,破裂传播速度由亚剪切破裂突然转化为超剪切破裂.与断层阶区促进超剪切破裂明显不同的是,这种超剪切破裂形态的转化不需要时间停顿.由于超剪切地震破裂的产生,破裂辐射的地震波发生相长干涉,形成马赫波,地震动被大大放大.计算给出的强地面运动峰值加速度在北川地区不仅数值高,而且高值区分布范围广,造成了北川地区的震害特别严重.模拟结果还表明,若高川右弯不存在（即映秀—北川断层面为平直）,则北川地区的断层不会形成超剪切破裂;如果高川右弯附近断层的不连续程度较大,则破裂会被终止;这些情况下,北川地区的震害都没有那么严重.所以,高川右弯的几何形状对于超剪切破裂的形成起着决定性作用,也是造成远离震中的北川地区震害特别严重的主要原因.因此,分析发震断层几何对于研究震源动力学过程及震害评估等有着非常重要的科学意义和实际应用价值.

     

断层自发破裂动力学过程的有限单元法模拟及其在地震研究中的应用

地震是断层的自发破裂动力学过程。数值模拟断层的自发破裂动力学过程对于认识地震的力学本质、减轻地震灾害等有着重要的科学意义及应用价值。本文首先对经典的滑移弱化摩擦关系进行了改进,然后对断层的破裂过程进行动态数值模拟。模拟结果表明,利用改进后的摩擦关系能够产生脉冲型(pulse-like)破裂模式。断层自发破裂过程受初始应力场及摩擦关系影响,若初始应力场中的剪应力水平较低或滑移弱化摩擦本构关系中的动摩擦系数较大,则容易产生脉冲型破裂;反之,则容易产生裂纹型(crack-like)破裂。另外,为了研究双材料(bimaterial)断层破裂对强地面运动的影响,我们采用正则化的速率-状态相关摩擦本构关系计算了破裂沿着双材料断层传播的二维有限元模型。模拟结果表明,双材料机制对地震破裂过程以及断层周边区域的强地面运动有显著影响。由断层破裂辐射出的地震波导致的强地面运动在整个空间上的分布是不对称的,其不对称性会随着断层两侧材料差异程度的增加而增加。断层破裂能否跨越断层阶区(stepover)继续传播,从而引发更大震级的地震,地震时断层是否发生超剪切破裂导致地震灾害加剧,都是震源动力学研究的重要内容。本文利用有限单元方法模拟断层阶区对地震破裂传播的控制作用以及对产生超剪切地震破裂的促进作用。研究结果表明:断层面上的摩擦系数减小、断层周边区域内初始剪应力增大以及较小的阶区间距等,都将增加断层破裂跳跃阶区传播的可能性;此外,这些物理因素都会对破裂的传播速度产生影响。在一定条件下,破裂传播速度会由在初始断层上的亚剪切波速度转为在次级断层上的超剪切波速度。结合以上在概念模型中对断层自发破裂过程的模拟研究结果,我们根据汶川地震和玉树地震发震断层的实际几何情况分别构建有限单元数值模型,研究了汶川地震单侧破裂过程的动力学机制以及玉树地震产生超剪切破裂过程的动力学机制。2008年汶川大地震的破裂过程极其复杂,向东北方向的破裂距离长达300 km,而向西南方向的破裂长度很小,呈现出单侧破裂的主要特征。文中模拟并分析了汶川地震的破裂过程,结果表明:龙门山断裂带两侧的物性差异是造成汶川大地震单侧传播的决定性因素。由于2010年玉树地震(Ms=7.1)产生了超剪切地震破裂,所以地震灾害特别严重。文中在模拟并分析玉树地震的破裂过程后认为:玉树地震发震断层走向与初始主应力方向之间的关系断层破裂是亚剪切转化为超剪切破裂的可能原因。     

断层破裂速度对地震动及其地震灾害影响的有限单元法模拟

断层破裂传播速度通常会影响强地面运动加速度的空间分布及地震的灾害程度,但究竟是如何影响的,目前未见全面的定量分析.为此,本研究设断层破裂速度从小到大发生改变(从亚剪切波速度到超剪切破裂速度),同时利用有限元方法计算在不同破裂传播速度的情况下,破裂产生的地震动及强地面运动峰值加速度(PGA)的空间分布,计算时保持在所有情况下其断层上的最终位错相等,即所有不同破裂情况下产生的地震矩震级保持不变,都为MW=7.0.计算结果显示：总体上,破裂传播的速度越快,PGA的幅值越大并且高值区分布的范围也越广,其地震灾害会越严重;超剪切破裂相对于亚剪切破裂会加剧地震灾害.但是,当破裂速度正好等于介质的剪切波速度时,断层面附近PGA峰值最大,断层附近的地震灾害最为严重.此外,研究中还发现,当破裂传播速度等于倍剪切波速度时(该速度称为Eshelby速度),不产生剪切马赫波;这与超剪切破裂产生剪切马赫波的普遍现象不一致;但此时PGA并没有大幅度减小,其震害也不会有显著改变.因此,本研究对于深入认识地震震源过程、地震动效应以及灾害评估有着重要的科学意义及实用价值.     

基于边界积分方程方法的弯折断层破裂传播过程控制因素分析

本文利用边界积分方程方法,以基于三角形网格的全空间格林函数及离散积分核计算为基础,进行了最常见的弯折断层的破裂传播过程模拟.为了去除边界积分方程方法中格林函数计算存在的高度奇异性,研究采用分部积分等方法对动力学方程进行了重整化和离散化处理.地震力学过程可以被视为断层由静摩擦转为动摩擦的过程,对于震源破裂过程的动力学模拟,摩擦准则起着重要作用,本研究采用常用的滑动弱化摩擦准则.计算引入Courant-Friedrich-Lewy比值来表达场点的影响,并控制计算的收敛性和稳定性.通过与典型算例的比对,检验了方法的正确性和有效性.地震破裂能否穿越断层弯折部位继续传播是震源动力学研究的重要内容,基于此,本文建立了多种理论弯折断层模型,模拟了断层弯折对地震破裂传播的控制作用,并通过改变断层周边初始应力场、断层弯折角度大小以及滑动弱化距离大小等来分析各个因素对破裂传播的影响.模拟结果表明:断层面上初始破裂区域内外的应力越高,破裂越容易越过断层弯折部位继续传播;初始破裂区域半径越大,或滑动弱化距离越小,破裂也越容易发生,并越过弯折部位继续传播.同样的初始条件,断层弯折角度越大,断层弯折作为障碍体,对破裂传播的阻碍作用越显著.小的弯折角,其破裂传播过程与平面断层差别不明显,基本仍以椭圆方式对称向两侧传播.     

不同摩擦本构关系对断层自发破裂动力学过程的影响

地震是断层的摩擦失稳过程.摩擦本构关系对断层的破裂成核、破裂传播、破裂速度、能量释放、破裂终止等起着至关重要的控制作用.为了比较不同摩擦关系在断层自发破裂动力学过程中的影响,文中引入目前应用最为广泛的4种典型摩擦本构关系,它们分别是：滑移弱化摩擦关系,速率弱化摩擦关系,以及速率-状态相依摩擦关系中的老化定律和滑动定律.研究中利用有限单元方法对上述4种摩擦关系控制的断层自发破裂过程分别进行模拟计算,模拟结果显示：当模型参数相同时,不同摩擦关系模拟的破裂行为总体上具有一致性,都可以产生亚剪切破裂或超剪切破裂,并且破裂传播速度的大小与摩擦本构关系的类型无关.此外,它们之间还存在着较大的差异：（1）速率弱化摩擦关系可以模拟脉冲型破裂;而其他3个摩擦关系只能模拟裂纹性破裂.（2）不同摩擦关系模拟的超剪切破裂转换长度不同,速率-状态相关摩擦关系的老化定律相比其他摩擦关系需要更大的转换长度才能实现亚剪切破裂转变为超剪切破裂;而速率弱化的摩擦关系的超剪切转换长度可以为0,即不需要转换距离就直接产生超剪切破裂.（3）速率弱化摩擦关系模拟的破裂速度自成核后很快就达到稳定值;而其他类型摩擦关系模拟的破裂传播则要经历由缓慢破裂到逐渐加速直至达到稳定破裂的过程.值得特别指出的是,本文所使用的4种摩擦关系都不能完整地反映实际大地震破裂过程的摩擦属性,需要进一步深入研究.

     

不同摩擦本构关系对断层自发破裂动力学过程的影响

地震是断层的摩擦失稳过程.摩擦本构关系对断层的破裂成核、破裂传播、破裂速度、能量释放、破裂终止等起着至关重要的控制作用.为了比较不同摩擦关系在断层自发破裂动力学过程中的影响,文中引入目前应用最为广泛的4种典型摩擦本构关系,它们分别是：滑移弱化摩擦关系,速率弱化摩擦关系,以及速率-状态相依摩擦关系中的老化定律和滑动定律.研究中利用有限单元方法对上述4种摩擦关系控制的断层自发破裂过程分别进行模拟计算,模拟结果显示：当模型参数相同时,不同摩擦关系模拟的破裂行为总体上具有一致性,都可以产生亚剪切破裂或超剪切破裂,并且破裂传播速度的大小与摩擦本构关系的类型无关.此外,它们之间还存在着较大的差异：（1）速率弱化摩擦关系可以模拟脉冲型破裂;而其他3个摩擦关系只能模拟裂纹性破裂.（2）不同摩擦关系模拟的超剪切破裂转换长度不同,速率-状态相关摩擦关系的老化定律相比其他摩擦关系需要更大的转换长度才能实现亚剪切破裂转变为超剪切破裂;而速率弱化的摩擦关系的超剪切转换长度可以为0,即不需要转换距离就直接产生超剪切破裂.（3）速率弱化摩擦关系模拟的破裂速度自成核后很快就达到稳定值;而其他类型摩擦关系模拟的破裂传播则要经历由缓慢破裂到逐渐加速直至达到稳定破裂的过程.值得特别指出的是,本文所使用的4种摩擦关系都不能完整地反映实际大地震破裂过程的摩擦属性,需要进一步深入研究.     

Imaging the rupture process of the 10 January 2018 MW7.5 Swan island,Honduras earthquake*

The 10 January 2018 M_W7.5 Swan island, Honduras earthquake occurred on the Swan island fault, which is a transform plate boundary between the North American and Caribbean plates. Here we back-project the rupture process of the earthquake using dense seismic stations in Alaska, and find that the earthquake ruptured at least three faults (three stages) for a duration of ~40 s. The rupture speed for the longest fault (stage 3) is as fast as 5 km/s, which is much faster than the local shear wave velocity of ~4 km/s. Supershear rupture was incidentally observed on long and straight strike-slip faults. This study shows a supershear rupture that occured on a strike-slip fault with moderate length, implying that supershear rupture might commonly occur on large strike-slip earthquakes. The common occurrence of supershear rupture on strike-slip earthquakes will challenge present understanding of crack physics, as well as strong ground motion evaluation in earthquake engineering.     

断层阶区对产生超剪切地震破裂的促进作用

地震时若断层发生超剪切破裂,地震灾害会显著加剧.因此研究超剪切破裂的形成机理有着非常重要的科学意义.本文利用动力有限单元方法,模拟断层破裂从初始断层跳跃传播到另一条平行的次级断层（断层阶区）时破裂速度的变化情况,并分析断层阶区几何特征等物理参数对产生超剪切地震破裂的促进作用.计算结果表明,断层阶区的诸多物理因素（如：重叠长度、相隔距离以及摩擦系数等）都会对破裂的传播速度产生影响.在一定条件下,破裂传播速度会由在初始断层上的亚剪切波速度,转换为在次级断层上的超剪切波速度.在破裂速度转换过程中,断层间隔起着重要作用,当断层阶区中两断层垂直间隔距离小到一定程度时,破裂跳跃阶区后,破裂速度不会发生变化.所以对于分段断层（可视为一种特殊的断层阶区）,由于其断层垂直间隔为0,也就不会出现破裂速度变化的现象,模拟结果对此也进行了证实.然而,若断层间隔太大,当其距离超过一定的限度后,破裂通常无法跨越断层阶区继续传播,而是终止在初始断层上.模拟结果还表明,初始断层与次级断层之间的重叠距离也十分重要,只有当断层阶区中两平行断层之间的重叠部分达到一定长度后,断层的破裂速度才有可能发生转换.此外,计算结果显示,破裂过程中断层面上的应力变化可能是破裂传播速度发生转换的直接原因.最后,模拟还发现,当破裂跨越断层阶区发生速度转换时,破裂需要停顿一定的时间,以便积聚足够的能量来实现破裂速度的增快.

     

1999年山西大同Ms 5.6地震的震源断层

大同震区先后在 1989、1991和 1999年发生MS >5地震 ,利用大同遥测地震台网的记录资料进行比较精确的地震序列震源定位 ,结合宏观烈度分布和震源机制解资料 ,详细地分析对比了 3次子序列的异同。结果显示 ,1999年MS5 .6地震的震源断层是走向NWW、长 16km、宽12km、埋深 5km以下、倾角近直立的左旋走滑断层。而前 2个子序列是NNE为主的右旋走滑断层活动所致 ,表明地震破裂方向发生了变化。这种 2个以上方向先后出现、并且强弱有别的地震破裂是普遍存在的 ,表明震源环境的复杂程度与地震序列的类型有关。虽然震区存在NE向的大王村断裂和NW向的团堡断裂 ,但目前没有证据说明震源断层和 2条构造断层连通。 3次子序列的震源断层都是走滑断层 ,也和 2条构造正断层有别。 1999年的子序列可能属于新破裂。     

2014年2月12日M_w6.9于田地震震源破裂过程及对周围断层的应力影响

2014年2月12日,在新疆于田县发生了里氏7.3级地震.在该地震震中附近,前人研究证明发育了大量规模不同的活动断层(如康西瓦断裂与贡嘎错断裂等).根据地震触发理论,地震发生后因地壳同震变形会导致其周边不同性质断裂破裂应力发生变化,进而影响其地震的潜在危险性.本文利用地震远场波形记录,反演了该地震滑动模型.之后,根据弹性无限半空间位错理论,计算了该地震在近场范围内活动断裂上的同震应力变化.其目的在于讨论于田地震引起的附近断裂上的库仑应力变化以及这些活动断裂可能潜在的地震危险性.在地震发生后,从国际地震学联合会(IRIS)地震数据中心,下载了震中距离介于30°~90°的地震远场波形记录,为保证台站方位角分布均匀,从中挑选了27个不同方位角的高信噪比地震记录参与理论地震图的生成和波形反演过程.我们采用广义射线理论计算生成远场理论地震波形数据.每个子断层参数的反演则利用基于全局化反演的快速模拟退火反演方法.在有限断层反演过程中,我们采用了强调波形拟合的相关误差函数作为待反演的目标函数,拟合的断层参数使目标函数为最小.之后,根据弹性无限半空间位错理论,以库仑破裂准则为基础,结合反演得到的地震震源机制解和地震位错模型,计算该地震引起的近场断层面上库仑应力的变化.由远场波形计算结果可以看到,于田地震的震源深度为10km,地震断层的倾角约71.9°,破裂面上最大的同震位移达到210cm,以左旋走滑为主并具有正倾滑分量,地震能量主要在前15s内释放.由此得到该地震的地震矩为2.91×1019 N·m,地震震级为Mw6.9.于田地震引发的余震,大致分布在三个区域内:普鲁断裂北部、康西瓦断裂东部和贡嘎错断裂中部.弹性应力计算结果表明,于田地震导致阿尔金断裂西段、普鲁断裂中段、康西瓦断裂东段和贡嘎错断裂中段的静态库仑应力明显增加,其中以康西瓦断裂东段和贡嘎错断裂中段应力增量为最大,分别达到了0.05 MPa和0.04 MPa.大量研究证明,当地震所导致的库仑应力变化大于0.01 MPa时将具有明显的地震触发作用.根据本文结果,2014年于田Mw6.9地震使普鲁断裂、贡嘎错断裂和康西瓦断裂上的库仑应力增量均超过了触发阈值,具有被触发出地震的潜在危险.因此,在以后的地震学研究中,应加强对该三条断裂地震危险性的研究和监测.此外,近6年以来,研究区域发生了3次6级以上的地震.这些地震均沿着贡嘎错断裂,由南西向北东迁移,逐步靠近阿尔金断裂,并且逐渐由正倾滑型地震转变为走滑型地震.阿尔金断裂的走滑速率达到了9mm·a-1,所以,尽管本次地震导致的阿尔金断裂库仑应力增量小于0.01 MPa,阿尔金的地震危险性也应该加强监测.     

基于时间域半空间Green函数建立的边界积分方程及其在动态破裂模拟中的应用

研究地震断层的自发破裂传播可以为通过震源过程了解地球介质应力状态提供理论基础,有助于理解震源破裂的机理和规律,在震源动力学研究中具有重要的地位.边界积分方程方法是研究地震断层破裂传播问题的常用方法之一,在处理复杂断层系统方面具有优势.本文从半无限空间Green函数的时间域积分形式解出发,推导出离散化的边界积分方程积分核,并将其分别应用到平面断层和弯折断层两种断层模型的动力学破裂过程的模拟中.模拟结果表明,当破裂传播到自由表面处时,会产生反射震相,反射震相与直达震相的耦合效应会导致自由表面处的破裂传播速率明显加快,同时这种耦合效应也可以促进破裂在断层弯折部位的传播.对于断层有一定埋深的情况,自由表面对破裂过程的影响较小,因而自由表面的效应仅在断层直接与地表相交或者断层上沿距离地表很近时需要纳入考虑范围.本文将自由表面的影响加入到断层的自发破裂传播研究中,有助于更真实地了解复杂断层系统中的动力学破裂过程.

     

基于非结构化网格的边界积分方程方法的断层自发破裂模拟

地震自发破裂模拟是震源动力学研究的重要内容,了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用,大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征,为了更为灵活地刻画断层几何复杂性,我们建立断层模型的三角形网格离散方案,通过精确的解析解形式来计算断层各个单元之间的应力格林函数,联立滑动弱化摩擦准则和非奇异边界积分方程,对断层的自发破裂过程进行了模拟.在简单的平面断层模型下,将计算结果与前人的结果进行了对比,验证了方法的正确性与有效性.对于几种常见的复杂断层模型,例如弯折、阶跃、含障碍体断层等,我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明,非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层,同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景.