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

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

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

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

     

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

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

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

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

     

2008年汶川大地震单侧破裂过程的动力学机制研究

2008年汶川大地震的破裂过程极其发杂,向东北方向的破裂距离长达300 km,而向西南方向的破裂长度很小,呈现出单侧破裂的主要特征.尽管汶川地震破裂呈单侧传播的现象引起许多地震学家的关注,但其物理机制至今还不是十分清楚.本文利用有限单元计算方法,模拟了汶川地震的破裂过程.模型中根据龙门山断裂带两侧（东南侧为四川盆地,西北侧为川西高原）实际的地震波速度来确定模型的介质物性参数,利用目前观测的应力环境来选定初始应力条件.模拟结果表明：破裂在汶川地震的震中处成核后,先向断层两侧自发传播,但向东北方向的传播距离明显大于向西南方向;断层面上的正应力在东北方向（破裂的正方向）随着传播距离的增大而不断减小,位错速率随着破裂的传播距离而越来越大,其脉冲变得越来越尖锐,即产生了Weertman脉冲.研究结果显示：由于这种脉冲的出现,破裂在正方向上（东北方向）能够自己放大、自己愈合、自行维持,摩擦热极小,所以破裂能够沿着东北方向一直传播,直到应力场方位发生变化,不利于破裂时才最后终止.但在西南方向,破裂过程中断层面上的正应力增大,阻碍破裂继续扩展.最后就出现了汶川地震中破裂朝东北方向单侧优势传播的基本格局.模拟结果还表明：若断层面两侧介质均匀,则破裂向两侧是对称传播,且破裂距离很短,因此这种情况无法产生像汶川大地震那样的特大地震.因此,文中的模拟结果表明龙门山断裂带两侧的物性差异是造成汶川大地震单侧传播的决定性因素.断层两侧物性差异（bimaterial contrast）影响断层破裂过程的研究对于深入认识地震动力学过程、地震灾害预测及评估等有重要的科学意义.     

为什么自然界中超剪切破裂的地震是如此之少?

超剪切破裂的地震由于比同震级的亚剪切地震造成更为严重的地震灾害,因而受到广泛重视.可是,自然界中超剪切破裂的地震数量极其稀少;到目前为止,全球只是发现了十几个超剪切破裂的震例.那么是哪些因素造成了超剪切破裂地震是如此之少?这个科学问题,至今未见前人有系统的分析和研究;为此,本文利用有限单元的计算方法,对触发产生超剪切破裂的几种常见力学因素进行定量分析.有限元模拟结果显示,地表作为触发产生超剪切破裂效率最高的因素,但当地表附近具有沉积层、未固结的断层泥等松软介质时,地表附近断层介质是由速度强化的摩擦本构关系所控制,此时亚剪切破裂无法转换为超剪切破裂,近地表的这层薄薄的摩擦速度强化层会有效抑制超剪切破裂的发生.此外,模拟结果还表明,断层上的障碍体、反凹凸体、断层阶区虽然可以促使亚剪切破裂转换为超剪切破裂,但在转化过程中,由于破裂出现时间上的停顿,这样整个断层上的破裂速度就被平均为亚剪切破裂,即远场地震台站接收的地震信息很可能无法辨别局部超剪切破裂的发生.若要产生一个目前可以识别的超剪切破裂地震事件,其条件十分苛刻,通常需要断层的长度足够长,断层几何要足够平直,初始应力水平要足够高,地表的摩擦速度强化层厚度要足够薄,近场要有足够多的地震台站,等.由此可见,上述这些因素导致了自然界中可以辨别的超剪切破裂地震非常稀少.本研究有助于我们深入理解超剪切破裂发生的力学机制,有助于我们更好地评估地震灾害.     

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

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

断层自发破裂动力过程的有限单元法模拟

断层自发破裂动力过程的研究对于认识地震过程及减轻地震灾害有着重要的科学意义.为合理地模拟断层的自发破裂过程,本文首先对经典的滑移弱化摩擦关系进行了改进,然后利用有限单元方法对破裂过程进行动态数值模拟.模拟结果表明,利用改进后的摩擦关系能够产生脉冲型(pulse-like)破裂模式,而经典的滑移弱化摩擦关系不能产生这种破裂形态.模拟结果还显示,断层自发破裂过程受初始应力场及摩擦关系影响,当初始应力场中剪应力水平较低时,容易产生脉冲型破裂;但当初始剪应力较高时,会产生裂纹型(crack-like)破裂.这个现象与在实验室里进行的岩石破裂实验结果是一致的.在相同的初始应力情况下,若滑移弱化摩擦本构关系中的动摩擦系数较大,断层将易于产生脉冲型破裂;若动摩擦系数较小,将倾向于产生裂纹型破裂.此外,本文也采用速率弱化摩擦关系对断层自发破裂过程进行了模拟,结果发现,在初始场及其他条件不变时,如果摩擦关系中的b-a值较小,容易产生脉冲型破裂;如果b-a值较大,会产生裂纹型破裂.     

断层自发破裂动力过程的有限单元法模拟

断层自发破裂动力过程的研究对于认识地震过程及减轻地震灾害有着重要的科学意义.为合理地模拟断层的自发破裂过程,本文首先对经典的滑移弱化摩擦关系进行了改进,然后利用有限单元方法对破裂过程进行动态数值模拟.模拟结果表明,利用改进后的摩擦关系能够产生脉冲型（pulse-like）破裂模式,而经典的滑移弱化摩擦关系不能产生这种破裂形态.模拟结果还显示,断层自发破裂过程受初始应力场及摩擦关系影响,当初始应力场中剪应力水平较低时,容易产生脉冲型破裂;但当初始剪应力较高时,会产生裂纹型（crack-like）破裂.这个现象与在实验室里进行的岩石破裂实验结果是一致的.在相同的初始应力情况下,若滑移弱化摩擦本构关系中的动摩擦系数较大,断层将易于产生脉冲型破裂;若动摩擦系数较小,将倾向于产生裂纹型破裂.此外,本文也采用速率弱化摩擦关系对断层自发破裂过程进行了模拟,结果发现,在初始场及其他条件不变时,如果摩擦关系中的b-a值较小,容易产生脉冲型破裂;如果b-a值较大,会产生裂纹型破裂.     

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

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

     

断层阶区对震源破裂传播过程的控制作用研究

地震破裂能否穿越断层阶区（stepover）引发更大震级的地震是震源动力学研究的重要内容.本文利用不连续变形体接触力学的动态有限单元方法,模拟断层阶区对地震破裂传播的控制作用.通过改变断层周边初始应力场、断层面上的摩擦本构关系以及断层阶区的间距大小来分析各个因素对破裂传播的影响,并定量分析产生这些影响的力学机制.模拟结果表明：断层面上的摩擦系数减小或断层周边区域内的初始剪应力增大,都将增加断层破裂跳跃阶区传播的可能性;此外,若断层阶区间距越小,断层破裂也越容易跳跃阶区传播.计算结果还显示：断层上的摩擦系数大、初始剪应力小、断层阶区间隔大,那么此阶区所在之处将可能是断层破裂的终止位置;相反,当断层面上的摩擦系数较小、初始剪应力较大、断层阶区间隔较小,破裂就容易穿越阶区而出现较大的地震.同时,从模拟结果可以看出,在发震断层破裂停止后,应力将继续向四周传播;当应力积累达到破裂极限时,触发断层阶区中的另一断层产生破裂,因此在破裂跳跃断层阶区的过程中存在一个时间延迟.最后,破裂能否跳跃断层阶区,可以利用库仑应力在空间的分布进行合理的解释.     

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

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

实际地形对2010年玉树MS7.1地震动力学破裂过程及地震动特征的影响

基于震源动力学模型,采用曲线网格有限差分法模拟了水平自由地表及实际地形条件下2010年玉树MS7.1地震的震源破裂过程及相应的地面地震动,获取了断层面上的破裂传播过程、最终滑移量分布及峰值地面速度分布,讨论了实际地形对玉树地震破裂过程及相应地震动的影响.基于本文设定的动力学模型,模拟结果显示:断层面上的高应力降是玉树地...     

新生破裂影响地震破裂跨越断层阶区传播过程的扩展有限单元法模拟

采用扩展有限元方法计算了断层阶区内介质产生的新生破裂对地震破裂跨越断层阶区传播过程的影响。模型中新生的断层扩展遵循最大剪应力破坏准则，当最大剪应力超过岩石的承受极限时，完整介质产生破裂形成新的断层，并且新断层的扩展方向为最大剪应力方向。扩展有限元法模拟结果表明，断层阶区内新生的断层改变了断层阶区的几何形态，同时也改变了断层破裂后的应力状态。新生破裂可以改变库仑应力在空间的分布格局，特别是可以提高断层上的应力水平，从而提高地震破裂跨越断层阶区的能力。模拟结果还显示，断层阶区内新生破裂的产生，可以使得地震破裂跨越10 km宽的断层阶区，若阶区内部介质没有产生新生破裂，则地震破裂无法跨越该断层阶区。本研究有助于进一步认识地震破裂跨越断层阶区的传播过程，特别是对地震震源过程分析及地震灾害评估等具有重要的科学意义。     

断层自发破裂传播的数值模拟:速率状态摩擦定律本构参数的影响

地震往往受控于滑动面的摩擦性质,这种摩擦性质可以由速率状态摩擦定律较好地描述.速率状态摩擦定律中的本构参数a和b与动态摩擦系数相关,从而影响着同震位移与剪切应力的时空演化.本文在前人工作的基础上,采用三维边界积分方程法模拟速率状态摩擦定律控制下均匀全空间中平面断层的自发破裂传播过程,并详细讨论了a和b对滑动速率、剪切应力和破裂传播速度的影响.数值结果表明a和b的不同取值将导致不同的破裂行为,b－a的值越大,断层越不稳定,这种不稳定性有利于裂纹的产生与扩展.但滑动速率的时空分布不只依赖b－a,而且还与a和b的具体取值有关,断层面上滑动速率峰值与剪切破裂强度均随着a的减小而增大,随着b的增大而增大.相关结果有助于加深对断层自发破裂传播的认识.     

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

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

汶川Ms8.0地震动力学背景的数值模拟实验

2008年5月12日汶川Ms8.0地震发生在龙门山断裂带.地球物理研究者利用多种方法对此次地震的动力学成因及过程进行了研究,本文从研究断层错动行为的角度出发来探索影响汶川地震断层初始破裂及破裂的传播过程的动力学因素.     

断层破裂传播速度与破坏性地震的烈度分布——浅源地震的断层不均匀性与短周期波动的多普勒效应

地震的等震线不仅反映了地震的大小,而且也反映了地震断层过程的类型和破裂速度.单侧破裂的断层过程呈现出蛋型的等震线,而双侧破裂的断层过程呈现出椭圆的等震线,等震线的长轴与短轴之比对于破裂传播速度十分敏感.本研究利用理论等震线与观测结果相比较的方法,确定了1964年日本新地震、1983年日本海中部地震、1975年中国海城地震和1976年中国唐山地震的破裂传播速度、断层类型以及断层走向.得到的断层破裂速度是剪切波的0.7—0.9倍.这些值比用长周期地震波所确定的要稍大.产生这种差异的原因是:对烈度起主要影响的短周期地震波强烈地依赖于小尺度的断层不均匀破裂过程,以及局部的破裂传播的开始和终止;而由长周期波所得到的破裂速度却反映了在整个断层上破裂传播的平均过程.根据等震线图所得出的断层类型及断层走向与其它独立方法的结果相一致.这意味着本方法可以应用于推断某些历史地震的断层类型、破裂传播速度及破裂传播方向.     

断层破裂传播速度与破坏性地震的烈度分布——浅源地震的断层不均匀性与短周期波动的多普勒效应

地震的等震线不仅反映了地震的大小,而且也反映了地震断层过程的类型和破裂速度.单侧破裂的断层过程呈现出蛋型的等震线,而双侧破裂的断层过程呈现出椭圆的等震线,等震线的长轴与短轴之比对于破裂传播速度十分敏感.本研究利用理论等震线与观测结果相比较的方法,确定了1964年日本新地震、1983年日本海中部地震、1975年中国海城地震和1976年中国唐山地震的破裂传播速度、断层类型以及断层走向.得到的断层破裂速度是剪切波的0.7——0.9倍.这些值比用长周期地震波所确定的要稍大.产生这种差异的原因是:对烈度起主要影响的短周期地震波强烈地依赖于小尺度的断层不均匀破裂过程,以及局部的破裂传播的开始和终止;而由长周期波所得到的破裂速度却反映了在整个断层上破裂传播的平均过程.根据等震线图所得出的断层类型及断层走向与其它独立方法的结果相一致.这意味着本方法可以应用于推断某些历史地震的断层类型、破裂传播速度及破裂传播方向.      

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

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