水源热泵工程集中区地下水温度场的数值模拟分析

《水源热泵工程集中区地下水温度场的数值模拟分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、水源热泵工程集中区地下水温度场的数值模拟分析摘要：为了避免地下水源热泵用户集中区用户间的热贯通，需要合理布设取回水井。文章采用有限单元法建立地下水温和水位数值模型，预测水温和水位的变化趋势，模拟结果区内地下水温冬季变化范围6.3～14.2℃，夏季变化范围11.5～21.2℃，满足水源热泵工程对取回水的要求，能够避免了用户间的热贯通。通过建立合理的数值模拟，能够定量地分析水源热泵工程抽回水对地下水温度场变化，为合理取回水井设计提供依据。关键词：水源热泵；地下水位和水温；数值模拟；热贯通地下水源热泵空调系统是

2、一种利用浅层地热能进行制冷和供暖的采能技术[1],近年来该系统的工程应用已比较广泛［2］。在应用过程中，盲目的利用水源热泵技术，在新建水源热泵工程设计过程中有时只考虑满足自身的使用，而忽略了对周围地下水环境的影响，从而造成热贯通[3-4]等现象，导致水源热泵机组换热效率降低。笔者结合实际项目建立地下水水温和水位的数值模型，预测水温和水位的变化趋势，为在水源热泵用户较多的地区的新建水源热泵工程项目提供参考依据。沈阳市某区A项目拟采用水源热泵工程，其采暖和制冷面积为59920.96m2。该项目采暖期用水量为8

3、50432m3，制冷期用水量为678442m3，设计抽水井5眼，回水井17眼。项目周围水源热泵用户比较多，如图1所示。B用户距离A100m，用水量26.48m3/h，有1眼抽水井和1眼回灌井；C用户距离A260.0m，用水量90.0m3／h，1眼抽水井，3眼回水井；D用户距离A区380.0m，用水量为253.0m3/h，2眼抽水井，6眼回灌井；E用户距离A100.0m，用水量为98.0m3/h，1眼抽水井，1眼回灌井；F用户距离A230.0m，用水量为240.0m3/h，抽水井2眼，回灌井7眼，上述用户均

4、利用水源热泵采暖和制冷。因此对地下水水温和水位进行数值模拟检验对已有用户B-F是否造成热贯通和取回水影响，成了该项目能否通过水行政专管部门审批的关键。各项目位置如图1所示。图1位置示意图1模型建立1.1水文地质概化模型①含水层类型的确定：含水层岩性为上部为中粗砂含卵砾石含有少量粘土，厚约21.5～22.5m，下部是由砂砾组成含粉质粘土，成半胶结状，砂砾石微风化，厚度约12.5m。将含水层视为潜水含水层，第三系风化砂砾岩被视为含水层底板。根据含水层的类型、岩性、厚度、导水特征及导温特征等，将模型概化为非均质

5、各向同性含水层，局部可视为均质。4②地下水特征.地下水位埋深为11.0～12.0m左右，地下水温度12℃.在天然状况下，场区内地下水流向为北东至西南,水力坡度为0.001-0.002左右.③地下水流的概化：计算区地下水位受枯丰水期影响有一定变化，水流呈非稳定状态，但总体上区域地下水为层流运动，水流符合达西定律，全区可视为非稳定二维平面流。④模拟区边界类型的划分：本次将模拟区4个方向侧面边界概化为定水位边界，其中东部边界水位为36.0m，西部边界水位为30.0m；模拟区地表入渗条件较差，大气降水入渗补给忽略

6、不计，将顶部边界概化为隔水边界；将下更新统地层定为隔水底板边界。源汇项的处理工程A场区内用户水源热泵工程共计5眼抽水井，17眼回灌井；工程B-F场区内用户水源热泵工程共计7眼抽水井，18眼回灌井。1.2地下水流与温度耦合数值模型（1）地下水流数值模型根据水文地质概念模型，建立研究区地下水水流数学模型，其具体表达式为（1）（2）（3）式中:n为孔隙度；ρ0为液体的参考密度，kg/m³；βP为水的压缩系数，Pa-1；p为地下水压力，Pa；t为时间，s；βT为水的热膨胀系数，℃-1；T为水和孔隙介质的温度，℃；

7、ρ为地下水密度，kg/m³；αb为空隙介质的压缩系数，Pa-1；▽为梯度算子；KP为渗透能力张量，㎡；μ为水的黏度，Pa·s；g为自由落体加速度，m2/s；D为计算区范围；Г1为第一类边界条件。（2）地下水温数值模型根据水文地质条件和热力学条件,建立地下水流与温度耦合热传递数值模型，如下：（4）（5）（6）式中:cƒ为水的比热容，J/（kg·℃）；ρS为空隙介质密度，kg/m³；cS为空隙介质比热容，J/（Kg·℃）；Kƒ为水的导热系数，W/(m·℃）；KS为空隙介质导热系数，W/(m·℃）；I为单位矩阵

8、；D为热机械弥散张量，W/(m·℃）；ν为水的流速，m/s。1.3区域剖分采用有限单元法进行地下水开采后水流数值模拟，模拟软件选用Feflow5.0，根据区内流场特征、系统的结构特征、水力特征和边界条件，将模拟区划分为的的长方形区域，模拟区的面积1.0km2。将模拟区剖分为6160个三角单元，含6248个节点。1.4水文地质参数确定计算所用的渗透系数、给水度、孔隙度等参数主要是根据试验和观测资料所确定，见表41。地下水密度、热