一种基于集合卡尔曼滤波的对无资料地区径流量推求方法技术

本发明专利技术公开了一种基于集合卡尔曼滤波的对无资料地区径流量推求方法，涉及水文模拟技术领域。所述方法：计算全流域的优化参数值；计算子流域A的出口点C点的估计流量及其分布规律；计算变量和参数的集合扰动值；计算下游子流域B的出口点即水文站D的径流量预报集合；对水文站D的实测径流值进行高斯扰动，并将得到高斯扰动值作为观测数据集合；融入观测数据集合进行同化更新；最终得到C点径流量QC的最优估计值。本发明专利技术简化了滤波同化过程，较好的提高了同化精度；同时也避免了常规方法对无资料地区求径流量的缺陷。

【技术实现步骤摘要】
一种基于集合卡尔曼滤波的对无资料地区径流量推求方法
本专利技术涉及水文模拟
，尤其涉及一种基于集合卡尔曼滤波的对无资料地区径流量推求方法。
技术介绍
自然界中有不少区域因雨量或水文站早期建设少、历史径流资料序列短，无法满足预报需求；同时，因人类活动等影响，导致流域产汇流特性发生较大变化，历史资料无法反映当前来水特性，历史资料不可用。因此对无资料地区的研究也是亟待解决的。目前，对无资料地区径流预报常用的方法为区域化方法，即通过某种途径，利用有资料流域的模型参数推求无资料流域的模型参数，从而对无资料流域进行预报。常用的区域化方法有属性相似法、空间相近法和回归法。其中，属性相似法指找出与研究流域属性如地形、土壤和气候等相似的流域，并把其参数作为研究流域的参数，属性相似法研究根据为同一区域的物理和气候属性相对一致，因此相邻流域的水文行为相似。空间相近法指找出与研究流域(如无资料流域)距离上相近的一个或者多个有资料流域，并把有资料流域参数作为研究流域的参数。回归法指根据有资料流域的模型参数和流域属性，建立二者之间的多元回归方程，从而利用无资料流域的流域属性推求其模型参数。但是，现有技术存在如下缺点：一、未考虑不同参数的区别：目前，无资料地区参数确定通常采用移用或回归分析的方法。现有，确定移用时通常先找到参照流域，而后所有参数统一移用或回归分析。可是实际上，水文模型中每个参数代表的物理意义不同，因每个流域都有其独特性，很难找到与模型中所有参数的相关特征均相似的流域，虽现今已有研究仅依据一种或几种相关特征寻找参照流域，而所有参数统一移用或采用同一些流域，但这种方式很明显无法反映参数的独特性。二、相似流域的选取缺乏全面性和客观性：现有研究通常选用一种或几种地形、地貌特征参数的综合指数如地形指数作为属性判别寻找相似流域，而后进行参数的移用或分析。因为不同水文模型其参数的物理意义不同，仅仅依据一种或几种流域特征参数进行判别是不够全面和客观的，很难反映真正的流域相似，并且采用所有流域特征进行分析存在信息量大、分类复杂以及无法判断的问题。三、主观判别和效率低：现有技术均采用人工判别处理的方式，解决无资料地区的洪水预报问题，此技术难以大量推广，并且其中包含人为的主观性和效率低下的问题，寻找其他有效简便的方法是目前急需解决的。
技术实现思路
本专利技术的目的在于提供一种基于集合卡尔曼滤波的对无资料地区径流量推求方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术所述基于集合卡尔曼滤波的对无资料地区径流量推求方法，所述方法包括：S1，计算全流域的优化参数值全流域中存在至少一个无水文资料子流域和至少一个有水文资料子流域，所述无水文资料的子流域的上游流域和下游流域均为全流域中子流域；对全流域进行参数优化，得到参数最优值，选取任意两个参数的最优值进行数据同化处理；两个参数分别为参数β和参数γ；S2，计算子流域A的出口的估计流量及其分布规律设定无水文资料子流域为子流域A，子流域A的出口点为C点，子流域A的下游流域为有水文资料子流域B，子流域B中的出口点为水文站D；子流域B出口对应的径流量是已知的，记为QD；然后，通过公式(1)计算子流域A的出口点C点的初始估计流量QC′；QC′＝QD×(SA/SA+B)(1)；其中，SA为子流域A的面积，SA+B为子流域A和子流域B的面积之和；S3，计算变量和参数的集合扰动值将初始估计流量QC′的均值和方差作为C点径流量QC初始设定的符合高斯分布的均值与方差；以参数β和参数γ的最优值为参考，对C点径流量QC、参数β和参数γ分别进行高斯扰动，分别得到变量、参数β和参数γ的高斯扰动值集合；S4，计算下游有水文资料水文站D的径流量预报集合在实测降水和蒸发数据的基础上，逐天运行水文模型，将S3中得到的C点径流量QC、参数β和参数γ的高斯扰动集合值逐一代入水文模型，并得到每个时间步长上的水文站D的径流量预报集合，每个时间步长上的水文站D的径流量预报集合数量大于等于1；S5，对水文站D实测径流值进行高斯扰动，并将得到的高斯扰动集合作为观测数据集合；S6，融入观测数据集合在已知观测算子的条件下，计算卡尔曼增益，将每个所述水文站D的径流量预报集合融入观测数据集合进行更新；S7，重复进行S4至S6，直到水文模型计算结束，最终得到C点径流量QC的最优估计值，所述最优估计值即为上游无资料子流域A出口点C点的径流量值。优选地，S1中，对全流域进行参数优化采用的优化算法为所述粒子群算法。优选地，所述水文模型为新安江模型。优选地，步骤S3中，所述高斯扰动，按照公式(2)实现：Ea＝ea+ε(2)；其中，Ea表示增加高斯扰动后的模型状态变量或者观测值；ea表示初始模型状态变量或者观测值；ε为符合高斯分布的扰动值，即ε～N(0，1)。优选地，步骤S6中，计算卡尔曼增益，具体按照下述计算：Kt＝PtHT(HPtHT+Rt)-1(9)；其中，m表示集合样本的总数；表示t时刻m个径流量预报集合的预报均值；表示t时刻水文模型第i个径流量预报集合中集合成员的一个预报状态向量；表示t时刻每一个径流量预报集合的集合成员与的扰动；yi,t表示t时刻第i个集合成员的观测数据向量；表示t时刻m个观测数据集合的观测均值；Δyi,t表示t时刻每一个观测数据的集合成员与的扰动；Pt表示t时刻的径流量预报集合的误差协方差；Rt表示t时刻观测数据集合的误差协方差矩阵；H为观测算子；Kt为卡尔曼增益，表示观测数据的权重；T表示观测算子的转置。优选地，步骤S6中，对每一个集合进行更新，采用公式(10)进行更新：其中，表示t时刻状态下的水文站D的径流量的最优估计集合，的均值即为t时刻状态下的最优估计值，表示t时刻的水文站D的径流量预报集合，yt表示t时刻的观测数据集合，Kt为卡尔曼增益，表示观测数据的权重。本专利技术的有益效果是：本专利技术利用集合卡尔曼滤波推求无资料地区径流量的研究，为无资料地区的研究展开了新思路和新方法。集合卡尔曼滤波的基本思想即是利用MonteCarlo方法设计预测状态的一个集合，该集合的平均可作为最佳估计，该集合的样本协方差即作为预测误差协方差的近似，该集合通过不断向前滤波，每个样本分别更新分析变量，而对变量的最佳估计即为各更新分析变量的样本平均。本专利技术提出的基于集合卡尔曼滤波的对无资料地区径流量推求方法，基于集合统计的思想，无须对非线性系统进行线性化，避免了Jacobian矩阵的繁冗计算，更简化了滤波同化过程。本专利技术利用一种扩展模型状态变量方法，即将模型参数和变量置于联合的向量中，将上游无资料地区径流量看作变量，利用下游已知流域的参数并将其径流量作为观测值，同时实现变量与参数的同化，其中变量同化的均值即为无资料地区的径流量。基于集合统计的思想，无须对非线性系统进行线性化，避免了Jacobian矩阵的繁冗计算，并且同化过程中同时考虑输入降雨误差、模型误差以及观测误差，较好的提高了同化精度。同时也避免了常规方法对无资料地区求径流量的缺陷，如相似流域选取缺乏全面、客观性，现有技术未考虑不同参数的区别，主观判别、效率低等，该方法将会对无资料地区预报完善洪水预报理论、提高预报精度、水库大坝建设和调度设计等具有理论价值和指导意义。附图说明图1是全本文档来自技高网...

【技术保护点】
一种基于集合卡尔曼滤波的对无资料地区径流量推求方法，其特征在于，所述方法包括：S1，计算全流域的优化参数值全流域中存在至少一个无水文资料子流域和至少一个有水文资料子流域，所述无水文资料的子流域的上游流域和下游流域均为全流域中子流域；对全流域进行参数优化，得到参数最优值，选取任意两个参数的最优值进行数据同化处理；两个参数分别为参数β和参数γ；S2，计算子流域A的出口的估计流量及其分布规律设定无水文资料子流域为子流域A，子流域A的出口点为C点，子流域A的下游流域为有水文资料子流域B，子流域B中的出口点为水文站D；子流域B出口对应的径流量是已知的，记为QD；然后，通过公式(1)计算子流域A的出口点C点的初始估计流量QC′；QC′＝QD×(SA/SA+B)            (1)；其中，SA为子流域A的面积，SA+B为子流域A和子流域B的面积之和；S3，计算变量和参数的集合扰动值将初始估计流量QC′的均值和方差作为C点径流量QC初始设定的符合高斯分布的均值与方差；以参数β和参数γ的最优值为参考，对C点径流量QC、参数β和参数γ分别进行高斯扰动，分别得到变量、参数β和参数γ的高斯扰动值集合；S4，计算下游有水文资料水文站D的径流量预报集合在实测降水和蒸发数据的基础上，逐天运行水文模型，将S3中得到的C点径流量QC、参数β和参数γ的高斯扰动集合值逐一代入水文模型，并得到每个时间步长上的水文站D的径流量预报集合，每个时间步长上的水文站D的径流量预报集合数量大于等于1；S5，对水文站D的实测径流值进行高斯扰动，并将得到的高斯扰动集合作为观测数据集合；S6，融入观测数据集合在已知观测算子的条件下，计算卡尔曼增益，将每个所述水文站D的径流量预报集合融入观测数据集合进行更新；S7，重复进行S4至S6，直到水文模型计算结束，最终得到C点径流量QC的最优估计值，所述最优估计值即为上游无资料子流域A出口点的径流量值。...

【技术特征摘要】
1.一种基于集合卡尔曼滤波的对无资料地区径流量推求方法，其特征在于，所述方法包括：S1，计算全流域的优化参数值全流域中存在至少一个无水文资料子流域和至少一个有水文资料子流域，所述无水文资料的子流域的上游流域和下游流域均为全流域中子流域；对全流域进行参数优化，得到参数最优值，选取任意两个参数的最优值进行数据同化处理；两个参数分别为参数β和参数γ；S2，计算子流域A的出口的估计流量及其分布规律设定无水文资料子流域为子流域A，子流域A的出口点为C点，子流域A的下游流域为有水文资料子流域B，子流域B中的出口点为水文站D；子流域B出口对应的径流量是已知的，记为QD；然后，通过公式(1)计算子流域A的出口点C点的初始估计流量QC′；QC′＝QD×(SA/SA+B)(1)；其中，SA为子流域A的面积，SA+B为子流域A和子流域B的面积之和；S3，计算变量和参数的集合扰动值将初始估计流量QC′的均值和方差作为C点径流量QC初始设定的符合高斯分布的均值与方差；以参数β和参数γ的最优值为参考，对C点径流量QC、参数β和参数γ分别进行高斯扰动，分别得到变量、参数β和参数γ的高斯扰动值集合；S4，计算下游有水文资料水文站D的径流量预报集合在实测降水和蒸发数据的基础上，逐天运行水文模型，将S3中得到的C点径流量QC、参数β和参数γ的高斯扰动集合值逐一代入水文模型，并得到每个时间步长上的水文站D的径流量预报集合，每个时间步长上的水文站D的径流量预报集合数量大于等于1；S5，对水文站D的实测径流值进行高斯扰动，并将得到的高斯扰动集合作为观测数据集合；S6，融入观测数据集合在已知观测算子的条件下，计算卡尔曼增益，将每个所述水文站D的径流量预报集合融入观测数据集合进行更新；S7，重复进行S4至S6，直到水文模型计算结束，最终得到C点径流量QC的最优估计值，所述最优估计值即为上游无资料子流域A出口点的径流量值。2.根据权利要求1所述方法，其特征在于，S1中，对全流域进行参数优化采用的优化算法为所述粒子群算法。3.根据权利要求1所述方法，其特征在于，所述水文模型为新安江模型。4.根据权利要求1所述方法，其特征在于，步骤S3中，所述高斯扰动，按照公式(2)...

【专利技术属性】
技术研发人员：雷晓辉，廖卫红，张苹苹，谢先红，王明元，殷兆凯，秦韬，张云辉，鲍淑君，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京,11