The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident...The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.展开更多
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.展开更多
Seismic wave reflection method is an advanced geophysical detection method in tunnel geological prediction.It is more sensitive and effective in detecting geological anomalies such as fault fracture zone and karst.In ...Seismic wave reflection method is an advanced geophysical detection method in tunnel geological prediction.It is more sensitive and effective in detecting geological anomalies such as fault fracture zone and karst.In order to verify the prediction efficacy and accuracy of the seismic wave reflection method with different instruments and equipment(tunnel geological prediction[TGP]/tunnel seismic prediction[TSP])and different vibration modes(hammering,explosives),a comparison test was carried out in Jinping Tunnel.The test results showed that the time-consumption of the hammering source was short,which can greatly reduce the impact on the construction site;different vibration sources methods of seismic wave reflection can predict the unfavorable geological sections accurately.展开更多
The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then sin...The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.展开更多
Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the...Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.展开更多
The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an anal...The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.展开更多
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...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.展开更多
This study quantifies seismic amplifications in near-shore arising from seaquakes. Within the Boundary Element Method, boundary elements are used to irradiate waves and force densities obtained for each element. Huyg...This study quantifies seismic amplifications in near-shore arising from seaquakes. Within the Boundary Element Method, boundary elements are used to irradiate waves and force densities obtained for each element. Huygens′ Principle is implemented since the diffracted waves are constructed at the boundary from which they are radiated, which is equivalent to Somigliana′s theorem. Application of boundary conditions leads to a system of integral equations of the Fredholm type of second kind and zero order. Several numerical configurations are analyzed： The first is used to verify the present formulation with ideal sea floor configurations to estimate seismic amplifications. With the formulation verified, simple slope configurations are studied to estimate spectra of seismic motions. It is found that P-waves can produce seismic amplifications from 1.2 to 3.9 times the amplitude of the incident wave. SV-waves can generate seismic amplifications up to 4.5 times the incident wave. Another relevant finding is that the highest amplifications are at the shore compared to the ones at the sea floor.展开更多
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.展开更多
Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using t...Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.展开更多
Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for mult...Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.展开更多
A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system o...A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.展开更多
基金National Natural Science Foundation of China under Grant No.51478247National Key Research and Development Program of China under Grant No.2016YFC1402800
文摘The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.
基金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.
文摘Seismic wave reflection method is an advanced geophysical detection method in tunnel geological prediction.It is more sensitive and effective in detecting geological anomalies such as fault fracture zone and karst.In order to verify the prediction efficacy and accuracy of the seismic wave reflection method with different instruments and equipment(tunnel geological prediction[TGP]/tunnel seismic prediction[TSP])and different vibration modes(hammering,explosives),a comparison test was carried out in Jinping Tunnel.The test results showed that the time-consumption of the hammering source was short,which can greatly reduce the impact on the construction site;different vibration sources methods of seismic wave reflection can predict the unfavorable geological sections accurately.
基金supported by the National Natural Science Foundation of China under grant Nos.40474012,40874020 and 40821062
文摘The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.
基金supported by National Natural Science Foundation of China (No. 41174085)Chinese Academy of Sciences (KZZD-EW-TZ-19)China Geological Survey (12120113101400)
文摘Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.
文摘The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.
基金supported by the National Natural Science Foundation of China(Grants 41374046 and41174034)
文摘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.
文摘This study quantifies seismic amplifications in near-shore arising from seaquakes. Within the Boundary Element Method, boundary elements are used to irradiate waves and force densities obtained for each element. Huygens′ Principle is implemented since the diffracted waves are constructed at the boundary from which they are radiated, which is equivalent to Somigliana′s theorem. Application of boundary conditions leads to a system of integral equations of the Fredholm type of second kind and zero order. Several numerical configurations are analyzed： The first is used to verify the present formulation with ideal sea floor configurations to estimate seismic amplifications. With the formulation verified, simple slope configurations are studied to estimate spectra of seismic motions. It is found that P-waves can produce seismic amplifications from 1.2 to 3.9 times the amplitude of the incident wave. SV-waves can generate seismic amplifications up to 4.5 times the incident wave. Another relevant finding is that the highest amplifications are at the shore compared to the ones at the sea floor.
基金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.
基金funded by the Earthquake Science Fund of Xinjiang (201404)Seismic Tracing Oriented Task of China Earthquake Administration(2015010122)
文摘Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.
文摘Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.
文摘A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.
