Sensitivity Analysis of Two-Phase Reservoir History Matching

1989 ◽  
Vol 4 (03) ◽  
pp. 319-324 ◽  
Author(s):  
A.T. Watson
1984 ◽  
Vol 24 (06) ◽  
pp. 697-706 ◽  
Author(s):  
A.T. Watson ◽  
G.R. Gavalas ◽  
J.H. Seinfeld

Abstract Since the number of parameters to be estimated in a reservoir history match is potentially quite large, it is important to determine which parameters can be estimated with reasonable accuracy from the available data. This aspect can be called determining the identifiability of the parameters. The identifiability of porosity and absolute parameters. The identifiability of porosity and absolute and relative permeabilities on the basis of flow and pressure data in a two-phase (oil/water) reservoir is pressure data in a two-phase (oil/water) reservoir is considered. The question posed is: How accurately can one expect to estimate spatially variable porosity and absolute permeability and relative permeabilities given typical permeability and relative permeabilities given typical production and pressure data" To gain insight into this production and pressure data" To gain insight into this question, analytical solutions for pressure and saturation in a one-dimensional (1D) waterflood are used. The following, conclusions are obtained.Only the average value of the porosity can be determined on the basis of water/oil flow measurements.The permeability distribution can be determined from pressure drop data with an accuracy depending on the pressure drop data with an accuracy depending on the mobility ratio.Exponents in a power function representation of the relative permeabilities can he determined from WOR data alone but not nearly so accurately as when pressure drop and flow data are used simultaneously. Introduction The utility of reservoir simulation in predicting reservoir behavior is limited by the accuracy with which reservoir properties can be estimated. Because of the high costs properties can be estimated. Because of the high costs associated with coring analysis, reservoir engineers must rely, on history matching as a means of estimating reservoir properties. In this process a history match is carried out by choosing the reservoir properties as those that result in simulated well pressure and flow data that match as closely as possible those measured during production. In general, reservoir properties at each gridblock in the simulator represent the unknown values to be determined. Although there are efficient methods for estimating such a large number of unknowns, it has long been recognized from the results of single phase history matching exercises that many different sets of parameter values may yield a nearly identical match of observed and predicted pressures. The conventional single phase predicted pressures. The conventional single phase history matching problem is in fact a mathematically illposed problem, which explains its nonunique behavior. Such a situation is, in short, the result of the large number of unknowns to be estimated on the basis of the available data and the lack of sensitivity of the simulator solutions to the parameters. Because of this lack of sensitivity, the need to reduce the number of unknown Parameters or to introduce some additional constraints, such as "smoothness" of the estimated parameters, has been recognized. A problem as important as that of choosing which minimization method to employ in history matching is that of choosing, on the basis of the available well data. which properties actually should be estimated. This selection properties actually should be estimated. This selection depends on the relationship of the unknown parameters to the simulated well data. Ideally one would want to knowwhich parameters can be determined uniquely if the measurements were exact, andgiven the expected level of error in the measurements, how accurately we can expect to be able to estimate the parameters. The first question, that of establishing uniqueness of the estimated parameters, is notoriously difficult to answer, and for a parameters, is notoriously difficult to answer, and for a problem as complicated as reservoir history matching, problem as complicated as reservoir history matching, there are virtually no general results available that allow one to establish uniqueness for permeability or porosity. Thus, it is not possible in general to base our choice of which parameters to estimate on rigorous mathematical uniqueness results. In lieu of an answer to Question 1, the selection of parameters to be estimated can be based on Question 2, parameters to be estimated can be based on Question 2, which is amenable to theoretical analysis. If the expected errors in estimation of any of the parameters, or any linear combination of the parameters, are extremely large, then that parameter or set of parameters can be judged as not identifiable. In such a case, steps may be taken to reduce the number of unknown parameters. In summary, the reservoir history matching problem is a difficult parameter estimation problem, and understanding the relationship between the unknown parameters and the measured data is essential to obtaining meaningful estimates of the reservoir properties. Quantitative studies regarding the accuracy of estimates for single-phase history matching problems have been reported by Shah et al. and Dogru et al. Shah et al,. investigated the optimal level of zonation for use with 1D single-phase (oil) situations. SPEJ P. 697


Energies ◽  
2020 ◽  
Vol 13 (17) ◽  
pp. 4290
Author(s):  
Dongmei Zhang ◽  
Yuyang Zhang ◽  
Bohou Jiang ◽  
Xinwei Jiang ◽  
Zhijiang Kang

Reservoir history matching is a well-known inverse problem for production prediction where enormous uncertain reservoir parameters of a reservoir numerical model are optimized by minimizing the misfit between the simulated and history production data. Gaussian Process (GP) has shown promising performance for assisted history matching due to the efficient nonparametric and nonlinear model with few model parameters to be tuned automatically. Recently introduced Gaussian Processes proxy models and Variogram Analysis of Response Surface-based sensitivity analysis (GP-VARS) uses forward and inverse Gaussian Processes (GP) based proxy models with the VARS-based sensitivity analysis to optimize the high-dimensional reservoir parameters. However, the inverse GP solution (GPIS) in GP-VARS are unsatisfactory especially for enormous reservoir parameters where the mapping from low-dimensional misfits to high-dimensional uncertain reservoir parameters could be poorly modeled by GP. To improve the performance of GP-VARS, in this paper we propose the Gaussian Processes proxy models with Latent Variable Models and VARS-based sensitivity analysis (GPLVM-VARS) where Gaussian Processes Latent Variable Model (GPLVM)-based inverse solution (GPLVMIS) instead of GP-based GPIS is provided with the inputs and outputs of GPIS reversed. The experimental results demonstrate the effectiveness of the proposed GPLVM-VARS in terms of accuracy and complexity. The source code of the proposed GPLVM-VARS is available at https://github.com/XinweiJiang/GPLVM-VARS.


SPE Journal ◽  
2017 ◽  
Vol 22 (05) ◽  
pp. 1506-1518 ◽  
Author(s):  
Pedram Mahzari ◽  
Mehran Sohrabi

Summary Three-phase flow in porous media during water-alternating-gas (WAG) injections and the associated cycle-dependent hysteresis have been subject of studies experimentally and theoretically. In spite of attempts to develop models and simulation methods for WAG injections and three-phase flow, current lack of a solid approach to handle hysteresis effects in simulating WAG-injection scenarios has resulted in misinterpretations of simulation outcomes in laboratory and field scales. In this work, by use of our improved methodology, the first cycle of the WAG experiments (first waterflood and the subsequent gasflood) was history matched to estimate the two-phase krs (oil/water and gas/oil). For subsequent cycles, pertinent parameters of the WAG hysteresis model are included in the automatic-history-matching process to reproduce all WAG cycles together. The results indicate that history matching the whole WAG experiment would lead to a significantly improved simulation outcome, which highlights the importance of two elements in evaluating WAG experiments: inclusion of the full WAG experiments in history matching and use of a more-representative set of two-phase krs, which was originated from our new methodology to estimate two-phase krs from the first cycle of a WAG experiment. Because WAG-related parameters should be able to model any three-phase flow irrespective of WAG scenarios, in another exercise, the tuned parameters obtained from a WAG experiment (starting with water) were used in a similar coreflood test (WAG starting with gas) to assess predictive capability for simulating three-phase flow in porous media. After identifying shortcomings of existing models, an improved methodology was used to history match multiple coreflood experiments simultaneously to estimate parameters that can reasonably capture processes taking place in WAG at different scenarios—that is, starting with water or gas. The comprehensive simulation study performed here would shed some light on a consolidated methodology to estimate saturation functions that can simulate WAG injections at different scenarios.


1980 ◽  
Vol 20 (06) ◽  
pp. 521-532 ◽  
Author(s):  
A.T. Watson ◽  
J.H. Seinfeld ◽  
G.R. Gavalas ◽  
P.T. Woo

Abstract An automatic history-matching algorithm based onan optimal control approach has been formulated forjoint estimation of spatially varying permeability andporosity and coefficients of relative permeabilityfunctions in two-phase reservoirs. The algorithm usespressure and production rate data simultaneously. The performance of the algorithm for thewaterflooding of one- and two-dimensional hypotheticalreservoirs is examined, and properties associatedwith the parameter estimation problem are discussed. Introduction There has been considerable interest in thedevelopment of automatic history-matchingalgorithms. Most of the published work to date onautomatic history matching has been devoted tosingle-phase reservoirs in which the unknownparameters to be estimated are often the reservoirporosity (or storage) and absolute permeability (ortransmissibility). In the single-phase problem, theobjective function usually consists of the deviationsbetween the predicted and measured reservoirpressures at the wells. Parameter estimation, orhistory matching, in multiphase reservoirs isfundamentally more difficult than in single-phasereservoirs. The multiphase equations are nonlinear, and in addition to the porosity and absolutepermeability, the relative permeabilities of each phasemay be unknown and subject to estimation. Measurements of the relative rates of flow of oil, water, and gas at the wells also may be available forthe objective function. The aspect of the reservoir history-matchingproblem that distinguishes it from other parameterestimation problems in science and engineering is thelarge dimensionality of both the system state and theunknown parameters. As a result of this largedimensionality, computational efficiency becomes aprime consideration in the implementation of anautomatic history-matching method. In all parameterestimation methods, a trade-off exists between theamount of computation performed per iteration andthe speed of convergence of the method. Animportant saving in computing time was realized insingle-phase automatic history matching through theintroduction of optimal control theory as a methodfor calculating the gradient of the objective functionwith respect to the unknown parameters. Thistechnique currently is limited to first-order gradientmethods. First-order gradient methods generallyconverge more slowly than those of higher order.Nevertheless, the amount of computation requiredper iteration is significantly less than that requiredfor higher-order optimization methods; thus, first-order methods are attractive for automatic historymatching. The optimal control algorithm forautomatic history matching has been shown toproduce excellent results when applied to field problems. Therefore, the first approach to thedevelopment of a general automatic history-matchingalgorithm for multiphase reservoirs wouldseem to proceed through the development of anoptimal control approach for calculating the gradientof the objective function with respect to theparameters for use in a first-order method. SPEJ P. 521^


2021 ◽  
Author(s):  
Omar Chaabi ◽  
Emad W. Al-Shalabi ◽  
Waleed Alameri

Abstract Low salinity polymer (LSP) flooding is getting more attention due to its potential of enhancing both displacement and sweep efficiencies. Modeling LSP flooding is challenging due to the complicated physical processes and the sensitivity of polymers to brine salinity. In this study, a coupled numerical model has been implemented to allow investigating the polymer-brine-rock geochemical interactions associated with LSP flooding along with the flow dynamics. MRST was coupled with the geochemical software IPhreeqc. The effects of polymer were captured by considering Todd-Longstaff mixing model, inaccessible pore volume, permeability reduction, polymer adsorption as well as salinity and shear rate effects on polymer viscosity. Regarding geochemistry, the presence of polymer in the aqueous phase was considered by adding a new solution specie and related chemical reactions to PHREEQC database files. Thus, allowing for modeling the geochemical interactions related to the presence of polymer. Coupling the two simulators was successfully performed, verified, and validated through several case studies. The coupled MRST-IPhreeqc simulator allows for modeling a wide variety of geochemical reactions including aqueous, mineral precipitation/dissolution, and ion exchange reactions. Capturing these reactions allows for real time tracking of the aqueous phase salinity and its effect on polymer rheological properties. The coupled simulator was verified against PHREEQC for a realistic reactive transport scenario. Furthermore, the coupled simulator was validated through history matching a single-phase LSP coreflood from the literature. This paper provides an insight into the geochemical interactions between partially hydrolyzed polyacrylamide (HPAM) and aqueous solution chemistry (salinity and hardness), and their related effect on polymer viscosity. This work is also considered as a base for future two-phase polymer solution and oil interactions, and their related effect on oil recovery.


Sign in / Sign up

Export Citation Format

Share Document