0 引言
随着世界经济的持续发展,人类对能源的消耗量日益增加。在过去的几十年里,人类的能源消耗量增加了一倍[1]。目前,化石能源依旧主导着世界能源结构,化石能源的大量消耗给环境带来了巨大压力,因此有必要进行节能以减少化石能源的消耗。目前普遍认为建筑节能是各种节能途径中最有效且潜力最大的方式[2]。地源热泵（ground source heat pump, GSHP）作为一种地热能利用形式,在世界范围内得到了广泛的应用,被应用于建筑的供暖以及提供生活热水等。地热换热器（borehole heat exchanger, BHE）是地源热泵系统的一个重要组成部分。最常见的BHE是将换热管置于钻孔内,并利用回填材料填充U型管与钻孔壁之间的环空来保持钻孔稳定以及强化传热[3,4,5,6,7]。除此之外,也有利用地下水填充钻孔的BHE形式（groundwater-filled borehole heat exchanger, GFBHE）,这种BHE通常在强岩层区域进行钻孔,地下水通过岩石的裂隙自然地流入钻孔中[8]。GFBHE结构简图如图1所示。在灌浆填充的钻孔中,主要的传热方式是导热;在由水进行填充的钻孔中,主要的传热方式是自然对流传热,有时还受到含水层地下水平流的影响[9]。利用地下水填充钻孔由于不需要填充材料可以减少初投资,相比水泥灌浆也能减少建设时间[10]。并且已有的研究也表明这种BHE具有可靠的传热性能[11,12],甚至还有人提出了在密封良好的钻孔中注入人造液体取代传统灌浆的新型BHE结构[13]。
Fig. 1 GFBHE structure diagram图1 GFBHE结构简图 |
作为一种有效的手段,数值模拟的方法也被应用到GFBHE的研究中[9,11-12,17-19]。GUSTAFSSON等[11]利用数值模拟软件FLUENT建立了一个3 m长的BHE模型,其研究重点是钻孔内部的自然对流,没有考虑含水层的渗透性。HOLMBERG等[17]针对GFBHE提出了一种简化的数值模型,发现钻孔内部的热阻在取热以及放热的不同运行模式下会有较大的不同。SPITLER等[9]得出了钻孔内部的努塞尔数与修正的瑞利数相关的结论。SHI等[19]对比研究了不同结构的换热器在地下水横流以及自然对流共同作用下的表现。JOHNSSON等[12]证明了钻孔内自然对流的存在,得出钻孔内水的有效导热系数的范围为1 ~ 2 W/(m∙K)。
如前所述,已经有许多研究者针对GFBHE利用热响应测试以及数值模拟的方法进行了研究。然而,热响应测试并不能真实地反映出GFBHE在冬季工况下的性能表现;已经建立的数值模型也大多数没有考虑发生在含水层内部的自然对流现象,而根据已有的研究,钻孔周围的对流现象会影响换热器的换热性能[20,21,22,23]。考虑到GFBHE的应用环境往往具有渗透性良好的条件,钻孔外部的对流传热也应考虑到模型之中。本文从GFBHE的实际应用环境出发,建立GFBHE的三维数值模型,并与灌浆填充的BHE（grout-grouted BHE, GGBHE）进行对比。研究两种BHE在相同环境和运行条件下的性能差异以及影响GFBHE性能的关键因素,分析在富水区域利用GFBHE取代传统灌浆BHE的可行性,让更多人认识到GFBHE的使用价值,为地热换热器的建设提供更多选择。
1 模型描述
1.1 模型假设
为了完整地描述GFBHE的传热过程,建立了耦合导热以及对流的三维非稳态数值模型。模型涉及U型管内外的对流传热以及钻孔与含水层之间的换热,现作出以下假设：
（1）含水层被视为均质且各向同性的饱和多孔介质;
（2）忽略地下水横流,仅考虑自然对流的作用;
（3）忽略地温梯度,认为含水层初始水温分布均匀。
上述假设对GGBHE模型同样适用。
1.2 数值模型
1.2.1 连续性方程
$$\frac{\partial \rho }{\partial t}+\nabla \cdot \left( \rho \vec{v} \right)=0 (1)$$
式中：$\rho $为流体的密度;$\vec{v}$为速度矢量。
1.2.2 动量守恒方程
流体的动量守恒方程如下：
$$\frac{\partial }{\partial t}\left( \rho \vec{v} \right)+\nabla \cdot \left( \rho \vec{v}\vec{v} \right)=-\nabla p+\nabla \cdot \bar{\bar{\tau }}+\rho \vec{g} (2)$$
式中：p为静压强;$\bar{\bar{\tau }}$为应力张量;$\rho \vec{g}$为浮升力项。
由于本文对自然对流的模拟采用Boussinesq假设,式（2）中除了浮力项之外其他项中的密度$\rho $为常数,对于浮力项中的密度$\rho $由式（3）近似表达。
$$\rho ={{\rho }_{0}}\left[ 1-\beta \left( T-{{T}_{0}} \right) \right] (3)$$
式中：${{\rho }_{0}}$为参考密度;${{T}_{0}}$为参考温度;β 为流体热膨胀系数;T为流体温度。
多孔介质区域的动量方程可以通过在式（2）的基础上加入动量源项的方式实现,由于本文不考虑地下水的横流对传热的影响,只考虑由BHE换热引起的自然对流,而由自然对流引起的水流速度不高,因此多孔介质中的水流可通过达西定律进行描述：
$$\nabla p=-\frac{\mu }{\alpha }\vec{v} (4)$$
式中：μ为流体黏度;α为渗透率。
1.2.3 能量方程
单一流体的不可压缩流能量方程有如下形式：
$$\frac{\partial }{\partial t}\left( \rho E \right)+\nabla \cdot \left[ \vec{v}\left( \rho E+p \right) \right]=\nabla \cdot \left( K\nabla T \right) (5)$$
式中：E为总能;$K$为导热系数;T为流体温度。
对于不可压缩流,式（5）中总能E由下式表达：
$$E=h+\frac{{{v}^{2}}}{2} (6)$$
固体区域能量方程形式如下：
$$\frac{\partial }{\partial t}\left( \rho h \right)=\nabla \cdot \left( k\nabla T \right) (7)$$
式（6）和式（7）中,h为比焓,其表达式为：
$$h=\int_{{{T}_{\text{ref}}}}^{T}{{{C}_{p}}}\text{d}T (8)$$
其中：Cp为相应介质的比热容;Tref为参考温度。
多孔介质中的传热过程采用热平衡假设模型[24]进行描述,即认为在传热过程中多孔介质固相以及液相的温度保持一致。多孔介质的能量方程如下：
$$\frac{\partial }{\partial t}\left[ \varepsilon {{\rho }_{\text{f}}}{{E}_{\text{f}}}+\left( 1-\varepsilon \right){{\rho }_{\text{s}}}{{h}_{\text{s}}} \right]+\nabla \cdot \left[ \vec{v}\left( {{\rho }_{\text{f}}}{{E}_{\text{f}}}+p \right) \right]=\nabla \cdot \left( {{k}_{\text{e}}}\nabla T \right) (9)$$
$${{K}_{\text{e}}}=\varepsilon \cdot {{K}_{\text{f}}}+\left( 1-\varepsilon \right)\cdot {{K}_{\text{s}}} (10)$$
式中：Ef为流体总能;hs为固体比焓;ρf为液体密度;ρs为固体密度;Kf为液体导热系数;Ks为固体导热系数;Ke为等效导热系数;ε为孔隙率。
2 基本模拟参数
本文U型管换热器的材料采用较为常见的高密度聚乙烯管。基于对比分析的出发点,两种BHE的基本结构参数相同。模型的整体直径为15 m,钻孔和含水层的长度均为80 m,直径为130 mm;U型管外径为25 mm,壁厚2.5 mm,管中心距60 mm。默认的含水层孔隙率为0.3,渗透系数为1 × 10-5 m/s。假设含水层的初始平均温度为20℃。换热器入口水流速度设置为0.5 m/s,入口水温为5℃。忽略模型上下面的传热,模型的侧面边界设置为恒温边界条件。其他参数如表1所示。
Table 1 Physical parameters used in the simulation表1 模拟用到的物性参数 |
| 材料 | 导热系数 / [W/(m2∙K)] | 密度 / (kg/m3) | 比热容 / [J/(kg∙K)] |
|---|---|---|---|
| HDPE | 0.44 | 950 | 1 500 |
| 含水层岩石 | 2.50 | 2 350 | 940 |
| 灌浆材料 | 0.93 | 1 800 | 840 |
| 水 | 0.60 | 999.7 | 4 182 |
3 结果与讨论
GFBHE换热过程复杂,涉及的影响因素较多,本文主要对比分析包括渗透系数、钻孔孔径、岩石导热系数在内的因素对GFBHE以及GGBHE性能的影响,并分析两种BHE的运行特性。模拟所涉及的含水层渗透系数的参数变化范围在1×10-5 ~ 1×10-3 m/s之间,含水层温度的变化范围在10 ~ 30℃之间,钻孔直径变化范围在110 ~ 190 mm之间。本文所进行的研究仅针对冬季工况展开。
3.1 两种BHE的运行特性分析
图2是GFBHE和GGBHE在默认条件下模拟120 h得到的出口温度变化曲线。从图中可以看到,在默认的条件下GFBHE的出口温度（Tout）始终高于GGBHE,这在一定程度上体现了GFBHE的可靠性能。在模拟的时间范围内,GFBHE和GGBHE的平均换热功率分别为3 270.8 W和2 915.8 W,GFBHE的换热性能要略优于GGBHE。此外,从图中还可以看出,两种BHE的出口温度变化趋势相同,在模拟进行的前20 h内出口温度都迅速降低,之后温度变化逐渐趋于平缓,120 h左右时,出口温度已经相对稳定。因此接下来将取模拟第120 h的结果来分析两种BHE的运行特性,并研究不同因素对两种BHE性能的影响。
Fig. 2 Water temperature changes with time at the outlet of two BHE图2 两种BHE出口水温随时间的变化 |
图3显示了GFBHE和GGBHE在渗透系数为1×10-3 m/s时运行120 h后钻孔周围不同深度位置含水层温度变化情况,其中h表示含水层深度位置,横坐标r表示距钻孔轴心的距离。可以看到,GGBHE钻孔周围不同深度的温度分布几乎相同,图中各点温度高度重合,而GFBHE钻孔周围不同深度温度分布则有极大的不同。显然,越靠近含水层底部,GFBHE钻孔周围的含水层温度变化越大,其热影响范围也越大,说明GFBHE的热提取过程更倾向在含水层底部进行。基于有关深井换热器的报道[25,26,27],推测两种BHE运行特性的不同是由于GFBHE钻孔内的水可以自由地在含水层与钻孔内进行交换造成的。钻孔内参与换热的冷流体在进入含水层之后受重力的影响向含水层底部流动,随着换热的不断进行,冷流体不断在底部聚集使含水层底部的热影响范围逐渐扩大。而GGBHE由于不存在GFBHE的这种特性难以在含水层建立起足够强度的自然对流,不能充分发挥含水层的自然对流换热作用。
Fig. 3 Temperature distribution at different depths of aquifer after 120 h operation of two BHE with permeability coefficient of 1×10-3 m/s图3 渗透系数为1×10-3 m/s时两种BHE运行120 h后含水层不同深度位置温度分布情况 |
为验证上述推测绘制了含水层和钻孔内的流线图与温度云图。从图4中可以看到,钻孔内的水自上而下流动,在靠近钻孔底部的位置进入含水层,而含水层内的水在靠近上方的位置进入钻孔,实现了钻孔内外的冷热流体交换;含水层只有在靠近钻孔下方的位置有比较明显的温度场变化。
Fig. 4 Streamline and temperature distribution in boreholes and aquifers图4 钻孔及含水层内部流线和温度分布图 |
3.2 含水层渗透系数对两种BHE性能的影响
不同渗透系数下,两种地热换热器平均功率的对比如图5显示。从图中可以看出,GFBHE的换热功率随着渗透系数的增大而增大,在渗透系数小于1×10-4 m/s时,渗透系数的改变对GFBHE换热性能影响相对较小,当渗透系数大于1×10-4 m/s时,渗透系数的改变对GFBHE的性能影响较大。此外,从图中还可以看到GGBHE换热性能受渗透率改变的影响并不明显。
Fig. 5 Changes of average heat transfer power with the permeability coefficient of two BHE图5 两种BHE的平均换热功率随渗透系数的变化 |
图6所示为两种BHE在不同渗透系数下进出口温差随时间的变化情况,图中GGBHE的进出口温差曲线几乎重合,而GFBHE的进出口温差曲线则明显受到渗透系数的影响,渗透系数越大,进出口温差越大。因此含水层渗透率越高的地区使用GFBHE得到的收益更大。含水层的渗透性越好,钻孔内外的原水交换就越容易,原水交换的同时,实现了对井内外换热效果的提升。而GGBHE没有钻孔内外的原水交换,对于含水层的渗透性变化的不敏感也从侧面说明了钻孔内外原水交换的重要性,因此在高渗透性含水层使用GGBHE可能无法充分利用含水层的渗透性来提升地热换热器的效果。当然,只有在渗透系数大于1×10-4 m/s时,GFBHE的换热效果提升才比较明显。
Fig. 6 Change of inlet and outlet temperature difference under different permeability coefficient图6 不同渗透系数下进出口温差变化情况 |
3.3 钻孔直径对两种BHE性能的影响
钻孔直径是BHE在建设过程中的一个重要参数,钻孔直径的选择不仅影响BHE的换热效果,也影响BHE的建设费用,因此本节对比分析两种BHE在不同孔径下的运行特性。在模拟过程中发现,GGBHE换热功率随孔径大小变化趋势与LUO等[28]的描述不符,分析认为这是由于灌浆材料导热能力的差异造成的。目前大多数灌浆材料的导热系数集中在0.8 ~ 2.4 W/(m∙K)之间[29],而导热系数的大小是灌浆材料的一个重要特性。因此,为了使对比结果更加可靠,人为地将灌浆的导热系数修改为1.5 W/(m∙K)和2.5 W/(m∙K)进行模拟,并将所有模拟结果绘制在图7中,GGBHE、GGBHE-1.5和GGBHE-2.5分别代表灌浆导热系数为0.93 W/(m∙K)、1.5 W/(m∙K)和2.5 W/(m∙K)的GGBHE。从图7中可以看到GFBHE的换热能力完全可以与GGBHE相比。而且随着孔径的增大,GFBHE的换热性能也相应提高。当孔径为190 mm时已经可以和灌浆导热系数最大的GGBHE性能非常接近。而GGBHE的性能随孔径的变化情况则随着灌浆导热系数的变化而有所不同,导热系数较小的GGBHE换热能力随孔径的增大而减小,而采用较大导热系数的GGBHE性能则随孔径的增大而有所提高。
Fig. 7 Changes of average heat transfer power with borehole diameter of two BHE图7 两种BHE的平均换热功率随孔径的变化情况 |
图8显示了GFBHE钻孔空间中水的平均流速随孔径的变化情况,可以看到钻孔内水的平均流速随着孔径增大而增大,而水流速度的提升有利于对流传热。对于GFBHE换热性能随孔径增大而提升的特性,认为是由于钻孔空间的增大使得钻孔内的自然对流作用增强引起的。因此,建议在建设GFBHE的过程中,在费用允许的前提下尽量使用较大的钻孔孔径。
Fig. 8 Change of average velocity with borehole diameter of water in GFBHE borehole图8 GFBHE钻孔内水的平均流速随孔径的变化 |
3.4 含水层温度对两种BHE性能的影响
GFBHE传热的关键在于自然对流传热过程,而自然对流由温差进行驱动,可以预测GFBHE的换热性能受含水层温度的影响较大。图9显示了两种BHE在不同的地温条件下的换热功率以及出口温度变化情况。可以看到两种BHE的换热功率以及出口温度都随着含水层温度的升高而升高。但两种BHE换热性能的增长速度并不相同,GFBHE的换热功率增长曲线斜率更大。在含水层温度为10℃时两种BHE的换热性能不相上下,随着含水层温度的升高,GFBHE相对GGBHE的性能提升越来越明显。因此温度越高的含水层,使用GFBHE可以获得更好的换热性能。
Fig. 9 Heat transfer power and outlet temperature of two BHE under different ground temperature conditions图9 两种BHE在不同地温条件下的换热功率及出口温度 |
4 结论
针对GFBHE和GGBHE建立了三维非稳态数值模型并进行了求解,对比分析了GFBHE和GGBHE的性能差异,研究了包括渗透系数、钻孔孔径、含水层温度在内的关键因素对两种BHE性能的影响。主要结论如下：
（1）含水层渗透性对GFBHE的影响较大,对GGBHE的影响不明显。当渗透系数大于1×10-4 m/s时,随着渗透系数的增大,GFBHE的换热性能明显提升,当渗透系数小于此值时,渗透系数的变化对GFBHE的影响相对较小。
（2）在冬季工况下,随着含水层温度的升高,GFBHE和GGBHE的换热功率都在提升。但是,GFBHE换热功率变化曲线斜率比GGBHE的要大,说明随着地温的升高GFBHE性能提升更大。
（3）GFBHE的性能随着钻孔孔径的增大而升高。GGBHE换热性能随孔径变化情况则与灌浆的导热能力有关,采用导热能力小的灌浆的GGBHE性能随孔径增大而下降,反之,则提升。
总之,在富水条件下建设GFBHE取代传统的灌浆BHE是可行的,并且由于节省了灌浆的费用,GFBHE的经济性可能会更好。此外,从各因素对两种BHE的影响来看,越是高渗透性以及高地温的地区使用GFBHE效果越好,建造GFBHE时也应该有意使用较大的孔径以获得更好的换热效果。
