Wide‐azimuth angle gathers for anisotropic wave‐equation migration

Extended common‐image‐point gathers (CIP) constructed by wide‐azimuth TI wave‐equation migration contain all the necessary information for angle decomposition as a function of the reflection and azimuth angles at selected locations in the subsurface. The aperture and azimuth angles are derived from the extended images using analytic relations between the space‐ and time‐lag extensions using information which is already available at the time of migration, i.e. the anisotropic model parameters. CIPs are cheap to compute because they can be distributed in the image at the most relevant positions, as indicated by the geologic structure. If the reflector dip is known at the CIP locations, then the computational cost can be reduced by evaluating only two components of the space‐lag vector. The transformation from extended images to angle gathers is a planar Radon transform which depends on the local medium parameters. This transformation allows us to separate all illumination directions for a given experiment, or between different experiments. We do not need to decompose the reconstructed wavefields or to choose the most energetic directions for decomposition. Applications of the method include illumination studies in complex areas where ray‐based methods fail, and assuming that the subsurface illumination is sufficiently dense, the study of amplitude variation with aperture and azimuth angles.     

基于倾角-偏移距域道集的绕射波成像

地震绕射波是地下非连续性地质体的地震响应,绕射波成像对地下断层、尖灭和小尺度绕射体的识别具有重要的意义.在倾角域共成像点道集中,反射波同相轴表现为一条下凸曲线,能量主要集中在菲涅耳带内,绕射波能量则比较发散.由于倾角域菲涅耳带随偏移距变化而存在差异,因此本文提出一种在倾角-偏移距域道集中精确估计菲涅耳带的方法,在各偏移距的倾角域共成像点道集中实现菲涅耳带的精确切除,从而压制反射波.在倾角-偏移距域道集中还可以分别实现绕射波增强,绕射波同相轴相位校正,因此能量弱的绕射波可以清晰地成像.在倾角域共成像点道集中,反射波同相轴的最低点对应于菲涅耳带估计所用的倾角,因此本文提出一种在倾角域共成像点道集中直接自动拾取倾角场的方法.理论与实际资料试算验证了本文绕射波成像方法的有效性.     

Migration to multiple offset and velocity analysis1

The stacking velocity best characterizes the normal moveout curves in a common-mid-point gather, while the migration velocity characterizes the diffraction curves in a zero-offset section as well as in a common-midpoint gather. For horizontally layered media, the two velocity types coincide due to the conformance of the normal and the image ray. In the case of dipping subsurface structures, stacking velocities depend on the dip of the reflector and relate to normal rays, but with a dip-dependent lateral smear of the reflection point. After dip-moveout correction, the stacking velocities are reduced while the reflection-point smear vanishes, focusing the rays on the common reflection points. For homogeneous media the dip-moveout correction is independent of the actual velocity and can be applied as a dip-moveout correction to multiple offset before velocity analysis. Migration to multiple offset is a prestack, time-migration technique, which presents data sets which mimic high-fold, bin-centre adjusted, common-midpoint gathers. This method is independent of velocity and can migrate any 2D or 3D data set with arbitrary acquisition geometry. The gathers generated can be analysed for normal-moveout velocities using traditional methods such as the interpretation of multivelocity-function stacks. These stacks, however, are equivalent to multi-velocity-function time migrations and the derived velocities are migration velocities.     

Estimation of geological dip and curvature from time‐migrated zero‐offset reflections in heterogeneous anisotropic media

Starting from a given time‐migrated zero‐offset data volume and time‐migration velocity, recent literature has shown that it is possible to simultaneously trace image rays in depth and reconstruct the depth‐velocity model along them. This, in turn, allows image‐ray migration, namely to map time‐migrated reflections into depth by tracing the image ray until half of the reflection time is consumed. As known since the 1980s, image‐ray migration can be made more complete if, besides reflection time, also estimates of its first and second derivatives with respect to the time‐migration datum coordinates are available. Such information provides, in addition to the location and dip of the reflectors in depth, also an estimation of their curvature. The expressions explicitly relate geological dip and curvature to first and second derivatives of reflection time with respect to time‐migration datum coordinates. Such quantitative relationships can provide useful constraints for improved construction of reflectors at depth in the presence of uncertainty. Furthermore, the results of image‐ray migration can be used to verify and improve time‐migration algorithms and can therefore be considered complementary to those of normal‐ray migration. So far, image‐ray migration algorithms have been restricted to layered models with isotropic smooth velocities within the layers. Using the methodology of surface‐to‐surface paraxial matrices, we obtain a natural extension to smooth or layered anisotropic media.     

Parsimonious 3D post-stack Kirchhoff depth migration

Parsimonious post‐stack migration is extended to three dimensions. By tracing single rays back along each incident wave direction (as determined by a local slant stack at the receivers), the ray tracing can be embedded in the migration. This approach significantly reduces the computer time and disk space needed because it is not necessary to build and save image time maps; 3D migration can be performed on a workstation or personal computer rather than using a supercomputer or cluster. The location of a reflector in the output image is defined by tracing a zero‐offset ray to the one‐way traveltime (the image condition); the orientation of the reflector is defined as a surface perpendicular to the raypath. The migration impulse response operator is confined to the first Fresnel zone around the estimated reflection point, which is much smaller than the large isochronic surface in traditional Kirchhoff depth migration. Additional efficiency is obtained by applying an amplitude threshold to reduce the amount of data to be migrated. Tests on synthetic data show that the proposed implementation of parsimonious 3D post‐stack Kirchhoff depth migration is at least two orders of magnitude faster than traditional Kirchhoff migration, at the expense of slightly degraded migration image coherence. The proposed migration is expected to be a useful complement to conventional time migrations for fast initial imaging of subsurface structures and for real‐time imaging of near‐offset sections during data acquisition for quality control.     

What is DMO coverage?

'Coverage' or 'fold' is defined as the multiplicity of common-midpoint (CMP) data. For CMP stacking the coverage is consistent with the number of traces sharing a common reflection point on flat subsurface reflectors. This relationship is not true for dipping reflectors. The deficiencies of CMP stacking with respect to imaging dipping events have long been overcome by the introduction of the dip-moveout (DMO) correction. However, the concept of coverage has not yet satisfactorily been updated to a 'DMO coverage' consistent with DMO stacking. A definition of constant-velocity DMO coverage will be proposed here. A subsurface reflector will be illuminated from a given source and receiver location if the time difference between the reflector zero-offset traveltime and the NMO- and DMO-corrected traveltime of the reflection event is less than half a dominant wavelength. Due to the fact that a subsurface reflector location is determined by its zero-offset traveltime, its strike and its dip, the DMO coverage also depends on these three parameters. For every surface location, the proposed DMO coverage consists of a 3D fold distribution over reflector strike, dip and zero-offset traveltime.     

3D description and inversion of reflection moveout of PS‐waves in anisotropic media

Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t²(x²) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (V_P0 and V_S0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (V_P0, V_S0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves.     

Passive seismic data primary estimation and noise removal via focal‐denoising closed‐loop surface‐related multiple elimination based on 3D L1‐norm sparse inversion

Passive seismic has recently attracted a great deal of attention because non‐artificial source is used in subsurface imaging. The utilization of passive source is low cost compared with artificial‐source exploration. In general, constructing virtual shot gathers by using cross‐correlation is a preliminary step in passive seismic data processing, which provides the basis for applying conventional seismic processing methods. However, the subsurface structure is not uniformly illuminated by passive sources, which leads to that the ray path of passive seismic does not fit the hyperbolic hypothesis. Thereby, travel time is incorrect in the virtual shot gathers. Besides, the cross‐correlation results are contaminated by incoherent noise since the passive sources are always natural. Such noise is kinematically similar to seismic events and challenging to be attenuated, which will inevitably reduce the accuracy in the subsequent process. Although primary estimation for transient‐source seismic data has already been proposed, it is not feasible to noise‐source seismic data due to the incoherent noise. To overcome the above problems, we proposed to combine focal transform and local similarity into a highly integrated operator and then added it into the closed‐loop surface‐related multiple elimination based on the 3D L1‐norm sparse inversion framework. Results proved that the method was capable of reliably estimating noise‐free primaries and correcting travel time at far offsets for a foresaid virtual shot gathers in a simultaneous closed‐loop inversion manner.     

频率-空间域有限差分法叠前深度偏移

为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果.     

Diffraction imaging by prestack reverse‐time migration in the dip‐angle domain

Prestack image volumes may be decomposed into specular and non‐specular parts by filters defined in the dip‐angle domain. For space‐shift extended image volumes, the dip‐angle decomposition is derived via local Radon transform in depth and midpoint coordinates, followed by an averaging over space‐shifts. We propose to employ prestack space‐shift extended reverse‐time migration and dip‐angle decomposition for imaging small‐scale structural elements, considered as seismic diffractors, in models with arbitrary complexity. A suitable design of a specularity filter in the dip‐angle domain rejects the dominant reflectors and enhances diffractors and other non‐specular image content. The filter exploits a clear discrimination in dip between specular reflections and diffractions. The former are stationary at the specular dip, whereas the latter are non‐stationary without a preferred dip direction. While the filtered image volume features other than the diffractor images (for example, noise and truncation artefacts are also present), synthetic and field data examples suggest that diffractors tend to dominate and are readily recognisable. Averaging over space‐shifts in the filter construction makes the reflectors? rejection robust against migration velocity errors. Another consequence of the space‐shift extension and its angle‐domain transforms is the possibility of exploring the image in a multiple set of common‐image gathers. The filtered diffractions may be analysed simultaneously in space‐shift, scattering‐angle, and dip‐angle image gathers by means of a single migration job. The deliverables of our method obviously enrich the processed material on the interpreter's desk. We expect them to further supplement our understanding of the Earth's interior.     

地震干涉技术被动源地震成像

被动源地震成像是基于地下反射波响应和透射波响应之间的关系,通过在地表接收由地下非人工震源发出的透射波响应,利用互相关运算来合成反射波响应的方法.合成的反射波响应,又称虚炮集记录等价于地表地震剖面记录.本文研究了由地下随机分布的被动源发出的透射波记录来合成虚炮集记录的基本原理,推导了反射波记录和透射波记录的关系公式;并通过数值模拟方法,分析了被动源成像的可行性.结果表明:利用被动源地震信息进行成像是可行的.     

Depth migration anisotropy analysis in the time domain

In areas of complex geology such as the Canadian Foothills, the effects of anisotropy are apparent in seismic data and estimation of anisotropic parameters for use in seismic imaging is not a trivial task. Here we explore the applicability of common‐focus point (CFP)‐based velocity analysis to estimate anisotropic parameters for the variably tilted shale thrust sheet in the Canadian Foothills model. To avoid the inherent velocity‐depth ambiguity, we assume that the elastic properties of thrust‐sheet with respect to transverse isotropy symmetry axis are homogeneous, the reflector below the thrust‐sheet is flat, and that the anisotropy is weak. In our CFP approach to velocity analysis, for a poorly imaged reflection point, a traveltime residual is obtained as the time difference between the focusing operator for an assumed subsurface velocity model and the corresponding CFP response obtained from the reflection data. We assume that this residual is due to unknown values for anisotropy, and we perform an iterative linear inversion to obtain new model parameters that minimize the residuals. Migration of the data using parameters obtained from our inversion results in a correctly positioned and better focused reflector below the thrust sheet. For traveltime computation we use a brute force mapping scheme that takes into account weakly tilted transverse isotropy media. For inversion, the problem is set up as a generalized Newton's equation where traveltime error (differential time shift) is linearly dependent on the parameter updates. The iterative updates of parameters are obtained by a least‐squares solution of Newton's equations. The significance of this work lies in its applicability to areas where transverse isotropy layers are heterogeneous laterally, and where transverse isotropy layers are overlain by complex structures that preclude a moveout curve fitting.     

偏移成像中的垂直分辨率

Kirchhoff真振幅偏移研究表明,偏移对地震子波具有拉伸作用,拉伸程度与地震波速度、反射角和地层倾角有关.子波拉伸将改变子波宽度(子波延续长度),对垂直分辨率有重要影响.本文通过对衰减Sinc子波的实际计算,分析了地震波速度、反射角和地层倾角分别与子波宽度的关系.从分析来看,陡倾角和大反射角（大炮检距）造成偏移成像垂直分辨率严重下降.高主频子波对各参数均有较低敏感度.     

3D3C陆地反射PS转换波CMP成像影响因素分析

三维三分量(3D3C)陆地反射PS转换波共中心点(CMP)叠加成像方法,虽然抽道集简单,但是对实际资料处理结果往往不理想.尤其当反射界面为三维倾斜界面时,其成像质量较差.本文提出有三个主要因素影响其成像质量:第一,转换点离散.运用实例计算得出,转换点离散度随着纵横波速度比、偏移距和界面倾角的增大而增大.相同界面倾角,不同测线方位的转换点离散度不同,视倾角的绝对值越大离散度也越大;第二,道集内静校正量差异增大.CMP道集中,由于转换点离散使得转换点横向跨度较大,经倾斜界面反射转换的S波出射到近地表地层时的角度差异也较大,导致静校突出;第三,加大动校叠加复杂性.三维倾斜界面PS波CMP道集近炮检距时距方程可表示为双曲形式,但是曲线的顶点位置和动校速度同时随测线方位变化,使得CMP道集同相轴很难校平,动校叠加过程很复杂.     

Pre-stack Kirchhoff depth migration local space-shift imaging condition: synthetic and data examples

Extracting accurate common image angle gathers from pre-stack depth migrations is important in the generation of any incremental uplift to the amplitude versus angle attributes and seismic inversions that can lead to significant impacts in exploration and development success. The commonly used Kirchhoff migration outputs surface common offset image gathers that require a transformation to angle gathers for amplitude versus angle analysis. The accuracy of this transformation is one of the factors that determine the robustness of the amplitude versus angle measurements. Here, we investigate the possibility of implementing an extended imaging condition, focusing on the space-lag condition, for generating subsurface reflection angle gathers within a Kirchhoff migration. The objective is to determine if exploiting the spatial local shift imaging condition can provide any increase in angle gather fidelity relative to the common offset image gathers. The same restrictions with a ray-based approach will apply using the extended imaging condition as both the offset and extended imaging condition method use travel times derived from solutions to an Eikonal equation. The aims are to offer an alternative ray-based method to generate subsurface angle gathers and to understand the impact on the amplitude versus angle response. To this end, the implementation of the space-shift imaging condition is discussed and results of three different data sets are presented. A layered three-dimensional model and a complex two-dimensional model are used to assess the space shift image gathers output from such a migration scheme and to evaluate the seismic attributes relative to the traditional surface offset common image gathers. The synthetic results show that the extended imaging condition clearly provides an uplift in the measured amplitude versus angle over the surface offset migration. The noise profile post-migration is also improved for the space-lag migration due to the double summation inside the migration. Finally, we show an example of a space-lag gather from deep marine data and compare the resultant angle gathers with those generated from an offset migration and a time-shift imaging condition Kirchhoff migration. The comparison of the real data with a well log shows that the space-lag result is a better match to the well compared to the time-lag extended imaging condition and the common offset Kirchhoff migration. Overall, the results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers with an incremental improvement in amplitude versus angle fidelity and lower noise but comes at a higher computational cost.     

Eliminating 3D pre‐stack migration artefacts by 6D filtering

State‐of‐the‐art 3D seismic acquisition geometries have poor sampling along at least one dimension. This results in coherent migration noise that always contaminates pre‐stack migrated data, including high‐fold surveys, if prior‐to‐migration interpolation was not applied. We present a method for effective noise suppression in migrated gathers, competing with data interpolation before pre‐stack migration. The proposed technique is based on a dip decomposition of common‐offset volumes and a semblance‐type measure computation via offset for all constant‐dip gathers. Thus the processing engages six dimensions: offset, inline, crossline, depth, inline dip, and crossline dip. To reduce computational costs, we apply a two‐pass (4D in each pass) noise suppression: inline processing and then crossline processing (or vice versa). Synthetic and real‐data examples verify that the technique preserves signal amplitudes, including amplitude‐versus‐offset dependence, and that faults are not smeared.     

Handling dip discrimination phenomenon in common‐reflection‐surface stack via combination of output‐imaging‐scheme and migration/demigration

In the application of a conventional common‐reflection‐surface (CRS) stack, it is well‐known that only one optimum stacking operator is determined for each zero‐offset sample to be simulated. As a result, the conflicting dip situations are not taken into account and only the most prominent event contributes to any a particular stack sample. In this paper, we name this phenomenon caused by conflicting dip problems as ‘dip discrimination phenomenon’. This phenomenon is not welcome because it not only leads to the loss of weak reflections and tips of diffractions in the final zero‐offset‐CRS stacked section but also to a deteriorated quality in subsequent migration. The common‐reflection‐surface stack with the output imaging scheme (CRS‐OIS) is a novel technique to implement a CRS stack based on a unified Kirchhoff imaging approach. As far as dealing with conflicting dip problems is concerned, the CRS‐OIS is a better option than a conventional CRS stack. However, we think the CRS‐OIS can do more in this aspect. In this paper, we propose a workflow to handle the dip discrimination phenomenon based on a cascaded implementation of prestack time migration, CRS‐OIS and prestack time demigration. Firstly, a common offset prestack time migration is implemented. Then, a CRS‐OIS is applied to the time‐migrated common offset gather. Afterwards, a prestack time demigration is performed to reconstruct each unmigrated common offset gather with its reflections being greatly enhanced and diffractions being well preserved. Compared with existing techniques dealing with conflicting dip problems, the technique presented in this paper preserves most of the diffractions and accounts for reflections from all possible dips properly. More importantly, both the post‐stacked data set and prestacked data set can be of much better quality after the implementation of the presented scheme. It serves as a promising alternative to other techniques except that it cannot provide the typical CRS wavefield attributes. The numerical tests on a synthetic Marmousi data set and a real 2D marine data set demonstrated its effectiveness and robustness.     

Enabling affordable omnidirectional subsurface extended image volumes via probing

Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix–vector formulation from which different image gathers can be formed and with which amplitude‐versus‐angle and wave‐equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two‐way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave‐equation‐based migration‐velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips.