3-D physical model in strong ground motion numerical simulation： A case study of Kunming basin

Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and projectconstruction,while numerical simulation of strong ground motion is an important method for seismic hazard pre-diction and assessment.A 3-D physical model in conformity with real strata configuration of(mainly)the Quater-nary is a prerequisite to ensure the reliability of the simulation results.In this paper,we give a detailed account ofthe technical scheme and process for creating a 3-D physical model in Kunming basin.The data used are synthe-sized from seismogeological data,borehole data,topographic data,digital elevation mode(DEM)data,seismicexploration results and wave velocity measurements.Stratigraphic division is based mainly on shear wave velocity,with strata sequence taken into consideration.The model construction is finally accomplished with ArcGIS andmany relevant programming techniques via layer-by-layer stacking(in depth direction)of the adjacent mediuminterfaces(meshes).Meanwhile,a database of 3-D physical models is set up,which provides model data and pa-rameters for strong ground motion simulation.Some processing methods and significant issues are also addressedin the paper in accordance with different types of exploration and experimental data.     

长周期地震动三维有限元数值模拟方法

采用数值方法对长周期地震动进行模拟已成为当前强地面运动预测的研究方向之一。本文提出了一种有限元并行计算的高效地震动数值模拟方法,可以同时考虑震源、传播途径和场地条件对地震动的影响,并根据昆明地区未来可能发生的设定地震对昆明市的影响地震动进行了预测。这一方法可以应用于国内已探测的发震活断层进行地震危害性评价工作。     

三维地震动场数值模拟并行计算系统

并行计算系统是进行三维复杂场地地震动场数值模拟的基础。本文作者利用10台微机构建了一套并行计算集群(常称为Beowulf系统),节点微机上运行Windows2000操作系统,采用MPI标准的MPICH1．2．4版本作为并行计算的支撑环境,并以100M bps高速交换式以太网作为互连网络。利用编制的并行计算程序对该并行计算系统性能进行了测试。结果表明：该系统能够达到非常高的并行加速比和并行效率。该系统将用于三维地震动场数值模拟研究。     

Basin and crustal velocity structure models for the simulation of strong ground motions in the Kinki area,Japan

We constructed a prototype of the basin and crustal structure model for the Kinki area, southwest of Japan, for the simulation of strong ground motions of hypothetical crustal and subduction earthquakes. We collected results of the deep seismic velocity profiles obtained by the reflection experiments and seismic imaging results, which were conducted in the Kinki area. The obtained profiles give underground velocity structures of the crust, from the surface to the subducting slab. We also gather the basin velocity structure information of the Osaka, Kyoto, Nara, and Ohmi basins. To examine the applicability of the constructed velocity structure model to the ground motion simulation, we simulated waveforms of an intermediate size event occurred near the source area of the hypothetical subduction earthquakes. Simulated ground motions using the basin and crustal velocity structure model are fairly well reproducing the observations at most of stations, and the constructed basin and crustal velocity structure model is applicable for the long-period ground motion simulations.     

Prediction of long-period ground motions from huge subduction earthquakes in Osaka,Japan

There is a high possibility of reoccurrence of the Tonankai and Nankai earthquakes along the Nankai Trough in Japan. It is very important to predict the long-period ground motions from the next Tonankai and Nankai earthquakes with moment magnitudes of 8.1 and 8.4, respectively, to mitigate their disastrous effects. In this study, long-period (>2.5 s) ground motions were predicted using an earthquake scenario proposed by the Headquarters for Earthquake Research Promotion in Japan. The calculations were performed using a fourth-order finite difference method with a variable spacing staggered-grid in the frequency range 0.05–0.4 Hz. The attenuation characteristics (Q) in the finite difference simulations were assumed to be proportional to frequency (f) and S-wave velocity (V _s) represented by Q = f · V _s / 2. Such optimum attenuation characteristic for the sedimentary layers in the Osaka basin was obtained empirically by comparing the observed motions during the actual M5.5 event with the modeling results. We used the velocity structure model of the Osaka basin consisting of three sedimentary layers on bedrock. The characteristics of the predicted long-period ground motions from the next Tonankai and Nankai earthquakes depend significantly on the complex thickness distribution of the sediments inside the basin. The duration of the predicted long-period ground motions in the city of Osaka is more than 4 min, and the largest peak ground velocities (PGVs) exceed 80 cm/s. The predominant period is 5 to 6 s. These preliminary results indicate the possibility of earthquake damage because of future subduction earthquakes in large-scale constructions such as tall buildings, long-span bridges, and oil storage tanks in the Osaka area.     

北京盆地结构对长周期地震动反应谱的影响

建立了北京盆地含起伏地形的地下介质三维速度结构模型, 以有限差分方法对盆地内的3次地震动进行了数值模拟． 计算了地表加速度反应谱, 论证了反应谱的有效周期范围, 从而获得了盆地结构对长周期段反应谱的平均放大系数, 并对不同周期时反应谱的平均放大系数进行了分区, 为今后城市的合理规划和长周期建筑物的建设提供了抗震设防参考． 最后对反应谱平均放大系数与等效沉积物厚度之间的关系进行了回归, 得出了反应谱放大系数随沉积物厚度增大而增强, 以及随周期增大而减弱等结论．      

3D crustal structure and long-period ground motions from a M9.0 megathrust earthquake in the Pacific Northwest region

We have developed a community velocity model for the Pacific Northwest region from northern California to southern Canada and carried out the first 3D simulation of a M_w 9.0 megathrust earthquake rupturing along the Cascadia subduction zone using a parallel supercomputer. A long-period (<0.5 Hz) source model was designed by mapping the inversion results for the December 26, 2004 Sumatra–Andaman earthquake (Han et al., Science 313(5787):658–662, 2006) onto the Cascadia subduction zone. Representative peak ground velocities for the metropolitan centers of the region include 42 cm/s in the Seattle area and 8–20 cm/s in the Tacoma, Olympia, Vancouver, and Portland areas. Combined with an extended duration of the shaking up to 5 min, these long-period ground motions may inflict significant damage on the built environment, in particular on the highrises in downtown Seattle.     

Prediction of ground motion in the Osaka sedimentary basin associated with the hypothetical Nankai earthquake

We studied the long-period ground motions in the Osaka sedimentary basin, Japan, which contains a 1- to 3-km thickness of sediments and is the site of many buildings or construction structures with long-natural period. We simulated the broadband ground motions likely to be produced by the hypothetical Nankai earthquake: the earthquake expected to give rise to the most severe long-period ground motion within the basin. For the simulation, we constructed multiscale heterogeneous source models based on the Central Disaster Management Council of Japan (CDMC) source model and adopted a hybrid computation method in which long-period motion and short-period motion are computed using a 3-D finite difference method and the stochastic Green’s function method, respectively. In computing long-period motions, we used a 3-D structure model of the crust and the Osaka sedimentary basin. The ground motions are estimated to have peak velocities of 50–90 cm/s, prolonged durations exceeding 300 s, and long predominant periods of 5–10 s in the area with great thickness of sediments. The predominant periods are in agreement with an approximate evaluation by 4 H/V _s where H and V _s are the thickness of the sediment and the average S wave velocity, respectively.     

长周期地震动衰减关系研究的迫切性

通过对长周期震害及长周期地震动研究的历史进行回顾指出,目前,大多数反应谱的周期还不够长,反应谱所反映的谱形与距离、震级和场地条件的关系还很有限,长周期地震动的研究仍需要不断完善,未来的长周期地震动的研究应注重软土场地的大震远场效应。     

An evaluation of 3-D velocity models of the Kanto basin for long-period ground motion simulations

We performed three-dimensional (3-D) finite difference simulations of long-period ground motions (2–10 s) in the Kanto basin using the Japan Seismic Hazard Information Station (J-SHIS 2009), Yamada and Yamanaka (Exploration Geophysics 65(3):139–150, 2012) (YY), and Head Quarter for Earthquake Research Promotion (HERP 2012) velocity models for two intermediate depth (68–80 km) moderate earthquakes (Mw 5.8–5.9), which occurred beneath the Kanto basin. The models primarily differ in the basic data set used in the construction of the velocity models. The J-SHIS and HERP models are the results of integration of mainly geological, geophysical, and earthquake data. On the other hand, the YY model is oriented towards the microtremor-array-observation data. We obtained a goodness of fit between the observed and synthetic data based on three parameters, peak ground velocities (PGVs), smoothed Fourier spectra (FFT), and cross-correlations, using an algorithm proposed by Olsen and Mayhew (Seism Res Lett 81:715–723, 2010). We found that the three models reproduced the PGVs and FFT satisfactorily at most sites. However, the models performed poorly in terms of cross-correlations especially at the basin edges. We found that the synthetics using the YY model overestimate the observed waveforms at several sites located in the areas having V _s 0.3 km/s in the top layer; on the other hand, the J-SHIS and HERP models explain the waveforms better at the sites and perform similarly at most sites. We also found that the J-SHIS and HERP models consist of thick sediments beneath some sites, where the YY model is preferable. Thus, we have concluded that the models require revisions for the reliable prediction of long-period ground motions from future large earthquakes.     

Lessons from the 2003 Tokachi-oki, Japan, earthquake for prediction of long-period strong ground motions and sloshing damage to oil storage tanks

The 2003 Tokachi-oki earthquake (M _w 8.0) in northern Japan generated large-amplitude long-period (4–8 s) ground motions in the Yufutsu sedimentary basin, causing severe damage to seven large oil storage tanks with floating roof structures because of severe sloshing of oil. The 30,000–40,000-m³ tanks having suffered the severe damage such as fires and sinking of floating roofs experienced the sloshing with large amplitudes exceeding 3 m in which the fundamental mode was predominant. The second mode of sloshing was also excited in the 110,000-m³ tanks in which their floating roofs sank into oil, indicating that the higher modes of sloshing as well as the fundamental mode should be considered in damage prediction. The strong ground motion recordings demonstrated the earthquake dependency of predominant periods and the substantial spatial variation of the long-period shaking observed within the Yufutsu basin, meaning the necessity of source- and site-specific prediction of long-period strong ground motions. The two-dimensional numerical modeling suggested the importance of detailed structures of soft near-surface sediments as well as deep basin structure for accurate prediction of long-period strong ground motions in deep sedimentary basins.     

特殊长周期地震动的参数特征研究

近断层脉冲型地震动和远场软土层场地类谐和地震动是两类特殊的长周期地震动,当前的规范均很少对这两类地震作用进行具体的规定。研究了近断层脉冲型和远场类谐和地震动的幅值、幅值比（V／A,D／V）、傅里叶幅值谱和反应谱的差别,分析了相位角和作用循环周期数对简单脉冲的影响,并用于解释两类特殊地震动的工程特征。以集集地震动为数据基础,分析了两类长周期地震动的傅里叶谱和反应谱特征;将平均加速度和位移规准反应谱分别与规范设计谱进行了比较。建议设计谱在长周期段考虑近断层作用和软土场地面波效应的影响。     

1999年台湾集集大地震强地面运动的模拟

利用有限断层模型对1999年台湾集集大地震所产生的近场强地面运动进行了数值模拟,尝试了将具有明确物理意义的震源动力模型的研究结果应用于描述震源过程。计算结果显示,基岩台或接近工程基岩台的计算模拟记录同观测记录符合得很好。进一步研究表明,当综合考虑场地效应对地震动影响时,模拟计算结果得到了明显的改善,同实际观测记录更为相符。本文的研究结果说明,为了科学合理地预测大地震所产生的强地面运动,应综合考虑影响地震所产生的地面运动的三要素,即:震源过程、传播路径及场地效应。     

Stochastic models for simulation of strong ground motion in Iceland

Two types of modelling approaches for simulating ground motion in Iceland are studied and compared. The first type of models, named discrete‐time series models (ARMA), are based solely on measured acceleration in earthquakes occurring in Iceland. The second type of models are based on a theoretical seismic source model called the extended Brune model. Based on measured acceleration in Iceland during the period 1986–1996, the parameters for the extended Brune models have been estimated. The seismic source models are presented here as ARMA models, which simplifies the simulation process. A single‐layer soil amplification model is used in conjunction with the extended Brune model to estimate local site amplification. Emphasis is put on the ground motion models representing the variability in the measured earthquakes, with respect to energy, duration and frequency content. Demonstration is made using these models for constructing linear and non‐linear probabilistic response spectra using a discretised version of the Bouc–Wen model for the hysteresis of the second‐order system. Copyright © 2001 John Wiley & Sons, Ltd.     

强地面运动数值模拟中三维物理模型的建立方法——以昆明盆地为例

基于活动断层的地震危害性综合评价可为城市规划和工程建设提供科学依据, 强地面运动数值模拟则是进行地震危害性预测和评价的重要方法, 而建立以第四系为主体的符合真实地层结构的三维物理模型是保证数值模拟结果可靠性的必要条件之一. 本文以昆明盆地为例, 综合利用地震地质、钻孔、地形地貌、DEM、地震勘探、波速测试等资料,以剪切波速为主要分层指标,参考地层层序, 结合ArcGIS等多种软件和相关的编程技术,采用将相邻介质分界面（层网）在深度方向上逐层叠加的方法建立三维物理模型;详细阐述了建立三维物理模型的主要技术思路和实现过程, 同时建立三维物理模型的数据库,为强地面运动数值模拟提供模型数据和参数. 文中还针对不同类型的探测、实验数据,总结了相应的技术处理方法和需要注意的问题.      

地震动加速度过程的小波模拟

本文应用基数B－样条小波模拟地震动加速度过程，将小波分析引入地震动加速度过程的研究，实现了地震动加速度的模拟，并可给出其解析表达式，便于设计地震动的调整及地震反应的求解。     

考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟

2008年5月12日四川汶川地区发生M_W7.9地震,震中位置103.4°E,31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡,且因为高山等地形复杂区域抢险救灾的艰巨性,为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响,建立了包含地形起伏影响及去除地形影响的两类模型.同时,依据震源破裂过程运动学反演结果,建立了包含障碍体破裂过程的震源滑动模型,实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法,通过对地震波传播过程的数值计算和后处理分析,模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用,距离断层破裂带越近,地形效应越明显,反之,距离越远,则地形效应越微弱;（2）因为地形高差与障碍体的影响,地震造成的峰值可能出现在震中区域之外;（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近;而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段,处清平、安县附近;对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义.

     

基于GP14.3运动学混合震源模型和SPECFEM 3D谱元法的宽频地震动模拟

基于确定性物理模型的全过程地震动模拟是现代地震工程的重要发展方向.然而受限于合理震源模型和计算资源需求,目前模拟的有效频率还多处于低频范围,难以满足工程结构敏感频带（5~10 Hz或更高）需求.本文即借助运动学混合震源模型能激发宽频地震波和谱元法空间高精度及计算收敛快的优势,首先将确定性的凹凸体震源模型与GP14.3随机震源模型结合得到有限断层运动学混合震源模型,进而将上述混合震源模型开发到SPECFEM 3D谱元法开源代码中,实现了基于谱元法和运动学混合震源模型的全过程宽频带地震动模拟.将方法首先应用于一维波速结构模型0~10 Hz地震动模拟,通过与频率波数域（FK）方法结果进行比较,验证了方法的精度;进而应用于2021年5月21日云南漾濞6.4级地震0.1~5 Hz地震动模拟,通过与4个台站的时程记录和相应反应谱的比较,以及与NGA-West2地震动衰减方程在频率0.1~5 Hz的反应谱的比较,检验了方法的适用性;最后给出了漾濞地区的地震动峰值加速度（PGA）和峰值速度（PGV）云图,分析了漾濞地震下近场强地面运动的空间分布特征.结果显示,震中PGA接近400 cm·s^-2,PGV达到45 cm·s^-1,烈度达到Ⅸ度,且受局部地形起伏影响,大理以及洱海西侧位置出现高烈度异常区.

     

基于希尔伯特变换的非平稳地震动模拟方法的验证

基于希尔伯特变换的非平稳地震动模拟方法能够生成一系列地震动样本,具有与给定天然地震动相似的非平稳特征. 该方法在模拟天然地震动时域波形的形态和能量的时频分布方面,效果非常显著且已得到验证. 本文以理想单自由度体系弹塑性地震反应为基础,深入验证该方法在模拟天然地震动对结构影响效应方面的效果. 研究结果表明,该方法能够较好地模拟天然地震动的工程特性,其生成的地震动样本的峰值特性及其引起的结构弹塑性地震响应具有与天然地震动相似的特征.