饱和含水冻土区埋地管道水热耦合数值模拟.pdfVIP

2011年 11月 西安石油大学学报 (自然科学版) NOV．20l1 第26卷第6期 JournalofXianShiyouUniversity(NaturalScienceEdition) Vol_26No．6 文章编号：1673-064X(2011)06-0066-04 饱和含水冻土区埋地管道水热耦合数值模拟 付俊鹏，马贵 阳 (辽宁石油化工大学 石油天然气工程学院，辽宁抚顺 113001) 摘要：通过对冻土区管道运行研究，提出针对管道安全运行的措施．采用有限容积法，得到多年冻土 多孔介质水热耦合数学模型，地表面温度采用随年周期性变化条件，应用SIMPLER算法对模型进 行数值求解．通过对无保温层和有相同厚度两种保温材质的管道在春季、夏季和冬季的土壤温度场 进行数值模拟，显示在地表温度波动的条件下，热流密度随土壤温度波动呈现周期性变化．在长期 运行管道中，无保温措施的管道周围冻土融化剧烈，管壁热流密度大且振幅波动大．使用厚度为40 mm的两种保温材料中，40mm聚氨酯保温效果较好．冻土区运行管道应加敷导热系数较小的保温 材料，可极大降低融化深度，保护管道安全运行． 关键词：冻土区；埋地管道；保温层；水热耦合；数值模拟 中图分类号：TE867 文献标识码：A 在温度场与水分场耦合问题的研究中，经常假 动量 怛 、能量万程守恒方程如 卜： 设土壤在空间上是分层均质的，土壤含水量较小，均 (1)质量守恒方程为 匀分布各向同性．冻土层作为隔水层，无地下水补给 +v ．( )．：0． (1) 及其他边界的水分补给和排泄作用，各土层含水量 稳定．但实际多年冻土具有较强的源汇作用，水分场 式中：lf为速度矢量，m／s；t为时间，s；p为流体密度 kg／m。 是不可忽略的因素 1]．以往只考虑其单一输运，未 ． (2)动量守恒方程为 讨论二者耦合效应．土壤中温度分布不均，不仅导致 热流，而且影响水分流，尤其在多年冻土区，土壤温 +V ·(p~uxu)= O +V ·(z／gra )+ x 度引起水的相变更加明显 j．水分迁移同样伴随热 量流．分析水热场耦合迁移运动，需建立符合实际的 ( 1II)+绛 (2) 水热耦合方程． +V．(p州 一 ap+ ·(x／gra )+ 1 数学模型 ( c1II)+ + 以下为饱和冻土多孔介质水热耦合控制方程． pfgot(T一 ef)． (3) 土壤孔隙是指孔隙中完全充满水时水的质量与 式中-U，U 分别为 U在 ，Y方向上的速度矢量， 固体颗粒质量之比．假设土体各相均质连续，相变过 程流体密度变化符合boussinesq假设，低速渗流时， Ⅱs；=150南 为多孔介质渗透率，m；c·= 水分迁移符合达西定律，忽略由相变融化引起的速 度变化 4。．根据有限容积理论，建立的质量守恒、 ×{ 为惯性损失系数，m一；8为孔隙率，无 收稿 日期：2011-06-20 基金项目：辽宁省 自然科学基金项 目“冻土区埋地管道周围