Estimation and application of near‐surface full waveform redatuming operators

Land seismic data quality can be severely affected by near‐surface anomalies. The imprint of a complex near‐surface can be removed by redatuming the data to a level below the surface, from where the subsurface structures are assumed to be relatively smooth. However, to derive a velocity‐depth model that explains the propagation effects of the near‐surface is a non‐trivial task. Therefore, an alternative approach has been proposed, where the redatuming operators are obtained in a data‐driven manner from the reflection event related to the datum. In the current implementation, the estimation of these redatuming operators is done in terms of traveltimes only, based on a high‐frequency approximation. The accompanying amplitudes are usually derived from a local homogeneous medium, which is obviously a simplification of reality. Such parametrization has produced encouraging results in the past but cannot completely remove the near‐surface complexities, leaving artefacts in the redatumed results. In this paper we propose a method that estimates the redatuming operators directly from the data, i.e., without using a velocity model, in a full waveform manner, such that detailed amplitude and phase variations are included. The method directly outputs the inverse propagation operators that are needed for true‐amplitude redatuming. Based on 2D synthetic data it is demonstrated that the resulting redatuming quality is improved and artefacts are reduced.     

利用CFP技术解决复杂近地表问题综述

复杂近地表条件会降低地震数据的质量,通常采用基于地表一致性的时移静校正消除其影响.但静校正与速度是密不可分的,而确定复杂近地表速度是非常困难的.基于CFP技术处理复杂近地表问题时避免了对速度的直接操作,使得静校正和速度的确定相互独立.首先根据叠前数据估算出波场的传播算子,然后依据等时原理在DTS模板中进行算子更新,再用这些更新的算子重建基准面和实现近地表单程时间成像.获得正确的算子振幅也是重建基准面的关键.     

3D sparse-data Kirchhoff redatuming

Wave‐equation redatuming can be a very efficient method of overcoming the overburden imprint on the target area. Owing to the growing amount of 3D data, it is increasingly important to develop a feasible method for the redatuming of 3D prestack data. Common 3D acquisition designs produce relatively sparse data sets, which cannot be redatumed successfully by applying conventional wave‐equation redatuming. We propose a redatuming approach that can be used to perform wave‐equation redatuming of sparse 3D data. In this new approach, additional information about the medium velocity below the new datum is included, i.e. redatumed root‐mean‐square (RMS) velocities, which can be extracted from the input data set by conventional velocity analysis, are used. Inclusion of this additional information has the following implications: (i) it becomes possible to simplify the 4D redatuming integral into a 2D integral such that the number of traces needed to calculate one output time sample and the computational effort are both reduced; (ii) the information about the subsurface enables an infill of traces which are needed for the integral calculation but which are missing in the sparse input data set. Two tests applying this new approach to fully sampled 2D data show satisfactory results, implying that this method can certainly be used for the redatuming of sparse 3D data sets.     

Preserved‐amplitude angle domain migration by shot‐receiver wavefield continuation

We present preserved‐amplitude downward continuation migration formulas in the aperture angle domain. Our approach is based on shot‐receiver wavefield continuation. Since source and receiver points are close to the image point, a local homogeneous reference velocity can be approximated after redatuming. We analyse this approach in the framework of linearized inversion of Kirchhoff and Born approximations. From our analysis, preserved‐amplitude Kirchhoff and Born inverse formulas can be derived for the 2D case. They involve slant stacks of filtered subsurface offset domain common image gathers followed by the application of the appropriate weighting factors. For the numerical implementation of these formulas, we develop an algorithm based on the true amplitude version of the one‐way paraxial approximation. Finally, we demonstrate the relevance of our approach with a set of applications on synthetic datasets and compare our results with those obtained on the Marmousi model by multi‐arrival ray‐based preserved‐amplitude migration. While results are similar, we observe that our results are less affected by artefacts.     

3D TABLE-DRIVEN MIGRATION1

An efficient full 3D wavefield extrapolation technique is presented. The method can be used for any type of subsurface structure and the degree of accuracy and dip-angle performance are user-defined. The extrapolation is performed in the space-frequency domain as a space-dependent spatial convolution with recursive Kirchhoff extrapolation operators. To get a high level of efficiency the operators are optimized such that they have the smallest possible size for a specified accuracy and dip-angle performance. As both accuracy and maximum dip-angle are input parameters for the operator calculation, the method offers the possibility of a trade-off between these quantities and efficiency. The operators are calculated in advance and stored in a table for a range of wavenumbers. Once they have been calculated they can be used many times. At the basis of the operator design is the well-known phase-shift operator. Although this operator is exact for homogeneous media only, it is assumed that it may be applied locally in case of inhomogeneities. Lateral velocity variations can then be handled by choosing the extrapolation operator according to the local value of the velocity. Optionally the operators can be designed such that they act as spatially variant high-cut filters. This means that the evanescent field can be suppressed in one pass with the extrapolation. The extrapolation method can be used both in prestack and post-stack applications. In this paper we use it in zero-offset migration. Tests on 2D and 3D synthetic and 2D real data show the excellent quality of the method. The full 3D result is much better then the result of two-pass migration, which has been applied to the same data. The implementation yields a code that is fully vectorizable, which makes the method very suitable for vector computers.     

EFFICIENT 2D AND 3D SHOT RECORD REDATUMING1

In order to make 3D prestack depth migration feasible on modern computers it is necessary to use a target-oriented migration scheme. By limiting the output of the migration to a specific depth interval (target zone), the efficiency of the scheme is improved considerably. The first step in such a target-oriented approach is redatuming of the shot records at the surface to the upper boundary of the target zone. For this purpose, efficient non-recursive wavefield extrapolation operators should be generated. We propose a ray tracing method or the Gaussian beam method. With both methods operators can be efficiently generated for any irregular shooting geometry at the surface. As expected, the amplitude behaviour of the Gaussian beam method is better than that of the ray tracing based operators. The redatuming algorithm is performed per shot record, which makes the data handling very efficient. From the shot records at the surface‘genuine zero-offset data’are generated at the upper boundary of the target zone. Particularly in situations with a complicated overburden, the quality of target-oriented zero-offset data is much better than can be reached with a CMP stacking method at the surface. The target-oriented zero-offset data can be used as input to a full 3D zero-offset depth migration scheme, in order to obtain a depth section of the target zone.     

The common focal point and common reflection surface methodologies: subsets or complements?

The common focal point (CFP) method and the common reflection surface (CRS) stack method are compared. The CRS method is a fast, highly automated procedure that provides high S/N ratio simulation of zero‐offset (ZO) images by combining, per image point, the reflection energy of an arc segment that is tangential to the reflector. It uses smooth parametrized two‐way stacking operators, based on a data‐driven triplet of attributes in 2D (eight parameters in 3D). As a spin‐off, the attributes can be used for several applications, such as the determination of the geometrical spreading factor, multiple prediction, and tomographic inversion into a smooth background velocity model. The CFP method aims at decomposing two‐way seismic reflection data into two full‐aperture one‐way propagation operators. By applying an iterative updating procedure in a half‐migrated domain, it provides non‐smooth focusing operators for prestack imaging using only the energy from one focal point at the reflector. The data‐driven operators inhibit all propagation effects of the overburden. The CFP method provides several spin‐offs, amongst which is the CFP matrix related to one focal point, which displays the reflection amplitudes as measured at the surface for each source–receiver pair. The CFP matrix can be used to determine the specular reflection source–receiver pairs and the Fresnel zone at the surface for reflection in one single focal point. Other spin‐offs are the prediction of internal multiples, the determination of reflectivity effects, velocity‐independent redatuming and tomographic inversion to obtain a velocity–depth model. The CFP method is less fast and less automated than the CRS method. From a pointwise comparison of features it is concluded that one method is not a subset of the other, but that both methods can be regarded as being to some extent complementary.     

An integrated method for resolving the seismic complex near-surface problem

We describe an integrated method for solving the complex near‐surface problem in land seismic imaging. This solution is based on an imaging approach and is obtained without deriving a complex near‐surface velocity model. We start by obtaining from the data the kinematics of the one‐way focusing operators (i.e. time‐reversed Green's functions) that describe propagation between the acquisition surface and a chosen datum reflector using the common‐focus‐point technology. The conventional statics solutions obtained from prior information about the near surface are integrated in the initial estimates of the focusing operators. The focusing operators are updated iteratively until the imaging principle of equal traveltime is fulfilled for each subsurface gridpoint of the datum reflector. Therefore, the seismic data is left intact without any application of time shifts, which makes this method an uncommitted statics solution. The focusing operators can be used directly for wave‐equation redatuming to the respective reflector or for prestack imaging if determined for multiple reflecting boundaries. The underlying velocity model is determined by tomographic inversion of the focusing operators while also integrating any hard prior information (e.g. well information). This velocity model can be used to perform prestack depth imaging or to calculate the depth of the new datum level. We demonstrate this approach on 2D seismic data acquired in Saudi Arabia in an area characterized by rugged topography and complex near‐surface geology.     

Incorporation of the Effects of Absorption into Kirchhoff Migration

—The effects of absorption are incorporated into Kirchhoff migration. The aim is to reconstruct the structures with true-amplitude seismic data and to increase the resolution of migrated data. A complex wave velocity is introduced into the solution of the Helmholtz equation, the starting point of Kirchhoff migration. This leads to an additional filter, the antidissipation operator, which is convolved with the wave field. The general structure of Kirchhoff migration remains unchanged. The effects of the antidissipation operator are illustrated on synthetic data. The new operator is valid for complex media with varying velocity and varying quality factor Q. Moreover there is no limitation to constant-Q, frequency-dependent Q can also be handled. The success of anelastic migration depends on how well the Q macro model is known.     

水成物对GNOS掩星弯曲角同化的影响评估

本文在无线电掩星弯曲角射线追踪正演算子中引入水成物的影响,针对台风个例,利用FY-3c GNOS弯曲角资料的同化展开研究.通过分析水成物对掩星弯曲角正演精度的影响,指出当掩星剖面跨越一定厚度的台风区云雨大气时,多相态水成物对GNOS弯曲角正演精度的影响不可忽略.进而提出一种考虑云雨影响的掩星折射率正演算法,将掩星折射率的正演分别在晴空区和云雨区进行,在云雨区正演算子中增加多相态水成物含量对正演掩星折射率的贡献,改进了FY-3c GNOS弯曲角资料在云雨大气环境的同化方案.针对2018年24号台风个例,进行了同化的参照试验、未考虑和考虑水成物影响时GNOS弯曲角的3DVAR同化试验,考量云雨环境下的GNOS弯曲角资料同化对台风模拟的影响差异.试验结果表明,两种同化方案皆能改善台风路径预报,台风中心海平面气压模拟都能接近实际观测,台风最大风速也不同程度增大.而考虑水成物含量的影响后,资料同化能更有效缩小观测空间与背景场空间之间的偏差,同化后观测与分析的偏差更接近高斯分布,台风外围动力场和热力场环境能够得到更优的调整,使得96 h的台风路径模拟平均距离误差较不考虑水成物影响的情形减小了约14...     

ELASTIC REDATUMING OF MULTICOMPONENT SEISMIC DATA1

Elastic redatuming can be carried out before or after decomposition of the multicomponent data into independent PP, PS, SP, and SS responses. We argue that from a practical point of view, elastic redatuming is preferably applied after decomposition. We review forward and inverse extrapolation of decomposed P- and S-wavefields. We use the forward extrapolation operators to derive a model of discrete multicomponent seismic data. This forward model is fully described in terms of matrix manipulations. By applying these matrix manipulations in reverse order we arrive at an elastic processing scheme for multicomponent data in which elastic redatuming plays an essential role. Finally, we illustrate elastic redatuming with a controlled 2D example, consisting of simulated multicomponent seismic data.     

单程波算子积分解的象征表示

单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法,他们代表了算子的不同表示方法,当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式,给出了算子象征的代数运算的头几阶表达式,这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式,对改善Kirchhoff积分的聚焦性能有重要意义.     

基于全波形反演的探地雷达数据逆时偏移成像

逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果.     

结合基准面重建的叠前时间偏移方法

提出了一种结合虚拟界面、瑞利积分和相移法的混合的基准面重建方法.通过与叠前时间偏移方法结合,形成了针对起伏地表采集数据的叠前时间偏移方法和新流程.该方法能正确考虑波在近地表传播的实际路径,克服了高速层出露时静校正方法的误差;它也能自己确定虚拟层速度,避免了现行基于波场延拓的基准面重建方法需要准确近地表速度的困难.文中分别用近地表存在明显低速层和近地表有高速层出露这两类模型的理论数据,验证了所发展方法和流程的有效性.     

双平方根方程三维叠前深度偏移

从双平方根（DSR）形式的波动方程出发，基于沉降观测概念和地震波扰动理论，介绍了深度域的DSR全偏移算子及共成像道集的生成方法. 根据三维地震数据的方位角特征，通过对全偏移算子的稳相近似，依次导出了适应于零方位角道集、Cross line共偏移距道集以及共偏移距矢量道集的偏移算子. 理论分析与合成数据的数值试验表明，DSR全偏移算子、共方位角偏移算子对介质速度变化的适应性很强，而其余两种偏移算子仅适用于缓变速情况.     

Residual dip moveout in VTI media

Dip‐moveout (DMO) correction is often applied to common‐offset sections of seismic data using a homogeneous isotropic medium assumption, which results in a fast execution. Velocity‐residual DMO is developed to correct for the medium‐treatment limitation of the fast DMO. For reasonable‐sized velocity perturbations, the residual DMO operator is small, and thus is an efficient means of applying a conventional Kirchhoff approach. However, the shape of the residual DMO operator is complicated and may form caustics. We use the Fourier domain for the operator development part of the residual DMO, while performing the convolution with common‐offset data in the space–time domain. Since the application is based on an integral (Kirchhoff) method, this residual DMO preserves all the flexibility features of an integral DMO. An application to synthetic and real data demonstrates effectiveness of the velocity‐residual DMO in data processing and velocity analysis.     

Application of the topographic datuming operator to a data setfrom the Eastern Arabian Peninsula

We apply a redatuming methodology, designed to handle rugged topography and the presence of high‐velocity layers near the acquisition surface, to a 2D land seismic data set acquired in Saudi Arabia. This methodology is based on a recently developed prestack operator, which we call the topographic datuming operator (TDO). The TDO, unlike static corrections, allows for the movement of reflections laterally with respect to their true locations, corresponding to the new datum level. Thus, it mitigates mispositioning of events and velocity bias introduced by the assumption of surface consistency and the time‐invariant time shifts brought about by static corrections. Using the shallow velocities estimated from refracted events, the TDO provides a superior continuity of reflections and better focusing than that obtained from conventional static corrections in most parts of the processed 2D line. The computational cost of applying the TDO is only slightly higher than static corrections. The marginal additional computational cost and the possibility of estimating, after TDO redatuming, stacking velocities that are not affected by a spurious positive bias, as in the case of static corrections, are further advantages of the proposed methodology. The likelihood of strong heterogeneities in the most complex part of the line limits the applicability of any approach based upon geometrical optics; however, the TDO produces results that are slightly better than those obtained from static corrections because of its ability to partially collapse diffractions generated in the near surface.     

基于直射线近似的一步法基准面校正方法研究

地形基准面校正算子（Topographic Datuming Operator,以下简称TDO）是一种基于直射线近似得到的基准面延拓算子.TDO可以视为是两步法波动方程基准面校正与常规静校正之间的一种过渡算法,该方法的最大特点在于它可以基于共炮点道集将炮点和检波点同时延拓到给定的水平基准面,因此相对于常规的两步法叠前波动方程基准面校正,TDO方法可以认为是一种更为高效的一步法基准面延拓方法.本文基于理论与实际数据论证了上述观点.     

Internal multiple attenuation by Kirchhoff extrapolation

Surface‐related multiple elimination is the leading methodology for surface multiple removal. This data‐driven approach can be extended to interbed multiple prediction at the expense of a huge increase of the computational burden. This cost makes model‐driven methods still attractive, especially for the three dimensional case. In this paper we present a methodology that extends Kirchhoff wavefield extrapolation to interbed multiple prediction. In Kirchhoff wavefield extrapolation for surface multiple prediction a single round trip to an interpreted reflector is added to the recorded data. Here we show that interbed multiples generated between two interpreted reflectors can be predicted by applying the Kirchhoff wavefield extrapolation operator twice. In the first extrapolation step Kirchhoff wavefield extrapolation propagates the data backward in time to simulate a round trip to the shallower reflector. In the second extrapolation step Kirchhoff wavefield extrapolation propagates the data forward in time to simulate a round trip to the deeper reflector. In the Kirchhoff extrapolation kernel we use asymptotic Green's functions. The prediction of multiples via Kirchhoff wavefield extrapolation is possibly sped up by computing the required traveltimes via a shifted hyperbola approximation. The effectiveness of the method is demonstrated by results on both synthetic and field data sets.     

Fourier finite-difference migration for steeply dipping reflectors with complex overburden1

The improvement in accuracy and efficiency of wave-equation migration techniques is an ongoing topic of research. The main problem is the correct imaging of steeply dipping reflectors in media with strong lateral velocity variations. We propose an improved migration method which is based on cascading phase-shift and finite-difference operators for downward continuation. Due to these cascaded operators we call this method‘Fourier finite-difference migration’(FFD migration). In our approach we try to generalize and improve the split-step Fourier migration method for strong lateral velocity variations using an additional finite-difference correction term. Like most of the current migration methods in use today, our method is based on the one-way wave equation. It is solved by first applying the square-root operator but using a constant velocity at each depth step which has to be the minimum velocity. In a second step, the approximate difference between the correct square-root operator and this constant-velocity squareroot operator (the error made in the first step) is implemented as an implicit FD migration scheme, part of which is the split-step Fourier correction term. Some practical aspects of the new FFD method are discussed. Its performance is compared with that of split-step and standard FD migration schemes. First applications to synthetic and real data sets are presented. They show that the superiority of FFD migration becomes evident by migrating steeply dipping reflectors with complex overburden having strong lateral velocity variations. If velocity is laterally constant, FFD migration has the accuracy of the phase-shift method. The maximum migration angle is velocity adaptive, in contrast to conventional FD migration schemes. It varies laterally depending on the local level of velocity variation. FFD migration is more efficient than higher-order implicit FD schemes. These schemes use two cascaded downward-continuation steps in order to attain comparable migration performance.