深井回灌式水源热泵井群运行的地下含水层传、蓄热性能模拟研究(1)

1、    深井回灌式水源热泵井群运行的地下含水层传、蓄热性能模拟研究(1)    目前国内进行深井回灌式水源热泵工程的井群设计和施工过程中，系统方案的可行性判据基本取决于单井出水量是否满足求，以及能否实现良好的人工回灌。然而良好的设计还需考虑井群当地地下含水层的水热运行与水文地质条件、环境气象因素和工程措施之间的关系。笔者通过比较目前流行的含水层流动传热模拟程序，选择利用了美国地质调查局编写的HST3D程序，对一典型双井承压含水层的温度场和流场进行了全年运行模拟，对该程序应用于此类问题的功能性和适用性作出评价，指出其

2、需完善之处。关键词 深井回灌 水源热泵 含水层水热运动 热贯通 建筑容积率 HST3D     问题的提出深井回灌式水源热泵技术作为一种有益于环境保护和可持续发展的冷热源形式，在国内外空调工程界已经得到了越来越多的应用1 2，文献3给出了其基本原理与相关技术经济分析。这一系统方式利用温度全年相对恒定的地下水作为水源热泵的水源，通过建造抽水及回灌井群，实现夏季抽冷水、灌热水，冬季抽热水、灌冷水的这一全年角色轮换的运行过程，地下含水层内部的热量或冷量被提取、蓄存和转移。井群是深井回灌式水源热泵系统的一个关键组成部分，其正常运行与否决定了应用水源热泵系统工程的

3、成败，井群的设计布局应当是慎之又慎的关键环节。目前国内进行此类工程的井群设计和施工过程中，系统方案的可行性判据基本取决于单井出水量是否满足求，以及能否实现良好的人工回灌。然而在进行该类工程井群部分的可行性分析和设计中，还需考虑以下几方面的问题：（1）当地含水层中的能量蓄存、转移过程。应用深井回灌方式，需在设计阶段知道当地含水层的能量提供能力有多大，也就是系统可负担的建筑容积率极限是多少。当建筑物全年冷热负荷不均匀时，系统对于含水层温度的常年影响效果如何，是否会造成含水层"背景温度"逐年降低或升高，从而导致系统运行失败，相关的应对策略如何制定？（2）"热贯通&quo

4、t;影响的避免由于回灌水与原始含水层温度存在的差异，在导热和对流等作用下，回灌井水"温度锋面"会导致近抽水井出水温度有不同程度的升高或降低，通常称为"热贯通"现象。如何确定适宜的井间距，如何确定井群的布局，避免"热贯通"的影响，是设计人员关心的主问题。对于高密度住宅小区或城区商用建筑应用深井回灌式水源热泵系统来说，由于可利用建筑用地的面积限制，如何优化井群布局及其各自对应的抽水或回灌角色，最大限度地避免"热贯通"的不利影响是尤为关键的。（3）水文地质条件的影响作用该问题的核心是如何考虑速度相对较大的当地地下水自然流

5、动的存在对于地下含水层温度场和井群布局的影响，如何在建筑物冷全年热负荷不均匀的情况下，利用自然地下水流场的存在，合理地优化各井的抽水、回灌角色和轮换方式，从而实现对于能源的最优利用。（4）地面机组和管道系统形式和运行模式的影响作用利用"小流量、大温差"的系统运行方式，能够实现对于含水层蓄能的最大利用，同时减少对于地下水资源最小程度的开采利用3。在"小流量、大温差"和传统的"大流量、小温差"两种运行工况下，地下含水量水层温度场全年变化过程，以及所导致对应的井群布局差异如何，也是研究设计人员所关心的问题。以上这些问题都需寻求对于井群部分的

6、含水层水热运动过程适用的计算分析工具，研究井群当地地下含水层的水热运动与水文地质条件、环境气象因素和工程措施之间的关系，为该类系统的相关设计与分析提供科学的依据，以推动集中式水源热泵机组应用的进一步推广。在此工程背景下，笔者通过比较目前工程学术界流行的含水层流动传热模拟程序，选择利用了美国地质调查局开发的地下水流动、传热、传质三维有限差分模拟程序HST3D，对一典型双井承压含水层进行了全年温度场和流场模拟，对该程序应用本问题的功能性和适用性作出评价，指出其需完善之处。2 含水层中的水热运动及相关数值模拟21 含水层中的水热运动自然界含水层中的地下水流动一般满足达西定律，而含水层内部的传热过程包

7、括：（1） 地下水的对流换热过程；（2） 地下水的导热过程；（3） 固体骨架的导热过程；（4） 由于通过多孔介质孔隙的不同流动通道液体的机械混合造成的局部热弥散，以及由于不同地质成分构造混合所造成的宏观热弥散过程；（5） 地下水与固体骨架之间的传热。在含水层传热中，当含水层骨架颗粒较小和流体流动雷诺数较低时，可以假定流体温度场与固体骨架的温度场时一致的，因此可以不考虑上面第5项的热量传递。第4项由含水层多孔介质热弥散机理所造成的热量传递，对于含水层内部的整体传热过程存在着不可忽略的影响，特别是单井周边含水层中由于相对较高的地下水孔隙流速，热弥散的影响作用尤为突出。22 含水层水热运动数值模拟含

8、水层水热运动的建模和模拟计算工作，在含水层季节性热蓄能等相关领域已经进行了较长时间的深入研究。含水层季节性热蓄能的应用和研究，早于上世纪七十年代中期已经相继在我国、北欧、北美等地区和国家开展。美国加州大学的Lawrence Berkeley Laboratory (LBL)建立了相应的单井含水层蓄能有限差分数值模拟程序CCC，并通过该程序对现场实验进行了模拟和分析4。国内陈兆祥5和薛禹群6等亦完成了相关模拟计算和现场实验工作。国际能源组织于1993年至2000年期间所完成的IEA ANNEX 8： Implementing Underground Thermal Energy Storage，

9、其中一个子项目就是关于地热蓄能的设计分析工具应用与评价，G?ran Helstr?m 7在其为该子项目所作的总结报告中列出了适用于含水层水热运动分析现行软件，包括AST、TWOW、SUTRA、Tradikon、HST3D等程序。Chiasson8的论文中亦列出了适用于含水层水热运动分析的现行软件，除以上所列以外，还包括SWIFT、AQUA3D、FEFLOW等。下面将对HST3D程序作原理和功能介绍。3 HST3D简介            摘 目前国内进行深井回灌式水源热泵工

10、程的井群设计和施工过程中，系统方案的可行性判据基本取决于单井出水         本篇论文是由3COME文档频道的网友为您在网络上收集整理饼投稿至本站的，论文版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否者后果自负，如果此文侵犯您的合法权益，请联系我们。HST3D910是英文Three-Dimensional flow，Heat, and Solute Transport model 的简称，它是美国地质调查局（USGS）于80年代末开发的一套开放型研究用程序。HST3D采用控制容积的能量平衡法对三维

11、流动、传热和传质微分方程进行离散求解，能够实现饱和含水层中流动、传热和传质过程的非稳态模拟，可用于饱和地下含水层相关流动、传热和传质问题的模拟，包括热田和土壤热、海水入侵、放射性核废料填埋等问题。HST3D具有很强的实用性，其性能包括程序设计结构的模块化、离散方法的简单化和求解方法的多样化、允许采用多种坐标系及不等距网络等，其开放性的模块化结构信纸科研人员可以根据需添加、修改或删除相应的模块。HST3D所求解的流动、传热以及物性方程11分别如下：饱和含水层的流动微分方程： （1）饱和含水层的传热微分方程： （2）假定密度为压力和温度的函数，其液体物性方程： （3）HST3D对于离散方程的系数矩

12、阵的求解方法包括：（1）三对角直接求解法（2）逐次超松弛迭代法（3）基于红黑排序的通用共轭梯度法（4）基于D4Z排序的通用共轭梯度法。HST3D能够处理第一、第二和第三类边界条件，能够处理点源和面源问题，在三维网格坐标方向允许设置不同的土壤传导参数及容积参数。能够处理承压含水层问题，以及存在自由水面的潜水含水层问题。我们选择HST3D作为问题计算分析工具，为进一步评价该程序应用于本问题的功能性和适用性，利用 HST3D对一典型近似工况下的双井承压含水层的全年温度场和流场进行了模拟。4 双井承压含水层模拟选取区域为长300m×宽200m×厚30m的具有上下不透水层的双井承压含

13、水量水层为计算模型（见图1），中部承压含水层以及上下不透水层的厚度均为10m，相关水力热力参数见表1。模型区域中央为一抽水井和一回灌井，两井相距100m，两井均为完整井（透水井壁空越整个承压含水层）。模型的初始温度为15，ABCD、EFGH边界面为15恒温边界条件，ABFE、CDHG边界面为定水头边界条件，ADHE、BCGF边界为不透水边界。模型设定为无自然水头条件。图1 承压含不层计算模型示意图模拟模型的含水层水力热力参数 表1         承压含水层    不透水层&

14、#160;   单位    渗透率    5.3×10-11    1×10-12    m2    孔隙度    0.25    0.35        固体骨架可压缩系数 &#

15、160;  4.6×10-4    4.6×10-4    Pa-1    固体骨架比热容    696    696    J/(kg·)    固体骨架的热传导系数    2600    

16、;2600    kg/m3    纵向弥散率    4    0    m    横向弥散率    1    0    m 为了尽可能模拟水源热泵系统的全年"大温差，小流量"的运行工况，首先进行持续50天的夏季工况运行

17、，抽水及回灌流量均为1200t/d，回灌水温恒为25；然后为持续50天的过渡季工况，两井停止运行；最后为持续50天的冬季运行，抽水及回灌井轮换角色，流量均为1200t/d，回灌水温恒为6。需强调的是，由于HST3D程序的输入功能限制，本算例不得不采用固定流量和固定回灌温度。该算例的模拟结果如下：由图2可以看出，夏季工况期间，回灌热水锋面已经到达抽水井，出现"热贯通"现象，同时部分回灌热量以导热为主的方式向上下不透水层传递。由图3可以看出，过渡季工况期间，热量传递过程以导热为主，大部分夏季回灌热量在含水层内部实现"跨季节"储存。由图4可以看出，冬季工况期间

18、，回灌冷水锋面同样已经到达抽水井，但是由于夏季回灌热水的存在，在一定程度上缓解了冷水锋面对于抽水井出水温度的影响。图2 25连续回灌50天（夏季工况）含水量水层中心剖面温度分布，A井回灌，B井抽水图3 停止运行连续50天（过渡季工况）含水层中心剖面温度分布，A和B井停止运行图4 6冷水回灌连续50天（冬季工况）含水层中心剖面温度分布， A井抽水，B井回灌由图6可以看出，夏季工况期间，抽水井出水温度逐步上升，由原始含水层温度15升至约19；冬季工况运行开始阶段，由于夏季回灌热量在含水层中的蓄存，抽水井出水温度远高于含水层原始水温，达24，可以看出在冬季运行期间实现了部分夏季回灌热量的"

19、热回收利用"。图5 双井运行期间含水层中心平面流速分布示意图图6 冬夏季抽水井出水水温变化曲线5HST3D适用性评价笔者认为，作为适用于深井回灌式水源热泵系统井群部分含水层水热运动的模拟分析的工具，需能够实现以下几方面的功能：（1）能够处理非稳态问题，能够读入动态的边界条件参数（2）能够实现对于井群的参数设置（3）能够处理热边界条件（4）能够实现压力场（水头）的计算（包括单井水头和远端边界的水头影响）（5）能够反映符合工程实际的真实的物理过程，能实现地面系统与井群两部分的联合运行工况分析。通过利用HST3D进行双井承压含水层的模拟，可以看出HST3D可以满足以上所列前4项功能，能够实

20、现一定水文地质条件下含水层水热运动的非稳态模拟，能够给出在一定的井群运行工况下能量在含水层中转移、蓄存的过程分析，能够给出抽水井的动态温度变化。但是HST3D还不能够实现地面以上部分（热泵机组）和地面以下部分（井群及共周边土壤）的全年联合运行分析，也就是说不能够根据建筑物全年动态负荷变化判断回灌温度和水量，从而对于符合工程实际上的过程进行模拟。由于HST3D具有一定的开放性和可拓展性，进一步的工作将为添加相关的建筑物负荷模块，完善其作为深井回灌式水源热泵系统井群运行的地下含水层传、蓄热性能模拟研究的计算分析工具。       

21、;     摘 目前国内进行深井回灌式水源热泵工程的井群设计和施工过程中，系统方案的可行性判据基本取决于单井出水         本篇论文是由3COME文档频道的网友为您在网络上收集整理饼投稿至本站的，论文版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否者后果自负，如果此文侵犯您的合法权益，请联系我们。    摘 本文主介绍和评析几种以铜管为通水部件、抗蚀能力强、能适应我国一般供暖水质使用求的新型铜铝散热器，适合于

22、在住宅供暖中使用，可供同行在工程设计和产品开发时参考。关键词 住宅 铜铝 散热器 评析     2，其中较大部分在采暖区，年需散热器约为1.5亿片（折合铸铁四柱-600型散热器）。加上旧房改建及老供暖系统更换散热器，数量会更多一些。由于建设部已明令自2000年10月1日起实施按户计热与控制，并且城市供热方式还是以城市集中供热为主，所以适合热水供暖使用的散热器仍然是我们研究的重点。已于2001年4月1日实施的国家采暖通风与空气调节设计规范采暖部分，对铝制散热器的内防腐及铸铁散热器内腔不含粘砂等问题提出了明确的求。在现有的各种散热器中，铸铁散热器推广和采用内

23、腔无粘砂新工艺；钢（板）制散热器沿需解决人防腐技术；铝制散热器按即将实施的建工行业标准求，必须采用人防腐措施，以能适用于水的pH12,氯离子含量120ppm的供暖系统。而新型铜铝散热器是以铜管为通水部件、充分发挥铜材耐腐蚀的优点；以铝材做成散热部件，充分发挥铝材密度小、易成型等优点，可以生产在型式、尺寸、重量等方面都能适合住宅使用求的产品。本文将重点介绍我国目前兴起的几种新型铜铝散热器，以供同行选用或开发产品时参考。2 住宅对散热器的求一般说来，住宅用散热器的求为"安全可靠、轻薄美新"。"安全可靠"包括热工性能及制造质量两大方面。求在规定的压力条件下不漏

24、水、在与普通钢管近似的有效使用年限内安全工作且散热量无大的衰减、以及不会对人体造成损伤等方面；"轻"指采用钢、铝、铜、塑等薄壁轻型材料生产散热器，以减轻散热器自重，以及在热工性能提高的基础上减少金属用量；"薄"指散热器的厚度（宽度）尽量小一些，以少占平面面积，希望厚度在40100mm范围内；"美"指散热器的外形、尺度、比例、色彩能与建筑装修协调且有良好的表面处理；"新"造型新颖及新材料、新工艺、新技术的采用。我国是发展中国这，住宅有高、中、低档的判别，所使用的散热器也应当有不同的档次以满足不同标准住宅的求。但上述的

25、一般求则是一致的。对此求，新型铜铝散热器有极大的适应性。面对我国的现实情况，将散热器的通水部件和散热部件区别对待，各用其材分而治之，可能是发展我国散热器生产、保证安全使用的新出路。3 新型铜铝散热器的主优点本产品的突出特点就是将散热器中通风部件与散热部件分别处理、发挥铜材耐蚀特长及铝材重量轻、导热好、易成型的特点，具体优点分述于后：31 以铜管作通水部件，耐蚀能力强，可以保证使用寿命高于供暖系统所用普通钢管的使用寿命。按国家标准工业循环冷却水处理规范求（GB50050-95），碳钢管壁的年腐蚀率应小于0.125mm/a；铜、铜合金及不锈钢管管壁的年腐蚀率应小于0.005mm/a。亦即铜管的耐蚀

26、能力约为碳钢管的25倍。即是用壁厚0.5mm的铜管也相当于壁厚12.5mm碳钢管的耐蚀年限。32 以铜管为通水部件，可以达到较高的承压能力，提高散热器的工作压力。按照国家标准铜及铜合金拉制管（GB/T 1527-1997）的规定，对T2、T3牌号的铜管，产品的壁厚范围为0.510mm，出厂液压试验压力的计算公式为： P=2·S·t/(D-0.8t) (MPa)式中，t为壁厚（mm）；D为管外径（mm）;S为允许应力（纯铜为41.2MPa）。按此公式计算求得不同管径及壁厚铜管的出厂试压值见表1，据此可以根据散热器工作压力的求选用不同壁厚的铜管。常用铜管出厂试验压力计算值（MP

27、a） 表1      外径    壁厚（mm）    0.5    0.6    0.7    0.8    0.9    1.0    1.1    15   

28、0;2.82    3.40    3.99                    16    2.64    3.19    3.74     &#

29、160;              19    2.22    2.67    2.13                    20

30、0;   2.10    2.53    2.97                    22    1.91    2.30    2.69  &#

31、160;                 25    1.67    2.02    2.36    2.71                27    1.55    1.86    2.18    2.50