Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive Data analytics”

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 505 records · Page 28

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗

Jet substructure from dark sector showers

We examine the robustness of collider phenomenology predictions for a dark sector scenario with QCD-like properties. Pair production of dark quarks at the LHC can result in a wide variety of signatures, depending on the details of the new physics model. A particularly challenging signal results when prompt production induces a parton shower that yields a high multiplicity of collimated dark hadrons with subsequent decays to Standard Model hadrons. The final states contain jets whose substructure encodes their non-QCD origin. This is a relatively subtle signature of strongly coupled beyond the Standard Model dynamics, and thus it is crucial that analyses incorporate systematic errors to account for the approximations that are being made when modeling the signal. We estimate theoretical uncertainties for a canonical substructure observable designed to be sensitive to the gauge structure of the underlying object, the two-point energy correlator $ {e}_2^{\left(\beta \right)} $ , by computing envelopes between resummed analytic distributions and numerical results from Pythia. We explore the separability against the QCD background as the confinement scale, number of colors, number of flavors, and dark quark masses are varied. Additionally, we investigate the uncertainties inherent to modeling dark sector hadronization. Simple estimates are provided that quantify one’s ability to distinguish these dark sector jets from the overwhelming QCD background. Such a search would benefit from theory advances to improve the predictions, and the increase in statistics using the data to be collected at the high luminosity LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Chemical Differences between Phenolic Secondary Organic Aerosol Formed through Gas-Phase and Aqueous-Phase Reactions

Phenolic compounds, which are significant emissions from biomass burning (BB), undergo rapid photochemical reactions in both gas and aqueous phases to form secondary organic aerosol, namely, gasSOA and aqSOA, respectively. The formation of gasSOA and aqSOA involves different reaction mechanisms, leading to different product distributions. In this study, we investigate the gaseous and aqueous reactions of guaiacol a representative BB phenol to elucidate the compositional differences between phenolic aqSOA and gasSOA. Aqueous-phase reactions of guaiacol produce higher SOA yields than gas-phase reactions (e.g., roughly 60 vs 30% at one half-life of guaiacol). These aqueous reactions involve more complex reaction mechanisms and exhibit a more gradual SOA evolution than their gaseous counterparts. Initially, gasSOA forms with high oxidation levels (O/C > 0.82), while aqSOA starts with lower O/C (0.55–0.75). However, prolonged aqueous-phase reactions substantially increase the oxidation state of aqSOA, making its bulk chemical composition closer to that of gasSOA. Additionally, aqueous reactions form a greater abundance of oligomers and high-molecular-weight compounds, alongside a more sustained production of carboxylic acids. AMS spectral signatures representative of phenolic gasSOA have been identified, which, together with tracer ions of aqSOA, can aid in the interpretation of field observation data on aerosol aging within BB smoke. The notable chemical differences between phenolic gasSOA and aqSOA highlighted in this study also underscore the importance of accurately representing both pathways in atmospheric models to better predict the aerosol properties and their environmental impacts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating magnetic van der Waals materials using data-driven approaches

In this work, we investigate magnetic monolayers of the form A i A ii B 4 X 8 based on the well-known intrinsic topological magnetic van der Waals (vdW) material MnBi 2 Te 4 (MBT) using first-principles calculations and machine learning techniques. We select an initial subset of structures to calculate the thermodynamic properties, electronic properties, such as the band gap, and magnetic properties, such as the magnetic moment and magnetic order using density functional theory (DFT). Data analytics approaches are used to gain insight into the microscopic origin of materials’ properties. The dependence of materials’ properties on chemical composition is also explored. For example, we find that the formation energy and magnetic moment depend largely on A and B sites whereas the band gap depends on all three sites. Finally, we employ machine learning tools to accelerate the search for novel vdW magnets in the MBT family with optimized properties. Finally, this study creates avenues for rapidly predicting novel materials with desirable properties that could enable applications in spintronics, optoelectronics, and quantum computing.

36 MATERIALS SCIENCE↗

Modeling of Additively Manufactured Ceramic Heat Exchangers with Semi-Elliptical Cross-Section Flow Channels

An approximate and easily applied analytical model was developed for heat transfer calculations of heat exchangers consisting of multiple rows and columns of heat transfer fluid flow channels with semi-elliptical cross sections. Heat exchangers of this type are being developed by using ceramic material and additive manufacturing for high temperature and pressure-concentrating solar electric power plants. Calculations using the model require only the geometrical dimensions and flow conditions of the heat exchanger. Comparisons of modeling predictions, with both simulation results and experimental data, were conducted to verify the viability of the model. The results showed good agreement where almost all the modeling predictions were within 20% of the simulation results or the experimental data. Finally, the proposed modeling approach is more generally applicable to heat transfer analysis of heat exchangers with similar flow channel configurations to those considered in this study.

Ceramic↗

Reexamining supercritical gas adsorption theories in nano-porous shales under geological conditions

Adsorption is an important phenomenon in surface chemistry, especially for nano-porous shales. In the shale gas-in-place, adsorbed gas could contribute up to 85%. However, adsorption is hard to quantitatively characterize due to different adsorption mechanisms, patterns, surfaces, and pore sizes. Moreover, key thermodynamic parameters, such as the enthalpy of adsorption, are challenging to determine due to uncertainties in adsorbed gas densities used in constructing absolute isotherms. In this paper, we revisit the Brunauer, Emmett, and Teller (BET) model, analytically simplify Ono-Kondo (OK) models for subsurface shales, and compare commonly used mono- and multilayer adsorption models (e.g., Langmuir, supercritical Dubinin-Radushkevich (SDR), supercritical BET (SBET), and simplified OK (OKs) models) with recently proposed pressure-dependent adsorption densities to develop a practical and reliable methodology that can be used in the supercritical state, typical for subsurface black shale conditions.Three independent data sets were used for nitrogen and methane adsorption isotherms at different temperatures. Furthermore, we demonstrate that adsorption predicted by the SDR model is comparable or lower than that predicted by the SBET model, but higher than the amounts predicted by Langmuir and OK models. The nitrogen BET method tends to underestimate the accessible SSA for methane. Measurement of isosteric heat of adsorption is suggested instead of the experimental fitting method, due to the significant difference between the calculated results by two commonly used methods. Experimental fitting and simulation methods are also briefly reviewed to guide future research on shale gas adsorption.

58 GEOSCIENCES↗

Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials

We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.

CREASE↗

Statistically-informed deep learning for gravitational wave parameter estimation

We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_\mathrm f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell = m = 2$ bar mode, $(\omega_\mathrm R, \omega_\mathrm I)$. Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_\mathrm f, \omega_\mathrm R, \omega_\mathrm I)$ of five binary black holes: GW150914, GW170104, GW170814, GW190521 and GW190630. We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. Furthermore, this neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science.

79 ASTRONOMY AND ASTROPHYSICS↗

Complete Round-Robin Hydrogen Gas Analysis Capability Comparison (Milestone 2.6)

The Department of Energy (DOE) is currently evaluating strategies for the extended dry storage of aluminum-clad spent nuclear fuel (ASNF). Part of this assessment concerns the extent of radiolytic molecular hydrogen (H2) generation from the aluminum cladding’s oxyhydroxide corrosion layers. Understanding this radiation-induced process and the factors effecting it (e.g., system conditions such as temperature and gaseous environment) are essential for the development of predictive computer models to support the Technical Considerations and Challenges for Extended (> 50 yrs) Dry Storage of ASNF program. To achieve this goal and ensure that the experimental data gathered by Task 2 (Oxyhydroxide Layer Radiolytic Gas Generation Resolution) research groups (Idaho National Laboratory and Savannah River National Laboratory) are consistent, a round-robin H2 analysis capability comparison was initiated. Here we present the results from said round-robin and conclude that despite differences in sample preparation, irradiation parameters, and analytical procedures, the measured data are sufficiently consistent between the two laboratories (= 15%).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Study of the deformed halo nucleus 31 Ne with Glauber model based on microscopic self-consistent structures

In this work, the exotic deformed nucleus 31 Ne is studied with an approach that combines self-consistent structure and reaction theories. The deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) is utilized to demonstrate that deformation and pairing correlations give rise to a halo structure with a large-amplitude p-wave configuration in 31 Ne. Then the valence nucleon wave functions and angle-averaged density distributions of 30 Ne from this theory are used as input for the Glauber reaction model to study the observables of neutron-rich neon isotopes to search halo signatures. With NL1 effective interaction, our predictions of the reaction cross sections for these exotic neon isotopes on a carbon target can better reproduce the experimental data than those from the relativistic mean field model for a spherical shape with resonances and pairing correlation contributions, and are roughly 3.3% (~ 50 mb) larger than those with Gaussian function-fitted densities of the core nuclei. The calculated one-neutron removal cross section at 240 MeV/nucleon, and the inclusive longitudinal momentum distribution of the 30 Ne residues from the 31Ne breakup reaction are largely improved over previous theoretical predictions and agree well with data. These reaction evaluations indicate a dilute density distribution in coordinate space and are a canonical signature of a halo structure. Moreover, our predictions with the NL3 and PK1 effective interactions give slightly better descriptions of reaction observables for exotic neon isotopes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unbundling Smart Meter Services Through Spatio-Temporal Decomposition Agents in DER-rich Environment

Smart meters and the advanced metering infrastructure (AMI) facilitate distribution system operators (DSOs) to gather information on energy consumption at the customer level. With the increasing penetration of building-level intermittent distributed energy resources (DERs) behind the meter, DER information is not available to DSOs. At the same time, smart meter enables users to participate in grid, with real-time information. Information for behind the meter is needed by user to coordinate building level assets for maximum benefits. The concept of unbundled smart meter (USM) needs agents to decompose smart meter measurements to provide service to DSO as well as customers. In this paper, we propose a Spatio-Temporal Decomposition Agent (STDA) for USM based on Artificial Intelligence (AI). STDA can help users optimize their energy usage, help DSO to utilize building assets for the grid operation. The energy usage strategy developed by STDA is suitable for different users, and can be customized by deep learning (DL) models according to the different energy consumption habits of each user. The power prediction performance results of various DL models and evaluation using a set of data from a Hawaii utility is presented. Furthermore, STDA integration with Home Energy management Systems (HEMS) to manage resources is presented and validated. STDA pre-processes the measurements before model training, and provides the spatio-temporal decomposed forecasting.

42 ENGINEERING↗

The role of detonation condensates on the performance of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) detonation

Thermochemical models of detonation are widely used to estimate energy delivery, but they are based on the assumption that the carbon-rich condensates (soot) formed during detonation are very similar to bulk carbon. Here, we present an analytic equation of state (EOS) based on experimental detonation data for soot formed during the detonation of triaminotrinitrobenzene (TATB)-based high explosives. X-ray photoelectron spectra of several detonation soots are used to determine the elemental nitrogen abundance, with surprisingly high values for TATB. The proposed TATB soot EOS is highly compressible at low pressures and shares some features of glassy carbon, exhibiting graphite- and diamond-like behavior as a function of pressure. We demonstrate the influence of formed soot on detonation performance, including a lowering of the detonation velocity at typical charge densities, and a more compressive product Hugoniot at overdriven conditions. The soot model improves the accuracy of thermochemical calculations for TATB-based explosives across a wide range of states. Detonation velocity predictions for HMX (cyclotetramethylene-tetranitramine)-TATB blends with 80% or more TATB content, as well predictions for 1,3-diamino-2,4,6-trinitrobenzene (DATB) and 3-nitro-1,2,4-triazol-5-one (NTO), which share some features with TATB, are also improved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Precision beam diagnostics at the NuMI facility using muon monitor observations

The Neutrinos at the Main Injector (NuMI) facility at Fermilab delivers an intense neutrino beam for multiple experiments by producing pions that decay into neutrinos, muons, and other particles. Magnetic horns—the primary pion focusing elements in the NuMI beamline—exhibit predominantly linear optics, enabling a predictable relationship between the proton beam and the resulting pion and muon phase spaces. This study has two primary objectives: first, to evaluate and confirm the linearity of the horn focusing mechanism using analytical models and numerical simulations; and second, to demonstrate that key beam parameters—such as proton beam intensity, beam position on target, and horn current—can be extracted from muon monitor observations within this linear optics framework. Using a machine learning model trained on spill-by-spill muon monitor data, we infer the horn current with a precision of ±0.05%, the beam intensity with ±0.1%, and the beam position on target with ±0.018⁢ mm horizontally and ±0.013⁢ mm vertically. This approach provides a reliable cross-check of beam parameters, helping to reduce systematic uncertainties that are critical for future experiments such as the Deep Underground Neutrino Experiment, which will rely on the neutrino beam produced by the Long-Baseline Neutrino Facility.

Beam control↗

Verification, Validation, and Uncertainty Quantification in Thermal Hydraulics, Freeman Scholar Lecture (2019)

Engineering problems are generally solved by analytical models or computer codes. These models, in addition to conservation equations, also include many empirical relationships and approximate numerical methods. Each of these components contributes to the uncertainty in the prediction. A systematic approach to judge the applicability of the code to the intended application is needed. It starts from verification of implementation of formulation in the code, identification of important phenomena, finding relevant tests with quantified uncertainty for these phenomena, and validation of the code by comparing predictions with the relevant test data. The relevant tests must address phenomena as expected in the intended application. In case of small size or limited condition tests, the scaling analyses are needed to assess the relevancy of the tests. Finally, a statement of uncertainty in the prediction is needed. Systematic approaches are described to aggregate uncertainties from different components of the code for intended application. Here, verification, validation, and uncertainty quantifications (VVUQs) are briefly described.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Universal diamond edge Raman scale to 0.5 terapascal and implications for the metallization of hydrogen

The recent progress in generating static pressures up to terapascal values opens opportunities for studying novel materials with unusual properties, such as metallization of hydrogen and high-temperature superconductivity. However, an evaluation of pressure above ~0.3 terapascal is a challenge. We report a universal high-pressure scale up to ~0.5 terapascal, which is based on the shift of the Raman edge of stressed diamond anvils correlated with the equation of state of Au and does not require an additional pressure sensor. According to the new scale, the pressure values are substantially lower by 20% at ~0.5 terapascal compared to the extrapolation of the existing scales. We compare the available data of H 2 at the highest static pressures. We show that the onset of the proposed metallization of molecular hydrogen reported by different groups is consistent when corrected with the new scale and can be compared with various theoretical predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Limits on Mode Coherence in Pulsating DA White Dwarfs Due to a Nonstatic Convection Zone

The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of convection zones in pulsating hydrogen-atmosphere (DAV) white dwarfs is very sensitive to temperature (i.e., M CZ α T –90 eff ), leading to the possibility that even low-amplitude pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is out of phase with what is required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. In this paper we show that: (1) whether a mode is calculated to propagate to the base of the convection zone is an accurate predictor of its width in the Fourier spectrum, (2) the phase shifts produced by reflection from the outer turning point are large enough to produce significant damping, and (3) amplitudes and periods are predicted to increase from the blue edge to the middle of the instability strip, and subsequently decrease as the red edge is approached. Furthermore, this amplitude decrease is in agreement with the observational data while the period decrease has not yet been systematically studied.

79 ASTRONOMY AND ASTROPHYSICS↗

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗