Modeling Infinite-Conductivity Vertical Fractures With Source and Sink Terms
Nghiem, Long X., SPE, Computer Modeling Group Abstract This paper describes a method for handling infinite-conductivity vertical fractures in reservoir simulation by using source and sink terms. It begins with a review of the concept of source/sink in reservoir simulation, and uses that concept to develop a method for computing the flow into or out of the fracture by assuming elliptical tow and by using the pressures of the blocks surrounding those containing the fracture. The assumption that the flow into the fracture is everywhere perpendicular to the fracture plane (i. e., linear flow) and the effect of the skin factor are also investigated. Test runs showed excellent agreement between computed results and those obtained by analytical and variational methods for single-phase systems. The formulation was extended to multiphase systems. and simulation of a waterflood yielded physically reasonable results. Introduction The simulation of fractured wells is of considerable interest because of the large number of wells that have been hydraulically fractured to increase the injectivity/ productivity in low-permeability formations. Hydraulic productivity in low-permeability formations. Hydraulic fracturing usually yields a vertical fracture plane that intersects the wellbore. Agarwal et al. showed that hydraulic fractures for which the dimensionless fracture flow capacity, FCD, defined as . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (1) is greater than 500 can be represented by infinite-conductivity vertical fractures. In Eq. 1 kf is the fracture permeability, vi, is the fracture width, k is the formation permeability, and x., is the fracture half-length. permeability, and x., is the fracture half-length. Conventional hydraulic fractures as opposed to massive hydraulic fractures usually fall into this category. Conceptually, an infinite-conductivity vertical fracture (ICVF) has no thickness and no pressure drop along the fracture plane and is represented by a finite line source or sink in an areal representation of the reservoir. Russell and Truitt simulated single-phase flow into an ICVF parallel to a boundary of a square reservoir by using finite differences. The fracture was located symmetrically within a no-flow boundary drainage area and was treated by these authors as a boundary condition. A more general problem was later solved by Gringarten et al., who used the analytical method of source/sink and Green's function proposed by Gringarten and Ramey. Gringarten et al. considered a fracture within a no-flow boundary rectangular reservoir. Using the same analytical method, Raghavan and Hadinoto obtained solutions for a fracture within a constant-pressure outer boundary. In both cases, the fracture was parallel to a boundary of the reservoir. Using Galerkin's method, Bennett et al. investigated the case of a fracture lying along one of the diagonals of a constant-pressure outer boundary square reservoir. The solutions for all these cases are reported in the form of tables and plots of dimensionless pressure drop vs. dimensionless time for different fracture-penetration ratios. This paper presents a method for modeling ICVF's with source/sink terms. The fracture is treated as a singularity and it is assumed that elliptical flow applies in the neighborhood of the fracture. The flow into or out of the fracture is computed from the fracture pressure and the pressures of the blocks surrounding those containing the fracture. SPEJ p. 633