scholarly journals Using constraints to address the instabilities of automated prestack velocity analysis

Geophysics ◽  
1992 ◽  
Vol 57 (3) ◽  
pp. 404-419 ◽  
Author(s):  
Christof Stork ◽  
Robert W. Clayton

Generalized prestack velocity analysis methods that use an automated approach to resolve laterally variable interval velocity fields are beset by a series of problems. The problem of resolving lateral velocity variations has inherent complications that prevent automated methods from being robust enough to be applied routinely to data from a variety of geologic provinces. The use of automated prestack velocity analysis methods will not eliminate the step of carefully producing an initial velocity model derived from regional geologic information and an interpretation of a conventionally processed section. For the methods to regularly produce useful additional information, the unique characteristics of each application must be input into the prestack velocity analysis with the use of inversion constraints. These constraints serve either to adapt the generalized prestack velocity analysis to a focused objective in a particular area or to provide iterative, interpretational tools that help the user produce a velocity model.

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. U21-U28 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan

A new migration velocity analysis is developed by combining the speed of parsimonious prestack depth migration with velocity adjustments estimated within and across common-reflection-point (CRP) gathers. The proposed approach is much more efficient than conventional tomographic velocity analysis because only the traces that contribute to a series of CRP gathers are depth migrated at each iteration. The local interval-velocity adjustments for each CRP are obtained by maximizing the stack amplitude over the predicted (nonhyperbolic) moveout in each CRP gather; this does not involve retracing rays. At every iteration, the velocity in each pixel is updated by averaging over all the predicted velocity updates. Finally, CRP positions and orientations are updated by parsimonious migration, and rays are retraced to define new CRP gathers for the next iteration; this ensures internal consistency between the updated velocity model and the CRP gather. Because the algorithm has a gridded-model parameterization, no explicit representation or fitting of reflectors is involved. Strong lateral-velocity variations, such as those found at salt flanks, can be handled. Application to synthetic and field data sets show that the proposed algorithm works effectively and efficiently.


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.


Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. B241-B252 ◽  
Author(s):  
Daniele Colombo ◽  
Diego Rovetta ◽  
Ersan Turkoglu

Seismic imaging in salt geology is complicated by highly contrasted velocity fields and irregular salt geometries, which cause complex seismic wavefield scattering. Although the imaging challenges can be addressed by advanced imaging algorithms, a fundamental problem remains in the determination of robust velocity fields in high-noise conditions. Conventional migration velocity analysis is often ineffective, and even the most advanced methods for depth-domain velocity analysis, such as full-waveform inversion, require starting from a good initial estimate of the velocity model to converge to a correct result. Nonseismic methods, such as electromagnetics, can help guide the generation of robust velocity models to be used for further processing. Using the multiphysics data acquired in the deepwater section of the Red Sea, we apply a controlled-source electromagnetic (CSEM) resistivity-regularized seismic velocity inversion for enhancing the velocity model in a complex area dominated by nappe-style salt tectonics. The integration is achieved by a rigorous approach of multiscaled inversions looping over model dimensions (1D first, followed by 3D), variable offsets and increasing frequencies, data-driven and interpretation-supported approaches, leading to a hierarchical inversion guided by a parameter sensitivity analysis. The final step of the integration consists of the inversion of seismic traveltimes subject to CSEM model constraints in which a common-structure coupling mechanism is used. Minimization is performed over the seismic data residuals and cross-gradient objective functions without inverting for the resistivity model, which is used as a reference for the seismic inversion (hierarchical approach). Results are demonstrated through depth imaging in which the velocity model derived through CSEM-regularized hierarchical inversion outperforms the results of a seismic-only derived velocity model.


1989 ◽  
Vol 20 (2) ◽  
pp. 301
Author(s):  
P.D. Grant

The Puffin Field is located within the Vulcan Sub-basin of the Timor sea, off the Northwest Coast of Australia. It lies within the offshore exploration permit AC/P2, operated by BHP Petroleum and its co-venturers. It is situated on the Ashmore Platform, an old Triassic horst which is normal faulted against the Swan Graben, a major Mesozoic depocentre and the regional source area. Three wells were drilled in the 1970's. Puffin-1 and Puffin-3 encountered oil in "FIT" tests from within the Maastrichtian 100 ft sand, and Puffin-2 flowed over 4000 barrels of oil per day from a slightly younger 4 m sand. On examination of the results of the Puffin wells, it was evident that there were severe velocity anomalies and differing oil water contacts in the Puffin field. The top of the 100 ft reservoir sand is at 2031.4 m subsea in Puffin-1, 2045 m subsea at Puffin-2 and 2074 m subsea at Puffin-3. The two way times to these events were 1392 ms, 1328 ms and 1398 ms respectively. The interpreted oil water contacts in Puffin-1 and Puffin-3 were 2033 and 2077 ms subsea respectively with no contact seen at Puffin-2. In an attempt to resolve these anomalies the AC/P2 joint venture undertook a detailed seismic reprocessing project of the 1980 data with special emphasis on detailed velocity analysis. This 1987 reprocessing effort involved two passes of velocity filtering and velocity analysis at every 600 m. Velocity analyses were picked on a horizon-consistent basis, such that variations in interval velocity for key horizons could be established for later use in depth conversion. Although sceptical in using stacking functions as the input velocities to depth conversion, they were used, as no viable alternative was feasible. Data quality was reliable to the top of the Palaeocene Calcilutite, and six horizons were picked with their respective velocities to this level. Analysis of the data indicated that the two major units exhibiting interval velocity variation were the Pliocene "low velocity layer" and the Eocene carbonates. Using the smoothed stacking velocity down to the Top Palaeocene Calcilutite the three wells tied the depth conversion with an accuracy of 0.5%. Below this horizon two constant interval velocities were used from well data as the quality of the seismic pick were not as reliable. To verify this model BHPP also undertook a "layer-cake" velocity approach which, although confirming the anomalous zones, could not be used laterally away from the three wells, which unfortunately all lay in a straight line. Two wells, Puffin-4 and Parry-1 were drilled in 1988 to test the resultant interpretation. The wells intersected the Top Palaeocene Calcilutite within 1% of prognosis at Puffin-4 and within 2.2% of prognosis at Parry-1, therefore confirming the stacking velocity model used in depth conversion. However, both wells came in deep to prognosis at the deeper, objective level as a result, in the case of Puffin-4, of being on the downthrown side of a small fault, and at Parry-1 due to a thickening of the Paleocene section and seismic mispicking of the Top Palaeocene Calcilutite. Had the mispick at Parry-1 been avoided then the tie would have been less than 1.0%. Both these mis-interpretations were made in the part of the section where the quality of seismic was poorest. These two results suggest that even though the depth conversion to the Top Paleocene Calcilutite is accurate to within 1%, the magnitude of the velocity variation is larger than the magnitude of the independent depth closure. The Puffin Field requires both better quality seismic below the Base Palaeocene Calcilutite, or the means to resolve the lateral extent and possible thickness of a 4 m sand away from Puffin-2. Until such a method of obtaining either better quality seismic to the objective level, or to be able to define the seismic resolution of the differing sand bodies of a minimum size of 4 m, the Puffin Field will remain a Geophysical enigma.


Geophysics ◽  
1989 ◽  
Vol 54 (2) ◽  
pp. 191-199 ◽  
Author(s):  
John L. Toldi

Conventionally, interval velocities are derived from picked stacking velocities. The velocity‐analysis algorithm proposed in this paper is also based on stacking velocities; however, it eliminates the conventional picking stage by always considering stacking velocities from the point of view of an interval‐velocity model. This view leads to a model‐based, automatic velocity‐analysis algorithm. The algorithm seeks to find an interval‐velocity model such that the stacking velocities calculated from that model give the most powerful stack. An additional penalty is incurred for models that differ in smoothness from an initial interval‐velocity model. The search for the best model is conducted by means of a conjugate‐gradient method. The connection between the interval‐velocity model and the stacking velocities plays an important role in the algorithm proposed in this paper. In the simplest case, stacking velocity is assumed to be equal to rms velocity. For the more general case, a linear theory is developed, connecting interval velocity and stacking velocity through the intermediary of traveltime. When applied to a field data set, the method produces an interval‐velocity model that explains the lateral variation in both stacking velocity and traveltime.


Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. U75-U88 ◽  
Author(s):  
Jintan Li ◽  
William W. Symes

The differential semblance method of velocity analysis flattens image gathers automatically by updating interval velocity to minimize the mean square difference of neighboring traces. We detail an implementation using hyperbolic normal moveout correction as the imaging method. The algorithm is fully automatic, accommodates arbitrary acquisition geometry, and outputs 1D, 2D, or 3D interval velocity models. This variant of differential semblance velocity analysis is effective within the limits of its imaging methodology: mild lateral heterogeneity and data dominated by primary events. Coherent noise events such as multiple reflections tend to degrade the quality of the velocity model estimated by differential semblance. We show how to combine differential semblance velocity analysis with dip filtering to suppress multiple reflections and thus improve considerably the accuracy of the velocity estimate. We illustrate this possibility using multiple-rich data from a 2D marine survey.


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


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE353-VE360 ◽  
Author(s):  
Moshe Reshef

When interval velocity analysis is conducted over complex geologic regions, scattering-angle gathers may cause significant inaccuracies. These inaccuracies are related to the loss of structural dip information when generating common-image gathers (CIGs). In this study, the idea of performing interval velocity analysis in the dip-angle domain was examined and demonstrated with synthetic and real data examples. The effects of migration velocity errors and their identification in this domain were analyzed in detail. Carrying the analysis directly on dip-angle gathers is practically impossible. The ability to perform a standard analysis based on flattening the events in the CIGs is achieved by replacing the dip-angle measure with an equivalent offset measure. This equivalent offset provides higher sensitivity to velocity errors and may improve the accuracy of the resultant velocity model.


Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1331-1339 ◽  
Author(s):  
Tariq Alkhalifah

Prestack migration velocity analysis in the time domain reduces the velocity‐depth ambiguity usually hampering the performance of prestack depth‐migration velocity analysis. In prestack τ migration velocity analysis, we keep the interval velocity model and the output images in vertical time. This allows us to avoid placing reflectors at erroneous depths during the velocity analysis process and, thus, avoid slowing down its convergence to the true velocity model. Using a 1D velocity update scheme, the prestack τ migration velocity analysis performed well on synthetic data from a model with a complex near‐surface velocity. Accurate velocity information and images were obtained using this time‐domain method. Problems occurred only in resolving a thin layer where the low resolution and fold of the synthetic data made it practically impossible to estimate velocity accurately in this layer. This 1D approach also provided us reasonable results for synthetic data from the Marmousi model. Despite the complexity of this model, the τ domain implementation of the prestack migration velocity analysis converged to a generally reasonable result, which includes properly imaging the elusive top‐of‐the‐reservoir layer.


Geophysics ◽  
1993 ◽  
Vol 58 (11) ◽  
pp. 1596-1606 ◽  
Author(s):  
Hans J. Tieman

Stacking velocities can be directly estimated from seismic data without recourse to a multivelocity stack and subsequent search techniques that many current procedures use. This is done as follows: (1) apply NMO to the data (over a window, for a particular common midpoint) using initial estimates for zero offset time and velocity; (2) produce two stacks by summing the data over offset after applying different weighting functions; (3) cross correlate the two stacks; and (4) translate the lag into velocity and time updates. The procedure is iterated until convergence has occurred. Referred to as ARAMVEL (U.S. Patent No. 4,813,027), the method is best implemented as an interactive continuous velocity analysis. Although very simple, both empirical studies and theoretical analysis have shown that it determines velocities more accurately than more traditional approaches based on a scan approach. Convergence is fast, with only one or two iterations usually necessary. The method is robust, as only approximate information is necessary initially. Results with real data show that the method can economically give the detailed velocity control necessary for processing data from areas with considerable lateral velocity variation, as well as provide traveltime information that can be used for sophisticated inversion into interval velocity and depth.


Sign in / Sign up

Export Citation Format

Share Document