Modeling of Matrix/Fracture Transfer with Nonuniform-Block Distributions in Low-Permeability Fractured Reservoirs

Unconventional shale-gas and tight oil reservoirs are commonly naturally fractured, and developing these kinds of reservoirs requires stimulation by means of hydraulic fracturing to create conductive fluid-flow paths through open-fracture networks for practical exploitation. The presence of the multiscale-fracture network, including hydraulic fractures, stimulated and nonstimulated natural fractures, and microfractures, increases the complexity of the reservoir simulation. The matrix-block sizes are not uniform and can vary in a very wide range, from several tens of centimeters to meters. In such a reservoir, the matrix provides most of the pore volume for storage but makes only a small contribution to the global flow; the fracture supplies the flow, but with negligible contributions to reservoir porosity. The hydrocarbon is mainly produced from matrix/fracture interaction. So, it is essential to accurately model the matrix/fracture transfers with a reservoir simulator. For the fluid-flow simulation in shale-gas and tight oil reservoirs, dual-porosity models are widely used. In a commonly used dual-porosity-reservoir simulator, fractures are homogenized from a discrete-fracture network, and a shape factor based on the homogenized-matrix-block size is applied to model the matrix/fracture transfer. Even for the embedded discrete-fracture model (EDFM), the matrix/fracture interaction is also commonly modeled using the dual-porosity concept with a constant shape factor (or matrix/fracture transmissibility). However, in real cases, the discrete-fracture networks are very complex and nonuniformly distributed. It is difficult to determine an equivalent shape factor to compute matrix/fracture transfer in a multiple-block system. So, a dual-porosity approach might not be accurate for the simulation of shale-gas and tight oil reservoirs because of the presence of complex multi-scalefracture networks. In this paper, we study the multiple-interacting-continua (MINC) method for flow modeling in fractured reservoirs. MINC is commonly considered as an improvement of the dual-porosity model. However, a standard MINC approach, using transmissibilities derived from the MINC proximity function, cannot always correctly handle the matrix/fracture transfers when the matrix-block sizes are not uniformly distributed. To overcome this insufficiency, some new approaches for the MINC subdivision and the transmissibility computations are presented in this paper. Several examples are presented to show that using the new approaches significantly improves the dual-porosity model and the standard MINC method for nonuniform-block-size distributions.