Decline curve analysis using a pseudo-pressure-based interporosity flow equation for naturally fractured gas reservoirs

Significant amounts of oil and gas are trapped in naturally fractured reservoirs, a phenomenon which has attracted growing attention as production from unconventional reservoirs starts to outpace production from conventional sources. Traditionally, the dual-porosity model has been used in modeling naturally fractured reservoirs. In a dual-porosity model, fluid flows through the fracture system in the reservoir, while matrix blocks are segregated by the fractures and act as fluid sources for them. This model was originally developed for liquid flow in naturally fractured systems and it is therefore inadequate for capturing pressure-dependent effects such as viscosity–compressibility changes in gas systems in its original form. This study presents a rigorous derivation of a gas interporosity flow equation that accounts for the effects of such pressure-sensitive properties. A numerical simulator using the gas interporosity flow equation is built and demonstrates a significant difference in system response from that of a simulator implementing a liquid-form interporosity flow equation. For this reason, rigorous modeling of interporosity flow is considered essential to decline curve analysis for naturally fractured gas reservoirs. In this study, we also show that the use of the proposed gas interporosity flow equation eliminates late-time decline discrepancies and enables rigorous decline curve analysis. The applicability of density-based approach in dual-porosity gas systems is investigated, and the approach reveals that gas production can be forecast in terms of a rescaled liquid solution that uses depletion-driven parameters, λ and β. Application of this approach demonstrated that, at the second decline stage, gas production profile shifted from its liquid counterpart is identical to gas numerical responses with gas interporosity flow equation in effects. The production rates from the pseudo-function approach and those from simulations implementing the gas interporosity flow equation for the synthetic reservoirs are compared against each other, which demonstrated good matches during decline.