Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLOPE”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigation of Interlayer Dielectric in BaTiO 3 /III–Nitride Transistors

In this paper, the impact of varying the thickness of the Al 2 O 3 interlayer dielectric on the electrical characteristics of BaTiO 3 /III-nitride transistors is investigated. In the findings, it is revealed that a minimum thickness of 8 nm for the Al 2 O 3 layer is crucial to maintain high device performance and protect against sputtering-induced damage during BaTiO 3 deposition. The fabricated BaTiO 3 /Al 2 O 3 /AlGaN/GaN high electron mobility transistors exhibit exceptional electrical properties, including a maximum current density of 700 mA mm –1 , an on-resistance of 5 Ω mm, an I ON /I OFF ratio of 10 7 , a subthreshold slope of 119 mV dec –1 , and significantly reduced gate leakage current. The devices with the optimal 8 nm Al 2 O 3 thickness demonstrates excellent agreement between theoretical and experimental values for effective mobility, achieving a value of 1188 cm 2 V –1 ·s at a 2D electron gas density of 10 13 cm –2 . Furthermore, in the study, it is confirmed that increasing the Al 2 O 3 thickness also improves the quality of interface charge density, as evidenced by the results obtained from capacitance–voltage measurements. In these findings, the critical role of controlling the Al 2 O 3 thickness in optimizing the electrical characteristics and overall performance of BaTiO 3 /III-nitride transistors are highlighted.

42 ENGINEERING↗
DOE OSTI2023

Watersheds and stream networks viewed longitudinally: Example insights from novel spatial portrayals of watershed characteristics

Abstract Longitudinal depictions of watershed structure and characteristics, including topography, stream networks, wetlands, ground water levels, and land use, can provide watershed knowledge and understanding unavailable from standard plan view maps. Three case studies provide examples of knowledge gained by applying longitudinal views of stream networks, watershed hydrologic behavior, and land use distributions. Longitudinal views of mountain stream networks show extreme variability in the slope‐area relationships of low Strahler order streams, large discontinuities in drainage area (large parts of drainage area space are absent in networks), and large variations in network curvature. Longitudinal views of a groundwater‐dominated headwater watershed increase the inference available from limited groundwater observations and clearly reveal how groundwater connections affect the permanence of surface water features and the distribution of vadose zone storage in the landscape. Plotting land uses longitudinally illuminates and allows a quantitative analysis of how land uses are distributed relative to topographic position. Viewing watersheds and stream networks longitudinally can provide new insights into watershed forms and processes and motivate new questions and research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Non-Equilibrium Synthesis of Highly Active Nanostructured, Oxygen-Incorporated Amorphous Molybdenum Sulfide HER Electrocatalyst

Molybdenum sulfide emerged as promising hydrogen evolution reaction (HER) electrocatalyst thanks to its high intrinsic activity, however its limited active sites exposure and low conductivity hamper its performance. To address these drawbacks, the non-equilibrium nature of pulsed laser deposition (PLD) is exploited to synthesize self-supported hierarchical nanoarchitectures by gas phase nucleation and sequential attachment of defective molybdenum sulfide clusters. Furthermore, the physics of the process are studied by in situ diagnostics and correlated to the properties of the resulting electrocatalyst. The as-synthesized architectures have a disordered nanocrystalline structure, with nanodomains of bent, defective S-Mo-S layers embedded in an amorphous matrix, with excess sulfur and segregated molybdenum particles. Oxygen incorporation in this structure fosters the creation of amorphous oxide/oxysulfide nanophases with high electrical conductivity, enabling fast electron transfer to the active sites. The combined effect of the nanocrystalline pristine structure and the surface oxidation enhances the performance leading to small overpotentials, very fast kinetics (35.1 mV dec -1 Tafel slope) and remarkable long-term stability for continuous operation up to -1 A cm -2 . This work shows possible new avenues in catalytic design arising from a non-equilibrium technique as PLD and the importance of structural and chemical control to improve the HER performance of MoS-based catalysts.

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

NiFe Nanoparticle Nest Supported on Graphene as Electrocatalyst for Highly Efficient Oxygen Evolution Reaction

Abstract Designing cost‐efffective electrocatalysts for the oxygen evolution reaction (OER) holds significant importance in the progression of clean energy generation and efficient energy storage technologies, such as water splitting and rechargeable metal–air batteries. In this work, an OER electrocatalyst is developed using Ni and Fe precursors in combination with different proportions of graphene oxide. The catalyst synthesis involved a rapid reduction process, facilitated by adding sodium borohydride, which successfully formed NiFe nanoparticle nests on graphene support (NiFe NNG). The incorporation of graphene support enhances the catalytic activity, electron transferability, and electrical conductivity of the NiFe‐based catalyst. The NiFe NNG catalyst exhibits outstanding performance, characterized by a low overpotential of 292.3 mV and a Tafel slope of 48 mV dec −1 , achieved at a current density of 10 mA cm − 2 . Moreover, the catalyst exhibits remarkable stability over extended durations. The OER performance of NiFe NNG is on par with that of commercial IrO 2 in alkaline media. Such superb OER catalytic performance can be attributed to the synergistic effect between the NiFe nanoparticle nests and graphene, which arises from their large surface area and outstanding intrinsic catalytic activity. The excellent electrochemical properties of NiFe NNG hold great promise for further applications in energy storage and conversion devices.

Chemistry↗
DOE OSTI2024

P‐Doping Modulated RuIr Nanoparticles Anchored on Co/N/C Catalysts with Improved Alkaline Hydrogen Evolution Activity and Stability

Achieving efficient and stable hydrogen evolution reactions in alkaline conditions is crucial for hydrogen production. In this study, a RuIr/Co (SA) NC-P catalyst featuring RuIr alloys alongside P-doping and CoNx sites is developed. RuIr alloying optimizes the electronic structure between Ru and Ir, promoting electron transfer from Ru to Ir. P-doping further modulates the electronic properties of RuIr alloys, optimizing hydrogen binding energy and weakening Ru─OH binding energy, facilitating rapid H 2 generation and OH ad transfer. Meanwhile, CoNx promotes water dissociation, providing a fast proton delivery path for RuIr alloys. The catalyst exhibits enhanced HER activity with a low overpotential of 20 mV at 10 mA cm −2 , a Tafel slope of 20.4 mV dec −1 , and a turnover frequency of 19.5 H 2 s −1 at 150 mV overpotential. Moreover, catalyst stability is improved 8 times by mitigating RuIr alloy dissolution/agglomeration via P-doping. In conclusion, this work introduces a promising approach for developing efficient and stable HER electrocatalysts.

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

Alternating salt and freshwater floods of coastal soils impact soil structure, hydraulic properties, and oxygen dynamics

Coastal soils are increasingly impacted by hydrologic intensification in the form of rising sea level and flooding from storm surges and precipitation. Alternating saltwater (SW) and freshwater (FW) exposure has the potential to disperse colloids, which can lead to disintegration of soil structure, clogging of pore spaces, and reduction in ecologically and biogeochemically important functions like infiltration or gas exchange. To investigate how hydrologic intensification affects soil structure and oxygen dynamics, we conducted a series of laboratory-based flood simulations. Intact soil cores of the A and B horizons (22 total) from the toe slope of an upland coastal forest along the western shore of Chesapeake Bay. We subjected the cores to 24 h of saturation with either FW or alternating brackish SW and FW, with a 24-h draining event in between flood events. Oxygen diffusion into soil during draining was reduced by up to 30% for soils flooded alternating SW-FW compared to soils flooded only with FW. We attributed the reduced oxygen diffusion to clogging of smaller pores by colloids, with colloid redistribution observed as an increase in specific surface area down the core profile of up to 59%. After three SW-FW floods (six floods total), there were significant changes in pore size distribution, significant redistribution of colloids, and the A horizon became sodic. We concluded that a small number of SW flooding events can induce a measurable change in soil physical properties that directly impacts the biogeochemical dynamics.

Rod, Kenton A. [Pacific Northwest National Laborat↗
DOE OSTI2023

Veneziano variations: how unique are string amplitudes?

String theory offers an elegant and concrete realization of how to consistently couple states of arbitrarily high spin. But how unique is this construction? In this paper we derive a novel, multi-parameter family of four-point scattering amplitudes exhibiting i) polynomially bounded high-energy behavior and ii) exchange of an infinite tower of high-spin modes, albeit with a finite number of states at each resonance. These amplitudes take an infinite-product form and, depending on parameters, exhibit mass spectra that are either unbounded or bounded, thus corresponding to generalizations of the Veneziano and Coon amplitudes, respectively. For the bounded case, masses converge to an accumulation point, a peculiar feature seen in the Coon amplitude but more recently understood to arise naturally in string theory [1]. Importantly, our amplitudes contain free parameters allowing for the customization of the slope and offset of the spin-dependence in the Regge trajectory. We compute all partial waves for this multi-parameter class of amplitudes and identify unitary regions of parameter space. For the unbounded case, we apply similar methods to derive new deformations of the Veneziano and Virasoro-Shapiro amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Measurement of the central exclusive production of charged particle pairs in proton-proton collisions at $\sqrt{s}$ = 200 GeV with the STAR detector at RHIC

We report on the measurement of the Central Exclusive Production of charged particle pairs h+h- (h = π, K, p) with the STAR detector at RHIC in proton-proton collisions at $\sqrt{s}$ = 200 GeV. The charged particle pairs produced in the reaction pp → p' + h + h - + p' are reconstructed from the tracks in the central detector and identified using the specific energy loss and the time of flight method, while the forward-scattered protons are measured in the Roman Pot system. Exclusivity of the event is guaranteed by requiring the transverse momentum balance of all four final-state particles. Differential cross sections are measured as functions of observables related to the central hadronic final state and to the forward-scattered protons. They are measured in a fiducial region corresponding to the acceptance of the STAR detector and determined by the central particles’ transverse momenta and pseudorapidities as well as by the forward-scattered protons’ momenta. This fiducial region roughly corresponds to the square of the four-momentum transfers at the proton vertices in the range 0.04 GeV 2 < -t 1 , -t 2 < 0.2 GeV 2 , invariant masses of the charged particle pairs up to a few GeV and pseudorapidities of the centrally-produced hadrons in the range |η| < 0.7. The measured cross sections are compared to phenomenological predictions based on the Double Pomeron Exchange (DPE) model. Structures observed in the mass spectra of π + π - and K + K - pairs are consistent with the DPE model, while angular distributions of pions suggest a dominant spin-0 contribution to π + π - production. For π + π - production, the fiducial cross section is extrapolated to the Lorentz-invariant region, which allows decomposition of the invariant mass spectrum into continuum and resonant contributions. The extrapolated cross section is well described by the continuum production and at least three resonances, the f 0 (980), f 2 (1270) and f 0 (1500), with a possible small contribution from the f 0 (1370). Fits to the extrapolated differential cross section as a function of t 1 and t 2 enable extraction of the exponential slope parameters in several bins of the invariant mass of π + π - pairs. These parameters are sensitive to the size of the interaction region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Study of η′ → π+π−l+l− decays at BESIII

Abstract With a sample of (10087±44)×10 6 J/ψevents accumulated with the BESIII detector, we analyze the decaysη′→ π + π − l + l − (l=e, μ) via the processJ/ψ → γη′. The branching fractions are measured to be$$ \mathcal{B} $$ B (η′→ π + π − e + e − ) = (2.45±0.02(stat.)±0.08(syst.))×10 −3 and$$ \mathcal{B} $$ B (η′→ π + π − μ + μ − ) = (2.16±0.12(stat.)±0.06(syst.))×10 −5 , and the ratio is$$ \frac{\mathcal{B}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{e}^{+}{e}^{-}\right)}{\mathcal{B}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{\mu}^{+}{\mu}^{-}\right)}=113.4\pm 0.9\left(\textrm{stat}.\right)\pm 3.7\left(\textrm{syst}.\right) $$ B η ′ → π + π − e + e − B η ′ → π + π − μ + μ − = 113.4 ± 0.9 stat . ± 3.7 syst . . In addition, by combining theη′ →π + π − e + e − andη′ →π + π − μ + μ − decays, the slope parameter of the electromagnetic transition form factor is measured to beb η′ = 1.30 ± 0.19 (GeV/c 2 ) −2 , which is consistent with previous measurements from BESIII and theoretical predictions from the VMD model. The asymmetry in the angle between theπ + π − andl + l − decay planes, which has the potential to reveal theCP-violation originating from an unconventional electric dipole transition, is also investigated. The asymmetry parameters are determined to be$$ {\mathcal{A}}_{CP}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{e}^{+}{e}^{-}\right)=\left(-0.21\pm 0.73\left(\textrm{stat}.\right)\pm 0.01\left(\textrm{syst}.\right)\right)\% $$ A CP η ′ → π + π − e + e − = − 0.21 ± 0.73 stat . ± 0.01 syst . % and$$ {\mathcal{A}}_{CP}\left({\eta}^{\prime}\to {\pi}^{+}{\pi}^{-}{\mu}^{+}{\mu}^{-}\right)=\left(0.62\pm 4.71\left(\textrm{stat}.\right)\pm 0.08\left(\textrm{syst}.\right)\right)\% $$ A CP η ′ → π + π − μ + μ − = 0.62 ± 4.71 stat . ± 0.08 syst . % , implying that no evidence ofCP-violation is observed at the present statistics. Finally, an axion-like particle is searched for via the decayη′ →π + π − a, a→e + e − , and upper limits of the branching fractions are presented for the mass assumptions of the axion-like particle in the range of 0−500 MeV/c 2 .

Physics↗
DOE OSTI2025

Light neutral-meson production in pp collisions at $\sqrt{\text{s}}$ = 13 TeV

The momentum-differential invariant cross sections of π 0 and η mesons are reported for pp collisions at $\sqrt{s}$ = 13 TeV at midrapidity (|y| < 0.8). The measurement is performed in a broad transverse-momentum range of 0.2 < p T < 200 GeV/c and 0.4 < p T < 60 GeV/c for the π 0 and η, respectively, extending the p T coverage of previous measurements. Transverse-mass-scaling violation of up to 60% at low transverse momentum has been observed, agreeing with measurements at lower collision energies. Transverse Bjorken x (x T ) scaling of the π 0 cross sections at LHC energies is fulfilled with a power-law exponent of n = 5.01 ± 0.05, consistent with values obtained for charged pions at similar collision energies. The data are compared to predictions from next-to-leading order perturbative QCD calculations, where the π 0 spectrum is best described using the CT18 parton distribution function and the NNFF1.0 or BDSS fragmentation function. Expectations from PYTHIA8 and EPOS LHC overestimate the spectrum for the π 0 and are not able to describe the shape and magnitude of the η spectrum. The charged-particle multiplicity dependent π 0 and η p T spectra show the expected change of the spectral shape, characterized by a flatter slope with increasing multiplicity. This is demonstrated across a broad transverse-momentum range and up to events with a charged-particle multiplicity exceeding five times the mean value in minimum bias collisions. The η/π 0 ratio depends on the charged-particle multiplicity for p T < 4 GeV/c. PYTHIA8 and EPOS LHC qualitatively explain this behavior with an increasing contribution from the feed-down of heavier particles to the π 0 spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Corners and islands in the S-matrix bootstrap of the open superstring

We bootstrap the Veneziano superstring amplitude in 10 dimensions from the bottom-up. Starting with the most general maximally supersymmetric Yang-Mills EFT, we input information about the lowest-lying massive states, which we assume contribute via tree-level exchanges to the 4-point amplitude. We show the following: (1) if there is only a single state at the lowest mass, it must be a scalar. (2) Assuming a string-inspired gap between the mass of this scalar and any other massive states, the allowed region of Wilson coefficients has a new sharp corner where the Veneziano amplitude is located. (3) Upon fixing the next massive state to be a vector, the EFT bounds have a one-parameter family of corners; these would correspond to models with linear Regge trajectories of varying slopes, one of which is the open superstring. (4) When the ratio between the massive scalar coupling and the tr F 4 coefficient is fixed to its string value, the spin and mass of the second massive state is determined by the bootstrap and the Veneziano amplitude is isolated on a small island in parameter space. Finally, we compare with other recent bootstraps approaches, both the pion model and imposing Regge-inspired maximal spin constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Relations between b → cτν decay modes in scalar models

As a consequence of the Ward identity for hadronic matrix elements, we find relations between the differential decay rates of semileptonic decay modes with the underlying quark-level transition b → cτν, which are valid in scalar models. The decay-mode dependent scalar form factor is the only necessary theoretical ingredient for the relations. Otherwise, they combine measurable decay rates as a function of the invariant mass-squared of the lepton pair $q^2$ in such a way that a universal decay-mode independent function is found for decays to vector and pseudoscalar mesons, respectively. This can be applied to the decays $B → D^{\ast}τv, B_s → D^\ast_sτv, B_c → J/ψτv$ and $B → Dτv, B_s → D_sτv, B_c → η_cτv$, with implications for $R(D^{(\ast)}), R(D^{(\ast)}_s), R(J/ψ), R(η_c)$, and $\mathcal{B}(B_c → τv)$. The slope and curvature of the characteristic $q^2$-dependence is proportional to scalar new physics parameters, facilitating their straight forward extraction, complementary to global fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗
DOE OSTI2022

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Prediction of gas hydrate saturation using machine learning and optimal set of well-logs

We report resistivity and acoustic logs are widely used to estimate gas hydrate saturation in various sedimentary systems using one of the two popular methods ((1) acoustic velocity and (2) electrical resistivity), but the limitations of these two methods are often overlooked, which include (i) well-specific calibration of empirical exponents in the electrical resistivity method, (ii) assumption of known pore morphology for gas hydrates in the acoustic velocity method, and (iii) presence of unknown mineralogy and bulk modulus terms in the acoustic velocity method. NMR-density porosity-derived gas hydrate saturation based on the analysis of the transverse magnetization relaxation time (T2) is considered the most precise method, but acquisition of NMR-based logs is limited at relatively recent drilled sites; additionally, its use in conventional oil and gas reservoirs is not that common due to higher cost and operational deployment limitations associated with acquiring NMR well-logs. This study proposes a new method that predicts gas hydrate saturation (S h ) for any well using porosity, bulk density, and compressional wave (P wave) velocity well-logs with neural network (or stochastic gradient descent regression) without any well-specific calibration and/or other aforementioned shortcomings of the existing methods. The method is developed by examining the underlying dependency between S h and different combinations of well-logs, chosen from 6 routine logs, with 12 different machine learning (ML) algorithms. The accuracy of the proposed method in predicting S h is ~ 84%, which is better than the accuracy of seismic and electrical resistivity methods (≤ 75%) per the results reported by three different studies. The robustness of the method in the specific case of permafrost-associated gas hydrates is demonstrated with well-log data from two wells drilled on the Alaska North Slope.

58 GEOSCIENCES↗