• 匿名盲审
  • 学术期刊非营利性
  • 全球免费开放获取全文
  • 最新科研成果提供绿色通道

留言板

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

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

考虑CVaR的机组组合和多场景备用决策联合优化

李建钊 谢敏 李舒佳 林盛振 黄彬彬

李建钊,谢敏,李舒佳等.考虑CVaR的机组组合和多场景备用决策联合优化[J].南方能源建设,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
引用本文: 李建钊,谢敏,李舒佳等.考虑CVaR的机组组合和多场景备用决策联合优化[J].南方能源建设,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
LI Jianzhao,XIE Min,LI Shujia,et al.Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR[J].Southern Energy Construction,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
Citation: LI Jianzhao,XIE Min,LI Shujia,et al.Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR[J].Southern Energy Construction,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008

考虑CVaR的机组组合和多场景备用决策联合优化

doi: 10.16516/j.gedi.issn2095-8676.2021.04.008
基金项目: 

广东省自然科学基金面上项目“虚拟电厂动态聚合演化理论与计算方法研究” 2021A1515012245

详细信息
    作者简介:

    李建钊(通信作者)1993-,男,广西北海人,华南理工大学电力系统及其自动化硕士研究生,主要从事电力市场、电力系统运行优化与控制(e-mail)178048863@qq.com

    谢敏1978-,女,湖南衡阳人,博士,副教授,主要研究方向为电力市场、电力系统运行优化与控制(e-mail)minxie@scut.edu.cn

    李舒佳1996-,女,江西萍乡人,华南理工大学电气工程硕士研究生,主要研究方向为电力市场、虚拟电厂运营机制(e-mail)417520690@qq.com

  • 中图分类号: TM7

Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR

  • LI Jianzhao,XIE Min,LI Shujia,et al.Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR[J].Southern Energy Construction,2021,08(04):50-65.
  • 摘要:   目的  为提高可再生能源消纳水平以及衡量其出力不确定性带来的风险,该文以特性各异电源互补协调调度作为提升系统运行经济性和促进可再生能源消纳的手段,并针对系统内可再生出力和负荷的波动性,引入条件风险价值(conditional value-at-risk,CVaR)理论,量化其不确定性对系统调度运行时产生的风险。  方法  利用多场景技术模拟系统不确定性,建立考虑条件风险价值(CVaR)的机组组合和多场景备用决策联合优化模型,将风险水平限制在可接受的前提下,追求系统运行成本最小。最后在通用代数建模系统(the general algebraic modeling system, GAMS)平台上编程,并调用混合整数线性规划CPLEX求解器进行求解。  结果  仿真计算表明,不同类型电源的机组组合和备用决策能有效提高系统运行经济性、可靠性和可再生能源消纳水平。  结论  基于改进IEEE 24和某省实际电网的算例分析验证了所提模型和方法的合理性和有效性。
  • 图  1  系统调度框架

    Fig.  1  Frame of system dispatching

    图  2  计算流程图

    Fig.  2  Flow chart of the calculation

    图  3  风电场景

    Fig.  3  Wind power scenarios

    图  4  负荷场景

    Fig.  4  Load scenarios

    图  5  方案1各时段机组出力

    Fig.  5  Unit output of scheme 1

    图  6  方案4各时段机组出力

    Fig.  6  Unit output of scheme 4

    图  7  方案1和方案4各时段非核机组的启动台数

    Fig.  7  Number of non-nuclear generators in operation in each period of scheme 1 and scheme 4

    图  8  不同风电渗透率下的总调度成本和CVaR的有效边界

    Fig.  8  Efficient frontier of Dispatching CVaR under different penetration of wind power

    图  9  不同风险系数下系统总调度成本和CVaR的有效边界

    Fig.  9  Efficient frontier of dispatching cost vs CVaR under different risk parameters

    图  10  负荷场景

    Fig.  10  Load scenarios

    图  11  风电1场景

    Fig.  11  Scenarios of No.1 wind farm

    图  12  风电2场景

    Fig.  12  Scenarios of No.2 wind farm

    表  1  风电和负荷场景概率

    Tab.  1.   Probability of wind power and load scenarios

    风电场景S概率ps负荷场景S概率ps
    S10.109S10.270
    S20.198S20.160
    S30.055S30.273
    S40.059S40.114
    S50.047S50.183
    S60.026
    S70.207
    S80.038
    S90.109
    S100.152
    下载: 导出CSV

    表  2  不同调度方案

    Tab.  2.   Dispatch of different schemes

    方案火电水电可中断负荷核电储能
    方案1+++
    方案2++++
    方案3++++
    方案4+++++

    注: +代表参与提供备用(核电参与调峰);-代表不参与。

    下载: 导出CSV

    表  3  不同调度方案下的系统运行指标

    Tab.  3.   System performance indices of different schemes

    运行指标方案1方案2方案3方案4
    总发电调度成本811 327.3783 272.3781 827.0778 702.2
    预想场景调度成本743 327.1761 252.3745 026.7757 102.1
    误差场景期望调度成本68 000.222 020.036 800.421 600.1
    备用容量期望成本51 160.920 682.633 544.621 600.1
    弃风损失期望16 466.01 337.43 255.80.0
    失负荷损失期望373.30.00.00.0
    弃风、失负荷损失期望16 839.31 337.43 255.80.0
    CVaR813 409.4783 334.2781 857.2778 724.3
    下载: 导出CSV

    表  4  方案1各时段机组启动计划及启动成本

    Tab.  4.   Unit commitment planning and cost of schemes 1

    机组启停状态机组编号时间段成本/$
    启动1t81 430.4
    3t81 725.0
    4t103 056.7
    6t8312.0
    9t11624.0
    10t102 298.0
    停机1t4715.2
    3t16862.5
    4t221 528.35
    6t18156.0
    7t24156.0
    9t3312.0
    10t11 149.0
    下载: 导出CSV

    表  5  方案4各时段机组启停计划及启停成本

    Tab.  5.   Unit commitment planning and cost of schemes 4

    机组启停状态机组编号时间段成本/$
    启动4t103 056.7
    6t9312.0
    停机4t221 528.35
    6t18156.0
    下载: 导出CSV

    表  6  不同风电渗透率下方案1和方案4计算结果

    Tab.  6.   Calculation results of scheme 1 and scheme 4 under different wind power penetration

    风电容量/MW渗透率/%方案预测场景调度成本/$备用容量期望成本/$弃风、失负荷损失期望/$
    3006.8方案1823 275.37 709.4100.7
    3006.8方案4822 871.95 242.80
    60013.6方案1798 662.39 723.40
    60013.6方案4796 863.39 395.80
    90020.3方案1773 252.818 031.81 523.8
    90020.3方案4768 538.319 214.02 327.8
    1 20027.1方案1743 327.151 160.916 839.3
    1 20027.1方案4757 102.121 600.10
    风电容量/MW渗透率/%方案误差场景期望调度成本/$总调度成本/$CVaR/$
    3006.8方案17 810.1831 085.5832 058.0
    3006.8方案45 242.8828 114.7828 165.1
    60013.6方案19 723.4808 385.7808 501.4
    60013.6方案49 395.8806 259.1806 360.7
    90020.3方案119 555.6792 808.4792 826.1
    90020.3方案421 541.8790 080.0790 080.0
    1 20027.1方案168 000.2811 327.3813 409.4
    1 20027.1方案421 600.1778 702.2778 724.3
    下载: 导出CSV

    表  7  方案4总调度成本和CVaR降幅

    Tab.  7.   Total dispatching cost and CVaR reduction of scheme 4

    渗透率/%总调度成本/%CVaR/%
    6.8
    13.62.642.63
    20.34.594.60
    27.15.975.97
    下载: 导出CSV

    表  8  不同风险系数下的计算结果

    Tab.  8.   Calculation results under different risk parameters

    风险系数β预想场景调度成本误差场景期望调度成本总发电调度成本CVaR
    0737 365.720 960.1758 325.8799 213.3
    0.1739 283.019 192.8758 475.7795 472.6
    0.2741 279.517 629.8758 909.4793 106.0
    0.3743 618.316 072.8759 691.1790 768.8
    0.4748 617.113 694.4762 311.5785 779.5
    0.5754 022.312 337.9766 360.2780 875.7
    0.6756 267.512 139.0768 406.5779 036.2
    0.7756 858.512 006.1768 864.6778 800.1
    0.8757 065.011 982.9769 047.8778 738.1
    0.9757 102.112 014.4769 116.5778 724.3
    1757 102.121 600.1778 702.2778 724.3
    下载: 导出CSV

    表  9  不同调度方案下的系统运行指标

    Tab.  9.   System performance indices of different schemes

    运行指标方案1方案4
    总发电调度成本599 682 595598 848 191
    预想场景调度成本592 098 821592 006 420
    误差场景期望调度成本7 583 7746 841 771
    备用容量期望成本7 583 7746 841 771
    弃风损失期望0.00.0
    失负荷损失期望0.00.0
    弃风、失负荷损失期望0.00.0
    CVaR599 974 173599 038 682
    下载: 导出CSV

    表  10  不同风险系数下的计算结果

    Tab.  10.   Calculation results under different risk parameters

    风险系数β预想场景调度成本误差场景期望调度成本总发电调度成本CVaR
    0591 506 8547 087 759598 594 613599 782 842
    0.1591 799 1656 925 843598 725 008599 548 623
    0.2591 874 1696 921 227598 795 396599 448 130
    0.3591 911 4556 902 821598 814 276599 116 590
    0.4591 952 0636 876 704598 828 767599 069 834
    0.5592 006 4206 841 771598 848 191599 038 682
    0.6592 056 1746 826 255598 882 429599 015 802
    0.7592 087 3246 802 164598 889 488599 010 729
    0.8592 139 4616 784 591598 924 052599 008 430
    0.9592 237 6616 764 004599 001 665598 856 317
    1592 237 9366 892 163599 130 099598 853 578
    下载: 导出CSV

    表  A1  发电机技术参数

    Tab.  A1.   Technical parameters of generating units

    Unit #12345678910
    Node12713151516222323
    Offer cost/[$·(MWh)-111.9611.9620.0320.3223.7810.2510.25010.2510.66
    Λup/$1 430.41 430.41 725.03 056.7437.0312.0312.00.0624.02 298.0
    Pmin/(MW)30.430.475206.851254.2554.25300108.5140
    Pmax/(MW)15215230059160155155300310350
    Rdn / Rup/(MW·h-115215230054060155155300310240
    Rdn / Rup/(MW·h-1121.6121.6225384.1548100.75100.75300201.5210
    Sdn / Sup/(MW·h-115215230054060155155300310240
    DT/h448102880848
    UT/h888124880824
    S0/h0021120000
    U0/h22220000102410300
    u0/h1100001111
    下载: 导出CSV

    表  A2  输电线路电抗和容量

    Tab.  A2.   Reactance and capacity of transmission lines

    Form nodeTo nodeReactance/(p.u.)Capacity/(MVA)FormnodeTonodeCapacity/(MVA)Capacity/(MVA)FormnodeCapacity/(MVA)Reactance/(p.u.)Capacity/(MVA)
    120.014 62458100.176 224515210.024 91 400
    130.225 32459110.084 056015240.052 9700
    150.090 72459120.084 056016170.026 3700
    240.135 624510110.084 056016190.023 4700
    260.205 024510120.084 056017180.014 3700
    390.127 124511130.048 870017220.106 9700
    3240.084 056011140.042 670018210.013 21 400
    490.111 024512130.048 870019200.020 31 400
    5100.094 024512230.098 570020230.011 21 400
    6100.064 224513230.088 470021220.069 2700
    780.065 224514160.059 4700
    890.176 224515160.017 2700
    下载: 导出CSV

    表  A3  节点位置和系统总需求分布

    Tab.  A3.   Node location and distribution of the total system demand

    Load #Node #% of system loadLoad #Node #% of system load
    113.810106.8
    223.411139.3
    336.312146.8
    442.6131511.1
    552.514163.5
    664.8151811.7
    774.416196.4
    886.017204.5
    996.1
    下载: 导出CSV

    表  A4  风电、核电和储能位置分布

    Tab.  A4.   Location and distribution of wind power, nuclear power and energy storage

    UnitNode #Rated power/MW
    Wind farm #120300
    Wind farm #28300
    Wind farm #36300
    Wind farm #47300
    Nuclear plant #118900
    Nuclear plant #2211 000
    Energy storage #17300
    下载: 导出CSV
  • [1] 舒印彪, 张智刚, 郭剑波, 等. 新能源消纳关键因素分析及解决措施研究 [J]. 中国电机工程学报, 2017, 37(1): 1-9.

    SHUY B, ZHANGZ G, GUOJ B, et al. Study on key factors and solution of renewable energy accommodation [J]. Proceedings of the CSEE, 2017, 37(1): 1-9.
    [2] 中国电力企业联合会. 中国电力行业年度发展报告2018 [R]. 北京: 中国电力企业联合会, 2018.

    China Electricity Council. Annual development report of China's power industry 2018 [R]. Beijing: China Electricity Council, 2018.
    [3] 中国电力企业联合会. 2017-2018年度全国电力供需形势分析预测报告 [J]. 电器工业, 2018(2): 11-15.

    China Electricity Council. Analysis and prediction report of national power supply and demand situation in 2017-2018 [J]. China Electrical Equipment Industry, 2018(2): 11-15.
    [4] RAU G H, WILLAUERH D, RENZ J. The global potential for converting renewable electricity to negative-CO2-emissions hydrogen [J]. Nature Climate Change, 2018, 8(7): 621-625.
    [5] CALDEIRAK, JAINA K, HOFFERTM I. Climate sensitivity uncertainty and the need for energy without CO2 emission [J]. Science, 2003, 299(5615): 2052-2054.
    [6] OMERA M. Energy, environment and sustainable development [J]. Renewable and Sustainable Energy Reviews, 2007, 12(9): 2265-2300.
    [7] 周玮, 彭昱, 孙辉, 等. 含风电场的电力系统动态经济调度 [J]. 中国电机工程学报, 2009, 29(25): 13-18.

    ZHOUW, PENGY, SUNH, et al. Dynamic economic dispatch in wind power integrated system [J]. Proceedings of the CSEE, 2009, 29(25): 13-18.
    [8] 戴蕾思, 叶承晋, 傅旭华, 等. 考虑概率分布约束的含高渗透率风电电力系统储能鲁棒优化方法 [J]. 电网技术, 2017, 41(3): 769-777.

    DAIL S, YEC J, FUX H, et al. Distributional robust joint chance constrained optimal capacity installment of energy storage in power system with high penetration of wind power [J]. Power System Technology, 2017, 41(3): 769-777.
    [9] 陈厚合, 王杨, 张儒峰, 等. 考虑源荷协调的风电并网系统旋转备用容量优化 [J]. 电力自动化设备, 2017, 37(8): 185-192.

    CHENH H, WANGY, ZHANGR F, et al. Spinning reserve capacity optimization considering coordination between source and load for power system with wind power [J]. Electric Power Automation Equipment, 2017, 37(8): 185-192.
    [10] XIAOX, ZHOUY H, HES Y, et al. Optimal joint operation of hydro-thermal-wind hybrid power system with cascaded hydro power [J]. Electric Power Automation Equipment, 2018, 38(2): 100-108.
    [11] 刘宇穗. 新型核电机组启停及给水系统控制策略 [J]. 南方能源建设, 2020, 7(2): 127-131.

    LIUY S. Control strategy of start/stop and feed water system of new type nuclear power plant [J]. Southern Energy Construction, 2020, 7(2): 127-131.
    [12] DONGZ, FENGJ T, HUANGX J. Nonlinear observer-based feedback dissipation load-following control for nuclear reactors [J]. IEEE Transactions on Nuclear Science, 2009, 56(1): 272-285.
    [13] 王骏, 赵洁, 刘涤尘, 等. 考虑核电参与的调峰优化运行模型 [J]. 中国电机工程学报, 2018, 38(6): 1665-1674+1903.

    WANGJ, ZHAOJ, LIUD C, et al. Optimal scheduling model of peak load regulation considering participation of nuclear power plant [J]. Proceedings of the CSEE, 2018, 38(6): 1665-1674+1903.
    [14] 隋鑫, 卢盛阳, 苏安龙, 等. 计及风电和柔性负荷的核电多目标优化调度研究 [J]. 中国电机工程学报, 2019, 39(24): 7232-7241+7496.

    SUIX, LUS Y, SUA L, et al. Research on multi-objective optimal scheduling of nuclear power considering wind power and flexible load [J]. Proceeding of the CSEE, 2019, 39(24): 7232-7241+7496.
    [15] GÖNSCHJ, HASSLERM. Optimizing the conditional value-at-risk in revenue management [J]. Review of Managerial Science, 2014, 8(4): 495-521.
    [16] ARTZNERP, DELBAENF, EBERJ M, et al. Coherent measures of risk [J]. Mathematical Finance, 1999, 9(3): 203-228.
    [17] WUJ L, ZHANGB H, DENGW S, et al. Application of cost-CVaR model in determining optimal spinning reserve for wind power penetrated system [J]. International Journal of Electrical Power and Energy Systems, 2015, 66(3): 110-115.
    [18] 王旭, 蒋传文, 刘玉娇, 等. 含大规模风电场的电力系统旋转备用风险决策模型和方法 [J]. 电力系统自动化, 2014, 38(13): 64-70.

    WANGX, JIANGC W, LIUY J, et al. Decisionmaking model and method for spinning reserve and risk of power systems integrated with large-scale wind farms [J]. Automation of Electric Power Systems, 2014, 38(13): 64-70.
    [19] 吕梦璇, 娄素华, 刘建琴, 等. 含高比例风电的虚拟电厂多类型备用协调优化 [J]. 中国电机工程学报, 2018, 38(10): 2874-2882+3138.

    LÜM X, LOUS H, LIUJ Q, et al. Coordinated optimization of multi-type reserve in virtual power plant accommodated high shares of wind power [J]. Proceedings of the CSEE, 2018, 38(10): 2874-2882+3138.
    [20] TANZ F, WANGG, JUL W, et al. Application of CVaR risk aversion approach in the dynamical scheduling optimization model for virtual power plant connected with wind-photovoltaic-energy storage system with uncertainties and demand response [J]. Energy, 2017(124): 198-213.
    [21] 王海冰, 王承民, 张庚午, 等. 考虑条件风险价值的两阶段发电调度随机规划模型和方法 [J]. 中国电机工程学报, 2016, 36(24): 6838-6848+6939.

    WANGH B, WANGC M, ZHANGG W, et al. Two-stage stochastic generation dispatching model and method considering conditional value-at-risk [J]. Proceedings of the CSEE, 2016, 36(24): 6838-6848+6939.
    [22] REZAEIN, KALANTARM. Stochastic frequency-security constrained energy and reserve management of an inverter interfaced islanded microgrid considering demand response programs [J]. International Journal of Electrical Power and Energy Systems, 2015(69): 273-286.
    [23] 周任军, 姚龙华, 童小娇, 等. 采用条件风险方法的含风电系统安全经济调度 [J]. 中国电机工程学报, 2012, 32(1): 56-63+18.

    ZHOUR J, YAOL H, TONGX J, et al. Security economic dispatch in wind power integrated systems using a conditional risk method [J]. Proceedings of the CSEE, 2012, 32(1): 56-63+18.
    [24] 易国伟, 童小娇, 周鹏, 等. CVaR和EVaR安全运行风险管理下的电力系统经济调度 [J]. 电力系统保护与控制, 2016, 44(8): 49-56.

    YIG W, TONGX J, ZHOUP, et al. Power system economic dispatch under CVaR and EVaR security operation risk management [J]. Power System Protection and Control, 2016, 44(8): 49-56.
    [25] 苏鹏, 刘天琪, 李兴源. 含风电的系统最优旋转备用的确定 [J]. 电网技术, 2010, 34(12): 158-162.

    SUP, LIUT Q, LIX Y. Determination of optimal spinning reserve of power grid containing wind [J]. Power System Technology, 2010, 34(12): 158-162.
    [26] LOH W L. On latin hypercube sampling [J]. The Annals of Statistics, 1996, 24(5): 2058-2080.
    [27] 吴雄, 王秀丽, 王建学, 等. 微网调度问题的混合整数规划方法 [J]. 中国电机工程学报, 2013, 33(28): 1-8.

    WUX, WANGX L, WANGJ X, et al. Economic generation scheduling of a microgrid using mixed integer programming [J]. Proceedings of the CSEE, 2013, 33(28): 1-8.
    [28] GROWEK N, HEITSCHH, ROMISCHW. Scenario reduction and scenario tree construction for power management problem [C]//IEEE. Power Tech Conference Proceedings, 2003 IEEEBologna, Bologna, 2003. Bologna: IEEE, 2003: 1-7.
    [29] 梁子鹏, 陈皓勇, 雷佳, 等. 考虑风电不确定度的风-火-水-气-核-抽水蓄能多源协同旋转备用优化 [J]. 电网技术, 2018, 42(7): 2111-2119+2123.

    LIANGZ P, CHENH Y, LEIJ, et al. A multi-source coordinated spinning reserve model considering wind power uncertainty [J]. Power System Technology, 2018, 42(7): 2111-2119+2123.
    [30] 周任军, 张浩, 范文帅, 等. 基于条件风险方法的含风电系统旋转备用优化调度 [J]. 电力科学与技术学报, 2015, 30(1): 3-9.

    ZHOUR J, ZHANGH, FANW S, et al. Spinning reserve optimal dispatch in wind power integrated system with CVaR method [J]. Journal of Electric Power Science and Technology, 2015, 30(1): 3-9.
    [31] URYASEVS. Conditional value-at-risk, optimization algorithms and applications [J]. Financial Engineering News, 2000, 2(3): 1-5.
    [32] ROCKAFELLARR T, URYASEVS. Optimization of conditional value-at-risk [J]. Journal of Risk, 2000, 2(3): 21-41.
    [33] 谢敏, 韦薇, 李建钊, 等. 储能电站参与能量-调频市场联合调度模式研究 [J]. 电力建设, 2019, 40(2): 20-28.

    XIEM, WEIW, LIJ Z, et al. Research on combined dispatching mode of energy storage market [J]. Electric Power Construction, 2019, 40(2): 20-28.
    [34] 侯文庭, 韦化. 考虑核电可调度性的风-光-核-水-火多源协调短期优化调度 [J]. 电工技术学报, 2018, 33(12): 2873-2882.

    HOUW T, WEIH. A multi-source coordinated short-term dispatch model considering the dispatchability of nuclear power plants [J]. Transactions of China Electrotechnical Society, 2018, 33(12): 2873-2882.
  • 加载中
图(12) / 表 (14)
计量
  • 文章访问数:  258
  • HTML全文浏览量:  57
  • PDF下载量:  39
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-13
  • 修回日期:  2021-06-09
  • 刊出日期:  2021-12-25

考虑CVaR的机组组合和多场景备用决策联合优化

doi: 10.16516/j.gedi.issn2095-8676.2021.04.008
    基金项目:

    广东省自然科学基金面上项目“虚拟电厂动态聚合演化理论与计算方法研究” 2021A1515012245

    作者简介: 李建钊(通信作者)1993-,男,广西北海人,华南理工大学电力系统及其自动化硕士研究生,主要从事电力市场、电力系统运行优化与控制(e-mail)178048863@qq.com

    谢敏1978-,女,湖南衡阳人,博士,副教授,主要研究方向为电力市场、电力系统运行优化与控制(e-mail)minxie@scut.edu.cn

    李舒佳1996-,女,江西萍乡人,华南理工大学电气工程硕士研究生,主要研究方向为电力市场、虚拟电厂运营机制(e-mail)417520690@qq.com

  • 中图分类号: TM7

摘要:   目的  为提高可再生能源消纳水平以及衡量其出力不确定性带来的风险,该文以特性各异电源互补协调调度作为提升系统运行经济性和促进可再生能源消纳的手段,并针对系统内可再生出力和负荷的波动性,引入条件风险价值(conditional value-at-risk,CVaR)理论,量化其不确定性对系统调度运行时产生的风险。  方法  利用多场景技术模拟系统不确定性,建立考虑条件风险价值(CVaR)的机组组合和多场景备用决策联合优化模型,将风险水平限制在可接受的前提下,追求系统运行成本最小。最后在通用代数建模系统(the general algebraic modeling system, GAMS)平台上编程,并调用混合整数线性规划CPLEX求解器进行求解。  结果  仿真计算表明,不同类型电源的机组组合和备用决策能有效提高系统运行经济性、可靠性和可再生能源消纳水平。  结论  基于改进IEEE 24和某省实际电网的算例分析验证了所提模型和方法的合理性和有效性。

English Abstract

李建钊,谢敏,李舒佳等.考虑CVaR的机组组合和多场景备用决策联合优化[J].南方能源建设,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
引用本文: 李建钊,谢敏,李舒佳等.考虑CVaR的机组组合和多场景备用决策联合优化[J].南方能源建设,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
LI Jianzhao,XIE Min,LI Shujia,et al.Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR[J].Southern Energy Construction,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
Citation: LI Jianzhao,XIE Min,LI Shujia,et al.Optimal Dispatch of Unit Commitment and Multi-Scenario Reserve Decision Considering CVaR[J].Southern Energy Construction,2021,08(04):50-65. doi:  10.16516/j.gedi.issn2095-8676.2021.04.008
  • OA:https://www.energychina.press/

    近年来,我国可再生能源的装机容量在电网中的占比逐年攀升1-3,高比例可再生能源并网发电已成为我国能源发展战略的必然趋势。大规模可再生能源并网发电,一方面能有效缓解能源短缺和改善气候环境问题14-6,另一方面其出力不确定性给电力系统的优化调度带来巨大的挑战。在此背景下,传统确定性动态经济调度7已不再适用,亟需新的优化调度方法。

    目前已有许多文献在考虑可再生能源出力不确定性的优化调度策略方面做了深入研究。文献[8]考虑风电功率预测误差不确定性的电力系统储能容量优化配置的计算方法,并采用概率分布鲁棒联合机会约束规划模型进行描述。文献[9]综合考虑风电、负荷等不确定性因素及需求侧资源互动,提出了一种考虑源-荷的多场景优化系统旋转备用的方法,并提高了风电的消纳水平。文献[10]构建了考虑旋转备用的水火风互补发电优化调度模型,并分析了不同情景下对系统发电量和旋转备用的影响。随着高比例风电并网,仅靠常规水、火电机组参与调节和提供旋转备用已不能满足系统可再生能源消纳的要求。因此,需要考虑系统其他可调节资源与火电机组共同参与调节,以提高系统调度灵活性,如核电,作为未来能源的重要组成部分11,其最大调峰深度达70%的负荷跟踪能力12,可参与调峰。文献[13-14]对核电调峰模型进行了详细建模,并分析了不确定性对系统调度及备用的影响,有效提高了风电利用率和系统调度灵活性。

    但上述优化问题中忽略了系统运行存在的风险,这样会导致在决策值最优的情况下产生最大的预期利润或最小预期成本,但在某些不利情况下则以极低利润或高额成本为代价。为避免这种情况,可在问题中引入风险度量,对与利润或成本相关的可变性风险进行建模,衡量与利润或成本分配相关的风险。目前,在电力市场和金融优化的随机规划问题中,最常见的风险度量方法有:方差(Variance)、预期短缺ES(Expected shortage)、风险值VaR(Value-at-Risk)、条件风险价值CVaR(Conditional Value-at-Risk)等15。其中,CVaR除了满足一致性风险度量平移不变性、次可加性、正齐次性、单调性16的要求外,还可以用线性公式表示,因此在电力市场的相关问题中得到广泛运用。文献[17-18]利用CVaR对系统不确定所带来的风险进行评估,提出了确定市场最优备用需求的成本-CVaR模型,但并未考虑系统的发电优化调度过程。在此基础上,文献[19-21]利用CVaR度量风电不确定性给系统带来的风险损失,并对发电优化调度过程进行了分析讨论。此外,不少文献将CVaR作为风险约束进行考虑,文献[22-24]计及系统不确定性,构建风险约束的随机优化模型,将CVaR作为风险约束条件,对系统面临的风险进行约束,进而提高系统运行经济性和可靠性。

    上述研究在促进可再生能源消纳,提高系统运行经济性上具有一定的理论研究和应用价值,但仍存在可进一步研究和探讨的问题:(1)多从配网、微网角度分析,或主网电源类型较单一,缺少有关策略对含高比例可再生能源的多类型电源系统风险管控模型进行研究和分析;(2)鲜有研究分析核电参与调峰对含高比例可再生能源的系统调度策略、灵活性和风险损失的影响。

    综上,本文建立了计及火、水、风、核、储多类型电源机组组合和备用决策的调度模型,该模型通过多场景分析法应对风电出力和负荷预测的不确定性,并计及不确定性给系统带来的风险损失,采用条件风险价值(CVaR)构造系统在进行策略调控时的风险决策模型。算例结果验证了所提模型和方法的合理性和有效性。

    • 在高渗透率可再生能源背景下,可中断负荷具备快速响应的能力,在系统备用不足时,可通过负荷中断减少容量缺额,相当于提供瞬时的上调备用。电池储能系统具备响应速度快、时间短,实时调节能力强等特点,能提供快速的旋转备用。在满足储能系统相关约束的前提下,通过低储高发,提供上/下调备用。核电带基荷运行时,常规机组需要承担较大的调峰压力,部分常规机组需要频繁启停以满足系统的调峰需求;核电参与调峰时,常规机组仅需降负荷运行即可满足系统调峰需求,能有效缓解系统调峰压力。针对系统各类型资源特点进行协调优化调度,可充分挖掘系统的灵活性,缓解系统灵活性调节容量的紧缺。

      本文基于多场景分析方法建立了机组组合和备用决策联合优化的模型,能够考虑风电出力和负荷波动的随机性对电网调度运行的影响。通过风电随机出力和负荷预测误差的预想场景和误差场景分析系统发电和备用需求,综合考虑系统机组组合、备用决策和风险损失之间的平衡,得到兼顾经济性和可靠性的最优调度计划。其中,为应对高比例可再生能源出力和负荷波动出现系统供需难以平衡的现象,误差场景考虑多类型电源提供备用、可再生能源削减出力及可中断负荷等调节措施,并在预想场景中考虑核电参与调峰来进一步缓解高比例可再生能系统电源备用紧缺和降低弃风和失负荷损失。

      该优化过程包含两类决策过程:(1)通过预想场景风电出力与负荷值确定各类型电源出力、启停计划及核电调峰计划;(2)综合考虑误差场景不确定性带来的调度风险,通过多场景分析反映各类型电源在各种运行场景下的调节能力,确定各类型电源的备用需求及备用经济分配。在基于预想场景和误差场景的优化决策过程中,预想场景与误差场景的决策量相互影响,通过预想场景各电源的出力计划与启停计划影响误差场景的备用计划。反之,误差场景的备用计划同样影响预想场景各电源的出力计划和启停计划。

    • 基于上述优化思路的调度模型框架如图1所示。首先,基于风电出力和负荷预测预想场景,采用场景法生成ω个误差场景。其中场景法包括场景生成和场景削减两个部分。在场景生成部分,假设风电出力、负荷预测误差均服从正态分布25,采用拉丁超立方采样法26模拟生成大量随机场景,并采用Cholesky分解法降低各随机变量之间的相依性;在场景削减部分,由于场景个数过多会增加计算量,降低计算效率,而场景个数过少又不具代表性,影响最后的模拟精度,为了同时兼顾求解复杂程度和结果精确性,本文采用以概率距离27-28为基础的场景削减方法对生成场景进行削减。

      图  1  系统调度框架

      Figure 1.  Frame of system dispatching

      其次,以系统总调度成本(预想场景调度成本和误差场景调度期望成本之和)与CVaR值之和最小为目标,在满足预想场景、误差场景、CVaR约束的前提下,优化求解得到满足系统风险水平的最优调度计划。

    • 本文的计算框图如图2所示。首先根据预想场景负荷、风电出力求解线性规划问题得到机组组合序列;然后考虑误差场景不确定性,确定各类型电源的备用需求及备用经济分配;最后判断其是否满足误差场景约束及CVaR约束,依此往复直到所有约束均满足为止。

      图  2  计算流程图

      Figure 2.  Flow chart of the calculation

    • 忽略风电运行成本,本文提出的考虑条件风险价值的机组组合和多场景备用决策联合优化调度模型的目标函数包含两部分,如式(1)所示。

      min F=(1-β)×EC+β×CVaR ((1))
      EC=CBAS+CERR ((2))

      第一部分为系统总调度成本EC,包括预想场景下的调度成本CBAS和误差场景下的期望调度成本CERR,其计算模型分别详见第2.1和2.2节;第二部分为条件风险价值CVaR,其计算模型详见第2.3节;β为风险系数,取值范围为,其大小反映了在进行发电调度决策时对风险的厌恶水平。当β为0时,问题为风险中性模型,即决策者只关注成本的期望值,而忽略了系统不确定性造成成本波动的风险。β值越大,则表明决策者是风险厌恶型。

    • 预想场景下的机组机组组合目标函数如式(3)所示:

      CBAS=t=1Tg=1NG(CgPgtBAS+CgtSU+CgtSD)+t=1Th=1NHChPhtBAS+t=1Tn=1NN[CnPntBAS+CnR(Pnmax-PntBAS)]+t=1Te=1NE(CcePcetBAS+CdePdetBAS)  ((3))

      式中:第一项为火电机组调度成本,包括机组运行成本、启动成本和停机成本;第二项为水电机组调度成本;第三项为核电机组调度成本,包括运行成本和调峰成本;第四项为储能调度成本,包括充电和放电成本。

      NG ——火电机组个数;

      NH ——水电机组个数;

      NN ——核电机组个数;

      NE ——储能机组个数;

      T ——调度周期,每天考虑24时段,则T为24;

      PgtBAS——火电机组g在时段t的机组出力;

      PhtBAS——水电机组h在时段t的机组出力;

      PntBAS——核电机组n在时段t的机组出力;

      PcetBAS——储能c在时段t的充电功率;

      PdetBAS——储能c在时段t的充电和放电功率;

      Pnmax ——核电机组n的额定功率;

      CgtSU ——火电机组g在时段t的启动成本;

      CgtSD ——火电机组g在时段t的停机成本;

      Cg ——火电机组g的边际发电成本报价;

      Ch ——水电机组h的边际发电成本报价;

      Cn ——核电机组n的边际发电成本报价;

      CnR ——核电机组n的调峰成本报价;

      Cce ——储能c在时段t的充电成本报价;

      Cde ——储能c在时段t的放电成本报价。

    • gGmPgtBAS+hHmPhtBAS+nNmPntBAS+eEm(PdetBAS-PcetBAS)+wNwPwtBAS-jJmLjtBAS-rΛmft(m,r)=0 ((4))
      ft(m,r)=B(m,r)(δmt-δrt) ((5))

      式中:

      gGm ——节点m上火电机组g

      hHm ——节点m上水电机组h

      nNm ——节点m上核电机组n

      eEm ——节点m上储能e

      wNw ——节点m上风电机组w

      jJm ——节点m上负荷j

      rΛm ——和节点m相连的节点r

      PwtBAS ——风电w在时段t的预测出力;

      LjtBAS ——负荷j在时段t的预测值;

      ft(m,r) ——时段t从节点m到节点r的线路潮流;

      B(m,r)——节点mr之间的电纳;

      δmt ——时段t节点m的电压相角;

      δrt ——时段t节点r的电压相角。

    • 火电机组的运行约束主要包括:

      P̲gtPgtBASP¯gtP¯gtPgmax(ugt-zg,t+1)+Sgdnzg,t+1P¯gtPg,t-1BAS+Rgupug,t-1+SgupygtP̲gtPgminugtP̲gtPg,t-1BAS-Rgdnug,t-Sgdnzgt ((6))
      ygt-zgt=ugt-ug,t-1ygt-zgt1ygt,zgt,ugt{0,1} ((7))
      t=1ζg[1-ug(t)]=0,ζg=min[T,(UTg-Ug0)ugt,t=0]t=kk+UTg-1ugtUTgygk,k=ζg+1T-UTg+1t=kT[ugt-ygt]0,k=T-UTg+2T ((8))
      t=1ξgugt=0,ξg=min{T,(DTg-Sg0)(1-ug,t=0)}t=kk+DTg-1(1-ugt)DTgzgk,k=ξg+1T-DTg+1t=kT(1-ugt-zgt)0,k=T-DTg+2T ((9))
      CgtSU=λgSUygtCgtSD=λgSDzgt ((10))

      其中,式(6)为火电机组出力和上/下爬坡约束,式(7)为火电机组启停0-1变量约束,式(8)~(9)为火电机组最小连续开机/停机时间约束,式(10)为火电机组启停机成本约束。

      式中:

      P̲gt/P¯gt ——火电机组出力上/下限;

      Pgmax/Pgmin ——火电机组的最大/最小出力;

      Rgup/Rgdn ——火电机组的向上/向下爬坡速率;

      Sgup/Sgdn ——火电机组的启动/停机爬坡速率;

      T ——调度周期;

      Ug0/Sg0 ——机组g在调度计划开始时(0小时结束)已经在线运行/离线的时间;

      UTg/DTg——机组的最小运行/停机时间;

      ugt ——0/1变量,如果机组g在时段t运行,则为1,否则为零;

      ygt ——启动0/1变量,如果机组g在时段t开始时启动,则等于1;

      zgt ——停机0/1变量,如果机组g在时段t开始时停机,则等于1;

      ug,t=0 ——初始运行状态,在线则为1,否则为0;

      λgSU/λgSD ——机组g的启动/停机成本。

    • 水电机组运行约束与火电机组类似,除要考虑机组出力约束、上/下爬坡约束外(参见火电机组运行约束),还要考虑水电转换约束、水量约束:

      PhtDA=0.009  81ηhQhtHht ((11))
      t=1TQhtdtQhmax ((12))

      式(11)为水电转换约束,ηh为水电机组h的水能转换效率,QhtHht分别为水电机组m在时段t的发电耗水量、水库水头高度;式(12)为水量约束,Qhmax为水电机组h的最大水库容量。

    • 考虑核电机组的运行特性,有2种运行模式: A模式和G模式。A模式下,机组带基荷,维持满功率运行29;G模式下,机组存在满功率和低功率两种运行方式13-14,可表示为:

      PntBAS=hntPnmax-i=13lntiPn,li+k=19sntkPn,sk ((13))
      hnt+i=13lnti+k=19sntk=1 ((14))
      hnxhnt-hn,t-1,  x=t,t+1,,t+TnH-1lnxlnt-ln,t-1,  x=t,t+1,,t+TnL-1 ((15))

      其中,式(13)为核电出力约束,式(14)为功率水平0/1变量约束,式(15)为机组满/低功率持续运行时间约束。

      式中:

      hnt ——核电机组n时段t的满功率运行标志,为0/1变量;

      lnti ——核电机组n时段t的低功率运行标志,为0/1变量;

      sntk ——核电机组n时段t的升/降功率运行标志,为0/1变量;

      Pnmax ——核电机组n的满功率水平;

      Pn,li ——不同调峰深度lnti对应的低功率水平,lnti对应3种调峰深度;

      Pn,sk ——运行标志sntk对应的过渡功率水平分别30%、50%、70%;

      TnH/TnL——机组满功率/低功率最小持续运行时间。

      另外,限于篇幅,升/降功率时间耦合约束未列出,可参照文献[13]。

    • ucetPceminPcetBASucetPcemaxudetPdeminPdetBASudetPdemaxucet+udet1SOCet=SOCe,t-1+(PcetBASηce-PdetBAS/ηde)ΔtSOCeminSOCetSOCemax ((16))

      式中:

      Pcemax ——储能e时段t的最大充电功率;

      Pcemin ——储能e时段t的最小充电功率;

      Pdemax ——储能e时段t的最大放电功率;

      Pdemin ——储能e时段t的最小放电功率;

      ucet ——充电0/1变量,充电则为1;

      udet ——放电0/1变量,放电则为1;

      ηce ——储能e充电效率;

      ηde ——储能e放电效率;

      SOCemax ——储能e的最大容量。

      SOCemin ——储能e的最小容量。

    • 0PwtDAPwmax ((17))

      式中:

      Pwmax——风电机组w的最大出力。

    • -fmax(m,r)ft(m,r)fmax(m,r) ((18))

      式中:

      fmax(m,r)——线路m-r的潮流传输极限。

    • 以误差场景模拟系统不确定性,考虑火电、水电和储能提供备用、风电削减出力、可中断负荷以应对系统不确定性。

    • 考虑多种不确定性的多场景备用决策模型目标函数如下式所示:

      CERR=ω=1NΩπω[t=1Tg=1NG(CgRURUgtω+CgRDRDgtω)+t=1Th=1NH(ChRURUhtω+ChRDRDhtω)+t=1Te=1NE(CeRURUetω+CeRDRDetω)+t=1Tw=1NWVwspillSwtω+t=1Tj=1NLVjLOLLjtωshed] ((19))

      式中:

      包含五项:第一项、第二项、第三项分别为火电、水电、储能的备用成本;第四项为弃风、失负荷损失;

      NΩ    ——场景集;

      NW    ——风电场个数;

      NL        ——负荷节点数;

      πω    ——场景ω的概率,下标ω表示误差场景编号,下同;

      RUgtω/RDgtω——火电机组g在时段t提供的上/下旋转备用;

      RUhtω/RDhtω——水电机组h在时段t提供的上/下旋转备用;

      RUetω/RDetω ——储能e在时段t提供的上/下旋转备用;

      CgRU/CgRD  ——火电机组g的上/下旋转备用报价;

      ChRU/ChRD  ——水电机组h的上/下旋转备用报价;

      CeRU/CeRD      ——储能e的上/下旋转备用报价;

      VjLOL          ——单位失负荷价值;

      Vwspill           ——单位弃风电量的惩罚因子;

      Ljtωshed            ——负荷j在时段t的中断功率;

      Swtω            ——风电场w在时段t的弃风功率。

    • gGmPgtω+hHmPhtω+nNmPntDA+eEm(Pdetω-Pcetω)+wNw(Pwtω-Swtω)-jJm(Ljtω-Ljtωshed)-rΛmftω(m,r)=0Pgtω=PgtDA+RUgtω-RDgtωPhtω=PhtDA+RUhtω-RDhtωPdetω=PdetDA+RUetωPcetω=PcetDA+RDetωftω(m,r)=B(m,r)(δmtω-δrtω) ((20))

      式中:

      Pgtω        ——火电机组g在时段t的机组出力;

      Phtω      ——水电机组h在时段t的机组出力;

      Pcetω/Pdetω    ——储能e在时段t的充电/放电功率;

      Pwtω     ——风电机组w在时段t的机组出力;

      Ljtω    ——在时段t节点j处的负荷值;

      ftω(m,r) ——时段t从节点m到节点r的线路潮流;

      B(m,r)   ——节点mr之间的电纳;

      δmtωδrtω   ——时段t节点mr的电压相角。

    • 机组出力约束包括火电、水电机组出力约束,储能充放电功率约束:

      ugtPgminPgtωugtPgmaxuhtPhminPhtωuhtPhmaxucetPceminPdetωucetPcemaxudetPdeminPcetωudetPdemax ((21))
    • 机组爬坡约束包括火电、水电机组爬坡约束:

      -RgdnPgtω-Pg,t-1,ωRgup-RhdnPhtω-Ph,t-1,ωRhup ((22))

      式中:

      Rhup——水电机组的向上爬坡速率;

      Rhdn——水电机组的向下爬坡速率。

    • 备用容量约束包括火电、水电机组和储能的备用约束:

      0RUgtωRupg0RDgtωRdng0RUhtωRuph0RDhtωRdnhRUetω0,RDetω0 ((23))

      式中:

      Rupg——火电机组的最大向上备用容量;

      Rdng——火电机组的最大向下备用容量;

      Ruph——水电机组的最大向上备用容量;

      Rdnh——水电机组的最大向下备用容量。

    • 0SwtωPwtω ((24))
    • 0LjtωshedLjtω ((25))
    • -fmax(m,r)ftω(m,r)fmax(m,r) ((26))
    • CVaR(条件风险价值)和VaR(风险价值)是常用的风险度量工具,但VaR存在不满足一致性公理,缺乏次可加性,尾部风险信息识别不准确,计算复杂等缺点30-31。CVaR在VaR的基础上考虑了超额风险损失的平均水平,除满足一致性风险度量外,还可以用线性公式表示,是电力市场相关问题中最常用的风险度量。基于CVaR的损失风险可转化为下面的线性优化问题求解32

      CVaR=miny+11-αs=1Nspszss.t.   zω[CDA+t=1NTg=1NG(CgRURUgtω+CgRDRDgtω)          +t=1NTh=1NHChRURUhtω+ChRDRDhtω          +t=1NTe=1NE(CeRURUetω+CeRDRDetω)          +t=1NTj=1NLVjLOLLjtωshed+t=1NTw=1NWVwspillSwtω]-y      zω0,ω=1,2,...,NΩ ((27))

      式中:

      yzω——辅助变量。

      当优化得到CVaR时,y即为所有场景下置信度为α时的经济损失风险价值。1-α为极端风险尾部风险的分布范围。

    • 本文模型目标函数是一次函数,约束条件是线性约束,且包含整数变量,故该优化模型为混合整数线性规划模型(MILP)。在GMAS 24.4环境下编写计算程序,通过调用GAMS混合混合整数线性规划(MILP)求解器CPLEX进行求解。

    • 本文采用改进IEEE 24节点系统作算例分析。该系统包含34条输电线路,2台核电机组,1个容量为320 MWh储能电站,10台非核电机组,其中8号机组为水电机组,其余为火电机组。4个风电场额定容量均为300 MW,并假设风电场预测出力均相同。非核电机组上调、下调备用为边际成本报价的2倍,停机成本为机组启动成本的0.5倍。储能电站数据取自文献[33],上调、下调备用单位成本取15.33 $/MWh19。核电机组参数取自文献[13],边际成本报价均为9.98 $/MWh,核电调峰单位成本均为13.76 $/MWh34。调度时段T取24小时。失负荷价值VjLOL取1 000 $/MWh,弃风电量惩罚因子Vwspill取200 $/MWh。置信水平α取95%。各电源、电网、负荷、储能的运行参数详见附录A,如无特别说明,算例分析均采用4个风电场同时并网,风险系数β均取1。

    • 根据风电和负荷预测场景,采用拉丁超立方采样法生成1 000个风电场景和1 000个负荷场景,并采用以概率距离为基础的场景缩减方法将风电场景缩减到10个场景,负荷场景缩减到5个场景。风电出力与负荷场景分别如图3图4所示。可以看出,在误差场景阶段的调整过程中,风电出力的波动较大,且具有明显的反调峰特性,而负荷的波动性较小。

      图  3  风电场景

      Figure 3.  Wind power scenarios

      图  4  负荷场景

      Figure 4.  Load scenarios

      本文只考虑风电和负荷不确定性,总共包含50个场景,场景概率为相应风电出力场景概率与负荷场景概率的乘积。风电出力场景概率和负荷场景概率如表1所示。

      表 1  风电和负荷场景概率

      Table 1.  Probability of wind power and load scenarios

      风电场景S概率ps负荷场景S概率ps
      S10.109S10.270
      S20.198S20.160
      S30.055S30.273
      S40.059S40.114
      S50.047S50.183
      S60.026
      S70.207
      S80.038
      S90.109
      S100.152
    • 为比较多类型电源协调调度对系统运行技术经济性的影响,本文设计了4种不同调度方案作对比分析。具体方案如表2所示。通过优化,得到4种不同调度方案下运行的各项指标,如表3所示。

      表 2  不同调度方案

      Table 2.  Dispatch of different schemes

      方案火电水电可中断负荷核电储能
      方案1+++
      方案2++++
      方案3++++
      方案4+++++

      表 3  不同调度方案下的系统运行指标

      Table 3.  System performance indices of different schemes

      运行指标方案1方案2方案3方案4
      总发电调度成本811 327.3783 272.3781 827.0778 702.2
      预想场景调度成本743 327.1761 252.3745 026.7757 102.1
      误差场景期望调度成本68 000.222 020.036 800.421 600.1
      备用容量期望成本51 160.920 682.633 544.621 600.1
      弃风损失期望16 466.01 337.43 255.80.0
      失负荷损失期望373.30.00.00.0
      弃风、失负荷损失期望16 839.31 337.43 255.80.0
      CVaR813 409.4783 334.2781 857.2778 724.3

      表3可以看出,与方案1相比,方案2增加核电机组参与调峰弥补了谷荷时段下调备用的不足,弃风损失减少了92%,失负荷损失减为0;方案3增加了储能提供上调、下调旋转备用,弃风损失减少了80%,失负荷损失减为0;方案4(本文方案)核电、储能参与调节,有效缓解了系统备用容量的紧缺,弃风、失负荷损失均降为0。

      方案1和方案4各时段机组出力分别如图5图6所示。可以看到,与方案1相比,在负荷低谷、风电多发时段(1:00-7:00),方案4中的核电参与调峰,相当于提供下调备用,增加了调峰费用和非核电机组出力,储能则通过低储高发提供上调、下调旋转备用;总体上,核电和储能参与调峰,使预想场景调度成本增加了2.4%,但备用容量期望成本、弃风和失负荷损失均大幅降低,即误差场景期望调度成本大幅降低(降低67.6%),总调度成本降低了3.5%;方案4的总发电调度成本和CVaR均最小,说明多类型电源调度能够大幅提高系统运行的经济性和降低系统面临的潜在风险损失,满足具有较高风险厌恶程度的决策者的调度需求。

      图  5  方案1各时段机组出力

      Figure 5.  Unit output of scheme 1

      图  6  方案4各时段机组出力

      Figure 6.  Unit output of scheme 4

      为进一步说明不同类型电源调度对系统运行的影响,方案1和方案4各时段非核电机组启停分别如表4表5所示,各时段非核机组的启动台数如图7所示。结合图4图7表4表5可以看出,两种方案的机组启动台数的变化趋势与负荷的变化趋势基本一致,在负荷的低谷期(1:00-7:00)启动台数减少,在负荷高峰期(10:00-17:00)启动台数增多。图7中,方案4曲线相较于方案1平缓,这是因为方案1中,核电带基荷运行、储能不参与调峰,大规模风电并网和负荷波动大幅增加了系统的不确定性,非核机组需要频繁启停以及时对风电和负荷的波动做出反应,进而大幅增加了非核机组的启动频次、启动成本和备用成本;相较于方案1,方案4核电和储能参与调峰,降低了非核电旋转备用预留的容量,避免非核电机组因预留容量的频繁启停,节约了旋转备容预留成本(降低57.8%)和机组启停成本(降低64.7%),但随着投入机组的增多,系统的运行成本会有所增加(预想场景调度成本增加1.9%),但总的调度成本降低了4.0%,结合表3得出上述结论。

      表 4  方案1各时段机组启动计划及启动成本

      Table 4.  Unit commitment planning and cost of schemes 1

      机组启停状态机组编号时间段成本/$
      启动1t81 430.4
      3t81 725.0
      4t103 056.7
      6t8312.0
      9t11624.0
      10t102 298.0
      停机1t4715.2
      3t16862.5
      4t221 528.35
      6t18156.0
      7t24156.0
      9t3312.0
      10t11 149.0

      表 5  方案4各时段机组启停计划及启停成本

      Table 5.  Unit commitment planning and cost of schemes 4

      机组启停状态机组编号时间段成本/$
      启动4t103 056.7
      6t9312.0
      停机4t221 528.35
      6t18156.0

      图  7  方案1和方案4各时段非核机组的启动台数

      Figure 7.  Number of non-nuclear generators in operation in each period of scheme 1 and scheme 4

      此外,结合表3~表5注意到,相较方案1,方案4 CVaR降低了4.3%,说明系统决策者通过减少机组频繁启停来规避系统不确定性带来的成本波动风险,显著提高了系统运行的经济性和可靠性。

    • 针对不同风电渗透率,方案1和方案4的计算结果如表6。结果表明,随着风电渗透率的增大,机组出力降低,预想场景调度成本呈逐步下降趋势,为应对系统的不确定性,需预留足够的备用容量,备用容量期望成本呈逐步上升趋势,同时风电实际出力与预测出力的偏差增大,弃风和失负荷损失逐步增大,误差场景期望调度成本呈上升趋势。同时,两方案的总调度成本均呈下降趋势,这表明风电渗透率的增大,系统运行经济性得到改善。根据表6计算结果,可以得到不同风电渗透率下的CVaR和总调度成本的有效边界如图7所示。由表6图7注意到:

      表 6  不同风电渗透率下方案1和方案4计算结果

      Table 6.  Calculation results of scheme 1 and scheme 4 under different wind power penetration

      风电容量/MW渗透率/%方案预测场景调度成本/$备用容量期望成本/$弃风、失负荷损失期望/$
      3006.8方案1823 275.37 709.4100.7
      3006.8方案4822 871.95 242.80
      60013.6方案1798 662.39 723.40
      60013.6方案4796 863.39 395.80
      90020.3方案1773 252.818 031.81 523.8
      90020.3方案4768 538.319 214.02 327.8
      1 20027.1方案1743 327.151 160.916 839.3
      1 20027.1方案4757 102.121 600.10
      风电容量/MW渗透率/%方案误差场景期望调度成本/$总调度成本/$CVaR/$
      3006.8方案17 810.1831 085.5832 058.0
      3006.8方案45 242.8828 114.7828 165.1
      60013.6方案19 723.4808 385.7808 501.4
      60013.6方案49 395.8806 259.1806 360.7
      90020.3方案119 555.6792 808.4792 826.1
      90020.3方案421 541.8790 080.0790 080.0
      1 20027.1方案168 000.2811 327.3813 409.4
      1 20027.1方案421 600.1778 702.2778 724.3

      方案1中,总调度成本和CVaR随着风电渗透率的增加呈下降趋势,但在风电渗透率为27%,总调度成本和CVaR发生跃变,这是由于大规模风电并网,系统备用期望成本大幅增加和在负荷低谷、风电多发时段(1:00-7:00)下调备用不足,弃风和失负荷大幅增加所引起的。且风电渗透率越大,弃风和失负荷损失越严重。

      方案4中,渗透率为20.3%产生了弃风和失负荷现象,这是由于核电参与调峰的成本较弃风、失负荷损失高,核电未参与调峰以及机组和储能提供的下调备用不足所引起的;当随着风电渗透率进一步增加(27.1%),核电参与调峰,系统备用充足,弃风和失负荷损失降为0,更为有效促进风电消纳;总调度成本和CVaR均呈下降趋势,说明多类型电源调度能提高系统的经济性,降低系统的风险损失。且由图8表7方案4各渗透率较6.8%渗透率总调度成本和CVaR下降幅度注意到,风电渗透率越高,下降趋势越明显。进一步表明多类型电源调度提高系统经济性和降低系统风险损失的效果随风电渗透率的增加愈加显著。

      图  8  不同风电渗透率下的总调度成本和CVaR的有效边界

      Figure 8.  Efficient frontier of Dispatching CVaR under different penetration of wind power

      表 7  方案4总调度成本和CVaR降幅

      Table 7.  Total dispatching cost and CVaR reduction of scheme 4

      渗透率/%总调度成本/%CVaR/%
      6.8
      13.62.642.63
      20.34.594.60
      27.15.975.97
    • 表8为不同风险系数的计算结果,可以看出,随着风险系数的逐渐增大,系统总调度成本呈上升趋势,变化范围从7.58×105 $到7.79×105 $,CVaR呈下降趋势,变化范围从7.99×105 $到7.79×105 $,CVaR数值减少2.6%导致系统总调度成本增加2.7%。根据表8的计算结果,可得到不同风险系数下总调度成本和CVaR的有效边界,如图9所示。可以看到,风险厌恶程度较低,即风险系数β较小时,随着β增大,系统总调度成本增加较平缓;当风险厌恶程度较高,即β较大,随着β增大,系统总调度成本快速增加。

      表 8  不同风险系数下的计算结果

      Table 8.  Calculation results under different risk parameters

      风险系数β预想场景调度成本误差场景期望调度成本总发电调度成本CVaR
      0737 365.720 960.1758 325.8799 213.3
      0.1739 283.019 192.8758 475.7795 472.6
      0.2741 279.517 629.8758 909.4793 106.0
      0.3743 618.316 072.8759 691.1790 768.8
      0.4748 617.113 694.4762 311.5785 779.5
      0.5754 022.312 337.9766 360.2780 875.7
      0.6756 267.512 139.0768 406.5779 036.2
      0.7756 858.512 006.1768 864.6778 800.1
      0.8757 065.011 982.9769 047.8778 738.1
      0.9757 102.112 014.4769 116.5778 724.3
      1757 102.121 600.1778 702.2778 724.3

      图  9  不同风险系数下系统总调度成本和CVaR的有效边界

      Figure 9.  Efficient frontier of dispatching cost vs CVaR under different risk parameters

    • 为进一步验证本文所提模型和方法的正确性和有效性,在某省级电网系统中进行算例分析。该系统等值后的电网数据包含89台发电机、1 151个母线节点以及1 485条输电线路。其中,89台发电机中包括43个水电机组、41个火电机组、2个风电场,2个核电厂,1个储能电站;1 151个母线节点包括202个负荷节点,其中,202个负荷节点中包含6个等值联络线负荷节点(10回联络线等值);1 485条输电线路中220 kV及以上电压等级线路有621条(已并联处理,去除变压器支路)。根据某省电网96点负荷曲线及风电场预测出力,采用拉丁超立方及场景缩减法得到风电出力与负荷场景分别如图10图11图12所示。

      图  10  负荷场景

      Figure 10.  Load scenarios

      图  11  风电1场景

      Figure 11.  Scenarios of No.1 wind farm

      图  12  风电2场景

      Figure 12.  Scenarios of No.2 wind farm

      表9为方案1和方案4风险系数为0.5时系统调度方案的运行指标。可以看到,由于系统电源较多,且风电占比较小,系统备用充足,各方案不存在弃风失负荷现象,但由于方案4中核电和储能参与调峰,降低了非核电旋转备用预留的容量,避免非核电机组因预留容量的频繁启停,节约了旋转备容预留成本,其总发电调度成本和CVaR均较小于方案1,说明多类型电源调度能够有效提高系统运行的经济性和降低系统面临的潜在风险损失。表10为不同风险系数的计算结果,可以看出,随着风险系数的逐渐增大,系统总调度成本呈上升趋势。

      表 9  不同调度方案下的系统运行指标

      Table 9.  System performance indices of different schemes

      运行指标方案1方案4
      总发电调度成本599 682 595598 848 191
      预想场景调度成本592 098 821592 006 420
      误差场景期望调度成本7 583 7746 841 771
      备用容量期望成本7 583 7746 841 771
      弃风损失期望0.00.0
      失负荷损失期望0.00.0
      弃风、失负荷损失期望0.00.0
      CVaR599 974 173599 038 682

      表 10  不同风险系数下的计算结果

      Table 10.  Calculation results under different risk parameters

      风险系数β预想场景调度成本误差场景期望调度成本总发电调度成本CVaR
      0591 506 8547 087 759598 594 613599 782 842
      0.1591 799 1656 925 843598 725 008599 548 623
      0.2591 874 1696 921 227598 795 396599 448 130
      0.3591 911 4556 902 821598 814 276599 116 590
      0.4591 952 0636 876 704598 828 767599 069 834
      0.5592 006 4206 841 771598 848 191599 038 682
      0.6592 056 1746 826 255598 882 429599 015 802
      0.7592 087 3246 802 164598 889 488599 010 729
      0.8592 139 4616 784 591598 924 052599 008 430
      0.9592 237 6616 764 004599 001 665598 856 317
      1592 237 9366 892 163599 130 099598 853 578
    • 大规模可再生能源并网以及负荷波动性显著增加了系统的不确定。此外,火、水、风、核、储能等多类型电源并网,增加了系统的运行复杂度,有必要考虑多类型电源互补优化模型的构建。基于此,本文提出了计及多类型电源调度和系统不确定性的机组组合和备用决策的联合优化模型和方法,主要结论如下:

      1)根据火、水、风、核、储能等多类型电源系统的特点,从电源侧、电网侧、负荷侧和储能四部分入手构建了含多类型电源的源-网-荷-储的互补协同优化模型。

      2)充分挖掘多类型电源系统的灵活性资源,提出了考虑可中断负荷和不同类型电源的机组组合和备用决策的联合优化模型。算例表明,随着风电渗透率提高,不同类型电源的机组组合和备用决策在提高系统运行经济性、可靠性和降低弃风电量方面的作用越明显,能够满足大规模可再生能源并网的调度决策要求。

      3)在多类型电源发电调度的基础上引入条件风险价值(CVaR),利用CVaR度量系统不确定性带来的风险损失,能够获得考虑系统运行风险的最优调度计划。通过不同风险厌恶程度下的算例仿真,验证了所提模型的合理性和有效性。同时,探讨了系统发电调度费用和风险水平之间的关系,为系统不同风险厌恶程度决策者提供有价值的参考信息。

  • 李建钊,谢敏,李舒佳等.考虑CVaR的机组组合和多场景备用决策联合优化[J].南方能源建设,2021,08(04):50-65.
  • 参考文献 (34)

    目录

      /

      返回文章
      返回