参考作物腾发量时间序列的长程相关性和多重分形分布 Long-term correlation and multi-fractality of reference crop evapotranspiration time series.pdfVIP

水 利 学 报 XUEBAO 2008年12月 SHUILI 第39卷第12期 文章编号：0559．9350(2008)12-1327—07 参考作物腾发量时间序列的长程相关性和多重分形分布 谢先红，崔远来，周玉桃 (武汉大学水资源与水电工程科学国家重点实验室，湖北武汉430072) 摘要：本文应用修正的Penman-Monteith模型计算了分布在全国各地的16个站点的参考作物腾发量时间序列，并采 用新近发展的多重分形非趋势波动分析法(MF-DFA)计算了序列的长程相关指数日和奇异谱f(a)。结果表明，潜 在蒸散序列围绕一定趋势的波动不是完全随机的现象，而是由内在自相似机制决定的长程相关过程；其不同层次 的波动是一个有序的多重分形分布，奇异谱可用二项倍增串级模型描述。 关键词：参考作物腾发量；波动；长程相关性；非趋势波动分析；多重分形 中图分类号：S161．4 文献标识码：A 1研究背景 蒸发蒸腾(亦称腾发或蒸散)是水文循环的重要环节，且在水循环中最直接受土地利用和气候变化 的影响¨】。参考作物腾发量是在供水充足，一定气象条件下植被地带最大的蒸发蒸腾量，常常是计算实 际腾发量的基础。因此，对参考作物腾发序列波动过程的研究，用以预报、模拟水文变量，乃至进一步认 识气候变化过程引起了广泛的关注。气候和水文变量，如气温、径流等序列具有多种周期，且波动不是 完全随机的，而是长程幂率(分形)相关的，相关性可以延续至几年乃至几十年心J1。蒸发蒸腾过程表现 出相似特征，随季节、年和年代际而表现出周期变化和非线性过程。杨辉等H1利用小波理论分析指出华 o发 北地区参考作物蒸发(蒸腾)的年代际周期主要是15年，而年际周期在不同年代有所不同；任玉敦b 现参考作物腾发量存在较为明显的5年、8一lO年及20年变化周期。刘丙军等m，1进一步发现了参考作 物腾发量序列在小尺度和大尺度上存在局部和整体的相似性，即具有分形特征。但是蒸发蒸腾的长程 (不同尺度)波动是否与其它气候水文变量一样具有长程相关性，而且其不同层次的波动是否可用单分 形结构完全描述，还需进一步研究。 就分析方法来说，传统的功率谱分析、相关分析等只适合研究平稳信号，并且单分形只用一个标度 指数也往往无法描述非平稳且含有噪声的时间序列不同层次波动性。因而，一些数学方法被提出以克 DetrendedFluctuation 服这些缺陷，其中结合了多重分形理论的非趋势波动分析法(Muhffractal Analysis， 简写MF—DFA)被认为是行之有效的方法，广泛应用于气候、水文序列研究中心J】。从动力学角度讲，该 方法滤除了原始序列中未知的趋势成份，剩下的离差序列保留着原始序列的波动成分，可避免长程相关 性的误判㈨。 本文应用MF．DFA方法探讨参考作物腾发量的长程相关性，并计算奇异谱来描述其不同层次的波 动。研究结果将有益于了解全球气候变化下蒸散过程的动力学机制和不同尺度的变化特征，及其未来 的变化趋势；同时为建立综合的气候、水文模型或灌区作物需水的预测提供新的思路。 收稿日期：2007-09．19 一1327— 万方数据 2资料和方法 2．1基础数据所采用数据由国家气象局资料中心提供，包括各站点的经纬度以及最低和最高气温、 风速、气压、相对湿度、日照时数的日值，时间范围从起始年(见表1)的1月1日至2005年12月31日。 个别缺测数据或异常数据采用前3日的平均值代替。经统计检验，所选站点数据可靠，适合参考作物腾 发量计算分析。根据所采用16个站点的多年平均降水和参考作物腾发量，计算得各站干燥度，由此评 判各站所在气候区¨0|。可见，它们分别处于湿润、半湿润、半干旱、干旱、极干旱等5个气候区，从而采 用数据具有一定