霍润科，钱美婷，李曙光，熊爱华，丁凡，张茹萍

(1. 西安建筑科技大学 a.土木工程学院；b. 陕西省岩土与地下空间工程重点实验室，西安 710055；2.中铁二十局集团有限公司 高原隧道施工技术及装备研发中心，西安 710016)

随着现代工业的迅速发展，酸性环境污染已经成为工程建设中不容忽视的破坏因素，作为一种边坡工程中常见的非均质材料，砂岩内部含有大量随机分布的微裂隙、孔洞、界面等缺陷，这些不连续的缺陷构成了砂岩多相复合的结构特征，从而影响砂岩的宏观物理性质，使其力学性能呈现各向异性、非线弹性和时效性[1]。随着CT技术的快速发展，其应用范围逐渐拓宽，学者们开始利用CT扫描技术探究受酸腐蚀砂岩截面的物质组成、孔隙结构等岩心参数，分析酸性环境变化对砂岩腐蚀过程的影响。

早期数字岩心技术多选用数值建模方法，主要基于电镜的二维扫描图，通过随机法或过程模拟等数值模拟算法实现三维重构[2]。Øren等[3]提出过程模拟法，考虑岩石基质的粒径分布，模拟真实岩心形成的过程，重现岩石真实的传导性质和几何属性；Hazlett在Metropolis和Kirkpatrick模拟退火算法的基础上，将目标函数定义为重建介质与模拟介质性质差值的平方和，通过体系更新，使其逐渐稳定，从而实现数字岩心；Okabe等[4]提出多点地质统计法，巧妙运用地质信息，利用二维岩心薄片，解决数字岩心重构过程中的连通性难题，随后孙建孟等[5-7]开展孔隙网络流动模拟研究。在Berea[8-10]数字岩心的理论基础上，刘洋等[11]将数字岩心技术应用于研究低渗储层岩石的渗流机理，分析流体在油湿和水湿两种环境下的饱和状态，并判断其饱和度的区别；……