Residual stereotomographic inversion

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. E35-E39 ◽  
Author(s):  
Dmitry Neckludov ◽  
Reda Baina ◽  
Evgeny Landa

Depth migration requires highly accurate knowledge of the subsurface velocity field. Different traveltime tomographic methods are used for this purpose. Stereotomography is a tomographic method that uses local dip estimates in addition to traveltimes for velocity model estimation. We present a new methodology for velocity model updating. It combines poststack stereotomography and residual moveout velocity inversion. The former is used for initial model construction and the latter for updating the velocity model. Residual inversion is a kind of stereotomographic inversion applied to common reflection point (CRP) gathers after model-based moveout correction. Velocity analysis can be made more efficient by preselecting the traces that contribute to a series of CRP gathers and using only these traces for inversion. The algorithm is defined in a two-step procedure. First, ray tracing from the reflection point for nonzero reflection offsets defines the source and receiver locations of the data traces in the CRS gather. Then these traces are moveout corrected according to the calculated traveltimes and residual moveout is estimated. The interval velocity model is updated by fitting the velocity that minimizes estimated residuals. Application of the proposed technique demonstrates its robustness and reliability for fast and automatic velocity model estimation.

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 1053-1070 ◽  
Author(s):  
Einar Iversen

The isochron, the name given to a surface of equal two‐way time, has a profound position in seismic imaging. In this paper, I introduce a framework for construction of isochrons for a given velocity model. The basic idea is to let trajectories called isochron rays be associated with iso chrons in an way analogous to the association of conventional rays with wavefronts. In the context of prestack depth migration, an isochron ray based on conventional ray theory represents a simultaneous downward continuation from both source and receiver. The isochron ray is a generalization of the normal ray for poststack map migration. I have organized the process with systems of ordinary differential equations appearing on two levels. The upper level is model‐independent, and the lower level consists of conventional one‐way ray tracing. An advantage of the new method is that interpolation in a ray domain using isochron rays is able to treat triplications (multiarrivals) accurately, as opposed to interpolation in the depth domain based on one‐way traveltime tables. Another nice property is that the Beylkin determinant, an important correction factor in amplitude‐preserving seismic imaging, is closely related to the geometric spreading of isochron rays. For these reasons, the isochron ray has the potential to become a core part of future implementations of prestack depth migration. In addition, isochron rays can be applied in many contexts of forward and inverse seismic modeling, e.g., generation of Fresnel volumes, map migration of prestack traveltime events, and generation of a depth‐domain–based cost function for velocity model updating.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB109-WB118 ◽  
Author(s):  
Jonathan Liu ◽  
Gopal Palacharla

Kirchhoff-type prestack depth migration is the method most popular for outputting offset gathers for velocity-model updating because of its flexibility and efficiency. However, conventional implementations of Kirchhoff migration use only single arrivals. This limits its ability to image complex structures such as subsalt areas. We use the beam methodology to develop a multiarrival Kirchhoff beam migration. The theory and algorithm of our beam migration are analogs to Gaussian beam migration, but we focus on attaining kinematic accuracy and implementation efficiency. The input wavefield of every common offset panel is decomposed into local plane waves at beam centers on the acquisition surface by local slant stacking. Each plane wave contributes a potential single-arrival in Kirchhoff migration. In this way, our method is able to handle multiarrivals caused by model complexity and, therefore, to overcome the limitation of conventional single-arrival Kirchhoff migration. The choice of the width of the beam is critical to the implementation of beam migration. We provide a formula for optimal beam width that achieves both accuracy and efficiency when the velocity model is reasonably smooth. The resulting structural imaging in subsalt and other structurally complex areas is of better quality than that from single-arrival Kirchhoff migration.


Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 533-546 ◽  
Author(s):  
Robert G. Clapp ◽  
Biondo L. Biondi ◽  
Jon F. Claerbout

In areas of complex geology, prestack depth migration is often necessary if we are to produce an accurate image of the subsurface. Prestack depth migration requires an accurate interval velocity model. With few exceptions, the subsurface velocities are not known beforehand and should be estimated. When the velocity structure is complex, with significant lateral variations, reflection‐tomography methods are often an effective tool for improving the velocity estimate. Unfortunately, reflection tomography often converges slowly, to a model that is geologically unreasonable, or it does not converge at all. The large null space of reflection‐tomography problems often forces us to add a sparse parameterization of the model and/or regularization criteria to the estimation. Standard tomography schemes tend to create isotropic features in velocity models that are inconsistent with geology. These isotropic features result, in large part, from using symmetric regularization operators or from choosing a poor model parameterization. If we replace the symmetric operators with nonstationary operators that tend to spread information along structural dips, the tomography will produce velocity models that are geologically more reasonable. In addition, by forming the operators in helical 1D space and performing polynomial division, we apply the inverse of these space‐varying anisotropic operators. The inverse operators can be used as a preconditioner to a standard tomography problem, thereby significantly improving the speed of convergence compared with the typical regularized inversion problem. Results from 2D synthetic and 2D field data are shown. In each case, the velocity obtained improves the focusing of the migrated image.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE269-VE280 ◽  
Author(s):  
Priyank Jaiswal ◽  
Colin A. Zelt

Imaging 2D multichannel land seismic data can be accomplished effectively by a combination of traveltime inversion and prestack depth migration (PSDM), referred to as unified imaging. Unified imaging begins by inverting the direct-arrival times to estimate a velocity model that is used in static corrections and stacking velocity analysis. The interval velocity model (from stacking velocities) is used for PSDM. The stacked data and the PSDM image are interpreted for common horizons, and the corresponding wide-aperture reflections are identified in the shot gathers. Using the interval velocity model, the stack interpretations are inverted as zero-offset reflections to constrain the corresponding interfaces in depth; the interval velocity model remains stationary. We define a coefficient of congruence [Formula: see text] that measures the discrepancy between horizons from the PSDM image andtheir counterparts from the zero-offset inversion. A value of unity for [Formula: see text] implies that the interpreted and inverted horizons are consistent to within the interpretational uncertainties, and the unified imaging is said to have converged. For [Formula: see text] greater than unity, the interval velocity model and the horizon depths are updated by jointly inverting the direct arrivals with the zero-offset and wide-aperture reflections. The updated interval velocity model is used again for both PSDM and a zero-offset inversion. Interpretations of the new PSDM image are the updated horizon depths. The unified imaging is applied to seismic data from the Naga Thrust and Fold Belt in India. Wide-aperture and zero-offset data from three geologically significant horizons are used. Three runs of joint inversion and PSDM are required in a cyclic manner for [Formula: see text] to converge to unity. A joint interpretation of the final velocity model and depth image reveals the presence of a triangle zone that could be promising for exploration.


2014 ◽  
Vol 69 (6) ◽  
Author(s):  
Sudra Irawan ◽  
Sismanto Sismanto ◽  
Adang Sukmatiawan

Seismic data processing is one of the three stages in the seismic method that has an important role in the exploration of oil and gas. Without good data processing, it is impossible to get seismic image cross section for good interpretation. A research using seismic data processing was done to update the velocity model by horizon based tomography method in SBI Field, North West Java Basin. This method reduces error of seismic wave travel time through the analyzed horizon because the existence velocity of high lateral variation in research area. There are three parameters used to determine the accuracy of the resulting interval velocity model, namely, flat depth gathers, semblance residual moveout that coincides with the axis zero residual moveout, and the correspondence between image depth (horizon) with wells marker  (well seismic tie). Pre Stack Depth Migration (PSDM) form interval velocity model and updating using horizon-based tomography method gives better imaging of under-surfaced structure results than PSDM before using tomography. There are three faults found in the research area, two normal faults have southwest-northeast strike and the other has northwest-southeast strike. The thickness of reservoir in SBI field, North West Java Basin, is predicted between 71 to 175 meters and the hydrocarbon (oil) reserve is predicted about  with 22.6% porosity and 70.7% water saturation. 


1999 ◽  
Author(s):  
Robert Bloor ◽  
Alfonso Gonzalez ◽  
Uwe Albertin ◽  
David Yingst

2016 ◽  
Vol 4 (01) ◽  
pp. 63
Author(s):  
Yuninggar Dwi Nugroho ◽  
Sudarmaji S

<span>The input data for pre stack time migration and pre stack depth migration is velocity model. <span>The exact velocity model can provide maximum result in seismic section. The best seismic <span>section can minimize possibility of errors during interpretation. Model based and grid based <span>tomography are used to refine the interval velocity model. The interval velocity will be used as <span>input in the pre stack depth migration. Initial interval velocity is obtained from RMS velocity<br /><span>using Dix formula. This velocity will be refined by global depth tomography method. The <span>global depth tomography method is divided into model based and grid based tomography. <span>Velocity analysis is performed along the horizon (depth model). Residual depth move out is <span>obtained from picking velocity. It is used as input in tomography method. The flat gather is <span>obtained at tenth iteration. The interval velocity that is obtained from tenth iteration has the <span>small errors. Tomography method can provide maximum result on velocity refinement. That is <span>shown by the result that the pre stack depth migration is much better than using initial interval <span>velocity. The pull up effect can be corrected by tomography method.</span></span></span></span></span></span></span></span></span></span></span></span><br /></span>


2016 ◽  
Vol 28 (2) ◽  
pp. 43
Author(s):  
Hagayudha Timotius ◽  
Yulinar Firdaus

The main goal of seismic exploration is to get an accurate image of subsurface section so it can be easily interpreted. Pre Stack Depth Migration (PSDM) is such a powerful imaging tool especially for complex area such an area where strong lateral velocity variations exist. The main challenge of PSDM is the need of accurate interval velocity model.In this research, Dix Transformation, coherency inversion, and tomography are used for initial interval velocity model, and then tomography is used for interval velocity model refinement. We compare also between seismic image resulted from PSDM and PSTM to determine the best method. The seismic data that processed in this paper is derived from north western part of Australian Waters. Kata kunci: Pre Stack Depth Migration, Dix Transformation, coherency inversion, tomography. Tujuan utama dari eksplorasi seismik adalah menghasilkan citra yang akurat dari penampang bawah permukaan sehingga diinterpretasi lebih mudah. Pre Stack Depth Migration (PSDM) merupakan suatu metode yang memberikan hasil peningkatan kualitas citra seismik pada daerah kompleks dimana terjadi variasi kecepatan lateral yang signifikan. Salah satu syarat penting yang harus dipenuhi agar hasil PSDM lebih optimal adalah model kecepatan interval yang akurat. Dalam penelitian ini Transformasi Dix, inversi koheren, dan tomografi digunakan untuk memenuhi syarat tersebut. Perbandingan hasil penampang seimik PSDM dan PSTM dilakukan untuk menentukan metode terbaik. Data seismik yang diolah dalam tulisan ini berasal dari wilayah Perairan Baratlaut Australia. Kata kunci: Pre Stack Depth Migration, Transformasi Dix, inversi koheren, tomografi


2010 ◽  
Vol 50 (2) ◽  
pp. 723
Author(s):  
Sergey Birdus ◽  
Erika Angerer ◽  
Iftikhar Abassi

Processing of multi and wide-azimuth seismic data faces some new challenges, and one of them is depth-velocity modelling and imaging with azimuthal velocity anisotropy. Analysis of multi-azimuth data very often reveals noticeable fluctuations in moveout between different acquisition directions. They can be caused by several factors: real azimuthal interval velocity anisotropy associated with quasi-vertical fractures or present day stress field within the sediments; short-wavelength velocity heterogeneities in the overburden; TTI (or VTI) anisotropy in the overburden; or, random distortions due to noise, multiples, irregularities in the acquisition geometry, etcetera. In order to build a velocity model for multi-azimuth pre-stack depth migration (MAZ PSDM) taking into account observed azimuthal anisotropy, we need to recognise, separate and estimate all the effects listed above during iterative depth-velocity modelling. Analysis of seismic data from a full azimuth 3D seismic land survey revealed the presence of strong spatially variable azimuthal velocity anisotropy that had to be taken into consideration. Using real data examples we discuss major steps in depth processing workflow that took such anisotropy into account: residual moveout estimation in azimuth sectors; separation of different effects causing apparent azimuthal anisotropy (see A–D above); iterative depth-velocity modelling with azimuthal anisotropy; and, subsequent MAZ anisotropic PSDM. The presented workflow solved problems with azimuthal anisotropy in our multi-azimuth dataset. Some of the lessons learned during this MAZ project are relevant to every standard narrow azimuth seismic survey recorded in complex geological settings.


Geophysics ◽  
1986 ◽  
Vol 51 (5) ◽  
pp. 1087-1109 ◽  
Author(s):  
N. D. Whitmore ◽  
Larry R. Lines

Vertical seismic profiles (VSPs) can supply information about both velocity and subsurface interface locations. Properly designed VSPs can be used to map steeply dipping interfaces such as salt dome flanks. Mapping subsurface interfaces with VSP data requires careful survey design, appropriate data processing, interval velocity estimation, and reflector mapping. The first of these four ingredients is satisfied, in most cases, by preacquisition modeling. The second is accomplished by careful data processing. Initial velocity estimates are provided by seismic tomography. Velocity‐model refinement is accomplished by a combination of iterative modeling and iterative least‐squares inversion. Finally, the resultant interval velocities are used in depth migration of the processed VSP. These four ingredients have been combined to map a salt dome flank.


Sign in / Sign up

Export Citation Format

Share Document