首页 / 专利库 / 信号处理 / 信号 / 一种配电网的自适应同步相量量测方法

一种配电网的自适应同步相量量测方法

阅读:738发布:2021-06-11

专利汇可以提供一种配电网的自适应同步相量量测方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种配 电网 的自适应同步相量量测方法,属于电 力 系统配电网同步相量测量技术领域。首先,设计一个4周波的稳态有限冲激响应 滤波器 和一个2周 波动 态有限冲激响应滤波器,这两个滤波器均考虑了基频附近难以滤除的带外干扰成分,且中心 频率 具有在线调整的能力。其次,基于现代谱估计方法,利用10周波 信号 采样 值分析信号的频率成分;进而,根据设计的暂态识别方法,实时识别信号由稳态到动态的瞬时变化过程,并在稳态同步相量计算结果和动态同步相量计算结果中进行选择输出。本发明所设计的同步相量量测方法,能够解决配电网同步相量计算过程中带外干扰成分导致的 精度 不高的问题,同时提高了动态变化场景下同步相量计算结果的响应速度。,下面是一种配电网的自适应同步相量量测方法专利的具体信息内容。

1.一种配电网的自适应同步相量量测方法,其特征在于该方法包括以下步骤:
(1)采集配电网的电压电流信号,包括以下步骤:
(1-1)以时间间隔T采集配电网的电压或电流信号,将当前采样数据点的前10个周波采样数据记为X;
(1-2)取采样数据X中当前采样数据点前的4个周波数据,记为Xs,采样数据的个数记为N1,将采样数据分别编号-K1,-(K1-1),……,0,1,……,(K1-1),K1;
(1-3)取采样数据X中当前采样数据点前的2个周波数据,记为Xd,采样数据个数记为N2,将采样数据分别编号-K2,-(K2-1),……,0,1,……,(K2-1),K2;
(2)采用现代谱估计方法,确定步骤(1)的采样数据X中电压或电流信号的频率在25Hz~150Hz范围内的谐波和/或间谐波的数目和频率值;
(3)根据步骤(1-2)周波数据Xs和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个稳态有限冲激响应滤波器,该滤波器的系数为Coeff_SteadyFilter:
Coeff_SteadyFilter=Fs(1,:)
其中,Fs为稳态滤波器系数矩阵,将稳态滤波器系数矩阵Fs中的第一行作为稳态有限冲激响应滤波器的系数Coeff_SteadyFilter,稳态滤波器系数矩阵Fs由稳态信号系数矩阵Cs和稳态滤波器权重矩阵W1计算得到,其中:
Fs=(CHsW1HW1Cs)-1CsW1
稳态信号系数矩阵Cs的计算公式为:
其中,上标H为矩阵共轭转置,e为自然对数的底数,j为虚数单位,K1为步骤(1-2)中电压或电流采样数据的编号值的最大值,T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准频率值,ωc等于配电网基准频率与2π的乘积,ωbi为电压或电流信号中频率在
25Hz~150Hz范围内的谐波和/或间谐波的角频率值,ωbi等于电压或电流信号中频率在
25Hz~150Hz范围内的谐波和/或间谐波的频率值与2π的乘积;
稳态滤波器权重矩阵W1为对角矩阵,即除了矩阵对角线上元素不为零外,矩阵其余元素均为零,计算公式为:
其中,w1,w2,……,wN1为稳态滤波器权重矩阵W1的对角元素,计算公式为:
n取1,2,…,N1
其中,N1为所设计的稳态有限冲激响应滤波器的阶数,该阶数等于4个周波电压或电流采样数据的点数,N1=2*K1+1,K1即为4周波电压或电流采样数据的编号值的最大值,cosh表示双曲余弦函数,ac和ρm分别是窗函数的特性调节系数,ac的取值为7.0,ρm的取值为1.2;
(4)根据步骤(1-3)周波数据Xd和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个动态有限冲激响应滤波器,该滤波器的系数为Coeff_DynamicFilter,计算公式为:
Coeff_DynamicFilter=Fd(1:3,:)
其中,Fd为动态滤波器系数矩阵,将动态滤波器系数矩阵Fd中的第一行至第三行作为动态有限冲激响应滤波器的系数CoeffDynamicFilter,动态滤波器系数矩阵Fd由动态信号系数矩阵Cd和权重矩阵W2计算得到,计算公式为:
Fd=(CHdW2HW2Cd)-1CdW2
其中,动态滤波器权重矩阵W2的计算公式如下:
其中,w1,w2,……,wN2为动态滤波器权重矩阵W2的对角元素,计算公式为:
n取1,2,…,N2
其中,N2为动态有限冲激响应滤波器的阶数,等于2个周波电压或电流采样值的点数,N2=2*K2+1,K2即为步骤(1-3)中电压或电流采样数据的编号值的最大值;
动态信号系数矩阵Cd的计算公式为:
其中,动态信号系数矩阵Cd的元素有两类,一类元素与信号的基准频率有关,其计算公式如下:
k取值-K2,-(K2-1),...,0,...,K2
其中,c(k,c)为基准频率信号的泰勒级数展开系数行向量,由电压或电流基准频率信号在编号值为k的采样数据点处的二阶泰勒展开计算得到, 为行向量c(k,c)的共轭行向量,其中每一个元素都是向量c(k,c)的共轭值,基准频率e为自然对数的底数,j为虚数单位,K2为步骤(1-3)中2周波电压或电流采样数据的编号值的最大值,K2=(N2-1)/2,N2为2周波电压或电流采样数据的点数,T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准角频率值,等于配电网基准频率与2π的乘积,
动态信号系数矩阵Cd的第二类元素,与电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波有关,其计算公式如下:
k取值-K2,-(K2-1),...,0,...,K2
其中,c(k,bi)为谐波和/或间谐波信号的泰勒级数展开系数行向量,由电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波信号在编号为k的采样数据点处的二阶泰勒展开计算得到, 为c(k,bi)的共轭行向量,其中每一个元素都是向量c(k,bi)的共轭值,(5)根据步骤(3)中设计的稳态有限冲激响应滤波器Coeff_SteadyFilter和步骤(1)的4周波采样数据Xs,计算配电网的稳态同步相量,计算公式如下:
ps=Coeff_SteadyFilter·Xs
其中,ps为计算得到的电压或电流信号的稳态相量值;
(6)根据步骤(4)中设计的动态有限冲激响应滤波器Coeff_DynamicFilter和步骤(1)的2周波采样数据Xd,计算配电网的动态同步相量pd,计算公式如下:
pd=Coeff_DynamicFilter·Xd
其中, 其中, 为动态相量值, 为动态相量一阶导数值,
为动态相量二阶导数值;
根据上述动态相量值 计算出配电网的电压或电流信号动态相量的幅值A1和相位计算公式如下:
其中, 表示动态相量的模值,即动态相量的幅值,a0为动态相量幅值的零阶导数,表示动态相量的相位,b0为动态相量相位的零阶导数值;
根据上述动态相量值 及一阶导数值 计算得到配电网电压或电流信号的频率值
fc,计算公式如下:
其中,f0为配电网电压或电流信号的额定频率值,b1为相位的一阶导数值:
其中,Im{}表示求取复数的虚部值,j为虚数单位;
根据上述动态相量值 一阶导数值 和二阶导数值 计算得到配电网电压或
电流信号的频率变化率值ROCOF,计算公式如下:
其中,b2为动态相量值 的相位的二阶导数值,计算公式如下:
其中,b1为动态相量的相位的一阶导数值,a1为动态相量的幅值的一阶导数,计算公式为:
其中,Re{}表示求取复数的实部值;
(7)设定配电网电压或电流信号状态判断规则:
(7-1)根据步骤(6)中计算得到的动态相量的幅值一阶导数a1和相位一阶导数b1,设定两个状态判断因子Index1和Index2,计算公式如下:
Index1=|a1|,Index2=|b1|
(7-2)根据步骤(5)中计算得到的电压或电流信号的稳态相量ps,设定两个状态判断因子Index3和Index4,计算公式如下:
其中,K2为电压或电流信号的采样数据Xd中编号的最大值,x(k)为步骤(1)中采样数据Xd中编号为k的采样值, 为采样数据Xd中编号为k的采样值相对应的估计值,计算公式如下:
其中,|ps|为稳态相量值ps的模值,∠ps为稳态相量值ps的相位值,fc为步骤(6)计算得到的配电网电压或电流信号的频率值;
(7-3)对配电网电压或电流信号当前的状态进行判断,若满足:
{Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=1
其中,σth1、σth2、σth2和σth2分别是Index1、Index2、Index3和Index4四个状态判断因子的边界值,∪为并集符号,表示“或”的关系,∩为取交集,表示“与”的关系,则判定配电网处于动态,选择步骤(6)中计算得到的动态相量pd作为最终得到的同步相量计算值P;
若满足:
{Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=0
则判定配电网处于稳态,将步骤(5)中计算得到的稳态相量ps作为最终得到的同步相量计算值P,实现配电网的自适应同步相量量测。

说明书全文

一种配电网的自适应同步相量量测方法

技术领域

[0001] 本发明涉及一种配电网的自适应同步相量量测方法,属于电系统配电网同步相量测量技术领域。

背景技术

[0002] 随着非线性负荷、电动汽车和以换流器为接口的分布式新能源的大量接入配电网,引入了大量谐波和间谐波;此外,由于配电网位于电网用户侧,受各种电力负荷动态行为的影响较大;配电网高频噪声干扰问题也十分严重,这些都为基于高频离散采样的同步相量高精度计算提出了新的挑战。
[0003] 目前,输电网同步相量计算方法根据方法对量测精度和响应性能的侧重不同,主要可以分为两类:稳态算法和动态算法。稳态算法包括:离散傅里叶变换以及其改进算法(补偿、加窗、插值、自适应采样率等)、有限和无限冲激响应滤波器方法以及最小二乘法等;动态算法包括:卡尔曼滤波、泰勒级数变换以及小波变换等。例如,徐全,陆超,刘映尚,等公开的基于泰勒级数和离散傅里叶变换的综合自适应相量算法,电力系统自动化,2016,40(19):37-43,将电压电流量测信号作为单一基波频率正弦波信号进行研究。然而,配电网量测信号中通常含有大量的谐波和间谐波等干扰成分。通过设计具有合适截断频率的低通滤波器,可以有效地滤除一部分谐波和间谐波干扰。但是,对于那些位于基波频率信号附近,且处于滤波器通带或者过渡带的干扰成分(称为带外干扰成分)衰减效果却不明显。因此,有必要在同步相量计算方法设计中,充分考虑这些谐波和间谐波干扰成分的影响。同时,所设计的同步相量计算方法要能够兼顾动态变化场景下的计算结果的响应速度。

发明内容

[0004] 本发明的目的是提出一种配电网的自适应同步相量量测方法,以解决配电网同步相量计算过程中带外干扰成分导致的精度不高的问题,同时提高在动态变化场景下的同步相量计算结果响应速度。
[0005] 本发明提出的配电网的自适应同步相量量测方法,包括以下步骤:
[0006] (1)采集配电网的电压或电流信号,包括以下步骤:
[0007] (1-1)以时间间隔T采集配电网的电压或电流信号,将当前采样数据点的前10个周波采样数据记为X;
[0008] (1-2)取采样数据X中当前采样数据点前的4个周波数据,记为Xs,采样数据的个数记为N1,将采样数据分别编号-K1,-(K1-1),……,0,1,……,(K1-1),K1;
[0009] (1-3)取采样数据X中当前采样数据点前的2个周波数据,记为Xd,采样数据个数记为N2,将采样数据分别编号-K2,-(K2-1),……,0,1,……,(K2-1),K2;
[0010] (2)采用现代谱估计方法,确定步骤(1)的采样数据X中电压或电流信号的频率在25Hz~150Hz范围内的谐波和/或间谐波的数目和频率值;
[0011] (3)根据步骤(1-2)周波数据Xs和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个稳态有限冲激响应滤波器,该滤波器的系数为Coeff_SteadyFilter:
[0012] Coeff_SteadyFilter=Fs(1,:)
[0013] 其中,Fs为稳态滤波器系数矩阵,将稳态滤波器系数矩阵Fs中的第一行作为稳态有限冲激响应滤波器的系数Coeff_SteadyFilter,稳态滤波器系数矩阵Fs由稳态信号系数矩阵Cs和稳态滤波器权重矩阵W1计算得到,其中:
[0014] Fs=(CHsW1HW1Cs)-1CsW1
[0015] 稳态信号系数矩阵Cs的计算公式为:
[0016]
[0017] 其中,上标H为矩阵共轭转置,e为自然对数的底数,j为虚数单位,K1为步骤(1-2)中电压或电流采样数据的编号值的最大值,T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准频率值,ωc等于配电网基准频率与2π的乘积,ωbi为电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波的角频率值,ωbi等于电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波的频率值与2π的乘积;
[0018] 稳态滤波器权重矩阵W1为对角矩阵,即除了矩阵对角线上元素不为零外,矩阵其余元素均为零,计算公式为:
[0019]
[0020] 其中,w1,w2,……,wN1为稳态滤波器权重矩阵W1的对角元素,计算公式为:
[0021] n取1,2,…,N1
[0022] 其中,N1为所设计的稳态有限冲激响应滤波器的阶数,该阶数等于4个周波电压或电流采样数据的点数,N1=2*K1+1,K1即为4周波电压或电流采样数据的编号值的最大值,Cosh表示双曲余弦函数,ac和ρm分别是窗函数的特性调节系数,ac的取值为7.0,ρm的取值为1.2;
[0023] (4)根据步骤(1-3)周波数据Xd和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个动态有限冲激响应滤波器,该滤波器的系数为Coeff_dynamicFilter,计算公式为:
[0024] Coeff_DynamicFilter=Fd(1:3,:)
[0025] 其中,Fd为动态滤波器系数矩阵,将动态滤波器系数矩阵Fd中的第一行至第三行作为动态有限冲激响应滤波器的系数Coeff_dynamicFilter,动态滤波器系数矩阵Fd由动态信号系数矩阵Cd和权重矩阵W2计算得到,计算公式为:
[0026] Fd=(CHdW2HW2Cd)-1CdW2
[0027] 其中,动态滤波器权重矩阵W2的计算公式如下:
[0028]
[0029] 其中,w1,w2,……,wN2为稳态滤波器权重矩阵W2的对角元素,计算公式为:
[0030] n取1,2,…,N2
[0031] 其中,N2为动态有限冲激响应滤波器的阶数,等于2个周波电压或电流采样值的点数,N2=2*K2+1,K2即为步骤(1-3)中电压或电流采样数据的编号值的最大值;
[0032] 动态信号系数矩阵Cd的计算公式为:
[0033]
[0034] 其中,动态信号系数矩阵Cd的元素有两类,一类元素与信号的基准频率有关,其计算公式如下:
[0035] k取值-K2,-(K2-1),…,0,…,K2
[0036] 其中,c(k,c)为基准频率信号的泰勒级数展开系数行向量,由电压或电流基准频率信号在编号值为k的采样数据点处的二阶泰勒展开计算得到, 为行向量c(k,c)的共轭行向量,其中每一个元素都是向量c(k,c)的共轭值,基准频率e为自然对数的底数,j为虚数单位,K2为步骤(1-3)中2周波电压或电流采样数据的编号值的最大值,K2=(N2-1)/2,N2为2周波电压或电流采样数据的点数,T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准角频率值,等于配电网基准频率与2π的乘积,
[0037] 动态信号系数矩阵Cd的第二类元素,与电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波有关,其计算公式如下:
[0038] k取值-K2,-(K2-1),…,0,…,K2
[0039] 其中,c(k,bi)为谐波和/或间谐波信号的泰勒级数展开系数行向量,由电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波信号在编号为k的采样数据点处的二阶泰勒展开计算得到, 为c(k,bi)的共轭行向量,其中每一个元素都是向量c(k,bi)的共轭值,
[0040] (5)根据步骤(4)中设计的稳态有限冲激响应滤波器Coeff_SteadyFilter和步骤(1)的4周波采样数据Xs,计算配电网的稳态同步相量,计算公式如下:
[0041] ps=Coeff_SteadyFilter·Xs
[0042] 其中,ps为计算得到的电压或电流信号的稳态相量值;
[0043] (6)根据步骤(4)中设计的动态有限冲激响应滤波器Coeff_dynamicFilter和步骤(1)的2周波采样数据Xd,计算配电网的动态同步相量pd,计算公式如下:
[0044] pd=Coeff_DynamicFilter·Xd
[0045] 其中, 其中, 为动态相量值, 为动态相量一阶导数值,为动态相量二阶导数值;
[0046] 根据上述动态相量值 计算出配电网的电压或电流信号动态相量的幅值A1和相位 计算公式如下:
[0047]
[0048]
[0049] 其中, 表示动态相量的模值,即动态相量的幅值,a0为动态相量幅值的零阶导数, 表示动态相量的相位,b0为动态相量相位的零阶导数值;
[0050] 根据上述动态相量值 及一阶导数值 计算得到配电网电压或电流信号的频率值fc,计算公式如下:
[0051]
[0052] 其中,f0为配电网电压或电流信号的额定频率值,b1为相位的一阶导数值:
[0053]
[0054] 其中,Im{}表示求取复数的虚部值,j为虚数单位;
[0055] 根据上述动态相量值 一阶导数值 和二阶导数值 计算得到配电网电压或电流信号的频率变化率值ROCOF,计算公式如下:
[0056]
[0057] 其中,b2为动态相量值 的相位的二阶导数值,计算公式如下:
[0058]
[0059] 其中,b1为动态相量的相位的一阶导数值,a1为动态相量的幅值的一阶导数,计算公式为:
[0060]
[0061] 其中,Re{}表示求取复数的实部值;
[0062] (7)设定配电网电压或电流信号状态判断规则:
[0063] (7-1)根据步骤(6)中计算得到的动态相量的幅值一阶导数a1和相位一阶导数b1,设定两个状态判断因子Index1和Index2,计算公式如下:
[0064] Index1=|a1|,Index2=|b1|
[0065] (7-2)根据步骤(5)中计算得到的电压或电流信号的稳态相量ps,设定两个状态判断因子Index3和Index4,计算公式如下:
[0066]
[0067]
[0068] 其中,K2为电压或电流信号的采样数据Xd中编号的最大值,x(k)为步骤(1)中采样数据Xd中编号为k的采样值, 为采样数据Xd中编号为k的采样值相对应的估计值,计算公式如下:
[0069]
[0070] 其中,|ps|为稳态相量值ps的模值,∠ps为稳态相量值ps的相位值,fc为步骤(6)计算得到的配电网电压或电流信号的频率值;
[0071] (7-3)对配电网电压或电流信号当前的状态进行判断,若满足:
[0072] {Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=1[0073] 其中,σth1、σth2、σth2和σth2分别是Index1、Index2、Index3和Index4四个状态判断因子的边界值,∪为并集符号,表示“或”的关系,∩为取交集,表示“与”的关系,则判定配电网处于动态,选择步骤(6)中计算得到的动态相量pd作为最终得到的同步相量计算值P;
[0074] 若满足:
[0075] {Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=0[0076] 则判定配电网处于稳态,将步骤(5)中计算得到的稳态相量ps作为最终得到的同步相量计算值P,实现配电网的自适应同步相量量测。
[0077] 本发明提出的配电网的自适应同步相量量测方法,其优点是:
[0078] 本发明方法采用稳态和动态相量计算两个通道并行计算的方式,并融合信号频率成分分析和配电网暂态变化快速识别技术,实现稳态和动态场景下同步相量的可靠计算,兼顾了同步相量计算的精确性和快速响应性。目前算法的同步相量在稳态基波、频率偏差±5Hz、含10%的2~25次谐波、含10%的25Hz~150Hz带外干扰的稳态量测场景下以及调制幅值为10%,调制频率为5Hz的动态场景下,幅值计算误差均小于0.2%,相角计算精度均小于0.05度,阶跃响应时间小于30毫秒。附图说明
[0079] 图1是本发明提出的配电网的自适应同步相量量测方法的流程框图

具体实施方式

[0080] 本发明提出的配电网的自适应同步相量量测方法,其流程框图如图1所示,包括以下步骤:
[0081] (1)采集配电网的电压或电流信号,包括以下步骤:
[0082] (1-1)以时间间隔T采集配电网的电压或电流信号,将当前采样数据点的前10个周波采样数据记为X;
[0083] (1-2)取采样数据X中当前采样数据点前的4个周波数据,记为Xs,采样数据的个数记为N1,将采样数据分别编号-K1,-(K1-1),……,0,1,……,(K1-1),K1;
[0084] (1-3)取采样数据X中当前采样数据点前的2个周波数据,记为Xd,采样数据个数记为N2,将采样数据分别编号-K2,-(K2-1),……,0,1,……,(K2-1),K2;
[0085] (2)采用现代谱估计方法,确定步骤(1)的采样数据X中电压或电流信号的频率在25Hz~150Hz范围内的谐波和/或间谐波的数目和频率值;
[0086] (3)根据步骤(1-2)周波数据Xs和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个稳态有限冲激响应滤波器,该滤波器的系数为Coeff_SteadyFilter:
[0087] Coeff_SteadyFilter=Fs(1,:)
[0088] 其中,Fs为稳态滤波器系数矩阵,将稳态滤波器系数矩阵Fs中的第一行作为稳态有限冲激响应滤波器的系数Coeff_SteadyFilter,稳态滤波器系数矩阵Fs由稳态信号系数矩阵Cs和稳态滤波器权重矩阵W1计算得到,其中:
[0089] Fs=(CHsW1HW1Cs)-1CsW1
[0090] 稳态信号系数矩阵Cs的计算公式为:
[0091]
[0092] 其中,上标H为矩阵共轭转置,e为自然对数的底数,j为虚数单位,K1为步骤(1-2)中电压或电流采样数据的编号值的最大值,T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准角频率值,ωc等于配电网基准频率与2π的乘积,ωbi为电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波的角频率值,ωbi等于电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波的频率值与2π的乘积;
[0093] 稳态滤波器权重矩阵W1为对角矩阵,即除了矩阵对角线上元素不为零外,矩阵其余元素均为零,计算公式为:
[0094]
[0095] 其中,w1,w2,……,wN1为稳态滤波器权重矩阵W1的对角元素,计算公式为:
[0096] n取1,2,…,N1
[0097] 其中,N1为所设计的稳态有限冲激响应滤波器的阶数,该阶数等于4个周波电压或电流采样数据的点数,N1=2*K1+1,K1即为4周波电压或电流采样数据的编号值的最大值,Cosh表示双曲余弦函数,ac和ρm分别是窗函数的特性调节系数,ac的取值为7.0,ρm的取值为1.2;
[0098] (4)根据步骤(1-3)周波数据Xd和步骤(2)的谐波和/或间谐波的数目和频率值,设计一个动态有限冲激响应滤波器,该滤波器的系数为Coeff_dynamicFilter,计算公式为:
[0099] Coeff_DynamicFilter=Fd(1:3,:)
[0100] 其中,Fd为动态滤波器系数矩阵,将动态滤波器系数矩阵Fd中的第一行至第三行作为动态有限冲激响应滤波器的系数Coeff_dynamicFilter,动态滤波器系数矩阵Fd由动态信号系数矩阵Cd和权重矩阵W2计算得到,计算公式为:
[0101] Fd=(CHdW2HW2Cd)-1CdW2
[0102] 其中,动态滤波器权重矩阵W2的计算公式如下:
[0103]
[0104] 其中,w1,w2,……,wN2为稳态滤波器权重矩阵W2的对角元素,计算公式为:
[0105] n取1,2,…,N2
[0106] 其中,N2为动态有限冲激响应滤波器的阶数,等于2个周波电压或电流采样值的点数,N2=2*K2+1,K2即为步骤(1-3)中电压或电流采样数据的编号值的最大值;
[0107] 动态信号系数矩阵Cd的计算公式为:
[0108]
[0109] 其中,动态信号系数矩阵Cd的元素有两类,一类元素与信号的基准频率有关,其计算公式如下:
[0110] k取值-K2,-(K2-1),…,0,…,K2
[0111] 其中,c(k,c)为基准频率信号的泰勒级数展开系数行向量,由电压或电流基准频率信号在编号值为k的采样数据点处的二阶泰勒展开计算得到。 为行向量c(k,c)的共轭行向量,其中每一个元素都是向量c(k,c)的共轭值。基准频率e为自然对数的底数,j为虚数单位,K2为步骤(1-3)中2周波电压或电流采样数据的编号值的最大值,K2=(N2-1)/2,N2为2周波电压或电流采样数据的点数。T为各采样数据之间的时间间隔,ωc为电压或电流信号的基准角频率值,等于配电网基准频率与2π的乘积。
[0112] 动态信号系数矩阵Cd的第二类元素,与电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波有关,其计算公式如下:
[0113] k取值-K2,-(K2-1),…,0,…,K2
[0114] 其中,c(k,bi)为谐波和/或间谐波信号的泰勒级数展开系数行向量,由电压或电流信号中频率在25Hz~150Hz范围内的谐波和/或间谐波信号在编号为k的采样数据点处的二阶泰勒展开计算得到。 为c(k,bi)的共轭行向量,其中每一个元素都是向量c(k,bi)的共轭值。
[0115] (5)根据步骤(4)中设计的稳态有限冲激响应滤波器Coeff_SteadyFilter和步骤(1)的4周波采样数据Xs,计算配电网的稳态同步相量,计算公式如下:
[0116] ps=Coeff_SteadyFilter·Xs
[0117] 其中,ps为计算得到的电压或电流信号的稳态相量值;
[0118] (6)根据步骤(4)中设计的动态有限冲激响应滤波器Coeff_dynamicFilter和步骤(1)的2周波采样数据Xd,计算配电网的动态同步相量pd,计算公式如下:
[0119] pd=Coeff_DynamicFilter·Xd
[0120] 其中, 其中, 为动态相量值, 为动态相量一阶导数值,为动态相量二阶导数值;
[0121] 根据上述动态相量值 计算出配电网的电压或电流信号动态相量的幅值A1和相位 计算公式如下:
[0122]
[0123]
[0124] 其中, 表示动态相量的模值,即动态相量的幅值,a0为动态相量幅值的零阶导数, 表示动态相量的相位,b0为动态相量相位的零阶导数值;
[0125] 根据上述动态相量值 及一阶导数值 计算得到配电网电压或电流信号的频率值fc,计算公式如下:
[0126]
[0127] 其中,f0为配电网电压或电流信号的额定频率值,b1为相位的一阶导数值:
[0128]
[0129] 其中,Im{}表示求取复数的虚部值,j为虚数单位;
[0130] 根据上述动态相量值 一阶导数值 和二阶导数值 计算得到配电网电压或电流信号的频率变化率值ROCOF,计算公式如下:
[0131]
[0132] 其中,b2为动态相量值 的相位的二阶导数值,计算公式如下:
[0133]
[0134] 其中,b1为动态相量的相位的一阶导数值,a1为动态相量的幅值的一阶导数,计算公式为:
[0135]
[0136] 其中,Re{}表示求取复数的实部值;
[0137] (7)设定配电网电压或电流信号状态判断规则:
[0138] (7-1)根据步骤(6)中计算得到的动态相量的幅值一阶导数a1和相位一阶导数b1,设定两个状态判断因子Index1和Index2,计算公式如下:
[0139] Index1=|a1|,Index2=|b1|
[0140] (7-2)根据步骤(5)中计算得到的电压或电流信号的稳态相量ps,设定两个状态判断因子Index3和Index4,计算公式如下:
[0141]
[0142]
[0143] 其中,K2为电压或电流信号的采样数据Xd中编号的最大值,x(k)为步骤(1)中采样数据Xd中编号为k的采样值, 为采样数据Xd中编号为k的采样值相对应的估计值,计算公式如下:
[0144]
[0145] 其中,|ps|为稳态相量值ps的模值,∠ps为稳态相量值ps的相位值,fc为步骤(6)计算得到的配电网电压或电流信号的频率值;
[0146] (7-3)对配电网电压或电流信号当前的状态进行判断,若满足:
[0147] {Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=1[0148] 其中,σth1、σth2、σth2和σth2分别是Index1、Index2、Index3和Index4四个状态判断因子的边界值,本发明方法的一个实施例中,边界值的取值分别为:0.001、0.001、0.00001、0.00001,∪为并集符号,表示“或”的关系,∩为取交集,表示“与”的关系,则判定配电网处于动态,选择步骤(6)中计算得到的动态相量pd作为最终得到的同步相量计算值P;
[0149] 若满足:
[0150] {Index1>σth1}∪{Index2>σth2}∩{Index3>σth3}∩{Index4>σth4}=0[0151] 则判定配电网处于稳态,将步骤(5)中计算得到的稳态相量ps作为最终得到的同步相量计算值P,实现配电网的自适应同步相量量测。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈