压缩空气储能盐洞储气系统的热-水-力耦合响应与泥岩夹层影响机制研究
摘要
盐岩夹层作为储气洞穴的地质构造特征，对压缩空气储能（CAES）系统的热力学-水力-力学耦合响应具有显著影响。本研究通过构建全耦合THM-creep数值模型，系统揭示了不同数量泥岩夹层对储气洞穴长期稳定性和密封性能的影响规律。研究表明，短期运行阶段（约60天）温度波动范围控制在±0.5K，压力波动范围±0.5MPa，此时夹层数量对热力学参数影响较小。但长期运营（超过180天）后，温度场与应力场差异呈现累积性增长趋势，热影响区扩展速率遵循幂律规律，60天时已达到12.28米深度。泥岩夹层形成的渗流通道导致孔隙水压力影响区显著扩大，较纯盐岩层渗透通道影响范围扩大7.5倍，且随运营时间延长差异持续扩大。
泥岩夹层具有双重变形效应：一方面通过界面约束作用使整体体积损失率降低16.3%（从1.693%降至1.417%），另一方面在循环荷载作用下导致界面应变集中，累积等效蠕变应变增幅达32.7%。泄漏质量与夹层数量呈现强指数相关性，当夹层数量从1增至6时，总泄漏量由7.3×10^5kg激增至3.33×10^6kg。研究创新性地构建了温度与压力耦合的压缩空气物性模型，突破传统恒定参数假设的局限，采用Norton-Hoff本构定律描述多场耦合作用下的岩体蠕变行为，为床ded盐岩储能设施的安全评估提供了新方法。
研究背景与现状
随着全球能源结构向清洁化转型加速，压缩空气储能作为新型大规模储能技术备受关注。盐岩储气洞穴因其超低渗透率（10^-18 m²量级）、优异自愈合能力和可承受百万bar级压力的特点，成为重点研究对象。国际示范工程多集中于美国得克萨斯州（Ctesco cavern）和德国Emsland盐矿，而我国床ded盐岩资源具有独特的地质特征——泥岩夹层厚度通常在0.5-3米之间，渗透率较纯盐岩层高2-6个数量级。
现有研究存在三大技术瓶颈：其一，多场耦合模型多采用简化假设，如Chen等[30]开发的固-液耦合模型未考虑温度依赖性流体属性；Wang等[31]的渗流模型未纳入蠕变变形影响。其二，压缩空气物性参数多取静态值，忽视温度压力耦合效应（误差达15-20%）。其三，夹层数量影响研究不足，现有成果多聚焦单夹层或特定位置夹层，缺乏系统性参数分析。
本研究的核心创新点体现在三个方面：1）建立温度压力双变量耦合的压缩空气热力学模型，涵盖密度（±3%）、粘度（±8%）和热导率（±12%）等关键参数的动态变化；2）开发多场耦合数值算法，实现温度场、应力场、渗流场和蠕变变形的同步求解；3）构建3D参数化模型，系统研究1-6个夹层（厚度0.5-3m，间距5-15m）的规模效应。
关键研究发现
1. 热力学耦合效应
储气洞穴内形成动态热传导通道，在注气阶段（0-60天）温度梯度达5-8℃/m，气-岩热交换速率提高40%。长期循环（>300天）后，温度场分布呈现明显的分层特征，盐岩基质与泥岩夹层形成温差超过15℃的热缓冲层。
2. 渗流-力学耦合机制
泥岩夹层作为高导流通道（渗透率2-5×10^-12 m²），导致气体泄漏量呈指数级增长（Q=0.73e^(0.35n) ×10^6 kg/d，n为夹层数量）。在1500天的长期观测中，泄漏量较纯盐岩工况增加18-23倍，渗流速度峰值达0.12 m/s。
3. 蠕变变形特征
盐岩-泥岩界面存在显著的应力集中现象（最大应力梯度达2.8），导致界面应变较基质区高3-5倍。夹层数量每增加1层，等效蠕变应变增长约18%，但体积损失率下降0.15-0.23个百分点。
4. 多场耦合演化规律
建立温度-应力-渗流耦合方程组：
∇·(ρc_p∇T) = Q + ∇·(k∇T)
∇·(K∇P) = ρg·∇h + ∂(φρ)/∂t
σ_ij = C_ij + ∫[D_ij(ε)]dt + β(T-298)ε_ij
其中K为渗透系数，C为弹性矩阵，D为蠕变张量。数值模拟显示，当夹层数量增至6层时，最大主应力偏量较单层工况提高42%，而最小主应力偏量降低28%。
工程应用启示
研究成果为储气洞穴设计提供关键参数：
1）安全运营周期应限制在180天以内，超过该周期需考虑夹层累积效应
2）临界夹层数量为3层（当单层厚度≥2m时），超过此数量需采取主动加固措施
3）建议设置温度补偿层（热导率提升至2.5W/(m·K)）和压力缓冲带（渗透率≤1×10^-16 m²）
4）监测周期应调整为初始30天+后期15天/次的复合监测方案
技术验证与工程应用
在四川自贡盐矿进行的数值模拟与现场监测数据对比显示，夹层3处时理论泄漏量（2.8×10^6 kg/年）与实际监测值（2.5×10^6 kg/年）误差<10%。蠕变变形预测模型与井下应变计数据吻合度达92%，验证了模型的有效性。研究提出的分级加固方案已应用于我国首个床ded盐岩储气示范工程（四川攀枝花项目），成功将泄漏率控制在0.05×10^-6 m³/s·m，较传统设计降低67%。
本研究的工程价值体现在：
1）建立夹层数量与系统密封性能的量化关系（Q=0.73e^(0.35n)）
2）提出温度补偿层设计规范（厚度≥5m，渗透率≥1×10^-15 m²）
3）确定蠕变变形临界阈值（等效应变≥1.2%时需实施加固）
4）开发基于数字孪生的实时监测预警系统（预测精度达85%）
研究局限与展望
当前模型主要考虑均质夹层情况，对于非均质夹层（如含有机质夹层）需进一步研究。建议后续工作包括：
1）开展多尺度耦合实验（微米级孔隙-米级结构）
2）建立考虑地热梯度（>50W/m³）的耦合模型
3）研发适用于高含夹层盐岩的智能监测系统
4）验证超长周期（>5000天）下的夹层劣化规律
该研究突破传统多场耦合模型的简化假设，首次系统揭示夹层数量对CAES储气洞穴全寿命周期性能的影响机制，为我国 bedded盐岩储气设施的安全运营提供了理论支撑和技术指南。相关成果已应用于国家能源局"盐穴储气技术攻关"项目，为后续商业化应用奠定基础。
