Migration for velocity and attenuation perturbations

Migration maps seismic data to reflectors in the Earth. Reflections are not only caused by small‐scale variations of the velocity and density but also of the quality factor that describes attenuation. We investigated scattering due to velocity and attenuation perturbations by computing the resolution function or point‐spread function in a homogeneous background model. The resolution function is the migration image of seismic reflection data generated by a point scatterer. We found that the resolution function mixes velocity and attenuation parameter perturbations to the extent that they cannot be reconstructed independently. This is true for a typical seismic setting with sources and receivers at the surface and a buried scatterer. As a result, it will be impossible to simultaneously invert for velocity and attenuation perturbations in the scattering approach, also known as the Born approximation. We proceeded to investigate other acquisition geometries that may resolve the ambiguity between velocity and attenuation perturbations. With sources and receivers on a circle around the scatterer, in 2D, the ambiguity disappears. It still shows up in a cross‐well setting, although the mixing of velocity and attenuation parameters is less severe than in the surface‐to‐surface case. We also consider illumination of the target by diving waves in a background model that has velocity increasing linearly with depth. The improvement in illumination is, however, still insufficient to remove the ambiguity.     

上月壳中的散射引起月震尾波的数值模拟研究

在阿波罗月震记录中普遍存在着强烈持久的尾波信号,这样的波形特征无法用均匀分层月球模型解释.一个普遍被接受的解释是月震尾波由月球浅层结构对月震波的散射引起.我们采用基于交错网格的伪谱和有限差分混合方法模拟研究非均匀上月壳对月震波的散射效应,在此基础上解释月震尾波的形成机制,并估计出上月壳速度扰动的强度.我们发现,在均匀分层模型基础上,进一步考虑上月壳中的非均匀结构对月震波的散射效应,能有效地解释月震信号中强烈持久的尾波.我们认为月震尾波可能是由上月壳中的低波速、低衰减和散射这三个因素的共同作用所引起.采用不同的扰动标准差模拟上月壳的非均匀性,并比较模拟波形与真实月震图的相似程度,我们发现上月壳中速度扰动的标准差应该在3%到5%之间,很可能接近于3%.     

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

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

接收函数的克希霍夫２D偏移方法

本文将地震反射资料成像的Kirchhoff偏移方法引入到接收函数成像研究，提出了一种接收函数的Kirchhoff偏移方法，以适应介质速度的横向变化，提高接收函数成像的精度和分辨率. 模型检验表明，与传统的水平分层的共转换点偏移叠加方法相比，Kirchhoff偏移方法能够实现转换波的有效归位，消除水平分层共转换点偏移叠加造成的假象. 一个很重要的现象是，倾斜或弯曲界面的成像明显受控于远震的入射方向，其成像精度的提高有赖于不同方向接收函数的联合成像.     

Generalized acoustic diffraction tomogra phy1

A generalized diffraction tomography algorithm is developed, which in principle can handle irregularly spaced data, curved acquisition lines and non-uniform background models. By direct comparison with medical diffraction tomography, it is shown that the generalized method involves the same two processing steps: data filtering and back-propagation. The filter handles the irregular sampling of the model space and the uneven energy coverage, while the back-propagation operator removes the wave propagation effects. Paraxial ray-tracing techniques are employed to compute both these quantities. In medical diffraction tomography, the resolution vector (i.e. the Fourier vector of the model space) is defined by the incident and scattered plane-wave directions. It is shown here that a similar relationship exists for a non-uniform background, where the resolution vector at a particular image point is defined by the incident and scattered ray directions. Consequently, the impulse response of the generalized algorithm becomes space variant. Finally, a general processing procedure for transmission mode seismic data, based on this generalized algorithm, is proposed. The potential of the method is demonstrated using synthetic cross-hole data.     

地层衰减对地震波速度逆散射反演的影响研究

在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想.     

基于逆散射成像条件的最小二乘逆时偏移

最小二乘逆时偏移（LSRTM）相对于常规逆时偏移（RTM）具有分辨率更高、振幅更准确、噪音更少等优势,可以对复杂的地质构造进行有效的成像.这种迭代更新反演成像方法十分依赖目标函数的梯度质量和计算效率.当地质模型中存在强反射界面或者记录中存在折射波时,基于常规互相关成像条件（CCC）的最小二乘逆时偏移梯度会包含很强的低频噪音,从而使反演的收敛速度和成像质量降低.为此,本文在最小二乘逆时偏移的梯度中引进了逆散射成像条件来压制这种低频噪音,并以此提出基于逆散射成像条件（ISC）的最小二乘逆时偏移方法.数值模拟结果表明,两者计算耗时基本一致,但逆散射成像条件能高效压制梯度中的低频噪音,从而使反演过程中收敛加速,成像质量得到显著提高.     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

Image‐domain wavefield tomography with extended common‐image‐point gathers

Waveform inversion is a velocity‐model‐building technique based on full waveforms as the input and seismic wavefields as the information carrier. Conventional waveform inversion is implemented in the data domain. However, similar techniques referred to as image‐domain wavefield tomography can be formulated in the image domain and use a seismic image as the input and seismic wavefields as the information carrier. The objective function for the image‐domain approach is designed to optimize the coherency of reflections in extended common‐image gathers. The function applies a penalty operator to the gathers, thus highlighting image inaccuracies arising from the velocity model error. Minimizing the objective function optimizes the model and improves the image quality. The gradient of the objective function is computed using the adjoint state method in a way similar to that in the analogous data‐domain implementation. We propose an image‐domain velocity‐model building method using extended common‐image‐point space‐ and time‐lag gathers constructed sparsely at reflections in the image. The gathers are effective in reconstructing the velocity model in complex geologic environments and can be used as an economical replacement for conventional common‐image gathers in wave‐equation tomography. A test on the Marmousi model illustrates successful updating of the velocity model using common‐image‐point gathers and resulting improved image quality.     

Wavefield tomography based on local image correlations

The estimation of a velocity model from seismic data is a crucial step for obtaining a high‐quality image of the subsurface. Velocity estimation is usually formulated as an optimization problem where an objective function measures the mismatch between synthetic and recorded wavefields and its gradient is used to update the model. The objective function can be defined in the data‐space (as in full‐waveform inversion) or in the image space (as in migration velocity analysis). In general, the latter leads to smooth objective functions, which are monomodal in a wider basin about the global minimum compared to the objective functions defined in the data‐space. Nonetheless, migration velocity analysis requires construction of common‐image gathers at fixed spatial locations and subsampling of the image in order to assess the consistency between the trial velocity model and the observed data. We present an objective function that extracts the velocity error information directly in the image domain without analysing the information in common‐image gathers. In order to include the full complexity of the wavefield in the velocity estimation algorithm, we consider a two‐way (as opposed to one‐way) wave operator, we do not linearize the imaging operator with respect to the model parameters (as in linearized wave‐equation migration velocity analysis) and compute the gradient of the objective function using the adjoint‐state method. We illustrate our methodology with a few synthetic examples and test it on a real 2D marine streamer data set.     

扩展成像条件下的最小二乘逆时偏移

逆时偏移(RTM)是复杂介质条件下地震成像的重要手段.因受观测系统限制、上覆地层影响以及波场带宽有限等因素的影响,现行的常规RTM所采用的互相关成像条件通常对地下构造进行模糊成像.最小二乘逆时偏移(LSRTM)通过最小化线性Born近似正演数据和采集数据之间的波形差异,采用梯度类反演算法优化反射系数模型,获得的成像结果具有更高的分辨率和更可靠的振幅保真度.然而,基于波形拟合的LSRTM对背景速度模型的依赖性很强.误差太大的速度模型容易产生周波跳跃现象,导致LSRTM难以获得全局最优解.为了克服这一问题,本文基于扩展模型的思想,在线性Born近似下,推导得到RTM扩展成像条件.并基于最小二乘反演理论,提出扩展成像条件下的LSRTM方法.理论模型试算表明,本文方法不仅可以提供分辨率更高、振幅属性更为可靠的成像结果,而且能够在一定程度上消除速度误差对反演成像的影响.     

基于成像射线的偏移剖面和速度时深转换评述

深度域地震波速度模型是地震勘探中最重要的参数之一,是获得精确深度域偏移剖面的基础.本文基于成像射线的概念,阐述了时间偏移域和深度域之间的关系,明确了时间偏移后的像点是通过成像射线与深度域的散射点联系起来的.在有横向变速的情况下,时间偏移剖面需要经过时深转换才能真实反映深度域层位,同样,从时间偏移剖面上提取的速度也需要经过时深转换才能得到深度域的速度.与前者不同的是,在速度的时深转换中,不仅速度的位置需要改变,且其数值也需要进行校正.该校正因子有几种等价的描述形式：度规量、速度扩散因子和几何扩散因子.     

地震勘探中的时间-暴光声学

本文提出用于地下成像无源地震方法的分析,在该方法中应用环境地震噪声作为地下散射体的照明源。该成像算法能够以递归方式将新的数据融入到成像中,这使成像背景噪声随时间减少。在空间域不相干环境背景噪声的假设前提下推导出成像算法的点-扩展函数。点-扩展函数表征成像的分辨率,即接收排列长度和环境带宽的函数。     

跨孔地震层析成像的级联方法

对厚度小于1/4波长的超薄波速干扰体进行高分辨率成像,作者曾采用走时反演和波场反演相结合的方法,体现了由低分辨率向高分辨率逐步逼近的思想.级联算法是这种思想的进一步发展,我们将具有不同分辨率的算法串联起来,以达到高分辨率成像的目的.本文介绍一种三级串联的算法,并进行了算法分析,数值计算的例子说明这种级联算法分辨率高、稳定性好,只需要地震资料而不要求其它先验信息,因此能较好地满足实际应用的要求.     

High‐resolution and super stacking of time‐reversal mirrors in locating seismic sources

Time reversal mirrors can be used to backpropagate and refocus incident wavefields to their actual source location, with the subsequent benefits of imaging with high‐resolution and super‐stacking properties. These benefits of time reversal mirrors have been previously verified with computer simulations and laboratory experiments but not with exploration‐scale seismic data. We now demonstrate the high‐resolution and the super‐stacking properties in locating seismic sources with field seismic data that include multiple scattering. Tests on both synthetic data and field data show that a time reversal mirror has the potential to exceed the Rayleigh resolution limit by factors of 4 or more. Results also show that a time reversal mirror has a significant resilience to strong Gaussian noise and that accurate imaging of source locations from passive seismic data can be accomplished with traces having signal‐to‐noise ratios as low as 0.001. Synthetic tests also demonstrate that time reversal mirrors can sometimes enhance the signal by a factor proportional to the square root of the product of the number of traces, denoted as N and the number of events in the traces. This enhancement property is denoted as super‐stacking and greatly exceeds the classical signal‐to‐noise enhancement factor of . High‐resolution and super‐stacking are properties also enjoyed by seismic interferometry and reverse‐time migration with the exact velocity model.     

Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Wide-angle linear forward modelling of synthetic seismograms

We address the issue of linearity and scale dependence in forward modelling of seismic data from well logs, for large ray parameters, wide angles or large offsets. We present a forward model, within the context of seismic‐to‐well matching, that is linearized in the elastic properties of the earth. This model preserves linearity at large ray parameters and can handle fine‐layering effects such as induced anisotropy. Starting from a low‐contrast small‐ray‐parameter model, we extend it to a large‐ray‐parameter model by fully linearizing the elastic‐property contrasts. Overall linearity of the forward model is extended by partitioning the compressional‐wave and shear‐wave velocity fields into two fundamental scales: a kinematic scale that governs wavefield propagation effects and a dynamic scale that governs wavefield scattering effects. This analysis reveals that the standard practice in forward modelling of strongly filtering the ratio of compressional‐wave velocity to shear‐wave velocity is well founded in the underlying physics. The partitioning of the velocity fields also leads naturally to forward modelling that accounts fully for stretch effects, to resolution of the angle‐of‐incidence versus ray‐parameter dichotomy in seismic‐amplitude analysis, and to full accounting for induced anisotropy and dispersion effects due to fine‐layering of isotropic media. With the onset of routine long‐offset acquisition and the compelling need to optimize asset management in order to maximize reserve recovery, this forward model recognizes the physics of seismic wave propagation and enables a more complete exploitation of amplitude information in pre‐critical seismic data.     

高斯波包反射走时速度反演方法

扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.     

基于Duffing振子混沌系统的地震速度分析方法

强随机噪声干扰是导致地震勘探资料低信噪比的主要原因,如何在强随机噪声干扰下获取有效的信息是值得关注的问题.Duffing振子混沌系统是一个非线性的动力学系统,其对强随机噪声具有免疫能力,而对特定的周期性信号具有敏感性.本文提出一种基于Duffing振子混沌系统的速度分析方法.对CMP道集按照时距曲线关系进行移动窗口截取,将所截取的信号构建为待测信号加入Duffing振子混沌系统,通过相图网格分割方法(GPM)判断系统状态的改变,从而在强随机噪声背景下获得高分辨率的速度谱.理论模型和实际资料的处理结果表明,与传统的水平叠加速度分析方法相比,本方法能够在强随机噪声背景下获得更准确的速度分析结果.     

Extended imaging conditions for wave‐equation migration

Wavefield‐based migration velocity analysis using the semblance principle requires computation of images in an extended space in which we can evaluate the imaging consistency as a function of overlapping experiments. Usual industry practice is to assemble those seismic images in common‐image gathers that represent reflectivity as a function of depth and extensions, e.g., reflection angles. We introduce extended common‐image point (CIP) gathers constructed only as a function of the space‐ and time‐lag extensions at sparse and irregularly distributed points in the image. Semblance analysis using CIP's constructed by this procedure is advantageous because we do not need to compute gathers at regular surface locations and we do not need to compute extensions at all depth levels. The CIP's also give us the flexibility to distribute them in the image at irregular locations aligned with the geologic structure. Furthermore, the CIP's remove the depth bias of common‐image gathers constructed as a function of the depth axis. An interpretation of the CIP's using the scattering theory shows that they are scattered wavefields associated with sources and receivers inside the subsurface. Thus, when the surface wavefields are correctly reconstructed, the extended CIP's are characterized by focused energy at the origin of the space‐ and time‐lag axes. Otherwise, the energy defocuses from the origin of the lag axes proportionally with the cumulative velocity error in the overburden. This information can be used for wavefield‐based tomographic updates of the velocity model, and if the velocity used for imaging is correct, the coordinate‐independent CIP's can be a decomposed as a function of the angles of incidence.