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.     

大型河谷场地地震动特性研究

采用有限差分方法,通过算例研究了大型河谷场地地震动特性分析中的人工边界的选取方法,对比分析了不同人工边界的选取对数值模拟结果的影响,确定了散射场地震反应分析输入边界的地震动输入方法,认为在进行有限差分动力计算时,模型两侧施加自由场边界的模拟效果要优于两侧施加粘性边界的模拟效果;同时,对FLAC计算软件进行二次开发,对2个地形差异较大的河谷场地,采用线性和非线性摩尔-库仑模型进行了地震反应对比分析,研究了河谷场地地震动幅值和频谱特性随地形变化的规律。模拟数据表明,河谷场地谷底处地震动基本无放大作用,地势凸起处放大作用则较为明显;当考虑土体非线性时,随着地震动强度的增加,放大作用逐渐减小;谱分析结果表明,地势凸起处受高频地震动的影响显著,而地势平坦的谷底则受低频地震动的影响显著。     

中国—欧洲共同体地震科学合作项目“北京西北延庆—怀来盆地地震学研究”进展

为探索强烈地震的前兆，在欧洲共同体的支持下，中国和英国的地震学家共同在北京西北的延庆—怀来盆地建立了一个数字地震台网．延庆—怀来数字地震台网由５个野外台和一个中心台构成。研究工作的重点是探索强烈地震前地震横波分裂现象中的参数变化．１９９５年４月３０日台网建成并投入运作．１９９５年７月下旬．怀来地区发生ＭＬ＝４．０地震及其余震，延庆—怀来数字地震台网很好地记录了这一地震序列。初步的分析表明．这一地震序列表明了怀来县城以南６—８ｋｍ处，怀来盆地南缘下部的断层活动，它与怀来地区地表有明显出露的主要断裂没有明显的直接联系。     

Estimation of Site Effects in Beijing City

— For the realistic modeling of the seismic ground motion in lateral heterogeneous anelastic media, the database of 3-D geophysical structures for Beijing City has been built up to model the seismic ground motion in the City, caused by the 1976 Tangshan and the 1998 Zhangbei earthquakes. The hybrid method, which combines the modal summation and the finite-difference algorithms, is used in the simulation. The modeling of the seismic ground motion, for both the Tangshan and the Zhangbei earthquakes, shows that the thick Quaternary sedimentary cover amplifies the peak values and increases the duration of the seismic ground motion in the northwestern part of the City. Therefore the thickness of the Quaternary sediments in Beijing City is the key factor controling the local ground effects. Four zones are defined on the base of the different thickness of the Quaternary sediments. The response spectra for each zone are computed, indicating that peak spectral values as high as 0.1 g are compatible with past seismicity and can be well exceeded if an event similar to the 1697 Sanhe-Pinggu occurs.     

2014年2月12日新疆于田地震强地面运动初步模拟及烈度预测

2014年2月12日新疆自治区于田县发生了M_S7.3级地震.在已知的三维介质模型、地形数据基础上,利用震源运动学初步反演模型（张勇等,2014）,作者采用曲线坐标网格三维曲线有限差分方法模拟了于田地震波场传播过程,并计算了模拟区域地震烈度分布.结果表明：地震最大烈度为7度,距离震中最近的于田县城烈度约为5度,断层西北侧地面震动略强于断层东南侧.另外模拟结果还显示地震动在山脊处具有较大的幅度.该地震本身主要能量释放区域在中地壳,主要滑动未及地表,因此对地表造成的破坏有限,这与目前尚无人员伤亡报告的情况相符合.     

Development of regional simulation of seismic ground‐motion and induced liquefaction enhanced by GPU computing

Much research has been conducted for physics‐based ground‐motion simulation to reproduce seismic response of soil and structures precisely and to mitigate damages caused by earthquakes. We aimed at enabling physics‐based ground‐motion simulations of complex three‐dimensional (3D) models with multiple materials, such as a digital twin (high‐fidelity 3D model of the physical world that is constructed in cyberspace). To perform one case of such simulation requires high computational cost and it is necessary to perform a number of simulations for the estimation of parameters or consideration of the uncertainty of underground soil structure data. To overcome this problem, we proposed a fast simulation method using graphics processing unit computing that enables a simulation with small computational resources. We developed a finite‐element‐based method for large‐scale 3D seismic response analysis with small programming effort and high maintainability by using OpenACC, a directive‐based parallel programming model. A lower precision variable format was introduced to achieve further speeding up of the simulation. For an example usage of the developed method, we applied the developed method to soil liquefaction analysis and conducted two sets of simulations that compared the effect of countermeasures against soil liquefaction: grid‐form ground improvement to strengthen the earthquake resistance of existing houses and replacement of liquefiable backfill soil of river wharves for seismic reinforcement of the wharf structure. The developed method accelerates the simulation and enables us to quantitatively estimate the effect of countermeasures using the high‐fidelity 3D soil‐structure models on a small cluster of computers.     

近断层速度脉冲地震动的三维有限差分模拟

根据台湾西部地质地貌特征和1999年集集M_W7.6地震的研究成果,建立三维速度结构模型和震源模型,并采用三维有限差分法对双冬断层可能产生的近断层脉冲型地震动进行数值模拟。结果表明,方向性效应引起的双向速度脉冲集中在垂直于断层滑动分量的方向上,而滑冲效应引起的单向速度脉冲则集中在平行于断层滑动分量的方向上。受方向性效应和上盘效应的共同调制,近断层脉冲型地震动反映出不对称带状分布的特征,速度脉冲主要分布在距离断层面约10 km的范围内。凹凸体的特性影响着地震动的时空分布,由地震波场显示南投和台中处于强地震动危险区。近场脉冲型地震动的研究对分析速度脉冲形成机理以及地震危险性有一定的参考意义。     

2D numerical simulations of earthquake ground motion: examples from the Marche Region,Italy

A detailed numerical simulation of the ground motion and a site response analysis for two towns in the Marche Region (Treia and Cagli) is carried out on the basis of structural models deduced from available geological and geophysical data. In both cases, the reference event is an M = 5.7 earthquake associated with a normal fault located beneath each town. The ground motion is computed using the 2D spectral element method (SPEM 2D). The method solves the propagation of the seismic field through complex geological structures and enables an estimate of the effects of deep crustal structure, superficial geology, and topography on ground motion. Numerical simulations of the seismic field are performed along 2D vertical planes containing the seismic source. Strong ground motion has not been yet recorded in the two towns; therefore, the numerical simulation of ground motion represents a way to overcome the lack of instrumental data. The simulations carried out for Treia show that ground motion is influenced by both source mechanism and effects due to propagation through the geological structure, while ground motion in Cagli features strong local effects, caused by the presence of alluvial deposits under a large area of the town.     

2014年8月3日云南鲁甸地震强地面运动初步模拟及烈度预测

2014年8月3日云南省昭通市鲁甸县发生了M_s6.5级地震.在已知的三维介质模型、地形数据基础上,利用震源运动学初步反演模型（张勇等,2014）,作者采用曲线坐标网格三维曲线有限差分方法模拟了鲁甸地震波场传播过程,并计算了模拟区域地震烈度分布.结果表明：地震最大烈度为Ⅶ度,破坏主要集中在鲁甸县以及巧家县、会泽县靠近鲁甸县的边界.另外,模拟结果还显示地震动在山峰、山脊处具有较大的幅度.计算表明断层东北侧的低速沉积盆地的波动放大效应加强了该地区的地震灾害.     

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.即唐山主震的发生使得其周边的应力场有一个瞬态的应力调整,唐山主震对后续余震的发生有促发作用.     

Broadband Ground Motion Reconstruction for the Kanto Basin during the 1923 Kanto Earthquake

Ground motions of the 1923 Kanto Earthquake inside the Kanto Basin are numerically simulated in a wide frequency range (0?C10?Hz) based on new knowledge of the earthquake??s source processes, the sedimentary structure of the basin, and techniques for generating broadband source models of great earthquakes. The Kanto Earthquake remains one of the most important exemplars for ground motion prediction in Japan due to its size, faulting geometry, and location beneath the densely populated Kanto sedimentary basin. We reconstruct a broadband source model of the 1923 Kanto Earthquake from inversion results by introducing small-scale heterogeneities. The corresponding ground motions are simulated using a hybrid technique comprising the following four calculations: (1) low-frequency ground motion of the engineering basement, modeled using a finite difference method; (2) high-frequency ground motion of the engineering basement, modeled using a stochastic Green??s function method; (3) total ground motion of the engineering basement (i.e. 1?+?2); and (4) ground motion at the surface in response to the total basement ground motion. We employ a recently developed three-dimensional (3D) velocity structure model of the Kanto Basin that incorporates prospecting data, microtremor observations and measurements derived from strong ground motion records. Our calculations reveal peak ground velocities (PGV) exceeding 50?cm/s in the area above the fault plane: to the south, where the fault plane is shallowest, PGV reaches 150?C200?cm/s at the engineering basement and 200?C250?cm/s at the surface. Intensity 7, the maximum value in the Japan Meteorological Agency??s intensity scale, is calculated to have occurred widely in Sagami Bay, which corresponds well with observed house-collapse rates due to the 1923 event. The modeling reveals a pronounced forward directivity effect for the area lying above the southern, shallow part of the fault plane. The high PGV and intensity seen above the southeastern corner of the fault plane and further east are largely due to this effect. Waveforms above the fault plane contain both short- and long-period components, but the short-period components are not observed further afield. Away from the fault, long-period waves (>2?s) dominate the ground motion, and in areas where the base of the third layer is relatively deep, the predominant period is >5?s. Levels of long-period ground motion in the southern part of the study area, around Sagami Bay and the southern parts of Boso Peninsula and Tokyo Bay, exceed that recorded at Tomakomai during the 2003 Tokachi-oki earthquake, when large oil storage tanks collapsed in response to sloshing generated by strong long-period motions.     

Efficient sampling and data reduction techniques for probabilistic seismic lifeline risk assessment

Probabilistic seismic risk assessment for spatially distributed lifelines is less straightforward than for individual structures. While procedures such as the ‘PEER framework’ have been developed for risk assessment of individual structures, these are not easily applicable to distributed lifeline systems, due to difficulties in describing ground‐motion intensity (e.g. spectral acceleration) over a region (in contrast to ground‐motion intensity at a single site, which is easily quantified using Probabilistic Seismic Hazard Analysis), and since the link between the ground‐motion intensities and lifeline performance is usually not available in closed form. As a result, Monte Carlo simulation (MCS) and its variants are well suited for characterizing ground motions and computing resulting losses to lifelines. This paper proposes a simulation‐based framework for developing a small but stochastically representative catalog of earthquake ground‐motion intensity maps that can be used for lifeline risk assessment. In this framework, Importance Sampling is used to preferentially sample ‘important’ ground‐motion intensity maps, and K‐Means Clustering is used to identify and combine redundant maps in order to obtain a small catalog. The effects of sampling and clustering are accounted for through a weighting on each remaining map, so that the resulting catalog is still a probabilistically correct representation. The feasibility of the proposed simulation framework is illustrated by using it to assess the seismic risk of a simplified model of the San Francisco Bay Area transportation network. A catalog of just 150 intensity maps is generated to represent hazard at 1038 sites from 10 regional fault segments causing earthquakes with magnitudes between five and eight. The risk estimates obtained using these maps are consistent with those obtained using conventional MCS utilizing many orders of magnitudes more ground‐motion intensity maps. Therefore, the proposed technique can be used to drastically reduce the computational expense of a simulation‐based risk assessment, without compromising the accuracy of the risk estimates. This will facilitate computationally intensive risk analysis of systems such as transportation networks. Finally, the study shows that the uncertainties in the ground‐motion intensities and the spatial correlations between ground‐motion intensities at various sites must be modeled in order to obtain unbiased estimates of lifeline risk. Copyright © 2010 John Wiley & Sons, Ltd.     

内马铁路一期工程局部不规则二维场地地震反应分析

基于内马铁路一期工程地质资料，利用ABAQUS有限元软件建立某局部不规则二维场地有限元计算模型，利用Python语言进行二次开发，编制二维等效线性化计算程序。开展50年超越概率63%、10%和2%水准下的土层地震反应分析计算，对不规则地形不同位置处的地震动峰值加速度和频谱特性进行对比分析。研究结果表明，不规则地形对地震动特性的影响显著。本文研究结果对内马铁路一期工程地震动参数确定具有指导意义，同时可为跨越不规则地形工程结构抗震设计提供参考。     

Site effects in 3D basins using 1D and 2D models: an evaluation of the differences based on simulations of the seismic response of Euroseistest

Site effects are one of the most predictable factors of destructive earthquake ground motion but results depend on the type of model chosen. We compare simulations of ground motion for a 3D model of the Mygdonian basin in northern Greece (Euroseistest) using different approximation for this basin. Site effects predicted using simple 1D models at many points inside the basin are compared to site effects predicted using four different 2D cross sections across the basin and with results for a full 3D simulation. Surface topography was neglected but anelastic attenuation was included in the simulations. We show that lateral heterogeneity may increase ground motion amplification by 100 %. Larger amplification is distributed in a wide frequency range, and amplification may occur at frequencies different from the expected resonant frequencies for the soil column. In contrast, on a different cross section, smaller conversion of incident energy into surface waves and larger dispersion leads to similar amplitudes of ground motion for 2D and 1D models. In general, results from 2D simulations are similar to those from a complete 3D model. 2D models may overestimate local surface wave amplitudes, especially when the boundaries of the basin are oblique to the selected cross section. However, the differences between 2D and 3D site effects are small, especially in regard of the difficulties and uncertainties associated to building a reliable 3D model for a large basin.     

地震动残差分析盆地附加放大效应:日本关东盆地为例

我国地震动预测及地震危险性分析通常仅考虑局部场地浅层岩土层对地震动的放大效应,不能考虑较大范围的地质条件影响,如沉积盆地厚沉积层对地震动的附加放大效应通常被忽略,造成盆地内地震动及地震危险性预测结果普遍被低估。本文以地震动观测记录数据充足的日本关东盆地为例,采用地震动残差分析方法评估盆地附加放大效应,分析覆盖层厚度、盆地内空间位置、震级、震源距对地震动放大效应的影响,建立关东盆地附加放大效应经验评估模型。分析表明:关东盆地附加放大效应与反应谱周期相关,整体上从短周期的1.0逐渐增大至周期为5s时的1.5,附加放大效应与覆盖层厚度相关性较小,主要受盆地空间位置和震源距的影响;盆地北部边缘及西北部地区附加放大效应更强烈,盆地南部附加放大效应较小,这可能与盆地边缘效应密切相关。本文建立的关东盆地附加放大效应经验模型略高于BSSA14和ASK14模型的放大效应预测。相关研究结果可用于我国地震动预测、下一代地震动区划图修订等。     

2017年九寨沟MS7.0地震场地特征分析

场地条件对地震动具有较大影响, 研究不同场地条件下的地震动特征对地震动的校正具有重要意义。 本文以九寨沟地震为例, 收集了66个台站的198条三分量强震观测记录和SRTM(Shuttle Radar Topography Mission)公里格网的DEM(Digital Elevation Model)数据, 从多个角度对场地特征进行分析。 首先根据坡度法使用DEM数据对九寨沟特征进行了场地分类, 然后讨论了不同场地类型下的加速度时程、 加速度反应谱、 地震动衰减等地震动特征。 研究表明, 缺少实测资料的情况下, 地形坡度可以作为Vs30(地表以下30 m范围的平均剪切波速度)的一种替代指标, 利用坡度法可以较快速地对场地进行分类; 利用强震观测记录能够从多角度对台站的场地特征进行分析, 不同场地类型对地震动影响不同, 其中, 土层对地震动具有明显的放大效应。 该结论可以为地震动结果的校正提供依据。     

9.21台湾集集地震中场地类别对地震动若干特性的影响

对9．21台湾集集地震主震的404组强震记录按照不同震中距和场地类别进行分组，然后对加速度水平分量和垂直分量的持时、峰值加速度、包络线、反应谱及其均方差作了，统计分析。结果表明，场地条件对不同分量加速度某些特征的影响较为显著，断层方向对地震动水平特性也有一定的影响。该结果对地震危险性分析以及结构抗震验算有一定参考价值。此外，对地震动三分量包络线的研究为三维地震动模拟提供了一定参考。     

Multi-Channel Satellite Image Analysis Using a Variational Approach

Currently, meteorological satellites provide multichannel image sequences including visible, temperature and water vapor channels. Based on a variational approach, we propose mathematical models to address some of the usual challenges in satellite image analysis such as: (i) the estimation and smoothing of the cloud structures by decoupling them into different layers depending on their altitudes, (ii) the estimation of the cloud structure motion by combining information from all the channels, and (iii) the 3D visualization of both the cloud structure and the estimated displacements. We include information of all the channels in a single variational motion estimation model. The associated Euler-Lagrange equations yield to a nonlinear system of partial differential equations that we solve numerically using finite-difference schemes. We illustrate the performance of the proposed models with numerical experiments on two multichannel satellite sequences of the North Atlantic, one of them from the Hurricane Vince. Based on a realistic synthetic ground truth motion, we show that our multichannel approach overcomes the single channel estimation for both the average Euclidean and angular errors.     

Parallel PSM/FDM Hybrid Simulation of Ground Motions from the 1999 Chi-Chi,Taiwan, Earthquake

-- A new technique for the parallel computing of 3-D seismic wave propagation simulation is developed by hybridizing the Fourier pseudospectral method (PSM) and the finite-difference method (FDM). This PSM/FDM hybrid offers a good speed-up rate using a large number of processors. To show the feasibility of the hybrid, a numerical 3-D simulation of strong ground motion was conducted for the 1999 Chi-Chi, Taiwan earthquake (M_w 7.6). Comparisons between the simulation results and observed waveforms from a dense strong ground motion network in Taiwan clearly demonstrate that the variation of the subsurface structure and the complex fault slip distribution greatly affect the damage during the Chi-Chi earthquake. The directivity effect of the fault rupture produced large S-wave pulses along the direction of the rupture propagation. Slips in the shallow part of the fault generate significant surface waves in Coastal Plain along the western coast. A large velocity gradient in the upper crust can propagate seismic waves to longer distances with minimum attenuation. The S waves and surface waves were finally amplified further by the site effect of low-velocity sediments in basins, and caused the significant disasters.