CN107239863A - 电网安全约束的鲁棒机组组合方法 - Google Patents
电网安全约束的鲁棒机组组合方法 Download PDFInfo
- Publication number
- CN107239863A CN107239863A CN201710515679.1A CN201710515679A CN107239863A CN 107239863 A CN107239863 A CN 107239863A CN 201710515679 A CN201710515679 A CN 201710515679A CN 107239863 A CN107239863 A CN 107239863A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- mtd
- msub
- mtr
- 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 26
- 238000012795 verification Methods 0.000 claims abstract description 30
- 238000003462 Bender reaction Methods 0.000 claims abstract description 22
- 238000005520 cutting process Methods 0.000 claims abstract description 20
- 230000005611 electricity Effects 0.000 claims abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 230000009467 reduction Effects 0.000 claims abstract description 10
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims description 14
- 239000000446 fuel Substances 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 9
- 238000009987 spinning Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 4
- 230000002040 relaxant effect Effects 0.000 claims description 4
- 208000030623 Low phospholipid-associated cholelithiasis Diseases 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010977 unit operation 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
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种电网安全约束的鲁棒机组组合方法,步骤如下:S1,预测风电出力,并对其进行拉丁超立方抽样,生成场景集,然后通过场景缩减技术对场景集缩减。S2,对缩减后的场景集按照其相应概率大小从大到小排序。S3,根据风电的预测出力求解主问题,即确定性考虑安全约束的机组组合问题。S4,对求得机组组合策略对场景逐一进行安全性校验。若通过校验,则计算该启停策略对其他场景集的切负荷和弃风所产生成本,并将该场景加入“鲁棒置信区间”,并则对下一个场景进行安全性校验。若不通过,则生成benders 割返回主问题,重复步骤S3。S5,当对所有场景进行过安全性校验后,总成本最小的即为最优策略。
Description
技术领域
本发明属于电力系统调度技术领域,特别是涉及一种电网安全约束的鲁棒机组组合方法。
背景技术
为应对化石能源枯竭和环境恶化问题,间歇性可再生能源如风电和太阳能在电力系统中的渗透率逐步提高,给系统安全与经济运行带来了问题甚至挑战。由于很难对间歇性可再生能源的发电出力进行准确预测,这给预先制定系统运行方式和实时调度都带来了困难。这样,针对如此不确定环境,如何在保证系统安全性的前提下确定兼顾系统经济性的运行和调度策略就是值得研究的重要问题。
针对不确定性参数的特征,可以采用概率、区间数、模糊数等进行模拟。概率模型一般采用蒙特卡洛仿真(Monte Carlo Simulation,MCS)生成场景和广义代数建模系统(General Algebraic Modeling System,GAMS)进行场景缩减,来描述和处理电力系统中的不确定性因素。区间数模型只需关注不确定性变量的上下界,通过区间优化方法得到乐观解和悲观解。在此背景下,重点研究考虑安全约束的情况下,使用场景模拟风电出力不确定来对电力系统机组组合鲁棒优化。
发明内容
基于此,本发明的目的在于提供一种电网安全约束的鲁棒机组组合方法,能够对电力系统机组组合鲁棒优化。
为实现上述目的,本发明采用如下技术方案:一种电网安全约束的鲁棒机组组合方法,包括以下步骤:
S1,预测风电出力,并对其进行拉丁超立方抽样,生成场景集,然后通过场景缩减技术对场景集缩减;
S2,对缩减后的场景集按照其相应概率大小从大到小排序;
S3,根据风电的预测出力求解主问题,即确定性考虑安全约束的机组组合问题;
S4,对求得机组组合策略对场景逐一进行安全性校验;若通过校验,则计算该启停策略对其他场景集的切负荷和弃风所产生成本,并将该场景加入“鲁棒置信区间”,并则对下一个场景进行安全性校验;若不通过,则生成benders割返回主问题,重复步骤S3;
S5,当对所有场景进行过安全性校验后,总成本最小的即为最优策略,其对应的“鲁棒置信区间”的概率即为鲁棒置信度。
进一步,步骤S3中,在鲁棒优化的框架下,建立能够满足给定鲁棒置信区间的UC模型,为提高计算速度,对燃料成本曲线进行分段线性化处理,所构造的UC目标函数为:
式中:Ng为发电机组总数;T为一个调度周期内所包括的时段总数;K为机组燃料成本的分段数目;li,k为机组燃料成本的分段线性系数(i=1,2,...,Ng;k=1,2,...,K);为发电机组i在时段t的状态,0和1分别表示停运和运行;表示火电机组i在时段t是否开机,其为1表示开机,否则为0;为火电机组i在时段t的开机成本;为火电机组i在时段t的有功出力;Fld为切负荷和弃风的损失;S为需要考虑的不确定性变量的场景集,SR为鲁棒置信区间所包含的场景集;ρs为场景s的概率,为场景s在时段t的切负荷功率,为场景s在时段t的弃风功率;
此外模型还需满足一定的约束条件,主要包括两部分:1)确定性调度模型中的约束条件;2)安全性校验中的约束条件,即鲁棒约束条件。
进一步,确定性调度的约束条件,
(1)功率平衡约束
式中:Nw为风电场数;为预测的风电场j在时段t的出力;为节点k在时段t的负荷;为时段t的线路损耗,可采用交流潮流模型求取;
(2)发电机组出力约束
式中:和分别为机组i的出力下限和上限;
(3)发电机组最小启停时间约束
式中:和分别为给定的机组i最小运行时间和最小停机时间;
(4)机组爬坡约束
式中:和分别为发电机组i的爬坡和滑坡速率;
(5)机组旋转备用约束
式中r为系统的旋转备用率;
(6)支路潮流约束
式中:Fij为支路ij可输送的有功功率极限;为支路ij在时段t输送的有功功率。
进一步,鲁棒约束条件,风电出力置信区间用场景集SR表示,需要满足的约束条件如下:
(1)功率平衡约束
式中:和分别为在场景s下时段t的常规机组出力、风电出力和网络损耗;
(2)发电机出力约束
(3)机组爬坡约束
(4)支路潮流约束
式中:为在场景s下支路ij的有功功率;
式(8)-(11)分别对应式(2)、(3)、(5)和(7),表示在场景s下的相应约束。
进一步,步骤S4中,采用以下求解方法:
(1)潮流线性化
在每个优化时段,采用LPAC方法对交流潮流方程做近似线性化处理;首先给定sinθij≈θij,Vi≈1.0,对cosθij采用多面松弛技术进行简化处理,这样,支路潮流方程可简化为:
Pij=gij-gij〈cosθij>-bijθij (12)
式中:gij和bij分别为线路ij的导纳和容纳,<cosθij>可由式(14)和式(15)表示:
式中:<cosθij>是cosθij的近似形式;是节点i的电压Vi的偏移量, 为预估的相角差最大取值范围;H为切线个数;d为相邻切线的角度差,
此外,每个节点需要满足功率平衡约束,发电机节点还需要满足无功出力上下限约束并考虑约束起作用时由PV节点向PQ节点的转换;
(2)Benders分解
经过上述简化处理后,可得到MILP模型,采用Benders分解将该MILP模型分解为主问题和子问题,主问题是确定性的UC问题,子问题则对主问题生成的UC方案进行安全性检验;若由主问题获得的UC方案满足安全性校验,则取其为最终结果;若不满足,则产生Benders割并反馈到主问题继续求解;
主问题:
其中:w(x)≤0为Benders割,w(x)的表达式在下述子问题中给出,即式(22)。
子问题:
对主问题求得的机组启停计划进行安全校验,引入非负松弛变量和构建如下优化模型:
由式(17)求得的vs如果为零,则表示满足所有约束,否则就有约束没能满足,
其中:为主问题得到的启停策略,和分别为相应约束条件的对偶变量。当子问题的目标函数值vs等于零,可以认为所有约束均满足,即通过安全性校验。否则必定有约束不满足。如果vs大于预先设定的阈值ε,则认为不满足安全性校验,需由式(22)求得的Benders割w(x)返回到主问题继续求解。
式中:为场景s下子问题的目标函数值,
由子问题模型可以看出最小弃风量和最小切负荷量至少有一个为零;当此时可以功率平衡约束中负荷过大,可通过切负荷实现功率平衡;而当此时可认为风电出力过高,需要通过弃风实现功率平衡;其最小弃风量和切负荷量可认为相应松弛变量取值,如式(23)所示:
本发明首先建立基于场景生成的鲁棒优化SCUC(Robust Security ConstrainedUnit Commitment,RSCUC)模型,由此获得的鲁棒机组组合策略满足给定的置信度,对处于置信区间之外的极端场景则采取弃风或切负荷等不得已措施来维持系统功率平衡,从而在系统运行的经济性和保守性之间合理折衷。之后,采用Benders分解法求解所建模型,将该问题分解为主问题和子问题。其中,主问题为确定性的SCUC问题;子问题则对考虑风电场出力随机变化时的系统状态进行安全性校验,若通过校验则表明所求得的SCUC策略满足鲁棒性约束,否则就生成相应的安全约束即Benders割并反馈给主问题。最后,采用修改的IEEE39节点系统来说明所提方法的有效性。
附图说明
图1是本发明一实施例的机组组合鲁棒优化流程示意图;
图2是本发明一实施例的采用IEEE 10机39节点系统的拓扑结构图;
图3是本发明风电场出力原始场景集和缩减后的场景集示意图;
图4为基准情形下的成本曲线图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明,不能理解为对本专利的限制。
请参照图1,在一个实施例中,本发明提出的机组组合鲁棒优化方法包括步骤如下:
S1,预测风电出力,并对其进行拉丁超立方抽样,生成场景集,然后通过场景缩减技术对场景集缩减;
S2,对缩减后的场景集按照其相应概率大小从大到小排序;
S3,根据风电的预测出力求解主问题,即确定性考虑安全约束的机组组合问题;
S4,对求得机组组合策略对场景逐一进行安全性校验;若通过校验,则计算该启停策略对其他场景集的切负荷和弃风所产生成本,并将该场景加入“鲁棒置信区间”,并则对下一个场景进行安全性校验;若不通过,则生成benders割返回主问题,重复步骤S3;
S5,当对所有场景进行过安全性校验后,总成本最小的即为最优策略,其对应的“鲁棒置信区间”的概率即为鲁棒置信度。
在传统的机组最优组合(UC)模型中,一般以给定调度时段内包括机组燃料成本和机组组合成本的系统运行总费用最小为优化目标。随着风电和太阳能等间歇性电源容量渗透率的不断提高,需要对传统UC模型进行改进,引入惩罚项来表征弃风和切负荷等措施所导致的损失。针对这一问题,在鲁棒优化的框架下,本节建立能够满足给定鲁棒置信区间的UC模型。为提高计算速度,对燃料成本曲线进行分段线性化处理。所构造的UC目标函数为:
式中:Ng为发电机组总数;T为一个调度周期内所包括的时段总数;K为机组燃料成本的分段数目;li,k为机组燃料成本的分段线性系数(i=1,2,...,Ng;k=1,2,...,K);为发电机组i在时段t的状态,0和1分别表示停运和运行;表示火电机组i在时段t是否开机,其为1表示开机,否则为0;为火电机组i在时段t的开机成本;为火电机组i在时段t的有功出力;Fld为切负荷和弃风的损失;S为需要考虑的不确定性变量的场景集,SR为鲁棒置信区间所包含的场景集;ρs为场景s的概率,为场景s在时段t的切负荷功率,为场景s在时段t的弃风功率;
此外模型还需满足一定的约束条件,主要包括两部分:1)确定性调度模型中的约束条件;2)安全性校验中的约束条件,即鲁棒约束条件。
一、确定性调度的约束条件
(1)功率平衡约束
式中:Nw为风电场数;为预测的风电场j在时段t的出力;为节点k在时段t的负荷;为时段t的线路损耗,可采用交流潮流模型求取,在此不再赘述;
(2)发电机组出力约束
式中:和分别为机组i的出力下限和上限;
(3)发电机组最小启停时间约束
式中:和分别为给定的机组i最小运行时间和最小停机时间;
(4)机组爬坡约束
式中:和分别为发电机组i的爬坡和滑坡速率;
(5)机组旋转备用约束
式中r为系统的旋转备用率;
(6)支路潮流约束
式中:Fij为支路ij可输送的有功功率极限;为支路ij在时段t输送的有功功率。
二、鲁棒约束条件
这里的鲁棒性指UC策略对于风电机组在其出力的鲁棒置信区间变化时均能找到相应的经济调度方式来消纳风电出力。风电出力置信区间用场景集SR表示,需要满足的约束条件如下:
(1)功率平衡约束
式中:和分别为在场景s下时段t的常规机组出力、风电出力和网络损耗;
(2)发电机出力约束
(3)机组爬坡约束
(4)支路潮流约束
式中:为在场景s下支路ij的有功功率;
式(8)-(11)分别对应式(2)、(3)、(5)和(7),表示在场景s下的相应约束。
求解方法:
(1)潮流线性化
这里采用多面松弛技术,将交流潮流松弛为线性模型,在提高求解效率的同时,结果也比较精确。在每个优化时段,可采用LPAC方法对交流潮流方程做近似线性化处理。首先给定sinθij≈θij,Vi≈1.0,对cosθij采用多面松弛技术进行简化处理,这样,支路潮流方程可简化为:
Pij=gij-gij<cosθij>-bijθij(12)
式中:gij和bij分别为线路ij的导纳和容纳,<cosθij>可由式(14)和式(15)表示:
式中:<cosθij>是cosθij的近似形式;是节点i的电压Vi的偏移量, 为预估的相角差最大取值范围;H为切线个数;d为相邻切线的角度差,
此外,每个节点需要满足功率平衡约束,发电机节点还需要满足无功出力上下限约束并考虑约束起作用时由PV节点向PQ节点的转换;
(2)Benders分解
经过上述简化处理后,可得到MILP模型,采用Benders分解将该MILP模型分解为主问题和子问题,主问题是确定性的UC问题,子问题则对主问题生成的UC方案进行安全性检验;若由主问题获得的UC方案满足安全性校验,则取其为最终结果;若不满足,则产生Benders割并反馈到主问题继续求解;
主问题:
其中:w(x)≤0为Benders割,w(x)的表达式在下述子问题中给出,即式(22)。
子问题:
对主问题求得的机组启停计划进行安全校验,引入非负松弛变量和构建如下优化模型:
由式(17)求得的vs如果为零,则表示满足所有约束,否则就有约束没能满足,
其中:为主问题得到的启停策略,和分别为相应约束条件的对偶变量。当子问题的目标函数值vs等于零,可以认为所有约束均满足,即通过安全性校验。否则必定有约束不满足。如果vs大于预先设定的阈值ε,则认为不满足安全性校验,需由式(22)求得的Benders割w(x)返回到主问题继续
求解,
式中:为场景s下子问题的目标函数值,
由子问题模型可以看出最小弃风量和最小切负荷量至少有一个为零;当此时可以功率平衡约束中负荷过大,可通过切负荷实现功率平衡;而当此时可认为风电出力过高,需要通过弃风实现功率平衡;其最小弃风量和切负荷量可认为相应松弛变量取值,如式(23)所示:
实施案例
1.数据处理
采用IEEE 10机39节点系统对所提出的方法进行说明,该系统的拓扑图见附图2。在具有双核3.2GHz CPU和4GB内存的个人计算机上进行了实现,并由商业求解器CPLEX12.5求解RSCUC模型中的混合整数规划问题。对机组燃料成本曲线进行分段线性化处理,分为3段。对于机组每小时的最大爬坡和滑坡功率分别给定为相应机组最大功率的1/3和1/6。Benders分解中阈值ε选取为10-3;目标函数中的Fld数值为800USD/MWh。给定与交流潮流线性化相关的参数为:h=10。风电场在节点29接入系统,每时段系统总负荷和风电预测出力如表1,各节点负荷占总负荷的比例如表2所示。
表1各时段系统负荷和风电出力的预测值
Table 1The forecasted values of the loads and wind power outputs
表2各负荷节点有功负荷占系统总有功负荷的百分比
Table 2Active power demand at each load bus over the total activeload in the system
假定风电出力服从正态分布,在时刻t的期望值μt为预测值,预测误差σt为μt的10%。为方便后续叙述,将预测误差σt为μt的10%这种情形称为“基准情形”。给定风电场出力的置信概率α为99.74%,对应的Kα=3。这样,风电出力置信区间为:
给定由5000个风电出力原始场景经场景缩减后得到100个场景。图2展示了风电场出力原始场景集和缩减后的场景集。把参数按上述取值时的运行场景作为基准情形。
2.基准情形下的计算结果
图3展示了基准情形下的计算结果。可以看出,随着鲁棒置信水平的提高,需要满足的场景增多,RSCUC成本随之增大。基准情形下一天24小时的最优RSCUC的机组成本为$411422.3,此时相应的切机和切负荷导致的赔偿费用的期望值为$280.3,总成本为$414702.6,鲁棒置信度为59%。表3列出了基准情形下的机组最优组合策略,1和0分别表示机组运行和停机。
当给定鲁棒置信度为100%时,总成本为$419923.9,比上述RSCUC最优策略的成本高$5221.3。而当给定鲁棒置信度为0时,即不计及风电出力预测误差时,RSCUC的成本为$417412.9,比最优解对应的成本高$2990.6。
表3基准情形下的机组最优组合策略
Table 3Optimal unit commitment strategy for the base case
3.风电出力预测误差对RSCUC结果的影响
本小节分析风电出力预测误差对RSCUC结果的影响,并对模型和常规鲁棒优化模型进行对比分析。定义两种RSCUC模式:
模式1:采用基于常规鲁棒优化模型的RSCUC,不考虑弃风和切负荷。这种模式为所提模型的极端形式,对应的鲁棒置信度为100%。
模式2:采用所提出模型。
用σ/μ表示风电出力预测误差,给定不同的σ/μ值,针对基态情形计算这两种模式下RSCUC优化结果对应的运行总成本。场景集取缩减后的100个场景。计算结果列于表4。
表4风电出力预测误差对RSCUC运行总成本的影响
Table 4Impacts of wind power prediction error on the overalloperation cost of RSCUC
由表4可以看出,当σ/μ=0即风电出力预测完全准确时,此时即为确定性的SCUC问题,两种模式求得的结果相同。只要σ/μ大于0,由模式2求得的运行成本总是小于模式1相关的成本,且随着σ/μ的增大,两种模式的优化结果相差增大。随着σ/μ增大,两种模式所求的运行总成本均增大,这符合预期,也说明了风电出力预测精度影响系统运行的经济性。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (5)
1.一种电网安全约束的鲁棒机组组合方法,其特征在于,包括以下步骤:
S1,预测风电出力,并对其进行拉丁超立方抽样,生成场景集,然后通过场景缩减技术对场景集缩减;
S2,对缩减后的场景集按照其相应概率大小从大到小排序;
S3,根据风电的预测出力求解主问题,即确定性考虑安全约束的机组组合问题;
S4,对求得机组组合策略对场景逐一进行安全性校验;若通过校验,则计算该启停策略对其他场景集的切负荷和弃风所产生成本,并将该场景加入“鲁棒置信区间”,并则对下一个场景进行安全性校验;若不通过,则生成benders割返回主问题,重复步骤S3;
S5,当对所有场景进行过安全性校验后,总成本最小的即为最优策略,其对应的“鲁棒置信区间”的概率即为鲁棒置信度。
2.根据权利要求1所述的电网安全约束的鲁棒机组组合方法,其特征在于,步骤S3中,在鲁棒优化的框架下,建立能够满足给定鲁棒置信区间的UC模型,为提高计算速度,对燃料成本曲线进行分段线性化处理,所构造的UC目标函数为:
<mrow>
<mi>min</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>k</mi>
<mi>K</mi>
</munderover>
<msub>
<mi>l</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>c</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>e</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&rsqb;</mo>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mo>&Element;</mo>
<mi>S</mi>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>F</mi>
<mrow>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>s</mi>
<mi>h</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:Ng为发电机组总数;T为一个调度周期内所包括的时段总数;K为机组燃料成本的分段数目;li,k为机组燃料成本的分段线性系数(i=1,2,...,Ng;k=1,2,...,K);为发电机组i在时段t的状态,0和1分别表示停运和运行;表示火电机组i在时段t是否开机,其为1表示开机,否则为0;为火电机组i在时段t的开机成本;为火电机组i在时段t的有功出力;Fld为切负荷和弃风的损失;S为需要考虑的不确定性变量的场景集,SR为鲁棒置信区间所包含的场景集;ρs为场景s的概率,为场景s在时段t的切负荷功率,为场景s在时段t的弃风功率;
此外模型还需满足一定的约束条件,主要包括两部分:1)确定性调度模型中的约束条件;2)安全性校验中的约束条件,即鲁棒约束条件。
3.根据权利要求2所述的电网安全约束的鲁棒机组组合方法,其特征在于,确定性调度的约束条件,
(1)功率平衡约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msubsup>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>k</mi>
<msub>
<mi>N</mi>
<mi>d</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:Nw为风电场数;为预测的风电场j在时段t的出力;为节点k在时段t的负荷;为时段t的线路损耗,可采用交流潮流模型求取;
(2)发电机组出力约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>min</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:和分别为机组i的出力下限和上限;
(3)发电机组最小启停时间约束
<mrow>
<mtable>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
1
式中:和分别为给定的机组i最小运行时间和最小停机时间;
(4)机组爬坡约束
<mrow>
<mtable>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>+</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>-</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:和分别为发电机组i的爬坡和滑坡速率;
(5)机组旋转备用约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>p</mi>
<mi>max</mi>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msubsup>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>k</mi>
<msub>
<mi>N</mi>
<mi>d</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中r为系统的旋转备用率;
(6)支路潮流约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:Fij为支路ij可输送的有功功率极限;为支路ij在时段t输送的有功功率。
4.根据权利要求3所述的电网安全约束的鲁棒机组组合方法,其特征在于,鲁棒约束条件,风电出力置信区间用场景集SR表示,需要满足的约束条件如下:
(1)功率平衡约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mi>k</mi>
<msub>
<mi>N</mi>
<mi>d</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>s</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:和分别为在场景s下时段t的常规机组出力、风电出力和网络损耗;
(2)发电机出力约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>min</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>s</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
(3)机组爬坡约束
<mrow>
<mtable>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>+</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>-</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>s</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
(4)支路潮流约束
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>s</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:为在场景s下支路ij的有功功率;
式(8)-(11)分别对应式(2)、(3)、(5)和(7),表示在场景s下的相应约束。
5.根据权利要求4所述的电网安全约束的鲁棒机组组合方法,其特征在于,步骤S4中,采用以下求解方法:
(1)潮流线性化
在每个优化时段,采用LPAC方法对交流潮流方程做近似线性化处理;首先给定sinθij≈θij,Vi≈1.0,对cosθij采用多面松弛技术进行简化处理,这样,支路潮流方程可简化为:
Pij=gij-gij<cosθij>-bijθij (12)
式中:gij和bij分别为线路ij的导纳和容纳,<cosθij〉可由式(14)和式(15)表示:
<mrow>
<mo><</mo>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>></mo>
<mo>&le;</mo>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mi>d</mi>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>h</mi>
<mi>d</mi>
<mo>+</mo>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mi>d</mi>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>h</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>H</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
2
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&le;</mo>
<mo><</mo>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>></mo>
<mo>&le;</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:〈cosθij〉是cosθij的近似形式;是节点i的电压Vi的偏移量, 为预估的相角差最大取值范围;H为切线个数;d为相邻切线的角度差,
此外,每个节点需要满足功率平衡约束,发电机节点还需要满足无功出力上下限约束并考虑约束起作用时由PV节点向PQ节点的转换;
(2)Benders分解
经过上述简化处理后,可得到MILP模型,采用Benders分解将该MILP模型分解为主问题和子问题,主问题是确定性的UC问题,子问题则对主问题生成的UC方案进行安全性检验;若由主问题获得的UC方案满足安全性校验,则取其为最终结果;若不满足,则产生Benders割并反馈到主问题继续求解;
主问题:
<mrow>
<munder>
<mi>min</mi>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</munder>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<munder>
<mi>&Sigma;</mi>
<mi>k</mi>
</munder>
<msub>
<mi>l</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>c</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>e</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&rsqb;</mo>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mo>&Element;</mo>
<mi>S</mi>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mi>T</mi>
</mrow>
</munderover>
<msub>
<mi>F</mi>
<mrow>
<mi>l</mi>
<mi>d</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>s</mi>
<mi>h</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>w</mi>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>&le;</mo>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中:w(x)≤0为Benders割,w(x)的表达式在下述子问题中给出,即式(22)。
子问题:
对主问题求得的机组启停计划进行安全校验,引入非负松弛变量和构建如下优化模型:
<mrow>
<msup>
<mi>v</mi>
<mi>s</mi>
</msup>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
由式(17)求得的vs如果为零,则表示满足所有约束,否则就有约束没能满足,
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<munderover>
<mi>&Sigma;</mi>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mi>j</mi>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mi>k</mi>
<msub>
<mi>N</mi>
<mi>d</mi>
</msub>
</munderover>
<msubsup>
<mi>p</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>p</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>max</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>max</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&RightArrow;</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>min</mi>
</msubsup>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>&RightArrow;</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>+</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<mo>&RightArrow;</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>3</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>-</mo>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<mo>&RightArrow;</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>4</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
其中:为主问题得到的启停策略,和分别为相应约束条件的对偶变量。当子问题的目标函数值vs等于零,可以认为所有约束均满足,即通过安全性校验。否则必定有约束不满足。如果vs大于预先设定的阈值ε,则认为不满足安全性校验,需由式(22)求得的Benders割w(x)返回到主问题继续求解,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mover>
<mi>v</mi>
<mo>^</mo>
</mover>
<mi>s</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mi>i</mi>
<msub>
<mi>N</mi>
<mi>g</mi>
</msub>
</munderover>
<mo>{</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>min</mi>
</msubsup>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>+</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>3</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mo>-</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mi>max</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>&theta;</mi>
<mrow>
<mn>4</mn>
<mi>i</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:为场景s下子问题的目标函数值,
由子问题模型可以看出最小弃风量和最小切负荷量至少有一个为零;当此时可以功率平衡约束中负荷过大,可通过切负荷实现功率平衡;而当此时可认为风电出力过高,需要通过弃风实现功率平衡;其最小弃风量和切负荷量可认为相应松弛变量取值,如式(23)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>s</mi>
<mi>h</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
4
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710235984 | 2017-04-12 | ||
CN2017102359845 | 2017-04-12 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107239863A true CN107239863A (zh) | 2017-10-10 |
CN107239863B CN107239863B (zh) | 2020-07-14 |
Family
ID=59989984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710515679.1A Active CN107239863B (zh) | 2017-04-12 | 2017-06-29 | 一种电网安全约束的鲁棒机组组合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107239863B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742901A (zh) * | 2017-11-14 | 2018-02-27 | 清华大学 | 考虑压缩空气储能的风电并网机组组合方法及装置 |
CN107947164A (zh) * | 2017-11-30 | 2018-04-20 | 三峡大学 | 一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
CN108054757A (zh) * | 2017-12-22 | 2018-05-18 | 清华大学 | 一种内嵌无功和电压的n-1闭环安全校核方法 |
CN108258695A (zh) * | 2018-02-28 | 2018-07-06 | 东南大学 | 一种交直流混联微网的随机鲁棒耦合型优化调度方法 |
CN109409595A (zh) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | 一种园区多能互补系统日前调度方法 |
WO2019128012A1 (zh) * | 2017-12-28 | 2019-07-04 | 东南大学 | 一种交直流混合微网鲁棒优化协调调度方法 |
CN110137955A (zh) * | 2019-05-21 | 2019-08-16 | 国网能源研究院有限公司 | 一种计及CVaR的鲁棒机组组合调度的决策方法 |
CN110311427A (zh) * | 2019-06-18 | 2019-10-08 | 华北电力大学 | 计及故障概率的两阶段n-k鲁棒故障约束机组组合方法 |
CN110533304A (zh) * | 2019-08-12 | 2019-12-03 | 湖南大学 | 一种电力系统负荷不确定性分析方法 |
CN112039058A (zh) * | 2020-08-18 | 2020-12-04 | 国网山东省电力公司临沂供电公司 | 基于风功率预测区间的机组组合方法、系统、介质及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102684224A (zh) * | 2012-05-25 | 2012-09-19 | 浙江大学 | 一种求解考虑风电波动性的机组组合方法 |
CN104182910A (zh) * | 2014-08-25 | 2014-12-03 | 国家电网公司 | 一种涉及相关性的风电出力场景构建方法 |
CN104716670A (zh) * | 2015-04-01 | 2015-06-17 | 国家电网公司 | 光伏并网下基于网络安全约束的机组组合方法 |
CN105024398A (zh) * | 2015-06-29 | 2015-11-04 | 国电南瑞科技股份有限公司 | 一种基于最优风电置信度的优化调度方法 |
-
2017
- 2017-06-29 CN CN201710515679.1A patent/CN107239863B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102684224A (zh) * | 2012-05-25 | 2012-09-19 | 浙江大学 | 一种求解考虑风电波动性的机组组合方法 |
CN104182910A (zh) * | 2014-08-25 | 2014-12-03 | 国家电网公司 | 一种涉及相关性的风电出力场景构建方法 |
CN104716670A (zh) * | 2015-04-01 | 2015-06-17 | 国家电网公司 | 光伏并网下基于网络安全约束的机组组合方法 |
CN105024398A (zh) * | 2015-06-29 | 2015-11-04 | 国电南瑞科技股份有限公司 | 一种基于最优风电置信度的优化调度方法 |
Non-Patent Citations (2)
Title |
---|
ANTONIO FRANGIONI等: "Tighter Approximated MILP Formulations for Unit Commitment Problems", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
王颖等: "计及风电最优置信度的机组组合优化方法", 《电网技术》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742901A (zh) * | 2017-11-14 | 2018-02-27 | 清华大学 | 考虑压缩空气储能的风电并网机组组合方法及装置 |
CN107947164A (zh) * | 2017-11-30 | 2018-04-20 | 三峡大学 | 一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
CN112186765A (zh) * | 2017-11-30 | 2021-01-05 | 三峡大学 | 一种机组组合决策的日前调度模型的建模方法 |
CN112186765B (zh) * | 2017-11-30 | 2022-06-17 | 三峡大学 | 一种机组组合决策的日前调度模型的建模方法 |
CN107947164B (zh) * | 2017-11-30 | 2020-11-17 | 三峡大学 | 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
CN112186766B (zh) * | 2017-11-30 | 2022-07-01 | 三峡大学 | 一种基于Benders分解的多阶段分解方法 |
CN112186766A (zh) * | 2017-11-30 | 2021-01-05 | 三峡大学 | 一种基于Benders分解的多阶段分解方法 |
CN108054757B (zh) * | 2017-12-22 | 2019-08-16 | 清华大学 | 一种内嵌无功和电压的n-1闭环安全校核方法 |
CN108054757A (zh) * | 2017-12-22 | 2018-05-18 | 清华大学 | 一种内嵌无功和电压的n-1闭环安全校核方法 |
WO2019128012A1 (zh) * | 2017-12-28 | 2019-07-04 | 东南大学 | 一种交直流混合微网鲁棒优化协调调度方法 |
CN108258695A (zh) * | 2018-02-28 | 2018-07-06 | 东南大学 | 一种交直流混联微网的随机鲁棒耦合型优化调度方法 |
CN109409595A (zh) * | 2018-10-19 | 2019-03-01 | 南京南瑞继保电气有限公司 | 一种园区多能互补系统日前调度方法 |
CN110137955A (zh) * | 2019-05-21 | 2019-08-16 | 国网能源研究院有限公司 | 一种计及CVaR的鲁棒机组组合调度的决策方法 |
CN110311427A (zh) * | 2019-06-18 | 2019-10-08 | 华北电力大学 | 计及故障概率的两阶段n-k鲁棒故障约束机组组合方法 |
CN110311427B (zh) * | 2019-06-18 | 2020-12-01 | 华北电力大学 | 计及故障概率的两阶段n-k鲁棒故障约束机组组合方法 |
CN110533304A (zh) * | 2019-08-12 | 2019-12-03 | 湖南大学 | 一种电力系统负荷不确定性分析方法 |
CN110533304B (zh) * | 2019-08-12 | 2022-12-09 | 湖南大学 | 一种电力系统负荷不确定性分析方法 |
CN112039058A (zh) * | 2020-08-18 | 2020-12-04 | 国网山东省电力公司临沂供电公司 | 基于风功率预测区间的机组组合方法、系统、介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107239863B (zh) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107239863A (zh) | 电网安全约束的鲁棒机组组合方法 | |
Akbari-Dibavar et al. | Economic-emission dispatch problem in power systems with carbon capture power plants | |
Lara et al. | Deterministic electric power infrastructure planning: Mixed-integer programming model and nested decomposition algorithm | |
CN105046395B (zh) | 一种含多类型新能源的电力系统日内滚动计划编制方法 | |
Guerra et al. | An optimization framework for the integrated planning of generation and transmission expansion in interconnected power systems | |
Palmintier | Incorporating operational flexibility into electric generation planning: Impacts and methods for system design and policy analysis | |
Fu et al. | Research on joint optimal dispatching method for hybrid power system considering system security | |
CN104779611B (zh) | 基于集中式和分布式双层优化策略的微电网经济调度方法 | |
CN106485352A (zh) | 一种多类电源发电计划组合决策方法 | |
CN108599268B (zh) | 一种考虑风电场时空关联约束的日前优化调度方法 | |
CN105322566A (zh) | 考虑预测误差时序分布的含风电机组组合模型建立方法 | |
CN103699941A (zh) | 一种电力系统调度运行年方案的制定方法 | |
CN107947244A (zh) | 新能源电站并网功率超前控制方法 | |
CN105305423A (zh) | 一种考虑间歇性能源不确定性的最优误差边界的确定方法 | |
CN106160000A (zh) | 一种分调‑省调两级协调的风电增发控制方法 | |
CN103956773A (zh) | 含风电系统机组的备用配置优化方法 | |
Wang et al. | A novel security stochastic unit commitment for wind-thermal system operation | |
CN111082466A (zh) | 考虑风电不确定性的新能源接入与网架扩建优化方法 | |
CN105528668A (zh) | 一种风电并网电力系统动态环境经济调度方法 | |
CN102622269A (zh) | 基于JADE的智能电网发电调度多Agent系统 | |
CN105305485A (zh) | 一种消纳大规模间歇性能源的安全约束经济调度方法 | |
Allen et al. | An optimization framework for solving integrated planning and scheduling problems for dense energy carriers | |
CN103633641A (zh) | 一种考虑风电接纳的中长期交易运营计划获取方法 | |
CN105226707A (zh) | 一种基于Shapley值风电并网系统输电固定成本的分摊方法 | |
CN112994087B (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 |