回填土质材料对地下换热器冻胀特性的影响研究.pdf

第35卷 第14期 农 业 工 程 学 报 Vol.35 No.14 2019年 7月 Transactions of the Chinese Society of Agricultural Engineering Jul. 2019 205 回填土质材料对地下换热器冻胀特性的影响研究 王有镗，郑 斌，王春光，李成宇，毛明明，刘旭阳 （山东理工大学交通与车辆工程学院，淄博 255049） 摘 要：以地源热泵技术在寒区设施农业中的应用为背景，开展地下换热器低温冻胀特性研究。通过岩土埋管冻胀试验，基于冻结半径追踪和管体应变测量，开展冻结区发展规律以及埋管变形收缩特性的研究，对比砂土基和黏土基回填料对冻胀的影响，并进一步考察换热管容积与流阻变化。结果表明，冻结区在进出水2管中心连线方向发展速度大于其垂直方向，受冻胀影响，U型换热管产生椭圆化变形与收缩，砂土基回填比黏土基回填具有更大的冻结范围，但前者换热管变形和收缩程度却小于后者。换热管在100 h内持续由0降至10 后，测得管容积减小率为0.4%，流阻增大率为6.5%，由此可知，埋管冻胀的变形与收缩成为循环流量减小和系统效率下降的原因之一。 关键词：热泵；传热；冻胀；地下换热器；回填料；埋管变形 doi：10.11975/j.issn.1002-6819.2019.14.026 中图分类号：S215; TK521 文献标志码：A 文章编号：1002-6819(2019)-14-0205-07 王有镗，郑 斌，王春光，李成宇，毛明明，刘旭阳. 回填土质材料对地下换热器冻胀特性的影响研究J. 农业工程学报，2019，35(14)：205211. doi：10.11975/j.issn.1002-6819.2019.14.026 http:/www.tcsae.org Wang Youtang, Zheng Bin, Wang Chunguang, Li Chengyu, Mao Mingming, Liu Xuyang. Effects of backfill materials on frost heave characteristic of ground heat exchangerJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(14): 205211. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.14.026 http:/www.tcsae.org 0 引 言 近年来，地源热泵技术在设施农业领域中得以广泛应用1，其中以日光温室的应用最为典型2-3。然而为了达到良好的采光效果，温室的外围护结构通常采用单层玻璃或塑料膜，因此在光照条件不佳时，地下换热器的运行温度受大气温度影响较大。对于中国广大的北方严寒地区而言，地下换热器（埋管）运行温度常持续低于0 ，由此致使换热岩土发生冻胀，冻胀作用不但会挤压埋管发生变形，而且也影响到系统运行效率及安全性4-8。 对于冻胀岩土挤压管道的研究，成果主要体现在油气运输和管渠输水等工程领域中。理论研究方面，学者们主要应用管-土模型计算管土作用9-11、预测管位移12-13、分析管变形应力和应变14-16。在试验研究方面，研究者已针对埋管冻胀变形，开展了许多大规模的工程测试，诸如加拿大Calgary冻胀试验17、法国Caen全尺寸管线冻胀测试12,18、美国Fairbanks高速路管线冻胀试验19-20以及美国UAF冻胀隆起管线测试21-22。此外，学者们也在农业设施的抗冻设计方面进行了一系列的理论探索23-25与验证26-27。近年来，尽管有学者对与地下换热器相似的垂直埋管28和管道温变影响29进行了管土冻胀研究，但由于应用背景不同，对地热热泵工程仍缺乏针对性。目前针对寒区地下换热器运行中易出现的冻胀问题开始收稿日期：2018-02-24 修订日期：2019-06-28 基金项目：国家自然科学基金资助项目（51806130）；山东省自然科学基金资助项目（ZR2017LEE031）；山东省重点研发计划资助项目（2019GHY112076） 作者简介：王有镗，讲师，博士，主要从事浅层地能利用研究。Email：wytsdut.edu.cn 得到关注，最初由加拿大工程师Lenarduzzi等7结合工程实际指出该问题的研究意义，此后中国学者郑平等30提出了考虑冻胀的土壤源热泵水热力耦合数值模型，笔者通过建立岩土冻胀试验系统探讨了换热管形态对其冻胀变形的影响31，同时基于孔隙增长率函数建立了地下换热器的冻胀变形模型32。 本文基于前述研究基础，针对工程中常用的砂土基与黏土基回填料，试验研究竖埋U型管周围土壤（以下简称围土）冻胀以及管截面变形特性；对比2种回填料的影响；考察管变形给循环系统带来的影响，为指导工程应用提供参考。 1 试验装置及方法 1.1 试验系统 本试验模拟地下竖埋U型换热器在饱和岩土环境中的低温运行过程，U型埋管内循环流动低温液体，土壤由此发生冻胀而使换热管受力变形。通过追踪土体内0 冻结锋面发展，研究冻结区发展特性，利用应变测试手段研究换热管截面变形特性。试验系统主要包括岩土槽、围土、U型换热管和冷源循环系统等，如图1所示。 1.2 试验及测试装置 圆柱形岩土槽的高为900 mm，内径为800 mm，槽身下部开4个补水口，控制槽内湿岩土水位，槽内底部铺垫一层150 mm厚砾石，以利于补水均匀渗入，砾石上部添置试验围土（包含地层土壤和回填料），为更好体现冻胀特性，试验采用均匀的饱和围土。用PE板制作环套，使岩土槽外围形成环腔空间，为槽内岩土提供4 农业工程学报（http:/www.tcsae.org） 2019年 206 恒定的边界温度，槽盖和槽底做绝热处理。U型换热管竖向埋置于岩土槽中央，壁厚3.5 mm，外径32 mm，U型间距为80 mm。循环液采用30%乙二醇溶液，由其将热量从岩土槽运送至冷源。 1.地层土壤 2.U型管 3.砂土基回填 4.黏土基回填 5.砾石层 6.排液管 7.岩土槽 8.补水口 9.温控单元 10.循环泵 11.流量计 12.冷源 13.水箱 14.球阀 15.软连接 16.压差计 17.槽盖 18.接线出口 19.环套 1.Ground soil 2.U-pipe 3.Sand-based backfill 4.Clay-based backfill 5.Gravel 6.Bleeder pipe 7.Soil tank 8.Water supply 9.Temperature control unit 10.Pump 11.Flowmeter 12.Cold source 13.Water tank 14.Ball valve 15.Flexible connection 16.Differential pressure gauge 17.Tank cover 18.Wiring channel 19.Annular sleeve a. 结构示意图 a. Structure diagram b. 试验布置 b. Experimental setup 图1 试验系统图 Fig.1 Experimental system 在U型管进出口布置热电偶测试循环液温度，在换热管埋深200 mm（砂土基回填层）处和600 mm（黏土基回填层）处的水平面上分别布置热电偶，以测试围土温度，如图2所示。以进、出水2管中心连线为X方向，与其垂直方向为Y方向，2管的中间位置为原点O，建立X-Y直角坐标系，在X和Y方向上各对称布置10个测点，相邻测点间距为40 mm。电阻应变片测试换热管外壁周向应变，利用应变变化判断管截面的变形，测试位置1、 1和2、2分别在换热管埋深180 mm（砂土基回填层）处和580 mm（黏土基回填层）处的进、出水管上，每处在X和Y方向对称布置、和 4个应变测点。U型压差计接U型管同一水平高度的进口与出口，测量其间压差变化。主要试验测试仪器型号及规格见表1。 注：1、1和2、2为应变测试位置，、和为应变测点，上角标和上角标分别表示进水管和出水管。 Note: 1, 1, 2 and 2 represent the strain test positions, , , and represent the strain test points, superscript and superscript represent the inlet pipe and the outlet pipe, respectively. 图2 应变和温度测试 Fig.2 Strain and temperature test 表1 仪器型号及规格 Table 1 Types and specifications of instruments 仪器Instrument 型号规格Type and specification 温度记录仪 Temperature recorder 横河 DX1012 YOKOGAWA DX1012 热电偶Thermocouple 欧米伽T型 OMEGA T-type 应变采集仪 Strain collector 无线静态型DH3819 Wireless static type DH3819 应变片 Strain gage 康铜箔电阻型 BE120-3AA Constantan foil resistance type BE120-3AA 压差计 Differential pressure gauge 测量范围0500 mm Measurement range 0-500 mm 温控仪Temperature controller Yudian AI-518P 1.3 试验方法 地源热泵工程技术规范中已经指出：细砂和膨润土的混合浆或专用灌浆材料宜作为回填材料，以实现保水防渗功能；然而，由于成本或施工便利等原因，工程中也存在回收使用钻孔屑黏土基回填料的情况，因此，本文选择以上2种回填料进行分析。地层土壤与黏土基回填料选用天然细粒黏土，砂土基回填料选用细砂和膨润土混合物（质量比为7:3），其基本特性参数见表2，模拟回填孔直径为150 mm。试验过程中，对围土分层作填充-夯实-填充处理，回填料的填充以回填高度的一半为界（埋深375 mm处），下半部为黏土基回填，上半部为砂土基回填，如图2所示。最后，通过补水口使土体达到饱和状态。 表2 回填料基本参数 Table 2 Basic parameters of backfills 不同粒径的组成Composition of grain size/% 填料类型 Types of backfill 0.005 mm (0.005 mm, 0.075 mm (0.075 mm, 0.25 mm (0.25 mm, 0.5 mm) 孔隙率Porosity 含水率 Moisture content/% 密度 Density/(kgm-3) 砂土基回填料 Sand-based backfill 4 9 51 36 0.32 17 2117 黏土基回填料 Clay-based backfill 31 42 18 9 0.33 18 2124 第14期 王有镗等：回填土质材料对地下换热器冻胀特性的影响研究 207 试验主要基于变形的应变测量，分析变形的基本形态。初始截面为圆形的换热管（此时应变=0）受到挤压后，如果其截面为椭圆化变形，则椭圆化向外凸出的管壁曲率变大时，该处外壁面将承受拉应力，即表现为拉应变（0）；椭圆化向内收缩的管壁曲率变小时，该处外壁面将承受压应力，即表现为压应变（0），如图3所示。 图3 变形与应变 Fig.3 Deformation and strain 试验过程中，直接测量参数的不确定度可由公式（1）确定 5022221 .dnddd uuuu （1） 式中ud为参数的不确定度，ud1, ud2, , udn为相互独立的不确定因素。 本试验中，温度测量的不确定度为2.02%，应变测量的不确定度为4.08%，液柱差测量的不确定度为4.54%。 2 结果及分析 2.1 冻结区发展分析 使系统连续运行100 h，换热管内流体温度由0 逐渐降低至10 ，如图4所示。由于U型换热管的长度仅为0.9 m，运行过程中其进、出口的温度差异不大（不超过0.5 ），进、出水2管及其轴向温度可视为均匀一致。 图4 换热管进出口温度 Fig.4 Inlet and outlet temperature of pipe 以O点为基准点，定义冻结锋面的0 点与基准点之间的距离为冻结半径，如图5a所示，以衡量冻结区范围。围土内X和Y方向的冻结半径发展如图5b所示，可见冻结半径在不同回填层内发展规律一致。以黏土基回填层为例，X和Y方向的冻结半径发展至56 mm所用时间分别为11和36 h，二者相差25 h，此后2方向冻结半径发展至96和136 mm位置的时间差分别为16和13 h，可见冻结区在X方向的发展速度始终大于Y方向，同时随着冻结范围的扩大，2方向的发展差异逐渐减小。显然，不对称的冻结区发展特性源于换热管的布置方式，进、出水2管位于X方向上，冻结范围较小时，冻结区形态受到管结构的影响较大，随着冻结区边缘（冻结锋面）距离换热管越远，其形态受管结构的影响逐渐减弱。 试验过程中，砂土基回填层的冻结半径发展平均速度（X向：1.64 mm/h；Y向：1.43 mm/h）要略快于黏土基回填层（X向：1.58 mm/h；Y向：1.37 mm/h）。事实上，除了砂土导热性能良好之外，在初始含水率相差不大的情况下，由于砂土基回填料内膨润土中亲水矿物的作用，使得回填料会吸收地层土壤中的水，这种水分迁移作用，可使正在冻结的区域冰晶含量不断增大33，冰的导热系数大于水和土颗粒。因此冰晶较多的砂土基回填层具有更大的冻结范围。 a. 冻结区 a. Freezing area b. 冻结半径对比 b. Comparison of freezing radius 图5 冻结区发展特性 Fig.5 Development characteristic of freezing area 2.2 管截面变形分析 取各应变测试位置处每25 h的应变变化，如图6所示。 图6 管体应变变化 Fig.6 Variations in pipe strains 由图6可以发现，各处测点均为负应变（压应变），且应变值不断增大，然而各测点变化规律却有所差异，农业工程学报（http:/www.tcsae.org） 2019年 208 和测点应变值普遍大于和测点。事实上，所测管体应变取决于2个因素：温度和变形，管壁不断的降温收缩使得压应变值不断增大，而变形使管壁附加了额外的拉压应力。由应变测试结果及前述变形与应变关系可知，相当于和测点附加了压应变，和测点附加了拉应变，管段发生了椭圆化变形，且椭圆长轴在Y向。 冻结区的不均匀发展导致了冻胀力的不均匀产生，由管截面的变形可知，来自X向的冻胀力大于Y向，这也与所测冻结区发展特性相一致。同时，随着运行时间增加，和测点压应变与和测点压应变的差距也在不断增大，表明管截面椭圆化变形程度不断加重。 对比可知，1和1在各个测试时间的变形都要小于2和2位置，处于砂土基回填层的管段，其椭圆化变形程度要小于黏土基回填层的管段。这表明相同条件下，黏土基回填层冻结所产生的冻胀力要比砂土基回填层的大。由于进出水2管温差不大，2管周围冻结区的发展较相似，因此相同埋深位置处，进出水2管段的变形程度差别不大。 2.3 管截面收缩分析 由各测试位置的管体应变（图6）可以看出，管截面在不断椭圆化的同时，也在不断收缩。管体平均应变（m）反映了截面周长的变化率，进而可反映管截面胀缩程度。试验中采用4个测点应变的平均值近似表示测试位置处的管体平均应变（m），如图7a所示。可见在前30 h内，4处测试位置的m基本维持在相同的水平，而大约运行至30 h时，2和2位置的m开始大于1和1位置，且差距逐渐增大，最终前者的m值达到2 50010-6左右，后者仅为1 90010-6左右。表明当冻结发展至一定程度时，黏土基回填层内的管截面收缩程度要比砂土基回填层内的大。 事实上，在低温运行环境下，管截面的收缩取决于温降和冻胀挤压作用的大小。温降引起的截面收缩应变T可由公式（2）求得。 0)T it t （ （2） 式中为HDPE管热胀系数，取值13010-6-134；0t 和it 分别为初始温度和第i时刻温度，。 因此，试验中各测试位置的m与T之差可体现冻胀挤压作用对管体应变的影响，如图7b所示。可见在前30 h内，4 处测试位置的冻胀应变值不断增大且增幅水平相当，从30 h至100 h，1和1处冻胀应变值由40010-6仅下降至50010-6左右，而在2和2处冻胀应变值却降至1 00010-6以下。表明当冻结区发展超过一定范围，黏土基回填对换热管的冻胀作用将大于砂土基回填，且差距逐渐增大。砂土基回填料中的膨润土存在遇水膨胀软化特性，从而使其强度降低35，冻胀力对换热管的作用相对减弱，因此在砂土中适量添加膨润土，可减轻冻胀产生的换热管变形与收缩程度。对于寒区地源热泵工程，从地下换热器安全运行角度而言，按照规范使用砂土基回填料显得尤为必要。 a. 管体平均应变m a. Mean strains in pipe surface m b. 冻胀作用应变m-T b. Frost heave strains m-T 图7 管体收缩应变 Fig.7 Contraction strains in pipe surface 2.4 管容积与流阻分析 在运行温度不断降低和冻胀作用力不断加大的情况下，换热管降温收缩和椭圆化变形会使截面面积减小，从而导致管容积减小。对于内径为25 mm的换热管，结合4个应变测点的平均值m，可求得冻胀变形后管体容积V，再根据初始管容积V0，可得容积变化量=V0V，如图8所示。可见随着运行时间的增加，管容积呈现不断减小的趋势，运行最终，换热管有约3.3 mL的容积减少量，与计算管段（管长750 mm）的初始容积V0=736 mL相比，减少了0.4%。 图8 换热管容积减少量 Fig.8 Volume decrease of heat exchange pipe 换热管的变形与收缩同样也会增大系统的流动阻力，系统流阻的变化可通过U型压差计的液柱差来反映。冻胀试验过程中的液柱差变化如图9所示，可见液柱差不断增大，且增幅随运行时间有增加的趋势，这与不断增大的管体变形程度有关。在100 h的运行时间内，液柱差由137 mm增至146 mm，表明U型换热管的流阻有6.5%的增加量。 事实上，换热管由于冻胀变形造成的容积减小以及第14期 王有镗等：回填土质材料对地下换热器冻胀特性的影响研究 209 流阻增大问题，在文献6和7中所阐述的工程问题中已有提及，由此所带来的诸如循环液溢流、系统流量减小以及COP（性能系数，coefficient of performance）下降等问题不容忽视，在工程设计中应加以考虑。 图9 液柱差增幅 Fig.9 Increase in liquid column difference 3 结论与讨论 本试验基于U型埋管地下换热器的持续降温运行过程，针对不同回填土质材料（砂土基和黏土基），开展围土冻结区发展以及埋管变形特性等基础研究，得出如下结论： 1）受U型结构影响，围土冻结区域不均匀发展，冻结区在进出水2管中心连线方向发展速度大于其垂直方向，该差距随冻结范围增大而减小，砂土基回填可较黏土基回填具有更大的冻结范围。 2）换热管截面产生椭圆化变形和收缩，椭圆短轴在2管中心连线方向，随着冻结区发展，变形与收缩程度不断加重，黏土基回填料中由冻胀产生的换热管变形与收缩程度更为严重，对于寒区地源热泵工程，使用添加膨润土的砂土基回填料显得尤为必要。 3）岩土冻胀所致的换热管截面变形及收缩，会引发管容积减小以及流阻增大问题，运行100 h、温降10 的试验条件下测得管容积减小率为0.4%，流阻增大率为6.5%，该问题应在工程设计中加以考虑。 事实上，冻胀特性主要取决于土壤、含水率和温度。工程中回填土质呈现多样化特点，土壤粒径和结构不尽相同，即使对于相同材质的回填料，土体强度也会由于回填深度不同而产生差异；孔隙含水率的变化可致使回填料热学和力学特性产生差异，同样，冻融变化也会对孔隙结构及含水率产生影响；此外，环境温度的变化也会影响地下换热器的低温运行模式。因此，上述因素均会影响地下换热器的结构冻变特性，进而作用于管土换热效率，限于篇幅，拟在后续对相关问题进行深入探讨。 参 考 文 献 1 方慧，杨其长，王柟，等. 浅层地热源节能技术及其在设施农业中的应用J. 农业工程学报，2008，24(10)：286290. Fang Hui, Yang Qichang, Wang Nan, et al. Geothermal technology and its applications in protected agricultureJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CSAE), 2008, 24(10): 286290. (in Chinese with English abstract) 2 孙维拓，张义，杨其长，等. 温室主动蓄放热-热泵联合加温系统热力学分析J. 农业工程学报，2014，30(14)：179188. Sun Weituo, Zhang Yi, Yang Qichang, et al. Thermodynamic analysis of active heat storage-release associated with heat pump heating system in greenhouseJ. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2014, 30(14): 179188. (in Chinese with English abstract) 3 石惠娴，任亦可，孟祥真，等. 植物工厂水蓄能型地下水源热泵供热系统节能运行特性研究J. 农业工程学报，2018，34(23)：157163. Shi Huixian, Ren Yike, Meng Xiangzhen, et al. Research on energy-saving operating characteristics of water storage groundwater source heat pump heating system in plant factoryJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(23): 157163. (in Chinese with English abstract) 4 李南生，谢利辉，陈薛浩. 寒区浅埋输油管线冻胀安全性分析J. 结构工程师，2008，24(1)：3540. Li Nansheng, Xie Lihui, Chen Xuehao. Frozen-heaving of frost-soil bed of shallow-buried oil-pipeline in cold regionJ. Structural Engineers, 2008, 24(1): 3540. (in Chinese with English abstract) 5 胡宗柳，吴明，陈杨，等. 输油管道冻胀安全性研究现状与趋势J. 油气储运，2011，30(12)：881883. Hu Zongliu, Wu Ming, Chen Yang, et al. Progress in the safety study of frost heaving of oil pipelineJ. Oil and Gas Storage and Transportation, 2011, 30(12): 881883. (in Chinese with English abstract) 6 Gabrielsson A, Lehtmets M, Moritz L, et al. Heat storage in soft clay field tests with heating (70) and freezing of the soilR. Sweden: Swedish Geotechnical Institute, 1997. 7 Lenarduzzi Frank J, Cragg Chris B H, Radhakrishna H S. Importance of grouting to enhance the performance of earth energy systemsJ. ASHRAE Transactions, 2000, 106: 424434. 8 Tafreshi S N M, Khalaj O. Laboratory tests of small-diameter HDPE pipes buried in reinforced sand under repeated-loadJ. Geotextiles and Geomembranes, 2008, 26: 145163. 9 Selvadurai A P S, Hu J, Konuk I. Computational modelling of frost heave induced soil-pipeline interaction I. Modelling of frost heaveJ. Cold Regions Science and Technology, 1999, 29(3): 215228. 10 Selvadurai A P S, Hu J, Konuk I. Computational modelling of frost heave induced soil-pipeline interaction II. Modelling of experiments at the Caen test facilityJ. Cold Regions Science and Technology, 1999, 29(3): 229257. 11 Michalowski R L, Zhu M. Frost heave modelling using porosity rate functionJ. International Journal for Numerical and Analytical Methods in Geomechanics, 2006, 30(8): 703722. 12 White T L. Pipelines in Permafrost and Freezing Ground, Engineering Resource Library and Database IndexesM. Canada: Permafrost Environmental Consulting Inc, 2006. 农业工程学报（http:/www.tcsae.org） 2019年 210 13 Susan A T, Ian D, Moore, Muge Balkaya. Parametric study of frost-induced bending moments in buried cast iron water pipesJ. Tunnelling and Underground Space Technology, 2016, 51: 291300. 14 Wu Yaping, Sheng Yu, Wang Yong. Stresses and deformations in a buried oil pipeline subject to differential frost heave in permafrost regionsJ. Cold Regions Science and Technology, 2010, 64(3): 256261. 15 Wang Jiabin, Niu Ditao, He Hui. Frost durability and stressstrain relationship of lining shotcrete in cold environmentJ. Construction and Building Materials, 2019, 198(2): 5869. 16 Razaqpur A G, Wang D. Frost-induced deformations and stresses in pipelinesJ. International Journal of Pressure Vessels and Piping, 1996, 69(2): 105118. 17 Carlson L E, Nixon, J F. Subsoil investigation of ice lensing at the Calgary, Canada, frost heave facilityJ. Canadian Geotechnical Journal, 1988, 25(2): 307319. 18 Williams P J, Riseborough D W, Smith M W. France-Canada joint study of deformation of an experimental pipe line by differential frost heaveJ. International Journal of Offshore and Polar Engineering, 1993, 3(1): 5660. 19 Nixon J F. The role of convective heat transport in the thawing of frozen soilsJ. Canadian Geotechnical Journal, 1975, 12(3): 425429. 20 Nixon J F. Discrete ice lens theory for frost heave in soilsJ. Canadian Geotechnical Journal, 1991, 28(6): 843859. 21 Huang S L, Bray M T, Akagawa S, et al. Field investigation of soil heave by a large diameter chilled gas pipeline experiment, Fairbanks, AlaskaJ. Journal of Cold Regions Engineering, 2004, 18(1): 234. 22 Kim K, Zhou W, Huang S L. Frost heave predictions of buried chilled gas pipelines with the effect of permafrostJ. Cold Regions Science and Technology, 2008, 53(3): 382396. 23 谢崇宝，白静，吴志琴，等. 季节性冻土区灌溉管道排空防冻模式设计J. 农业工程学报，2018，34(21)：8288. Xie Chongbao, Bai Jing, Wu Zhiqin, et al. Design of irrigation pipeline emptying anti-freezing mode in seasonal frozen soil regionJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(21): 8288. (in Chinese with English abstract) 24 唐少容，王红雨，潘 鑫，等. U 形混凝土衬砌结构冻胀性能离心模型试验研究J. 农业工程学报，2019，35(1)：157163. Tang Shaorong, Wang Hongyu, Pan Xin, et al. Frost heave performance of U-shaped canal concrete lining based on centrifuge model testJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 157163. (in Chinese with English abstract) 25 白静，谢崇宝，吴志琴. 季节性冻土区管道浅埋换填防冻模式研究J. 水利学报，2018，49(5)：588597. Bai Jing, Xie Chongbao, Wu Zhiqin. Research on anti-freezing pattern of shallow buried pipe by replacement filling of pipe trench in seasonal frozen areaJ. Journal of Hydraulic Engineering, 2018, 49(5): 588597. (in Chinese with English abstract) 26 戴佳琦. 寒区村镇饮用水输配管道的