基于互相关相移的主动源地震面波频散成像方法

以多道面波分析技术(MASW)为代表的主动源面波方法已经在地表勘探领域获得越来越多的应用.准确的面波相速度频散成像是利用主动源面波法进行近地表横波速度反演的必要条件.本文提出了一种新的主动源面波频散成像方法—CCPS(Cross-Correlation and Phase-Shifting)法.首先,将单炮面波记录由时空域变换到频率空间域并消去振幅谱的影响,得到归一化频谱.然后,依次将每道信号与整个记录进行互相关的同时进行相速度扫描叠加.最后经过一次求和过程即得到最终的频散图像.通过理论数据、实际工程地震面波资料和石油勘探面波资料比较了普通相移法和CCPS法的频散成像效果.研究表明:新方法相对于普通相移法而言,明显提高了频散图像的质量,相速度聚焦性能和背景噪声压制能力有显著改善.     

用面波联合勘探技术探测浅部速度结构

瑞雷面波勘探技术以其快速经济、受场地条件限制小等优点广泛应用于浅部横波速度结构探测.人工源面波勘探方法对浅部地层的探测精度高,但探测深度较浅;天然源面波勘探方法探测深度较深,但对浅部速度结构的探测精度不高.本文在夏垫和玉溪两地分别开展了人工源和天然源面波联合探测试验,尝试采用不同排列和相同排列两种方法,采集人工源和天然源面波信息,联合处理数据并提取频散曲线,反演得到浅部地层的横波速度结构.探测结果表明:人工源和天然源面波联合勘探,尤其是采用相同排列的方法,可以在几乎不增加常规面波勘探工作量的条件下,既能保证浅部地层的探测精度,又明显拓展探测深度,大大提升了面波勘探能力,有望在工程勘察领域中推广应用.     

基于面波频散曲线聚类分析的近地表横波速度反演

获取准确的近地表横波速度对复杂地表条件下弹性波地震数据处理和成像非常重要.在浅层面波工程勘探中通过反演提取的频散曲线可以获得近地表横波速度结构.在多道面波频散曲线分析中,频散关系拾取的精度直接影响速度反演结果的可靠性.本文在多道面波叠加及自动拾取频散曲线基础上,提出了基于面波频散曲线聚类分析的近地表横波速度反演方法.该方法充分考虑了低信噪比条件下面波频散曲线的不确定性,通过在频散曲线拾取中引入曼哈顿距离K-Means聚类算法提高频散曲线拾取的准确性.采用多道多窗口叠加技术提高了面波反演对横向速度变化的适应性,通过聚类算法和多窗口叠加提高反演的可靠性,聚类算法获得较准确的频散曲线更利于后续横波速度反演过程.模拟数据算例对比表明本文提出的方法比常规算法效果更好,精度更高.将提出的方法应用于工程勘探和油气勘探的面波数据反演中,结果也验证了该方法的有效性.     

埋入源多道面波分析(MASW)中最小偏移距的估计方法

多道面波分析（MASW）法是一种被广泛使用的浅层地震勘探方法.为解决埋入源近场效应对MASW法应用的影响,本文分析了埋入源产生的Rayleigh波传播机制,探讨了MASW法的排列参数对地震记录的影响,提出了确定最小偏移距的经验公式.通过对实际地质模型的有限元模拟,验证了经验公式的合理性.在以经验公式为指导的野外实测中,利用多种方法提取了相对可靠的频散曲线,采用基阶和高阶面波联合反演,所得的一维速度结构与钻孔资料具有很强的相关性,表明本文发展的最小偏移距估计方法具有一定的实际应用价值.     

面波频散谱多模式高分辨率成像的多道信号比较法

高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,同时极大地提高了频散谱的成像分辨率.然而,在研究中我们发现,仅利用两个台站的地震记录,并不能将面波的频散特性完整地考虑在内,导致LSC和NLSC方法对高阶模式的成像存在较大的误差.由于主动源面波勘探多道采集的方式,基于信号比较理论的多道信号比较法(MSC)充分利用多道地震信号,可以获得准确的多模式成像,然而该方法需要计算任意两道的频散谱并叠加,存在冗余的计算,导致计算效率较低.因此,本研究对MSC方法进行了相应的改进,通过追踪炮集记录上的面波波组提高了原方法的计算效率,同时,利用理论频散曲线进行叠加分析,验证了改进的MSC方法的正确性和有效性.通过与相移法、LSC和NLSC方法的对比分析,展示了MSC方法是一种准确的、高分辨率的面波多模式频散谱成像方法.实际地震资料的应用,揭示了MSC方法在浅层地震勘探中用于提取面波多模式频散信息的巨大潜力.     

基于矢量波数变换法(VWTM)的多道Rayleigh波分析方法

在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息.     

基于改进Frequency-Bessel变换法的主动源Love波频散分析

近年来,Love波多道分析技术在浅地表结构探测中受到越来越广泛的应用.与Rayleigh波相比,Love波相速度不受纵波速度影响,其反演参数较少,可使反演过程更加稳定、求解的横波速度模型更加可靠.高分辨率和多模式面波频散分析是面波多道分析技术中至关重要的一环.本文利用改进Frequency-Bessel变换法对主动源Love波进行频散分析.通过改进的0阶Frequency-Bessel变换,将时间-空间域多道Love波记录变换到频率-波数域,获得其频散能量谱.通过公式推导和数值模拟,验证了该方法的有效性.实例测试表明,该方法具有较高的分辨率和多模式分辨能力,为多分量面波勘探提供了一种有效的频散分析手段.     

长三角地区地壳P、S波速度结构模型

利用收集到的瑞利面波频散数据和区域地震台网震相数据,分别采用面波层析成像和体波走时层析成像方法反演了长三角地区地下速度异常分布,获得了该区地壳P、S波速度结构.层析成像结果表明速度结构模型与地质特征基本一致,研究结果有助于深化对长三角地区沉积盖层和地壳结构以及强震活动深部环境的认识.     

多道面波分析方法及其应用研究

正浅层地表介质的剪切波速是工程抗震设计、场地地震危险性分析、场地地震效应研究等研究的重要基础。面波勘探方法是一种利用人工震源激发产生的瑞雷面波来反演地下介质弹性参数特别是剪切波速的现场原位测试方法,工作效率高且经济成本低。本文针对面波勘探中相对较新的多道面波分析方法(MASW)在地下介质剪切波速探测中的应用进行了研究。主要研究内容及取得的成果如下:(1)利用传递矩阵法(Dunkin)计算了水平层位条件下的瑞雷波相速度频散曲线,并总结     

SASW及MASW方法在隐伏断层探测中的应用

针对地壳浅部介质横向变化大的特点,用瞬态瑞利波SASW与MASW方法对河北夏垫断裂进行了探测与数据处理分析.分别用时间域互相关法和频率域f-k分析法提取频散曲线,然后进行反演得到横波速度结构.实验结果表明瞬态瑞利波法在应用中便捷、灵活,对于探测地表介质中非均匀异常体及断裂向浅部的延伸情况具有独特的应用前景.     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.     

Generation of a pseudo-2D shear-wave velocity section by inversion of a series of 1D dispersion curves

Multichannel Analysis of Surface Waves utilizes a multichannel recording system to estimate near-surface shear (S)-wave velocities from high-frequency Rayleigh waves. A pseudo-2D S-wave velocity (v_S) section is constructed by aligning 1D models at the midpoint of each receiver spread and using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. The receiver spread length sets the theoretical lower limit and any v_S structure with its lateral dimension smaller than this length will not be properly resolved in the final v_S section. A source interval smaller than the spread length will not improve the horizontal resolution because spatial smearing has already been introduced by the receiver spread.In this paper, we first analyze the horizontal resolution of a pair of synthetic traces. Resolution analysis shows that (1) a pair of traces with a smaller receiver spacing achieves higher horizontal resolution of inverted S-wave velocities but results in a larger relative error; (2) the relative error of the phase velocity at a high frequency is smaller than at a low frequency; and (3) a relative error of the inverted S-wave velocity is affected by the signal-to-noise ratio of data. These results provide us with a guideline to balance the trade-off between receiver spacing (horizontal resolution) and accuracy of the inverted S-wave velocity. We then present a scheme to generate a pseudo-2D S-wave velocity section with high horizontal resolution using multichannel records by inverting high-frequency surface-wave dispersion curves calculated through cross-correlation combined with a phase-shift scanning method. This method chooses only a pair of consecutive traces within a shot gather to calculate a dispersion curve. We finally invert surface-wave dispersion curves of synthetic and real-world data. Inversion results of both synthetic and real-world data demonstrate that inverting high-frequency surface-wave dispersion curves – by a pair of traces through cross-correlation with phase-shift scanning method and with the damped least-square method and the singular-value decomposition technique – can feasibly achieve a reliable pseudo-2D S-wave velocity section with relatively high horizontal resolution.     

Research on the middle-of-receiver-spread assumption of the MASW method

The multichannel analysis of surface wave (MASW) method has been effectively used to determine near-surface shear- (S-) wave velocity. Estimating the S-wave velocity profile from Rayleigh-wave measurements is straightforward. A three-step process is required to obtain S-wave velocity profiles: acquisition of a multiple number of multichannel records along a linear survey line by use of the roll-along mode, extraction of dispersion curves of Rayleigh waves, and inversion of dispersion curves for an S-wave velocity profile for each shot gather. A pseudo-2D S-wave velocity section can be generated by aligning 1D S-wave velocity models. In this process, it is very important to understand where the inverted 1D S-wave velocity profile should be located: the midpoint of each spread (a middle-of-receiver-spread assumption) or somewhere between the source and the last receiver. In other words, the extracted dispersion curve is determined by the geophysical structure within the geophone spread or strongly affected by the source geophysical structure. In this paper, dispersion curves of synthetic datasets and a real-world example are calculated by fixing the receiver spread and changing the source location. Results demonstrate that the dispersion curves are mainly determined by structures within a receiver spread.     

Multichannel analysis of surface wave method with the autojuggie

The shear (S)-wave velocity of near-surface materials and its effect on seismic-wave propagation are of fundamental interest in many engineering, environmental, and groundwater studies. The multichannel analysis of surface wave (MASW) method provides a robust, efficient, and accurate tool to observe near-surface S-wave velocity. A recently developed device used to place large numbers of closely spaced geophones simultaneously and automatically (the ‘autojuggie’) is shown here to be applicable to the collection of MASW data. In order to demonstrate the use of the autojuggie in the MASW method, we compared high-frequency surface-wave data acquired from conventionally planted geophones (control line) to data collected in parallel with the automatically planted geophones attached to steel bars (test line). The results demonstrate that the autojuggie can be applied in the MASW method. Implementation of the autojuggie in very shallow MASW surveys could drastically reduce the time required and costs incurred in such surveys.     

Application of surface-wave tomography to mineral exploration: a case study from Siilinjärvi,Finland

In order to assess the feasibility and validity of surface-wave tomography as a tool for mineral exploration, we present an active seismic three-dimensional case study from the Siilinjärvi mine in Eastern Finland. The aim of the survey is to identify the formation carrying the mineralization in an area south of the main pit, which will be mined in the future. Before acquiring the data, we performed an accurate survey design to maximize data coverage and minimize the time for deployment and recollection of the equipment. We extract path-averaged Rayleigh-wave phase-velocity dispersion curves by means of a two-station method. We invert them using a computationally efficient tomographic code which does not require the computation of phase-velocity maps and inverts directly for one-dimensional S-wave velocity models. The retrieved velocities are in good agreement with the data from a borehole in the vicinity, and the pseudo three–dimensional S-wave velocity volume allows us to identify the geological contact between the formation hosting most of the mineralization and the surrounding rock. We conclude that the proposed method is a valid tool, given the small amount of equipment used and the acceptable amount of time required to process the data.     

Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements

Recent field tests illustrate the accuracy and consistency of calculating near-surface shear (S)-wave velocities using multichannel analysis of surface waves (MASW). S-wave velocity profiles (S-wave velocity vs. depth) derived from MASW compared favorably to direct borehole measurements at sites in Kansas, British Columbia, and Wyoming. Effects of changing the total number of recording channels, sampling interval, source offset, and receiver spacing on the inverted S-wave velocity were studied at a test site in Lawrence, Kansas. On the average, the difference between MASW calculated V_s and borehole measured V_s in eight wells along the Fraser River in Vancouver, Canada was less than 15%. One of the eight wells was a blind test well with the calculated overall difference between MASW and borehole measurements less than 9%. No systematic differences were observed in derived V_s values from any of the eight test sites. Surface wave analysis performed on surface data from Wyoming provided S-wave velocities in near-surface materials. Velocity profiles from MASW were confirmed by measurements based on suspension log analysis.     

Apparent phase velocities and fundamental-mode phase velocities of Rayleigh waves

Rayleigh wave dispersion data usually appear in the form of multimodal spectra for a layered model structure. The number of dispersion modal curves depends on the number of layers in the model. The measured dispersion velocities from the frequency–wavenumber (F–K) space, however, may not represent the true phase velocities of the fundamental-mode dispersion curve, but apparent phase velocities. The present study discusses how multimode curves are generated in the multichannel analysis of surface waves (MASW) method and the cause of the apparent velocity produced by the F–K method. Results from a field trial demonstrate that apparent phase velocities fail to reveal thin layers with low velocities. A better agreement of the inverted model with the geotechnical record is obtained by using the data points extracted from the fundamental-mode curve of the MASW spectral image.     

The Selection of Field Acquisition Parameters for Dispersion Images from Multichannel Surface Wave Data

The accuracy and resolution of surface wave dispersion results depend on the parameters used for acquiring data in the field. The optimized field parameters for acquiring multichannel analysis of surface wave (MASW) dispersion images can be determined if preliminary information on the phase velocity range and interface depth is available. In a case study on a fill slope in Hong Kong, the optimal acquisition parameters were first determined from a preliminary seismic survey prior to a MASW survey. Field tests using different sets of receiver distances and array lengths showed that the most consistent and useful dispersion images were obtained from the optimal acquisition parameters predicted. The inverted S-wave velocities from the dispersion curve obtained at the optimal offset distance range also agreed with those obtained by using direct refraction survey.     

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

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