scholarly journals The competition between fracture nucleation, propagation, and coalescence in dry and water-saturated crystalline rock

Solid Earth ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 375-387
Author(s):  
Jessica A. McBeck ◽  
Wenlu Zhu ◽  
François Renard

Abstract. The continuum of behavior that emerges during fracture network development in crystalline rock may be categorized into three end-member modes: fracture nucleation, isolated fracture propagation, and fracture coalescence. These different modes of fracture growth produce fracture networks with distinctive geometric attributes, such as clustering and connectivity, that exert important controls on permeability and the extent of fluid–rock interactions. To track how these modes of fracture development vary in dominance throughout loading toward failure and thus how the geometric attributes of fracture networks may vary under these conditions, we perform in situ X-ray tomography triaxial compression experiments on low-porosity crystalline rock (monzonite) under upper-crustal stress conditions. To examine the influence of pore fluid on the varying dominance of the three modes of growth, we perform two experiments under nominally dry conditions and one under water-saturated conditions with 5 MPa of pore fluid pressure. We impose a confining pressure of 20–35 MPa and then increase the differential stress in steps until the rock fails macroscopically. After each stress step of 1–5 MPa we acquire a three-dimensional (3D) X-ray adsorption coefficient field from which we extract the 3D fracture network. We develop a novel method of tracking individual fractures between subsequent tomographic scans that identifies whether fractures grow from the coalescence and linkage of several fractures or from the propagation of a single fracture. Throughout loading in all of the experiments, the volume of preexisting fractures is larger than that of nucleating fractures, indicating that the growth of preexisting fractures dominates the nucleation of new fractures. Throughout loading until close to macroscopic failure in all of the experiments, the volume of coalescing fractures is smaller than the volume of propagating fractures, indicating that fracture propagation dominates coalescence. Immediately preceding failure, however, the volume of coalescing fractures is at least double the volume of propagating fractures in the experiments performed at nominally dry conditions. In the water-saturated sample, in contrast, although the volume of coalescing fractures increases during the stage preceding failure, the volume of propagating fractures remains dominant. The influence of stress corrosion cracking associated with hydration reactions at fracture tips and/or dilatant hardening may explain the observed difference in fracture development under dry and water-saturated conditions.

2020 ◽  
Author(s):  
Jessica A. McBeck ◽  
Wenlu Zhu ◽  
François Renard

Abstract. Different modes of fracture growth produce fracture networks with distinctive geometric attributes that exert important controls on the extent of fluid-rock interactions. We perform in situ X-ray tomography triaxial compression experiments on monzonite to investigate the influence of fracture nucleation, preexisting fracture propagation, and coalescence on fracture network development in crystalline rocks under crustal conditions. We impose a confining pressure of 20–35 MPa and then increase the differential stress in steps until the rock fails macroscopically. After each stress step we acquire a three-dimensional (3D) X-ray adsorption coefficient field from which we extract the 3D fracture network. To examine the influence of pore fluid on fracture network development, we perform two experiments under nominally-dry conditions and one under water-saturated conditions with 5 MPa pore fluid pressure. We develop a method of tracking individual fractures between subsequent tomographic scans that identifies whether fractures grow from the coalescence and linkage of several fractures or from the propagation of a single fracture. Throughout loading until shortly before failure in all of the experiments, the volume of coalescing fractures is smaller than the volume of propagating fractures, indicating that fracture propagation dominates coalescence. Immediately preceding failure, however, the volume of coalescing fractures is at least double the volume of propagating fractures in the experiments deformed at nominally dry conditions. In the water-saturated sample, although the volume of coalescing fractures increases during this stage, the volume of propagating fractures remains dominant. The influence of stress corrosion cracking associated with hydration reactions at fracture tips and/or dilatant hardening may explain the observed difference in fracture development under dry and water-saturated conditions. Our experimental data on fracture growth at different conditions provide new constraints in assessing fluid flow in subsurface fracture networks that are central to energy and environmental engineering practices.


Author(s):  
Hannes Hofmann ◽  
Tayfun Babadagli ◽  
Günter Zimmermann

The creation of large complex fracture networks by hydraulic fracturing is imperative for enhanced oil recovery from tight sand or shale reservoirs, tight gas extraction, and Hot-Dry-Rock (HDR) geothermal systems to improve the contact area to the rock matrix. Although conventional fracturing treatments may result in bi-wing fractures, there is evidence by microseismic mapping that fracture networks can develop in many unconventional reservoirs, especially when natural fracture systems are present and the differences between the principle stresses are low. However, not much insight is gained about fracture development as well as fluid and proppant transport in naturally fractured tight formations. In order to clarify the relationship between rock and treatment parameters, and resulting fracture properties, numerical simulations were performed using a commercial Discrete Fracture Network (DFN) simulator. A comprehensive sensitivity analysis is presented to identify typical fracture network patterns resulting from massive water fracturing treatments in different geological conditions. It is shown how the treatment parameters influence the fracture development and what type of fracture patterns may result from different treatment designs. The focus of this study is on complex fracture network development in different natural fracture systems. Additionally, the applicability of the DFN simulator for modeling shale gas stimulation and HDR stimulation is critically discussed. The approach stated above gives an insight into the relationships between rock properties (specifically matrix properties and characteristics of natural fracture systems) and the properties of developed fracture networks. Various simulated scenarios show typical conditions under which different complex fracture patterns can develop and prescribe efficient treatment designs to generate these fracture systems. Hydraulic stimulation is essential for the production of oil, gas, or heat from ultratight formations like shales and basement rocks (mainly granite). If natural fracture systems are present, the fracturing process becomes more complex to simulate. Our simulation results reveal valuable information about main parameters influencing fracture network properties, major factors leading to complex fracture network development, and differences between HDR and shale gas/oil shale stimulations.


Processes ◽  
2018 ◽  
Vol 6 (8) ◽  
pp. 113 ◽  
Author(s):  
Shen Wang ◽  
Huamin Li ◽  
Dongyin Li

To investigate the mechanism of hydraulic fracture propagation in coal seams with discontinuous natural fractures, an innovative finite element meshing scheme for modeling hydraulic fracturing was proposed. Hydraulic fracture propagation and interaction with discontinuous natural fracture networks in coal seams were modeled based on the cohesive element method. The hydraulic fracture network characteristics, the growth process of the secondary hydraulic fractures, the pore pressure distribution and the variation of bottomhole pressure were analyzed. The improved cohesive element method, which considers the leak-off and seepage behaviors of fracturing liquid, is capable of modeling hydraulic fracturing in naturally fractured formations. The results indicate that under high stress difference conditions, the hydraulic fracture network is spindle-shaped, and shows a multi-level branch structure. The ratio of secondary fracture total length to main fracture total length was 2.11~3.62, suggesting that the secondary fractures are an important part of the hydraulic fracture network in coal seams. In deep coal seams, the break pressure of discontinuous natural fractures mainly depends on the in-situ stress field and the direction of natural fractures. The mechanism of hydraulic fracture propagation in deep coal seams is significantly different from that in hard and tight rock layers.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-15 ◽  
Author(s):  
Yongxiang Zheng ◽  
Jianjun Liu ◽  
Yun Lei

The formation of the fracture network in shale hydraulic fracturing is the key to the successful development of shale gas. In order to analyze the mechanism of hydraulic fracturing fracture propagation in cemented fractured formations, a numerical simulation about fracture behavior in cemented joints was conducted based firstly on the block discrete element. And the critical pressure of three fracture propagation modes under the intersection of hydraulic fracturing fracture and closed natural fracture is derived, and the parameter analysis is carried out by univariate analysis and the response surface method (RSM). The results show that at a low intersecting angle, hydraulic fractures will turn and move forward at the same time, forming intersecting fractures. At medium angles, the cracks only turn. At high angles, the crack will expand directly forward without turning. In conclusion, low-angle intersecting fractures are more likely to form complex fracture networks, followed by medium-angle intersecting fractures, and high-angle intersecting fractures have more difficulty in forming fracture networks. The research results have important theoretical guiding significance for the hydraulic fracturing design.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Hannes Hofmann ◽  
Tayfun Babadagli ◽  
Günter Zimmermann

The creation of large complex fracture networks by hydraulic fracturing is imperative for enhanced oil recovery from tight sand or shale reservoirs, tight gas extraction, and hot-dry-rock (HDR) geothermal systems to improve the contact area to the rock matrix. Although conventional fracturing treatments may result in biwing fractures, there is evidence by microseismic mapping that fracture networks can develop in many unconventional reservoirs, especially when natural fracture systems are present and the differences between the principle stresses are low. However, not much insight is gained about fracture development as well as fluid and proppant transport in naturally fractured tight formations. In order to clarify the relationship between rock and treatment parameters, and resulting fracture properties, numerical simulations were performed using a commercial discrete fracture network (DFN) simulator. A comprehensive sensitivity analysis is presented to identify typical fracture network patterns resulting from massive water fracturing treatments in different geological conditions. It is shown how the treatment parameters influence the fracture development and what type of fracture patterns may result from different treatment designs. The focus of this study is on complex fracture network development in different natural fracture systems. Additionally, the applicability of the DFN simulator for modeling shale gas stimulation and HDR stimulation is critically discussed. The approach stated above gives an insight into the relationships between rock properties (specifically matrix properties and characteristics of natural fracture systems) and the properties of developed fracture networks. Various simulated scenarios show typical conditions under which different complex fracture patterns can develop and prescribe efficient treatment designs to generate these fracture systems. Hydraulic stimulation is essential for the production of oil, gas, or heat from ultratight formations like shales and basement rocks (mainly granite). If natural fracture systems are present, the fracturing process becomes more complex to simulate. Our simulations suggest that stress state, in situ fracture networks, and fluid type are the main parameters influencing hydraulic fracture network development. Major factors leading to more complex fracture networks are an extensive pre-existing natural fracture network, small fracture spacings, low differences between the principle stresses, well contained formations, high tensile strength, high Young’s modulus, low viscosity fracturing fluid, and large fluid volumes. The differences between 5 km deep granitic HDR and 2.5 km deep shale gas stimulations are the following: (1) the reservoir temperature in granites is higher, (2) the pressures and stresses in granites are higher, (3) surface treatment pressures in granites are higher, (4) the fluid leak-off in granites is less, and (5) the mechanical parameters tensile strength and Young’s modulus of granites are usually higher than those of shales.


2021 ◽  
Author(s):  
Rahul Prabhakaran ◽  
Giovanni Bertotti ◽  
Janos Urai ◽  
David Smeulders

Abstract. We investigate the spatial variation of 2D fracture networks digitized from the well-known Lilstock limestone pavements, Bristol Channel, UK. By treating fracture networks as spatial graphs, we utilize a novel approach combining graph similarity measures and hierarchical clustering to identify spatial clusters within fracture networks and quantify spatial variation. We use four graph similarity measures: fingerprint distance, D-measure, NetLSD, and portrait divergence to compare fracture graphs. The technique takes into account both topological relationship and geometry of the networks and is applied to three large fractured regions consisting of nearly 300,000 fractures spread over 14,200 sq.m. The results indicates presence of spatial clusters within fracture networks with that vary gradually over distances of tens of metres. One region is not influenced by faulting but still displays variation in background fracturing. Variation in fracture development in the other two regions are interpreted to be primarily influenced by proximity to faults that gradually gives way to background fracturing. Comparative analysis of the graph similarity-derived clusters with fracture persistence measures indicate that there is a general correspondence between patterns; however, additional variations are highlighted that is not obvious from fracture intensity and density plots. The proposed method provides a quantitative way to identify spatial variations in fracture networks which can be used to guide stochastic and geostatistical approaches to fracture network modelling.


Materials ◽  
2020 ◽  
Vol 13 (6) ◽  
pp. 1405 ◽  
Author(s):  
Petr Koudelka ◽  
Tomas Fila ◽  
Vaclav Rada ◽  
Petr Zlamal ◽  
Jan Sleichrt ◽  
...  

Several methods, including X-ray radiography, have been developed for the investigation of the characteristics of water-saturated quasi-brittle materials. Here, the water content is one of the most important factors influencing their strength and fracture properties, in particular, as regards to porous building materials. However, the research concentrated on the three-dimensional fracture propagation characteristics is still significantly limited due to the problems encountered with the instrumentation requirements and the size effect. In this paper, we study the influence of the water content in a natural quasi-brittle material on its mechanical characteristics and fracture development during in-situ four-point bending by employing high-resolution X-ray differential micro-tomography. The cylindrical samples with a chevron notch were loaded using an in-house designed four-point bending loading device with the vertical orientation of the sample. The in-house designed modular micro-CT scanner was used for the visualisation of the specimen’s behaviour during the loading experiments. Several tomographic scans were performed throughout the force-displacement diagrams of the samples. The reconstructed 3D images were processed using an in-house developed differential tomography and digital volume correlation algorithms. The apparent reduction in the ultimate strength was observed due to the moisture content. The crack growth process in the water-saturated specimens was identified to be different in comparison with the dry specimens.


2021 ◽  
Author(s):  
Zhao Hui ◽  
Sheng Guanglong ◽  
Huang Luoyi ◽  
Zhong Xun ◽  
Fu Jingang ◽  
...  

Abstract Accurately characterizing fracture network morphology is necessary for flow simulation and fracturing evaluation. The complex natural fractures and reservoir heterogeneity in unconventional reservoirs make the induced fracture network resulting from hydraulic fracturing more difficult to describe. Existing fracture propagation simulation and fracture network inversion methods cannot accurately match actual fracture network morphology. Considering the lightning breakdown similar as fracture propagation, a new efficient approach for inversion of fracture network morphology is proposed. Based on the dielectric breakdown model (DBM) for lightning breakdown simulation and similarity principle, an induced fracture propagation algorithm integrating reservoir in-situ stress, rock mechanical parameters, and stress shadow effect is proposed. The fractal index and random function are coupled to quantitatively characterize the probability distribution of induced fracture propagation path. At the same time, a matching rate function is proposed to quantitatively evaluate the fitting between fracture network morphology and the micro seismic data. Combined with automatic history matching method, the actual fracture network morphology can be inverted with the matching rate as objective function. The proposed approach is applied to fracture network simulation of mult-fractured horizontal wells of shale oil reservoir in China, and the fracture networks from inversion fit well with the micro seismic data. A simulation of 94 fractures in the 32 section of Well X2 shows that the well propagates more obvious branch fractures. The single-wing fracture network communicates approximately 200m horizontally and approximately 10m vertically. In single fracture flow simulation, it is necessary to consider the influence of complex fracture network morphology, but when simulating fluid flow for a single well or even a reservoir, only the main fracture needs to be considered. This paper proposes an induced fracture propagation algorithm that integrates reservoir in-situ stress, rock mechanical parameters, and stress shadowing effects. This algorithm greatly improves the calculation efficiency on the premise of ensuring the accuracy of induced fracture network morphology. The approach in this paper provides a theoretical basis for flow simulation of stimulated reservoirs and optimization of fracture networks.


Author(s):  
Stephanie G. Zihms ◽  
Helen Lewis ◽  
Tiago Siqueira de Miranda ◽  
Stephen A. Hall ◽  
James M. Somerville

Abstract: Comparing outcrop data to laboratory results is important to verify and validate experiments of analogue and reservoir materials especially regarding conditions for deformation experiments. This is important better understand highly complex carbonate reservoir strata and their response to changes in subsurface conditions, reducing subsurface uncertainty. This study develops methods to allow for a more straightforward comparison of outcrop data (m-scale) with experimentally created fracture arrays developed in cylindrical samples (cm-scale). The main objective is to assess usefulness of experimentally-produced fracture networks as analogues for subsurface structures, typically at the meter and above scale by developing new techniques to use the lab deformation. It analyses key characteristics of laboratory-induced fracture networks by adapting scanline methods to use with x-ray tomography (XRT) images to allow for comparison with outcrop and field data. To test and verify these new methods two low permeability carbonate samples were used for deformation testing and analysis. Applying the different scanline methods we show that they can be used to analyse lab induced fractures (mm to cm-scale) identified in XRT images for comparison with outcrop data (m-scale). In addition, these methods also allow for quantification of fracture network attributes e.g. fracture spacing, fracture apertures, orientation. This new data bridges the gap between micro-scanlines using thin sections and outcrop scanlines.


1981 ◽  
Vol 11 ◽  
Author(s):  
M.H. Bradbury ◽  
D. Lever ◽  
D. Kinsey

One of the options being considered for the disposal of radioactive waste is deep burial in crystalline rocks such as granite. It is generally recognised that in such rocks groundwater flows mainly through the fracture networks so that these will be the “highways” for the return of radionuclides to the biosphere. The main factors retarding the radionuclide transport have been considered to be the slow water movement in the fissures over the long distances involved together with sorption both in man-made barriers surrounding the waste, and onto rock surfaces and degradation products in the fissures.


Sign in / Sign up

Export Citation Format

Share Document