申请/专利权人：广西电网有限责任公司电力科学研究院

申请日：2018-12-17

公开（公告）日：2020-11-24

公开（公告）号：CN109993639B

主分类号：G06Q30/08(20120101)

分类号：G06Q30/08(20120101);G06Q10/04(20120101);G06Q50/06(20120101)

优先权：

专利状态码：有效-授权

法律状态：2020.11.24#授权;2019.08.02#实质审查的生效;2019.07.09#公开

摘要：本发明涉及电力系统及其自动化技术领域，具体涉及一种虚拟电厂参与日前电力市场的优化投标策略。本发明针对虚拟电厂参与日前电力市场投标研究，建立了优化投标模型，提出了优化投标策略，为今后虚拟电厂参与日前电力市场提供了科学理论支持。所建立的成本容量模型，充分考虑了虚拟电厂内部多类型不同需求侧资源出力特性，生成的优化投标策略计及了投标主体内部不同需求侧资源在日前市场投标电量和投标电价，体现了投标策略的优越性。

主权项：1.一种虚拟电厂参与日前电力市场的优化投标方法，其特征在于：包括以下步骤：S1：在各虚拟电厂提交其需求侧资源的投标电量和电价后，以最小成本为优化目标构建调度中心市场出清模型，获取市场出清价格和出清电量；所述步骤S1中的构建调度中心市场出清模型的具体步骤如下：S11：在各虚拟电厂提交其需求侧资源i的投标电量Pi和投标电价λi后，调度中心进行经济调度获取出清电价πd和出清电量gi，假设出清电量gi≥0，其目标函数如下： 式中，I表示所有需求侧资源集合；S12：负荷供需平衡约束： 式中，d为负荷；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；S13：出清电量受到投标电量约束如下：gi≤Pi,i∈Iπgi3；式中，πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益；S2：根据步骤S1得到的市场出清价格，结合虚拟电厂内部需求侧资源的变动运行成本和投标电量，获得该虚拟电厂净收益函数；获得该虚拟电厂净收益函数的计算方式为：设该虚拟电厂为虚拟电厂A，虚拟电厂A的总收入减去运行成本即得净收益R，具体如下： 式中，ci表示需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量；S3：以虚拟电厂最大净收益为优化目标，结合其内部各需求侧资源最大容量约束，建立虚拟电厂的投标优化模型；建立虚拟电厂的投标优化模型的步骤具体如下：S31：设该虚拟电厂为虚拟电厂A，虚拟电厂A投标优化目标是其净收益最大化，具体如下： 式中，ci表示虚拟电厂A的需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量；S32：虚拟电厂A的需求侧资源i投标电量具有出力约束，具体如下：Pi≤Pimaxi∈A6；式中，Pi为虚拟电厂A的需求侧资源i的投标电量，Pimax表示虚拟电厂A的需求侧资源i的最大容量；S4：根据步骤S1获得的市场出清模型的最优性条件，将其作为虚拟电厂投标优化的部分约束条件，从而将步骤S3建立的投标优化模型转化成具有平衡约束的数学规划模型；所述步骤S4中的将步骤S3建立的投标优化模型转化成平衡约束的数学规划模型的步骤具体如下：S41：根据对偶理论，获取市场出清模型的最优性条件，具体如式7-16，将其加入到步骤S3建立的虚拟电厂A的投标优化模型的约束条件中，得到虚拟电厂A具有平衡约束的数学规划模型； gi-Pi≤0,i∈A8； πd+πgi-λi≤0,i∈A10； πgi≤0,i∈I12；πd+πgi-λigi＝0,i∈A13； Pi-giπgi＝0,i∈A15； 其中，式7-9是步骤S1的构建的调度中心市场出清模型的约束条件，式10-12是对偶问题的约束条件，式13-16是互补条件；式中，为虚拟电厂A的预测市场负荷；表示虚拟电厂A预测其余需求侧资源的投标电量-电价；Pi为虚拟电厂A的需求侧资源i的投标电量，gi为出清电量；λi为投标电价；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益；表示市场其余需求侧资源集合；S42：根据gi≥0,i∈I，对公式10和11进行变形，得到：πd+πgi-λigi≤0,i∈A17； S43：根据公式12，对公式8和9进行变形，得到：gi-Piπgi≤0,i∈A19； S44：根据公式17-20的非正性，得到互补条件13-16的等效公式： S45：使用公式7:重新安排公式21各项，得到等效公式22，即原对偶等式条件为： S46：至此，生成步骤S3建立的虚拟电厂A的投标优化模型的目标函数及其约束条件，步骤S3建立的虚拟电厂A的投标优化模型的目标函数如式5所示，其约束条件为式6-12和式22；S5：采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项，即连续变量乘积，形成混合整数规划模型；所述步骤S5中采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项的具体步骤如下：S51：应用式13和15的互补条件，用λigi和πgiPi重写目标函数5中的非线性项πdgi：πdgi＝λigi-πgiPi23；S52：采用二进制扩充方法处理约束条件22中的非线性项λigi：首先，通过一组离散值{λim,m＝0,1,...,M1}近似连续决策变量λi，其近似表达式为： 这里，λi的取值区间为λimin,λimax，Δλi＝λimax-λiminM1，xki为二进制变量；然后，将gi乘以等式24两侧，得到新变量zki＝xkigi，从而得到： 同时，引入相应约束：0≤gi-zki≤G11-xki26；0≤zki≤G1xki27；式中，G1是一个常数，S53：采用二进制扩充方法，处理约束条件22中的非线性项πgiPi：首先，通过一组离散值{Pim,m＝0,1,...,M2}近似连续决策变量Pi，其近似表达式为 这里，Pi的取值区间为Pimin,Pimax，ΔPi＝Pimax-PiminM2，yki为二进制变量；然后，将πgi乘以等式28两侧，得到新变量wki＝πgiyki，从而得到： 同时，引入相应约束0≤wki-πgi≤G21-yki30；0≤-wki≤G2yki31；式中，G2是一个常数，S54：在目标函数5，约束条件6-12和22应用上述步骤S51-S53的线性转化形式，得到混合整数规划模型： πgi≤0,i∈I40； 0≤gi-zki≤G11-xki,i∈A42；0≤zki≤G1xki,i∈A43；0≤wki-πgi≤G21-yki,i∈A44；0≤-wki≤G2yki,i∈A45；xki,yki∈{0,1},i∈A46；S6：结合虚拟电厂内部需求侧资源的出力特性，确定虚拟电厂参与日前市场的优化投标模型；具体步骤如下：S61：确定虚拟电厂A所参与日前市场的时间间隔为Δt，时段数S62：根据虚拟电厂A和其他虚拟电厂内部需求侧资源组成情况，确定市场中所有需求侧资源种类、数目以及所属虚拟电厂；虚拟电厂A的需求侧资源包括风电机组、燃气轮机、电动汽车和柔性负荷；S63：根据历史信息和市场信息，虚拟电厂A预测得到市场负荷和其余需求侧资源投标电量和电价；S64：根据虚拟电厂A内部需求侧资源特性，构建其在t时段容量成本模型CA,t；S641：虚拟电厂A在t时段容量成本模型CA,t具体包括风电机组WT成本CA,WT,t、燃气轮机GT成本CA,GT,t、柔性负荷PSDR成本CA,PSDR,t和电动汽车EV成本CA,EV,t，即：CA,t＝CA,WT,t+CA,GT,t+CA,PSDR,t+CA,EV,t；47①风电机组成本由于风电机组发电成本为零，这里风电机组成本主要是弃风成本：CA,WT,t＝cWTPA,WT,t48； 式中，cWT表示弃风电惩罚费用；PA,WT,t表示在t时刻弃风电；表示虚拟电厂A预测风电机组在时刻t发电功率；②燃气轮机成本燃气轮机成本由燃料费用、启停成本和碳排放惩罚成本组成，即： 其中：ηA,GT,t＝aA,GTPA,GT,t2+bA,GTPA,GT,t+cA,GT51； 式中，FA,GT,t为GT燃气耗量；LHV为低位热值；ηA,GT,t分别为GT发电效率；分别为GT启停成本；DK、YK、VK分别为第k类污染气体的排放量及单位排放量对应的罚款和环境价值；PA,GT,t为GT发电功率；aA,GT、bA,GT、cA,GT为GT效率系数；分别为GT每次开机、停机费用；uA,GT,t为二进制变量，分别为GT机组在时段t工作状态；GT的相关约束条件为： 式中，表示GT在时刻t的最大出力；③柔性负荷成本柔性负荷成本是需求响应负荷削减或增加产生的费用，即：CA,PSDR,t＝aA,PSDRPA,PSDR,t2+bA,PSDRPA,PSDR,t+cA,PSDR,t55；式中，PA,PSDR,t表示PSDR在时刻t负荷削减量或增加量；aA,PSDR、bA,PSDR、cA,PSDR,t表示与价格弹性系数及电价折扣率有关系数；PSDR的相关约束条件如下： 式中：表示PSDR在时刻t最大削减或增加负荷；和分别表示PSDR在时刻t削减负荷和增加负荷；和分别表示PSDR在时刻t是否处于削减负荷状态和增加负荷状态；④电动汽车成本电动汽车成本是电池损耗成本，即： 式中：nA,EV表示EV数量；cA,EV,i、EA,EV,i、NA,EV,i、DODA,EV,i分别表示第i辆EV的电池购买成本、电池容量、电池可用循环次数及电池可用放电限度；为EV的放电效率；LA,EV,i,t为第i辆EV在时刻t的行驶距离；SA,EV,i表示第i辆EV能量需求；为决策变量，表示EV的放电功率；EV的相关约束条件具体如下： 式中，EA,EV,i,t表示第i辆EV在时刻t的蓄电量；分别表示第i辆EV蓄电量上、下限；表示第i辆EV在时刻t的充电功率；分别表示第i辆EV充放功率上限；分别表示第i辆EV在时刻t是否处于充放状态；表示第i辆EV在时刻t是否处于接入电网状态；分别表示第i辆EV始、末时刻蓄电量；分别表示第i辆EV充电效率；S642：应用分段线性化方法，处理GT容量成本模型及约束条件；具体是将GT在t时刻燃料费用与出力函数分成nGT段，每一段引入一个状态变量τi,GT,t和一个连续变量pi,GT,t，其燃料费用的线性化表达式如下： 相关约束条件为： pi,GTτi,t≤pi,GT,t≤pi+1,GTτi,GT,t72； 对于GT启停成本模型，可以重写为线性形式： S643：应用分段线性化方法，处理PSDR成本容量模型及相关约束条件；具体是将PSDR在t时刻增加负荷费用与出力函数分成段，每一段引入一个状态变量和一个连续变量其增加负荷费用的线性化表达式如下： 相关约束条件为： 将PSDR在t时刻削减负荷费用与出力函数分成段，每一段引入一个状态变量和一个连续变量其削减负荷费用的线性化表达式如下： 相关约束条件为： S65：根据成本容量模型及相关约束，结合预测信息，构建虚拟电厂A参加日前电力市场优化投标模型；S7：采用YamipGurobi软件计算虚拟电厂参与日前电力市场的最大收益投标模型，获取投标电量和投标电价。

全文数据：一种虚拟电厂参与日前电力市场的优化投标策略技术领域本发明涉及电力系统及其自动化技术领域，具体涉及一种虚拟电厂参与日前电力市场的优化投标策略。背景技术近年来，风电、太阳能等间歇性能源发展迅猛，随着大规模可再生能源接入，其随机性和间歇性特点使得电力系统安全稳定运行带来一定挑战。燃气轮机、电动汽车、柔性负荷等需求侧资源因其灵活、经济、可靠、环保特点，能有效解决上述问题，因此得到迅速发展。但需求侧资源一般容量较小，且多以分散形式存在，在电力市场上不具有竞争力，亟需协调多类型需求侧资源管理技术出现。虚拟电厂能有效聚合分散的不同容量等级需求侧资源，通过协调优化控制，能够提高系统整体的稳定性，减轻调度负担，同时作为整体参与电力市场，具有更强的市场竞争力，为用户带来可观收益，有利于调动电力用户参与系统运行的积极性实现资源的高效利用。在需求响应项目中，虚拟电厂聚合的广义需求侧资源主要包括分布式发电distributedgeneration，DG资源、储能资源和负荷资源。其中，DG资源包括2类：一类式是诸如风能、太阳能等间歇式可再生能源，另一类是可控DG资源，如燃料电池、微型燃气轮机、柴油发电机等。储能资源包括静态储能和电动汽车。负荷资源包括不可平移负荷、可平移负荷和可削减负荷等。现在研究关于虚拟电厂参与日前电力市场投标优化，缺少对协调多类型不同需求侧资源出力特性的刻画，没有充分考虑投标主体内部不同需求侧资源在日前市场投标电量和投标电价的获取，没有实现合理协调需求侧资源出力，达到利益最大化。发明内容为了解决上述问题，本发明提供了一种虚拟电厂参与日前电力市场的优化投标策略，具体技术方案如下：一种虚拟电厂参与日前电力市场的优化投标策略包括以下步骤：S1：在各虚拟电厂提交其需求侧资源的投标电量和电价后，以最小成本为优化目标构建调度中心市场出清模型，获取市场出清价格和出清电量；S2：根据步骤S1得到的市场出清价格，结合虚拟电厂内部需求侧资源的变动运行成本和投标电量，获得该虚拟电厂净收益函数；S3：以虚拟电厂最大净收益为优化目标，结合其内部各需求侧资源最大容量约束，建立虚拟电厂的投标优化模型；S4：根据步骤S1获得的市场出清模型的最优性条件，将其作为虚拟电厂投标优化的部分约束条件，从而将步骤S3建立的投标优化模型转化成具有平衡约束的数学规划模型；S5：采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项，即连续变量乘积，形成混合整数规划模型；S6：结合虚拟电厂内部需求侧资源的出力特性，确定虚拟电厂参与日前市场的优化投标模型；S7：采用YamipGurobi软件计算虚拟电厂参与日前电力市场的最大收益投标模型，获取投标电量和投标电价。优选地，所述步骤S1中的构建调度中心市场出清模型的具体步骤如下：S11：在各虚拟电厂提交其需求侧资源i的投标电量Pi和投标电价λi后，调度中心进行经济调度获取出清电价πd和出清电量gi，假设出清电量gi≥0，其目标函数如下：式中，I表示所有需求侧资源集合；S12：负荷供需平衡约束：式中，d为负荷；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；S13：出清电量受到投标电量约束如下：gi≤Pi,i∈Iπgi3；式中，πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益。优选地，所述步骤S2中的获得该虚拟电厂净收益函数的计算方式为：设该虚拟电厂为虚拟电厂A，虚拟电厂A的总收入减去运行成本即得净收益R，具体如下：式中，ci表示需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。优选地，所述步骤S3中建立虚拟电厂的投标优化模型的步骤具体如下：S31：设该细腻电厂为虚拟电厂A，虚拟电厂A投标优化目标是其净收益最大化，具体如下：式中，ci表示虚拟电厂A的需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。S32：虚拟电厂A的需求侧资源i投标电量具有出力约束，具体如下：Pi≤Pimaxi∈A6；式中，Pi为虚拟电厂A的需求侧资源i的投标电量，Pimax表示虚拟电厂A的需求侧资源i的最大容量。优选地，所述步骤S4中的将步骤S3建立的投标优化模型转化成平衡约束的数学规划模型的步骤具体如下：S41：根据对偶理论，获取市场出清模型的最优性条件，具体如式7-16，将其加入到步骤S3建立的虚拟电厂A的投标优化模型的约束条件中，得到虚拟电厂A具有平衡约束的数学规划模型；gi-Pi≤0,i∈A8；πd+πgi-λi≤0,i∈A10；πgi≤0,i∈I12；πd+πgi-λigi＝0,i∈A13；Pi-giπgi＝0,i∈A15；其中，式7-9是步骤S1的构建的调度中心市场出清模型的约束条件，式10-12是对偶问题的约束条件，式13-16是互补条件；式中，为虚拟电厂A的预测市场负荷；表示虚拟电厂A预测其余需求侧资源的投标电量-电价；Pi为虚拟电厂A的需求侧资源i的投标电量，gi为出清电量；λi为投标电价；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益；表示市场其余需求侧资源集合；S42：根据gi≥0,i∈I，对公式10和11进行变形，得到：πd+πgi-λigi≤0,i∈A17；S43：根据公式12，对公式8和9进行变形，得到：gi-Piπgi≤0,i∈A19；S44：根据公式17-20的非正性，得到互补条件13-16的等效公式：S45：使用公式7:重新安排公式21各项，得到等效公式22，即原对偶等式条件为：S46：至此，生成步骤S3建立的虚拟电厂A的投标优化模型的目标函数及其约束条件，步骤S3建立的虚拟电厂A的投标优化模型的目标函数如式5所示，其约束条件为式6-12和式22。优选地，所述步骤S5中采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项的具体步骤如下：S51：应用式13和15的互补条件，用λigi和πgiPi重写目标函数5中的非线性项πdgi：πdgi＝λigi-πgiPi23；S52：采用二进制扩充方法处理约束条件22中的非线性项λigi：首先，通过一组离散值{λim,m＝0,1,...,M1}近似连续决策变量λi，其近似表达式为：这里，λi的取值区间为xki为二进制变量；然后，将gi乘以等式24两侧，得到新变量zki＝xkigi，从而得到：同时，引入相应约束：0≤gi-zki≤G11-xki26；0≤zki≤G1xki27；式中，G1是一个常数，通常S53：采用二进制扩充方法，处理约束条件22中的非线性项πgiPi：首先，通过一组离散值{Pim,m＝0,1,...,M2}近似连续决策变量Pi，其近似表达式为这里，Pi的取值区间为Pimin,Pimax，ΔPi＝Pimax-PiminM2，yki为二进制变量；然后，将πgi乘以等式28两侧，得到新变量wki＝πgiyki，从而得到：同时，引入相应约束0≤wki-πgi≤G21-yki30；0≤-wki≤G2yki31；式中，G2是一个常数，通常S54：在目标函数5，约束条件6-12和22应用上述步骤S51-S53的线性转化形式，得到混合整数规划模型：πgi≤0,i∈I40；0≤gi-zki≤G11-xki,i∈A42；0≤zki≤G1xki,i∈A43；0≤wki-πgi≤G21-yki,i∈A44；0≤-wki≤G2yki,i∈A45；xki,yki∈{0,1},i∈A46。优选地，所述步骤S6中确定虚拟电厂参与日前市场的优化投标模型的具体步骤如下：S61：确定虚拟电厂A所参与日前市场的时间间隔为Δt，时段数S62：根据虚拟电厂A和其他虚拟电厂内部需求侧资源组成情况，确定市场中所有需求侧资源种类、数目以及所属虚拟电厂；S63：根据历史信息和市场信息，虚拟电厂A预测得到市场负荷和其余需求侧资源投标电量和电价；S64：根据虚拟电厂A内部需求侧资源特性，构建其在t时段容量成本模型CA,t；S65：根据成本容量模型及相关约束，结合预测信息，构建虚拟电厂A参加日前电力市场优化投标模型。优选地，所述Δt设置为1h，则时段数＝241＝24时段。本发明的有益效果为：本发明针对虚拟电厂参与日前电力市场投标研究，建立了优化投标模型，提出了优化投标策略，为今后虚拟电厂参与日前电力市场提供了科学理论支持。所建立的成本容量模型，充分考虑了虚拟电厂内部多类型不同需求侧资源出力特性，生成的优化投标策略计及了投标主体内部不同需求侧资源在日前市场投标电量和投标电价，体现了投标策略的优越性。具体实施方式为了更好的理解本发明，下面具体实施例对本发明作进一步说明：一种虚拟电厂参与日前电力市场的优化投标策略包括以下步骤：S1：在各虚拟电厂提交其需求侧资源的投标电量和电价后，以最小成本为优化目标构建调度中心市场出清模型，获取市场出清价格和出清电量；构建调度中心市场出清模型的具体步骤如下：S11：在各虚拟电厂提交其需求侧资源i的投标电量Pi和投标电价λi后，调度中心进行经济调度获取出清电价πd和出清电量gi，假设出清电量gi≥0，其目标函数如下：式中，I表示所有需求侧资源集合；S12：负荷供需平衡约束：式中，d为负荷；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；S13：出清电量受到投标电量约束如下：gi≤Pi,i∈Iπgi3；式中，πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益。S2：根据步骤S1得到的市场出清价格，结合虚拟电厂内部需求侧资源的变动运行成本和投标电量，获得该虚拟电厂净收益函数；获得该虚拟电厂净收益函数的计算方式为：设该虚拟电厂为虚拟电厂A，虚拟电厂A的总收入减去运行成本即得净收益R，具体如下：式中，ci表示需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。S3：以虚拟电厂最大净收益为优化目标，结合其内部各需求侧资源最大容量约束，建立虚拟电厂的投标优化模型；建立虚拟电厂的投标优化模型的步骤具体如下：S31：设该细腻电厂为虚拟电厂A，虚拟电厂A投标优化目标是其净收益最大化，具体如下：式中，ci表示虚拟电厂A的需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。S32：虚拟电厂A的需求侧资源i投标电量具有出力约束，具体如下：Pi≤Pimaxi∈A6；式中，Pi为虚拟电厂A的需求侧资源i的投标电量，Pimax表示虚拟电厂A的需求侧资源i的最大容量。S4：根据步骤S1获得的市场出清模型的最优性条件，将其作为虚拟电厂投标优化的部分约束条件，从而将步骤S3建立的投标优化模型转化成具有平衡约束的数学规划模型；具体如下：S41：根据对偶理论，获取市场出清模型的最优性条件，具体如式7-16，将其加入到步骤S3建立的虚拟电厂A的投标优化模型的约束条件中，得到虚拟电厂A具有平衡约束的数学规划模型；gi-Pi≤0,i∈A8；πd+πgi-λi≤0,i∈A10；πgi≤0,i∈I12；πd+πgi-λigi＝0,i∈A13；Pi-giπgi＝0,i∈A15；其中，式7-9是步骤S1的构建的调度中心市场出清模型的约束条件，式10-12是对偶问题的约束条件，式13-16是互补条件；式中，为虚拟电厂A的预测市场负荷；表示虚拟电厂A预测其余需求侧资源的投标电量-电价；Pi为虚拟电厂A的需求侧资源i的投标电量，gi为出清电量；λi为投标电价；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益；表示市场其余需求侧资源集合；S42：根据gi≥0,i∈I，对公式10和11进行变形，得到：πd+πgi-λigi≤0,i∈A17；S43：根据公式12，对公式8和9进行变形，得到：gi-Piπgi≤0,i∈A19；S44：根据公式17-20的非正性，得到互补条件13-16的等效公式：S45：使用公式7:重新安排公式21各项，得到等效公式22，即原对偶等式条件为：S46：至此，生成步骤S3建立的虚拟电厂A的投标优化模型的目标函数及其约束条件，步骤S3建立的虚拟电厂A的投标优化模型的目标函数如式5所示，其约束条件为式6-12和式22。S5：采用二进制扩充方法abinaryexpansionapproach处理步骤S4中的具有平衡约束的数学规划模型中的非线性项，即连续变量乘积，形成混合整数规划模型；具体步骤如下：S51：应用式13和15的互补条件，用λigi和πgiPi重写目标函数5中的非线性项πdgi：πdgi＝λigi-πgiPi23；S52：采用二进制扩充方法处理约束条件22中的非线性项λigi：首先，通过一组离散值{λim,m＝0,1,...,M1}近似连续决策变量λi，其近似表达式为：这里，λi的取值区间为xki为二进制变量；然后，将gi乘以等式24两侧，得到新变量zki＝xkigi，从而得到：同时，引入相应约束：0≤gi-zki≤G11-xki26；0≤zki≤G1xki27；式中，G1是一个常数，通常S53：采用二进制扩充方法，处理约束条件22中的非线性项πgiPi：首先，通过一组离散值{Pim,m＝0,1,...,M2}近似连续决策变量Pi，其近似表达式为这里，Pi的取值区间为Pimin,Pimax，ΔPi＝Pimax-PiminM2，yki为二进制变量；然后，将πgi乘以等式28两侧，得到新变量wki＝πgiyki，从而得到：同时，引入相应约束0≤wki-πgi≤G21-yki30；0≤-wki≤G2yki31；式中，G2是一个常数，通常S54：在目标函数5，约束条件6-12和22应用上述步骤S51-S53的线性转化形式，得到混合整数规划模型：πgi≤0,i∈I40；0≤gi-zki≤G11-xki,i∈A42；0≤zki≤G1xki,i∈A43；0≤wki-πgi≤G21-yki,i∈A44；0≤-wki≤G2yki,i∈A45；xki,yki∈{0,1},i∈A46。S6：结合虚拟电厂内部需求侧资源的出力特性，确定虚拟电厂参与日前市场的优化投标模型；具体步骤如下：S61：确定虚拟电厂A所参与日前市场的时间间隔为Δt，时段数本实施例中Δt设置为1h，则时段数＝241＝24时段。S62：根据虚拟电厂A和其他虚拟电厂内部需求侧资源组成情况，确定市场中所有需求侧资源种类、数目以及所属虚拟电厂。以虚拟电厂A、B和C参加区域日前市场投标为例进行介绍。虚拟电厂A由4类需求侧资源组成，包括风电机组windturbine,PT、燃气轮机gasturbine，GT、电动汽车electricvehicle，EV和柔性负荷price-sensitivedemandresponse，PSDR。虚拟电厂B由GT和PSDR组成。虚拟电厂C由GT和EV组成。此处PSDR是指基于激励的需求响应负荷。S63：根据历史信息和市场信息，虚拟电厂A预测得到市场负荷和其余需求侧资源投标电量和电价。包括虚拟电厂B的投标信息和虚拟电厂C的投标信息此处PSDR+表示增加需求响应负荷，而PSDR-表示减少需求响应负荷。S64：根据虚拟电厂A内部需求侧资源特性，构建其在t时段容量成本模型CA,t。S641：虚拟电厂A在t时段容量成本模型CA,t具体包括WT成本CA,WT,t、GT成本CA,GT,t、PSDR成本CA,PSDR,t和EV成本CA,EV,t，即：CA,t＝CA,WT,t+CA,GT,t+CA,PSDR,t+CA,EV,t；47①WT成本由于WT发电成本为零，这里WT成本主要是弃风成本：CA,WT,t＝cWTPA,WT,t48；式中，cWT表示弃风电惩罚费用；PA,WT,t表示在t时刻弃风电；表示虚拟电厂A预测风电机组在时刻t发电功率。②GT成本GT成本由燃料费用、启停成本和碳排放惩罚成本组成，即：其中：ηA,GT,t＝aA,GTPA,GT,t2+bA,GTPA,GT,t+cA,GT51；式中，FA,GT,t为GT燃气耗量；LHV为低位热值；ηA,GT,t分别为GT发电效率；分别为GT启停成本；DK、YK、VK分别为第k类污染气体的排放量及单位排放量对应的罚款和环境价值；PA,GT,t为GT发电功率；aA,GT、bA,GT、cA,GT为GT效率系数；分别为GT每次开机、停机费用；uA,GT,t为二进制变量，分别为GT机组在时段t工作状态。GT的相关约束条件为：式中，表示GT在时刻t的最大出力。③PSDR成本PSDR成本是需求响应负荷削减或增加产生的费用，即：CA,PSDR,t＝aA,PSDRPA,PSDR,t2+bA,PSDRPA,PSDR,t+cA,PSDR,t55；式中，PA,PSDR,t表示PSDR在时刻t负荷削减量或增加量；aA,PSDR、bA,PSDR、cA,PSDR,t表示与价格弹性系数及电价折扣率有关系数。PSDR的相关约束条件如下：式中：表示PSDR在时刻t最大削减或增加负荷；和分别表示PSDR在时刻t削减负荷和增加负荷；和分别表示PSDR在时刻t是否处于削减负荷状态和增加负荷状态。④EV成本EV成本是电池损耗成本，即：式中：nA,EV表示EV数量；cA,EV,i、EA,EV,i、NA,EV,i、DODA,EV,i分别表示第i辆EV的电池购买成本、电池容量、电池可用循环次数及电池可用放电限度；为EV的放电效率；LA,EV,i,t为第i辆EV在时刻t的行驶距离；SA,EV,i表示第i辆EV能量需求；为决策变量，表示EV的放电功率。EV的相关约束条件具体如下：式中，EA,EV,i,t表示第i辆EV在时刻t的蓄电量；分别表示第i辆EV蓄电量上、下限；表示第i辆EV在时刻t的充电功率；分别表示第i辆EV充放功率上限；分别表示第i辆EV在时刻t是否处于充放状态；表示第i辆EV在时刻t是否处于接入电网状态；分别表示第i辆EV始、末时刻蓄电量；分别表示第i辆EV充电效率。S642：应用分段线性化方法，处理GT容量成本模型及约束条件。具体是将GT在t时刻燃料费用与出力函数分成nGT段，每一段引入一个状态变量τi,GT,t和一个连续变量pi,GT,t，其燃料费用的线性化表达式如下：相关约束条件为：pi,GTτi,t≤pi,GT,t≤pi+1,GTτi,GT,t72；对于GT启停成本模型，可以重写为线性形式：S643：应用分段线性化方法，处理PSDR成本容量模型及相关约束条件。具体是将PSDR在t时刻增加负荷费用与出力函数分成段，每一段引入一个状态变量和一个连续变量其增加负荷费用的线性化表达式如下：相关约束条件为：将PSDR在t时刻削减负荷费用与出力函数分成段，每一段引入一个状态变量和一个连续变量其削减负荷费用的线性化表达式如下：。相关约束条件为：S65：根据成本容量模型及相关约束，结合预测信息，构建虚拟电厂A参加日前电力市场优化投标模型。虚拟电厂A出清电量包括弃风电gA,WT,t、GT出力gA,GT,t、需求响应增加负荷需求响应减少负荷EV充电功率和EV放电功率市场所有需求侧资源用集合I表示。虚拟电厂A在时段t总出清电量可表示为：虚拟电厂B总出清电量可表示为：虚拟电厂C总出清电量表示为：构建的虚拟电厂A参加日前电力市场优化投标模型为：πgi,t≤0,i∈I；1020≤gi,t-zki,t≤G11-xki,t,i∈A；1040≤zki,t≤G1xki,t,i∈A；1050≤wki,t-πgi,t≤G21-yki,t,i∈A；1060≤-wki,t≤G2yki,t,i∈A；107ui,t,xki,t,yki,t∈{0,1},i∈A；108式中，系数δi,t取+1或-1，即公式93gi,t系数为正数时取+1，系数为负数时取-1。此外，虚拟电厂A出清电量满足GT有关约束70-78，PSDR相关约束56-60、80-84、86-90，EV有关约束62-68。S7：采用YamipGurobi软件计算虚拟电厂参与日前电力市场的最大收益投标模型，获取投标电量和投标电价，具体步骤如下：S71：应用蒙特卡罗方法对虚拟电厂A的每辆电动汽车行驶距离和并网时间进行模拟，生成nEV辆电动汽车单一时刻驶入和驶离的电动汽车数量及行驶距离输入数据，并采用K-means聚类算法抽取NEV类典型电动汽车，得到每类电动汽车具有相近并网时间和行驶距离。S72：根据NEV类典型电动汽车并网时间、行驶距离和数量，累加得到每类EV充、放电功率上限、并网容量和行驶距离，用NEV类典型EV参数更新投标优化模型中涉及EV模型参数。S73：采用YamipGurobi商用软件计算虚拟电厂参加日前市场的优化投标模型，获取虚拟电厂A内部各类需求侧资源投标电量和电价。本发明不局限于以上所述的具体实施方式，以上所述仅为本发明的较佳实施案例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

权利要求：1.一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：包括以下步骤：S1：在各虚拟电厂提交其需求侧资源的投标电量和电价后，以最小成本为优化目标构建调度中心市场出清模型，获取市场出清价格和出清电量；S2：根据步骤S1得到的市场出清价格，结合虚拟电厂内部需求侧资源的变动运行成本和投标电量，获得该虚拟电厂净收益函数；S3：以虚拟电厂最大净收益为优化目标，结合其内部各需求侧资源最大容量约束，建立虚拟电厂的投标优化模型；S4：根据步骤S1获得的市场出清模型的最优性条件，将其作为虚拟电厂投标优化的部分约束条件，从而将步骤S3建立的投标优化模型转化成具有平衡约束的数学规划模型；S5：采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项，即连续变量乘积，形成混合整数规划模型；S6：结合虚拟电厂内部需求侧资源的出力特性，确定虚拟电厂参与日前市场的优化投标模型；S7：采用YamipGurobi软件计算虚拟电厂参与日前电力市场的最大收益投标模型，获取投标电量和投标电价。2.根据权利要求1所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S1中的构建调度中心市场出清模型的具体步骤如下：S11：在各虚拟电厂提交其需求侧资源i的投标电量Pi和投标电价λi后，调度中心进行经济调度获取出清电价πd和出清电量gi，假设出清电量gi≥0，其目标函数如下：式中，I表示所有需求侧资源集合；S12：负荷供需平衡约束：式中，d为负荷；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；S13：出清电量受到投标电量约束如下：gi≤Pi,i∈Iπgi3；式中，πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益。3.根据权利要求1所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S2中的获得该虚拟电厂净收益函数的计算方式为：设该虚拟电厂为虚拟电厂A，虚拟电厂A的总收入减去运行成本即得净收益R，具体如下：式中，ci表示需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。4.根据权利要求1所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S3中建立虚拟电厂的投标优化模型的步骤具体如下：S31：设该细腻电厂为虚拟电厂A，虚拟电厂A投标优化目标是其净收益最大化，具体如下：式中，ci表示虚拟电厂A的需求侧资源i的变动运行成本；i∈A表示属于虚拟电厂A的需求侧资源集合；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；gi为出清电量。S32：虚拟电厂A的需求侧资源i投标电量具有出力约束，具体如下：Pi≤Pimaxi∈A6；式中，Pi为虚拟电厂A的需求侧资源i的投标电量，Pimax表示虚拟电厂A的需求侧资源i的最大容量。5.根据权利要求4所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S4中的将步骤S3建立的投标优化模型转化成平衡约束的数学规划模型的步骤具体如下：S41：根据对偶理论，获取市场出清模型的最优性条件，具体如式7-16，将其加入到步骤S3建立的虚拟电厂A的投标优化模型的约束条件中，得到虚拟电厂A具有平衡约束的数学规划模型；gi-Pi≤0,i∈A8；πd+πgi-λi≤0,i∈A10；πgi≤0,i∈I12；πd+πgi-λigi＝0,i∈A13；Pi-giπgi＝0,i∈A15；其中，式7-9是步骤S1的构建的调度中心市场出清模型的约束条件，式10-12是对偶问题的约束条件，式13-16是互补条件；式中，为虚拟电厂A的预测市场负荷；表示虚拟电厂A预测其余需求侧资源的投标电量-电价；Pi为虚拟电厂A的需求侧资源i的投标电量，gi为出清电量；λi为投标电价；πd是对偶变量，表示现货价格，即增加单位负荷的边际成本；πgi是对偶变量，表示需求侧资源i增加单位投标电量的边际效益；表示市场其余需求侧资源集合；S42：根据gi≥0,i∈I，对公式10和11进行变形，得到：πd+πgi-λigi≤0,i∈A17；S43：根据公式12，对公式8和9进行变形，得到：gi-Piπgi≤0,i∈A19；S44：根据公式17-20的非正性，得到互补条件13-16的等效公式：S45：使用公式7:重新安排公式21各项，得到等效公式22，即原对偶等式条件为：S46：至此，生成步骤S3建立的虚拟电厂A的投标优化模型的目标函数及其约束条件，步骤S3建立的虚拟电厂A的投标优化模型的目标函数如式5所示，其约束条件为式6-12和式22。6.根据权利要求5所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S5中采用二进制扩充方法处理步骤S4中的具有平衡约束的数学规划模型中的非线性项的具体步骤如下：S51：应用式13和15的互补条件，用λigi和πgiPi重写目标函数5中的非线性项πdgi：πdgi＝λigi-πgiPi23；S52：采用二进制扩充方法处理约束条件22中的非线性项λigi：首先，通过一组离散值{λim,m＝0,1,...,M1}近似连续决策变量λi，其近似表达式为：这里，λi的取值区间为xki为二进制变量；然后，将gi乘以等式24两侧，得到新变量zki＝xkigi，从而得到：同时，引入相应约束：0≤gi-zki≤G11-xki26；0≤zki≤G1xki27；式中，G1是一个常数，通常S53：采用二进制扩充方法，处理约束条件22中的非线性项πgiPi：首先，通过一组离散值{Pim,m＝0,1,...,M2}近似连续决策变量Pi，其近似表达式为这里，Pi的取值区间为Pimin,Pimax，ΔPi＝Pimax-PiminM2，yki为二进制变量；然后，将πgi乘以等式28两侧，得到新变量wki＝πgiyki，从而得到：同时，引入相应约束0≤wki-πgi≤G21-yki30；0≤-wki≤G2yki31；式中，G2是一个常数，通常S54：在目标函数5，约束条件6-12和22应用上述步骤S51-S53的线性转化形式，得到混合整数规划模型：πgi≤0,i∈I40；0≤gi-zki≤G11-xki,i∈A42；0≤zki≤G1xki,i∈A43；0≤wki-πgi≤G21-yki,i∈A44；0≤-wki≤G2yki,i∈A45；xki,yki∈{0,1},i∈A46。7.根据权利要求1所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述步骤S6中确定虚拟电厂参与日前市场的优化投标模型的具体步骤如下：S61：确定虚拟电厂A所参与日前市场的时间间隔为Δt，时段数S62：根据虚拟电厂A和其他虚拟电厂内部需求侧资源组成情况，确定市场中所有需求侧资源种类、数目以及所属虚拟电厂；S63：根据历史信息和市场信息，虚拟电厂A预测得到市场负荷和其余需求侧资源投标电量和电价；S64：根据虚拟电厂A内部需求侧资源特性，构建其在t时段容量成本模型CA,t；S65：根据成本容量模型及相关约束，结合预测信息，构建虚拟电厂A参加日前电力市场优化投标模型。8.根据权利要求7所述的一种虚拟电厂参与日前电力市场的优化投标策略，其特征在于：所述Δt设置为1h，则时段数＝241＝24时段。

百度查询： 广西电网有限责任公司电力科学研究院 一种虚拟电厂参与日前电力市场的优化投标策略

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD