结合人工特征和残差网络的电力系统暂态稳定预测方法

1.一种结合人工特征和残差网络的电力系统暂态稳定预测方法，其特征在于，该方法包括以下步骤：(1)对具有N台发电机的电力系统，根据电力系统历史运行情况和运行人员经验，对s种运行工况在f种故障下的暂态稳定性进行时域仿真计算，得到s×f种运行场景的特征向量X k 和暂态稳定性y k ，其中，上标k表示第k种运行场景，k＝1,2,…,s×f，y k ＝(0,1)表示电力系统在故障清除后能够保持暂态稳定，y k ＝(1,0)表示电力系统在故障清除后不能保持暂态稳定，故障清除时间根据人工经验设定，第k种运行场景中故障清除后n个采样点的发电机有功功率P Gi k 、发电机转子角度δ i k 、发电机转子角速度ω i k 、发电机母线的电压幅值V Gi k 和发电机母线的电压相角θ Gi k 构成特征向量X k ：X k ＝[P Gi k (t),δ i k (t),ω i k (t),V Gi k (t),θ Gi k (t)]其中，下标i表示电力系统中的第i台发电机，i＝1,2,…,N，t表示第t个采样点，t＝1,2,…,n，n为人为设定的采样点数，采样频率选为电力系统的额定频率；(2)根据步骤(1)中的特征向量X k ＝[P Gi k (t),δ i k (t),ω i k (t),V Gi k (t),θ Gi k (t)]，计算如下30×n个人为定义的人工特征：人工特征 其中0 - 表示故障发生前的最后一个采样值人工特征 人工特征 人工特征 人工特征Y 5 k (t)＝Y 3 k (t)-Y 4 k (t)人工特征Y 6 k (t)＝Y 2 k (t)/Y 1 k (t)人工特征 人工特征 人工特征 人工特征 人工特征Y 11 k (t)＝Y 9 k (t)-Y 10 k (t)人工特征Y 12 k (t)＝Y 8 k (t)/Y 7 k (t)人工特征 人工特征 人工特征 人工特征 人工特征Y 17 k (t)＝Y 15 k (t)-Y 16 k (t)人工特征Y 18 k (t)＝Y 14 k (t)/Y 13 k (t)人工特征 人工特征 人工特征 人工特征 人工特征Y 23 k (t)＝Y 21 k (t)-Y 22 k (t)人工特征Y 24 k (t)＝Y 20 k (t)/Y 19 k (t)人工特征 人工特征 人工特征 人工特征 人工特征Y 29 k (t)＝Y 27 k (t)-Y 28 k (t)人工特征Y 30 k (t)＝Y 26 k (t)/Y 25 k (t)将上述30×n个人工特征Y p k (t)进行最大最小归一化，其中，下标p＝1,…,30，得到归一化后的人工特征 归一化的公式为：(3)将每种场景的特征向量X k ＝[P Gi k (t),δ i k (t),ω i k (t),V Gi k (t),θ Gi k (t)]进行最大最小归一化，按照发电机维度、变量维度和时间维度排列成三维数据 然后设置卷积层、池化层、残差单元和全连接层，结合步骤(2)中的30×n个归一化后的人工特征 得到暂态稳定预测模型M的结构，具体包括以下步骤：(3-1)将步骤(1)得到的每一种运行场景下的特征向量X k ＝[P Gi k (t),δ i k (t),ω i k (t),V Gi k (t),θ Gi k (t)]进行最大最小归一化，归一化的公式为：然后，将归一化后的数据 按照发电机维度、时间维度和变量维度排列成三维数据 该三维数据 的维度为N×n×5；(3-2)结合步骤(2)中的人工特征和深度学习中的残差网络，设计得到暂态稳定预测模型M的结构，M的输入数据是步骤(2)得到的归一化后的人工特征 和步骤(3-1)得到的三维数据 M的输出是O k 2 ，当O k 2 ＝(0,1)时，表示电力系统在第k种运行场景下能够保持暂态稳定，当O k 2 ＝(1,0)时，表示电力系统在第k种运行场景下不能保持暂态稳定，模型M由如下多个单元堆叠而成：(3-2-1)卷积层：利用c个卷积核w l 和一个偏置矩阵V 0 对步骤(3-1)中第k种运行场景的三维数据 进行卷积操作，得到特征向量O k ，其中，l＝1,…,c，卷积核w l 和偏置矩阵V 0 是步骤(4)的待求参数，w l ∈R a×d ，R a×d 表示a×d维矩阵，且矩阵中每一个元素都是实数，a和d的值取为奇数，且满足a≤N，d≤n，卷积核的个数c≥5；(3-2-2)池化层：对特征向量O k 进行最大池化，获得池化后的特征A k 0 ；(3-2-3)堆叠的m个残差单元：利用深度残差网络中堆叠的m个残差单元，对步骤(3-2-3)的池化后特征A k 0 进行特征提取，其中，第f个残差单元的输出为：A k f ＝σ(σ(A k f-1 *J l f,1 +V f,1 )*J l f,2 +V f,2 +A k f-1 )其中，上标f＝1,…,m，m为残差单元的总数，m的值由人为设定，σ()是ReLU激活函数，A k f 是第f个残差单元的输出，A k f-1 是第f-1个残差单元的输出，J l f,1 是第f级残差单元使用的第一层卷积层的第l个卷积核，J l f,2 是第f级残差单元使用的第二层卷积层的第l个卷积核，l＝1,…,c，V f,1 是第f级残差单元使用的第一层卷积层的偏置矩阵，V f,2 是第f级残差单元使用的第二层卷积层的偏置矩阵，卷积核J l f,1 和J l f,2 、偏置矩阵V l f,1 和V l f,2 是步骤(4)的待求参数；(3-2-4)池化层：对第m级残差单元的输出A k m 进行最大池化，获得池化后特征Q k ；(3-2-5)批归一化层：利用批归一化方法，对步骤(3-2-4)的池化后特征Q k 进行归一化处理，得到归一化后的特征U k ；(3-2-6)平铺层：利用平铺函数，将步骤(3-2-5)的归一化特征U k 平铺为h×1维特征向量V k ，其中h的大小由归一化特征U k 的维度决定；(3-2-7)第一层全连接层：将步骤(2)得到的30×n个归一化后的人工特征 与步骤(3-2-6)的h×1维特征向量V k 合并为Z k ，Z k 是(h+30×n)×1维向量，然后将Z k 输入到第一层全连接层中，得到第一层全连接层的输出为O k 1 ：O k 1 ＝σ(GZ k +b 1 )其中，上标1表示第一层全连接层，权重矩阵G∈R g×h ，R g×h 表示g×h维矩阵，且矩阵中每一个元素都是实数，第一层全连接层的偏置向量b 1 ∈R g×1 ，R g×1 表示g×1维向量，且向量中每一个元素都是实数，g表示全连接层的输出维度，第一层全连接层的输出维度由人为设定，取值范围设为g∈(2,h+30×n)，权重矩阵G和偏置向量b 1 是步骤(4)的待求参数；(3-2-8)第二层全连接层：将步骤(3-2-7)的输出O k 1 输入到第二层全连接层中，得到第二层全连接层的输出为O k 2 ：O k 2 ＝Softmax(HO k 1 +b 2 )其中，上标2表示第二层全连接层，权重矩阵H∈R 2×g ，R 2×g 表示2×g维矩阵，且矩阵中每一个元素都是实数，第二层全连接层的偏置向量b 2 ∈R 2×1 ，R 2×1 表示2×1维向量，且向量中每一个元素都是实数，Softmax()是Softmax激活函数，权重矩阵H和偏置向量b 2 均为步骤(4)的待求参数；(4)根据步骤(1)得到的s×f个样本和基于自适应矩估计的梯度下降算法，即Adam算法，迭代计算M中的待求参数，得到最终的暂态稳定预测模型，具体包括如下步骤：(4-1)从步骤(1)得到的s×f个样本中随机抽取 个样本作为训练集，剩余 个样本作为验证集，其中 表示对0.8×s×f向下取整；(4-2)设定一个集合S＝{e max ,A max ,M max }，其中A max 是迭代过程中得到暂态稳定预测模型的最高预测准确率，e max 是取得最高预测准确率A max 时的迭代次数，M max 是第e max 次迭代得到的暂态稳定预测模型，记迭代次数为r，最大迭代次数为r max ，最小迭代次数为r min ，其中r max 和r min 的值由人为设定，并满足r max ＞r min ≥10，设迭代次数r的初值为0，e max 初值为0，A max 的初值为-1，模型M max 设为空；(4-3)将迭代次数r与最大迭代次数r max 进行比较：(4-3-1)若r≥r max ，则集合S中的M max 作为最终的电力系统暂态稳定预测模型；(4-3-2)若r＜r max ，则令r:＝r+1，转入步骤(4-4)；(4-4)利用步骤(4-1)的训练集和Adam算法，计算步骤(3)中模型M的所有待求参数，包括w l 、V 0 、J l f,1 、J l f,2 、V l f,1 、V l f,2 、G、H、b 1 和b 2 ，得到第r次迭代当前参数对应的暂态稳定预测模型M r ；(4-5)利用M r 对步骤(4-1)验证集中所有样本的暂态稳定性进行预测，得到预测准确率，记为A r ，将A r 的值与A max 进行比较：(4-5-1)若A r >A max ，则令e max ＝r，A max ＝A r ，M max ＝M r ，更新得到新的集合S，然后转入步骤(4-3)；(4-5-2)若A r ≤A max ，则将e max 的值与r和r-r min 的值进行比较，如果满足r-r min ≤e max ≤r，则转入步骤(4-3)，如果不满足r-r min ≤e max ≤r，则停止迭代，取集合S中的M max 作为最终的电力系统暂态稳定预测模型；(5)获取故障切除后发电机的有功功率、转子角度、转子角速度、电压幅值和电压相角，经过计算并输入至步骤(4)得到的电力系统暂态稳定预测模型中，得到暂态稳定预测结果，具体包括如下步骤：(5-1)利用离线时域仿真计算或直接采集电力系统广域测量信息系统的量测数据，得到故障清除后n个采样点的发电机有功功率P Gi (t)、发电机转子角度δ i (t)、发电机转子角速度ω i (t)、发电机母线的电压幅值V Gi (t)和发电机母线的电压相角θ Gi (t)，构成初始输入特征，其中t＝1,…,n；(5-2)对步骤(5-1)的初始输入特征计算得到步骤(2)定义的30×n个归一化后的人工特征；(5-3)利用步骤(3-1)的最大最小归一化对初始输入特征进行归一化处理后并排列成N×n×5的三维数据；(5-4)将步骤(5-2)得到的30×n个归一化后的人工特征和步骤(5-3)得到的N×n×5的三维数据共同输入至步骤(4)得到的电力系统暂态稳定预测模型中，得到电力系统的暂态稳定预测结果。