Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plumes”

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 253 records · Page 14

An adaptable technique for comparative image assessment: Application to crosswell electromagnetic survey design for fluid monitoring

Reservoir integrity stewardship accompanying carbon capture and sequestration considers fluid extraction and reinjection as a risk-mitigating method against overpressuring that could lead to caprock damage and ensuing [Formula: see text] leakage. Crosswell electromagnetics offers a technically viable monitoring method with the spatial volume coverage necessary for reservoir-encompassing pressure management. However, a certain logistic dilemma for deep gas sequestration into saline and thus electrically conductive aquifers is that crosswell magnetic-field measurements underperform in the imaging of more resistive plume bodies, further exacerbated when vertical arrays intersect, as opposed to surround, plumes. Furthermore, comparative synthetic-data plume imaging of such scenarios rates the information content of magnetic-field versus electric-field 3D crosswell layouts for reservoir and infrastructure conditions of a representative pilot site in a coastal area in Florida. The image quality of the resulting plume replications can be ranked numerically through a newly proposed semblance qualifier, appraising the model goodness of fit to a given reference. In contrast to common least-squares measures for goodness of fit, the semblance formulation uses classifying logistic function types, thus enabling a better distinction of predefined anomaly features.

58 GEOSCIENCES↗

Numerical Simulations of Geologic Storage Reservoir Management to Support Risk Mitigation Evaluation

This report provides a detailed description of a set of numerical simulations that represent reservoir behavior over time in response to different operational decision scenarios for detection of potential leakage and reduction or avoidance of leakage impact at a hypothetical geological carbon storage (GCS) site. These simulations serve as the basis for a series of GCS leakage risk forecasts that are to be developed using the National Risk Assessment Partnership’s Open-Source Integrated Assessment Model (NRAP-Open-IAM), and a demonstration of a simple decision support workflow for evaluation of mitigation strategies based on results of those system model forecasts. This risk assessment and decision support study is forthcoming. Four injection scenarios were considered: a constant rate carbon dioxide (CO 2 ) injection case (base case), a case with CO 2 injection rate adjustment, a case with early termination of injection operations, and a case with brine extraction. CO 2 injection operations were controlled to ensure that the pressure transient remains below the defined manageable reservoir fracture pressure, with consideration shown to hypothetical locations within the modeled spatial domain where the overburden was weaker and lower transient pressure increases were allowable. Additionally, a brine extraction alternative was considered as a reservoir management and risk mitigation option to reduce reservoir pressure, steer the plume away from any hypothetical geohazard such as fault as needed, and enhance storage capacity. Such operational actions contribute to risk management overtime. This study explores the potential utility of reservoir management for risk reduction at GCS sites. This study shows that the injection design may modify the time to CO 2 breakthrough at a legacy well; in particular, these results show that brine extraction can add value for mitigating risk both by delaying leakage and reducing pressure build-up. For the scenario considered, both pressure plots and pressure distributions demonstrate that pressure build-up was decreased by 3% with brine extraction. Additionally, extraction of brine afforded enhancement of CO 2 storage capacity by 5% compared to the base case. These findings suggest that brine extraction has substantial potential to steer the risk-related reservoir effects away from known geohazards (e.g., faults and legacy wells) by conducting pressure transient effects and CO 2 plume movement toward the production well. Injection rate adjustment scenarios considered in this study show potential value for managing both reservoir pressure transients and CO 2 plume behavior. Operational actions for reducing injection and/or early termination of injection (as compared to the base case), however, require careful design; tailoring both the extent and timing of injection rate adjustment over the injection and post-injection operational period must be thoroughly planned to balance maximizing storage and minimizing subsurface environmental risk. The study also gives preliminary consideration to the effectiveness that monitoring strategy may play in providing useful information to inform reservoir management decisions for risk reduction. Two types of monitoring were considered: 1) pressure build-up or pressure transient; and 2) potential leakage detection from a CO 2 mass or plume. Four hypothetical legacy wells, two plugged and two abandoned, were placed in the model domain. Monitoring along these wells was measured over time in individual stacked reservoir formations and shale formations to support risk mitigation decisions, especially operational decisions that assisted in risk reduction. These simulations will serve as the basis for a series of GCS leakage risk forecasts that are to be developed using NRAP’s Open-IAM, and demonstration of a simple decision support workflow for comparative assessment of mitigation alternatives based on results of those system model forecasts. This risk assessment and decision support study is forthcoming.

54 ENVIRONMENTAL SCIENCES↗

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗

A Review of Offshore Methane Quantification Methodologies

Since pre-industrial times, anthropogenic methane emissions have increased and are partly responsible for a changing global climate. Natural gas and oil extraction activities are one significant source of anthropogenic methane. While methods have been developed and refined to quantify onshore methane emissions, the ability of methods to directly quantify emissions from offshore production facilities remains largely unknown. Here, we review recent studies that have directly measured emissions from offshore production facilities and critically evaluate the suitability of these measurement strategies for emission quantification in a marine environment. The average methane emissions from production platforms measured using downwind dispersion methods were 32 kg h−1 from 188 platforms; 118 kg h−1 from 104 platforms using mass balance methods; 284 kg h−1 from 151 platforms using aircraft remote sensing; and 19,088 kg h−1 from 10 platforms using satellite remote sensing. Upon review of the methods, we suggest the unusually large emissions, or zero emissions observed could be caused by the effects of a decoupling of the marine boundary layer (MBL). Decoupling can happen when the MBL becomes too deep or when there is cloud cover and results in a stratified MBL with air layers of different depths moving at different speeds. Decoupling could cause: some aircraft remote sensing observations to be biased high (lower wind speed at the height of the plume); the mass balance measurements to be biased high (narrow plume being extrapolated too far vertically) or low (transects miss the plume); and the downwind dispersion measurements much lower than the other methods or zero (plume lofting in a decoupled section of the boundary layer). To date, there has been little research on the marine boundary layer, and guidance on when decoupling happens is not currently available. We suggest an offshore controlled release program could provide a better understanding of these results by explaining how and when stratification happens in the MBL and how this affects quantification methodologies.

Riddick, Stuart N. (ORCID:0000000316841843)↗

Aerosols from anthropogenic and biogenic sources and their interactions – modeling aerosol formation, optical properties, and impacts over the central Amazon basin

The Green Ocean Amazon experiment – GoAmazon 2014–2015 – explored the interactions between natural biogenic forest emissions from central Amazonia and urban air pollution from Manaus. Previous GoAmazon 2014–2015 studies showed that nitrogen oxide (NO x = NO + NO 2 ) and sulfur oxide (SO x ) emissions from Manaus strongly interact with biogenic volatile organic compounds (BVOCs), affecting secondary organic aerosol (SOA) formation. In previous studies, ground-based and aircraft measurements provided evidence of SOA formation and strong changes in aerosol composition and properties. Aerosol optical properties also evolve, and their impacts on the Amazonian ecosystem can be significant. As particles age, some processes, such as SOA production, black carbon (BC) deposition, particle growth and the BC lensing effect change the aerosol optical properties, affecting the solar radiation flux at the surface. This study analyzes data and models SOA formation using the Weather Research and Forecasting with Chemistry (WRF-Chem) model to assess the spatial variability in aerosol optical properties as the Manaus plumes interact with the natural atmosphere. The following aerosol optical properties are investigated: single scattering albedo (SSA), asymmetry parameter (g aer ), absorption Ångström exponent (AAE) and scattering Ångström exponent (SAE). These simulations were validated using ground-based measurements at three experimental sites, namely the Amazon Tall Tower Observatory – ATTO (T0a), downtown Manaus (T1), Tiwa Hotel (T2) and Manacapuru (T3), as well as the U.S. Department of Energy (DOE) Gulfstream 1 (G-1) aircraft flights. WRF-Chem simulations were performed over 7 d during March 2014. Results show a mean biogenic SOA (BSOA) mass enrichment of 512 % at the T1 site, 450 % in regions downwind of Manaus, such as the T3 site, and 850 % in areas north of the T3 site in simulations with anthropogenic emissions. The SOA formation is rather fast, with about 80 % of the SOA mass produced in 3–4 h. Comparing the plume from simulations with and without anthropogenic emissions, SSA shows a downwind reduction of approximately 10 %, 11 % and 6 % at the T1, T2 and T3 sites, respectively. Other regions, such as those further downwind of the T3 site, are also affected. The g aer values increased from 0.62 to 0.74 at the T1 site and from 0.67 to 0.72 at the T3 site when anthropogenic emissions are active. During the Manaus plume-aging process, a plume tracking analysis shows an increase in SSA from 0.91 close to Manaus to 0.98 160 km downwind of Manaus as a result of SOA production and BC deposition.

54 ENVIRONMENTAL SCIENCES↗

On the differences in the vertical distribution of modeled aerosol optical depth over the southeastern Atlantic

Abstract. The southeastern Atlantic is home to an expansive smoke aerosol plume overlying a large cloud deck for approximately a third of the year. The aerosol plume is mainly attributed to the extensive biomass burning activities that occur in southern Africa. Current Earth system models (ESMs) reveal significant differences in their estimates of regional aerosol radiative effects over this region. Such large differences partially stem from uncertainties in the vertical distribution of aerosols in the troposphere. These uncertainties translate into different aerosol optical depths (AODs) in the planetary boundary layer (PBL) and the free troposphere (FT). This study examines differences of AOD fraction in the FT and AOD differences among ESMs (WRF-CAM5, WRF-FINN, GEOS-Chem, EAM-E3SM, ALADIN, GEOS-FP, and MERRA-2) and aircraft-based measurements from the NASA ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) field campaign. Models frequently define the PBL as the well-mixed surface-based layer, but this definition misses the upper parts of decoupled PBLs, in which most low-level clouds occur. To account for the presence of decoupled boundary layers in the models, the height of maximum vertical gradient of specific humidity profiles from each model is used to define PBL heights. Results indicate that the monthly mean contribution of AOD in the FT to the total-column AOD ranges from 44 % to 74 % in September 2016 and from 54 % to 71 % in August 2017 within the region bounded by 25∘ S–0∘ N–S and 15∘ W–15∘ E (excluding land) among the ESMs. ALADIN and GEOS-Chem show similar aerosol plume patterns to a derived above-cloud aerosol product from the Moderate Resolution Imaging Spectroradiometer (MODIS) during September 2016, but none of the models show a similar above-cloud plume pattern to MODIS in August 2017. Using the second-generation High Spectral Resolution Lidar (HSRL-2) to derive an aircraft-based constraint on the AOD and the fractional AOD, we found that WRF-CAM5 produces 40 % less AOD than those from the HSRL-2 measurements, but it performs well at separating AOD fraction between the FT and the PBL. AOD fractions in the FT for GEOS-Chem and EAM-E3SM are, respectively, 10 % and 15 % lower than the AOD fractions from the HSRL-2. Their similar mean AODs reflect a cancellation of high and low AOD biases. Compared with aircraft-based observations, GEOS-FP, MERRA-2, and ALADIN produce 24 %–36 % less AOD and tend to misplace more aerosols in the PBL. The models generally underestimate AODs for measured AODs that are above 0.8, indicating their limitations at reproducing high AODs. The differences in the absolute AOD, FT AOD, and the vertical apportioning of AOD in different models highlight the need to continue improving the accuracy of modeled AOD distributions. These differences affect the sign and magnitude of the net aerosol radiative forcing, especially when aerosols are in contact with clouds.

54 ENVIRONMENTAL SCIENCES↗

Clean Water Production in Cooling Towers

This project developed and demonstrated a novel technology that produces clean water from cooling tower recirculating water by using the natural evaporation and condensation cycle inside cooling towers. The system captures the escaping plume and converts blowdown quality water into high purity water suitable for on-site reuse such as boiler feed. The technology uses electric fields to ionize exhaust plumes, charge the entrained droplets, and direct them toward collection electrodes where they coalesce and flow downward. This allows water recovery at a low energy cost while reducing visible plume emissions. In addition, we developed a complementary software platform that improves overall cooling tower performance. The system uses wireless sensors and physics-based machine learning algorithms to optimize key parameters of the cooling process. For power generation facilities, this increases the thermal efficiency of the cooling loop and condenser, resulting in measurable cycle efficiency gains. Improvements of one percent or more can deliver significant increases in electricity production for the same fuel input.

01 COAL, LIGNITE, AND PEAT↗

Pulsed laser synthesis of mesoporous metal chalcogenide thin films

Mesoporous films of the metal chalcogenide B-FeSe were grown on MgO substrates by KrF pulsed laser deposition (PLD) in an argon background. At 100 mTorr, gated intensified charge-coupled device imaging and ion probe measurements showed that the plasma plume responsible for crystal growth initially comprised three components, with distinct expansion velocities. Plume interactions with the substrate heater and ablation target gave rise to complex dynamics, including collisions between the charged leading edge—rebounding between the substrate and the target—and slower-moving species in the plume interior. Film growth was dominated by species with kinetic energies ≤0.5 eV/atom. X-ray reflectivity revealed that films grown in this environment—with a substrate temperature of 350 ° C, a laser fluence of 1.0 J cm −2 , and a 7.5 mm 2 spot area—formed a porous framework with 15% porosity. Atomic force microscopy showed surface features that suggest pore sizes below 100 nm. X-ray diffraction indicated that the porous films were epitaxial with respect to the substrate and likely grew by oriented-attachment of gas-phase molecular clusters or very small nanoparticles, in contrast to the conventional epitaxy of vacuum films from atomic constituents. The in-plane orientation of the mesoporous films was B-FeSe [100]||[110] MgO, attributed to the soft landing of pre-formed crystallites on the MgO substrates, where protruding Se rows of B-FeSe aligned with corrugations of the MgO surface. In conclusion, this work implies that growth of candidate electrocatalyst materials by PLD in inert gas background may allow mesoporous frameworks with a single crystallographic orientation that expose specific crystal facets for electrochemical reactions and active site engineering.

cluster epitaxy↗

Well–based monitoring of CO 2 geological sequestration operations in saline aquifers: Critical insights into key questions

Geological carbon sequestration in saline aquifers is one of the most promising strategies to help mitigate emissions of CO 2 to the atmosphere. Significant challenges in ensuring the security of the sequestration process rest in the evolution and expansion of the CO 2 plume in the subsurface. The ability to track the movement of the injected CO 2 poses another challenge. Critical questions related to the integrity of the sequestration operations in saline aquifers relate to plume characteristics that can we monitor using well-based variables. We addressed this and related questions using an integrated modeling framework through a numerical investigation of carbon sequestration in saline aquifers during long-term and post-injection periods. This modeling paradigm incorporates the effect of structural, geological, and petrophysical characteristics. That way, we can account for critical physicochemical processes, rock-fluid interactions, and lithology dependencies. The well fluid variables investigated include fluid composition, pH, fluid density, and ion activity. We learned that fluid property analytics and diagnostics can be powerful tools to estimate the movement of CO 2 and its storage in different trapping mechanisms. These analytics can help optimize operational aspects and simplify reservoir-scale models while still reflecting the complex nature of the CO 2 interactions underground and offering insights into plume evolution.

58 GEOSCIENCES↗

Challenges of designing, operating, and managing brine extraction at a carbon dioxide storage site

Developing and deploying commercial-scale carbon dioxide (CO 2 ) storage at a site requires operational strategies to manage pressure and CO 2 plume distribution. In the natural gas storage and oil and gas industries, fluid extraction has been successfully implemented as a strategy to manage pressure and control fluid distribution (i.e., hydrocarbon recovery and natural gas storage efficiency, respectively). The technical challenges of integrating brine extraction at CO 2 storage sites are expected to be similar, if not identical, to those in the natural gas storage and petroleum industries. Major technical challenges of using brine extraction to manage storage in a storage unit include (1) designing and completing the brine extraction wells, (2) locating the brine extraction well, (3) assessing the efficacy of brine extraction, and (4) handling the extracted brine. The design and completion of the brine extraction well include the well type (i.e., vertical or horizontal) and perforated interval. The location of the well is influenced by the timing of extraction operations (i.e., prior or after CO 2 injection starts) and project objectives (i.e., desired directional movement of the CO 2 plume and avoid or minimize CO 2 extraction). The efficacy of brine extraction requires a baseline projection of the pressure and CO 2 plume distribution. Options for handling extracted brine include subsurface disposal, surface disposal, and industrial use. Here this paper addresses means of overcoming technical challenges that may likely be encountered before or after CO 2 injection starts.

25 ENERGY STORAGE↗

Evaluating the state-of-the-art in remote volcanic eruption characterization Part I: Raikoke volcano, Kuril Islands

Raikoke, a small, unmonitored volcano in the Kuril Islands, erupted in June 2019. We integrate data from satellites (including Sentinel-2, TROPOMI, MODIS, Himawari-8), the International Monitoring System (IMS) infrasound network, and global lightning detection network (GLD360) with information from local authorities and social media to retrospectively characterize the eruptive sequence and improve understanding of the pre-, syn- and post- eruptive behavior. In this work, we observe six infrasound pulses beginning on 21 June at 17:49:55 UTC as well as the main Plinian phase on 21 June at 22:29 UTC. Each pulse is tracked in space and time using lightning and satellite imagery as the plumes drift eastward. Post-eruption visible satellite imagery shows expansion of the island's surface area, an increase in crater size, and a possibly-linked algal bloom south of the island. We use thermal satellite imagery and plume modeling to estimate plume height at 10–12 km asl and 1.5–2 × 10 6 kg/s mass eruption rate. Remote infrasound data provide insight into syn-eruptive changes in eruption intensity. Our analysis illustrates the value of interdisciplinary analyses of remote data to illuminate eruptive processes. However, our inability to identify deformation, pre-eruptive outgassing, and thermal signals, which may reflect the relatively short duration (~12 h) of the eruption and minimal land area around the volcano and/or the character of closed-system eruptions, highlights current limitations in the application of remote sensing for eruption detection and characterization.

58 GEOSCIENCES↗

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

The composition of gases from a diffusion flame above longleaf pine needle fuel beds

The gas and tar composition of a diffusion flame from longleaf pine needles is currently poorly understood and more data are needed to fill in the gap between pyrolysis data and smoke plume data, thus improving physical and chemical modeling of wildland smoke formation. A pilot experiment to measure light gas and tar composition of such a flame is described for three flame regions: persistent flame (flame base), intermittent flame, and smoke plume. Flame gases from 24 experimental fires were collected in canisters and analyzed using EPA method TO-14A for CO 2 , CO, H 2 , CH 4 , and C 2 to C 7 hydrocarbon gases. Condensed gas (tar) samples were collected and analyzed using GC/MS. Other light gases were measured using FTIR spectroscopy. Results from compositional data analysis suggest significant differences in (relative) concentration of compounds detected in the three regions of the flame. Statistical tests for differences in flame zones were performed using the canister data: Concentration of hydrocarbons relative to CO and CO 2 decreased from the persistent flame zone above the pyrolyzing needles through the intermittent flame region into the flame-free plume. This was likely due to both chemical reactions (oxidation) occurring in the flame as well as the introduction of air into the flame/plume by entrainment.

Biomass↗

Impacts on mechanical strength of chemical reactions induced by hydrous supercritical CO2 in Boise Sandstone

Geomechanics experiments were used to assess mechanical alteration of Boise Sandstone promoted by reactions with supercritical carbon dioxide (scCO 2 ) and water vapor. During geologic carbon storage, scCO 2 is injected into subsurface reservoirs, forming buoyant plumes. At brine-plume interfaces, scCO 2 can dissolve into native brines, and water from brines can partition into scCO 2 , forming hydrous scCO 2 . This study investigates the effect of hydrous scCO 2 on the strength of Boise Sandstone. Samples are first exposed to recirculating hydrous scCO 2 for 24 h at 70 °C and 13.8 MPa scCO 2 pressure. Samples are reacted with scCO 2 with added water contents up to 500 mL. After scCO 2 exposure, samples are deformed at room temperature under confining pressures of 3.4, 6.9, and 10.3 MPa. The results demonstrate that hydrous scCO 2 induces chemical reactions in Boise Sandstone, with ions migrating from the solid into the hydrous scCO 2 phase. At the longer time-scales, these reactions could lead to mechanical weakening in the samples; however, on the scale of our experiments, the strength changes are within sample variability. Because the solubility of water in scCO 2 is extremely low (0.008 mol H 2 O per 1 mol CO 2 ), the mineral dissolution of Boise Sandstone was under 0.002 wt.%. Additionally, mineral grains and pore throats in Boise Sandstone are cemented with quartz, which is not susceptible to dissolution at these conditions. Our results indicate that humidity in scCO 2 plumes is unlikely to sustain chemical reactions and induce long term strength changes in quartz cemented sandstones due to resistant mineralogies and low water solubility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗