Search

Abstract In this paper, we investigate the effects of poroelastic properties and permeability on hydraulic fracturing in an anisotropic medium. A hybrid approach combining the EDFM (embedded discrete fracture model) and XFEM (extended finite element method) is proposed. Fractures are embedded into a nonconforming grid used to determine the fluid flow in a porous matrix by the mimetic finite difference method, accounting for an anisotropic permeability tensor. The stress-strain problem is solved by the extended finite element method with the same grid. The fluid flow and rock deformation as well as the fracture propagation are iteratively coupled. The proposed approach is validated against the analytical solutions for Mandel's problem and the KDG model. A series of calculations are performed, and the obtained results are analyzed to investigate the effects of the anisotropy of permeability, that of the elastic modulus and Biot's coefficient on the hydraulic fracturing process. It is found that the anisotropy of permeability has a significant influence on the geometrical parameter of a fracture, while the anisotropy of the elastic modulus has a dominating influence on the propagation direction of a fracture. Biot's coefficient also has an influence on the fracture propagation kinetics. [ABSTRACT FROM AUTHOR]

In view of fluid flow in cracks and rock deformation, the mathematical model for shale hydraulic crack propagation was established. The crack flow field and the rock stress field were solved with the finite element method and the extended finite element method respectively, and the two fields were coupled through the Picard iteration. The presented model gave results consistent with those of the classic model, which verified correctness of the former. Based on the model, the effects of the rock elasticity modulus, Poisson's ratio and injection rate on the crack geometry, and the dynamic process of a hydraulic crack approaching a natural crack at an arbitrary angle, were simulated. The numerical results show that the elasticity modulus and injection rate have significant influence on the crack geometry, while the Poisson's ratio has little effect; the more brittle the shale is, the longer and narrower the hydraulic crack will grow; the greater the principal stress difference and the approaching angle are, the easier the hydraulic crack crosses the natural crack; there is a relatively large decrease in the crack width at the intersection between a hydraulic crack and a natural crack; the extended finite element method avoids mesh reconstruction and refinement during computation, and reduces the computing time. The presented model provides an effective theoretical tool for the shale fracturing design. [ABSTRACT FROM AUTHOR]

This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing. The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method (DDM) and the Picard iterative method. The shale gas flow considers multiple transport mechanisms, and the flow in the fracture network is handled by the embedded discrete fracture model (EDFM). A series of numerical simulations are conducted to analyze the effects of the cluster number, stage spacing, stress difference coefficient, and natural fracture distribution on the stimulated fracture area, fractal dimension, and cumulative gas production, and their correlation coefficients are obtained. The results show that the most influential factors to the stimulated fracture area are the stress difference ratio, stage spacing, and natural fracture density, while those to the cumulative gas production are the stress difference ratio, natural fracture density, and cluster number. This indicates that the stress condition dominates the gas production, and employing intensive volume fracturing (by properly increasing the cluster number) is beneficial for improving the final cumulative gas production. [ABSTRACT FROM AUTHOR]

Rock yielding may well take place during hydraulic fracturing in deep reservoirs. The prevailing models based on the linear elastic fracture mechanics (LEFM) are incapable of describing the evolution process of hydraulic fractures accurately. In this paper, a hydro-elasto-plastic model is proposed to investigate the hydraulic fracture propagation in deep reservoirs. The Drucker–Prager plasticity model, Darcy's law, cubic law and cohesive zone model are employed to describe the plastic deformation, matrix flow, fracture flow and evolution of hydraulic fractures, respectively. Combining the embedded discrete fracture model (EDFM), extended finite element method (XFEM) and finite volume method, a hybrid numerical scheme is presented to carry out simulations. A dual-layer iterative procedure is developed based on the fixed-stress split method, Picard iterative method and Newton–Raphson iterative method. The iterative procedure is used to deal with the coupling between nonlinear deformation with fracture extension and fluid flow. The proposed model is verified against analytical solutions and other numerical simulation results. A series of numerical cases are performed to investigate the influences of rock plasticity, internal friction angle, dilatancy angle and permeability on hydraulic fracture propagation. Finally, the proposed model is extended to simulate multiple hydraulic fracture propagation. The result shows that plastic deformation can enhance the stress-shadowing effect. [ABSTRACT FROM AUTHOR]

• A 3D hydro-mechanical coupled model is presented to study hydraulic fracturing interaction with bedding planes. • Cohesive zone method is employed to represent the bedding planes explicitly. • Fracture propagation behaviors are analyzed under different patterns. Bedding planes are abundant in the continental shale oil reservoirs, but their impact on hydraulic fracturing is still not fully understood. To investigate the effects of stress condition, bedding plane's properties and layer lithology on hydraulic fracture propagation in layered formation, a 3D hydro-mechanical coupled model is presented. The bedding planes are explicitly represented by cohesive elements with degraded properties compared to the rock. The results show that increasing the stress difference ratio of the vertical principal stress to horizontal principal stress promotes the hydraulic fracture propagating across bedding planes, and the critical value of vertical principal stress required for hydraulic fracture crossing rises with the decrease of bedding plane's critical traction. Furthermore, when the critical traction ratio of bedding plane to rock drops below 0.7, the fracture propagation pattern changes from crossing to deflecting. Additionally, the extension pressure increases slightly, the fracture length decreases and the area of opening bedding planes rises. As the permeability of bedding plane increases, the fracture cannot cross the bedding planes and even stops extension. The layer lithology contrast also affects fracture propagation patterns. This study examines the fracture propagation behaviors under different patterns, including the opening morphology of bedding plane and the evolution of injection pressure and strain energy. The simulation results may provide useful insights for hydraulic fracturing design in layered formations. [ABSTRACT FROM AUTHOR]

Fracture penetration or deflection may occur when the hydraulic fracture approaches weak interface in multi-layered formation. However, the intrinsic mechanism of fracture propagation pattern remains unclear. A model was presented to investigate hydraulic fracture propagation in multi-layered formation via the phase field method. The interface is explicitly represented by a weak thin layer. Fluid flow in the fractured media was modeled by using the indicator function of phase field. The iterative algorithm was described and validated against the benchmark cases. Simulation results show that the evolution of injection pressure is more complex when the fracture deflects than that when the fracture penetrates the interface. The fracture is prone to penetrate the softer layer the perforated layer while deflects to the interface if the adjacent layer is stiffer. Decreasing the modulus ratio of interface to layers promotes the fracture deflecting to the interface, while increasing the stress ratio contributes to the fracture penetrating the interface. The stacking mode of layers greatly affects the fracture propagation morphology in three or more layered formations. Rock heterogeneity doesn't alter the fracture propagation pattern as a whole, but it affects the local fracture morphology. The results provide a cloud chart of the relationship between fracture propagation pattern and controlling factors for the potential of field application. [ABSTRACT FROM AUTHOR]