水文地质模拟分析及在地下水源热泵中的应用.doc

沈阳城区水文地质模拟分析及在地下水源热泵中的应用 摘要: 用数字模拟方法分析了沈阳城区地下水系统的水文地质条件。基于分析结果结合解析方法提出了沈阳城区地下水源热泵抽水井和回灌井的取水条件以及抽水井之间、回灌井之间、抽水井与回灌井之间的合理间距。对于指导沈阳城区地下水源热泵项目的建设具有直接的指导作用。 关键词: 沈阳城区; 地下水水源热泵; 地下水; 影响半径 中图分类号: P641. 8?文献标识码: A 文章编号: 1007-2284(2011) 04-0014- 04 0 引言 地下水源热泵空调系统是一种利用浅层地热能进行制冷和供暖的采能技术，最早出现于美国20 世纪40 年代，由于其环境效益显著，运行效率较高，在欧洲及美国迅速发展[1- 3]。我国于20 世纪末引入该技术，在全国各地得到了推广应用。 水源热泵理论上可以利用一切的水资源，但在实际工程中，不同的水资源利用的成本差异是相当大的。所以在不同的地区是否有合适的水源成为水源热泵应用是限制水源热泵技术的一个重要限制，而且还受含水层结构的限制。水源热泵若利用地下水，必须考虑水源的回灌，对于回灌技术，必须结合当地的地质情况，确保可以在经济条件下打井找到合适的水源，保证地下水源热泵的利用效率，才能使地下水源热泵得到广泛的推广应用。 但是在地下水源热泵的相关研究中，国内外文献多见于单井地下水源热泵对地下水的影响变化研究，对大区域大规模多井群地下水源热泵应用研究目前很少，缺少理论指导[4- 8]。本文在考虑地下含水层水及其地质条件影响基础上针对沈阳城区地下水源热泵井间距问题统进行研究探索，为地下水源热泵可持续健康发展提供参考[3- 6]。 1 研究区域概况 研究区沈阳处于东经123°16′- 123°35′、北纬41°42′-41°55′之间。区域总面积为468 km2 。地处新华夏第二隆起带与第二沉降带之交界地带，是辽宁省东部山地和中部平原的过渡带，断裂带比较多，基岩出露面积较小，由东北走向西南，但地壳较稳定。其地表为第四纪冲积层所覆盖，黏土、沙砾、砾石交互成层，堆积很厚。地层由下更新统、上更新统及全新统冲洪积堆积物构成。地貌形态包括构造剥蚀丘陵，山前冲洪积扇和冲积平原。沈阳地区以平原为主，地势平坦。北部部分地区为辽北丘陵，地势由东北向西南缓缓倾斜，西部是辽河、浑河冲积平原，东部为低山丘陵，属于辽宁东部山区的一部分。 据统计，截至2007 年底，沈阳市地源热泵技术应用面积已达1 840. 41 万m2 ，已建地源热泵项目中地下水地源热泵项目占到了总项目数的90%，占总供热面积的98. 4%。沈阳城区2007 年可利用地下水资源总量为17. 33 亿m3 ，能充分满足地下水源热泵的发展需求。 2 水文地质数学模型的建立 2. 1?沈阳城区水文地质概念模型 对沈阳城区地下水系统及环境分析，可将其地下水系统的概念模型做如下归纳: 区域内虽然有全全新冲洪积砂砾石空隙潜压系统、上更新统冲洪积砂砾石空隙潜水和微承压水压系统及其中下更新统冰水沉积、局部承压水压系统。但考虑3 个亚系统之间并无连续稳定的隔水系统，相互间具有密切的水力联系，且区域内地下水开采几乎皆为统一开采，故可将其视为一个统一的潜水含水系统。系统内含水介质其透水性为空间坐标有关，为非均质; 但任一点各方向的渗透性能相同，二者综合可将其含水层归为非均质各项同性含水层。 2. 2?沈阳城区地下水系统数学模型 通过分析沈阳市的地质条件和水文地质条件，含水层厚度与水文地质参数随空间变化而变化，但地下水流水平运动流速缓慢，符合达西定律。结合以上结构的概化，将区域内的地下水流系统概念化为非均质各项同性、二维平面、非稳定地下水流系统，建立如下数学模型: 式中: K 为含水层渗透系数，m/ d; Ω为计算区域; h 为水位，m; t为时间，d;m 为含水层厚度，m; S 为含水层的给水度; W(x ，y ，t) 为垂向补排强度，m/ d; Q(x j ，y j ，t) 为t 时第j 号井抽水量，m3/d; h0 初始时刻水位值，m; h1 (x ，y ，t) 为第一轮边界水位值，m; q(x ，y ，t) 第二类边界的单位长度流量，m3/d; Γ1 第一类边界; Γ2 第二类边界;边界的法线方向。 3 参数分区与水流模拟模型验证 ? 根据水文地质概念模型，并结合本水文地质资料，利用GMS 软件中的TINS 模块和solid 模块对研究区进行三维网格剖分，共剖分三角网4 391 个，形成数据节点2 280 个，每个三角网的面积为0. 10 km2 。根据工程勘察资料，在平面上将研究区域划分为8 个参数分区(图1、表1) 。 图1 沈阳城区渗透系数及给水度分区图 Fig. 1 The sectorized picture of osm