Stochastic inversion of pressure and saturation changes from time-lapse AVO data

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. C81-C92 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Hilde Grude Borgos ◽  
Martin Landrø
Keyword(s):  
Stochastic Model ◽  
Seismic Data ◽  
Synthetic Data ◽  
Time Lapse ◽  
Bayesian Framework ◽  
Reservoir Modeling ◽  
Data Set ◽  
The North ◽  
Fluid Saturation ◽  
Likelihood Model

Effects of pressure and fluid saturation can have the same degree of impact on seismic amplitudes and differential traveltimes in the reservoir interval; thus, they are often inseparable by analysis of a single stacked seismic data set. In such cases, time-lapse AVO analysis offers an opportunity to discriminate between the two effects. We quantify the uncertainty in estimations to utilize information about pressure- and saturation-related changes in reservoir modeling and simulation. One way of analyzing uncertainties is to formulate the problem in a Bayesian framework. Here, the solution of the problem will be represented by a probability density function (PDF), providing estimations of uncertainties as well as direct estimations of the properties. A stochastic model for estimation of pressure and saturation changes from time-lapse seismic AVO data is investigated within a Bayesian framework. Well-known rock physical relationships are used to set up a prior stochastic model. PP reflection coefficient differences are used to establish a likelihood model for linking reservoir variables and time-lapse seismic data. The methodology incorporates correlation between different variables of the model as well as spatial dependencies for each of the variables. In addition, information about possible bottlenecks causing large uncertainties in the estimations can be identified through sensitivity analysis of the system. The method has been tested on 1D synthetic data and on field time-lapse seismic AVO data from the Gullfaks Field in the North Sea.

Discrimination between pressure and fluid saturation changes from marine multicomponent time‐lapse seismic data

Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1592-1599 ◽  
Author(s):  
Martin Landrø ◽  
Helene Hafslund Veire ◽  
Kenneth Duffaut ◽  
Nazih Najjar
Keyword(s):  
Pore Pressure ◽  
Seismic Data ◽  
Water Saturation ◽  
Synthetic Data ◽  
Time Lapse ◽  
Seismic Survey ◽  
Data Sets ◽  
Amplitude Variation With Offset ◽  
Fluid Saturation ◽  
Pressure Changes

Explicit expressions for computation of saturation and pressure‐related changes from marine multicomponent time‐lapse seismic data are presented. Necessary input is PP and PS stacked data for the baseline seismic survey and the repeat survey. Compared to earlier methods based on PP data only, this method is expected to be more robust since two independent measurements are used in the computation. Due to a lack of real marine multicomponent time‐lapse seismic data sets, the methodology is tested on synthetic data sets, illustrating strengths and weaknesses of the proposed technique. Testing ten scenarios for various changes in pore pressure and fluid saturation, we find that it is more robust for most cases to use the proposed 4D PP/PS technique instead of a 4D PP amplitude variation with offset (AVO) technique. The fit between estimated and “real” changes in water saturation and pore pressure were good for most cases. On the average, we find that the deviation in estimated saturation changes is 8% and 0.3 MPa for the estimated pore pressure changes. For PP AVO, we find that the corresponding average errors are 9% and 1.0 MPa. In the present method, only 4D PP and PS amplitude changes are used in the calculations. It is straightforward to include use of 4D traveltime shifts in the algorithm and, if reliable time shifts can be measured, this will most likely further stabilize the presented method.

Prediction of residual oil saturation by using the ratio of amplitude attributes of time-lapse seismic data

Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. IM1-IM12 ◽  
Author(s):  
Meng Li ◽  
Zhen Liu ◽  
Minzhu Liu ◽  
Huilai Zhang
Keyword(s):  
Seismic Data ◽  
Time Lapse ◽  
Seismic Interpretation ◽  
Impedance Change ◽  
Residual Oil ◽  
Data Set ◽  
Oil Saturation ◽  
Residual Oil Saturation ◽  
Amplitude Difference ◽  
Fluid Saturation

Subtraction of baseline and monitoring seismic data is a common step in highlighting reservoir changes in time-lapse seismic interpretation. However, ambiguity exists in the interpretation of the amplitude difference, which is controlled by fluid change and reservoir thickness. To estimate the residual oil saturation quantitatively, we have developed a time-lapse seismic interpretation method that uses the ratio of amplitude attributes extracted from the baseline and monitoring seismic data. The relationship between impedance change and the ratio of the baseline and monitoring amplitude attributes is determined to avoid the influence of reservoir thickness. Subsequently, the fluid saturation is calculated from the impedance change by using a proper petrophysical relationship. We have tested our new method on a real time-lapse seismic data set from a water-flooded reservoir in the deepwater area of West Africa. The water-flooded area determined from the amplitude difference does not completely match the production logs because of the influence of variations in the reservoir thickness. However, the residual oil distribution calculated with the proposed method matches the production logs well. The connectivity of sandstone bodies is also evaluated based on an integrated interpretation of estimated oil saturation. With its simple principles and easy accessibility, our method improves the accuracy of time-lapse seismic data interpretation in water-flooded oil reservoirs. Furthermore, the quantitative interpretation of fluid change enables the time-lapse seismic technology to guide reservoir development directly.

Discrimination between pressure and fluid saturation changes from time‐lapse seismic data

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 836-844 ◽  
Author(s):  
Martin Landrø
Keyword(s):  
Seismic Data ◽  
Time Lapse ◽  
Seismic Survey ◽  
Baseline Survey ◽  
Water Contact ◽  
Data Set ◽  
Seismic Surveys ◽  
Fluid Saturation ◽  
Oil Water ◽  
Pressure Changes

Explicit expressions for computing saturation‐ and pressure‐related changes from time‐lapse seismic data have been derived and tested on a real time‐lapse seismic data set. Necessary input is near‐and far‐offset stacks for the baseline seismic survey and the repeat survey. The method has been tested successfully in a segment where pressure measurements in two wells verify a pore‐pressure increase of 5 to 6 MPa between the baseline survey and the monitor survey. Estimated pressure changes using the proposed relationships fit very well with observations. Between the baseline and monitor seismic surveys, 27% of the estimated recoverable hydrocarbon reserves were produced from this segment. The estimated saturation changes also agree well with observed changes, apart from some areas in the water zone that are mapped as being exposed to saturation changes (which is unlikely). Saturation changes in other segments close to the original oil‐water contact and the top reservoir interface are also estimated and confirmed by observations in various wells.

Regularization and datuming of seismic data by weighted, damped least squares

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Synthetic Data ◽  
Real Data ◽  
Structural Data ◽  
Irregular Sampling ◽  
Data Set ◽  
Near Surface ◽  
Damped Least Squares ◽  
Wavefield Extrapolation

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.

Sparse reflectivity inversion for nonstationary seismic data with surface-related multiples: Numerical and field-data experiments

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. R199-R217 ◽  
Author(s):  
Xintao Chai ◽  
Shangxu Wang ◽  
Genyang Tang
Keyword(s):  
Seismic Data ◽  
Resolution Enhancement ◽  
Synthetic Data ◽  
Data Sets ◽  
Data Set ◽  
Anelastic Attenuation ◽  
Seismic Resolution ◽  
Text Filtering ◽  
The Stability ◽  
Reflectivity Inversion

Seismic data are nonstationary due to subsurface anelastic attenuation and dispersion effects. These effects, also referred to as the earth’s [Formula: see text]-filtering effects, can diminish seismic resolution. We previously developed a method of nonstationary sparse reflectivity inversion (NSRI) for resolution enhancement, which avoids the intrinsic instability associated with inverse [Formula: see text] filtering and generates superior [Formula: see text] compensation results. Applying NSRI to data sets that contain multiples (addressing surface-related multiples only) requires a demultiple preprocessing step because NSRI cannot distinguish primaries from multiples and will treat them as interference convolved with incorrect [Formula: see text] values. However, multiples contain information about subsurface properties. To use information carried by multiples, with the feedback model and NSRI theory, we adapt NSRI to the context of nonstationary seismic data with surface-related multiples. Consequently, not only are the benefits of NSRI (e.g., circumventing the intrinsic instability associated with inverse [Formula: see text] filtering) extended, but also multiples are considered. Our method is limited to be a 1D implementation. Theoretical and numerical analyses verify that given a wavelet, the input [Formula: see text] values primarily affect the inverted reflectivities and exert little effect on the estimated multiples; i.e., multiple estimation need not consider [Formula: see text] filtering effects explicitly. However, there are benefits for NSRI considering multiples. The periodicity and amplitude of the multiples imply the position of the reflectivities and amplitude of the wavelet. Multiples assist in overcoming scaling and shifting ambiguities of conventional problems in which multiples are not considered. Experiments using a 1D algorithm on a synthetic data set, the publicly available Pluto 1.5 data set, and a marine data set support the aforementioned findings and reveal the stability, capabilities, and limitations of the proposed method.

scholarly journals Time-Lapse Seismic Monitoring of Onshore Reservoirs in Niger Delta, Field ‘K’ as a Case Study

2017 ◽  
Vol 1 (1) ◽  
pp. 23-27 ◽  
Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Niger Delta ◽  
Rock Physics ◽  
Time Lapse ◽  
Seismic Survey ◽  
Impedance Change ◽  
Data Set ◽  
Reservoir Fluid ◽  
4D Seismic

Field ‘K’ is a mature field in the coastal swamp onshore Niger delta, which has been producing since 1960. As a huge producing field with some potential for further sustainable production, field monitoring is therefore important in the identification of areas of unproduced hydrocarbon. This can be achieved by comparing production data with the corresponding changes in acoustic impedance observed in the maps generated from base survey (initial 3D seismic) and monitor seismic survey (4D seismic) across the field. This will enable the 4D seismic data set to be used for mapping reservoir details such as advancing water front and un-swept zones. The availability of good quality onshore time-lapse seismic data for Field ‘K’ acquired in 1987 and 2002 provided the opportunity to evaluate the effect of changes in reservoir fluid saturations on time-lapse amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic impedance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir fluid saturations were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where reservoir porosity is 0.25%. In this paper, we demonstrated the use of repeat Seismic to delineate swept zones and areas hit with water override in a producing onshore reservoir.

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.

Marlim R3D: A realistic model for controlled-source electromagnetic simulations — Phase 2: The controlled-source electromagnetic data set

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. E293-E299
Author(s):  
Jorlivan L. Correa ◽  
Paulo T. L. Menezes
Keyword(s):  
A Priori ◽  
Synthetic Data ◽  
Realistic Model ◽  
Earth Model ◽  
Data Sets ◽  
Data Set ◽  
Geoelectric Model ◽  
Controlled Source ◽  
The North ◽  
Electromagnetic Simulations

Synthetic data provided by geoelectric earth models are a powerful tool to evaluate a priori a controlled-source electromagnetic (CSEM) workflow effectiveness. Marlim R3D (MR3D) is an open-source complex and realistic geoelectric model for CSEM simulations of the postsalt turbiditic reservoirs at the Brazilian offshore margin. We have developed a 3D CSEM finite-difference time-domain forward study to generate the full-azimuth CSEM data set for the MR3D earth model. To that end, we fabricated a full-azimuth survey with 45 towlines striking the north–south and east–west directions over a total of 500 receivers evenly spaced at 1 km intervals along the rugged seafloor of the MR3D model. To correctly represent the thin, disconnected, and complex geometries of the studied reservoirs, we have built a finely discretized mesh of [Formula: see text] cells leading to a large mesh with a total of approximately 90 million cells. We computed the six electromagnetic field components (Ex, Ey, Ez, Hx, Hy, and Hz) at six frequencies in the range of 0.125–1.25 Hz. In our efforts to mimic noise in real CSEM data, we summed to the data a multiplicative noise with a 1% standard deviation. Both CSEM data sets (noise free and noise added), with inline and broadside geometries, are distributed for research or commercial use, under the Creative Common License, at the Zenodo platform.

Seismic characterization of a blocky mass-transport deposit in the Trealla Limestone Formation, North Carnarvon Basin, Australia

2020 ◽  
Vol 8 (4) ◽  
pp. SR53-SR58
Author(s):  
Laura Ortiz-Sanguino ◽  
Javier Tellez ◽  
Heather Bedle ◽  
Dilan Martinez-Sanchez
Keyword(s):  
Mass Transport ◽  
Seismic Data ◽  
Information Management System ◽  
Data Set ◽  
Differential Dissolution ◽  
Passive Margins ◽  
The North ◽  
Seismic Characterization ◽  
The Individual ◽  
Carnarvon Basin

The deepwater Cenozoic strata in the North Carnarvon Basin, Australia, represent an interval of interest for stratigraphic studies in passive margins settings of mixed siliciclastic-carbonate environments. We have explored the geomorphological characteristics of a mass-transport deposit (MTD) within the Trealla Limestone Formation to describe in detail the differences among the blocks. To characterize the individual geometry and structural configuration of the blocks within the MTD, we used geometric seismic attributes such as coherence, curvature, dip azimuth, and dip magnitude using horizon slices and vertical profiles. The evaluation finds two types of blocks: remnant and glide (or rafted) blocks. Remnant blocks are in situ and stratigraphically continuous fragments with the underlying strata. This type of block is frequently fault-bounded and displays low deformation evidence. Glide blocks are part of the transported material detached from a paleoslope. These blocks are deformed and occasionally appear as “floating” fragments embedded within a chaotic matrix in the MTD. Glide blocks are used as kinematic indicators of the direction of deposition of MTDs. We evaluate these elements in a modern continental analog that resembles a similar setting for a better understanding of the slide occurrence. Geological feature: Glide blocks, North Carnarvon Basin, Australia Seismic appearance: Discrete angular blocks with internal reflectors Alternative interpretations: Differential dissolution in a mixed siliciclastic-carbonate environment Features with a similar appearance: Carbonate buildups, differential dissolution blocks Formation: Trealla Limestone Formation, North Carnarvon Basin Age: Early-Middle Miocene Location: Offshore Northwest Australia, North Carnarvon Basin Seismic data: Obtained from Western Australian Petroleum and Geothermal Information Management System, Draeck 3D seismic data set Analysis tools: Visualization software (Petrel 2019) and attribute performance software (AASPI 6.0)

Sign in / Sign up

Export Citation Format

Share Document