CN116306051B - 一种河湖关系水量交换信息确定方法、装置及电子设备 - Google Patents
一种河湖关系水量交换信息确定方法、装置及电子设备 Download PDFInfo
- Publication number
- CN116306051B CN116306051B CN202310588557.0A CN202310588557A CN116306051B CN 116306051 B CN116306051 B CN 116306051B CN 202310588557 A CN202310588557 A CN 202310588557A CN 116306051 B CN116306051 B CN 116306051B
- Authority
- CN
- China
- Prior art keywords
- target
- lake
- river
- water level
- period
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 321
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000008569 process Effects 0.000 claims abstract description 26
- 238000004088 simulation Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 16
- 238000011144 upstream manufacturing Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 19
- 230000015654 memory Effects 0.000 claims description 19
- 238000001556 precipitation Methods 0.000 claims description 17
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000005484 gravity Effects 0.000 claims description 10
- 230000008020 evaporation Effects 0.000 claims description 7
- 238000001704 evaporation Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 description 7
- 230000008595 infiltration Effects 0.000 description 4
- 238000001764 infiltration Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000029142 excretion Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003204 osmotic effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种河湖关系水量交换信息确定方法、装置及电子设备,包括:获取目标模型和分布式水文模型;利用分布式水文模型对水流在目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;将目标时段内目标断面的河道径流量信息输入到目标模型中,得到目标时段内的目标断面的水位;利用分布式水文模型对与目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据模拟结果计算目标时段内的湖泊水位;根据目标时段内的湖泊水位以及目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息。本发明提供的方法,解决了现有技术中无法准确评估河湖关系的问题。
Description
技术领域
本发明涉及水循环模拟技术领域,具体涉及一种河湖关系水量交换信息确定方法、装置及电子设备。
背景技术
对于一个典型的河湖系统来讲,当河流水位较高,顶托和倒灌作用可使得湖泊水位升高,使得汛期湖泊防洪压力增加,非汛期湖区旱情可得到缓解;湖区水位较高时,湖水大量涌入河流,一方面可以缓解湖泊防洪压力,同时也可导致湖泊的枯水期延长,旱情加剧。因此,定量揭示通湖泊与河流之间的水量相互作用关系,对流域防洪抗旱和经济社会可持续发展具有重要意义。
传统河湖关系的评估方法多是从局部入手,基于实测数据和水动力模型对河湖关系进行评估,无法实现考虑多影响因素,使得评估结果不够准确。
发明内容
因此,本发明要解决的技术问题在于克服现有河湖关系的评估不够准确的缺陷,从而提供一种河湖关系水量交换信息确定方法、装置及电子设备。
第一方面,本发明实施例公开了一种河湖关系水量交换信息确定方法,所述方法包括:获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息;利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位;利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位;根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息。
本发明提供的河湖关系水量交换信息确定方法,基于分布式水文模型确定目标河流与目标湖泊之间的水流量交换信息,与现有技术相比,能够从整个流域水文情势出发,更加全面系统地模拟评估河湖关系,同时基于分布式水文模型的径流预测能力,能够实现未来情景下河湖关系的预测评估,解决了现有技术中无法准确评估河湖关系的问题。
结合第一方面,在第一方面的一种可能的实施方式中,所述获取目标模型之前,所述方法还包括:获取目标河流的中目标断面的历史水位和流量数据;根据所述目标断面的历史水位和流量数据建立目标模型。
本实施方式提供的方法,通过目标断面的历史水位和流量数据建立目标模型,便于后续通过目标模型计算目标断面水位。
结合第一方面,在第一方面的一种可能的实施方式中,所述利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位,包括:基于所述分布式水文模型计算目标时段与所述目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量;获取所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息、目标时段湖泊向下游河流流出的水量以及人工干预的水量信息;根据所述目标时段与所述目标湖泊相连的上游河流流入湖泊的水量、目标时段内湖泊与含水层之间的总渗透流量、所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息以及人工干预的水量信息,计算目标时段内的湖泊水位。
本实施方式提供的方法,通过分布式水文模型确定目标时段与所述目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量,并根据目标时段与所述目标湖泊相连的上游河流流入湖泊的水量、目标时段内目标湖泊向下游河流流出的水量、目标湖泊的历史水位、目标时段内目标湖泊接受的降水量信息以及人工干预的水量信息,计算目标时段内的湖泊水位,在计算过程中有效考虑了影响湖泊水位的多个因素,使得计算得到的湖泊水位更加准确。
结合第一方面,在第一方面的一种可能的实施方式中,计算目标时段内的湖泊水位包括:
根据第一关系式计算所述目标时段内的湖泊水位,所述第一关系式为:
其中,hl,n为目标时段的湖泊水位,hl,n-1为前一时段湖泊的水位;为目标时段的时间步长;/>为目标时段湖泊接受的降水量;/>为目标时段湖泊的蒸发量;/>为目标时段湖泊的人工干预水量;/>为目标时段与湖泊相连的上游河流流入湖泊的水量;/>为目标时段湖泊向下游河流流出的水量;As,n-1表示前一时段湖泊水面面积;/>表示目标时段内湖泊与含水层之间的总渗流量。
本实施方式提供的方法,在计算过程中有效考虑了影响湖泊水位的多个因素,使得计算得到的湖泊水位更加准确。
结合第一方面,在第一方面的一种可能的实施方式中,所述根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息,包括:将所述目标时段内的湖泊水位与所述目标时段内目标断面的水位进行比对;根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息。
本实施方式提供的方法,基于目标时段内的湖泊水位与目标时段内目标断面的水位的比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息,使得计算得到的目标河流与目标湖泊之间的水流量交换信息更加准确。
结合第一方面,在第一方面的一种可能的实施方式中,根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息,包括:
当所述目标时段内的湖泊水位小于所述目标时段内目标断面的水位,根据第二关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第二关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度,hr为河流水位,hl,n为湖泊水位;
当所述目标时段内的湖泊水位大于所述目标时段内目标断面的水位,通过第三关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第三关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度;hr为河流水位;hl,n为湖泊水位。
本实施方式提供的方法,基于目标时段内的湖泊水位与目标时段内目标断面的水位的比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息,使得计算得到的目标河流与目标湖泊之间的水流量交换信息更加准确。
第二方面,本发明实施例还公开了一种河湖关系水量交换信息确定装置,所述装置包括:第一获取模块,用于获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息;第一确定模块,用于利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;第二确定模块,用于将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位;计算模块,用于利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位;第三确定模块,用于根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息。
结合第二方面,在第二方面的一种可能的实施方式中,所述装置还包括;第二获取模块,用于获取目标河流的中目标断面的历史水位和流量数据;建立模块,用于根据所述目标断面的历史水位和流量数据建立目标模型。
第三方面,本发明实施例还公开了一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如第一方面或第一方面任一可选实施方式所述的河湖关系水量交换信息确定方法。
第四方面,本发明实施方式还公开了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如第一方面或第一方面任一可选实施方式所述的河湖关系水量交换信息确定方法。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中河湖关系水量交换信息确定方法的一个具体示例的流程图;
图2为本发明实施例中河湖关系水量交换信息确定装置的一个具体示例的原理框图;
图3为本发明实施例中电子设备的一个具体示例图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
正如本申请背景技术所述,传统河湖关系的评估方法多是从局部入手,基于实测数据和水动力模型,无法实现考虑多影响因素,更无法对未来情景下河湖关系进行预测,因此,本发明提供了一种河湖关系水量交换信息确定方法,可以应用于任一处理器,或集成有处理器的电子设备,本方法基于分布式水文模型,与传统方法相比,能够从整个流域水文情势出发,更加全面系统地模拟评估河湖关系,定量揭示了湖泊与河流之间的水量相互作用关系,对流域防洪抗旱和经济社会可持续发展具有重要意义。
本发明实施例公开了一种河湖关系水量交换信息确定方法,如图1所示,该方法包括如下步骤:
步骤101,获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息。
示例性地,目标模型是用于表征目标河流的中目标断面的水位和流量关系的模型,目标河流可以是任一河湖系统中的河流,目标断面可以是目标河流中的任一断面,本申请实施例中,目标模型可以是一个数学模型,可以表征目标断面处水位和流量的关系;分布式用于模拟河湖系统中的水资源的变化信息,本申请实施例中,分布式水文模型是进行区域产汇流模拟的有效工具,其输入数据包括基本空间数据和水循环驱动数据。基本空间数据包括目标区域的DEM数据、土地利用数据、土壤类型数据;水循环驱动数据主要包括气象数据和取用水数据。
步骤102,利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息。
示例性地,本申请实施例中,分布式水文模型对目标河流中的流动过程进行模拟过程中,首先模拟各个计算单元的地表产流/入渗过程,产流量将形成地表径流进入河道,然后通过对水量在河道中的演进过程进行模拟得到不同断面的河道径流量。其中地表产流/入渗计算采用Green-Ampt(入渗模型)方法,指定断面处的河道径流量的计算采用马斯京根法,计算公式可以如下式(1)所示:
(1)
其中:Qout,2为计算时段末(目标时段末)目标断面处的河道径流量;Qin,1为计算时段初河道首端入口流量;Qin,2为计算时段末河道首端入口流量;Qout,1为计算时段初目标断面处的河道径流量;C1,C2,C3为马斯京根法的计算系数,对应计算公式如下式(2)所示:
(2)
其中,△t为时间步长;K为河道单元的存储时间常数;X为权重因子,它的值是楔形水体的函数,对于河流,取值在0~0.3之间,平均值约为0.2。
步骤103,将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位。
示例性地,本申请实施例中,目标时段内目标断面的河道径流量信息可以为Qout,2,目标模型可以用如下式(3)的关系式表示:
(3)
其中,hr表示河道指定断面处的水位;f(Qout,2)表示目标断面处的水位流量函数关系,Qout,2可以是输入的指定断面处的河道径流量信息。
步骤104,利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位。
示例性地,本申请实施例中,分布式水文模型可以实现对与目标河流所连接的目标湖泊的湖泊水量变化过程的模拟操作,根据模拟结果可以计算得到目标河流上游流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗流量,根据计算结果并结合目标湖泊的降水量、目标湖泊向目标河流下游流出的水量以及蒸发量等已知数据,可以准确计算目标时段内的湖泊水位。
步骤105,根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息。
示例性地,根据目标时段内的湖泊水位以及目标时段内目标断面的水位可以确定水流是从湖泊流向河流还是从河流流向湖泊,进而可以确定目标时段内目标河流与目标湖泊之间的水流量交换信息。
本发明提供的河湖关系水量交换信息确定方法,基于分布式水文模型确定目标河流与目标湖泊之间的水流量交换信息,与现有技术相比,能够从整个流域水文情势出发,更加全面系统地模拟评估河湖关系,同时基于分布式水文模型的径流预测能力,能够实现未来情景下河湖关系的预测评估,解决了现有技术中无法准确评估河湖关系的问题。
作为本发明一个可选实施方式,步骤101之前,该方法还包括:
获取目标河流的中目标断面的历史水位和流量数据。
示例性地,本申请实施例中,目标断面的历史水位和流量数据可以通过历史数据库得到,历史数据库中可以存储有目标河流中多个断面的历史水位和流量数据。
根据所述目标断面的历史水位和流量数据建立目标模型。
示例性地,根据目标断面的历史水位和流量数据之间的对应关系,建立表征目标断面水位和流量关系的数学模型。
作为本发明一个可选实施方式,步骤104,包括:
基于所述分布式水文模型计算目标时段与所述目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量。
示例性地,基于分布式水文模型可以计算目标时段与目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量;分布式水文模型的输入数据为基本的空间数据和水循环驱动数据,通过产流汇流过程的模拟计算可以得到上游河道末端出流量以及目标时段内湖泊与含水层之间的总渗透流量。
获取所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息、为目标时段湖泊向下游河流流出的水量以及人工干预的水量信息。
示例性地,目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息、为目标时段湖泊向下游河流流出的水量以及人工干预的水量信息都为已知输入数据,可以通过历史数据库或者相关资料获取。
根据所述目标时段与所述目标湖泊相连的上游河流流入湖泊的水量、目标时段内湖泊与含水层之间的总渗透流量、所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息以及人工干预的水量信息,计算目标时段内的湖泊水位。
示例性地,通过目标时段与目标湖泊相连的上游河流流入湖泊的水量、目标时段内湖泊与含水层之间的总渗透流量、目标湖泊的历史水位、目标时段内目标湖泊接受的降水量信息以及人工干预的水量信息,可以准确计算目标时段内的湖泊水位。
作为本发明一个可选实施方式,计算目标时段内的湖泊水位包括:
根据第一关系式计算所述目标时段内的湖泊水位,所述第一关系式为:
(4)
其中,hl,n为目标时段的湖泊水位,hl,n-1为前一时段湖泊的水位;为目标时段的时间步长;/>为目标时段湖泊接受的降水量;/>为目标时段湖泊的蒸发量;/>为目标时段湖泊的人工干预水量;/>为目标时段与湖泊相连的上游河流流入湖泊的水量;/>为目标时段湖泊向下游河流流出的水量;As,n-1表示前一时段湖泊水面面积;/>表示目标时段内湖泊与含水层之间的总渗流量。
示例性地,上式(4)中,hl,n为目标时段的湖泊水位(L),hl,n-1为前一时段湖泊的水位(L);为目标时段的时间步长(T);/>为目标时段湖泊接受的降水量(L3/T),为已知输入数据;/>为目标时段湖泊的蒸发量(L3/T),可基于湖泊上的参考作物腾发量和水面蒸发修正因子计算得到;/>为目标时段湖泊的人工干预水量(L3/T),包括人工补水或抽排量,正表示抽排,负表示补水,为已知输入数据;/>为目标时段与湖泊相连的上游河流流入湖泊的水量(L3/T),为上游河道的末端出流量;/>为目标时段湖泊向下游河流流出的水量(L3/T),当下游河道有实测数据时,通过实测值进行下泄;As,n-1表示前一时段湖泊水面面积(L2),是湖泊水位的函数;/>表示目标时段内湖泊与含水层之间的总渗流量(L3/T),正表示湖泊向含水层净排泄,负表示含水层向湖泊净排泄,可基于分布式水文模型计算得到。
作为本发明一个可选实施方式,步骤105,包括:
将所述目标时段内的湖泊水位与所述目标时段内目标断面的水位进行比对。
示例性地,本申请实施例中,比对过程可以包括但不限于对目标时段内的湖泊水位与目标时段内目标断面的水位之间大小关系的比对。
根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息。
示例性地,根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息。
作为本发明一个可选实施方式,根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息,包括:
当所述目标时段内的湖泊水位小于所述目标时段内目标断面的水位,根据第二关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第二关系式为:
(5)
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度,hr为河流水位,hl,n为湖泊水位;
当所述目标时段内的湖泊水位大于所述目标时段内目标断面的水位,通过第三关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第三关系式为:
(6)
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度;hr为河流水位;hl,n为湖泊水位。
示例性地,当目标时段内的湖泊水位小于所述目标时段内目标断面的水位时,水流从河流流向湖泊,根据第二关系式(5)计算从河流流向湖泊的水流流量;当目标时段内的湖泊水位大于所述目标时段内目标断面的水位时,水流从湖泊流向河流,根据第三关系式(6)计算从湖泊流向河流的水流流量;
本发明实施例还公开了一种河湖关系水量交换信息确定装置,如图2所示,该装置包括:第一获取模块201,用于获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息;第一确定模块202,用于利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;第二确定模块203,用于将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位;计算模块204,用于利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位;第三确定模块205,用于根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息。
本发明提供的河湖关系水量交换信息确定装置,基于分布式水文模型确定目标河流与目标湖泊之间的水流量交换信息,与现有技术相比,能够从整个流域水文情势出发,更加全面系统地模拟评估河湖关系,同时基于分布式水文模型的径流预测能力,能够实现未来情景下河湖关系的预测评估,解决了现有技术中无法准确评估河湖关系的问题。
作为本发明一个可选实施方式,该装置还包括;第二获取模块,用于获取目标河流的中目标断面的历史水位和流量数据;建立模块,用于根据所述目标断面的历史水位和流量数据建立目标模型。
作为本发明一个可选实施方式,计算模块,包括:第一计算子模块,用于基于所述分布式水文模型计算目标时段与所述目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量;获取子模块,用于获取所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息、目标时段湖泊向下游河流流出的水量以及人工干预的水量信息;第二计算子模块,用于根据所述目标时段与所述目标湖泊相连的上游河流流入湖泊的水量、目标时段内湖泊与含水层之间的总渗透流量、所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息以及人工干预的水量信息,计算目标时段内的湖泊水位。
作为本发明一个可选实施方式,第二计算子模块,包括:第三计算子模块,用于根据第一关系式计算所述目标时段内的湖泊水位,所述第一关系式为:
其中,hl,n为目标时段的湖泊水位,hl,n-1为前一时段湖泊的水位;:目标时段的时间步长;/>为目标时段湖泊接受的降水量;/>为目标时段湖泊的蒸发量;/>为目标时段湖泊的人工干预水量;/>:目标时段与湖泊相连的上游河流流入湖泊的水量;/>为目标时段湖泊向下游河流流出的水量;As,n-1:前一时段湖泊水面面积;/>表示目标时段内湖泊与含水层之间的总渗流量。
作为本发明一个可选实施方式,第三确定模块,包括:比对子模块,用于将所述目标时段内的湖泊水位与所述目标时段内目标断面的水位进行比对;第四计算子模块,用于根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息。
作为本发明一个可选实施方式,第四计算子模块,包括:第五计算子模块,用于当所述目标时段内的湖泊水位小于所述目标时段内目标断面的水位,根据第二关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第二关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度,hr为河流水位,hl,n为湖泊水位;
第六计算子模块,用于当所述目标时段内的湖泊水位大于所述目标时段内目标断面的水位,通过第三关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第三关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度;hr为河流水位;hl,n为湖泊水位。
本发明实施例还提供了一种电子设备,如图3所示,该电子设备可以包括处理器401和存储器402,其中处理器401和存储器402可以通过总线或者其他方式连接,图3中以通过总线连接为例。
处理器401可以为中央处理器(Central Processing Unit,CPU)。处理器401还可以为其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片,或者上述各类芯片的组合。
存储器402作为一种非暂态计算机可读存储介质,可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块,如本发明实施例中的河湖关系水量交换信息确定方法对应的程序指令/模块。处理器401通过运行存储在存储器402中的非暂态软件程序、指令以及模块,从而执行处理器的各种功能应用以及数据处理,即实现上述方法实施例中的河湖关系水量交换信息确定方法。
存储器402可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储处理器401所创建的数据等。此外,存储器402可以包括高速随机存取存储器,还可以包括非暂态存储器,例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中,存储器402可选包括相对于处理器401远程设置的存储器,这些远程存储器可以通过网络连接至处理器401。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个模块存储在所述存储器402中,当被所述处理器401执行时,执行如图1所示实施例中的河湖关系水量交换信息确定方法。
上述电子设备具体细节可以对应参阅图1所示的实施例中对应的相关描述和效果进行理解,此处不再赘述。
本领域技术人员可以理解,实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)、随机存储记忆体(Random AccessMemory,RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive,缩写:HDD)或固态硬盘(Solid-State Drive,SSD)等;所述存储介质还可以包括上述种类的存储器的组合。
虽然结合附图描述了本发明的实施例,但是本领域技术人员可以在不脱离本发明的精神和范围的情况下作出各种修改和变型,这样的修改和变型均落入由所附权利要求所限定的范围之内。
Claims (8)
1.一种河湖关系水量交换信息确定方法,其特征在于,所述方法包括:
获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息;
利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;
将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位;
利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位;
根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息;
所述根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息,包括:
将所述目标时段内的湖泊水位与所述目标时段内目标断面的水位进行比对;
根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息;
根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息,包括:
当所述目标时段内的湖泊水位小于所述目标时段内目标断面的水位,根据第二关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第二关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度,hr为河流水位,hl,n为湖泊水位;
当所述目标时段内的湖泊水位大于所述目标时段内目标断面的水位,通过第三关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第三关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度;hr为河流水位;hl,n为湖泊水位。
2.根据权利要求1所述的方法,其特征在于,所述获取目标模型之前,所述方法还包括:
获取目标河流的中目标断面的历史水位和流量数据;
根据所述目标断面的历史水位和流量数据建立目标模型。
3.根据权利要求1所述的方法,其特征在于,所述利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位,包括:
基于所述分布式水文模型计算目标时段与所述目标湖泊相连的上游河流流入湖泊的水量以及目标时段内湖泊与含水层之间的总渗透流量;
获取所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息、目标时段湖泊向下游河流流出的水量以及人工干预的水量信息;
根据所述目标时段与所述目标湖泊相连的上游河流流入湖泊的水量、目标时段内湖泊与含水层之间的总渗透流量、所述目标湖泊的历史水位、目标时段内所述目标湖泊接受的降水量信息以及人工干预的水量信息,计算目标时段内的湖泊水位。
4.根据权利要求3所述的方法,其特征在于,计算目标时段内的湖泊水位包括:
根据第一关系式计算所述目标时段内的湖泊水位,所述第一关系式为:
其中,hl,n为目标时段的湖泊水位,hl,n-1为前一时段湖泊的水位;为目标时段的时间步长;/>为目标时段湖泊接受的降水量;/>为目标时段湖泊的蒸发量;/>为目标时段湖泊的人工干预水量;/>为目标时段与湖泊相连的上游河流流入湖泊的水量;/>为目标时段湖泊向下游河流流出的水量;As,n-1表示前一时段湖泊水面面积;/> 表示目标时段内湖泊与含水层之间的总渗流量。
5.一种河湖关系水量交换信息确定装置,其特征在于,所述装置包括:
第一获取模块,用于获取目标模型和分布式水文模型,所述目标模型用于表征目标河流的中目标断面的水位和流量关系,所述分布式水文模型用于模拟河湖系统中的水资源的变化信息;
第一确定模块,用于利用所述分布式水文模型对水流在所述目标河流中的流动过程进行模拟,确定目标时段内目标断面的河道径流量信息;
第二确定模块,用于将所述目标时段内目标断面的河道径流量信息输入到所述目标模型中,得到所述目标时段内的目标断面的水位;
计算模块,用于利用所述分布式水文模型对与所述目标河流所连接的目标湖泊的湖泊水量变化过程进行模拟,得到模拟结果,并根据所述模拟结果计算目标时段内的湖泊水位;
第三确定模块,用于根据所述目标时段内的湖泊水位以及所述目标时段内目标断面的水位确定目标时段内目标河流与目标湖泊之间的水流量交换信息;
第三确定模块,包括:
比对子模块,用于将所述目标时段内的湖泊水位与所述目标时段内目标断面的水位进行比对;
第四计算子模块,用于根据比对结果计算目标时段内目标河流与目标湖泊之间的水流量交换信息;
第四计算子模块,包括:
第五计算子模块,用于当所述目标时段内的湖泊水位小于所述目标时段内目标断面的水位,根据第二关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第二关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度,hr为河流水位,hl,n为湖泊水位;
第六计算子模块,用于当所述目标时段内的湖泊水位大于所述目标时段内目标断面的水位,通过第三关系式计算所述目标时段内目标河流与目标湖泊之间的水流量交换信息,所述第三关系式为:
其中,Q e为从目标河流流向目标湖泊的水流流量,为工程参数固定值;g为重力加速度;hr为河流水位;hl,n为湖泊水位。
6.根据权利要求5所述的装置,其特征在于,所述装置还包括;
第二获取模块,用于获取目标河流的中目标断面的历史水位和流量数据;
建立模块,用于根据所述目标断面的历史水位和流量数据建立目标模型。
7.一种电子设备,其特征在于,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如权利要求1-4任一所述的河湖关系水量交换信息确定方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1-4中任一项所述的河湖关系水量交换信息确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310588557.0A CN116306051B (zh) | 2023-05-24 | 2023-05-24 | 一种河湖关系水量交换信息确定方法、装置及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310588557.0A CN116306051B (zh) | 2023-05-24 | 2023-05-24 | 一种河湖关系水量交换信息确定方法、装置及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116306051A CN116306051A (zh) | 2023-06-23 |
CN116306051B true CN116306051B (zh) | 2023-08-18 |
Family
ID=86829134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310588557.0A Active CN116306051B (zh) | 2023-05-24 | 2023-05-24 | 一种河湖关系水量交换信息确定方法、装置及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116306051B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06117888A (ja) * | 1992-10-09 | 1994-04-28 | Fujitsu Ltd | 堰の放流量演算方法 |
CN109858666A (zh) * | 2018-12-12 | 2019-06-07 | 中国水利水电科学研究院 | 一种湖泊水资源量评估与预测方法 |
CN113106916A (zh) * | 2021-04-15 | 2021-07-13 | 河海大学 | 一种基于湖泊水位-流量与水位-面积关系的江湖交汇河道江湖交互作用量化方法 |
-
2023
- 2023-05-24 CN CN202310588557.0A patent/CN116306051B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06117888A (ja) * | 1992-10-09 | 1994-04-28 | Fujitsu Ltd | 堰の放流量演算方法 |
CN109858666A (zh) * | 2018-12-12 | 2019-06-07 | 中国水利水电科学研究院 | 一种湖泊水资源量评估与预测方法 |
CN113106916A (zh) * | 2021-04-15 | 2021-07-13 | 河海大学 | 一种基于湖泊水位-流量与水位-面积关系的江湖交汇河道江湖交互作用量化方法 |
Non-Patent Citations (1)
Title |
---|
采煤沉陷区分布式水循环模型研究;李慧;中国博士学位论文全文数据库 工程科技Ⅱ辑;第3章 * |
Also Published As
Publication number | Publication date |
---|---|
CN116306051A (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114997541B (zh) | 一种基于数字孪生技术的城市内涝预测方法及预警平台 | |
CN114254526B (zh) | 一种湖泊水-气界面二氧化碳交换量的评估方法及系统 | |
CN112784502B (zh) | 一种水文水力学动态双向耦合的洪水预测方法 | |
CN112632871B (zh) | 基于遥感的无资料自由溢流水库出流过程动态估算方法 | |
CN115983483A (zh) | 耦合循环神经网络与河道水位流量关系反算建库前来水情景下河道水位的预测方法 | |
WO2018078674A1 (ja) | シミュレーション装置、シミュレーション方法及び記録媒体 | |
CN115375198A (zh) | 区域河湖水系连通联合调度和水质安全保障方法及系统 | |
CN112528563A (zh) | 一种基于svm算法的城市内涝预警方法 | |
CN111985129A (zh) | 一种城市暴雨内涝精细化模拟方法 | |
KR101758116B1 (ko) | 실측 유량 전이를 이용한 미계측 유역 장기유출량 산정방법 | |
CN113869804B (zh) | 一种洪涝灾害下的电网设备风险预警方法及系统 | |
CN115688246A (zh) | 一种局部坐标系下的水库库容模拟方法及装置 | |
CN115290148A (zh) | 一种复杂水流条件下的流量测量方法及系统 | |
CN115510771A (zh) | 一种基于水文水动力耦合模型的社区雨涝数值模拟方法 | |
CN116306051B (zh) | 一种河湖关系水量交换信息确定方法、装置及电子设备 | |
CN113946964A (zh) | 一种草原河流河道流量汇流计算方法 | |
CN103870699A (zh) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 | |
KR101703809B1 (ko) | 기지 유량 전이를 이용한 미계측 유역 장기유출량 산정방법 | |
CN112964295B (zh) | 一种基于5g技术的闸泵站智能测控终端 | |
CN115511291A (zh) | 基于内涝淹没水深监测的岩溶洼地过水能力分析方法 | |
KR20240000915A (ko) | 유역 상세화를 통한 실적자료 기반 물수급 분석 방법 및 그 방법을 구현하기 위한 컴퓨터로 판독 가능한 기록매체에 저장된 프로그램 | |
CN114861571A (zh) | 一种河道型水库动态边界计算方法、装置及存储介质 | |
Wang et al. | R-statistic based predictor variables selection and vine structure determination approach for stochastic streamflow generation considering temporal and spatial dependence | |
CN113987767A (zh) | 自动校准的洪水滚动预报方法 | |
Basri et al. | Hydrological Modelling of Surface Runoff for Temengor Reservoir Using GR4H Model. |
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 |