径流的混沌特性分析及非线性预测方法技术

本发明专利技术公开了径流的混沌特性分析及非线性预测方法，采集的径流数据大于39个；对每个径流数据进行小波变换，生成若干个子序列；对每个子序列采用小数量法或其他方法分别计算最大lyapunov指数；对最大lyapunov指数为正数的子序列判定为混沌时间序列，分别计算混沌时间序列的延迟时间和嵌入维数；根据计算出来的延迟时间和嵌入维数，对相应的具有混沌特性的子序列进行相空间重构，将重构结果作为预测的数据输入神经网络模型完成最终预测。本发明专利技术的有益效果是江河海域的径流记录数据大于39个时，即可进行混沌特性分析及非线性预测。

【技术实现步骤摘要】

本专利技术属于水利
，设及径流的混浊特性分析及非线性预测方法。
技术介绍
非线性现象是自然界的普遍现象，非线性系统具有超越不同学科领域的性质（相 似性）。该种复杂系统的相似性，并非在于定性描述，而在于实际测验或理论细节，如孤立 子、相干结构等。近年来，非线性科学取得了长足的进展，尤其是混浊理论的建立，为各个学 科领域展现了新的视野。混浊是当今的前沿课题及学术热点，他揭示了自然界及人类社会 普遍存在的复杂性、有序性与无序性的统一、确定性与随机性的统一，拓展了人们的视野， 加深了对客观世界的认识。混浊涵盖了自然界及社会科学等各个领域，其研究的进展将有 力地促进几乎所有学科和
的发展川。二十世纪80年代，混浊科学在得到广泛发展 的同时，又与其他学科相互渗透，与时序分析的研究相交叉的结果是产生了"混浊时间序列 分析"。混浊时间序列分析在水文系统中的应用研究是一项极有意义的开创性工作。对水 文气象科学来讲，运用混浊理论中新的思想和方法，为水文科学的研究注入了新的活力，尤 其是实验分析研究，可W说已经从过去的半经验理论和统计理论研究步入了系统动力学理 论研究的阶段。 混浊是出于确定映射的似随机，物理学上，混浊通常被认为是确定的、耗散的非线 性动力系统中无序的、不可预知的行为。该里应该强调的是：混浊不能被简单地等同于无 序，更确切地说它是一种非线性的有序。混浊学研究的是无序中的有序，许多现象即使遵 循严格的确定性规则，但大体上仍是无法预测的，比如大气中的端流、中长期天气等，混浊 现象在不同的时间标度下表现出相似的模式，该与分形在空间标度下表现的相似性十分相 像。 混浊主要讨论非线性动力系统的不稳定、发散的过程，但系统状态在相空间中总 是收敛于一定的吸引子，该与分形的生成过程十分相像径流过程是一种复杂的水文现象， 其表现为强烈的非线性特征。要完整而准确地描述该样一种复杂的非线性过程，传统的欧 式几何显得力不从屯、，W混浊为核屯、的非线性科学，成了描述该类过程的有力工具。作为水 文系统一个重要的分支，河川径流情势的时空变化，是一个十分复杂的非线性与不确定性 干扰过程。从大尺度水文问题来看，河川径流变化受到全球气候变化影响，并且对气候有十 分明显的非线性反馈作用。从中尺度水文问题来看，河川径流变化不仅与气候因子（降水、 气温等）有关，而且受下垫面、地形、地貌W及人类活动影响。河道洪水波运动所表征的高 度非线性，使得径流描述需要考虑不同来源的非线性动力学影响。随着现代科学技术的不 断发展，高效率、大面积范围的收集水文信息已经成为可能。与此同时，水文非线性系统识 别与参数估计理论的研究也极大地丰富了系统水文学分析途径。正是由于径流系统内部因 素相互作用的复杂性W及影响其演变过程的径流过程是一种复杂的水文现象诸多外部因 素的难W处理，使得径流系统的演变特征难W理解和刻画。于是，构成径流过程不同因素的 相互作用产生的可观察结果一河川径流时间序列，就成为进行径流特性提取的有效研究对 象。 水文系统是一个开放的、复杂的巨系统，同时又是一个动态的非线性复合系统。一 方面，它是地球大气圈环境内相互作用和依赖的若干水文要素组成的具有水文循环与演化 功能的整体；另一方面，它又受地球及宇宙自然力的作用及来自人类的不同程度的生产活 动的影响，从而形成了水文系统复杂的演化规律。然而，由于哲学观和科学技术方法论的限 审IJ，长期W来，人们一直用传统的确定性方法或随机性方法，或将二者结合的方法来描述水 文过程，得W揭示的是水文系统的确定性规律。根据水文要素变化的非线性特点，引进新的 分析途径是十分必要的。 混浊理论和水文科学的结合则产生了一个新的研究领域，混浊理论的引入可W丰 富水文学的研究内容，推动水文科学的发展。混浊理论开启了探索水文现象变化的新途径， 通过应用混浊理论中的相空间重构技术，把水文时间序列嵌入到重构的相空间中，便可W 在相空间中揭示出水文动力系统复杂运动特征，该样就可能从复杂水文系统运动中发现其 内在的、有序的、确定性规律。由于该方面的研究尚处于初级阶段，认识还不够深入，需要研 究和解决的问题很多。因此，开展混浊理论及其应用的研究，从混浊动力学的角度去认识水 文系统的演变规律具有重要的现实意义和科学价值。 现有的混浊分析方法，一般要求被分析对象应包含较多数据，如至少上千，或更 多。但由于历史原因，相当数量的江河海域的径流数据记载只有上百或几十年。该对于按 年测量的径流量（如径流总量等）来说，就仅有上百甚至几十个数据。实验证明，该些数据 量不大的年径流量也同样具有混浊特性。此外，怎么将径流的混浊特性分析与预测技术结 合W提高径流预测模型的性能也鲜有研究。
技术实现思路
[000引本专利技术的目的在于提供混浊特性分析及非线性预测方法，解决了现有方法对于数 据量不大的江河海域的径流记录数据，无法进行混浊特性分析及非线性预测的问题。 本专利技术所采用的技术方案是按照W下步骤进行： 步骤1;采集径流数据，包括影响径流演化的所有相关数据，收集径流数据，径流 数据为每年采集一次，采集的径流数据大于39个； 步骤2 ;对每个径流数据进行小波变换，生成若干个子序列；[001引步骤3 ;对每个子序列采用小数量法或其他方法分别计算最大lyapunov指数对 最大lyapunov指数为正数的子序列判定为混浊时间序列，分别计算混浊时间序列的延迟 时间和嵌入维数； 步骤4;根据计算出来的延迟时间和嵌入维数，对相应的具有混浊特性的子序列 进行相空间重构，相空间重构表后的信息可W表示为维数为m的时延向量： X(t) = (x(t), x(t-S),…，x(t-(m-l)巧）T(2) 其中m为嵌入维数，S为延迟时间； 步骤5 ;将重构结果即式（2)作为预测的数据输入神经网络模型完成最终预测。 进一步，所述步骤3中计算Lyapunov指数过程如下；[001引对时间序列{x^X2,…，xj进行FFT变换，计算出平均周期P，进行相空间重构，寻 找相空间中每个点Xi的最近邻点X并限制短暂分离，即 其中Ij-n(j)I>p，p为混浊时间序列平均周期，对相空间中每个点X斯最近邻 点对的i个离散时间步后的距离dj.(i)，对每个i，求出所有j的Indj.(i)平均值，并用最小 二乘法作出加归直线，该直线的斜率就是最大Lyapunov指数入。 进一步，所述步骤3中混浊时间序列的延迟时间和嵌入维数计算方法为： 使用互信息法确定延迟时间： 绘制T~I(T)曲线，将曲线第一次降低到极小值时的延迟时间Tm。乃作为相 空间重构的延迟时间； 使用Cao算法计算嵌入维数： 该里用U范数，记所有a(i，m)关于i的均值为 [002引E(m)只依赖于嵌入维数m和延迟时间T，为研究嵌入维数从m变为m+1时相空间 的变化情况，定义[003U如果当m大于某个m。时，Ei(m)停止变化，则就是重构相空间的最小嵌入维 数，Cao定义量； 本专利技术的有益效果是江河海域的径流记录数据大于39个时，即可进行混浊特性 分析及非线性预测。【附图说明】 图1是本专利技术基于径流的混浊特性分析及非线性预测系统仿真实验中柳江径流 总量不意图； 图2是本专利技术仿真实验四种模型拟合结果绝对误差比较；[0本文档来自技高网...

【技术保护点】
径流的混沌特性分析及非线性预测方法，其特征在于：按照以下步骤进行：步骤1：采集径流数据，包括影响径流演化的所有相关数据，收集径流数据，径流数据为每年采集一次，采集的径流数据大于39个；步骤2：对每个径流数据进行小波变换，生成若干个子序列；步骤3：对每个子序列采用小数量法或其他方法分别计算最大lyapunov指数；对最大lyapunov指数为正数的子序列判定为混沌时间序列，分别计算混沌时间序列的延迟时间和嵌入维数；步骤4：根据计算出来的延迟时间和嵌入维数，对相应的具有混沌特性的子序列进行相空间重构，相空间重构表后的信息表示为维数为m的时延向量：X(t)＝(x(t),x(t‑S),…,x(t‑(m‑1)S))T   (2)其中m为嵌入维数，S为延迟时间；步骤5：将重构结果即式(2)作为预测的数据输入神经网络模型完成最终预测。

【技术特征摘要】

【专利技术属性】
技术研发人员：丁红，
申请(专利权)人：柳州师范高等专科学校，
类型：发明
国别省市：广西;45