SEARCH · Engineering Papers
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.
Co-Author Voon Hui Lai, Australian National University, will be attending IUGG and presenting the work. The session of interested is "S05 - Advances in Earthquake and Explosion Monitoring Using Distributed Acoustic Sensing"
We examine the ability of three different inversion methods: 1) first motion (FM) inversions, 2) amplitude inversions, and 3) full waveform (FW) time variable moment tensor (TVMT) inversions, to recover an accurate source mechanism for simulated data for an earthquake as well as an explosion. The ability of inversion methods, such as the ones described above, to recover accurate models representative of the data depends on both a priori information as well as the quality of data being inverted. Therefore, we examine the effect that station geometry, geologic model, and noise level has on the inversion and estimated source mechanism. We find that FM data can provide more robust solutions than amplitude data and are not as sensitive to inaccurate earth models, especially in low-noise cases, and overall FW inversions are the most accurate out of the methods examined, but can still be biased by inaccurate earth models.
As part of the Cyclotron Road program, METAseismic sought to investigate a multi-scale approach, both in terms of manufacturing and design, for seismic metamaterials with applications in the energy sector. Specifically, this project explored the feasibility of using nanofabrication to enable the incorporation of 3D nanostructured components into a mechanical metamaterial. These components were explored to increase the seismic resilience of our infrastructure to earthquakes, with the final aim of providing customized levels of protection for its components. This was envisioned as the base of a multi-scale design approach that consisted of using a multiplicity of metamaterial architectures ranging from the nano- to the meso-scale to address the stringent seismic design objectives of electric-power equipment.
During natural hazard events (hurricanes, wildfires, earthquakes, etc.) and recent man-made events (e.g., cyber attacks), the exchange of near real-time, spatially refined data within the response community is critical. The EAGLE-I$^{TM}$ platform is one tool that facilitates this data for decision makers within the energy sector. While much information can be collected and integrated into the system directly, other pertinent data must be augmented by other derived data products to enhance the information and allow for a consistent evaluation of on-the-ground conditions. One such data set that requires the addition of other derived data is the electric utility customer outage data that is aggregated to the county level within the EAGLE-I application. Without a county customer data set, outages can only be compared on total counts, which gives greater importance to higher population outages. Including an electric utility customer data set at the county level allows for these outage counts to be converted to percent outages and brings a consistent classification of outages and equal importance to all outages. To achieve this, several available data sets were combined and spatial disaggregation techniques were employed to model customer estimates at the county scale. This paper presents the approach to produce this data for the United States and lessons learned from working with these disparate data sets. Data validation is provided, where possible, and limitations of the model and possible improvements are discussed.
We present a preliminary version of the East Asia Tomography (EAT) model, an adjoint waveform tomography model of East and Southeast Asia. We used SPiRaL (Simmons et al., 2021) as our starting model and source parameters for 250 earthquakes from the Global Centroid Moment Tensor catalogue (Ekström et al., 2012). Over 198 iterations, we have iterated the EAT model down to a minimum period of 35 seconds. We plan on continuing our iteration technique down to 30 seconds period before updating our misfit function to use a normalized cross correlation-based misfit functions (e.g., Tao et al., 2018) to better constrain Earth structure. We hope to iterate the current extent of the model to 25 seconds minimum period before iterating to shorter periods for a smaller subregion of the full model.
Fifteen years have passed since the March 11, 2011, the Great East Japan Earthquake (GEJE) triggered a tsunami that inundated the Fukushima Daiichi Nuclear Power Station (Daiichi or 1F) and led to severe accidents in multiple units at the plant. During these fifteen years since these events, Daiichi forensics examinations and evaluations have revealed many insights for enhancing the safety of nuclear power plants worldwide. However, uncertainties remain regarding the events after the GEJE. Additional Daiichi forensics information is required to inform future Decontamination and Decommissioning (D&D) activi- ties conducted by the plant owner, Tokyo Electric Power Company Holdings, Incorporated (TEPCO Hold- ings). This forensics information also has important implications for the safety and operation of U.S. commercial nuclear power plants. A collaborative U.S. and Japanese effort, which includes approximately 50 experts in reactor safety and plant operations, was initiated in 2014 by the Department of Energy Office of Nuclear Energy (U.S. DOE-NE) to identify Daiichi examination needs and evaluate recent Daiichi examination data to address these needs. This collaborative effort, which is beneficial to the U.S. and Japan, continues as Japan moves forward with D&D activities.
The U.S. Department of Energy’s (DOE) Energy Transitions Initiative Partnership Project (ETIPP) connects remote and island communities, regional partners, and the DOE national laboratories to support communities as they seek to build resilience in their energy systems. The Makah Tribe faces several energy challenges, including frequent power outages and the potential for an extended outage due to an earthquake or tsunami. The Tribe joined ETIPP in 2022 to address those challenges, seeking to build energy resilience and sovereignty in the community. The Makah Tribe, Spark Northwest, the Pacific Northwest National Laboratory (PNNL), and the National Renewable Energy Laboratory (NREL) collaborated to develop a strategic energy plan as part of the second cohort of ETIPP communities. The long-term energy vision of the Makah Tribe includes increasing energy efficiency in the community, improving energy management capacity, and developing the renewable energy generation and storage sufficient to independently power the Reservation for one year. Additionally, the ETIPP team worked with Makah leadership, staff, and community members to identify a set of community priorities, values, and goals to guide energy development as the Tribe takes the incremental steps toward their vision for energy sovereignty. Those energy values include ecosystem-based management, energy sovereignty and project ownership, workforce development and capacity, economic opportunity, community wellbeing and priorities, and emergency disaster resilience. To understand what would be needed for a year for energy independence, the PNNL team conducted an assessment to determine the current energy usage of the Tribe and also modeled several scenarios for future energy use. Using the energy usage values, the team estimated how two types of renewable energy technology, specifically locally deployed solar and small-scale wind, could contribute towards the energy independence goal. The energy baseline and resource assessment produced the following key findings:
Energy Resilience Analysis for Distribution Power System (ERAD) is a free, open-source Python toolkit for estimating the energy and service impacts of hazards like earthquakes and flooding. It uses a graph-based approach to capture high resolution connectivity among the grid, critical services, and customers and rapidly compute household level metrics and aggregated statistics across large distribution systems. It uses asset fragility curves that relate hazard severity to survival probability for power system equipment including cables, transformers, substations, etc. The tool is designed to be modular and extensible, allowing it to interface with third-party hazard simulators and integrate into broader resilience analysis workflows. ERAD enables researchers, students, communities, distribution utilities, and other stakeholders to understand hazard impacts and evaluate the effectiveness of different programs to improve energy resilience. The webinar was hosted by NLR researcher Aadil Latif.
A dramatic increase in frequency of minor to major flooding since 2000 has caused significant economic losses across the world. To mitigate and recover from these losses, actions have been taken to build resilient communities and infrastructures, specifically, by providing situational awareness in near real-time about flood impacts to enhance response and recovery efforts. Several hydrologic and hydraulic flood models are available at various spatial and temporal resolutions to forecast flood events at regional to global scale. Given the global coverage of two operational flood models – GloFAS (Global Flood Awareness System) and GFMS (Global Flood Monitoring System), the purpose of this project is to implement a Model of Models (MoM) approach to integrate the outputs from these two models to classify flood severity at watershed level worldwide, and send alerts based on severity similar to the USGS PAGER (used for severity alerting and impact analysis for earthquakes) to flood impacted communities. The alerts containing flood impacts and severity information will be disseminated through the DisasterAWARE platform, operated by the Pacific Disaster Center (PDC), that provides global multi-hazard alerting and Situational Awareness information to the emergency management community and public. The current version of the MoM approach was implemented for a case study flood event that occurred during January and February of 2020 in South and Central Africa. The findings of the case study event reveal that the approach is effective in identifying potential flood impact areas and the spatio-temporal variation of flood severity, flood depth and extent at watershed level, which will be used to assess infrastructure and societal impacts using earth-observation data and for alerting.
Natural hazard-induced nuclear accidents, such as the Fukushima Daiichi Accident that occurred in Japan in 2011, have significantly increased reactor safety studies in understanding nuclear power plant (NPP) responses to external hazard events such as earthquakes and floods. Natural hazards could cause the loss of offsite power in nuclear power plants, potentially leading to a Station Blackout (SBO) accident that significantly contributes to the overall risk of nuclear power plant accidents. Despite the fact that extensive research has been conducted on the station blackout accident for nuclear power plant, further understanding of these events is needed, particularly in the context of the dynamic nature of external hazards such as external flooding. This paper estimates the progression of station blackout events for a generic pressurized water reactor (PWR) in response to external flooding events. The original RELAP5-3D model of the Westinghouse four-loop design pressurized water reactor was adopted and modified to simulate the external flood-induced station blackout accident, including the short-term and long-term station blackout scenarios. A sensitivity analysis of long-term station blackout, examining reactor operation times and analyzing key parameters over time, was also conducted in this work. The results of the analyses, especially the critical timing parameters of key event sequences, provide useful insights about the time during the external flooding event, which is important for plant operators to make timely decisions to prevent potential core damage. This paper represents significant progress toward developing an integrated risk assessment framework for further identifying and assessing the effects of the critical sources of uncertainties of nuclear power plant under external hazard-induced events.
The field of fracture mechanics was developed during the throes of World War II, and since then, it has been a very active area of research. This is a very energetic and vibrant field because fracture processes are relevant for a vast range of engineering applications spanning different temporal and length scales. Depending on the application, material fracture and fragmentation may be either a desirable feature (for example, for the purpose of hydraulic fracturing operations), or something that is to be avoided at all costs (for example, in key components of mechanical systems). In some other instances, fractures are already present and must be accommodated through novel design (for example, in the case of earthquake ruptures occurring in the Earth’s crust). Because of this wide spectra of motivation, there is a continuous need for the fracture mechanics research community to gain more insight on how fractures are created, evolve, and interact with themselves and with material microstructures. In addition, there are still many remaining questions about how fracture propagation, growth, and interaction impact overall material behavior, and how to develop accurate models to predict this behavior.
The use of uncrewed aerial systems (UAS) increases the opportunities for detecting surface changes in remote areas and in challenging terrain. Detecting surface topographic changes offers an important constraint for understanding earthquake damage, groundwater depletion, effects of mining, and other events. For these purposes, changes on the order of 5–10 cm are readily detected, but sometimes it is necessary to detect smaller changes. An example is the surface changes that result from underground explosions, which can be as small as 3 cm. Previous studies that described change detection methodologies were generally not aimed at detecting sub-5-cm changes. Additionally, studies focused on high-fidelity accuracy were either computationally modeled or did not fully provide the necessary examples to highlight the usability of these workflows. Detecting changes at this threshold may be critical in certain applications, such as global security research and monitoring for high-consequence natural hazards, including landslides. Here we provide a detailed description of the methodology we used to detect 2–3 cm changes in an important applied research setting—surface changes related to underground explosions. This methodology improves the accuracy of change detection data collection and analysis through the optimization of pre-field planning, surveying, flight operations, and post-processing the collected data, all of which are critical to obtaining the highest output data resolution possible. We applied this methodology to a field study location, collecting 1.4 Tb of images over the course of 30 flights, and location data for 239 ground control points (GCPs). We independently verified changes with orthoimagery, and found that structure-from-motion, software-reported root mean square errors (RMSEs) for both control and check points underestimated the actual error. We found that 3 cm changes are detectable with this methodology, thereby improving our knowledge of a rock’s response to underground explosions.
The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations.
Extreme, tsunami-like wave runup events in the absence of earthquakes or landslides have been attributed to trapped waves over shallow bathymetry, long waves created by atmospheric disturbances, and long waves generated by abrupt breaking. These runup events are associated with inland excursions of hundreds of meters and periods of minutes. While the theory of radiation stress implies that nearshore energy transfer from the carrier waves to the infragravity waves can also lead to very large runup, there have not been observations of runup events induced by this process with magnitudes and periods comparable to the other three mechanisms. This work presents observations of several runup events in the US Pacific Northwest that are comparable to extreme runup events related to the other three mechanisms. It also discusses possible generation mechanisms and shows that energy transfer from carrier waves to bound infragravity waves is a plausible generation mechanism. In addition, a method to predict and forecast extreme runup events with similar characteristics is presented.
Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.
Porosity reduction in rocks from a fault core can cause elevated pore fluid pressures and consequently influence the recurrence time of earthquakes. We investigated the porosity distribution in the New Zealand's Alpine Fault core in samples recovered during the first phase of the Deep Fault Drilling Project (DFDP-1B) by using two-dimensional nanoscale and three-dimensional microscale imaging. Synchrotron X-ray microtomography-derived analyses of open pore spaces show total microscale porosities in the range of 0.1 %–0.24 %. These pores have mainly non-spherical, elongated, flat shapes and show subtle bipolar orientation. Scanning and transmission electron microscopy reveal the samples' microstructural organization, where nanoscale pores ornament grain boundaries of the gouge material, especially clay minerals. Our data imply that (i) the porosity of the fault core is very small and not connected; (ii) the distribution of clay minerals controls the shape and orientation of the associated pores; (iii) porosity was reduced due to pressure solution processes; and (iv) mineral precipitation in fluid-filled pores can affect the mechanical behavior of the Alpine Fault by decreasing the already critically low total porosity of the fault core, causing elevated pore fluid pressures and/or introducing weak mineral phases, and thus lowering the overall fault frictional strength. We conclude that the current state of very low porosity in the Alpine Fault core is likely to play a key role in the initiation of the next fault rupture.
This paper presents a methodology for freight traffic assignment in a large-scale road-rail intermodal network under uncertainty. Network uncertainties caused by natural disasters have dramatically increased in recent years. Several of these disasters (e.g., Hurricane Sandy, Mississippi River Flooding, and Hurricane Harvey) severely disrupted the U.S. freight transportation network, and consequently, the supply chain. To account for these network uncertainties, a stochastic freight traffic assignment model is formulated. An algorithmic framework, involving the sample average approximation and gradient projection algorithm, is proposed to solve this challenging problem. The developed methodology is tested on the U.S. intermodal network with freight flow data from the Freight Analysis Framework. The experiments consider three types of natural disasters that have different risks and impacts on transportation networks: earthquakes, hurricanes, and floods. It is found that for all disaster scenarios, freight ton-miles are higher compared to the base case without uncertainty. The increase in freight ton-miles is the highest under the flooding scenario; this is because there are more states in the flood-risk areas, and they are scattered throughout the U.S.
Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.