地埋管换热仿真从入门到进阶：单管/管群/U型管和螺旋管工程案例

01

如下图（截自参考文献）所示包含九个地埋管的井群（管群）模型，其中地埋管均为双U管（即每孔两根U型管）。

材料参数可参考下表（实际制作时取值与表中略有区别）。

实际仿真使用的模型图（地埋管埋入深度140m，回填土深度142m；计算域为长方体，尺寸为35mX35mX142m）。土壤初始温度为17℃，水初始温度16℃。进口流速0.4m/s，温度5℃，出口回流温度9.6℃。

计算时，考虑土壤的分层和渗流，模型中对各层土壤进行了划分（即将模型中土壤部分切成若干层，如上图）。各层土壤参数可参考下表（根据实际取值）。

地下水渗流以及管孔布置如下图所示，渗流方向与X轴呈30度夹角。

相关信息（摘自文献）：不考虑地下水渗流时，模型中1#、3#、7#、9#处于井群的顶角位置，称为“角井”; 2#、4#、6#、8# 处于边缘部位，称为“边井”; 5#处于井群的中间部位，称为“中井”。

考虑渗流时，各层土壤简化为均一的多孔介质，相关参数可参考下表（根据实际取值）。

1 地下水渗流速度很大时的结果

如第一层的渗流速度为-5次方量级

各埋管间的相互影响很小。

地下3米水平面温度分布

地下30米水平面温度分布

立面温度分布

各层土壤性质及渗流速度不同，温度分布不同。

监测点的温度变化曲线

2 地下水渗流速度很小时的结果

各层渗流速度分别比前一种情况小100-500倍。

从图上可见，渗流速度小时，各埋管间的相互影响很大。

地下3米水平面温度分布

地下30米水平面温度分布

立面温度分布

分层面位置很明显。

这些图形近似对称（有渗流时肯定不完全对称），其实是因为工况时间比较短。若计算的工况时间较长，比如640天结果图如下：

使用软件处理一下比较漂亮

3 假定土壤为各向异性多孔介质

此为假定情况，实际来说肯定不会是这样的。

此处假定土壤中沿渗流速度方向（与X轴正向呈30度夹角方向）渗流阻力较小（与前两种情况两同），其余两个方向上渗流阻力为该方向的100倍。

由如下结果图可见，若假设沿渗流速度方向阻力较小时，云图中的小尾巴更加明显了。