scholarly journals Streamline-Based Time-Lapse-Seismic-Data Integration Incorporating Pressure and Saturation Effects

SPE Journal ◽  
2017 ◽  
Vol 22 (04) ◽  
pp. 1261-1279 ◽  
Author(s):  
Shingo Watanabe ◽  
Jichao Han ◽  
Gill Hetz ◽  
Akhil Datta-Gupta ◽  
Michael J. King ◽  
...  

Summary We present an efficient history-matching technique that simultaneously integrates 4D repeat seismic surveys with well-production data. This approach is particularly well-suited for the calibration of the reservoir properties of high-resolution geologic models because the seismic data are areally dense but sparse in time, whereas the production data are finely sampled in time but spatially averaged. The joint history matching is performed by use of streamline-based sensitivities derived from either finite-difference or streamline-based flow simulation. For the most part, earlier approaches have focused on the role of saturation changes, but the effects of pressure have largely been ignored. Here, we present a streamline-based semianalytic approach for computing model-parameter sensitivities, accounting for both pressure and saturation effects. The novelty of the method lies in the semianalytic sensitivity computations, making it computationally efficient for high-resolution geologic models. The approach is implemented by use of a finite-difference simulator incorporating the detailed physics. Its efficacy is demonstrated by use of both synthetic and field applications. For both the synthetic and the field cases, the advantages of incorporating the time-lapse variations are clear, seen through the improved estimation of the permeability distribution, the pressure profile, the evolution of the fluid saturation, and the swept volumes.

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. M73-M87 ◽  
Author(s):  
Alvaro Rey ◽  
Eric Bhark ◽  
Kai Gao ◽  
Akhil Datta-Gupta ◽  
Richard Gibson

We have developed an efficient approach of petroleum reservoir model calibration that integrates 4D seismic surveys together with well-production data. The approach is particularly well-suited for the calibration of high-resolution reservoir properties (permeability) because the field-scale seismic data are areally dense, whereas the production data are effectively averaged over interwell spacing. The joint calibration procedure is performed using streamline-based sensitivities derived from finite-difference flow simulation. The inverted seismic data (i.e., changes in elastic impedance or fluid saturations) are distributed as a 3D high-resolution grid cell property. The sensitivities of the seismic and production surveillance data to perturbations in absolute permeability at individual grid cells are efficiently computed via semianalytical streamline techniques. We generalize previous formulations of streamline-based seismic inversion to incorporate realistic field situations such as changing boundary conditions due to infill drilling, pattern conversion, etc. A commercial finite-difference flow simulator is used for reservoir simulation and to generate the time-dependent velocity fields through which streamlines are traced and the sensitivity coefficients are computed. The commercial simulator allows us to incorporate detailed physical processes including compressibility and nonconvective forces, e.g., capillary pressure effects, while the streamline trajectories provide a rapid evaluation of the sensitivities. The efficacy of our proposed approach was tested with synthetic and field applications. The synthetic example was the Society of Petroleum Engineers benchmark Brugge field case. The field example involves waterflooding of a North Sea reservoir with multiple seismic surveys. In both cases, the advantages of incorporating the time-lapse variations were clearly demonstrated through improved estimation of the permeability heterogeneity, fluid saturation evolution, and swept and drained volumes. The value of the seismic data integration was in particular proven through the identification of the continuity in reservoir sands and barriers, and by the preservation of geologic realism in the calibrated model.


Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. M15-M31 ◽  
Author(s):  
Mingliang Liu ◽  
Dario Grana

We have developed a time-lapse seismic history matching framework to assimilate production data and time-lapse seismic data for the prediction of static reservoir models. An iterative data assimilation method, the ensemble smoother with multiple data assimilation is adopted to iteratively update an ensemble of reservoir models until their predicted observations match the actual production and seismic measurements and to quantify the model uncertainty of the posterior reservoir models. To address computational and numerical challenges when applying ensemble-based optimization methods on large seismic data volumes, we develop a deep representation learning method, namely, the deep convolutional autoencoder. Such a method is used to reduce the data dimensionality by sparsely and approximately representing the seismic data with a set of hidden features to capture the nonlinear and spatial correlations in the data space. Instead of using the entire seismic data set, which would require an extremely large number of models, the ensemble of reservoir models is iteratively updated by conditioning the reservoir realizations on the production data and the low-dimensional hidden features extracted from the seismic measurements. We test our methodology on two synthetic data sets: a simplified 2D reservoir used for method validation and a 3D application with multiple channelized reservoirs. The results indicate that the deep convolutional autoencoder is extremely efficient in sparsely representing the seismic data and that the reservoir models can be accurately updated according to production data and the reparameterized time-lapse seismic data.


2011 ◽  
Vol 14 (05) ◽  
pp. 621-633 ◽  
Author(s):  
Alireza Kazemi ◽  
Karl D. Stephen ◽  
Asghar Shams

Summary History matching of a reservoir model is always a difficult task. In some fields, we can use time-lapse (4D) seismic data to detect production-induced changes as a complement to more conventional production data. In seismic history matching, we predict these data and compare to observations. Observed time-lapse data often consist of relative measures of change, which require normalization. We investigate different normalization approaches, based on predicted 4D data, and assess their impact on history matching. We apply the approach to the Nelson field in which four surveys are available over 9 years of production. We normalize the 4D signature in a number of ways. First, we use predictions of 4D signature from vertical wells that match production, and we derive a normalization function. As an alternative, we use crossplots of the full-field prediction against observation. Normalized observations are used in an automatic-history-matching process, in which the model is updated. We analyze the results of the two normalization approaches and compare against the case of just using production data. The result shows that when we use 4D data normalized to wells, we obtain 49% reduced misfit along with 36% improvement in predictions. Also over the whole reservoir, 8 and 7% reduction of misfits for 4D seismic are obtained in history and prediction periods, respectively. When we use only production data, the production history match is improved to a similar degree (45%), but in predictions, the improvement is only 25% and the 4D seismic misfit is 10% worse. Finding the unswept areas in the reservoir is always a challenge in reservoir management. By using 4D data in history matching, we can better predict reservoir behavior and identify regions of remaining oil.


SPE Journal ◽  
2012 ◽  
Vol 18 (01) ◽  
pp. 159-171 ◽  
Author(s):  
Mario Trani ◽  
Rob Arts ◽  
Olwijn Leeuwenburgh

Summary Time-lapse seismic data provide information on the dynamics of multiphase reservoir fluid flow in places where no production data from wells are available. This information, in principle, could be used to estimate unknown reservoir properties. However, the amount, resolution, and character of the data have long posed significant challenges for quantitative use in assisted-history-matching workflows. Previous studies, therefore, have generally investigated methods for updating single models with reduced parameter-uncertainty space. Recent developments in ensemble-based history-matching methods have shown the feasibility of multimodel history and matching of production data while maintaining a full uncertainty description. Here, we introduce a robust and flexible reparameterization for interpreted fluid fronts or seismic attribute isolines that extends these developments to seismic history matching. The seismic data set is reparameterized, in terms of arrival times, at observed front positions, thereby significantly reducing the number of data while retaining essential information. A simple 1D example is used to introduce the concepts of the approach. A synthetic 3D example, with spatial complexity that is typical for many waterfloods, is examined in detail. History-matching cases based on both separate and combined use of production and seismic data are examined. It is shown that consistent multimodel history matches can be obtained without the need for reduction of the parameter space or for localization of the impact of observations. The quality of forecasts based on the history-matched models is evaluated by simulating both expected production and saturation changes throughout the field for a fixed operating strategy. It is shown that bias and uncertainty in the forecasts of production both at existing wells and in the flooded area are reduced considerably when both production and seismic data are incorporated. The proposed workflow, therefore, enables better decisions on field developments that require optimal placement of infill wells.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. C81-C92 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Hilde Grude Borgos ◽  
Martin Landrø

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.


2021 ◽  
Author(s):  
Rick Schrynemeeckers

Abstract Current offshore hydrocarbon detection methods employ vessels to collect cores along transects over structures defined by seismic imaging which are then analyzed by standard geochemical methods. Due to the cost of core collection, the sample density over these structures is often insufficient to map hydrocarbon accumulation boundaries. Traditional offshore geochemical methods cannot define reservoir sweet spots (i.e. areas of enhanced porosity, pressure, or net pay thickness) or measure light oil or gas condensate in the C7 – C15 carbon range. Thus, conventional geochemical methods are limited in their ability to help optimize offshore field development production. The capability to attach ultrasensitive geochemical modules to Ocean Bottom Seismic (OBS) nodes provides a new capability to the industry which allows these modules to be deployed in very dense grid patterns that provide extensive coverage both on structure and off structure. Thus, both high resolution seismic data and high-resolution hydrocarbon data can be captured simultaneously. Field trials were performed in offshore Ghana. The trial was not intended to duplicate normal field operations, but rather provide a pilot study to assess the viability of passive hydrocarbon modules to function properly in real world conditions in deep waters at elevated pressures. Water depth for the pilot survey ranged from 1500 – 1700 meters. Positive thermogenic signatures were detected in the Gabon samples. A baseline (i.e. non-thermogenic) signature was also detected. The results indicated the positive signatures were thermogenic and could easily be differentiated from baseline or non-thermogenic signatures. The ability to deploy geochemical modules with OBS nodes for reoccurring surveys in repetitive locations provides the ability to map the movement of hydrocarbons over time as well as discern depletion affects (i.e. time lapse geochemistry). The combined technologies will also be able to: Identify compartmentalization, maximize production and profitability by mapping reservoir sweet spots (i.e. areas of higher porosity, pressure, & hydrocarbon richness), rank prospects, reduce risk by identifying poor prospectivity areas, accurately map hydrocarbon charge in pre-salt sequences, augment seismic data in highly thrusted and faulted areas.


2021 ◽  
Author(s):  
S Al Naqbi ◽  
J Ahmed ◽  
J Vargas Rios ◽  
Y Utami ◽  
A Elila ◽  
...  

Abstract The Thamama group of reservoirs consist of porous carbonates laminated with tight carbonates, with pronounced lateral heterogeneities in porosity, permeability, and reservoir thickness. The main objective of our study was mapping variations and reservoir quality prediction away from well control. As the reservoirs were thin and beyond seismic resolution, it was vital that the facies and porosity be mapped in high resolution, with a high predictability, for successful placement of horizontal wells for future development of the field. We established a unified workflow of geostatistical inversion and rock physics to characterize the reservoirs. Geostatistical inversion was run in static models that were converted from depth to time domain. A robust two-way velocity model was built to map the depth grid and its zones on the time seismic data. This ensured correct placement of the predicted high-resolution elastic attributes in the depth static model. Rock physics modeling and Bayesian classification were used to convert the elastic properties into porosity and lithology (static rock-type (SRT)), which were validated in blind wells and used to rank the multiple realizations. In the geostatistical pre-stack inversion, the elastic property prediction was constrained by the seismic data and controlled by variograms, probability distributions and a guide model. The deterministic inversion was used as a guide or prior model and served as a laterally varying mean. Initially, unconstrained inversion was tested by keeping all wells as blind and the predictions were optimized by updating the input parameters. The stochastic inversion results were also frequency filtered in several frequency bands, to understand the impact of seismic data and variograms on the prediction. Finally, 30 wells were used as input, to generate 80 realizations of P-impedance, S-impedance, Vp/Vs, and density. After converting back to depth, 30 additional blind wells were used to validate the predicted porosity, with a high correlation of more than 0.8. The realizations were ranked based on the porosity predictability in blind wells combined with the pore volume histograms. Realizations with high predictability and close to the P10, P50 and P90 cases (of pore volume) were selected for further use. Based on the rock physics analysis, the predicted lithology classes were associated with the geological rock-types (SRT) for incorporation in the static model. The study presents an innovative approach to successfully integrate geostatistical inversion and rock physics with static modeling. This workflow will generate seismically constrained high-resolution reservoir properties for thin reservoirs, such as porosity and lithology, which are seamlessly mapped in the depth domain for optimized development of the field. It will also account for the uncertainties in the reservoir model through the generation of multiple equiprobable realizations or scenarios.


2020 ◽  
Vol 39 (7) ◽  
pp. 480-487
Author(s):  
Patrick Smith ◽  
Brandon Mattox

The P-Cable high-resolution 3D marine acquisition system tows many short, closely separated streamers behind a small source. It can provide 3D seismic data of very high temporal and spatial resolution. Since the system is containerized and has small dimensions, it can be deployed at short notice and relatively low cost, making it attractive for time-lapse seismic reservoir monitoring. During acquisition of a 3D high-resolution survey in the Gulf of Mexico in 2014, a pair of sail lines were repeated to form a time-lapse seismic test. We processed these in 2019 to evaluate their geometric and seismic repeatability. Geometric repetition accuracy was excellent, with source repositioning errors below 10 m and bin-based receiver positioning errors below 6.25 m. Seismic data comparisons showed normalized root-mean-square difference values below 10% between 40 and 150 Hz. Refinements to the acquisition system since 2014 are expected to further improve repeatability of the low-frequency components. Residual energy on 4D difference seismic data was low, and timing stability was good. We conclude that the acquisition system is well suited to time-lapse seismic surveying in areas where the reservoir and time-lapse seismic signal can be adequately imaged by small-source, short-offset, low-fold data.


Sign in / Sign up

Export Citation Format

Share Document