Articles | Volume 58
https://doi.org/10.5194/adgeo-58-1-2022
https://doi.org/10.5194/adgeo-58-1-2022
11 Oct 2022
 | 11 Oct 2022

Hysteresis in permeability evolution simulated for a sandstone by mineral precipitation and dissolution

Maria Wetzel, Thomas Kempka, and Michael Kühn
Abstract

Mineral dissolution and precipitation can substantially affect rock permeability, which is a critical parameter for a broad range of geological subsurface applications. Virtual experiments on digital pore-scale samples represent a powerful and flexible approach to understand the impact of microstructural alterations on evolving hydraulic rock behaviour and quantify trends in permeability.

In the present study, porosity-permeability relations are simulated for a precipitation-dissolution cycle within a typical reservoir sandstone. A hysteresis in permeability is observed depending on the geochemical process and dominating reaction regime, whereby permeability of the six investigated reaction paths varies by more than two orders of magnitude at a porosity of 17 %. Controlling parameters for this hysteresis phenomenon are the closure and re-opening of micro-scale flow channels, derived from changes in pore throat diameter and connectivity of the pore network. In general, a transport-limited regime exhibits a stronger impact on permeability than a reaction-limited regime, which uniformly alters the pore space. In case of mineral precipitation, higher permeability reduction results from successive clogging of pore throats, whereas in case of dissolution, permeability significantly increases due to a widening of existing flow paths. Both, the geochemical process and dominating reaction regime govern characteristic microstructural alterations, which cannot be simply reversed by the inversion of the geochemical processes itself. Hence, permeability evolution clearly depends on the hydrogeochemical history of the sample.

1 Introduction

Permeability is a crucial hydrogeological parameter for the characterisation of fluid flow in porous media and highly relevant for a broad range of geological subsurface applications, including geothermal energy recovery, carbon dioxide, natural gas and hydrogen storage or hydrocarbon production, as well as subsurface contaminant transport (De Lucia et al.2011; Jacquey et al.2015; Hofmann et al.2014; Kleinitz et al.2001; Regenspurg et al.2015). Analytical or empirical methods are commonly applied in reservoir models to estimate evolving permeability trends. These approaches generally describe permeability as a function of porosity (Bernabé et al.2003; Ma2015), such as the Kozeny-Carman relationship (Carman1937) or the power law (Hommel et al.2018). However, their predictive capabilities are limited, since they have been originally developed for idealised microstructures. In general, mineral precipitation reduces the porosity of a rock (Fig. 1), but its effect on transport properties strongly depends on the specific spatial alterations of the microstructure. Thus, porosity-permeability relationships can be more complex especially for geochemical processes, because mineral precipitation and/or dissolution locally change the pore structure, and thus the flow field, e.g. by closing or opening preferential flow paths. The resulting non-uniform permeability response cannot be captured by common analytical approaches (Beckingham2017; Lesueur et al.2020).

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f01

Figure 1The geological subsurface is of essential relevance for a broad range of applications. Geochemical processes alter the microstructure of rocks, which in turn affects their hydraulic behaviour on the macro scale, and thereby governs the evolution of geological subsurface systems.

Download

Numerical simulations on pore-scale samples represent a flexible method to determine fundamental rock property relations (Berg et al.2017). In this context, the research field of digital rock physics comprises the numerical simulation of hydraulic (Blunt et al.2013; Herring et al.2017), elastic (Andrä et al.2013; Shulakova et al.2017), electric (Keehm et al.2001; Wu et al.2020) and/or thermal properties (Ettemeyer et al.2020; Schepp et al.2020), which are directly calculated on digital samples, mainly derived from micro-CT scans. These virtual experiments offer a non-destructive method with direct access to rock parameters, which is further applied to compute evolving trends in rock properties based on microstructural changes (Niu and Zhang2019; Noiriel et al.2016; Sain et al.2016; Wetzel et al.2020b).

Within the presented study, a cycle of secondary mineral precipitation and subsequent dissolution is simulated for a synthetic sample of a Fontainebleau sandstone. Trends in permeability are computed for two contrasting geochemical reaction regimes: (1) a transport-limited and (2) a surface reaction-limited regime. Moreover, morphometric parameters of the pore space are computed and analysed to quantify arising changes in the microstructure and understand hydraulic impacts depending on the prevailing geochemical processes and reaction regimes.

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f02

Figure 2(a) Virtual sample of the unaltered synthetic sandstone including the flow field, represented by streamlines coloured by the local fluid flow velocity. Frequency distribution of the (b) pore throat size and (c) connectivity of the pore network with their respective average value.

Download

2 Materials and methods

In the present study, a digital sample of a highly porous sandstone is synthetically generated to allow for an assessment of the reaction-induced permeability evolution over a preferably wide range of porosities. Therefore, a process-based approach is applied, combining the gravity-driven deposition of irregularly shaped grains, and their diagenetic cementation. A more detailed description of this concept is provided in Wetzel et al. (2021), demonstrating that the constructed homogenous sandstone samples show microstructural and physical properties comparable to their natural equivalents. Grain size distribution and grain shape of a Fontainebleau sandstone are applied to generate the synthetic sample, since it represents a comprehensively examined reference rock for reservoir engineering applications (Gomez et al.2010; Walderhaug et al.2012). The virtual sandstone investigated here has a resolution of 4503 voxels with an edge length of 5 µm, whereby the total porosity of the unaltered sample is 24.7 % (Fig. 2a).

The morphometric parameters of the pore space are calculated using the Python package PoreSpy (Gostick et al.2019). In order to quantify changes within the pore space, the individual pores are initially extracted for the unaltered microstructure by means of a watershed segmentation. Then, labelling of the pores is maintained and overlaid by the altered microstructure, so that variations in pore morphology due to the simulated alterations can be calculated. Thus, median pore throat diameter and mean connectivity of the pore network are determined, whereby the latter represents the number of throats connected to a single pore (Fig. 2b and c).

Permeability is calculated from the flow field by solving the steady-state Stokes equation using the OpenFOAM software package (Weller et al.1998). The virtual experiment simulates a cycle of (1) precipitation of an arbitrary secondary mineral until the sample is nearly impermeable, and (2) the subsequent dissolution of the previously precipitated mineral. Corresponding trends in permeability are further determined for two contrasting geochemical reaction regimes, presupposing that the system is generally advection‐dominated (high Péclet numbers). In case of a surface reaction-limited process, the chemical reaction rate is the limiting factor, and the pore space is altered uniformly around the grains (Damköhler number ≪1). For the transport-limited alteration, the chemical reaction is limited by the availability of reactants, which are transported to the fluid-rock interface (Damköhler number ≫1). In this case, the local velocity is used as a proxy for solvent flux, and precipitation as well as dissolution are defined to occur in regions of high flow velocities (>75th percentile). A detailed description of this iterative approach can be found in Wetzel et al. (2020a). The aim of this simplification is to approximate characteristic precipitation and dissolution patterns instead of performing fully-coupled reactive transport simulations. Thus, a bandwidth of hydraulic property changes can be investigated for, e.g., sandstones or carbonates (Niu and Zhang2019; Miller et al.2017), whereby the simulated porosity-permeability curves are comparable to experimentally determined relations (Wetzel2021).

3 Results

3.1 Precipitation paths

The precipitation of an arbitrary secondary mineral by the two contrasting geochemical reaction regimes leads to distinctive alteration patterns of the pore space. The surface reaction-limited precipitation uniformly alters the pore structure. All pores are equally affected, which illustrates the continuous reduction in pore throat diameter and connectivity of the pore network. At a porosity of 6 %, the virtual sample becomes nearly impermeable and precipitation is stopped. The resulting decrease in permeability for the reaction-limited regime can be approximated by the power law with an exponent of four (Fig. 3a).

By contrast, the transport-limited regime shows by definition a localised precipitation in areas of high fluid flow velocities. This leads to a preferential clogging of pore throats, indicated by the morphometric evolution of the pore space (Fig. 3b): (I) First, fluid flow paths become narrower, since the median pore throat diameter is strongly reduced, while simultaneously the connectivity of the pore network is changed only slightly. (II) Successive mineral growth leads to a closure of fluid flow paths, whereby connectivity of the pore network is drastically reduced, while the throat diameter is less affected. The described morphological alteration of the pore space is the reason for the steep permeability decrease, whereby the porosity-permeability relation can be described by the power law with an exponent above ten. Consequently, low permeabilities comparable to the reaction-limited alteration are reached at a considerably higher porosity of 15 %, and the amount of precipitated minerals is significantly lower.

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f03

Figure 3(a) Simulated porosity-permeability trends for reaction- and transport-limited precipitation regimes compared to analytical relations based on the power law with different exponents (η). (b) Evolution of median throat diameter and mean connectivity for both chemical regimes with the characteristic phases of (I) flow path narrowing and subsequent their complete clogging (II).

Download

3.2 Dissolution paths

The subsequent mineral dissolution also shows characteristic morphometric alterations of the pore space, depending on the dominant geochemical reaction regime. The initial microstructure and permeability (Fig. 4a) will be achieved by definition, since the quartz grains are considered as inert and only the previously precipitated minerals can be dissolved. In case of reaction-limited dissolution, (I) the pore structure is initially uniformly affected, so that the pore throat diameter and connectivity of the pore network constantly decrease (Fig. 4c). Since the fluid-mineral contact area is substantially lower in narrow throats, which are filled with the previously precipitated minerals, it takes additional iterations for their uniform dissolution. Hence, a substantial amount of 53 % of the pore throats is still filled with secondary minerals at the high porosities of 19 % (Fig. 4d). Thus, (II) dissolution of these narrow throats leads to a drastic increase in connectivity of the pore network, and simultaneously reduces the median throat diameter. Nevertheless, the permeability curve of the sample is still continuous and does not exhibit any sharp changes. The porosity-permeability relation of the reaction-limited dissolution regime can be approximated by the power law with an exponent of five.

The transport-limited mineral dissolution (I) initially exhibits a widening of the existing main flow paths, indicated by a sharp increase in the throat diameter, while the pore network connectivity remains nearly constant (Fig. 4c). Fluid flow velocities substantially increase in these flow paths, since flow is channelled and all of the previously precipitated minerals in the respective pores are dissolved (Fig. 4d). (II) Successive mineral dissolution leads to the development of new flow paths, where connectivity of the pore network constantly increases and throat diameters slightly widen. (III) As for the reaction-limited regime, secondary minerals in narrow throats exhibiting a lower fluid-to-mineral interface are dissolved, what demonstrates the decrease in the median throat diameter resulting from the increasing amount of newly dissolved small throats. The previously described steps of the morphological alteration result in a non-uniform porosity-permeability relation, comprising a steep initial increase due to the fluid flow path widening, followed by a comparably lower permeability increase until the initial state is reached again (Fig. 4a).

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f04

Figure 4(a) Simulated porosity-permeability trends for reaction- and transport-limited dissolution compared to analytical power law relations with different exponents (η). Evolution of median throat diameter and mean connectivity for the (b) transport- and (c) reaction-limited regimes with phases of characteristic morphometric evolution. (d) Differences in the flow field between both dissolution regimes represented by streamlines (xy projection of the 3D sample) as well as the microstructure (2D slice).

Download

3.3 Permeability hysteresis of the precipitation-dissolution cycle

The precipitation-dissolution cycle under the assumption of two contrasting geochemical reaction regimes was investigated by means of six simulation scenarios with different porosity-permeability relationships (Fig. 5a). A hysteresis in the evolution of hydraulic properties depending on the geochemical process and dominating reaction regime can be observed. Controlling parameters for this phenomenon are characteristic changes of the pore space, as the closing or opening of flow channels described in Sect. 3.1 and 3.2.

Reaction-limited precipitation (Scenario 1) and subsequent dissolution (Scenario 1b) uniformly alter the pore structure, and thus largely show continuous changes in pore throat diameters and connectivity (Fig. 5b). Thus, permeability is constantly changed in both cases. Permeabilities of the dissolution path (Scenario 1b) are consistently lower due to the smaller fluid-to-mineral interface of the narrow throats, filled with the previously precipitated minerals. Also for Scenario 2b, the reaction-limited dissolution mostly leads to uniform changes in throat diameters and pore network connectivity (Fig. 5b). Nevertheless, the parameter evolution depends on the initial distribution of secondary minerals. Since previous precipitation preferentially occurs in pore throats (Fig. 6b), uniform dissolution leads to a comparatively gentle increase in the hydraulic properties, exhibiting the lowest permeabilities of all investigated scenarios (Fig. 5a).

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f05

Figure 5(a) Porosity-permeability curves for the six investigated reaction paths of the simulated precipitation-dissolution cycle. (b) Characteristic evolution of the median throat diameter (green) and mean connectivity (blue) depend on the geochemical process and dominating reaction regime. At a porosity of 17 %, permeability varies by more than two orders of magnitude.

Download

https://adgeo.copernicus.org/articles/58/1/2022/adgeo-58-1-2022-f06

Figure 6(a) Characteristic differences in the fluid flow regime represented by streamlines (xy projection of the 3D sample) at a porosity of 17.4 %, where permeability varies by more than two orders in magnitude. (b) Comparison of characteristic morphometric differences and the local flow velocities within the pore space for the six simulation scenarios (1a–2b).

Download

Transport-limited geochemical processes exhibit characteristic non-uniform alteration patterns, where flow paths are systematically clogged or dissolved. The narrowing of pore throats and their successive closure are typical for the precipitation processes occurring in this reaction regime (Scenario 2) and result in a comparably steep permeability decrease (Fig. 5a). The transport-limited dissolution of minerals (Scenarios 1a and 2a) in turn initially lead to a predominant widening of existing fluid flow paths and force channelised flow. This exhibits the sharp increase in the throat diameter, while connectivity remains nearly constant, resulting in a steep porosity-permeability curve (Fig. 5b). Successive mineral dissolution then induces new fluid flow paths, what increases the pore network connectivity and leads to a less pronounced increase in permeability. In general, this characteristic evolution can be observed for both transport-limited dissolution scenarios, but the intensity of each phase depends on the spatial distribution pattern of the previously precipitated minerals. In Scenario 2a, these dissolution steps are less prominent, since secondary minerals are initially localised around pore throats, what reduces the effect of flow path widening and the related intensity of the permeability increase (Fig. 5a and b).

Thus, a permeability hysteresis depending on the geochemical process and dominating reaction regime can be observed. In general, reaction-limited precipitation results in higher permeabilities than transport-limited precipitation. Analogously, transport-limited dissolution processes exhibit consistently higher permeabilities than the respective precipitation path and the reaction-limited dissolution (Fig. 5a). For the same porosity of 17 %, permeability varies by more than two orders in magnitude between the six simulated scenarios. Reaction-limited precipitation followed by a transport-limited dissolution (Scenario 1a) exhibits the highest permeability of 2.9 Darcy (1 Darcy =9.869×10-13 m2) due to the initially uniform precipitation and fluid flow channelling, resulting from preferential mineral dissolution inside the main fluid flow paths (Fig. 6a and b). Lowest permeabilities of 30 Millidarcy occur in Scenario 2b, since the previous transport-limited precipitation preferentially clogged pore throats, while the uniform reaction-limited dissolution requires more iterations to extend narrow throats due to the smaller fluid-to-mineral interface (Fig. 6b). Hence, the results clearly illustrate that permeability evolution depends on the geochemical history of a rock sample, where both the geochemical process and dominating reaction regime govern characteristic microstructural alterations within the pore space. These cannot be simply reversed by the inversion of the geochemical processes itself, so that the respective paths require specific assessments.

4 Discussion

The observed permeability hysteresis of the investigated geochemical cycle of secondary mineral precipitation and its subsequent dissolution can be explained by characteristic morphometric alterations of the pore space. Controlling parameters for the observed changes in hydraulic parameters are the clogging and re-opening of pore-scale fluid flow channels, which can be quantified by the pore throat diameter and connectivity of the pore network. The observed predominant clogging of pore throats as well as the related drastic permeability decrease in the transport-limited precipitation regime have been reported by several studies (Tenthorey and Scholz2002; Crandell et al.2012; Niu and Zhang2019). Also the selective widening of fluid flow paths and a flow channelling by preferential pathways are well aligned with observations from experiments, especially for evolving wormholes in carbonatic rocks (Lebedev et al.2017; Luquot et al.2016; Menke et al.2014; Miller et al.2017), and have been also observed for strongly cemented sandstones (Kühn2004).

Lesueur et al. (2020) modelled a dissolution-precipitation cycle by uniformly altering digital rock samples and reported a persisting permeability increase. For a sandstone sample, they demonstrated the presence of a hysteresis phenomenon, where the final permeability is less than 3 % above the initial one. In contrast, a permeability increase of 110 % has been observed for a carbonate sample. In the models presented in this study, the initial permeability is achieved by definition, since only the previously precipitated minerals can be dissolved, while the granular structure is considered as inert. Nevertheless, there is a clear tendency showing that reaction-limited dissolution leads to lower permeabilities compared to the previous precipitation path, since dissolution cannot restore all pore thoats previously filled by mineral precipitation. This is in agreement with the work of Lesueur et al. (2020), where for the reverse cycle of dissolution followed by precipitation permeabilities increase, since dissolved new fluid flow channels remained open after mineral precipitation. Hence, considering a reaction-limited regime, a dissolution-precipitation cycle would result in consistently higher permeabilities compared to the initial state, whereas a precipitation-dissolution cycle would induce a decrease in permeability. Moreover, even if not investigated here, the hysteresis phenomenon and associated permeability changes will be more pronounced in carbonates or rocks in which the entire microstructure can be chemically altered. For example, the bend within the porosity-permeability curve due to the change from flow path widening to the development of new flow paths is associated with the underlying microstructure of the transport-limited dissolution Scenarios 1a and 2a. Since only secondary minerals within an inert granular structure are allowed to be dissolved, the remaining flow paths will be dissolved by definition. In contrast to this, flow path widening and flow channelling would continue in carbonates.

Nevertheless, the investigated geochemical reaction regimes are end-members and reactive transport is more complex and not exclusively depending on fluid flow velocity. Chemical reactions are controlled by various other factors like fluid chemistry, mineralogy, temperature, pore morphology, pore and lithological pressures as well as transport properties (Cil et al.2017; Fazeli et al.2020; Schepers and Milsch2013; Beckingham2017). Moreover, the presented approach presumes advection-dominated species transport and implicitly neglects diffusion. This is considered as a valid assumption for many applications in geological subsurface systems. Despite the fact that the field of physical pore-scale simulations is quickly evolving due to increasing computational efficiency, there are up to now no studies known which explicitly perform reactive transport simulations for mineral precipitation and dissolution on realistic three-dimensional pore space models at the given spatial discretisation. In this regard, the presented approach and investigated scenarios represent a flexible and moreover feasible method to quantify permeability trends and support the understanding of the hydraulic impact of microstructural alterations.

5 Conclusions

In this study, porosity-permeability trends are computed for a precipitation-dissolution cycle within a typical reservoir sandstone. Virtual experiments are performed, simulating secondary mineral precipitation within the pore space until the sample is nearly impermeable, and the subsequent dissolution of the previously precipitated minerals. Two contrasting geochemical reaction regimes are assessed, which govern the preferential location of the alteration in the pore space.

The resulting six scenarios of the investigated precipitation-dissolution cycle clearly show a permeability hysteresis, depending on the geochemical process and dominating reaction regime. Hereby, permeability varies by more than two orders of magnitude for the same porosity of 17 %. Controlling parameters for this hysteresis phenomenon are the closure and re-opening of pore-scale fluid flow channels, which can be quantified by the morphometric parameters pore throat diameter and connectivity of the pore network. In general, a transport-limited regime exhibits a stronger impact on permeability compared to the reaction-limited geochemical regime, which uniformly alters the pore space. Transport-limited precipitation is characterised by the narrowing of pore throats and their successive closure, resulting in a comparably high permeability reduction. Transport-limited dissolution in turn leads to a predominant widening of existing flow paths, and thereby a significant increase in permeability. Reaction-limited processes uniformly alter the pore structure, whereby permeability is continuously changed. Nevertheless, permeabilities of the reaction-limited dissolution paths are consistently lower compared to the precipitation paths, due to a smaller fluid-to-mineral interface of the narrow throats filled with the previously precipitated minerals.

The results clearly illustrate that the permeability evolution is history-dependent, where both the geochemical process and the dominating reaction regime govern the characteristic microstructural alterations of the pore space, which cannot be simply reversed by an inversion of the geochemical processes itself. The presented modelling approach enables to simulate cycles of precipitation and dissolution on highly resolved three-dimensional pore space models without the requirement of implementing complex and computationally expensive reactive transport simulations. Both evolving trends in permeability and causal changes in pore morphology can be quantified, and thereby explain the impact of microstructural alterations on hydraulic rock properties. Hence, the proposed method is of essential importance for a wide range of natural and engineered subsurface applications, since it improves process understanding, and thereby the predictive capabilities of reservoir models.

Code and data availability

The data and software will be only made available on specific user requests.

Author contributions

MW, TK, and MK conceived and designed the simulations; MW performed the research; MW and TK analysed the data; MW, TK, and MK wrote the paper. All authors read and agreed to the published version of the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “European Geosciences Union General Assembly 2022, EGU Division Energy, Resources & Environment (ERE)”. It is a result of the EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022.

Acknowledgements

We appreciate the constructive comments from Elco Luijendijk and an anonymous reviewer.

Financial support

This publication has been supported by the funding programme “Open Access Publikationskosten” Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project Number 491075472.

The article processing charges for this open-access publication were covered by the Helmholtz Centre Potsdam – GFZ German Research Centre for Geosciences.

Review statement

This paper was edited by Christopher Juhlin and reviewed by Elco Luijendijk and one anonymous referee.

References

Andrä, H., Combaret, N., Dvorkin, J., Glatt, E., Han, J., Kabel, M., Keehm, Y., Krzikalla, F., Lee, M., Madonna, C., Marsh, M., Mukerji, T., Saenger, E., Sain, R., Saxena, N., Ricker, S., Wiegmann, A., and Zhan, X.: Digital rock physics benchmarks – Part II: Computing effective properties, Comput. Geosci., 50, 33–43, https://doi.org/10.1016/j.cageo.2012.09.008, 2013. a

Beckingham, L.: Evaluation of macroscopic porosity-permeability relationships in heterogeneous mineral dissolution and precipitation scenarios, Water Resour. Res., 53, 10217–10230, https://doi.org/10.1002/2017WR021306, 2017. a, b

Berg, C. F., Lopez, O., and Berland, H.: Industrial applications of digital rock technology, J. Petr. Sci. Eng., 157, 131–147, https://doi.org/10.1016/j.petrol.2017.06.074, 2017. a

Bernabé, Y., Mok, U., and Evans, B.: Permeability-porosity relationships in rocks subjected to various evolution processes, Thermo-Hydro-Mechanical Coupling in Fractured Rock, 160, 937–960, https://doi.org/10.1007/978-3-0348-8083-1_9, 2003. a

Blunt, M., Bijeljic, B., Dong, H., Gharbi, O., Iglauer, S., Mostaghimi, P., Paluszny, A., and Pentland, C.: Pore-scale imaging and modelling, Adv. Water Resour., 51, 197–216, https://doi.org/10.1016/j.advwatres.2012.03.003, 2013. a

Carman, P.: Fluid flow through granular beds, Transactions of the Institution of Chemical Engineers, 15, 150–166, 1937. a

Cil, M., Xie, M., Packman, A., and Buscarnera, G.: Solute mixing regulates heterogeneity of mineral precipitation in porous media, Geophys. Res. Lett., 44, 6658–6666, https://doi.org/10.1002/2017GL073999, 2017. a

Crandell, L., Peters, C., Um, W., Jones, K., and Lindquist, W.: Changes in the pore network structure of Hanford sediment after reaction with caustic tank wastes, J. Contam. Hydrol., 131, 89–99, https://doi.org/10.1016/j.jconhyd.2012.02.002, 2012. a

De Lucia, M., Lagneau, V., de Fouquet, C., and Bruno, R.: The influence of spatial variability on 2D reactive transport simulations, Comptes Rendus – Geoscience, 343, 406–416, https://doi.org/10.1016/j.crte.2011.04.003, 2011. a

Ettemeyer, F., Lechner, P., Hofmann, T., Andrä, H., Schneider, M., Grund, D., Volk, W., and Günther, D.: Digital sand core physics: Predicting physical properties of sand cores by simulations on digital microstructures, Int. J. Solids Struct., 188, 155–168, https://doi.org/10.1016/j.ijsolstr.2019.09.014, 2020. a

Fazeli, H., Masoudi, M., Patel, R., Aagaard, P., and Hellevang, H.: Pore-scale modeling of nucleation and growth in porous media, ACS Earth Space Chem., 4, 249–260, https://doi.org/10.1021/acsearthspacechem.9b00290, 2020. a

Gomez, C. T., Dvorkin, J., and Vanorio, T.: Laboratory measurements of resistivity and velocity for Fontainebleau sandstones, Society of Exploration Geophysicists International Exposition and 80th Annual Meeting 2010, SEG 2010, 75, 2452–2458, https://doi.org/10.1190/1.3513347, 2010. a

Gostick, J., Khan, Z., Tranter, T., Kok, M., Agnaou, M., Sadeghi, M., and Jervis, R.: PoreSpy: A python toolkit for quantitative analysis of porous media images, J. Open Source Softw., 4, 1296, https://doi.org/10.21105/joss.01296, 2019. a

Herring, A., Middleton, J., Walsh, R., Kingston, A., and Sheppard, A.: Flow rate impacts on capillary pressure and interface curvature of connected and disconnected fluid phases during multiphase flow in sandstone, Adv. Water Resour., 107, 460–469, https://doi.org/10.1016/j.advwatres.2017.05.011, 2017. a

Hofmann, H., Blöcher, G., Börsing, N., Maronde, N., Pastrik, N., and Zimmermann, G.: Potential for enhanced geothermal systems in low permeability limestones – stimulation strategies for the Western Malm karst (Bavaria), Geothermics, 51, 351–367, https://doi.org/10.1016/j.geothermics.2014.03.003, 2014. a

Hommel, J., Coltman, E., and Class, H.: Porosity-permeability relations for evolving pore space: A review with a focus on (bio-)geochemically altered porous media, Transport Porous M., 124, 589–629, https://doi.org/10.1007/s11242-018-1086-2, 2018. a

Jacquey, A., Cacace, M., Blöcher, G., Watanabe, N., and Scheck-Wenderoth, M.: Hydro-mechanical evolution of transport properties in porous media: Constraints for numerical simulations, Transport Porous M., 110, 409–428, https://doi.org/10.1007/s11242-015-0564-z, 2015. a

Keehm, Y., Mukerji, T., and Nur, A.: Computational rock physics at the pore scale: Transport properties and diagenesis in realistic pore geometries, Leading Edge, 20, 180–183, https://doi.org/10.1190/1.1438904, 2001. a

Kleinitz, W., Koehler, M., and Dietzsch, G.: The precipitation of salt in gas producing wells, in: SPE – European Formation Damage Conference, Society of Petroleum Engineers, 295–301, https://doi.org/10.2118/68953-MS, 2001. a

Kühn, M.: Reactive flow modeling of hydrothermal systems, vol. 103, Lecture Notes in Earth Sciences, Springer-Verlag Berlin Heidelberg, https://doi.org/10.1007/b13902, 2004. a

Lebedev, M., Zhang, Y., Sarmadivaleh, M., Barifcani, A., Al-Khdheeawi, E., and Iglauer, S.: Carbon geosequestration in limestone: Pore-scale dissolution and geomechanical weakening, Int. J. Greenh. Gas Contr., 66, 106–119, https://doi.org/10.1016/j.ijggc.2017.09.016, 2017. a

Lesueur, M., Poulet, T., and Veveakis, M.: Permeability hysteresis from microchannels opening during dissolution/reprecipitation cycle, Geophys. Res. Lett., 47, 1–8, https://doi.org/10.1029/2020GL088674, 2020. a, b, c

Luquot, L., Hebert, V., and Rodriguez, O.: Calculating structural and geometrical parameters by laboratory measurements and X-ray microtomography: a comparative study applied to a limestone sample before and after a dissolution experiment, Solid Earth, 7, 441–456, https://doi.org/10.5194/se-7-441-2016, 2016. a

Ma, J.: Review of permeability evolution model for fractured porous media, J. Rock Mech. Geotech. Eng., 7, 351–357, https://doi.org/10.1016/j.jrmge.2014.12.003, 2015. a

Menke, H., Bijeljic, B., Andrew, M., and Blunt, M. J.: Dynamic pore-scale imaging of reactive transport in heterogeneous carbonates at reservoir conditions, Energy Proced., 63, 5503–5511, https://doi.org/10.1016/j.egypro.2014.11.583, 2014. a

Miller, K., Vanorio, T., and Keehm, Y.: Evolution of permeability and microstructure of tight carbonates due to numerical simulation of calcite dissolution, J. Geophys. Res.-Sol. Ea., 122, 4460–4474, https://doi.org/10.1002/2017JB013972, 2017. a, b

Niu, Q. and Zhang, C.: Permeability prediction in rocks experiencing mineral precipitation and dissolution: A numerical study, Water Resour. Res., 55, 3107–3121, https://doi.org/10.1029/2018WR024174, 2019. a, b, c

Noiriel, C., Steefel, C., Yang, L., and Bernard, D.: Effects of pore-scale precipitation on permeability and flow, Adv. Water Resour., 95, 125–137, https://doi.org/10.1016/j.advwatres.2015.11.013, 2016. a

Regenspurg, S., Feldbusch, E., Byrne, J., Deon, F., Driba, D., Henninges, J., Kappler, A., Naumann, R., Reinsch, T., and Schubert, C.: Mineral precipitation during production of geothermal fluid from a Permian Rotliegend reservoir, Geothermics, 54, 122–135, https://doi.org/10.1016/j.geothermics.2015.01.003, 2015. a

Sain, R., Mukerji, T., and Mavko, G.: On microscale heterogeneity in granular media and its impact on elastic property estimation, Geophysics, 81, D561–D571, https://doi.org/10.1190/geo2016-0152.1, 2016. a

Schepers, A. and Milsch, H.: Dissolution-precipitation reactions in hydrothermal experiments with quartz-feldspar aggregates, Contr. Mineral. Petrol., 165, 83–101, https://doi.org/10.1007/s00410-012-0793-x, 2013. a

Schepp, L., Ahrens, B., Balcewicz, M., Duda, M., Nehler, M., Osorno, M., Uribe, D., Steeb, H., Nigon, B., Stöckhert, F., Swanson, D., Siegert, M., Gurris, M., and Saenger, E.: Digital rock physics and laboratory considerations on a high-porosity volcanic rock, Sci. Rep.-UK, 10, 1–16, https://doi.org/10.1038/s41598-020-62741-1, 2020. a

Shulakova, V., Sarout, J., Pimienta, L., Lebedev, M., Mayo, S., Clennell, M., and Pervukhina, M.: Effect of supercritical CO2 on carbonates: Savonnières sample case study, Geophys. Prospect., 65, 251–265, https://doi.org/10.1111/1365-2478.12369, 2017. a

Tenthorey, E. and Scholz, C.: Mapping secondary mineral formation in porous media using heavy metal tracers, J. Geophys. Res.-Sol. Ea., 107, 1–12, https://doi.org/10.1029/2000jb000109, 2002. a

Walderhaug, O., Eliassen, A., and Aase, N. E.: Prediction of permeability in quartz-rich sandstones: examples from the Norwegian continental shelf and the Fontainebleau sandstone, J. Sedimen. Res., 82, 899–912, https://doi.org/10.2110/jsr.2012.79, 2012.  a

Weller, H., Tabor, G., Jasak, H., and Fureby, C.: A tensorial approach to computational continuum mechanics using object-oriented techniques, Comput. Phys., 12, 620–631, https://doi.org/10.1063/1.168744, 1998. a

Wetzel, M.: Pore space alterations and their impact on hydraulic and mechanical rock properties quantified by numerical simulations, PhD thesis, Universität Potsdam, https://doi.org/10.25932/publishup-51206, 2021. a

Wetzel, M., Kempka, T., and Kühn, M.: Hydraulic and mechanical impacts of pore space alterations within a sandstone quantified by a flow velocity-dependent precipitation approach, Materials, 13, 3100, https://doi.org/10.3390/ma13143100, 2020a. a

Wetzel, M., Kempka, T., and Kühn, M.: Digital rock physics approach to simulate hydraulic effects of anhydrite cement in Bentheim sandstone, Adv. Geosci., 54, 33–39, https://doi.org/10.5194/adgeo-54-33-2020, 2020b. a

Wetzel, M., Kempka, T., and Kühn, M.: Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics, Minerals, 11, 151, https://doi.org/10.3390/min11020151, 2021. a

Wu, Y., Lin, C., Yan, W., Liu, Q., Zhao, P., and Ren, L.: Pore-scale simulations of electrical and elastic properties of shale samples based on multicomponent and multiscale digital rocks, Marine Petr. Geol., 117, 104369, https://doi.org/10.1016/j.marpetgeo.2020.104369, 2020. a

Download
Short summary
Porosity-permeability relations are simulated for a precipitation-dissolution cycle in a virtual sandstone. A hysteresis in permeability is observed depending on the geochemical process and dominating reaction regime, whereby permeability varies by more than two orders of magnitude. Controlling parameters for this hysteresis phenomenon are the closure and re-opening of micro-scale flow channels, derived from changes in pore throat diameter and connectivity of the pore network.