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.

Determination of Direction of Arrival of Seismic Wave by a Single Tri-axial Fiber Optic Geophone

A fiber Bragg grating (FBG) geophone and a surface seismic wave-based algorithm for detecting the direction of arrival (DOA) are described. The operational principle of FBG geophone is introduced and illustrated with systematic experimental data, demonstrating an improved FBG geophone with many advantages over the conventional geophones. An innovative, robust, and simple algorithm is developed for obtaining the bearing information on the seismic events, such as people walking, or vehicles moving. Such DOA estimate is based on the interactions and projections of surface-propagating seismic waves generated by the moving personnel or vehicles with a single tri-axial seismic sensor based on FBGs. Of particular interest is the case when the distance between the source of the seismic wave and the detector is less than or comparable to one wavelength (less than 100 m), corresponding to near-field detection, where an effective method of DOA finding lacks.

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 of the Response Features of Apparent Seismic Wave Velocity Ratio in a Horizontal Layered Medium

This paper deals with the response features of AR(apparent ratio of seismic wave velocities to the changes of TR(true ratio of wave velocities)in the horizontal layered model by mathematical modeling.The results show that:(1)the response features of AR are associated with the parameters of the structure and its dynamic changes,and the relative position between the hypocenters and the monitoring networks,showing complicated patterns strongly related to the concrete paths of propagation of seismic waves from the source to the receiver in the observatories of the network;(2)the depth of the seismic source would have important influence on the response features of AR,especially the capacity to carry the anomalous information in the condition of the earth media,being in the anomalous state would be greater for those earthquakes which occur inside the anomalous layers than those underneath the anomalous layers;(3)the response features of AR are clearly related to the changes of TR(true ratio of wave velocities)instead of changes of wave velocities themselves,i.e.the response could be small as the changes in TR is small even in the case of large changes in the wave velocities.It is suggested that more attention must be paid to all these features in combination with detailed investigation of the velocity structure of the earth media in the study region and best fitting of precise hypocenter locations when one wants to obtain the reliable precursors from the changes in AR.

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.

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.

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.

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.

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.

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.

Propagation simulation and attenuation law analysis of seismic wave in elastoplastic tunnel rock medium

The energy caused by the dynamic impact in mining engineering forth release and spread by the way of seismic waves, monitoring is an effective way for forecasting mine dynamical disasters, such as rockburst and coal and gas outburst. Three-dimensional dynamic model was built to simulate the propagating progress of seismic waves in the elastoplastic tunnel rock and analyzed the propagating law of perturbation acceleration around tunnel, based on the finite element dynamic analysis software ANSYS/L S-DYNA. The simulation results indicate that： （1） The propagation attenuation of seismic wave is a negative index relationship; （2） The acceleration amplitude of seismic wave decays rapidly in near-field and decays slowly in far-field; （3） When the perturbation is generated in the dead ahead of tunnel, the acceleration of seismic wave become smaller and smaller away from the roadway-rib;（4） The elastic and plastic stress state of tunnel rock is also an important factor for propagation process of wave, the energy of seismic wave is mainly consumed for geometric spreading and plastic deformation in propagation in the elastoplastic medium model.

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.