Adaptive parameterization for solving of thermal/compositional nonlinear flow and transport with buoyancy

The nonlinear nature of flow and transport in porous media requires a linearization of the governing numerical-model equations. We propose a new linearization approach and apply it to complex thermal/compositional problems. The key idea of the approach is the transformation of discretized mass- and energy-conservation equations to an operator...

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...

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...

Influential factors on the development of a low-enthalpy geothermal reservoir: A sensitivity study of a realistic field

A realistic deep low-enthalpy geothermal reservoir based on real data with high detail and complicated sedimentary structure is utilized to perform sensitivity analyses of the geological features influencing reservoir properties. We perform simulations using the Delft Advanced Research Terra Simulator (DARTS). Compelling numerical performance...

Operator-based linearization approach for modeling of multiphase flow with buoyancy and capillarity

Numerical simulation of coupled multiphase multicomponent flow and transport in porous media is a crucial tool for understanding and forecasting of complex industrial applications related to the subsurface. The discretized governing equations are highly nonlinear and usually need to be solved with Newton's method, which corresponds with high...

Discrete well affinity data-driven proxy model for production forecast

A physics-based data-driven model is proposed for forecasting of subsurface energy production. The model fully relies on production data and does not require any in-depth knowledge of reservoir geology or governing physics. In the proposed approach, we use the Delft Advanced Reservoir Terra Simulator (DARTS) as a workhorse for data-driven...

A scalable collocated finite volume scheme for simulation of induced fault slip

We present a scalable collocated Finite Volume Method (FVM) to simulate induced seismicity as a result of pore pressure changes. A discrete system is obtained based on a fully-implicit fully-coupled description of flow, elastic deformation, and contact mechanics at fault surfaces on a flexible unstructured mesh. The cell-centered collocated...

Influence of process-based, stochastic and deterministic methods for representing heterogeneity in fluvial geothermal systems

Focus is on comparing stochastic, process-based and deterministic methods for modelling heterogeneity in hydraulic properties of fluvial geothermal reservoirs. Models are considered a generalized representation of a fluvial sequence in the upper part of the Gassum Formation in northern Denmark. Two ensemble realizations of process-based and...

Uncertainty quantification in a heterogeneous fluvial sandstone reservoir using GPU-based Monte Carlo simulation

The efficient operation and management of a geothermal project can be largely affected by geological, physical, operational and economic uncertainties. Systematic uncertainty quantification (UQ) involving these parameters helps to determine the probability of the focused outputs, e.g., energy production, Net Present Value (NPV), etc. However,...

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...