历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原MS8.5大地震为例

大地震的发生会引起区域位移场和应力场发生变化,进而改变区域内及临近断层的应力状态和地震活动性.目前,研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化,同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而,历史上曾经发生过不少大地震,例如,1920年的海原M_S8.5大地震,是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料,前人在研究历史地震的影响时往往给出一些简单的断层滑动模型,将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响,基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料,以海原M_S8.5大地震为例,探讨了如何建立海原大地震断层滑动模型,并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算,发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时,进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.

     

2001年昆仑山口西MS8.1地震对周围断层的应力影响数值分析

大地震的发生往往会引起周围区域形变场和应力场变化,且对临近断层上的应力状态也有影响.2001年11月4日,昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型,建立了三维含地形高程的横向不均匀性椭球型地球有限元模型,采用等效体力方法,分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比,该模型考虑了地形、Moho面起伏和地球介质横向不均匀性;与一般的有限元数值模拟相比,该模型考虑了地球曲率和椭率,合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数,计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加,对这些地震的发生起促进作用;而造成玉树和德令哈地震发震断层上的库仑应力变化为负值,在一定程度上抑制了这些断层的地震活动性.此外,计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响,其中地形和椭率造成的同震位移场相对误差约10%.

     

2017年九寨沟MS7.0地震对周围地区的静态应力影响

基于九寨沟M_S7.0地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在周围主要活动断层上产生的库仑应力变化、在周围地区产生的应力场和位移场和同震库仑应力变化对余震的触发.结果表明：（1）九寨沟地震造成虎牙断裂中段库仑破裂应力有较大增加,已经超过0.01 MPa的阈值,虎牙断裂北段、塔藏断裂中段和岷江断裂北段北部的库仑破裂应力有较大降低,因此尤其要注意虎牙断裂中段的危险性.（2）水平面应力场在该地震震中东西两侧增加（拉张）,张应力起主要作用.在震中南北两侧降低（压缩）,压应力起主要作用.从水平主压和主张应力方向来看,均呈现出条形磁铁的磁场形态.从剖面上的应力场来看,在上盘的面膨胀区域内,大部分点的主张应力方向与地表是垂直的,在其他区域内,主张应力和主压应力均以震中为中心,向外呈辐射状.（3）从地表水平位移场来看,震中东西两侧物质朝震中位置汇聚,南北两侧物质向外流出,在震中处的最大水平位移量达43 mm.从地表垂直位移场来看,震中南北两侧出现明显的隆升,隆升最大值达56.8 mm.震中东西两侧出现明显的沉降,沉降最大值达74.5 mm.从剖面的位移场来看,九寨沟地震为左旋走滑地震,且有一定的正断成分.由分析可以推测该断层破裂在大致22~26 km的深度上就截止了.并推测下盘物质运动的动力来自震源北东东方向（四川块体）深度在6~30 km的下盘下层物质,上盘物质运动的动力来自震源北西西方向（巴颜喀拉块体）深度在0~6 km的上盘上层物质.（4）通过计算不同深度上主震对余震的触发作用可知,主震后的最大余震受到了主震的触发作用,多数其他余震也受到主震的触发作用.主震的发生促使了库仑应力增加地区余震的发生,抑制了一部分库仑应力减少地区余震的发生.

     

均匀弹性分层介质模型中的同震位移、应变以及应力

地震发生时产生的同震位移、应变以及应力变化,特别是同震应力的变化,在地震触发等问题的研究中有着重要的意义.本文进一步发展了基于均匀弹性水平层状介质,利用广义反射透射系数矩阵和离散波数计算同震位移的方法,使之可以计算相应的应变、应力以及同震库仑应力变化.可适用于多种情况,接收点可以位于地表以及地表以下,震源类型可以是剪切位错源以及拉张位错源.通过与半无限介质的解析解相比较,结果一致,验证了方法的可靠性以及计算精度,可以用于计算地震之后库仑应力变化,为判断余震分布提供参考.在计算同震位移时,使用了梯形积分与Filon积分相结合的积分方式,即提高了同震位移计算的速度,又保证了计算精度,有利于反演问题研究.     

福建沿海边缘陆域的原地应力测量研究

通过在福建沿海边缘陆域地区活动断裂带附近进行水压致裂原地应力测量, 得到了地表浅部构造应力的大小、方向和分布特征. 根据实测的地应力资料,用库仑摩擦滑动准则分析研究了断裂带的性质和活动性. 结果表明：① 沿海岸线边缘自北向南,其最大水平主压应力方向为北西向,与测区北西向断裂带走向近于平行,与地质构造、跨断层形变测量等方法反映的主压应力方向大体一致,与震源机制解反映的主压应力方向夹角20deg;左右; ② 水平主应力值随深度增加而增大,与垂直应力值的关系为SHSVSh和SHasymp;SV Sh, 应力状态有利于走滑断层和正断层活动;③ 利用库仑摩擦滑动准则,摩擦强度取0.6～1.0分析, 断层附近的现今地应力状态达到或超过产生正断层摩擦滑动的临界值,表明测区现今构造活动以正断层为主; ④ 反映出本区现今构造应力场的力源主要来自太平洋板块向西和菲律宾板块向北西对欧亚板块的水平挤压作用.     

周边地区强震对新疆境内中强地震活动的影响分析

该文基于1950年以来新疆周边强震与境内中强地震活动的统计分析结果, 解算了1976—2007年新疆周边强震活跃、 境内相对平静时段, 周边3组强震产生的库仑应力变化, 并结合其后3年新疆地区中强地震活动特征, 分析探讨了周边强震活动对境内地震活动的影响作用。 结果表明, 新疆周边这3组强震活动对新疆境内具有一定程度的减震作用。     

Stress triggering of normal aftershocks due to strike slip earthquakes in compressive regime

A few cases of occurrence of normal aftershocks after strike slip earthquakes in compressive regime have been reported in the literature. Occurrence of such aftershocks is intriguing as they occurred despite the apparent stabilizing influence of compressive plate tectonic stresses on the normal faults. To investigate the occurrence processes of such earthquakes, we calculate change in static stress on optimally oriented normal and reverse faults in the dilational and compressional step over zones, respectively, due to slip on a vertical strike slip fault under compressive regime. We find that change in static stress is much more pronounced on normal faults as compared to that on reverse faults, for all values of fault friction. Change in static stress on reverse fault is marginally positive only when the fault friction is low, whereas for normal faults it is positive for all values of fault friction, and is maximum for high fault friction. We suggest that strike slip faulting in compressive regime creates a localized tensile environment in the dilational step over zone, which causes normal faulting in that region. The aftershocks on such normal faults are considered to have occurred as an almost instantaneous response of stress transfer due to strike slip motion.     

Elastic and visco-elastic stress triggering in the South Iceland Seismic Zone due to large earthquakes since 1706

Damaging earthquakes in the South Iceland Seismic Zone (SISZ) occur fairly regularly and often as a series of events with a few days only between individual events. Tolerably reliable information on epicentre locations and mechanisms are available for 13 M ≥ 6 events between 1706 and 2000. For these events, we computed the co- and post-seismic stress fields, hereby approximating the SISZ by a mixed elastic/visco-elastic layered half-space. The horizontal shear stress and the Coulomb stress changes were analysed to detect possible trigger mechanisms, which may aid future earthquake mitigation efforts. We tested several criteria but must conclude that the start of an earthquake series in the SISZ cannot be explained by triggering through previous events. Inside an individual series, however, one may infer triggering. Our results are in contradiction with the findings in other regions of the world. The reason might be related to the fact that the SISZ is not a mature fault zone, in which old faults are re-activated if a certain stress level threshold is passed. In addition, uncertainties in the model parameters as well as the neglect of horizontal variations in the model and of possible stress transfer due to volcanic activity further complicate the evaluation of our results and need to be taken into account in future studies.     

Failure analysis of the granite for a dam foundation

This paper presents the numerical and experimental investigation of a dam foundation. The site consists of granite and especially weathered granite. Numerical and analytical analyses have been conducted to evaluate suitability of the granite for a dam foundation. Mohr–Coulomb failure criterion has been used to calculate the compressive strength of the rock mass. Finite element implementation of Mohr–Coulomb criterion has been used for deformation and settlement analysis. Packer testing and in-situ estimation of rock mass deformability using the Goodman jack have been performed. The results of the evaluations show that very small normal and differential settlements can be expected in structural regions due to dam construction, and the compacted concrete dam can safely be constructed on either the competent granite or the weathered granite. However, the dam foundation on the weathered granite has high conductivity potential problems for the reservoir.     

Spatio-temporal evolution of Coulomb stress in the Pacific slab inverted from the seismicity rate change and its tectonic interpretation in Hokkaido, Northern Japan

We estimated spatio-temporal evolution of Coulomb stress within the subducted Pacific slab in Hokkaido from the analysis of seismicity rate change. For this purpose we used earthquake catalog from the Institute of Seismology and Volcanology (ISV), Hokkaido University for the period 1993/4/1–2006/12/31 after relocating to compensate location errors due to the heterogeneous P- and S-wave structure beneath Hokkaido. We found that spatial pattern of Coulomb stress change inverted from the seismicity rate change is comparable with static change in Coulomb stress estimated from dislocation models. Our results and analyses reveal important insights on spatio-temporal pattern of deformation of the subducted Pacific slab in terms of Coulomb stress change. We found that the 2003 Tokachi Oki earthquake (Mw = 8.0) pervasively perturbed Coulomb stress in a regional scale with a significant impact to trigger the 2004 Kushiro Oki earthquake. The 2004 Kushiro Oki earthquake (Mw = 7.0) is another significant stressing event that changed the pattern of Coulomb stress in the area. We found that stressing events with magnitude smaller than 7.0 has minimal impact on Coulomb stress change in the Pacific slab. Similarly, comparatively deep focused large earthquakes could not change Coulomb stress significantly. Further the pattern of Coulomb stress change after the 2003 Tokachi Oki earthquake correlates the pattern of afterslip distribution in Hokkaido.