CN113792953A - 一种虚拟电厂优化调度方法及*** - Google Patents

一种虚拟电厂优化调度方法及*** Download PDF

Info

Publication number
CN113792953A
CN113792953A CN202110763812.1A CN202110763812A CN113792953A CN 113792953 A CN113792953 A CN 113792953A CN 202110763812 A CN202110763812 A CN 202110763812A CN 113792953 A CN113792953 A CN 113792953A
Authority
CN
China
Prior art keywords
load
scene
power plant
virtual power
probability
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.)
Granted
Application number
CN202110763812.1A
Other languages
English (en)
Other versions
CN113792953B (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
NARI Group Corp
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
State Grid Electric Power Research Institute
Beijing State Grid Purui UHV Transmission Technology Co Ltd
Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
NARI Group Corp
Nari Technology Co Ltd
NARI Nanjing Control System Co Ltd
State Grid Electric Power Research Institute
Beijing State Grid Purui UHV Transmission Technology Co Ltd
Economic and Technological Research Institute of State Grid Shanxi 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, NARI Group Corp, Nari Technology Co Ltd, NARI Nanjing Control System Co Ltd, State Grid Electric Power Research Institute, Beijing State Grid Purui UHV Transmission Technology Co Ltd, Economic and Technological Research Institute of State Grid Shanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110763812.1A priority Critical patent/CN113792953B/zh
Publication of CN113792953A publication Critical patent/CN113792953A/zh
Application granted granted Critical
Publication of CN113792953B publication Critical patent/CN113792953B/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/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • 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
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • 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)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种虚拟电厂优化调度方法及***,包括:根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;将上述数据输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;求解所述优化调度模型,得到虚拟电厂的优化调度方案。优点:通过将柔性负荷参与到虚拟电厂激励型需求响应中去,充分发挥虚拟电厂中分布式能源、储能***、电动汽车和柔性负荷的削峰填谷作用,减少热电互联机组的污染物的排放量。

Description

一种虚拟电厂优化调度方法及***
技术领域
本发明涉及一种虚拟电厂优化调度方法及***,属于电力***优化调度技术领域。
背景技术
面对传统电力生产带来的环境问题,以太阳能、风能和燃料电池为能源的分布式发电已逐渐被公认为是满足未来能源需求的可靠、清洁的发电方式。而随着新能源发电在电网中的渗透率的逐渐增加,新能源发电具备的不确定性和波动性对电力***带来了很大的干扰,而且由于区域内部消纳能力不足,弃风弃光现象严重。在2019年,泛在电力物联网建设概念被国家电网有限公司提出,重点通过建设虚拟电厂和多能互补等措施,提高分布式新能源的友好并网水平,促进清洁能源消纳。
虚拟电厂是指将某一区域内的发电机组、储能***、分布式能源以及可控柔性负荷有机结合在一起,通过先进的控制、计量以及通信技术,在满足区域内部的电力负荷的前提下,合并为一个整体参与电网的运行。
虚拟电厂包含区域内各种分布式能源、储能***、热电联产机组、电动汽车以及柔性负荷,虚拟电厂通过对其进行调控,获得利润,并对参加激励型需求响应负荷进行经济补偿。因此如何平衡虚拟电厂中分布式能源出力和负荷需求之间的不平衡性,整合分布式能源出力,以及用户侧柔性负荷积极参与激励型需求响应,提高虚拟电厂的经济效益,减少虚拟电厂污染物排放量,是虚拟电厂优化调度需要解决的重要问题。
发明内容
本发明所要解决的技术问题是克服现有技术的缺陷,提供一种虚拟电厂优化调度方法及***。
为解决上述技术问题,本发明提供一种虚拟电厂优化调度方法,包括:
根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
将次日风力、光伏发电所发功率场景以及概率和虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;
求解所述优化调度模型,得到虚拟电厂中的各发电资产每小时的出力情况,根据出力情况确定虚拟电厂的优化调度方案。
进一步的,所述根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率的过程包括:
通过下式建立边缘分布函数的对数似然函数,
Figure BDA0003150049800000021
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数;
通过最大化式(1)得到f(xt,δ)的边缘分布函数
Figure BDA0003150049800000022
并利用
Figure BDA0003150049800000023
求取边缘概率密度函数估计参数
Figure BDA0003150049800000024
利用最大化式(2)求取Copula函数估计参数
Figure BDA0003150049800000025
Figure BDA0003150049800000026
式中:c为Copula概率密度函数,θ表示Copula函数参数,根据求取的Copula函数估计参数代入到Student-T Copula生成风速、光照强度场景集;
利用同步回代技术,缩减Student-T Copula生成的风速、光照强度场景的数量,得到缩减后的风速、光照强度场景集;
利用下式根据缩减后的风速、光照场景集得到预测的次日的风力功率、光伏功率出力场景;
Figure BDA0003150049800000031
PPV,t=ηPVSPVεt
式中:gWPP,t为风力机组在时刻t的发电功率;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;εt为预测光伏辐射强度。
进一步的,所述利用同步回代技术,缩减场景数量的过程包括:
S1,将缩减后的风速、光照场景集设置为C,迭代次数k=1,设迭代过程的场景集为Ck=C;
S2,确定缩减后的风速、光照场景集中任意两个场景的概率距离,概率距离公式如下
Figure BDA0003150049800000032
式中:║·║2为范数表达式;i,j∈Ck,且i≠j,
Figure BDA0003150049800000033
对应任意两种情景,其对应的概率分别为πi、πj
S3,找到最接近场景i的场景,确定该场景与其他场景之间最小距离为
Figure BDA0003150049800000034
然后找到一个与这个最小距离J相匹配的场景;
S4,将每个场景的最小距离与其对应的概率相乘,得到一个与上一步中确定的最小值相匹配的场景,然后根据
Figure BDA0003150049800000035
确定场景C1,ZC1=minZi,i∈Ck
S5,消除场景C1,同时将场景C1的概率转移到最接近的概率,其余场景设定为Ck+1=Ck-C1
S6,记录剩余场景的数量,如果其余剩余场景数量满足计算需求,则继续下一步骤,如果不满足,则返回步骤S2;
S7,保留集合Ck+1的每个场景以及所对应的概率。
进一步的,所述激励型需求响应柔性负荷包括可转移负荷和可中断负荷;
其中,可转移负荷的模型为
Figure BDA0003150049800000041
式中:
Figure BDA0003150049800000042
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;
当用户侧可转移负荷参与到虚拟电厂需求响应中,虚拟电厂对转移负荷进行经济补偿,可转移负荷的经济补偿成本表达式为
Figure BDA0003150049800000043
式中:σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用;
可中断负荷用于直接消减负荷总量,可中断负荷的数学模型以及补偿费用如下式表达
Figure BDA0003150049800000044
Figure BDA0003150049800000045
式中:
Figure BDA0003150049800000046
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure BDA0003150049800000047
为第d级负荷的中断量补偿电价;nd表示电力负荷等级。
进一步的,所述优化调度模型的表达式为:
Figure BDA0003150049800000048
式中:π(ω)为场景ω发生的概率;W为总场景个数;F为虚拟电厂的净收益;
Figure BDA0003150049800000051
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure BDA0003150049800000052
为虚拟电厂在场景ω下的售热收益;
Figure BDA0003150049800000053
为热电联产机组发电发热成本;
Figure BDA0003150049800000054
为虚拟电厂运行环保代价;
Figure BDA0003150049800000055
为电动汽车运行成本;
Figure BDA0003150049800000056
为储能***运行成本;
Figure BDA0003150049800000057
为电锅炉出力成本;
Figure BDA0003150049800000058
为偏离出力计划的惩罚成本;
Figure BDA0003150049800000059
分别为转移负荷补偿成本、中断负荷补偿成本;
优化调度模型的约束条件包括:电功率平衡约束、热功率平衡约束、分布式电源出力约束、热电联产***出力与爬坡约束、储能***约束、电动汽车荷电状态约束、激励型需求响应约束。
一种虚拟电厂优化调度***,包括:
预测模块,用于根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
统计模块,用于统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
输入模块,用于将次日风力、光伏发电所发功率场景以及概率和虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;
求解模块,用于求解所述优化调度模型,得到虚拟电厂中的各发电资产每小时的出力情况,根据出力情况确定虚拟电厂的优化调度方案。
进一步的,所述预测模块包括:
第一计算模块,用于通过下式建立边缘分布函数的对数似然函数,
Figure BDA00031500498000000510
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数;
第二计算模块,用于通过最大化式(1)得到f(xt,δ)的边缘分布函数
Figure BDA0003150049800000061
并利用
Figure BDA0003150049800000062
求取边缘概率密度函数估计参数
Figure BDA0003150049800000063
第三计算模块,用于利用最大化式(2)求取Copula函数估计参数
Figure BDA0003150049800000064
Figure BDA0003150049800000065
式中:c为Copula概率密度函数,θ表示Copula函数参数,根据求取的Copula函数估计参数代入到Student-T Copula生成风速、光照强度场景集;
缩减模块,用于利用同步回代技术,缩减Student-T Copula生成的风速、光照强度场景的数量,得到缩减后的风速、光照强度场景集;
第四计算模块,用于利用下式根据缩减后的风速、光照场景集得到预测的次日的风力功率、光伏功率出力场景;
Figure BDA0003150049800000066
PPV,t=ηPVSPVεt
式中:gWPP,t为风力机组在时刻t的发电功率;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;εt为预测光伏辐射强度。
进一步的,所述缩减模块包括:
设置模块,用于将缩减后的风速、光照场景集设置为C,迭代次数k=1,设迭代过程的场景集为Ck=C;
第一确定模块,用于确定缩减后的风速、光照场景集中任意两个场景的概率距离,概率距离公式如下
Figure BDA0003150049800000067
式中:║·║2为范数表达式;i,j∈Ck,且i≠j,
Figure BDA0003150049800000068
对应任意两种情景,其对应的概率分别为πi、πj
第二确定模块,用于找到最接近场景i的场景,确定该场景与其他场景之间最小距离为
Figure BDA0003150049800000071
然后找到一个与这个最小距离J相匹配的场景;
第三确定模块,用于将每个场景的最小距离与其对应的概率相乘,得到一个与上一步中确定的最小值相匹配的场景,然后根据
Figure BDA0003150049800000072
确定场景C1,
Figure BDA0003150049800000073
消除模块,用于消除场景C1,同时将场景C1的概率转移到最接近的概率,其余场景设定为Ck+1=Ck-C1
记录模块,用于记录剩余场景的数量,如果其余剩余场景数量满足计算需求,则继续下一步骤,如果不满足,则依次通过第一确定模块、第二确定模块、第三确定模块和消除模块重新确定其余场景;
保留模块,用于保留集合Ck+1的每个场景以及所对应的概率。
进一步的,所述统计模块包括:
可转移负荷计算模块,用于计算可转移负荷:
Figure BDA0003150049800000074
式中:
Figure BDA0003150049800000075
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;
当用户侧可转移负荷参与到虚拟电厂需求响应中,虚拟电厂对转移负荷进行经济补偿,可转移负荷的经济补偿成本表达式为
Figure BDA0003150049800000076
式中:σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用;
可转移负荷计算模块,用于计算可中断负荷以及可中断负荷的补偿费用:
Figure BDA0003150049800000081
Figure BDA0003150049800000082
式中:
Figure BDA0003150049800000083
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure BDA0003150049800000084
为第d级负荷的中断量补偿电价;nd表示电力负荷等级。
进一步的,所述优化调度模型的表达式为:
Figure BDA0003150049800000085
式中:π(ω)为场景ω发生的概率;W为总场景个数;F为虚拟电厂的净收益;
Figure BDA0003150049800000086
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure BDA0003150049800000087
为虚拟电厂在场景ω下的售热收益;
Figure BDA0003150049800000088
为热电联产机组发电发热成本;
Figure BDA0003150049800000089
为虚拟电厂运行环保代价;
Figure BDA00031500498000000810
为电动汽车运行成本;
Figure BDA00031500498000000811
为储能***运行成本;
Figure BDA00031500498000000812
为电锅炉出力成本;
Figure BDA00031500498000000813
为偏离出力计划的惩罚成本;
Figure BDA00031500498000000814
分别为转移负荷补偿成本、中断负荷补偿成本;
优化调度模型的约束条件包括:电功率平衡约束、热功率平衡约束、分布式电源出力约束、热电联产***出力与爬坡约束、储能***约束、电动汽车荷电状态约束、激励型需求响应约束。
本发明所达到的有益效果:
本发明通过将柔性负荷参与到虚拟电厂激励型需求响应中去,充分发挥虚拟电厂中分布式能源、储能***、电动汽车和柔性负荷的削峰填谷作用,平缓负荷曲线,提升虚拟电厂整体运行经济,减少了热电互联机组的污染物的排放量。
附图说明
图1为本发明实例中虚拟电厂优化调度方法的流程图;
图2为本发明实例中虚拟电厂整体框架。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1-2所示,本发明提出一种虚拟电厂优化调度方法,实施流程包含如下详细步骤:
步骤1、根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
在已知历史风速,光照强度的前提下,判断历史风速、光照强度在不同时刻具有上下尾依赖性;
通过下式建立边缘分布函数的对数似然函数,
Figure BDA0003150049800000091
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数。
通过最大化上式得到f(xt)的边缘分布函数F(xt),并可以利用
Figure BDA0003150049800000092
求取边缘概率密度函数估计参数
Figure BDA0003150049800000093
利用最大化下式求取Copula函数估计参数
Figure BDA0003150049800000094
Figure BDA0003150049800000095
式中:c为Copula概率密度函数,θ表示Copula函数参数。根据求取的Copula函数参数代入到Student-T Copula生成风速、光照强度场景集。
利用下式根据上步生成的风速、光照场景得到风力功率、光伏功率出力场景;
Figure BDA0003150049800000096
PPV,t=ηPVSPVεt
式中:gWPP,t为风力发电机组在时刻t的发电;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;ε为预测光伏辐射强度;
由于经过Student-T Copula函数生成场景数据量庞大,各场景之间相似度很高,影响计算速度,因此为更有效的将相近场景合并,加快计算速度,利用同步回代技术,缩减场景数量。
步骤2、统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
在虚拟电厂中,电动汽车作为一种移动式储能,不仅可以借助充电桩从电网上获取电能,同时还可以在激励机制下向电网输送电能,实现与电网之间的V2G(Vehicle-to-Grid)互动,从而达到平抑分布式发电的不确定性以及波动性的作用。电动汽车的荷电状态模型表示为
Figure BDA0003150049800000101
式中:SOCt为t时刻电动汽车荷电状态;PEV,t为电动汽车等效出力值,大于0表示充电,小于0表示放电,等于0表示离网状态;CEV为电动汽车电池容量;λ为电动汽车单位里程耗电量;L为该电动汽车行驶里程。
柔性负荷包含可转移负荷与可中断负荷,两者均可以参与需求响应调度,根据供需情况进行转移与中断,可控程度较高。
可转移负荷是在保证整个调度周期内总负荷量不变的条件下,可以改变其使用时段的负荷。当用户侧的可转移负荷参与到需求响应中时,虚拟电厂需要对转移负荷进行经济补偿。其补偿费用不仅与转移负荷量的大小有关,还应与负荷转移的时间间隔有关,因为优化前后可转移负荷调度间隔越大,则会严重影响用户的用电舒适度。可转移负荷模型与补偿费用表示如下:
Figure BDA0003150049800000111
Figure BDA0003150049800000112
式中,
Figure BDA0003150049800000113
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用。
与可转移负荷不改变负荷总量不同,可中断负荷可以直接消减负荷总量。可中断负荷的数学模型以及补偿费用分别为
Figure BDA0003150049800000114
Figure BDA0003150049800000115
式中,
Figure BDA0003150049800000116
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure BDA0003150049800000117
为第d级负荷的中断量补偿电价,nd表示电力负荷等级。
步骤3、在满足用户电负荷需求以及热负荷需求的前提下,以最大化虚拟电厂收益为目标建立优化模型,进行虚拟电厂经济优化调度;
3.1虚拟电厂整体框架
虚拟电厂在发电侧聚合了热电联产机组、风力、光伏、储能***、电动汽车、参与电力市场交易,并与电动汽车用户、一级用户和二级用户签订双边交易协议,如附图2所示,其中一级用户由电量需求较大的工业和商业负荷组成,其用电高峰相对分散且电价高,二级用户由电量需求较小的居民负荷组成,用电高峰集中且电价低。一级用户和二级用户均由传统负荷和柔性负荷组成,传统负荷不参与虚拟电厂调度,必须优先满足,柔性负荷可以参与到激励型需求响应中,调整自身用电需求。
3.2虚拟电厂调度目标函数
虚拟电厂调度目标函数为整体利润最大,目标函数主要分为总收入与运行成本,目标函数表示如下:
Figure BDA0003150049800000121
式中:F为虚拟电厂的净收益;π(ω)为场景ω发生的概率;W为总场景个数;
Figure BDA0003150049800000122
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure BDA0003150049800000123
为虚拟电厂在场景ω下的售热收益;
Figure BDA0003150049800000124
为热电联产机组发电发热成本;
Figure BDA0003150049800000125
为虚拟电厂运行环保代价;
Figure BDA0003150049800000126
为电动汽车运行成本;
Figure BDA0003150049800000127
为储能***运行成本;
Figure BDA0003150049800000128
为电锅炉出力成本;
Figure BDA0003150049800000129
为偏离出力计划的惩罚成本;CDR1为可转移负荷总补偿费用;CDR2为中断负荷补偿成本。
3.3虚拟电厂调度约束条件
(1)电功率平衡约束
Figure BDA00031500498000001210
式中:PCHP为热电联产机组电功率出力;P′WPP、P′PV分别为风力、光伏出力;
Figure BDA00031500498000001211
Figure BDA00031500498000001212
分别为电动汽车放电、充电功率;PE1,dis、PE1,chr分别为储能***放电、充电功率;P′load1、P′load2分别为一级、二级用户电负荷需求;Peb为电转热设备消耗功率;
(2)热平衡约束
HCHP+Heb+HE2,dis=Hload+HE2,chr
式中:HCHP为热电联产机组热功率出力;Heb为电转热设备输出热功率;HE2,chr、HE2,dis分别为储热罐充、放热功率;Hload为用户热负荷需求。
(3)电动汽车约束
PEV,t|≤PEV,max
SOCmin≤SOCt≤SOCmax
SOCout,t≥SOCdemand
式中:PEV,max为电动汽车最大充放电功率;SOCmin、SOCmax为电动汽车最小、最大荷电状态;SOCdemand为电动汽车车主用车期望荷电状态。
(4)可转移负荷约束
可转移负荷量约束
Figure BDA0003150049800000131
Figure BDA0003150049800000132
式中:
Figure BDA0003150049800000133
为t时刻最大允许转入负荷量;
Figure BDA0003150049800000134
为t时刻最大允许转出负荷量。
可转移负荷转移时间约束
p(t,t′)=0 t′∈T1
p(t′,t)=0 t∈T2
式中:T1、T2分别为不允许转出时间、不允许转入时间。
转移前后负荷总量不变约束
Figure BDA0003150049800000135
(5)可中断负荷约束
PDR,d,min,t≤PDR,d,t≤PDR,d,max,t
PDR,d,t-1+PDR,d,t≤PDR,d,max
式中:PDR,d,min,t为中断负荷功率的最小值;PDR,d,max,t为中断负荷功率的最大值;PDR,d,max为第d级负荷连续时间内可调控负荷的最大响应功率。
步骤4、求解虚拟电厂优化调度模型,得到虚拟电厂的优化调度方案。
相应的本发明还提供一种虚拟电厂优化调度***,包括:
预测模块,用于根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
统计模块,用于统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
输入模块,用于将次日风力、光伏发电所发功率场景以及概率和虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;
求解模块,用于求解所述优化调度模型,得到虚拟电厂中的各发电资产每小时的出力情况,根据出力情况确定虚拟电厂的优化调度方案。
进一步的,所述预测模块包括:
第一计算模块,用于通过下式建立边缘分布函数的对数似然函数,
Figure BDA0003150049800000141
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数;
第二计算模块,用于通过最大化式(1)得到f(xt,δ)的边缘分布函数
Figure BDA0003150049800000142
并利用
Figure BDA0003150049800000143
求取边缘概率密度函数估计参数
Figure BDA0003150049800000144
第三计算模块,用于利用最大化式(2)求取Copula函数估计参数
Figure BDA0003150049800000145
Figure BDA0003150049800000146
式中:c为Copula概率密度函数,θ表示Copula函数参数,根据求取的Copula函数估计参数代入到Student-T Copula生成风速、光照强度场景集;
缩减模块,用于利用同步回代技术,缩减Student-T Copula生成的风速、光照强度场景的数量,得到缩减后的风速、光照强度场景集;
第四计算模块,用于利用下式根据缩减后的风速、光照场景集得到预测的次日的风力功率、光伏功率出力场景;
Figure BDA0003150049800000151
PPV,t=ηPVSPVεt
式中:gWPP,t为风力机组在时刻t的发电功率;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;εt为预测光伏辐射强度。
进一步的,所述缩减模块包括:
设置模块,用于将缩减后的风速、光照场景集设置为C,迭代次数k=1,设迭代过程的场景集为Ck=C;
第一确定模块,用于确定缩减后的风速、光照场景集中任意两个场景的概率距离,概率距离公式如下
Figure BDA0003150049800000152
式中:║·║2为范数表达式;i,j∈Ck,且i≠j,
Figure BDA0003150049800000153
对应任意两种情景,其对应的概率分别为πi、πj
第二确定模块,用于找到最接近场景i的场景,确定该场景与其他场景之间最小距离为
Figure BDA0003150049800000154
然后找到一个与这个最小距离J相匹配的场景;
第三确定模块,用于将每个场景的最小距离与其对应的概率相乘,得到一个与上一步中确定的最小值相匹配的场景,然后根据
Figure BDA0003150049800000155
确定场景C1,
Figure BDA0003150049800000161
消除模块,用于消除场景C1,同时将场景C1的概率转移到最接近的概率,其余场景设定为Ck+1=Ck-C1
记录模块,用于记录剩余场景的数量,如果其余剩余场景数量满足计算需求,则继续下一步骤,如果不满足,则依次通过第一确定模块、第二确定模块、第三确定模块和消除模块重新确定其余场景;
保留模块,用于保留集合Ck+1的每个场景以及所对应的概率。
进一步的,所述统计模块包括:
可转移负荷计算模块,用于计算可转移负荷:
Figure BDA0003150049800000162
式中:
Figure BDA0003150049800000163
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;
当用户侧可转移负荷参与到虚拟电厂需求响应中,虚拟电厂对转移负荷进行经济补偿,可转移负荷的经济补偿成本表达式为
Figure BDA0003150049800000164
式中:σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用;
可转移负荷计算模块,用于计算可中断负荷以及可中断负荷的补偿费用:
Figure BDA0003150049800000165
Figure BDA0003150049800000166
式中:
Figure BDA0003150049800000167
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure BDA0003150049800000171
为第d级负荷的中断量补偿电价;nd表示电力负荷等级。
进一步的,所述优化调度模型的表达式为:
Figure BDA0003150049800000172
式中:π(ω)为场景ω发生的概率;W为总场景个数;F为虚拟电厂的净收益;
Figure BDA0003150049800000173
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure BDA0003150049800000174
为虚拟电厂在场景ω下的售热收益;
Figure BDA0003150049800000175
为热电联产机组发电发热成本;
Figure BDA0003150049800000176
为虚拟电厂运行环保代价;
Figure BDA0003150049800000177
为电动汽车运行成本;
Figure BDA0003150049800000178
为储能***运行成本;
Figure BDA0003150049800000179
为电锅炉出力成本;
Figure BDA00031500498000001710
为偏离出力计划的惩罚成本;
Figure BDA00031500498000001711
分别为转移负荷补偿成本、中断负荷补偿成本;
优化调度模型的约束条件包括:电功率平衡约束、热功率平衡约束、分布式电源出力约束、热电联产***出力与爬坡约束、储能***约束、电动汽车荷电状态约束、激励型需求响应约束。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (10)

1.一种虚拟电厂优化调度方法,其特征在于,包括:
根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
将次日风力、光伏发电所发功率场景以及概率和虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;
求解所述优化调度模型,得到虚拟电厂中的各发电资产每小时的出力情况,根据出力情况确定虚拟电厂的优化调度方案。
2.根据权利要求1所述的虚拟电厂优化调度方法,其特征在于,所述根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率的过程包括:
通过下式建立边缘分布函数的对数似然函数,
Figure FDA0003150049790000011
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数;
通过最大化式(1)得到f(xt,δ)的边缘分布函数
Figure FDA0003150049790000012
并利用
Figure FDA0003150049790000013
求取边缘概率密度函数估计参数
Figure FDA0003150049790000014
利用最大化式(2)求取Copula函数估计参数
Figure FDA0003150049790000015
Figure FDA0003150049790000016
式中:c为Copula概率密度函数,θ表示Copula函数参数,根据求取的Copula函数估计参数代入到Student-T Copula生成风速、光照强度场景集;
利用同步回代技术,缩减Student-T Copula生成的风速、光照强度场景的数量,得到缩减后的风速、光照强度场景集;
利用下式根据缩减后的风速、光照场景集得到预测的次日的风力功率、光伏功率出力场景;
Figure FDA0003150049790000021
PPV,t=ηPVSPVεt
式中:gWPP,t为风力机组在时刻t的发电功率;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;εt为预测光伏辐射强度。
3.根据权利要求2所述的虚拟电厂优化调度方法,其特征在于,所述利用同步回代技术,缩减场景数量的过程包括:
S1,将缩减后的风速、光照场景集设置为C,迭代次数k=1,设迭代过程的场景集为Ck=C;
S2,确定缩减后的风速、光照场景集中任意两个场景的概率距离,概率距离公式如下
Figure FDA0003150049790000022
式中:║·║2为范数表达式;i,j∈Ck,且i≠j,
Figure FDA0003150049790000023
对应任意两种情景,其对应的概率分别为πi、πj
S3,找到最接近场景i的场景,确定该场景与其他场景之间最小距离为
Figure FDA0003150049790000024
然后找到一个与这个最小距离J相匹配的场景;
S4,将每个场景的最小距离与其对应的概率相乘,得到一个与上一步中确定的最小值相匹配的场景,然后根据
Figure FDA0003150049790000031
确定场景C1,
Figure FDA0003150049790000032
S5,消除场景C1,同时将场景C1的概率转移到最接近的概率,其余场景设定为Ck+1=Ck-C1
S6,记录剩余场景的数量,如果其余剩余场景数量满足计算需求,则继续下一步骤,如果不满足,则返回步骤S2;
S7,保留集合Ck+1的每个场景以及所对应的概率。
4.根据权利要求3所述的虚拟电厂优化调度方法,其特征在于,所述激励型需求响应柔性负荷包括可转移负荷和可中断负荷;
其中,可转移负荷的模型为
Figure FDA0003150049790000033
式中:
Figure FDA0003150049790000034
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;
当用户侧可转移负荷参与到虚拟电厂需求响应中,虚拟电厂对转移负荷进行经济补偿,可转移负荷的经济补偿成本表达式为
Figure FDA0003150049790000035
式中:σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用;
可中断负荷用于直接消减负荷总量,可中断负荷的数学模型以及补偿费用如下式表达
Figure FDA0003150049790000036
Figure FDA0003150049790000041
式中:
Figure FDA0003150049790000042
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure FDA0003150049790000043
为第d级负荷的中断量补偿电价;nd表示电力负荷等级。
5.根据权利要求4所述的虚拟电厂优化调度方法,其特征在于,所述优化调度模型的表达式为:
Figure FDA0003150049790000044
式中:π(ω)为场景ω发生的概率;W为总场景个数;F为虚拟电厂的净收益;
Figure FDA0003150049790000045
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure FDA0003150049790000046
为虚拟电厂在场景ω下的售热收益
Figure FDA0003150049790000047
为热电联产机组发电发热成本;
Figure FDA0003150049790000048
为虚拟电厂运行环保代价;
Figure FDA0003150049790000049
为电动汽车运行成本;
Figure FDA00031500497900000410
为储能***运行成本;
Figure FDA00031500497900000411
为电锅炉出力成本;
Figure FDA00031500497900000412
为偏离出力计划的惩罚成本;
Figure FDA00031500497900000413
分别为转移负荷补偿成本、中断负荷补偿成本;
优化调度模型的约束条件包括:电功率平衡约束、热功率平衡约束、分布式电源出力约束、热电联产***出力与爬坡约束、储能***约束、电动汽车荷电状态约束、激励型需求响应约束。
6.一种虚拟电厂优化调度***,其特征在于,包括:
预测模块,用于根据虚拟电厂内部可再生分布式能源所在区域风速、温度、光照辐射强度历史信息,预测次日风力、光伏发电所发功率场景以及概率;
统计模块,用于统计虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量;
输入模块,用于将次日风力、光伏发电所发功率场景以及概率和虚拟电厂内部电动汽车数量、储能***容量、热电联产机组数量以及参与激励型需求响应柔性负荷数量输入到预先构建的满足用户电负荷需求以及热负荷需求的,以最大化虚拟电厂收益为目标的优化调度模型;
求解模块,用于求解所述优化调度模型,得到虚拟电厂中的各发电资产每小时的出力情况,根据出力情况确定虚拟电厂的优化调度方案。
7.根据权利要求6所述的虚拟电厂优化调度***,其特征在于,所述预测模块包括:
第一计算模块,用于通过下式建立边缘分布函数的对数似然函数,
Figure FDA0003150049790000051
式中:L(δ)为对数似然函数,f(xt,δ)为边缘概率密度函数,xt为t时刻的风速/光照强度,δ为边缘概率密度函数参数,T表示为一天总小时数;
第二计算模块,用于通过最大化式(1)得到f(xt,δ)的边缘分布函数
Figure FDA0003150049790000052
并利用
Figure FDA0003150049790000053
求取边缘概率密度函数估计参数
Figure FDA0003150049790000054
第三计算模块,用于利用最大化式(2)求取Copula函数估计参数
Figure FDA0003150049790000055
Figure FDA0003150049790000056
式中:c为Copula概率密度函数,θ表示Copula函数参数,根据求取的Copula函数估计参数代入到Student-T Copula生成风速、光照强度场景集;
缩减模块,用于利用同步回代技术,缩减Student-T Copula生成的风速、光照强度场景的数量,得到缩减后的风速、光照强度场景集;
第四计算模块,用于利用下式根据缩减后的风速、光照场景集得到预测的次日的风力功率、光伏功率出力场景;
Figure FDA0003150049790000057
PPV,t=ηPVSPVεt
式中:gWPP,t为风力机组在时刻t的发电功率;vt为时刻t的预测风速;vin和vout为切入风速和切出风速;vr为额定风速;gr为额定输出功率;PPV,t为光伏发电机组在时刻t的出力;ηPV为光伏转换效率;SPV为光伏面积;εt为预测光伏辐射强度。
8.根据权利要求7所述的虚拟电厂优化调度***,其特征在于,所述缩减模块包括:
设置模块,用于将缩减后的风速、光照场景集设置为C,迭代次数k=1,设迭代过程的场景集为Ck=C;
第一确定模块,用于确定缩减后的风速、光照场景集中任意两个场景的概率距离,概率距离公式如下
Figure FDA0003150049790000061
式中:║·║2为范数表达式;i,j∈Ck,且i≠j,
Figure FDA0003150049790000062
对应任意两种情景,其对应的概率分别为πi、πj
第二确定模块,用于找到最接近场景i的场景,确定该场景与其他场景之间最小距离为
Figure FDA0003150049790000063
然后找到一个与这个最小距离J相匹配的场景;
第三确定模块,用于将每个场景的最小距离与其对应的概率相乘,得到一个与上一步中确定的最小值相匹配的场景,然后根据
Figure FDA0003150049790000064
确定场景C1,
Figure FDA0003150049790000065
消除模块,用于消除场景C1,同时将场景C1的概率转移到最接近的概率,其余场景设定为Ck+1=Ck-C1
记录模块,用于记录剩余场景的数量,如果其余剩余场景数量满足计算需求,则继续下一步骤,如果不满足,则依次通过第一确定模块、第二确定模块、第三确定模块和消除模块重新确定其余场景;
保留模块,用于保留集合Ck+1的每个场景以及所对应的概率。
9.根据权利要求8所述的虚拟电厂优化调度***,其特征在于,所述统计模块包括:
可转移负荷计算模块,用于计算可转移负荷:
Figure FDA0003150049790000071
式中:
Figure FDA0003150049790000072
分别为d级负荷在时刻t的转移负荷前、后的负荷;p(t,t′)为从t时刻转移到t′时刻的负荷量;
当用户侧可转移负荷参与到虚拟电厂需求响应中,虚拟电厂对转移负荷进行经济补偿,可转移负荷的经济补偿成本表达式为
Figure FDA0003150049790000073
式中:σ1、σ2为可转移负荷补偿费用的基本补偿费用与用户舒适度补偿费用;Tmax为最大转移时间间隔;CDR1为可转移负荷总补偿费用;
可转移负荷计算模块,用于计算可中断负荷以及可中断负荷的补偿费用:
Figure FDA0003150049790000074
Figure FDA0003150049790000075
式中:
Figure FDA0003150049790000076
P′load,d,t分别为d级负荷中断负荷前、后的负荷总量;CDR2,d为中断负荷补偿成本;PDR,d,t为第d级负荷中断量;
Figure FDA0003150049790000077
为第d级负荷的中断量补偿电价;nd表示电力负荷等级。
10.根据权利要求9所述的虚拟电厂优化调度***,其特征在于,所述优化调度模型的表达式为:
Figure FDA0003150049790000078
式中:π(ω)为场景ω发生的概率;W为总场景个数;F为虚拟电厂的净收益;
Figure FDA0003150049790000079
为虚拟电厂在场景ω下为总负荷供电的收益;
Figure FDA00031500497900000710
为虚拟电厂在场景ω下的售热收益
Figure FDA00031500497900000711
为热电联产机组发电发热成本;
Figure FDA00031500497900000712
为虚拟电厂运行环保代价;
Figure FDA00031500497900000713
为电动汽车运行成本;
Figure FDA00031500497900000714
为储能***运行成本;
Figure FDA00031500497900000715
为电锅炉出力成本;
Figure FDA00031500497900000716
为偏离出力计划的惩罚成本;
Figure FDA00031500497900000717
分别为转移负荷补偿成本、中断负荷补偿成本;
优化调度模型的约束条件包括:电功率平衡约束、热功率平衡约束、分布式电源出力约束、热电联产***出力与爬坡约束、储能***约束、电动汽车荷电状态约束、激励型需求响应约束。
CN202110763812.1A 2021-07-06 2021-07-06 一种虚拟电厂优化调度方法及*** Active CN113792953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110763812.1A CN113792953B (zh) 2021-07-06 2021-07-06 一种虚拟电厂优化调度方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110763812.1A CN113792953B (zh) 2021-07-06 2021-07-06 一种虚拟电厂优化调度方法及***

Publications (2)

Publication Number Publication Date
CN113792953A true CN113792953A (zh) 2021-12-14
CN113792953B CN113792953B (zh) 2024-03-19

Family

ID=79180992

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110763812.1A Active CN113792953B (zh) 2021-07-06 2021-07-06 一种虚拟电厂优化调度方法及***

Country Status (1)

Country Link
CN (1) CN113792953B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114626603A (zh) * 2022-03-11 2022-06-14 国网冀北电力有限公司承德供电公司 一种配电网需求响应调度运行方案的优化处理方法和装置
CN115358519A (zh) * 2022-07-13 2022-11-18 上海嘉柒智能科技有限公司 一种虚拟电厂优化调度方法及装置
CN116031886A (zh) * 2023-01-30 2023-04-28 天津大学 热泵虚拟电厂日前-日内两阶段灵活性爬坡能力控制方法
CN116031954A (zh) * 2023-02-21 2023-04-28 国网浙江省电力有限公司金华供电公司 模糊机会约束下的虚拟电厂优化调度方法
CN116562575A (zh) * 2023-05-16 2023-08-08 中国电力工程顾问集团有限公司 一种综合能源***的优化调度方法
CN117094453A (zh) * 2023-10-20 2023-11-21 国网安徽省电力有限公司合肥供电公司 一种用于虚拟电厂的调度优化***及其方法
CN117791627A (zh) * 2024-02-26 2024-03-29 国网山东省电力公司东营供电公司 考虑虚拟电厂不确定性的柔性负荷动态聚合方法及***
CN117977599A (zh) * 2024-03-28 2024-05-03 国网浙江新兴科技有限公司 含有电动汽车的虚拟电厂的调度方法、***、设备和介质
CN118096224A (zh) * 2024-04-24 2024-05-28 国网安徽省电力有限公司经济技术研究院 基于虚拟电厂优化运行的电力市场出清方法及***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150117085A (ko) * 2014-04-09 2015-10-19 전자부품연구원 가상 발전 관리 플랫폼 시스템 및 에너지 관리 방법
CN109523052A (zh) * 2018-09-18 2019-03-26 国网浙江省电力有限公司经济技术研究院 一种考虑需求响应和碳交易的虚拟电厂优化调度方法
KR20190105173A (ko) * 2018-02-22 2019-09-16 (주)솔라라이트 에너지 저장 장치를 이용한 가상 발전소
CN110516851A (zh) * 2019-08-05 2019-11-29 华北电力大学 一种基于虚拟电厂的源荷双侧热电联合随机优化调度方法
CN112072640A (zh) * 2020-08-13 2020-12-11 清华大学 一种虚拟电厂聚合资源的容量优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20150117085A (ko) * 2014-04-09 2015-10-19 전자부품연구원 가상 발전 관리 플랫폼 시스템 및 에너지 관리 방법
KR20190105173A (ko) * 2018-02-22 2019-09-16 (주)솔라라이트 에너지 저장 장치를 이용한 가상 발전소
CN109523052A (zh) * 2018-09-18 2019-03-26 国网浙江省电力有限公司经济技术研究院 一种考虑需求响应和碳交易的虚拟电厂优化调度方法
CN110516851A (zh) * 2019-08-05 2019-11-29 华北电力大学 一种基于虚拟电厂的源荷双侧热电联合随机优化调度方法
CN112072640A (zh) * 2020-08-13 2020-12-11 清华大学 一种虚拟电厂聚合资源的容量优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YI WANG ETAL: ""Control of PMSG-Based Wind Turbines for System Inertial Response and Power Oscillation Damping"", 《CONTROL OF PMSG-BASED WIND TURBINES FOR SYSTEM INERTIAL RESPONSE AND POWER OSCILLATION DAMPING》, vol. 6, no. 2, pages 565 - 574, XP011576028, DOI: 10.1109/TSTE.2015.2394363 *
王天旺等: ""虚拟电厂下计及分布式风电与储能***的电力***优化调度_王天旺"", 《电力建设》, vol. 37, no. 11, pages 108 - 114 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114626603A (zh) * 2022-03-11 2022-06-14 国网冀北电力有限公司承德供电公司 一种配电网需求响应调度运行方案的优化处理方法和装置
CN115358519A (zh) * 2022-07-13 2022-11-18 上海嘉柒智能科技有限公司 一种虚拟电厂优化调度方法及装置
CN116031886B (zh) * 2023-01-30 2024-02-13 天津大学 热泵虚拟电厂日前-日内两阶段灵活性爬坡能力控制方法
CN116031886A (zh) * 2023-01-30 2023-04-28 天津大学 热泵虚拟电厂日前-日内两阶段灵活性爬坡能力控制方法
CN116031954A (zh) * 2023-02-21 2023-04-28 国网浙江省电力有限公司金华供电公司 模糊机会约束下的虚拟电厂优化调度方法
CN116562575A (zh) * 2023-05-16 2023-08-08 中国电力工程顾问集团有限公司 一种综合能源***的优化调度方法
CN116562575B (zh) * 2023-05-16 2023-10-31 中国电力工程顾问集团有限公司 一种综合能源***的优化调度方法
CN117094453A (zh) * 2023-10-20 2023-11-21 国网安徽省电力有限公司合肥供电公司 一种用于虚拟电厂的调度优化***及其方法
CN117094453B (zh) * 2023-10-20 2024-01-30 国网安徽省电力有限公司合肥供电公司 一种用于虚拟电厂的调度优化***及其方法
CN117791627A (zh) * 2024-02-26 2024-03-29 国网山东省电力公司东营供电公司 考虑虚拟电厂不确定性的柔性负荷动态聚合方法及***
CN117791627B (zh) * 2024-02-26 2024-05-14 国网山东省电力公司东营供电公司 考虑虚拟电厂不确定性的柔性负荷动态聚合方法及***
CN117977599A (zh) * 2024-03-28 2024-05-03 国网浙江新兴科技有限公司 含有电动汽车的虚拟电厂的调度方法、***、设备和介质
CN117977599B (zh) * 2024-03-28 2024-06-11 国网浙江新兴科技有限公司 含有电动汽车的虚拟电厂的调度方法、***、设备和介质
CN118096224A (zh) * 2024-04-24 2024-05-28 国网安徽省电力有限公司经济技术研究院 基于虚拟电厂优化运行的电力市场出清方法及***
CN118096224B (zh) * 2024-04-24 2024-08-13 国网安徽省电力有限公司经济技术研究院 基于虚拟电厂优化运行的电力市场出清方法及***

Also Published As

Publication number Publication date
CN113792953B (zh) 2024-03-19

Similar Documents

Publication Publication Date Title
CN113792953B (zh) 一种虚拟电厂优化调度方法及***
Hou et al. Multi-objective economic dispatch of a microgrid considering electric vehicle and transferable load
Wang et al. Multi-stage optimal energy management of multi-energy microgrid in deregulated electricity markets
Moradi et al. Operational strategy optimization in an optimal sized smart microgrid
CN108206543B (zh) 一种基于能源梯级利用的能源路由器及其运行优化方法
CN111340274A (zh) 一种基于虚拟电厂参与的综合能源***优化方法和***
Ju et al. A Tri-dimensional Equilibrium-based stochastic optimal dispatching model for a novel virtual power plant incorporating carbon Capture, Power-to-Gas and electric vehicle aggregator
Wang et al. Synergistic planning of an integrated energy system containing hydrogen storage with the coupled use of electric-thermal energy
CN103151797A (zh) 基于多目标调度模型的并网运行方式下微网能量控制方法
CN110957717A (zh) 一种多电源电力***多目标日前优化调度方法
Wang et al. Energy management strategy of hybrid energy storage based on Pareto optimality
CN116316567A (zh) 一种阶梯碳交易机制下综合能源需求响应优化方法
Zhao et al. Day-ahead robust optimal dispatch of integrated energy station considering battery exchange service
CN111293718A (zh) 基于场景分析的交直流混合微网分区二层优化运行方法
CN115759610A (zh) 一种电力***源网荷储协同的多目标规划方法及其应用
CN111030101B (zh) 一种基于多元化大数据清洁能源消纳联动调控方法及***
CN115271264A (zh) 一种工业园区能源***调配方法及计算设备
CN116167573A (zh) 一种高铁站综合能源需求响应优化调度方法及相关装置
CN115689166A (zh) 一种区域分布式能源资源聚合利用方法及***
CN111126675A (zh) 多能互补微网***优化方法
CN112541612B (zh) 一种基于多网耦合的虚拟电厂经济优化调度方法
CN112883630A (zh) 用于风电消纳的多微网***日前优化经济调度方法
Mahoor et al. Smart energy management for a micro-grid with consideration of demand response plans
CN117039862A (zh) 基于氢共享储能服务的微网群双层优化配置方法和装置
CN116128101A (zh) 虚拟电厂优化调度模型构建方法

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