CN108681783A - 一种水库实时多目标随机优化调度和风险评估方法 - Google Patents
一种水库实时多目标随机优化调度和风险评估方法 Download PDFInfo
- Publication number
- CN108681783A CN108681783A CN201810293796.2A CN201810293796A CN108681783A CN 108681783 A CN108681783 A CN 108681783A CN 201810293796 A CN201810293796 A CN 201810293796A CN 108681783 A CN108681783 A CN 108681783A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- formula
- flood
- real
- risk
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 238000005457 optimization Methods 0.000 title claims abstract description 38
- 238000012502 risk assessment Methods 0.000 title claims abstract description 20
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000013210 evaluation model Methods 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 4
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000013459 approach Methods 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 6
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000002948 stochastic simulation Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 2
- 238000013508 migration Methods 0.000 abstract description 2
- 230000005012 migration Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 241001672694 Citrus reticulata Species 0.000 description 2
- 230000006854 communication Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000034994 death Effects 0.000 description 1
- 231100000517 death Toxicity 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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/0635—Risk analysis of enterprise or organisation activities
-
- 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
-
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision 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
技术领域
本发明涉及水库优化调度与风险决策方法,尤其涉及一种水库实时多目标随机优化调度和风险评估方法。
背景技术
流域洪水具有巨大的破坏强度,是对人类社会危害最大的自然灾害之一。洪水往往会造成严重的人员伤亡、经济损失和生态破坏。自2006年至2016年,全球范围内的洪水灾害共造成了61730名人员死亡和3750亿美元的经济损失。为了抵御洪水,我国建造了大量的水库大坝,截止目前,共建成水库约98000余座。作为防洪工程措施和非工程措施有机结合的典范,水库优化调度对于削减洪峰、减轻洪灾损失具有重要作用。
水文预报提供了关于未来水文形势的关键信息,对于指导水库实时优化调度具有十分重要的意义。现有的技术往往局限于确定性环境之中,即采用确定性的水文预报成果和确定性的水库优化调度模型,忽略了水文预报和优化调度结果的不确定性,使得水库调度决策普遍缺乏科学性和可靠性。实际上,由于流域水文模型的内在误差、不精确的初始和边界条件、以及对实际物理水文过程的概化,水文预报具有固有的不确定性。制约水库优化调度决策的主要挑战之一就在于调度所依据的水文预报信息的不确定性。不确定性信息并不利于简化决策过程,但忽略不确定性就是忽略现实的正确规律。因此,如何在水库优化调度中充分地考虑水文预报不确定性,追踪不确定性的来源及其迁移转化规律,建立水库实时多目标随机优化调度和风险评估方法是目前亟待解决的技术难题。
发明内容
发明目的:针对现有技术的不足,本发明综合考虑水文预报不确定性及其动态演进特性,结合随机优化思想,提出了一种水库实时多目标随机优化调度和风险评估方法。
技术方案:本发明提供了一种水库实时多目标随机优化调度和风险评估方法,包括以下步骤:
(1)建立水文预报演进模型,模拟水文预报不确定性的动态演进过程,合成入库水文集合预报并构建洪水情景树;
(2)建立基于实时修正策略的水库实时多目标随机优化调度模型;
(3)采用约束法和两步求解过程推求不确定环境下的帕累托前沿;
(4)建立风险评估模型,量化各调度指标的不确定性和调度风险;
(5)采用随机多属性决策模型从帕累托前沿上确定最理想的水库调度方案,并定量评估决策结果的不确定性;
(6)建立基于实时信息修正的滚动预报-调度-决策流程。
进一步的,所述步骤(1)中定义单时段的预报更新uw,t为相邻两个时刻之间水纹预报误差的差,uw,t计算公式如下:
式中,qt为t时刻的流量实测值;fw,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报,t=w+i,i=0,1,…,H,H为预见期;ew,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报误差;
假定面临时段(即当前时段t)的水文预报是完美的,即ft,t=qt,推求水文预报不确定性的动态演进过程,计算公式如下:
单时段的水文预报更新uw,w+i-1,i=1,…,H,构成如下所示的水文预报更新向量Uw,-:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1] (3);
采用以下方差-协方差矩阵VCV来度量水文预报更新向量的统计特性:
式中,vari为uw,w+i-1的方差,covi,j为uw,w+i-1和uw,w+j-1之间的协方差,i=1,…,H;
采用Cholesky(乔列斯基)分解将VCV转换成一个下三角矩阵和其转置的乘积,即VCV=V·VT,然后,采用以下公式对水文预报更新向量Uw,-进行随机模拟:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1]=[θ1,θ2,…θH]·VT (5);
式中,θi(i=1,2,…,H)为服从独立正态同分布的随机变量;
将水文预报更新向量Uw,-代入预报演进过程计算公式(2)中即可合成入库水文集合预报;
基于随机生成的水文集合预报,构建洪水情景树为:
其中,I为情景树规模;ωt i代表洪水情景树中的第i个洪水情景,即洪水情景树的第i个树节点;该洪水情景树由从当前时段w到预见期末w+H的一系列树节点构成。
进一步的,所述步骤(2)中水库实时多目标随机优化调度模型的第一个优化目标为水库期望最高水位最低,表达式为:
式中,t为调度时序;w为水文集合预报发布的时刻;Zt i为t时刻第i个洪水情景下的水库水位(m);P(ωt i)为第i个洪水情景出现的概率;
第二个目标为期望最大出库流量最小,表达式为:
式中,Qt i是t时刻第i个洪水情景下的水库出库流量(m3/s);
水库实时多目标随机优化调度模型包括以下约束条件:
1)水量平衡约束:
Vt i=Vt-1 i+(It i-Qt i)·△t (8);
式中,Vt i和Vt-1 i分别为第i个洪水情景下t-1时刻的水库蓄水量;It i和Qt i分别为t时刻第i个洪水情景下的水库入库流量和出库流量(m3/s);△t为调度时间间隔(s);
2)水库水位约束:
式中,Zt 和分别为t时刻第i个洪水情景下水库水位的最低和最高限值;
3)水库出库流量约束:
式中,Qt 和分别为t时刻第i个洪水情景下水库出库流量的最低和最高限值;
4)水库泄流能力约束:
Qt i≤Qmax(Zt i) (11);
式中,Qmax(Zt i)为t时刻第i个洪水情景下的水库水位对应的水库泄流能力;
5)水库出库流量变幅约束:
|Qt i-Qt-1 i|≤△Qm (12);
式中,△Qm为允许的水库出库流量最大变幅,该约束主要用于将相邻时段的水库出库流量限定在指定范围以内;
6)初始和边界条件:
Z1 i=ZINI;ZT i→ZFL (13);
式中,ZINI分别为ZFL初始水库水位和目标期末水位。
进一步的,所述步骤(3)包括:
(31)确定每个优化目标的上下界;先优化第一个目标而不考虑第二个目标,再优化第二个目标而不考虑第一个目标;然后,采用约束法将多目标问题转换为一系列单目标优化问题;选择第二个目标作为当前的单个优化目标,将第一个目标作为额外的约束条件加到约束集中;在上下界之间参数化摄动第一个目标值即可获得一系列非劣解集(记为K个);选取的决策变量为所有洪水情景下每个时段的水库出库流量;采用非线性规划求解上述转换后的单目标优化问题;期望最大出库流量为
(32)将步骤(31)中求得的期望最大出库流量作为固定的约束值,仅优化第一个目标,从而量化水库最高水位的不确定性;期望最高水位为
通过采用上述约束法结合两步求解过程,可以获得不确定环境下的帕累托前沿,并将水库最高水位和最大出库流量描述为服从一定分布的随机变量。
进一步的,所述步骤(4)风险评估模型包括水库漫坝风险和下游洪灾风险,将水库漫坝风险定义为最高水位超过安全水位的概率,表达式为:
式中,f(Zmax)为最高水位Zmax的概率密度函数;
将下游洪灾风险定义为最大出库流量超过安全泄量的概率,表达式为:
式中,f(Qmax)为最大出库流量Qmax的概率密度函数。
进一步的,所述步骤(5)中利用层次分析法将决策者的主观偏好信息量化为指标权重,首先从帕累托前沿的概率分布中随机抽样生成N个指标值样本;然后将该样本和指标权重作为TOPSIS方法的输入,通过N次蒙特卡罗迭代计算贴近度系数值;最后通过统计分析推求贴近度系数的概率分布,进而确定各方案的概率排序,建立随机多属性决策模型,从帕累托前沿上确定最理想的水库调度方案;采用排序不确定度RUD定量评估决策结果的整体不确定性程度,排序不确定度定义为:各方案获得除各自确定性排序以外所有排序的概率和,计算公式如下:
式中,bi r表示确定性排序为i的方案在随机因素的作用下获得排序r的概率,表示确定性排序为i的方案在随机因素的作用下获得排序si的概率;si为第i个方案的确定性排序。
进一步的,所述步骤(6)包括:
(61)随机生成水文集合预报并构建洪水情景树,作为水库实时多目标随机优化调度模型的输入;
(62)基于面临时段w发布的水文集合预报[fw,w i,fw,w+1 i,…,fw,T i](上标i指第i个洪水情景),采用约束法、两阶段求解过程和非线性规划相耦合的方法确定期望意义上的非劣解集
(63)对每个非劣解进行风险评估,计算水库漫坝风险以及下游洪灾风险;
(64)基于决策者提供的主观偏好信息,采用随机多属性决策方法从帕累托前沿上选择令决策者最为满意的期望非劣解
(65)对于所选择的期望决策序列仅实施面临时刻w的期望泄流策略后续时段的泄流策略将随着逐时段更新的水文集合预报而动态更新;
(66)随着时间推移至w+1时段,采用水量平衡公式更新水库状态,并将水文集合预报更新为[fw+1,w+1 i,fw+1,w+2 i,…,fw+1,T i]。
重复上述步骤(61)至步骤(66)直至水库调度期末,即获得每个时段的水库泄流策略和相应的风险信息。
有益效果:与现有技术相比,本发明具有以下优点和有益效果:
1、本发明方法在水库优化建模过程中明确地考虑了水文预报不确定性及其动态演进特性,并定量评估了不确定性沿预报-优化-决策过程链的发展、迁移和传播规律;
2、本发明方法适用于水库实时调度,能够提高调度决策的可靠性,为不确定性环境下的水库调度提供必要的风险信息和决策支撑。
附图说明
图1是本发明方法的流程图;
图2为不确定环境下的帕累托前沿示意图。
具体实施方式
下面通过实施例并结合附图,对本发明的技术方案进行详细说明。
如图1所示,一种水库实时多目标随机优化调度和风险评估方法,包括以下步骤:
(1)建立水文预报演进模型,模拟水文预报不确定性的动态演进过程,合成入库水文集合预报并构建洪水情景树。通过定义单时段的预报更新,推求水文预报不确定性演进公式,并基于预报更新的方差协方差矩阵随机模拟水文集合预报并据此构建入库洪水情景树,单时段的预报更新(从w-1时刻到w时刻)定义为相邻两个时刻之间预报误差的差。
按照预见期的长短划分,流域水文预报可以分为短期、中期和长期预报。短期水文预报主要依据实测的“落地”降雨信息,通常采用站点观测网实时获取(例如雨量站和雷达测雨等)。虽然短期水文预报比中长期水文预报更具可靠性,但其预见期有限,主要取决于流域汇流时间(一般为几小时至几天)。中期水文预报(预见期为2-15天)通常需要耦合数值天气预报模式,尤其适用于上游流量数据缺失或者在极端洪水条件下数据监测和传输设备损坏的情形。长期水文预报的预见期通常大于15天,一般采用多因子统计分析和遥相关的方法进行建模。尽管雷达测雨技术、数值天气预报模式、流域水文模型和水气遥相关研究取得了较大的进展和突破,但是不确定性始终存在于水文预报之中,给水库优化调度决策带来了较大的挑战。
中短期水文预报一般采取实时滚动预报与更新策略,因此,预报不确定性也随时间的推移呈现出动态演进规律。在一场洪水事件的开始阶段,未来时段水文预报的发布一般基于当前所获取的水文信息。随着时间的推移,未来某时刻的水文预报将根据不断更新的水文信息而同步更新。令fw,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报(t=w+i,i=0,1,…,H,H为预见期),令ew,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报误差;定义单时段的预报更新uw,t(从w-1时刻到w时刻)为相邻两个时刻之间水文预报误差的差,uw,t计算公式如下:
式中,qt为t时刻的流量实测值。上述公式表明:uw,t代表了fw,t(当前时段对qt预报)在fw-1,t(前一个时段对qt预报)基础上的更新。
假定面临时段的水文预报是完美的(即ft,t=qt),推求水文预报不确定性的演进过程,计算公式如下:
上述公式表明水文预报不确定性是时序相关的,并且呈现出随时间推移而动态演进的规律。
单时段的预报更新uw,w+i-1(i=1,…,H)构成如下所示的预报更新向量Uw,-:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1] (3);
本发明方法采用以下方差-协方差矩阵VCV来度量预报更新向量的统计特性:
式中,vari为uw,w+i-1的方差,covi,j为uw,w+i-1和uw,w+j-1之间的协方差,i=1,…,H。
当历史水文预报资料充分时,可以采用上述VCV矩阵分析水文预报不确定性特征,并量化不确定性的动态演进过程;此外,也可以根据上述公式合成水文集合预报,表征未来一定预见期内可能发生的各种洪水情景。由于VCV矩阵是一个半正定矩阵,可以采用Cholesky(乔列斯基)分解将其转换成一个下三角矩阵和其转置的乘积VCV=V·VT。然后,采用以下公式对预报更新向量Uw,-进行随机模拟:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1]=[θ1,θ2,…θH]·VT (5);
式中,θi(i=1,2,…,H)为服从独立正态同分布的随机变量。
根据历史水文预报资料样本进行统计分析,确定VCV矩阵的参数;采用上述公式随机模拟预报更新向量;将其代入预报演进过程计算公式(2)中即可合成入库水文集合预报。基于随机生成的水文集合预报,构建洪水情景树为:
其中,I为情景树规模;ωt i代表洪水情景树中的第i个洪水情景,即洪水情景树的第i个树节点;该洪水情景树由从当前时段w到预见期末w+H的一系列树节点构成。本步骤所构建的洪水情景树为水库多目标随机优化调度模型提供输入条件。
(2)建立基于实时修正策略的水库实时多目标随机优化调度模型。
汛期水库优化调度的目标是在保证水库大坝安全的前提下最小化水库出库流量。本发明方法考虑两个冲突目标:(1)大坝自身及其上游安全;(2)下游防洪控制断面安全。第一个调度目标希望尽可能地加大水库出库流量、降低水库水位从而以大坝防洪安全。相反地,第二个调度目标则希望尽可能多地拦蓄洪水、降低水库出库流量以避免下游防洪断面的洪灾损失。本发明方法采用洪水情景树来量化水文预报不确定性,建立基于实时修正策略的水库实时多目标随机优化调度模型,考虑的优化目标为:期望最高水位最低和期望最大出库流量最小。与现有技术不同的是,现有的水库优化调度模型均为确定性模型,而本发明方法构建的是考虑水文预报不确定性的随机优化调度模型。
第一个优化目标为水库期望最高水位最低,以保证水库自身和库区的防洪安全,表达式为:
式中,t为调度时序;w为水文集合预报发布的时刻;Zt i为t时刻第i个洪水情景下的水库水位(m);P(ωt i)为第i个洪水情景出现的概率;表明上述目标函数值与洪水情景树相关。Zt i为一个随机过程,与具体的洪水情景及其出现的概率相关。
第二个目标为期望最大出库流量最小,以保证下游防洪断面的安全,表达式为:
式中,Qt i是t时刻第i个洪水情景下的水库出库流量(m3/s);表明上述目标函数值与洪水情景树相关。同样地,Qt i也是一个随机过程,与具体的洪水情景及其出现的概率相关。
所建立的多目标随机优化调度模型包括以下约束条件:
1)水量平衡约束:
Vt i=Vt-1 i+(It i-Qt i)·△t (8);
式中,Vt i和Vt-1 i分别为第i个洪水情景下t-1时刻的水库蓄水量;It i和Qt i分别为t时刻第i个洪水情景下的水库入库流量和出库流量(m3/s);△t为调度时间间隔(s)。
2)水库水位约束:
式中,Zt 和分别为t时刻第i个洪水情景下水库水位的最低和最高限值。
3)水库出库流量约束:
式中,Qt 和分别为t时刻第i个洪水情景下水库出库流量的最低和最高限值。
4)水库泄流能力约束:
Qt i≤Qmax(Zt i) (11);
式中,Qmax(Zt i)为t时刻第i个洪水情景下的水库水位对应的水库泄流能力。
5)水库出库流量变幅约束:
|Qt i-Qt-1i|≤△Qm (12);
式中,△Qm为允许的水库出库流量最大变幅,该约束主要用于将相邻时段的水库出库流量限定在指定范围以内。
6)初始和边界条件:
Z1 i=ZINI;ZT i→ZFL (13);
式中,ZINI分别为ZFL初始水库水位和目标期末水位。
(3)采用约束法和两步求解过程推求不确定环境下的帕累托前沿。
如步骤(1)所述,本发明方法采用一个由一系列入流情景组成的洪水情景树来刻画水文预报的不确定性。与现有技术不同的是,为了应对水文预报的不确定性,本发明方法在步骤(2)中建立了多目标随机优化调度模型用于调度决策。为了增强对预报不确定性的适应性,该随机优化模型针对每一个水库入流情景制定一个不同的最优决策,以确保目标函数在平均意义上达到最优。因此,在考虑预报不确定性的情况下,每个非劣解将不再表现为帕累托前沿上的一个单点,而是使得两个目标函数值为服从一定统计分布的随机变量。不确定环境下的帕累托前沿示意图如图2所示,图中圆圈代表期望值,实线代表固定最高水位,虚线代表固定最大出库流量。
本发明提出约束法结合两步求解过程推求不确定环境下的帕累托前沿,进一步包括以下子步骤:
(31)首先,确定每个优化目标的上下界。
先优化第一个目标而不考虑第二个目标,再优化第二个目标而不考虑第一个目标。然后,采用约束法将多目标问题转换为一系列单目标优化问题。选择第二个目标作为当前的单个优化目标,将第一个目标作为额外的约束条件加到约束集中。在上下界之间参数化摄动第一个目标值即可获得一系列非劣解集。选取的决策变量为所有洪水情景下每个时段的水库出库流量。采用非线性规划求解上述转换后的单目标优化问题。在本步骤中,每个洪水情景下的水库最高水位等于固定的约束值,而最大出库流量则互不相同且呈现出一定的统计分布特性,期望最大出库流量为
(32)其次,将步骤(31)中求得的期望最大出库流量作为固定的约束值,仅优化第一个目标,从而量化水库最高水位的不确定性。在本步骤中,每个洪水情景下的水库最大出库流量等于固定的约束值。因此,期望最高水位为
通过采用上述约束法结合两步求解过程,可以获得不确定环境下的帕累托前沿,并将水库最高水位和最大出库流量描述为服从一定分布的随机变量。
(4)建立风险评估模型,量化各调度指标的不确定性和调度风险。
现有的风险评估方法主要基于固定的水库调度规则,未考虑优化调度。本发明方法与现有技术的不同之处在于:本发明直接将将水文预报不确定性融入水库优化调度过程之中。在所建立的水库多目标随机优化调度模型中,风险来源于不同洪水情景下最优调度策略的差异。为此,本步骤定量评估从水文预报到优化调度的风险传播过程,所构建的风险评估模型主要考虑两类风险,分别为:水库漫坝风险和下游洪灾风险。
风险Risk定义为荷载L超过抗力R的超越概率,即Risk=P(L>R)。对于水库漫坝风险而言,荷载为水库最高水位Zmax,抗力为水库安全水位Zs。本发明将水库漫坝风险定义为最高水位超过安全水位的概率,表达式为:
式中,f(Zmax)为最高水位Zmax的概率密度函数。
下游防洪断面的荷载为最大出库流量Qmax,抗力为安全泄量Qs。因此,本发明将下游洪灾风险定义为最大出库流量超过安全泄量的概率,表达式为:
式中,f(Qmax)为最大出库流量Qmax的概率密度函数。
(5)采用随机多属性决策模型从帕累托前沿上确定最理想的水库调度方案,并定量评估决策结果的不确定性。
利用层次分析法将决策者的主观偏好信息量化为指标权重。基于理想解法(TOPSIS,Technique for Order Preference by Similarity to an Ideal Solution)和蒙特卡罗模拟建立随机多属性决策模型,从帕累托前沿上确定最理想的水库调度方案。基本思路是首先从帕累托前沿的概率分布中随机抽样生成N个指标值样本;然后将该样本和指标权重作为TOPSIS方法的输入,通过N次蒙特卡罗迭代计算贴近度系数值;最后通过统计分析推求贴近度系数的概率分布,进而确定各方案的概率排序。本步骤中的层次分析法和TOPSIS方法属于已有技术,在本实施例中不作赘述。
在确定性环境下,TOPSIS方法依据各方案的确定性贴近度系数进行方案排序。然而,由于水文预报不确定性因素的影响,各方案的贴近度系数不再是常量,而是服从一定分布的随机变量。现有的技术方法一般只考虑各方案的确定性排序,没有对决策结果的不确定性进行定量评估。为此,本发明方法提出排序不确定度RUD,据此定量评估决策结果的整体不确定性程度。排序不确定度定义为:各方案获得除各自确定性排序以外所有排序的概率和,计算公式如下:
式中,bi r表示确定性排序为i的方案在随机因素的作用下获得排序r的概率,表示确定性排序为i的方案在随机因素的作用下获得排序si的概率;si为第i个方案的确定性排序。
(6)建立基于实时信息修正的滚动预报-调度-决策流程。
为了确定随机优化调度的实时修正策略,本发明方法将逐时段更新的水文预报与随机优化调度模型相结合,并明确地考虑预报不确定性的动态演进过程,建立基于实时信息修正的滚动预报-调度-决策流程。本步骤进一步包括以下子步骤,以确定每个时段的水库泄流策略和相应的风险:
(61)随机生成水文集合预报并构建洪水情景树,作为水库实时多目标随机优化调度模型的输入;
(62)基于面临时段w发布的水文集合预报[fw,w i,fw,w+1 i,…,fw,T i](上标i指第i个洪水情景),采用约束法、两阶段求解过程和非线性规划相耦合的方法确定期望意义上的非劣解集
(63)对每个非劣解进行风险评估,计算水库漫坝风险以及下游洪灾风险;
(64)基于决策者提供的主观偏好信息,采用随机多属性决策方法从帕累托前沿上选择令决策者最为满意的期望非劣解
(65)对于所选择的期望决策序列仅实施面临时刻w的期望泄流策略后续时段的泄流策略将随着逐时段更新的水文集合预报而动态更新;
(66)随着时间推移至w+1时段,采用水量平衡公式更新水库状态,并将水文集合预报更新为[fw+1,w+1 i,fw+1,w+2 i,…,fw+1,T i]。
实施上述技术方案可以得到不确定环境下的帕累托前沿、每个时段最优的水库泄流策略、水库漫坝风险和下游洪灾风险。与现有技术相比,上述技术方案在水库优化建模过程中明确地考虑了水文预报不确定性及其动态演进特性,并定量评估了不确定性沿预报-调度-决策过程链的动态传播和是时程演化规律,能够显著提高调度决策的可靠性,为不确定性环境下的水库调度提供决策支持。
以上详细描述了本发明的优选实施方式。但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种等同变换,这些等同变换均属于本发明的保护范围。
Claims (7)
1.一种水库实时多目标随机优化调度和风险评估方法,其特征在于,包括以下步骤:
(1)建立水文预报演进模型,模拟水文预报不确定性的动态演进过程,合成入库水文集合预报并构建洪水情景树;
(2)建立基于实时修正策略的水库实时多目标随机优化调度模型;
(3)采用约束法和两步求解过程推求不确定环境下的帕累托前沿;
(4)建立风险评估模型,量化各调度指标的不确定性和调度风险;
(5)采用随机多属性决策模型从帕累托前沿上确定最理想的水库调度方案,并定量评估决策结果的不确定性;
(6)建立基于实时信息修正的滚动预报-调度-决策流程。
2.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(1)中定义单时段的预报更新uw,t为相邻两个时刻之间水纹预报误差的差,uw,t计算公式如下:
式中,qt为t时刻的流量实测值;fw,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报,t=w+i,i=0,1,…,H,H为预见期;ew,t(m3/s)表示w时刻发布的关于未来时刻t的水文预报误差;
假定面临时段的水文预报是完美的,即ft,t=qt,推求水文预报不确定性的动态演进过程,计算公式如下:
单时段的水文预报更新uw,w+i-1,i=1,…,H,构成如下所示的水文预报更新向量Uw,-:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1] (3);
采用以下方差-协方差矩阵VCV来度量水文预报更新向量的统计特性:
式中,vari为uw,w+i-1的方差,covi,j为uw,w+i-1和uw,w+j-1之间的协方差,i=1,…,H;
采用Cholesky分解将VCV转换成一个下三角矩阵和其转置的乘积,即VCV=V·VT,然后,采用以下公式对水文预报更新向量Uw,-进行随机模拟:
Uw,-=[uw,w,uw,w+1,…,uw,w+H-1]=[θ1,θ2,…θH]·VT (5);
式中,θi(i=1,2,…,H)为服从独立正态同分布的随机变量;
将水文预报更新向量Uw,-代入预报演进过程计算公式(2)中即可合成入库水文集合预报;
基于随机生成的水文集合预报,构建洪水情景树为:
其中,I为情景树规模;ωt i代表洪水情景树中的第i个洪水情景,即洪水情景树的第i个树节点;该洪水情景树由从当前时段w到预见期末w+H的一系列树节点构成。
3.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(2)中水库实时多目标随机优化调度模型的第一个优化目标为水库期望最高水位最低,表达式为:
式中,t为调度时序;w为水文集合预报发布的时刻;Zt i为t时刻第i个洪水情景下的水库水位(m);P(ωt i)为第i个洪水情景出现的概率;
第二个目标为期望最大出库流量最小,表达式为:
式中,Qt i是t时刻第i个洪水情景下的水库出库流量(m3/s);
水库实时多目标随机优化调度模型包括以下约束条件:
1)水量平衡约束:
Vt i=Vt-1 i+(It i-Qt i)·△t (8);
式中,Vt i和Vt-1 i分别为第i个洪水情景下t-1时刻的水库蓄水量;It i和Qt i分别为t时刻第i个洪水情景下的水库入库流量和出库流量(m3/s);△t为调度时间间隔(s);
2)水库水位约束:
式中,Zt 和分别为t时刻第i个洪水情景下水库水位的最低和最高限值;
3)水库出库流量约束:
式中,Qt 和分别为t时刻第i个洪水情景下水库出库流量的最低和最高限值;
4)水库泄流能力约束:
Qt i≤Qmax(Zt i) (11);
式中,Qmax(Zt i)为t时刻第i个洪水情景下的水库水位对应的水库泄流能力;
5)水库出库流量变幅约束:
|Qt i-Qt-1 i|≤△Qm (12);
式中,△Qm为允许的水库出库流量最大变幅,该约束主要用于将相邻时段的水库出库流量限定在指定范围以内;
6)初始和边界条件:
Z1 i=ZINI;ZT i→ZFL (13);
式中,ZINI分别为ZFL初始水库水位和目标期末水位。
4.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(3)包括:
(31)确定每个优化目标的上下界;先优化第一个目标而不考虑第二个目标,再优化第二个目标而不考虑第一个目标;然后,采用约束法将多目标问题转换为一系列单目标优化问题;选择第二个目标作为当前的单个优化目标,将第一个目标作为额外的约束条件加到约束集中;在上下界之间参数化摄动第一个目标值即可获得一系列非劣解集;选取的决策变量为所有洪水情景下每个时段的水库出库流量;采用非线性规划求解上述转换后的单目标优化问题;期望最大出库流量为
(32)将步骤(31)中求得的期望最大出库流量作为固定的约束值,仅优化第一个目标,从而量化水库最高水位的不确定性;期望最高水位为
通过采用上述约束法结合两步求解过程,可以获得不确定环境下的帕累托前沿,并将水库最高水位和最大出库流量描述为服从一定分布的随机变量。
5.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(4)风险评估模型包括水库漫坝风险和下游洪灾风险,将水库漫坝风险定义为最高水位超过安全水位的概率,表达式为:
式中,f(Zmax)为最高水位Zmax的概率密度函数;
将下游洪灾风险定义为最大出库流量超过安全泄量的概率,表达式为:
式中,f(Qmax)为最大出库流量Qmax的概率密度函数。
6.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(5)中利用层次分析法将决策者的主观偏好信息量化为指标权重,首先从帕累托前沿的概率分布中随机抽样生成N个指标值样本;然后将该样本和指标权重作为TOPSIS方法的输入,通过N次蒙特卡罗迭代计算贴近度系数值;最后通过统计分析推求贴近度系数的概率分布,进而确定各方案的概率排序,建立随机多属性决策模型,从帕累托前沿上确定最理想的水库调度方案;采用排序不确定度RUD定量评估决策结果的整体不确定性程度,排序不确定度定义为:各方案获得除各自确定性排序以外所有排序的概率和,计算公式如下:
式中,bi r表示确定性排序为i的方案在随机因素的作用下获得排序r的概率,表示确定性排序为i的方案在随机因素的作用下获得排序si的概率;si为第i个方案的确定性排序。
7.根据权利要求1所述的一种水库实时多目标随机优化调度和风险评估方法,其特征在于,所述步骤(6)包括:
(61)随机生成水文集合预报并构建洪水情景树,作为水库实时多目标随机优化调度模型的输入;
(62)基于面临时段w发布的水文集合预报[fw,w i,fw,w+1 i,…,fw,T i],采用约束法、两阶段求解过程和非线性规划相耦合的方法确定期望意义上的非劣解集
(63)对每个非劣解进行风险评估,计算水库漫坝风险以及下游洪灾风险;
(64)基于决策者提供的主观偏好信息,采用随机多属性决策方法从帕累托前沿上选择令决策者最为满意的期望非劣解
(65)对于所选择的期望决策序列仅实施面临时刻w的期望泄流策略后续时段的泄流策略将随着逐时段更新的水文集合预报而动态更新;
(66)随着时间推移至w+1时段,采用水量平衡公式更新水库状态,并将水文集合预报更新为[fw+1,w+1 i,fw+1,w+2 i,…,fw+1,T i]。
重复上述步骤(61)至步骤(66)直至水库调度期末,即获得每个时段的水库泄流策略和相应的风险信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810293796.2A CN108681783A (zh) | 2018-04-04 | 2018-04-04 | 一种水库实时多目标随机优化调度和风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810293796.2A CN108681783A (zh) | 2018-04-04 | 2018-04-04 | 一种水库实时多目标随机优化调度和风险评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108681783A true CN108681783A (zh) | 2018-10-19 |
Family
ID=63799461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810293796.2A Pending CN108681783A (zh) | 2018-04-04 | 2018-04-04 | 一种水库实时多目标随机优化调度和风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108681783A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109636098A (zh) * | 2018-10-31 | 2019-04-16 | 华中科技大学 | 一种基于风险熵的防洪调度风险分析方法 |
CN110033164A (zh) * | 2019-03-04 | 2019-07-19 | 华中科技大学 | 一种水库群联合防洪调度的风险评估与决策方法 |
CN110348701A (zh) * | 2019-06-21 | 2019-10-18 | 华中科技大学 | 一种水库群防洪调度风险传递规律分析方法 |
CN110969308A (zh) * | 2019-12-10 | 2020-04-07 | 浙江工业大学 | 基于动态关键任务的启发式资源受限项目调度优化方法 |
CN111461421A (zh) * | 2020-03-28 | 2020-07-28 | 华中科技大学 | 一种基于互馈关系解析的梯级水库风险评估方法及系统 |
CN111626592A (zh) * | 2020-05-22 | 2020-09-04 | 江苏建筑职业技术学院 | 一种工程进度中不确定风险的分析方法 |
CN112950033A (zh) * | 2021-03-04 | 2021-06-11 | 吴统明 | 一种基于水库调度规则合成的水库调度决策方法及系统 |
CN113191543A (zh) * | 2021-04-27 | 2021-07-30 | 西安理工大学 | 考虑水资源供需过程关联性与随机性的水库优化调度方法 |
CN113326632A (zh) * | 2021-06-19 | 2021-08-31 | 南昌工程学院 | 一种反推水库入库流量过程的优化修正方法 |
CN113433896A (zh) * | 2021-05-24 | 2021-09-24 | 深圳大学 | 一种动态生产调度方法及系统 |
CN113469528A (zh) * | 2021-06-30 | 2021-10-01 | 河海大学 | 一种考虑时空关联多维不确定性的水库群多目标防洪调度风险分析方法 |
CN113516305A (zh) * | 2021-06-29 | 2021-10-19 | 太湖流域管理局水利发展研究中心 | 一种情景-目标交互的水网区水资源调度智能决策方法及系统 |
CN113688542A (zh) * | 2021-10-26 | 2021-11-23 | 长江水利委员会长江科学院 | 智能优化水资源配置方法、装置、计算机设备及存储介质 |
CN114022024A (zh) * | 2021-11-17 | 2022-02-08 | 中山大学 | 一种基于信息论指标的水库调度模块评价方法及系统 |
CN117910711A (zh) * | 2024-03-20 | 2024-04-19 | 长江水利委员会长江科学院 | 均衡风险与效益的汛期水位智能动态控制模型的构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248015A (zh) * | 2017-07-07 | 2017-10-13 | 武汉大学 | 基于集合预报的水库实时蓄水调度方法 |
-
2018
- 2018-04-04 CN CN201810293796.2A patent/CN108681783A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107248015A (zh) * | 2017-07-07 | 2017-10-13 | 武汉大学 | 基于集合预报的水库实时蓄水调度方法 |
Non-Patent Citations (5)
Title |
---|
FEILIN ZHU 等: "《Real-Time Optimal Flood Control Decision Making and Risk Propagation Under Multiple Uncertainties》", 《WATER RESOURCES RESEARCH》 * |
孙兆峰: "《基于风险理论的水库汛限水位控制研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
欧阳硕等: "《金沙江下游梯级与三峡梯级多目标联合防洪优化调度研究》", 《水力发电学报》 * |
赵铜铁钢: "《考虑水文预报不确定性的水库优化调度研究》", 《中国博士学位论文全文数据库》 * |
陈璐等: "《水文预报不确定性对水库防洪调度的影响分析》", 《水利学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109636098A (zh) * | 2018-10-31 | 2019-04-16 | 华中科技大学 | 一种基于风险熵的防洪调度风险分析方法 |
CN110033164A (zh) * | 2019-03-04 | 2019-07-19 | 华中科技大学 | 一种水库群联合防洪调度的风险评估与决策方法 |
CN110348701A (zh) * | 2019-06-21 | 2019-10-18 | 华中科技大学 | 一种水库群防洪调度风险传递规律分析方法 |
CN110348701B (zh) * | 2019-06-21 | 2021-10-08 | 华中科技大学 | 一种水库群防洪调度风险传递规律分析方法 |
CN110969308A (zh) * | 2019-12-10 | 2020-04-07 | 浙江工业大学 | 基于动态关键任务的启发式资源受限项目调度优化方法 |
CN111461421B (zh) * | 2020-03-28 | 2021-02-09 | 华中科技大学 | 一种基于互馈关系解析的梯级水库风险评估方法及系统 |
CN111461421A (zh) * | 2020-03-28 | 2020-07-28 | 华中科技大学 | 一种基于互馈关系解析的梯级水库风险评估方法及系统 |
CN111626592A (zh) * | 2020-05-22 | 2020-09-04 | 江苏建筑职业技术学院 | 一种工程进度中不确定风险的分析方法 |
CN111626592B (zh) * | 2020-05-22 | 2023-09-05 | 江苏建筑职业技术学院 | 一种工程进度中不确定风险的分析方法 |
CN112950033A (zh) * | 2021-03-04 | 2021-06-11 | 吴统明 | 一种基于水库调度规则合成的水库调度决策方法及系统 |
CN112950033B (zh) * | 2021-03-04 | 2024-03-29 | 吴统明 | 一种基于水库调度规则合成的水库调度决策方法及系统 |
CN113191543B (zh) * | 2021-04-27 | 2022-07-29 | 西安理工大学 | 考虑水资源供需过程关联性与随机性的水库优化调度方法 |
CN113191543A (zh) * | 2021-04-27 | 2021-07-30 | 西安理工大学 | 考虑水资源供需过程关联性与随机性的水库优化调度方法 |
CN113433896A (zh) * | 2021-05-24 | 2021-09-24 | 深圳大学 | 一种动态生产调度方法及系统 |
CN113326632A (zh) * | 2021-06-19 | 2021-08-31 | 南昌工程学院 | 一种反推水库入库流量过程的优化修正方法 |
CN113516305B (zh) * | 2021-06-29 | 2022-05-03 | 太湖流域管理局水利发展研究中心 | 一种情景-目标交互的水网区水资源调度智能决策方法及系统 |
CN113516305A (zh) * | 2021-06-29 | 2021-10-19 | 太湖流域管理局水利发展研究中心 | 一种情景-目标交互的水网区水资源调度智能决策方法及系统 |
CN113469528A (zh) * | 2021-06-30 | 2021-10-01 | 河海大学 | 一种考虑时空关联多维不确定性的水库群多目标防洪调度风险分析方法 |
CN113469528B (zh) * | 2021-06-30 | 2023-09-12 | 河海大学 | 一种考虑时空关联多维不确定性的水库群多目标防洪调度风险分析方法 |
CN113688542A (zh) * | 2021-10-26 | 2021-11-23 | 长江水利委员会长江科学院 | 智能优化水资源配置方法、装置、计算机设备及存储介质 |
CN114022024A (zh) * | 2021-11-17 | 2022-02-08 | 中山大学 | 一种基于信息论指标的水库调度模块评价方法及系统 |
CN117910711A (zh) * | 2024-03-20 | 2024-04-19 | 长江水利委员会长江科学院 | 均衡风险与效益的汛期水位智能动态控制模型的构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108681783A (zh) | 一种水库实时多目标随机优化调度和风险评估方法 | |
Hung et al. | An artificial neural network model for rainfall forecasting in Bangkok, Thailand | |
Coulibaly et al. | Downscaling precipitation and temperature with temporal neural networks | |
Luk et al. | A study of optimal model lag and spatial inputs to artificial neural network for rainfall forecasting | |
JP2019194424A (ja) | 河川流域に関する水文解析および管理のためのプロセスおよびシステム | |
Meng et al. | Flood control operation coupled with risk assessment for cascade reservoirs | |
Lai et al. | Flood risk assessment and regionalization from past and future perspectives at basin scale | |
Chen et al. | Groundwater level prediction using SOM-RBFN multisite model | |
Adeyemo et al. | River flow forecasting using an improved artificial neural network | |
CN113050430B (zh) | 一种基于鲁棒强化学习的排水系统控制方法 | |
Donyaii et al. | Water reservoir multiobjective optimal operation using Grey Wolf optimizer | |
Ferreira et al. | Optimal and adaptive operation of a hydropower system with unit commitment and water quality constraints | |
Liu et al. | A risk assessment method based on RBF artificial neural network-cloud model for urban water hazard | |
CN114492233B (zh) | 一种基于webGIS平台及考虑综合利用要求的流域水模拟方法 | |
Cameron et al. | An evaluation of a traditional and a neural net modelling approach to flood forecasting for an upland catchment | |
CN116050628A (zh) | 一种基于数字孪生的考虑预报不确定性的水库防洪优化调度方法 | |
Huang et al. | Identification and prediction of urban traffic congestion via cyber-physical link optimization | |
Dariane et al. | Long term streamflow forecasting using a hybrid entropy model | |
Evsukoff et al. | A multi-model approach for long-term runoff modeling using rainfall forecasts | |
CN114399235B (zh) | 一种基于雨情数据判定灾害风险等级的方法及系统 | |
Wu et al. | Assessment of waterlogging risk in the deep foundation pit projects based on projection pursuit model | |
Nozari et al. | Simulation and optimization of control system operation and surface water allocation based on system dynamics modeling | |
Goyal et al. | Flash flood risk management modeling in indian cities using IoT based reinforcement learning | |
Jain et al. | State-of-the-art review: Operation of multi-purpose reservoirs during flood season | |
Zhijun | RBF neural networks optimization algorithm and application on tax forecasting |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |