A mathematical programming model for the optimal utilization of deep saline aquifers for CO2 storage

• Published in: Computers & chemical engineering Vol. 203; p. 109343
• Main Authors: Trucco, Diego J., Presser, Demian J., Cafaro, Diego C., Grossmann, Ignacio E., Usgaonkar, Saurabh Shenvi, Zhang, Qi, Misra, Pratik, Binagia, Heather, Rowe, Wayne, Mehta, Sanjay
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.12.2025
• Subjects: CCUS; Containment; Injection strategy; Injectivity; Optimization; Saline aquifers; Injection strategy; CCUS; Injectivity; Optimization; Saline aquifers; Containment
• ISSN: 0098-1354
• DOI: 10.1016/j.compchemeng.2025.109343

•Technically sound, safe, and efficient storage of CO2 in deep saline aquifers.•Maximization of CO2 storage subject to injectivity and containment constraints.•Application of mathematical programming, a novel approach to optimize CO2 injection.•Model and results validated through reservoir engineering and geological simulation. This work presents a novel nonlinear programming (NLP) formulation aimed at maximizing the overall amount of CO2 stored into deep saline aquifers in the long term. The goal is to optimally determine CO2 injection rates into vertical wells while properly managing bottom-hole pressures over time. The reservoir may comprise several layers with heterogeneous physical properties. The injection plan should meet the subsurface engineering policies for safe operations along with existing technical constraints. The major challenge is to track the CO2 migration across the reservoir to ensure containment during the injection periods and also in the long term. The NLP formulation is based on a discrete space and time representation of the reservoir, comprising pressure propagation and mass balance equations between every pair of adjacent blocks in the grid. Results for several illustrative case studies in two dimensions show the potential of the model to find optimal solutions in few seconds. Injection plans suggested by the optimization model are efficient and have been validated by accurate simulation runs. Based on these findings, the model has the potential to be extended to three dimensions and adapted to real-world cases. [Display omitted]