scholarly journals Semi-Implicit Finite Volume Procedure for Compositional Subsurface Flow Simulation in Highly Anisotropic Porous Media

Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 341
Author(s):  
Sebastián Echavarría-Montaña ◽  
Steven Velásquez ◽  
Nicolás Bueno ◽  
Juan David Valencia ◽  
Hillmert Alexander Solano ◽  
...  

Subsurface multiphase flow in porous media simulation is extensively used in many disciplines. Large meshes with non-orthogonalities (e.g., corner point geometries) and full tensor highly anisotropy ratios are usually required for subsurface flow applications. Nonetheless, simulations using two-point flux approximations (TPFA) fail to accurately calculate fluxes in these types of meshes. Several simulators account for non-orthogonal meshes, but their discretization method is usually non-conservative. In this work, we propose a semi-implicit procedure for general compositional flow simulation in highly anisotropic porous media as an extension of TPFA. This procedure accounts for non-orthogonalities by adding corrections to residual in the Newton-Raphson method. Our semi-implicit formulation poses the guideline for FlowTraM (Flow and Transport Modeller ) implementation for research and industry subsurface purposes. We validated FlowTraM with a non-orthogonal variation of the Third SPE Comparative Solution Project case. Our model is used to successfully simulating a real Colombian oil field.

2006 ◽  
Author(s):  
Hong Wang ◽  
Richard Ewing ◽  
Guan Qin ◽  
Stephen Lincoln Lyons

2019 ◽  
Vol 46 (12) ◽  
pp. 1151-1159 ◽  
Author(s):  
Naser Moosavian

In this paper, a new matrix framework has been developed for the simulation of flow and pressure in porous media. In this framework, the pressure gradient formulation in Darcy’s law is considered as the head-loss equation in pipe network modeling. Then, an artificial pipe network has been constructed to find the pressure head profile in porous media. Two explicit and implicit formulations have been advanced for linear and nonlinear analysis, which the latter is an implementation of the Newton–Raphson algorithm. Both formulations iteratively solve a linear system of equations for calculating the nodal heads and apply a matrix multiplication for updating the flow vector. While the explicit method needs few iterations, the implicit method requires at least 20 iterations to converge with acceptable accuracy. For testing these formulations, four different types of network configurations were tested. The analysis of three laboratory tests showed that the application of the implicit method provides reliable and accurate results.


Water ◽  
2019 ◽  
Vol 11 (12) ◽  
pp. 2647 ◽  
Author(s):  
Dominique Guérillot ◽  
Jérémie Bruyelle

One critical point of modeling of flow in porous media is the capacity to consider parameters that are highly variable in space. It is then very challenging to simulate numerically fluid flow on such heterogeneous porous media. The continuous increase in computing power makes it possible to integrate smaller and smaller heterogeneities into geological models of up to tens of millions of cells. On such meshes, despite computer performance, multi-phase flow equations cannot be solved in an acceptable time for hydrogeologists and reservoir engineers, especially when the modeling considers several components in each fluid and when taking into account rock-fluid interactions. Taking average reservoir properties is a common approach to reducing mesh size. During the last decades, many authors studied the upscaling topic. Two different ways have been investigated to upscale the absolute permeability: (1) an average of the permeability for each cell, which is then used for standard transmissibility calculation, or (2) computing directly the upscaled transmissibility values using the high-resolution permeability values. This paper is related to the second approach. The proposed method uses the half-block approach and combines the finite volume principles with algebraic methods to provide an upper and a lower bound of the upscaled transmissibility values. An application on an extracted map of the SPE10 model shows that this approach is more accurate and faster than the classical transmissibility upscaling method based on flow simulation. This approach keeps the contrast of transmissibility values observed at the high-resolution geological scale and improves the accuracy of field-scale flow simulation for highly heterogeneous reservoirs. Moreover, the upper and lower bounds delivered by the algebraic method allow checking the quality of the upscaling and the gridding.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-12
Author(s):  
Zhuangzhuang Wang ◽  
Zhaomin Li ◽  
Hailong Chen ◽  
Fei Wang ◽  
Dawei Hou ◽  
...  

Foam is widely used as a selective blocking agent through mobility control in oil field development. Its flow behavior in porous media has been investigated sufficiently, but few studies were carried out to understand the change of foam texture in flow. In this work, sandpack and micromodel experiments were conducted simultaneously to analyze foam flow behavior from the perspective of foam texture. Based on the measured flowing pressure and the observed foam image, the correlation between blocking pressure and foam texture was quantitatively investigated. The blocking pressure has a strong correlation with average diameter (-0.906) and variation coefficient (-0.78) and has a positive correlation with the filling ratio (0.84). These indicate that the blocking performance of foam is influenced by its texture closely. But path analysis shows only that the average diameter and variation coefficient have a significant direct effect on blocking pressure (-0.624 and -0.404). These show that the blocking capacity of foam is mainly influenced by the size and uniformity of bubbles. Tiny, dense, and homogeneous foam has a stronger blocking capacity. This study provides a deep insight of foam flow in porous media.


2016 ◽  
Vol 19 (4) ◽  
pp. 998-1014 ◽  
Author(s):  
Lei Zhang ◽  
Qinjun Kang ◽  
Li Chen ◽  
Jun Yao

AbstractThe unified lattice Boltzmann model is extended to the quadtree grids for simulation of fluid flow through porous media. The unified lattice Boltzmann model is capable of simulating flow in porous media at various scales or in systems where multiple length scales coexist. The quadtree grid is able to provide a high-resolution approximation to complex geometries, with great flexibility to control local grid density. The combination of the unified lattice Boltzmann model and the quadtree grids results in an efficient numerical model for calculating permeability of multi-scale porous media. The model is used for permeability calculation for three systems, including a fractured system used in a previous study, a Voronoi tessellation system, and a computationally-generated pore structure of fractured shale. The results are compared with those obtained using the conventional lattice Boltzmann model or the unified lattice Boltzmann model on rectangular or uniform square grid. It is shown that the proposed model is an accurate and efficient tool for flow simulation in multi-scale porous media. In addition, for the fractured shale, the contribution of flow in matrix and fractures to the overall permeability of the fractured shale is studied systematically.


Sign in / Sign up

Export Citation Format

Share Document