CN109376939A - 一种基于自适应神经网络的电网稳定性实时预测方法 - Google Patents

Abstract

一种基于自适应人工神经网络的发电机转子角预测方法，包括以下步骤：采用数据自动生成算法生成数据并存储在数据库中；将数据库中数据进行规范化；通过反向传播算法训练数据；最后通过对发电机转子角的稳定性预测来实现电网实时预测。本发明为电力***发生故障时实时预测电网稳定性提供了决策支持，对提高大规模电网中长期电压稳定裕度，改善电网运行的经济性与质量，均有着重大意义。

Description

一种基于自适应神经网络的电网稳定性实时预测方法

技术领域

本发明属于电力***安全分析评估领域，具体涉及一种基于自适应神经网络的电网稳定性实时预测方法。

背景技术

暂态稳定性是指***在不失去稳定性的情况下，能够承受从发电到负载的最大潮流转移的能力。此外，当网络条件发生突然变化时，如故障或负载突然增加时，就会出现这种现象，在这种情况下，***受到发电机转子角的较大偏移，并不受非线性功率角关系的影响。影响***稳定性的最重要参数是初始运行状态和故障的严重程度，而扰动前的稳态与扰动后的稳态不同，因此需要通过转子角度的变化进行动态分析来检查其稳定性。

发电机作为电力***中的重要组成部分，其故障会给电力***带来严重危害，有必要早期预警以防止事故的发生或进一步扩大。尤其是同步发电机的转子绕组由许多串联线圈组成，而转子具有较大的转动惯量，由于离心力的作用，在运行中线匝绝缘的移动，转子绕组端部的热变形，线匝端部垫块松动或护环绝缘衬垫老化，小的导电粒子或碎物进入转子线匝端部和转子通风沟导致转子绕组匝间短路发生。

在线动态安全评估仍然是大多数电网关注的主要问题，因为负载过重的发电机在长传输互连链路上运行，导致负阻尼、振荡和失去同步。在这样的***中，需要进行长期的模型分析来检测这些意外事件的影响，而对发电机转子角的稳定性预测存在很多局限性，且准确度不理想。因此应用神经网络解决电力***动态安全评价问题具有很大意义，现有方法主要存在以下缺陷：①作为神经网络的训练样本，都应用暂态能量函数进行稳定分析的结果，而这种分析方法受***数学模型详细程度和计算结果准确度的限制，使训练样本本身就不够准确，从而影响神经网络的训练结果；②在进行动态安全评价时，通常的已知量为***在故障前的运行参数，它们来自状态估计或在线潮流的计算结果，而要获得故障发生后的暂态量，将需要大量的附加计算，从而使它们不便于在线应用；③传统神经网络算法每次迭代的步长效率较低，从而导致算法收敛速度较慢，神经网络的灵敏度也相对不高，从而会影响神经网络应用程序的性能。

发明内容

针对上述问题与不足，本发明提供一种基于自适应神经网络的电网稳定性实时预测方法，该方法利用自动数据生成与自适应神经网络算法，训练出规范化数据并计算其输出输入的误差，从而对电力***稳定高效的在线动态安全评估，并预测出发电机转子角的稳定性，从而达到对电网稳定性实时预测的作用。该方法有效的处理在电力***安全评估以及预测转子角稳定性中所遇到的问题，有利于解决电网中许多安全问题，提高了电网安全运行的水平。

为了实现上述发明目的，本发明采取如下技术方案：

一种基于自适应神经网络的电网稳定性实时预测方法，包括以下步骤：

步骤1：采用动态安全评估(Dynamic Security Assessment，DSA)数据生成算法选择各发电机的初始条件、负载需求和故障编号作为神经网络预测***的输入参数，将这些数据存储在数据库中；

步骤2：将数据库里的数据转换成神经函数可以处理的形式，需要在自适应人工神经网络(Adaptive Artificial Neural Network，AANN)开始时进行规范化处理；

步骤3：将规范化的数据进行训练，通过反向传播算法，计算出各个数据输出值，并估计误差，达到训练目的；

步骤4：在神经网络完成了训练过程并达到目标后，进行发电机转子角的稳定性预测。

进一步地，在步骤1中，提出一种数据自动生成的算法，根据发电机之间的差异，采用自动数据生成模型确定了***中所有电机的转子角稳定性。如果角差达到峰值，然后减小，***就处于稳定状态。如果任何角度的差异无限增加，会导致至少一台机器失去同步，***就会在一个不稳定的条件下工作。

在此基础上，必须为稳定和不稳定条件选择目标，当***在扰动后稳定时，输出值设置为0，当***不稳定时，输出值设置为1。重复此过程，直到***发生故障前达到最大负载水平，就可以将数据进行存储。

进一步地，在所述步骤2中具体包含以下步骤：

步骤2-1：对输入的矢量数据进行归一化，以便将所有数据x_i都放在-1到1之间，以防止网络权重的任何波动，归一化的过程基于：

r＝x_i(max)-x_i(min) (2)

其中：x_i为输入的数据，x_i(max)、x_i(min)为输入值的最大值和最小值，r是输入值的取值范围，r₀是所有输入值的中点值。

步骤2-2：提出了一种优化方法来处理由于隐藏层中没有指定数目的神经元这一问题。

通过讨论输入与输出的电导率，可以得出输入层的神经元数目与输入的神经元数目相同，输出层的神经元数目与输出的神经元数目相同。此外，影响人工神经网络(Artificial Neural Network，ANN)输出的另外两个因素是：一是学习率系数，它改变了权重调整的大小；二是动量项，添加动量项可以提高收敛速度，当x_i为常数时，输出值被认为是零。

进一步地，在所述步骤3中包括如下步骤：

步骤3-1：网络体系结构第一层中，输入被应用到它们的输出目标上，并通过网络层传播来计算误差之和E_sum，通过使用均方根误差(Root Mean Square Errors，RMSE)来计算该误差之和，均方根误差定义如下所示：

其中：N为事故的总数，t_O为目标值，O_O为AANN输出值，P＝(1，2，...，N)。

步骤3-2：网络体系结构第二层中，有两个停止时刻的标准：第一个时刻是当一个阈值误差达到一个性能目标值时会停止，第二个时刻是当训练达到一定的迭代次数，网络不能达到阈值误差值时会停止。第二层的网络输出计算定义为：

net_k＝∑w_kiO_i+b_k (5)

其中：i为输入节点，k为隐藏节点，O为输出节点，O_i为第一层的输出，w_ki和w_ok是输入层和隐藏层之间，以及隐藏层和输出层之间的权重矩阵。将b_k设为1，加快收敛过程。

步骤3-3：网络体系结构第三层，网络输出的计算由以下定义：

net_O＝∑w_okO_k (6)

其中：w_ok是隐藏层和输出层之间的权重矩阵，O_k为隐藏层的输出。

进一步地，在所述步骤3-3中隐藏神经元和输出神经元采用的激活函数是双曲正切函数，激活函数由以下定义：

其中：f(net_k)为隐藏神经元的激活函数，f(net_O)为输出神经元的激活函数。

误差的计算定义为：

E＝t_O-O_O (9)

两个网络必须根据两个错误信号进行更新，通过两个错误信号更新，以获得更好的收敛性。

第一个错误信号δ_O在输出(O)和隐藏层(k)之间，δ_O定义为：

δ_O＝E×O_O(1-O_O) (10)

第二个错误信号δ_k在隐藏层(k)和输入层(i)之间，δ_k定义为：

δ_k＝O_k(1-O_k)∑δ_Ow_Ok (11)

进一步地，在所述步骤4中，一般情况下，扰动后的暂态稳定现象持续时间为3-5秒，对于具有主要振荡的区域间模态的大***，扰动持续时间可延长至10秒，这个时间足以确定转子的振荡是否稳定。

***中所有机器的转子角都是根据发电机之间的差异来确定的，转子角定义为：

Δδ(g_i,g_k)＝δg_i-δg_k (12)

其中：Δδ为转子角，g_i、g_k是发电机数量，i＝1，2，...(n-1)；k＝1，2，...n，n为监测的发电机总数。

当扰动后转子的振荡稳定时，输出值为0；不稳定时输出值为1，此时将发电机记入不稳定列表中。

当输出为1，不稳定列表中出现发电机时，***显示出不稳定状态，此时计算信号N_S：

当AANN处于训练时，新的操作条件遵循以下定义：

f(x)-S＝0 (14)

其中：x为输入参数，f(x)为用于适应新操作条件的静态训练，S为每个情况选择的安全状态。

独立动力***具有一般的非线性微分关系，该微分关系定义为：

f'(x,Δ)＝f'(x)-Δ×S' (15)

其中：S’为转子角的预测安全状态，在AANN训练完成后，神经网络会跟随该式进行***的模拟，并得出预测结果。

与现有技术比，本发明一种基于自适应神经网络的电网稳定性实时预测方法，有益效果在于：

①、传统方法讨论了用时间步长预测方法预测***行为的神经网络应用，然而该方法不依赖于集成时间步长进行安全预测，并且新的操作条件自动地适合于所训练的网络。

②、本发明采用了一种DSA自动数据生成方法，能快速的自动生成所需数据。DSA模型能够根据转子角度、负载水平和扰动类型的初始条件对发电机的稳定状态进行分类和预测，该方法的响应速度快，可以对***的每个序列进行稳定性评估，该方法十分便于在线应用。

③、本发明在利用传播算法的同时还增加了RMSE方程的应用，RMSE方程在误差达到最小且迭代次数更少的情况下，还能提高神经网络的灵敏度，这些因素增强了神经网络应用程序的性能，使其比标准的反向传播方法更稳定。

附图说明

图1是本发明***流程图；

图2是本发明DSA数据自动生成算法框图；

图3是本发明实例中不稳定发电机总数对分类误差的影响情况的散点图。

具体实施方式

下面结合附图对本发明作进一步说明。

本发明提供一种基于自适应神经网络的电网稳定性实时预测方法，所述方法包括以下步骤：

图1中步骤1：采用动态安全评估DSA数据生成算法选择各发电机的初始条件、负载需求和故障编号作为神经网络预测***的输入参数，将这些数据存储在数据库中。

在该步骤中提出一种数据自动生成的算法，根据发电机之间的差异，采用自动数据生成模型确定了***中所有电机的转子角稳定性。如果角差达到峰值，然后减小，***就处于稳定状态。如果任何角度的差异无限增加，会导致至少一台机器失去同步，***就会在一个不稳定的条件下工作。

其中，该数据自动生成算法如图2所示。

图1中步骤2：将数据库里的数据转换成神经函数可以处理的形式，需要在AANN开始时进行规范化处理。

所述步骤2包括以下步骤：

步骤2-1：对输入的矢量数据进行归一化，以便将所有数据x_i都放在-1到1之间，以防止网络权重的任何波动，归一化的过程基于：

r＝x_i(max)-x_i(min) (2)

其中：x_i为输入的数据，x_i(max)、x_i(min)为输入值的最大值和最小值，r是输入值的取值范围，r₀是所有输入值的中点值。

步骤2-2：提出了一种优化方法来处理由于隐藏层中没有指定数目的神经元这一问题。

通过讨论输入与输出的电导率，可以得出输入层的神经元数目与输入的神经元数目相同，输出层的神经元数目与输出的神经元数目相同。此外，影响ANN输出的另外两个因素是：一是学习率系数，它改变了权重调整的大小；二是动量项，添加动量项可以提高收敛速度，当x_i为常数时，输出值被认为是零。

图1中步骤3：将规范化的数据进行训练，通过反向传播算法，计算出各个数据输出值，并估计误差，达到训练目的。

所述步骤3包括以下步骤：

步骤3-1：网络体系结构第一层中，输入被应用到它们的输出目标上，并通过网络层传播来计算误差之和E_sum，通过使用均方根误差(RMSE)来计算该误差之和，均方根误差定义如下所示：

其中：N为事故的总数，t_O为目标值，O_O为AANN输出值，P＝(1，2，...，N)。

步骤3-2：网络体系结构第二层中，有两个停止时刻的标准：第一个时刻是当一个阈值误差达到一个性能目标值时会停止，第二个时刻是当训练达到一定的迭代次数，网络不能达到阈值误差值时会停止。第二层的网络输出计算定义为：

net_k＝∑w_kiO_i+b_k (5)

其中：i为输入节点，k为隐藏节点，O为输出节点，O_i为第一层的输出，w_ki和w_ok是输入层和隐藏层之间，以及隐藏层和输出层之间的权重矩阵。将b_k设为1，加快收敛过程。

步骤3-3：网络体系结构第三层，网络输出的计算由以下定义：

net_O＝∑w_okO_k (6)

其中：w_ok是隐藏层和输出层之间的权重矩阵，O_k为隐藏层的输出。

在所述步骤3-3中隐藏神经元和输出神经元采用的激活函数是双曲正切函数，激活函数由以下定义：

其中：f(net_k)为隐藏神经元的激活函数，f(net_O)为输出神经元的激活函数。

误差的计算定义为：

E＝t_O-O_O (9)

两个网络必须根据两个错误信号进行更新，通过两个错误信号更新，以获得更好的收敛性。

第一个错误信号δ_O在输出(O)和隐藏层(k)之间，δ_O定义为：

δ_O＝E×O_O(1-O_O) (10)

第二个错误信号δ_k在隐藏层(k)和输入层(i)之间，δ_k定义为：

δ_k＝O_k(1-O_k)∑δ_Ow_Ok (11)

图1中步骤4：在神经网络完成了训练过程并达到目标后，进行发电机转子角的稳定性预测。

在所述步骤4中，一般情况下，扰动后的暂态稳定现象持续时间为3-5秒，对于具有主要振荡的区域间模态的大***，扰动持续时间可延长至10秒，这个时间足以确定转子的振荡是否稳定。

***中所有机器的转子角都是根据发电机之间的差异来确定的，转子角定义为：

Δδ(g_i,g_k)＝δg_i-δg_k (12)

其中：Δδ为转子角，g_i、g_k是发电机数量，i＝1，2，...(n-1)；k＝1，2，...n，n为监测的发电机总数。

当扰动后转子的振荡稳定时，输出值为0；不稳定时输出值为1，此时将发电机记入不稳定列表中。

当输出为1，不稳定列表中出现发电机时，***显示出不稳定状态，此时计算信号N_S：

当AANN处于训练时，新的操作条件遵循以下定义：

f(x)-S＝0 (14)

其中：x为输入参数，f(x)为用于适应新操作条件的静态训练，S为每个情况选择的安全状态。

独立动力***具有一般的非线性微分关系，该微分关系定义为：

f'(x,Δ)＝f'(x)-Δ×S' (15)

其中：S’为转子角的预测安全状态，在AANN训练完成后，神经网络会跟随该式进行***的模拟，并得出预测结果。

实施例：

本发明提出的方法在一个87节点算例***中进行了测试，该***分为四个区域，即北部、东部、中部和南部。该***包含171条传输线及23台同步发电机。

对该***采用基于此发明的发电机转子角稳定性预测方法，为了检验使用自适应神经网络应用工具进行动态安全评估的能力，对神经网络进行了不同负载水平下的测试针对不同的负载水平生成了三组测试数据集，以检验AANN的精度和性能，得出此方法AANN的精度和性能比传统方法高的多。

图3验证了不稳定发电机总数对误差分类错误率的影响，可以看出发电机不稳定数量和分类误差率的上升比例趋于稳定，并且在有近50台不稳定发电机的情况下错误率也能维持在百分之五以下。

此外，当许多发电机不稳定时，分类误差会增加；因此不稳定的发电机应该分成一组，其余的发电机应该分成另一组，以保持安全的运行状态。

通过对训练后的神经网络进行RMSE检验，验证了该方法的有效性。见如下表1：

表1

区域	事故数量	均方根误差(％)	T1	T2
					北部	63	2.370	0.140	6.599
东部	27	0.463	0.022	0.554
					中部	53	0.021	0.061	3.213
南部	47	0.635	0.074	3.478

T1：每个事故的平均训练时间(秒)；T2：所有事故的训练时间(秒)。

可以看出，AANN可以在6.6秒内预测最大区域的稳定性状态，说明本发明的AANN工具通过转子角稳定性模拟能够在极短的时间内准确地预测***的状态。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。

Claims (10)

1.一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于包括以下步骤：

步骤1：采用动态安全评估数据生成算法，选择各发电机的初始条件、负载需求和故障编号作为神经网络预测***的输入参数，将这些数据存储在数据库中；

步骤2：将数据库里的数据转换成神经函数可以处理的形式，需要在自适应人工神经网络开始时进行规范化处理；

步骤3：将规范化的数据进行训练，通过反向传播算法，计算出各个数据输出值，并估计误差，达到训练目的；

步骤4：在神经网络完成了训练过程并达到目标后，进行发电机转子角的稳定性预测。

2.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤1中，提出一种数据自动生成的算法，根据发电机之间的差异，采用自动数据生成模型确定了***中所有电机的转子角稳定性，如果角差达到峰值，然后减小，***就处于稳定状态；如果任何角度的差异无限增加，会导致至少一台机器失去同步，***就会在一个不稳定的条件下工作。

3.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤1中，为稳定和不稳定条件选择目标；当***在扰动后稳定时，输出值设置为0，当***不稳定时，输出值设置为1；重复此过程，直到***发生故障前达到最大负载水平，就将数据进行存储。

4.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤2中具体包含以下步骤：

步骤2-1：对输入的矢量数据进行归一化，以便将所有数据x_i都放在-1到1之间，以防止网络权重的任何波动，归一化的过程基于：

r＝x_i(max)-x_i(min) (2)

其中：x_i为输入的数据，x_i(max)、x_i(min)为输入值的最大值和最小值，r是输入值的取值范围，r₀是所有输入值的中点值；

步骤2-2：提出了一种优化方法来处理由于隐藏层中没有指定数目的神经元这一问题；

通过讨论输入与输出的电导率，可以得出输入层的神经元数目与输入的神经元数目相同，输出层的神经元数目与输出的神经元数目相同；此外，影响人工神经网络输出的另外两个因素是：一是学习率系数，它改变了权重调整的大小；二是动量项，添加动量项可以提高收敛速度，当x_i为常数时，输出值被认为是零。

5.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤3中，在训练过程中采用反向传播算法，反向传播算法与感知器网络算法相似，有多个层组成；该算法中网络体系结构由三层组成：第一层与输入连接；第二层包含激活函数；第三层是网络的输出。

6.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤3中具体包含以下步骤：

步骤3-1：网络体系结构第一层，输入被应用到它们的输出目标上，并通过网络层传播来计算误差之和E_sum，通过使用均方根误差来计算该误差之和(Root Mean Square Errors，RMSE)，均方根误差定义如下所示：

其中：N为事故的总数，t_O为目标值，O_O为AANN输出值，P＝(1，2，...，N)；

步骤3-2：网络体系结构第二层，有两个停止时刻：第一个时刻是当一个阈值误差达到性能目标值时；第二个时刻是当训练达到一定的迭代次数，网络不能达到阈值误差值时；第二层的网络输出计算定义为：

net_k＝∑w_kiO_i+b_k (5)

其中：i为输入节点，k为隐藏节点，O为输出节点，O_i为第一层的输出，w_ki和w_ok是输入层和隐藏层之间，以及隐藏层和输出层之间的权重矩阵；将b_k设为1，加快收敛过程；

步骤3-3：网络体系结构第三层，网络输出的计算由以下定义：

net_O＝∑w_okO_k (6)

其中：w_ok是隐藏层和输出层之间的权重矩阵，O_k为隐藏层的输出。

7.根据权利要求1所述一种基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤3-3中隐藏神经元和输出神经元采用的激活函数是双曲正切函数，该激活函数由以下定义：

其中：f(net_k)为隐藏神经元的激活函数；f(net_O)为输出神经元的激活函数。

8.根据权利要求6所述的基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤3-3中误差的计算，即输出值与目标值进行比较，该误差计算定义为：

E＝t_O-O_O (9)。

9.根据权利要求6所述的基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤3-3中两个网络必须根据两个错误信号进行更新，通过两个错误信号更新，以获得更好的收敛性；

第一个错误信号δ_O在输出(O)和隐藏层(k)之间，δ_O定义为：

δ_O＝E×O_O(1-O_O) (10)

第二个错误信号δ_k在隐藏层(k)和输入层(i)之间，δ_k定义为：

δ_k＝O_k(1-O_k)∑δ_Ow_Ok (11)。

10.根据权利要求1所述的基于自适应神经网络的电网稳定性实时预测方法，其特征在于：在所述步骤4中，一般情况下，扰动后的暂态稳定现象持续时间为3-5秒，对于具有主要振荡的区域间模态的大***，扰动持续时间可延长至10秒，这个时间足以确定转子的振荡是否稳定；

***中所有机器的转子角都是根据发电机之间的差异来确定的，转子角定义为：

Δδ(g_i,g_k)＝δg_i-δg_k (12)

其中：Δδ为转子角，g_i、g_k是发电机数量，i＝1，2，...(n-1)；k＝1，2，...n，n为监测的发电机总数；

当扰动后转子的振荡稳定时，输出值为0；不稳定时输出值为1，此时将发电机记入不稳定列表中；

当输出为1，不稳定列表中出现发电机时，***显示出不稳定状态，此时计算信号N_S：

当AANN处于训练时，新的操作条件遵循以下定义：

f(x)-S＝0 (14)

其中：x为输入参数，f(x)为用于适应新操作条件的静态训练，S为每个情况选择的安全状态；

独立动力***具有一般的非线性微分关系，该微分关系定义为：

f'(x,Δ)＝f'(x)-Δ×S' (15)

其中：S’为转子角的预测安全状态,在AANN训练完成后，神经网络会跟随该式进行***的模拟，并得出预测结果。