CN104638672A - 一种考虑变量相关性的光伏穿透功率极限的确定方法 - Google Patents
一种考虑变量相关性的光伏穿透功率极限的确定方法 Download PDFInfo
- Publication number
- CN104638672A CN104638672A CN201510039704.4A CN201510039704A CN104638672A CN 104638672 A CN104638672 A CN 104638672A CN 201510039704 A CN201510039704 A CN 201510039704A CN 104638672 A CN104638672 A CN 104638672A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- population
- load
- individual
- correlation
- 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 44
- 230000005540 biological transmission Effects 0.000 title claims abstract description 12
- 238000005070 sampling Methods 0.000 claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 230000001580 bacterial effect Effects 0.000 claims abstract description 25
- 230000002431 foraging effect Effects 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000035515 penetration Effects 0.000 claims description 28
- 238000009987 spinning Methods 0.000 claims description 15
- 238000012546 transfer Methods 0.000 claims description 12
- 230000009194 climbing Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 230000006978 adaptation Effects 0.000 claims description 6
- 238000003646 Spearman's rank correlation coefficient Methods 0.000 claims description 4
- 238000005206 flow analysis Methods 0.000 claims description 4
- 230000008929 regeneration Effects 0.000 claims description 3
- 238000011069 regeneration method Methods 0.000 claims description 3
- 238000004064 recycling Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000010248 power generation Methods 0.000 abstract description 2
- 230000032683 aging Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000001617 migratory effect Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H02J3/383—
-
- 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/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种考虑变量相关性的光伏穿透功率极限的确定方法,通过建立随机规划问题设定光伏穿透功率极限最大化为目标函数,考虑光伏电站内部、区域负荷之间以及光伏发电与节点负荷变化的相关性,基于拉丁超立方采样方法对光伏出力、传统机组出力和负荷水平进行抽样来检验约束条件,最后利用改进的细菌觅食算法求解规划问题,计算符合条件的光伏穿透功率极限。本发明方法解决了因光伏穿透功率极限不符合实际生产实践所带来的光伏入网后的安全隐患问题,能够更好地指导光伏电站的规划工作,增加了采样覆盖率,在不影响精度的情况下,减小了计算量,在提高初期寻优速度的同时也确保了后期寻优精度。
Description
技术领域
本发明属于新能源并网的规划领域,尤其涉及一种考虑变量相关性的光伏穿透功率极限的确定方法。
背景技术
近年来,由于能源危机的爆发和环境污染的日趋严重,新能源的应用越来越广泛。光伏作为新能源的一种,具有的随机性和不可控性,对电网的可靠性、稳定性、电能质量以及系统调峰特性等均产生影响。所以当光伏并入电力系统时,应当考虑该规模的光伏发电是否对电力系统的安全稳定运行造成隐患。参考风电穿透功率极限的定义,将光伏穿透功率极限定义为系统能够接受的最大光伏电场装机容量占系统最大负荷的百分比。研究电力系统的光伏穿透功率极限对光伏并网的规划与控制有着非常重要的现实意义。
确定光伏穿透功率极限的方法一般有三种:数字仿真法、数学优化法、频率约束法。然而现有的计算方法中,一般将光伏出力、节点负荷等的波动视为相互独立的随机变量,忽视了其相互影响的关系。事实上,同时接入电网的光伏电场间会相互影响,甚至其波动也会带来节点负荷的波动。此外,相同区域间的节点负荷由于身处相同的环境,其波动特性有着相同的特征。所以,现有技术中没有考虑这些随机变量间的相关特性,所确定的光伏穿透功率极限不符合实际生产实践,根据这样的光伏穿透功率极限进行光伏并网,会对电力系统的安全稳定运行造成隐患。
发明内容
发明目的:为了解决现有技术中所存在的因没有考虑光伏出力、节点负荷等之间的相关特性,导致所确定的光伏穿透功率极限不符合实际生产实践,带来光伏入网后的安全隐患问题,本发明公开了一种考虑变量相关性的光伏穿透功率极限的确定方法,通过建立随机规划模型得到光伏穿透功率极限的最优解,实现更好地指导光伏电站的规划工作。
技术方案:为了实现上述目的,本发明提供的考虑变量相关性的光伏穿透功率极限的确定方法,包括以下步骤:
(1)提取传统电网网络数据及光伏、传统机组和负荷各随机变量的相关参数,以光伏穿透功率极限最大为目标函数,将系统安全运行及稳定性水平作为约束条件,采用随机规划方法确定规划问题;
(2)初始化细菌觅食算法的各参量,随机产生一组初代种群个体的位置,所述种群个体的位置指代光伏电场的装机容量;
(3)利用拉丁超立方采样方法生成具有相关关系的光伏、负荷随机样本,并利用概率潮流分析对当前种群个体进行可行性验证,系统状态满足所述约束条件则进入步骤(4),否则,重新产生种群个体进行验证,直到达到预设的种群规模;
(4)对系统状态满足约束条件的初始种群个体进行细菌觅食算法的趋向、复制操作得到新的种群个体,采用步骤(3)对所述新的种群个体进行可行性验证,系统状态满足所述约束条件则更新种群个体的位置并进行步骤(5)中的迁徙操作,否则继续进行趋向、复制操作计算新的种群个体的位置;
(5)种群中的个体在完成指定次数的趋向、复制操作后,以预先设定的迁徙概率Ped进行迁徙操作,重新分布到求解空间;
(6)重复步骤(4)、(5),直到完成指定代数,根据细菌适应性评估函数确定符合条件的光伏穿透功率极限。
其中,步骤(1)中所述规划问题为:
目标函数:
约束条件:
Pr{-Pmax≤Pij≤Pmax}≥α1
Pr{-Δfmax≤Δf≤Δfmax}≥α2
Pr{η1ΣPL,i≤Ps,up}≥α3
Pr{η2ΣPL,i≥Ps,down}≥α4
Pr{ΣγiPG,i≥ΣδiPs,i}≥α5
其中,Ps,max为光伏电场最大装机容量,Pl,max为系统最大负荷水平,α1、α2、α3、α4、α5分别表示线路有功潮流极限、系统频率偏移、旋转备用容量、机组爬坡能力的置信水平,Pij为线路有功潮流,Pmax为线路传输的有功潮流上限,Δf为系统频率偏移,Δfmax为其上限,PL,i为负荷出力,,η1为系统上旋转备用占负荷的预设值,η2为系统下旋转备用占负荷的预设值,Ps,up、Ps,down分别为系统上、下旋转备用容量,PG,i为常规机组出力,γi为常规机组爬坡能力系数,Ps,i为光伏出力,δi为光伏功率变化系数。
其中,步骤(3)中所述具有相关关系的光伏、负荷随机样本的生成,包括以下步骤:
S1.根据光伏出力变量与节点负荷变化量的历史数据并利用Spearman秩相关系数描述所述光伏出力变量与所述节点负荷变化量之间的相关关系,计算光伏负荷秩相关矩阵Rs;
S2.利用中值拉丁超立方对光伏电站出力随机变量、传统机组出力随机变量、节点负荷变化随机变量进行采样获取光伏负荷随机样本,得到原始采样矩阵X;
S3.通过Chloskey分解光伏负荷秩相关矩阵Rs得到矩阵C,再对所述原始采样矩阵X进行转换得到具有相关关系的采样矩阵Y,表达式如下:
Y=CX;
S4.通过等概率原则将具有相关关系的采样矩阵Y转换为符合正态分布的矩阵X’;
S5.采用Gram-Schmidt序列正交化方法进行排序,通过反复利用正向、反向迭代交替来降低矩阵X’中变量间的相关性,直至完成预先设定的迭代次数,得到最终样本矩阵。
其中,步骤(3)中所述利用概率潮流分析对种群个体进行可行性验证,包括以下步骤:
根据光伏、传统机组出力和节点负荷的抽样样本,对最终样本矩阵采用直流潮流算法计算得到系统线路传输的有功功率,检验系统的线路有功潮流极限、系统频率偏移、旋转备用容量、机组爬坡能力的概率是否满足置信水平,系统状态满足所述约束条件则认为种群个体可行,不满足则认为种群个体不可行。
其中,步骤(4)中所述细菌觅食算法的趋向步长为:
其中,S为变步长,S0为初始步长,rand(0,1)产生(0,1)间的随机数,D为种群代数,D=l+k,l表示第l次趋向操作,k表示第k次复制操作。
其中,步骤(5)中进行迁徙操作,包括以下步骤:
通过细菌适应性评估函数对当前种群个体降序排列得到最优个体,若最优个体以迁徙概率Ped被选中,则不进行迁徙,重新选择被迁徙个体。
有益效果:1.本发明方法考虑了光伏并网后光伏电场内部和电力系统内部以及它们之间的相互关系,计及了光伏电场出力和节点负荷水平等随机变量间的相关性,不再将各状态变量视为互相独立的随机变量,更加符合生产实际;且光伏负荷的相关关系对电力系统起到调峰作用,有利于光伏并网,提高光伏接入电网的容量;此外,通过本发明方法指导光伏并网的规划与设计工作,调节光伏负荷的相关关系,可以尽可能的增加光伏接入容量。
2.为了得到系统的潮流分布,验证种群个体的可行性,需要对系统的各随机变量进行采样。本发明在考虑位置相近的光伏电站、节点负荷变化之间均存在一定相关关系的基础上,计及光伏电站内部、节点负荷之间、光伏电站和负荷之间的相关性等,对光伏负荷随机变量采样过程中进行相关性处理。光伏发电主要受气象条件的影响,其输出功率具有随机性和不可控性的特点,但是,对于处于相邻地区或相同气象条件的多个光伏电站来说,其输出功率之间具有较强的相关关系;类似地,处于同一地区的同一类负荷,由于受到相同环境和社会等因素的影响,具有较强的相关关系;此外,由于电网负荷与光伏发电一样,白天高夜间低,两种功率波动也具有一定的相关关系,负荷高峰期光伏系统能够较好的提供电力,起到调峰作用,有利于光伏发电并网。所以在对随机变量进行采样时,必须计及光伏负荷变化的相关关系,才能正确确定光伏并网的光伏穿透功率极限,指导优化光伏阵列的朝向和倾角,分散光伏电站的建设,来调节光伏负荷间的相关关系,有利于尽可能大的容纳光伏电源,充分利用光伏能源。
本发明采用拉丁超立方采样方法,比较于一般穿透功率极限求解问题中用到的蒙特卡罗随机模拟法,增加了采样覆盖率,在不影响精度的情况下,大大减少了计算量,提高了计算时间。
3.由于细菌觅食算法中主要是复制和迁徙行为更新,因此原理简单,容易实现,收敛速度更快。本发明采用的改进的细菌觅食算法相对于标准细菌觅食算法,首先改进了趋向操作中的步长参量,将其设定为一变量,随着种群老化而趋向减小,在提高初期寻优速度的同时也确保了后期寻优精度;另外,对迁徙行为的改进,保证了最优个体不被迁徙,确保了算法的收敛速度。
附图说明
图1为本发明方法的流程图;
图2为实施例中基于拉丁超立方采样产生相关性随机变量样本的流程图;
图3为实施例中改进的细菌算法的流程图。
具体实施方式
下面结合实施例对本发明作更进一步的说明。
图1中,利用本发明方法确定光伏穿透功率极限,包括以下步骤:
步骤一:提取传统电网网络数据及光伏、传统机组和负荷各随机变量的相关参数,采用随机规划方法确定规划问题,具体过程如下:
采用随机规划方法将光伏穿透功率极限的最大化作为目标函数:
其中,Ps,max为光伏电场最大装机容量,Pl,max为系统最大负荷水平。由于光伏穿透功率极限的确定是相对于指定的电力系统,所以系统的最大负荷水平为一常量,目标函数等价于Ps,max的最大化。将系统潮流方程作为等式约束,考虑到光伏接入主要是对系统有功功率的影响较大,所以本发明简化潮流方程式为线性关系式,采用直流潮流算法。将系统安全运行及稳定性水平作为不等式约束,预设相应的置信水平α1、α2、α3、α4、α5,主要包括以下几个约束:
1.线路有功潮流上下限约束
Pr{-Pmax≤Pij≤Pmax}≥α1
其中Pij为线路有功潮流,Pmax为线路传输的有功潮流上限。
2.系统频率偏移约束
Pr{-Δfmax≤Δf≤Δfmax}≥α2
其中Δf为系统频率偏移,Δfmax为其上限。系统的频率偏移可以通过系统的静态频率特性计算:
其中ΔP为系统有功功率不平衡量,KS系统单位调节功率,包括发电机和负荷单位调节功率。
3.旋转备用容量约束
其中PL,i为负荷出力,Ps,up、Ps,down分别为系统上、下旋转备用容量,η1为系统上旋转备用占负荷的预设值,η2为系统下旋转备用占负荷的预设值。
4.机组爬坡能力约束
Pr{ΣγiPG,i≥ΣδiPs,i}≥α5
其中PG,i为常规机组出力,γi为常规机组爬坡能力系数,Ps,i为光伏出力,δi为光伏功率变化系数。
步骤二:根据预先设定的细菌觅食算法各参量,随机产生一组初代种群个体的位置,此处的种群个体指代光伏电场,个体的位置指代光伏电场的装机容量。
步骤三:通过拉丁超立方采样产生光伏、发电机和负荷随机样本,由相关性处理方法生成具有相关关系的光伏、负荷随机样本,进行概率潮流计算,验证种群个体是否满足预设的约束条件,否则,重新产生个体进行检验,直到达到预设的种群规模。
如图1和图2所示,本发明对种群个体进行可行性验证包括以下步骤:
1.根据光伏出力变量与节点负荷变化关系的历史数据计算光伏负荷秩相关矩阵Rs;
本实施例用Spearman秩相关系数描述变量间的相互关系,Spearman秩相关系数的表达式为:
其中Rx,Ry分别为变量x,y的秩,cov(·)为随机变量的标准差,var(·)为随机变量的协方差。将光伏电站出力变量与节点负荷变化量之间的相关系数用矩阵形式排列,则得到光伏负荷秩相关矩阵Rs;
2.利用中值拉丁超立方采样获取光伏负荷随机样本,得到原始采样矩阵X,具体的采样过程如下:设
x1,x2,...,xn1,xn1+1,xn1+2,...,xn1+n2,...,xn1+n2+n3(n1+n2+n3=N)
是待求解的N个随机变量,其中(x1,x2,...,xn1)为光伏电场中每个光伏电站出力随机变量,(xn1+1,xn1+2,...,xn1+n2)传统机组出力随机变量,(xn1+n2+1,...,xn1+n2+n3)为节点负荷变化随机变量。预设K为采样数,将N个随机变量中每一个xi取值空间等分为K段,设xi的采样点为每个区间的中点,则第k个采样区间对应xi的值可以确定为:
通过上述方法对N个随机变量进行采样,可以得到N×K维采样矩阵X;
3.通过Chloskey分解光伏负荷秩相关矩阵Rs得到矩阵C将原始采样矩阵X转换成具有相关关系的采样矩阵Y,经过此步骤的处理,使得采样样本更符合实际情况,表达式如下:
Y=CX
4.通过等概率原则将具有相关关系的采样矩阵Y转换为符合正态分布的矩阵X’;
5.采用Gram-Schmidt序列正交化方法进行排序,反复利用正向、反向迭代交替,达到预先设定的迭代次数,以达到降低变量间相关性的目的,得到最终样本矩阵;
6.获取到最终样本矩阵之后,基于最终采样矩阵进行系统概率潮流分析以验证种群个体的可行性,具体如下:
光伏电池的输出功率认为服从贝塔分布,即
其中,PM为光伏电池输出功率,最大输出功率PM,max=η·A·rmax,η为光伏转换效率(这里认为是一常数),A为光伏阵列面积,α、β为贝塔分布形状参数,由光伏电站安装当地光照强度统计数据决定;对于节点负荷,本发明采用正态分布模型,均视为在期望值附近按照正态分布规律变化的随机变量;对于传统机组出力,本发明将其视为均匀分布的连续性随机变量分布,取值范围即为发电机的有功上下限;
根据光伏、传统机组出力和节点负荷的抽样样本进行系统概率潮流分析可以计算得到系统线路传输的有功功率,从而通过代入系统光伏电站出力、传统机组出力、节点负荷水平样本以及对应的线路传输有功功率,检验系统的线路有功潮流上下限、系统频率偏移、旋转备用容量、机组爬坡能力约束的置信水平,检验是否满足步骤一中的约束条件,满足则认为该种群个体可行,否则认为该种群个体不可行。
如图1和图3所示,本发明采用细菌觅食算法进行步骤一中所确定的规划问题求解,体现在以下几个步骤:
步骤四:对经过可行性验证的初始种群个体利用细菌觅食算法进行趋向、复制操作,采用步骤二中检验种群个体位置可行性的步骤对新的种群位置进行验证,满足可行性则更新种群位置并进行步骤五中的迁徙操作,否则继续进行趋向、复制操作计算新的种群位置。
本发明方法相对于标准细菌觅食算法,改进了趋向操作中的步长参量。标准细菌觅食算法的趋向步长为一常数,然而事实上,对于新生种群,步长较大有利于提高寻优速度,但对于老化群体,步长较小有利于保证寻优精度。所以本发明采用的改进细菌觅食算法中设定步长为一变量随着种群代数的增加而减小,表达式如下:
其中,S为变步长,S0为初始步长,D为种群代数,rand(0,1)产生(0,1)间的随机数。这里种群代数定义为趋向次数与复制次数的和,主要是考虑到保持种群信息的一致性:
D=l+k
其中l表示第l次趋向操作,k表示第k次复制操作。由此,在提高初期寻优速度的同时也确保了后期寻优精度。
由于算法的目的是实现光伏装机容量的最大化,所以复制操作中的细菌适应性评估函数选取为:
其中xi为光伏出力随机变量,l,k,j分别为第l次趋向操作、第k次复制操作、第j次迁徙操作的次数,每一步复制操作都根据评估函数对种群中个体进行降序排列,选择其中最优的种群个体。
步骤五:种群中个体在完成一次复制操作后,以预先设定的迁徙概率Ped进行迁徙操作,重新分布到求解空间,以提高算法搜索全局的能力。本发明方法对标准细菌觅食算法中的迁徙操作做了改进,每一次迁徙操作前通过步骤四中的细菌适应性评估函数,降序排列得到最优个体,若最优个体以迁徙概率Ped被选中,则不进行迁徙,重新选择被迁徙个体。当前种群中的最优个体具有良好的适应性,最接近全局最优解,标准细菌算法中的迁徙操作很可能使当前取得的成果遭到破坏,造成种群的退化,改进后的细菌觅食算法在保证原有算法的全局搜索能力基础上,加快了算法的收敛速度。
步骤六:重复步骤四、五,直到指定代数,确定符合条件的光伏穿透功率极限。
以上详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种等同变换,这些等同变换均属于本发明的保护范围。
Claims (6)
1.一种考虑变量相关性的光伏穿透功率极限的确定方法,其特征在于,该方法包括以下步骤:
(1)提取传统电网网络数据及光伏、传统机组和负荷各随机变量的相关参数,以光伏穿透功率极限最大为目标函数,将系统安全运行及稳定性水平作为约束条件,采用随机规划方法确定规划问题;
(2)初始化细菌觅食算法的各参量,随机产生初代种群个体的位置,所述种群个体的位置指代光伏电场的装机容量;
(3)利用拉丁超立方采样方法生成具有相关关系的光伏、负荷随机样本,并利用概率潮流分析对当前种群个体进行可行性验证,系统状态满足所述约束条件则进入步骤(4),否则,重新产生种群个体进行验证,直到达到预设的种群规模;
(4)对可行的初始种群个体进行细菌觅食算法的趋向、复制操作得到新的种群个体,采用步骤(3)对所述新的种群个体进行可行性验证,系统状态满足所述约束条件则更新种群个体的位置并进行步骤(5)中的迁徙操作,否则继续进行趋向、复制操作计算新的种群个体的位置;
(5)种群中的个体在完成指定次数的趋向、复制操作后,以预先设定的迁徙概率Ped进行迁徙操作,重新分布求解空间;
(6)重复步骤(4)、(5),直到完成指定代数,根据细菌适应性评估函数确定符合条件的光伏穿透功率极限。
2.根据权利要求1所述的考虑变量相关性的光伏穿透功率极限的确定方法,其特征在于,步骤(1)中所述规划问题为:
目标函数:
约束条件:
Pr{-Pmax≤Pij≤Pmax}≥α1
Pr{-Δfmax≤Δf≤Δfmax}≥α2
Pr{η1ΣPL,i≤Ps,up}≥α3
Pr{η2ΣPL,i≥Ps,down}≥α4
Pr{ΣγiPG,i≥ΣδiPs,i}≥α5
其中,Ps,max为光伏电场最大装机容量,Pl,max为系统最大负荷水平,α1、α2、α3、α4、α5分别表示线路有功潮流极限、系统频率偏移、旋转备用容量、机组爬坡能力的置信水 平,Pij为线路有功潮流,Pmax为线路传输的有功潮流上限,Δf为系统频率偏移,Δfmax为其上限,PL,i为负荷出力,,η1为系统上旋转备用占负荷的预设值,η2为系统下旋转备用占负荷的预设值,Ps,up、Ps,down分别为系统上、下旋转备用容量,PG,i为常规机组出力,γi为常规机组爬坡能力系数,Ps,i为光伏出力,δi为光伏功率变化系数。
3.根据权利要求1所述的考虑变量相关性的光伏穿透功率极限的确定方法,其特征在于,步骤(3)中所述具有相关关系的光伏、负荷随机样本的生成,包括以下步骤:
(1)根据光伏出力变量与节点负荷变化量的历史数据并利用Spearman秩相关系数描述所述光伏出力变量与所述节点负荷变化量之间的相关关系,计算光伏负荷秩相关矩阵Rs;
(2)利用中值拉丁超立方对光伏电站出力随机变量、传统机组出力随机变量、节点负荷变化随机变量进行采样获取光伏负荷随机样本,得到原始采样矩阵X;
(3)通过Chloskey分解光伏负荷秩相关矩阵Rs得到矩阵C,再对所述原始采样矩阵X进行转换得到具有相关关系的采样矩阵Y,表达式如下:
Y=CX;
(4)通过等概率原则将具有相关关系的采样矩阵Y转换为符合正态分布的矩阵X’;
(5)采用Gram-Schmidt序列正交化方法进行排序,通过反复利用正向、反向迭代交替来降低矩阵X’中变量间的相关性,直至完成预先设定的迭代次数,得到最终样本矩阵。
4.根据权利要求1所述的考虑变量相关性的光伏穿透功率极限的确定方法,其特征在于,步骤(3)中所述利用概率潮流分析对种群个体进行可行性验证,包括以下步骤:
根据光伏、传统机组出力和节点负荷的概率分布,对最终样本矩阵采用直流潮流算法计算得到系统线路传输的有功功率,检验系统的线路有功潮流极限、系统频率偏移、旋转备用容量、机组爬坡能力的概率是否满足置信水平,系统状态满足所述约束条件则认为种群个体可行,不满足则认为种群个体不可行。
5.根据权利要求1所述的考虑变量相关性的光伏穿透功率极限的确定方法,其特征在于,步骤(4)中所述细菌觅食算法的趋向步长为:
其中,S为变步长,S0为初始步长,rand(0,1)产生(0,1)间的随机数,D为种群代数,D=l+k,l表示第l次趋向操作,k表示第k次复制操作。
6.根据权利要求1所述的考虑变量相关性的光伏穿透功率极限的确定方法,其特 征在于,步骤(5)中进行迁徙操作,包括以下步骤:
通过所述细菌适应性评估函数对当前种群个体降序排列得到最优个体,若最优个体以迁徙概率Ped被选中,则不进行迁徙,重新选择被迁徙个体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510039704.4A CN104638672B (zh) | 2015-01-26 | 2015-01-26 | 一种考虑变量相关性的光伏穿透功率极限的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510039704.4A CN104638672B (zh) | 2015-01-26 | 2015-01-26 | 一种考虑变量相关性的光伏穿透功率极限的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104638672A true CN104638672A (zh) | 2015-05-20 |
CN104638672B CN104638672B (zh) | 2017-02-01 |
Family
ID=53217122
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510039704.4A Expired - Fee Related CN104638672B (zh) | 2015-01-26 | 2015-01-26 | 一种考虑变量相关性的光伏穿透功率极限的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104638672B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104882910A (zh) * | 2015-06-19 | 2015-09-02 | 国网四川省电力公司经济技术研究院 | 基于调峰平衡约束的风电穿透功率极限确定方法 |
CN105226730A (zh) * | 2015-10-29 | 2016-01-06 | 东北电力大学 | 一种提高风电可调度性的风储联合系统调度方法及装置 |
CN106954255A (zh) * | 2017-03-28 | 2017-07-14 | 吉林化工学院 | 基于自适应细菌觅食的认知无线电功率控制方法及系统 |
CN107681685A (zh) * | 2017-09-13 | 2018-02-09 | 国网甘肃省电力公司电力科学研究院 | 一种考虑光伏非线性相关性的电力系统概率潮流计算方法 |
CN107769645A (zh) * | 2016-08-22 | 2018-03-06 | 中车株洲电力机车研究所有限公司 | 基于电驱系统的状态估计的协同控制方法 |
CN108539781A (zh) * | 2018-03-29 | 2018-09-14 | 国网江苏省电力有限公司电力科学研究院 | 提高恢复过程安全性的扩展黑启动方案二层规划优化方法 |
CN109586280A (zh) * | 2018-11-21 | 2019-04-05 | 国网冀北电力有限公司 | 一种机组交易极限电量确定方法及系统 |
CN110555628A (zh) * | 2019-09-11 | 2019-12-10 | 国网能源研究院有限公司 | 一种新能源发电出力与负荷特性匹配程度的综合评价方法 |
CN113128071A (zh) * | 2021-05-08 | 2021-07-16 | 南京工程学院 | 一种含光伏发电的发电系统可靠性评估方法 |
CN113507134A (zh) * | 2021-06-24 | 2021-10-15 | 东北电力大学 | 一种规划目标年新能源电源装机容量的优化规划方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441535A (zh) * | 2013-08-01 | 2013-12-11 | 国电南瑞科技股份有限公司 | 基于sced的日前发电计划光伏发电接纳能力分析方法 |
-
2015
- 2015-01-26 CN CN201510039704.4A patent/CN104638672B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441535A (zh) * | 2013-08-01 | 2013-12-11 | 国电南瑞科技股份有限公司 | 基于sced的日前发电计划光伏发电接纳能力分析方法 |
Non-Patent Citations (3)
Title |
---|
于晗等: "采用拉丁超立方采样的电力系统概率潮流计算方法", 《电力系统自动化》 * |
蒋平等: "考虑风光互补的间歇性能源准入功率极限研究", 《电网技术》 * |
马溪原等: "采用改进细菌觅食算法的风/光/储混合微电网电源优化配置", 《中国电机工程学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104882910A (zh) * | 2015-06-19 | 2015-09-02 | 国网四川省电力公司经济技术研究院 | 基于调峰平衡约束的风电穿透功率极限确定方法 |
CN105226730A (zh) * | 2015-10-29 | 2016-01-06 | 东北电力大学 | 一种提高风电可调度性的风储联合系统调度方法及装置 |
CN107769645B (zh) * | 2016-08-22 | 2020-04-21 | 中车株洲电力机车研究所有限公司 | 基于电驱系统的状态估计的协同控制方法 |
CN107769645A (zh) * | 2016-08-22 | 2018-03-06 | 中车株洲电力机车研究所有限公司 | 基于电驱系统的状态估计的协同控制方法 |
CN106954255A (zh) * | 2017-03-28 | 2017-07-14 | 吉林化工学院 | 基于自适应细菌觅食的认知无线电功率控制方法及系统 |
CN107681685A (zh) * | 2017-09-13 | 2018-02-09 | 国网甘肃省电力公司电力科学研究院 | 一种考虑光伏非线性相关性的电力系统概率潮流计算方法 |
CN108539781A (zh) * | 2018-03-29 | 2018-09-14 | 国网江苏省电力有限公司电力科学研究院 | 提高恢复过程安全性的扩展黑启动方案二层规划优化方法 |
CN109586280A (zh) * | 2018-11-21 | 2019-04-05 | 国网冀北电力有限公司 | 一种机组交易极限电量确定方法及系统 |
CN110555628A (zh) * | 2019-09-11 | 2019-12-10 | 国网能源研究院有限公司 | 一种新能源发电出力与负荷特性匹配程度的综合评价方法 |
CN110555628B (zh) * | 2019-09-11 | 2022-07-22 | 国网能源研究院有限公司 | 一种新能源发电出力与负荷特性匹配程度的综合评价方法 |
CN113128071A (zh) * | 2021-05-08 | 2021-07-16 | 南京工程学院 | 一种含光伏发电的发电系统可靠性评估方法 |
CN113128071B (zh) * | 2021-05-08 | 2024-02-09 | 南京工程学院 | 一种含光伏发电的发电系统可靠性评估方法 |
CN113507134A (zh) * | 2021-06-24 | 2021-10-15 | 东北电力大学 | 一种规划目标年新能源电源装机容量的优化规划方法 |
CN113507134B (zh) * | 2021-06-24 | 2022-04-05 | 东北电力大学 | 一种规划目标年新能源电源装机容量的优化规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104638672B (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104638672A (zh) | 一种考虑变量相关性的光伏穿透功率极限的确定方法 | |
Lin et al. | An improved moth-flame optimization algorithm for support vector machine prediction of photovoltaic power generation | |
CN103390116B (zh) | 采用分步方式的光伏电站发电功率预测方法 | |
Kardakos et al. | Application of time series and artificial neural network models in short-term forecasting of PV power generation | |
CN110474339B (zh) | 一种基于深度发电负荷预测的电网无功控制方法 | |
CN105184678A (zh) | 基于多种神经网络组合算法的光伏电站发电量短期预测模型的构建方法 | |
CN103106544B (zh) | 一种基于t-s型模糊神经网络的光伏发电预测系统 | |
CN113496311A (zh) | 光伏电站发电功率预测方法及系统 | |
CN104009494B (zh) | 一种环境经济发电调度方法 | |
CN105449713A (zh) | 考虑分布式电源特性的有源配电网智能软开关规划方法 | |
CN107480825B (zh) | 一种计及容量可信度的光伏电站优化规划方法 | |
CN107069814A (zh) | 配网分布式电源容量布点的模糊机会约束规划方法与系统 | |
CN110048456A (zh) | 一种基于大规模风电接入输电网的源网联合规划方法 | |
CN103441535A (zh) | 基于sced的日前发电计划光伏发电接纳能力分析方法 | |
CN103986193B (zh) | 一种最大风电并网容量获取的方法 | |
CN112149264A (zh) | 一种基于改进哈里斯鹰优化算法的主动配电网规划方法 | |
CN104218578A (zh) | 一种分布式供电系统的规划方法及装置 | |
CN103473393A (zh) | 一种考虑随机概率的输电裕度控制模型建模方法 | |
CN113191091B (zh) | 基于混合深度学习机制的风速预测方法及系统、设备 | |
CN110445127A (zh) | 一种面向多重随机不确定性的配电网无功优化方法及系统 | |
CN110084430A (zh) | 一种考虑时空特性设计分布式光伏出力预测模型的方法 | |
CN114595891A (zh) | 一种配电网电压及潮流越界风险评估方法、系统及设备 | |
CN104915788B (zh) | 一种考虑多风场相关性的电力系统动态经济调度的方法 | |
CN106026200A (zh) | 风电场的电力系统无功优化方法 | |
Yang et al. | Optimal placement and sizing of distributed generators based on multiobjective particle swarm optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170201 Termination date: 20210126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |