A Nonlinear Automatic History Matching Technique for Reservoir Simulation Models

1972 ◽  
Vol 12 (06) ◽  
pp. 508-514 ◽  
Author(s):  
L. Kent Thomas ◽  
L.J. Hellums ◽  
G.M. Reheis

Abstract This paper presents a nonlinear optimization technique that automatically varies reservoir parameters to obtain a history match of held parameters to obtain a history match of held performance. The method is based on the classical performance. The method is based on the classical Gauss-Newton least-squares procedure. The range of each parameter is restricted by a box-type constraint and special provisions are included to handle highly nonlinear cases. Any combination of reservoir parameters may be used as the optimization variables and any set or sets of held data may be included in the match. Several history matches are presented, including examples from previous papers for comparison. In each of these examples, the technique presented here resulted in equivalent history matches in as few or fewer simulation runs. Introduction The history matching phase of reservoir simulations usually requires a trial-and-error procedure of adjusting various reservoir parameters procedure of adjusting various reservoir parameters and then calculating field performance. This procedure is continued until an acceptable match procedure is continued until an acceptable match between field and calculated performance has been obtained and can become quite tedious and time consuming, even with a small number of reservoir parameters, because of the interaction between the parameters, because of the interaction between the parameters and calculated performance. parameters and calculated performance. Recently various automatic or semiautomatic history-matching techniques have been introduced. Jacquard and Jain presented a technique based on a version of the method of steepest descent. They did not consider their method to be fully operational, however, due to the lack of experience with convergence. Jahns presented a method based on the Gauss-Newton equation with a stepwise solution for speeding convergence; but his procedure still required a large number of reservoir simulations to proceed to a solution. Coats et al. presented a proceed to a solution. Coats et al. presented a workable automatic history-matching procedure based on least-squares and linear programming. The method presented by Slater and Durrer is based on a gradient method and linear programming. In their paper they mention the difficulty of choosing a step paper they mention the difficulty of choosing a step size for their gradient method, especially for problems involving low values of porosity and problems involving low values of porosity and permeability. They also point out the need for a permeability. They also point out the need for a fairly small range on their reservoir description parameters for highly nonlinear problems. Thus, parameters for highly nonlinear problems. Thus, work in this area to date has resulted either in techniques based on a linear parameter-error dependence or in nonlinear techniques which require a considerable number of simulation runs. The method presented here is a nonlinear algorithm that will match both linear and nonlinear systems in a reasonable number of simulations. HISTORY MATCHING In a reservoir simulation, various performance data for the field, such as well pressures, gas-oil ratios, and water-oil ratios, are used as the basis for the match. During the matching of these performance data certain reservoir and fluid performance data certain reservoir and fluid parameters are assumed to be known while other parameters are assumed to be known while other less reliable data, forming the set (x1, x2...xn), are varied to achieve a match. The objective of the history-matching procedure presented in this paper is to minimize, in a presented in this paper is to minimize, in a least-squares sense, the error between the set of observed and calculated performance data, Fk(x1, x2... xn). SPEJ P. 508

2012 ◽  
Vol 518-523 ◽  
pp. 4376-4379
Author(s):  
Bao Yi Jiang ◽  
Zhi Ping Li

With the increase in computational capability, numerical reservoir simulation has become an essential tool for reservoir engineering. To minimize the objective function involved in the history matching procedure, we need to apply the optimization algorithms. This paper is based on the optimization algorithms used in automatic history matching.


2013 ◽  
Vol 748 ◽  
pp. 614-618
Author(s):  
Bao Yi Jiang ◽  
Zhi Ping Li ◽  
Cheng Wen Zhang ◽  
Xi Gang Wang

Numerical reservoir models are constructed from limited available static and dynamic data, and history matching is a process of changing model parameters to find a set of values that will yield a reservoir simulation prediction of data that matches the observed historical production data. To minimize the objective function involved in the history matching procedure, we need to apply the optimization algorithms. This paper is based on the optimization algorithms used in automatic history matching. Several optimization algorithms will be compared in this paper.


1983 ◽  
Vol 23 (05) ◽  
pp. 717-726 ◽  
Author(s):  
Bertrand A. See ◽  
Roland N. Horne

Abstract This paper addresses the optimization of field operations under a given set of technical and economic constraints and demonstrates that an optimal control procedure may be applied to any industrial oil and/or gas reservoir at reasonable cost and at acceptable accuracy level. Both primary and secondary recovery processes can be considered. The method proposed here is presented in two main sections. First, the modeling phase provides an approximately and locally linear model of the reservoir. A previously calibrated reservoir simulator model is used to perform a series of experiments. and a multiple variable regression analysis is used to fit the experimental data. The experimental design was one of the key issues in this work. Second, the optimization phase is performed with a linear programming algorithm. Nonlinear effects, such as performed with a linear programming algorithm. Nonlinear effects, such as those generated by the presence of gas, are approximated by several procedural iterations. procedural iterations. The application of this method to the case of a hypothetical reservoir demonstrates the validity of the optimal control procedure and shows convergence within an acceptable number of iterations. Introduction This investigation demonstrates the application of linear programming to a set of behavior equations derived from reservoir simulation results by use of a least-squares inversion procedure. The method is intended to optimize the production schedule of any reservoir for which the producer or injector well locations have already been fixed. The accuracy of this optimization procedure depends on many factors. the most important being the approximate procedure depends on many factors. the most important being the approximate linearization of the nonlinear system. Also important is reducing the required number of simulation runs until a satisfactory cost/accuracy compromise is obtained. It appears, that the reservoir engineer may contribute in reducing the experimental and calculation cost by properly selecting the series of simulation experiments. Simulator experiments, multivariable regression, least-squares inversion, the simplex algorithm. and validation are the major steps of this theoretical optimization procedure. The application of the process is demonstrated by working out a hypothetical practical example. The optimization of field operations has been explored by many authors, and many approaches have been suggested. Lee and Aronofsky established the first principles of this type of procedure by designing a time-discretized optimization process applied to a set of single-well reservoirs. Linear programming was used, and the objective function considered was the net profit. This first approach was improved by Aronofsky and Williams, who reduced the assumptions concerning the reservoir. Attra et al. refined the Lee and Aronofsky production model by introducing additional economic and technical factors, such as sales contract requirements or gas compressor limitations. For all these methods, the linear equations constituting the reservoir model were derived from material-balance considerations, and the reservoirs generally were assumed uniform and single-phase. SPEJ p. 717


1970 ◽  
Vol 10 (01) ◽  
pp. 66-74 ◽  
Author(s):  
K.H. Coats ◽  
J.R. Dempsey ◽  
J.H. Henderson

Abstract Reservoir description data largely determine the validity of simulated reservoir performance. This paper presents a method that employs the least paper presents a method that employs the least squares and linear programming techniques to determine a reservoir description from given performance data. The method bandies multiphase performance data. The method bandies multiphase as well as single-phase flow Problems. The description parameters determined by the method may be any physical properties that influence calculated field performance. We believe The technique offers considerably greater efficiency than previously reported techniques. Example applications presented include cases of single-phase gas flow, single-phase oil flow and two-phase gas-water flow. In these particular applications the method gave accurate results with a large range of uncertainty in the reservoir parameters, and with a small number of simulation parameters, and with a small number of simulation runs. Introduction The purpose of reservoir simulation is estimation of future reservoir performance under alternative well configurations or operating conditions. This estimation is increasingly being performed using rather complex, numerical reservoir models. Reservoir description data constitute the bulk of the required input data for these models, and the accuracy of these data largely determine the validity of the calculated results. Thus an obvious problem is the determination of an accurate problem is the determination of an accurate reservoir description. We treat the problem of determining a reservoir description that, when used as input data to a reservoir simulator, results in close agreement between calculated and observed field performance. Field history or performance data are presumed available for some period of time designated the "match period". The available field history may reflect single- or multiphase, multidimensional flow, and the performance data to be matched may be any mix of observed pressures, producing rates, gas-oil and/or water-oil producing ratios. The observed field performance may correspond to a period of depletion and/or injection, or to an period of depletion and/or injection, or to an interference test. Our method for determining a viable reservoir description requires a number of runs using a reservoir simulator, each run using a reservoir description that is random within limits specified by the engineer. We then use a second, small program, that utilizes least squares and linear program, that utilizes least squares and linear programming; techniques, to process the data output programming; techniques, to process the data output from those runs to determine a reservoir description. To illustrate and test this new method, we constructed three example reservoirs experiencing single-phase gas, single-phase oil and two-phase (gas-water) flow, respectively, in two spatial dimensions. Simulator runs were made using a given set of reservoir description parameters. The results of these runs were then treated as "data" and the description parameters considered unknown. The automatic history matching method described in this paper was applied to back out description parameter values from the performance "data". parameter values from the performance "data". The agreement between these values and the true parameter values is given below. parameter values is given below. Reed et al. present an actual field case where the manual approach to matching production history proved prohibitive in both man and machine time. proved prohibitive in both man and machine time. Our least squares, linear programming technique was then used to achieve a satisfactory and economical match of the reservoir performance data. SPEJ P. 66


Author(s):  
Po Ting Lin ◽  
Wei-Hao Lu ◽  
Shu-Ping Lin

In the past few years, researchers have begun to investigate the existence of arbitrary uncertainties in the design optimization problems. Most traditional reliability-based design optimization (RBDO) methods transform the design space to the standard normal space for reliability analysis but may not work well when the random variables are arbitrarily distributed. It is because that the transformation to the standard normal space cannot be determined or the distribution type is unknown. The methods of Ensemble of Gaussian-based Reliability Analyses (EoGRA) and Ensemble of Gradient-based Transformed Reliability Analyses (EGTRA) have been developed to estimate the joint probability density function using the ensemble of kernel functions. EoGRA performs a series of Gaussian-based kernel reliability analyses and merged them together to compute the reliability of the design point. EGTRA transforms the design space to the single-variate design space toward the constraint gradient, where the kernel reliability analyses become much less costly. In this paper, a series of comprehensive investigations were performed to study the similarities and differences between EoGRA and EGTRA. The results showed that EGTRA performs accurate and effective reliability analyses for both linear and nonlinear problems. When the constraints are highly nonlinear, EGTRA may have little problem but still can be effective in terms of starting from deterministic optimal points. On the other hands, the sensitivity analyses of EoGRA may be ineffective when the random distribution is completely inside the feasible space or infeasible space. However, EoGRA can find acceptable design points when starting from deterministic optimal points. Moreover, EoGRA is capable of delivering estimated failure probability of each constraint during the optimization processes, which may be convenient for some applications.


Sign in / Sign up

Export Citation Format

Share Document