CN116485163B - 一种考虑输水间歇的输水系统水量调度方案编制方法 - Google Patents
一种考虑输水间歇的输水系统水量调度方案编制方法 Download PDFInfo
- Publication number
- CN116485163B CN116485163B CN202310745222.5A CN202310745222A CN116485163B CN 116485163 B CN116485163 B CN 116485163B CN 202310745222 A CN202310745222 A CN 202310745222A CN 116485163 B CN116485163 B CN 116485163B
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 288
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims description 16
- 238000005086 pumping Methods 0.000 claims description 10
- 238000005265 energy consumption Methods 0.000 claims description 5
- 230000035772 mutation Effects 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 3
- 230000001276 controlling effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 abstract description 10
- 230000008901 benefit Effects 0.000 abstract description 8
- 238000003911 water pollution Methods 0.000 abstract description 7
- 230000007812 deficiency Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000383 hazardous chemical Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 241001092142 Molina Species 0.000 description 1
- 241001454523 Quillaja saponaria Species 0.000 description 1
- 235000009001 Quillaja saponaria Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000344 soap Substances 0.000 description 1
- 238000006467 substitution reaction 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/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。
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 CN116485163A (zh) | 2023-07-25 |
CN116485163B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117252349B (zh) * | 2023-11-20 | 2024-03-01 | 中国一冶集团有限公司 | 基于水质水量双控的城市河流智能调控方法 |
Citations (5)
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 | 中国水利水电科学研究院 | 一种考虑水力响应特性的跨流域调水工程旬水量优化调度计划编制方法 |
-
2023
- 2023-06-25 CN CN202310745222.5A patent/CN116485163B/zh active Active
Patent Citations (5)
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)
Title |
---|
南水北调东线工程江苏段多目标优化调度研究;郭玉雪;张劲松;郑在洲;方国华;薛刘宇;闻昕;;水利学报(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116485163A (zh) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sarker et al. | Progress on the demand side management in smart grid and optimization approaches | |
CN111762051B (zh) | 基于聚合商的电动汽车参与受端电网低谷调峰需求响应调控方法 | |
CN110417048B (zh) | 一种考虑源网荷约束的直流电网送受端联合调峰优化方法 | |
CN116485163B (zh) | 一种考虑输水间歇的输水系统水量调度方案编制方法 | |
CN108460482A (zh) | 考虑可中断负荷和日前小时电价的风电消纳方法 | |
CN103904664B (zh) | 一种基于有效静态安全域的agc机组实时调度方法 | |
CN110417061B (zh) | 一种基于改进蛙跳算法的电热联合系统调度方法 | |
CN113393038B (zh) | 一种综合能源站网及需求响应协同的规划方法 | |
CN102621941B (zh) | 一种冷轧电镀锌机组钢卷排序方法及系统 | |
Saravanan | DSM in an area consisting of residential, commercial and industrial load in smart grid | |
CN109494810A (zh) | 间歇性能源电站调度评价方法及该评价的限电分配方法 | |
CN109359830B (zh) | 一种水电站电力交易计划电量分解方法 | |
CN108875207A (zh) | 一种核反应堆优化设计方法及系统 | |
CN114862457A (zh) | 梯级水电站群参与日前现货市场的竞价电量边界确定方法 | |
CN114243736A (zh) | 电网优化调度方法、电子设备、存储介质及程序产品 | |
CN110729721B (zh) | 电力系统全局备用容量计算方法 | |
CN116739370A (zh) | 一种综合能源系统成本计算方法、装置及综合能源系统 | |
CN113452079B (zh) | 一种计及水电流域梯级调度的发电自驾驶系统及方法 | |
CN112186813B (zh) | 一种电力系统区域电网调度方法、系统、装置及存储介质 | |
CN113240183A (zh) | 商业建筑电热负荷日前优化调度方法和系统 | |
CN113394820A (zh) | 一种新能源并网电力系统的优化调度方法 | |
CN113205250A (zh) | 一种考虑水库调度风险的风光水互补发电计划编制方法 | |
CN115378016B (zh) | 一种多电动汽车集群日前充电计划生成方法及系统 | |
CN112257994B (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 |