Sensitivity Analysis of Water Flooding Optimization by Dynamic Optimization

Document Type: Research Article


Research Institute of Petroleum Industry (RIPI), P.O. Box 18745-4163 Tehran, I.R. IRAN


This study concerns the scope to improve water flooding in heterogeneous reservoirs. We used an existing, in-house developed, optimization program consisting of a reservoir simulator in combination with an adjoint-based optimal control algorithm. In particular we aimed to examine the scope for optimization in a two-dimensional horizontal reservoir containing a single high permeable streak, as a function of reservoir and fluid parameters, which we combined in the form of 10 dimensionless parameters. We defined the parameter NPVimprovement to indicate the improvement in net present value (NPV) that can be achieved through optimization. For initial screening of the effect of the dimensionless parameters, a two-level D-optimal design of experiments (DOE) technique was used to obtain a linear response surface model with the aid of 11 water-flooding simulations. As a result 8 dimensionless groups were selected for more detailed analysis, and a full quadratic NPVimprovement model was constructed using a three-level D-optimal design using 50 simulations. It should be reminded that all the D-optimal matrix designs were generated by using commands of statistics Toolbox of MATLAB software. Finally, Pareto charts were plotted to visualize the sensitivity of the model as a function of the dimensionless parameters. Based on the present model we can draw the conclusion that the parameters (relative streak length),  (relative streak permeability) and the ratio of water cost and oil price have the largest effect on the scope for obtaining a high value of NPVimprovement.


Main Subjects

[1] Asheim, H., Maximization of Water Sweep Effi-ciency by Controlling Production and Injection Rates, paper SPE 18365,(1991).

[2] Brouwer, D.,R., Jansen, J.,D., van der Starre, S., and Van Kruijsdijk, C.P.J.W., Recovery Increase through Water Flooding using Smart Well Technology, paper SPE 68979, (2000).

[3] Brouwer, D., R., and Jansen, J.D., Dynamic Optimization of Water Flooding with Smart Wells Using Optimal Control, paper SPE 78278, (2001).

[4] Sudaryanto, B., and Yortsos, Y., C., Optimization of Fluid Front Dynamics in Porous Media Using Rate Control,. I. Equal Mobility Fluids, Phys. of Fluids, 12 (7), 1656, (2000).

[5] Yeten, B., Durlofsky, L.J. and Aziz, K., Optimization of Smart Well Control, paper SPE 79031, (2001).

[6] Box, G.E.P., Hunter, W.G. and Hunter, J.S., “Statistics for Experimenters”, John Wiley & Sons, (1978).

[7] Friedmann, F., Chawathe, A. and Larue, D., Assessing Uncertainty in Channelized Reservoirs using Experimental Designs, paper SPE 71622,(2000).

[8] Atkinson, A. C. and Donev A. N., Optimum Experi-mental Designs, Oxford University Press, Oxford, UK, (1992).

[9] Cook, R.D. and Nachtsheim, C.J., A Comparison of Algorithms for Constructing exact D-Optimal Designs, Technometrics, 22 (3), August (1980).

[10] Narayanan,  K.,  Applications for Response Surfaces in Reservoir Engineering, MS Thesis, University of Texas, Austin, (1999).

[11] Draper, N. R. and Lin, D. K. J., Small Response-Surface Designs, Technometrics, 32,187 (1990).

[12] Box, G.E.P. and Draper, N.R., “Empirical Model-Building and Response Surface”, John Wiley and Sons, (1987).

[13] Montgomery, D. C. and Peck, E. A., “Introduction to Linear Regression Analysis”, John Wiley and Sons, (1982).