CN107506854B - 一种考虑差异化场景的220kV电网网架结构规划方法 - Google Patents
一种考虑差异化场景的220kV电网网架结构规划方法 Download PDFInfo
- Publication number
- CN107506854B CN107506854B CN201710659970.6A CN201710659970A CN107506854B CN 107506854 B CN107506854 B CN 107506854B CN 201710659970 A CN201710659970 A CN 201710659970A CN 107506854 B CN107506854 B CN 107506854B
- Authority
- CN
- China
- Prior art keywords
- planning
- power
- representing
- wind
- nodes
- 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
- 238000013439 planning Methods 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000004069 differentiation Effects 0.000 title claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 49
- 230000002068 genetic effect Effects 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 41
- 238000012795 verification Methods 0.000 claims description 27
- 238000005457 optimization Methods 0.000 claims description 18
- 239000002245 particle Substances 0.000 claims description 17
- 238000004891 communication Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 241000039077 Copula Species 0.000 claims description 5
- 238000003491 array Methods 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008602 contraction Effects 0.000 claims description 3
- 238000003780 insertion Methods 0.000 claims description 2
- 230000037431 insertion Effects 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 description 14
- 230000035699 permeability Effects 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 11
- 238000013461 design Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000007142 ring opening reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种考虑差异化场景的220kV电网网架结构规划方法,涉及一种网架结构规划方法。本发明包括以下步骤:1)根据电网可靠性水平差异和新能源接入比例差异,生成多个规划场景;2)以综合投资最小为规划的目标函数,对规划的目标函数进行数学建模;3)对于不同规划方案,生成不同的导纳矩阵;4)获取最优网架序列所述220kV电网规划方法中,根据电网可靠性水平差异、新能源接入比例差异生成多个规划场景。本技术方案以综合投资最小为目标函数进行数学建模,对于不同规划方案,生成不同的导纳矩阵,并应用遗传算法进行求解。与现有技术相比,本发明综合考虑了电力系统规划中约束条件的差异性。在不同场景各项约束指标波动范围较大时,本发明有更好的适应性。
Description
技术领域
本发明涉及一种网架规划方法,尤其涉及一种考虑差异化场景的 220kV电网网架规划方法。
背景技术
首先,我国大部分地区的220kV电网处于电压等级序列的中间位置,介于输电与配电功能的中间地带,对输配电网起着重要的承上启下作用。当前,在我国负荷密度较高的大城市,220kV电网功能已由传统的输送电能的输电网转变为主要作用为分配电能的高压配电网;但在负荷密度不高的城镇和广大乡村地区,220kV电网仍然是骨干输电网。正是由于220kV 电网的承上启下作用,其结构研究就显得尤为重要。一方面,在负荷密度较高地区,220kV电网层面的短路电流问题已日益突出,成为影响电网安全运行的重要隐患,要求220kV电网结构能够尽量简单,分区运行;另一方面,在负荷密度不高地区,220kV电网作为重要供电电源,其供电可靠性问题又较为突出,要求220kV电网结构能尽量保证双电源,满足电网运行可靠性要求。220kV电网的可靠性可以用不同的指标来进行衡量。通常发输电系统可靠性评估是对发输电系统的充裕度(Adequacy)和安全性 (Security)进行评估。由此可见,对于220kV电网进行规划,所考虑的场景参数范围取值较大。从供电半径,负荷密度和电网安全等级,以及新能源渗透率来看,涉及的场景较多。对于参数取值范围较大的220kv电网规划问题,目前还缺少差异化的多场景规划建模方法。
其次,在规划220kV电压等级的电网网架结构时,出于限制短路电流和避免电磁环网的目的,广泛采用的原则是分区分层,电网解环。虽然这样的做法牺牲了电网一定的安全裕度,但在理论设计的过程中,极大地减少了所要考虑系统的节点个数,从而规避了电力系统运算当中极易出现的维数灾问题。同时,考虑分层和区域化的220kV电网结构,也可以从设计上更加明确其与500kV系统衔接时的具体分工分配。现有220kV规划尚不能做到对短路电流在方案形成阶段进行先校验,因此在寻优方向上不够准确。
再次,就目前而言,国内外的许多电力部门广泛采用确定性方法对发输电系统可靠性进行分析和评价,例如广为使用的N-1准则。N-1准则由于其使用的简单和明确性,在进行电网可靠性分析时,可以直接迅速地产生各种故障状态,并在此基础上进行分析(Contingency Analysis)。由于 N-1准则产生的故障状态是对电力系统中任意主要元件非计划退出运行的遍历,其故障状态在进行选择(Contingency Filtering)后,可以较好地代表电力系统中可能出现的重大典型故障。可以认为,N-1准则是判断电力系统可靠性的基本手段,在220kV电网的典型网架设计过程当中,无论是考虑送受端差异,还是考虑发展水平差异,N-1以至于N-2准则都是基础的,也是重要的规划原则。同时需要注意的是,在规划220kV等级网架结构时,N-1,N-2准则只是提供了故障态的便利手段,在考虑故障概率,在进一步进行差异化分析的过程当中,需要更多的量化标准来指导网架拓扑的设计与主要参数的选择。与此同时,日益增高的新能源渗透率使得N-1 校验的计算必须并入更多的机会约束,而传统的建模和求解都相对复杂耗时。现阶段在进行电网规划时,对于考虑不确定因素的机会约束处理还没有成熟的手段,先校验的N-1原则并没有得到充分的应用。
针对背景技术中提到的问题以及现有技术的不足,具体有下述问题亟待解决:
1)对于220kV电网规划参数取值波动较大,送受端差异明显的特点,缺少差异化的网架生成方法。
2)对于220kV电网的短路电流水平限制的问题上,缺少先校验的规划方法,方案求解的准确性和寻优方向的精度较差。
3)对于220kV电网的N‐1原则校验上,缺少以LOLP为指标,以先校验进行机会约束的停电事故风险校验。大规模风电并网规划中网架生成的灵活性不能得到保证。
发明内容
本发明要解决的技术问题和提出的技术任务是对现有技术方案进行完善与改进,提供一种考虑差异化场景的220kV电网网架规划方法,以达到差异化的网架生成目的。为此,本发明采取以下技术方案。
一种考虑差异化场景的220kV电网网架规划方法,其特征在于包括以下步骤:
1)根据电网可靠性水平差异和新能源接入比例差异,生成多个规划场景,包括低失负荷概率高渗透率规划场景、低失负荷概率低渗透率规划场景、高失负荷概率高渗透率规划场景、高失负荷概率低渗场景透率规划场景;
2)以综合投资最小为规划的目标函数,对规划的目标函数进行数学建模;
3)对于不同规划方案,生成不同的导纳矩阵;并应用遗传算法进行求解,排除不满足安全约束的网架结构;对满足安全约束的网架结构进行检验;其具体包括:
301)输入负荷数据,发电数据和算法参数;
302)应用Floyd法产生初始网架规划方案;
303)应用Warshall法校验系统连通性;当校验得到系统不连通时,则返回步骤302),当校验系统连通时,进入步骤304);
304)将产生的规划方案划入待选规划集合;
305)应用二进制粒子群算法对网架结构求解;
306)应用凸优化算法对LOLP约束和N-1约束进行校验;
307)若校验不通过,则返回至步骤302);若校验通过则进入下一步;
4)获取最优网架序列。
本技术方案在考虑规划参数取值范围变化较大的复杂综合规划问题时,本发明可以将不同的因素以分场景的方式进行讨论,所用规划求解计算时间随典型场景数量线性增加,在中型和大型系统中均具有良好的适应性。
在考虑网架规划模型时,对于不同的求解步骤,应用不同的算法进行计算,使得综合计算成本较低。对于网架规划模型,采用粒子群算法进行求解。对于风电出力,采用蒙特卡罗法模拟,并由边缘概率分布结合Copula 函数求得其联合概率分布。对于LOLP校验,采用凸优化算法进行快速求解。
在生成网架结构时,采用二进制粒子群算法,通过限制生成的规划方案数量来减少求解时间。同时采用Warshall算法检验系统的连通性,防止出现电力孤岛,预防节点导纳矩阵奇异导致计算不收敛的情况。
作为对上述技术方案的进一步完善和补充,本发明还包括以下附加技术特征。
进一步的,所述规划的目标函数设定为投资CI,其包括一次或二次设备投入成本IC;运行成本CO;故障成本CF。
进一步的,在步骤2)中,在确定目标函数后,模型的数学表达为:
min F(n(i,j))=CLCC
其中:N表示待选线路集合;C表示预想故障集合;L表示系统线路集合;B表示系统节点集合;T表示规划考虑连续时间集合; n(i,j)表示支路(i,j)间的规划待建线路条路;n(i,j),max表示支路(i,j)间的规划待建线路条路上限;n(i,j),max表示支路(i,j)间的原有线路条数;Pd为预测得到的负荷有功功率列向量;Pg,max为发电机有功出力上限列向量; S为节点支路关联矩阵;f表示正常情况下的支路有功功率列向量; f(i,j)表示支路(i,j)间的有功功率;f(i,j),max表示支路(i,j)间的有功功率上限;Pg表示发电机有功出力列向量;θi、θj为节点i、j的相角;在短路电流约束当中,Id表示短路电流,Id,max表示短路电流上限, Um表示故障处端口电压,Zmm表示故障处对应节点阻抗矩阵的自阻抗;在N-1约束当中,fc为故障态下的线路潮流向量,Pg,c为故障态下的发电机有功出力向量,Sc为故障态下的节点-支路关联矩阵,f(i,j),c为故障态下的支路(i,j)间的有功功率,γ(i,j)表示支路(i,j)间的电纳。
进一步的,步骤302)具体步骤包括:假设网络图中有N个顶点,通过有向图的带权矩阵中插入顶点的方法来构造矩阵,最后得到的矩阵为网络图的距离矩阵,同时记录插入点以便得到两点间的最短路径;即定义一矩阵,其元素为其值根据网络图中各顶点间的连接关系确定;插入一个顶点,两点之间距离根据下式由确定;以此类推,直到确定为止
进一步的,步骤303)中,包括:
A)通过Warshall法从图的邻接矩阵A推出可达矩阵P;得到可达矩阵 P后通过查询矩阵中对应元素的值,求得子图的连通性;子图中节点数为 N,线路数为L,则其邻接矩阵可以用两个数组IP和P2来表示;IP是一个N+1维的数组,第一个元素是0,将子图节点按节点号升序排列,之后的每一个元素代表子图的每一个节点在子图中相连的节点数,并将所相连的节点号按IP数组中节点的顺序存放在P2数组中,P2的维数则为2L;
B)对产生的备选方案进行连通性校验,计算连通性时,任意选择子图中的任意一个节点,查找其在图中的连通片,查找过程包括:a)将连通片节点自身放入连通片中,并查找其相连的节点,存放在连通片;b)查找连通片中的节点,若有其相连的节点不在连通片,将其加入连通片中,当没有新节点加入时,查找结束;c)若连通片的节点为子图的所有节点,子图连通,否则不连通。
进一步的,步骤305)中,二进制量子粒子群算法为:
更新过程简化为:
mbest表示平均最优位置;pbest表示局部最优位置;mbest表示全局最优位置;
M表示种群个数;pi表示位于局部最优pbest和全局最优gbest之间的随机位置;zid k为第i个种群第k次迭代的更新距离;φ和μ为属于[0,1]之间的随机数;
α为收缩扩张系数,用来控制算法的收敛速度;Zij k+1表示第i个种群的第j位的更新位置;prij表示其相对于随机位置的变异概率;dH(*) 为两个序列间的汉明距离。
进一步的,对于不确定性的模拟,考虑风电出力随机性带来的有功出力涨落;对于单风电场出力模拟,设风速服从韦伯分布:
其中,λ和k为对应节点风速韦伯分布的范围参数和形状参数;
设风机出力特性为:
上式中忽略风电场的尾流效应,风机可进行最大风能追踪,则可将整个风电场等效为等值风机;
其中:A为所有风机叶片迎风扫掠总面积;CP为风能转换效率,是桨距角和叶尖速比的函数;ρ为空气密度;PN为风机额定功率;vc为切入风速;vN为额定风速;vmax为切出风速。
对于多风电场接入的规划问题,采用Copula函数法,将已知的同时序下单风电场出力概率分布作为边缘分布,得到多风电场出力的联合概率密度分布。
进一步的,规划目标需满足多场景条件下的约束,具体包括:系统热备用约束、机组发电功率约束、机群出力约束、可靠性约束及线路功率约束。
有益效果:1)在考虑规划参数取值范围变化较大的复杂综合规划问题时,本发明可以将不同的因素以分场景的方式进行讨论,所用规划求解计算时间随典型场景数量线性增加,在中型和大型系统中均具有良好的适应性。
2)在考虑网架规划模型时,对于不同的求解步骤,应用不同的算法进行计算,使得综合计算成本较低。对于网架规划模型,采用粒子群算法进行求解。对于风电出力,采用蒙特卡罗法模拟,并由边缘概率分布结合 Copula函数求得其联合概率分布。对于LOLP校验,采用凸优化算法进行快速求解。
3)在生成网架结构时,采用二进制粒子群算法,通过限制生成的规划方案数量来减少求解时间。同时采用Warshall算法检验系统的连通性,防止出现电力孤岛,预防节点导纳矩阵奇异导致计算不收敛的情况。
附图说明
图1为本发明的规划方法求解流程图;
图2为低失负荷概率高渗透率场景规划系统图;
图3为低失负荷概率低渗透率场景规划系统图;
图4为高失负荷概率高渗透率场景规划系统图;
图5为高失负荷概率低渗透率场景规划系统图。
具体实施方式
以下结合说明书附图对本发明的技术方案做进一步的详细说明。
本发明包括以下步骤:
1)根据电网可靠性水平差异和新能源接入比例差异,生成多个规划场景,包括低失负荷概率高渗透率规划场景、低失负荷概率低渗透率规划场景、高失负荷概率高渗透率规划场景、高失负荷概率低渗场景透率规划场景;
2)以综合投资最小为规划的目标函数,对规划的目标函数进行数学建模;
3)对于不同规划方案,生成不同的导纳矩阵,并应用遗传算法进行求解,排除不满足安全约束的网架结构;对满足安全约束的网架结构进行检验;如图1所示,其具体包括:
301)输入负荷数据,发电数据和算法参数;
302)应用Floyd法产生初始网架规划方案;
303)应用Warshall法校验系统连通性;当校验得到系统不连通时,则返回步骤302),当校验系统连通时,进入步骤304);
304)将产生的规划方案划入待选规划集合;
305)应用二进制粒子群算法对网架结构求解;
306)应用凸优化算法对LOLP约束和N-1约束进行校验;
307)若校验不通过,则返回至步骤302);若校验通过则进入下一步输出最优网架序列;
4)获取最优网架序列。
本发明综合考虑了电力系统规划中约束条件的差异性。在不同场景各项约束指标波动范围较大时,本发明有更好的适应性。同时考虑了先校验的N-1原则和短路电流约束,规划方案可信度更高,在新能源大规模并网的背景下保证了网架生成的灵活性。
本技术方案考虑电网可靠性水平差异和新能源接入比例差异,以生成不同对应场景。对于电网可靠性水平差异,通过失负荷概率校验和N-1原则校验来保证;对于新能源渗透率差异,通过蒙特卡洛模拟方法来进行模拟。该差异化网架规划方法可以在约束条件取值较为宽泛的情况下给出典型的220kV电网网架序列,为系统规划提供参考。
规划的全周期成本目标函数为:
CLCC=m1+m2+m3+m4
m1为CI,即总投资成本,包括设备的置购成本和安装调试费用、方案整体的研究设计费用、工程前期准备费用如工程地块改造和购买费用、为满足环保要求而需额外花费的费用或者是支付的环保罚款,以及与上述投入成本有关的年度管理费用如运输费、监理费、公积金等。
m2为CO,即运行成本,包括设备运行过程中不可避免的损耗和网损。
m3为直接故障成本。直接故障成本需要从全网缺电成本的角度进行考虑。缺电成本与许多因素有关,包括停电发生的时间、停电量、停电持续时间、停电频率以及用户类型等。
m4为间接故障成本。其中间接故障成本包括赔偿费用、对社会造成的不良影响以及公司信誉受损等,这一部分成本可以通过对历史数据的统计获得一个合理的比例。设定间接故障成本与直接故障成本的比例为a, a通常是取较大量值以突出故障所带来的社会影响,那么有m4=am3。则故障总成本CF为:CF=m3+m4。
在本实施例中,规划的目标函数设定为投资CI(Comprehensive Investment),其包括一次或二次设备投入成本IC(Investment Costs);运行成本CO(Operation Costs);故障成本CF(Outage or Failure Costs)。
在确定目标函数后,模型的数学表达为:
min F(n(i,j))=CLCC
上述模型当中,
N表示待选线路集合;
C表示预想故障集合;
L表示系统线路集合;
B表示系统节点集合;
T表示规划考虑连续时间集合;
n(i,j)表示支路(i,j)间的规划待建线路条路;
n(i,j),max表示支路(i,j)间的规划待建线路条路上限;
n(i,j),max表示支路(i,j)间的原有线路条数;
Pd为预测得到的负荷有功功率列向量;
Pg,max为发电机有功出力上限列向量;
S为节点支路关联矩阵;
f表示正常情况下的支路有功功率列向量;
f(i,j)表示支路(i,j)间的有功功率;
f(i,j),max表示支路(i,j)间的有功功率上限;
Pg表示发电机有功出力列向量;
θi、θj为节点i、j的相角;
在短路电流约束当中,
Id表示短路电流;
Id,max表示短路电流上限;
Um表示故障处端口电压;
Zmm表示故障处对应节点阻抗矩阵的自阻抗;
在N-1约束当中,
fc为故障态下的线路潮流向量;
Pg,c为故障态下的发电机有功出力向量;
Sc为故障态下的节点-支路关联矩阵;
f(i,j),c为故障态下的支路(i,j)间的有功功率;
γ(i,j)表示支路(i,j)间的电纳。
规划目标需满足多场景条件下的约束。具体包括:
系统热备用约束:
∑iu(t)Pg≥Pd(t)+Ph(t)
机组发电功率约束:
Pi,min≤Pg(t)≤Pi,max
机群出力约束:
可靠性约束:
LOLP≤β
上式中失负荷概率LOLP可以表达为:
LOLP=∑t∈τ∑i∈B∑l∈L∑c∈CE(r,t,i,l,c)
其中,E(r,t,i,l,c)表示系统失负荷概率的期望值在连续时间,选定节点和选定线路以及选定故障待选集中的迭代之和。
其线路功率约束为:
0≤n(i,j)≤n(i,j),max
Sf+Pg=Pd
其中,L表示线路集合,n(i,j)为决策变量,表示待选线路走廊的架线数目; n(i,j),max为架线数目最大值;n(i,j),0表示支路(i,j)间原有线路的条数;γ(i,j)表示支路(i,j)间单条线路的导纳;f(i,j),max表示支路(i,j)间单条线路的有功传输极限。
最重要的约束条件即为线路不过载,即网络能够满足安全输送电力的要求,为此模型采取逐条线路开断后的传输容量校验。当任意一条线路断开后能够引起系统其它线路出现过负荷或系统解列时,说明网络没有满足 N-1检验。在这种情况下,必须应用某种电网 规划方法进行网络扩展,直到满足N-1检验为止。
随着场景数量增多,规划方案求解难度只随场景数线性增长。
本实施例提供一种考虑差异化场景的220kV电网网架规划方法,在进行网架结构规划时,绝大多数情况下是在已有的网架结构上进行增改。即已经满足连通性,甚至满足N-1条件下连通性的拓扑结构在规划初期就已经由实际的电网基础所决定了。在规划新的送电线路时,或者舍弃已有的冗余线路时,对电力系统网络的邻接矩阵进行连通性判断,即而已得到满足连通性校验的拓扑结构。在基于上述情况下,具体的方法步骤为:
1)在生成待选拓扑结构时,本专利采用Floyd最短路径算法。
Floyd最短路径算法的具体实施步骤为:假设网络图中有N个顶点,它通过有向图的带权矩阵中插入顶点的方法来构造矩阵,最后得到的矩阵为网络图的距离矩阵,同时记录插入点以便得到两点间的最短路径。即定义一矩阵,其元素为其值根据网络图中各顶点间的连接关系可以确定;插入一个顶点,两点之间距离根据下式由确定;以此类推,直到确定为止
连通性算法采用warshall算法。通过该算法可以从图的邻接矩阵A推出可达矩阵P。得到可达矩阵P后通过查询矩阵中对应元素的值,求得子图的连通性。
子图中节点数为N(Nunn),线路数为L(L≤m),则其邻接矩阵可以用两个数组IP和P2来表示。IP是一个N+1维的数组,第一个元素是0,将子图节点按节点号升序排列,之后的每一个元素代表子图的每一个节点在子图中相连的节点数,并将所相连的节点号按IP数组中节点的顺序存放在 P2数组中,P2的维数则为2L。
接下来对产生的备选方案进行连通性校验,计算连通性时,任意选择子图中的任意一个节点,查找其在图中的连通片,查找过程如下:(1)将连通片节点自身放入连通片中,并查找其相连的节点,存放在连通片;(2) 查找连通片中的节点,若有其相连的节点不在连通片,将其加入连通片中,当没有新节点加入时,查找结束;(3)若连通片的节点为子图的所有节点,子图连通,否则不连通。
在电力系统领域,连通性检验算法主要应用于输配电网络规划方法中,用于对输配电网规划方案的连通性进行检验,从而确定该方案是否可行。连通性检验算法可以和智能优化算法相结合,作为电网规划中连通性约束对优化问题进行求解。
差异化规划搜索的核心骨干网架,其数学模型同样需要满足基于图论的连通性约束,因此连通性检验算法可以应用于核心骨干网架的搜索,作为核心骨干网架中连通性约束的求解算法,与二进制量子粒子群优化算法和生物地理学优化算法相结合,实现核心骨干网架的搜索与构建。
2)由于本专利所涉及研究的优化问题包含约束条件较多,为了保证运算速度和最优性,优化算法采用遗传算法中的粒子群算法进行求解,用凸优化算法进行后校验。
二进制量子粒子群算法的更新表达式如下所示:
更新过程可简化为:
mbest表示平均最优位置;
pbest表示局部最优位置;
mbest表示全局最优位置;
M表示种群个数;
pi表示位于局部最优pbest和全局最优gbest之间的随机位置;
zid k为第i个种群第k次迭代的更新距离;
φ和μ为属于[0,1]之间的随机数;
α为收缩扩张系数,用来控制算法的收敛速度;
Zij k+1表示第i个种群的第j位的更新位置;
prij表示其相对于随机位置的变异概率;
dH(*)为两个序列间的汉明距离。
二进制量子粒子群算法在0-1整数规划问题上拥有良好的收敛性,已经在多个领域得到应用并取得理想效果。在电力系统领域,二进制量子粒子群算法在智能电网多目标优化规划研究、电力系统无功优化及电力系统故障诊断等方面都得到了广泛应用,体现了计算简单、复杂性低、收敛速度快和稳定性好等优点。
差异化规划搜索的核心骨干网架,有别于单纯地筛选支路或支路重要性评估,更加强调其对整个系统的作用。针对差异化规划设计的要求,要在大电网中搜索出一个符合电网特性及图论连通性的核心骨干网架是一个NP难题。因此,可以采用二进制量子粒子群算法进行核心骨干网架的搜索,建立了核心骨干网架的整体构建模型和流程。
3)对于不确定性的模拟,考虑风电出力随机性带来的有功出力涨落。对于单风电场出力模拟,设风速服从韦伯分布:
其中,λ和k为对应节点风速韦伯分布的范围参数和形状参数。
设风机出力特性为:
上式中忽略风电场的尾流效应,风机可进行最大风能追踪,则可将整个风电场等效为等值风机。
A为所有风机叶片迎风扫掠总面积;
CP为风能转换效率,是桨距角和叶尖速比的函数;
ρ为空气密度;
PN为风机额定功率;
vc为切入风速;
vN为额定风速;
vmax为切出风速。
对于多风电场接入的规划问题,采用Copula函数法,将已知的同时序下单风电场出力概率分布作为边缘分布,可得到多风电场出力的联合概率密度分布。
下面结合具体实施方案进一步详细说明本实施例,本实施例采用IEEE118节点系统进行仿真。此系统有118个节点,186条线路。其中有 54个发电机,99个节点与负荷相连。系统包含9个变压器。总装机容量为9966.2MW,平均发电量4374.9MW,平均负载4242.0MW,网损 132.86MW。对于此系统,分别针对失负荷概率为β1=5.7078e-06和β2=1.9026e-05进行电网规划。同时对新能源渗透率分别为k1=5%和k2=15%进行讨论,生成4个具体差异化算例。
Case1:低失负荷概率高渗透率
Case2:低失负荷概率低渗透率
Case3:高失负荷概率高渗透率
Case4:高失负荷概率低渗透率
将节点24、节点25、节点26和节点49替换为风电接入。待选线路集不再设定为全部已有线路,而改为15条已有线路和5条待选线路。
表1 IEEE118节点系统Case1规划结果
Case1规划结果如上表所示,规划新增线路共计12条,目标函数取值 62.5p.u.。短路电流水平为42.1kA,小于限定值50kA,如图2所示。
表2 IEEE118节点系统Case2规划结果
Case2规划结果如上表所示,规划新增线路共计8条,目标函数取值 51.7p.u.。短路电流水平为35.7kA,小于限定值50kA,如图3所示。
表3 IEEE118节点系统Case3规划结果
Case3规划结果如上表所示,规划新增线路共计6条,目标函数取值 43.3p.u.。短路电流水平为33.8kA,小于限定值50kA,如图4所示。
表4 IEEE118节点系统Case4规划结果
Case4规划结果如上表所示,规划新增线路共计3条,目标函数取值 27.4p.u.,短路电流水平31.9kA,小于限定值50kA,如图5所示。
对于实施例具体考虑的四个规划场景,本专利所提出的网架规划方法均给出了满足相应LOLP概率的规划方案。二进制粒子群算法对于网架结构的生成起到了筛选作用,对于不满足安全约束的网架结构进行了预先排除,极大地节省了计算时间。对于N-1约束校验,应用凸优化算法。对于电力系统中各个元器件的强迫停运情况,形成N-1的事故预想集。首先对于事故预想集进行筛选。对于N-1故障时,风电取极端值依旧可以满足电力系统稳态运行约束的情况,将N-1故障从中剔除。最终形成满足N-1要求的电网运行方案。
由以上四个具体实施例可以看出,风电渗透率和失负荷概率对于网架规划的保守性均有直接影响。相比较之下,Case2所对应的高失负荷概率高渗透率场景相较于Case3所对应的低失负荷概率低渗透率场景,其所需要的电网可靠性更高。其直接表现为规划所得的目标函数值次优,即最小投资成本更大。此结果同时表明,在本实施例所考虑的具体场景当中,考虑设定的风电渗透率等级k1=5%和k2=15%,失负荷概率等级β1=5.7078e-06 和β2=1.9026e-05。失负荷概率对于网架规划的方案形成影响更大。
同时需要指出的是,风电典型场景的生成精度对于规划结果同样有较大影响。在考虑系统失负荷概率LOLP水平较低,即电网比较坚强时,需要对风电生成的聚类取较高的精度,以保证在极端情况下的小概率风电出力场景不被忽略,而小概率的极端场景恰恰会对考虑可靠性的电网规划有较大影响。例如多个电场无风会使得部分线路重载,或者多个风电场出力由于风速达到切出风速,突然由额定值转为零时,对于系统的可靠性要求高于一般情况。或者在考虑小概率的不同N-1场景时,风电出力取较高或者较低值时,系统部分线路的传输容量极易越限。在生成风电出力典型场景时,如果考虑等概率场景分布,则为了达到足够的场景精度,会造成生成的场景数量过多。本方法在生成典型风电出力场景时,考虑极端情况不等概率聚类,以保证系统的失负荷概率在极端值影响下有较好的响应。
综上,本专利提出的方法对于参数波动较大的多场景220kV网架规划问题,可以生成差异化的网架序列。同时生成的规划方案满足N-1校验原则和短路电流约束,为大规模新能源并网的灵活性规划提供了新的应对策略。
以上图1所示的一种考虑差异化场景的220kV电网网架规划方法是本发明的具体实施例,已经体现出本发明实质性特点和进步,可根据实际的使用需要,在本发明的启示下,对其进行形状、结构等方面的等同修改,均在本方案的保护范围之列。
Claims (8)
1.一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于包括以下步骤:
1)根据电网可靠性水平差异和新能源接入比例差异,生成多个规划场景,包括低失负荷概率高渗透率规划场景、低失负荷概率低渗透率规划场景、高失负荷概率高渗透率规划场景、高失负荷概率低渗透率规划场景;
2)以综合投资最小为规划的目标函数,对规划的目标函数进行数学建模;
3)对于不同规划方案,生成不同的导纳矩阵,并应用遗传算法进行求解,排除不满足安全约束的网架结构;对满足安全约束的网架结构进行检验;其具体步骤包括:
301)输入负荷数据,发电数据和算法参数;
302)应用Floyd法产生初始网架规划方案;
303)应用Warshall法校验系统连通性;当校验得到系统不连通时,则返回步骤302),当校验系统连通时,进入步骤304);
304)将产生的规划方案划入待选规划集合;
305)应用二进制量子粒子群算法对网架结构求解;
306)应用凸优化算法对LOLP约束和N-1约束进行校验;
307)若校验不通过,则返回至步骤302);若校验通过则进入下一步;
4)获取最优网架序列。
2.根据权利要求1所述的一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于:所述规划的目标函数设定为投资CI,其包括一次或二次设备投入成本IC;运行成本CO;故障成本CF。
3.根据权利要求2所述的一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于:在步骤2)中,在确定目标函数后,模型的数学表达为:
min F(n(i,j))=CLCC
其中:表示待选线路集合;表示预想故障集合;表示系统线路集合;表示系统节点集合;表示规划考虑连续时间集合;n(i,j)表示支路(i,j)间的规划待建线路条路;n(i,j),max表示支路(i,j)间的规划待建线路条路上限;表示支路(i,j)间的原有线路条数;Pd为预测得到的负荷有功功率列向量;Pg,max为发电机有功出力上限列向量;S为节点支路关联矩阵;f表示正常情况下的支路有功功率列向量;f(i,j)表示支路(i,j)间的有功功率;f(i,j),max表示支路(i,j)间的有功功率上限;Pg表示发电机有功出力列向量;θi、θj为节点i、j的相角;在短路电流约束当中,Id表示短路电流,Id,max表示短路电流上限,Um表示故障处端口电压,Zmm表示故障处对应节点阻抗矩阵的自阻抗;在N-1约束当中,fc为故障态下的线路潮流向量,Pg,c为故障态下的发电机有功出力向量,Sc为故障态下的节点-支路关联矩阵,f(i,j),c为故障态下的支路(i,j)间的有功功率,γ(i,j)表示支路(i,j)间的电纳。
5.根据权利要求4所述的一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于:步骤303)中,包括:
A)通过Warshall法从图的邻接矩阵A推出可达矩阵P;得到可达矩阵P后通过查询矩阵中对应元素的值,求得子图的连通性;子图中节点数为N,线路数为L,则其邻接矩阵用两个数组IP和P2来表示;IP是一个N+1维的数组,第一个元素是0,将子图节点按节点号升序排列,之后的每一个元素代表子图的每一个节点在子图中相连的节点数,并将所相连的节点号按IP数组中节点的顺序存放在P2数组中,P2的维数则为2L;
B)对产生的备选方案进行连通性校验,计算连通性时,任意选择子图中的任意一个节点,查找其在图中的连通片,查找过程包括:a)将连通片节点自身放入连通片中,并查找其相连的节点,存放在连通片;b)查找连通片中的节点,若有其相连的节点不在连通片,将其加入连通片中,当没有新节点加入时,查找结束;c)若连通片的节点为子图的所有节点,子图连通,否则不连通。
7.根据权利要求6所述的一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于:对于不确定性的模拟,考虑风电出力随机性带来的有功出力涨落;对于单风电场出力模拟,设风速服从韦伯分布:
其中,λ和k为对应节点风速韦伯分布的范围参数和形状参数;
设风机出力特性为:
上式中忽略风电场的尾流效应,风机进行最大风能追踪,则将整个风电场等效为等值风机;
其中:A为所有风机叶片迎风扫掠总面积;CP为风能转换效率,是桨距角和叶尖速比的函数;ρ为空气密度;PN为风机额定功率;vc为切入风速;vN为额定风速;vmax为切出风速;
对于多风电场接入的规划问题,采用Copula函数法,将已知的同时序下单风电场出力概率分布作为边缘分布,得到多风电场出力的联合概率密度分布。
8.根据权利要求7所述的一种考虑差异化场景的220kV电网网架结构规划方法,其特征在于:规划目标需满足多场景条件下的约束,具体包括:系统热备用约束、机组发电功率约束、机群出力约束、可靠性约束及线路功率约束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710659970.6A CN107506854B (zh) | 2017-08-04 | 2017-08-04 | 一种考虑差异化场景的220kV电网网架结构规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710659970.6A CN107506854B (zh) | 2017-08-04 | 2017-08-04 | 一种考虑差异化场景的220kV电网网架结构规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107506854A CN107506854A (zh) | 2017-12-22 |
CN107506854B true CN107506854B (zh) | 2021-04-20 |
Family
ID=60690498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710659970.6A Active CN107506854B (zh) | 2017-08-04 | 2017-08-04 | 一种考虑差异化场景的220kV电网网架结构规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107506854B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL2031630A (en) * | 2021-04-28 | 2022-11-04 | Univ Qingdao Technology | Power grid wiring method applicable to peer-to-peer power delivery mode |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520330B (zh) * | 2018-02-13 | 2021-08-10 | 三峡大学 | 一种考虑负荷预测误差不确定性的配电网中压线路差异化规划方法 |
CN108667027B (zh) * | 2018-05-15 | 2020-04-24 | 四川大学 | 含柔性直流的交直流系统潮流转移搜索与量化分析方法 |
CN109659972B (zh) * | 2018-11-30 | 2020-10-09 | 国网浙江省电力有限公司经济技术研究院 | 基于Pair-copula理论的多点光伏出力概率预测方法和系统 |
CN109449951A (zh) * | 2018-12-07 | 2019-03-08 | 广东电网有限责任公司 | 一种电力市场背景下电力系统无功优化的方法及相关装置 |
CN109768543B (zh) * | 2018-12-18 | 2022-09-20 | 广西电网有限责任公司电力科学研究院 | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 |
CN110334391B (zh) * | 2019-05-23 | 2023-03-28 | 明阳智慧能源集团股份公司 | 一种多维度约束风电场集电线路自动规划方法 |
CN110460043B (zh) * | 2019-08-08 | 2020-11-24 | 武汉理工大学 | 基于多目标改进粒子群算法的配电网网架重构方法 |
CN110852565B (zh) * | 2019-10-10 | 2023-05-12 | 国家电网有限公司 | 一种考虑不同功能属性的输电网网架规划方法 |
CN110991727A (zh) * | 2019-11-28 | 2020-04-10 | 海南电网有限责任公司 | 一种基于潮流网损模型和线路约束模型的电网规划方法 |
CN112417631B (zh) * | 2020-12-01 | 2024-03-29 | 中国能源建设集团浙江省电力设计院有限公司 | 一种全电压等级电网防灾抗灾规划方法 |
CN113704934B (zh) * | 2021-07-28 | 2023-05-26 | 长江勘测规划设计研究有限责任公司 | 基于图论的多电缆路径规划方法 |
CN113887003B (zh) * | 2021-09-26 | 2023-05-19 | 南方电网能源发展研究院有限责任公司 | 坚强局部电网网架结构确定方法、装置、设备和存储介质 |
CN114362133B (zh) * | 2021-10-20 | 2024-05-14 | 国网内蒙古东部电力有限公司经济技术研究院 | 同质化条件下电网稳定控制方法 |
CN114511158B (zh) * | 2022-04-19 | 2022-06-21 | 南京理工大学 | 基于尾流偏转效应和2DJensen模型的风力机功率预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102545258A (zh) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | 一种大规模风电并网的电网优化规划方法 |
CN106815657A (zh) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | 一种考虑时序性和可靠性的配电网双层规划方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011088313A1 (de) * | 2011-12-12 | 2013-06-13 | Wobben Properties Gmbh | Verfahren zum Betrieb einer Windenergieanlage bzw. eines Windparks |
-
2017
- 2017-08-04 CN CN201710659970.6A patent/CN107506854B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102545258A (zh) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | 一种大规模风电并网的电网优化规划方法 |
CN106815657A (zh) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | 一种考虑时序性和可靠性的配电网双层规划方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL2031630A (en) * | 2021-04-28 | 2022-11-04 | Univ Qingdao Technology | Power grid wiring method applicable to peer-to-peer power delivery mode |
Also Published As
Publication number | Publication date |
---|---|
CN107506854A (zh) | 2017-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107506854B (zh) | 一种考虑差异化场景的220kV电网网架结构规划方法 | |
CN110969347B (zh) | 一种输电网结构形态评估方法 | |
Liang et al. | Probability-driven transmission expansion planning with high-penetration renewable power generation: A case study in northwestern China | |
CN109598377B (zh) | 一种基于故障约束的交直流混合配电网鲁棒规划方法 | |
Huang et al. | An efficient probabilistic approach based on area grey incidence decision making for optimal distributed generation planning | |
Gupta et al. | A review on optimal sizing and siting of distributed generation system: Integrating distributed generation into the grid | |
CN115640963A (zh) | 一种考虑投资运营模式的海上风电接入系统鲁棒规划方法 | |
CN103972927A (zh) | 含光伏/储能发电系统微电网的并网转孤网综合控制方法 | |
CN115017854A (zh) | 一种基于多维评估指标体系的配电网dg最大准入容量计算方法 | |
Lin et al. | Scenario generation and reduction methods for power flow examination of transmission expansion planning | |
CN113690930B (zh) | 基于nsga-iii算法的分布式光伏电源中长期选址定容方法 | |
Jaleel et al. | Reliability Evaluation of Electric Distribution Network with Distributed Generation Integrated | |
Lu et al. | Clean generation mix transition: Large-scale displacement of fossil fuel-fired units to cut emissions | |
Shurong et al. | Optimization of power collector system for large-scale offshore wind farm based on topological redundancy assessment | |
CN116415708B (zh) | 考虑置信水平的电网鲁棒规划方法 | |
Zhang et al. | Multi–objective cluster partition method for distribution network considering uncertainties of distributed generations and loads | |
CN104578057B (zh) | 一种考虑源流路径的含风电电网电压主导节点识别方法 | |
CN111293687A (zh) | 一种基于三维粒子群算法的分布式电源选址定容的方法 | |
Matute et al. | Transmission expansion planning considering the impact of distributed generation | |
CN114709825A (zh) | 基于电力电子变压器的交直流配电网动态防雷方法和系统 | |
Li et al. | Analysis of distribution network planning with distributed power supply based on ELM algorithm | |
CN109193617B (zh) | 基于紧约束识别的电力系统脆弱点评价方法及系统 | |
CN106169755A (zh) | 一种基于负荷分配均衡的配电网三双接线联络开关优化方法 | |
CN117455033B (zh) | 基于信息融合的双段式综合能源系统优化方法 | |
Bai et al. | Correlation analysis and prediction of power network loss based on mutual information and artificial neural network |
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 |