Although many fracture models are based on two-dimensional plane strain approximations, accurately predicting fracture propagation geometry requires accounting for the three-dimensional aspects of fractures. In this study, we implemented 3-D displacement discontinuity (DD) boundary element modeling to investigate the following intrinsically 3-D natural or hydraulic fracture propagation problems: the effect of fracture height on lateral propagation of vertical natural fractures, joint development in the vicinity of normal faults, and hydraulic fracture height growth and non-planar propagation paths. Fracture propagation is controlled by stress intensity factor (SIF) and its determination plays a central role in LEFM. The DD modeling is used to evaluate SIF in Mode I, II and III at the tip of an arbitrarily-shaped embedded crack by using crack-tip element displacement discontinuity. We examine the accuracy of SIF calculation is for rectangular, penny-shaped, and elliptical planar cracks. Using the aforementioned model for lateral propagation of overlapping fractures shows that the curving path of overlapping fractures is strongly influenced by the spacing-to-height ratio of fractures, as well as the differential stress magnitude. We show that the angle of intersection between two non-coincident but parallel en-echelon fractures depends strongly on the fracture height-to-spacing ratio, with intersection angles being asymptotic for "tall" fractures (large height-to-spacing ratios) and nearly orthogonal for "short" fractures. Stress perturbation around normal faults is three-dimensionally heterogeneous. That perturbation can result in joint development at the vicinity of normal faults. We examine the geometrical relationship between genetically related normal faults and joints in various geologic environments by considering a published case study of fault-related joints in the Arches National Park region, Utah. The results show that joint orientation is dependent on vertical position with respect to the normal fault, the spacing-to-height ratio of sub-parallel normal faults, and Poisson's ratio of the media. Our calculations represent a more physically reasonable match to measured field data than previously published, and we also identify a new mechanism to explain the driving stress for opening mode fracture propagation upon burial of quasi-elastic rocks. Hydraulic fractures may not necessarily start perpendicular to the minimum horizontal remote stress. We use the developed fracture propagation model to explain abnormality in the geometry of fracturing from misaligned horizontal wellbores. Results show that the misalignment causes non-planar lateral propagation and restriction in fracture height and fracture width in wellbore part. / text

Successful creations of multiple hydraulic fractures in horizontal wells are critical for economic development of unconventional reservoirs. The recent advances in diagnostic techniques suggest that multi-fracturing stimulation in unconventional reservoirs has often caused complex fracture geometry. The most important factors that might be responsible for the fracture complexity are fracture interaction and the intersection of the hydraulic and natural fracture. The complexity of fracture geometry results in significant uncertainty in fracturing treatment designs and production optimization. Modeling complex fracture propagation can provide a vital link between fracture geometry and stimulation treatments and play a significant role in economically developing unconventional reservoirs. In this research, a novel fracture propagation model was developed to simulate complex hydraulic fracture propagation in unconventional reservoirs. The model coupled rock deformation with fluid flow in the fractures and the horizontal wellbore. A Simplified Three Dimensional Displacement Discontinuity Method (S3D DDM) was proposed to describe rock deformation, calculating fracture opening and shearing as well as fracture interaction. This simplified 3D method is much more accurate than faster pseudo-3D methods for describing multiple fracture propagation but requires significantly less computational effort than fully three-dimensional methods. The mechanical interaction can enhance opening or induce closing of certain crack elements or non-planar propagation. Fluid flow in the fracture and the associated pressure drop were based on the lubrication theory. Fluid flow in the horizontal wellbore was treated as an electrical circuit network to compute the partition of flow rate between multiple fractures and maintain pressure compatibility between the horizontal wellbore and multiple fractures. Iteratively and fully coupled procedures were employed to couple rock deformation and fluid flow by the Newton-Raphson method and the Picard iteration method. The numerical model was applied to understand physical mechanisms of complex fracture geometry and offer insights for operators to design fracturing treatments and optimize the production. Modeling results suggested that non-planar fracture geometry could be generated by an initial fracture with an angle deviating from the direction of the maximum horizontal stress, or by multiple fracture propagation in closed spacing. Stress shadow effects are induced by opening fractures and affect multiple fracture propagation. For closely spaced multiple fractures growing simultaneously, width of the interior fractures are usually significantly restricted, and length of the exterior fractures are much longer than that of the interior fractures. The exterior fractures receive most of fluid and dominate propagation, resulting in immature development of the interior fractures. Natural fractures could further complicate fracture geometry. When a hydraulic fracture encounters a natural fracture and propagates along the pre-existing path of the natural fracture, fracture width on the natural fracture segment will be restricted and injection pressure will increase, as a result of stress shadow effects from hydraulic fracture segments and additional closing stresses from in-situ stress field. When multiple fractures propagate in naturally fracture reservoirs, complex fracture networks could be induced, which are affected by perforation cluster spacing, differential stress and natural fracture patterns. Combination of our numerical model and diagnostic methods (e.g. Microseismicity, DTS and DAS) is an effective approach to accurately characterize the complex fracture geometry. Furthermore, the physics-based complex fracture geometry provided by our model can be imported into reservoir simulation models for production analysis. / text

Large amounts of hydrocarbon reserves are trapped in fractured reservoirs where fluid flux is far more rapid along fractures than through the porous matrix, even though the volume of the pore space may be a hundred times greater than the volume of the fractures. These are considered extremely challenging in terms of accurate recovery prediction because of their complexity and heterogeneity. Conventional reservoir simulators are generally not suited to naturally fractured reservoirs’ production history simulation, especially when production processes are associated with large pressure and temperature changes that lead to large redistribution of effective stresses, causing natural fracture aperture alterations. In this case, all the effective processes, i.e. hydraulic, thermal and geomechanical, should be considered simultaneously to explain and evaluate the behavior of stress-sensitive reservoirs over the production period. This is called thermo-hydro-mechanical (THM) coupling. In this study, a fully coupled thermo-hydro-mechanical approach is developed to simulate the physical behavior of fractures in a plane strain thermo-poroelastic medium. A hybrid numerical method, which implements both the finite difference method (FDM) and the displacement discontinuity method (DDM), is established to study the pressure, temperature, deformation and stress variations of fractures and surrounding rocks during production processes. This method is straightforward and can be implemented in conventional reservoir simulators to update fracture conductivity as it uses the same grid block as the reservoir grids and requires only discretization of fractures. The hybrid model is then verified with couple of analytical solutions for the fracture aperture variation under different conditions. This model is implemented for some examples to present the behavior of fracture network as well as its surrounding rock under thermal injection and production. The results of this work clearly show the importance of rate, aspect ratio (i.e. geometry) and the coupling effects among fracture flow rate and aperture changes arising from coupled stress, pressure and temperature changes. The outcomes of this approach can be used to study the behavior of hydraulic injection for induced fracturing and promoting of shearing such as hydraulic fracturing of shale gas or shale oil reservoirs as well as massive waste disposal in the porous carbonate rocks. Furthermore, implementation of this technique should be able to lead to a better understanding of induced seismicity in injection projects of all kinds, whether it is for waste water disposal, or for the extraction of geothermal energy.

There are two big challenges which restrict the extensive application of fully coupled geomechanics-reservoir modeling. The first challenge is computational effort. Consider a 3-D simulation combining pressure and heat diffusion, elastoplastic mechanical response, and saturation changes; each node has at least 5 degrees of freedom, each leading to a separate equation. Furthermore, regions of large p, T and σ′ gradients require small-scale discretization for accurate solutions, greatly increasing the number of equations. When the rock mass surrounding the reservoir region is included, it is represented by many elements or nodes. These factors mean that accurate analysis of realistic 3-D problems is challenging, and will so remain as we seek to solve larger and larger coupled problems involving nonlinear responses. To overcome the first challenge, the displacement discontinuity method is introduced wherein a large-scale 3-D case is divided into a reservoir region where Δp, ΔT and non-linear effects are critical and analyzed using FEM, and an outside region in which the reservoir is encased where Δp and ΔT effects are inconsequential and the rock may be treated as elastic, analyzed with a 3D displacement discontinuity formulation. This scheme leads to a tremendous reduction in the degrees of freedom, yet allows for reasonably rigorous incorporation of the reactions of the surrounding rock. The second challenge arises from some forms of numerical instability. There are actually two types of sharp gradients implied in the transient advection-diffusion problem: one is caused by the high Peclet numbers, the other by the sharp gradient which appears during the small time steps due to the transient solution. The way to eliminate the spurious oscillations is different when the sharp gradients are induced by the transient evolution than when they are produced by the advective terms, and existing literature focuses mainly on eliminating the spurious spatial temperature oscillations caused by advection-dominated flow. To overcome the second challenge, numerical instability sources are addressed by introducing a new stabilized finite element method, the subgrid scale/gradient subgrid scale (SGS/GSGS) method.

There are two big challenges which restrict the extensive application of fully coupled geomechanics-reservoir modeling. The first challenge is computational effort. Consider a 3-D simulation combining pressure and heat diffusion, elastoplastic mechanical response, and saturation changes; each node has at least 5 degrees of freedom, each leading to a separate equation. Furthermore, regions of large p, T and σ′ gradients require small-scale discretization for accurate solutions, greatly increasing the number of equations. When the rock mass surrounding the reservoir region is included, it is represented by many elements or nodes. These factors mean that accurate analysis of realistic 3-D problems is challenging, and will so remain as we seek to solve larger and larger coupled problems involving nonlinear responses. To overcome the first challenge, the displacement discontinuity method is introduced wherein a large-scale 3-D case is divided into a reservoir region where Δp, ΔT and non-linear effects are critical and analyzed using FEM, and an outside region in which the reservoir is encased where Δp and ΔT effects are inconsequential and the rock may be treated as elastic, analyzed with a 3D displacement discontinuity formulation. This scheme leads to a tremendous reduction in the degrees of freedom, yet allows for reasonably rigorous incorporation of the reactions of the surrounding rock. The second challenge arises from some forms of numerical instability. There are actually two types of sharp gradients implied in the transient advection-diffusion problem: one is caused by the high Peclet numbers, the other by the sharp gradient which appears during the small time steps due to the transient solution. The way to eliminate the spurious oscillations is different when the sharp gradients are induced by the transient evolution than when they are produced by the advective terms, and existing literature focuses mainly on eliminating the spurious spatial temperature oscillations caused by advection-dominated flow. To overcome the second challenge, numerical instability sources are addressed by introducing a new stabilized finite element method, the subgrid scale/gradient subgrid scale (SGS/GSGS) method.

In geothermal reservoirs and unconventional gas reservoirs with very low matrix permeability, fractures are the main routes of fluid flow and heat transport, so the fracture permeability change is important. In fact, reservoir development under this circumstance relies on generation and stimulation of a fracture network. This thesis presents numerical simulation of the response of a fractured rock to injection and extraction considering the role of poro-thermoelasticity and joint deformation. Fluid flow and heat transport in the fracture are treated using a finite difference method while the fracture and rock matrix deformation are determined using the displacement discontinuity method (DDM). The fractures response to fluid injection and extraction is affected both by the induced stresses as well as by the initial far-field stress. The latter is accounted for using the non-equilibrium condition, i.e., relaxing the assumption that the rock joints are in equilibrium with the in-situ stress state. The fully coupled DDM simulation has been used to carry out several case studies to model the fracture response under different injection/extractions, in-situ stresses, joint geometries and properties, for both equilibrium and non-equilibrium conditions. The following observations are made: i) Fluid injection increases the pressure causing the joint to open. For non-isothermal injection, cooling increases the fracture aperture drastically by inducing tensile stresses. Higher fracture aperture means higher conductivity. ii) In a single fracture under constant anisotropic in-situ stress (non-equilibrium condition), permanent shear slip is encountered on all fracture segments when the shear strength is overcome by shear stress in response to fluid injection. With cooling operation, the fracture segments in the vicinity of the injection point are opened due to cooling-induced tensile stress and injection pressure, and all the fracture segments experience slip. iii) Fluid pressure in fractures increases in response to compression. The fluid compressibility and joint stiffness play a role. iv) When there are injection and extraction in fractured reservoirs, the cooler fluid flows through the fracture channels from the injection point to extraction well extracting heat from the warmer reservoir matrix. As the matrix cools, the resulting thermal stress increases the fracture apertures and thus increases the fracture conductivity. v) Injection decreases the amount of effective stress due to pressure increase in fracture and matrix near a well. In contrast, extraction increases the amount of effective stress due to pressure drop in fracture and matrix.

No description available.