应力腐蚀和压力溶解引起的裂隙开度-刚度变化对thmm耦合过程影响的数值模拟

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