Abstract
This paper presents an innovative numerical procedure for flow simulation and estimation of relative permeability of naturally fractured carbonate rocks by history matching laboratory derived multi-phase production data. The multi-phase flow simulator accounts for flow in discrete fractures with variable asperities. These discrete fractures are characterised by the fracture surface topography which is imaged with the use of surface laser scanning and micro computed tomography (mu-CT). Next, the production data for naturally fractured carbonate samples (matrix and fracture permeability), samples with no fracture (matrix permeability only) and samples with fracture only (near zero matrix permeability) for unsteady state conditions are obtained. Forward simulation performed by coupling both the matrix and fracture relative permeability to estimate production for a porous fractured system. This is validated by comparing to the laboratory obtained production data. This has allowed us to validate the proposed numerical procedure and subsequently history match the experimental production data for the porous fractured system to obtain overall porous fractured relative permeability curves. Results of this study show that there exists a strong discontinuity of phase saturation. It was also observed that Corey's parameters are best suited for history matching of the production data through porous matrix, and cubic parameters are best suited for fracture.