An upscaling procedure for fractured reservoirs with embedded grids

Upscaling of geological models for reservoir simulation is an active and important area of research. In particular, we are interested in reservoirs where the rock matrix exhibits an intricate network of fractures, which usually acts as a preferential path to the flow. Accounting for fractures' contribution in the simulation of a reservoir is of paramount importance. Here we have focused on obtaining effective parameters (e.g., transmissibility) on a 3-D computational grid on the reservoir scale, which account for the presence, at a finer spatial scale, of fractures and a network of fractures. We have essentially followed the idea illustrated in Karimi-Fard et al. (2006), yet this work has some notable aspects of innovation in the way the procedure has been implemented, and in its capability to consider rather general corner-point grids, like the ones normally used in reservoir simulations in the industry, and complex and realistic fracture networks, possibly not fully connected inside the coarse cells. In particular, novel contribution is the employment of an Embedded Discrete Fracture Model (EDFM) for computing fracture-fracture and matrix-fracture transmissibilities, with a remarkable gain in speedup. The output is in the form of transmissibility that, although obtained by considering single-phase flow, can be used for coarse-scale multiphase reservoir simulations, also via industrial software, such as Eclipse, Intersect, or GPRS. The results demonstrate the effectiveness and computational efficiency of the numerical procedure which is now ready for further testing and industrialization