Double-difference elastic-waveform inversion with prior information for time-lapse monitoring

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. R259-R273 ◽  
Author(s):  
Zhigang Zhang ◽  
Lianjie Huang

Quantitative time-lapse seismic monitoring can provide crucial information for enhanced oil recovery, geologic carbon storage, and enhanced geothermal systems. Recently developed double-difference elastic-waveform inversion has the potential to quantitatively monitor reservoirs using seismic reflection data. Because the approximate location of a reservoir or a target monitoring region is usually known, we incorporated this knowledge as prior information into double-difference elastic-waveform inversion. Using numerical examples of synthetic time-lapse models, we found that our new method can quantitatively monitor the changes of elastic properties within reservoirs. Therefore, the double-difference elastic-waveform inversion with prior information on the location of a monitoring region is a promising tool for quantitatively monitoring reservoir properties’ changes.

Geophysics ◽  
2021 ◽  
pp. 1-51
Author(s):  
Yanhua Liu ◽  
Ilya Tsvankin

Time-lapse full-waveform inversion can provide high-resolution information about changes in the reservoir properties during hydrocarbon production and CO2 injection. However, the accuracy of the estimated source wavelet, which is critically important for time-lapse FWI, is often insufficient for field-data applications. The so-called “source-independent” FWI is designed to reduce the influence of the source wavelet on the inversion results. We incorporate the convolution-based source-independent technique into a time-lapse FWI algorithm for VTI (transversely isotropic with a vertical symmetry axis) media. The gradient of the modified FWI objective function is obtained from the adjoint-state method. The algorithm is tested on a model with a graben structure and the modified VTI Marmousi model using three time-lapse strategies (the parallel-difference, sequential-difference, and double-difference methods). The results confirm the ability of the developed methodology to reconstruct the localized time-lapse parameter variations even for a strongly distorted source wavelet. The algorithm remains robust in the presence of moderate noise in the input data but the accuracy of the estimated time-lapse changes depends on the model complexity.


Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 845-860 ◽  
Author(s):  
François Clément ◽  
Guy Chavent ◽  
Susana Gómez

Migration‐based traveltime (MBTT) formulation provides algorithms for automatically determining background velocities from full‐waveform surface seismic reflection data using local optimization methods. In particular, it addresses the difficulty of the nonconvexity of the least‐squares data misfit function. The method consists of parameterizing the reflectivity in the time domain through a migration step and providing a multiscale representation for the smooth background velocity. We present an implementation of the MBTT approach for a 2-D finite‐difference (FD) full‐wave acoustic model. Numerical analysis on a 2-D synthetic example shows the ability of the method to find much more reliable estimates of both long and short wavelengths of the velocity than the classical least‐squares approach, even when starting from very poor initial guesses. This enlargement of the domain of attraction for the global minima of the least‐squares misfit has a price: each evaluation of the new objective function requires, besides the usual FD full‐wave forward modeling, an additional full‐wave prestack migration. Hence, the FD implementation of the MBTT approach presented in this paper is expected to provide a useful tool for the inversion of data sets of moderate size.


Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. R31-R41 ◽  
Author(s):  
Andrea Zunino ◽  
Klaus Mosegaard ◽  
Katrine Lange ◽  
Yulia Melnikova ◽  
Thomas Mejer Hansen

Determination of a petroleum reservoir structure and rock bulk properties relies extensively on inference from reflection seismology. However, classic deterministic methods to invert seismic data for reservoir properties suffer from some limitations, among which are the difficulty of handling complex, possibly nonlinear forward models, and the lack of robust uncertainty estimations. To overcome these limitations, we studied a methodology to invert seismic reflection data in the framework of the probabilistic approach to inverse problems, using a Markov chain Monte Carlo (McMC) algorithm with the goal to directly infer the rock facies and porosity of a target reservoir zone. We thus combined a rock-physics model with seismic data in a single inversion algorithm. For large data sets, the McMC method may become computationally impractical, so we relied on multiple-point-based a priori information to quantify geologically plausible models. We tested this methodology on a synthetic reservoir model. The solution of the inverse problem was then represented by a collection of facies and porosity reservoir models, which were samples of the posterior distribution. The final product included probability maps of the reservoir properties in obtained by performing statistical analysis on the collection of solutions.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE25-VE34 ◽  
Author(s):  
Gilles Lambaré

Stereotomography was proposed [Formula: see text] ago for estimating velocity macromodels from seismic reflection data. Initially, the goal was to retain the advantages of standard traveltime tomography while providing an alternative to difficult interpretive traveltime picking. Stereotomography relies on the concept of locally coherent events characterized by their local slopes in the prestack data cube. Currently, stereotomography has been developed in two and three dimensions, and precious experience has been gained. The expected advantages have been demonstrated fully (in particular, the efficiency and reliability of the semiautomatic stereotomographic picking strategies), and further studies have increased the method’s potential and flexibility. For example, stereotomographic picking can now be done in either the prestack or poststack domain, in either the time (migrated or unmigrated) or depth domain. It appears that the theoretical frame of stereotomography can reconcile, very satisfactorily and efficiently, most methods proposed for velocity-macromodel estimation for depth imaging. Moreover, an extension of the method to full-waveform inversion already exists and opens the way for very interesting developments.


Geophysics ◽  
2011 ◽  
Vol 76 (1) ◽  
pp. SA19-SA34 ◽  
Author(s):  
Shohei Minato ◽  
Toshifumi Matsuoka ◽  
Takeshi Tsuji ◽  
Deyan Draganov ◽  
Jürg Hunziker ◽  
...  

Crosswell reflection method is a high-resolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to time-lapse surveys. An alternative is to use surface sources in combination with seismic interferometry. Seismic interferometry (SI) could retrieve the reflection response at one of the boreholes as if from a source inside the other borehole. We investigate the applicability of SI for the retrieval of the reflection response between two boreholes using numerically modeled field data. We compare two SI approaches — crosscorrelation (CC) and multidimensional deconvolution (MDD). SI by MDD is less sensitive to underillumination from the source distribution, but requires inversion of the recordings at one of the receiver arrays from all the available sources. We find that the inversion problem is ill-posed, and propose to stabilize it using singular-value decomposition. The results show that the reflections from deep boundaries are retrieved very well using both the CC and MDD methods. Furthermore, the MDD results exhibit more realistic amplitudes than those from the CC method for downgoing reflections from shallow boundaries. We find that the results retrieved from the application of both methods to field data agree well with crosswell seismic-reflection data using borehole sources and with the logged P-wave velocity.


Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. P13-P25
Author(s):  
Michael J. Faggetter ◽  
Mark E. Vardy ◽  
Justin K. Dix ◽  
Jonathan M. Bull ◽  
Timothy J. Henstock

Time-lapse (4D) seismic imaging is now widely used as a tool to map and interpret changes in deep reservoirs as well as investigate dynamic, shallow hydrological processes in the near surface. However, there are very few examples of time-lapse analysis using ultra-high-frequency (UHF; kHz range) marine seismic reflection data. Exacting requirements for navigation can be prohibitive for acquiring coherent, true-3D volumes. Variable environmental noise can also lead to poor amplitude repeatability and make it difficult to identify differences that are related to real physical changes. Overcoming these challenges opens up a range of potential applications for monitoring the subsurface at decimetric resolution, including geohazards, geologic structures, as well as the bed-level and subsurface response to anthropogenic activities. Navigation postprocessing was incorporated to improve the acquisition and processing workflow for the 3D Chirp subbottom profiler and provide stable, centimeter-level absolute positioning, resulting in well-matched 3D data and mitigating 4D noise for data stacked into [Formula: see text] common-midpoint bins. Within an example 4D data set acquired on the south coast of the UK, interpretable differences are recorded within a shallow gas blanket. Reflections from the top and bottom of a gas pocket are imaged at low tide, whereas at high tide only the upper reflection is imaged. This case study demonstrates the viability of time-lapse UHF 3D seismic reflection for quantitative mapping of decimeter-scale changes within the shallow marine subsurface.


Sign in / Sign up

Export Citation Format

Share Document