Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computing Methodologies”

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 307 records · Page 17

Intelligent Reservoir Generation for Liquid State Machines using Evolutionary Optimization

Neuromorphic Computing is a burgeoning field of research. Many groups are exploring hardware architectures and theoretical ideas about spiking recurrent neural networks. The overarching goal is to exploit the low power promise of these neuromorphic systems. However, it is difficult to train spiking recurrent neural networks (SRNNs) to perform tasks and make efficient use of neuromorphic hardware. Reservoir Computing is an attractive methodology because it requires no tuning of weights for the reservoir itself. Yet, to find optimal reservoirs, manual tuning of hyperparameters such as hidden neurons, synaptic density, and natural structure is still required. Because of this, researchers often have to generate and evaluate many networks, which can result in non-trivial amounts of computation. This paper employs the reservoir computing technique (specifically liquid state machines) and genetic algorithms in order to develop useful networks that can be deployed on neuromorphic hardware. We build on past work in reservoir computing and genetic algorithms to demonstrate the power of combining these two techniques and the advantage it can provide over manually tuning reservoirs for use on classification tasks. We discuss the complexities of determining whether or not to use the genetic algorithms approach for liquid state machine generation.

Reynolds, John↗

SU-E-T-559: Monte Carlo Simulation of Cobalt-60 Teletherapy Unit Modeling In-Field and Out-Of-Field Doses for Applications in Computational Radiation Dosimetry

Purpose: Mean organ doses from structures located in field and outside of field boundaries during radiotherapy treatment must be considered when looking at secondary effects. Treatment planning in patients with 40 years of follow-up does not include 3-D treatment planning images and did not estimate dose to structures out of the direct field. Therefore, it is of interest to correlate actual clinical events with doses received. Methods: Accurate models of radiotherapy machines combined with whole body computational phantoms using Monte Carlo methods allow for dose reconstructions intended for studies on late radiation effects. The Theratron-780 radiotherapy unit and anatomically realistic hybrid computational phantoms are modeled in the Monte Carlo radiation transport code MCNPX. The major components of the machine including the source capsule, lead in the unit-head, collimators (fixed/adjustable), and trimmer bars are simulated. The MCNPX transport code is used to compare calculated values in a water phantom with published data from BJR suppl. 25 for in-field doses and experimental data from AAPM Task Group No. 36 for out-of-field doses. Next, the validated cobalt-60 teletherapy model is combined with the UF/NCI Family of Reference Hybrid Computational Phantoms as a methodology for estimating organ doses. Results: The model of Theratron-780 has shown to be agree with percentage depth dose data within approximately 1% and for out of field doses the machine is shown to agree within 8.8%. Organ doses are reported for reference hybrid phantoms. Conclusion: Combining the UF/NCI Family of Reference Hybrid Computational Phantoms along with a validated model of the Theratron-780 allows for organ dose estimates of both in-field and out-of-field organs. By changing field size, position, and adding patient-specific blocking more complicated treatment set-ups can be recreated for patients treated historically, particularly those who lack both 2D/3D image sets.

61 RADIATION PROTECTION AND DOSIMETRY↗

Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models

The document is the Final Technical Progress Report for the Nuclear Energy University Program’s Integrated Research Project (IRP) on “Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models” that supports the LWR Sustainability Program’s RISMC R&D Pathway. The project goal is to develop and demonstrate a data-driven methodology for validation of advanced computer models used in nuclear power plant safety analysis. Specifically, the advanced computer models are those in the toolkit developed to support risk-informed safety margin characterization (RISMC), an integrated deterministic/probabilistic safety analysis methodology developed in the Department of Energy’s Light Water Reactor Sustainability (LWR-S) program. The report reflects the progress made towards the project’s stated goal by contributions by researchers and graduate students from universities, with support from researchers from national laboratories and industry companies. The project organization, effort coordination and technical implementation are summarized, followed by discussion of main findings, issues, and path forward. Selected chapters provide a more detailed description of tasks, approaches and respective findings and recommendations. Noteworthy are contributions that serve as guidelines for methodology development. It is also noted that this report is complemented by other milestone reports (as stand-alone deliverables) that provide detailed discussion of the technical developments. The project results have been documented in a number (12) dissertations and these, 50+ peer-reviewed publications in technical journals and conference proceedings.

42 ENGINEERING↗

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗

Comparison of neutron organ and effective dose coefficients for PIMAL stylized phantom in bent postures in standard irradiation geometries

Neutron dose coefficients for standard irradiation geometries have been reported in International Commission on Radiological Protection (ICRP) Publication 116 for the ICRP Publication 110 adult reference phantoms. In the present paper, organ and effective dose coefficients have been calculated for a receptor in both upright and articulated (bent) postures representing more realistic working postures exposed to a mono-energetic neutron radiation field. This work builds upon prior work by Dewji and co-workers comparing upright and bent postures for exposure to mono-energetic photon fields. Simulations were conducted using the Oak Ridge National Laboratory’s articulated stylized adult phantom, “Phantom wIth Moving Arms and Legs” (PIMAL) software package, and the Monte Carlo N-Particle (MCNP) version 6.1.1 radiation transport code. Organ doses were compared for the upright and bent (45° and 90°) phantom postures for neutron energies ranging from 1 × 10 -9 to 20 MeV for the ICRP Publication 116 external exposure geometries—antero-posterior (AP), postero-anterior (PA), and left and right lateral (LLAT, RLAT). Using both male and female phantoms, effective dose coefficients were computed using ICRP Publication 103 methodology. The resulting coefficients for articulated phantoms were compared to those of the upright phantom. Computed organ and effective dose coefficients are discussed as a function of neutron energy, phantom posture, and source irradiation geometry. For example, it is shown here that for the AP and PA irradiation geometries, the differences in the organ coefficients between the upright and bent posture become more pronounced with increasing bending angle. In the AP geometry, the brain dose coefficients are expectedly higher in the bent postures than in the upright posture, while all other organs have lower dose coefficients, with the thyroid showing the greatest difference. Finally and overall, the effective dose estimated for the upright phantom is more conservative than that for the articulated phantom, which may have ramifications in the estimation or reconstruction of radiation doses.

61 RADIATION PROTECTION AND DOSIMETRY↗

Road Lidar Dataset for the TxDOT Austin District

This is a road lidar data collection for developing road elevation models and road inundation mapping methodologies, a joint work between ORNL and The University of Texas at Austin. This dataset is generated as part of the flood transportation infrastructure, partly funded by the NOAA CIROH project. ORNL is a project partner for high-performance computing-empowered flood inundation mapping methodology R&D. The dataset is computed using a GPU-accelerated lidar data processing workflow developed at ORNL. The lidar data source is from TxGIO, the state lidar data collection site. The output dataset is in two formats: laz and copc. It is organized by TxDOT's maintenance sections, covering the Austin District. Data size: 3.86 billion road lidar points, 1.67% of the entire lidar data input Projection: EPSG:32614 (WGS84/UTM zone 14N) Website: https://web.corral.tacc.utexas.edu/nfiedata/road3d/austin_district/AustinMaintenanceSections_H_epsg6343_V_epsg5703/ LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

Earth Sciences Applications for Energy and Global Security: A Few Vignettes [Slides]

Understanding the physics of flow and transport in the subsurface environment is crucial to several National Security applications. In the Energy Security realm, some examples include characterizing radioactive waste disposal, exploring natural gas extraction and sequestering CO2. An example in the Global Security arena is to be able to distinguish between anthropogenic and naturally occurring seismic disturbances, the former being of particular concern as underground nuclear testing activity has seen an increase in recent years. A common theme across these applications is the ability to accurately model fluid flow in porous and/or fractured subsurface media. For decades, researchers have grappled with how to include the effects of uncertainties in these models, including but not limited to heterogeneity, unknown initial and boundary conditions and issues with scalability. More recently advances in computing and machine learning methodologies have allowed a broader exploration of the uncertainty space as we strive for real time decision making. I will present an overview of the field and discuss a few examples of how my own research in these topics has evolved over nearly two decades.

58 GEOSCIENCES↗

Building and testing correlations for the estimation of one‐electron reduction potentials of a diverse set of organic molecules

We describe and evaluate a method for computationally predicting reduction potentials of a diverse group of organic molecules by linearly correlating calculated lowest unoccupied molecular orbital energies with ground state reduction potentials measured in acetonitrile. The approach is shown to provide a unique combination of extreme computational simplicity and excellent accuracy across a diverse range of organic structures and a wide window of reduction potentials. A disparate set of molecules (74 compounds belonging to six distinct structural families, comprised of molecules containing C, H, N, O, F, Cl, and Br, with functional groups including esters, ketones, halides, nitriles, quinones, alkenes, arenes, heteroarenes, and pyridinium and higher benzologs, all containing conjugated pi systems, spanning a 3.5‐V range of reduction potentials) was used to build the correlations. This methodology was found to be computationally inexpensive compared with other approaches and to permit the useful prediction of reduction potentials of additional molecules of diverse structural types not included in the families used to determine the correlation parameters. The effects of varying the basis set used in the B3LYP electronic structure calculations and including solvent (compared with calculations in gas phase) were also examined. It was found that the inclusion of a continuum solvent model in the calculations was required for accurate results, particularly when including cationic species in the correlations (although when only neutral molecules were examined, reasonable results could even be obtained in vacuo ). Several of the best correlations were used to predict the reduction potentials of seven much larger and structurally diverse chromophores that were not included in the correlation data set. Strong correlations ( r 2 values > 0.99) with very good predictive abilities (root mean square deviation < 60 mV) were found. This extremely simple and computationally efficient entirely closed‐shell methodology is proven robust and useful for the design of new molecules capable of participating in redox processes, including electron transfer reactions. Copyright © 2015 John Wiley & Sons, Ltd.

Méndez‐Hernández, Dalvin D.↗

PNNL Strategic Goods Testbed: A Data Library for Illicit Nuclear Trafficking

Pacific Northwest National Laboratory (PNNL) has put significant effort into nonproliferation activities as an institution, both in terms of the classical nuclear material focused approach and in the examination of other strategic goods necessary to implement a nuclear program. To assist in these efforts, several projects in the Analysis in Motion (AIM) and Signature Discovery (SDI) Initiatives at PNNL are developing machine learning methodology for human-computer interaction in real time environments to assist analysts in this domain. All of these technical projects require access to data – whether it is in terms of detector data, shipping records, financial information, company relations, or other communications. The first question that mathematical and computational researchers come up with when asked to build analyst assist or automated tools is “What does the data look like? ” They become frustrated when basic questions like this can not be easily answered and this can have the effect of pushing researchers away from the nuclear trafficking domain, especially in strategic commodity and export control areas where data sets can not easily be generated through standard experimental techniques. For small projects that are building a proof of concept for their methodology, obtaining this data can be arduous and expensive. To relieve the burden of data collection from these projects and grow a lab-wide capability, the Strategic Goods Testbed Team has taken over data collection and placed subscriptions and access to flat data files in a centralized location so that all projects can benefit from these items. We have collected shipping data in the form of PIERS records, judicial information about export control cases, NAC data on the nuclear fuel industry, and financial data from Dun and Bradstreet and our data sets are continuing to expand. With a single access agreement, researchers in data-mining and other fields can utilize all of the records that have been downloaded, make requests through subscription services, and interact with other researchers through our interface. Our testbed team provides more than a simply static repository by working with researchers to refine their data needs and insure data quality as well as quantity. We are currently working with laboratory and initiative specific management to examine effective ways for continuing data growth and sustainability.

data capability↗

TU-A-12A-08: Computing Longitudinal Material Changes in Bone Metastases Using Dual Energy Computed Tomography

Purpose: This study demonstrates a methodology for tracking changes in metastatic bone disease using trajectories in material basis space in serial dual energy computed tomography (DECT) studies. Methods: This study includes patients with bone metastases from breast cancer that had clinical surveillance CT scans using a General Electric CT750HD in dual energy mode. A radiologist defined regions-of-interested (ROI) for bone metastasis, normal bone, and marrow across the serial DECT scans. Our approach employs a Radon transform to forward-projection the basis images, namely, water and iodine, into sinogram space. This data is then repartitioned into fat/bone and effective density/Z image pairs using assumed energy spectrums for the x-ray energies. This approach both helps remove negative material densities and avoids adding spectrum-hardening artifacts. These new basis data sets were then reconstructed via filtered back-projection to create new material basis pair images. The trajectories of these pairs were then plotted in the new basis space providing a means to both visualize and quantitatively measure changes in the material properties of the tumors. Results: ROI containing radiologist defined metastatic bone disease showed well-defined trajectories in both fat/bone and effective density/Z space. ROI that contained radiologist defined normal bone and marrow did not exhibit any discernible trajectories and were stable from scan to scan. Conclusions: The preliminary results show that changes in material composition and effective density/Z image pairs were seen primarily in metastasis and not in normal tissue. This study indicates that by using routine clinical DECT it may be possible to monitor therapy response of bone metastases because healing or worsening bone metastases change material composition of bone. Additional studies are needed to further validate these results and to test for their correlation with outcome.

60 APPLIED LIFE SCIENCES↗

Matrix-Based Process Modeling in Microsoft Excel

The purpose of this document is to describe the theory and methodology of matrix-based process modeling and demonstrate its use through a generic process model within Microsoft Excel. This type of model allows a modeler to quickly query various metrics, such as total production time, for a specified production run. These matrix models tend to be quick to build and use, but come with certain limitations and are more suited for smaller processes.

97 MATHEMATICS AND COMPUTING↗

Activity Theory Literature Review

Complex challenges across Sandia National Laboratories' (SNL) mission areas underscore the need for systems level thinking, resulting in a better understanding of the organizational work systems and environments in which our hardware and software will be used. SNL researchers have successfully used Activity Theory (AT) as a framework to clarify work systems, informing product design, delivery, acceptance, and use. To increase familiarity with AT, a working group assembled to select key resources on the topic and generate an annotated bibliography. The resources in this bibliography are arranged in six categories: 1) An introduction to AT; 2) Advanced readings in AT; 3) AT and human computer interaction (HCI); 4) Methodological resources for practitioners; 5) Case studies; and 6) Related frameworks that have been used to study work systems. This annotated bibliography is expected to improve the reader's understanding of AT and enable more efficient and effective application of it.

42 ENGINEERING↗

Toward Interactive Scenario Analysis and Exploration

As Modeling and Simulation (M&S) tools have matured, their applicability and importance have increased across many national security challenges. In particular, they provide a way to test how something may behave without the need to do real world testing. However, current and future changes across several factors including capabilities, policy, and funding are driving a need for rapid response or evaluation in ways that many M&S tools cannot address. Issues around large data, computational requirements, delivery mechanisms, and analyst involvement already exist and pose significant challenges. Furthermore, rising expectations, rising input complexity, and increasing depth of analysis will only increase the difficulty of these challenges. In this study we examine whether innovations in M&S software coupled with advances in ''cloud'' computing and ''big-data'' methodologies can overcome many of these challenges. In particular, we propose a simple, horizontally-scalable distributed computing environment that could provide the foundation (i.e. ''cloud'') for next-generation M&S-based applications based on the notion of ''parallel multi-simulation''. In our context, the goal of parallel multi- simulation is to consider as many simultaneous paths of execution as possible. Therefore, with sufficient resources, the complexity is dominated by the cost of single scenario runs as opposed to the number of runs required. We show the feasibility of this architecture through a stable prototype implementation coupled with the Umbra Simulation Framework [6]. Finally, we highlight the utility through multiple novel analysis tools and by showing the performance improvement compared to existing tools.

97 MATHEMATICS AND COMPUTING↗

Development and Integration of Predictive Models for Manufacturing and Structural Performance of Carbon Fiber Composites in Automotive Applications (Final Report)

The goal of this project is to develop integrated, state-of-the-art, computational modeling tools based on an integrated computational materials engineering (ICME) methodology that are critically needed to enable structural carbon fiber (CF) applications in automobiles. These tools are designed to predict the manufacturing and structural performance of CF composites, including stochastic effects. During the first phase of the project, the manufacturing and structural performance tools, including a stochastic driver, were developed, calibrated, and validated against coupon and component level test results. For this project, from the manufacturing side, the development efforts were focused on high-pressure resin transfer molding (HP-RTM), a potentially game-changing manufacturing technology that is capable of reaching the 3-5-minute cycle time that is necessary for volume manufacturing of composites. On the structural performance side, the crashworthiness of CF structures was studied. Accordingly, computational tools were developed and validated. The resulting differences between the numerical predictions and the experimental results in the first phase of the project were required to be less than 15%. During the second phase of the project, the manufacturing and performance tools were integrated by mapping the manufacturing outcomes (e.g., fiber angles, residual stresses, degree of cure, and defects) into the structural models. Also, using the tools developed in this project, an automotive assembly currently manufactured in steel (2016 GM Malibu) was redesigned to be manufactured using advanced CF composites with the objective of comparing the performance of the CF assembly with the developed ICME model predictions. The weight, performance and cost of the CF assembly were compared with the baseline steel assembly as well. The CF underbody design was shown to be 30% lighter than the corresponding steel design. The CF underbody assembly withstood the side pole impact with less than half the intrusion for the steel assembly. Excellent correlations within 10% for both peak impact load and intrusion were observed between the numerical predictions and experimental results. Based on the cost models developed in the project, the cost increase per Kg saved was determined to be $22. The project team has prepared a total of 28 publications in various journals and conferences and broadly communicated the project results for the benefit of the composite industry. The project team has also prepared a total of 10 patent applications that were submitted to the US Patent and Trademark office.

36 MATERIALS SCIENCE↗

A multifractal approach to space-filling recovery for PET quantification

Purpose: A new image-based methodology is developed for estimating the apparent space-filling properties of an object of interest in PET imaging without need for a robust segmentation step and used to recover accurate estimates of total lesion activity (TLA). Methods: A multifractal approach and the fractal dimension are proposed to recover the apparent space-filling index of a lesion (tumor volume, TV) embedded in nonzero background. A practical implementation is proposed, and the index is subsequently used with mean standardized uptake value (SUV {sub mean}) to correct TLA estimates obtained from approximate lesion contours. The methodology is illustrated on fractal and synthetic objects contaminated by partial volume effects (PVEs), validated on realistic {sup 18}F-fluorodeoxyglucose PET simulations and tested for its robustness using a clinical {sup 18}F-fluorothymidine PET test–retest dataset. Results: TLA estimates were stable for a range of resolutions typical in PET oncology (4–6 mm). By contrast, the space-filling index and intensity estimates were resolution dependent. TLA was generally recovered within 15% of ground truth on postfiltered PET images affected by PVEs. Volumes were recovered within 15% variability in the repeatability study. Results indicated that TLA is a more robust index than other traditional metrics such as SUV {sub mean} or TV measurements across imaging protocols. Conclusions: The fractal procedure reported here is proposed as a simple and effective computational alternative to existing methodologies which require the incorporation of image preprocessing steps (i.e., partial volume correction and automatic segmentation) prior to quantification.

60 APPLIED LIFE SCIENCES↗

Automated Programmable Logic Controller Memory Forensics Using RGB Image Analysis and Deep Learning

The introduction of Industry 4.0 and Internet-based technologies has enhanced industrial control system operations but have inadvertently increased their vulnerabilities to cyber attacks. When an industrial control system is compromised, security analysts need to identify the root cause quickly to start the recovery process and develop mitigation strategies. Memory forensics is critical in the incident analysis process to ascertain what occurred. Approaches for analyzing the persistent memory in industrial control devices are limited and almost nonexistent for volatile memory. This chapter proposes an automated methodology for programmable logic controller memory dump analysis using computer vision and deep learning techniques. The methodology converts the sequences of bytes in a programmable logic controller memory dump to red-green-blue pixels and employs a deep learning model that learns the underlying patterns and features of pre-labeled forensic artifacts in images and segments them into distinct regions. The trained model is employed to automatically segment new memory images and identify forensic artifacts. Evaluation of the methodology on a Schneider Electric Modicon M221 programmable logic controller under code injection and code modification attacks demonstrates its ability to detect attack artifacts in memory dumps.

Asmar Awad, Rima [ORNL] (ORCID:0000000233407742)↗

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↗