Finite Element Model Simulations Associated With Hydraulic Fracturing

Abstract

Abstract Recently emphasis has been placed on the development and testing of innovative well stimulation techniques for the recovery of unconventional gas resources. The design of optimal hydraulic fracturing treatments for specified reservoir conditions requires sophisticated models for predicting the induced fracture geometry and interpreting governing mechanisms. This paper presents methodology and results pertinent to hydraulic fracture modeling for the U.S. DOE's Eastern Gas Shales Program (EGSP). The presented finite-element model simulations extend available modeling efforts and provide a unified framework for evaluation of fracture dimensions and associated responses. Examples illustrating the role of multilayering, in-situ stress, joint interaction, and branched cracks are given. Selected comparisons and applications also are discussed. Introduction Selection and design of stimulation treatments for Devonian shale wells has received considerable attention in recent years1-3. The production of natural gas from such tight eastern petroliferous basins is dependent on the vertical thickness of the organically rich shale matrix, its inherent fracture system density, anisotropy, and extent, and the communication-link characteristics of the induced fracture system(s). The investigation of stimulation techniques based on resource characterization, reservoir property evaluation, theoretical and laboratory model simulations, and field testing is a logical step toward the development of commercial technology for optimizing gas production and related costs. This paper reports formulations, methodology, and results associated with analytical simulations of hydraulic fracturing for EGSP. The presented model extends work reported by Perkins and Kern,4 Nordgren,5 Geertsma and DeKlerk,6 and Geertsma and Haafkens.7 The simulations provide a finite-element model framework for studying vertically induced fracture responses with the effects of multilayering and in-situ stress considered. In this context, Brechtel et al.,8 Daneshy,9 Cleary,10 and Anderson et al.11 have done recent studies addressing specific aspects of this problem. The use of finite-element model techniques for studying mixed-mode fracture problems encountered in dendritic fracturing and vertical fracture/joint interaction also is illustrated along with application of suitable failure criteria. Vertical Hydraulic Fracture Model Formulations Coupled structural fracture mechanics and fracture fluid response models for predicting hydraulically induced fracture responses have been reported previously.12,13 These simulations incorporate specified reservoir properties, in-situ stress conditions, and stimulation treatment parameters. One shortcoming of this modeling effort is that finite-element techniques are used for the structural and stress intensity simulations, while a finite-difference approach is used to evaluate the leakoff and fracture-fluid response in the vertical crack. A consistent framework for conducting all simulations using finite-element modeling is formulated here.