Application of transient Rayleigh wave in detection of tunnel lining void area

Transient Rayleigh wave detection is a high-precision nondestructive detection method.At present,it has been widely used in shallow exploration,but rarely used in tunnel lining quality detection.Through the tunnel lining physical model experiment,the layout defects of the double-layer reinforcement lining area were detected and the Rayleigh wave velocity profile and dispersion curve were analyzed after data process-ing,which finally verified the feasibility and accuracy of Rayleigh wave method in detecting the tunnel lining void area.The results show that the method is not affected by the reinforcement inside the lining,the shallow detection is less disturbed and the accuracy is higher,and the data will fluctuate slightly with the deepening of the detection depth.At the same time,this method responds quite accurately to the thickness of the concrete,allowing for the assessment of the tunnel lining’s lack of compactness.This method has high efficiency,good reliability,and simple data processing,and is suitable for nondestructive detection of internal defects of tun-nel lining structure.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Rayleigh wave nonlinear inversion based on the Firefly algorithm

Rayleigh waves have high amplitude, low frequency, and low velocity, which are treated as strong noise to be attenuated in reflected seismic surveys. This study addresses how to identify useful shear wave velocity profile and stratigraphic information from Rayleigh waves. We choose the Firefly algorithm for inversion of surface waves. The Firefly algorithm, a new type of particle swarm optimization, has the advantages of being robust, highly effective, and allows global searching. This algorithm is feasible and has advantages for use in Rayleigh wave inversion with both synthetic models and field data. The results show that the Firefly algorithm, which is a robust and practical method, can achieve nonlinear inversion of surface waves with high resolution.

Dispersion function of Rayleigh waves in porous layered half-space system

Rayleigh wave exploration is based on an elastic layered half-space model. If practical formations contain porous layers, these layers need to be simplified as an elastic medium. We studied the effects of this simplification on the results of Rayleigh wave exploration. Using a half-space model with coexisting porous and elastic layers, we derived the dispersion functions of Rayleigh waves in a porous layered half-space system with porous layers at different depths, and the problem of transferring variables to matrices of different orders is solved. To solve the significant digit overflow in the multiplication of transfer matrices, we propose a simple, effective method. Results suggest that dispersion curves differ in a low- frequency region when a porous layer is at the surface; otherwise, the difference is small.

Dynamic response of utility tunnel during the passage of Rayleigh waves

The modeling methodologies and calculation of dynamic response of underground structure under Rayleigh waves is investigated in this paper. First the free field responses under Rayleigh waves are analyzed and the numerical results agree well with the theoretical results. Then, the approximate Rayleigh waves are put forward based on the preliminary re- search, and Rayleigh wave field is obtained through fast Fourier transform technique. Taking a utility tunnel as an example, its dynamic responses under Rayleigh waves is calculated by ABAQUS. The results demonstrate that bending deformation is the main component of structural deformation and the deformation at the top of the structure is about twice as much as that at bottom of the structure. The effect of soil-structure interface and the buried depth of underground structure are also investi- gated via parameter analysis. For the shallow buffed underground structures, Rayleigh waves can be the key factor to control the responses and damage of the structure.

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.

High-resolution Rayleigh wave phase velocity maps from ambient noise tomography in North China

We presented high-resolution Rayleigh wave phase velocity maps at periods ranging from 5 s to 30 s in the northeast part of the North China Craton (NNCC). Continuous time-series of vertical component between October 2006 and December 2008, recorded by 187 broadband stations temporarily deployed in the NNCC region, have been cross-correlated to obtain estimated fundamental mode Rayleigh wave Green’s functions. Using the frequency and time analysis technique based on continuous wavelet transformation, we measured 3 667 Rayleigh wave phase velocity dispersion curves. High-resolution phase velocity maps at periods of 5, 10, 20 and 30 s were reconstructed with grid size 0.25°× 0.25°, which reveal lateral heterogeneity of shear wave structure in the crust and upper mantle of NNCC. For periods shorter than 10 s, the phase velocity variations are well correlated with the principal geological units in the NNCC, with low-speed anomalies corresponding to the major sedimentary basins and high-speed anomalies coinciding with the main mountain ranges. Within the period range from 20 s to 30 s, high phase velocity observed in eastern NCC is coincident with the thin crust, whereas low phase velocities imaged in central NCC is correlated to the thick crust. However, the low-velocity anomaly in the Beijing-Tianjin-Tangshan region displayed in the 20 s and 30 s phase maps may be associated with fluids.

Nonlinear Rayleigh wave inversion based on the shuffled frog-leaping algorithm

At present, near-surface shear wave velocities are mainly calculated through Rayleigh wave dispersion-curve inversions in engineering surface investigations, but the required calculations pose a highly nonlinear global optimization problem. In order to alleviate the risk of falling into a local optimal solution, this paper introduces a new global optimization method, the shuffle frog-leaping algorithm （SFLA）, into the Rayleigh wave dispersion-curve inversion process. SFLA is a swarm-intelligence-based algorithm that simulates a group of frogs searching for food. It uses a few parameters, achieves rapid convergence, and is capability of effective global searching. In order to test the reliability and calculation performance of SFLA, noise-free and noisy synthetic datasets were inverted. We conducted a comparative analysis with other established algorithms using the noise-free dataset, and then tested the ability of SFLA to cope with data noise. Finally, we inverted a real-world example to examine the applicability of SFLA. Results from both synthetic and field data demonstrated the effectiveness of SFLA in the interpretation of Rayleigh wave dispersion curves. We found that SFLA is superior to the established methods in terms of both reliability and computational efficiency, so it offers great potential to improve our ability to solve geophysical inversion problems.

Group velocity distribution of Rayleigh waves and crustal and upper mantle velocity structure of the Chinese mainland and its vicinity

Based on the long period digital surface wave data recorded by 11 CDSN stations and 11 IRIS stations, the dispersion curves of the group velocities of fundamental mode Rayleigh waves along 647 paths, with the periods from 10 s to 92 s, were measured by multi-filter. Their distribution at 25 central periods within the region of 18~54N, 70~140E was inverted by Dimtar-Yanovskaya method. Within the period from 10 s to 15.9 s, the group velocity distribution is laterally inhomogeneous and is closely related to geotectonic units, with two low velocity zones located in the Tarim basin and the East China Sea and its north regions, respectively. From 21 s to 33 s, the framework of tectonic blocks is revealed. From 36.6 s to 40 s, the lithospheric subdivision of the Chinese mainland is obviously uncovered, with distinct boundaries among the South-North seismic belt, the Tibetan plateau, the North China, the South China and the Northeast China. Four cross-sections of group velocity distribution with period along 30N, 38N, 90E and 120E, are discussed, respectively, which display the basic features of the crust and upper mantle of the Chinese mainland and its neighboring regions. There are distinguished velocity differences among the different tectonic blocks. There are low-velocity-zones (LVZ) in the middle crust of the eastern Tibetan plateau, high velocity featured as stable platform in the Tarim basin and the Yangtze platform, shallow and thick low-velocity-zone in the upper mantle of the North China. The upper mantle LVZ in the East China Sea and the Japan Sea is related to the frictional heat from the subduction of the Philippine slab and the strong extension since the Himalayan orogenic period.

A hybrid inversion method of damped least squares with simulated annealing used for Rayleigh wave dispersion curve inversion

Surface wave methods are becoming increasingly popular in many geotechnical applications and in earthquake seismology due to their noninvasive characteristics.Inverse surface wave dispersion curves are a crucial step in most surface wave methods.Many inversion methods have been applied to surface wave dispersion curve inversion,including linearized inversion and nonlinearized inversion methods.In this study,a hybrid inversion method of Damped Least Squares(DLS) with Very Fast Simulated Annealing(VFSA) is developed for multi-mode Rayleigh wave dispersion curve inversion.Both synthetic and in situ fi eld data were used to verify the validity of the proposed method.The results show that the proposed method is superior to the conventional VFSA method in aiming at global minimum,especially when parameter searching space is adjacent to real values of the parameters.The advantage of the new method is that it retains both the merits of VFSA for global search and DLS for local search.At high temperatures,the global search dominates the runs,while at a low temperatures,the local search dominates the runs.Thus,at low temperatures,the proposed method can almost directly approach the actual model.

A method for inversion of layered shear wavespeed azimuthal anisotropy from Rayleigh wave dispersion using the Neighborhood Algorithm

Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm （NA） for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depth- dependent Vsv （or the elastic parameter L） as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polariza- tion azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.

All-parameters Rayleigh wave inversion

Since S-wave velocity of the subsurface is an important parameter in near surface applications,many studies have been conducted for its estimation.Among the various methods that use surface waves or body waves,Rayleigh wave inversion is the most popular.In practice,the densities and P-wave velocities of different layers are usually assumed to be known to avoid ill-posed problems,as they have less influence on the dispersion curves.However,improper assignment of these two groups of parameters leads to inaccurate estimation of the S-wave velocity profile.In order to address this problem,the all-parameters Rayleigh wave inversion strategy is proposed in which the S-wave velocities,layer thicknesses,densities and P-wave velocities of different layers are included as the unknown parameters for inversion.Meanwhile,the transitional Markov Chain Monte Carlo(TMCMC)algorithm is applied for the implementation of all-parameters Rayleigh wave inversion.One simulated example and two real-test applications are demonstrated to verify the capability of the proposed method in the estimation of the S-wave velocity profile,the densities and the P-wave velocities.Furthermore,it is verified that the proposed method achieved more accurate S-wave velocity profile estimation than the traditional approach.

Bayesian Rayleigh wave inversion with an unknown number of layers

Surface wave methods have received much attention due to their efficient, flexible and convenient characteristics. However, there are still critical issues regarding a key step in surface wave inversion. In most existing methods, the number of layers is assumed to be known prior to the process of inversion. However, improper assignment of this parameter leads to erroneous inversion results. A Bayesian nonparametric method for Rayleigh wave inversion is proposed herein to address this problem. In this method, each model class represents a particular number of layers with unknown S-wave velocity and thickness of each layer. As a result, determination of the number of layers is equivalent to selection of the most applicable model class. Regarding each model class, the optimization search of S-wave velocity and thickness of each layer is implemented by using a genetic algorithm. Then, each model class is assessed in view of its efficiency under the Bayesian framework and the most efficient class is selected. Simulated and actual examples verify that the proposed Bayesian nonparametric approach is reliable and efficient for Rayleigh wave inversion, especially for its capability to determine the number of layers.

The acoustoelastic effect on Rayleigh waves in elastic-plastic deformed layered rocks

On the basis of the acoustoelastic theory for elastic-plastic materials, the influence of statically deformed states including both the elastic and plastic deformations induced by applied uniaxial stresses on the Rayleigh wave in layered rocks is investigated by using a transfer matrix method. The acoustoelastic effects of elastic plastic strains in rocks caused by static deformations, are discussed in detail. The Rayleigh-type and Sezawa modes exhibit similar trends in acoustoelastic effect： the acoustoelastic effect increasing rapidly with the frequency-thickness product and the phase velocity change approaching a constant value for thick layer and high frequency limit. Elastic-plastic deformations in the Castlegate layered rock obviously modify the phase velocity of the Rayleigh wave and the cutoff points for the Sezawa modes. The investigation may be useful for seismic exploration, geotechnical engineering and ultrasonic detection.

An analytical study on the Rayleigh wave generation in a stratified structure

In this paper,an analytical method is used to investigate the Rayleigh wave generation in a stratified structure and the wave generation in a dry sandy layer constrained between the couple stress and inhomogeneous orthotropic half-spaces.This study is devoted to analyzing the impact of various effective parameters associated with the media on the phase velocities of the wave.The displacement components for each medium are derived by implementing the separable variable method.The frequency equa-tion is secured by using the displacement components in the boundary conditions,imposed at the interfaces between the layer and half-spaces.Moreover,the secured equation is the relation between the phase velocity and the wave number.Numerical computations are performed,and graphical representations are demonstrated between the phase velocity and the wave number for both phase velocities with different values of the parameters.The comparison between the phase velocities is observed for the same value of each pa-rameter.

Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials

Dispersion of the generalized Rayleigh waves propagating in a covered halfspace made of viscoelastic materials is investigated by utilizing the exact equations of the theory of linear viscoelasticity.The dispersion equation is obtained for an arbitrary type of hereditary operator of the materials of the constituents and a solution algorithm is developed for obtaining numerical results on the dispersion of the waves under consideration.Dispersion curves are presented for certain attenuation cases and the influence of the viscosity of the materials is studied through three rheological parameters of the viscoelastic materials which characterize the characteristic creep time,long-term values and the mechanical behaviour of the viscoelastic material around the initial state of the deformation.Numerical results are presented and discussed for the case where the viscoelasticity of the materials is described through fractional-exponential operators by Rabotnov.As the result of this discussion,in particular,how the rheological parameters influence the dispersion of the generalized Rayleigh waves propagating in the covered half-space under consideration is established.

Active-source Rayleigh wave dispersion by the Aki spectral formulation

Rayleigh wave imaging is efficient in estimating the shear- （S） wave velocity in near-surface exploration. The key is to accurately extract the dispersion of Rayleigh wave. We propose a method to calculate the dispersion of the active-source Rayleigh wavefield by using the Aki formulation. The spectrum after the cross correlation of two-channel records in the frequency domain is expressed by the Bessel function. Using the corresponding relation between the zero point of the spectrum real part and the Bessel function root, the phase velocity at the discrete frequency point is obtained and the dispersion curve is extracted. First, the theoretical basis and calculation method used in the active-source Rayleigh wave data are introduced. Then, three sets of theoretical models are calculated by this method and the results are consistent with the theoretical dispersion. Finally, we process a group of real Rayleigh wave data and obtain the phase velocity profiles and compared them with the results obtained by the multichannel surface wave analysis method. The effectiveness and applicability of the Aki method in active-source data processing are verified. Compared with multichannel wave processing, the advantage of the Aki method lies in the use of two- channel data in a single-shot record. When the number of acquisition channels in a shot gathers is insufficient or there is a bad channel, the quality of the extracted dispersion is guaranteed.

Effect of lateral heterogeneity on 2-D Rayleigh wave ZH ratio sensitivity kernels based on the adjoint method: Synthetic and inversion examples

The ratio between vertical and radial amplitudes of Rayleigh waves(hereafter,the Rayleigh wave ZH ratio)is an important parameter used to constrain structures beneath seismic stations.Some previous studies have explored crust and upper mantle structures by joint inversion of the Rayleigh wave ZH ratio and surface wave dispersion.However,all these studies have used a 1-D depth sensitivity kernel,and this kernel may lack precision when the structure varies a great deal laterally.Here,we present a systematic investigation of the two-dimensional(2-D)Rayleigh wave ZH ratio kernel based on the adjoint-wavefield method and perform two synthetic tests using the new kernel.The 2-D ZH ratio kernel is consistent with the traditional 1-D sensitivity kernel but has an asymmetric pattern with a preferred orientation toward the source.The predominant effect caused by heterogeneity can clearly be seen from kernels calculated from models with 2-D heterogeneities,which confirms the necessity of using the new 2-D kernel in some complex regions.Inversion tests using synthetic data show that the 2-D ZH ratio kernel has the potential to resolve small anomalies as well as complex lateral structures.

Numerical simulation for recognition of coalfield fire areas by Rayleigh waves

Effective recognition of a coalfield fire area improves fire-fighting efficiency and helps avoid potential geological hazards. Coalfield fire areas are hard to detect accurately using general geophysical methods. This paper describes simulations of shallow, buried coalfield fires based on real geological conditions. Recognizing the coalfield fire by Rayleigh wave is proposed. Four representative geological models are constructed, namely; the non-burning model, the pseudo-burning model, the real-burning model, and the hidden-burning model. Numerical simulation using these models shows many markedly different characteristics between them in terms of Rayleigh wave dispersion and Eigen displacement. These characteristics, as well as the shear wave velocity obtained by inverting the fundamental dispersion, make it possible to distinguish the type of the coalfield fire area and indentify the real and serious coalfield fire area. The results are very helpful for future application of Rayleigh waves for the detection of coalfield fire area.

Methods to increase the depth and precision of transient Rayleigh wave exploration

In order to increase the exploration depth of Rayleigh wave, new idea that dif-ferent from the former principles in data acquisition was applied. Suitable data acquisition parameter was given out on the basis of large amount of experiments. By reducing the group interval, the low frequency signal are enhanced instead of been attenuated. Fur-thermore, to solve the problem that the precision of Rayleigh wave exploration method count much to the signal-to-noise ratio, some preprocessing methods were put forward. By using zero shift rectifying, digital F-K filtering and cutting, noises can be effectively eliminated.