Initial velocity model construction of seismic tomography in near-surface applications

The impact of initial velocity models on final image reconstruction results and how to construct a proper initial velocity model in near-surface tomography studies are investigated on a two-layer synthetic model with gradually increasing velocity with depth. Refraction initial velocity models and linear velocity function models are tested on both synthetic and field data to obtain images close to reality. It is concluded that velocity function type initial models should be preferred in soft alluvial deposits that exist within the investigated depths, whereas refraction initial models should be preferred in the groundwater table or with strong refractors' existence within the investigated depths to obtain optimum subsurface images in refraction–diving wave seismic tomography.     

The perturbation method in elastic wave scattering

Methods of theoretical study in seismic wave scattering are reviewed with the emphasis on the perturbation method. Detailed analysis for weak scattering using Born approximation is given. For elastic random media, the mean square amplitudes of scattered waves are derived using a new approach by working directly in the spectrum domain. The conditions for the scalar wave approximation are obtained. The problem of sensitivity of fore- and backscattering to heterogeneities with different scales and properties (velocity or impedance) is discussed.     

局域双曲线Radon变换压制散射波研究与应用

一般而言,由于地下非均质体的存在所产生的二次波源,由它再生成新的波场,叫散射波场。目前,将散射波作为有效波来成像,已开始在溶洞和裂缝等特殊地质体的识别中得到应用。但对于野外资料采集来说,地表复杂地区,如戈壁、砾石区和山前带,大量存在的散射波却是干扰波,它们的存在会严重影响资料的品质,而其研究与实际应用国内外还很少。因此,通过正演模拟,分析散射波的基本特征,在此基础上研究散射波的去噪方法显得十分必要。本文从地震波运动学时距关系出发,研究了反射波和散射波的几何特征;然后用有限差分正演,模拟了散射波场,用理论模型研究并测试了局域双曲线Radon变换散射波去噪新方法。对于实际炮集资料,分析了F-K滤波方法压制散射噪声的局限,采用局域双曲Radon变换有效地去除了炮集中存在的散射噪声,取得了较好的应用效果。     

二阶Born近似有限频率走时敏感核

有限频率层析成像考虑了非均匀介质中波的散射、衍射、波前愈合等物理性质,使得其对速度异常体的分辨能力远大于射线层析成像.推导和计算有限频率敏感核是进行有限频率层析成像的关键,当前推导有限频率敏感核多借助一阶Born近似,但这只适用于弱散射介质的情况.本文基于二阶Born近似并利用傅里叶变换推导了三维均匀介质情况下有限频率敏感核的解析表达式,并将其推广到非均匀介质中得到了三维非均匀介质中有限频率敏感核.研究表明:当介质中速度扰动小于2%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核差别很小,可近似认为相同;当介质中速度扰动大于5%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核有较大不同,表明此时已不能忽略二次散射.     

地震多级散射波正演模拟方法

地震波在穿越地下散射体群时会产生多级散射波,分析其地震响应特征,可推断散射体的分布情况和性质。本文从二维标量波动方程出发,结合地震散射理论和波恩近似理论,推导了多级散射波方程。在此基础上,采用高阶有限差分法对双点散射体模型和复杂散射体模型进行数值模拟,分析了多级散射波的传播规律和波场特征,并通过抽取多级散射记录和各级散射记录的单道记录与参考单道记录的对比,验证了本文推导散射波方程的准确性。      

瞬变电磁虚拟波场二阶Born近似成像算法

依据瞬变电磁扩散场与虚拟波动场之间的数学表达式,通过波场反变换得到虚拟波场,进而借用地震勘探上成熟的成像方法定位地质界面的位置以及形态.本文首先将时间域波动方程转换为频率域波动方程,将总波动场分为背景场和散射场、总波速分为背景速度和扰动速度,然后运用Green定理得到了散射场关于地质界面速度扰动量的表达式,并使用Bor...     

True amplitude imaging by inverse generalized Radon transform based on Gaussian beam decomposition of the acoustic Green's function

True amplitude migration is one of the most important procedures of seismic data processing. As a rule it is based on the decomposition of the velocity model of the medium into a known macrovelocity component and its sharp local perturbations to be determined. Under this decomposition the wavefield can be considered as the superposition of an incident and reflected/scattered waves. The single scattering approximation introduces the linear integral operator that connects the sharp local perturbations of the macrovelocity model with the multishot/multioffset data formed from reflected/scattered waves. We develop the pseudoinverse of this operator using the Gaussian beam based decomposition of acoustic Green's functions. The computation of this pseudoinverse operator is done pointwise by shooting Gaussian beams from the target area towards the acquisition system. The numerical implementation of the pseudoinverse operator was applied to the synthetic data Sigsbee2A. The results obtained demonstrate the high quality of the true amplitude images computed both in the smooth part of the model and under the salt body.     

Vector-scattering of elastic waves by directional structural space gradients

A nonstochastic and noniterative theory of vector scattering in inhomogeneous media is presented. The elastodynamic vector wave-equation for 3D inhomogeneous media is solved for a weak heterogeneity at the high-frequency region. It is shown that there exists a forward scattered field which decays slowly along the source-receiver path. Its rate of attenuation depends on the azimuth of the path relative to the direction of the inhomogeneity, but is independent of frequency. The Green's tensor for the above regime is derived in closed form and leads to the quantification of fields of dipolar sources in weak inhomogeneous media. The inhomogeneity at the source creates a source-induced scattering (in addition to path-scattering) having a radiation-pattern that bears the signature of the source. The availability of the analytic Green's tensor, in conjunction with the Huygens-Kirchhoff-Helmholtz formalism, opens new ways to calculate the scattered fields due to various structural inhomogeneities applicable to exploration and earthquake seismology. The theoretical results of this study point to the conclusion that the scalar wave approximation may not always be valid for the propagation of seismic waves in the earth's lithosphere.     

Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media

Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.     

地震干涉测量法近地表散射波分离技术

针对山地地震勘探数据低信噪比问题,近地表散射波分离意义显得尤为突出,地震干涉测量法为此提供了一种技术手段.本文将地震干涉测量理论和散射理论结合起来,导出了近地表散射波地震干涉测量表达式,分为互相关型和褶积型表达式,它们由实际波场和背景波场干涉测量构成.根据近地表散射波分离理论,结合陆上地震勘探实际观测系统,采用褶积和反褶积混合型地震干涉测量配置,用实际地震资料展示了近地表散射波分离技术的应用效果.经过理论分析和砾石区实际资料试验,表明地震干涉测量不仅能分离测线上散射源产生的散射波,而且能分离部分侧面散射波.该技术的优点在于它适应于起伏地形和不均匀近地表结构,并且不需要起伏地形和近地表速度信息.为了从实际资料中消除近地表散射波,本文采用多道匹配滤波自适应减法,在砾石区见到较好效果.     

地震波散射理论及应用研究综述

回顾了地震波散射的概念、研究历史及分类.对散射理论及应用研究成果进行了系统评述，特别对前人关于散射波数值模拟和应用研究的不足进行了总结：(1)观测系统设计不合理；(2)对散射波的特征研究缺乏系统性；(3)对复杂散射波场的识别方法没有进行系统研究；(4)适合散射成像的处理软件还有待开发.并介绍了作者关于地震散射波场特征数值模拟研究的初步成果.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.     

Azimuthally pre-stack seismic inversion for orthorhombic anisotropy driven by rock physics

Based on the long-wavelength approximation, a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study, we propose a reliable method for predicting the Thomsen’s weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen’s weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen’s weak anisotropic parameters and fracture weaknesses. Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth (AVAZ) inversion for Thomsen’s weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares (IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen’s weak anisotropic parameters and fracture weaknesses.     

Prestack data enhancement with phase corrections in time-frequency domain guided by local multidimensional stacking

We present a new approach to enhancing weak prestack reflection signals without sacrificing higher frequencies. As a first step, we employ known multidimensional local stacking to obtain an approximate ‘model of the signal’. Guided by phase spectra from this model, we can detect very weak signals and make them visible and coherent by ‘repairing’ corrupted phase of original data. Both presented approaches – phase substitution and phase sign corrections – show good performance on complex synthetic and field data suffering from severe near-surface scattering where conventional processing methods are rendered ineffective. The methods are mathematically formulated as a special case of time-frequency masking (common in speech processing) combined with the signal model from local stacking. This powerful combination opens the avenue for a completely new family of approaches for multi-channel seismic processing that can address seismic processing of land data with nodes and single sensors in the desert environment.     

地震散射波的高精度数值模拟与振幅分析

随着地震勘探要求和技术的不断提高,散射波地震勘探方法逐渐引起人们的重视,尤其在构造复杂地区的高精度勘探中,散射波理论具有广泛的应用前景.本文在讨论了现阶段几种散射波数值模拟方法的利弊后,设计了一种基于微扰论的FK域积分法,在散射场的二次震源和空间能量衰减处理两方面做了改进,其计算精度和效率均得到了提高.此外,采用该算法针对散射波动力学进行了正演模拟试算,从扰动介质空间分布和扰动量两个角度定量地分析了散射波振幅动力学方面的特征.在此基础上得到了五个经验常数,并给出了对应的计算公式和物理意义,通过这五个常数可以划分散射波不同形态的动力学特征.     

时间二阶积分波场的全波形反演

通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异.     

煤层陷落柱散射波数值模拟与成像

煤层陷落柱是煤田勘探开发中常见的一种典型的非均匀地质体.由于来自陷落柱的反射信号少、反射能量弱,使得基于反射波原理的常规地震成像方法难以有效识别陷落柱.本文以散射波理论为基础,采用数值模拟方法,研究了陷落柱的散射波场特征,研究表明地面接收的波场中含有来自陷落柱陡倾角界面的散射波场.通过共散射点道集波场的模拟,可以清晰地识别散射波,获得地下散射点和非均匀地质体的信息,判断散射点的位置,从而勾画出不均匀地质体的形态.采用等效偏移距假设抽取共散射点道集,在此基础上进行叠前偏移,对陷落柱成像;模拟与实际数据成像结果对比表明此方法能够合理地提取散射点的散射波场信息,对陷落柱形态及内部结构准确成像,是一种有效的煤田陷落柱成像方法.     

逆散射理论的发展及应用前景

逆散射理论作为量子力学发展的产物,目前在地震数据处理方面得到了较为广泛的应用.逆散射成像是通过散射体外部场的探测来获得散射体内部的信息,它是近年来发展比较迅速的一种成像方法,其突出的应用是运用逆散射序列中的一阶近似项进行深度成像.本文对逆散射理论在国内外的发展历程进行较为系统的概括和总结,在此基础上,对逆散射理论的应用前景进行了分析.逆散射理论在多次波的压制、界面奇性反演、大扰动深度成像等方面具有广阔的应用前景.     

Model‐based attenuation for scattered dispersive waves

Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data.     

基于卷积神经网络的地震偏移剖面中散射体的定位和成像

本文提出了一种基于深度学习卷积神经网络(CNN)的全波形反演方法,可对地震散射波场中的散射体进行成像和定位.本文的灵感来自如下猜想:在散射波场剖面上的每个点附近的局部波场与该点到各散射体之间的最小距离有关系,并且这个关系可以被CNN网络所识别.我们将该最小距离定义为散射距离场,并将散射距离场的类别(即大小等级)作为CNN网络的预期输出,而输入就是该点附近的局部波场.最后用上述CNN网络对散射波场进行逐点训练和识别.计算结果证实了我们的灵感猜想,即上述CNN网络能够在复杂散射波场中对散射体进行成像.只通过一个训练模型的学习,CNN网络即可反演多种散射模型的偏移剖面,最后得到"类别函数预测值"和"滤波剖面"两种成像结果,由此可以辨识出在复杂的偏移剖面中各散射体的位置.