一种电力系统高维静态安全域边界拟合方法与流程

技术特征：
1.一种电力系统高维静态安全域边界拟合方法，其特征在于，包括以下步骤：步骤1、初始样本集的生成；步骤2、为初始样本集中每一个样本x
i
寻找一个处于安全域内的样本x
s,i
，一个安全域外样本x
u,i
，并满足||x
u,i-x
s,i
||2≤ξ，ξ是距离阈值，工程需求设定；形成边界样本集s
b
；步骤3、将边界样本集s
b
中数据所处的原始功率注入空间x，通过深度神经网络模型转换到新的三维特征空间z中；步骤4、通过基于信息增益比的加权倾斜决策树算法在特征空间z中进行边界的提取，并评估边界性能，选出最优边界。2.根据权利要求1所述电力系统高维静态安全域边界拟合方法，其特征在于，步骤1的实现包括以下步骤：步骤1.1、基于拉丁超立方抽样算法生成发电机出力数据；设一个n节点的电力系统中，含n
g
个发电机节点，n
l
个负荷节点，发电机出力上下限为：[p
g_min,j
,p
g_max,j
],j∈1,2,...,n
g
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)式中，p
g_min,j
为第j个发电机的出力下限，p
g_max,j
为第j个发电机的出力上限；设置初始样本集采样数量为n
init
，则基于拉丁超立方抽样发电机出力数据为：式中，p
g,i,j
为第i个样本中第j个发电机的出力，且：步骤1.2、发电机出力约束的负荷节点采样方法包括以下步骤：步骤1.2.1、针对第i个样本，发电机组总出力采样值为p
g_out,i
，则系统平衡机出力上下限为[p
s_min
,p
s_max
]，设n
l
个负荷节点总需求采样值为不计网损情况下，为满足系统平衡机出力不越限，p
l,i
需满足：式(4)中的约束在n
l
维空间中围成了一个凸超多面体；将问题建模为在式(5)中的超多面体中对n
l
维向量进行采样：s.t.a
×
p
l
≤b
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)式中，a∈(n
l
2)
×
n
l
，p
l
∈n
l
×
1，b∈(n
l
2)
×
1，均为式(4)中的约束系数；步骤1.2.2、构造线性目标函数：式中，f(p
l
)是变量为p
l
的线性函数；
依据线性规划方法求得最优解p
l,0
，p
l,0
为式(5)所表示的超多面体的一个顶点；步骤1.2.3、以p
l,0
为初始值，求解超多面体的中心；构建一个非线性优化问题：式中，a
i,:
为矩阵a中第i行元素，p
l,cen
为上述凸超多面体几何中心；步骤1.2.4、以p
l,cen
为出发点，随机初始化采样方向u∈r
nl
×1，求解经过点p
l,cen
并以u为方向的直线与多面体的两个交点o1,o2：式中，z，c为中间计算量，c
p
、c
n
为该直线与多面体各个平面、包含延伸面的交点，o1、o2为该直线与多面体的实际交点；步1.2.5、在点o1,o2间进行均匀采样，随机生成采样系数λ∈[0,1]，得到采样点：p
l,sam
＝p
l,cen
λ
·
(ο1 (o
2-o1))
×
u
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)式中，p
l,sam
为最后采样点；步骤1.3、依据步骤1.2对n
init
个样本中负荷节点值进行采样，并通过交流潮流计算分析，筛选出处于安全域内的样本点，形成初始样本集合s＝{(x
i
,y
i
)},i∈1,2,...,n
s
，y
i
∈{0}，其中，s为初始样本集合，x
i
为第i个样本中节点注入功率向量，由发电机出力向量p
g,i
与节点负荷向量p
l,i
组成，n
s
为集合s中样本数，y
i
∈{0}为样本标签，0的含义为该样本处于安全域内。3.根据权利要求2所述电力系统高维静态安全域边界拟合方法，其特征在于，步骤2的实现包括以下步骤：步骤2.1、设定边界距离阈值ξ，搜索初始步长λ
r,set
；步骤2.2：依据样本集合s的样本数n
s
，基于拉丁超立方抽样算法采样n
s
个n
p
维功率增长方向向量，并对每一个方向向量进行归一化，形成搜索方向向量矩阵步骤2.3：对每一个处于样本集合s中的x
i
，执行如下算法得到边界样本集合s
b
，其中x
s,i
为样本对中处于安全域的边界样本，x
u,i
为样本对中处于安全域外的边界样本，集合s
b
中共2n
s
个样本；
4.根据权利要求3所述电力系统高维静态安全域边界拟合方法，其特征在于，步骤3的实现包括以下步骤：步骤3.1、依据边界样本集数据量2n
s
及数据维度n
p
，选择神经网络进行特征非线性变换；若样本集数量小于10000及数据维度小于100，选择全连接神经网络；否则选择卷积神经网络或循环神经网络；步骤3.2、选择神经网络优化器为自适应矩估计算法，设定神经网络超参数，进行模型训练，并选出最优模型，将边界样本集s
b
中数据所处的原始功率注入空间x，通过深度神经网络模型转换到新的三维特征空间z中，得到特征空间下的数据集s
z
。5.根据权利要求4所述电力系统高维静态安全域边界拟合方法，其特征在于，步骤4的实现包括以下步骤：步骤4.1、设定信息增益比的加权倾斜决策树模型中每一非叶子节点下，其输入条件为训练数据集s
z
，(z
i
,y
i
)∈s
z
,i＝1,2,3,...,2n
s
；设定模型最大深度为d，当前节点深度为d；其中，根节点下训练数据集为步骤3.2生成的数据集s
z
，其他非叶子结点下的训练数据集为其父节点的训练集划分得到的左子集s'
l
、右子集s'
r
；步骤4.2、基于数据集s
z
创建模型根节点g，并令当前节点深度d＝0；步骤4.3、如果当前节点深度d>模型最大深度d，则节点g设为叶子结点，其标签为数据集s
z
内样本数量最多的对应标签k；否则，转步骤4.4；步骤4.4、如果数据集s
z
内所有样本属于同一类别k，则节点g设为叶子结点，其标签为k；否则，转步骤4.5；步骤4.5、随机初始化模型参数θ0；步骤4.6、基于拟牛顿算法和初始值θ0求解模型目标函数最小值式(10)，得到模型最优
参数θ
best
；其中，l(θ)为模型目标函数，λ为l2正则化项系数，θ为模型各节点下的待训练参数，||θ||2为θ的二范数，h(s
z
)为样本集合s的经验熵，h(s
z
|θ)为样本集合s
z
在θ下的条件经验熵；其中，k是样本总类别数；k表示样本中第k类样本标签；|s
k
|是样本集合s
z
中第k类样本数量，|s
z
|是样本集合s
z
总样本数量；其中，w
l
为所有样本分属左子集的权重之和，w
r
为所有样本分属右子集的权重之和，h
l
为左子节点加权信息熵，h
r
为右子节点加权信息熵，m是该节点下样本集合s
z
总样本数量，θ为模型各节点下的待训练参数；其中，为样本(x
i
,y
i
)属于左子节点的权重，s
l
各样本关联属于左子节点权重信息的集合；其中，为样本(x
i
,y
i
)属于右子节点的权重，s
r
各样本关联属于右子节点权重信息的集合；其中，k是样本总类别数；k表示样本中第k类样本标签，为样本集合中k类别下样本属于左子集的权重之和，w
l
为样本集合s
z
中所有样本属于左子集的权重之和；其中，k是样本总类别数；k表示样本中第k类样本标签，w
rk
为样本集合中k类别下样本属于右子集的权重之和；w
r
为样本集合s
z
中所有样本属于右子集的权重之和；其中，(x
i
,y
i
)表示样本集合s
z
中第i个样本，s
l
为各样本关联分属左子集权重信息的集合，s
r
为各样本关联分属右子集权重信息的集合；
其中，为第i个样本分属左子节点的权重，为第i个样本分属右子节点的权重，σ(
·
)为sigmoid函数；其中，k是样本总类别数；k表示样本中第k类样本标签，为样本集合中k类别下样本属于左子集的权重之和；其中，k是样本总类别数；k表示样本中第k类样本标签，为样本集合中k类别下样本属于右子集的权重之和；步骤4.7、基于最优参数θ
best
得到左子集、右子集，s'
l
,s'
r
；其中，其中，(x
i
,y
i
)表示样本集合s
z
中第i个样本，为第i个样本分属左子节点的权重，s
l
为各样本关联分属左子集权重信息的集合，为第i个样本分属右子节点的权重，s
r
为各样本关联分属右子集权重信息的集合；步骤4.8、如果s'
l
样本数量小于阈值n
y
，置左子节点为叶节点，标签为s'
l
中数量最多的类别，否则，构建左子节点g
l
，并令该节点下训练集为s'
l
，d＝d 1，转步骤4.3；步骤4.9、如果s'
r
样本数量小于阈值n
y
，置右子节点为叶节点，标签为s'
r
中数量最多的类别，否则，构建右子节点g
r
，并令该节点下训练集为s'
r
，d＝d 1，转步骤4.3；步骤4.10、输出决策树模型，并转为边界表达式。

技术总结
本发明涉及电力行业安全域边界拟合技术，具体涉及一种电力系统高维静态安全域边界拟合方法，考虑系统发电机出力约束，进行节点负荷值采样，通过潮流计算，筛选出处于安全域内的样本，形成初始样本集。提出边界样本搜索算法，为初始样本集中每一个样本寻找一个处于安全域内的样本，一个安全域外样本，并满足二者距离小于设定的距离阈值。将所处原始功率注入空间的边界样本集中数据通过深度神经网络模型转换到新的三维特征空间中。通过基于信息增益比的加权倾斜决策树算法在特征空间中进行边界的提取，并评估边界性能，选出最优边界。该方法可实现电力系统安全域高维边界的拟合，并缩小与实际边界的误差。缩小与实际边界的误差。缩小与实际边界的误差。


技术研发人员：张俊 高天露 戴宇欣 张科 许沛东 陈思远 周自强 陈晓刚 章姝俊 李乃一
受保护的技术使用者：国网浙江省电力有限公司
技术研发日：2022.10.28
技术公布日：2023/2/3