Production optimization of thermodynamically rigorous isothermal and compositional models

Research output: Chapter in Book/Report/Conference proceedingArticle in proceedingsResearchpeer-review


In this work, we consider algorithms for solving production optimization problems that involve isothermal (constant temperature) and compositional oil production processes. The purpose of production optimization is to compute a long-term production strategy that is economically optimal. We present a thermodynamically rigorous model of isothermal oil production processes. We derive the model from first principles by applying a number of assumptions including the assumption of constant temperature. The model is based on two key principles, namely phase equilibrium and conservation of mass and energy. The conservation equations are expressed as partial differential equations, and we model the phase equilibrium as a VT flash process. It is common to formulate the phase equilibrium conditions in oil reservoir flow models as the fugacities being equal. We describe how to derive that condition from the phase equilibrium conditions from the VT flash problem. The VT flash is an adaption of the second law of thermodynamics, i.e. the entropy of a closed system in equilibrium is maximal, to isothermal systems. The VT flash can therefore be formulated as an inner optimization problem that needs to be solved for each grid cell in the discretized reservoir in the forward simulation of the oil production process. We demonstrate that it is natural to model such isothermal production processes with differential-algebraic equations in a semi-explicit index-1 form. We describe a single-shooting algorithm for solving the production optimization problem efficiently. It is key to the efficiency of such algorithms to compute gradients. For that purpose, we use an adjoint algorithm. We implement the single-shooting algorithm in C/C++ using the open-source software DUNE, the open-source thermodynamic software ThermoLib, and the numerical optimization software KNITRO. Finally, we present a numerical example that involves optimal waterflooding.
Original languageEnglish
Title of host publicationProceedings of 16th European Conference on the Mathematics of Oil Recovery
PublisherEuropean Association of Geoscientists and Engineers
Publication date2018
ISBN (Print)9789462822603
Publication statusPublished - 2018
Event16th European Conference on the Mathematics of Oil Recovery - World Trade Center Barcelona, Barcelona, Spain
Duration: 3 Sept 20186 Sept 2018
Conference number: 16


Conference16th European Conference on the Mathematics of Oil Recovery
LocationWorld Trade Center Barcelona
Internet address


Dive into the research topics of 'Production optimization of thermodynamically rigorous isothermal and compositional models'. Together they form a unique fingerprint.

Cite this