第一作者简介: 张乐(1992—),女,助理工程师,主要从事水工环地质工作与研究。Email: 765213301@qq.com。
通信作者简介: 周阳(1986—),男,工程师,主要从事水工环地质工作与研究。Email: 77196410@qq.com。
近年来,浅层地热能受到广泛关注,但由于在开发利用过程中对地层导热性能分析不足,现有开发利用项目失败率较高。为解决这一问题,对陕西省不同地貌单元已开展的46组热响应试验结果进行整理,以线热源模型为基础,采用叠加法对陕西省典型地貌单元的导热系数进行了分析统计。结果表明: 山间盆地、沙漠高原和断陷盆地3种地貌单元的导热能力依次减小; 在同一地貌单元中,第四系黄土覆盖层越厚,钻孔的导热能力越差; 黄土塬地区导热系数大小差距较大,可进一步深入研究该地貌单元内导热系数的影响因素。通过对不同地貌单元岩土导热系数分布特征的研究,提出了各地区浅层地热能开发利用的相关建议,从而为进一步高效、高质开发浅层地热能提供思路和技术支持。
In recent years, shallow geothermal energy has received extensive attention, but the failure rate of existing projects is high due to insufficient analysis of the thermal conductivity of the strata during the development and utilization process. In order to solve this problem, the authors collated 46 sets of thermal response test results that have been carried out in different geomorphological units in Shaanxi Province. Based on the linear heat source model, the thermal conductivity of typical geomorphological units in Shaanxi Province was analyzed by superposition method. The results show that the thermal conductivity of three geomorphological units(mountain basin, desert plateau and fault basin) decreases gradually. In the same geomorphological unit, the thicker the loess cover of Quaternary, the worse the thermal conductivity of the borehole. The difference of thermal conductivity in loess tableland is quite large, and the influencing factors of thermal conductivity in the geomorphological unit can be further studied. Based on the study of the distribution characteristics of thermal conductivity in different geomorphological units, some suggestions for the development and utilization of shallow geothermal energy in different regions were put forward, thus providing some ideas and technical supports for further efficient and high-quality development of shallow geothermal energy.
浅层地热能资源(Shallow Geothermal Resources)是指蕴藏在浅层岩土体、地下水、地表水中的地热资源[1], 一般低于25 ℃[2]。最新数据表明, 我国287个地级以上城市浅层地热资源量为每年2.87× 1020 J, 相当于95亿t标准煤, 至少可减少二氧化碳排放量6.52亿t[3]。据调查, 陕西省共有浅层地热能开发利用工程234处, 运行效果良好的仅占39.7%。对岩土体导热性能分析的不足, 是导致浅层地热资源开发工程失败率高的主要原因。
为设计浅层地热能开发利用工程, 要对岩土导热系数进行测试分析。对于岩土热物性参数的测量, 一般采用热响应试验法, 这种方法常用的数学模型有线热源模型和柱热源模型。通过热响应试验可较精确地得到岩土的热物性参数, 但这种方法的时间成本较高。为解决时间成本与工程成功率之间的矛盾, 本文基于大量热响应实验数据, 采用叠加法进行分析整理, 发现岩土导热系数与其所在地貌单元存在一定的相关性, 初步得到了适宜修建浅层地热能开发利用工程的地貌条件, 为浅层地热能的开发、利用提供了适宜性区划。
陕西省地貌总体呈SN走向, 北部为高原, 中部为盆地, 南部为山区, 地貌结构复杂, 由北向南大体可分为4种地貌类型, 分别为陕北沙漠高原、陕北黄土高原、关中断陷盆地和陕南山间盆地(图1)。
陕西省北部地区属大陆性季风气候, 年平均气温8~12 ℃[4], 年降水量300~500 mm[5]; 中部地区属温带湿润、半湿润气候, 年平均气温为11~13 ℃[6], 年降水量500~700 mm; 南部地区属亚热带湿润、半湿润气候, 年平均气温达13 ℃以上, 年降水量大于700 mm。
陕北沙漠高原区主要分布于长城以北, 是毛乌素沙漠的南缘。构造单元属鄂尔多斯盆地东翼北部的陕北单斜[7], 地势平坦, 以活动沙丘、沙垄及片沙为主。沙丘、沙地之间湖泊、海子众多, 间有滩地。地势总体西高东低, 海拔900~1 400 m。沟壑不发育, 风蚀、风积地貌发育。地下水埋深一般为10~30 m, 个别地方埋深低于5 m, 渗透系数多小于5 m/d。
陕北黄土高原是在新近纪末起伏和缓的准平原基础上历经第四纪以来多次黄土堆积和侵蚀作用形成的。地貌单元以延安为界: 北部地形破碎, 沟壑发育, 水土流失严重, 生态环境脆弱; 南部中低山夹黄土塬。黄土梁峁地区地形破碎, 不利于地下水的存储, 导致水量贫乏; 黄土塬地区的降水存储条件较好, 洼地富水更为丰富。
关中断陷是新生代断陷盆地, 盆地南依秦岭, 北连黄土高原, 西狭东阔, 盆地两侧地形向渭河倾斜, 呈阶梯状地貌景观。盆地中共发育13条张性或张扭性断裂, 走向呈EW向或NW向, 倾角介于50° ~80° 之间, 控盆断裂6条, 控单元断裂7条。水文地质方面, 该地区具有第四系松散岩类典型盆地的水文地质特征[8]。
陕南山间盆地区是经断陷作用和堆积作用形成的、由宽阔阶地、坝子、丘陵、河谷等构成的地貌单元, 区内主要有汉中盆地、安康盆地、商丹盆地、洛南盆地和石门盆地。山间盆地的含水岩性相互交接, 并被各类地质构造穿插。由于地形、地貌、岩性和构造的不同, 该地区富水性分布极不均匀。汉中盆地和安康盆地富水性极好, 商丹盆地富水性中等。总体来说, 地势越高, 富水性越差。
欧美许多国家在20世纪中叶就已经开始了对地埋管地源热泵换热过程的研究[9], Ingersoll和Plass于1948年利用Kelvin线热源理论首次提出了线热源模型, 并将其应用到无限大土壤中[10, 11, 12]。时至今日, 线热源模型仍广泛应用于地埋管地源热泵设计、浅层地热能分析、岩土导热系数计算等多个方面。
线热源模型应用于地埋管换热器的换热过程中时, 需做如下简化假定[12, 13]:
(1)地埋管周围土体无限大, 且具有均匀性和各向同性, 初始温度和热物性参数相同且不随温度的变化而变化;
(2)只考虑径向导热, 忽略轴向导热;
(3)不考虑钻孔的几何尺寸, 近似将其视为轴心上的线热源。
载热流体平均温度
式中:
经整理可得
式中:
热响应试验中, 流体平均温度
式中:
经整理可得
由上式可知如何绘制
常规的热响应试验分析方法是将2次热响应试验看作2次孤立的热负荷, 这种计算方法认为第2次热响应试验开始时的地温条件与第1次热响应试验开始时的地温条件完全一致, 但其实地温场已经受到干扰, 很可能2次结果的初始条件并不完全相同, 这就导致2次热响应试验计算出的导热系数结果相去甚远。因此, 高宽等[14]及周世玉[15]提出了线热源的叠加法。叠加法认为岩土的导热系数、热阻、扩散率等都是岩土的固有属性(不随温度的变化而变化), 在此基础上将2次热响应试验等效成3个持续热负荷叠加的结果, 从而避免了地温恢复这一未知过程所带来的误差。
热响应试验过程为: 0至
第1次热响应试验稳定后流体平均温度为
第2次热响应试验稳定后流体平均温度为
式中:
第1次热响应试验稳定后流体平均温度为
第2次热响应试验稳定后流体平均温度为
由上式可知, 基于叠加法的原理, 在连续2次热响应试验中, 第1次热响应试验稳定后流体平均温度
由此, 通过叠加法可计算出岩土体的导热系数。
本文以西安市一处热响应试验点为例, 分别进行了线热源法和叠加法计算, 并将2者结果做了对比。该点位于黄土台塬, 底层为第四系全新统冲积物, 岩性为粉质黏土、砂层, 钻孔深120 m, 地下水位埋深21.6 m。
2.3.1 线热源法
分别将第1次和第2次热响应试验过程中流体平均温度t与试验历时的对数
2.3.2 叠加法
利用第1次热响应试验数据绘制
可见, 利用叠加法计算导热系数可有效减少因2次热响应试验之间的相互干扰而对导热系数计算结果产生的影响, 能更准确地模拟热响应试验过程中的实际情况, 因此通过叠加法计算出的导热系数也更加精确可信。
本文所采用数据来源于陕西省水工环地质调查中心2013— 2018年多个浅层地热能相关项目, 包括“ 陕西省大中型城市浅层地热能调查评价” “ 陕西省主要城市浅层地温能开发区1∶ 5万水文地质调查” “ 陕西省浅层地热能示范基地建设” 和“ 陕西省浅层地热能赋存规路及开发利用关键技术” 等。项目组在省内4大典型地貌单元中选取具有代表性的地段, 开展了46组热响应试验, 试验数据比较丰富。
本文利用叠加法对已有的46组热响应试验数据(其中断陷盆地28组, 山间盆地9组, 黄土高原4组, 沙漠高原5组)进行分析计算, 并对计算结果进行处理, 筛选出2次热响应试验计算所得导热系数相对误差小于10%的试验点共计39组。后将这39组综合导热系数计算结果分类绘制于图中(图4)。
由图4可知: 总体来说, 山间盆地、沙漠高原、断陷盆地3种地貌单元, 岩土体综合导热系数依次减小; 黄土高原地貌单元中, 不同地貌类型导热系数差距较大。总之, 各地貌单元导热系数具有不同特征。
断陷盆地地貌单元包括低级阶地、高级阶地、黄土塬3种类型。其中: 低级阶地地区导热系数值分布最为集中, 集中分布于1.49~1.74 W/(m· K)之间, 均值为1.62 W/(m· K); 高级阶地地区导热系数值分布于1.09~1.83 W/(m· K)之间, 均值为1.44 W/(m· K), 各点导热系数值几乎均匀地分布于均值周围; 黄土塬地区导热系数值分布于1.34~2.30 W/(m· K)之间, 均值为1.62 W/(m· K), 偏向较小的边缘, 存在某个点计算所得导热系数过大, 从而拉高了平均水平的现象。就均值而言, 高级阶地< 低级阶地=黄土塬; 就区间长度而言, 黄土塬> 高级阶地> 低级阶地。这说明黄土塬地区热响应试验钻孔的导热能力分布并不均匀, 可能受局部地貌或地质条件影响较大。
山间盆地地貌单元包括低级阶地、高级阶地2种类型。低级阶地地区和高级阶地地区导热系数值分布都较为集中: 低级阶地地区导热系数值分布于2.04~2.36 W/(m· K)之间, 均值为2.19 W/(m· K); 高级阶地地区导热系数值分布于1.93~2.37 W/(m· K)之间, 均值为2.11 W/(m· K)。就区间长度上和均值而言, 两者均相差不大, 高级阶地区间长度略长于低级阶地地区, 高级阶地地区均值略低于低级阶地地区。
本次热响应试验涉及黄土高原地貌单元中的河谷阶地地区和黄土塬地区。河谷阶地地区导热系数值分布于2.01~2.32 W/(m· K)之间, 均值为2.16 W/(m· K); 黄土塬地区导热系数值分布于0.81~1.84 W/(m· K)之间, 均值为1.33 W/(m· K)。两者对比, 黄土塬地区导热系数值的不稳定性高于河谷阶地地区, 黄土塬地区导热系数均值明显低于河谷阶地地区。
沙漠高原地貌单元包括风沙河谷和风沙滩地2种类型。风沙河谷地区导热系数值分布于2.02~2.17 W/(m· K)之间, 均值为2.10 W/(m· K); 风沙滩地地区导热系数值分布于1.74~1.99 W/(m· K)之间, 均值为1.83 W/(m· K)。两者导热系数区间长度均较小, 风沙滩地地区导热系数略高于风沙河谷地区。
总的来说, 山间盆地、沙漠高原、断陷盆地3种地貌单元的导热能力依次减小; 同一地貌单元中, 第四系黄土覆盖层越厚, 钻孔导热能力越差; 黄土塬地区导热系数大小差距较大。
研究表明, 陕西省岩土体导热系数与其所处地貌类型有一定关系。山间盆地、沙漠高原、断陷盆地3种地貌单元的导热能力依次减小。在同一地貌单元中, 第四系黄土覆盖层越厚, 钻孔的导热能力越差。对断陷盆地低级阶地地区、山间盆地、黄土高原河谷阶地地区和沙漠高原而言, 地貌因素对岩土导热能力影响较大, 在相同地貌条件下, 岩土体导热系数相近。
本研究对陕西省浅层地热能开发、利用及选址具有指导意义。陕西省山间盆地、沙漠高原风沙河谷地区以及黄土高原河谷阶地地区岩土体导热性能良好, 且受地域影响不大, 可进行浅层地热能的推广利用; 断陷盆地低级阶地地区和沙漠高原风沙滩地地区的岩土导热性能受地域影响不大, 但导热能力一般, 可适当开发、利用浅层地热能资源; 在黄土高原黄土塬地区、断陷盆地高级阶地和黄土塬地区, 不同地域岩土导热系数相差较大, 应因地制宜, 针对不同地区、不同地质条件测定岩土热物性参数, 分析其是否适宜进行浅层地热能开发, 同时进一步有针对性地探讨影响该地区岩土导热系数的主要因素, 从而指导在该地貌单元上进行浅层地热能的开发、利用。
致谢: 参与本次野外调查工作的同事还有赵智强、宋巧莲、朱红玉、张卉、王克、桂忠强、王鸽、侯娟等, 在此表示衷心感谢!
(责任编辑: 刘丹)
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|