在使用Yasuhara等建立的裂隙开度的应力腐蚀和压力溶解模型的基础上,考虑裂隙闭合量对裂隙刚度的影响,针对一个假设的位于非饱和双重孔隙-裂隙岩体中且有核素泄漏的高放废物地质处置模型,拟定两种计算工况:(1)裂隙的刚度系数是裂隙闭合量的线性函数;(2)裂隙刚度为常数,进行热-水-应力-迁移耦合的二维有限元数值模拟,考察了岩体中的温度、裂隙开度的闭合速率、闭合量、孔(裂)隙水压力、地下水流速、核素浓度、裂隙刚度和正应力的变化、分布等情况。结果表明,两种工况岩体中的温度场、孔(裂)隙水压力、地下水流速、核素浓度无明显差别;裂隙闭合基本由应力腐蚀产生;在相同计算时间内两种工况的裂隙闭合量较为接近,工况1略大;工况1中离玻璃固化体越近,裂隙刚度值越高,并且在玻璃固化体附近的应力值较大,且集中程度较高。