Time-lapse image-domain tomography using adjoint-state methods

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. A29-A33 ◽  
Author(s):  
Jeffrey Shragge ◽  
Tongning Yang ◽  
Paul Sava
Keyword(s):  
Seismic Velocity ◽  
Time Lapse ◽  
Velocity Estimation ◽  
Data Sets ◽  
Velocity Change ◽  
Image Domain ◽  
Adjoint State ◽  
The Difference ◽  
Monitor Image ◽  
4D Seismic

Adjoint-state methods (ASMs) have proven successful for calculating the gradients of the functionals commonly found in geophysical inverse problems. The 3D ASM image-domain tomography (IDT) formulation of the seismic velocity estimation problem highlights imperfections in migrated image volumes and, using appropriate penalty functions (e.g., differential semblance), forms an objective function that can be minimized using standard optimization approaches. For time-lapse (4D) seismic scenarios, we show that the 3D ASM-IDT approach can be extended to multiple (e.g., baseline and monitor) data sets and offers high-quality estimates of subsurface velocity change. We discuss two different penalty operators that lead to absolute and relative 4D inversion strategies. The absolute approach uses the difference of two independent 3D inversions to estimate a 4D model perturbation (i.e., slowness squared). The relative approach inverts for the model perturbation that optimally matches the monitor image to the baseline image — even if migrated energy is imperfectly focused. Both approaches yield good 4D slowness estimates; however, we assert that the relative approach is more robust given the ubiquitous presence of nonrepeatable 4D acquisition noise and imperfect baseline model estimates.

Common-focus point-based target-oriented imaging approach for continuous seismic reservoir monitoring

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. M41-M48 ◽  
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali
Keyword(s):  
Seismic Data ◽  
Rock Physics ◽  
Time Lapse ◽  
Velocity Estimation ◽  
Estimation Methods ◽  
Data Sets ◽  
Data Set ◽  
Velocity Models ◽  
Reservoir Monitoring ◽  
Focus Point

The ideal approach for continuous reservoir monitoring allows generation of fast and accurate images to cope with the massive data sets acquired for such a task. Conventionally, rigorous depth-oriented velocity-estimation methods are performed to produce sufficiently accurate velocity models. Unlike the traditional way, the target-oriented imaging technology based on the common-focus point (CFP) theory can be an alternative for continuous reservoir monitoring. The solution is based on a robust data-driven iterative operator updating strategy without deriving a detailed velocity model. The same focusing operator is applied on successive 3D seismic data sets for the first time to generate efficient and accurate 4D target-oriented seismic stacked images from time-lapse field seismic data sets acquired in a [Formula: see text] injection project in Saudi Arabia. Using the focusing operator, target-oriented prestack angle domain common-image gathers (ADCIGs) could be derived to perform amplitude-versus-angle analysis. To preserve the amplitude information in the ADCIGs, an amplitude-balancing factor is applied by embedding a synthetic data set using the real acquisition geometry to remove the geometry imprint artifact. Applying the CFP-based target-oriented imaging to time-lapse data sets revealed changes at the reservoir level in the poststack and prestack time-lapse signals, which is consistent with the [Formula: see text] injection history and rock physics.

3D DAS-VSP Illumination Modeling for CO2 Plume Migration Monitoring in Offshore Sarawak, Malaysia

2021 ◽  
Author(s):  
Pankaj Kumar Tiwari ◽  
Zoann Low ◽  
Parimal Arjun Patil ◽  
Debasis Priyadarshan Das ◽  
Prasanna Chidambaram ◽  
...  
Keyword(s):  
Dynamic Simulation ◽  
Seismic Velocity ◽  
Time Lapse ◽  
Carbonate Reservoir ◽  
Fiber Optic Sensor ◽  
Co2 Injection ◽  
Storage Site ◽  
Plume Migration ◽  
Co2 Plume ◽  
4D Seismic

Abstract Monitoring of CO2 plume migration in a depleted carbonate reservoir is challenging and demand comprehensive and trailblazing monitoring technologies. 4D time-lapse seismic exhibits the migration of CO2 plume within geological storage but in the area affected by gas chimney due to poor signal-to-noise ratio (SNR), uncertainty in identifying and interpretation of CO2 plume gets exaggerated. High resolution 3D vertical seismic profile (VSP) survey using distributed acoustic sensor (DAS) technology fulfil the objective of obtaining the detailed subsurface image which include CO2 plume migration, reservoir architecture, sub-seismic faults and fracture networks as well as the caprock. Integration of quantitative geophysics and dynamic simulation with illumination modelling dignify the capabilities of 3D DAS-VSP for CO2 plume migration monitoring. The storage site has been studied in detailed and an integrated coupled dynamic simulation were performed and results were integrated with seismic forward modeling to demonstrate the CO2 plume migration with in reservoir and its impact on seismic amplitude. 3D VSP illumination modelling was carried out by integrating reservoir and overburden interpretations, acoustic logs and seismic velocity to illustrate the subsurface coverage area at top of reservoir. Several acquisition survey geometries were simulated based on different source carpet size for effective surface source contribution for subsurface illumination and results were analyzed to design the 3D VSP survey for early CO2 plume migration monitoring. The illumination simulation was integrated with dynamic simulation for fullfield CO2 plume migration monitoring with 3D DAS-VSP by incorporating Pseudo wells illumination analysis. Results of integrated coupled dynamic simulation and 4D seismic feasibility were analyzed for selection of best well location to deploy the multi fiber optic sensor system (M-FOSS) technology. Amplitude response of synthetic AVO (amplitude vs offsets) gathers at the top of carbonate reservoir were analyzed for near, mid and far angle stacks with respect to pre-production as well as pre-injection reservoir conditions. Observed promising results of distinguishable 25-30% of CO2 saturation in depleted reservoir from 4D time-lapse seismic envisage the application of 3D DAS-VSP acquisition. The source patch analysis of 3D VSP illumination modelling results indicate that a source carpet of 6km×6km would be cos-effectively sufficient to produce a maximum of approximately 2km in diameter subsurface illumination at the top of the reservoir. The Pseudo wells illumination analysis results show that current planned injection wells would probably able to monitor early CO2 injection but for the fullfield monitoring additional monitoring wells or a hybrid survey of VSP and surface seismic would be required. The integrated modeling approach ensures that 4D Seismic in subsurface CO2 plume monitoring is robust. Monitoring pressure build-ups from 3D DAS-VSP will reduce the associated risks.

Assessment of 4D seismic repeatability and CO2 detection limits using a sparse permanent land array at the Aquistore CO2 storage site

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. WA1-WA13 ◽  
Author(s):  
Lisa A. N. Roach ◽  
Donald J. White ◽  
Brian Roberts
Keyword(s):  
Co2 Storage ◽  
Time Lapse ◽  
Data Sets ◽  
Storage Site ◽  
Processing Sequence ◽  
Time Migration ◽  
Receiver Array ◽  
Co2 Detection ◽  
4D Seismic ◽  
The Impact

Two 3D time-lapse seismic surveys were acquired in 2012 and 2013 at the Aquistore [Formula: see text] storage site prior to the start of [Formula: see text] injection. Using these surveys, we determined the background time-lapse noise at the site and assessed the feasibility of using a sparse areal permanent receiver array as a monitoring tool. Applying a standard processing sequence to these data, we adequately imaged the reservoir at 3150–3350 m depth. Evaluation of the impact of each processing step on the repeatability revealed a general monotonic increase in similarity between the data sets as a function of processing. The prestack processing sequence reduced the normalized root mean squared difference (nrms) from 1.13 between the raw stacks to 0.13 after poststack time migration. The postmigration cross-equalization sequence further reduced the global nrms to 0.07. A simulation of the changes in seismic response due to a range of [Formula: see text] injection scenarios suggested that [Formula: see text] was detectable within the reservoir at the Aquistore site provided that zones of greater thickness than 6–13 m have reached [Formula: see text] saturations of greater than 5%.

Time-lapse seismic signal analysis for enhanced oil recovery at Cranfield CO2 sequestration site, Cranfield field, Mississippi

2013 ◽  
Vol 1 (2) ◽  
pp. T157-T166 ◽  
Author(s):  
Julie Ditkof ◽  
Eva Caspari ◽  
Roman Pevzner ◽  
Milovan Urosevic ◽  
Timothy A. Meckel ◽  
...  
Keyword(s):  
Oil Recovery ◽  
Time Lapse ◽  
Seismic Signal ◽  
Injection Well ◽  
Saturation Curve ◽  
Seismic Surveys ◽  
Seismic Amplitude ◽  
The Difference ◽  
Seismic Volumes ◽  
4D Seismic

The Cranfield field in southwest Mississippi has been under continuous [Formula: see text] injection by Denbury Onshore LLC since 2008. Two 3D seismic surveys were collected in 2007 and 2010. An initial 4D seismic response was characterized after three years of injection, where more than three million tons of [Formula: see text] remain in the subsurface. This interpretation showed coherent seismic amplitude anomalies in some areas that received large amounts of [Formula: see text] but not in others. To understand these effects better, we performed Gassmann substitution modeling at two wells: the 31F-2 observation well and the 28-1 injection well. We aimed to predict a postinjection saturation curve and acoustic impedance (AI) change through the reservoir. Seismic volumes were cross-equalized, well ties to seismic were performed, and AI inversions were subsequently carried out. Inversion results showed that the change in AI is higher than Gassmann substitution predicted for the 28-1 injection well. The time-lapse AI difference predicted by the inversion is similar in magnitude to the difference inferred from a time delay along a marker horizon below the reservoir.

scholarly journals Time-lapse difference static correction using prestack crosscorrelations: 4D seismic image enhancement case from Ketzin

Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. B243-B252 ◽  
Author(s):  
Peter Bergmann ◽  
Artem Kashubin ◽  
Monika Ivandic ◽  
Stefan Lüth ◽  
Christopher Juhlin
Keyword(s):  
Seismic Data ◽  
Time Lapse ◽  
Data Sets ◽  
Storage Site ◽  
Data Set ◽  
Near Surface ◽  
Static Correction ◽  
Seismic Image ◽  
Static Corrections ◽  
4D Seismic

A method for static correction of time-lapse differences in reflection arrival times of time-lapse prestack seismic data is presented. These arrival-time differences are typically caused by changes in the near-surface velocities between the acquisitions and had a detrimental impact on time-lapse seismic imaging. Trace-to-trace time shifts of the data sets from different vintages are determined by crosscorrelations. The time shifts are decomposed in a surface-consistent manner, which yields static corrections that tie the repeat data to the baseline data. Hence, this approach implies that new refraction static corrections for the repeat data sets are unnecessary. The approach is demonstrated on a 4D seismic data set from the Ketzin [Formula: see text] pilot storage site, Germany, and is compared with the result of an initial processing that was based on separate refraction static corrections. It is shown that the time-lapse difference static correction approach reduces 4D noise more effectively than separate refraction static corrections and is significantly less labor intensive.

Quantitative estimation of compaction and velocity changes using 4D impedance and traveltime changes

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 949-957 ◽  
Author(s):  
Martin Landrø ◽  
Jan Stammeijer
Keyword(s):  
Seismic Data ◽  
Quantitative Estimation ◽  
Time Lapse ◽  
Seismic Velocities ◽  
Compaction Process ◽  
Empirical Relationships ◽  
Velocity Change ◽  
Reservoir Rocks ◽  
Velocity Changes ◽  
4D Seismic

In some hydrocarbon reservoirs, severe compaction of the reservoir rocks is observed. This compaction is caused by production, and it is often associated with changes in the overburden. Time‐lapse (or 4D) seismic data are used to monitor this compaction process. Since the compaction causes changes in both layer thickness and seismic velocities, it is crucial to distinguish between the two effects. Two new seismic methods for monitoring compacting reservoirs are introduced, one based on measured seismic prestack traveltime changes, and the other based on poststack traveltime and amplitude changes. In contrast to earlier methods, these methods do not require additional empirical relationships, such as, for instance, a velocity‐porosity relationship. The uncertainties in estimates for compaction and velocity change are expressed in terms of errors in the traveltime and amplitude measurements. These errors are directly related to the quality and repeatability of time‐lapse seismic data. For a reservoir at 3000‐m depth with 9 m of compaction, and assuming a 4D timeshift error of 0.5 ms at near offset and 2 ms at far offset, we find relative uncertainty in the compaction estimate of approximately 50–60% using traveltime information only.

scholarly journals Time-lapse image-domain tomography using adjoint-state theory

2013 ◽  
Vol 2013 (1) ◽  
pp. 1-4
Author(s):  
Jeffrey Shragge ◽  
Tongning Yang ◽  
Paul Sava
Keyword(s):  
Time Lapse ◽  
State Theory ◽  
Image Domain ◽  
Adjoint State

scholarly journals Application of time-lapse ERT imaging to watershed characterization

Geophysics ◽  
2008 ◽  
Vol 73 (3) ◽  
pp. G7-G17 ◽  
Author(s):  
Carlyle R. Miller ◽  
Partha S. Routh ◽  
Troy R. Brosten ◽  
James P. McNamara
Keyword(s):  
Reference Model ◽  
Time Lapse ◽  
Data Sets ◽  
Data Set ◽  
Resistivity Tomography ◽  
Practical Applications ◽  
Data Inversion ◽  
Time Periods ◽  
The Difference ◽  
Base Data

Time-lapse electrical resistivity tomography (ERT) has many practical applications to the study of subsurface properties and processes. When inverting time-lapse ERT data, it is useful to proceed beyond straightforward inversion of data differences and take advantage of the time-lapse nature of the data. We assess various approaches for inverting and interpreting time-lapse ERT data and determine that two approaches work well. The first approach is model subtraction after separate inversion of the data from two time periods, and the second approach is to use the inverted model from a base data set as the reference model or prior information for subsequent time periods. We prefer this second approach. Data inversion methodology should be consideredwhen designing data acquisition; i.e., to utilize the second approach, it is important to collect one or more data sets for which the bulk of the subsurface is in a background or relatively unperturbed state. A third and commonly used approach to time-lapse inversion, inverting the difference between two data sets, localizes the regions of the model in which change has occurred; however, varying noise levels between the two data sets can be problematic. To further assess the various time-lapse inversion approaches, we acquired field data from a catchment within the Dry Creek Experimental Watershed near Boise, Idaho, U.S.A. We combined the complimentary information from individual static ERT inversions, time-lapse ERT images, and available hydrologic data in a robust interpretation scheme to aid in quantifying seasonal variations in subsurface moisture content.

scholarly journals Initial 4D seismic results after CO2 injection start-up at the Aquistore storage site

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. B95-B107 ◽  
Author(s):  
Lisa A. N. Roach ◽  
Donald J. White ◽  
Brian Roberts ◽  
Doug Angus
Keyword(s):  
Root Mean Square ◽  
Time Lapse ◽  
Seismic Survey ◽  
Co2 Injection ◽  
Data Sets ◽  
Storage Site ◽  
Mean Square ◽  
Simultaneous Processing ◽  
Start Up ◽  
4D Seismic

The first post-[Formula: see text]-injection 3D time-lapse seismic survey was conducted at the Aquistore [Formula: see text] storage site in February 2016 using the same permanent array of buried geophones used for acquisition of three previous pre-[Formula: see text]-injection surveys from March 2012 to November 2013. By February 2016, 36 kilotons of [Formula: see text] have been injected within the reservoir between 3170 and 3370 m depth. We have developed time-lapse results from analysis of the first post-[Formula: see text]-injection data and three pre-[Formula: see text]-injection data sets. The objective of our analysis was to evaluate the ability of the permanent array to detect the injected [Formula: see text]. A “4D-friendly simultaneous” processing flow was applied to the data in an effort to maximize the repeatability between the pre- and post-[Formula: see text]-injection volumes while optimizing the final subsurface image including the reservoir. Excellent repeatability was achieved among all surveys with global normalized root-mean-square (Gnrms) values of 1.13–1.19 for the raw prestack data relative to the baseline data, which decreased during processing to Gnrms values of approximately 0.10 for the final crossequalized migrated data volumes. A zone of high normalized root-mean-square (nrms) values (0.11–0.25 as compared with background values of 0.05–0.10) is identified within the upper Deadwood unit of the storage reservoir, which likely corresponds to approximately 18 kilotons of [Formula: see text]. No significant nrms anomalies are observed within the other reservoir units due to a combination of reduced seismic sensitivity, higher background nrms values, and/or small quantities of [Formula: see text] residing within these zones.

Sign in / Sign up

Export Citation Format

Share Document