Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Well Log”

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 739 records · Page 41

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗

Dependence of Z-R Relations on the Rain Type Classification Scheme

The TRMM Global Validation Program is giving us a unique opportunity to compare radar datasets from different sites since they are analyzed in a relatively uniform procedure. Monthly Ze-R relations for four different sites (i.e, Melbourne Florida, Houston Texas, Darwin Australia and Kwajalein Atoll) were derived. The relations were obtained using the Window Probability Matching Method (WPMM). This version of the PMM relies on matching unconditional probabilities of rain rates, R, and radar reflectivity, Ze, using rain gauge and radar data, respectively. This procedure was done separately for convective and stratiform rain type using the Steiner classification procedure. The radar and gauge data from all sites were quality controlled using the same algorithms, which include also an automatic procedure to filter unreliable rain gauge data upon comparison to radar data. An adjusted power law Z-R for each rain type was also derived by comparing the radar-gauge coincident pairs in order to adjust the total monthly rainfall to match the gauges. The obtained PMM based Ze-R relations are found to be curved lines in log-log space rather than any straight line power law. While the PMM based Ze-R curves were always distinctly different between the convective and stratiform rain, the power law based Z-R, in few cases, was found to be the same for both types. In general, a given reflectivity was matched to a much lower rain intensity in the convective rainfall as compared to that in stratiform rainfall. These findings are inherently contradictory to previous findings based on disdrometer data and suggest some precaution for using the latter Z-R relations on radar data when the partition of stratiform and convective rainfall amount is in concern. The inverse trends in the relations might be caused by effects such as partial beam fillings, the use of different classification schemes, as well as having a distinct difference in the Z-R relations between the initial convective and the trailing transition regions as suggested by recent findings.

Amitai, Eyal↗

Orbital debris and meteoroid population as estimated from LDEF impact data

Examination of LDEF's various surfaces shows numerous craters and holes due to hypervelocity impacts of meteoroids and man-made orbital debris. In this paper, the crater numbers as reported by Humes have been analyzed in an effort to understand the orbital debris and natural meteoroid environment in LEO. To determine the fraction of man-made to natural impacts, the side to top ratio of impacts and results of the Chemistry of Micrometeoroids Experiment are used. For craters in the 100 micron to 500 micron size range, about 25 percent to 30 percent of the impacts on the forward-facing surfaces and about 10 percent of the impacts on the trailing surfaces were estimated due to man-made orbital debris. A technique has been developed to convert crater numbers to particle fluxes, taking the fact into account that the distributions of impact velocity and incidence angle vary over the different surfaces of LDEF, as well as the ratio of the surface area flux to the cross-sectional area flux. Applying this technique, Humes' data concerning craters with limiting lip diameters of 100 micron, 200 micron and 500 micron have been converted into orbital debris and meteoroid fluxes ranging from about 20 micron to 200 micron particle diameter. The results exhibit good agreement with orbital debris model and meteoroid model. The converted meteoroid flux is slightly larger than Grun's model (by 40 to 70 percent). The converted orbital debris flux is slightly lower than Kessler's model for particle diameter smaller than about 30 micron and slightly larger than the model for particle diameter larger than about 40 micron. Taking also into account the IDE data point at about 0.8 micron particle diameter, it suggests to change the slope log (flux) versus log (diameter) of orbital debris flux in the 1 micron to 100 micron particle diameter range from 2.5 to 1.9.

Zhang, Jingchang↗

Is a data set distributed as a power law? A test, with application to gamma-ray burst brightnesses

We present a method to determine whether an observed sample of data is drawn from a parent distribution that is pure power law. The method starts from a class of statistics which have zero expectation value under the null hypothesis, H(sub 0), that the distribution is a pure power law: F(x) varies as x(exp -alpha). We study one simple member of the class, named the `bending statistic' B, in detail. It is most effective for detection a type of deviation from a power law where the power-law slope varies slowly and monotonically as a function of x. Our estimator of B has a distribution under H(sub 0) that depends only on the size of the sample, not on the parameters of the parent population, and is approximated well by a normal distribution even for modest sample sizes. The bending statistic can therefore be used to test a set of numbers is drawn from any power-law parent population. Since many measurable quantities in astrophysics have distriibutions that are approximately power laws, and since deviations from the ideal power law often provide interesting information about the object of study (e.g., a `bend' or `break' in a luminosity function, a line in an X- or gamma-ray spectrum), we believe that a test of this type will be useful in many different contexts. In the present paper, we apply our test to various subsamples of gamma-ray burst brightness from the first-year Burst and Transient Source Experiment (BATSE) catalog and show that we can only marginally detect the expected steepening of the log (N (greater than C(sub max))) - log (C(sub max)) distribution.

Wijers, Ralph A. M. J.↗

The SEDs and Host Galaxies of the Dustiest GRB Afterglows

The afterglows and host galaxies of long gamma-ray bursts (GRBs) offer unique opportunities to study star-forming galaxies in the high-z Universe, Until recently, however. the information inferred from GRB follow-up observations was mostly limited to optically bright afterglows. biasing all demographic studies against sight-lines that contain large amounts of dust. Aims. Here we present afterglow and host observations for a sample of bursts that are exemplary of previously missed ones because of high visual extinction (A(sub v) (Sup GRB) approx > 1 mag) along the sight-line. This facilitates an investigation of the properties, geometry and location of the absorbing dust of these poorly-explored host galaxies. and a comparison to hosts from optically-selected samples. Methods. This work is based on GROND optical/NIR and Swift/XRT X-ray observations of the afterglows, and multi-color imaging for eight GRB hosts. The afterglow and galaxy spectral energy distributions yield detailed insight into physical properties such as the dust and metal content along the GRB sight-line as well as galaxy-integrated characteristics like the host's stellar mass, luminosity. color-excess and star-formation rate. Results. For the eight afterglows considered in this study we report for the first time the redshift of GRBs 081109 (z = 0.97S7 +/- 0.0005). and the visual extinction towards GRBs 0801109 (A(sub v) (Sup GRB) = 3.4(sup +0.4) (sub -0.3) mag) and l00621A (A(sub v) (Sup GRB) = 3.8 +/- 0.2 mag), which are among the largest ever derived for GRB afterglows. Combined with non-extinguished GRBs. there is a strong anti-correlation between the afterglow's metals-to-dust ratio and visual extinction. The hosts of the dustiest afterglows are diverse in their properties, but on average redder(((R - K)(sub AB)) approximates 1.6 mag), more luminous ( approximates 0.9 L (sup *)) and massive ((log M(sup *) [M(solar]) approximates 9.8) than the hosts of optically-bright events. We hence probe a different galaxy population. suggesting that previous host samples miss most of the massive. chemically-evolved and metal-rich members. This also indicates that the dust along the sight-line is often related to host properties, and thus probably located in the diffuse ISM or interstellar clouds and not in the immediate GRB environment. Some of the hosts in our sample. are blue, young or of small stellar mass illustrating that even apparently non-extinguished galaxies possess very dusty sight-lines due to a patchy dust distribution. Conclusions. The afterglows and host galaxies of the dustiest GRBs provide evidence for a complex dust geometry in star-forming galaxies. In addition, they establish a population of luminous. massive and correspondingly chemically-evolved GRB hosts. This suggests that GRBs trace the global star-formation rate better than studies based on optically-selected host samples indicate, and the previously-claimed deficiency of high-mass host galaxies was at least partially a selection effect.

Kruhler, T.↗

Quantifying Distribution System Resilience From Utility Data: Large Event Risk and Benefits of Investments

We focus on blackouts in electric distribution systems that have a large cost to customers. To quantify resilience to these events, we show how to calculate risk metrics from the historical outage data routinely collected by utilities' outage management systems. Risk is defined using a customer cost exceedance curve. The exceedance curve has a heavy tail that implies large fluctuations in large blackout costs, and this makes estimating the mean large cost in the usual way impractical. To avoid this problem, we use new resilience metrics describing the large event risk; these metrics are the probability of a large cost event, the annual log cost resilience index, and the average of the logarithm of the cost of large-cost events or the slope magnitude of the tail on a log–log exceedance curve. Resilience can be improved by planned investments to upgrade system components or speed up restoration. The benefits that these investments would have had if they had been made in the past can be quantified by “rerunning history” with the effects of the investment included, and then recalculating the large event risk to find the improvement in resilience. An example using utility data shows a 2% reduction in the probability of a large cost event due to 10% wind hardening and 6%–7% reduction due to 10% faster restoration in two different areas of a distribution utility. This new data-driven approach to quantify resilience and resilience investments is realistic and much easier to apply than complicated approaches based on modeling all the phases of resilience. Moreover, an appeal to improvements to past lived experience may well be persuasive to customers and regulators in making the case for resilience investments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Core-Mantle Fractionation of Carbon in Earth and Mars: the Effects of Sulfur

Constraining carbon (C) fractionation between silicate magma ocean (MO) and core-forming alloy liquid during early differentiation is essential to understand the origin and early distribution of C between reservoirs such as the crust-atmosphere, mantle, and core of Earth and other terrestrial planets. Yet experimental data at high pressure (P)-temperature (T) on the effect of other light elements such as sulfur (S) in alloy liquid on alloy-silicate partitioning of C and C solubility in Fe-alloy compositions relevant for core formation is lacking. Here we have performed multi-anvil experiments at 6–13 GPa and 1800–2000 °C to examine the effects of S and Ni on the solubility limit of C in Fe-rich alloy liquid as well as partitioning behavior of C between alloy liquid and silicate melt (𝐷𝐷(sub C)(sup alloy/silicate)). The results show that C solubility in the alloy liquid as well as 𝐷𝐷(sub C)(sup alloy/silicate) decreases with increasing in S content in the alloy liquid. Empirical regression on C solubility in alloy liquid using our new experimental data and previous experiments demonstrates that C solubility significantly increases with increasing temperature, whereas unlike in S-poor or S-free alloy compositions, there is no discernible effect of Ni on C solubility in S-rich alloy liquid. Our modelling results confirm previous findings that in order to satisfy the C budget of BSE, the bulk Earth C undergoing alloy-silicate fractionation needs to be as high as those of CI-type carbonaceous chondrite, i.e., not leaving any room for volatility-induced loss of carbon during accretion. For Mars, on the other hand, an average single-stage core formation at relatively oxidized conditions (1.0 log unit below IW buffer) with 10-16 wt.% S in the core could yield a Martian mantle with a C budget similar to that of Earth’s BSE for a bulk C content of ~0.25-0.9 wt.%. For the scenario where C was delivered to the proto-Earth by a S-rich differentiated impactor at a later stage, our model calculations predict that bulk C content in the impactor can be as low as ~0.5 wt.% for an impactor mass that lies between 9-20% of present day Earth’s mass. This value is much higher than 0.05-0.1 wt.% bulk C in the impactor predicted by Li et al. (2016) because C-solubility limit of 0.3 wt.% in a S-rich alloy predicted by their models is significantly lower than the experimentally derived C-solubility of ~1.6 wt. % for the relevant S-content in the core of the impactor.

Carbon↗

North American Vegetation Dynamics Observed with Multi-Resolution Satellite Data

North American vegetation has been discovered to be a net carbon sink, with atypical behavior of drawing down more carbon from the atmosphere during the past century. It has been suggested that the Northern Hemisphere will respond favorably to climate warming by enhancing productivity and reducing the impact of fossil fuel emissions into the atmosphere. Many investigations are currently underway to understand and identify mechanisms of storage so they might be actively managed to offset carbon emissions which have detrimental consequences to the functioning of ecosystems and human well being. This paper used a time series of satellite data from multiple sensors at multiple resolutions over the past thlrty years to identify and understand mechanisms of change to vegetation productivity throughout North America. We found that humans had a marked impact to vegetation growth in half of the six selected study regions which cover greater than two million km2. We found climatic influences of increasing temperatures, and longer growing seasons with reduced snow cover in the northern regions of North America with forest fire recovery in the Northern coniferous forests of Canada. The Mid-latitudes had more direct land cover changes induced by humans coupled with climatic influences such as severe drought and altered production strategies of rain-fed agriculture in the upper Midwest, expansion of irrigated agriculture in the lower Midwest, and insect outbreaks followed by subsequent logging in the upper Northeast. Vegetation growth over long time periods (20+ years) in North America appears to be associated with long term climate change but most of the marked changes appear to be associated with climate variability on decadal and shorter time scales along with direct human land cover conversions. Our results document regional land cover land use change and climatic influences that have altered continental scale vegetation dynamics in North America.

Neigh, Christopher S. R.↗

H-Canyon Jumper Gasket Sealing Evaluation

H-Canyon Documented Safety Analysis approved December 2019 credits the leak integrity of H-Canyon (HCA) process tanks, piping, and jumpers in preventing the release of radioactive material in a post seismic environment. Demonstration of leak tightness capability for this equipment was achieved through substantial analytical effort. However, one potential weak point that affects leak integrity exists at the gasketed joint between Hanford Connectors (HC) and nozzles. Review of the HCA Crane Log over a three-year period indicates the HCs periodically leak. Savannah River National Laboratory (SRNL) was requested to perform gasket compression tests that evaluate gasket relaxation characteristics, i.e., creep, and provide improvements to the HC installation process that can reduce or eliminate the frequency of leaks at HCs. A HC torque value of 312 ft-lbf is the target value used because it is the minimum torque value, with 95% confidence, when an unlubricated HC is impacted for standard sealing period of 9 seconds. In the HCA pipe-jumper connections are installed remotely by using HCs. HCs inherently add to the uncertainty of a gasketed joint because a calibrated torque value cannot be applied. The gasket material is a special material made of asbestos-fiber impregnated with Teflon. This material has been used for decades and generally performs well but does leak intermittently as gaskets creep and preloaded torque is lost. Gasket types are either snap-ring (SR) or lampshade (LS). Both gasket types cover the same seating area. The SR gasket is designed to remain with the jumper when the jumper is removed leaving the nozzle faces clean of any debris. SR gaskets are initially installed when the jumper is fabricated. Changing gaskets is time consuming and the radioactive environment presents risks to personnel. When jumpers are removed LS gaskets could be used to better seal the connection to reduce down time and ALARA concerns. That is, an LS gasket would be remotely positioned over the nozzle and the HC reinstalled. Testing also included two situations not typical of plant operation. One situation was the use of the non-standard double-layered gaskets (DG), and the other situation was the use of a LS without removal of the installed SR gasket. The DG was used when the supply of LS gaskets ended, so that testing could continue. It was of interest to demonstrate this gasket’s creep response because the DG, which was two gaskets glued together, was thicker than the LS gasket and the thicker material was expected to elicit a greater creep characteristics. The second situation of using two gaskets together by installing an LS gasket without removing the installed gasket is not allowed by the ASME Codes. However, by demonstrating the creep history of such a gasket arrangement it could be considered as an option that could possibly seal a connection while significantly reducing operational down time and personnel risk. Forty tests were conducted, and the overall results demonstrate the initial applied force and how the gasket changes with time. The change in sealing force, in terms of torque, was then measured to obtain a force history as the gaskets creep. This report documents the results of those tests as well as the method used to perform the tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Implementation of a Time-domain Cosmic-Ray-Muon Tagger for the NEXUS Low-background Cryogenic Facility

Quantifying the effects of radiation on the operation of qubits both as quantum information systems as well as particle detectors has emerged as a pressing issue in quantum science in recent years. We present an overview of the design, operation, and deployment of a 90-$\mathrm{in}^2$ three-panel muon detector for use in the NEXUS experimental facility at Fermilab to temporally isolate correlated errors in qubits and determine if they possess an astrophysical origin. Constructed with three scintillator-attached PMTs read out with NIM modules in a triply-coincident logic scheme, we measure a surface-level muon luminosity of 9.7425 muons per second---consistent with an average surface-level cosmic-ray flux of approximately one muon per square centimeter per minute. Integrating the NIM modules with an MCC 128 DAQ HAT and Raspberry Pi, a Python script records exactly when a muon struck the detector and writes a timestamp to a log file for follow-up cross referencing. This experiment will broadly contribute to further studies aimed at understanding the source and mitigation of information loss in qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Colorado River Basin Water Resources: Pairing Phenology Citizen Science Observations with NASA Earth Observations to Monitor Riparian Vegetation

Native plant species in the Colorado River Basin are currently threatened by non-native species that are changing resource competition dynamics. Monitoring native species is essential for understanding changes in plant health throughout the year as well as how species’ distributions change inter-annually. However, current monitoring of species is limited to local efforts and cannot be efficiently conducted across more extensive regions of the basin. Project partners, including the USA National Phenology Network, United States Geological Survey, and Educating Children Outdoors, provided data and on-the-ground expertise to assist the project. NASA Earth observations, including Landsat 8 Operational Land Imager (OLI), Aqua and Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and Sentinel-2 MultiSpectral Instrument (MSI) from the European Space Agency, provided remote sensing data that covers the large spatial extent of the basin. This study combined both data sources over the study period of 2013 – 2019 and explored the relationship between on-the-ground phenology observations and remotely sensed phenology signals. The team found that pairing the datasets required spatial coordinates to be accurately logged and that it was necessary to have consistent observations throughout the year by citizen scientists. The team found promising relationships between phenology observations and the two vegetation indices, normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), when study locations were in an area with a high proportion of vegetation within a satellite pixel. The results and methods outlined in this study create new opportunities to leverage remote sensing alongside citizen science data for widespread monitoring of non-native species.

Tyler Pantle↗

Implementation of a Time-domain Cosmic-ray-muon Tagger for the NEXUS Low-background Cryogenic Facility

Quantifying the effects of radiation on the operation of qubits both as quantum information systems as well as particle detectors has emerged as a pressing issue in quantum science in recent years. We present an overview of the design, operation, and deployment of a 90-$\mathrm{in}^2$ three-panel muon detector for use in the NEXUS experimental facility at Fermilab to temporally isolate correlated errors in qubits and determine if they possess an astrophysical origin. Constructed with three scintillator-attached PMTs read out with NIM modules in a triply-coincident logic scheme, we measure a surface-level muon luminosity of 9.7425 muons per second---consistent with an average surface-level cosmic-ray flux of approximately one muon per square centimeter per minute. Integrating the NIM modules with an MCC 128 DAQ HAT and Raspberry Pi, a Python script records exactly when a muon struck the detector and writes a timestamp to a log file for follow-up cross referencing. This experiment will broadly contribute to further studies aimed at understanding the source and mitigation of information loss in qubits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗

Aging jets from low-mass stars

An extended faint optical jet is associated with the compact emission region plus faint star known as HH 55. HH 55 is located in the Lupus 2 cloud 2 min SW of the well studied T Tauri star RU Lupi. The HH 55 jet extends 55 sec N and 35 sec S in PA 160 deg. The HH 55 star is an emission line star of spectral type M3.5. Its image in the emission lines of H-alpha and (S II) is slightly elongated by 2 sec - 3 sec to the S but in continuum light is symmetrical and pointlike ((full width at half maximum) (FWHM) = 1.7 sec). The star and jet have several features in common with the star and jet known as Sz 102 = Th 28 in the nearby Lupus 3 cloud. We suggest that these objects are representative of the late evolutionary stage of the HH jet-outflow phenomenon and point out that such objects may be quite common although difficult to detect. With L(sub bol) approximately = 0.005 solar luminosity, and log T(sub e) approximately = 3.5, the HH 55 star is close to the main sequence and evolutionary tracks suggest an age of 3 x 10(exp 7) yr.

Graham, J. A.↗

A spectrophotometric study of the impact of pH and metal-to-ligand ratio on the speciation of the Pu( vi )-oxalate system

The oxalate ligand is prevalent throughout the nuclear fuel cycle. While the Pu(III)- and Pu(IV)-oxalate systems are well studied due to their use in plutonium metal and PuO 2 production, the effect of oxalate on Pu(VI) remains understudied. Absorption spectroscopy was employed to probe the solution behavior of the Pu(VI)-oxalate system as a function of pH (1, 3, 7) and metal-to-ligand ratio (M/L; 10 : 1–1 : 10). Peak changes in the UV-vis-NIR spectra were associated with the formation of multiple Pu(VI)-oxalate species with increasing oxalate concentration. Some insight into identification of species present in solution was gained from the limited Pu(VI)-oxalate literature and comparisons with the assumed isostructural U(VI)-oxalate system. A peak in the UV-vis-NIR spectrum at 839 nm, which corresponds to the formation of a 1 : 1 PuO 2 (C 2 O 4 )(aq) complex, was observed and used to determine the formation constant (log β° = 4.64 ± 0.06). A higher coordinated Pu(VI)-oxalate peak at 846 nm was tentatively assigned as the 1 : 2 complex PuO 2 (C 2 O 4 ) 2 2₋ and a preliminary formation constant was determined (log β° = 9.30 ± 0.08). In conclusion, the predominance of both complexes was shown in speciation diagrams calculated from the formation constants, illustrating the importance of considering the Pu(VI)-oxalate system in the nuclear fuel cycle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deterministic Discrete Fracture Network (DFN) Model for the EGS Collab Project on the 4850 Level of the Sanford Underground Research Facility (SURF)

The EGS Collab is conducting hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF) located in Lead, South Dakota. A total of eight ~60m-long subhorizontal boreholes were drilled from the 4850 Level (~1.5 km below the ground surface) into the crystalline rock of this former mine. Six of these holes are used for geophysical monitoring, one is used for hydraulic fracture stimulation, and the remaining hole was designed as a production borehole that receives water from the injection well via the induced and natural fracture system. The primary goal of creating the discrete fracture network model is to show that these modeling methods are critical for the development of enhanced geothermal systems (EGS). This includes the prediction of rock behavior during fracturing and during an extended period of water flow between the parallel injection and production boreholes. Understanding the results from the induced fracturing and flow is complicated by the presence of significant natural fractures that interact with the stimulation and/or flow pathways. The delineation and characterization of natural fractures is thus an important part of the project, and therefore a model of the Discrete Fracture Network (DFN) was developed on a deterministic basis. The DFN was populated using observations and interpretations integrated from drift (horizontal passageways that allow access in the underground) fracture mapping, analysis of core recovered from the eight boreholes, borehole televiewer logs and videos, and observations of flow between and within boreholes and in the drift. The natural fracture system is dominated by a pervasive northwest-trending, steeply dipping shear system that is identifiable in the drifts and the core. Hydraulic fracture stimulation, flow/tracer circulation tests, and geophysical monitoring revealed that the behavior of the injected water, and perhaps the growth of induced fractures, has been significantly influenced by the existing fractures identified in the DFN.

40 EE - Geothermal Technologies Office (EE-4G)↗

Utah FORGE 2439: Well 16B(78)-32 Field-Test Data from Mini-Frac Tests

This submittal includes the field-test data collected during stress tests conducted in the Utah FORGE 16B(78)-32 wellbore to measure/characterize the stresses in the geothermal reservoir. The type of stress test performed is referred to as a mini-frac test or a micro-frac test. The test is a hydraulic fracture test that involves injecting a small volume of water into a short interval of the reservoir that is isolated by a straddle packer to create a fracture. The test provides information about the minimum and maximum horizontal stress at the test depths. A total of seven mini-frac tests were performed within the upper (vertical or nearly vertical) section of the 9-5/8-inch diameter uncased wellbore, between depths of 5,202 ft MD [5,202 ft TVD] and 5,980 ft MD [5966.18 ft TVD]. Two types of data were generated as part of the mini-frac tests. The first type of data is time-series data for each mini-frac test (e.g., injection-rate, test-interval pressure), and are contained in the "Mini-Frac Test Data.zip" file below. The second type is geophysical log data, attached below in the "Geophysical Logs.zip" file. Here, two types of image logs and acoustic logs were obtained before and after the mini-frac tests to examine each test depth for visual evidence of induced fracture(s) and to measure the azimuth (orientation) of any observed fractures. The attached Readme file includes the test depth information and contents of each data file included here. Analysis of the field-test data was conducted and is documented in the attached report below.

15 GEOTHERMAL ENERGY↗