Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulated flooding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Heavy Oil Polymer EOR in the Challenging Alaskan Arctic - It Works!

Under co-sponsorship of the US Department of Energy and Hilcorp Alaska LLC the first ever polymer field pilot commenced on 8/28/2018 in the Schrader Bluff heavy oil reservoir at the Milne Point Field on Alaska North Slope (ANS). The primary objective of the pilot is to prove the efficacy of polymer Enhanced Oil Recovery (EOR) to unlock the vast heavy oil resources on ANS. More than two and half years after startup, the polymer injection, supporting laboratory experiments and simulation studies steadily continue. The pilot started injecting hydrolyzed polyacrylamide (HPAM), at a concentration of 1,750 ppm to achieve a target viscosity of 45 cP, into the two horizontal injectors in the flood pattern. Production is monitored in the two horizontal producers. Based on laboratory measurements of polymer viscosity at reservoir conditions, the team decided to reduce polymer concentration to 1,200 ppm since July 2020 in an effort to control injection pressure and optimize polymer utilization. Quality control (QC) on the field ensures uniform polymer solution properties. Representative rock and fluid systems and test conditions are utilized in the corefloods on polymer retention, effect of injection water salinity, polymer loading, and their combinations on oil recovery. A history matched reservoir simulation model for forecasting oil recovery was developed on the basis of all the available field data. Field concerns related to the post-polymer breakthrough impact on flow assurance is addressed via specialized laboratory tests. Notwithstanding early operational disruptions and hydration issues, continuous polymer injection in both injectors has been achieved. To date, 950,000 lbs of polymer or 2 million barrels of polymer solution, equating to 13% of total pore volume (PV), has been placed in the flood pattern, serving as an effective indicator of adequate polymer injectivity. So far, the success of polymer EOR is evident from drastically reduced water cut in the producers, an estimated incremental 1,000 bopd over waterflood, and a favorable polymer utilization of 1.7 lbs/barrel of incremental oil. Polymer breakthrough was observed 26 months after the start of polymer injection. Main observations from corefloods are unusually high polymer retention values in some cases and a positive response to low salinity water. Although the heterogeneity in the flood pattern and exceptionally low water cut pose some challenges, persistent novel and justifiable simulation approaches have resulted in a robust history matched model. Experimental results on produced fluids treatment provide operational guidance for improved oil-water separation and mitigation of heater tube fouling. The ongoing success of the pilot is a key indicator of bringing the team ever closer to meeting the project’s primary objective. The encouraging results of the pilot is one of the drivers that has provided the impetus to apply polymer EOR throughout the Milne Point Field, which would increase oil recovery and extend the economic life of the Trans Alaska Pipeline System. The scientific knowledge, including the many lessons learned from this pilot also has referential value for other potential heavy oil EOR projects throughout the world. The collected data, operational lessons learned, and the overall success of the pilot are summarized in the paper.

None, None↗
DOE OSTI2021

Maximizing oil production from water alternating gas (CO 2 ) injection into residual oil zones: The impact of oil saturation and heterogeneity

Residual oil zones (ROZs) are widespread reservoirs, characterized by oil at residual saturation, either underlying oil fields (brownfield) or lateral (greenfield) to such fields. These reservoirs have the potential to produce volumes of oil sufficiently significant to make appreciable impacts on the US’s oil reserves and associated incidental CO 2 sequestration. The objective of this study is to improve our understanding the impact of heterogeneous and low oil saturations, in brownfield ROZs, on the effectiveness of water alternating gas (WAG) injection strategies. ROZs occur in the Permian Basin and elsewhere, and operators are using CO 2 injection for enhanced oil recovery (EOR) in these zones. The consensus model for the formation of ROZs is that they were formed by the effect of faster regional aquifer flow, acting over millions of years. Both the magnitude of oil saturation and the spatial distribution of oil differ from water-flooded main pay zones (MPZs). To explore the most effective injection strategies, we conducted simulations of CO 2 injection into synthetic geologic reservoirs. These simulations focused on injection into reservoirs subject to either man-made waterflooding or long-term natural waterflooding. By exploring the impact of varying: oil saturation; well patterns; reservoir heterogeneity; and permeability anisotropy, we attempt to quantify the factors that most influence the effectiveness of WAG injection. WAG ratios (the ratio of injected water and CO 2 , in reservoir volumes) of interest are those that either minimize the net CO 2 utilization ratios or maximize oil production rates. In general, the most effective WAG ratios for ROZs, are consistently less than those observed undergoing CO 2 injection in the same geologic reservoir models after traditional (man-made) waterflooding. This work demonstrates that most favorable WAG ratios for oil production in ROZs are different from those in traditional MPZs because of oil saturation differences. Therefore, CO 2 injection into both zones or directly copying WAG injection designs from MPZs to ROZs might not maximize oil production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

The Role of Groundwater Withdrawals on River Regulation: Example From the Columbia River Basin

The Columbia River Basin (CRB) is heavily regulated by more than 250 dams on its river system while depending significantly on groundwater withdrawals in certain sub-basins. Neglecting groundwater withdrawals in hydrologic models of the basin could result in inaccurate predictions of its water budget and thus mislead water management decisions in the basin. This work aims to understand the impacts of groundwater pumping on the spatiotemporal patterns of modeling regulated streamflow in the CRB using a modified version of the Variable Infiltration Capacity (VIC) model integrated with a water management component that accounts for groundwater withdrawals, irrigation demands, and reservoir operation (VIC-GIRR). The VIC-GIRR simulations showed that considering additional groundwater withdrawals would alleviate the stress of irrigation water deficit in the Snake River Basin with an average reduction of 10 km 3 /year. Such a reduction in water deficit resulted in slight streamflow increase over the CRB with maximum increase up to 40% during dry period in certain locations. We also note that the implementation of groundwater withdrawal does not, however, improve the overall model performance in long-term averaged streamflow and storage predictions. Our results highlight the efforts needed to examine additional important processes in representing the interactions between water withdrawals and reservoir operations. Such efforts will aid in better simulation of multi-reservoir system and improve effectiveness for agricultural productivity, power generation, flood control, and navigation purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

AgroIBIS VSF Farmed Pothole Ponding Model Runs

This data set includes modeled (AgroIBIS-VSF) output of crop yield for miscanthus and a corn/soybean rotation in .nc format. Input includes the precipitation scenarios (.csv) used to generate the runs, and supplemental observed yield monitor data (.csv) from the research area of interest.

Agricultural Hydrology (Drainage, Flooding, Irriga↗
DOE OSTI2021

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhancing the representation of water management in global hydrological models

Abstract. This study enhances an existing global hydrological model (GHM), Xanthos, by adding a new water management module that distinguishes between the operational characteristics of irrigation, hydropower, and flood control reservoirs. We remapped reservoirs in the Global Reservoir and Dam (GRanD) database to the 0.5∘ spatial resolution in Xanthos so that a single lumped reservoir exists per grid cell, which yielded 3790 large reservoirs. We implemented unique operation rules for each reservoir type, based on their primary purposes. In particular, hydropower reservoirs have been treated as flood control reservoirs in previous GHM studies, while here, we determined the operation rules for hydropower reservoirs via optimization that maximizes long-term hydropower production. We conducted global simulations using the enhanced Xanthos and validated monthly streamflow for 91 large river basins, where high-quality observed streamflow data were available. A total of 1878 (296 hydropower, 486 irrigation, and 1096 flood control and others) out of the 3790 reservoirs are located in the 91 basins and are part of our reported results. The Kling–Gupta efficiency (KGE) value (after adding the new water management) is ≥ 0.5 and ≥ 0.0 in 39 and 81 basins, respectively. After adding the new water management module, model performance improved for 75 out of 91 basins and worsened for only 7. To measure the relative difference between explicitly representing hydropower reservoirs and representing hydropower reservoirs as flood control reservoirs (as is commonly done in other GHMs), we use the normalized root mean square error (NRMSE) and the coefficient of determination (R2). Out of the 296 hydropower reservoirs, the NRMSE is > 0.25 (i.e., considering 0.25 to represent a moderate difference) for over 44 % of the 296 reservoirs when comparing both the simulated reservoir releases and storage time series between the two simulations. We suggest that correctly representing hydropower reservoirs in GHMs could have important implications for our understanding and management of freshwater resource challenges at regional-to-global scales. This enhanced global water management modeling framework will allow the analysis of future global reservoir development and management from a coupled human–earth system perspective.

13 HYDRO ENERGY↗
DOE OSTI2025

DOE-ICoM/RIFT

Rapid Infrastructure Flood Tool (RIFT) is a two-dimensional hydrodynamic model based on the complete shallow water equations. RIFT has specifically been designed with rapid simulation in mind by utilizing commodity high performance computing technology and best-available nation-wide data. RIFT is used to predict the movement of water over land and resolve the spatial and temporal variability of flood depths, extent, and velocity. RIFT can be applied to many flood situations and has primarily been used to quantify flood extents from dam/levee failure or inland rainfall flooding.

Perkins, Bill [Pacific Northwest National Laborato↗
DOE OSTI2022

Improved Simulations of Atmospheric River Climatology and Variability in High-Resolution CESM

Atmospheric rivers (ARs) can trigger extreme rainfall and snowfall events, causing high risk of flooding. Meanwhile, they are sometimes responsible for the occurrence and break of drought events. So, it is of great importance for weather and climate forecast models to accurately represent ARs. Studies have shown that AR variability is linked to various climate modes, such as El Niño/Southern Oscillation, Pacific Decadal Oscillation, Pacific-North America pattern, providing the potential for predicting ARs at subseasonal-to-decadal timescales. In this study, we examine the impact of climate model resolution on AR simulations and their relationship with major climate modes by comparing a pair of high-resolution and a low-resolution (LR) Community Earth System Model (CESM) simulations against the observations. The results reveal major deficiencies in LR that include underestimated climatological AR strength and associated precipitation, and inaccurate representations of ARs relationships with climate modes. Significant improvements are observed by increasing CESM horizontal resolution, suggesting a path for more accurate prediction of ARs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data for: Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events

This dataset contains data and code associated with the paper Patel et al. "Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events".The coastal terrestrial-aquatic interface (TAI) is a highly dynamic system characterized by strong physical, chemical, and biological gradients. In particular, shifting soil redox conditions and consumption of terminal electron acceptors, due in part to dynamic hydrologic conditions, is a strong driver of carbon availability and transformations across TAIs. However, while redox dynamics are well described, our ability to quantitatively forecast rates of oxic to anoxic shifts in soils with different characteristics and inundation regimes is limited. We integrated field measurements, laboratory incubations, and model simulations to improve mechanistic understanding of oxygen consumption dynamics in coastal soils. Continuous in situ monitoring unexpectedly revealed that flooding caused temporary spikes in subsurface dissolved oxygen followed by rapid consumption in the wetlands. To further investigate these mechanisms in a controlled setting, we performed laboratory incubations using surface and subsurface soils from a TAI gradient (defined here as upland forest to transitional forest to wetland) in Western Lake Erie to measure oxygen consumption rates in TAI soils during flood events. In our experiments, wetland soils reached anoxia the fastest, in ∼ 9 h on average, whereas upland soils turned anoxic in ∼ 18 h. Subsurface upland soils did not turn anoxic even after two weeks of saturation in the lab, and their oxygen consumption patterns suggested carbon and/or nutrient limitation. These results are consistent with in-situ groundwater redox and oxygen measurements in the field, where wetland soils exhibited the highest rates of oxygen consumption along the TAI. Model simulations of oxygen consumption suggested that oxygen consumption had stronger abiotic controls in wetland soils but stronger biotic controls in upland soils, providing a useful framework for future incubation experiments. Microbial activity is a strong driver of oxygen consumption in TAI soils, although it is constrained by the availability of dissolved carbon in subsurface soils.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/COMPASS-DOE/cmps-time_to_anoxia (DOI: 10.5281/zenodo.10815374).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing climate-change-induced flood risk in the Conasauga River watershed: an application of ensemble hydrodynamic inundation modeling

Abstract. This study evaluates the impact of potential future climate change on flood regimes, floodplain protection, and electricity infrastructures across the Conasauga River watershed in the southeastern United States through ensemble hydrodynamic inundation modeling. The ensemble streamflow scenarios were simulated by the Distributed Hydrology Soil Vegetation Model (DHSVM) driven by (1) 1981–2012 Daymet meteorological observations and (2) 11 sets of downscaled global climate models (GCMs) during the 1966–2005 historical and 2011–2050 future periods. Surface inundation was simulated using a GPU-accelerated Two-dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) hydrodynamic model. A total of 9 out of the 11 GCMs exhibit an increase in the mean ensemble flood inundation areas. Moreover, at the 1 % annual exceedance probability level, the flood inundation frequency curves indicate a ∼ 16 km2 increase in floodplain area. The assessment also shows that even after flood-proofing, four of the substations could still be affected in the projected future period. The increase in floodplain area and substation vulnerability highlights the need to account for climate change in floodplain management. Overall, this study provides a proof-of-concept demonstration of how the computationally intensive hydrodynamic inundation modeling can be used to enhance flood frequency maps and vulnerability assessment under the changing climatic conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An Integrated Modeling Framework for Sediment Dynamics During Urban Flooding: Application to Hurricane Harvey in Houston

Floodwater can mobilize and redistribute large volumes of sediment from upland to downstream urban areas, threatening infrastructure, water quality, and ecosystem health. However, existing modeling approaches often fail to capture sediment dynamics in urban floodplains due to the lack of integration between upland hydrological processes and riverine sediment transport. This study presents the first integrated modeling framework that couples the Energy Exascale Earth System Model (E3SM) land component, which simulates runoff and hillslope erosion, with TELEMAC-GAIA, a two-dimensional hydrodynamic and sediment transport model. This framework enables the fully distributed, process-based simulation of high-resolution (as fine as 30 m) sediment dynamics from hillslopes to floodplains. Applied to a highly urbanized watershed in Houston during Hurricane Harvey, this framework reproduced observed water levels at 16 USGS gauges (median R 2 = 0.83 and KGE = 0.78), key sediment dynamics such as sediment transport and deposition processes, and reproduced spatial deposition patterns consistent with LiDAR-derived data. Based on the simulation, we estimate 8.0 million m 3 of event-scale sediment deposition, including 5.7 million m 3 trapped in the flood-control reservoirs and 2.3 million m 3 deposited along major channels and floodplains. Using a representative unit removal cost, this corresponds to an estimated dredging cost of $581 million for total deposition. These results provide a first-order, physically based quantification of Harvey-scale sediment impacts. This study provides a valuable tool for the holistic analysis of sediment dynamics triggered by extreme urban flooding, supporting flood-resilience planning. More broadly, it highlights the importance of integrating physically based hydrological processes for urban flooding and sediment research.

Hurricane Harvey↗
DOE OSTI2026

A hybrid CNN-LSTM surrogate model for hyper-resolution spatiotemporal flood forecasting in Norfolk, Virginia

Study region: Norfolk, Virginia, United States Study focus: Accurate and timely flood forecasting is essential for enhancing resilience in coastal urban areas in the context of increasing frequency and intensity of rainfall, sea level rise and urbanization. This study presents a hybrid deep learning-based surrogate model that integrates Convolutional Neural Networks (CNN) and Long Short-Term Memory (LSTM) networks to enable real-time spatiotemporal flood forecasting. The model leverages CNN to capture spatial features from inputs such as elevation and Topographic Wetness Index (TWI), while LSTM processes time-series inputs of rainfall and tide data to capture temporal features. New hydrologic insights for the region: The hybrid CNN-LSTM model was trained using the physics-based hydrodynamic model simulations obtained from the Two-dimensional Unsteady FLOW (TUFLOW) model for Norfolk, Virginia, and achieved high predictive accuracy across diverse flood-prone areas. The reduced computational time from four to six hours using TUFLOW to 3.2 min per event using CNN-LSTM enables rapid flood inundation mapping and early warning applications. The model effectively captured both spatial flood extents and their temporal evolution across different flooding scenarios, providing forecasts at a 2.5-m spatial resolution and 15-min temporal resolution and a one-hour-ahead prediction horizon. While challenges remain in terms of transferability to new regions and real-time data assimilation, this approach demonstrates strong potential for supporting operational flood risk management in coastal urban environments.

Coastal urban flooding↗
DOE OSTI2024

Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics

Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.

58 GEOSCIENCES↗
DOE OSTI2024

Enhancing Extreme Precipitation Predictions With Dynamical Downscaling: A Convection‐Permitting Modeling Study in Texas and Oklahoma

Precipitation in the Southern Plains of the United States is relatively well depicted by the Community Earth System Model (CESM). However, despite its ability to capture seasonal mean precipitation anomalies, CESM consistently underestimates extreme pluvial and drought events, rendering it an insufficient tool for extending simulation lead times for exceptional events, such as the abnormally dry May 2011, which helped drive Texas into its worst period of drought in more than a century, and the abnormally wet May 2015, which led to widespread flooding in that state. Ensemble‐based regional climate experiments are completed for the two extreme years using Weather Research and Forecasting model (WRF) and downscaled from CESM. WRF simulations are at convection‐permitting grid resolution for improved physical representation of simulated precipitation over the Southern Great Plains. By integrating convection‐permitting models (CPMs) into each individual member of a CESM ocean data assimilation ensemble, this study demonstrates that high‐resolution dynamical downscaling can improve model skillfulness at capturing these two events and is thus a potentially useful tool for forecasting extremely high and extremely low precipitation events at subseasonal or even seasonal lead times.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2026

The Value of Forecasters‐in‐the‐Loop in Real‐Time Flood Forecasting in the Age of Machine Learning

Machine learning (ML) applications in hydrological forecasting are increasingly prevalent and show great potential. However, many previous studies have only evaluated performance through reanalysis or retrospective simulations compared to simplified baselines. This study provides the first assessment of ML performance against actual operational forecasting systems operated by the California Nevada River Forecast Center (CNRFC), which combines the Community Hydrologic Prediction System (CHPS) with forecasters-in-the-loop. Results demonstrate that forecasters-in-the-loop systems consistently outperform ML models in both general forecasts and flood alerting across lead times up to 96 hr, even when ML models use observed forcings, while CNRFC operational process relies on biased weather forecasts. Our analysis reveals that forecaster expertise maintains forecast reliability despite inaccurate precipitation inputs, with human-guided systems showing superior performance degradation characteristics at extended lead times. These findings highlight the irreplaceable value of human expertise in operational forecasting and caution against overstating current ML capabilities in real-world applications.

Tran, Vinh Ngoc [Univ. of Michigan, Ann Arbor, MI ↗
DOE OSTI2022

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Polymer Retention Evaluation in a Heavy Oil Sand for a Polymer Flooding Application on Alaska’s North Slope

For a polymer flooding field trial in a heavy oil reservoir on Alaska’s North Slope, polymer retention is a key parameter. Because of the economic impact of retention, this parameter was extensively studied using field core material and conditions. Several important findings were noted. Polymer retention based on effluent viscosity measurements can be overestimated unless the correct (non-linear) relation between polymer concentration and viscosity is used. Polymer degradation (either mechanical or oxidative) can also lead viscosity-based measurements to overestimate retention. Inaccessible pore volume (IAPV) can be overestimated if insufficient brine is flushed through the sand between polymer banks. Around 100 pore volumes (PV) of brine may be needed to displace mobile polymer to approach a true residual resistance factor and properly measure IAPV. Even for a sand pack with kwsor = 20 mD, IAPV was zero for HPAM with molecular weight (Mw) of 18 million g/mol. Fine-grain particles (<20 μm) strongly impacted polymer retention values. Native NB#1 sand with a significant component of particles <20 μm exhibited 290 μg/g, while the same sand exhibited 28 μg/g after these small particles were removed. Polymer retention did not necessarily correlate with mineral composition. The NB#1, NB#3, and OA sands had similar elemental and clay compositions, but the NB#1 sand exhibited ~10 times higher retention than the NB#3 sand. Polymer retention did not necessarily correlate with permeability. NB#1 sand exhibited much higher retention than OA sand, even though NB#1 sand was twice as permeable as OA sand. No evidence of chromatographic separation of HPAM molecular weights was found in our experiments. Although retention tended to be greater without a residual oil saturation (than at Sor), the effect was not strong. Aging a core (with high oil saturation) at 60°C reduced HPAM retention by a factor of two. As expected, under similar conditions, polymer retention was greater for a higher Mw HPAM (18 million g/mol) than for a lower Mw HPAM (10-12 million g/mol). In many cases with high polymer retention values (e.g., 240 μg/g), polymer arrival at the end of the core was relatively quick, but achieving the injected concentration occurred gradually over many pore volumes. This effect was not due chromatographic separation of polymer molecular weights. Results from modeling of this behavior was consistent with concentrationdependent polymer retention. The form assumed for the retention function in a simulator can have an important impact on the timing and magnitude of the oil response from a polymer flood. Field-based observations can underestimate polymer retention, depending on when the tracer and polymer concentrations were measured and the assumptions made about reservoir heterogeneity.

Wang, Dongmei↗