核废料处置库近场的热-水-力耦合分析.pdfVIP

第 27卷第 4期 安 徽 工 程 大 学 学 报 Vo1．27．NO．4 2012年 12月 JournalofAnhuiPolytechnicUniversity Dec．，2O12 文章编号 ：1672—2477(2012)04—0075—05 核废料处置库近场的热一水一力耦合分析 张大斌 ，李子善。 (1．安徽理工大学 土木建筑学 院，安徽 淮南 232001； 2．芜湖市建筑工程管理处，安徽 芜湖 241000；3．芜湖市审计局，安徽 芜湖 241000) 摘要 ：为研究高放射性核废料地质处置 的热一水一力耦合过程 ，以FEBEX原位试验为计算模型，利用有限元软 件 code-bright进行数值模拟分析 ，得到热一水一力耦合作用下处置库关闭后近场膨润土和岩石内温度 、饱和度 、 吸力、应力及位移 的变化规律 ，其结果可为核废料处置库 的规划 、设计 以及缓冲／回填材料 的选取提供参考 ． 关 键 词 ：膨润土 ；核废料处置库 ；热一水一力耦合 ；数值模拟 中图分类 号 ：TB12 文献标识码 ：A 在高放废物地质处置库中，一方面库 内核废物衰变放热导致处置库周围缓冲 ／回填材料一膨润土产生 温度应力 ，使膨润土内水受热转化成水蒸汽，并且发生 向外 的扩散 ；另一方面 围岩高压力水头的作用使周 围岩石内水向膨润土 内发生渗透 ，使得高压实的、富含蒙脱石 的膨润土 的饱和度持续上升，并产生膨胀应 力。这样隔离屏障在热量和水作用下要产生非常复杂的 “热一水一力 (T—H—M)”现象，并且它们之 间相互关 联耦合．为了搞清楚这种耦合现象的作用机理 以及保障隔离屏障的长期安全性 ，张玉军 提出热一水一应力 耦合模型并编制了相应 的二维的有 限元程序 ，以FEBEX原位试验为模拟对象进行 了数值分析 ；田辉 。利 用高放核废料处置概念库 的模型，应用有限元软件 ABAQUS对缓冲材料一非饱和膨润土的均匀化过程进 行模拟研究，并对非饱和膨润土进行热一水一力耦合分析 ，得出不同时刻 的温度 、饱和度、孔隙水压、应力和 位移 的空间分布云图和随时间变化曲线 图；XavierPintado等 。以实验 (基准实验)和数值模拟的方法研 究了不同热力梯度条件下膨润土 内的热一水一力耦合现象；M．T．Zandarin等 采用有限元软件 CODE— BRIGHT针对瑞典废物 回取试验 (CRT)进行数值模拟分析 ，并预测 了缓冲材料与铜罐和缓冲材料与岩石 相接触处 的热一水一力耦合现象等． 本文选用适当的热、水、力本构方程，采用 CODE—BRIGHT软件嘲对处置库近场进行热一水一力 (T—H— M)数值模拟 ，得到处置库关闭前后近场膨润土和岩石的温度、饱和度、吸力 、应力和位移的变化规律． 1 本构方程 1．1 热本构方程 本研究假定热传导服从 Fourier定律[6 ．i一一 T，其 中，i为介质导热量 ； T为介质温度 的微 分 ；为热传导系数 ，表征材料导热性能的物性参数 (越大 ，导热性能越好)，显然和三相热传导系数 、 、 l-sz 、 空隙率 以及饱和度有关 ，具有以下关系式 ：— s卜 (1S 一 (一 ) (一 ：)卜 一一 ( ·A — ， 其中，S 为土的饱和度 ， 、 分别为完全干燥和完全饱和情况下土的热传导系数． 1．2 流体本构方程 本研究假定孔隙中液相 、气相遵循 Darcy定律 。．q一一K (vp一pzg)，q一一K (vp 一 g)，其 中，qz、q分别为液相和气相的流量；K 为液或气相渗透系数张量 (K 一忌K ／u，i— l，g)；U 为液相或气 相的动力粘滞系数k一为液相或气相相对渗透系数 ；g为重力 向量 ；K为介质固有渗透系数 ；VP、 P 分 别为水压力和气压力的梯度 ；、Pg分别为液相和气相的密度． 其 中与以上计算有关 的参数表达式如