共查询到20条相似文献，搜索用时 31 毫秒
1.
Effect of lateral heterogeneity on surface wave testing: Numerical simulations and a countermeasure 总被引：1，自引：0，他引：1
The surface wave (S-wave) method has gained popularity in engineering practice for determining S-wave velocity depth profiles. A growing trend is towards the application of S-wave testing for spatially 2-D S-wave velocity tomography, ignoring the assumption of horizontally layered medium. A fourth-order velocity-stress finite difference method is used to perform numerical simulations of S-wave testing in earth models with lateral variation. Results show that the lateral heterogeneity induces a non-stationary property in the space domain, resulting in false depth-related dispersion or higher modes if conventional approach based on stationary assumption is used for the dispersion analysis. Artifacts maybe introduced in spatially 2-D S-wave velocity imaging if the effect of lateral heterogeneity is not accounted for. As a potential countermeasure, a high-lateral-resolution S-wave method is proposed to reduce the effect of lateral heterogeneity while maintaining the resolution and depth range of dispersion analysis. It consists of a walk-away survey and a phase-seaming procedure when synthesizing seismograms with different nearest source-to-receiver offset, allowing wide-wavelength dispersion analysis within a small spatial range. The proof of concept is given with several numerical examples. They show that the high-resolution S-wave method can greatly alleviate the effect of lateral heterogeneity and increase spatial resolution. 相似文献
2.
由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质. 相似文献
3.
在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息. 相似文献
4.
A modified approach to surface wave dispersion analysis using active sources is proposed. The method is based on continuous recordings, and uses the continuous wavelet transform to analyze the phase velocity dispersion of surface waves. This gives the possibility to accurately localize the phase information in time, and to isolate the most significant contribution of the surface waves. To extract the dispersion information, then, a hybrid technique is applied to the narrowband filtered seismic recordings. The technique combines the flexibility of the slant stack method in identifying waves that propagate in space and time, with the resolution of f–k approaches. This is particularly beneficial for higher mode identification in cases of high noise levels. To process the continuous wavelet transform, a new mother wavelet is presented and compared to the classical and widely used Morlet type. The proposed wavelet is obtained from a raised-cosine envelope function (Hanning type). The proposed approach is particularly suitable when using continuous recordings (e.g., from seismological-like equipment) since it does not require any hardware-based source triggering. This can be subsequently done with the proposed method. Estimation of the surface wave phase delay is performed in the frequency domain by means of a covariance matrix averaging procedure over successive wave field excitations. Thus, no record stacking is necessary in the time domain and a large number of consecutive shots can be used. This leads to a certain simplification of the field procedures. To demonstrate the effectiveness of the method, we tested it on synthetics as well on real field data. For the real case we also combine dispersion curves from ambient vibrations and active measurements. 相似文献
5.
高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,同时极大地提高了频散谱的成像分辨率.然而,在研究中我们发现,仅利用两个台站的地震记录,并不能将面波的频散特性完整地考虑在内,导致LSC和NLSC方法对高阶模式的成像存在较大的误差.由于主动源面波勘探多道采集的方式,基于信号比较理论的多道信号比较法(MSC)充分利用多道地震信号,可以获得准确的多模式成像,然而该方法需要计算任意两道的频散谱并叠加,存在冗余的计算,导致计算效率较低.因此,本研究对MSC方法进行了相应的改进,通过追踪炮集记录上的面波波组提高了原方法的计算效率,同时,利用理论频散曲线进行叠加分析,验证了改进的MSC方法的正确性和有效性.通过与相移法、LSC和NLSC方法的对比分析,展示了MSC方法是一种准确的、高分辨率的面波多模式频散谱成像方法.实际地震资料的应用,揭示了MSC方法在浅层地震勘探中用于提取面波多模式频散信息的巨大潜力. 相似文献
6.
W.A. Mulder 《Geophysical Prospecting》2014,62(1):17-33
Migration velocity analysis with the constant‐density acoustic wave equation can be accomplished by the focusing of extended migration images, obtained by introducing a subsurface shift in the imaging condition. A reflector in a wrong velocity model will show up as a curve in the extended image. In the correct model, it should collapse to a point. The usual approach to obtain a focused image involves a cost functional that penalizes energy in the extended image at non‐zero shift. Its minimization by a gradient‐based method should then produce the correct velocity model. Here, asymptotic analysis and numerical examples show that this method may be too sensitive to amplitude peaks at large shifts at the wrong depth and to artefacts. A more robust alternative is proposed that can be interpreted as a generalization of stack power and maximizes the energy at zero‐subsurface shift. A real‐data example is included. 相似文献
7.
Reverse‐time migration has become an industry standard for imaging in complex geological areas. We present an approach for increasing its imaging resolution by employing time‐shift gathers. The method consists of two steps: (i) migrating seismic data with the extended imaging condition to get time‐shift gathers and (ii) accumulating the information from time‐shift gathers after they are transformed to zero‐lag time‐shift by a post‐stack depth migration on a finer grid. The final image is generated on a grid, which is denser than that of the original image, thus improving the resolution of the migrated images. Our method is based on the observation that non‐zero‐lag time‐shift images recorded on the regular computing grid contain the information of zero‐lag time‐shift image on a denser grid, and such information can be continued to zero‐lag time‐shift and refocused at the correct locations on the denser grid. The extra computational cost of the proposed method amounts to the computational cost of zero‐offset migration and is almost negligible compared with the cost of pre‐stack shot‐record reverse‐time migration. Numerical tests on synthetic models demonstrate that the method can effectively improve reverse‐time migration resolution. It can also be regarded as an approach to improve the efficiency of reverse‐time migration by performing wavefield extrapolation on a coarse grid and by generating the final image on the desired fine grid. 相似文献
8.
高斯束逆时偏移结合了射线类偏移的高计算效率和波动方程逆时偏移的高精度,能很好地处理焦散点、大倾角成像问题,并且具有面向目标成像的能力.多分量地震资料的偏移技术可以对地下复杂构造进行更准确的成像,由于实际地下介质具有黏滞性,研究黏弹性叠前逆时偏移具有一定的现实意义.本文采用高斯束逆时偏移方法对多分量地震数据进行吸收衰减补偿,首先分别给出纵波和转换波共炮域高斯束叠前逆时偏移方法原理,在此基础上推导补偿吸收衰减的表达式,校正Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的多分量高斯束叠前逆时偏移.数值模型的测试结果显示,在考虑地下介质的黏滞性时,本文方法具有更高的成像分辨率. 相似文献
9.
介质的黏滞性是普遍存在的.黏滞性介质中的真振幅成像需要校正由介质的黏滞性引起的振幅衰减与速度频散,然而常规的反Q偏移方法存在不稳定问题.本文在反演的框架下求解黏声介质成像问题,在有效避开不稳定的同时实现真振幅成像.首先将波动方程线性化,并依此建立黏声介质最小平方逆时偏移（LSRTM）的目标函数;然后推导波动方程伴随算子,并在此基础上借助伴随状态法推导迭代求解的具体算法;最后通过引入动态相位编码技术将计算量降至与常规逆时偏移相同的数量级.该方法在真振幅成像过程中考虑了介质黏滞性的影响,更接近实际情况,具有更好的振幅保持性.相对于常规逆时偏移,该方法能够自动压制成像噪声,具有更高的成像分辨率及精度.通过模型试算验证了方法的正确性. 相似文献
10.
The spectral analysis of surface waves (SASW) method is an in situ, seismic method for determining the shear wave velocity (or maximum shear modulus) profile of a site. The SASW test consists of three steps: field testing, evaluation of dispersion curve by phase unwrapping method, and determination of shear modulus profile by inversion process. In general, field testing and dispersion curve evaluation are regarded as simple work. However, because of characteristic of Fourier transform used in the conventional phase unwrapping method, dispersion curve is sensitive to background noise and body waves in the low frequency range. Furthermore, under some field conditions such as pavement site, the usual phase unwrapping method can lead to erroneous dispersion curve. To overcome problem of the usual phase unwrapping method, in this paper, a new method of determining dispersion curve for SASW method was applied using time–frequency analysis based on harmonic wavelet transform as an alternative method of a current phase unwrapping method. To estimate the applicability of proposed method to SASW method, numerical simulations at various layered soil and pavement profiles were performed and the dispersion curves by proposed method are more reliable than those by the usual phase unwrapping method. 相似文献
11.
Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models 
下载免费PDF全文
下载免费PDF全文 In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information. 相似文献
12.
时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质. 相似文献
13.
A deep-learning-based approach for seismic surface-wave dispersion inversion (SfNet) with application to the Chinese mainland 下载免费PDF全文
下载免费PDF全文 Surface-wave tomography is an important and widely used method for imaging the crust and upper mantle velocity structure of the Earth. In this study, we proposed a deep learning (DL) method based on convolutional neural network (CNN), named SfNet, to derive the vS model from the Rayleigh wave phase and group velocity dispersion curves. Training a network model usually requires large amount of training datasets, which is labor-intensive and expensive to acquire. Here we relied on synthetics generated automatically from various spline-based vS models instead of directly using the existing vS models of an area to build the training dataset, which enhances the generalization of the DL method. In addition, we used a random sampling strategy of the dispersion periods in the training dataset, which alleviates the problem that the real data used must be sampled strictly according to the periods of training dataset. Tests using synthetic data demonstrate that the proposed method is much faster, and the results for the vS model are more accurate and robust than those of conventional methods. We applied our method to a dataset for the Chinese mainland and obtained a new reference velocity model of the Chinese continent (ChinaVs-DL1.0), which has smaller dispersion misfits than those from the traditional method. The high accuracy and efficiency of our DL approach makes it an important method for vS model inversions from large amounts of surface-wave dispersion data. 相似文献
14.
Rayleigh-wave phase velocities have been utilized to determine shear (S)-wave velocities in near-surface geophysics since early 1980s. One of the key steps is to calculate theoretical dispersion curves of an earth model. When the S-wave velocity of the surface layer is higher than some of the layers below, however, the Rayleigh-wave phase velocity in a high-frequency range calculated by existing algorithms approaches the lowest S-wave velocity among the layers above the half-space, rather than a value related to the S-wave velocity of the surface layer. According to our numerical modeling results based on wave equation, trends of the Rayleigh-wave dispersive energy approach about a 91% of the S-wave velocity of the surface layer at a high-frequency range when its wavelength is much shorter than the thickness of the surface layer, which cannot be fitted by a dispersion curve calculated by existing algorithms. We propose a method to calculate Rayleigh-wave phase velocities of models with a high-velocity surface layer by considering its penetration depth. We build a substituted model that only contains the layer with the lowest S-wave velocity among the layers above the half-space and the layers above it. We use the substituted model to replace the original model to calculate phase velocities when the Rayleigh-wave wavelength is not long enough to penetrate the lowest S-wave velocity layer. Several synthetic models are used to verify fitness between the dispersion curve calculated by our proposed method and the trend of the highest dispersive energy. Examples of inversion also demonstrate high accuracy of using our method as the forward calculation method during the inversions. 相似文献
15.
作为近地表横波速度结构成像的主要手段之一,面波多道分析法的正问题研究对现场观测系统设计及后续反演计算具有重要意义.目前面波频散曲线的正演主要分为两类:一是对水平层状介质中面波的本征值问题进行求解,该类方法计算效率高但较难考虑地下介质在横向上的不均匀性;二是基于波动方程的全波场模拟,该类方法在理论上可考虑任意复杂的地质模型但计算成本相对较高.本文基于振幅归一化加权的聚束分析,提出了一种适用于横向非均匀介质模型的多道瑞雷波频散曲线正演方法.首先,基于聚束分析的计算公式推导得到了经振幅归一化加权后输出功率谱中相速度与局部相速度之间的关系,然后通过黄金分割极值搜索算法计算得到了多道瑞雷波数据的理论频散曲线.数值分析结果表明,该算法能够快速地实现横向非均匀介质中多道瑞雷波频散曲线的正演计算,所求取的频散曲线与采用二维弹性波时间域有限差分模拟分析得到的结果误差较小,这在一定程度上说明了该计算方法的可靠性,从而可为面波多道分析法中的观测系统快速优化设计以及横向非均匀介质中频散曲线的反演解释提供理论支撑. 相似文献
16.
多次波压制方法的研究一直都是地震数据处理中非常重要的一个课题.由于常用的多次波匹配方法主要针对多次波模型和实际多次波存在的振幅或相位的差异进行匹配校正,而无法直接校正多次波模型和实际多次波存在时移误差.本文构建了一种复曲波变换的算法,利用复曲波变换的时移不变性质,通过调整复曲波系数的振幅和相位实现对多次波模型振幅和时移误差的校正.为了更好地保护有效信号,在一次波和多次波分离前,引入一个非线性屏蔽滤波器,可以事先分离出大部分有效波,然后再将剩余部分数据作为输入数据,在复曲波域进行剩余一次波和多次波分离.最后通过模型试算和实际资料处理验证了本文提出的一次波和多次波分离方法的有效性. 相似文献
17.
Xiao-Ping Li 《Geophysical Prospecting》1997,45(6):945-961
A deterministic pure phase shift filter (PPSF) is developed to extract the fundamental mode from multimode surface Love waves. Because of different phase velocities of modes and hence different phase traveltimes for a fixed travel distance between source and receiver, the deterministic PPSF can be computed, provided that the dispersion relation of the medium is estimated from the existing transmission data. The process consists of (a) applying the deterministic PPSF to the multimode wave (this step of the filtering process results in a time series in which the amplitudes of the fundamental mode appear at acausal times and the amplitudes of higher modes appear at causal times); (b) setting amplitude values equal to zero for positive times; (c) applying the inverse PPSF to the filtered signal. By using such a deterministic PPSF process, the higher modes almost disappear. The method is applied to synthetic multimode data computed by the normal-mode summation method. 相似文献
18.
波场延拓是地震偏移成像的基础. 快速进行目标区波场延拓对石油勘探中急需发展的深部地震勘探和无组合海量地震数据的成像有重要意义. 在目标区成像中,目前已有的波场延拓方法,包括基于走时计算的Dix方法和射线追踪方法,以及基于小步长波场递推的方法,在适应复杂介质、计算精度和计算效率的某一方面还不能完全满足实际需要. 本文提出一种基于“算子相位”李代数积分的快速计算延拓算子的方法,称为大步长波场延拓方法. 在该方法中,指向目标区的波场延拓算子象征的复相位被表示成波数的线性组合. 线性组合的系数是层速度函数及其导数的深度积分,计算和存储较为方便. 波场延拓算子通过相移算子加校正的方法,利用快速Fourier变换在空间域和波数域予以实现. 利用动力学等价关系导出了便于计算的表达式. 本文比较了算子主象征函数用一步法展开和用两步法展开的精度,从而说明大步长方法的精度要高于递推方法. 在横向和纵向线性变化介质中,将大步长方法的脉冲响应与递推法做了比较,说明大步长延拓算子的走时精度主要取决于相移因子中的横向变速校正项;且在各种近似下,大步长算子发生的频散都非常小. 相似文献
19.
Surface wave methods are based on the inversion of observed Rayleigh wave phase‐velocity dispersion curves. The goal is to estimate mainly the shear‐wave velocity profile of the investigated site. The model used for the interpretation is 1D, hence results obtained wherever lateral variations are present cannot be considered reliable. In this paper, we study four synthetic models, all with a lateral heterogeneity. When we process the entire corresponding seismograms with traditional f‐k approach, the resulting 1D profiles are representative of the subsurface properties averaged over the whole length of the receivers lines. These results show that classical analysis disregards evidences of sharp lateral velocity changes even when they show up in the raw seismograms. In our research, we implement and test over the same synthetic models, a novel robust automated method to check the appropriateness of 1D model assumption and locate the discontinuities. This new approach is a development of the recent multi‐offset phase analysis with the following further advantages: it does not need previous noise evaluation and more than one shot. Only once the discontinuities are clearly identified, we confidently perform classical f‐k dispersion curve extraction and inversion separately on both sides of the discontinuity. Thus the final results, obtained by putting side by side the 1D profiles, are correct 2D reconstructions of the discontinuous S‐wave distributions obtained without any additional ad‐hoc hypotheses. 相似文献
20.
Juan Diego Jaramillo Juan David Gómez Doriam Restrepo Santiago Rivera 《地震工程与工程振动(英文版)》2014,13(3):535-543
A method is proposed to estimate the seismic soil-structure-interaction (SSI) effects for use in engineering practice. It is applicable to 2D structures subjected to vertically incident shear waves supported by homogenous half-spaces. The method is attractive since it keeps the simplicity of the spectral approach, overcomes some of the difficulties and inaccuracies of existing classical techniques and yet it considers a physically consistent excitation. This level of simplicity is achieved through a response spectra modification factor that can be applied to the free-field 5%-damped response spectra to yield design spectral ordinates that take into account the scattered motions introduced by the interaction effects. The modification factor is representative of the Transfer Function (TF) between the structural relative displacements and the free-field motion, which is described in terms of its maximum amplitude and associated frequency. Expressions to compute the modification factor by practicing engineers are proposed based upon a parametric study using 576 cases representative of actual structures. The method is tested in 10 cases spanning a wide range of common fundamental vibration periods. 相似文献