1.2 国内外研究现状与进展
1.2.1 河川径流年内、年际特征分析状况
径流年内各月流量分配的不均匀程度取决于时间和空间的变化情况。径流年内变化情况影响着河流管理中的优化调度和合理配置,也决定了流域内水利工程的建设规模。河流流量随时间变化的基本特征主要从两个方面分析,即年内变化和年际变化。其中,年内变化分析的主要内容是河流流量在一年内各月的分布特征,分析径流量随时间的变化特征。年际变化基本特征主要是对河流历年的流量变化特征进行分析,是进一步分析河流流量变化规律的基础。涂新军等分析了东江流域年内径流量不均匀性分配的年际变异特性,对多年的年内变化特征进行了对比分析;王双银等在冯家山水库河流流量特性研究中从全年各月的流量分配情况,集中度及集中期值计算等方面对流量的年内变化特征进行了分析;李艳等通过计算径流年内分配不均匀性、集中程度、变化幅度等指标对北江流域径流序列年内分配特性进行了分析;胡彩霞等提出了将基尼系数作为分析河流流量年内变化不均匀性特征的指标,以东江流域为例验证了其适用性;徐东霞等在分析嫩江的年内流量的变化情况时,分析的主要指标有流量年内各月的分布情况、每年的完全调节系数、集中度、集中期、年内分配变化幅度等。目前分析河流年内各月流量变化特征的指标主要有变差系数Cv、年内完全调节系数Ct、年内分配变化幅度、集中度、集中期和基尼系数等。年径流变差系数Cv可用于衡量河川径流流量的年际变化,其值反映了河流多年流量的波动变化情况,Cv值越大,河流多年的流量变化越剧烈,不利于流域内水资源的开发利用和防汛抗旱;反之,径流的年际变化平缓。
1.2.2 河川径流序列趋势分析和变异点检验状况
变量的演变主要存在两种形式:一种形式是连续性变化,另一种形式是不连续性变化。径流连续的演变形式是指河川径流的变化趋势,具体指河流的流量随时间的增加呈现连续性的变化,包括持续递减、持续递增和波动情况。马颖等在分析海河水系的径流时间序列的趋势变化时,采用了曼-肯德尔(Mann-Kendall)法;吕继强等采用累计均值法、滑动均值法等分析了新疆和田河流域的时域变化特征;杨帆等基于多种趋势分析方法的理论基础,提出了一种基于插值法的新的趋势分析方法;王生雄等以渭河华县站1956—2000年的多年流量资料为例,采用Mann-Kendall法和中值检验法对径流时间序列进行了趋势分析;刘建梅等在研究杂谷脑河径流变化情况时,引入了小波理论方法。趋势分析的常见方法有滑动平均、线性估计、二次滑动、累计距平和五点三次平滑等。
突变理论的数学基础知识是常微分方程中关于奇点的分析,其关键点是分析某一过程持续稳定状态的快速切换。基于统计学分析,突变的产生是从一个样本统计特征值到另一个样本统计特征值的快速切换。目前,对突变问题的分析还没有形成完善的理论体系,在处理常见的突变现象时,人们根据统计理论、概率理论等提出了一些比较准确的检验方法,这些方法分析数据的3个方面:数据的均值和方差有无突变改变;数据的回归系数是否存在突然改变;由数据资料统计出的事件的概率是否存在突然变化。在突变的研究中,引起广大学者关注的焦点问题是关于该理论在实际生产生活中的应用问题。由于理论的知识以及研究的局限性,到目前为止,突变现象的发生还不能给予非常明确的诠释,正是由于这种问题的存在,可能会得到错误的研究结论。针对该理论存在的问题可以采用以下3种途径进行完善:在检验时间序列的突变问题时,采用多种检验方法来检验,通过各种方法的检验结果对比,提高结果的准确率;对于检验出来的结论,再用多种严格的显著性水平的检验方法进行校核;结合实际,利用专业的水文学知识对突变现象进行判断,以使理论结果与实际现象相一致。在径流序列的突变现象检验中,现在常用的方法主要有:有序聚类法、Mann-Kendall法、山本(Yamamoto)法、滑动t检验、勒帕热(Lepage)法、克拉默(Cramer)法、佩蒂特(Pettitt)法和启发式分割算法(BG算法)等。
1.2.3 河川径流周期分析状况
水文时间序列一般含有概率意义上的周期成分,如何有效地识别、判定水文时间序列的周期成分,已经成为水文学者十分关注的一个研究课题。因为多种因素影响下的水文过程具有不确定性和复杂性,所以,对于水文时间序列的周期分析,学者们分别经历了从傅里叶分析到小波分析,再到希尔伯特-黄变换的研究历程。
(1)以傅里叶变换为基础的传统谱分析。傅里叶变换于1807年由法国数学家傅里叶(Fourier)最早提出,由此出现了与时域相对应的频域;1965年快速傅里叶变换的出现,使频域分析走向实用并迅速拓展。例如,方差谱估计是通过傅里叶变换将水文序列从时间域转为频率域进行周期分析的工具;最大熵谱法(maximum entropy method)是在数据系列熵达到最大的基础上,利用相关函数将假设存在的未知的那一部分数据通过迭代方法推算出来,从而得到功率谱。
奇异谱分析(singular spectrum analysis,SSA)是从时间序列的动力结构出发,并与经验正交函数(empirical orthogonal function,EOF)相联系的一种统计技术,它已广泛应用于时间范围上的信号处理中。SSA的优点主要表现在两方面:①基于自身数据的一种趋势估计方法,不需要其他指标数据,也不存在过度拟合等问题,十分适合于对非线性时间序列的变化进行分析;②对嵌套空间维数的限定,可以有效地对振荡的转换进行时间定位。SSA非常善于识别隐藏在时间序列中的弱信号,是近年来出现的研究周期振荡现象的一种新统计技术,使得提取水文序列周期的技术有了新的飞跃。
(2)小波分析。小波分析(wavelet analysis)由法国地质学家J Morlet最早于1984年提出,核心是小波变换。之后Grossman和Meyer等对小波进行的一系列深入研究使小波分析有了坚实的数学基础。1993年,Kumar等运用正交小波(Harr小波)变换分析了空间降水的尺度和振荡特征,这是小波变换的概念首次应用于水文学领域;自此之后,对于小波变换在水文时间序列分析和预测的方面,国内外的水文学者先后展开了研究。卢文喜等采用Morlet小波分析方法研究了吉林省大安地区年降水量序列的变化特征;翟劭燚等采用Morlet小波研究了海河流域近50年来降水变化的多时间尺度特征;李占玲等采用Morlet小波分析方法研究了雅鲁藏布江流域6个站点的径流序列;傅朝和王毅荣采用EOF和小波变换对黄土高原40年的月降水变化特征进行了研究;李远平和杨太保采用墨西哥帽小波函数研究了柴达木盆地近50年的气温和降水时间序列;邵晓梅等采用墨西哥帽小波函数研究了黄河流域近40年来降水的季节变化和年际变化时间序列;王澄海和崔洋将小波变换和SSA两种分析方法相结合,研究了西北地区26个站近50年的降水周期随时间的变化。
(3)EMD-HHT及EEMD-HHT分析。希尔伯特-黄变换(Hilbert-Huang transform,HHT)是由Huang等于1998年提出来的一种新型的非线性、非平稳的时频分析方法,主要由经验模态分解(empirical mode decomposition,EMD)方法和希尔伯特(Hilbert)变换两部分组成。与傅里叶变换、小波变换等时频分析方法相比,EMD方法具有自适应的优点,因此它能很好地处理非线性、非平稳过程;然而,传统EMD具有模态混叠的问题。对此,Huang曾经提出了中断检验的方法,这种方法的原理是先直接对信号进行EMD,然后再对分解出来的结果进行观察,如果出现模态混叠现象则重新进行分解,可以看出这种方法需要人为后判检验;重庆大学的谭善文提出了多分辨的EMD思想,该方法的原理是对每一个IMF分量都规定一个尺度范围之后再进行EMD,以此来解决模态混叠,可以看出这种方法牺牲了EMD方法良好的自适应性。最后,Huang等又提出了集合经验模态分解(ensemble empirical mode decomposition,EEMD)方法,EEMD方法是在传统EMD方法的基础上进行的改进,它是将白噪声加入信号来补充一些缺失的尺度,在信号分解中具有良好的表现,EEMD方法不仅保留了EMD方法自适应的优点,还有效抑制了其模态混叠的问题。EEMD与Hilbert变换相结合的方法,已经被广泛应用到机电工程、应用力学、机械故障诊断及地震信号检测等领域。
(4)其他分析方法。除了上述的分析方法之外,可以对水文时间序列进行周期分析的方法还有很多,如鲁棒PCA聚类分析法、连续功率谱分析法、R/S分析法、多因子逐步回归分析法以及模糊分析法、灰色系统法、混沌理论分析方法等。其中,主成分分析法(principle component analysis,PCA)对系统具有很高的线性依赖性,为此Xu Lei等提出了在目标函数中加入修正项来提高其鲁棒性的做法,邓红霞将这种改进的鲁棒PCA方法与聚类分析、模式识别等方法相结合,提出了一种特别应用于径流的周期特性研究的鲁棒PCA聚类分析法。
1.2.4 河川径流预测方法状况
国外对于径流预测研究起步较早,研究技术发展速度较快,已初步建立起了较完善的径流预测研究体系。受资料条件和科学技术水平的限制,在很长一段时间里,国内的预测研究一直处在一个与实践经验相关的发展水平上。一直以来,径流中长期预测都是难度比较大的研究课题之一,限于目前的科学技术水平,国内外径流的中长期预测尚无比较成熟有效的方法可以借鉴,该方面的研究仍处于不断探索与发展之中。
中长期径流预测作为水文水资源科学研究的热点和难点之一。国外19世纪末、20世纪初开始对中长期水文预报进行研究。“世界天气法”最先将洪水预报应用于尼罗河下游流域,之后又将其推广到北美和欧洲的一些地区和国家。20世纪30年代,涂长望依据前期气候变化特征对后期长江流域水文情势进行了预报。20世纪50年代,内蒙古水文总站依据“历史演变法”对黄河流域的长期洪水预报进行了研究。张家诚等以东亚大气环流为基础,对长江流域的中长期预测进行了研究。20世纪60年代之后,随着先进探测手段的出现及海洋学、气象学等相关学科的发展,中长期水文预报得到了进一步发展。1992年,美国国家气象中心首先将组合预报的方法引进中长期水文预报领域。近年来,各国不仅注重预报模式分辨率的提高、预报技巧的创新,而且开始强调数值预报产品的综合开发利用。
虽然河川径流预报是一门复杂的非线性、非平稳动力学,但因其具有重要的地位和广阔的应用前景,很多专家学者对此进行了大量的探讨研究工作。当前国内外有关径流中长期预测方面的研究仍处在积极的探索与不断的发展阶段,常用的水文预报方法有基于传统预报技术的数理统计法和成因分析法,以及基于现代人工智能技术的模糊分析法、小波分析法、分形理论、混沌理论、灰色理论、人工神经网络、支持向量机等方法,随着计算机科学的不断发展,预报方法与思路日益增多。根据预报模式的不同,可将水文预报模型粗略地分为数理统计模型和物理成因分析模型两大类。这两大类预测模型各有其适用条件和优缺点。物理成因分析法主要是借助大气环流、天气、海温、太阳黑子等指标进行预报。其主要思路是通过建立后期水文要素(如径流)与前期大气环流特性及与这些特性有关的各种气象因子之间的定量分析关系,来进行水文预报。然而由于水文气象要素长期演变规律的复杂性及其不确定性,其物理机制尚未完全被识别,对物理成因还很难做到客观的定性或定量评价,再加上该方法对数据资料要求较高,因此利用物理成因分析法对径流进行预测难度仍较大,目前物理成因分析法在水文预报部门应用还较少。成因分析法有实际的物理基础作支撑,是今后水文中长期预报一个重要的发展方向之一。数理统计方法是在中长期水文预报中被广泛应用的一种方法,它借助数理统计的理论和方法,从大量的历史资料中识别水文要素自身统计演变的规律性或者寻找预报对象与预报因子之间的内在联系,通过建立预报模型来进行水文预报。数理统计法主要包括单因素预报(时间序列分析法)和多因素预报。时间序列分析法是揭示各种水文现象自身演变规律性的一种有效而可靠的预测方法。作为近代统计学重要分支,多因素分析法在中长期预报中应用领域最多、范围最广。基于现代人工智能的预报技术是当前水文中长期预报研究的热点之一。
传统水文预报技术的逐步成熟和完善,以及理论、方法的不断创新,使上述预报方法的内容和手段不断得到充实,预报方法向多元化、多样化方向发展,从而出现了各种各样的组合预测方法。组合预测方法着眼于单项预测方法的近似性和局限性,通过将多个不同的水文预报方法进行适当组合(如线性加权组合、变化权重组合、熵值法、最优加权组合等),得到组合预报方法,其目的是尽量挖掘各个预测方法所提供的信息,充分发挥各个预测方法的最大优势,从而使预测精度得到提高、预测结果变得更加可靠。组合预报模型如小波分析与人工神经网络组合预报模型、模糊理论与人工神经网络组合预报模型、混沌理论与人工神经网络组合预报模型以及其他组合预报模型,这些组合模型的出现和发展,使径流预测精度得到了一定程度提高,但神经网络与混沌模型在提高模型的泛化能力上缺乏理论依据,在实际应用中受到限制;小波分析法需要事先人为地选择小波基函数以及分解尺度,由于其信息带有明显的主观性,并且这种主观性对分解结果会产生较大的影响,分解结果很难做到对序列自身特点的真实反映,因此,该方法在实际应用中同样有所限制。对于组合模型的研究目前还仅仅处于探索性研究阶段,组合预测的主要问题在于:当前在选择组合预测方法上,尚没有明确的指导原则;对于组合预测模型实际应用的实时性与有效性研究不够。所以,今后关于组合预测研究方面还需投入较大的精力,组合模型的探索研究将成为今后水文预报研究的热点之一。
1976年,Refsgaard和Hansen将AR模型与回归模型进行组合对径流进行预测,结果显示组合预测模型的预测精度明显高于单一模型的预测精度。1983年,Makridakis和Winkler对近10种不同的预测模型进行了研究,并针对不同的组合预测方法给出相应的评价,结果显示在4种组合预测模型中,时变误差平方和方法的预测效果最好,而包含误差协方差的组合方法的预测效果不甚理想。1987年,WcLeod将变换函数噪音模型、自回归模型和概念模型进行组合预测,结果显示组合预测的效果优于单独预测模型的预测效果。1996年,Donaldson和Kamstra将多种基于ANN的组合预测方法引入到经济学研究领域,结果表明组合方法比单项预测方法的预测更精确。1997年,Shamseldin等将ANN组合方法引入到水文学研究领域。1997年,Schreider等将概念性降雨-径流模型与自适应线性筛选方法加以耦合对墨累河的日径流进行了预测。2001年,Shamseldin和O'Commor从线性模型和非线性模型的角度对ANN模型进行了深入的研究,为人工神经网络在水文学中的应用打下了良好的基础。2006年,Kim等对均方差、常系数回归、简单平均和ANN组合预测方法等进行了研究,通过对比各组合模型的预测效果,得出组合方法使预测精度得到了一定程度的提高。2009年,Jeong和Kim将水利经验和解析推导作为参考,对预测组合技术加以选择,在此基础上,从理论和实证角度对组合模型的预测效果做出了评价。
国内对径流预测组合方面的研究开始的比较晚。1996年,黄伟军等借助最优组合手段将门限自回归模型与分级退水模型进行耦合,并利用贝叶斯方法对预测结果做出了评价。2004年,段召辉等将时间序列模型与径流响应线性模型加以耦合对日径流量进行了预测,预测结果显示耦合预测模型对径流预测精度的提高是有效的。2008年,殷峻暹等以丹江口水库为研究对象,首先研究单一时间序列分析模型和ANN模型的预测效果,然后建立二者的组合预测模型并对径流进行预测,结果显示组合预测模型的预测效果明显优于单一预测模型。2009年,傅新忠等通过建立时间序列预测模型(ARIMA)与BP神经网络模型的组合预测模型对径流量进行了预测研究,结果显示组合预测模型具有较强的容错能力,其预测精度也有明显提高,该耦合预测模型对提高径流中长期预测精度是行之有效的。2009年,黄志强等研究了各种水文预测的组合模型,并将其应用到长潭水库流域的洪水预报中,结果显示组合模型的实际应用效果良好。2011年,孙惠子等将差分自回归移动平均模型、人工神经网络和多元线性回归进行简单平均组合及最优加权组合,对枯季径流进行了预测研究,并将组合预报结果与单独预报模型的结果进行了比较,结果显示:与简单平均组合模型相比,通过最优加权法建立的组合预测模型在提高预测精度方面具有明显的优势,并且最优加权组合模型精度不仅取决于各单项预测模型精度,而且与预报精度和各单项模型间的相关性有关。2012年,郭华等将灰色GM(1,1)模型、BP人工神经网络与马尔科夫链进行组合对入库径流量进行预测,结果显示组合模型的预测精度高于单独模型的预测精度。2014年,黄生志等将EMD与SVM进行组合对渭河月径流进行预测,并将预测结果与单一SVM模型的预测效果进行比较,结果证明了EMD与SVM组合模型的可靠性。
国内外对径流中长期预测已做了大量的研究,但对于径流观测数据的处理多数是在时间序列平稳的假定下完成的,这样会导致预测精度大大降低,因此对径流序列的平稳化处理显得尤为重要。1998年,Huang等提出了一种信号处理方法——EMD,本质上该方法是对信号序列进行平稳化处理,本书首先利用EMD或EEMD法对径流序列进行平稳化处理,然后结合各种时间序列预测模型对河川径流序列进行分析预测,提高预测精度,为流域的水资源合理规划与有效管理、水量优化调度等工作提供可靠的水量时空分布信息,进而通过不断提高水资源的利用效率,使有限的水资源社会经济与生态环境功能得到充分的发挥。