基于数值模拟的设定地震场地谱时程拟合方法

本文将确定性数值模拟方法与地震动预测方程相结合,提出了一种重大水电工程场址设定地震的地震动时程生成方法。该方法基于场址设定地震,首先采用地震动预测方程确定场址的场地相关反应谱;其次建立包含震源和场址的场地模型,通过确定性数值模拟方法生成场址地震动时程;最后对生成的场址地震动时程进行调整,使其反应谱与设计谱相一致,用于工程抗震分析。这一方法生成的地震动时程既考虑了震源机制、传播路径以及局部场地效应等物理背景,又与场地相关的设计地震反应谱保持一致,为重大工程抗震分析与评价提供了一种新的思路。     

A numerical approach for modeling near-fault ground motion and its application in the 1994 Northridge earthquake

An approach for simulating near-fault ground motion was presented by combining the finite fault model with a numerical algorithm, named investigated lump method presented in this paper for wave propagation. The investigated lumps are constructed from the auxiliary quadrilateral grids. The dynamic equilibrium equations of a typical investigated lump have been derived and obtained by integrating the stresses along the contour of the investigated lump. The stresses are calculated using the constitutive relations and the interpolation techniques. The investigated lump method is then implemented using the equilibrium equations of investigated lumps and the calculations of stresses alternately in time domain. The stability criterion of the algorithm has been given. Comparisons with the discrete wave-number method solutions for predicting the ground motions at the Pacoima Dam during the San Fernando earthquake show the validity of the method presented in this paper for simulating near-field ground motions. A finite fault source model has been implemented in the algorithm here. The source parameters given by Wald et al. (1996) [18] are applied to synthesize the ground motions at three stations during the 1994 Northridge earthquake. The simulating results qualitatively match to the corresponding ground motion records. The studies demonstrated that the approach presented in this paper is an effective tool for the numerical simulation of near-fault ground motion.     

Response of an arch dam to non-uniform excitation generated by a seismic wave scattering model

Non-uniform ground motions are generated based on a single record available at a site and seismic wave scattering analysis. The Chino Hills 2008 earthquake records at the Pacoima Dam site are used to indicate the accuracy of the method. Dynamic analysis of the Pacoima dam-reservoir-foundation under uniform and non-uniform ground motions is carried out using the EACD-3D2008 software, and the results are compared to recorded responses at different locations on the dam. There is good agreement between computed and recorded displacements of the dam for non-uniform excitation. For uniform excitation, the displacements are underestimated in comparison with those obtained from recorded excitation. Significant intensification of stresses, especially near the foundation, and different patterns of stress distribution are observed for non-uniform excitation in comparison with uniform excitation. For uniform excitation maximum stresses occur in the crown cantilever near the crest, but for non-uniform excitation the maximum stresses occur along the sides and near the foundation.     

Determination and application of path duration of seismic ground motions based on the K-NET data in Sagami Bay,Japan

Duration models are one of the important parameters in ground-motion simulations. This model varies in different study areas, and plays a critical role in nonlinear structural response analysis. Currently, available empirical models are being globally used in ground-motion simulations, with limited research focusing on path duration in specific regions. In this study, we collected 6,486 sets of three-component strong-motion records from 29 K-NET stations in the Sagami Bay, Japan, and its surrounding areas between January 2000 to October 2018. We extracted the effective duration of 386 pieces of ground-motion records by manually picking up the S-wave arrival time and calculating the significant duration. We then obtained the path duration model of the study area based on the empirical equation of dynamic corner frequency and source duration of [7]. Compared with the results of the available empirical models, the Fourier spectrum of the simulated ground motion from our effective duration model showed higher accuracy in the long-term range, with less fitting residuals. This path duration model was then applied to simulate two earthquakes of M_W5.4 and M_W6.2, respectively, in the region using the stochastic finite-fault method with a set of reliable source, path, and site parameters determined for the study area. The simulation results of most stations fit well with observation records in the 0–30 Hz frequency band. For the M_W5.4 earthquake, the simulated ground motions at KNG005/KNG010/SZO008 stations were relatively weak in the mid to high frequency band (1–30 Hz) because the quality factor and geometric diffusion model used in the simulation were the averages of the entire Sagami Bay region, causing a bias in the results of a few stations owing to local crustal velocity anomalies and topographic effects. For the M_W6.2 earthquake, the simulated ground motions were relatively weak at all SZO and TKY stations, mainly because of the close distance from these stations to the epicenter and the complex seismic-wave propagation paths. The analysis suggests that the differences between the simulation results of the two earthquakes were mainly related to complex geological conditions and seismic-wave propagation paths.     

基于震源-河谷波场数值模拟的坝址地震动参数确定方法

针对目前高坝坝址地震动参数确定方法的一些不足,本文尝试了一种基于地震学理论的坝址区三维地震动场生成方法.基本思路是将地震学和工程学结合,针对设定地震,建立震源-传播介质-坝址峡谷场地数值模型,通过超大规模的数值计算,模拟地震波从发震断层破裂开始到坝址场地的物理传播过程,生成坝址区的三维地震动参数.与传统基于衰减关系的地震动参数确定方法相比,这一方法可以考虑震源机制、传播介质和坝址峡谷场地效应等三大要素的影响.对于特定的坝址,可以生成符合实际地质构造、区域岩体动力特性以及坝址峡谷地质地形条件的地震动荷载分布,具有针对性,为重大高坝枢纽工程遭遇极端地震荷载作用时的抗震安全性提供了一种新的分析手段.     

用随机模拟方法研究设定地震地面运动

用随机模拟方法研究设定地震地面运动可分两步进行：（1）采用中小地震的数字观测资料确定研究区路径、场地参数;（2）把确定的参数应用到研究区的强地面运动随机模拟中. 为此,本文首先采用了一种可行的方法,分离软基岩场地效应和非弹性衰减的影响,并把二者分别确定出来. 研究区北部中软基岩场地的平均效应在2～4Hz频率范围为15倍左右;研究区的S波品质因子为QS＝278f[KG*2]0.346. 把这些参数用于研究区的场地和路径模型中,并选择单拐角频率震源谱模型,随机模拟了研究区未来中强地震可能在北天山中段可能造成的地面运动,模拟加速度时程和反应谱可以服务于本地区的地震灾害预测和建筑物可靠性验算.     

System identification of a concrete arch dam and calibration of its finite element model

A magnitude 4.3 earthquake occurred near Pacoima Dam on 13 January 2001. An accelerometer array that had been upgraded after the Northridge earthquake recorded the motion with 17 channels on the dam and the dam–foundation interface. Using this data, properties of the first two modes are found from a system identification study. Modal properties are also determined from a forced vibration experiment performed in 2002 and indicate a significantly stiffer system than is estimated from the 2001 earthquake records. The 2001 earthquake, although small, must have induced temporary nonlinearity. This has implications for structural health monitoring. The source of the nonlinear behaviour is believed to be loss of stiffness in the foundation rock. A finite element model of Pacoima Dam is constructed and calibrated to match modal properties determined from the system identification study. A dynamic simulation of the 2001 earthquake response produces computed motions that agree fairly well with the recorded ones. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of local site effects by a generalized inversion scheme using observed records of ‘Small-Titan’

A simple-extended array of observation for strong motions, named ‘Small-Titan’, was installed in Sendai City, Japan in 1997. This array system consists of 20 observation sites, which have an average of about 4 km for the station-to-station distance. Each observation station was located at various kinds of soils to effectively obtain local site effects. Since its completion, the system has obtained about 160 motion records. As is well known, ground motions are influenced by the three main factors such as source, path and local site effects. An accurate prediction of ground motions should thus be performed based on the three factors. In this context, theoretical methods used commonly lack enough information related with source and path effects. A generalized inversion scheme to separate the effects of source, propagation path and local site from S-wave spectra was applied to the records of Small-Titan obtained in the Sendai region, Japan. The objective of this paper is to estimate local site effects in Sendai City and to obtain meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. This paper also describes the applicability of the generalized inversion scheme to seismic zoning.

Site amplification factors estimated from the generalized inversion scheme vary significantly, reflecting the effects of soil characteristics at each site. Strong site-dependency was found in the site amplification factors. In particular, the local site amplification factors in the western side are smaller than those in the eastern side, across the Nagamachi–Rifu Fault that exists in the central part of the Sendai area. Results from the generalized inversion scheme are compared with those obtained by the direct spectral ratio scheme, the H/V method and the bootstrap inversion. There is a good agreement in the amplification factors between the generalized inversion scheme, the direct spectral ratio scheme and the bootstrap inversion scheme. The H/V method for estimating site amplification is found to extract the same dominant periods as the other methods. This paper especially emphasizes a good agreement in not only dominant period, but also amplification factor.

The generalized inversion scheme is a useful method in separating the effects of source, propagation path and local site and in revealing the site amplification. In addition, the generalized inversion scheme provides meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. Moreover, this paper indicates that the generalized inversion scheme is applicable to the seismic zoning toward disaster mitigation.     

Ground motion observations and simulation for local earthquakes in the Campi Flegrei volcanic area

Ground motion produced by low magnitude earthquakes can be used to predict peak values in high seismic risk areas where large earthquakes data are not available. In the present work 20 local earthquakes (M_D∈[?0.3, 2.2]) occurred in the Campi Flegrei caldera during the last decade were analyzed. We followed this strategy: empirical relations were used to calibrate synthetic modeling, accounting for the source features and wave propagation effects. Once the source and path parameters of ground motion simulation were obtained from the reference data set, we extrapolated scenarios for stronger earthquakes for which real data are not available. The procedure is structured in two steps: (1) evaluation of ground motion prediction equation for Campi Flegrei area and assessment of input parameters for the source, path and site effects in order to use the finite fault stochastic approach (EXSIM code); (2) simulation of two moderate-to-large earthquake scenarios for which only historical data or partial information are available (M_w4.2 and M_w5.4). The results show that the investigated area is characterized by high attenuation of peak amplitude and not negligible site effects. The stochastic approach has revealed a good tool to calibrate source, path and site parameters on small earthquakes and to generate large earthquake scenario. The investigated magnitude range represents a lower limit to apply the stochastic method as a calibration tool, due to the small size of involved faults (fault length around 200/300 m).     

Estimating ground motion incoherence through finite source simulation: a case study of the 1980 El-Asnam Earthquake

Spatial variability of ground motions has significant influence on dynamic response of extended structures such as bridges and tunnels. In this study, the widely used finite-source ground motion simulation approach, the so-called Empirical Green’s Function (EGF) method, is extended to synthesize seismic motions across an array of stations located at bedrock in the epicentral region of the 1980 El-Asnam region (North-West Algeria). The target event being simulated is the October 10 1980 \( M_{s} = 7.2 \) Earthquake, and the EGF is obtained from the ground motion recorded at Sogedia Factory station during the 8 November 1980 \( M_{L} = 5.6 \) aftershock. Coherency functions are then estimated from the simulated ground accelerations. A parametric study investigating the influence of shear wave velocity, earthquake magnitude, and epicentral distance is conducted by simulating ground acceleration for different scenarios using the Hybrid Green’s Function method. The main finding of the study is that finite source effects can cause significant loss in coherency at bedrock in the near-field. In the far-field, the source effect alone does not seem to produce incoherent motion, which implies that scattering and local site effects could be dominating there. Furthermore, coherency functions are found to be more sensitive to inter-station separation in the near-field than in the far-field. Increasing shear wave velocity seems to increase coherency functions, and larger earthquakes seem to produce more incoherent motion than smaller ones. The simulation method presented here produces incoherent motion mainly due to the finite source effect, while path effects are partially accounted for through the EGF, and local site effects are not considered. In this sense, the estimated coherency functions represent that of plane waves. A parametric model of plane wave coherency is calibrated and presented based on the simulation results. The results indicate that the parametric model can be used as a first approximation, and at least an upper bound of lagged coherency in the near-field region of the El-Asnam Earthquake scenario. This model could be useful in random vibration analysis or generation of spatially variable ground motion for time history analysis of lifeline structures in the study area.     

Generation of spatially nonuniform ground motion for nonlinear analysis of a concrete arch dam

An accelerometer array at Pacoima Dam with three locations along the base and abutments recorded ground motion from a magnitude 4.3 earthquake on 13 January 2001. These records present an opportunity to study spatial nonuniformity for the motion in a canyon. Topographic amplification is characterized by ratios of response spectral displacement between locations, and seismic wave travel times are studied using cross‐correlation functions to obtain delays. Results of the analysis of the 2001 earthquake records are used to generate ground motion for the 1994 Northridge earthquake to replace records that were not able to be fully digitized. The ground motion generated for the Northridge earthquake is used as input to a finite element model of Pacoima Dam. The response of the model is consistent with observations of Pacoima Dam after the Northridge earthquake. Comparison of the response due to nonuniform input with the response due to uniform input demonstrates the importance of accounting for spatial nonuniformity because of the significance that the pseudostatic component has for the response to nonuniform input. Copyright © 2006 John Wiley & Sons, Ltd.     

一种改进的近断层脉冲型地震动模拟方法

本文基于小波包技术的随机地震动模拟方法,提出一种改进的参数化随机近断层脉冲型地震动模拟方法。然后,通过识别和提取近断层脉冲型地震动数据库中脉冲型地震动的特征参数,建立了基于震源、传播路径和场地特征等参数的脉冲模型参数预测方程。最后,通过模拟实际记录和误差分析检验了改进的模拟方法的有效性。结果表明:应用改进的模拟方法得到的地震动时程无论在波形、频率特性还是峰值上均与实际记录具有较好的一致性。改进的模拟方法在保留地震动时频非平稳性的基础上,能够有效地提高近断层脉冲型地震动的模拟效果,并且能够很好地体现脉冲型地震动的主要特征。     

Velocity structure building and ground motion simulation of the 2014 Ludian Ms 6.5 Earthquake

This study constructs a 3 D velocity structure model of the Ludian region in the Yunnan province, southwestern China, and simulates ground motion propagation of the 2014 Ludian Ms 6.5 earthquake. It aims to construct the local velocity structure of the Ludian region in three dimensions and with high precision. The simulation, using the spectral element method, is validated by field data from the Ludian earthquake records. Thus, it demonstrates that the adopted key parameters, such as the seismic source mechanism, propagation medium and geographical features of the engineering site, are appropriated for the simulation. Meanwhile, the simulation generates the ground motion distribution of the study region with an earthquakeinduced landslide in Ludian earthquake.     

Velocity structure building and ground motion simulation of the 2014 Ludian <Emphasis Type="Italic">M</Emphasis>s 6.5 Earthquake

This study constructs a 3D velocity structure model of the Ludian region in the Yunnan province, southwestern China, and simulates ground motion propagation of the 2014 Ludian Ms 6.5 earthquake. It aims to construct the local velocity structure of the Ludian region in three dimensions and with high precision. The simulation, using the spectral element method, is validated by field data from the Ludian earthquake records. Thus, it demonstrates that the adopted key parameters, such as the seismic source mechanism, propagation medium and geographical features of the engineering site, are appropriated for the simulation. Meanwhile, the simulation generates the ground motion distribution of the study region with an earthquakeinduced landslide in Ludian earthquake.     

汶川地震强地面运动模拟

基于确定性震源模型的方法主要用于计算低频(11 Hz)地震动常用经验格林函数法或随机方法,对低频地震动模拟不够准确．本文在确定性震源模型方法基础上,尝试采用分解给定的震源模型的方法来模拟宽频带(0.1——10 Hz)强地面运动,即采用分级离散断层面和分解断层面破裂单元上升时间的方法,增加震源时间函数中的高频信号,从而避免了对地震记录丰富程度和准确性的依赖．文中模拟计算了汶川MS8.0地震在8个地震观测台的地震动,将模拟结果和观测记录进行了加速度时程曲线和傅里叶振幅谱的对比分析．对比结果显示,模拟估计的地震动峰值加速度和持续时间与观测记录的数据基本在plusmn;50%的精度范围内相同,傅里叶振幅谱显示模拟结果有得到10 Hz左右的高频成分. 四川盆地中的台站模拟结果高频衰减比观测记录要快,原因是模拟过程没有考虑场地效应.对强地震动模拟还是要综合考虑震源、传播路径和场地的影响．研究结果表明,此改变震源输入的确定性方法可应用于模拟近断层宽频强地面运动．      

Strong ground motion simulation during the November 1759 Earthquake along Serghaya Fault in the metropolitan of Damascus, Syria

The seismic hazard potential for metropolitan of Damascus, Syria is mainly controlled by earthquakes along Serghaya Fault which is a branch of Dead Sea Fault System. In this study, strong ground motion due to the November 1759 Earthquake along the fault of Serghaya was estimated with a numerical simulation technique. In the simulation, the Kostrov-like slip-velocity function was used as an input to the discrete wave number method to simulate the strong ground motions in a broadband frequency range. In order to model the incoherent rupture propagation which can excite large high-frequency waves, random numbers are added to arrival time of circular rupture front. MMI intensities calculated from the synthetic ground motions are compared with the observed values by Ambraseys and Barazangi (J Geophys Res 94:4007-4013, 1989). The calculated intensities are in good agreement with the observed ones at the most sites that validate appropriateness of the proposed source model. The PGA and PGV in the eastern region of Damascus city are higher than those in the western region due to the effects of local site amplification. The simulated high-frequency (1.0–6.0 Hz) ground motions for the sites in the Damascus city are higher than the design requirements defined by the Syrian building code. Furthermore, the simulated high-frequency ground motions for sites in the focal region are bigger than the design requirements in the case of the near-fault factors and are not considered. That demonstrates the appropriateness of considering the near-fault factors for a site near the focal region as introduced by the new building code.     

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.     

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.     

Canyon topography effects on ground motion at Feitsui damsite

The free-field accelerograms along Feitsui Canyon are analyzed and modeled by a numerical scheme to study the effect of canyon topography. Since six strong-motion accelerometers (SC1–SC6) were deployed along the Feitsui Canyon in 1991; there are 14 earthquakes (4.9≤M_L≤6.6) recorded by these stations until June 1996, but only five triggered all six stations. The maximum PGA value is 68.6 cm s⁻² recorded at station SC1. According to the present data, the effect of the dam on the ground motions at canyon stations can be negligible. The amplitude of ground motion on the slopes of the canyon is bigger than that at its trough. The integral equation method is applied to a two dimensional model of Feitsui Canyon to study the effects of the canyon topography. We choose the ground motion of SC3 or SC4 station at the trough of the canyon as the input motion for the model, which is then used to predict the ground motion at the other five stations. Apart from the earthquake close to the damsite, the simple model can reproduce the observed accelerations at all frequencies below 4 Hz. Overall, the numerical method can well predict the ground motion along the canyon, although the high-frequency simulation is underestimated.     

Response spectral relationships for rock sites derived from the component attenuation model

The seismological model was developed initially from the fundamental relationship between earthquake ground motion properties and the seismic moment generated at the source of the earthquake. Following two decades of continuous seismological research in the United States, seismological models which realistically account for both the source and path effects on the seismic shear waves have been developed and their accuracy rigorously verified (particularly in the long and medium period ranges). An important finding from the seismological research by Atkinson and Boore and their co‐investigators is the similarity of the average frequency characteristics of seismic waves generated at the source between the seemingly very different seismic environments of Eastern and Western North America (ENA and WNA, respectively). A generic definition of the average source properties of earthquakes has therefore been postulated, referred to herein as the generic source model. Further, the generic ‘hard rock’ crustal model which is characteristic of ENA and the generic ‘rock’ crustal model characteristic of WNA have been developed to combine with the generic source model, hence enabling simulations to be made of the important path‐related modifications to ground motions arising from different types of crustal rock materials. It has been found that the anelastic contribution to whole path attenuation is consistent between the ENA and WNA models, for earthquake ground motions (response spectral velocities and displacements) in the near and medium fields, indicating that differences in the ENA and WNA motions arise principally from the other forms of path‐related modifications, namely the mid‐crust amplification and the combined effect of the upper‐crust amplification and attenuation, both of which are significant only for the generic WNA ‘rock’ earthquake ground motions. This paper aims to demonstrate the effective utilization of the latest seismological model, comprising the generic source and crustal models, to develop a response spectral attenuation model for direct engineering applications. The developed attenuation model also comprises a source factor and several crustal (wave‐path modification) component factors, and thus has also been termed herein the component attenuation model (CAM). Generic attenuation relationships in CAM, which embrace both ENA and WNA conditions, have been developed using stochastic simulations. The crustal classification of a region outside North America can be based upon regional seismological and geological information. CAM is particularly useful for areas where local strong motion data are lacking for satisfactory empirical modelling. In the companion paper entitled ‘response spectrum modelling for rock sites in low and moderate seismicity regions combining velocity, displacement and acceleration predictions’, the CAM procedure has been incorporated into a response spectrum model which can be used to effectively define the seismic hazard of bedrock sites in low and moderate seismicity regions. This paper and the companion paper constitute the basis of a long‐term objective of the authors, to develop and effectively utilize the seismological model for engineering applications worldwide.