Algebraic multiscale method for flow in heterogeneous porous media with embedded discrete fractures (F-AMS)

Matei Tene,Mohammed Saad Al Kobaisi,Hadi Hajibeygi

DOI: https://doi.org/10.1016/j.jcp.2016.06.012

2017-05-15

Abstract:This paper introduces an Algebraic MultiScale method for simulation of flow in heterogeneous porous media with embedded discrete Fractures (F-AMS). First, multiscale coarse grids are independently constructed for both porous matrix and fracture networks. Then, a map between coarse- and fine-scale is obtained by algebraically computing basis functions with local support. In order to extend the localization assumption to the fractured media, four types of basis functions are investigated: (1) Decoupled-AMS, in which the two media are completely decoupled, (2) Frac-AMS and (3) Rock-AMS, which take into account only one-way transmissibilities, and (4) Coupled-AMS, in which the matrix and fracture interpolators are fully coupled. In order to ensure scalability, the F-AMS framework permits full flexibility in terms of the resolution of the fracture coarse grids. Numerical results are presented for two- and three-dimensional heterogeneous test cases. During these experiments, the performance of F-AMS, paired with ILU(0) as second-stage smoother in a convergent iterative procedure, is studied by monitoring CPU times and convergence rates. Finally, in order to investigate the scalability of the method, an extensive benchmark study is conducted, where a commercial algebraic multigrid solver is used as reference. The results show that, given an appropriate coarsening strategy, F-AMS is insensitive to severe fracture and matrix conductivity contrasts, as well as the length of the fracture networks. Its unique feature is that a fine-scale mass conservative flux field can be reconstructed after any iteration, providing efficient approximate solutions in time-dependent simulations.

Numerical Analysis,Computational Engineering, Finance, and Science,Computational Physics,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the flow simulation problem of embedding discrete fractures in heterogeneous porous media. Specifically, the paper proposes an Algebraic Multiscale Method (F - AMS) for flow simulation in heterogeneous porous media with embedded discrete fractures. This method aims to overcome the computational complexity and efficiency problems encountered by traditional simulation methods when dealing with high - resolution fracture networks.

### Background and Challenges
In many geoscience applications, such as hydrocarbon production and geothermal energy development, the target formations are often naturally fractured. These formations are usually composed of highly heterogeneous bedrocks, which are interspersed with multiple low - dimensional and high - conductivity fracture networks at different scales. This poses important challenges to flow simulation and prompts researchers to develop advanced modeling and numerical solution techniques.

### Limitations of Existing Methods
The existing Multiscale Finite Element (MSFE) and Multiscale Finite Volume (MSFV) methods are mainly applicable to non - fractured heterogeneous porous media. These methods have limitations when dealing with fracture networks, especially when the conductivity differences between fractures and bedrocks are significant and the lengths of fracture networks are long, the convergence will be affected.

### Innovations of the F - AMS Method
The main innovations of the F - AMS method include:
1. **Independently Construct Multiscale Coarse Grids**: Independently construct multiscale coarse grids for the bedrock and the fracture network.
2. **Locally - Supported Basis Functions**: Calculate basis functions with local support through algebraic methods and establish the mapping relationship between coarse and fine scales.
3. **Multiple Coupling Strategies**: Four types of basis functions are studied, namely fully decoupled (Decoupled - AMS), unidirectional coupling (Frac - AMS and Rock - AMS), and fully coupled (Coupled - AMS).
4. **Flexibility**: The F - AMS framework allows complete flexibility in the resolution of the fracture coarse grid, thus adapting to fracture networks at different scales.
5. **Efficiency and Scalability**: Through systematic benchmark tests, the F - AMS method shows good performance and scalability when dealing with high - conductivity - ratio and long - fracture - network situations.

### Numerical Results and Verification
The paper verifies the performance of the F - AMS method through two - dimensional and three - dimensional numerical experiments. The experimental results show that the F - AMS method has high efficiency and accuracy when dealing with high - conductivity - ratio and long - fracture - network situations, and can provide effective approximate solutions in time - dependent simulations.

### Conclusion
The F - AMS method provides an efficient and scalable solution to the problem of flow simulation of embedding discrete fractures in heterogeneous porous media. The development of this method marks an important step in integrating multiscale methods as "black - box" pressure solvers into existing reservoir simulators, and is expected to be extended to more complex physical scenarios in the future.