CN116702658A - 一种基于预见期和形状参数的流域单位线模拟方法及应用 - Google Patents

一种基于预见期和形状参数的流域单位线模拟方法及应用 Download PDF

Info

Publication number
CN116702658A
CN116702658A CN202310981626.4A CN202310981626A CN116702658A CN 116702658 A CN116702658 A CN 116702658A CN 202310981626 A CN202310981626 A CN 202310981626A CN 116702658 A CN116702658 A CN 116702658A
Authority
CN
China
Prior art keywords
unit line
period
parameter
time
foresight
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.)
Granted
Application number
CN202310981626.4A
Other languages
English (en)
Other versions
CN116702658B (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.)
Ministry Of Water Resources Information Center Ministry Of Water Resources Monitoring And Forecasting Center For Hydrology And Water Resources
Tianjin Research Institute for Water Transport Engineering MOT
Original Assignee
Ministry Of Water Resources Information Center Ministry Of Water Resources Monitoring And Forecasting Center For Hydrology And Water Resources
Tianjin Research Institute for Water Transport Engineering MOT
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 Ministry Of Water Resources Information Center Ministry Of Water Resources Monitoring And Forecasting Center For Hydrology And Water Resources, Tianjin Research Institute for Water Transport Engineering MOT filed Critical Ministry Of Water Resources Information Center Ministry Of Water Resources Monitoring And Forecasting Center For Hydrology And Water Resources
Priority to CN202310981626.4A priority Critical patent/CN116702658B/zh
Publication of CN116702658A publication Critical patent/CN116702658A/zh
Application granted granted Critical
Publication of CN116702658B publication Critical patent/CN116702658B/zh
Priority to JP2024083398A priority patent/JP7557732B1/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种基于预见期和形状参数的流域单位线模拟方法及应用,采用基于预见期和形状参数的瞬时单位线公式计算流域的瞬时单位线。本发明有益效果:本发明提出的基于预见期和形状参数的瞬时单位线公式无论参数Tp、C3取值如何都具有以下几个特征:与时间轴包围面积等于1;函数形态先增加后减小,为单峰函数;在t≥0的范围内,函数值均大于等于0。以上性质使函数可以天然符合单位线的几个基本定义;避免了传统推求单位线的诸多困难,可以快速率定参数,提高了本单位线方法的率定、计算效率,且相对于Nash瞬时单位线,参数Tp、C3具有明确物理意义,方便所提方法的参数率定和应用推广。

Description

一种基于预见期和形状参数的流域单位线模拟方法及应用
技术领域
本发明属于水利工程数值计算技术领域,尤其是涉及一种基于预见期和形状参数的流域单位线模拟方法及应用。
背景技术
流域出口的流量过程是洪水预报关注的重点对象,也是高精度河道水动力模型计算中的关键输入参数,利用流域单位线进行短期洪水预报是目前普遍采用的技术手段。流域单位线的定义可以表述为:对任意一个具体的流域,其在单位时段内均匀分布的单位净雨量所形成的流域出口站的地面径流过程线称为单位线。流域单位线是一个概念性黑箱模型,其应用符合倍比叠加原理。
但是传统意义上的流域单位线推求困难,往往根据单位线的定义设置许多约束条件(例如:单位线先增大后减小、单位线与时间轴包围的面积等于流域降水量、单位线上的值均大于等于零等等)利用分析法、试错法、最小二乘法、线性规划等最优化算法进行求解。即便如此,有些推求得到的单位线仍然存在锯齿形,不符合一般自然规律。另一类方法是应用Nash瞬时单位线公式,但是Nash瞬时单位线实质是将流域等效成n个具有相同调蓄作用的线性水库串联,这在数学上获得了成功,但在物理概念上十分牵强,也导致了Nash瞬时单位线的参数(即n和K)不具有直观的物理解释。
发明内容
有鉴于此,本发明旨在提出一种基于预见期和形状参数的流域单位线模拟方法及应用,本流域单位线模型简单、实用,且模型参数具有明确的物理意义,从而克服传统流域单位线推求困难、Nash瞬时单位线参数物理意义不明确的困难,给流域洪水预报提供解决方案。
为达到上述目的,本发明的技术方案是这样实现的:
一种基于预见期和形状参数的流域单位线模拟方法,包括基于预见期和形状参数的瞬时单位线公式,采用基于预见期和形状参数的瞬时单位线公式计算流域的瞬时单位线,基于预见期和形状参数的瞬时单位线公式为:
式中,为流域的瞬时单位线;/>为计算地面径流;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>表示指数函数;Gamma()表示伽马函数。
进一步的,包括以下步骤:
S1:收集流域水文资料;
S2:应用基于预见期和形状参数的瞬时单位线公式;
S3:将基于预见期和形状参数的瞬时单位线转换为时段单位线;
S4:根据流域历史资料率定时段单位线的参数Tp、C3
S5:应用推求的时段单位线进行短期洪水预报。
进一步的,在步骤S1中,收集单位线洪水预报所需的流域水文资料,所述流域水文资料包括流域面积F、历史洪水的时间序列、净雨量序列、实测地面径流量序列。
进一步的,在步骤S3中,将基于预见期和形状参数的瞬时单位线转换为时段单位线,包括以下步骤:
S31、计算瞬时单位线关于时间t的积分:
其中,S(t)为流域瞬时单位线关于时间的积分函数;为流域瞬时单位线;/>为积分的时间变量;
得到:
式中,S(t)为流域瞬时单位线关于时间的积分函数;C3为形状参数;Tp为预见期;双参数的Gamma函数为上不完全伽马函数;
S32、通过两个S曲线相减,得到时段单位线:
其中,为时段单位线;/>为时间间隔;S(t)为流域瞬时单位线关于时间的积分函数;/>为时间轴上延迟/>的S曲线;
得到:
式中,为时段单位线;Gamma()表示伽马函数;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>为时间间隔;max()为取最大值函数;
S33、特别的,对于流域面积为F,单位为km2,时间间隔为,单位为h,净雨为10mm的单位线:
式中,为10mm净雨对应的单位线;F为流域面积;/>为时间间隔;/>为时段单位线。
进一步的,在步骤S4中,根据流域历史资料率定时段单位线的参数Tp、C3,包括以下率定过程:
将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程;
将各时段净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水模拟结果;公式为:
其中,为第i个时段序号的历史洪水模拟结果;/>为第j个时段净雨量;/>为模拟的时段序号;/>为单位线计算的时段序号,且/>;/>为10mm净雨对应的单位线;/>为时间间隔。
进一步的,在步骤S4中,参数率定采用历史洪水模拟误差最小原则,确定单位线的参数Tp和C3,参数率定过程使用参数率定算法,参数率定算法为遍历法或试错法或SCE-UA法。
进一步的,在步骤S5中,应用推求的时段单位线进行短期洪水预报,包括以下内容:
应用推求的时段单位线进行短期洪水预报时,将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程;
将各时段净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水预报结果;公式为:
其中,为第i个时段序号的洪水预报结果;/>为第j个时段净雨量;/>为预报的时段序号;/>为单位线计算的时段序号,且/>;/>为10mm净雨对应的单位线;/>为时间间隔。
进一步的,一种电子设备,包括处理器以及与处理器通信连接,且用于存储所述处理器可执行指令的存储器,所述处理器用于执行上述的短期洪水预报的方法。
进一步的,一种服务器,包括至少一个处理器,以及与所述处理器通信连接的存储器,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述处理器执行,以使所述至少一个处理器执行上述的短期洪水预报的方法。
进一步的,一种计算机可读取存储介质,存储有计算机程序,所述计算机程序被处理器执行时实现上述的短期洪水预报的方法。
相对于现有技术,本发明所述的一种基于预见期和形状参数的流域单位线模拟方法及应用具有以下优势:
(1)本发明所述的一种基于预见期和形状参数的流域单位线模拟方法及应用,本发明提出的基于预见期和形状参数的瞬时单位线公式无论参数Tp、C3取值如何都具有以下几个特征:与时间轴包围面积等于1;函数形态先增加后减小,为单峰函数;在的范围内,函数值均大于等于0。以上性质使函数可以天然符合单位线的几个基本定义;避免了传统推求单位线的诸多困难,可以快速率定参数,提高了本单位线方法的率定、计算效率。
(2)本发明所述的一种基于预见期和形状参数的流域单位线模拟方法及应用,本方法的函数峰值出现在t=Tp处,即:Tp直观表示流域预见期。在实际应用中用很少的水文资料,就可以确定参数Tp的大致范围。
(3)本发明所述的一种基于预见期和形状参数的流域单位线模拟方法及应用,与Nash瞬时单位线比较,函数峰值出现时间仅由参数Tp决定,在Tp确定后函数形状由参数C3决定。物理意义明确,方便所提方法的参数率定和应用推广。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明实施例所述的基于预见期和形状参数的瞬时单位线公式的函数图形示意图;
图2为本发明实施例所述的当参数C3确定,所提瞬时单位线参数Tp与函数结果的关系示意图;
图3为本发明实施例所述的当参数Tp确定,所提瞬时单位线参数C3与函数结果的关系示意图;
图4为本发明实施例所述的两个S曲线相减,得到时段单位线的曲线示意图;
图5为本发明实施例所述的参数率定过程中评判指标NSE的等值线示意图;
图6为本发明实施例所述的参数率定后的最优结果示意图;
图7为本发明实施例所述的率定得到的时段单位线示意图;
图8为本发明实施例所述的应用推求的时段单位线进行短期洪水预报示意图;
图9为本发明实施例所述的短期洪水预报方法流程示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以通过具体情况理解上述术语在本发明中的具体含义。
下面将参考附图并结合实施例来详细说明本发明。
如图1至图9所示,一种基于预见期和形状参数的流域单位线模拟方法,包括基于预见期和形状参数的瞬时单位线公式,采用基于预见期和形状参数的瞬时单位线公式计算流域的瞬时单位线,基于预见期和形状参数的瞬时单位线公式为:
式中,为流域的瞬时单位线;/>为计算地面径流;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>表示指数函数;Gamma()表示伽马函数。
如图9所示,在本发明一种优选的实施方式中,提出一种短期洪水预报的方法,应用上述基于预见期和形状参数的流域单位线模拟方法,包括以下步骤:
S1:收集流域水文资料;
收集单位线洪水预报所需的流域水文资料,所述流域水文资料包括流域面积F、历史洪水的时间序列、净雨量序列、实测地面径流量序列。
S2:应用基于预见期和形状参数的瞬时单位线公式;
假设流域瞬时单位线可表示成:
式中为计算地面径流,Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>表示指数函数;Gamma()表示伽马函数。
S3:将基于预见期和形状参数的瞬时单位线转换为时段单位线;
计算瞬时单位线关于时间t的积分:
其中,S(t)为流域瞬时单位线关于时间的积分函数(也称S曲线);为流域瞬时单位线;/>为积分的时间变量;
得到:
式中,S(t)为流域瞬时单位线关于时间的积分函数(也称S曲线);C3为形状参数;Tp为预见期;双参数的Gamma函数为上不完全伽马函数。
通过两个S曲线相减,得到时段单位线:
其中,为时段单位线;/>为时间间隔;S(t)为流域瞬时单位线关于时间的积分函数(也称S曲线);/>为时间轴上延迟/>的S曲线;
得到:
式中,为时段单位线;Gamma()表示伽马函数;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>为时间间隔;max()为取最大值函数。
特别的,对于流域面积为F(km2),时间间隔为 (h),净雨为10(mm)的单位线:
式中,为10mm净雨对应的单位线;F为流域面积;/>为时间间隔;/>为时段单位线。
S4:根据流域历史资料率定时段单位线的参数(Tp、C3);
根据流域洪水预报的实际需求确定时段单位线的单位时段长,给定参数Tp、C3,即可运用本发明所提时段单位线进行历史洪水模拟。具体过程为:将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程,将各个净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水预报结果,公式为:
其中,为第i个时段序号的历史洪水模拟结果;/>为第j个时段净雨量(单位mm);为模拟的时段序号;/>为单位线计算的时段序号(/>);/>为10mm净雨对应的单位线;/>为时间间隔。
本方法的参数率定采用历史洪水模拟误差最小原则,确定单位线的参数Tp和C3,参数具体过程可以使用任意成熟的参数率定算法,推荐的算法有遍历法、试错法、SCE-UA法等。
S5:应用推求的时段单位线进行短期洪水预报;
在确定应时段单位线的单位时段长,率定得到最优的参数Tp、C3后,即可用推求的时段单位线进行短期洪水预报。预报计算时,将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程,将各时段净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水预报结果;公式为:
其中,为第i个时段序号的洪水预报结果;/>为第j个时段净雨量(单位mm);/>为预报的时段序号;/>为单位线计算的时段序号(/>);/>为10mm净雨对应的单位线;为时间间隔。
公式推导过程:
S3推导过程如下:
将瞬时单位线转换为时段单位线时,采用S曲线,按照S曲线的定义:
得到S曲线的表达式:
其积分后的表达式为:
式中双参数的Gamma函数为上不完全伽马函数。
通过两个S曲线相减,得到时段单位线:
得到:
式中max()为取最大值函数。
对于流域面积为F(km2),时间间隔为 (h),净雨为10(mm)的单位线:
关于不完全Gamma函数:
在Mathematica中Gamma()表示上不完全Gamma函数,而在Matlab中,gammainc(z,a)表示正规化的下不完全Gamma函数,本发明说明书中的不完全Gamma函数与数学软件Mathematica中的一致,在数学软件Matlab中尝试本方法的读者需注意此区别。
本发明的优势:
(1)本发明提出的基于预见期和形状参数的瞬时单位线公式无论参数Tp、C3取值如何都具有以下几个特征:与时间轴包围面积等于1;函数形态先增加后减小,为单峰函数;在的范围内,函数值均大于等于0。以上性质使函数可以天然符合单位线的几个基本定义。避免了传统推求单位线的诸多困难,可以快速率定参数,提高了本单位线方法的率定、计算效率。
(2)函数峰值出现在t=Tp处,即:Tp直观表示流域预见期。在实际应用中用很少的水文资料,就可以确定参数Tp的大致范围。
(3)与Nash瞬时单位线比较,函数峰值出现时间仅由参数Tp决定,在Tp确定后函数形状由参数C3决定。物理意义明确,方便所提方法的参数率定和应用推广。
实施例1
1)收集流域水文资料;
以中国水利水电出版社《水文预报(第4版)》P71页计算实例为例,收集到南河开峰谷站控制的流域历史洪水资料,包括流域面积、时间、时段净雨量、实测地面径流量,其中流域面积为5239.02km2,其他资料如下表所示。
2)应用基于预见期和形状参数的瞬时单位线公式;
实施例过程同S2描述完全一致,在此不再赘述。
为直观展示本发明所提公式,图1展示了基于预见期和形状参数的瞬时单位线公式的函数图形;图2展示了当参数C3确定,所提瞬时单位线参数Tp与函数结果的关系;图3展示了当参数Tp确定,所提瞬时单位线参数C3与函数结果的关系。
3)将基于预见期和形状参数的瞬时单位线转换为时段单位线;
实施例过程与S3描述完全一致,在此不再赘述。
为直观展示本发明所提方法,图4展示了两个S曲线相减,得到时段单位线的曲线图。
4)根据流域历史资料率定时段单位线的参数;
在本实施例中,根据收集到的资料,确定单位线洪水预报的单位时段长 =6h。使用遍历法对单位线的参数Tp和C3进行参数寻优,利用纳什效率系数(NSE)评价模拟过程的优劣,即:对给定的Tp、C3参数寻优区间、寻优步长,遍历所有可能的参数组合,计算每组参数组合条件下模拟结果的NSE值(/>),NSE值约接近1说明模拟结果越接近观测,取NSE最大的参数组合作为最优参数。
本实施例中,参数Tp的寻优区间是,寻优步长是1h;参数C3的寻优区间是/>,寻优步长是0.2。在以上区间所有参数组合进行遍历计算,得到的NSE等值线图如图5所示,最大NSE=0.94最优参数为:Tp=9,C3=1.1,最优参数时,历史资料模拟结果如图6所示。
5)应用推求的时段单位线进行短期洪水预报;
利用率定得到最优参数,确定时段单位线,本实施例中流域面积F=5239.02km2,时间间隔为=6h,采用普遍接受的净雨10mm单位线,得到单位线如图7所示。
假设该流域未来发生一场降雨,时段净雨量如下表所示。
利用时段单位线进行短期洪水预报,得到的结果如图8所示。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于预见期和形状参数的流域单位线模拟方法,其特征在于:包括基于预见期和形状参数的瞬时单位线公式,采用基于预见期和形状参数的瞬时单位线公式计算流域的瞬时单位线,基于预见期和形状参数的瞬时单位线公式为:
式中,为流域的瞬时单位线;/>为计算地面径流;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>表示指数函数;Gamma()表示伽马函数。
2.一种短期洪水预报的方法,应用权利要求1所述的基于预见期和形状参数的流域单位线模拟方法,其特征在于:包括以下步骤:
S1:收集流域水文资料;
S2:应用基于预见期和形状参数的瞬时单位线公式;
S3:将基于预见期和形状参数的瞬时单位线转换为时段单位线;
S4:根据流域历史资料率定时段单位线的参数Tp、C3
S5:应用推求的时段单位线进行短期洪水预报。
3.根据权利要求2所述的一种短期洪水预报的方法,其特征在于:在步骤S1中,收集单位线洪水预报所需的流域水文资料,所述流域水文资料包括流域面积F、历史洪水的时间序列、净雨量序列、实测地面径流量序列。
4.根据权利要求3所述的一种短期洪水预报的方法,其特征在于:在步骤S3中,将基于预见期和形状参数的瞬时单位线转换为时段单位线,包括以下步骤:
S31、计算瞬时单位线关于时间t的积分:
其中,S(t)为流域瞬时单位线关于时间的积分函数;为流域瞬时单位线;/>为积分的时间变量;
得到:
式中,S(t)为流域瞬时单位线关于时间的积分函数;C3为形状参数;Tp为预见期;双参数的Gamma函数为上不完全伽马函数;
S32、通过两个S曲线相减,得到时段单位线:
其中,为时段单位线;/>为时间间隔;S(t)为流域瞬时单位线关于时间的积分函数;/>为时间轴上延迟/>的S曲线;
得到:
式中,为时段单位线;Gamma()表示伽马函数;Tp、C3为瞬时单位线参数,其中Tp为预见期,C3为形状参数;/>为时间间隔;max()为取最大值函数;
S33、特别的,对于流域面积为F,单位为km2,时间间隔为,单位为h,净雨为10mm的单位线:
式中,为10mm净雨对应的单位线;F为流域面积;/>为时间间隔;/>为时段单位线。
5.根据权利要求4所述的一种短期洪水预报的方法,其特征在于:在步骤S4中,根据流域历史资料率定时段单位线的参数Tp、C3,包括以下率定过程:
将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程;
将各时段净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水模拟结果;公式为:
其中,为第i个时段序号的历史洪水模拟结果;/>为第j个时段净雨量;/>为模拟的时段序号;/>为单位线计算的时段序号,且/>;/>为10mm净雨对应的单位线;/>为时间间隔。
6.根据权利要求5所述的一种短期洪水预报的方法,其特征在于:在步骤S4中,参数率定采用历史洪水模拟误差最小原则,确定单位线的参数Tp和C3,参数率定过程使用参数率定算法,参数率定算法为遍历法或试错法或SCE-UA法。
7.根据权利要求5所述的一种短期洪水预报的方法,其特征在于:在步骤S5中,应用推求的时段单位线进行短期洪水预报,包括以下内容:
应用推求的时段单位线进行短期洪水预报时,将每个时段的净雨量除以10mm,再乘以净雨10mm的时段单位线得到对应这个时段净雨量的出流过程;
将各时段净雨量产生的出流过程在时间轴上按顺序前后排序,在时间轴上对各出流过程进行叠加得到洪水预报结果;公式为:
其中,为第i个时段序号的洪水预报结果;/>为第j个时段净雨量;/>为预报的时段序号;/>为单位线计算的时段序号,且/>;/>为10mm净雨对应的单位线;/>为时间间隔。
8.一种电子设备,包括处理器以及与处理器通信连接,且用于存储所述处理器可执行指令的存储器,其特征在于:所述处理器用于执行上述权利要求2-7任一所述的短期洪水预报的方法。
9.一种服务器,其特征在于:包括至少一个处理器,以及与所述处理器通信连接的存储器,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述处理器执行,以使所述至少一个处理器执行如权利要求2-7任一所述的短期洪水预报的方法。
10.一种计算机可读取存储介质,存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现权利要求2-7任一所述的短期洪水预报的方法。
CN202310981626.4A 2023-08-07 2023-08-07 一种基于预见期和形状参数的流域单位线模拟方法及应用 Active CN116702658B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202310981626.4A CN116702658B (zh) 2023-08-07 2023-08-07 一种基于预见期和形状参数的流域单位线模拟方法及应用
JP2024083398A JP7557732B1 (ja) 2023-08-07 2024-05-22 洪水ピークディレイタイム及び形状係数に基づく流域単位図のシミュレーション方法及び応用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310981626.4A CN116702658B (zh) 2023-08-07 2023-08-07 一种基于预见期和形状参数的流域单位线模拟方法及应用

Publications (2)

Publication Number Publication Date
CN116702658A true CN116702658A (zh) 2023-09-05
CN116702658B CN116702658B (zh) 2023-10-03

Family

ID=87826210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310981626.4A Active CN116702658B (zh) 2023-08-07 2023-08-07 一种基于预见期和形状参数的流域单位线模拟方法及应用

Country Status (2)

Country Link
JP (1) JP7557732B1 (zh)
CN (1) CN116702658B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108874936A (zh) * 2018-06-01 2018-11-23 河海大学 一种基于改进新安江模型的适用于山丘区的水文预报方法
CN111460686A (zh) * 2020-04-23 2020-07-28 中国水利水电科学研究院 一种大气、陆面与水文三者双向耦合方法
CN112329257A (zh) * 2020-11-19 2021-02-05 四川大学 适用于山区小流域暴雨山洪洪水预报预警的水文模型分段筛选方法
WO2022032872A1 (zh) * 2020-08-14 2022-02-17 贵州东方世纪科技股份有限公司 一种基于大数据的水文预报方法
CN115391712A (zh) * 2022-05-18 2022-11-25 河海大学 一种城市洪涝风险预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4323565B1 (ja) 2009-03-30 2009-09-02 学校法人 君が淵学園 崇城大学 降雨による河川氾濫予測情報を導出する端末及びプログラム
KR101978351B1 (ko) 2018-11-13 2019-05-15 주식회사 하이드로셈 Cctv 영상 기반의 실시간 자동 유량계측 시스템 및 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108874936A (zh) * 2018-06-01 2018-11-23 河海大学 一种基于改进新安江模型的适用于山丘区的水文预报方法
CN111460686A (zh) * 2020-04-23 2020-07-28 中国水利水电科学研究院 一种大气、陆面与水文三者双向耦合方法
WO2022032872A1 (zh) * 2020-08-14 2022-02-17 贵州东方世纪科技股份有限公司 一种基于大数据的水文预报方法
CN112329257A (zh) * 2020-11-19 2021-02-05 四川大学 适用于山区小流域暴雨山洪洪水预报预警的水文模型分段筛选方法
CN115391712A (zh) * 2022-05-18 2022-11-25 河海大学 一种城市洪涝风险预测方法

Also Published As

Publication number Publication date
CN116702658B (zh) 2023-10-03
JP7557732B1 (ja) 2024-09-30

Similar Documents

Publication Publication Date Title
JP7337361B2 (ja) 多点並行修正に基づく分布型水文モデルパラメータの較正方法
CN105678067B (zh) 城市河涌防洪排涝动态预报控制方法及系统
CN105887751B (zh) 一种考虑水库调洪演算的多变量洪水重现期计算方法
CN107633141B (zh) 一种一维水动力模型泵站特性曲线系数的辨识方法
CN106647289A (zh) 明渠输水调度闸门延时启闭时间的前馈控制方法
CN101899820B (zh) 面向河流生态系统保护的流域地表水可利用量确定方法
CN110377989B (zh) 基于水热耦合平衡的非一致性两变量设计洪水推求方法
CN110445130B (zh) 考虑最优无功支撑的静态电压稳定裕度计算装置
CN114548680B (zh) 城市暴雨洪水管理模型参数自动率定方法及系统
CN112966221A (zh) 一种南方丰水地区考核年份用水总量折算方法
CN116401327A (zh) 无资料地区中小流域设计暴雨洪水计算辅助系统
CN116702658B (zh) 一种基于预见期和形状参数的流域单位线模拟方法及应用
CN108763829A (zh) 基于水动力相似的区域降雨模拟系统降雨过程设计方法
CN103530700A (zh) 城区配电网饱和负荷综合预测方法
CN111709147A (zh) 一种基于水文失事机理的设计洪水地区组成计算方法
CN104239716A (zh) 一种基于参数偏差灵敏度的电网设备参数识别与估计方法
CN113869804A (zh) 一种洪涝灾害下的电网设备风险预警方法及系统
CN114840989B (zh) 一种栅格尺度考虑水利工程拦蓄的河道汇流演算方法
CN108824349B (zh) 基于水动力过程相似的模型流域单位线的构建方法
CN109165402B (zh) 一种求解通用幂函数形明渠水力最优断面的方法
CN116611361A (zh) 一种水库峰值流量抵达坝前历时计算方法及装置
CN116488155A (zh) 电量预测方法及其装置、计算机可读存储介质
CN106557614B (zh) 一种基于Halphen IB分布的洪水频率分析方法
CN111861014B (zh) 一种分布式单位线洪水预报方法
CN108920571A (zh) 基于分布式水文模型的室内水文实验模型比例尺确定方法

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