Numerical simulations of the impacts of CO₂-induced calcite dissolution on the permeability and elastic parameters of rocks

Authors

YANG Bo, Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China; Jilin Provincial Key Laboratory of Water Resources and Environment, Jilin University, Changchun 130021, ChinaFollow
XU Tianfu, Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China; Jilin Provincial Key Laboratory of Water Resources and Environment, Jilin University, Changchun 130021, China
FENG Guanhong, Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China; Jilin Provincial Key Laboratory of Water Resources and Environment, Jilin University, Changchun 130021, China
ZHU Huixing, Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China; Jilin Provincial Key Laboratory of Water Resources and Environment, Jilin University, Changchun 130021, ChinaFollow

Abstract

Objective CO₂ dissolution in water tends to induce mineral reactions, thus significantly changing the pore structure of rocks and further affecting their seepage and mechanical properties. This phenomenon poses a major impact on the suitability and long-term safety of geologic CO₂ sequestration engineering. Methods Using a simulation technology that combines the lattice Boltzmann method (LBM) and the finite element method (FEM), this study systematically investigated the calcite dissolution characteristics of calcite-bearing rock samples under varying injection rates of a saturated solution of CO₂. Furthermore, the dynamic evolution patterns of the permeability and elastic parameters (bulk and shear moduli) of the rock samples were revealed. Results and Conclusions At lower injection rates, the chemical reaction products became enriched in the mid-to-distal part of the samples, inhibiting calcite dissolution in this part. Consequently, the chemical reactions occurred merely near the injection port. As the injection rate increased, the saturated solution of CO₂ exhibited an enhanced capacity to penetrate samples and then diluted reactants. As a result, calcites in the mid-to-distal part could participate in the reactions, resulting in a more uniform spatial distribution of calcite dissolution positions across the samples. Calcite dissolution significantly enhanced the rock permeability, especially at higher injection rates. At lower injection rates (30‒150 m/a), the solute transport was primarily achieved by diffusion, leading to an insignificant increase in rock permeability. When the injection rate increased to 750‒18 750 m/a, the solute transport mechanism gradually transitioned into advection, leading to greater differences in the rock matrix structure and more pronounced permeability enhancement. For the fitting of the permeability-porosity relationship, the power function model generally exhibited fit index n ranging from 2.8 to 6.5, while the Carman-Kozeny model showed n values varying from 0.9 to 4.6. Notably, the n values of both models showed a monotonically increasing trend with an increase in the injection rate. Comparative analysis reveals that the Carman-Kozeny model outperformed the power function model in terms of prediction accuracy. Calcite dissolution significantly reduced the elastic properties of rocks. As the porosity increased from 0.44 to 0.56, both the shear and bulk moduli decreased by approximately 20%. Additionally, both moduli declined rapidly initially and then decreased slowly in the late stage. Such nonlinear evolution renders the power function model more advantageous than the traditional linear model. This study also revealed that the elastic properties of rocks experienced exacerbated degradation with an increase in the injection rate. Compared to the bulk modulus, the shear modulus was more sensitive to variations in the injection rate. This study determines the key mechanism where the fluid injection rate dominates the evolution of the permeability and mechanical properties of rocks by governing the spatial distribution of dissolution, providing a theoretical basis for assessing the long-term safety of geologic CO₂ sequestration sites.

Keywords

saline aquifer, geologic CO2 sequestration, LBM-FEM coupling simulation, calcite dissolution, permeability, evolutionary pattern of elastic parameters, evolutionary laws

DOI

10.12363/issn.1001-1986.25.07.0548

Recommended Citation

YANG Bo, XU Tianfu, FENG Guanhong, et al. (2026) "Numerical simulations of the impacts of CO₂-induced calcite dissolution on the permeability and elastic parameters of rocks," Coal Geology & Exploration: Vol. 54: Iss. 1, Article 18.
DOI: 10.12363/issn.1001-1986.25.07.0548
Available at: https://cge.researchcommons.org/journal/vol54/iss1/18

Reference

[1] 邓晨光，孙小龙，刘可禹，等. 储层孔隙结构对CO₂地质封存井周围盐沉淀作用的影响机制[J]. 煤田地质与勘探，2025，53(10)：128−140.

DENG Chenguang，SUN Xiaolong，LIU Keyu，et al. Mechanisms underlying the impacts of reservoir pore structures on near–well salt precipitation in CO₂ geological storage[J]. Coal Geology & Exploration，2025，53(10)：128−140.

[2] KONG Xiangguo，LIN Xi，CAI Yuchu，et al. Effects of pore characteristics on CO₂ adsorption performance of coal slime with different metamorphic degrees[J]. Journal of Environmental Management，2024，368：122094.

[3] KONG Xiangguo，HU Jie，CAI Yuchu，et al. Pore structure and mineral composition characteristics of coal slime before and after ashing and the effects on CO₂ adsorption[J]. Environmental Science and Pollution Research，2024，31(43)：55597−55609.

[4] BASHIR A，ALI M，PATIL S，et al. Comprehensive review of CO₂ geological storage：Exploring principles，mechanisms，and prospects[J]. Earth–Science Reviews，2024，249：104672.

[5] 卢帅屹. 咸水层CO₂地质封存泥岩盖层力学稳定性研究：以吉林油田大情字井为例[D]. 长春：吉林大学，2024.

LU Shuaiyi. Study on mechanical stability of CO₂ geological storage mudstone cap in saltwater reservoir：A case study of Daqing Zi Well in Jilin Oilfield[D]. Changchun：Jilin University，2024.

[6] KHORAMIAN R，SALAUDEEN I，POURAFSHARY P，et al. Impact of relative permeability hysteresis on CO₂ storage in saline aquifers[J]. Greenhouse Gases：Science and Technology，2025，15(1)：3−12.

[7] XU Tianfu，APPS J A，PRUESS K. Reactive geochemical transport simulation to study mineral trapping for CO₂ disposal in deep arenaceous formations[J]. Journal of Geophysical Research：Solid Earth，2003，108(B2)：2071.

[8] WANG Zhiqiang，LIU Shuyang，LI Hangyu，et al. CO₂ storage in saline aquifers：A simulation on quantifying the impact of permeability heterogeneity[J]. Journal of Cleaner Production，2024，471：143415.

[9] DU Jiapei，ZHOU Annan，LIN Xiaoshan，et al. Modeling microstructural mechanical behavior of expansive soil at various water contents and dry densities by molecular dynamics simulation[J]. Computers and Geotechnics，2023，158：105371.

[10] SMITH M M，HAO Y，CARROLL S A. Development and calibration of a reactive transport model for carbonate reservoir porosity and permeability changes based on CO₂ core–flood experiments[J]. International Journal of Greenhouse Gas Control，2017，57：73−88.

[11] YANG Yongfei，LI Yingwen，YAO Jun，et al. Dynamic pore–scale dissolution by CO₂–saturated brine in carbonates：Impact of homogeneous versus fractured versus vuggy pore structure[J]. Water Resources Research，2020，56(4)：e2019WR026112.

[12] ZHANG Hao，ZHONG Ying，ZHANG Jiang，et al. Experimental research on deterioration of mechanical properties of carbonate rocks under acidified conditions[J]. Journal of Petroleum Science and Engineering，2020，185：106612.

[13] YANG Bing，WANG Haizhu，WANG Bin，et al. Effect of supercritical CO₂–water/brine–rock interaction on microstructures and mechanical properties of tight sandstone[J]. Transport in Porous Media，2023，149(1)：87−115.

[14] LI Hao，ZHONG Zuliang，LIU Xinrong，et al. Micro–damage evolution and macro–mechanical property degradation of limestone due to chemical effects[J]. International Journal of Rock Mechanics and Mining Sciences，2018，110：257−265.

[15] 李华焜，郑刘根，陈永春，等. 基于CT扫描的重构土壤孔隙结构及其对水盐运移影响[J]. 煤田地质与勘探，2024，52(4)：120−127.

LI Huakun，ZHENG Liugen，CHEN Yongchun，et al. Exploring the pore structure of reconstructed soils and its effects on water and salt transport based on CT scanning[J]. Coal Geology & Exploration，2024，52(4)：120−127.

[16] 张永超，刘昌岭，吴能友，等. 含水合物沉积物孔隙结构特征与微观渗流模拟研究[J]. 海洋地质前沿，2020，36(9)：23−33.

ZHANG Yongchao，LIU Changling，WU Nengyou，et al. Advances in the pore–structure characteristics and micro–scopic seepage numerical simulation of the hydrate–bearing sediments[J]. Marine Geology Frontiers，2020，36(9)：23−33.

[17] 林伟. 致密油储层数字岩心建模及自发渗吸研究[D]. 北京：中国科学院大学，2020.

LIN Wei. Study on the construction of digital rocks and spontaneous imbibition in tight oil reservoirs[D]. Beijing：University of Chinese Academy of Sciences，2020.

[18] YANG Bo，ZHU Huixing，TIAN Hailong，et al. Pore–scale direct simulation of CO₂，water，and oil three–phase relative permeabilities considering mineral wettability and dissolution[J]. Physics of Fluids，2025，37(2)：026630.

[19] LI Xiaoliang，LIU Zhentang，FENG Xiaojun，et al. Effects of acid sulfate and chloride ion on the pore structure and mechanical properties of sandstone under dynamic loading[J]. Rock Mechanics and Rock Engineering，2021，54(12)：6105−6121.

[20] CHEN Li，HE An，ZHAO Jianlin，et al. Pore–scale modeling of complex transport phenomena in porous media[J]. Progress in Energy and Combustion Science，2022，88：100968.

[21] 李衡. CO₂地质封存过程中孔隙演化的多尺度反应运移模型研究及应用[D]. 长春：吉林大学，2024.

LI Heng. Research and application of multi–scale reactive transport model for pore evolution during CO₂ geological storage[D]. Changchun：Jilin University，2024.

[22] AGRAWAL P，MASCINI A，BULTREYS T，et al. The impact of pore–throat shape evolution during dissolution on carbonate rock permeability：Pore network modeling and experiments[J]. Advances in Water Resources，2021，155：103991.

[23] CHEN Li，WANG Mengyi，KANG Qinjun，et al. Pore scale study of multiphase multicomponent reactive transport during CO₂ dissolution trapping[J]. Advances in Water Resources，2018，116：208−218.

[24] CHEN Yu，VALOCCHI A J，KANG Qinjun，et al. Inertial effects during the process of supercritical CO₂ displacing brine in a sandstone：Lattice Boltzmann simulations based on the continuum–surface–force and geometrical wetting models[J]. Water Resources Research，2019，55(12)：11144−11165.

[25] KANG Qinjun，CHEN Li，VALOCCHI A J，et al. Pore–scale study of dissolution–induced changes in permeability and porosity of porous media[J]. Journal of Hydrology，2014，517：1049−1055.

[26] AN Senyou，ERFANI H，HELLEVANG H，et al. Lattice–Boltzmann simulation of dissolution of carbonate rock during CO₂–saturated brine injection[J]. Chemical Engineering Journal，2021，408：127235.

[27] XIE Qiuheng，WANG Wendong，SU Yuliang，et al. Pore–scale study of calcite dissolution during CO₂–saturated brine injection for sequestration in carbonate aquifers[J]. Gas Science and Engineering，2023，114：204978.

[28] GARBOCZI E J. Finite element and finite difference programs for computing the linear electric and elastic properties of digital images of random materials[M]. Gaithersburg：Department of Commerce，1998.

[29] MEILLE S，GARBOCZI E J. Linear elastic properties of 2D and 3D models of porous materials made from elongated objects[J]. Modelling and Simulation in Materials Science and Engineering，2001，9(5)：371−390.

[30] KANG Qinjun，ZHANG Dongxiao，LICHTNER P C，et al. Lattice Boltzmann model for crystal growth from supersaturated solution[J]. Geophysical Research Letters，2004，31(21)：L21604.

[31] KANG Qinjun，LICHTNER P C，ZHANG Dongxiao. Lattice Boltzmann pore–scale model for multicomponent reactive transport in porous media[J]. Journal of Geophysical Research：Solid Earth，2006，111(B5)：B05203.

[32] YANG Bo，XU Tianfu，DU Yiling，et al. Numerical investigation on the influence of CO₂–induced mineral dissolution on hydrogeological and mechanical properties of sandstone using coupled lattice Boltzmann and finite element model[J]. Journal of Hydrology，2024，639：131616.

[33] CHEN Li，KANG Qinjun，TANG Qing，et al. Pore–scale simulation of multicomponent multiphase reactive transport with dissolution and precipitation[J]. International Journal of Heat and Mass Transfer，2015，85：935−949.

[34] QIAN Yuehong，D’HUMIÈRES D，LALLEMAND P. Lattice BGK models for navier–stokes equation[J]. Europhysics Letters (EPL)，1992，17(6)：479−484.

[35] KANG Qinjun，LICHTNER P C，VISWANATHAN H S，et al. Pore scale modeling of reactive transport involved in geologic CO₂ sequestration[J]. Transport in Porous Media，2010，82(1)：197−213.

[36] GRAY F，ANABARAONYE B，SHAH S，et al. Chemical mechanisms of dissolution of calcite by HCl in porous media：Simulations and experiment[J]. Advances in Water Resources，2018，121：369−387.

[37] GRAY F，ANABARAONYE B U，CRAWSHAW J P，et al. Pore–scale dissolution mechanisms in calcite–CO₂–brine systems：The impact of non–linear reaction kinetics and coupled ion transport[J]. Geochimica et Cosmochimica Acta，2021，305：323−338.

[38] WANG Zi，CHEN Li，WEI Hangkai，et al. Pore–scale study of mineral dissolution in heterogeneous structures and deep learning prediction of permeability[J]. Physics of Fluids，2022，34(11)：116609.

[39] TSCHACHE S，VINJE V，LIE J E，et al. Quantifying amplitude–variation–with–offset uncertainties related to calcite–cemented beds using a Monte Carlo simulation[J]. Interpretation，2023，11(2)：T313−T327.

[40] AL–ABRI A，HIWA S，ROBERT A. Experimental investigation of the velocity–dependent relative permeability and sweep efficiency of supercritical CO₂ injection into gas condensate reservoirs[J]. Journal of Natural Gas Science and Engineering，2009，1(4/5)：158−164.

[41] ZOU Qisu，HE Xiaoyi. On pressure and velocity boundary conditions for the lattice Boltzmann BGK model[J]. Physics of Fluids，1997，9(6)：1591−1598.

[42] YANG Bo，ZHU Huixing，YUAN Yilong，et al. Prediction of CO₂ breakthrough pressure across diverse lithologies：From direct numerical simulation to a new efficient estimation method[J]. Engineering Geology，2025，355：108263.

[43] LOU Qin，GUO Zhaoli，SHI Baochang. Evaluation of outflow boundary conditions for two–phase lattice Boltzmann equation[J]. Physical Review E，2013，87(6)：063301.

[44] MAVKO G，MUKERJI T，DVORKIN J. The rock physics handbook：Tools for seismic analysis of porous media[M]. Cambridge：Cambridge University Press，2003.

[45] SUN Yue，YU Hang，YANG Bo. Impact of wettability on CO₂ dynamic dissolution in three–dimensional porous media：Pore–scale simulation using the lattice Boltzmann method[J]. Langmuir，2024，40(43)：22658−22672.

[46] LIU Min，MOSTAGHIMI P. Pore–scale simulation of dissolution–induced variations in rock mechanical properties[J]. International Journal of Heat and Mass Transfer，2017，111：842−851.