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 73 records · Page 4

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗

Effects of Supercritical CO2 on Matrix Permeability of Unconventional Formations

We studied the effects of supercritical carbon dioxide (scCO2) on the matrix permeability of reservoir rocks from the Eagle Ford, Utica, and Wolfcamp formations. We measured permeability using argon before exposure of the samples to scCO2 over time periods ranging from days to weeks. We measured permeability (and the change of permeability with confining pressure) when both argon and scCO2 were the pore fluids. In all three formations, we generally observe a negative correlation between initial permeability and carbonate content—the higher the carbonate content, the lower the initial permeability. In clay- and organic-rich samples, swelling of the matrix resulting from adsorption decreased the permeability by about 50% when the pore fluid was scCO2 although this permeability change is largely reversible. In carbonate-rich samples, dissolution of carbonate minerals by carbonic acid irreversibly increased matrix permeability, in some cases by more than one order of magnitude. This dissolution also increases the pressure dependence of permeability apparently due to enhanced mechanical compaction. Despite these trends, we observed no general correlation between mineralogy and the magnitude of the change in permeability with argon before and after exposure to scCO2. Flow of scCO2 through μm-scale cracks appears to play an important role in determining matrix permeability and the pressure dependence of permeability. Extended permeability measurements show that while adsorption is nearly instantaneous and reversible, dissolution is time-dependent, probably owing to reaction kinetics. Our results indicate that the composition and microstructure of matrix flow pathways control both the initial permeability and how permeability changes after interaction with scCO2. Electron microscopy images with Back-Scattered Electron (BSE) and Energy Dispersive Spectroscopy (EDS) revealed dissolution and etching of calcite minerals and precipitation of calcium sulfide resulting from exposure to scCO2.

58 GEOSCIENCES↗

Long-term permeability evolution of shale seal rocks with argon and scCO 2

In this work, we conducted a suite of experiments to evaluate the long-term permeability evolution of shale rocks under constant effective stress, before, during, and after interaction with supercritical carbon dioxide (scCO 2 ). To do so, we measured the time-dependent evolution of argon permeabilities for the pre- and post-reacted samples (i.e., before and after long-term interaction with scCO 2 ). In addition, we obtained permeability evolution during long-term interaction with scCO 2 . The samples showed either relatively constant permeabilities or a moderate decrease during pre reaction long-term argon tests. The permeability evolution during long-term CO 2 tests showed continuous increase, continuous decrease, or cycles of increase/decrease in permeability. The long term response of the samples to CO 2 included phenomena such as 20 (i) salt precipitation, (ii) swelling-induced cracks, and (iii) carbonate dissolution. While it is obvious that salt precipitation and swelling-induced cracks decrease and increase the permeability, respectively, the sample response to carbonate dissolution proved to be more complex and may increase or decrease the permeability. The permeability evolution during post-reaction long-term argon injection is also affected by the contribution of each of these three phenomena, during long-term interaction with CO 2 . We observe increase, decrease, and constant permeability evolution during post-reaction argon tests. Our experiments reveal that the initial permeability of the samples plays a significant role on the long-term permeability response of shales in the presence argon and scCO 2 fluids. This study shows, when shales are hydraulically fractured with CO 2 , their initial permeability has a more significant role than their permeability evolution over time.

58 GEOSCIENCES↗

Permeability evolution of shale during spontaneous imbibition

Shales have small pore and throat sizes ranging from nano to micron scales, low porosity and limited permeability. The poor permeability and complex pore connectivity of shales pose technical challenges to (a) understanding flow and transport mechanisms in such systems and, (b) in predicting permeability changes under dynamic saturation conditions. This paper presents quantitative experimental evidence of the migration of water through a generic shale core plug using micro CT imaging. In addition, in-situ measurements of gas permeability were performed during counter-current spontaneous imbibition of water in nano-darcy permeability Marcellus and Haynesville core plugs. It was seen that water blocks severely reduced the effective permeability of the core plugs, leading to losses of up to 99.5% of the initial permeability in experiments lasting 30 days. There was also evidence of clay swelling which further hindered gas flow. When results from this study were compared with similar counter-current gas permeability experiments reported in the literature, the initial (base) permeability of the rock was found to be a key factor in determining the time evolution of effective gas permeability during spontaneous imbibition. With time, a recovery of effective permeability was seen in the higher permeability rocks, while becoming progressively detrimental and irreversible in tighter rocks. Finally, these results suggest that matrix permeability of ultra-tight rocks is susceptible to water damage following hydraulic fracturing stimulation and, while shut-in/soaking time helps clearing-up fractures from resident fluid, its effect on the adjacent matrix permeability could be detrimental.

04 OIL SHALES AND TAR SANDS↗

Discontinuities in effective permeability due to fracture percolation

Motivated by a triaxial coreflood experiment with a sample of Utica shale where an abrupt jump in permeability was observed, possibly due to the creation of a percolating fracture network through the sample, we perform numerical simulations based on the experiment to characterize how the effective permeability of otherwise low-permeability porous media depends on fracture formation, connectivity, and the contrast between the fracture and matrix permeabilities. While a change in effective permeability due to fracture formation is expected, the dependence of its magnitude upon the contrast between the matrix permeability and fracture permeability and the fracture network structure is poorly characterized. We use two different high-fidelity fracture network models to characterize how effective permeability changes as percolation occurs. The first is a dynamic two-dimensional fracture propagation model designed to mimic the laboratory settings of the experiment. The second is a static three-dimensional discrete fracture network (DFN) model, whose fracture and network statistics are based on the fractured sample of Utica shale. Once the network connects the inflow and outflow boundaries, the effective permeability increases non-linearly with network density. In most networks considered, a jump in the effective permeability was observed when the embedded fracture network percolated. We characterize how the magnitude of the jump, should it occur, depends on the contrast between the fracture and matrix permeabilities. For small contrasts between the matrix and fracture permeabilities the change is insignificant. However, for larger contrasts, there is a substantial jump whose magnitude depends non-linearly on the difference between matrix and fracture permeabilities. A power-law relationship between the size of the jump and the difference between the matrix and fracture permeabilities is observed. In conclusion, the presented results underscore the importance of fracture network topology on the upscaled properties of the porous medium in which it is embedded.

36 MATERIALS SCIENCE↗

Relative permeability of gas and water flow in hydrate-bearing porous media: A micro-scale study by lattice Boltzmann simulation

The water-gas relative permeability is an important parameter to characterize multiphase flow in sediments. To study the water-gas relative permeability of hydrate-bearing porous media, multiphase flow simulations were carried out at the pore scale using the lattice Boltzmann method. In this work, the effects of hydrate saturation and hydrate-growth habits on the water-gas relative permeability, which is scaled by the relative permeability considering the hydrate only, were evaluated in a two-dimensional porous medium. Results show that the increase of hydrate saturation causes the decrease of water-gas effective permeability as expected. However, the effect of hydrate saturation on the water-gas relative permeability is different from that of hydrate saturation on the water-gas effective permeability. The water-gas relative permeability increases with the increase of hydrate saturation in the pore-filling case. The water-gas relative permeability decreases with the increase of hydrate saturation in the grain-coating case. The wettability of solid phase has a different effect on the relative permeability of wetting phase and nonwetting phase. The Jamin effect (phase blocking) was observed and may exist in the production of gas from natural gas hydrate reservoirs. This seriously affects the multiphase flow characteristics. The changes of microscale fluid distribution effect the changes of water-gas relative permeability. The relationship between the water-gas relative permeability and the characterization parameters of microscale fluid distribution was analyzed.

58 GEOSCIENCES↗

From core-scale experiment to reservoir-scale modeling: A scale-up approach to investigate reaction-induced permeability evolution of CO2 storage reservoir and caprock at a U.S. CO2 storage site

Mineral dissolution and secondary mineral precipitation can cause porosity and permeability changes of CO 2 storage reservoirs and caprocks after injection of CO 2 . In this paper, a 3-step approach (core-scale experiment →core-scale modeling →reservoir-scale modeling) is developed to simulate reservoir-scale porosity and permeability evolution of CO 2 storage formation and caprock at a model CO 2 storage site. The model site is based on characteristics of a real site in Mississippi, USA. Important chemical and permeability modeling parameters in the reservoir-scale model are validated by core-scale experimental and reactive transport modeling results. The reservoir-scale model predicts a maximum 3.2% permeability increase of the CO 2 storage formation and a maximum 1.1% permeability increase of the caprock after 1000 years of exposure to CO 2 -rich brine, while the core-scale model predicts 7% permeability decrease for a small CO 2 storage formation core and 296% permeability increase for a small caprock core after 180-day exposure to CO 2 -rich brine. The discrepancy between permeability results of reservoir-scale model and core-scale model is attributed to strong pH buffering effect of CO 2 storage formation with large mass of H+-consuming minerals. Therefore, using core-scale experiments/models only is not sufficient to elucidate reservoir-scale permeability evolution. Variations of key model parameters have a small effect on permeability evolution of both CO 2 storage formation and caprock, except for variations of K eq (SiO 2 (am)) and the exponent n in permeability-porosity correlation. SiO 2 (am) is a key mineral that governs permeability evolution of CO 2 storage formation and caprock, given the characteristics of the model CO 2 storage site.

58 GEOSCIENCES↗

Effect of medium permeability anisotropy on the morphological evolution of two non-uniformities in a geochemical dissolution system

The morphological evolutions of chemical dissolution fronts have attracted increasing interest in the field of the geological sciences and in industrial applications. Extensive research based on numerical simulations has been conducted to understand how various mechanisms and processes influence the morphological evolution of chemical dissolution fronts within geological media. Most researchers in previous studies have assumed the medium permeability to be isotropic for developing numerical models, despite isotropic geological media being uncommon in the real world. This study investigates the effect of medium permeability anisotropy on the morphological evolutions of two non-uniformities with higher permeability in a geochemical dissolution system. A series of numerical simulations are performed to evaluate the effect of medium permeability anisotropy on the morphological evolution of a chemical dissolution front. The simulation results indicate that the patterns of the dissolution reaction front are substantially affected by medium permeability anisotropy. An increase in the permeability anisotropy ratio, which is defined as the ratio of the permeability in the transverse direction to that in the longitudinal direction, enhances the dominance of the flow-focusing effect over the stabilizing or merging effect induced by diffusion/dispersion mechanism. Therefore, an increase in the permeability anisotropy ratio can increase the fingering length of the dissolution front or cause the dissolution front to have a more unstable pattern. By contrast, a reduction in the permeability anisotropy ratio will weaken the flow-focusing effect, thereby reducing the fingering length of the dissolution front or changing the front morphology such that it has a more stable status. The effect of the permeability anisotropy ratio on the morphological evolution tends to decrease when the Zhao number (negative dimensionless upstream pressure gradient) of the system increases. The consideration of medium permeability anisotropy in the geochemical dissolution model renders the simulation of the morphological evolutions of dissolution reaction fronts more realistic.

58 GEOSCIENCES↗

A Mechanistic Model for Relative Permeability of Gas and Water Flow in Hydrate-Bearing Porous Media With Capillarity

Lack of mechanistic models to describe petrophysical properties of mobile phases in gas hydrate–bearing sediments is one of the key challenges in accurately predicting gas production. The major drawback of empirical models that are used to fit a relative permeability curve for mobile phases in gas hydrate–bearing sediments is that they are based on experimental data, which are limited, and not able to account for hydrate morphology in pore space. This study proposes a relative permeability model that is mechanistic in nature and developed to account capillarity by building on the original nonempirical relative permeability model that assumes negligible capillary pressure. The proposed model implicitly accounts for capillarity with the help of four empirical parameters (two for each mobile phase) that incorporate each mobile phase pressure. It is shown that the proposed model provides an improved match to relative permeability data (derived from pore–scale simulation of gas hydrate–bearing sediments) than nonempirical relative permeability by accounting for the effect of capillary pressure. Additionally, unlike fully empirical models that can predict relative permeabilities reliably only at gas hydrate saturations for which experimental data are available, the proposed model only requires fitting the empirical parameters once with experimental data at any single gas hydrate saturation and then it can then be used to predict relative permeability at any gas hydrate saturation. Here, the mechanistic nature of the proposed model allows studying relative permeability of hydrate–bearing sediments as a function of hydrate morphology and wettability (fluid phase distribution) besides other physical parameters of the model (e.g., porosity, gas, and water residual saturation). Based on the sensitivity analysis of different hydrate morphologies on gas/water relative permeability, it is found that gas relative permeability is sensitive to hydrate morphologies, while water relative permeability shows little dependency on hydrate localization in pore space.

58 GEOSCIENCES↗

Porosity-Permeability Evolution During Simultaneous Mineral Dissolution and Precipitation

Permeability in reactive porous media evolves in complex ways that are hard to predict. Macroscopic empirical equations are often used to estimate permeability evolution but fail to reflect the impact of pore scale reactions on permeability. This work aims to understand the evolution of permeability in systems where mineral dissolution and precipitation occur simultaneously. Pore network models are used for permeability simulation of a Paluxy sandstone using pore and pore-throat size distributions from the analysis of X-ray CT images. Pore and pore-throat radii are increased and decreased to reflect the effects of dissolution and precipitation reactions. Varying spatial distributions of reactions are simulated and the resulting porosity and permeability values compared with commonly used macroscopic equations. It is observed that when dissolution occurs at the inlet and precipitation at the outlet, porosity increases while permeability decreases. Similar results are also observed in the respective opposite scenario where precipitation occurs at the inlet and dissolution at the outlet. When dissolution and precipitation reactions are randomly distributed throughout the network, porosity increases with little change in permeability that is not captured by the empirical equations. When reactions are controlled by pore and pore-throat sizes, such as dissolution occurring in pores and pore-throats of small size and precipitation in pores and pore-throats of larger size, porosity and permeability decrease. When precipitation occurs in pores and pore-throats of smaller size and dissolution in pores and pore-throats of larger size, porosity and permeability increase and is described relatively well by the Verma-Pruess porosity-permeability equation.

54 ENVIRONMENTAL SCIENCES↗

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

Coupling fracture facies with in-situ permeability measurements to generate stochastic simulations of tight carbonate aquifer properties: Example from the Lower Cretaceous aquifer, Northern Provence, SE France

The relationships between fracture facies and permeability distribution have been explored using a well-constrained 3D implicit structural model of hemipelagic low porosity/low permeability carbonate from the Northern Provence region, SE France. Fracture and permeability facies were determined using an extensive dataset of 99 hydrogeological wells. Data processing was undertaken using a step-by-step approach, involving: (i) identification of the fracture facies based on well logs detailing structure and stratigraphy; (ii) determination of permeability facies from the a priori correlation between the dimension of the hydraulic radius of influence (deduced from slug test analyses) and the type of reservoir heterogeneity (fissure, fracture, fault zone, etc.); (iii) three dimensional plot of fracture and permeability facies in the geological model using a variographic analysis of data. Thirty-three sequential indicator simulations (SIS) based on geostatistic analyses were realised on both fracture and permeability facies. Finally, a connectivity algorithm was developed to compute the probability of connection between selected infiltration areas and the major aquifer springs via moderate- to high-permeability geological bodies. Key results are: (i) fault zones have a greater role in controlling permeability facies distribution than on fracture facies repartition; (ii) there is little correlation between permeability and fracture facies distributions; (iii) connectivity results highlight the compartmentalisation of aquifers in the Cadarache area, the extensions of permeable geological bodies being limited by the N130 faults.

58 GEOSCIENCES↗

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗

Silver diamine fluoride differentially affects dentin and hypomineralized enamel permeabilities

OBJECTIVES: To investigate the physicochemical effect of silver diamine fluoride (SDF) by correlating permeability with mineral density and elemental composition of hypomineralized enamel and carious dentin. METHODS: Enamel and dentin from human carious primary teeth with and without SDF treatment in-vivo, and hypomineralized enamel from permanent molars with and without SDF treatment in-vitro were scanned using micro X-ray computed tomography. Spatial maps of biometals (calcium, zinc), phosphorus, and silver were generated using X-ray fluorescence microprobe. Permeabilities were computed using Porous Microstructure Analysis software. RESULTS: The intrinsic permeability of SDF-treated carious dentin was 14.3 % lower than untreated sound dentin (6.39e-15 ± 3.01e-15 m² vs 7.46e-15 ± 1.82e-15 m²; P < 0.0001), while untreated carious dentin was 98.4 % higher (1.48e-14 ± 7.11e-15 m²; P < 0.0001). SDF-treated and untreated transparent dentin showed similar reduced permeabilities (75.6 % and 78.4 % lower than untreated sound dentin, respectively; P = 0.93). Severely hypomineralized enamel showed permeability reaching 108.1 % of adjacent sound dentin (5.71e-15 ± 2.04e-15 m² vs 5.28e-15 ± 1.30e-15 m²; P = 0.1409) and was significantly higher than mildly hypomineralized enamel (1.39e-15 ± 1.04e-15 m²; P < 0.0001). SDF treatment did not significantly impact the permeability of severely hypomineralized enamel (12.4 % reduction; P = 0.07). Principal component regression identified Zn level as a significant effector of tissue permeabilities in carious primary teeth (P < 0.0001). SIGNIFICANCE: This study introduces a computational method to measure dental tissue permeability, and demonstrates that SDF significantly reduces permeability in carious dentin but not intact hypomineralized enamel. The study reveals biometal Zn localization can alter dentin and enamel permeabilities, providing new insights into pathobiological mechanisms underlying caries and hypomineralization.

Chou, Conrad↗