CN108306303A - 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 - Google Patents
一种考虑负荷增长和新能源出力随机的电压稳定评估方法 Download PDFInfo
- Publication number
- CN108306303A CN108306303A CN201810045708.7A CN201810045708A CN108306303A CN 108306303 A CN108306303 A CN 108306303A CN 201810045708 A CN201810045708 A CN 201810045708A CN 108306303 A CN108306303 A CN 108306303A
- Authority
- CN
- China
- Prior art keywords
- load
- power
- growth
- node
- new energy
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 9
- 230000012010 growth Effects 0.000 claims abstract description 83
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000009826 distribution Methods 0.000 claims abstract description 37
- 230000006641 stabilisation Effects 0.000 claims abstract description 17
- 238000013097 stability assessment Methods 0.000 claims abstract description 17
- 238000011105 stabilization Methods 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000010248 power generation Methods 0.000 claims description 29
- 230000005611 electricity Effects 0.000 claims description 17
- 230000007773 growth pattern Effects 0.000 claims description 17
- 239000003102 growth factor Substances 0.000 claims description 16
- 238000005516 engineering process Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 240000002853 Nelumbo nucifera Species 0.000 claims description 4
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 4
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 4
- 230000005283 ground state Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 2
- 244000131316 Panax pseudoginseng Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 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
- 230000005684 electric field Effects 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 238000005192 partition Methods 0.000 claims description 2
- 238000003064 k means clustering Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 9
- 241000196324 Embryophyta Species 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/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)建立负荷增长概率分布模型
对于聚类后的各负荷节点群，将各节点群有功负荷增长系数γj列为随机变量，假定负荷增长波动满足正态分布，即，正态分布的期望μj＝ηj,t，标准差为σj，则γj的概率密度函数为
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为区域基态有功负荷，本发明考虑的负荷裕度统计特征如下：
1)区域相对负荷裕度：
2)期望：
3)标准差：
4)最大值LMmin与最小值LMmax及其各自对应的负荷增长方向Dmin，Dmax。
5)概率密度分布信息。
连续潮流通过在常规潮流方程中引入负荷变化参数并增加一维扩展方程，采用预测-校正的方法解决常规潮流方程在临界点附近雅克比矩阵奇异的问题。连续潮流主要过程包括预测、校正、参数化策略和步长控制四部分，其计算示意图如图3所示。
下面结合图4至图8，以IEEE 118节点***为具体实例，说明本发明所属模型和方法的可行性和有效性。
IEEE118节点***接线图如图4所示，将整个***分为两个区域，分别为区域A和区域B。风电场和光伏发电***数据如表1和表2所示。在IEEE118节点***中，节点23、39和114分别接入风电场，节点44和118分别接入光伏发电***。
表1风电场相关参数
| 风电场 | Nw | Pr/MW | vci/(m/s) | vr/(m/s) | kw | cw |
| 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 | η | rmax | α | β |
| 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 (6)
1.一种考虑负荷增长和新能源出力随机的电压稳定评估方法，其特征在于，包括以下具体步骤：
1)建立负荷增长和新能源出力随机的电压稳定负荷裕度计算模型；
2)根据评估意愿和先验知识指定负荷和发电增长范围，获取增长范围内各节点相关参数；
3)对指定负荷增长范围内负荷节点的历史负荷数据采用改进K均值聚类方法进行分类，得到各负荷节点群，根据负荷分类和负荷预测计算负荷增长范围内节点群负荷增长方向，基于负荷增长方向和新能源出力的随机性分别建立概率分布模型；
4)根据各概率分布通过拉丁超立方采样技术得到负荷增长方向和新能源出力采样样本；
5)根据各样本的负荷增长方向和新能源出力样本，利用连续潮流计算有功负荷裕度。直至计算完所有采样样本，得到各样本下的负荷裕度值，求取负荷裕度的统计特征，分析电压稳定概率评估结果。
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所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法，其特征在于，所述步骤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)建立负荷增长概率分布模型
对于聚类后的各负荷节点群，将各节点群有功负荷增长系数γj列为随机变量，假定负荷增长波动满足正态分布，即正态分布的期望μj＝ηj,t，标准差为σj，则γj的概率密度函数为
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分别列为随机变量，则的新能源机组随机出力可表示如下：
5.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法，其特征在于，所述步骤4)的拉丁超立方采样方法包括如下过程：
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 -1(·)是Fk(·)的反函数。
42)经过步骤41)采样过程得到初始采样矩阵XN×M相关性较高，用cholesky分解法进行排序，通过改变采样值的排列顺序来减小相关性。
43)经过41)和42)排序两个步骤得到满足各随机变量概率分布且相关性较低的负荷方向和新能源出力样本集。
6.根据权利要求1所述的一种考虑负荷增长和新能源出力随机的电压稳定评估方法，其特征在于，所述步骤5)包括如下过程：
51)针对每一个采样样本中的新能源出力PRi和QRi和负荷增长系数kPi和kQi，采用连续潮流求解电压稳定负荷裕度值。
52)重复步骤51)，直至计算完所有采样样本，得到不同样本下各个负荷裕度值。
53)根据各负荷裕度值求取负荷裕度统计特征，设LM(Xi)为第i个样本Xi对应的负荷裕度，M为样本总数，P0为区域基态有功负荷，考虑的负荷裕度统计特征如下：
1)区域相对负荷裕度：
2)期望：
3)标准差：
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 true CN108306303A (zh) | 2018-07-20 |
| CN108306303B 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) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108920868A (zh) * | 2018-07-23 | 2018-11-30 | 国网江西省电力有限公司电力科学研究院 | 基于智能电表数据的负荷辨识模型参数修正方法 |
| CN109119987A (zh) * | 2018-09-05 | 2019-01-01 | 国网天津市电力公司 | 基于电网灵活性的城市电网不平衡功率波动范围评估方法 |
| CN109450000A (zh) * | 2017-11-09 | 2019-03-08 | 广东电网有限责任公司电力调度控制中心 | 一种基于负荷率调整方向的发电计划偏差电量分配方法 |
| CN109617083A (zh) * | 2018-12-26 | 2019-04-12 | 贵州电网有限责任公司 | 基于粒子群算法的配电网故障修复方法 |
| CN110096730A (zh) * | 2019-03-15 | 2019-08-06 | 国网辽宁省电力有限公司 | 一种电网电压快速评估方法及*** |
| CN110460085A (zh) * | 2019-08-28 | 2019-11-15 | 华北电力大学（保定） | 一种考虑风电和负荷特性对电力***影响的方法 |
| CN110531271A (zh) * | 2019-10-14 | 2019-12-03 | 深圳市蓝度汽车电控技术有限公司 | 一种汽车电瓶静态电压的检测方法、***及终端 |
| CN110932277A (zh) * | 2019-12-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | 一种电网静态电压稳定负荷裕度分析方法、装置和设备 |
| CN111463796A (zh) * | 2020-02-28 | 2020-07-28 | 广东电网有限责任公司 | 一种考虑风电出力的电压稳定性评估方法及装置 |
| CN112366713A (zh) * | 2020-11-30 | 2021-02-12 | 南方电网科学研究院有限责任公司 | 交直流混联电网静态电压稳定计算方法、装置及存储介质 |
| CN112564181A (zh) * | 2020-11-25 | 2021-03-26 | 云南电网有限责任公司 | 考虑新能源出力不确定性的在线概率安全监测方法 |
| CN113328467A (zh) * | 2021-05-26 | 2021-08-31 | 广东电网有限责任公司 | 一种概率电压稳定性评估方法、***、终端设备及介质 |
| CN113452028A (zh) * | 2021-07-29 | 2021-09-28 | 广东电网有限责任公司 | 低压配电网概率潮流计算方法、***、终端和存储介质 |
| 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 |
|---|
| H. YU等: "Probabilistic load flow evaluation with hybrid Latin hypercube sampling and Cholesky decomposition", 《IEEE》 * |
| 伍利等: "实际负荷变化模式下的静态电压稳定分析", 《电力***及其自动化学报》 * |
| 李俊芳等: "基于进化算法改进拉丁超立方抽样的概率潮流计算", 《中国电机工程学报》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109450000B (zh) * | 2017-11-09 | 2021-07-30 | 广东电网有限责任公司电力调度控制中心 | 一种基于负荷率调整方向的发电计划偏差电量分配方法 |
| CN109450000A (zh) * | 2017-11-09 | 2019-03-08 | 广东电网有限责任公司电力调度控制中心 | 一种基于负荷率调整方向的发电计划偏差电量分配方法 |
| CN108920868A (zh) * | 2018-07-23 | 2018-11-30 | 国网江西省电力有限公司电力科学研究院 | 基于智能电表数据的负荷辨识模型参数修正方法 |
| CN109119987A (zh) * | 2018-09-05 | 2019-01-01 | 国网天津市电力公司 | 基于电网灵活性的城市电网不平衡功率波动范围评估方法 |
| CN109617083A (zh) * | 2018-12-26 | 2019-04-12 | 贵州电网有限责任公司 | 基于粒子群算法的配电网故障修复方法 |
| CN110096730A (zh) * | 2019-03-15 | 2019-08-06 | 国网辽宁省电力有限公司 | 一种电网电压快速评估方法及*** |
| CN110096730B (zh) * | 2019-03-15 | 2023-04-07 | 国网辽宁省电力有限公司 | 一种电网电压快速评估方法及*** |
| CN110460085A (zh) * | 2019-08-28 | 2019-11-15 | 华北电力大学（保定） | 一种考虑风电和负荷特性对电力***影响的方法 |
| CN110531271A (zh) * | 2019-10-14 | 2019-12-03 | 深圳市蓝度汽车电控技术有限公司 | 一种汽车电瓶静态电压的检测方法、***及终端 |
| CN110932277A (zh) * | 2019-12-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | 一种电网静态电压稳定负荷裕度分析方法、装置和设备 |
| CN111463796A (zh) * | 2020-02-28 | 2020-07-28 | 广东电网有限责任公司 | 一种考虑风电出力的电压稳定性评估方法及装置 |
| CN111463796B (zh) * | 2020-02-28 | 2021-10-08 | 广东电网有限责任公司 | 一种考虑风电出力的电压稳定性评估方法及装置 |
| CN112564181A (zh) * | 2020-11-25 | 2021-03-26 | 云南电网有限责任公司 | 考虑新能源出力不确定性的在线概率安全监测方法 |
| CN112564181B (zh) * | 2020-11-25 | 2023-05-26 | 云南电网有限责任公司 | 考虑新能源出力不确定性的在线概率安全监测方法 |
| CN112366713B (zh) * | 2020-11-30 | 2023-02-28 | 南方电网科学研究院有限责任公司 | 交直流混联电网静态电压稳定计算方法、装置及存储介质 |
| CN112366713A (zh) * | 2020-11-30 | 2021-02-12 | 南方电网科学研究院有限责任公司 | 交直流混联电网静态电压稳定计算方法、装置及存储介质 |
| CN113328467A (zh) * | 2021-05-26 | 2021-08-31 | 广东电网有限责任公司 | 一种概率电压稳定性评估方法、***、终端设备及介质 |
| CN113452028A (zh) * | 2021-07-29 | 2021-09-28 | 广东电网有限责任公司 | 低压配电网概率潮流计算方法、***、终端和存储介质 |
| CN113887809A (zh) * | 2021-10-11 | 2022-01-04 | 国网新疆电力有限公司巴州供电公司 | 双碳目标下配电网供需平衡方法、***、介质及计算设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108306303B (zh) | 2021-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108306303A (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
| CN107301472B (zh) | 基于场景分析法和电压调节策略的分布式光伏规划方法 | |
| Ganguly et al. | Distributed generation allocation on radial distribution networks under uncertainties of load and generation using genetic algorithm | |
| Liang et al. | Probability-driven transmission expansion planning with high-penetration renewable power generation: A case study in northwestern China | |
| CN107591844A (zh) | 考虑节点注入功率不确定性的主动配电网鲁棒重构方法 | |
| CN107910863A (zh) | 综合考虑光伏出力和负荷需求预测区间的配电网调度方法 | |
| Capizzi et al. | Recurrent neural network-based control strategy for battery energy storage in generation systems with intermittent renewable energy sources | |
| CN106130004B (zh) | 一种综合考虑稳定特性的并网点新能源接纳能力的评估方法 | |
| CN107069814A (zh) | 配网分布式电源容量布点的模糊机会约束规划方法与*** | |
| CN109560574A (zh) | 一种考虑不确定因素的智能配电网网架规划方法 | |
| CN105186499A (zh) | 一种配电网多目标概率最优潮流模糊建模与求解方法 | |
| CN110380444A (zh) | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 | |
| CN112418496B (zh) | 一种基于深度学习的配电台区储能配置方法 | |
| CN116822719A (zh) | 一种配电网的多目标规划方法及装置 | |
| Han et al. | Optimal sizing considering power uncertainty and power supply reliability based on LSTM and MOPSO for SWPBMs | |
| Dasan et al. | Optimal siting and sizing of hybrid distributed generation using EP | |
| CN108694475B (zh) | 基于混合模型的短时间尺度光伏电池发电量预测方法 | |
| CN113344283B (zh) | 基于边缘智能的能源互联网新能源消纳能力评估方法 | |
| CN109726416B (zh) | 一种基于新能源集群预测及潮流计算的调度决策方法 | |
| Ni et al. | A review of line loss analysis of the low-voltage distribution system | |
| CN113609699A (zh) | 辐射状配电网交流潮流模型的计算方法和*** | |
| CN106655266B (zh) | 一种接入新能源的地区电网可灵活调节电源配置方法 | |
| CN109586309B (zh) | 基于大数据自由熵理论和场景匹配的配电网无功优化方法 | |
| CN112039051A (zh) | 一种含双馈风力发电机接入变电站母线负荷的实时建模方法 | |
| Fu et al. | A spatial forecasting method for photovoltaic power generation combined of improved similar historical days and dynamic weights allocation |
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 |