Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media
Gypsum dissolution leads to the development of karstic features within much shorter timescales than in other sedimentary rocks, potentially leading to rapid deterioration of groundwater quality and increasing the risk of catastrophes caused by subsidence. Here, we present a 2-D reactive transport mo...
Main Authors: | , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
MDPI AG
2022-01-01
|
Series: | Energies |
Subjects: | |
Online Access: | https://www.mdpi.com/1996-1073/15/3/761 |
_version_ | 1797488199006683136 |
---|---|
author | Jesús F. Águila Javier Samper Belén Buil Paloma Gómez Luis Montenegro |
author_facet | Jesús F. Águila Javier Samper Belén Buil Paloma Gómez Luis Montenegro |
author_sort | Jesús F. Águila |
collection | DOAJ |
description | Gypsum dissolution leads to the development of karstic features within much shorter timescales than in other sedimentary rocks, potentially leading to rapid deterioration of groundwater quality and increasing the risk of catastrophes caused by subsidence. Here, we present a 2-D reactive transport model to evaluate gypsum karstification in physically and chemically heterogeneous systems. The model considers a low-permeability rock matrix composed mainly of gypsum and a discontinuity (fracture), which acts as a preferential water pathway. Several scenarios are analyzed and simulated to investigate the relevance for gypsum karstification of: (1) the dynamic update of flow and transport parameters due to porosity changes; (2) the spatial distribution of minerals in the rock matrix; (3) the time evolution of water inflows through the boundaries of the model; (4) the functions relating permeability, <i>k</i>, to porosity, <i>ϕ</i>. The average porosity of the matrix after 1000 years of simulation increases from 0.045 to 0.29 when flow, transport, and chemical parameters and the water inflows through the boundary are dynamically updated according to the porosity changes. On the contrary, the porosity of the matrix hardly changes when the porosity feedback effect is not considered, while its average increases to 0.13 if the water inflow occurs through the discontinuity. Moreover, the dissolution of small amounts of highly soluble sulfate minerals plays a major role in the development of additional fractures. The increase in hydraulic conductivity is largest for the power law with an exponent of <i>n</i> = 5, as well as the Kozeny-Carman and the modified Fair-atch <i>k-ϕ</i> relationships. The gypsum dissolution front propagates into the matrix faster when the power law with <i>n</i> = 2 and 3 and the Verma–Pruess <i>k-ϕ</i> relationships are used. |
first_indexed | 2024-03-09T23:59:37Z |
format | Article |
id | doaj.art-b6180eadebd54b31a6479d6c7306156c |
institution | Directory Open Access Journal |
issn | 1996-1073 |
language | English |
last_indexed | 2024-03-09T23:59:37Z |
publishDate | 2022-01-01 |
publisher | MDPI AG |
record_format | Article |
series | Energies |
spelling | doaj.art-b6180eadebd54b31a6479d6c7306156c2023-11-23T16:19:27ZengMDPI AGEnergies1996-10732022-01-0115376110.3390/en15030761Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured MediaJesús F. Águila0Javier Samper1Belén Buil2Paloma Gómez3Luis Montenegro4Center for Advanced Scientific Research (CICA), Civil Engineering School, Campus de Elviña, University of A Coruña, 15071 A Coruña, SpainCenter for Advanced Scientific Research (CICA), Civil Engineering School, Campus de Elviña, University of A Coruña, 15071 A Coruña, SpainCenter for Energy, Environment and Technology Research (CIEMAT), 28040 Madrid, SpainCenter for Energy, Environment and Technology Research (CIEMAT), 28040 Madrid, SpainCenter for Advanced Scientific Research (CICA), Civil Engineering School, Campus de Elviña, University of A Coruña, 15071 A Coruña, SpainGypsum dissolution leads to the development of karstic features within much shorter timescales than in other sedimentary rocks, potentially leading to rapid deterioration of groundwater quality and increasing the risk of catastrophes caused by subsidence. Here, we present a 2-D reactive transport model to evaluate gypsum karstification in physically and chemically heterogeneous systems. The model considers a low-permeability rock matrix composed mainly of gypsum and a discontinuity (fracture), which acts as a preferential water pathway. Several scenarios are analyzed and simulated to investigate the relevance for gypsum karstification of: (1) the dynamic update of flow and transport parameters due to porosity changes; (2) the spatial distribution of minerals in the rock matrix; (3) the time evolution of water inflows through the boundaries of the model; (4) the functions relating permeability, <i>k</i>, to porosity, <i>ϕ</i>. The average porosity of the matrix after 1000 years of simulation increases from 0.045 to 0.29 when flow, transport, and chemical parameters and the water inflows through the boundary are dynamically updated according to the porosity changes. On the contrary, the porosity of the matrix hardly changes when the porosity feedback effect is not considered, while its average increases to 0.13 if the water inflow occurs through the discontinuity. Moreover, the dissolution of small amounts of highly soluble sulfate minerals plays a major role in the development of additional fractures. The increase in hydraulic conductivity is largest for the power law with an exponent of <i>n</i> = 5, as well as the Kozeny-Carman and the modified Fair-atch <i>k-ϕ</i> relationships. The gypsum dissolution front propagates into the matrix faster when the power law with <i>n</i> = 2 and 3 and the Verma–Pruess <i>k-ϕ</i> relationships are used.https://www.mdpi.com/1996-1073/15/3/761gypsum karstificationreactive transport modeldissolutionporosity changespermeability-porosity relationshipfractured media |
spellingShingle | Jesús F. Águila Javier Samper Belén Buil Paloma Gómez Luis Montenegro Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media Energies gypsum karstification reactive transport model dissolution porosity changes permeability-porosity relationship fractured media |
title | Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media |
title_full | Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media |
title_fullStr | Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media |
title_full_unstemmed | Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media |
title_short | Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media |
title_sort | reactive transport model of gypsum karstification in physically and chemically heterogeneous fractured media |
topic | gypsum karstification reactive transport model dissolution porosity changes permeability-porosity relationship fractured media |
url | https://www.mdpi.com/1996-1073/15/3/761 |
work_keys_str_mv | AT jesusfaguila reactivetransportmodelofgypsumkarstificationinphysicallyandchemicallyheterogeneousfracturedmedia AT javiersamper reactivetransportmodelofgypsumkarstificationinphysicallyandchemicallyheterogeneousfracturedmedia AT belenbuil reactivetransportmodelofgypsumkarstificationinphysicallyandchemicallyheterogeneousfracturedmedia AT palomagomez reactivetransportmodelofgypsumkarstificationinphysicallyandchemicallyheterogeneousfracturedmedia AT luismontenegro reactivetransportmodelofgypsumkarstificationinphysicallyandchemicallyheterogeneousfracturedmedia |