Analysis of effects of active sources on observed phase velocity based on the thin layer method

In the wave field induced by active sources, the observed phase velocity of surface waves is influenced by both mode incompatibility (i.e. non-planar spread of surface waves is idealized as plane waves) and body waves. Effects of sources are usually investigated based on numerical simulations and physical models. Several methods have been proposed to mitigate the effects. In application, however, these methods may also have difficulties since the energy of the body waves depends on soil stratification and parameters. There are multiple modes of surface waves in layered media, among which the higher modes dominate the wave field for soils with the irregular shear velocity profiles. Considering the mode incompatibility and the higher modes, we derive analytical expressions for the effective phase velocity of the surface waves based on the thin layer stiffness method, and investigate the effects of the body waves on the observed phase velocity through the phase analysis of the vibrations of both the surface waves and the body waves. The results indicate that the effective phase velocity of the surface waves in layered media varies with the frequency and the spread distance, and is underestimated compared to that of the plane surface waves in the spread range less than about one wavelength. The oscillations that appeared in the observed phase velocity are due to the involvement of the body waves. The mode incompatibility can be ignored in the range beyond one wavelength, while the influence range of the body waves is far beyond one wavelength. The body waves have a significant influence on the observed phase velocity of the surface waves in soils with a soft layer trapped between the first and the second layers because of strong reflections.     

Signal reconstruction of surface waves on SASW measurement using Gaussian Derivative wavelet transform

Surface wave method consists of measurement and processing of the dispersive Rayleigh waves recorded from two or more vertical transducers. The dispersive phase data are inverted and the shear wave velocity versus depth is obtained. However, in case of residual soil, the reliable phase spectrum curve is difficult to be produced. Noises from nature and other human-made sources disturb the generated surface wave data. In this paper, a continuous wavelet transform based on mother wavelet of Gaussian Derivative was used to analyze seismic waves in different frequency and time. Time-frequency wavelet spectrum was employed to localize the interested seismic response spectrum of generated surface waves. It can also distinguish the fundamental mode of the surface wave from the higher modes of reflected body waves. The results presented in this paper showed that the wavelet analysis is able to determine reliable surface wave spectrum of sandy clayey residual soil.     

Dispersion Energy Analysis of Rayleigh and Love Waves in the Presence of Low-Velocity Layers in Near-Surface Seismic Surveys

High-frequency surface-wave analysis methods have been effectively and widely used to determine near-surface shear (S) wave velocity. To image the dispersion energy and identify different dispersive modes of surface waves accurately is one of key steps of using surface-wave methods. We analyzed the dispersion energy characteristics of Rayleigh and Love waves in near-surface layered models based on numerical simulations. It has been found that if there is a low-velocity layer (LVL) in the half-space, the dispersion energy of Rayleigh or Love waves is discontinuous and ‘‘jumping’’ appears from the fundamental mode to higher modes on dispersive images. We introduce the guided waves generated in an LVL (LVL-guided waves, a trapped wave mode) to clarify the complexity of the dispersion energy. We confirm the LVL-guided waves by analyzing the snapshots of SH and P–SV wavefield and comparing the dispersive energy with theoretical values of phase velocities. Results demonstrate that LVL-guided waves possess energy on dispersive images, which can interfere with the normal dispersion energy of Rayleigh or Love waves. Each mode of LVL-guided waves having lack of energy at the free surface in some high frequency range causes the discontinuity of dispersive energy on dispersive images, which is because shorter wavelengths (generally with lower phase velocities and higher frequencies) of LVL-guided waves cannot penetrate to the free surface. If the S wave velocity of the LVL is higher than that of the surface layer, the energy of LVL-guided waves only contaminates higher mode energy of surface waves and there is no interlacement with the fundamental mode of surface waves, while if the S wave velocity of the LVL is lower than that of the surface layer, the energy of LVL-guided waves may interlace with the fundamental mode of surface waves. Both of the interlacements with the fundamental mode or higher mode energy may cause misidentification for the dispersion curves of surface waves.     

Plane-layered models for the analysis of wave propagation in reservoir environments1

The long-wavelength propagation and attenuation characteristics of three geological structures that frequently occur in reservoir environments are investigated using a theoretical model that consists of a stack of fine and viscoelastic plane layers, with the layers being either solid or fluid. Backus theory properly describes fine layering and a set of fluid-filled microfractures, under the assumption that interfaces between different materials are bonded. The effects of saturation on wave attenuation are modelled by the relative values of the bulk and shear quality factors. The anisotropic quality factor in a fine-layered system shows a variety of behaviours depending on the saturation and velocities of the single constituents. The wave is less attenuated along the layering direction when the quality factors are proportional to velocity, and vice versa when inversely proportional to velocity. Fractured rocks have very anisotropic wavefronts and quality factors, in particular for the shear modes which are strongly dependent on the characteristics of the fluid filling the microfractures. When the size of the boundary layer is much smaller than the thickness of the fluid layer, the stack of solid-fluid layers becomes a layered porous media of the Biot type. This behaviour is caused by the slip-wall condition at the interface between the solid and the fluid. As in Biot theory, there are two compressional waves, but here the medium is anisotropic and the slow wave does not propagate perpendicular to the layers. Moreover, this wave shows pronounced cusps along the layering direction, like shear waves in a very anisotropic single-phase medium.     

倾斜边坡多道面波分析中的最小偏移距估计方法

为减小倾斜边坡中能量较强的近场体波对面波识别的干扰,确保多道面波分析方法接收到的瑞雷波分量具有较强的能量,本文通过分析地下震源在倾斜边坡产生的瑞雷波及其传播规律,基于几何地震学提出了在倾斜地表生成瑞雷波的最小偏移距的经验公式,建立了界面起伏的层状倾斜边坡模型,从而获得模拟共炮点记录,并将基于共炮点记录得到的地表质点运动...     

Techniques for mode separation in Rayleigh wave testing

The spectral analysis of surface wave (SASW) developed in the early eighties has opened the way to the use of surface waves for the definition of shear wave velocity profiles in soil deposits or pavement structures without the need of any borings or intrusion. The SASW testing procedure was designed to minimize the contribution of higher modes and thus assumes that the Rayleigh waves which propagate at the surface belong only to the fundamental mode. Several studies have however demonstrated that, in some conditions, higher Rayleigh modes can contribute significantly to the dispersion curve. Different tests configurations exist today to deal with Rayleigh mode problem by the use of an array of receivers. In spite of that, the SASW configuration remains attractive due to the limited number of receivers, as well as, the Rayleigh modes contributing in SASW records configuration can be identified by multiple-filter technique and isolated using time-variable filters. The proposed techniques are first validated by simulated records and then applied to SASW records obtained in the field. The study confirms that higher modes can participate and even dominate in SASW records. An important contribution of higher Rayleigh modes can also exist, even if the shear wave velocity increases regularly with depth. The higher Rayleigh modes can significantly affect the accuracy of the shear wave velocity profile if they are not properly identified and separated. A multi-mode inversion process is shown to be important to have an accurate soil characterization.     

Forward modeling of Rayleigh surface waves for analytical characterization of dominant dispersion trends

Forward modeling is of critical importance for inversion analysis of surface wave methods to obtain shear-wave velocity (V_S) profiles of soil sites. The dynamic stiffness matrix (DSM) method can provide forward modeling of Rayleigh surface waves to simulate complex wave propagation in layered soil sites. However, contamination from body waves and interference of multiple Rayleigh wave modes can reduce the accuracy of theoretical dispersion curves, especially at irregular soil sites with embedded low-velocity or high-velocity layers. An analytical method is developed herein to combine the techniques of the multichannel analysis of surface waves method with the DSM method to improve the accuracy of the theoretical dispersion analysis for soil sites. The proposed method implements multichannel analysis of the analytical displacement responses to capture dominant dispersion trends. Comparison of the results obtained with the new method against those from the transfer matrix method and the literature indicates that the new method can (1) effectively minimize the effects of contamination caused by body waves and interference from several Rayleigh wave modes, and (2) generate accurate dominant dispersion trends for soil sites with various stiffness profiles, especially for the high-frequency dispersion characteristics of the profiles with embedded low-velocity layers.     

浅层有限频率面波成像中的3D灵敏度核分析

本文利用面波散射的模式耦合方法,基于波恩近似和远场假设,研究了有限频率面波三维灵敏度核,针对面波在工程应用中常遇到的水平分层的背景介质模型,计算了介质扰动引起的面波相位和幅度扰动的三维灵敏度核,分析了模式耦合对三维灵敏度核的影响.结果表明,仅考虑模式自身耦合的JWKB近似,介质密度和波速扰动引起的三维灵敏度核可以蜕化为面波相速度扰动的二维灵敏度核,有限频率面波成像和传统基于射线的方法类似,可以分两步进行,只需在纯路径反演过程中将基于大圆路径假设下的一维灵敏度核用二维灵敏度核取代.如果交叉模式耦合的影响不能忽略,在反演时必须引入三维面波灵敏度核,直接对介质参数进行反演.     

裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2—3kHz）和长源距（3—4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2-3kHz）和长源距（3-4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

On the propagation characteristics of tunnel surface‐waves for seismic prediction

To increase the safety and efficiency of tunnel constructions, online seismic exploration ahead of a tunnel has become a valuable tool. One recent successful forward looking approach is based on the excitation and registration of tunnel surface‐waves. For further development and for finding optimal acquisition geometries it is important to study the propagation characteristics of tunnel surface‐waves. 3D seismic finite difference modelling and analytic solutions of the wave equation in cylindrical coordinates reveal that at higher frequencies, i.e., if the tunnel‐diameter is significantly larger than the wavelength of surface‐waves, these surface‐waves can be regarded as Rayleigh‐waves confined to the tunnel wall and following helical paths along the tunnel axis. For lower frequencies, i.e., when the tunnel surface‐wavelength approaches the tunnel‐diameter, the propagation characteristics of these surface‐waves are similar to S‐waves. We define the surface‐wave wavelength‐to‐tunnel diameter ratio w to be a gauge for separating Rayleigh‐ from S‐wave excitation. For w > 1.2 tunnel surface‐waves behave like S‐waves, i.e. their velocity approaches the S‐wave velocity and the particle motion is linear and perpendicular to the ray direction. For w < 0.6 they behave like Rayleigh‐waves, i.e., their velocity approaches the Rayleigh‐wave velocity and they exhibit elliptical particle motion. For 0.6 < w < 1.2 a mixture of both types is observed. Field data from the Gotthard Base Tunnel (Switzerland) show both types of tunnel surface‐waves and S‐waves propagating along the tunnel.     

Effect of surface stress on magneto-elastic surface waves in finitely conducting media

This paper deals with the investigation of the effect of surface stress and conductivity on the propagation of surface wave in isotropic, homogeneous, elastic media under the action of a primary magnetic field. Formulation of the general surface wave propagation problem has been made, and the corresponding frequency equation has been derived. Frequency equations for Rayleigh wave, surface shear wave and Stoneley wave have been deduced from that of general surface wave as special cases. The effects of surface stress and magnetic field on the wave velocities and attenuation factors of Rayleigh wave and surface shear wave are shown by numerical calculations and graphs. Some important wave velocity equations, as obtained by other authors, have been deduced as special cases from the wave velocity equation for Stoneley wave. It is found that the combined effect of surface stress and magnetic field modulates the wave velocity ratios and attenuation factors of Rayleigh wave and surface shear wave to a considerable extent.     

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity–stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10–16% when the frequency is above 10 Hz due to the velocity dispersion of P and S waves.     

Simple equations guide high-frequency surface-wave investigation techniques

We discuss five useful equations related to high-frequency surface-wave techniques and their implications in practice. These equations are theoretical results from published literature regarding source selection, data-acquisition parameters, resolution of a dispersion curve image in the frequency–velocity domain, and the cut-off frequency of high modes. The first equation suggests Rayleigh waves appear in the shortest offset when a source is located on the ground surface, which supports our observations that surface impact sources are the best source for surface-wave techniques. The second and third equations, based on the layered earth model, reveal a relationship between the optimal nearest offset in Rayleigh-wave data acquisition and seismic setting—the observed maximum and minimum phase velocities, and the maximum wavelength. Comparison among data acquired with different offsets at one test site confirms the better data were acquired with the suggested optimal nearest offset. The fourth equation illustrates that resolution of a dispersion curve image at a given frequency is directly proportional to the product of a length of a geophone array and the frequency. We used real-world data to verify the fourth equation. The last equation shows that the cut-off frequency of high modes of Love waves for a two-layer model is determined by shear-wave velocities and the thickness of the top layer. We applied this equation to Rayleigh waves and multi-layer models with the average velocity and obtained encouraging results. This equation not only endows with a criterion to distinguish high modes from numerical artifacts but also provides a straightforward means to resolve the depth to the half space of a layered earth model.     

高频面波方法的若干新进展

面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中, 由于该方法具有非侵入性、无损、高效及价格低的特点, 越来越受到浅地表地球物理和地质工程学界的重视, 视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质, 尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量, 从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少, 这使得勒夫波的频散曲线比瑞雷波的简单.因此, 勒夫波反演更稳定, 非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰, 并具有更高的分辨率, 从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言, 瑞雷波的探测深度是勒夫波的1.3~1.4倍; 而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体, 利用反卷积消除了地震子波对勒夫波波形的影响, 通过更新每个块体的S波速度来拟合勒夫波波形, 从而获得地下S波速度模型.该方法不基于水平层状模型假设, 适用于任意二维介质模型.     

基于萤火虫算法的雷瑞波非线性反演(英文)

雷瑞波具有强振幅、低频和低速的特点,在反射地震勘探中通常是需要被压制的强噪声。本文研究如何利用雷瑞波获取近地表地层的横波速度和地下结构,选取萤火虫优化算法进行面波的反演,萤火虫优化算法是一种新的粒子群算法理论,具有稳定、快捷、全局搜索等特点。针对萤火虫优化算法优缺点进行了讨论和改进,通过对理论模型和野外数据的测试应用,将提取的瑞利面波频散曲线反演得到横波速度信息。结果表明萤火虫优化算法能实现面波非线性反演,并具有分辨率高、抗干扰能力强等优点和实际使用前景。     

Analytical expressions for the borehole seismic signatures excited by an explosive source buried in an elastic half‐space

Borehole guided waves that are excited by explosive sources outside of the borehole are important for interpreting borehole seismic surveys and for rock property inversion workflows. Borehole seismograms are typically modelled using numerical methods of wave propagation. In order to benchmark such numerical algorithms and partially to interpret the results of modelling, an analytical methodology is presented here to compute synthetic seismograms. The specific setup is a wavefield emanating from a monopole point source embedded within a homogeneous elastic medium that interacts with a fluid‐filled borehole and a free surface. The methodology assumes that the wavelength of the seismic signal is much larger than the borehole radius. In this paper, it is supposed that there is no poroelastic coupling between the formation and the borehole. The total wavefield solution consists of P, PP, and PS body waves; the surface Rayleigh wave; and the low‐frequency guided Stoneley wave (often referred as the tube wave) within the borehole. In its turn, the tube wave consists of the partial responses generated by the incident P‐wave and the reflected PP and PS body waves at the borehole mouth and by the Rayleigh wave, as well as the Stoneley wave eigenmode. The Mach tube wave, which is a conic tube wave, additionally appears in the Mach cone in a slow formation with the tube‐wave velocity greater than the shear one. The conditions of appearance of the Mach wave in a slow formation are formulated. It is shown that the amplitude of the Mach tube wave strongly depends on Poisson's ratio of the slow surrounding formation. The amplitude of the Mach tube wave exponentially decreases when the source depth grows for weakly compressible elastic media with Poisson's ratio close to 0.5 (i.e., saturated clays and saturated clay soils). Asymptotic expressions are also provided to compute the wavefield amplitudes for different combinations of source depth and source‐well offset. These expressions allow an approximate solution of the wavefield to be computed much faster (within several seconds) than directly computing the implicit integrals arising from the analytical formulation.     

土中Love波弥散特性

利用有限单元法及解析法建立和求解了土中Love波特征方程以及位移计算公式．计算结果表明,这一计算方法比纯解析法优越,可以用来分析均质和非均质上中Love波弥散性．本文利用这一方法详细讨论了Love波在上软下硬地基及软夹层地基中的传播特性和弥散特性．上软下硬地基Love波具有弥散性,土层的剪切波及厚度对Love波弥散曲线影响较大,而质量密度的相对变化对Love彼弥散曲线影响较小．软夹层地基中低频时Love波以第一模态波为主,现场所测为第一模态波波速;高频时存在多个高模态波,土中传播的波为这几个高模态波的叠加波,现场所测波速随两传感器的位置不同而有波动．     

黏弹性介质瑞雷波有限差分模拟与特性分析

本文通过数值模拟研究了介质黏弹性对瑞雷波传播的影响.模拟采用结合了交错Adams-Bashforth时间积分法、应力镜像法和多轴完美匹配层的标准交错网格高阶有限差分方案.通过模拟结果和理论结果对比,测试了方法的精度,验证了结果的正确性.在均匀半空间模型中,分别从波场快照、波形曲线及频散能量图三个角度,对黏弹性介质瑞雷波衰减和频散特性进行了详细分析.两层速度递增模型被用于进一步分析瑞雷波在黏弹性层状介质中的特性.结果表明：由于介质的黏弹性,瑞雷波振幅发生衰减,高频成分比低频成分衰减更剧烈,衰减程度随偏移距增大而增强;瑞雷波相速度发生频散,且随频率增大而增大,频散能量的分辨率有所降低;黏弹性波动方程中的参考频率,不会影响瑞雷波振幅衰减和相速度频散的程度,但决定了黏弹性和弹性介质瑞雷波相速度相等的频率位置.本研究有助于人们更好地理解地球介质中瑞雷波的行为,并为瑞雷波勘探的应用和研究提供了科学和有价值的参考.