第43卷 第8期 2010年8月 天 津 大 学 学 报 Journal of Tianjin University
V ol.43 No.8.Aug. 2010
收稿日期：2009-03-02；修回日期：2009-09-25.
基金项目：国家“十一五”科技支撑计划资助项目(2006BAJ04B04)；国家科技支撑计划资助项目(2008BAJ12B03). 作者简介：张于峰（1954— ），男，博士，教授． 通讯作者：张于峰，yufengfa@https://www.doczj.com/doc/4d17789675.html,.
U 型埋管系统地下传热数值模拟
张于峰1，陈成敏1，聂金哲1，周小珠1,2，胡晓微1，马洪亭1
(1. 天津大学环境科学与工程学院，天津 300072；2. 钦州学院物理与电子工程系，钦州 535000)
摘 要：采用柱热源模型，建立了无限大区域内U 型换热器与土壤间非稳态传热的二维数学模型．以天津地区的U 型垂直埋管实验得到的数据作为参数，利用FLUENT 软件进行模拟，得出该地区土壤源热泵间歇运行6年的土壤温度分布规律：距离地埋管最近的1m 点其变化速度最快、幅度最大；随着典型点与U 管的距离增加，其变化速度及幅度都将减小；距离地埋管最远的7 m 点其温度几乎不受热泵系统运行的影响，始终维持在初始温度14.5℃左右． 关键词：U 型埋管换热器；土壤源热泵；数值模拟
中图分类号：X57 文献标志码：A 文章编号：0493-2137(2010)08-0717-05
Numerical Simulation of Heat Transfer of Soil Around
U -Tube Undergroud Heat Exchanger
ZHANG Yu-feng 1，CHEN Cheng-min 1，NIE Jin-zhe 1，ZHOU Xiao-zhu 1,2，
HU Xiao-wei 1，MA Hong-ting 1
(1. School of Environmental Science and Engineering ，Tianjin University ，Tianjin 300072，China ； 2. Department of Physics and Electronic Engineering ，Qinzhou University ，Qinzhou 535000，China )
Abstract ：A two-dimensional unsteady conduction model of heat transfer between a U-tube heat exchanger and the
soil in an infinite zone was established by employing the column heat reservoir model. The data obtained from a U-tube performance experiment conducted in Tianjin were used as input ，and the temperature distribution of six-year soil in the same area was simulated by FLUENT software. The temperature of the nearest point ，which is 1 m away from the pipe ，drops the fastest in velocity and greatest in amplitude. The velocity and amplitude decrease as the dis-tance between the typical point and the U-tube increases. The temperature of the farthest point ，which is 7 m away from the pipe ，is free from the impact of the heat exchanger ，remaining at the initial temperature of about 14.5 ℃.
Keywords ：U-tube underground heat exchanger ；ground source heat pump ；numerical simulation
地源热泵系统是一种利用地下浅层资源的既可以供热又可制冷的高效节能空调系统．其工作原理是节能系统通过地源热泵机将地下的热能提取出来对建筑物供暖，或者将建筑物中的热能释放到地下从而实现对建筑物的制冷[1]． 夏季，可将建筑内的热能储存于地层中以备冬用，同样，冬季可以将富余的冷量储存于地层以备夏用．这样，通过利用地层自身的热工特点实现对建筑物和环境的能量交换．地源热泵通过输入少量的高品位能源(如电能)，实现低温位热能向高温位转移．理论上，地源热泵消耗1kW 能量，用户可以得到4kW 以上的热量或冷量[2]．比电
锅炉加热节省2/3以上电能；比燃料锅炉节省1/2以
上的能量．由于地源热泵的热源温度全年较为稳定，天津地区一般为10～15℃左右[3]，其制冷、制热系数可达3.5～4.4，与传统的空气源热源相比，要高出40%左右，其运行费用为普通中央空调系统的50%～60%．因此近年来，地源热泵空调系统在北美国家取得了较快发展，我国的地源热泵市场也趋活跃[4]．地源热泵分为地下水源热泵、土壤源热泵和地表水源热泵等[5]．土壤源热泵技术能否被广泛地推广应用，很大程度上取决于精确、可靠的系统设计方法和计算工具的有效性，地下埋管换热器长期运行性能研究是整
·718·天津大学学报第43卷 第8期个系统的核心部件．土壤源热泵系统运行过程中对地
下土壤温度产生的影响需要进行更进一步的研究．
U型垂直埋管属于土壤源热泵的一种，具有良好
的节能、环保等特性，而且经济效益显著[6]，适用于城
乡居民住所及办公楼等的采暖、制冷需求．笔者以U
型垂直埋管(井埋)为研究对象，建立传热系统数学模
型，利用FLUENT对其进行求解，并模拟出土壤温度
变化情况，为地源热泵在天津地区的应用推广提供必
要的理论指导．
1 U型埋管间歇性运行系统及其原理
以夏季工作原理为例，冬季工作原理只需将图中
冷凝器和蒸发器的位置互换即可．
本实验以国家“十一五”规划新农村过程中采暖
空调技术在农村的应用为基础，以天津静海县独流镇
的一建筑物作为工程对象，进行地源热泵空调系统设
计与安装，其原理如图1所示．系统运行模式为：制
冷-休息-供热-休息，以12个月为周期，每个过程运行
时间均为3个月．由于农村地区的特殊条件：院落可
利用面积较大，故实验中U形管布置间距较大，为
8m，可认为2个U形管之间无传热影响．
图1夏季土壤源热泵工作原理
Fig.1Principle of ground source heat pump in summer 大地换热器数据采集主要采用PT1000型铂电阻温度传感器对地下温度进行测量，采用安捷伦34970A数据采集仪读取并记录测量值，测试得到天津地区地下土壤初温度为14.5℃．其他数据采集及处理见文献[7]．
2 垂直埋管系统传热数学模型
2.1 基本设定
U型管换热器用一个等效单管换热器代替，等效
管管径为1/2
ed0
()
D n D
=[8-9]，0D为U型管外径，对单管，n＝2；对双U型管，n＝4；土壤与埋管之间的传热方式为纯导热，且土壤按深度分层，每层土壤的导热系数为常数，土壤与埋管接触良好，忽略接触热阻；忽略土壤中因水分迁移而引起的热迁移；等价管的管壁与土壤的换热量各处是均匀的；与垂直埋管底部具有相同深度的土壤层上下之间无热交换；土壤被视为无限大、均匀、各向同性介质；没有邻近管井的热干扰；不考虑灌浆材料的热影响；地源热泵机组间歇性运行；埋管内有对流换热热阻．
2.2 传热数学模型
对于常热流密度的圆柱热源[10]，埋管井周围无限土壤介质中温度场分布的表达式为
()
g r0
s
,
q
T T T G z p
k L
∞
Δ=?=(1)
2
01
22
11
10
1e1
(,)[()()
()()
1
()()]d
z
G z p J p Y
J Y
J Y p
β
ββ
ββ
βββ
?
∞
2
?
=?
π+
∫
式中：s
2
t
z
r
α
=；
r
p
r
=；如果假设土壤的热物性是不变的，则()
G z p
，仅是时间t和半径r的函数．直接计算()
G z p
，十分困难，可利用计算线图直接查得．
对于U型埋管，采用当量直径
eq
D，且有
eq0
D=(2)
埋管管壁热阻为
00i
p
p
ln(/)
r r r
R
k
=(3)埋管内对流换热热阻(光滑圆管，紊流)为
f
i0i
1
(/)
R
h r r
=(4)
式中：f
i i
i
2
k
h Nu
r
=，0.8
i i
0.023n
Nu Re Pr
=，当0.4
n=时用于制热工况，当0.3
n=时用于制冷工况．
从埋管外壁到管内流体的当量导热热阻为
eq
p f
1
h
R R
=
+
(5)
将
p
R和f R代入式(5)得
1
00
i
eq
0i p i
ln
r r
r
h
r h k r
?
??
??
=+
??
??
??
??
??
(6)
埋管外壁与管内流体的温差为
p
eq0
q
T
h A
Δ=(7)
00
2
A r L
=π
U 型埋管系统地下传热数值模拟 ·719·
2.3 初始及边界条件
埋管所处区域不同度深度土壤初始温度一致且不考虑地面传热，初始温度为14.5℃；在地面与空气接触的地方，认为计算区域的土壤与外界空气不进行换热，使用绝热边界条件；模型两边边界设置为距U 型管中心20m 处、下部边界150m 处不受土壤传热的影响．土壤远边界温度随深度和以年为单位的时间发生变化，该温度可使用解析方程组确定，再使用加权平均计算整个耦合的平均远边界温度． 2.4 结果修正
对于单U 型埋管，由于两腿间距离很近，通过管壁的热流密度不均匀，对式(7)进行修正得
p 0t eq 0q
T C N h A Δ=，其中，0C 为热流不均匀分布修正系
数，单U 型管取0.85，双U 型管取0.6；t N 为钻孔内埋管根数，单U 型管取2，双U 型管取4． 埋管内水的平均温度为
wa g p T T T T ∞=+Δ+Δ (8) 由能量平衡得流体进出口温差为 wo wi p
q
T T mc ?=
(9) 定义流体平均温度为 wo wi
wa 2
T T T +=
(10) 代入式(9)得流体出口温度为 wo wa 2p
q
T T mc =+
(11) 考虑U 型埋管两腿间热短路的修正则有 wo wa sc 2p
q
T T T mc =++Δ (12) 式中 ()sc 2
sc
2p q
T mc R ?Δ=
sc f pw s 2R R R R ′=++ f 4
3i i R r h L
′=
π ()
0pw p 4ln /3i r r R k L
=π
[]100s s cos (2)/2x r r R k L
?+=π
考虑热泵机组的间歇运行，将时间分隔使用不同的热流密度，计算式为
{
g sc 111s 1
(,)(,)n n T R F q G z p G z p k L
???Δ=?+?? sc 2212(,)(,)n n R F q G z p G z p ?????++?? sc 11(2)(1)(,)(,)n n n n n n R F q G z p G z p ?????????+??
}
sc 1(,)n n R F q G z p ???? (13)
3 U 型地下埋管循环系统参数
U 型管外径为30mm ，根据等效管径计算公式eq 0D =可得到等效直径为42mm 的圆柱；U 型管
采用高密度聚乙烯管(即PE 管)，等效直径为
42mm ，管中循环液为水；工程所在地区的土壤为重土-潮湿类型，密度为20943kg /m ，比热为962J /kg K ?（
），导热系数为1.3W/(m K)?；管中循环水温度取平均值：制冷时24.08℃，供热时5.85℃，流量为2.03m /h ．循环水与U 管接触的壁面按对流换热进行计算，循环水与管壁之间的对流换热系数的计算为
0.8f 0.023/n h Re Pr d λ= (14) 式中：Re 为雷诺数，地源热泵系统制冷时取2.98×104，供热时取1.75×104；Pr 为普朗特数，地源热泵系统制冷时取6.22，供热时取11.595；n 为地源热泵系统制冷时取0.4，供热时取0.3；f λ为导热系数，地源热泵系统制冷时取0.6049W/(m K)?，供热时取0.5625W/(m K)?；d 为等效直径，d ＝0.042m ．
4 模拟及其结果分析
网格划分采用GMBIT 软件建立等效圆柱体的几何模型，并以0.5m 为步长划分长方形结构格图，所有边界条件设置为WALL 类型；选择2DMesh 类型输出网格．计算时间步长采用3 600 s ，地源热泵系统运行模式为：制冷-休息-供热-休息，以12个月为周期，每个过程运行时间均为3个月，共对地源热泵系统运行6年间对地下温度场产生的影响进行模拟．利用FLUENT 进行仿真计算，求解器采用非耦合求解法和隐式算法，按非稳态传热设置[11]，选择能量方程，其他设置按照默认值．地下温度场的模拟结果，如图2～图7所示．
由图2～图7可以看出，除第1年到第2年期间受地埋管传热影响的地下土壤范围增大之外，其余几年温度场几乎没有变化．
·720· 天 津 大 学 学 报 第43卷 第8期
图2 系统运行1年后温度场 (单位：K)
Fig.2
Temperature field after one year ′s operation (unit: K)
图3 系统运行2年后温度场(单位：K)
Fig.3 Temperature field after two years ′operation (unit: K)
图4 系统运行3年后温度场(单位：K) Fig.4 Temperature field after three years ′
operation (unit:K)
图5 系统运行4年后温度场(单位：K)
Fig.5 Temperature field after four years ′operation (unit: K)
图8是以曲线的形式描述地下温度场各典型点的温度变化情况图．由图中可以看出，在土壤耦合热泵系统第1个运行周期之内，地下温度场随热泵系统运行工况不同呈现出不同的变化趋势；同时，在各典型点温度变化过程中，距离地埋管最近的1m 点其变
化速度最快、幅度最大，随着典型点与双U 管的距离增加，其变化速度及幅度都将减小，距离地埋管最远的7 m 点其温度几乎不受热泵系统运行的影响，始终维持在初始温度14.5 ℃左右．
图6 系统运行5年后温度场(单位：K)
Fig.6 Temperature field after five years ′operation (unit: K)
图7 系统运行6年后温度场(单位：K)
Fig.7 Temperature field after six years ′operation (unit: K)
图8 系统运行第1个周期内地下温度场各典型点温度变化 Fig.8 Temperature change of typical point in under - ground temperature field during the first operation cycle
5 结 语
计算结果显示，土壤耦合热泵系统运行的每个周
期之内，地下温度场随热泵系统运行工况不同呈现出不同的变化趋势，系统制冷时地下温度升高，系统制热时地下温度降低，系统停止运行时地下温度呈现出向初始温度回归的趋势，而距离地埋管最近的点受热泵系统运行的影响最强烈也最直接，随着距离地增加，影响将逐渐减小．计算结果同时显示，在系统每个运行周期结束之后，地下温度相对周期开始时刻都
2010年8月 张于峰等：U 型埋管系统地下传热数值模拟 ·721·
是降低的．而随着温度场中各点与地埋管距离地增加，其温度降低幅度是逐渐减小的，并且发现在系统运行的第1个周期，各点温度降低速度最快，而后随着运行周期地增多，温度变化速度明显减缓． 符号说明：
0A —埋管与土壤间的传热面积，2
m ；
p c —循环流体定压比热，J/(kg K)?； n F —制冷或制热运行份额； f h —对流换热系数，2W/(m K)?；
eq h —从埋管外壁到管内流体的当量热传递系数，
2W/(m K)?；
0J ,1J —零阶和一阶第一类贝塞尔函数；
p k —管壁导热系数，W/(m K)?；
f k —流体导热系数，W/(m K)?； s k —土壤导热系数，W/(m K)?；
L —竖井深度，
m ； m —循环流体流量，kg/s ；
n q —地下埋管换热量，
W ； q —地下埋管换热量，W ； i r —埋管内半径，
m ； 0r —地下埋管外壁的半径，m ；
r —距离原点的半径，
m ； p R —管壁的导热热阻，
K/W ； f R —流体与内管壁对流换热热阻，
2m K /W ?（）； f R ′—考虑热短路的流体与内管壁的对流换热热阻，
2
m K /W ?（）；
pw R —考虑热短路的管壁的导热热阻，K/W ； s R —考虑热短路的土壤的导热热阻，
K/W ； t —运行时间，
s ； wa T —埋管内平均水温，K ；
wi T —埋管进口水温，K ； wo T —埋管出口水温，K ； T ∞—远边界地温，K ；
r 0T —地下埋管外壁的温度，
K ； g T Δ—远边界地温(原始温度)与地下埋管外壁的温度的差值，K ；
g T Δ—埋管外壁温度与远边界处大地温度的差值，K ； sc T Δ—U 型埋管两臂间热短路的温度修正，K ； p T Δ—埋管外壁温度与管内流体平均温度的差值，K ； 0Y ，1Y —零阶和一阶第二类贝塞尔函数； Nu —努塞尔数；
Re —雷诺数； Pr —普朗特数；
x —U 型管两臂的间距，m ；
s α—土壤的导温系数，2m /s ．
参考文献：
［1］ Huttrer W. Geothermal heat pumps ：An increasingly
successful technology [J ]. Renewable Energy ，1997，10(2/3)：481-488.
［2］ Omer A M. Ground-source heat pumps systems and ap-plications [J ]. Renewable and Sustainable Energy Re-views ，2008，12(2)：344-371.
［3］ 马最良，姚 杨，杨自强. 水环热泵空调系统设计
[M ]. 北京：化学工业出版社，2005.
Ma Zuiliang ，Yao Yang ，Yang Ziqiang. System Design for Pendular Water Ring Heat Pump Air-Conditioning [M ]. Beijing ：Chemical Industry Press ，2005(in Chi-nese ).
［4］ 徐 伟. 中国地源热泵情况调查与分析[J ]. 工程建设
与设计，2006(12)：16-19.
Xu Wei. Survey of domestic GSHP application [J ]. Con-struction and Design for Project China ，2006(12)：16-19.
［5］ 刁乃仁，方肇洪. 地埋管地源热泵技术[M ]. 北京：
高等教育出版社，2006.
Diao N airen ，Fang Zhaohong. The Technology of Ground-couple Heat Pump [M ]. Beijing ：Higher Educa-tion Press ，2006(in Chinese ).
［6］ Hepbasli A ，Kalinci Y . A review of heat pump water
heating systems [J ]. Renewable and Sustainable Energy Review ，2009，13(6/7)：1211-1229.
［7］ 聂金哲. 村镇地下储能直接冷却技术研究[D ]. 天津：
天津大学环境科学与工程学院，2008.
N ie Jinzhe. Villages and Small Towns Underground Stored Energy Direct Cooling Engineering Research [D ]. Tianjin ：School of Environment Science and Engineer-ing ，Tianjin University ，2008(in Chinese ).
［8］ 刘正华，陈汝东，李 芃，等. 土壤源热泵系统埋地
换热器换热性能研究[J ]. 流体机械，2007，35(3)：63-67.
Liu Zhenghua ，Chen Rudong ，Li Peng ，et al. Per-formance study of underground heat exchanger in soil source heat pump system [J ]. Fluid Machinery ，2007，35(3)：63-67(in Chinese ).
［9］ 张开黎，王如竹，于立强. 垂直埋管土壤源热泵供热
供冷实验与分析[J ]. 流体机械，2001，29(9)：57-60. Zhang Kaili ，Wang Ruzhu ，Yu Liqiang. Theoretical re-search of vertical buried ground source heat pump [J ]. Fluid Machinery ，2001，29(9)：57-60(in Chinese ).
［10］ Deerman J D ，Kavanaugh S P. Simulation of vertical u-tube ground-coupled heat pump systems using the cylin-drical heat source solution [J ]. ASHRAE Transactions ，1991，97(1)：287-295.
［11］ 陶文铨. 数值传热学[M ]. 2版. 西安：西安交通大学
出版社，2001.
Tao Wenquan. Numerica Heat Transfer [M ]. 2nd ed. Xi'an ：Xi'an Jiaotong University Press ，2001(in Chi-nese ).