Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale”

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 847 records · Page 47

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗

Method and Device for Extraction of Liquids from a Solid Particle Material

A method, system, and device for separating oil from oil sands or oil shale is disclosed. The method includes heating the oil sands, spinning the heated oil sands, confining the sand particles mechanically, and recovering the oil substantially free of the sand. The method can be used without the addition of chemical extraction agents. The system includes a source of centrifugal force, a heat source, a separation device, and a recovery device. The separation device includes a method of confining the sands while allowing the oil to escape, such as through an aperture.

deMayo, Benjamin↗

Observations of VOC Emissions and Photochemical Products over US Oil- and Gas-Producing Regions Using High-Resolution H3O+ CIMS (PTR-ToF-MS)

VOCs (Volatile Organic Compounds) related to oil and gas extraction operations in the United States were measured by H3O (sup plus) chemical ionization time-of-flight mass spectrometry (H3O (sup plus) ToFCIMS/PTR-ToF-MS (Time of Flight Chemical Ionization Mass Spectrometry/Proton Transfer Reaction-Time of Flight-Mass Spectroscopy) from aircraft during the Shale Oil and Natural Gas Nexus (SONGNEX) campaign in March-April 2015. This work presents an overview of major VOC species measured in nine oil- and gas-producing regions, and a more detailed analysis of H3O (sup plus) ToF-CIMS measurements in the Permian Basin within Texas and New Mexico. Mass spectra are dominated by small photochemically produced oxygenates and compounds typically found in crude oil: aromatics, cyclic alkanes, and alkanes. Mixing ratios of aromatics were frequently as high as those measured downwind of large urban areas. In the Permian, the H3O (sup plus) ToF-CIMS measured a number of underexplored or previously unreported species, including aromatic and cycloalkane oxidation products, nitrogen heterocycles including pyrrole (C4H5N) and pyrroline (C4H7N), H2S, and a diamondoid (adamantane) or unusual monoterpene. We additionally assess the specificity of a number of ion masses resulting from H3O (sup plus) ion chemistry previously reported in the literature, including several new or alternate interpretations.

mass spectrometr↗

Zircon, Baddeleyite, and Reidite Found in Ries Crater Suevite

Introduction: Impact events can generate superheated impact melts and even vapor [1]. In the past decade, evidence of the high temperatures and high pressures of the impact process has been found in the impact melt from terrestrial craters (e.g., Mistastin Lake [2], Meteor Crater [3]) and even on the Moon [4–5]. Often, these studies involve zircon or zirconium-bearing phases. Zircon is a particularly useful mineral, due to its robustness and durability against weathering. It is used principally for chronology, but it has a multitude of geologic applications, including geothermometry and fingerprinting magma sources. Here we focus on its capability of recording impact conditions. Grains of zircon (ZrSiO4) are converted at high temperature and/or pressure during an impact event to reidite (a high-pressure polymorph of ZrSiO4) or to tetragonal-ZrO2 + SiO2 [6]. Reidite and tertragonal-ZrO2 leave identifying relicts in the rocks, markers of the high pressures and/or temperatures these rocks underwent. Ries Crater is a 26-km-diameter peak ring crater in southern Germany that formed approximately 15 Ma [7–9]. Ries is the type locality of the polymict impact breccia known as suevite. Here, we report a microanalytical study of suevite from the Ries impact structure. These analyses are used to inform our understanding of the pressure and temperature conditions involved in creating a polymict impact breccia such as suevite in a peak ring crater. Sample Description: The sample, denoted 16RS08, originates from Otting Quarry, at 48.8777° N, 10.7921° E, approximately 17 km from the center of the Ries Crater in Germany (approximately 4 km outside the crater rim). The material recorded in sample 16RS08 is therefore considered an outer suevite. The outer suevite is a discontinuous layer of polymict impact breccia that occurs outside the central ring of Ries, up to 22 km from the center of the crater [9]. Lithic clasts in the outer suevite consist primarily of crystalline basement rocks (gneiss, granite, amphibolite), with less than 5% of lithic clasts being overlying sedimentary rocks (limestone, sandstone, shale) [9]. The thin section studied contains variably shocked lithic and mineral clasts, impact glass, and interstitial minerals that make up the matrix of the breccia (Fig. 1) [9]. Methods: We used a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona to obtain 15 elemental X-ray maps of 16RS08. We next used the JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC) to obtain electron backscatter diffraction (EBSD) maps and energy dispersive X-ray spectroscopy (EDS) maps of select portions of the section. The EBSD data were collected under beam conditions of 20 kV, and ~9 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We used the elemental X-ray maps to identify the areas of interest in the section, particularly phosphates and Zr-bearing grains. These areas of interest were then targeted for follow up EBSD and EDS analyses. We have identified, through combined EBSD and EDS analysis, the presence of zircon, reidite, and monoclinic-ZrO2 (baddeleyite) in 16RS08. In 16RS08, we have found singular grains of zircon, zircon with a vermicular baddeleyite halo (Fig. 2a), and granular zircon with reidite (Fig. 2b). The different Zr-rich phases and their corresponding textures signifies that this sample underwent a broad spectrum of pressure and temperature conditions during the impact event. For example, zircons surrounded by a vermicular baddeleyite and SiO2 intergrowth (i.e., Fig. 2a) have been shown to preserve evidence of the extremely high temperatures of impact melt, upwards of 2370 °C [2, 6]. Similarly, reidite and granular zircon (i.e., Fig. 2b) have been shown to preserve evidence of high pressure, as the transition to reidite occurs >30 GPa [3, 6]. Future Work: Next, we will process the EBSD and EDS data for these Zr-rich grains, specifically looking for indicators of cubic- or tetragonal-ZrO2 in the baddeleyite remnants and of shock-precursors to the reidite. To further inform our work, we will also obtain BSE images of these grains using a Hitachi S-4800 SEM in KMICF at the University of Arizona, as well as geochemical spot analyses via EPMA. The data collected will be used to constrain the formation conditions of the Ries Crater outer suevite. Acknowledgments: We thank Ken Domanik and Jerry Chang for their support with data collection. This work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. TME thanks A. Cavosie and N. Timms for assistance during field sampling. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Melosh H. J. (1989) Oxf. U. Press. [2] Timms et al. (2017) EPSL 477, 52–58. [3] Cavosie et al. (2016) Geology 44:9, 703–706. [4] White et al. (2020) Nature Astr. 4, 974–978. [5] Crow C. A. et al. (2017) GCA 202, 264–284. [6] Timms et al. (2017) Earth-Sci. Rev. 165, 185–202. [7] Schmieder M. et al. (2018) GCA 220, 146–157. [8] Schwarz W. H. et al. (2020) M&PS 55:2, 312–325. [9] Stöffler et al. (2013) M&PS 43:4, 515–589.

Zircon↗

Anoxic photochemical weathering of pyrite on Archean continents

Sulfur is an essential element of life that is assimilated by Earth’s biosphere through the chemical breakdown of pyrite. On the early Earth, pyrite weathering by atmospheric oxygen was severely limited, and low marine sulfate concentrations persisted for much of the Archean eon. Here, we show an anoxic photochemical mechanism of pyrite weathering that could have provided substantial amounts of sulfate to the oceans as continents formed in the late Archean. Pyrite grains suspended in anoxic ferrous iron solutions produced millimolar sulfate concentrations when irradiated with ultraviolet light. The Fe2+(aq) was photooxidized, which, in turn, led to the chemical oxidation of pyritic sulfur. Additional experiments conducted with 2.68 Ga shale demonstrated that photochemically derived ferric iron oxidizes and dissolves sedimentary pyrite during chemical weathering. The results suggest that before the rise of atmospheric oxygen, oxidative pyrite weathering on Archean continents was controlled by the exposure of land to sunlight.

Jihua Hao↗

Potentially Underestimated Gas Flaring Activities—A New Approach to Detect Combustion Using Machine Learning and NASA’s Black Marble Product Suite

Monitoring changes in greenhouse gas (GHG) emission is critical for assessing climate mitigation efforts towards the Paris Agreement goal. A crucial aspect of science-based GHG monitoring is to provide objective information for quality assurance and uncertainty assessment of the reported emissions. Emission estimates from combustion events (gas flaring and biomass burning) are often calculated based on activity data (AD) from satellite observations, such as those detected from the visible infrared imaging radiometer suite (VIIRS) onboard the Suomi-NPP and NOAA-20 satellites. These estimates are often incorporated into carbon models for calculating emissions and removals. Consequently, errors and uncertainties associated with AD propagate into these models and impact emission estimates. Deriving uncertainty of AD is therefore crucial for transparency of emission estimates but remains a challenge due to the lack of evaluation data or alternate estimates. This work proposes a new approach using machine learning (ML) for combustion detection from NASA's Black Marble product suite and explores the assessment of potential uncertainties through comparison with existing detections. We jointly characterize combustion using thermal and light emission signals, with the latter improving detection of probable weaker combustion with less distinct thermal signatures. Being methodologically independent, the differences in ML-derived estimates with existing approaches can indicate the potential uncertainties in detection. The approach was applied to detect gas flares over the Eagle Ford Shale, Texas. We analyzed the spatio-temporal variations in detections and found that approximately 79.04% and 72.14% of the light emission-based detections are missed by ML-derived detections from VIIRS thermal bands and existing datasets, respectively. This improvement in combustion detection and scope for uncertainty assessment is essential for comprehensive monitoring of resulting emissions and we discuss the steps for extending this globally.

gas flaring↗

Opportunities for Clean Energy in Natural Gas Well Operations: Preprint

The oil and gas industry is increasingly seeking operational improvements to reduce both costs and emissions. Currently, oil and gas directly and indirectly contributes forty-two percent of global greenhouse gas emissions, with over twenty percent of the industry’s emissions coming from operations. Given the opportunity for emissions reductions, this study describes techno-economic analysis evaluating opportunities for distributed energy generation and storage technologies – including solar photovoltaics (PV), distributed wind energy, and battery energy storage – to support companies’ energy cost savings targets, clean energy goals, and energy resiliency needs at hypothetical upstream well sites in the Marcellus Shale in Pennsylvania, both grid-connected and off-grid. These technologies reduce the site’s consumption of grid electricity and natural gas and thus help reduce Scope 1 and 2 emissions associated with electricity and natural gas. For each scenario, a cost of avoided emissions was calculated; these values can be compared to internal organizational value placed on emissions reductions, compared to other emissions reduction strategies such as energy efficiency, reducing flaring, and direct carbon capture and sequestration, and compared to existing (albeit limited) U.S. carbon markets such as California’s Low Carbon Fuel Standard. The study also explores the ability of these electric clean energy technologies to support site resiliency against utility outages.

42 ENGINEERING↗

Opportunities for Clean Energy in Natural Gas Well Operations

The oil and gas industry is increasingly seeking operational improvements to reduce both costs and emissions while improving resilience against electric grid outages. This study describes techno-economic analysis of opportunities for distributed energy generation and storage technologies to support companies' energy cost savings, clean energy, and energy resiliency goals. Specifically, the analysis evaluates solar photovoltaics (PV), distributed wind energy, and battery energy storage at hypothetical upstream well sites in the Marcellus Shale in Pennsylvania, both grid-connected and off-grid. Results indicate opportunity for solar PV to reduce operational costs. Additionally, these technologies reduce the site's consumption of grid electricity and natural gas and thus can help reduce Scope 1 and 2 emissions associated with electricity and natural gas consumption. For each emissions reduction scenario, a cost of avoided emissions was calculated; these values can be compared to internal organizational value placed on emissions reductions, compared to other emissions reduction strategies such as energy efficiency, reducing flaring, and direct carbon capture and sequestration, and compared to existing (albeit limited) U.S. carbon markets such as California's Low Carbon Fuel Standard. Results indicate that the associated costs of emissions reductions via distributed renewables are competitive with these options and markets. The study also explores the ability of these electric clean energy technologies to support site resiliency against utility outages.

42 ENGINEERING↗

Wellbore Leakage Mitigation Using Advanced Mineral Precipitation Strategies (Final Scientific/Technical Report)

This research project, entitled Wellbore Leakage Mitigation using Advanced Mineral Precipitation Strategies, was directed toward advancing mineral precipitation technologies to seal wellbore leaks and other potential carbon dioxide (CO 2 ) leakage pathways, including compromised wellbore cement and fracture networks of varying geometries throughout the wellbore and surrounding formation. The advanced mineral sealing techniques were evaluated in the laboratory and the field to expand the range of temperature, pressure, and chemical conditions where mineral precipitation could be used. Rates of enzymatically or thermally induced calcium carbonate precipitation (EICP, TICP) were determined as well as rates of enzyme inactivation at elevated temperatures (>60 °C). Immobilization was shown to enhance the thermal stability of urease enzyme and provide protection from inhibition in low pH conditions. Mineral seals promoted by EICP and TICP were used to seal fractured composite (cement/sandstone) and shale cores in high pressure reactors operated at temperatures greater than 55 °C and 130 °C, respectively. A computational model was updated to model processes of EICP and TICP. The work resulted in a successful field test, which showed that EICP in the presence of low pH fluids (CO 2 -impacted brine) could promote permeability reduction and seal a channel in the near-wellbore environment.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulations of Gas Production from Gas Hydrate Reservoirs at the Prudhoe Bay Unit 7-11-12 Pad on Alaska North Slope

In December 2018, a partnership between the U.S. Department of Energy National Energy Technology Laboratory (DOE NETL), the Japan Oil, Gas and Metals National Corporation (JOGMEC), and the U.S. Geological Survey (USGS) successfully drilled and logged the Hydrate-01 Stratigraphic Test Well (STW) in the greater Prudhoe Bay oil field on Alaska North Slope. The logging-while-drilling (LWD) data confirmed the presence of gas hydrate-bearing reservoirs within two sand reservoirs (Unit B and Unit D) that are suitable targets for future testing. The interpreted log data and core sample measurements were used to create reservoir models for the Prudhoe Bay Unit (PBU) Kuparuk 7-11-12 Site. The models combine both gas hydrate-bearing sections in Unit B and Unit D together with the intermediate Unit C and over- and underburden sands and shales. The vertical heterogeneity in porosity, gas hydrate saturation, and permeability distributions for reservoir and non-reservoir units was implemented using fine mesh discretization. The depressurization method was applied to Unit B to induce gas hydrate destabilization at constant bottom hole pressure (BHP) values. The results of the numerical simulations support the development of production scenarios, well design, surface facilities, and field test procedures with the main goal to perform efficient and safe scientific production testing. Project results contribute to the knowledge base of permafrost-associated gas hydrate accumulations as a future energy source.

Myshakin, Evgeniy↗

Development of Methane Emissions Model to Assess Fuel Recovery Potential at Gas Well Sites Using On Site Compression

The U.S. natural gas production and consumption has increased 85.5% since 2005 primarily due to the unconventional production methods of horizontal drilling and hydraulic fracturing. Natural gas used as a fuel has a lower greenhouse gas (GHG) footprint than coal and petroleum due to lower Carbon Dioxide (CO2) emissions when combusted. However, the “greener” benefit to natural gas may be negated by leaks in production and transmission systems. Methane (CH4), the primary hydrocarbon in natural gas, has an estimated Global Warming Potential (GWP) of 28-36 over 100 years, meaning it can absorb 28-36 more energy than CO2 which has a GWP of 1.0. Natural gas well sites are prone to methane emissions, or leaks and irregular gas releases, vented to atmosphere throughout production and transmission. The U.S. Department of Energy (DOE) and the National Energy Technology Laboratory (NETL) has recently granted West Virginia University (WVU) funding under agreement DE-FOA-0002005, to “Advance technologies to mitigate methane emissions and increase the efficiency of the natural gas transportation infrastructure”. As part of this funding WVU was tasked with identifying and quantifying sources of methane emissions at unconventional well sites, processing this data, and developing a system to recapture these emissions. A 0-D Simulink model was developed, utilizing standardized methodologies, data from previously conducted studies, as well as collected data from well sites in the Marcellus shale play region. The model was developed to predict emission rates from various components at natural gas well sites as well as the potential to utilize these emissions as fuel for the natural gas powered compressor engines on-site. This model was utilized to run high, medium, and low cases for four identified emission sources, engine size, pneumatic controller count, liquid level production which dictates tank emissions, and compressor packing vent emissions. Due to discrepancies in transient tank emission data, a high and low emission factor for tanks was used, resulting in two sets of 81 executed cases, and 162 unique cases of total site emissions and potential for fuel consumption. Each of the cases were run over 86,400 seconds at a 1 Hz, representative of a full 24 hour day of operation. The fuel consumption offset an average of 557% of fuel consumption on an energy density basis across all 81 cases with the high tank emission factor with a maximum offset of 2334%. The fuel consumption offset was an average of 82.9% for all 81 cases with the low tank emission factor with a maximum offset of 337%. This study highlights flaws in the use of publicly available methane number calculations to determine natural gas’s suitability as an engine fuel as well as the lack of public data for transient liquid storage tank emissions.

03 NATURAL GAS↗

Methane Emissions Produced by Pneumatic Devices and Produced Water Tanks on Natural Gas Wellsites

In recent years, more scrutiny has been placed on the release of greenhouse gasses (GHGs) by the oil and gas industry in North America as climate concerns increase. In the coming years, stricter regulations proposed by the Environmental Protection Agency (EPA) concerning the release of GHGs will further direct attention to GHG emissions quantification, modeling, and overall reduction. Active natural gas wellsites employ production equipment that routinely emits methane during normal operation. Production equipment often includes, but is not limited to, gas production units (GPUs), pneumatic controllers (PCs) and actuators powered by produced natural gas, and produced water storage tanks. This work details direct methane emission measurements taken at active natural gas wellsites originating from intermittent pneumatic controllers and produced water storage tanks. Efforts to model methane flows originating from these sources are also detailed, and those results are compared to direct measurement. The direct methane emissions measurement campaign collected flow data from two active dry wellsites in the Marcellus shale region. These measurement campaigns captured, in total, several weeks of continuous data from these sites. Data included multiple pneumatic device vent flowrate channels, a single channel monitoring total emissions vented from the produced water storage tanks, and corresponding site weather data. Periodic composition measurements were taken from flows vented from produced water tanks. Computational modeling efforts attempted to mimic a given site’s methane emissions behavior given limited input parameters. Direct measurement data and simulation data yielded, through reduction and analysis, apparent methane emissions factors, statistics on pneumatic device actuation events, important relationships between GPU parameters and pneumatic device actuation events, produced water storage tank vent flows and flow compositions, and parameters most contributing to the magnitude and frequency of these fugitive methane emissions. Measurements, considering data from a 72hr period, produced PC emissions factors from 0.008 to 0.03 SCFH per device with an average device emitting 0.026 SCFH at the first site. The second site produced PC emissions factors spanning from 0.05 to 25.47 SCFH with an average device emitting 1.93 SCFH not including the 25.47 SCFH device, which was assumed to be malfunctioning. Two unique water tank emission measurements were made, 618 g CH4 per hour and 32.6 g CH4 per hour. On a water basis, the water tank emissions factors were 0.66 and 0.042 kg CH4 per bbl of water, respectively.

03 NATURAL GAS↗

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

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

Benge, Margaret↗

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗