In this work, we consider a single-phase flow and heat transfer problem in fractured geothermal reservoirs. Mixed dimensional problems are considered, where the temperature and pressure equations are solved for porous matrix and fracture networks with transfer term between them. For the fine-grid approximation, a finite volume method with embedded fracture model is employed. To...

This paper presents a non-intrusive subdomain POD-TPWL (SD POD-TPWL) for reservoir history matching through integrating domain decomposition (DD), proper orthogonal decomposition (POD), radial basis function (RBF) interpolation, and the trajectory piecewise linearization (TPWL). It is an efficient approach for model reduction and linearization of general non-linear time-dependent...

Field-scale simulation of flow in porous media in presence of incomplete mixing demands for high-resolution computational grids, much beyond the scope of state-of-the-art simulators. Hence, the upscaling-based Todd and Longstaff (TL) approach is typically used, where coarse grid cells are employed with effective mixing fluid properties and parameters found by matching results...

In this work, we model the biofilm growth at the microscale using a rectangular pore network model in 2D and a cubic network in 3D. For the 2D network, we study the effects of bioclogging on porosity and permeability when we change parameters like the number of nodes in the network, the network size, and the concentration of nutrients at the inlet. We use a 3D cubic network to...

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to...

Reactive transport in porous media with dissolution and precipitation has important applications in oil and gas industry and groundwater remediation. In this work, we present a simulation method for reactive flow in porous media of two salts that share an ion. The method consists of a front-tracking solver that uses the Riemann solutions of the underlying set of hyperbolic...

Three sets of synthetic images were created from two original datasets. A suite exhibiting greyscale contrast was produced from an 8.96-μm voxel size 3D X-ray microscopy image of a sandstone rock and a two suites (one showing greyscale contrast and one showing both greyscale and textural contrast) were produced from a 5 × 5 × 5 nm voxel size FIB-SEM image of a shale rock. The...

Flows of multiple fluid phases are common in many subsurface reservoirs. Numerical simulation of these flows can be challenging and computationally expensive. Dynamic adaptive mesh optimisation and related approaches, such as adaptive grid refinement can increase solution accuracy at reduced computational cost. However, in models or parts of the model domain, where the local...

There are several issues to consider when we use ensemble smoothers to condition reservoir models on rate data. The values in a time series of rate data contain redundant information that may lead to poorly conditioned inversions and thereby influence the stability of the numerical computation of the update. A time series of rate data typically has correlated measurement errors...

Efficient heat exploitation strategies from geothermal systems demand for accurate and efficient simulation of coupled flow-heat equations on large-scale heterogeneous fractured formations. While the accuracy depends on honouring high-resolution discrete fractures and rock heterogeneities, specially avoiding excessive upscaled quantities, the efficiency can be maintained if...

This paper examines the properties of the Iterated Ensemble Smoother (IES) and the Multiple Data Assimilation Ensemble Smoother (ES–MDA) for solving the history matching problem. The iterative methods are compared with the standard Ensemble Smoother (ES) to improve the understanding of the similarities and differences between them. We derive the three smoothers from Bayes...

Hydrological models are necessary tools for simulating the water cycle and for understanding changes in water resources. To achieve realistic model simulation results, real-world observations are used to determine model parameters within a “calibration” procedure. Optimization techniques are usually applied in the model calibration step, which assures a maximum similarity between...

We present an analytical solution to estimate the minimum polymer slug size needed to ensure that viscous fingering of chase water does not cause its breakdown during secondary oil recovery. Polymer flooding is typically used to improve oil recovery from more viscous oil reservoirs. The polymer is injected as a slug followed by chase water to reduce costs; however, the water is...

We investigate the accuracy and robustness of one of the most common methods used in glaciology for finite element discretization of the ?-Stokes equations: linear equal order finite elements with Galerkin least-squares (GLS) stabilization on anisotropic meshes. Furthermore, we compare the results to other stabilized methods. We find that the vertical velocity component is more...

In this paper, we apply the enhanced local pressure (ELP) model to study crack interaction in hydraulic fracturing. The method is based on the extended finite element method (X-FEM) where the pressure and the displacement fields are assumed to be discontinuous over the fracture exploiting the partition of unity property of finite element shape functions. The material is fully...

In simulation of fluid injection in fractured geothermal reservoirs, the characteristics of the physical processes are severely affected by the local occurence of connected fractures. To resolve these structurally dominated processes, there is a need to develop discretization strategies that also limit computational effort. In this paper, we present an upscaling methodology for...

Finite difference is the most widely used method for seismic wavefield modeling. However, most finite-difference implementations discretize the Earth model over a fixed grid interval. This can lead to irregular model geometries being represented by ‘staircase’ discretization, and potentially causes mispositioning of interfaces within the media. This misrepresentation is a major...

Due to an oversight, some author’s corrections were not carried out during Performing proof corrections stage. The Publisher apologizes for these mistakes. The original article was updated.

Grid generation for reservoir simulation must honor classical key constraints and be boundary aligned such that control-volume boundaries are aligned with geological features such as layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. An unstructured grid generation procedure is proposed that automates control-volume and/or control point boundary...

Novel cell-centred finite-volume formulations are presented for incompressible and immiscible two-phase flow with both gravity and capillary pressure effects on structured and unstructured grids. The Darcy-flux is approximated by a control-volume distributed multipoint flux approximation (CVD-MPFA) coupled with a higher resolution approximation for convective transport. The CVD...