发明名称 一种含抽水蓄能电站的电力系统发电计划制定方法
摘要 本发明专利公开了一种含抽水蓄能电站的电力系统发电计划制定方法,其特征在于包括以下步骤:首先制定检修计划,获取电力系统负荷曲线;然后在负荷平衡、机组出力与爬坡限制、系统旋转备用和抽水蓄能电站库容约束条件下,以调度周期内所有火电机组的发电成本与启停费用之和最小为目标建立模型,用混沌控制的基本粒子群算法求解。本发明有效克服了求解电力系统发电计划模型时收敛性不佳、容易陷入局部最优的问题,制定更合理的电力系统发电计划。
申请公布号 CN103699938A 申请公布日期 2014.04.02
申请号 CN201310469126.9 申请日期 2013.10.10
申请人 华北电力大学(保定);中国南方电网有限责任公司 发明人 谢红玲;李燕青;孙凯航;董驰;李翔;王坚;梁志飞;傅志伟;付妍;沈博一
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I;G06N3/00(2006.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 石家庄新世纪专利商标事务所有限公司 13100 代理人 董金国
主权项 1.一种含抽水蓄能电站的电力系统发电计划制定方法,其特征在于包括以下步骤:A. 制定检修计划:确定电力系统参数、火力发电机组数量及参数、抽水蓄能电站数量及参数、各抽水蓄能电站中的抽水蓄能发电机组数量及参数;所述电力系统参数包括最低旋转备用容量;所述火力发电机组参数包括单位发电成本、启停费用、额定容量、最小技术出力百分数、爬坡速率;所述抽水蓄能电站参数包括最小、最大库容电量;所述抽水蓄能发电机组参数包括单位发电成本、启停费用、额定容量;B. 获取电力系统负荷曲线:所述电力系统负荷预测数据时间间隔固定为Δ<i>t</i>,单位为小时;在时间-电力系统负荷预测值坐标系下将所述电力系统负荷预测数据相邻的两点用直线连接,形成电力系统负荷曲线;所述电力系统负荷预测数据将电力系统负荷曲线划分为T个调度区间;C. 建立抽水蓄能电站发电计划模型:              目标函数为:<img file="301512DEST_PATH_IMAGE001.GIF" wi="526" he="41" />(1)其中,T为调度时段数;<img file="256960DEST_PATH_IMAGE002.GIF" wi="18" he="24" />为第n火力发电机组的单位发电成本,MW,<img file="190281DEST_PATH_IMAGE003.GIF" wi="17" he="24" />为第n火力发电机组的启停费用,元/kWh;<img file="796843DEST_PATH_IMAGE004.GIF" wi="24" he="24" />为第t调度时段,第n火力发电机组的发电功率,MW,<img file="242737DEST_PATH_IMAGE005.GIF" wi="35" he="24" />为第t时段,第n火力发电机组的状态,取值为0时停运,为1时投运;M为抽水蓄能发电机组台数;<img file="731487DEST_PATH_IMAGE006.GIF" wi="20" he="24" />为第m抽水蓄能发电机组的单位发电成本,元/kWh,<img file="671761DEST_PATH_IMAGE007.GIF" wi="24" he="24" />为第t调度时段,第m抽水蓄能发电机组的发电功率,MW;<img file="880632DEST_PATH_IMAGE008.GIF" wi="36" he="24" />为第t时段,第m抽水蓄能发电机组的状态,取值为0时停运,为1时投运;<img file="44897DEST_PATH_IMAGE009.GIF" wi="18" he="24" />为第m抽水蓄能发电机组的启停费用,万元/次;       约束条件包括:       负荷平衡约束:<img file="224206DEST_PATH_IMAGE010.GIF" wi="173" he="45" />(2)   其中,<i>p</i><sub><i>L,t</i></sub>为第t时段电力系统负荷预测值,MW;      火力发电机组出力约束<img file="951859DEST_PATH_IMAGE011.GIF" wi="300" he="24" />(3)其中,<img file="329751DEST_PATH_IMAGE012.GIF" wi="36" he="24" />为第n火力发电机组的最大技术出力,单位;<img file="868180DEST_PATH_IMAGE013.GIF" wi="36" he="24" />为第n火力发电机组的最小技术出力,MW;      抽水蓄能发电机组出力约束:<img file="19938DEST_PATH_IMAGE014.GIF" wi="299" he="24" />(3)<img file="302015DEST_PATH_IMAGE015.GIF" wi="39" he="24" />为第m抽水蓄能发电机组的最大技术出力,MW;<img file="799992DEST_PATH_IMAGE013.GIF" wi="36" he="24" />为第m抽水蓄能发电机组的最小技术出力,MW;火力发电机组爬坡限制约束:<img file="493010DEST_PATH_IMAGE016.GIF" wi="339" he="24" />(4)其中,<img file="270080DEST_PATH_IMAGE017.GIF" wi="27" he="24" />为第n火力发电机组的爬坡速度,计算方法为:<img file="339536DEST_PATH_IMAGE018.GIF" wi="92" he="24" />(5)所述(5)式中<img file="629703DEST_PATH_IMAGE019.GIF" wi="17" he="24" />为第n火力发电机组的爬坡速率,%;       系统旋转备用约束:<img file="995088DEST_PATH_IMAGE020.GIF" wi="465" he="46" />(6)       其中,backup<sub>t</sub>为第t调度时段的最低旋转备用容量,MW;                       抽水蓄能电站库容约束:<img file="433022DEST_PATH_IMAGE021.GIF" wi="183" he="120" />(7)其中:<img file="56902DEST_PATH_IMAGE022.GIF" wi="27" he="24" />为第t调度时段末第k抽水蓄能电站的库容电量,MW;<img file="450843DEST_PATH_IMAGE023.GIF" wi="27" he="24" />为第k抽水蓄能电站库容电量的初始值,MW;<img file="564292DEST_PATH_IMAGE024.GIF" wi="39" he="24" />、<img file="427206DEST_PATH_IMAGE025.GIF" wi="39" he="24" />为第k抽水蓄能电站的最小、最大库容电量,MW;       D:采用改进的粒子群优化算法求解所述抽水蓄能电站发电计划模型的最优解,包括以下分步骤;D1:确定初始粒子群TN:所述粒子为T*NN二维矩阵,NN=N+M,其元素为火力发电机组或抽水蓄能发电机组在各调度区间的发电功率,MW,所述初始粒子群由从多于T*NN粒子中选择目标函数值最小的T*NN个粒子组成,所述各粒子的速度和位置在混沌序列控制下产生;<img file="337000DEST_PATH_IMAGE026.GIF" wi="186" he="99" />(8)D2:依次更新所述各粒子的速度和位置:按(8)式更新所述各粒子的速度:<img file="664076DEST_PATH_IMAGE027.GIF" wi="281" he="24" />(8)其中,<img file="886110DEST_PATH_IMAGE028.GIF" wi="17" he="24" />为更新前第k粒子的速度,<img file="220008DEST_PATH_IMAGE029.GIF" wi="18" he="24" />为更新前第k粒子的位置,MW,<img file="185690DEST_PATH_IMAGE030.GIF" wi="27" he="24" />为更新后第k粒子的速度,<img file="367273DEST_PATH_IMAGE031.GIF" wi="36" he="24" />为第k粒子的最优个体极值,MW,<img file="510940DEST_PATH_IMAGE032.GIF" wi="35" he="24" />为所有粒子的最优全局极值,MW,r为0~1的随机数;c<sub>1</sub>和c<sub>2</sub>为加速因子;w<sub>k</sub>为惯性系数;按(9)式更新所述各粒子的位置:<img file="348446DEST_PATH_IMAGE033.GIF" wi="95" he="24" />(9)其中,<img file="914557DEST_PATH_IMAGE034.GIF" wi="27" he="24" />为更新后第k粒子的位置,MW;D3:更新各粒子的个体最优位置:依次计算更新后各粒子的目标函数值,若小于所述粒子的个体极值,用其更新所述粒子的个体极值;D4:更新全局最优位置:目标函数值最小粒子个体最优位置作为更新后的全局最优位置;D5:混沌扰动更新各粒子位置,如果所述更新后粒子位置对应的目标函数值小于更新前粒子位置对应的目标函数值,则更新所述粒子位置,否则所述粒子位置不变;D6:对新的粒子进行约束处理,如出现越限后,超过最大值取上边界,超过最小值取下边界,使所有的粒子都满足约束条件;D7:判断是否达到预定迭代次数,如果没有,转向步骤D2,否则,输出最优全局极值及对应的最小目标函数值。
地址 071003 河北省保定市北市区永华北大街619号