Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geologic Map”

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.

571 records · Page 32

Changing seasonality of moderate and extreme precipitation events in the Alps

The intensity of precipitation events is expected to increase in the future.The rate of increase depends on the strength or rarity of the events; very strong and rare events tend to follow the Clausius–Clapeyron relation,whereas weaker events or precipitation averages increase at a smaller rate than expected from the Clausius–Clapeyron relation. An often overlooked aspect is seasonal occurrence of such events, which might change in the future. To address the impact of seasonality, we use a large ensemble of regional and global climate model simulations, comprising tens of thousands of model years of daily temperature and precipitation for the past, present,and future. In order to make the data comparable, they are quantile mapped to observation-based time series representative of the Aare catchment in Switzerland. Model simulations show no increase in annual maximum 1-day precipitation events (Rx1day) over the last 400 years and an increase of 10 %–20 % until the end of the century for a strong (RCP8.5) forcing scenario. This fits with a Clausius–Clapeyron scaling of temperature at the event day, which increases less than annual mean temperature. An important reason for this is a shift in seasonality. Rx1day events become less frequent in late summer and more frequent in early summer and early autumn, when it is cooler. The seasonality shift is shown to be related to summer drying. Models with decreasing annual mean or summer mean precipitation show this behaviour more strongly. The highest Rx1day per decade, in contrast, shows no change in seasonality in the future. This discrepancy implies that decadal-scale extremes are thermodynamically limited; conditions conducive to strong events still occur during the hottest time of the year on a decadal scale. In contrast, Rx1day events are also limited by other factors. Conducive conditions are not reached every summer in the present, and even less so in the future. Finally, results suggest that changes in the seasonal cycle need to be accounted for when preparing for moderately extreme precipitation events and assessing their socio-economic impacts.

54 ENVIRONMENTAL SCIENCES↗

STAR FORMATION ACTIVITY IN THE GALACTIC H II COMPLEX S255-S257

We present results on the star formation activity of an optically obscured region containing an embedded cluster (S255-IR) and molecular gas between two evolved H II regions, S255 and S257. We have studied the complex using optical and near-infrared (NIR) imaging, optical spectroscopy, and radio continuum mapping at 15 GHz, along with Spitzer-IRAC results. We found that the main exciting sources of the evolved H II regions S255 and S257 and the compact H II regions associated with S255-IR are of O9.5-B3 V nature, consistent with previous observations. Our NIR observations reveal 109 likely young stellar object (YSO) candidates in an area of {approx}4.'9 x 4.'9 centered on S255-IR, which include 69 new YSO candidates. To see the global star formation, we constructed the V - I/V diagram for 51 optically identified IRAC YSOs in an area of {approx}13' x 13' centered on S255-IR. We suggest that these YSOs have an approximate age between 0.1 and 4 Myr, indicating a non-coeval star formation. Using spectral energy distribution models, we constrained physical properties and evolutionary status of 31 and 16 YSO candidates outside and inside the gas ridge, respectively. The models suggest that the sources associated with the gas ridge are younger (mean age {approx}1.2 Myr) than the sources outside the gas ridge (mean age {approx}2.5 Myr). The positions of the young sources inside the gas ridge at the interface of the H II regions S255 and S257 favor a site of induced star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Metals in Past Societies: A Global Perspective on Indigenous African Metallurgy Shadreck Chirikure

This slim book (166 pages) shines a spotlight on pre-industrial African metallurgy, its global connections, and anthropological implications. It integrates seemingly disparate disciplines, such as history, geology, ethnography, archeology, and metallurgy, to illustrate the diversity and innovation in metallurgy across Africa and the role of metals in the rise of socio-economic inequalities and political power. The book has 7 chapters and the focus on metals as enablers of human needs and wants is evident in each chapter. The first chapter presents the context of the work and data sources. The second chapter focuses on the origin and development of mining and metallurgy in pre-industrial Africa. Chapter 3 is dedicated to the interaction of nature and culture in the process of mining. Chapter 4 deals with the transformation of the ore into metal by smelting and the sociocultural aspects of this process. Chapter 5 explores the social and cultural roles acquired by metals as a result of fabrication into objects. Chapter 6 examines the social role of metals, trade in metals, cultural contact, proto-globalization, and technology transfer. Finally, Chapter 7 draws lessons for global anthropology from the African experience. The sources of information are adequately cited and the long list of references at the end of each chapter will be a boon to researchers in this field. The author highlights the cultural aspects and social context of the adoption of metallurgy in Africa while drawing parallels between practices in pre-industrial Africa and those in other parts of the world. The book is peppered with delightful vignettes that offer insights into the process of transforming nature into culturally significant objects. For instance, African miners, like their counterparts in Nepal and Latin America, called upon deities, spirits and ancestors to mediate between nature and humans. Women had distinct roles in this process, but there were variations in these roles and in the caste status of metal workers across Africa. Taboos, rituals and magic were very much a part of the development of metallurgical technology. The smelting of metal was considered analogous to conception, gestation, and child birth. Power, fertility and metallurgy were intertwined as shown by the decoration of furnaces with female anatomical features. The book concludes with a warning against broad generalizations and stereotypes about African practices. This volume is well written and illustrated with photos, micrographs, colorful maps and illustrations that make it eminently readable.

36 MATERIALS SCIENCE↗

Recover Act. Verification of Geothermal Tracer Methods in Highly Constrained Field Experiments

The prediction of the geothermal system efficiency is strong linked to the character of the flow system that connects injector and producer wells. If water flow develops channels or “short circuiting” between injection and extraction wells thermal sweep is poor and much of the reservoir is left untapped. The purpose of this project was to understand how channelized flow develops in fracture geothermal reservoirs and how it can be measured in the field. We explored two methods of assessing channelization: hydraulic connectivity tests and tracer tests. These methods were tested at a field site using two verification methods: ground penetrating radar (GPR) images of saline tracer and heat transfer measurements using distributed temperature sensing (DTS). The field site for these studies was the Altona Flat Fractured Rock Research Site located in northeastern New York State. Altona Flat Rock is an experimental site considered a geologic analog for some geothermal reservoirs given its low matrix porosity. Because soil overburden is thin, it provided unique access to saturated bedrock fractures and the ability image using GPR which does not effectively penetrate most soils. Five boreholes were drilled in a “five spot” pattern covering 100 m2 and hydraulically isolated in a single bedding plane fracture. This simple system allowed a complete characterization of the fracture. Nine small diameter boreholes were drilled from the surface to just above the fracture to allow the measurement of heat transfer between the fracture and the rock matrix. The focus of the hydraulic investigation was periodic hydraulic testing. In such tests, rather than pumping or injection in a well at a constant rate, flow is varied to produce an oscillating pressure signal. This pressure signal is sensed in other wells and the attenuation and phase lag between the source and receptor is an indication of hydraulic connection. We found that these tests were much more effective than constant pumping tests in identifying a poorly connected well. As a result, we were able to predict which well pairs would demonstrate channelized flow. The focus of the tracer investigation was multi-ionic tests. In multi-ionic tests several ionic tracers are injected simultaneously and the detected in a nearby pumping well. The time history of concentration, or breakthrough curve, will show a separation of the tracers. Anionic tracers travel with the water but cationic tracer undergo chemical exchange with cations on the surface of the rock. The degree of separation is indicative of the surface area exposed to the tracer. Consequently, flow channelization will tend to decrease the separation in the breakthrough. Estimation of specific surface area (the ration of fracture surface area to formation volume) is performed through matching the breakthrough curve with a transport model. We found that the tracer estimates of surface area were confirmed the prediction of channelized flow between well pairs produced by the periodic hydraulic tests. To confirm that the hydraulic and tracer tests were correctly predicting channelize flow, we imaged the flow field using surface GPR. Saline water was injected between the well pairs which produced a change in the amplitude and phase of the reflected radar signal. A map was produced of the migration of saline tracer from these tests which qualitatively confirmed the flow channelization predicted by the hydraulic and tracer tests. The resolution of the GPR was insufficient to quantitatively estimate swept surface area, however. Surface GPR is not applicable in typical geothermal fields because the penetration depths do not exceed 10’s of meters. Nevertheless, the method of using of phase to measure electrical conductivity and the assessment of antennae polarization represent a significant advancement in the field of surface GPR. The effect of flow character on fracture / rock thermal exchange was evaluated using heated water as a tracer. Water elevated 30 degrees C above the formation water was circulated between two wells pairs. One well pair had been identified in hydraulic and tracer testing as well connected and the other poorly connected. Temperature rise was measured in the adjacent rock matrix using coiled fiber optic cable interrogated for temperature using a DTS. This experimental design produced over 4000 temperature measurements every hour. We found that heat transfer between the fracture and the rock matrix was highly impacted by the character of the flow field. The strongly connected wells which had demonstrated flow channelization produced heat rise in a much more limited area than the more poorly connected wells. In addition, the heat increase followed the natural permeability of the fracture rather than the induced flow field. The primary findings of this work are (1) even in a single relatively planar fracture, the flow field can be highly heterogeneous and exhibit flow channeling, (2) channeling results from a combination of fracture permeability structure and the induced flow field, and (3) flow channeling leads to reduced heat transfer. Multi-ionic tracers effectively estimate relative surface area but an estimate of ion-exchange coefficients are necessary to provide an absolute measure of specific surface area. Periodic hydraulic tests also proved a relative indicator of connectivity but cannot prove an absolute measure of specific surface area.

15 GEOTHERMAL ENERGY↗

Low-temperature Stirling Engine for Geothermal Electricity Generation

Up to 2700 terawatt-hours per year of geothermal electricity generation capacity has been shown to be available within North America, typically with wells drilled into geologically active regions of the earth's crust where this energy is concentrated (Huttrer, 2001). Of this potential, about half is considered to have temperatures high enough for conventional (steam-based) power production, while the other half requires unconventional power conversion approaches, such as organic Rankine cycle systems or Stirling engines. If captured and converted effectively, geothermal power generation could replace up to 100GW of fossil fuel electric power generation, leading to a significant reduction of US power sector emissions. In addition, with the rapid growth of hydro-fracking in oil and gas production, there are smaller-scale distributed power generation opportunities in heated liquids that are co-produced with the main products. Since 2006, Cool Energy, Inc. (CEI) has designed, fabricated and tested four generations of low-temperature (100°C to 300°C) Stirling engine power conversion equipment. The electric power output of these engines has been demonstrated at over 2kWe and over 16% thermal conversion efficiency for an input temperature of 215°C and a rejection temperature of 15°C. Initial pilot units have been shipped to development partners for further testing and validation, and significantly larger engines (20+ kWe) have been shown to be feasible and conceptually designed. Originally intended for waste heat recovery (WHR) applications, these engines are easily adaptable to geothermal heat sources, as the heat supply temperatures are similar. Both the current and the 20+ kWe designs use novel approaches of self-lubricating, low-wear-rate bearing surfaces, non-metallic regenerators, and high-effectiveness heat exchangers. By extending CEI's current 3 kWe SolarHeart® Engine into the tens of kWe range, many additional applications are possible, as one 20 kWe design produces nearly seven times the power output of the 3 kWe unit but at only 2.5 times the estimated fabrication cost. Phase I of the proposed SBIR program will therefore study the feasibility of generating electricity with one or more 20 kWe or larger Stirling engines, powered by geothermal heat produced by current and possibly some forward-looking borehole extraction methods, and from producing oil and gas wells. The feasibility study will include full analysis of the thermodynamic and heat transfer processes within the engine (necessary to produce optimum theoretical designs and performance maps), the cost of pumping the geothermal heat recovery fluid, and how the system tradeoffs impact the overall system economics. The goal is a geothermal system design that could be demonstrated during a Phase II follow-on program at a field test site.

15 GEOTHERMAL ENERGY↗

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗

Risk Assessment and Monitoring of Stored CO2 in Organic Rocks Under Non-Equilibrium Conditions

The USA is embarking upon tackling the serious environmental challenges posed to the world by greenhouse gases, especially carbon dioxide (CO2). The dimension of the problem is daunting. In fact, according to the Energy Information Agency, nearly 6 billion metric tons of CO2 were produced in the USA in 2007 with coal-burning power plants contributing about 2 billion metric tons. To mitigate the concerns associated with CO2 emission, geological sequestration holds promise. Among the potential geological storage sites, unmineable coal seams and shale formations in particular show promise because of the probability of methane recovery while sequestering the CO2. However. the success of large-scale sequestration of CO2 in coal and shale would hinge on a thorough understanding of CO2's interactions with host reservoirs. An important parameter for successful storage of CO2 reservoirs would be whether the pressurized CO2 would remain invariant in coal and shale formations under reasonable internal and/or external perturbations. Recent research has brought to the fore the potential of induced seismicity, which may result in caprock compromise. Therefore, to evaluate the potential risks involved in sequestering CO2 in Illinois bituminous coal seams and shale, we studied: (i) the mechanical behavior of Murphysboro (Illinois) and Houchin Creek (Illinois) coals, (ii) thermodynamic behavior of Illinois bituminous coal at - 100oC ≤ T ≤ 300oC, (iii) how high pressure CO2 (up to 20.7 MPa) modifies the viscosity of the host, (iv) the rate of emission of CO2 from Illinois bituminous coal and shale cores if the cores, which were pressurized with high pressure (≤ 20.7 MPa) CO2, were exposed to an atmospheric pressure, simulating the development of leakage pathways, (v) whether there are any fractions of CO2 stored in these hosts which are resistance to emission by simply exposing the cores to atmospheric pressure, and (vi) how compressive shockwaves applied to the coal and shale cores, which were pressurized with high pressure CO2, determine the fate of sequestered CO2 in these cores. Our results suggested that Illinois bituminous coal in its unperturbed state, i.e., when not pressurized with CO2, showed large variations in the mechanical properties. Modulus varied from 0.7 GPa to 3.4 GPa even though samples were extracted from a single large chunk of coal. We did not observe any glass transition for Illinois bituminous coal at - 100oC ≤ T ≤ 300oC, however, when the coal was pressurized with CO2 at ambient ≤ P ≤ 20.7 MPa, the viscosity of the coal decreased and inversely scaled with the CO2 pressure. The decrease in viscosity as a function of pressure could pose CO2 injection problems for coal as lower viscosity would allow the solid coal to flow to plug the fractures, fissures, and cleats. Our experiments also showed a very small fraction of CO2 was absorbed in coal; and when CO2 pressurized coals were exposed to atmospheric conditions, the loss of CO2 from coals was massive. Half of the sequestered gas from the coal cores was lost in less than 20 minutes. Our shockwave experiments on Illinois bituminous coal, New Albany shale (Illinois), Devonian shale (Ohio), and Utica shale (Ohio) presented clear evidence that the significant emission of the sequestered CO2 from these formations cannot be discounted during seismic activity, especially if caprock is compromised. It is argued that additional shockwave studies, both compressive and transverse, would be required for successfully mapping the risks associated with sequestering high pressure CO2 in coal and shale formations.

Malhotra, Vivak↗

Global rules for translating land-use change (LUH2) to land-cover change for CMIP6 using GLM2

Anthropogenic land-use and land-cover change activities play a critical role in Earth system dynamics through significant alterations to biogeophysical and biogeochemical properties at local to global scales. To quantify the magnitude of these impacts, climate models need consistent land-cover change time series at a global scale, based on land-use information from observations or dedicated land-use change models. However, a specific land-use change cannot be unambiguously mapped to a specific land-cover change. Here, nine translation rules are evaluated based on assumptions about the way land-use change could potentially impact land cover. Utilizing the Global Land-use Model 2 (GLM2), the model underlying the latest Land-Use Harmonization dataset (LUH2), the land-cover dynamics resulting from land-use change were simulated based on multiple alternative translation rules from 850 to 2015 globally. For each rule, the resulting forest cover, carbon density and carbon emissions were compared with independent estimates from remote sensing observations, U.N. Food and Agricultural Organization reports, and other studies. The translation rule previously suggested by the authors of the HYDE 3.2 dataset, that underlies LUH2, is consistent with the results of our examinations at global, country and grid scales. This rule recommends that for CMIP6 simulations, models should (1) completely clear vegetation in land-use changes from primary and secondary land (including both forested and non-forested) to cropland, urban land and managed pasture; (2) completely clear vegetation in land-use changes from primary forest and/or secondary forest to rangeland; (3) keep vegetation inland-use changes from primary non-forest and/or secondary non-forest to rangeland. Our analysis shows that this rule is one of three (out of nine)rules that produce comparable estimates of forest cover, vegetation carbon and emissions to independent estimates and also mitigate the anomalously high carbon emissions from land-use change observed in previous studies in the 1950s. According to the three translation rules, contemporary global forest area is estimated to be 37.42×106 km2, within the range derived from remote sensing products. Likewise, the estimated carbon stock is in close agreement with reference biomass datasets, particularly over regions with more than 50 % forest cover

58 GEOSCIENCES↗

Geometry and spatial distribution of lenticulae on Europa

The surface of Europa contains several types of roughly elliptical features, collectively called lenticulae. Lenticulae may have positive relief (domes) or negative relief (pits), may disrupt the crust (chaos), or discolor the surface (spots); some lenticulae have attributes of both domes and chaos (dome/chaos). In this study, we map the location, dimensions and shapes of all lenticulae and their interactions with other lenticulae and lineaments. We find that (1) pits and domes have similar sizes; (2) chaos are larger than the other lenticulae; (3) pits are clustered within the trailing antijovian quadrant and the leading subjovian quadrant whereas domes, dome/chaos, and chaos terrains are more uniformly distributed; (4) the areal density for all lenticulae is not uniform; (5) lenticulae do not divert the path of younger lineaments such as ridges. Taken together, these observations are consistent with conceptual models in which lenticulae are created by intrusion of liquid water bodies, or convection within, the ice shell. Additionally, the observations are consistent with the notion that each type of lenticula is a surface expression of dynamics within the ice shell at a different stage of lenticulae evolution. The similar size and shape of pits and domes suggests that one may evolve into the other. Because domes are more numerous and more uniformly distributed than pits, they are more likely to represent the end stage of this evolution, assuming the end-stage leaves a longer-lasting surface expression. Models also predict that larger features are more likely to disrupt the crust, which is consistent with dome/chaos and chaos being larger than pits and domes. We find no examples of lineaments offsetting pits but lineaments do cross some chaos. Pits also have a preferred northwest-southeast elongation, whereas domes, dome/chaos, and chaos do not have a preferred orientation. If lenticulae orientation is influenced by crustal stress, then pits may have formed during a shorter time interval than the other features. As a result, pits may sample a shorter, more recent time period than domes, dome/chaos, and chaos, consistent with pits being the earliest stage in the evolution of lenticulae. In conclusion, we find no strong evidence that lineaments are deflected by lenticulae, implying either that the stresses created by lenticulae are too small to influence lineaments, or that the complete evolution of lenticulae occurs on a time scale that is short compared to the time between the formation of lineaments at a given location.

58 GEOSCIENCES↗

Life Cycle Water Consumption and Water Resource Assessment for Utility-Scale Geothermal Systems: An In-Depth Analysis of Historical and Forthcoming EGS Projects

This report is the third in a series of reports sponsored by the U.S. Department of Energy Geothermal Technologies Program in which a range of water-related issues surrounding geothermal power production are evaluated. The first report made an initial attempt at quantifying the life cycle fresh water requirements of geothermal power-generating systems and explored operational and environmental concerns related to the geochemical composition of geothermal fluids. The initial analysis of life cycle fresh water consumption of geothermal power-generating systems identified that operational water requirements consumed the vast majority of water across the life cycle. However, it relied upon limited operational water consumption data and did not account for belowground operational losses for enhanced geothermal systems (EGSs). A second report presented an initial assessment of fresh water demand for future growth in utility-scale geothermal power generation. The current analysis builds upon this work to improve life cycle fresh water consumption estimates and incorporates regional water availability into the resource assessment to improve the identification of areas where future growth in geothermal electricity generation may encounter water challenges. This report is divided into nine chapters. Chapter 1 gives the background of the project and its purpose, which is to assess the water consumption of geothermal technologies and identify areas where water availability may present a challenge to utility-scale geothermal development. Water consumption refers to the water that is withdrawn from a resource such as a river, lake, or nongeothermal aquifer that is not returned to that resource. The geothermal electricity generation technologies evaluated in this study include conventional hydrothermal flash and binary systems, as well as EGSs that rely on engineering a productive reservoir where heat exists, but where water availability or permeability may be limited. Chapter 2 describes the approach and methods for this work and identifies the four power plant scenarios evaluated: a 20-MW EGS binary plant, a 50-MW EGS binary plant, a 10-MW hydrothermal binary plant, and a 50-MW hydrothermal flash plant. The methods focus on (1) the collection of data to improve estimation of EGS stimulation volumes, aboveground operational consumption for all geothermal technologies, and belowground operational consumption for EGS; and (2) the mapping of the geothermal and water resources of the western United States to assist in the identification of potential water challenges to geothermal growth. Chapters 3 and 4 present the water requirements for the power plant life cycle. Chapter 3 presents the results of the current data collection effort, and Chapter 4 presents the normalized volume of fresh water consumed at each life cycle stage per lifetime energy output for the power plant scenarios evaluated. Over the life cycle of a geothermal power plant, from construction through 30 years of operation, the majority of water is consumed by plant operations. For the EGS binary scenarios, where dry cooling was assumed, belowground operational water loss is the greatest contributor depending upon the physical and operational conditions of the reservoir. Total life cycle water consumption requirements for air-cooled EGS binary scenarios vary between 0.22 and 1.85 gal/kWh, depending upon the extent of belowground operational water consumption. The air-cooled hydrothermal binary and flash plants experience far less fresh water consumption over the life cycle, at 0.04 gal/kWh. Fresh water requirements associated with air- cooled binary operations are primarily from aboveground water needs, including dust control, maintenance, and domestic use. Although wet-cooled hydrothermal flash systems require water for cooling, these plants generally rely upon the geofluid, fluid from the geothermal reservoir, which typically has high salinity and total dissolved solids concentration and is much warmer than normal groundwater sources, for their cooling water needs; thus, while there is considerable geofluid loss at 2.7 gal/kWh, fresh water consumption during operations is similar to that of aircooled binary systems. Chapter 5 presents the assessment of water demand for future growth in deployment of utility-scale geothermal power generation. The approach combines the life cycle analysis of geothermal water consumption with a geothermal supply curve according to resource type, levelized cost of electricity (LCOE), and potential growth scenarios. A total of 17 growth scenarios were evaluated. In general, the scenarios that assumed lower costs for EGSs as a result of learning and technological improvements resulted in greater geothermal potential, but also significantly greater water demand due to the higher water consumption by EGSs. It was shown, however, that this effect could be largely mitigated if nonpotable water sources were used for belowground operational water demands. The geographical areas that showed the highest water demand for most growth scenarios were southern and northern California, as well as most of Nevada. In addition to water demand by geothermal power production, Chapter 5 includes data on water availability for geothermal development areas. A qualitative analysis is included that identifies some of the basins where the limited availability of water is most likely to affect the development of geothermal resources. The data indicate that water availability is fairly limited, especially under drought conditions, in most of the areas with significant near- and medium-term geothermal potential. Southern California was found to have the greatest potential for water-related challenges with its combination of high geothermal potential and limited water availability. The results of this work are summarized in Chapter 6. Overall, this work highlights the importance of utilizing dry cooling systems for binary and EGS systems and minimizing fresh water consumption throughout the life cycle of geothermal power development. The large resource base for EGSs represents a major opportunity for the geothermal industry; however, depending upon geology, these systems can require large quantities of makeup water due to belowground reservoir losses. Identifying potential sources of compatible degraded or low-quality water for use for makeup injection for EGS and flash systems represents an important opportunity to reduce the impacts of geothermal development on fresh water resources. The importance of identifying alternative water sources for geothermal systems is heightened by the fact that a large fraction of the geothermal resource is located in areas already experiencing water stress. Chapter 7 is a glossary of the technical terms used in the report, and Chapters 8 and 9 provide references and a bibliography, respectively.

15 GEOTHERMAL ENERGY↗

Assessment of Brine Management for Geologic Carbon Sequestration

Geologic carbon sequestration (GCS) is the injection of carbon dioxide (CO 2 ), typically captured from stationary emission sources, into deep geologic formations to prevent its entry into the atmosphere. Active pilot facilities run by regional United States (US) carbon sequestration partnerships inject on the order of one million metric tonnes (mt) CO 2 annually while the US electric power sector emits over 2000 million mt-CO 2 annually. GCS is likely to play an increasing role in US carbon mitigation initiatives, but scaling up GCS poses several challenges. Injecting CO 2 into sedimentary basins raises fluid pressure in the pore space, which is typically already occupied by naturally occurring, or native, brine. The resulting elevated pore pressures increase the likelihood of induced seismicity, of brine or CO 2 escaping into potable groundwater resources, and of CO 2 escaping into the atmosphere. Brine extraction is one method for pressure management, in which brine in the injection formation is brought to the surface through extraction wells. Removal of the brine makes room for the CO 2 and decreases pressurization. Although the technology required for brine extraction is mature, this form of pressure management will only be applicable if there are cost-­effective and sustainable methods of disposing of the extracted brine. Brine extraction, treatment, and disposal may increase the already substantial capital, energy, and water demands of Carbon dioxide Capture and Sequestration (CCS). But, regionally specific brine management strategies may be able to treat the extracted water as a source of revenue, energy, and water to subsidize CCS costs, while minimizing environmental impacts. By this approach, value from the extracted water would be recovered before disposing of any resulting byproducts. Until a price is placed on carbon, we expect that utilities and other CO 2 sources will be reluctant to invest in capital intensive, high risk GCS projects; early technical, economic, and environmental assessments of brine management are extremely valuable for determining the potential role of GCS in the US. We performed a first order feasibility and economic assessment, at three different locations in the US, of twelve GCS extracted-­water management options, including: geothermal energy extraction, desalination, salt and mineral harvesting, rare-­earth element harvesting, aquaculture, algae biodiesel production, road de-­icing, enhanced geothermal system (EGS) recharge, underground reinjection, landfill disposal, ocean disposal, and evaporation pond disposal. Three saline aquifers from different regions of the US were selected as hypothetical GCS project sites to encompass variation in parameters that are relevant to the feasibility and economics of brine disposal. The three aquifers are the southern Mt. Simon Sandstone Formation in the Illinois Basin, IL; the Vedder Formation in the southern San Joaquin Basin, CA; and the Jasper Interval in the eastern Texas Gulf Basin, TX. These aquifers are candidates for GCS due to their physical characteristics and their close proximity to large CO 2 emission sources. Feasibility and impacts were calculated using one mt-­CO 2 injected as the functional unit of brine management. Scenarios were performed for typical 1000MW coal-­fired power plants (CFPP) that incurred an assumed 24 percent carbon capture energy penalty (EP), injected 90 percent of CO 2 emissions (~9 million mt-­ CO 2 injected annually), and treated extracted water onsite. Net present value (NPV), land requirements, laws and regulations, and technological limits were determined for each stage of disposal, and used to estimate feasibility. The boundary of the assessment began once extracted water was brought to the surface, and ended once the water evaporated, was injected underground, or was discharged into surface water bodies. Results of the assessment were generated, stored, and analyzed using Microsoft Excel spreadsheets and ESRI Geographical Information System (GIS) maps. Conclusions about the relative benefits and impacts of alternative brine-­management strategies were highly sensitive to local climate and weather, and aquifer water chemistry. The NPV of certain scenarios ranged from -­$50/mt-­CO 2 (a cost) to +$10/mt-­CO 2 (revenue). The land footprint of the scenarios in this study ranged from <1 km 2 to 100 km 2 . Brine extraction as a pressure management tool for GCS has potential for improving the economics and for minimizing the environmental impacts of CCS. In order to maximize this potential, careful analysis of each saline aquifer and region must be conducted to determine a regionally appropriate brine use sequence (BUS) at the time of site selection. Models that use GIS will be essential tools in determining such sequences for individual CFPP. Future studies that perform risk and life cycle assessments (LCA) of BUS scenarios, incorporate additional impact metrics into the BUS model, and enhance the temporal sensitivity of the model would improve the robustness of this regional assessment method.

54 ENVIRONMENTAL SCIENCES↗

Potosi Reservoir Modeling; History and Recommendations

As a part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO₂) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from two projects: the US DOE-funded Illinois Basin–Decatur Project being conducted by the Midwest Geological Sequestration Consortium in Macon County, Illinois, as well as data from the Illinois – Industrial Carbon Capture and Sequestration (IL-ICCS) project funded through the American Recovery and Reinvestment Act. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the Verification Well 1 (VW1) and the Injection Well (CCS1), structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for the Potosi Formation. The intention was for two million tonnes per annum (MTPA) of CO₂ to be injected for 20 years into the Potosi Formation. In 2013, updated reservoir models for the Cambrian Potosi Formation were evaluated. The data included formation tops from mud logs, well logs from the CCS1, VW1, and Verification Well 2 (VW2) wells, structural and stratigraphic formation from a larger 3D seismic survey, and field data from several waste water injection wells for Potosi Formation. The objective is to simulate the injection of CO₂ at a rate 3.5 million tons per annum (3.2 million tonnes per annum [MTPA]) for 30 years 106 million tons (96 MT total) into the Potosi Formation. The Potosi geomodeling efforts have evolved from using data from a single well in 2010 to the inclusion of data from three wells in 2013 which largely leverage the porosity and permeability logs plus knowledge of lost circulation zones. The first Potosi model (Potosi Geobody Model 2010) attempted to use the available seismic inversion data to inform the geomodel and predict vugular zones in advance of drilling VW1. Lost circulation zones in VW1 came in as the geologists predicted. The model was not implemented in subsequent simulation work. To date, reservoir models used for flow simulation work have relied predominately on Gaussian distributed properties (porosity and permeability) and have employed a single injection well. Potosi Model 2013b incorporated the new VW2 logs, and exhibited an extra level of sophistication by delineating the vugular intervals. This method added further realism that likely represents the best reservoir approximation to date. Where the 2010 reservoir models were 10 by 10 mi (16 by 16 km) in area, the 2013 models were expanded in size to 30 by 30 mi (48 by 48 km). The latest reservoir simulations show that a minimum of four injectors might be required to meet target injection rates. Still, there is data that requires further scrutiny and modeling methodologies that require testing for the Potosi Formation. This work is currently ongoing, and the next phase of the reservoir modeling intends to implement valuable data like porosity derived from seismic inversion, seismically derived geobodies, or a combination of both to further define vugular zones and the porosity distribution within the Potosi Formation. Understanding the dual porosity, dual permeability character of the Potosi remains the greatest challenge in representing this formation. Further analysis of the FMI* fullbore formation microimager data may aid in assessing this uncertainty. The Potosi Formation is indeed an interesting formation, and recommendations to further characterize it are included in the following list: - Data acquisition to identify the vugs permeability, distribution, and interconnectivity could be considered to perform a more rigorous evaluation of the Potosi Formation injectivity and capacity. This could be achieved by performing an injection test on a vugular interval to determine the vugs permeability, and an interference test between wells to evaluate the local vugs extent and interconnectivity. - A thorough study of the available FMI data may reveal specifics on estimating the vug to matrix ratio. This estimate could be used to further condition the porosity distribution. Porosity logs alone might underestimate the formation’s porosity associated with vugs. Porosity mapping derived from the seismic inversion could also be used in the succeeding task to characterize the lateral porosity distribution within the reservoir. This could involve the geobody methodology previously attempted in 2010. With or without seismic inversion porosity mapping, it is worth exploring whether increased lateral heterogeneity plays a significant role in Potosi injectivity. Investigations on vugular, dolomitic outcrops suggest that there may be significantly greater lateral heterogeneity than what has been modeled here. - The FMI data also reveals the presence of and helps describe open fractures. The presence of fractures will further enhance the formation’s permeability. The task of leveraging this data in the geomodeling effort still remains. Under the best of circumstances, this data describing open fractures may be combined with seismic attributes to delineate fracture corridors. Fracture modeling would certainly add another layer of sophistication to the model. Its contribution and applicability remain to be explored. - Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. - Currently there is no fracture gradient data available for the Potosi in the Decatur project area. The acquisition of the fracture pressure data could be considered to determine an appropriate maximum allowable bottomhole injection pressure. This would allow the evaluation of injectivity and the required number of wells in a more precise manner. - Special core analysis (SCAL) to determine the relative permeability and capillary pressure of the vugs and matrix could be considered to have a better estimation of the reservoir injectivity and plume extent. - Formation water sampling and analysis could be considered for the Potosi to estimate the water salinity and properties. A vertical flow performance evaluation could be considered for the succeeding task to determine the appropriate tubing size, the required injection tubing head pressure, and to investigate whether the corresponding well injection rate falls within the tubing erosional velocity limit. - A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.

Smith, Valerie↗

The Potosi Reservoir Model 2013

As a part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium in Macon County, Illinois. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the VW1 and the CCS1 wells, structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. Intention was for two million tons per annum (MTPA) of CO2 to be injected for 20 years. In the preceding, the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010" in this topical report) was re-run using a new injection scenario; 3.2 MTPA for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. It was not sufficiently large enough to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 miles x 30 miles (48.3km x48.3km), while preserving all property modeling workflows and layering. This model was retained as the base case of Potosi Dynamic Model 2013a. The Potosi reservoir model was updated to take into account the new data from the verification well VW2 which was drilled in 2012. The new porosity and permeability modeling was performed to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.2 MTPA for 30 years. This new dynamic model was named Potosi Dynamic Model 2013b. Due to the major uncertainties on the vugs permeability, two models were built; the Pessimistic and Optimistic Cases. The Optimistic Case assumes vugs permeability of 9,000 mD, which is analog to the vugs permeability identified in the pressure fall off test of a waste water injector in the Tuscola site, approx. 40 miles (64.4km) away from the IBDP area. The Pessimistic Case assumes that the vugs permeability is equal to the log data, which does not take into account the permeability from secondary porosity. The probability of such case is deemed low and could be treated as the worst case scenario, since the contribution of secondary porosity to the permeability is neglected and the loss circulation events might correspond to a much higher permeability. It is considered important, however, to identify the range of possible reservoir performance since there are no rigorous data available for the vugs permeability. The Optimistic Case gives an average CO2 injection rate of 0.8 MTPA and cumulative injection of 26 MT in 30 years, which corresponds to 27% of the injection target. The injection rate is approx. 3.2 MTPA in the first year as the well is injecting into the surrounding vugs, and declines rapidly to 0.8 MTPA in year 4 once the surrounding vugs are full and the CO2 start to reach the matrix. This implies that according to this preliminary model, a minimum of four (4) wells could be required to achieve the injection target. This result is lower than the injectivity estimated in the Potosi Dynamic Model 2013a (43 MT in 30 years), since the permeability model applied in the Potosi Dynamic Model 2013b is more conservative. This revision was deemed necessary to treat the uncertainty in a more appropriate manner. As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. For the Optimistic Case, the plume extends 17 miles (27.4km) in E-W and 14 miles (22.5km) in N-S directions after 30 years. After injection is completed, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 100 years post injection, the plume extends 20 miles (32.2km) in E-W and 15.5 miles (24.9km) in N-S directions. Should the targeted cumulative injection of 96 MT be achieved; a much larger plume extent could be expected. For the Optimistic Case, the increase of reservoir pressure at the end of injection is approximately 1200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 30 psia (206.8 kPa) beyond 14 miles (22.5km) away from injector. Should the targeted cumulative injection of 96 MT be achieved; a much larger areal pressure increase could be expected. The initial reservoir pressure is nearly restored after approximately 100 years post injection. The presence of matrix slows down the pressure dissipations. The Pessimistic Case gives an average CO2 injection rate of 0.2 MTPA and cumulative injection of 7 MT in 30 years, which corresponds to 7% of the injection target. This implies that in the worst case scenario, a minimum of sixteen (16) wells could be required to achieve the injection target. The present evaluation is mainly associated with uncertainty on the vugs permeability, distribution, and interconnectivity. The different results indicated by the Optimistic and Pessimistic Cases signify the importance of vugs permeability characterization. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Porosity mapping derived from the seismic inversion could also be used in the succeeding task to characterize the lateral porosity distribution within the reservoir. With or without seismic inversion porosity mapping, it is worth exploring whether increased lateral heterogeneity plays a significant role in Potosi injectivity. Investigations on vugular, dolomitic outcrops suggest that there may be significantly greater lateral heterogeneity than what has been modeled here. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗