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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Integrated characterization of CO 2 storage reservoirs on the Rock Springs Uplift combining geomechanics, geochemistry, and flow modeling

This project improves understanding of the effects CO 2 injection and storage on geomechanical, petrophysical, and other properties on two potential reservoirs, the Madison Limestone and Weber Sandstone, at the Rock Springs Uplift, southwest Wyoming. Rock physics studies, laboratory experiments, and computer simulations were employed. Two groups of rock samples were prepared, one aged with formation brine and the second aged with CO 2 -saturated formation brine and subsequently used for geomechanical, geochemical, and petrophysical testing. Rock physics studies provided a facies classification to discriminate among different lithologies and to investigate the feasibility of predicting petrophysical properties using seismic inversion. By applying a zone-wise regression relation established in this facies analysis, a 3D reservoir model was developed. It was determined that time-lapse seismic data are not feasible for monitoring CO 2 displacement at the Rock Springs Uplift or to mitigate leakage-risks situations because pressure and fluid saturation have small effects on elastic velocities and synthetic seismograms. The small variations that do occur are due to the low porosity of the sandstone and the high stiffness of the dolomite. In the Madison Limestone, mineral dissolution due to reaction with CO 2 -saturated brine was observed in highly permeable pathways while mineral precipitation happened in intergranular or micropores. Changes in pore size were quantified via Time-Domain Nuclear Magnetic Resonance transverse relaxation time (T 2 ) and diffusion coefficient distributions. As opposed to continuous injection of CO 2 -saturated brine, a subtle increase in porosity after static ‘soaking’ with CO 2 -saturated brine increases permeability due to the heterogeneity of changes in pore and pore throat sizes. Porosity-permeability relationships are fitted using a power-law function that can be adopted to model time-evolved porosity and permeability in a transport-limited CO 2 injection process. The effect of CO 2 -saturated brine on the linear elastic properties of Weber Sandstone and Madison Limestone cannot be generally concluded. Considering horizontal and vertical orientations and different confining pressures, no consistent relationship between stress (or strain) data and reaction with CO 2 -saturated brine was observed in the nonlinear plastic and post-failure regime. The change of elastic constants due to reaction with CO 2 -saturated brine is more significant in Madison Limestone than in Weber Sandstone; however, no consistent trend was observed. Reservoir simulations suggest that injectivity for Weber Sandstone is adversely impacted by low permeability. Changes in confining pressure do not impact injection rate and can be ignored while the effect of changes to pore pressure cannot. The effect of confining pressure and pore pressure on the size of the CO 2 plume in the Madison Limestone is negligible because of higher permeability. Geomechanical and petrophysical properties induced by CO 2 -brine-rock reactions were time-dependent. Simulations considered these time-dependent properties and indicate that injectivity, porosity distribution, pressure distribution and CO 2 saturation were affected.

58 GEOSCIENCES↗

DOE OSTI2023

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗

DOE OSTI2022

Experiments and Modeling of Proppant Embedment and Fracture Conductivity for the Caney Shale, Oklahoma, USA

ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.

Katende, A.↗

DOE OSTI2023

Produced Fluid Induced Mineralogy and Elemental Alterations of Caney Shale, Southern Oklahoma

ABSTRACT This study involves batch reactor experiments and subsequent analyses of samples from Caney Shale in the Ardmore Basin of South-Central Oklahoma. Samples include mainly rock cores and cuttings recovered from two wells respectively drilled vertically through and horizontally across the Caney Shale. Mineralogical compositions are obtained by X-Ray Diffraction (XRD) measurements whilst microstructure and elemental distribution are acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS) respectively. Batch experiments are then conducted using selected rock samples and produced fluid from the Caney Formation. Deionized water is also reacted with some samples to serve as standard. Experiments are conducted at 95°C and ambient pressure for 7 and 30 days to assess the geochemical rock-fluid interactions. Results show rock mineralogical compositions are predominantly quartz, feldspar, carbonates, and clay with minor pyrite. Post-experimental mineralogical changes observed in samples include increased amorphous entities especially within the clay portions of XRD plots and dissolution of feldspar and carbonate minerals and formation of new mineral phases, mostly clays and salts. These are corroborated by EDS elemental analyses which show decreased elemental compositions. The implications of reactions mentioned above include but not limited to, scale formation, clay fines migration and shale softening all of which pose significant permeability impairment on formation over time. INTRODUCTION Shale reservoirs account for a large share of unconventional reservoirs in the world (Lyu et al., 2015). However, ultra-low permeability and high clay compositions pose significant challenges when producing from these reservoirs (Dawuda and Srinivasan, 2022, 2023). Producing from these reservoirs therefore requires horizontal drilling and hydraulic fracturing technologies which have proven their efficacy in generating substantial permeability in reservoirs to ensure production (Fujian et al., 2019; Liu et al., 2018). Even after expensive horizontal drilling and hydraulic fracturing, geochemical reactions between engineered fluids and formation leads to fracture constriction and adversely impact petrophysical properties (permeability and porosity) of the reservoir. These technologies are therefore under constant development and improvement in various aspects to ensure fine tuning for specific reservoirs. Under present conditions, much of the hydrocarbon reserves in unconventional shale reservoirs are left unproduced due to rapid decline in permeability following resumption to production after hydraulic fracturing. It is therefore essential to understand the range of geochemical reactions that cause rapid depletion of permeability after hydraulic fracturing and apply these to each shale reservoir to ensure substantial recovery rates.

Awejori, G. A.↗

NASA NTRS1981

Porous tooling process for manufacture of graphite/polyimide composites

A porous tooling system was selected for the processing of Graphite/PMR-15 Polyimide laminates in thickness up to 3.2 mm. (0.125 inch). This tool system must have a reasonable strength, permeability dimensional stability, and thermal conductivity to accomplish curing at 600 F and 200 psi and 200 psi autoclave temperature and pressure. A permeability measuring apparatus was constructed and permeability vs. casting water level determined to produce tools at three different permeability levels. On these tools, laminates of 5, 11, and 22 plies (.027, .060, and 0.121 inch) were produced and evaluated by ultrasonic, mechanical, and thermal tests to determine the effect of the tool permeability on the cured laminates. All tools produced acceptable laminates at 5 and 11 plies but only the highest permeability produced acceptable clear ultrasonic C-Scans. Recommendations are made for future investigations of design geometry, and strengthening techniques for porous ceramic tooling.

Smiser, L. W.↗

NASA NTRS2003

Intracellular localization of a group II chaperonin indicates a membrane-related function

Chaperonins are protein complexes that are believed to function as part of a protein folding system in the cytoplasm of the cell. We observed, however, that the group II chaperonins known as rosettasomes in the hyperthermophilic archaeon Sulfolobus shibatae, are not cytoplasmic but membrane associated. This association was observed in cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C by using immunofluorescence microscopy and in thick sections of rapidly frozen cells grown at 76 degrees C by using immunogold electron microscopy. We observed that increased abundance of rosettasomes after heat shock correlated with decreased membrane permeability at lethal temperature (92 degrees C). This change in permeability was not seen in cells heat-shocked in the presence of the amino acid analogue azetidine 2-carboxylic acid, indicating functional protein synthesis influences permeability. Azetidine experiments also indicated that observed heat-induced changes in lipid composition in S. shibatae could not account for changes in membrane permeability. Rosettasomes purified from cultures grown at 60 degrees C and 76 degrees C or heat-shocked at 85 degrees C bind to liposomes made from either the bipolar tetraether lipids of Sulfolobus or a variety of artificial lipid mixtures. The presence of rosettasomes did not significantly change the transition temperature of liposomes, as indicated by differential scanning calorimetry, or the proton permeability of liposomes, as indicated by pyranine fluorescence. We propose that these group II chaperonins function as a structural element in the natural membrane based on their intracellular location, the correlation between their functional abundance and membrane permeability, and their potential distribution on the membrane surface.

Cell Membrane/physiology/ultrastructure↗

NASA NTRS2009

Amorphous Silk Fibroin Membranes for Separation of CO2

Amorphous silk fibroin has shown promise as a polymeric material derivable from natural sources for making membranes for use in removing CO2 from mixed-gas streams. For most applications of silk fibroin, for purposes other than gas separation, this material is used in its highly crystalline, nearly natural form because this form has uncommonly high tensile strength. However, the crystalline phase of silk fibroin is impermeable, making it necessary to convert the material to amorphous form to obtain the high permeability needed for gas separation. Accordingly, one aspect of the present development is a process for generating amorphous silk fibroin by treating native silk fibroin in an aqueous methanol/salt solution. The resulting material remains self-standing and can be prepared as thin film suitable for permeation testing. The permeability of this material by pure CO2 has been found to be highly improved, and its mixed-gas permeability has been found to exceed the mixed-gas permeabilities of several ultrahigh-CO2-permeable synthetic polymers. Only one of the synthetic polymers poly(trimethylsilylpropyne) [PTMSP] may be more highly permeable by CO2. PTMSP becomes unstable with time, whereas amorphous silk should not, although at the time of this reporting this has not been conclusively proven.

Aberg, Christopher M.↗

DOE OSTI2011

Application of real rock pore-threat statistics to a regular pore network model

This work reports the application of real rock statistical data to a previously developed regular pore network model in an attempt to produce an accurate simulation tool with low computational overhead. A core plug from the St. Peter Sandstone formation in Indiana was scanned with a high resolution micro CT scanner. The pore-throat statistics of the three-dimensional reconstructed rock were extracted and the distribution of the pore-throat sizes was applied to the regular pore network model. In order to keep the equivalent model regular, only the throat area or the throat radius was varied. Ten realizations of randomly distributed throat sizes were generated to simulate the drainage process and relative permeability was calculated and compared with the experimentally determined values of the original rock sample. The numerical and experimental procedures are explained in detail and the performance of the model in relation to the experimental data is discussed and analyzed. Petrophysical properties such as relative permeability are important in many applied fields such as production of petroleum fluids, enhanced oil recovery, carbon dioxide sequestration, ground water flow, etc. Relative permeability data are used for a wide range of conventional reservoir engineering calculations and in numerical reservoir simulation. Two-phase oil water relative permeability data are generated on the same core plug from both pore network model and experimental procedure. The shape and size of the relative permeability curves were compared and analyzed and good match has been observed for wetting phase relative permeability but for non-wetting phase, simulation results were found to be deviated from the experimental ones. Efforts to determine petrophysical properties of rocks using numerical techniques are to eliminate the necessity of regular core analysis, which can be time consuming and expensive. So a numerical technique is expected to be fast and to produce reliable results. In applied engineering, sometimes quick result with reasonable accuracy is acceptable than the more time consuming results. Present work is an effort to check the accuracy and validity of a previously developed pore network model for obtaining important petrophysical properties of rocks based on cutting-sized sample data.

58 GEOSCIENCES↗

DOE OSTI2011

Application of real rock pore-throat statistics to a regular pore network model

42 ENGINEERING↗

DOE OSTI2014

The Potosi Reservoir Model 2013c, Property Modeling Update

As 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 this project as well as two other separately funded projects: the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (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 Potosi Formation. The intention was for 2.2 million tons per annum (2 million tonnes per annum [MTPA]) of CO2 to be injected for 20 years. In the Task Error! Reference source not found., the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010") was re-run using a new injection scenario of 3.5 million tons per annum (3.2 MTPA) for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. The models size was insufficient 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 by 30 mi (48 by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. In the preceding Task [1], the Potosi reservoir model was updated to take into account the new data from the Verification Well #2 (VW2) which was drilled in 2012. The porosity and permeability modeling was revised 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.5 million tons per annum (3.2 MTPA) for 30 years. This dynamic model was named Potosi Dynamic Model 2013b. In this Task, a new property modeling workflow was applied, where seismic inversion data guided the porosity mapping and geobody extraction. The static reservoir model was fully guided by PorosityCube interpretations and derivations coupled with petrophysical logs from three wells. The two main assumptions are: porosity features in the PorosityCube that correlate with lost circulation zones represent vugular zones, and that these vugular zones are laterally continuous. Extrapolation was done carefully to populate the vugular facies and their corresponding properties outside the seismic footprint up to the boundary of the 30 by 30 mi (48 by 48 km) model. Dynamic simulations were also run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This new dynamic model was named Potosi Dynamic Model 2013c. Reservoir simulation with the latest model gives a cumulative injection of 43 million tons (39 MT) in 30 years with a single well, which corresponds to 40% of the injection target. The injection rate is approx. 3.2 MTPA in the first six months as the well is injecting into the surrounding vugs, and declines rapidly to 1.8 million tons per annum (1.6 MTPA) in year 3 once the surrounding vugs are full and the CO2 start to reach the matrix. After, the injection rate declines gradually to 1.2 million tons per annum (1.1 MTPA) in year 18 and stays constant. This implies that a minimum of three (3) wells could be required in the Potosi to reach the injection target. The injectivity evaluated in this Task was higher compared to the preceding Task, since the current facies modeling (guided by the porosity map from the seismic inversion) indicated a higher density of vugs within the vugular zones. 5 As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. After 30 years of injection, the plume extends 13.7 mi (22 km) in E-W and 9.7 mi (16 km) in N-S directions. After injection finishes, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 60 years post-injection, the plume extends 14.2 mi (22.8 km) in E-W and 10 mi (16 km) in N-S directions, and remains constant as the remaining pressure gradient has become very low. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger plume extent could be expected. The increase of reservoir pressure at the end of injection is approximately 1,200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 10 psia (69 kPa) beyond 14 mi (23 km) away from injector. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger areal pressure increase could be expected. The reservoir pressure declines rapidly during the first 30 years post injection and the initial reservoir pressure is nearly restored after 100 years post-injection. The present evaluation is mainly associated with uncertainty on the vugs permeability and interconnectivity. The use of porosity mapping from seismic inversion might have reduced the uncertainty on the lateral vugs body distributions. However, major uncertainties on the Potosi vugs permeability remains. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. 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 and appropriate spacing to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗

DOE OSTI2012

Review: Some low-frequency electrical methods for subsurface characterization and monitoring in hydrogeology

Low-frequency geoelectrical methods include mainly self-potential, resistivity, and induced polarization techniques, which have potential inmany environmental and hydrogeological applications. They provide complementary information to each other and to in-situ measurements. The self-potential method is a passive measurement of the electrical response associated with the in-situ generation of electrical current due to the flow of pore water in porous media, a salinity gradient, and/or the concentration of redoxactive species. Under some conditions, this method can be used to visualize groundwater flow, to determine permeability, and to detect preferential flow paths. Electrical resistivity is dependent on the water content, the temperature, the salinity of the pore water, and the clay content and mineralogy. Time-lapse resistivity can be used to assess the permeability and dispersivity distributions and to monitor contaminant plumes. Induced polarization characterizes the ability of rocks to reversibly store electrical energy. It can be used to image permeability and to monitor chemistry of the pore water-minerals interface. These geophysical methods, reviewed in this paper, should always be used in concert with additional in-situ measurements (e.g. in-situ pumping tests, chemical measurements of the pore water), for instance through joint inversion schemes, which is an area of fertile on-going research.

Revil, Andre↗

DOE OSTI2018

Effect of Magnetic Inclusions on the Effective Magnetostriction of Bulk Superconductors

A simple model is presented based on the Kim–Anderson model to further investigate the dependence of the effective magnetostriction of magnetic inclusion-superconducting matrix system on both the elastic and magnetic parameters including the elastic modulus, permeability, and volume fraction. The effect of the permeability on the magnetostriction is also obtained by implementing the continuity conditions of displacement and strain at the interface between the inclusion and the matrix through the magnetostriction loop. The results indicate that a stiffer inclusion can decrease the effective magnetostriction no matter whether the inclusion is magnetic or not and a larger effective magnetostriction can be obtained by choosing the matrix with a higher permeability, which gives an explanation about why the composite made from a matrix with a high permeability but a negligibly small magnetostriction yields unexpectedly low magnetostriction. Of particular interest is that in a certain range the effective magnetostriction of composites can be enhanced until it is saturated by increasing the permeability of matrix.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

Experimental and numerical analysis of parallel reactant flow and transverse mixing with mineral precipitation in homogeneous and heterogeneous porous media

Formation of mineral precipitates in the mixing interface between two reactant solutions flowing in parallel in porous media is governed by reactant mixing by diffusion and dispersion and is coupled to changes in porosity/permeability due to precipitation. The spatial and temporal distribution of mixing-dependent precipitation of barium sulfate in porous media was investigated with side-by-side injection of barium chloride and sodium sulfate solutions in thin rectangular flow cells packed with quartz sand. The results for homogeneous sand beds were compared to beds with higher or lower permeability inclusions positioned in the path of the mixing zone. In the homogeneous and high permeability inclusion experiments, BaSO 4 precipitate (barite) formed in a narrow deposit along the length and in the center of the solution–solution mixing zone even though dispersion was enhanced within, and downstream of, the high permeability inclusion. In the low permeability inclusion experiment, the deflected BaSO 4 precipitation zone broadened around one side and downstream of the inclusion and was observed to migrate laterally toward the sulfate solution. A continuum-scale fully coupled reactive transport model that simultaneously solves the nonlinear governing equations for fluid flow, transport of reactants and geochemical reactions was used to simulate the experiments and provide insight into mechanisms underlying the experimental observations. Lastly, migration of the precipitation zone in the low permeability inclusion experiment could be explained by the coupling effects among fluid flow, reactant transport and localized mineral precipitation reaction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Diffusion and Flow of CO{sub 2} in Carbon Anode for Aluminium Smelting

Using a well-made carbon anode with the right porosity characteristics is essential to the successful operation of a Hall–Héroult cell. Bubble generation at this carbon anode and its contribution to the overall voltage drop in aluminum production holds significant potential for reducing overpotential and improving energy savings. This voltage drop is believed to be greatly influenced by the gas diffusion in the anode carbon, for which there are a limited number of measurements to correlate with the carbon anode properties. In the present study, the CO{sub 2} gas diffusion characteristics were first investigated via measuring anode porosity using mercury intrusion porosimetry (MIP) for different samples with different permeability. Second, diffusion experiments were conducted by flowing CO{sub 2} gas into the anode sample at elevated temperatures (up to 960 °C). The impact of temperature, average pore size, and permeability on the diffusion coefficient was investigated. The diffusion coefficient was consequently calculated by curve fitting using diffusion theories. The value obtained varied from 1.38 × 10{sup −6} to 7.89 × 10{sup −6} m{sup 2}/s. A high diffusion coefficient was measured in the carbon sample with larger average pore size and higher permeability. It was found that anomalous diffusion behavior in the temperature range 600 °C to 960 °C was caused by convective flow effects, different diffusion mechanisms, and the Boudouard reaction for these carbon anode samples. Furthermore, by plotting the measured diffusion coefficient against average pore size and permeability, an increasing trend of diffusivity with an increase in pore size and permeability was observed.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Insights from in-situ X-ray computed tomography during axial impregnation of unidirectional fiber beds

In-situ X-ray computed tomography during axial impregnation of unidirectional fiber beds is used to study coupled effects of fluid velocity, fiber movement and preferred flow channeling on permeability. Here, in order to interpret the experimental measurements, a new computational tool for predicting axial permeability of very large 2D arrays of non-uniformly packed fibers is developed. The results show that, when the impregnation velocity is high, full saturation is attained behind the flow front and the fibers rearrange into a less uniform configuration with higher permeability. In contrast, when the velocity is low, fluid flows preferentially in the narrowest channels between fibers, yielding unsaturated permeabilities that are lower than those in the saturated state. Lastly, these insights combined with a new computational tool will enable improved prediction of permeability, ultimately for use in optimization of composite manufacturing via liquid impregnation.

36 MATERIALS SCIENCE↗

DOE OSTI2012

Characterization of Tight Gas Reservoir Pore Structure Using USANS/SANS and Gas Adsorption Analysis

Small-angle and ultra-small-angle neutron scattering (SANS and USANS) measurements were performed on samples from the Triassic Montney tight gas reservoir in Western Canada in order to determine the applicability of these techniques for characterizing the full pore size spectrum and to gain insight into the nature of the pore structure and its control on permeability. The subject tight gas reservoir consists of a finely laminated siltstone sequence; extensive cementation and moderate clay content are the primary causes of low permeability. SANS/USANS experiments run at ambient pressure and temperature conditions on lithologically-diverse sub-samples of three core plugs demonstrated that a broad pore size distribution could be interpreted from the data. Two interpretation methods were used to evaluate total porosity, pore size distribution and surface area and the results were compared to independent estimates derived from helium porosimetry (connected porosity) and low-pressure N{sub 2} and CO{sub 2} adsorption (accessible surface area and pore size distribution). The pore structure of the three samples as interpreted from SANS/USANS is fairly uniform, with small differences in the small-pore range (< 2000 {angstrom}), possibly related to differences in degree of cementation, and mineralogy, in particular clay content. Total porosity interpreted from USANS/SANS is similar to (but systematically higher than) helium porosities measured on the whole core plug. Both methods were used to estimate the percentage of open porosity expressed here as a ratio of connected porosity, as established from helium adsorption, to the total porosity, as estimated from SANS/USANS techniques. Open porosity appears to control permeability (determined using pressure and pulse-decay techniques), with the highest permeability sample also having the highest percentage of open porosity. Surface area, as calculated from low-pressure N{sub 2} and CO{sub 2} adsorption, is significantly less than surface area estimates from SANS/USANS, which is due in part to limited accessibility of the gases to all pores. The similarity between N{sub 2} and CO{sub 2}-accessible surface area suggests an absence of microporosity in these samples, which is in agreement with SANS analysis. A core gamma ray profile run on the same core from which the core plug samples were taken correlates to profile permeability measurements run on the slabbed core. This correlation is related to clay content, which possibly controls the percentage of open porosity. Continued study of these effects will prove useful in log-core calibration efforts for tight gas.

03 NATURAL GAS↗

DOE OSTI2013

Geomechanical effects on CO 2 leakage through fault zones during large-scale underground injection

The importance of geomechanics—including the potential for faults to reactivate during large-scale geologic carbon sequestration operations—has recently become more widely recognized. However, notwithstanding the potential for triggering notable (felt) seismic events, the potential for buoyancy-driven CO 2 to reach potable groundwater and the ground surface is actually more important from public safety and storage-efficiency perspectives. In this context, this paper extends the previous studies on the geomechanical modeling of fault responses during underground carbon dioxide injection, focusing on the short-term integrity of the sealing caprock, and hence on the potential for leakage of either brine or CO 2 to reach the shallow groundwater aquifers during active injection. We consider stress/strain-dependent permeability and study the leakage through the fault zone as its permeability changes during a reactivation, also causing seismicity. We analyze several scenarios related to the volume of CO 2 injected (and hence as a function of the overpressure), involving both minor and major faults, and analyze the profile risks of leakage for different stress/strain-permeability coupling functions. We conclude that whereas it is very difficult to predict how much fault permeability could change upon reactivation, this process can have a significant impact on the leakage rate. Moreover, our analysis shows that induced seismicity associated with fault reactivation may not necessarily open up a new flow path for leakage. Results show a poor correlation between magnitude and amount of fluid leakage, meaning that a single event is generally not enough to substantially change the permeability along the entire fault length. Finally, and consequently, even if some changes in permeability occur, this does not mean that the CO 2 will migrate up along the entire fault, breaking through the caprock to enter the overlying aquifer.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

A rock physics and seismic reservoir characterization study of the Rock Springs Uplift, a carbon dioxide sequestration site in Southwestern Wyoming

We present a reservoir geophysics study, including rock physics modeling and seismic inversion, of a carbon dioxide sequestration site in Southwestern Wyoming, namely the Rock Springs Uplift, and build a petrophysical model for the potential injection reservoirs for carbon dioxide sequestration. Our objectives include the facies classification and the estimation of the spatial model of porosity and permeability for two sequestration targets of interest, the Madison Limestone and the Weber Sandstone. The available dataset includes a complete set of well logs at the location of the borehole available in the area, a set of 110 core samples, and a seismic survey acquired in the area around the well. The proposed study includes a formation evaluation analysis and facies classification at the well location, the calibration of a rock physics model to link petrophysical properties and elastic attributes using well log data and core samples, the elastic inversion of the pre-stack seismic data, and the estimation of the reservoir model of facies, porosity and permeability conditioned by seismic inverted elastic attributes and well log data. In particular, the rock physics relations are facies-dependent and include granular media equations for clean and shaley sandstone, and inclusion models for the dolomitized limestone. The permeability model has been computed by applying a facies-dependent porosity-permeability relation calibrated using core sample measurements. Finally, the study shows that both formations show good storage capabilities. The Madison Limestone includes a homogeneous layer of high-porosity high-permeability dolomite; the Weber Sandstone is characterized by a lower average porosity but the layer is thicker than the Madison Limestone.

58 GEOSCIENCES↗