Moveout Correction Research Articles

Goos–Hänchen induced normal moveout correction for wide-angle reflections

The Goos–Hanchen (GH) lateral shift has been theoretically simulated and observed in lab. GH lateral shift introduces additional traveltime and distance when the incidence angles are larger than the critical angle. For seismic wave, this GH shift is caused by the total reflection of an incident beam of P-wave from low to high impedance medium at near and post-critical angles. Because of its large influences on traveltime and lateral shift displacement, the GH shift should be corrected in normal moveout (NMO) correction for wide-angle reflections in seismic data processing. In this paper, we derive the partial derivatives of reflection coefficients (PP- and PSV-wave) with respect to circular frequency using the Zoeppritz equations. Then, the delay time and NMO correction term with the behavior of GH lateral shift is derived. The characteristics of delay time and GH induced time differences are analyzed. The results show that this GH shift could be either positive or negative and the delay on time has large influences on seismic reflections when the incidence angles are larger than the critical angles. The efficiency of GH induced NMO correction is tested using synthetic seismic data. The GH induced NMO correction should be done for wide-angle reflections during the progress of seismic data processing.

PRE STACK DEPTH MIGRATION UNTUK KOREKSI EFEK PULL UP DENGAN MENGGUNAKAN METODE HORIZON BASED DEPTH TOMOGRAPHY PADA LAPANGAN ‘A1 DAN A2’

In the case of seismic data processing with sandstone lithology such as shale and carbonate often get the result of data processing which have pull up effect especially on the time domain migration result. Pre stack depth migration is a processing based on focusing the amplitude according to the actual depth by using the input interval velocity. Migration is performed using kirchhoff pre stack depth migration algorithm. Pre stack depth migration is done with modeling of horizontal based depth tomography method. This method uses residual moveout correction applied along the horizon-picking line. This research uses two field data that is A1 and A2 Field. A1field has characteristics of carbonate rock that produce pull up shaped similar to carbonate layer. A2 field has a pull-up effect that is not very clear but has build up because of the layer above it. Stages performed starting from the processing of pre stack time migration in the form of velocity picking, generate rms velocity and migration time domain. The pre stack depth migration process begins with a velocity transformation with the dix transformation equation to generate interval velocity, migrate Pre stack depth migration, perform horizon interpretations and perform velocity modeling using the horizon based depth tomography method. The iteration is done 4 times and resulted in the final section of pre stack depth migration which has been corrected by pull up effect.

Upper mantle and mantle transition zone thermal and water content anomalies beneath NE Asia: Constraints from receiver function imaging of the 410 and 660 km discontinuities

The 410 and 660 km discontinuities (d410 and d660) bordering the mantle transition zone (MTZ) beneath NE Asia, including NE China, Eastern Mongolia, and southern Siberia, are imaged in successive circular bins with a radius of 1 degree by stacking a total of 274,413 P-to-s radial receiver functions recorded by 799 broadband seismic stations. After moveout correction based on the 1-D IASP91 Earth model, the resulting apparent depths of the discontinuities exhibit significant and spatially systematic variations. Three approximately N-S elongated narrow zones with significantly thickened MTZ are observed, which may be associated with the thermal effect and dehydration of subducted slabs. The major volcanoes in NE China are underlain by a d660 that is apparently depressed by ∼19 km, which can be interpreted by the presence of an anomalously high water concentration in the lower MTZ released from the stagnated slabs. Low wavespeed anomalies above the d410 west of the Datong volcanic fields are underlain by an MTZ with normal thickness, and are attributable by dehydration of the leading portion of the stagnant Pacific slab in the MTZ that is revealed in an N-S oriented narrow zone east of this area. The lateral shift of the upper mantle low wavespeed zone and the area with thickened MTZ may suggest a westward drift of the upper mantle relative to the subducted slab. An abnormally thin MTZ is observed beneath the Hangay Dome in central Mongolia, suggesting the possible existence of thermal upwelling from the lower mantle through the MTZ.

Stacking Seismic Data Based on Ramanujan Sums

Real seismic data consist of multiple traces captured by an array of receivers. These multiple traces are sorted by the common midpoint between the source and the receiver, and the time-lag between different traces is corrected by a normal move-out correction process. After these preprocessing steps, the sorted traces contain the same information about the earth sublayers and are corrupted by noise. The next step, termed stacking in seismic signal processing, involves the construction of an optimum trace with an improved signal-to-noise ratio (SNR) from these sorted traces. In this letter, we present an improved method for weighted stacking, where each trace is weighed in accordance with the noise variance. Using the first-order derivative property of Ramanujan sums, we perform the estimation of noise variance in each trace. Numerical results demonstrate that the method presented in this letter has better SNR for the trace obtained after stacking in comparison with existing methods.

Estimating normal moveout velocity using the recurrent neural network

Machine learning (ML) has recently gained immense popularity because of its successful application in complex problems. It develops an abstract relation between the input and output. We have evaluated the application of ML to the most basic seismic processing of normal moveout (NMO) correction. The arrival times of reflection events in a common midpoint (CMP) gather follow a hyperbolic trajectory; thus, they require a correction term to flatten the CMP gather before stacking. This correction term depends on an rms velocity, also referred to as the NMO velocity. In general, NMO velocity is estimated using the semblance measures and picking the peaks in the velocity panel. This process requires a lot of human intervention and computation time. We have developed a novel method using one of the tools based on an ML- approach and applied to the NMO velocity estimation problem. We use the recurrent neural network (RNN) to estimate the NMO velocity directly from the seismic data. The input to the network is a seismic gather and corresponding precalculated NMO velocity (as prelabeled data set) to flatten the gather. We first train the network to develop a relationship between the input gathers (before NMO correction) and the corresponding NMO velocities for a few CMPs as a supervised learning process. Adam optimization algorithm is used to train the RNN. The output from the network is then compared against the correct NMO velocity. The error between the two velocities is then used to update the weight of the neurons and to minimize the mean-squared error between the two velocities. After the network is trained, it can be used to calculate the NMO velocity for the rest of the seismic gathers. We evaluate our method on a noisy data set from Poland. We used only 10% of the CMPs to train the network, and then we used the trained network to predict NMO velocity for the remaining CMP locations. The stack section obtained by using RNN-generated NMO velocities is nearly identical to that obtained by the conventional semblance method.

Automated static and moveout corrections of high‐resolution seismic data from the Baltic Sea

ABSTRACTHigh‐frequency multichannel seismic systems provide detailed images of the shallow marine subsurface. In order to exploit the redundancy inherent in such data optimally, traveltime corrections need to account for normal moveout and static effects due to vertical source and receiver variations. Misalignment of reflections in common‐midpoint gathers will significantly lower the frequency content in the final stack, making this correction particularly important for very high‐frequency seismic data. Traditionally, normal moveout correction involves labour‐intensive picking of stacking velocities, while static corrections can be, by some techniques, performed automatically. In this paper, we present a high‐frequency seismic case study from the Baltic Sea, using seismic image matching as a novel, fully automated technique to perform joint moveout and static corrections. Our multichannel test profiles were acquired offshore Rügen island for wind farm development. Owing to the regular passage of up to 1.5 m high ocean waves during data acquisition, these boomer profiles suffer from strong static effects. We perform joint normal moveout and static corrections by defining the nearest common offset section as a fixed reference frame and minimizing its difference in traveltime with respect to all available common offset sections. Time shifts are computed independent of a pre‐defined traveltime curve, using the normalized cross‐correlation as a measure of data similarity while penalizing irregular displacements by a regularization term. Time shifts are converted to stacking velocities based on the traditional hyperbolic traveltime equation. Our results are compared with those derived by conventional manual velocity analysis and subsequent trim static corrections. We find that image matching produces stacks of similar quality and stacking velocity models of similar to slightly better quality compared with the conventionally derived ones, revealing the potential of this technique to automatize and significantly speed up this first part of the seismic processing chain.

A field experiment of walkaway distributed acoustic sensing vertical seismic profile in a deep and deviated onshore well in Japan using a fibre optic cable deployed inside coiled tubing

ABSTRACTA two‐dimensional walkaway vertical seismic profiling survey using distributed acoustic sensing was conducted at an onshore site in Japan. The maximum depth and the deviation of the observation well were more than 4,000 m and 81 degrees, respectively. Among the several methods for installing fibre optic cables, we adopted the inside coiled tubing method, in which coiled tubing containing a fibre optic cable is deployed. The signal‐to‐noise ratio of the raw shot gather was low, possibly due to poor coupling between the fibre optic cable and the subsurface formation resulting from the fibre optic cable deployment method and the existence of considerable tubewave noise. Nevertheless, direct P‐wave arrivals, P–P reflections and P–S converted waves exhibited acceptable signal‐to‐noise ratios after careful optimization of gauge length for distributed acoustic sensing optical processing and the application of carefully parameterized tubewave noise suppression. One of the challenges in current distributed acoustic sensing vertical seismic profile data processing is the separation of P‐ and S‐waves using only one‐component measurements. Hence, we applied moveout correction using two‐dimensional ray tracing. This process effectively highlights only reflected P‐waves, which are used in subsequent subsurface imaging. Comparison with synthetic well seismograms and two‐dimensional surface seismic data confirms that the final imaging result has a sufficiently high quality for subsurface monitoring. We acquired distributed acoustic sensing vertical seismic profile data under both flowing conditions and closed conditions, in which the well was shut off and no fluid flow was allowed. The two imaging results are comparable and suggest the possibility of subsurface imaging and time‐lapse monitoring using data acquired under flowing conditions. The results of this study suggest that, by adopting the inside coiled tubing method without drilling a new observation well, more affordable distributed acoustic sensing vertical seismic profile monitoring can be achieved in fields such as CO2 capture and storage and unconventional shale projects, where monitoring costs have to be minimized.

Spherical multifocusing method for irregular topography

Moveout correction for irregular topography has been a longstanding challenge in processing seismic exploration data. Irregular topography usually results in large moveout among traces, a low signal-to-noise ratio (S/N), and difficulty in modeling near-surface velocities. Conventional normal moveout (NMO) corrections and elevation static methods are imprecise and tend to introduce significant errors for large offsets. Over the past two decades, several multiparameter time corrections and stacking techniques to reduce noise and improve resolution have been proposed in place of the classic NMO and common-midpoint stack. These include the common-reflection-surface (CRS), common-offset CRS, nonhyperbolic CRS, implicit CRS, multifocusing (MF), irregular surface MF (IS-MF), spherical MF (SMF), and common-offset MF methods. Various CRS-type operators that consider the top-surface topography have been proposed. For MF-type operators, only IS-MF can be applied directly to the irregular topography with no elevation statics required. In this study, we have developed a new MF formulation, modifying the SMF method to consider nonzero elevations of sources and receivers and we corrected moveout of nonplanar data directly without prior elevation static corrections. The proposed extension combines the sensitivity to spherical reflectors of SMF with the applicability of the IS-MF method to irregular topography. We investigated the behavior of the new operator using a physical model data set and compared the results with those from the conventional IS-MF method. The results revealed that the new operator is more robust over a wide range of source and receiver elevations and has advantages on strongly curved interfaces. We also confirmed the potential of the proposed approach by comparing stacking results for a real-land data set with a low S/N.

Goal-Oriented Inversion-Based NMO Correction Using a Convex $l_{2,1}$ -Norm

Normal moveout (NMO) correction is a routine step in seismic data processing, which has an important impact on other seismic processing procedures, seismic inversion, and interpretation. We propose a goal-oriented inversion-based NMO correction method using a convex $l_{2,1}$ -norm. The proposed method corrects the superresolution multichannel offset-dependent reflectivity rather than the bandlimited data itself sample by sample, block by block, or wavelet by wavelet. Therefore, the proposed method can essentially reduce the amplitude and even phase distortion introduced by data-based NMO correction methods in the presence of strong wavelet interference. We impose two goal-oriented constraints including both the temporal sparsity and the horizontal continuity of reflectivity, which are approximately represented by a convex $l_{2,1}$ -norm, on the geometric moveout relationship from offset-dependent trajectories to zero offsets to build a new objective function for NMO correction. The goal-desired temporal sparsity of reflectivity can induce the superresolution solution; meanwhile, the goal-desired horizontal continuity introduces a reasonable intrinsic structure to further limit the solution space and is particularly suitable to processing interfering reflections. Attributing to these two additional constraints, the new NMO correction method can flatten the interfering events and the intersecting events with favorable offset-dependent amplitude and phase variations even in the presence of noise. Synthetic and real data examples are adopted to verify the performance of our method. The results show that goal-oriented inversion-based NMO correction using the $l_{2,1}$ -norm is a potentially effective, stable, and high-quality NMO correction tool, especially for strong wavelet interference and at far offsets.

Three-parameter normal moveout correction in layered anisotropic media: A stretch-free approach

In conventional normal moveout (NMO) correction, some parts of the recorded data at larger offsets are discarded because of NMO distortions. Deviation from the true traveltime of reflections due to the anisotropy and heterogeneity of the earth, and wavelet stretching are two reasons of these distortions. The magnitudes of both problems increase with increasing the offset to depth ratio. Therefore, to be able to keep larger offsets of shallower reflections, both problems should be obviated. Accordingly, first, we have studied different traveltime approximations being in use, alongside new parameterizations for two classical functional equations, to select suitable equations for NMO correction. We numerically quantify the fitting accuracy and uncertainty of known nonhyperbolic traveltime approximations for P-waves in transversely isotropic media with vertical symmetry axis (VTI). We select three suitable three-parameter approximations for NMO in layered VTI media as the VTI generalized moveout approximation, a double-square-root approximation, and a perturbation-based approximation. Second, we have developed an extension of the earlier proposed stretch-free NMO method, using the selected moveout approximations. This method involves an automatic modification of the input parameters in anisotropic NMO correction, for selected reflections. Our anisotropic stretch-free NMO method is tested on synthetic and three real data sets from Gulf of Mexico and Iranian oil fields. The results verify the success of the method in extending the usable offsets, by generating flat and stretch-free NMO corrected reflections.

Automatic migration velocity estimation for prestack time migration

Migration velocity analysis is a labor-intensive part of the iterative prestack time migration (PSTM) process. We have developed a velocity estimation scheme to improve the efficiency of the velocity analysis process using an automatic approach. Our scheme is the numerical implementation of the conventional velocity analysis process based on residual moveout analysis. The key aspect of this scheme is the automatic event picking in the common-reflection point (CRP) gathers, which is implemented by semblance scanning trace by trace. With the picked traveltime curves, we estimate the velocities at discrete grids in the velocity model using the least-squares method, and build the final root-mean-square (rms) velocity model by spatial interpolation. The main advantage of our method is that it can generate an appropriate rms velocity model for PSTM in just a few iterations without manual manipulations. In addition, using the fitting curves of the picked events in a range of offsets to estimate the velocity model, which is fitting to a normal moveout correction, can prevent our scheme from the local minima issue. The Sigsbee2B model and a field data set are used to verify the feasibility of our scheme. High-quality velocity model and imaging results are obtained. Compared with the computational cost to generate the CRP gathers, the cost of our scheme can be neglected, and the quality of the initial velocity is not critical.

Reflection‐constrained 2D and 3D non‐hyperbolic moveout analysis using particle swarm optimization

ABSTRACTThe objective of moveout parameter inversion is to derive sets of parameter models that can be used for moveout correction and stacking at each common midpoint location to increase the signal‐to‐noise ratio of the data and to provide insights into the kinematic characteristics of the data amongst other things. In this paper, we introduce a data‐driven user‐constrained optimization scheme that utilizes manual picks at a point on each reflector within a common midpoint gather to constrain the search space in which an optimization procedure can search for the optimal parameter sets at each reflection. The picks are used to create boundary curves which can be derived approximately via an optimization technique or analytically via the derivation of an analytical bounds function. In this paper, we derive analytical forms of bounds functions for four different moveout cases. These are normal moveout, non‐hyperbolic moveout, azimuthally dependent normal moveout and azimuthally dependent non‐hyperbolic moveout. The optimization procedure utilized here to search for the optimal moveout parameters is the particle swarm optimization technique. However, any metaheuristic optimization procedure could be modified to account for the constraints introduced in this paper. The technique is tested on two‐layer synthetic models based on three of the four moveout cases discussed in this paper. It is also applied to an elastic forward modelled synthetic model called the HESS model, and finally to real 2D land data from Alaska. The resultant stacks show a marked improvement in the signal‐to‐noise ratio compared to the raw stacks. The results for the normal moveout, non‐hyperbolic moveout and azimuthally dependent normal moveout tests suggest that the method is viable for said models. Results demonstrate that our method offers potential as an alternative to conventional parameter picking and inversion schemes, particularly for some cases where the number of parameters in the moveout approximation is 2 or greater.

Superresolution stacking based on sparse Radon transform

Normal moveout (NMO) correction is routinely applied to each common-midpoint (CMP) gather before forming a stack section. Trace stretching is inherent in the conventional NMO correction, causing frequency distortion and temporal resolution degradation in the stacked image. We have developed a novel method to compute a superresolution stack directly from CMP gathers by solving an optimization problem that leads to an [Formula: see text]-norm minimization of the model based on compressive sensing principles. Our algorithm is able to recover the superposition of signals in a specialized parallel-shifted hyperbolic Radon-transform domain. The superresolution stacking result is then extracted directly from the transform domain without inversely transforming back to the time-offset domain. The effectiveness of the developed methodology in producing superresolution stack images was tested on a synthetic data set and a land data set from central Saudi Arabia. Test results clearly indicate that our method reduces the stretch effect, enhances the signal, attenuates multiples and noise, and results in a substantial improvement on the temporal and spatial resolutions of the stack image especially in regions dominated by noise.

Reprocessing legacy three-dimensional seismic data from the Halfmile Lake and Brunswick No. 6 volcanogenic massive sulphide deposits, New Brunswick, Canada

We reprocessed legacy three-dimensional (3D) seismic data from the Halfmile Lake and Brunswick areas, both of which were acquired for mineral exploration in the Bathurst Mining Camp, New Brunswick. Each 3D seismic survey was acquired over known volcanogenic massive sulphide deposits and covered areas with strong mineral potential. Most improvements resulted from a reduction of coherent and random noise on prestack gathers and from an improved velocity model, combined with re-imaging with dip moveout corrections and poststack migration or prestack time migration. At Halfmile Lake, the new imaging results show the Deep zone and a possible extension of the sulphide mineralization at greater depth. True amplitude processing has shown that this anomaly has strong amplitudes and is offset from the Deep zone by a shallowly dipping fault (<15°). With the clearer geological context provided by our results, this anomaly, which appears as a stand-alone anomaly on an original image obtained by Noranda Exploration Ltd., becomes a defendable exploration target. Nonorthogonal acquisition geometry and receiver patches of the Brunswick No. 6 3D seismic survey generated artefacts after dip moveout processing that reduced the overall quality of the seismic volumes. By using a filtering approach based on the application of a weighted Laplacian-Gaussian filter in the Kx–Ky domain, we reduced the noise and improved the continuity of reflections. We also imaged the short and flat reflections observed previously only in the shallow part of prestack time migrated data. These short reflections appear as diffractions on the filtered stacked section with dip moveout corrections, indicating that they originate from small geological bodies or discontinuities in the subsurface.

Stretch‐free generalized normal moveout correction

ABSTRACTThe effective application of normal moveout correction processes mainly depends on four factors: the chosen traveltime approximation, the stretching associated with the given traveltime, crossing events and phase changes, the last two being inherent to the seismic data. In this context, we conduct a quantitative analysis on stretching considering a general traveltime expression depending on half‐offset and midpoint coordinates. Through this analysis, we propose a mathematically proven procedure to eliminate stretching, which can be applied to any traveltime approximation. The proposed method is applied to synthetic and real data sets, considering different traveltime approximations and achieved complete elimination of stretching.

Absence of thermal influence from the African Superswell and cratonic keels on the mantle transition zone beneath southern Africa: Evidence from receiver function imaging

The depths of the 410 km (d410) and 660 km (d660) discontinuities beneath southern Africa, which is presumably underlain by the lower-mantle African Superswell, are imaged in 1∘ radius consecutive circular bins using over 6400 P-to-S receiver functions (RFs) recorded by 130 seismic stations over a 27 yr period. When the IASP91 standard Earth model is utilized for moveout correction and time-depth conversion, a normal mantle transition zone (MTZ) thickness of 246 ± 7 km is observed, suggesting that the Superswell has no discernible effect on mantle transition zone temperature. Based on the negligible disparity of the mean MTZ thicknesses between on (246 ± 6 km) and off (246 ± 8 km) cratonic regions, we conclude that the deep Archean cratonic keels possess limited influence on MTZ thermal structure. The apparently shallower-than-normal MTZ discontinuities and the parallelism between the d410 and d660 are mostly the results of upper mantle high wave speed anomalies probably corresponding to a thick lithosphere with a mean thickness of about 245 km beneath the Kaapvaal and 215 km beneath the Zimbabwe cratons. In contradiction to conclusions from some of the previous studies, the resulting spatial distribution of the stacking amplitudes of the P-to-S converted phases at the discontinuities is inconsistent with the presence of an excessive amount of water in the MTZ and atop the d410.

Seismic Data Processing of Subba Oil Field in South Iraq

Evaluation study was conducted for seismic interpretation using two-dimensional seismic data for Subba oil field, which is located in the southern Iraq. The Subba oil field was discovered in 1973 through the results of the seismic surveys and the digging of the first exploratory well SU-1 in 1975 to the south of the Subba oil field. The entire length of the field is 35 km and its width is about 10 km. The Subba oil field contains 15 wells most of them distributed in the central of the field. This study is dealing with the field data and how to process it for the purpose of interpretation; the processes included conversion of field data format, compensation of lost data and noise disposal, as well as the application of normal moveout correction for the purpose of improving the quality of seismic sections, using Omiga software at the processing department/ Petroleum Exploration Company After processing RMS velocity map was done also constructed from the processing data and showed a low velocity compatibility north of the study area and an increase in velocity to the south of the oil field

Topography of the Mantle Transition Zone Discontinuities Beneath Alaska and Its Geodynamic Implications: Constraints From Receiver Function Stacking

AbstractThe 410 and 660 km discontinuities (d410 and d660, respectively) beneath Alaska and adjacent areas are imaged by stacking 75,296 radial receiver functions recorded by 438 broadband seismic stations with up to 30 years of recording period. When the 1‐D IASP91 Earth model is used for moveout correction and time depth conversion, significant and spatially systematic variations in the apparent depths of the d410 and d660 are observed. The mean apparent depth of the d410 and d660 for the entire study area is 417 ± 12 km and 665 ± 12 km, respectively, and the mean mantle transition zone (MTZ) thickness is 248 ± 8 km which is statistically identical to the global average. For most of the areas, the undulations of the apparent depths of the d410 and d660 are highly correlated, indicating that lateral velocity variations in the upper mantle above the d410 contribute to the bulk of the observed apparent depth variations by affecting the traveltimes of the P‐to‐S converted phases from both discontinuities. Beneath central Alaska, a broad zone with greater than normal MTZ thicknesses and shallower than normal d410 is imaged, implying that the subducting Pacific slab has reached the MTZ and is fragmented or significantly thickened. Within the proposed Northern Cordilleran slab window, an overall thinner than normal MTZ is observed and is most likely the result of a depressed d410. This observation, when combined with results from seismic tomography investigations, may indicate advective thermal upwelling from the upper MTZ through the slab window.

Nonstretching normal-moveout correction using a dynamic time warping algorithm

Normal-moveout (NMO) correction is one of the most important routines in seismic processing. NMO is usually implemented by a sample-by-sample procedure; unfortunately, such implementation not only decreases the frequency content but also distorts the amplitude of seismic waveforms resulting from the well-known stretch. The degree of stretch increases with increasing offset. To minimize severe stretch associated with far offset, we use a dynamic time warping (DTW) algorithm to achieve an automatic dynamic matching NMO nonstretch correction, which does not handle crossing events and convoluted events such as thin layers. Our algorithm minimizes the stretch through an automatic static temporal correction of seismic wavelets. The local static time shifts are obtained using a DTW algorithm, which is a nonlinear optimization method. To mitigate the influence of noise, we evaluated a multitrace window strategy to improve the signal-to-noise ratio of seismic data by obtaining a more precise moveout correction at far-offset traces. To illustrate the effectiveness of our algorithm, we first applied our method to synthetic data and then to field seismic data. Both tests illustrate that our algorithm minimizes the stretch associated with far offsets, and the method preserves the amplitude fidelity.

Oriented time-domain dip moveout correction for planar reflectors in common-source domain

Oriented time-domain imaging can be orders of magnitude faster than the routine techniques, which rely on velocity analysis. The term “oriented” refers to those techniques that use the information carried by local slopes. Time-domain dip moveout (DMO) correction, which has often been ignored by the seismic imaging community, has been coming back to attention within the last few years. I have developed an oriented time-domain DMO correction workflow that does not face the problematic loop between the dip-dependent and/or dip-independent velocities existing in the classic DMO correction algorithms. Use of the proposed approach is also advantageous over the previous oriented techniques; the proposed technique is independent of the wavefront curvature, and the input seismic data do not need to be sorted in two different domains. The application of the technique is limited to reflectors with a small curvature. The theory of the proposed technique is investigated on a simple synthetic data example and then applied to a 2D marine data set.