近场强地震动模拟中对破裂的方向性效应和上盘效应的表达.

为了模拟近场强地震动，采用了基于有限断层模型的一种随机合成方法．震源破裂面被剖分成一定数量的子源，总的地震矩分成数量更多的子震的矩．一给定子源中一次子震引起的场地地震动的傅氏谱，通过考虑点源的谱、随着距离的衰减、能量耗散及近地表效应等导出．据此幅值谱和一个随机相位谱，并与一个时程包络函数结合，合成一个子震时程．将各子源中各个子震引起的所有时程叠加，得出场地的地震动时程．叠加中，各子震时程之间的时滞，据子震发震时差和子源至场地的距离差别引起的时差确定．对一个设定地震，选用4个断层面倾角，计算了近场21个地点的地震动．结果表明，本文方法可以很好地表达破裂的方向性效应和上盘效应．为了验证方法的可靠性，对1994年美国北岭地震中3个近断层台站MCN，LV3和PCD模拟的地震动与实际记录的加速度反应谱和时程作了比较．      

近场强地震动数值模拟的简化计算方法

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性     

Broadband Kinematic Stochastic Simulation of an Earthquake Source: a Refined Procedure for Application in Seismic Hazard Studies

To carry out a realistic simulation of earthquake strong ground motion for applied studies, one needs an earthquake fault/source simulator that can integrate most relevant features of observed earthquake ruptures. A procedure of this kind is proposed that creates a broadband kinematic source model. At lower frequencies, the source is described as propagating slip pulse with locally variable velocity. The final slip is assumed to be a two-dimensional (2D) random function. At higher frequencies, radiation from the same running strip is assumed to be random and incoherent in space. The model is discretized in space as a grid of point subsources with certain time histories. At lower frequencies, a realistic shape of source spectrum is generated implicitly by simulated kinematics of slip pulse propagation. At higher frequencies, the original approach is used to generate signals with spectra that plausibly approximate the prescribed smooth far-field source spectrum. This spectrum is set on the basis of the assumedly known regional empirical spectral scaling law, and subsource moment rate time histories are conditioned so as to fit this expected spectrum. For the random function that describes final slip over the fault area, lognormal probability distribution of amplitudes is assumed, on the basis of exploratory analysis of inverted slip distributions. Similarly, random functions that describe local slip rate time histories are assumed to have lognormal distribution of envelope amplitudes. In this way one can effectively emulate expressed ??asperities?? of final slip and occasional occurrence of large spikes on near-source accelerograms. A special procedure is proposed to simulate the spatial coherence of high-frequency fault motion. This approach permits the simulation of fault motion plausibly at high spatial resolution, fulfilling the prerequisite for simulation of strong motion in the vicinity of a fault. A particular realization (sample) of a source created in a simulation run depends on several random seeds, and also on a considerable number of parameters. Their values can be selected so as to take into account expected source features; they can also be perturbed to examine the source-related component of uncertainty of strong motion. The proposed approach to earthquake source specification is well adapted to the study of deterministic seismic hazard: it may be used for simulation of individual scenario events, or suites of such events, as well as for analysis of uncertainty for expected ground motion parameters from a particular class of events. Examples are given of application of the proposed approach to strong motion simulations and related uncertainty estimation.     

近场地震动格林函数的解析法与数值法对比研究

分析了在均匀弹性介质中，分别用解析方法和三维波动有限元数值方法计算无限全空间近场地震动的格林函数的理论与方法。将震源处理分为：(1)计算单一位错点源单位脉冲引起的格林函数；(2)用有限断层模型将断面划分为若干子源并有时间延迟，叠加所有子源的格林函数。本以1994年美国北岭地震为例，用解析法和数值方法分别计算了LV3、PCD、MCN三个场地的格林函数，并进行了对比分析。     

近断层地震动场预测——以长春市为例

断层带附近地震动场分布的研究,是当前地震工程领域研究的热点问题之一。近断层地震动场的分布对在断层附近进行抗震结构设计时,不仅是提供地震动输入,也是确定建设场地避让范围的重要依据之一。以区域地震构造背景分析、目标断层活动性鉴定、地震危险性评价为基础,结合断层探测结果,利用统计经验关系等最终确定发震断层,并建立相应的震源模型。采用显式有限元和并行计算技术计算目标区域场地的长周期地震动。利用有限断层随机合成的方法,计算高频地震动。将低频和高频地震动合成为目标区域内的宽频带地震动时程。对局部特殊场地条件地区,基于场地调查和勘探的数据,利用等效线性化等方法进行一维土层的非线性反应计算,给出这些特殊场地的宽频带地震动时程。最后,根据地震动时程获得设定地震发生时,目标区域的峰值加速度分布预测图和相应的反应谱。以长春市为例预测了在设定地震发生时,近断层地震动场的分布情况。当长春尖山子—卡伦断层发生6.0级地震时,潜在破坏性地震动的影响范围集中在附近,沿断层走向分布。加速度峰值沿断层垂直变化,主要为90 Gal～140 Gal。只是在长春市南部加速度峰值达到200 Gal。本研究的预测结果具备断层附近地震动的一些最基本的特征,符合当前对断层附近地震动的基本认识。     

利用PGA快速确定汶川地震破裂特征

快速确定断层破裂特征是烈度速报的一项重要技术,断层破裂特征可为烈度速报提供震源模型,提升烈度速报准确性。通过汶川地震加速度记录,提出一种快速计算震源破裂参数的方法。 假定断层为线源模型,以一定间距将断层离散化为若干子源,以震中为不动点,通过旋转获得所有断层可能的走向,通过每次移动一个子源,获得断层所有可能的破裂方式,将二者结合即可给出断层所有可能的空间分布;计算每种断层空间分布与每个台站的断层距,利用加速度记录峰值和断层距统计回归衰减关系,分析每个衰减关系的拟合残差,残差最小拟合效果最好的衰减关系所对应断层参数,能够对该次地震的地震动场有最合理的解释,也最有可能是实际地震中的断层空间分布。     

2010年4月14日玉树Ms7.1地震加速度场预测

基于有限断层震源、且使用动力学拐角频率的地震动随机模拟方法预测玉树地震近断层的加速度场.首先,基于有限断层震源建模方法建立该次地震的震源模型;然后,基于上述地震动模拟方法预测玉树地震近断层191个节点的加速度时程.在此基础上,取每个结点的加速度峰值绘制该次地震的近断层加速度场.结果表明:(1)近断层加速度场主要受震源破裂过程和断层面上滑动分布的影响.断层面上凹凸体投影到地表的区域附近,加速度峰值最大,也是震害最严重的区域;(2)对于走滑地震,断层沿线附近的场地并非均会发生破裂方向性效应;发生破裂方向性效应的场地与凹凸体在断层面上的位置有关.     

Temporal and spatial rupture process of the great Kunlun Mountain Pass earthquake of November 14, 2001 from the GDSN long period waveform data

The temporal and spatial rupture process of the 14 November 2001 Kunlun Mountain Pass earthquake (KMPE) is obtained by inverting the high signal-to-noise-ratio P-waveform data of vertical components of 20 stations with epicentral distances less than 90°, which are of Global Digital Seismogragh Network (GDSN). The inverted results indicate that the KMPE consists of 3 sub-events. The rupture of the first sub-event initiated at the instrumental epicenter (35.97°N, 90.59°E) and then propagated both westwards and eastwards, extending 140 km westwards at the speed of 4.0 km/s and 80 km eastwards at the speed of 2.2 km/s, which appeared to be an asymmetrical bilateral rupture dominantly from east to west. This sub-event formed a 220-km-long fault. Fifty-two seconds after initiation of the first sub-event, at which time the first sub-event was not over but in its healing phase, the rupture of the second sub-event initiated 220 km west of the epicenter and propagated both westwards and eastwards, extending 50     

关于有限断层计算模型的研究——考虑位错时、空不均匀分布的滑动时间函数

断层面上某一点的滑动时间函数需要用三个参数来表达:断层的平均位错量、上升时间和破裂传播的时间延迟。本文基于有限断层模型,根据断层面上每个子源的位错量时、空不均匀分布特征,采用Brune模型系统确定了相应滑动时间函数的三个参数,构建了位错量呈时、空不均匀分布的有限断层滑动模型。该方法可以表现破裂面上断层位错量在时间和空间上分布的不均匀特征,既能够表现震源的复杂性同时又能够简化震源模型的构建过程,为考虑断层附近的强地震动数值模拟计算提供了有效的途径。     

Prediction of near-field strong ground motions for scenario earthquakes on active fault

1 Introduction Earthquake disaster investigations show that numerous strong earthquakes were caused by remobilization of active faults. Major casualties and severe damage to buildings as well as signi?cant economic losses resulted from the ground motions of strongearthquakescausedbyactivefaultslocatedbeneath urban areas. Recently, the potential hazard prediction of and its mitigation against active faults located beneath urban areas have become an important research topic for seismologists and…     

Modelling of Strong Ground Motions from 1991 Uttarkashi,India, Earthquake Using a Hybrid Technique

We present a simple and efficient hybrid technique for simulating earthquake strong ground motion. This procedure is the combination of the techniques of envelope function (Midorikawa et al. Tectonophysics 218:287–295, 1993) and composite source model (Zeng et al. Geophys Res Lett 21:725–728, 1994). The first step of the technique is based on the construction of the envelope function of the large earthquake by superposition of envelope functions for smaller earthquakes. The smaller earthquakes (sub-events) of varying sizes are distributed randomly, instead of uniform distribution of same size sub-events, on the fault plane. The accelerogram of large event is then obtained by combining the envelope function with a band-limited white noise. The low-cut frequency of the band-limited white noise is chosen to correspond to the corner frequency for the target earthquake magnitude and the high-cut to the Boore’s f _max or a desired frequency for the simulation. Below the low-cut frequency, the fall-off slope is 2 in accordance with the ω² earthquake source model. The technique requires the parameters such as fault area, orientation of the fault, hypocenter, size of the sub-events, stress drop, rupture velocity, duration, source–site distance and attenuation parameter. The fidelity of the technique has been demonstrated by successful modeling of the 1991 Uttarkashi, Himalaya earthquake (Ms 7). The acceptable locations of the sub-events on the fault plane have been determined using a genetic algorithm. The main characteristics of the simulated accelerograms, comprised of the duration of strong ground shaking, peak ground acceleration and Fourier and response spectra, are, in general, in good agreement with those observed at most of the sites. At some of the sites the simulated accelerograms differ from observed ones by a factor of 2–3. The local site geology and topography may cause such a difference, as these effects have not been considered in the present technique. The advantage of the technique lies in the fact that detailed parameters such as velocity-Q structures and empirical Green’s functions are not required or the records of the actual time history from the past earthquakes are not available. This method may find its application in preparing a wide range of scenarios based on simulation. This provides information that is complementary to the information available in probabilistic hazard maps.     

Modeling of the ground motion for the Petropavlovsk earthquake of November 24, 1971 (M = 7.6)

Based on the example of the strong earthquake of November 24, 1971, with the earthquake source near the Petropavlovsk-Kamchatskii, the new modeling technique of the strong ground motions within the broadband is tested. In this technique, the seismologically-substantiated models of the radiation source and elastic medium are used. The source is represented by an array of point radiation sources-dislocations with the random seismic moments (amplitudes) and with the random time functions. The new method of calculation of the Green functions is developed to describe the propagation of waves and residual displacements of a layered medium. The method is used for the simulation of the horizontal ground motion, recorded by a S5S-ISO instrument for the strong earthquake that took place on November 24, 1971 with the source near the Petropavlovsk-Kamchatskii (with a depth of 105 km). The position of the hypocenter, the sizes and the position of the fault, and “fault mechanism” were considered to be known a priori. By a trial and error method of the duration of the source process and only two spectral parameters, it was possible to simulate successfully the fundamental characteristics of the ground vibrations: the amplitude of acceleration, the velocity and displacement of the ground, their Fourier spectrum, the duration of the vibrations, and the response spectrum. The surprisingly high level of high-frequency radiation, probably connected with the intraplate position of the source, is a specific feature of the source.     

Simulation of orthogonal horizontal components of near‐fault ground motion for specified earthquake source and site characteristics

A procedure to generate horizontal pairs of synthetic near‐fault ground motion components for specified earthquake source and site characteristics is presented. Some near‐fault ground motions contain a forward directivity pulse; others do not, even when the conditions for such a pulse are favorable. The proposed procedure generates pulse‐like and non‐pulse‐like motions in appropriate proportions. We use our recent stochastic models of pulse‐like and non‐pulse‐like near‐fault ground motions that are formulated in terms of physically meaningful parameters. The parameters of these models are fitted to databases of recorded pulse‐like and non‐pulse‐like motions. Using these empirical “observations,” predictive relations are developed for the model parameters in terms of the earthquake source and site characteristics (type of faulting, earthquake magnitude, depth to top of rupture plane, source‐to‐site distance, site characteristics, and directivity parameters). The correlation coefficients between the model parameters are also estimated. For a given earthquake scenario, the probability of occurrence of a directivity pulse is first computed; pulse‐like and non‐pulse‐like motions are then simulated according to the predicted proportions using the empirical predictive models. The resulting time series are realistic and reproduce important features of recorded near‐fault ground motions, including the natural variability. Moreover, the statistics of their elastic response spectra agree with those of the NGA‐West2 dataset, with the additional feature of distinguishing between pulse‐like and non‐pulse‐like cases and between forward and backward directivity scenarios. The synthetic motions can be used in addition to or in place of recorded motions in performance‐based earthquake engineering, particularly when recorded motions are scarce.     

近断层强地震动场预测

以1997年4月11日新疆伽师地震（M_w6.1）为例,详细介绍了近断层强地震动场的预测方法.首先,用有限断层震源建模方法建立了该次地震的震源模型;然后,基于动力学拐角频率的地震动随机模拟方法,模拟了该次地震仅有主震加速度记录、且位于巨厚土层上的三个台站的加速度时程,并用实际地震记录进行了验证.在此基础上,基于预测的近断层77个节点的加速度时程的峰值绘制了该次地震的加速度场.结果表明,上述方法模拟的加速度时程在0.5 Hz以上的高频段是可行的、实用的;预测的近断层加速度场具有非常明显的上盘效应.地表最大加速度的范围与断层面上最大凹凸体位置相对应,说明与断层面上凹凸体相对应的地面上的建（构）筑物将会遭受到较为严重的震害.     

2013年四川芦山MS7.0地震近断层强地面运动模拟及烈度分布估计

2013年4月20日在我国四川省发生了芦山M_S7.0地震,地震给当地群众的生命财产安全带来了巨大的损失,其中最严重的破坏发生在震中附近的芦山、宝兴等地区.根据地震发生的快速反演结果,及发震断层面上滑动位移的分布情况,构建有限断层模型,对近断层区域的强地面运动进行初步计算,并基于强地面运动的模拟结果给出震区烈度分布的初步估计.模拟结果显示：模拟烈度图显示极震区的烈度在IX级左右,VI级烈度影响范围大致为16000 km²,该结果与中国地震局于4月27日给出的震区实测烈度图一致程度较高.强烈地震发生后,基于近断层强地面运动模拟计算的结果,可以给出相对合理的地震烈度分布情况估计,对震区震情判定及救灾工作具备较高的现实意义.     

Earthquake fault‐induced surface rupture—A hybrid strong ground motion simulation technique and discussion for structural design

Unique to the near‐source region of a large earthquake is the occurrence of strong impulsive ground motion and surface faulting referred to as ‘fling‐step’ motion. The objective of this study is to synthesize broad‐band time histories over a wide range of frequencies by characterizing rupture directivity and fling effects from the comprehensive strong motion database of the near‐fault Chi‐Chi event. To aid in the generation of these special types of ground motions, a hybrid modeling technique is introduced based on the stochastic finite‐fault radiation method and an efficient analytical approach to incorporate the observed low‐frequency features in the records close to the ruptured fault. The results show that the overall agreement among the developed hybrid methodology and recorded waveforms and response spectra is quite satisfying. A brief discussion on the design of infrastructures near seismic fault is also included. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of ground acceleration in the near source range: the case of 1976, Friuli earthquake (M = 6.5), northern Italy

A mixed statistical-deterministic model of earthquake rupture is developed for evaluating the strong ground motion in the near source range (receiver distance comparable to the fault length). The source parametrization is based on the k-square model and the propagation is computed by asymptotic Green's functions. The method is applied to the case of 1976, Friuli earthquake (M = 6.5) in northern Italy which occurred on a low-dip thrusting fault. Acceleration records at 29 stations are computed for 100 simulations of rupture histories. The mean value map of peak ground accelerations shows clearly a maximum to the south due to the inner geometry and directivity of the source. The variation of the estimated PGA versus the epicentral distance is strongly dependent on azimuth and is not decreasing monotonically. The comparison of these curves with those predicted by empirical acceleration–distance relationships shows discrepancies in the near source distance range. This study shows the importance of considering the complexity of the source rupture process for strong motion estimate in the near source range.     

漾濞MS6.4地震长周期强地震动的有限差分初步模拟

根据我国川滇地区的地质地貌特征、前人地质资料和2021年漾濞M_S6.4地震的震源参数初步反演结果,建立三维速度结构模型,采用三维有限差分法对漾濞地震的长周期地震动进行研究.对实际地震烈度和模拟烈度进行对比,同时在区域附近布置了36个观测点,给出其中6个观测点的三分量速度时程,并给出所有观测点阻尼比为5%的速度反应谱,结果表明:(1)模拟的速度峰值已超过25 cm/s,与实际漾濞地区附近的烈度为8度吻合;(2)受断层滑动分布影响,速度水平分量在西北方向分布范围更广,衰减速度明显慢于东南方向,并且在极震区,东北—西南方向的地震动衰减也较慢.方向性效应对峰值大小及其分布范围的影响显著;(3)观测点速度反应谱的特征周期分布在1~3 s范围内,可能会因为共振效应对区域附近的大型建筑物产生较为严重的破坏.需要对漾濞地震进行进一步研究,以期对漾濞地震的长周期地震动有更为细致的认识.     

Seismological asperities from the point of view of dynamic rupture modeling: the 2007 Mw6.6 Chuetsu-Oki,Japan, earthquake

We study the ground motion simulations based on three finite-source models for the 2007 Mw6.6 Niigata Chuetsu-oki, Japan, earthquake in order to discuss the performance of the input ground motion estimations for the near-field seismic hazard analysis. The three models include a kinematic source inverted from the regional accelerations, a dynamic source on a planar fault with three asperities inferred from the very-near-field ground motion particle motions, and another dynamic source model with conjugate fault segments. The ground motions are calculated for an available 3D geological model using a finite-difference method. For the comparison, we apply a goodness-of-fit score to the ground motion parameters at different stations, including the nearest one that is almost directly above the ruptured fault segments. The dynamic rupture models show good performance. We find that seismologically inferred earthquake asperities on a single fault plane can be expressed with two conjugate segments. The rupture transfer from one segment to another can generate a significant radiation; this could be interpreted as an asperity projected onto a single fault plane. This example illustrates the importance of the fault geometry that has to be taken into account when estimating the very-near-field ground motion.     

相交早第四纪断裂近断层地震动场预测——以枣庄市新城区为例

近断层地震动场的预测不仅可以为抗震结构设计提供更加丰富的地震动输入,也可以作为土地合理规划的参考依据。在一些地震活动相对较弱的地区,地震往往发生在两个早第四纪断裂相交的位置处。因此对于这种相交早第四纪断裂区的近断层地震动场预测具有实际的意义。本文以枣庄市新城区为例,预测了在设定地震发生时,相交早第四纪断裂近断层地震动场的分布情况。先分别计算了近断层区低频和高频的地震动,然后将低频和高频成分合成宽频带的地震动场。在计算低频地震动场时,建立了枣庄市新城区地下三维传播介质模型和震源模型。计算的过程是模拟震源发生破裂产生地震波后,地震波在地下介质中的传播过程。使用了三维显式有限元、快速拉格朗日有限元法等数值计算方法以及并行计算技术来提高计算的速度。由于地震动的频率成分是连续的,因此对于高频地震动使用了随机震动理论来进行合成。本文研究的预测结果具备断层附近地震动的一些最基本的特征,符合当前对断层附近地震动的基本认识。