Numerical Simulation in Hydraulic Fracturing: Multiphysics Theory and Applications


Book Description

The expansion of unconventional petroleum resources in the recent decade and the rapid development of computational technology have provided the opportunity to develop and apply 3D numerical modeling technology to simulate the hydraulic fracturing of shale and tight sand formations. This book presents 3D numerical modeling technologies for hydraulic fracturing developed in recent years, and introduces solutions to various 3D geomechanical problems related to hydraulic fracturing. In the solution processes of the case studies included in the book, fully coupled multi-physics modeling has been adopted, along with innovative computational techniques, such as submodeling. In practice, hydraulic fracturing is an essential project component in shale gas/oil development and tight sand oil, and provides an essential measure in the process of drilling cuttings reinjection (CRI). It is also an essential measure for widened mud weight window (MWW) when drilling through naturally fractured formations; the process of hydraulic plugging is a typical application of hydraulic fracturing. 3D modeling and numerical analysis of hydraulic fracturing is essential for the successful development of tight oil/gas formations: it provides accurate solutions for optimized stage intervals in a multistage fracking job. It also provides optimized well-spacing for the design of zipper-frac wells. Numerical estimation of casing integrity under stimulation injection in the hydraulic fracturing process is one of major concerns in the successful development of unconventional resources. This topic is also investigated numerically in this book. Numerical solutions to several other typical geomechanics problems related to hydraulic fracturing, such as fluid migration caused by fault reactivation and seismic activities, are also presented. This book can be used as a reference textbook to petroleum, geotechnical and geothermal engineers, to senior undergraduate, graduate and postgraduate students, and to geologists, hydrogeologists, geophysicists and applied mathematicians working in this field. This book is also a synthetic compendium of both the fundamentals and some of the most advanced aspects of hydraulic fracturing technology.







Numerical Modeling of Complex Hydraulic Fracture Propagation in Layered Reservoirs with Auto-optimization


Book Description

Hydraulic fracturing brings economic unconventional reservoir developments, and multi-cluster completion designs result in complex hydraulic fracture geometries. Therefore, accurate yet efficient modeling of the propagation of multiple non-planar hydraulic fractures is desired to study the mechanisms of hydraulic fracture propagation and optimize field completion designs. In this research, a novel hydraulic fracture model is developed to simulate the propagation of multiple hydraulic fractures with proppant transport in layered and naturally fractured reservoirs. The simplified three-dimensional displacement discontinuity method (S3D DDM) is enhanced to compute the hydraulic fracture deformation and propagation with analytical fracture height growth and vertical width variation. Using a single row of DDM elements, the enhanced S3D DDM hydraulic fracture model computes the fully 3D geometries with a similar computational intensity to a 2D model. Then an Eulerian-Lagrangian proppant transport model is developed, where the slurry flow rate and pressure are solved within the Eulerian regime, and the movement of solid proppant particles is solved within the Lagrangian regime. The adaptive proppant gridding scheme in the model allows a smaller grid size at the earlier fracturing stage for higher resolution and a larger grid size at the later fracturing stage for higher efficiency. Besides the physical model, an optimization module that utilizes advanced optimization algorithms such as genetic algorithm (GA) and pattern search algorithm (PSA) is proposed to automatically optimize the completion designs according to the preset targets. Numerical results show that hydraulic fracture propagation is under the combined influence of the in-situ stress, pumping schedule, natural fractures, and cluster placement. Hence, numerical simulation is needed to predict complex hydraulic fracture geometries under various geologic and completion settings. The complex hydraulic fracture geometries, together with fracturing fluid and proppant properties, also affect proppant placement. Moreover, the stress contrast at layer interfaces can cause proppant bridging and form barriers on the proppant transport path. The optimized completion designs increase effective hydraulic and propped areas, but they vary depending on the optimization targets. The developed hydraulic fracture model provides insights into the hydraulic fracturing process and benefits unconventional reservoir development




New numerical approaches to model hydraulic fracturing in tight reservoirs with consideration of hydro-mechanical coupling effects


Book Description

In this dissertation, two new numerical approaches for hydraulic fracturing in tight reservoir were developed. A more physical-based numerical 3D-model was developed for simulating the whole hydraulic fracturing process including fracture propagation, closure and contact as well as proppant transport and settling. In this approach rock formation, pore and fracture systems were assembled together, in which hydro-mechanical coupling effect, proppant transport and settling as well as their influences on fracture closure and contact were fully considered. A combined FDM and FVM schema was used to solve the problem. Three applications by using the new approach were presented. The results illustrated the whole hydraulic fracturing process well and seemed to be logical, which confirmed the ability of the developed approach to model the in-situ hydraulic fracturing operation from injection start till fully closure. In order to investigate the orientation problem of hydraulic fracturing in tight reservoir, a new approach for simulating arbitrary fracture propagation and orientation in 2D was developed. It was solved by a hybrid schema of XFEM and FVM. Three numerical studies were illustrated, which proved the ability of the developed approach to solve the orientation problem in field cases.







Optimization of Multistage Hydraulic Fracturing Treatment for Maximization of the Tight Gas Productivity


Book Description

Hydraulic fracturing is essential technology for the development of unconventional resources such as tight gas. So far, there are no numerical tools which can optimize the whole process from geological modeling, hydraulic fracturing until production simulation with the same 3D model with consideration of the thermo-hydro-mechanical coupling. In this dissertation, a workflow and a numerical tool chain were developed for design and optimization of multistage hydraulic fracturing in horizontal well regarding a maximum productivity of the tight gas wellbore. After the verification a full 3D reservoir model is generated based on a real tight gas field in the North German Basin. Through analysis of simulation results, a new calculation formula of FCD was proposed, which takes the proppant position and concentration into account and can predict the gas production rate more accurately. However, not only FCD but also proppant distribution and hydraulic connection of stimulated fractures to the well, geological structure and the interaction between fractures are determinant for the gas production volume. Through analysis the numerical results of sensitivity analysis and optimization variations, there is no unique criterion to determine the optimal number and spacing of the fractures, it should be analyzed firstly in detail to the actual situation and decided then from case to case.




Numerical Investigation of Water Loss Mechanisms During Hydraulic Fracturing Flow-back Operation in Tight Oil Reservoirs


Book Description

Multi-stage hydraulic fracturing is widely applied in tight reservoir exploitation. Production is enhanced significantly if hydraulic fractures can connect to regions with enhanced permeability due to the presence of micro (and induced) fractures. However, less than 50% of fracturing fluids are typically recovered. This study models the mechanisms of water loss and retention in fracture-matrix system. The effects of capillarity and geomechanics are systematically investigated, and the time scale of water imbibition under different reservoir conditions is tested. During the shut-in (soaking) and flow-back periods, the fracture conductivity decreases as effective stress increases due to imbibition. Previous works have addressed fracture closure during the production phase; however, the coupling of imbibition due to multiphase flow and stress-dependent fracture properties during shut-in is less understood. Numerical simulation results indicate the circumstances under which this phenomenon might be beneficial or detrimental to subsequent on tight oil production. A series of mechanistic simulation models consisting of both hydraulic fractures and stochastically distributed micro fractures are constructed to simulate fluid distribution during shut-in and flow-back. Three systems: matrix, hydraulic fracture and micro fractures are explicitly represented in the computational domain. Fluid loss and retention mechanisms are systematically investigated in this study by subjecting mechanistic model to different reservoir conditions. Water imbibition into the matrix would help to displace hydrocarbons into nearby micro and hydraulic fractures, and this process could lead to an increase in initial rate. Although other water loss mechanisms including water loss in desiccated matrix and water trapping in induced micro fractures were proposed in literature, detailed understanding of the roles of water trapping in these systems is still lacking. Impacts of secondary fracture distributions and properties, matrix permeability, multiphase flow functions, wettability, initial saturation, water injection rate and shut-in duration on fluid retention and the associated time scales are assessed. Increase in short-term oil production as a result of imbibition could be counteracted by the reduction in flow capability due to fracture closure. Therefore, the coupling of stress-dependent fracture conductivity and imbibition are studied next. Our results indicate that fracture compaction can enhance imbibition and water loss, which in turn leads to further reduction in fracture pressure and conductivity. Spatial variability in micro-fracture properties is also modeled probabilistically to investigate whether it is possible for fracturing fluid to be trapped in the micro fractures, or conversely, the micro fractures could provide alternate pathways for fluids to access the hydraulic fracture systems. This work presents a quantitative study of the controlling factors of water retention due to fluid-rock properties and geomechanics. It investigates the roles of multi-scale fractures in flow-back behavior and ensuing recovery performance. The results highlight 1) the crucial interplay between shut-in duration and properties of connected fractures in short- and long-term production performances; 2) the critical interaction between imbibition and geomechanics in short- and long-term production performances. The results would have considerable impacts on understanding and improving current industry practice on fracturing design and assessment of stimulated reservoir volume.







Numerical Investigation of Interaction Between Hydraulic Fractures and Natural Fractures


Book Description

Hydraulic fracturing of a naturally-fractured reservoir is a challenge for industry, as fractures can have complex growth patterns when propagating in systems of natural fractures in the reservoir. Fracture propagation near a natural fracture (NF) considering interaction between a hydraulic fracture (HF) and a pre-existing NF, has been investigated comprehensively using a two dimensional Displacement Discontinuity Method (DDM) Model in this thesis. The rock is first considered as an elastic impermeable medium (with no leakoff), and then the effects of pore pressure change as a result of leakoff of fracturing fluid are considered. A uniform pressure fluid model and a Newtonian fluid flow model are used to calculate the fluid flow, fluid pressure and width distribution along the fracture. Joint elements are implemented to describe different NF contact modes (stick, slip, and open mode). The structural criterion is used for predicting the direction and mode of fracture propagation. The numerical model was used to first examine the mechanical response of the NF to predict potential reactivation of the NF and the resultant probable location for fracture re-initiation. Results demonstrate that: 1) Before the HF reaches a NF, the possibility of fracture re-initiation across the NF and with an offset is enhanced when the NF has weaker interfaces; 2) During the stage of fluid infiltration along the NF, a maximum tensile stress peak can be generated at the end of the opening zone along the NF ahead of the fluid front; 3) Poroelastic effects, arising from fluid diffusion into the rock deformation can induce closure and compressive stress at the center of the NF ahead of the HF tip before HF arrival. Upon coalescence when fluid flows along the NF, the poroelastic effects tend to reduce the value of the HF aperture and this decreases the tension peak and the possibility of fracture re-initiation with time. Next, HF trajectories near a NF were examined prior to coalesce with the NF using different joint, rock and fluid properties. Our analysis shows that: 1) Hydraulic fracture trajectories near a NF may bend and deviate from the direction of the maximum horizontal stress when using a joint model that includes initial joint deformation; 2) Hydraulic fractures propagating with higher injection rate or fracturing fluid of higher viscosity propagate longer distance when turning to the direction of maximum horizontal stress; 3) Fracture trajectories are less dependent on injection rate or fluid viscosity when using a joint model that includes initial joint deformation; whereas, they are more dominated by injection rate and fluid viscosity when using a joint model that excludes initial joint deformation.