研究与开发 基于Gurobi的抽水蓄能电站优化运行分析 张津玮 顾秀芳 (内蒙古工业大学, 内蒙古0101O0) 摘要 抽水蓄能机组加入系统后,可以利用自身调峰填谷的能力,提高系统的经济性,减少 煤耗量。本文基于某电网的目前数据,以系统运行成本最小为目标,考虑系统中各项约束条件建 立优化模型,探讨了抽水蓄能电站为系统带来的好处。实验以Matlab为平台,编译程序调用Gurobi 软件进行求解。优化结果表明,抽水蓄能电站可以很好的跟踪负荷的变化,降低煤耗,提高系统 的经济性。同时验证了Gurobi优化软件在电力系统运行优化方面的可行性和实用性。 关键词:Gurobi;Matlab;抽水蓄能;优化运行 Analysis on optimized operation of pumped storage power station based on gurobi ZhangJinwei GuXiufang (Inner Mongolia University of Technology,Inner Mongolia 0 1 0 1 O0) Abstract After pumped storage units join the system,it can use its own ability of adjusting the peak valley to improve the system economy,and reduce coal consumption.Based on the data of a certain power grid,this paper discusses the advantages of the pumped storage power station with the minimum operating cost of the system and the optimized model is built considering the constraints of the system. The experiment based on Matlab,and the compiler calls Gurobi software to solve.The optimization results show that the pumped storage power station can track the change of the load well,reduce the coal consumption and improve the economy of the system.At the same time,it verifies the feasibility and practicability of Gurobi optimization software in the optimal operation of power system. Keywords:gurobi;matlab;pumped storage;optimized operation 抽水蓄能机组是一种既能发电又能储存电能的 装置,起停迅速,运行灵活可靠,具有削峰填谷和 快速跟踪…等功能,还具有显著的节能环保效益『2]。 除此之外,抽水蓄能机组还是系统中承担调峰、调 频、调相、调压、旋转备用、事故备用和黑起动【3] 的重要技术手段,维护系统的安全稳定运行[4 】。 在国内外抽水蓄能机组经济运行的研究中,绝 模拟算法,通过实例验证了模型可以准确的模拟、 分析和优化系统的运行状况。文献[8]构建了考虑经 济成本和环保成本的数学模型,并且在一般引力搜 索算法的基础上引入混沌算法,应用这种改进的引 力搜索算法进行实例计算。文献[9】以系统收益最大 为目标分别对抽水蓄能电站的日调节和周调节方式 建立了混合整数线性规划模型,使用动态规划法求 解,对两种调节方式分别进行规划。 大多数学者都以系统的成本最/Ix[ 罐]或者收益最大【 。 ] 为目标进行规划。文献[6]最先以一个日调度期内的 现有文献对抽水蓄能电站的运行优化大都着重 考虑系统的经济效益,本文借鉴已有文献,以全网 电力系统运行费用最小为目标,包括燃料费和固定 运行费,使用动态规划法求解,讨论了抽蓄电站容 量在承担旋转备用和调峰填谷任务中该如何划分的 问题。文献[7]基于随机生产模拟中的负荷曲线分解 运行成本最低为目标,考虑系统各项约束条件建立 数学模型。针对所建立的模型,由于模型规模较大, 包含的变量和约束条件数量较多,将Matlab作为平 台,调用Gurobi软件进行抽蓄电站优化运行的规划。 技术和动态规划思想,提出优化模型和相应的概率 22 l胃I.I藏柬2018 ̄2 研究与开发 1 抽蓄电站优化运行建模 1.1 目标函数 以上时,机组爬坡约束可不考虑【H]。 (1)抽水蓄能电站库容约束 Wmi ≤ ≤ 。 ∑ + 以全网运行成本最小为目标,包括火电机组煤 耗成本和抽蓄机组的起停成本,即 (6) (7) i ” ≤ 4一 ≤8m ∑ minF= ∑∑ ( ,,)+ , K (1) 一 式(6)表示任意时段水库的库容在最小和最大 库容之间, i 、 分别为水库的储能上下限; = 式中, 是煤价; 为调度的时段数;t为调度时段 式(7)了水库库容在每目的始末变化范围, 的索引号;n为火电机组的数目;i为火电机组的索 引号;P ,是火电机组i在f时刻的出力;S“是抽蓄 机组的起停成本;K是起停次数; (户“)是机组f在 时段t的发电成本函数,一般由出力的二次函数表 示,即 f,(Pi,,)=airj+ , +Ci (2) 式中,系数a 、b 、C 分别为火电机组i的煤耗量特 性系数。 1.2模型约束条件 1)系统需求及相关运行约束 (1)功率平衡约束 要求每个时段火电机组出力和抽蓄机组出力之 和减去抽蓄机组抽水的功率与此时段的负荷保持 平衡,即 式中, 、 分别为抽蓄机组 在t时刻的发电 和抽水出力; .,为时段,的负荷大小。 (2)火电机组最大、最小出力约束 . in≤ .f≤ . ax (4)式中,Pi, i Pf‘m 分别为电站最小、最大出力。 (3)火电机组爬坡速度约束 At≤ f+1一 ,f≤UiAt (5) ,式中,D 、 分别为第f台火电机组有功出力的下 降速率和上升速率; 为一个时间段的时长。 2)抽水蓄能机组相关约束 抽蓄机组受自身固有的特性以及电站运行方面 的条件和【12-13],具备自己特定的运行特性和约 束。抽蓄机组响应速度较快,机组起停及工况转换 等能在很短时间内完成,当优化的时段问隔为0.5h 8mi 、 分别为每天首末时段最小和最大变动库 容,大小取可调用库容的5%【 1。 (2)抽水蓄能电站上水库电量平衡约束 +I= 一 + ‘ir3p (8) 式中,rig、rip分别为发电效率和抽水效率。式(8) 表示相邻时段水库库容的变化关系。 (3)抽水蓄能机组发电与抽水出力约束 xjf≤ ≤ gxj(9) ,,f Yj,f≤P P≤ p (10) ,f 两式中, 、 s分别为抽水蓄能机组 的最 小和最大发电出力; 、 舣p是最小和最 大抽水出力; ,、 ,,为两个布尔变量,表示机 组在t时段的运行状态,其值为0表示停机,值为1 表示开机。 考虑到抽水蓄能电站运行的经济性,要求机组 不可同时发生发电和抽水的情况,即 +Y,f≤1 (11) ,(4)起停次数约束 在工作中,过多的起停机组对抽水蓄能机组会 造成损害,且起停时会增加水头的损失,从而增加 成本,因此需要对每日的机组起停次数进行。 设 表示每台机组每日最大起停次数,则机组每 曰起停次数 应不大于 即 ≤ (12) 3)系统备用容量约束 考虑到抽水蓄能机组可以实现发电/抽水间的 快速转换,因而抽水状态下可以提供的备用容量为 抽水功率和最大发电功率的和。 故当抽蓄机组处于发电工况时,有 2018 ̄2 嘲|_I技玳I 23 研究与开发 ∑( m 一 , )+∑( g一 )≥Rt(13) i=1 j=l 当抽蓄机组处于抽水工况时,有 ∑( …一 ,,)+∑( + )≥Rt (14) i=l j=l 式中,R,为旋转备用率,按系统总负荷的7%[ 来 考虑。 2 求解方法 建立的数学模型属于MIP(混合整数规划模 型),通常研究人员都使用Matlab自带的求解器或 者Lingo软件来解决规模较小的整数规划模型,但 是此模型的变量和约束条件相对较多,规模较大, 用传统的方法难以计算。本文采用在Matlab中使用 Yalmip工具箱调用Gurobi优化软件求解。 Yalmip是由Lofbeig开发的免费Matlab工具箱, 其最大的优点在于可以调用外部的优化求解器,并 且有自己的一套简洁直观的编程语言,可以一定程 度的减少程序编译的难度。Gurobi是由美国Gurobi 公司开发的新一代大规模数学规划优化器,在 Decision Tree for Optimization Software网站举行的 第三方优化器评估中,展示出更快的优化速度和精 度,是目前世界上顶尖的求解线性规划、整数规划 和某些非线性规划的软件包之一,在很多行业都有 广泛应用。 3 算例分析 3.1参数设置 以某地区电网作为研究对象,系统中有4台火 电机组,参数见表1;抽水蓄能电站的参数见表2, 其起动和停机的收费标准来自瑞典学者的调研报 告,抽蓄电站每次起动收费约2000元fI 7J;典型日一 天24h每个时段负荷数据见表3。 表1火电机组参数 参数 PGl PG2 PG3 PG4 P …/MW 570 425 700 6l0 Pi,mi /MW 228 l70 280 244 日f 0.000176 0.000230 0.000117 0.000l51 b 0.11O6 O.1509 0.0704 0.1308 C 3.02 5.O3 7.04 3.Ol 爬坡约束/MW l7l 128 21O 183 24 l电囊一I技,代 2o18年第2期 表2抽水蓄能机组参数 参数 取值 水库蓄能最大/最小值/MW·h 900/0 水泵抽水最大/最小值/MW l20/0 水力发电最大/最小值/MW 120/0 抽水效率、发电效率 0.9 每日最大起停次数 2 表3典型日负荷数据 时刻 负荷/MW 时刻 负荷/MW 时刻 负荷/MW 1:O0 1179 9:O0 l772 17:00 1653 2:OO 1149 1O:OO 1831 18:00 1980 3:0O 1060 11:O0 1861 19:00 2109 4:O0 l030 12:00 2068 20:00 1772 5:00 l179 13:0O l920 2l:O0 1327 6:00 1327 l4:O0 1772 22:00 1149 7:O0 l624 15:00 l624 23:O0 1090 8:00 1713 l6:O0 1476 24:00 1031 3.2 结果分析 为验证抽蓄对系统的调节作用,设置两种方案: ①无抽蓄的优化运行;②抽蓄电站加入后的最优运 行。通过Matlab编译程序调用Gurobi软件,计算 得到抽水蓄能电站优化运行结果(见图1)。图1为 抽蓄机组的优化出力曲线图,大于零为发电工况, 小于零为抽水工况(为了便于对比,将抽蓄出力计 算结果放大l0倍)。可以看到抽水蓄能机组发挥了其 调峰填谷的作用,在负荷较低的第0:00—5:o0以及 第21:00—24:00时段进行抽水,在10:00—19:00 点负荷较高时进行发电,并且在负荷最高的第12:00、 13:00、18:00和19:00时段出力明显最高。 时刻 图1抽水蓄能机组出力 表4是在两种运行方案下,系统的总煤耗,总 成本以及常规机组发电量的对比。 研究与开发 表4 两种运行方案计算结果 运行方案 总煤耗/t 总成本/元 常规机组发电量/(MW·h) A B 6686.0 6649.7 1203480 l200954 36696 36836 可以得到,在抽蓄电站投入的情况下,因调峰 填谷的功能,可以每天减少煤耗量0.6%(36.3t); 火电机组总发电量增大,这是抽水蓄能机组在发挥 调峰填谷作用时,因为其能量转换时产生的损耗, 导致抽水工况时吸收的火电机组电量小于火电机组 发出的电量,并且在发电工况时进一步产生损耗, 最终导致了火电机组总发电量增大,具体各火电机 组发电情况见表5。 表5火电机组出力比较 机组 方案A 方案B 1 8407 8340 2 兰 得怒 5069 4867 O O O O O O O O 0 O ∞勰 加 3 l4683 l5388 4 8537 8241 火电机组发电量增大但煤耗减少的原因分析 如下: 1)4台火电机组的煤耗系数各不相同,可以看 到只有煤耗量系数最小的机组3的日出力增大,其 余3台机组的出力均有不同程度的减少。此消彼长 之下使得系统的总煤耗减少。 2)4台机组在两种方案的出力对比如图2所示, 图中虚线为方案A的数据,实线为方案B的数据。 抽水蓄能机组调峰填谷的作用使得其在负荷低谷时 段吸收煤耗系数较小的机组1和机组3的电量;在 负荷高峰时段发电,减少火电机组的发电压力,令 除去煤耗系数最小的机组3之外,其他3台机组的 出力均减小,从而降低火电机组的运行成本。 (a)机组1 羹 怒 (d)机组4 图2火电机组出力 4结论 本文以某地区电网的实际数据为基础,以系统 运行成本最低为目标建立了抽水蓄能电站的运行优 化模型。算例分析表明,抽水蓄能电站可以很好的 跟踪负荷的变化,提高了系统的经济性。同时,实 验验证了Gurobi优化软件在电力系统运行优化规划 方面的实用性,结果表明Gurobi计算效率高,速度 快,且在Yalmip工具箱的作用下,程序编译更加便 捷,为电力系统运行规划的计算方法提供了新的 思路。 参考文献 [1]Wang H,Jun L I,Jiang H,et a1.Design ofthe Control 2018年第2期霸 麓柬I 25 研究与开发 Sequence for Monitoring and Control System of 水电自动化与大坝监测,2005,29(6):6-10. Pushihe Pumped Storage Power Station[J].Hydropower Automation&Dam Monitoring,2010,1:329-33 1. [11]唐海华,黄春雷,丁杰.混合式抽水蓄能电站优化 调度策略[J].电力系统自动化,2011,35(21):40—45. [12】杨柳青,林舜江,刘明波,等.考虑风电接入的大型 电力系统多目标动态优化调度[J].电工技术学报, 2014,29(10):286—295. [2]何永秀,关雷,蔡琪,等.抽水蓄能电站在电网中的 保安功能与效益分析[J]_电网技术,2004,28(20): 54—57.67. [3] 肖白,丛晶,高晓峰,等.风电一抽水蓄能联合系统 综合效益评价方法[J].电网技术,2014,38(2):400— 404. [13]王伟华,王红宇,许国瑞,等.基于时步有限元的抽 水蓄能电机瞬态参数计算方法的对比[J].电工技术 学报,2015,30(11:89.97. [4】程诗明,郭明宇,贺儒飞.基于实测信号的抽水蓄 能机组频率自适应算法研究[J].电气技术,2017, 18(4):106—109. [14】邹金,赖旭,汪宁渤.以减少电网弃风为目标的风 电与抽水蓄能协调运行[J].电网技术,2015,39(9): 2472.2477. [5] 冯凌云,邱小波,吴吴.抽水蓄能机组调速系统动 态仿真平台的设计与开发 .电气技术,2015, 16(12):98—101. [15】胡泽春,丁华杰,孔涛.风电一抽水蓄能联合日运行 优化调度模型[J].电力系统自动化,2012,36(2): 36—41.57. [6】 徐得潜,翟国寿.抽水蓄能电站与火电站配合运行 优化模型研究[J].水力发电学报,1996(4):l1.20. [16]徐飞,陈磊,金和平,等.抽水蓄能电站与风电的联 合优化运行建模及应用分析[J].电力系统自动化, 2013,37(1):149-154. 【7】 丘文干.抽水蓄能电站及其系统运行优化的概率模 拟算法[J].电网技术,2006,30(15):95.100. [8】 柴大鹏,成欢,葛宇剑,等.火电一抽水蓄能联合运 营系统的效益分析及优化算法研究[J].华东电力, 2014,42(5):1012-1019. [9]丘文千.抽水蓄能电站日调节与周调节方式运行优 [1 7]Nilsson O,sjelvgren D.Hydro unit start—up costs and their impact on the short term scheduling strategies of Swedish power producers[J].IEEE Transactions on Power Systems,1997,12(1):38—44. 化模型[J].水电自动化与大坝监测,2004,28(5): 51.54. 收稿日期:2017-07.12 作者简介 张津玮(1993一),男,硕士研究生,研究方向为电力系统运行分析。 【10]丘文干.抽水蓄能电站运行优化的动态规划模型[J]. (上接第21页) [8]刘杰.表贴式永磁同步电机无位置传感器驱动技术 [4】汤长波,顾敏明,潘海鹏.无位置传感器无刷直流 电机启动方法研究[J].浙江理工大学学报,2011, 28(6):886—889. 研究[D].哈尔滨:哈尔滨工业大学,2014. [9] 王子辉.永磁同步电机全速度范围无位置传感器控 制策略研究[D].杭州:浙江大学,2012. [10】王力,黄科元,黄守道,等.改进的永磁同步电机无 传感器快速启动方法[J].电力电子技术,2015,49(7) 62—64. [5] 万晓凤,康利平,郑博嘉.无位置传感器BLDCM的 软件升频启动策略[J].工程设计学报,2014,21(6): 596—602. 【6】刘颖,周波,冯瑛,等.基于脉振高频电流注入 SPMSM低速无位置传感器控制[J].电工技术学报, 2012,27(7):139—145. 收稿日期:2017-08.22 作者简介 林卫鹏(1991一),男,福建省漳州市龙海市人,硕士研究生,研究 方向为新型电机理论与控制技术。 [7] 李翠萍,王新生,高阳.无位置传感器永磁同步电 动机控制系统[J].微特电机,2007,35(12):39.41,45. 26 I瞎一|l藏硪2018年第2期