Search NASA⌕ Search

SEARCH · Search NASA

Results for “well fractures”

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 181 records · Page 10

Stress-sensitivity of fracture conductivity of Tuscaloosa Marine Shale cores

Production decline from wells producing from shale rocks is typically quite rapid. One possible factor responsible for such decline is loss in fracture conductivity. This is true in the case of the Tuscaloosa Marine shale reservoir. In this work, the stress-dependent fracture conductivity of TMS core samples is investigated using five propped cores under varying stress conditions. In addition, the effects of static rock mechanical properties and rock mineralogy on conductivity decline are investigated. The mineralogy of each core sample was determined by X-ray Diffraction analysis. Fracture conductivity was measured using a Hassler-type pressure core holder with a pressure rating of 15,000 psi (103.42 MPa). Estimated rock mechanical properties showed that shale anisotropy could affect fracture conductivity. The time decline in conductivity revealed two possible decline trends, a) decline from the reduction in induced fracture width, and b) decline from the healing of developed micro-cracks. An exponential decline in fracture conductivity with increasing confinement pressure was observed with an average decline rate constant of 3.15 x 10 -10 mPa -1 . While rock mineralogy did not have any clear relationship with conductivity decline, the effects of rock mechanical properties were only important at 10.34 MPa confinement pressure and above. The observation showed that a high Young's modulus correlate to high conductivity, whereas a high Poisson's ratio correlate to low conductivity. This study provides useful observations for optimizing fracture design to enhance well productivity in the TMS.

02 PETROLEUM↗

Changes in environmental and engineered conditions alter the plasma membrane lipidome of fractured shale bacteria

ABSTRACT Microorganisms that persist in fractured shale reservoirs cause several problems including secreting foul gases and forming biofilms. Current biocontrol measures often fail due to limited knowledge of their in situ activities. The plasma membrane protects the cell, mediates many of its critical functions, and responds to intracellular cues and ecological perturbations through physicochemical modifications. As such, it provides valuable insight into the physiological adaptation of microorganisms in disturbed environmental systems. Here, we (i) demonstrate how changes in salinity and hydraulic retention time (HRT) influence the plasma membrane intact polar lipid (IPL) chemistry of model bacterium, Halanaerobium congolense WG10, and mixed microbial consortia enriched from shale-produced fluids and (ii) elucidate adjustments in membrane IPL chemistry during biofilm growth relative to planktonic cells. We incubated H. congolense WG10 in chemostats under three salinities (7%, 13%, and 20% NaCl), operated under three HRTs (19.2, 24, and 48 h), and in drip flow biofilm reactors under the same salinity gradients. Also, mixed microbial consortia in produced fluids were enriched in triplicate chemostat vessels under three HRTs (19.2, 24, and 72 h) and biofilm reactors. Lipids were analyzed by ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Our results show that phosphatidylglycerols, cardiolipins, and phosphatidylethanolamines were predominantly enriched in planktonic H. congolense WG10 cells grown at hypersalinity (20%) compared to optimum (13%). In addition, several zwitterionic phosphatidylcholines and phosphatidylethanolamines were higher in abundance during biofilm growth. These observations suggest that microbial adaptation and biofilm formation in fractured shale are enabled by strategic plasma membrane IPL chemistry adjustments. IMPORTANCE Microorganisms inadvertently introduced into the shale reservoir during fracturing face multiple stressors including brine-level salinities and starvation. However, some anaerobic halotolerant bacteria adapt and persist for long periods of time. They produce hydrogen sulfide, which sours the reservoir and corrodes engineering infrastructure. In addition, they form biofilms on rock matrices, which decrease shale permeability and clog fracture networks. These reduce well productivity and increase extraction costs. Under stress, microbes remodel their plasma membrane to optimize its roles in protection and mediating cellular processes such as signaling, transport, and energy metabolism. Hence, by observing changes in the membrane lipidome of model shale bacteria, Halanaerobium congolense WG10, and mixed consortia enriched from produced fluids under varying subsurface conditions and growth modes, we provide insight that advances our knowledge of the fractured shale biosystem. We also offer data-driven recommendations for improving biocontrol efficacy and the efficiency of energy recovery from unconventional formations.

03 NATURAL GAS↗

Fallon, NV FORGE well 21-31 Lithology, Mineral, Image Log, and Injection Test data

Attached is 5 datasets collected from Fallon FORGE well 21-31 located in Churchill County, Nevada collected and interpreted between Feburary 2018-January 2020. This submission includes 1) new lithology interpretation derived from petrographic thin sections of sidewall cores taken from well 21-31, 2) X-Ray Diffraction (XRD) interpretation of sidewall cores and cuttings from well 21-31, 3) Hyperspectral (Short-Wave Infrared and Long-Wave Infrared) interpretation of cuttings and sidewall cores from well 21-31, 4) Interpretation of Formation MicroImager (FMI) and Borehole Televiewer (BHTV) image logs for natural fractures and induced structures for well 21-31, 5) Injection test data and associated analytical model parameters from well testing of well 21-31.

15 GEOTHERMAL ENERGY↗

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗

Optimal pumping schedule with high-viscosity gel for uniform distribution of proppant in unconventional reservoirs

Hydrocarbon recovery from a hydraulic-fractured unconventional reservoir can be enhanced by replacing slickwater with non-Newtonian high-viscosity gel as a fracturing fluid. However, in the literature, the role of gel rheology on hydrocarbon production from a fractured reservoir is not very well understood. This lack of understanding makes it difficult for the industry to choose appropriate viscosity parameters of fracturing gels. Motivated by these limitations, first, we have developed a high-fidelity model of non-Newtonian fluid flow to simulate fracture propagation, fluid leak-off, and gel flowback processes to predict the oil production from a hydraulic fractured reservoir. Simulation results of the high-fidelity model show how fracture geometry, formation damage, and gel cleanup process varies with the fracturing fluid rheology. Therefore, the developed model is used to understand the effect of gel rheology on hydrocarbon production via a sensitivity analysis, which provides a set of rheological parameters that maximize the cumulative amount of oil production. Finally, utilizing this optimal high-viscosity gel, a model predictive controller is designed to obtain an optimal pumping schedule necessary for producing a target fracture geometry and proppant concentration. The closed-loop simulation results demonstrate that the obtained pumping schedule successfully achieved the desired fracture geometry and proppant concentration at the end of pumping, leading to the maximum oil production from an unconventional reservoir.

13 HYDRO ENERGY↗

Utah FORGE: Discrete Fracture Network (DFN) Data

The FORGE team is making these fracture models available to researchers wanting a set of natural fractures in the FORGE reservoir for use in their own modeling work. They have been used to predict stimulation distances during hydraulic stimulation at the open toe section of well 16A(78)-32. These fracture sets are fully stochastic and do not contain the deterministic set that matches the pilot well 58-32 FMI data. Well 58-32 has been completed and 16A(78)-32 is to be drilled as part of Phase 3. The original .fab files are not included due to redundancy. The *.fabgz data for the 800m and 1200m depth areas are in the native FracMan format and have been compressed using Gzip. Filtered data for the 800m depth area includes .csv spreadsheets, native FracMan (.fab), and GOCAD (.ts) files that are in a compressed zip format. The file titled "SGW 2020 Finnila and Podgorney DFN fracture files on GDR.pdf" is a description of the data and should be reviewed prior to data use.

15 GEOTHERMAL ENERGY↗

"Hybrid fracture/matrix modeling for well completion options evaluation"

Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Here we present the development of a methodology to rapidly and efficiently numerically simulate mixed fracture-matrix flow systems for evaluation of well design and completion options. The methodology is based on a loose coupling framework, allowing the fracture and matrix systems to be meshed separately. The fracture system includes the integration of fracture growth and aperture data from well stimulation simulations of stochastically generated fracture networks. Automatic mesh refinement is used in the matrix simulation to resolve heat transfer near the fracture network. This simulation framework is used to efficiently determine optimal production and injection well placement using adaptive sampling.

15 GEOTHERMAL ENERGY↗

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗

Revisiting the first Sandia Fracture Challenge with transient deformation heating and strain localization considerations

In this study, we present a systematic experimental study and accompanying theoretical analysis of the dependence of ductile fracture on strain localization, strain hardening rates, deformation-induced thermal softening, and transient heat conduction under spatially uniform as well as spatially heterogeneous deformation. Spatially uniform cases are studied via standard dogbone shaped specimens subject to uniaxial tension. Spatially heterogeneous deformation is studied via the so-called Sandia Fracture Challenge (SFC) specimen, which is a standard compact tension specimen modified with three machined holes in front of a blunt notch. We utilize the same precipitation hardened martensitic stainless steel used in the first SFC experiments, i.e. 15-5 PH with an H1075 heat treatment. We also study 15-5 PH in Condition A and with the H900 heat treatment, each of which has a different hardening behavior and ductility. We find that the ductility does not correlate with the deformation to first crack initiation in the SFC specimen. Instead, hardening rates are better correlated. Moreover, a re-examination finds that an accurate calibration of the hardening rates is strongly correlated with accurate blind predictions of ductile fracture in the SFC specimens. Given that thermal softening can greatly affect hardening rates, we provided a thermographic analysis of deformation-induced heating and transient heat transfer in both the dogbone shaped samples and the SFC specimens. Transient heating in stainless steels is found to have first-order effects on the strain at the onset of necking, strain to fracture, as well as strain localization and ductile fracture in complex geometries even under so-called quasi-static loading rates.

42 ENGINEERING↗

Microstructure of irradiated AGR TRISO particle buffer layers as measured by X-ray computed tomography

Shrinkage of the initially low-density buffer layer in tristructural isotropic (TRISO) coated fuel particles during irradiation is a well-known phenomenon with potential implications for fission product and actinide transport, as well as potential pyrocarbon fracture that in some cases has been observed to impact particle performance. During post-irradiation examination, the buffer layer's structure is commonly determined using 2D microscopy of a particle cross section or a series of particle cross sections at staggered depths. Although these methods provide a general idea of the irradiated buffer microstructure, they do not provide a full picture of the TRISO coating microstructure. By contrast, x-ray computed tomography (XCT) provides full 3D imaging of the TRISO particle. Particles from various compacts irradiated for the Advanced Gas Reactor Fuel Development and Qualification Program were imaged using XCT. Some of these particles were selected because of abnormal fission product inventories (e.g., low 137 Cs), whereas others were randomly selected from the center of the fission product inventory distribution. Qualitative and quantitative analysis techniques were applied to the randomly selected particles as representatives of typical TRISO behavior to study the post-irradiation structure of the buffer layer. These results showed that, while separation of the buffer and inner pyrolytic carbon layers was a common behavior during the radiation-induced shrinkage of the buffer, a portion of the original buffer/inner pyrolytic carbon interface remained intact throughout irradiation in nearly all cases. Furthermore, these results also indicated clear trends in the degree of buffer densification with irradiation temperature and fluence.

36 MATERIALS SCIENCE↗

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗

Water chemistry data and calculations for Fracture Sustainability Tests

Measured cation, anion, and silica are presented as well as calculations evaluating mass balance performed for 4 fracture sustainability tests under EGS conditions: Stripa granite, Metamudstone, Rhyolitic Ashflow Tuff, and Silicified Rhyolitic Tuff. Calculations of charge balance are included for water samples in the 4 tests.

15 GEOTHERMAL ENERGY↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Evaluating sealability of blended smart polymer and fiber additive for geothermal drilling with the effect of fracture opening size

Geothermal formations often contain extensive fracture networks. These fracture networks contribute to the significant loss of drilling fluids during geothermal drilling. Multiple loss circulation materials (LCM) such as fiber, granules, and pills have been proposed to tackle this problem but with only limited success. Recent advances in materials science have led to the development of thermoset shape memory polymers (SMP) to address the lost circulation problem. In this paper, we evaluate a thermoset SMP performance in sealing near wellbore fractures of different sizes in geothermal wells. The SMP performance was assessed using granite disks and cylindrical granite cores having fracture sizes of 1000 μm and 3000 μm. A static filtration test was performed using cedar fiber, CaCO 3 , and SMP. Results showed cedar fiber performed better than the CaCO 3 ., reducing fluid loss by 89% and improving sealing pressure by 200 psi. A novel dynamic testing unit that allows for high-temperature testing under flowing conditions was used in this study. The analysis showed that 3% by weight SMP and fiber blends could bridge and plug the 1000 μm fracture. For a larger fracture of 3000 μm width, there was a need to increase the weight concentration of the SMP to 6% to plug the fracture opening effectively. We showed the influence of key parameters such as the type of LCM, concentration, and particle size distribution in optimizing the performance of drilling fluid loss treatment.

02 PETROLEUM↗