3D velocity-independent elliptically anisotropic moveout correction

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB129-WB136 ◽  
Author(s):  
William Burnett ◽  
Sergey Fomel

Azimuthal anisotropy or lateral velocity variations cause azimuthal variations in moveout velocity, which can degrade seismic images if handled improperly. In cases in which apparent azimuthally anisotropic moveout is present, a single picked velocity is inadequate to flatten an event on a 3D CMP gather. Conventional velocity-analysis techniques require a significant amount of time and effort, especially in areas where apparent anisotropy is observed. We propose a velocity-independent imaging approach to perform an elliptically anisotropic moveout correction in three dimensions. The velocity-independent approach relies on volumetric local traveltime slopes rather than aggregate velocities and therefore provides an azimuthally flexible description of traveltime geometries throughout the gather. We derive theoretical expressions for extracting the moveout slowness matrix and the angle between the symmetry and acquisition axes as volumetric local attributes. A practical inversion scheme to extract the same parameters is also developed. These parameters are used to solve for moveout slowness as a function of azimuth. Tests on a synthetic common-midpoint (CMP) gather show accurate results for the automatic moveout correction and the inversion scheme. A field data example from west Texas illustrates the application of the automatic moveout correction as a residual moveout.

Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 872-882 ◽  
Author(s):  
Biondo Biondi

I present an unconditionally stable, implicit finite‐difference operator that corrects the constant‐velocity phase‐shift operator for lateral velocity variations. The method is based on the Fourier finite‐difference (FFD) method. Contrary to previous results, my correction operator is stable even when the medium velocity has sharp discontinuities, and the reference velocity is higher than the medium velocity. The stability of the new correction enables the definition of a new downward‐continuation method based on the interpolation of two wavefields: the first wavefield is obtained by applying the FFD correction starting from a reference velocity lower than the medium velocity; the second wavefield is obtained by applying the FFD correction starting from a reference velocity higher than the medium velocity. The proposed Fourier finite‐difference plus interpolation (FFDPI) method combines the advantages of the FFD technique with the advantages of interpolation. A simple and economical procedure for defining frequency‐dependent interpolation weight is presented. When the interpolation step is performed using these frequency‐dependent interpolation weights, it significantly reduces the residual phase error after interpolation, the frequency dispersion caused by the discretization of the Laplacian operator, and the azimuthal anisotropy caused by splitting. Tests on zero‐offset data from the SEG‐EAGE salt data set show that the FFDPI method improves the imaging of a fault reflection with respect to a similar interpolation scheme that uses a split‐step correction for adapting to lateral velocity variations.


Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 164-175 ◽  
Author(s):  
Hans J. Tieman

The most common migration velocity analysis algorithms, iterative profile migration, focusing analysis, and stack power analysis are based on restrictive subsurface assumptions that may cause the methods to break down in the presence of lateral velocity variations. Typically, the subsurface is assumed to have either a constant velocity, or at most, a depth variable velocity. Recent innovations include the use of traveltime inversion philosophy to invert migration measurements of the curvature as a function of offset exhibited by an event following migration. Traveltime inversion makes few assumptions regarding the subsurface, but is a rather unstable process. Thus, an important question is, “Under what conditions do the traditional methods break down and make the use of traveltime inversion methods mandatory?” Marrying the common migration analysis with tomography results in a set of equations that, while useful for generating updates from migration measurements, are too complex for answering the above question. However, by restricting the subsurface to low relief structures and assuming small angle wave propagation, these updating equations can be approximated by forms that are identical to the traditional updating equations, except for a factor that is dependent upon the magnitude, position, and spatial wavelength of potential lateral velocity variations. These simpler equations indicate that even relatively long wavelength anomalies can cause updates to point in the wrong direction and iterative procedures to diverge, and under certain conditions, cause the accuracy of these updates to decrease dramatically.


2017 ◽  
Vol 5 (3) ◽  
pp. SN25-SN32 ◽  
Author(s):  
Ping Wang ◽  
Shouting Huang ◽  
Ming Wang

Complex overburdens often distort reservoir images in terms of structural positioning, stratigraphic resolution, and amplitude fidelity. One prime example of a complex overburden is in the deepwater Gulf of Mexico, where thick and irregular layers of remobilized (i.e., allochthonous) salt are situated above prospective reservoir intervals. The highly variant salt layers create large lateral velocity variations that distort wave propagation and the illumination of deeper reservoir targets. In subsalt imaging, tools such as reflection tomography, full-waveform inversion, and detailed salt interpretation are needed to derive a high-resolution velocity model that captures the lateral velocity variations. Once a velocity field is obtained, reverse time migration (RTM) can be applied to restore structural positioning of events below and around the salt. However, RTM by nature is unable to fully recover the reflectivity for desired amplitudes and resolution. This shortcoming is well-recognized by the imaging community, and it has propelled the emergence of least-squares RTM (LSRTM) in recent years. We have investigated how current LSRTM methods perform on subsalt images. First, we compared the formulation of data-domain versus image-domain least-squares migration, as well as methods using single-iteration approximation versus iterative inversion. Then, we examined the resulting subsalt images of several LSRTM methods applied on the synthetic and field data. Among our tests, we found that image-domain single-iteration LSRTM methods, including an extension of an approximate inverse Hessian method in the curvelet domain, not only compensated for amplitude loss due to poor illumination caused by complex salt bodies, but it also produced subsalt images with fewer migration artifacts in the field data. In contrast, an iterative inversion method showed its potential for broadening the bandwidth in the subsalt, but it was less effective in reducing migration artifacts and noise. Based on our understanding, we evaluated the current state of LSRTM for subsalt imaging.


Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1662-1672 ◽  
Author(s):  
Sergey Fomel

Time‐migration velocity analysis can be performed by velocity continuation, an incremental process that transforms migrated seismic sections according to changes in the migration velocity. Velocity continuation enhances residual normal moveout correction by properly taking into account both vertical and lateral movements of events on seismic images. Finite‐difference and spectral algorithms provide efficient practical implementations for velocity continuation. Synthetic and field data examples demonstrate the performance of the method and confirm theoretical expectations.


Geophysics ◽  
1989 ◽  
Vol 54 (6) ◽  
pp. 718-729 ◽  
Author(s):  
Kamal Al‐Yahya

In conventional seismic processing, velocity analysis is performed by using the normal moveout (NMO) equation which is based on the assumption of flat, horizontal reflectors. Imaging by migration (either before or after stack) is done normally in a subsequent step using these velocities. In this paper, velocity analysis and imaging are combined in one step, and migration itself is used as a velocity indicator. Because, unlike NMO, migration can be formulated for any velocity function, migration‐based velocity analysis methods are capable of handling arbitrary structures, i.e., those with lateral velocity variations. In the proposed scheme, each shot gather (profile) is migrated with an initial depth‐velocity model. Profile migration is implemented in the (x, ω) domain, but the actual implementation of profile migration is not critical, as long as it is not done in a spatial‐wavenumber domain, which would preclude handling of lateral velocity variations. After migration with an initial velocity model, the velocity error is estimated, and the initial velocity model is updated; the process is repeated until convergence is achieved. The velocity analysis is based on the principle that after prestack migration with the correct velocity model, an image in a common‐receiver gather (CRG) is aligned horizontally regardless of structure. The deviation from horizontal alignment is therefore a measure of the error in velocity. If the migration velocity is lower than the velocity of the medium, events curve upward, whereas if the migration velocity is higher than the velocity of the medium, events curve downward.


1999 ◽  
Vol 2 (03) ◽  
pp. 271-280 ◽  
Author(s):  
Ekrem Kasap ◽  
Kun Huang ◽  
Than Shwe ◽  
Dan Georgi

Summary The formation-rate-analysis (FRASM) technique is introduced. The technique is based on the calculated formation rate by correcting the piston rate with fluid compressibility. A geometric factor is used to account for irregular flow geometry caused by probe drawdown. The technique focuses on the flow from formation, is applicable to both drawdown and buildup data simultaneously, does not require long buildup periods, and can be implemented with a multilinear regression, from which near-wellbore permeability, p * and formation fluid compressibility are readily determined. The field data applications indicate that FRA is much less amenable to data quality because it utilizes the entire data set. Introduction A wireline formation test (WFT) is initiated when a probe from the tool is set against the formation. A measured volume of fluid is then withdrawn from the formation through the probe. The test continues with a buildup period until pressure in the tool reaches formation pressure. WFTs provide formation fluid samples and produce high-precision vertical pressure profiles, which, in turn, can be used to identify formation fluid types and locate fluid contacts. Wireline formation testing is much faster compared with the regular pressure transient testing. Total drawdown time for a formation test is just a few seconds and buildup times vary from less than a second (for permeability of hundreds of millidarcy) to half a minute (for permeability of less than 0.1 md), depending on system volume, drawdown rate, and formation permeability. Because WFT tested volume can be small (a few cubic centimeters), the details of reservoir heterogeneity on a fine scale are given with better spatial resolution than is possible with conventional pressure transient tests. Furthermore, WFTs may be preferable to laboratory core permeability measurements since WFTs are conducted at in-situ reservoir stress and temperature. Various conventional analysis techniques are used in the industry. Spherical-flow analysis utilizes early-time buildup data and usually gives permeability that is within an order of magnitude of the true permeability. For p* determination, cylindrical-flow analysis is preferred because it focuses on late-time buildup data. However, both the cylindrical- and spherical-flow analyses have their drawbacks. Early-time data in spherical-flow analysis results in erroneous p* estimation. Late-time data are obtained after long testing times, especially in low-permeability formations; however, long testing periods are not desirable because of potential tool "sticking" problems. Even after extended testing times, the cylindrical-flow period may not occur or may not be detectable on WFTs. When it does occur, permeability estimates derived from the cylindrical-flow period may be incorrect and their validity is difficult to judge. New concepts and analysis techniques, combined with 3-D numerical studies, have recently been reported in the literature.1–7 Three-dimensional numerical simulation studies1–6 have contributed to the diagnosis of WFT-related problems and the improved analysis of WFT data. The experimental studies7 showed that the geometric factor concept is valid for unsteady state probe pressure tests. This study presents the FRA technique8 that can be applied to the entire WFT where a plot for both drawdown and buildup periods renders straight lines with identical slopes. Numerical simulation studies were used to generate data to test both the conventional and the FRA techniques. The numerical simulation data are ideally suited for such studies because the correct answer is known (e.g., the input data). The new technique and the conventional analysis techniques are also applied to the field data and the results are compared. We first review the theory of conventional analysis techniques, then present the FRA technique for combined drawdown and buildup data. A discussion of the numerical results and the field data applications are followed by the conclusions. Analysis Techniques It has been industry practice to use three conventional techniques, i.e., pseudo-steady-state drawdown (PSSDD), spherical and cylindrical-flow analyses, to calculate permeability and p* Conventional Techniques Pseudo-Steady-State Drawdown (PSSDD). When drawdown data are analyzed, it is assumed that late in the drawdown period the pressure drop stabilizes and the system approaches to a pseudo-steady state when the formation flow rate is equal to the drawdown rate. PSSDD permeability is calculated from Darcy's equation with the stabilized (maximum) pressure drop and the flowrate resulting from the piston withdrawal:9–11 $$k {d}=1754.5\left({q\mu \over r {i}\Delta p {{\rm max}}}\right),\eqno ({\rm 1})$$where kd=PSSDD permeability, md. The other parameters are given in Nomenclature.


Author(s):  
Debbie Care ◽  
Shirley Nichols ◽  
Derek Woodfield

The use of low-ionic-strength hydroponic culture and image analysis techniques to discriminate and isolate morphologically distinct, genetically differentiated root types within white clover is described. Advantages of this method include the ability to view the genetic expression of the root systems without the modifying effects of growth in soil, to examine the growth and structure of roots over time, and to store the images for further examination. It is recognised that although the root systems grow in three dimensions, they are constrained to two dimensions by the flatbed scanner. However, the morphological parameters determined by image analysis would not be altered whether this analysis was measured in two or three dimensions. Keywords: image analysis, root morphology, solution culture, Trifolium repens


Author(s):  
D. Amazonas ◽  
R. Aleixo ◽  
J. Schleicher ◽  
J. Costa ◽  
A. Novais ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document