Adjoint based optimal control of partially miscible two-phase flow in porous media with applications to CO2 sequestration in underground reservoirs
KAUST Grant NumberUK-C0020
Online Publication Date2014-11-14
Print Publication Date2015-03
Permanent link to this recordhttp://hdl.handle.net/10754/597465
MetadataShow full item record
Abstract© 2014, Springer Science+Business Media New York. With the target of optimizing CO2 sequestration in underground reservoirs, we investigate constrained optimal control problems with partially miscible two-phase flow in porous media. Our objective is to maximize the amount of trapped CO2 in an underground reservoir after a fixed period of CO2 injection, while time-dependent injection rates in multiple wells are used as control parameters. We describe the governing two-phase two-component Darcy flow PDE system, formulate the optimal control problem and derive the continuous adjoint equations. For the discretization we apply a variant of the so-called BOX method, a locally conservative control-volume FE method that we further stabilize by a periodic averaging feature to reduce oscillations. The timestep-wise Lagrange function of the control problem is implemented as a variational form in Sundance, a toolbox for rapid development of parallel FE simulations, which is part of the HPC software Trilinos. We discuss the BOX method and our implementation in Sundance. The MPI parallelized Sundance state and adjoint solvers are linked to the interior point optimization package IPOPT, using limited-memory BFGS updates for approximating second derivatives. Finally, we present and discuss different types of optimal control results.
CitationSimon M, Ulbrich M (2014) Adjoint based optimal control of partially miscible two-phase flow in porous media with applications to CO2 sequestration in underground reservoirs. Optim Eng 16: 103–130. Available: http://dx.doi.org/10.1007/s11081-014-9270-x.
SponsorsThe support from Award No. UK-C0020, made by King Abdullah University of Science and Technology (KAUST) is gratefully acknowledged. This work was conducted as part of the MAC-KAUST project K1 "Simulating CO<INF>2</INF> Sequestration" within the Munich Centre of Advanced Computing (MAC) at TUM. The computations were performed on a compute cluster that was partially funded by DFG INST 95/919-1 FUGG. Finally, we thank the referees for their valuable suggestions that helped us to improve the quality of the article.
JournalOptimization and Engineering