Scattering of seismic waves by three-dimensional large-scale hill topography simulated by a fast parallel IBEM

To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual(GMRES) algorithm and constructing an Open MP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.

Seismic Waves in Layered Media

This paper is organized as follows. After a discussion of the differential equations for wave propagation in the horizontally stratified medium and of the initial and boundary conditions, the displacements are derived on the free surface of the layered medium for plane waves when a point source is located on the s-th imaginary boundary at the depth -s （physical parameters of the layers s and （s ＋ 1） are put to be identical）. Then, the source will be represented as a single force of arbitrary orientation and a general moment tensor point source. Further, ＂a primary field＂ for a point source will be introduced. Matrix method for the solution of the direct seismic problem is considered based on the matrix method of Thomson-Haskell and its modifications.

Response of Surface Soft Soil and Fine-Grained Sand Layers to Seismic Waves

Shanghai is located in eastern China and is built on overburden soil layers. It can be seen from the Mexico M S=8.1 earthquake on September 19, 1985 and the Hanshin M S=7.4 earthquake on January 17,1995 that heavy casualties and property losses have a direct relationship with overburden soil layers. Ground motions caused by earthquakes are significantly amplified when passing through the soil layers. Under the influence of these amplified motions, building structures, whose nature frequency is within the frequency band of soil amplification response, will experience more severe damage than those built on bedrock. Therefore, engineering seismologists have paid considerable attention the amplification responses in the Shanghai overburden soil layers. The amplification responses of soil and sand layers in this paper are given by the M L=4.1 earthquake in Nantong, Jiangsu Province on December 25, 2001 at 31.8° N, 120.9° E. It can be seen that the responses of soil and sand layers are very different. That is important.

A study of damping factors in perfectly matched layers for the numerical simulation of seismic waves

When simulating seismic wave propagation in free space, it is essential to introduce absorbing boundary conditions to eliminate reflections from artificially trtmcated boundaries. In this paper, a damping factor referred to as the Gaussian dmping factor is proposed. The Gaussian damping factor is based on the idea of perfectly matched layers （PMLs）. This work presents a detailed analysis of the theoretical foundations and advantages of the Gaussian damping factor. Additionally, numerical experiments for the simulation of seismic waves are presented based on two numerical models： a homogeneous model and a multi-layer model. The results show that the proposed factor works better. The Gaussian damping factor achieves a higher Signal-to-Noise Ratio （SNR） than previously used factors when using same number of PMLs, and requires less PMLs than other methods to achieve an identical SNR.

Numerical simulation on radiation and energy of blast-induced seismic waves in deep rock masses

With regard to blasting in deep rock masses,it is commonly thought that an increase in the in-situ stress will change the blast-induced rock crack propagation and ultimately affect rock fragmentation.However,little attention has been given to the change in seismic wave radiation when the fractured zone changes with the in-situ stress.In this study,the influences of in-situ stress on blast-induced rock fracture and seismic wave radiation are numerically investigated by a coupled SPH-FEM simulation method.The results show that the change in blast-induced rock fracture with in-situ stress has a considerable effect on the seismic wave energy and composition.As the in-situ stress level increases,the size of the fractured zone is significantly reduced,and more explosion energy is transformed into seismic energy.A reduction in the size of the fractured zone(seismic wave source zone)results in a higher frequency content of the seismic waves.In a nonhydrostatic in-situ stress field,blast-induced cracks are most suppressed in the direction of the minimum in-situ stress,and thus the seismic waves generated in this direction have the highest energy density.In addition to P-waves,Swaves are also generated when a circular explosive is detonated in a nonhydrostatic in-situ stress field.The S-waves result from the asymmetrical release of rock strain energy due to the anisotropic blast-induced fracture pattern.

Nonlinear finite element analysis of effect of seismic waves on dynamic response of Shiziping dam

Many high earth-rockfill dams are constructed in the west of China. The seismic intensity at the dam site is usually very high, thus it is of great importance to ensure the safety of the dam in meizoseismal area. A 3D FEM model is established to analyze the seismic responses of Shiziping earth-rockfill dam. The nonlinear elastic Duncan-Chang constitutive model and the equivalent viscoelastic constitutive model are used to simulate the static and dynamic stress strain relationships of the dam materials, respectively. Four groups of seismic waves are inputted from the top of the bedrock to analyze the dynamic responses of the dam. The numerical results show that the calculated dynamic magnification factors display a good consistency with the specification values. The site spectrum results in larger acceleration response than the specification spectrum. The analysis of relative dynamic displacement indicates that the displacement at the downstream side of the dam is larger than that at the upstream side. The displacement response reduces from the center of river valley to two banks. The displacement responses corresponding to the specification spectrum are a little smaller than those corresponding to the site spectrum. The analysis of shear stress indicates that a large shear stress area appears in the upstream overburden layer, where the shear stress caused by site waves is larger than that caused by specification waves. The analysis of dynamic principal stress indicates that the minimum dynamic stresses in corridor caused by specification and site waves have little difference. The maximum and minimum dynamic stresses are relatively large at two sides. The largest tensile stress occurs at two sides of the floor of grouting corridor, which may result in the crack near the corridor side. The numerical results present good consistency with the observation data of the grouting corridor in Wenchuan earthquake.

A laboratory method to simulate seismic waves induced by underground explosions

The seismic waves induced by underground explosions generate geological hazards affecting deep buried tunnels such as rockbursts and engineering-induced earthquakes. This issue is difficult to study through full-scale testing due to the expense and unpredictable danger. To solve this problem, the authors developed experimental apparatus and presented a laboratory method to simulate seismic waves induced by underground explosions. In this apparatus, a combined structure of a diffusive-shaped water capsule and a special-shaped oil capsule was designed. This structure can provide an applied confining stress and freely transmit the stress wave generated by external impact. Therefore, the coupled loading of in situ stress and seismic waves induced by underground explosions in the deep rock mass was simulated. The positive pressure time and peak value of the stress wave could be adjusted by changing the pulse-shaper and the initial impact energy. The obtained stress waves in the experiments correspond to that generated by 0.15-120 kt of TNT equivalent explosion at a scaled distance of 89.9-207.44 m/kt.

Monitoring and Prediction of the Vibration Intensity of Seismic Waves Induced in Underwater Rock by Underwater Drilling and Blasting

All underwater drilling and blasting operations generate seismic waves.However,due to a lack of suitable vibration sensing instruments,most studies on the propagation of seismic waves have been limited to shorelines near construction areas or wharfs,whereas comparatively few studies have beerconducted on the larger seafloor itself.To address this gap,a seafloor vibration sensor system was developed and applied in this study that consists of an autonomous acquisition storage terminal,soft-ware platform,and hole-plugging device that was designed to record the blasting vibration intensities received through submarine rocks at a given measurement point.Additionally,dimensional analyses were used to derive a predictive equation for the strength of blast vibrations that considered the in fluence of the water depth.By combining reliable vibration data obtained using the sensor system in submarine rock and the developed predictive equation,it was determined that the water depth was ar important factor influencing the measured vibration strength.The results using the newly derivedequation were compared to those determined using the Sadowski equation,which is commonly used on land,and it was found that predictions using the derived equation were closer to the experimental values with an average error of less than 10%,representing a significant improvement.Based on these results the developed sensor system and preliminary theoretical basis was deemed suitable for studying the propagation behavior of submarine seismic waves generated by underwater drilling and blasting operations.

Modeling of low-frequency seismic waves in a shallow sea using the staggered grid difference method

Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships. To obtain the propagation characteristics of ship seismic waves, an algorithm for calculating Seismic waves at the seafloor is presented based on the staggered-grid finite difference method. The accuracy of the algorithm was tested by comparison with analytical solutions. Numerical simulation of seismic waves generated by a low-frequency point sotmd source in a typical shallow sea environment was carried out. Using various source frequencies and locations in the numerical simulation, we show that the seismic waves in the near field are composed mostly of transmitted S-waves and interface waves while transmitted P-waves are weak near the seafloor. However, in the far field, the wave components of the seismic wave are mainly normal modes and interface waves, with the latter being relatively strong in the waveforms, As the source frequency decreases, the normal modes become smaller and the interface waves dominate the time series of the seismic waves.

Studies on seismic waves

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.

Seismic Waves in a Layered Half-Space from an Arbitrary Buried Source

In this paper, a transfer matrix and a three-dimensional dynamic response of a layered half-space to an arbitrary buried source are derived with the aid of a technique which combines the Laplace and two-dimensional Fourier transforms in a rectangular coordinate system. This method is clear in concept, and the corresponding formulas given in the paper are simple and convenient for marine seismic prospecting and other fields' applications. An example is presented and the calculated results are in good agreement with those of the finite element method (FEM).

Energy flux characteristics of seismic waves at the interface between soil layers with different saturations

Based on the multiphase poroelasticity theory describing the propagation of waves in the unsaturated fluid-saturated porous medium,the reflection and transmission coefficients of the seismic waves at the interface between soil layers with different saturations are obtained.Our unsaturated model consists of a deformable skeleton in which two compressible and viscous fluids(i.e.,water and gas)flow in the interstices.Three compressional waves(i.e.,P1,P2,and P3 waves)and one shear(i.e.,S wave)wave exist in the unsaturated soils.The expressions for the energy ratios of the various reflected and transmitted waves at the interface during the transmission and reflection processes are presented in explicit forms accordingly.At last,numerical computations are performed and the results obtained are respectively depicted graphically.The variation of the energy ratios with the incident angle,wave frequency and saturation degrees of the upper and lower soil layers is illustrated in detail.The calculation results show that the allocation of incident seismic waves at the interface is influenced not only by the angle and frequency of the incident seismic waves,but also by the saturations of the upper and lower soil layers.It is also verified that,at the interface,the sum of energy ratios of the reflected and transmitted waves is approximately equal to unity as was expected.This study is of importance to several fields such as geotechnical engineering,seismology,and geophysics.

Reflection and transmission of seismic waves at an interface betwen two saturated soils

Based on the modified Biot model for asturated soils, taking the compressibilities of the grains and the pore fluid as well as the viscous coupling into account, the reflection and transmission of seismic aves at an interface between two saturated soils are studied in this paper. A formula is derived for calculation of the amplitude reflection and transmission coefficients of various waves. A aumerical investigation of the dependence of the coefficients on the angle of incidence and the frequency is performed. This study is of a value for seismological studies and geophysical exploration.

Influence of long period seismic waves induced by the Taiwan large earthquake on Shanghai area

In this paper a group of long period seismic waves in Shanghai area induced by Taiwan large earthquake is presented by using the method of semi empirical Green function, the period is up to 20 s. Such seismic waves can be used as a reference curve to test the strength of long period structures and their aseismic design. The long period part of seismic influence curve presented in 'Architecture Aseismic Design Code'(GBJ11 89) is less than 3 s, and uncertainties exist in the effects of earthquake safety evaluation. This research will be able to eliminate these shortages.

A real time processing system of seismic waves using personal computers-Function and characteristics

Seismological Bureau of Sichuan Province, Chengdu 610041, China2) Center for Analysis and Prediction, State Seismological Bureau, Beijing 100036, China3) Observation Center for Prediction of Earthquakes and Volcanic Eruptions, Faculty of Sciences, Tohoku University, Sendai 98077, Japan

A comparative analysis of seismic response of shallow buried underground structure under incident P,SV and Rayleigh waves

In this study,A time-domain seismic response analysis method and a calculation model of the underground structure that can realize the input of seismic P,SV and Rayleigh waves are established,based on the viscoelastic artificial boundary elements and the boundary substructure method for seismic wave input.After verifying the calculation accuracy,a comparative study on seismic response of a shallow-buried,double-deck,double-span subway station structure under incident P,SV and Rayleigh waves is conducted.The research results show that there are certain differences in the cross-sectional internal force distribution characteristics of underground structures under different types of seismic waves.The research results show that there are certain differences in the internal force distribution characteristics of underground structures under different types of seismic waves.At the bottom of the side wall,the top and bottom of the center pillar of the underground structure,the section bending moments of the underground structure under the incidences of SV wave and Rayleigh wave are relatively close,and are significantly larger than the calculation result under the incidence of P wave.At the center of the side wall and the top floor of the structure,the peak value of the cross-sectional internal force under the incident Rayleigh wave is larger than the calculation result under SV wave.In addition,the floor of the underground structure under Rayleigh waves vibrates in both the horizontal and vertical directions,and the magnification effect in the vertical direction is more significant.Considering that the current seismic research of underground structures mainly considers the effect of body waves such as the shear waves,sufficient attention should be paid to the incidence of Rayleigh waves in the future seismic design of shallow underground structures.

Rayleigh waves from correlation of seismic noise in Great Island of Tierra del Fuego,Argentina:Constraints on upper crustal structure

In this study, the ambient seismic noise cross-correlation technique is applied to estimate the upper structure of the crust beneath Great Island of Tierra del Fuego（TdF）, Argentina, by the analysis of shortperiod Rayleigh wave group velocities. The island, situated in the southernmost South America, is a key area of investigation among the interaction between the South American and Scotia plates and is considered as a very seismically active one. Through cross-correlating the vertical components of ambient seismic noise registered at four broadband stations in TdF, we were able to extract Rayleigh waves which were used to estimate group velocities in the period band of 2.5-16 s using a timefrequency analysis. Although ambient noise sources are distributed in homogeneously, robust empirical Green＇s functions could be recovered from the cross-correlation of 12 months of ambient noise, The observed group velocities were inverted considering a non-linear iterative damped least-squares inversion procedure and several 1-D shear wave velocity models of the upper crust were obtained.According to the inversion results, the S-wave velocity ranges between 1.75 and 3,7 km/s in the first10 km of crust, depending on the pair of stations considered. These results are in agreement to the major known surface and sub-surface geological and tectonic features known in the area. This study represents the first ambient seismic noise analysis in TdF in order to constraint the upper crust beneath this region.It can also be considered as a successful feasibility study for future analyses with a denser station deployment for a more detailed imaging of structure.

Nonlinear Time History Analysis for the Different Column Orientations under Seismic Wave Synthetic Approach

The significant impact of earthquakes on human lives and the built environment underscores the extensive human and economic losses caused by structural collapses. Over the years, researchers have focused on improving seismic design to mitigate earthquake-induced damages and enhance structural performance. In this study, a specific reinforced concrete (RC) frame structure at Kyungpook National University, designed for educational purposes, is analyzed as a representative case. Utilizing SAP 2000, the research conducts a nonlinear time history analysis to assess the structural performance under seismic conditions. The primary objective is to evaluate the influence of different column section designs, while maintaining identical column section areas, on structural behavior. The study employs two distinct seismic waves from Abeno (ABN) and Takatori (TKT) for the analysis, comparing the structural performance under varying seismic conditions. Key aspects examined include displacement, base shear force, base moment, joint radians, and layer displacement angle. This research is anticipated to serve as a valuable reference for seismic restraint reinforcement work on RC buildings, enriching the methods used for evaluating structures through nonlinear time history analysis based on the synthetic seismic wave approach.

Experimental Research on Amplitude Change of Blasting Seismic Wave with Topography

The propagation characteristics of the amplitude of the blasting seismic wave under the conditions of various topographies are approached by means of experiments. Some factors affecting the effects of quake insulation groove, such as the size, the depth and the position of the quake insulation groove, are studied. The amplitudes of the blasting seismic waves under the conditions of the different sizes of the quake insulation groove are measured. According to the experiments, the effects of the quake insulation groove are related to the position, the distance, the energy of the explosion source and the size of the quake insulation groove itself. The farther it is from the explosion source, the smaller the energy is. The lower its position is and the larger its size is, the more remarkable the effects of the quake insulation groove are.

Research on the Seismic Wave Field of Karst Cavern Reservoirs near Deep Carbonate Weathered Crusts

Fracture and cavern hydrocarbon reservoirs in carbonates are an important pool type worldwide. The karst cavern reservoirs are easiest to identify on seismic reflection data. The prediction, exploration, and development of this type of reservoir require theoretical research on seismic wave fields reflected from complex inhomogeneous media. We compute synthetic seismic sections for fluidfilled cavern reservoirs of various heights and widths using random media models and inhomogeneous media elastic wave equations. Results indicate that even caverns significantly smaller than 1/ 4 wavelength are detectible on conventional band-width seismic sections as diffractions migrated into bead-type events. Diffraction amplitude is a function of cavern height and width. We introduce a width-amplitude factor which can be used to calculate the diffraction amplitude of a cavern with a limited width from the diffraction amplitude computed for an infinitely wide cavern.