CN108306303B - 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 - Google Patents
一种考虑负荷增长和新能源出力随机的电压稳定评估方法 Download PDFInfo
- Publication number
- CN108306303B CN108306303B CN201810045708.7A CN201810045708A CN108306303B CN 108306303 B CN108306303 B CN 108306303B CN 201810045708 A CN201810045708 A CN 201810045708A CN 108306303 B CN108306303 B CN 108306303B
- Authority
- CN
- China
- Prior art keywords
- load
- output
- node
- power
- active
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 30
- 230000012010 growth Effects 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000009826 distribution Methods 0.000 claims abstract description 38
- 238000005070 sampling Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000006641 stabilisation Effects 0.000 claims abstract description 15
- 238000011105 stabilization Methods 0.000 claims abstract description 15
- 238000013097 stability assessment Methods 0.000 claims abstract description 10
- 238000003064 k means clustering Methods 0.000 claims abstract description 8
- 238000010248 power generation Methods 0.000 claims description 36
- 230000005684 electric field Effects 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 7
- 238000005286 illumination Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 4
- 230000005283 ground state Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000005315 distribution function Methods 0.000 claims description 2
- 230000005855 radiation Effects 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims 1
- 230000000087 stabilizing effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 8
- 238000012937 correction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000007773 growth pattern Effects 0.000 description 1
- 238000002789 length control Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- 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
-
- 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
- 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
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- 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
- 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
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑负荷增长和新能源出力随机的电压稳定评估方法,属于电力系统稳定分析与控制的技术领域。本发明考虑负荷增长和新能源出力随机对电压稳定评估的影响,提出了一种电压稳定评估方法,将不确定问题转化为若干确定性问题求解。该方法基于历史负荷数据采用改进K均值聚类算法将负荷节点分类,根据负荷预测结果定义负荷随机增长的期望。根据负荷增长方向、风电和光伏出力的概率分布应用拉丁超立方采样获得负荷增长方向和电源出力样本,在此基础上运用连续潮流对各样本进行确定性的电压稳定负荷裕度计算,进而获得负荷裕度的统计特征,得到概率评估结果。
Description
技术领域
本发明涉及一种考虑负荷增长和新能源出力随机的电压稳定评估方法,属于电力系统稳定分析与控制的技术领域。
背景技术
目前基于确定性的电压稳定评估体系已经趋于成熟,但电网在实际运行中存在诸多随机因素,特别是近些年来风电、光伏等新能源机组并网,使得确定性研究无法满足系统随机变化的要求。因此,研究计及随机因素的电压稳定评估方法对系统的安全稳定运行具有重要意义。
计及负荷变化、电源出力和元件故障的电压稳定评估方法得到了充分研究。研究计及负荷增长和新能源出力随机的电压稳定概率评估具有重要工程实际意义。目前对于负荷增长不确定性的研究主要问题是难以定义负荷增长的期望,现有的定义方法缺乏实际意义。文献一《基于负荷不确定性建模的静态电压稳定性风险评估》(中国电机工程学报,2016年第36卷第13期第3471页)考虑负荷增长的随机性,用正态分布来模拟随机增长,通过将增负荷节点分类来考虑各类节点增长方式的不同,但各节点随机增长的期望取该节点的基态负荷值缺乏实际意义。文献二《考虑不均衡区域负荷增长的在线电压稳定评估》(电力自动化设备,2011年第31卷第3期第58页)在确定性电压稳定评估中负荷增长量的给定是根据下一时段负荷预测结果与当前负荷水平的差值,给负荷随机增长期望值的给定提供参考。电源出力的随机性主要考虑大规模新能源机组出力的间歇性和不确定性,主要有风电和光伏出力,需要根据评估地区的气象因素、时间因素等,其出力随风速和太阳辐照度的变化而变化。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足,提出了一种考虑负荷增长和新能源出力随机的电压稳定评估方法,该方法定义负荷增长方向和新能源出力随机时,计及了历史负荷数据、负荷预测结果以及风电场、光伏电场出力的不确定性,引入概率分析方法,使得电压稳定评估结果更具有参考意义。
本发明采用如下技术方案实现:本发明考虑负荷增长和新能源出力随机的电压稳定评估方法,包括以下具体步骤:
1)建立负荷增长和新能源出力随机的电压稳定负荷裕度计算模型;
2)根据评估意愿和先验知识指定负荷和发电增长范围,获取增长范围内各节点相关参数;
3)对指定负荷增长范围内负荷节点的历史负荷数据采用改进K均值聚类方法进行分类,得到各负荷节点群,根据负荷分类和负荷预测计算负荷增长范围内节点群负荷增长方向,基于负荷增长方向和新能源出力的随机性分别建立概率分布模型;
4)根据各概率分布通过拉丁超立方采样技术得到负荷增长方向和新能源出力采样样本;
5)根据各样本的负荷增长方向和新能源出力样本,利用连续潮流计算有功负荷裕度。直至计算完所有采样样本,得到各样本下的负荷裕度值,求取负荷裕度的统计特征,分析电压稳定概率评估结果。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)对负荷节点进行聚类,考虑不同类型负荷增长方式的不同,计及同一类负荷变化的相关性,不同类型负荷变化的独立性,同时将各负荷节点群增长系数列为随机变量,降低了随机变量的维度,从而减少了计算量;
(2)基于负荷预测结果定义负荷增长随机变化的期望,相较于其它方法更能真实反映电力系统实际负荷增长方式。
(3)与确定性电压稳定评估方法相比,本发明计及了负荷增长和新能源出力随机性,用概率分布来体现不确定性,概率评估结果有助于系统运行人员全面了解电网电压稳定负荷裕度的概率分布情况,对电网运行方式的制定和电力系统调度有一定的指导意义。
附图说明
图1为本发明方法的流程图;
图2为拉丁超立方采样示意图;
图3为连续潮流的计算过程示意图;
图4为IEEE118节点系统接线图;
图5为区域B负荷-发电增长模式下的负荷裕度概率密度分布示意图;
图6为区域A和区域B负荷-发电增长模式下的负荷裕度概率密度分布示意图;
图7为区域A发电增长区域B负荷增长模式下的负荷裕度概率密度分布示意图;
图8为区域B发电增长区域A负荷增长模式下的负荷裕度概率密度分布示意图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明公开了一种考虑负荷增长和新能源出力随机的电压稳定评估方法,属于电力系统稳定分析与控制的技术领域。本发明考虑负荷增长和新能源出力随机对电压稳定评估的影响,提出了一种随机电压稳定负荷裕度计算方法,将不确定问题转化为确定性问题求解。该方法考虑了负荷增长和新能源出力的随机性,根据历史负荷数据采用改进K均值聚类算法将负荷节点分类,基于负荷预测结果定义负荷随机增长的期望,根据负荷增长方向、风电和光伏出力的概率分布应用拉丁超立方采样获得负荷增长方向和电源出力样本,在此基础上运用连续潮流对各样本进行确定性的电压稳定负荷裕度计算,进而获得负荷裕度的统计特征。
本发明所涉及的考虑负荷增长和新能源出力随机的电压稳定评估方法计算流程图如图1所示。具体按照以下步骤:
1)建立负荷增长和新能源出力随机的电压稳定负荷裕度计算模型。
电压稳定负荷裕度计算通常转化为电压稳定临界点的计算,以负荷裕度为目标函数,考虑负荷增长方向向量和新能源出力的随机性以及系统静态安全约束,构成如下随机电压稳定负荷裕度计算模型:
式中:λ为负荷参数,标量,无单位;PLi0、QLi0分别为节点i负荷的有功和无功功率;PGi0、QGi0和分别为节点i常规发电机有功、无功出力和新能源发电机有功、无功出力,其中上标“~”表示随机变量;PGi,max、PGi,min和QGi,max、QGi,min分别为节点i常规发电机有功出力上、下限和无功出力上、下限;PRi,max、PRi,min和QRi,max、QRi,min分别为节点i新能源发电机有功出力上、下限和无功出力上、下限;和kGi分别为节点i的有功、无功负荷增长系数和常规发电机有功出力增长系数;Vi、Vj分别为节点i与j的电压幅值;Gij、Bij分别为导纳矩阵中第i行第j列元素的实部和虚部;θij为节点i与j之间的相角差。ΩN为所有节点集合;ΩG为常规发电机节点集合;ΩR为新能源发电机节点集合;ΩL为负荷节点集合。
由于各节点的负荷增长系数kPi、kQi和新能源出力PRi、QRi为随机变量,则计算得到的负荷裕度也是一个随机量,其概率分布特征与kPi、kQi和PRi、QRi的随机分布特征密切相关。
2)根据评估意愿和先验知识指定负荷和发电增长范围,获取增长范围内各节点相关参数。
在进行电压稳定评估之前,系统运行人员需要根据评估意愿指定负荷和发电增长的范围,如按区域或地区增长,并得到各负荷节点和发电机节点的相关参数信息,如负荷节点的历史负荷和负荷预测值、发电机节点出力上下限等等。启动计算的系统规划和运行人员具有电网运行特征的先验知识,如规划人员了解未来一段时间哪些地区规划有新的电源点,其发电会增长,哪些地区规划有新的负荷点,其负荷会大量增长;而调度人员了解短期内哪些分区发电有增加的裕度,哪些分区负荷会增长。
3)对指定负荷增长范围内负荷节点的历史负荷数据采用改进K均值聚类方法进行分类,得到各负荷节点群,根据负荷分类和负荷预测计算负荷增长范围内节点群负荷增长方向,基于负荷增长方向和新能源出力的随机性分别建立概率分布模型。
31)对于系统运行人员指定的负荷增长范围内的负荷,其增长模式也是与其负荷类型紧密相关的,在负荷增长方向的定义中考虑其分类,使得负荷裕度计算结果更加具有实际意义。本发明基于各负荷节点历史负荷数据采用改进K均值聚类算法将负荷变化相似的节点归到同一类,实现负荷分类的改进K均值聚类算法具体步骤如下:
311)已知N个负荷节点的历史负荷数据(1天24点),得到各负荷样本PLi=[PLi,1,PLi,2…PLi,24],(i=1,2…N),指定聚类数目K;
312)先对各负荷样本进行一次层次聚类,根据层次聚类的结果得到各类的初始均值向量uj=[uj,1,uj,2…uj,24],(j=1,2…K);
313)计算各负荷样本PLi与各均值向量uj的欧氏距离dist(PLi,uj):
将距离某个均值向量最小的样本划入该均值向量所在的类,对所有样本划分后得到K个负荷节点群Cj={Pj,1,Pj,2…Pj,nj},其中nj为第j个节点群的负荷节点个数;
314)根据步骤313)的聚类结果更新各节点群均值向量uj:
315)重复步骤313)和步骤314),直到聚类结果不再变化时停止计算,得到负荷聚类结果。
32)得到各负荷节点群后,根据各节点群负荷预测结果,定义第j类节点群第t时刻的负荷增长系数ηj,t为:
式中:Pj,k,tmax表示第j类负荷中的第k个负荷节点日尖峰时刻tmax的有功负荷预测值;Pj,k,t表示第j类负荷中的第k个负荷节点第t时刻的实际有功负荷值。
33)由于负荷聚类是依据负荷变化的相似程度,定义负荷增长方向时,假定不同节点群之间负荷增长相互独立,同一节点群内各节点有功负荷增长系数与节点群增长系数相同,即群内各负荷同调增长,且各节点无功负荷按恒定功率因数比例增长,负荷增长带来的网损增量由全网的松弛节点来承担,则对于第j个节点群内的第i个负荷节点的增长系数有:
kPi,j=kQi,j=γj
式中:kPi,j和kQi,j分别表示第j个节点群内的第i个节点的有功增长系数和无功增长系数。
34)定义常规发电机组的发电增长方式为各发电机按有功备用容量比例分配有功增长功率,其中有功备用容量PRES,i定义为
PRES,i=PGi,max-PGi0
则发电机i出力的有功增长方向为
该种增长方式下各发电机出力同时达到上限,此时发电增长范围内无可用的有功储备,为计算至电压稳定临界点,功率缺额由指定发电增长范围外的发电机以某种方式承担。
35)建立负荷增长概率分布模型
36)建立风电场出力概率模型
对风电场而言,其输出功率取决于该地区的风速,风速的随机变化近似服从双参数的威布尔分布,其概率密度函数为:
式中:v为风速;kw为威布尔分布的形状参数;cw为尺度参数。
单个风机有功输出功率PWg随风速v的变化用如下函数表示:
式中:vci为切入风速;vr为额定风速;vco为切出风速;Pr为风机的额定功率;k1和k2为常数,其中k1=Pr/(vr-vci),k2=-k1vci。
经相关统计,大部分时间内风速基本维持在vci与vr之间,则单个风机出力的概率密度函数为
风电场输出功率PWi为:
PWi=NWiPWg
式中:NWi为风电场的风机台数。
风电场运行方式按恒功率因数运行,则无功功率出力QWi为
37)建立光伏发电系统出力概率模型
对于光伏发电系统,光照辐照度的变化在一定时间段内近似满足贝塔分布,其概率密度函数为
式中:r太阳光照辐照度;rmax为最大辐射度;α和β为贝塔分布形状参数。
光伏电场的输出有功功率PSi与光照辐照度r的关系用如下函数表示:
PSi=rAη
式中:A为太阳能方阵的总面积;η为光电转换效率。
则光伏电场出力的概率密度函数为
式中:Pmax=rmaxAη。
通常光伏发电系统只向电网提供有功功率,本发明不考虑光伏发电系统的无功出力。
38)新能源机组并网后在潮流计算中,将每一个风电场或光伏电场等值为一个PQ节点,其有功和无功出力分别为PRi和QRi。计及了新能源机组的随机出力后,将PRi和QRi分别列为随机变量,则的新能源机组随机出力可表示如下:
4)根据各概率分布通过拉丁超立方采样技术得到负荷增长方向和新能源出力采样样本。
根据步骤35)、步骤36)和步骤37)得到的各个节点群的负荷增长系数和新能源出力概率分布,利用拉丁超立方采样方法得到负荷增长方向和新能源出力样本集,其步骤如下:
41)假设系统有N个随机变量,输入随机变量X为各负荷节点群的有功负荷增长方向、各风电场和光伏电场的有功出力,设为X=[x1,x2,…,xN],采样规模为M,其中一个随机变量xk(k=1,2,…,N)的概率分布函数为Yk=Fk(xk),值域为[0,1]。采样时将值域分为M个等间距不重叠的子区间,选取子区间的中点作为采样值,则随机变量xk的第i(i=1,2,…,M)个采样值为其中是Fk(·)的反函数。采样过程如图2所示。
42)经过步骤41)采样过程得到初始采样矩阵XN×M相关性较高,用cholesky分解法进行排序,通过改变采样值的排列顺序来减小相关性。
43)经过步骤41)和步骤42)两个步骤得到满足各随机变量概率分布且相关性较低的负荷方向和新能源出力样本集。
5)根据各样本的负荷增长方向和新能源出力样本,利用连续潮流计算有功负荷裕度。直至计算完所有采样样本,得到各样本下的负荷裕度值,求取负荷裕度的统计特征,分析电压稳定概率评估结果。
51)针对每一个采样样本中的新能源出力PRi、QRi以及负荷增长系数kPi、kQi采用连续潮流求解电压稳定负荷裕度值。考虑新能源出力以及负荷增长和发电增长的参数化连续潮流方程为:
52)重复步骤51),直至计算完所有负荷方向和新能源出力样本集,得到不同样本下各个负荷裕度值。
53)根据各负荷裕度值求取负荷裕度统计特征,设LM(Xi)为第i个样本Xi对应的负荷裕度,M为样本总数,P0为区域基态有功负荷,本发明考虑的负荷裕度统计特征如下:
4)最大值LMmin与最小值LMmax及其各自对应的负荷增长方向Dmin,Dmax。
5)概率密度分布信息。
连续潮流通过在常规潮流方程中引入负荷变化参数并增加一维扩展方程,采用预测-校正的方法解决常规潮流方程在临界点附近雅克比矩阵奇异的问题。连续潮流主要过程包括预测、校正、参数化策略和步长控制四部分,其计算示意图如图3所示。
下面结合图4至图8,以IEEE 118节点系统为具体实例,说明本发明所属模型和方法的可行性和有效性。
IEEE118节点系统接线图如图4所示,将整个系统分为两个区域,分别为区域A和区域B。风电场和光伏发电系统数据如表1和表2所示。在IEEE118节点系统中,节点23、39和114分别接入风电场,节点44和118分别接入光伏发电系统。
表1风电场相关参数
风电场 | N<sub>w</sub> | P<sub>r</sub>/MW | v<sub>ci</sub>/(m/s) | v<sub>r</sub>/(m/s) | k<sub>w</sub> | c<sub>w</sub> |
1 | 100 | 0.75 | 4.0 | 15.0 | 1.4 | 6.0 |
2 | 50 | 1.50 | 3.0 | 14.0 | 1.8 | 7.0 |
3 | 40 | 2.00 | 3.0 | 14.0 | 1.6 | 6.5 |
表2光伏发电系统相关参数
光伏电站 | A | η | r<sub>max</sub> | α | β |
1 | 2000 | 14 | 700 | 0.95 | 0.95 |
2 | 1800 | 14 | 700 | 0.90 | 0.90 |
本实例模拟两种负荷-发电增长模式:1)区域负荷-发电增长;2)区域间负荷-发电增长。将指定负荷增长区域内的负荷分为3类,基于各节点群的负荷预测峰值定义负荷增长系数随机变化的期望,为表明所提算法的有效性,定义两种运行场景:场景1:确定性评估。指定负荷增长区域内所有负荷确定性增长,功率因数保持恒定,各节点群增长系数取期望值,不计及新能源机组的出力。场景2:本发明所提概率性评估方法。正态分布标准差取期望值的5%,拉丁超立方采样次数取200。分别计算两种运行场景下的各负荷-发电增长模式得到负荷裕度的统计特征和系统薄弱母线信息(负荷裕度单位:MW)。
(1)区域负荷-发电增长
表3区域负荷-发电增长模式统计特征
区域B负荷-发电增长和区域A和区域B负荷-发电增长模式下的负荷裕度概率密度分布分别如图5和图6所示。
(2)各区域间负荷-发电增长
表4各区域间负荷-发电增长统计特征
区域A发电增长区域B负荷增长和区域A发电增长区域B负荷增长模式下的负荷裕度概率密度分布分别如图7和图8所示。
从以上各计算结果可知,考虑负荷增长和新能源出力的电压稳定概率评估相较于电压稳定确定性评估的结果更全面,后者是前者取样样本中的一个特例,且前者可给出各负荷增长方向和各新能源出力水平下的各负荷裕度值和负荷裕度的整体概率分布信息。综合看来,考虑了负荷增长和新能源出力随机的电压稳定概率评估结果能给系统的安全稳定运行提供有价值的参考。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种考虑负荷增长和新能源出力随机的电压稳定评估方法,其特征在于,包括以下具体步骤:
1)建立负荷增长和新能源出力随机的电压稳定负荷裕度计算模型;
2)根据评估要求和先验知识指定负荷和发电增长范围,获取增长范围内各节点相关参数;
3)对指定负荷增长范围内负荷节点的历史负荷数据采用改进K均值聚类方法进行分类,得到各负荷节点群,根据负荷分类和负荷预测计算负荷增长范围内节点群负荷增长方向,基于负荷增长方向和新能源出力的随机性分别建立概率分布模型;
4)根据各概率分布通过拉丁超立方采样技术得到负荷增长方向和新能源出力采样样本;
5)根据各样本的负荷增长方向和新能源出力样本,利用连续潮流计算有功负荷裕度,直至计算完所有采样样本,得到各样本下的负荷裕度值,求取负荷裕度的统计特征,分析电压稳定概率评估结果;
所述步骤3)包括如下过程:
31)对于系统运行人员指定的负荷增长范围内的负荷,其增长模式也是与其负荷类型紧密相关的,在负荷增长方向的定义中考虑其分类,使得负荷裕度计算结果更加具有实际意义;基于各负荷节点历史负荷数据采用改进K均值聚类算法将负荷变化相似的节点归到同一类,得到各负荷节点群;
32)得到各负荷节点群后,再根据各节点群负荷预测结果,定义第j类节点群的有功负荷增长系数γj取值为ηj,t:
式中:ηj,t为第j类节点群第t时刻的增长系数;Pj,k,tmax表示第j类负荷中的第k个负荷节点日尖峰时刻tmax的有功负荷预测值;Pj,k,t表示第j类负荷中的第k个负荷节点第t时刻的实际有功负荷值;
33)由于负荷聚类是依据负荷变化的相似程度,定义负荷增长方向时,假定不同节点群之间负荷增长相互独立,同一节点群内各节点有功负荷增长系数与节点群增长系数相同,即群内各负荷同调增长,且各节点无功负荷按恒定功率因数比例增长,负荷增长带来的网损增量由全网的松弛节点来承担,则对于第j个节点群内的第i个负荷节点的增长系数有:
kPi,j=kQi,j=γj
式中:kPi,j和kQi,j分别表示第j个节点群内的第i个节点的有功增长系数和无功增长系数;
34)定义常规发电机组的发电增长方式为各发电机按有功备用容量比例分配有功增长功率,其中有功备用容量PRES,i为
PRES,i=PGi,max-PGi0
则发电机i出力的有功增长方向为
该增长方向下各发电机出力同时达到上限,此时发电增长范围内无可用的有功储备,为计算至电压稳定临界点,功率缺额由指定发电增长范围外的发电机以某种方式承担;
35)建立负荷增长概率分布模型
36)建立风电场出力概率模型
对风电场而言,其输出功率取决于风电场的风速,风速的随机变化近似服从双参数的威布尔分布,其概率密度函数为:
式中:v为风速;kw为威布尔分布的形状参数;cw为尺度参数;
单个风机有功输出功率PWg随风速v的变化用如下函数表示:
式中:vci为切入风速;vr为额定风速;vco为切出风速;Pr为风机的额定功率;k1和k2为常数,其中k1=Pr/(vr-vci),k2=-k1vci;
经相关统计,大部分时间内风速基本维持在vci与vr之间,则单个风机出力的概率密度函数为
风电场输出功率PWi为
PWi=NWiPWg
式中:NWi为风电场的风机台数;
风电场运行方式按恒功率因数运行,则无功功率出力QWi为
37)建立光伏发电系统出力概率模型
对于光伏发电系统,光照辐照度的变化在一定时间段内近似满足贝塔分布,其概率密度函数为
式中:r太阳光照辐照度;rmax为最大辐射度;α和β为贝塔分布形状参数;
光伏电场的输出有功功率PSi与光照辐照度r的关系用如下函数表示:
PSi=rAη
式中:A为太阳能方阵的总面积;η为光电转换效率;
则光伏电场出力的概率密度函数为
式中:Pmax=rmaxAη;
通常光伏发电系统只向电网提供有功功率,不考虑光伏发电系统的无功出力;
38)新能源机组并网后在潮流计算中,将每一个风电场或光伏电场等值为一个PQ节点,其有功和无功出力分别为PRi和QRi,计及了新能源机组的随机出力后,将PRi和QRi分别列为随机变量,则的新能源机组随机出力可表示如下:
2.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法,其特征在于,所述步骤1)包括如下内容:
电压稳定负荷裕度计算通常转化为电压稳定临界点的计算,以负荷裕度为目标函数,考虑负荷增长方向向量和新能源出力的随机性以及系统静态安全约束,构成如下随机电压稳定负荷裕度计算模型:
式中:λ为负荷参数,标量,无单位;PLi0、QLi0分别为节点i负荷的有功和无功功率;PGi0、QGi0和分别为节点i常规发电机有功、无功出力和新能源发电机有功、无功出力,其中上标“~”表示随机变量;PGi,max、PGi,min和QGi,max、QGi,min分别为节点i常规发电机有功出力上、下限和无功出力上、下限;PRi,max、PRi,min和QRi,max、QRi,min分别为节点i新能源发电机有功出力上、下限和无功出力上、下限;和kGi分别为节点i的有功、无功负荷增长系数和常规发电机有功出力增长系数;Vi、Vj分别为节点i与j的电压幅值;Gij、Bij分别为导纳矩阵中第i行第j列元素的实部和虚部;θij为节点i与j之间的相角差,ΩN为所有节点集合;ΩG为常规发电机节点集合;ΩR为新能源发电机节点集合;ΩL为负荷节点集合。
3.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法,其特征在于,所述步骤2)包括如下内容:
系统运行人员进行电压稳定评估时,需要根据评估意愿指定负荷和发电增长范围,按区域或地区增长,并得到增长范围内各负荷节点和发电机节点的相关参数信息,启动计算的系统运行人员具有电网运行特征的先验知识,规划人员了解未来一段时间哪些地区规划有新的电源点,其发电会大量增长,哪些地区规划有新的负荷点,其负荷会大量增长;而调度人员了解短期内哪些分区发电有增加的裕度,哪些分区负荷会增长。
4.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法,其特征在于,所述步骤4)的拉丁超立方采样方法包括如下过程:
41)假设系统有N个随机变量,输入随机变量X为各负荷节点群的有功负荷增长方向、各风电场和光伏电场的有功出力,设为X=[x1,x2,…,xN],采样规模为M,其中一个随机变量xk的概率分布函数为Yk=Fk(xk),值域为[0,1],采样时将值域分为M个等间距不重叠的子区间,选取子区间的中点作为采样值,则随机变量xk的第i个采样值为其中Fk -1(·)是Fk(·)的反函数,k=1,2,L,N,i=1,2,L,M;
42)经过步骤41)采样过程得到初始采样矩阵XN×M相关性较高,用cholesky分解法进行排序,通过改变采样值的排列顺序来减小相关性;
43)经过41)和42)排序两个步骤得到满足各随机变量概率分布且相关性较低的负荷方向和新能源出力样本集。
5.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法,其特征在于,所述步骤5)包括如下过程:
51)针对每一个采样样本中的新能源出力PRi和QRi和负荷增长系数kPi和kQi,采用连续潮流求解电压稳定负荷裕度值;
52)重复步骤51),直至计算完所有采样样本,得到不同样本下各个负荷裕度值;
53)根据各负荷裕度值求取负荷裕度统计特征,设LM(Xi)为第i个样本Xi对应的负荷裕度,M为样本总数,P0为区域基态有功负荷,考虑的负荷裕度统计特征如下:
4)最大值LMmin与最小值LMmax及其各自对应的负荷增长方向Dmin,Dmax;
5)概率密度分布信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810045708.7A CN108306303B (zh) | 2018-01-17 | 2018-01-17 | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810045708.7A CN108306303B (zh) | 2018-01-17 | 2018-01-17 | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108306303A CN108306303A (zh) | 2018-07-20 |
CN108306303B true CN108306303B (zh) | 2021-06-15 |
Family
ID=62865800
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810045708.7A Active CN108306303B (zh) | 2018-01-17 | 2018-01-17 | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108306303B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109450000B (zh) * | 2017-11-09 | 2021-07-30 | 广东电网有限责任公司电力调度控制中心 | 一种基于负荷率调整方向的发电计划偏差电量分配方法 |
CN108920868A (zh) * | 2018-07-23 | 2018-11-30 | 国网江西省电力有限公司电力科学研究院 | 基于智能电表数据的负荷辨识模型参数修正方法 |
CN109119987B (zh) * | 2018-09-05 | 2022-11-01 | 国网天津市电力公司 | 基于电网灵活性的城市电网不平衡功率波动范围评估方法 |
CN109617083A (zh) * | 2018-12-26 | 2019-04-12 | 贵州电网有限责任公司 | 基于粒子群算法的配电网故障修复方法 |
CN110096730B (zh) * | 2019-03-15 | 2023-04-07 | 国网辽宁省电力有限公司 | 一种电网电压快速评估方法及系统 |
CN110460085B (zh) * | 2019-08-28 | 2023-04-07 | 华北电力大学(保定) | 一种考虑风电和负荷特性对电力系统影响的方法 |
CN110531271A (zh) * | 2019-10-14 | 2019-12-03 | 深圳市蓝度汽车电控技术有限公司 | 一种汽车电瓶静态电压的检测方法、系统及终端 |
CN110932277A (zh) * | 2019-12-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | 一种电网静态电压稳定负荷裕度分析方法、装置和设备 |
CN111463796B (zh) * | 2020-02-28 | 2021-10-08 | 广东电网有限责任公司 | 一种考虑风电出力的电压稳定性评估方法及装置 |
CN112564181B (zh) * | 2020-11-25 | 2023-05-26 | 云南电网有限责任公司 | 考虑新能源出力不确定性的在线概率安全监测方法 |
CN112366713B (zh) * | 2020-11-30 | 2023-02-28 | 南方电网科学研究院有限责任公司 | 交直流混联电网静态电压稳定计算方法、装置及存储介质 |
CN113328467B (zh) * | 2021-05-26 | 2023-01-20 | 广东电网有限责任公司 | 一种概率电压稳定性评估方法、系统、终端设备及介质 |
CN113452028B (zh) * | 2021-07-29 | 2022-05-10 | 广东电网有限责任公司 | 低压配电网概率潮流计算方法、系统、终端和存储介质 |
CN113887809A (zh) * | 2021-10-11 | 2022-01-04 | 国网新疆电力有限公司巴州供电公司 | 双碳目标下配电网供需平衡方法、系统、介质及计算设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318489A (zh) * | 2014-10-21 | 2015-01-28 | 广东电网有限责任公司电力科学研究院 | 一种基于负载特性分析的变压器分群方法 |
CN105305488A (zh) * | 2015-10-13 | 2016-02-03 | 国家电网公司 | 一种考虑新能源并网对输电网利用率影响的评价方法 |
CN105654245A (zh) * | 2015-12-31 | 2016-06-08 | 太原理工大学 | 基于负荷不确定性建模的电网静态电压稳定风险评估方法 |
CN106712080A (zh) * | 2016-12-29 | 2017-05-24 | 中国电力科学研究院 | 一种基于表计配置的有源配电网可观测性分析方法及系统 |
-
2018
- 2018-01-17 CN CN201810045708.7A patent/CN108306303B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104318489A (zh) * | 2014-10-21 | 2015-01-28 | 广东电网有限责任公司电力科学研究院 | 一种基于负载特性分析的变压器分群方法 |
CN105305488A (zh) * | 2015-10-13 | 2016-02-03 | 国家电网公司 | 一种考虑新能源并网对输电网利用率影响的评价方法 |
CN105654245A (zh) * | 2015-12-31 | 2016-06-08 | 太原理工大学 | 基于负荷不确定性建模的电网静态电压稳定风险评估方法 |
CN106712080A (zh) * | 2016-12-29 | 2017-05-24 | 中国电力科学研究院 | 一种基于表计配置的有源配电网可观测性分析方法及系统 |
Non-Patent Citations (3)
Title |
---|
Probabilistic load flow evaluation with hybrid Latin hypercube sampling and Cholesky decomposition;H. Yu等;《IEEE》;20090531;第24卷(第2期);第661页-第667页 * |
基于进化算法改进拉丁超立方抽样的概率潮流计算;李俊芳等;《中国电机工程学报》;20110905;第31卷(第25期);第90页-第96页 * |
实际负荷变化模式下的静态电压稳定分析;伍利等;《电力系统及其自动化学报》;20140531;第26卷(第5期);第50页-第54页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108306303A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108306303B (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
CN110048420B (zh) | 一种配电网随机优化调度的方法、装置和介质 | |
Ganguly et al. | Distributed generation allocation on radial distribution networks under uncertainties of load and generation using genetic algorithm | |
Syahputra et al. | Performance improvement of radial distribution network with distributed generation integration using extended particle swarm optimization algorithm | |
Amjady et al. | Short-term load forecast of microgrids by a new bilevel prediction strategy | |
CN103903073B (zh) | 一种含分布式电源及储能的微电网优化规划方法及系统 | |
CN109508823B (zh) | 一种基于场景分析法的配电网中分布式电源规划方法 | |
CN105186525B (zh) | 风电接入下无功电压控制分区方法 | |
CN107909211B (zh) | 基于模糊c均值聚类算法的风场等值建模及优化控制方法 | |
CN110298138A (zh) | 一种综合能源系统优化方法、装置、设备及可读存储介质 | |
CN106130004B (zh) | 一种综合考虑稳定特性的并网点新能源接纳能力的评估方法 | |
CN105207573A (zh) | 基于离散概率模型的风光互补供电系统定量优化配置方法 | |
CN114336785B (zh) | 一种基于网格聚类的分布式电源群控群调控制方法及装置 | |
CN110445127A (zh) | 一种面向多重随机不确定性的配电网无功优化方法及系统 | |
CN113783224A (zh) | 一种考虑多种分布式能源运行的配电网双层优化规划方法 | |
CN112651576A (zh) | 长期风电功率预测方法及装置 | |
Azizivahed et al. | Multi-area dynamic economic dispatch considering water consumption minimization, wind generation, and energy storage system | |
Zhao et al. | Probabilistic voltage stability assessment considering stochastic load growth direction and renewable energy generation | |
CN111985805A (zh) | 综合能源系统动态需求响应的方法及系统 | |
CN113609699A (zh) | 辐射状配电网交流潮流模型的计算方法和系统 | |
CN117439090A (zh) | 以灵活调节系数为指标的灵活性资源配置或调度方法 | |
Phonrattanasak et al. | Optimal location and sizing of solar farm on Japan east power system using multiobjective Bees algorithm | |
Zhen et al. | Multi-level wavelet decomposition based day-ahead solar irradiance forecasting | |
Ymeri et al. | Impact of Photovoltaic Systems Placement, Sizing on Power Quality in Distribution Network. | |
Chenxu et al. | An improved cumulant method for probabilistic load flow calculation |
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 |