CN109085530B - 一种基于游程特性的动态负荷信号分析方法 - Google Patents

一种基于游程特性的动态负荷信号分析方法 Download PDF

Info

Publication number
CN109085530B
CN109085530B CN201811124662.4A CN201811124662A CN109085530B CN 109085530 B CN109085530 B CN 109085530B CN 201811124662 A CN201811124662 A CN 201811124662A CN 109085530 B CN109085530 B CN 109085530B
Authority
CN
China
Prior art keywords
dynamic load
load current
run
current amplitude
run length
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
Application number
CN201811124662.4A
Other languages
English (en)
Other versions
CN109085530A (zh
Inventor
李博
曹敏
王恩
朱梦梦
唐标
刘斯扬
沈鑫
翟少磊
韩彤
雷岩团哏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Yunnan Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Yunnan Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Yunnan Power Grid Co Ltd filed Critical Electric Power Research Institute of Yunnan Power Grid Co Ltd
Priority to CN201811124662.4A priority Critical patent/CN109085530B/zh
Publication of CN109085530A publication Critical patent/CN109085530A/zh
Application granted granted Critical
Publication of CN109085530B publication Critical patent/CN109085530B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/04Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本申请提供了一种基于游程特性的动态负荷信号分析方法,利用动态负荷电流幅值向量,得到动态负荷电流幅值的二元游程序列,通过对二元游程序列的处理,并利用动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值,工作人员可利用概率密度值和自相关函数值对动态负荷电流进行随机特性分析,进而提高电能表误差测试的准确性。

Description

一种基于游程特性的动态负荷信号分析方法
技术领域
本申请涉及电力技术领域,尤其涉及一种基于游程特性的动态负荷信号分析方法。
背景技术
近年来,随着工业技术的发展和智能电网的不断提升,高电压、特大功率的冲击负荷设备,如电弧炉、电气化铁路等已广泛接入电网,这类设备在运行过程中频繁开启与关断,呈现出负荷功率短时冲击性、以及负荷电流快速随机波动的复杂特性,引起电能表出现较大动态误差,导致电能计量不准确。
因此,分析这类动态负荷的随机特性是确定模型结构和参数、建立电能表动态误差测试信号模型的基础。从而,研究特大功率动态负荷的随机变化特性,分析动态负荷的主要特征对电能表误差的影响,对提高系统运行经济性方面,提升智能电网中电能准确计量水平意义十分重大。
但是,现有的研究基本都是采用长时间间隔负荷信号测量数据的分析与建模方法,在大时间尺度上分析得到动态负荷特性和模型,需要对负荷信号长时间的测量和分析,目前缺少对动态负荷小时间尺度上随机特性的分析,以及相应的动态负荷信号随机特性研究,不能满足分析动态负荷的主要特征对电能表动态误差影响的要求。
发明内容
本申请提供了一种基于游程特性的动态负载信号分析方法,以解决目前缺少对动态负荷小时间尺度上随机特性的分析,以及相应的动态负荷信号随机特性研究,不能满足分析动态负荷的主要特征对电能表动态误差影响的要求的问题。
本申请提供了一种基于游程特性的动态负荷信号分析方法,所述方法包括:
创建动态负荷电流幅值向量A={α(0),α(1),…,α(n)},其中,α(n)为第n个工频周期的动态电流幅值,n=0,1,2,…,N-1;N为选取的动态负荷电流信号随机样本的工频周期数;
对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列;
对所述二元游程序列进行进行处理,得到正游程G1(τ)和负游程G-1(τ)的数量,τ为游程长度,且为大于或等于1的正整数;
根据所述正游程G1(τ)和负游程G-1(τ)的数量,利用如下公式,得到游程长度为τ的游程的总数量,
G(τ)=G-1(τ)+G1(τ),其中,G(τ)为游程长度为τ的游程的总数量,τ为大于或等于1的正整数;
利用所述游程长度为τ的游程的总数量,计算得到所述动态负荷电流幅值游程长度的期望值和方差;
根据所述期望值和方差,利用所述动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值。
进一步地,所述根据预设的条件,对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列包括:
利用下面的二值化条件,对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列,
Figure BDA0001812063900000021
其中,α(n)为第n个工频周期的动态电流幅值,n=0,1,2,…,N-1;N为选取的动态负荷电流信号随机样本的工频周期数;X(n)动态负荷电流幅值的游程阈值,其取决于所述动态负荷电流基波的上下包络线的平均值。
进一步地,所述对所述二元游程序列进行进行处理,得到正游程G1(τ)和负游程G-1(τ)的数量包括:
在所述二元游程序列中,将连续τ个α(n)值为1的动态电流幅值所构成的游程序列确定为正游程G1(τ),τ为大于或等于1的正整数,即游程长度;
将连续将连续τ个α(n)值为-1的动态电流幅值所构成的游程序列确定为负游程G-1(τ),τ为大于或等于1的正整数,即游程长度;
统计所述正游程G1(τ)和负游程G-1(τ)的数量。
进一步地,所述动态负荷电流幅值游程长度的概率密度函数为
Figure BDA0001812063900000022
其中,P(τ)为动态负荷电流幅值游程长度的概率密度函数;G(τ)为游程长度为τ的游程的总数量;τ为游程长度,且为大于或等于1的正整数,N为选取的动态负荷电流信号随机样本的工频周期数;。
进一步地,所述动态负荷电流幅值游程长度的自相关函数为
Figure BDA0001812063900000031
其中,σ2为所述动态负荷电流幅值游程长度的期望值方差,E和
Figure BDA0001812063900000032
为所述动态负荷电流幅值游程长度的期望值,η为游程时延;τ为游程长度,且为大于或等于1的正整数;G(τ)为游程长度为τ的游程的总数量。
由以上技术方案可知,本申请提供了一种基于游程特性的动态负荷信号分析方法,利用动态负荷电流幅值向量,得到动态负荷电流幅值的二元游程序列,通过对二元游程序列的处理,并利用动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值,工作人员可利用概率密度值和自相关函数值对动态负荷电流进行随机特性分析,进而提高电能表误差测试的准确性。
附图说明
为了更清楚地说明本申请的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请提供了一种基于游程特性的动态负荷信号分析方法的流程图。
具体实施方式
参见图1,本申请提供了一种基于游程特性的动态负荷信号分析方法,包括如下步骤:
步骤11:创建动态负荷电流幅值向量A={α(0),α(1),…,α(n)},其中,α(n)为第n个工频周期的动态电流幅值,n=0,1,2,…,N-1。
步骤12:对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列。
具体地,利用下面的二值化条件,对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列,
Figure BDA0001812063900000033
其中,α(n)为第n个工频周期的动态电流幅值,n=0,1,2,…,N-1;N为选取的动态负荷电流信号随机样本的工频周期数;X(n)动态负荷电流幅值的游程阈值,其取决于所述动态负荷电流基波的上下包络线的平均值。
X(n)是一个由动态负荷信号电流幅度决定的函数,随着电流包络的波动而不断变化。
步骤13:对所述二元游程序列进行进行处理,得到正游程G1(τ)和负游程G-1(τ)的数量,τ为游程长度,且为大于或等于1的正整数。
具体地,在所述二元游程序列中,将连续τ个α(n)值为1的动态电流幅值所构成的游程序列确定为正游程G1(τ),τ为大于或等于1的正整数,即游程长度;将连续将连续τ个α(n)值为-1的动态电流幅值所构成的游程序列确定为负游程G-1(τ),τ为大于或等于1的正整数,即游程长度;统计所述正游程G1(τ)和负游程G-1(τ)的数量。
步骤14:根据所述正游程G1(τ)和负游程G-1(τ)的数量,利用如下公式,得到游程长度为τ的游程的总数量;
G(τ)=G-1(τ)+G1(τ),其中,G(τ)为游程长度为τ的游程的总数量,τ为大于或等于1的正整数。
步骤15:利用所述游程长度为τ的游程的总数量,计算得到所述动态负荷电流幅值游程长度的期望值和方差。
具体地,动态负荷电流幅值游程长度的期望计算公式如下:
Figure BDA0001812063900000041
其中,E(τ)或
Figure BDA0001812063900000042
为动态负荷电流幅值游程长度的期望值;G(τ)为游程长度为τ的游程的总数量,τ为大于或等于1的正整数;N为选取的动态负荷电流信号随机样本的工频周期数。
动态负荷电流幅值游程长度的期望值是游程长度落在某点或某个区间的概率所构成的曲线,反映动态负荷信号电流幅值游程长度的集中趋势数。
动态负荷电流幅值游程长度的方差计算公式如下:
Figure BDA0001812063900000043
其中,σ2为所述动态负荷电流幅值游程长度的期望值方差,
Figure BDA0001812063900000044
为动态负荷电流幅值游程长度的期望值;G(τ)为游程长度为τ的游程的总数量,τ为大于或等于1的正整数;N为选取的动态负荷电流信号随机样本的工频周期数。
动态负荷电流幅值游程长度的期望值方差越大,表明电流动态项幅度游程长度的波动范围越大,动态负荷越不稳定。
步骤16:根据所述期望值和方差,利用所述动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值。
进一步地,所述动态负荷电流幅值游程长度的概率密度函数为
Figure BDA0001812063900000051
其中,P(τ)为动态负荷电流幅值游程长度的概率密度函数;G(τ)为游程长度为τ的游程的总数量;τ为游程长度,且为大于或等于1的正整数,N为选取的动态负荷电流信号随机样本的工频周期数。
动态负荷电流幅值游程长度的概率密度函数用来描述电流动态项幅度游程长度在某个确定取值点附近的概率。
进一步地,所述动态负荷电流幅值游程长度的自相关函数为
Figure BDA0001812063900000052
其中,σ2为所述动态负荷电流幅值游程长度的期望值方差,E和
Figure BDA0001812063900000053
为所述动态负荷电流幅值游程长度的期望值,η为游程时延,即游程长度的推移数;τ为游程长度,且为大于或等于1的正整数;G(τ)为游程长度为τ的游程的总数量。
自相关函数表示动态负荷电流幅度的游程长度随机起伏变化的快慢。
由以上技术方案可知,本申请提供了一种基于游程特性的动态负荷信号分析方法,利用动态负荷电流幅值向量,得到动态负荷电流幅值的二元游程序列,通过对二元游程序列的处理,并利用动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值,工作人员可利用概率密度值和自相关函数值对动态负荷电流进行随机特性分析,进而提高电能表误差测试的准确性。

Claims (5)

1.一种基于游程特性的动态负荷信号分析方法,其特征在于,所述方法包括:
创建动态负荷电流幅值向量A={α(0),α(1),…,α(n)},其中,α(n)为第n个工频周期的动态负荷电流幅值,n=0,1,2,…,N-1;N为选取的动态负荷电流信号随机样本的工频周期数;
对所述动态负荷电流幅值向量进行二值化处理,得到所述动态负荷电流幅值向量的二元游程序列;
对所述二元游程序列进行处理,得到正游程G1(τ)和负游程G-1(τ)的数量,τ为游程长度,且为大于或等于1的正整数;
根据所述正游程G1(τ)和负游程G-1(τ)的数量,利用如下公式,得到游程长度为τ的游程的总数量,
G(τ)=G-1(τ)+G1(τ),其中,G(τ)为游程长度为τ的游程的总数量,τ为大于或等于1的正整数;
利用所述游程长度为τ的游程的总数量,计算得到所述动态负荷电流幅值游程长度的期望值和方差;
根据所述期望值和方差,利用所述动态负荷电流幅值游程长度的概率密度函数和自相关函数,计算得到所述动态负荷电流幅值游程长度的概率密度值和自相关函数值。
2.如权利要求1所述的方法,其特征在于,根据预设的条件,对所述动态负荷电流幅值向量进行二值化处理,得到所述动态负荷电流幅值的二元游程序列包括:
利用下面的二值化条件,对所述动态负荷电流幅值向量进行二值化处理,得所述动态负荷电流幅值的二元游程序列,
Figure FDA0002739258880000011
其中,α(n)为第n个工频周期的动态负荷电流幅值,n=0,1,2,…,N-1;N为选取的动态负荷电流信号随机样本的工频周期数;X(n)动态负荷电流幅值的游程阈值,其取决于动态负荷电流基波的上下包络线的平均值。
3.如权利要求2所述的方法,其特征在于,所述对所述二元游程序列进行处理,得到正游程G1(τ)和负游程G-1(τ)的数量包括:
在所述二元游程序列中,将连续τ个α(n)值为1的动态负荷电流幅值所构成的游程序列确定为正游程G1(τ),τ为大于或等于1的正整数,即游程长度;
将连续τ个α(n)值为-1的动态负荷电流幅值所构成的游程序列确定为负游程G-1(τ),τ为大于或等于1的正整数,即游程长度;
统计所述正游程G1(τ)和负游程G-1(τ)的数量。
4.如权利要求1所述的方法,其特征在于,所述动态负荷电流幅值游程长度的概率密度函数为
Figure FDA0002739258880000021
其中,P(τ)为动态负荷电流幅值游程长度的概率密度函数;G(τ)为游程长度为τ的游程的总数量;τ为游程长度,且为大于或等于1的正整数,N为选取的动态负荷电流信号随机样本的工频周期数。
5.如权利要求1所述的方法,其特征在于,所述动态负荷电流幅值游程长度的自相关函数为
Figure FDA0002739258880000022
其中,σ2为所述动态负荷电流幅值游程长度的期望值方差,E和
Figure FDA0002739258880000023
为所述动态负荷电流幅值游程长度的期望值,η为游程时延;τ为游程长度,且为大于或等于1的正整数;G(τ)为游程长度为τ的游程的总数量。
CN201811124662.4A 2018-09-26 2018-09-26 一种基于游程特性的动态负荷信号分析方法 Active CN109085530B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811124662.4A CN109085530B (zh) 2018-09-26 2018-09-26 一种基于游程特性的动态负荷信号分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811124662.4A CN109085530B (zh) 2018-09-26 2018-09-26 一种基于游程特性的动态负荷信号分析方法

Publications (2)

Publication Number Publication Date
CN109085530A CN109085530A (zh) 2018-12-25
CN109085530B true CN109085530B (zh) 2021-02-02

Family

ID=64842576

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811124662.4A Active CN109085530B (zh) 2018-09-26 2018-09-26 一种基于游程特性的动态负荷信号分析方法

Country Status (1)

Country Link
CN (1) CN109085530B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113609648B (zh) * 2021-07-12 2024-04-19 湖南省计量检测研究院 基于动态负荷的电能误差计量方法、装置及电子设备
CN116008903A (zh) * 2023-01-09 2023-04-25 国网山西省电力公司太原供电公司 复杂用电工况对智能电能表运行状态影响的评价方法
CN116148750A (zh) * 2023-01-09 2023-05-23 国网山西省电力公司太原供电公司 智能电能表动态误差仿真实验系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1289480A (zh) * 1998-10-01 2001-03-28 皇家菲利浦电子有限公司 游程受限制的数字信息信号的产生
CN1801629A (zh) * 2004-12-31 2006-07-12 建兴电子科技股份有限公司 游程长度限制码错误更正解码方法
US7436331B1 (en) * 2006-04-12 2008-10-14 Marvell International Ltd. Rate-28/30 DC-free RLL code
CN104569900A (zh) * 2013-10-14 2015-04-29 北京化工大学 智能电能表动态误差测试方法与装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1289480A (zh) * 1998-10-01 2001-03-28 皇家菲利浦电子有限公司 游程受限制的数字信息信号的产生
CN1801629A (zh) * 2004-12-31 2006-07-12 建兴电子科技股份有限公司 游程长度限制码错误更正解码方法
US7436331B1 (en) * 2006-04-12 2008-10-14 Marvell International Ltd. Rate-28/30 DC-free RLL code
CN104569900A (zh) * 2013-10-14 2015-04-29 北京化工大学 智能电能表动态误差测试方法与装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
非线性负荷随机特性分析与动态测试信号建模;陈碧希;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180415;C042-727 *

Also Published As

Publication number Publication date
CN109085530A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
CN109085530B (zh) 一种基于游程特性的动态负荷信号分析方法
US11852666B2 (en) Circuit and method for width measurement of digital pulse signals
CN109255159B (zh) 一种基于机器学习的电路路径延时波动预测方法
CN109752584B (zh) 一种周期信号有效值测量方法
CN112633427A (zh) 一种基于离群点检测的超高次谐波发射信号检测方法
CN105223482B (zh) 局部放电特高频信号波形的小波分解二值去噪方法
CN106127192A (zh) 一种基于相似性的轴承剩余使用寿命预测方法
CN105699738B (zh) 一种基于pwm的交流信号有效值测量方法
CN111178587A (zh) 一种基于spark框架的短期电力负荷快速预测方法
CN105223481A (zh) 基于差值能量函数的局部放电特高频信号起始时刻确定方法
CN105242231B (zh) 数字化电能计量性能检测与溯源方法及闭环系统
CN103472333A (zh) 风电并网电能质量综合性能检测方法
CN108693084B (zh) 产生源分析装置以及产生源分析系统
CN102508210A (zh) 一种雷电脉冲信号检测及抗干扰装置
CN103969508A (zh) 一种实时高精密的电力谐波分析方法及装置
CN117748507A (zh) 基于高斯回归模型的配网谐波接入不确定性评估方法
CN110363677B (zh) 一种覆冰输电线路找形方法
CN102346809A (zh) 一种将爆破振动加速度转换为速度的方法
CN201548612U (zh) 高精度宽带峰值检波装置
CN108132399B (zh) 一种提高数字化变电站电能质量分析精度的简化插值方法
DE50005988D1 (de) Verfahren zur seismischen datenverarbeitung
CN111736014B (zh) 一种频率信号采集和数据处理方法及装置
CN107729980B (zh) 一种波形信号自适应计数方法
CN104122492A (zh) 一种预测半导体器件寿命的工作电压的方法
CN114034966A (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