CN107947164B - 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 - Google Patents
考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 Download PDFInfo
- Publication number
- CN107947164B CN107947164B CN201711238480.5A CN201711238480A CN107947164B CN 107947164 B CN107947164 B CN 107947164B CN 201711238480 A CN201711238480 A CN 201711238480A CN 107947164 B CN107947164 B CN 107947164B
- Authority
- CN
- China
- Prior art keywords
- day
- matrix
- ahead
- correlation
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 19
- 238000003462 Bender reaction Methods 0.000 claims abstract description 13
- 230000006870 function Effects 0.000 claims description 48
- 239000011159 matrix material Substances 0.000 claims description 45
- 238000005070 sampling Methods 0.000 claims description 29
- 230000009466 transformation Effects 0.000 claims description 9
- 238000010248 power generation Methods 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- 238000005457 optimization Methods 0.000 abstract description 6
- 230000009194 climbing Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000009977 dual effect Effects 0.000 description 3
- 108010074506 Transfer Factor Proteins 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 208000011580 syndromic disease Diseases 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000006386 memory function Effects 0.000 description 1
Images
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/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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/067—Enterprise or organisation modelling
-
- 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
-
- 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
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Entrepreneurship & Innovation (AREA)
- Power Engineering (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明属于电网调度领域,具体是一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,用来求解含新能源的电力系统动态经济调度问题。提出了一种综合考虑负荷、风电以及光伏出力不确定性及概率相关性的日前鲁棒调度方法。首先构建考虑多重不确定性因素及概率相关性的改进鲁棒优化调度模型;然后利用Cholesky分解法将具有相关性的随机样本转换为相互独立的随机样本,从而基于样本特征直接确定最坏场景;最后利用Benders分解法对模型进行求解。本发明所提方法可以在多重不确定性因素下,保证日前调度计划鲁棒性的同时,有效提升其经济性,而基于Cholesky分解的最坏场景确定方法也有效提升了鲁棒调度模型的紧凑性,使其计算效率得到显著提升。
Description
技术领域
本发明一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,涉及电力系统调度领域。
背景技术
风电及光伏是无污染、绿色的可再生能源,其分布广泛、能量密度高,适合大规模开发,因此,风力和光伏发电技术已受到了世界各国的高度重视。然而,由于其出力具有随机性与波动性的特点,大规模接入电网将会给传统的调度方法带来极大的挑战,因此,研究多种新能源大规模接入下的电力系统日前调度方法具有重要的理论价值与现实意义。
目前很多专家学者从不同角度研究了新能源接入下的电力系统日前调度问题,但普遍仅考虑了单一不确定性变量,然而实际电力系统中包含诸如风电出力、光伏出力、负荷预测误差等多重不确定性因素,现有仅考虑单一不确定性的日前调度显然难以保证其决策的有效性并影响系统运行的经济性。因此在调度问题中考虑多重不确定性因素的影响已经成为近年来专家学者们研究的热点。
发明内容
针对现有方法的不足,本发明提出了一种综合考虑负荷、风电以及光伏不确定性及相关性的改进鲁棒调度方法,该方法首先构造考虑多种不确定性因素的日前鲁棒优化调度模型,然后引入Cholesky分解法将具有相关性的随机样本转换为相互独立的随机样本,并基于样本特征直接确定最坏场景,最后进行鲁棒机组组合求解。
本发明采取的技术方案为:
一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,包括以下步骤:
1)考虑多重随机因素的日前鲁棒调度建模;
2)最坏场景求解;
3)通过Benders分解法来对模型进行求解;
在步骤1)中,首先考虑多重随机因素的日前鲁棒调度建模;然后进行基本场景下的日前调度建模;再进行不确定场景下的日前调度建模;
在步骤2)中首先进行随机因素的概率密度函数建模,然后样本抽样,再进行正交转换矩阵推导;接着进行最坏场景求取;
在步骤3)中,首先进行UC主问题建模,然后进行基本场景下的安全子问题建模;再进行考虑不确定性的安全子问题建模。
在步骤1)中,多重随机因素包括风电、光伏及负荷预测误差因素。
在步骤1)中,将机组组合决策分为基本场景和最坏场景分别建模。
上述基本场景以不确定性因素功率预测值为基础,以系统总运行成本最小为目标,同时考虑系统在确定性环境下的各种常规约束条件。
上述最坏场景以不确定性电源出力的最大波动出力为基础,考虑不确定性约束条件。
在步骤2)中,先寻求最坏场景,然后直接利用最坏场景进行鲁棒性校核。
首先利用非参数核密度估计分别构建风电、光伏以及负荷的概率密度函数,然后利用拉丁超立方抽样生成样本,最后采用Cholesky分解法将上述具有相关性的随机样本转换为相互独立的随机样本,并以此为基础确定最坏场景。
已知基于历史数据的负荷样本个数为n,则基于非参数核密度估计方法构建负荷的概率密度模型为:
式中:φ(Pd)为负荷的概率密度函数;K(Pd,l)为核函数;Pdm为负荷样本中的第m个样本值, l为带宽。
选择高斯函数作为负荷概率密度模型的核函数,并对带宽l进行求解,从而得到系统负荷的概率密度函数φ(Pd)。
求得风电出力的概率密度函数φ(Pw)以及光伏出力概率密度函数φ(Pv)。如式(2)和式 (3)。
式中:φ(Pw)为负荷的概率密度函数;K(Pw,l)为核函数;Pwm为风电出力样本中的第m个样本值,l为带宽。
式中:φ(Pv)为负荷的概率密度函数;K(Pv,l)为核函数;Pvm为光伏出力样本中的第m个样本值,l为带宽。
采用拉丁超立方抽样进行分层采样。
设采样规模为N,Ym=Fm(Xm)表示第m个随机变量Xm的概率密度函数,其具体抽样过程如下:将区间[0,1]平均分为N等分,选取每个子区间的中间值,通过其反函数得到采样值所有随机变量采样完成后,则得到其样本矩阵。
利用相关系数矩阵来描述负荷、风电以及光伏发电间的相关性,设通过拉丁超立方抽样得到的样本矩阵为W=[w1,w2,…wl]T,其相关系数矩阵为Cw:
该矩阵各元素可由式(5)求得:
相关系数矩阵Cw是正定矩阵,则可以对系数矩阵进行Cholesky分解:
Cw=GGT (6)
式中:G为下三角矩阵,其中元素可由式(7)求得:
设存在一正交矩阵B,可将具有相关性的输入随机变量W转换为不相关的随机变量Y:
Y=BW (8)
由于不相关随机变量Y的相关系数矩阵CY为单位矩阵I,因而有:
CY=ρ(Y,YT)=ρ(BW,WTBT)=Bρ(W,WT)BT=BCWBT=I (9)
又由式(6)可得:
CY=BCWBT=BGGTBT=(BG)(BG)T=I (10)
由上式推导得到:
B=G-1 (11)
在已知有相关性的输入不确定量W的前提下,通过正交变换矩阵,可以将其变为不相关的随机变量Y。
通过得到的正交变换矩阵,将存在相关性的负荷、风电及光伏出力矩阵转换为独立矩阵,消除其间的相关性,进而通过最坏场景线性叠加的方法求得最坏场景。
采用Benders分解法对模型进行求解;将原问题分解为一个主问题和两个子问题,主问题为基本场景下的UC决策主问题,两个子问题分别为,基本场景下的网络安全校核子问题和最坏场景下的鲁棒校核子问题;
包括以下步骤:首先,进行UC主问题建模;然后,进行基本场景下的安全子问题建模;接着,考虑不确定性的安全子问题建模;
具体步骤如下:
1):UC主问题建模
UC主问题包括目标函数式(1)及约束条件式(3-9)以及所有生成的Benders割。
2):基本场景下的安全子问题建模
基本场景下的网络安全校核子问题如式(27),其通过校验UC主问题的潮流越限情况来确保调度方案的网络安全。
式中:λ1,l,t,λ2,l,t是网络安全约束的对偶变量;vl,t为松弛变量。
子问题中引入松弛变量vl,t的作用是当于约束条件不能满足时,用松弛变量暂时缓解网络安全约束,以保证子问题始终有解。若最终优化得出的vl,t大于给定的安全阈值,则表示主问题求得的最优机组组合方案不能满足网络安全约束,因此需要返回Benders割如下:
3):考虑不确定性的安全子问题建模
S.t.
0≤v1,lt,v2t,v3t (29)
如果不能满足安全约束,则返回Benders割式(30)到主问题,其作为约束条件,使得机组组合和出力方案在最坏场景上自适应调整满足网络安全约束和功率平衡。
通过上述步骤,完成考虑多重不确定性及相关性的电力系统日前鲁棒调度。
本发明的技术效果如下:
1)、与传统方法相比,本发明一方面可以有效计及多重不确定性因素及其相关性,因而可以在保证调度决策方案鲁棒性的同时提高系统运行的经济性;另一方面,改进的鲁棒调度模型有效避免了在海量场景下求解机组组合问题,具有较高的求解效率。
2)、本发明通过引入Cholesky分解法对鲁棒优化方法进行改进,从而无需通过多场景潮流计算直接确定最坏场景,有效提升了鲁棒优化模型的适用性和求解效率。
附图说明
图1本发明整体建模框图;
图2本发明最坏场景求取流程图;
图3本发明算法总体思路框图;
图4本发明实施例风电及光伏有功出力曲线;
图5为时刻1下的最坏场景对比表格;
图6为日前调度成本对比结果表格;
图7为计算效率对比表格。
具体实施方式
一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,包括以下步骤:
步骤1:基本场景下的日前调度建模
本发明基于鲁棒优化思想,构建考虑风电、光伏及负荷预测误差不确定性的日前调度模型。本发明将机组组合决策分为基本场景和最坏场景分别建模,基本场景以不确定性因素功率预测值为基础,以系统总运行成本最小为目标,同时考虑系统在确定性环境下的各种常规约束条件,从而保障调度决策的经济性;最坏场景以不确定性电源出力的最大波动出力为基础,考虑不确定性约束条件,从而保证基本场景下的决策方案在不确定环境下的鲁棒性。模型框图如图1所示。
步骤1.1:基本场景下的日前调度建模:包含目标函数及约束条件
1)确定包含了常规机组的燃料费用和机组启停成本的目标函数:
其中:
2)确立该模型中含有的约束条件,如:
(1)系统功率平衡约束,不计网损的情况下:
(2)常规机组出力约束:
(3)机组爬坡速率约束
上升爬坡速率约束
下降爬坡速率约束
式中:URi,DRi为常规机组爬坡功率限制。
(4)常规机组启停时间约束
(5)网络安全约束
步骤1.2:不确定场景下的日前调度建模
基本场景下求得的日前机组组合和机组出力计划应当保证系统在不确定性环境下的鲁棒性,因此,需要利用不确定性环境下的系统约束条件对日前计划进行校核。需要指出的是,在实际计算中可以先寻求最坏场景,然后直接利用最坏场景进行鲁棒性校核。
确立模型中的约束条件:
(1)系统功率平衡约束,不计网损的情况下:
(2)常规机组出力约束:
(3)机组爬坡速率约束
上升爬坡速率约束
下降爬坡速率约束
(4)旋转备用容量约束
(5)网络安全约束
步骤2:最坏场景求解
由步骤1可知,在实际计算中可以先寻求最坏场景,然后直接利用最坏场景进行鲁棒性校核。因此本发明基于Cholesky分解理论,提出了一种适用于多重相关随机性因素的最坏场景快速推求方法。首先利用非参数核密度估计分别构建风电、光伏以及负荷的概率密度函数,然后利用拉丁超立方抽样生成样本,最后采用Cholesky分解法将上述具有相关性的随机样本转换为相互独立的随机样本,并以此为基础确定最坏场景。
步骤2.1:随机因素的概率密度函数建模
已知基于历史数据的负荷样本个数为n,则基于非参数核密度估计方法构建负荷的概率密度模型为:
式中:φ(Pd)为负荷的概率密度函数;K(Pd,l)为核函数;Pdm为负荷样本中的第m个样本值, l为带宽。
本发明选择高斯函数作为负荷概率密度模型的核函数,并对带宽l进行求解,从而得到系统负荷的概率密度函数φ(Pd)。
同理可以求得风电出力的概率密度函数φ(Pw)以及光伏出力概率密度函数φ(Pv)。如式 (17)和式(18)。
式中:φ(Pw)为负荷的概率密度函数;K(Pw,l)为核函数;Pwm为风电出力样本中的第m个样本值,l为带宽。
式中:φ(Pv)为负荷的概率密度函数;K(Pv,l)为核函数;Pvm为光伏出力样本中的第m个样本值,l为带宽。
步骤2.2:样本抽样
拉丁超立方抽样是一种分层采样方法,其具有样本记忆功能,可避免抽取已经出现的样本。设采样规模为N,Ym=Fm(Xm)表示第m个随机变量Xm的概率密度函数。其具体抽样过程如下:将区间[0,1]平均分为N等分,选取每个子区间的中间值,通过其反函数得到采样值所有随机变量采样完成后,则得到其样本矩阵。
步骤2.3:正交转换矩阵推导
利用相关系数矩阵来描述负荷、风电以及光伏发电间的相关性,设通过拉丁超立方抽样得到的样本矩阵为W=[w1,w2,…wl]T,其相关系数矩阵为Cw:
该矩阵各元素可由式(20)求得:
由定义不难看出,本发明的相关系数矩阵Cw是正定矩阵,则可以对系数矩阵进行Cholesky分解:
Cw=GGT (6)
式中:G为下三角矩阵,其中元素可由式(22)求得。
设存在一正交矩阵B,可将具有相关性的输入随机变量W转换为不相关的随机变量Y:
Y=BW (8)
由于不相关随机变量Y的相关系数矩阵CY为单位矩阵I,因而有:
CY=ρ(Y,YT)=ρ(BW,WTBT)=Bρ(W,WT)BT=BCWBT=I (9)
又由式(21)可得:
CY=BCWBT=BGGTBT=(BG)(BG)T=I (10)
由上式推导得到:
B=G-1 (11)
在已知有相关性的输入不确定量W的前提下,通过正交变换矩阵,可以将其变为不相关的随机变量Y。
步骤2.4:最坏场景求取;
通过步骤2.3得到的正交变换矩阵,将存在相关性的负荷、风电及光伏出力矩阵转换为独立矩阵,消除其间的相关性,进而通过最坏场景线性叠加的方法求得最坏场景。最坏场景求取流程图如图2所示。
步骤3:模型求解
本发明提出了一种基于Benders分解的多阶段分解算法来对模型进行求解。
传统Benders分解法一般是将原问题分解为主次两个问题,而考虑到本发明模型的特殊性,将原问题分解为一个主问题和两个子问题,主问题为基本场景下的UC决策主问题,两个子问题分别为:基本场景下的网络安全校核子问题和最坏场景下的鲁棒校核子问题。其整体框架如图1所示。
步骤3.1:UC主问题建模
UC主问题包括目标函数式(1)及约束条件式(3-9)以及所有生成的Benders割。
步骤3.2:基本场景下的安全子问题建模
基本场景下的网络安全校核子问题如式(27),其通过校验UC主问题的潮流越限情况来确保调度方案的网络安全。
式中:λ1,l,t,λ2,l,t是网络安全约束的对偶变量;vl,t为松弛变量。
子问题中引入松弛变量vl,t的作用是当于约束条件不能满足时,用松弛变量暂时缓解网络安全约束,以保证子问题始终有解。若最终优化得出的vl,t大于给定的安全阈值,则表示主问题求得的最优机组组合方案不能满足网络安全约束,因此需要返回Benders割如下:
步骤3.3:考虑不确定性的安全子问题建模
S.t.
0≤v1,lt,v2t,v3t (29)
如果不能满足安全约束,则返回Benders割式(30)到主问题,其作为约束条件,使得机组组合和出力方案在最坏场景上自适应调整满足网络安全约束和功率平衡。
步骤4:确立运行模式;
为对比分析本发明所建调度模型的有效性与正确性,确立以下两种运行模式:
模式1:考虑多重不确定性但不考虑相关性的鲁棒日前调度。
模式2:考虑多重不确定性及其相关性的鲁棒日前调度模型。
具体的,本发明以修改的IEEE-118节点系统来验证所提模型的正确性。该系统包含了 54台常规火电机组、3个风电场以及1个光伏电站。其中风电场额定功率分别为100MW、200MW、250MW,位于5、9、48号节点;光伏电站的容量为300MW,位于20号节点。风电与光伏的有功出力曲线如图5所示。系统中常规机组正旋转备用需求为系统最大负荷的8%,负旋转备用需求为系统最小负荷的2%。安全校核子问题校验值的阈值均取 10-3MWh。采用拉丁超立方抽样对负荷、风电以及光伏出力每小时抽样100次,共计2400 组样本。相关计算均在英特尔酷睿i5-4460处理器3.20GHz,8G内存计算机上完成,采用 Matlab和Cplex 12.5对算例进行编程求解。
根据本发明所提最坏场景求解方法,求出模式1和模式2的不确定性因素最坏场景对比结果如图5所示。
由图5可知,如果不考虑不确定性因素间的概率相关性,通过简单线性叠加求出的最坏场景较为保守,这将会影响到日前调度决策的经济性,通过本文方法的精确计算可以看出,由于不确定性因素之间存在相关性,这种最坏场景其实是不可能出现的。
分别计算模式1和模式2下的日前调度方案,其结果如图6所示。
由图6可知,相较于不考虑多重不确定性因素概率相关性的传统日前调度方法,本发明提出的方法使本算例的日前调度成本降低了约38650$。其原因是,本发明方法充分考虑了风电、光伏与负荷间的概率相关性,有效避免了在鲁棒性校核过程中计及不可能发生的极端场景,从而使日前调度决策在保证系统鲁棒性的同时,降低其运行成本。
为比较本发明方法相较于传统方法的优势,分别采用标准Benders分解方法和本发明方法对同一算例进行仿真,并对比其计算效率。对比结果如图7所示。
由图7可知:本发明所提方法较计算效率得到了非常显著的提升,这是由于本发明提出的基于Cholesky分解的最坏场景求取方法有效避免了在大量场景下进行机组组合和潮流计算的求解问题,因而计算效率提升了411.5%。
本发明按照优选实施例进行了说明,但上述实施例不以任何形式限定本发明,凡采用等同替换或等效变换的形式所获得的技术方案,均落在本发明技术方案的范围内。
Claims (9)
1.一种考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,其特征在于:包括以下步骤:
1)考虑多重随机因素的日前鲁棒调度建模;
2)最坏场景求解;
3)通过Benders分解法来对模型进行求解;
在步骤1)中,首先考虑多重随机因素的日前鲁棒调度建模;然后进行基本场景下的日前调度建模;再进行不确定场景下的日前调度建模;
在步骤2)中首先进行随机因素的概率密度函数建模,然后样本抽样,再进行正交转换矩阵推导;接着进行最坏场景求取;
在步骤3)中,首先进行UC主问题建模,然后进行基本场景下的安全子问题建模;再进行考虑不确定性的安全子问题建模;
在步骤2)中,首先利用非参数核密度估计分别构建风电、光伏以及负荷的概率密度函数,然后利用拉丁超立方抽样生成样本,最后采用Cholesky分解法将上述具有相关性的随机样本转换为相互独立的随机样本,并以此为基础确定最坏场景。
3.根据权利要求2所述的考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,其特征在于:选择高斯函数作为负荷概率密度模型的核函数,并对带宽l进行求解,从而得到系统负荷的概率密度函数φ(Pd)。
5.根据权利要求1所述的考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,其特征在于:采用拉丁超立方抽样进行分层采样。
8.根据权利要求7所述的考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,其特征在于:相关系数矩阵Cw是正定矩阵,则可以对系数矩阵进行Cholesky分解:
Cw=GGT (6)
式中:G为下三角矩阵,其中元素可由式(7)求得:
设存在一正交矩阵B,可将具有相关性的输入随机变量W转换为不相关的随机变量Y:
Y=BW (8)
由于不相关随机变量Y的相关系数矩阵CY为单位矩阵I,因而有:
CY=ρ(Y,YT)=ρ(BW,WTBT)=Bρ(W,WT)BT=BCWBT=I (9)
又由式(6)可得:
CY=BCWBT=BGGTBT=(BG)(BG)T=I (10)
由上式推导得到:
B=G-1 (11)
在已知有相关性的输入不确定量W的前提下,通过正交变换矩阵,可以将其变为不相关的随机变量Y。
9.根据权利要求8所述的考虑多重不确定性及相关性的电力系统日前鲁棒调度方法,其特征在于:通过得到的正交变换矩阵,将存在相关性的负荷、风电及光伏出力矩阵转换为独立矩阵,消除其间的相关性,进而通过最坏场景线性叠加的方法求得最坏场景。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011110581.6A CN112186765B (zh) | 2017-11-30 | 2017-11-30 | 一种机组组合决策的日前调度模型的建模方法 |
CN202011112195.0A CN112186766B (zh) | 2017-11-30 | 2017-11-30 | 一种基于Benders分解的多阶段分解方法 |
CN201711238480.5A CN107947164B (zh) | 2017-11-30 | 2017-11-30 | 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711238480.5A CN107947164B (zh) | 2017-11-30 | 2017-11-30 | 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
Related Child Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011112195.0A Division CN112186766B (zh) | 2017-11-30 | 2017-11-30 | 一种基于Benders分解的多阶段分解方法 |
CN202011110581.6A Division CN112186765B (zh) | 2017-11-30 | 2017-11-30 | 一种机组组合决策的日前调度模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107947164A CN107947164A (zh) | 2018-04-20 |
CN107947164B true CN107947164B (zh) | 2020-11-17 |
Family
ID=61947028
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011110581.6A Active CN112186765B (zh) | 2017-11-30 | 2017-11-30 | 一种机组组合决策的日前调度模型的建模方法 |
CN201711238480.5A Active CN107947164B (zh) | 2017-11-30 | 2017-11-30 | 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 |
CN202011112195.0A Active CN112186766B (zh) | 2017-11-30 | 2017-11-30 | 一种基于Benders分解的多阶段分解方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011110581.6A Active CN112186765B (zh) | 2017-11-30 | 2017-11-30 | 一种机组组合决策的日前调度模型的建模方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011112195.0A Active CN112186766B (zh) | 2017-11-30 | 2017-11-30 | 一种基于Benders分解的多阶段分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (3) | CN112186765B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108667010B (zh) * | 2018-05-04 | 2019-08-16 | 清华大学 | 一种基于分布鲁棒优化的配电网经济调度方法 |
CN109831381B (zh) * | 2018-12-15 | 2020-07-10 | 清华大学 | 一种考虑信息物理耦合的电力通信网鲁棒路由优化方法 |
CN109657385A (zh) * | 2018-12-26 | 2019-04-19 | 国网江苏省电力有限公司 | 一种适用于输电线路随机状态生成的仿真计算方法 |
CN109785183A (zh) * | 2018-12-26 | 2019-05-21 | 国网山西省电力公司电力科学研究院 | 一种考虑风电与负荷预测不确定性的鲁棒调度方法 |
CN110535120A (zh) * | 2019-08-16 | 2019-12-03 | 华南理工大学 | 考虑大气污染扩散的电-气互联系统分布式优化调度方法 |
CN111985566B (zh) * | 2019-11-21 | 2021-09-21 | 国网江苏省电力有限公司南通供电分公司 | 一种基于Copula函数的分布式电源时序联合出力典型场景生成方法 |
CN111932048B (zh) * | 2020-04-16 | 2022-08-30 | 贵州电网有限责任公司 | 基于极端场景的两阶段月度交易安全校核方法 |
CN111711184B (zh) * | 2020-05-25 | 2022-02-11 | 国网青海省电力公司 | 一种基于最坏场景辨识的电力系统鲁棒经济调度方法 |
CN112310959A (zh) * | 2020-10-13 | 2021-02-02 | 国网湖北省电力有限公司宜昌供电公司 | 考虑不确定性因素及其相关性的配电网低电压综合治理方法 |
CN112668751B (zh) * | 2020-11-26 | 2022-06-17 | 广西大学 | 一种机组优化调度模型的建立方法及装置 |
CN112989536A (zh) * | 2021-03-11 | 2021-06-18 | 南方电网科学研究院有限责任公司 | 基于场景分解的电气多能流系统优化调度方法 |
CN113642792A (zh) * | 2021-08-12 | 2021-11-12 | 中国南方电网有限责任公司 | 全面考虑复杂大电网运行约束条件的中长期机组组合精确建模技术方法 |
CN114862130B (zh) * | 2022-04-14 | 2023-03-24 | 山东大学 | 一种基于双层分解算法的多场景随机机组组合方法及系统 |
CN114759602B (zh) * | 2022-04-24 | 2024-04-05 | 国网山东省电力公司潍坊供电公司 | 考虑光伏极端场景的配电网接纳能力评估方法 |
CN115099467A (zh) * | 2022-06-01 | 2022-09-23 | 中国矿业大学 | 一种针对大气污染物扩散不确定性的电力系统分布鲁棒调度方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102832614A (zh) * | 2012-08-03 | 2012-12-19 | 国电南瑞科技股份有限公司 | 不确定性环境下发电计划的鲁棒优化方法 |
KR20150037357A (ko) * | 2013-09-30 | 2015-04-08 | 한국전력공사 | 유연 송전 시스템의 운용 장치 |
CN104809327A (zh) * | 2014-09-02 | 2015-07-29 | 长沙理工大学 | 含新能源电力调度矩不确定分布鲁棒优化方法 |
CN106655190A (zh) * | 2016-10-19 | 2017-05-10 | 云南电网有限责任公司电力科学研究院 | 一种求解风电场概率最优潮流的方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
CN107239863A (zh) * | 2017-04-12 | 2017-10-10 | 广东电网有限责任公司电力调度控制中心 | 电网安全约束的鲁棒机组组合方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8606420B2 (en) * | 2011-07-15 | 2013-12-10 | International Business Machines Corporation | Solving large-scale security-constrained economic dispatch problem in real-time |
CN102684224B (zh) * | 2012-05-25 | 2014-04-16 | 浙江大学 | 一种求解考虑风电波动性的机组组合方法 |
CN102983573B (zh) * | 2012-11-09 | 2014-10-15 | 天津大学 | 基于安全域的安全约束经济调度方法 |
CN104485665B (zh) * | 2014-12-17 | 2016-08-31 | 河海大学 | 计及风速预测误差时空相关性的动态概率潮流计算方法 |
CN104716670B (zh) * | 2015-04-01 | 2017-03-29 | 国家电网公司 | 光伏并网下基于网络安全约束的机组组合方法 |
CN106953354B (zh) * | 2017-03-10 | 2019-11-08 | 国网山东省电力公司经济技术研究院 | 考虑电压支撑的含风电电力系统机组组合方法 |
CN107039977B (zh) * | 2017-06-03 | 2020-03-31 | 广东博慎智库能源科技发展有限公司 | 电力系统鲁棒调度不确定集构建方法 |
CN107330568A (zh) * | 2017-08-16 | 2017-11-07 | 东南大学 | 基于Benders解耦的储能、分布式电源与配电网协调规划方法 |
-
2017
- 2017-11-30 CN CN202011110581.6A patent/CN112186765B/zh active Active
- 2017-11-30 CN CN201711238480.5A patent/CN107947164B/zh active Active
- 2017-11-30 CN CN202011112195.0A patent/CN112186766B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102832614A (zh) * | 2012-08-03 | 2012-12-19 | 国电南瑞科技股份有限公司 | 不确定性环境下发电计划的鲁棒优化方法 |
KR20150037357A (ko) * | 2013-09-30 | 2015-04-08 | 한국전력공사 | 유연 송전 시스템의 운용 장치 |
CN104809327A (zh) * | 2014-09-02 | 2015-07-29 | 长沙理工大学 | 含新能源电力调度矩不确定分布鲁棒优化方法 |
CN106655190A (zh) * | 2016-10-19 | 2017-05-10 | 云南电网有限责任公司电力科学研究院 | 一种求解风电场概率最优潮流的方法 |
CN107038292A (zh) * | 2017-04-01 | 2017-08-11 | 三峡大学 | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 |
CN107239863A (zh) * | 2017-04-12 | 2017-10-10 | 广东电网有限责任公司电力调度控制中心 | 电网安全约束的鲁棒机组组合方法 |
Non-Patent Citations (2)
Title |
---|
含风电场电力系统的模糊鲁棒优化调度;季峰;《中国电机工程学报》;20141005;第34卷(第28期);全文 * |
基于IGDT的有源配电网鲁棒恢复决策方法;陈珂宁;《中国电机工程学报》;20140705;第34卷(第19期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112186766A (zh) | 2021-01-05 |
CN112186765A (zh) | 2021-01-05 |
CN112186765B (zh) | 2022-06-17 |
CN107947164A (zh) | 2018-04-20 |
CN112186766B (zh) | 2022-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107947164B (zh) | 考虑多重不确定性及相关性的电力系统日前鲁棒调度方法 | |
CN107977744B (zh) | 一种基于传统Benders分解法的电力系统日前鲁棒调度方法 | |
CN112467807B (zh) | 一种多能源电力系统日前优化调度方法及系统 | |
CN110298138B (zh) | 一种综合能源系统优化方法、装置、设备及可读存储介质 | |
CN107239863B (zh) | 一种电网安全约束的鲁棒机组组合方法 | |
Fu et al. | Multiobjective stochastic economic dispatch with variable wind generation using scenario-based decomposition and asynchronous block iteration | |
Bai et al. | A two-level approach to AC optimal transmission switching with an accelerating technique | |
Yang et al. | An improved robust SCUC approach considering multiple uncertainty and correlation | |
Liu et al. | Computational study of security constrained economic dispatch with multi-stage rescheduling | |
CN107316113B (zh) | 一种输电网规划方法及系统 | |
CN109886446B (zh) | 基于改进混沌粒子群算法的电力系统动态经济调度方法 | |
CN110912205A (zh) | 基于场景集的含风电电力系统运行备用动态调度优化方法 | |
CN111313475A (zh) | 一种功率平衡约束考虑预测误差不确定变量的电力系统调度方法 | |
Huang et al. | Joint generation and reserve scheduling of wind‐solar‐pumped storage power systems under multiple uncertainties | |
CN110009141A (zh) | 基于sdae特征提取和svm分类模型的爬坡事件预测方法及系统 | |
CN115425668A (zh) | 一种基于电力系统时序生产模拟的储能容量优化配置方法 | |
CN115438964A (zh) | 一种输电网随机规划方法、设备及存储介质 | |
CN105305485A (zh) | 一种消纳大规模间歇性能源的安全约束经济调度方法 | |
Yang | Multi‐objective optimization of integrated gas–electricity energy system based on improved multi‐object cuckoo algorithm | |
Zhang et al. | Short‐Term Power Prediction of Wind Power Generation System Based on Logistic Chaos Atom Search Optimization BP Neural Network | |
Wang et al. | Optimal recovery strategy of DERs integrated distribution network based on scheduling rationality | |
CN113778681B (zh) | 基于云计算的数据处理方法、装置及存储介质 | |
CN112422520A (zh) | 一种考虑通信时滞的微电网分布式经济调度算法及装置 | |
Liu et al. | Two-stage robust optimal dispatch method considering wind power and load correlation | |
Ji et al. | A scenario probability based method to solve unit commitment of large scale energy storage system and thermal generation in high wind power penetration level system |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20180420 Assignee: Hubei Yunzhihang Drone Technology Co.,Ltd. Assignor: CHINA THREE GORGES University Contract record no.: X2023980044730 Denomination of invention: A robust day-ahead scheduling method for power systems considering multiple uncertainties and correlations Granted publication date: 20201117 License type: Common License Record date: 20231027 |
|
EE01 | Entry into force of recordation of patent licensing contract |