REVIEW 4 major objections 5 minor 30 references
Reactive Transport Simulation of Silicate-Rich Shale Rocks when Exposed to CO2 Saturated Brine Under High Pressure and High Temperature
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper claims that a 1D CrunchTope model of Permian shale in CO2-saturated brine reproduces measured mineral-reaction depths (1080 and 1480 μm vs.
desk verdict A transparent but overclaimed CrunchTope validation study: the model reproduces reaction depths and qualitative trends, but not the two dominant mechanisms (quartz stress corrosion, clay swelling) that the experiments flag. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is the coupled reaction-transport loop implemented in CrunchTope, an open-source reactive transport code. Aqueous species diffuse through a 1D brine-rock grid under a porosity-dependent effective diffusion coefficient $D_i^e = \phi^m D_i$, while mineral dissolution and precipitation obey the transition-state-theory rate law $R = A k a_{H^+}^n (1 - IAP/K_{sp})$. Nucleation of secondary phases (amorphous SiO2, kaolinite, gibbsite) is added through the classical nucleation barrier $\Delta G^* = 16\pi v^2 a^3 / [3 k_B^2 T^2 (\ln(IAP/K_{sp}))^2]$, and porosity is updated from the evolving mineral volume fractions. This machinery converts the batch-reaction chemistry into spatial phase-concentration profiles that can be compared directly with the SEM-EDS measurements.
What would settle it
At 14 days, the paper reports the simulated quartz-rich phase concentration at the brine-rock interface as 48.30%, against about 13.90% measured; a decisive check would be to add a stress-corrosion-cracking rate law for quartz and see whether the interface concentration drops toward the measured value while still matching the reaction depths and porosity profile. If the quartz concentration remains near its initial 54.9% while the measured value is near 13.9%, the claim that the simulation replicates chemo-mechanical alteration is falsified.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that a continuum-scale, one-dimensional reactive transport simulation can replicate the key pore-scale chemo-mechanical response of silicate-rich shale to CO2-saturated brine: feldspar dissolves, quartz dissolves and then partly reprecipitates, clay precipitates near the reacted surface, and the reaction fronts advance at 1080 μm (14 days) and 1480 μm (28 days), matching the experimentally measured 1100 μm and 1500 μm. The authors conclude that CrunchTope successfully simulated the pore-scale transport processes and chemo-mechanical reactions in the shale, including the nucleation and growth of secondary minerals such as amorphous SiO2 and kaolinite, and that the simulation results are validated by the experimental findings. The paper is explicit that two experimentally observed processes—stress-corrosion cracking of quartz and clay swelling—are not captured by the code, which is why the simulated magnitudes of quartz dissolution and clay precipitation fall short of the measured values.
Load-bearing premise
The load-bearing premise is that quartz is essentially inert and highly weathering-resistant, which the experiments contradict by showing large quartz loss through stress-corrosion cracking, and if that premise fails the simulated mineral balance and chemo-mechanical predictions understate reality; the model also cannot simulate clay swelling, further limiting clay-phase magnitudes.
Editorial extensions
If this is right
- Reaction depth grows with the square root of exposure time, from 1080 μm at 14 days to 1480 μm at 28 days, matching the measured 1100 μm and 1500 μm; this means shorter simulations can be extrapolated to longer exposures by Fickian scaling.
- Porosity near the brine-rock interface rises from 5% to about 15–18%, with only about 0.89–1.19% mineral precipitation, so CO2-saturated brine does not seal shale pores by precipitation and further brine diffusion remains possible.
- The model nucleates and grows secondary kaolinite and amorphous SiO2, reproducing the experimentally inferred pathway in which feldspar transforms to clay and quartz precipitates.
- Because the simulated reaction fronts align with measured depths, 1D reactive transport runs can estimate the thickness of the altered zone in shale caprocks without a full experimental campaign.
Reading between the lines
- A natural next test is to replace the inert-quartz assumption with a stress-corrosion-cracking rate law; if calibrated to the measured interface concentration, the model could predict the experimentally observed >50% drop in quartz-rich indentation modulus and thus caprock mechanical integrity, not just mineralogy.
- The same 1D setup could be extended to 2D or 3D with explicit fractures to ask whether the thin precipitation bands (0.89–1.19%) change fracture permeability differently than matrix porosity.
- Because the authors validate only phase concentrations and reaction depths, a stronger test of the model would compare simulated effluent chemistry (Na+, SiO2(aq), pH evolution) and nanoindentation maps against the experimental record; the paper reports only pH at the interface, not the full solute record.
- If the porosity feedback is right, the main risk to storage integrity in silicate-rich shales may be stress-corrosion weakening rather than pore clogging, which would shift monitoring and caprock-failure assessments toward mechanical damage indicators.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a one-dimensional reactive transport simulation (CrunchTope) of a Permian shale sample exposed to CO2-saturated brine at 100 °C and 12.4 MPa for 14 and 28 days, with the stated aim of replicating the experimental chemo-mechanical alteration profiles reported by Prakash et al. (2024). The authors report feldspar dissolution, dissolution followed by precipitation of quartz-rich phases, clay precipitation, and reaction depths of 1080 μm and 1480 μm at 14 and 28 days, respectively, compared with experimental depths of 1100 μm and 1500 μm. They conclude that CrunchTope successfully simulated the pore-scale transport processes and chemo-mechanical reactions, while acknowledging in Section 4 two key limitations: the model cannot represent quartz dissolution by stress corrosion cracking or clay swelling.
Significance. If convincingly validated, the simulation would offer a faster and cheaper complement to batch experiments for screening shale-CO2 interactions. The paper has notable strengths: it states the governing equations and parameter values clearly, it compares against a specific experimental dataset, and it explicitly enumerates model limitations. However, as presented, the validation is largely qualitative: the two chemo-mechanical mechanisms identified by the experimental study as dominant—quartz stress-corrosion cracking and clay swelling—are absent from the model, and the quantitative phase-concentration agreement is poor. The reaction-depth agreement is consistent with a diffusion-controlled front under generic reactive chemistry and is not sufficient to establish predictive skill. The manuscript therefore needs substantial revision to temper its central claim and to provide a more rigorous validation.
major comments (4)
- [Section 4; Section 3.1] The concluding claim that "CrunchTope successfully simulated the pore-scale transport processes and chemo-mechanical reactions" is not supported by the quantitative comparison in Section 3.1. The simulated quartz-rich phase concentration falls only from 54.9% at 30 μm to 48.3% at the brine-rock interface, while the experimental profile falls from about 44.6% at 132 μm to about 13.9% at the interface. Section 3.1 also states that simulated clay precipitation is lower than measured because CrunchTope cannot represent clay swelling. Since quartz dissolution and clay swelling are the two dominant chemo-mechanical mechanisms identified in the experiments, the conclusion should be revised to claim reproduction of the qualitative reaction sequence and penetration depth, not successful simulation of chemo-mechanical reactions.
- [Section 2.2] The model is configured to reproduce the experimental outcomes in a circular manner. The authors write that "for actual replication of the experimental study, the nucleation of amorphous silicon dioxide, kaolinite, and gibbsite was considered," and they add a thin inert substrate layer in the brine domain to trigger this nucleation. The simulated precipitation of these phases is therefore an input assumption rather than an independent model prediction. Please separate calibrated inputs from predicted outputs and provide a sensitivity analysis (e.g., varying J0, interfacial energy a, and the presence of the substrate) to demonstrate that the main conclusions do not depend on these ad hoc choices.
- [Section 3.1; Eq. (3)] The reaction-depth agreement (1080 vs 1100 μm at 14 days, and 1480 vs 1500 μm at 28 days) is presented as validation, but a diffusion-controlled reaction front advances as the square root of time under almost any reactive chemistry, so this agreement is weak evidence. The front location is controlled by the effective diffusion coefficient D_i^e = φ^m D_i (Eq. 3) with m = 2.40 and D = 1.48×10^−9 m²/s, values that were selected without independent constraint, and the inert nucleation substrate introduced in Section 2.2 also influences the front. Please report parameter sensitivity or uncertainty bounds (e.g., variations in m, D, J0, and interfacial energy) to show that the depth match is not merely a tuning outcome.
- [Section 3.4] The porosity profile in Section 3.4 is presented without any comparison to experimental porosity measurements. The statement that porosity increases from 5.0% to 15.3% at 14 days and to 17.9% at 28 days is a model output, not a validated result. Since porosity evolution is a central component of the claimed "transport processes," either provide the corresponding experimental data or present this as a qualitative model prediction rather than a validated outcome.
minor comments (5)
- [Table 1] The mineral proportions in Table 1 sum to 105.14% (54.59 + 23.31 + 16.61 + 2.54 + 1.67 + 1.27 + 5.15), which is internally inconsistent; please clarify whether these are weight fractions, volume fractions, or unnormalized values and correct the normalization.
- [Section 2.1, Eq. (1)] The notation in Eq. (1) is confusing because the summation index k is used both for the denominator and for the running index k = 1, ..., n, while the left-hand side uses a different index i; please rewrite the equation with distinct indices and define φ, w_i, and ρ_i explicitly.
- [Abstract; Section 2.2] The software name is spelled inconsistently as "Crunch Tope" in the Abstract and "CrunchTope" in the body; please use a single consistent spelling throughout the manuscript.
- [References] Reference 27 (Varanasi, n.d.) is incomplete; please provide the full dissertation title, year, and institution, or remove the citation if it is not essential.
- [Figures] The manuscript refers to Figures 1–6, but the figures are not embedded in the text of the submission; please ensure that all figures are included with sufficient resolution so that the experimental data points can be distinguished from the simulated curves.
Circularity Check
Partial circularity: precipitation of SiO2(am), kaolinite, and gibbsite is inserted as model input to replicate the experiments, then reported as a successful simulation.
-
fitted input called prediction
[Section 2.2 (Reactive Transport Simulation Modeling); see also Section 4 (Conclusion)]
"The model was also updated to include the nucleation of the minerals, such as amorphous silicon dioxide (SiO2(am)), kaolinite, and gibbsite. To simulate the nucleation of these minerals, a thin inert substrate layer was included in the brine domain closer to the brine rock interface."
The Introduction states this was done "for actual replication of the experimental study" (Section 1), so the phases that the experiments had already shown precipitating are inserted into the model as nucleating phases with an artificial inert substrate.
full rationale
This is a validation-style study rather than a first-principles derivation, and the self-citation to Prakash et al. (2024) is not itself circular in the harmful sense: the cited nanoindentation/SEM-EDS data are external, falsifiable experimental measurements, not outputs of the model. The main circularity is narrower and located at the nucleation setup. The authors add nucleation of SiO2(am), kaolinite, and gibbsite "for actual replication of the experimental study" (Section 1) and install a thin inert substrate layer to make those phases nucleate near the interface (Section 2.2), then treat the resulting precipitation as evidence that CrunchTope "successfully simulated" the chemo-mechanical reactions (Section 4). That specific validated outcome is therefore an input assumption rather than a prediction. The paper's own limitations are correctness concerns, not additional circularity: Section 3.1 admits "CrunchTope simulated results do not seem to demonstrate significant dissolution" for quartz and "CrunchTope is its inability to simulate the swelling of the clay phase," which undercut the strength of the validation but do not make the derivation self-referential. The reaction-depth agreement is also weak evidence because a diffusion front scales with the square root of time under Fickian transport, but this is not counted as circular since the paper does not state that the diffusion coefficient or cementation exponent were fitted to those depths. Overall, one central "prediction" reduces to the model input, giving partial circularity (score 6) rather than a fully forced derivation.
Assumptions & free parameters
free parameters (5)
- Cementation exponent m =
2.40
- Effective diffusion coefficient D =
1.48e-9 m2/s
- Nucleation kinetic factor J0 =
1.00e-8 mol/m2/s
- Interfacial energy a =
47 ± 1 mJ/m2
- Inclusion of SiO2(am), kaolinite, and gibbsite nucleation =
n/a
assumptions (5)
- domain assumption Quartz is a highly stable and insoluble mineral with high resistance to weathering
- domain assumption Clay swelling and grain detachment need not be simulated
- standard math Transition state theory (TST) rate law describes mineral dissolution and precipitation kinetics
- domain assumption A single uniform diffusion coefficient represents all aqueous species
- domain assumption A 1D domain captures the pore-scale chemo-mechanical response of the heterogeneous rock
invented entities (1)
-
Thin inert substrate layer for mineral nucleation
Cite this review
Pith. "Pith review of Reactive Transport Simulation of Silicate-Rich Shale Rocks when Exposed to CO2 Saturated Brine Under High Pressure and High Temperature." pith.science (2026). https://pith.science/paper/ZQ525EZM
@misc{pith2026250605122,
author = {Pith},
title = {Pith review of: Reactive Transport Simulation of Silicate-Rich Shale Rocks when Exposed to CO2 Saturated Brine Under High Pressure and High Temperature},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZQ525EZM}},
note = {Machine review of arXiv:2506.05122}
}
read the original abstract
This study examines the feasibility of carbon dioxide storage in shale rocks and the reliability of reactive transport models in achieving accurate replication of the chemo-mechanical interactions and transport processes transpiring in these rocks when subjected to CO2 saturated brine. Owing to the heterogeneity of rocks, experimental testing for adequate deductions and findings, could be an expensive and time-intensive process. Therefore, this study proposes utilization of reactive transport modeling to replicate the pore-scale chemo-mechanical reactions and transport processes occurring in silicate-rich shale rocks in the presence of CO2 saturated brine under high pressure and high temperature. For this study, Crunch Tope has been adopted to simulate a one-dimensional reactive transport model of a Permian rock specimen exposed to the acidic brine at a temperature of 100 {\deg}C and pressure of 12.40 MPa (1800 psi) for a period of 14 and 28 days. The results demonstrated significant dissolution followed by precipitation of quartz rich phases, precipitation and swelling of clay rich phases, and dissolution of feldspar rich phases closer to the acidic brine-rock interface. Moreover, porosity against reaction depth curve showed nearly 1.00% mineral precipitation occur at 14 and 28 days, which is insufficient to completely fill the pore spaces.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Abedi, S., Slim, M., Hofmann, R., Bryndzia, T., & Ulm, F. J. (2016). Nanochemo -mechanical signature of organic -rich shales: A coupled indentation -EDX analysis. Acta Geotechnica, 11(3), 559–572. https://doi.org/10.1007/s11440-015-0426-4
-
[2]
Abedi, S., Slim, M., & Ulm, F. J. (2016). Nanomechanics of organic-rich shales: the role of thermal maturity and organic matter content on texture. Acta Geotechnica , 11(4), 775 –787. https://doi.org/10.1007/s11440-016-0476-2
-
[3]
S., Raza, A., Zhang, Y., Mahmoud, M., & Arif, M
Adila, A. S., Raza, A., Zhang, Y., Mahmoud, M., & Arif, M. (2023). Geochemical Interactions Among Rock/CO2/Brine Systems: Implications for CO2 Geo -Storage. Society of Petroleum 13 Engineers - Gas and Oil Technology Showcase and Conference, GOTS 2023 , 1 –11. https://doi.org/10.2118/214029-MS
-
[4]
Aman, M., Espinoza, D. N., Ilgen, A. G., Major, J. R., Eichhubl, P., & Dewers, T. A. (2018). CO2- induced chemo-mechanical alteration in reservoir rocks assessed via batch reaction experiments and scratch testing. Greenhouse Gases: Science and Technology , 8(1), 133 –149. https://doi.org/10.1002/ghg.1726
-
[5]
Asadi, P., Salek, M. F., & Beckingham, L. E. (2024). Modeling the Spatial and Temporal Evolution of Fractured Systems with a Heterogenous Mineralogy. ACS Earth and Space Chemistry , 8(6), 1133–1145. https://doi.org/10.1021/acsearthspacechem.3c00322
-
[6]
Berrezueta, E., Moita, P., Pedro, J., Abdoulghafour, H., Mirão, J., Beltrame, M., Barrulas, P., Araújo, A., Caeiro, M. H., Luís, L., & Carneiro, J. (2023). Laboratory experiments and modelling of the geochemical interaction of a gabbro -anorthosite with seawater and supercritical CO2: A mineral carbonation study. Geoenergy Science and Engineering , 228(Ju...
-
[7]
Clark, A. C., & Vanorio, T. (2016). The rock physics and geochemistry of carbonates exposed to reactive brines. Journal of Geophysical Research: Solid Earth , 121, 1497 –1513. https://doi.org/10.1002/2014JB011037.Received
-
[8]
Dávila, G., Dalton, L., Crandall, D. M., Garing, C., Werth, C. J., & Druhan, J. L. (2020). Reactive alteration of a Mt. Simon Sandstone due to CO2 -rich brine displacement. Geochimica et Cosmochimica Acta, 271, 227–247. https://doi.org/10.1016/j.gca.2019.12.015
Show all 30 references
-
[9]
Ben, Bennour, Z., Gholami, R., & Hossain, M
Fatah, A., Mahmud, H. Ben, Bennour, Z., Gholami, R., & Hossain, M. (2022). Geochemical modelling of CO2 interactions with shale: Kinetics of mineral dissolution and precipitation on geological time scales. Chemical Geology , 592, 120742. https://doi.org/10.1016/j.chemgeo.2022.120742
2022
-
[10]
G., L iu, Q., & Etsell, T
Geramian, M., Osacky, M., Ivey, D. G., L iu, Q., & Etsell, T. H. (2016). Effect of Swelling Clay Minerals (Montmorillonite and Illite-Smectite) on Nonaqueous Bitumen Extraction from Alberta Oil Sands. Energy and Fuels, 30(10), 8083–8090. https://doi.org/10.1021/acs.energyfuels.6b01026
2016 doi
-
[11]
Houhou, R., & Laloui, L. (2022). Geomechanics for energy and the environment: Current developments. Geomechanics for Energy and the Environment , 32, 100345. https://doi.org/10.3390/en16186588
2022 doi
-
[12]
N., Siddiqui, S., & Thakur, G
Khan, M. N., Siddiqui, S., & Thakur, G. C. (2024). Recent Advances in Geoch emical and Mineralogical Studies on CO2 –Brine–Rock Interaction for CO2 Sequestration: Laboratory and Simulation Studies. Energies, 17(13), 3346. https://doi.org/10.3390/en17133346
2024 doi
-
[13]
I., & Jun, Y
Li, Q., Steefel, C. I., & Jun, Y. S. (2017). Incorporating Nanoscale Effects into a Continuum-Scale Reactive Transport Model for CO2-Deteriorated Cement. Environmental Science and Technology, 51(18), 10861–10871. https://doi.org/10.1021/acs.est.7b00594
2017 doi
-
[14]
W., Soong, Y., Hellevang, H., & Zhu, C
Liu, F., Lu, P., Griffith, C., Hedges, S. W., Soong, Y., Hellevang, H., & Zhu, C. (2012). CO2-brine- caprock interaction: Reactivity experiments on Eau Claire shale and a review of relevant literature. International Journal of Greenhouse Gas Control , 7, 153 –167. https://doi....
2012 doi
-
[15]
Martogi, D., & Abedi, S. (2020). Microscale approximation of the elastic mechanical properties of randomly oriented rock cuttings. Acta Geotechnica , 15(12), 3511 –3524. https://doi.org/10.1007/s11440-020-01020-9
2020 doi
-
[16]
Murugesu, M. P. (2024). Imaging Reactive Tran sport in Multiscale Geologcal Media . Stanford University
2024
-
[17]
Nguene, P. C. K. (2019). Chemo-Mechanical Characterization of Co2 -Rock-Fluid Interactions- Eagle Ford and Wolfcamp Shales. Texas A&M University
2019
-
[18]
Prakash, R., & Abedi, S. (2023). Experimentally Informed Simulation of Creep Behavior in Shale Rocks Induced by Chemo -mechanical Loading. Rock Mechanics and Rock Engineering , 56(9), 6631–6645. https://doi.org/10.1007/s00603-023-03413-0
2023 doi
-
[19]
D., Noshadravan, A., & Abedi, S
Prakash, R., Kana Ng uene, P., Seers, T. D., Noshadravan, A., & Abedi, S. (2019). Chemo - 14 mechanical investigation of CO2 -fluid-rock interaction in CO2 storage and CO2 -EOR processes in unconventional reservoirs. American Rock Mechanics Association (ARMA), 1–10
2019
-
[20]
A., & Abedi, S
Prakash, R., Mahgoub, S. A., & Abedi, S. (2024). Chemo -mechanical Alteration of Silicate-Rich Shale Rock after Exposure to CO2-Rich Brine at High Temperature and Pressure. Rock Mechanics and Rock Engineering, 57(8), 5317–5333. https://doi.org/10.1007/s00603-023-03664-x
2024 doi
-
[21]
Prakash, R., Nguene, P. C. K., Noshadravan, A., & Abedi, S. (2022). Chemical reactions of carbonate-rich mudstones with aqueous CO2 and their impacts on rock’s local microstructural and chemo-mechanical properties. Journal of Natural Gas Science and Enginee ring, 103, 104587. ...
2022
-
[22]
Sharma, P., Prakash, R., & Abedi, S. (2019). Effect of temperature on nano- and microscale creep properties of organic-rich shales. Journal of Petroleum Science and Engineering , 175(December 2018), 375–388. https://doi.org/10.1016/j.petrol.2018.12.039
2019 doi
-
[23]
Sobhbidari, F., & Hu, Q. (2021). Recent advances in the mechanical characterization of shales at nano-to micro -scales: A review. Mechanics of Materials , 162, 104043. https://doi.org/10.1016/j.mechmat.2021.104043
2021
-
[24]
I., & Lichtner, P
Steefel, C. I., & Lichtner, P. C. (1998). Multicomponent reactive transport in discrete fractures: I. Controls on reaction front geometry. Journal of Hydrology , 209(1–4), 186 –199. https://doi.org/10.1016/S0022-1694(98)00146-2
1998 doi
-
[25]
Tournassat, C., & Steefel, C. I. (2019). Reactive Transport Modeling of Coupled Processes in Nanoporous Media. Reviews in Mineralogy and Geochemistry , 85, 75 –109. https://doi.org/10.2138/rmg.2019.85.4
2019 doi
-
[26]
Ulmer-Scholle, D., Scholle, P., Schieber, J., & Raine, R. (2014). Diagenesis: iron sulfide, oxide & hydroxide cements. 347–360
2014
-
[27]
Varanasi, V. R. S. B. (n.d.). Effect of Pressure on CO2 -Rock-Fluid Interactions in Shale: Experimental and Reactive Tansport Modeling Study. Texas A&M University, College Station
-
[28]
Vialle, S., & Vanorio, T. (2011). Laboratory measurements of elastic properties of carbonate rocks during injection of reactive CO2 -saturated water. Geophysical Research Letters, 38(1), L01302. https://doi.org/10.1029/2010GL045606
2011 doi
-
[29]
Wu, Y., Li, Y., Luo, S., Lu, M., Zhou, N., Wang, D., & Zhang, G. (2020). Multiscale elastic anisotropy of a shale characterized by cross-scale big data nanoindentation. International Journal of Rock Mechanics and Mining Sciences , 134, 104458. https://doi.org/10.1016/j.ijrmms....
2020
-
[30]
Yuan, G., Cao, Y., Schulz, H.-M., Hao, F., Gluyas, J., Liu, K., Yang, T., Wang, Y., Xi, K., & Li, F. (2019). A review of feldspar alteration and its geological significance in sedimentary basins: From shallow aquifers to deep hydrocarbon reservoirs. Earth-Science Reviews, 191(...
2019 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.