Study on simulating strong ground motion by fractal stochastic method

The time history of strong ground motion can be synthesized by empirical Green's function (EGF) method.Firstly a large seismic event is discretized into a series of subevents; secondly recordings of earthquakes with proper size and spatial distribution are chosen as time history (EGF) of those subevents; finally the EGFs are summated to get the time history of ground motion caused by the large event.     

Seismic ground motion relationship in southern China based on stochastic finite-fault model

The characteristics of seismic ground motions in southern China are difficult to determine statistically due to a lack of strong ground motion data. In this study, a stochastic finite-fault ground motion model was adopted to simulate the seismic ground motions at bedrock for southern China, based on parameters derived from small and medium earthquakes that have occurred in the region. From these, the response spectra was estimated. A set of ground motion attenuation relations hipswas then developed based on simulated peak ground motions and response spectral parameters through regression, which would be applicable for use in engineering practice. Through comparisons, it was demonstrated that the proposed ground motion relationships are generally consistent with those obtained from other reported ground motion attenuation models for southern China.     

中强地震随机有限断层模型应力降参数的确定方法

随机有限断层模型是模拟地震动加速度时程的一个重要工具.但将其应用于中强地震时,由于震源信息的准确性较差从而使模型参数具有较大的不确定性.尤其针对其中最为关键的应力降参数,目前相关研究较为缺乏且尚未形成系统的确定方法.本文基于美国Little Skull Mountain Mw5.6级地震2个近场台站记录的地震动模拟,详细研究了采用随机有限断层法拟合中强地震地震动伪加速度反应谱（PSA）来确定应力降参数值的方法,并在计算应力降时引入了其它震源参数的不确定性,随后对此方法的可行性进行了验证.研究表明：采用不同频段反应谱残差和计算得到的应力降值差别较大,确定中强地震应力降较为合适的反应谱频段是中高频,采用该频段确定的应力降参数值模拟的反应谱和峰值加速度与实际记录较为符合;脉冲子断层百分比、断层长宽、倾角和深度等震源参数按截断的正态分布或均匀分布随机抽样,拟合得到的应力降参数值与通过实际震源模型参数得到的值相近.以上研究结果对确定一个区域中强地震应力降或中强地震近场强震动模拟研究提供了更进一步的研究方法和研究方向.     

三危山断层近场地震动对敦煌莫高窟的影响研究

三危山断层距敦煌莫高窟仅约3km。本建立横向均匀成层的地震断层模型，利用随机有限断层方法，预测三危山断层发生5．5级、6．0级、6．5级地震时地震动对莫高窟的影响。研究表明逆断层使三危山断层辐射的地震动分布形态呈现出明显的上下盘效应。同时，将计算PGA值与已有的根据其它地区地震记录推出的PGA衰减关系进行了对比分析，认为该结果可供莫高窟抗震防护加固设计时参考。     

龙门山地区的kappa（κ0）模型及汶川MS8.0地震的强地震动模拟

在地震动模拟中,高频衰减参数（κ₀）是一个重要的参数,它控制了傅里叶谱高频部分的衰减特性.本文利用汶川M_S8.0地震和芦山M_S7.0地震主余震的1671组强震动观测记录,计算了龙门山地区50个断层距小于150 km的强震动台站的κ₀,基于随机有限断层法模拟了汶川地震中这些台站的加速度时程、傅里叶振幅谱和反应谱,并与前人的研究进行比较.结果表明,合理计算的κ₀可以有效地改善加速度时程振幅和高频谱（>1 Hz）的模拟结果.另一方面,本文基于κ₀与地形高程的相关性,建立了龙门山地区的κ₀经验模型.分别采用该经验模型和κ₀= 0.04 s模拟了汶川地震的峰值加速度分布,并与观测记录进行比较.结果表明,采用本文提出的κ₀模型可以更好地重现汶川地震的峰值加速度分布,特别是在断层破裂方向的反方向区域和山区.综上,汶川地震中山区的峰值加速度明显大于盆地地区的现象,不仅与断层破裂产状有关,还与山区和盆地地区的κ₀之间显著的差异有关.     

Calibration of stochastic finite-fault ground motion simulations for the 1997 Umbria-Marche, Central Italy, earthquake sequence

The recent 1997 Umbria-Marche, Central Italy, earthquake sequence allowed us to model recorded ground motions using a method developed by Beresnev and Atkinson [Bull Seism Soc Am 87 (1997) 67–84; Seism Res Lett, 69 (1998) 27–32; Bull Seism Soc Am 88 (1998) 1392–1401]. The method generalizes the stochastic ground-motion simulation technique, developed for point sources, to the case of finite faults. It subdivides the fault plane into subfaults and assumes each subfault to be a point source with a ω² spectrum. Geometric spreading and regional anelastic attenuation are included in the model. The data include horizontal acceleration recordings from the SSN and ENEL databases of the 1997 Umbria-Marche events on 26 September, at 00:33 GMT, with M_w=5.7, and at 09:40 GMT, with M_w=6.0; and on 14 October at 15:23 GMT, with M_w=5.6. The strong motion simulations are performed using model parameters based on the results of previous studies, and adjusting the subfault size to calibrate the simulation model against recorded ground motions. Local site response is considered to account for observed amplification effects at specific recording sites (e.g. Nocera Umbra). A good agreement is found between the simulated response spectra and the recorded data, concluding that this method reproduces the salient ground-motion characteristics at different distances and azimuths.     

震后近断层震动图的快速产出研究—以2022年1月8日青海门源地震为例

震后近断层震动图的快速产出对于政府相关部门快速确定地震影响区、评估震害及其导致的经济损失和人员伤亡、科学决策应急救援方案和措施以减轻人员伤亡和财产损失具有重要意义。以2022年门源地震为例,利用滑动分布、应力降均不同的两个震源模型（即,王卫民等反演的震源模型和本研究生成的随机滑动震源模型）以及相同的路径、场地模型和其它输入参数,使用动力学拐角频率的随机有限断层方法和本研究建立的峰值地面速度（PGV）、水平向最大谱烈度（SI）和中国仪器地震烈度（I_I）的经验模型研究了快速产出近断层震动图（包括峰值地面加速度图、PGV图、SI图、I_I图和中国地震烈度图）的实效性。结果表明:（1）上述方法和经验模型可用于震后震动图的快速产出,其实效性主要取决于震源、路径和场地模型的可靠性;（2）基于上述两种震源模型产出的地震烈度图与中国地震局发布的该次地震的烈度图在总体上具有高度一致性,均可用于确定地震影响区,但基于反演震源模型产出的地震烈度图可以给出极震区,而基于随机滑动震源模型产出的地震烈度图则需要根据其最大等震线和发震断层的位置大致估计极震区的位置。     

Strong ground motion simulation for the 2013 Lushan M_W6.6 earthquake,Sichuan, China,based on the inverted and synthetic slip models

It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. （Chin J Geophys 56：1412-1417, 2013） and Wang et al. （Chin J Geophys 56：1408-1411,2013） published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models： the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang＇s model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model （k 2 model II） with 2 times value of comer wave number of the original k-2 model （k 2 model I） are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity （MMI） show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of sli     

2010年新西兰Mw7.0主震与2011年Mw6.1余震近场强地面运动特征比较分析

2010年9月3日16时35分46秒新西兰南岛Greendale附近发生了Mw7.0地震,震源深度约10.0 km.2011年2月21日新西兰南岛又发生了Mw6.1地震,为2010年Mw7.0主震后最大的一次余震,震源深度约5.0 km,发震断层为Christchurch南约9 km一条近东西走向逆冲的隐伏断层,该地震造成Christchurch城内多处建筑物严重损毁.本文分析了2010年新西兰地震事件Mw7.0主震与Mw6.1余震强地面运动的特征.新西兰Mw6.1余震近场强地面运动整体高于Mw7.0主震.将主震和余震的强震观测记录分别与新一代衰减关系(NGA)进行对比,发现余震强震观测数据整体高于其震级对应的NGA.分别选取距离主震和余震震中最近且强震观测记录最高的两个台站(GDLC台站和HVSC台站)作为参照台站,建立动态复合震源模型(DCSM)及有限断层随机振动模型(SFFM)进行强地面运动的模拟计算,分析两种模型的模拟结果并对比二者的优势及局限,以便在未来工作中更好地通过模型计算强地面运动特征,实现区域化特征快速、实时分析及局部重点、细致分析相结合的目标.     

Stochastic strong ground motion simulations in sparsely-monitored regions: A validation and sensitivity study on the 13 March 1992 Erzincan (Turkey) earthquake

Stochastic simulations have recently become quite popular for estimating synthetic ground motion time histories. For seismically active regions that are not well-monitored or studied extensively, input parameters of the simulations should be carefully selected as the reliability of the simulation results directly depends on the accuracy of the input parameters. In the first part of this study, 13 March 1992 Erzincan (eastern Turkey) earthquake (Mw=6.6), which is recorded at only three strong ground motion stations, is simulated using the stochastic finite-fault method. The source and regional path parameters for this event are adopted from previously validated studies whereas the local site parameters are derived herein. In the second part of the paper, sensitivity of the simulation results with respect to small changes in selected input seismic parameters is investigated. The parameters for which sensitivities are computed include stress drop, crustal shear-wave quality factor and kappa operator. A change of 20% in stress drop value results in 14% change in PGA, whereas a 20% difference in the Q₀ value causes 17% change in PGA, and a 20% variation in kappa leads to 15% difference in PGA. Numerical experiments presented in this study prove that the ground motion simulations are prone to trade-off between the source, path and site filters. Hence, input models must be implemented carefully for reliable synthetic ground motions.