A Unified Thermal-Reactive Compositional Simulation Framework for Modeling CO₂ Sequestration at Various Scales

In this work, we present a unified framework for the simulation of CO₂ sequestration problems at various time and space scales. The parametrization technique utilizes thermodynamic state-dependent operators expressing the governing equations for the thermal-compositional-reactive system to solve the nonlinear problem. This approach...

Coupled modeling of well and reservoir for geo-energy applications

The energy transition is inevitable since approximately two-thirds of the current global GHG emissions are related to energy production. Subsurface can provide a great opportunity for innovative low-carbon energy solutions such as geothermal energy production, hydrogen storage, carbon capture, and sequestration, etc. Well and borehole...

Applying Massively Parallel Interface for MPFA scheme with advanced linearization for fluid flow in porous media

We apply Massively Parallel Interface for MPFA-O scheme with state-of-the-art Operator-Based Linearization (OBL) approach for multiphase flow in porous media. The implementation of MPFA-O scheme enhances the modelling capabilities for non-K-orthogonal mesh. A fully implicit scheme is applied to guarantee the stability of solutions when a mass...

Nonlinear solver based on trust region approximation for CO₂ utilization and storage in subsurface reservoir

Simulation of CO₂ utilization and storage (CCUS) in subsurface reservoirs with complex heterogeneous structures requires a model that captures multiphase compositional flow and transport. Accurate simulation of these processes necessitates the use of stable numerical methods that are based on an implicit treatment of the flux term...

FluidFlower Benchmark: Lessons Learned from the Perspective of Subsurface Simulation

In this work, we describe our decisions made to perform the FluidFlower simulation study and discuss various aspects of the benchmark that are different from our usual subsurface simulation practice. We will discuss the impact of various modeling choices on the outcomes of the simulation models, such as gridding, discretization, and solver...

Towards the Uncertainty Quantification of Fractured Karst Systems: Reactive Transport and Fracture Networks: Where Numerical Modeling Meets Outcrop Observations

Society relies on large amounts of energy to progress and allow for a high standard of living. The recent severe climate changes require advanced technologies related to cleaner energy resources. One such technology beneficial for accelerating this current energy transition is geothermal energy. This type of energy is often found in fractured...

Operator-Based Modeling and Inversion: An Operator Approach to the Forward and Inverse Scattering Problems

The seismic method has many applications. It is important in the critical sector of energy. Besides being used in imaging oil and gas reservoirs, it is also utilized in other sectors of energy such as geothermal energy exploration and development. It also plays a role in extracting other resources such as minerals or in the process of monitoring...

Numerical investigations of foam-assisted CO₂ storage in saline aquifers

CO₂-foam injection is a promising technology for reducing gas mobility and increasing trapping within the swept region in deep brine aquifers. In this work, a consistent thermodynamic model based on a combination of the Peng-Robinson equation of state (PR EOS) for gas components with an activity model for the aqueous phase is...

An efficient numerical simulator for geothermal simulation: A benchmark study

Accurate prediction of temperature and pressure distribution is essential for geothermal reservoir exploitation with cold water re-injection. Depending on our knowledge about the heterogeneous structure of the subsurface, the reservoir development scheme can be optimized and the overall lifetime of the geothermal field can be extended. In...

Delft Advanced Research Terra Simulator: General Purpose Reservoir Simulator with Operator-Based Linearization

Numerical simulation is based on space and time discretization providing a trade-off between accuracy and computational performance. The operator-based Linearization (OBL) approach introduces an additional discretization domain for the physical description of fluid and rock. Delft Advanced Research Terra Simulator (DARTS) is a general-purpose...

Element balance formulation in reactive compositional flow and transport with parameterization technique

We present a novel nonlinear formulation for modeling reactive-compositional flow and transport in the presence of complex phase behavior due to a combination of thermodynamic and chemical equilibria in multi-phase systems. We apply this formulation to model precipitation/dissolution of minerals in reactive multiphase flow in subsurface...

Optimization of CO₂ injection using multi-scale reconstruction of composition transport

The current situation with green gas emission requires the development of low-carbon energy solutions. However, a significant part of the modern energy industry still relies on fossil fuels. To combine these two contradictory targets, we investigate a strategy based on a combination of CO₂ sequestration with enhanced oil recovery ...

Operator-based linearization for efficient modeling of geothermal processes

Numerical simulation is one of the most important tools required for financial and operational management of geothermal reservoirs. The modern geothermal industry is challenged to run large ensembles of numerical models for uncertainty analysis, causing simulation performance to become a critical issue. Geothermal reservoir modeling requires...

Operator-based linearization approach for modeling of multiphase multi-component flow in porous media

A new approach for the linearization of governing equations that describe flow and transport in porous media is proposed in this work. It is based on an approximate representation of the exact physics of the problem, which is similar to an approximate representation of space and time discretization performed in conventional simulation. The...