CN109558618B - 流域流量的获取方法、装置、设备及可读存储介质 - Google Patents

流域流量的获取方法、装置、设备及可读存储介质 Download PDF

Info

Publication number
CN109558618B
CN109558618B CN201810095189.5A CN201810095189A CN109558618B CN 109558618 B CN109558618 B CN 109558618B CN 201810095189 A CN201810095189 A CN 201810095189A CN 109558618 B CN109558618 B CN 109558618B
Authority
CN
China
Prior art keywords
basin
flow
water outlet
outlet section
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810095189.5A
Other languages
English (en)
Other versions
CN109558618A (zh
Inventor
龙笛
杜明达
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201810095189.5A priority Critical patent/CN109558618B/zh
Publication of CN109558618A publication Critical patent/CN109558618A/zh
Application granted granted Critical
Publication of CN109558618B publication Critical patent/CN109558618B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/52—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring the height of the fluid level due to the lifting power of the fluid flow

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明提供一种流域流量的获取方法、装置、设备及可读存储介质。该方法包括：获取待测流域的出水口断面的水面高程；根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域的出水口断面的水位和流量关系的流量相关参数。本发明通过获取待测流域的出水口断面的水面高程，并通过预设的水文模拟模型实现了对待测流域的流域流量进行测算，提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。

Description

流域流量的获取方法、装置、设备及可读存储介质
技术领域
本发明涉及水利领域，特别是涉及一种流域流量的获取方法、装置、设备及可读存储介质。
背景技术
随着科学技术的发展，为了更好的了解人类的生存环境，以使人类的生产活动更适应自然的规律，人们对于全球气候的预测需求越来越多。因此随着人类在水利领域和气象领域的深入研究，现有的观测站已经无法完全覆盖人们所需的观测范围，因而导致部分流域的水文数据无法通过观测站进行精确观测。
传统技术中，对于一些少资料流域，通常需要借用临近流域的水文模型来测算少资料流域的水文数据，其中，少资料流域指的是地面观测站点能够从该流域获取的数据量有限的区域。
然而，利用临近流域的水文模型测算少资料流域的水文数据时，其测算的结果误差较大，无法准确反映少资料流域的流域特性。
发明内容
基于此，有必要针对借用临近流域的水文模型测算少资料流域的水文数据所导致的误差较大的问题，提供一种流域流量的获取方法、装置、设备及可读存储介质。
第一方面，本发明实施例提供一种流域流量的获取方法，包括：
获取待测流域的出水口断面的水面高程；
根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域内的出水口断面的水位和流量关系的流量相关参数。
在其中一个实施例中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
在其中一个实施例中，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；
其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述的出水口断面水面高程的影响。
在其中一个实施例中，所述根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，还包括：
根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；
将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
在其中一个实施例中，所述根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值，包括：
执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程；
若所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件，则调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
在其中一个实施例中，所述判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件，包括：
根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量；
判断所述函数输出值是否位于预设的收敛范围内；
若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。
在其中一个实施例中，所述收敛函数包括：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量；所述所述
第二方面，本发明实施例提供一种流域流量的获取装置，包括：
获取模块，用于获取待测流域的出水口断面的水面高程；
处理模块，用于根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域内的出水口断面的水位和流量关系的流量相关参数。
在其中一个实施例中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的的出水口断面的水面高程。
在其中一个实施例中，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；
其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述出水口断面的水面高程的影响。
在其中一个实施例中，所述流域流量的获取装置，还包括：
确定模块，用于在所述处理模块根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；
替换模块，用于将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
在其中一个实施例中，所述确定模块，具体包括：
执行单元，用于执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程；
循环处理单元，用于在所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件时，调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并指示所述执行单元返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
在其中一个实施例中，所述执行单元，具体用于根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；并判断所述函数输出值是否位于预设的收敛范围内；若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量。
在其中一个实施例中，所述收敛函数包括：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量；所述所述
第三方面，本发明实施例提供的一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器用于执行所述计算机程序时实现上述任一实施例中所述方法的步骤。
第四方面，本发明实施例提供的一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一实施例中所述方法的步骤。
本发明提供的流域流量的获取方法、装置、设备以及可读存储介质，通过获取待测流域的出水口断面的水面高程，并根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，以确定所述待测流域的流域流量，由于该水文模拟模型中所用的参数包含了能够表征当前待测流域的流量相关参数以及表征当前待测流域的出水口断面的水底高程的参数，因此，该水文模拟模型是与当前待测流域的特性相匹配的模型，其相较于传统技术中采用临近流域的模型测算当前流域的水文数据，本实施例的水文模拟模型测算的水文数据准确度更高，进而利用该水文模拟模型计算流域流量的时候，所得到的待测流域的流域流量与待测流域的实际流域流量误差较小，其大大提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。
附图说明
图1为一个实施例提供的流域流量的获取方法流程示意图；
图2为另一个实施例提供的流域流量的获取方法流程示意图；
图3为又一个实施例提供的流域流量的获取方法流程示意图；
图4为又一个实施例所示的流域流量的获取方法流程示意图；
图5为又一个实施例所示的流域流量的获取方法流程示意图；
图6为一个实施例所示的流域流量的获取装置的结构示意图；
图7为另一个实施例所示的流域流量的获取装置的结构示意图；
图8为又一个实施例所示的流域流量的获取装置的结构示意图。
具体实施方式
随着科学技术的发展，人们针对水文领域的研究越来越深入。在水文领域中，由于观测站所获取的水文数据量有限，常常导致部分流域的观测数据较少，影响该区域的水文数据的测算。通常在研究过程中，人们通过利用临近流域的水文模型对少资料流域的水文数据进行测算。然而，该方法无法准确反映少资料流域的流域特性，进而对该地区的气候特征不能准确评估。
本发明实施例提供的流域流量的获取方法、装置、设备及可存储介质，通过所建立的精准水文模型对待测流域的流域流量进行测算，大大提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。
需要说明的是，本发明实施例提供的流域流量的获取方法，其执行主体可以是流域流量的获取装置，该装置可以通过软件、硬件或者软硬件结合的方式实现成为计算机设备的部分或者全部。可选的，该计算机设备可以为PC、便携式设备、服务器等具有数据处理功能的电子设备，本实施例对计算机设备的具体形式并不做限定。下述方法实施例的执行主体以计算机设备为例来进行说明。
为了使本发明的目的、技术方案及优点更加清楚明白，通过下述实施例并结合附图，对本发明实施例中的技术方案的进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
图1为一个实施例提供的流域流量的获取方法流程示意图，本实施例涉及的是计算机设备根据所获取的待测流域出水口断面的水面高程和预设的待测流域的水文模拟模型确定流域流量的具体过程。如图1所示，该方法包括：
S101、获取待测流域出水口断面的水面高程。
具体的，计算机设备可以通过测高卫星获取待测流域出水口断面的水面高程，该待测流域可以为由分水岭所包围的河流集水区。其中，出水口断面的水面高程为表征待测流域的出水口断面的水面高度的参数，该待测流域的出水口断面的水面高度可以为待测流域的水面的海拔高度，也可以为待测流域的水面相对于某个固定基准的高度，例如可以为相对于待测流域的地面观测站点的高度，对此本实施例不做限定。通常，选取与待测流域出水口最接近的测高点处进行出水口断面的水面高程的采集。
S102、根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域内的出水口断面的水位和流量关系的流量相关参数。
具体的，计算机设备将获取的出水口断面的水面高程作为预设的待测流域的水文模拟模型的参照数据，输入该水文模拟模型中，并通过该水文模拟模型的计算，输出相应的待测流域的流域流量。可选的，上述流量相关参数可以包括水位流量关系综合参数和流域出水口断面形态系数，可选的，还可以包括其他的参数。该水位流量关系综合系数可以用于表征待测流域的出水口断面的水位和流量之间的部分关系，上述流域出水口断面形态系数为表征待测流域的横截面形状的参数，该水位流量关系综合系数和流域出水口断面形态系数共同用于表征待测流域的出水口断面的水位和流量之间的关系。可选的，上述出水口断面的水底高程可以为获取到的水底的高度数据，该水底的高度数据可以为流域底部的海拔高度，当上述流域为河流时，上述出水口断面的水底高程可以表征该河流的河床的海拔高度。
本实施例中，由于上述水文模拟模型中所用的参数包含了能够表征当前待测流域的流量相关参数以及表征当前待测流域的出水口断面的水底高程的参数，因此，该水文模拟模型是与当前待测流域的特性相匹配的模型，其相较于于传统技术中采用临近流域的模型测算当前流域的水文数据，本实施例的水文模拟模型测算的水文数据准确度更高，进而利用该水文模拟模型计算流域流量的时候，所得到的待测流域的流域流量与待测流域的实际流域流量误差较小，其大大提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。
作为上述实施例的一种可选的实现方式，上述水文模拟模型可以为包含a(H-z)b的表达式；其中，a为水位流量关系综合参数、b为流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
本实现方式中，水文模拟模型可以是Qalt＝a(H-z)b，其中，Qalt为流域流量，或者可以包括a(H-z)b的表达式，本实施例对此并不做限定。其中，上述出水口断面的水面高程H，表征流域的水面的高度数据，该水面的高度数据可以为流域水面的海拔高度，当上述流域为河流时，上述出水口断面的水面高程可以表征该河流的河面的海拔高度。计算机设备可以根据获取到的待测流域的出水口断面的水面高程H，通过公式Qalt＝a(H-z)b得到待测流域的流域流量Qalt，或者还可以通过包含a(H-z)b的其他公式得到待测流域的流域流量Qalt。可选的，这里的“其他公式”例如可以是包含了能够计算冰冻圈的冰雪融化产生的流水量以及冰雪层所持有的液态水对流域流量的影响的公式，该部分公式可以作为流域流量Qalt的校正或者辅助公式，用于对特定气候圈中的流域流量进行更准确的测算。可选的，a的取值范围可以为0-250，b的取值范围可以为0.1-3，z的取值范围可以为卫星观测的河底的高度序列中的最低值到最低值减去10米的区间。
可选的，上述出水口断面的水面高程可以通过卫星测算得到，该卫星可以是测高卫星，还可以是其他集成了测高功能的卫星。例如，该测高卫星可以采用Jason-2卫星，也可以采用Jason-3卫星或者其他测高卫星，对此本实施例并不做限定，只要能够得到待测流域的出水口断面的水面高程即可。
可选的，上述待测流域的出水口断面的水面高程H可以满足H＝hsat-ρ-Δh的关系式，其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项。需要说明的是，该出水口断面的水面高程H可以直接通过H＝hsat-ρ-Δh获得，还可以通过包含hsat-ρ-Δh的其他公式获得，可选的，这里的“其他公式”可以是用于作为待测流域的出水口断面的水面高程H的校正或者辅助公式，例如可以为包含大气层中的云层、降雨以及降雪等气候现象对出水口断面的水面高程的影响。
针对上述H＝hsat-ρ-Δh，hsat可以通过测高卫星测量卫星与地面的距离得到，该距离hsat为测高卫星在测高瞬间的卫星距离地面测高点的高度；另外，卫星获取的重定后的星地距离ρ表征卫星到地面的距离，其可以通过测高卫星Jason-2的三级产品传感器地球物理数据(Sensor Geophys ica l data record，简称SGDR)获取，例如通过卫星的SGDR产品中的波形得到重定后的星地距离ρ。可选的，对星地距离进行重定可以为：计算机设备通过卫星发射具有特定波形的信号，该信号经由地面反射，再由卫星接收，然后，计算机设备根据接收到的反射信号的变化计算得到修正后的的卫星到地面的距离，即为重定后的星地距离。并且，上述环境修正项Δh为表征地表以及大气层中的对流层和电离层等空间因素对卫星测量数据的影响，通过该环境修正项Δh，对测高卫星获取的卫星高度数据进行补偿，可选的，该环境修正项Δh可以由数据产品SGDR提供。当使用Jason-2测高卫星时，由于Jason-2有10天的重访周期，因此测高卫星获取的数据为每隔十天一个的出水口断面的水面高程，通过对该间隔十天的出水口断面的水面高程按天进行线性插值，则可以得到出水口断面的水面高程的日连续序列。
由于出水口断面的水面高程H的关系式中包含了通过卫星产品对星地距离进行重定的修正，并且加入了地表以及大气层中的对流层和电离层等空间因素对卫星数据的修正，从而使得的出水口断面的水面高程H更为准确，进而提高了测算的流域流量的精度。
上述实现方式通过计算机设备将卫星获取的出水口断面的水面高程输入到包含表达式a(H-z)b的水文模拟模型中，从而结合了待测流域自身的水位流量关系综合参数a，流域出水口断面形态系数b以及待测流域的出水口断面的水底高程z，使得该水文模拟模型结合了待测流域自身的特性参数，其相较于传统技术中采用临近流域的水文模型测算当前流域的水文数据，该实现方式的水文模拟模型测算的数据准确度更高，进而利用该水文模拟模型计算流域流量的时候，所得到的待测流域的流域流量与待测流域的实际流域流量误差较小，其大大提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。
上述各实施例分别列举了计算机设备如何根据待测流域的驱动数据和预设的水文模拟模型得到待测流域的流域流量的过程，下面通过下述的实施例描述上述预设的水文模拟模型的训练过程。需要说明的是，以下方法仅用于解释本发明而不用于对本发明进行限定。
图2为另一个实施例提供的流域流量的获取方法流程示意图，本实施例涉及的是计算机设备如何确定上述水文模拟模型的具体过程。在上述施例的基础上，可选的，S101之前还可以包括：
S201、根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；所述流量相关参数包括水位流量关系综合参数和流域出水口断面形态系数。
需要说明的是，上述流量相关参数的初始设定值以及训练流域的出水口断面的水底高程的初始设定值为上述CREST模型设定的初始值，该初始设定值可以是参考以往模型的经验值，也可以是借鉴其他相关模型得到的，对此本实施例不做限定。
具体的，上述模型训练数据可以包括已知的输入数据和已知的输出数据，已知的输入数据可以为训练流域的出水口断面的水面高程，已知的输出数据可以为通过卫星或者地面观测站点实际观测到的流域流量，可选的，该实际观测到的流域流量可以为观测站点在连续三年所测得的数据中所选取的每一年中流域流量的最高值和最低值所形成的序列。该已知的输入数据和已知的输出数据之间具有一一对应关系。
当计算机设备得到模型训练数据之后，将模型训练数据中的训练流域的出水口断面的水面高程输入预设的CREST模型，通过该CREST模型进行计算，并输出待收敛的流域流量。进一步的，计算机设备通过将该CREST模型输出的待收敛的流域流量与实测的流域流量之间的相关性进行收敛限定，可以得到流量相关参数的实际值以及出水口断面的水底高程的实际值。其中，CREST模型中所包括的流量相关参数的初始设定值以及出水口断面的水底高程的初始设定值为上述CREST模型中的已知数据。需要说明的是，上述训练流域即为需要测算其流域流量的待测流域，且上述训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值的计算过程，其可能的实现方式可以参见下述图3和图4所示实施例的描述。
S202、将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
本实施例中，通过根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值，并且将CREST模型中的流量相关参数的初始设定值替换为流量相关参数的实际值，以及将出水口断面的水底高程的初始设定值替换为出水口断面的水底高程的实际值，从而使得该水文模拟模型结合了待测流域自身的特性参数，其相较于传统技术中采用临近流域的模型测算当前流域的水文数据，本实施例所得到的水文模拟模型测算的水文数据准确度更高，进而利用该水文模拟模型计算的流域流量，其与待测流域的实际流域流量误差较小，其大大提高了少资料流域的流域流量的测算精度，进而提高了少资料地区的气候估算的准确性。
图3为又一个实施例所示的流域流量的获取方法流程示意图，本实施例涉及的是计算机设备根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值的具体过程。在上述施例的基础上，可选的，如图3所示，S201具体可以包括：
S301、执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件。
S302、若所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件，则调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
具体的，计算机设备首先根据该CREST模型以及该模型中的流量相关参数的初始设定值以及训练流域的出水口断面的水底高程的初始设定值，通过公式Qalt'＝a'(H'-z')b'或者该公式的变形，从而确定待收敛的流域流量Qalt'，该训练流域的出水口断面的水面高程通常可以由测高卫星获取得到。然后，计算机设备可以进一步判断所确定的待收敛的流域流量Qalt'与地面站实测的流域流量之间的关系是否满足预设的收敛条件。若待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件，则将当前CREST模型中的流量相关参数的值以及出水口断面的水底高程的值，确定为上述训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；若待收敛的流域流量与实测的流域流量之间的关系不满足预设收敛条件，则调整CREST模型中的流量相关参数的值和出水口断面的水底高程的值，并返回执行上述处理操作，即再次根据公式Qalt'＝a'(H'-z')b'或者该公式的变形，确定新的待收敛的流域流量，并判断该新的待收敛的流域流量与实测的流域流量之间的关系是否满足预设收敛条件，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
可选的，上述S301中计算机设备判断待收敛的流域流量与实测的流域流量之间的关系是否满足预设收敛条件，可以通过图4所示的方法实现，具体为：
S401、根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量。
可选的，上述收敛函数可以是BIAS(Qsim，Qobs)和NSCE(Qsim，Qobs)，还可以是
和
其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量，所 述所述
S402、判断所述函数输出值是否位于预设的收敛范围内。
可选的，该收敛范围可以是满足预设条件数值范围，还可以是满足输出值的坐标落在收敛图形内的区域范围。
可选的，上述收敛函数为BIAS(Qsim，Qobs)和NSCE(Qsim，Qobs)时，其收敛范围可以设置为：BIAS(Qsim，Qobs)满足0±0.2之间的范围，以及NSCE(Qsim，Qobs)满足0.8～1之间的范围。
可选的，上述收敛函数为上述公式(1)和公式(2)时，其收敛范围可以设置为：BIAS(Qsim，Qobs)满足0±0.2之间的范围，以及OBJV满足0.8～1之间的范围。可选的，当需要测算的数据精度高时，可以将收敛范围设置的较小，例如OBJV的收敛范围可以设置为0.9～1之间，BIAS(Qsim，Qobs)的收敛范围可以设置为0±0.1之间，以满足高精度数据的测算需求；当需要测算的数据精度低时，可以将收敛范围设置的较大，例如OBJV的收敛范围可以设置为0.6～1之间，BIAS(Qsim，Qobs)的收敛范围可以设置为0±0.3之间，在满足数据精度的情况下提高数据测算的效率。
需要说明的是，OBJV为将Qsim与Qobs二者的纳什系数NSCE(Qsim，Qobs)和二者对数的纳什系数NSCE(logQsim，logQobs)取平均值，通过将logQsim和logQobs加入到收敛函数中，可以使得收敛函数能够表征Qsim与Qobs之间的相关性的同时，能够避免由于所获取的数据异常，例如获取到误差较大的数据时，使得二者相关性产生较大偏差的不良影响。
S403、若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。
为了便于理解上述收敛函数在训练过程中的具体用法，下面对该收敛函数进行举例说明：假设第一次根据预设的CREST模型和已知的模型训练数据中的出水口断面的水面高程H1得到第一次待收敛的流域流量Qalt'的值(设为A，该A与H1对应)，计算机设备将上述收敛函数中的Qsim替换为A，然后将Qobs替换为H1对应的卫星实测的流域流量，以计算得到上述收敛函数的值，同时判断上述收敛函数的值是否分别达到各自预设的收敛范围。如果没有达到，则调整上述CREST模型中的流量相关参数a、b以及出水口断面的水底高程z的值，继续通过调整后的CREST模型和H1得到第二次待收敛的流域流量Qalt'的值(设为B)，计算机设备将上述收敛函数中的Qsim替换为B，然后将Qobs替换为H1对应的卫星实测的流域流量，再次计算得到收敛函数的值，同时判断上述收敛函数的值是否分别达到各自的预设的收敛范围；若达到，则将上述调整后的流量相关参数以及出水口断面的水底高程的值分别作为流量相关参数的实际值以及出水口断面的水底高程的实际值，从而得到训练上述预设的水文模拟模型。
可选的，通常能够满足预设收敛条件的参数组合不止一组，则在调整流量相关参数的值以及出水口断面的水底高程的值时，利用多目标遗传算法，从算法收敛的结果，即帕累托前沿解集中选择一组参数组合作为最终结果，并重新计算与实际观测流量的系数进行验证。
本实施例中通过判断收敛函数是否满足预设的收敛范围，以确定当前模型中的流量相关参数以及出水口断面的水底高程能否满足预设的精度要求，从而确保训练的模型测算的流域流量的精度达到使用需求，进而保证了少资料地区的气候估算的准确性。
图5为又一个实施例所示的流域流量的获取方法流程示意图。如图5所示，该方法具体包括：
S501、根据已知的模型训练数据和预设的CREST模型Qalt'＝a'(H'-z')b'，执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，将该待收敛的流域流量Qalt'作为公式中的Qsim输入公式，并根据所述待收敛的流域流量和和确定函数输出值；判断函数输出值是否位于预设的收敛范围内；若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量。
S502、调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
S503、将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型，该水文模拟模型包含a(H-z)b的表达式。
其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
S504、根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量。
以上S501至S504的执行过程具体可以参见上述实施例的描述，其实现原理和技术效果类似，在此不再赘述。
图6为一个实施例提供的一种流域流量的获取装置的结构示意图。如图6所示，该装置包括：获取模块11和处理模块12。
具体的，获取模块11，用于获取待测流域的出水口断面的水面高程；
处理模块12，用于根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括用于表征待测流域内的流量相关参数以及用于表征所述待测流域的出水口断面的水底高程的参数。
本实施例提供的流域流量的获取装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
可选的，上述水文模拟模型可以满足Qalt＝a(H-z)b的关系式，也可以包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数、b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
可选的，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项。
图7为另一个实施例所示的流域流量的获取装置的结构示意图。在上述各实施例的基础上，如图7所示，上述装置还可以包括：确定模块13和替换模块14。
具体的，确定模块13，用于在所述处理模块12根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；所述流量相关参数包括水位流量关系综合参数和流域出水口断面形态系数。
替换模块14，用于将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
本实施例提供的流域流量的获取装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
图8为又一个实施例所示的流域流量的获取装置的结构示意图。在上述图7实施例的基础上，如图8所示，上述确定模块13具体可以包括：执行单元131和循环处理单元132。
具体的，执行单元131，用于执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程。
循环处理单元132，用于在所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件时，调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并指示所述执行单元返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
本实施例提供的流域流量的获取装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在其中的一个实施例中，上述执行单元131具体用于根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；并判断所述函数输出值是否位于预设的收敛范围内；若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量。
本实施例提供的流域流量的获取装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在其中的一个实施例中，上述收敛函数可以包括如下函数：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量，所述所述
本实施例提供的流域流量的获取装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在一个实施例中，还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序在被所述处理器执行时实现如下步骤：
获取待测流域的出水口断面的水面高程；
根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域内的出水口断面的水位和流量关系的流量相关参数。
在一个实施例中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
在一个实施例中，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述出水口断面的水面高程的影响。
在一个实施例中，处理器执行计算机程序时，在所述根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，还实现以下步骤：
根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程；若所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件，则调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量；判断所述函数输出值是否位于预设的收敛范围内；若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。
在一个实施例中，所述收敛函数包括：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量；所述所述
上述实施例提供的计算机设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，还提供了一种计算机可读存储介质，所述可读存储介质中存储有计算机程序，当所述计算机程序在计算机上运行时，使得计算机执行如下步骤：
获取待测流域的出水口断面的水面高程；
根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域内的出水口断面的水位和流量关系的流量相关参数。
在一个实施例中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
在一个实施例中，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述出水口断面的水面高程的影响。
在一个实施例中，处理器执行计算机程序时，在所述根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，还实现以下步骤：
根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
执行处理操作，所述处理操作包括：根据公式Qalt'＝a'(H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程；若所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件，则调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量；判断所述函数输出值是否位于预设的收敛范围内；若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。
在一个实施例中，，所述收敛函数包括：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量；所述所述
上述实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
上述实施例可以全部或者部分的通过软件、硬件、固件或者其任意组合来实现。当使用软件程序实现时，可以全部或者部分地以计算机程序产品的形式实现，所述计算机程序产品可以包括一个或者多个计算机指令。在计算机上加载和执行所述计算机程序指令时，全部或者部分地产生按照本发明实施例所述的流程或者功能。其中，所述计算机可以是通用的计算机、专用计算机、计算机网络或者其他可编程装置；所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输。例如，所述计算机指令可以从一个网站站点、计算机、服务器或者数据中心通过有线(例如同轴电缆、光纤数字用户线DSL)或者无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或者数据中心进行传输；所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或者多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，例如软盘、硬盘、磁带、光介质(例如DVD)、或者半导体介质(例如固态硬盘SSD)等。
本领域内的技术人员应明白，本申请的实施例可提供为方法、装置、计算机设备或可读存储介质等计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置、设备和方法，在没有超过本申请的范围内，可以通过其他的方式实现。例如，以上所描述的实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种流域流量的获取方法，其特征在于，包括：
获取待测流域的出水口断面的水面高程；
根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域的出水口断面的水位和流量关系的流量相关参数；
其中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
2.根据权利要求1所述的方法，其特征在于，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；
其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述出水口断面的水面高程的影响。
3.根据权利要求1-2任一项所述的方法，其特征在于，所述根据所述待测流域的出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量之前，还包括：
根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的的出水口断面的水底高程的实际值；其中，所述CREST模型中包括训练流域的流量相关参数的初始设定值和训练流域的出水口断面的水底高程的初始设定值；
将所述CREST模型中的流量相关参数的初始设定值替换为所述流量相关参数的实际值，以及将所述出水口断面的水底高程的初始设定值替换为所述出水口断面的水底高程的实际值，得到所述待测流域的水文模拟模型；其中，所述模型训练数据包括：训练流域的出水口断面的水面高程和实测的流域流量。
4.根据权利要求3所述的方法，其特征在于，所述根据已知的模型训练数据和CREST模型，确定训练流域的流量相关参数的实际值以及训练流域的出水口断面的水底高程的实际值，包括：
执行处理操作，所述处理操作包括：根据公式Qalt'＝a('H'-z')b'确定待收敛的流域流量Qalt'，并判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件；其中，H'为所述训练流域的出水口断面的水面高程；
若所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件，则调整所述CREST模型中的流量相关参数的值和所述CREST模型中的出水口断面的水底高程的值，并返回执行所述处理操作，直至所得到的新的待收敛的流域流量与实测的流域流量之间的关系满足预设收敛条件为止，并将调整后的流量相关参数的值作为所述流量相关参数的实际值以及将调整后的出水口断面的水底高程的值作为所述出水口断面的水底高程的实际值。
5.根据权利要求4所述的方法，其特征在于，所述判断所述待收敛的流域流量与所述实测的流域流量之间的关系是否满足预设收敛条件，包括：
根据所述待收敛的流域流量和收敛函数，确定所述待收敛的流域流量对应的函数输出值；其中，所述收敛函数的因变量包括所述待收敛的流域流量与所述实测的流域流量；
判断所述函数输出值是否位于预设的收敛范围内；
若否，则确定所述待收敛的流域流量与所述实测的流域流量之间的关系不满足预设收敛条件。
6.根据权利要求5所述的方法，其特征在于，所述收敛函数包括：和其中，Qsim为所述CREST模型输出的待收敛的流域流量，Qobs为实测的流域流量；所述所述
7.一种流域流量的获取装置，其特征在于，包括：
获取模块，用于获取待测流域的出水口断面的水面高程；
处理模块，用于根据所述出水口断面的水面高程和预设的待测流域的水文模拟模型，确定所述待测流域的流域流量；其中，所述水文模拟模型包括所述待测流域的出水口断面的水底高程、以及用于表征待测流域的出水口断面的水位和流量关系的流量相关参数；
其中，所述水文模拟模型包含a(H-z)b的表达式；其中，a为所述流量相关参数中的水位流量关系综合参数，b为所述流量相关参数中的流域出水口断面形态系数，z为所述待测流域的出水口断面的水底高程，H为所述待测流域的出水口断面的水面高程。
8.根据权利要求7所述的装置，其特征在于，所述待测流域的出水口断面的水面高程H满足H＝hsat-ρ-Δh的关系式；
其中，hsat为卫星高度，ρ为重定后的星地距离，Δh为环境修正项，所述环境修正项Δh用于表征所述待测流域对应的大气层对所述出水口断面的水面高程的影响。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器用于执行所述计算机程序时实现权利要求1-6中任一项方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-6中任一项方法的步骤。
CN201810095189.5A 2018-01-31 2018-01-31 流域流量的获取方法、装置、设备及可读存储介质 Active CN109558618B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810095189.5A CN109558618B (zh) 2018-01-31 2018-01-31 流域流量的获取方法、装置、设备及可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810095189.5A CN109558618B (zh) 2018-01-31 2018-01-31 流域流量的获取方法、装置、设备及可读存储介质

Publications (2)

Publication Number Publication Date
CN109558618A CN109558618A (zh) 2019-04-02
CN109558618B true CN109558618B (zh) 2019-09-20

Family

ID=65864473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810095189.5A Active CN109558618B (zh) 2018-01-31 2018-01-31 流域流量的获取方法、装置、设备及可读存储介质

Country Status (1)

Country Link
CN (1) CN109558618B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111061985B (zh) * 2019-11-05 2020-11-06 清华大学 无资料地区水库调节径流的计算方法、装置和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289570A (zh) * 2011-07-23 2011-12-21 浙江大学 基于降雨-径流-洪水演进计算的洪水预报方法
CN103093114A (zh) * 2013-02-05 2013-05-08 河海大学 一种基于地形和土壤特性的分布式流域缺水量测算方法
CN103473809A (zh) * 2013-09-27 2013-12-25 武汉大学 基于3d打印技术的流域水文模拟方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102034001A (zh) * 2010-12-16 2011-04-27 南京大学 一种以栅格为模拟单元的分布式水文模型设计方法
CN103106625B (zh) * 2013-03-08 2016-01-20 珠江水利委员会珠江水利科学研究院 水库、闸泵群联合抑咸调度方法
CN104408318A (zh) * 2014-12-03 2015-03-11 中国水利水电科学研究院 一种河流型水功能区季节性设计流量计算方法
CN106777775B (zh) * 2017-01-10 2020-02-14 清华大学 一种基于多断面水位预测河流流量的神经网络方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102289570A (zh) * 2011-07-23 2011-12-21 浙江大学 基于降雨-径流-洪水演进计算的洪水预报方法
CN103093114A (zh) * 2013-02-05 2013-05-08 河海大学 一种基于地形和土壤特性的分布式流域缺水量测算方法
CN103473809A (zh) * 2013-09-27 2013-12-25 武汉大学 基于3d打印技术的流域水文模拟方法

Also Published As

Publication number Publication date
CN109558618A (zh) 2019-04-02

Similar Documents

Publication Publication Date Title
Lauri et al. Using reanalysis and remotely sensed temperature and precipitation data for hydrological modeling in monsoon climate: Mekong River case study
Despotovic et al. Review and statistical analysis of different global solar radiation sunshine models
Anderson et al. Real-time hydraulic and hydrodynamic model of the St. Clair River, Lake St. Clair, Detroit river system
CN109543353B (zh) 三维水汽反演方法、装置、设备和计算机可读存储介质
Berre et al. Filtering of background error variances and correlations by local spatial averaging: A review
CN109035105A (zh) 一种月尺度蒸散发量的定量估算方法
CN104656108B (zh) 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法
Liang et al. Application of Bayesian model averaging approach to multimodel ensemble hydrologic forecasting
CN113848577A (zh) 一种基于动态分区的大规模gnss网并行解算方法及***
Zhao et al. High-precision ZTD model of altitude-related correction
CN102175213B (zh) Gps水准高程传递方法
Briggs et al. A glacial systems model configured for large ensemble analysis of Antarctic deglaciation
CN103473809B (zh) 基于3d打印技术的流域水文模拟方法
CN105842692A (zh) 一种insar测量中的大气校正方法
CN104040378B (zh) 气象预测装置以及气象预测方法
Shen et al. Long-term variability of the Kuroshio transport east of Taiwan and the climate it conveys
CN109558617A (zh) 径流数据的获取方法、装置、计算机设备及可读存储介质
CN105069295A (zh) 基于卡尔曼滤波的卫星以及地面降水测量值同化方法
CN115544785B (zh) 一种无资料梯级水库流域水文模拟方法和***
CN113009531A (zh) 一种小尺度高精度的低空对流层大气折射率模型
CN109558618B (zh) 流域流量的获取方法、装置、设备及可读存储介质
Malamataris et al. Integrated hydrological modelling of surface water and groundwater under climate change: the case of the Mygdonia basin in Greece
Yang et al. Space‐time evolution of Greenland Ice Sheet elevation and mass from Envisat and GRACE data
Jacob et al. Experimental and numerical study of the hydrodynamics of the western sector of Ria Formosa
Hoteit et al. Assessing 4D-VAR for dynamical mapping of coastal high-frequency radar in San Diego

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant