黄河流域参考作物蒸散发量的混沌特性研究

　　摘要：为了更好的研究流域参考作物蒸散发量在时间上和空间上的变化特征，并合理利用参数来预测未来参考作物蒸散发量的变化规律，本文以黄河流域内玛曲站为代表，通过混沌理论和相空间重构技术，在重构相空间的基础上分别计算出饱和关联维数D2、最大Lyapunov指数来表征系统的混沌性质，为进一步预测黄河流域参考作物蒸散量的未来变化趋势提供依据，也为缺乏气象资料的地区预测参考作物的变化提供新的思路。

　　关键词：黄河流域；蒸散发量；混沌特性

　　蒸散发量是水文循环的重要环节，准确的预算流域蒸散发量能深入了解流域内陆—气之间水热交换和水分循环过程，对于流域水资源规划管理和利用极为重要。通过了解流域参考作物蒸散发量的非线性变化特征及其特征的变化规律，对预测黄河流域参考作物蒸散发量未来几年可能的非线性变化、指导黄河流域绿色农业参考作物种植区域生态结构的变化进行科学良性调整、合理开发和利用可调配的种植地区水土资源、促进种植区域绿色农业和种植地区经济与生态环境的良性平衡发展均有重要意义[1]。

　　1研究背景

　　近年来，在我国的水文科学中已经开始了广泛应用混沌理论。在我国，丁晶最先深入研究混沌理论，他论述了将混沌运动特性分析的理论和方法，并把它广泛应用于暴雨和洪水的时间序列，证明了在洪水过程中的混沌现象[2]。本次研究以黄河流域为研究对象，利用相关气象数据，通过Penman-Monteith公式计算出参考作物蒸散发量，进而分析得出黄河流域参考作物蒸散发量的时间和空间变化特征，最后运用量化方法描述出不同时空黄河流域参考作物蒸散发量的混沌特性。

　　2研究区域概况

　　黄河流域总面积752443km2，涉及9个省（区），黄河流域自然和生态条件复杂多样，东西走向的跨越构成我国地形的三大阶梯，由于区域内地形以及气候条件的复杂多样性，造成了水资源的时空分布不均匀，洪涝灾害频发，水资源的供需不平衡，人和自然的矛盾日益突出，严重影响了流域内经济的可持续发展。

　　3研究的方法

　　3.1Penman-Monteith公式法

　　采用FAO1998年推荐的Penman-Monteith公式法计算逐日参考作物蒸散发量[3]：ETO=0.408（Rn-G）+γ900T+273μ（2ls-la）△+γ（1+0.34μ2）（1）式中，ET0参考作物蒸散发量为逐日参考作物蒸散发量（mm/d）；Rn为冠层表面净辐射[MJ/（m2·d）]；G为土壤热通量[MJ/（m2·d）]，一般情况下为0；T为平均温度（℃）；Δ为饱和水汽压温度曲线的斜率（KPa/℃）；μ2为高度2.0m处风速（m/s）；es为饱和水汽压（KPa）；ea为实际水汽压（KPa）；γ为湿度常数（KPa/℃）。

　　3.2相空间重构理论

　　一般时间序列主要是在时间域中进行模型研究，混沌时间序列必须在相空间中模型研究，因此说相空间重构是分析时间序列的关键环节[4]。

　　4实例应用及结果分析

　　4.1参考作物蒸散发量的时间序列变化特征

　　利用相关气象数据通过Penman-Monteith公式计算出参考作物蒸散发量，以此作为混沌识别的基础序列。以玛曲站为代表站绘制时序图，从图1可以看出，参考作物蒸散发量的时间序列的过程图表现出强烈的季节周期性和波动性，但是仅从图像并不能判断出此序列是否为混沌序列，需要借助于其他方法来判断参考作物蒸散发量的时间序列中是否存在混沌现象[6]。

　　4.2参考作物蒸散发量混沌特性识别

　　由于混沌理论在水文领域的发展历程较短，没有任何一位学者明确提出某一种混沌特性识别的方法，因此使用单一的方法进行时间序列的识别可靠性不高[7]。本文对序列的识别主要采用的是定量分析中的饱和关联维数法和最大Lyapunov指数法进行对比识别。

　　4.2.1饱和关联维数法混沌特性识别由公式（2），计算出黄河流域玛曲站点1970—2009年的参考作物蒸散发量年月均值的自相关函数值，并且作出自相关函数值随延迟时间变化的图像，选择自相关函数值第一次经过零点的时间作为延迟时间。如图2所示。

　　4.2.2最大Lyapunov指数法混沌特性识别

　　由于本文数据量很小，所以选择了易于操作且计算精度相对较高的小数据量法。首先利用快速傅里叶变换估算出时间序列的平均周期P，取最高谱线所对应的下标求出平均周期。运用MATLAB中的相关程序计算最大Lyapunov指数，计算结果见图5。最后利用最小二乘法拟合得到的直线的斜率就是参考作物蒸散发量时间序列的最大Lyapunov指数λ。玛曲站的最大Lyapunov指数λ为0.007。

　　4.3参考作物蒸散发量混沌特性分析

　　黄河流域玛曲站点的混沌特性的分析方法采取以下两种方法，一是饱和关联维数法，二是最大Lyapunov指数法，计算成果如表1所示。

　　分析表1可知，黄河流域玛曲站均存在饱和关联维数且都为分数维，并且有最大Lyapunov指数λ＞0，充分说明黄河流域参考作物蒸散发量时间序列存在混沌现象。

　　5结论

　　本文选取了黄河流域内的玛曲气象站点1970—2009年的逐日气象资料，采用Penman-Monteith公式计算出参考作物蒸散发量的时间序列。对参考作物蒸散发量时间序列的混沌识别，主要在相空间重构的基础上，通过饱和关联维数法和最大Lyapunov指数法两种量化方法，得出了时间序列的关联维数随着嵌入维数的增加而达到饱和状态且为分数维，并且均有最大Lyapunov指数λ＞0，计算结果均表明了黄河流域参考作物蒸散发量时间序列存在混沌特征。这不仅可以了解黄河流域参考作物蒸散发量的非线性变化规律，而且为预测流域内参考作物蒸散发量的未来变化趋势提供依据。

　　参考文献

　　[1]杜加强,舒俭民,刘成程,等.黄河上游参考作物蒸散量变化特征及其对气候变化的响应[J].农业工程学报,2012(12):92-100.

　　[2]傅军,丁晶,邓育仁.洪水混沌特性初步研究[J].水科学进展,1996(3):43-47.

　　[3]张艳兰,栾元重,尹燕运,等.混沌时间序列相空间重构及特性识别[J].测绘科学,2016(4):15-18.

　　[4]马建华,楚纯洁.黄河流域动力系统泥沙时序混沌特征分析——地理系统综合研究的一种尝试[J].地理研究,2006(6):949-958.

　　[5]战国隆,马孝义,刘继龙,等.蒸发蒸腾量时间序列的混沌特征分析[J].农村水利水电,2010(3):22-24+28.

　　[6]李国良,付强,冯艳,等.混沌理论及其在水文水资源中的应用[J].农业系统科学与综合研究,2006(3):173-176.

　　[7]刘道文,忽海娜.非线性时间序列混沌特性分析及短期预测[J].计算机仿真,2012(10):370-373.

　　[8]陈引锋,马长铃.黑河日径流量混沌变化特性的研究——Ⅲ最大李雅普诺夫(Lyapunov)指数的确定[J].地下水,2011(3):136-138.

　　张勇

《黄河流域参考作物蒸散发量的混沌特性研究》

上一篇：刍议质量管理在食品行业的应用

下一篇：浅谈交通运输设施与农村经济发展