技术特征：
1.一种谐波发射水平评估方法，其特征在于，包括如下步骤：测量用户侧停机时系统侧谐波电压v
s
，使用扩散核密度估计算法计算系统侧谐波电压的概率密度函数f(v
s
)，将系统谐波阻抗z
s
的实部r和虚部x视为随机变量，利用jeffrey先验算法估计r和x的先验概率p(r,x)，结合r和x的先验概率p(r,x)和系统侧谐波电压的概率密度函数f(v
s
)计算条件概率p(v
sx
|r,x)；根据先验概率p(r,x)、概率密度函数f(v
s
)和条件概率p(v
sx
|r,x)计算得到r和x的后验概率p(r,x|v
sx
)，根据r和x的后验概率p(r,x|v
sx
)计算得到系统谐波阻抗z
s
的估计值；根据系统谐波阻抗z
s
的估计值计算用户侧和系统侧的谐波发射水平。2.如权利要求1所述的谐波发射水平评估方法，其特征在于：使用扩散核密度估计算法，计算系统侧谐波电压的概率密度函数f(v
s
)，具体包括：以如下核密度估计表达式计算v
s
的概率密度函数：式中，v
sxi
为测得的系统侧谐波电压实部的样本值，n为样本容量，h为带宽，k为核函数。3.如权利要求2所述的谐波发射水平评估方法，其特征在于：带宽h的最优带宽为h*，采用如下公式计算：式中，k
d
(
·
)为扩散核函数，n代表高斯核密度估计的初始条件系数，x和t为k
d
(
·
)中的随机变量，a(x)、b(x)为在定义域中具有二阶导数的任何正函数。4.如权利要求1所述的谐波发射水平评估方法，其特征在于：利用jeffrey先验算法估计r和x的先验概率p(r,x)，具体包括：将z
s
的先验分布视为正态分布，参数为(μ,σ2)，利用jeffrey先验方法估计参数(μ,σ2)，(μ,σ2)的分布如下：根据贝叶斯公式，(μ,σ2)的联合后验分布由如下公式获得式中，和分别是μ和σ2的数学期望，根据标准化熵损失函数的定义，参数μ和σ2的估计取数学期望值，得到μ和σ2的估计值如下：
式中，e1代表后验概率的数学期望，和分别为μ和σ2的估计值，z
s
的先验分布为对先验分布函数求导得到r和x的先验概率p(r
i
,x
j
)。5.如权利要求1所述的谐波发射水平评估方法，其特征在于：结合r和x的先验概率p(r,x)和系统侧谐波电压的概率密度函数f(v
s
)计算条件概率p(v
sx
|r,x)，具体包括：将随机变量r，x的先验分布记为：通过下式计算条件概率：式中，n为用于计算条件概率样本数据的长度，p(v
sx
(t|r,x))代表v
sx
(r,x)中第t个元素的概率值。6.如权利要求1所述的谐波发射水平评估方法，其特征在于：根据先验概率p(r,x)、概率密度函数f(v
s
)和条件概率p(v
sx
|r,x)，使用贝叶斯公式计算r和x的后验概率p(r,x|v
sx
)：计算r,x的后验概率值，f(v
sx
)视为常数，上式写为：p(r,x|v
sx
)＝λp(v
sx
|r,x)p(r,x)式中，λ表示任一非零实数。7.如权利要求5或6所述的谐波发射水平评估方法，其特征在于：根据r和x的后验概率p(r,x|v
sx
)计算得到系统谐波阻抗z
s
的估计值具体包括：根据r和x的后验概率p(r,x|v
sx
)得到r和x的联合概率密度函数π(r,x|v
sx
)，根据r和x的联合概率密度函数得到边缘概率密度函数π(r|v
sx
)和π(x|v
sx
)：式中，π(r,x|v
sx
)为r和x的联合概率密度函数，π(r|v
sx
)和π(x|v
sx
)分别为r和x的边缘概率密度函数；通过计算边缘概率密度函数的数学期望得到系统谐波阻抗实部和虚部的估计值：式中，z
sx
和z
sy
分别为系统谐波阻抗实部和虚部的估计值，根据z
sx
和z
sy
计算系统谐波阻抗z
s
的估计值，即z
s
＝z
sx
+jz
sy
。8.如权利要求1所述的谐波发射水平评估方法，其特征在于：根据系统谐波阻抗z
s
的估计值计算用户侧和系统侧的谐波发射水平，具体包括：根据如下公式计算用户侧和系统侧谐波发射水平：v
pcc_c
＝z
s
i
pcc
v
pcc_u
＝v
pcc-z
s
i
pcc
其中v
pcc_c
为用户侧的谐波发射水平，v
pcc_u
为系统侧的谐波发射水平，i
pcc
为共同耦合点处的谐波电流，v
pcc
为共同耦合点处的谐波电压。9.一种谐波发射水平评估装置，其特征在于：包括条件概率计算模块，用于测量用户侧停机时系统侧谐波电压v
s
，使用扩散核密度估计算法计算系统侧谐波电压的概率密度函数f(v
s
)，将系统谐波阻抗z
s
的实部r和虚部x视为随机变量，利用jeffrey先验算法估计r和x的先验概率p(r,x)，结合r和x的先验概率p(r,x)和系统侧谐波电压的概率密度函数f(v
s
)计算条件概率p(v
sx
|r,x)；系统谐波阻抗估计值计算模块，用于根据先验概率p(r,x)、概率密度函数f(v
s
)和条件概率p(v
sx
|r,x)计算得到r和x的后验概率p(r,x|v
sx
)，根据r和x的后验概率p(r,x|v
sx
)计算得到系统谐波阻抗z
s
的估计值；系统谐波阻抗计算模块，用于根据系统谐波阻抗z
s
的估计值计算用户侧和系统侧的谐波发射水平。10.如权利要求9所述的谐波发射水平评估装置，其特征在于：所述系统谐波阻抗计算模块根据系统谐波阻抗z
s
的估计值计算用户侧和系统侧的谐波发射水平，具体包括：根据如下公式计算用户侧和系统侧谐波发射水平：v
pcc_c
＝z
s
i
pccvpcc_u
＝v
pcc-z
s
i
pcc
其中v
pcc_c
为用户侧的谐波发射水平，v
pcc_u
为系统侧的谐波发射水平，i
pcc
为共同耦合点处的谐波电流，v
pcc
为共同耦合点处的谐波电压。11.一种电子设备，其特征在于：包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现如权利要求1-8中任一项所述的谐波发射水平评估方法。12.一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如权利要求1-8中任一项所述的谐波发射水平评估方法。
技术总结
本发明提供一种谐波发射水平评估方法及装置，所述方法包括：测量用户侧停机时系统侧谐波电压V
技术研发人员：柳丹 冀肖彤 李小平 王伟 梅欣 邓万婷 陈孝明 曹侃 叶畅 康逸群 江克证 熊平 胡畔 肖繁 谭道军
受保护的技术使用者：国网湖北省电力有限公司
技术研发日：2022.04.07
技术公布日：2022/8/5