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.

Seismic wave propagation in viscoelastic and saturated rock

onsidering engineering practice, the viscoelastic twophase model is adopted, seismic wave propagation in saturated rock is studied. Not only the effect of the viscosity of rock skeleton but also the effect of ground water on the propagation of the seismic wave can be considered by this model, the propagation characteristics of seismic wave in saturated rock can be understood comprehensively and the model is more reasonable than other model by which seismic wave propagation is studied. The effect of frequency, water content and viscosity constant on the wave velocity and attenuation are studied by numerical examples and some valuable conclusions are drawn.

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

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.

Forward prediction for tunnel geology and classification of surrounding rock based on seismic wave velocity layered tomography

Excavation under complex geological conditions requires effective and accurate geological forward-prospecting to detect the unfavorable geological structure and estimate the classification of surround-ing rock in front of the tunnel face.In this work,a forward-prediction method for tunnel geology and classification of surrounding rock is developed based on seismic wave velocity layered tomography.In particular,for the problem of strong multi-solution of wave velocity inversion caused by few ray paths in the narrow space of the tunnel,a layered inversion based on regularization is proposed.By reducing the inversion area of each iteration step and applying straight-line interface assumption,the convergence and accuracy of wave velocity inversion are effectively improved.Furthermore,a surrounding rock classification network based on autoencoder is constructed.The mapping relationship between wave velocity and classification of surrounding rock is established with density,Poisson’s ratio and elastic modulus as links.Two numerical examples with geological conditions similar to that in the field tunnel and a field case study in an urban subway tunnel verify the potential of the proposed method for practical application.

Comparative Analysis of Application of Seismic Wave Reflection Method in Advanced Geological Prediction

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.

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.

Characterization Method of Explosion Seismic Wave

A new characterization method of explosion seismic wave is suggested on the basis of the analysis of experimental measured results. The seismic wave function is resolved into amplitude modulation part and random one. For the latter, the fractal dimension and the relevant characterization parameters are yielded by using the Weirstrass Mandelbrot (W M) fractal function. In contrast with conventional statistical parameters, the new set of parameters is independent of the chosen time length scales and the measuring instruments. A modeling example is presented which shows that the theoretical results are in agreement with the experimental results.

Numerical simulation of seismic wave propagation in complex media by convolu-tional differentiator

We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.

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.

Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries

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.

Study on Dynamic Response of the “Dualistic” Structure Rock Slope with Seismic Wave Theory

Currently, scant attention has been paid to the theoretical analysis on dynamic response mechanism of the ＂Dualistic＂ structure roek slope. The analysis presented here provides insight into the dynamic response of the ＂Dualistie＂ structure rock slope. By investigating the principle of energy distribution, it is shown that the effect of a joint plays a significant role in slope stability analysis. A dynamic reflection and transmission model （RTM） for the ＂Dualistic＂ structure rock slope and explicit dynamic equations are established to analyze the dynamic response of a slope, based on the theory of elastic mechanics and the principle of seismic wave propagation. The theoretical simulation solutions show that the dynamic response of the ＂Dualistic＂ structure rock slope （soft-hard） model is greater than that of the ＂Dualistic＂ strueture rock slope （hard-soft） model, especially in the slope crest. The magnifying effect of rigid foundation on the dynamic response is more obvious than that of soft foundation. With the amplitude increasing, the cracks could be found in the right slope （soft-hard） crest. The crest failure is firstly observed in the right slope （soft-hard） during the experimental process. The reliability of theoretical model is also investigated by experiment analysis. The conclusions derived in this paper could also be used in future evaluations of Multi-layer rock slopes.

Physics informed machine learning: Seismic wave equation

Similar to many fields of sciences,recent deep learning advances have been applied extensively in geosciences for both small-and large-scale problems.However,the necessity of using large training data and the’black box’nature of learning have limited them in practice and difficult to interpret.Furthermore,including the governing equations and physical facts in such methods is also another challenge,which entails either ignoring the physics or simplifying them using unrealistic data.To address such issues,physics informed machine learning methods have been developed which can integrate the governing physics law into the learning process.In this work,a 1-dimensional(1 D)time-dependent seismic wave equation is considered and solved using two methods,namely Gaussian process(GP)and physics informed neural networks.We show that these meshless methods are trained by smaller amount of data and can predict the solution of the equation with even high accuracy.They are also capable of inverting any parameter involved in the governing equation such as wave velocity in our case.Results show that the GP can predict the solution of the seismic wave equation with a lower level of error,while our developed neural network is more accurate for velocity(P-and S-wave)and density inversion.

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.

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

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.

Stochastic seismic wave interaction with a rock joint having Coulomb-slip behavior

Seismic wave interaction with a slippery rock joint with an arbitrary impinging angle is analytically studied based on the conservation of momentum on the wave fronts. Based on the displacement discontinuity method, the wave propagation equations are derived for incident P- and S-waves. By comparison, the calculated transmission and reflection coefficients for normal incident waves are the same as the existing results, which proves the wave propagation equation obtained in the paper is correct. The wave propagation derived in the context can be applied to incident waves with different waveforms. Stochastic seismic waves are then used to analyze the seismic wave interaction with the slippery rock joint, where the stochastic seismic waves are generated from frequency spectra. The parametric studies are carried out to investigate the effect of type, intensity and impinging angle of the incident seismic waves on the wave propagation across the slippery rock joint.

Influence of seismic wave type and incident direction on the dynamic response of tall concrete-faced rockfill dams

Owing to the stochastic behavior of earthquakes and complex crustal structure,wave type and incident direction are uncertain when seismic waves arrive at a structure.In addition,because of the different types of the structures and terrains,the traveling wave effects have different influences on the dynamic response of the structures.For the tall concrete-faced rockfill dam(CFRD),it is not only built in the complex terrain such as river valley,but also its height has reached 300 m level,which puts forward higher requirements for the seismic safety of the anti-seepage system mainly comprising concrete face slabs,especially the accurate location of the weak area in seism.Considering the limitations of the traditional uniform vibration analysis method,we implemented an efficient dynamic interaction analysis between a tall CFRD and its foundation using a non-uniform wave input method with a viscous-spring artificial boundary and equivalent nodal loads.This method was then applied to investigate the dynamic stress distribution on the concrete face slabs for different seismic wave types and incident directions.The results indicate that dam-foundation interactions behave differently at different wave incident angles,and that the traveling wave effect becomes more evident in valley topography.Seismic wave type and incident direction dramatically influenced stress in the face slab,and the extreme stress values and distribution law will vary under oblique wave incidence.The influence of the incident direction on slab stress was particularly apparent when SH-waves arrived from the left bank.Specifically,the extreme stress values in the face slab increased with an increasing incident angle.Interestingly,the locations of the extreme stress values changed mainly along the axis of the dam,and did not exhibit large changes in height.The seismic safety of CFRDs is therefore lower at higher incident angles from an anti-seepage perspective.Therefore,it is necessary to consider both the seismic wave type and incident direction during seismic capacity evaluations of tall CFRDs.