Seismic anisotropy and upper mantle dynamics in Alaska:A review of shear wave splitting analyses

Shear wave splitting(SWS)is regarded as the most effective geophysical method to delineate mantle flow fields by detecting seismic azimuthal anisotropy in the earth's upper mantle,especially in tectonically active regions such as subduction zones.The Aleutian-Alaska subduction zone has a convergence rate of approximately 50 mm/yr,with a trench length reaching nearly 2800 km.Such a long subduction zone has led to intensive continental deformation and numerous strong earthquakes in southern and central Alaska,while northern Alaska is relatively inactive.The sharp contrast makes Alaska a favorable locale to investigate the impact of subduction on mantle dynamics.Moreover,the uniqueness of this subduction zone,including the unusual subducting type,varying slab geometry,and atypical magmatic activity and composition,has intrigued the curiosity of many geoscientists.To identify different sources of seismic anisotropy beneath the Alaska region and probe the influence of a geometrically varying subducting slab on mantle dynamics,extensive SWS analyses have been conducted in the past decades.However,the insufficient station and azimuthal coverage,especially in early studies,not only led to some conflicting results but also strongly limited the in-depth investigation of layered anisotropy and the estimation of anisotropy depth.With the completion of the Transportable Array project in Alaska,recent studies have revealed more detailed mantle structures and characteristics based on the dense station coverage and newly collected massive seismic data.In this study,we review significant regional-and continental-scale SWS studies in the Alaska region and conclude the mantle flow fields therein,to understand how a geometrically varying subducting slab alters the regional mantle dynamics.The summarized mantle flow mechanisms are believed to be conducive to the understanding of seismic anisotropy patterns in other subduction zones with a complicated tectonic setting.

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.

Accurate simulations of pure-viscoacoustic wave propagation in tilted transversely isotropic media

Forward modeling of seismic wave propagation is crucial for the realization of reverse time migration(RTM) and full waveform inversion(FWI) in attenuating transversely isotropic media. To describe the attenuation and anisotropy properties of subsurface media, the pure-viscoacoustic anisotropic wave equations are established for wavefield simulations, because they can provide clear and stable wavefields. However, due to the use of several approximations in deriving the wave equation and the introduction of a fractional Laplacian approximation in solving the derived equation, the wavefields simulated by the previous pure-viscoacoustic tilted transversely isotropic(TTI) wave equations has low accuracy. To accurately simulate wavefields in media with velocity anisotropy and attenuation anisotropy, we first derive a new pure-viscoacoustic TTI wave equation from the exact complex-valued dispersion formula in viscoelastic vertical transversely isotropic(VTI) media. Then, we present the hybrid finite-difference and low-rank decomposition(HFDLRD) method to accurately solve our proposed pure-viscoacoustic TTI wave equation. Theoretical analysis and numerical examples suggest that our pure-viscoacoustic TTI wave equation has higher accuracy than previous pure-viscoacoustic TTI wave equations in describing q P-wave kinematic and attenuation characteristics. Additionally, the numerical experiment in a simple two-layer model shows that the HFDLRD technique outperforms the hybrid finite-difference and pseudo-spectral(HFDPS) method in terms of accuracy of wavefield modeling.

A review of the wave gradiometry method for seismic imaging

As dense seismic arrays at different scales are deployed,the techniques to make full use of array data with low computing cost become increasingly needed.The wave gradiometry method(WGM)is a new branch in seismic tomography,which utilizes the spatial gradients of the wavefield to determine the phase velocity,wave propagation direction,geometrical spreading,and radiation pattern.Seismic wave propagation parameters obtained using the WGM can be further applied to invert 3D velocity models,Q values,and anisotropy at lithospheric(crust and/or mantle)and smaller scales(e.g.,industrial oilfield or fault zone).Herein,we review the theoretical foundation,technical development,and major applications of the WGM,and compared the WGM with other commonly used major array imaging methods.Future development of the WGM is also discussed.

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.

Weighted-elastic-wave interferometric imaging of microseismic source location

Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging method.

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.

Separation of P. and SV-wavefields from multi-componen seismic data

In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.

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.

Contribution of tuned liquid column gas dampers to the performance of offshore wind turbines under wind, wave, and seismic excitations

The main intention of the present study is to reduce wind, wave, and seismic induced vibrations of jacket- type offshore wind turbines （JOWTs） through a newly developed vibration absorber, called tuned liquid column gas damper （TLCGD）. Using a Simulink-based model, an analytical model is developed to simulate global behavior of JOWTs under different dynamic excitations. The study is followed by a parametric study to explore efficiency of the TLCGD in terms of nacelle acceleration reduction under wind, wave, and earthquake loads. Study results indicate that optimum frequency of the TLCGD is rather insensitive to excitation type. In addition, while the gain in vibration control from TLCGDs with higher mass ratios is generally more pronounced, heavy TLCGDs are more sensitive to their tuned frequency such that ill-regulated TLCGD with high mass ratio can lead to destructive results. It is revealed that a well regulated TLCGD has noticeable contribution to the dynamic response of the JOWT under any excitation.

Seismic stability of jointed rock slopes under obliquely incident earthquake waves

Seismic stability of slopes has been traditionally analyzed with vertically propagated earthquake waves.However,for rock slopes,the earthquake waves might approach the outcrop still with a evidently oblique direction.To investigate the impact of obliquely incident earthquake excitations,the input method for SV and P waves with arbitrary incident angles is conducted,respectively,by adopting the equivalent nodal force method together with a viscous-spring boundary.Then,the input method is introduced within the framework of ABAQUS software and verified by a numerical example.Both SV and P waves input are considered herein for a 2 D jointed rock slope.For the jointed rock mass,the jointed material model in ABAQUS software is employed to simulate its behavior as a continuum.Results of the study show that the earthquake incident angles have significance on the seismic stability of jointed rock slopes.The larger the incident angle,the greater the risk of slope instability.Furthermore,the stability of the jointed rock slopes also is affected by wave types of earthquakes heavily.P waves induce weaker responses and SV waves are shown to be more critical.

Characteristics of the Seismic Waves from a New Active Source Based on Methane Gaseous Detonation

Active seismic sources are critical for obtaining high resolution images of the subsurface.For active imaging in urban areas,environment friendly and green seismic sources are required.In present work,we introduce a new type of green active source based on the gaseous detonation of methane and oxygen.When fired in a closed container,the chemical reaction,i.e.gaseous detonation,will produce high pressure air over 150 MPa.Seismic waves are produced when high pressure air is quickly released to impact the surroundings.The first field experiment of this active source was carried out in December,2017 in Jingdezhen,Jiangxi Province,where a series of active sources were excited to explore their potential in mine exploration.In current work,we analyzed the seismic waves recorded by near-field accelerators and a dense short-period seismic array and compared them with those from a mobile airgun source,another kind of active source by releasing high pressure air into water.The results demonstrate that it can be used for high resolution near surface imaging.Firstly,the gaseous detonation productions are harmless CO2 and water,making it a green explosive source.Secondly,the dominant seismic frequencies are 10-80 Hz and a single shot can be recorded up to 15 km,making it suitable for local structure investigations.Thirdly,it can be excited in vertical wells,similar to traditional powder explosive sources.It can also act as an additional on-land active source to airgun sources,which requires a suitable water body as intermediate media to generate repeating signals.Moreover,the short duration and high frequency signature of the source signals make it safe with no damage to nearby buildings.These make it convenient to excite in urban areas.As a new explosive source,the excitation equipment and conditions,such as gas ratio,sink depth and air-releasing directions,need further investigation to improve seismic wave generation efficiency.

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.

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.