SWAT 2000气象模拟器的随机模拟原理、验证及改进.pdfVIP

26卷第4期 资 源科学 V01．26。No．4 2004年7月 REsouRcEsscIENcE July，2004 SWAT 2000气象模拟器的随机 模拟原理、验证及改进 张 东1’2，张万昌1 (1．南京大学国际地球系统科学研究所(ESSI)，南京210093； 2．南京师范大学地理科学学院，南京 210097) 摘要：该文详细介绍了swAT 2000气象模拟器(最高、最低气温、太阳辐射)的随机模拟原理，利用甘肃 河西地区酒泉、张掖等7个站的完整气象资料，提出了气候条件相似的区域气象台站月平均气象数据的计 算方法，并计算出模拟器需要的输入参数，对气象模拟结果进行了验证。在此基础上利用逐日平均气温修 正模拟生成的最高、最低气温并求5日滑动平均，对模拟器进行了改进。通过模拟结果与气象台站实测资 料的对比可以看出，改进后最高气温相关系数从0．85提高到o．90，最低气温相关系数从O．89提高到o．94， 气温数据模拟精度有一定提高。最后利用改进后的气象模拟器对河西黑河地区祁连、札马什克、莺落峡3 个数据短缺站点进行了气温和太阳辐射的模拟，通过swAT模型的应用证明模拟结果能够满足分布式水 文模型的气象输入数据精度要求。 关键词：气象模拟器；随机模拟；气温修正；5日滑动平均；分布式水文模型 中圈分类号：s11；s16文献标识码：A 文章编号：1007—7588(2004)04—0028一09 1 引言 近年来，国内外发展了许多综合考虑气候、土壤、植被和水文过程的分布式水文模型n。]，其中国 外比较有代表性的有swAT(Amold等，1995)，VIC(xu 生练等b。提出了基于DEM的分布式流域水文模型。在这些分布式水文模型中，气象数据作为影响水 文循环的重要驱动因素，是模型不可缺少的输入数据，通过它们可以分析水文过程中的地表能量平衡 和水分循环等过程。但是在国内，由于气象站点分布不均匀或者站点气象数据记录不完整，往往得不 到满足模型需要的研究区域气象数据，从而影响了模型的运行和模型效率的发挥。swAT模型…21内 嵌的气象模拟器利用数理统计原理，通过气象资料月平均数据(如气温和太阳辐射的月平均值、标准 差等)的统计规律加上随机噪声矩阵来模拟周期性的气象变化过程，可以实现站点气象资料不全时的 填补或气象数据缺测时的模拟。但是随机噪声矩阵的引入会带来模拟的最高气温低于最低气温的逻 辑错误以及数据产生异常振荡等问题，而且如何在数据资料不全的情况下得到气象模拟器所需要的 月平均气象输入数据也是进行气象模拟前必须解决的问题。本文主要针对解决以上问题以及提高气 象模拟器的模拟精度进行了探讨，以期为今后高精度的气象模拟提供参考。 2气温、太阳辐射随机模拟原理 气温、太阳辐射的变化是一个自然过程，具有季节性和地带性变化的特征。通过对大量气象资料 的数理统计分析可以看出，气温、太阳辐射的周期性变化具有变量与时间呈自相关，变量与变量呈互 收稿日期：2003—08一∞；修订日期：2003—09一Ol 术重点项目资助(编号：2001)。 作者简介：张东(1975一)，男，江苏通州人，博士生，主要从事RS＆GIs与分布式水文模型集成研究。 E·mail：珏柚窖dfr@8ina．c锄． 26卷第4期 张东等：swAT2000气象模拟器的随机模拟原理、验证及改进 29 相关的特点心’63。它们的变化可以看作是一种弱稳定多元随机过程，并随月平均值和标准差的变化而 数据。残差元序列计算方法如下： 尺Si(-，)=A·R|si一。(_，)+B·E；(_，)(-，=1，2，3) (1) 式中：J．依次分别代表最高、最低气温和太阳辐射3个变量序列；船；(，)和麟㈠(J)为考虑了时滞效应 后的第i天和第i一1天3个变量的残差元矩阵；E；(，)为随机噪声矩阵，该随机噪声服从正态分布并