Integration of time-lapse seismic data into dynamic reservoir model is an efficient process in calibrating reservoir parameters update. The choice of the metric which will measure the misfit between observed data and simulated model has a considerable effect on the history matching process, and then on the optimal ensemble model acquired. History matching using 4D seismic and production data simultaneously is still a challenge due to the nature of the two different type of data (time-series and maps or volumes based).
Conventionally, the formulation used for the misfit is least square, which is widely used for production data matching. Distance measurement based objective functions designed for 4D image comparison have been explored in recent years and has been proven to be reliable. This study explores history matching process by introducing a merged objective function, between the production and the 4D seismic data. The proposed approach in this paper is to make comparable this two type of data (well and seismic) in a unique objective function, which will be optimised, avoiding by then the question of weights. An adaptive evolutionary optimisation algorithm has been used for the history matching loop. Local and global reservoir parameters are perturbed in this process, which include porosity, permeability, net-to-gross, and fault transmissibility.
This production and seismic history matching has been applied on a UKCS field, it shows that a acceptalbe production data matching is achieved while honouring saturation information obtained from 4D seismic surveys.
An ensemble-based history-matching framework is proposed to enhance the characterization of petroleum reservoirs through the assimilation of crosswell electromagnetic (EM) data. As one of advanced technologies in reservoir surveillance, crosswell EM tomography can provide a cross-sectional conductivity map and hence saturation profile at an interwell scale by exploiting the sharp contrast in conductivity between hydrocarbons and saline water. Incorporating this new information into reservoir simulation in combination with other available observations is therefore expected to enhance the forecasting capability of reservoir models and to lead to better quantification of uncertainty.
The proposed approach applies ensemble-based data-assimilation methods to build a robust and flexible framework under which various sources of available measurements can be readily integrated. Because the assimilation of crosswell EM data can be implemented in different ways (e.g., components of EM fields or inverted conductivity), a comparative study is conducted. The first approach integrates crosswell EM data in its original form which entails establishing a forward model simulating observed EM responses. In this work, the forward model is based on Archie's law that provides a link between fluid properties and formation conductivity, and Maxwell’s equations that describe how EM fields behave given the spatial distribution of conductivity. Alternatively, formation conductivity can be used for history matching, which is obtained from the original EM data through inversion using an adjoint gradient-based optimization method. Because the inverted conductivity is usually of high dimension and very noisy, an image-oriented distance parameterization utilizing fluid front information is applied aiming to assimilate the conductivity field efficiently and robustly. Numerical experiments for different test cases with increasing complexity are carried out to examine the performance of the proposed integration schemes and potential of crosswell EM data for improving the estimation of relevant model parameters. The results demonstrate the efficiency of the developed history-matching workflow and added value of crosswell EM data in enhancing the characterization of reservoir models and reliability of model forecasts.
We present a novel sampling algorithm for characterization and uncertainty quantification of heterogeneous multiple facies reservoirs. The method implements a Bayesian inversion framework to estimate physically plausible porosity distributions. This inversion process incorporates data matching at the well locations and constrains the model space by adding
The proposed workflow uses an ensemble-based Markov Chain Monte Carlo approach combined with sampling probability distributions that are physically meaningful. Moreover, the method targets geostatistical modeling to specific zones in the reservoir. Accordingly, it improves fulfilling the inherent stationarity assumption in geostatistical simulation techniques. Parameter sampling and geostatistical simulations are calculated through an inversion process. In other words, the models fit the known porosity field at the well locations and are structurally consistent within main reservoir compartments, zones, and layers obtained from the seismic impedance volume. The new sampling algorithm ensures that the automated history matching algorithm maintains diversity among ensemble members avoiding underestimation of the uncertainty in the posterior probability distribution.
We evaluate the efficiency of the sampling methodology on a synthetic model of a waterflooding field. The predictive capability of the assimilated ensemble is assessed by using production data and dynamic measurements. Also, the qualities of the results are examined by comparing the geological realism of the assimilated ensemble with the reference probability distribution of the model parameters and computing the predicted dynamic data mismatch. Our numerical examples show that incorporating the seismically constrained models as prior information results in an efficient model update scheme and favorable history matching.
Chen, Hongquan (Texas A&M University) | Yang, Changdong (Texas A&M University) | Datta-Gupta, Akhil (Texas A&M University) | Zhang, Jianye (Institution of Exploration and Development of Tarim Oilfield Company-Petro China) | Chen, Liqun (Institution of Exploration and Development of Tarim Oilfield Company-Petro China) | Liu, Lei (Institution of Exploration and Development of Tarim Oilfield Company-Petro China) | Chen, Baoxin (Institution of Exploration and Development of Tarim Oilfield Company-Petro China) | Cui, Xiaofei (Optimization Petroleum Technology, Inc.) | Shi, Fashun (Optimization Petroleum Technology, Inc.) | Bahar, Asnul (Kelkar and Associates, Inc.)
History matching of million-cell reservoir models still remains an outstanding challenge for the industry. This paper presents a hierarchical multi-scale approach to history matching high resolution dual porosity reservoir models using a combination of evolutionary algorithm and streamline method. The efficacy of the approach is demonstrated through application to a high pressure high temperature (HPHT) fractured gas reservoir in the Tarim basin, China with wells located at an average depth of 7500 meters.
Our proposed multi-scale history matching approach consists of two-stages: global and local. For the global stage, we calibrate coarse-scale static and dynamic parameters using an evolutionary algorithm. The global calibration uses coarse-scale simulations and applies regional multipliers to match RFT data, well bottom hole pressures, and field average pressure. For the local stage, we calibrate fracture permeability using streamline based sensitivities to further match well bottom-hole pressures. The streamlines are derived from the fracture cell fluxes and the sensitivities are analytically computed for highly compressible flow. The sensitivities are validated by comparison with the pertubation method.
The proposed hierarchical multiscale history matching workflow is applied to a faulted and highly fractured deep gas reservoir in the Tarim basin, China. The excessive well cost arising from the large well depth (7500 meters) and high pressure (18000 psi) necessitates optimal field development with limited number of wells. The fracture properties of dual porosity model are upscaled from a highly dense discrete fracture network model generated based on well data and seismic attributes. The history matching includes RFT data, static pressure data and flowing bottom-hole pressure data in producing wells. Field average pressure and RFT (static pressure) data were well matched during the global stage using coarse scale models while flowing bottom-hole pressure is further matched during the local stage calibration using fine scale models. Streamline method has been applied previously mainly to incompressible or slightly compressible flow. However in this application, the results show that the modified streamline-based sensitivity can also significantly reduce data misfit for highly compressible flow. The history matched models are used to visualize well drainage volumes using streamlines. The well drainage volumes in conjunction with static reservoir properties are used to define a ‘depletion capacity map’ which is then used for optimal infill well placement.
The novelty of our approach lies in the application of streamlines derived from dual porosity finite-difference simulation to facilitate history matching and well placement optimization in a tight gas reservoir. The newly developed streamline-based analytical sensitivities are suitable for highly compressible flow. To our knowledge, this is the first time streamlines have been used to facilitate history matching and optimal well placement for gas reservoirs.
Standard history matching workflows use qualitative 4D seismic observations to assist in reservoir modeling and simulation. However, such workflows lack a robust framework for quantitatively integrating 4D seismic interpretations. 4D or time-lapse seismic interpretations provide valuable inter-well saturation and pressure information and quantitatively integrating this inter-well data can help to constrain simulation parameters and improve the reliability of production modeling. This paper outlines technologies aimed at leveraging the value of 4D for reducing uncertainty in the range of history matched models and improving the production forecast.
The proposed 4D Assisted History Match (4DAHM) workflows utilize interpretations of 4D seismic anomalies for improving the reservoir simulation models. Design of Experiments (DOE) is initially used to generate the probabilistic history match simulations by varying the range of uncertain parameters. Saturation maps are extracted from the Production History Matched (PHM) simulations and then compared with 4D predicted swept anomalies. An automated extraction method was created and is used to reconcile spatial sampling differences between 4D data and simulation output. Interpreted 4D data is compared with simulation output, and the mismatch generated is used as a 4D filter to refine the suite of reservoir simulation models. The selected models are used to identify reservoir simulation parameters that are sensitive for generating a good match.
The application of 4DAHM workflows has resulted in reduced uncertainty in volumetric predictions of oil fields, probabilistic saturation S-curves at target locations, and fundamental changes to the dynamic model needed to improve the match to production data. Results from adopting this workflow in two different deep-water reservoirs are discussed. They not only resulted in reduced uncertainty, but also provided information on key performance indicators that are critical in obtaining a robust history match. In the first case study presented, the deep-water oil field 4DAHM resulted in a reduction of uncertainty by 20% in OOIP and by 25% in EUR in the P90-P10 range estimates. In the second case study, 4DAHM workflow exploited discrepancies between 4D seismic and simulation data to identify features necessary to be included in the dynamic model. Connectivity was increased through newly interpreted inter-channel erosional contacts, and sub-seismic faults. Moreover, the workflow provided an improved drilling location which has the higher probability of tapping unswept oil and better EUR. The 4D filters constrained the suite of reservoir simulation models and helped to identify 4 out of 24 simulation parameters critical for success. The updated PHM models honor both the production data and 4D interpretations, resulting in reduced uncertainty across the S-curve and, in this case, an increased P50 OOIP of 24% for a proposed infill drilling location, plus a significant cycle-time savings.
Several authors have investigated quantitative approaches for incorporating a large number of seismic data into history-matching work flows for production data. The work flows adopted in most of these studies require significant reduction of the uncertainty space or produce on a single history-matched model. In recent years, an algorithm that has gained increased popularity is the ensemble Kalman filter (EnKF). The EnKF and its derived algorithms are called ensemble methods, and their most notable characteristic is being computationally feasible for large systems and being relatively easy to implement. Several authors have investigated the effect of assimilating 3D- or 4D-seismic data with EnKF and have had to face the problem of building an accurate petroelastic model. To circumvent this problem, some authors have reparameterized 4D-seismic anomaly-front data into arrival times. Despite the advantage of eliminating seismic inversion, this method presents the disadvantage of at least doubling the simulation time for accurate ensemble arrival-time predictions. This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 183901, “Ensemble-Based Assisted History Matching With 4D-Seismic Fluid-Front Parameterization,” by Mario Trani, Konrad Wojnar, Arthur Moncorgé, and Philippe Berthet, Total, prepared for the 2017 SPE Middle East Oil and Gas Show and Conference, Manama, Bahrain, 6–9 March. The paper has not been peer reviewed.
Khani, Hojjat (University of Calgary) | Hamdi, Hamidreza (University of Calgary) | Nghiem, Long (Computer Modelling Group Ltd.) | Chen, Zhangxing (University of Calgary) | Sousa, Mario Costa (University of Calgary)
The overall objective of reservoir modeling is to reduce the uncertainty in production forecasts by utilizing all available data to construct a calibrated reservoir model. Geological heterogeneities have a fundamental impact on the growth of a steam chamber and the performance of a SAGD (steam assisted gravity drainage) process. The objective of this work is to incorporate geological heterogeneities into the history matching process using a probability perturbation method (PPM) to preserve the geological consistency of a reservoir model.
A PPM is a geological data integration framework which employs a multiple-point geostatistics (MPS) algorithm. The heart of this method is to systematically perturb the underlying probabilities used to generate the reservoir facies. A PPM generally consists of two loops: an outer loop which is responsible for randomly generating a global configuration of the facies and an inner loop which systematically perturbs the generated facies to match the dynamic data. The combination of these two iterations creates a set of realizations that preserve the geological information.
In this paper, a training image is built based on a 3D outcrop description of a meandering channelized reservoir that is analogous to some of the Canadian heavy oil reservoirs. All other available data including reservoir properties at well locations, trends and production data are also incorporated into the PPM framework for this history matching process. The reservoir model is characterized by three facies: clean sands, medium-grained sandstones and silts, which have different porosity, horizontal permeability and vertical permeability. The SAGD performance is a function of steam chamber development, which depends on the level of heterogeneity in the reservoir. The results show that the heterogeneity distribution has a large impact on the fluid flow at different stages of production. The results show that such complexities can be well preserved during the history matching process using the PPM by generating the geological patterns depicted in a training image. The PPM is shown to be an efficient approach for history matching in presence of complex reservoir geology.
Nguyen, Ngoc T. B. (University of Calgary) | Dang, Cuong T. Q. (Computer Modeling Group Ltd.) | Yang, Chaodong (Computer Modeling Group Ltd.) | Nghiem, Long X. (Computer Modeling Group Ltd.) | Chen, Zhangxin (University of Calgary)
Steam Assisted Gravity Drainage (SAGD) has been widely applied to unlock hydrocarbon resources in oil sands reservoirs. This method uses steam, which is generated at the surface, to heat a formation and create a steam chamber around an injector. Past studies have indicated that reservoir heterogeneity is one of the crucial factors that directly affectthe performance of the SAGD process. This paper presents an innovative integrated modeling approach for evaluating, assisted history matching, and production forecasting of the SAGD process with the presence of a complex shale barrier system in oil sands reservoirs.
As SAGD is a strongly geological dependent recovery process and, unfortunately, there are many uncertainties associated with reservoir geology in reality. Therefore, it requires generating a large number of geological realizations to capture the critical effects of geology, especially with the presence of shale barriers, in history matching and field development planning of the SAGD process. To quantify the impact and improve the quality of history matching compared with the traditional method, an efficient integrated workflow has been developed in which geological information generated from a geological modeling package is automatically updated for a reservoir simulator and controlled by an intelligent optimizer in a big-loop modeling approach.
A detailed workflow on the integrated modeling approach that includes shale barriers for a typical oil sands reservoir is described in the first section of this paper. Shale bodies are geostatistically distributed in the geological models. A comprehensive parametric study was conducted with numerous geological realizations to identify the critical role of shale barriers in SAGD performance including shale geometry, shale length and thickness, shale distribution and proportions. Then the Bayesian algorithm with a Proxy-based Acceptance-Rejection sampling method is employed for assisted history matching of SAGD production profiles. With the presence of complex shale barriers, it requires simultaneous updating of both geological and reservoir engineering parameters. Using the proposed approach, the global history matching errors were drastically reduced in all production wells. Validation results indicate that the integrated modeling approach effectively helps to update the properties and distribution of shale barriers to find the closest geological distribution compared to the true solution. Finally, an ensemble of the best-matched simulation models is used to perform a probabilistic forecasting to capture the uncertainties in future production profiles.
Not limited to history matching ofthe SAGD process, the proposed approach can be also applied to different complex problems such as robust optimization for various recovery methods from conventional to unconventional reservoirs.
In this work, a Bayesian data assimilation methodology for simultaneous estimation of channelized facies and petrophysical properties (e.g., permeability fields) is explored. Based on the work of
This work focuses on the improvement of an integrated methodology for the automatic history matching of compartmentalised reservoirs using 4D seismic results, stochastic initialization and the Ensemble Kalman Filter method. We show the comparison of two different history matching approaches using the Ensemble Kalman Filter (EnKF) to update the Fault Transmissibility Multipliers (FTM) initially estimated with and without considering the 4D seismic results. In this study, the parameters updated during the history matching are two-phase fault transmissibility multipliers (FTM), absolute permeability and effective porosity of a synthetic realistic 3D reservoir. The true impedance map and the changes in reservoir pressure and saturation were previously computed from 4D seismic results. The systematic estimation of two-phase fault transmissibility multipliers is based on the integration of the collected 4D seismic results and an established method validated in our previous work based on a deterministic model, using the gradient-based History Matching, Levenberg Marquardt method (LM). We present the history matching of a synthetic reservoir using the Ensemble Kalman Filter (EnKF) considering the 4D seismic results to update the models and geostatistical techniques to produce the initial geological models. The stochastic method used is the Sequential Gaussian Simulation (SGS) technique to generate 100 initial models. During history matching using the EnKF, the saturation distributions are computed from the forward modelling of a two-phase system (oil-water). The impedance maps are then estimated using the Gassmann equation and compared with the true impedance map as part of the History Matching process. To validate the results, the cost function consisting of two components is calculated, the first is the structural similarity index of the two reconstructed impedance images to the real impedance image and the second is the RMS cost function value,