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

留言板

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

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

基于连续隐马尔可夫模型的风水火联合低碳检修优化

何知纯 谢敏 黄莹 李弋升 张世平

何知纯, 谢敏, 黄莹, 李弋升, 张世平. 基于连续隐马尔可夫模型的风水火联合低碳检修优化[J]. 南方能源建设, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
引用本文: 何知纯, 谢敏, 黄莹, 李弋升, 张世平. 基于连续隐马尔可夫模型的风水火联合低碳检修优化[J]. 南方能源建设, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
HE Zhichun, XIE Min, HUANG Ying, LI Yisheng, ZHANG Shiping. Wind Power, Hydropower and Thermal Power Combined Low-Carbon Maintenance Optimization Based on Continuous Hidden Markov Model[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
Citation: HE Zhichun, XIE Min, HUANG Ying, LI Yisheng, ZHANG Shiping. Wind Power, Hydropower and Thermal Power Combined Low-Carbon Maintenance Optimization Based on Continuous Hidden Markov Model[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005

基于连续隐马尔可夫模型的风水火联合低碳检修优化

doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
基金项目: 广东省基础与应用基础研究基金面上项目“虚拟电厂动态聚合演化理论与计算方法研究”(2021A1515012245);中央高校基本科研业务费重大产学研合作扶持专项“主动配电网多网多体协同优化决策技术研究”(x2dlD2201280)
详细信息
    作者简介:

    何知纯,1998-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1329029641@qq.com

    谢敏,1978-,女,博士,副教授,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)minxie@scut.edu.cn

    黄莹,1999-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1419928558@qq.com

    李弋升,1998-,男,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)907353249@qq.com

    张世平,1998-,男,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1632522154@qq.com

    通讯作者:

    何知纯,1998-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1329029641@qq.com

  • 中图分类号: TK89;TM315

Wind Power, Hydropower and Thermal Power Combined Low-Carbon Maintenance Optimization Based on Continuous Hidden Markov Model

  • 摘要:   目的  新型电力系统背景下,风机的低碳检修与常规机组的协同检修等问题亟待解决。文章兼顾多属性气象因子影响和低碳性、经济性需求,建立基于连续隐马尔可夫的风水火联合低碳检修优化模型。  方法  首先,以降雨量、风速、雷电危险度为观测序列,以检修容量为隐状态序列,利用连续隐马尔可夫(Continuous Hidden Markov Model, CHMM)过程实现对风电场检修容量的动态跟踪。然后,以最优检修容量为决策依据、以总成本最小为优化目标,统筹考虑检修约束、系统控制约束等,构建风水火联合低碳检修优化模型。最后,以IEEE30节点系统进行算例展开研究。  结果  结果表明,所提模型具有更显著的经济效益和低碳特性。  结论  文章所做研究对风机的运行维护具有较高的理论价值,工程适用性较强。
  • 图  1  风电场日检修容量值的确定步骤

    Fig.  1  Steps to determine the daily maintenance capacity of wind farm

    图  2  基于CHMM的风电场检修容量跟踪模型的流程

    Fig.  2  Flow chart of wind farm maintenance capacity tracking model based on CHMM

    图  3  风水火联合低碳检修优化模型的求解步骤

    Fig.  3  Solving steps of the optimization model for wind power, hydropower and thermal power combined low-carbon maintenance

    图  4  最优观测值-隐状态分布

    Fig.  4  Hidden state distribution of optimal observation value

    图  5  IEEE30节点系统拓扑图

    Fig.  5  IEEE 30-node system topology diagram

    图  6  风电场预测功率曲线

    Fig.  6  Forecast power of wind farm

    图  7  风电场日负荷预测曲线

    Fig.  7  Daily load forecasting curve of wind farm

    图  8  水火电/风电场单位容量的检修费用曲线

    Fig.  8  Maintenance cost curve of unit capacity of hydropower station, thermal power station/wind farm

    图  9  规划期内各机组检修计划

    Fig.  9  Maintenance plan of each WTGS in the planning period

    图  10  各类型机组出力变化

    Fig.  10  Output change of each type of WTGS

    图  11  各风电场的出力曲线

    Fig.  11  Output curve of each wind farm

    图  12  优化过程中各成本变化曲线

    Fig.  12  Cost change curves in the optimization process

    图  13  两模型求解检修区间对比图

    Fig.  13  Comparison diagram of maintenance interval solved by two models

    表  1  各隐状态物理含义

    Tab.  1.   Physical meaning of each hidden state

    状态检修容量物理含义
    1$(0.00 \sim 0.50){P_{\max ,w} }$少部分机组需要检修
    2$(0.50 \sim 0.75){P_{\max ,w} }$
    3$(0.75 \sim 0.90){P_{\max ,w} }$大部分机组需要检修
    4$(0.90 \sim 0.95){P_{\max ,w} }$
    5$(0.95 \sim 1.00){P_{\max ,w} }$几乎所有的机组需要检修
    下载: 导出CSV

    表  2  观测序列数据集形式

    Tab.  2.   Observation sequence data set form

    日期降雨量/mm风速/(m·s−1)雷电危险度
    D38151.680.3
    D391679.710.5
    ……………………
    下载: 导出CSV

    表  3  初始状态概率矩阵${\boldsymbol{\pi}} $

    Tab.  3.   Initial state probability matrix ${\boldsymbol{\pi}} $

    初始概率 ${s_1}$${s_2}$${s_3}$${s_4}$${s_5}$
    ${\boldsymbol{\pi}}$10000
    下载: 导出CSV

    表  4  状态转移概率矩阵${\boldsymbol{A}}$

    Tab.  4.   State transition probability matrix ${\boldsymbol{A}}$

    隐状态 ${s_1}$${s_2}$${s_3}$${s_4}$${s_5}$
    ${s_1}$0.263 00.483 90.161 90.035 00.056 2
    ${s_2}$0.385 70.261 30.205 00.039 00.109 1
    ${s_3}$0.301 90.122 50.373 20.177 10.025 3
    ${s_4}$0.263 60.000 00.139 90.458 00.138 5
    ${s_5}$0.211 00.156 20.288 50.055 60.288 7
    下载: 导出CSV

    表  5  隐状态概率矩阵${\boldsymbol{B}}$

    Tab.  5.   Hidden state probability matrix ${\boldsymbol{B}}$

    日期s1s2s3s4s5
    D1 0.791 624 7 0.207 005 4 0.001 272 4 0.000 086 1 0.000 011 3
    D2 0.038 939 9 0.960 987 4 0.000 016 6 0.000 052 3 0.000 003 8
    D3 0.001 347 6 0.998 619 9 0.000 000 4 0.000 031 2 0.000 000 9
    D4 0.791 687 4 0.203 395 1 0.004 805 6 0.000 099 1 0.000 012 7
    D5 0.785 881 3 0.212 917 9 0.001 102 6 0.000 086 8 0.000 011 4
    …… …… …… …… …… ……
    下载: 导出CSV

    表  6  未来一月的最优隐状态序列

    Tab.  6.   Optimal hidden state sequence for the next month

    隐状态天数/d天数占比/%
    11240.0
    21033.3
    3516.7
    400.0
    5310.0
    下载: 导出CSV

    表  7  各机组特性参数和检修工期

    Tab.  7.   Characteristic parameters and maintenance duration of each WTGS

    机组$ {P_{i,\max }} $/
    MW
    $ {P_{i,\min }} $/
    MW
    c2/
    (元·MW−1)
    c1/
    (元·MW−1)
    c0/
    (元·MW−1)
    检修工期/
    d
    i1 80.0 0 0.140 0 14.00 0 10
    i2 80.0 0 0.121 7 12.17 0 7
    i3 50.0 0 4
    i4 55.0 0 0.058 0 22.60 0 5
    i5 30.0 0 0.173 9 20.86 0 3
    i6 40.0 0 0.173 9 20.86 0 4
    w1 15.0 0 3
    w2 22.5 0 3
    w3 30.0 0 3
    w4 30.0 0 3
    下载: 导出CSV

    表  8  规划周期内各风电场检修安排

    Tab.  8.   Maintenance arrangement of each wind farm within the planning period

    风电场开始时间/d持续时间/d结束时间/d检修容量/MW
    w1283306.0
    w22232513.5
    w32832016.5
    w473922.5
    下载: 导出CSV

    表  9  不同模型求解结果对比

    Tab.  9.   Comparison of solution results of different models

    参数优化模型对照模型对比
    总成本/元 2 628 384 2 629 746 −0.05%
    碳排放成本/元 268 731.3 269 106 −0.14%
    运行成本/元 1 881 674 1 892 588 −0.58%
    检修成本/元 477 979 468 052.1 2.12%
    水火电检修成本/元 450 653.2 447 661.8 0.67%
    风电检修成本/元 27 325.8 20 390.33 34.01%
    水火电出力/MW 5 893.229 5 901.448 −0.14%
    风电出力/MW 922.654 2 914.435 9 0.90%
    下载: 导出CSV
  • [1] 张智刚, 康重庆. 碳中和目标下构建新型电力系统的挑战与展望 [J]. 中国电机工程学报, 2022, 42(8): 2806-2818. DOI:  10.13334/j.0258-8013.pcsee.220467.

    ZHANG Z G, KANG C Q. Challenges and prospects for constructing the new-type power system towards a carbon neutrality future [J]. Proceedings of the CSEE, 2022, 42(8): 2806-2818. DOI:  10.13334/j.0258-8013.pcsee.220467.
    [2] 舒印彪, 张丽英, 张运洲, 等. 我国电力碳达峰、碳中和路径研究 [J]. 中国工程科学, 2021, 23(6): 1-14. DOI:  10.15302/J-SSCAE-2021.06.001.

    SHU Y B, ZHANG L Y, ZHANG Y Z, et al. Carbon peak and carbon neutrality path for China′s power industry [J]. Strategic study of CAE, 2021, 23(6): 1-14. DOI:  10.15302/J-SSCAE-2021.06.001.
    [3] 葛畅, 阎洁, 刘永前, 等. 海上风电场运行控制维护关键技术综述 [J]. 中国电机工程学报, 2022, 42(12): 4278-4291. DOI:  10.13334/j.0258-8013.pcsee.220637.

    GE C, YAN J, LIU Y Q, et al. Review of key technologies for operation control and maintenance of offshore wind farm [J]. Proceedings of the CSEE, 2022, 42(12): 4278-4291. DOI:  10.13334/j.0258-8013.pcsee.220637.
    [4] 官嫣嫣, 庄佳才. 基于风险管理的海上风电进度管理方法研究 [J]. 南方能源建设, 2022, 9(1): 34-39. DOI:  10.16516/j.gedi.issn2095-8676.2022.01.005.

    GUAN Y Y, ZHUANG J C. Research on schedule management method of offshore wind power based on risk management [J]. Southern energy construction, 2022, 9(1): 34-39. DOI:  10.16516/j.gedi.issn2095-8676.2022.01.005.
    [5] 田帅. 新能源风险投资全周期风险管理体系的构建 [J]. 南方能源建设, 2022, 9(1): 40-46. DOI:  10.16516/j.gedi.issn2095-8676.2022.01.006.

    TIAN S. Construction of the whole cycle risk management system of new energy venture capital [J]. Southern energy construction, 2022, 9(1): 40-46. DOI:  10.16516/j.gedi.issn2095-8676.2022.01.006.
    [6] 晁化伟. 计及恶劣天气影响的海上风电场可靠性评估与维修/维护资源优化 [D]. 重庆: 重庆大学, 2019. DOI: 10.27670/d.cnki.gcqdu.2019.002117.

    CHAO H W. Offshore wind farms reliability evaluation and repair/maintenance resources optimization considering severe weather conditions [D]. Chongqing: Chongqing University, 2019. DOI: 10.27670/d.cnki.gcqdu.2019.002117.
    [7] 马广璐. 风电机组机会维修策略研究 [D]. 北京: 华北电力大学(北京), 2021. DOI: 10.27140/d.cnki.ghbbu.2021.000562.

    MA G L. Research on opportunity maintenance strategy of wind turbine [D]. Beijing: North China Electric Power University (Beijing), 2021. DOI: 10.27140/d.cnki.ghbbu.2021.000562.
    [8] 黄玲玲, 曹家麟, 张开华, 等. 海上风电机组运行维护现状研究与展望 [J]. 中国电机工程学报, 2016, 36(3): 729-738. DOI:  10.13334/j.0258-8013.pcsee.2016.03.017.

    HUANG L L, CAO J L, ZHANG K H, et al. Status and prospects on operation and maintenance of offshore wind turbines [J]. Proceedings of the CSEE, 2016, 36(3): 729-738. DOI:  10.13334/j.0258-8013.pcsee.2016.03.017.
    [9] 符杨, 黄路遥, 刘璐洁, 等. 基于状态自适应评估的海上风电机组预防性维护策略 [J]. 电力自动化设备, 2022, 42(1): 1-9. DOI:  10.16081/j.epae.202110009.

    FU Y, HUANG L Y, LIU L J, et al. Preventive maintenance strategy for offshore wind turbine based on state adaptive assessment [J]. Electric power automation equipment, 2022, 42(1): 1-9. DOI:  10.16081/j.epae.202110009.
    [10] SI X S, LI T M, ZHANG Q, et al. An optimal condition-based replacement method for systems with observed degradation signals [J]. IEEE transactions on reliability, 2018, 67(3): 1281-1293. DOI:  10.1109/TR.2018.2830188.
    [11] 李锁, 黄玲玲, 刘阳, 等. 基于风电机组状态信息的海上风电场维护策略 [J]. 现代电力, 2022, 39(1): 26-35. DOI:  10.19725/j.cnki.1007-2322.2021.0004.

    LI S, HUANG L L, LIU Y, et al. An offshore wind farm maintenance strategy based on wind turbine condition information [J]. Modern electric power, 2022, 39(1): 26-35. DOI:  10.19725/j.cnki.1007-2322.2021.0004.
    [12] 黄海悦. 基于零件可靠度的海上风电机组机会维护策略 [D]. 北京: 华北电力大学(北京), 2018.

    HUANG H Y. Offshore wind turbine opportunity maintenance strategy based on component reliability [D]. Beijing: North China Electric Power University (Beijing), 2018.
    [13] 苏春, 胡照勇, 郑玉巧. 基于可用度约束的风力机单部件顺序维修优化 [J]. 东南大学学报(自然科学版), 2019, 49(1): 110-115. DOI:  10.3969/j.issn.1001-0505.2019.01.016.

    SU C, HU Z Y, ZHENG Y Q. Single part sequential maintenance optimization for wind turbines based on availability constraint [J]. Journal of southeast university (natural science edition), 2019, 49(1): 110-115. DOI:  10.3969/j.issn.1001-0505.2019.01.016.
    [14] 逯红霞, 张蕊萍, 董海鹰. 考虑故障相关性的风电机组维修策略 [J]. 可再生能源, 2020, 38(4): 477-483. DOI:  10.13941/j.cnki.21-1469/tk.2020.04.009.

    LU H X, ZHANG R P, DONG H Y. Considering maintenance strategy of wind turbines with fault correlation [J]. Renewable energy resources, 2020, 38(4): 477-483. DOI:  10.13941/j.cnki.21-1469/tk.2020.04.009.
    [15] 符杨, 许伟欣, 刘璐洁, 等. 考虑天气因素的海上风电机组预防性机会维护策略优化方法 [J]. 中国电机工程学报, 2018, 38(20): 5947-5956. DOI:  10.13334/j.0258-8013.pcsee.171695.

    FU Y, XU W X, LIU L J, et al. Optimization of preventive opportunistic maintenance strategy for offshore wind turbine considering weather conditions [J]. Proceedings of the CSEE, 2018, 38(20): 5947-5956. DOI:  10.13334/j.0258-8013.pcsee.171695.
    [16] 符杨, 杨凡, 刘璐洁, 等. 考虑部件相关性的海上风电机组预防性维护策略 [J]. 电网技术, 2019, 43(11): 4057-4063. DOI:  10.13335/j.1000-3673.pst.2019.0178.

    FU Y, YANG F, LIU L J, et al. Preventive maintenance strategy for offshore wind turbines considering component correlation [J]. Power system technology, 2019, 43(11): 4057-4063. DOI:  10.13335/j.1000-3673.pst.2019.0178.
    [17] IRAWAN C A, OUELHADJ D, JONES D, et al. Optimisation of maintenance routing and scheduling for offshore wind farms [J]. European journal of operational research, 2017, 256(1): 76-89. DOI:  10.1016/j.ejor.2016.05.059.
    [18] 王东. 多重相关性影响下的风机系统维护决策优化 [D]. 上海: 上海海洋大学, 2022. DOI: 10.27314/d.cnki.gsscu.2022.000433.

    WANG D. Optimization of wind turbine system maintenance decision under the influence of multiple correlations [D]. Shanghai: Shanghai Ocean University, 2022. DOI: 10.27314/d.cnki.gsscu.2022.000433.
    [19] 鞠冠章, 王靖然, 崔琛, 等. 极端天气事件对新能源发电和电网运行影响研究 [J]. 智慧电力, 2022, 50(11): 77-83. DOI:  10.3969/j.issn.1673-7598.2022.11.013.

    JU G Z, WANG J R, CUI C, et al. Impact of extreme weather events on new energy power generation and power grid operation [J]. Smart power, 2022, 50(11): 77-83. DOI:  10.3969/j.issn.1673-7598.2022.11.013.
    [20] 程江洲, 冯馨以, 冯梦婷, 等. 计及气象因素的风电机组故障诊断与风险预测 [J]. 科学技术与工程, 2022, 22(22): 9645-9651. DOI:  10.3969/j.issn.1671-1815.2022.22.025.

    CHENG J Z, FENG X Y, FENG M T, et al. Wind turbine fault diagnosis and risk prediction considering meteorological factors [J]. Science technology and engineering, 2022, 22(22): 9645-9651. DOI:  10.3969/j.issn.1671-1815.2022.22.025.
    [21] 李凌飞, 孙悦, 黄莹, 等. 考虑恶劣天气影响的海上风电场及柔性直流并网系统可靠性评估 [J]. 南方电网技术, 2020, 14(12): 32-42. DOI:  10.13648/j.cnki.issn1674-0629.2020.12.005.

    LI L F, SUN Y, HUANG Y, et al. Reliability evaluation of offshore wind farm and VSC-HVDC integrated system considering the influence of extreme weather [J]. Southern power system technology, 2020, 14(12): 32-42. DOI:  10.13648/j.cnki.issn1674-0629.2020.12.005.
    [22] 李欣, 段詠程. 基于改进隐马尔可夫模型的网络安全态势评估方法 [J]. 计算机科学, 2020, 47(7): 287-291. DOI:  10.11896/jsjkx.190300045.

    LI X, DUAN Y C. Network security situation assessment method based on improved hidden Markov model [J]. Computer science, 2020, 47(7): 287-291. DOI:  10.11896/jsjkx.190300045.
    [23] 朱熀秋, 樊帅. 基于改进连续隐马尔可夫模型的六极径向主动磁轴承转子位移软测量 [J]. 中国电机工程学报, 2021, 41(11): 3933-3943. DOI:  10.13334/j.0258-8013.pcsee.201007.

    ZHU H Q, FAN S. Soft-sensing modeling for rotor displacements of six-pole radial active magnetic bearing using improved continuous hidden Markov model [J]. Proceedings of the CSEE, 2021, 41(11): 3933-3943. DOI:  10.13334/j.0258-8013.pcsee.201007.
    [24] 郑守红, 毕果, 苏史博, 等. 基于连续隐马尔可夫模型的砂轮磨削性能退化评估 [J]. 厦门大学学报(自然科学版), 2021, 60(6): 1064-1070. DOI:  10.6043/j.issn.0438-0479.202007015.

    ZHENG S H, BI G, SU S B, et al. Evaluation of grinding wheel degradation performance based on continuous hidden Markov model [J]. Journal of Xiamen University (natural science edition), 2021, 60(6): 1064-1070. DOI:  10.6043/j.issn.0438-0479.202007015.
    [25] 沈云云. 基于连续隐马尔可夫模型的语音识别抗噪问题研究 [D]. 哈尔滨: 哈尔滨工业大学, 2020. DOI: 10.27061/d.cnki.ghgdu.2020.003973.

    SHEN Y Y. Research on anti-noise of speech recognition based on continuous hidden Markov model [D]. Harbin: Harbin Institute of Technology, 2020. DOI: 10.27061/d.cnki.ghgdu.2020.003973.
    [26] 刘嘉宁, 潮铸, 钟华赞, 等. 电网调度操作过程中的天气风险源建模及分析 [J]. 广东技术师范学院学报(自然科学), 2016, 37(11): 8-12,46. DOI:  10.13408/j.cnki.gjsxb.2016.11.003.

    LIU J N, CHAO Z, ZHONG H Z, et al. Weather risk sources quantified model and analysis for power grid dispatching [J]. Journal of Guangdong Polytechnic Normal University, 2016, 37(11): 8-12,46. DOI:  10.13408/j.cnki.gjsxb.2016.11.003.
    [27] 熊小伏, 王尉军, 于洋, 等. 多气象因素组合的输电线路风险分析 [J]. 电力系统及其自动化学报, 2011, 23(6): 11-15,28. DOI:  10.3969/j.issn.1003-8930.2011.06.003.

    XIONG X F, WANG W J, YU Y, et al. Risk analysis method for transmission line combining of various meteorological factors [J]. Proceedings of the CSU-EPSA, 2011, 23(6): 11-15,28. DOI:  10.3969/j.issn.1003-8930.2011.06.003.
    [28] 王建, 姚江宁, 刘泽青, 等. 恶劣天气下配电网故障统计分析及其概率分布拟合 [J]. 电力系统保护与控制, 2022, 50(17): 143-153. DOI:  10.19783/j.cnki.pspc.211536.

    WANG J, YAO J N, LIU Z Q, et al. Fault statistical analysis and probability distribution fitting for a power distribution network in adverse weather conditions [J]. Power system protection and control, 2022, 50(17): 143-153. DOI:  10.19783/j.cnki.pspc.211536.
    [29] 高晨, 赵勇, 汪德良, 等. 海上风电机组电气设备状态检修技术研究现状与展望 [J]. 电工技术学报, 2022, 37(增刊1): 30-42. DOI:  10.19595/j.cnki.1000-6753.tces.L90275.

    GAO C, ZHAO Y, WANG D L, et al. Research status and prospect of condition based maintenance technology for offshore wind turbine electrical equipment [J]. Transactions of China electrotechnical society, 2022, 37(Suppl. 1): 30-42. DOI:  10.19595/j.cnki.1000-6753.tces.L90275.
    [30] 罗森森. 基于故障数据和监测数据的风电机组优化维修策略研究 [D]. 北京: 华北电力大学(北京), 2019. DOI: 10.27140/d.cnki.ghbbu.2019.000377.

    LUO S S. Research on optimization maintenance strategy of wind turbine based on fault data and monitoring data [D]. Beijing: North China Electric Power University (Beijing), 2019. DOI: 10.27140/d.cnki.ghbbu.2019.000377.
    [31] 邓杰, 姜飞, 王文烨, 等. 考虑电热柔性负荷与氢能精细化建模的综合能源系统低碳运行 [J]. 电网技术, 2022, 46(5): 1692-1702. DOI:  10.13335/j.1000-3673.pst.2021.1373.

    DENG J, JIANG F, WANG W Y, et al. Low-carbon optimized operation of integrated energy system considering electric-heat flexible load and hydrogen energy refined modeling [J]. Power system technology, 2022, 46(5): 1692-1702. DOI:  10.13335/j.1000-3673.pst.2021.1373.
    [32] 马志侠, 张林鍹, 郑兴, 等. 基于PEMFC-P2G与风光不确定的综合能源系统优化调度 [J]. 太阳能学报, 2022, 43(6): 441-447. DOI:  10.19912/j.0254-0096.tynxb.2022-0345.

    MA Z X, ZHANG L X, ZHENG X, et al. Optimal scheduling of integrated energy system based on PEMFC-P2G and inpact of wind power and photovoltaic uncertainty [J]. Acta energiae solaris sinica, 2022, 43(6): 441-447. DOI:  10.19912/j.0254-0096.tynxb.2022-0345.
    [33] 蓝静, 朱继忠, 李盛林, 等. 考虑碳惩罚的电化学储能消纳风光与调峰研究 [J]. 综合智慧能源, 2022, 44(1): 9-17. DOI:  10.3969/j.issn.2097-0706.2022.01.002.

    LAN J, ZHU J Z, LI S L, et al. Research on electrochemical energy storage to assist new energy consumption and peak load regulation considering carbon penalty [J]. Integrated intelligent energy, 2022, 44(1): 9-17. DOI:  10.3969/j.issn.2097-0706.2022.01.002.
    [34] 谢敏, 何润泉, 刘明波, 等. 考虑复杂振动区特性的梯级水火电低碳机组组合模型与算法研究 [J]. 电网技术, 2023, 47(2): 645-657. DOI:  10.13335/j.1000-3673.pst.2022.0900.

    XIE M, HE R Q, LIU M B, et al. Research on the model and algorithm for the hydro-thermal low-carbon unit commitment problem considering forbidden zone constraints [J]. Power system technology, 2023, 47(2): 645-657. DOI:  10.13335/j.1000-3673.pst.2022.0900.
    [35] 李建涛, 姚鸿韦, 梅德文. 碳中和目标下我国碳市场定价机制研究 [J]. 环境保护, 2021, 49(14): 24-29. DOI:  10.14026/j.cnki.0253-9705.2021.14.009.

    LI J T, YAO H W, MEI D W. Research on the pricing mechanism of China's carbon market under carbon neutralization target [J]. Environmental protection, 2021, 49(14): 24-29. DOI:  10.14026/j.cnki.0253-9705.2021.14.009.
    [36] 汪超群, 韦化, 吴思缘. 计及潮流约束的水火电力系统机组组合问题的分解–协调算法 [J]. 中国电机工程学报, 2017, 37(11): 3148-3161. DOI:  10.13334/j.0258-8013.pcsee.160677.

    WANG C Q, WEI H, WU S Y. A decomposition-coordination algorithm applied to hydro-thermal unit commitment problems with power flow constraints [J]. Proceedings of the CSEE, 2017, 37(11): 3148-3161. DOI:  10.13334/j.0258-8013.pcsee.160677.
    [37] 代江, 田年杰, 姜有泉, 等. 考虑梯级耦合的水火电检修计划与机组组合协同优化 [J]. 电力工程技术, 2022, 41(3): 83-91. DOI:  10.12158/j.2096-3203.2022.03.010.

    DAI J, TIAN N J, JIANG Y Q, et al. Collaborative maintenance scheduling and unit commitment for hydropower and thermal power systems considering cascade hydropower coupling [J]. Electric power engineering technology, 2022, 41(3): 83-91. DOI:  10.12158/j.2096-3203.2022.03.010.
    [38] 付聪, 王砚平, 刘俊磊, 等. 基于辅助优化问题的安全约束机组组合约束削减方法 [J]. 电力系统保护与控制, 2021, 49(21): 9-17. DOI:  10.19783/j.cnki.pspc.210152.

    FU C, WANG Y P, LIU J L, et al. Constraint reduction method for security-constrained unit commitment based on an auxiliary optimization problem [J]. Power system protection and control, 2021, 49(21): 9-17. DOI:  10.19783/j.cnki.pspc.210152.
    [39] 代江, 田年杰, 姜有泉, 等. 考虑天然来水随机性的水火电系统机组检修计划 [J]. 电力系统保护与控制, 2022, 50(12): 44-53. DOI:  10.19783/j.cnki.pspc.211035.

    DAI J, TIAN N J, JIANG Y Q, et al. Generator maintenance schedule of hydro-thermal power systems considering randomness of natural water inflow [J]. Power system protection and control, 2022, 50(12): 44-53. DOI:  10.19783/j.cnki.pspc.211035.
    [40] 梅竞成, 齐冬莲, 张建良, 等. 一种考虑预测电价和碳排放成本的大规模机组检修决策方法 [J]. 电子与信息学报, 2022, 44(11): 3767-3776. DOI:  10.11999/JEIT220491.

    MEI J C, QI D L, ZHANG J L, et al. A decision method of the large-scale unit maintenance scheduling considering predicted electricity price and carbon emission cost [J]. Journal of electronics & information technology, 2022, 44(11): 3767-3776. DOI:  10.11999/JEIT220491.
    [41] 俞晨玺, 孔维禄, 俞柏红, 等. 考虑需求响应的多目标机组检修调度优化 [J]. 电力系统保护与控制, 2020, 48(11): 110-118. DOI:  10.19783/j.cnki.pspc.190920.

    YU C X, KONG W L, YU B H, et al. Multi-objective optimization of generation maintenance scheduling considering demand response [J]. Power system protection and control, 2020, 48(11): 110-118. DOI:  10.19783/j.cnki.pspc.190920.
    [42] 代江, 姜有泉, 田年杰, 等. 考虑发电权交易的水火电机组检修双层优化 [J]. 高电压技术, 2022, 48(10): 4143-4153. DOI:  10.13336/j.1003-6520.hve.20210778.

    DAI J, JIANG Y Q, TIAN N J, et al. Bi-level optimization for hydrothermal unit maintenance considering generation right transaction [J]. High voltage engineering, 2022, 48(10): 4143-4153. DOI:  10.13336/j.1003-6520.hve.20210778.
    [43] WANG Y, KIRSCHEN D S, ZHONG H W, et al. Coordination of generation maintenance scheduling in electricity markets [J]. IEEE transactions on power systems, 2016, 31(6): 4565-4574. DOI:  10.1109/TPWRS.2016.2514527.
  • 加载中
图(13) / 表 (9)
计量
  • 文章访问数:  243
  • HTML全文浏览量:  147
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-12-08
  • 修回日期:  2023-02-06
  • 网络出版日期:  2023-05-26
  • 刊出日期:  2023-07-10

基于连续隐马尔可夫模型的风水火联合低碳检修优化

doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
    基金项目:  广东省基础与应用基础研究基金面上项目“虚拟电厂动态聚合演化理论与计算方法研究”(2021A1515012245);中央高校基本科研业务费重大产学研合作扶持专项“主动配电网多网多体协同优化决策技术研究”(x2dlD2201280)
    作者简介:

    何知纯,1998-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1329029641@qq.com

    谢敏,1978-,女,博士,副教授,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)minxie@scut.edu.cn

    黄莹,1999-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1419928558@qq.com

    李弋升,1998-,男,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)907353249@qq.com

    张世平,1998-,男,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1632522154@qq.com

    通讯作者: 何知纯,1998-,女,华南理工大学电气工程硕士,主要从事电力系统智能化调度与控制、电力市场等方向研究工作(e-mail)1329029641@qq.com
  • 中图分类号: TK89;TM315

摘要:   目的  新型电力系统背景下,风机的低碳检修与常规机组的协同检修等问题亟待解决。文章兼顾多属性气象因子影响和低碳性、经济性需求,建立基于连续隐马尔可夫的风水火联合低碳检修优化模型。  方法  首先,以降雨量、风速、雷电危险度为观测序列,以检修容量为隐状态序列,利用连续隐马尔可夫(Continuous Hidden Markov Model, CHMM)过程实现对风电场检修容量的动态跟踪。然后,以最优检修容量为决策依据、以总成本最小为优化目标,统筹考虑检修约束、系统控制约束等,构建风水火联合低碳检修优化模型。最后,以IEEE30节点系统进行算例展开研究。  结果  结果表明,所提模型具有更显著的经济效益和低碳特性。  结论  文章所做研究对风机的运行维护具有较高的理论价值,工程适用性较强。

English Abstract

何知纯, 谢敏, 黄莹, 李弋升, 张世平. 基于连续隐马尔可夫模型的风水火联合低碳检修优化[J]. 南方能源建设, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
引用本文: 何知纯, 谢敏, 黄莹, 李弋升, 张世平. 基于连续隐马尔可夫模型的风水火联合低碳检修优化[J]. 南方能源建设, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
HE Zhichun, XIE Min, HUANG Ying, LI Yisheng, ZHANG Shiping. Wind Power, Hydropower and Thermal Power Combined Low-Carbon Maintenance Optimization Based on Continuous Hidden Markov Model[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
Citation: HE Zhichun, XIE Min, HUANG Ying, LI Yisheng, ZHANG Shiping. Wind Power, Hydropower and Thermal Power Combined Low-Carbon Maintenance Optimization Based on Continuous Hidden Markov Model[J]. SOUTHERN ENERGY CONSTRUCTION, 2023, 10(4): 43-56. doi: 10.16516/j.gedi.issn2095-8676.2023.04.005
    • 双碳“战略”下,构建以新能源为主体的新型电力系统已大势所趋,风电因技术成熟、潜力巨大呈现迅速发展的态势,在新能源领域占据主导地位[1-2]。然而,复杂的设备结构和恶劣的运行环境不仅加剧风机故障,给风机的检修维护带来诸多不便,更引发系统失稳、经济受损等现象。从长远考虑,风电的检修维护问题亟待解决[3-5]

      目前风机检修体系为“故障检修为主、状态检修和计划检修为辅”[6-7]。故障检修又称事后检修,是发生在设备故障后的被动检修。受限于检修资源的可及性、快速恢复运行的必要性,检修事项往往流于表面,仅适用于重要性低的机组[8]。状态检修以多源异质数据监测、故障诊断为基础,评估机组运行状态和剩余寿命。在实际应用中,受恶劣天气、监测设备性能、多部件层级式耦合的影响,机组状态特征提取困难,监测数据呈现波动剧烈、灵敏不足等特点,易引起故障定位错误、检修过剩等问题[9-10]。计划检修站在全生命周期视角,统筹考虑机组劣化过程和历史故障规律,提前规划检修事项,被认为是最经济、最可行的检修方式[11-12]。目前,计划检修的研究聚焦在检修策略的制定,苏青等[13]引入有效年龄的概念,充分考虑检修次数、检修程度对机组故障演化的影响,求解风机部件的顺序维修。逯红霞等[14]基于故障有向图理论,推演风机维护阈值函数,以最小化维护成本为目标优化检修周期。符杨等[15]同样以维护成本为优化目标,综合考虑天气、设备劣化对检修策略的影响。符杨等[16]借助Copula函数进行多部件联合风险评估,以此为纽带构建双层优化模型,探寻兼顾经济性和安全性的检修策略。IRAWAN C A等[17]考虑检修资源的可及性,构建风电场检修优化模型。王东[18]将风电机组建模为串联系统,利用马尔可夫链和威尔分布描述机组寿命进程,优化检修策略。

      上述研究主要存在2个方面的问题:(1)现有研究局限于对风机内部故障演化规律的探讨,外界影响因子等对检修策略的综合影响还有待研究。作为经济调度的前提,风机的检修应在挖掘自身运行特性、气象因素的影响上,系统地实现与其他常规机组的检修协同;(2)检修维护过程忽视了低碳运行这一现实需求,对模型的优化目标考虑不充分。风机的检修引发常规机组出力增加、启停频繁等一系列现象,造成不必要的资源浪费,合理的检修计划应促使调度过程中各类资源实现多能互补。

      鉴于此,本文计及多属性气象因子(降雨量、风速、雷电危险度)对检修容量的影响,基于连续隐马尔可夫(Continuous Hidden Markov Model,CHMM)理论构建风电场检修容量跟踪模型,预估契合运行环境的检修容量。并以此为决策依据,构建风水火联合低碳检修优化模型,在兼顾系统低碳性、经济性目标的同时,实现风电与水火电的协同检修优化。IEEE30节点系统算例表明,本文所提模型的降本效益与低碳特性显著,能为风水火的联合检修维护提供有效指导。

    • 相较于水火电机组,风机的运行环境较为恶劣。大量调研表明[19-21],气象因素成为机组故障与检修的重要诱因,其中降雨量、风速、雷电危险度占据主导地位。本节以风电场为研究对象,以降雨量、风速、雷电危险度为多属性气象因子,引入CHMM理论,跟踪求解多属性气象因子波动下的风电场检修容量。

    • CHMM经马尔可夫发展而来,它以两个关联的马尔可夫随机过程表征隐状态序列(不可观测)与观测序列的联系[22-23]。其中,马尔科夫链描述隐状态的变化,随机过程则描述两序列映射关系的概率分布,观测序列的概率密度函数通常以高斯混合模型表示[24]

      风电场检修容量无法直接观测得到,但可以通过分析多属性气象因子(降雨量、风速、雷电危险度)的时变过程来间接获得,且多属性气象因子与检修容量之间存在复杂的统计关系。据此,本文将风电场检修容量跟踪问题映射为CHMM。隐状态序列(以天为时间尺度,下同)为检修容量的连续时变过程,观测序列为时变过程内获取的多属性气象因子数值。

      风电场检修容量的CHMM由一个五元组$\lambda = (N,M,{\boldsymbol{\pi}} ,A,B)$描述[25],具体如下。

      1)隐状态数目N。记隐状态集合$S = \left\{ {{s_1},{s_2}, \cdots {s_N}} \right\}$,观测序列$Q = \left\{ {{q_1},{q_2}, \cdots ,{q_T}} \right\}$,则t时刻所处状态${q_t} \in S$。本文认为风电场检修事项中,每个状态为特定检修容量区间。根据历史检修数据、恶劣天气设定标准等进行风电场检修容量的区间划分[26-28],设定5种隐状态,即N=5,其定义如表1所示。

      表 1  各隐状态物理含义

      Table 1.  Physical meaning of each hidden state

      状态检修容量物理含义
      1$(0.00 \sim 0.50){P_{\max ,w} }$少部分机组需要检修
      2$(0.50 \sim 0.75){P_{\max ,w} }$
      3$(0.75 \sim 0.90){P_{\max ,w} }$大部分机组需要检修
      4$(0.90 \sim 0.95){P_{\max ,w} }$
      5$(0.95 \sim 1.00){P_{\max ,w} }$几乎所有的机组需要检修

      2)观测值数目M:每个隐状态对应的高斯混合数。记观测值集合$V = \left\{ {{v_1},{v_2}, \cdots {v_M}} \right\}$,则t时刻观测值${o_t} \in V$,观测序列$O = \left\{ {{o_1},{o_2}, \cdots ,{o_T}} \right\}$。本文以相互关联的降雨量序列、风速序列、雷电危险度序列作为观测序列,采用3个高斯分布的线性组合无限逼近,即M=3。

      3)初始状态概率矩阵${\boldsymbol{\pi}}$${\boldsymbol{\pi}}$表示初始t=1时的概率矩阵,${\boldsymbol{\pi}} = \left\{ {{{\boldsymbol{\pi}} _i}} \right\}$,其中${{ {π}} _i}$描述状态${q_1}$${S_i}$的概率,即:

      $$ {{\boldsymbol{\pi}}_i}= P({q_1} = {S_ i})1 \leqslant i \leqslant N $$ (1)

      4)状态转移概率矩阵${\boldsymbol{A}}$${\boldsymbol{A}}$描述从状态${S_i}$(t时刻)转移到状态${S_j}$(t+1时刻)的概率,即${\boldsymbol{A}} = {({a_{ij}})_{N \times N}}$

      $$ {a_{ij}} = P({q_{t + 1}} = {S_ j}\left| {{q_t} = {S_i}} \right.),{\text{ 1}} \leqslant i,j \leqslant N $$ (2)

      5)观测概率矩阵${\boldsymbol{B}}$。CHMM采用高斯混合模型拟合各隐状态对应的观测值的概率密度函数,即:

      $$ {b_j}({o_t}) = \displaystyle \sum\limits_{m = 1}^M {{w_{j,m}}{b_{j,m}}({o_t}) = } \displaystyle \sum\limits_{m = 1}^M {{w_{j,m}}G\left( {{o_t},{\mu _{j,m}},{\displaystyle \sum\limits_{j,m}}} \right)} $$ (3)

      式中:

      M ——高斯分量个数;

      ${w_{j,m}}$${\mu _{j,m}}$${\sum _{j,m}}$ ——状态${S_j}$下第m个高斯分布的权值、均值向量、方差矩阵;

      $G({o_t},{\mu _{j,m}},{\sum _{j,m}})$ ——混合高斯概率密度函数。状态j下的第m个高斯分布$ {b_{j,m}}({o_t}) $为:

      $$ \begin{split} &{b_{j,m}}({o_t}) = G({o_t},{\mu _{j,m}},{\displaystyle \sum\limits_{ j,m}})= \\& \dfrac{1}{{\sqrt {{{(2{\boldsymbol{\pi}})}^D}\left| {{\displaystyle \sum\limits _{j,m}}} \right|} }}\exp ( - \dfrac{1}{2}{\left( {{o_t} - {\mu _{j,m}}} \right)^T}{\displaystyle \sum\limits_{j,m}}^{ - 1}({o_t} - {\mu _{j,m}})) \end{split} $$ (4)

      权重${w_{j,m}}$受以下约束限制:

      $$ {w_{j,m}} > 0,\displaystyle \sum\limits_{m = 1}^M {{w_{j,m}} = 1,1 \leqslant j \leqslant N,1 \leqslant m \leqslant {M_j}} $$ (5)

      鉴于五元组各参数联系紧密,参数${\boldsymbol{B}}$可等效为:${\boldsymbol{B}} \to \left( {w,\mu ,\displaystyle \sum\limits } \right)$,五元组$\lambda $可等效为:$\lambda \to \left( {{\boldsymbol{\pi}} ,{\boldsymbol{A}},w,\mu ,\displaystyle \sum\limits } \right)$

    • CHMM的训练本质是参数学习问题,即通过循环迭代求解最优参数$\lambda = \left( {{\boldsymbol{\pi}} ,{\boldsymbol{A}},w,\mu ,\displaystyle \sum\limits } \right)$,使得$P(O\left| \lambda \right.)$最大。本文以历史观测序列$O$为训练数据,估计最优参数$\lambda $步骤如下:

      1)基于$O$,采用Baum-Welch算法评估参数(第k次循环时求解)。

      2)采用Viterbi算法计算参数${\hat \lambda _\kappa }$下的$\hat P(O{\left| {\hat \lambda } \right._k})$

      3)判断收敛条件$\left| {\hat P(O{{\left| {\hat \lambda } \right.}_k}) - \hat P(O{{\left| {\hat \lambda } \right.}_{k - 1}})} \right| < \delta $$\delta $为设定阈值)是否成立;若成立,将当前求解的${\hat \lambda _\kappa }$作为最优参数$\lambda $;若不成立,重复前两步,直至收敛。

    • 风电场检修容量的确定可视为解码问题,即在给定观测序列$O'$和最优参数$\lambda $(1.2节得到)的基础上,采用Viterbi算法求解最贴合观测序列时变特性的隐状态序列。其核心思想如下:

      假定在最优参数和观测序列下,${q_t}$描述$ P({q_t} = {S_i}\left| {O',\lambda } \right.) $达最大时的t时刻状态,${q_{1:T}}$描述$P({q_T} = {S_i}\left| {O',\lambda } \right.)$达最大时的状态序列,即:

      $$ {q_t} = \arg \mathop {\max }\limits_{1 \leqslant i \leqslant N} [P({q_t} = {S_i}\left| {O',\lambda } \right.)] $$ (6)
      $$ {q_{1:T}} = \arg \mathop {\max }\limits_{1 \leqslant i \leqslant N} [P({q_T} = {S_i}\left| {O',\lambda } \right.)] $$ (7)
      $$ P({q_T} = {S_i}\left| {O',\lambda } \right.) = \dfrac{{P(O',{q_T} = {S_i}\left| \lambda \right.)}}{{P(O'\left| \lambda \right.)}} $$ (8)

      求解隐状态序列$Q'$即得到了风电场的日检修容量区间,依据风电场的运行特点与历史检修维护数据[29-30],求解此区间内的日检修容量值$ {\beta '_w}(t) $。即日检修容量值的确定分两步进行:(1)计及多属性气象因子构建CHMM模型,求解式(7)、式(8)得到隐状态(即风电场待检修的容量区间);(2)依据风电场的检修特性等,从式(7)、式(8)求解的容量区间内选取具体的检修容量值。两步骤的逻辑关系如图1所示,具体公式表述如下:

      图  1  风电场日检修容量值的确定步骤

      Figure 1.  Steps to determine the daily maintenance capacity of wind farm

      $$ {\beta '_w}(t) = {\gamma _t}{P_{\max ,w}} $$ (9)

      式中:

      $ {\gamma _t} $ ——t时刻隐状态所系区间内的随机系数(抽样法得到),描述风电场检修容量占比,假定t时刻隐状态为${s_1}$,则${\gamma _{_t}} \in \left( {0 \sim 0.5} \right)$

      $ {P_{\max ,w}} $ ——风电场w的最大出力。

    • 基于CHMM的风电场检修容量跟踪模型的整体流程如图2所示。

      图  2  基于CHMM的风电场检修容量跟踪模型的流程

      Figure 2.  Flow chart of wind farm maintenance capacity tracking model based on CHMM

    • 依据风机检修维护存在极易受外界因素影响、故障频繁、即坏即检等特点,基于简化模型的原则,设定如下假设:

      1)考虑到同一地区下各风机类似的系统特性和外部气象特性,风电与水火电机组共同参与检修优化。

      2)为确保多属性气象因子的时效性和精准性,根据天气预报分类原则,以月度为检修规划周期,以天为时间尺度。

      3)通常风机检修维护周期较短,参考文献,本文设定风机检修工期为3 d。

      4)在检修工期内,某一风电场检修容量值须为固定值,本文取检修区间内的最小值,具体见式(17)。

    • 本文在构建模型时以最小化经济损失、助力电力系统低碳运行为目标,统筹考虑机组检修成本、运行成本、碳排放成本等目标。以下采用各子目标(成本)加权求和的方式构建种优化目标,其中各项权值系数可以体系对该项成本值的惩罚强度[31-33]

      构建总优化目标:

      $$ \min F = \min ({\varpi _1}{F_1} + {\varpi _2}{F_2} + {\varpi _3}{F_3}) $$ (10)

      机组检修成本$ {F_1} $表示为:

      $$ {F_1} = \sum\limits_{t = 1}^T {\sum\limits_{i = 1}^Z {{C_{jx,Z,t}}{P_{\max ,i}}{a_i}(t)} } + \sum\limits_{t = 1}^T {\sum\limits_{w = 1}^W {{C_{jx,w,t}}{\beta _w}(t)} } $$ (11)

      机组运行成本$ {F_2} $表示为:

      $$ {F_2} = \sum\limits_{t = 1}^T {\sum\limits_{i = 1}^{Th} {[{c_{0,i}}{P_i}{{(t)}^2} + {c_{1,i}}{P_i}(t) + {c_{2,i}}(1 - {a_i}(t)]} } $$ (12)

      碳排放成本$ {F_3} $表示为:

      $$ {F_3} = \lambda \sum\limits_{t = 1}^T {\sum\limits_{i = 1}^{T_{\rm{h}}} {{P_i}(t){\delta _t}} } $$ (13)

      式中:

      T ——检修时段(本文以月为检修周期,T=30);

      i ——机组编号;

      Z ——水火电机组总数,Z=Th+TmThTm分别为火电机组数目、水电机组数目);

      W ——风电场总数;

      $ {C_{jx,Z,t}} $$ {C_{jx,w,t}} $ ——单位容量水火电机组、风机的检修费用,受实际情况下人力物力资源的限制,该项费用月内随时间变化,人力物力资源缺乏时,其值较高;

      $ {P_{\max ,i}} $ ——水火电机组的最大出力;

      $ {a_i}(t) $ ——机组i的检修0~1变量,$ {a_i}(t) = 1 $表示t时刻机组i处于检修状态,$ {a_i}(t) = 0 $表示t时刻机组i处于非检修状态;

      $ {\beta _w}(t) $ ——风电场w的检修连续变量,$ {\beta _w}(t) $值反映风电场中实际需要检修的容量值;

      $ {P_i}(t) $ ——火电机组it时刻的出力;

      $ {c_{0,i}} $$ {c_{1,i}} $$ {c_{2,i}} $ ——火电机组的煤耗系数。

      根据文献[34-35]$\lambda $为度电碳成本,取值为0.019元/kWh;$ {\delta _t} $为时间间隔,即每天的小时数(本文检修模型以天为尺度),取24 h;$ {\varpi }_{1}、{\varpi }_{2}、{\varpi }_{3} $分别为机组检修成本、运行成本、碳排放成本的权值,根据实际系统的碳惩罚需求而设定。

    • 优化过程考虑机组检修约束、机组特性约束、系统运行控制约束等[36-38],为风水火联合系统制定合理的检修优化策略[39]

    • 1)水火电机组检修持续性约束

      $$ {a}_{i}(t)=\left\{ {\begin{array}{c}0\text{,}t\notin [{{t}^{\prime }}_{\text{sta},i},{{t}^{\prime }}_{\text{end},i}]\\ 1\text{,}t\in [{{t}^{\prime }}_{\text{sta},i},{{t}^{\prime }}_{\text{end},i}]\end{array}} \right. $$ (14)

      式中:

      $ {t'_{{\text{sta}},i}} $$ {t'_{{\text{end}},i}} $ ——机组i实际检修计划的开始时间和终止时间。

      2)水火电机组工期约束

      $$ \displaystyle \sum\limits_{t = 1}^T {{a_i}(t) = {T_{{\text{cx}},i}}} $$ (15)

      式中:

      $ {T_{{\rm{cx}},i}} $ ——机组i的检修工期。

    • 1)检修容量平衡方程

      风电场的实际检修容量由检修0~1变量和预测检修容量决定,即系统地兼顾模型优化过程与外界多属性气象因子的共同作用。

      $$ {\beta _w}(t) = {\beta '_w}(t) \cdot {b_w}(t) $$ (16)

      式中:

      $ {\beta '_w}(t) $ ——t时刻风电场w的预测检修容量,经第1节模型求得,即所述CHMM的隐状态对应值(式9);

      $ {b_w}(t) $ ——风电场检修0~1变量;$ {b_w}(t) = 1 $表示t时刻风电场w处于检修状态,$ {b_w}(t) = 0 $表示t时刻风电场w处于非检修状态。

      2)检修容量一致性约束

      依据2.1节所述假设4,在优化得到的检修区间内,检修容量应保持固定且为区间最小值:

      $$ \begin{gathered} {\beta _w}({{t'}_{{\text{sta}},w}}) = {\beta _w}({{t'}_{{\text{sta}},w}} + 1) = {\beta _w}({{t'}_{{\text{end}},w}}) \\ {\text{ }} = \min \left\{ {{\beta _w}({{t'}_{{\text{sta}},w}}),{\beta _w}({{t'}_{{\text{sta}},w}} + 1),{\beta _w}({{t'}_{{\text{end}},w}})} \right\} \\ \end{gathered} $$ (17)

      式中:

      $ {t'_{{\text{sta}},w}} $$ {t'_{{\text{end}},w}} $——风电场w实际检修计划的开始时间和终止时间。

      3)检修持续性约束

      $$ {b}_{w}(t)=\left\{ {\begin{array}{c}0\text{,}t\notin [{{t}^{\prime }}_{\text{sta},w},{{t}^{\prime }}_{\text{end},w}]\\ 1\text{,}t\in [{{t}^{\prime }}_{\text{sta},w},{{t}^{\prime }}_{\text{end},w}]\end{array}} \right. $$ (18)

      4)检修工期约束

      $$ \sum\limits_{t = 1}^T {{b_w}(t) = {T_{{\text{cx}},w}}} $$ (19)

      式中:

      $ {T_{{\text{cx}},w}} $——风电场w的检修工期,根据假设3,$ {T_{{\text{cx}},w}} = 3 $

    • 1)火电机组功率上下限约束

      $$ [1 - {a_i}(t)]{P_{i,\min }} \leqslant {P_i}(t) \leqslant [1 - {a_i}(t)]{P_{i,\max }} $$ (20)

      式中:

      $ {P_{i,\min }} $$ {P_{i,\max }} $($i \in {\varOmega _{T_{\rm{h}}}}$)——火电机组i的最小技术出力、最大技术出力。

      2)水电机组功率上下限约束

      $$ 0 \leqslant {P_i}(t) \leqslant [1 - {a_i}(t)]{P_{i,\max }} $$ (21)

      式中:

      $ {P_{i,\max }} $($i \in {\varOmega _{T_{\rm{m}}}}$)——水电机组i的最大技术出力。

      3)风电场功率约束

      $$ {P_w}(t) = {P'_w}(t) \cdot [1 - {\gamma _t}{b_w}(t)] $$ (22)

      式中:

      $ {P_w}(t) $$ {P'_w}(t) $ ——风电场w的调度出力、预测出力;

      $ {\gamma _t} $ ——风电场检修容量占比,由第1节模型求得(式9)。

    • 1)功率平衡约束

      $$ \sum\limits_{i = 1}^Z {{P_i}(t)} + \sum\limits_{w = 1}^W {{P_w}(t)} \pm \sum\limits_{h = 1}^H {{P_h}(t)} = {D_t} $$ (23)

      式中:

      $ {P_h}(t) $——t时刻外区域与本区域电网的交换功率;

      H ——联络线总数;

      Dt ——本区域的日最大负荷,为系统总负荷值。

      2)系统备用约束

      国内外学者在机组检修优化问题中很少详细考虑失负荷量,通常将失负荷量的考虑计入系统备用约束中[40-42],即:

      $$ \sum\limits_{i = 1}^Z {{P_{i,\max }} - \sum\limits_{i = 1}^Z {{P_{i,\max }}} {a_i}(t) - (1 + \theta ){D_t} \geqslant 0} $$ (24)

      式中:

      $ \theta $——依据系统运行特点设定的最佳备用系数,本文取$ \theta = 0.05 $

      3)系统潮流约束

      $$ - {P_{l,\max }} \leqslant \sum\limits_{l \in L}^{} {[{P_i}(t)} + {P_w}(t) + {P_h}(t) - {D_t}] \leqslant {P_{l,\max }} $$ (25)

      式中:

      $ P_{l, {\rm{max}}} $ ——输电线路$ l $的最大传输功率;

      L ——线路总数。

    • 将第1节求出的风电场检修容量参数$[{\beta '_w}(t),{\gamma _t}]$作为优化模型的边界条件,受制于同一检修工期内风电场的检修容量一致,即假设4和式(17)必须满足的要求下,本文迭代求解模型,求解流程如图3所示。具体步骤如下:

      图  3  风水火联合低碳检修优化模型的求解步骤

      Figure 3.  Solving steps of the optimization model for wind power, hydropower and thermal power combined low-carbon maintenance

      1)将$[{\beta '_w}(t),{\gamma _t}]$作为初始边界条件,求解风水火联合低碳检修优化模型,得到风电场初始检修计划$ {b_w}(t) $

      2)按照式(17)修改检修工期$ ({t'_{{\text{sta}},w}},{t'_{{\text{end}},w}}) $内的$ {\beta '_w}(t) $$ {\gamma _t} $值,动态调整初始边界条件${[{\beta '_w}(t),{\gamma _t}]^k }$(第$ k $次调整)。

      3)以调整后的${[{\beta '_w}(t),{\gamma _t}]^k }$为边界条件循环求解优化模型,重复两阶段步骤,直至$ {b_w}{(t)^k } ={b_w}{(t)^{k - 1}} $($ k \geqslant 2 $)。

    • 基于模型特性和时效性拆分观测序列,以过去两年的数据为训练集,以未来一个月的数据为预测集。原始观测序列包含降雨量、风速、雷电危险度3个维度,前2个维度可由气象机构直接监测得到,雷电危险度受所在片区的气象数据、雷电天气强度和风机特性共同作用,设定原则参考文献[28]。原始观测序列数据集形式如表2所示。

      表 2  观测序列数据集形式

      Table 2.  Observation sequence data set form

      日期降雨量/mm风速/(m·s−1)雷电危险度
      D38151.680.3
      D391679.710.5
      ……………………
    • 经CHMM训练,采用Baum-Welch算法求解模型参数$\lambda $表3~表5所示。

      表 3  初始状态概率矩阵${\boldsymbol{\pi}} $

      Table 3.  Initial state probability matrix ${\boldsymbol{\pi}} $

      初始概率 ${s_1}$${s_2}$${s_3}$${s_4}$${s_5}$
      ${\boldsymbol{\pi}}$10000

      表 4  状态转移概率矩阵${\boldsymbol{A}}$

      Table 4.  State transition probability matrix ${\boldsymbol{A}}$

      隐状态 ${s_1}$${s_2}$${s_3}$${s_4}$${s_5}$
      ${s_1}$0.263 00.483 90.161 90.035 00.056 2
      ${s_2}$0.385 70.261 30.205 00.039 00.109 1
      ${s_3}$0.301 90.122 50.373 20.177 10.025 3
      ${s_4}$0.263 60.000 00.139 90.458 00.138 5
      ${s_5}$0.211 00.156 20.288 50.055 60.288 7
    • 将3.1.2节求出的参数作为CHMM的最优参数$\lambda $,采用Viterbi算法求解未来一个月与观测序列关联的最优隐状态序列,最优观测值-隐状态分布如图4所示,结果如表6所示。

      图  4  最优观测值-隐状态分布

      Figure 4.  Hidden state distribution of optimal observation value

      表 5  隐状态概率矩阵${\boldsymbol{B}}$

      Table 5.  Hidden state probability matrix ${\boldsymbol{B}}$

      日期s1s2s3s4s5
      D1 0.791 624 7 0.207 005 4 0.001 272 4 0.000 086 1 0.000 011 3
      D2 0.038 939 9 0.960 987 4 0.000 016 6 0.000 052 3 0.000 003 8
      D3 0.001 347 6 0.998 619 9 0.000 000 4 0.000 031 2 0.000 000 9
      D4 0.791 687 4 0.203 395 1 0.004 805 6 0.000 099 1 0.000 012 7
      D5 0.785 881 3 0.212 917 9 0.001 102 6 0.000 086 8 0.000 011 4
      …… …… …… …… …… ……

      表 6  未来一月的最优隐状态序列

      Table 6.  Optimal hidden state sequence for the next month

      隐状态天数/d天数占比/%
      11240.0
      21033.3
      3516.7
      400.0
      5310.0

      由上述图表可知,未来一个月内,检修容量处于隐状态1的天数最多,占比40%,大都分布在降雨量和雷电危险度较低的时刻,此时外界天气状况尚佳,有少部分(0%~50%)机组需要安排检修。处于隐状态4、隐状态5的天数较少,占比10%,即实际运行中安排90%以上机组检修的情况极为罕见,这与有限的检修资源(人力、预算等)、亟待满足的负荷需求密切相关。同时,隐状态2、隐状态3出现的频数高于隐状态5,且相应的天气状况整体更为恶劣,这是因为极其严峻的外部气象环境在加剧机组的损害,同时也给检修维护工作带来了极大的不便,此时,安排检修的机组比例相对更低。经CHMM求出的风电场检修容量充分考虑了外界气象因素对检修容量的影响,为风电场参与系统整体检修优化提供有利参照。

    • 本文通过IEEE30节点系统算例验证模型($ {\varpi }_{1}、 {\varpi }_{2}、{\varpi }_{3} $均取为1),系统包含30个节点、6台常规水火电机组、41条线路,在5、8、11、27节点接有4座风场。系统拓扑图如图5所示(W为风电,Th为火电,Tm为水电),4座风电场的预测功率曲线如图6所示,日负荷预测曲线如图7所示,以BPA运行数据中各节点负荷数据为基础,等比分配总负荷得到各节点负荷,各机组特性参数和检修工期如表7所示,水火电机组/风电场单位容量的检修费用曲线如图8所示。所采用的计算机为Inter(R) Core(TM) i7-9700 CPU @ 3.00GHz(8 CPUs)。

      图  5  IEEE30节点系统拓扑图

      Figure 5.  IEEE 30-node system topology diagram

      图  6  风电场预测功率曲线

      Figure 6.  Forecast power of wind farm

      图  7  风电场日负荷预测曲线

      Figure 7.  Daily load forecasting curve of wind farm

      表 7  各机组特性参数和检修工期

      Table 7.  Characteristic parameters and maintenance duration of each WTGS

      机组$ {P_{i,\max }} $/
      MW
      $ {P_{i,\min }} $/
      MW
      c2/
      (元·MW−1)
      c1/
      (元·MW−1)
      c0/
      (元·MW−1)
      检修工期/
      d
      i1 80.0 0 0.140 0 14.00 0 10
      i2 80.0 0 0.121 7 12.17 0 7
      i3 50.0 0 4
      i4 55.0 0 0.058 0 22.60 0 5
      i5 30.0 0 0.173 9 20.86 0 3
      i6 40.0 0 0.173 9 20.86 0 4
      w1 15.0 0 3
      w2 22.5 0 3
      w3 30.0 0 3
      w4 30.0 0 3

      图  8  水火电/风电场单位容量的检修费用曲线

      Figure 8.  Maintenance cost curve of unit capacity of hydropower station, thermal power station/wind farm

    • 经GAMS仿真计算,规划周期内优化的各风电场检修安排如表8所示,各机组检修优化决策结果如图9所示。

      表 8  规划周期内各风电场检修安排

      Table 8.  Maintenance arrangement of each wind farm within the planning period

      风电场开始时间/d持续时间/d结束时间/d检修容量/MW
      w1283306.0
      w22232513.5
      w32832016.5
      w473922.5

      图  9  规划期内各机组检修计划

      Figure 9.  Maintenance plan of each WTGS in the planning period

      由上述图表可知,风电场大都安排在月末检修,此时段处于负荷低谷期、单位容量检修费用较低时段,此检修安排兼顾系统运行稳定性与经济性。大部分水火电机组的检修安排在中下旬(D15~D25),此时段负荷需求呈现大幅度下降的趋势,如此可保证负荷高峰期内大部分机组处于持续运行状态,满足系统备用需求。整个规划周期内,实现了风水火的协同检修、错峰检修,保障了机组出力平稳,有力地促进系统“等备用运行”。

      检修优化过程中,系统的水火电出力、风电出力随时间的变化曲线如图10所示。

      图  10  各类型机组出力变化

      Figure 10.  Output change of each type of WTGS

      图10可知,水火电出力曲线与负荷曲线的变化趋势较为相近,负荷高峰期即水火电机组出力充裕时期。D5D15D24前后,受限于部分水火电机组的检修,呈现水火电出力欠缺、风电出力充足的运行特性。月末,风电场内部分机组需要检修,风电出力不足,此时段负荷需求几乎由水火电供给。整体上,各机组出力变化特性反映了优化模型对机组检修和系统运行的有效协同。

      具体地,各风电场的出力曲线如图11所示,差值部分表征检修引起的出力缺额,各缺额时段对应检修区间。w1w3检修工期内风电出力较低,出力缺额较少,检修对整个风电场出力的影响较小;w4检修工期内出力较高、检修容量最大,此时出力缺额最多,检修对出力的影响最为显著。综合分析,检修计划的优化充分考虑了风电场出力特性,符合将检修安排在出力低谷时段这一现实要求。

      图  11  各风电场的出力曲线

      Figure 11.  Output curve of each wind farm

    • 整个优化过程中,各成本变化曲线如图12所示。运行成本与火电出力紧密相关,两曲线呈现相同的变化趋势。D19~D22内检修成本曲线处于峰部,分析可知,此时容量较大的机组i1i4处于检修状态。碳排放成本曲线较为平滑缓和,表明规划期内火电机组出力保持平稳。

      图  12  优化过程中各成本变化曲线

      Figure 12.  Cost change curves in the optimization process

      文献[43]模型在检修优化时忽略了对低碳的要求,为方便阐述,称文献[43]模型为“对照模型”。为验证本文模型的优越性,应用本文算例求解两模型,两模型求解检修区间对比图如图13所示,得到的各项子结果如表9所示。

      图  13  两模型求解检修区间对比图

      Figure 13.  Comparison diagram of maintenance interval solved by two models

      表 9  不同模型求解结果对比

      Table 9.  Comparison of solution results of different models

      参数优化模型对照模型对比
      总成本/元 2 628 384 2 629 746 −0.05%
      碳排放成本/元 268 731.3 269 106 −0.14%
      运行成本/元 1 881 674 1 892 588 −0.58%
      检修成本/元 477 979 468 052.1 2.12%
      水火电检修成本/元 450 653.2 447 661.8 0.67%
      风电检修成本/元 27 325.8 20 390.33 34.01%
      水火电出力/MW 5 893.229 5 901.448 −0.14%
      风电出力/MW 922.654 2 914.435 9 0.90%

      表9可知,相比于对照模型,优化模型在目标函数中考虑了碳排放成本,求解得到的水火电出力、和碳排放成本均降低0.14%,运行成本降低15.8%,而风电出力增加0.9%,表明优化模型以促进风电出力来平抑水火电出力变化的影响。优化模型求解的检修区间内风机出力较小,有利于风电消纳、降低水火电出力,但此时段单位检修容量成本和检修容量值均更高(见图8),故总检修成本有所增加,由于优化模型优化的是3个子目标的加权求和,系统总成本降低0.05%,增加的检修成本得到有效补偿,表明优化模型以牺牲检修成本满足对经济效益的追求。对比各项指标,优化模型在改善风电出力、助力低碳、降低成本方面整体上优于对照模型,即优化模型在合理安排检修计划的同时,实现了系统的低碳、经济运行,为“双碳”战略下风水火的联合低碳检修优化提供有力支撑。

    • 本文计及多属性气象因子,建立基于CHMM的风电场检修容量跟踪模型。以跟踪模型求解结果为边界,以低碳、经济运行为优化目标,统筹考虑机组运行特性、检修控制约束等,建立风水火联合低碳检修优化模型,通过IEEE30节点系统进行验证,得到以下结论:

      1)本文充分计及多属性气象因子对风机检修的影响,制定合理的检修容量配置方案,为优化模型提供契合现实的数据。

      2)本文兼顾多属性气象因子和系统运行优化的共同作用,风电场的检修大都被安排在出力低谷期,有力地促进了风电消纳。

      3)本文模型综合考虑系统运行特性和低碳目标,实现了风水火协同检修、错峰检修,有效平抑机组出力波动性,助力系统平稳运行。

      4)与对照模型相比,本文模型的优化程度更佳,更好地实现了“双碳”目标下的低碳检修和经济运行。

      但本文也存在一些不足,受制于系统内风机性能、算例规模,模型的优化空间有限,下一步将完善模型,计及多类型电源,实现多场景下的机组联合低碳检修优化。

参考文献 (43)

目录

    /

    返回文章
    返回