We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units...We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional 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.展开更多
The analysis of seismic wave propagation and amplification in complex geological structures requires efficient numerical methods.In this article,following up on recent studies devoted to the formulation,implementation...The analysis of seismic wave propagation and amplification in complex geological structures requires efficient numerical methods.In this article,following up on recent studies devoted to the formulation,implementation and evaluation of 3-D single-and multi-region elastodynamic fast multipole boundary element methods(FM-BEMs),a simple preconditioning strategy is proposed.Its efficiency is demonstrated on both the single-andmulti-region versions using benchmark examples(scattering of plane waves by canyons and basins).Finally,the preconditioned FM-BEM is applied to the scattering of plane seismic waves in an actual configuration(alpine basin of Grenoble,France),for which the high velocity contrast is seen to significantly affect the overall efficiency of the multi-region FM-BEM.展开更多
We apply the spectral-element method(SEM),a high-order finite-element method(FEM) to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometri...We apply the spectral-element method(SEM),a high-order finite-element method(FEM) to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities,and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots,synthetic seismograms,and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems.展开更多
The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogen...The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous media, anisotropic and porous media, surface wave and seismic wave inversion, and seismic wave study in prospecting and logging problems. Important projects in the current studies on seismic wave is suggested as the development of high efficient numerical methods, and applying them to the studies of excitation and propagation of seismic waves in complex media and strong ground motion, which will form a foundation for refined earthquake hazard analysis and prediction.展开更多
Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomo...Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomography,seismic-wave scattering and its applications to studying Earth’s heterogeneity,oneway wave propagation and one-return wave modeling,beamlet and dreamlet applications,strong non-linear full-waveform inversion,and direct envelop inversion.展开更多
A stable and explicit second order accurate finite difference method for the elastic wave equation in curvilinear coordinates is presented.The discretization of the spatial operators in the method is shown to be self-...A stable and explicit second order accurate finite difference method for the elastic wave equation in curvilinear coordinates is presented.The discretization of the spatial operators in the method is shown to be self-adjoint for free-surface,Dirichlet and periodic boundary conditions.The fully discrete version of the method conserves a discrete energy to machine precision.展开更多
基金supported by the School of Energy Resources at the University of WyomingThe GPU hardware used in this study was purchased using the NSF Grant EAR-0930040
文摘We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional 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.
文摘The analysis of seismic wave propagation and amplification in complex geological structures requires efficient numerical methods.In this article,following up on recent studies devoted to the formulation,implementation and evaluation of 3-D single-and multi-region elastodynamic fast multipole boundary element methods(FM-BEMs),a simple preconditioning strategy is proposed.Its efficiency is demonstrated on both the single-andmulti-region versions using benchmark examples(scattering of plane waves by canyons and basins).Finally,the preconditioned FM-BEM is applied to the scattering of plane seismic waves in an actual configuration(alpine basin of Grenoble,France),for which the high velocity contrast is seen to significantly affect the overall efficiency of the multi-region FM-BEM.
基金supported by the Natural Sciences and Engineering Research Council of Canada (NSERC)Center for Excellence in Mining Innovations (CEMI,through SUMIT project)+2 种基金Computations for this study were performed on hardwares purchased through the combined funding of Canada Foundation for Innovation (CFI)Ontario Research Fund (ORF)University of Toronto Startup Fund
文摘We apply the spectral-element method(SEM),a high-order finite-element method(FEM) to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities,and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots,synthetic seismograms,and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems.
基金State National Science Foundation of China (grant No. 40134010).
文摘The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous media, anisotropic and porous media, surface wave and seismic wave inversion, and seismic wave study in prospecting and logging problems. Important projects in the current studies on seismic wave is suggested as the development of high efficient numerical methods, and applying them to the studies of excitation and propagation of seismic waves in complex media and strong ground motion, which will form a foundation for refined earthquake hazard analysis and prediction.
文摘Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomography,seismic-wave scattering and its applications to studying Earth’s heterogeneity,oneway wave propagation and one-return wave modeling,beamlet and dreamlet applications,strong non-linear full-waveform inversion,and direct envelop inversion.
基金This work performed under the auspices of the U.S.Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
文摘A stable and explicit second order accurate finite difference method for the elastic wave equation in curvilinear coordinates is presented.The discretization of the spatial operators in the method is shown to be self-adjoint for free-surface,Dirichlet and periodic boundary conditions.The fully discrete version of the method conserves a discrete energy to machine precision.