SEARCH · Engineering Papers
Results for “SANDSTONES”
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.
CO 2 charged brines changed rock strength and stiffness at Crystal Geyser, Utah: Implications for leaking subsurface CO 2 storage reservoirs
Explore the source record for details and available documents.
Assignment and calibration of relative permeability by hydrostratigraphic units for multiphase flow analysis, case study: CO 2 -EOR operations at the Farnsworth Unit, Texas
Included in the most critical factors for geological CO 2 storage site screening, selection, and operation is effective simulations of multiphase flow and transport. Relative permeability is probably the greatest source of potential uncertainty in multiphase flow simulation, second only to intrinsic permeability heterogeneity. The specific relative permeability relationship assigned greatly impacts forecasts of CO 2 trapping mechanisms, phase behavior, and long-term plume movement. A primary goal of this study is to evaluate the impacts and implications of different methods of assigning relative permeability relationships for CO 2 -EOR model forecasts. Most simulation studies published in the literature base selection of relative permeability functions on the geologic formation or rock type alone. In this study, we initially implemented reservoir model grids with previously-identified hydrostratigraphic units based on porosity and permeability relationship of the Morrow ‘B’ Sandstone, then assigned relative permeability functions for those hydrostratigraphic units. Specific, constrained relative permeability relationships were created and assigned to each hydrostratigraphic unit using petrophysical data and Mercury Intrusion Capillary Pressure (MICP) measurements, from core samples of each hydrostratigraphic unit. Results of forward simulations with the newly-calibrated models will be compared to those of previous models as well as to simulation results for a range of different relative permeability relationships. The study site is the Farnsworth Unit (FWU) in the northeast Texas Panhandle, an active CO 2 -EOR operation. The target formation is the Morrow ‘B’ Sandstone, a clastic formation composed of medium to course sands.
Low Pore Connectivity in Natural Rock
As repositories for CO₂ and radioactive waste, as oil and gas reservoirs, and as contaminated sites needing remediation, rock formations play a central role in energy and environmental management. The connectivity of the rock's porespace strongly affects fluid flow and solute transport. This work examines pore connectivity and its implications for fluid flow and chemical transport. Three experimental approaches (imbibition, tracer concentration profiles, and imaging) were used in combination with network modeling. In the imbibition results, three types of imbibition slope [log (cumulative imbibition) vs. log (imbibition time)] were found: the classical 0.5, plus 0.26, and 0.26 transitioning to 0.5. The imbibition slope of 0.26 seen in Indiana sandstone, metagraywacke, and Barnett shale indicates low pore connectivity, in contrast to the slope of 0.5 seen in the well-connected Berea sandstone. In the tracer profile work, rocks exhibited different distances to the plateau porosity, consistent with the pore connectivity from the imbibition tests. Injection of a molten metal into connected pore spaces, followed by 2-D imaging of the solidified alloy in polished thin sections, allowed direct assessment of pore structure and lateral connection in the rock samples. Pore-scale network modeling gave results consistent with measurements, confirming pore connectivity as the underlying cause of both anomalous behaviors: imbibition slope not having the classical value of 0.5, and accessible porosity being a function of distance from the edge. A poorly connected porespace will exhibit anomalous behavior in fluid flow and chemical transport, such as a lower imbibition slope (in air–water system) and diffusion rate than expected from classical behavior.
An experimental feasibility study on the use of CO 2 -soluble polyfluoroacrylates for CO 2 mobility and conformance control applications
We report high molecular weight polyfluoroacrylate (PFA) is an amorphous, sticky, hydrophobic, oil-phobic, CO 2 -soluble homopolymer. 1 wt% PFA-in-CO 2 solutions are four times as viscous as pure CO 2 . Although these solutions provided modest mobility control, unexpectedly large increases in pressure drop occurred because a portion of the PFA adsorbed onto the pore surfaces, altering wettability and significantly permeability. Further, PFA exhibited diminished solubility in CO 2 in the presence of light extracted hydrocarbons. Although these effects rendered CO 2 -PFA solutions inappropriate for mobility control, excellent conformance control for dual parallel brine-saturated cores was attained when a CO 2 -PFA solution was first injected into the isolated high permeability sandstone core. This core was then placed in parallel with a low permeability sandstone core, and all of the subsequently injected CO 2 was diverted to the low permeability core. Conformance control was not as effective in limestone, possibly due the low pH-induced erosion of the limestone that mitigated PFA adsorption.
Machine learning-assisted upscaling analysis of reservoir rock core properties based on micro-computed tomography imagery
Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.
Effects of Mineral Surface Properties on Supercritical CO 2 Wettability in a Siliciclastic Reservoir
Wettability is a key reservoir characteristic influencing geological carbon sequestration (GCS) processes, such as CO 2 transport and storage capacity. Wettability is often determined on a limited number of reservoir samples by measuring the contact angle at the CO 2 /brine/mineral interface, but the ability to predict this value remains a challenge. In this work, minerals comprising a natural reservoir sample were identified, and the influence of their surface roughness and mineralogy on the contact angle was quantified to evaluate predictive models and controlling mechanisms. The natural sample was obtained from the Mount Simon formation, a representative siliciclastic reservoir that is the site of a United States Department of Energy CO 2 injection project. A thin section of the Mount Simon sandstone was examined with compound light microscopy and environmental scanning electron microscopy (ESEM) coupled with energy-dispersive X-ray spectroscopy (EDS). Quartz and feldspar were identified as dominant minerals and were coated with various reddish black precipitates consistent with illite clay and iron oxide hematite. Contact angle (θ) measurements were conducted for the four representative minerals and the Mount Simon sample over a range of pressures (2–25 MPa) at 40 °C. At supercritical conditions, all samples are strongly water-wet, with contact angles between 27° and 45°. Several predictive models for contact angle were evaluated for the mineral and Mount Simon samples, including the Wenzel and Cassie–Baxter models, plus newly proposed modifications of these that account for the fraction of different minerals comprising the reservoir sample surface, the surface roughness, and the extent that roughness pits are filled with brine. Modeling results suggest that the fraction of mineral surfaces containing roughness pits filled with brine is the most important reservoir characteristic that controls wettability in the Mount Simon sandstone, followed by surface mineralogy. Finally, to our knowledge, this is one of the few studies to investigate the effects of individual minerals on the wettability of a natural reservoir sample.
Investigation of Mass Transfer and Sorption in CO 2 /Brine/Rock Systems via In Situ FT-IR
CO 2 geological storage in deep saline formations is considered a promising method to mitigate anthropogenic CO 2 emissions and, thereby, minimize changes to the Earth’s atmosphere. A fundamental understanding of CO 2 mass transfer and sorption phenomena in brine-saturated reservoir formations is necessary to understand the long-term fate of injected CO 2 as it is subjected to different (physical, dissolution, and mineral) trapping mechanisms. In this work, we investigate CO 2 sorption in brine-saturated and dry Mt. Simon sandstone samples via in-situ Fourier Transform infrared spectroscopy (FT-IR) at elevated pressures, ranging from 0.3 MPa to 8.3 MPa, at a temperature of 50 ⁰C. The FT-IR spectra of bulk-phase CO 2 were simultaneously recorded under the same conditions. For bulk-phase CO 2 , we observed, in agreement with past studies, a doublet peak at 2361 cm -1 and 2336 cm -1 and another peak (ν 2 bending mode) at 667 cm -1 . With increasing pressure, the position of the peak at 667 cm -1 remains invariant, however, when crossing into the supercritical region the doublet peak degenerates onto a single peak at 2336 cm -1 with a barely visible shoulder at 2361 cm -1 . The bulk CO 2 data provide a perfect fit for Beer’s law for the whole range of pressure conditions. For the dry sample, the IR spectrum is experimentally indistinguishable from the bulk CO 2 spectrum, signifying that if physical adsorption occurs to any significant extent, the adsorbed CO 2 molecules are not substantially more rotationally constrained than the dense bulk CO 2 molecules. For the brine-saturated sample, we observe a strong band centered at 2342 cm -1 and a small companion peak at 2360-2361 cm -1 that degenerates into a barely visible shoulder peak at the higher pressures. The 2342 cm -1 band has been previously observed by other investigators for CO 2 dissolved in bulk water/brine as well during its adsorption on a variety of other wet natural porous media. We observe no peaks corresponding to bicarbonate or carbonate bulk species, which correlates well with the prior literature on similar low-pH aqueous solutions. The integrated peak area for the CO 2 sorbed in the brine-saturated sample correlates linearly with its solubility in the same bulk brine, as measured separately via a PVT-cell approach. Here, this validates the accuracy of both techniques, and the potential of the FT-IR method to be used in the study of mass transfer and adsorption in such systems. To that effect, a simple mathematical model is presented to analyze the FT-IR data to determine the CO 2 effective diffusivity in the brine-saturated sandstone sample.
A Quantum Mechanical Study of Quartz (101) Interfacial Boundaries
Density functional theory (DFT) calculations were performed on periodic and molecular models to explore the energetics of bonding modes at quartz–quartz interfaces, common grain contacts in sandstones. Four interface types were modeled: H 2 O-mediated, silanol H-bond mediated, O–O peroxide bonds, and siloxane bonds. Each type of interface may exist in the subsurface of the Earth, depending upon water activity, temperature, and pressure, as interfacial structures are stable within DFT molecular dynamics simulations up to 473 K. The results predict interfacial energies SiOSi > SiOOSi > SiOH-6H 2 O > 2SiOH, as expected based on similar bond strengths in other systems. The proposed SiOOSi interface is not generally accounted for in chemomechanical models but could be important in sandstones at depth in the Earth. The implications for intergrain fracturing and H2 generation are discussed.
Time-Dependent Consolidation in Porous Geomaterials at In Situ Conditions of Temperature and Pressure
Analysis of quartz sandstones shows that grain-scale crushing (fracture and rearrangement) and associated sealing of fractures contribute significantly to consolidation. The crushing strength (P*) for granular material is defined by laboratory experiments conducted at strain rates of 10 -4 to 10 -5 s -1 and room temperature. Based on experiments, many sandstones would require burial depths in excess of the actual maximum burial depth to create observed microstructure and density. In this paper, we use experiments and soil mechanics principles to determine rate laws for brittle consolidation of fine-grained quartz sand to better estimate in situ failure conditions of porous geomaterials. Experiments were conducted on St. Peter sand utilizing different isostatic consolidation and creep load paths at temperatures to 200 °C and at strain rates of 10 -4 to 10 -10 s -1 . Experiment results are consistent with observed rate dependence of consolidation in soils, and P* for sand can be identified by the change in the dependence of consolidation rate with stress, allowing the extrapolation of P* determined in the laboratory to geologic rates and temperatures. Additionally, normalized P* values can be described by a polynomial function to quantify temperature, stress, and strain-rate relationships for the consolidation of porous geomaterials by subcritical cracking. At geologic loading rates, P* for fine-grained quartz sand is achieved within ~3-km burial depth, and thus, shear-enhanced compaction under nonisostatic stress can occur at even shallower depths. Finally, these results demonstrate that time and temperature effects must be considered for predicting the brittle consolidation of sediments in depositional basins, petroleum reservoirs, and engineering applications.
Experimental Poroviscoelasticity of Common Sedimentary Rocks
The success of geoenergy applications such as petroleum recovery or geological storage of CO 2 depends on properly addressing the physical coupling between the pore fluid diffusion and mechanical deformation of the subsurface rock. Constitutive models should include short-term hydromechanical interactions and long-term behavior and should incorporate the principles behind the mathematical models for poroelastic and poroviscoelastic responses. However, the viscous parameters in constitutive relationships still need to be validated and estimated. In this work, we experimentally quantify the time-dependent response of fluid-filled sedimentary rocks at room temperature and isotropic stress states. Drained, undrained, and unjacketed geomechanical tests are performed to measure the poroelastic parameters for Berea sandstone, Apulian limestone, clay-rich material, and Opalinus clay (shale). A poroviscous model parameter, the bulk viscosity, is included in the constitutive relationships. The bulk viscosity is estimated under constant isotropic stress conditions from time-dependent deformation of rock in the drained regime for timescales ~10 5 s and from observations of the pore pressure growth under undrained conditions at timescales of ~10 4 s. The bulk viscosity is on the order of 10 15 –10 16 Pa s for sandstone, limestone, and shale and ~10 13 Pa s for clay-rich material, and it decreases with an increase in pore pressure despite a corresponding decrease in the effective stress. In the long term, fluid pressure can asymptotically approach minimum principal stress, which in natural reservoirs may lead to liquefaction or rock embrittlement, causing slip instabilities and earthquakes and creating high-permeability channels in low-permeable rock.
Grain Size Variations in the Murray Formation: Stratigraphic Evidence for Changing Depositional Environments in Gale Crater, Mars
The lowermost exposure of the Murray formation in Gale crater, Mars, was interpreted as sediment deposited in an ancient lake based on data collected by the Curiosity rover. Constraining the stratigraphic extent and duration of this environment has important implications for the paleohydrology of Gale. Insights into early Martian environments and paleofluid flow velocity can be obtained from grain size in rocks. Visual inspection of grain size is not always available for rocks investigated at field sites on Mars due to limited image coverage. But grain sizes can also be estimated from the Gini Index Mean Score, a grain-size proxy that uses point-to-point chemical variations in ChemCam Laser Induced Breakdown Spectroscopy data. Here, new Gini Index Mean Score results indicate that the Murray formation is dominated by mudstones with grains smaller than the spatial resolution of all rover cameras. Intervals of fine to coarse sandstone also are present, some of which are verified using observations of grain size and sedimentary structures in associated images. Overall, results demonstrate that most of the Murray consists of mudstone, suggesting settling of grains from suspension in low energy depositional environments such as lakes. Some of the mudstones contain desiccation cracks indicating periods of drying with a lowering of lake water level. However, beds and lenses of cross-bedded sandstones are common at specific intervals, suggesting episodes of fluvial and possibly eolian deposition. The persistence of lacustrine deposits interspersed with fluvial deposits suggests that liquid water was sustained on the Martian surface for tens of thousands to millions of years.
Dynamic Processes of CO 2 Storage in the Field: 1. Multiscale and Multipath Channeling of CO 2 Flow in the Hierarchical Fluvial Reservoir at Cranfield, Mississippi
A consistent picture of dynamic channeling, invasion, spreading, and breakthrough (CISB) of supercritical CO 2 in the hierarchical fluvial reservoir at Cranfield, Mississippi is presented after 10 years of integration and analysis of complementary field monitoring and characterization data. The dynamic CISB with small-scale CO 2 -flow channels in the F1-F2-F3 cross section (F1, F2, and F3 are one injection and two monitoring wells) was imaged by daily electrical resistance tomography (ERT) and time-lapse crosswell seismic surveys. One, three, and four CO 2 flow channels logged at F1, F2, and F3, respectively, were dynamically connected with strong temporal variations in CO 2 saturation during 221 days of drainage with injection rate doubling twice and 81 days of imbibition. Three intermediate-scale CO 2 flow channels (with highest CO 2 saturation) normal to the cross section were ERT-imaged during late-time drainage. A large-scale, sinuous fluvial CO 2 flow channel was imaged by repeat surface seismic survey at the end of the imbibition. The fluvial sandstone channel sinuously bypasses the F1-F2-F3 cross section in a point bar, but the channel is connected to the cross section through an intermediate-scale sandstone channel, forming a complicated flow channel network. The multiscale flow channel network (in the fluvial channel-point bar system) revealed from the observed CISB enables us to consistently interpret the hydrological monitoring data of three tracer tests, each conducted during an injection rate step, and preinjection hydraulic-thermal-tracer tests. Finally, this interpretation of the CISB and flow channel network can guide future modeling and data inversion to best understand the effects of natural heterogeneity on CO 2 storage efficiency and residual trapping.
The Curiosity Rover's Exploration of Glen Torridon, Gale Crater, Mars: An Overview of the Campaign and Scientific Results
The Mars Science Laboratory rover, Curiosity, explored the clay mineral-bearing Glen Torridon region for 1 Martian year between January 2019 and January 2021, including a short campaign onto the Greenheugh pediment. The Glen Torridon campaign sought to characterize the geology of the area, seek evidence of habitable environments, and document the onset of a potentially global climatic transition during the Hesperian era. Curiosity roved 5 km in total throughout Glen Torridon, from the Vera Rubin ridge to the northern margin of the Greenheugh pediment. Curiosity acquired samples from 11 drill holes during this campaign and conducted the first Martian thermochemolytic-based organics detection experiment with the Sample Analysis at Mars instrument suite. The lowest elevations within Glen Torridon represent a continuation of lacustrine Murray formation deposits, but overlying widespread cross bedded sandstones indicate an interval of more energetic fluvial environments and prompted the definition of a new stratigraphic formation in the Mount Sharp group called the Carolyn Shoemaker formation. Glen Torridon hosts abundant phyllosilicates yet remains compositionally and mineralogically comparable to the rest of the Mount Sharp group. Glen Torridon samples have a great diversity and abundance of sulfur-bearing organic molecules, which are consistent with the presence of ancient refractory organic matter. The Glen Torridon region experienced heterogeneous diagenesis, with the most striking alteration occurring just below the Siccar Point unconformity at the Greenheugh pediment. Results from the pediment campaign show that the capping sandstone formed within the Stimson Hesperian aeolian sand sea that experienced seasonal variations in wind direction.
Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin
Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.
Constraining Erosion Rates and Landscape Evolution With In Situ 10 Be and 26 Al Cosmogenic Nuclides at Table Mountain, Antarctica
Abstract This study investigates surface weathering and sediment preservation at Table Mountain, a high‐elevation, hyperarid, polar landscape in the Transantarctic Mountains. We report cosmogenic nuclide concentrations ( 10 Be and 26 Al) in quartz from bedrock surfaces, erratic boulder lag, and cobbles embedded within Sirius Group sediments to quantify erosion rates. In situ 10 Be and 26 Al depth profiles from a 2.95 m permafrost core in the Sirius Group further constrain surface erosion rates and elucidate landscape stability. Measured 10 Be and 26 Al concentrations from two sandstone bedrock surfaces adjacent to Sirius Group sediments give erosion rates of 0.18–0.28 m/Myr. An erratic sandstone boulder within the lag above the Sirius Group yields erosion rates of ∼0.42 ± 0.03 m/Myr, whereas two cobbles embedded within the Sirius Group yield higher rates of 0.81–1.12 m/Myr. Depth profiles of in situ 10 Be and 26 Al indicate no vertical mixing of Sirius Group permafrost since deposition. Depth profile models are best explained by erosion rates of 0.53 +0.13 / −0.12 m/Myr, and an exposure age of 0.78 +0.06 / −0.08 Ma. We view the model “age” to represent the ∼0.8‐million‐year time‐scale for surface lowering equivalent to one attenuation length of cosmic ray production to achieve steady‐state conditions. Continual exhumation of embedded clasts from within the Sirius Group results in an accumulation of clasts forming the observed erosional lag deposit covering the landscape. Our erosion rates of the Sirius Group surface based on in situ 10 Be and 26 Al depth profiles are an order‐of‐magnitude larger than those based on meteoric 10 Be infiltration and further clarification is required.
Quantification of Dissolved O 2 in Bulk Aqueous Solutions and Porous Media Using NMR Relaxometry
Effects of dissolved paramagnetic oxygen (O 2 ) in water on 1 H nuclear magnetic resonance (NMR) Carr-Purcell-Meiboom-Gill (CPMG) experiments is evaluated at a 1 H Larmor frequency of 2 MHz. Dissolution of O 2 into water significantly reduces the 1 H transverse relaxation coefficient (T 2 ). For deoxygenated water, T 2 is 3388 ms, water at ambient atmospheric conditions (7.4 mg/L O2) exhibits a T2 of 2465 ms, and dissolution of 2710 mg/L O 2 further reduces T 2 to 36 ms. The results were fit with an empirical model to facilitate prediction of T 2 times for bulk water as a function of paramagnetic oxygen concentrations in solution. Dissolved O 2 also greatly influences 1 H NMR CPMG experiments of confined water in a model system composed of Berea sandstone. For this system, 90 mg/L O 2 in H 2 O enhances T 2 relaxation of bulk water such that the relaxation time is comparable to physically confined water in the sandstone pores. Given the sensitivity of NMR T 2 coefficients to paramagnetic oxygen, low-field NMR-based characterization of fluid and porous media structure requires control of dissolved oxygen, as geospatial variation in the partial pressure of O 2 alone is expected to perturb fluid and pore relaxation times by up to 60 and 36%, respectively.
Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks
Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.