Search NASA⌕ Search

Engineering topics

Cole, David R.

Publications and source records attributed to Cole, David R..

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics

Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.

58 GEOSCIENCES↗
DOE OSTI2022

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗
DOE OSTI2022

Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA

For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.

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

Sodium diffusion in heterogeneous porous media: Connecting laboratory experiments and simulations

Sodium has been suggested as a tracer for brine in reservoir formations where a significant amount of sodium ion is found mainly in the aqueous phase. Fortunately, the most abundant sodium isotope, 23 Na, is Nuclear Magnetic Resonance (NMR) active, making it possible to study the structure and dynamical properties of both bulk and pore confined fluid systems. We investigated the diffusion of key dissolved sodium species in bulk solution and porous media as a function of composition, ionic strength, and permeability of the carbonate samples using NMR relaxometry and pulsed-field gradient (PFG) NMR techniques. We use Indiana limestone as an example of natural porous media and water as a freely diffusible tracer and carrier. We demonstrate apparent diffusion measurement of sodium based on changes in spin–spin relaxation time (T 2 ) signal in pore confinement. The diffusion rate of sodium decreases with increasing counter-ion size. This effect is greater at higher ionic strengths and lower chemical potentials in porous media. The reactive transport code, CrunchFlow, was used to complement the NMR experiments to simulate diffusion behavior in porous media. Combining these two methods provides a powerful approach to estimating effective diffusion coefficients in heterogeneous matrices. The modeling considers the influence of physical properties (porosity and tortuosity) and chemical properties (geochemical composition and chemical gradients). The results presented in this work highlight the advantage of measuring apparent diffusivity using NMR T 2 relaxometry in conjunction with numerical simulation to derive effective diffusivity and the corresponding matrix properties (i.e., tortuosity) of the system.

58 GEOSCIENCES↗
DOE OSTI2022

Effects of pore connectivity and tortuosity on the dynamics of fluids confined in sub-nanometer pores

Dynamical behavior of fluids under nano-pore confinement is studied extensively as it has important implications for several industrial as well as geological processes. Pore network in many porous materials exhibits a varied degree of inter connections. The extent of this pore connectivity may affect the structural and dynamical behavior of the confined fluid. However, studies of fluid confinement addressing these effects systematically are lacking. In this paper we report molecular dynamics simulation studies addressing the effects of pore connectivity on the dynamics of two representative fluids – CO 2 and ethane in silicalite by systematically varying the degree of pore connectivity through selectively blocking some pore space with immobile methane molecules. By selectively turning off the pore spaces in the shape of straight, or tortuous zigzag channels, we also probe the effects of pore tortuosity. In general, pore connectivity is found to facilitate both the translational as well as rotational dynamics of both fluids, while the intermolecular modes of vibration in both fluids remain largely unaffected. The effects of providing connections between a set of straight or zigzag channel-like pores are however more nuanced. Pore tortuosity facilitates the rotational motion, but suppresses the translational motion of CO 2 , while its effects on the rotational and translational motion of ethane are less pronounced. The intermolecular vibrational modes of both fluids shift to higher energies with an increase in the number of tortuous pores. The results reported here provide a detailed molecular level understanding of the effects of pore connectivity on the dynamics of fluids and thus have implications for applications like fluid separation.

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

Molecular Structure of Adsorbed Water Phases in Silica Nanopores

The adsorption of water vapor in silica nanopores with different pore morphologies and surface hydrophilicities was studied to quantify the densities and thicknesses of the water sorption layers and deduce their molecular structures. Water adsorption to surface hydroxyls is described by a multilayer sorption model. At low pressure, the water adsorption isotherms are largely independent of pore size and the adsorbed amounts scale with the surface hydroxyl density. Adsorbed-phase densities corresponding to the adsorption of two water molecules per surface hydroxyl group are found in the first adsorbed water layer for a wide range of surface hydroxyl densities. The densities and layer thickness values found in narrow pores indicate that patchy adsorbed layers form if not enough water molecules exist for a full layer, which coexist with dry pore surface regions. This behavior indicates cooperative adsorption effects, i.e., a preference for the formation of hydrogen bonds between water molecules bound to surface hydroxyls. In narrow pores, pore condensation limits further growth of the sorption layer, and in larger pores and at the planar quartz surface, a second adsorbed water layer is formed, which can hold up to approximately 4 additional water molecules per surface hydroxyl group. The water molecules in these thicker adsorption layers arrange such that the sorption layer density is similar to the bulk water density. Pore confinement limitations on the sorption layers are observed in pores with radii of as large as 8 nm. Molecular dynamics modeling reveals two preferential orientations for water molecules adsorbing to surface hydroxyl groups and suggests an intralayer structuring in the adsorbed monolayer. Furthermore, adsorbed water molecules in the sorption layer are bonded to the surface hydroxyl group via the donation and acceptance of hydrogen bonds.

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

Characterization of Nanoscale Pores in Tight Gas Sandstones Using Complex Techniques: A Case Study of a Linxing Tight Gas Sandstone Reservoir

Pore structures with rich nanopores and permeability in tight gas reservoirs are poorly understood up to date. Advanced techniques are needed to be employed to accurately characterize pore structures, especially tiny pores which include micron and nanopores. In this study, various experimental techniques such as scanning electron microscopy (SEM), nuclear magnetic resonance (NMR) T 2 , nitrogen adsorption method, and NMR cryoporometry (NMRC) are combined to interrogate the complex pore systems of the tight gas reservoir in the Linxing formation, Ordos Basin, China. Results show that tight gas sandstones are primarily comprised of residual interparticle and clay-dominated pores. Clay and quartz are two dominate minerals while pyrite occupies a nontrivial amount as well. The permeability of tight gas sandstones is very low, exhibiting an extremely poor positive correlation with porosity. While pore types and relative pore contents are more influential factors on the permeability, accurate characterization of pore size distribution is critical for the permeability of tight gas sandstones. Therefore, complementary characterization methods are carried out, indicating that neither small pores with radii < 100 nm (around peak 1 in NMR T 2 distribution) nor large pores with radii > 5 μ m (around peak 3 in NMR T 2 distribution) control the permeability by analyzing the connectivity of the pores in various size ranges, but rather pores averaging approximately 350 ± X nm (around peak 2 in NMR T 2 distribution) have sufficient connectivity to host and transmit hydrocarbons. The pore size of tight gas sandstones is dominated by the clay-rich mineral assemblage. The study shows that the NMRC technique can be a very promising method, especially when referred to as a promising “roadmap” on how to interrogate tight formations such as the tight gas sands or even shale especially for the nanopore characterization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Effects of Pore Connectivity on the Sorption of Fluids in Nanoporous Material: Ethane and CO2 Sorption in Silicalite

Adsorption of fluids in nanoporous materials is important for several applications including gas storage and catalysis. The pore network in natural, as well as engineered, materials can exhibit different degrees of connectivity between pores. While this might have important implications for the sorption of fluids, the effects of pore connectivity are seldom addressed in the studies of fluid sorption. We have carried out Monte Carlo simulations of the sorption of ethane and CO2 in silicalite, a nanoporous material characterized by sub-nanometer pores of different geometries (straight and zigzag channel like pores), with varied degrees of pore connectivity. The variation in pore connectivity is achieved by selectively blocking some pores by loading them with methane molecules that are treated as a part of the rigid nanoporous matrix in the simulations. Normalized to the pore space available for adsorption, the magnitude of sorption increases with a decrease in pore connectivity. The increased adsorption in the systems where pore connections are removed by blocking them is because of additional, albeit weaker, adsorption sites provided by the blocker molecules. By selectively blocking all straight or zigzag channels, we find differences in the absorption behavior of guest molecules in these channels.

Gautam, Siddharth (ORCID:0000000314435382)↗
DOE OSTI2021

Nanoscale Structure and Dynamics in Geochemical Systems

Neutron scattering is a powerful tool to elucidate the structure and dynamics of systems that are important to geochemists, including ion association in complex aqueous solutions, solvent-exchange reactions at mineral–water interfaces, and reaction and transport of fluids in nanoporous materials. This article focusses on three techniques: neutron diffraction, which can reveal the atomic-level structure of aqueous solutions and solids; quasi-elastic neutron scattering, which measures the diffusional dynamics at mineral–water interfaces; and small-angle neutron scattering, which can show how properties of nanoporous systems change during gas, liquid, and solute imbibition and reaction. Furthermore, the usefulness and applicability of the experimental results are extended by rigorous comparison to computational simulations.

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

Comparative geochemistry of flowback chemistry from the Utica/Point Pleasant and Marcellus formations

Flowback/Produced fluid samples were collected from several wells from two Utica/Point Pleasant (UPP) sites (UPPW and UPPS) in Ohio, and one Marcellus (Marcellus Shale Energy and Environment Laboratory (MSEEL)) site in West Virginia over a period of approximately two years. Although these formations have different ages, depositional environments, diagenetic histories, and geochemical and mineralogical compositions (i.e. the UPP is significantly more carbonate rich than the Marcellus which is more siliceous), analysis of trends in fluid species over time shows that, overall, the TDS and major solubilized elements (Na, Ca, Cl) in the UPP and Marcellus brines are remarkably similar. Total dissolved solutes (TDS) in these brines ranged from approximately 40 to 250 g/L salt, and in general, concentrations increased with time elapsed since natural gas well completion and stimulation. The behavior of Na, Br, and Cl suggests that the produced water signatures from these formations are largely derived from the native formational brines which display evidence of originating from evaporated seawater. There is a strong correlation between Cl and Br, indicating that both species behave conservatively, and the similarity among each of these brines suggests no appreciable contribution of salt from halite dissolution because Br is excluded from the halite structure. Cl/Br ratios in the brines range from ~80 to 120 (mg/L/mg/L). Other elements, such as K, which readily reacts between fluids and ion exchange sites on clays, generally exhibit conservative behavior for an individual site, but show significant variations among each of the different well pads. The concentrations of Sr and Ba vary dramatically among well sites, and increase with respect to Cl- over time, suggesting increasing solubilization, presumably from desorption from clay minerals or dissolution of carbonates or sulfates from the source formation(s). The UPPW well site has very low Ba due to high-sulfate input fluid, which resulted in precipitation of barite/celestite in the brines. In contrast the UPPS well site had elevated Sr (~ 3500 mg/L), presumably due to the use of Sr-rich recycled brine used in hydraulic fracturing. The Marcellus site had the highest Ba concentrations (up to 10 g/L) and highest Ba/Sr ratios in the fluids, due to the high concentration of barium in the Marcellus target (~ 1000 ppm, as compared to ~200 ppm in the UPP). These observations suggest that solutes in the FP fluids are derived from native brines, water-rock interactions that have occurred over geologic time scales, as well as some contribution from contemporaneous reactions in the subsurface. The results also show that the composition of the injected fluid can influence flowback fluid chemistry and possibly production efficiency.

58 GEOSCIENCES↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗