Predicting productivity index of hydraulically fractured formations

Abstract Productivity index and inflow performance of horizontal wells intersecting multiple hydraulic fractures are of great importance. This importance comes from the fact that the fracturing process has become a common stimulation technique in the petroleum industry. However, few models for the productivity index and inflow performance have been presented in the literatures due to the complexity governing this topic. This paper introduces a new technique for estimating the pseudo-steady state productivity index of horizontal wells intersecting multiple hydraulic fractures. Based on the instantaneous source solutions for the diffusivity equation, seven analytical models have been derived for different source solutions. Four of them represent the effect of the formation height and fracture height (the vertical direction), while the other three represent the solution for the horizontal plane. For vertical hydraulic fractures, the four solutions of the vertical direction, representing the pseudo-skin factor, are almost neglected. The three horizontal plane solutions are the main parameters that control the productivity index and inflow performance of the fractured formations. In this technique, the horizontal wells are acting in finite reservoirs where the pseudo-steady state flow is expected to develop. Reservoir geometry, reservoir properties, and fracture dimensions were considered in this technique. The number of fractures and the spacing between them were also investigated in this study. A new analytical model for estimating the required number of hydraulic fractures has been introduced in this study based on the reservoir drainage area and the surface area of fractures. The models have been used to establish several plots to estimate the shape factor group based on the number of fractures and the half fracture length. This group is one of the main terms in the productivity index model. Several plots for the shape factor of fractured formations have been introduced in this study. The results obtained from the new technique have been compared with the results from previous models. Several numerical examples will be included in the paper.