基于数据驱动升维灵敏度分析的电力系统电压控制方法技术方案

本发明专利技术公开了基于数据驱动升维灵敏度分析的电力系统电压控制方法，包括如下步骤：S1、初始化模型参数；S2、以节点电压和无功历史量测数据作为训练样本输入构建线性潮流方程矩阵；S3、基于所述线性潮流方程矩阵，计算节点电压与无功灵敏度系数；S4、构建电力系统电压优化控制模型；S5、采用基于矩阵分裂的分布式算法求解步骤S4中得到的含综合能源路由电力系统电压优化控制模型，直至满足收敛精度；S6、输出步骤S5的求解结果，包括目标函数值、各节点电压幅值、综合能源路由无功注入调整量。方案设计电力系统电压优化控制模型，收敛速度快、精度高，具有更强的工程适用性。具有更强的工程适用性。具有更强的工程适用性。

【技术实现步骤摘要】
基于数据驱动升维灵敏度分析的电力系统电压控制方法


[0001]本专利技术涉及电力系统运行控制
，具体的，涉及基于数据驱动升维灵敏度分析 的电力系统电压控制方法。

技术介绍

[0002]随着分布式新能源的大量开发利用，分布式能源以集群化形式接入电力系统，给电力 系统安全稳定运行、供电可靠性带来挑战。分布式能源的不确定性加剧了电力系统运行状态 的波动，大量分布式能源进行协调增加了电力系统电压控制的复杂性。
[0003]与此同时，多类型能源的接入提升了电力系统运行的灵活性，也增加了系统运行的复 杂度。具体而言，综合能源以综合能源路由的形式接入，能够灵活控制电网侧功率，参与电 网调节。另外，综合能源路由的功率注入改变了电网传统单一的潮流方向，导致电压波动加 剧，影响用电设备寿命的同时造成脱网风险。综合能源路由自身具有快速的无功调节能力， 随着综合能源路由接入量的逐渐提升，其电压控制作用成为电力系统中无法忽视的调节手段， 由于无需配置额外的调节装置，其调节成本更低，且对出力随机波动具有更快响应速度，对 于保障系统安全高效运行具有关键意义。
[0004]集中式电压控制方法需基于集中式通信网络对全局模型进行统一维护，并实时采集大 量运行数据，不仅通信时延较长，优化模型的求解更可能遭遇“维数灾”的问题，收敛速度 较慢甚至难以收敛。分布式控制更适合含综合能源路由的电力系统电压控制，基于多代理系 统的应用。但这传统分布式控制方法仅具备一阶收敛速度，当综合能源路由数量较大时，电 压控制收敛速度将极其缓慢，无法满足在线控制的需求。近年来，基于二阶梯度的分布式方 法得到了广泛研究应用，由于采用了二阶梯度，其迭代过程具有超线性收敛速度，但其矩阵 的求逆过程一般需要全局信息，因此对分布式迭代控制提出了较高的要求。对于综合能源路 由广泛接入的电力系统而言，调度中心往往不掌握其精确模型参数，将对分布式电压控制效 果及收敛速度造成不良影响。为解决上述问题，需要发展一种具有更快的收敛速度和更优的 控制收敛结果，且不受参数不精确问题影响，具有更强的工程适用性的综合能源路由电压控 制方法。

技术实现思路

[0005]本专利技术的目的是提出了基于数据驱动升维灵敏度分析的含综合能源路由电力系统电 压控制方法，通过利用矩阵分裂方法，实现了矩阵的分布式求逆，从而将分布式控制敛速提 升至超线性收敛，更适合在线应用，基于Koopman升维方法，利用综合能源路由历史运行数 据作为训练样本，构建高维线性精确潮流模型，从而推导得到电压
‑
无功全局灵敏度，以此校 正分布式控制中的迭代方向使得模型具备更快的收敛速度和更优的控制收敛结果，且不受参 数不精确问题影响，具有更强的工程适用性。
[0006]本专利技术实施例中提供的一种技术方案是基于数据驱动升维灵敏度分析的电力系统电 压控制方法，包括如下步骤：
S1、输入综合能源路由的接入位置和容量、综合能源路由无功调节量的下限和上限值，含综 合能源路由电力系统基准电压和基准功率的初值、各节点安全电压约束、节点电压和无功历 史量测数据、训练样本数量、升维维度、节点电压自然分布向量、本地二阶梯度矩阵初值、 修正参数、迭代步长、收敛精度；S2、以节点电压和无功历史量测数据作为训练样本输入，采用基于Koopman状态空间升维变 换方法构建线性潮流方程矩阵；S3、基于所述线性潮流方程矩阵，计算节点电压与无功灵敏度系数；S4、在线性化电压
‑
无功方程的基础上，构建含综合能源路由电力系统电压优化控制模型，包 括：设定含综合能源路由电力系统电压偏差之和最小为目标函数，考虑综合能源路由无功调 节能力受逆变器容量的约束；S5、采用基于矩阵分裂的分布式算法求解步骤S4中得到的含综合能源路由电力系统电压优化 控制模型，直至满足收敛精度；S6、输出步骤S5的求解结果。
[0007]作为优选，S2，包括如下步骤：S21、对输入变量进行升维变换得到升维变换扩充出的第t维输入变量；S22、利用最小二乘方法估计线性潮流方程矩阵L的值。
[0008]作为优选，S21包括：对输入变量进行升维变换，升维后的输入变量公式表示如下：式中，x为输入变量，包括节点注入有功和无功功率，ψ(x)表示升维扩充出的m维输入变量， x
lift
为升维后的输入变量；升维变换扩充出的第t维输入变量，公式表示如下：式中，ψ
t
(x)为ψ(x)中的第t个元素，f
lift
为升维函数，x
v
为x中的第v个元素，K为x的维 度，c
t
为随机生成的K维基底向量，c
tv
为c
t
中的第v个元素。
[0009]作为优选，S22包括：利用最小二乘方法估计线性统一潮流方程矩阵L的值，公式如下：式中，为输入变量样本矩阵，包含S个时间断面的节点相角及电压幅值向量，为升维 后的输入变量样本矩阵，包含S个时间断面的升维后的输入变量，S为训练样本数量，为 状态变量样本矩阵，包含S个时间断面的节点相角及电压幅值向量，[
·
]+
为矩阵的 Moore
‑
Penrose逆。
[0010]作为优选，S3中，计算节点电压与无功灵敏度系数，公式如下：
式中，Z
ji
为节点电压与无功灵敏度矩阵Z中第j行第i列的元素，L
ji
为矩阵L中对应电压V
j
和无功调整量
△
q
i
的元素，L
j,(K+t)
为矩阵L中对应电压V
j
和升维变换扩充出的第t维输入变量 ψ
t
(x)的值，m为升维维度。
[0011]作为优选，S4中，所述的线性化电压
‑
无功方程表示为：其中，
△
q为节点综合能源路由无功注入调整量，为综合能源路由无功不进行调节时的节点 电压自然分布向量。
[0012]作为优选，S4中，电力系统电压优化控制模型，公式如下：作为优选，S4中，电力系统电压优化控制模型，公式如下：式中，f(Δq)为目标函数，代表矩阵元素的平方和，v为节点额定电压所构成的列向量，为额定电压与自然电压分布之间的偏差，Δq和为节点综合能源路由无功调节量 的下限和上限值。
[0013]作为优选，S5中，S51、计算各个节点自身的目标函数f
i
(
△
q(k))：式中，f
i
(
△
q(k))为节点i自身的目标函数，
△
q(k)为第k步迭代时综合能源路由无功注入调 整量构成的向量，N
i
为与节点i通过一条支路直接相连，且相对于参考节点位于下游的节点 所构成的集合，Z
j
为Z中的第j行，为中第j个元素；S52、计算节点i自身目标函数的梯度：其中，g
i
(k)为g(k)中的第i个元素，V...

【技术保护点】

【技术特征摘要】
1.基于数据驱动升维灵敏度分析的电力系统电压控制方法，其特征在于，包括如下步骤：S1、输入综合能源路由的接入位置和容量、综合能源路由无功调节量的下限和上限值，含综合能源路由电力系统基准电压和基准功率的初值、各节点安全电压约束、节点电压和无功历史量测数据、训练样本数量、升维维度、节点电压自然分布向量、本地二阶梯度矩阵初值、修正参数、迭代步长、收敛精度；S2、以节点电压和无功历史量测数据作为训练样本输入，采用基于Koopman状态空间升维变换方法构建线性潮流方程矩阵；S3、基于所述线性潮流方程矩阵，计算节点电压与无功灵敏度系数；S4、在线性化电压
‑
无功方程的基础上，构建含综合能源路由电力系统电压优化控制模型；S5、采用基于矩阵分裂的分布式算法求解步骤S4中得到的含综合能源路由电力系统电压优化控制模型，直至满足收敛精度；S6、输出步骤S5的求解结果。2.根据权利要求1所述的基于数据驱动升维灵敏度分析的电力系统电压控制方法，其特征在于，S2，包括如下步骤：S21、对输入变量进行升维变换得到升维变换扩充出的第t维输入变量；S22、利用最小二乘方法估计线性潮流方程矩阵L的值。3.根据权利要求2所述的基于数据驱动升维灵敏度分析的电力系统电压控制方法，其特征在于，S21包括：对输入变量进行升维变换，升维后的输入变量公式表示如下：式中，x为输入变量，包括节点注入有功和无功功率，ψ(x)表示升维扩充出的m维输入变量，x
lift
为升维后的输入变量；升维变换扩充出的第t维输入变量，公式表示如下：式中，ψ
t
(x)为ψ(x)中的第t个元素，f
lift
为升维函数，x
v
为x中的第v个元素，K为x的维度，c
t
为随机生成的K维基底向量，c
tv
为c
t
中的第v个元素。4.根据权利要求2所述的基于数据驱动升维灵敏度分析的电力系统电压控制方法，其特征在于，S22包括：利用最小二乘方法估计线性统一潮流方程矩阵L的值，公式如下：
式中，为输入变量样本矩阵，包含S个时间断面的节点相角及电压幅值向量，为升维后的输入变量样本矩阵，包含S个时间断面的升维后的输入变量，S为训练样本数量，为状态变量样本矩阵，包含S个时间断面的节点相角及电压幅值向量，[
·
]
+
为矩阵的Moore
‑
Penrose逆。5.根据权利要求1所述的基于数据驱动升维灵敏度分析的电力系统电压控制方法，其特征在于，S3中，计算节点电压与无功灵敏度系数，公式如下：式中，Z
ji
为节点电压与无功灵敏度矩阵Z中第j行第i列的元素，L
ji
为矩阵L中对应电压V
j
和无功调整量Δq
i
的元素，L
j,(K+t)
为矩阵L中对应电压V
...

【专利技术属性】
技术研发人员：黄天恩，莫雅俊，吴振杰，王源涛，廖培，谢佳烨，唐剑，王中冠，李祥，徐双蝶，王艳，李城达，张超，夏衍，董航，周依希，孙思聪，张洁，许鹏，陈嘉宁，苏熀兴，杨兴超，李跃华，黄佳斌，来益博，胡浩，
申请(专利权)人：天津大学，
类型：发明
国别省市：