CN113627673B - 一种调水工程多水力设施联合优化调控方案制定方法 - Google Patents
一种调水工程多水力设施联合优化调控方案制定方法 Download PDFInfo
- Publication number
- CN113627673B CN113627673B CN202110920075.1A CN202110920075A CN113627673B CN 113627673 B CN113627673 B CN 113627673B CN 202110920075 A CN202110920075 A CN 202110920075A CN 113627673 B CN113627673 B CN 113627673B
- Authority
- CN
- China
- Prior art keywords
- model
- regulation
- hydraulic
- joint optimization
- scheme set
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 135
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 82
- 230000001052 transient effect Effects 0.000 claims abstract description 74
- 230000008569 process Effects 0.000 claims abstract description 53
- 238000012795 verification Methods 0.000 claims abstract description 20
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000010168 coupling process Methods 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 238000012544 monitoring process Methods 0.000 claims abstract description 6
- 238000006467 substitution reaction Methods 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000012546 transfer Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 230000001105 regulatory effect Effects 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 239000004973 liquid crystal related substance Substances 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
- 238000009834 vaporization 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
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- 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/06—Energy or water supply
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
Abstract
一种调水工程多水力设施联合优化调控方案制定方法,获取调水工程沿线水情及运行状态监测数据,率定并验证得到符合调水工程实际运行的水力瞬变模型;建立多水力设施联合优化调控模型,与所述水力瞬变模型进行耦合,得到父代调控方案集;并构建所述水力瞬变模型的替代模型,通过所述替代模型和所述多水力设施联合优化调控模型的耦合,不断对所述父代调控方案集进行迭代调整,得到多水力设施联合优化调控的最优调控方案集。本发明可以有效提高水力瞬变模型模拟的准确性,提高多水力设施联合优化调控模型的效率和实用性,以及调控方案的切实性和可操作性,从而及时指导实际工程的调控运行,减少调控过程风险,避免反复调控,提高调控运行品质。
Description
技术领域
本发明涉及调水工程领域,尤其涉及一种调水工程多水力设施联合优化调控方案制定方法。
背景技术
为了均衡水资源空间配置,缓解缺水地区水资源短缺,促进社会、经济以及生态等方面的协调发展,使得调水工程的建设成为必然。但由于地形、地质、城市建设、调水规模等因素的限制,调水方式、结构、水力设施日趋复杂,水力响应也越来越复杂,对于多水力设施的协调联合调控需求越来越高。而多水力设施间调控相互干扰、相互影响,若联合调控不当,往往需要反复调控,机械损耗大,人力消耗大,决策效率低,还可能造成管道系统运行安全风险的叠加,增加最大瞬时压力超过管道的承受能力,进而引起爆管、阀门损坏等事故的风险;增加最小瞬时压力降至对应环境温度下的水体汽化压力,使得管道局部呈现真空状态,进而造成管道压瘪、坍塌,或引发水柱分离和弥合现象的风险;增加沿线调节池出现漏空或漫溢、隧洞局部断面水位出现急剧下降或出现明满交替现象、泵站前池水位出现低于最低运行水位或漫溢,水泵出现空蚀空化等的风险,运行调控难度大。
因此,综合考虑各安全风险,寻求多水力设施联合优化调控方案制定方法,对有效指导调水工程的调度运行具有重要的应用价值和科学意义。而现有技术中利用优化算法结合水力瞬变模型制定多水力设施联合优化调控方案,往往需要进行成百上千次的瞬变计算,耗时长,难以满足决策的需求。
发明内容
本发明的目的在于提供一种调水工程多水力设施联合优化调控方案制定方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明采用的技术方案如下:
一种调水工程多水力设施联合优化调控方案制定方法,包括以下步骤:
S1、整理调水工程沿线水情及运行状态监测数据,进行率定并验证得到符合调水工程实际运行的水力瞬变模型;
S2、根据调水工程瞬变过程中安全风险的需求,建立多水力设施联合优化调控模型,并与所述水力瞬变模型进行耦合;
S3、当迭代次数为1时,对耦合所述水力瞬变模型的所述多水力设施联合优化调控模型进行求解,得到原始的决策变量对应的目标函数值,形成所述多水力设施联合优化调控的父代调控方案集;
S4、根据所述父代调控方案集中各方案的决策变量和对应的目标函数值,建立所述水力瞬变模型的替代模型,进而将所述替代模型和所述多水力设施联合优化模型进行耦合;
S5、当迭代次数大于1时,将所述父代调控方案集中的所述决策变量通过优化算法生成第二代决策变量,通过耦合所述替代模型的所述水力设施联合优化调控模型进行求解,得到所述第二代决策变量对应的目标函数值,形成所述多水力设施联合优化调控的子代调控方案集;
S6、将所述父代调控方案集与所述子代调控方案集进行结合,并按照非支配等级层筛选得到新父代调控方案集;
S7、将所述父代调控方案集与所述子代调控方案集进行结合并选择出非支配方案集,将所述非支配方案集代入所述水力瞬变模型进行校验,得到校验非支配方案集,将得到的所述校验非支配方案集并入步骤S3形成的所述父代调控方案集作为训练方案集,对所述替代模型进行训练,实现所述替代模型的第一次迭代;
S8、利用所述新父代调控方案集代替所述父代调控方案集,并利用耦合迭代后替代模型的所述水力设施联合优化调控模型求解,重复步骤S5—S7;直至迭代次数达到设置的最大迭代次数,得到所述多水力设施联合优化调控模型的最优调控方案集。
优选的,所述替代模型的输入值为所述训练方案集中各方案的决策变量,输出值为所述训练方案集中各方案的目标函数;所述替代模型的高斯曲面为:
其中,X,y分别表示输入和输出矩阵;X.表示预测输入矩阵,即用来检测替代模型的样本;f.表示预测输出向量。
优选的,所述多水力设施联合优化调控模型与所述水力瞬变模型之间的耦合方式为:将所述多水力设施联合优化调控模型生成的决策变量作为所述水力瞬变模型的输入,根据所述水力瞬变模型生成计算所述多水力设施联合优化调控模型的目标函数所需变量。
优选的,耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤:
S31、根据所述多水力设施联合优化调控模型建立需要进行调控的水力设施决策变量种群,并对所述水力设施决策变量种群中的水力设施设置约束条件,通过所述水力设施的约束条件对所述多水力设施联合优化调控模型的决策变量进行处理;
S32、驱动所述水力瞬变模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S33、根据所述目标函数的所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值。
优选的,耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程并行处理,作为所述多水力设施联合优化调控模型的第一并行域。
优选的,耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤:
S51、根据水力设施的约束条件,对所述多水力设施联合优化调控模型的第二代决策变量进行处理;
S52、驱动所述替代模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S53、根据步骤S52得到的所述第二代决策变量对应的目标函数的所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值,形成非支配方案集。
优选的,所述第二代决策变量的生成过程与耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程并行处理,作为所述多水力设施联合优化调控模型的第二并行域。
优选的,所述父代调控方案集、所述子代调控方案集、所述非支配方案集和所述校验非支配方案集均为调控方案集,所述调控方案集包括所述多水力设施联合优化调控模型的决策变量及其对应的目标函数值。
优选的,所述父代调控方案集首先通过NSGA-II优化算法进行非支配等级排序,再将所述父代调控方案集中的决策变量通过优化算法生成第二代决策变量。
优选的,驱动所述水力瞬变模型对所述非支配方案集进行的校验过程并行处理,作为所述多水力设施联合优化调控模型的第三并行域。
本发明的有益效果是:本发明公开了一种调水工程多水力设施联合优化调控方案制定方法,将水力瞬变模型和多水力设施联合优化调控模型进行耦合,得到满足约束条件的调控方案,并通过建立所述瞬变模型的替代模型对得到的调控方案不断进行迭代更新,得到多水力设施联合优化调控的最优调控方案。本发明可以有效提高水力瞬变模型模拟的准确性,提高多水力设施联合优化调控的效率和实用性,提高多水力设施联合调控方案的切实性、可操作性,从而可以及时的指导实际工程的调控运行,减少调控过程风险、避免反复调控、提高调控运行品质。
附图说明
图1是多水力设施联合优化调控模型方案制定流程图;
图2是隧洞特征线网格示意图;
图3是多水力设施联合优化调控模型与水力瞬变模型及替代模型耦合关系图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
一种调水工程多水力设施联合优化调控方案制定方法,涉及调水工程多水力设施联合优化调控技术领域,如图1所示,包括以下步骤:
S1、整理调水工程沿线水情及运行状态监测数据,进行率定并验证得到符合调水工程实际运行的水力瞬变模型;
所述调水工程沿线水情及运行状态监测数据包括管道流量、管道压力,隧洞流量、隧洞水位,调节池水位,泵站水泵前后压力、水泵过流量、水泵调频等运行状态数据,沿线阀门前后压力、开度、过流量等运行状态数据;
S11、结合调水工程系统拓扑图,利用特征线法建立调水工程管道内部节点数学模型、隧洞内部节点数学模型以及管道连接边界、隧洞连接边界、管道与隧洞连接边界、水泵边界、空气阀边界、阀门边界、水池边界、定流量、定水位边界等设备边界模型;
管道内部节点的特征线方程为:
其中, i表示管道上的结点号;j表示时层号;Q表示管道流量,代替控制方程中的平均流速;A表示管道横截面积;Δt表示时间步长;Δx表示管道分段长,即空间步长;g表示重力加速度,f表示水力摩阻系数,H表示液体压头;
如图2隧洞特征线网格所示,隧洞内部节点的特征线方程为:
其中,分别为f在P与L两点、P与R两点上的平均值;/> 分别为Bc-在P与L两点,Bc+在P与R两点的平均值;对于R、L点上断面的水深、流量等信息可通过同一时层已知的网格结点A、B、C三点的值的二次插值得到;
该所述隧洞内部节点的特征线满足库朗稳定条件:
各设备边界模型可通过结合边界自身的物理特性和管道、隧洞的特征线方程得到;
S12、分别选取两组包含工况变化的持续120分钟的步长为0.01秒的水情、运行状态监测数据,以其中一组数据率定管道糙率、隧洞糙率、阀门以及水泵特性的参数,另外一组作为验证,控制验证误差在10%以内,以得到符合调水工程实际运行的水力瞬变模型。
S2、根据调水工程瞬变过程各安全风险的需求,建立多水力设施联合优化调控模型,并利用调用外部程序的方法,与所述水力瞬变模型进行耦合;
S21、调水工程多水力设施联合调控主要包括对水泵、泵后阀以及沿线多个阀门的联合调控,其中决策变量主要包括每一个管道的m个水泵中相邻水泵之间关机的间隔时长tde,d表示管道的编码,e表示同一个管道中相邻水泵之间间隔的编号,e=1,2,…,m-1;泵后阀快关过程或慢关过程的关闭时长t1;泵后阀快关过程或慢关过程相对应的开度θ;泵后阀的开度从全开到全关的总关闭时长t2;以及沿线n个阀门的由于水泵依次停机的调控开度βlk,l表示沿线阀门的编号,k表示水泵的编号,l=1,2,…,n,k=1,2,…,m;
S22、根据调水工程瞬变过程各安全的需求,确立目标函数,包括:调水系统在瞬变过程中沿线的最大压力最小化min(Hmax)、最小压力最大化max(Hmin)、机组最大倒转速最小化max(nmin)、控制建筑物最低水位最大化max(zqmin)和控制建筑物最高水位最小化min(zqmax);
其中,Hmax,Hmin分别表示调水工程系统在瞬变过程中沿线最大压力和最小压力,nmin表示泵站机组事故停电时水泵的最大倒转速,倒转速为负值,q表示调节池、隧洞进出口等控制建筑物的编号,zqmin、zqmax表示控制建筑物瞬变过程中的最低、最高水位;
S23、根据水泵、泵后阀、沿线阀门等各水力设施的调控要求以及系统安全风险的要求,设置约束条件:
根据泵后阀的制造要求、泵站事故停机、阀门拒动工况下0流量出现的时刻,得到泵后阀关闭规律中决策变量泵后阀快关过程或慢关过程的开度、快关过程或慢关过程的时长以及总关闭时间的约束条件;所述多水力设施联合优化调控模型中对于其中的多个水力设施中泵后阀的约束条件为:
其中,tc_min表示多水力设施联合调控系统中相应的泵后阀可操控的最短时长,tc_max表示多水力设施联合调控系统中相应的泵后阀可操纵的最长时长;
且当所述泵后阀的关闭顺序为先快关后慢关时,所述泵后阀的关闭规律为:
当所述泵后阀的关闭顺序为先慢关后快关时,所述泵后阀的关闭规律为:
对于所述多水力设施联合优化调控模型中水泵的约束条件为:nmin>-1.2nr;
对于所述多水力设施联合优化调控模型中管道的约束条件为:
Hmax≤Hs_max=1.3Hw_max
Hmin>-7.5m
其中,Hw_max表示沿线稳态时的最大压力值,Hs_max表示瞬变过程允许出现的最大压力,nr表示水泵的额定转速,Hmax、Hmin分别表示系统在瞬变过程中沿线最大压力和最小压力,nmin表示泵站机组事故停电时水泵的最大倒转速;
对于所述多水力设施联合优化调控模型中控制断面运行水位的约束条件为:
其中,Zqs max和Zqs min分别表示控制建筑物设计规定的最高及最低运行水位,Zqmax和Zqmin分别表示瞬变过程中建筑物出现的最高及最低水位;
S24、根据所述多水力设施联合优化调控模型的决策变量、目标函数和约束条件,建立多水力设施联合优化调控模型;将所述多水力设施联合优化调控模型和所述水力瞬变模型进行耦合,如图3所示具体的耦合方式为:将所述多水力设施联合优化调控模型生成的决策变量传递给所述水力瞬变模型,作为所述水力瞬变模型的输入,利用所述水力瞬变模型进行瞬变过程计算,并将计算管道最大、最小压力,隧洞进出口、调节池最高、最低水位,机组最大倒转速等多水力设施联合优化调控模型的目标函数所需的变量传递至多水力设施联合优化调控模型,实现模型间的耦合。
所述多水力设施联合优化调控模型的求解流程如图1所示:依次进行种群初始化、种群约束条件并行处理、并行驱动水力瞬变模型进行瞬变过程计算、并行目标函数值的计算、非支配等级排序、拥挤距离计算、选择计算、交叉计算、变异过程并行计算、并行驱动所述替代模型进行瞬变过程计算、并行目标函数值的计算、生成子代、父代子代融合、新父代的生成;针对上述求解过程中的并行性进行分析,所述多水力设施联合优化调控模型在求解过程中的并行域包括种群约束条件处理、驱动水力瞬变模型进行瞬变过程计算、目标函数值的计算、驱动替代模型进行瞬变过程计算以及NSGA-II优化算法生成子代决策变量中变异计算过程。
S3、当迭代次数为1时,利用并行算法求解耦合所述水力瞬变模型的多水力设施联合优化调控模型,得到原始的决策变量对应的目标函数值,形成所述多水力设施联合优化调控的父代调控方案集;
上述耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程,具体包括以下步骤:
S31、根据所述多水力设施联合优化调控模型建立需要进行调控的水力设施决策变量种群,对所述水力设施决策变量种群进行初始化;并对所述水力设施决策变量种群中的水力设施设置约束条件,通过所述水力设施的约束条件对所述多水力设施联合优化调控模型的决策变量进行处理;
S32、驱动所述水力瞬变模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S33、根据所述目标函数所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值;
其中步骤S31—步骤S33作为耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程的并行域,作为所述多水力设施联合优化调控模型的第一并行域。
S4、根据所述父代调控方案集中各方案的决策变量和对应的目标函数值,作为水力瞬变模型的输入值和输出值,利用matlab采用高斯曲面,建立所述水力瞬变模型的替代模型,进而采用所述水力瞬变模型与所述多水力设施联合优化模型之间的耦合方式,将所述替代模型和所述多水力设施联合优化模型进行耦合;
所述替代模型的输入值为训练方案集中各方案的决策变量,输出值为训练方案集中各方案的目标函数;所述替代模型的高斯曲面为:
其中,X,y表示输入和输出矩阵;X.表示预测输入矩阵,即用来检测替代模型的样本;f.表示预测输出向量。
S5、针对所述父代调控方案集通过NSGA-II优化算法进行非支配等级排序和拥挤距离计算,当迭代此处大于1时,再将所述父代调控方案集中的决策变量通过选择计算、交叉计算、变异过程等优化算法得到第二代决策变量;通过耦合所述替代模型的所述水力设施联合优化调控模型进行求解,得到所述第二代决策变量对应的目标函数值,形成所述多水力设施联合优化调控的子代调控方案集;具体包括一下步骤:
S51、根据水力设施的约束条件,对所述多水力设施联合优化调控模型的第二代决策变量进行处理;
S52、驱动所述替代模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S53、根据步骤S52得到的所述第二代决策变量对应的目标函数的所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值,形成非支配方案集;
其中,所述第二代决策变量的生成过程与步骤S51—步骤S53作为耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程的并行域,作为所述多水力设施联合优化调控模型的第二并行域。
S6、将所述父代调控方案集与所述子代调控方案集进行结合,并通过NSGA-II优化算法进行非支配等级排序和拥挤距离计算,按照非支配等级层筛选得到与所述父代调控方案集中元素数量相同的新父代调控方案集。
S7、所述父代调控方案集与所述子代调控方案集进行结合,并通过NSGA-II优化算法进行非支配等级排序和拥挤距离计算,按照非支配等级层筛选出非支配方案集,将所述非支配方案集代入所述水力瞬变模型进行校验,得到校验非支配方案集,将得到的所述校验非支配方案集并入步骤S3形成的所述父代调控方案集作为所述训练方案集,对所述替代模型进行训练,实现所述替代模型的第一次迭代;
其中,驱动所述水力瞬变模型对所述非支配方案集进行的校验过程并行处理,作为所述多水力设施联合优化调控模型的第三并行域。
S8、利用所述新父代调控方案集代替所述父代调控方案集,并利用耦合迭代后替代模型的水力设施联合优化调控模型求解,重复步骤S5—S7;直至迭代次达到设置的最大迭代次数,得到的所述多水力设施联合优化调控模型的最优调控方案集,通过所述调控方案集对所述调水工程多水力设施的调控进行指导。
所述父代调控方案集、所述子代调控方案集、所述非支配方案集和所述校验非支配方案集均为调控方案集,所述调控方案集包括所述多水力设施联合优化调控模型的决策变量及其对应的目标函数值。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明公开了一种调水工程多水力设施联合优化调控方案制定方法,将水力瞬变模型和多水力设施联合优化调控模型进行耦合,得到满足约束条件的调控方案,并通过建立所述瞬变模型的替代模型对得到的调控方案不断进行迭代更新,得到多水力设施联合优化调控的最优调控方案。本发明可以有效提高水力瞬变模型模拟的准确性,提高多水力设施联合优化调控的效率和实用性,提高多水力设施联合调控方案的切实性、可操作性,从而可以及时的指导实际工程的调控运行,减少调控过程风险、避免反复调控、提高调控运行品质。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (10)
1.一种调水工程多水力设施联合优化调控方案制定方法,其特征在于,包括以下步骤:
S1、整理调水工程沿线水情及运行状态监测数据,进行率定并验证得到符合调水工程实际运行的水力瞬变模型;
S2、根据调水工程瞬变过程中安全风险的需求,建立多水力设施联合优化调控模型,并与所述水力瞬变模型进行耦合;
S3、当迭代次数为1时,对耦合所述水力瞬变模型的所述多水力设施联合优化调控模型进行求解,得到原始的决策变量对应的目标函数值,形成所述多水力设施联合优化调控的父代调控方案集;
S4、根据所述父代调控方案集中各方案的决策变量和对应的目标函数值,建立所述水力瞬变模型的替代模型,进而将所述替代模型和所述多水力设施联合优化模型进行耦合;
S5、当迭代次数大于1时,将所述父代调控方案集中的所述决策变量通过优化算法生成第二代决策变量;通过耦合所述替代模型的所述水力设施联合优化调控模型进行求解,得到所述第二代决策变量对应的目标函数值,形成所述多水力设施联合优化调控的子代调控方案集;
S6、将所述父代调控方案集与所述子代调控方案集进行结合,并按照非支配等级层筛选得到新父代调控方案集;
S7、将所述父代调控方案集与所述子代调控方案集进行结合并选择出非支配方案集,将所述非支配方案集代入所述水力瞬变模型进行校验,得到校验非支配方案集,将得到的所述校验非支配方案集并入步骤S3形成的所述父代调控方案集作为训练方案集,对所述替代模型进行训练,实现所述替代模型的第一次迭代;
S8、利用所述新父代调控方案集代替所述父代调控方案集,并利用耦合迭代后替代模型的所述水力设施联合优化调控模型求解,重复步骤S5—S7;直至迭代次数达到设置的最大迭代次数,得到所述多水力设施联合优化调控模型的最优调控方案集。
2.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,所述替代模型的输入值为所述训练方案集中各方案的决策变量,输出值为所述训练方案集中各方案的目标函数;所述替代模型的高斯曲面为:
其中,X,y分别表示输入和输出矩阵;X.表示预测输入矩阵,即用来检测替代模型的样本;f.表示预测输出向量。
3.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,所述多水力设施联合优化调控模型与所述水力瞬变模型之间的耦合方式为:将所述多水力设施联合优化调控模型生成的决策变量作为所述水力瞬变模型的输入,根据所述水力瞬变模型生成计算所述多水力设施联合优化调控模型的目标函数所需变量。
4.根据权利要求3所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤:
S31、根据所述多水力设施联合优化调控模型建立需要进行调控的水力设施决策变量种群,并对所述水力设施决策变量种群中的水力设施设置约束条件,通过所述水力设施的约束条件对所述多水力设施联合优化调控模型的决策变量进行处理;
S32、驱动所述水力瞬变模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S33、根据所述目标函数的所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值。
5.根据权利要求4所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程并行处理,作为所述多水力设施联合优化调控模型的第一并行域。
6.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤:
S51、根据水力设施的约束条件,对所述多水力设施联合优化调控模型的第二代决策变量进行处理;
S52、驱动所述替代模型进行瞬变过程中的计算,得到计算所述多水力设施联合优化调控模型的目标函数所需变量;
S53、根据步骤S52得到的所述第二代决策变量对应的目标函数的所需变量,通过所述多水力设施联合优化模型的计算得到对应的目标函数值,形成非支配方案集。
7.根据权利要求6所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,所述第二代决策变量的生成过程与耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程并行处理,作为所述多水力设施联合优化调控模型的第二并行域。
8.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,所述父代调控方案集、所述子代调控方案集、所述非支配方案集和所述校验非支配方案集均为调控方案集,所述调控方案集包括所述多水力设施联合优化调控模型的决策变量及其对应的目标函数值。
9.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,所述父代调控方案集首先通过NSGA-II优化算法进行非支配等级排序,再将所述父代调控方案集中的决策变量通过优化算法生成第二代决策变量。
10.根据权利要求1所述的调水工程多水力设施联合优化调控方案制定方法,其特征在于,驱动所述水力瞬变模型对所述非支配方案集进行的校验过程并行处理,作为所述多水力设施联合优化调控模型的第三并行域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110920075.1A CN113627673B (zh) | 2021-08-11 | 2021-08-11 | 一种调水工程多水力设施联合优化调控方案制定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110920075.1A CN113627673B (zh) | 2021-08-11 | 2021-08-11 | 一种调水工程多水力设施联合优化调控方案制定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113627673A CN113627673A (zh) | 2021-11-09 |
CN113627673B true CN113627673B (zh) | 2023-08-04 |
Family
ID=78384583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110920075.1A Active CN113627673B (zh) | 2021-08-11 | 2021-08-11 | 一种调水工程多水力设施联合优化调控方案制定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113627673B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006353068A (ja) * | 2005-06-20 | 2006-12-28 | Chugoku Electric Power Co Inc:The | 流込式水力発電設備における発電出力制御装置、水位調整装置および、発電出力制御方法 |
JP2015074971A (ja) * | 2013-10-07 | 2015-04-20 | 駿 横山 | 砂防ダム水力発電システム |
CN107355684A (zh) * | 2017-07-19 | 2017-11-17 | 中国水利水电科学研究院 | 一种管网故障水力监测实验系统及其实现故障辨识的方法 |
CN111783369A (zh) * | 2020-07-22 | 2020-10-16 | 中国水利水电科学研究院 | 一种多闸群明渠调水工程的短期多目标优化调度方法 |
CN111985166A (zh) * | 2020-07-31 | 2020-11-24 | 河海大学 | 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质 |
-
2021
- 2021-08-11 CN CN202110920075.1A patent/CN113627673B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006353068A (ja) * | 2005-06-20 | 2006-12-28 | Chugoku Electric Power Co Inc:The | 流込式水力発電設備における発電出力制御装置、水位調整装置および、発電出力制御方法 |
JP2015074971A (ja) * | 2013-10-07 | 2015-04-20 | 駿 横山 | 砂防ダム水力発電システム |
CN107355684A (zh) * | 2017-07-19 | 2017-11-17 | 中国水利水电科学研究院 | 一种管网故障水力监测实验系统及其实现故障辨识的方法 |
CN111783369A (zh) * | 2020-07-22 | 2020-10-16 | 中国水利水电科学研究院 | 一种多闸群明渠调水工程的短期多目标优化调度方法 |
CN111985166A (zh) * | 2020-07-31 | 2020-11-24 | 河海大学 | 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113627673A (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2024109103A1 (zh) | 一种闸门联动控制方法、装置及并联供水发电系统 | |
CN114777192B (zh) | 基于数据关联和深度学习的二级网供热自主优化调控方法 | |
He et al. | A methodology for the optimal design of gathering pipeline system in old oilfield during its phased development process | |
CN115456422B (zh) | 一种基于计算水动力学的灌区配水计划动态预演修正方法 | |
CN112113146B (zh) | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 | |
Jung et al. | Fluid transients and pipeline optimization using GA and PSO: The diameter connection | |
Zargar et al. | Optimization of gated spillways operation for flood risk management in multi-reservoir systems | |
Kong et al. | Predictive control for the operation of cascade pumping stations in water supply canal systems considering energy consumption and costs | |
CN113361214A (zh) | 一种基于水位流量数据的明渠控制模型参数辨识方法 | |
CN113627673B (zh) | 一种调水工程多水力设施联合优化调控方案制定方法 | |
Liu et al. | An adaptive predictive control algorithm for comprehensive dendritic canal systems | |
US20180355589A1 (en) | Management of liquid conduit systems | |
Abbasi et al. | Ant-colony-based simulation–optimization modeling for the design of a forced water pipeline system considering the effects of dynamic pressures | |
CN116717839A (zh) | 供热控制方法、控制装置和供热系统 | |
CN117348400A (zh) | 一种水网多通道并联联动控制系统 | |
Bao-Feng et al. | Optimization of tree pipe networks layout and size, using particle swarm optimization | |
CN116956744A (zh) | 基于改进粒子群算法的多回路沟槽电缆稳态温升预测方法 | |
Carpitella et al. | Multi-objective and multi-criteria analysis for optimal pump scheduling in water systems | |
Colmenares et al. | MLD Systems: Modeling and Control—Experience With a Pilot Process | |
Haddad et al. | Dynamic penalty function as a strategy in solving water resources combinatorial optimization problems with honey-bee mating optimization (HBMO) algorithm | |
Arsene et al. | A study on modeling and simulation of water distribution systems based on loop corrective flows and containing controlling hydraulics elements | |
Tolson et al. | Battle of background leakage assessment for water networks (BBLAWN): An incremental savings approach | |
Jeong et al. | Implementation of simplified sequential stochastic model predictive control for operation of hydropower system under uncertainty | |
Sarbu et al. | Hydraulic simulation and optimisation of water transmission and distribution systems | |
Gangwani et al. | Optimal pipe-sizing design of water distribution networks using Modified Rao-II Algorithm |
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 |