冰川冻土, 2021, 43(6): 1617-1627 doi: DOI：10.7522/j.issn.1000-0240.2021.0108

第二次青藏高原综合科学考察研究

高山区雪堆储雪的实验和模拟计算

王兴,1,2, 王飞腾,1, 任贾文1,2, 秦大河1,2

1.中国科学院 西北生态环境资源研究院 冰冻圈科学国家重点实验室，甘肃 兰州 730000

2.中国科学院大学，北京 100049

Experiment and numerical simulation of snow storage in alpine region

WANG Xing,1,2, WANG Feiteng,1, REN Jiawen1,2, QIN Dahe1,2

1.State Key Laboratory of Cryospheric Science，Northwest Institute of Eco-Environment and Resources，Chinese Academy of Sciences，Lanzhou 730000，China

2.University of Chinese Academy of Sciences，Beijing 100049，China

通讯作者: 王飞腾，研究员，主要从事冰冻圈科学研究. E-mail： wangfeiteng@lzb.ac.cn

编委: 庞瑜

收稿日期: 2020-07-01   修回日期: 2021-12-12  

基金资助: 中国科学院重点部署项目.  KFZD-SW-323
中国科学院战略性先导科技专项.  XDA20060201
第二次青藏高原综合科学考察研究项目.  2019QZKK0201

Received: 2020-07-01   Revised: 2021-12-12  

作者简介 About authors

王兴，硕士研究生，主要从事冰雪物理研究.E-mail：wangxing181@mails.ucas.ac.cn , E-mail：wangxing181@mails.ucas.ac.cn

摘要

随着全球气候变暖，滑雪运动呈现出对气候变化高度的敏感性和依赖性，各地滑雪季将不同程度地缩短。储雪作为一种应对气候变化的方法在雪务保障中逐渐得到研究和应用。在新疆阿勒泰地区吉木乃县开展储雪实验，应用谐波反应法研究了雪堆融化量与外界空气温度、太阳辐射和绝热保温结构热学性能之间的定量关系。实验中，覆盖两层绝热保温材料雪堆1的平均融化量为18.0 kg·d-1（相当于初始质量的0.85%），未覆盖绝热保温材料雪堆2的平均融化量为120.8 kg·d-1（相当于初始质量的6.67%）。模拟期间观测到的雪堆1质量减少了1 438.0 kg，对应的模拟值为1 520.8 kg。谐波反应法可以为评价绝热保温结构的热学性能和提前确定储雪量提供重要的参考依据。由于许多物理过程都未考虑，导致模拟的雪堆融化量与实际融化量之间存在不确定性。雪堆所用绝热保温结构的性能可以用反射率、总传热系数、衰减度和延迟时间来衡量。

关键词： 储雪 ; 综合空气温度 ; 绝热保温结构 ; 谐波反应法

Abstract

With global warming， skiing sports show sensitivity and dependence on climate change； the skiing season will be shortened to different degrees worldwide. As a coping method， snow storage has been gradually studied and applied in the skiing industry. In this paper， a snow storage experiment was carried out in Jiminay County （47°08′ N，85°35′ E，2 915 m a.s.l.）， Altay， Xinjiang in May 2019. The harmonic response method studied the quantitative relationship between the amount of snowpack melting， external air temperature， solar radiation， and thermal properties of the insulation structure. The average melting rate of Snowpack 1 with two layers insulating materials was 18.0 kg·d-1 （equal to 0.85% of the initial weight）. In contrast， the average melting rate of Snowpack 2 without insulating material was 120.8 kg·d-1 （equal to 6.67% of the initial weight）. Glass wool blanket and aluminum foil reflective film on the upper surface of Snowpack 1 can effectively reduce the amount of snowpack melting. The external air temperature and solar radiation have good periodic variation； the 12 order Fourier series can simulate the change of solar and external air temperature. The weight of Snowpack 1 observed decreased by 1 438.0 kg， and the corresponding simulated value was 1 520.8 kg. The harmonic response method can provide momentous significance for evaluating the performance of the thermal insulation structure and determining the amount of snow storage in advance. Because many physical processes are not considered， there is still uncertainty between the amount of simulated melting and the amount of actual melting. The performance of thermal insulation structures used in snow storage can be measured by reflectivity， total heat transfer coefficient， amplitude decrement and time lag.

Keywords： snow storage ; solar and external air temperature ; insulation structure ; harmonic response method

PDF (4333KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王兴, 王飞腾, 任贾文, 秦大河. 高山区雪堆储雪的实验和模拟计算[J]. 冰川冻土, 2021, 43(6): 1617-1627 doi:DOI：10.7522/j.issn.1000-0240.2021.0108

WANG Xing, WANG Feiteng, REN Jiawen, QIN Dahe. Experiment and numerical simulation of snow storage in alpine region[J]. Journal of Glaciology and Geocryology, 2021, 43(6): 1617-1627 doi:DOI：10.7522/j.issn.1000-0240.2021.0108

0 引言

储雪是指按照低环境影响和经济性原则，将雪储存起来以度过一年中温暖的季节［1-2］。随着人类历史的发展，储雪的应用目的不断发生变化。早期储雪被广泛应用于食品储存和房间制冷［3-4］，尤其是随着气候变暖［5］，对于建筑制冷的需求不断增加［6-7］；近年来，随着冰雪运动的发展，储雪对滑雪产业的作用越发重要［8-9］，能够增加雪场运营时间、降低投资成本和节能环保等。许多冬季雪上运动赛事都会提前进行储雪作业，尤其是近几十年，储雪成为雪上运动赛事规划必不可少的一部分，例如索契和平昌冬奥会都提前进行了储雪［10-11］。

目前最常见的储雪方法是将雪储存在地面上并覆盖绝热保温材料［12］。为了保证储雪结束时的用雪安全，需要对储雪融化量进行模拟，并指导储雪方案的设计。Olefs等［9］在奥地利高海拔地区使用SNOWPACK模型模拟了绝热保温层对积雪能量平衡的影响。Grünewald等［1］发现SNOWPACK模型在点尺度上可以很好地模拟雪堆的高度变化，但是对雪堆三维变化的模拟误差较大。Lintzén等［2］在瑞士使用稳态方法计算了由于地面温度、降雨和空气温度所导致雪堆融化的量。总的来说，目前储雪实验研究较少，主要在高纬度地区开展，缺乏中低纬度地区的储雪技术和经验；储雪材料的选择多依靠个人经验，缺乏理论指导；储雪结果关注质量变化，而忽视雪质的好坏；另外，国外在储雪技术方面对国内进行技术封锁。

储雪所用绝热保温结构的作用是减少由外界环境进入雪堆的热量。绝热保温结构的传热是一个复杂的非稳态过程，其温度场在不断变化，因此在计算绝热保温结构的传热量时需要使用非稳态的导热方法，常见的非稳态导热方法有反应系数法［13-15］、谐波反应法［14-15］、有限差分法［15-17］和有限元法［18-19］。谐波反应法的物理意义明显，不需要反复迭代运算，是国内计算室内空调冷负荷问题的主要方法［14］。外界气象条件具有不同时间尺度的变化特征，但从逐日来看，其大致以24小时为周期发生变化。因此本文基于平均逐时气象数据，以24小时为周期，使用谐波反应法计算雪堆的融化量并分析储雪过程中绝热保温结构的传热特征。

1 方案与数据

1.1　实验地点

本次储雪实验于2019年5月25日开始，地点位于阿勒泰山南麓，萨吾尔山北部，距离木斯岛冰川北侧6 km，新疆阿勒泰地区吉木乃县正南39 km处的高山区气象综合观测场附近（47°08′ N， 85°35′ E，海拔2 915 m）（图1）。该地区属于中国三大稳定积雪区之一，冬季降雪丰富，为本次实验提供了充足的物质基础。其降水主要受西风气流的控制，其次受北冰洋冷湿气流的影响，属于典型的北温带大陆性寒冷气候［20-21］。实验地点所在的阿勒泰地区滑雪历史悠久，滑雪产业发达，到2013年阿勒泰地区滑雪场已达30多家［22］。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   储雪位置

Fig.1   Location of snow storage


1.2　气象及雪堆数据

气象数据由高山区气象综合观测场的自动气象站提供，模拟所使用的气象要素包括气温、入射和出射的短波辐射。数据采集器（CR1000）每半小时存储一次气象数据。雪堆的初始密度及含水率数据使用SLF积雪传感器进行测量。雪堆质量变化数据使用5个压力传感器为一组进行测量。由于传感器技术故障，气象数据于2019年6月22日开始获取至8月25日结束。具体传感器类型和参数见表1。

表1   传感器类型和参数

Table 1  Sensor information and technical specifications

参数传感器名称及型号精度单位架设高度
气温温湿度传感器HC2S3±0.1 ℃℃2.1 m
短波辐射四分量辐射传感器CNR4±10%W·m-22 m
雪堆密度SLF积雪传感器—kg·m-3—
雪堆含水率—vol%—
雪堆质量悬臂梁式压力传感器HZC-30A±0.1%kg0.3 m

新窗口打开| 下载CSV


1.3　实验方案

储雪实验自2017年就已经在河北延庆石京龙滑雪场开展，其主储雪堆上表面绝热保温结构为表层土工织布（0.5 cm厚度）+多层铝箔PE保温层（0.5~3.0 cm厚度）+底部土工织布的三元结构。2018年又在河北崇礼万龙滑雪场开展了四种储雪堆覆盖方案的对比研究实验［23］。本研究基于前期的经验，选择了较优的储雪方案，在吉木乃县开展实验。

本次实验设计有两个雪堆（图2），其形状为正四棱锥，底面边长为2.8 m，棱锥的高为1.5 m，棱锥每个面的坡度为49°，四个面的朝向分别为南、西、北和东。雪堆1上表面覆盖绝热保温材料，雪堆2上表面则不覆盖任何材料。为了评价雪堆上表面覆盖材料的绝热保温效果，两个雪堆底部均置于压力传感器之上并与地面分离，其分别由铁架、铁皮、挤塑聚苯乙烯泡沫塑料保温板和铝箔反射膜组成。针对以往储雪实验中使用硬性绝热保温材料所带来的施工困难，以及与雪堆贴合性差、易滑落等问题［24］，本次实验雪堆1上表面使用玻璃棉毡软性材料。高山地区，太阳辐射是积雪融化的重要能量来源［8，25］，因此在玻璃棉外侧覆盖铝箔反射膜以降低雪堆太阳辐射得热量。铝箔反射膜外覆遮阳网，起到固定内侧绝热保温材料和降低内侧太阳辐射的作用。各材料具体的参数见表2。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   高山区气象观测场附近的储雪实验

Fig.2   Snow storage near the meteorologicalobservation site in alpine region


表2   材料的物理特征

Table 2  Physical characteristics of materials in this study

材料厚度d/m密度ρ/（kg·m-3）导热系数λ/（W·m-1·K-1）比热容c/（kJ·kg-1·K-1）
铁板0.005007 85062.8000.46
挤塑聚苯乙烯泡沫塑料保温板0.10000280.0301.50
玻璃棉毡0.06000120.0581.34
铝箔反射膜0.000042 700237.0000.88

新窗口打开| 下载CSV


2 方法

2.1　模拟方案

融雪期，自然状态下的积雪与周围环境的热量交换包括雪面的净辐射通量、感热通量、潜热通量、地热通量和雨水热通量［26-27］。储雪所用绝热保温结构的特性直接或间接影响着雪堆的能量平衡。绝热保温结构的辐射特性影响雪面的净辐射通量；绝热保温结构的总传热系数影响雪面的感热通量和地热通量；绝热保温结构的渗透率影响雪面的潜热通量和雨水热通量。雪堆上部热量和水分的传输过程如图3所示。影响雪堆融化速率的外界气象要素主要有两项，太阳辐射和外界空气温度［28-29］。雪堆融化过程中，外界环境通过绝热保温结构向内部的传热可以分为三个过程：绝热保温结构外表面的吸热过程，绝热保温结构的导热过程和绝热保温结构内表面的放热过程［30］。外界环境通过绝热保温结构向雪堆传递热量主要通过覆盖层的热传导。雪堆吸热后，冷储减少，雪堆温度升高。本实验因为雪堆完成覆盖一段时间后，才获得各种有效观测数据，且模拟初期雪堆质量就已经开始减少，亦即雪堆内部已经有融化且通过融水流失使雪堆质量减少，所以假定雪堆各层和内部空气层的温度均为0 ℃，从而使绝热保温结构的传热量全部用于雪堆的融化。当雪堆温度升高为0 ℃时，雪堆开始融化，大部分融水自上而下进行水分的迁移，但有时间滞后，而且一部分融水被雪堆截留。融水从雪堆中流出会引起雪堆质量的变化，本实验中假定融水能够全部流出雪堆。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   雪堆上部热量和水分的传输

Fig.3   Transportation of heat and water on theupper part of the snowpack


2.2　斜面太阳辐射总量逐时计算模型

斜面太阳辐射总量I（W·m-2）由三部分组成［15］，直接辐射量IDβ（W·m-2），散射辐射量Idβ（W·m-2），地面反射辐射量Irβ（W·m-2）。

I=IDβ+Idβ+Irβ=IhRb+DhRd+Rh1-cos β2

式中：Ih为水平面的直接辐射量（W·m-2）；Rb为直接辐射转换系数；Dh为水平面的散射辐射量（W·m-2）；Rd为散射辐射转换系数；Rh为地面反射的辐射量（W·m-2）；β为斜面的倾角。

直接辐射转换系数Rb为：

Rb=cos θcos z

式中：θ为斜面太阳入射角（°）；z为太阳天顶角（°）。

散射辐射转换系数Rd的计算采用各向异性模型［31］：

Rd=IhIERb+121-IhIE1+cos β

式中：IE为水平面天文辐射量（W·m-2）。

2.3　非稳态热传导模型

谐波反应法是建立在绝热保温结构导热方程经典求解的基础上，将外界综合空气温度视为以24小时为周期的不规则周期函数（图4），使用傅立叶级数对外界综合温度进行表达，并引入衰减和延迟的概念。谐波反应法计算绝热保温结构的传热量分为两部分，一部分为外界平均综合温度与绝热保温结构内部温度之差所引起的传热量；另一部分是由于外界综合温度相对于平均温度的波动所引起的绝热保温结构内表面波动而产生的附加传热量。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   雪堆1上表面各材料层表面蓄热系数和蓄热系数的计算

Fig.4   Calculation of surface heat storage coefficient and heatstorage coefficient of each material layer on theupper surface of Snowpack 1


热量在固体中的传播服从傅里叶基本导热微分方程。一维无内热源非稳态导热微分方程为［32］：

∂T∂t=a∂2T∂x2

式中：a为热扩散率。

外界综合空气温度Te（℃）可表述为［30］：

Te=Tair+ρ(1-τ)Iαe

式中：Tair为外界空气温度（℃）；ρ为铝箔反射膜的吸收率（0.2）；τ为遮阳网的遮阳率（0.8）；αe为绝热保温结构的外表面换热系数，取其工程上的计算数值为23.26 W·m-2·K-1［33］。

外界综合空气温度Te展开为傅里叶级数的表达式为［14］：

Te=T¯e+∑n=1kAncos (ωnt-φn)

式中：T¯e为外界综合温度的平均值（℃）；An为第n阶余弦波外扰的振幅；ωn为第n阶余弦波外扰的频率；φn为第n阶余弦波外扰的初相位（弧度，rad）。

通过绝热保温结构的导热量计算公式为［14］：

Q=KFT¯e-Ti+αiF∑n=1kAnvncos (ωnt-φn-εn)

式中：K为绝热保温结构总传热系数（W·m-2·K-1）；F为绝热保温结构的面积（m2）；αi为绝热保温结构的内表面换热系数，取其在工程应用上的数值为 7.5 W·m-2·K-1［34］；Ti为绝热保温结构内部空气层的温度（0 ℃），vn和εn分别为绝热保温结构的衰减倍数和延迟时间（弧度，rad）。

雪堆的融化量计算公式为：

G=QL

式中：G为雪堆融化量（kg）；L为雪的融化潜热（MJ·kg-1）。

3 结果

3.1　气象数据

高山区综合观测场中的自动气象站观测结果显示，2019年6月22日至8月25日的平均气温为7.5 ℃，逐日平均气温和气温日较差变化均较大，其范围分别为1.2~13.4 ℃和1.8~10.8 ℃，原因可能与高海拔山区地带大气逆辐射的增温效果较弱和山区小气候复杂多变有关。水平面太阳辐射总量和地面反射辐射量的平均值分别为250.9 W·m-2和 33.5 W·m-2。2019年6月22日至8月25日，水平面太阳辐射总量为1 408.9 MJ·m-2（表3），地面反射的太阳辐射量为188.2 MJ·m-2。模拟期间（2019年6月22日至8月25日）所需要的逐时气象数据见图5。

表3   太阳辐射总量

Table 3  The global solar radiation

辐射面直接辐射量/（MJ·m-2）散射辐射量/（MJ·m-2）地面反射辐射量/（MJ·m-2）太阳辐射总量/（MJ·m-2）
南向斜面739.6512.2124.51 376.3
西向斜面576.0481.9124.51 182.4
北向斜面-344.7124.5469.2
东向斜面747.5504.9124.51 376.9
水平面810.4598.5-1 408.9

新窗口打开| 下载CSV


图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   逐时气象数据

Fig.5   Hourly meteorological data ［external air temperature （a）， global solar radiation on the horizontal surface （b）， reflected radiation on the ground （c）］


3.2　模拟结果与分析

3.2.1　各朝向斜面太阳辐射总量

通过斜面太阳辐射总量逐时模型计算得到雪堆南（S）、西（W）、北（N）和东（E）四个不同朝向斜面的太阳辐射总量（直接辐射量、散射辐射量、地面反射的辐射量）的逐时值（图6）。模拟结果显示，雪堆各朝向斜面太阳辐射总量具有很好的周期性变化特征，可以发现各朝向斜面的辐射变化规律与地球绕太阳的运行规律相对应［图6（a）、6（b）、6（c）］。北向斜面无太阳直接辐射［图6（b）］，北向斜面的太阳辐射总量仅由散射辐射量和地面反射辐射量两部分组成。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   各朝向斜面平均日太阳辐射量的逐时值 （S：南向斜面；W：西向斜面；N：北向斜面；E：东向斜面）

Fig.6   Hourly average daily solar radiation of each inclined plane （S： south-facing surface tilted； W： west-facing surface tilted； N： north-facing surface tilted； E： east-facing surface tilted； B： bottom surface） ［hourly global solar radiation of each inclined plane （a）， hourly direct solar radiation of each inclined plane （b）， hourly diffuse solar radiation of each inclined plane （c）， hourly reflected radiation of each inclined plane （d）］


实验期间雪堆各朝向斜面的太阳辐射总量都未超过水平面的太阳辐射总量（表3）。斜面太阳辐射总量最大的方位角受地区纬度、气象条件［35］和日期的影响，在本实验中，四个朝向斜面以东向斜面的太阳辐射总量最多，南向斜面的太阳辐射总量与东向斜面相接近。从各斜面太阳辐射总量的组成比例来看，南、西和东向斜面太阳直接辐射总量占比最大，分别占太阳总辐射量的53.7%、48.7%和54.3%。北向斜面散射辐射量的占比最大，达到73.5%。雪堆四个朝向斜面中地面反射辐射量的占比最小。

3.2.2　外界综合空气温度

外界综合空气温度的大小是由外界空气温度和太阳辐射共同决定的，各朝向斜面平均日外界逐时综合空气温度的计算结果见表4。00：00—06：00和22：00—23：00（北京时间，下同）各朝向斜面的太阳辐射总量为0，因此这两个时间段的各斜面综合空气温度相等，并在06：00达到最低值。受太阳辐射的影响，南、西、北和东向斜面外界综合空气温度有不同幅度的提升，其外界平均综合空气温度分别为7.9 ℃、7.8 ℃、7.6 ℃和7.9 ℃（外界平均空气温度为7.5 ℃）。由于外界综合空气温度的变化规律不是一个完美的谐波曲线，同时为了提高傅里叶级数拟合的精度，因此在本文中将四个朝向斜面的外界综合空气温度以及外界空气温度展开为12阶傅里叶级数来进行叠加拟合（图7），拟合后的温度与原温度之间的差值小于0.01 ℃。

表4   各朝向斜面平均日外界逐时综合空气温度

Table 4  Hourly average daily solar and external air temperature of each inclined plane

北京时间外界综合空气温度/℃北京时间外界综合空气温度/℃
南向斜面西向斜面北向斜面东向斜面南向斜面西向斜面北向斜面东向斜面
00：006.56.56.56.512：0010.39.69.410.4
01：006.36.36.36.313：0010.59.99.510.3
02：006.16.16.16.114：0010.410.19.610.1
03：005.95.95.95.915：0010.210.19.49.7
04：005.75.75.75.716：009.910.09.49.5
05：005.65.65.65.617：009.79.99.29.2
06：005.45.45.45.418：009.19.58.98.9
07：006.06.06.06.119：008.69.18.48.4
08：007.17.07.07.720：008.08.68.08.0
09：008.38.18.19.021：007.47.67.47.4
10：009.28.68.69.822：006.86.86.86.8
11：009.99.19.110.223：006.66.66.66.6

新窗口打开| 下载CSV


图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   各朝向斜面外界综合空气温度傅里叶级数拟合（A：外界空气温度；S：南向斜面；W：西向斜面；N：北向斜面；E：东向斜面）

Fig.7   Fourier series fitting of the solar and external air temperature in each inclined plane （A： external air temperature； S： south-facing surface tilted； W： west-facing surface tilted； N： north-facing surface tilted； E： east-facing surface tilted）


3.2.3　逐时传热量和雪堆融化量

外界环境通过绝热保温结构的热量可以分为两个部分，雪堆1上表面绝热保温结构的导热量（第一部分）和雪堆底部绝热保温结构的导热量（第二部分）。非稳态方法计算第一部分使用的外界环境参数为各朝向斜面的外界综合空气温度，第二部分使用的外界环境参数为外界空气温度。雪堆上表面绝热保温结构各朝向斜面的面积均相等，为 2.8 m2，底部面积为7.8 m2。图8显示有两个时间段雪堆底部绝热保温结构的导热量大于各朝向斜面的导热量，分别为00：00—07：00和19：00—23：00。由于雪堆底部和上表面绝热保温结构的延迟时间不同，导致第一部分与第二部分逐时曲线的波谷位置不同（图8）。绝热保温结构南、西、北、东朝向斜面逐时导热量分别在14：00、15：00、14：00和13：00达到峰值，绝热保温结构底部逐时导热量在18：00达到峰值。实验期间外界环境通过绝热保温结构的热量中第一部分占比为79.7%，第二部分占比为20.3%，表明通过绝热保温结构的热量中第一部分是影响雪堆融化量的主要因素，这与其他储雪实验结果一致［2］，经分析发现，原因主要有两个：雪堆上表面绝热保温结构的面积大于底部的面积；雪堆上表面绝热保温结构的边界条件比底部更有利于增加传热量。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   绝热保温结构的逐时导热量（S：南向斜面；W：西向斜面；N：北向斜面；E：东向斜面；B：底部）

Fig.8   Hourly conduction of thermal insulation structure（S： south-facing surface tilted； W： west-facing surface tilted； N： north-facing surface tilted； E： east-facing surface tilted； B： bottom surface）


储雪实验于2019年5月25日开始，雪堆1初始质量为2 115.0 kg，雪堆2的初始质量为1 812.0 kg，雪堆2在储雪开始15天后于6月8日全部融化（图9）。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   雪堆的状态

Fig.9   The state of the snowpack


本次实验中，上表面覆盖绝热保温材料的雪堆1平均融化量为18 kg·d-1（相当于初始质量的0.85%），未覆盖绝热保温材料的雪堆2平均融化量为120.8 kg·d-1（相当于初始质量的6.67%）。由于传感器技术故障，模拟从6月22日开始，模拟期间雪堆1质量减少了1 438.0 kg，对应的模拟值为 1 520.8 kg（6月22日至8月25日）。

4 讨论

非稳态导热问题的求解，本质上就是在定解条件下求解导热微分方程［36］。谐波反应法将边界条件的离散数据转化为傅里叶级数，引入衰减倍数和延迟时间参数，物理意义明确，但该方法受限于周期性边界假设条件的限制，只适用于计算平均日导热微分方程，而无法计算全年任意时间的导热微分方程。在本实验中，使用谐波反应法计算平均日绝热保温结构的传热量，可为分析储雪所用绝热保温结构的传热特性提供重要的参考。储雪过程中能量和物质的传输是一个受多要素综合影响的复杂过程，具体的模拟和量化的研究较少［1-2，37-38］。通过对模拟和观测结果的比较证明，谐波反应法可以为建立外界空气温度、太阳辐射、绝热保温结构的热学性质和雪堆融化量的关系提供科学依据，对提前确定储雪量具有重要的参考价值，因此在储雪工程上求解平均日一维非稳态导热问题时，可以使用谐波反应法。

4.1　雪堆融化速率敏感性分析

储雪期间，雪堆1的平均融化量为18.0 kg·d-1（相当于初始质量的0.85%），雪堆2的平均融化量为120.8 kg·d-1（相当于初始质量的6.67%），说明在本次实验中雪堆1上表面所使用的覆盖方案可有效减少外界气象环境对雪堆融化的影响。雪堆的融化速率主要受外界气象条件、绝热保温结构的特性和雪堆几何特征的影响［39］。储雪地点可以选择在遮阳、低风速的位置，以降低外界气象条件对雪堆融化的影响。各材料层的热阻是绝热保温结构热阻的主要贡献者，因此在符合经济性原则的基础上可以选择导热系数较小的材料或增加各材料层的厚度来减少雪堆融化的速率。当外界气象条件一定时，相同体积下，越小的雪堆表面积意味着越少的热量进入雪堆，即雪堆表面积与体积的比例越小，雪堆的融化速率越小。

4.2　误差分析

模拟计算雪堆融化量的误差主要包括以下五部分。

（1）融水全部流出雪堆。雪堆的融化外流过程可以分为两个阶段［40］：雪堆升温为0 ℃时，雪堆表面开始融化并渗入雪堆内部的停蓄阶段；当雪堆内部的液态水达到最大持水量时，融水才能从雪堆中流出的外流阶段。实际融化的雪水并没有全部流出雪堆，有一部分被雪堆截留，保存在雪堆内部。雪堆质量变化反映的是流出雪堆的那部分融水损失。融水的外流时刻和雪堆截留的融水量主要与雪堆的融化速率、最大持水量、雪堆的高度和体积有关。融化速率主要受能量输入大小的影响。雪堆最大持水量的大小主要受雪堆密度的影响［28-29］，储雪过程中，由于压力和变质作用等，导致雪堆密度增大，持水能力降低，促进了融水的外流。雪堆高度的下降将导致融水在雪堆中迁移路径减小，同时高度和体积的减少，使得雪堆容纳融水的空间变小，促进了融水的外流。因此雪堆质量的变化并不能完全代表雪堆的融化量，此误差是由截留在雪堆内部的融水量造成的。雪堆截留水量对融化量模拟的影响在储雪初期较大，导致模拟的雪堆融化量大于雪堆质量的变化。

（2）非稳态方法计算传热量所用绝热保温结构的几何参数。模拟中所使用绝热保温结构的几何参数不随时间的变化而变化，但实际过程中随着雪堆的融化，绝热保温结构的几何参数必然会发生相应变化，此误差将导致模拟的雪堆融化量大于实际的融化量，但在实际的储雪工程中对提前确定储雪量是相对安全的。

（3）雪堆太阳辐射热量。铝箔反射膜与遮阳网之间存在一系列的反射和透射过程［41］，因此在计算雪堆融化量的过程中可能严重低估了太阳辐射对于雪堆融化的影响。

（4）模拟计算所用各材料的基础参数均是通过查表获得，可能与材料实际的参数不符，导致计算结果存在偏差，同时模拟计算结果仅以雪堆质量的变化作为验证标准，缺少雪堆表面热通量等实测数据，使得模拟计算的绝热保温结构的传热量与实际值的偏差未知。

（5）没有考虑长波辐射项。涉及夏季绝热保温领域的研究中，长波辐射项对绝热保温结构外表面的热作用通常被忽略，这是因为在夏季忽略此项对于计算结果是相对安全的［15］，此误差将导致模拟的雪堆融化量大于实际的融化量。

5 结论

通过开展储雪实验，以及定量、定性分析，结果表明雪堆1的平均融化量为18.0 kg·d-1（相当于初始质量的0.85%），雪堆2的平均融化量为120.8 kg·d-1（相当于初始质量的6.67%）。雪堆1上表面绝热保温结构可以有效减少雪堆融化量，证明在新疆阿勒泰地区进行人工储雪是可行的。使用本实验的储雪覆盖方案（反射层和保温层）可以有效减缓雪堆的融化，为今后阿勒泰地区的储雪工程提供技术和方案等方面的借鉴。通过对外界综合空气温度的傅里叶级数描述发现，对于不规则的谐波变化曲线，使用高阶傅里叶级数叠加拟合可以完美描述这一曲线变化。虽然本文对雪堆融化量的计算是一个比较粗糙的过程，但模拟结果与观测值的结果比较仍可以证明，使用数学方法计算雪堆的融化量是可行的。谐波反应法可以为建立外界空气温度、太阳辐射、绝热保温结构的热学性质和雪堆融化量的关系提供科学依据。雪堆覆盖的绝热保温结构可以用反射率、总传热系数、衰减度和延迟时间来进行评价。

虽然雪堆融化量的模拟可以较好地反映雪堆质量的变化，但许多物理过程都未考虑，尤其是未考虑储雪过程中雪堆形状和截留水量的变化，使得模拟结果与实际雪堆融化量的关系仍存在不确定性。下一步储雪实验研究可使用地面激光扫描仪定期获取雪堆随时间变化的高精度几何数据，并在雪堆表面和内部分别布设热通量和含水量传感器，从而进一步验证模拟结果与实际融化量的关系。

参考文献

Grünewald T， Wolfsperger F， Lehning M.

Snow farming： conserving snow over the summer season

［J］. The Cryosphere， 2018， 12（1）： 385-400.

[本文引用: 3]

Lintzén N， Knutsson S.

Snow storage-Modelling， theory and some new research

［J］. Cold Regions Science and Technology， 2018， 153： 45-54.

[本文引用: 4]

Skogsberg K.

Seasonal snow storage for space and process cooling

［D］. Luleå： Luleå University of Technology， 2005.

[本文引用: 1]

Taylor T B.

Ice ponds

［C］//AIP Conference Proceedings. AIP， 1985， 135（1）： 562-575.

[本文引用: 1]

IPCC. Climate change 2013： The physical science basis： working group I.

Contribution to the fifth assessment report of the intergovernmental panel on climate change

［M］. New York： Cambridge University Press， 2013： 187-201.

[本文引用: 1]

Moe J M.

Using stored snow as cooling at Oslo Airport， Norway

［J］. Proceedings of the Institution of Civil Engineers-Civil Engineering， 2018， 171（5）： 11-16.

[本文引用: 1]

Skogsberg K， Nordell B.

The Sundsvall hospital snow storage

［J］. Cold Regions Science and Technology， 2001， 32（1）： 63-70.

[本文引用: 1]

Olefs M， Fischer A.

Comparative study of technical measures to reduce snow and ice ablation in Alpine glacier ski resorts

［J］. Cold Regions Science and Technology， 2008， 52（3）： 371-384.

[本文引用: 2]

Olefs M， Lehning M.

Textile protection of snow and ice： measured and simulated effects on the energy and mass balance

［J］. Cold Regions Science and Technology， 2010， 62（2/3）： 126-141.

[本文引用: 2]

Pestereva N M.

Modern engineering technology to adapt to the adverse weather and climatic conditions at mountain ski resorts

［J］. Life Science Journal， 2014， 11（9）： 800-804.

[本文引用: 1]

Yang Zhanwu.

Research on the snow-making for Beijing 2022 winter olympics and winter Paralympics

［J］. China Winter Sports， 2017， 39（1）： 1-8.

[本文引用: 1]

杨占武.

北京冬奥会和冬残奥会人工造雪的研究

［J］. 冰雪运动， 2017， 39（1）： 1-8.

[本文引用: 1]

Paksoy H Ö. Thermal energy storage for sustainable energy consumption［M］. Dordrecht： Springer Netherlands， 2007：349-353.

[本文引用: 1]

He Jiapeng， Wang Dongfang.

The application study for unsteady heat transmission of wall-method of reaction coefficient

［J］. Journal of Engineering Thermophysics， 2000， 21（1）： 93-96.

[本文引用: 1]

何嘉鹏， 王东方.

对墙体不稳定传热——反应系数法的应用研究

［J］. 工程热物理学报， 2000， 21（1）： 93-96.

[本文引用: 1]

Xu Xiangbo， Hu Yixiong. Air-conditioning loads of building and vehicle［M］. Changsha： National University of Defense Technology Press， 1997： 39-65.

[本文引用: 4]

徐湘波， 胡益雄. 建筑物及汽车空调负荷［M］. 长沙： 国防科技大学出版社， 1997： 39-65.

[本文引用: 4]

Yan Qisen， Zhao Qingzhu. Process of buildings［M］. Beijing： China Architecture & Building Press， 1986： 26-163.

[本文引用: 5]

彦启森， 赵庆珠. 建筑热过程［M］. 北京： 中国建筑工业出版社， 1986： 26-163.

[本文引用: 5]

Tian Qingchao， Wu Jiansheng， Fan Xinmin， et al.

Finite difference technique for heat conduction in multi-layer Gun barrels

［J］. Acta Armamentarii， 2000， 21（4）： 297-300.

田青超， 吴建生， 樊新民， 等.

有限差分法在复合枪管传热中的应用

［J］. 兵工学报， 2000， 21（4）： 297-300.

Ozziek M N. Heat conduction［M］. Beijing： Higher Education Press， 1983： 523-540.

[本文引用: 1]

奥齐西克. 热传导［M］. 北京： 高等教育出版社， 1983： 523-540.

[本文引用: 1]

Kong Xiangqian. Applying of finite element method in heat transfer［M］. 2nd ed. Beijing： Science Press， 1986： 57-142.

[本文引用: 1]

孔祥谦. 有限单元法在传热学中的应用［M］. 2版. 北京： 科学出版社， 1986： 57-142.

[本文引用: 1]

Li Weili， Zhou Feng， Hou Yunpeng， et al.

Calculation of rotor temperature field for hydro-generator as well as the analysis on relevant factors

［J］. Proceedings of the CSEE， 2002， 22（10）： 85-90.

[本文引用: 1]

李伟力， 周封， 侯云鹏， 等.

大型水轮发电机转子温度场的有限元计算及相关因素的分析

［J］. 中国电机工程学报， 2002， 22（10）： 85-90.

[本文引用: 1]

Wang Rongxiao.

Climate change and its influence on main crop planting in Altay region

［D］. Urumqi： Xinjiang Agricultural University， 2014：10.

[本文引用: 1]

王荣晓.

阿勒泰地区气候变化及对主要农作物种植的影响

［D］. 乌鲁木齐： 新疆农业大学， 2014：10.

[本文引用: 1]

Hu Jian.

Reconstruction and analysis of climate based on the width of tree rings in Altay region

［D］. Urumqi： Xinjiang Normal University， 2014：18-20.

[本文引用: 1]

胡建.

基于树轮宽度的阿勒泰地区气候重建与分析

［D］. 乌鲁木齐： 新疆师范大学， 2014：18-20.

[本文引用: 1]

Yusup Abula， Hu Jinming， Tursun Ayshemgul.

Xinjiang Altay skiing

［J］. Shandong Sports Science & Technology， 2015， 37（2）： 36-38.

[本文引用: 1]

阿不拉·玉素甫， 胡金明， 阿依夏木古丽·吐尔逊.

新疆阿勒泰滑雪运动研究

［J］. 山东体育科技， 2015， 37（2）： 36-38.

[本文引用: 1]

Yang Shuping.

Study on the optimization and simulation of cross-season snow storage

［D］. Beijing： University of Chinese Academy of Sciences， 2020：17-19.

[本文引用: 1]

杨淑萍.

跨季储雪的优化方案与模拟研究

［D］. 北京： 中国科学院大学， 2020：17-19.

[本文引用: 1]

Weiss H， Bierman P， Dubief Y， et al.

Optimization of over-summer snow storage at midlatitudes and low elevation

［J］. The Cryosphere， 2019， 13（12）： 3367-3382.

[本文引用: 1]

Yu Jinghua， Yang Changzhi， Tian Liwei， et al.

A study on optimum insulation thicknesses of external walls in hot summer and cold winter zone of China

［J］. Applied Energy， 2009， 86（11）： 2520-2529.

[本文引用: 1]

Anderson E A.

A point energy and mass balance model of a snow cover

［C］// NOAA Technical report NWS 19. U.S. Department of Commence （Washington D.C.

[本文引用: 1]

）， 1976： 150.

[本文引用: 1]

Kondo J， Yamazaki T.

A prediction model for snowmelt， snow surface temperature and freezing depth using a heat balance method

［J］. Journal of Applied Meteorology， 1990， 29（5）： 375-384.

[本文引用: 1]

Wei Zhaocai.

The study on snow melt processes simulation and snow characteristics analysis

［D］. Urumqi： Xinjiang University， 2010：13-46.

[本文引用: 2]

魏召才.

融雪过程模拟及积雪特性分析研究

［D］. 乌鲁木齐： 新疆大学， 2010：13-46.

[本文引用: 2]

Zhang Na.

Research on the effect of snow cover characteristic and radiation energy on snow melting process

［D］. Shenyang： Shenyang Agricultural University， 2017：46.

[本文引用: 2]

张娜.

积雪特性与辐射能量对积雪融化过程的影响研究

［D］. 沈阳： 沈阳农业大学， 2017：46.

[本文引用: 2]

Liu Nianxiong， Qin Youguo. Building thermal environment［M］. Beijing： Tsinghua University Press， 2005： 98-107.

[本文引用: 2]

刘念雄， 秦佑国. 建筑热环境［M］. 北京： 清华大学出版社， 2005： 98-107.

[本文引用: 2]

Hay J E.

Calculation of monthly mean solar radiation for horizontal and inclined surfaces

［J］. Solar Energy， 1979， 23（4）： 301-307.

[本文引用: 1]

Zhang Hongji. Heat conduction［M］. Beijing： Higher Education Press， 1992： 13-25.

[本文引用: 1]

张洪济. 热传导［M］. 北京： 高等教育出版社， 1992： 13-25.

[本文引用: 1]

Lu Yaoqing. Practical heating and air conditioning design manual［M］. 2nd ed. Beijing： China Architecture & Building Press， 2008： 251-254.

[本文引用: 1]

陆耀庆. 实用供热空调设计手册［M］. 2版. 北京： 中国建筑工业出版社， 2008： 251-254.

[本文引用: 1]

Mcquiston F C， Parker J D， Spitler J D. Heating， ventilating， and air conditioning： analysis and design［M］. New York： John Wiley & Sons， Inc， 2004： 129-131.

[本文引用: 1]

Zhao Shuxing.

Optimal tilt angle and installation azimuth of flat plate solar collector for heating

［J］. China Construction Dynamics （Solar & Renewable Energy Sources）， 2005（2）： 56-59.

[本文引用: 1]

赵树兴.

采暖用平板式太阳能集热器的最佳倾角和最佳方位角

［J］. 中国建设动态（阳光能源）， 2005（2）： 56-59.

[本文引用: 1]

Yang Shiming， Tao Wenquan. Heat transfer［M］. 4th ed. Beijing： Higher Education Press， 2006： 44-45.

[本文引用: 1]

杨世铭， 陶文铨. 传热学［M］. 4版. 北京： 高等教育出版社， 2006： 44-45.

[本文引用: 1]

Chen Xiaofei， Tian Jing， Zhang Xueping， et al.

An overview of the researches on the coupled transport of water， heat and solute in the accumulating and melting processes of snowpack

［J］. Journal of Glaciology and Geocryology， 2006， 28（2）： 288-292.

[本文引用: 1]

陈晓飞， 田静， 张雪萍， 等.

积雪融雪过程中水、热、溶质耦合运移规律的研究进展

［J］. 冰川冻土， 2006， 28（2）： 288-292.

[本文引用: 1]

Skogsberg K， Lundberg A.

Wood chips as thermal insulation of snow

［J］. Cold Regions Science and Technology， 2005， 43（3）： 207-218.

[本文引用: 1]

Wang Xing， Qin Dahe， Ren Jiawen， et al.

Numerical estimation of thermal insulation performance of different coverage schemes at three places for snow storage

［J］. Advances in Climate Change Research， 2021， 12（6）： 903-912.

[本文引用: 1]

Shen Bing， Huang Honghu. Principle of hydrology［M］. Beijing： China Water Power Press， 2008： 225-228.

[本文引用: 1]

沈冰， 黄红虎. 水文学原理［M］. 北京： 中国水利水电出版社， 2008： 225-228.

[本文引用: 1]

Shen Mingwei， Hao Feilin.

Investigation of effects of outer shading screens on illumination inside multi-span greenhouse

［J］. Transactions of the Chinese Society of Agricultural Engineering， 2003， 19（6）： 245-247.

[本文引用: 1]

沈明卫， 郝飞麟.

外遮阳对连栋塑料温室内光环境的影响研究

［J］. 农业工程学报， 2003， 19（6）： 245-247.

[本文引用: 1]

/

〈 〉

陇ICP备2021001824号-10
　 甘公网安备 62010202000676号
版权所有 © 《冰川冻土》编辑部
地址：兰州市天水中路8号　
邮编：730000　电话：0931-8260767　
E-mail：edjgg@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司