CN109193668A - 一种基于分布鲁棒优化的合同电量分解方法 - Google Patents
一种基于分布鲁棒优化的合同电量分解方法 Download PDFInfo
- Publication number
- CN109193668A CN109193668A CN201811282374.1A CN201811282374A CN109193668A CN 109193668 A CN109193668 A CN 109193668A CN 201811282374 A CN201811282374 A CN 201811282374A CN 109193668 A CN109193668 A CN 109193668A
- Authority
- CN
- China
- Prior art keywords
- unit
- constraint
- contract
- formula
- probability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005096 rolling process Methods 0.000 title claims abstract description 44
- 230000005611 electricity Effects 0.000 claims abstract description 41
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 23
- 239000000243 solution Substances 0.000 claims description 22
- 238000001803 electron scattering Methods 0.000 claims description 21
- 238000004146 energy storage Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000007599 discharging Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000010977 unit operation Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000009194 climbing Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000002402 nanowire electron scattering Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- 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
- 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
- Y02E70/00—Other energy conversion or management systems reducing GHG emissions
- Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
-
- 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
Abstract
本发明公开了一种基于分布鲁棒优化的合同电量分解方法,包括如下步骤:S1、采集成本信息、月发电量与合同电量的差额信息、机组参数以及功率数据,并根据其建立合同电量分解模型;S2、根据合同电量分解模型,建立考虑风电不确定性的分布鲁棒优化模型;S3、根据分布鲁棒优化模型,采用列与约束生成算法,得到风电不确定场景发生的概率值结果;S4、根据概率值结果,采用分布鲁棒优化模型,得到合同电量分解结果;本发明解决了风电以及负荷的不确定性在合同电量分解中的问题,减少了风电不确定性所带来的波动。
Description
技术领域
本发明属于电力市场领域,具体涉及一种基于分布鲁棒优化的合同电量分解方法。
背景技术
合同电量的分解是电网公司电力交易中心的主要任务之一。例如,省电网在年初制定发电计划时,需要根据月负荷预测曲线、机组检修计划、热电联产等因素将年合同电量划分为月份。同时,每月发电计划也需要分解为一天甚至一段时间。高质量的合同电能分解方案可以大大降低后续日调度的难度,是电网经济调度和节能发电调度的基础,因此具有重要的意义。此外,在发电计划的实际实施中,需要根据一天内各时段负荷的不确定性,对所有设备进行调整,对分解方案进行修正,以使电力合同偏差最小。
月度机组组合作为一种中长期电力的资源优化,可以较长时间考虑电网的运行效率。这种能源结构也决定了机组不应该采用频繁启动和关闭的最佳运行模式,从而提高了月度机组组合在实际生产中的重要性。因此,可以考虑基于月度机组组合的月度发电计划来解决合同功率分解问题。
然而,近年来,随着风力发电等分布式发电系统的接入,电力系统变得越来越复杂。合同电量的分解可能会面临一个巨大的挑战,以适应不确定的DG输出和负载。对于DG输出的不确定性建模,常用的方法通常包括随机优化(SO)和鲁棒优化(RO)。虽然这两种不确定性方法在一定程度上达到了一定的处理效果,但它们也有各自的缺陷。分布式鲁棒优化(DRO)作为一种很有前途的不确定性建模方法,近年来得到了广泛的研究。它把随机优化和鲁棒优化结合起来,在给定的置信集范围内寻找最坏概率分布。然而,DRO方法通常通过利用Wasserstein度量或不确定的一阶/二阶矩信息来构造概率分布集,因此通过使用对偶性和刻画矩信息,求解过程相当复杂。
发明内容
针对现有技术中的上述不足,本发明提供的一种基于分布鲁棒优化的合同电量分解方法,用于解决风电以及负荷的不确定性在合同电量分解中的问题,减少风电不确定性所带来的波动。
为了达到上述发明目的,本发明采用的技术方案为:
一种基于分布鲁棒优化的合同电量分解方法,包括如下步骤:
S1:采集成本信息、月发电量与合同电量的差额信息、机组参数以及功率数据,并根据其建立合同电量分解模型;
S2:根据合同电量分解模型,建立考虑风电不确定性的分布鲁棒优化模型;
S3:根据分布鲁棒优化模型,采用列与约束生成算法,得到风电不确定场景发生的概率值结果;
S4:根据概率值结果,采用分布鲁棒优化模型,得到合同电量分解结果。
进一步地,步骤S1中,建立合同电量分解模型,包括如下步骤:
S1-1:根据成本信息和月发电量与合同电量的差额信息,构建合同电量分解模型的目标函数;
S1-2:根据机组参数和功率数据,制定合同电量分解模型的约束条件。
进一步地,步骤S1-1中,目标函数的公式为:
F=min(F1+F2+F3)
式中,F为总体的优化目标函数;F1为机组起停成本;F2为机组月发电量与合同电量的差额;F3为弃风成本。
进一步地,机组起停成本的计算公式为:
式中,F1为机组的起停成本;T为总天数;td为一天的总时间段数;NG为总机组数;STi、SDi分别为i机组的开始和停止成本;ui,t、ui,t-1分别为i机组在t和t-1时间段的启动和停止状态值,停止为0,启动为1;t为时间段变量;i为机组变量;
机组月发电量与合同电量的差额的计算公式为:
式中,F2为机组月发电量与合同电量的差额;Li每月第i机组功率调度偏差的增量;
弃风成本的计算公式为:
式中,F3为弃风成本;NW为风力机数量;δ为弃风惩罚系数;为第一台风力机在t时间段内的预测和实际调度输出。
进一步地,步骤S1-2中,约束条件包括电力平衡约束、机组运行约束、ESS约束以及潮流约束;
电力平衡约束的公式为:
式中,Pi,t为i机组在t时间段内的有功输出;分别为i机组在t时间段内的充电和放电有功注入;为t时间段内的负载总有功功率;ESS为储能系统的机组节点集合;NESS为总储能系统ESS数量;为i风电机组在t时间段内的有功功率;
机组运行约束的公式为:
式中,Pi,min、Pi,max分别为i机组有功功率输出的下限和上限;ui,t分别为i机组在t时间段内的启动和停止状态值,停止为0,启动为1;yi,t为i机组在t时间段内的是否有停机到开机状态变化的标志值,没有为0,有为1;zi,t为i机组在t时间段内的从开始到关闭状态的变化的标志值,关闭为0,开始为1;rUi、rDi分别为i机组的爬坡上限和下限;TDi、TUi分别为最小启动时间和最小停机时间;t为时间段变量;i为机组变量;
ESS约束的公式为:
式中,ESS为储能系统的机组节点集合;Ei,t为i机组在t时间段内所储存的电量;分别为充电效率和放电效率;分别为i机组的充电和放电状态;为i机组在t时间段内的充电和放电有功功率值;为i机组在t时间段内的充电和放电有功功率的上限值;SOCj.max、SOCj.min为存储有功功率比的上限和下限;分别为i机组上存储的有功功率的上限和下限;
潮流约束的公式为:
式中,Pt、分别为各常规机组、风电机组、储能装置和负荷需求在t时间段内有功功率在系统总节点维度下的向量表示形式;Pline为支路功率;为支路功率上限和下限;xn为支路n的电抗,n=1,2,...,NL;NL为总支路数;L为系统支路节点的连接矩阵;B为B系数矩阵;Bdiag为B系数对角矩阵;B-1为B系数逆矩阵。
进一步地,步骤S2包括如下步骤:
S2-1:在合同电量分解模型中添加风电不确定性,得到不确定模型;
S2-2:对不确定模型的概率进行线性化处理,得到分布鲁棒优化模型。
进一步地,步骤S2-1中,不确定模型的公式为:
式中,zs为离散型决策变量即第一阶段变量;ys为连续型决策变量即第二阶段变量;us为风电的不确定性出力及负荷的预测值;ψ为概率分布置信区间;ps为s场景下的概率值;a为关于目标函数的约束系数的向量形式表示;A为关于第一阶段变量的约束系数的向量形式表示;b为关于第一阶段变量的约束的常数的向量形式表示;c为二阶锥约束中松弛外变量系数的向量形式表示;C为关于第二阶段变量的约束系数的向量形式表示;d为一、二阶段变量耦合约束中第二阶段变量的系数矩阵;D为一、二阶段变量耦合约束中第二阶段变量的系数的向量形式表示;Q为二阶锥约束中松弛内变量系数的向量形式表示;G为一、二阶段耦合约束中第一阶段变量的系数;E为风电不确定出力约束的系数;Ys为关于第二阶段变量的所有约束的向量形式表示;q为二阶锥约束中松弛内的常数的向量形式表示;g为一、二阶段耦合约束中的常数;T为表示矩阵的转置;s为场景变量;Ns为典型场景个数。
进一步地,概率分布置信区间的公式为:
式中,ψ为概率分布置信区间;ψ1、ψ∞分别为由1范数和无穷范数所限制的概率分布置信区间;p0为初始概率值;p、p0分别为概率和初始概率值的向量表示形式;ps为场景s的概率;θ1、θ∞分别为一范数和无穷范数所限制的概率波动范围;表示场景正数概率。
进一步地,步骤S2-2中,线性化处理的公式为:
式中,为概率变大的标志变量,变大为1;为概率变大量;为概率减小量;为场景s的初始概率值;为概率减小的标志变量,变小为1。
进一步地,步骤S3中,使用列与约束生成算法进行求解计算,包括如下步骤:
S3-1:根据分布鲁棒优化模型,初始化参数,包括上界、下界以及迭代参数;
S3-2:进行第一阶段求解,即求解列与约束生成算法主问题,获得最优决策结果,并更新下界值;
S3-3:固定第一阶段变量的决策值,进行第二阶段求解,即求解列与约束生成算法子问题,获得最优解以及最优目标函数值;
最优解为风电不确定场景发生的概率值;
S3-4:更新上界值,判断上、下界值之差是否小于等于预先设定的允许偏差值,若是则停止迭代,返回最优解,否则更新主问题恶劣概率分布,在主问题中更新第二阶段变量并添加新的第二阶段变量约束,更新迭代参数,返回步骤S3-2。
本发明的有益效果为:
(1)本发明提出的一种基于分布鲁棒优化的合同电量分解方法,通过改变网络潮流,储能系统ESS可以大大减少机组月发电量和合同电量的差额,从而减少风电不确定性所带来的波动;
(2)分布式鲁棒优化模型在经济性和稳健性方面取得了较好的均衡效果,充分考虑了场景概率分布的不确定性,获得了最低的期望值,显示最优经济性能和不确定性适应性;
(3)利用数据驱动下的DRO方法解决风电以及负荷的不确定性在合同电量分解中的问题。
附图说明
图1为基于分布鲁棒优化的合同电量分解方法流程图;
图2为建立合同电量分解模型的方法流程图;
图3为建立考虑风电不确定性的分布鲁棒优化模型的方法流程图;
图4为使用列与约束生成算法进行求解计算的方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
一种基于分布鲁棒优化的合同电量分解方法,如图1所示,包括如下步骤:
S1:采集成本信息、月发电量与合同电量的差额信息、机组参数以及功率数据,并根据其建立合同电量分解模型,合同电量分解模型的优化目标是在保证安全约束的前提下,通过调节ESS和WT,使机组月发电量与合同电量之差最小,同时,为了减少启停机组的数量,在优化目标中增加了机组的启动成本,保证了机组的连续运行,步骤S1中,建立合同电量分解模型,如图2所示,包括如下步骤:
S1-1:根据成本信息和月发电量与合同电量的差额信息,构建合同电量分解模型的目标函数,目标函数的公式为:
F=min(F1+F2+F3)
式中,F为总体的优化目标函数;F1为机组起停成本;F2为机组月发电量与合同电量的差额;F3为弃风成本;
机组起停成本的计算公式为:
式中,F1为机组的起停成本;T为总天数;td为一天的总时间段数;NG为总机组数;STi、SDi分别为i机组的开始和停止成本;ui,t、ui,t-1分别为i机组在t和t-1时间段的启动和停止状态值,停止为0,启动为1;t为时间段变量;i为机组变量;
机组月发电量与合同电量的差额的计算公式为:
式中,F2为机组月发电量与合同电量的差额;Li每月第i机组功率调度偏差的增量;
机组功率调度偏差的增量的公式为:
式中,为i机组在tc天的实际发电量;E0,i为i机组的合同月电量;tc为天数变量;
弃风成本的计算公式为:
式中,F3为弃风成本;NW为风力机数量;δ为弃风惩罚系数;为第一台风力机在t时间段内的预测和实际调度输出;
S1-2:根据机组参数和功率数据,制定合同电量分解模型的约束条件;
约束条件包括电力平衡约束、机组运行约束、ESS约束以及潮流约束;
电力平衡约束的公式为:
式中,Pi,t为i机组在t时间段内的有功输出;分别为i机组在t时间段内的充电和放电有功注入;为t时间段内的负载总有功功率;ESS为储能系统的机组节点集合;NESS为总储能系统ESS数量;为i风电机组在t时间段内的有功功率;
机组运行约束的公式为:
式中,Pi,min、Pi,max分别为i机组有功功率输出的下限和上限;ui,t分别为i机组在t时间段内的启动和停止状态值,停止为0,启动为1;yi,t为i机组在t时间段内的是否有停机到开机状态变化的标志值,没有为0,有为1;zi,t为i机组在t时间段内的从开始到关闭状态的变化的标志值,关闭为0,开始为1;rUi、rDi分别为i机组的爬坡上限和下限;TDi、TUi分别为最小启动时间和最小停机时间;t为时间段变量;i为机组变量;
ESS约束的公式为:
式中,ESS为储能系统的机组节点集合;Ei,t为i机组在t时间段内所储存的电量;分别为充电效率和放电效率;分别为i机组的充电和放电状态;为i机组在t时间段内的充电和放电有功功率值;为i机组在t时间段内的充电和放电有功功率的上限值;SOCj.max、SOCj.min为存储有功功率比的上限和下限;分别为i机组上存储的有功功率的上限和下限;
潮流约束的公式为:
式中,Pt、分别为各常规机组、风电机组、储能装置和负荷需求在t时间段内有功功率在系统总节点维度下的向量表示形式;Pline为支路功率;为支路功率上限和下限;xn为支路n的电抗,n=1,2,...,NL;NL为总支路数;L为系统支路节点的连接矩阵;B为B系数矩阵;Bdiag为B系数对角矩阵;B-1为B系数逆矩阵;
S2:根据合同电量分解模型,建立考虑风电不确定性的分布鲁棒优化模型,采用数据驱动的分布式鲁棒方法来分析风电和负荷的不确定性,根据各设备调节的灵活性,对启停状态设计第一相位变量,将常规机组的蓄电相关变量设计成第一阶段变量,而其他变量设计成第二阶段变量,可以根据实际风机出力和负荷波动灵活调整,假设K个实际场景是通过历史数据获得的,并且选择有限的离散场景,在每个场景中的概率分布是P0,很明显,现实中的概率分布不一定是P0,从而得到基于数据驱动的分布式鲁棒无功优化模型,如图3所示,包括如下步骤:
S2-1:在合同电量分解模型中添加风电不确定性,得到不确定模型;
不确定模型的公式为:
式中,zs为离散型决策变量即第一阶段变量;ys为连续型决策变量即第二阶段变量;us为风电的不确定性出力及负荷的预测值;ψ为概率分布置信区间;ps为s场景下的概率值;a为关于目标函数的约束系数的向量形式表示;A为关于第一阶段变量的约束系数的向量形式表示;b为关于第一阶段变量的约束的常数的向量形式表示;c为二阶锥约束中松弛外变量系数的向量形式表示;C为关于第二阶段变量的约束系数的向量形式表示;d为一、二阶段变量耦合约束中第二阶段变量的系数矩阵;D为一、二阶段变量耦合约束中第二阶段变量的系数的向量形式表示;Q为二阶锥约束中松弛内变量系数的向量形式表示;G为一、二阶段耦合约束中第一阶段变量的系数;E为风电不确定出力约束的系数;Ys为关于第二阶段变量的所有约束的向量形式表示;q为二阶锥约束中松弛内的常数的向量形式表示;g为一、二阶段耦合约束中的常数;T为表示矩阵的转置;s为场景变量;Ns为典型场景个数;
由于历史数据的信息有限,不确定性的概率分布无法精确确定。因此,我们允许不确定性的概率分布在由历史数据构造的预定义置信集内是任意的,概率分布置信区间的公式为:
式中,ψ为概率分布置信区间;ψ1、ψ∞分别为由1范数和无穷范数所限制的概率分布置信区间;p0为初始概率值;p、p0分别为概率和初始概率值的向量表示形式;ps为场景s的概率;θ1、θ∞分别为一范数和无穷范数所限制的概率波动范围;表示场景正数概率;
将概率分布置信区间公式进行变换,得到求解公式为:
式中,Pr{·}为概率的波动范围的求和;K为生成的总场景数目;α1、α∞为1范数和无穷范数对应的置信度;
S2-2:对不确定模型的概率进行线性化处理,得到分布鲁棒优化模型;
线性化处理的公式为:
式中,为概率变大的标志变量,变大为1;为概率变大量;为概率减小量;为场景s的初始概率值;为概率减小的标志变量,变小为1;
S3:根据分布鲁棒优化模型,采用列与约束生成算法,得到风电不确定场景发生的概率值结果,如图4所示,包括如下步骤:
S3-1:根据分布鲁棒优化模型,初始化参数,包括上界、下界以及迭代参数;
设置下界LB=0,上界UB=+∞,迭代次数m=0;
S3-2:进行第一阶段求解,即求解列与约束生成算法主问题,获得最优决策结果(z*,η*),并更新下界值LB=max{LB,η*};
其中,η*为第一阶段优化目标值,z*为第一阶段变量的决策值;
S3-3:固定第一阶段变量的决策值z*,进行第二阶段求解,即求解列与约束生成算法子问题,获得最优解以及最优目标函数值
S3-4:更新上界值判断上、下界值之差是否小于等于预先设定的允许偏差值,即(UB-LB)≤ε,若是则停止迭代,返回最优解即第m次迭代中风电不确定场景发生的概率值,否则更新主问题恶劣概率分布在主问题中更新第二阶段变量ys,m+1并添加新的第二阶段变量约束Ys,m+1,更新迭代参数m=m+1,返回步骤S3-2;
S4:根据概率值结果,采用分布鲁棒优化模型,得到合同电量分解结果。
算例分析:对改进的IEEE 39节点进行了算例分析,以证明所提出的模型,实验用MATLAB R2014a在个人计算机上实现,ESS安装在母线16和33上,它们的充电或放电功率限制均为0.3毫瓦,SOCj.min=0.8,SOCj.max=1.2,常规单元的参数如表1发电机组参数表所示。
表1
此外,假设各母线的负荷和风电预测基础数据相同,利用均值为0、方差为0.25的正态分布,基于基础数据生成10000个情景,将其聚类为10个典型情景。
(1)储能系统的影响
将对比有无储能系统的效果,而其他约束保持不变,如表2有无ESS时机组月发电量与合同电量偏差比较表所示。
表2
如表2所示,机组的月发电量与合同电量之间的差额将大大减少,而储能系统可以改变网络功率流动以减小这种差异。机组的月发电量与一天内的合同电量之差可以显著降低,采用储能系统的优化效果更好。
(2)与随机优化、鲁棒方法和分布式鲁棒方法的比较分析
将分布式鲁棒优化(DRO)与传统随机优化(SO)和鲁棒优化(RO)方法进行比较,随机方法基于分布式鲁棒方法已知的10个场景,每个场景的概率值为0.1,鲁棒优化采用预测值的0.25倍作为波动范围,并采用区域优化试验,随机生成10000个随机概率分布组合,以获得概率分布网络损耗的平均值和最大值,如表3的概率分布表所示,在DRO法中α∞=0.8,α1∈{0.2,0.99}。
表3
从表3中看出,RO通常考虑最坏的风场和负载场景信息,并且相应的优化结果最大,而SO仅考虑10个已知的离散性,并且优化结果最小;DRO模型对经济和保守主义具有较好的均衡效应,并得到了最小概率分布期望,并显示了最优的经济性能和不确定性的适应性。
本发明的有益效果为:
(1)本发明提出的一种基于分布鲁棒优化的合同电量分解方法,通过改变网络潮流,储能系统ESS可以大大减少机组月发电量和合同电量的差额,从而减少风电不确定性所带来的波动;
(2)分布式鲁棒优化模型在经济性和稳健性方面取得了较好的均衡效果,充分考虑了场景概率分布的不确定性,获得了最低的期望值,显示最优经济性能和不确定性适应性;
(3)利用数据驱动下的DRO方法解决风电以及负荷的不确定性在合同电量分解中的问题。
Claims (10)
1.一种基于分布鲁棒优化的合同电量分解方法,其特征在于,包括如下步骤:
S1:采集成本信息、月发电量与合同电量的差额信息、机组参数以及功率数据,并根据其建立合同电量分解模型;
S2:根据合同电量分解模型,建立考虑风电不确定性的分布鲁棒优化模型;
S3:根据分布鲁棒优化模型,采用列与约束生成算法,得到风电不确定场景发生的概率值结果;
S4:根据概率值结果,采用分布鲁棒优化模型,得到合同电量分解结果。
2.根据权利要求1所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S1中,建立合同电量分解模型,包括如下步骤:
S1-1:根据成本信息和月发电量与合同电量的差额信息,构建合同电量分解模型的目标函数;
S1-2:根据机组参数和功率数据,制定合同电量分解模型的约束条件。
3.根据权利要求2所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S1-1中,目标函数的公式为:
F=min(F1+F2+F3)
式中,F为总体的优化目标函数;F1为机组起停成本;F2为机组月发电量与合同电量的差额;F3为弃风成本。
4.根据权利要求3所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,机组起停成本的计算公式为:
式中,F1为机组的起停成本;T为总天数;td为一天的总时间段数;NG为总机组数;STi、SDi分别为i机组的开始和停止成本;ui,t、ui,t-1分别为i机组在t和t-1时间段的启动和停止状态值,停止为0,启动为1;t为时间段变量;i为机组变量;
机组月发电量与合同电量的差额的计算公式为:
式中,F2为机组月发电量与合同电量的差额;Li每月第i机组功率调度偏差的增量;
弃风成本的计算公式为:
式中,F3为弃风成本;NW为风力机数量;δ为弃风惩罚系数;为第一台风力机在t时间段内的预测和实际调度输出。
5.根据权利要求2所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S1-2中,所述约束条件包括电力平衡约束、机组运行约束、ESS约束以及潮流约束;
电力平衡约束的公式为:
式中,Pi,t为i机组在t时间段内的有功输出;分别为i机组在t时间段内的充电和放电有功注入;Pt L为t时间段内的负载总有功功率;ESS为储能系统的机组节点集合;NESS为总储能系统ESS数量;为i风电机组在t时间段内的有功功率;
机组运行约束的公式为:
式中,Pi,min、Pi,max分别为i机组有功功率输出的下限和上限;ui,t分别为i机组在t时间段内的启动和停止状态值,停止为0,启动为1;yi,t为i机组在t时间段内的是否有停机到开机状态变化的标志值,没有为0,有为1;zi,t为i机组在t时间段内的从开始到关闭状态的变化的标志值,关闭为0,开始为1;rUi、rDi分别为i机组的爬坡上限和下限;TDi、TUi分别为最小启动时间和最小停机时间;t为时间段变量;i为机组变量;
ESS约束的公式为:
式中,ESS为储能系统的机组节点集合;Ei,t为i机组在t时间段内所储存的电量;分别为充电效率和放电效率;分别为i机组的充电和放电状态;为i机组在t时间段内的充电和放电有功功率值;为i机组在t时间段内的充电和放电有功功率的上限值;SOCj.max、SOCj.min为存储有功功率比的上限和下限;分别为i机组上存储的有功功率的上限和下限;
潮流约束的公式为:
式中,Pt、Pt W、Pt ESS、Pt L分别为各常规机组、风电机组、储能装置和负荷需求在t时间段内有功功率在系统总节点维度下的向量表示形式;Pline为支路功率;为支路功率上限和下限;xn为支路n的电抗,n=1,2,...,NL;NL为总支路数;L为系统支路节点的连接矩阵;B为B系数矩阵;Bdiag为B系数对角矩阵;B-1为B系数逆矩阵。
6.根据权利要求1所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S2包括如下步骤:
S2-1:在合同电量分解模型中添加风电不确定性,得到不确定模型;
S2-2:对不确定模型的概率进行线性化处理,得到分布鲁棒优化模型。
7.根据权利要求6所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S2-1中,不确定模型的公式为:
式中,zs为离散型决策变量即第一阶段变量;ys为连续型决策变量即第二阶段变量;us为风电的不确定性出力及负荷的预测值;ψ为概率分布置信区间;ps为s场景下的概率值;a为关于目标函数的约束系数的向量形式表示;A为关于第一阶段变量的约束系数的向量形式表示;b为关于第一阶段变量的约束的常数的向量形式表示;c为二阶锥约束中松弛外变量系数的向量形式表示;C为关于第二阶段变量的约束系数的向量形式表示;d为一、二阶段变量耦合约束中第二阶段变量的系数矩阵;D为一、二阶段变量耦合约束中第二阶段变量的系数的向量形式表示;Q为二阶锥约束中松弛内变量系数的向量形式表示;G为一、二阶段耦合约束中第一阶段变量的系数;E为风电不确定出力约束的系数;Ys为关于第二阶段变量的所有约束的向量形式表示;q为二阶锥约束中松弛内的常数的向量形式表示;g为一、二阶段耦合约束中的常数;T为表示矩阵的转置;s为场景变量;Ns为典型场景个数。
8.根据权利要求7所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,概率分布置信区间的公式为:
式中,ψ为概率分布置信区间;ψ1、ψ∞分别为由1范数和无穷范数所限制的概率分布置信区间;p0为初始概率值;p、p0分别为概率和初始概率值的向量表示形式;ps为场景s的概率;θ1、θ∞分别为一范数和无穷范数所限制的概率波动范围;表示场景正数概率。
9.根据权利要求8所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S2-2中,线性化处理的公式为:
式中,为概率变大的标志变量,变大为1;为概率变大量;为概率减小量;为场景s的初始概率值;为概率减小的标志变量,变小为1。
10.根据权利要求1所述的基于分布鲁棒优化的合同电量分解方法,其特征在于,所述步骤S3中,使用列与约束生成算法进行求解计算,包括如下步骤:
S3-1:根据分布鲁棒优化模型,初始化参数,包括上界、下界以及迭代参数;
S3-2:进行第一阶段求解,即求解列与约束生成算法主问题,获得最优决策结果,并更新下界值;
S3-3:固定第一阶段变量的决策值,进行第二阶段求解,即求解列与约束生成算法子问题,获得最优解以及最优目标函数值;
所述最优解为风电不确定场景发生的概率值;
S3-4:更新上界值,判断上、下界值之差是否小于等于预先设定的允许偏差值,若是则停止迭代,返回最优解,否则更新主问题恶劣概率分布,在主问题中更新第二阶段变量并添加新的第二阶段变量约束,更新迭代参数,返回步骤S3-2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811282374.1A CN109193668B (zh) | 2018-10-31 | 2018-10-31 | 一种基于分布鲁棒优化的合同电量分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811282374.1A CN109193668B (zh) | 2018-10-31 | 2018-10-31 | 一种基于分布鲁棒优化的合同电量分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109193668A true CN109193668A (zh) | 2019-01-11 |
CN109193668B CN109193668B (zh) | 2021-07-20 |
Family
ID=64941030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811282374.1A Expired - Fee Related CN109193668B (zh) | 2018-10-31 | 2018-10-31 | 一种基于分布鲁棒优化的合同电量分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109193668B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法 |
CN110543966A (zh) * | 2019-07-23 | 2019-12-06 | 四川大学 | 一种电热气联供的微能源网的日前调度优化方法 |
CN110826763A (zh) * | 2019-09-19 | 2020-02-21 | 国网浙江省电力有限公司台州供电公司 | 一种基于引导学习策略的中长期合约电量分解方法 |
CN111080082A (zh) * | 2019-11-27 | 2020-04-28 | 国网山东省电力公司经济技术研究院 | 一种适应低碳电源发展的电网规划方法 |
CN112053034A (zh) * | 2020-08-11 | 2020-12-08 | 云南电网有限责任公司 | 考虑风电不确定性分布特征的电网可调节鲁棒优化调度方法 |
CN113869927A (zh) * | 2021-03-29 | 2021-12-31 | 四川大学 | 一种促进社区微网内部多产消者能源共享的分时定价方法 |
CN114123316A (zh) * | 2021-11-10 | 2022-03-01 | 云南电网有限责任公司 | 计及风电不确定性的scuc模型优化方法及设备 |
CN115688394A (zh) * | 2022-10-18 | 2023-02-03 | 上海科技大学 | 考虑电网多重不确定性的v2g分布鲁棒优化方法 |
CN117422183A (zh) * | 2023-12-18 | 2024-01-19 | 国网四川省电力公司 | 机组检修优化方法、装置、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017028869A (ja) * | 2015-07-22 | 2017-02-02 | 富士電機株式会社 | 需給計画作成装置、プログラム |
CN107785885A (zh) * | 2016-08-26 | 2018-03-09 | 中国电力科学研究院 | 一种双边电量合同参与日前电力电量平衡的协调优化方法 |
CN108631328A (zh) * | 2018-07-04 | 2018-10-09 | 四川大学 | 一种考虑dg无功支撑和开关重构的主动配电网分布鲁棒无功优化方法 |
-
2018
- 2018-10-31 CN CN201811282374.1A patent/CN109193668B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017028869A (ja) * | 2015-07-22 | 2017-02-02 | 富士電機株式会社 | 需給計画作成装置、プログラム |
CN107785885A (zh) * | 2016-08-26 | 2018-03-09 | 中国电力科学研究院 | 一种双边电量合同参与日前电力电量平衡的协调优化方法 |
CN108631328A (zh) * | 2018-07-04 | 2018-10-09 | 四川大学 | 一种考虑dg无功支撑和开关重构的主动配电网分布鲁棒无功优化方法 |
Non-Patent Citations (2)
Title |
---|
ROBERT MIETH: "Data-Driven Distributionally Robust Optimal Power Flow for Distribution Systems", 《IEEE CONTROL SYSTEMS LETTERS》 * |
税月: "考虑风电不确定性的电气能源系统两阶段分布鲁棒协同调度", 《电力系统自动化》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009152A (zh) * | 2019-04-03 | 2019-07-12 | 东南大学 | 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法 |
CN110543966A (zh) * | 2019-07-23 | 2019-12-06 | 四川大学 | 一种电热气联供的微能源网的日前调度优化方法 |
CN110826763B (zh) * | 2019-09-19 | 2022-07-26 | 国网浙江省电力有限公司台州供电公司 | 一种基于引导学习策略的中长期合约电量分解方法 |
CN110826763A (zh) * | 2019-09-19 | 2020-02-21 | 国网浙江省电力有限公司台州供电公司 | 一种基于引导学习策略的中长期合约电量分解方法 |
CN111080082A (zh) * | 2019-11-27 | 2020-04-28 | 国网山东省电力公司经济技术研究院 | 一种适应低碳电源发展的电网规划方法 |
CN111080082B (zh) * | 2019-11-27 | 2023-08-18 | 国网山东省电力公司经济技术研究院 | 一种适应低碳电源发展的电网规划方法 |
CN112053034B (zh) * | 2020-08-11 | 2022-08-05 | 云南电网有限责任公司 | 一种考虑风电不确定性分布特征的电网可调节鲁棒优化调度方法 |
CN112053034A (zh) * | 2020-08-11 | 2020-12-08 | 云南电网有限责任公司 | 考虑风电不确定性分布特征的电网可调节鲁棒优化调度方法 |
CN113869927A (zh) * | 2021-03-29 | 2021-12-31 | 四川大学 | 一种促进社区微网内部多产消者能源共享的分时定价方法 |
CN113869927B (zh) * | 2021-03-29 | 2023-10-24 | 四川大学 | 一种促进社区微网内部多产消者能源共享的分时定价方法 |
CN114123316A (zh) * | 2021-11-10 | 2022-03-01 | 云南电网有限责任公司 | 计及风电不确定性的scuc模型优化方法及设备 |
CN114123316B (zh) * | 2021-11-10 | 2023-06-13 | 云南电网有限责任公司 | 计及风电不确定性的scuc模型优化方法及设备 |
CN115688394A (zh) * | 2022-10-18 | 2023-02-03 | 上海科技大学 | 考虑电网多重不确定性的v2g分布鲁棒优化方法 |
CN115688394B (zh) * | 2022-10-18 | 2023-12-26 | 上海科技大学 | 考虑电网多重不确定性的v2g分布鲁棒优化方法 |
CN117422183A (zh) * | 2023-12-18 | 2024-01-19 | 国网四川省电力公司 | 机组检修优化方法、装置、设备及存储介质 |
CN117422183B (zh) * | 2023-12-18 | 2024-03-19 | 国网四川省电力公司 | 机组检修优化方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109193668B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109193668A (zh) | 一种基于分布鲁棒优化的合同电量分解方法 | |
Khatod et al. | Evolutionary programming based optimal placement of renewable distributed generators | |
CN111355265B (zh) | 一种微电网能量两阶段鲁棒优化方法及系统 | |
CN109149571A (zh) | 一种考虑系统燃气和火电机组特性的储能优化配置方法 | |
CN108306303A (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
CN109636056A (zh) | 一种基于多智能体技术的多能源微网去中心化优化调度方法 | |
CN107623337B (zh) | 一种微电网能量管理方法 | |
CN107017656A (zh) | 一种考虑电能质量提升的主动配电系统优化调度方法 | |
CN109034587B (zh) | 一种协调多种可控单元的主动配电系统优化调度方法 | |
CN105932722A (zh) | 计及风电出力不确定性的抽水蓄能-风电联合调度方法 | |
CN113705962A (zh) | 一种基于分布鲁棒优化的虚拟电厂日前调度方法 | |
He et al. | Hierarchical optimal energy management strategy of hybrid energy storage considering uncertainty for a 100% clean energy town | |
CN115622146A (zh) | 一种梯级水光蓄互补系统调度决策方法 | |
CN105305485A (zh) | 一种消纳大规模间歇性能源的安全约束经济调度方法 | |
CN109038653A (zh) | 一种追求分散式风电高效消纳的配电网优化运行方法 | |
Chang et al. | Model predictive control based energy collaborative optimization management for energy storage system of virtual power plant | |
CN111523947A (zh) | 一种虚拟电厂发电成本生成方法 | |
CN110717694B (zh) | 基于新能源消纳期望值的储能配置随机决策方法及装置 | |
CN113013909A (zh) | 一种基于平抑牵引的储能容量改进方法 | |
Fan et al. | Research on Optimal Operation and Capacity Configuration of Energy Storage Based on Pumped Storage Station | |
Ma et al. | Consensus-based hierachical demand side management in microgrid | |
Wang et al. | An integrated energy distribution system planning method with multiple energy storage systems | |
Zhou et al. | Wind Power Penetration Limit Calculation of Black-Start Based on Copula Theory | |
Gerlach et al. | Experts versus Algorithms? Optimized Fuzzy Logic Energy Management of Autonomous PV Hybrid Systems with Battery and H2 Storage. Energies 2021, 14, 1777 | |
Xingyou et al. | The Utility Analysis of LHSS Combined SMES System Considering Capacity Configuration |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210720 Termination date: 20211031 |