中深层地埋管群的传热分析.docxVIP

方法开发计算

程序，采用有限差分法求解。算例考虑的系统是由两个钻孔组成的套管式中深层地埋管换热器，两个钻孔的几何参数、物理条件完全相同，埋设在由

3层水平地层组成的岩土中。系统的主要参数如下:钻孔直径d=0.28m、深度H=2000m、钻孔间距D=10m、大地热流qg=0.

075W/m2、单孔循环水流量M=12kg/s.单孔取热速率恒定Q=200kWo在取热负荷不变的条件下模拟了系统运行8000h的传

热工况。采用文章提出的降维算法，可以把三维问题分解为若干柱坐标中的二维问题求解。柱坐标中的二维导热问题在径向采用变步长，纵向由于温

度梯度很小，可以采用较大的步长。整个问题(包括钻孔内)采用的节点数为10752个，完全可以保证计算这样的大型三维传热问题的精度。

经测试，在台式电脑英特尔酷睿(i7-8700CPU@320GHz3.19GHz)上完成以上运算用时10.412s。对算例进行计

算的结果如图4所示。其中，系统初始时刻温度场温度分布如图4⑻所示，系统结束运行时(8000h)岩土中的温度场分布如图4(b)

所示。随着地埋管换热器不断取热，钻孔周边岩土层的温度逐渐降低。靠近埋管区域的岩土温度降低得最为剧烈，远离地埋管换热器的岩土温度逐渐

趋于初始地温。当运行8000h时，双孔之间的岩土温度明显低于外侧岩土温度，说明二者之间存在热干扰现象，这对埋管换热器的取热效率有不

利的影响。图4双钻孔地埋管换热器的岩土温度响应图单孔取热

及双孔取热两个不同工况下，流体进出口温度随运行时间的变化如图5所示。在恒

定取热工况的条件下，单孔及双孔的进出口水温均随着运行时间的增大而逐渐降低。持续取热的时间小于约半年时，两者的进出水温度几乎完全相同

,说明在运行时间足够短时，双孔之间的热影响可以忽略不计。计算结果表明，在第1年的运行期间，两工况的进水和出水温度基本一致，说明在运

行时间足够短的情况下，可以忽略相邻钻孔间的热干扰。这充分表明了降维算法的可信度。随着运行时间的增加，两者之间的差别逐渐加大。在两孔

的间距为10m的情况下，持续运行10年后，双孔和单孔换热器的进出口水温可以相差约6.1℃o这表明两个相邻地埋管换热器之间的热干扰

已经严重影响到了地埋管换热器的效能。（双孔间距10m，持续负荷200kWo）图5双孔及单孔地埋管换热器进出口水温随时间

的变化趋势图5、结论文章基于线性叠加原理提出了多孔地埋管换热器的降维计算方法及地埋管群的传热模型，摒弃了传统地源热泵传热模型中关于钻孔内外区域连接条件的均匀温度和均匀热流的简化假定，因而更准确地描绘了中深层地埋管中的传热过程。主要结论如下：（1）地下岩层中初始温度分布的影响可以由单独的一项叠加到零初始温度的温度响应中。因此，叠加原理可以扩展到考虑地温梯度的场合。（2）单个钻孔地埋管换热器钻孔外的导热通常可以简化为柱坐标中的二维问题;而多个钻孔地埋管换热器钻孔外的导热必须要按三维问题处理。在中深层地埋管群的传热问题中，线性叠加原理可以适用于钻孔外区域的导热问题。因此，地埋管群的传热计算可以转化为多个单孔地埋管换热器传热（二维问题）的叠加来求解，从而大大提高计算的效率。（3）基于有限差分法的单孔中深层地埋管换热器的二维数值计算模型，采用了追赶法求解相应的线性代数方程组，具有计算高效快捷的特点，特别适合在用叠加原理求解多孔地埋管换热器的场合使用。（4）相邻的地埋管换热器在长期运行过程中存在热干扰现象，其影响的大小与钻孔间距、运行时间以及取热速率有关。关于中深层地埋管群的参数分析以及设计和优化运行将在后续的研究中给出。

长（至少为数

月至数年），特别是由多个地埋管换热器组成的地埋管群中的传热问题都必须按三维问题来考虑。因此，中深层地埋管换热器传热的研究是地源热

泵空调系统的技术难点，同时也是该项技术研究的核心和应用的

基础。地埋管换热器的传热分析方法大体上可以分为两大类，即解析解方法和数值解

方法，其各有长处和局限性［3-4］。解析解方法的基本出发点是研究无限大或半无限大介质中单个钻孔埋管的传热规律，此时地埋管换热器简化

为线热源或圆柱面热源，并得到其在恒定热流作用下温度响应的解析解。20世纪80年代，ESKILSON［5］提出了一种叠加原理的新思路

,基于描述无限大介质中单个钻孔在恒定热流加热条件下的温度响应，再利用叠加原理得到多个钻孔组成的地埋管换热器在变化负荷作用下的实际温

度响应。这种方法的物理概念清晰、计算简捷，考虑到了地埋管换热器的复杂几何配置和负荷随时间的变化，已经广泛应用于实际的工程设计计算和

能耗分析［6-8］。这种解析解方法的一个不足是需要假设地下岩土为均匀