Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Snow distribution patterns revisited: A physics-based and machine learning hybrid approach to snow distribution mapping in the sub-Arctic Supporting Data

Snow in the Arctic and sub-Arctic is highly variable at fine scales, with deep drifts and shallow scoured areas creating a complex pattern of snow on the landscape. This fine-scale variation in snow is driven primarily by landscape and vegetation properties. Some landscape features, such as river beds, will rapidly fill in with snow during the wintertime due to high winds, while shrubs will trap blowing snow, resulting in drifts. Meanwhile, snow will blow off of exposed areas, resulting in abnormally shallow snow. These complex interactions between wind, vegetation, and terrain are difficult to represent well with physically-based models, but machine learning techniques have shown promise in the past. Here, we propose a hybrid modeling approach, where we use machine learning derived snow pattern maps to inform SnowModel, a physically-based snow process model. We develop and test this technique at the Teller 27 Seward Peninsula NGEE-Arctic study site. This dataset includes 5 *.nc files of model inputs and outputs plus one user guide (*.pdf). We present the data we use to drive SnowModel (vegetation type, digital elevation model), the snow pattern maps used to inform SnowModel (Standardized Depth Values maps, Machine Learning Snow Distribution Pattern), and our machine learning, SnowModel, and Hybrid snow depth and snow water equivalent results. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa

Abstract Background Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials. Results Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions. Conclusion Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structuring Nutrient Yields throughout Mississippi/Atchafalaya River Basin Using Machine Learning Approaches

To minimize the eutrophication pressure along the Gulf of Mexico or reduce the size of the hypoxic zone in the Gulf of Mexico, it is important to understand the underlying temporal and spatial variations and correlations in excess nutrient loads, which are strongly associated with the formation of hypoxia. This study’s objective was to reveal and visualize structures in high-dimensional datasets of nutrient yield distributions throughout the Mississippi/Atchafalaya River Basin (MARB). For this purpose, the annual mean nutrient concentrations were collected from thirty-three US Geological Survey (USGS) water stations scattered in the upper and lower MARB from 1996 to 2020. Eight surface water quality indicators were selected to make comparisons among water stations along the MARB over the past two decades. Principal component analysis (PCA) was used to comprehensively evaluate the nutrient yields across thirty-three USGS monitoring stations and identify the major contributing nutrient loads. The results showed that all samples could be analyzed using two main components, which accounted for 81.6% of the total variance. The PCA results showed that yields of orthophosphate (OP), silica (SI), nitrate–nitrites (NO 3 -NO 2 ), and total suspended sediment (TSS) are major contributors to nutrient yields. It also showed that land-planted crops, density of population, domestic and industrial discharges, and precipitation are fundamental causes of excess nutrient loads in MARB. These factors are of great significance for the excess nutrient load management and pollution control of the Mississippi River. It was found that the average nutrient yields were stable within the sub-MARB area, but the large nitrogen yields in the upper MARB and the large phosphorus yields in the lower MARB were of great concern. t-distributed stochastic neighbor embedding (t-SNE) revealed interesting nonlinear and local structures in nutrient yield distributions. Clustering analysis (CA) showed the detailed development of similarities in the nutrient yield distribution. Moreover, PCA, t-SNE, and CA showed consistent clustering results. This study demonstrated that the integration of dimension reduction techniques, PCA, and t-SNE with CA techniques in machine learning are effective tools for the visualization of the structures of the correlations in high-dimensional datasets of nutrient yields and provide a comprehensive understanding of the correlations in the distributions of nutrient loads across the MARB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Revealing Nanoscale Chemical Heterogeneities in Polycrystalline Mo-BiVO 4 Thin Films

The activity of polycrystalline thin film photoelectrodes is impacted by local variations of the material properties due to the exposure of different crystal facets and the presence of grain/domain boundaries. Here a multi-modal approach is applied to correlate nanoscale heterogeneities in chemical composition and electronic structure with nanoscale morphology in polycrystalline Mo-BiVO 4 . By using scanning transmission X-ray microscopy, the characteristic structure of polycrystalline film is used to disentangle the different X-ray absorption spectra corresponding to grain centers and grain boundaries. Comparing both spectra reveals phase segregation of V 2 O 5 at grain boundaries of Mo-BiVO 4 thin films, which is further supported by X-ray photoelectron spectroscopy and many-body density functional theory calculations. Theoretical calculations also enable to predict the X-ray absorption spectral fingerprint of polarons in Mo-BiVO 4 . After photo-electrochemical operation, the degraded Mo-BiVO 4 films show similar grain center and grain boundary spectra indicating V 2 O 5 dissolution in the course of the reaction. Overall, these findings provide valuable insights into the degradation mechanism and the impact of material heterogeneities on the material performance and stability of polycrystalline photoelectrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Finite element modeling of electropolishing of radio frequency dipole Nb crab cavity in hydrofluoric-sulfuric acid electrolyte

The superior performance of niobium superconducting radio frequency (SRF) cavities is achieved through meticulous surface treatment, notably via chemical electropolishing, ensuring exceptionally smooth surfaces. While this technique has been extensively employed for cylindrically symmetric structures, addressing more intricate geometries poses significant challenges in achieving uniform polishing and controlled material removal, especially when moving away from retractable cathodes. This presents a multifaceted electrochemical, thermal, and fluid dynamics conundrum. A prime example is the 197 MHz radio frequency dipole (RFD) crabbing cavity proposed for the Electron Ion Collider (EIC) project, exemplifying such complex structures. Our groundbreaking work integrates the localized oxide thickness variation, considering its impact on current distribution and Joule heating, within the framework of multi-physics modeling using the COMSOL® simulation suite. This approach was applied to a comprehensive model of the RFD cavity, allowing us to investigate optimal external cooling water flow conditions necessary for achieving desirable outcomes. In conclusion, this illustrates the increasing utility of such multi-physics codes to greatly aid in designing solutions to challenging processing requirements for increasingly complex accelerator cavities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

High precision control of laser energy for laser–matter interaction studies

Precise, highly reproducible control of the laser energy is required for high confidence laser–matter interaction research such as in dynamic compression science and high energy density physics. The energy must be adjustable without affecting the pulse shape (time varying intensity) or beam smoothness. We have developed a convenient two-stage energy tuning method for a nominal 100 J, 351 nm (UV) laser. The energy is adjusted in 10 J (10%) increments by operating the laser at full energy and inserting a beam splitter in the laser output. As the splitter is located after the final frequency tripling optics, the UV pulse shape is unchanged. The energy is varied by substituting a splitter of different reflectivity. For finer 3 J (3%) increments, the infrared pulse is attenuated inside the laser before the final amplifier. This requires modest tuning to preserve the pulse shape. The demonstrated variation in shot-to-shot reproducibility is less than +/–2.5 J (5% of the full energy), irrespective of the laser output energy. These approaches can be adapted to most ~100 J class lasers. We describe these techniques and show two examples where they have elucidated the underlying physics in laser shock compression experiments. One used only the beam splitters to establish the pressure for melting in iron. The other combined both techniques to finely increment the peak stress (~2 GPa steps) in germanium to precisely determine the onset and completion of melting—including the melting kinetics. Furthermore, these unambiguous results would not be possible without the developments described here.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Photoinduced electric effects in various plasmonic materials

Photoinduced voltages associated with surface plasmon polariton excitations are studied both theoretically and experimentally in various plasmonic systems as the function of material, wavelength, and type of structure. Experimental photovoltage normalized to the absorbed power shows a general decrease upon an increase in the wavelength, enhancement in the nanostructured samples, and a strong variation in the magnitude as a function of the material, which are not in line with the theoretical predictions of the simple plasmonic pressure approach. Furthermore, the results can be used for clarification of the mechanisms and further development of an adequate theoretcial approach to the plasmon drag effect.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Learning nuclear cross sections across the chart of nuclides with graph neural networks

We explore the use of deep learning techniques to learn how nuclear cross sections change as we add or remove protons and neutrons. As a proof of principle, we focus on the neutron-induced reactions in the fast energy regime. Our approach follows a two-stage learning framework. First, we apply representation learning to encode cross section data into a latent space using either variational autoencoders (VAEs) or implicit neural representations (INRs). Then, we train graph neural networks (GNNs) on the resulting embeddings to predict missing values across the nuclear chart by leveraging the topological structure of neighboring isotopes. We demonstrate accurate cross section predictions within a 9 × 9 block of missing nuclei. We also find that the optimal GNN training strategy depends on the type of latent representation used, with VAE embeddings performing best under end-to-end optimization in the original space, while INR embeddings achieve better results when the GNN is trained only in the latent space. Furthermore, using clustering algorithms, we map groups of latent vectors into regions of the nuclear chart and show that VAEs and INRs can discover some of the neutron magic numbers. These findings suggest that deep-learning models based on the representation encoding of cross sections combined with graph neural networks hold significant potential in augmenting nuclear theory models, e.g., by providing reliable estimates of covariances of cross sections, including cross-material covariances.

Machine learning↗
DOE OSTI2020

Using Machine Learning Tools to Predict Compressor Stall

Clean energy has become an increasingly important consideration in today’s power systems. As the push for clean energy continues, many coal-fired power plants are being decommissioned in favor of renewable power sources such as wind and solar. However, the intermittent nature of renewables means that dynamic load following traditional power systems is crucial to grid stability. With high flexibility and fast response at a wide range of operating conditions, gas turbine systems are poised to become the main load following component in the power grid. Yet, rapid changes in load can lead to fluid flow instabilities in gas turbine power systems. These instabilities often lead to compressor surge and stall, which are some of the most critical problems facing the safe and efficient operation of compressors in turbomachinery today. Although the topic of compressor surge and stall has been extensively researched, no methods for early prediction have been proven effective. This study explores the utilization of machine learning tools to predict compressor stall. The long short-term memory (LSTM) model, a form of recurrent neural network (RNN), was trained using real compressor stall datasets from a 100 kW recuperated gas turbine power system designed for hybrid configuration. Two variations of the LSTM model, classification and regression, were tested to determine optimal performance. The regression scheme was determined to be the most accurate approach, and a tool for predicting compressor stall was developed using this configuration. Overall, results show that the tool is capable of predicting stalls 5–20 ms before they occur. With a high-speed controller capable of 5 ms time-steps, mitigating action could be taken to prevent compressor stall before it occurs.

42 ENGINEERING↗
DOE OSTI2020

Coupling between plasmonic and photonic crystal modes in suspended three-dimensional meta-films

A complementary metal oxide semiconductor (CMOS) compatible fabrication method for creating three-dimensional (3D) meta-films is presented. In contrast to metasurfaces, meta-films possess structural variation throughout the thickness of the film and can possess a sub-wavelength scale structure in all three dimensions. Here we use this approach to create 2D arrays of cubic silicon nitride unit cells with plasmonic inclusions of elliptical metallic disks in horizontal and vertical orientations with lateral array-dimensions on the order of millimeters. Fourier transform infrared (FTIR) spectroscopy is used to measure the infrared transmission of meta-films with either horizontally or vertically oriented ellipses with varying eccentricity. Shape effects due to the ellipse eccentricity, as well as localized surface plasmon resonance (LSPR) effects due to the effective plasmonic wavelength are observed in the scattering response. The structures were modeled using rigorous coupled wave analysis (RCWA), finite difference time domain (Lumerical), and frequency domain finite element (COMSOL). The silicon nitride support structure possesses a complex in-plane photonic crystal slab band structure due to the periodicity of the unit cells. We show that adjustments to the physical dimensions of the ellipses can be used to control the coupling to this band structure. The horizontally oriented ellipses show narrow, distinct plasmonic resonances while the vertically oriented ellipses possess broader resonances, with lower overall transmission amplitude for a given ellipse geometry. We attribute this difference in resonance behavior to retardation effects. The ability to couple photonic slab modes with plasmonic inclusions enables a richer space of optical functionality for design of metamaterial-inspired optical components.

42 ENGINEERING↗
DOE OSTI2021

Gaussian Kernel Methods for Seismic Fragility and Risk Assessment of Mid-Rise Buildings

Seismic fragility functions can be evaluated using the cloud analysis method with linear regression which makes three fundamental assumptions about the relation between structural response and seismic intensity: log-linear median relationship, constant standard deviation, and Gaussian distributed errors. While cloud analysis with linear regression is a popular method, the degree to which these individual and compounded assumptions affect the fragility and the risk of mid-rise buildings needs to be systematically studied. This paper conducts such a study considering three building archetypes that make up a bulk of the building stock: RC moment frame, steel moment frame, and wood shear wall. Gaussian kernel methods are employed to capture the data-driven variations in the median structural response and standard deviation and the distributions of residuals with the intensity level. With reference to the Gaussian kernels approach, it is found that while the linear regression assumptions may not affect the fragility functions of lower damage states, this conclusion does not hold for the higher damage states (such as the Complete state). In addition, the effects of linear regression assumptions on the seismic risk are evaluated. For predicting the demand hazard, it is found that the linear regression assumptions can impact the computed risk for larger structural response values. However, for predicting the loss hazard with downtime as the decision variable, linear regression can be considered adequate for all practical purposes.

58 GEOSCIENCES↗
DOE OSTI2020

Forward Modeling the Orbits of Companions to Pulsating Stars from Their Light Travel Time Variations

Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to undersampling near periapsis. We introduce a new approach that simultaneously forward models each sample in the light curve and show that this method improves upon current sensitivity limits—especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the signal-to-noise ratio of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short-period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗