An official website of the United States government
Abstract The flow‐induced dissolution of porous rocks governs many important subsurface processes and applications. Solute mixing, which determines pore‐scale concentration fields, is a key process that affects dissolution. Despite its importance, the effects of pore‐scale mixing on large‐scale dissolution patterns have not been investigated. Here, we use a pore network model to elucidate the mixing effects on macroscopic dissolution patterns and solute transport. We consider two mixing rules at pore intersections that represent two end members in terms of the mixing intensity. We observe that the mixing effect on dissolution is the strongest at moderate Damköhler number, when the reactive and advective time scales are comparable. This is the regime where wormholes spontaneously appear. Incomplete mixing is shown to enhance flow focusing at the tips of the dissolution channels, which results in thinner wormholes and shorter breakthrough times. These effects on passive solute transport are evident independent of initial network heterogeneity.
more »
« less
Persistent Homology as a Heterogeneity Metric for Predicting Pore Size Change in Dissolving Carbonates
Abstract Accurate prediction of physical alterations in carbonate reservoirs under dissolution is critical for development of subsurface energy technologies. The impact of mineral dissolution on flow characteristics depends on the connectivity and tortuosity of the pore network. Persistent homology is a tool from algebraic topology that describes the size and connectivity of topological features. When applied to 3D X‐ray computed tomography (XCT) imagery of rock cores, it provides a novel metric of pore network heterogeneity. Prior works have demonstrated the efficacy of persistent homology in predicting flow properties in numerical simulations of flow through porous media. Its ability to combine size, spatial distribution, and connectivity information make it a promising tool for understanding reactive transport in complex pore networks, yet limited work has been done to apply persistence analysis to experimental studies on natural rocks. In this study, three limestone cores were imaged by XCT before and after acid‐driven dissolution flow‐through experiments. Each XCT scan was analyzed using persistent homology. In all three rocks, permeability increase was driven by the growth of large, connected pore bodies. The two most homogenous samples saw an increased effect nearer to the flow inlet, suggesting emerging preferential flow paths as the reaction front progresses. The most heterogeneous sample showed an increase in along‐core homogeneity during reaction. Variability of persistence showed moderate positive correlation with pore body size increase. Persistence heterogeneity analysis could be used to anticipate where greatest pore size evolution may occur in a reservoir targeted for subsurface development, improving confidence in project viability.
more »
« less
- Award ID(s):
- 1943726
- PAR ID:
- 10465056
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 59
- Issue:
- 9
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time‐resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid‐solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction‐front propagation and provide a foundation for controlling dissolution hotspots.more » « less
-
Pore-scale modeling is essential in understanding and predicting flow and transport properties of rocks. Generally, pore-scale modeling is dependent on imaging technologies such as Micro Computed Tomography (micro-CT), which provides visual confirmation into the pore microstructures of rocks at a representative scale. However, this technique is limited in the ability to provide high resolution images showing the pore-throats connecting pore bodies. Pore scale simulations of flow and transport properties of rocks are generally done on a single 3D pore microstructure image. As such, the simulated properties are only representative of the simulated pore-scale rock volume. These are the technological and computational limitations which we address here by using a stochastic pore-scale simulation approach. This approach consists of constructing hundreds of 3D pore microstructures of the same pore size distribution and overall porosity but different pore connectivity. The construction of the 3D pore microstructures incorporates the use of Mercury Injection Capillary Pressure (MICP) data to account for pore throat size distribution, and micro-CT images to account for pore body size distribution. The approach requires a small micro-CT image volume (7–19 mm3) to reveal key pore microstructural features that control flow and transport properties of highly heterogeneous rocks at the core-scale. Four carbonate rock samples were used to test the proposed approach. Permeability calculations from the introduced approach were validated by comparing laboratory measured permeability of rock cores and permeability estimated using five well-known core-scale empirical model equations. The results show that accounting for the stochastic connectivity of pores results in a probabilistic distribution of flow properties which can be used to upscale pore-scale simulated flow properties to the core-scale. The use of the introduced stochastic pore-scale simulation approach is more beneficial when there is a higher degree of heterogeneity in pore size distribution. This is shown to be the case with permeability and hydraulic tortuosity which are key controls of flow and transport processes in rocks.more » « less
-
How does the physical and chemical structure of the Critical Zone (CZ), defined as the zone from treetops to the bottom of groundwater, govern its hydro-biogeochemical functioning? Multiple lines of evidence from past and newly emerging research have prompted the shallow and deep partitioning concentration-discharge (C-Q) hypothesis. The hypothesis states that in-stream C-Q relationships are shaped by distinct source waters from flow paths at different depths. Base flows are often dominated by deep groundwater and mostly reflect groundwater chemistry, whereas high flows are often dominated by shallow soil water and thus mostly reflect soil water chemistry. The contrasts between shallow soil water versus deeper groundwater chemistry shape stream solute export patterns. In this context, the vertical connectivity that regulates the shallow and deep flow partitioning is essential in determining chemical contrasts, biogeochemical reaction rates in soils and parent rocks, and ultimately solute export patterns. This talk will highlight insights gleaned from multiple lines of recent studies that include collation of water chemistry data from soils, rocks, and streams in intensively monitored watersheds, meta-analysis of stream chemistry data at the continental scale, and integrated reactive transport modeling at the hillslope and watershed scales. The hypothesis underscores the importance of subsurface vertical structure and connectivity relative to the extensively studied horizontal connectivity. It also alludes to the potential of using streams as mirrors for subsurface water chemistry, and the potential of using C-Q relationships to infer flow paths and biogeochemical reaction rates and the response of earth’s subsurface to climate and human perturbations. Broadly, this simple conceptual framework links CZ subsurface structure to its functioning under diverse climate, geology, and land cover conditions.more » « less
-
Abstract Intermittent streams currently constitute >50% of the global river network, and the number of intermittent streams is expected to increase due to changes in land use and climate. Surface flows are known to expand and contract within the headwater channel network due to changes in the water table driven by climate, often changing seasonally. However, the underlying causes of disconnections and reconnections throughout the stream network remain poorly understood and may reflect subsurface flow capacity. We assess how 3D subsurface flowpaths control local surface flows at Gibson Jack Creek in the Rocky Mountains, Idaho, USA. Water table dynamics, hydraulic gradients, and hyporheic exchange were monitored along a 200‐m section of the stream throughout the seasonal recession in WY2018. Shallow lateral hillslope‐riparian‐stream connectivity was more frequent in transects spanning perennially flowing stream reaches than intermittent reaches. During low‐flow periods, larger losing vertical hydraulic gradients were observed in paired piezometers in intermittent reaches than in adjacent perennial reaches. Contrary to dominant conceptual models, longitudinal measurements of hydrologic exchange in both intermittent and perennial reaches were seasonally variable except for one perennial reach that showed consistent significant gains. Observed drying dynamics, as well as subsurface pathways, were highly variable even over short distances (30 m). Flow probability and subsurface flow capacity at upstream locations can be assessed with an outlet hydrograph and upstream flow measurements. Accurate characterization of subsurface storage, discharge, and connection is critical to understanding the drivers of drying cycles in intermittent streams and their likely responses to future change.more » « less
