Numerical Simulation of Hydraulically Induced Fracture Network Propagation in Shale Formation

Abstract

Hydraulic fracturing treatment has been proven to be one of the key technologies for shale gas development. Micro-seismic mapping data has shown that hydraulic fracturing stimulation has often resulted in complex fracture network due to the geological complexity of shale formations. Hydraulic fracturing is a coupled process of shale formation deformation and flow of engineered fluid that includes water, proppant and other chemicals. Moreover, the pre-existing natural fractures in shale formation may complicate hydraulic fracture propagation process and alter its Young's modulus. In this paper, we have developed a numerical model for modeling hydraulic fracture propagation in highly fractured shale formation. The proposed numerical model has integrated turbulent flow, rock stress response, interactions of hydraulic fracture propagation with natural fractures, and influence of natural fractures on formation Young's modulus. Mixed finite element method is employed for numerical solution of the nonlinear partial differential equations. The proposed model has been validated with bi-wing hydraulic fracture model through regression tests. The preliminary numerical results illustrate the significant differences in modeling hydraulic fracture propagation in comparison with current models that assume laminar flow in hydraulic fracture process. Two simulated cases with different initial natural fracture mapping are given. The preliminary numerical results show that length and density of natural fracture have significant impact on formation Young's modulus, and interactions between hydraulic fracture and natural fractures create the complex fracture network. This model provides an opportunity to optimize hydraulic fracturing stimulation design through numerical simulations, which is vital in shale reservoir development. Introduction Hydraulic fracturing stimulation is one of the key technologies in shale gas development and has often, as shown in micro seismic mapping, resulted in complex fracture network geometry due to the interaction between hydraulic fracture and pre-existing natural fractures (Fisher et al. 2002; Maxwell et al. 2002; Daniels et al. 2007; Le Calvez et al. 2007; X. Wang et al. 2011). Though it has been widely used in oil and gas industry for many years, it remains a great challenge to quantitatively characterize the hydraulically induced fracture network in shale gas development due to the complexity of the network and lack of the observational data. Hydraulic fracturing is a stimulation process, in which engineered fluids, a mixture of water, sand and chemicals, are pumped into shale formation at high pressure and rate to create an induced fracture network. In the stimulation process, high flow rate and fluid pressure initiate and propagate the hydraulically induced fracture. The induced fracture can open the natural fracture or cross it, or can be rested by the natural fracture once it intersects with a natural fracture (Blanton, 1982; Warpinski, N.R, 1987; Renshaw, C.E., 1995, and N. Potluri, 2005). Consequently, it is critical to model the coupling effects of the flow process and the fracture propagation.

Numerical Simulation of Shale Gas Development Based on Complex Fracture Network Growth

“…The basic data is shown in Table 2, and multi-stage hydraulic fracture distribution is shown is Figure 3. Using the above well data well producing rate is calculated based on bi-wing hydraulic fracture model and our previous model (Li, 2012(Li, , 2013. As shown in Figure 4, compared to bi-wing model our numerical result have a better agreement with real-world measured well producing rate.…”

“…(2) Hydraulic fracture network prediction Based on the above natural fracture distribution, hydraulic fracture network is proposed using our previous model (Li, 2012(Li, , 2013. During this process interaction between hydraulic fracture and large-scale natural fracture, turbulent flow in fracture, and influence of natural fracture on rock property (i.e., Young's modulus and Possion's ratio) are considered.…”

“…This model contains quite non-linear equations, which needs numerical method to solve and requires enough basic data, such as reservoir formation information data, logging data, fluid property, and hydraulic fracture treatment data and so on. Integrated of turbulence flow, interaction of hydraulic fracture and natural fractures, and influence of natural fracture on rock properties, one numerical model is proposed to predict fracture network propagation in shale formation (Li, 2012(Li, , 2013. This model provide one method for large scale fracture network prediction.…”

“…In this paper, we used our previous proposed model (Li, 2012(Li, , 2013 to generate hydraulic fracture network distribution in shale formation, taking into account interaction between hydraulic fracture and actual large-scale natural fractures. Integrating hydraulic fracture network results and natural fractures in non-stimulated area, highly constrained unstructured gridding and a connection list are constructed, using Mixed Discrete Fracturing Modeling (MDFM) method.…”

Diagnostics of Casing Deformation in Multi-Stage Hydraulic Fracturing Stimulation in Lower Silurian Marine Shale Play in Southwestern China

“…The highly non-linear equations in this model was solved numerically based on various data, such as formation information data, logging data, fluid property, and hydraulic fracture treatment data and so on. Combining turbulence flow, interaction of hydraulic and natural fractures, and influence of natural fractures on rock properties, Li et al proposed a model to predict fracture network propagation in shale formations (Li, 2012(Li, , 2013. This model also predicts the large scale fracture network.…”