Two-dimensional models for waterflooding induced hydraulic fracture accounting for the poroelastic effects on a reservoir scale

Abstract

This work focuses on the developing of a computationally efficient numerical model for a waterflooding induced planar hydraulic fracture in a poroelastic medium. In particular, under certain assumptions it is possible to decouple the fully coupled model into poroelastic reservoir and purely elastic fracture subproblems. Further, the general 3D problem formulation for a planar fracture is reduced to two dimensions for early time plane strain fracture and late time constant height fracture. Noticeably, the division of the problem into two subproblems allows us to use the displacement discontinuity approach for the fracture related part, and finite element method for the reservoir related part. The latter is especially important for the constant height fracture, which cannot be modeled using exclusively the finite element approach in two dimensions. The algorithm is benchmarked against the available analytical solutions, as well as compared to other numerical codes. The developed numerical method is used to study the influence of the heterogeneous pore pressure and the corresponding stress for different waterflooding fracturing situations. We consider the cases with both injection and production wells operating nearby. With the presence of such wells, the fracture growth becomes asymmetrical and both the well type and location are important. The qualitatively different fracture behavior is observed for various cases. Other numerical examples include the simulation of a step-rate test as well as investigation of the possibility of fracture control by altering the injection rate.