Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Oxidation of soil organic carbon during an anoxic-oxic transition

Redox reactions of iron (Fe) can play an important role in controlling the stability and transformation of organic carbon (OC) in soils. However, there is limited knowledge about the dynamics of Fe and OC in soils during anoxic-oxic transitions, which are common in temperate and tropical biomes. In this study, we investigated the fate and transformation of Fe and OC during a 5-day oxic incubation of four pre-reduced (8 days anoxic) forest soils. Up to 3.7% of total OC in the pre-reduced soils was oxidized to CO 2 , in conjunction with the oxidation of 31 –84% of the extractable Fe(II). OC oxidation in pre-reduced soils was lower than non-reducedcontrol soils, indicating that recently anoxically-incubated soils support lower OC metabolism upon re-oxidation than consistently oxic soils. Oxic OC oxidation was negatively correlated with the fraction of Fe-bound OC and also the fraction of short-range-ordered (SRO) Fe oxides. These results suggest that association with Fe oxides—especially SRO Fe oxides—can inhibit the availability of OC for oxidation. However, the difference between the oxidation of OC in pre40 reduced soils and consistently oxic controls could not be linked to the moderate changes in the Fe mineral composition during the anoxic-oxic transition, indicating the importance of longer term OC association with Fe oxides, as well as other processes, in regulating the oxidation of OC during the redox fluctuations. OC association with Fe oxides should be considered during simulations of biogeochemical cycling of OC in soils during the redox fluctuations, for which other processes such as the anoxic transformation of OC and response of the microbial community also need to be incorporated.

58 GEOSCIENCES↗
DOE OSTI2021

Field-scale estimation of soil properties from spectral induced polarization tomography

Estimates of soil properties such as Cation Exchange Capacity (CEC), water content, grain size characteristics, and permeability are important in geotechnical engineering, water resources, and agriculture. We develop a non-intrusive approach to estimate these properties in the field using spectral induced polarization (SIP) tomography. This geophysical method provides information about the frequency dependence of the complex electrical conductivity of porous media. Using 18 soil samples collected from a Bordeaux vineyard, we first conducted a laboratory study using SIP over the frequency range 10 mHz-45 kHz. The laboratory data were used to confirm the accuracy of a recently developed dynamic Stern layer petrophysical model. The results are consistent with published values from previous works using soils. A comparison was made by comparing the field complex conductivity spectra and the experimental data at two locations where core samples were obtained. The model was then used in concert with field data to image the spatial distribution of CEC, water content, permeability, and mean grain size along a vineyard transect. For clay and sandy textures found in the field, measured and estimated CEC agree rather well (from 6 to 40% discrepancy). Furthermore, our approach provides an efficient way to estimate important soil properties in a non-invasive manner, in high resolution, and over field-relevant scales of the critical zone of the Earth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Soil organic matter is principally root derived in an Ultisol under oak forest

To a large degree, the sources and stability of soil organic carbon remain poorly constrained. A clear understanding of links among the components of the soil C cycle is hampered by the complexity of the system as well as challenges associated with partitioning bulk soil C into meaningful fractions. A large accidental 14 CO 2 release at the Oak Ridge Reservation in Tennessee, USA provided a strong label pulse into adjacent, well-studied oak forests, resulting in highly elevated Δ 14 C values in leaf litter (~1000‰) and roots (~260–450‰). A four-year manipulative study was conducted to determine the relative contribution of litter versus roots to the bulk mineral soil C pool, as well as to free light, occluded light and heavy fractions. The heavy fraction was further split into fractions with densities of 1.7–2.4 g cm –3 and >2.4 g cm –3 to test the homogeneity of the mineral-associated fraction of C. Substantial concentrations of label were detected in all soil fractions within a year of the 14 CO 2 release, indicating rapid incorporation of newly fixed photosynthates in all fractions of soil organic C, regardless of differences in stability inferred by previous work. This rapid incorporation of label occurred only in treatments where roots were labeled, indicating that roots are the major source of inputs to mineral soil C stocks at these sites. Separation of the heavy fraction into subfractions of intermediate (1.7–2.4 g cm –3 ) and high (>2.4 g cm –3 ) density indicated that both subfractions incorporated label at similar rates, despite significant differences in degree of microbial processing. In general, the rate of label incorporation suggested a much faster turnover for all fractions than indicated by natural radiocarbon abundance values. This study suggests that within each soil fraction there are portions of slow-cycling and fast-cycling materials, and the determination of an average turnover time or mean age is dependent on experimental approach. The rapid incorporation of label into all fractions within a year regardless of inferred stability implies a high degree of heterogeneity in all fractions regardless of how finely the soils are partitioned. Further refinement of the nature and drivers of this heterogeneity could yield important insights into the soil C cycle.

58 GEOSCIENCES↗
DOE OSTI2021

Quantifying the influence of different biocrust community states and their responses to warming temperatures on soil biogeochemistry in field and mesocosm studies

Biocrusts influence soil biogeochemistry by fixing carbon (C) and nitrogen (N) and through leachate inputs to soils. Functional rates can vary among biocrust community states and in response to edaphic properties, heterotrophic microbial activity, and global change. Using soils and biocrusts from the Colorado Plateau, Utah, USA, we aimed to quantify the influence of early-successional (ES) and late-successional (LS) biocrust community states on soil biogeochemistry. In a field setting, we found soil was less “fertile” under ES than LS biocrust, but ES biocrust had a relative influence 1.3 times greater than LS biocrust on soil fertility. Leachate collected from LS biocrust had, on average, 6 times more organic C and 1.7 times more dissolved N than ES biocrust, but concentrations of phosphorus (P) and inorganic N did not differ among the two biocrust types. Here, to disentangle influences of biocrusts and soil properties on biogeochemical pools, we constructed mesocosms from homogenized soil and left the surface bare or covered with ES or LS biocrust, before assignment to ambient or warmed (+5 °C) temperature treatments for 3 months. Multivariate biogeochemical properties differed among cover types, yet all exhibited losses of P, N, and organic C and nearly half of the biogeochemical variables considered did not differ among cover types. Mesocosms with LS biocrust retained more dissolved N, supported 8 additional, significant correlations among biogeochemical pools of the biocrust and mineral soil layer on average, and lost fewer of these correlations under warming. Overall, while soils under LS biocrusts were more fertile (i.e., had higher nutrient concentrations) than under ES, we did not find evidence implicating leachate as the primary driver of this difference. Biocrust influences on soil fertility were greater when mineral soil nutrients were in lower concentrations, highlighting the value of even incipient biocrusts for dryland functioning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Chronological changes in soil biogeochemical properties of the glacier foreland of Midtre Lovénbreen, Svalbard, attributed to soil-forming factors

We report that Glacier forelands provide an excellent opportunity to investigate vegetation succession and soil development along the chronosequence; however, there are few studies on soil biogeochemical changes from environmental factors, aside from time. This study aimed to investigate soil development and biogeochemical changes in the glacier foreland of Midtre Lovenbreen, Svalbard, by considering various factors, including time. Eighteen vegetation and soil variables were measured at 38 different sampling sites of varying soil age, depth, and glacio-fluvial activity. Soil organic matter (SOM) was quantitatively measured, and the compositional changes in SOM were determined following size-density fractionation. In the topsoil, the soil organic carbon (SOC) and total nitrogen (N) content was found to increase along the soil chronosequence and were highly correlated with vegetation-associated variables. These findings suggest that plant-derived material was the main driver of the light fraction of SOM accumulation in the topsoil. The heavy fractions of SOM were composed of microbially transformed organic compounds, eventually contributing to SOM stabilization within short 90-yr deglaciation under harsh climatic conditions. In addition to time, the soil vertical profiles showed that other environmental parameters, also affected the soil biogeochemical properties. The high total phosphorous (P) content and electrical conductivity in the topsoil were attributed to unweathered subglacial materials and a considerable amount of inorganic ions from subglacial meltwater. The high P and magnesium content in the subsoil were attributed to parent materials, while the high sodium and potassium content in the surface soil were a result of sea-salt deposition. Glacio-fluvial runoff hampered ecosystem development by inhibiting vegetation development and SOM accumulation. This study emphasizes the importance of considering various soil-forming factors, including parent/subglacial materials, aeolian deposition, and glacio-fluvial runoff, as well as soil age, to obtain a comprehensive understanding of the ecosystem development in glacier forelands.

glacio-fluvial runoff↗
DOE OSTI2022

Contributions of anoxic microsites to soil carbon protection across soil textures

Anoxic microsites, zones of oxygen depletion in otherwise oxic soils, may slow soil C turnover. However, the abundance of anoxic microsites and their contribution to soil C protection is yet undefined. In this study, we determine the contribution of anoxic microsites to soil C protection in soils of three distinct textures (clay loam, loam, sandy loam) across a range of soil moistures. We examined the influence of soil oxygen supply by increasing oxygen content in the incubation atmosphere (oxygen enrichment) and through disaggregation. We attributed increases in CO 2 efflux to the aeration of anoxic microsites. The contribution of anoxic microsites to soil C protection increased with decreasing clay content. Clay loam CO 2 efflux was relatively unaffected by aeration. Moderately moist, loam soils had CO 2 effluxes that did not increase with oxygen enrichment but increased by 375% upon disaggregation. Sandy loam soil CO 2 efflux increased by 50–75% with oxygen enrichment and 250% with disaggregation. Geochemical and microbial data reveal that anoxic microsite abundance also increased with decreasing clay content. The proportion of acid extractable Fe present as Fe(II) increased with decreasing clay content, and methanogens were more abundant in loam and sandy loam soils. Our results suggest that oxygen demand, rather than supply, can regulate anoxic microsite formation and that anoxic protection of soil C can be diminished through physical disruption of soil structure. Our findings further illustrate that anoxic microsites should be included in conceptual models of soil C protection to avoid soil C loss and improve predictions of soil C response to disturbance.

58 GEOSCIENCES↗
DOE OSTI2023

An interlaboratory comparison of mid-infrared spectra acquisition: Instruments and procedures matter

Diffuse reflectance spectroscopy has been extensively employed to deliver timely and cost-effective predictions of a number of soil properties. However, although several soil spectral laboratories have been established worldwide, the distinct characteristics of instruments and operations still hamper further integration and interoperability across mid-infrared (MIR) soil spectral libraries. In this study, we conducted a large-scale ring trial experiment to understand the lab-to-lab variability of multiple MIR instruments. By developing a systematic evaluation of different mathematical treatments with modeling algorithms, including regular preprocessing and spectral standardization, we quantified and evaluated instruments' dissimilarity and how this impacts internal and shared model performance. We found that all instruments delivered good predictions when calibrated internally using the same instruments' characteristics and standard operating procedures by solely relying on regular spectral preprocessing that accounts for light scattering and multiplicative/additive effects, e.g., using standard normal variate (SNV). When performing model transfer from a large public library (the USDA NSSCKSSL MIR library) to secondary instruments, good performance was also achieved by regular preprocessing (e. g., SNV) if both instruments shared the same manufacturer. However, significant differences between the KSSL MIR library and contrasting ring trial instruments responses were evident and confirmed by a semi-unsupervised spectral clustering. For heavily contrasting setups, spectral standardization was necessary before transferring prediction models. Non-linear model types like Cubist and memory-based learning delivered more precise estimates because they seemed to be less sensitive to spectral variations than global partial least square regression. In summary, the results from this study can assist new laboratories in building spectroscopy capacity utilizing existing MIR spectral libraries and support the recent global efforts to make soil spectroscopy universally accessible with centralized or shared operating procedures.

58 GEOSCIENCES↗
DOE OSTI2025

The role of Ca-bridged organic matter in an alkaline soil, as revealed by multimodal chemical imaging

Mineral–organic matter (OM) studies have predominantly focused on acidic soils that are abundant in iron (Fe) oxides and aluminum (Al) oxides. We have probed mineral–OM interactions in an alkaline or calcareous soil of the Aridisols class. Unlike the role of Fe and Al, the role of Ca-minerals (particularly calcite), which are ubiquitous in alkaline soils, in OM sequestration is not well understood. Multiple recent model studies with aqueous Ca2+ or synthetic calcite and a suite of OM compounds have shown Ca-OM assemblages to be spatially correlated with calcite at the microscale. To study the chemical state of both Ca and Fe and their competing role in soil organic matter (SOM) stabilization, we performed laboratory characterization using x-ray diffraction, Mössbauer spectroscopy, x-ray photoelectron spectroscopy, scanning electron microscopy, and scanning transmission electron microscopy, alongside synchrotron-based microscale chemical imaging using scanning transmission x-ray microscopy combined with near-edge x-ray absorption fine structure. Ca mineral–organic associations were found to be ubiquitous in this system and are likely critical for understanding SOM stabilization/degradation in alkaline soils. From our findings on mineralogy, speciation, and the nature of Ca-OM bridging, we identified differences in C and Ca chemistry based on the relative location of OM to Ca minerals. The OM near the calcite crystal was enriched in lipid and protein moieties, Ca-OM next to Fe minerals displayed a strong contribution from aromatic compounds, while on the surface of microbes, the carbonate was believed to be of microbial in origin, as also suggested by preliminary works reporting on the formation of amorphous calcite or nano-calcite. In Ca-OM admixed with carbonate, it was difficult to distinguish Ca-associated OM from amorphous calcite or nano-calcite.

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

Composition and persistence of soil organic matter along eroding and depositional transects in buried vs. modern soil layers: A case of the Brady paleosol at Wauneta, Nebraska

Paleosols form when soils are buried through deposition by aeolian, colluvial, alluvial or other processes. Burial of former topsoil isolates soil organic matter (SOM) from surface conditions, allowing carbon to accumulate and potentially remain stable for millennia. In this study, SOM composition, distribution, and persistence were analyzed in the Brady Soil of Nebraska, USA to compare SOM spatial variability in modern and buried soils, as well as the impact of erosional exposure on SOM stability. The Brady Soil, formed as a surface soil during the Pleistocene-Holocene transition and now a paleosol buried up to 6 m deep (or more) by loess deposition during the Holocene, was sampled along burial (up to 5.8 m depth) and erosional (up to 1.8 m depth) transects to compare SOM dynamics in different geomorphic settings. Fourier Transform Infrared Spectroscopy (FTIR) and Fourier Transform ion cyclotron resonance mass spectrometry (FTICR-MS) were used to analyze SOM composition, while δ 13 C isotope analyses identified SOM sources and radiocarbon values were used to estimate turnover rates. Results confirmed a vegetation shift from C3 to C4 plants after Brady Soil formation, reflecting warming climatic conditions. Increasing SOM age and decreasing δ 13 C and δ 15 N values with depth indicated slowing of decomposition rate in buried soils. Higher pH in the Brady Soil suggested greater base cation content, supporting SOM stabilization through organo-mineral associations and aggregate formation. However, exposure of the Brady Soil due to surface erosion caused faster SOM turnover. This result suggested susceptibility of buried SOM to losses via decomposition upon erosional exposure, possibly accelerated by priming in response to modern SOM inputs. These findings highlight the potential loss of carbon stocks in buried soils under future climate change, as shifts in soil physicochemical properties may destabilize long-preserved SOM.

Environmental sciences↗
DOE OSTI2023

Economic and operational investigation of CO 2 sequestration through enhanced oil recovery in unconventional reservoirs in Colorado, USA

The ongoing CCUS commercial projects are highly relied on the support of government incentives due to massive capital investment. Here, this study analyzes the economics of carbon capture utilization and sequestration (CCUS) projects, shows a state-wide CCUS deployment exercise, followed by simulation results of enhanced oil recovery (EOR) based CO 2 storage in unconventional reservoirs. The comprehensive economic analysis of capture, transportation, sequestration costs, enhanced 45Q tax credits, and EOR revenue implies the practicality of CO 2 -EOR to offset the high CCUS costs. With the economics understanding, we study the top CO 2 sources, existing CO 2 pipelines, and sequestration sinks in the state of Colorado, USA. This paper next presents results from EOR simulation in one section of the unconventional Denver-Julesburg (DJ) Basin Niobrara and Codell reservoirs. The simulation model is based on a geological static model, incorporated with hydraulic fracture stimulation, history matched to production, and calibrated to the microseismic and time-lapse surface seismic data. The CO 2 -EOR simulation results show that oil production can be increased and more CO 2 stored with: a longer primary production period; the presence of a shut-in period; higher injection rates; and multi-well injectors. The modeling results show that about 7–10 Mscf of CO 2 will be stored when recovering 1 stb of EOR oil. By adding the enhanced oil revenue and the carbon credits together, it is estimated that the most economic case can generate $\$13$ MM when oil price is assumed to be $\$80$/stb, and the EOR oil revenue is 3.4 times greater than that generated from 45Q incentives. It corresponds to the scenario that a five-year primary production is followed by CO 2 injection into four wells with the sequence of injection (4 MMscf/day for 6 months), shut-in (6 months) and production (12 months). The best practices in this study will provide valuable insights for similar CCUS projects in other unconventional fields. Furthermore, this study defines a term named “Carbon Neutrality Index (CNI)” by comparing the amount of CO 2 stored with that burned by EOR oil. The CNI value of 0 indicates the enhanced oil is carbon neutral; a negative CNI value implies there is a net reduction in carbon emission. The 4-year huff-n-puff (HnP) simulation leads to a positive CNI value, indicating that the EOR oil generated in this process is not carbon neutral yet.

03 NATURAL GAS↗
DOE OSTI2023

Investigation of methane mass transfer and sorption in Marcellus shale under variable net-stress

Natural gas in shale exists as free and adsorbed gas, subject to prevailing pore pressures and stress conditions. Accordingly, to accurately estimate/predict the shale gas recovery potential, a central requirement is to represent gas transport and sorption behavior under varying stress conditions. The objective of this work is to facilitate the interpretation of laboratory-scale experiments, at relevant conditions, in an attempt to bridge the gap in scales between laboratory- and field-scale observations. We have conducted a series of high-pressure experiments on a full-diameter core sample from the Marcellus shale. These include gas loading (pressure-decay) and depletion (production) experiments with pure methane (CH 4 ) at variable stress conditions to characterize transport and sorption behavior under reservoir-relevant conditions. Here we have formulated and applied a novel integral model for mass transfer and storage in multi-porosity shale systems that allows us to effectively investigate transport and sorption phenomena: We delineate gas transport by interpreting helium (He) pressure-decay experiments and demonstrate how to use the information gained to calculate the relevant transport coefficients of CH 4 and other gases. A separate measurement of the CH 4 sorption isotherm on a smaller sample (a shale cube) was interpreted and combined with the transport description to predict the production behavior of CH 4 from the experiments with the full-diameter core. Our experiments demonstrate that the representation of sorption hysteresis is crucial for predicting and guiding shale gas production: At the end of both gas production experiments, approximately 20% of the initial gas in place remained in the core. Without accounting for sorption hysteresis, our modeling demonstrates that the CH 4 production could be overestimated by 10%. We demonstrate that our integral, triple-porosity model provides an effective approach for the interpretation and prediction of gas transport and sorption behavior during loading and production experiments on shale cores under variable net-stress conditions. In summary, our work combines measurements and modeling of mass transfer and sorption in shales at different scales to validate a characterization approach that facilitates an improved understanding of shale gas production. Furthermore, the triple-porosity model utilized in our work defines a potential pathway for the translation of laboratory-scale experimentation to larger-scale applications.

58 GEOSCIENCES↗
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗
DOE OSTI2023

Using high pressure solutions of polyfluoroacrylate and CO 2 to Seal cement cracks for improved wellbore integrity

Polyfluoroacrylate (PFA) is a hydrophobic and oleophobic polymer that is soluble in high pressure carbon dioxide (CO 2 ). In this study, the ability of PFA-CO 2 solutions to greatly reduce the apparent permeability of split or cracked Portland cement cylindrical samples is assessed. The apparent permeability values of confined samples were determined before and after treatment with PFA-CO 2 solutions. In four tests, PFA-CO 2 solutions were continuously displacing pure CO 2 from the cracked cement and the decrease in apparent permeability due to PFA adsorption and wettability alteration was monitored. The lowest apparent permeability cracked cement sample (81 nD) was completely sealed with a very small amount of solution. Samples with initial apparent permeabilities of 89 μD and 29.4 mD exhibited 92% and 99% reductions in permeability, respectively, before the experiments had to be stopped because of the excessively large increase in pressure drop. A 50% reduction in apparent permeability was observed with a 3.80 mD sample. Four other split cement samples (bound together with tape) with an initial apparent permeability in the 9.0–70 mD range were removed from the core holder and immersed in a PFA-CO 2 solution for 24 h to allow for PFA adsorption. Then the PFA-CO 2 solution was depressurized, allowing for the deposition of additional PFA from the solution within the crack as the pressure fell below the cloud point pressure of the PFA-CO 2 solution. These four samples were then confined again in a core holder and apparent permeability reductions of 29–93% were observed. Finally, results from these eight experiments indicates that the more substantial reductions in the nD – mD apparent permeability of the cracked cement correlated to lower initial crack permeability, higher PFA concentration, and slower injection rate of the PFA-CO 2 solution into the crack.

58 GEOSCIENCES↗
DOE OSTI2023

Uncertainty quantification of the convolutional neural networks on permeability estimation from micro-CT scanned sandstone and carbonate rock images

Rock permeability is one of the most crucial properties affecting subsurface fluid flow behaviors. To accurately and robustly estimate the permeability, Digital Rock Physics, including micro-CT scanning technology and direct flow simulations on scanned images, has prevailed in recent years. Besides, machine learning techniques such as convolutional neural networks (CNNs) have been widely adopted and achieved success in permeability estimations directly from rock images. However, existing ML methods used for permeability estimation from rock images lack uncertainty quantification that causes unreliable predictions and overconfident estimations on out-of-distribution (OOD) samples. Here, in this work, we propose a PI3NN-CNN framework to address this problem. PI3NN-CNN consists of a CNN model for absolute permeability estimation and a PI3NN method to quantify the estimation uncertainty. It is able to quantify the uncertainty for in-distribution (InD) data with a desired confidence level, and identify OOD samples to avoid overconfident predictions. We demonstrate the method using micro-CT scanned images from two sandstone and two carbonate rocks. We found that PI3NN-CNN generates accurate predictions for InD samples, while producing high-quality prediction uncertainties regardless of the prediction accuracy. Meanwhile, PI3NN-CNN identifies OOD samples using its special network initialization scheme. The unique feature of PI3NN-CNN makes it applicable to more complex real-world image-based data for robust learning and predictions without overconfident estimations when the ground-truth information is unavailable.

58 GEOSCIENCES↗
DOE OSTI2024

Stochastic reconstruction of fracture network pattern using spatial point processes

Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗