1976年M_S7.8唐山地震断层动态破裂及近断层强地面运动特征

采用美国南加州地震委员会(SCEC)Steven Day博士提供的三维有限差分断层瞬态破裂动力学模型(3D-FDM),以1976年唐山M_S7.8地震为例,从简化的断层双侧破裂模式出发,对该地震发震断层的动态破裂过程及近断层地表运动特征进行了仿真模拟和计算.研究区域为围绕发震断层200 km×140 km×40 km(深度)的长方形块体组成,模拟计算的空间分辨率和时间分辨率分别为200 m和0.012 s,形成的空间网格节点数为1051×701×201.在DELL小型上作站上,我们实现了对源程序的移植和并行计算.同时,通过引进计算机可视化技术,对模拟数据进行了3D/4D解释分析.另外,在对源程序修改过程中,实现了对京津唐地区三维地壳速度结构的嵌入,在一定程度上增强了对地震波传播以及地面运动模拟的真实性,并讨论了地震破裂的方向性对近断层地表运动的影响.最后根据初步研究结果结合京津唐地区活动断层构造特征,对唐山M_S7.8级主震后随之而来的1976滦县M_S7.1级余震及宁河M_S6.9级余震的动态触发机制提出了新的解释.由于受主震破裂方向性作用的影响,使得主震对后续两个较大余震产生的动态应力变化的峰值在断层的走滑方向上较大,为2～3 MPa,在逆冲方向上较小,为0.1～0.2 MPa.即唐山主震的发生使得其周边的应力场有一个瞬态的应力调整,唐山主震对后续余震的发生有促发作用.     

1976年MS7.8唐山地震断层动态破裂及近断层强地面运动特征

采用美国南加州地震委员会（SCEC）Steven Day博士提供的三维有限差分断层瞬态破裂动力学模型（3D-FDM）,以1976年唐山M_S7.8地震为例,从简化的断层双侧破裂模式出发,对该地震发震断层的动态破裂过程及近断层地表运动特征进行了仿真模拟和计算.研究区域为围绕发震断层200 km×140 km×40 km（深度）的长方形块体组成,模拟计算的空间分辨率和时间分辨率分别为200 m和0.012 s,形成的空间网格节点数为1051×701×201.在DELL小型工作站上,我们实现了对源程序的移植和并行计算.同时,通过引进计算机可视化技术,对模拟数据进行了3D/4D解释分析.另外,在对源程序修改过程中,实现了对京津唐地区三维地壳速度结构的嵌入,在一定程度上增强了对地震波传播以及地面运动模拟的真实性,并讨论了地震破裂的方向性对近断层地表运动的影响.最后根据初步研究结果结合京津唐地区活动断层构造特征,对唐山M_S7.8级主震后随之而来的1976滦县M_S7.1级余震及宁河M_S6.9级余震的动态触发机制提出了新的解释.由于受主震破裂方向性作用的影响,使得主震对后续两个较大余震产生的动态应力变化的峰值在断层的走滑方向上较大,为2~3 MPa,在逆冲方向上较小,为0.1~0.2 MPa.即唐山主震的发生使得其周边的应力场有一个瞬态的应力调整,唐山主震对后续余震的发生有促发作用.     

Rupture directivity effect on strong ground motion during the 12 May 2008 MW7.9 Wenchuan earthquake

We focus here on the rupture directivity effect on the spatial distribution and attenuation characteristics of near-field ground motions during the 2008 MW7.9 Wenchuan earthquake. We examine the difference between the observed ground motions in and opposite the rupture directions and compare them with Next Generation Attenuation-West2 (NGA-West2) ground motion prediction models. The isochrone directivity predictor is used to quantify the band-limited nature of the rupture directivity effect on strong ground motion. Our results show that the observed peak ground velocity (PGV) and spectral accelerations of periods of 1.0 s and longer are significantly amplified in the rupture direction, but de-amplified in the opposite direction affected by rupture directivity effect of this event. In contrast, the effect of rupture directivity on the observed peak ground acceleration (PGA) and periods of shorter than 1.0 s are relatively weak. The rupture directivity of this event shows clear period dependent and band limited characteristics with the strongest effect occurring around the period of 7.5 s.     

Development of attenuation relation for the near fault ground motion from the characteristic earthquake

A composite source model has been used to simulate a broadband strong ground motion with an associated fault rupture process. A scenario earthquake fault model has been used to generate 1 000 earthquake events with a magni-tude of Mw8.0. The simulated results show that, for the characteristic event with a strike-slip faulting, the character istics of near fault ground motion is strongly dependent on the rupture directivity. If the distance between the sites and fault was given, the ground motion in the forward direction (Site A) is much larger than that in the backward direction (Site C) and that close to the fault (Site B). The SH waves radiated from the fault, which corresponds to the fault-normal component plays a key role in the ground motion amplification. Corresponding to the sites A, B, and C, the statistical analysis shows that the ratio of their aPG is 2.15:1.5:1 and their standard deviations are about 0.12, 0.11, and 0.13, respectively. If these results are applied in the current probabilistic seismic hazard analysis (PSHA), then, for the lower annual frequency of exceedance of peak ground acceleration, the predicted aPG from the hazard curve could reduce by 30% or more compared with the current PSHA model used in the developing of seismic hazard map in the USA. Therefore, with a consideration of near fault ground motion caused by the rupture directivity, the regression model used in the development of the regional attenuation relation should be modified accordingly.     

2010年玉树超剪切破裂地震破裂过程反演

为探讨玉树地震是否为超剪切破裂事件,利用Yagi的方法重新进行了此次地震震源破裂过程的反演。通过给定不同的破裂速度进行对比,发现当破裂速度为4.7km/s时,理论与实际观测结果拟合残差最小,且反演结果更符合实际。另据P波高频辐射能量包络反演,破裂传播速度在4.7～5.8km/s范围内,而该地区的剪切波速度则为3.0～3.6km/s,从而证明了此次地震超剪切破裂现象的存在。反演结果表明,此次地震在玉树段的NW和SE段形成了19和31km的地表破裂,而震中所在的中段由于隆宝湖拉分盆地的存在,造成了15km的未破裂区。超剪切破裂是造成SE段玉树县城遭受严重破坏的原因之一。     

2015年尼泊尔Gorkha地震强地面运动记录分析

2015年4月25日在尼泊尔Gorkha地区发生M_W7.8地震,距离发震断层约11 km的KATNP台站完整记录了主震的加速度时程.本文根据KATNP台站记录的加速度数据分析了Gorkha地震的地震动特征.结果表明Gorkha地震在KATNP台站处产生的水平向峰值加速度为0.17 g,竖直向峰值加速度为0.19 g,该数值小于科学家们对如此大规模地震产生的地震动的预期,初步推测这可能是由加德满都山谷产生的非线性响应造成的（Dixit et al.,2015）;地震在KATNP台站处产生了地表永久位移,其中竖向永久位移为131.9 cm,水平向永久位移的绝对值为159.2 cm,方向为南偏西19°（199°）,据此可简单推算出断层走向约为289°（109°）.地震产生了脉冲型地震动,影响因素有盆地效应、地震破裂的向前的方向性效应以及滑冲效应,其中盆地效应的周期约为5 s左右,方向性效应产生的速度脉冲的周期约为8 s左右.加速度反应谱显示在0.5 s和5.0 s左右各有一个峰值,前者是由地震破裂的脉冲式滑移产生的大量高频地震动造成的,后者是由于盆地效应和地震破裂的方向性效应造成的.基于阿里亚斯烈度计算的地震动持时约在36~46 s之间,小于与其规模相当的地震产生的地震动持时,并且不同方向上的地震动持时可能与地震破裂方向有关.阿里亚斯烈度随时间的变化比较简单,也反映了Gorkha地震是一次连续的、能量释放相对简单的地震事件.     

特征地震强地面运动参数衰减公式拟合

计算机仿真模拟设定地震断层动态破裂传播和近断层强地表运动响应的结果表明， 对于特征地震而言，近断层附近的地表运动特征与断层破裂传播的方向性有着强烈的依赖关系. 当场地（观测点）至断层的距离给定时，正向于破裂传播方向的场地（场地A）的地表质点运动（位移、速度、加速度），远远大于震中附近（场地B）和反向于破裂传播方向的场地（场地C）的地表质点运动，而且沿断层垂直分量所辐射的SH波的传播起到了主导作用. 对应于场地A，B和C，统计分析结果表明，峰值加速度的几何平均值之比为2.15:1.5:1， 而且各自的均方差分别为0.12, 0.11和0.13. 如果将所得的研究结果应用于概率地震危险性分析中，对于较低的年超越频度，近断层附近的地表峰值加速度的估算值可下降15%~30%. 因此，考虑到断层破裂传播方向性对地表运动的影响，区域衰减曲线的回归分析模型应该给予恰当的修正.      

震源破裂方式和断层性质对近场强地震动特征的影响

基于数值格林函数法的近场强地震动数值模拟方法,以1994年Northridge地震断层面上位错量的不均匀分布模式和该地区的地层剪切波速度结构为震源模型和计算模型,做了两个方面的模拟研究:(1)直立走滑断层(断层倾角为90°)情况下,模拟分析了有限断层单侧破裂模式和双侧破裂模式对强地震动特征——破裂方向性和上盘效应的影响;(2)对于倾斜断层(倾角为45°),模拟分析了正断层、逆断层和走滑断层情况下,单侧破裂模式对其强地震动主要特征——破裂方向性效应和上盘效应的影响.结果表明:断层的破裂方式直接影响着地表地震动峰值和矢量分布;在近场区无论直立断层还是倾斜断层,其地表地震动峰值分布所表达的破裂方向性效应显著,位于破裂传播前方的地震动强度大,反映了波前被压缩的趋势,破裂后方地震强度明显变小;倾斜断层引起的上盘效应明显,NS向分量和竖向分量的地表地震动峰值的最大值出现在上盘靠近断层迹线处,EW向分量的峰值在断层迹线两侧呈不对称分布,且逆断层引起的地震动峰值最大,走滑断层的次之,正断层的最小.     

Distant‐earthquake simulations considering source rupture propagation: refining the seismic hazard of Hong Kong

Deterministic and probabilistic seismic hazard analyses should be complementary, in the sense that probabilistic analysis may be used to identify the controlling deterministic design‐level earthquake events, and more sophisticated models of these events may then be developed to account for uncertainties that could not have been included directly in the probabilistic analysis. De‐aggregation of the tentative uniform hazard spectra (UHS) in Hong Kong resulting from a probabilistic seismic hazard assessment (PSHA) indicates that strong and major distant earthquakes, rather than moderate local earthquakes, make the largest contribution to the seismic hazard level within the natural‐period range longer than 0.3 s. Ground‐motion simulations of controlling events located 90 and 340 km from Hong Kong, taking into account uncertainties in the rupture process, reveal that the tentative UHS resulting from the PSHA may have significantly underestimated the mid‐to‐long period components. This is attributed mainly to the adoption of double‐corner source‐spectrum models in the attenuation relationships employed in the PSHA. The results of the simulations indicate clearly that rupture directivity and rupture velocity can significantly affect the characteristics of ground motions, even from such distant earthquakes. The rupture‐directivity effects have profound implications in elongating the second corner period where the constant velocity intersects the constant displacement, thus increasing the associated displacement demand. However, demands for acceleration and velocity are found to be not sensitive to the presence of the directivity pulses. Ground pulses resulting from forward rupture directivity of distant earthquakes have longer predominant periods than the usual near‐fault directivity pulses. These long‐period pulses may have profound implications for metropolises, such as Hong Kong and others in Southeast Asia, having large concentration of high‐rise buildings. Copyright © 2005 John Wiley & Sons, Ltd.     

Response of buildings to near-field pulse-like ground motions

Ground motions affected by directivity focusing at near-field stations contain distinct pulses in acceleration, velocity, and displacement histories. For the same Peak Ground Acceleration (PGA) and duration of shaking, ground motions with directivity pulses can generate much higher base shears, inter-storey drifts, and roof displacements in high-rise buildings as compared to the 1940 El Centro ground motion which does not contain these pulses. Also, the ductility demand can be much higher and the effectiveness of supplemental damping lower for pulse-like ground motions. This paper presents a simple interpretation of the response characteristics of three recorded and one synthetic near-field ground motions. It is seen that for pulse-like ground motions—similar to any other ground motion—the Peak values of Ground Acceleration, Velocity, and Displacement (PGA, PGV and PGD) are the key response parameters. Near-field ground motions with directivity effects tend to have high PGV/PGA ratio, which dramatically influences their response characteristics. Copyright © 1999 John Wiley & Sons, Ltd.     

Attenuation of peak ground accelerations from the great Wenchuan earthquake

Over 800 accelerograms recorded by 272 ground-level stations during the Wenchuan earthquake are used to analyze the influence of rupture distance, local site conditions and azimuth on peak ground accelerations (PGAs). To achieve a better understanding of the characteristics of ground motions, the spatial distributions of the EW, NS and UD components of PGAs arc obtained. Comparisons between the EW and NS components, the fault-normal and fault-parallel components,and the vertical and horizontal components of PGAs are performed, and the regression formula of the vertical-to-horizontal ratio of PGAs is developed. The attenuation relationship of peak horizontal accelerations (PHAs) is compared with several contemporary attenuation relationships. In addition, an analysis of residuals is conducted to identify the potential effects of rupture distance, azimuth and site conditions on the observed values of PHAs. The analysis focuses on medium-hard soil site conditions, as they provided most of the data used in this study.     

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

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

近断层地震动的基本特征

本文对近断层地震动的基本特征作了详细分析,这些基本特征主要包括近断层强地震动的集中性、地表破裂、地面永久变形、破裂的方向性效应、近断层速度大脉冲和上盘效应。这些特征虽然在一次地震中不一定同时出现,但它们是已经被强震观测资料和数值模拟证实了的,在模拟和预测近断层地震动时,必须充分考虑这些特点,合理的近断层地震动模拟或预测模型和方法应当在结果中再现这些特点。     

A Condition for Super-shear Rupture Propagation in a Heterogeneous Stress Field

— We have used numerical simulations with the boundary integral equation method to investigate a mechanism to excite super-shear rupture velocities in a homogeneous stress field including an asperity of increased initial stress. When the rupture, with the slip-weakening distance selected to generate sub-Rayleigh speed, encounters the asperity it either accelerates to super-shear velocities or maintains the sub-Rayleigh speed, dependent on the size and amplitude of the asperity. Three classes of rupture propagation are identified: the velocity (a) for the most narrow asperities increases slowly towards the Rayleigh wave speed, (b) for intermediate width of the asperities jumps to super-shear values for a short distance but then decreases to sub-Rayleigh wave speeds, and (c) for the widest asperities jumps to super-shear values and pertains to values between the S- and P-wave velocities. The transitions between the three classes of rupture propagation are characterized by very narrow (critical) ranges of rupture resistance. If the size of the initial asperity is smaller than critical, it becomes difficult for rupture to propagate with super-shear velocities even if the initial stress level is high. Our results suggest that stress variation along the rupture path helps homogenize the rupture velocity and propagate with sub-Rayleigh wave speeds.     

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

基于震源动力学模型,采用曲线网格有限差分法模拟了水平自由地表及实际地形条件下2010年玉树M_S7.1地震的震源破裂过程及相应的地面地震动,获取了断层面上的破裂传播过程、最终滑移量分布及峰值地面速度分布,讨论了实际地形对玉树地震破裂过程及相应地震动的影响。基于本文设定的动力学模型,模拟结果显示:断层面上的高应力降是玉树地震出现超剪切破裂传播现象的主要原因;计算区域的实际地形阻碍了由自由地表引起的超剪切破裂的产生,对断层面滑移量的分布特征及滑动速率影响较大,进而在一定程度上降低了地震动峰值,同时对地震动的分布特征产生影响,且地震动平行断层面的水平分量相对受影响更大。      

台湾双冬断层近场脉冲型地震动的数值模拟

根据我国台湾地区西部的地质地貌特征和1999年集集M_W7.6地震的震源参数,建立了三维速度结构模型和两类震源模型。基于地壳中断层的位错积累量和岩石破裂后应力应变的传播特性,采用三维有限差分法对双冬断层活动可能产生的近场脉冲型地震动进行了模拟研究。结果表明:走滑断层垂直于断层走向的水平分量和逆断层垂直分量的峰值速度较大;由方向性效应所产生的双向速度脉冲主要集中在垂直于断层滑动分量方向,而由滑冲效应所产生的单向速度脉冲则主要集中在平行于断层滑动分量的方向;受方向性效应和上盘效应的共同制约,近场脉冲型地震动呈不对称带状分布,速度脉冲多分布在距离走滑断层迹线15 km和逆断层迹线10 km的范围内;速度反应谱在断层面的覆盖范围内沿破裂方向逐渐增大,且速度脉冲可能会对大型建筑物产生严重的剪切破坏。受凹凸体特性的影响,地震波场显示南投、台中和苗栗处于强地震动危险区。      

Numerical simulation of strong ground motion for the M_s8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground mo-tion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the north-western margin of the Sichuan Basin and caused by both the directivity of fault rupture and the ampli-fication in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard as-sessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008.     

2018年2月6日花莲M_W6.4地震近场地震动方向性效应

利用中国台湾省内222个强震动台站以及Palert地震预警系统520个台站所观测的三分量加速度记录,研究此次花莲M_W6.4地震近场强地震动空间分布和衰减特征,将观测结果与美国NGA-West2地震动经验预测模型进行对比,揭示此次台湾花莲地震近场地震动的长周期特点,基于回归残差分析研究地震动峰值加速度(PGA)、峰值速度(PGV)和不同周期地震动的空间分布差异,定量考察近场地震动的方向性效应.研究结果表明:(1)整体上此次地震的近场PGV观测值和周期1.0s以上的长周期加速度谱值与美国NGA-West2地震动预测模型结果接近,PGA观测值和周期小于1.0s的加速度反应谱略低于预测模型结果.从空间分布来看,周期1.0s以上的长周期地震动在断层的不同方位有系统性差异,在破裂传播前方(震中西南方位),周期大于1.0s时的反应谱明显高于美国NGA-West2地震动经验预测模型,在破裂传播后方(震中东北方位),周期大于1.0s时的反应谱低于经验预测模型,表明此次地震近场地震动具有显著的方向性效应.(2)破裂传播的方向性效应主要影响周期超过1.0s的长周期,而对PGA以及周期小于1.0s的短周期地震动影响较弱.在破裂传播前方,周期1.0～10.0s的加速度反应谱值被增强到整体观测平均水平的1.16～1.52倍;在破裂传播后方,周期1.0～10.0s的加速度反应谱值被减弱到整体观测平均水平的0.36～0.70倍.(3)此次地震破裂方向性效应的影响表现出明显的窄带效应,破裂方向性的影响(包括破裂传播前方的增强作用和破裂传播后方的减弱作用)在周期T=3.0s时达到最大,在该周期破裂传播前方的增强系数为1.52,破裂传播后方的减弱系数为0.36.从周期T=3.0s到10.0s,破裂方向性效应的影响随周期增大总体上呈减弱趋势,这与2016年日本熊本M_W7.0地震破裂方向性效应的影响特点显著不同.     

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

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

2016年2月6日美浓6.4级地震P-alert台网强震动记录初步分析

P-alert台网实时数据对地震预警及烈度速报和工程地震研究都是重要的补充,处理分析这些数据对客观衡量P-alert台网数据质量和数据用途有重要意义。对2016年2月6日台湾美浓ML6.4地震P-alert台网获取的记录进行了处理和初步分析,统计分析显示100gal以上加速度记录有112条,200gal以上加速度记录有33条,400gal以上加速度记录有7条,东西向最大峰值加速度为466.4gal,南北向最大峰值加速度为498.4gal,竖直向最大峰值加速度为258.6gal,最大仪器地震烈度为9.5度。竖直向峰值加速度和峰值速度比水平向峰值加速度和峰值速度衰减快。峰值加速度比峰值速度衰减快,观测峰值加速度和峰值速度与台湾西南地区峰值加速度和峰值速度衰减公式比较一致。计算得到了近场台站的永久位移,显示P-alert台网绝大多数永久位移在1cm到5cm之间,最大永久位移达8cm。