Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “log response”

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.

504 records · Page 28

Retrieval of Ozone Column Content from Airborne Sun Photometer Measurements During SOLVE II: Comparison with SAGE III, POAM III,THOMAS and GOME Measurements

During the Second SAGE 111 Ozone Loss and Validation Experiment (SOLVE II), the 14- channel NASA Ames Airborne Trackmg Sunphotometer (AATS-14) was mounted on the NASA DC-8 and successfully measured spectra of total and aerosol optical depth (TOD and AOD) during the sunlit portions of eight science flights. Values of ozone column content above the aircraft have been derived from the AATS-14 data by using a linear least squares method. For each AATS-14 measured TOD spectrum, this method iteratively finds the ozone column content that yields the best match between measured and calculated TOD. The calculations assume the known Chappuis ozone band shape and a three-parameter AOD shape (quadratic in log-log space). Seven of the AATS-14 channels (each employing an interference filter with a nominal full-width at half maximum bandpass of -5 nm) are within the Chappuis band, with center wavelengths between 452.9 nm and 864.5 nm. One channel (604.4 nm) is near the peak, and three channels (499.4, 519.4 and 675.1 nm) have ozone absorption within 30-40% of that at the peak. For the typical DC-8 SOLVE II cruising altitudes of approx. 8-12 km and the background stratospheric aerosol conditions that prevailed during SOLVE 11, absorption of incoming solar radiation by ozone comprised a significant fraction of the aerosol-plus-ozone optical depth measured in the four AATS-14 channels centered between 499.4 and 675.1 nm. Typical AODs above the DC-8 ranged from 0.003-0.008 in these channels. For comparison, an ozone overburden of 0.3 atm-cm (300 DU) translates to ozone optical depths of 0.009,0.014, 0.041, and 0.012, respectively, at these same wavelengths. In this paper, we compare AATS-14 values of ozone column content with temporally and spatially near-coincident values derived from measurements acquired by the Stratospheric Aerosol and Gas Experiment III (SAGE III) and the Polar Ozone and Aerosol Measurement 111 (POAM III) satellite sensors. We also compare AATS-14 ozone retrievals during selected DC-8 latitudinal and longitudinal transects with total column ozone data acquired by the Total Ozone Mapping Spectrometer (TOMS) and the Global Ozone Monitoring Experiment (GOME) satellite sensors. To enable this comparison, the amount of ozone in the column below the aircraft is estimated by combining SAGE and/or POAM data with high resolution, fast response in-situ ozone measurements acquired during the DC-8 ascent at the start of each science flight.

Livingston, J.↗

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗

Multi-Scale Mass Transfer Processes Controlling Natural Attenuation and Engineered Remediation: An IFRC Focused on Hanford’s 300 Area Uranium Plume

The Integrated Field Research Challenge (IFRC) at the Hanford Site 300 Area uranium (U) plume addresses multi-scale mass transfer processes in a complex subsurface biogeochemical setting where groundwater and riverwater interact. A series of forefront science questions on reactive mass transfer motivates research. These questions relate to the effect of spatial heterogeneities; the importance of scale; coupled interactions between biogeochemical, hydrologic, and mass transfer processes; and measurements and approaches needed to characterize and model a mass-transfer dominated biogeochemical system. The project was initiated in February 2007, with CY 2007, CY 2008, CY 2009, and CY 2010 progress summarized in preceding reports. A project peer review was held in March 2010, and the IFRC project acted upon all suggestions and recommendations made in consequence by reviewers and SBR/DOE. These responses have included the development of 'Modeling' and 'Well-Field Mitigation' plans that are now posted on the Hanford IFRC web-site, and modifications to the IFRC well-field completed in CY 2011. The site has 35 instrumented wells, and an extensive monitoring system. It includes a deep borehole for microbiologic and biogeochemical research that sampled the entire thickness of the unconfined 300 A aquifer. Significant, impactful progress has been made in CY 2011 including: (i) well modifications to eliminate well-bore flows, (ii) hydrologic testing of the modified well-field and upper aquifer, (iii) geophysical monitoring of winter precipitation infiltration through the U-contaminated vadose zone and spring river water intrusion to the IFRC, (iv) injection experimentation to probe the lower vadose zone and to evaluate the transport behavior of high U concentrations, (v) extended passive monitoring during the period of water table rise and fall, and (vi) collaborative down-hole experimentation with the PNNL SFA on the biogeochemistry of the 300 A Hanford-Ringold contact and the underlying redox transition zone. The modified well-field has functioned superbly without any evidence for well-bore flows. Beyond these experimental efforts, our site-wide reactive transport models (PFLOTRAN and eSTOMP) have been updated to include site geostatistical models of both hydrologic properties and adsorbed U distribution; and new hydrologic characterization measurements of the upper aquifer. These increasingly robust models are being used to simulate past and recent U desorption-adsorption experiments performed under different hydrologic conditions, and heuristic modeling to understand the complex functioning of the smear zone. We continued efforts to assimilate geophysical logging and 3D ERT characterization data into our site wide geophysical model, with significant and positive progress in 2011 that will enable publication in 2012. Our increasingly comprehensive field experimental results and robust reactive transport simulators, along with the field and laboratory characterization, are leading to a new conceptual model of U(VI) flow and transport in the IFRC footprint and the 300 Area in general, and insights on the microbiological community and associated biogeochemical processes influencing N, S, C, Mn, and Fe. Collectively these findings and higher scale models are providing a unique and unparalleled system-scale understanding of the biogeochemical function of the groundwater-river interaction zone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of Ex Vivo Ionizing Radiation on Static and Fatigue Properties of Mouse Vertebral Bodies

For a variety of medical and scientific reasons, human bones can be exposed to a wide range of ionizing radiation levels. In vivo radiation therapy (0.05 kGy) is used in cancer treatment, and ex vivo irradiation (25-35 kGy) is used to sterilize bone allografts. Ionizing radiation in these applications has been shown to increase risk of fracture, decrease bone quality and degrade collagen integrity. Past studies have investigated the deleterious effects of radiation on cortical or trabecular bone specimens individually, but to date no studies have examined whole bones containing both cortical and trabecular tissue. Furthermore, a clear relationship between the dose and the mechanical and biochemical response of bone's extracellular matrix has yet to be established for doses ranging from cancer therapy to allograft sterilization (0.05-35 kGy). To gain insight into these issues, we conducted an ex vivo radiation study to investigate non-cellular (i.e. matrix) effects of ionizing radiation dose on vertebral whole bone mechanical properties, over a range of radiation doses (0.05-35 kGy), with a focus on any radiation-induced changes in collagen. With underlying mechanisms of action in mind, we hypothesized that any induced reductions in mechanical properties would be associated with changes in collagen integrity. METHODS: 20-week old female mice were euthanized and the lumbar spine was dissected using IACUC approved protocols. The lumbar vertebrae (L1- S1) were extracted from the spine via cuts through adjacent intervertebral discs, and the endplates, posterior processes, surrounding musculature, and soft tissues were removed (approx. 1.5mm diameter, approx. 2mm height). Specimens were randomly assigned to one of five groups for ex vivo radiation exposure: x-ray irradiation at 0.05, 1, 17, or 35 kGy, or a 0 kGy control. Following irradiation, the vertebrae were imaged using microcomputed tomography (micro-CT) and then subjected to either monotonic compressive loading to failure or uniform cyclic compressive loading. During cyclic testing, samples were loaded in force control to a force level that corresponded to a strain of 0.46%, as determined in advance by a linearly elastic micro-CT-based finite element analysis for each specimen. Tests were stopped at imminent fracture, defined as a rapid increase in strain. The main outcome for the monotonic test was the strength (maximum force); for cyclic testing it was the fatigue life (log of the number of cycles of loading at imminent failure). A fluorometric assay was used on the S1 vertebrae to measure the number of non-enzymatic collagen crosslinks[4]. A one-way ANOVA was performed on mechanical properties and collagen crosslinks; means were compared with controls using Dunnett's method, with a Tukey-Kramer post-hoc analysis when significance was found (p < or = 0.05). RESULTS: Compared to the unirradiated control group, the concentration of non-enzymatic collagen crosslinks was significantly increased for all irradiated groups (p < 0.0001), and being higher by at least 50% (Figure 1a). By contrast, the radiation effects on the collagen were only evident at the higher doses. For irradiation exposures of 17 kGy or more, strength decreased substantially as the radiation level was increased, but no effect was evident below 17 kGy (Figure 1b). There was no significant change in the stiffness or maximum displacement for any radiation dose (p>0.05). The finite element analysis prescribed force level for cyclic loading exceeded the measured (monotonic) strength of the 17 and 35 kGy irradiated groups (mean +/- SD, 20.6 +/- 5.6 N; 13.2 +/- 3.7 N, respectively) and therefore these groups were eliminated from the fatigue study. The fatigue life for the 0.05 and 1 kGy groups were similar to each other and were not statistically significantly different from the control group (Figure 1c).

radiation↗

Laser Optogalvanic Spectroscopy pf Neon and Argon in a Discharge Plasma and its Significance for Microgravity Combustion

A detailed study of combustion mechanisms in flames, employing laser-based diagnostics, has provided good knowledge and understanding of the physical phenomena, and led to better characterization of the dynamical and chemical combustion processes, both under low-gravity (in space) and normal gravity (in ground based facilities, e.g. drop towers). Laser induced fluorescence (LIF), laser-induced incandescence (LII) and LIF thermometry have been widely used to perform nonintrusive measurements and to better understand combustion phenomena. Laser optogalvanic (LOG) spectroscopy has well-established applications in ion mobility measurements, atomic and molecular spectroscopy, ionization rates, recombination rates, velocity measurements and as a combustion probe for trace element detection. Absorption spectra of atomic and molecular species in flames can be obtained via LOG spectroscopy by measuring the voltage and current changes induced by laser irradiation. There are different kinds of processes that contribute to a discharge current, namely: (1) electron impact ionization, (2) collisions among the excited atoms of the discharge species and (3) Penning ionization. In general, at higher discharge currents, the mechanism of electron impact ionization dominates over Penning ionization, whereby the latter is hardly noticeable. In a plasma, whenever the wavelength of a laser coincides with the absorption of an atomic or molecular species, the rate of ionization of the species momentarily increases or decreases due to laser-assisted acceleration of collisional ionization. Such a rate of change in the ionization is monitored as a variation in the transient current by inserting a high voltage electrode into the plasma. Optogalvanic spectroscopy in discharges has been useful for characterizing laser line-widths and for providing convenient calibration lines for tunable dye lasers in the ultraviolet, visible and infrared wavelength regions. Different kinds of quantitative information, such as the electron collisional ionization rate, can be extracted from the complex processes occurring within the discharge. In the optogalvanic effect (OGE), there is no problem of overlap from background emissions, and hence even weak signals can be detected with a high signal-to-noise ratio, which makes the optogalvanic effect sensitive enough to resolve vibrational changes in molecular bonds and differences in energy levels brought about by different electron spins. For calibration purposes, neon and argon gaseous discharges have been employed most extensively, because these gases are commonly used as buffer gases within hollow-cathode lamps and provide an acceptable density of calibration lines. In the present work, our main aim has been to understand the dominant physical processes responsible for the production of the OGE signal, based on the extensive time resolved optogalvanic waveforms recorded, and also to extract quantitative information on the rates of excited state collisional processes.

Misra, Prabhakar↗

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS↗

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗

Architectural Analysis of Complex Evolving Systems of Systems

The goal of this collaborative project between FC-MD, APL, and GSFC and supported by NASA IV&V Software Assurance Research Program (SARP), was to develop a tool, Dynamic SAVE, or Dyn-SAVE for short, for analyzing architectures of systems of systems. The project team was comprised of the principal investigator (PI) from FC-MD and four other FC-MD scientists (part time) and several FC-MD students (full time), as well as, two APL software architects (part time), and one NASA POC (part time). The PI and FC-MD scientists together with APL architects were responsible for requirements analysis, and for applying and evaluating the Dyn-SAVE tool and method. The PI and a group of FC-MD scientists were responsible for improving the method and conducting outreach activities, while another group of FC-MD scientists were responsible for development and improvement of the tool. Oversight and reporting was conducted by the PI and NASA POC. The project team produced many results including several prototypes of the Dyn-SAVE tool and method, several case studies documenting how the tool and method was applied to APL s software systems, and several published papers in highly respected conferences and journals. Dyn-SAVE as developed and enhanced throughout this research period, is a software tool intended for software developers and architects, software integration testers, and persons who need to analyze software systems from the point of view of how it communicates with other systems. Using the tool, the user specifies the planned communication behavior of the system modeled as a sequence diagram. The user then captures and imports the actual communication behavior of the system, which is then converted and visualized as a sequence diagram by Dyn-SAVE. After mapping the planned to the actual and specifying parameter and timing constraints, Dyn-SAVE detects and highlights deviations between the planned and the actual behavior. Requirements based on the need to analyze two inter-system communication protocols that are representative of protocols used in the Aerospace industry have been specified. The protocols are related: APL s Common Ground System (CGS) as used in the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) and the Radiation Belt Space Probes (RBSP) missions. The analyzed communications were implementations of the Telemetry protocol and the CCSDS File Delivery Protocol (CFDP) protocol. Based on these requirements, three prototypes of Dyn-SAVE were developed and applied to these protocols. The application of Dyn-SAVE to these protocols resulted in the detection of several issues. Dyn-SAVE was also applied to several Testbeds that have previously been used for experimentation earlier on this project, as well as, to other protocols and logs for testing its broader applicability. For example, Dyn-SAVE was used to analyze 1) the communication pattern between a web browser and a web server, 2) the system log of a computer in order to detect offnominal computer shut-down behavior, and 3) the actual test cases of NASA Goddard s Core Flight System (CFS) and automatically generated test cases in order to determine the overlap between the two sets of test cases. In all cases, Dyn-SAVE assisted in providing insightful conclusions about each of the cases identified above.

Lindvall, Mikael↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

H-Canyon Jumper Gasket Sealing Evaluation

H-Canyon Documented Safety Analysis approved December 2019 credits the leak integrity of H-Canyon (HCA) process tanks, piping, and jumpers in preventing the release of radioactive material in a post seismic environment. Demonstration of leak tightness capability for this equipment was achieved through substantial analytical effort. However, one potential weak point that affects leak integrity exists at the gasketed joint between Hanford Connectors (HC) and nozzles. Review of the HCA Crane Log over a three-year period indicates the HCs periodically leak. Savannah River National Laboratory (SRNL) was requested to perform gasket compression tests that evaluate gasket relaxation characteristics, i.e., creep, and provide improvements to the HC installation process that can reduce or eliminate the frequency of leaks at HCs. A HC torque value of 312 ft-lbf is the target value used because it is the minimum torque value, with 95% confidence, when an unlubricated HC is impacted for standard sealing period of 9 seconds. In the HCA pipe-jumper connections are installed remotely by using HCs. HCs inherently add to the uncertainty of a gasketed joint because a calibrated torque value cannot be applied. The gasket material is a special material made of asbestos-fiber impregnated with Teflon. This material has been used for decades and generally performs well but does leak intermittently as gaskets creep and preloaded torque is lost. Gasket types are either snap-ring (SR) or lampshade (LS). Both gasket types cover the same seating area. The SR gasket is designed to remain with the jumper when the jumper is removed leaving the nozzle faces clean of any debris. SR gaskets are initially installed when the jumper is fabricated. Changing gaskets is time consuming and the radioactive environment presents risks to personnel. When jumpers are removed LS gaskets could be used to better seal the connection to reduce down time and ALARA concerns. That is, an LS gasket would be remotely positioned over the nozzle and the HC reinstalled. Testing also included two situations not typical of plant operation. One situation was the use of the non-standard double-layered gaskets (DG), and the other situation was the use of a LS without removal of the installed SR gasket. The DG was used when the supply of LS gaskets ended, so that testing could continue. It was of interest to demonstrate this gasket’s creep response because the DG, which was two gaskets glued together, was thicker than the LS gasket and the thicker material was expected to elicit a greater creep characteristics. The second situation of using two gaskets together by installing an LS gasket without removing the installed gasket is not allowed by the ASME Codes. However, by demonstrating the creep history of such a gasket arrangement it could be considered as an option that could possibly seal a connection while significantly reducing operational down time and personnel risk. Forty tests were conducted, and the overall results demonstrate the initial applied force and how the gasket changes with time. The change in sealing force, in terms of torque, was then measured to obtain a force history as the gaskets creep. This report documents the results of those tests as well as the method used to perform the tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Adaptive Changes In Postural Equilibrium And Motion Sickness Following Repeated Exposures To Virtual Environments

Virtual environments offer unique training opportunities, particularly for training astronauts and preadapting them to the novel sensory conditions of microgravity. Two unresolved human factors issues in virtual reality (VR) systems are: 1) potential "cybersickness", and 2) maladaptive sensorimotor performance following exposure to VR systems. Interestingly, these aftereffects are often quite similar to adaptive sensorimotor responses observed in astronauts during and/or following space flight. Changes in the environmental sensory stimulus conditions and the way we interact with the new stimuli may result in motion sickness, and perceptual, spatial orientation and sensorimotor disturbances. Initial interpretation of novel sensory information may be inappropriate and result in perceptual errors. Active exploratory behavior in a new environment, with resulting feedback and the formation of new associations between sensory inputs and response outputs, promotes appropriate perception and motor control in the new environment. Thus, people adapt to consistent, sustained alterations of sensory input such as those produced by microgravity, unilateral labyrinthectomy and experimentally produced stimulus rearrangements. Adaptation is revealed by aftereffects including perceptual disturbances and sensorimotor control disturbances. The purpose of the current study was to compare disturbances in postural control produced by dome and head-mounted virtual environment displays, and to examine the effects of exposure duration, and repeated exposures to VR systems. Forty-one subjects (21 men, 20 women) participated in the study with an age range of 21-49 years old. One training session was completed in order to achieve stable performance on the posture and VR tasks before participating in the experimental sessions. Three experimental sessions were performed each separated by one day. The subjects performed a navigation and pick and place task in either a dome or head-mounted display (HMD) VR system for either 30 or 60 min. The environment was a square room with 15 pedestals on two opposite walls. The objects appeared on one set of pedestals and the subject s objective was to move the objects to the other set of pedestals. After the subject picked up an object, a pathway appeared and they were required to follow the pathway to the other side of the room. The subject was instructed to perform the task as quickly and accurately as possible, avoiding hitting walls and other any obstacles and placing the object on the center of the pedestal. Postural equilibrium was measured (using the Equitest CDP balance system, Neurocom, International) before, immediately after, and at 1 hr, 2 hr, 4 hr and 6 hr following exposure to VR. Postural equilibrium was measured during quiet stance with eyes open, eyes closed and vision and/or ankle proprioceptive inputs selectively altered by servo-controlling the visual surround and/or support surface to the subject s center of mass sway. Posture data was normalized using a log transformation and motion sickness data were normalized using the square root. In general, we found that exposure to VR resulted in decrements in postural stability. The largest decrements were observed in the tests performed immediately following exposure to VR and showed a fairly rapid recovery across the remaining test sessions. In addition, subjects generally showed improvement across days. We found significant main effects for day and time for the composite equilibrium score and for sensory organization tests (SOT) 1, 2 and 6. Significant main effects were observed for day for SOT 3 and 5. Although we found no significant main effects for gender (when center of gravity was used as a covariate), we did observe significant gender X time interaction effects for composite equilibrium and for SOT 1, 3, 4 and 5. Women appeared to show larger decrements in postural stability immediately after exposure to VR than men, but recover more quickly than n. Finally, we found no significant main effects for type of VR device or for exposure duration, however, these factors did interact with other factors during some of the SOTs. Subjects exhibited rapid recovery of motion sickness symptoms across time following exposure to VR and significantly less severe symptoms across days. We did not observe main effects for gender, type of device or duration of exposure. Individuals recovered from the detrimental effects of exposure to virtual reality on postural control and motion sickness within one hour. Sickness severity and initial decrements in postural equilibrium decreases over days, which suggests that subjects become dual-adapted over time. These findings provide some direction for developing training schedules for VR users that facilitate adaptation, and support the idea that preflight training of astronauts may serve as useful countermeasure for the sensorimotor effects of space flight.

Harm, D. L.↗

Hydroacoustic Evaluation of Juvenile Salmonid Passage and Distribution at Lookout Point Dam, 2010

Pacific Northwest National Laboratory evaluated juvenile salmonid passage and distribution at Lookout Point Dam (LOP) on the Middle Fork Willamette River for the U.S. Army Corps of Engineers, Portland District (USACE), to provide data to support decisions on long-term measures to enhance downstream passage at LOP and others dams in USACE's Willamette Valley Project. This study was conducted in response to the listing of Upper Willamette River Spring Chinook salmon (Oncorhynchus tshawytscha) and Upper Willamette River steelhead (O. mykiss) as threatened under the Endangered Species Act. We conducted a hydroacoustic evaluation of juvenile salmonid passage and distribution at LOP during February 2010 through January 2011. Findings from this 1 year of study should be applied carefully because annual variation can be expected due to variability in adult salmon escapement, egg-to-fry and fry-to-smolt survival rates, reservoir rearing and predation, dam operations, and weather. Fish passage rates for smolt-size fish (> ~90 mm and < 300 mm) were highest during December-January and lowest in mid-summer through early fall. Passage peaks were also evident in early spring, early summer, and late fall. During the entire study period, an estimated total of 142,463 fish ± 4,444 (95% confidence interval) smolt-size fish passed through turbine penstock intakes. Of this total, 84% passed during December-January. Run timing for small-size fish (~65-90 mm) peaked (702 fish) on December 18. Diel periodicity of smolt-size fish showing crepuscular peaks was evident in fish passage into turbine penstock intakes. Relatively few fish passed into the Regulating Outlets (ROs) when they were open in summer (2 fish/d) and winter (8 fish/d). Overall, when the ROs were open, RO efficiency (RO passage divided by total project passage) was 0.004. In linear regression analyses, daily fish passage (turbines and ROs combined) for smolt-size fish was significantly related to project discharge (P<0.001). This relationship was positive, but there was no relationship between total project passage and forebay elevation (P=0.48) or forebay elevation delta, i.e., day-to-day change in forebay elevation (P=0.16). In multiple regression analyses, a relatively parsimonious model was selected that predicted the observed data well. The multiple regression model indicates a positive trend between expected daily fish passage and each of the three variables in the model-Julian day, log(discharge), and log(abs(forebay delta)); i.e., as any of the environmental variables increase, expected daily fish passage increases. For vertical distribution of fish at the face of the dam, fish were surface-oriented with 62%-80% occurring above 10 m deep. The highest percentage of fish (30%-60%) was found between 5-10-m-deep. During spring and summer, mean target strengths for the analysis periods ranged from -44.2 to -42.1 dB. These values are indicative of yearling-sized juvenile salmon. In contrast, mean target strengths in fall and winter were about -49.0 dB, which are representative of subyearling-sized fish. The high-resolution spatial and temporal data reported herein provide detailed information about vertical, horizontal, diel, daily, and seasonal fish passage rates and distributions at LOP from March 2010 through January 2011. This information will support management decisions on design and development of surface passage and collection devices to help restore Chinook salmon populations in the Middle Fork Willamette River watershed above LOP.

13 HYDRO ENERGY↗

Genesis Solar Wind – Capture, Return, Curate and Analyze: Looking Backward and Creating a Timeline

Introduction: In 1997 NASA’S Discovery Program selected the Genesis mission proposal to return solar wind samples to Earth for laboratory analyses. Principal Investigator Donald S. Burnett and the science team defined the purity of collector materials and ability to analyze solar wind composition to the precision required for planetary science. As a small mission, focused on a well-defined science goal, yet needing careful attention to engineering details, the communication among scientists and engineers, nurtured by Don Burnett, was exceptional. Genesis Mission and Curation Legacy: Genesis, as the first U. S. spacecraft to return astromaterial samples since Apollo, not only integrated the mission planning and flight teams, but also the science and sample curation teams during the mission development period. Since Genesis is a sample return mission, the Science Team was essential in certifying the collectors (sample containers for solar atoms). From inception, Genesis established mission funding for returned sample curation. JSC was lead in contamination control during mission preparation, including establishment of an ISO 4 cleanroom facility and use of ultrapure water (UPW) for cleaning flight hardware (and, as it turned out, for cleaning collectors after the mishap). Reliable, fast communication among scientists, engineers and curators at the hands-on level established deep respect among team members and efficient decision-making. JSC’s 50-years of astromaterial sample curation provided experienced sample processors onsite during recovery in Utah (a deep bench for emergency response). Post-recovery curation included iterative collaboration with science sample users to clean or verify cleanliness of samples. The science legacy from Genesis is addressed by Burnett and Jurewicz, this volume. In The Beginning: After Apollo sample return, Burnett and Marcia Neugebauer at JPL began discussing a solar wind sample return, with Neugebauer arguing that separate collection of solar wind regimes was essential science. By 1992 a solar wind sample return mission was presented at a workshop, and by 1994 a mission was proposed named Suess-Urey. The mission was re-proposed under a new name GENESIS and selected in 1997. Susan Niebur captured the Genesis mission history and stories, from high level management documents and from many interviews with participants [2]. Her account lets readers glimpse personality of participants in quotations from interviews. Need and Scope for Detailed Technical Timeline: A timeline constructed from lower level task documents has been initiated to document the resources and skills actually used, as well as task sequence or concurrency. Timelines for high level mission events are captured in two documents [1] [2] and for detailed re-entry events in [3]. A detailed technical timeline for Genesis mission and curation activities will provide data points for lower level tasks, such as ISO 4 curation facility construction time, preparation for nominal sample field recovery, mishap recovery, and UPW expansion. Changes in technology context 1990-2024: Semiconductor technologies were easily accessible in the U.S.A. (1990-1999), and the Genesis team used those resources for cleanroom design and UPW system expansion. Image documentation was changing from film to digital during cleanroom construction and payload cleaning (1997-2001). Engineering design was done using computer aided design proprietary software, making more difficult the archiving of payload configuration and materials. Email of documents, tracked delivery service and virtual meeting capability greatly improved communication efficiency. Information sources – Pre-launch mission preparation: Examples of mission science, engineering and contamination control are collector purity testing, payload design/fabrication and ISO 4 cleanroom construction. Information on timing of these activities comes from facility readiness reviews, management reviews, shipping documents, procurement documents, test reports, travel documents, laboratory logs, Quality Assurance documents, dates on images, participant notebooks and emails. Information sources – Sample return re-entry and field recovery activities: Information comes from event timelines produced by Mid-Air Recovery team, Lockheed team lead notes and from chase video, JPL Quality Assurance. Information also comes from images and logbooks from UTTR cleanroom operations and from curatorial documents. Information sources – Resulting science and sample cleaning processes: Agendas from the annual gatherings of the science team initially trace testing for collector purity/cleanliness, and after sample recovery, include collector cleaning and cleanliness assessment. Post-recovery documents include curatorial orders and procedures, sample allocation documents and LPSC abstracts. Timeline Objectives: A simple spreadsheet timeline with headers DATE, EVENT, PEOPLE, COMMENT, INFORMATION SOURCE has been initiated and currently has over 90 entries. While this is not definitive historical research, it is a quick look at the evolution of Genesis curation with pointers to documents or people with information. Engineers for future missions may find useful points of comparison for development of facilities. References:[1] Genesis Mission Reference Document, (2011) JPL D-62382.[2] Niebur S. M., edited by Brown D. W. (2023) NASA’s Discovery Program: The First 20 Years of Competitive Planetary Exploration, NASA-SP-2023-4238.[3] Genesis Mishap Investigation Board Report, Vol. 1 (July 2005).

solar wind↗

Analysis of Cavern and Well Stability at the West Hackberry SPR Site Using a Full-Dome Model

This report presents computational analyses that simulate the structural response of caverns at the Strategic Petroleum Reserve (SPR) West Hackberry site. The cavern field comprises 22 caverns. Five caverns (6, 7, 8, 9, 11) were acquired from industry and have unusual shapes and a history dating back to 1946. The other 17 caverns (101-117) were leached according to SPR standards in the mid-1980s and have tall cylindrical shapes. The history of the caverns and their shapes are simulated in a three-dimensional geomechanics model of the site that predicts deformations, strains, and stresses. Future leaching scenarios corresponding to oil drawdowns using fresh water are also simulated by increasing the volume of the caverns. Cavern pressures are varied in the model to capture operational practices in the field. The results of the finite element model are interpreted to provide information on the current and future status of subsidence, well integrity, and cavern stability. The most significant results in this report are relevant to Cavern 6. The cavern is shaped like a bowl with a large ceiling span and is in close proximity to Cavern 9. The analyses predict tensile stresses at the edge of the ceiling during repressurization of Cavern 6 following workover conditions. During a workover the cavern is at low pressure to service a well. The wellhead pressures are atmospheric. When the workover is complete, the cavern is repressurized. The resulting elastic stresses are sufficient to cause tension around the edge of the large ceiling span. With time, these stresses relax to a compressive state because of salt creep. However, the potential for salt fracture and propagation exists, particularly towards Cavern 9. With only 200 feet of salt between the caverns, the operational consequences must be examined if the two caverns become connected. A critical time may be during a workover of Cavern 9 in part because of the operational vulnerabilities, but also because dilatant damage is predicted under the ledge that forms the lower lobe in the cavern. The remaining caverns have no significant issues regarding cavern stability and may be safely enlarged during subsequent oil drawdowns. Predicted well strains and subsidence are significant and consequently future remedial actions may be necessary. These predicted well strains certainly suggest appropriate monitoring through a well-logging program. Subsidence is currently being monitored.

02 PETROLEUM↗

Assessment of CO 2 Storage Resources in Depleted Oil and Gas Fields in the Ship Shoal Area, Gulf of Mexico

The Gulf of Mexico is one of the most important regions in the United States for energy resources and infrastructure. Gulf of Mexico federal offshore oil and gas production accounts for 17% of total U.S. crude oil production and 5% of total U.S. gas production (EIA Gulf of Mexico Fact Sheet). This region presents an excellent combination of high need and significant opportunity for large scale geologic storage of CO 2 . The Ship Shoal Area is located about 20 miles offshore Louisiana within the Gulf Coast federal waters. Miocene and Pliocene sediments in the Ship Shoal Area are proven to provide excellent and secure traps for oil and gas. The Ship Shoal Area contains a large number of depleted oil and gas fields either currently abandoned, or planned for abandonment by 2025, which may provide very significant potential CO 2 storage capacity. GeoMechanics Technologies conducted a comprehensive research project to better characterize Neogene sediments in the Ship Shoal Area for high volume CO 2 storage. The data generated from literature, well data, well log, and formation evaluation were input into Rockwork 16 geologic software to create a detailed geologic model for the Ship Shoal study area. The geologic grids were then fed into TOUGH2, the gas migration model software and FLAC 3D, the geomechanical model software. The fluid flow modeling for both the Miocene and Pliocene simulation from both Block 107 and Block 84 indicated that Pliocene and Miocene are a good reservoir for the CO 2 sequestration. Thirty years of CO 2 injection (at a rate of 1 million metric tons CO 2 per year) and 30 years of observation simulations were run with the fluid flow modeling. The 3D geomechanical modeling for Ship Shoal Block 107 and Block 84 were developed to evaluate induced surface deformation and potential fault reactivation after 30 years CO 2 injection at base of Pliocene and at upper Miocene injection locations. Evaluation results indicate low to no risks for fault slips or fault reactivation. Twelve well bores and 76 well bores for Ship Shoal Blocks 84 and 107 respectively were reviewed for their cement history. Most wells (57 out of 76 and 7 out of 12 wells in Ship Shoal Block 107 and 84 respectively) have good integrity. Nineteen and 5 wells (in Ship Shoal Block 107 and 84 respectively) with no top plug, incomplete cement or Plug and Abandonment information are given yellow cautionary indicators. These cautionary wells may provide leakage paths of CO 2 through the well bores to the USDWs. Using our Quantitative Risk & Decision Analysis Tool (QRDAT) for caprock integrity evaluation, we compared Ship Shoal’s risk to that of In Salah, Sleipner, Kevin Dome, Loudon, Illinois Industrial CCS and Wilmington Graben. We found the risk at the Ship Shoal Blocks 84 and 107 fields are similar to the known CO 2 active sequestration sites, but lower than the Wilmington Graben turbidities offshore California studied site. The risk of natural seismicity in the Gulf of Mexico is relatively low. Geomechanics Technologies has documented the top 25 CO 2 emission sources within the close proximity of the Ship Shoal Block 84 and 107 fields. All the offshore and onshore pipelines have been digitized and can be viewed in an interactive website (http://www.geomechanicstech.com/shipshoal.html). We also performed a feasibility study on the potential for converting existing oil or gas pipelines for CO 2 transport. There are abandoned, idled and retired onshore and offshore pipelines. However, there is no standard specification for maximum pipeline pressures needed for CO 2 transport. It is a function of design, materials, and testing provided; thus it will be the responsibility of the pipeline operator to correctly determine, maintain and operate within the limits of the pipelines. There are a few transit corridors extending from onshore Louisiana (Cailou Bay) to offshore trunk-lines. The cost for constructing a pipeline has increased about 46% from 2015 to an average cost of $5,064,046 per mile in 2016. We use the NETL CO 2 storage resource mass estimate to calculate the potential resources for Ship Shoal Block 84 and 107 fields using the depleted oil and gas reservoir volume and sand volume generated from our geologic modeling. The estimated storage resource results are greater when using the sand volumes derived through geologic modeling than the oil and gas reservoir data. The difference is due to the depleted oil and gas reservoir information not accounting for the water-flooded sand located below the oil/gas-water contact or unproductive sand units. Resource calculation using the sand volume obtained through geologic modeling overestimate the storage capacity as the model accounts for all sand within the formation, not just the interconnected sand.

02 PETROLEUM↗

Reservoir Characterization of Bridgeport and Cypress Sandstones in Lawrence Field Illinois to Improve Petroleum Recovery by Alkaline-Surfactant-Polymer Flood

Within the Illinois Basin, most of the oilfields are mature and have been extensively waterflooded with water cuts that range up to 99% in many of the larger fields. In order to maximize production of significant remaining mobile oil from these fields, new recovery techniques need to be researched and applied. The purpose of this project was to conduct reservoir characterization studies supporting Alkaline-Surfactant-Polymer Floods in two distinct sandstone reservoirs in Lawrence Field, Lawrence County, Illinois. A project using alkaline-surfactantpolymer (ASP) has been established in the century old Lawrence Field in southeastern Illinois where original oil in place (OOIP) is estimated at over a billion barrels and 400 million barrels have been recovered leaving more than 600 million barrels as an EOR target. Radial core flood analysis using core from the field demonstrated recoveries greater than 20% of OOIP. While the lab results are likely optimistic to actual field performance, the ASP tests indicate that substantial reserves could be recovered even if the field results are 5 to 10% of OOIP. Reservoir characterization is a key factor in the success of any EOR application. Reservoirs within the Illinois Basin are frequently characterized as being highly compartmentalized resulting in multiple flow unit configurations. The research conducted on Lawrence Field focused on characteristics that define reservoir compartmentalization in order to delineate preferred target areas so that the chemical flood can be designed and implemented for the greatest recovery potential. Along with traditional facies mapping, core analyses and petrographic analyses, conceptual geological models were constructed and used to develop 3D geocellular models, a valuable tool for visualizing reservoir architecture and also a prerequisite for reservoir simulation modeling. Cores were described and potential permeability barriers were correlated using geophysical logs. Petrographic analyses were used to better understand porosity and permeability trends in the region and to characterize barriers and define flow units. Diagenetic alterations that impact porosity and permeability include development of quartz overgrowths, sutured quartz grains, dissolution of feldspar grains, formation of clay mineral coatings on grains, and calcite cementation. Many of these alterations are controlled by facies. Mapping efforts identified distinct flow units in the northern part of the field showing that the Pennsylvanian Bridgeport consists of a series of thick incised channel fill sequences. The sandstones are about 75-150 feet thick and typically consist of medium grained and poorly sorted fluvial to distributary channel fill deposits at the base. The sandstones become indistinctly bedded distributary channel deposits in the main part of the reservoir before fining upwards and becoming more tidally influenced near their top. These channel deposits have core permeabilities ranging from 20 md to well over 1000 md. The tidally influenced deposits are more compartmentalized compared to the thicker and more continuous basal fluvial deposits. Fine grained sandstones that are laterally equivalent to the thicker channel type deposits have permeabilities rarely reaching above 250 md. Most of the unrecovered oil in Lawrence Field is contained in Pennsylvanian Age Bridgeport sandstones and Mississippian Age Cypress sandstones. These reservoirs are highly complex and compartmentalized. Detailed reservoir characterization including the development of 3-D geologic and geocellular models of target areas in the field were completed to identify areas with the best potential to recover remaining reserves including unswept and by-passed oil. This project consisted of tasks designed to compile, interpret, and analyze the data required to conduct reservoir characterization for the Bridgeport and Cypress sandstones in pilot areas in anticipation of expanded implementation of ASP flooding in Lawrence Field. Geologic and geocellular modeling needed for reservoir characterization and reservoir simulation were completed as prerequisites for design of efficient ASP flood patterns. Characterizing the complex reservoir geology that identifies the geologic conditions that will optimize oil recoveries for expansion of the ASP pilots in the Bridgeport and Cypress sandstones to other areas of Lawrence Field is the primary objective of this project. It will permit evaluation of efficiency of oil recovery from Bridgeport and Cypress sandstone reservoirs using ASP technology. Additionally evaluation of similar Pennsylvanian and Chesterian reservoirs shows that it is likely that ASP flood technology can be successfully applied to similar reservoirs in the Illinois Basin as well as to other U.S. reservoirs. Chemical flooding was introduced in stages with the first flood initiated in 2010 and a second offset pilot project initiated during 2011. Rex Energy Corporation is reporting a positive response on its ASP Middagh pilot project in the Pennsylvanian Bridgeport B reservoir, Lawrence Field. Oil response in the 15 acre flood has continued to show an increase in oil cut from 1% to 12%. Total pattern production increased from 16 BOPD and stabilized at a range of 65-75 BOPD in the last three months of 2011. Peak production rose to 100 + BOPD. Oil cut in the pilot increased for 1.0% to ~ 12.0% with an individual well showing oil cuts greater than 20%. A second, 58 acre pilot (Perkins-Smith) adjacent to and likely in communication with the Middagh pilot has been initiated. Preliminary brine injection has been implemented and ASP injection was initiated in mid-2012. Response is expected by mid-2013 with peak recovery expected by late 2013. Rex Energy is projecting full scale expansion with the next step of development being a 351 acre project scheduled to begin in mid-2013. Preliminary development has been initiated in this Delta Unit area located in the south half of section 32, T4N, R12W.

02 PETROLEUM↗

Oilwell Conversion (Well API 121913310501) to Geothermal Heat Storage Well for Flexible Electricity Storage

Geothermal growth is limited by a lack of geographically dispersed high-temperature thermal resources and high initial upfront investment in characterization and well construction. This project intended to address the challenges of energy supply intermittency and enhance grid resilience, reliability, and energy security by storing energy provided from currently available renewable resources in the subsurface to harvest it a later time during at-peak energy demand. This project intended to improve geothermal adoption, reduce initial project risk, and improve price competitiveness through utilizing existing oil and gas infrastructure such as non-productive wells, non-economic fields, dry holes, and orphaned wells. The project also intended to address the lack of geographically dispersed thermal resources and enhance grid resilience, reliability, and energy security by introducing an economical method for storing energy from currently available renewable resources in the subsurface for usage during at-peak energy demand. During this research, the project furthered the understanding of the feasibility of utilizing abandoned oil and gas wells as geothermal heat storage wells. The project team investigated the heat storage and hydrogeological characteristics of subsurface reservoirs in the Illinois Basin to evaluate their response to heat injection for determining the evolution of temperature profiles and heat losses over time using existing and available data sets. The project team then performed modeling and simulation to evaluate the heat losses of returning fluids during heat extraction. The outputs were used to select an optimal candidate reservoir and location in Southern Illinois. The team designed and performed a small-scale field test in an existing oil well to refine the model and to demonstrate the permitting and regulatory pathways necessary for the conversion of oil and gas assets to geothermal use. The field test also serves as a proof of concept and can guide the procedures for future research and implementation. Additionally, the project team, conducted initial market research and customer discovery to develop a go to market strategy for an Advanced Geothermal Energy Storage (AGES) system. The project team in this research also identified the parameters to be refined in future research, to improve the current go to market strategy economic model. To this end several subject matter experts were also identified to assist in future research with geothermal infrastructure setup, energy storage policy and law, energy storage market demand, potential siting based on demand etc. Future research will involve further sophistication of the site commercial modeling, implementing a larger-scale test, and further refinement of the thermodynamic modeling/simulation process. The output will be lifecycle costs and economics suitable for comparison to alternative approaches from a validated full-scale demonstration for venture capital investment into this technology. The project successfully demonstrated the ability to leverage existing oilfield infrastructure, permits, and land access and leasing agreements, to enable geothermal storage projects to come online faster and cheaper than a greenfield development could. This technology could allow for greater energy independence and security through long-term energy storage solutions. The longer duration allows for greater storage for renewables currently limited by hours-long storage durations of lithium-ion. The AGES system would support the growth of renewable energy farms, and provide greater opportunities for a cleaner energy infrastructure.

15 GEOTHERMAL ENERGY↗

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↗