Data-driven Prestack Waveform Inversion Using Genetic Algorithm- Methodology and Examples

2021 ◽  
pp. 1-97
Author(s):  
Lingxiao Jia ◽  
Subhashis Mallick ◽  
Cheng Wang

The choice of an initial model for seismic waveform inversion is important. In matured exploration areas with adequate well control, we can generate a suitable initial model using well information. However, in new areas where well control is sparse or unavailable, such an initial model is compromised and/or biased by the regions with more well controls. Even in matured exploration areas, if we use time-lapse seismic data to predict dynamic reservoir properties, an initial model, that we obtain from the existing preproduction wells could be incorrect. In this work, we outline a new methodology and workflow for a nonlinear prestack isotropic elastic waveform inversion. We call this method a data driven inversion, meaning that we derive the initial model entirely from the seismic data without using any well information. By assuming a locally horizonal stratification for every common midpoint and starting from the interval P-wave velocity, estimated entirely from seismic data, our method generates pseudo wells by running a two-pass one-dimensional isotropic elastic prestack waveform inversion that uses the reflectivity method for forward modeling and genetic algorithm for optimization. We then use the estimated pseudo wells to build the initial model for seismic inversion. By applying this methodology to real seismic data from two different geological settings, we demonstrate the usefulness of our method. We believe that our new method is potentially applicable for subsurface characterization in areas where well information is sparse or unavailable. Additional research is however necessary to improve the compute-efficiency of the methodology.

2014 ◽  
Vol 54 (1) ◽  
pp. 69
Author(s):  
Andrew Long ◽  
Cyrille Reiser

Ultra-low seismic frequencies less than about 7 Hz cannot be produced by conventional air gun arrays, for any configuration and for any towing depth. There is a profound difference between improving low-frequency recovery by removing source and receiver ghosts (achievable) and improving low-frequency injection on the source side (an unrealised dream). If 1–7 Hz amplitudes could be usefully injected into the earth, it would be possible to facilitate much sharper seismic representation of geological contacts and internal features, and seismic inversion would yield robust and precise predictions of reservoir properties—without well control. The net result is fewer exploration and appraisal wells, greatly reduced exploration and development risks, and optimised recoverable reserves. Furthermore, an emerging seismic pursuit known as full waveform inversion (FWI) makes the bold promise that raw seismic field gathers can be directly used to invert for the highest achievable velocity models, almost without any human intervention. These models will bypass the traditional lack of low-frequency information in band-limited seismic data, and facilitate the aforementioned ambition of seismic inversion without well control. FWI, however, is confronted by the paradox that ultra-low-frequency seismic gathers are the necessary input for stable results. This paper describes new technologies that may enable the injection of strong 2–7 Hz amplitudes into the earth, and explains in simple terms how FWI can already be pursued as a robust complement to the prediction of accurate reservoir properties. The low-frequency revolution is already here.


Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. B1-B12 ◽  
Author(s):  
Josiane Pafeng ◽  
Subhashis Mallick ◽  
Hema Sharma

Applying seismic inversion to estimate subsurface elastic earth properties for reservoir characterization is a challenge in exploration seismology. In recent years, waveform-based seismic inversions have gained popularity, but due to high computational costs, their applications are limited, and amplitude-variation-with-offset/angle inversion is still the current state-of-the-art. We have developed a genetic-algorithm-based prestack seismic waveform inversion methodology. By parallelizing at multiple levels and assuming a locally 1D structure such that forward computation of wave equation synthetics is computationally efficient, this method is capable of inverting 3D prestack seismic data on parallel computers. Applying this inversion to a real prestack seismic data volume from the Rock Springs Uplift (RSU) located in Wyoming, USA, we determined that our method is capable of inverting the data in a reasonable runtime and producing much higher quality results than amplitude-variation-with-offset/angle inversion. Because the primary purpose for seismic data acquisition at the RSU was to characterize the subsurface for potential targets for carbon dioxide sequestration, we also identified and analyzed some potential primary and secondary storage formations and their associated sealing lithologies from our inversion results.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. B229-B240 ◽  
Author(s):  
Rajive Kumar ◽  
Prashant Bansal ◽  
Bader S. Al-Mal ◽  
Sagnik Dasgupta ◽  
Colin Sayers ◽  
...  

Optimization of production from unconventional reservoirs requires estimates of reservoir properties such as porosity, total organic carbon (TOC) content, clay content, fluid saturation, and fracture intensity. The porosity and TOC content help to determine reservoir quality, and the natural fracture intensity provides information important for the completion strategy. Because shale reservoirs display intrinsic anisotropy due to layering and the partial alignment of clay minerals and kerogen with the bedding plane, the minimum acceptable representation of the anisotropy of naturally fractured shale-gas reservoirs is orthotropy, in which a set of vertical compliant fractures is embedded in a vertical transverse isotropic (VTI) background medium. Full-azimuth seismic data are required to characterize such reservoirs and to invert for the anisotropic elastic properties. Orthotropic inversion uses azimuthally sectored seismic data stacked according to the incident angle. Even for high-fold acquisition, this azimuth/angle grouping can result in low-fold angle stacks. Orthotropic amplitude-variation-with-offset-and-azimuth (AVOAz) inversion requires seismic preconditioning techniques that ensure proper primary amplitude preservation, noise attenuation, and data alignment, and a workflow implemented for the construction of an orthotropic rock-physics model. This model integrates well and core data to estimate reservoir properties using the results of the AVOAz inversion. The seismic inversion results include the P- and S-impedance and parameters quantifying the azimuthal anisotropy. The rock model assumes a VTI kerogen-rich layer, containing aligned vertical fractures, and it uses prestack orthotropic AVOAz inversion results to predict porosity, TOC, and fracture intensity.


2022 ◽  
Vol 9 ◽  
Author(s):  
Zhonghan Liu ◽  
Yingcai Zheng ◽  
Hua-Wei Zhou

To better interpret the subsurface structures and characterize the reservoir, a depth model quantifying P-wave velocity together with additional rock’s physical parameters such as density, the S-wave velocity, and anisotropy is always preferred by geologists and engineers. Tradeoffs among different parameters can bring extra challenges to the seismic inversion process. In this study, we propose and test the Direct Waveform Inversion (DWI) scheme to simultaneously invert for 1D layered velocity and density profiles, using reflection seismic waveforms recorded on the surface. The recorded data includes primary reflections and interbed multiples. DWI is implemented in the time-space domain then followed by a wavefield extrapolation to downward continue the source and receiver. By explicitly enforcing the wavefield time-space causality, DWI can recursively determine the subsurface seismic structure in a local layer-by-layer fashion for both sharp interfaces and the properties of the layers, from shallow to deep depths. DWI is different from the layer stripping methods in the frequency domain. By not requiring a global initial model, DWI also avoids many nonlinear optimization problems, such as the local minima or the need for an accurate initial model in most waveform inversion schemes. Two numerical tests show the validity of this DWI scheme serving as a new strategy for multi-parameter seismic inversion.


Author(s):  
Michael Gineste ◽  
Jo Eidsvik

AbstractAn ensemble-based method for seismic inversion to estimate elastic attributes is considered, namely the iterative ensemble Kalman smoother. The main focus of this work is the challenge associated with ensemble-based inversion of seismic waveform data. The amount of seismic data is large and, depending on ensemble size, it cannot be processed in a single batch. Instead a solution strategy of partitioning the data recordings in time windows and processing these sequentially is suggested. This work demonstrates how this partitioning can be done adaptively, with a focus on reliable and efficient estimation. The adaptivity relies on an analysis of the update direction used in the iterative procedure, and an interpretation of contributions from prior and likelihood to this update. The idea is that these must balance; if the prior dominates, the estimation process is inefficient while the estimation is likely to overfit and diverge if data dominates. Two approaches to meet this balance are formulated and evaluated. One is based on an interpretation of eigenvalue distributions and how this enters and affects weighting of prior and likelihood contributions. The other is based on balancing the norm magnitude of prior and likelihood vector components in the update. Only the latter is found to sufficiently regularize the data window. Although no guarantees for avoiding ensemble divergence are provided in the paper, the results of the adaptive procedure indicate that robust estimation performance can be achieved for ensemble-based inversion of seismic waveform data.


2019 ◽  
Vol 38 (2) ◽  
pp. 106-115 ◽  
Author(s):  
Phuong Hoang ◽  
Arcangelo Sena ◽  
Benjamin Lascaud

The characterization of shale plays involves an understanding of tectonic history, geologic settings, reservoir properties, and the in-situ stresses of the potential producing zones in the subsurface. The associated hydrocarbons are generally recovered by horizontal drilling and hydraulic fracturing. Historically, seismic data have been used mainly for structural interpretation of the shale reservoirs. A primary benefit of surface seismic has been the ability to locate and avoid drilling into shallow carbonate karsting zones, salt structures, and basement-related major faults which adversely affect the ability to drill and complete the well effectively. More recent advances in prestack seismic data analysis yield attributes that appear to be correlated to formation lithology, rock strength, and stress fields. From these, we may infer preferential drilling locations or sweet spots. Knowledge and proper utilization of these attributes may prove valuable in the optimization of drilling and completion activities. In recent years, geophysical data have played an increasing role in supporting well planning, hydraulic fracturing, well stacking, and spacing. We have implemented an integrated workflow combining prestack seismic inversion and multiattribute analysis, microseismic data, well-log data, and geologic modeling to demonstrate key applications of quantitative seismic analysis utilized in developing ConocoPhillips' acreage in the Delaware Basin located in Texas. These applications range from reservoir characterization to well planning/execution, stacking/spacing optimization, and saltwater disposal. We show that multidisciplinary technology integration is the key for success in unconventional play exploration and development.


2001 ◽  
Vol 41 (2) ◽  
pp. 131
Author(s):  
A.G. Sena ◽  
T.M. Smith

The successful exploration for new reservoirs in mature areas, as well as the optimal development of existing fields, requires the integration of unconventional geological and geophysical techniques. In particular, the calibration of 3D seismic data to well log information is crucial to obtain a quantitative understanding of reservoir properties. The advent of new technology for prestack seismic data analysis and 3D visualisation has resulted in improved fluid and lithology predictions prior to expensive drilling. Increased reservoir resolution has been achieved by combining seismic inversion with AVO analysis to minimise exploration risk.In this paper we present an integrated and systematic approach to prospect evaluation in an oil/gas field. We will show how petrophysical analysis of well log data can be used as a feasibility tool to determine the fluid and lithology discrimination capabilities of AVO and inversion techniques. Then, a description of effective AVO and prestack inversion tools for reservoir property quantification will be discussed. Finally, the incorporation of the geological interpretation and the use of 3D visualisation will be presented as a key integration tool for the discovery of new plays.


2019 ◽  
Author(s):  
Clàudia Gras ◽  
Valentí Sallarès ◽  
Daniel Dagnino ◽  
C. Estela Jiménez ◽  
Adrià Meléndez ◽  
...  

Abstract. We present a high-resolution P-wave velocity model of the sedimentary cover and the uppermost basement until ~ 3 km depth obtained by full-waveform inversion of multichannel seismic data acquired with a 6 km-long streamer in the Alboran Sea (SE Iberia). The inherent non-linearity of the method, especially for short-offset, band-limited seismic data as this one, is circumvented by applying a data processing/modeling sequence consisting of three steps: (1) data re-datuming by back-propagation of the recorded seismograms to the seafloor; (2) joint refraction and reflection travel-time tomography combining the original and the re-datumed shot gathers; and (3) FWI of the original shot gathers using the model obtained by travel-time tomography as initial reference. The final velocity model shows a number of geological structures that cannot be identified in the travel-time tomography models or easily interpreted from seismic reflection images alone. A sharp strong velocity contrast accurately defines the geometry of the top of the basement. Several low-velocity zones that may correspond to the abrupt velocity change across steeply dipping normal faults are observed at the flanks of the basin. A 200–300 m thick, high-velocity layer embedded within lower velocity sediment may correspond to evaporites deposited during the Messinian crisis. The results confirm that the combination of data re-datuming and joint refraction and reflection travel-time inversion provides reference models that are accurate enough to apply full-waveform inversion to relatively short offset streamer data in deep water settings starting at field-data standard low frequency content of 6 Hz.


Sign in / Sign up

Export Citation Format

Share Document