留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于仿射鲁棒优化的多站融合协同优化调度

张笑演 郭创新 金国胜 尹康 高亚栋 周盈

张笑演, 郭创新, 金国胜, 等. 基于仿射鲁棒优化的多站融合协同优化调度[J]. 现代电力, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
引用本文: 张笑演, 郭创新, 金国胜, 等. 基于仿射鲁棒优化的多站融合协同优化调度[J]. 现代电力, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
ZHANG Xiaoyan, GUO Chuangxin, JIN Guosheng, etc. Collaborative Optimal Dispatch of Multi-Station Integration Based on Affine Robust Optimization[J]. Modern Electric Power, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
Citation: ZHANG Xiaoyan, GUO Chuangxin, JIN Guosheng, etc. Collaborative Optimal Dispatch of Multi-Station Integration Based on Affine Robust Optimization[J]. Modern Electric Power, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139

基于仿射鲁棒优化的多站融合协同优化调度

doi: 10.19725/j.cnki.1007-2322.2021.0139
详细信息
    作者简介:

    张笑演(1998),男,硕士研究生,研究方向为综合能源系统经济调度与优化运行,E-mail:673107293@qq.com

    郭创新(1969),男,博士,教授,通信作者,研究方向为能源互联网运行与规划、风险调度、新能源接入,E-mail:guochuangxin@zju.edu.cn

  • 中图分类号: TM715

Collaborative Optimal Dispatch of Multi-Station Integration Based on Affine Robust Optimization

  • 摘要: 为解决新能源出力不确定性的多能电站(数据中心、电动汽车充电站、分布式光伏、风电、储能站)优化运行和收益分配问题,提出基于仿射鲁棒优化的多站融合协同优化策略。考虑数据中心与其余子系统分属不同运营商,对多站融合系统中各子系统进行建模,其中数据中心的负载具有可平移负荷的特性;通过盒式不确定集对风电和光伏出力的不确定性进行建模;以多站融合协同总收益最大化为目标函数,采用仿射鲁棒优化方法进行求解。提出了基于影子价格的能源站与数据中心站的收益分配方法。最后,通过仿真算例说明了所提模型的有效性和所提收益分配方法的合理性。
  • 图  1  多站融合系统拓扑

    Figure  1.  Topology of multi-station integration

    图  2  数据中心负荷预测

    Figure  2.  Load prediction of data center

    图  3  典型日光伏、风机预测出力曲线

    Figure  3.  Forecasted output of photovoltaic generation and wind generator in a typical day

    图  4  购电价格

    Figure  4.  Purchase/sell electricity price

    图  5  各个时刻电动汽车并网数量

    Figure  5.  The number of electric vehicles connected to the grid at each moment

    图  6  能源站电、冷优化调度结果

    Figure  6.  Optimal dispatch result of power, heating and cooling load in energy station

    图  7  电动汽车充放电情况

    Figure  7.  Charging and discharging conditions of electric vehicles

    图  8  设备参与因子优化结果

    Figure  8.  Optimization results of equipment participating factors

    图  9  数据中心电负荷调度前后对比

    Figure  9.  Comparison of power load of data center before and after dispatching

    图  10  数据中心冷负荷调度前后对比

    Figure  10.  Comparison of cooling load of data center before and after dispatching

    图  11  2种场景下与电网交换功率情况

    Figure  11.  Power exchange with the grid in two scenarios

    A1  设备参数

    A1.   Parameters of equipment

    设备名称燃气轮机燃气锅炉吸收式制冷机电制冷机电储能光伏风机电动汽车充电桩
    投资费用/(元/kW)282662812569501800753680001000
    维修费用/(元/kW·h)0.0940.0030.0140.0080.0060.0120.015
    发电效率0.43
    热电比1.03
    制热效率0.92
    制冷效率1.13.5
    充电效率0.950.9
    放电效率0.950.9
    容量150050010001000150040040060
    soc下限/并网soc0.47≤0.4
    soc上限/离网soc0.951
    最大功率/kW3003030
    下载: 导出CSV

    表  1  “多站融合”成本分析

    Table  1.   Cost analysis on multi-station integration

    模型
    编号
    购电
    成本
    /(元/日)
    购气
    成本
    /(元/日)
    运维
    成本
    /(元/日)
    充电站
    收益
    /(元/日)
    数据中
    心收益
    /(元/日)
    净收益
    /(元/日)
    1−30120153263805497
    210593463213202703202936451
    310307461813072712201846664
    下载: 导出CSV

    表  2  各个设备的可用度

    Table  2.   Availability of each device

    设备可用度
    燃气轮机 0.991293
    吸收式制冷机0.995732
    电制冷机0.996898
    电储能0.999830
    光伏0.9999
    风机0.9999
    断路器0.999998
    下载: 导出CSV
  • [1] 韩晓平. 关于国网革命的几点看法[J]. 上海节能, 2019(4): 241−244.

    HAN Xiaoping. Some views on the revolution of state grid[J]. Shanghai Energy Conservation, 2019(4): 241−244(in Chinese).
    [2] 刘成运, 孟超, 景锐, 等. A级数据中心综合能源系统多目标优化设计和调度[J]. 电力系统自动化, 2019, 43(14): 136−142. doi:  10.7500/AEPS20180731007

    LIU Chengyun, MENG Chao, JING Rui, et al. Multi-objective optimization design and dispatch of A level data center through integrated energy system[J]. Automation of Electric Power Systems, 2019, 43(14): 136−142(in Chinese). doi:  10.7500/AEPS20180731007
    [3] 杜先波, 陈国琳, 唐一铭, 等. 多站融合型数据中心交直流供电模式对比分析[J]. 电力需求侧管理, 2021, 23(1): 84−89. doi:  10.3969/j.issn.1009-1831.2021.01.016

    DU Xianbo, CHEN Guolin, TANG Yiming, et al. Comparison and analysis of AC/DC power supply mode of data center under multi-station integration[J]. Power Demand Side Management, 2021, 23(1): 84−89(in Chinese). doi:  10.3969/j.issn.1009-1831.2021.01.016
    [4] 白中华, 李强, 陈晶, 等. 储能电站在多站融合场景下的运行策略优化[J]. 中国电力, 2021, 54(6): 136−144.

    BAI Zhonghua, LI Qiang, CHEN Jing, et al. Operation strategy optimization of energy storage power station under multi-station integration[J]. Electric Power, 2021, 54(6): 136−144(in Chinese).
    [5] 王伯伊, 张越, 刘明波, 等. “多站融合”运营模式的探索研究[J]. 电力信息与通信技术, 2019, 17(7): 41−45.

    WANG Boyi, ZHANG Yue, LIU Mingbo, et al. Exploration and research on the operation mode of Multi-station Integration[J]. Electric Power Information and Communication Technology, 2019, 17(7): 41−45(in Chinese).
    [6] 杨彭, 顾颖, 周旭. 多站融合业务应用场景与建设运营模式探索研究[J]. 数据通信, 2020(3): 1−3.

    YANG Peng, GU Ying, ZHOU Xu. Research on business application scenarios and construction operation mode of multi-station integration[J]. Shuju Tongxin, 2020(3): 1−3(in Chinese).
    [7] 中华人民共和国工业和信息化部. 数据中心设计规范GB 50174—2017[S]. 北京: 中国计划出版社, 2017: 20
    [8] 郝建华, 索龙, 李红艳. 数据中心网络基于路径关键度的拥塞避免重路由方法[J]. 电力信息与通信技术, 2021, 19(5): 6−13.

    HAO Jianhua, SUO Long, LI Hongyan. Path criticality based congestion avoidance rerouting method for data center networks[J]. Electric Power Information and Communication Technology, 2021, 19(5): 6−13(in Chinese).
    [9] 许刚, 张丙旭, 张广超. 电动汽车集群并网的分布式鲁棒优化调度模型[J]. 电工技术学报, 2021, 36(3): 565−578.

    XU Gang, ZHANG Bingxu, ZHANG Guangchao. Distributed robust optimal dispatch model for electric vehicles connected to grid[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 565−578(in Chinese).
    [10] 施云辉, 郭创新, 丁筱. 基于仿射可调鲁棒优化的园区综合能源系统经济调度[J]. 发电技术, 2020, 41(2): 118−125. doi:  10.12096/j.2096-4528.pgt.19139

    SHI Yunhui, GUO Chuangxin, DING Xiao. Integrated energy system economic dispatch based on affine adjustable robust optimization[J]. Power Generation Technology, 2020, 41(2): 118−125(in Chinese). doi:  10.12096/j.2096-4528.pgt.19139
    [11] SHIMIZU K, AIYOSHI E. Necessary conditions for Min-Max problems and algorithms by a relaxation procedure.[J]. IEEE Transactions on Automatic Control, 1980, 25(1): 62−67. doi:  10.1109/TAC.1980.1102226
    [12] 李利利, 丁恰, 涂孟夫, 等. 机组组合问题的仿射可调整鲁棒优化模型与算法[J]. 电力工程技术, 2017, 36(3): 31−35.

    LI Lili, DING Qia, TU Mengfu, et al. Affinely adjustable robust optimization model and algorithm for unit commitment problem[J]. Electric Power Engineering Technology, 2017, 36(3): 31−35(in Chinese).
  • [1] 潘华, 姚正, 林顺富, 方静, 颜静汝, 高旭.  基于信息间隙决策理论的含光热电站及热泵的综合能源系统低碳调度优化 . 现代电力, 2022, 39(2): 169-181. doi: 10.19725/j.cnki.1007-2322.2021.0058
    [2] 李彬, 杜亚彬, 曹望璋, 祁兵, 孙毅, 陈宋宋.  考虑风光储互补与工作负载分配的数据中心优化调度 . 现代电力, 2022, 39(3): 356-362. doi: 10.19725/j.cnki.1007-2322.2021.0129
    [3] 仲志强, 成海生, 蔡华, 张辉辉, 江御龙.  基于分时电价和需求响应的家庭微电网系统协同控制策略及其实现 . 现代电力, 2021, 38(1): 69-78. doi: 10.19725/j.cnki.1007-2322.2020.0101
    [4] 魏震波, 都成.  计及天气对电动汽车接入家庭型微电网供需侧影响的经济调度模型 . 现代电力, 2020, 37(1): 10-18. doi: 10.19725/j.cnki.1007-2322.2019.0031
    [5] 蒋怡静, 于艾清, 屠亚南.  基于IBQPSO算法的电动汽车时空双尺度有序充电引导策略 . 现代电力, 2019, 36(6): 1-8.
    [6] 林俐, 张玉, 王世谦.  计及微电源效率约束的多源微网优化调度 . 现代电力, 2019, 36(4): 17-23.
    [7] 何阳, 张宇, 王育飞, 金商鹤.  考虑负荷优化的电动汽车光伏充电站储能容量配置 . 现代电力, 2019, 36(5): 76-81.
    [8] 赵琦玮, 王昕, 王鑫, 郎永波, 贾立凯.  微电网环境下考虑日前预测误差的电动汽车多时间尺度优化调度模型 . 现代电力, 2019, 36(5): 47-53.
    [9] 陈小闽, 王钢, 汪隆君.  基于车联网框架的电动汽车有序充电策略 . 现代电力, 2018, 35(4): 1-7.
    [10] 张禄, 孙舟, 王伟贤, 李香龙, 周杨, 黄其进, 陈雁.  应用曲线分群预测的电动汽车充电设施规划方法 . 现代电力, 2018, 35(4): 21-26.
    [11] 苏小林, 张艳娟, 武中, 阎晓霞.  规模化电动汽车充电负荷的预测及其对电网的影响 . 现代电力, 2018, 35(1): 45-54.
    [12] 满洲, 郭春林, 裴鑫, 胡德权.  基于V2G的电动汽车无功功率补偿技术 . 现代电力, 2018, 35(4): 15-20.
    [13] 刘广, 曾成碧, 苗虹.  基于粗糙集的电动汽车充电站规划综合评价 . 现代电力, 2018, 35(1): 55-60.
    [14] 周凌锋, 王杰.  基于时空分布负荷预测的电动汽车充电优化 . 现代电力, 2018, 35(5): 10-16.
    [15] 胡文平, 巫伟南, 郝婉梦, 于腾凯, 杨 军.  含电动汽车的配电网运行风险评估 . 现代电力, 2016, 33(2): 45-50.
    [16] 张乐平, 张明明, 林伟斌.  电动汽车与电网统一互动架构设计与探讨 . 现代电力, 2014, 31(1): 34-39.
    [17] 刘人玮, 尹忠东, 曹松伟, 姜喆
    基于直接功率控制的电流型PWM电动汽车充电装置 . 现代电力, 2014, 31(2): 27-30.
    [18] 许璐, 袁越, 王敏, 刘冠群.  含电动汽车充电站的配电网谐波潮流计算 . 现代电力, 2014, 31(2): 21-26.
    [19] 寇凌峰, 刘自发, 周 欢.  区域电动汽车充电站规划的模型与算法 . 现代电力, 2010, 27(4): 44-48.
    [20] 杨 婷, 景占荣, 高 田.  电动汽车车载智能快速充电器的研究与设计 . 现代电力, 2010, 27(5): 62-66.
  • 加载中
图(11) / 表 (3)
计量
  • 文章访问数:  182
  • HTML全文浏览量:  52
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-05-03
  • 刊出日期:  2022-07-27

基于仿射鲁棒优化的多站融合协同优化调度

doi: 10.19725/j.cnki.1007-2322.2021.0139
    作者简介:

    张笑演(1998),男,硕士研究生,研究方向为综合能源系统经济调度与优化运行,E-mail:673107293@qq.com

    郭创新(1969),男,博士,教授,通信作者,研究方向为能源互联网运行与规划、风险调度、新能源接入,E-mail:guochuangxin@zju.edu.cn

  • 中图分类号: TM715

摘要: 为解决新能源出力不确定性的多能电站(数据中心、电动汽车充电站、分布式光伏、风电、储能站)优化运行和收益分配问题,提出基于仿射鲁棒优化的多站融合协同优化策略。考虑数据中心与其余子系统分属不同运营商,对多站融合系统中各子系统进行建模,其中数据中心的负载具有可平移负荷的特性;通过盒式不确定集对风电和光伏出力的不确定性进行建模;以多站融合协同总收益最大化为目标函数,采用仿射鲁棒优化方法进行求解。提出了基于影子价格的能源站与数据中心站的收益分配方法。最后,通过仿真算例说明了所提模型的有效性和所提收益分配方法的合理性。

English Abstract

张笑演, 郭创新, 金国胜, 等. 基于仿射鲁棒优化的多站融合协同优化调度[J]. 现代电力, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
引用本文: 张笑演, 郭创新, 金国胜, 等. 基于仿射鲁棒优化的多站融合协同优化调度[J]. 现代电力, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
ZHANG Xiaoyan, GUO Chuangxin, JIN Guosheng, etc. Collaborative Optimal Dispatch of Multi-Station Integration Based on Affine Robust Optimization[J]. Modern Electric Power, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
Citation: ZHANG Xiaoyan, GUO Chuangxin, JIN Guosheng, etc. Collaborative Optimal Dispatch of Multi-Station Integration Based on Affine Robust Optimization[J]. Modern Electric Power, 2022, 39(4): 379-387. doi: 10.19725/j.cnki.1007-2322.2021.0139
    • 国家电网公司于2019年提出打造枢纽型、平台型、共享型企业,建设运营好坚强智能电网和泛在电力物联网(“三型两网”)[1]。并且随着中国科技的发展,海量的数据需求和边缘计算技术需求也呈现出高速增长的态势。基于此背景,人们提出了多站融合的概念,即挖掘利用变电站资源,使之融合光伏、风机等分布式发电站,边缘数据中心站,储能站以及充电站等。对于数据中心,文献[2]研究了在综合能源系统下的优化运行方法,考虑了不同季节下的用能需求,具有一定参考价值;文献[3]研究了交直流供电方式对于数据中心供电可靠性的影响。对于储能站,文献[4]研究了在多站融合场景下储能站的优化运行方法,采用储能站年净收益的年值为目标函数,对其容量进行优化配置,并以节点收益减去电池损耗为目标函数进行运行策略的优化,但是对于充电站的建模只是采用了电动汽车典型负荷曲线,没有考虑并网时刻的随机性和充电桩的数量对于充电站负荷的约束。对于多站融合的运营模式,文献[5]分别阐述了数据中心站、充电站和储能站的运营模式,并分为面向电网、面向企业2个方向进行分析,最后结合实际案例对运营模式可行性进行探究,但是只罗列了最后的收益情况,缺乏中间计算过程;文献[6]针对多站融合业务特点,创造性地提出了“电网独资”、“建设–移交”、“建设–经营–移交”和“委托运营”等建设运营模式。以上对于多站融合的研究中,鲜有涉及对新能源出力不确定性的考虑。

      本文在上述研究的基础上,在多站融合场景下,考虑数据中心的负载存在可平移的特性、3类不同电动汽车的充电需求以及电动汽车并网时间的随机性、风机和光伏出力的不确定性,通过盒式不确定集进行建模,并且考虑充电站中充电桩的数量对其负荷的约束情况,以多站融合协同总收益最大化为目标函数,采用仿射鲁棒优化方法进行求解。提出基于影子价格的能源站与数据中心站的收益分配方法。最后,通过仿真算例说明所提模型的有效性和所提收益分配方法的合理性。

    • 典型的多站融合系统拓扑如图1所示。风机、光伏和燃气轮机构成了分布式发电站,为了给数据中心提供相应的冷负荷,又加入了电制冷机、吸收式制冷机和燃气锅炉。本文考虑的数据中心为B级数据中心,为保证数据中心供电可靠性,接了两路电网线路[7],其中一路在正常运行时不向数据中心供电,只有在发生事故时才会进行供电。

      图  1  多站融合系统拓扑

      Figure 1.  Topology of multi-station integration

    • 数据中心通常利用虚拟机技术实现服务器集群的大规模配置,并将实体服务器以虚拟机的形式租赁给用户使用[8]。工作负载即定义为由虚拟机在完整的软件环境下处理的、连续且不可中断的批处理计算作业请求[8],每一个工作负载都包含一定的CPU利用率、调度等级、允许执行时限等参数。本文设定数据中心的工作负载分为实时性负载W0和3种非实时性的负载W1W2W3,其中W0需要在任务到达后1 h内完成,CPU占用率为g0W1W2W3则可在任务到达后2~4 h内完成,计算资源占用率分别为g1g2g3

      在常规模式下,所有计算任务将被平均分配到所允许的各个时段。然而,运营商可以提前非实时性负载的完成时间来提升用户满意度,也可以灵活调整各时段的任务分配来实现电负荷的削峰填谷。以调度因子$x_{i,s}^t$为决策变量,代表$t$时段抵达数据中心的任务${W_i}$$t + s$时段执行的比例。因此,$t$时段调度后数据中心的总任务量可以表示为

      $$ {W^t} = W_0^t + \displaystyle\sum\limits_{i = {\text{1}}}^{\text{3}} {\displaystyle\sum\limits_{s = {\text{0}}}^i {x_{i,s}^{t - s}W_i^{t - s}} } $$ (1)

      $t$时段调度后各非实时性任务的任务量分别为

      $$ W_{{\text{d}}i}^t = \displaystyle\sum\limits_{s = {\text{0}}}^i {x_{i,s}^{t - s}} W_i^{t - s}\quad i = {\text{1,2,3}} $$ (2)

      $t$时段调度后计算中心的CPU占用率为

      $$ {g^t} = {g_{\text{0}}}W_{\text{0}}^t + \displaystyle\sum\limits_{i = {\text{1}}}^{\text{3}} {{g_i}} W_{{\text{d}}i}^t $$ (3)

      计算中心电负荷为服务器集群的功耗,等于CPU实时功耗和静态功耗的总和,因此数据中心电负荷$E_{{\text{DC}}}^t$的表达式为

      $$ E_{{\text{DC}}}^t = k{g_t} + bN $$ (4)

      式中:k为服务器集群CPU功耗系数;b为静态功耗系数;N为开机总数。由于本文假定各服务器同构,因此开机总数又通过公式(5)约束

      $$ {\text{0}}{\text{.9}}N \geqslant {g_t}\quad t = {\text{1,}}\cdots{\text{,}}T $$ (5)

      式中:0.9为冗余系数,代表单台服务器CPU占用率不超过90%,从而避免拥塞的发生。

      数据中心运行过程中会产生大量的热量,需要制冷系统将机房的温度维持在规定的范围内,因此,引入数据中心冷负荷$C_{{\text{DC}}}^t$,其表达式为

      $$ C_{{\text{DC}}}^t = E_{{\text{DC}}}^t/{F_{\text{C}}} $$ (6)

      式中:$ {F_{\text{C}}} $为负荷系数,与服务器电源效率、CPU效率等因素相关。

    • 储能站可以平抑光伏和风机出力的波动,提升新能源的消纳水平。此外,储能站还具有削峰填谷的作用,在低谷时进行充电,高峰时进行放电。储能站的充放电模型如式(7)所示

      $$\begin{cases} \begin{array}{l} {S}_{\text{cha}}+{S}_{\text{dis}}\leqslant\text{1}\\ -{S}_{\text{cha}}\cdot {P}_{\text{r}}\leqslant{P}_{\text{cha,}t}\leqslant\text{0}\\ \text{0}\leqslant{P}_{\text{dis},t}\leqslant{S}_{\text{dis}}\cdot {P}_{\text{r}}\end{array} \end{cases}$$ (7)

      式中: $ {S_{{\text{cha}}}} $$ {S_{{\text{dis}}}} $为二进制变量,表征储能站的充放电状态,避免储能站既充电又放电;$ {P_{\text{r}}} $为储能电站的额定充放电功率;$ {P_{{\text{cha,}}t}} $$ {P_{{\text{dis}},t}} $分别表示t时刻的充电、放电功率。储能站的荷电状态${S_{{\text{oc}}}}$(state of charge)模型如式(8)所示

      $$ \begin{cases}\begin{array}{l}{S}_{\text{oc}}\text{(}t+\text{1)}={S}_{\text{oc}}\text{(}t\text{)}-\left(\text{1/}{E}_{\text{es}}\right)\cdot \left({P}_{\text{cha ,}t}{\eta }_{\text{s}}\text+{P}_{\text{dis,}t}\text{/}{\eta }_{\text{s}}\right)\cdot \Delta t\\ {S}_{\text{ocmin}}\leqslant{S}_{\text{oc}}\text{(}t\text{)}\leqslant{S}_{\text{oc}}{}_{\text{max}},\forall t\end{array}\end{cases} $$ (8)

      式中:$ {E_{{\text{es}}}} $为储能站的额定容量;$ \Delta t $为调度时间间隔;$ {S_{{\text{oc}}}}{\text{(}}t{\text{)}} $为储能站在$t$时刻的荷电状态;$ {S_{{\text{ocmin}}}} $$ {S_{{\text{oc}}}}_{{\text{max}}} $分别为储能站荷电状态的下限和上限;$ {\eta _{\text{s}}} $为储能站的充放电效率。因为储能站要保证数据中心站的供电可靠性,在任何时候都要保证数据中心站最大负荷运行30min的备用容量[7],因此在多站融合场景下考虑调度问题时储能电站允许的SOC下限值会比在一般情况下高。

    • 对于充电站的建模主要在于对充电汽车的建模。本文仿照文献[9]的做法,将电动汽车按需求分为3类。第1类电动汽车一直以额定功率充电,对应实际中希望以最快速度充满电的需求情况,对应的充电费用较高。第2类电动汽车允许充电功率小于额定功率,但不能放电,并且必须在4h内充满电,充电费用较低,对应实际中用车不迫切,希望能减少充电费用的需求。第3类电动车允许充电功率低于额定功率,并且允许放电,但必须在6h内完成充电,充电费用最低。对应实际中有充裕的时间可以进行充电,希望最小化充电成本的需求。第1类电动汽车的充电模型如式(9)所示,对于$\forall i \in{C_1} $

      $$ \begin{cases}\begin{array}{l}{P}_{i,t}={P}_{\text{cr}}\\ {S}_{\text{oc}\text{,}}{}_{i}\text{(}t+\text{1)}={S}_{\text{oc}\text{,}}{}_{i}\text{(}t\text{)}+{P}_{i,t}\cdot {\eta }_{\text{c}}\cdot \Delta t/{E}_{\text{c}}\end{array}\end{cases} $$ (9)

      式中:${C_{\text{1}}}$表示第1类电动汽车集合;$ {P_{i,t}} $表示电动汽车在$ t $时刻的充电功率;$ {P_{{\text{cr}}}} $表示充电桩额定充电功率;$ {\eta _{\text{c}}} $表示电动汽车充放电效率;$ {E_{\text{c}}} $为电动汽车额定容量。

      第2类电动汽车的充电模型如式(10)所示,对于$\forall i \in {C_{\text{2}}}$

      $$ \begin{cases} \begin{array}{l}\text{0}\leqslant{P}_{i,t}\leqslant{P}_{\text{cr}}\\ {S}_{\text{oc}\text{,}}{}_{i}\text{(}t+\text{1)}={S}_{\text{oc}\text{,}}{}_{i}\text{(}t\text{)}+{P}_{i,t}\cdot {\eta }_{\text{c}}\cdot \Delta t/{E}_{\text{c}}\\ {S}_{\text{oc}\text{,}}{}_{i}\left({t}_{i,\text{0}}+\text{4}\right)=\text{1}\end{array}\end{cases} $$ (10)

      式中:${C_{\text{2}}}$表示第2类电动汽车的集合;$ {t_{i,{\text{0}}}} $表示电动汽车$i$的并网时刻。

      第3类电动汽车的充放电模型如式(11)所示,对于$\forall i \in {C_{\text{3}}}$

      $$ \begin{cases}\begin{array}{l}-{P}_{\text{cr}}\leqslant{P}_{i,t}\leqslant{P}_{\text{cr}}\\ {S}_{\text{oc}\text{,}}{}_{i}\text{(}t+\text{1)}={S}_{\text{oc}\text{,}}{}_{i}\text{(}t\text{)}+{P}_{i,t}\cdot {\eta }_{\text{c}}\cdot \Delta t/{E}_{\text{c}}\\ {S}_{\text{oc}\text{,}}{}_{i}\left({t}_{i,\text{0}}+\text{6}\right)=\text{1}\end{array} \end{cases}$$ (11)
    • 本文假定能源站与数据中心的利益不存在明显冲突,两者采取合作盈利的模式。能源站的优化目标为提高充电站的收益,并减少总购能成本和运维费用,数据中心站的目标为减少用电和用冷的费用。两者首先进行协同优化,接着用能站按照对偶变量所得的边际电价和边际冷价向数据中心征收费用。协同优化的目标函数为

      $$ {\text{max}}\left( {{A_{{\text{car}}}} - {C_{{\text{gas}}}} - {C_{\text{e}}} - {C_{\text{m}}}} \right) $$ (12)

      式中:$ {A_{{\text{car}}}} $表示给电动汽车充电获得的收益,其计算公式为

      $$ \begin{split} &\sum\limits_{t = {\text{1}}}^T {\sum\limits_{i = {\text{1}}}^{{n_{\text{1}}}} {{\rho _{{\text{c,1}}}}} } \cdot {E_{\text{c}}} \cdot \left( {{\text{1-}}{S_{{\text{oc}}}}\left( {{t_{i,{\text{0}}}}} \right)} \right) + \\ {A_{\rm{{car}}}} = &\sum\limits_{t = {\text{1}}}^T {\sum\limits_{i = {\text{1}}}^{{n_{\text{2}}}} {{\rho _{c,{\text{2}}}}} } \cdot {E_{\text{c}}} \cdot \left( {{\text{1-}}{S_{{\text{oc}}}}\left( {{t_{i,{\text{0}}}}} \right)} \right) + \\ &\sum\limits_{t = {\text{1}}}^T {\sum\limits_{i = {\text{1}}}^{{n_{\text{3}}}} {{\rho _{c,{\text{3}}}}} } \cdot {E_{\text{c}}} \cdot \left( {{\text{1-}}{S_{{\text{oc}}}}\left( {{t_{i,{\text{0}}}}} \right)} \right) \\ \end{split} $$ (13)

      式中:$ {n}_{\text{1}}、{n}_{\text{2}}、{n}_{\text{3}} $分别表示第1类电动汽车、第2类电动汽车、第3类电动汽车的数量;$\,{\rho }_{\text{c,1、}} {\rho }_{\text{c,2、}} \,{\rho }_{\text{c,3}}$分别表示第1类电动汽车、第2类电动汽车、第3类电动汽车的单位充电价格。

      ${C_{{\text{gas}}}}$为购买天然气的费用,其计算公式为

      $$ {C_{{\text{gas}}}} = \sum\limits_{t = 1}^T {{\rho _{{\text{gas}}}}} \cdot \left( {\frac{{{P_{{\text{gt}},t}}}}{{{\eta _{{\text{gt}}}}}} + \frac{{{H_{{\text{gb}},t}}}}{{{\eta _{{\text{gb}}}}}}} \right) $$ (14)

      式中:$ {\rho _{{\text{gas}}}} $为天然气的单位热值价格;$ {P_{{\text{gt}},t}} $$ {\eta _{{\text{gt}}}} $分别为燃气轮机$t$时刻的发电功率和发电效率; $ {H_{{\text{gb}},t}} $$ {\eta _{{\text{gb}}}} $分别为燃气锅炉$ t$时刻的发热功率和制热效率。

      ${C_{\text{e}}}$为向电网购电的费用,其计算公式为

      $$ {C_{\text{e}}} = \sum\limits_{t = 1}^T {{P_{{\text{net}},t}}} \cdot {\rho _{{\text{net}},t}} $$ (15)

      式中:$ {P_{{\text{net}},t}} $$t$时刻向电网购电的功率;$ {\rho _{{\text{net}},t}} $$t$ 时刻的购电价格。

      ${C_{\text{m}}}$为设备年维护费用,其计算公式为

      $$ \sum\limits_{i \in I} {\sum\limits_{t = 1}^T {\left| {{P_{i,t}}} \right|} } \cdot {\rho _{i,{\text{m}}}} $$ (16)

      式中:$ {P_{i,t}} $为设备$i$$t$时刻的出力值;$ {\rho _{i,{\text{m}}}} $为设备$i$的单位功率维修费用。

    • 1)电能平衡约束。

      $$ \begin{gathered} {P_{{\text{v}},t}} + {P_{{\text{w}},t}} + {P_{{\text{dis}},t}} + {P_{{\text{cha}},t}} + {P_{{\text{net}},t}} + {P_{{\text{gt}},t}} = \\ E_{{\text{DC}}}^t + {P_{{\text{car}},t}} + {P_{{\text{ec}},t}} \\ \end{gathered} $$ (17)

      式中:$ {P_{{\text{v}},t}} $$ {P_{{\text{w}},t}} $分别表示$t$时刻的光伏出力和风机出力;$ {P_{{\text{car}},t}} $表示$t$时刻3类电动汽车的总负荷;$ {P_{{\text{ec}},t}} $表示$t$时刻电制冷机消耗的电功率。

      2)冷功率平衡约束。

      $$ {C_{{\text{ec}},t}} + {C_{{\text{ac}},t}} = C_{{\text{DC}}}^t $$ (18)

      式中:$ {C_{{\text{ec}},t}} $$ {C_{{\text{ac}},t}} $分别表示$t$时刻电制冷机的制冷功率和吸收式制冷机的制冷功率。

      3)电动汽车数量和功率约束。

      $$ \begin{cases}\begin{array}{l}{n}_{1,t}+{n}_{2,t}+{n}_{3,t}\leqslant{n}_{\text{cha}}\\ {P}_{\text{car},t}\leqslant{n}_{\text{cha}}\cdot {P}_{\text{cr}}\end{array}\end{cases} $$ (19)

      式中:$ {n_{1,t}} $$ {n_{2,t}} $$ {n_{3,t}} $分别表示在$t$时刻并网的1类、2类、3类电动汽车数量;$ {n_{{\text{cha}}}} $为充电桩的数量。

      4)燃气轮机约束。

      $$ \left\{ \begin{gathered} {P_{{\text{gt}},t}} = {P_{{\text{gas}},{\text{gt}},t}} \cdot {\eta _{{\text{gt}}}} \\ {H_{{\text{gt}},t}} = {P_{{\text{gt}},t}} \cdot {\zeta _{{\text{gt}}}} \\ \end{gathered} \right. $$ (20)

      式中:$ {P_{{\text{gas}},{\text{gt}},t}} $$t$时刻燃气轮机消耗的天然气功率;$ {\zeta _{{\text{gt}}}} $为燃气轮机的热电比。

      5)吸收式制冷机约束

      $$ {C_{{\text{ac}},t}} = {H_{{\text{gt}},t}} \cdot {\eta _{{\text{ac}}}} $$ (21)

      式中:$ {\eta _{{\text{ac}}}} $表示吸收式制冷机的制冷效率;$ {H_{{\text{gt}},t}} $表示燃气轮机的余热。

      6)电制冷机约束。

      $$ {C_{{\text{ec}},t}} = {P_{{\text{ec}},t}} \cdot {\eta _{{\text{ec}}}} $$ (22)

      式中:$ {\eta _{{\text{ec}}}} $表示电制冷机的制冷效率。

      此外,优化过程所要满足的约束条件还包含式(1)—(11)。

    • 以上模型的目标是能源站与数据中心的总福利最大化,并未涉及费用的分摊。本文采用边际价格模型来对数据中心的费用进行结算,因其可以客观地反映数据中心的可平移负荷资源在协同优化中起到的降低总成本的作用。能源站应向数据中心征收的费用为

      $$ {C_{{\text{DC}}}} = \sum\limits_{t = 1}^T {{\pi _{e,t}}} \widehat {E_{{\text{DC}}}^t} + \sum\limits_{t = 1}^T {{\pi _{\phi ,t}}} \widehat {C_{{\text{DC}}}^t} $$ (23)

      式中:$ {\pi _{e,t}} $为协同优化模型中固定变量$ E_{{\text{DC}}}^t $为其最优解$ \widehat {E_{{\text{DC}}}^t} $所得的影子价格,代表t时段的边际电价;$ {\pi _{\phi ,t}} $为协同优化模型中固定$ C_{{\text{DC}}}^t $$ \widehat {C_{{\text{DC}}}^t} $所得的影子价格,代表t时段的边际冷价。

    • 公式(1)—(23)可以建模为两阶段鲁棒优化模型:

      $$ \mathop {\min }\limits_{x,y} \left( {{{\mathbf{c}}^{\text{T}}}x + \mathop {\max }\limits_{{\xi } \in U} \;{{\mathbf{b}}^{\text{T}}}y} \right) $$ (24)

      约束条件为

      $$ \begin{array}{c}x\in X\\ y\in Y\\ {\bf{A}}x+{\bf{B}}y\leqslant{\bf{g}}\\ {\bf{T}}\xi +{\bf{C}}x+{\bf{D}}y={\bf{h}}\end{array} $$ (25)

      式中:$ x $$ y $分别为第1阶段和第2阶段的决策变量,对应的可行域分别是$ X $$ Y $,第2阶段的决策变量$ y $与不确定量$\xi $以及第1阶段决策变量$ x $有关。$ {c} $$ {b} $为成本系数矩阵,$ {\mathbf{A}} $$ {\mathbf{B}} $$ {\mathbf{g}} $$ {\mathbf{T}} $$ {\mathbf{C}} $$ {\mathbf{D}} $$ {\mathbf{h}} $是与约束相关的常数矩阵,均可通过式(12)—(23)推得。公式(24)中的max代表的含义是不确定量$ \xi $在不确定集$U$所能取到的最恶劣的情况。假设第2阶段变量$ y $与不确定量$ {\mathbf{\xi }} $呈线性关系,即$y = W\xi $,其中$W$即为设备的参与因子。可以将上述min-max问题转化为:

      $$ \mathop {\min }\limits_{_{{\mathbf{x,y}}}} \left( {{{\mathbf{c}}^{\text{T}}}x + Q} \right) $$ (26)

      约束条件为:

      $$ x \in X $$ (27)
      $$ y \in Y $$ (28)
      $$ \underset{{}_{\xi \in \text{U}}}{\mathrm{max}}\;{\bf{b}}^{\text{T}}(W\xi )\leqslant Q $$ (29)
      $$ \underset{{}_{\xi \in \text{U}}}{\mathrm{max}}\;{\bf{A}}x+{\bf{B}}W\xi \leqslant{\bf{g}} $$ (30)
      $$ {\bf{T}}\xi {\mathbf{ + }}{\bf{C}}x{\mathbf{ + }}{\bf{D}}W\xi {\mathbf{ = }}{\bf{h}} $$ (31)

      公式(24)—(31)即为原优化问题公式(12)—(23)的紧凑形式。其中$Q$为辅助变量。

    • 紧凑模型公式(12)—(23)中有双线性项$ W\xi $,无法直接求解,需要将双线性项进行一些数学处理,如引入松弛变量[10];或者通过列生成逐步排除不符合约束的解 [11],但需要求解另一个非线性问题;本文采用对偶理论,对约束进行转换[12]。假设不确定集$U$为:

      $$ -\Delta \xi \leqslant\xi \leqslant\Delta \xi $$ (32)

      式中:$ {\bf{\xi }} $表示风电或者光伏出力的不确定量,而$ \Delta {\bf{\xi }} $表示与基准预测出力的最大偏差量。根据对偶理论,式(29)可以转换为:

      $$ \underset{{}_{r,s}}{\mathrm{min}}\;{\bf{b}}^{\text{T}}[(r+s)\Delta \xi ]\leqslant Q $$ (33)
      $$ {\text{ s}}{\text{.t}}{\text{. }}\quad r - s = {{\bf{b}}^{\text{T}}}W $$ (34)

      式中:$ {\mathbf{r}} $$ {\mathbf{s}} $为非负的对偶变量。而大于等于一个最小值,即等价于存在性问题,因此可以把式(29)转化为

      $$ {\bf{b}}^{\text{T}}[(r+s)\Delta \xi ]\leqslant Q $$ (35)
      $$ r - s = {{\bf{b}}^{\text{T}}}W $$ (36)

      式(30)的转化与上述步骤相同,不再赘述。而对于式(31)的等式约束,可以转化为:

      $$ {\bf{T}} + {\bf{D}}W = 0 $$ (37)
      $$ {\bf{C}}x - {\bf{h}} = 0 $$ (38)

      经过处理之后模型公式(12)—(23)成为一个线性模型,可以调用成熟的商业求解器进行求解。

    • 以某数据中心的负荷需求为参考[2],选取典型日的负荷数据作为多站融合场景下数据中心的冷、电能需求数据,如图2所示。并且选取典型日的光伏、风机出力曲线,如图3所示。而投资方向电网购电的价格如图4所示。各个设备的参数如附录A表A1所示,其中储能站的SOC下限值已经经过计算,保证可以维持数据中心以最大负荷正常运行30 min。在充电站方面,第1类、第2类、第3类电动汽车的数量分别设置为12辆、20辆、30辆,充电价格分别为1.2元/kWh、1元/kWh和0.8元/kWh。考虑到电动汽车并网的随机性,其并网时刻由计算机随机生成。但是考虑实际情况,第1类电动汽车对应车主用车较急的情况,因此第1类电动汽车的并网时刻设置在7:00—21:00之间;而第2、第3类电动汽车的并网时刻则是在1:00—23:00之间。在并网时刻生成时,还要考虑充电桩数量约束的限制。本文设置充电桩的数量为10。而电动汽车并网时的SOC值也是随机生成,考虑到实际情况,规定并网时的SOC在0.1~0.4之间。

      图  2  数据中心负荷预测

      Figure 2.  Load prediction of data center

      图  3  典型日光伏、风机预测出力曲线

      Figure 3.  Forecasted output of photovoltaic generation and wind generator in a typical day

      图  4  购电价格

      Figure 4.  Purchase/sell electricity price

      表 A1  设备参数

      Table A1.  Parameters of equipment

      设备名称燃气轮机燃气锅炉吸收式制冷机电制冷机电储能光伏风机电动汽车充电桩
      投资费用/(元/kW)282662812569501800753680001000
      维修费用/(元/kW·h)0.0940.0030.0140.0080.0060.0120.015
      发电效率0.43
      热电比1.03
      制热效率0.92
      制冷效率1.13.5
      充电效率0.950.9
      放电效率0.950.9
      容量150050010001000150040040060
      soc下限/并网soc0.47≤0.4
      soc上限/离网soc0.951
      最大功率/kW3003030

      算例中对比以下3种模型,以验证所提多站融合模式和仿射鲁棒优化的有效性。

      模型1 :数据中心与能源站独立运行,采用确定性的调度方法。

      模型2:数据中心和能源站协同优化,采用确定性的调度方法。

      模型3:数据中心和能源站协同优化,采用仿射鲁棒优化方法。

    • 首先基于前一节的输入数据验证模型3,即采用仿射鲁棒优化方法的能源站与数据中心协同优化,考虑目标函数式(12)与约束条件式(7)—(11)以及式(17)—(22),在Python中调用Docplex包进行计算,得到最优调度结果。图5展示了某一次随机生成的电动汽车并网数量情况,图6 (a)和(b)分别展示了能源站的最优供电和供冷计划,在电价低谷时段,储能充电,燃气轮机不发电,电制冷机开启,而在电价平时段和高峰时段,储能放电,燃气轮机发电,电制冷机关闭,从而实现峰谷套利。

      图  5  各个时刻电动汽车并网数量

      Figure 5.  The number of electric vehicles connected to the grid at each moment

      图  6  能源站电、冷优化调度结果

      Figure 6.  Optimal dispatch result of power, heating and cooling load in energy station

      在充电站方面,由于电动汽车的并网时刻和并网SOC是随机生成的,每次优化计算的结果都不相同,图7展示了平均意义下优化后的电动汽车充放电情况,正表示充电,负表示放电。电动汽车的充放电总体与并网数量相关,这是因为电动汽车的充电总时长存在限制,无法实现电储能的长时间范围充放电。然而,在短时间范围内,如低谷时段0:00—6:00,电动汽车可以利用电价低谷尽量多充电,而在某些高峰时段,如8:00—9:00、15:00、21:00,电动汽车甚至可以反向送电,从而降低该时段向电网购电的费用。

      图  7  电动汽车充放电情况

      Figure 7.  Charging and discharging conditions of electric vehicles

      在应对新能源出力波动方面,各类设备均参与了调节,其参与因子如图8所示。吸收式制冷机、燃气轮机、储能、电网的参与因子为负,这是因为当新能源出力大于预测值时,需要减少这些设备的出力来维持功率平衡;电制冷机的参与因子为正,这是因为当新能源出力大于预测值时,需要增加电制冷机的供冷量来消纳富余的新能源。

      图  8  设备参与因子优化结果

      Figure 8.  Optimization results of equipment participating factors

    • 图9图10展示了数据中心在调度前后电负荷和冷负荷的对比。在电价低谷时段,如4:00—7:00、12:00—14:00和23:00—24:00,数据中心将相邻时段的计算任务调节至这些时段以增加电负荷和冷负荷。在电价高峰时段,如8:00—11:00和20:00—22:00,数据中心计算任务调节至相邻时段以减少电负荷和冷负荷。

      图  9  数据中心电负荷调度前后对比

      Figure 9.  Comparison of power load of data center before and after dispatching

      图  10  数据中心冷负荷调度前后对比

      Figure 10.  Comparison of cooling load of data center before and after dispatching

      对于数据中心来说,由于从能源站购能价格较低,降低了其运行成本。假设数据中心直接通过向电网购电满足其电负荷和制冷需求,则单日的购电成本为23675元。而在多站融合场景下,通过向设备投资方购电和购冷的成本为20184元,降低了17.3%。此外,因为有储能站保证了数据中心的备用容量,以及多种能源供应方式带来的可靠性的提高,数据中心不需要再投资UPS,节省了其投资费用。

    • 表1展示了在10000组随机生成的新能源出力情况时不同模型下能源站的收益情况,其中假设新能源出力服从标准差为预测值的10%的正态分布。模型1采用了确定性优化方法,能源站与数据中心独立运行,由于充电站负荷相对较低,能源站有富余的电能向主网出售,因此购电成本小于0,但因为缺少了数据中心的收入,因此总收益较低。在能源站与数据中心协同优化后(模型2),增加了购电成本和购气成本,但因此获得了更多的数据中心收入,因此总收益较高。进一步采用仿射鲁棒优化方法后(模型3),由于系统应对不确定性的能力增强,相对减少了购电和购气成本,因此总收益进一步提高。

      表 1  “多站融合”成本分析

      Table 1.  Cost analysis on multi-station integration

      模型
      编号
      购电
      成本
      /(元/日)
      购气
      成本
      /(元/日)
      运维
      成本
      /(元/日)
      充电站
      收益
      /(元/日)
      数据中
      心收益
      /(元/日)
      净收益
      /(元/日)
      1−30120153263805497
      210593463213202703202936451
      310307461813072712201846664

      图11展示了电网公司直接向数据中心供电时的交换功率和多站融合场景下与电网公司的交换功率情况。可以看到,设备投资方为了使得自身运行成本尽可能降低,其与电网交换功率大于0的时刻都集中在平时段或者谷时段。并且因为储能站的存在,在谷时段的交换功率会大于数据中心的电负荷,把电储存起来,在峰时段释放。因此如果把多站融合作为一个负荷整体来看,其负荷特性恰好发挥了需求侧响应的效果,具有削峰填谷的作用。可以使电网公司避免应对峰时段的高负荷而增加备用投资,节约了成本。

      图  11  2种场景下与电网交换功率情况

      Figure 11.  Power exchange with the grid in two scenarios

      从环境的角度来说,因为每个时刻尽量利用新能源出力可以节约购电费用和购气费用,所以弃风弃光率为0,光伏、风电的消纳率很高,避免了资源浪费。因此多站融合的模式对于电网公司、设备投资主体以及数据中心都有益处,具有一定的发展前景。

      此外,为了验证在多站融合的框架下,数据中心的可靠性要求是否能达到要求,本文固定了各个设备的可用度,如表2所示[3]。通过蒙特卡洛模拟方法,验证数据中心的供电可靠性。注意到供电的两路电网线路,其中一路在正常运行时不向数据中心供电,只有在发生事故时才会进行供电;储能站SOC在正常优化运行时的下限为0.47,以保证数据中心以最大负荷正常运行30 min,但是在计算可靠性时,当设备发生故障导致供能受到影响时,应该将储能站的SOC下限值设为0,以此来计算数据中心供电可靠性。本文选取蒙特卡洛模拟的次数为3000次,计算得到其切负荷的概率为0.999563,满足其可靠性要求[7]。因此不必投资建设柴油发电机和不间断电源,进一步降低投资成本。

      表 2  各个设备的可用度

      Table 2.  Availability of each device

      设备可用度
      燃气轮机 0.991293
      吸收式制冷机0.995732
      电制冷机0.996898
      电储能0.999830
      光伏0.9999
      风机0.9999
      断路器0.999998
    • 本文以设备投资方的利润最大为目标函数,考虑风电和光伏的不确定性建立了多能电站的仿射鲁棒最优调度模型。仿真结果表明,多站融合的模式对于电网、设备投资方以及数据中心都有益处,在满足数据中心供电可靠性的要求的同时,在新能源消纳方面也具有较好的表现。

      本文只是考虑了风电和光伏出力的不确定性。在对电动汽车负荷的建模上也只考虑了随机性,没有考虑不确定性。因此如何对电动汽车不确定性进行建模,并与风电、光伏出力的不确定性相结合,并在此基础上进行优化计算是下一步的研究方向。

      (本刊附录请见网络版,印刷版略)

参考文献 (12)

目录

    /

    返回文章
    返回