[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN109523052B - 一种考虑需求响应和碳交易的虚拟电厂优化调度方法 - Google Patents

一种考虑需求响应和碳交易的虚拟电厂优化调度方法 Download PDF

Info

Publication number
CN109523052B
CN109523052B CN201811088161.5A CN201811088161A CN109523052B CN 109523052 B CN109523052 B CN 109523052B CN 201811088161 A CN201811088161 A CN 201811088161A CN 109523052 B CN109523052 B CN 109523052B
Authority
CN
China
Prior art keywords
representing
power plant
time
virtual power
electric
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
CN201811088161.5A
Other languages
English (en)
Other versions
CN109523052A (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.)
Hangzhou Guodian Electric Power Technology Development Co ltd
Economic and Technological Research Institute of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Hangzhou Guodian Electric Power Technology Development Co ltd
Economic and Technological Research Institute of State Grid Zhejiang 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 Hangzhou Guodian Electric Power Technology Development Co ltd, Economic and Technological Research Institute of State Grid Zhejiang Electric Power Co Ltd filed Critical Hangzhou Guodian Electric Power Technology Development Co ltd
Priority to CN201811088161.5A priority Critical patent/CN109523052B/zh
Publication of CN109523052A publication Critical patent/CN109523052A/zh
Application granted granted Critical
Publication of CN109523052B publication Critical patent/CN109523052B/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/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/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
    • 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems
    • 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)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (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

一种考虑需求响应和碳交易的虚拟电厂优化调度方法
技术领域
本发明属于电力系统优化调度领域,特别是一种考虑需求响应和碳交易的虚拟电厂优化调度方法。
背景技术
为解决化石能源紧缺和化石燃料燃烧带来的环境污染问题,清洁可再生能源发电机组得到了越来越广泛的应用;为缓解电力系统负荷侧与日俱增的需求压力,可控负荷的应用逐渐受到重视。可再生能源发电机组和可控负荷的不断增加,诸如发电机组出力的波动性、源荷功率平衡的不确定性等新问题也逐渐出现,为电力系统的安全经济运行带来了新的挑战。基于此,“虚拟电厂”的概念应运而生。
虚拟电厂通过先进的通信手段和软件系统,将燃气轮机、间歇性可再生能源发电机组、可控负荷等分布式能源聚合在一起,形成一个整体参与到电力系统运行和电力市场运营之中。目前,国内外已有一些虚拟电厂项目,比如:国电云南小中甸风光水分布式电源示范项目、欧洲FENIX虚拟电厂项目等。
现有关于虚拟电厂优化调度已展开了一定的研究,主要集中于虚拟电厂利用储能平抑间歇性可再生能源发电机组出力的波动性、参与包含日前市场和实时市场等多阶段电力市场以获取经济收益方面。除此之外,虚拟电厂仍有其他优良特性可以在优化调度中发挥出来。虚拟电厂中含有的清洁可再生能源发电机组和电动汽车,由于其优良的环保特性,可以有效减少电力系统的碳排放,带来显著的环境效益;同时,虚拟电厂中的可控负荷,可以作为需求响应资源,为用户带来经济收益的同时优化电力系统负荷曲线,提升电力系统的经济性与安全性。目前,对于虚拟电厂在优化调度中同时考虑需求响应和碳减排效益方面,鲜有有效的科学研究。
发明内容
本发明所要解决的技术问题是克服上述现有技术存在的缺陷,提供一种考虑需求响应和碳交易的虚拟电厂优化调度方法,以提升电力系统的减排效益,提升虚拟电厂的经济收益以及电力系统的运行安全性。
为此,本发明采用以下的技术方案:一种考虑需求响应和碳交易的虚拟电厂优化调度方法,其包括步骤:
1)统计虚拟电厂内分布式电源的数量,以及电动汽车、水蓄冷空调系统参与需求响应的可控负荷数量,所述的分布式电源包括燃气轮机、风电机组和光伏发电机组;
2)对次日虚拟电厂内风电、光伏机组的出力和电力系统内刚性负荷进行预测;
3)虚拟电厂的能量管理系统接受电动汽车车主的日前申报信息,包括各电动汽车的入网及离网时刻、预估的初始电池荷电状态和期望的离网荷电状态;
4)虚拟电厂的能量管理系统根据日前预测信息、电动汽车申报信息以及电力市场和碳交易市场的价格信息,建立虚拟电厂同时参与电力市场和碳交易市场的优化调度模型;
5)求解虚拟电厂的优化调度模型,得到虚拟电厂的日前优化调度方案。
作为上述技术方案的补充,虚拟电厂优化调度模型的目标函数为虚拟电厂所获得的经济收益最大,其数学表达式如下所示:
Figure BDA0001803679040000021
Figure BDA0001803679040000022
Figure BDA0001803679040000023
Figure BDA0001803679040000024
Figure BDA0001803679040000025
Figure BDA0001803679040000026
式中:T表示一天内的时段数;
Figure BDA0001803679040000027
表示在时刻t的分时电价;Pt s表示在时刻t虚拟电厂与配电系统的交互电功率;
Figure BDA0001803679040000028
Figure BDA0001803679040000029
分别表示在时刻t的燃气轮机运行成本、碳交易成本和电动汽车充放电成本;NG表示虚拟电厂内燃气轮机的数量;kg,0表示燃气轮机g的固定运行成本;ug,t表示在时刻t燃气轮机g的开关状态变量;nj表示燃气轮机运行成本的线性化段数;kg,j表示第j段线性化运行成本函数的斜率;Pg,j,t表示在时刻t燃气轮机g的第j段出力大小;Pg,t表示在时刻t燃气轮机g的总出力大小;
Figure BDA00018036790400000210
Figure BDA00018036790400000211
分别表示在时刻t燃气轮机g的启动和停机成本;
Figure BDA00018036790400000212
Figure BDA00018036790400000213
表示燃气轮机g启动和停机一次所消耗的成本;Nv表示虚拟电厂内电动汽车的数量;
Figure BDA00018036790400000214
Figure BDA00018036790400000215
分别表示电动汽车v在时刻t的充放电功率;
Figure BDA0001803679040000031
Figure BDA0001803679040000032
分别表示在时刻t的充放电电价。
作为上述技术方案的补充,虚拟电厂参与的碳交易机制表示为:
Figure BDA0001803679040000033
Figure BDA0001803679040000034
ΔMt=MC,t-MD,t
Figure BDA0001803679040000035
式中:N表示虚拟电厂中所有发电机组的数量;ε表示虚拟电厂的单位电量碳排放分配系数;Pn,t表示在时刻t发电机组n的有功出力;KC表示单位碳排放量交易价格;MC,t表示在时刻t虚拟电厂的实际碳排放量;MD,t表示在时刻t虚拟电厂分配到的碳排放配额;Qg表示燃气轮机g的单位电量碳排放强度。
如果在时刻t虚拟电厂的实际碳排放量低于分配得到的碳排放配额,虚拟电厂通过将多余的碳排放配额出售至碳交易市场以获取相应的利润;否则,虚拟电厂必须在碳交易市场中购买相应的碳排放权。
作为上述技术方案的补充,作为需求响应资源的电动汽车充放电负荷需满足的约束条件如下:
Figure BDA0001803679040000036
Figure BDA0001803679040000037
Figure BDA0001803679040000038
Sv,min≤Sv,t≤Sv,max
Figure BDA0001803679040000039
Figure BDA00018036790400000310
Figure BDA00018036790400000311
式中:
Figure BDA00018036790400000312
Figure BDA00018036790400000313
分别表示电动汽车v在时刻t的充电和放电状态变量;
Figure BDA00018036790400000314
Figure BDA00018036790400000315
分别表示电动汽车v的额定充放电功率;
Figure BDA00018036790400000316
表示电动汽车v在时刻t的调度状态;Sv,t表示电动汽车v在时刻t的电池荷电状态;Sv,min和Sv,max分别表示电动汽车v的最小和最大荷电状态;
Figure BDA0001803679040000041
Figure BDA0001803679040000042
分别表示电动汽车v的入网和离网时刻;
Figure BDA0001803679040000043
Figure BDA0001803679040000044
分别表示电动汽车v在入网时刻和离网时刻对应的电池荷电状态;Sv,a表示电动汽车v在入网时刻的电池荷电状态;Sv,d表示电动汽车v在离网时刻的目标电池荷电状态;
Figure BDA0001803679040000045
Figure BDA0001803679040000046
分别表示电动汽车v的充放电效率;Emax表示电动汽车电池容量;Δt表示时段长度。
作为上述技术方案的补充,作为需求响应资源的水蓄冷空调负荷建模如下:
1)建筑物热动态平衡方程
Figure BDA0001803679040000047
式中:θout,t和θin,t分别表示在时刻t的户外和室内温度;B是建筑物的热传导系数;Ca表示空气的比热容系数;V表示室内容积;ρa表示空气密度;Qt表示在时刻t建筑物的瞬时得热量;CD,t表示在时刻t建筑物内制冷装置提供的冷量;Δt表示时段长度;
2)热舒适度建模
采用PMV方法衡量建筑物内的热舒适度指标,PMV值和温度之间的关系为:
Figure BDA0001803679040000048
式中:IPMV表示室内的PMV值;θ表示温度,℃;
3)水蓄冷空调系统建模
单位时间内水蓄冷空调系统提供的冷量CD,i,t为:
CD,i,t=Cp,i,t-Cs,i,t+Cr,i,t
式中:Cp,i,t,Cs,i,t和Cr,i,t分别表示在时刻t水蓄冷空调i的制冷量、蓄冷量和放冷量;
水蓄冷空调需满足的运行约束为:
0≤Is,i,t+Ir,i,t≤1,
Figure BDA0001803679040000049
Figure BDA00018036790400000410
Figure BDA00018036790400000411
Figure BDA00018036790400000412
Figure BDA00018036790400000413
式中:Is,i,t和Ir,i,t分别表示在时刻t水蓄冷空调i的蓄冷和放冷运行状态变量;
Figure BDA0001803679040000051
Figure BDA0001803679040000052
分别表示在单位时间内水蓄冷空调的最大制冷量、蓄冷量和放冷量;
Figure BDA0001803679040000053
表示在时刻t水蓄冷空调i储存的冷量大小;Sc,max表示水蓄冷空调的蓄冷水箱容量;ηs和ηr分别表示蓄冷效率和放冷效率;Δt表示时段长度;
水蓄冷空调i在时刻t消耗的电功率Pi,t为:
Figure BDA0001803679040000054
式中:μp表示能耗比,用来描述水蓄冷空调的制冷机的制冷效果和消耗的电功率之间的关系;μs和μr分别表示蓄冷量和放冷量与消耗电能的转换系数。
作为上述技术方案的补充,需满足的燃气轮机运行约束和电力系统能量平衡约束如下:
Figure BDA0001803679040000055
Figure BDA0001803679040000056
Figure BDA0001803679040000057
式中:
Figure BDA0001803679040000058
Figure BDA0001803679040000059
分别表示燃气轮机g的最小和最大输出有功功率;
Figure BDA00018036790400000510
Figure BDA00018036790400000511
分别表示燃气轮机g的爬坡和滑坡功率;Δt表示时段长度;Nq表示虚拟电厂内间歇性可再生能源发电机组(主要指风电机组和光伏发电机组)的数量;Pk,t表示可再生能源发电机组k在时刻t的有功出力;Nc表示虚拟电厂内水蓄冷空调系统的数量;Pt R表示虚拟电厂在时刻t的刚性负荷需求。
本发明具有的有益效果如下:
本发明提供的考虑需求响应和碳交易的虚拟电厂优化调度方法,通过将碳交易机制引入到虚拟电厂优化调度中,可以充分发挥虚拟电厂中大量风电、光伏发电机组和低碳排放量的燃气轮机的环保优势,提升电力系统的减排效益;水蓄冷空调系统和电动汽车通过参与需求响应,提升虚拟电厂的经济收益以及电力系统的运行安全性。
附图说明
图1为本发明实施例中考虑需求响应和碳交易的虚拟电厂优化调度方法的流程图;
图2为本发明应用例中一个典型水蓄冷空调系统的电功率和蓄冷量随时间变化情况图;
图3为本发明应用例中电动汽车充放电负荷随时间变化情况图。
具体实施方式
为更好地理解本发明的目的、技术方案以及技术效果,以下结合附图对本发明进行进一步的说明。
实施例
本发明提出了一种考虑需求响应和碳交易的虚拟电厂优化调度方法,实施流程包括如下详细步骤:
步骤1、统计虚拟电厂内分布式电源(如燃气轮机、风电机组、光伏发电机组)的数量,以及电动汽车、水蓄冷空调系统等可参与需求响应的可控负荷数量;
步骤2、对次日虚拟电厂内风电、光伏机组出力和电力系统内刚性负荷进行预测;
步骤3、虚拟电厂的能量管理系统接受电动汽车车主的日前申报信息,具体包括:各电动汽车的入网及离网时刻、预估的初始电池荷电状态和期望的离网荷电状态;
步骤4、建立虚拟电厂同时考虑需求响应和碳交易的优化调度模型;
具体地,以虚拟电厂所获得的经济收益最大为目标函数建立优化模型,其数学表达式如下所示:
Figure BDA0001803679040000061
Figure BDA0001803679040000062
Figure BDA0001803679040000063
Figure BDA0001803679040000064
Figure BDA0001803679040000065
Figure BDA0001803679040000066
式中:T表示一天内的时段数;
Figure BDA0001803679040000067
表示在时刻t的分时电价;Pt s表示在时刻t虚拟电厂与配电系统的交互电功率;
Figure BDA0001803679040000068
Figure BDA0001803679040000069
分别表示在时刻t的燃气轮机运行成本、碳交易成本和电动汽车充放电成本;NG表示虚拟电厂内燃气轮机的数量;kg,0表示燃气轮机g的固定运行成本;ug,t表示在时刻t燃气轮机g的开关状态变量;nj表示燃气轮机运行成本的线性化段数;kg,j表示第j段线性运行成本函数的斜率;Pg,j,t表示在时刻t燃气轮机g的第j段出力大小;Pg,t表示在时刻t燃气轮机g的总出力大小;
Figure BDA0001803679040000071
Figure BDA0001803679040000072
分别表示在时刻t燃气轮机g的启动和停机成本;
Figure BDA0001803679040000073
表示燃气轮机g启动和停机一次所消耗的成本;Nv表示虚拟电厂内电动汽车的数量;
Figure BDA0001803679040000074
Figure BDA0001803679040000075
分别表示电动汽车v在时刻t的充放电功率;
Figure BDA0001803679040000076
分别表示在时刻t的充放电电价。
虚拟电厂参与的碳交易机制,可表示为:
Figure BDA0001803679040000077
Figure BDA0001803679040000078
Figure BDA0001803679040000079
式中:N表示虚拟电厂中所有发电机组的数量;ε表示虚拟电厂的单位电量碳排放分配系数,由国家发改委发布的“区域电网基准线排放因子”确定;KC表示单位碳排放量交易价格;MC,t表示在时刻t虚拟电厂的实际碳排放量;MD,t表示在时刻t虚拟电厂分配到的碳排放配额;Qg表示燃气轮机g的单位电量碳排放强度;Pn,t表示在时刻t发电机组n的有功出力。
如果在时刻t虚拟电厂的实际碳排放量低于分配得到的碳排放配额,虚拟电厂可以通过将多余的碳排放配额出售至碳交易市场以获取相应的利润;否则,虚拟电厂必须在碳交易市场中购买相应的碳排放权。
电动汽车和水蓄冷空调负荷,是虚拟电厂中典型的需求响应资源。其中,电动汽车充放电负荷参与需求响应,需满足如下约束条件:
Figure BDA00018036790400000710
Figure BDA00018036790400000711
Figure BDA00018036790400000712
Sv,min≤Sv,t≤Sv,max
Figure BDA00018036790400000713
Figure BDA00018036790400000714
Figure BDA0001803679040000081
式中:
Figure BDA0001803679040000082
Figure BDA0001803679040000083
分别表示电动汽车v在时刻t的充电和放电状态变量;
Figure BDA0001803679040000084
Figure BDA0001803679040000085
分别表示电动汽车v的额定充放电功率;
Figure BDA0001803679040000086
表示电动汽车v在时刻t的调度状态;Sv,t表示电动汽车v在时刻t的电池荷电状态;Sv,min和Sv,max分别表示电动汽车v的最小和最大荷电状态;
Figure BDA0001803679040000087
Figure BDA0001803679040000088
分别表示电动汽车v的入网和离网时刻;
Figure BDA0001803679040000089
Figure BDA00018036790400000810
分别表示电动汽车v在入网时刻和离网时刻对应的电池荷电状态;Sv,a表示电动汽车v在入网时刻的电池荷电状态;Sv,d表示电动汽车v在离网时刻的目标电池荷电状态;
Figure BDA00018036790400000811
Figure BDA00018036790400000812
分别表示电动汽车v的充放电效率;Emax表示电动汽车电池容量;Δt表示时段长度。
作为另一个重要需求响应资源,水蓄冷空调系统建模如下:
(1)建筑物热动态平衡方程:
Figure BDA00018036790400000813
式中:θout,t和θin,t分别表示在时刻t的户外和室内温度;B是建筑物的热传导系数;Ca表示空气的比热容系数;V表示室内容积;ρa表示空气密度;Qt表示在时刻t建筑物的瞬时得热量,其主要和太阳辐射、建筑物内装置散热等有关;CD,t表示在时刻t建筑物内制冷装置提供的冷量。
(2)热舒适度建模:
采用PMV方法衡量建筑物内的热舒适度指标。人体可接受的PMV值在-0.5至0.5之间。PMV值和温度之间的关系为:
Figure BDA00018036790400000814
式中:IPMV表示室内的PMV值;θ表示温度,℃。
(3)水蓄冷空调系统建模:
单位时间内水蓄冷空调系统提供的冷量为:
CD,i,t=Cp,i,t-Cs,i,t+Cr,i,t
式中:Cp,i,t,Cs,i,t和Cr,i,t分别表示在时刻t水蓄冷空调i的制冷量、蓄冷量和放冷量。
水蓄冷空调需满足的运行约束为:
0≤Is,i,t+Ir,i,t≤1,
Figure BDA0001803679040000091
Figure BDA0001803679040000092
Figure BDA0001803679040000093
Figure BDA0001803679040000094
Figure BDA0001803679040000095
式中:Is,i,t和Ir,i,t分别表示在时刻t水蓄冷空调i的蓄冷和放冷运行状态变量;
Figure BDA0001803679040000096
Figure BDA0001803679040000097
分别表示在单位时间内水蓄冷空调的最大制冷量、蓄冷量和放冷量;
Figure BDA0001803679040000098
表示在时刻t水蓄冷空调i储存的冷量大小;Sc,max表示水蓄冷空调的蓄冷水箱容量;ηs和ηr分别表示蓄冷效率和放冷效率。
水蓄冷空调i在时刻t消耗的电功率为:
Figure BDA0001803679040000099
式中:μp表示能耗比,用来描述水蓄冷空调的制冷机的制冷效果和消耗的电功率之间的关系;μs和μr分别表示蓄冷量和放冷量与消耗电能的转换系数。
另外,优化调度模型中,下列燃气轮机运行约束和电力系统能量平衡约束仍需满足:
Figure BDA00018036790400000910
Figure BDA00018036790400000911
Figure BDA00018036790400000912
式中:
Figure BDA00018036790400000913
Figure BDA00018036790400000914
分别表示燃气轮机g的最小和最大输出有功功率;
Figure BDA00018036790400000915
Figure BDA00018036790400000916
分别表示燃气轮机g的爬坡和滑坡功率;Δt表示时段长度;Nq表示虚拟电厂内间歇性可再生能源发电机组(注:主要指风电机组和光伏发电机组)数量;Pk,t表示可再生能源发电机组k在时刻t的有功出力;Nc表示虚拟电厂内水蓄冷空调系统的数量;Pt R表示虚拟电厂在时刻t的刚性负荷需求。
步骤5、求解上述优化调度模型,得到虚拟电厂的日前优化调度方案,如:次日各个时刻燃气轮机的出力大小、各电动汽车的充放电功率、其他需求响应资源的用电情况等。
应用例
为了进一步理解本发明,以下以一个简单虚拟电厂优化调度问题为例,来解释本发明的实际应用。
假设虚拟电厂中含有2台燃气轮机、1台风电机组、1台光伏发电机组和50个智能用户。假设每一个智能用户都拥有1辆电动汽车和1台水蓄冷空调。电动汽车的初始电池荷电状态取在[0.15,0.35]内均匀分布的随机数。电动汽车的离网时刻和入网时刻分别服从以下两式所描述的概率密度函数。
Figure BDA0001803679040000101
Figure BDA0001803679040000102
式中:μd=8.92,σd=3.24,μa=17.47,σa=3.41。
水蓄冷空调系统的参数如表1所示;电力系统的分时电价如表2所示。
表1水蓄冷空调系统参数
Figure BDA0001803679040000103
表2电力系统分时电价
时段 电价($/kWh)
7:00-8:15 0.103
8:30-10:15 0.164
10:30-11:30 0.174
11:45-17:45 0.103
18:00-18:45 0.164
19:00-20:45 0.174
21:00-22:45 0.164
23:00-6:45 0.041
根据本发明所提出的方法,求解虚拟电厂优化调度模型。首先,得到一个典型水蓄冷空调系统的电功率和蓄冷量随时间变化情况如图2所示。可以看出,水蓄冷空调系统在需求响应机制的引导下,主要在负荷低谷时段消耗大量电功率,可以优化负荷曲线轮廓。
电动汽车充放电负荷随时间变化情况如图3所示。同水蓄冷空调系统类似,电动汽车负荷通过响应分时电价,可以起到削峰填谷的作用。
改变碳排放价格,虚拟电厂的碳排放量以及通过碳交易获得的收益变化如表3所示。可以看出,通过设计适当的碳交易机制以及合理安排碳排放价格,可以减少虚拟电厂的碳排放量,充分发挥虚拟电厂在碳减排方面的突出效益。
表3三种碳排放价格下虚拟电厂碳排放量及碳交易收益
Figure BDA0001803679040000111
通过以上分析,可见本发明提出的考虑需求响应和碳交易的虚拟电厂优化调度方法,不仅可以充分利用虚拟电厂内的需求响应资源以达到增加经济收益、优化电力系统负荷曲线的目的,还可以使虚拟电厂通过参与适当设计的碳交易机制以减少电力系统的碳排放量。这种优化调度方法,对电力系统的经济性、安全性和环保性均有一定的提升作用。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (1)

1.一种考虑需求响应和碳交易的虚拟电厂优化调度方法,其特征在于,包括步骤:
1)统计虚拟电厂内分布式电源的数量,以及电动汽车、水蓄冷空调系统参与需求响应的可控负荷数量,所述的分布式电源包括燃气轮机、风电机组和光伏发电机组;
2)对次日虚拟电厂内风电、光伏机组的出力和电力系统内刚性负荷进行预测;
3)虚拟电厂的能量管理系统接受电动汽车车主的日前申报信息,包括各电动汽车的入网及离网时刻、预估的初始电池荷电状态和期望的离网荷电状态;
4)虚拟电厂的能量管理系统根据日前预测信息、电动汽车申报信息以及电力市场和碳交易市场的价格信息,建立虚拟电厂同时参与电力市场和碳交易市场的优化调度模型;
5)求解虚拟电厂的优化调度模型,得到虚拟电厂的日前优化调度方案;
虚拟电厂优化调度模型的目标函数为虚拟电厂所获得的经济收益最大,其数学表达式如下所示:
Figure FDA0003143137000000011
Figure FDA0003143137000000012
Figure FDA0003143137000000013
Figure FDA0003143137000000014
Figure FDA0003143137000000015
Figure FDA0003143137000000016
式中:T表示一天内的时段数;
Figure FDA0003143137000000017
表示在时刻t的分时电价;
Figure FDA0003143137000000018
表示在时刻t虚拟电厂与配电系统的交互电功率;
Figure FDA0003143137000000019
Figure FDA00031431370000000110
分别表示在时刻t的燃气轮机运行成本、碳交易成本和电动汽车充放电成本;NG表示虚拟电厂内燃气轮机的数量;kg,0表示燃气轮机g的固定运行成本;ug,t表示在时刻t燃气轮机g的开关状态变量;nj表示燃气轮机运行成本的线性化段数;kg,j表示第j段线性化运行成本函数的斜率;Pg,j,t表示在时刻t燃气轮机g的第j段出力大小;Pg,t表示在时刻t燃气轮机g的总出力大小;
Figure FDA00031431370000000111
Figure FDA00031431370000000112
分别表示在时刻t燃气轮机g的启动和停机成本;
Figure FDA00031431370000000113
Figure FDA00031431370000000114
表示燃气轮机g启动和停机一次所消耗的成本;Nv表示虚拟电厂内电动汽车的数量;
Figure FDA00031431370000000115
Figure FDA00031431370000000116
分别表示电动汽车v在时刻t的充放电功率;
Figure FDA00031431370000000117
Figure FDA00031431370000000118
分别表示在时刻t的充放电电价;
虚拟电厂参与的碳交易机制表示为:
Figure FDA0003143137000000021
Figure FDA0003143137000000022
ΔMt=MC,t-MD,t
Figure FDA0003143137000000023
式中:N表示虚拟电厂中所有发电机组的数量;ε表示虚拟电厂的单位电量碳排放分配系数;Pn,t表示在时刻t发电机组n的有功出力;KC表示单位碳排放量交易价格;MC,t表示在时刻t虚拟电厂的实际碳排放量;MD,t表示在时刻t虚拟电厂分配到的碳排放配额;Qg表示燃气轮机g的单位电量碳排放强度;
如果在时刻t虚拟电厂的实际碳排放量低于分配得到的碳排放配额,虚拟电厂通过将多余的碳排放配额出售至碳交易市场以获取相应的利润;否则,虚拟电厂必须在碳交易市场中购买相应的碳排放权;
作为需求响应资源的电动汽车充放电负荷需满足的约束条件如下:
Figure FDA0003143137000000024
Figure FDA0003143137000000025
Figure FDA0003143137000000026
Sv,min≤Sv,t≤Sv,max
Figure FDA0003143137000000027
Figure FDA0003143137000000028
Figure FDA0003143137000000029
式中:
Figure FDA00031431370000000210
Figure FDA00031431370000000211
分别表示电动汽车v在时刻t的充电和放电状态变量;
Figure FDA00031431370000000212
Figure FDA00031431370000000213
分别表示电动汽车v的额定充放电功率;
Figure FDA00031431370000000214
表示电动汽车v在时刻t的调度状态;Sv,t表示电动汽车v在时刻t的电池荷电状态;Sv,min和Sv,max分别表示电动汽车v的最小和最大荷电状态;
Figure FDA00031431370000000215
Figure FDA00031431370000000216
分别表示电动汽车v的入网和离网时刻;
Figure FDA00031431370000000217
Figure FDA00031431370000000218
分别表示电动汽车v在入网时刻和离网时刻对应的电池荷电状态;Sv,a表示电动汽车v在入网时刻的电池荷电状态;Sv,d表示电动汽车v在离网时刻的目标电池荷电状态;
Figure FDA0003143137000000031
Figure FDA0003143137000000032
分别表示电动汽车v的充放电效率;Emax表示电动汽车电池容量;Δt表示时段长度;
作为需求响应资源的水蓄冷空调负荷建模如下:
1)建筑物热动态平衡方程
Figure FDA0003143137000000033
式中:θout,t和θin,t分别表示在时刻t的户外和室内温度;B是建筑物的热传导系数;Ca表示空气的比热容系数;V表示室内容积;ρa表示空气密度;Qt表示在时刻t建筑物的瞬时得热量;CD,t表示在时刻t建筑物内制冷装置提供的冷量;Δt表示时段长度;
2)热舒适度建模
采用PMV方法衡量建筑物内的热舒适度指标,PMV值和温度之间的关系为:
Figure FDA0003143137000000034
式中:IPMV表示室内的PMV值;θ表示温度,℃;
3)水蓄冷空调系统建模
单位时间内水蓄冷空调系统提供的冷量CD,i,t为:
CD,i,t=Cp,i,t-Cs,i,t+Cr,i,t
式中:Cp,i,t,Cs,i,t和Cr,i,t分别表示在时刻t水蓄冷空调i的制冷量、蓄冷量和放冷量;
水蓄冷空调需满足的运行约束为:
0≤Is,i,t+Ir,i,t≤1,
Figure FDA0003143137000000035
Figure FDA0003143137000000036
Figure FDA0003143137000000037
Figure FDA0003143137000000038
Figure FDA0003143137000000039
式中:Is,i,t和Ir,i,t分别表示在时刻t水蓄冷空调i的蓄冷和放冷运行状态变量;
Figure FDA00031431370000000310
Figure FDA00031431370000000311
分别表示在单位时间内水蓄冷空调的最大制冷量、蓄冷量和放冷量;
Figure FDA00031431370000000312
表示在时刻t水蓄冷空调i储存的冷量大小;Sc,max表示水蓄冷空调的蓄冷水箱容量;ηs和ηr分别表示蓄冷效率和放冷效率;Δt表示时段长度;
水蓄冷空调i在时刻t消耗的电功率Pi,t为:
Figure FDA0003143137000000041
式中:μp表示能耗比,用来描述水蓄冷空调的制冷机的制冷效果和消耗的电功率之间的关系;μs和μr分别表示蓄冷量和放冷量与消耗电能的转换系数;
需满足的燃气轮机运行约束和电力系统能量平衡约束如下:
Figure FDA0003143137000000042
Figure FDA0003143137000000043
Figure FDA0003143137000000044
式中:
Figure FDA0003143137000000045
Figure FDA0003143137000000046
分别表示燃气轮机g的最小和最大输出有功功率;
Figure FDA0003143137000000047
Figure FDA0003143137000000048
分别表示燃气轮机g的爬坡和滑坡功率;Δt表示时段长度;Nq表示虚拟电厂内间歇性可再生能源发电机组的数量;Pk,t表示可再生能源发电机组k在时刻t的有功出力;Nc表示虚拟电厂内水蓄冷空调系统的数量;
Figure FDA0003143137000000049
表示虚拟电厂在时刻t的刚性负荷需求。
CN201811088161.5A 2018-09-18 2018-09-18 一种考虑需求响应和碳交易的虚拟电厂优化调度方法 Active CN109523052B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811088161.5A CN109523052B (zh) 2018-09-18 2018-09-18 一种考虑需求响应和碳交易的虚拟电厂优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811088161.5A CN109523052B (zh) 2018-09-18 2018-09-18 一种考虑需求响应和碳交易的虚拟电厂优化调度方法

Publications (2)

Publication Number Publication Date
CN109523052A CN109523052A (zh) 2019-03-26
CN109523052B true CN109523052B (zh) 2021-09-10

Family

ID=65771535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811088161.5A Active CN109523052B (zh) 2018-09-18 2018-09-18 一种考虑需求响应和碳交易的虚拟电厂优化调度方法

Country Status (1)

Country Link
CN (1) CN109523052B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110165715B (zh) * 2019-05-31 2022-07-26 电子科技大学 一种将电动汽车储能式充电站接入虚拟电厂的方法
CN111179110B (zh) * 2019-12-06 2023-09-05 清华-伯克利深圳学院筹备办公室 虚拟电厂可变阶聚合等效鲁棒动态模型建模方法及装置
CN111784055B (zh) * 2020-07-02 2024-09-20 长沙理工大学 一种考虑电力需求侧响应特性的虚拟电厂优化调度方法
CN112348283B (zh) * 2020-11-26 2022-06-17 国网天津市电力公司电力科学研究院 蓄热式电采暖虚拟电厂日前可调度潜力评估方法及装置
CN113327124B (zh) * 2021-04-13 2024-06-28 东南大学 一种多能源-碳排放联合市场上的低碳p2p用能交易方法
CN113077111A (zh) * 2021-04-26 2021-07-06 上海电机学院 基于插电混动汽车v2g技术的虚拟电厂优化调度方法
CN113792953B (zh) * 2021-07-06 2024-03-19 国电南瑞南京控制系统有限公司 一种虚拟电厂优化调度方法及系统
CN114066046B (zh) * 2021-11-12 2022-11-11 国网江苏省电力有限公司镇江供电分公司 一种面向深度调峰的虚拟电厂内光储荷优化调度方法
CN114421540B (zh) * 2022-03-02 2024-09-20 国网浙江省电力有限公司经济技术研究院 一种基于虚拟电厂的分布式抽水蓄能调度方法
CN114971154A (zh) * 2022-03-31 2022-08-30 华南理工大学 一种包含碳交易机制的可再生能源消纳方法
CN114912943B (zh) * 2022-04-14 2024-07-19 国网江苏省电力有限公司盐城供电分公司 一种考虑碳排放交易的虚拟电厂中长期联合调度方法
CN114844100A (zh) * 2022-04-26 2022-08-02 东南大学 一种虚拟电厂参与配电网多目标优化调度模型与方法
CN115081715B (zh) * 2022-06-24 2024-07-16 东北电力大学 考虑碳交易与绿证交易的虚拟电厂优化调度方法
CN115907981B (zh) * 2022-11-08 2023-09-19 国网冀北张家口风光储输新能源有限公司 一种虚拟电厂低碳联合交易优化方法
DE202022106297U1 (de) 2022-11-09 2022-11-21 Kawther Ali Dhlan Aldhlan Intelligentes System für Energiebedarfsprognosen und nachhaltiges Energiemanagement durch maschinelles Lernen und künstliche Intelligenz
CN117318166B (zh) * 2023-08-09 2024-07-02 河海大学 考虑故障风险下基于合作博弈的多虚拟电厂低碳调度方法
CN117239810B (zh) * 2023-11-09 2024-03-26 南方电网数字电网研究院有限公司 虚拟电厂电能调度方案获取方法、装置和设备
CN117688520A (zh) * 2024-02-04 2024-03-12 国网安徽省电力有限公司经济技术研究院 一种基于电碳市场关联的电-碳价格传导数据分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105117805A (zh) * 2015-09-15 2015-12-02 武汉大学 一种利用电动汽车的虚拟发电厂优化调度方法及系统
CN105741027A (zh) * 2016-01-27 2016-07-06 河海大学 一种含电动汽车的虚拟电厂优化调度方法
CN105761109A (zh) * 2016-02-19 2016-07-13 刘隽琦 虚拟电厂能源管理和电力交易的智能管理系统及其优化运行方法
CN106096790A (zh) * 2016-06-22 2016-11-09 东南大学 基于变频空调虚拟机组建模的虚拟电厂日前和实时市场优化调度方法
CN107482688A (zh) * 2017-07-26 2017-12-15 长沙理工大学 一种碳捕集虚拟电厂调度优化方法
CN108510212A (zh) * 2018-04-17 2018-09-07 香港中文大学(深圳) 一种交互式能源系统的分布式能源调度方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105117805A (zh) * 2015-09-15 2015-12-02 武汉大学 一种利用电动汽车的虚拟发电厂优化调度方法及系统
CN105741027A (zh) * 2016-01-27 2016-07-06 河海大学 一种含电动汽车的虚拟电厂优化调度方法
CN105761109A (zh) * 2016-02-19 2016-07-13 刘隽琦 虚拟电厂能源管理和电力交易的智能管理系统及其优化运行方法
CN106096790A (zh) * 2016-06-22 2016-11-09 东南大学 基于变频空调虚拟机组建模的虚拟电厂日前和实时市场优化调度方法
CN107482688A (zh) * 2017-07-26 2017-12-15 长沙理工大学 一种碳捕集虚拟电厂调度优化方法
CN108510212A (zh) * 2018-04-17 2018-09-07 香港中文大学(深圳) 一种交互式能源系统的分布式能源调度方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"计及风险与源–网–荷双层协调的电力实时平衡调度";朱建全 等;《中国电机工程学报》;20150705;全文 *

Also Published As

Publication number Publication date
CN109523052A (zh) 2019-03-26

Similar Documents

Publication Publication Date Title
CN109523052B (zh) 一种考虑需求响应和碳交易的虚拟电厂优化调度方法
Liu et al. Two-phase collaborative optimization and operation strategy for a new distributed energy system that combines multi-energy storage for a nearly zero energy community
CN108154309B (zh) 计及冷热电多负荷动态响应的能源互联网经济调度方法
CN113962828B (zh) 考虑碳消耗的综合能源系统协调调度方法
CN110807588B (zh) 一种多能源耦合综合能源系统的优化调度方法
CN111339689B (zh) 建筑综合能源调度方法、系统、存储介质及计算机设备
CN111640044A (zh) 一种考虑虚拟储能的综合能源系统dg消纳策略研究方法
Yang et al. Research on the contribution of regional Energy Internet emission reduction considering time-of-use tariff
CN112186755A (zh) 一种区域综合能源系统柔性负荷储能化建模方法
CN112465236B (zh) 一种考虑综合满意度的社区综合能源系统调度方法
Yi et al. Low carbon economy scheduling of integrated energy system considering the mutual response of supply and demand
CN115912419A (zh) 一种基于负荷聚合商的储热-储电协同调度方法
CN115730747A (zh) 一种综合能源系统多主体利益分配方法及其应用
CN111967728B (zh) 计及用能舒适度时变性的商场楼宇调峰容量评估方法
CN113313329A (zh) 一种含综合能源系统的配电网优化调度方法
CN113128868A (zh) 一种区域综合能源系统调度优化方法及装置
CN113128799A (zh) 能源管控方法、装置、电子设备及计算机存储介质
CN110783917A (zh) 含有新能源消纳的多能源能量枢纽配置方法
CN115470564A (zh) 一种公共建筑能源系统协调控制方法及其控制组件
Shi et al. Economic operation of industrial microgrids with multiple kinds of flexible loads
Ma et al. Research on collaborative planning of distributed resources and energy storage in regional integrated energy system
CN107528352A (zh) 一种基于可再生能源高渗透率的配电网有功优化方法
CN113036751A (zh) 一种考虑虚拟储能的可再生能源微电网优化调度方法
Jiarui et al. Research on Demand Response Strategy of Electricity Market Based on Intelligent Power Consumption
Ling et al. Optimal dispatch of multi energy system using power-to-gas technology considering flexible load on user side

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