Applying orthogonal collocation for rapid and reliable solutions of transient flow in naturally fractured reservoirs

The results of OC are verified with exact analytical solutions using the Laplace transform. Thereafter, the effects of three parameters (storativity ratio, interporosity flow coefficient, and radial distance of the outer boundary from the wellbore) are studied using pressure-derivative plots. Employing the superposition theorem, the pressure response of NFR to the multi-rate production from the wellbore is calculated. The results indicate that this approximate approach can predict the pressure transient response of NFRs with an acceptable accuracy, i.e., absolute average relative deviation (AARD) of 0.16%. The equations can be adapted to deal with other flow conditions, boundary conditions and fracture stimulated reservoirs. The proposed method provides a straightforward, faster and more accurate solution in comparison with other methods commonly used. As the OC method is simpler to execute than analytical solutions (e.g. Laplace Transform), it is easier to modify and expand for addressing more complex reservoir scenarios. Likewise, the numerical methods (e.g. finite difference and finite elements) commonly applied to model fractured reservoirs involve significantly more time-consuming calculations, making them more cumbersome to use in reservoir simulations than the proposed orthogonal collocation method.