.地下水源热泵取水引起含水层土体颗粒运动的力学特性分析.docxVIP

地下水源热泵取水引起含水层土体颗粒运动的力学特性分析来源：全国地源热泵网　日期：2011-07-27　点击：次王松庆 张旭（同济大学 暖通空调及燃气研究所，上海 200092）摘? 要：本文对地下水源热泵地下取水时引起含水层土体颗粒运动的力学特性进行了理论分析，为研究和解决地下水源热泵取水井长期运行后可能出现的井壁坍塌和井壁周围颗粒重组等实际问题提供了理论依据。关键词：地下水源热泵；含水层；颗粒运动1 前言地下水源热泵是地源热泵的一个分支，它把地下水作为热泵的低位热源。由于较深地层中的地下水不会受到大气温度变化的影响，故能常年保持较为恒定的温度，既高于冬季的室外空气温度，也低于夏季的室外空气温度，且具有较大的热容量，因此地下水源热泵系统的效率比空气源热泵高，且不存在结霜等问题。所以近几年地下水源热泵在国内如山东、河南、湖北、辽宁、黑龙江、北京和河北等地得到了广泛的应用[1]。但是地下水源热泵在使用的过程中也相继出现了一些问题。如文献[2]中指出由于地下水源热泵的取水井常常连续取水，而连续取水所引起的最显著问题是井水含沙问题，如把含沙的地下水长期不断的抽到空调管路系统中，可能造成管路阻塞及损坏，更有甚者因长期抽水而形成取水井的井壁坍塌，造成地下水源热泵系统无法继续使用的严重后果。文献[3]也指出当注水井又兼抽水井时，反复的抽、注水会引起井壁周围的细颗粒介质的重组。这些情况都说明在研究地下水源热泵时，研究地下取水对含水层土体颗粒所产生的影响是十分必要的。本文通过地下水源热泵地下取水时引起含水层土体颗粒运动的力学特性进行分析，为研究和解决地下水源热泵取水井长期运行后可能出现的井壁坍塌和井壁周围颗粒重组等实际问题提供了理论依据。2 土体颗粒受力状况分析2.1 土体颗粒受力状况分析在取水井附近取水，土体颗粒同时受到水流的拖曳力，粘着力，水的升力以及自身的重力等多种作用力。2.2 拖曳力对于半径为r的球形颗粒，在无限大来流速度为 的液体中的受到的拖曳力为 ，可以表示为[4]：Fd=Cd(1/2)ρfv2Sp? （1）式中：ρf ——水的密度；Sp——颗粒的迎风面积，当其为完整的圆形时，Sp=πr2；Cd——阻力系数。对于阻力系数，不同的学者给出了不同的计算公式。Stokes认为，在流速和颗粒雷诺数很低的情况下，可忽略N-S方程式中的惯性项，阻力系数可表示为：Cd=(24/Rep) (Rep＜1)? （2）Ossen近似的考虑了惯性项而得到的阻力系数计算公式为：Cd=(24/Rep) [1+(3/16) Rep](Rep＜5) （3）对于颗粒雷诺数较大的情况，也可以采用下面经验公式进行计算。Cd=(24/Rep) [1+0.15Rep0.687] （4）同时，土体颗粒受到的拖曳力与壁面切应力有关，其计算式也可以表示为[5]：FH=32 r2τ? （5）当时，临界拖曳力为FH=32 r2τcr式中：τ——切应力；τcr——临界切应力。对处于流速较慢的层流中的颗粒来说，其拖曳力计算公式还可以表示为[6]：FH=1.7(6π)μrU（6）式中： U——与壁面距离为r处的水流速度；μ——水的粘滞系数。2.3 粘着力粘着力是研究颗粒从固体骨架脱离过程中的一个重要的组成部分。Sharma认为粘着力主要是作用在颗粒上的van der Waals力，与平衡分离距离有关[5]。Ingles和Young则认为粘着力还与静电力（其大小与土壤、水流中的离子特性有关）、水流的毛细管作用等因素有关[7-8]。对于半径为r的可变形球形颗粒，当它与光滑表面接触时，其受到的van der Waals力可通过下式计算[9]：Fv=(Ar/6h2)[1+(a2/hr)]（7）式中： Fv——van der Waal力；A——Hamaker 常数；h——分离距离；a——接触半径。我们可以近似的认为粘着力等于van der Waals力。因此，粘着力可以表示为：FA≈Fv（8）2.4升力颗粒所受到的升力是指使得颗粒脱离表面的垂直方向的力。根据文献[5]和[10]可知，升力的大小与粘着力相比，相差几个数量级，因此通常在考虑了粘着力的情况下，可以不考虑升力对颗粒的影响。同时文献[10]还指出，颗粒受到的重力小于其受到的升力，因此在一般条件下，我们可以不考虑颗粒的自重。3 土体颗粒运动状况分析土体颗粒从含水层固体骨架脱离这一过程受到了很多因素的影响。如受到颗粒自身性质的影响（几何形状、尺寸大小等）；水流的物理和化学性质的影响（如流速、温度、pH值等）[11]。本文不考虑化学作用对颗粒脱离过程的影响，仅从物理方面（力学角度）对脱离过程进行分析。本文假设颗粒仅受到由水流引起的拖曳力和粘着力，忽略升力和颗粒自身重力及其它外力对颗粒脱离的影响。文献[5]中指出颗粒从固体骨架的脱离可能是通过滑动的