CN115660309A - 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法 - Google Patents

考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法 Download PDF

Info

Publication number
CN115660309A
CN115660309A CN202210991160.1A CN202210991160A CN115660309A CN 115660309 A CN115660309 A CN 115660309A CN 202210991160 A CN202210991160 A CN 202210991160A CN 115660309 A CN115660309 A CN 115660309A
Authority
CN
China
Prior art keywords
cost
carbon
power
vpp
load
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.)
Pending
Application number
CN202210991160.1A
Other languages
English (en)
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.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
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 China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN202210991160.1A priority Critical patent/CN115660309A/zh
Publication of CN115660309A publication Critical patent/CN115660309A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • Y—GENERAL 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

考虑电‑碳‑绿证联合交易的虚拟电厂日前优化调度方法，包括：构建电‑碳‑绿色证书联合交易框架；构建联合交易机制下VPP日前优化调度模型；构建优化调度模型的约束条件；对电‑碳‑绿证联合交易下VPP调度方案和市场交易策略进行风险度量；以VPP总成本最小为目标函数进行求解。首先对VPP分别在电‑碳市场、电‑绿证市场、碳‑绿证市场的市场交易策略和调度方案进行分析；其次，通过设置四种不同场景，对其在可再生能源消纳效果和二氧化碳排放的削减效果进行对比分析；最后，对于风光的不确定性，采用信息间隙决策理论对不确定参数进行处理，得出了风险规避模型和风险偏好模型约束下VPP的最优调度方案及市场交易策略，能够为不同投资者提供决策参考。

Description

考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法
技术领域
本发明涉及园区虚拟电厂目前优化调度技术领域，具体涉及一种考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法。
背景技术
当前，聚合大量可再生能源的虚拟电厂(virtual power plant，VPP)成为配网侧节能减排的重要载体，并可以通过配置诸如碳捕集系统、P2G等节能减排装置进一步降低C02的排放量。然而，目前有关促进VPP节能减排的市场交易机制不够完善，相关的补贴政策往往过于单一，无法充分发挥出VPP聚合资源的灵活性，导致VPP绿色化、低碳化发展水平较低。因此，亟需开展综合考虑电力市场交易、绿色证书交易以及碳排放权交易的VPP运行优化模型研究。
目前，国内外学者针对VPP在不同市场交易机制下的优化调度模型开展了一系列的研究，并得出了不同交易机制下VPP的最优调度策略。然而，研究中所涉及的碳交易机制大多为基于基准线法的碳排放权交易方法，碳价为单一固定的数值，无法体现对VPP节能减排的激励作用。此外，研究中碳捕集电厂的能耗来源大多为燃气轮机或火电机组，而VPP中的风光等可再生能源则全部用来上网，无法体现VPP在电-碳市场的选择性和交互性，且受限于联络线传输功率，其可再生能源的消纳水平也普遍不高。
在此背景下，可再生能源配额制以及绿色证书交易机制为提高其消纳率提供了新的途径。考虑到碳交易与绿色证书交易的互补特性，在减少碳排放量同时，尽可能的提高清洁能源的消纳比例。绿色证书交易与风光等可再生能源的出力密切相关，然而目前的研究并未考虑由风光出力的不确定性导致的绿色证书交易风险，且未涉及电-碳-绿色证书三者联合交易，因此需要进一步在风光不确定性下对VPP参与多市场的交易风险进行研究。
发明内容
为了解决上述技术问题，本发明提出了一种考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，该方法首先建立了电-碳-绿色证书联合交易的市场框架。在此基础上，搭建了VPP内部风-光-碳捕集系统协调调度框架，并进一步建立了联合交易机制下VPP日前调度模型；与此同时，为了准确衡量VPP在不同市场所面临的交易风险，引入信息间隙决策理论对VPP内部风电、光伏等可再生能源出力的不确定性进行处理，并得出了鲁棒模型和机会模型约束下VPP的最优调度方案以及市场交易策略，能够为不同投资者提供决策参考。
本发明采取的技术方案为：
考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，包括以下步骤：
步骤1：构建电-碳-绿色证书联合交易框架；
步骤2：构建联合交易框架下VPP日前优化调度模型；
步骤3：设置步骤2中优化调度模型的约束条件；
步骤4：对电-碳-绿证联合交易下VPP调度方案和市场交易策略进行风险度量；
步骤5：在MATLAB环境下，以VPP调度总成本最小为目标函数，通过YALMIP工具箱调用CPLEX求解器进行求解；
步骤6：首先，对VPP分别在电-碳市场、电-绿证市场、碳-绿证市场的市场交易策略和调度方案进行分析；其次，由于VPP聚合风-光-火碳捕集系统，通过设置四种不同场景，对其在可再生能源消纳效果和二氧化碳排放的削减效果进行对比分析；最后，对于风光的不确定性，采用信息间隙决策理论对不确定参数进行处理，得出了风险规避模型和风险偏好模型约束下VPP的最优调度方案以及市场交易策略。
所述步骤1中，
考虑VPP能够同时参与电-碳-绿色证书市场交易，构建了电-碳-绿色证书联合交易框架。VPP参与多个市场的交易流程为：首先，VPP向政府监管部门提出申请，监管部门审核后为VPP制定可再生能源消纳配额与碳排放配额；其次，VPP根据相关配额，在满足要求前提下在CET平台和GCT平台进行交易；最后，不同市场交易量确定后，按照相关结算规则进行结算，获取收益；
电-碳-绿色证书市场之间的耦合关系主要有以下两个方面：
首先，由于VPP内部聚合风-光-火碳捕集系统，当电力市场价格较高时，则风-光-火大部分电能用于上网以从电力市场中获取较高收益，此时用于处理CO2排放的能耗较低，因此VPP需要从碳排放市场购买较多的碳排放配额，即电-碳市场存在一定的耦合关系；其次，可再生能源的消纳比重与VPP内部的等效负荷相关，等效负荷的变化会影响相应购售电量的变化，因此，绿色证书市场交易量也与电力市场密切相关；
综上所述，在电-碳-绿色证书联合交易框架下，VPP将根据其在不同市场获取的潜在效益，合理制定内部聚合单元的运行计划，从而获取更高的经济效益。
所述步骤2中，
VPP整体参与市场包括：碳交易市场、电力市场、绿色证书交易市场；
VPP内部聚合资源包括以下部分：
电源测包括：风电、光伏、传统火电基础上改装而成的碳捕集电厂；
负荷侧包括：基础负荷、可转移负荷、可中断负荷、中央空调负荷；
储能设备包括：蓄电池；
在电-碳-绿色证书联合交易框架下，以VPP运行成本最小为目标函数，构建联合交易机制下VPP日前优化调度模型，如式(1)所示，具体包含以下六个部分，分别为碳捕集电厂燃料成本、储能电池损耗成本、可控负荷成本、碳排放权交易成本、绿色证书交易成本和弃风、弃光惩罚成本：
Figure BDA0003751864390000031
式(1)中，F为VPP总成本，
Figure BDA0003751864390000032
为t时刻VPP的碳捕集电厂燃料成本、
Figure BDA0003751864390000033
为储能电池损耗成本、
Figure BDA0003751864390000034
为可控负荷成本、
Figure BDA0003751864390000035
为碳排放权交易成本、
Figure BDA0003751864390000036
为绿色证书交易成本、
Figure BDA0003751864390000037
为弃风、弃光惩罚费用。T表示调度周期，取24小时；t表示时间。
(i)碳捕集电厂燃料成本
Figure BDA0003751864390000038
为简化计算，碳捕集电厂成本不考虑阀点效应，故总燃料费用计算式为：
Figure BDA0003751864390000039
式(2)中，
Figure BDA00037518643900000310
为t时刻VPP的碳捕集电厂燃料成本，a、b、c分别为火电机组的燃料成本系数，Pt G为火电机组t时刻的总出力。
(ii)储能系统的电池损耗费用
Figure BDA00037518643900000311
计算如下：
Figure BDA00037518643900000312
式(3)中，
Figure BDA00037518643900000313
为储能电池损耗成本，cess为电池损耗成本系数，Pt esc、Pt esd分别表示储能的充、放电功率。
(iii)可控负荷的调控费用
Figure BDA00037518643900000314
可控负荷的调控费用主要包括可转移负荷费用
Figure BDA00037518643900000315
可中断负荷的补偿费用
Figure BDA00037518643900000316
以及空调负荷调控引起的用户不舒适度补偿
Figure BDA00037518643900000317
Figure BDA00037518643900000318
Figure BDA0003751864390000041
Figure BDA0003751864390000042
Figure BDA0003751864390000043
式(4)-式(7)中，
Figure BDA0003751864390000044
可控负荷成本，
Figure BDA0003751864390000045
分别为t时刻中断负荷、转移负荷以及空调负荷的调控费用；cm、cs、cres分别为中断负荷、转移负荷及空调负荷的成本系数；
Figure BDA0003751864390000046
表示第m级中断负荷的响应量；Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷，Tt in表示室内温度，Tt ref表示参考室温。
(iv)碳交易成本
Figure BDA0003751864390000047
Figure BDA0003751864390000048
式(8)中，
Figure BDA0003751864390000049
为阶梯式碳交易成本；λ为碳交易基础价格，
Figure BDA00037518643900000410
为VPP实际的碳排放量，l为碳排放区间长度，α为价格增长率。
(v)绿色证书交易成本
Figure BDA00037518643900000411
Figure BDA00037518643900000412
式(9)中，
Figure BDA00037518643900000413
为绿色证书交易成本，Pt gc为绿色证书交易价格，
Figure BDA00037518643900000414
为VPP在t时刻购买或出售的绿色证书量，若为正值表示VPP此刻购买绿色证书，若为负值表示此刻出售绿色证书。
(vi)弃风、弃光惩罚成本
Figure BDA00037518643900000415
Figure BDA00037518643900000416
式(10)中，
Figure BDA00037518643900000417
为弃风、弃光惩罚成本，cpv、cpw分别为弃光、弃风惩罚费用，Pt ppv、Pt ppw分别为光伏、风电出力的预测值，Pt pv、Pt pw为实际消纳的光伏、风电功率。
所述步骤3中，VPP日前优化调度模型的约束条件主要为聚合单元的运行约束、功率平衡约束以及市场交易约束等，具体包括以下部分：
(a)风-光-火碳捕集系统运行约束：
Pt GN＝Pt G-Pt GC (11)；
Pt GC＝Pt B+Pt OP (12)；
式(11)、式(12)中，Pt GN为碳捕集电厂的实际出力；Pt G为火电总出力；Pt GC为碳捕集的总能耗，包括固定能耗Pt B及运行能耗Pt OP；
将风光机组与碳捕集机组联合运行，风光的部分功率可以作为碳捕集系统的运行能耗使用，具体如下：
Figure BDA0003751864390000055
Figure BDA0003751864390000051
式(13)、式(14)中，Pt OP为运行能耗；Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率；wc为处理单位CO2需要的功耗；
Figure BDA0003751864390000052
为再生塔处理的CO2量；
Figure BDA0003751864390000053
为实际火电碳排放量；eg为单位火电机组出力产生的CO2量；Pt G为火电总出力。
风光机组除了供给碳捕集系统，余电依然可以由VPP向外部电网出售，表达式如下：
Pt pv＝Pt vc+Pt vw (15)；
Pt pw＝Pt wc+Pt ww (16)；
式(15)、式(16)中，Pt pv、Pt pw为实际消纳的光伏、风电功率；Pt vw、Pt ww分别为光伏以及风电的上网功率；Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率。
风、光、火电机组需满足最大出力约束：
0≤Pt pw≤Pt ppw (17)；
0≤Pt pv≤Pt ppv (18)；
Pgmin≤Pt G≤Pgmax (19)；
式(17)、式(18)、式(19)中，Pgmin、Pgmax分别为火电机组出力的技术最小值和最大值；Pt G为火电总出力；Pt pv、Pt pw为实际消纳的光伏、风电功率；Pt ppv、Pt ppw分别为光伏、风电出力的预测值。此外，火电机组还须满足爬坡率约束：
Figure BDA0003751864390000054
式(20)中，Pt G为t时刻火电出力，
Figure BDA0003751864390000061
为t-1时刻火电出力，Pramp为火电机组的爬坡功率。
(b)储能装置运行约束：
Figure BDA0003751864390000062
式(21)中，
Figure BDA0003751864390000063
分别为t时刻、t-1储能蓄电量值；ηc、ηd为储能系统的充放电效率；Pt esc、Pt esd分别表示储能的充、放电功率。
(c)可控负荷运行约束：
Figure BDA0003751864390000064
Figure BDA0003751864390000065
Figure BDA0003751864390000066
式(22)、式(23)、式(24)中，
Figure BDA0003751864390000067
为第m级削减负荷的最大削减量，
Figure BDA0003751864390000068
为第m级可削减负荷，Pt iL为t时刻的总削减负荷，
Figure BDA0003751864390000069
为t-1时刻总削减负荷；M表示削减负荷的级别。Pcmax为最大连续削减负荷量；式(24)表示避免连续负荷削减导致用户舒适度下降。
Figure BDA00037518643900000610
Figure BDA00037518643900000611
Figure BDA00037518643900000612
Figure BDA00037518643900000613
式(25)、式(26)、式(27)、式(28)中，Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；Pt f max为t时段最大转移负荷；
Figure BDA00037518643900000614
分别表示t时刻是否有负荷转移到别处或从别处转移过来，为布尔变量。
式(27)表示负荷的转移在某个确定的时段式单向进行的，式(28)表明负荷转移前后总负荷功率应保持不变。
Figure BDA00037518643900000615
Tin min≤Tt in≤Tin max (30)
式(29)、式(30)中，γ、αt、β为与建筑物蓄热特性相关的参数；σ为空调的制冷能效比；Pt cold为中央空调的制冷功率；Tt in为t时刻温度，
Figure BDA0003751864390000071
为t-1时刻温度；Tin min、Tin max分别为满足用户舒适度的温度区间上、下限。
(d)系统功率平衡约束：
Pt cold+Pt loads+Pt esc+Pt mgs＝Pt vw+Pt ww+Pt GN+Pt mgb+Pt esd (31)
Pt loads＝Pt load-Pt sfl+Pt sfq-Pt iL (32)
式(31)、式(32)中，Pt cold为中央空调的制冷功率；Pt esc、Pt esd分别表示储能的充、放电功率；Pt loads为t时刻的等效负荷；Pt load为t时刻的基础负荷；Pt vw、Pt ww分别为光伏以及风电的上网功率；Pt mgs、Pt mgb分别为t时刻VPP向外部电网出售以及购入的电量；Pt GN为碳捕集电厂的实际出力；Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；Pt iL为t时刻的总削减负荷。
(e)联络线传输功率约束：
0≤Pt mgs≤Pmgmax (33)
0≤Pt mgb≤Pmgmax (34)
式(33)、式(34)中，Pt mgmax表示联络线最大传输功率，Pt mgs为售电量，Pt mgb为购电量。
(f)碳排放市场运行约束：
采用基准线法确定VPP的碳排放配额：
Et D＝ηbPt mgb+ηgPt G (35)
式(35)中，
Figure BDA0003751864390000072
为VPP核发的碳配额总量，ηb、ηg为从市场购电及火电机组的配额系数，Pt mgb为购电量，Pt G为火电总出力。
VPP的实际碳排放量计算如下：
Et P＝ap(Pt mgb)2+bpPt mgb+cP+Et N (36)
式(36)中，
Figure BDA0003751864390000073
为VPP的实际碳排放量，ap、bp、cP为实际排放系数，Pt mgb为购电量，
Figure BDA0003751864390000074
为碳捕集电厂净排放。
碳排放市场交易需求计算如下：
Figure BDA0003751864390000075
式(37)中，
Figure BDA0003751864390000081
为VPP碳排放市场交易量，
Figure BDA0003751864390000082
为VPP的实际碳排放量，
Figure BDA0003751864390000083
为VPP核发的碳配额总量。
(g)VPP消纳权责及绿证交易约束：
VPP内可再生能源的消纳量与绿色证书交易量总和应不小于消纳权责规定的电力消纳总量：
Figure BDA0003751864390000084
Tt total＝Pt loads+Pt cold+Pt esc-Pt esd (39)
Figure BDA0003751864390000085
式(38)、式(39)、式(40)中，
Figure BDA0003751864390000086
表示可再生能源消纳总量；w为消纳权重；
Figure BDA0003751864390000087
为总电量；
Figure BDA0003751864390000088
为VPP在t时刻购买或出售的绿色证书量；Pt loads为t时刻的等效负荷；Pt cold为中央空调的制冷功率；Pt esc、Pt esd分别表示储能的充、放电功率；Pt pv、Pt pw分别为实际消纳的光伏、风电功率。
所述步骤4中，由于风光出力具有较强的随机性与间歇性，因此VPP在日前制定的各类型市场申报策略和调度方案可能会与实际日内申报结果产生偏差，从而遭受高额的不平衡惩罚。采用信息间隙决策理论(Information Gap Decision Theory，IGDT)处理不确定性参数；
按照预设目标的好坏，不确定性参数对目标函数的影响分为积极和消极两个方面，相应的IGDT模型又可分为风险偏好型(机会)以及风险规避型(鲁棒)两种，其反应了不同投资者在面对风险时的价值取向。一种认为不确定性参数将会有利于目标函数实现，另一种则认为不确定性参数将会对目标函数产生负面影响，二者对应的数学模型分别如下：
Figure BDA0003751864390000089
Figure BDA0003751864390000091
式(41)、式(42)中，B0为预设目标值，通常由确定性的模型求解得出；Bc、Ba为不同风险偏向模型下的期望目标函数；βc、βa为偏差因子，表示预期成本与预设目标的偏差程度；H(X，d)＝0、G(X，d)≥0分别为优化模型的等式与不等式约束；α为不确定性参数的波动程度；
Figure BDA0003751864390000092
为不确定性参数的预测值；X为不确定参数实际值；d为决策变量，
Figure BDA0003751864390000093
为不确定参数波动范围函数，s.t.表示约束条件。
式(41)对应IGDT鲁棒模型，为保守者的决策模型，在最大化不确定性扰动下以期实现最低的预设目标；
式(42)对应IGDT机会模型，为激进者决策模型，在最小化不确定性扰动下以期实现最高的预设目标。
本发明一种考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，技术效果如下：
1)本发明能够让虚拟电厂在电-碳-绿色证书联合交易框架下，根据其在不同市场获取的潜在收益，合理制定内部聚合单元的运行计划，从而获取更高的经济效益；
2)风-光-火联合碳捕集电厂为VPP的调度策略提供了更高的灵活性和可靠性。通过碳捕集系统可以优先将无法消纳的可再生能源消纳，其余不足的能耗再由边际成本较高的火电机组供给，因此整体运行策略更加具有经济性；
3)针对风光出力具有较强的随机性与间歇性，导致VPP在目前制定的各类型市场申报策略可能会与实际日内申报结果产生偏差，从而遭受高额的不平衡惩罚。本发明引入信息间隙决策理论对不确定性参数进行处理，根据决策者偏好，得出了鲁棒模型和机会模型约束下VPP的最优调度方案及市场交易策略。
附图说明
图1是电-碳-绿色证书联合交易框架图。
图2为本发明所建立电-碳-绿色证书联合交易下VPP日前调度模型。
图3为风电、光伏、负荷日前预测曲线图。
图4为分时电价下碳市场交易图。
图5为碳捕集电厂C02处理于实际排放图。
图6为碳捕集电厂电源供给图。
图7为GCT交易量及等效负荷图。
图8为联合交易下VPP调度成本三维曲面图。
图9为扰动水平随调度成本变化曲线图。
图10为不同模型下VPP市场交易策略对比图。
具体实施方式
如图1所示，为电-碳-绿色证书联合交易框架图，电-碳-绿色证书市场之间的耦合关系主要有以下几个方面。首先，电-碳市场存在一定的耦合关系，由于VPP内部聚合风-光-火碳捕集系统，在电价峰时段，VPP会选择减少碳捕集系统能耗，在满足负荷需求情况下，更多的参与电市场以获利；在电价谷时段，此时不能在电市场获得更高收益，则会选择通过风、光、火更多为捕集供给能耗，减少二氧化碳的排放，参与碳市场以获得更高利润，风、光为碳捕集提供能耗的同时也意味着会优先对可能无法消纳的可再生能源进行了消纳，进而减少了弃风、弃光量；其次，可再生能源的消纳比重与VPP内部的等效负荷相关，等效负荷的变化会影响相应购售电量的变化，因此，绿色证书市场交易量也与电力市场密切相关；此外，碳市场也与绿证市场有着紧密的联系，随着火电出力的增加，二氧化碳排放量也会不断增加，在考虑碳价随碳排放量递增的阶梯型碳交易市场下，碳排放成本会越来越高，因此通过风-光-火碳捕集电厂协调运行，可以在促进可再生能源消纳进而达到可再生能源消纳权责的同时，减少火电机组的出力，进而减少碳排放成本。在电-碳-绿色证书联合交易框架基础上，搭建了VPP内部风-光-碳捕集系统协调调度框架，并进一步建立了联合交易机制下VPP日前调度模型。
由于风光出力具有较强的随机性与间歇性，因此VPP在日前制定的各类型市场的申报策略和调度方案可能会与实际日内申报结果产生偏差，从而遭受高额的不平衡惩罚。为了应对风光出力不确定性带来的潜藏风险，采用信息间隙决策理论(Information GapDecision Theory，IGDT)处理不确定性参数，并得出了风险规避模型和风险偏好模型约束下VPP的最优调度策略及市场交易策略，最后通过算例验证了所构模型和所用算法的有效性。具体包括以下步骤：
步骤1、构建电-碳-绿色证书联合交易框架：
鉴于不同市场之间良好的互动能力和互补潜力，考虑VPP可同时参与电-碳-绿色证书市场交易，构建了电-碳-绿色证书联合交易框架。VPP参与多个市场的交易流程为：首先，VPP向政府监管部门提出申请，监管部门审核后为VPP制定可再生能源消纳配额与碳排放配额；其次，VPP根据相关配额，在满足要求前提下在CET平台和GCT平台进行交易；最后，不同市场交易量确定后，按照相关结算规则进行结算，获取收益。
此外，电-碳-绿色证书市场之间的耦合关系主要有以下两个方面：首先，由于VPP内部聚合风-光-火碳捕集系统，当电力市场价格较高时，则风-光-火大部分电能用于上网以从电力市场中获取较高收益，此时用于处理CO2排放的能耗较低，因此VPP需要从碳排放市场购买较多的碳排放配额，即电-碳市场存在一定的耦合关系；其次，可再生能源的消纳比重与VPP内部的等效负荷相关，等效负荷的变化会影响相应购售电量的变化，因此，绿色证书市场交易量也与电力市场密切相关。综上所述，在电-碳-绿色证书联合交易框架下，VPP将根据其在不同市场获取的潜在效益，合理制定内部聚合单元的运行计划，从而获取更高的经济效益。
步骤2、构建联合交易框架下VPP日前优化调度模型：
VPP整体参与的市场交易包括：碳交易市场、电力市场、绿色证书交易市场；
VPP内部聚合资源包括以下部分：
电源测包括：风电、光伏、传统火电基础上改装而成的碳捕集电厂；
负荷侧包括：基础负荷、可转移负荷、可中断负荷、中央空调负荷；
储能设备包括：蓄电池；
在电-碳-绿色证书联合交易机制下，以VPP运行成本最小为目标函数，如式(1)所示，具体包含以下六个部分，分别为碳捕集电厂燃料成本、储能电池损耗成本、可控负荷成本、碳排放权交易成本、绿色证书交易成本和弃风、弃光惩罚成本。
Figure BDA0003751864390000111
式(1)中，F为VPP总成本；
Figure BDA0003751864390000112
为t时刻VPP的碳捕集电厂燃料成本；
Figure BDA0003751864390000113
为储能电池损耗成本；
Figure BDA0003751864390000114
为可控负荷成本；
Figure BDA0003751864390000115
为碳排放权交易成本；
Figure BDA0003751864390000116
为绿色证书交易成本；
Figure BDA0003751864390000117
为弃风、弃光惩罚费用。
(i)碳捕集电厂燃料成本
Figure BDA0003751864390000121
为简化计算，碳捕集电厂成本不考虑阀点效应，故总燃料费用计算式为：
Figure BDA0003751864390000122
式(2)中，
Figure BDA0003751864390000123
为t时刻VPP的碳捕集电厂燃料成本；
a、b、c分别为火电机组的燃料成本系数；
Pt G为火电机组t时刻的总出力。
(ii)储能系统的电池损耗费用
Figure BDA0003751864390000124
计算如下：
Figure BDA0003751864390000125
式(3)中，
Figure BDA0003751864390000126
为储能电池损耗成本；
cess为电池损耗成本系数；
Pt esc、Pt esd分别表示储能的充、放电功率。
(iii)可控负荷的调控费用
Figure BDA0003751864390000127
可控负荷的调控费用主要包括可转移负荷费用
Figure BDA0003751864390000128
可中断负荷的补偿费用
Figure BDA0003751864390000129
以及空调负荷调控引起的用户不舒适度补偿
Figure BDA00037518643900001210
Figure BDA00037518643900001211
Figure BDA00037518643900001212
Figure BDA00037518643900001213
Figure BDA00037518643900001214
式(4)-式(7)中，
Figure BDA00037518643900001215
可控负荷成本；
Figure BDA00037518643900001216
分别为t时刻中断负荷、转移负荷以及空调负荷的调控费用；
cm、cs、cres分别为中断负荷、转移负荷以及空调负荷的成本系数；
Figure BDA00037518643900001217
表示第m级中断负荷的响应量；
Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；
Tt in表示室内温度；
Tt ref表示参考室温。
(iv)碳交易成本
Figure BDA0003751864390000131
Figure BDA0003751864390000132
式(8)中，
Figure BDA0003751864390000133
为阶梯式碳交易成本；
λ为碳交易基础价格；
Et VPP为VPP实际的碳排放量；
l为碳排放区间长度；
α为价格增长率。
(v)绿色证书交易成本
Figure BDA0003751864390000134
Figure BDA0003751864390000135
式(9)中，
Figure BDA0003751864390000136
为绿色证书交易成本；
Pt gc为绿色证书交易价格；
Figure BDA0003751864390000137
为VPP在t时刻购买或出售的绿色证书量
(vi)弃风、弃光惩罚成本
Figure BDA0003751864390000138
Figure BDA0003751864390000139
式(10)中，
Figure BDA00037518643900001310
为弃风、弃光惩罚成本；
cpv、cpw分别为弃光、弃风惩罚费用；
Pt ppv、Pt ppw分别为光伏、风电出力的预测值；
Pt pv、Pt pw为实际消纳的光伏、风电功率。
步骤3、构建优化调度模型的约束条件：
VPP日前调度模型的约束条件主要为聚合单元的运行约束、功率平衡约束以及市场交易约束等，具体包括以下部分：
(a)风-光-火碳捕集系统运行约束：
Pt GN＝Pt G-Pt GC (11)
Pt GC＝Pt B+Pt OP (12)
式(11)-式(12)中，Pt GN为碳捕集电厂的实际出力；
Pt G为火电总出力；
Pt GC为碳捕集的总能耗；
固定能耗Pt B；
运行能耗Pt OP。
将风光机组与碳捕集机组联合运行，风光的部分功率可以作为碳捕集系统的运行能耗使用，具体如下：
Figure BDA0003751864390000141
Figure BDA0003751864390000142
式(13)-式(14)中，Pt OP为运行能耗，
Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率；
wc为处理单位CO2需要的功耗；
Figure BDA0003751864390000143
为再生塔处理的CO2量；
Figure BDA0003751864390000144
为实际火电碳排放量；
eg为单位火电机组出力产生的CO2量；
Pt G为火电总出力。
风光机组除了供给碳捕集系统，余电依然可以由VPP向外部电网出售，表达式如下：
Pt pv＝Pt vc+Pt vw (15)
Pt pw＝Pt wc+Pt ww (16)
式(15)-式(16)中，Pt pv、Pt pw为实际消纳的光伏、风电功率；
Pt vw、Pt ww分别为光伏以及风电的上网功率；
Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率。
风、光、火电机组需满足最大出力约束：
0≤Pt pw≤Pt ppw (17)
0≤Pt pv≤Pt ppv (18)
Pgmin≤Pt G≤pgmax (19)
式(17)-式(19)中，Pgmin、Pgmax分别为火电机组出力的技术最小值和最大值；
Pt G为火电总出力；
Pt pv、Pt pw为实际消纳的光伏、风电功率；
Pt ppv、Pt ppw分别为光伏、风电出力的预测值。
此外，火电机组还须满足爬坡率约束：
Figure BDA0003751864390000151
式(20)中，Pt G为t时刻火电总出力；
Figure BDA0003751864390000152
为(t-1)时刻；
Pramp为火电机组的爬坡功率。
(b)储能装置运行约束：
Figure BDA0003751864390000153
式(21)中，
Figure BDA0003751864390000154
分别为t时刻、t-1储能蓄电量值，
ηc、ηd为储能系统的充放电效率，
Pt esc、Pt esd分别表示储能的充、放电功率。
(c)可控负荷运行约束：
Figure BDA0003751864390000155
Figure BDA0003751864390000156
Figure BDA0003751864390000157
式(22)-式(24)中，
Figure BDA0003751864390000158
为第m级削减负荷的最大削减量；
Figure BDA0003751864390000159
为第m级可削减负荷；
Pt iL为t时刻的总削减负荷；
Figure BDA00037518643900001510
为(t-1)时刻总削减负荷；
Pcmax为最大连续削减负荷量。
Figure BDA0003751864390000161
Figure BDA0003751864390000162
Figure BDA0003751864390000163
Figure BDA0003751864390000164
式(25)-式(28)中，Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；
Pt f max为t时段最大转移负荷；
Figure BDA0003751864390000165
分别表示t时刻是否有负荷转移到别处或从别处转移过来。
Figure BDA0003751864390000166
Tin min≤Tt in≤Tin max (30)
式(29)-式(30)中，γ、αt、β为与建筑物蓄热特性相关的参数；
σ为空调的制冷能效比；
Pt cold为中央空调的制冷功率；
Tt in为t时刻温度；
Tin min、Tin max为满足用户舒适度的温度区间上下限。
(d)系统功率平衡约束：
Pt cold+Pt loads+Pt esc+Pt mgs＝Pt vw+Pt ww+Pt GN+Pt mgb+Pt esd (31)
Pt loads＝Pt load-Pt sfl+Pt sfq-Pt iL (32)
式(31)-式(32)中，Pt cold为中央空调的制冷功率；
Pt esc、Pt esd分别表示储能的充、放电功率；
Pt loads为t时刻的等效负荷；
Pt load为t时刻的基础负荷；
Pt vw、Pt ww分别为光伏以及风电的上网功率；
Pt mgs、Pt mgb分别为t时刻VPP向外部电网出售以及购入的电量；
Pt GN为碳捕集电厂的实际出力；
Pt iL为t时刻的总削减负荷；
Pt load为基础负荷；
Pt sfl、Pt sfq分别表示转移到别处的负荷和从别处转移来的负荷。
(e)联络线传输功率约束：
0≤Pt mgs≤Pmgmax (33)
0≤Pt mgb≤Pmgmax (34)
式(33)-式(34)中，Pt mgmax表示联络线最大传输功率；
Pt mgs为售电量；
Pt mgb为购电量。
(f)碳排放市场运行约束：
采用基准线法确定VPP的碳排放配额：
Figure BDA0003751864390000171
式(35)中，
Figure BDA0003751864390000172
为VPP核发的碳配额总量；
ηb、ηg为从市场购电及火电机组的配额系数；
Pt mgb为购电量；
Pt G为火电总出力。
VPP的实际碳排放量计算如下：
Figure BDA0003751864390000173
式中，
Figure BDA0003751864390000174
为VPP的实际碳排放量；
ap、bp、cP为实际排放系数；
Pt mgb为购电量；
Figure BDA0003751864390000175
为碳捕集电厂净排放；
碳排放市场交易需求计算如下：
Figure BDA0003751864390000176
式(37)中，
Figure BDA0003751864390000177
为VPP碳排放市场交易量；
Figure BDA0003751864390000178
为VPP的实际碳排放量；
Figure BDA0003751864390000181
为VPP核发的碳配额总量。
(g)VPP消纳权责及绿证交易约束：
VPP内可再生能源的消纳量与绿色证书交易量总和应不小于消纳权责规定的电力消纳总量：
Figure BDA0003751864390000182
Et total＝Pt loads+Pt cold+Pt esc-Pt esd (39)
Figure BDA0003751864390000183
式(38)-式(40)中，
Figure BDA0003751864390000184
表示可再生能源消纳总量；
w为消纳权重；
Figure BDA0003751864390000185
为总电量；
Figure BDA0003751864390000186
为VPP在t时刻购买或出售的绿色证书量；
Pt loads为t时刻的等效负荷；
Pt cold为中央空调的制冷功率；
Pt esc、Pt esd分别表示储能的充、放电功率；
Pt pv、Pt pw为实际消纳的光伏、风电功率。
步骤4、对电-碳-绿证联合交易下VPP调度方案和市场交易策略的风险度量：
由于风光出力具有较强的随机性与间歇性，因此VPP在日前制定的各类型市场申报策略和调度方案可能会与实际日内申报结果产生偏差，从而遭受高额的不平衡惩罚。采用信息间隙决策理论(Information Gap Decision Theory，IGDT)处理不确定性参数。
按照预设目标的好坏，不确定性参数对目标函数的影响分为积极和消极两个方面，相应的IGDT模型又可分为风险偏好型(机会)以及风险规避型(鲁棒)两种，其反应了不同投资者在面对风险时的价值取向。一种认为不确定性参数将会有利于目标函数实现，另一种则认为不确定性参数将会对目标函数产生负面影响，二者对应的数学模型分别如下：
Figure BDA0003751864390000191
Figure BDA0003751864390000192
式(41)-式(42)中，B0为预设目标值，通常由确定性的模型求解得出；
Bc、Ba为风险规避模型和风险偏好模型下的期望目标函数；
βc、βa为偏差因子，表示预期成本与预设目标的偏差程度；
H(X，d)＝0、G(X，d)≥0分别为优化模型的等式与不等式约束；
α为不确定性参数的波动程度；
Figure BDA0003751864390000193
为不确定性参数的预测值；
X为不确定参数实际值；
d为决策变量；
Figure BDA0003751864390000194
为不确定参数波动范围函数；
S.t.表示约束条件。
式(41)对应IGDT鲁棒模型，为保守者的决策模型，在最大化不确定性扰动下以期实现最低的预i设目标；式(42)对应IGDT机会模型，为激进者决策模型，在最小化不确定性扰动下以期实现最高的预设目标。
步骤5、在MATLAB环境中进行求解：
本发明在MATLAB环境下，以VPP总成本最小为目标，通过YALMIP工具箱调用CPLEX求解器进行求解。
步骤6、VPP参与多种市场交易策略和调度方案风险度量并进行算例分析：
本发明以图2虚拟电厂为例，验证所构模型以及所提算法的有效性。风、光以及负荷的日前预测值如图3所示，VPP与上级电网的交易电价如表1所示，储能参数方面，最大充放电功率均取为1MW，储能容量为4MW，初始蓄电量为0，充放电效率取为95％；可削减负荷的各级最大削减比例为15％、10％、8％，各级削减负荷的补偿费用为500/700/800元/MW，联络线的传输功率取为10MW，转移负荷的补偿费用为50元/MW，单时段最大转移负荷比例为20％，碳排放配额取火电为0.385，市场购电为0.798，基础碳价取为0.25元/kg，火电最大出力3.56MW，取最低消纳权重为12.5％，激励型消纳权责为13.8％。
(1)：首先对VPP参与电-碳市场所获经济效益进行分析，仿真结果如图4所示。
从图4中可以明显看出，碳排放权交易与分时电价之间具有显著的相关性，即说明电力市场与碳市场存在相关耦合特性。具体而言，当电价较低，VPP不能从电力市场获取高额收益时，VPP增加碳市场的交易量，以期从碳排放市场获取更高的收益。相反，当电力市场价格较高时，VPP从电力市场可以获取较高的收益，此时VPP的碳市场交易量较低，甚至需要大量从外部购入碳排放权。因此，通过电-碳市场联动，不仅提高了VPP市场交易策略的灵活性，也提高了VPP整体的经济效益。
表1 VPP电力市场交易价格
Figure BDA0003751864390000201
VPP之所以能够实现电-碳市场的友好互动，很大程度上依赖于自身聚合的风-光-碳捕集系统的协调调度，通过碳捕集系统控制VPP实际排放的CO2，从而能够实现对电价以及碳价的准确响应，VPP调度周期内再生塔处理的CO2量和碳捕集电厂实际排放的CO2量如图5所示，碳捕集电厂的能耗组成如图6所示。
从图5可以看出，在电价较低的时期，VPP通过增加碳捕集能耗，增加再生塔对CO2的吸收，整体碳捕集系统对外实现了零碳排放，从而可以对外出售更多的CO2排放权。
从图6中可以看出，风-光-火联合碳捕集系统为VPP的调度策略提供了更高的灵活性和可靠性。通过碳捕集系统可以优先将无法消纳的可再生能源消纳，其余不足的能耗再由边际成本较高的火电机组供给，因此整体运行策略更加具有经济性，且充足的能耗供给来源使得VPP可以精准响应外部电价和碳价。
为了进一步证明本发明所构模型的经济性和促进可再生能源的消纳，分别设置以下方案：
方案1：风-光-火碳捕集系统加阶梯式碳交易；
方案2：风-光-火碳捕集系统加传统碳交易；
方案3：仅火电碳捕集系统加阶梯式碳交易；
方案4：仅火电碳捕集系统加传统碳交易；
四种方案下可再生能源的消纳率以及整体VPP的目标函数如表2所示，对比表2可以看出，四种场景下VPP的最终调度结果有较为明显的区别。对比方案1和方案2可知，阶梯式碳交易与传统碳交易相比能够进一步引导VPP参与碳市场交易，降低VPP参与碳市场交易的成本，且阶梯式碳交易下可再生能源的消纳率略有提升。对比方案1以及方案3可知，风-光-火碳捕集系统能够充分利用风光等可再生能源，降低碳捕集的成本，降低VPP整体调度成本，同时能最大限度的提高风光等可再生能源的消纳率。
表2不同方案下VPP调度结果
Figure BDA0003751864390000211
(2)VPP参与电-绿证市场的交互分析：
为体现VPP电-绿证市场交互关系，图7给出了VPP等效负荷与GCT交易量之间的关系图。主要有以下特点：图中可以看出，由于高电价时段10-15h可控负荷大量响应，造成等效负荷大幅度降低，此时消纳权责较低，且在满足负荷需求下VPP可以尽可能的参与电市场，故此时绿色证书交易量以及电力市场交易量有较为明显的增加。而在15-18h，由于等效负荷此时较峰值电价时期有所升高，因此绿色证书交易量开始较少。由此可见，绿色证书交易量与电力市场交易量近似呈正相关，而与等效负荷近似呈负相关。
(3)碳-绿证市场灵敏度分析：
在电-碳-绿色证书联合交易机制下，VPP的总成本与碳排放单价、绿色证书单价密切相关，为了进一步明确VPP的总调度成本随多种市场交易价格变化规律，图8绘制了碳价-绿色证书价格变化下的VPP调度成本三维曲面图。随着碳价的增加，VPP成本呈现先增后减的趋势，这是因为在碳价较低的时段，VPP会优先选择从外部购入碳排放权，因此随着碳价的增加，VPP的成本也不断增加；当碳价达到某一临界值，即VPP从外部购买碳排放权的成本高于碳捕集的成本时，VPP选择利用碳捕集系统减少自身碳排放量，并可以对外出售部分碳排放权，因此随着碳价的继续增加，VPP碳市场的收益也不断增加，故总运行成本不断降低。
此外，随着绿色证书交易单价的不断增加，VPP在绿色证书交易市场上的收益也将不断增加，故VPP的调度成本也呈线性趋势下降。随着碳市场和绿色证书市场的不断开放，VPP可参考图8所得的三维曲面图，灵活选择报价合理的交易对象，从而降低自身运行成本，提高其运行效益。
(4)基于IGDT算法的VPP调度策略分析：
假设偏差因子最大值为20％，依次计算偏差因子从20％～0变化过程中扰动水平的数值，得到图9所示轨迹图。可以看到，随着VPP调度成本的逐渐增加，机会模型的扰动水平αa逐步降低，鲁棒模型的扰动水平αc逐渐增加，变化呈现相反的趋势。可以理解为，αc越高，风光的出力越偏离预测值，决策者对风光波动的看法越悲观，在此情况下所做出的决策就更加偏向保守，但同时风光容许波动范围也在不断扩大，这意味着系统调度方案和市场交易方案对不确定量的变化的敏感程度大大降低，因此系统鲁棒性会相对较高；相反，随着αa增加，风光的出力越高于所作预测，可再生能源能够更加多的发电，决策者对风光的波动更加的乐观，认为能够带来更多的收益，但所面临风险也会随之增加。
图10给出了偏差因子为10％时，电-碳-绿色证书市场交易结果。在鲁棒模型下，无论是哪一种类型的市场，VPP的调度策略几乎总是最为保守的，电、碳、绿色证书的交易出售量较低，但是购入值却较高。与之相反，机会模型下，VPP调度策略则愈发激进，表现为电、碳、绿色证书的出售量更高，且购入值更低。由此可见，本发明提出的IGDT模型能够准确对由风光出力不确定性引起的不同类型市场的交易策略的风险进行度量，从而为不同投资者提供决策参考。
本发明构建了电-碳-绿色证书联合交易下含风光、碳捕集、储能及需求响应等的VPP日前调度模型，并考虑可再生能源出力不确定性造成的潜在市场运行风险，建立了基于IGDT算法的VPP优化调度模型，经过算例分析，得出以下主要结论：
(1)通过电-碳市场联动、电-绿色证书市场联动，能够提高VPP市场交易策略的灵活性和整体的经济性；
(2)VPP的调度成本与基础碳价近似呈二次函数关系，而与绿色证书价格近似成本线性关系，随着碳-绿色证书市场的不断开放，VPP可灵活选择报价合理的交易对象，从而降低自身运行成本；
(3)随着VPP调度成本的增加，机会模型的扰动水平呈线性降低，鲁棒模型的扰动水平呈线性增加，VPP可以根据自身风险的偏好程度，灵活合理选择其在电-碳-绿证市场的交易策略。

Claims (6)

1.考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，其特征在于包括以下步骤：
步骤1：构建电-碳-绿色证书联合交易框架；
步骤2：构建联合交易框架下VPP目前优化调度模型；
步骤3：设置步骤2中优化调度模型的约束条件；
步骤4：对电-碳-绿证联合交易下VPP调度方案和市场交易策略进行风险度量；
步骤5：以VPP调度总成本最小为目标函数进行求解。
2.根据权利要求1所述考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，其特征在于：还包括步骤6：首先，对VPP分别在电-碳市场、电-绿证市场、碳-绿证市场的市场交易策略和调度方案进行分析；其次，由于VPP聚合风-光-火碳捕集系统，通过设置四种不同场景，对其在可再生能源消纳效果和二氧化碳排放的削减效果进行对比分析；最后，对于风光的不确定性，采用信息间隙决策理论对不确定参数进行处理，得出了风险规避模型和风险偏好模型约束下VPP的最优调度方案以及市场交易策略。
3.根据权利要求1所述考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，其特征在于：所述步骤2中，在电-碳-绿色证书联合交易框架下，以VPP运行成本最小为目标函数，构建联合交易机制下VPP目前优化调度模型，如式(1)所示，具体包含碳捕集电厂燃料成本、储能电池损耗成本、可控负荷成本、碳排放权交易成本、绿色证书交易成本和弃风、弃光惩罚成本：
Figure FDA0003751864380000011
式(1)中，F为VPP总成本，
Figure FDA0003751864380000012
为t时刻VPP的碳捕集电厂燃料成本、
Figure FDA0003751864380000013
为储能电池损耗成本、
Figure FDA0003751864380000014
为可控负荷成本、
Figure FDA0003751864380000015
为碳排放权交易成本、
Figure FDA0003751864380000016
为绿色证书交易成本、
Figure FDA0003751864380000017
为弃风、弃光惩罚费用；T表示调度周期；t表示时间；
(i)碳捕集电厂燃料成本
Figure FDA0003751864380000018
碳捕集电厂成本不考虑阀点效应，故总燃料费用计算式为：
Figure FDA0003751864380000019
式(2)中，
Figure FDA00037518643800000110
为t时刻VPP的碳捕集电厂燃料成本，a、b、c分别为火电机组的燃料成本系数，Pt G为火电机组t时刻的总出力；
(ii)储能系统的电池损耗费用
Figure FDA0003751864380000021
计算如下：
Figure FDA0003751864380000022
式(3)中，
Figure FDA0003751864380000023
为储能电池损耗成本，cess为电池损耗成本系数，Pt esc、Pt esd分别表示储能的充、放电功率；
(iii)可控负荷的调控费用
Figure FDA0003751864380000024
可控负荷的调控费用包括可转移负荷费用
Figure FDA0003751864380000025
可中断负荷的补偿费用
Figure FDA0003751864380000026
以及空调负荷调控引起的用户不舒适度补偿
Figure FDA0003751864380000027
Figure FDA0003751864380000028
Figure FDA0003751864380000029
Figure FDA00037518643800000210
Figure FDA00037518643800000211
式(4)～式(7)中，
Figure FDA00037518643800000212
可控负荷成本，
Figure FDA00037518643800000213
分别为t时刻中断负荷、转移负荷以及空调负荷的调控费用；cm、cs、cres分别为中断负荷、转移负荷及空调负荷的成本系数；
Figure FDA00037518643800000214
表示第m级中断负荷的响应量；Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷，Tt in表示室内温度，Tt ref表示参考室温；
(iv)碳交易成本
Figure FDA00037518643800000215
Figure FDA00037518643800000216
式(8)中，
Figure FDA00037518643800000217
为阶梯式碳交易成本；λ为碳交易基础价格，
Figure FDA00037518643800000218
为VPP实际的碳排放量，l为碳排放区间长度，α为价格增长率；
(v)绿色证书交易成本
Figure FDA00037518643800000219
Figure FDA00037518643800000220
式(9)中，
Figure FDA0003751864380000031
为绿色证书交易成本，Pt gc为绿色证书交易价格，
Figure FDA0003751864380000032
为VPP在t时刻购买或出售的绿色证书量，若为正值表示VPP此刻购买绿色证书，若为负值表示此刻出售绿色证书；
(vi)弃风、弃光惩罚成本
Figure FDA0003751864380000033
Figure FDA0003751864380000034
式(10)中，
Figure FDA0003751864380000035
为弃风、弃光惩罚成本，cpv、cpw分别为弃光、弃风惩罚费用，Pt ppv、Pt ppw分别为光伏、风电出力的预测值，Pt pv、Pt pw为实际消纳的光伏、风电功率。
4.根据权利要求1所述考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，其特征在于：所述步骤3中，VPP日前优化调度模型的约束条件包括：聚合单元的运行约束、功率平衡约束以及市场交易约束，具体包括以下部分：
(a)风-光-火碳捕集系统运行约束：
Pt GN＝Pt G-Pt GC (11)；
Pt GC＝Pt B+Pt OP (12)；
式(11)、式(12)中，Pt GN为碳捕集电厂的实际出力；Pt G为火电总出力；Pt GC为碳捕集的总能耗，包括固定能耗Pt B及运行能耗Pt OP；
将风光机组与碳捕集机组联合运行，风光的部分功率可以作为碳捕集系统的运行能耗使用，具体如下：
Figure FDA0003751864380000036
Figure FDA0003751864380000037
式(13)、式(14)中，Pt OP为运行能耗；Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率；wc为处理单位CO2需要的功耗；
Figure FDA0003751864380000038
为再生塔处理的CO2量；Et N为实际火电碳排放量；eg为单位火电机组出力产生的CO2量；Pt G为火电总出力；
风光机组除了供给碳捕集系统，余电依然能够由VPP向外部电网出售，表达式如下：
Pt pv＝Pt vc+Pt vw (15)；
Pt pw＝Pt wc+Pt ww (16)；
式(15)、式(16)中，Pt pv、Pt pw为实际消纳的光伏、风电功率；Pt vw、Pt ww分别为光伏以及风电的上网功率；Pt wc、Pt vc分别为风电、光伏用于碳捕集能耗的功率；
风、光、火电机组需满足最大出力约束：
0≤Pt pw≤Pt ppw (17)；
0≤Pt pv≤Pt ppv (18)；
Pgmin≤Pt G≤Pgmax (19)；
式(17)、式(18)、式(19)中，Pgmin、Pgmax分别为火电机组出力的技术最小值和最大值；PtG为火电总出力；Pt pv、Pt pw为实际消纳的光伏、风电功率；Pt ppv、Pt ppw分别为光伏、风电出力的预测值；此外，火电机组还须满足爬坡率约束：
Figure FDA0003751864380000041
式(20)中，Pt G为t时刻火电出力，
Figure FDA0003751864380000042
为t-1时刻火电出力，Pramp为火电机组的爬坡功率；
(b)储能装置运行约束：
Figure FDA0003751864380000043
式(21)中，
Figure FDA0003751864380000044
分别为t时刻、t-1储能蓄电量值；ηc、ηd为储能系统的充放电效率；Pt esc、Pt esd分别表示储能的充、放电功率；
(c)可控负荷运行约束：
Figure FDA0003751864380000045
Figure FDA0003751864380000046
Figure FDA0003751864380000047
式(22)、式(23)、式(24)中，
Figure FDA0003751864380000048
为第m级削减负荷的最大削减量，
Figure FDA0003751864380000049
为第m级可削减负荷，Pt iL为t时刻的总削减负荷，
Figure FDA00037518643800000410
为t-1时刻总削减负荷；M表示削减负荷的级别；Pcmax为最大连续削减负荷量；式(24)表示避免连续负荷削减导致用户舒适度下降；
Figure FDA00037518643800000411
Figure FDA00037518643800000412
Figure FDA0003751864380000051
Figure FDA0003751864380000052
式(25)、式(26)、式(27)、式(28)中，Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；Pt fmax为t时段最大转移负荷；
Figure FDA0003751864380000053
分别表示t时刻是否有负荷转移到别处或从别处转移过来，为布尔变量；
Figure FDA0003751864380000054
Tinmin≤Tt in≤Tinmax (30)
式(29)、式(30)中，γ、αt、β为与建筑物蓄热特性相关的参数；σ为空调的制冷能效比；Pt cold为中央空调的制冷功率；Tt in为t时刻温度，
Figure FDA0003751864380000056
为t-1时刻温度；Tinmin、Tinmax分别为满足用户舒适度的温度区间上、下限；
(d)系统功率平衡约束：
Pt cold+Pt loads+Pt esc+Pt mgs＝Pt vw+Pt ww+Pt GN+Pt mgb+Pt esd (31)
Pt loads＝Pt load-Pt sfl+Pt sfq-Pt iL (32)
式(31)、式(32)中，Pt cold为中央空调的制冷功率；Pt esc、Pt esd分别表示储能的充、放电功率；Pt loads为t时刻的等效负荷；Pt load为t时刻的基础负荷；Pt vw、Pt ww分别为光伏以及风电的上网功率；Pt mgs、Pt mgb分别为t时刻VPP向外部电网出售以及购入的电量；Pt GN为碳捕集电厂的实际出力；Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷；Pt iL为t时刻的总削减负荷；
(e)联络线传输功率约束：
0≤Pt mgs≤Pmgmax (33)
0≤Pt mgb≤Pmgmax (34)
式(33)、式(34)中，Pt mgmax表示联络线最大传输功率，Pt mgs为售电量，Pt mgb为购电量；
(f)碳排放市场运行约束：
采用基准线法确定VPP的碳排放配额：
Figure FDA0003751864380000055
式(35)中，
Figure FDA0003751864380000061
为VPP核发的碳配额总量，ηb、ηg为从市场购电及火电机组的配额系数，Pt mgb为购电量，Pt G为火电总出力；
VPP的实际碳排放量计算如下：
Figure FDA0003751864380000062
式(36)中，
Figure FDA0003751864380000063
为VPP的实际碳排放量，ap、bp、cP为实际排放系数，Pt mgb为购电量，
Figure FDA0003751864380000064
为碳捕集电厂净排放；
碳排放市场交易需求计算如下：
Figure FDA0003751864380000065
式(37)中，
Figure FDA0003751864380000066
为VPP碳排放市场交易量，
Figure FDA0003751864380000067
为VPP的实际碳排放量，
Figure FDA0003751864380000068
为VPP核发的碳配额总量；
(g)VPP消纳权责及绿证交易约束：
VPP内可再生能源的消纳量与绿色证书交易量总和应不小于消纳权责规定的电力消纳总量：
Figure FDA0003751864380000069
Figure FDA00037518643800000610
Figure FDA00037518643800000611
式(38)、式(39)、式(40)中，
Figure FDA00037518643800000612
表示可再生能源消纳总量；w为消纳权重；
Figure FDA00037518643800000613
为总电量；
Figure FDA00037518643800000614
为VPP在t时刻购买或出售的绿色证书量；Pt loads为t时刻的等效负荷；Pt cold为中央空调的制冷功率；Pt esc、Pt esd分别表示储能的充、放电功率；Pt pv、Pt pw分别为实际消纳的光伏、风电功率。
5.根据权利要求1所述考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法，其特征在于：所述步骤4中，采用信息间隙决策理论处理不确定性参数，具体如下：
按照预设目标的好坏，不确定性参数对目标函数的影响分为积极和消极两个方面，相应的IGDT模型又可分为机会以及鲁棒两种，其反应了不同投资者在面对风险时的价值取向；一种认为不确定性参数将会有利于目标函数实现，另一种则认为不确定性参数将会对目标函数产生负面影响，二者对应的数学模型分别如下：
Figure FDA0003751864380000071
Figure FDA0003751864380000072
式(41)、式(42)中，B0为预设目标值，由确定性的模型求解得出；Bc、Ba为不同风险偏向模型下的期望目标函数；βc、βa为偏差因子，表示预期成本与预设目标的偏差程度；H(X，d)＝0、G(X，d)≥0分别为优化模型的等式与不等式约束；α为不确定性参数的波动程度；
Figure FDA0003751864380000073
为不确定性参数的预测值；X为不确定参数实际值；d为决策变量，
Figure FDA0003751864380000074
为不确定参数波动范围函数，s.t.表示约束条件；
式(41)对应IGDT鲁棒模型，在最大化不确定性扰动下以期实现最低的预设目标；
式(42)对应IGDT机会模型，在最小化不确定性扰动下以期实现最高的预设目标。
6.一种联合交易机制下VPP目前优化调度模型，其特征在于：该模型以VPP运行成本最小为目标函数，包含碳捕集电厂燃料成本、储能电池损耗成本、可控负荷成本、碳排放权交易成本、绿色证书交易成本和弃风、弃光惩罚成本，如式(1)所示：
Figure FDA0003751864380000075
式(1)中，F为VPP总成本，
Figure FDA0003751864380000076
为t时刻VPP的碳捕集电厂燃料成本、
Figure FDA0003751864380000077
为储能电池损耗成本、
Figure FDA0003751864380000078
为可控负荷成本、
Figure FDA0003751864380000079
为碳排放权交易成本、
Figure FDA00037518643800000710
为绿色证书交易成本、
Figure FDA00037518643800000711
为弃风、弃光惩罚费用；T表示调度周期；t表示时间；
(i)碳捕集电厂燃料成本
Figure FDA00037518643800000712
碳捕集电厂成本不考虑阀点效应，故总燃料费用计算式为：
Figure FDA00037518643800000713
式(2)中，
Figure FDA0003751864380000081
为t时刻VPP的碳捕集电厂燃料成本，a、b、c分别为火电机组的燃料成本系数，Pt G为火电机组t时刻的总出力；
(ii)储能系统的电池损耗费用
Figure FDA0003751864380000082
计算如下：
Figure FDA0003751864380000083
式(3)中，
Figure FDA0003751864380000084
为储能电池损耗成本，cess为电池损耗成本系数，Pt esc、Pt esd分别表示储能的充、放电功率；
(iii)可控负荷的调控费用
Figure FDA0003751864380000085
可控负荷的调控费用包括可转移负荷费用
Figure FDA0003751864380000086
可中断负荷的补偿费用
Figure FDA0003751864380000087
以及空调负荷调控引起的用户不舒适度补偿
Figure FDA0003751864380000088
Figure FDA0003751864380000089
Figure FDA00037518643800000810
Figure FDA00037518643800000811
Figure FDA00037518643800000812
式(4)～式(7)中，
Figure FDA00037518643800000813
可控负荷成本，
Figure FDA00037518643800000814
分别为t时刻中断负荷、转移负荷以及空调负荷的调控费用；cm、cs、cres分别为中断负荷、转移负荷及空调负荷的成本系数；
Figure FDA00037518643800000815
表示第m级中断负荷的响应量；Pt sfl、Pt sfq分别表示转移至别的时段的负荷和从别的时段转移过来的负荷，Tt in表示室内温度，Tt ref表示参考室温；
(iv)碳交易成本
Figure FDA00037518643800000816
Figure FDA00037518643800000817
式(8)中，
Figure FDA00037518643800000818
为阶梯式碳交易成本；λ为碳交易基础价格，
Figure FDA00037518643800000819
为VPP实际的碳排放量，l为碳排放区间长度，α为价格增长率；
(v)绿色证书交易成本
Figure FDA0003751864380000091
Figure FDA0003751864380000092
式(9)中，
Figure FDA0003751864380000093
为绿色证书交易成本，Pt gc为绿色证书交易价格，
Figure FDA0003751864380000094
为VPP在t时刻购买或出售的绿色证书量，若为正值表示VPP此刻购买绿色证书，若为负值表示此刻出售绿色证书；
(vi)弃风、弃光惩罚成本
Figure FDA0003751864380000095
Figure FDA0003751864380000096
式(10)中，
Figure FDA0003751864380000097
为弃风、弃光惩罚成本，cpv、cpw分别为弃光、弃风惩罚费用，Pt ppv、Pt ppw分别为光伏、风电出力的预测值，Pt pv、Pt pw为实际消纳的光伏、风电功率。
CN202210991160.1A 2022-07-18 2022-07-18 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法 Pending CN115660309A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210991160.1A CN115660309A (zh) 2022-07-18 2022-07-18 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210991160.1A CN115660309A (zh) 2022-07-18 2022-07-18 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法

Publications (1)

Publication Number Publication Date
CN115660309A true CN115660309A (zh) 2023-01-31

Family

ID=84983809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210991160.1A Pending CN115660309A (zh) 2022-07-18 2022-07-18 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法

Country Status (1)

Country Link
CN (1) CN115660309A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116109335A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司 一种基于数字孪生的新能源动态电价数据管理系统
CN116128553A (zh) * 2023-04-19 2023-05-16 南京师范大学 基于绿证和碳交易交互的综合能源调度方法、系统
CN117422249A (zh) * 2023-10-20 2024-01-19 国网甘肃省电力公司超高压公司 综合能源系统日前优化调度方法
CN117439080A (zh) * 2023-12-20 2024-01-23 国网山东省电力公司营销服务中心(计量中心) 一种虚拟电厂的调度方法、系统、存储介质及设备
CN117458482A (zh) * 2023-11-21 2024-01-26 华北电力大学 一种虚拟电厂与风光储供需匹配的绿电交易方法
CN117689234A (zh) * 2024-02-04 2024-03-12 山东科技大学 基于多主体双层博弈的园区综合能源系统调度方法及系统
CN117933445A (zh) * 2023-11-27 2024-04-26 国网青海省电力公司清洁能源发展研究院 一种农村虚拟电厂售电商系统的双层优化方法和装置
CN118539435A (zh) * 2024-07-24 2024-08-23 国网山东省电力公司莱芜供电公司 计及电碳耦合约束的电力系统可用灵活性提升优化方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116109335A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司 一种基于数字孪生的新能源动态电价数据管理系统
CN116109335B (zh) * 2023-04-10 2023-09-08 国网浙江省电力有限公司 一种基于数字孪生的新能源动态电价数据管理系统
CN116128553A (zh) * 2023-04-19 2023-05-16 南京师范大学 基于绿证和碳交易交互的综合能源调度方法、系统
CN117422249A (zh) * 2023-10-20 2024-01-19 国网甘肃省电力公司超高压公司 综合能源系统日前优化调度方法
CN117458482A (zh) * 2023-11-21 2024-01-26 华北电力大学 一种虚拟电厂与风光储供需匹配的绿电交易方法
CN117933445A (zh) * 2023-11-27 2024-04-26 国网青海省电力公司清洁能源发展研究院 一种农村虚拟电厂售电商系统的双层优化方法和装置
CN117439080A (zh) * 2023-12-20 2024-01-23 国网山东省电力公司营销服务中心(计量中心) 一种虚拟电厂的调度方法、系统、存储介质及设备
CN117439080B (zh) * 2023-12-20 2024-04-12 国网山东省电力公司营销服务中心(计量中心) 一种虚拟电厂的调度方法、系统、存储介质及设备
CN117689234A (zh) * 2024-02-04 2024-03-12 山东科技大学 基于多主体双层博弈的园区综合能源系统调度方法及系统
CN117689234B (zh) * 2024-02-04 2024-05-03 山东科技大学 基于多主体双层博弈的园区综合能源系统调度方法及系统
CN118539435A (zh) * 2024-07-24 2024-08-23 国网山东省电力公司莱芜供电公司 计及电碳耦合约束的电力系统可用灵活性提升优化方法

Similar Documents

Publication Publication Date Title
CN115660309A (zh) 考虑电-碳-绿证联合交易的虚拟电厂日前优化调度方法
Wang et al. Demand response comprehensive incentive mechanism-based multi-time scale optimization scheduling for park integrated energy system
Wang et al. Optimal scheduling of the regional integrated energy system considering economy and environment
CN111738497B (zh) 计及需求侧响应的虚拟电厂双层优化调度方法
Liu et al. Energy management method of integrated energy system based on collaborative optimization of distributed flexible resources
CN108009693A (zh) 基于两级需求响应的并网微电网双层优化方法
CN104376385A (zh) 一种微电网电价优化方法
CN107194514B (zh) 一种针对风电功率预测误差的需求响应多时间尺度调度方法
CN115115277A (zh) 园区综合能源调度方法、装置、计算机设备及存储介质
CN113746089B (zh) 面向多用户的多时间尺度电力套餐与家庭能量优化方法
CN113869927B (zh) 一种促进社区微网内部多产消者能源共享的分时定价方法
CN111210079A (zh) 一种分布式能源虚拟电厂运行优化方法及系统
CN115293485A (zh) 考虑电动汽车和需求响应的综合能源系统低碳调度方法
CN116706920A (zh) 虚拟电厂参与电能量调峰市场的最优策略与优化调度方法
CN116961044A (zh) 基于模糊机会约束规划广义共享储能优化配置方法及系统
CN115689166A (zh) 一种区域分布式能源资源聚合利用方法及系统
CN116187648A (zh) 基于热电解耦的虚拟电厂热电联合优化调度方法
CN113690875B (zh) 一种微网实时互动等值模型建立方法
Guo et al. How to realize the power demand side actively matching the supply side?——A virtual real-time electricity prices optimization model based on credit mechanism
CN117526451A (zh) 考虑柔性负荷下的区域综合能源系统配置优化方法
CN115438906A (zh) 一种多虚拟电厂点对点交易方法、电子设备及存储介质
CN114282708A (zh) 考虑多尺度需求响应的跨区域综合能源系统优化运行方法及系统
Zhong et al. Two Stage Optimization Model Considering Demand Response Dispatch Value of Thermoelectric Load
CN118100175B (zh) 一种基于分布式优化的能源系统能量管理方法
CN115563816B (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