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.

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.

SegNet-based first-break picking via seismic waveform classification directly from shot gathers with sparsely distributed traces

Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to address this issue and present a fast automatic seismic waveform classification method to pick densely-sampled FBs directly from common-shot gathers with sparsely distributed traces.Through feeding a large number of representative shot gathers with missing traces and the corresponding binary labels segmented by manually interpreted fully-sampled FBs,we can obtain a welltrained Seg Net model.When any unseen gather including the one with irregular trace spacing is inputted,the Seg Net can output the probability distribution of different categories for waveform classification.Then FBs can be picked by locating the boundaries between one class on post-FBs data and the other on pre-FBs background.Two land datasets with each over 2000 shots are adopted to illustrate that one well-trained 25-layer Seg Net can favorably classify waveform and further pick fully-sampled FBs verified by the manually-derived ones,even when the proportion of randomly missing traces reaches50%,21 traces are missing consecutively,or traces are missing regularly.

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.

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.

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.

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.

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.

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.

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

在多分量地震勘探资料中,水平和垂直分量都记录有P波和SV波。本文研究了P波和SV波波场分解的方法。当上行P波和S波分别入射时,通过研究它们引起的水平和垂直位移,把地震记录的水平和垂直分量进行分解,得到纵波和转换波;在已知海底介质中纵横波速度和介质密度的情况下,可以在τ-P域内实现波场分解,然后变换到时间域,得到时间域内的纵波和转换波剖面。把本文研究的波场分解方法应用于合成资料,能够有效地分解得到纵波波场和转换波波场。最后,海上多分量实际资料的实例应用表明,本文研究的波场分解方法是可行有效的。本方法也适用于自由表面的资料。

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.

Seismic Response of a Typical Fixed Jacket-Type Offshore Platform (SPD1) Under Sea Waves

Offshore platforms in seismically active areas should be designed to service severe earthquake excitations with no global structural failure. In seismic design of offshore platforms, it is often necessary to perform a dynamic analysis that accounts for nonlinear pile soil structures interaction effects. This paper summarizes the nonlinear dynamic analysis of a 3-D model of a typical Jacket-Type platform which is installed in Persian Gulf (SPD1), under simultaneously wave and earthquake loading has been conducted. It is assumed that they act in the same and different directions. The interaction between soil and piles is modeled by equivalent pile length theory. The structure is modeled by finite element method (Ansys Inc.). It be concluded that when the longitudinal components of the earthquake and wave are in different directions, an increase on the response of platform can be seen.

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 scheme to treat the singularity in global seismic wavefield simulation using pseudospectral method with staggered grids

The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.

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 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.

Seismic Wave Propagation and Excitation in Multi-layered Media with Irregular Interfaces Part (Ⅰ): SH Case

Study of seismic wave excitation and propagation in laterally heterogeneous media was an active and important subject in seismology in the past two decades, numerous analytical and numerical efforts have been made in this research field. In this article, I have, first, made a brief review on those developments and then introduced and summarized a unified and efficient method, global generalized reflection-transmission (abbreviated to R/T thereafter) matrices method, for synthetic seismograms in multi-layered media with irregular interfaces developed by the author [24~26]. As demonstrated in this article, this method could be regarded as an extension of the generalized R/T coefficients method for the horizontally layered case [2,5] to the layered media with irregularly shaped interfaces by incorporating the T matrices technique [27,28]. Because of the use of a recursive scheme in computing the global generalized R/T matrices, this method is efficient, particularly for the case with a large number of

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 rock slope.The analysis presented here provides insight into the dynamic response of the "Dualistic" 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" structure 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.