CN113128768A - 考虑风电不确定性的水风火短期优化调度方法 - Google Patents

考虑风电不确定性的水风火短期优化调度方法 Download PDF

Info

Publication number
CN113128768A
CN113128768A CN202110435073.3A CN202110435073A CN113128768A CN 113128768 A CN113128768 A CN 113128768A CN 202110435073 A CN202110435073 A CN 202110435073A CN 113128768 A CN113128768 A CN 113128768A
Authority
CN
China
Prior art keywords
wind
power
output
thermal power
water
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
Application number
CN202110435073.3A
Other languages
English (en)
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.)
PowerChina Zhongnan Engineering Corp Ltd
Original Assignee
PowerChina Zhongnan Engineering Corp Ltd
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 PowerChina Zhongnan Engineering Corp Ltd filed Critical PowerChina Zhongnan Engineering Corp Ltd
Priority to CN202110435073.3A priority Critical patent/CN113128768A/zh
Publication of CN113128768A publication Critical patent/CN113128768A/zh
Pending legal-status Critical Current

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
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • 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
    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

本发明公开了一种考虑风电不确定性的水风火短期优化调度方法，在满足水风火约束的基础上，构建了清洁能源利用率最大和火电机组出力波动最小的多目标优化调度模型，该模型利用梯级水电启停灵活、爬坡速度快等优点平抑风电出力波动。由于梯级水电、风电、火电联合优化调度问题是个高度复杂的混合整数非线性规划问题，整体求解难度较大。本发明基于分层求解思想，将模型分解为风电调度层、水电调度层与火电调度层，层与层之间以余荷相关联，并提出了集成神经网络、区间估计、改进的遗传算法与启发式算法相结合的求解框架，实现了模型的快速求解。

Description

考虑风电不确定性的水风火短期优化调度方法
技术领域
本发明涉及水风火电调度领域，特别是一种考虑风电不确定性的水风火短期优化调度方法。
背景技术
风电等可再生能源具有随机性、波动性与反调峰性，风电等可在生能源的大规模并网严重影响电力系统的安全性、可靠性与稳定性。在可再生能源发展的新趋势下，减少弃风，提高电网对风电等清洁能源的消纳，是目前电力系统发展的关键任务。
弃风的主要原因是电力系统的灵活性电源不足，电力系统没有充足的灵活性电源平抑风电波动，部分地区仅仅依靠火电应对风电出力的不确定性。随着风电等清洁性能源大规模政策并网，部分地区仅仅依靠火电平抑风电波动，火电机组将频繁启停或运行在深度调峰区，严重影响火电机组的寿命，严重影响电力系统的安全。随着我国西北、西南梯级水电站群的并网，如何充分梯级水电等灵活性电源的调节能力，制定合理的水风火等多电源联合调度方案，以缓解风电等间歇性电源的大规模馈入带来的运行风险，是当前亟需解决的问题之一。
目前已有的梯级水电、火电与风电联合调度，大都考虑确定性调度，并不符合风电随机波动的特点。
本发明中的“短期”是指调度周期为一天，每15min一个点，共计96个点进行计算。
发明内容
本发明所要解决的技术问题是，针对现有技术不足，提供一种考虑风电不确定性的水风火短期优化调度方法，充分利用梯级水电的灵活性平抑风电波动，实现风电、水电多发，火电少发，计算速度快、时效性强，有效减少系统煤耗，增加系统效益。
为解决上述技术问题，本发明所采用的技术方案是：一种考虑风电不确定性的水风火短期优化调度方法，包括以下步骤：
S1、以t时刻风速预报数据、历史实际风电功率数据为神经网络的输入，预测t时刻风电出力
Figure BDA0003032860840000021
t＝1,2，…，T；
S2、根据下式平抑风电所需要的最小出力
Figure BDA0003032860840000022
Figure BDA0003032860840000023
其中，
Figure BDA0003032860840000024
Figure BDA0003032860840000025
为
Figure BDA0003032860840000026
的反函数，
Figure BDA0003032860840000027
Ωt(j)＝{bi|i＝t,t+y,...,t+ny}，n表示误差集合Ωt(j)中的元素个数；ξ为给定误差水平；bi为排序后风电出力偏差值；num表示bi≤ξ的个数；y表示时间间隔(y＝15min)，β为误差；
S3、将火电机组总负荷减去平抑风电所需要的最小出力和梯级水电站出力，得到的结果分配到火电机组，并对火电机组开机进行优化。
本发明在满足水电、风电、火电约束的基础上，充分利用梯级水电的灵活性平抑风电波动，实现风电、水电多发，火电少发，计算速度快、时效性强，有效减少系统煤耗，增加系统效益。
步骤S3中，所述梯级水电站出力的计算过程包括：
I)选取M组水库运行过程中水位序列，利用所述水位序列初始化种群，得到M个个体，每个个体表示为Zr＝[z1r,z2r,...,zTr],zTr为个体Zr在T时刻的染色体，表示个体Zr在T时刻的水位值；T表示总的时段长度(T＝96,每个时段长度为15min)；r表示个体编号1≤r≤M；
II)按发电量最大为目标进行评价，计算每个个体的适应度；
III)对所有个体的适应度排序，利用确定的选择概率Ps(取Ps＝0.6)，选择Ps×M个适应度最高的个体，直接进入下一代；被淘汰的(1-Ps)×M个个体由交叉或变异而产生的新的个体所替代；
IV)重复步骤I)-III)，直到最优个体的适应度不发生改变时，结束；最优个体对应的适应度值，即为梯级水电站出力。
步骤I)～IV)提出的遗传算法具有群体搜索的特性。它的搜索过程是从一个具有多个个体的初始群体开始的，一方面可以有效地避免搜索一些不必搜索的点，另一方面基于概率规则，而不是确定性规则，这使得搜索灵活，效率高。
由于遗传算法易受初始解的影响，步骤I)中，采用Logistic映射初始化种群，即采用公式xr+1＝μxr(1-xr)进行映射，xr表示采用Logistic映射后生成的第r个个体；μ为正向可调参数，μ可取3.8。
为了使Logistic生成的个体在可行域区内，需要将生成的序列，映射到原可行域区间内，优选地，按式
Figure BDA0003032860840000031
进行映射，
Figure BDA0003032860840000032
z分别表示水位变量上限、下限；Zr(1)表示映射后个体的初始水位序列。
为了进一步改善遗传算法的性能，步骤III)中，对种群中的所有个体进行两两随机交叉配对，采用随机数的方式生成交叉点，依据设定的交叉概率Pc在其交叉点处相互交换两个个体的染色体，从而产生两个新的个体；交叉概率按照如下公式进行计算：
Figure BDA0003032860840000033
其中，Pc为交叉概率，fc为要交叉的两个个体中最大的适应度值，fmax为种群中最大适应度值，fave为每代群体的平均适应度值，Pc1、Pc2为正向可调参数。
步骤III)中，为了保证交叉过程中产生的新个体与原个体差异性，在随机交叉配对过程中，当公式
Figure BDA0003032860840000034
成立时，个体Z1与Z2执行交叉操作，否则重新配对；Z1、Z2为两个不同的个体，cn为总的迭代次数，cj为当前迭代的次数，b为阈值(本发明中，可取b＝0.8)。
为了进一步改善遗传算法的性能，步骤III)中，对交叉后的所有个体的染色体，以概率为Pm进行突变，突变后的染色体值大小满足水位最高与最低约束；Pm按照如下公式进行计算：
Figure BDA0003032860840000035
其中，Pm为变异概率，fmax为种群中最大的适应度值，fave为每代群体的平均适应度值，fm要变异个体的适应度值，Pm1与Pm2为正向可调参数。
步骤S3的火电机组开机优化具体实现过程包括：
A)在不考虑启停成本的基础上，计算火电机组的单位煤耗成本，并根据单位煤耗成本的大小，建立机组开机顺序表；
B)剔除所述开机顺序表中检修、尚未达到开机时间的机组，得到更新后的开机顺序表；
C)根据更新后的开机顺序表，选择单位电量消耗煤量最小的机组，并累加其最大出力，直到其大于系统净负荷(电网负荷减去风电与水电出力)最大值；
D)将步骤C)中各火电机组最小出力累加，并验证其是否小于系统净负荷最小值，若不满足则根据步骤A)的开机顺序表，调整机组，并验证步骤C)；是否成立，若不成立重新安排火电站机组，并调整火电机组承担负荷。
基于步骤A)-D)的火电机组开停机方法，能在系统负荷一定的情况下尽可能减少火电机组的系统煤耗，增加火电机组发电效益。
还包括：
E)计算各个时段是否满足火电机组爬坡约束，若不满足，则重新选择火电机组，重复步骤A)～步骤E)，直到满足火电站爬坡约束；所述火电爬坡约束为：
Figure BDA0003032860840000041
Figure BDA0003032860840000042
为火电机组j在t时刻的出力，
Figure BDA0003032860840000043
Δj 分别表示向上、向下爬坡速度。
基于步骤E)的火电机组爬坡验证方法，能在一定的条件内快速、准确识别火电机组出力是否满足爬坡约束，使得火电机组出力满足实际需求。
还包括：
S4、判断公式
Figure BDA0003032860840000044
是否成立，若是，则表明梯级水电能平抑风电波动，水电出力满足
Figure BDA0003032860840000045
否则，梯级水电不足以抵消风电出力的不确定性，火电机组在满足各种约束的条件下，承担平抑风电波动出力ζt；
S5、判断
Figure BDA0003032860840000051
是否成立，若成立则按照步骤S3对梯级水电站负荷进行分配；否则，在满足各种约束条件的前提下以步长Δε增加ζt，直到
Figure BDA0003032860840000052
成立，按照步骤S3对梯级水电站负荷进行分配，并安排火电机组，Nsd表示水电站的个数，
Figure BDA0003032860840000053
表示水电站m在t时刻的出力。
本发明还提供了一种考虑风电不确定性的水风火短期优化调度系统，其包括计算机设备；所述计算机设备被配置或编程为用于执行上述方法的步骤。
与现有技术相比，本发明所具有的有益效果为：本发明在满足水风火约束的基础上，构建了清洁能源利用率最大和火电机组出力波动最小的多目标优化调度模型，该模型利用梯级水电启停灵活、爬坡速度快等优点平抑风电出力波动。其技术方案为：由于梯级水电、风电、火电联合优化调度问题是个高度复杂的混合整数非线性规划问题，整体求解难度较大。基于分层求解思想，将模型分解为风电调度层、水电调度层与火电调度层，层与层之间以余荷相关联，并提出了集成神经网络、区间估计、改进的遗传算法与启发式算法相结合的求解框架，实现了模型的快速求解。以某区域为例进行模拟计算，算例表明该方法时效性强，仅需89.2s，即可得出结果，实现了风电、水电等清洁能源多发、火电少发，增加了系统效益。
附图说明
图1为本发明实施例总体计算原理；
图2为本发明实施例实际系统负荷与实际风电出力；
图3为本发明实施例置信区间为90％风电出力；
图4(a)和图4(b)分别为梯级水电优化前后实验结果图；
图5为优化前后梯级水电站出力结果；
图6为优化前后火电出力结果图；
图7(a)优化前调度结果；
图7(b)优化后调度结果。
具体实施方式
本发明实施例构建了清洁能源消纳最大与火电机组出力波动最小的多目标优化调度模型。为了实现模型的快速求解，基于分层求解的思想将模型分解为三层，层与层之间以余荷相关联(如图1所示)，并根据梯级水电、风电与火电的不同特性，制定了不同的调度策略。首先根据风速预报数据、历史风电功率数据与历史风速数据对风电出力进行预测，并给定置信水平，采用经验法计算风电出力的置信区间；根据风电出力区间，安排梯级水电出力过程使得水电出力尽可能提供充足的灵活性，以应对置信区间内的风电波动，同时提供更多的调峰电量，尽可能的保证火电出力平稳；火电机组调节能力差，因此应尽量保持余留给火电的净负荷尽可能平稳，使之承担基荷。在求解过程中，采用Elman神经网络与非参数法对风电运行层置信区间进行计算；采用改进的遗传算法对梯级水电站的负荷进行分配；采用启发式算法安排火电机组出力。需要说明的是采用该分层策略求解得到的解，虽然可能不是全局最优解，但求解速度较快，在可接受的时间内求得较优解。步骤如下：
步骤1：水-风-火联合优化调度模型的总体目标是在保证弃风、弃水最小的前提下，尽量保证火电出力平稳。具体目标函数如下：
Figure BDA0003032860840000061
式中f1用来反映火电出力的平稳情况，若f1越大表示火电出力越不平稳；
Figure BDA0003032860840000062
表示火电机组在t时刻的出力；f2为调度模型中弃风、弃水量；
Figure BDA0003032860840000063
分别表示弃风、弃水总功率；T为调度周期时段数。
在求解过程中，应当满足以下约束：
1)电力平衡约束
在求解过程中应保证电力系统供需平衡：
Figure BDA0003032860840000064
其中，
Figure BDA0003032860840000065
表示t时刻系统预测负荷需求；
Figure BDA0003032860840000066
与
Figure BDA0003032860840000067
分别表示梯级水电与风电机组预测出力。
2)梯级水电站约束
水力约束：
梯级水电上下游水力联系很强，上游出库直接影响下游入库；此外梯级水电还应当满足水位库容约束、出库流量约束。
Figure BDA0003032860840000071
其中：
Figure BDA0003032860840000072
表示水电站m在t时刻的库容，
Figure BDA0003032860840000073
表示m的上游电站在t-τ时刻的出库流量；τ为上游出库到达下游时间；
Figure BDA0003032860840000074
与
Figure BDA0003032860840000075
分别表示梯级水电站m在时刻t的区间入流、出库、蒸发、渗漏损失的水量；
Figure BDA0003032860840000076
Figure BDA0003032860840000077
Zm,t 分别表示电站m出库下限、出库上限、水位上限、水位下限。
运行约束：
水电站在运行过程中需满足机组过流上下限约束、出力大小上下限约束、最小启停时间约束。
Figure BDA0003032860840000078
式中
Figure BDA0003032860840000079
表示机组i最小、最大过流；
Figure BDA00030328608400000710
为机组i在t时刻开停机状态参数开机为0关机为1；
Figure BDA00030328608400000711
表示水电机组i的最大出力；
Figure BDA00030328608400000712
表示机组开机、关机持续时间；
Figure BDA00030328608400000713
分别表示最小运行时长与停机时长。
3)火电机组约束
火电机组出力上下限约束：
Figure BDA00030328608400000714
其中，
Figure BDA00030328608400000715
当启停状态改变时为1其他时段为0；
Figure BDA00030328608400000716
表示火电机组j在t时刻的出力；
Figure BDA00030328608400000717
表示火电机组j的最小与最大出力；Nsh表示火电机组个数。
火电机组启停与爬坡约束：
由于火电机组启停时间较长，不考虑火电机组的日内启停，即火电机组发电状态一旦确定，在日内不会发生变化。
Figure BDA0003032860840000081
火电机组爬坡约束：
Figure BDA0003032860840000082
式中
Figure BDA0003032860840000083
Δh,j 分别表示向上、向下爬坡速度。
4)备用约束
假设风电出力预测误差区间为
Figure BDA0003032860840000084
负荷功率预测误差为
Figure BDA0003032860840000085
则系统备用需求为：
Figure BDA0003032860840000086
式中
Figure BDA0003032860840000087
分别表示t时刻的正、负旋转备用。
步骤2：基于Elman神经网络与非参数法的风电运行层
采用Elman神经网络与非参数法对风电运行层的出力区间置信上下限进行求解，具体步骤如下：
1)以t时刻风速预报数据、历史风电功率数据为输入，选取动态Elman神经网络(Elman Neural Network、Elman)模型作为基础预测工具，预测t时刻风电出力
Figure BDA0003032860840000088
t＝1,2,...,T，T＝96，t与t+1的时间间隔为15min；
2)：假设t时刻风电出力预测值与实际值的偏差为bt，按照下式求得误差分布函数
Figure BDA0003032860840000089
Figure BDA00030328608400000810
Ωt(j)＝{bi|i＝t,t+y,...,t+ny} (10)
其中：n表示误差集合Ωt(j)中的所有元素个数；ξ为给定误差水平；bi为排序后风电出力偏差值；num表示bi≤ξ的个数；其中t表示时间，y表示时间间隔。
3)给定ξ，根据
Figure BDA0003032860840000091
即可采用经验法求得相应误差水平下的风电功率预测区间；假设误差为β对应的一个1-β的概率预测区间为：
Figure BDA0003032860840000092
其中
Figure BDA0003032860840000093
为
Figure BDA0003032860840000094
的反函数，
Figure BDA0003032860840000095
为t时刻风电功率预测值，为了方便下文表示，记区间下限
Figure BDA0003032860840000096
区间上限
Figure BDA0003032860840000097
步骤3：根据风电运行层风电出力上限与出力下限计算平抑风电所需要的最小出力：
Figure BDA0003032860840000098
步骤4：将步骤3)所得的
Figure BDA0003032860840000099
分配到梯级水电调度层，采用梯级水电出力平抑风电波动，梯级水电出力计算方法，见步骤5)。
步骤5：采用遗传算法对梯级水电出力进行计算，使得每个时段的水电出力
Figure BDA00030328608400000910
均大于
Figure BDA00030328608400000911
采用遗传算法对梯级电站出力计算步骤如下：
1)初始解(个体)的生成：选取M(取M＝50)组水库运行过程中水位序列，利用所述水位序列初始化种群，得到M个个体，每个个体可表示为Zr＝[z1r,z2r,...,zTr],式中zTr为个体Zr的染色体，表示在T时刻的水位值；T表示总的时段长度(T＝96,每个时段长度为15min)；r表示个体编号1≤r≤M。由于遗传算法易受初始解的影响，采用Logistic映射初始化种群，见式(13)，从而提升算法搜索效率。
xr+1＝μxr(1-xr) (13)
式中：xr表示采用Logistic映射后生成的第r个个体；μ为正向可调参数。
研究表明，Logistic映射的参数在一定范围内是混沌映射，可取u＝3.8。为了使Logistic生成的个体在可行域区内，需要将生成的序列，映射到原可行域区间内，可按式
Figure BDA00030328608400000912
进行映射(
Figure BDA00030328608400000913
y分别表示原优化变量上限、下限；Zr(1)表示映射后个体的初始位置)。
2)适应度评价：按发电量最大为目标进行评价(在保证
Figure BDA0003032860840000101
与期末水位约束的基础上尽可能的增加梯级水电出力)，计算每个个体的适应度。
3)选择运算：采用保存最优个体策略，假定群体中的个体规模为M(M＝50)，首先根据个体的适应度值进行排序；然后利用确定的选择概率Ps(取Ps＝0.6)，从母体种群中选择Ps×M个适应较高的个体,直接进入下一代；而被淘汰的(1-Ps)×M个体由交叉或变异而产生的新个体所替代。
4)交叉操作：对种群中的所有个体进行两两随机配对，采用随机数的方式生成交叉点，依据设定的交叉概率Pc在其交叉点处相互交换两个个体的染色体，从而产生两个新的个体。Pc值对产生新个体至关重要，Pc过大导致遗传算法的适应度高的个体结果很快被破坏、偏小就不容易产生新的个体结构导致搜索停滞不前。交叉概率按照如下公式进行计算：
Figure BDA0003032860840000102
式中fc为要交叉的两个个体适应度较大的适应度值；fmax为种群中最大的适应度值；fave为每代群体的平均适应度值；Pc1、Pc2为正向可调参数，可取Pc1＝0.9，Pc2＝0.6。
为了保证交叉过程中产生的新个体与原个体差异性，在随机交叉配对过程中，假设个体Z1与Z2可用向量Z1＝[z11,z21...,zT1]与Z2＝[z12,z22...,zT2]表示，在设置交叉中应采取如下措施：
Figure BDA0003032860840000103
式中：cn为计算过程中总的迭代次数，cj为当前迭代的次数，b为阈值(取b＝0.8)。
5)变异操作：对交叉后的所有个体的染色体，以概率为Pm进行突变，突变后的染色体值大小，应满足水位最高与最低约束。Pm值的大小的选取对算法性能至关重要，Pm偏小不容易产生新的个体，偏大就变成了随机搜索。Pm按照如下公式进行：
Figure BDA0003032860840000111
fm要变异个体的适应度值Pm1与Pm2为正向可调参数，可取Pm1＝0.1，Pm2＝0.001
6)重复步骤1)-5)直到最优个体的适应度不发生改变为止。最优个体对应的适应度值，即为梯级水电站出力值。
将余荷(总负荷减去风电与梯级水电出力)分配到火电机组，并按照步骤6)对火电机组开机进行优化。
步骤6：基于启发式算法的火电调度层
采用启发式算法寻求机组的最优开机方式，具体步骤如下：
1)：在不考虑启停成本的基础上，计算火电机组的单位煤耗成本，并根据煤耗成本的大小，建立机组开机顺序表。单位电量煤耗成本可按下式计算：
Figure BDA0003032860840000112
式中βj为火电机组j生产单位电量所消耗的煤耗成本；
Figure BDA0003032860840000113
表示火电机组j的当前出力；aj、bj、cj为火电机组j的运行参数为固定值。
2)：剔除开机顺序表中检修、尚未达到开机时间的机组。
3)：根据开停机顺序表，选择单位电量消耗煤量最小的机组，并累加其最大出力，直到其大于系统净负荷(电网负荷减去风电与水电出力)的最大值。
4)：将3)中各火电机组最小出力累加，并验证其是否小于系统净负荷最小值，若不满足则根据1)的开机顺序表，调整机组，并验证3)是否成立，若不成立重新安排火电站机组，并调整火电机组承担负荷。
5)：计算各个时段是否满足火电站爬坡约束，若不满足，则重新选择火电机组执1)-5)，直到满足要求。
步骤7：在满足各种约束的条件下，按照式(18)计算梯级水电是否能平抑风电波动。
Figure BDA0003032860840000121
步骤8：若式(18)成立，表明梯级水电能平抑风电波动，在负荷低谷时段，为了减少风电的随机波动对电力系统的影响同时在负荷高峰期预留更多的水电出力应保证
Figure BDA0003032860840000122
步骤9：若不满足式(18)，则梯级水电不足以抵消风电出力的不确定性，此时火电机组在满足各种约束的条件下，应承担平抑风电波动出力ζt。
步骤10：计算
Figure BDA0003032860840000123
是否成立。
a)若成立则按照步骤5对梯级水电站负荷进行分配，按步骤6安排火电机组。
b)若不成立则在满足各种约束条件的前提下以较小的步长Δε＝0.1增加ζt，直到
Figure BDA0003032860840000124
成立。并按照步骤5对梯级水电站负荷进行分配，按步骤6安排火电机组。
为了验证本发明方法的有效性，以我国西南某地区冬季典型日短期调度为例进行计算。该地区冬季风电出力波动性较大，以冬季为某日进行计算更能反映风电出力不确定性对电力系统的影响。算例包含距离不足20km的二级梯级水电站，上游水电站为多年调节电站装机为1200MW，下游电站为日调节电站装机为1320MW，装机总容量为13500MW的风电场群，以及30台装机总容量为13650MW的火电机组群。选取调度时段为15min，调度周期取1天，在调度周期内实际风电出力与系统实际负荷如图2所示，可以看出风电出力具有反调峰性，增加了电力系统的调度难度，图3为置信区间为90％时风电场群实际出力与预测出力过程。
梯级水电在调度过程中各个参数如表1所示，其中水电站1处于上游，为多年调节电站；电站2处于下游电站，为径流式电站。采用改进的遗传算法迭代500次，初始种群为50。所有仿真程序均用pathon语言编写，运行环境为4核CUP、主频3.2GHZ、内存为16GB、硬盘500GB、Windows系统的联想计算机。仿真计算100次，最长计算时间为89.2s，可有效满足短期调度时效性的需求，取最优计算结果进行分析。
表1梯级水电站运行参数
Figure BDA0003032860840000131
优化前与优化后各个电源发电量如表2所示，从表中可以看出优化前与优化后置信区间内的风电均全额消纳，但优化后火电少发45.37kWh，实现了风电、水电等清洁能源多发、火电少发。
表2优化前后各电源出力结果
Figure BDA0003032860840000132
优化前(实际调度过程)与优化后梯级水电站水位如图4(a)与4(b)所示，可以看出优化前后上游水电站由于是多年调节电站水位变化不大。而下游水电站是日调节水电站，调度期初水位被抬升，耗水率更小，更有利于发电。优化前后梯级水电站的出力如图5所示，可以看出，优化前与优化后梯级水电在任何时段都能“平抑”风电波动，即任意时段梯级水电站出力满足
Figure BDA0003032860840000133
优化前梯级水电站出力较大，后期可供电站支配的调峰水量较小，后期水电调峰能力不足，导致火电机组深度调峰，调峰压力较大；优化后梯级水电在前期出力较小，后期可供支配的水量较多，调峰出力相对较大，使得火电机组出力更加平稳。与优化前相比水电多发45.56万kWh，显著提高了发电效益。
通过计算火电机组的开机顺序表如表3所示，火电出力优化前与优化后对比如图4所示，可以看出优化后火电机组出力峰谷差为2426.12MW，与优化前相比峰谷差减小了682.15MW；优化后标准差为790.43MW，与优化前相比减小了226.75MW，火电机组出力更加平稳；优化后火电机组开机台数为8台，与优化前相比较小三台，煤耗成本减小23.33万元。火电机组优化前后出力结果如图6所示。
表3启发式算法确定火电机组经济性顺序表
Figure BDA0003032860840000141
表4优化前后火电出力对比
Figure BDA0003032860840000142
图7(a)与图7(b)展示了各个电源优化前后出力情况，可以看出无论是优化前还是优后，火电机组出力占比较大，优化后火电出力占比为70.42％，减少了0.11％。优化后梯级水电在负荷高峰期(40-52,64-88)出力加大，在负荷低谷时段出力较小；火电在负荷高峰期调峰负荷减小，负荷低谷时出力增加，火电出力更加平稳。从优化结果来看，梯级水电、风电、火电联合调度能有效减小风电出力不确定性对电力系统的影响，使得火电机组出力更加平稳。实现了清洁能源多发、火电少发的目的，获得了较好的调度结果。

Claims (10)

1.一种考虑风电不确定性的水风火短期优化调度方法，其特征在于，
包括以下步骤：
S1、以t时刻风速预报数据、历史实际风电功率数据为神经网络的输入，预测t时刻风电出力Pt pred；
S2、根据下式平抑风电所需要的最小出力Pt py：
Figure FDA0003032860830000011
其中，
Figure FDA0003032860830000012
Figure FDA0003032860830000013
为Ft *(ξ)的反函数，
Figure FDA0003032860830000014
Ωt(j)＝{bi|i＝t,t+y,...,t+ny}，n表示误差集合Ωt(j)中的元素个数；ξ为给定误差水平；bi为排序后风电出力偏差值；num表示bi≤ξ的个数；y表示时间间隔，β为误差；
S3、将火电机组总负荷减去平抑风电所需要的最小出力和梯级水电站出力，得到的结果分配到火电机组，并对火电机组开机进行优化。
2.根据权利要求1所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤S3中，所述梯级水电站出力的计算过程包括：
I)选取M组水库运行过程中水位序列，利用所述水位序列初始化种群，得到M个个体，每个个体表示为Zr＝[z1r,z2r,...,zTr]；zTr为个体Zr在T时刻的染色体，表示个体Zr在T时刻的水位值；T表示总的时段长度；r表示个体编号1≤r≤M；
II)按发电量最大为目标进行评价，计算每个个体的适应度；
III)对所有个体的适应度排序，利用确定的选择概率Ps，选择Ps×M个适应度最高的个体，直接进入下一代；被淘汰的(1-Ps)×M个个体由交叉或变异而产生的新的个体所替代；
IV)重复步骤I)-III)，直到最优个体的适应度不发生改变时，结束；最优个体对应的适应度值，即为梯级水电站出力。
3.根据权利要求2所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤I)中，采用Logistic映射初始化种群，即采用公式xr+1＝μxr(1-xr)进行初始化种群中的个体，xr表示采用Logistic映射后生成的第r个个体；μ为正向可调参数；优选地，按式
Figure FDA0003032860830000021
进行映射，
Figure FDA0003032860830000022
z分别表示水位变量上限、下限；Zr(1)表示映射后个体Zr初始水位序列。
4.根据权利要求2所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤III)中，对种群中的所有个体进行两两随机交叉配对，采用随机数的方式生成交叉点，依据设定的交叉概率Pc在其交叉点处相互交换两个个体的染色体，从而产生两个新的个体；交叉概率按照如下公式进行计算：
Figure FDA0003032860830000023
其中，Pc为交叉概率，fc为要交叉的两个个体中最大的适应度值，fmax为种群中最大适应度值，fave为每代群体的平均适应度值，Pc1、Pc2为正向可调参数。
5.根据权利要求4所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤III)中，在随机交叉配对过程中，当公式
Figure FDA0003032860830000024
成立时，个体Z1与Z2执行交叉操作，否则重新配对；Z1、Z2为两个不同的个体，cn为总的迭代次数，cj为当前迭代的次数，b为阈值。
6.根据权利要求2所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤III)中，对交叉后的所有个体的染色体，以概率为Pm进行突变，突变后的染色体值大小满足水位最高与最低约束；Pm按照如下公式进行计算：
Figure FDA0003032860830000031
其中，Pm为变异概率，fmax为种群中最大的适应度值，fave为每代群体的平均适应度值，fm要变异个体的适应度值，Pm1与Pm2为正向可调参数。
7.根据权利要求1所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，步骤S3的具体实现过程包括：
A)在不考虑启停成本的基础上，计算火电机组的单位煤耗成本，并根据单位煤耗成本的大小，建立机组开机顺序表；
B)剔除所述开机顺序表中检修、尚未达到开机时间的机组，得到更新后的开机顺序表；
C)根据更新后的开机顺序表，选择单位电量消耗煤量最小的机组，并累加其最大出力，直到其大于系统净负荷(电网负荷减去风电与水电出力)最大值；
D)将步骤C)中各火电机组最小出力累加，并验证其是否小于系统净负荷最小值，若不满足则根据步骤A)的开机顺序表，调整机组，并验证步骤C)是否成立，若不成立重新安排火电站机组，并调整火电机组承担负荷。
8.根据权利要求7所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，还包括：
E)计算各个时段是否满足火电站爬坡约束，若不满足，则重新选择火电机组，重复步骤A)～步骤E)，直到满足火电站爬坡约束；所述火电爬坡约束为：
Figure FDA0003032860830000032
Figure FDA0003032860830000033
为火电机组j在t时刻的出力，
Figure FDA0003032860830000034
Δj 分别表示向上、向下爬坡速度。
9.根据权利要求1所述的考虑风电不确定性的水风火短期优化调度方法，其特征在于，还包括：
S4、判断公式
Figure FDA0003032860830000041
是否成立，若是，则表明梯级水电能平抑风电波动，水电出力满足
Figure FDA0003032860830000042
否则，梯级水电不足以抵消风电出力的不确定性，火电机组在满足各种约束的条件下，承担平抑风电波动出力ζt；
S5、判断
Figure FDA0003032860830000043
是否成立，若成立则按照步骤S3对梯级水电站负荷进行分配；否则，在满足各种约束条件的前提下以步长Δε增加ζt，直到
Figure FDA0003032860830000044
成立，并按照步骤S3对梯级水电站负荷进行分配，并安排火电机组，Nsd表示水电站的个数，
Figure FDA0003032860830000045
表示水电站m在t时刻的出力。
10.一种考虑风电不确定性的水风火短期优化调度系统，其特征在于，包括计算机设备；所述计算机设备被配置或编程为用于执行权利要求1～9之一所述方法的步骤。
CN202110435073.3A 2021-04-22 2021-04-22 考虑风电不确定性的水风火短期优化调度方法 Pending CN113128768A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110435073.3A CN113128768A (zh) 2021-04-22 2021-04-22 考虑风电不确定性的水风火短期优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110435073.3A CN113128768A (zh) 2021-04-22 2021-04-22 考虑风电不确定性的水风火短期优化调度方法

Publications (1)

Publication Number Publication Date
CN113128768A true CN113128768A (zh) 2021-07-16

Family

ID=76779063

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110435073.3A Pending CN113128768A (zh) 2021-04-22 2021-04-22 考虑风电不确定性的水风火短期优化调度方法

Country Status (1)

Country Link
CN (1) CN113128768A (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003050987A (ja) * 1994-11-24 2003-02-21 Matsushita Electric Ind Co Ltd 最適化調整方法と最適化調整装置
CN101599870A (zh) * 2009-06-30 2009-12-09 西北工业大学 网络链路性能测量方法
CN104504229A (zh) * 2014-09-19 2015-04-08 杭州电子科技大学 一种基于混合启发式算法的智能公交调度方法
CN105303450A (zh) * 2015-11-09 2016-02-03 北京工业大学 基于谱聚类改进交叉的复杂网络社区发现方法
CN106130079A (zh) * 2016-07-08 2016-11-16 华北电力大学 一种风水火短期联合优化调度方法
CN107358231A (zh) * 2017-06-09 2017-11-17 东北电力大学 一种基于sift算子和混沌遗传算法的车牌字符识别方法
CN110212579A (zh) * 2019-06-17 2019-09-06 国网山西省电力公司电力科学研究院 一种风-水-火联合鲁棒机组组合方法
CN110472825A (zh) * 2019-07-09 2019-11-19 贵州黔源电力股份有限公司 一种多级调度机构协调的梯级水电站实时调度弃水削减方法
CN111476474A (zh) * 2020-04-01 2020-07-31 贵州黔源电力股份有限公司 梯级水电站减少弃水量的调度方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003050987A (ja) * 1994-11-24 2003-02-21 Matsushita Electric Ind Co Ltd 最適化調整方法と最適化調整装置
CN101599870A (zh) * 2009-06-30 2009-12-09 西北工业大学 网络链路性能测量方法
CN104504229A (zh) * 2014-09-19 2015-04-08 杭州电子科技大学 一种基于混合启发式算法的智能公交调度方法
CN105303450A (zh) * 2015-11-09 2016-02-03 北京工业大学 基于谱聚类改进交叉的复杂网络社区发现方法
CN106130079A (zh) * 2016-07-08 2016-11-16 华北电力大学 一种风水火短期联合优化调度方法
CN107358231A (zh) * 2017-06-09 2017-11-17 东北电力大学 一种基于sift算子和混沌遗传算法的车牌字符识别方法
CN110212579A (zh) * 2019-06-17 2019-09-06 国网山西省电力公司电力科学研究院 一种风-水-火联合鲁棒机组组合方法
CN110472825A (zh) * 2019-07-09 2019-11-19 贵州黔源电力股份有限公司 一种多级调度机构协调的梯级水电站实时调度弃水削减方法
CN111476474A (zh) * 2020-04-01 2020-07-31 贵州黔源电力股份有限公司 梯级水电站减少弃水量的调度方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
严太山;杨堤;张立波;江海升;郑晓琼;柴宏博;张超;: "基于改进小生境遗传算法的输电网安全经济规划模型", 水电能源科学, no. 07, pages 206 - 210 *
刘建文;丁洁玉;潘坤;张晓强;: "基于个体相似度的改进自适应遗传算法研究", 青岛大学学报(工程技术版), no. 01, pages 16 - 19 *
周艳聪;孙晓晨;余伟翔;: "基于改进遗传算法的物流配送路径优化研究", 计算机工程与科学, no. 10, pages 118 - 122 *
尹君;薛福珍;: "一种综合改进型遗传算法及其在多变量解耦中的应用", 化工自动化及仪表, no. 02, pages 16 - 20 *
柳伯超;秦茂玲;刘弘;: "基于遗传算法的多约束背包问题求解方案", 山东师范大学学报(自然科学版), no. 01, pages 19 - 22 *
邱柳钦;: "改进遗传算法在MF-TDMA资源规划中的研究", 通信技术, no. 04, pages 117 - 120 *
金新峰: "大规模风电馈入下富水电地区多能源短期互补协调研究", 中国优秀硕士学位论文全文数据库（电子期刊）工程科技Ⅱ辑C042-1101, pages 105 - 106 *
闻骥骏,张学明,沈成武: "避免近亲繁殖的遗传算法", 交通科技, no. 05, pages 60 - 62 *

Similar Documents

Publication Publication Date Title
Wei et al. Short-term optimal operation of hydro-wind-solar hybrid system with improved generative adversarial networks
CN109347151B (zh) 一种新能源参与调峰的送端电网电源结构优化方法
CN103151803B (zh) 一种含风电系统机组及备用配置的优化方法
Zhang et al. A hybrid particle swarm optimization with small population size to solve the optimal short-term hydro-thermal unit commitment problem
CN105790265B (zh) 一种考虑交流潮流约束的不确定性机组组合模型及求解方法
CN103580061B (zh) 微电网运行方法
CN108879787B (zh) 一种含风电的电网随机调度优化模型与方法
CN111555281A (zh) 一种电力系统灵活性资源配置仿真方法及装置
CN112381375B (zh) 一种基于潮流分配矩阵的电网经济运行域快速生成方法
CN109103929A (zh) 基于改进动态克里丝金模型的配电网经济优化调度方法
CN105244870B (zh) 一种快速计算电网风电场弃风率及机组发电量的方法
CN105305423A (zh) 一种考虑间歇性能源不确定性的最优误差边界的确定方法
CN109546689B (zh) 一种适用于大规模系统的两阶段机组组合运行控制方法
CN115940292A (zh) 基于圆圈搜索算法的含风储电力系统优化调度方法及系统
Zhou et al. Multiobjective scheduling method for short-term peak shaving operation of cascade hydro plants
CN115114854A (zh) 一种虚拟电厂分布式资源两级自组织优化聚合方法及系统
Najafi-Ghalelou et al. Risk-Constrained Scheduling of Energy Hubs: A Stochastic $ p $-Robust Optimization Approach
CN115051388A (zh) 一种基于分布鲁棒的“源-网-荷-储”两阶段调度优化方法
CN106451424A (zh) 含大规模光伏发电并网的输电网随机规划方法
CN117674290A (zh) 一种基于多场景的水电平抑分布鲁棒优化方法
Lopes et al. Metaheuristic methods applied to the pumps and turbines configuration design of water pumped storage systems
CN113128768A (zh) 考虑风电不确定性的水风火短期优化调度方法
CN115271244A (zh) 一种基于两阶段分布鲁棒优化的梯级水电站短期调峰模型
CN111967646A (zh) 一种虚拟电厂可再生能源优化配置方法
Miracle et al. Hybrid Metaheuristic Model for Optimal Economic Load Dispatch in Renewable Hybrid Energy System

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210716