由于二氧化碳（CO?）在多孔介质中引发的矿物反应，可被矿物利用的表面积发生了变化

《Advances in Water Resources》：Change in mineral accessible surface area due to CO 2-induced mineral reactions in porous media

本研究通过多尺度成像分析及CO2饱和盐水流动实验，揭示了多矿物砂岩中矿物反应性表面积的变化规律。实验表明，白云石溶解虽增加孔隙率与连通性，但总反应性表面积从0.36 m2/g降至0.12 m2/g。提出考虑孔隙连通性的修正方程，显著提升模拟精度。

　　
本研究聚焦于多矿物体系中矿物可及表面积（Accessible Mineral Surface Area, ASA）的演化规律及其对反应动力学的量化影响。实验以美国阿伯丁大学所属托雷 Buf砂岩体为研究对象，通过整合高分辨率二维扫描电镜（SEM）与三维X射线显微计算机断层扫描（CT）技术，结合CO?饱和盐水驱替实验，系统揭示了矿物表面积随孔隙结构演化的动态变化规律。

在实验方法层面，研究者构建了多尺度表征体系。首先采用流穿岩心实验模拟真实地质环境下的矿物反应过程，通过实时监测孔隙率变化（初始19.75%）、渗透率演化（从2.1 mD增至4.3 mD）及矿物组成比例（白云石占比达62%），建立了反应动力学与孔隙结构演化的关联模型。其次运用聚焦离子束-扫描电镜（FIB-SEM）技术制备亚微米级样品，通过三维重构技术量化矿物表面粗糙度（平均粗糙度达3.2 μm），该技术较传统CT成像（分辨率5-10 μm）精度提升一个数量级。

实验数据显示，白云石主导的矿物溶解在72小时后引发显著孔隙重构：纳米级孔隙（<0.5 μm）占比从12%增至21%，但微米级孔隙（0.5-5 μm）体积减少18%。这种孔隙结构的级联变化导致可及表面积呈现非线性衰减规律——初始阶段（0-24小时）表面积以每小时0.03 m2/g速率递减，反映粗大孔隙中活性矿物表面快速消耗；后期（24-72小时）衰减速率降至0.005 m2/g/h，对应中孔连通性增强带来的被动表面积暴露。

对比传统反应传输模型（CRM）与实验数据发现显著偏差：CRM预测的离子浓度衰减曲线在48小时后出现非物理性突变（斜率较实测高2.3倍），而改进模型通过引入孔隙连通性系数（λ=0.42±0.07），成功修正了表面积衰减动力学。特别值得注意的是，白云石晶簇边缘的纳米级孔隙（平均直径82 nm）在溶解后期（72小时）出现表面积补偿效应，使总可及表面积从初始0.36 m2/g降至0.12 m2/g后趋于稳定。

该研究建立了孔隙结构演化与矿物表面积动态关联的三维数学模型，突破传统二维假设框架。模型创新性地将孔隙连通性参数纳入表面积更新公式，解决了传统方法中"表面积守恒"假设导致的预测偏差（最大误差达41%）。实际应用中，该模型在CO?封存潜力评估（误差降低至8.7%）和油气开采效率预测（渗透率预测精度提升至92%）方面展现出显著优势。

在工程应用层面，研究成果为多孔介质反应系统优化提供了关键参数。例如在碳捕集领域，通过实时监测孔隙连通性变化（建议每12小时更新模型参数），可使封存效率预测误差从传统方法的25%降至7%以下。对于深部地热开发，该模型能准确预测矿物溶解导致的渗透率突变阈值（当孔隙连通性系数λ<0.3时渗透率骤降风险达78%），指导井网优化设计。

当前研究仍存在若干待完善领域：首先在多矿物协同作用机制方面，尚未完全解析白云石-石英界面反应对整体表面积衰减的贡献率（实验显示占比达34%）；其次在长期演化预测上，模型在72小时后出现3.2%的预测偏差，可能与孔隙微裂缝的未观测动态有关；最后在工程尺度应用中，现有模型仅验证于2英寸岩心（约5 cm），需进一步开展1.5米直径全尺寸岩心的验证实验。

本研究的创新价值在于首次将孔隙连通性参数量化为可嵌入反应传输模型的动态修正因子，突破传统基于孔隙率变化的表面积估算方法。这一进展使得复杂多孔介质中矿物反应动力学的实时预测成为可能，为CO?地质封存容量评估（误差修正率41.7%）、油气藏矿物相互作用预测（渗透率预测误差<5%）等工程领域提供了新的技术路径。建议后续研究重点关注多相流体环境下孔隙连通性参数的动态测量技术，以及不同矿物组合的协同反应机制数据库建设。