Time-lapse velocity analysis — Application to onshore continuous reservoir monitoring

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. B105-B117 ◽  
Author(s):  
Julien Cotton ◽  
Hervé Chauris ◽  
Eric Forgues ◽  
Paul Hardouin

In 4D seismic, the velocity model used for imaging and reservoir characterization can change as production from the reservoir progresses. This is particularly true for heavy oil reservoirs stimulated by steam injection. In the context of sparse and low-fold seismic acquisitions, conventional migration velocity analyses can be inadequate because of a poorly and irregularly sampled offset dimension. We update the velocity model in the context of daily acquisitions with buried sources and receivers. The main objective is to demonstrate that subtle time-lapse effects can be detected over the calendar time on onshore sparse acquisitions. We develop a modified version of the conventional prestack time migration to detect velocity changes obtained after crosscorrelation of the base and monitor surveys. This technique is applied on a heavy oil real data set from the Netherlands and reveals how the steam diffuses over time within the reservoir.

Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. S105-S111 ◽  
Author(s):  
Sheng Xu ◽  
Feng Chen ◽  
Bing Tang ◽  
Gilles Lambare

When using seismic data to image complex structures, the reverse time migration (RTM) algorithm generally provides the best results when the velocity model is accurate. With an inexact model, moveouts appear in common image gathers (CIGs), which are either in the surface offset domain or in subsurface angle domain; thus, the stacked image is not well focused. In extended image gathers, the strongest energy of a seismic event may occur at non-zero-lag in time-shift or offset-shift gathers. Based on the operation of RTM images produced by the time-shift imaging condition, the non-zero-lag time-shift images exhibit a spatial shift; we propose an approach to correct them by a second pass of migration similar to zero-offset depth migration; the proposed approach is based on the local poststack depth migration assumption. After the proposed second-pass migration, the time-shift CIGs appear to be flat and can be stacked. The stack enhances the energy of seismic events that are defocused at zero time lag due to the inaccuracy of the model, even though the new focused events stay at the previous positions, which might deviate from the true positions of seismic reflection. With the stack, our proposed approach is also able to attenuate the long-wavelength RTM artifacts. In the case of tilted transverse isotropic migration, we propose a scheme to defocus the coherent noise, such as migration artifacts from residual multiples, by applying the original migration velocity model along the symmetry axis but with different anisotropic parameters in the second pass of migration. We demonstrate that our approach is effective to attenuate the coherent noise at subsalt area with two synthetic data sets and one real data set from the Gulf of Mexico.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 631-650 ◽  
Author(s):  
Mark E. Mathisen ◽  
Anthony A. Vasiliou ◽  
Paul Cunningham ◽  
J. Shaw ◽  
J. H. Justice ◽  
...  

Time‐lapse crosswell seismic data acquired with a cemented receiver cable have been processed into P‐ and S‐wave tomograms which image heavy oil sand lithofacies and changes as a result of steam injection. Twenty‐seven crosswell surveys were acquired between two wells over a 3.5 month period before, during, and after a 34‐day, 30 MBBL [Formula: see text] steam injection cycle. Interpretation was based on correlations with reservoir data and models, observation well data, and engineering documentation of the production history and steam cycle. Baseline S‐ and P‐wave tomograms image reservoir sand flow units and areas affected by past cyclic steam injection. S‐wave tomograms define lithology and porosity contrasts between the excellent reservoir quality, “high flow” turbidite channel facies and the interbedded “low to moderate flow” bioturbated levee facies. The reservoir dip of approximately 20° is defined by the velocity contrast between lithofacies. P‐wave baseline tomograms image lithology, porosity, structure, and several low velocity zones caused by past steam injection. Previous steam‐heat injection caused the formation of gas which reduced velocities as much as several thousand ft/s (600 m/s), an amount which obscures the velocity contrast between lithofacies and smaller velocity reductions as a result of temperature alone. Time‐lapse and difference P‐wave tomograms document several areas with small decreases in velocity during steam injection and larger decreases after cyclic steam injection. Velocity reductions range from 300 to 900 ft/s (90 to 270 m/s) adjacent to and above injectors located 20 to 50 feet (6 to 15 m) from the tomogram cross‐section. Poisson’s ratio tomograms show a significant decrease (.10) in the same area, and include low values indicative of gas saturation. Continuous injectors located 50 to 350 feet (15 to 100 m) from the survey area also caused a progressive decrease in velocity of the “high flow” channel sands during the time‐lapse survey. Interdisciplinary interpretation indicates that tomograms not only complement other borehole‐derived reservoir characterization and temperature monitoring data but can be used to quantitatively characterize interwell reservoir properties and monitor changes as a result of the thermal recovery process. Monitoring results over 3.5 months confirms that stratification has controlled the flow of steam, in contrast to gravity override. This suggests that tomographic images of reservoir flow‐units and gas‐bearing high temperature zones should be useful for positioning wells and optimizing injection intervals, steam volumes, and producing well completions.


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. U29-U38 ◽  
Author(s):  
Andreas Kjelsrud Evensen ◽  
Martin Landrø

Most seismic studies of changes in traveltimes are of a qualitative nature and a major challenge in four dimensions is to use the information contained in time shifts to quantify the nature of velocity changes in the subsurface layers. We propose a 4D tomographic inversion method that uses time shifts from prestack seismic data to estimate parameters describing the 2D velocity field after changes have occurred. Prestack data allow for the usage of many offsets, thus increasing the information input for the inversion. The velocity changes are parameterized by a chosen number of Gaussian functions in two dimensions and weighted least-squares inversion is used to estimate the parameters describing these functions. We have found that the parameters describing the position and shape of the Gaussian velocity anomalies can be estimated with this method for simple synthetic cases. For more complex cases with overlapping Gaussian functions, resolution of the parameters can be difficult and in these cases our recommendation is to find the best fit for a simple smooth anomaly to a more complex real world. The method is tested on a real data set from a [Formula: see text] injection project above the Sleipner field in the North Sea, where quantification of changes is important for monitoring purposes. We have found that the noise levels in prestack traveltime data are on the high side for large-scale analysis; however, we estimate reasonable [Formula: see text] layer thickness and velocity compared to previous work in a nearby area.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. C141-C152
Author(s):  
Nicolas Mastio ◽  
Pierre Thore ◽  
Marianne Conin ◽  
Guillaume Caumon

In the petroleum industry, time-lapse (4D) studies are commonly used for reservoir monitoring, but they are also useful to perform risk assessment for potential overburden deformations (e.g., well shearing, cap-rock integrity). Although complex anisotropic velocity changes are predicted in the overburden by geomechanical studies, conventional time-lapse inversion workflows only deal with vertical velocity changes. To retrieve the geomechanically induced anisotropy, we have adopted a reflection traveltime tomography method coupled with a time-shift estimation algorithm of prestack data of the baseline and monitor simultaneously. For the 2D approach, we parameterize the anisotropy using five coefficients, enough to cover any type of anisotropy. Before applying the workflow to a real data set, we first study a synthetic data set based on the real data set and include velocity variations between baseline and monitor found in the literature (vertical P-wave velocity decrease in the cap rock and isotropic P-wave velocity change in the reservoir). On the synthetics, we measure the angular ray coverage necessary to retrieve the target anisotropy and observe that the retrieved anisotropies depend on the offset range. Based on a synthetic experiment, we believe that the acquisition of the real case study is suitable for performing tomographic inversion. The anisotropic velocity changes obtained on three inlines separated by 375 m are consistent and show a strong positive anomaly in the cap rock along the 45° direction (the [Formula: see text] parameter in Thomsen notation), whereas the vertical velocity change is surprisingly almost negligible. We adopt a rock-physics explanation compatible with these observations and geologic considerations: a reactivation of water-filled subvertical cracks.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.


Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.


Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.


Author(s):  
S. Yuh ◽  
M. Le Ravalec ◽  
C. Hubans ◽  
P.-O. Lys ◽  
D. Foulon

Sign in / Sign up

Export Citation Format

Share Document