地源热泵系统已经成为了应用最广泛的“绿色”空调系统，在恰当的管理下，能达到高效和节能的效果，估计全世界有超过110万个已经运行的地源热泵系统 [13]。但多数建筑冬夏冷热负荷不同，从而引起地埋管换热器的吸热和放热不平衡，多余的热量或冷量堆积在地下，土壤温度偏离初始温度。目前，对地源热泵系统热不平衡的研究多数集中引起热不平衡的原因及解决措施上，即复合式系统及运行控制策略的研究[46]。
摘 要：以夏热冬冷地区某实际地埋管地源热泵系统为分析对象，对夏季采用冷却塔供冷而仅冬季采用地埋管和消防水池联合供热的系统运行参数进行了6 a测试。建立了地埋管三维管群模型，通过数值计算方法对地埋管周围岩土的温度分布进行了热平衡分析和计算。通过测试数据与理论计算结果进行对比分析，得到了影响大地自调节能力的影响因素。
关键词：核心期刊征稿投稿,地源热泵,大地自调节能力,热平衡,三维模型
受太阳辐射、大气循环、降水等因素的作用，土壤与大气进行着能量的交换导致土壤温度周期性的日变化和季节变化，其趋势呈现为竖向上的梯度变化[7]。地源热泵系统利用土壤作为系统的低位冷热源，即土壤本身具备一定的能量，可以称为大地能。大地能通过与外界的能量交换调节土壤温度，使得土壤温度按照一定的规律变化。对于地埋管地源热泵系统，建筑负荷通过地埋管向土壤释热或吸热，对原始土壤温度的变化形成一定的外扰，使得土壤温度升高或降低。当输入到土壤中的热量控制在一定的范围内，则土壤温度变化值可以控制在一定范围内。当一个供冷(供热)季结束后，下一个供热(供冷)季开始前，通过大地土壤上下层热量的扩散、地表与大气之间热量传递和辐射换热等自然换热条件下，使得土壤温度能够逐渐基本恢复到初始温度。这种大地抗外扰的能力可以称为大地的自调节能力。
在大地蓄热能力分析中，大地的自然调节能力是不能忽略的[8]。为此，必须正确了解地源热泵系统热平衡率的本质。地源热泵系统利用地埋管进行取热和放热，释放到大地中的热量或冷量因大地的蓄能能力而被储存在大地中。由于蓄冷或蓄热的不平衡可能会导致多年地温偏移初始温度，而影响系统的正常运行。但这种不平衡是以超过大地的自调节能力体现出来的。若不能正确理解热平衡的意义，在工程实施中可能会增加埋管量或其他技术措施来保证大地的热平衡，但这会增加地源热泵系统的初投资，同时也不能保证地源热泵的高效运行。由于地埋管换热是一个复杂的非稳态过程，土壤温度受诸多因素的影响，并不能简单的以建筑冷热负荷的差异来判断土壤热平衡 [9]。大多学者以土壤温度的恢复情况来判断土壤的热平衡[1011]，即地源热泵系统运行一个循环周期后土壤温度能够恢复到初始温度则表示系统是热平衡的。影响土壤温度恢复的因素很多，包括对土壤的物性参数、室外气候条件、埋管布置形式、系统的运行方式以及建筑负荷强度等[1213]。由于多种因素的存在，对热平衡的研究显得较为复杂。但研究大地的自调节能力是热平衡研究的基础，因此有必要对大地的自调节能力进行研究。
王 勇，等：冬季工况下地埋管地源热泵系统中大地的自调节能力分析
在夏热冬冷地区，夏季冷负荷大于冬季热负荷，地埋管地源热泵系统是夏季向土壤放热而冬季向土壤取热，研究的重点集中在夏季地埋管换热研究和复合式地源热泵系统的优化与控制研究[1416]，对于只有冬季向土壤取热的地源热泵系统的运行情况研究较少。本文主要是针对夏热冬冷地区某实际工程，在夏季采用冷却塔供冷，而仅冬季采用地埋管取热的地源热泵系统供热的运行模式下，对系统运行参数进行了分析。通过该系统多年的冬季运行数据分析和土壤温度实验测试，并建立了地埋管换热器三维管群模型，利用FLUENT软件计算，讨论和分析大地对释冷量的自调节能力范围。
1 地源热泵工程概况
某医院用地源热泵系统位于夏热冬冷地区的重庆市，该建筑有11层，1～10层主要为病房及医生办公室，11层主要为手术室和医生办公室。夏季全楼最大冷负荷1 650 kW，其中手术室冷负荷208 kW;冬季全楼最大热负荷为759 kW，其中手术室热负荷138 kW。地埋管系统采用DN32聚乙烯管，孔深80 m，设计流量128 L/s，流速0.66 m/s，共240个换热孔，孔间距3 m。
该地源热泵系统2006年冬季开始运行至今。图1为医院地源热泵系统图，11层手术室采用水水热泵机组，热负荷138 kW，其余均采用水空气热泵机组，其中消防水池起减压作用。从负荷大小可以看出，手术室夏季冷负荷相对较小，对地埋管换热的影响较小，可以忽略。根据建筑使用情况，可以认为夏季采用冷却塔供冷，冬季采用地源热泵系统供热。图2为地埋管换热器各测点布置图，地埋管外壁面温度测点布置在进出水管-1、-9、 -29 m，以及出水管-59 m处。根据该建筑地源热泵系统的实际使用情况，只有冬季使用地源热泵系统，建筑通过地埋管与土壤的热交换将冷量释放到大地，系统只向大地取热，这与一般夏热冬冷地区地源热泵系统的使用情况不同。为此，本文主要讨论该运行模式下土壤温度的变化，从而分析大地对冷量的自身调节能力。根据DEST软件建模计算该建筑逐时负荷，冬季累积热负荷902 532.92 kWh，累积单位孔深热负荷47 kWh/m，该建筑的的运行热负荷大部分集中在40%～80%的范围。
2 实验数据分析
2.1 多年运行数据分析
项目于2008年建立了较为完整的系统运行数据资料，包括系统的启停状态和末端用户供水温度。根据原始记录数据整理，从2008―2011年冬季供暖期间(12月、1月、2月)，冷却水泵一直连续运行，而地埋管循环水泵则间断开启。冬季1月份热负荷最大，选取1月份典型日循环水泵启停情况以及末端用户机组供水温度，见图3。 从图3可以看出，冬季运行期间，末端用户供水温度均在15℃以上，22℃以下，大部分时间保持在17～20℃。地埋管侧循环水泵间歇运行，随着循环水泵的开启，供水温度上升，且可保持在一个较高的温度范围，但在不同时间段内，其提供的最高温度有所不同，这主要受系统形式以及室外空气温度的影响。
从系统形式上看，由于系统中增加了消防水池，其830 m3的水容量有一定的蓄热作用，地埋管环路水与大地之间交换得到的热量则先储存在消防水池中。当消防水池温度达到20℃后，该温度已经接近土壤的温度，若继续开启地埋管循环水泵，地埋管换热效果差，同时增加了水泵能耗，该段时间循环水泵停止运行。当末端用户供水温度低于18℃，即消防水池中的水温低于 18℃，此时启动循环水泵，经过地埋管换热，消防水池中的水温升高，升高到20℃后循环水泵又停止运行。与常规的地源热泵系统相比，加入消防水池后，这种系统运行方式可以提高末端机组的供水温度，使得机组的效率大大提高，但进入地埋管水的水质要求提高。同时由于系统形式的特殊性，循环水泵间歇运行，也为土壤温度提供了恢复时间，管壁周围土壤温度不至于过低而导致系统效率降低。
为了分析室外空气温度的影响，典型年重庆市1月1日―1月22日室外逐时空气干球温度变化见图4。如1月16日―1月22日，17日开始降温，到21日早上降到最低5.0℃。系统运行情况如图3(c)，17日地埋管循环水泵停止运行后，因为气温降低，建筑热负荷增大，消防水池中储存的热量很快使用完，用户进水温度直接降到17℃，而后几天内，地埋管循环水泵均开启以维持用户进水温度16～18℃，同时因为昼夜温差的不同，用户进水温度也有所波动。由此看出当室外空气温度降低时，系统所需提供的热量随着建筑热负荷的增大而增大，由于消防水池的蓄能能力有限[17]，消防水池储存的热量在较短时间内被末端用户置换完毕，需要开启地埋管以提高用户供水温度。
2.2 土壤温度分析
根据测试2006年8月土壤初始自然温度为213℃。2012年1月7日系统运行了1个月后，管壁周围平均温度18.9℃，见图5。此时地埋管侧流量为 201 m3/h，进出水温差约1～1.5℃，即地埋管提供的热量约233～350 kW，约为DEST软件计算出的建筑热负荷的31%～50%。地埋管承担的热负荷比设计热负荷小，根据实际使用情况调查，主要有以下几个原因：1)手术室的使用时间根据手术安排，并不是24 h连续开启;2)实际医院入住率基本保证在100%，大部分科室走道中也增加了床位，同时白天探病的亲人较多，人数大大的增多，从而增加了人体散热量，使得热负荷减小;3)医院新风主要靠卫士间排风和自然通风，冬季病人一般不喜欢开窗，则新风热负荷减小;4)根据不同病人的情况，一些病人手术后因为供暖时空气干燥而容易缺水，则并未使用。
2012年3月1日早上8点地埋管循环水泵停止运行，管壁周围温度开始恢复，测试是从中午11：50开始，各点管壁温度恢复情况见图6。由于土壤温度的恢复在初始时刻变化最大，在早8：00地埋管循环水泵停止运行时刻管壁温度约18.5℃。从图中可以看出，各点温度缓慢上升，12 h后，到晚上20：00各点温度约为19.5～19.8℃，比岩土初始自然温度低1.7～1.5℃。
根据多年系统运行数据分析和土壤温度分析，系统运行一个循环周期后，管壁周围土壤温度能够基本恢复到初始温度，冷量没有堆积在管壁周围，且能保持系统长期高效运行。这是因为冬季建筑通过地埋管向土壤吸收的热量没有破坏大地的自调节能力。由此看出，在这种运行模式下，当地埋管向土壤的取热量在一定范围内时，土壤温度可以通过大地的自调节能力，使其管壁周围温度恢复到初始温度。地埋管周围岩土的自调节能力的研究可以通过数值计算方法进行分析确定。
3 大地自调节能力分析
3.1 地下换热器传热模型建立
3.1.1 物理模型 地埋管与大地的传热是一个复杂的非稳态传热过程，线热源模型最简单而又较实用[1819]。但线热源模型忽略了地表对土壤温度的影响以及竖向温度的变化，为了分析大地对温度的自平衡能力的大时，需要考虑上表面的换热、土壤远边界条件的精确设置以及短时间内计算的准确性，需要采用三维模型[2021]。本文采用了三维管群模型，为了简化，作如下假设：
1)忽略地表面辐射换热对岩土温度的影响，认为岩土初始温度均匀一致，为当地的冬季岩土初始温度。
2)在整个传热过程中岩土的物理成分、热物性参数不变。
3)仅考虑纯导热，忽略各种接触热阻。
4)不考虑水分迁移对热量传递的影响。
5)忽略底部U型管与土壤间的换热。
管群三维模型土壤长×宽×高为24 m×20 m×80 m，共12个钻孔，钻孔直径0.130 m，孔间距4 m。U型管外径为0.032 m，内径0.026 m，管道中心距离为0.06 m，如图7。由于U型管底部弯管处换热量较小，此处忽略不计。采用软件Gambit软件进行模型的建立及划分网格，同时因为管壁周围对换热影响较大，对地埋管周围网格进行了局部加密，总网格个数为1 800 160，如图8。同时各种材料的物性参数见表1。
3.1.2 数学模型
1)控制方程
土壤的传热模型是三维的非稳态传热，式(1)为直角坐标系中的传热方程。
U型管内流体流动是紊流流动，kε双方程模型是紊流粘性系数模型中应用最广泛和最成功的一种模型。管中流动水为粘性常数的不可压缩流体，且不考虑重力作用，式(2)～(6)为简化后的控制方程。
2)边界条件和初始条件
进口：设置为速度进口，进口速度为0.44 m/s，进口温度根据DeST软件计算得到的逐时负荷Q和每一步迭代得到的出水温度t2，由式(9)计算出相应进水温度t1。 Q(1-1COP)=CPm(t1-t2)(9)
出口：设置为压力出口，出口温度由每一步迭代求得。
进出水管底部：进水管底部设置为压力出口;出水管底部设置为速度进口，速度为0.44 m/s，且其进水温度等于进水管底部出水温度。
上表面：设置为固定壁面，采用第2类边界条件，边界面周围空气温度为重庆市典型年室外逐时干球温度。边界面与空气之间的表面传热系数根据式(10)[11]计算，v为低空附近空气速度。
h=5.8+3.7v(10)
下表面及远边界：设置为固定壁面，采用第1类边界条件，温度采用测试得到的冬季岩土初始温度19.5℃。
钻孔孔壁及U型管管壁：钻孔孔壁设置为固定壁面，壁厚度为0 mm;U型管管壁设置为耦合壁面，管壁厚度为3 mm，从而进行U型管内流体、回填材料以及岩土的耦合传热计算。
3)计算工况及时间步长
为了分析土壤温度自身调节能力的大小，讨论不同累积热负荷大小和不同负荷强度情况下土壤温度的变化情况，计算工况见表2。
工况累积热负荷比例负荷强度比例累积热负荷/kWh累积单位孔深热负荷/(kWh・m-1)平均单位孔深换热量/(W・m-1)
冬季运行12月、1月和2月，1 d运行12 h，停机12 h，恢复期为3月和4月，时间步长为12 h，共302步，采用FLUENT软件计算。
3.2 计算结果和讨论
在分析土壤温度变化时，由于在短期内埋管竖直方向温度变化较小，且地埋管换热主要集中在地埋管上部，则以地下30 m处土壤温度为代表，不同工况下冬季运行3个月后和恢复2个月后土壤温度分布见图9。
从图9可以看出，工况1下，运行3个月后钻孔之间基本没有相互影响，且经过2个月的恢复期后，埋管区域温度基本恢复到了初始温度。而工况3下，运行期间钻孔之间相互影响大且恢复2个月后，埋管中间区域温度仍然较低，即工况3下所需土壤温度的恢复时间更长。为了分析土壤温度随运行时间的变化情况，根据模型边界条件的设置选取管群中心最不利的钻孔1和远边界换热最有利的钻孔2为分析对象。钻孔1和钻孔2地下30 m处孔壁平均温度随运行时间变化见图10。
从图10可以看出：工况1，冬季地源热泵系统运行3个月后孔壁温度缓慢下降到18.0℃，比初始温度仅降低了1.5℃，同时恢复速度也较慢，恢复2个月后孔壁温度到达18.6～18.8℃，与初始温度相差不到1℃，钻孔1和钻孔2温度接近，这是因为承担的热负荷小，冷扩散半径小，钻孔之间没有受到相互热干扰;工况3，承担负荷较大，孔壁温度下降速度快，在运行3个月后，孔壁温度降低到了14℃，比初始温度降低了5.5℃，但其恢复速度也较快，恢复2个月后两个钻孔分别恢复到了16.2℃和16.9℃;工况2孔壁温度介于工况1和工况3之间，运行3个月后，孔壁温度降低到16.6℃，而恢复2个月后能分别恢复到18.1、18.5℃，此时钻孔2的温度已经接近工况1下孔壁温度。
经过以上分析看出热负荷大于工况2后，土壤温度下降速度大。在工况3下，冬季运行仅20 d后，孔壁温度就降低到了16.6℃，此时累积热负荷为182 193.32 kWh，还小于工况1的3个月累积热负荷值。这是因为当热负荷强度增大到一定值后，需要向土壤吸收的热量增大，破坏了大地的自调节能力，使得冷量堆积在管壁周围，即在短时间内土壤不能温度恢复到接近初始温度。但对于工况3，大地的自调节能力有所破坏，在供热末期孔壁温度在13～14℃，此时地埋管出口水温在11～12℃，仍然能保证地源热泵系统的正常运行。
对比系统实际运行工况和计算结果，实际运行工况介于工况2和工况1之间，在这种间歇运行的情况下，可认为工况2为大地自调节能力的临界点，平均单孔吸收的热量21.76 W/m，大地的自调节能力能够得到有效的保证。因此，从以上分析可以看出，自调节能力大小不仅与累积热负荷大小有关，而且与逐时的负荷强度有关。
4 结 论
1)间歇运行模式为地源热泵系统正常运行条件中较重要的影响参数。保证间歇运行的技术方案较多，而利用消防水池的蓄热能力来降低地埋管的运行时间，是值到考虑的技术方案。通过该医院住院部地源热泵系统的测试和分析发现，在冬季，利用消防水池蓄热能力来保证土壤温度恢复时间，从而提高地埋管的换热效率，该方案有利系统长期高效运行。
2)大地的自调节能力是研究地源热泵系统高效运行中不可忽略的一个重要方面。忽略该参数会导致盲目采用辅热或辅冷等其他技术措施来保证其正常运行，其不必要的技术措施会造成初投资或运行费增加。因此，大地的自调节能力分析的准确性有利保障地源热泵的技术体系。从该医院多年测试数据以及理论计算结果分析，当负荷强度以及累积热负荷限定在一定的区间内时，不会造成热平衡失调而导致的地温过度变化。对于该项目的具体情况，其单孔承担的热量控制在21.76 W/m左右，大地的自调节能力能够保证埋管周围的岩土温度在季节周期内恢复到初始温度，保证大地温度的恒定。
3)从计算分析以及测试结果可以看到，大地对释冷量的自调节能力，取决于负荷特征以及土壤的地质结构以及系统的运行模式。在负荷特征中，累积负荷以及逐时负荷强度是其重要影响参数。本文通过案例定义和分析了大地的自调节能力，在不同条件下，大地的自调节能力是不同的，应做具体分析。
参考文献：
[1]Spitler J D.Groundsource heat pump system researchpast， present and future [J].HVAC&R Research，2005，11(2)：165167.
[2]Younis M，Bolisetti T，Ting D S K.Ground source heat pump systems： current status [J].International Journal of Environmental Studies，2010，67(3)：405415. [3]Abdeen M O.Groundsource heat pumps systems and applications [J].Renewable and Sustainable Energy Reviews，2008，12(2)：344371.
[4]於仲义，胡平放，王彬，等.混合式地源热泵系统优化设计[J].暖通空调，2007，39(9)：105109.
Yu Z Y，Hu P F，Wang B，et al.Optimizing design of hybrid groundsource heat pump systems [J].Heating Ventilating & Air Conditioning，2007，39(9)：105109.
[5]Man Y，Yang H X，Fang Z H.Study on hybrid groundcoupled heat pump systems [J].Energy and Buildings，2008，40(11)：20282036.
[6]Hackel S，Neillis G，Sanford K.Optimization of cooling dominated hybrid ground coupled heat pump systems [J].ASHRAE Transactions， 2009：565580.
[7]张慧智，史学正，于东升，等.中国土壤温度的季节性变化及其区域分异研究[J].土壤学报，2009，46(2)：227234.
