CN116485163A - 一种考虑输水间歇的输水系统水量调度方案编制方法 - Google Patents

一种考虑输水间歇的输水系统水量调度方案编制方法 Download PDF

Info

Publication number
CN116485163A
CN116485163A CN202310745222.5A CN202310745222A CN116485163A CN 116485163 A CN116485163 A CN 116485163A CN 202310745222 A CN202310745222 A CN 202310745222A CN 116485163 A CN116485163 A CN 116485163A
Authority
CN
China
Prior art keywords
period
water delivery
water
scheduling
delivery system
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
CN202310745222.5A
Other languages
English (en)
Other versions
CN116485163B (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.)
Jiangsu Water Source Co ltd Of East Line Of South To North Water Transfer
Jiangsu University
China Institute of Water Resources and Hydropower Research
Original Assignee
Jiangsu Water Source Co ltd Of East Line Of South To North Water Transfer
Jiangsu University
China Institute of Water Resources and Hydropower Research
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 Jiangsu Water Source Co ltd Of East Line Of South To North Water Transfer, Jiangsu University, China Institute of Water Resources and Hydropower Research filed Critical Jiangsu Water Source Co ltd Of East Line Of South To North Water Transfer
Priority to CN202310745222.5A priority Critical patent/CN116485163B/zh
Publication of CN116485163A publication Critical patent/CN116485163A/zh
Application granted granted Critical
Publication of CN116485163B publication Critical patent/CN116485163B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/06314—Calendaring for a resource
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/12—Computing arrangements based on biological models using genetic models
    • G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • 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/06313—Resource planning in a project environment
    • 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

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Quality & Reliability (AREA)
  • Biophysics (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Educational Administration (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Physiology (AREA)
  • Molecular Biology (AREA)
  • Genetics & Genomics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)

Abstract

本发明涉及一种考虑输水间歇的输水系统水量调度方案编制方法，该方法考虑输水系统输水的间歇特性，将调度时段划分为前间歇期、连续期和后间歇期，根据该间歇特性构建调度期时间违反值计算模型，并以调度时段连续期的开始时间、结束时间和平均输水量为决策变量建立决策空间，然后以时间违反值、输水量任务违反值及总运行成本最小为目标，采用差分进化算法求解最优决策变量作为输水系统水量调度方案。该方法不仅满足了输水系统危化品运输及突发水污染治理的时间需求，还在保证完成供水任务的同时降低了调度运行成本，有效提高了水量调度的综合效益。

Description

一种考虑输水间歇的输水系统水量调度方案编制方法
技术领域
本发明涉及水资源调度技术领域，尤其涉及一种考虑输水间歇的输水系统水量调度方案编制方法。
背景技术
水量调度方案编制是输水系统开展水量优化调度的重要工作之一，但由于梯级泵站群的输水系统耦联关系复杂难以统筹调度，且梯级泵站群的联合调度需要充分考虑泵站的运行效益，导致现有的水量调度方案编制方法应用于实际调度中往往存在综合效益不高的问题，难以满足实际水量调度需求。虽然国内外专家学者已经引入、改进和提出了大量的水量优化调度模型。但是，这些水量优化调度模型尚缺少对河渠输水间歇特性的考虑，编制出的水量调度方案无法给危化品运输或突发水污染治理安排合理时间，因而难以在实际调度决策中给予水质安全层面的保障，并且经济性不好，从而导致这些方法难以有效落地，需要进一步改善。
发明内容
本发明所要解决的技术问题是提供一种可为危化品运输和水质污染处理预留时间并能提高水量调度综合效益的考虑输水间歇的输水系统水量调度方案编制方法。
本发明解决上述技术问题所采用的技术方案为：一种考虑输水间歇的输水系统水量调度方案编制方法，包括如下步骤：
S1，将调度期划分为n个调度时段，n为大于1的自然数，并将每个调度时段按时间顺序依次划分为前间歇期、连续期和后间歇期，在调度时段内的前间歇期和后间歇期内输水系统不输水而在连续期内输水；
S2，以每一个调度时段内连续期的开始时间、结束时间、输水系统的平均输水流量为决策变量构建决策空间，其中，连续期的开始时间的取值范围为对应调度时段的起始时间节点至末时间节点之间的时段，连续期的结束时间的取值范围为对应调度时段的开始时间至末时间节点之间的时段，输水系统的平均输水流量的范围为零至该输水系统的输水上限；
S3，构建考虑输水间歇特性的水量调度时间违反值计算模型：
两个相邻调度时段中，将在前调度时段的前间歇期、连续期、后间歇期分别定义为、、，将在后调度时段的前间歇期、连续期、后间歇期分别定义为、、，其中i表示调度时段序号，则：
当且，或，且，或，且时，
，上式中，为调度期内总的间歇期时长下限违反值，为重复计入的间歇期时间，为连续不 输水时长下限；
当和均为0，且，和不同时为0时，
，上式中，为调度期内总的连续输水期时长上限违反值，所述连续输水期为调度期内连续的输水 时段，为重复计入的连续期时间，为连续输水时长上限；
为调度期内时间违反值，、分别为和的权重因子；
S4，构建输水系统水量调度目标函数：
式中，为调度期内最小时间违反值，为调度期内最小的供 水量任务违反值，为调度期内实际供水量，为调度期内供水任务，为 调度期内最小的输水系统运行成本，为调度期内调度时段内的输水系统总 运行成本，n为调度期内调度时段总数；
约束条件：
输水系统的输水能力约束：
式中，为第j渠段的输水流量，为第j渠段的输水能力；
梯级泵站的输水能力约束：
式中，为第k座泵站的输水流量，为第k座泵站的输水能力；
输水系统的水量平衡：
式中，表示第j渠段的入流流量，为第j渠段的损失流量， 为第j渠段的出流流量；
S5，根据水量调度目标函数采用差分进化算法求解调度期内最优的各调度时段决策变量，即求解最优的各调度时段的连续期开始时间、结束时间以及输水系统的平均输水流量，从而形成输水系统的水量调度方案。
进一步的，所述S5包括如下步骤：
S51，随机生成决策变量的初始解并根据种群规模构建初始化种群，计算初始化种 群中个体的时间违反值、供水量任务违反值、输水系统总运行成本；
S52，通过差分进化算法的种群更新机制不断通过变异、交叉、选择操作生成新个 体，并计算新个体的、、值；
S53，将新个体加入上一代种群，按分层思想逐一对、、进行排序，即 对初始化种群中的个体先按值大小进行排序，如果存在值相同的个体，则对值相同 的个体再按值大小进行排序，如果存在值相同的个体，则对值相同的 个体再按值大小进行排序，并生成个体的适应度值，淘汰适应度相对大的个体以生成 新一代种群；
S54，判断种群代数是否达到种群终止代数，若是，则执行S55，若否，则返回执行S52；
S55，从最新一代种群中选取适应度最小的个体作为最优个体，最优个体记录的每个调度时段的连续期开始时间、结束时间及输水系统平均输水量即为输水系统的水量调度方案。
进一步的，所述步骤S4中调度期内第i调度时段内的输水系统总运行成本的计算方式如下：
上式中，表示输水系统中第k个泵站在第i调度时段内总的运行成本，G表示 输水系统中的泵站个数，、、分别表示第k个泵站在第i调度 时段的运行成本、额外成本、捞草成本，表示第k个泵站在第i调度时段的平均抽调流 量，表示第k个泵站在泵站扬程为时泵站的历史平均能源单耗，为泵站预估 扬程，为第k个泵站的电价，表示第k个泵站抽调单方水产生的额外费用，表示第k个泵站的平均捞草成本，为第i调度时段的时长。
进一步的，所述S1中划分调度时段的方式为将调度期匹配至自然月的旬，调度期匹配在每旬的天数即为一个调度时段的长度。
进一步的，所述权重因子P1、P2均设为1000。
进一步的，所述种群规模为100。
进一步的，所述种群终止代数为20000。
与现有技术相比，本发明的优点在于：本发明提供了一种考虑了输水间歇的水量调度方法，该方法将调度时段按河渠输水间歇特性划分为前间歇期、连续期和后间歇期，根据该间歇特性构建调度期时间违反值计算模型，并以调度时段连续期的开始时间、结束时间和平均输水量为决策变量建立决策空间，然后以时间违反值、输水量任务违反值及总运行成本最小为目标，采用差分进化算法求解最优决策变量作为输水系统水量调度方案。该方法不仅满足了输水系统危化品运输及突发水污染治理的时间需求，还在保证完成供水任务的同时降低了调度运行成本，有效提高了水量调度的综合效益。
附图说明
图1为调度时段分期基本形式及变式示意图。
图2为差分进化算法求解决策变量的流程示意图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
在本申请的描述中，需要说明的是，对于方位词，如有术语“中心”、“横向”、“纵向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示方位和位置关系为基于附图所示的方位或位置关系，仅是为了便于叙述本申请和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定方位构造和操作，不能理解为限制本申请的具体保护范围。术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
本发明涉及一种考虑输水间歇的输水系统水量调度方案编制方法，包括如下步骤：
S1，将调度期划分为n个调度时段，n为大于1的自然数，并将每个调度时段按时间顺序依次划分为前间歇期、连续期和后间歇期，在调度时段内的前间歇期和后间歇期内输水系统不输水而在连续期内输水；
S2，以每一个调度时段内连续期的开始时间、结束时间、输水系统的平均输水流量为决策变量构建决策空间，其中，连续期的开始时间的取值范围为对应调度时段的起始时间节点至末时间节点之间的时段，连续期的结束时间的取值范围为对应调度时段的开始时间至末时间节点之间的时段，输水系统的平均输水流量的范围为零至该输水系统的输水上限，其中输水系统的输水上限根据输水系统中的输水建筑物的输水能力确定；
S3，构建考虑输水间歇特性的水量调度时间违反值计算模型：
两个相邻调度时段中，将在前调度时段的前间歇期、连续期、后间歇期分别定义为、、，将在后调度时段的前间歇期、连续期、后间歇期分别定义为、、，其中i表示调度期内从前往后的调度时段的序号，则：
当且，或，且，或，且时，
，上式中， 为调度期内总的间歇期时长下限违反值，为重复计入的间歇期时间，为连续不输水 时长下限；需要说明的是，因为计算总的间歇期的时长下限违反值时，可能会存在对间歇期 重复计算的情况，因此需要扣除重复计入的间歇期时间，并非一个固定值，是根据相邻 调度时段中重复计入的间歇期时间确定的，一般由输水系统的运输环境、危化品运输 要求等根据经验确定；如果一个调度时段的连续期为零，则该调度时段的前间歇期时长等 于后间歇期时长并等于该调度时段时长；
当和均为0，且，和不同时为0时，
，上式中，为 调度期内总的连续输水期时长上限违反值，所述连续输水期为调度期内连续的输水时段，为重复计入的连续期时间，为连续输水时长上限；需要说明的是，因为计算总的连 续期的时长上限违反值时，可能会存在对连续期重复计算的情况，因此需要扣除重复计入 的连续期时间，并非一个固定值，是根据相邻调度时段中重复计入的连续期时间确定 的，一般由输水系统的运输环境、危化品运输要求等根据经验确定；
为调度期内时间违反值，、分别为和的权重因子；由于和具有同等重要性，因此、二者的取值相等，具体取值根据经验得到，比如可以均 取1000；
S4，构建输水系统水量调度目标函数：
式中，为调度期内最小时间违反值，为调度期内最小的供 水量任务违反值，为调度期内实际供水量，为调度期内供水任务，为 调度期内最小的输水系统运行成本，为调度期内调度时段内的输水系统总 运行成本，n为调度期内调度时段总数；
约束条件：
输水系统的输水能力约束：
式中，为第j渠段的输水流量，为第j渠段的输水能力，渠段输水能力 可根据渠段的水力参数确定；
梯级泵站的输水能力约束：
式中，为第k座泵站的输水流量，为第k座泵站的输水能力，泵站的 输水能力可根据泵站的设计参数确定；
输水系统的水量平衡：
式中，表示第j渠段的入流流量，为第j渠段的损失流量， 为第j渠段的出流流量；
S5，根据水量调度目标函数采用差分进化算法求解调度期内最优的各调度时段决策变量，即求解最优的各调度时段的连续期开始时间、结束时间以及输水系统的平均输水流量，从而形成输水系统的水量调度方案。
需要说明的是，在进行调度过程中，调度时段的前间歇期、连续期、后间歇期都可能为零，但不会同时为零，因此可以得到调度时段的基本形式和变式，如图1所示。构建计算时间违反值计算模型的目的在于，使得输水系统为危化品运输提供足够的运输时间以及在应对突发水质污染能够有足够的时间治理。因此，从这个角度看，需要输水系统的连续输水时间越短越好，而不输水的时间越长越好，也即需要使得间歇期时间越长越好，连续期时间越短越好，为此，在编制水量调度方案时，考虑对调度方案的时间违反值进行计算。
需要进一步说明的是，在调度期最后调度时段时，如果没有危化品运输任务而又 需要迫切地完成输水任务时，可以将最后调度时段的时长设置大一些，并同时将最后调度 时段的连续输水时长上限值设置大一些，由此可以针对调度期最后调度时段计算得出 对应的时间违反值，从而得到调度期内的时间违反值为于、、与对应权重 因子乘积之和，各权重因子相同并根据经验进行具体取值。
需要进一步说明的是，上述S1中划分调度时段的方式为将调度期匹配至自然月的旬，调度期匹配在每旬的天数即为一个调度时段。比如调度期为1月5日至2月12日，则划分的调度时段分别为1月5日-1月10日、1月11日-1月20日、1月21日至1月31日、2月1日至2月10日、2月11日-2月12日。
在本实施例中，如图2所示，上述S5中根据水量调度目标函数采用差分进化算法求解调度期内最优的各调度时段决策变量，具体包括如下步骤：
S51，随机生成决策变量的初始解并根据种群规模构建初始化种群，种群规模可以 根据实际计算效率的需要根据经验进行选择，比如100等，计算初始化种群中个体的时间违 反值、供水量任务违反值、输水系统总运行成本；
S52，通过差分进化算法的种群更新机制不断通过变异、交叉、选择操作生成新个 体，并计算新个体的、、值；
S53，将新个体加入上一代种群，按分层思想逐一对、、进行排序，即 对初始化种群中的个体先按值大小进行排序，如果存在值相同的个体，则对值相同 的个体再按值大小进行排序，如果存在值相同的个体，则对值相同的 个体再按值大小进行排序，并生成个体的适应度值，淘汰适应度相对大的个体以生成 新一代种群；
S54，判断种群代数是否达到种群终止代数，种群代数可以根据全局收敛速度的要求根据经验进行确定，比如20000代，若是，则执行S55，若否，则返回执行S52；
S55，从最新一代种群中选取适应度最小的个体作为最优个体，最优个体记录的每个调度时段的连续期开始时间、结束时间及输水系统平均输水量即为输水系统的水量调度方案。
由于差分进化算法及变异、交叉、选择等操作均属于现有技术，本实施例中仅简要列出计算步骤，本领域技术人员在知晓采用该方法求解连续期开始时间、结束时间、及输水系统平均输水量的决策变量后，应当知道具体的实施步骤，不在此详述。
在本实施例中，步骤S4中调度期内第i调度时段内的输水系统总运行成本的计算方式如下：
上式中，表示输水系统中第k个泵站在第i调度时段内总的运行成本，G表示 输水系统中的泵站个数，、、分别表示第k个泵站在第i调度 时段的运行成本、额外成本、捞草成本，表示第k个泵站在第i调度时段的平均抽调流 量，表示第k个泵站在泵站扬程为时泵站的历史平均能源单耗，为泵站预估 扬程，为第k个泵站的电价，表示第k个泵站抽调单方水产生的额外费用，表示第k个泵站的平均捞草成本，为第i调度时段的时长，由于泵站的运行成本计 算属于现有技术，本申请并未对其进行改进，不在此进行详细描述。
下面以南水北调东线工程江苏段洪泽湖-骆马湖段输水系统为例，对本发明的水量调度方案编制方法进行进一步说明。南水北调东线工程江苏段洪泽湖-骆马湖段输水系统包括徐洪河、中运河两条输水线路，其中徐洪河输水线路包括泗洪站、睢宁二站、沙集站、邳州站四个泵站，中运河输水线路包括泗阳站、刘老涧二站、皂河二站三个泵站，各泵站参数如表1和表2所示。为体现本发明考虑输水间歇特性的水量调度方案的优越性，选取未考虑河渠输水间歇的水量调度方案进行对比，通过比较两种方法的水量调度方案的效果，对本发明的实施过程和效果进行陈述。
表1 徐洪河输水线路各梯级泵站参数表
梯级泵站 预估扬程（m） 能源单耗（kw/(kt*m)） 电费单价元（KW.h） 单方水额外费（元/m³） 捞草费用单价（元/天） 抽调流量上限（m³/s） 至下游泵站渠段损失流量（m³/s）
泗洪站 2.35 4.56 0.6164 0 619.2263 120 4.97
睢宁二站 6.7 4.41 0.6164 0 1535.525 60 3.46
沙集站 6.7 4.6 0.484 0.006 1872.948 50 3.46
邳州站 1.7 4 0.6164 0 784.3957 100 \
表2 中运河输水线路各梯级泵站参数表
梯级泵站 预估扬程（m） 能源单耗（kw/(kt*m)） 电费单价元（KW.h） 单方水额外费（元/m³） 捞草费用单价（元/天） 抽调流量上限（m³/s） 至下游泵站渠段损失流量（m³/s）
泗阳站 5.55 3.83 0.6164 0 216.6667 164 6.24
刘老涧二站 3.7 4.06 0.6164 0 1413.483 80 8.45
皂河二站 4.45 3.91 0.6164 0 0 75 \
水量调度计划按0.5亿m³、1亿m³、1.5亿m³、2亿m³选取4组工况，分别定义为工况一至工况四，调度期为2022年11月13日至2023年1月13日，调度时段按旬划分。假设连续不输水时长下限为5天，连续输水时长上限为15天，权重因子p1、p2均为1000，差分进化算法种群规模100，变异系数F=0.4，交叉系数CR=0.3，终止代数20000代。
将调度期按旬划分为7个调度时段即7旬，各调度时段的天数为8、10、10、10、11、10、3，此即是未考虑河渠输水间歇的运行天数。采用本发明的水量调度方案编制方法得出最优的调度方案，针对不同工况可得如表3至表6中第二列所示的指标数据，未考虑河渠输水间歇方法进行编制的不同工况下的水量调度方案对应的指标数据如表3至表6第三列所示。需要说明的是，表3至表6中的指标计算方式如下：输水效率=（实际供水量-水量损失）/实际供水量，缺水率=（计划供水量-实际供水量）/计划供水量，这些指标计算属于现有技术，不在此详细展开。
表3 两种方法在工况一下的计算结果对比表
项目 本发明方法 无间歇方法
计划供水量（亿m³） 0.5 0.5
实际供水量（亿m³） 0.5003 0.5003
水量损失（亿m³） 0.0801 0.4443
总成本（元） 488029 870431
缺水率(%) -0.06 -0.06
输水效率(%) 84 11
表4 两种方法在工况二下的计算结果对比表
项目 本发明方法 无间歇方法
计划供水量（亿m³） 1 1
实际供水量（亿m³） 0.9996 0.9996
水量损失（亿m³） 0.1457 0.4443
总成本（元） 971690 1272853
缺水率(%) 0.04 0.04
输水效率(%) 85 56
表5 两种方法在工况三下的计算结果对比表
项目 本发明方法 无间歇方法
计划供水量（亿m³） 1.5 1.5
实际供水量（亿m³） 1.4999 1.4999
水量损失（亿m³） 0.2331 0.4443
总成本（元） 1496646 1675971
缺水率(%) 0.007 0.007
输水效率(%) 84 70
表6 两种方法在工况四下的计算结果对比表
项目 本发明方法 无间歇方法
计划供水量（亿m³） 2 2
实际供水量（亿m³） 2.002 2.002
水量损失（亿m³） 0.2404 0.4443
总成本（元） 2061908 2079089
缺水率(%) -0.01 -0.01
输水效率(%) 88 78
由表3至表6可以看出，本发明方法、无间歇方法在4组工况下的缺水率基本相同；而本发明方法的水量损失普遍更小，因而本发明方法的输水效率更高；本发明方法在四组工况中的供水成本均低于无间歇方法，而且无间歇方法无法给危化品运输及水质污染治理预留时间，而本发明方法可以利用输水间歇为危化品运输和突发水污染治理预留时间，具有更好的综合效益。
因此，本发明方法与无间歇方法相比，能够在保障供水任务完成的同时，大幅度提升输水效率，整体上降低供水成本，并能够为危化品运输及突发水污染治理预留足够的时间，有利于提升实际水量调度的综合效益，为输水系统的水量调度方案编制提供了一种新的思路。
尽管已经示出和描述了本发明的实施例，本领域技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变形，本发明的范围由权利要求及其等同物限定。

Claims (7)

1.一种考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：包括，
S1，将调度期划分为n个调度时段，n为大于1的自然数，并将每个调度时段按时间顺序依次划分为前间歇期、连续期和后间歇期，在调度时段内的前间歇期和后间歇期内输水系统不输水而在连续期内输水；
S2，以每一个调度时段内连续期的开始时间、结束时间、输水系统的平均输水流量为决策变量构建决策空间，其中，连续期的开始时间的取值范围为对应调度时段的起始时间节点至末时间节点之间的时段，连续期的结束时间的取值范围为对应调度时段的开始时间至末时间节点之间的时段，输水系统的平均输水流量的范围为零至该输水系统的输水上限；
S3，构建考虑输水间歇特性的水量调度时间违反值计算模型：
两个相邻调度时段中，将在前调度时段的前间歇期、连续期、后间歇期分别定义为、、，将在后调度时段的前间歇期、连续期、后间歇期分别定义为、、，其中i表示调度时段序号，则：
当且，或，且，或，且时，
，式中，为调度期内总 的间歇期时长下限违反值，为重复计入的间歇期时间，为连续不输水时长下限；
当和均为0，且，和不同时为0时，
，式中，为调度期 内总的连续输水期时长上限违反值，所述连续输水期为调度期内连续的输水时段，为重 复计入的连续期时间，为连续输水时长上限；
为调度期内时间违反值，、分别为和的权重因子；
S4，构建输水系统水量调度目标函数：
式中，为调度期内最小时间违反值，为调度期内最小的供水量 任务违反值，为调度期内实际供水量，为调度期内供水任务，为调度 期内最小的输水系统运行成本，为调度期内调度时段内的输水系统总运行 成本，n为调度期内调度时段总数；
约束条件：
输水系统的输水能力约束：
式中，为第j渠段的输水流量，为第j渠段的输水能力；
梯级泵站的输水能力约束：
式中，为第k座泵站的输水流量，为第k座泵站的输水能力；
输水系统的水量平衡：
式中，表示第j渠段的入流流量，为第j渠段的损失流量，为第 j渠段的出流流量；
S5，根据水量调度目标函数采用差分进化算法求解调度期内最优的各调度时段决策变量，即求解最优的各调度时段的连续期开始时间、结束时间以及输水系统的平均输水流量，从而形成输水系统的水量调度方案。
2.根据权利要求1所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：所述S5包括如下步骤：
S51，随机生成决策变量的初始解并根据种群规模构建初始化种群，计算初始化种群中 个体的时间违反值、供水量任务违反值、输水系统总运行成本；
S52，通过差分进化算法的种群更新机制不断通过变异、交叉、选择操作生成新个体，并 计算新个体的、、值；
S53，将新个体加入上一代种群，按分层思想逐一对、、进行排序，即对初 始化种群中的个体先按值大小进行排序，如果存在值相同的个体，则对值相同的个 体再按值大小进行排序，如果存在值相同的个体，则对值相同的个体 再按值大小进行排序，并生成个体的适应度值，淘汰适应度相对大的个体以生成新一 代种群；
S54，判断种群代数是否达到种群终止代数，若是，则执行S55，若否，则返回执行S52；
S55，从最新一代种群中选取适应度最小的个体作为最优个体，最优个体记录的每个调度时段的连续期开始时间、结束时间及输水系统平均输水量即为输水系统的水量调度方案。
3.根据权利要求2所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在 于：所述步骤S4中调度期内第i调度时段内的输水系统总运行成本的计算方式 如下：
上式中，表示输水系统中第k个泵站在第i调度时段内总的运行成本，G表示输水 系统中的泵站个数，、、分别表示第k个泵站在第i调度时段 的运行成本、额外成本、捞草成本，表示第k个泵站在第i调度时段的平均抽调流量，表示第k个泵站在泵站扬程为时泵站的历史平均能源单耗，为泵站预估扬 程，为第k个泵站的电价，表示第k个泵站抽调单方水产生的额外费用，表示第k个泵站的平均捞草成本，为第i调度时段的时长。
4.根据权利要求1所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：所述S1中划分调度时段的方式为将调度期匹配至自然月的旬，调度期匹配在每旬的天数即为一个调度时段的长度。
5.根据权利要求3所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：所述权重因子P1、P2均设为1000。
6.根据权利要求2所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：所述种群规模为100。
7.根据权利要求2所述的考虑输水间歇的输水系统水量调度方案编制方法，其特征在于：所述种群终止代数为20000。
CN202310745222.5A 2023-06-25 2023-06-25 一种考虑输水间歇的输水系统水量调度方案编制方法 Active CN116485163B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310745222.5A CN116485163B (zh) 2023-06-25 2023-06-25 一种考虑输水间歇的输水系统水量调度方案编制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310745222.5A CN116485163B (zh) 2023-06-25 2023-06-25 一种考虑输水间歇的输水系统水量调度方案编制方法

Publications (2)

Publication Number Publication Date
CN116485163A true CN116485163A (zh) 2023-07-25
CN116485163B CN116485163B (zh) 2023-08-18

Family

ID=87212321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310745222.5A Active CN116485163B (zh) 2023-06-25 2023-06-25 一种考虑输水间歇的输水系统水量调度方案编制方法

Country Status (1)

Country Link
CN (1) CN116485163B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117252349A (zh) * 2023-11-20 2023-12-19 中国一冶集团有限公司 基于水质水量双控的城市河流智能调控方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110097342A (ko) * 2010-02-25 2011-08-31 (주)씨에스 씨엔아이 배수지의 지능화 운전 시스템 및 방법
CN109325610A (zh) * 2018-08-15 2019-02-12 东莞理工学院 一种面向多水源梯级泵站调水系统的分段线性调度方法
CN109598408A (zh) * 2018-10-29 2019-04-09 华中科技大学 一种兼顾用水公平性和重要性的年水量调度计划编制方法
CN111709134A (zh) * 2020-06-12 2020-09-25 河海大学 一种基于闸泵站群控制的多水源跨流域调度方法及装置
CN112785087A (zh) * 2021-02-22 2021-05-11 中国水利水电科学研究院 一种考虑水力响应特性的跨流域调水工程旬水量优化调度计划编制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110097342A (ko) * 2010-02-25 2011-08-31 (주)씨에스 씨엔아이 배수지의 지능화 운전 시스템 및 방법
CN109325610A (zh) * 2018-08-15 2019-02-12 东莞理工学院 一种面向多水源梯级泵站调水系统的分段线性调度方法
CN109598408A (zh) * 2018-10-29 2019-04-09 华中科技大学 一种兼顾用水公平性和重要性的年水量调度计划编制方法
CN111709134A (zh) * 2020-06-12 2020-09-25 河海大学 一种基于闸泵站群控制的多水源跨流域调度方法及装置
CN112785087A (zh) * 2021-02-22 2021-05-11 中国水利水电科学研究院 一种考虑水力响应特性的跨流域调水工程旬水量优化调度计划编制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭玉雪;张劲松;郑在洲;方国华;薛刘宇;闻昕;: "南水北调东线工程江苏段多目标优化调度研究", 水利学报, no. 11 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117252349A (zh) * 2023-11-20 2023-12-19 中国一冶集团有限公司 基于水质水量双控的城市河流智能调控方法
CN117252349B (zh) * 2023-11-20 2024-03-01 中国一冶集团有限公司 基于水质水量双控的城市河流智能调控方法

Also Published As

Publication number Publication date
CN116485163B (zh) 2023-08-18

Similar Documents

Publication Publication Date Title
Sarker et al. Progress on the demand side management in smart grid and optimization approaches
CN111762051B (zh) 基于聚合商的电动汽车参与受端电网低谷调峰需求响应调控方法
CN116485163B (zh) 一种考虑输水间歇的输水系统水量调度方案编制方法
CN110417048B (zh) 一种考虑源网荷约束的直流电网送受端联合调峰优化方法
CN108460482A (zh) 考虑可中断负荷和日前小时电价的风电消纳方法
CN106650105A (zh) 一种混流泵叶轮的设计方法
CN112785087B (zh) 一种考虑水力响应特性的跨流域调水工程旬水量优化调度计划编制方法
CN107248015A (zh) 基于集合预报的水库实时蓄水调度方法
CN114548546A (zh) 一种调水工程水量的优化调度方法
CN115115277A (zh) 园区综合能源调度方法、装置、计算机设备及存储介质
CN110417061B (zh) 一种基于改进蛙跳算法的电热联合系统调度方法
Skworcow et al. Pump schedules optimisation with pressure aspects in complex large-scale water distribution systems
CN103699938A (zh) 一种含抽水蓄能电站的电力系统发电计划制定方法
Saravanan DSM in an area consisting of residential, commercial and industrial load in smart grid
CN115313378A (zh) 一种风光储发电站日前有功出力优化调度方法及系统
CN109359830B (zh) 一种水电站电力交易计划电量分解方法
CN104463414A (zh) 一种分时电价下热轧生产经济负荷调度方法
CN109508813B (zh) 一种兼顾调峰和弃水的水电电量分解模型建立方法
CN108365637A (zh) 一种送电计划与抽水储能发电计划优化方法及系统
CN113393038B (zh) 一种综合能源站网及需求响应协同的规划方法
Haimes et al. Coordination of regional water resource supply and demand planning models
CN107748930B (zh) 一种充分灌溉条件下直接补渠的单库-多站系统水资源优化配置方法
CN112186813B (zh) 一种电力系统区域电网调度方法、系统、装置及存储介质
CN115496378A (zh) 一种计入风能减排效益的电力系统经济调度方法
Du et al. Optimization of pump scheduling in waterworks considering load balancing using improved genetic 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