The distributions of pore pressure and stresses around a fracture are of interest in conventional hydraulic fracturing operations, fracturing during water-flooding of petroleum reservoirs, shale gas, and injection/extraction operations in a geothermal reservoir. During the operations, the pore pressure will increase with fluid injection into the fracture and leak off to surround the formation. The pore pressure increase will induce the stress variations around the fracture surface. This can cause the slip of weakness planes in the formation and cause the variation of the permeability in the reservoir. Therefore, the investigation on the pore pressure and stress variations around a hydraulic fracture in petroleum and geothermal reservoirs has practical applications. The stress and pore pressure fields around a fracture are affected by: poroelastic, thermoelastic phenomena as well as by fracture opening under the combined action of applied pressure and in-situ stress. In our study, we built up two models. One is a model (WFPSD model) of water-flood induced fracturing from a single well in an infinite reservoir. WFPSD model calculates the length of a water flood fracture and the extent of the cooled and flooded zones. The second model (FracJStim model) calculates the stress and pore pressure distribution around a fracture of a given length under the action of applied internal pressure and in-situ stresses as well as their variation due to cooling and pore pressure changes. In our FracJStim model, the Structural Permeability Diagram is used to estimate the required additional pore pressure to reactivate the joints in the rock formations of the reservoir. By estimating the failed reservoir volume and comparing with the actual stimulated reservoir volume, the enhanced reservoir permeability in the stimulated zone can be estimated. In our research, the traditional two dimensional hydraulic fracturing propagation models are reviewed, the propagation and recession of a poroelastic PKN hydraulic fracturing model are studied, and the pore pressure and stress distributions around a hydraulically induced fracture are calculated and plotted at a specific time. The pore pressure and stress distributions are used to estimate the failure potentials of the joints in rock formations around the hydraulic fracture. The joint slips and rock failure result in permeability change which can be calculated under certain conditions. As a case study and verification step, the failure of rock mass around a hydraulic fracture for the stimulation of Barnett Shale is considered. With the simulations using our models, the pore pressure and poro-induced stresses around a hydraulic fracture are elliptically distributed near the fracture. From the case study on Barnett Shale, the required additional pore pressure is about 0.06 psi/ft. With the given treatment pressure, the enhanced permeability after the stimulation of hydraulic fracture is calculated and plotted. And the results can be verified by previous work by Palmer, Moschovidis and Cameron in 2007.

Multiple transverse fracturing along a horizontal well is a relatively new technology that is designed to increase well productivity by increasing the contact between the reservoir and the wellbore. For multiple transverse fractures, the performance of the well system is determined by three aspects: the inflow from the reservoir to the fracture, the flow from the fracture to the wellbore, and the inflow from the reservoir to the horizontal wellbore. These three aspects influence each other and combined, influence the wellbore outflow. In this study, we develop a model to effectively formulate the inter-relationships of a multi-fracture system. This model includes a reservoir model and a wellbore model. The reservoir model is established to calculate both independent and inter-fracture productivity index to quantify the contribution from all fractures on pressure drop of each fracture, by using the source functions to solve the single-phase gas reservoir flow model. The wellbore model is used to calculate the pressure distribution along the wellbore and the relationship of pressure between neighboring fractures, based on the basic pressure drop model derived from the mechanical energy balance. A set of equations with exactly the same number of fractures will be formed to model the system by integrating the two models. Because the equations are nonlinear, iteration method is used to solve them. With our integrated reservoir and wellbore model, we conduct a field study to find the best strategy to develop the field by hydraulic fracturing. The influence of reservoir size, horizontal and vertical permeability, well placement, and fracture orientation, type (longitudinal and transverse), number and distribution are completely examined in this study. For any specific field, a rigorous step-by-step procedure is proposed to optimize the field.

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.

This thesis presents a methodology based on Proppant Number approach for optimal fracture treatment design of natural gas wells considering non-Darcy flow effects in the design process. Closure stress is taken into account, by default, because it is the first factor decreasing propped pack permeability at in-situ conditions. Gel damage was also considered in order to evaluate the impact of incorporating more damaging factors on ultimate well performance and optimal geometry. Effective fracture permeability and optimal fracture geometry are calculated through an iterative process. This approach was implemented in a spreadsheet. Non-Darcy flow is described by the β factor. All β factor correlations available in the literature were evaluated. It is recommended to use the correlation developed specifically for the given type of proppant and mesh size, if available. Otherwise, the Pursell et al. or the Martins et al. equations are recommended as across the board reliable correlations for predicting non-Darcy flow effects in the propped pack. The proposed methodology was implemented in the design of 11 fracture treatments of 3 natural tight gas wells in South Texas. Results show that optimal fracture design might increase expected production in 9.64 MMscf with respect to design that assumes Darcy flow through the propped pack. The basic finding is that for a given amount of proppant shorter and wider fractures compensate the non-Darcy and/or gel damage effect. Dynamic programming technique was implemented in design of multistage fractures for one of the wells under study for maximizing total gas production. Results show it is a powerful and simple technique for this application. It is recommended to expand its use in multistage fracture designs.

The complex interrelations among the different physical processes involved in acid fracturing make it difficult to design, and later, to predict the outcome of stimulation jobs. Actual tendencies require the use of computational models to deal with the dynamic interaction of variables. This thesis presents a new study of acidized surface textures by means of a laser profilometer to improve our understanding of the remaining etched surface topography and its hydraulic response. Visualization plots generated by the profilometer identified hydrodynamic channels that could not be identified by the naked eye in acidized surfaces. The plots clarified the existence of rock heterogeneities and revealed how the processes of dissolution function in chalk rock. Experimental data showed clearly that the effect of dissolution depends on the type of rock and the fluid system; dolomite, for example, dissolves more rapidly but more roughly than limestone. Fluid leakoff rate and temperature also affect the dissolution. Further research is necessary to clarify the effects of conductivity.

Hydraulic fracturing of shale has combined traditional oil and gas industry techniques to create significant new reserves in the United States. Poor science, incomplete media coverage and politicization of the issues threaten broad understanding of issues of genuine concern while overstating others. The Environmental Protection Agency should focus on science-based regulation prior to enumerating new rules and should continue to cede primacy to the states where traditional regimes have proven successful in regulating oil and gas. The most critical issues associated with hydraulic fracturing tend to be regional and predicated on local hydrogeology. Surface water disposal and emissions standards need revision and strengthening. Scarce resources should be dedicated to better understanding regional water availability and to heightened awareness of the energy-water nexus. / text

Large volumes of natural gas exist in tight fissured reservoirs. Hydraulic fracturing is one of the main stimulating techniques to enhance recovery from these fractured reservoirs. Although hydraulic fracturing has been used for decades for the stimulation of tight gas reservoirs, a thorough understanding of the interaction between induced hydraulic fractures and natural fractures is still lacking. Recent examples of hydraulic fracture diagnostic data suggest complex, multi-stranded hydraulic fracture geometry is a common occurrence. The interaction between pre-existing natural fractures and the advancing hydraulic fracture is a key condition leading to complex fracture patterns. Large populations of natural fractures that exist in formations such as the Barnett shale are sealed by precipitated cements which could be quartz, calcite, etc. Even though there is no porosity in the sealed fractures, they may still serve as weak paths for fracture initiation and/or for diverting the path of the growing hydraulic fractures. Performing hydraulic fracture design calculations under these complex conditions requires modeling of fracture intersections and tracking fluid fronts in the network of reactivated fissures. In this dissertation, the effect of the cohesiveness of the sealed natural fractures and the intact rock toughness in hydraulic fracturing are studied. Accordingly, the role of the pre-existing fracture geometry is also investigated. The results provide some explanations for significant differences in hydraulic fracturing in naturally fractured reservoirs from non-fractured reservoirs. For the purpose of this research, an extended finite element method (XFEM) code is developed to simulate fracture propagation, initiation and intersection. The motivation behind applying XFEM are the desire to avoid remeshing in each step of the fracture propagation, being able to consider arbitrary varying geometry of natural fractures and the insensitivity of fracture propagation to mesh geometry. New modifications are introduced into XFEM to improve stress intensity factor calculations, including fracture intersection criteria into the model and improving accuracy of the solution in near crack tip regions. The presented coupled fluid flow-fracture mechanics simulations extend available modeling efforts and provide a unified framework for evaluating fracture design parameters and their consequences. Results demonstrate that fracture pattern complexity is strongly controlled by the magnitude of in situ stress anisotropy, the rock toughness, the natural fracture cement strength, and the approach angle of the hydraulic fracture to the natural fracture. Previous studies (mostly based on frictional fault stability analysis) have concentrated on predicting the onset of natural fracture failure. However, the use of fracture mechanics and XFEM makes it possible to evaluate the progression of fracture growth over time as fluid is diverted into the natural fractures. Analysis shows that the growing hydraulic fracture may exert enough tensile and/or shear stresses on cemented natural fractures that they may be opened or slip in advance of hydraulic fracture tip arrival, while under some conditions, natural fractures will be unaffected by the hydraulic fracture. A threshold is defined for the fracture energy of cements where, for cases below this threshold, hydraulic fractures divert into the natural fractures. The value of this threshold is calculated for different fracture set orientations. Finally, detailed pressure profile and aperture distributions at the intersection between fracture segments show the potential for difficulty in proppant transport under complex fracture propagation conditions. Whether a hydraulic fracture crosses or is arrested by a pre-existing natural fracture is controlled by shear strength and potential slippage at the fracture intersections, as well as potential debonding of sealed cracks in the near-tip region of a propagating hydraulic fracture. We introduce a new more general criterion for fracture propagation at the intersections. We present a complex hydraulic fracture pattern propagation model based on the Extended Finite Element Method as a design tool that can be used to optimize treatment parameters under complex propagation conditions. / text

Micro seismic data and coring studies suggest that hydraulic fractures interact heavily with natural fractures creating complex fracture networks in naturally fractured reservoirs such as the Barnett shale, the Eagle Ford shale, and the Marcellus shale. However, since direct observations of subsurface hydraulic fracture geometries are incomplete or nonexistent, we look to properly scaled experimental research and computer modeling based on realistic assumptions to help us understand fracture intersection geometries. Most experimental analysis of this problem has focused on natural fractures with frictional interfaces. However, core observations from the Barnett and other shale plays suggest that natural fractures are largely cemented. To examine hydraulic fracture interactions with cemented natural fractures, we performed 9 hydraulic fracturing experiments in gypsum cement blocks that contained embedded planar glass, sandstone, and plaster discontinuities which acted as proxies for cemented natural fractures. There were three main fracture intersection geometries observed in our experimental program. 1) A hydraulic fracture is diverted into a different propagation path(s) along a natural fracture. 2) A taller hydraulic fracture bypasses a shorter natural fracture by propagating around it via height growth while also separating the weakly bonded interface between the natural fracture and the host rock. 3) A hydraulic fracture bypasses a natural fracture and also diverts down it to form separate fractures. The three main factors that seemed to have the strongest influence on fracture intersection geometry were the angle of intersection, the ratio of hydraulic fracture height to natural fracture height, and the differential stress. Our results show that bypass, separation of weakly bonded interfaces, diversion, and mixed mode propagation are likely in hydraulic fracture intersections with cemented natural fractures. The impact of this finding is that we need fully 3D computer models capable of accounting for bypass and mixed mode I-III fracture propagation in order to realistically simulate subsurface hydraulic fracture geometries. / text

This dissertation presents three inter-related studies. Chapter 2 presents a study of scaling of crack propagation in rubber sheets. Two different scaling laws for supersonic and subsonic cracks were discovered. Experiments and numerical simulations have been conducted to investigate subsonic and supersonic cracks. The experiments are performed at 85 °C to suppress strain-induced crystallites that complicate experiments at lower temperature. Calibration experiments were performed to obtain the parameters needed to compare with a theory including viscous dissipation. Both experiments and numerical simulations support supersonic cracks, and a transition from subsonic to supersonic is discovered in the plot of experimental crack speed curves versus extension ratio for different sized samples. Both experiments and simulations show two different scaling regimes: the speed of subsonic cracks scales with the elastic energy density while the speed of supersonic cracks scales with the extension ratio. Crack openings have qualitatively different shapes in the two scaling regimes. Chapter 3 describes a theory of oscillating cracks. Oscillating cracks are not seen very widely, but observed in rubber and gels. A theory has been proposed for the onset of oscillation in gels, but the oscillation of cracks in rubber has not been explained. This study provides a theory able to describe both rubber and gels and recover the experimental phase diagram for oscillating cracks in rubber. The main new idea is that the oscillations of cracks follow from basic features of fracture mechanics and are independent of details of the crack equation of motion. From the fact that oscillations exist, one can deduce some conditions on forms that equations of motion can take. A discrete model of hydraulic fracture is mentioned in Chapter 4. Hydraulic fracturing is a stimulation treatment wherein fluids are injected into reservoirs under high pressure to generate fractures in reservoirs. In this study, a lattice-based pseduo-3D model is developed to simulate hydraulic fracturing. This mode has been validated via a comparison with the KGD model. A series of pilot simulations was systematically tested for complex geometries under more realistic operation conditions, including flexible boundary conditions, randomness in elastic properties of shales and perforations. The simulation results confirm that perforation is likely to increase the complexity of fracture networks; the results also suggest that the interference between neighboring fractures is key to fracture network formation. / text

Thesis (Ph.D.)--Memorial University of Newfoundland. / Typescript. Restricted until November 1997. Bibliography: leaves [125]-133. Also available online.