CN113627673A - 一种调水工程多水力设施联合优化调控方案制定方法 - Google Patents

一种调水工程多水力设施联合优化调控方案制定方法 Download PDF

Info

Publication number
CN113627673A
CN113627673A CN202110920075.1A CN202110920075A CN113627673A CN 113627673 A CN113627673 A CN 113627673A CN 202110920075 A CN202110920075 A CN 202110920075A CN 113627673 A CN113627673 A CN 113627673A
Authority
CN
China
Prior art keywords
hydraulic
regulation
model
facility
control
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
CN202110920075.1A
Other languages
English (en)
Other versions
CN113627673B (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.)
Taiyuan University of Technology
Original Assignee
Taiyuan University of Technology
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 Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202110920075.1A priority Critical patent/CN113627673B/zh
Publication of CN113627673A publication Critical patent/CN113627673A/zh
Application granted granted Critical
Publication of CN113627673B publication Critical patent/CN113627673B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Feedback Control In General (AREA)

Abstract

一种调水工程多水力设施联合优化调控方案制定方法，获取调水工程沿线水情及运行状态监测数据，率定并验证得到符合调水工程实际运行的水力瞬变模型；建立多水力设施联合优化调控模型，与所述水力瞬变模型进行耦合，得到父代调控方案集；并构建所述水力瞬变模型的替代模型，通过所述替代模型和所述多水力设施联合优化调控模型的耦合，不断对所述父代调控方案集进行迭代调整，得到多水力设施联合优化调控的最优调控方案集。本发明可以有效提高水力瞬变模型模拟的准确性，提高多水力设施联合优化调控模型的效率和实用性，以及调控方案的切实性和可操作性，从而及时指导实际工程的调控运行，减少调控过程风险，避免反复调控，提高调控运行品质。

Description

一种调水工程多水力设施联合优化调控方案制定方法
技术领域
本发明涉及调水工程领域，尤其涉及一种调水工程多水力设施联合优化调控方案制定方法。
背景技术
为了均衡水资源空间配置，缓解缺水地区水资源短缺，促进社会、经济以及生态等方面的协调发展，使得调水工程的建设成为必然。但由于地形、地质、城市建设、调水规模等因素的限制，调水方式、结构、水力设施日趋复杂，水力响应也越来越复杂，对于多水力设施的协调联合调控需求越来越高。而多水力设施间调控相互干扰、相互影响，若联合调控不当，往往需要反复调控，机械损耗大，人力消耗大，决策效率低，还可能造成管道系统运行安全风险的叠加，增加最大瞬时压力超过管道的承受能力，进而引起爆管、阀门损坏等事故的风险；增加最小瞬时压力降至对应环境温度下的水体汽化压力，使得管道局部呈现真空状态，进而造成管道压瘪、坍塌，或引发水柱分离和弥合现象的风险；增加沿线调节池出现漏空或漫溢、隧洞局部断面水位出现急剧下降或出现明满交替现象、泵站前池水位出现低于最低运行水位或漫溢，水泵出现空蚀空化等的风险，运行调控难度大。
因此，综合考虑各安全风险，寻求多水力设施联合优化调控方案制定方法，对有效指导调水工程的调度运行具有重要的应用价值和科学意义。而现有技术中利用优化算法结合水力瞬变模型制定多水力设施联合优化调控方案，往往需要进行成百上千次的瞬变计算，耗时长，难以满足决策的需求。
发明内容
本发明的目的在于提供一种调水工程多水力设施联合优化调控方案制定方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明采用的技术方案如下：
一种调水工程多水力设施联合优化调控方案制定方法，包括以下步骤：
S1、整理调水工程沿线水情及运行状态监测数据，进行率定并验证得到符合调水工程实际运行的水力瞬变模型；
S2、根据调水工程瞬变过程中安全风险的需求，建立多水力设施联合优化调控模型，并与所述水力瞬变模型进行耦合；
S3、当迭代次数为1时，对耦合所述水力瞬变模型的所述多水力设施联合优化调控模型进行求解，得到原始的决策变量对应的目标函数值，形成所述多水力设施联合优化调控的父代调控方案集；
S4、根据所述父代调控方案集中各方案的决策变量和对应的目标函数值，建立所述水力瞬变模型的替代模型，进而将所述替代模型和所述多水力设施联合优化模型进行耦合；
S5、当迭代次数大于1时，将所述父代调控方案集中的所述决策变量通过优化算法生成第二代决策变量，通过耦合所述替代模型的所述水力设施联合优化调控模型进行求解，得到所述第二代决策变量对应的目标函数值，形成所述多水力设施联合优化调控的子代调控方案集；
S6、将所述父代调控方案集与所述子代调控方案集进行结合，并按照非支配等级层筛选得到新父代调控方案集；
S7、将所述父代调控方案集与所述子代调控方案集进行结合并选择出非支配方案集，将所述非支配方案集代入所述水力瞬变模型进行校验，得到校验非支配方案集，将得到的所述校验非支配方案集并入步骤S3形成的所述父代调控方案集作为训练方案集，对所述替代模型进行训练，实现所述替代模型的第一次迭代；
S8、利用所述新父代调控方案集代替所述父代调控方案集，并利用耦合迭代后替代模型的所述水力设施联合优化调控模型求解，重复步骤S5—S7；直至迭代次数达到设置的最大迭代次数，得到所述多水力设施联合优化调控模型的最优调控方案集。
优选的，所述替代模型的输入值为所述训练方案集中各方案的决策变量，输出值为所述训练方案集中各方案的目标函数；所述替代模型的高斯曲面为：
Figure BDA0003207066090000031
其中，X，y分别表示输入和输出矩阵；X.表示预测输入矩阵，即用来检测替代模型的样本；f.表示预测输出向量。
优选的，所述多水力设施联合优化调控模型与所述水力瞬变模型之间的耦合方式为：将所述多水力设施联合优化调控模型生成的决策变量作为所述水力瞬变模型的输入，根据所述水力瞬变模型生成计算所述多水力设施联合优化调控模型的目标函数所需变量。
优选的，耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤：
S31、根据所述多水力设施联合优化调控模型建立需要进行调控的水力设施决策变量种群，并对所述水力设施决策变量种群中的水力设施设置约束条件，通过所述水力设施的约束条件对所述多水力设施联合优化调控模型的决策变量进行处理；
S32、驱动所述水力瞬变模型进行瞬变过程中的计算，得到计算所述多水力设施联合优化调控模型的目标函数所需变量；
S33、根据所述目标函数的所需变量，通过所述多水力设施联合优化模型的计算得到对应的目标函数值。
优选的，耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程并行处理，作为所述多水力设施联合优化调控模型的第一并行域。
优选的，耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程包括以下步骤：
S51、根据水力设施的约束条件，对所述多水力设施联合优化调控模型的第二代决策变量进行处理；
S52、驱动所述替代模型进行瞬变过程中的计算，得到计算所述多水力设施联合优化调控模型的目标函数所需变量；
S53、根据步骤S52得到的所述第二代决策变量对应的目标函数的所需变量，通过所述多水力设施联合优化模型的计算得到对应的目标函数值，形成非支配方案集。
优选的，所述第二代决策变量的生成过程与耦合所述替代模型的所述多水力设施联合优化调控模型的求解过程并行处理，作为所述多水力设施联合优化调控模型的第二并行域。
优选的，所述父代调控方案集、所述子代调控方案集、所述非支配方案集和所述校验非支配方案集均为调控方案集，所述调控方案集包括所述多水力设施联合优化调控模型的决策变量及其对应的目标函数值。
优选的，所述父代调控方案集首先通过NSGA-II优化算法进行非支配等级排序，再将所述父代调控方案集中的决策变量通过优化算法生成第二代决策变量。
优选的，驱动所述水力瞬变模型对所述非支配方案集进行的校验过程并行处理，作为所述多水力设施联合优化调控模型的第三并行域。
本发明的有益效果是：本发明公开了一种调水工程多水力设施联合优化调控方案制定方法，将水力瞬变模型和多水力设施联合优化调控模型进行耦合，得到满足约束条件的调控方案，并通过建立所述瞬变模型的替代模型对得到的调控方案不断进行迭代更新，得到多水力设施联合优化调控的最优调控方案。本发明可以有效提高水力瞬变模型模拟的准确性，提高多水力设施联合优化调控的效率和实用性，提高多水力设施联合调控方案的切实性、可操作性，从而可以及时的指导实际工程的调控运行，减少调控过程风险、避免反复调控、提高调控运行品质。
附图说明
图1是多水力设施联合优化调控模型方案制定流程图；
图2是隧洞特征线网格示意图；
图3是多水力设施联合优化调控模型与水力瞬变模型及替代模型耦合关系图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
一种调水工程多水力设施联合优化调控方案制定方法，涉及调水工程多水力设施联合优化调控技术领域，如图1所示，包括以下步骤：
S1、整理调水工程沿线水情及运行状态监测数据，进行率定并验证得到符合调水工程实际运行的水力瞬变模型；
所述调水工程沿线水情及运行状态监测数据包括管道流量、管道压力，隧洞流量、隧洞水位，调节池水位，泵站水泵前后压力、水泵过流量、水泵调频等运行状态数据，沿线阀门前后压力、开度、过流量等运行状态数据；
S11、结合调水工程系统拓扑图，利用特征线法建立调水工程管道内部节点数学模型、隧洞内部节点数学模型以及管道连接边界、隧洞连接边界、管道与隧洞连接边界、水泵边界、空气阀边界、阀门边界、水池边界、定流量、定水位边界等设备边界模型；
管道内部节点的特征线方程为：
Figure BDA0003207066090000051
Figure BDA0003207066090000052
其中，
Figure BDA0003207066090000053
Figure BDA0003207066090000054
i表示管道上的结点号；j表示时层号；Q表示管道流量，代替控制方程中的平均流速；A表示管道横截面积；Δt表示时间步长；Δx表示管道分段长，即空间步长；g表示重力加速度，f表示水力摩阻系数，H表示液体压头；
如图2隧洞特征线网格所示，隧洞内部节点的特征线方程为：
Figure BDA0003207066090000055
Figure BDA0003207066090000061
其中，
Figure BDA0003207066090000062
分别为f在P与L两点、P与R两点上的平均值；
Figure BDA0003207066090000063
Figure BDA0003207066090000064
分别为Bc-在P与L两点，Bc+在P与R两点的平均值；对于R、L点上断面的水深、流量等信息可通过同一时层已知的网格结点A、B、C三点的值的二次插值得到；
该所述隧洞内部节点的特征线满足库朗稳定条件：
Figure BDA0003207066090000065
各设备边界模型可通过结合边界自身的物理特性和管道、隧洞的特征线方程得到；
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流量出现的时刻，得到泵后阀关闭规律中决策变量泵后阀快关过程或慢关过程的开度、快关过程或慢关过程的时长以及总关闭时间的约束条件；所述多水力设施联合优化调控模型中对于其中的多个水力设施中泵后阀的约束条件为：
Figure BDA0003207066090000071
其中，tc_min表示多水力设施联合调控系统中相应的泵后阀可操控的最短时长，tc_max表示多水力设施联合调控系统中相应的泵后阀可操纵的最长时长；
且当所述泵后阀的关闭顺序为先快关后慢关时，所述泵后阀的关闭规律为：
Figure BDA0003207066090000072
当所述泵后阀的关闭顺序为先慢关后快关时，所述泵后阀的关闭规律为：
Figure BDA0003207066090000073
对于所述多水力设施联合优化调控模型中水泵的约束条件为：nmin>-1.2nr；
对于所述多水力设施联合优化调控模型中管道的约束条件为：
Hmax≤Hs_max＝1.3Hw_max
Hmin>-7.5m
其中，Hw_max表示沿线稳态时的最大压力值，Hs_max表示瞬变过程允许出现的最大压力，nr表示水泵的额定转速，Hmax、Hmin分别表示系统在瞬变过程中沿线最大压力和最小压力，nmin表示泵站机组事故停电时水泵的最大倒转速；
对于所述多水力设施联合优化调控模型中控制断面运行水位的约束条件为：
Figure BDA0003207066090000081
其中，Zqs max和Zqs min分别表示控制建筑物设计规定的最高及最低运行水位，Zqmax和Zqmin分别表示瞬变过程中建筑物出现的最高及最低水位；
S24、根据所述多水力设施联合优化调控模型的决策变量、目标函数和约束条件，建立多水力设施联合优化调控模型；将所述多水力设施联合优化调控模型和所述水力瞬变模型进行耦合，如图3所示具体的耦合方式为：将所述多水力设施联合优化调控模型生成的决策变量传递给所述水力瞬变模型，作为所述水力瞬变模型的输入，利用所述水力瞬变模型进行瞬变过程计算，并将计算管道最大、最小压力，隧洞进出口、调节池最高、最低水位，机组最大倒转速等多水力设施联合优化调控模型的目标函数所需的变量传递至多水力设施联合优化调控模型，实现模型间的耦合。
所述多水力设施联合优化调控模型的求解流程如图1所示：依次进行种群初始化、种群约束条件并行处理、并行驱动水力瞬变模型进行瞬变过程计算、并行目标函数值的计算、非支配等级排序、拥挤距离计算、选择计算、交叉计算、变异过程并行计算、并行驱动所述替代模型进行瞬变过程计算、并行目标函数值的计算、生成子代、父代子代融合、新父代的生成；针对上述求解过程中的并行性进行分析，所述多水力设施联合优化调控模型在求解过程中的并行域包括种群约束条件处理、驱动水力瞬变模型进行瞬变过程计算、目标函数值的计算、驱动替代模型进行瞬变过程计算以及NSGA-II优化算法生成子代决策变量中变异计算过程。
S3、当迭代次数为1时，利用并行算法求解耦合所述水力瞬变模型的多水力设施联合优化调控模型，得到原始的决策变量对应的目标函数值，形成所述多水力设施联合优化调控的父代调控方案集；
上述耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程，具体包括以下步骤：
S31、根据所述多水力设施联合优化调控模型建立需要进行调控的水力设施决策变量种群，对所述水力设施决策变量种群进行初始化；并对所述水力设施决策变量种群中的水力设施设置约束条件，通过所述水力设施的约束条件对所述多水力设施联合优化调控模型的决策变量进行处理；
S32、驱动所述水力瞬变模型进行瞬变过程中的计算，得到计算所述多水力设施联合优化调控模型的目标函数所需变量；
S33、根据所述目标函数所需变量，通过所述多水力设施联合优化模型的计算得到对应的目标函数值；
其中步骤S31—步骤S33作为耦合所述水力瞬变模型的所述多水力设施联合优化调控模型的求解过程的并行域，作为所述多水力设施联合优化调控模型的第一并行域。
S4、根据所述父代调控方案集中各方案的决策变量和对应的目标函数值，作为水力瞬变模型的输入值和输出值，利用matlab采用高斯曲面，建立所述水力瞬变模型的替代模型，进而采用所述水力瞬变模型与所述多水力设施联合优化模型之间的耦合方式，将所述替代模型和所述多水力设施联合优化模型进行耦合；
所述替代模型的输入值为训练方案集中各方案的决策变量，输出值为训练方案集中各方案的目标函数；所述替代模型的高斯曲面为：
Figure BDA0003207066090000091
其中，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所述的调水工程多水力设施联合优化调控方案制定方法，其特征在于，所述替代模型的输入值为所述训练方案集中各方案的决策变量，输出值为所述训练方案集中各方案的目标函数；所述替代模型的高斯曲面为：
Figure FDA0003207066080000021
其中，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所述的调水工程多水力设施联合优化调控方案制定方法，其特征在于，驱动所述水力瞬变模型对所述非支配方案集进行的校验过程并行处理，作为所述多水力设施联合优化调控模型的第三并行域。
CN202110920075.1A 2021-08-11 2021-08-11 一种调水工程多水力设施联合优化调控方案制定方法 Active CN113627673B (zh)

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 true CN113627673A (zh) 2021-11-09
CN113627673B 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)

* Cited by examiner, † Cited by third party
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 河海大学 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
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
CN113627673B (zh) 2023-08-04

Similar Documents

Publication Publication Date Title
Saloux et al. Model-based predictive control to minimize primary energy use in a solar district heating system with seasonal thermal energy storage
CN107591844A (zh) 考虑节点注入功率不确定性的主动配电网鲁棒重构方法
Li et al. Identifying explicit formulation of operating rules for multi-reservoir systems using genetic programming
CN104616060A (zh) 基于bp神经网络和模糊逻辑的绝缘子污秽度预测方法
CN111080483A (zh) 基于图划分的供水管网dma自动分区方法
CN104050517A (zh) 基于grnn神经网络的光伏发电预测方法
WO2024109103A1 (zh) 一种闸门联动控制方法、装置及并联供水发电系统
Yi et al. Intelligent prediction of transmission line project cost based on least squares support vector machine optimized by particle swarm optimization
CN113850692A (zh) 一种基于深度学习的城市供水系统闸泵群优化调度方法
Jung et al. Fluid transients and pipeline optimization using GA and PSO: The diameter connection
CN113361214A (zh) 一种基于水位流量数据的明渠控制模型参数辨识方法
CN108717584B (zh) 一种供水管网多目标分区方法
Askari et al. Forecasting semi-dynamic response of natural gas networks to nodal gas consumptions using genetic fuzzy systems
CN114777192A (zh) 基于数据关联和深度学习的二级网供热自主优化调控方法
CN105389442A (zh) 一种耦合遗传算法、神经网络与数值模拟的反向设计方法
CN114219334A (zh) 一种基于遗传算法的贝叶斯网络天然气管道泄漏概率计算方法
Zhang et al. A universal stability quantification method for grid-connected hydropower plant considering FOPI controller and complex nonlinear characteristics based on improved GWO
CN105389614A (zh) 一种神经网络自更新过程的实现方法
Ferdowsi et al. A survey of PSO contributions to water and environmental sciences
CN116717839A (zh) 供热控制方法、控制装置和供热系统
CN113627673A (zh) 一种调水工程多水力设施联合优化调控方案制定方法
Pérez-Rúa et al. Closed-loop two-stage stochastic optimization of offshore wind farm collection system
Mukesh et al. Influence of search algorithms on aerodynamic design optimisation of aircraft wings
CN114909706A (zh) 一种基于强化学习算法和压差控制的二级网平衡调控方法
CN114357874A (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