基于增强型时序模式注意力的PMU数据恢复方法及系统

基于增强型时序模式注意力的pmu数据恢复方法及系统
技术领域
1.本发明涉及同步相量测量装置(phasor measurement unit，pmu)的数据恢复技术领域，具体涉及一种基于增强型时序模式注意力的pmu数据恢复方法及系统。


背景技术：

2.同步相量测量装置是电网广域测量系统的一种重要监测工具，其采集的同步相量数据质量对电网势态感知和动态监测方面起着关键性作用。由于电网的扩张和新型电力电子设备的接入，通信基础设施传输的数据量正在急剧增加。因此，受通信问题、网络攻击与电磁干扰等会出现数据丢失等问题。同步相量测量单元(phasor measurement unit，pmu)因其同步性、快速性和准确性成为电网动态实时监测最有效的手段之一，并得到迅速发展。但由于受通信问题、网络攻击与电磁干扰等诸多因素影响，pmu数据存在不同程度的缺失问题。pmu数据丢失直接影响其在电力系统中的应用，甚至威胁到系统的安全稳定运行。
3.总体而言，目前的pmu数据恢复方法存在以下几个缺点：1，pmu数据丢失存在随机性，且实际量测数据易受到系统噪声、谐波、振荡等影响，因此如何建立一种精准且高效的方法是pmu数据恢复的难点。2，pmu丢失的数据类型丰富，包括系统正常运行和系统处于故障状态的幅值、相角、频率数据，因此对适应于不同类型数据的特征提取方法提出了更高的要求。3，pmu采集的数据量庞大，因此所用方法应适用于大数据处理，并且有好的实时性。总而言之，现有方法对pmu数据恢复能在一定程度上实现量测数据的恢复，但恢复精确度和速度有待提升。


技术实现要素：

4.本发明要解决的技术问题：针对现有技术的上述问题，提供一种基于增强型时序模式注意力的pmu数据恢复方法及系统，本发明旨在提高挖掘pmu数据样本的数据间相关性的能力，实现更精确的pmu数据恢复。
5.为了解决上述技术问题，本发明采用的技术方案为：一种基于增强型时序模式注意力的pmu数据恢复方法，包括：s101，针对待恢复的pmu数据采用最优变分模态分解提取特征；s102，将特征归一化后输入结合增强型时序模式注意力和长短期记忆的网络模型，在结合增强型时序模式注意力和长短期记忆的网络模型中的处理包括：利用长短期记忆网络针对输入的特征计算隐状态向量以获取隐状态向量，采用滑动窗口获取lstm网络的隐状态向量构建行列的隐状态矩阵，其中是滑动窗口长度；是隐状态特征数；对隐状态矩阵用个大小的卷积核做卷积得到行列的时序模式矩阵；对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵；根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权
重，并对权重进行归一化得到注意力权重；基于注意力权重对增强时序模式矩阵的每一行进行加权求和得到加权和；将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出，并将输出输入全连接层得到恢复后的pmu数据。
6.可选地，步骤s102中对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵包括：s201，对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量；其中空间范围(,1)中的是指矩阵的行数，空间范围(1,)中的是指矩阵的列数；s202，将张量和张量融合成大小的张量并进行卷积、利用非线性激励函数得到用于对水平方向和垂直方向上的空间信息进行编码的中间特征图；s203，将中间特征图沿着空间维度拆分为张量和；将张量和分别用函数激活得到张量和，再将张量和进行张量相乘得到增强时序模式矩阵。
7.可选地，步骤s201中对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量的函数表达式为：，，上式中，为张量的第个通道的分量，为第个通道第行第个分量，为张量的第个通道的分量，为第个通道第列第个分量，为矩阵的列数，为矩阵的行数。
8.可选地，步骤s202中生成中间特征图的函数表达式为：
上式中，表示将张量和张量级联，w,b分别表示一个滤波器核的权值和偏置项，relu表示relu激励函数。
9.可选地，步骤s203中将张量和进行张量相乘得到增强时序模式矩阵的函数表达式为：，上式中，为增强时序模式矩阵的第行列元素，为张量的第个通道第行第个分量，为张量的第个通道第列第个分量。
10.可选地，步骤s102中根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权重的函数表达式为：，上式中，表示强时序模式矩阵中的第行，为权重矩阵，为当前时刻的lstm网络的隐状态向量；步骤s102中对权重进行归一化得到注意力权重的函数表达式为：，上式中，表示函数。
11.可选地，步骤s102中基于注意力权重对增强特征的每一行进行加权求和得到加权和的函数表达式为：，上式中，为隐状态矩阵的行数。
12.可选地，步骤s102中将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出的函数表达式为：，=，上式中，和为融合矩阵权重，为输出矩阵权重。
13.此外，本发明还提供一种基于增强型时序模式注意力的pmu数据恢复系统，包括相互连接的微处理器和存储器，所述微处理器被编程或配置以执行所述基于增强型时序模式注意力的pmu数据恢复方法。
14.此外，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序，所述计算机程序用于被微处理器编程或配置以执行所述基于增强型时序模式注意力的pmu数据恢复方法。
15.和现有技术相比，本发明主要具有下述优点：本发明提出了一种最优变分模态（ovmd）分解-增强型时序模式注意力（enhanced temporal pattern attention，etpa）-长短期记忆（long short-term memory，lstm）结合的pmu数据恢复技术，通过将上述技术的结
合，在时间模式注意力机制的基础上，利用得到的时间模式矩阵进行不同空间维度池化、加权变换，能够提高挖掘pmu数据样本的数据间相关性的能力，通过较少参数的增加来实现恢复精度大幅度的提升，实现更精确的pmu数据恢复。
附图说明
16.图1为本发明实施例方法的基本流程示意图。
17.图2为本发明实施例中etpa-lstm网络模型的原理示意图。
18.图3为本发明实施例中增强时序模式矩阵的变换生成原理示意图。
19.图4为本发明实施例方法和现有方法对静态数据进行恢复的实验结果对比。
20.图5为本发明实施例方法和现有方法对动态数据进行恢复的实验结果对比。
21.图6为本发明实施例方法和现有方法的静态数据恢复误差实验结果对比。
22.图7为本发明实施例方法和现有方法的动态数据恢复误差实验结果对比。
具体实施方式
23.如图1所示，本实施例基于增强型时序模式注意力的pmu数据恢复方法包括：s101，针对待恢复的pmu数据采用最优变分模态分解（optimal variational mode decomposition，ovmd）提取特征；s102，将特征归一化后输入结合增强型时序模式注意力（enhanced temporal pattern attention，etpa）和长短期记忆（long short-term memory，lstm）的网络模型（本实施例中记为etpa-lstm网络模型），在结合增强型时序模式注意力和长短期记忆的网络模型中的处理包括：利用长短期记忆网络针对输入的特征计算隐状态向量以获取隐状态向量，采用滑动窗口获取lstm网络的隐状态向量构建行列的隐状态矩阵，其中是滑动窗口长度；是隐状态特征数；对隐状态矩阵用个大小的卷积核做卷积得到行列的时序模式矩阵；对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵；根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权重，并对权重进行归一化得到注意力权重；基于注意力权重对增强时序模式矩阵的每一行进行加权求和得到加权和；将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出，并将输出输入全连接层得到恢复后的pmu数据。
24.本实施例步骤s101中采用的最优变分模态分解（ovmd）是对变分模态分解（vmd）的改进。对于给定的信号，变分模态分解（vmd）通过求解约束变分问题将数据分解为多个模态函数之和，各模态函数信号具有有限带宽，约束变分问题描述如下：
，（1）式中，为分解得到的模态总数，与分别为分解得到第个模态的时域信号和中心频率，为时刻的第个模态的时域信号，表示的偏导数，表示脉冲函数，为虚数单位，为给定的信号。为求解上述约束变分优化问题，引入二次惩罚因子和lagrange乘子，其中二次惩罚因子能有效降低高斯噪声干扰，而lagrange乘子则能增强约束严格性，得到增广变分问题：，（2）上式中，为lagrange乘子，为时刻的lagrange乘子。
25.利用基于对偶分解和lagrange法的交替方向乘子方法求解式（2）的增广变分问题，对、与进行交替迭代寻优，可得如下迭代公式：，（3），（4），（5）上式中，为第个模态的时域信号傅里叶变化的次迭代结果，为频率，为待分解信号的傅里叶变化，为第个模态的时域信号的傅里叶变化，为lagrange乘子的傅里叶变化，为第个模态中心频率的次迭代结果，为第个模态的时域信号的傅里叶变化，为lagrange乘子傅里叶变化的次迭代结果，为lagrange乘子傅里叶变化的次迭代结果，为步长。
26.对于给定的求解精度，满足下式时停止迭代：，（6）上式中，为第个模态的时域信号的次迭代结果，为第个模态
的时域信号的次迭代结果。
27.变分模态分解(vmd)的步骤包括：a）初始化、、与最大迭代次数，迭代次数；b）根据式（3）、（4）更新、；c）根据式（5）更新，；d）根据式（6）判断收敛性，若不收敛且，则跳转步骤b），否则停止迭代，得到最终模态函数和中心频率。
28.变分模态分解（vmd）的性能主要受模态总数和步长的影响。为了提高pmu数据恢复的精度，最优变分模态分解（ovmd）是对变分模态分解（vmd）增加了基于模态总数和步长双参数优化，其中模态总数采用基于中心频率观察法确定，步长采用残差的均方根误差最小化确定。双参数优化的步骤包括：第一步：初始化，，其中表示正整数集。
29.第二步：计算并分析不同值下分解模式的。若，此时的模态总数被确定为用于分解的最优的模态总数。
30.第三步：设置步长的取值范围，如式（7）所示。
31.，（7）上式中，为因变量步长的自变量，表示定义域从0-100，为自然数集，表示定义域从0-100的离散的101个点。
32.第四步：将信号，最优的模态总数，步长从小到大依次代入式（8）进行计算：，（8）上式中，~分别为第1~个模态的时域信号，ovmd为在确定最优的模态总数和对步长寻优过程中的最优变分模态分解。
33.第五步：将每个对应步长分解的信号进行重构，表达式为：，（9）上式中，为重构信号。
34.第六步：计算重构信号与原始信号的均方根误差(root mean squared error，rmse)，表达式如下：，（10）
上式中，为分解信号的长度，为第个重构信号，为第个原始信号。
35.第七步：选出rmse最小时对应的最优的步长。
36.在此基础上，即可利用最优的模态总数和最优的步长，即可针对pmu数据进行ovmd分解，完成特征提取。
37.本实施例的etpa-lstm网络模型中利用长短期记忆网络（lstm网络）针对输入的特征计算隐状态向量以获取隐状态向量。长短期记忆(long short-term memory，lstm)网络，与全连接神经网络相比，lstm隐藏层的值不仅由该时刻的输入层决定，还由上时刻的隐藏层决定。同时，lstm网络加入门结构，以控制每一时刻信息的记忆与遗忘，有选择地加以保留。记忆模块是lstm区别于一般神经网络的核心。lstm记忆模块包括记忆细胞和三个门控单元：输入门、输出门和遗忘门。其中，输入门决定信息是否输入至记忆细胞，遗忘门决定记忆细胞的信息是否删除，输出门决定记忆细胞的信息是否输出。lstm的三个门均采用函数，激励函数采用tanh函数。
38.本实施例的etpa-lstm网络模型中，增强型时序模式注意力（enhanced temporal pattern attention，etpa）是利用lstm网络得到隐藏状态矩阵，对该矩阵先用卷积提取特征行向量，然后对特征行向量进行不同空间维度的变换。由于特征行向量包含特征在多时间步上的时序信息，因此对特征行向量进行进一步变换，使增强型时间模式注意力机制(enhanced temporal pattern attention，etpa)能捕获到特征时间序列的更长期的相关性，以此来提高恢复精度。etpa-lstm网络模型的实现包括下述步骤：第一步：lstm网络中前一时刻的隐状态向量和当前输入信息同时传递给遗忘门，遗忘门的计算函数表达式为：，（11），（12）第二步：lstm网络的输入门决定输入记忆单元的哪些信息应该加入中，输入门的计算函数表达式为：，（13），（14），（15），（16）第三步：lstm网络的输出门判断哪些信息输出到时刻的隐藏层状态。输出门的计算函数表达式为：，（17），（18）式（11）-（18）中：为时刻的元胞状态；、、和分别为遗忘门、输入
门、元胞状态和输出门的输入权值矩阵；、、和分别为遗忘门、输入门、元胞状态和输出门的输出权值矩阵；、、和分别为遗忘门、输入门、元胞状态和输出门的偏置；、、和分别为遗忘门、输入门、元胞状态和输出门的输出向量；和为中间产生的过渡向量；为激活函数。
39.第四步：求隐状态矩阵。通过lstm网络每个时刻输出一个隐状态向量，个组合得到隐状态矩阵，大小为行列。隐状态矩阵的列向量表示时间步，行向量表示隐状态特征。其中是滑动窗口长度；是隐状态特征数。
40.第五步，卷积提取特征。对隐状态矩阵用个1
ⅹ
大小的卷积核从上往下沿着个特征做卷积得到时序模式矩阵，大小为行列，时序模式矩阵的行向量是特征行向量，列向量表示卷积核。这一步进行卷积操作可以提取隐状态矩阵的时序模式信息，增强模型的学习能力。时序模式矩阵的计算函数表达式为：，（19）上式中，为时序模式矩阵的第行列的元素，为隐状态矩阵中的第行列的元素，其中为总的时间步长，为代表滑窗的窗口长度，为与隐状态矩阵所包含的矩阵个数的相关参数，为第个卷积核的第个权重和偏置项，为第个卷积核，为隐状态矩阵的第行。
41.第六步，对时序模式矩阵进行变换，通过此变换，可以实现时序模式矩阵跨通道的交互和信息整合，进一步地提取时序模式矩阵的特征，这些特征能有效地提高pmu数据恢复的精度，其变换过程如图3所示。如图3所示，步骤s102中对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵包括：s201，对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量；其中空间范围(,1)中的是指矩阵的行数，空间范围(1,)中的是指矩阵的列数；s202，将张量和张量融合成大小的张量并进行卷积、利用非线性激励函数得到用于对水平方向和垂直方向上的空间信息进行编码的中间特征图；s203，将中间特征图沿着空间维度拆分为张量和；将张量和分别用
函数激活得到张量和，再将张量和进行张量相乘得到增强时序模式矩阵。
42.本实施例步骤s201中对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量的函数表达式为：，（20），（21）上式中，为张量的第个通道的分量，为第个通道第行第个分量，为张量的第个通道的分量，为第个通道第列第个分量，为矩阵的列数，为矩阵的行数。
43.本实施例步骤s202中生成中间特征图的函数表达式为：，（22）上式中，表示将张量和张量级联，w,b分别表示一个滤波器核的权值和偏置项，relu表示relu激励函数。
44.本实施例步骤s203中将张量和分别用函数激活得到张量和可表示为：，（23），（24）本实施例步骤s203中将张量和进行张量相乘得到增强时序模式矩阵的函数表达式为：，（25）上式中，为增强时序模式矩阵的第行列元素，为张量的第个通道第行第个分量，为张量的第个通道第列第个分量。
45.第七步，加权求和。根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权重，并对权重进行归一化得到注意力权重；基于注意力权重对增强时序模式矩阵的每一行进行加权求和得到加权和；本实施例中，步骤s102中根据增强时序模式矩阵和当前时刻的lstm网络的隐
状态向量确定增强时序模式矩阵中每一行的权重的函数表达式为：，（26）上式中，表示强时序模式矩阵中的第行，为权重矩阵，为当前时刻的lstm网络的隐状态向量；步骤s102中对权重进行归一化得到注意力权重的函数表达式为：，（27）上式中，表示函数。
46.紧接着，基于注意力权重对增强特征的每一行进行加权求和得到加权和，加权和是特征行向量的加权和，而特征行向量包含到多个时间步上的信息，因此它能捕获到特征时间序列的更长期的相关性。本实施例中，步骤s102中基于注意力权重对增强特征的每一行进行加权求和得到加权和的函数表达式为：，（28）上式中，为隐状态矩阵的行数。
47.第八步，融合函数。将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出。
48.具体地本实施例步骤s102中将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出的函数表达式为：，（29）=，（30）上式中，和为融合矩阵权重，为输出矩阵权重。
49.最终，通过全连接层即可输出恢复后的pmu数据。
50.为了验证本实施例基于增强型时序模式注意力的pmu数据恢复方法的性能，本实施例中将ovmd-tpa-lstm方法和tpa-lstm方法作为本实施例方法（ovmd-etpa-lstm方法）的对比，采用了静态和动态仿真数据进行实验。其中ovmd-tpa-lstm方法为ovmd+基于时序模式注意力机制（tpa）+长短时记忆（lstm）网络的方法，tpa-lstm方法即为进一步精简了ovmd的对比方法。动态数据采用基于wecc-179节点模型仿真的振荡案例作为测试案例。此案例的pmu上送频率为，即60s内可上送1800个同步相量数据。选取该案例连续180s的同步相量数据，即5400个同步相量数据作为实验样本。静态数据采用如下表达式生成5400个数据点作为实验样本：，（31）上式中，为pmu数据（静态数据），为相量幅值；为工频；为初相角；
为噪声部分；且本实施例中；；；的信噪比为60db。
51.首先确定静态和动态数据的ovmd分解的最优模态总数和最优步长；其次构建数据集，时间段的数据作为ovmd-etpa-lstm网络的输入，时刻的同步相量数据作为输出，滑动窗口为一个数据点，即下一个样本的输入为时间段的数据，输出为时刻的数据。假设连续丢失20个数据点，丢失时间段为179.333s-180s，分别采用tpa-lstm、ovmd-tpa-lstm和本实施例方法对丢失数据进行恢复，实验结果如图4和图5所示，绝对值误差如下图6和图7所示。从图4中可见，tpa-lstm方法的电压数据恢复值是一条直线，而采用了ovmd进行特征提取的ovmd-tpa-lstm方法和本实施例方法（ovmd-etpa-lstm方法）的恢复值与原始数据有相同的趋势，说明了ovmd进行特征提取对pmu数据恢复的重要性。从图4和图5可见，相较于tpa-lstm方法和ovmd-tpa-lstm方法，本实施例方法（ovmd-etpa-lstm方法）都更趋近于原始数据；同时从图6和图7中可见，本实施例方法（ovmd-etpa-lstm方法）对于静态和动态数据的恢复，其绝对值误差都明显小于tpa-lstm方法和ovmd-tpa-lstm方法。由此可以说明本实施例方法（ovmd-etpa-lstm方法）对pmu数据恢复更具优势。
52.为了更加精准的对不同算法的结果进行对比，采用综合矢量误差(total vector error，tve)来衡量恢复精度，各方法综合矢量误差tve如表1所示，各方法数据恢复时间如表2所示。
53.表1：不同方法的综合矢量误差tve对比。
54.数据类型tpa-lstmovmd-tpa-lstm本实施例方法静态数据1.20
×
10-52.57
×
10-61.55
×
10-6动态数据1.32
×
10-41.04
×
10-44.29
×
10-5表2：不同方法的数据恢复时间对比。
55.数据类型tpa-lstmovmd-tpa-lstm本实施例方法静态数据0.057s0.068s0.106s动态数据0.062s0.073s0.109s从表1可知，对于静态和动态数据恢复的综合矢量误差，本实施例方法（ovmd-etpa-lstm方法）都明显小于tpa-lstm方法和ovmd-tpa-lstm方法。由此说明本实施例方法（ovmd-etpa-lstm方法）的恢复精度更优。从表2可知，由于本实施例方法（ovmd-etpa-lstm方法）对时间模式矩阵进行了变换，该过程增加了少量参数，使得恢复时间有一定提高，但该方法在时间的整体数量级别上仍有较好的实时性。可见，本实施例方法（ovmd-etpa-lstm方法）在增加少量参数的情况下，能实现数据恢复精度的大幅度提高。
56.此外，本实施例还提供一种基于增强型时序模式注意力的pmu数据恢复系统，包括相互连接的微处理器和存储器，所述微处理器被编程或配置以执行所述基于增强型时序模式注意力的pmu数据恢复方法。此外，本实施例还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序，所述计算机程序用于被微处理器编程或配置以执行所述基于增强型时序模式注意力的pmu数据恢复方法。
57.本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可读存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。本技术是参照根据本技术实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
58.以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

技术特征：
1.一种基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，包括：s101，针对待恢复的pmu数据采用最优变分模态分解提取特征；s102，将特征归一化后输入结合增强型时序模式注意力和长短期记忆的网络模型，在结合增强型时序模式注意力和长短期记忆的网络模型中的处理包括：利用长短期记忆网络针对输入的特征计算隐状态向量以获取隐状态向量，采用滑动窗口获取lstm网络的隐状态向量构建行列的隐状态矩阵，其中是滑动窗口长度；是隐状态特征数；对隐状态矩阵用个大小的卷积核做卷积得到行列的时序模式矩阵；对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵；根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权重，并对权重进行归一化得到注意力权重；基于注意力权重对增强时序模式矩阵的每一行进行加权求和得到加权和；将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出，并将输出输入全连接层得到恢复后的pmu数据。2.根据权利要求1所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s102中对时序模式矩阵分别沿水平方向池化得到水平方向张量、沿竖直方向池化得到竖直方向张量，将水平方向张量、竖直方向张量相乘得到增强时序模式矩阵包括：s201，对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量；其中空间范围(,1)中的是指矩阵的行数，空间范围(1,)中的是指矩阵的列数；s202，将张量和张量融合成大小的张量并进行卷积、利用非线性激励函数得到用于对水平方向和垂直方向上的空间信息进行编码的中间特征图；s203，将中间特征图沿着空间维度拆分为张量和；将张量和分别用函数激活得到张量和，再将张量和进行张量相乘得到增强时序模式矩阵。3.根据权利要求2所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s201中对时序模式矩阵使用池化内核的空间范围(,1)沿水平坐标对每个通道进行编码得到大小的张量，对时序模式矩阵使用池化内核的空间范围(1,)沿垂直坐标对每个通道进行编码得到大小的张量的函数表达式为：
，，上式中，为张量的第c个通道的分量，为第c个通道第行第个分量，为张量的第c个通道的分量，为第c个通道第列第个分量，w为矩阵的列数，h为矩阵的行数。4.根据权利要求2所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s202中生成中间特征图的函数表达式为：上式中，表示将张量和张量级联，w,b分别表示一个滤波器核的权值和偏置项，relu表示relu激励函数。5.根据权利要求2所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s203中将张量和进行张量相乘得到增强时序模式矩阵的函数表达式为：，上式中，为增强时序模式矩阵的第行列元素，为张量的第个通道第行第个分量，为张量的第个通道第列第个分量。6.根据权利要求1所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s102中根据增强时序模式矩阵和当前时刻的lstm网络的隐状态向量确定增强时序模式矩阵中每一行的权重的函数表达式为：，上式中，表示强时序模式矩阵中的第行，为权重矩阵，为当前时刻的lstm网络的隐状态向量；步骤s102中对权重进行归一化得到注意力权重的函数表达式为：，上式中，表示函数。7.根据权利要求1所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s102中基于注意力权重对增强特征的每一行进行加权求和得到加权和的函数表达式为：，
上式中，为隐状态矩阵的行数。8.根据权利要求1所述的基于增强型时序模式注意力的pmu数据恢复方法，其特征在于，步骤s102中将加权和和当前时刻的lstm网络的隐状态向量进行加权求和得到加权后的隐状态向量，并根据隐状态向量计算出输出的函数表达式为：，=，上式中，和为融合矩阵权重，为输出矩阵权重。9.一种基于增强型时序模式注意力的pmu数据恢复系统，包括相互连接的微处理器和存储器，其特征在于，所述微处理器被编程或配置以执行权利要求1～8中任意一项所述基于增强型时序模式注意力的pmu数据恢复方法。10.一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序，其特征在于，所述计算机程序用于被微处理器编程或配置以执行权利要求1～8中任意一项所述基于增强型时序模式注意力的pmu数据恢复方法。

技术总结
本发明公开了一种基于增强型时序模式注意力的PMU数据恢复方法及系统，本发明方法包括针对待恢复的PMU数据采用最优变分模态分解提取特征再归一化后输入结合ETPA和LSTM的网络模型，利用LSTM获取隐状态向量并构建隐状态矩阵，对隐状态矩阵卷积得到时序模式矩阵并对该矩阵进行变换得到增强时序模式矩阵；基于注意力权重对增强时序模式矩阵的每一行进行加权求和得到加权和；将加权和和隐状态向量加权为隐状态向量并计算输出，将输出输入全连接层得到恢复后的PMU数据。本发明旨在提高挖掘PMU数据样本的数据间相关性的能力，实现更精确的PMU数据恢复。PMU数据恢复。PMU数据恢复。


技术研发人员：姚文轩 杨秋爽 邱伟 黄琴 马钰浩
受保护的技术使用者：湖南大学
技术研发日：2023.06.12
技术公布日：2023/7/20