Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fracturing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Computational Modeling of Hydraulic Properties of a Sheared Single Rock Fracture

Flow through a single mechanically sheared Marcellus shale fracture was investigated computationally and experimentally. To provide a better understanding of the variation of hydraulic and geometrical characteristics of a fracture subjected to shearing, coupled shear flow tests on the fracture for four shearing displacement steps under constant normal stress were performed. At the end of each shearing step, computed tomography (CT) scans with resolution 26.8 μm were obtained and the corresponding fracture geometries were evaluated. The CT images were used to generate full aperture maps of the fracture configuration. In addition, average aperture maps were also created by averaging the full-resolution data over 10 × 10 pixels, smoothing out fine structural details. Computational modeling of water flow through the fractures at different shearing steps was performed using a modified local cubic law approach and the 3D full Navier–Stokes equations with the use of the ANSYS-Fluent software. Both the average aperture maps and full maps were used in these simulations. The experimental pressure drops of the fracture at shearing step 1, which has very small apertures, poorly matched the numerical results, quite likely because the fracture structure was inadequately captured by the scanning resolution. Shearing typically increased the aperture height of the fracture, whose features were then better captured by the CT scan. Good agreement between the experimental data and the numerical results of the full map for shearing step 2 was observed. The simulations were performed for both full and average aperture maps, and the effects of scan resolution and surface roughness on the accuracy of the results were studied. The modified local cubic law and full Navier–Stokes simulations of the averaged map fracture were found to be in good agreement. It was conjectured that this was because the nonlinear losses were insignificant for the smoothed out averaged map fracture. Similar comparisons with those of the full map showed agreement in trends, but there were some quantitative differences. The averaged fracture map simulations also predicted lower pressure drops compared to the full map, particularly for high flow rates. These differences were due to the fine-scale geometrical complexity (surface roughness) of fracture geometry that affects the fluid flow in the fracture. In conclusion, an improved cubic law model was also proposed, and its accuracy was verified by comparing its predictions with those of the Navier–Stokes simulations.

58 GEOSCIENCES↗
DOE OSTI2023

Using in-situ strain measurements to evaluate the accuracy of stress estimation procedures from fracture injection/shut-in tests

Fracture injection/shut-in tests are commonly used to measure the state of stress. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of Hayashi and Haimson. In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. In conclusion, the calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

58 GEOSCIENCES↗
DOE OSTI2018

Injection‐Induced Shear Slip and Permeability Enhancement in Granite Fractures

Abstract Shear stimulation, so‐called hydroshearing , is believed to reactivate critical or near‐critical preexisting fractures by pressurized fluid injection, causing them to slip and dilate for economic production from geothermal and unconventional petroleum reservoirs. However, little or no experiments have been conducted to directly probe the coupled hydromechanical responses and permeability evolution during fracture shearing. In this paper, we present the results of novel injection‐induced shear tests on cylindrical granite samples each containing a single tensile or saw cut fracture. The shear test on rough fractures demonstrates that retainable permeability enhancement can be achieved through dilatant shear slip. The characteristic stick‐slip behavior of the fracture during shearing is also observed. Fracture aperture controlled by effective normal stress dominates the fluid flow in stick state, while in slip state, the major contributor to production increase is the irreversible normal dilation caused by fracture shear slip. The experimental data show that depending on the fracture roughness, the fracture slip process includes two quasi‐static slip intervals with a slip velocity of ~10 −9 to ~10 −6 m/s and a dynamic slip interval with a slip velocity of ~10 −7 to ~10 −5 m/s. The fracture slip correlates well with the associated stress relaxation: a faster fracture slip induces a quicker stress relaxation, and vice versa. The influence of surface roughness on the hydromechanical responses during fracture shearing and the shear‐induced asperity degradation are also studied. The various experimental observations on rough and smooth fractures clearly indicate that the key component to shear stimulation is fracture self‐propping by asperities.

Ye, Zhi↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in theMidwest United States

This report presents final technical results for the project Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in the Midwest United States (DE-FE0023330). The project was a three-year effort consisting of seven technical tasks focused on defining geomechanical factors for CO 2 storage applications in deep saline rock formations in Ohio and the Midwest United States, because geomechancial issues have been identified as a significant risk factor for large-scale CO 2 storage applications. A basin-scale stress-strain analysis was completed to describe the geomechanical setting for rock formations of Ordovician-Cambrian age in Ohio and adjacent areas of the Midwest United States in relation to geologic CO 2 storage applications. The tectonic setting, stress orientation-magnitude, and geomechanical and petrophysical parameters for CO 2 storage zones and caprocks in the region were cataloged. Ten geophysical image logs were analyzed for natural fractures, borehole breakouts, and drilling-induced fractures. The logs indicated mostly less than 10 fractures per 100 vertical feet in the borehole, with mostly N65E principal stress orientation through the section. Geophysical image logs and other logs were obtained for three wells located near the sites where specific models were developed for geomechanical simulations: Arches site in Boone County, Kentucky; Northern Appalachian Basin site in Chautauqua County, New York; and E-Central Appalachian Basin site in Tuscarawas County, Ohio. For these three wells, 9,700 feet of image logs were processed and interpreted to provide a systematic review of the distribution within each well of natural fractures, wellbore breakouts, faults, and drilling induced fractures. There were many borehole breakouts and drilling-induced tensile fractures but few natural fractures. Concentrated fractures were present at the Rome-basal sandstone and basal sandstone-Precambrian contacts at the Arches and East-Central Appalachian Basin sites. Geophysical logs were utilized to develop local-scale geologic models by determining geomechanical and petrophysical parameters within the geologic formations. These data were ported to coupled fluid-flow and reservoir geomechanics multi-phase CO 2 injection simulations. The models were developed to emphasize the geomechanical layers within the CO 2 storage zones and caprocks. A series of simulations were completed for each site to evaluate whether commercial-scale CO 2 could be safely injected into each site, given site-specific geologic and geomechanical controls. This involved analyzing the simulation results for the integrity of the caprock, intermediate, and reservoir zones, as well quantifying the areal uplift at the surface. Simulation results were also examined to ensure that the stress-stress perturbations were isolated within the subsurface, and that there was only limited upward migration of the CO 2 . Simulations showed capacity to inject more than 10 million metric tons of CO 2 in a single well at the Arches and East Central Appalachian Basin sites without excessive geomechanical risks. Low-permeability rock layers at the Northern Appalachian Basin study area well resulted in very low CO 2 injection capacity. Fracture models developed for the sites suggests that the sites have sparse fracture network in the deeper Cambrian rocks. However, there were indicators in image logs of a moderate fracture matrix in the Rose Run Sandstone at the Northern Appalachian Basin site. Dual permeability fracture matrix simulations suggest the much higher injection rates may be feasible in the fractured interval. Guidance was developed for geomechanical site characterization in the areas of geophysical logging, rock core testing, well testing, and site monitoring. The guidance demonstrates that there is a suitable array of options for addressing geomechanical issues at CO 2 storage sites. Finally, a review of Marcellus and Utica-Point Pleasant shale gas wells and CO 2 storage intervals indicates that these items are vertically separated, except for the Oriskany sandstone and Marcellus wells in southwest Pennsylvania and northern West Virginia. Together, project results present a more realistic portrayal of geomechanical risk factors related to CO 2 storage for existing and future coal-fired power plants in Ohio.

58 GEOSCIENCES↗
DOE OSTI2017

Report on Developing a US-China Joint Project on CO2-Based Fracturing Techniques

There are several types of hydraulic fracturing methods that have been investigated in the past decades that use little or no water. The criteria for selecting an optimal hydraulic fracturing fluids are many, and no single fluid can meet all the desired characteristics. Developing a set of technologies by using different fracturing fluids can better meet individual geologic, economic, and environmental needs. CO2 as an alternative of water for hydraulic fracturing has its unique place among all other fracturing fluids, and can meet both economic and environmental objectives under certain circumstances. There are three main categories of CO2-based fracturing fluids, simple liquid CO2 (LCO2), directly thickened LCO2, and CO2-foams/emulsions. Although LCO2 has many advantages, critical disadvantages arise due to its high operating pressure, low viscosity, poor proppant-carrying capacity, and generation of shallower fracture networks. To date, techniques for direct thickening of LCO2 to overcome these limitations have proven ineffective. To increase the apparent viscosity of CO2, we propose to further develop CO2-foam (and emulsion) fracturing techniques. While CO2-foams have long been used in EOR, the required properties for hydraulic fracturing are different, and relatively few peer-reviewed publications and field test reports are available regarding their applications in fracturing, and the surfactants and additives used are often proprietary This report describes a brief state-of-the-art analysis of CO2-based hydraulic fracturing methods and then discusses a joint plan for a long-term R&D collaboration to advance the most promising CO2-based hydraulic fracturing techniques as part of the U.S.-China Clean Energy Research Center, Water Energy Technology (CERC-WET). The objective of this U.S.-China joint research project is to design and laboratory test a CO2-foam fracturing fluid system, and then conduct pilot field-tests on them in China’s Jilin Oilfield. This joint research project will benefit both U.S. and China for advancing their shared goal of developing better technologies for CO2-based fracturing. The Chinese CERC-WET has two teams joining the project, one from RIPED, responsible for aboveground infrastructure design and construction. Another team from the Jilin Oilfield will be responsible for the field test design and conducting the test. The U.S.CERC-WET team at LBNL will be responsible for providing the CO2-foam formulations and their rheology parameters to be tested in the field (a challenging task and a key to success); participating in the equipment design and field test design; assisting in analyzing field data; and reporting to U.S. DOE and CERC-WET.

58 GEOSCIENCES↗
NASA NTRS2017

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗
DOE OSTI2020

Alternative, Less Water-Intensive, Fracturing Fluids for Enhanced Geothermal Systems

Four fracturing fluids; water, CO2, CO2 with water, and CO2 with an aqueous solution of a CO2-reactive polymer, poly(allylamine) (PAA, 1wt%) were evaluated using a high-temperature true-triaxial fracturing apparatus and ½ foot side granite cubic samples. All three CO2-based fracturing fluids, CO2, CO2 with water, and CO2 with aqueous PAA fractured granite at higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When partially saturating the rock sample with PAA solution followed by fracturing with CO2, the volume expansion caused by CO2-induced cross-linking of PAA leads to a faster pressure increase due to the associated stress generated and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. It is concluded that CO2 (when injected in HDR) and CO2/PAA are the fluids of choice for stimulation of granitic rock samples under the studied geothermal P/T conditions. CO2/PAA, however, offer the following three additional advantages; 1) it requires significantly lower volumes of CO2, 2) fracture permeability is independent of injection flow rate, and 3) the reversible (previously reported) viscosity increase is beneficial to transport proppants if they are ever developed enhanced geothermal systems.

Jian, Guoqing↗
DOE OSTI2018

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.

58 GEOSCIENCES↗
DOE OSTI2015

Rapid imbibition of water in fractures within unsaturated sedimentary rock

The spontaneous imbibition of water and other liquids into gas-filled fractures in variably-saturated porous media is important in a variety of engineering and geological contexts. However, surprisingly few studies have investigated this phenomenon. In this paper, we present a theoretical framework for predicting the 1-dimensional movement of water into air-filled fractures within a porous medium based on early-time capillary dynamics and spreading over the rough surfaces of fracture faces. The theory permits estimation of sorptivity values for the matrix and fracture zone, as well as a dispersion parameter which quantifies the extent of spreading of the wetting front. Quantitative data on spontaneous imbibition of water in unsaturated Berea sandstone cores were acquired to evaluate the proposed model. The cores with different permeability classes ranging from 50 to 500 mD and were fractured using the Brazilian method. Spontaneous imbibition in the fractured cores was measured by dynamic neutron radiography at the Neutron Imaging Prototype Facility (beam line CG-1D, HFIR), Oak Ridge National Laboratory. Water uptake into both the matrix and the fracture zone exhibited square-root-of-time behavior. The matrix sorptivities ranged from 2.9 to 4.6 mm s -0.5 , and increased linearly as the permeability class increased. The sorptivities of the fracture zones ranged from 17.9 to 27.1 mm s -0.5 , and increased linearly with increasing fracture aperture width. The dispersion coefficients ranged from 23.7 to 66.7 mm 2 s -1 and increased linearly with increasing fracture aperture width and damage zone width. Both theory and observations indicate that fractures can significantly increase spontaneous imbibition in unsaturated sedimentary rock by capillary action and surface spreading on rough fracture faces. Fractures also increase the dispersion of the wetting front. In conclusion, further research is needed to investigate this phenomenon in other natural and engineered porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Numerical simulations of vertical growth of hydraulic fractures and brine migration in geological formations above the Marcellus shale

One of the critical environmental questions about hydraulic fracturing in shales is the potential for contamination of ground and surface water. There are two specific concerns arising from hydraulic treatments: 1) whether hydraulic fractures extend upward through overlying strata to reach overlying aquifers containing drinking water, and 2) whether injected fluids push native fluids upward into these overlying aquifers. In this work, the extent of likely fracture growth through overlying layers during hydraulic treatment of the Marcellus shale was estimated using a hydraulic fracture model. A wide range of material and fluid flow properties in a multi-layered geologic model was considered. The model was based on conditions and characteristics applicable to the Marcellus shale in that part of the Appalachian basin within southwestern Pennsylvania. Predictions of vertical termination frequencies for hydraulic fractures were used in a multi-layer model of the strata and natural fractures for studying brine migration through the natural and induced fracture network. NFFLOW, the software for explicitly modeling flow within networks of fractures, was utilized to compute transient flow rates according to the schedule of injected fluid during hydraulic fracturing. To aid our analysis, the modeled sequence of geologic strata was capped with a fictitious unfractured, but moderately-permeable layer, which serves as a monitoring zone. The analysis assumes one well lateral was placed in the middle of the Marcellus shale with hydraulic fractures penetrating layers in the model. The newly-developed geomechanical module within NFFLOW was used to represent stress-sensitivity of the fractures. This allows the opening and closing of fracture apertures with changes in fluid pressures within fracture segments. Pressure increases in the formations overlying the Tully limestone, indicating fluid flow, was observed due to the hydraulic stimulation; and the impact of these increased pressures on brine migration towards the surface was considered.

03 NATURAL GAS↗
DOE OSTI2016

The direct-current response of electrically conducting fractures excited by a grounded current source

Hydraulic fracture stimulation of low permeability reservoir rocks is an established and cross–cutting technology for enhancing hydrocarbon production in sedimentary formations and increasing heat exchange in crystalline geothermal systems. Whereas the primary measure of success is the ability to keep the newly generated fractures sufficiently open, long–term reservoir management requires a knowledge of the spatial extent, morphology, and distribution of the fractures — knowledge primarily informed by microseismic and ground deformation monitoring. To minimize the uncertainty associated with interpreting such data, we investigate through numerical simulation the usefulness of direct-current (DC) resistivity data for characterizing subsurface fractures with elevated electrical conductivity by considering a geophysical experiment consisting of a grounded current source deployed in a steel cased borehole. In doing so, the casing efficiently energizes the fractures with steady current. Finite element simulations of this experiment for a horizontal well intersecting a small set of vertical fractures indicate that the fractures manifest electrically in (at least) two ways: (1) a local perturbation in electric potential proximal to the fracture set, with limited farfield expression and (2) an overall reduction in the electric potential along the borehole casing due to enhanced current flow through the fractures into the surrounding formation. The change in casing potential results in a measurable effect that can be observed far from fractures themselves. Under these conditions, our results suggest that farfield, timelapse measurements of DC potentials can be interpreted by simple, linear inversion for a Coulomb charge distribution along the borehole path, including a local charge perturbation due to the fractures. As a result, this approach offers an inexpensive method for detecting and monitoring the time-evolution of electrically conducting fractures while ultimately providing an estimate of their effective conductivity — the latter providing an important measure independent of seismic methods on fracture shape, size, and hydraulic connectivity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗
DOE OSTI2015

Testing Novel CR-39 Detector Deployment System for Identification of Subsurface Fractures, Soda Springs, ID

The Environmental Protection Agency (EPA) has teamed with Battelle Energy Alliance, LLC (BEA) at Idaho National Laboratory (INL) to facilitate further testing of geologic-fracture-identification methodology at a field site near the Monsanto Superfund Site located in Soda Springs, Idaho. INL has the necessary testing and technological expertise to perform this work. Battelle Memorial Institute (BMI) has engaged INL to perform this work through a Work for Others (WFO) Agreement. This study continues a multi-year collaborative effort between INL and EPA to test the efficacy of using field deployed Cr-39 radon in soil portals. This research enables identification of active fractures capable of transporting contaminants at sites where fractures are suspected pathways into the subsurface. Current state of the art methods for mapping fracture networks are exceedingly expensive and notoriously inaccurate. The proposed WFO will evaluate the applicability of using cheap, readily available, passive radon detectors to identify conductive geologic structures (i.e. fractures, and fracture networks) in the subsurface that control the transport of contaminants at fracture-dominated sites. The proposed WFO utilizes proven off-the-shelf technology in the form of CR-39 radon detectors, which have been widely deployed to detect radon levels in homes and businesses. In an existing collaborative EPA/INL study outside of this workscope, CR-39 detectors are being utilized to determine the location of active transport fractures in a fractured granitic upland adjacent to a landfill site at the Fort Devens, MA that EPA-designated as National Priorities List (NPL) site. The innovative concept of using an easily deployed port that allows the CR-39 to measure the Rn-222 in the soil or alluvium above the fractured rock, while restricting atmospheric Rn-222 and soil sourced Ra from contaminating the detector is unique to INL and EPA approach previously developed. By deploying a series of these inexpensive detector-casing combinations statistical samples of the Rn-222 flux can be measured, elucidating the most communicative fractures (i.e. fractures that are actively transporting water and gasses). The Rn-222 measurements can then be used as an input to create a more accurate conceptual model to be used for transport modeling and related cleanup activities. If the team’s approach is demonstrated to be applicable to a wide variety of rock types and soil conditions it might potentially offer significant cost saving without a reduction in data quality at Monsanto Superfund and other sites underlain by fracture-dominated bedrock.

58 GEOSCIENCES↗
DOE OSTI2018

Fracture Diagnostics Using Low Frequency Electromagnetic Induction and Electrically Conductive Proppants

The monitoring and diagnostics of induced fractures are important for the real-time performance evaluation of hydraulic fracturing operations. Previous electromagnetic (EM) based studies show that single backbone tri-axial induction logging tools are promising candidates for real-time monitoring and diagnosis of fractures in open-hole wells. To support the development of field deployable tools, the concept must be tested in experiments, in a controllable environment, before it is tested under field-like conditions. To this end, we have developed numerical tools which can simulate any wellbore environment while logging hydraulic fractures with the induction tool. We have designed and built a prototype induction tool and performed two sets of tests to compare with numerical simulation results. The computational and experimental setup consists of tri-axial transmitter- and receiver-coils in co-axial, co-planar and cross-polarized configurations. Both lab and shallow earth measurements are shown to be in good agreement with simulations for all examined cases. The average relative and maximum discrepancies of the measured signals from the simulated ones were lower than 3% and 10%, respectively. With the prototype tool, strong signals sensitive to the fracture’s surface area and dip-angle were measured in co-axial coil configuration, while weaker signals sensitive to the fracture’s aspect ratio were observed in the co-planar configuration. Cross-polarized signals are also shown to be strong and sensitive to the fracture’s dip. Lastly, we resolved the detectable components of the signal tensor to obtain parameters of simplified fracture geometries. The inversion algorithm, a derivative free directional search model, uses an objective function defined as a combination of co-axial and cross-polarized signals from different tool spacings to provide a well behaved global minimum. The robustness of inversion algorithm is tested on synthetic data for single cluster fractures in homogeneous and heterogeneous background electrical conductivity. All the effective model parameters for different cases, fracture conductivity (width), propped area (length) and dip-angle, are shown to be recovered with sufficient accuracy. We also evaluated the effect of multiple neighboring fractures and suggested a multi-cluster inversion algorithm that can accurately recover the proppant distribution in each fracture. Based on these results we provide a downhole tool design, with detailed specifications, that can effectively recover far-field proppant distribution in hydraulic fractures.

02 PETROLEUM↗
DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗
DOE OSTI2021

Role of Mineralogy in Controlling Fracture Formation

The presence of fractures in caprocks can pose increased risks in subsurface energy systems and processes like CO 2 sequestration by introducing high-permeability leakage paths. Fracture apertures and permeability can be altered through mineral dissolution and precipitation reactions, but the reactive evolution of fractures is not well understood. In fractures, minerals that are otherwise inaccessible to reactive fluids can become exposed, resulting in mineral reactions unpredicted by bulk formation data. This work seeks to understand the relationship between mineralogy and fracture formation to enhance our understanding of reactive fracture evolution and CO 2 leakage potential. Here, the mineral compositions of mechanically induced fracture surfaces in samples of the Mancos and Marcellus shales have been quantified and compared to those of the near-fracture matrices using imaging and bulk X-ray diffraction (XRD) data. In the Mancos shale, the concentrations of clay minerals are enhanced along fracture surfaces with respect to the bulk, and the fracture is most likely to form at kaolinite–kaolinite interfaces. Further, evaluation of the mineralogical spatial variability through cross-correlation analysis of the surrounding matrix in images of samples cut perpendicular to the fracture shows that clay is 16.7 times more likely to be present than carbonate minerals near the fracture surface. The high correlation persists roughly 200 μm into the surrounding matrix for the Mancos sample and implies that the fracture formed within a defined clay-rich lithofacies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hydraulic Modeling of Induced and Propagated Fractures: Analysis of Flow and Pressure Data From Hydromechanical Experiments in the COSC-1 Deep Borehole in Crystalline Rock Near Åre, Sweden

To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.

58 GEOSCIENCES↗