Hybrid Coupled Discrete-Fracture/Matrix and Multicontinuum Models for Unconventional-Reservoir Simulation

Abstract

Summary Unconventional reservoirs are the focus of considerable attention as a primary energy source. Numerical simulation is a core kernel of reservoir-engineering work flows for reservoir evaluation, optimization, and management. Accurate and efficient numerical simulation of unconventional reservoirs is challenging. There is substantial physical complexity involving a number of tightly coupled mechanisms in the modeling of these reservoirs. The complexity is further amplified by the multicontinuum nature of the stimulated formation, and the complex fracture networks with a wide range of fracture-length scales and topologies. To adequately capture the effects of the multiscaled fracture system, we develop two alternative hybrid approaches that are aimed at combining the advantages of multicontinuum and discrete-fracture/matrix (DFM) representations. During the development of unconventional resources, geological and geophysical information may be available in some cases to suggest a prior characterization, whereas in many other cases, this prior model may be incomplete and limited to hydraulic fractures. The two hybrid approaches could be used for different applications depending on the available characterization data and different requirements for efficiency and accuracy considerations. The first hybrid model couples an embedded-discrete-fracture model (EDFM) with multiple interacting continua (MINC) into EDFM/MINC, which simulates the fracture network characterized by stimulated-reservoir-volume (SRV) concept. This optimized model can reduce the computational cost that is associated with the widely applied logarithmically spaced/locally refined (LS/LR) DFM technique, while improving the flexibility to model the complex geometry of hydraulic fractures. The MINC concept allows the hybrid model to handle the extreme contrast in conductivity between the small-scale fracture network and the ultratight matrix that results in steep potential gradients. For the second class of hybrid model (unstructured DFM/continuum), the primary fractures are described by use of DFM with unstructured gridding, and the small-scale fractures are simulated by continuum-type approaches in a fully coupled manner. Optimized local-grid refinement is used to accurately handle the transient-flow regime around primary fractures. An upscaling technique that applies EDFM on the detailed realization of the discrete-fracture network by use of the target unstructured grid to generate an appropriate dual-permeability model is also developed. The upscaling technique is suitable for cases where a detailed prior model for the complete fracture network is available. Simulation studies demonstrate the applicability of the developed hybrid-fracture models. Model verification is conducted against several reference solutions.

A Comprehensive Model for Real Gas Transport in Shale Formations with Complex Non-planar Fracture Networks

“…(2012) integrated a shale-gas transport model with extended finite element method (XFEM) to study the main flow gas mechanism of shale in complex fracture network18. Jiang and Younis (2015) proposed two hybrid approaches: one is the coupling EDFM with multiple interacting continua (MINC), the other is the coupling of unstructured DFM with continuum-type approaches19. However, these numerical methods are still challenging to apply due to complicated gridding issues, an expensive computational cost, and complexities in development of computational codes.…”

Flow and Transport in Tight and Shale Formations: A Review

“…The embedded discrete fracture models have been employed by Li and Lee [128] and Moinfar et al [129] discretized the complex fractures into a number of segments while the matrix is treated as structured grids. Jiang and Younis [130] used the two hybrid methods, namely, embedded discrete fracture models with multiple interacting continua and the coupling of unstructured DFM with continuum type. The fracture networks' complexity has important effects on the well performance.…”

Simulation of a multistage fractured horizontal well in a tight oil reservoir using an embedded discrete fracture model

“…Based on the above advantages, 23,29-40 many scholars have performed numerous simulations using the EDFM, demonstrating its good adaptability for simulating fractured reservoirs. 23,[29][30][31][32][33][34][35][36][37][38][39][40] The EDFM can describe the interporosity flow between different media, so each cell is connected with different numbers and types of cells. To understand and record these connections 30 proposed the non-neighboring connection (NNC) method, and Jiang et al 31 proposed the connection list, but they did not study their methods further.…”

“…The EDFM uses an orthogonal structured grid for the mesh matrix, and the fractures are divided into a series of segments based on the boundary of the matrix cell. Based on the above advantages, 23,[29][30][31][32][33][34][35][36][37][38][39][40] many scholars have performed numerous simulations using the EDFM, demonstrating its good adaptability for simulating fractured reservoirs. The EDFM is very flexible for handing various kinds of intersections, such as fracture-fracture and fracture-well intersections 28 .…”

“…These segments are the fracture grids in the EDFM. 23,[29][30][31][32][33][34][35][36][37][38][39][40] The EDFM can describe the interporosity flow between different media, so each cell is connected with different numbers and types of cells. The interporosity between two systems is calculated using the Peaceman-like formulation, which appears as the source/sink term in the governing equations.…”