Filtering coherent noise during prestack depth migration

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1054-1066 ◽  
Author(s):  
Bertrand Duquet ◽  
Kurt J. Marfurt

We can often suppress short‐period multiples by predictive deconvolution. We can often suppress coherent noise with significantly different moveout by time‐invariant dip filtering on common‐shot, common‐receiver or NMO-corrected common‐midpoint gathers. Unfortunately, even time variant dip filtering on NMO-corrected data breaks down in the presence of strong lateral velocity variation where the underlying NMO correction breaks down. Underattenuated multiples, converted waves, and diffracted head waves can significantly impede and/or degrade prestack migration‐driven velocity analysis and amplitude variation with offset analysis as well as the quality of the final stacked image. Generalization of time‐variant dip filtering based on conventional NMO corrections of common‐midpoint gathers also breaks down for less conventional data processing situations where we wish to enhance data having nonhyperbolic moveout, such as converted wave energy or long‐offset P-wave reflections in structurally deformed anisotropic media. We present a methodology that defines a depth‐variant velocity filter based on an approximation to the true velocity/depth structure of the earth developed by the interpreter/processor during the normal course of their prestack imaging work flow. Velocity filtering in the depth domain requires the design and calibration of two new least‐squares transforms: a constrained least‐squares common offset Kirchhoff depth migration transform and a transform in residual migration‐velocity moveout space. Each of these new least‐squares transforms can be considered to be generalizations of the well‐known discrete Radon transform commonly used in the oil and gas exploration industry.

Geophysics ◽  
1982 ◽  
Vol 47 (6) ◽  
pp. 884-897 ◽  
Author(s):  
Walter S. Lynn ◽  
Jon F. Claerbout

In areas of large lateral variations in velocity, stacking velocities computed on the basis of hyperbolic moveout can differ substantially from the actual root mean square (rms) velocities. This paper addresses the problem of obtaining rms or migration velocities from stacking velocities in such areas. The first‐order difference between the stacking and the vertical rms velocities due to lateral variations in velocity are shown to be related to the second lateral derivative of the rms slowness [Formula: see text]. Approximations leading to this relation are straight raypaths and that the vertical rms slowness to a given interface can be expressed as a second‐order Taylor series expansion in the midpoint direction. Under these approximations, the effect of the first lateral derivative of the slowness on the traveltime is negligible. The linearization of the equation relating the stacking and true velocities results in a set of equations whose inversion is unstable. Stability is achieved, however, by adding a nonphysical fourth derivative term which affects only the higher spatial wavenumbers, those beyond the lateral resolution of the lateral derivative method (LDM). Thus, given the stacking velocities and the zero‐offset traveltime to a given event as a function of midpoint, the LDM provides an estimate of the true vertical rms velocity to that event with a lateral resolution of about two mute zones or cable lengths. The LDM is applicable when lateral variations of velocity greater than 2 percent occur over the mute zone. At variations of 30 percent or greater, the internal assumptions of the LDM begin to break down. Synthetic models designed to test the LDM when the different assumptions are violated show that, in all cases, the results are not seriously affected. A test of the LDM on field data having a lateral velocity variation caused by sea floor topography gives a result which is supported by depth migration.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. U87-U96 ◽  
Author(s):  
Mamoru Takanashi ◽  
Ilya Tsvankin

One of the most serious problems in anisotropic velocity analysis is the trade-off between anisotropy and lateral heterogeneity, especially if velocity varies on a scale smaller than the maximum offset. We have developed a P-wave MVA (migration velocity analysis) algorithm for transversely isotropic (TI) models that include layers with small-scale lateral heterogeneity. Each layer is described by constant Thomsen parameters [Formula: see text] and [Formula: see text] and the symmetry-direction velocity [Formula: see text] that varies as a quadratic function of the distance along the layer boundaries. For tilted TI media (TTI), the symmetry axis is taken orthogonal to the reflectors. We analyzed the influence of lateral heterogeneity on image gathers obtained after prestack depth migration and found that quadratic lateral velocity variation in the overburden can significantly distort the moveout of the target reflection. Consequently, medium parameters beneath the heterogeneous layer(s) are estimated with substantial error, even when borehole information (e.g., check shots or sonic logs) is available. Because residual moveout in the image gathers is highly sensitive to lateral heterogeneity in the overburden, our algorithm simultaneously inverts for the interval parameters of all layers. Synthetic tests for models with a gently dipping overburden demonstrate that if the vertical profile of the symmetry-direction velocity [Formula: see text] is known at one location, the algorithm can reconstruct the other relevant parameters of TI models. The proposed approach helps increase the robustness of anisotropic velocity model-building and enhance image quality in the presence of small-scale lateral heterogeneity in the overburden.


2017 ◽  
Vol 5 (3) ◽  
pp. SK179-SK187 ◽  
Author(s):  
Thang Ha ◽  
Kurt Marfurt

The Panhandle-Hugoton field, of Texas, Oklahoma, and Kansas, is a giant oil field and is the largest conventional gas field in North America. Most hydrocarbon production in this field comes from the Wichita Uplift area, where the basement is the most shallow. Although the field has been extensively produced, many local hydrocarbon accumulations have not been fully exploited. Recent drilling activity in the survey indicates that some wells produce directly from basement fractures, suggesting a new play type for the area. Because the target is shallow, the seismic data are heavily contaminated by coherent noise, such as ground roll and head waves, creating challenges for seismic processing. To improve the seismic interpretation, we carefully reprocessed the field gathers resulting in improved correlation within the sedimentary and the basement sections. Correlating well control to seismic attribute volumes indicates that a fractured basement gives rise to lower P-wave impedance and strong amplitude versus azimuth anomalies. The azimuthal anisotropy is strongest in a direction parallel to the regional maximum horizontal stress, suggesting that these fractures are open. Coherence anomalies indicate a rugose basement surface, whereas curvature shows two lineament sets, consistent with the weathering and fractured exposure of basement in the Wichita Mountains to the southeast.


1970 ◽  
Vol 19 (4) ◽  
pp. 439-445 ◽  
Author(s):  
G. Willey ◽  
J. R. Cleary ◽  
P. D. Marshall
Keyword(s):  
P Wave ◽  

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 1052-1059 ◽  
Author(s):  
Richard S. Lu ◽  
Dennis E. Willen ◽  
Ian A. Watson

The large velocity contrast between salt and the surrounding sediments generates strong conversions between P‐ and S‐wave energy. The resulting converted events can be noise on P‐wave migrated images and should be identified and removed to facilitate interpretation. On the other hand, they can also be used to image a salt body and its adjacent sediments when the P‐wave image is inadequate. The converted waves with smaller reflection and transmission angles and much larger critical angles generate substantially different illumination than does the P‐wave. In areas where time migration is valid, the ratio between salt thickness in time and the time interval between the P‐wave and the converted‐wave salt base on a time‐migrated image is about 2.6 or 1.3, depending upon whether the seismic wave propagates along one or both of the downgoing and upcoming raypaths in salt as the S‐wave, respectively. These ratios can be used together with forward seismic modeling and 2D prestack depth migration to identify the converted‐wave base‐of‐salt (BOS) events in time and depth and to correctly interpret the subsalt sediments. It is possible to mute converted‐wave events from prestack traces according to their computed arrival times. Prestack depth migration of the muted data extends the updip continuation of subsalt sedimentary beds, and improves the salt–sediment terminations in the P‐wave image. Prestack and poststack depth‐migrated examples illustrate that the P‐wave and the three modes of converted waves preferentially image different parts of the base of salt. In some areas, the P‐wave BOS can be very weak, obscured by noise, or completely absent. Converted‐wave imaging complements P‐wave imaging in delineating the BOS for velocity model building.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1202-1218 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin

Reflection moveout recorded over an azimuthally anisotropic medium (e.g., caused by vertical or dipping fractures) varies with the azimuth of the source‐receiver line. Normal‐moveout (NMO) velocity, responsible for the reflection traveltimes on conventional‐length spreads, forms an elliptical curve in the horizontal plane. While this result remains valid in the presence of arbitrary anisotropy and heterogeneity, the inversion of the NMO ellipse for the medium parameters has been discussed so far only for horizontally homogeneous models above a horizontal or dipping reflector. Here, we develop an analytic moveout correction for weak lateral velocity variation in horizontally layered azimuthally anisotropic media. The correction term is proportional to the curvature of the zero‐offset traveltime surface at the common midpoint and, therefore, can be estimated from surface seismic data. After the influence of lateral velocity variation on the effective NMO ellipses has been stripped, the generalized Dix equation can be used to compute the interval ellipses and evaluate the magnitude of azimuthal anisotropy (measured by P-wave NMO velocity) within the layer of interest. This methodology was applied to a 3-D “wide‐azimuth” data set acquired over a fractured reservoir in the Powder River Basin, Wyoming. The processing sequence included 3-D semblance analysis (based on the elliptical NMO equation) for a grid of common‐midpoint “supergathers,” spatial smoothing of the effective NMO ellipses and zero‐offset traveltimes, correction for lateral velocity variation, and generalized Dix differentiation. Our estimates of depth‐varying fracture trends in the survey area, based on the interval P-wave NMO ellipses, are in good agreement with the results of outcrop and borehole measurements and the rotational analysis of four‐ component S-wave data.


2014 ◽  
Vol 54 (2) ◽  
pp. 504
Author(s):  
Sanjeev Rajput ◽  
Michael Ring

For the past two decades, most of the shear-wave (S-wave) or converted wave (P-S) acquisitions were performed with P-wave source by making the use of downgoing P-waves converting to upgoing S-waves at the mode conversion boundaries. The processing of converted waves requires studying asymmetric reflection at the conversion point, difference in geometries and conditions of source and receiver, and the partitioning of energy into orthogonally polarised components. Interpretation of P-S sections incorporates the identification of P-S waves, full waveform modeling, correlation with P-wave sections and depth migration. The main applications of P-S wave imaging are to obtain a measure of subsurface S-wave properties relating to rock type and fluid saturation (in addition to the P-wave values), imaging through gas clouds and shale diapers, and imaging interfaces with low P-wave contrast but significant S-wave changes. This study examines the major differences in processing of P and P-S wave surveys and the feasibility of identifying converted mode reflections by P-wave sources in anisotropic media. Two-dimensional synthetic seismograms for a realistic rocky mountain foothills model were studied. A Kirchhoff-based technique that includes anisotropic velocities is used for depth migration of converted waves. The results from depth imaging show that P-S section help in distinguishing amplitude associated with hydrocarbons from those caused by localised stratigraphic changes. In addition, the full waveform elastic modeling is useful in finding an appropriate balance between capturing high-quality P-wave data and P-S data challenges in a survey.


Geophysics ◽  
1988 ◽  
Vol 53 (1) ◽  
pp. 32-49 ◽  
Author(s):  
John A. Dickinson

The theoretically correct way to perform a three‐dimensional (3-D) migration of seismic data requires large amounts of data manipulation on the computer. In order to alleviate this problem, a true, one‐pass 3-D migration is commonly replaced with an approximate technique in which a series of two‐dimensional (2-D) migrations is performed in orthogonal directions. This two‐pass algorithm produces the correct answer when the velocity is constant, both horizontally and vertically. Here I analyze the error due to this algorithm when the velocities vary vertically. The analysis has two parts: first, a theoretical analysis is performed in which a formula for the error is derived; and second, a field data comparison between one‐pass and two‐pass migrations is shown. My conclusion is that two‐pass 3-D migration is, in general, a very good approximation. Its errors are usually small, the exceptions being when both the reflector dip is large (in practice this typically means greater than about 25 to 40 degrees) and the orientation of the reflector is in neither the inline nor the crossline direction. Even then the error is the same order of magnitude as that due to the uncertainty in the migration velocities. These conclusions are still valid when there is lateral velocity variation, as long as this variation is accounted for by trace stretching. The analysis presented here deals with time migration; no claims are made regarding depth migration.


Geophysics ◽  
1994 ◽  
Vol 59 (9) ◽  
pp. 1419-1434 ◽  
Author(s):  
James L. Black ◽  
Matthew A. Brzostowski

Even if the correct velocity is used, time migration mispositions events whenever the velocity changes laterally. These errors increase with lateral velocity variation, depth of burial, and dip angle θ. Our analyses of two model types, one with an implicit gradient and one with an explicit gradient, yield simple “rules of thumb” for these errors to first order in the lateral gradient. The x error is [Formula: see text], and the z error is [Formula: see text], where the quantity A = A(x, z) contains the information about depth of burial and magnitude of lateral gradient. These rules can be used to determine when depth migration is needed. Further analysis also shows that the image‐ray correction to time migration is accurate only at small dip. For dipping events, the image‐ray correction must be supplemented by a shift in x of the form [Formula: see text] and a shift in z given by [Formula: see text]. These time‐migration corrections take the same form for both the models we have studied, suggesting a general scheme for correcting time migration, which we call “remedial migration.”


2021 ◽  
pp. 1-52
Author(s):  
Youfang Liu ◽  
James Simmons

Several P-wave azimuthal anisotropy studies have been conducted for the SEAM II Barrett model data. However, these analyses provide fracture property estimation that is inconsistent with the actual model properties. Therefore, we perform a feasibility study to understand the influence of the overburden and reservoir properties, and the processing and inversion steps, which together determine the success of the fracture interpretation from seismic data. 1D model properties (orthorhombic for both overburden and reservoir) are first extracted from the actual Barrett model properties at two locations. Anisotropic prestack reflectivity modeling exposes the true orthorhombic response of the 1D medium in the form of Common Offset and Common Azimuth (COCA) gathers. The true anisotropic response is obscured in the Barrett data (generated by finite element modeling) due to the mild lateral velocity variations and orthorhombic anisotropy in the overburden. We then expose the reservoir anisotropic response by using an isotropic overburden in the reflectivity modeling. This shows that the P-wave VVAZ responses generated by the reservoir itself are weak, which leads to an unstable VVAZ inversion to estimate the interval NMO velocity anisotropy. The reservoir thickness (125m or 65ms TWT) or NMO velocity anisotropy (6-7%) needs to be at least doubled to obtain a stable VVAZ inversion. Anisotropic geometrical-spreading correction improves the amplitude-versus-azimuth (AVAZ) inversion results when reflectivity modeling models orthorhombic overburden. The converted wave ( C-wave) has a stronger VVAZ response compared to the P-wave. We suggest that the C-wave data could be useful to constrain fracture interpretation in the Barrett model. We conclude that the results of previous studies are due to the combination of the residual influence of overburden after processing and imaging, and the weak anisotropy responses from the reservoir.


Sign in / Sign up

Export Citation Format

Share Document