CN103544546A - 一种在线预测电力系统静态电压稳定极限的方法 - Google Patents

Abstract

本发明公开了一种在线预测电力系统静态电压稳定极限的方法，该方法由基于Lasso的样本降维方法、基于自组织特征映射网络的样本筛选方法和误差反向传播型神经网络三部分组成。本发明通过自组织特征映射网络和Lasso方法对电力系统离线仿真生成的大量随机训练样本进行预处理，Lasso方法用于减少输入变量的维数，自组织特征映射网络用于减少训练样本的数量，之后采用误差反向传播型神经网络对电力系统静态电压稳定极限进行离线训练和在线预测，本发明的方法可以有效提高误差反向传播型神经网络的离线训练效率和在线预测效果。

Description

一种在线预测电力系统静态电压稳定极限的方法

技术领域

本发明属于电力系统稳定性分析技术领域，具体涉及一种在线预测电力系统静态电压稳定极限的方法。

背景技术

而随着现代电力系统的发展，系统的规模不断扩大，网络结构日趋复杂，环境和经济因素的制约使电力系统的运行情况更加接近极限条件。在线评估电压稳定性是预防大停电事故的重要手段之一。

快速求解系统静态电压稳定极限是在线评估电力系统电压稳定性最基本的要求。电力系统在运行过程中，随着电力系统中的负荷不断增加，系统的潮流雅克比矩阵将趋向于奇异，当潮流雅克比奇异时，系统达到潮流临界解，该状态下系统的负荷总有功功率即为系统的静态电压稳定极限（或称极限传输功率、极限负荷能力），系统静态电压稳定极限越大，则可认为系统在该情景下的静态电压稳定性越好。连续潮流法是计算电力系统静态电压稳定极限的经典方法，但由于其计算量较大、耗时较多，难于在线应用。由于电力系统静态电压稳定极限无法用解析表达式描述，因此该领域研究主要集中在采用神经网络对静态电压稳定极限进行在线估计。

李妍、程时杰等在标题为电压稳定极限附近神经网络模拟及静态稳定裕度判定（电力系统自动化，1999(21)，37-40）的文献中提出了一种用误差反向传播神经网络来模拟节点电压静态稳定临界状态附近的过渡过程的方法，该方法中可用于分析在某种运行方式电网各节点的电压静态稳定储备系数。崔峰、齐占庆、姜萌在标题为基于模糊神经网络的电力系统电压稳定评估（电力系统保护与控制，2009(11)，40-44）的文献中设计了一个多输入单输出的模糊神经网络，通过计算带静止无功补偿器电力系统的负荷能力极限，对系统的电压稳定性进行评估，有功与无功负载的不确定性、有功与无功发生器、母线电压以及无功补偿装置的参数都被考虑进来，该方法采用Kohonen自组织映射网络聚类所有节点的有功和无功负荷来减少输入量个数。陈爱军、刘爱国在标题为基于神经网络模型的在线电压稳定极限评估（电力科学与工程，2010(08)，19-23）的文献将回归分析和人工神经网络模型同时应用于电力系统电压稳定评估，用回归分析法求取负荷对电压稳定裕度的灵敏度，根据预先设置好的灵敏度阀值来进行特征选择，从而减少输入变量的维数。

上述三种采用神经网络来离线训练和在线估计静态电压稳定极限时存在以下两个问题：（1）电力系统可用的输入变量过多；（2）训练样本数量多但分布并不理想。这两个问题将降低神经网络的离线训练效率并影响其预测效果，在系统规模较大时其影响十分明显。

发明内容

针对现有技术所存在的上述技术问题，本发明提供了一种在线预测电力系统静态电压稳定极限的方法，能够有效提高神经网络的离线训练效率和在线预测效果。

一种在线预测电力系统静态电压稳定极限的方法，包括如下步骤：

（1）通过对电力系统进行离线仿真，生成系统在随机多个工况下对应的多个训练样本；所述的训练样本包括输入变量和输出变量，所述的输出变量为系统的静态电压稳定极限（又称为负荷能力极限、极限传输功率）；

每个样本的系统的静态电压稳定极限的求取方法如下：电网在某一随机工况下，不断增加全网负荷和发电机的功率，负荷功率增长方式为全网各负荷保持恒功率因数按各自的初始有功功率比例同步增长，发电机功率增长方式为除平衡机以外的全体发电机按各自的初始有功出力比例共同承担全网增长的负荷总有功功率。随着全网负荷和发电机功率的不断增加，系统的潮流雅克比矩阵趋向奇异，当系统达到潮流临界解时，该状态下系统的全网负荷总有功功率即为系统的静态电压稳定极限。

（2）对所述的训练样本进行降维，然后通过自组织特征映射网络对所有训练样本进行分类筛选；

（3）利用筛选保留下来的训练样本对误差反向传播型神经网络进行训练，得到用于预测电力系统静态电压稳定极限的神经网络模型；

（4）实时采集电力系统的状态信息，通过所述的神经网络模型计算出系统的静态电压稳定极限。

所述的输入变量和状态信息均包括系统的全网负荷总有功功率、全网负荷总无功功率、全网发电机总有功出力、全网发电机总无功出力、系统中各节点的节点电压和节点相角、各发电机的有功出力和无功出力、各负荷的有功功率和无功功率。

所述的步骤（2）中对训练样本进行降维的具体过程如下：

A1.通过Lasso（Least Absolute Shrinkage and Select Operator，最小绝对值收缩选择）方法确定训练样本输入变量中各特征量对应的回归系数；

A2.对于任一训练样本，删除其输入变量中回归系数为0对应的特征量；

A3.根据步骤A2，遍历所有训练样本。

所述的步骤（2）中对训练样本进行分类筛选的具体过程如下：

B1.利用自组织特征映射网络对所有训练样本进行分类；

B2.对于任一类训练样本，计算该类所有训练样本关于静态电压稳定极限的均方差，若求得的均方差小于预设的阈值，则保留该类中静态电压稳定极限最大和最小的两个训练样本并删除其他训练样本；否则，保留该类所有训练样本；

B3.根据步骤B2，遍历每一类训练样本。

所述的步骤（3）中，采用单隐含层的误差反向传播型神经网络进行训练。

本发明的方法通过自组织特征映射网络和Lasso方法用于对离线训练样本集的输入数据进行预处理，采用误差反向传播型神经网络对电力系统静态电压稳定极限进行离线训练和在线预测，可以有效提高误差反向传播型神经网络的离线训练效率和在线预测效果。

附图说明

图1为本发明方法的流程示意图。

图2为新英格兰电力试验系统的单线结构示意图。

图3为测试样本的误差分布示意图。

具体实施方式

为了更为具体地描述本发明，下面结合附图及具体实施方式对本发明的技术方案进行详细说明。

如图1所示，一种在线预测电力系统静态电压稳定极限的方法，包括如下步骤：

步骤一：

通过离线仿真产生训练和测试所需的样本；测试样本的产生方式和训练样本的产生方式完全相同。

本实施方式的电力系统为10机39节点新英格兰电力试验系统，如图2所示，该系统包含有39个节点（其中19个负荷节点）、10台发电机，基本运行状态下，系统总负荷为6097.1MW，其余的随机初始运行状态的获取方法如下：各负荷在基准值的±30%以内各自随机变化（负荷功率因数保持不变），各发电机的有功出力在基准值的±30%以内各自随机变化，各发电机的机端电压在基准值的±3%以内各自随机变化。负荷功率增长方式为全网各负荷保持恒功率因数按各自的初始有功功率比例同步增长，发电机功率增长方式为除平衡机以外的全体发电机按各自的初始有功出力比例共同承担全网增长的负荷总有功功率，考虑发电机的无功限制。本实施共采样3000组训练样本和500组测试样本。每个样本需要记录的状态量包括：全网负荷总有功功率、全网负荷总无功功率、全网发电机总有功出力、全网发电机总无功出力、各节点电压、各节点相角、各发电机有功出力、各发电机无功出力、各负荷有功功率、各负荷有功功率无功功率，共139维输入变量和1维输出变量。

步骤二：

基于Lasso方法的对训练样本维数进行降阶：

A1.根据训练样本通过Lasso方法确定训练样本输入变量中各特征量对应的回归系数。

Lasso方法是Robert Tibshirani于1996年提出的变量选择技术，Lasso方法用模型系数的绝对值函数做为惩罚来压缩模型系数，使一些回归系数变小，甚至还使一些绝对值较小的系数直接变为0，这就使得这种方法同时据有了子集选择和岭回归的优点。与传统的变量选择方法相比，Lasso方法很好的克服了传统方法在选择变量上的不足，因此该方法在统计领域受到了极大的重视。标准Lasso的目标函数和约束条件如下：

$\begin{array}{ccc}\underset{b}{\min }{\left|\right|\mathrm{Ab}-y\left|\right|}_2^2& s.t.\underset{i=1}{\overset{n_1}{\mathrm{\Σ}}}\left|b_i\right|\≤t& A=\left[\begin{array}{c}{x}_1\\ .\\ .\\ .\\ x_{3000}\end{array}\right]\end{array}$

上式的目标函数中，矩阵A为所有训练样本的输入变量的集合，x_i表示第i个训练样本的输入变量（1×139维行向量），i为自然数且1≤i≤3000，A的维数为3000×139，向量b为每个输入变量对应的回归系数，即b=[b₁，b₂，…，b₃₀₀₀]^T，向量y为所有训练样本的输出变量的集合，即y=[y₁，y₂，…，y₃₀₀₀]^T，目标函数的意义是通过确定向量b的取值使表达式的2-范数最小；约束条件中b_i为向量b的第i个元素，t为调和参数，约束条件的意义为调和参数t的控制使回归系数总体变小。

对于Lasso的计算，主要就是调和参数t的确定和解二次规划问题。确定调和参数t的方法有交叉核实、广义交叉核实等。确定调和参数t后，对Lasso的数学模型进行二次规划问题求解，首先把带有绝对值的不等式约束化条件转化为简单的线性约束，通过公开文献披露的引入顺序不等式约束的方法可以找到满足龙格—库塔条件的可行解。

A2.对于任一训练样本，删除其输入变量中回归系数为0对应的特征量；

通过求解标准Lasso方法的数学模型，可以得到回归系数向量b，b中零元素对应的输入状态量对系统静态电压稳定极限的影响非常小，因此在后续的误差反向传播型网络训练时可以删去，反之，b中非零元素对应的输入状态量应予以保留。

A3.根据步骤A2，遍历所有训练样本。

本实施例中每个训练样本共有139维输入变量和1维输出变量，本实施例中Lasso方法的惩罚系数λ取50，算法的终止判据为迭代次数达到10000次或收敛判据小于1×10^-7。Lasso方法解出的回归系数向量b中含零元素61个，其对应的系统状态量在后续误差反向传播型神经网络的训练样本和预测样本中可以删去，共剩余78维系统状态量作为后续误差反向传播型神经网络的输入。

步骤三：

基于自组织特征映射网络的对训练样本进行分类筛选：

B1.利用自组织特征映射网络对所有训练样本进行分类；

通过自组织特征映射网络对所有训练样本进行分类。自组织特征映射网络是由Tuevo Kohonen在1979年和1982年提出和发展的，自组织特征映射网络采用无监督的学习方式，能够根据网络的输入模式而组织自身去学习其内在的相似性并进行分类，采用最小不相似性度量来确定胜出单元。对于同一电网的任意两个运行状态，在两者的静态电压稳定极限和运行状态量均很接近时，则可认为这两个运行状态属于同一类。本方法采用自组织特征映射网络对系统运行状态进行分类，本实施例中自组织特征映射网络的输入样本组成如下： $X_{s\left(i\right)}=$ $[P_{G1\left(i\right)}^{*},...,P_{\mathrm{Gp}\left(i\right)}^{*},V_{1\left(i\right)}^{*},...,V_{\mathrm{nb}\left(i\right)}^{*},P_{L1\left(i\right)}^{*},...,P_{\mathrm{Lq}\left(i\right)}^{*},Q_{L1\left(i\right)}^{*},...,Q_{\mathrm{Lq}\left(i\right)}^{*},\mathrm{\ω}{P}_{M\left(i\right)}^{*}],$ x_s(i)为第i个输入样本，p为电网中发电机母线的数量，nb为电网中母线的数量，q为电网中负荷母线的数量，x_s(i)中变量的上标“*”表示该变量已经归一化。

为第i个训练样本对应系统工况下第k个节点的发电机有功功率，

为第i个训练样本对应系统工况下第k个节点的节点电压，

为第i个训练样本对应系统工况下第k个节点的负荷有功功率，为第i个训练样本对应系统工况下第k个节点的负荷无功功率，ω为静态电压稳定极限的权重系数，

为第i个训练样本对应系统工况下的系统静态电压稳定极限。i为自然数，本实施例中1≤i≤3000。

B2.对于任一类训练样本，计算该类所有训练样本关于静态电压稳定极限的均方差，若求得的均方差小于预设的阈值，则保留该类中静态电压稳定极限最大和最小的两个训练样本并删除其他训练样本；否则，保留该类所有训练样本；

本实施例采用6×6的二维网格状自组织特征映射网络对样本进行分类，自组织特征映射网络隐含层的形状为六边形，神经元间的距离用欧氏距离计算，训练次数为20次。本实施例的自组织特征映射网络输入变量包括系统中各节点的节点电压、各负荷的有功功率和无功功率，共58维，静态电压稳定极限的权重系数ω取10，样本集合的均方差筛选阈值取1MW。采用3000组训练样本训练自组织特征映射网络，训练完成后样本被分为36类。

B3.根据步骤B2，遍历每一类训练样本。

步骤三中3000组训练样本被自组织特征映射网络分为36类，其中有9类样本集合(共含1179个训练样本)的均方差小于1MW，在这9类样本集合中各选择1个静态电压稳定极限最小和1个静态电压稳定极限最大的样本作为代表，从而用18个样本代替1179个样本作为后续误差反向传播型神经网络的训练样本。该步骤共减少样本数量1161个，占训练样本总数的38.7%。

步骤四：

用简化后的训练样本训练误差反向传播型神经网络。

误差反向传播型神经网络的训练样本为1839个78维的样本。本算例的误差反向传播型神经网络采用单隐含层结构，隐含层含10个神经元，活化函数采用Sigmoid函数，输入层和输出层均为线性Purelin函数，训练采取Levenberg-Marquardt算法，初始学习率为1，学习率增大因子为10，学习率减小因子为0.1，学习的终止条件为训练周期达到100次或性能指标小于1×10^-5。本发明方法的数据预处理和训练误差反向传播型神经网络的用时共为45.01s，而直接采用误差反向传播型神经网络训练初始3000个139维的训练样本需用时96.24s，可以看出本发明方法能够明显提高误差反向传播型神经网络的训练速度。

步骤五：

测试训练好的误差反向传播型神经网络是否满足在线预测静态电压稳定极限的要求。用以下两个指标测试误差反向传播型神经网络的预测性能，第一个性能指标为测试样本的平均相对误差，第二个性能指标为测试样本的最大相对误差。

第一个性能指标测试样本的平均相对误差e_Bmean的计算公式如下：

$e_{\mathrm{Bmean}}=\frac{1}{n_2}\underset{i=1}{\overset{n_2}{\mathrm{\Σ}}}\left|\frac{y_{\mathrm{BPi}}-y_i}{y_i}\right|$

第二个性能指标测试样本的最大相对误差e_Bmax的计算公式如下：

$e_{B\max }=\max \left|\frac{y_{\mathrm{BPi}}-y_i}{y_i}\right|,(i=\mathrm{1,2},...,n_2)$

式中：y_BPi为误差反向传播型神经网络对第i个测试样本的静态电压稳定极限预测值，y_i为第i个测试样本的静态电压稳定极限实际值，n₂为测试样本的数量。

通过500个测试样本对训练好的误差反向传播型神经网络进行测试，测试结果为：测试样本的平均相对误差为0.30%，最大相对误差为2.59%，测试样本的误差分布如图3所示。测试结果表明所有样本下静态电压稳定极限预测值和静态电压稳定极限实际值的误差均较小，能够满足在线预测静态电压稳定极限的要求。

Claims (5)

1.一种在线预测电力系统静态电压稳定极限的方法，包括如下步骤：

（1）通过对电力系统进行离线仿真，生成系统在随机多个工况下对应的多个训练样本；所述的训练样本包括输入变量和输出变量，所述的输出变量为系统的静态电压稳定极限；

（2）对所述的训练样本进行降维，然后通过自组织特征映射网络对所有训练样本进行分类筛选；

（3）利用筛选保留下来的训练样本对误差反向传播型神经网络进行训练，得到用于预测电力系统静态电压稳定极限的神经网络模型；

（4）实时采集电力系统的状态信息，通过所述的神经网络模型计算出系统的静态电压稳定极限。

2.根据权利要求1所述的方法，其特征在于：所述的输入变量和状态信息均包括系统的全网负荷总有功功率、全网负荷总无功功率、全网发电机总有功出力、全网发电机总无功出力、系统中各节点的节点电压和节点相角、各发电机的有功出力和无功出力、各负荷的有功功率和无功功率。

3.根据权利要求1所述的方法，其特征在于：所述的步骤（2）中对训练样本进行降维的具体过程如下：

A1.通过Lasso方法确定训练样本输入变量中各特征量对应的回归系数；

A2.对于任一训练样本，删除其输入变量中回归系数为0对应的特征量；

A3.根据步骤A2，遍历所有训练样本。

4.根据权利要求1所述的方法，其特征在于：所述的步骤（2）中对训练样本进行分类筛选的具体过程如下：

B1.利用自组织特征映射网络对所有训练样本进行分类；

B2.对于任一类训练样本，计算该类所有训练样本关于静态电压稳定极限的均方差，若求得的均方差小于预设的阈值，则保留该类中静态电压稳定极限最大和最小的两个训练样本并删除其他训练样本；否则，保留该类所有训练样本；

B3.根据步骤B2，遍历每一类训练样本。

5.根据权利要求1所述的方法，其特征在于：所述的步骤（3）中，采用单隐含层的误差反向传播型神经网络进行训练。

Images