CN109993445B - 一种考虑光伏预测误差的综合能源系统脆弱性评估方法 - Google Patents

一种考虑光伏预测误差的综合能源系统脆弱性评估方法 Download PDF

Info

Publication number
CN109993445B
CN109993445B CN201910271051.0A CN201910271051A CN109993445B CN 109993445 B CN109993445 B CN 109993445B CN 201910271051 A CN201910271051 A CN 201910271051A CN 109993445 B CN109993445 B CN 109993445B
Authority
CN
China
Prior art keywords
prediction error
energy system
load
power generation
photovoltaic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910271051.0A
Other languages
English (en)
Other versions
CN109993445A (zh
Inventor
丁凯
钱一民
胡畔
李伟
王易
陈乔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910271051.0A priority Critical patent/CN109993445B/zh
Publication of CN109993445A publication Critical patent/CN109993445A/zh
Application granted granted Critical
Publication of CN109993445B publication Critical patent/CN109993445B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003Load forecast, e.g. methods or systems for forecasting future load demand
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80Management or planning
    • Y02P90/82Energy audits or management systems therefor
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Power Engineering (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种考虑光伏预测误差的综合能源系统脆弱性评估方法,属于电力系统与新能源技术领域。所述方法包括以下步骤:获取综合能源系统的设备和负荷数据;构建设备运行状态及负荷需求和光伏发电预测误差的累积分布函数;利用蒙特卡洛模拟对设备是否故障和预测误差大小进行采样;计算综合能源系统的脆弱性指标;判断脆弱性指标收敛精度。用于评估综合能源系统运行时一旦发生风险所造成结果的严重程度,帮助综合能源系统的调度与维护人员合理分析系统运行状态并及时发现脆弱环节进行重点关注,实现综合能源系统的精准运维。

Description

一种考虑光伏预测误差的综合能源系统脆弱性评估方法
技术领域
本发明涉及电力系统与新能源技术领域,具体是一种考虑光伏预测误差的综合能源系统脆弱性评估方法。
背景技术
随着经济、社会的不断发展,对能源的需求迅速增长。煤、石油等化石燃料储量有限,能源问题日渐紧迫。这些问题迫使研究人员寻找一种更加高效、可持续的能源利用方式。同时,伴随着分布式发电技术、微电网技术的发展,以电能作为能量传输和利用枢纽的综合能源系统架构被研究人员提出,并作为能量梯级利用、综合利用、清洁利用的方式展开相关研究和应用。发电厂与用户构成的热电联产系统作为最基本综合能源系统,已广泛应用于北方冬季的供暖和供电需求。燃气轮机、光伏发电、压缩机、锅炉、集热热泵等设备的加入与一体化设计,使得气-电、光-电、电-冷、电-热、热-电之间的转化更为灵活;系统内能量的梯级利用使得能源利用效率得到提高。综合能源系统作为能源和环境问题的解决手段之一,已得到世界范围内的普遍认可。
但相比于传统的电力系统,综合能源系统存在冷、热、电、气多种能源形式,设备多、技术杂、能量转化过程复杂多样,管理起来更为困难,一种能源供给出现问题往往会牵扯到其他能源的供给。而对于存在分布式光伏发电的综合能源系统,光伏发电的间歇性与不确定性更成为了综合能源系统供能环节的一大风险因素。因此,对综合能源系统运行风险的预判尤为重要,合理开展综合能源系统的脆弱性分析对降低系统运行风险、提高供能可靠性具有实际意义。
发明内容
本发明的目的是提供一种综合能源系统脆弱性评估的方法,用于评估综合能源系统运行时一旦发生风险所造成结果的严重程度,帮助综合能源系统的调度与维护人员合理分析系统运行状态并及时发现脆弱环节进行重点关注,实现综合能源系统的精准运维。
一种考虑光伏预测误差的综合能源系统脆弱性评估方法,具体包含以下步骤:
步骤1、获取综合能源系统的设备数据和负荷数据,根据所获取的设备数据和负荷数据构建综合能源系统功率交换模型,所述综合能源系统功率交换模型用于计算电负荷需求Le、热负荷需求Lh和冷负荷需求量Lc
步骤2、构建设备运行状态及负荷需求和光伏发电预测误差的累积分布函数,其中设备运行状态的概率服从均匀分布,负荷需求和光伏发电预测误差发生的概率服从高斯分布;
步骤3、利用蒙特卡洛模拟,基于步骤2所得累积分布函数对设备运行状态和预测误差大小进行随机抽样,获得各设备采样后新的设备运行状态、以及考虑预测误差后新的负荷需求量电负荷需求
Figure BDA0002018394020000021
热负荷需求
Figure BDA0002018394020000022
冷负荷需求
Figure BDA0002018394020000023
和光伏发电量
Figure BDA0002018394020000024
步骤4、计算步骤3中蒙特卡洛模拟后考虑光伏预测误差和设备运行状态的实际供给量
Figure BDA0002018394020000025
根据考虑预测误差后新的负荷需求量
Figure BDA0002018394020000026
以及考虑光伏预测误差和设备运行状态的实际供给量
Figure BDA0002018394020000027
计算电负荷损失值ΔLe、热负荷损失值ΔLh、冷负荷损失值ΔL;
步骤5、重复步骤2-4多次得到三组数据,分别计算三组数据的均值和方差得到能源系统脆弱性指标的收敛精度,在收敛精度满足预设条件时输出能源系统脆弱性评估的最终结果,即电负荷、热负荷和冷负荷损失值。
进一步的,所述步骤1中,所述设备数据包括能量分配系数v、设备运行状态S、设备效率η和设备数量N,所述负荷数据包括能量供给量P、能量需求量L。
进一步的,所述步骤1构建的综合能源系统功率交换模型可用下式进行描述:
Figure BDA0002018394020000031
式中,Le、Lh、Lc分别表示电负荷需求、热负荷需求和冷负荷需求量;vDN表示从配电网传输至综合能源系统中电能用以供给电负荷需求的分配系数,
Figure BDA0002018394020000032
表示第i台光伏发电设备产生的电能用以供给电负荷需求的分配系数,
Figure BDA0002018394020000033
表示第j台燃气轮机产生的电能用以供给电负荷需求的分配系数;SDN表示配电网运行状态,
Figure BDA0002018394020000034
表示第i台光伏发电设备的运行状态,
Figure BDA0002018394020000035
表示第j台燃气轮机的运行状态,
Figure BDA0002018394020000036
表示第k台中央空调的运行状态,数值0表示故障,1表示正常运行;
Figure BDA0002018394020000037
Figure BDA0002018394020000038
分别表示第j台燃气轮机的气-电转化效率和气-热转化效率;
Figure BDA0002018394020000039
表示第k台中央空调的电制冷效率;
Figure BDA00020183940200000310
表示配电网传输至综合能源系统的电能总量,
Figure BDA0002018394020000041
表示第i台光伏发电设备所发电能总量,Pj,g表示天然气网传输至第j台燃气轮机所包含的能量大小。
进一步的,所述步骤2中,对于任意设备的运行状态,其均匀分布的累积分布函数均满足:
Figure BDA0002018394020000042
对于任意负荷和光伏发电量预测误差的累积分布函数均满足:
Figure BDA0002018394020000043
式中,预测误差的均值μ、方差σ通过统计预测误差,利用距估计或点估计法得到。
进一步的,所述步骤3中,采样后通过下式计算各设备采样后新的运行状态:
Figure BDA0002018394020000044
式中,
Figure BDA0002018394020000045
表示第i台光伏发电设备新的运行状态,
Figure BDA0002018394020000046
表示第j台燃气轮机新的运行状态,
Figure BDA0002018394020000047
表示第k台中央空调新的运行状态;
Figure BDA0002018394020000048
Figure BDA0002018394020000049
分别表示光伏发电设备、燃气轮机、中央空调发生故障的概率。
进一步的,所述步骤3中,考虑预测误差后,新的负荷需求量
Figure BDA0002018394020000051
和光伏发电量
Figure BDA0002018394020000052
变为:
Figure BDA0002018394020000053
式中,
Figure BDA0002018394020000054
为第i台光伏发电设备发电量的预测值,
Figure BDA0002018394020000055
Figure BDA0002018394020000056
分别为第i台光伏发电设备发电量预测误差的均值和方差;Le,f为电负荷需求的预测值,μL和σL分别为电负荷需求预测误差的均值和方差;其他L、μ、σ下标h和c同理表示为热负荷、冷负荷需求的相关变量。
进一步的,所述步骤4具体为:
将步骤3中求得的考虑预测误差后新的光伏发电量
Figure BDA0002018394020000057
光伏发电设备、燃气轮机、中央空调采样后新的运行状态
Figure BDA0002018394020000058
Figure BDA0002018394020000059
代入下式计算得出考虑光伏预测误差和设备运行状态的实际供给量
Figure BDA00020183940200000510
Figure BDA00020183940200000511
Figure BDA00020183940200000512
其次,根据考虑预测误差后新的负荷需求量
Figure BDA00020183940200000513
以及考虑光伏预测误差和设备运行状态的实际供给量
Figure BDA00020183940200000514
计算电负荷损失值ΔLe、热负荷损失值ΔLh、冷负荷损失值ΔLc
Figure BDA0002018394020000061
进一步的,所述步骤5中具体为:
重复步骤2至步骤4共m次,m为大于等于3的自然数,得到三组数据:
Figure BDA0002018394020000062
分别计算三组数据的均值μe、μh、μc和方差σe、σh、σc,代入下式计算脆弱性指标的收敛精度θ:
Figure BDA0002018394020000063
进一步的,所述预设条件为:脆弱性指标的收敛精度θ<0.03。
本发明具有的有益的效果是:能够以具体数值的形式展示综合能源系统的脆弱性大小,并能分别展示电、冷、热子系统的脆弱性大小,同时考虑了负荷预测误差和光伏发电预测误差对评估结果的影响,相比于没有考虑预测误差的方法,本方法的评估结果更能反映系统运行的实际情况,帮助运维人员更好的安排设备维护。
附图说明
图1本发明考虑光伏预测误差的综合能源系统脆弱性评估方法的流程图;
图2是综合能源系统架构图。
具体实施方式
下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述。
下面将结合附图,对本发明实施例中的技术方案进行清楚、完整地描述。
综合能源系统如图1所示,结合附图2进行说明,一种考虑光伏预测误差的综合能源系统脆弱性评估方法包括以下步骤:
步骤1.获取综合能源系统的设备和负荷数据,根据所获取的设备和负荷数据构建综合能源系统功率交换模型,所述综合能源系统功率交换模型用于计算电负荷需求Le、热负荷需求Lh和冷负荷需求量Lc,其可用于步骤4中负荷损失指标的计算。需要获取的具体数据包括设备数据:能量分配系数v,设备运行状态S,设备效率η和设备数量N,以及负荷数据:能量供给量P,能量需求量L,这些数据可以通过读取计量表、设备铭牌和综合能源系统的控制中心得到。所构建的综合能源系统功率交换模型可用下式进行描述:
Figure BDA0002018394020000071
式中,Le、Lh、Lc分别表示电负荷需求、热负荷需求和冷负荷需求量;vDN表示从配电网传输至综合能源系统中电能用以供给电负荷需求的分配系数,
Figure BDA0002018394020000072
表示第i台光伏发电设备产生的电能用以供给电负荷需求的分配系数,
Figure BDA0002018394020000073
表示第j台燃气轮机产生的电能用以供给电负荷需求的分配系数;SDN表示配电网运行状态,
Figure BDA0002018394020000074
表示第i台光伏发电设备的运行状态,
Figure BDA0002018394020000081
表示第j台燃气轮机的运行状态,
Figure BDA0002018394020000082
表示第k台中央空调的运行状态,数值0表示故障,1表示正常运行;
Figure BDA0002018394020000089
Figure BDA00020183940200000810
分别表示第j台燃气轮机的气-电转化效率和气-热转化效率;
Figure BDA0002018394020000084
表示第k台中央空调的电制冷效率;
Figure BDA0002018394020000085
表示配电网传输至综合能源系统的电能总量,
Figure BDA0002018394020000086
表示第i台光伏发电设备所发电能总量,Pg表示天然气网传输至综合能源系统所包含的能量大小。
步骤2.构建设备运行状态及负荷需求和光伏发电预测误差的累积分布函数,其中光伏发电设备、燃气轮机、中央空调发生故障的概率服从均匀分布;负荷需求和光伏发电预测误差发生的概率服从高斯分布。
对于任意设备的运行状态,其均匀分布的累积分布函数均满足:
Figure BDA0002018394020000087
对于任意负荷和光伏发电量预测误差的累积分布函数均满足:
Figure BDA0002018394020000088
式中,预测误差的均值μ、方差σ通过统计预测误差,利用距估计或点估计法易于得到。
步骤3.利用蒙特卡洛模拟,基于步骤2所得累积分布函数对设备运行状态(例如是否发生故障)和预测误差大小进行随机抽样,获得设备运行状态、以及负荷需求和光伏发电量的预测误差。
采样后,可通过下式计算各设备采样后新的运行状态,区别于步骤1中的运行状态变量,新的运行状态变量加入了上标“’”作为区别。由于我国城市配电网供电可靠率高达99.9%,配电网完全故障至无法提供电能的概率可以忽略不计,因此令步骤1中的SDN恒等于1。
Figure BDA0002018394020000091
式中,
Figure BDA0002018394020000092
表示第i台光伏发电设备新的运行状态,
Figure BDA0002018394020000093
表示第j台燃气轮机新的运行状态,
Figure BDA0002018394020000094
表示第k台中央空调新的运行状态;
Figure BDA0002018394020000095
Figure BDA0002018394020000096
分别表示光伏发电设备、燃气轮机、中央空调发生故障的概率。
考虑预测误差后,新的负荷需求量
Figure BDA0002018394020000097
和光伏发电量
Figure BDA0002018394020000098
变为:
Figure BDA0002018394020000099
式中,
Figure BDA00020183940200000910
为第i台光伏发电设备发电量的预测值,
Figure BDA00020183940200000911
Figure BDA00020183940200000912
分别为第i台光伏发电设备发电量预测误差的均值和方差;Le,f为电负荷需求的预测值,μL和σL分别为电负荷需求预测误差的均值和方差;其他L、μ、σ下标h和c同理表示为热负荷、冷负荷需求的相关变量。
步骤4.计算综合能源系统的脆弱性指标,具体指标包括:电负荷损失指标、热负荷损失指标、冷负荷损失指标。
首先需要计算步骤3中蒙特卡洛模拟后考虑光伏预测误差和设备运行状态(有无故障)的实际供给量,具体的,将步骤3中求得的考虑预测误差后新的光伏发电量
Figure BDA0002018394020000101
光伏发电设备、燃气轮机、中央空调采样后新的设备运行状态
Figure BDA0002018394020000102
Figure BDA0002018394020000103
代入下式计算得出考虑光伏预测误差和设备运行状态(有无故障)的实际供给量
Figure BDA0002018394020000104
Figure BDA0002018394020000105
Figure BDA0002018394020000106
其次,根据考虑预测误差后新的负荷需求量
Figure BDA0002018394020000107
以及考虑光伏预测误差和设备运行状态(有无故障)的实际供给量
Figure BDA0002018394020000108
Figure BDA0002018394020000109
计算电负荷损失值ΔLe、热负荷损失值ΔLh、冷负荷损失值ΔLc
Figure BDA00020183940200001010
步骤5.重复步骤2至步骤4共m次,得到三组数据:
Figure BDA00020183940200001011
分别计算三组数据的均值μe、μh、μc和方差σe、σh、σc,代入下式计算收敛精度:
Figure BDA0002018394020000111
若脆弱性指标的收敛精度满足预设条件(例如θ<0.03),则输出最后一次指标结果,即脆弱性评估的最终结果为电负荷、热负荷和冷负荷损失值分别为μe、μh和μc;反之,则重复步骤5直到获得到一组收敛精度满足预设条件(例如θ<0.03)的值,由此可得到三组脆弱性评估值。
经过以上步骤即完成了考虑光伏预测误差的综合能源系统脆弱性评估方法,损失值越大,则系统越脆弱,对于脆弱性高的电、冷、热部分需要设备运行维护人员格外关注。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何属于本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (8)

1.一种考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于包括以下步骤:
步骤1、获取综合能源系统的设备数据和负荷数据,根据所获取的设备数据和负荷数据构建综合能源系统功率交换模型,所述综合能源系统功率交换模型用于计算电负荷需求Le、热负荷需求Lh和冷负荷需求量Lc
步骤2、构建设备运行状态及负荷需求和光伏发电预测误差的累积分布函数,其中设备运行状态的概率服从均匀分布,负荷需求和光伏发电预测误差发生的概率服从高斯分布;
步骤3、利用蒙特卡洛模拟,基于步骤2所得累积分布函数对设备运行状态和预测误差大小进行随机抽样,获得各设备采样后新的设备运行状态、以及考虑预测误差后新的负荷需求量电负荷需求
Figure FDA0003129988760000011
热负荷需求
Figure FDA0003129988760000012
冷负荷需求
Figure FDA0003129988760000013
和光伏发电量
Figure FDA0003129988760000014
步骤4、计算步骤3中蒙特卡洛模拟后考虑光伏预测误差和设备运行状态的实际供给量
Figure FDA0003129988760000015
根据考虑预测误差后新的负荷需求量
Figure FDA0003129988760000016
以及考虑光伏预测误差和设备运行状态的实际供给量
Figure FDA0003129988760000017
计算电负荷损失值ΔLe、热负荷损失值ΔLh、冷负荷损失值ΔLc
步骤5、重复步骤2-4多次得到三组数据,分别计算三组数据的均值和方差得到能源系统脆弱性指标的收敛精度,在收敛精度满足预设条件时输出能源系统脆弱性评估的最终结果,即电负荷、热负荷和冷负荷损失值;
所述步骤4具体为:
将步骤3中求得的考虑预测误差后新的光伏发电量
Figure FDA0003129988760000018
光伏发电设备、燃气轮机、中央空调采样后新的运行状态
Figure FDA0003129988760000019
Figure FDA00031299887600000110
代入下式计算得出考虑光伏预测误差和设备运行状态的实际供给量
Figure FDA00031299887600000111
Figure FDA00031299887600000112
Figure FDA0003129988760000021
其中,
Figure FDA0003129988760000022
分别表示考虑光伏预测误差和设备运行状态的电负荷、热负荷和冷负荷的实际供给量;vDN表示从配电网传输至综合能源系统中电能用以供给电负荷需求的分配系数,
Figure FDA0003129988760000023
表示第i台光伏发电设备产生的电能用以供给电负荷需求的分配系数,
Figure FDA0003129988760000024
表示第j台燃气轮机产生的电能用以供给电负荷需求的分配系数;
Figure FDA0003129988760000025
Figure FDA0003129988760000026
分别表示第j台燃气轮机的气-电转化效率和气-热转化效率;
Figure FDA0003129988760000027
表示第k台中央空调的电制冷效率;
Figure FDA0003129988760000028
表示配电网传输至综合能源系统的电能总量,Pj,g表示天然气网传输至第j台燃气轮机所包含的能量大小,SDN’表示新的配电网运行状态;
其次,根据考虑预测误差后新的负荷需求量
Figure FDA0003129988760000029
以及考虑光伏预测误差和设备运行状态的实际供给量
Figure FDA00031299887600000210
计算电负荷损失值ΔLe、热负荷损失值ΔLh、冷负荷损失值ΔLc
Figure FDA00031299887600000211
2.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤1中,所述设备数据包括能量分配系数v、设备运行状态S、设备效率η和设备数量N,所述负荷数据包括能量供给量P、能量需求量L。
3.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤1构建的综合能源系统功率交换模型可用下式进行描述:
Figure FDA0003129988760000031
式中,Le、Lh、Lc分别表示电负荷需求、热负荷需求和冷负荷需求量;vDN表示从配电网传输至综合能源系统中电能用以供给电负荷需求的分配系数,
Figure FDA0003129988760000032
表示第i台光伏发电设备产生的电能用以供给电负荷需求的分配系数,
Figure FDA0003129988760000033
表示第j台燃气轮机产生的电能用以供给电负荷需求的分配系数;SDN表示配电网运行状态,
Figure FDA0003129988760000034
表示第i台光伏发电设备的运行状态,
Figure FDA0003129988760000035
表示第j台燃气轮机的运行状态,
Figure FDA0003129988760000036
表示第k台中央空调的运行状态,数值0表示故障,1表示正常运行;
Figure FDA0003129988760000037
Figure FDA0003129988760000038
Figure FDA0003129988760000039
分别表示第j台燃气轮机的气-电转化效率和气-热转化效率;
Figure FDA00031299887600000310
表示第k台中央空调的电制冷效率;
Figure FDA00031299887600000311
表示配电网传输至综合能源系统的电能总量,
Figure FDA00031299887600000312
表示第i台光伏发电设备所发电能总量,Pj,g表示天然气网传输至第j台燃气轮机所包含的能量大小。
4.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤2中,对于任意设备的运行状态,其均匀分布的累积分布函数均满足:
Figure FDA00031299887600000313
对于任意负荷和光伏发电量预测误差的累积分布函数均满足:
Figure FDA00031299887600000314
式中,预测误差的均值μ、方差σ通过统计预测误差,利用距估计或点估计法得到。
5.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤3中,采样后通过下式计算各设备采样后新的运行状态:
Figure FDA0003129988760000041
式中,
Figure FDA0003129988760000042
表示第i台光伏发电设备新的运行状态,
Figure FDA0003129988760000043
表示第j台燃气轮机新的运行状态,
Figure FDA0003129988760000044
表示第k台中央空调新的运行状态;
Figure FDA0003129988760000045
Figure FDA0003129988760000046
Figure FDA0003129988760000047
分别表示光伏发电设备、燃气轮机、中央空调发生故障的概率。
6.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤3中,考虑预测误差后,新的负荷需求量
Figure FDA0003129988760000048
和光伏发电量
Figure FDA0003129988760000049
变为:
Figure FDA00031299887600000410
式中,
Figure FDA00031299887600000411
为第i台光伏发电设备发电量的预测值,
Figure FDA00031299887600000412
Figure FDA00031299887600000413
分别为第i台光伏发电设备发电量预测误差的均值和方差;Le,f为电负荷需求的预测值,
Figure FDA00031299887600000414
Figure FDA00031299887600000415
分别为电负荷需求预测误差的均值和方差;其他L、μ、σ下标h和c同理表示为热负荷、冷负荷需求的相关变量。
7.如权利要求1所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述步骤5中具体为:
重复步骤2至步骤4共m次,m为大于等于3的自然数,得到三组数据:
Figure FDA0003129988760000051
分别计算三组数据的均值μe、μh、μc和方差σe、σh、σc,代入下式计算脆弱性指标的收敛精度θ:
Figure FDA0003129988760000052
8.如权利要求1或7所述的考虑光伏预测误差的综合能源系统脆弱性评估方法,其特征在于:所述预设条件为:脆弱性指标的收敛精度θ<0.03。
CN201910271051.0A 2019-04-04 2019-04-04 一种考虑光伏预测误差的综合能源系统脆弱性评估方法 Active CN109993445B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910271051.0A CN109993445B (zh) 2019-04-04 2019-04-04 一种考虑光伏预测误差的综合能源系统脆弱性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910271051.0A CN109993445B (zh) 2019-04-04 2019-04-04 一种考虑光伏预测误差的综合能源系统脆弱性评估方法

Publications (2)

Publication Number Publication Date
CN109993445A CN109993445A (zh) 2019-07-09
CN109993445B true CN109993445B (zh) 2021-07-30

Family

ID=67132413

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910271051.0A Active CN109993445B (zh) 2019-04-04 2019-04-04 一种考虑光伏预测误差的综合能源系统脆弱性评估方法

Country Status (1)

Country Link
CN (1) CN109993445B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109494816B (zh) * 2018-12-28 2020-07-31 清华-伯克利深圳学院筹备办公室 电-热耦合的多能流系统的风险评估方法和装置
CN110705866B (zh) * 2019-09-27 2022-04-29 新奥数能科技有限公司 评估综合能源系统能源生产设备配置与负荷匹配度的方法
CN110852631B (zh) * 2019-11-14 2022-05-24 沈阳工业大学 基于负荷预测误差的多能源系统储能容量指标计算方法
CN111857066A (zh) * 2020-06-08 2020-10-30 北京燃气能源发展有限公司 区域性能源集中管控优化方法
CN113077175A (zh) * 2021-04-21 2021-07-06 广东电网有限责任公司 一种综合能源系统运行状态评估方法
CN113420937B (zh) * 2021-07-05 2022-12-30 山西大学 一种考虑输送路径脆弱性的综合能源切换方法及系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109102132A (zh) * 2017-06-20 2018-12-28 中国电力科学研究院 一种基于概率论的微电网脆弱性评估方法及装置
CN107730129B (zh) * 2017-10-24 2021-06-22 重庆大学 考虑光热热电联产与电锅炉的电-气-热互联系统风险评估方法
CN107732970B (zh) * 2017-11-10 2020-03-17 国网甘肃省电力公司经济技术研究院 一种新能源并网电力系统的静态安全概率评估方法
CN108446865B (zh) * 2018-04-17 2019-08-16 北京清大高科系统控制有限公司 基于区间方法的热-电耦合多能流系统功率风险评估方法
CN108494021B (zh) * 2018-04-20 2021-06-01 东北大学 电-热-气综合能源系统的稳定评估与静态控制方法
CN108898287A (zh) * 2018-06-11 2018-11-27 国网江西省电力有限公司电力科学研究院 大规模光伏并网的配电网运行风险评估方法
CN109064062B (zh) * 2018-09-11 2021-08-06 浙江大学 一种考虑多能耦合交互的用户侧综合能源系统运行风险评估方法

Also Published As

Publication number Publication date
CN109993445A (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN109993445B (zh) 一种考虑光伏预测误差的综合能源系统脆弱性评估方法
Teh et al. Reliability impacts of the dynamic thermal rating and battery energy storage systems on wind-integrated power networks
CN102831321B (zh) 一种基于蒙特卡洛方法的风电场风险评估方法
He et al. Reliability modeling and evaluation of urban multi-energy systems: A review of the state of the art and future challenges
CN109377008B (zh) 电热耦合的综合能源系统风险评估方法
CN110112790B (zh) 有源配电网运行效率的评价方法、装置、设备及存储介质
CN111030094B (zh) 一种多能互补园区供能可靠性评估方法
Wang et al. Optimizing for clean-heating improvements in a district energy system with high penetration of wind power
CN103825272A (zh) 一种基于解析法的含分散风电的配电网可靠性确定方法
CN111859812B (zh) 天气影响下的海上风电场及柔直并网系统可靠性评估方法
CN112580927A (zh) 一种基于物联网的工商业园区综合能效评估方法
CN109524959B (zh) 一种计及天然气网故障的发输电系统充裕性评估方法
Meng et al. Reliability evaluation of electricity-gas-heat multi-energy consumption based on user experience
CN113610316B (zh) 不确定环境下考虑综合需求响应的园区综合能源系统优化调度方法
Bogatenkov et al. Secure state power system models building
CN113158547B (zh) 计及经济性和可靠性的区域综合能源系统优化配置方法
CN114240011A (zh) 一种综合能源系统多能流可靠性评估方法
Li et al. Reliability modeling and assessment for integrated energy system: a review of the research status and future prospects
Ge et al. Analysis method and empirical research on economic benefit of large‐scale consumptive power grid investment
AU2021209785B2 (en) Certification of at least one facility-specific amount of energy of an energy system
CN114765451A (zh) 基于故障知识图谱的光伏电站故障检测方法
Chi et al. Dynamic Reliability Assessment Framework for Integrated Energy Systems Based on the Improved Universal Generating Function
Yao et al. Static stability analysis based on probabilistic power flow calculation considering p2g technology
Chen et al. The Influence of FDIAs on Integrated Power Flow of the Integrated Energy System
Zhu et al. Reliability Evaluation of Integrated Energy System Considering Integrated Demand Response

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant