Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable selection”

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

Fe-rich X-ray amorphous material records past climate and persistence of water on Mars

X-ray amorphous material comprises 15-73 wt.% of sedimentary rocks and eolian sediments in Gale crater. This material is variably siliceous and iron rich but aluminum poor. The presence of volatiles is consistent with the existence of incipient weathering products. To better understand the implications of this material for past aqueous conditions on Mars, here we investigate X-ray amorphous material formation and longevity within terrestrial iron rich soils with varying ages and environmental conditions using bulk and selective dissolution methods, X-ray diffraction, and transmission electron microscopy. Results indicate that in situ aqueous alteration is required to concentrate iron into clay-size fraction material. Cooler climates promote the formation and persistence of X-ray amorphous material whereas warmer climates promote the formation of crystalline secondary phases. Iron rich X-ray amorphous material formation and persistence on Mars are therefore consistent with past cool and relatively wet environments followed by long-term cold and dry conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A device for volatile organic compound (VOC) analysis from skin using heated dynamic headspace sampling

Abstract Human skin is an important source of volatile organic compounds (VOCs) offering noninvasive methods to gain clinical metabolite information. This work was focused on the development of a skin sampling device based on a dynamic headspace sampling method with the addition of temperature to increase VOC metabolite recovery. The device preconcentrates skin VOC emissions onto a sorbent substrate, which can either be preserved for offline analysis or attached to a real time sensor downstream. In this work, skin VOC samples were analyzed offline using thermal desorption-gas chromatography-mass spectrometry. A list of 10 common skin VOCs was pre-selected to optimize parameters of sampling time, sampling temperature, and sorbent selection. Overall, this study highlights an effective skin VOC sampling technology with a heating dimension (40 °C, rather than 30 °C or no heating) with a sampling time of 15 min (rather than 5 or 30 mins) and onto Tenax TA sorbent (rather than PDMS), which collectively increases the recovery of compounds with lower vapor pressure and decreases the observed variability in skin VOC measurements. Finally, a list of 79 skin VOC compounds were detected and identified within a cohort of 20 young, healthy volunteers.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Cosmic topology. Part IIa. Eigenmodes, correlation matrices, and detectability of orientable Euclidean manifolds

Abstract If the Universe has non-trivial spatial topology, observables depend on both the parameters of the spatial manifold and the position and orientation of the observer. In infinite Euclidean space, most cosmological observables arise from the amplitudes of Fourier modes of primordial scalar curvature perturbations. Topological boundary conditions replace the full set of Fourier modes with specific linear combinations of selected Fourier modes as the eigenmodes of the scalar Laplacian. We present formulas for eigenmodes in orientable Euclidean manifolds with the topologiesE 1 –E 6 ,E 11 ,E 12 ,E 16 , andE 18 that encompass the full range of manifold parameters and observer positions, generalizing previous treatments. Under the assumption that the amplitudes of primordial scalar curvature eigenmodes are independent random variables, for each topology we obtain the correlation matrices of Fourier-mode amplitudes (of scalar fields linearly related to the scalar curvature) and the correlation matrices of spherical-harmonic coefficients of such fields sampled on a sphere, such as the temperature of the cosmic microwave background (CMB). We evaluate the detectability of these correlations given the cosmic variance of the observed CMB sky. We find that topologies where the distance to our nearest clone is less than about 1.2 times the diameter of the last scattering surface of the CMB give a correlation signal that is larger than cosmic variance noise in the CMB. This implies that if cosmic topology is the explanation of large-angle anomalies in the CMB, then the distance to our nearest clone is not much larger than the diameter of the last scattering surface. We argue that the topological information is likely to be better preserved in three-dimensional data, such as will eventually be available from large-scale structure surveys.

Astronomy & Astrophysics↗
DOE OSTI2023

Geospatial and Information Substitution and Anonymization Tool (GISA)

The Geospatial and Information Substitution and Anonymization Tool (GISA) incorporates techniques for obfuscating identifiable information from point data or documents, while simultaneously maintaining chosen variables to enable future use and meaningful analysis. This approach promotes collaboration and data sharing while also reducing the risk of exposure to sensitive information. GISA can be used in a number of different ways, including the anonymization of point spatial data, batch replacement/removal of user-specified terms from file names and from within file content, and aid with the selection and redaction of images and terms based on recommendations using natural language processing. Version 1 of the tool, published here, has updated functionality and enhanced capabilities to the beta version published in 2023. Please see User Documentation for further information on capabilities, as well as a guide for how to download and use the tool. If there are any feedback you would like to provide for the tool, please reach out with your feedback to edxsupport@netl.doe.gov. Disclaimer: This project was funded by the United States Department of Energy, National Energy Technology Laboratory, in part, through a site support contract. Neither the United States Government nor any agency thereof, nor any of their employees, nor the support contractor, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. The Geospatial and Information Substitution and Anonymization Tool (GISA) was developed jointly through the U.S. DOE Office of Fossil Energy and Carbon Management’s EDX4CCS Project, in part, from the Bipartisan Infrastructure Law.

Bipartisan Infrastructure Law↗
DOE OSTI2022

First Extraction of Single Differential Cross-Sections on Argon for CC1$\mu$2p0$\pi$ Event Topologies in the MicroBooNE Detector

The precise measurement of cross-sections for a variety of interactions is critical to the success of forthcoming flagship neutrino experiments. Of special interest are neutrino interactions that leave the nucleus in a 2-particle 2-hole state (2p2h). This note will present the first cross-section measurement for the production of 2p2h states on argon. Using MicroBooNE data collected from the first three years of operation, we select events corresponding to a charged-current $ν_µ$ interaction that left the argon nucleus in a 2p2h state. These interactions produce a topology with one muon and two protons in the final state (CC1$\mu$2p0$\pi$). This analysis targets a single differential cross-section measurement for CC1$\mu$2p0$\pi$ event topologies in a variety of kinematic variables. Comparisons are made to a set of theoretical models that explore different cross-section modeling configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Performance Assessment of a Photonic Radiative Cooling System for Office Buildings

Spectrally selective materials show high technical potential for use as radiative cooling heat exchangers for buildings. Recent work by Stanford University has applied a photonic approach to tailor the optical properties of a coating material. In a prototype test, a radiative cooler based on that photonic material has demonstrated the ability to maintain radiator surfaces at below-ambient temperatures in the presence of intense, direct sunlight. This paper presents the simulation and parametric analysis of photonic radiative cooling system that integrates the use of the photonic radiative cooler and radiant floor cooling via the whole energy simulation program EnergyPlus. The simulation was made for a medium office building with three floors and a total floor area of 5,000 m2. Three key design parameters including the storage tank size, the radiator area, and the number of floors are expected to have a large impact on system performance were investigated in this paper. The results show that the percentage of cooling electricity savings from the photonic radiative cooling system relative to a reference variable-air-volume system has a linear relationship with the radiator area, and a quadratic correlation with both tank size and the number of floors. The findings from this parametric analysis are valuable to understand how component sizing affects system performance.

Wang, Weimin↗
DOE OSTI2024

Impact of ultraviolet-B radiation on early-season morpho-physiological traits of indica and japonica rice genotypes

Ultraviolet (UV)-B radiation is considered one of the major detrimental rays coming from the Sun. UV-B radiation has a harmful impact on plant growth and development. The effect of UV-B radiation was studied on 64 rice (Oryza sativa L.) genotypes during the vegetative season. An equal number of genotypes from the japonica (50%) and indica (50%) subspecies were phenotyped using the Soil-Plant-Atmosphere-Research (SPAR) units. The 10 kJ UV-B was imposed 12 days after planting (DAP) and continued for three weeks (21 d). Based on the combined ultraviolet-B radiation response index (CUVBRI) for each genotype, the 64 rice genotypes were classified into sensitive, moderately sensitive, moderately tolerant, and tolerant. Various shoot traits, such as plant height, tiller, and leaf numbers, were measured. We also studied critical root phenological traits like root volume, diameter, tips, and forks. Out of all the studied shoot traits, leaf area showed maximum reduction for both indica (54%) and japonica (48%). Among the root traits, root length decreased by negligible (1%) for indica as compared to japonica (5%), while root crossing and forks showed a maximum decline for japonica (37 and 42%), respectively. This study is timely, meaningful, and required because it will help breeders select a tolerant or sensitive rice line for better yield and production under abiotic stresses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Smart culture medium optimization for recombinant protein production: Experimental, modeling, and AI/ML-driven strategies

Recombinant protein production (RPP) is central to biotechnology, where recombinant proteins are used as either end products or catalysts in the synthesis of chemicals, fuels, and materials. Among the major cost drivers, culture medium plays a pivotal role in determining protein yield and quality. This review presents a comprehensive perspective on the critical stages of “smart” culture medium optimization: planning, screening, modeling, optimization, and validation. In the planning stage, we examine the nutritional and energetic roles of medium components, including carbon, nitrogen, amino acids, salts, and trace metals, and their impacts on culture parameters such as pH, oxidative state, and osmolality. We highlight the variability in trace metal content due to water sources, culture vessels, and raw materials, which can substantially influence RPP. The screening stage covers Design of Experiments (DoE) approaches, assessing their theoretical basis, implementation, and limitations. For modeling, we describe methods that integrate experimental data to develop predictive models for smart medium formulation. Model-based optimization strategies can then be employed to select optimal media compositions for a given application. The validation stage aims to evaluate model predictions and provide feedback for model training and refinement. Finally, we survey mechanistic and artificial intelligence/machine learning (AI/ML)-driven models as integrated, transformational tools for predictive modeling of bioprocess conditions, nutrient availability, cellular metabolism, and protein quality, with the goal of optimizing culture media to enhance protein yields while reducing costs and environmental impact. We conclude by addressing the challenges of translating laboratory-scale medium optimization to industrial-scale settings and exploring future AI/ML-driven approaches that may overcome current bottlenecks and accelerate medium design for RPP. Overall, this review provides a unified framework for advancing smart medium design in RPP.

Artificial Intelligence/Machine Learning (AI/ML)↗
DOE OSTI2021

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Results of the International Energy Agency Bioenergy Round Robin on the Analysis of Heteroatoms in Biomass Liquefaction Oils

A round robin study evaluating the analysis of biomass liquefaction oils (BLOs) from fast pyrolysis and hydrothermal liquefaction (HTL) was performed, involving fifteen laboratories in seven countries in order to assess the current status of analytical techniques for the determination of nitrogen, sulfur, and chlorine content in BLOs and evaluate potential differences in origin (i.e. fast pyrolysis vs HTL). The BLOs were produced from a range of feedstocks including pine, mixed softwoods, forest residues, micro-algae, miscanthus, and wheat straw to cover a variety in nitrogen, sulfur, and chlorine content and speciation. Nine samples were distributed, comprised of eight separate BLOs and one blind duplicate produced by five producers. The samples were analyzed for water, carbon, hydrogen, nitrogen, sulfur, and chlorine content. No analytical test method was mandated; laboratories were encouraged to utilize whichever method they determined would be most applicable, relying on the existing body of BLO literature as a guide. The results of this round robin study are presented in this paper. The results of the carbon, hydrogen, and water measurements as reference analyses had relative standard deviations (2.9%, 3.5, and 5.6%, respectively) that were comparable to those found in past round robin studies on fast pyrolysis bio-oil. The analysis of nitrogen, sulfur, and chlorine showed higher levels of variability. Laboratories mostly chose the same method for water, carbon, hydrogen and nitrogen determination whereas there were a variety of methods chosen for sulfur and chlorine determination. The results suggest that specific analytical methods for the determination of nitrogen, sulfur, and chlorine should be further refined to ensure reproducible and accurate results for BLO analysis due to their importance on emissions, material selection, and catalyst activity.

Biomass, IEA Bioenergy, pyrolysis, hydrothermal li↗
DOE OSTI2017

Multireference Density Functional Theory with Generalized Auxiliary Systems for Ground and Excited States

To describe static correlation, we create a new approach to density functional theory (DFT), which uses a generalized auxiliary system that is of a different symmetry, such as particle number or spin, from that of the physical system. The total energy of the physical system consists of two parts: the energy of the auxiliary system, which is determined with a chosen density functional approximation (DFA), and the excitation energy from an approximate linear response theory that restores the symmetry to that of the physical system, thus rigorously leading to a multideterminant description of the physical system. The electron density of the physical system is different from that of the auxiliary system and is uniquely determined from the functional derivative of the total energy with respect to the external potential. Our energy functional is thus an implicit functional of the physical system density, but an explicit functional of the auxiliary system density. We demonstrate that the total energy minimum and stationary states, describing the ground and excited states of the physical system, can be obtained by a self-consistent optimization with respect to the explicit variable, the generalized Kohn–Sham noninteracting density matrix. We have developed the generalized optimized effective potential method for the self-consistent optimization. Among options of the auxiliary system and the associated linear response theory, reformulated versions of the particle–particle random phase approximation (pp-RPA) and the spin-flip time-dependent density functional theory (SF-TDDFT) are selected for illustration of principle. Numerical findings show that our multireference DFT successfully describes static correlation in bond dissociation and double bond rotation.

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

Exact and efficient hybrid Monte Carlo algorithm for accelerated Bayesian inference of gene expression models from snapshots of single-cell transcripts

Single cells exhibit a significant amount of variability in transcript levels, which arises from slow, stochastic transitions between gene expression states. Elucidating the nature of these states and understanding how transition rates are affected by different regulatory mechanisms require state-of-the-art methods to infer underlying models of gene expression from single cell data. A Bayesian approach to statistical inference is the most suitable method for model selection and uncertainty quantification of kinetic parameters using small data sets. However, this approach is impractical because current algorithms are too slow to handle typical models of gene expression. To solve this problem, we first show that time-dependent mRNA distributions of discrete-state models of gene expression are dynamic Poisson mixtures, whose mixing kernels are characterized by a piecewise deterministic Markov process. Here, we combined this analytical result with a kinetic Monte Carlo algorithm to create a hybrid numerical method that accelerates the calculation of time-dependent mRNA distributions by 1000-fold compared to current methods. We then integrated the hybrid algorithm into an existing Monte Carlo sampler to estimate the Bayesian posterior distribution of many different, competing models in a reasonable amount of time. We demonstrate that kinetic parameters can be reasonably constrained for modestly sampled data sets if the model is known a priori. If there are many competing models, Bayesian evidence can rigorously quantify the likelihood of a model relative to other models from the data. We demonstrate that Bayesian evidence selects the true model and outperforms approximate metrics typically used for model selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Signal-Compensation-Based Adaptive PID Control for Fused Magnesia Smelting Processes

For fused magnesia smelting processes, the proportional-integral-derivative (PID) controller is difficult to regulate the smelting currents within desired ranges due to unexpected variations of the dynamics caused by material feeding operations. To solve this problem, we first proposed its dynamic model composed of a linear model and an uncertain nonlinear dynamic system that describes the unknown time-varying dynamics of nonlinearity and strong couplings among controlled variables. Then, using the accurately calculated value of the uncertain part at the previous sampling time instant and the tracking error that reflects the influence of its change rate, we developed a feedforward and a one-step optimal control to obtain two signal compensators to eliminate the impact of frequent variations of the uncertainties. For the compensators, we proposed a weight selection method to minimize the sum of the squared tracking error with ensured closed-loop stability. Based on the abovementioned modeling and control structure, we finally developed a novel adaptive PID controller augmented with signal compensators to control the whole process. Industrial applications have shown that the proposed method can always keep the tracking error and the control input within their targeted ranges, leading to a significantly reduced energy consumption per ton and increased fused magnesia production.

Wang, Weizhou↗
DOE OSTI2023

Performance assessment for climate intervention (PACI): preliminary application to a stratospheric aerosol injection scenario

As the prospect of exceeding global temperature targets set forth in the Paris Agreement becomes more likely, methods of climate intervention are increasingly being explored. With this increased interest there is a need for an assessment process to understand the range of impacts across different scenarios against a set of performance goals in order to support policy decisions. The methodology and tools developed for Performance Assessment (PA) for nuclear waste repositories shares many similarities with the needs and requirements for a framework for climate intervention. Using PA, we outline and test an evaluation framework for climate intervention, called Performance Assessment for Climate Intervention (PACI) with a focus on Stratospheric Aerosol Injection (SAI). We define a set of key technical components for the example PACI framework which include identifying performance goals, the extent of the system, and identifying which features, events, and processes are relevant and impactful to calculating model output for the system given the performance goals. Having identified a set of performance goals, the performance of the system, including uncertainty, can then be evaluated against these goals. Using the Geoengineering Large Ensemble (GLENS) scenario, we develop a set of performance goals for monthly temperature, precipitation, drought index, soil water, solar flux, and surface runoff. The assessment assumes that targets may be framed in the context of risk-risk via a risk ratio, or the ratio of the risk of exceeding the performance goal for the SAI scenario against the risk of exceeding the performance goal for the emissions scenario. From regional responses, across multiple climate variables, it is then possible to assess which pathway carries lower risk relative to the goals. The assessment is not comprehensive but rather a demonstration of the evaluation of an SAI scenario. Future work is needed to develop a more complete assessment that would provide additional simulations to cover parametric and aleatory uncertainty and enable a deeper understanding of impacts, informed scenario selection, and allow further refinements to the approach.

climate intervention↗
DOE OSTI2019

The Effect of Calcium on the Cohesive Strength and Flexural Properties of Low-Methoxyl Pectin Biopolymers

Pectin binds the mesothelial glycocalyx of visceral organs, suggesting its potential role as a mesothelial sealant. To assess the mechanical properties of pectin films, we compared pectin films with a less than 50% degree of methyl esterification (low-methoxyl pectin, LMP) to films with greater than 50% methyl esterification (high-methoxyl pectin, HMP). LMP and HMP polymers were prepared by step-wise dissolution and high-shear mixing. Both LMP and HMP films demonstrated a comparable clear appearance. Fracture mechanics demonstrated that the LMP films had a lower burst strength than HMP films at a variety of calcium concentrations and hydration states. The water content also influenced the extensibility of the LMP films with increased extensibility (probe distance) with an increasing water content. Similar to the burst strength, the extensibility of the LMP films was less than that of HMP films. Flexural properties, demonstrated with the 3-point bend test, showed that the force required to displace the LMP films increased with an increased calcium concentration ( p < 0.01). Toughness, here reflecting deformability (ductility), was variable, but increased with an increased calcium concentration. Similarly, titrations of calcium concentrations demonstrated LMP films with a decreased cohesive strength and increased stiffness. We conclude that LMP films, particularly with the addition of calcium up to 10 mM concentrations, demonstrate lower strength and toughness than comparable HMP films. These physical properties suggest that HMP has superior physical properties to LMP for selected biomedical applications.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Hydrogen Leak Modeling for Development of Smart Distributed Monitoring Under Unintended Releases

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower. Hydrogen has the potential to play a crucial role in decarbonizing industrial processes that are currently reliant on fossil fuels and provide long-duration and/or seasonal energy storage to enable electricity decarbonization. Hydrogen can also be used as a fuel for fuel cell vehicles, providing a zero-emission alternative to traditional internal combustion engines. DOE launched the Hydrogen Energy Earthshot (Hydrogen Shot) in June 2021 to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). While promising, Hydrogen is highly-flammable, and in the presence of oxygen, it can form explosive mixtures. . Therefore, understanding leak scenarios is essential to evaluate and mitigate the safety risks associated with potential hydrogen leaks. An increased understanding of leak behavior, and having tools to model leaks, can help assess how hydrogen would disperse in different environments, influencing emergency response plans and safety measures, and identify potential issues with materials and design systems that can withstand the challenges posed by hydrogen. Recently, researchers have attempted to study hydrogen leaks for development of risk management strategies. However, the focus has been on closed or semi-closed spaces like storage rooms, vehicles, garages, and fueling stations - all promising locations for future hydrogen infrastructure. In this presentation, the modeling environment extends the span of research further by modeling hydrogen leak in an outdoor, open space. We will present the key challenges with modeling hydrogen leaks in an uncontrollable environment, how they were handled, and how modeling results informed sensor selection and placement. A Hydrogen research facility at the National Renewable Energy Laboratory (NREL) was used as a case study to model hydrogen leaks. In the future, Hydrogen wide area detection methodologies will be developed and tested at this site to monitor for unintended and operational hydrogen releases. The data generated from modeling will be used to develop a predictive model to detect hydrogen leak location based on concentration measured by sensors in this open space. Furthermore, the facility was also chosen because controlled hydrogen releases can be performed. A computational fluid dynamics (CFD) based modeling approach was taken to model hydrogen leak. The full-scale hydrogen facility was modeled with a large ambient domain. The electrolyzer at the facility can produce a controlled release rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind direction, wind speed at various altitudes, and temperature were used as inputs to the model. To capture the variability of weather conditions, a subset of the weather conditions experienced during daytime hours without precipitation over the course of three months was generated; using established data clustering techniques, a total of 100 condition sets were chosen. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. These distributions are compared to experimental data from a constant mass flow, controlled hydrogen release at the facility. The stochastic wind conditions of the release make direct validation difficult, therefore, statistical comparison approaches were used. Wind conditions are found to significantly impact the release behavior, including direction and concentration. Sensor selection and placement is proposed for the facility and is now based on release behavior predicted for the facility given its weather patterns; this is much more informed than without the modeling results. The methodology and analysis procedure can be translated to other facilities using modified geometries and site-specific weather conditions. Hydrogen holds great promise as a renewable energy fuel, but ensuring safety in its production, storage, and use is paramount. Studying potential leak scenarios in an open space will help develop sensors to detect hydrogen on a large spectrum of concentration and eventually build a smart distributed monitoring system.

CFD↗
DOE OSTI2019

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Numerical Simulation of Reservoir Response to Depressurization

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Boswell, Ray↗