第ｌ３卷Ｏ２期 南水北调与水利科技 Ｓｏｕｔｈ－ｔｏ－Ｎｏｒｔｈ Ｗａｔｅｒ Ｔｒａｎｓｆｅｒｓ ａｎｄ Ｗａｔｅｒ Ｓｃｉｅｎｃｅ ８Ｌ Ｔｅｃｈｎｏｌｏｇｙ ＶｏＬ１３ Ｎｏ．０２ ２０１５年６月 ＪｕｒＬ ２０１５ 西安市降水量 的时 序特性分析 乔 亮 ，钱会 ，高盼盼 （１．长安大学环境科学与工程学院，西安７１００５４；２．旱区地下水文与生态效应教育部重点实验室，西安７１００５４） 摘要：利用西安市１９５１年－－２００８年（２０Ｏ６年缺测）的月降水量观测资料，采用线性拟合、云模型及Ｍａｎｎ－Ｋｅｎｄａｌｌ突 变检验３种方法，分析了该地区降水量在丰水期、枯水期及全年三个不同时间尺度的时序特性。通过线性趋势拟合 分析可知，西安市近５７年来降水量总体呈减少趋势，其减少速率为１３．６１３ ｍｒｎ／１０ａｌ用云模型研究降水量分布的不 确定性，结果表明其分布的不均匀性为全年＞丰水期＞枯水期，不均匀性的稳定程度则与之相反；用Ｍａｎｎ－Ｋｅｎｄａｌｌ 法分析降水量的突变情况，发现各时期的变化大都未达到显著性水平，只有丰水期和全年极少时域内存在明显 突变。 关键词：降水量；线性拟合；云模型；Ｍａｎｎ－Ｋｅｎｄａｌｌ检验；西安市 中图分类号：Ｐ４２６．６ 文献标志码：Ａ 文章编号：１６７２—１６８３（２Ｏ１５）ＯＯ２—００６３—０３ Ａｎａｌｙｓｉｓ ｏｎ ｔｈｅ ｓｅｑｕｅｎｃｉａｌ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｐｒｅｃｉｐｉｔａｔｉｏｎ ｉｎ Ｘｉ ａｎ ＱＩＡＯ Ｌｉａｎｇ ．一．ＱＩＡＮ Ｈｕｉ 一．ＧＡＯ Ｐａｎ－ｐａｎ ， （１．Ｓｃｈｏｏｌ ｏｆＥｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｃｈａｎｇ￣ａｎ Ｕｎｉｖｅｒｓｉｔｙ，ＸｉＪａｎ ７１００５４，Ｃｈｉｎａ； ２．Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｓｕｂｓｕｒｆａｃｅ Ｈｙｄｒｏｌｏｇｙ ａｎｄ Ｅｃｏｌｏｇｙ ｉｎ Ａｒｉｄ Ａｒｅａｓ，Ｍｉｎｉｓｔｒｙ ｏｆ Ｅｄｕｃａｔｉｏｎ，Ｘｉ＇ａｎ ７１００５４，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｗｅｔ ｓｅａｓｏｎ，ｄｒｙ ｓｅａｓｏｎ ａｎｄ ａｎｎｕａｌ ｐｒｅｃｉｐｉｔａｔｉｏｎ ｒｅｃｏｒｄｓ，ｆｒｏｍ １９５１ ｔｏ ２００８（１ａｃｋ ｍｅａｓｕｒｅｍｅｎｔ ｏｆ ２００６）ｉｎ Ｘｉ ａｎ ａｒｅ ａｎａｌｙｚｅｄ ｂｙ ｔｈｅ ｌｉｎｅａｒ ｆｉｔｔｉｎｇ，ｃｌｏｕｄ ｍｏｄｅｌ ａｎｄ Ｍａｎｎ－Ｋｅｎｄａｌｌ ａｂｒｕｐｔ ｃｈａｎｇｅ ｔｅｓｔ．Ｔｈｅ ｌｉｎｅａｒ ｆｉｔｔｉｎｇ ｍｅｔｈｏｄ ｉｓ ｕｓｅｄ ｔｏ ａｎａｌｙｚｅ ｔｈｅ ｔｒｅｎｄ ｃｈａｎｇｅｓ．Ｄｅｃｒｅａｓｉｎｇ ｔｅｎｄ— ｅｎｃｙ ｏｆ ｐｒｅｃｉｐｉｔａｔｉｏｎ ｉｓ ｋｎｏｗｎ ｉｎ ｌａｓｔ ５７ ｙｅａｒｓ ｗｉｔｈ １３．６１３ ｍｍ ｐｅｒ ｔｅｎ ｙｅａｒｓ．Ｔｈｅ ｃｌｏｕｄ ｍｏｄｅｌ ｉｓ ｕｓｅｄ ｔｏ ｒｅｓｅａｒｃｈ ｔｈｅ ｕｎｃｅｒｔａｉｎｔｙ ｏｆ ｐｒｅｃｉｐｉｔａｔｉｏｎ ｄｉｓｔｒｉｂｕｔｉｏｎ，ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ａｎｎｕａ１＞ｗｅｔ ｓｅａｓｏｎ＞ｄｒｙ ｓｅａｓｏｎ ｉｎ ｎｏｎ－ｕｎｉｆｏｒｍｉｔｙ ｏｆ ｄｉｓｔｒｉｂｕｔｉｏｎ，ｗｈｉｌｅ ｉｔｓ ｓｔａｂｉｌｉｔｙ ｉｓ ｏｎ ｔｈｅ ｃｏｎｔｒａｒｙ． Ｔｈｅ Ｍａｎｎ－Ｋｅｎｄａｌｌ ｉｓ ｕｓｅｄ ｔｏ ｔｅｓｔ ｔｈｅ ａｂｒｕｐｔ ｃｈａｎｇｅｓ ｏｆ ｐｒｅｃｉｐｉｔａｔｉｏｎ，ｆｉｎｄｉｎｇ ｔｈａｔ ｐｒｅｃｉｐｉｔａｔｉｏｎ ａｂｒｕｐｔ ｃｈａｎｇｅｓ ａｍｏｎｇ ｔｈｅｓｅ ｐｅｒｉｏｄｓ ｄｏ ｎｏｔ ｒｅａｃｈ ｓｉｇｎｉｆｉｃａｎｔ ｌｅｖｅｌ，ｂｕｔ ｏｎｌｙ ｉｎ ｒａｒｅｒ ｔｉｍｅ ｏｆ ｗｅｔ ｓｅａｓｏｎ ａｎｄ ａｎｎｕａｌ ｅｘｉｓｔ ｏｂ￣ｏｕｓ ｃｈａｎｇｅ￣ Ｋｅｙ ｗｏｒｄｓ：ｐｒｅｃｉｐｉｔａｔｉｏｎ ｌｉｎｅａｒ ｆｉｔｔｉｎｇ￣ｃｌｏｕｄ ｍｏｄｅｌ；Ｍａｎｎ－Ｋｅｎｄａｌｌ ｔｅｓｔ；Ｘｉ ａｎ 西安市位于中国腹地黄河流域中部的关中盆地，南 丰水期、枯水期及全年３个不同时间尺度进行分析，其变化 趋势采用线性拟合、不确定性采用云模型、突变检验则采用 Ｍａｎｎ－Ｋｅｎｄａｌｌ法。 靠北秦岭主脊，北跨渭河平原，地貌结构多样，海拔高差位居 全国各市之冠。受暖温带半湿润季风气候影响，其夏秋季高 温多雨，春冬季降水稀少。据统计，该地区２００８年人均占有 地表水资源量不足３８０ ，仅相当于全国人均占有量的六分 之一；地下水资源严重超采，引起地下水位下降、地面下沉， 且近年来水质污染状况愈发严峻，污染正从地表向地下渗 透，这又加剧了水资源的短缺程度ｌ＿１］。大气降水作为水资源 １材料及方法 １．１ 降水资料 西安市１９５１年一２００８年（２００６年缺测）的月降水量观 测数据收集自中国气象科学数据共享服务网。采用统计学 方法得出历年丰、枯水期及全年的降水量作为时序特性分析 的主要补给源，其丰枯程度将直接影响到水资源总量，进而 制约现代社会经济的可持续发展，因此对西安市降水量的时 序特性进行准确分析对于水资源的合理利用及降水的长期 预测具有十分重要的指导意义。 的目标，其中降水的丰枯期按４月一９月为丰水期，１０月至 翌年３月为枯水期划分。 １．２研究方法 已有许多学者对降水量的时序特性做过分析，采用的方 法也多种多样，如线性回归、滑动平均 ］、累积距平、小波分 １．２．１ 线性拟合 也称线性回归 ，是通过研究两个变量问的相关关系 得到某一趋势的数学方法。假设降水量Ｑ与时间ｔ满足一 析＿＿３ ］、Ｍａｎｎ－Ｋｅｎｄａｌｌ检验等。本文主要从西安市降水量在 收稿日期：２０１４—１２—０９ 修回日期：２０１５－０３—２３ 基金项目：水利部公益性行业科研专项经费项目（２０１３０１０８４） 作者简介：乔亮（１９８３一），男，甘肃临夏人，博士，主要从事水资源方面的研究。Ｅ－ｍａｉｌ：ｑｉａｎｇ．＿ｌｉａｎｇ＠ｆｏｘｍａｉｌ．ｃｏｒｎ ｜水｜ｌ丈ｌ采资ｇ濂 ．６３． 第１３卷Ｏ２期・南水北调与水利科技・２０１５年６月 元线性回归方程： Ｑ＝ａｔ＋ｂ （￡一１，２，３，…， ） 式中：参数。为降水量的变化趋势，当ａ＞０时，降水量随时 间变化呈增加趋势，反之则减少，其绝对值的大小反映了此 种变化程度的显著与否。 １．２．２云模型分析 在实际生活中，人们通常习惯于用自然语言来划分、评 判事物所属的类别、等级，如降水的“多”与“少”。云模型 “］ 即可实现这种不确定性语言与定量数值之间的自然转化，它 年份 图ｌ丰、枯水期变化趋势 １０００ ９００ 舶０ ７００ 是在传统模糊数学和概率统计交叉渗透的基础上建立起来 的，通常根据无确定度的逆向正态云发生器（算法１）统计出 降水数据所对应的数字特征（Ｅｘ，Ｅｎ，Ｈｅ）作为云模型的定性 概念（在数据未知情况下可以人为给定），随后采用正向正态 云发生器（算法２）获得一系列云滴并构成云图，从云图便可 直观的反映出降水量随时间变化的不确定性程度。 算法１ 算法２ 输入：样本点ｚ 磊 特征 如， ， 生成云 （Ｅ 输 ： 字特征 输出：，西 ，ＨＰ） 。 个云滴（“” ，… ） ｆ一ｂ ｕｎｃｔ ｉｏｎ ， ａｃｋｃｌｏｕｄ（ｍ］ｆｌ ） ｕｎｃｔｉ。ｎｏｎ ＥＬｘ，ｕ￣ｕｊ＝ｃｌ１ ｄ（ｏｕｄ（ＥｘＥｘ，’ Ｅ—ｎ，’ 胁，’ Ｊｎ） Ｅｘ＝ｍｅａｎ（ｘ）； 求平Ｅｎｎ＝ｎｏｒｍｒｎｄ（Ｅｎ，Ｈｅ＇２，１， ）； 生成以 均值 Ｅｎ为期望，Ｈｅ ２为方差的正态随机数 ｆｏｒ ｋ一１ １ ｃ 一 一， ：” 以Ｅ（ ）一ｎｏｒｍｒｎｄ（Ｅｘ，Ｅ ｎ ２，１，ｎ）；ｘ为期望，西玎２为方差的正态随 薮‘ 生成 Ｅｎ＝ｍｅａｎ（ａｂｓ（ｘ（ｉ）一Ｕ（ ）一ｅｘｐ（一（ｚ（ｋ）一Ｅｘ） ２／（２＊Ｅｎｎ Ｅｘ））＊ｓｑｒｔ（ｐｉ／２）； （ ）一２）））； 计算隶属度 ＥＨ 叩（ ａｒ（ （ｉ） ｎ ２）： ｅｎｄ …“ ｅｎｄ ｐｌｏｔ（ｘ，“，’．’）； 绘制云图 １．２．３ Ｍａｎｎ－Ｋｅｎｄａｌ１突变检验 Ｍａｎｎ－Ｋｅｎｄａｌｌ突变分析法ｌ＿６ 是一种非参数统计检验 方法，其优点在于计算简便，不需要样本遵从一定的分布，且 不受个别异常值的干扰［１３ １５］而受到许多学者的青睐。许多 文献中对该方法均有详细介绍，此处不再赘述。 ２结果与分析 ２．１ 变化趋势的线性分析 分别绘制丰、枯期及全年降水折线图并拟合一元线性回 归曲线（见图１、图２）。由图可以看出：西安市近５７年来降 水量在不同时期的年际变化均呈减少趋势，其中全年降水量 的减少趋势最为显著，减少速率为１３．６１３ ｍｍ／１０ａ，丰水期 略为次之，减少速率为１０．８２４ ｍｍ／ｌＯａ，而枯水期降水量减 少趋势则并不明显（速率２．７８９ ｍｍ／１０ａ）。从降水折线图中 亦可看出，降水量在丰水期和全年的波动幅度较大，在枯水 期的波动则较为平稳，且历年来丰水期降水量均明显大于枯 水期，丰水期降水量的大小及其变化对全年降水量的影响较 为显著。 ２．２分布的不确定性分析 由算法１分别计算各时期降水量云模型的数字特征（见 表】），其中Ｅｘ为多年来降水量大小的平均值，反映多年降 ．６４． 水 丈水 餐｜ 善咖 ５００ ４００ ３００ ２００ １９５１ １９５６ １９６１ ｌ９６６ １９７Ｉ １９７６ １９８ｌ ｌ９８６ １９９１ １９９６ ２００１ ２００７ 年份 图２全年变化趋势 水的总体情况；Ｅｎ是熵，表示降水量大小相对于Ｅｘ的确定 度，亦可反映降水量相对于Ｅｘ的离散度，其值愈大，降水量 分布愈不均匀；Ｈｅ是超熵，即熵的熵，体现降水量离散度的 稳定性。通过算法２绘制出各期对应的隶属云图（见图３、图 ４），从图中可看出：枯水期降水量最少，但分布却最为集中， 稳定性也最高；丰水期降水量虽远大于枯水期，但分布的不 均匀程度及不均匀程度的稳定性较大；全年降水量的分布最 为分散，稳定性也最小。通过云图并结合表１各时期降水的 云模型即可对其在时间尺度上分布的不确定性有更为直观 的认识。 表１ 丰、枯水期及全年降水量的云模型 时期 （Ｅｘ，Ｅｎ，Ｈｅ） 丰水期 （４３３．９６，１０９．２５，１５．９７） 枯水期 （１３８．３７，４１．０９，３．１５） 全年 （５７３．３３，１２０．８７，３９．６９） １ ０．９ ０．８ ０．７ ０．６ 魁 噬０．５ 惴 ０．４ ３ ０．２ ｎ１ ０ ２００ ４００ ６００ ８００ １０（Ｄ ｌ２００ 降水量／ａｒｍ 图３丰水期、枯水期云图 ２．３ Ｍａｎｎ－Ｋｅｎｄａｌｌ突变分析 利用Ｍａｔｌａｂ ７．１软件平台编程分别计算降水量在丰、枯 水期及全年的统计量ＵＦ 、ＵＢ ，同时获得检验图，见图５。 从图５看，降水量在１９５１年一１９６１年间呈增加趋势，且 于１９５８年达到显著性，１９６１年后则呈减少趋势；降水量在 １９６２年和１９６５年前后发生突变并以１９６５年的突变较为显 乔 亮等・西安市降水量的时序特性分析 图４全年云图 著。图６中，ＵＦ、ＵＢ曲线均处于临界区域内，说明降水在枯 水期的变化并不明显，但在１９５２年前后、１９５５年、１９５７年、 １９６７年、１９６８年、１９７１年末、２００３年、２００４年均有突变现象， 其中１９７１年末发生的突变最为显著；降水量在１９５３年一 １９５４年、１９６１年一１９６６年及１９６８年左右呈增加趋势，其余 年份呈则减少。图７中，降水量在１９５１年一１９６２年、１９６４ 年一１９６５年两个时间段呈增加趋势，其余年份表现为减少趋 势；于１９６６年开始突变（降水量由多变少），在１９７９年至 １９８０年末、１９８２、１９９５、１９９７年各年前后及２００１年一２００３年 间突变均达到显著性水平。综合图５一图７，各时期大部分 的统计量位于９５ 置信区间［一１．９６，１．９６］内，且ＵＦ曲线 也大多位于零值以下，说明西安市１９５１年一２００８年的降水 量整体呈减少趋势，这与线性拟合分析的结论一致。 ４ ——ＵＦ统计量 …・ ３ ． 昀 计量 ……－０ ０５显著水平 ２ 懈ｔ 七 ／／ 螺０ 一Ｉ ／‘ ：也 ＿３ １９５５ ｌ９５０ １９６５ １９７０ ｌ９７５ １９８０ １９８５ １９９０ １９９５ ２０Ｏ０ ２００５ 年份 图５丰水期突变检验 ——ＵＦ统计量 …・． ， “蹴计量 ……－０ ０５显著水平 姆 Ａ： 年份 图６枯水期突变检验 ３结论 线性拟合、云模型及Ｍａｎｎ－Ｋｅｎｄａｌｌ突变检验三种方法 一ＵＦ统计量 …・． “蹴计量 ……一００５显著水平 咖 士 蠕 ． Ｖ… Ｖ一 年份 图７全年突变检验 均能较好的反映西安市１９５１年一２００８年（２００６年缺测）降 水量的时序特性，但仅通过其中一种难以做出全面的分析。 综合上述三种方法可获得以下结论。 （１）西安市近５７年来降水量总体呈减少趋势，其中全年 表现最为显著，减少速率为１３．６１３ ｍｍ／１０ａ，丰水期略为次 之，为１０．８２４ ｍｍ／１０ａ，枯水期为２．７８９ ｍｍ／１０ａ。 （２）降水量的不确定性表现为分布的不均匀性及不均匀 性的稳定程度，其全年降水量的分布最为分散，稳定程度最 低，其次为丰水期，枯水期则相对集中，稳定性也高。 （３）降水量在各时期的变化大都未达到显著性水平，只 有丰水期和全年极少时域内发生显著突变，枯水期则突变点 较多，但变化幅度均较小，呈上下波动状。 （４）从降水量的线性趋势、不确定性及突变情况看，丰水 期和全年的变化状况均较相似，说明丰水期降水量的大小对 全年降水量的影响较为显著。 参考文献： ［１］樊春贤．西安市水资源可持续利用的探讨ＥＪ］．地下水，２００７，２９ （１）：１７—１９． ［２］吴建华，李培月，钱会．西安市气象要素变化特征及可利用降雨 量预测模型口］．南水北调与水利科技，２０１３，１１（１）：５０—５４． ［３］张允，赵景波．西安市近５５年来降水的多时间尺度分析［Ｊ］．中 国农业气象，２００８，２９（４）：４０６—４１０． ［４］易淑珍，王钊．水文时间系列周期分析方法探讨［Ｊ］．水文， ２００５，２５（４）：２６—２８． ［５］衡彤，王文圣，丁晶．降水量时间序列变化的小波特征［Ｊ］．长江 流域资源与环境，２００２，１１（５）：１０３—１０５． Ｅ６］赵华荣，夏北成，江学顶．桂林市近５０年降雨变化特征分ＥＪ］． 桂林理工大学学报，２０１１，３１（２）：２０２—２０６． Ｅ７］王颖华，张鑫．西宁市降水量特征及变化趋势分析＿Ｊ］．水土保 持研究，２０１１，１８（５）：１５６—１６０． ｒ８］ 陈敏．三门峡市近５Ｏ年降水变化趋势及周期分析［Ｊ］．气象与 环境科学，２０１３，３６（１）：４４—４９． ［９］李德毅，杜鹚．不确定性人工智能［Ｍ］．北京：国防工业出版社， ２００５． ｒ１Ｏ］高文荣，段春青，陈晓楠，等．基于云推理的年降水预测模型 ［Ｊ］．中国农村水利水电，２００９，（３）：３２—３５． ［１１］尹国定，卫红．云计算一实现概念计算的方法［Ｊ］．东南大学学 报（自然科学版），２００３，３３（４）：５０２—５０６． Ｅ１２３ Ｋｅｎｄａｌｌ，眦Ｇ．Ｒａｎｋ Ｃｏｒｒｅｌａｔｉｏｎ Ｍｅｔｈｏｄｓ，４ｔｈ ｅｄ．Ｃｈａｒｌｅｓ Ｇｒｉｆｉｎ，Ｌｏｎｄｏｎ，１９７５． （下转第９６页） ｌ永 采 蒜｜穗辫 ・ ・ 第１３卷０２期・南水北调与水利科技・２０１５年６月 行对比，结果见图９。可见在２０１１年下半年，Ｂ７孔１５ ｎｌ深 验，初步获得了窟顶覆盖层以及山体砂岩包气带不同部位的 土壤湿度、气压、气温和空气湿度数据，并对其时空变化特征 的成因进行了讨论。对监测数据的分析表明石窟山体包气 带空气相对湿度则接近饱和状态。随深度增加温度变幅相 对减小，埋深１ ＩＴＩ以下的砂砾层气温稳定在１１．６℃～１１．７ 处气压与大气压都经历了总体抬升的过程，但包气带气压抬 升的幅度低于大气压，夏季大气压低于包气带，而冬季则相 反。１ ５ ｒｎ深处温度约为９．５℃，几乎保持不变。有效气压 差与气温的变化关系见图１０，容易看出在９月份气温开始快 速下降，而１５ ｎｌ深处包气带的气压也从高于大气压逐渐变 为低于大气压。这说明，云冈石窟包气带夏季存在向上的气 压梯度，驱动空气和湿气向上运移，而冬季存在向下的气压 梯度，驱动空气和湿气下沉。这种随空气流动产生的湿气运 移对石窟文物有什么影响，还需要更多的调查研究。 ℃左右，反应了云冈地区的多年平均气温。气压观测发现包 气带存在垂向上的气压梯度，夏季驱动空气和湿气向上运 移，而冬季驱动空气和湿气下沉，这样的气压梯度为空气和 湿气的运移创造了客观条件。 参考文献： ［１］黄继忠．世界文化遗产云冈石窟的防水保护ＦＪ－１．文物保护与考 古科学，２００８（ｓ１）：１１４—１２１． 占 壬 目 ３ ［２］甘向群．云冈石窟凝结水研究与防治［Ｄ］．中国地质大学（ｊ 京），２００６． 趔 Ｉ－３］尚剑宝，上官学兵．山西大同云冈石窟防水保护方案探讨ｒＪ］． 山西建筑，２００３，６（２９）：７５￣７６． Ｅ４］张红梅，马国栋，速宝玉，等．大同云冈石窟文物渗水病害防治 方案探讨［Ｊ］．水文地质工程质，２００４（５）：６４—６７． ［５］万力，曹文炳，王旭升，等．云冈石窟水汽转化特征的初步研究 ［Ｊ］．工程勘察，２０１２（１１）：６－１１． ［６］Ｃｈａｒｏｌａ Ａ Ｅ．Ｓａｌｔｓ ｉｎ ｔｈｅ ｄｅｔｅｒｉｏｒａｔｉｏｎ ｏｆ ｐｏｒｏｕｓ ｍａｔｅｒｉａｌｓ：ａｎ ｏｖｅｒｖｉｅｗ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ａｍｅｒｉｃａｎ Ｉｎｓｔｉｔｕｔｅ ｆｏｒ Ｃｏｎｓｅｒｖａｔｉｏｎ， ２０００，３９（３）：３２７—３４３． ｒ７］Ａｍｏｒｏｓｏ Ｇ Ｇ，Ｆａｓｓｉｎａ Ｓｔｏｎｅ ｄｅｃａｙ ａｎｄ ｃｏｎｓｅｒｖａｔｉｏｎ［Ｍ］． Ａｍｓｔｅｒｄａｍ ａｎｄ Ｎｅｗ Ｙｏｒｋ：Ｅｌｓｅｖｉｅｒ，１９８３． ［８］傅志斌，彭涛．云冈石窟防水保护工程勘察报告ＥＲ］．建设综合 勘察研究设计院，２００３． ［９］黄继忠，张俊芳．粉尘对云冈石窟石雕影响的研究［Ｊ］．文物保 护与考古科学，２００４，１（１６）：１—８． ［１Ｏ］董佩，王旭升，万力，等．云冈石窟覆盖层特征及其对石窟渗水 的影响口］．科技导报，２０１３，３１（１Ｏ）：６２—６５． ［１１］王旭升，万力，彭涛，等．云冈石窟人渗水的形成和运移［Ｊ］．工 程勘察，２０１２（１１）：１２—１５． 图ｉ０ ２０１１下半年Ｂ７子Ｌ气压差和温度变化曲线 （有效气压差是用大气压减去１５ ｍ深度 ［１２］曹文炳，万力，周训，等．西北地区沙丘凝结水形成机制及其对 生态环境影响初步探讨［Ｊ］．水文地质工程地质，２００３（２）：６－１０． ［１３］万力，曹文炳，周训，等．包气带中温度变化对水分分布影响的 实验研究［Ｊ］．水文地质工程地质，２００４（３）：２５—２８． ［１４３吴宝燕，张爱国，朱丽君．云冈石窟水害及其治理措施［Ｊ］．地 下水，２００８，５（３０）：１１６一】１８． 补偿气压得到的）［５］ ４结论 本文特别 测‘云『Ｘ］石窟围岩和覆盖层所在包气带的湿 热 子进行了覆盖层监ＮＳＬ、砂岩垂直和水平监测孔观测试 ［１５］黄继忠，万力，彭涛，等．云冈石窟水分来源探查工程及若干成 果［Ｊ］．工程勘察，２０１２（１１）：１－４． （上接第６５页） ｒ１ ３］Ｚｈａｎｇ Ｑ．Ｘｕ Ｃ Ｙ，Ｚｈａｎｇ ｚ， ａＩ．ｓｐａｔｉａｌ ａｎｄ ｔｅｍｐｏｒａｌ ｖａｒｉａ ｂｉｌｉｔｙ ｏｆ ｐｒｅｃｉｐｉｔａｔｉｏｎ ｏｖｅｒ Ｃｈｉｎａ，１９５１ ２００５［Ｊ］．Ｔｈｅｏｒｅｔｉｃａｌ ａｎｄ Ａｐｐｌｉｅｄ Ｃｌｉｍａｔｏｌｏｇｙ，２００９，９５（１ ２）：５３—６８． ［ｊ］．Ｔｈｅｏｒｅｔｉｃａｌ ａｎｄ Ａｐｐｌｉｅｄ Ｃｌｉｍａｔｏｌｏｇｙ，２００８，９３（１－２）：３５— ４３． ［１５］Ｚｈａｎｇ Ｙ Ｊ，Ｋａｎｇ Ｓ Ｃ，Ｑｉｎ Ｄ Ｈ，ｅｔ ａ１．Ｃｈａｎｇｅｓ ｉｎ ａｎｎｕ— ａｌ ａｃｃｕｍｕｌａｔｉｏｎ ｒｅｔｒｉｅｖｅｄ ｆｒｏｍ Ｇｅｌａ ｄａｉｎｄｏｎｇ ｉｃｅ ｃｏｒｅ ａｎｄ ｉｔｓ ｒｅｌａｔｉｏｎｓｈｉｐ ｔｏ ａｔｍｏｓｐｈｅｒｉｃ ｃｉｒｃｕｌａｔｉｏｎ ｏｖｅｒ ｔｈｅ Ｔｉｂｅｔａｎ Ｐｌａｔ— ［１４］ Ｚｈａｎｇ Ｑ，Ｘｕ Ｃ Ｙ，Ｚｈａｎｇ Ｚ，ｅｔ ａ１．Ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｏｒ ｖａｒｉａｂｉｌｉｔｙ．Ｔｈｅ ｃａｓｅ ｏｆ Ｙｅｌｌｏｗ ｒｉｖｅｒ ａｓ ｉｎｄｉｃａｔｅｄ ｂｙ ｅｘｔｒｅｍｅ ｍａｘｉｍｕｍ ａｎｄ ｍｉｎｉｎｍｍ ａｉｒ ｔｅｍｐｅｒａｔｕｒｅ ｄｕｒｉｎｇ １９６０—２００４ ｅａｕ［Ｊ］．Ｃｈｉｎｅｓｅ Ｓｃｉｅｎｃｅ Ｂｕｌｌｅｔｉｎ，２００７，５２（２３）：３２６１—３２６６． ・９６・ 生态与环境
