Gilles Lambaré

Antileakage Fourier transform for seismic data regularization

Seismic data regularization, which spatially transforms irregularly sampled acquired data to regularly sampled data, is a long-standing problem in seismic data processing. Data regularization can be implemented using Fourier theory by using a method that estimates the spatial frequency content on an irregularly sampled grid. The data can then be reconstructed on any desired grid. Difficulties arise from the nonorthogonality of the global Fourier basis functions on an irregular grid, which results in the problem of “spectral leakage”: energy from one Fourier coefficient leaks onto others. We investigate the nonorthogonality of the Fourier basis on an irregularly sampled grid and propose a technique called “antileakage Fourier transform” to overcome the spectral leakage. In the antileakage Fourier transform, we first solve for the most energetic Fourier coefficient, assuming that it causes the most severe leakage. To attenuate all aliases and the leakage of this component onto other Fourier coefficients, the data component corresponding to this most energetic Fourier coefficient is subtracted from the original input on the irregular grid. We then use this new input to solve for the next Fourier coefficient, repeating the procedure until all Fourier coefficients are estimated. This procedure is equivalent to “reorthogonalizing” the global Fourier basis on an irregularly sampled grid. We demonstrate the robustness and effectiveness of this technique with successful applications to both synthetic and real data examples.

Common-angle migration: A strategy for imaging complex media

Complex velocity models characterized by strong lateral variations are certainly a great motivation, but also a great challenge, for depth imaging. In this context, some unexpected results can occur when using depth imaging algorithms. In general, after a common shot or common offset migration, the resulting depth images are sorted into common‐image gathers (CIG), for further processing such as migration‐based velocity analysis or amplitude‐variation‐with‐offset analysis. In this paper, we show that CIGs calculated by common‐shot or common‐offset migration can be strongly affected by artifacts, even when a correct velocity model is used for the migration. The CIGs are simply not flat, due to unexpected curved events (kinematic artifacts) and strong lateral variations of the amplitude (dynamic artifacts). Kinematic artifacts do not depend on the migration algorithm provided it can take into account lateral variations of the velocity model. This can be observed when migrating the 2‐D Marmousi dataset either with a wave‐equation migration or with a multivalued Kirchhoff migration/inversion. On the contrary, dynamic artifacts are specific to multi‐arrival ray‐based migration/inversion. This approach, which should provide a quantitative estimation of the reflectivity of the model, provides in this context dramatic results. In this paper, we propose an analysis of these artifacts through the study of the ray‐based migration/inversion operator. The artifacts appear when migrating a single‐fold subdata set with multivalued ray fields. They are due to the ambiguous focusing of individual reflected events at different locations in the image. No information is a priori available in the single‐fold data set for selecting the focusing position, while migration of multifold data would provide this information and remove the artifacts by the stack of the CIGs. Analysis of the migration/inversion operator provides a physical condition, the imaging condition, for insuring artifact free CIGs. The specific cases of common‐shot and common‐offset single‐fold gathers are studied. It appears clearly that the imaging condition generally breaks down in complex velocity models for both these configurations. For artifact free CIGs, we propose a novel strategy: compute CIGs versus the diffracting/reflecting angle. Working in the angle domain seems the natural way for unfolding multivalued ray fields, and it can be demonstrated theoretically and practically that common‐angle imaging satisfies the imaging condition in the great majority of cases. Practically, the sorting into angle gathers can not be done a priori over the data set, but is done in the inner depth migration loop. Depth‐migrated images are obtained for each angle range. A canonical example is used for illustrating the theoretical derivations. Finally, an application to the Marmousi model is presented, demonstrating the relevance of the approach.

Practical aspects and applications of 2D stereotomography

Stereotomography is a new velocity estimation method. This tomographic approach aims at retrieving subsurface velocities from prestack seismic data. In addition to traveltimes, the slope of locally coherent events are picked simultaneously in common offset, common source, common receiver, and common midpoint gathers. As the picking is realized on locally coherent events, they do not need to be interpreted in terms of reflection on given interfaces, but may represent diffractions or reflections from anywhere in the image. In the high‐frequency approximation, each one of these events corresponds to a ray trajectory in the subsurface. Stereotomography consists of picking and analyzing these events to update both the associated ray paths and velocity model. In this paper, we describe the implementation of two critical features needed to put stereotomography into practice: an automatic picking tool and a robust multiscale iterative inversion technique. Applications to 2D reflection seismic are presented on syn...

Can we quantitatively image complex structures with rays

Ray‐based prestack depth migration fails to image quantitatively complex structures when a single arrival—for example, the first or the strongest one—is taken into account. In this paper, we investigate whether accounting for multiple arrivals in ray‐based preserved amplitude prestack depth migration allows one to improve quantitative imaging of complex media. The asymptotic ray‐Born migration/inversion, originally designed to process one single arrival, is extended to the case of multiple arrivals by accounting for the cross‐contributions of all the source and receiver raypaths. Multiple arrivals in the folded ray fields are computed by a dynamic ray tracing based on a wavefront construction technique. With an application to the complex Marmousi model, we demonstrate that ray‐Born inversion can provide a reliable quantitative estimation of the relative impedance perturbation even in the complex deep part of the model, for which the amplitudes were underestimated drastically when only a single arrival was...

Migration velocity analysis from locally coherent events in 2‐D laterally heterogeneous media, Part I: Theoretical aspects

We present a new method based on migration velocity analysis (MVA) to estimate 2‐D velocity models from seismic reflection data with no assumption on reflector geometry or the background velocity field. Classical approaches using picking on common image gathers (CIGs) must consider continuous events over the whole panel. This interpretive step may be difficult—particularly for applications on real data sets. We propose to overcome the limiting factor by considering locally coherent events. A locally coherent event can be defined whenever the imaged reflectivity locally shows lateral coherency at some location in the image cube.In the prestack depth‐migrated volume obtained for an a priori velocity model, locally coherent events are picked automatically, without interpretation, and are characterized by their positions and slopes (tangent to the event). Even a single locally coherent event has information on the unknown velocity model, carried by the value of the slope measured in the CIG. The velocity is e...

Migration velocity analysis from locally coherent events in 2‐D laterally heterogeneous media, Part II: Applications on synthetic and real data

We demonstrate a method for estimating 2-D velocity models from synthetic and real seismic reflection data in the framework of migration velocity analysis (MVA). No assumption is required on the reflector geometry or on the unknown background velocity field, provided that the data only contain primary reflections/diffractions. In the prestack depth-migrated volume, locations where the reflectivity exhibits local coherency are automatically picked without interpretation in two panels: common image gathers (CIGs) and common offset gathers (COGs). They are characterized by both their positions and two slopes. The velocity is estimated by minimizing all slopes picked in the CIGs. We test the applicability of the method on a real data set, showing the possibility of an efficient inversion using (1) the migration of selected CIGs and COGs, (2) automatic picking on prior uncorrelated locally coherent events, (3) efficient computation of the gradient of the cost function via paraxial ray tracing from the picked events to the surface, and (4) a gradient-type optimization algorithm for convergence.

3D ray+Born migration/inversion—Part 1: Theory

Prestack ray+Born migration/inversion can be split in two steps : the computation of common image gathers (CIGs) and their weighted stack (the migration stack). The choice of the domain for the CIGs (shot, offset, angle, etc.) has a direct impact on the resolution of the migration stack. This resolution can be studied easily in the frame of ray+Born migration/inversion theory resulting into improved migration/inversion formulas according to the acquisition geometry. This paper is devoted to this analysis in the cases of a simple 2D acquisition and of a 3D swath acquisition, both corresponding to classical data sets from the SEG/EAGE 3D overthrust experiment. We show that the migration formula originally designed for 3D marine acquisition is not adaptable to the 3D swath acquisition. Finally, we propose a new formula for this specific acquisition, which improves the resolution of the final migrated image. The relevance of this new formula is illustrated in the frame of the SEG/EAGE experiment in the companion paper.

Fast migration/inversion with multivalued rayfields: Part I—Method, validation test, and application in 2D to Marmousi

We address the problem of building an efficient algorithm for ray-based migration/inversion in complex media characterized by multivalued traveltime. Since the multivalued ray maps can be computed by wavefront construction, the key point lies in the interpolation of these multivalued coarse maps. For this goal we propose two new strategies. The first one is based on an interpolation of the rayfield after identifying the ray branches on the elementary cells of the coarse grid; the second is based on an extrapolation of the arrivals from the coarse grid points. Both strategies are evaluated in two dimensions with the Marmousi data set, demonstrating the CPU efficiency and accuracy of both approaches. Since they both consider multiarrivals locally, the associated overhead cost with respect to a single-arrival implementation is quite low (twice the computing time of the Marmousi test case). The application also reveals that the new interpolation strategy is more efficient in terms of CPU time compared to the extrapolation strategy but at the expense of a more complex numerical implementation and less accurate results.

3D ray+Born migration/inversion—Part 2: Application to the SEG/EAGE overthrust experiment

The SEG/EAGE overthrust model is a synthetic onshore velocity model that was used to generate several large synthetic seismic data sets using acoustic finite-difference modeling. From this database, several realistic subdata sets were extracted and made available for testing 3D processing methods. For example, classic onshore-type data-acquisition geometries are available such as a swath acquisition, which is characterized by a nonuniform distribution of long offsets with azimuth and midpoints. In this paper, we present an application of 2.5D and 3D ray+Born migration/inversion to several classical data sets from the SEG/EAGE overthrust experiment. The method is formulated as a linearized inversion of the scattered wavefield. The method allows quantitative estimates of short wavelength components of the velocity model. First, we apply a 3D migration/inversion formula formerly developed for marine acquisitions to the swath data set. The migrated sections exhibit significant amplitude artifacts and acquisition footprints, also revealed by the shape of the local spatial resolution filters. From the analysis of these spatial resolution filters, we propose a new formula significantly improving the migrated dip section. We also present 3D migrated results for the strike section and a small 3D target containing a channel. Finally, the applications demonstrate, that the ray+ Born migration formula must be adapted to the acquisition geometry to obtain reliable estimates of the true amplitude of the model perturbations. This adaptation is relatively straightforward in the frame of the ray+Born formalism and can be guided by the analysis of the resolution operator.

Stereotomography – Past, Present and Future

Stereotomography was proposed six years ago for estimating velocity macro-models from seismic reflection data. Initially we aimed at defining a method keeping the advantages of standard travel time tomography but remedying to the difficult interpretative picking.