1.考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于包括以下步骤：步骤1：以多二元表判据为基础，提出暂态电压安全裕度指标，用于评估系统暂态电压稳定性；

步骤2：基于暂态电压安全裕度指标，提出动态无功补偿设备的布点方法；

步骤3：建立差异化动态无功补偿优化模型；

步骤4：通过改进的熵权优劣解距离法，筛选出各场景下的最佳方案，继而确定动态无功补偿设备在系统中的最终配置方案。

2.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤1中，提出暂态电压安全裕度指标前，先基于风电场景概率理论，构建出风电典型运行场景，如公式(1)～(6)所示：式(1)中，f(v)为风速的概率密度函数；v为风速大小；k为风速分布的形状参数；c为尺度参数，其值可由风速均值μ和标准差σ计算出；e为自然常数；

式(2)中，Pw为风电机组的输出功率；Pr为风电机组的额定容量；vci、vr、vco分别为风电机组的切入风速、额定风速、切出风速；

式(3)～(5)中，P1、P2、P3分别为风电机组运行于零输出场景、欠额输出场景和额定输出场景下的概率；

为风电机组在欠额输出场景下的输出功率大小；

零输出场景下的风电输出功率：

为风电机组在零输出场景下的输出功率大小；

额定输出场景下的风电输出功率：

为风电机组在额定输出场景下的输出功率大小。

3.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤1中，多二元表判据指的是：通过设置电压二元表(Vcr.1，Tcr.1)、(Vcr.2，Tcr.2)、…、(Vcr.n，Tcr.n)来要求某母线电压Vi低于各预设门槛值Vcr.1、Vcr.2、…、Vcr.n的最长持续时间Tb分别不超过规定时间Tcr.1、Tcr.2、…、Tcr.n，当某母线的暂态电压满足此条件时，则认为该母线的暂态电压稳定，否则暂态电压失稳。

4.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤1中，暂态电压安全裕度指标包括：

母线的暂态电压安全裕度指标：

式(7)中，σ为式(12)所示的阶跃因子；ηi为n个低电压多二元表和m个过电压二元表约束下的母线i的暂态电压安全裕度指标；ηi.d.n为n个低电压多二元表约束下的母线i的电压暂降安全裕度指标，具体如公式(8)所示；ηi.r.m为m个过电压多二元表约束下母线i的电压暂升指标，具体如公式(10)所示：式(8)中，ViN为额定电压参考值；tk和t'k分别为母线i电压在跌落过程中低于低电压二元表中的电压阈值Vcr.d.k和恢复过程中高于电压阈值Vcr.d.k的时刻；tk+1和t'k+1分别为母线i电压在跌落过程中低于低电压二元表中的电压阈值Vcr.d.k+1和恢复过程中高于电压阈值Vcr.d.k+1的时刻；Vi(t)为母线i的暂态电压响应曲线；Kn为第n个阈值区间的权重系数；

Kk为各阈值区间的权重系数，可由公式(9)逐级求解；

式(9)中，Tcr.d.1为第1个低电压二元表的时间阈值；Tcr.d.2为第2个低电压二元表的时间阈值；Tcr.d.n为第n个低电压二元表的时间阈值；Vcr.d.1为第1个低电压二元表的电压阈值；

Vcr.d.2为第2个低电压二元表的电压阈值；Vcr.d.n为第n个低电压二元表的电压阈值；K1为第1个阈值区间的权重系数；K2为第2个阈值区间的权重系数；Kn为第n个阈值区间的权重系数；

ηi.r.m为m个过电压多二元表约束下母线i的电压暂升指标，具体如公式(10)所示：式(10)中，Vcr.r.k+1为第k+1个过电压二元表的电压阈值；ti.r.k为母线i电压在扰动后的上升过程中初次高于第k个过电压二元表电压阈值的时刻；ti.r.k+1为母线i电压在扰动后的上升过程中初次高于第k+1个过电压二元表电压阈值的时刻；t'i.r.k为母线i电压在上升结束后的下降过程中初次等于第k个过电压二元表电压阈值的时刻；m为正整数；t'i.r.k+1为母线i电压在上升结束后的下降过程中初次等于第k+1个过电压二元表电压阈值的时刻；ti.r.m为母线i电压在扰动后的上升过程中初次高于第m个过电压二元表电压阈值的时刻；t'i.r.m为母线i电压在上升结束后的下降过程中初次等于第m个过电压二元表电压阈值的时刻；

Kr.m为{0≤Vi(t)≤Vcr.r.m+1}∩{ti.r.m≤t≤t'i.r.m+1}之间区域的权重系数；

Kr.k为{0≤Vi(t)≤Vcr.r.k+1}∩{ti.r.k≤t≤t'i.r.k+1}之间区域的权重系数，其数值大小可由公式(11)解得；

式(11)中，k为正整数；Tcr.r.1为第1个过电压二元表中的时间阈值；Tcr.r.2为第2个过电压二元表中的时间阈值；Tcr.r.k为第k个过电压二元表中的时间阈值；Tcr.r.k+1为第k+1个过电压二元表中的时间阈值；Kr.1为{0≤Vi(t)≤Vcr.r.2}∩{ti.r.1≤t≤t'i.r.2}之间区域的权重系数；Kr.k为{0≤Vi(t)≤Vcr.r.k+1}∩{ti.r.k≤t≤t'i.r.k+1}之间区域的权重系数；Vcr.r.1为第1个过电压二元表的电压阈值；Vcr.r.2为第2个过电压二元表的电压阈值；Vcr.r.k为第k个过电压二元表的电压阈值；Vcr.r.k+1为第k+1个过电压二元表的电压阈值；

式(12)中，e为自然常数；t为时间；ω为峰度参数，用来表征函数的陡峭程度；

区域电压合格率指标：

式(13)中，Pa为考虑系统S种运行方式、T个典型的风电场景、M个预选故障集和整个系统所有N个负荷母线后区域a的电压合格率指标；Pl为系统以运行方式l运行的概率；Na.l.v.b.d为在l运行方式及风电场景v下，负荷母线b处发生d类型故障时，区域a内电压合格母线的数量，当母线i的ηi<1时，即可认为母线i为电压合格母线；Na为区域a内的所有的负荷母线数目；PWT.v为风电以场景v运行的概率；δl.v.d为在l运行方式及风电场景v下，故障d的权重系数，数值上等于其发生的概率，由于各典型故障相互独立，故有 为母线b发生故障d的概率，这里设故障d发生在系统各母线处的概率相等，即

区域电压稳定裕度指标：

式(14)中，ηa为考虑系统S种运行方式、T个典型的风电场景、M个预选故障集和整个系统所有N个负荷母线后区域a的电压稳定裕度指标；ηi.l.v.d.b为考虑系统在l运行方式、风电场景v后母线b发生故障d时，区域a内母线i的暂态电压安全裕度指标；Pl为系统以运行方式l运行的概率；PWT.v为风电以场景v运行的概率；δl.v.d为在l运行方式及风电场景v下，故障d的权重系数； 为在l运行方式及风电场景v下，母线b发生故障d的概率。

5.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤2中，动态无功补偿设备的布点方法包括以下步骤：步骤2.1：根据公式(15)～(18)筛选出系统中的关键母线：

Bi.zs＝λiη'i.risk (15)；

式(15)中，Bi.zs为母线i的中枢值；λi为母线i的权重系数，其值由公式(16)定义；

式(16)中，Di为母线i的度，反映与母线i相关联的边的条数； 为权重系数，满足Si为母线i的实际注入功率；Sbase为系统功率的基准值； 表征母线i在系统中传输或分配功率的大小，其值越大，则代表该母线在系统中传输和分配的功率越多，在系统中的重要性程度越大；

η'i.risk为母线i的电压失稳风险因子，其值由公式(17)定义：式(17)中，Nl.v.d.i为在l运行方式及风电场景v下，母线i发生故障d时，其动态分区内电压失稳母线的总数；η'l.v.g.d为在l运行方式及风电场景v下，母线i发生故障d时，其动态分区内电压失稳母线g的暂态电压安全裕度，当η'l.v.g.d>1时，即可认为母线g暂态电压失稳；Pl为系统以运行方式l运行的概率；PWT.v为风电以场景v运行的概率；δl.v.d为在l运行方式及风电场景v下，故障d的权重系数； 为在l运行方式及风电场景v下，母线i发生故障d的概率；η'i.risk反映了母线i发生故障后，其动态分区内所有失稳母线的安全裕度指标平均大小，其值越大，则代表该母线受故障的影响程度越大；

通过计算各母线的中枢值，并根据公式(18)对其进行降序排列，从而确定出系统的关键母线集：Pbus＝{i|sort{Bi.zs},i∈{1,2,...,N}} (18)；

式(18)中，sort{Bi.zs}为各母线按照Bi.zs值大小进行降序排列组成的集合；Bi.zs为母线i的中枢值；i为母线的编号；1,2,...,N为系统内的N个负荷母线编号；

步骤2.2：由公式(19)～(23)构造出待补偿的候选母线集合；

式(19)中，SI1.i为母线i基于母线的暂态电压安全裕度的灵敏度指标；Nl.v.d.i为在l运行方式及风电场景v下，母线i发生故障d时，其动态分区内电压失稳母线的总数；ηg0.d为在母线i处安装动态无功补偿装置前，考虑n个低电压多二元表和m个过电压二元表约束后，系统发生故障d时，母线g的暂态电压安全裕度指标；ηg.d为在母线i处安装某动态无功补偿装置后，系统发生故障d时，n个低电压多二元表和m个过电压二元表约束下的母线g的暂态电压安全裕度指标；△Qc.i为在母线i处安装的动态无功补偿设备的容量；

若在暂态过程期间，对动态无功补偿设备快速电压支撑能力提出更高要求时，可用公式(20)定义的动态无功补偿设备i的系统动态无功备用来对SI1.i进一步修正；

‑t

QRTSi＝∫e (Qi‑Qi0)dt (20)；

式(20)中，QRTSi是表征接于母线i处的动态无功补偿设备的动态无功备用的量化值；Qi为接于母线i处的动态无功补偿设备在暂态过程中增发的无功功率；Qi0为稳态运行时，接于‑t母线i处的动态无功补偿设备的初始无功功率；e 为引入的衰减因子，用来量化无功补偿设备在暂态过程的各个时刻增发的无功功率，动态无功补偿设备增发无功越快对系统的暂态电压稳定越有利；

基于公式(19)、公式(20)，定义如公式(21)所示的基于母线的暂态电压安全裕度及动态无功响应速率的灵敏度指标；

式(21)中，SI2.i为基于母线的暂态电压安全裕度及动态无功响应速率的灵敏度指标；

max{QRTSi,i∈{1,2,…,N}}为各QRTSi值中的最大者；SI1.i为母线i基于母线的暂态电压安全裕度的灵敏度指标；

进一步构造出如公式(22)所示的基于SI1.i的候选母线集合和如公式(23)所示的基于SI2.i的候选母线集合；

CSI1.bus＝{i|sort{SI1.i},i∈{1,2,...,N}} (22)；

式(22)中，CSI1.bus为基于SI1.i下的候选母线集；sort{SI1.i}为各母线按照SI1.i值大小进行降序排列组成的集合；

CSI2.bus＝{i|sort{SI2.i},i∈{1,2,...,N}} (23)；

式(23)中，CSI2.bus为基于SI2.i下的候选母线集；sort{SI2.i}为各母线按照SI2.i值大小进行降序排列组成的集合。

6.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤3中，差异化动态无功补偿优化模型的目标函数如公式(24)～(26)所示:f1(x)＝ω1[(1‑Pa)+ηa] (24)；

式(24)中，Pa为区域a的电压合格率指标；ηa为区域a的电压稳定裕度指标；

min f＝{f1(x),f2(x)} (26)；

式(24)～(26)中，f1(x)和f2(x)为待优化的两个子目标函数，分别表征动态无功补偿效果和动态无功补偿经济成本；ω1、ω2为子目标函数的优化权重，满足ω1+ω2＝1；1‑Pa为系统电压失稳率指标；T2为SVC的运行年限；Csvc.u为SVC的无功补偿单价；Qsvc.u为SVC的无功补偿容量；Fsvc.u为SVC的安装费用；T1为STATCOM的运行年限；CSTATCOM.h为STATCOM的无功补偿单价；QSTATCOM.h为STATCOM的无功补偿容量；FSTATCOM.h为STATCOM的安装费用；Zl.v为在l运行方式及风电场景v下,安装SVC的补偿节点个数；Hl.v为在l运行方式及风电场景v下安装STATCOM的补偿节点个数；Ce为电价；ζ为年最大负荷小时数；Pl.v.loss为在l运行方式及风电场景v下，系统的网络损耗；minf表示两个子目标函数最小。

7.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤4中，改进的熵权TOPSIS法具体如下：

对于s个评估方案，w个评估指标，其评估步骤如下所示：

步骤(1)、运用极值法将评价矩阵A＝[axy]s×w标准化，得到标准化后的评价矩阵A'＝[a'xy]s×w，其中：a'xy取值如下：式(27)中，x为当前评估方案；y为当前评估指标；axy为指标的初始值；a'xy为指标的新值； 分别对应评估指标y下axy的最大值和最小值；

步骤(2)、根据公式(28)，计算变异系数Vy；

式(28)中， Sy分别为如公式(29)、公式(30)所示的a'xy的平均值和标准差；

步骤(3)、计算基于变异系数Vy下各指标的权重值W1y；

步骤(4)、由a'xy的离散分布情况计算其信息熵Ey；

式(32)中，Ey为信息熵；s为评估方案的总数；x为当前评估方案；a'xy为指标的新值；ln为以自然常数为底数的对数函数；

步骤(5)、计算基于信息熵Ey下各指标的权重值W2y；

步骤(6)、计算各指标的最终权重Wy，构建其加权矩阵R；

R＝(Wy×a'xy)s×w＝(rxy)s×w (35)；

式(35)中，R为加权矩阵；Wy为指标的最终权重；a'xy为指标的新值；rxy为加权矩阵中的指标值；s为评估方案的总数；w为评估指标的总数；

步骤(7)、由加权矩阵R确定各指标下的最优方案 和最劣方案

式(36)中， 为最优方案；rxy为加权矩阵中的指标值；x为当前评估方案；y为当前评估指标；w为评估指标的总数； 为指标rxy中的最大值；

式(37)中， 为最劣方案；rxy为加权矩阵中的指标值；x为当前评估方案；y为当前评估指标；w为评估指标的总数； 为指标rxy中的最小值；

步骤(8)、依次计算各待评估方案与最优方案 及最劣方案 之间的欧式距离式(38)中， 为各待评估方案与最优方案 之间的欧式距离；y为当前评估指标；w为评估指标的总数； 为最优方案；rxy为加权矩阵中的指标值；

式(39)中， 为各待评估方案与最劣方案 之间的欧式距离；y为当前评估指标；w为评估指标的总数； 为最劣方案；rxy为加权矩阵中的指标值；

步骤(9)、计算各待评估方案与理想方案的贴近度Fx；

8.根据权利要求1所述考虑暂态电压稳定的含高渗透率风电的受端电网无功规划方法，其特征在于：所述步骤4中，最终配置方案由公式(41)确定：

式(41)中，PROfin为最终的动态无功补偿设备配置方案； 为在l运行方式及风电场景v下筛选出的最佳方案的经济成本；l为系统的运行方式；S为系统运行方式的总数；v为风电的运行场景；T为典型风电场景的总数。

9.差异化动态无功补偿优化模型，其特征在于：该模型的目标函数如公式(24)～(26)所示：f1(x)＝ω1[(1‑Pa)+ηa] (24)；

式(24)中，Pa为区域a的电压合格率指标；ηa为区域a的电压稳定裕度指标；

min f＝{f1(x),f2(x)} (26)；

式(24)～(26)中，f1(x)和f2(x)为待优化的两个子目标函数，分别表征动态无功补偿效果和动态无功补偿经济成本；ω1、ω2为子目标函数的优化权重，满足ω1+ω2＝1；1‑Pa为系统电压失稳率指标；T2为SVC的运行年限；Csvc.u为SVC的无功补偿单价；Qsvc.u为SVC的无功补偿容量；Fsvc.u为SVC的安装费用；T1为STATCOM的运行年限；CSTATCOM.h为STATCOM的无功补偿单价；QSTATCOM.h为STATCOM的无功补偿容量；FSTATCOM.h为STATCOM的安装费用；Zl.v为在l运行方式及风电场景v下安装SVC的补偿节点个数；Hl.v为在l运行方式及风电场景v下安装STATCOM的补偿节点个数；Ce为电价；ζ为年最大负荷小时数；Pl.v.loss为在l运行方式及风电场景v下，系统的网络损耗；minf表示两个子目标函数最小。

10.改进的熵权TOPSIS法，其特征在于：该方法对s个评估方案，w个评估指标进行评价，包括如下步骤：步骤1)、运用极值法将评价矩阵A＝[axy]s×w标准化，得到标准化后的评价矩阵A'＝[a'xy]s×w，其中a'xy取值如下：式(27)中，x为当前评估方案；y为当前评估指标；axy为指标的初始值；a'xy为指标的新值； 分别对应评估指标y下axy的最大值和最小值；

步骤2)、根据公式(28)，计算变异系数Vy；

式(28)中， Sy分别为如公式(29)、公式(30)所示的a'xy的平均值和标准差；

步骤3)、计算基于变异系数Vy下各指标的权重值W1y；

步骤4)、由a'xy的离散分布情况计算其信息熵Ey；

式(32)中，Ey为信息熵；s为评估方案的总数；x为当前评估方案；a'xy为指标的新值；ln为以自然常数为底数的对数函数；

步骤5)、计算基于信息熵Ey下各指标的权重值W2y；

步骤6)、计算各指标的最终权重Wy，构建其加权矩阵R；

R＝(Wy×a'xy)s×w＝(rxy)s×w (35)；

式(35)中，R为加权矩阵；Wy为指标的最终权重；a'xy为指标的新值；rxy为加权矩阵中的指标值；s为评估方案的总数；w为评估指标的总数；

步骤7)、由加权矩阵R确定各指标下的最优方案 和最劣方案

式(36)中， 为最优方案；rxy为加权矩阵中的指标值；x为当前评估方案；y为当前评估指标；w为评估指标的总数； 为指标rxy中的最大值；

式(37)中， 为最劣方案；rxy为加权矩阵中的指标值；x为当前评估方案；y为当前评估指标；w为评估指标的总数； 为指标rxy中的最小值；

步骤8)、依次计算各待评估方案与最优方案 及最劣方案 之间的欧式距离式(38)中， 为各待评估方案与最优方案 之间的欧式距离；y为当前评估指标；w为评估指标的总数； 为最优方案；rxy为加权矩阵中的指标值；

式(39)中， 为各待评估方案与最劣方案 之间的欧式距离；y为当前评估指标；w为评估指标的总数； 为最劣方案；rxy为加权矩阵中的指标值；

步骤9)、计算各待评估方案与理想方案的贴近度Fx；

由公式(41)确定动态无功补偿设备在系统中的最终配置方案：

式(41)中，PROfin为最终的动态无功补偿设备配置方案； 为在l运行方式及风电场景v下筛选出的最佳方案的经济成本；l为系统的运行方式；S为系统运行方式的总数；v为风电的运行场景；T为典型风电场景的总数。