CN108614154B - 一种基于频谱分析的电力用户相位识别方法 - Google Patents
一种基于频谱分析的电力用户相位识别方法 Download PDFInfo
- Publication number
- CN108614154B CN108614154B CN201810498584.8A CN201810498584A CN108614154B CN 108614154 B CN108614154 B CN 108614154B CN 201810498584 A CN201810498584 A CN 201810498584A CN 108614154 B CN108614154 B CN 108614154B
- Authority
- CN
- China
- Prior art keywords
- phase
- htm
- user
- meter
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种基于频谱分析的电力用户相位识别方法,采用频谱分析技术,如傅利叶变换和反变换或小波分解和重构,对变压器台区下的用户电表用电量数据和变压器每个相位用电量数据进行分析计算,获得其时域数据高频分量的差值,或用电量差值的时域数据高频分量,再计算用户电表的特征值,排序后取前若干个,提取对应周期的时域数据高频分量的差值或用电量差值的时域数据高频分量构成特征集合,最后根据特征集合利用相关性分析方法获得超过指定阈值的最大的相似度值,以此判断用户电表所属相位。本发明无需对用户智能电表和抄表方式进行任何改造,只需用电量数据即可较为准确的判断用户所属相位,且用电量数据越多,识别率越高,便于推广应用。
Description
技术领域
本发明属于电气工程(电力***)领域,特别涉及一种基于频谱分析的电力用户相位识别方法。
背景技术
发电厂产生三相交流电压,为向家庭或建筑物供电,需通过变压器逐步降压。当最后分配电力到低压用户时,每个低压变压器供电约几十户到几百户。然而,低压用户使用的电力是单相的,每个低压用户都连接到低压变压器的三相之一:A相,B相或C相。
低压变压器三个相位上连接的低压用户负荷随时间变化不断波动,如果三相负荷差异较大,即三相负载不对称,会造成变压器损耗增大,严重时会危及用电安全。因此,需要定期检查变压器三相用户的负荷情况,必要时需要对低压用户的相位进行调整。
对低压用户相位进行调整的前提是知道低压用户所属的相位。目前国内的低压用户电表抄表方式包括有线和无线两大类,无线抄表的用户占低压用户总数的三分之一以上。无线抄表的低压用户无法获得其所属的变压器相位,如果将无线抄表方式更换为有线抄表方式,成本巨大。因此,如何识别无线抄表低压用户的相位,对变压器的稳定运行至关重要。
现有的技术方案,多采用人力现场观察确定或采用相位识别仪等专业仪器。使用人力现场观察确定,存在工程量巨大,耗时耗力,检查的准确性和效率较差等问题;而用专业的相位识别仪器多数基于载波通信技术和电流脉冲技术。对于载波通信技术而言,由于载波信号会在台变或电缆内跨相耦合,使得在某一相加载信息后,其他两相也可以接收到,即存在“信号串相”问题。对于电流脉冲技术而言,需要在变压器出线端安装电流互感器,操作上存在一定安全隐患,且可控制性差。而且使用专业仪器会增加电力监控管理设备的投入成本,增加电网的运维管理投入。所以目前亟需一种操作简单,且行之有效的相位识别方法。
公开号为107147430A的专利提供了一种基于载波通信的电力相位识别方法。方法要求一报文发送设备,存在硬件及其维护成本,同时也存在部署工程量大等问题。
公开号为107271946A的专利提供了一种利用智能电能表采集的电压与一致相位电能表的几个时刻的电压值做线性相关运算来确定相位的电能表相位识别方法。该方法要求用户的智能电表必须采集电压数据,并且对时间的精度要求较高,即需要同时采集变压器关口表和用户电表的电压读数,通过电压的相关性进行判断。国内低压用户的智能电表虽然具有采集电压读数的功能,但在电力***的实际运行中,该数据并未采集和保存,低压用户智能电表采集并保存的数据通常只有用电量(电表读数)。并且,大多数低压用户智能电表也无法保证数据采集时间的准确性。总体来说,该方法缺乏实用性。
发明内容
为克服上述现有技术的缺点,本发明的目的在于提供一种基于频谱分析的电力用户相位识别方法,该方法只需采集变压器关口表每个相位和低压用户智能电表的用电量,即可较为准确的识别低压用户所属变压器相位。采集的用电量数据越多,识别率越高。
为了实现上述目的,本发明采用的技术方案是:
一种基于频谱分析的电力用户相位识别方法,包括以下步骤:
1)采集第i周期变压器关口表电表每一相位p的电量数据tmi,p,p∈{A,B,C},和每一用户电表j的电量数据umi,j,j=1,2,...,K,按T个周期汇总得如下集合:
TMp={tmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
UMj={umi,j,i∈N且i∈[1,T]},j=1,2,...,K
其中,i表示周期序号,A、B、C表示变压器关口表的三个相位,K为用户电表的数量,T表示周期个数,TMp表示变压器关口表p相位的电量数据集合,UMj表示用户电表j的电量数据集合,N表示自然数;
2)逐一对每个相位p按如下公式计算△tmi,p:
△tmi,p=tmi,p-tmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△TMp={△tmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△umi,j:
△umi,j=umi,j-umi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△UMj={△umi,j,i∈N且i∈[2,T]},j=1,2,...,K
3)逐一对每个相位p,用频谱分析方法处理△TMp,获得如下形式的时域数据高频分量集合△TMHp:
△TMHp={△tmhi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理△UMj,获得如下形式的时域数据高频分量集合△UMHj:
△UMHj={△umhi,j,i∈N且i∈[2,T]},j=1,2,...,K
4)按如下公式逐一计算用户电表j的特征值fi,j:
逐一对每个用户电表j,按fi,j的大小降序排序并保留前X个,分别在集合△TMHA、△TMHB、△TMHC和△UMHj中依次找出特征值所对应周期序号的元素并构成特征集合△CA,j、△CB,j、△CC,j和△Cj;
其中,所述X为提取的特征值的个数;
5)逐一对每个用户电表j,采用相似度计算方法分别对△CA,j与△Cj、△CB,j与△Cj、△CC,j与△Cj进行计算,获得三个相位上的相似度值SA,j、SB,j和SC,j,取SA,j、SB,j和SC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位;如否,则无法判断所属相位。
在上述方案中,所述步骤2),步骤3),步骤4)和步骤5)可以如下步骤2a),步骤3a),步骤4a)和步骤5a)分别代替:
步骤2a)
逐一对每个相位p,用频谱分析方法处理TMp,获得如下形式的时域数据高频分量集合HTMp:
HTMp={htmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理UMj,获得如下形式的时域数据高频分量集合HUMj:
HUMj={humi,j,i∈N且i∈[1,T]},j=1,2,...,K
步骤3a)
逐一对每个相位p按如下公式计算△htmi,p:
△htmi,p=htmi,p-htmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HTMp={△htmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△humi,j:
△humi,j=humi,j-humi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HUMj={△humi,j,i∈N且i∈[2,T]},j=1,2,...,K
步骤4a)
按如下公式逐一计算用户电表j的特征值hfi,j:
逐一对每个用户电表j,按hfi,j的大小降序排序并保留前HX个,分别在集合△HTMA、△HTMB、△HTMC和△HUMj中依次找出特征值所对应周期序号的元素并构成特征集合△HCA,j、△HCB,j、△HCC,j和△HCj;
其中,所述HX为提取的特征值的个数;
步骤5a)
逐一对每个用户电表j,采用相似度计算方法分别对△HCA,j与△HCj、△HCB,j与△HCj、△HCC,j与△HCj进行计算,获得三个相位上的相似度值HSA,j、HSB,j和HSC,j,取HSA,j、HSB,j和HSC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位,如否,则无法判断所属相位。
所述频谱分析的方法为傅里叶变换和反变换。
所述频谱分析的方法为小波分解和重构。
所述相似度计算方法为Pearson相关性分析方法或余弦相似计算方法或hash相似度计算方法,对应的相似度为相关系数。
与现有技术相比,本发明无需对用户智能电表和抄表方式进行任何改造,不存在任何改造成本;本发明所提供的方法只需采集变压器关口表每个相位和低压用户智能电表的用电量,这些数据都是电力***肯定会采集和保存的,因此实用性极高;与传统方式相比较,本发明提供的方法只需要通过计算即可对低压用户所属相位进行较为准确的识别,降低了相位识别的工作量和成本;本发明提供的方法只需在现有数据的基础上进行计算即可获得结果,可迅速见效,效率极高并便于推广应用。
具体实施方式
下面结合实施例详细说明本发明的实施方式,本发明一共有四类实施例。
一、实施例(一)
本发明的具体实施例(一)步骤如下:
步骤1),采集第i周期变压器关口表电表每一相位p的电量数据tmi,p,p∈{A,B,C},和每一用户电表j的电量数据umi,j,j=1,2,...,K,按T个周期汇总可得如下集合:
TMp={tmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
UMj={umi,j,i∈N且i∈[1,T]},j=1,2,...,K
举实例如下:
采集总周期T=10000的某变压器关口表和该关口表下10个用户电表的电量数据,电量数据单位为kwh,周期时间间隔为半小时,用从1到10的数字编号代表10个用户电表,用A,B,C代表关口表的三个相位。得到数据如下表1-1所示;
表1-1 采集数据部分举例表(受限于篇幅未能全部罗列)
在实际中,关口表不一定直接记录每个相位的电量数据,但会记录每个相位的电压值、电流值,因此某一相位的电量数据可以通过对该相位的电压值、电流值、周期时间进行计算获得近似值。
关口表A相位集合TMA={tmi,A,i∈N且i∈[1,10000]},其余相位以此类推;
用户电表1的电量数据集合UM1={umi,1,i∈N且i∈[1,10000]},其余用户电表以此类推。
步骤2)逐一对每个相位p按如下公式计算△tmi,p:
△tmi,p=tmi,p-tmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△TMp={△tmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△umi,j:
△umi,j=umi,j-umi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△UMj={△umi,j,i∈N且i∈[2,T]},j=1,2,...,K
以步骤1)例中的TMA和UM1为例,按相应的公式进行计算可得如下表1-2所示数据。
表1-2 以TMA和UM1为例的差值数据
关口表A相位集合△TMA={△tmi,A,i∈N且i∈[2,10000]},其余相位以此类推;
用户电表1的电量数据集合△UM1={△umi,1,i∈N且i∈[2,10000]},其余用户电表以此类推。
步骤3)逐一对每个相位p,用频谱分析方法处理△TMp,获得如下形式的时域数据高频分量集合△TMHp:
△TMHp={△tmhi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理△UMj,获得如下形式的时域数据高频分量集合△UMHj:
△UMHj={△umhi,j,i∈N且i∈[2,T]},j=1,2,...,K
以步骤2)中的△TMA为例,用频谱分析方法处理△TMA:
首先,对△TMA进行离散傅里叶变换,可以得到下表1-3的数据;
表1-3 对△TMA进行傅里叶变换的数据
随后,可设定截止频率fc=1000,将变换后的数据中高于该频率位置的幅值置零,如下表1-4;
表1-4 对△TMA进行傅里叶变换后再滤波的数据
然后,根据表1-4将滤波后的数据进行傅里叶反变换,取实数部分,此时得到的是低频信息,如下表1-5;
表1-5 将滤波后的数据反变换后的数据
i=2 | i=3 | ... | i=10000 | |
低频信息数据 | -0.7454 | -0.9833 | -0.4877 |
最后,由原始数据减去低频信息数据即为时域数据高频分量,得下表1-6:
表1-6 △TMA经频谱分析后的时域数据高频分量
其余的相位同理可以获得时域数据高频分量,对用户电表也是同样的处理,在整个相位识别过程中的截止频率应当保持一致。对所有相位数据和用户电表数据都处理后,得到下表1-7;
表1-7 傅里叶方法的时域数据高频分量表
汇总可得相位A的时域数据高频分量集合,如下所示,其余相位以此类推。
△TMHA={△tmhi,A,i∈N且i∈[2,10000]}
汇总可得用户电表1的时域数据高频分量集合,如下所示,其余用户电表以此类推。
△UMH1={△umhi,1,i∈N且i∈[2,10000]}
步骤4),按如下公式逐一计算用户电表j的特征值fi,j:
逐一对每个用户电表j,按fi,j的大小降序排序并保留前X个,分别在集合△TMHA、△TMHB、△TMHC和△UMHj中依次找出特征值所对应周期序号的元素并构成特征集合△CA,j、△CB,j、△CC,j和△Cj;
对表1-7中用户电表的有关数据,按下述公式逐一计算用户电表j的特征值fi,j,得到下表1-8:
表1-8 用户电表的特征值表
以表1-8中j=1为例,首先按fi,1的大小降序排序并保留前X=100个,分别在集合△TMHA、△TMHB、△TMHC和△UMH1中依次找出对应周期序号的元素值构成集合△CA,1、△CB,1、△CC,1和△C1。如下表1-9所示:
表1-9 j=1时对应周期序号的△CA,1、△CB,1、△CC,1和△C1
表1-9中,第三行的数据△tmh4056,A、△tmh8758,A、...、△tmh7921,A构成△CA,1,以此类推。
步骤5),逐一对每个用户电表j,采用相似度计算方法分别对△CA,j与△Cj、△CB,j与△Cj、△CC,j与△Cj进行计算,获得三个相位上的相似度值SA,j、SB,j和SC,j,取SA,j、SB,j和SC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位,如否,则无法判断所属相位。
以用户电表1为例,取表1-9中数据分别计算三个Pearson相似度SA,1、SB,1和SC,1,以计算SA,1为例:
分别计算△C1和△CA,1的平均值和方差,见表1-10:
表1-10 平均值方差表
其中,△ck,1和△ck,A,1分别表示△C1和△CA,1中的第k个元素。计算Pearson相似度:
同理获得:
SB,1=0.0438
SC,1=-0.0077
最大相似度为SA,1=0.1850>SB,1>SC,1,且SA,1>LS=0.1,其中LS为设定的阈值,所以判别用户电表1属于相位A。其余用户电表同理判别。
二、实施例(二)
本实施例与实施例(一)的不同之处在于步骤3),将实施例(一)中步骤3)的频谱分析采用小波分解和重构作为实施方法。
以△UM1为例:基波函数采用dB2函数,分解层数为4层,得到的数据如下表2-1:
表2-1 对△UM1用小波分解和重构后的数据
取D1+D2+D3+D4作为重构后的△UM1的时域数据高频分量。
同理应用于其他行数据可以得到下表2-2:
表2-2 用小波分解和重构方法的时域数据高频分量表
用表2-2代替表1-7参与实施例(一)之后的步骤3),步骤4)和步骤5)的运算即可。
三、实施例(三)
本实施例与实施例(一)比较而言,其步骤1)相同,都为采集数据,但之后对数据的处理不同,以表1-1的数据为例。
步骤2a)逐一对每个相位p,用频谱分析方法处理TMp,获得如下形式的时域数据高频分量集合HTMp:
HTMp={htmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理UMj,获得如下形式的时域数据高频分量集合HUMj:
HUMj={humi,j,i∈N且i∈[1,T]},j=1,2,...,K
对表1-1中的每行数据进行傅里叶变换,以TMA为例,进行傅里叶变换,获得其频域的值,如下表3-1:
表3-1 对TMA进行傅里叶变换的数据
在表3-1中,可设置截止频率fc=1000,将傅里叶变换后数据中高于该频率位置的幅值置零,得到如下表3-2所示:
表3-2 对TMA进行傅里叶变换后再滤波的数据
然后,根据表3-2将滤波后的数据进行傅里叶反变换,取实数部分,此时得到的是低频信息,由原始数据减去低频信息数据即为时域数据高频分量,如下表3-3:
表3-3 经频谱分析后的时域数据高频分量
其余的相位同理可以获得时域数据高频分量,对用户电表也是同样的处。对所有相位数据和用户电表数据都处理后。可以得到下表3-4:
表3-4 时域数据高频分量表
汇总可得相位A的时域数据高频分量集合,如下所示,其余相位以此类推。
HTMA={htmi,A,i∈N且i∈[1,10000]}
汇总可得用户电表1的时域数据高频分量集合,如下所示,其余用户电表以此类推。
HUM1={humi,1,i∈N且i∈[1,10000]}
步骤3a)逐一对每个相位p按如下公式计算△htmi,p:
△htmi,p=htmi,p-htmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HTMp={△htmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△humi,j:
△humi,j=humi,j-humi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HUMj={△humi,j,i∈N且i∈[2,T]},j=1,2,...,K
对于步骤2a)中的结果,以HTMA和HUM1为例可得到下表3-5。
表3-5 以HTMA和HUM1为例得到的差值数据
关口表A相位集合△HTMA={△htmi,A,i∈N且i∈[2,10000]},其余相位以此类推;
用户电表1的电量数据集合△HUM1={△humi,1,i∈N且i∈[2,10000]},其余用户电表以此类推。
下表3-6列出关口表相位和用户电表的部分的差值数据。
表3-6 差值数据表
步骤4a)按如下公式逐一计算用户电表j的特征值hfi,j:
逐一对每个用户电表j,按hfi,j的大小降序排序并保留前HX个,分别在集合△HTMA、△HTMB、△HTMC和△HUMj中依次找出特征值所对应周期序号的元素并构成特征集合△HCA,j、△HCB,j、△HCC,j和△HCj;
其中,所述HX为提取的特征值的个数。
对步骤3a)中用户电表的有关数据,按下述公式逐一计算用户电表j的特征值hfi,j,得到下表3-7:
表3-7 用户电表的特征值表
以表3-7中j=1为例,首先按hfi,1的大小降序排序并保留前HX=100个,分别在集合△HTMA、△HTMB、△HTMC和△HUM1中依次找出对应周期序号的元素值构成集合△HCA,1、△HCB,1、△HCC,1和△HC1。如下表3-8所示:
表3-8 j=1时对应周期的△HCA,1、△HCB,1、△HCC,1和△HC1
步骤5a)逐一对每个用户电表j,采用相似度计算方法分别对△HCA,j与△HCj、△HCB,j与△HCj、△HCC,j与△HCj进行计算,获得三个相位上的相似度值HSA,j、HSB,j和HSC,j,取HSA,j、HSB,j和HSC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位,如否,则无法判断所属相位。
以用户电表1为例,取表3-8中数据分别计算三个Pearson相似度HSA,1、HSB,1和HSC,1,以计算HSA,1为例:
分别计算△HC1和△HCA,1的平均值和方差,见表3-9:
表3-9 平均值方差表
其中,△hck,1和△hck,A,1分别表示△HC1和△HCA,1中的第k个元素。计算Pearson相似度:
同理获得
HSB,1=0.0437
HSC,1=-0.0079
最大相似度为HSA,1=0.1848>HSB,1>HSC,1,且HSA,1>LS=0.1,其中LS为设定的阈值,所以用户电表1属于相位A。
四、实施例(四)
本实施例与实施例(三)的不同之处在于步骤2a),将实施例(三)中步骤2a)的频谱分析采用小波分解与重构作为实施方法。
以UM1为例:基波函数采用dB2函数,分解层数为4层,得到的数据如下表4-1:
表4-1 对UM1用小波分解和重构后的数据
取D1+D2+D3+D4作为重构后的UM1的时域数据高频分量。
同理应用于其他相位和用户电表可以得到下表4-2:
表4-2 用小波分解和重构方法的时域数据高频分量表
用表4-2代替表3-4参与实施例(三)之后的步骤3a),步骤4a)和步骤5a)的运算即可。
Claims (5)
1.一种基于频谱分析的电力用户相位识别方法,其特征在于,包括以下步骤:
1)采集第i周期变压器关口表电表每一相位p的电量数据tmi,p,p∈{A,B,C},和每一用户电表j的电量数据umi,j,j=1,2,...,K,按T个周期汇总得如下集合:
TMp={tmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
UMj={umi,j,i∈N且i∈[1,T]},j=1,2,...,K
其中,i表示周期序号,A、B、C表示变压器关口表的三个相位,K为用户电表的数量,T表示周期个数,TMp表示变压器关口表p相位的电量数据集合,UMj表示用户电表j的电量数据集合,N表示自然数;
2)逐一对每个相位p按如下公式计算△tmi,p:
△tmi,p=tmi,p-tmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△TMp={△tmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△umi,j:
△umi,j=umi,j-umi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△UMj={△umi,j,i∈N且i∈[2,T]},j=1,2,...,K
3)逐一对每个相位p,用频谱分析方法处理△TMp,获得如下形式的时域数据高频分量集合△TMHp:
△TMHp={△tmhi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理△UMj,获得如下形式的时域数据高频分量集合△UMHj:
△UMHj={△umhi,j,i∈N且i∈[2,T]},j=1,2,...,K
4)按如下公式逐一计算用户电表j的特征值fi,j:
逐一对每个用户电表j,按fi,j的大小降序排序并保留前X个,分别在集合△TMHA、△TMHB、△TMHC和△UMHj中依次找出特征值所对应周期序号的元素并构成特征集合△CA,j、△CB,j、△CC,j和△Cj;
其中,所述X为提取的特征值的个数;
5)逐一对每个用户电表j,采用相似度计算方法分别对△CA,j与△Cj、△CB,j与△Cj、△CC,j与△Cj进行计算,获得三个相位上的相似度值SA,j、SB,j和SC,j,取SA,j、SB,j和SC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位;如否,则无法判断所属相位。
2.根据权利要求1所述基于频谱分析的电力用户相位识别方法,其特征在于,利用如下步骤2a),步骤3a),步骤4a)和步骤5a)分别代替所述步骤2),步骤3),步骤4)和步骤5)
步骤2a)
逐一对每个相位p,用频谱分析方法处理TMp,获得如下形式的时域数据高频分量集合HTMp:
HTMp={htmi,p,i∈N且i∈[1,T]},p∈{A,B,C}
逐一对每个用户电表j,用频谱分析方法处理UMj,获得如下形式的时域数据高频分量集合HUMj:
HUMj={humi,j,i∈N且i∈[1,T]},j=1,2,...,K
步骤3a)
逐一对每个相位p按如下公式计算△htmi,p:
△htmi,p=htmi,p-htmi-1,p,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HTMp={△htmi,p,i∈N且i∈[2,T]},p∈{A,B,C}
逐一对每个用户电表j按如下公式计算△humi,j:
△humi,j=humi,j-humi-1,j,i∈N且i∈[2,T]
按T-1个周期汇总可得如下集合:
△HUMj={△humi,j,i∈N且i∈[2,T]},j=1,2,...,K
步骤4a)
按如下公式逐一计算用户电表j的特征值hfi,j:
i∈N且i∈[2,T]
逐一对每个用户电表j,按hfi,j的大小降序排序并保留前HX个,分别在集合△HTMA、△HTMB、△HTMC和△HUMj中依次找出特征值所对应周期序号的元素并构成特征集合△HCA,j、△HCB,j、△HCC,j和△HCj;
其中,所述HX为提取的特征值的个数;
步骤5a)
逐一对每个用户电表j,采用相似度计算方法分别对△HCA,j与△HCj、△HCB,j与△HCj、△HCC,j与△HCj进行计算,获得三个相位上的相似度值HSA,j、HSB,j和HSC,j,取HSA,j、HSB,j和HSC,j中最大的值并判断是否超过阈值;如是,则该最大值所对应相位即为用户电表j所属相位,如否,则无法判断所属相位。
3.根据权利要求1或2所述基于频谱分析的电力用户相位识别方法,其特征在于,所述频谱分析的方法为傅里叶变换和反变换。
4.根据权利要求1或2所述基于频谱分析的电力用户相位识别方法,其特征在于,所述频谱分析的方法为小波分解和重构。
5.根据权利要求1或2所述基于频谱分析的电力用户相位识别方法,其特征在于,所述相似度计算方法为Pearson相关性分析方法或余弦相似计算方法或hash相似度计算方法,对应的相似度为相关系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810498584.8A CN108614154B (zh) | 2018-05-23 | 2018-05-23 | 一种基于频谱分析的电力用户相位识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810498584.8A CN108614154B (zh) | 2018-05-23 | 2018-05-23 | 一种基于频谱分析的电力用户相位识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108614154A CN108614154A (zh) | 2018-10-02 |
CN108614154B true CN108614154B (zh) | 2020-01-21 |
Family
ID=63663731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810498584.8A Active CN108614154B (zh) | 2018-05-23 | 2018-05-23 | 一种基于频谱分析的电力用户相位识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108614154B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146743B (zh) * | 2019-04-30 | 2021-05-11 | 国网江苏省电力有限公司电力科学研究院 | 一种台区相位识别方法 |
CN112527863A (zh) * | 2020-12-10 | 2021-03-19 | 国网陕西省电力公司咸阳供电公司 | 一种基于电表数据相关性分析的线变关系识别方法 |
CN112730984B (zh) * | 2021-03-31 | 2021-07-13 | 国网江西省电力有限公司供电服务管理中心 | 一种基于智能电表的低压配电网相位识别方法 |
CN114118304B (zh) * | 2022-01-25 | 2022-06-24 | 佰聆数据股份有限公司 | 用户电能表相位识别方法、***及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1633368A1 (ru) * | 1988-06-15 | 1991-03-07 | Запорожский машиностроительный институт им.В.Я.Чубаря | Цифровой измеритель показателей качества электрической энергии трехфазной сети |
CN101788614A (zh) * | 2009-12-28 | 2010-07-28 | 江苏省电力公司苏州供电公司 | 高低压综合相序相位仪及其检测方法 |
CN104092481A (zh) * | 2014-07-17 | 2014-10-08 | 江苏林洋电子股份有限公司 | 一种通过电压特征区分台区和相别的方法 |
CN107271946A (zh) * | 2017-06-01 | 2017-10-20 | 宁波迦南智能电气股份有限公司 | 一种电能表相位识别方法 |
-
2018
- 2018-05-23 CN CN201810498584.8A patent/CN108614154B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1633368A1 (ru) * | 1988-06-15 | 1991-03-07 | Запорожский машиностроительный институт им.В.Я.Чубаря | Цифровой измеритель показателей качества электрической энергии трехфазной сети |
CN101788614A (zh) * | 2009-12-28 | 2010-07-28 | 江苏省电力公司苏州供电公司 | 高低压综合相序相位仪及其检测方法 |
CN104092481A (zh) * | 2014-07-17 | 2014-10-08 | 江苏林洋电子股份有限公司 | 一种通过电压特征区分台区和相别的方法 |
CN107271946A (zh) * | 2017-06-01 | 2017-10-20 | 宁波迦南智能电气股份有限公司 | 一种电能表相位识别方法 |
Non-Patent Citations (2)
Title |
---|
一种基于载波通信的电力相位识别方法;郑杰等;《电测与仪表》;20150818;第52卷(第16A期);77-79 * |
电力相位识别技术研究与实现;孙乐乐等;《机电信息》;20171231(第30期);90-91 * |
Also Published As
Publication number | Publication date |
---|---|
CN108614154A (zh) | 2018-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108614154B (zh) | 一种基于频谱分析的电力用户相位识别方法 | |
CN110956220B (zh) | 一种非侵入式家用电器负荷识别方法 | |
CN107025365B (zh) | 一种用于用户侧的非侵入负荷辨识方法 | |
CN111242391B (zh) | 用于电力负荷识别的机器学习模型训练方法及*** | |
CN105676028B (zh) | 一种基于模板匹配滤波的居民负荷用电识别方法 | |
CN109299134B (zh) | 一种非侵入式监测***下电器的识别方法 | |
CN106154040B (zh) | 一种用于计算电网单点等值时等效谐波阻抗方法 | |
CN104638664A (zh) | 风电厂次同步振荡保护方法及保护装置 | |
CN105137177A (zh) | 一种用于配电网单点监测的谐波电压责任计算告警方法 | |
CN110571806B (zh) | 一种用于配电网负荷类别的特征提取与辨识方法 | |
CN112819649A (zh) | 确定台区户变关系的方法及装置 | |
CN110703006B (zh) | 一种基于卷积神经网络的三相电能质量扰动检测方法 | |
CN105305437B (zh) | 电力负荷三重可信度匹配辨识方法 | |
CN111027886A (zh) | 一种考虑单位成本效用的低电压治理方案的评估方法 | |
WO2024104037A1 (zh) | 基于数学形态学和模式识别的直流电弧检测方法 | |
CN110231503A (zh) | 高损台区基于格兰杰因果检验的窃电用户识别定位方法 | |
Gao et al. | High-impedance arc fault modeling for distribution networks based on dynamic geometry dimension | |
CN113466548A (zh) | 一种基于相量测量技术的智能电表台区识别方法 | |
CN205786952U (zh) | 一种变压器绕组变形综合诊断*** | |
CN112182499A (zh) | 一种基于时序电量数据的低压配电网拓扑结构辨识方法 | |
CN112730984B (zh) | 一种基于智能电表的低压配电网相位识别方法 | |
CN113054747B (zh) | 一种配电台区非线性负荷的非侵入式检测方法 | |
CN108982958B (zh) | 一种基于小波变换的单通道变频用电信号盲源分离方法 | |
CN112949563B (zh) | 基于变点检测和改进knn算法的非侵入式负荷识别方法 | |
CN110277835B (zh) | 一种基于用电信息采集***的低压用户过电压风险监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |