计及时序相关性的数据驱动电-气综合能源系统态势感知方法与模型

1.本发明涉及综合能源系统态势感知技术领域，尤其指一种计及时序相关性的 数据驱动电-气综合能源系统态势感知方法。


背景技术：

2.电-气综合能源系统的多源异构特性为系统经济、高效、安全、可靠运行带 来了问题与挑战，如何对电-气综合能源系统进行高效准确态势感知，对于保障 系统安全可靠运行具有重大意义。
3.近年来，出于电力系统可靠性与安全性需求，国内外专家学者对态势感知技 术在电力领域的应用进行了探索与研究。lin jin等人对主动配电网中态势感知 的路线、技术与瓶颈进行了分析与展望。李大虎等人建立了元件级与系统级的运 行风险态势指标，对台风气象下的电网运行风险进行了评估与预测。肖白等人通 过降噪自编码器对电力负荷噪声进行处理，基于奇异谱分析方法对电力负荷高/ 低频分量进行分解，最后利用长短期记忆神经网络实现了空间电力负荷的准确感 知与预测。wang weikang等人提出了一种基于相角转移量与频率变化率的电力 系统频率扰动感知识别方法。jing hongyan等人提出了一种基于深度q网络的 电、热、冷负荷态势感知方法。刘柳等人针对以能量枢纽为多能供应源的区域综 合能源系统，提出了安全域概念与模型，实现了系统安全区间准确感知。以上研 究促进了态势感知技术在电力领域的推广与应用。但是电-气综合能源系统场景 下的态势感知技术研究仍较少且不够深入，缺乏对于海量多源异构运行量测数据 隐含信息的有效挖掘，为电-气综合能源系统的安全可靠运行带来了问题与挑战。


技术实现要素：

4.本发明所要解决的技术问题是提供一种计及时序相关性的数据驱动电-气综 合能源系统态势感知方法，该方法考虑到偏离量的时序相关性特征，并基于数据 驱动理论建模，可有效挖掘历史数据的隐含信息，提高系统态势感知准确率。
5.为了解决上述技术问题，本发明采用如下技术方法：一种计及时序相关性的 数据驱动电-气综合能源系统态势感知方法，包括以下步骤：
6.s1、采集数据：采集电-气综合能源系统中量测装置与设备的历史数据；
7.s2、处理数据：利用状态估计对步骤s1中采集的量测噪声滤波处理，进行 电-气综合能源系统状态量的求解并察觉系统偏离量，将得到的系统偏离量划分 为历史偏离与未来偏离，基于偏离量在时域上的变化情况得到历史时段偏离量对 应的特征图像以及未来时段偏离量对应的分类标签；
8.s3、构建并训练模型：构建卷积神经网络模型，使用步骤s2中得到的特征 图像和分类标签对所述卷积神经网络模型进行训练，以得到模型输入-输出的映 射关系；
9.s4、预测当前时段的未来运行趋势：采用步骤s1-s2中的方法采集并处理电
ꢀ‑
气综
合能源系统中量测装置与设备在当前时段的量测数据，以得到当前时段偏 离量对应的图像特征，将该图像特征输入训练好的卷积神经网络模型中，得到系 统未来运行趋势的分类标签。
10.进一步地，在步骤s1中，采集数据时：
11.采集电-气综合能源系统中相量测量单元、广域测量系统、压力变送计器、 涡轮流量计等量测装置与设备的量测数据，所述量测数据包括：电网电压幅值 ui、电网电压相角θ、电网节点有功功率pi、电网节点无功功率qi、电网线路有 功负荷p
ij
、电网线路无功负荷q
ij
、气网节点压强pi、气网节点流量fi、气网管 道流量f
ij
。
12.进一步地，在步骤s2中：
13.采用基于加权最小二乘法建立的电-气综合能源系统状态估计模型对量测噪 声滤波处理，该状态估计模型如下：
[0014][0015]
式中：ze与zg分别为电网与气网量测量向量；he与hg分别为电网与气网量 测模型；xe与xg分别为电网与气网状态量向量；re与rg分别为电网与气网量测协 方差矩阵；c(x)为零注入节点伪量测约束。
[0016]
采用电-气综合能源系统偏离量计算模型计算电-气综合能源系统偏离量，该 偏离量计算模型如下：
[0017][0018]
式中：d为电-气综合能源系统状态偏离量；ρ为权重系数；me与mg分 别为电网与气网节点数；ui与pj分别为电网节点i电压幅值与气网节点j压 强；un与pn分别为电网额定电压与气网额定压强。
[0019]
更进一步地，在步骤s2中，将得到的系统偏离量划分为历史偏离与未来偏 离时：假设所述系统偏离量为p～f时段内获取的系统偏离量，t为当前时刻， t∈(p,f)，则p～t时段偏离量d
p
～
t
为历史偏离，t～f时段偏离量d
t
～f为未来偏 离。
[0020]
更进一步地，在步骤s2中，利用格拉姆角差场理论对历史偏离量时序相关 性进行分析得到历史时段偏离量对应的特征图像，利用格拉姆角差场理论对未来 偏离量运行趋势进行定性分类分析得到未来时段偏离量对应的分类标签。
[0021]
再进一步地，在步骤s2中，得到历史时段偏离量对应的特征图像的具体方 法如下：
[0022]
1)p～t时段系统偏离量归一化：对p～t时段偏离量d
p～t
进行归一化操作， 将数据区间缩放至[-1,1]，如下式：
[0023][0024]
式中：为n时刻归一化系统偏离量；dn为n时刻系统偏离量；max(d
p～t
) 与min
(d
p～t
)分别为p～t时段系统偏离量最大值与最小值；
[0025]
2)极坐标系形式转换：将归一化处理后的系统状态偏离量数值通过反 余弦计算变换为角度，对应时间戳变换为半径，实现系统偏离从直角坐标系 到极坐标系的转换与表达，如下式：
[0026][0027]
式中：θn为偏离量在极坐标系中角度编码；rn为偏离量在极坐标系中半 径编码；
[0028]
3)格拉姆角差场矩阵生成：在得到偏离量时间序列的极坐标系形式后， 利用各时间断面偏离量角度差的正弦运算生成格拉姆角差场矩阵，以识别多 时间断面间时序自相关性，如下式：
[0029][0030]
式中：g为p～t时段偏离量dp～t对应格拉姆角差场矩阵；
[0031]
4)将格拉姆角差场变换所得格拉姆角差场矩阵g进行图像编码，得到 p～t时段偏离量dp～t对应的特征图像。
[0032]
再进一步地，在步骤s2中，得到未来时段偏离量对应的分类标签的具体方 法如下：
[0033]
依据未来系统偏离量变化将系统未来运行趋势划分为健康、恢复、临界、 紧急四种状态：
[0034]
1)判断max(d
t～f
)≤dh是否成立，其中，dh为健康状态偏离量阈值，若是， 则未来运行趋势为健康状态；若否，转到2)；
[0035]
2)判断dh＜max(d
t～f
)≤dc是否成立，其中，dc为临界状态偏离量阈值，若 是，转到3)；若否，转到4)；
[0036]
3)判断是否成立，若是，则未来运行趋势为临界状态；若 否，则未来运行趋势为恢复状态；
[0037]
4)判断是否成立，若是，则未来运行趋势为紧急状态；若 否，则未来运行趋势为恢复状态；
[0038]
经分类所得的未来运行趋势即为t～f时段偏离量dt～f对应的分类标签。
[0039]
优选地，在步骤s3中构建的卷积神经网络模型包括1层输入层、2层卷积 层、2层池化层、2层激活层、1层扁平层、1层全连接层与1层输出层，所述卷 积神经网络模型的卷积公式如下：
[0040]
[0041]
式中：x为卷积计算输入矩阵；w代表大小为m
×
n的卷积核；y为卷积计算 输出矩阵。
[0042]
优选地，在步骤s4中，将当前时段偏离量对应的图像特征输入至卷积神经 网络模型中以预测未来运行趋势时：
[0043]
所述卷积神经网络模型先利用其输入层输入当前时段格拉姆角差场图像并 将其转换为32
×
32
×
3的像素矩阵；再通过两次“卷积-激活-池化”处理得到5
×
5 大小的128个特征图，其中激活层均采用relu作为激活函数，池化层均采用最 大池化方式；接着利用扁平层将三维特征图压缩得到长度为3200一维数组；最 后通过全连接层将扁平层信息传递至输出层，利用softmax分类器得出当前时段 的未来运行趋势的分类标签。
[0044]
作为本发明的另一面，一种计及时序相关性的数据驱动电-气综合能源系统 态势感知模型，包括前述计及时序相关性的数据驱动电-气综合能源系统态势感 知方法中的卷积神经网络模型。
[0045]
本发明提供的计及时序相关性的数据驱动电-气综合能源系统态势感知方法 与模型是态势感知技术在多源异构数据系统中的一次突破，促进了态势感知技术 的发展，也实现了电-气综合能源系统运行态势感知准确率的提高。具体而言， 本发明通过巧妙的运用偏离量在时域上的自相关性，并基于数据驱动理论建立态 势预测模型，可有效挖掘电-气综合能源系统中隐含的信息，提高电-气综合能源 系统运行态势感知准确率，为电-气综合能源系统提供更高效准确的风险预警信 息，保障系统安全可靠的运行。
附图说明
[0046]
图1为本发明所涉计及时序相关性的数据驱动电-气综合能源系统态势感知 方法的流程图；
[0047]
图2是本发明所涉计及时序相关性的数据驱动电-气综合能源系统态势感知 方法中格拉姆角差场变化示意图；
[0048]
图3是本发明所涉计及时序相关性的数据驱动电-气综合能源系统态势感知 方法与模型中的卷积神经网络模型架构图；
[0049]
图4是本发明实施例中电-气综合能源系统气网节点1-7的压强状态数据图；
[0050]
图5是本发明实施例中电-气综合能源系统电网节点1-14的电压状态数据图。
具体实施方式
[0051]
为了便于本领域技术人员的理解，下面结合实施例与附图对本发明作进一步 的说明，实施方式提及的内容并非对本发明的限定。
[0052]
如图1所示，一种计及时序相关性的数据驱动电-气综合能源系统态势感知 方法，包括四大步骤，具体如下。
[0053]
步骤s1、采集历史量测数据；
[0054]
采集电-气综合能源系统中相量测量单元、广域测量系统、压力变送计器、 涡轮流量计等量测装置与设备在p～f时段的历史量测数据，该历史量测数据包 括：电网电压幅值ui、电网电压相角θ、电网节点有功功率pi、电网节点无功功 率qi、电网线路有功负荷p
ij
、电网线路无功负荷q
ij
、气网节点压强pi、气网节 点流量fi、气网管道流量f
ij
。
[0055]
步骤s2、处理历史量测数据；
[0056]
s21，采用基于加权最小二乘法建立的电-气综合能源系统状态估计模型对量 测噪声滤波处理，该状态估计模型如下：
[0057][0058]
式中：ze与zg分别为电网与气网量测量向量；he与hg分别为电网与气网量 测模型；xe与xg分别为电网与气网状态量向量；re与rg分别为电网与气网量测协 方差矩阵；c(x)为零注入节点伪量测约束。
[0059]
s22，采用电-气综合能源系统偏离量计算模型计算电-气综合能源系统在p～f 时段的偏离量，该偏离量计算模型如下：
[0060][0061]
式中：d为电-气综合能源系统状态偏离量；ρ为权重系数；me与mg分 别为电网与气网节点数；ui与pj分别为电网节点i电压幅值与气网节点j压 强；un与pn分别为电网额定电压与气网额定压强。
[0062]
s23，基于步骤s23获取了系统p～f时段偏离量，假设t为当前时刻，t∈(p,f)， 则电-气综合能源系统在p～t时段偏离量d
p
～
t
为历史偏离，在t～f时段偏离量 d
t
～f为未来偏离。
[0063]
s24，利用格拉姆角差场理论对历史偏离量时序相关性进行分析得到历史时 段偏离量对应的特征图像，利用格拉姆角差场理论对未来偏离量运行趋势进行定 性分类分析得到未来时段偏离量对应的分类标签，具体步骤如下：
[0064]
s241，获取历史时段偏离量对应的特征图像；
[0065]
1)如图2所示，p～t时段系统偏离量归一化：对p～t时段偏离量d
p～t
进 行归一化操作，将数据区间缩放至[-1,1]，如下式：
[0066][0067]
式中：为n时刻归一化系统偏离量；dn为n时刻系统偏离量；max(d
p～t
) 与min(d
p～t
)分别为p～t时段系统偏离量最大值与最小值；
[0068]
2)如图2所示，极坐标系形式转换：将归一化处理后的系统状态偏离 量数值通过反余弦计算变换为角度，对应时间戳变换为半径，实现系统偏离 从直角坐标系到极坐标系的转换与表达，如下式：
[0069][0070]
式中：θn为偏离量在极坐标系中角度编码；rn为偏离量在极坐标系中半 径编码；
[0071]
3)如图2所示，格拉姆角差场矩阵生成：在得到偏离量时间序列的极 坐标系形式后，利用各时间断面偏离量角度差的正弦运算生成格拉姆角差场 矩阵，以识别多时间断面间时序自相关性，如下式：
[0072][0073]
式中：g为p～t时段偏离量dp～t对应格拉姆角差场矩阵；
[0074]
4)如图2所示，将格拉姆角差场变换所得格拉姆角差场矩阵g进行图 像编码，得到p～t时段偏离量dp～t对应的特征图像。
[0075]
s242，确定未来时段偏离量对应的分类标签；
[0076]
采取定量与定性相结合的方式，依据未来系统偏离量变化将系统未来运行 趋势划分为健康、恢复、临界、紧急四种状态：
[0077]
1)判断max(d
t～f
)≤dh是否成立，其中，dh为健康状态偏离量阈值，若是， 则未来运行趋势为健康状态；若否，转到2)；
[0078]
2)判断dh＜max(d
t～f
)≤dc是否成立，其中，dc为临界状态偏离量阈值，若 是，转到3)；若否，转到4)；
[0079]
3)判断是否成立，若是，则未来运行趋势为临界状态；若 否，则未来运行趋势为恢复状态；
[0080]
4)判断是否成立，若是，则未来运行趋势为紧急状态；若 否，则未来运行趋势为恢复状态；
[0081]
经分类所得的未来运行趋势即为t～f时段偏离量dt～f对应的分类标签。
[0082]
步骤s3、构建并训练模型；
[0083]
s31，构建卷积神经网络模型；
[0084]
如图3所示，构建包括1层输入层、2层卷积层、2层池化层、2层激活层、 1层扁平层、1层全连接层与1层输出层的卷积神经网络模型，该卷积神经网络 模型的卷积公式如下：
[0085][0086]
式中：x为卷积计算输入矩阵；w代表大小为m
×
n的卷积核；y为卷积计算 输出矩阵。
[0087]
s32，将步骤s241中得到的特征图像和步骤s242中得到的分类标签作为样 本集对卷积神经网络模型进行训练，具体而言，将步骤s241中得到的特征图像 作为卷积神经网络模型的输入，将步骤s242中得到的分类标签作为卷积神经网 络模型的输出，训练得到卷积神经网络模型输入-输出的映射关系。
[0088]
步骤s4，预测当前时段(下文以a-b时段表示)的未来运行趋势；
[0089]
s41，参照步骤s1，采集电-气综合能源系统a-b时段的量测数据。
[0090]
s42，参照步骤s21对步骤s41中得到的量测数据噪声滤波处理。
[0091]
s43，参照步骤s22，计算电-气综合能源系统a-b时段偏离量。
[0092]
s44，参照步骤s24，得到a-b时段偏离量对应的特征图像。
[0093]
s45，将步骤s44中得到的图像特征输入训练好的卷积神经网络模型中，卷 积神经网络模型先利用其输入层输入该某时段格拉姆角差场图像并将其转换为 32
×
32
×
3的像素矩阵；再通过两次“卷积-激活-池化”处理得到5
×
5大小的128 个特征图，其中激活层均采用relu作为激活函数，池化层均采用最大池化方式； 接着利用扁平层将三维特征图压缩得到长度为3200一维数组；最后通过全连接 层将扁平层信息传递至输出层，利用softmax分类器得出系统未来运行趋势的分 类标签。
[0094]
作为本发明的另一面，一种计及时序相关性的数据驱动电-气综合能源系统 态势感知模型，包括前述计及时序相关性的数据驱动电-气综合能源系统态势感 知方法中的卷积神经网络模型。
[0095]
为测试本发明所涉计及时序相关性的数据驱动电-气综合能源系统态势感知 方法的性能，下面以14节点电力系统与7节点天然气系统耦合的电-气综合能源 系统为例进行说明。在本实施例中，电-气综合能源系统的参数如表1、表2所示， 卷积神经网络各层级参数与操作如表3所示。
[0096]
表1电力系统参数表
[0097][0098]
表2天然气系统参数表
[0099][0100]
表3卷积神经网络各层级参数、操作与输入输出
[0101][0102]
以上述电-气综合能源系统为例，采用本发明所涉计及时序相关性的数据驱 动电-气综合能源系统态势感知方法开展系统运行趋势态势感知，在该过程中， 量测数据一共采集2000组90分钟尺度数据，采集时间尺度设置为5分钟，p～t 间隔设置为6，t～f间隔设置为3。卷积神经网络模型训练时样本集中的训练集与 验证集比值设置为7：3。最终本实施例得到的态势感知方法的性能指标如表4 所示。
[0103]
表4态势感知性能指标
[0104][0105]
另外，在本实施例中，图1与图2展示了单次态势感知对应的状态量数据， 此次预测的未来运行趋势为健康状态，与实际相符。
[0106]
由表4以及图4、5展示的结果看，本发明所涉计及时序相关性的数据驱动 电-气综合能源系统态势感知方法在态势感知方面的准确率很高，这可有力证明 本发明的可行性与准确性。
[0107]
上述实施例为本发明较佳的实现方案，除此之外，本发明还可以其它方式现 实，在不脱离本技术方案构思的前提下任何显而易见的替换均在本发明的保护范 围之内。
[0108]
为了让本领域普通技术人员更方便地理解本发明相对于现有技术的改进之 处，本发明的一些附图和描述已经被简化，并且为了清楚起见，本技术文件还省 略了一些其他元素，本领域普通技术人员应该意识到这些省略的元素也可构成本 发明的内容。