Adjoint based optimal control of partially miscible two-phase flow in porous media with applications to CO2 sequestration in underground reservoirs
Type
ArticleAuthors
Simon, MoritzUlbrich, Michael
KAUST Grant Number
UK-C0020Date
2014-11-14Online Publication Date
2014-11-14Print Publication Date
2015-03Permanent link to this record
http://hdl.handle.net/10754/597465
Metadata
Show full item recordAbstract
© 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.Citation
Simon 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.Sponsors
The 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.Publisher
Springer NatureJournal
Optimization and Engineeringae974a485f413a2113503eed53cd6c53
10.1007/s11081-014-9270-x