不同气象信息插值算法对SWAT模型模拟影响研究.pdfVIP

不同气象信息插值算法对SWAT 模型模拟的影响研究 雷晓辉 王海潮 陈宁 蒋云钟 (中国水利水电科学研究院水资源研究所 北京 100044) 【摘 要】 本文对影响SWAT模型模拟精度的气象信息插值算法进行研究，认为其现有方法比较 适合气象数据缺失的地区。通过雨量／气象插值方法的研究，选取4种有代表性的降水量插值方法 和1种温度插值方法对雨量／气象站较多，观测系列较长的地区进行雨量／气象插值。以北京市为 例，建立SWAT模型，并用多线程并行计算对1956～2000年系列雨量／气象信息进行插值，修改 一——《≈===童◆葛==∞—～—————— 了利用统计方法计算已知信息，较为精确地模拟了雨量／气象信息，提高了SWAT模型模拟精度。 【关键词】 降水量插值 温度插值 SWAT l SwAT模型现有方法 一致：各个子流域寻找距离每个子流域中心最近的雨量站，然后用该雨量站信息来代表该子 流域的雨量信息。降水量文件(．pcp)的存储所有雨量站的历史系列，其中每一列只存储一 个雨量站的降水过程。每个子流域参数文件(．sub)只存储当前子流域对应降雨量文件中的 某列降水量过程，即．pcp文件的哪一列。气象信息的处理方式和降水量基本一样。其中 MWSWAT为了进一步减少计算时间，无论用户提供多少雨量站，只是在一个指定的雨量 站数量(例如30)中进行搜索。另外，这种方法只是针对雨量／气象站来选择，而很多情况 下有些雨量／气象站的记录会出现缺测或者丢失等，SWAT模型通过气象发生器来自动弥补 这部分缺失的数据。这种处理方式比较适合气象资料缺乏地区，这种情况下其他插值方法不 好使用，只有选择距离每个子流域最近的水文站／气象站来作为代表站。而事实上这些缺失 的雨量／气象站的信息完全可以通过周边的雨量／气象站信息插值得到。 SWAT模型中可以对降水和温度进行高程修正，把子流域划分为多个等高带(小于10个)， 然后根据该子流域中某个等高带与该子流域的代表站的高程，通过高程修正分别计算每个等高带 的雨量／气象信息，然后把各个等高带的信息进行面积加权平均得到整个子流域的信息。但是修 正过程需要修改每个子流域的．sub文件，手动添加等高带等等，实现过程比较繁琐。 2雨量／气象插值方法介绍 区域降水量是水文模拟中的重要参数之一，对于研究区域流域水资源合理规划与利用具 有重要科学价值。目前，降水量都是来自固定气象站点的独立观测资料，但是，当研究扩展 到景观和区域尺度时，仅仅基于站点的观测资料远远不能满足要求，降水的连续性空间分布数 据更为重要。因此，对降水量进行插值和模拟的研究越来越多L1剖，由于降水通常随海拔高度 升高而升高[5]，因而大量早期研究都是通过一定的插值算法或趋势分析模型，将样点的地理位 置和海拔高度作为空间变量，利用空间变量与降水的关系进行插值并生成降水空间分布图。本 文选用泰森多边形、Kriging、距离反比法、修正距离反比法计算研究区降水量分布。 气温是由太阳辐射到达大气形成的，在地球表面的分布与地表的地形特征密不可分，特 别是在山区，地形复杂多样，影响程度不一，温度的空间变化很大。关于山区平均温度随海 22 中国水利学会第四届青年科技论坛论文集 拔高度的分布问题，许多人进行过详细的研究睁¨。山区气温的推算大致分为两类：～类是 分解综合法，另一类是函数相关法(主要包括统计模拟法、回归余项法和物候法等)。本文 选用多气象站线性插值法计算研究区气温分布。 2．1泰森多边形法 泰森多边形法又称最近距离法，是降水资料空间展布和计算流域面雨量时最常用的方 法。它的基本思想十分简单，无观测地点的雨量等于离其最近观测站点的雨量。泰森多边形 法的插值步骤包括： (1)定义一组控制站点。 (2)相邻三个站点间连接成三角形。 (3)找出所有控制站点的外边界凸体多边形。 (4)对每个三角形的每条边画垂直平分线。 (5)这些垂直平分线构成的多边形便是泰森多边形。 泰森多边形内所有位置的雨量值等于该多边形内雨量站点的雨量值。 2．2 Kriging法 Kriging方法是一种特定的滑动平均加权法。这种方法假定采样点问的距离和方向可反 映一定的空间关联，并利用它们来解释空间的变异性。它主要是利用一定的数学函数对待估 点或者是给定搜索半径内的所有点