Hydraulic fracturing is an indispensable technique for developing unconventional resources such as shale gas and tight oil. When hydraulic fractures interact with pre-existing natural fractures, it can result in a complex fracture network. The interaction depends on in-situ stresses, rock and natural fracture mechanical properties, approach angle and hydraulic fracture treatment parameters. Most simulation studies treat natural fractures as frictional interfaces with cohesive properties. However, from core observation, partially cemented and fully cemented natural fractures are widely present and it is not clear whether they would fit the common description. In this study, semi-circular bending test is utilized to examine the propagation paths and strength of samples with pre-existing cemented fractures. Synthetic hydrostone samples are used to represent the rock and different inclusion slices with different mechanical properties are used to mimic cemented natural fractures. In a series of experiments, we assess the influence of the fracture approach angle, inclusion strength, and inclusion thickness on fracture propagation. Current results show that fractures tend to cross the inclusion when the approach angle is high and divert into the inclusion when the approach angle is low. The crossing surface is not a clean cut, but often has a jog distance. The thickness of the inclusion does not change the crossing/diverting behavior for orthogonal approaching samples, however it does change the jog distance along the interface. Preliminary simulation results using finite element software, ABAQUS, are presented better to analyze the experimental observations. The assessments of fracture interaction in this study are in good agreement with previous work and theories. / text

Natural hydraulic fractures (NHFs) are fractures whose growths are driven by fluid loading. The fluid flow properties of the host rock have a primary, but hitherto little appreciated control on the NHF propagation rates. This study focuses on investigating the impacts of host rock fluid flow on the propagation and pattern development of multiple NHF in a poroelastic media. A realistic geomechanical model is developed to combine both the fluid flow and mechanical interactions between multiple fractures. The natural hydraulic fracture propagation is observed to consist of a series of crack-seal processes indicating incremental stop-start growth. Growth timing is on the scale of millions of years based on recent natural fracture growth reconstructions. These time scales are compatible with some model scenarios. My newly developed numerical model captures the crack-seal process for multiple NHF propagation. A sensitivity study conducted to investigate the impacts of different fluid flow properties on NHF propagation shows that permeability is a predominate influence on the timescale of NHF development. In low-permeability rocks, fractures have more stable initiation and much longer propagation timing compared to those in high-permeability rocks. Another aspect of great interest is the influence of fluid flow on fracture spacing and pattern development for multiple NHFs propagation in a poroelastic environment. My new poroelastic geomechanial model combines the natural hydraulic fracturing mechanism with the mechanical interactions between fractures. The numerical results show that as host rock permeability decreases, more fractures can propagate and a much smaller spacing is reached for a given fracture set. The low permeability slows down the propagation of long fractures and prevents them from dominating the fracture pattern. As a result, more fractures are able to grow at a similar speed and a more closely spaced fracture pattern is achieved for either regularly spaced or randomly distributed multiple fractures in low-permeability rocks. Investigation is also conducted in analyzing the distributions of fracture attributes (length, aperture and spacing) in low- and high-permeability rocks. For shales with high subcritical index, low permeability helps the fractures propagate more closely spaced instead of clustering. Meanwhile, in low-permeability rocks, factures have relatively smaller apertures, which lead to a slower fracture opening rate. The competition between the slow fracture opening rate and quartz precipitation rate will affect the effective permeability and porosity of the naturally fractured reservoir. However, the competition is trivial in high-permeability rocks. Other factors, such as reservoir boundary condition, layer thickness, subcritical index and pattern development stage, all have considerable impact on fracture pattern development and attribute distribution in a poroelastic media. / text

Modeling the stimulated reservoir volume during hydraulic fracturing is important to geothermal and petroleum reservoir stimulation. The interaction between a hydraulic fracture and pre-existing natural fractures exerts significant control on stimulated volume and fracture network complexity. This thesis presents a boundary element and finite difference based method for modeling this interaction during hydraulic fracturing process. In addition, an improved boundary element model is developed to more accurately calculate the total stimulated reservoir volume. The improved boundary element model incorporates a patch to calculate the tangential stresses on fracture walls accurately, and includes a special crack tip element at the fracture end to capture the correct stress singularity the tips The fracture propagation model couples fluid flow to fracture deformation, and accounts for fracture propagation including the transition of a mechanically-closed natural fractures to a hydraulic fracture. The numerical model is used to analyze a number of stimulation scenarios and to study the resulting hydraulic fracture trajectory, fracture aperture, and pressures as a function of injection time. The injection pressure, fracture aperture profiles shows the complexity of the propagation process and its impact on stimulation design and proppant placement. The injection pressure is observed to decrease initially as hydraulic fracture propagates and then it either increases or decreases depending on the factors such as distance between hydraulic fracture and natural fracture, viscosity of the injected fluid, injection rate and also other factor that are discussed in detail in below sections. Also, the influence of flaws on natural fracture in its opening is modeled. Results shows flaws that are very small in length will not propagate but are influencing the opening of natural fracture. If the flaw is located near to one end tip the other end tip will likely propagate first and vice versa. This behavior is observed due to the stress shadowing effect of flaw on the natural fracture. In addition, sequential and simultaneous injection and propagation of multiple fractures is modeled. Results show that for sequential injection, the pressure needed to initiate the later fractures increases but the geometry of the fractures is less complicated than that obtained from simultaneous injection under the same fracture spacing and injection. It is also observed that when mechanical interaction is present, the fractures in sequential fracturing have a higher width reduction as the later fractures are formed

Sets of natural opening-mode fractures in sedimentary rocks may show a variety of types of aperture-size distributions. A frequently documented size distribution type, in the literature and in data presented here, is the power law. The emergence of power-law distributions of fracture aperture and length sizes has been simulated using various quasi-mechanical fracture-growth routines but models based on linear-elastic fracture mechanics rarely produce such patterns. I collected a fracture-size dataset of unprecedented size and resolution using core and field methods and scanning electron microscope-based cathodoluminescence (SEM-CL) images. This dataset confirms the prevalence of power laws with a narrow range of power-law exponents among fractures that contain synkinematic cement. Organized microfractures are ubiquitous in sandstones. A fracture-growth simulation I devised reproduces observed size-scaling patterns by distributing fracture-opening increments among actively growing fractures. The simulated opening increments have a uniform size, which can be specified; uniform opening size is consistent with observations of narrow ranges of micron-scale widths of opening increments within crack-seal texture in natural fractures. Thus power-law size scaling of natural fractures can be explained using non-power-law (uniform-sized) opening increments, arranged using rules designed to simulate the effects of cement precipitation during fracture opening. A fundamental shortcoming of previous models of fracture-set evolution is the absence of a test because only natural fracture end states, not growth histories, could be measured. Using a technique to constrain fracture timing based on fluid inclusion microthermometry and thermal history modeling, I tested growth models by reconstructing the opening history of a set of natural fractures in the Triassic El Alamar Formation in northeast Mexico. The natural-fracture data show that, consistent with simulations, new microscopic fractures are continually introduced during natural fracture pattern evolution. As well, larger fractures represent sites of concentrated reactivation, although smaller fractures may be reactivated after long periods of quiescence. The pattern likely arises through feedback between fracture growth and the mechanically adhesive effects of contemporaneous fracture cement deposition. The narrow range in power-law exponents documented among fractures can help improve estimates of meter-scale large-fracture spacing where limited fracture samples are available. / text

It is essential to understand natural fracture systems embedded in shale-gas reservoirs and the stress fields that influence how induced fractures form in targeted shale units. Multicomponent seismic technology and elastic seismic stratigraphy allow geologic formations to be better images through analysis of different S-wave modes as well as the P-wave mode. Significant amounts of energy produced by P-wave sources radiate through the Earth as downgoing SV-wave energy. A vertical-force source is an effective source for direct SV radiation and provides a pure shear-wave mode (SV-SV) that should reveal crucial information about geologic surfaces located in anisotropic media. SV-SV shear wave modes should carry important information about petrophysical characteristics of hydrocarbon systems that cannot be obtained using other elastic-wave modes. Regardless of the difficulties of extracting good-quality SV-SV signal, direct shear waves as well as direct P and converted S energy should be accounted for in 3C seismic studies. Acquisition of full-azimuth seismic data and sampling data at small intervals over long offsets are required for detailed anisotropy analysis. If 3C3D data can be acquired with improved signal-to-noise ratio, more uniform illumination of targets, increased lateral resolution, more accurate amplitude attributes, and better multiple attenuation, such data will have strong interest by the industry. The objectives of this research are: (1) determine the feasibility of extracting direct SV-SV common-mid-point sections from 3-C seismic surveys, (2) improve the exploration for stratigraphic traps by developing systematic relationship between petrophysical properties and combinations of P and S wave modes, (3) create compelling examples illustrating how hydrocarbon-bearing reservoirs in low-permeable rocks (particularly anisotropic shale formations) can be better characterized using different S-wave modes (P-SV, SV-SV) in addition to the conventional P-P modes, and (4) analyze P and S radiation patterns produced by a variety of seismic sources. The research done in this study has contributed to understanding the physics involved in direct-S radiation from vertical-force source stations. A U.S. Patent issued to the Board of Regents of the University of Texas System now protects the intellectual property the Exploration Geophysics Laboratory has developed related to S-wave generation by vertical-force sources. The University’s Office of Technology Commercialization is actively engaged in commercializing this new S-wave reflection seismic technology on behalf of the Board of Regents. / text

Drilling is crucial to many industries, including hydrocarbon extraction, CO2 sequestration, geothermal energy, and others. During penetrating the subsurface rocks, drilling fluid (mud) is used for drilling bit cooling, lubrication, removing rock cuttings, and providing wellbore mechanical stability. Significant mud loss from the wellbore into the surrounding formation causes fluid lost-circulation incidents. This phenomenon leads to cost overrun, environmental pollution, delays project time and causes safety issues. Although lost-circulation exacerbates wellbore conditions, prediction of the characteristics of subsurface formations can be obtained. Generally, four formation types cause lost-circulation: natural fractures, and induced fractures, vugs and caves, and porous/permeable medium. The focus in this work is on naturally fractured formations, which is the most common cause of lost circulation. In this work, a novel prediction tool is developed based on analytical solutions and type-curves (TC). Type-curves are derived from the Cauchy equation of motion and mass conservation for non-Newtonian fluid model, corresponding to Herschel-Bulkley model (HB). Experimental setup from literature mimicking a deformed fracture supports the establishment of the tool. Upscaling the model of a natural fracture at subsurface conditions is implemented into the equations to achieve a group of mud type-curves (MTC) alongside another set of derivative-based mud type-curves (DMTC). The developed approach is verified with numerical simulations. Further, verification is performed with other analytical solutions. This proposed tool serves various functionalities; It predicts the volume loss as a function of time, based on wellbore operating conditions. The time-dependent fluid loss penetration from the wellbore into the surrounding formation can be computed. Additionally, the hydraulic aperture of the fracture in the surrounding formation can be estimated. Due to the non-Newtonian behavior of the drilling mud, the tool can be used to assess the fluid loss stopping time. Validation of the tool is performed by using actual field datasets and published experimental measurements. Machine-Learning is finally investigated as a complementary approach to determine the flow behavior of mud loss and the corresponding fracture properties.

Understanding the mechanisms that govern hydraulic fracturing applications in unconventional formations, such as gas-bearing shales, is of increasing interest to the petroleum upstream industry. Among such mechanisms, the geomechanical interactions between hydraulic fractures and pre-existing fractures on one hand, and simultaneous multiple hydraulic fractures on the other hand are seen of high importance. Although the petroleum engineering and related literature contains a number of studies that discusses such topics of hydraulic fracture interactions, there still remain some aspects that require answers, validations, or further supporting data. Particularly, experimental evidence is fairly scarce and keenly needed to solidify the understanding of such complex applications. In this work, the investigation methodology uses a series of hydraulic fracturing laboratory tests performed on synthetic rocks made of gypsum-based cements such as hydrostone and plaster in various experimental set ups. Those laboratory tests aim to closely investigate hydraulic fracture intersection with pre-existing fractures by assessing some factors that govern its outcomes. Specifically, the roles of the pre-existing fracture cementation, aperture, and relative height on the intersection mode are examined. The results show dominant effect of the cement-fill type relative to the host-rock matrix in determining whether hydraulic fracture crossing the pre-existing interface may occur. Similarly, hydraulic fracture height relative to the height of the pre-existing fracture may dictate the intersection results. However, the intersection mode seems to be insensitive of the pre-existing fracture aperture. Moreover, simultaneous multi-fracture propagation is examined and found to be impacted by the interference of the stresses induced from each fracturing source on neighboring fracturing sources. Such stress interference increases as the number of the propagating hydraulic fractures increase. While hydraulic fractures initiating from fracturing sources located in the middle of the fracturing stage seem to have inhibited propagation, outer hydraulic fractures may continue propagating with outward curvatures. Overall, the experimental results and analyses offer more insights for understanding hydraulic fracture complexity in unconventional formations. / text

La fracturation hydraulique est au cœur d'un certain nombre de phénomènes naturels et induits et est cruciale pour un développement durable de la production de ressources énergétiques. Compte tenu de son rôle crucial, ce phénomène a été pris en compte au cours des trois dernières décennies par le monde académique. Néanmoins, un certain nombre d'aspects très importants de ce processus ont été systématiquement négligés par la communauté. Deux des plus remarquables sont l'incapacité de la grande majorité des modèles existants à aborder la propagation des fractures hydrauliques dans les massifs rocheux fracturés où l'injection de fluide peut à la fois conduire à la fracturation de la roche intacte et à la réactivation de fractures préexistantes. Un autre aspect essentiel de ce processus est qu'il est intrinsèquement tridimensionnel, ce qui est souvent négligé par les modèles actuellement disponibles. Pour aborder ce problème essentiel, un modèle hydro-mécanique couplé basé sur la méthode des éléments discrets a été développé. La masse rocheuse est ici représentée par un ensemble d'éléments discrets interagissant à travers des lois de contact cohésifs qui peuvent se casser pour former des fissures à l'intérieur du milieu simulé. Ces fissures peuvent se coalescer pour former des fractures. Une méthode de volume fini est utilisée pour simuler l'écoulement de fluide entre les éléments discrets. L'écoulement est calculé en fonction de la déformation de l'espace poreux dans le milieu intact et de l'ouverture des fissures dans les fractures. De plus, les fractures naturelles sont modélisées explicitement de sorte qu'elles peuvent présentées des comportements mécanique et hydraulique différents de ceux de la matrice rocheuse intacte. La simulation des processus de fracturation hydraulique dans un milieu initialement intact en considérant plusieurs points d'injection plus ou moins espacés a permis de mettre en évidence l'évolution spatio-temporelle des fractures hydrauliques et de quantifier l'impact des différentes stratégies d'injection sur des indices représentatifs du volume fracturé, de l'intensité et de la densité des fractures ou encore sur la pression de fluide au niveau du puits. De plus, l'injection dans une fente de perforation non alignée sur le plan de contrainte minimum a génère des fractures hydrauliques non planaires percolantes si la connectivité est faible, ce qui peut être gênant pour la mise en place du proppant. En outre, des interactions fortes prennent place entre des fractures hydrauliques étroitement espacées ont été mises en évidence grâce au le suivi de la orientation de contrainte principale locale et ont révélé l'importance des effets d'ombre de contrainte. Des solutions sont proposées pour optimiser les traitements multiples à partir d'un puits de forage non parfaitement aligné. Enfin, l'interaction entre une seule fracture hydraulique et une seule fracture naturelle de propriétés et d'orientations variables a été étudiée à l'aide du modèle proposé. L'évolution de la fracture hydraulique et la réponse globale de l'échantillon ont été enregistrées d'une manière comparable aux données expérimentales existantes pour établir un pont entre les résultats expérimentaux et numériques. Les fractures naturelles persistantes semblent être des barrières pour la fracture hydraulique si leur conductance est élevée par apport a celle de la matrice ou si leur raideur est faible par rapport a la rigidité du milieu environnant. D'autre part, une faible rigidité dans les discontinuités non persistantes pourrait provoquer une bifurcation de la fracture hydraulique principale. De plus, des angles d'approche élevés et des contraintes différentielles fortes semblent favoriser le croisement de la fracture naturelle alors que des angles faibles engendrent plutôt un glissement ou une dilatation par cisaillement de la partie du plan qui n'est pas affectée par la perturbation de la contrainte. / Hydraulic fracturing is at the core of a number of naturally occurring and induced phenomena and crucial for a sustainable development of energy resource production. Given its crucial role this process has been given increasing attention in the last three decades from the academic world. Nonetheless a number of very significant aspects of this process have been systematically overlooked by the community. Two of the most notable ones are the inability of the vast majority of existing models to tackle at once the propagation of hydraulic fractures in realistic, fractured rocks-masses where hydraulic fracturing is a competing dipole mechanism between fracturing of the intact rock and re-activation of exiting fracture networks. Another essential aspect of this process is that it is intrinsically three-dimensional which is neglected by most models. To tackle this vital problem taking into account these pivotal aspects, a fully coupled hydro-mechanical model based on the discrete element method has been developed. The rock mass is here represented by a set of discrete elements interacting through elastic-brittle bonds that can break to form cracks inside the simulated medium. Theses cracks can coalesce to form fractures. A finite volume scheme is used to simulate the fluid flow in between these discrete elements. The flow is computed as a function of the pore space deformation in the intact medium and of the cracks' aperture in the fractures. Furthermore, the natural fractures are modelled explicitly and present mechanical and hydraulic properties different from the rock matrix. Employing this model in an intact numerical specimen, single fluid injection and multiple closely spaced sequential injections, enabled the description the full spatio-temporal evolution of HF propagation and its impact on quantitative indexes used in description of hydraulic fracturing treatments, such as fractured volume, fracture intensity and down-the-hole pressure for different control parameters and in-situ stress-fields. Moreover, injections from perforation slots which are not well aligned to the minimum stress plane showed possible creation of percolating non-planar hydraulic fractures of low connectivity, which can be troublesome for proppant placement. Also, strong interactions between closely spaced HF were highlighted by tracking the local principal stress rotation around the injection zones, emphasizing the importance of stress shadow effects. Optimization solutions are proposed for multiple treatments from a non-perfectly aligned wellbore. Finally, interaction between a single hydraulic fracture and a single natural fracture of varying properties and orientations was studied using the proposed model. The evolution of the hydraulic fracture and the global response of the specimen were recorded in a way comparable to existing experimental data to bridge the experimental and numerical findings. Persistent natural fractures appeared to be barriers for the hydraulic fracture if their conductance is high compared to the matrix conductivity or if their stiffness is significantly low compared to the rock matrix rigidity. Low stiffness in non-persistent defects might also cause a bifurcation of the main hydraulic fracture due to the local stress field perturbation around the defect and ahead of the hydraulic fracture tip. Furthermore, high approach angles and differential stresses seemed to favour crossing of the natural fracture while low angles enable shear slippage or dilation on the part of the plane which is not affected by the local stress perturbation.