Search NASA⌕ Search

SEARCH · Search NASA

Results for “solvation structure”

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.

213 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

The Transition from Hydrogen Bonding to Ionization in (HCI)n(NH3)n and (HCI)n(H2O)n Clusters: Consequences for Anharmonic Vibrational Spectroscopy

Anharmonic vibrational frequencies and intensities are calculated for 1:1 and 2:2 (HCl)(sub n)(NH3)(sub n) and (HCl)(sub n)(H2O)(sub n) complexes, employing the correlation-corrected vibrational self-consistent field method with ab initio potential surfaces at the MP2/TZP computational level. In this method, the anharmonic coupling between all vibrational modes is included, which is found to be important for the systems studied. For the 4:4 (HCl)(sub n)(H2O)(sub n) complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Just as the (HCl)(sub n)(NH3)(sub n) structure switches from hydrogen-bonded to ionic for n=2, the (HCl)(sub n)(H2O)(sub n) switches to ionic structure for n=4. For (HCl)2(H2O)2, the lowest energy structure corresponds to the hydrogen-bonded form. However, configurations of the ionic form are separated from this minimum by a barrier of less than an O-H stretching quantum. This suggests the possibility of experiments on ionization dynamics using infrared excitation of the hydrogen-bonded form. The strong cooperative effects on the hydrogen bonding, and concomitant transition to ionic bonding, makes an accurate estimate of the large anharmonicity crucial for understanding the infrared spectra of these systems. The anharmonicity is typically of the order of several hundred wave numbers for the proton stretching motions involved in hydrogen or ionic bonding, and can also be quite large for the intramolecular modes. In addition, the large cooperative effects in the 2:2 and higher order (HCl(sub n)(H2O)(sub n) complexes may have interesting implications for solvation of hydrogen halides at ice surfaces.

Chaban, Galina M.↗
DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2024

A Fluorinated Lewis Acidic Organoboron Tunes Polysulfide Complex Structure for High–Performance Lithium–Sulfur Batteries

Many challenges in lithium-sulfur (Li–S) batteries are associated with the radical change in lithium polysulfide (LPS) solubility during cycling, but chemical approaches to address such inconsistency are still lacking. Here, the use of a strong Lewis acidic fluorinated organoboron, tri(2,2,2-trifluoroethyl) borate (TFEB), is reported as a multi-functional mediator to simultaneously overcome multiple technical barriers in practical Li–S batteries. TFEB acts as an anion acceptor and forms strong molecular complexes with Lewis basic LPS. The TFEB-LPS complexes have consistent solubility across the full polysulfide spectrum and deliver several times improved better redox kinetics, unlocking a true redox catalytic mechanism that covers the majority of redox events in thick sulfur cathodes. As a result, Li–S batteries evaluated under practical conditions exhibit significantly improved discharge capacity, rate capability, and cycling stability with the addition of the TFEB additive. More importantly, TFEB also contributes to the stabilization of lithium anode in the presence of polysulfides by generating strong interfacial film. These attributes significantly improve the cycling stability of practical Li–S pouch cells, which are assembled with a unit energy density of 219 Wh kg –1 . Finally, the results provide new molecular insights on the design of unlocking solvation networks of practical Li–S systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

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

In-Situ FTIR Detection of Transition Metal (TM)-Ion Dissolution From Cathodes in Li-Ion Batteries

Transition metal (TM) ions, commonly Ni and Mn, play a crucial role in Li-ion battery cathodes as the reaction centers for rapid redox reactions. A major challenge with TM-based cathodes is capacity degradation, particularly at higher operating voltages. This degradation is closely linked to the dissolution of TMs from the cathode materials and their subsequent deposition on the anode. This process not only modifies the surface structure of the cathode but, more significantly, alters the SEI composition on the anode [1-2]. The dissolution of TMs cations into a liquid electrolyte from cathode materials, such as Mn-ion dissolution from Mn-rich cathode (LMR), is detrimental to the cycling performance of Li-ion batteries [3-4]. Much attention has been paid to this issue but there remains a lack of characterization techniques which can detect the TM-ion dissolution from the cathode during electrochemical measurements. In our study, we use in-situ ATR-FTIR as an effective technique to probe the TM-ion dissolution from the cathode. We have first demonstrated the detrimental effects of TM ions on the electrochemical performance of Li-ion batteries by adding a small amount of TM salt (50 mM Mn(PF6)) to the electrolyte of a Li-ion coin cell with LFP and graphite electrode. We observed a rapid capacity fade after the first delithiation cycle. To investigate TM ion dissolution, we established a baseline IR spectrum for various TM solvation states (such as Mn and Ni) by measuring concentration-dependent IR spectra. This baseline spectrum helps us detect TM ion dissolution during battery cycling. In this work, we discuss in detail the effect of TM ions on the electrochemical performance of Li-ion batteries and the detection of TM ions during battery cycling using in-situ FTIR spectroscopy. We will compare TM dissolution between coated and uncoated cathodes to examine the effect of cathode coatings to mitigate degradation due to TM dissolution and cross-over from cathode to anode. References: (1) Zhan, C.; Wu, T.; Lu, J.; Amine, K. Dissolution, migration, anddeposition of transition metal ions in Li-ion batteries exemplified byMn-based cathodes - a critical review. Energy Environ. Sci. 2018, 11,243-257. (2) Jung, R.; Linsenmann, F.; Thomas, R.; Wandt, J.; Solchenbach,S.; Maglia, F.; Stinner, C.; Tromp, M.; Gasteiger, H. A. Nickel,Manganese, and Cobalt Dissolution from Ni-Rich NMC and TheirEffects on NMC622-Graphite Cells. J. Electrochem. Soc. 2019, 166,A378-A389. (3) Zhao, L.; Chenard, E.; Capraz, O. O.; Sottos, N. R.; White, S.R. Direct Detection of Manganese Ions in Organic Electrolyte by UV-Vis Spectroscopy. J. Electrochem. Soc. 2018, 165, A345-A348 (4) Zhang, Y.; Hu, A.; Xia, D.; Hwang, S.; Sainio, S.; Nordlund, D.;Michel, F. M.; Moore, R. B.; Li, L.; Lin, F. Operando characterization and regulation of metal dissolution and redeposition dynamics nearbattery electrode surface. Nat. Nanotechnol. 2023, 18, 790.

25 ENERGY STORAGE↗
DOE OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
DOE OSTI2025

Coordination Chemistry and Photoluminescence of Sm(II) Dibenzo-24-crown-8 Complexes

Three Sm(II) dibenzo-24-crown-8 (db24c8) complexes were synthesized in anhydrous, air-free conditions via the reaction of SmI 2 with db24c8 and tetrabutylammonium tetraphenylborate ([TBA][BPh 4 ]; where needed) in acetonitrile (CH 3 CN), dimethoxyethane (DME), and tetrahydrofuran (THF) to yield [Sm(db24c8)(CH 3 CN) 2 ][BPh 4 ][I]·CH 3 CN, [Sm(db24c8)(DME)]I 2 , and [Sm(db24c8)(THF) 2 ]I 2 , respectively. In each case, a 10-coordinate, staggered dodecahedral (2:6:2) environment is formed around the Sm 2+ center that is completed by either two solvent molecules (CH 3 CN or THF) or one bidentate solvent molecule (DME). Inner-sphere solvent molecules can be excluded by reacting SmI 2 with db24c8 in 1:3 THF:toluene to yield Sm(db24c8)I 2 . Here, this molecule features a distorted, eight-coordinate, hexagonal pyramidal Sm 2+ metal center, where the coordinated db24c8 molecule shows a torsion angle unexpectedly close to the 180° antiperiplanar arrangement and two uncoordinated db24c8 oxygen atoms. Solution UV–vis–NIR measurements demonstrate that Sm 2+ is a good size match for the cavity of various db24c8 conformations and that Eu 2+ and Yb 2+ exhibit competition between acetonitrile solvation and the Eu 2+ and Yb 2+ /db24c8 complexes in solution. During excitation by 546 nm light, both [Sm(db24c8)(DME)]I 2 and [Sm(db24c8)(THF) 2 ]I 2 exhibit mixed 5d → 4f and 4f → 4f emission at 20 °C and exclusively 4f → 4f at −180 °C, whereas Sm(db24c8)I 2 only shows 5d → 4f emission regardless of temperature. Photoluminescence from [Sm(db24c8)(CH 3 CN) 2 ][BPh 4 ][I]·CH 3 CN is quenched.

Cations↗
DOE OSTI2025

Environmental photosensitizers can exhibit enhanced actinic absorption in microhydrated clusters compared to solution

Brown carbon chromophores at environmental air-water interfaces often act as photosensitizers that absorb sunlight and subsequently transfer energy to nearby molecules, initiating a wide variety of chemical reactions. Despite their importance to understanding daytime chemistry at these airwater interfaces, little is known about the role of the solvation environment on the photophysical properties of these photosensitizers. In this work, we present a joint experimental-theoretical study of the vibrational and photophysical properties of microhydrated protonated and deprotonated 4-benzoylbenzoic acid (4-BBA), a key model system for environmental photosensitizers. We find that for protonated 4-BBAH + ∙(H 2 O) 0-1 , representing photosensitizers in very acidic conditions, a single bright state dominates the UV-Vis spectrum between 280 and 400 nm. Comparing the experimental UV-Vis spectra and quantum chemistry-predicted spectra of 4-BBA + ∙(H 2 O) 0-2 , we find that the degree of microhydration has little effect on the UV-Vis spectra or the orbitals of the dominant feature. For deprotonated 4-BBA ‒ , representing photosensitizers in basic conditions, quantum chemistry calculations predict that the UV-Vis spectra are ~3x weaker in intensity than the brightest 4-BBAH + ∙(H 2 O) 0-1 features and were not observed experimentally. Quantum chemistry calculations predict one dominant UV-Vis feature is present in the spectra of 4-BBA ‒ ∙(H 2 O) 0-2 , which exhibit minor shifts with degree of microhydration. While 4-BBA in bulk solution over a range of pH values has relatively weak absorption within the solar actinic region, we show that microhydrated 4-BBA has bright transitions within the actinic region. This indicates that the complex structure of environmental air-water interfaces can shift the absorption maximum of photosensitizers into the actinic region for enhanced absorption of sunlight and subsequent enhancement of photosensitizer-driven reactions.

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

Exploring Local Vibrational Structure in Protein-Bound Chlorophyll a: Isotope-Enrichment Experiments and Electrostatic Analysis

Local chlorophyll (Chl) vibrations play an essential role in biological photosynthesis by facilitating rapid energy transfer between pigments. In research studies, they also provide a useful spectroscopic probe of the local protein environment that surrounds each pigment. However, measuring the complete vibrational absorption spectrum of a protein-bound Chl molecule is much more difficult than, for example, Chl in neat solvent due to overlap with protein vibrations that typically drown out Chl vibrational features. Resonance Raman and fluorescence spectroscopies provide a way around this problem for Franck–Condon active vibrations, but these often rely on cryogenic measurement conditions and fail to capture vibrational signatures from, e.g., ester-group vibrations that lack coupling to a convenient electronic transition. In the present contribution, we use 13 C-enrichment of the protein backbone to shift protein background signals to lower frequency, providing a largely clean spectral window in which to study local Chl a C═O stretch modes. A room-temperature absorption spectrum for Chl a in the water-soluble chlorophyll protein (WSCP) of Lepidium virginicum is thus extracted as a difference between protein-plus-pigment and protein-only vibrational spectra. Excellent agreement in the molecular fingerprint region with the vibrational spectrum of Chl a in organic solvents confirms that the resulting spectrum represents the response of the protein-bound Chl a molecule. Furthermore, the ester group resonance is observed to shift in response to the S53P mutation that eliminates a 17 3 ester group hydrogen bond from the protein environment. Finally, we analyze these experimental results using MD-based electrostatic analysis, finding that electric-field mapping at the C atom of the ester group provides a satisfactory explanation for the observed frequency shifts between organic solvent and protein environment. MD analysis further suggests that a red-shifted ester peak observed experimentally for the S53P mutant results from solvation of the 17 3 Chl ester group due increase water penetration into the Chl-binding pocket relative to wild-type WSCP. Here, we anticipate that these results will prove useful both for benchmarking future simulation work and as a reference for interpreting Chl vibrational spectra as a probe of pigment–protein interactions.

14 SOLAR ENERGY↗
DOE OSTI2025

Hydration and Restructuring of Polar Polymer Interfaces: Implications in Antifouling and Responsive Materials

Manipulating polymer interfaces is crucial for understanding how structure influences function in applications spanning biofouling prevention to energy storage. Moreover, observing how polymers adapt their microscopic structure to changes in their local environment can reveal essential properties that govern their performance in such applications, providing key insights into how to design more effective interfaces. Here, in this study, a series of “grafting-from” polymer brushes with side chains varying in charge, are probed by sum frequency generation (SFG) and modeled using all-atom molecular dynamics (MD) simulations to elucidate how chemical makeup and charge mediate interfacial restructuring in dry versus hydrated states. Results show that charge, in progressing from nonpolar to cationic to zwitterionic, results in dramatic changes in interfacial structure and overall hydration. While net neutral systems, regardless of bulk phase polarity, show minimal interfacial water structuring, the cationic species exhibits strong bulk water signals from the surface potential. Meanwhile, the polymer brushes themselves restructure in water differently independent of charge, impacting the functional groups that are presented to the aqueous phase. Nonpolar and cationic species for instance undergo a change in alkyl group orientations to accommodate hydrating water molecules, whereas the zwitterionic polymer becomes completely disordered in water. Overall, the structure-based behavior trends presented herein have implications in antifouling applications and responsive material interfaces.

biointerface↗
DOE OSTI2025

Humidity Disrupts Structural and Chiroptical Properties of Chiral 2D Perovskites

Chiral two-dimensional (2D) hybrid organic–inorganic metal halide perovskite semiconductors have emerged as an exceptional material platform with many design opportunities for spintronic applications. However, a comprehensive understanding of changes to the crystal structure and chiroptical properties upon exposure to atmospheric humidity has not been established. We demonstrate phase degradation to the 1D (MBA)PbI 3 (MBA = methylbenzylammonium) and the hypothetical (MBA) 3 PbI 5 ·H 2 O hydrate phases, accompanied by a reduction and disappearance of the chiroptical response. Firstprinciple simulations show that water molecules preferentially locate at the interface between the organic cations and the inorganic framework, thereby disrupting the hydrogen bonding, impacting both the structural chirality and stability of the material. These findings provide critical insights into phase degradation mechanisms and their impact on chiroptical activity in chiral 2D perovskites.

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

Pathways to High-Performance Salt Hydrate Thermochemical Energy Storage Materials and Systems

Thermochemical materials (TCMs) based on salt hydrates are promising for thermal energy storage as they combine high energy densities with low reaction temperatures. However, their adoption is hindered by poor structural integrity and degradation under hygrothermal cycling. Storage performance is governed not only by the chemical reaction, but also by the coupled thermo-chemo-mechanical behavior that evolves with cycling. Understanding and controlling this coupling across length scales (material-to-reactor) is necessary to improve TCM stability and lifetime. In this perspective, we discuss the shortcomings of current characterization approaches and emphasize the need for measuring transport properties and structural transformations using in situ techniques that capture the dynamic evolution of these materials. We also outline opportunities for multiscale modeling frameworks that link thermodynamics and mechanics, enabling predictive evaluation of composite architectures designed for cycling stability. We conclude by identifying research questions that must be addressed to transform TCMs into viable energy storage technologies.

Composites↗
NASA NTRS2016

Computational and Experimental Study of Li-doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Psi). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Psi. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Psi-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Psi and suggest a source of impedance for Li+ to reach the electrode.

Haskins, Justin B.↗
DOE OSTI2025

Encapsulation of Ba in InSb Framework Introduces Chirality in Clathrate-Like BaIn 4 Sb 4

The discovery of Zintl compounds remains a powerful strategy for identifying materials with tunable electronic and thermal transport properties. During a concerted search for new inorganic clathrates with In–Sb frameworks, we discovered BaIn 4 Sb 4 . The composition of this phase deviates from that expected for a type-I clathrate with tetrahedral coordination of all In and Sb atoms (Ba 8 In 31 Sb 15 ). Instead, in the chiral structure of BaIn 4 Sb 4 (space group P 3 1 21, No. 152), a part of the In atoms have a trigonal planar coordination of 1In + 2Sb, forming Sb 2 –In–In–Sb 2 nonplanar fragments isostructural to diborane(4) B 2 H 4 with D 2 d symmetry. Ba atoms are located inside 16-vertex In 8 Sb 8 polyhedra, which share vertices and edges to form a chiral framework around the 3 1 screw axes. The title compound is electron-balanced, [Ba 2+ ][In 2+ ] 2 [In 3+ ] 2 [Sb3–] 4 , which was confirmed by characterization of the charge and heat transport properties. BaIn 4 Sb 4 exhibits a low thermal conductivity and high Seebeck coefficient, suggesting its untapped potential for thermoelectric applications. Finally, density functional theory (DFT) calculations indicate that chemical doping may enhance carrier concentration and improve the originally low electrical conductivity, thus enhancing thermoelectric performance.

electrical properties↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗