基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统可靠性评估方法
【专利摘要】本发明公开了一种基于几何优化?最小方差化方法的重要抽样蒙特卡洛电力系统可靠性评估方法，首次将规划领域的几何优化方法和重要抽样中的方差最小化模型结合在一起，利用几何优化求解可靠性评估中的方差最小化模型，然后利用求解得到的重要抽样参数进行可靠性评估，本发明相比于传统的基于重要抽样的可靠性评估方法，具有收敛方差更小、收敛速度更快和评估精度更高的特点，通过提高系统可靠性评估的评估速度和评估精度，尤其是捕捉那些小概率/高影响的稀有事件，使之能够快速评估电网规划不同方案之间的可靠性的细微差别，为电网规划方案的选择提供定量准确的辅助参考依据。
【专利说明】
基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系 统可靠性评估方法
技术领域
[0001] 本发明属于电力系统规划评估领域，设及能够快速提高可靠性评估的速度且考虑 稀有事件的基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统可靠性评估方 法。
【背景技术】
[0002] 近年来，由于电力系统的网架日益坚强，人们对电力系统中的极端事件越发关注， 而运些极端的小概率/高影响事件一般很难抽样，面对运些事件，传统的可靠性评估方法常 常显得过于无力。另一方面，计算运些小概率的事件对于电力系统的规划十分有意义，由于 现在的电网网架一般都非常坚强，因此能够准确计算出运些方案在可靠性指标上的差异就 显得尤其重要。运样的方法可W对电力系统的扩张规划提供合理的参考建议。目前，重要抽 样的方法十分流行于电力系统的可靠性评估中，虽然它能够加快整个评估过程的收敛速 度，但是评估的精度问题却很少得到讨论。
[0003] 比如，目前十分流行的交叉赌方法虽然能够显著加快可靠性评估的速度，但是其 精度究竟怎样还很少得到讨论。与此同时，最小方差技术一直被认为具有更好的收敛特性 和估计误差，因此被建议使用在电力系统的可靠性评估中。然而，方差最小化模型一般比较 难W求解，因此合理的求解该问题是一个难点。
【发明内容】
[0004] 针对上面提到的问题，本发明的目的在于提出一种基于几何优化-最小方差化方 法的重要抽样蒙特卡洛电力系统可靠性评估方法。
[0005] 为达到上述目的，本发明采用了 W下技术方案：
[0006] 1)预抽样阶段
[0007] 1.1)确定预抽样阶段需要的参数，包括系统的峰值负荷水平、线路容量、线路原始 强迫停运率UL、发电机组发电额定容量W及发电机组原始强迫停运率UG，定义预抽样的样本 数No;
[000引1.2)置模拟迭代次数k = 0，置所有状态的状态变量为0，即：
[0009] H(IiNo)=O
[0010] 1.3)令k = k+l;然后抽样系统状态Xk,令系统中的元件为两状态模型，记录第k次 抽样元件的状态，元件的状态由发电机组状态Xg和线路状态Xl组成，对元件的状态进行评估 分析后计算系统是否会失负荷，然后统计评估分析的结果，如果失负荷，贝化(Xk) = I;否则， H(Xk) =0;H( ?)表示数组；
[0011] 1.4)判断k是否达到No，如果没有达到，则转至步骤1.3)，如果达到，则利用几何优 化求解基于重要抽样的电力系统可靠性评估的方差最小化模型，从而得到新的发电机组强 迫停运率VG和线路强迫停运率VL ;
[0012] 2)主抽样阶段
[0013] 2.1)定义主抽样的样本数化，对可靠性指标进行初始化，并置模拟迭代次数n = 0;
[0014] 2.2)令n = n+l;然后利用新的元件强迫停运率VG和VL抽样系统的状态，对抽样的状 态进行状态评估分析，并根据是否发生失负荷更新可靠性指标；
[001引2.3)判断n是否达到化，如果没有达到，则转至步骤2.2)，否则，结束可靠性评估。
[0016] 所述利用几何优化求解基于重要抽样的电力系统可靠性评估的方差最小化模型 具体包括W下步骤:将建立的基于重要抽样的电力系统可靠性评估的方差最小化模型松弛 为满足几何优化模型的几何优化-方差最小化模型，然后求解该模型，得到新的重要抽样参 数，即新的元件强迫停运率。
[0017] 所述基于重要抽样的电力系统可靠性评估的方差最小化模型的建立具体包括W 下步骤:按照原始强迫停运率预抽样No次后，根据得到的评估结果按W下公式形成基于重 要抽样的电力系统可靠性评估的方差最小化模型：
[001 引
[0019]其中，Ng表示系统中发电机组的数目，Nl表示系统中线路的数目。[0020] 所述几何优化-方差最小化模型的获取方法包括W下步骤：
[0021] 引入 虚拟变j
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
[002引
[0029]
[0030] 通过对数变换和内点法求解所述几何优化-方差最小化模型。
[0031] 所述可靠性指标包括电力不足概率LOLP和电力电量不足期望值邸NS,可靠性指标 按照W下公式进行更新：
[0032] LOLP = LOLP+! XW(Xn；u,V)
[0033] 邸 NS =邸 NS+APnXW(Xn;u，v)
[0034] 其中，A Pn表示第n次抽样导致的负荷损失;W(Xn;u，v)表示第n次抽样得到的系统 状态的似然率比，该似然率比表示如下：
[0035]
[0036] 其中，Ng表示系统中发电机组的数目，化表示系统中线路的数目。
[0037] 本发明的有益效果体现在：
[0038] 本发明利用几何优化方法高效求解电力系统可靠性评估中的方差最小模型得到 新的重要抽样参数，从而使用该参数完成可靠性评估，得到系统的可靠性评估指标。相比于 流行的交叉赌等重要抽样方法，因为求解的新的重要抽样参数更加接近理想的重要抽样参 数，因此本发明能够在评估收敛速度和评估精度上取得优势。相比于传统的蒙特卡洛方法， 本发明能够极大的提高评估效率，通过提高系统可靠性评估的评估速度和评估精度，尤其 是捕捉那些小概率/高影响的稀有事件(该方法对于稀有事件的捕捉更为有益，主要是重要 抽样方法的性质决定了运一点），使之能够快速评估电网规划不同方案之间的可靠性的细 微差别，为电网规划方案的选择提供定量准确的辅助参考依据，尤其是输电网的扩展规划。 本发明可W有效提高可靠性评估的计算效率，是一种极有潜力应用在电网规划中的可靠性 评估方法。
【附图说明】
[0039] 图1为基于几何优化-最小方差化方法的重要抽样蒙特卡洛可靠性评估方法的流 程框图。
[0040] 图2为状态评估流程框图。
【具体实施方式】
[0041 ]下面结合附图和实施例对本发明做详细说明。
[0042] 本发明是一种基于几何优化求解方差最小模型的重要抽样可靠性评估方法，该方 法利用几何优化(GP)高效求解可靠性评估中的方差最小化(VM)模型，得到需要的重要抽样 参数，然后利用运些参数进行相关的可靠性评估。尤其是针对那些高度可靠的系统展开可 靠性评估。其具体分为预抽样中利用几何优化求解方差最小模型，从而求解可靠性评估参 数，和主抽样中进行可靠性评估两阶段：
[0043] Level 1预抽样(求解参数）：首先预抽样产生初始样本，然后在运些样本的基础 上，利用几何优化求解可靠性评估VM模型，从而得到需要的重要抽样可靠性评估参数。
[0044] Level 2主抽样(评估）：在得到的重要抽样可靠性评估参数的基础上，对系统的 可靠性展开评估，得到最终的可靠性指标。
[0045] 本发明基于几何优化-最小方差(GP-VM)的可靠性评估的流程如图1所示，具体包 括两个阶段(预抽样阶段和主抽样阶段）：
[0046] Level 1预抽样阶段
[0047] 1)输入预抽样阶段需要的参数，包含系统的峰值负荷水平、线路(包含线路容量、 线路原始强迫停运率UL)、发电机组(包含发电额定容量、发电机组原始强迫停运率UG)信息 等原始数据，定义预抽样的样本数No;下标G表示发电机，下标L表示线路；
[004引2)置模拟迭代次数k初值为0，即k = 0，置所有状态的状态变量为0，即H( 1: No)=化
[0049] 3)k = k+l;
[0050] 4)抽样系统状态Xk，一般认为系统中的元件为两状态模型，记录第k次抽样元件的 状态(发电机状态Xek和线路状态XlO,对该状态开展评估分析，计算系统是否会失负荷，最 后统计此次可靠性评估的结果，如果失负荷H(Xk) = 1;否则，H(Xk) =0。步骤4)具体包括状态 抽样、状态评估和结果统计：
[0051] a)状态抽样:考虑所有的元件均为两状态模型，即对于第i个元件有：
[0052]
口-^
[0053] 0是一个[0，1]空间均匀分布的随机数，Ui是第i个元件的强迫停运率。
[0054] Xi是第i个元件的状态，1表示正常工作，0表示故障。一个综合系统中如果有输电 线路化条，W及发电机Ng台，其状态分别用Xli和Xg康示，其中，对发电机的状态进行抽样：
[0化5]
(2-2)
[0056] UCi是第i个发电机的强迫停运率。
[0057] 对输电线路的状态进行抽样：
[0058]
C--3)
[0059] ULi是第i条线路的强迫停运率。
[0060] b)状态评估:状态评估主要是对抽样的系统状态进行可靠性评估，主要是分析研 究该状态是否导致切负荷。
[0061] 如图2所示，状态评估的主要流程如下：
[0062] BI,判断线路是否停运(故障），若有线路停运，则进入解列分析模块;若无线路停 运，则输出结果，结束评估。
[0063] B2,分析系统是否解列，若系统解列，则在各个解列块内进行功率平衡分析，统计 各解列块内的失负荷量A P，累加后得到第k次抽样系统总的失负荷量A Pk,结束评估;如果 没有解列，则进入潮流分析模块B3。
[0064] B3,判断系统是否有孤立节点，若有孤立节点，去除孤立节点再形成导纳矩阵，然 后计算线路潮流;若无过负荷结束评估，若有过负荷进入状态校正模块评估是否有失负荷， 统计第k次系统的失负荷量APk,结束评估。其中的状态校正模型如下：
[0065]
(2-4)
[0066] 式中i是节点的编号，D是节点的负荷功率向量，P是节点的发电机注入功率向量，N 是节点总数，是各节点的负荷需求向量，pmax是各发电机的最大发电量向量，PTDF是功率 分布因子矩阵，歹是线路潮流容量上限向量，E是线路潮流容量下限向量。Pi和di分别是各 节点的注入功率和负荷功率，是决策变量。该状态校正模型的目标函数是系统总的失负荷 量最小。
[0067] C)结果统计:运里主要计算用于计算新的抽样参数所需要的状态变量：
[006引如果第k次抽样没有导致切负荷，即APk是0,那么记录此时的状态：
[0069] H(Xk)=O (2-5)
[0070] 如果第k次抽样导致切负荷，那么APk大于0,更新系统的状态：
[0071] H(Xk) = I (2-6)
[0072] 并记录此时系统中各元件的状态向量拍谢X化具体设置情况如下：
[0073]
(2-巧
[0074] 5)判断k是否到达No,如果没有到达，转向步骤3)，否则转向下面求解计算新的重 要抽样参数，即利用几何优化(GP)求解电力系统可靠性评估中的方差最小化(VM)模型，从 而得到新的元件强迫停运率VG和VL。
[0075] 对于基于重要抽样的电力系统可靠性评估，其一般可W表示为：
[0076]
[0077] 其中f (Xk)是系统的原始抽样密度函数，g(Xk)是系统的重要抽样密度函数，N代表 抽样次数；
[0078] 那么，方差最小化橫巧就是将上面的估计进行方差最小化处理，即：
[0079]
[0080] ，上面的式子可W进一步表示为：
[0081]
(2-8)
[0082] E表示求期望;对于求解新的重要抽样参数v(v表示新的强迫停运率），就是求解式 2-8的模型。
[0083] 对该模型积分形式进行离散化后得到：
[0084]
(2-9)
[0085] 式中的f( ?)是系统的概率密度函数。
[0086] 根据预抽样得到的评估结果和电力系统元件服从两状态的伯努利分布的假设，并 展开其（指公式2-9)中的f(Xk;u)和f(Xk;v)可W得到下面的电力系统可靠性评估方差最小 化模型：
[0087]
(2-10)
[0088] 也就是按照原始强迫停运率U预抽样No次后，根据得到的评估结果可按上面的式 子形成电力系统可靠性评估的方差最小化模型。
[0089] 将上述电力系统可靠性评估的方差最小化模型松弛为满足几何优化模型的几何 优化-方差最小化模型，并给出了证明。其主要思路如下：
[0090] 、虚拟变 量tj和Sj
[0091] (2-11)
[0092] 然而几何优化模型如下：
[0093] mjn !]\x)
[0094] s.t.fi(x)《l ,i = l ,2, . . . ,m
[0095] gi(x) = l,i = l,2,...,p
[0096] 其中的fi是多项式，gi是单项式，x=(xi，…，xn)是n个非负变量。其中的单项式具 有如下特点：
C三0并且ai G R。式中的多项式可W表示为许多单项式 的和
[0097] 易得出，模型（式2-11)并不满足几何优化模型的形式，因为其多项式约束是等式 约束而非几何优化要求的不等式约束，所W模型(式2-11)并不是几何优化模型，而是一个 N-P hard问题。故对该模型（式2-11)进行等价松弛，得到下面的满足几何优化模型条件的 几何优化-最小方差(GP-VM)模型：
[009引
(2-12)
[0099] 实际上，松弛得到的几何优化-最小方差(GP-VM)模型和原方差最小化(VM)模型是 等价的，下面给出证明。
[0100] 证明：首先式2-11中的目标函数是求其最小值，而它是t神日Sj的一个减函数，因此 它们应该被尽量取大。同时，t詞日Sj的上界被限定在了 1-1%和1-。因此，运个优化解必须 被取在护1-?和矿1-?。运就意味着，上面的不等式约束必须达到边界才能成立。因此 式2-11中的模型可W被松弛为式2-12的形式。证毕！
[0101] 利用几何优化对上面的模型（式2-12)进行求解得到最终需要的新的重要抽样参 数，即元件新的强迫停运率VG和VL。几何优化模型通过一个对数变换和内点法可W被有效的 求解。几何优化模型可W被视为一个黑盒，任何相关的模型可W放入其中求解。因此，通过 对数变换和内点法就可W高效求解上面的模型(式2-12)。
[0102] Level 2主抽样阶段(使用通过求解几何优化模型得到的新的重要抽样参数V来 进行评估）
[0103] 6)定义主抽样的样本数化，对可靠性指标进行初始化，具体为LOLP(电力不足概 率）=0，EENS(电力电量不足期望值）=0;
[0104] 7)置模拟迭代次数n初值为0，即n = 0，置所有状态的状态变量为0，即H( 1: No) = 0; [010 引 8)n = n+l;
[0106] 9)利用新的元件强迫停运率VG和VL抽样系统的状态，对抽样的状态进行状态评估 分析，并根据是否发生失负荷更新可靠性指标。
[0107] 步骤9)具体包括状态抽样、状态评估和结果统计：
[0108] a)状态抽样:考虑所有的元件均为两状态模型，即对于第i个元件有：
[0109]
(2-13)
[0110] 0是一个[0，1]空间均匀分布的随机数，Vi是第i个元件的强迫停运率。
[0111] Xi是第i个元件的状态，1表示正常工作，0表示故障。一个综合系统中如果有输电 线路化条，发电机Ng台，其状态分别用Xli和拍1表示，则对发电机进行抽样：
[0112]
(2-14)
[01 1 3] VGi是第i个发电机的强迫停运率。
[0114]对输电线路的状态进行抽样：
[0115]
焊-：1 匀
[0116] VLi是第i条线路的强迫停运率。
[0117] b)状态评估:状态评估主要是对抽样的系统状态进行可靠性评估，主要是分析研 究该状态是否导致切负荷。
[0118] 如图2所示，状态评估的主要流程如下：
[0119] bl首先判断线路是否停运(故障），若有线路停运，则进入解列分析模块;若无线路 停运，则输出结果，结束评估。
[0120] b2分析系统是否解列，若系统解列，则在各个解列块内进行功率平衡分析，统计各 解列块内的失负荷量A P，累加各解列块的负荷损失量得到第n次的失负荷量(负荷损失量） A Pn，结束评估;如果没有解列，则进入潮流分析模块b3。
[0121] b3判断系统是否有孤立节点，有孤立节点，则去除孤立节点再形成导纳矩阵，然后 计算线路潮流;若无过负荷，结束评估，若有过负荷，进入状态校正环节评估是否有失负荷， 统计第n次系统的失负荷量APn,结束评估。其中的状态校正模型如下：
[012。
（2-峭
[0123] 式中i是节点的编号，D是节点的负荷功率向量，P是节点的发电机注入功率向量，N 是节点总数，是各节点的负荷需求向量，pmax是各发电机的最大发电量向量，PTDF是功率 分布因子矩阵，示是线路潮流容量上限向量，g是线路潮流容量下限向量。Pi和di分别是各 节点的注入功率和负荷功率，是决策变量。该状态校正模型的目标函数是系统总的失负荷 量最小。
[0124] C)结果统计:运里主要计算用于统计此时的可靠性评估得到的指标。具体如下：
[0125] 如果负荷损失A Pn大于0,按照下式分别更新LOLP和邸NS指标，具体计算如下：
[0126]
(2-17)
[0127] 式中的W(Xn;U，V)是似然率比，在重要抽样中用于指标的修正。具体按下式计算：
[012 引
(2-lS)
[0129] 式中W(Xn;u，v)是第n次抽样得到的状态的似然率比。j表示元件编号，Ng表示发电 机数目，化表示线路数目，Xnj表示第n次抽样得到的系统状态中第j个元件的状态，如果该变 量为1，表示该元件正常工作；否则，表示该元件故障。UW是第j条线路的原始强迫停运率， VW是第j条线路的新的强迫停运率。UU是第j条线路的原始强迫停运率，VU是第j条线路的 新的强迫停运率。
[0130] 10)判断n是否到达化，如果没有到达，转向步骤8)，否则结束计算。
[0131] 算例分析
[0132] 采用标准IE邸-RTS 79测试系统作为算例说明，该系统有33台发电机，38条线路， 下面的预抽样尺寸分别取为3000、4000、5000、6000。主抽样尺寸取为10000。并将6?-￥1结果 与CE(交叉赌是目前一种主流的重要抽样可靠性评估方法)结果进行对比，COV是收敛方差 系数，该系数越小，表明收敛越快，计算精度越高，越接近真实的可靠性指标，具体评估结果 表1所示：
[0133] 表1 GP-VM和CE可靠性评估结果对比 [01341
[0135] 一般而言，对于实际系统，GP-VM的预抽样样本数取到25000,主抽样样本依据需要 的指标精度来定，一般取为100000。
【主权项】
1. 基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统可靠性评估方法，其 特征在于:包括以下步骤： 1) 预抽样阶段 1.1) 确定预抽样阶段需要的参数，包括系统的峰值负荷水平、线路容量、线路原始强迫 停运率UL、发电机组发电额定容量以及发电机组原始强迫停运率U C，定义预抽样的样本数 No； 1.2) 置模拟迭代次数k=0; 1.3) 令k=k+1;然后抽样系统状态Xk，令系统中的元件为两状态模型，记录第k次抽样元 件的状态，元件的状态由发电机组状态Xc和线路状态Xl组成，对元件的状态进行评估分析后 计算系统是否会失负荷，然后统计评估分析的结果，如果失负荷，则H(X k) = I;否则，H(Xk) = 〇;H( ·)表示数组； 1.4) 判断k是否达到No，如果没有达到，则转至步骤1.3)，如果达到，则利用几何优化求 解基于重要抽样的电力系统可靠性评估的方差最小化模型，从而得到新的发电机组强迫停 运率 VG和线路强迫停运率VL; 2) 主抽样阶段 2.1) 定义主抽样的样本数N1，对可靠性指标进行初始化，并置模拟迭代次数n = 0; 2.2) 令n = n+l;然后利用新的元件强迫停运率％和1抽样系统的状态，对抽样的状态进 行状态评估分析，并根据是否发生失负荷更新可靠性指标； 2.3) 判断η是否达到见，如果没有达到，则转至步骤2.2)，否则，结束可靠性评估。2. 根据权利要求1所述基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统 可靠性评估方法，其特征在于:所述利用几何优化求解基于重要抽样的电力系统可靠性评 估的方差最小化模型具体包括以下步骤:将建立的基于重要抽样的电力系统可靠性评估的 方差最小化模型松弛为满足几何优化模型的几何优化-方差最小化模型，然后求解该模型， 得到新的重要抽样参数，即新的元件强迫停运率。3. 根据权利要求2所述基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统 可靠性评估方法，其特征在于:所述基于重要抽样的电力系统可靠性评估的方差最小化模 型的建立具体包括以下步骤:按照原始强迫停运率预抽样No次后，根据得到的评估结果按 以下公式形成基于重要抽样的电力系统可靠性评估的方差最小化模型：其中，Ng表示系统中发电机组的数目，Nl表示系统中线路的数目。4. 根据权利要求3所述基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统 可靠性评估方法，其特征在于:所述几何优化-方差最小化模型的获取方法包括以下步骤：，其中1?=1，2，···，Ν〇;引入虚拟变 量b和以，满足浐1-4和$广1-匕,，得到以下的优化模型：5. 根据权利要求2所述基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统 可靠性评估方法，其特征在于：通过对数变换和内点法求解所述几何优化-方差最小化模 型。6. 根据权利要求1所述基于几何优化-最小方差化方法的重要抽样蒙特卡洛电力系统 可靠性评估方法，其特征在于:所述可靠性指标包括电力不足概率LOLP和电力电量不足期 望值EENS，可靠性指标按照以下公式进行更新： L0LP = L0LP+lXff(Xn；u,v) EENS = EENs+Δ PnXff(Xn；u,v) 其中，△ Pn表示第η次抽样导致的负荷损失;W(Xn;u，v)表示第η次抽样得到的系统状态 的似然率比，该似然率比表示如下：其中，Ng表示系统中发电机组的数目，Nl表示系统中线路的数目。
【文档编号】G06Q50/06GK106022581SQ201610316648
【公开日】2016年10月12日
【申请日】2016年5月12日
【发明人】别朝红, 严超, 丁涛, 王灿, 胡源
【申请人】西安交通大学