阮传侠，冯树友，牟双喜，程万庆，赵苏民

(1.中国地质大学(北京)水资源与环境学院，北京 100083； 2.天津地热勘查开发设计院，天津 300250)

天津地区地层热物性特征及影响因素分析

阮传侠1,2，冯树友2，牟双喜2，程万庆2，赵苏民2

(1.中国地质大学(北京)水资源与环境学院，北京 100083； 2.天津地热勘查开发设计院，天津 300250)

为分析天津地区地层热物性参数的地区分布特征和平面分布规律，布置88个勘查孔、现场采集1 076个岩土样、室内分析热导率、比热容、热扩散系数等地层热物性实验，对测试数据进行统计分析。结果表明：天津市岩土体热导率在1.26～1.70 W/(m·K)，比热容在2 050～2 090 J/(kg·K)，热扩散系数在0.45×10-6～0.74×10-6m2/s。同一地区不同岩性的比热容，黏土最大、粉砂最小，热导率刚好相反，热扩散系数与比热容规律相同；不同地区同一岩性的热物性参数差别不大。天津市比热容与热导率呈现大致相反的趋势，热导率高值区位于蓟县、宝坻和宁河的东部、武清西部、静海南部以及滨海新区的中部地区；比热容高值区位于蓟县、宝坻和宁河西部、武清东部、静海北部以及滨海新区大部分地区。

地热地质；热导率；比热容；热扩散系数；水文地质；地层热物性；天津地区

E-mail:rchx12@sohu.com

天津地区夏季炎热、冬季寒冷，气候条件决定了天津夏季需要制冷，冬季需要供暖。优越的水文地质、基础地质条件为开发利用地源热泵提供了基础条件，因此地源热泵在天津得到了快速发展。它是可再生能源开发利用的新方式[1～3]。地源热泵系统设计以地层热物性参数为基础[4～5]，如果热物性参数不准确，则可能导致地源热泵系统设计不当[6]，给后期的应用带来诸多问题，因此地层热物性参数对于地源系统工程设计非常重要。

2010年天津第一次开展了天津市浅层地热能调查工作，取得了一批天津市岩土体的热物性参数，为天津市开发利用浅层地热能提供了基础数据。孙宝成等[7]进行了天津市浅层地热能资源评价及开发利用条件分析的研究，主要开展浅层地热能适宜性分区和资源评价等。于建水等[8]开展了天津市浅层地热能开发利用经济和环境效益评价工作，为后期浅层地热能开发利用提供了基础依据。唐永香等[9～10]对天津地区浅层地热能赋存条件及前景进行了分析，并针对天津滨海新区浅层地热能资源评价及开发利用对策进行了研究。刘九龙等[11]开展了专项地埋管地源热泵系统工程的高效利用技术研究，主要从热泵利用效率方面开展了研究。曾梅香等[15]从地质环境监测系统设计和施工要求、监测数据采集传输及整理分析、监测系统运行与维护等方面进行研究，为浅层地热能动态监测提供依据。

据调查统计，截至2016年底天津市浅层地热能开发利用工程达到306个，应用建筑面积为933×104m2。其中地埋管地源热泵系统多达231处，可见地埋管地源热泵工程利用越来越广泛，对岩土体热物性特征的研究显得尤其重要。

为了保证浅层地热能资源高效、合理、可持续开发利用，在以往工作成果的基础上，对天津市不同区域相同地层和同一区域不同地层的热物性参数进行对比分析，研究热物性参数平面分布特征。结果表明天津市岩土体热导率为1.26～1.70 W/(m·K)，比热容为2 050～2 090 J/(kg·K)，热扩散系数为0.45～0.74×10-6m2/s。同一地区不同岩性的比热容，黏土最大、粉砂最小，热导率刚好相反，热扩散系数与比热容规律相同；不同地区同一岩性的热物性参数差别不大，为浅层地热能工程设计提供了基础参数。

1 地热地质条件

天津市地处华北平原东北部，环渤海的中心，东南临渤海，北依燕山山脉，南北距离为189 km，东西距离为117 km。天津市国土面积11 946 km2，北部山区面积755 km2，南部面积11 191 km2。天津市地貌主要有山地、丘陵、平原、洼地、滩涂等，地貌特征是从蓟县北部山区向南部平原逐级下降，西面从武清区永定河冲积扇尾部向东缓缓倾斜，南面从静海县南运河大堤向海河河口渐渐降低，形成北高南低、西高东低的形态。最高处为蓟县九山顶，海拔1 085.5 m，最低处为塘沽区大沽口，海拔为-2 m。

1.1地质条件

天津市地质构造复杂，新构造运动控制了新生代地层的沉积。

以宁河—宝坻断裂为界，北部为阴山纬向构造体系，南部为以新华夏构造体系，为一系列北北东向的隆起坳陷和构造断裂，自东向西的Ⅲ级构造单元依次有黄骅坳陷、沧县隆起、冀中坳陷和一系列北东向断裂。北部山区受二级燕山沉降带的控制，在蓟县山前断裂与宁河—宝坻断裂之间形成北部平原，为Ⅲ级构造单元蓟宝隆褶。南部为南部平原。

天津市第四系在北部山区零星分布，主要广布于平原区，厚度由近山前不足百米向南变厚，坳陷区己超过400 m，多为冲洪积、冲积、沼、湖积和海积，形成了砂、砂性土和黏性土不规则互层，无论从平面上还是从剖面上其岩性组合较单一(表1)。

表1 天津市第四系地层岩性特征Table 1 Lithology of the Quaternary sediments in Tianjin

1.2水文地质条件

天津市第四系孔隙水广泛分布于广大平原区，按埋藏条件和水文地质特征可分为四种类型：全淡水、浅层淡水、咸水和深层淡水。根据以往的地质、水文勘查成果，第四系孔隙水划分为四个含水组，即第Ⅰ含水组、第Ⅱ含水组、第Ⅲ含水组和第Ⅳ含水组。各含水层特征见表2。

表2 含水组特征Table 2 Characteristics of the aquifers

1.3地温场特征

地壳浅层地温场在垂向上的变化特征受地层岩性、结构、孔隙度、当地气候、地下水活动等多种因素的影响。从上至下依次分为：变温层、恒温层和增温带，天津市恒温带深度一般在30 m左右，“恒温带”之上称为“变温带”，“恒温带”之下称为“增温带”，其温度主要受地壳传导热的影响，地温随深度的增加而增高，不同地区具有不同的地热增温率。

不同地区地温梯度也不同，这主要和控热的地质构造、热储层结构、岩浆和断裂活动及水文地质环境等有密切关系。天津市200 m以浅地层的地温梯度值一般为1.0～4.0 ℃/100 m。选取了不同构造区的增温带(200 m以浅)的地温梯度特征值，见表3。

表3 地温梯度分布特征值Table 3 Distribution characteristics of the geothermal gradient

2 岩土体物性特征

地层热物性特征，主要是岩土体的热导率、比热容和热扩散系数，反映了岩土体的蓄热和导热能力，是影响浅层地热能资源赋存的重要影响因素，是浅层地热能资源开发利用适宜性分区及资源潜力评价的重要参数。

岩土体热物性特征主要通过实验室对岩土体样品测试取得。共完成88个勘查孔，勘查孔的设计深度为120 m，直径为110 mm，根据岩性采集岩土样，对勘查孔进行岩性分析，每个勘查孔采集10～12个岩土样，每个样品长度为25 cm，采样间隔为10 m或者根据地层埋深进行调整，共采集1 076个岩土样。现场按照试验室对岩土样要求进行包装，送往实验室，测试其热导率、比热容以及热扩散系数。

2.1测试方法

岩土热导率和热扩散系数采用微细热探针法进行测量。热探针法作为一种非稳态测量方法，已经广泛应用于多孔介质传热介质研究领域。试验热探针主体采用Φ2 mm的不锈钢管，长度为200 mm。热源采用Φ0.1 mm的绕制漆包康铜丝，外接WYJ-15V型恒流源保证稳定的加热功率。测温传感器采用自制的微细铜-康铜热电偶，通过Julabo型高精度恒温水浴标定，测试数据通过Aglient 34970A型温度采集仪记录和储存，温度测量精度为±0.1 ℃。试验中，岩土热探针的测试结果采用丙三醇和筛分干燥细砂进行标定，从而确定平均修正因子。经过标定后，岩土样品的热导率和热扩散系数测试误差可以控制在±2%之内[13～14]。

比热容采用混合法(绝热恒容槽)，岩土体样品温度测量采用Φ0.05 mm铜-康铜热电偶，标定精度±0.1 ℃，温度数据记录采用HP/Agilent 34970自动采集系统，恒温水浴采用德国Julabo型实验室级加热/冷却循环装置(±0.02 ℃)。测试结果进行修订，修订后的测试结果误差在±5%之内，满足工程设计的需求[13～14]。

2.2岩土体物理性质统计分析

按照岩土体的岩性、物理性质分类，在数理统计(采取方差剔除异常值)方法进行全孔段(120 m)的加权平均，根据勘查孔岩性资料，把天津市岩土体主要概化为黏土、粉质黏土、粉土、粉砂四种，对1 076个热物性测试成果进行统计，分析岩土体热物性参数分布规律，并且初步解释其影响因素。

2.2.1不同地点同一岩性热物性参数对比分析

依据采样的地点，岩土体热物性参数测试结果按照天津市行政区进行统计分析，为了直观表示出不同地点、同一岩性热物性变化，将统计分析的结果绘制曲线图，见图1。

图1 不同地点同一岩性热导率、比热容和热扩散系数对比分析图Fig.1 Comparison of the thermal conductivity, specific heat and thermal diffusion coefficient of different rocks in different sities

通过对比分析图1，能明显看出，热导率总体上呈现黏土>粉质黏土>粉土>粉砂的分布特征，但是对于同一种岩性而言，不同地点热导率变化不大，因此可以认为热导率的大小仅受地层岩性的影响而与分布的位置相关性不大。从图1看，不同岩性、不同地区的比热容和热扩散系数变化不大，只是因所处的区域和岩性不同存在差异，总体上呈现粉土>粉质黏土>粉土>粉砂，应更一步分析其影响因素。

2.2.2同一地点不同岩性热物性参数对比分析

为了分析同一测试点不同岩性之间各物理性质的差异，对蓟县、中心城区及静海勘查孔的岩土体测试数据进行对比分析，分别代表天津市北部、中部和南部地区的地层物理性质特征，见表4[15]。

通过同一地区不同岩性之间热物性对比分析，结果表明：岩土体的热物性参数主要受地层岩性的影响，对比热容而言，黏土最大，粉砂最小，但是热导率刚好相反，黏土最小，而粉砂最大，热扩散系数与热导率规律相同。

2.2.3岩土体热物性平面特征

为了直观看出综合热导率和比热容的平面变化趋势，绘制平面等值线图(图2)。 图2热导率高值区位于蓟县、宝坻和宁河的东部、武清西部、静海南部以及滨海新区的中部，这正好反映水文地质条件与热导率之间的联系。热导率高的地区一般都是地下水动力条件好的地区，该地区单井涌水量也比较大[16]。

天津市热导率的变化规律是：东西方向上两边高中间低。比热容与热导率呈现大致相反的趋势，比热容的高值区主要分布在蓟县、宝坻和宁河西部、武清东部、静海北部以及滨海新区大部分地区，低值区位于宁河、天津市区西部和津南等地。

表4 三个地区不同岩性的热物性参数Table 4 Thermal properties of different rocks in three location

图2 天津市岩土体热导率和比热容平面分布图Fig.2 Plane distribution of the thermal conductivity and specific heat capacity of rock and soil in Tianjin City

3 分析与讨论

岩土体热物性参数影响因素主要是以下几个方面：

受地层岩性的影响，同一种岩性地层热物性参数与所处的位置关系不大(表7)。由表可知，热导率与热扩散系数变化规律一致，热导率大的岩性，热扩散系数也相对较大。

从平面分布规律分析(图2)，不同地区的地层热物性存在明显的差异。蓟县地区热导率在1.54～1.70 W/(m·K)之间，比热容在2 030～2 070 J/(kg·K)之间，热扩散系数在0.52×10-6～0.74×10-6m2/s之间。中心城区热导率在1.54～1.70 W/(m·K)之间，比热容在2 030～2 070 J/(kg·K)之间，热扩散系数在0.52×10-6～0.74×10-6m2/s之间。静海地区热导率在1.34～1.39 W/(m·K)之间，比热容在2 050～2 090 J/(kg·K)之间，热扩散系数在0.52×10-6～0.74×10-6m2/s之间。不同地区差值的区别主要受地下水径流条件影响，在热导率值相对比较的高位置，比热容值相对较小，如热导率为1.70 W/(m·K)，地区比热容值为2 010 J/(kg·K)，热扩散系数为0.52×10-6～0.74×10-6m2/s。

表7 热物性参数对比表Table 7 Comparison of thermal physical properties

4 结论

(1)通过对88个勘查孔，1 076个岩土样测试结果统计分析，结果表明：同一地区不同岩性的比热容，黏土最大、粉砂最小，热导率刚好相反，热扩散系数与热导率规律相同；不同地区同一岩性的热物性参数差别不大。

(2)对比分析热导率和比热容平面分布规律，天津市比热容与热导率呈现大致相反的趋势，热导率高值区位于蓟县、宝坻和宁河的东部、武清西部、静海南部以及滨海新区的中部地区；比热容的高值区位于蓟县、宝坻和宁河西部、武清东部、静海北部以及滨海新区大部分地区。

(3)综合分析天津地区地热地质条件、水文地质条件、第四系岩性特征以及岩土体热物性参数，天津市适宜开发利用地埋管地源热泵系统，适宜的地埋管埋深为100～120 m。

致谢：非常感谢天津市国土资源和房屋管理提供支撑，得以顺利开展项目，采集了大量的样品，完成实验室测试！

[1] 韩再生, 冉伟彦，佟红兵,等. 浅层地热能勘查评价[J]. 中国地质，2007, 34(6):1115-1121.[HAN Z S, RAN W Y, TONG H B,etal. Exploration and Evaluation of shallow geothermal energy exploration[J]. Geology in China, 2007, 34(6):1115-1121. (in Chinese)]

[2] 王秉忱. 科学运用地源热泵技术开发地热能[J]. 水文地质工程地质，2008，35(3): I.[WANG B C. The use of ground source heat pump in the development of shallow geothermal energy[J]. Hydrogeology & Engineering Geology, 2008，35(3): I. (in Chinese)]

[3] 王秉忱. 我国利用水源、地源热泵技术开发浅层地热能呈现好形势[J]. 工程设计与建设，2008 (1):8-11.[WANG B C. The utilization of water source and ground source heat pump technology in the development of shallow geothermal energy can be a good situation[J]. Construction & design for project, 2008 (1):8-11. (in Chinese)]

[4] 赵军，张春雷，李新国,等. U型管埋地换热器三维传热模型及实验对比分析[J].太阳能学报, 2006, 27 (1):63-66.[ZHAO J, ZHANG C L, LI X G,etal. Three dimensional heat transfer model and experimental contrast analysis of U pipe buried heat exchanger[J]. Acta Energiae Solaris Sinica, 2006, 27 (1) :63-66. (in Chinese)]

[5] Pribnow D, Sass J H. Determination of thermal conductivity from deep boreholes[J]. Journal of geophysical research, 1995,100(B6)：9981-9994.

[6] 毕文明，岳丽燕，韩再生,等. 地埋管换热性能综合微缩试验研究[J]. 水文地质工程地质，2014,41(1)：144-148.[BI W M, YUE L Y, HAN Z S,etal. A comprehensive miniature test study of buried tube heat exchange performance[J]. Hydrogeology & Engineering Geology, 2014,41(1)：144-148.(in Chinese)]

[7] 孙宝成，林黎, 张云霞,等. 天津市浅层地热能调查报告[R]. 天津:天津地热勘查开发设计院, 2010.[SUN B C, LIN L, ZHANG Y X,etal. Investigation report of shallow geothermal energy in Tianjin City[R]. Tianjin:Tianjin geothermal exploration and development designing institute, 2010. (in Chinese)]

[8] 于建水，贾佩，郑晓菲,等. 天津市浅层地热能示范工程经济环境效益评价[R]. 天津:天津地热勘查开发设计院，2012.[YU J S, JIA P, ZHENG X F,etal. Tianjin shallow geothermal energy demonstration project economy and environment benefit evaluation[R]. Tianjin: Tianjin geothermal exploration and development designing institute, 2012. (in Chinese)]

[9] 唐永香，李嫄嫄，俞礽安,等. 天津地区浅层地热能赋存条件浅析及前景展望[J]. 地质调查与研究, 2013, 36(3): 213-220.[TANG Y X, LI Y Y, YU R G,etal. Analysis on occurrence conditions and prospect of the shallow geothermal energy in Tianjin area[J]. Geological Survey and Research, 2013, 36(3): 213-220. (in Chinese)]

[10] 唐永香，李嫄嫄，俞礽安,等. 天津滨海新区浅层地热能资源评价及开发利用对策分析[J]. 地质找矿论坛, 2014, 29(4): 622-627.[TANG Y X, LI Y Y, YU R G,etal. Resource assessment, exploitation and utilization measures of the shallow geothermal energy in Tianjin Binhai New area[J]. Contributions to Geology and Mineral Resources Research, 2014, 29(4): 622-627. (in Chinese)]

[11] 刘九龙，郑晓菲，边宗斌,等. 宝坻国土局浅层地热能示范项目系统高效运行与节能减排研究[M]. 天津:天津地热勘查开发设计院，2015.[LIU J L, ZHENG X F, BIAN Z B,etal. The system high efficiently running, saving energy and reducing emission of shallow geothermal energy demonstration project in Baodi land and resources bureau[R]. Tianjin: Tianjin Geothermal Exploration and Development Designing Institute,2015. (in Chinese)]

[12] DB 12/T 608—2015. 天津市浅层地热能地质环境监测技术规程[S].[DB 12/T 608—2015. The regulation for shallow geothermal energy geological environment monitoring in Tianjin[S]. (in Chinese)]

[13] 胡芃,陈则韶. 量热技术和热物性测定[R]. 合肥：中国科学技术大学出版社，2009.[HU P, CHEN Z S. Calorimetry techniques and thermal physical properties determination[R]. Hefei: University of Science and Technology of China Press, 2009. (in Chinese)]

[14] 杨世铭，陶文栓. 传热学[R]. 北京：高等教育出版社，2006.[YANG S M, TAO W S. Heat Transmission Science[R]. Beijing: Higher Education Press, 2006. (in Chinese)]

[15] 韩斯东，金光，毕文明,等. 某地源热泵场地浅层土壤分层热物性响应实验[J]. 水文地质工程地质，2015,42(4)：144-149.[HAN S D, JIN G, BI W M,etal. Test of shallow soil thermal properties of layered thermal response of a ground source heat pump system in a certain place[J]. Hydrogeology & Engineering Geology, 2015,42(4)：144-149. (in Chinese)]

[16] 刘爱华，佟红兵，冉伟彦. 北京某垂直地埋管区地温场变化规律研究[J]. 水文地质工程地质，2016,43(4)：165-170.[LIU A H, TONG H B, RAN W Y. A study of ground temperature changes in a vertical heat exchanger area of Beijing[J]. Hydrogeology & Engineering Geology, 2016,43(4)：165-170. (in Chinese)]

责任编辑：张若琳

AnanalysisofthecharacteristicsofthermalphysicalpropertiesandtheirinfluencingfactorsintheTianjinarea

RUAN Chuanxia1，2, FENG Shuyou2, MOU Shuangxi2, CHENG Wanqing2, ZHAO Sumin2

(1.SchoolofWaterResourcesandEnvironment,ChinaUniversityofGeosciences(Beijing),Beijing100083,China; 2.TianjinGeothermalExplorationandDevelopmentDesigningInstitute,Tianjin300250China)

In order to analyze the regional distribution and plane distribution regulations of thermal physical parameters of the strata in the Tianjin area, the authors use 88 exploration holes, collect 1 076 soil samples, analyze thermal conductivity, specific heat and thermal diffusion coefficient of formation in the laboratory and carries out statistical analyses of the test data. The results show that the heat conductivity of rock and soil body in the Tianjin area is 1.26～1.70 W/(m·K), the specific heat capacity ranges from 2 050 to 2 090 J/(kg·K), and the thermal diffusion coefficient varies from 0.45×10-6to 0.74×10-6m2/s. The specific heat of different rocks in the same region is different: clay is the smallest and silt is the biggest. The situation is contrary for the thermal conductivity. The thermal diffusivity and the specific heat capacity are of the same changes. Difference in thermal physical parameters in different parts of the same rocks is not large. In the Tianjin region the heat capacity and the thermal conductivity are roughly in an opposite trend, i.e., the areas of high value of thermal conductivity is located in the easttern parts of Jixian, Baodi and Ninghe county, the western part of Wuqing, the southern part of Jinghai and the middle part of the Binhai new area; the areas of high value of the specific heat capacity lie mainly in the western part of Jixian, Baodi and Ninghe county, the eastern part of Wuqing, the northern part of Jinghai and the most part of the Binnhai new area.

geothermal geology; thermal conductivity; specific heat; thermal diffusion coefficient; hydrogeology; thermal physical properties; Tianjin area

P314.1

A

1000-3665(2017)05-0158-06

10.16030/j.cnki.issn.1000-3665.2017.05.24

2016-10-10;

2016-12-08

国家自然科学基金项目(41502256)

阮传侠(1978-)，女，硕士，高级工程师，主要从事地热和浅层地热能资源评价，动态监测和回灌等研究。