盐岩非线性流变性态及地下复杂溶腔变形数值分析

针对平顶山盐田盐岩层数多、单层厚度薄、泥岩夹层多的特点,建立等效复合单元代替互层状盐岩体,利用ABAQUS软件对马庄井田的地下溶腔进行三维数值模拟研究,探讨不同内压下,腔周塑性区围岩的非线性流变性态以及溶腔体积缩减规律。结果表明:(1)等效复合单元在稳态蠕变阶段具有较好的数值稳定性,泥岩能有效地限制黏塑性区范围以及围岩向洞内位移的发展;(2)在围压(8 MPa)作用下,互层状盐岩单元体基本呈现线性小变形特征,溶腔体积收缩率与蠕变持续时间均呈现最低值,分别为3%与2.5 a,随着内压增大(8~12 MPa),互层状盐岩单元体开始呈现非线性蠕变,但经过较短时间后(2.5 a)体积收缩率逐渐稳定至一较小值,约为3.5%,随着内压的迅速减小(≤4 MPa),腔周盐岩逐渐呈现非线性大变形特征,表现为稳态蠕变时间增长,蠕变速率增加,5 a后围岩蠕变仍在稳定增加;(3)塑性区围岩内的黏滞系数随时间普遍加大,体现了互层状盐岩在储气运行过程中越来越稳定的变形趋势,因此,实际运行时储气压力应控制在合理范围(4~12 MPa),尤其要减少在较小内压下的蓄气运行时间。