本发明专利技术公开了一种流域洪水模拟计算的实时校正方法、系统、装置及介质，其方法包括：确定模拟计算的开始时刻和目标时段；初始化模拟计算的开始时刻的土壤含水量，并输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值；计算目标时段内洪量的模拟计算值；获取目标时段内洪量的实测值，并计算目标时段内洪量的模拟计算值与实测值之间的误差值；判断误差值是否小于等于误差阈值，若是，则输出土壤含水量最优值；若否，则对土壤含水量进行校正并进行迭代；将土壤含水量最优值输入流域水文模型获取模拟计算时段内流量的模拟计算值最优值；从中获取当前时刻之后的预报时段内流量的模拟计算最优值；本发明专利技术能够提高实时洪水预报的成果精度。报的成果精度。报的成果精度。
【技术实现步骤摘要】
流域洪水模拟计算的实时校正方法、系统、装置及介质
[0001]本专利技术涉及一种流域洪水模拟计算的实时校正方法、系统、装置及介质，属于洪水预报
技术介绍
[0002]降雨产生洪水的量级与降雨量大小、前期土壤含水量和流域特征有关。水文工作者根据流域内雨量站的实际降雨量，使用流域水文模型进行流域产汇流计算，模拟计算预报流域出口断面位置未来一段时间内的流量过程。我国南方流域多使用新安江蓄满产流模型，只有土壤蓄满部分的流域产流，而预报时刻的土壤含水量是未知的，通常对模型进行长时间演算推求预报时刻的土壤含水量，也因流域蒸散发不确定性，容易造成径流计算误差。实时预报系统往往使用实时校正技术来减小误差，例如自回归模型、卡尔曼滤波、神经网络算法等，使用当前的误差去预测未来的误差，校正模型预报计算值。将误差作为时间序列，采用数值统计和状态估计方法，预测未来误差，但上述方法没有抓住误差成因或者是算法太复杂，往往效果不好。
技术实现思路
[0003]本专利技术的目的在于克服现有技术中的不足，提供一种流域洪水模拟计算的实时校正方法、系统、装置及介质，解决由于土壤含水量误差导致洪水模拟计算不准确的技术问题。
[0004]为达到上述目的，本专利技术是采用下述技术方案实现的：
[0005]第一方面，本专利技术提供了一种流域洪水模拟计算的实时校正方法，包括：
[0006]获取模拟计算的当前时刻，并基于当前时刻确定模拟计算的开始时刻；
[0007]获取当前时刻的洪峰状态，并结合开始时刻确定模拟计算的目标时段；
[0008]初始化模拟计算的开始时刻的土壤含水量，并输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值；
[0009]根据模拟计算时段内各时刻流量的模拟计算值计算目标时段内洪量的模拟计算值；
[0010]获取目标时段内洪量的实测值，并计算目标时段内洪量的模拟计算值与实测值之间的误差值；
[0011]判断误差值是否小于等于误差阈值，若是，则输出土壤含水量最优值；若否，则对土壤含水量进行校正，并将校正后的土壤含水量替换初始化的土壤含水量且返回步骤：初始化模拟计算的开始时刻的土壤含水量；
[0012]将土壤含水量最优值输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值最优值；
[0013]从模拟计算时段内各时刻流量的模拟计算值最优值中获取当前时刻之后的预报时段内各时刻流量的模拟计算最优值。
[0014]优选的，所述模拟计算的开始时刻为：模拟计算的当前时刻之前时长T内的最小流量出现的时刻，所述时长T＝流域预见期+48h。
[0015]优选的，所述模拟计算的目标时段为：
[0016]若当前时刻的洪峰状态为洪峰未出现，则模拟计算的目标时段为当前时刻之前一个流域预见期；
[0017]若当前时刻的洪峰状态为洪峰已出现，但洪峰出现时刻至当前时刻未超过半个流域预见期，则模拟计算的目标时段为当前时刻之前一个流域预见期；
[0018]若当前时刻的洪峰状态为洪峰已出现，但洪峰出现时刻至当前时刻超过半个流域预见期，则模拟计算的目标时段为洪峰出现时刻前后各半个流域预见期。
[0019]优选的，所述初始化模拟计算的开始时刻的土壤含水量θ为：
[0020]θ＝a*Q
b
+c
[0021]式中，a、b、c为常数参数，a的取值为1
‑
5，b的取值为0.15
‑
0.65，c的取值为30
‑
60，Q为模拟计算的开始时刻的断面流量。
[0022]优选的，所述对土壤含水量进行校正的表达式为：
[0023][0024]式中，θ
t+1
、θ
t
分别为第t+1、t次迭代时模拟计算的开始时刻的土壤含水量，W
m
、W
s
分别为目标时段内洪量的模拟计算值、实测值。
[0025]第二方面，本专利技术提供了一种流域洪水模拟计算的实时校正系统，所述系统包括：
[0026]开始时刻确定模块，用于获取模拟计算的当前时刻，并基于当前时刻确定模拟计算的开始时刻；
[0027]目标时段确定模块，用于获取当前时刻的洪峰状态，并结合开始时刻确定模拟计算的目标时段；
[0028]第一模拟计算模块，用于初始化模拟计算的开始时刻的土壤含水量，并输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值；
[0029]第一获取模块，用于根据模拟计算时段内各时刻流量的模拟计算值计算目标时段内洪量的模拟计算值；
[0030]误差计算模块，用于获取目标时段内洪量的实测值，并计算目标时段内洪量的模拟计算值与实测值之间的误差值；
[0031]土壤含水量校正模块，用于判断误差值是否小于等于误差阈值，若是，则输出土壤含水量最优值；若否，则对土壤含水量进行校正，并将校正后的土壤含水量替换初始化的土壤含水量且返回第一模拟计算模块中循环迭代；
[0032]第二模拟计算模块，用于将土壤含水量最优值输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值最优值；
[0033]第二获取模块，用于从模拟计算时段内各时刻流量的模拟计算值最优值中获取当前时刻之后的预报时段内各时刻流量的模拟计算最优值。
[0034]第三方面，本专利技术提供了一种流域洪水模拟计算的实时校正装置，包括处理器及存储介质；
[0035]所述存储介质用于存储指令；
[0036]所述处理器用于根据所述指令进行操作以执行根据上述方法的步骤。
[0037]第四方面，本专利技术提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述方法的步骤。
[0038]与现有技术相比，本专利技术所达到的有益效果：
[0039]本专利技术提供的一种流域洪水模拟计算的实时校正方法、系统、装置及介质，通过确定合适的模型计算开始时刻，保证开始时刻的土壤含水量对后期预报时段内的模拟计算产生影响；提出了初始土壤含水量的逐步迭代过程，根据当前时刻与洪峰出现时刻分情况选定目标时段，使目标时段内的模拟成果接近实测成果，从而优化土壤含水量；通过优化后的土壤含水量提高了实时洪水预报的成果精度。
附图说明
[0040]图1是本专利技术实施例一提供的一种流域洪水模拟计算的实时校正方法的流程图；
[0041]图2是本专利技术实施例一提供的洪峰未出现的流量示意图；
[0042]图3是本专利技术实施例一提供的洪峰出现且洪峰出现时刻至当前时刻未超过半个流域预见期的流量示意图；
[0043]图4是本专利技术实施例一提供的洪峰出现且洪峰出现时刻至当前时刻超过半个流域预见期的流量示意图。
具体实施方式
[0044]下面结合附图对本专利技术作进一步描述。以下实施例仅用于更加清楚地说明本专利技术的技术方案，而不能以此来限制本专利技术的保护范围。
[0045]实施例一：
[0046]如图1所示，本专利技术本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种流域洪水模拟计算的实时校正方法，其特征在于，包括：获取模拟计算的当前时刻，并基于当前时刻确定模拟计算的开始时刻；获取当前时刻的洪峰状态，并结合开始时刻确定模拟计算的目标时段；初始化模拟计算的开始时刻的土壤含水量，并输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值；根据模拟计算时段内各时刻流量的模拟计算值计算目标时段内洪量的模拟计算值；获取目标时段内洪量的实测值，并计算目标时段内洪量的模拟计算值与实测值之间的误差值；判断误差值是否小于等于误差阈值，若是，则输出土壤含水量最优值；若否，则对土壤含水量进行校正，并将校正后的土壤含水量替换初始化的土壤含水量且返回步骤：初始化模拟计算的开始时刻的土壤含水量；将土壤含水量最优值输入流域水文模型获取模拟计算时段内各时刻流量的模拟计算值最优值；从模拟计算时段内各时刻流量的模拟计算值最优值中获取当前时刻之后的预报时段内各时刻流量的模拟计算最优值。2.根据权利要求1所述的一种流域洪水模拟计算的实时校正方法，其特征在于，所述模拟计算的开始时刻为：模拟计算的当前时刻之前时长T内的最小流量出现的时刻，所述时长T＝流域预见期+48h。3.根据权利要求1所述的一种流域洪水模拟计算的实时校正方法，其特征在于，所述模拟计算的目标时段为：若当前时刻的洪峰状态为洪峰未出现，则模拟计算的目标时段为当前时刻之前一个流域预见期；若当前时刻的洪峰状态为洪峰已出现，但洪峰出现时刻至当前时刻未超过半个流域预见期，则模拟计算的目标时段为当前时刻之前一个流域预见期；若当前时刻的洪峰状态为洪峰已出现，但洪峰出现时刻至当前时刻超过半个流域预见期，则模拟计算的目标时段为洪峰出现时刻前后各半个流域预见期。4.根据权利要求1所述的一种流域洪水模拟计算的实时校正方法，其特征在于，所述初始化模拟计算的开始时刻的土壤含水量θ为：θ＝a*Q
b
+c式中，a、b、c为常数参数，a的取值为1
‑
5，b的取值为0.15
‑
0.65，c的取值为30
...
【专利技术属性】
技术研发人员：舒大兴，
申请(专利权)人：河海大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。