Premium
An upscaling procedure for fractured reservoirs with embedded grids
Author(s) -
Fumagalli Alessio,
Pasquale Luca,
Zonca Stefano,
Micheletti Stefano
Publication year - 2016
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2015wr017729
Subject(s) - fracture (geology) , grid , computer science , matrix (chemical analysis) , scale (ratio) , flow (mathematics) , regular grid , speedup , transmissibility (structural dynamics) , petroleum engineering , geology , simulation , geotechnical engineering , mathematics , geometry , materials science , parallel computing , geodesy , physics , quantum mechanics , vibration isolation , vibration , composite material
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.