应力-渗流-溶蚀耦合作用下三维岩石裂隙渗透特性数值计算研究.pdfVIP

第47卷第2期岩土工程学报Vol.47No.2

2025年.2月ChineseJournalofGeotechnicalEngineeringFeb.2025

DOI：10.11779/CJG

应力-渗流-溶蚀耦合作用下三维岩石裂隙渗透特性

数值计算研究

112*111

申林方，吕倩文，刘文连，张家明，杨鸿忠，李泽

（1.昆明理工大学建筑工程学院，云南昆明650500；2.中国有色金属工业昆明勘察设计研究院，云南昆明650051）

摘要：基于格子Boltzmann方法采用双分布函数分别模拟渗流速度场与溶质浓度场的演化过程，建立了三维岩石裂隙

应力-渗流-溶蚀耦合作用机制的数值计算模型，并讨论了渗流流速、法向应力、溶蚀反应速率等因素对裂隙渗透特性

演化规律的影响。结果表明：在渗流流速较低时，壁面溶蚀出来的离子得不到及时输运，使得出口处浓度较高溶蚀速

度慢，裂隙结构呈“喇叭口”状。增大法向应力会减小裂隙开度，减慢溶质的运移速率，使得裂隙出口处的溶蚀速率

显著降低，从而限制了其渗透率的发展。当壁面溶蚀反应速率较小时，裂隙渗透率呈持续缓慢增长的状态；随着溶蚀

反应速率增加，出口处的溶蚀量会明显小于入口处，导致出口处壁面发生显著溶蚀之前，裂隙渗透率发展缓慢，此后

渗透率便呈急速突变增长趋势。研究成果能够为酸蚀作用下岩石裂隙渗透能力的定量评价提供重要理论支撑。

关键词：岩石裂隙；应力-渗流-溶蚀耦合；渗透特性；格子Boltzmann方法；数值计算

中图分类号：TU45文献标识码：A文章编号：1000-4548(2025)02-0428-10

作者简介：申林方(1982—)，女，博士，教授，主要从事水岩多场耦合方面的研究工作。E-mail:linfangshen@126.com。

Numericalstudyonpermeabilitypropertiesofthree-dimensionalrockfracture

undercoupledstress-seedissolutionprocess

112111

SHENLinfang,LÜQianwen,LIUWenlian,ZHANGJiaming,YANGHongzhong,LIZe

(1.FacultyofCivilEngineeringandMechanics,KunmingUniversityofScienceandTechnology,Kunming650500,China;2.Kunming

ProspectingDesignInstituteofChinaNonferrousMetalsIndustry,Kunming650051，China)