Global SH-wave propagation in a 2D whole Moon model using the parallel hybrid PSM/FDM method

We present numerical modeling of SH-wave propagation for the recently proposed whole Moon model and try to improve our understanding of lunar seismic wave propagation. We use a hybrid PSM/FDM method on staggered grids to solve the wave equations and implement the calculation on a parallel PC cluster to improve the computing efficiency. Features of global SH-wave propagation are firstly discussed for a 100-km shallow and900-km deep moonquakes, respectively. Effects of frequency range and lateral variation of crust thickness are then investigated with various models. Our synthetic waveforms are finally compared with observed Apollo data to show the features of wave propagation that were produced by our model and those not reproduced by our models. Our numerical modeling show that the low-velocity upper crust plays significant role in the development of reverberating wave trains. Increasing frequency enhances the strength and duration of the reverberations.Surface multiples dominate wavefields for shallow event.Core–mantle reflections can be clearly identified for deep event at low frequency. The layered whole Moon model and the low-velocity upper crust produce the reverberating wave trains following each phases consistent with observation. However, more realistic Moon model should be considered in order to explain the strong and slow decay scattering between various phases shown on observation data.     

Global SH-wavefield calculation for a two-dimensional whole- Earth model with the parallel hybrid PSM/FDM algorithm

We present a parallel hybrid algorithm based on pseudospectral method （PSM） and finite difference method （FDM） for two-dimensional （2-D） global SH- wavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corre- sponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.     

Global SH-wavefield calculation for a two-dimensional whole-Earth model with the parallel hybrid PSM/FDM algorithm

We present a parallel hybrid algorithm based on pseudospectral method (PSM) and finite difference method (FDM) for two-dimensional (2-D) global SH-wavefield simulation. The whole-Earth model is taken as a cross section of spherical Earth, and corresponding wave equations are defined in 2-D cylindrical coordinates. Spatial derivatives in the wave equations are approximated with efficient and high accuracy PSM in the lateral and high-order FDM in the radial direction on staggered grids. This algorithm allows us to divide the whole-Earth into sub-domains in radial direction and implement efficient parallel computing on PC cluster, while retains high accuracy and efficiency of PSM in lateral direction. A transformation of moment tensor between 3-D spherical Earth and our 2-D model was proposed to give corresponding moment tensor components used in 2-D modeling. Comparison of modeling results with those obtained by direct solution method shows very good accuracy of our algorithm. We also demonstrate its feasibility with a lateral heterogeneous whole-Earth model with localized velocity perturbation.     

地震波传播的三维伪谱和高阶有限差分混合方法并行模拟

基于交错网格伪谱法和高阶精度有限差分方法,发展了模拟非均匀介质地震波传播的三维伪谱和有限差分混合算法.该方法在两个水平方向利用交错网格伪谱算子计算空间微分,保留了该方法高效、高精度的优势,在垂直方向采用交错网格高阶精度有限差分算子实现空间微分计算．利用有限差分方法的局部性特征,将三维计算区域在垂直方向上划分为一系列子区域,并分配给不同的处理器,实现了在并行计算机集群上的三维并行计算．通过模拟算例,与离散波数法比较,检验了该算法的精度．为了检验该方法的实用性,在64个处理器上,对三维沉积盆地模型进行了67108864个网格点的并行计算,模拟的波场主频率为1.25Hz,讨论了沉积盆地深度对三维沉积盆地地面运动的影响．      

基于图形处理器的伪谱和高阶有限 差分混合方法地震波数值模拟

伪谱和高阶有限差分混合方法, 在垂直方向采用交错网格有限差分算子, 利用其并行程度高的特点, 在水平方向采用伪谱算子, 保留其高精度的优势, 是计算地震波场的有效方法. 图形处理器(graphic processing unit, 简写为GPU) 由于其高度并行性, 在计算此类问题中有显著的优势. 由英伟达(NVIDIA)公司推出的统一计算设备架构(compute unified device architecture, 简写为CUDA)平台极大地简化了GPU编程的难度. 为提高计算效率, 本文实现了基于CUDA 平台的混合方法二维地震波场模拟. 然后基于二维均匀介质模型将CPU与GPU版本的运行时间进行对比. 实际测试结果表明, 基于CUDA 的并行模拟方法在保证计算精度的同时显著地提高了计算速度, 为开展大规模非均匀地球介质地震波传播数值模拟提供了一种可选的方法.      

基于伪谱和有限差分混合方法的兰州盆地强地面运动二维数值模拟

利用伪谱和有限差分混合方法对兰州盆地进行二维强地面运动模拟.假设3个不同的震源深度5、10、20 km的同一地震下,以一个5层的二维剖面作为模型,以此来研究地震波传播过程并分析复杂的二维非均质路径效应,以及局部地下速度构造对强地面运动振幅加乘的影响,从而了解路径效应对强地面运动最大峰值位移的影响程度.研究结果显示:伪谱和有限差分混合方法模拟地震波场结合了有限差分法和伪谱法的优点,弥补了二者的不足,能较好地处理介质不连续面的计算,同时保证了和伪谱法相当的计算精度.通过地震波场模拟可知沉积盆地的强地面运动较基岩相比具有放大效应,震源深度的不同对兰州盆地垂直分量地面最大峰值位移影响较大,在震源深度为5 km时最小,为0.06 cm;震源深度10 km时在水平剖面30～55 km范围内峰值位移最大,达到0.14 cm;但对水平分量的峰值位移影响较小.     

Accelerating the discontinuous Galerkin method for seismic wave propagation simulations using multiple GPUs with CUDA and MPI

We have successfully ported an arbitrary high-order discontinuous Galerkin method for solving the three-dimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) of NVIDIA and Message Passing Interface (MPI) and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix–matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.     

非均匀介质全局广义R/T递推传播矩阵法二维SH地震波模拟研究

地球深部圈层及沉积盆地是一种分区非均匀介质系统,其中不规则地层边界（含起伏地表）对地震波的主要特征有显著影响,而地层的随机非均匀性则主要影响地震波的散射和衰减特征.为了精确刻画不规则地层边界对地震波的反射、透射效应以及非均质体散射引起的地震波衰减效应,全局广义R/T递推传播矩阵法（GGRTM）被提出并逐步发展成为继有限元和有限差分方法之后的另一种复杂介质高精度地震波传播半解析求解方法.在已有的此类方法中,不规则边界均匀地层GGRTM法的优势在于对不规则地层边界的反射和透射效应的准确模拟,而非均质地层薄板化GGRTM法则能准确描述非均质体散射对地震波衰减的影响.本文吸收这两种已有方法的优势,提出了一种考虑非均匀介质、不规则边界的全局广义R/T递推传播矩阵混合方法,并将其用于对边界不规则、层内非均质的复杂模型的二维SH波场模拟.随后在本文方法与边界元法对比研究的基础上讨论了方法的模拟精度.研究结果表明本文提出的混合法是一种解决复杂模型高精度地震模拟的有效方法.

     

Pseudo‐spectral method using rotated staggered grid for elastic wave propagation in 3D arbitrary anisotropic media

Staggering grid is a very effective way to reduce the Nyquist errors and to suppress the non‐causal ringing artefacts in the pseudo‐spectral solution of first‐order elastic wave equations. However, the straightforward use of a staggered‐grid pseudo‐spectral method is problematic for simulating wave propagation when the anisotropy level is greater than orthorhombic or when the anisotropic symmetries are not aligned with the computational grids. Inspired by the idea of rotated staggered‐grid finite‐difference method, we propose a modified pseudo‐spectral method for wave propagation in arbitrary anisotropic media. Compared with an existing remedy of staggered‐grid pseudo‐spectral method based on stiffness matrix decomposition and a possible alternative using the Lebedev grids, the rotated staggered‐grid‐based pseudo‐spectral method possesses the best balance between the mitigation of artefacts and efficiency. A 2D example on a transversely isotropic model with tilted symmetry axis verifies its effectiveness to suppress the ringing artefacts. Two 3D examples of increasing anisotropy levels demonstrate that the rotated staggered‐grid‐based pseudo‐spectral method can successfully simulate complex wavefields in such anisotropic formations.     

流体饱和半空间二维地形三分量弹性波散射间接边界元模拟

无限长局部地形地震波斜入射响应问题称为二维三分量问题,在计算量远小于三维的情况下,一定程度上反映了近地表场地的三维动力响应特征.基于天然土体的成层性及固液两相耦合特性,以层状多孔介质内部移动线荷载(孔隙水压)动力格林函数作为基本解,开展流体饱和半空间二维地形三分量弹性波散射的2.5维间接边界元模拟研究.总场响应由自由波...     

Estimation of uncertainty propagation in flood inundation mapping using a 1‐D hydraulic model

A need for more accurate flood inundation maps has recently arisen because of the increasing frequency and extremity of flood events. The accuracy of flood inundation maps is determined by the uncertainty propagated from all of the variables involved in the overall process of flood inundation modelling. Despite our advanced understanding of flood progression, it is impossible to eliminate the uncertainty because of the constraints involving cost, time, knowledge, and technology. Nevertheless, uncertainty analysis in flood inundation mapping can provide useful information for flood risk management. The twin objectives of this study were firstly to estimate the propagated uncertainty rates of key variables in flood inundation mapping by using the first‐order approximation method and secondly to evaluate the relative sensitivities of the model variables by using the Hornberger–Spear–Young (HSY) method. Monte Carlo simulations using the Hydrologic Engineering Center's River Analysis System and triangle‐based interpolation were performed to investigate the uncertainty arising from discharge, topography, and Manning's n in the East Fork of the White River near Seymour, Indiana, and in Strouds Creek in Orange County, North Carolina. We found that the uncertainty of a single variable is propagated differently to the flood inundation area depending on the effects of other variables in the overall process. The uncertainty was linearly/nonlinearly propagated corresponding to valley shapes of the reaches. In addition, the HSY sensitivity analysis revealed the topography of Seymour reach and the discharge of Strouds Creek to be major contributors to the change of flood inundation area. Copyright © 2014 John Wiley & Sons, Ltd.     

Seismic site effects by an optimized 2D BE/FE method II. Quantification of site effects in two-dimensional sedimentary valleys

This paper deals with the evaluation of seismic site effects due to the local topographical and geotechnical characteristics. The amplification of surface motions is calculated by a numerical method combining finite elements in the near field and boundary elements in the far field (FEM/BEM). The numerical technique is improved by time truncation. In the first part of this article, the accuracy and the relevance of this optimized method are presented. Moreover, parametric studies are done on slopes, ridges and canyons to characterize topographical site effects. The second part deals with sedimentary valleys. The complexity of the combination of geometrical and sedimentary effects is underlined. Extensive parametrical studies are done to discriminate the topographical and geotechnical effects on seismic ground movement amplifications in two-dimensional irregular configurations. Characteristic coefficients are defined to predict the amplifications of horizontal displacements. The accuracy of this quantitative evaluation technique is tested and discussed.     

Parameter identification for on‐line model updating in hybrid simulations using a gradient‐based method

To improve the efficiency of model fitting, parameter identification techniques have been actively investigated. Recently, the applications of parameter identification migrated from off‐line model fitting to on‐line model updating. The objective of this study is to develop a gradient‐based method for model updating to advance hybrid simulation also called hybrid test. A novel modification of the proposed method, which can reduce the number of design variables to improve the identification efficiency, is illustrated in detail. To investigate the model updating, simulated hybrid tests were conducted with a 5‐story steel frame equipped with buckling‐restrained braces (BRBs) utilized in the shaking table tests conducted in E‐Defense in Japan in 2009. The calibrated analytical model that was verified with the test results can serve as the reference model. In the simulated hybrid tests, the physical BRB substructure is numerically simulated by utilizing a truss element with the 2‐surface model identical to the part of the reference model. Such numerical verification allows simulation of measurement errors for investigation on the performance of the proposed method. Moreover, the feasibility of sharing the identified parameter values, which were obtained from the physical substructure responses, with the relevant numerical models is also verified with the artificial component responses derived from the physical experiments.     

Numerical modelling of wave propagation in anisotropic soil using a displacement unit-impulse-response-based formulation of the scaled boundary finite element method

An efficient method for modelling the propagation of elastic waves in unbounded domains is developed. It is applicable to soil–structure interaction problems involving scalar and vector waves, unbounded domains of arbitrary geometry and anisotropic soil. The scaled boundary finite element method is employed to derive a novel equation for the displacement unit-impulse response matrix on the soil–structure interface. The proposed method is based on a piecewise linear approximation of the first derivative of the displacement unit-impulse response matrix and on the introduction of an extrapolation parameter in order to improve the numerical stability. In combination, these two ideas allow for the choice of significantly larger time steps compared to conventional methods, and thus lead to increased efficiency. As the displacement unit-impulse response approaches zero, the convolution integral representing the force–displacement relationship can be truncated. After the truncation the computational effort only increases linearly with time. Thus, a considerable reduction of computational effort is achieved in a time domain analysis. Numerical examples demonstrate the accuracy and high efficiency of the new method for two-dimensional soil–structure interaction problems.     

3D finite-difference simulations of strong ground motions in the Yanhuai area, Beijing, China during the 1720 Shacheng earthquake (Ms 7.0) using a stochastic finite-fault model

Three-dimensional finite-difference (3D-FD) simulations of earthquake wave propagations in the Yanhuai area were performed for the 1720 Shacheng earthquake (M_s 7.0) using a stochastic finite-fault model, running on a parallel supercomputer Hitachi-SR8000 (http://www.lrz-muenchen.de). A stochastic finite-fault model was implemented into the 3D-FD program. The basic idea of the stochastic finite-fault model is that the fault plane can be subdivided into several subfaults (or elements, sources). Radiation from a large earthquake is the sum of contributions from all subfaults with proper time delays, each of which acts as a small independent double-couple point source. Heterogeneity of the fault rupture process was modeled by randomizing the location of initial rupture (hypocenter), slip vectors (slip, rake), and rise-times of subfaults in this study. A 3D velocity model of the Yanhuai area was constructed based on studies that analyzed available geological and geophysical information. A grid increment of 75 m in three directions was used in the 3D-FD simulation, which made it possible to capture the short period information with a resolution as low as 0.5 s in sediment regions. The uncertainties of simulated results caused by the stochastic finite-fault model were studied with a homogeneous 3D model. We found that the effects of the randomness of source on simulated ground motions are only limited in near-fault-region including the surface exposure of the fault and its vicinities, which occupies about 5% of the whole study area. This article presents an integrated approach for simulating the strong ground motions for engineering purpose using the 3D-FD method. Such simulations would be useful for hazard mapping, land using planning, insurance rate assessment, particularly in planning, preparedness, and coordinating emergency response, which must be based on realistic situations induced by concrete (historic or scenario) earthquakes.     

迭代法在三维网络模型最短路径射线追踪中的应用

目前有许多基于二维网格模型的射线追踪算法,但应用与三维网格模型的射线追踪算法较少,主要是因为三维网格射线追踪算法在算法实现、编程以及图形显示方面都有一定难度.本文提出了一种适用于三维网格射线追踪的迭代算法,并应用于几种理论模型,同时实现了三维绘图显示,通过研究证明迭代算法能在一定程度上提高三维网格射线追踪的精度和速度.     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明：新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.

     

Modeling the potential distribution of shallow-seated landslides using the weights of evidence method and a logistic regression model： a case study of the Sabae Area, Japan

A number of statistical methods are typically used to effectively predict potential landslide distributions. In this study two multivariate statistical analysis methods were used （weights of evidence and logistic regression） to predict the potential distribution of shallow-seated landslides in the Kamikawachi area of Sabae City, Fukui Prefecture, Japan. First, the dependent variable （shallow-seated landslides） was divided into presence and absence, and the independent variables （environmental factors such as slope and altitude） were categorized according to their characteristics. Then, using the weights of evidence （WE） method, the weights of pairs comprising presence （w^＋（i）） or absence （w^-（i））, and the contrast values for each category of independent variable （evidence）, were calculated, Using the method that integrated the weights of evidence method and a logistic regression model, score values were calculated for each category of independent variable. Based on these contrast values, three models were selected to sum the score values of every gird in the study area. According to a receiver operating characteristic curve analysis （ROC）, model 2 yielded the best fit for predicting the potential distribution of shallow-seated landslide hazards, with 89% correctness and a 54.5% hit ratio when the occurrence probability （OP） of landslides was 70%. The model was tested using data from an area close to the study region, and showed 94% correctness and a hit ratio of 45.7% when the OP of landslides was 70%. Finally, the potential distribution of shallow-seated landslides, based on the OP, was mapped using a geographical information system.     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.