Effects of obstacles on Rayleigh wave dispersion obtained from the SASW test

This study investigates the effects of underground obstacles on a Rayleigh Wave dispersion curve obtained by spectral analysis of surface waves (SASW) testing. The SASW test is simulated by a transient response analysis of an axisymmetric finite element method. Both rigid obstacles and cavities with a variety of shapes and embedment depths are considered. Results show strong fluctuations in the dispersion curve obtained from the signals recorded by receivers in the vicinity of an obstacle. The fluctuations are primarily a result of wave reflections from the near and far faces of an obstacle. Results obtained from a finite element model are explained through a comparison with a simple theoretical model and experimental results are published. Fluctuations shown in the dispersion curve can result in misinterpretation of the shear wave velocity profile. Wherever obstacles are known to exist, an SASW test should be performed in a way that minimizes these effects. The parameters that will reduce or aggravate these effects were studied and are discussed.     

Inversion of SASW dispersion curves based on maximum flexibility coefficients in the wave number domain

Spectral analysis of surface waves (SASW) is a nondestructive in-situ testing method that is used to determine stiffness profiles of soil and pavement sites based on dispersion characteristics of Rayleigh-type surface waves.Inversion of the Rayleigh wave dispersion curve of a site provides information on the variation of shear-wave velocity with depth. In the inversion procedures currently used for SASW tests, the field dispersion curve is matched with a theoretical dispersion curve obtained for the fundamental mode of surface wave propagation.In order to overcome difficulties associated with the presence of multi-modes in SASW signals, a new inversion method based on the maximum vertical flexibility coefficient is introduced in this paper. Unlike root-searching methods, the new method easily identifies the predominant propagation modes. In this new approach, the simplex method is used to match field and theoretical dispersion curves automatically. The purpose of this paper is to present the details of the new method and to demonstrate its advantages.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

煤巷小构造Rayleigh型槽波超前探测数值模拟

对煤巷小构造地震波场进行了数值模拟研究,分析了层状煤层中地震波的传播特征.研究表明:(1)在煤巷迎头前方煤层内以纵波震源激发的Rayleigh型槽波相对于体波能量较强,波列较长,波速较低.(2)沿煤层传播的Rayleigh型槽波在小构造面上产生Rayleigh型槽波反射波,反射Rayleigh型槽波垂直分量相对于水平分量能量较强.沿煤层反向传播的反射Rayleigh型槽波在煤巷迎头面上转换为沿煤巷底板传播的Rayleigh面波.沿煤巷底板可以接收到能量较强的反射Rayleigh型槽波产生的Rayleigh面波,其可以作为超前探测小构造面的特征波.在地震记录上反射Rayleigh型槽波产生的Rayleigh面波波至最迟,在时间域与其他波列时间间隔较大,其垂直分量能量相对于水平分量较强,在地震记录上容易识别.(3)在相同的地质条件下应用反射地震超前探测方法,标志煤巷迎头前方存在小构造面的反射地震波能量较弱,受煤巷顶、底板界面和采煤迎头面的强反射波干扰,在地震记录中难以识别.     

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.     

On arrangement of source and receivers in SASW testing

This study investigates the effects of source and receivers arrangement on the Rayleigh wave dispersion curve in SASW testing. Analytical studies and numerical simulations with coupled finite and infinite elements are presented in this paper. It is shown that arrangement of source and receivers has a significant effect on test results, especially for soils with high Poisson's ratio or saturated soils. Larger source-to-receiver distance and receiver spacing usually give better results, and it is unnecessary to keep them equal. To satisfy the error control requirement in Rayleigh wave phase velocity measurement, source-to-receiver distance and receiver spacing should meet corresponding minimum values, which are proposed for different Poisson's ratios of soil in this paper.     

微动多阶瑞雷波SPAC系数反演方法及应用研究

为充分利用微动信号中基阶和高阶模式瑞雷波,本文研究了基于多阶瑞雷波SPAC系数直接反演的方法.该方法首先基于地层介质响应计算多阶瑞雷波的能量占比,考虑实际观测台阵有限台站个数对SPAC系数影响,正演计算多阶瑞雷波SPAC系数,再采用快速模拟退火算法对其反演以获得地下介质横波速度结构.在此基础上,本文通过数值模拟验证该方法的可靠性,分别选取三种典型地质模型,基于模式叠加算法合成理论微动信号,采用本文方法计算其理论多阶瑞雷波SPAC系数并反演,给出反演结果与真实模型对比.我们将该方法应用于上海中心城区的地质调查中,通过与钻探结果对比,进一步验证该方法的有效性.本文理论与实际应用研究表明,基于多阶瑞雷波SPAC系数直接反演的微动探测方法有助于提高反演结果的可靠性,尤其对含软硬夹层的复杂地层介质,可提高探测精度.     

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

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

Rayleigh波频散曲线“交叉”及多模式耦合作用研究

Rayleigh波可以用来反演近地表结构,在工程物探、石油物探、地球内部结构探测中均有重要意义.数值计算得到的含低速层的层状介质对应的Rayleigh波频散曲线会出现看似“交叉”的现象,但是对于这种现象目前还没有进行系统的研究.事实上可以验证,有些看似交叉的频散曲线实际上不相交.改变低速层的厚度和横波速度发现低速层越明显（即低速层速度越低或层厚越厚）频散曲线越不容易相交.凡友华等在2007年提出频散曲线对应着四种基本模式,在频散曲线发生“交叉”现象的区域实际上存在两个以上模式的频散曲线.本文主要研究了存在R模和S₂模的区域内频散曲线的“交叉”现象.首先利用竖直本征振动曲线研究R模和S2模Rayleigh 波的振动特点,发现R模对应的本征振动主要集中在地表,随着深度变化能量快速衰减,S₂模对应的本征振动主要集中在第2层.研究“交叉点”附近频散点对应的本征振动曲线发现这一区域有些Rayleigh波同时具有R模和S₂模的振动特点,对应着一种耦合模式.通过对实例的研究发现,在“交叉点”附近,若两条频散曲线不发生交叉,则每条曲线对应的模式会发生R模和S₂模之间经由耦合模式的转变,本文称这种现象为两种模式发生耦合;若两条频散曲线相交,则同一条频散曲线上的Rayleigh波模式几乎相同,只是在离交点很近的区域会存在一些耦合模式,本文称此时两种模式不发生耦合.本文研究结果主要供Rayleigh波对低速层结构的反演研究参考.     

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.     

场地瑞利波频散特性的测试方法研究

瑞利波频散特性测试是波速结构的面波探测方法中的关键环节。本文不进行速度结构的比较,而是直接与表面波谱分析方法比较所提取的频散曲线,从而将反演方法可能引入的不确定性因素剥离了出去,能够更客观地说明地脉动台阵方法的有效性和精度。在唐山26个场地进行的两种方法中美联合观测和分析结果表明。地脉动方法在83％的场地达到了与表面波谱分析方法相同的精度,是一种经济有效的场地瑞利波频散特性和波速结构测试方法。     

基阶与高阶瑞利波联合反演研究

研究了六层层状介质模型瑞利波基阶和高模式波相速度对横波速度、深度的敏感性,结果表明：基阶波较高模式波对7 m以内浅部地层的横波速度更敏感,敏感性频带在10～25 Hz范围内,峰值频带集中在18 Hz左右;高模式波较基阶波对深部地层的横波速度更敏感,敏感性频带宽,峰值分散.基阶波对浅层的敏感性和高模式波穿透深度更深的特点为近地表岩土层二维横波速度结构的联合反演提供了前提条件.利用阻尼最小二乘SVD（Singular Value Decomposition）算法联合基阶与高模式波对理论模型和实例数据进行横波速度反演,反演结果表明联合反演增强了反演的稳定性,提高了反演的精度.     

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

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

Simultaneous measurement and inversion of surface wave dispersion and attenuation curves

Surface wave tests are non-invasive seismic techniques that have traditionally been used to determine the shear wave velocity (i.e. shear modulus) profile of soil deposits and pavement systems. Recently, Rix et al. [J. Geotech. Geoenviron. Engng 126 (2000) 472] developed a procedure to obtain near-surface values of material damping ratio from measurements of the spatial attenuation of Rayleigh waves. To date, however, the shear wave velocity and shear damping ratio profiles have been determined separately. This practice neglects the coupling between surface wave phase velocity and attenuation that arises from material dispersion in dissipative media. This paper presents a procedure to measure and invert surface wave dispersion and attenuation data simultaneously and, thus, account for the close coupling between the two quantities. The methodology also introduces consistency between phase velocity and attenuation measurements by using the same experimental configuration for both. The new approach has been applied at a site in Memphis, TN and the results obtained are compared with independent measurements.     

Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm

Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity.     

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.     

瑞利波测试技术用于剪切波速勘探之研究

本文根据瑞利波测量所得波形的变化特征，分析研究了地基参数对瑞利波弥散曲线的影响规律，提出了利用瑞利波测试剪切波速的方法。在研究过程中，对瑞利波测试技术进行了探讨，并引入了相关分析原理，以消除噪声、系统的非线性以及其它一些人为因素的干扰。最后结合实例，论证该技术用于实际勘探的可行性并指出其存在的问题和有待进一步完善的工作。     

黏弹性与弹性介质中Rayleigh面波特性对比研究

Rayleigh面波的频散特性可以用来研究地表浅层结构. 本文使用时域有限差分法来模拟复杂黏弹性介质中的Rayleigh面波,研究了Q值对面波频散特性的影响.文中采用旋转交错网格有限差分,以非分裂卷积形式的完全匹配层为吸收边界,推出了求解二阶位移-应力各向同性黏弹性波动方程的数值方法.为了检验数值解的精度,首先将简单模型的正演结果与解析解对比,验证了方法的正确性;然后模拟了横向缓变层状介质和含有洞穴的介质中的面波,对弹性和黏弹性介质中的面波的频散特性进行对比分析.模拟结果表明浅层Q值对面波的频散特性有显著的影响;强吸收情况下,高阶面波的能量相对低阶面波能量显著增强.     

Analysis of dispersed multi-mode signals of the SASW method using the multiple filter/crosscorrelation technique

Spectral analysis of surface waves (SASW) is a nondestructive in-situ testing method used for determining the thickness and elastic properties of pavement and soil sites using the dispersion characteristics of surface waves. In previous studies, it has been demonstrated that for some sites errors may arise in experimental dispersion curves when the usual SASW test and data analysis procedures are followed, in particular the phase unwrapping procedure and source-to-near-receiver spacing distance. These errors occur due to the participation of more than one surface wave mode in SASW signals. In this study, the multiple filter/crosscorrelation technique often used in earthquake seismology for the analysis of multi-mode earthquake records is presented and applied for calculating phase velocities from SASW signals. It is demonstrated that this technique produces results that are generally more accurate than those produced by the usual phase unwrapping procedure.     

2015年8月12日天津化学爆炸产生的多模式面波分析及其应用研究

2015年8月12日天津滨海新区发生的强烈化学品爆炸造成了巨大的经济损失和社会影响.天津爆炸产生了清晰的大振幅面波信号,分析结果表明这组信号由基阶和高阶面波组成,可以追踪到约135 km外的远处台站.利用这组面波信号分别开展了以下研究：（1）利用附近三个台站记录的四个单频基阶Rayleigh波信号对爆破事件的绝对位置进行了网格搜索,结果与利用GPS测量的位置相差仅0.498 km;（2）分别利用网格搜索和主事件定位法,对两次子事件的相对位置进行了确定,距离约75 m左右,与前人研究结果吻合;（3）从面波记录中测量到36条基阶Rayleigh波、49条第一高阶Rayleigh波、9条基阶Love波和29条第一高阶Love波的频散曲线,并进一步反演获得研究区域地下4 km内的S波速度结构.反演结果显示地表处S波速度低至0.375 km·s^-1,在小于1 km的浅地表速度梯度较大,符合典型的盆地结构特征.本文的研究结果为类似爆炸等突发事件快速定位提供了新的思路,有助于灾后救援的迅速展开;同时得到天津滨海新区及周边浅层精细的速度结构,对于地震灾害评估有很大帮助.