CN117236216B - 水库水动力模型糙率率定方法、装置、计算机设备及介质 - Google Patents

Abstract

本发明涉及模型参数率定技术领域，公开了水库水动力模型糙率率定方法、装置、计算机设备及介质，方法基于各断面分区的初始糙率值以及预设糙率率定步长确定各断面分区分别对应的第一糙率值序列，基于水位计算边界信息以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。通过各断面分区的第一糙率值序列计算各糙率值对应的第二水位信息，基于各糙率值对应的第二水位信息以及第一水位信息确定各断面分区的糙率率定值，糙率率定结果更合理。

Description

水库水动力模型糙率率定方法、装置、计算机设备及介质

技术领域

本发明涉及模型参数率定技术领域，具体涉及水库水动力模型糙率率定方法、装置、计算机设备及介质。

背景技术

糙率是水动力模型计算中一个相当重要参数，对河道水流及其冲淤变化的计算结果影响很大。通过对水动力模型的糙率率定可以有效提高模型的计算精度。

现有水动力模型的糙率率定方案一般通过迭代试算各监测点水位的模型计算值与实测值之间的误差，判断水位误差是否小于预设阈值，若水位误差大于预设阈值，对糙率值进行更新，根据更新后的糙率值计算监测点水位的模型计算值，直至监测点水位的模型计算值与实测值之间的误差小于预设阈值，将当前糙率值确定为糙率率定结果。由于预设阈值需要人为确定，使得模型糙率率定结果不具备科学性。

发明内容

有鉴于此，本发明提供了一种水库水动力模型糙率率定方法、装置、计算机设备及介质，以解决相关技术中通过人为确定的预设阈值率定模型糙率导致糙率率定结果不具备科学性的问题。

第一方面，本发明提供了一种水库水动力模型糙率率定方法，该方法包括：获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，第一水位信息为实测水位信息；基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值；基于水位计算边界信息、水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。

本发明提供的水库水动力模型糙率率定方法，基于各断面分区的初始糙率值以及预设糙率率定步长确定各断面分区分别对应的第一糙率值序列，基于水位计算边界信息以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。本发明提供的方法，通过各断面分区的第一糙率值序列计算各糙率值对应的第二水位信息，基于各糙率值对应的第二水位信息以及第一水位信息确定各断面分区的糙率率定值，使得糙率率定结果更合理，解决了相关技术中通过人为确定的预设阈值率定模型糙率导致糙率率定结果不具备科学性的问题。

在一种可选的实施方式中，根据各糙率值对应的第二水位与第一水位的误差，确定对应断面分区的糙率率定值的步骤，包括：获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序；令N＝1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值；判断N是否小于M，若N小于M，令N＝N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值；返回判断N是否小于M的步骤，直至N＝M，得到各断面的第一糙率率定值。

本可选实施方式提供的方法，通过各断面分区的第一糙率率定值序列以及第1至N个断面的第一糙率率定值计算待率定的第N+1个分区的第一糙率率定值，有效减少了由于已率定分区的糙率率定值误差的累积，对后续待率定分区第一糙率率定值的影响，保证了各断面分区第一糙率率定值的准确性。

在一种可选的实施方式中，根据各糙率值对应的第二水位与第一水位的误差，确定对应断面分区的糙率率定值的步骤，还包括：将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值；若第二目标数量大于第一目标数量，将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量小于第一目标数量，得到各断面分区的第二糙率率定值，第二目标数量用于表征所述第一差异值大于第一预设值的断面分区的数量。

本可选实施方式提供的方法，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值，当第一差异值大于第一预设值的断面分区的数量大于第一目标数量，将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值再次对各断面分区进行率定，保证了各断面分区糙率率定值的准确性。

在一种可选的实施方式中，断面分区是根据水库的监测站点分布信息对水库沿程断面进行划分得到的。

在一种可选的实施方式中，每个分区的糙率率定顺序通过如下步骤获取：获取水库的沿程断面的水流特征信息以及每个断面分区的位置信息；基于水流特征信息确定每个断面分区的分区类型；基于每个分区的位置信息以及分区类型确定对应分区的糙率率定顺序。

本可选实施方式提供的方法，基于水库的沿程断面的水流特征信息以及每个断面分区的位置信息确定各断面分区的糙率率定顺序，使得分区糙率率定更符合物理规律。

在一种可选的实施方式中，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值的步骤，包括：将第1个断面分区的各糙率值对应的第二水位信息和对应第一水位信息进行误差计算，得到第一糙率值序列中各糙率值分别对应的第一误差计算结果；基于第一误差计算结果，将第一糙率值序列中误差最小的糙率值作为第1个断面的第一糙率率定值。

在一种可选的实施方式中，将第1个断面分区在每个糙率值条件下计算得到的第二水位变化信息和对应第一水位信息进行误差计算，得到第一误差计算结果的步骤，包括：基于第1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；基于第1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第1个断面分区在每个糙率值条件下的第二水位变化曲线；基于第1个断面分区的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第1个断面分区在每个糙率值条件下水位的相位误差，得到第一误差计算结果。

本可选实施方式提供的方法，基于第二水位变化曲线以及第一水位变化曲线对断面分区的糙率进行率定，提高了糙率率定的计算速度和计算结果的准确性。

在一种可选的实施方式中，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值的步骤，包括：分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果；基于第二误差计算结果，将第一糙率值序列中误差最小的糙率值作为第N+1个断面分区的糙率率定值。

在一种可选的实施方式中，分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果的步骤，包括：基于第N+1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；基于第N+1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第N+1个断面分区在每个糙率值条件下的第二水位变化曲线；基于第N+1个断面的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第N+1个断面在每个糙率值条件下水位的相位误差，得到第二误差计算结果。

第二方面，本发明提供了一种水库水动力模型糙率率定装置，该装置包括：获取模块，用于获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，第一水位信息为实测水位信息；第一确定模块，用于基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值；计算模块，用于基于水位计算边界信息、水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；第二确定模块，用于根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。

在一种可选的实施方式中，第二确定模块，包括：排序子模块，用于获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序；第一确定子模块，用于令N＝1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值；第二确定子模块，用于判断N是否小于M，若N小于M，令N＝N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值；第三确定子模块，用于返回判断N是否小于M的步骤，直至N＝M，得到各断面的第一糙率率定值。

在一种可选的实施方式中，第二确定模块，还包括：比对子模块，用于将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值；第四确定子模块，用于若第二目标数量大于第一目标数量，将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量小于第一目标数量，得到各断面分区的第二糙率率定值，第二目标数量用于表征第一差异值大于第一预设值的断面分区的数量。

第三方面，本发明提供了一种计算机设备，包括：存储器和处理器，存储器和处理器之间互相通信连接，存储器中存储有计算机指令，处理器通过执行计算机指令，从而执行上述第一方面或其对应的任一实施方式的水库水动力模型糙率率定方法。

第四方面，本发明提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机指令，计算机指令用于使计算机执行上述第一方面或其对应的任一实施方式的水库水动力模型糙率率定方法。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是根据本发明实施例的水库水动力模型糙率率定方法的流程示意图；

图2是根据本发明实施例的另一水库水动力模型糙率率定方法的流程示意图；

图3为A河段站点实测水位与糙率率定后的计算水位对比图；

图4为I河段站点实测水位与糙率率定后的计算水位对比图；

图5为P河段站点实测水位与糙率率定后的计算水位对比图；

图6是根据本发明实施例的水库水动力模型糙率率定装置的结构框图；

图7是本发明实施例的计算机设备的硬件结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

相关技术中，通过迭代试算各监测点水位的模型计算值与实测值之间的误差，判断水位误差是否小于预设阈值，若水位误差大于预设阈值，对糙率值进行更新，根据更新后的糙率值计算监测点水位的模型计算值，直至监测点水位的模型计算值与实测值之间的误差小于预设阈值，将当前糙率值确定为糙率率定结果。由于预设阈值需要人为确定，使得模型糙率率定结果不具备科学性。

有鉴于此，本发明实施例提供了一种水动力模型糙率率定方法，过各断面分区的第一糙率值序列计算各糙率值对应的第二水位信息，基于各糙率值对应的第二水位信息以及第一水位信息确定各断面分区的糙率率定值，糙率率定结果更合理，解决了相关技术中通过人为确定的预设阈值率定模型糙率导致糙率率定结果不具备科学性的问题。

根据本发明实施例，提供了一种水库水动力模型糙率率定方法实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

在本实施例中提供了一种水库水动力模型糙率率定方法，可用于上述的处理器，图1是根据本发明实施例的水库水动力模型糙率率定方法的流程图，如图1所示，该流程包括如下步骤：

步骤S101，获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，第一水位信息为实测水位信息。

示例性地，水库可以包括但不限于河道型水库，非恒定流水动力模型是基于河道型水库的实测来流流量、出流流量、各监测站水位、地形断面资料等可信数据构建得到的；不同断面分区的初始糙率值可以是基于预设的规则确定的，本申请实施例中对不同断面分区的初始糙率值的确定方法不做限定，只要符合基本物理意义即可，糙率具有物理意义，反映的是河段的阻力特性，经验总结天然河道的阻力特性在0.025～0.07之间，一维水动力中用糙率可以表征河道由于各类原因(断面开形状突变、河流流向突变等)导致的水头损失，因此在取值时选定范围大于天然河道，所有河段可参考同类计算结果选择一致；水位计算边界信息可以包括但不限于预设时段内水库的实测来流流量、出流流量、坝前水位以及各监测站水位等数据；不同断面分区的第一水位信息是实测水位信息，不同断面分区的第一水位信息可以由设置在不同断面分区内的监测站点采集得到；糙率率定范围可以是糙率率定值的上限和下限，本申请实施例中，糙率率定值的上限可以包括但不限于0.1，下限可以包括但不限于0.02；预设糙率率定步长可以基于实际需求确定，本申请实施例中，预设糙率率定步长应不小于0.0025。

步骤S102，基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值。

示例性地，本申请实施例中，对于任一断面分区，将该断面分区的初始糙率值作为中值，按照预设糙率率定步长上下调整，调整后得到的糙率值序列不超过糙率率定范围。确定断面分区的第一糙率值序列时，首先确定当前糙率率定的断面分区，该分区需要进行计算比较的糙率取值有哪些，根据计算机的不同配置有不同的计算情况，若为多核并行计算，建议根据核心数设置计算工况数目，如32核，则设置不多于32个计算工况；若为单核，则建议采取二分法处理；现电脑CPU均为多核，建议采用多核并行计算；另外此处计算方法近似穷举，并不是最优的计算方法，但对于一维水动力计算而言，其糙率调整相互影响，且计算速度较快，最主要的是能快速生成计算工况，在这种条件下采用，每一次循环步采用固定步长，较为合适。

步骤S103，基于水位计算边界信息、水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息。

示例性地，本申请实施例中，水位计算边界信息作为边界条件，依次将任一断面分区对应第一糙率值序列中的多个糙率值依次输入到非恒定流水动力模型中，对模型求解，得到各糙率值对应的第二水位信息

步骤S104，根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。

示例性地，本申请实施例中，基于各糙率值对应的第二水位信息以及第一水位信息可以确定各糙率值的计算误差，将误差最小的糙率值作为糙率率定值。

本实施例提供的水库水动力模型糙率率定方法，基于各断面分区的初始糙率值以及预设糙率率定步长确定各断面分区分别对应的第一糙率值序列，基于水位计算边界信息以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。本发明提供的方法，通过各断面分区的第一糙率值序列计算各糙率值对应的第二水位信息，基于各糙率值对应的第二水位信息以及第一水位信息确定各断面分区的糙率率定值，糙率率定结果更合理，解决了相关技术中通过人为确定的预设阈值率定模型糙率导致糙率率定结果不具备科学性的问题。

在本实施例中提供了一种水库水动力模型糙率率定方法，可用于上述的处理器，图2是根据本发明实施例的水库水动力模型糙率率定方法的流程图，如图2所示，该流程包括如下步骤：

步骤S201，获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，第一水位信息为实测水位信息。详细请参见图1所示实施例的步骤S101，在此不再赘述。

步骤S202，基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值。详细请参见图1所示实施例的步骤S102，在此不再赘述。

步骤S203，基于水位计算边界信息、水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息。详细请参见图1所示实施例的步骤S103，在此不再赘述。

步骤S204，根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。详细请参见图1所示实施例的步骤S104，在此不再赘述。

具体地，上述步骤S204包括：

步骤S2041，获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序。

示例性地，每个断面分区的糙率率定顺序可以是基于预设需求确定的，本申请实施例对每个断面分区的糙率率定顺序不做限定，本领域技术人员可以根据需求确定。

步骤S2042，令N＝1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值。

示例性地，N可以是当前待率定断面分区的编号，当N＝1时，可以确定当前要率定的断面分区为第1个断面分区，基于该断面分区的第一糙率序列中各糙率值对应的第二水位信息和第一水位信息确定计算误差做小的糙率值，并将该误差最小的糙率值糙率第1个断面分区的第一糙率率定值。

步骤S2043，判断N是否小于M，若N小于M，令N＝N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值。

示例性地，M为断面分区的数量，当N小于M，令N＝N+1，并将第1至N个断面的第一糙率率定值输入到非恒定流水动力模型中，通过求解非恒定流水动力模型，得到第N+1个断面的第一糙率值序列中各糙率值对应的第二水位信息，基于各糙率值对应的第二水位信息和第一水位信息可以确定计算误差最小的糙率值，并将该计算误差最小的糙率值作为第N+1个断面的第一糙率率定值。

步骤S2044，返回判断N是否小于M的步骤，直至N＝M，得到各断面的第一糙率率定值。

示例性地，返回判断N是否小于M的步骤，直至N＝M，可以确定所有断面分区的第一糙率率定值。

步骤S2045，将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值。

示例性地，将每个断面分区的糙率率定值与对应断面分区的初始糙率值进行比对，可以确定各断面分区的第一糙率值与对应的初始糙率值之间的第一差异值。

步骤S2046，若第二目标数量大于第一目标数量，将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量小于第一目标数量，得到各断面分区的第二糙率率定值，第二目标数量用于表征第一差异值大于第一预设值的断面分区的数量。

示例性地，本申请实施例对第一预设值和第一目标数量的具体内容不做限定，本领域技术人员可以根据需求确定。当第二目标数量大于第一目标数量，将各断面分区的第一糙率率定值作为初始糙率值，并返回基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值的步骤，对各断面分区的糙率进行多次率定，直至第二目标数量小于第一目标数量，得到各断面分区的第二糙率率定值。

在一些可选的实施方式中，断面分区是根据水库的监测站点分布信息对水库沿程断面进行划分得到的。

在一些可选的实施方式中，每个分区的糙率率定顺序通过如下步骤获取：

步骤a1，获取水库的沿程断面的水流特征信息以及每个断面分区的位置信息。示例性地，本申请实施例中，沿程断面的水流特征信息可以是沿程断面的水流流动方向。

步骤a2，基于水流特征信息确定每个断面分区的分区类型。示例性地，本申请实施例中，基于沿程断面的水流特征信息可以将各断面分区划分为河道区和湖泊区

步骤a3，基于每个分区的位置信息以及分区类型确定对应分区的糙率率定顺序。示例性地，本申请实施例中，断面分区的率定顺序可以为先进行河道区的糙率率定，基于每个分区的位置信息依次从上游向下游率定，然后进行湖泊区对应的断面分区的糙率率定，依次从坝前向上游率定。此外，也可以任意指定不参与糙率率定的断面分区。

本发明实施例提供的方法，在进行分段糙率率定时考虑了率定顺序的影响。大型河道型水库水流特性差异极大，根据水流特性可大致分为河道区和湖泊区，在进行糙率率定时必须分别考虑。河道区的水位主要由上游来流流量控制，因而水位的变化随流量传递是自上而下的。湖泊区水位主要受坝前水位控制，因而水位传递是自下而上的。这样的处理方式符合物理规律，在计算时糙率率定收敛速度更快、精度更高，也可以避免数值震荡问题，得到的率定结果更为准确。

在一些可选的实施方式中，上述步骤S2042包括：

步骤b1，将第1个断面分区的各糙率值对应的第二水位信息和对应第一水位信息进行误差计算，得到第一糙率值序列中各糙率值分别对应的第一误差计算结果。示例性地，本申请实施中，可以将第1个断面分区的各糙率值的第二水位信息与对应的第一水位信息之间的差异值作为第一误差计算结果。

步骤b2，基于第一误差计算结果，将第一糙率值序列中误差最小的糙率值作为第1个断面的第一糙率率定值。示例性地，将计算误差值最小的糙率值作为第一糙率率定值。

在一些可选的实施方式中，上述步骤b1包括：

步骤b11，基于第1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线。示例性地，第一水位信息可以是水位随时间的变化信息，基于第一水位信息可以确定第一水位变化曲线，第一水位变化曲线用于表征实测水位与时间之间的对应关系。

步骤b12，基于第1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第1个断面分区在每个糙率值条件下的第二水位变化曲线。示例性地，每个糙率值条件下计算得到的第二水位信息为对应的计算水位随时间的变化信息。

步骤b13，基于第1个断面分区的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第1个断面分区在每个糙率值条件下水位的相位误差，得到第一误差计算结果。示例性地，本申请实施例中，确定第1个断面分区在每个糙率值条件下水位的相位误差，在确定相位误差时，一般会确定每个糙率值条件下的第二水位变化曲线和第二水位变化曲线中的水位峰值点(若有多个，则同时记录，包括最大值和最小值)，确定出各峰值点出现的时刻，即峰现时刻，然后确定每个第二水位变化曲线的峰现时刻与第一变化曲线的峰现时刻的误差的平方，基于计算得到的误差平方确定的误差最小的糙率值，即可作为该断面分区的糙率的率定值，若率定值在确定糙率率定值的上下边界上且不在糙率的上下限上，重新确定该断面分区的第一糙率值序列，直到率定值不在糙率的上下限上，即完成该步在该断面分区的糙率率定，并将该值作为该断面分区的率定值。

对于坝前湖泊区，水位整体变幅主要受坝前水位控制，坝前水位波动传播时间主要受水深影响，糙率影响不大，且受波浪、风速等影响，水位采集信息不一定完全准确，会有一定误差，采用水位误差平方和最小进行糙率率定时，特别容易出现对于不同糙率值，水位的计算误差变化不大，这时候就难以确定相对准确的糙率值。而采用水位波动的相位误差控制，就可得到相对合理的取值，糙率率定的计算速度和稳定性更好，特别适用于坝前水位较高、坝前水位波动，湖泊区范围大的水库库区河段糙率的率定

在一些可选的实施方式中，上述步骤S2043包括：

步骤c1，分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果。

步骤c2，基于第二误差计算结果，将第一糙率值序列中误差最小的糙率值作为第N+1个断面分区的糙率率定值。

在一些可选的实施方式中，上述步骤c1包括：

步骤c11，基于第N+1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线。

步骤c12，基于第N+1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第N+1个断面分区在每个糙率值条件下的第二水位变化曲线。

步骤c13，基于第N+1个断面的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第N+1个断面在每个糙率值条件下水位的相位误差，得到第二误差计算结果。

下面通过具体的实施例对本发明提供的一种水库水动力模型糙率率定方法进行说明。

实施例1：

对于某一大型水库，库区长约600公里，沿程糙率未知，采用一维水动力模型计算时段为2020年6月20日至9月10日的洪水在整个库区的传播过程，即水位随时间变化过程。可同步获取该时段内上游来流、支流汇入、坝前水位及坝前出流信息。沿河道从上游到坝前共布置14个水位测站，其中2个站点位于河道区，12个站点位于湖泊区，各测站在该时段内的水位信息均可获得。根据已知信息率定沿程各河段(断面分区)的糙率，各河段的糙率值如表1所示。根据计算方法，确定糙率率定方案如下：

表1

糙率率定的步长在本计算中选定为0.0025。

首先确定A断面分区的计算工况，糙率初始取值为0.03，步长为0.0025，根据电脑性能，确定并行计算核心数为9，相应断面A糙率取值数列为：(0.0400、0.0375、0.0350、0.0325、0.0300、0.0275、0.0250、0.0225、0.0200)。

将各河段的糙率值及A河段的糙率序列分别代入水动力模型计算，得到相对应的9组计算结果，经比较，当河段糙率取为0.03时，水位误差的平方和最小，吻合度最高，因此可以确定A河段糙率值为0.03，并将其作为率定值进行下一河段计算。下一河段为B河段。

对于B河段，采用和A断面相似计算方式，经计算确定B河段糙率值为0.03；依次序对各河段逐步计算直到C河段，得到各河段在第一次循环后的糙率取值如下表2所示：

表2

经比较，第一次循环得到糙率与初始值差别较大，需要以第一次循环得到的糙率值为新的初始值，重新从A河段开始率定，进行第二轮循环，计算结果如下表3所示：

表3

经比较，第二轮循环得到结果与第一轮完全一致，各河段糙率率定完成。

率定后的沿程典型站点水位实测值与计算值对比见下图3、图4和图5所示。图3为A河段站点实测水位与糙率率定后的计算水位对比图，图4为I河段站点实测水位与糙率率定后的计算水位对比图，图5为P河段站点实测水位与糙率率定后的计算水位对比图。基于对比图可知，本率定方法得到的结果与实测值吻合较好。

实施例2：

对于某一大型水库，库区长约600公里，拟采用一维水动力模型计算其洪水在库区的传播过程，但沿程糙率未知，需要根据已知信息率定沿程各河段的糙率。拟研究时段为2021年9月1日至10月2日，可同步获取该时段内上游来流、支流汇入、坝前水位及坝前出流信息。沿河道从上游到坝前共布置14个水位测站，其中1个站点位于河道区，10个站点位于湖泊区，3个站点位于过渡区，各测站在该时段内的水位信息均可获得。根据计算方法，确定其糙率率定顺序如下表4所示，沿程各河段糙率初始值如下表5所示：

表4

表5

河段名称	A	B	C	D	E	F	G
								糙率取值	0.03	0.05	0.05	0.06	0.05	0.03	0.05
河段名称	H	I	J	K	L	M	N
								糙率取值	0.05	0.05	0.06	0.0825	0.06	0.05	0.05

糙率率定的步长在本计算中选定为0.0025。

首先确定A断面的计算工况，糙率初始取值为0.03，步长为0.0025，根据电脑性能，确定并行计算核心数为15，相应断面A糙率取值数列为(0.0475、0.0450、0.0425、0.0400、0.0375、0.0350、0.0325、0.0300、0.0275、0.0250、0.0225、0.0200、0.0175、0.0150、0.0125)。

将各河段的糙率值及A河段的糙率序列分别代入水动力模型计算，得到相对应的15组计算结果，经比较，当河段糙率取为0.03时，计算水位峰值出现时刻与实测水位峰值出现时刻相位差最小，吻合度最高，因此可以确定A河段糙率值为0.03，并将其作为率定值进行下一断面计算。下一断面为N河段，对于N河段，采用和A断面相似计算方式，经计算确定N河段糙率值为0.065；依次序对各河段逐步计算直到B河段，得到各河段在第一次循环中的糙率取值如下表6所示：

表6

河道名称	A	B	C	D	E	F	G
								糙率取值	0.0300	0.0500	0.0475	0.0600	0.0500	0.0400	0.0375
河道名称	H	I	J	K	L	M	N
								糙率取值	0.0525	0.0300	0.0625	0.0675	0.0775	0.0750	0.0650

经比较，第一次循环得到计算与初始值差别较大，需要从A河段开始重新率定，进行第二轮循环，计算结果如下表7所示：

表7

河道名称	A	B	C	D	E	F	G
								糙率取值	0.0300	0.0500	0.0475	0.0600	0.0500	0.0400	0.0375
河道名称	H	I	J	K	L	M	N
								糙率取值	0.0525	0.0300	0.0625	0.0675	0.0775	0.0750	0.0650

经比较，第二轮循环得到结果与第一轮完全一致。各河段糙率率定结束。

率定后的沿程各典型站点水位实测值与计算值进行对比可以确定，本率定方法得到的结果与实测值吻合较好，水位峰值出现时机完全对应，可准确预测洪水在库区的峰现时间，提高预报精度。

在本实施例中还提供了一种水库水动力模型糙率率定装置，该装置用于实现上述实施例及优选实施方式，已经进行过说明的不再赘述。如以下所使用的，术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

本实施例提供一种水库水动力模型糙率率定装置，如图6所示，包括：

获取模块601，用于获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，第一水位信息为实测水位信息；

第一确定模块602，用于基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值；

计算模块603，用于基于水位计算边界信息、水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算第一糙率值序列中各糙率值对应的第二水位信息；

第二确定模块604，用于根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值。

在一些可选的实施方式中，第二确定模块604包括：

排序子模块，用于获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序；

第一确定子模块，用于令N＝1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值；

第二确定子模块，用于判断N是否小于M，若N小于M，令N＝N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值；

第三确定子模块，用于返回判断N是否小于M的步骤，直至N＝M，得到各断面的第一糙率率定值。

在一些可选的实施方式中，第二确定模块604，还包括：

比对子模块，用于将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值；

第四确定子模块，用于若第二目标数量大于第一目标数量，

将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回基于预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量小于第一目标数量，得到各断面分区的第二糙率率定值，第二目标数量用于表征第一差异值大于第一预设值的断面分区的数量。

在一些可选的实施方式中，断面分区是根据水库的监测站点分布信息对水库沿程断面进行划分得到的。

在一些可选的实施方式中，每个分区的糙率率定顺序通过如下步骤获取：

获取水库的沿程断面的水流特征信息以及每个断面分区的位置信息；

基于水流特征信息确定每个断面分区的分区类型；

基于每个分区的位置信息以及分区类型确定对应分区的糙率率定顺序。

在一些可选的实施方式中，第一确定子模块包括：

第一误差计算单元，用于将第1个断面分区的各糙率值对应的第二水位信息和对应第一水位信息进行误差计算，得到第一糙率值序列中各糙率值分别对应的第一误差计算结果；

第一确定单元，用于基于第一误差计算结果，将第一糙率值序列中误差最小的糙率值作为第1个断面的第一糙率率定值。

在一些可选的实施方式中，第一误差计算单元包括：

第一确定子单元，用于基于第1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；

第二确定子单元，用于基于第1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第1个断面分区在每个糙率值条件下的第二水位变化曲线；

第三确定子单元，用于基于第1个断面分区的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第1个断面分区在每个糙率值条件下水位的相位误差，得到第一误差计算结果。

在一些可选的实施例中，第二确定子模块包括：

第二误差计算单元，用于分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果；

第二确定单元，用于基于第二误差计算结果，将第一糙率值序列中误差最小的糙率值作为第N+1个断面分区的糙率率定值。

在一些可选的实施方式中，第二误差计算单元包括：

第四确定子单元，用于基于第N+1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；

第五确定子单元，用于基于第N+1个断面分区在每个糙率值条件下计算得到的第二水位信息确定第N+1个断面分区在每个糙率值条件下的第二水位变化曲线；

第六确定子单元，用于基于第N+1个断面的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第N+1个断面在每个糙率值条件下水位的相位误差，得到第二误差计算结果。

上述各个模块和单元的更进一步的功能描述与上述对应实施例相同，在此不再赘述。

本实施例中的水库水动力模型糙率率定装置是以功能单元的形式来呈现，这里的单元是指ASIC(Application Specific Integrated Circuit，专用集成电路)电路，执行一个或多个软件或固定程序的处理器和存储器，和/或其他可以提供上述功能的器件。

本发明实施例还提供一种计算机设备，具有上述图6所示的水库水动力模型糙率率定装置。

请参阅图7，图7是本发明可选实施例提供的一种计算机设备的结构示意图，如图7所示，该计算机设备包括：一个或多个处理器10、存储器20，以及用于连接各部件的接口，包括高速接口和低速接口。各个部件利用不同的总线互相通信连接，并且可以被安装在公共主板上或者根据需要以其它方式安装。处理器可以对在计算机设备内执行的指令进行处理，包括存储在存储器中或者存储器上以在外部输入/输出装置(诸如，耦合至接口的显示设备)上显示GUI的图形信息的指令。在一些可选的实施方式中，若需要，可以将多个处理器和/或多条总线与多个存储器和多个存储器一起使用。同样，可以连接多个计算机设备，各个设备提供部分必要的操作(例如，作为服务器阵列、一组刀片式服务器、或者多处理器系统)。图7中以一个处理器10为例。

处理器10可以是中央处理器，网络处理器或其组合。其中，处理器10还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路，可编程逻辑器件或其组合。上述可编程逻辑器件可以是复杂可编程逻辑器件，现场可编程逻辑门阵列，通用阵列逻辑或其任意组合。

其中，存储器20存储有可由至少一个处理器10执行的指令，以使所述至少一个处理器10执行实现上述实施例示出的方法。

存储器20可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储根据计算机设备的使用所创建的数据等。此外，存储器20可以包括高速随机存取存储器，还可以包括非瞬时存储器，例如至少一个磁盘存储器件、闪存器件、或其他非瞬时固态存储器件。在一些可选的实施方式中，存储器20可选包括相对于处理器10远程设置的存储器，这些远程存储器可以通过网络连接至该计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

存储器20可以包括易失性存储器，例如，随机存取存储器；存储器也可以包括非易失性存储器，例如，快闪存储器，硬盘或固态硬盘；存储器20还可以包括上述种类的存储器的组合。

该计算机设备还包括通信接口30，用于该计算机设备与其他设备或通信网络通信。

本发明实施例还提供了一种计算机可读存储介质，上述根据本发明实施例的方法可在硬件、固件中实现，或者被实现为可记录在存储介质，或者被实现通过网络下载的原始存储在远程存储介质或非暂时机器可读存储介质中并将被存储在本地存储介质中的计算机代码，从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件的存储介质上的这样的软件处理。其中，存储介质可为磁碟、光盘、只读存储记忆体、随机存储记忆体、快闪存储器、硬盘或固态硬盘等；进一步地，存储介质还可以包括上述种类的存储器的组合。可以理解，计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件，当软件或计算机代码被计算机、处理器或硬件访问且执行时，实现上述实施例示出的方法。

虽然结合附图描述了本发明的实施例，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (10)

1.一种水库水动力模型糙率率定方法，其特征在于，所述方法包括：

获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，所述第一水位信息为实测水位信息；

基于所述预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，所述第一糙率值序列中包含不同的糙率值；

基于所述水位计算边界信息、所述水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算所述第一糙率值序列中各糙率值对应的第二水位信息；

根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值；

根据各糙率值对应的第二水位与第一水位的误差，确定对应断面分区的糙率率定值的步骤，包括：

获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序；

令N=1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值；

判断N是否小于M，若N小于M，令N=N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值；

返回判断N是否小于M的步骤，直至N=M，得到各断面的第一糙率率定值；

所述根据各糙率值对应的第二水位与第一水位的误差，确定对应断面分区的糙率率定值的步骤，还包括：

将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值；

若第二目标数量大于第一目标数量，

将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回所述基于所述预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，所述第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量小于所述第一目标数量，得到各断面分区的第二糙率率定值，所述第二目标数量用于表征所述第一差异值大于第一预设值的断面分区的数量；

基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值的步骤，包括：

将第1个断面分区的各糙率值对应的第二水位信息和对应第一水位信息进行误差计算，得到第一糙率值序列中各糙率值分别对应的第一误差计算结果；

基于所述第一误差计算结果，将所述第一糙率值序列中误差最小的糙率值作为第1个断面的第一糙率率定值。

2.根据权利要求1所述的方法，其特征在于，所述断面分区是根据水库的监测站点分布信息对水库沿程断面进行划分得到的。

3.根据权利要求1所述的方法，其特征在于，所述每个分区的糙率率定顺序通过如下步骤获取：

获取水库的沿程断面的水流特征信息以及每个断面分区的位置信息；

基于所述水流特征信息确定每个断面分区的分区类型；

基于所述每个分区的位置信息以及分区类型确定对应分区的糙率率定顺序。

4.根据权利要求1所述的方法，其特征在于，将第1个断面分区在每个糙率值条件下计算得到的第二水位变化信息和对应第一水位信息进行误差计算，得到第一误差计算结果的步骤，包括：

基于第1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；

基于第1个断面分区在每个糙率值条件下计算得到的第二水位信息确定所述第1个断面分区在每个糙率值条件下的第二水位变化曲线；

基于第1个断面分区的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第1个断面分区在每个糙率值条件下水位的相位误差，得到第一误差计算结果。

5.根据权利要求1所述的方法，其特征在于，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值的步骤，包括：

分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果；

基于第二误差计算结果，将所述第一糙率值序列中误差最小的糙率值作为第N+1个断面分区的糙率率定值。

6.根据权利要求5所述的方法，其特征在于，分别计算第N+1个断面的各糙率值分别对应的第二水位信息与对应的第一水位信息的误差，得到第一糙率值序列中各糙率值分别对应的第二误差计算结果的步骤，包括：

基于第N+1个断面分区的第一水位信息确定第1个断面分区的第一水位变化曲线；

基于所述第N+1个断面分区在每个糙率值条件下计算得到的第二水位信息确定所述第N+1个断面分区在每个糙率值条件下的第二水位变化曲线；

基于第N+1个断面的第一水位变化曲线以及对应每个糙率值条件下的第二水位变化曲线计算得到于第N+1个断面在每个糙率值条件下水位的相位误差，得到第二误差计算结果。

7.一种水库水动力模型糙率率定装置，其特征在于，所述装置包括：

获取模块，用于获取水库的非恒定流水动力模型、水库不同断面分区的初始糙率值、水位计算边界信息、不同断面分区第一水位信息、糙率率定范围以及预设糙率率定步长，所述第一水位信息为实测水位信息；

第一确定模块，用于基于所述预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，所述第一糙率值序列中包含不同的糙率值；

计算模块，用于基于所述水位计算边界信息、所述水库的非恒定流水动力模型以及各断面分区分别对应的第一糙率值序列，分别计算所述第一糙率值序列中各糙率值对应的第二水位信息；

第二确定模块，用于根据各糙率值对应的第二水位信息与第一水位信息的误差，确定对应断面分区的糙率率定值；

所述第二确定模块，包括：

排序子模块，用于获取每个断面分区的糙率率定顺序，按照率定顺序对各断面分区进行排序；

第一确定子模块，用于令N=1，基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值；

第二确定子模块，用于判断N是否小于M，若N小于M，令N=N+1，利用第1至N个断面的第一糙率率定值、第N+1个断面的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第N+1个断面进行糙率率定，得到第N+1个断面的第一糙率率定值；

第三确定子模块，用于返回判断N是否小于M的步骤，直至N=M，得到各断面的第一糙率率定值；

所述第二确定模块，还包括：

比对子模块，用于将每个断面分区的第一糙率率定值与初始糙率值进行比对，确定各对应断面分区的第一糙率率定值与初始糙率值之间的第一差异值；

第四确定子模块，用于若第二目标数量大于第一目标数量，

将各断面分区的第一糙率率定值分别作为各断面分区的初始糙率值，返回所述基于所述预设糙率率定步长、糙率率定范围以及水库不同断面分区的初始糙率值，确定各断面分区分别对应的第一糙率值序列，所述第一糙率值序列中包含不同的糙率值的步骤，直至第二目标数量数量小于所述第一目标数量，得到各断面分区的第二糙率率定值，所述第二目标数量用于表征所述第一差异值大于第一预设值的断面分区的数量；

所述基于第1个断面分区的第一糙率值序列、对应断面分区的第一水位信息、各糙率值对应的第二水位信息以及所述非恒定流水动力模型对第1个断面分区进行糙率率定，得到第1个断面分区的第一糙率率定值，包括：

将第1个断面分区的各糙率值对应的第二水位信息和对应第一水位信息进行误差计算，得到第一糙率值序列中各糙率值分别对应的第一误差计算结果；

基于所述第一误差计算结果，将所述第一糙率值序列中误差最小的糙率值作为第1个断面的第一糙率率定值。

8.根据权利要求7所述的装置，其特征在于，所述断面分区是根据水库的监测站点分布信息对水库沿程断面进行划分得到的。

9.一种计算机设备，其特征在于，包括：

存储器和处理器，所述存储器和所述处理器之间互相通信连接，所述存储器中存储有计算机指令，所述处理器通过执行所述计算机指令，从而执行权利要求1至6中任一项所述的水库水动力模型糙率率定方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机指令，所述计算机指令用于使计算机执行权利要求1至6中任一项所述的水库水动力模型糙率率定方法。