v2一种考虑地下水渗流的地埋管三维传热模型.docVIP

考虑地下水渗流的地 埋管地源热泵三维数值模拟模型及其应用 胡平放1﹡胡磊2 雷飞1 孙启明1 徐菱虹1 1、华中科技大学环境科学与工程学院 2、中信建筑设计研究总院有限公司 摘要： 本文建立了一种考虑地下水渗流的地源热泵地埋管三维传热模型，并利用模型在Matlab环境下模拟了在无渗流、有渗流及改变地下水渗流速度时，埋管出水温度的变化以及埋管周围温度分布的变化，说明地下水渗流有利于埋管换热，在系统设计中是应该考虑的因素；借助所建模型，通过反演方法，由埋管出水温度推导出土壤综合导热系数，该方法考虑地下水渗流因素，可数值求解综合导热系数，结果可用于指导实际工程设计。 关键词： 地下埋管换热器； Matlab程序； 地下水渗流； 土壤综合导热系数 如何使地源热泵系统设计更为合理，运行更为高效，成为研究者关注的最重要的问题。地源热泵的地下埋管深度一般在60~200m，在这么深的地层内，或多或少都存在地下水的渗流。地下水渗流无疑会对地埋管换热器的传热产生影响。 日本有相关文献用二维数值模型检验了井内温度随抽热速度和地下水流速变化的关系。计算结果表明：对于岩土中地下水渗透流速大于10-5m/s以上者，与假定只有导热的系统设计比较，可以达到节省埋管长度的目的[1]。 美国的Eskilson[2]利用了Carslaw[3]等给出的移动线热源问题的稳态解析解，讨论了在达到稳定状态以后渗流对地热换热器的影响。另外Chiasson[4]等人也利用有限元法数值求解了二维的渗流问题，对一些相关问题进行了计算和比较，但未能找到较一般的规律和结论。 Gehlin 等 讨论了坚硬岩石垂直缝隙中地下水渗流对岩土热物性测试结果的影响，指出即使较小的地下水渗流也会明显加强地下的传热[5]。 哈尔滨工业大学的范蕊、马最良指出目前国内外关于垂直U 型埋管换热器的传热模型都是基于纯导热的模型，主要是由于考虑地下水渗流时情况复杂，分析比较困难，但认为地下水渗流利于埋管的换热[6，7]。王庆鹏采用整体求解方法求得管内流体、地下埋管换热器及周围岩土的温度场数值解。并且通过分析求解结果，研究了地下水渗流对传热过程的影响，表明地下水流动对原温度场的影响明显[8]。 本文建立了一个考虑地下水渗流的地埋管三维传热模拟模型。应用该模型对地下水渗流因素对地埋管流体温度、土壤温度，土壤导热系数计算结果的影响等进行分析。 1 考虑地下水渗流影响的地下埋管三维传热模型 1.1 传热模型的建立 热渗耦合传热的能量控制方程 对于多孔介质土壤，考虑地下水渗流的热渗耦合传热能量控制方程已经由前人推导得出了，孔隙率为，单相流体非等温渗流的能量方程为[6]： (1) 而且，; ; ; ; . 其中，下标表示液相参数；下标表示固相参数；下标表示两相综合参数；表示密度；表示比热容；表示源项能量；表示热容比；表示热扩散率；代表控制体内能量变化项；代表地下水渗流能量变化项；代表导热项；代表源项。将上述能量控制方程进行离散处理，然后对地下埋管及周围土壤进行网格划分，针对每个节点建立能量平衡方程。 对于微元体而言，不考虑内热源的情况下，在微元时间段⊿τ 内，根据控制体单位时间的内能变化量⊿U 与地下水渗流引起的能量变化⊿Qu等于控制体在空间上相邻节点单位时间传入的热量之和，能量平衡关系式建立如式2： 图1 水平截面的网格划分 图2 垂直方向的网格划分 模型中网格用直角坐标（x，y，z）来表示，以埋管支管间节点为中心，分进出水口端节点，两端对称节点坐标相同。为了适合整体网格的划分方式，水平网格划分时，圆形的埋管被等价于截面积相等的方形网格，管段壁厚按照原壁厚选取，管道沿中心线剖开，下端面为对称绝热边界，网格从埋管节点开始向外以1.3倍的比例因子逐步加宽；垂直方向网格按等距划分，在埋管底部 K=Nver-1 的节点层以下，随运行时间的延长，埋管周围土壤的热量积聚程度逐渐加重，考虑到该层节点与其下方的土壤之间有温差的存在，会发生垂直方向的换热，因此在该层以下加设了一层土壤节点（即图中的 K=Ne的节点层）。 1.3 模型假设条件 （1）、将土壤看成一个均匀的、均值的、刚性、各向同性的饱和多孔介质，具有相同的导热系数值；流体与固体瞬间达到局部热平衡；假设水和砂土的热平衡是瞬间发生的，即砂粒骨架和周围流水具有相同的温度；液相通过对流作用输送热量； （2）、液相流动在层流范围内，符合Darcy定律；只考虑单方向的地下水渗流情况； （3）、钻井区域回填材料与原状土壤