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

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

波动方程偏移成像阴影的照明补偿

受地下复杂构造和地震数据采集系统的影响,使地震波对地下目标的照明出现不均匀性,地震采集系统难以有效地获取地下某些目标的反射信息,进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开,利用角度域波动方程偏移成像和波动方程照明分析,并结合波动方程反演理论,提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响,消除复杂构造区的偏移成像阴影,改进波动方程叠前深度偏移成像在复杂构造区的成像效果.     

Extending illumination using all multiples: application to 3D acquisition geometry analysis

Recent advances in survey design have led to conventional common‐midpoint‐based analysis being replaced by subsurface‐based seismic acquisition analysis, with emphasis on advanced techniques of illumination analysis. Among them is the so‐called focal beam method, which is a wave‐equation‐based seismic illumination analysis method. The objective of the focal beam method is to provide a quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on the resolution and angle‐dependent amplitude fidelity of the image. The method distinguishes between illumination and sensing capability of a particular acquisition geometry by computing the focal source beam and the focal detector beam, respectively. Sensing is related to the detection properties of a detector configuration, whereas illumination is related to the emission properties of a source configuration. The focal source beam analyses the incident wavefield at a specific subsurface grid point from all available sources, whereas the focal detector beam analyses the sensing wavefield reaching at the detector locations from the same subsurface grid point. In the past, this method could only address illumination by primary reflections. In this paper, we will extend the concept of the focal beam method to incorporate the illumination due to the surface and internal multiples. This in fact complies with the trend of including multiples in the imaging process. Multiple reflections can illuminate a target location from other angles compared with primary reflections, resulting in a higher resolution and an improved illumination. We demonstrate how an acquisition‐related footprint can be corrected using both the surface and the internal multiples.     

基于平面波照明的偏移成像补偿

受地下复杂构造和地震数据采集系统的影响,地震波对地下目标的照明出现不均匀,在地震数据的偏移成像中出现成像阴影.根据地震数据最小二乘偏移/反演理论,和把地震波场照明结果作为最小二乘偏移/反演中的Hessian矩阵的近似对偏移成像进行补偿的原理,提出一种应用平面波照明结果对平面波偏移成像结果进行补偿以消除偏移成像阴影的方法.这种基于平面波照明的偏移成像补偿方法相对于局部角度域的照明偏移成像补偿方法具有计算效率上的优势.     

Stacking angle‐domain common‐image gathers for normalization of illumination

Unequal illumination of the subsurface highly impacts the quality of seismic imaging. Different image points receive different folds of reflection‐angle illumination, which can be caused by irregular acquisition or by wave propagation in complex media. Illumination problems can deteriorate amplitudes in migrated images. To address this problem, we present a method of stacking angle‐domain common‐image gathers, in which we use local similarity with soft thresholding to determine the folds of local illumination. Normalization by local similarity regularizes local illumination of reflection angles for each image point of the subsurface model. This approach compensates for irregular illumination by selective stacking in the image space, regardless of the cause of acquisition or propagation irregularities. Additional migration is not required because the methodology is implemented in the reflection angle domain after migration. We use two synthetic examples to demonstrate that our method can normalize migration amplitudes and effectively suppress migration artefacts.     

Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.     

利用井间地震层析成像方法探测隐伏断层

针对目前隐伏断层探测中存在的技术难题,提出了在地面反射地震勘探初步确定的断点附近,用井间地震层析成像方法探测隐伏断层及其上断点的埋深。首先介绍井间地震技术的应用现状和井间地震层析成像的基本原理,并结合实例讨论了井间初至走时层析成像的野外数据采集和室内资料处理方法。芦花台隐伏断层的井间地震层析成像结果表明,采用弯曲射线追踪方法反演得到的跨芦花台隐伏断层的井间纵波速度结构图像与钻孔联合剖面揭示的断层特征相吻合。      

复杂地区观测系统分析评价的照明度函数设计

在复杂构造地区进行照明度分析,对观测系统的优化设计有着重要意义.人们对如何进行正演照明做了大量的研究工作,但是如何利用照明结果进行观测系统的优化设计并没有一个统一的标准.为了使照明度分析有一定的规范,通过文字和图形定义了炮点入射照明、检波点接收照明、共反射点面元入射照明、共反射点面元成像照明和炮点-共反射点面元-检波点...     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.     

Seismic imaging of deep hydrocarbon reservoirs

This paper examines those aspects of reflection seismology which require special consideration when imaging deeper hydrocarbon reservoirs, including the constraints imposed by vertical resolution, lateral resolution, and velocity analysis. We derive quantitative expressions relating the uncertainties in stacking velocities and in interval velocities derived from stacking velocities to acquisition parameters, as well as expressions for the lateral resolution which can theoretically be achieved for migrated seismic images. This analysis shows that the most significant limitations of seismic imaging at depth involve the finite lateral resolution of the seismic method, and the proper lateral positioning of seismic images. These difficulties are overcome in large measure through the proper migration of a seismic dataset, which becomes more critical as deeper horizons are imaged. If these horizons are suspected of having significant 3-D structure, a strong argument may be made for acquiring a 3-D seismic survey over the prospect. Migration of this dataset will then generate an image of the subsurface with good lateral resolution in both the X and Y directions.     

Increasing illumination and sensitivity of reverse‐time migration with internal multiples

Reverse‐time migration is a two‐way time‐domain finite‐frequency technique that accurately handles the propagation of complex scattered waves and produces a band‐limited representation of the subsurface structure that is conventionally assumed to be linear in the contrasts in model parameters. Because of this underlying linear single‐scattering assumption, most implementations of this method do not satisfy the energy conservation principle and do not optimally use illumination and model sensitivity of multiply scattered waves. Migrating multiply scattered waves requires preserving the non‐linear relation between the image and perturbation of model parameters. I modify the extrapolation of source and receiver wavefields to more accurately handle multiply scattered waves. I extend the concept of the imaging condition in order to map into the subsurface structurally coherent seismic events that correspond to the interaction of both singly and multiply scattered waves. This results in an imaging process referred to here as non‐linear reverse‐time migration. It includes a strategy that analyses separated contributions of singly and multiply scattered waves to a final non‐linear image. The goal is to provide a tool suitable for seismic interpretation and potentially migration velocity analysis that benefits from increased illumination and sensitivity from multiply scattered seismic waves. It is noteworthy that this method can migrate internal multiples, a clear advantage for imaging challenging complex subsurface features, e.g., in salt and basalt environments. The results of synthetic seismic imaging experiments, including a subsalt imaging example, illustrate the technique.     

基于目标区的局部角度域方向照明分析及其在偏移成像中的应用

叠前深度偏移已成为近年来地震偏移成像领域广泛应用的技术之一,其中复杂盐丘下方目标区域结构的成像始终是偏移成像中较难解决的问题.局部角度域波场分解可以提供在空间和方向上的双重局部化信息,因而这一技术被广泛地应用在方向照明分析、成像振幅校正等方面.本文在采用局部指数标架小波束进行角度域方向照明分析的基础上,研究采集系统对复杂盐丘下部层状结构成像质量的影响;同时通过分析目标区对应于相对地表采集系统分布的照明能量,确定叠前炮集数据中对盐下目标区成像起重要贡献的部分,并将这些数据用来对目标区成像,从而达到提高盐下结构成像质量的目的.本文将二维SEG-EAGE盐丘模型中盐下反射体作为目标区,分别计算其照明分析和由照明能量分布确定的炮集数据对目标区的偏移成像结果,同时通过与成像振幅校正前后的成像结果对比,说明了该方法对提高目标区结构成像质量的有效性.     

Source side seismic tomography(3STomo):A novelmethod to image the subsurface structure beneath seismically active region

We propose a novel seismic tomography method,Source Side Seismic Tomography(3STomo),which is designed particularly to image the subsurface structure beneath seismically active regions.Unlike the teleseismic tomography,in which the data are relative traveltime residuals between closely spaced stations for each teleseismic event,3STomo uses relative traveltime shifts between earthquakes within the study region for each distant station.Given the relatively evener distribution of global seismic stations,this method has unique advantages for imaging the structure beneath regions that have numerous earthquakes but lack of dense seismic stations,for example,some subduction zones and spreading ridges in the ocean.In addition,3STomo has potentially better vertical resolution at shallow depths than the traditional teleseismic tomography.The effect of the inaccurate source parameters on its resolution can be minimized by using depth phases and the technique of joint source and structure inversion.Numerical experiments and application to Luzon Island,Philippines show that 3STomo can be a valuable tool to investigate the subsurface structure beneath some areas where the traditional method cannot be applied to,or at least it can be used as a complementary component of conventional teleseismic tomography to obtain better back-azimuth coverage and achieve higher resolution at shallow depths in the inversion.     

地震照明叠前深度偏移方法综述

对地震照明叠前深度偏移的基本概念、实现方法进行了分类和阐述.地震照明叠前深度偏移是通过对震源和炮记录进行合理的选择和合成,从而进行地震照明成像的一种有效方法.其可以分为平面波偏移和小束波偏移.理论模型的处理效果证明地震照明叠前深度偏移成像技术有很高的计算效率,并且还可以提高地下特定目标的成像质量.     

Migration of multiples from the South China Sea

The South China Sea where water depth is up to 5000 m is the most promising oil and gas exploration area in China in the future.The seismic data acquired in the South China Sea contain various types of multiples that need to be removed before imaging can be developed.However,compared with the conventional reflection migration,multiples carry more information of the underground structure that helps provide better subsurface imaging.This paper presents a method to modify the conventional reverse time migration so that multiple reflections can migrate to their correct locations in the subsurface.This approach replaces the numerical impulsive source with the recorded data including primaries and multiples on the surface,and replaces the recorded primary reflection data with multiples.In the reverse time migration process,multiples recorded on the surface are extrapolated backward in time to each depth level,while primaries and multiples recorded on the surface are extrapolated forward in time to the same depth levels.By matching the difference between the primary and multiple images using an objective function,this algorithm improves the primary resultant image.Synthetic tests on Sigsbee2 B show that the proposed method can obtain a greater range and better underground illumination.Images of deep water in the South China Sea are obtained using multiples and their matching with primaries.They demonstrate that multiples can make up for the reflection illumination and the migration of multiples is an important research direction in the future.     

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

基于吸收衰减补偿的多分量高斯束逆时偏移

高斯束逆时偏移结合了射线类偏移的高计算效率和波动方程逆时偏移的高精度,能很好地处理焦散点、大倾角成像问题,并且具有面向目标成像的能力.多分量地震资料的偏移技术可以对地下复杂构造进行更准确的成像,由于实际地下介质具有黏滞性,研究黏弹性叠前逆时偏移具有一定的现实意义.本文采用高斯束逆时偏移方法对多分量地震数据进行吸收衰减补偿,首先分别给出纵波和转换波共炮域高斯束叠前逆时偏移方法原理,在此基础上推导补偿吸收衰减的表达式,校正Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的多分量高斯束叠前逆时偏移.数值模型的测试结果显示,在考虑地下介质的黏滞性时,本文方法具有更高的成像分辨率.     

Converted‐wave beam migration with sparse sources or receivers

Gaussian beam depth migration overcomes the single‐wavefront limitation of most implementations of Kirchhoff migration and provides a cost‐effective alternative to full‐wavefield imaging methods such as reverse‐time migration. Common‐offset beam migration was originally derived to exploit symmetries available in marine towed‐streamer acquisition. However, sparse acquisition geometries, such as cross‐spread and ocean bottom, do not easily accommodate requirements for common‐offset, common‐azimuth (or common‐offset‐vector) migration. Seismic data interpolation or regularization can be used to mitigate this problem by forming well‐populated common‐offset‐vector volumes. This procedure is computationally intensive and can, in the case of converted‐wave imaging with sparse receivers, compromise the final image resolution. As an alternative, we introduce a common‐shot (or common‐receiver) beam migration implementation, which allows migration of datasets rich in azimuth, without any regularization pre‐processing required. Using analytic, synthetic, and field data examples, we demonstrate that converted‐wave imaging of ocean‐bottom‐node data benefits from this formulation, particularly in the shallow subsurface where regularization for common‐offset‐vector migration is both necessary and difficult.     

Super-resolution least-squares prestack Kirchhoff depth migration using the <Emphasis Type="Italic">L</Emphasis><Subscript>0</Subscript>-norm

Least-squares migration (LSM) is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L₂-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L₀-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.     

地震照明分析及其在地震采集设计中的应用

地震照明分析是在给定地下背景结构的情况下用模拟方法研究采集系统对地下结构的探测能力.对目标照明的影响主要来自观测系统的设置,上覆地质结构的复杂性以及目标倾角等三项因素.我们将上述因素的影响统一纳入到对照明的定量计算中.本文阐述了关于地震照明的一系列基本概念,并基于波动理论给出对不同类型照明的计算方法.首先将入射和散射波分解到局部角度域并导出局部照明矩阵,它包括了在目标附近所有可能的入射和散射方向对照明的贡献.从照明矩阵出发,研究入射和反射波在角度域中与反射面的相互作用,并由此获得对地下结构的不同照明描述.作为例子,具体给出了对局部照明矩阵,空间体照明,反射面定向照明,反射面照明角度覆盖,以及对反射面可视性的计算方法,并给出了相应的数值计算结果. 最后,讨论了照明分析方法在地震采集设计中的一些可能应用.