一种基于随机最优潮流的风电穿透功率极限分析方法技术

本发明专利技术公开了一种基于随机最优潮流的风电穿透功率极限分析方法。方法为：首先构建基于概率最优潮流的风电穿透功率极限计算模型，并设置系统信息以及粒子群算法初始参数；随机生成控制变量初始值和粒子初始速度；然后随机生成风电场风速样本，计算粒子在各个样本下的系统潮流，判断是否满足机会约束；接着将机会约束加入到适应度函数中，计算各粒子适应度值，获取个体最优值和全局最优值，并更新各粒子速度和位置；最后判断是否满足最大迭代次数，若不满足，则重新随机生成风电场风速样本进行计算，若满足输出最优决策变量以及风电穿透功率极限值。本发明专利技术能够满足风电随机性和潮流等式约束，从而直接计算风电接入能力。

【技术实现步骤摘要】
一种基于随机最优潮流的风电穿透功率极限分析方法
本专利技术涉及电力系统运行与控制
，特别是一种基于随机最优潮流的风电穿透功率极限分析方法。
技术介绍
作为一种可再生、无污染的绿色能源，风能的大规模开发利用在缓解能源危机、保护环境、促进可持续发展等方面的作用日益突出，受到了世界各国的高度重视，风力发电已经成为目前最成熟、最具大规模开发和商业化潜力的新能源发电技术。随着国家在政策上对可再生能源发电的大力支持，我国的风电建设也进入了一个快速发展的时期。然而，由于风能固有的随机性、间歇性和波动性，大规模风电接入后必将对电力系统的安全稳定运行造成严重的影响。因此，评估一个系统可以接纳的最大风电容量，即风电接入能力，是风电场规划阶段迫切需要解决的关键问题，也是对大规模风电接入后的系统进行分析时所要解决的首要问题。由于风电接入对系统的影响涉及因素较多，范围较广，分析计算十分复杂，因此至今尚没有统一的求解风电穿透功率极限的方法。现有的方法一般仅针对某一种特定的制约因素，无法综合考虑多个影响风电接入的制约因素，无法全面考虑系统的各种运行方式和风电场风速条件。
技术实现思路
本专利技术的目的在于提供一种能够满足风电随机性和潮流等式约束的基于随机最优潮流的风电穿透功率极限分析方法，从而直接计算风电接入能力。实现本专利技术目的的技术解决方案为：一种基于随机最优潮流的风电穿透功率极限分析方法，包括以下步骤：步骤1、构建基于概率最优潮流的风电穿透功率极限计算模型，并设置系统信息以及粒子群算法初始参数；步骤2、随机生成控制变量初始值和粒子初始速度；步骤3、随机生成风电场风速样本，计算粒子在各个样本下的系统潮流，判断是否满足机会约束；步骤4、将机会约束加入到适应度函数中，计算各粒子适应度值，获取个体最优值和全局最优值，并更新各粒子速度和位置；步骤5、判断是否满足最大迭代次数，若不满足，则返回步骤3；若满足，则进行步骤6；步骤6、输出最优决策变量以及风电穿透功率极限值。进一步地，步骤1中所述的构建基于概率最优潮流的风电穿透功率极限计算模型，具体如下：该模型以系统可接纳的各风电场装机容量之和最大化作为目标，目标函数为：式中，m为风电场个数，ni为第i个风电场中风机的个数，PNWi为第i个风电场中风机的额定功率；等式约束为系统的潮流方程，约束为：式中，PGi、QGi分别为节点i处的常规发电机组的有功和无功出力，PWi、QWi分别为节点i处的风电场的有功和无功出力，PLi、QLi分别为节点i处的有功和无功负荷，Ui、Uj、θij分别为节点i和节点j的电压幅值和相角差，Gij、Bij分别为系统导纳矩阵中的实部和虚部，CPQ、CPV分别为PQ、PV节点的集合；不等式约束包括决策变量和状态变量的约束，其中决策变量为风机个数和常规机组出力，约束为：式中，为第i个风电场中风机的最大配备个数，N*为正整数，分别为第i台常规发电机组的最小和最大有功出力，CG为常规发电机组的集合；状态变量包括节点电压幅值、常规发电机组无功出力、线路潮流以及系统的上下旋转备用，约束为：式中，分别为节点i的最小和最大电压幅值，分别为第i台常规发电机组的最小和最大无功出力，PLi为第i条线路上的潮流，为第i条线路上的潮流最大限值，CL为线路的集合，分别为系统的上、下旋转备用，取值为系统总负荷的5％，α1～α5分别为相应状态变量所对应不等式约束的置信水平。进一步地，步骤3中所述的随机生成风电场风速样本，具体如下：风速近似服从双参数的威布尔分布，概率密度函数f(v)为：式中，v为风电场的风速；k为威布尔分布的形状系数，取值范围为1.8～2.3；c为威布尔分布的尺度系数，表示某一时间段内该地区的平均风速；风电机组有功出力与风速的关系如下：式中，vin、vout、vN分别为风电机组的切入风速、切出风速和额定风速，p、pN分别为风电机组的实际输出功率和额定输出功率；采用随机模拟算法处理带有随机变量ξ的机会约束：Pr{g(x,ξ)≤0}≥α(19)随机模拟算法如下：a)置计数器N′＝0；b)根据随机变量ξ的概率分布Φ(ξ)生成随机变量ξ；c)如果g(x,ξ)≤0成立，则N′＝N′+1；d)重复步骤b)和c)共N次；e)如果N′/N≥α，则机会约束成立，否则不成立。进一步地，步骤4所述的将机会约束加入到适应度函数中，计算各粒子适应度值，获取个体最优值和全局最优值，并更新各粒子速度和位置，具体如下：计算每个初始粒子的风速样本是否满足机会约束条件，并计算适应度函数为：式中，x＝1时表示满足机会约束，x＝0则相反；设定粒子群共由N个粒子组成，每个粒子定义为D维空间，则根据下式进行更新：式中，i＝1,2....,M为粒子的个数；d＝1,2....,D为粒子的维数，即待优化问题的解的维数；c1、c2为学习因子；r1、r2为(0,1)上均匀分布的随机数；ω为惯性权重；分别为粒子i在第k次迭代的速度和位置；分别为粒子i的个体历史最优值和全部粒子的全局历史最优值；其中惯性权重ω采用非线性递减策略，以凹函数递减：ω＝(ωstart-ωend)(t/tmax)2+(ωend-ωstart)(2t/tmax)+ωstart(21)式中，ωstart、ωend分别为初始惯性权重和终止惯性权重；t、tmax分别为当前迭代次数和最大迭代次数；整数型决策变量ni速度更新的公式为：式中，int表示取整函数；表示区间上均匀分布的随机数，当时，取取当时，取取在每次速度更新后，判断速度是否越限，如果越限，则根据下式对速度进行修正：对于搜索空间限制在[Xmin,Xmax]的粒子，设定最大速度vmax为：vmax＝λ(Xmax-Xmin)/2,0.1≤λ≤1(24)本专利技术与现有技术相比，其显著优点在于：(1)针对风电的随机性，能够直接计算风电接入能力；(2)有效地计算出考虑风速随机性情况下各接入点的风电穿透功率极限值，并对风电场的规划设计进行指导。附图说明图1是本专利技术基于随机最优潮流的风电穿透功率极限分析方法的流程示意图。具体实施方式本专利技术基于随机最优潮流的风电穿透功率极限分析方法，针对风电的随机性，解决了传统最优潮流无法直接运用于风电接入能力计算的问题。本方法首先构建了基于概率最优潮流的风电穿透功率极限计算模型，针对风速服从威布尔分布，建立机会约束规划模型，并采用随机模拟技术处理机会约束，最后采用粒子群算法求解整体模型，输出最优决策变量以及风电穿透功率极限值。结合图1，本专利技术提出的基于随机最优潮流的风电穿透功率极限分析方法，具体步骤如下：步骤1、构建基于概率最优潮流的风电穿透功率极限计算模型，并设置系统信息以及粒子群算法初始参数，具体如下：基于概率最优潮流的风电穿透功率极限计算模型结合机会约束规划理论，在满足系统潮流等式约束和一系列系统可靠安全运行的不等式约束前提下，以系统可接纳的各风电场装机容量之和最大化作为目标，选取各风电场风机的配备个数和常规机组的有功出力作为决策变量进行优化调整。从本质上讲，此问题是一个含随机变量的多变量、多约束、非线性混合整数规划问题。该模型以系统可接纳的各风电场装机容量之和最大化作为目标，目标函数为：式中，m为风电场个数，ni为第i个风电场中风机的个数，PNWi为第i个风电场中风机的额定功率；等式约本文档来自技高网...

【技术保护点】
1.一种基于随机最优潮流的风电穿透功率极限分析方法，其特征在于，包括以下步骤：步骤1、构建基于概率最优潮流的风电穿透功率极限计算模型，并设置系统信息以及粒子群算法初始参数；步骤2、随机生成控制变量初始值和粒子初始速度；步骤3、随机生成风电场风速样本，计算粒子在各个样本下的系统潮流，判断是否满足机会约束；步骤4、将机会约束加入到适应度函数中，计算各粒子适应度值，获取个体最优值和全局最优值，并更新各粒子速度和位置；步骤5、判断是否满足最大迭代次数，若不满足，则返回步骤3；若满足，则进行步骤6；步骤6、输出最优决策变量以及风电穿透功率极限值。

【技术特征摘要】
1.一种基于随机最优潮流的风电穿透功率极限分析方法，其特征在于，包括以下步骤：步骤1、构建基于概率最优潮流的风电穿透功率极限计算模型，并设置系统信息以及粒子群算法初始参数；步骤2、随机生成控制变量初始值和粒子初始速度；步骤3、随机生成风电场风速样本，计算粒子在各个样本下的系统潮流，判断是否满足机会约束；步骤4、将机会约束加入到适应度函数中，计算各粒子适应度值，获取个体最优值和全局最优值，并更新各粒子速度和位置；步骤5、判断是否满足最大迭代次数，若不满足，则返回步骤3；若满足，则进行步骤6；步骤6、输出最优决策变量以及风电穿透功率极限值。2.根据权利要求1所述的基于随机最优潮流的风电穿透功率极限分析方法，其特征在于，步骤1中所述的构建基于概率最优潮流的风电穿透功率极限计算模型，具体如下：该模型以系统可接纳的各风电场装机容量之和最大化作为目标，目标函数为：式中，m为风电场个数，ni为第i个风电场中风机的个数，PNWi为第i个风电场中风机的额定功率；等式约束为系统的潮流方程，约束为：式中，PGi、QGi分别为节点i处的常规发电机组的有功和无功出力，PWi、QWi分别为节点i处的风电场的有功和无功出力，PLi、QLi分别为节点i处的有功和无功负荷，Ui、Uj、θij分别为节点i和节点j的电压幅值和相角差，Gij、Bij分别为系统导纳矩阵中的实部和虚部，CPQ、CPV分别为PQ、PV节点的集合；不等式约束包括决策变量和状态变量的约束，其中决策变量为风机个数和常规机组出力，约束为：式中，为第i个风电场中风机的最大配备个数，N*为正整数，分别为第i台常规发电机组的最小和最大有功出力，CG为常规发电机组的集合；状态变量包括节点电压幅值、常规发电机组无功出力、线路潮流以及系统的上下旋转备用，约束为：式中，分别为节点i的最小和最大电压幅值，分别为第i台常规发电机组的最小和最大无功出力，PLi为第i条线路上的潮流，为第i条线路上的潮流最大限值，CL为线路的集合，分别为系统的上、下旋转备用，取值为系统总负荷的5％，α1～α5分别为相应状态变量所对应不等式约束的置信水平。3.根据权利要求1所述的基于随机最优潮流的风电穿透功率极限分析方法，其特征在于，步骤3中所述的随机生成风电场风速样本，具体如下：...

【专利技术属性】
技术研发人员：涂杉杉，卜京，郑文哲，焦永辉，张宁宇，周前，殷明慧，谢云云，姚娟，
申请(专利权)人：南京理工大学，国网江苏省电力有限公司电力科学研究院，
类型：发明
国别省市：江苏,32