Due to numerical difficulties in conducting high fidelity simulation of recovery mechanisms in complex natural fracture systems, there are no published studies that address the impact of preserving details of the fracture networks. We used highly refined grids to conduct fine scale simulations of various recovery mechanisms in different complex fracture settings and compared the results to those obtained on simplified dual porosity dual permeability (DPDK) representations created by applying a consistent upscaling procedure.

Our study considers densely connected, sparsely connected, and isolated fracture networks that are extracted from a field-scale fractured carbonate reservoir model. Discrete fracture-matrix (DFM) models were constructed using an unstructured grid with refinement of the matrix rock near fractures. High-resolution simulations of spontaneous imbibition, gravity drainage, and viscous displacement recovery mechanisms were conducted on these DFM models. We also built equivalent DPDK models by using single phase flow-based upscaling and actual fracture geometry and distribution. The recovery mechanisms were simulated on these DPDK models and compared to high-resolution DFM models.

The fine scale simulations revealed that lateral viscous displacement recovery depends on the details of the fracture networks and can be significantly higher than those predicted from equivalent DPDK models. The DPDK models all predict the same recovery. For spontaneous imbibition, both fine scale and equivalent DPDK models show dependence on fracture geometry, but the DPDK models predict much higher rates. Fine scale and equivalent DPDK models agree reasonably for gravity drainage. These findings are explained by analyzing the matrix-fracture flows, and implications on efforts to improve shape factors in DPDK models and upscaling efforts in DFM models are discussed.

You can access this article if you purchase or spend a download.