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.

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.

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.

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.

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.

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.

Studies on seismic waves

The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous media, anisotropic and porous media, surface wave and seismic wave inversion, and seismic wave study in prospecting and logging problems. Important projects in the current studies on seismic wave is suggested as the development of high efficient numerical methods, and applying them to the studies of excitation and propagation of seismic waves in complex media and strong ground motion, which will form a foundation for refined earthquake hazard analysis and prediction.

Present State of Explosion Seismic Wave Research and Primary Investigation on Its Characteristics

The present state and the significance of research on explosion seismic waves are discussed, and meanwhile the main contents and the basic problems to be solved in the study of explosion seismic waves are analyzed. The spectra characteristics of explosion seismic waves, functions of the isolated-seismic grooves and influences of the detonating methods on explosion seismic waves are investigated by experiments. The experimental method is introduced. Some experimental results are presented which are concerned with the influences of topographical conditions, explosive charges, ignition patterns, isolated-seismic grooves and the other related factors on the characteristics of seismic waves.

A laboratory method to simulate seismic waves induced by underground explosions

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

SIMULATION SEISMIC WAVE PROPAGATION IN TOPOGRAPHIC STRUCTURES USING ASYMMETRIC STAGGERED GRIDS

A new 3D finite-difference(FD) method of spatially asymmetric staggered grids was presented to simulate elastic wave propagation in topographic structures.The method approximated the first-order elastic wave equations by irregular grids finite difference operator with second-order time precise and fourth-order spatial precise.Additional introduced finite difference formula solved the asymmetric problem arisen in non-uniform staggered grid scheme.The method had no interpolation between the fine and coarse grids.All grids were computed at the same spatial iteration.Complicated geometrical structures like rough submarine interface,fault and nonplanar interfaces were treated with fine irregular grids.Theoretical analysis and numerical simulations show that this method saves considerable memory and computing time,at the same time,has satisfactory stability and accuracy.

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.

A micro-kinetic model of enhanced foam stability under artificial seismic wave

To get a deeper understanding on the synergistic enhancement effect of low frequency artificial seismic wave on foam stability,a micro-kinetic model of enhanced foam stability under low frequency artificial seismic wave is established based on a vertical liquid film drainage model and elastic wave theory.The model is solved by non-dimensional transformation of the high order partial differential equations and a compound solution of implicit and explicit differences and is verified to be accurate.The foam film thickness,surfactant concentration distribution and drainage velocity under the action of low frequency artificial seismic wave are quantitatively analyzed.The research shows that low-frequency vibration can reduce the difference between the maximum and minimum concentrations of surfactant in the foam liquid film at the later stage of drainage,enhance the effect of Marangoni effect,and improve the stability of the foam liquid film.When the vibration frequency is close to the natural frequency of the foam liquid film,the vibration effect is the best,and the best vibration frequency is about 50 Hz.The higher the vibration acceleration,the faster the recovery rate of surfactant concentration in the foam liquid film is.The higher the vibration acceleration,the stronger the ability of Marangoni effect to delay the drainage of foam liquid film and the better the foam stability is.It is not the higher the vibration acceleration,the better.The best vibration acceleration is about 0.5 times of gravity acceleration.Reasonable vibration parameters would greatly enhance the effect of Marangoni effect.The smaller the initial concentration of surfactant,the better the vibration works in enhancing Marangoni effect.

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.

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

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.