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]

As an effective approach to evaluate shale gas production, macroscale numerical simulations have been conducted; however, multiscale storage space and complex gas–water two-phase transport mechanisms are not systematically considered. To deal with this, a rigorous cross-scale gas–water two-phase transport simulation method that considers the pore, core, and field scale is proposed herein. First, to investigate the fluid storage state in different pore types during the gas–water two-phase transport process, a shale pore network gas–water transport simulation is conducted. Subsequently, a core-scale model that considers the distribution of organic matter is constructed, and the gas–water two-phase transport behaviors are upscaled from the pore to the core by incorporating nanomicroscale effects. Next, a shale gas–water two-phase macroscale numerical simulation is implemented based on the core-scale simulation results. Multiple interacting continua are used to describe the shale matrix and microfractures, whereas the embedded discrete fracture method is employed to represent hydraulic fractures. Finally, the manner by which nanomicroscale effects, stress sensitivity, and SRV area affect the production of shale gas reservoirs is discussed comprehensively. This study provides a practical method to estimate shale gas production by rigorously considering gas–water nanomicroscale effects to reduce uncertainties during productivity evaluation. [ABSTRACT FROM AUTHOR]

In this study, we present an efficient numerical solution for studying hydraulic fracturing under coupled thermal–hydraulic conditions in an elastic–plastic porous medium. The propagation of macroscopic fracture is described by using an extended finite element method. Both the fluid flow through the porous medium and the exchange between the medium and fracture are taken into account. It is the same for the heat transfer. An efficient iterative scheme is then proposed to deal with the coupling between material deformation with fracture growth, fluid flow and heat transfer. The proposed method is assessed through comparisons with analytical solutions for a number of well-established problems. A series of numerical calculations are further performed in order to investigate the effect of plastic deformation and temperature change on the process of hydraulic fracture propagation. [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]

In order to investigate pressure performance of multiple fractured horizontal wells (MFHWs) penetrating heterogeneous unconventional reservoir and avoid the high computational cost of numerical simulation, a semi-analytical model for MFHWs combining Green function solution and boundary element method has been obtained, where the reservoir is divided into different homogeneous substructures and coupled at interface boundaries by plane source function in a closed rectangular parallelepiped. Hydraulic fractures are assumed uniform flux and dual porosity model is used for natural fractures system. Then the model is validated by compared with analytical solution of MFHWs in a homogeneous reservoir and trilinear flow model, which shows that this model can achieve high accuracy even with a small interface discretization number, and it can consider the radial flow around each hydraulic fractures. Finally, the pressure responses with heterogeneous parameters of reservoirs are discussed including heterogeneous permeability, non-uniform block-length and fracture half-length distribution as well as dual porosity parameters like elastic storage ratio and crossflow ratio. [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]

Nonuniform proppant distribution is fairly common in hydraulic fractures, and different closure behaviors of the propped and unpropped fractures have been observed in lots of physical experiments. However, the modeling of partially propped fracture closure is rarely performed, and its effect on gas production is not well understood as a result of previous studies. In this paper, a fully coupled fluid flow and geomechanics model is developed to simulate partially propped fracture closure, and to examine its effect on gas production in fractured shale reservoirs. Specifically, an efficient hybrid model, which consists of a single porosity model, a multiple porosity model and the embedded discrete fracture model (EDFM), is adopted to model the hydro-mechanical coupling process in fractured shale reservoirs. In flow equations, the Klinkenberg effect is considered in gas apparent permeability, and adsorption/desorption is treated as an additional source term. In the geomechanical domain, the closure behaviors of propped and unpropped fractures are described through two different constitutive models. Then, a stabilized extended finite element method (XFEM) iterative formulation, which is based on the polynomial pressure projection (PPP) technique, is developed to simulate a partially propped fracture closure with the consideration of displacement discontinuity at the fracture interfaces. After that, the sequential implicit method is applied to solve the coupled problem, in which the finite volume method (FVM) and stabilized XFEM are applied to discretize the flow and geomechanics equations, respectively. Finally, the proposed method is validated through some numerical examples, and then it is further used to study the effect of partially propped fracture closures on gas production in 3D fractured shale reservoir simulation models. This work will contribute to a better understanding of the dynamic behaviors of fractured shale reservoirs during gas production, and will provide more realistic production forecasts. [ABSTRACT FROM AUTHOR]