CN108053057A - 一种基于CVaR的虚拟电厂优化调度建模方法 - Google Patents
一种基于CVaR的虚拟电厂优化调度建模方法 Download PDFInfo
- Publication number
- CN108053057A CN108053057A CN201711222793.1A CN201711222793A CN108053057A CN 108053057 A CN108053057 A CN 108053057A CN 201711222793 A CN201711222793 A CN 201711222793A CN 108053057 A CN108053057 A CN 108053057A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- vpp
- risk
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001172 regenerating effect Effects 0.000 claims abstract description 17
- 238000005553 drilling Methods 0.000 claims abstract description 10
- 230000005611 electricity Effects 0.000 claims description 35
- 238000004146 energy storage Methods 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000007726 management method Methods 0.000 claims description 12
- 238000010248 power generation Methods 0.000 claims description 8
- 239000002283 diesel fuel Substances 0.000 claims description 7
- 238000005094 computer simulation Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000002253 acid Substances 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- 239000000446 fuel Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009194 climbing Effects 0.000 claims description 2
- 238000005315 distribution function Methods 0.000 claims description 2
- 238000005286 illumination Methods 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims 1
- 240000002853 Nelumbo nucifera Species 0.000 claims 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 229920000620 organic polymer Polymers 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Secondary Cells (AREA)
Abstract
本发明提出了一种基于条件风险价值CVaR的VPP优化调度建模方法,用风险的理论去探讨VPP调度优化问题,VPP中含有大量风光等可再生能源,其出力与常规能源不同,具有很强的间隙性及波动性,因此VPP优化调度问题是一个具有不确定性的问题,从而使得调度VPP有可能面临风险,而利用CVaR理论可以准确地度量VPP在调度运行中的风险,从而权衡经济性与风险性。本发明在传统VPP调度运行的目标函数中,加入了度量风险的CVaR项,使得该问题数学模型的目标函数不仅考虑总成本最小,还要考虑风险最小,将CVaR乘上一个权重系数加入目标函数,从而将多目标问题转化为单目标问题求解。本发明还对不同风险偏好的运营商进行分级,为不同风险偏好的VPP运营商提供具有指导性的方案。
Description
技术领域
本发明属于电力系统运行与调度技术领域,特别涉及一种基于CVaR的虚拟电厂优化调 度建模方法。
背景技术
传统发电带来的环境污染以及当今世界可持续发展的理念促使包括风力和光伏发电在内 的可再生能源迅速发展。但是,可再生能源发电具有很强的随机性和不可控性,风电和光伏 的单独并网会影响到电力系统的安全运行。为了克服上述缺点,国内外学者提出了虚拟发电 厂(virtual power plant,VPP)[3-4]的概念。
VPP通过对各类分布式发电机组(Distributed Generator,DG)、可控负荷(dispatchable load, DL)和分布式储能设施(Ditributed Energy Storage,DES)等设备的有机聚合,不但可以减少可再 生能源发电间歇性和随机性的影响,并且可以获得规模经济的效益。VPP的特点符合世界的 发展需求,有着广阔的应用前景,截止到2014年底,全球VPP总容量增至4.8GW,预计到 2023年,其规模将增至28GW,年收入可达53亿美元[6]。
然而在虚拟电厂调度运行过程中,由于各分布式电源发电具有间隙性和随机性,导致VPP 运营商的收益具有波动性,如何合理地分配各分布式电源发电量,不仅影响到收益和成本等 经济问题,对于资源能否充分利用、环境问题等也具有深刻的影响。本发明从降低运行商的 经济风险角度出发,用投资组合理论中的条件风险价值为衡量指标,以此做出合理的日前调 度策略。
发明内容
本发明本发明提供了一种基于CVaR的虚拟电厂优化调度模型,能够更好的为决策者提 供最优策略方案。
本发明具体为一种基于CVaR的虚拟电厂优化调度建模方法,所述方法包括以下步骤:
步骤1:确定VPP组成部分与各部分功能,并对各个部分进行建模及分析;
步骤2:为仿真实际运行下VPP各分布式电源出力不确定性,以场景技术通过综合考虑 各场景下电源出力来模拟不确定性;
步骤3:综合考虑各组成部分的实际运行约束,确立以系统总运行成本最优为目标的VPP 优化调度模型目标函数;
步骤4:在上述模型的目标函数中加入CVaR条件风险价值作为理论基础,以此衡量运行 调度成本的经济性与风险性;
步骤5:对风险系数进行分级,以此来表示不同风险偏好的VPP运营商的要求;
步骤6:综合以上步骤,建立考虑CVaR的VPP优化调度模型。
所述步骤1具体包括以下步骤:
(1)光伏发电系统建模分析:
光伏发电量视为服从Beta分布,具有强随机性,假设给定一组光伏电池阵列,其面积和 光电转换效率分别为A和η,则这组光伏电池阵列的输出功率为:
Ppv=rAη
式中r为这一时间段内的实际光照强度;
(2)风力发电系统建模分析:
风力发电量视为服从威布尔分布,Pr为风力发电机的额定功率,vci,vr,vco分别为该风力 发电机的切入风速、额定风速和切出风速,则该风机的输出功率为:
(3)柴油发电机建模分析:柴油发电机的作用是平抑VPP中可再生能源出力的不确定 性,使得各组成部分能够形成一个整体向外作为一个电厂的角色售电;
(4)储能系统建模分析:在VPP中配置储能系统的作用也是为了平抑可再生能源出力 的波动性,储能系统既能对外放电也能向内充电,其存在能够确保VPP每日负荷的供应;
(5)电力负荷建模分析:VPP作为一个电厂自然要承担电厂应有的电力负荷供应需求。
所述步骤2具体为:
由于可再生能源出力的波动不确定性及市场电价的波动性,导致VPP调度运行收益也存 在不确定性,使用场景集的方法处理不确定性问题,在多个场景下使得将不确定性问题转化 为确定性问题处理,分别选取风电出力、光伏出力及电价场景集w、s、p个,各场景概率 分别为π(w)、π(s)、π(p),则VPP在实时运行中的变量均为在各场景下的形式。
所述步骤3具体为:
基于CVaR的VPP优化调度模型的目标函数包括日调度运行总成本及CVaR风险项,其 中日调度总成本包括柴油机燃料成本、可中断负荷成本、VPP运行售电收益;
1)柴油机燃料成本
式中:表示单台柴油发电机年燃料成本,通过下式计算得到:
式中:pf为柴油的单价;为t时刻柴油机出力;Q为柴油机的燃油消耗率;
2)可中断负荷补偿成本
式中:λcurt表示VPP中断用户负荷时给予的中断补偿费用;表示t时刻VPP中断负荷 功率;
3)VPP运行收益
VPP的运行收益来自于在日前市场的收益:
式中:分别为t时刻向电网售电量和购电量;λp(t)表示t时刻在场景p下日 前能量市场统一出清价格MCP;kpurchase表示购电价格在原电价基础上乘以的比例系数。
所述步骤4具体包括以下步骤:
(1)采用历史模拟法来计算风险价值,VaR反映投资组合在给定的置信度β上潜在的最 大损失,设f(x,y)为损失函数,x为决策变量,y为随机变量,假设ρ(y)是y的概率密度函 数,则损失函数f(x,y)不大于边界值α的分布函数为:
对于给定的置信度β,VaR可由下式得到:
VaRβ=min{α∈R:ψ(x,α)≥β}
(2)由于VaR只是某个置信度下的分位点,而分位点以后的风险信息未被考虑,即存 在“尾部风险”,因此采用CVaR风险度量方法,其作为VaR的一种补充风险度量能够更好地 反映“尾部风险”:
由于VaRβ(x)解析式难以求出,所以推导出变换函数Fβ(x,α):
式中:[f(x,y)-α]+表示max{f(x,y)-α,0},α即为VaR的值,
用以下估计式来计算变换函数Fβ(x,α):
式中:y1,...,yq为y的q个样本,易得
所述步骤5具体包括以下步骤:
基于CVaR的VPP调度优化有两个优化目标,一是使调度运行的总成本最小化,二是使 运行VPP不确定因素导致的风险最小化,使用设置权重系数的方法将双目标函数转化为单目 标函数处理,具体方法是在CVaR项上乘以一个权重因子,称之为风险系数L,以此衡量经 济性与风险性,最终目标函数形式如下:
该目标函数分为两个部分,第一部分为VPP调度运行发电成本,第二部分为风险度量项; 定义权重系数L为风险偏好系数,表示运营商对于风险的态度,其取值范围为L≥0,当L小 于0.1时,表示运营商为风险偏好型,当L大于0.5时,表示运营商为风险规避型;VPP运营 商通过对风险系数的设定来衡量经济与风险,按其对于风险的承受程度来选择合适的风险系 数。
所述步骤6具体包括以下步骤:
VPP中各组成部分在运行时具有如下约束:
1)常规机组出力约束
式中:Pc.max为柴油机输出额定功率;ramp为柴油机的爬坡率;为t时段各场景下柴油 机组的总出力,为决策变量;
2)蓄电池电量及充放电约束
式中:分别为t时段各场景下储电量和充、放电量,为决策变量;Wb.max为单组铅酸蓄电池的额定容量;CDOD为蓄电池的最大放电深度;Pbc.max、Pbd.max分别表示单组 铅酸蓄电池的最大充、放电功率;
3)可中断负荷约束
式中:表示t时段各场景下中断负荷量,为决策变量;kcurt为中断系数;Pel(t)表示t时 段下的电负荷;
4)功率平衡约束
式中:分别表示t时刻在场景w和场景s下单台风机和单组光伏阵列出力;
5)CVaR风险约束
为便于求解引入虚拟变量令表示超过VaR的损失,损失函数 f(x,y)取收益的负值;为了便于计算将其松弛为下面两个不等式;
则
式中:α即为VaR的值;β为投资商设定的置信度。
本发明中优化问题的决策变量包括:
本发明具有如下优点和技术效果:
(1)本发明区别于传统VPP优化调度模型只考虑运行成本最小,创新性地提出了考虑 以CVaR为理论基础的风险项,能够权衡经济性与风险性,同时也给不同风险偏好的VPP运 行商提供定量的依据;
(2)通过对本发明所建数学模型的算例分析,得到了一下结论:不同风险偏好的VPP 运营商调度策略有很大不同,当运行策略保守时柴油机及储能发电以平抑可再生能源出力波 动,降低风险性;当运行策略激进时灵活调度柴油机及储能从而增加售电收益,提高经济性。
附图说明
图1为本发明一种基于CVaR的虚拟电厂优化调度建模方法的流程示意图;
图2为VPP系统调度运行总成本关于风险值CVaR的有效边界示意图;
图3为不同风险系数下的优化结果示意图(从左到右依次为激进型、较激进型、较保守 型和保守型);
图4VPP日出力计划调度图。
具体实施方式
下面结合附图对本发明一种基于CVaR的虚拟电厂优化调度建模方法的具体实施方式做 详细阐述。
如图1所示,本发明的一种基于CVaR的虚拟电厂优化调度建模方法包括以下步骤:
步骤1:确定VPP组成部分与各部分功能,并对各个部分进行建模及分析;
步骤2:为仿真实际运行下VPP各分布式电源出力不确定性,以场景技术通过综合考虑 各场景下电源出力来模拟不确定性;
步骤3:综合考虑各组成部分的实际运行约束,确立以系统总运行成本最优为目标的VPP 优化调度模型目标函数;
步骤4:在上述模型的目标函数中加入CVaR条件风险价值作为理论基础,以此衡量运行 调度成本的经济性与风险性;
步骤5:对风险系数进行分级,以此来表示不同风险偏好的VPP运营商的要求;
步骤6:综合以上步骤,建立考虑CVaR的VPP优化调度模型。
下面以运营一个VPP为例,介绍本发明:
本发明算例分析选择不同地区的风、光等自然资源数据及电价、电力负荷数据进行模拟, 选择时间步长1h,调度运行时间24h。分别选择不同的64种场景进行不确定性模拟。并设置 不同的风险偏好系数以表述不同风险偏好的VPP运营商,得到不同的最优调度结果,根据该 计算结果可以画出运行成本与风险值之间的有效前沿曲线,如图2所示。并且在图中将不同 的风险系数划分为4个部分以表示大致的四种不同类型VPP运营商,分别分别代表运营商为 保守型、较保守型、较激进型和激进型。
研究结果表明:在风险系数较小时,运营商倾向于获得收益,调度策略比较激进,随着 CVaR的减小总成本增加缓慢;在风险系数较大时,运营商极力规避风险,调度策略趋向保守, 总成本迅速增加。
为了给不同风险偏好类型的VPP运营商提供具有指导意义的结果,分别选取4种具有代 表性的风险系数:L=0.05代表运营商为激进型;L=0.1代表运营商为较激进型;L=0.5代 表运营商为较保守型;L=2代表运营商为保守型。分析在这4种不同风险系数下,VPP的调 度运行结果就发生怎样的变化。具体数据见图3。
由图中结果可以看见,当风险系数变化时,主要影响VPP调度策略中的柴油机及储能出 力。当风险系数较小时,柴油机出力不仅用于平抑可再生能源出力的波动和负荷供应,还用 于将剩余的柴发出力在电价合适的时候向市场出售,因此此时VPP调度策略也面临着较大的 风险。当风险系数较大时,为了规避风险,主要依靠从电网购电及中断负荷来维持VPP各部 分相对稳定地运作,因此总购电量、总中断负荷量有所上升。
VPP出力计划分析:
VPP作为一个电厂,其日前调度出力计划必不可少,当投资商确定风险偏好水平后,本 文研究模型能够得到VPP每日出力计划。
选取一个典型日,投资商风险系数选值为1。VPP日出力计划如图4所示:
0~5h:电价较低,可再生能源出力小于负荷需求;VPP优先调度储能进行放电,当储能 放电量不足以弥补功率缺额时,启动柴油机进行发电,由于电价很低,选择购入部分电量及 中断部分负荷以满足负荷需求。
6~10h:电价较低,可再生能源出力略大于负荷需求;VPP调度可再生能源和柴油机剩余 出力以及购入部分电量对储能进行充电,当储能电量充满时,停止柴油机和购电。不在电价 低时向市场售电。
11~20h:电价较高,可再生能源出力大于负荷需求;VPP调度可再生能源剩余出力向市 场出售,并调度储能适当放电向市场出售。
21~24h:电价高,可再生能源出力小于负荷需求;VPP调度储能放电及柴油机发电,将 满足负荷需求之外的剩余电力全部向市场出售。
以上仿真结果验证了本发明所构模型有效性和实用性。
条件风险价值的引入不仅能够处理该问题由于风光等可再生能源出力及电价波动带来的 不确定型=性,也能够给出具体有形的风险数据,从而使得该模型能够为不同风险偏好的VPP 运营商在处理经济性与风险之间权衡问题时提供定量的VPP每日调度运行参考依据,更够给 决策者更大的选择空间。
最后应该说明的是,结合上述实施例仅说明本发明的技术方案而非对其限制。所属领域 的普通技术人员应当理解到,本领域技术人员可以对本发明的具体实施方式进行修改或者等 同替换,但这些修改或变更均在申请待批的权利要求保护范围之中。
Claims (7)
1.一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于,所述方法包括以下步骤:
步骤1:确定VPP组成部分与各部分功能,并对各个部分进行建模及分析;
步骤2:为仿真实际运行下VPP各分布式电源出力不确定性,以场景技术通过综合考虑各场景下电源出力来模拟不确定性;
步骤3:综合考虑各组成部分的实际运行约束,确立以系统总运行成本最优为目标的VPP优化调度模型目标函数;
步骤4:在上述模型的目标函数中加入CVaR条件风险价值作为理论基础,以此衡量运行调度成本的经济性与风险性;
步骤5:对风险系数进行分级,以此来表示不同风险偏好的VPP运营商的要求;
步骤6:综合以上步骤,建立考虑CVaR的VPP优化调度模型。
2.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤1具体包括以下步骤:
(1)光伏发电系统建模分析:
光伏发电量视为服从Beta分布,具有强随机性,假设给定一组光伏电池阵列,其面积和光电转换效率分别为A和η,则这组光伏电池阵列的输出功率为:
Ppv=rAη
式中r为这一时间段内的实际光照强度;
(2)风力发电系统建模分析:
风力发电量视为服从威布尔分布,Pr为风力发电机的额定功率,vci,vr,vco分别为该风力发电机的切入风速、额定风速和切出风速,则该风机的输出功率为:
<mrow>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mi>v</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo><</mo>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
<mo><</mo>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mo>></mo>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
k2=-k1vci;
(3)柴油发电机建模分析:柴油发电机的作用是平抑VPP中可再生能源出力的不确定性,使得各组成部分能够形成一个整体向外作为一个电厂的角色售电;
(4)储能系统建模分析:在VPP中配置储能系统的作用也是为了平抑可再生能源出力的波动性,储能系统既能对外放电也能向内充电,其存在能够确保VPP每日负荷的供应;
(5)电力负荷建模分析:VPP作为一个电厂自然要承担电厂应有的电力负荷供应需求。
3.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤2具体为:
由于可再生能源出力的波动不确定性及市场电价的波动性,导致VPP调度运行收益也存在不确定性,使用场景集的方法处理不确定性问题,在多个场景下使得将不确定性问题转化为确定性问题处理,分别选取风电出力、光伏出力及电价场景集w、s、p个,各场景概率分别为π(w)、π(s)、π(p),则VPP在实时运行中的变量均为在各场景下的形式。
4.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤3具体为:
基于CVaR的VPP优化调度模型的目标函数包括日调度运行总成本及CVaR风险项,其中日调度总成本包括柴油机燃料成本、可中断负荷成本、VPP运行售电收益;
1)柴油机燃料成本
<mrow>
<msubsup>
<mi>C</mi>
<mi>F</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>C</mi>
<mi>f</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
</mrow>
式中:表示单台柴油发电机年燃料成本,通过下式计算得到:
<mrow>
<msubsup>
<mi>C</mi>
<mi>f</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>p</mi>
<mi>f</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>v</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>Q</mi>
</mrow>
式中:pf为柴油的单价;为t时刻柴油机出力;Q为柴油机的燃油消耗率;
2)可中断负荷补偿成本
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:λcurt表示VPP中断用户负荷时给予的中断补偿费用;表示t时刻VPP中断负荷功率;
3)VPP运行收益
VPP的运行收益来自于在日前市场的收益:
<mrow>
<msup>
<mi>R</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>&lambda;</mi>
<mi>p</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>E</mi>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>p</mi>
<mi>u</mi>
<mi>r</mi>
<mi>c</mi>
<mi>h</mi>
<mi>a</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>&lambda;</mi>
<mi>p</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:分别为t时刻向电网售电量和购电量;λp(t)表示t时刻在场景p下日前能量市场统一出清价格MCP;kpurchase表示购电价格在原电价基础上乘以的比例系数。
5.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤4具体包括以下步骤:
(1)采用历史模拟法来计算风险价值,VaR反映投资组合在给定的置信度β上潜在的最大损失,设f(x,y)为损失函数,x为决策变量,y为随机变量,假设ρ(y)是y的概率密度函数,则损失函数f(x,y)不大于边界值α的分布函数为:
<mrow>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Integral;</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>&alpha;</mi>
</mrow>
</munder>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>y</mi>
</mrow>
对于给定的置信度β,VaR可由下式得到:
VaRβ=min{α∈R:ψ(x,α)≥β}
(2)由于VaR只是某个置信度下的分位点,而分位点以后的风险信息未被考虑,即存在“尾部风险”,因此采用CVaR风险度量方法,其作为VaR的一种补充风险度量能够更好地反映“尾部风险”:
<mrow>
<msub>
<mi>CVaR</mi>
<mi>&beta;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
</mfrac>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>VaR</mi>
<mi>&beta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>y</mi>
</mrow>
由于VaRβ(x)解析式难以求出,所以推导出变换函数Fβ(x,α):
<mrow>
<msub>
<mi>F</mi>
<mi>&beta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
</mfrac>
<munder>
<mo>&Integral;</mo>
<mrow>
<mi>y</mi>
<mo>&Element;</mo>
<msup>
<mi>R</mi>
<mi>m</mi>
</msup>
</mrow>
</munder>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>&rsqb;</mo>
</mrow>
<mo>+</mo>
</msup>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>y</mi>
</mrow>
式中:[f(x,y)-α]+表示max{f(x,y)-α,0},α即为VaR的值,
用以下估计式来计算变换函数Fβ(x,α):
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>~</mo>
</mover>
<mi>&beta;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>q</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>f</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>&rsqb;</mo>
</mrow>
<mo>+</mo>
</msup>
</mrow>
式中:y1,...,yq为y的q个样本,易得
6.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤5具体包括以下步骤:
基于CVaR的VPP调度优化有两个优化目标,一是使调度运行的总成本最小化,二是使运行VPP不确定因素导致的风险最小化,使用设置权重系数的方法将双目标函数转化为单目标函数处理,具体方法是在CVaR项上乘以一个权重因子,称之为风险系数L,以此衡量经济性与风险性,最终目标函数形式如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>n</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>{</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>C</mi>
<mi>F</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msup>
<mi>R</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>+</mo>
<mi>L</mi>
<mi>C</mi>
<mi>V</mi>
<mi>a</mi>
<mi>R</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
该目标函数分为两个部分,第一部分为VPP调度运行发电成本,第二部分为风险度量项;定义权重系数L为风险偏好系数,表示运营商对于风险的态度,其取值范围为L≥0,当L小于0.1时,表示运营商为风险偏好型,当L大于0.5时,表示运营商为风险规避型;VPP运营商通过对风险系数的设定来衡量经济与风险,按其对于风险的承受程度来选择合适的风险系数。
7.根据权利要求1所述的一种基于CVaR的虚拟电厂优化调度建模方法,其特征在于:所述步骤6具体包括以下步骤:
VPP中各组成部分在运行时具有如下约束:
1)常规机组出力约束
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>v</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>N</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mo>.</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mi>p</mi>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>v</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>v</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>r</mi>
<mi>a</mi>
<mi>m</mi>
<mi>p</mi>
</mrow>
式中:Pc.max为柴油机输出额定功率;ramp为柴油机的爬坡率;为t时段各场景下柴油机组的总出力,为决策变量;
2)蓄电池电量及充放电约束
<mrow>
<msubsup>
<mi>W</mi>
<mi>b</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>W</mi>
<mi>b</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&eta;</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>&eta;</mi>
</mrow>
(1-CDOD)·Nb·Wb.max≤Wb wsp(t)≤Nb·Wb.max
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>N</mi>
<mi>b</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
<mo>.</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>N</mi>
<mi>b</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mi>d</mi>
<mo>.</mo>
<mi>max</mi>
</mrow>
</msub>
</mrow>
式中:分别为t时段各场景下储电量和充、放电量,为决策变量;Wb.max为单组铅酸蓄电池的额定容量;CDOD为蓄电池的最大放电深度;Pbc.max、Pbd.max分别表示单组铅酸蓄电池的最大充、放电功率;
3)可中断负荷约束
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:表示t时段各场景下中断负荷量,为决策变量;kcurt为中断系数;Pel(t)表示t时段下的电负荷;
4)功率平衡约束
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
<mi>w</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>v</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:分别表示t时刻在场景w和场景s下单台风机和单组光伏阵列出力;
5)CVaR风险约束
为便于求解引入虚拟变量令表示超过VaR的损失,损失函数f(x,y)取收益的负值;为了便于计算将其松弛为下面两个不等式;
<mrow>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
<mrow>
<msubsup>
<mi>z</mi>
<mi>k</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mo>{</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msup>
<mi>&lambda;</mi>
<mi>p</mi>
</msup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msubsup>
<mi>E</mi>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&CenterDot;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>p</mi>
<mi>u</mi>
<mi>r</mi>
<mi>c</mi>
<mi>h</mi>
<mi>a</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>&lambda;</mi>
<mi>p</mi>
</msup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mi>F</mi>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>s</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>}</mo>
<mo>-</mo>
<mi>&alpha;</mi>
</mrow>
则
式中:α即为VaR的值;β为投资商设定的置信度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711222793.1A CN108053057B (zh) | 2017-11-29 | 2017-11-29 | 一种基于CVaR的虚拟电厂优化调度建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711222793.1A CN108053057B (zh) | 2017-11-29 | 2017-11-29 | 一种基于CVaR的虚拟电厂优化调度建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108053057A true CN108053057A (zh) | 2018-05-18 |
CN108053057B CN108053057B (zh) | 2021-12-14 |
Family
ID=62120903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711222793.1A Active CN108053057B (zh) | 2017-11-29 | 2017-11-29 | 一种基于CVaR的虚拟电厂优化调度建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108053057B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108960510A (zh) * | 2018-07-04 | 2018-12-07 | 四川大学 | 一种基于两阶段随机规划的虚拟电厂优化交易策略模型 |
CN109066769A (zh) * | 2018-07-20 | 2018-12-21 | 国网四川省电力公司经济技术研究院 | 风电全消纳下虚拟电厂内部资源调度控制方法 |
CN109638886A (zh) * | 2019-01-10 | 2019-04-16 | 国网江苏省电力有限公司宿迁供电分公司 | 一种基于CVaR的风电日前出力申报优化方法 |
CN110075533A (zh) * | 2019-06-05 | 2019-08-02 | 国网北京市电力公司 | 模拟虚拟电厂的实现方法及装置 |
CN110808615A (zh) * | 2019-12-07 | 2020-02-18 | 国家电网有限公司 | 考虑不确定性的气电虚拟电厂调度优化方法 |
CN111402015A (zh) * | 2020-03-17 | 2020-07-10 | 南京工业大学 | 一种基于购售风险的虚拟电厂双层竞标方法及系统 |
CN111582751A (zh) * | 2020-05-19 | 2020-08-25 | 国网吉林省电力有限公司 | 一种时间加权的购电风险预警方法 |
CN112994087A (zh) * | 2021-02-08 | 2021-06-18 | 四川大学 | 基于条件风险约束的多源电力系统中期优化调度方法 |
CN113435714A (zh) * | 2021-06-10 | 2021-09-24 | 华中科技大学 | 一种基于WCVaR风险分析的电网调度优化方法和系统 |
-
2017
- 2017-11-29 CN CN201711222793.1A patent/CN108053057B/zh active Active
Non-Patent Citations (2)
Title |
---|
王海冰: ""考虑条件风险价值的两阶段发电调度随机规划模型和方法"", 《中国电机工程学报》 * |
谭忠富: ""CVAR风险规避方法在具有不确定性和需求响应的风光储能系统连接虚拟电厂动态调度优化模型中的应用"", 《SCIENCEDIRECT》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108960510B (zh) * | 2018-07-04 | 2021-11-26 | 四川大学 | 一种基于两阶段随机规划的虚拟电厂优化交易策略装置 |
CN108960510A (zh) * | 2018-07-04 | 2018-12-07 | 四川大学 | 一种基于两阶段随机规划的虚拟电厂优化交易策略模型 |
CN109066769A (zh) * | 2018-07-20 | 2018-12-21 | 国网四川省电力公司经济技术研究院 | 风电全消纳下虚拟电厂内部资源调度控制方法 |
CN109066769B (zh) * | 2018-07-20 | 2020-03-27 | 国网四川省电力公司经济技术研究院 | 风电全消纳下虚拟电厂内部资源调度控制方法 |
CN109638886A (zh) * | 2019-01-10 | 2019-04-16 | 国网江苏省电力有限公司宿迁供电分公司 | 一种基于CVaR的风电日前出力申报优化方法 |
CN110075533A (zh) * | 2019-06-05 | 2019-08-02 | 国网北京市电力公司 | 模拟虚拟电厂的实现方法及装置 |
CN110808615A (zh) * | 2019-12-07 | 2020-02-18 | 国家电网有限公司 | 考虑不确定性的气电虚拟电厂调度优化方法 |
CN111402015A (zh) * | 2020-03-17 | 2020-07-10 | 南京工业大学 | 一种基于购售风险的虚拟电厂双层竞标方法及系统 |
CN111582751A (zh) * | 2020-05-19 | 2020-08-25 | 国网吉林省电力有限公司 | 一种时间加权的购电风险预警方法 |
CN111582751B (zh) * | 2020-05-19 | 2022-06-14 | 国网吉林省电力有限公司 | 一种时间加权的购电风险预警方法 |
CN112994087A (zh) * | 2021-02-08 | 2021-06-18 | 四川大学 | 基于条件风险约束的多源电力系统中期优化调度方法 |
CN112994087B (zh) * | 2021-02-08 | 2023-06-30 | 四川大学 | 基于条件风险约束的多源电力系统中期优化调度方法 |
CN113435714A (zh) * | 2021-06-10 | 2021-09-24 | 华中科技大学 | 一种基于WCVaR风险分析的电网调度优化方法和系统 |
CN113435714B (zh) * | 2021-06-10 | 2022-08-02 | 华中科技大学 | 一种基于WCVaR风险分析的电网调度优化方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108053057B (zh) | 2021-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108053057A (zh) | 一种基于CVaR的虚拟电厂优化调度建模方法 | |
CN108960510B (zh) | 一种基于两阶段随机规划的虚拟电厂优化交易策略装置 | |
Yousri et al. | Integrated model for optimal energy management and demand response of microgrids considering hybrid hydrogen-battery storage systems | |
Wang et al. | Reliable-economical equilibrium based short-term scheduling towards hybrid hydro-photovoltaic generation systems: Case study from China | |
CN111738497B (zh) | 计及需求侧响应的虚拟电厂双层优化调度方法 | |
CN107464010A (zh) | 一种虚拟电厂容量优化配置方法 | |
Aaslid et al. | Stochastic optimization of microgrid operation with renewable generation and energy storages | |
Hamilton et al. | Dispatch optimization of concentrating solar power with utility-scale photovoltaics | |
Bernal-Agustín et al. | Multi-objective design and control of hybrid systems minimizing costs and unmet load | |
CN104362677A (zh) | 一种主动配电网优化配置结构及其配置方法 | |
CN107565607A (zh) | 一种基于实时电价机制的微电网多时间尺度能量调度方法 | |
CN114069678A (zh) | 一种考虑储能退化成本的光储直流微电网能量调度方法 | |
CN113888204A (zh) | 多主体博弈的虚拟电厂容量优化配置法 | |
Allen et al. | Infrastructure planning and operational scheduling for power generating systems: An energy-water nexus approach | |
Kumar et al. | A New Approach to Design and Optimize Sizing of Hybrid Microgrids in Deregulated Electricity Environment | |
Ali et al. | Development and planning of a hybrid power system based on advance optimization approach | |
Hasan et al. | Electricity cost optimization for large loads through energy storage and renewable energy | |
Meydani et al. | Hybrid Renewable Energy Sources: Their Potential to Meet Electricity Demand | |
TWI639962B (zh) | 一種應用於智慧電網之粒群最佳化模糊邏輯控制充電法 | |
CN116896072A (zh) | 远海风电与制氢系统的优化调控方法和装置 | |
CN116681294A (zh) | 采用改进WCVaR的可调节风险偏好综合能源系统容量控制方法 | |
Azaroual et al. | Model predictive control-based energy management strategy for grid-connected residential photovoltaic–wind–battery system | |
CN115986812A (zh) | 一种考虑储能和需求响应的微电网经济规划方法及装置 | |
CN113255957A (zh) | 综合服务站不确定因素的定量优化分析方法及系统 | |
CN114037331B (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 |