Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permeability”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Coupled Hydrodynamic-Geomechanical Modelling of CO2-WAG Field Development at Farnsworth Unit: A Case Study

Abstract The SWP project is located in a mature waterflood undergoing conversion to CO2-WAG operations at Farnsworth, Texas, USA. Utilized CO2 is anthropogenic, sourced from a fertilizer and an ethanol plant. Major project goals are optimizing the storage/production balance, ensuring storage permanence, and developing best practices for CCUS. This paper provides a review of work performed toward development of a 3D coupled Mechanical Earth Model (MEM) for use in assessment of caprock integrity, fault reactivation potential, and evaluation of stress dependent permeability in reservoir forecasting. Mechanical property estimates computed from geophysical logs at selected wellbores were integrated with 3D seismic elastic inversion products to create a 3D "static" mechanical property model sharing the same geological framework as the existing reservoir simulation model including 3 major faults. Stresses in the MEM were initialized from wellbore stress estimates and reservoir simulation pore pressures. One way and two way coupled simulations were performed using a compositional hydrodynamic flow model and geomechanical solvers. Coupled simulations were performed on history matched primary, secondary (waterflood), and tertiary (CO2 WAG) recovery periods, as well as an optimized WAG prediction period. These simulations suggest that the field has been operating at conditions which are not conducive to either caprock failure or fault reactivation. Two way coupled simulations were performed in which permeability was periodically updated as a function of volumetric strain using the Kozeny-Carmen porosity-permeability relationship. These simulations illustrate the importance of frequent permeability updating when recovery scenarios result in large pressure changes such as in field re-pressurization through waterflood after a long primary depletion recovery period. Conversely, production forecasting results are less sensitive to permeability update frequency when pressure cycles are short and shallow as in WAG cycles. This paper describes initial work on development of a mechanical earth model for use in assessment of geomechanical risks associated with CCUS operations at FWU. The emphasis of this work is on integration of available geomechanical data for creation of the static mechanical property model. Preliminary coupled hydro-mechanical simulations are presented to illustrate some of the key diagnostic output from coupled simulations which will be used in later work for in depth evaluation of specific risk factors such as induced seismicity and caprock integrity.

02 PETROLEUM↗
DOE OSTI2022

Coreflood Tests to Evaluate Enhanced Oil Recovery Potential of Wettability-Altering Surfactants for Oil-Wet Heterogeneous Carbonate Reservoirs

Summary Oil-wetness and heterogeneity are two main factors that result in low oil recovery (OR) by waterflood in carbonate reservoirs. The injected water is likely to flow through high-permeability regions and bypass the oil in the low-permeability matrix. In this study, systematic coreflood tests were carried out in both “homogeneous” cores and “heterogeneous” cores with a wettability-altering surfactant. The homogeneous coreflood tests were conducted to evaluate surfactant retention, as well as to compare tertiary surfactant flooding with secondary surfactant flooding. The heterogeneous coreflood test was proposed to model bypassing in low-permeability matrix during waterfloods, and dynamic imbibition of surfactant into the low-permeability matrix. Surfactant retention results suggest that retention increases as initial oil saturation decreases. The retention of selected surfactant in the target reservoir cores was measured to be within a range of 0.07–0.12 mg/g-rock, which is economically acceptable. The results of homogeneous coreflood tests showed that both secondary waterflood and secondary surfactant flood can achieve high OR (>50%) from relatively homogeneous oil-wet cores. A shut-in phase after the surfactant injection resulted in a surge in oil production, which suggests that enough time should be given for wettability alteration by surfactants. The results of heterogeneous coreflood tests showed that more oil is bypassed in the tighter matrix by waterflood if the permeability is higher in the flooded layer and this bypassed oil is the target for the wettability-altering surfactant floods. Slow wettability-altering surfactant injection leads to imbibition into bypassed regions. When the oil-wet carbonate reservoirs have large unswept regions after waterflood, wettability-altering surfactants can significantly improve OR if enough time is given for imbibition.

Engineering↗
DOE OSTI2021

Experimental Study of Transport Behavior of Swellable Microgel Particles in Superpermeable Channels for Conformance Control

Summary Gel treatment is an effective way to attack excessive water production during oil development. The transport behavior of gel materials in reservoirs is of crucial importance to the effectiveness of gel treatments. The aim of this paper is investigating the transport behavior of swellable micrometer-sized preformed particle gels (PPGs, or microgels) through superpermeable (super-K) channels. Sandpacks with permeabilities ranging from 27 to 221 darcies were used to mimic the super-K channels. Multiple pressure sensors were applied along the sandpack models to monitor the propagation behavior of the microgels. The tested microgel particles could transport through the super-K channels, and a higher driving pressure gradient was required when the particle/pore size ratio was larger. The pressure gradient distribution along the super-K channels was relatively uniform when the particle/pore ratio was low (less than 1.3). However, the inlet section would show increasingly higher pressure gradients as the particle/pore ratio was increased, indicating increased difficulty in propagation. The propagation of the gel particles was significantly slower compared with the carrying fluid. The delayed propagation behavior was more pronounced when the particle/pore ratio was higher. The injection pressure was much less sensitive to the injection flow rate compared with a Newtonian fluid. The gel dispersion exhibited an apparent shear thinning (pseudoplastic) behavior when transporting through the porous channels. Breakage of the gel particles was observed especially at high superficial velocities. The particle breakage was partially responsible for the apparent shear thinning behavior. The breakage phenomenon was in favor of deep placement of the gel particles. The channel permeabilities were significantly reduced by the microgels, bringing sufficient resistance to subsequent waterflooding (more than 99.5%). At given matching size conditions, softer gels were more likely to establish in-depth placement and uniform water blocking capacity in the channels. The microgel particles exhibited salinity-responsive behavior to the post-brine flush. The gel particles could shrink and reswell according to the salinity of the injected water. Possibilities were discussed to use this salinity-responsive behavior. Also, the microgels exhibited a particular disproportionate permeability reduction (DPR) effect. After gel injection, the channel permeability to water flow was reduced by more than 20 to 92 times of the permeability to oil flow. This work provides important support to understand the transport behavior of gel particles in super-K channels. The achievements are helpful for gel product selection and gel treatment design.

Engineering↗
DOE OSTI2022

Impact of Polymer Rheology on Gel Treatment Performance of Horizontal Wells with Severe Channeling

Summary Gel treatment has been a cost-effective method to control the conformance of a reservoir with severe heterogeneity problems. The water channels in such reservoirs can be classified as open fractures or high permeability porous media with pore-throat network. Many simulation studies have been conducted to discuss gel treatment performance for conformance control. However, nobody considered the polymer rheology difference in open fractures and porous media in simulation. Previous simulation studies also ignored the residual resistance factor as a function of rock permeability rather than a constant parameter. In this study, a conceptual simulation model was established to simulate the linear flow system for the reservoir with horizontal wells considering the two factors mentioned above. The results demonstrate that the gel treatment always provides the better placement results in the open fracture type channel than pore-throat network type channel. Moreover, it is very necessary to consider residual resistance factor as a function of permeability, which is based on the experimental results and can provide much greater plugging efficiency in the higher permeable channels than constant residual resistance factor. Sensitivity analysis studies and multifactor analysis indicate that increasing oil viscosity and permeability ratio has a strong positive influence on conformance control results, which indicate in-situ gel treatment can be better applied in heavy or viscous oil reservoirs with fracture-like channels. Besides, the results also indicate that in reservoirs with severe channeling problem where channel velocity was high, the differences of gelant placement and profile improvement in models with two different types of channels could be enlarged greatly.

Engineering↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗
DOE OSTI2012

Effects of Porous Medium Heterogeneity on Vadose Zone Desiccation: Intermediate-scale Laboratory Experiments and Simulations

Soil desiccation (drying), involving water evaporation induced by dry gas injection, is a potentially robust vadose zone remediation process to limit contaminant transport through the vadose zone. A series of four intermediate-scale flow cell experiments was conducted in homogeneous and simple layered heterogeneous porous medium systems to investigate the effects of heterogeneity on desiccation of unsaturated porous media. The permeability ratios of porous medium layers ranged from about five to almost two orders of magnitude. The insulated flow cell was equipped with twenty humidity and temperature sensors and a dual-energy gamma system was used to determine water saturations at various times. The multiphase code STOMP was used to simulate the desiccation process. Results show that injected dry gas flowed predominantly in the higher permeability layer and delayed water removal from the lower permeability material. For the configurations tested, water vapor diffusion from the lower to the higher permeability zone was considerable over the duration of the experiments, resulting in much larger relative humidity values of the outgoing air than based on permeability ratios alone. Acceptable numerical matches with the experimental data were obtained when an extension of the saturation-capillary pressure relation below the residual water saturation was used. The agreements between numerical and experimental results suggest that the correct physics are implemented in the simulator and that the thermal and hydraulic properties of the porous media, flow cell wall and insulation materials were properly represented.

contaminant transport↗
DOE OSTI2012

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

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

02 PETROLEUM↗
DOE OSTI2014

Multiphase Fluid Flow in Deformable Variable-Aperture Fractures (Final Report, February 1, 2009 – January 31, 2014)

Fractures provide flow paths that can potentially lead to fast migration of fluids or contaminants. A number of energy-­related applications involve fluid injections that significantly perturb both the pressures and chemical composition of subsurface fluids. These perturbations can cause both mechanical deformation and chemical alteration of host rocks with potential for significant changes in permeability. In fractured rock subjected to coupled chemical and mechanical stresses, it can be difficult to predict the sign of permeability changes, let alone the magnitude. This project integrated experimental and computational studies to improve mechanistic understanding of these coupled processes and develop and test predictive models and monitoring techniques. The project involved three major components: (1) study of two-­phase flow processes involving mass transfer between phases and dissolution of minerals along fracture surfaces (Detwiler et al., 2009; Detwiler, 2010); (2) study of fracture dissolution in fractures subjected to normal stresses using experimental techniques (Ameli, et al., 2013; Elkhoury et al., 2013; Elkhoury et al., 2014) and newly developed computational models (Ameli, et al., 2014); (3) evaluation of electrical resistivity tomography (ERT) as a method to detect and quantify gas leakage through a fractured caprock (Breen et al., 2012; Lochbuhler et al., 2014). The project provided support for one PhD student (Dr. Pasha Ameli; 2009-­2013) and partially supported a post-­doctoral scholar (Dr. Jean Elkhoury; 2010-­2013). In addition, the project provided supplemental funding to support collaboration with Dr. Charles Carrigan at Lawrence Livermore National Laboratory in connection with (3) and supported one MS student (Stephen Breen; 2011-­2013). Major results from each component of the project include the following: (1) Mineral dissolution in fractures occupied by two fluid phases (e.g., oil-­water or water-­CO 2 ) causes changes in local capillary forces and redistribution of fluids. These coupled processes enhance channel formation and the potential for development of fast flow paths through fractures. (2) Dissolution in fractures subjected to normal stress can result in behaviors ranging from development of dissolution channels and rapid permeability increases to fracture healing and significant permeability decreases. The timescales associated with advective transport of dissolved ions in the fracture, mineral dissolution rates, and diffusion within the adjacent porous matrix dictate the sign and magnitude of the resulting permeability changes. Furthermore, a high--resolution mechanistic model that couples elastic deformation of contacts and aperture-­dependent dissolution rates predicts the range of observed behaviors reasonably well. (3) ERT has potential as a tool for monitoring gas leakage in deep formations. Using probabilistic inversion methods further enhances the results by providing uncertainty estimates of inverted parameters.

02 PETROLEUM↗
DOE OSTI2014

Investigation of Coupled Processes and Impact of High Temperature Limits in Argillite Rock

Shale and clay-rich geological formations have been considered as potential host rock for geological disposal of high-level radioactive waste throughout the world, because of their low permeability, low diffusion coefficient, high retention capacity for radionuclides, and capability to self-seal fractures. Permeabilities arising from the primary porosity of these rock types are low because of the extremely fine-grained constituents and the very small pores associated with the fine-grained material. Permeability values are typically less than 10 microdarcies (10-17 m2 or about 10-10 m/s equivalent hydraulic conductivity) and are often two to four orders of magnitude lower. While fractures can occur, clay and shale often demonstrate the tendency to self-seal fractures, which reduces the effects of fractures on bulk permeability. This occurs as a result of swelling from the increased water potential in fractures, the generation of fracture infilling materials and mineral precipitation, and plastic deformation of the rock (Mazurek et al., 2003). The low permeability of clay and shale rock are well-known in the hydrogeology community where these rock types represent aquitards that severely limit groundwater movement, and in petroleum geology, where they act as caprocks limiting the rise of buoyant petroleum fluids. Other favorable characteristics of clay/shale rock are the strong sorptive behavior for many radionuclides and (in saturated systems) the low flow rates, which typically lead to reducing conditions because of the lack of oxygen transport from the surface. Clay and shale rock also act to chemically buffer the effects of materials introduced through repository construction, operation, and emplaced materials (Arcos et al., 2008).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Quantification of Wellbore Leakage Risk Using Non-destructive Borehole Logging Techniques

Well integrity is important at all potential CCS locations and may play a crucial role establishing leakage risk in areas where there is a high density of existing wells that could be impacted by the storage operations including depleted petroleum fields where EOR or CCS will occur. To address a need for risk quantification methods that can be directly applied to individual wells using borehole logging tools a study was conducted using data from five wells in Wyoming. The objectives of the study were: Objective 1: Develop methods to establish the baseline flow parameters (porosity and permeability or mobility) from individual measurements of the material properties and defects in a well. Objective 2: Develop a correlation between field flow-property data and cement logs that can be used to establish the flow-properties of well materials and well features using cement mapping tools. Objective 3: Establish a method that uses the flow-property model (Objective 2) to analyze the statistical uncertainties associated with individual well leakage that can provide basis for uncertainty in risk calculations. The project objectives were met through the logging of five wells in Carbon and Natrona County Wyoming to collect data that was used to estimate individual and average well flow properties and model the results using ultrasonic data collected during the logging. Three of the five wells provided data on point and average flow properties for well annuli. Data from the other two wells were used to create models of cement permeability and test whether information collected in one well could be used to characterize another well. The results of the in-situ point measurements were confirmed by the lab measurements sidewall cores collected near the same depths Objective 1 was met using the data collected through logging, testing, and sampling. The methods were developed that can establish baseline flow parameters of wells by both point and average test methods. The methods to estimate the flow properties modeling of point pressure tests, modeling of vertical interference tests, and laboratory measurement of cased-hole sidewall cores The wells were in sufficiently good shape to allow the development of the characterization methods while still having enough defects to study differences in results as they relate to well integrity. Samples and tests analyzed from three of five wells studied in showed the cements were largely intact and had not degraded from exposure native brines. Log results taken in conjunction with the core measurements indicate that interfaces and/or problems with cement placement due to eccentering provide preferential flow paths for fluids, which can increase the effective permeability of the barrier several orders of magnitude above the permeability of intact cement. The results of the maps created using logging tools indicating that the cement condition and bond are generally good identify a need for more research to understand how logs can be used to predicteffective well permeabilities such as those measured by the VITs in this study.

42 ENGINEERING↗
DOE OSTI2015

Conducting a 3D Converted Shear Wave Project to Reduce Exploration Risk at Wister, CA

Ormat sited 2 full-size exploration wells based on 3D seismic interpretation of fractures, prior drilling results, and temperature anomaly. The wells indicated commercial temperatures (>300 F), but almost no permeability, despite one of the wells being drilled within 820 ft of an older exploration well with reported indications of permeability. Following completion of the second well in 2012, Ormat undertook a lengthy program to 1) evaluate the lack of observed permeability, 2) estimate the likelihood of finding permeability with additional drilling, and 3) estimate resource size based on an anticipated extent of permeability.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

The Coal-Seq III Consortium. Advancing the Science of CO 2 Sequestration in Coal Seam and Gas Shale Reservoirs

The Coal-Seq consortium is a government-industry collaborative that was initially launched in 2000 as a U.S. Department of Energy sponsored investigation into CO2 sequestration in deep, unmineable coal seams. The consortium’s objective aimed to advancing industry’s understanding of complex coalbed methane and gas shale reservoir behavior in the presence of multi-component gases via laboratory experiments, theoretical model development and field validation studies. Research from this collaborative effort was utilized to produce modules to enhance reservoir simulation and modeling capabilities to assess the technical and economic potential for CO2 storage and enhanced coalbed methane recovery in coal basins. Coal-Seq Phase 3 expands upon the learnings garnered from Phase 1 & 2, which has led to further investigation into refined model development related to multicomponent equations-of-state, sorption and diffusion behavior, geomechanical and permeability studies, technical and economic feasibility studies for major international coal basins the extension of the work to gas shale reservoirs, and continued global technology exchange. The first research objective assesses changes in coal and shale properties with exposure to CO2 under field replicated conditions. Results indicate that no significant weakening occurs when coal and shale were exposed to CO2, therefore, there was no need to account for mechanical weakening of coal due to the injection of CO2 for modeling. The second major research objective evaluates cleat, Cp, and matrix, Cm, swelling/shrinkage compressibility under field replicated conditions. The experimental studies found that both Cp and Cm vary due to changes in reservoir pressure during injection and depletion under field replicated conditions. Using laboratory data from this study, a compressibility model was developed to predict the pore-volume compressibility, Cp, and the matrix compressibility, Cm, of coal and shale, which was applied to modeling software to enhance model robustness. Research was also conducted to improve algorithms and generalized adsorption models to facilitate realistic simulation of CO2 sequestration in coal seams and shale gas reservoirs. The interaction among water and the adsorbed gases, carbon dioxide (CO2), methane (CH4), and nitrogen (N2) in coalbeds is examined using experimental in situ laboratory techniques to comprehensively model CBM production and CO2 sequestration in coals. An equation of state (EOS) module was developed which is capable of predicting the density of pure components and mixtures involving the wet CBM gases CH4, CO2, and N2 at typical reservoir condition, and is used to inform CO2 injection models. The final research objective examined the effects adsorbed CO2 has on coal strength and permeability. This research studied the weakening or failure of coal by the adsorption of CO2 from empirically derived gas production data to develop models for advanced modeling of permeability changes during CO2 sequestration. The results of this research effort have been used to construct a new and improved model for assessing changes in permeability of coal reservoirs due CO2 injection. The modules developed from these studies and knowledge learned are applied to field validation and basin assessment studies. These data were used to assess the flow and storage of CO2 in a shale reservoir, test newly developed code against large-scale projects, and conduct a basin-oriented review of coal storage potential in the San Juan Basin. The storage potential and flow of CO2 was modeled for shale sequestration of a proprietary Marcellus Shale horizontal gas production well using COMET3 simulation software. Simulation results from five model runs indicate that stored CO2 quantities are linked to the duration of primary production preceding injection. Matrix CO2 saturation is observed to increase in each shale zone after injection with an increase in primary production, and the size of the CO2 plume is also observed to increase in size the longer initial production is sustained. The simulation modules developed around the Coal-Seq experimental work are also incorporated into a pre-existing large-scale numerical simulation model of the Pump Canyon CO2-ECBM pilot in the San Juan Basin. The new model was applied to re-history match the data set to explore the improvements made in permeability prediction against previously published data sets and to validate this module. The assessment of the new data, however, indicates that the impact of the variable Cp is negligible on the overall behavior of the coal for CO2 storage purposes. Applying these new modules, the San Juan Basin and the Marcellus Shale are assessed for their technical ECBM/AGR and CO2 storage potential and the economic potential of these operations. The San Juan Basin was divided into 4 unique geographic zones based on production history, and the Marcellus was divided into nine. Each was assessed based upon each zone’s properties, and simulations were run to assess the potential of full Basin development. Models of a fully developed San Juan Basin suggest the potential for up to 104 Tcf of CO2 storage, and 12.3 Tcf of methane recovery. The Marcellus models suggest 1,248 Tcf of CO2 storage and 924 Tcf of AGR. The economics are deemed favorable where credits cover the cost of CO2 in the San Juan Basin, and in many cases in the Marcellus, but to maximize storage potential, credits need to extend to pay the operator to store CO2.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2017

A detailed study of CO 2 -brine capillary trapping mechanisms as applied to geologic carbon storage. Final Report

The proposed research focuses on improved fundamental understanding of the efficiency of physical trapping mechanisms, and as such will provide the basis for subsequent upscaling efforts. The overarching hypothesis of the proposed research is that capillary pressure plays a significant role in capillary trapping of CO 2 , especially during the water imbibition stage of the sequestration process. We posit that the relevant physics of the sequestration process is more complex than is currently captured in relative permeability models, which are often based on so-called trapping models to represent relative permeability hysteresis. Our 4 main questions, guiding the 4 main tasks of the proposed research, are as follows: (1) What is the morphology of capillary trapped CO 2 at the pore scale as a function of temperature, pressure, brine concentration, interfacial tension, and pore-space morphology under injection and subsequent imbibition? (2) Is it possible to describe the capillary trapping process using formation-dependent, but otherwise unique continuum-scale functions in permeability-capillary pressure, interfacial area and saturation space, rather than hysteretic functions in permeability-saturation or capillary pressure-saturation space? (3) How do continuum-scale relationships between kr-Pc-S-Anw developed based on pore-scale observations compare with traditional models incorporating relative permeability hysteresis (such as Land’s and other models,) and with observations at the core (5-10cm) scale? (4) How can trapped CO 2 volume be optimized via engineered injection and sweep strategies, and as a function of formation type (incl. heterogeneity)?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Innovative computational tools for reducing exploration risk through integration of water-rock interactions and magnetotelluric surveys

Mapping permeability distributions in geothermal reservoirs is essential for reducing the cost of geothermal development. To avoid the cost and sampling bias of measuring permeability directly through drilling, we require remote methods of imaging permeability such as geophysics. Electrical resistivity (or its inverse, conductivity) is one of the most sensitive geophysical properties known to reflect long range fluid interconnection and thus the likelihood of permeability. Perhaps the most widely applied geophysical methods for imaging subsurface resistivity is magnetotellurics (MT) due to its relatively great penetration depths. A primary goal of this project is to confirm through ground truthing at existing geothermal systems that MT resistivity structure interpreted integratively is capable of revealing permeable fluid pathways into geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Geomechanics-Based Stochastic Analysis of Injection- Induced Seismicity

The production of geothermal energy from dry and low permeability reservoirs is achieved by water circulation in natural and/or man-made fractures, and is referred to as enhanced or engineered geothermal systems (EGS). Often, the permeable zones have to be created by stimulation, a process which involves fracture initiation and/or activation of discontinuities such as faults and joints due to pore pressure and the in-situ stress perturbations. The stimulation of a rock mass is often accompanied by multiple microseismic events. Micro-seismic events associated with rock failure in shear, and shear slip on new or pre-existing fracture planes and possibly their propagations. The microseismic signals contain information about the sources of energy that can be used for understanding the hydraulic fracturing process and the created reservoir properties. Detection and interpretation of microseismic events is useful for estimating the stimulated zone, created reservoir permeability and fracture growth, and geometry of the geological structures and the in-situ stress state. The process commonly is referred to as seismicity-based reservoir characterization (SBRC). Although, progress has been made by scientific & geothermal communities for quantitative and qualitative analysis of reservoir stimulation using SBRC several key questions remain unresolved in the analysis of micro-seismicity namely, variation of seismic activity with injection rate, delayed micro-seismicity, and the relation of stimulated zone to the injected volume and its rate, and the resulting reservoir permeability. In addition, the current approach to SBRC does not consider the full range of relevant poroelastic and thermoelastic phenomena and neglects the uncertainty in rock properties and in-situ stress in the data inversion process. The objective of this research and technology developments was to develop a 3D SBRC model that addresses these shortcomings by taking into account hydro-thermo-poro-mechanical mechanisms associated with injection and utilizing a state-of-the-art stochastic inversion procedure. The approach proposed herein is innovative and significantly improves the existing SBCR technology (e.g., Shapiro et al. 2003) for geothermal reservoirs in several ways. First, the current scope of the SBRC is limited with respect to the physical processes considered and the rock properties used. Usually, the geomechanics analyses within SBRC is limited to the pore pressure diffusion in the rock mass, which is modeled using a time-dependent parabolic equation and solved using a finite element algorithm with either a line or a point source. However, water injection induces both poroelastic and thermoelastic stresses in the rock mass which affect the stress state. In fact, it has been suggested that thermoelastic stresses can play a dominant role in reservoir seismicity (Ghassemi et al., 2007). We include these important effects by using a fully-coupled poro-thermoelastic constitutive equations for the rock mass which will be solved using a 3D finite element model with more realistic injection geometries such as multiple injection/extraction sources (and in fractures), uncertainty in the material parameters and the in-situ stress distribution to better reflect the pore pressure and stress distributions. In addition, we developed a 3D stochastic fracture network model to study MEQ generation in fracture rocks. The model was verified using laboratory experiments, and calibrated and applied to Newberry EGS stimulation. In previous SBRC approaches, the triggering of micro-seismicity is modeled base on the assumption that the prior stochastic criticality model of the rock mass is a valid and adequate description. However, this assumption often does not hold in the field. Thus, we improved upon the current SBRC approach by using the micro-seismic responses to estimate the hydraulic diffusivity as well as the criticality distribution itself within the field. In this way, instead of relying on our a priori knowledge of criticality distribution, we combine an initial probabilistic description of criticality with the information contained in microseismic measurements to arrive at criticality solutions that are conditioned on both field data and our prior knowledge. Previous SBRC have relied upon a deterministic inversion approach to estimate the permeability, and the extent of the stimulated zone, whereas a stochastic inversion algorithm that recognizes and quantifies the uncertainties in the prior model, the time evolution of pore pressure distributions (modeling errors), and the observed seismic events is developed and used herein to realistically assess the quality of the solution. Finally, we developed a technique for processing discrete MEQ data to estimate fracture network properties such as dip and dip directions. The approach was successfully applied to the Fenton Hill HRD experiment and the Newberry EGS with results in good agreement with field observations.

15 GEOTHERMAL ENERGY↗