CN111625773B - 基于可信性测度的计量回路电流测量异常评估方法及装置 - Google Patents

基于可信性测度的计量回路电流测量异常评估方法及装置 Download PDF

Info

Publication number
CN111625773B
CN111625773B CN202010448059.2A CN202010448059A CN111625773B CN 111625773 B CN111625773 B CN 111625773B CN 202010448059 A CN202010448059 A CN 202010448059A CN 111625773 B CN111625773 B CN 111625773B
Authority
CN
China
Prior art keywords
current
current amplitude
user group
calculating
dynamic load
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
CN202010448059.2A
Other languages
English (en)
Other versions
CN111625773A (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power 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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010448059.2A priority Critical patent/CN111625773B/zh
Publication of CN111625773A publication Critical patent/CN111625773A/zh
Application granted granted Critical
Publication of CN111625773B publication Critical patent/CN111625773B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/02Testing or calibrating of apparatus covered by the other groups of this subclass of auxiliary devices, e.g. of instrument transformers according to prescribed transformation ratio, phase angle, or wattage rating
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明提供了基于可信性测度的计量回路电流测量异常评估方法及装置,包括以下步骤:步骤一:收集平均电流幅值数据,计算偏移特征量;步骤二:计算不同电流幅值偏移特征量的电流无异常概率;步骤三:将不同电流幅值偏移特征量与对应电流无异常的概率拟合得到电流隶属度函数,计算可信性测度;步骤四:设定用户群电流幅值测量预警值,并计算动态负荷电流幅值越限严重度;步骤五:应用可信性测度与动态负荷电流幅值越限严重度判断台区是否异常,本发明运用模糊数学理论,能在复杂场景下评估出台区电流测量异常的可能性,方法灵活,限制少;采用改进的可信性测度方法,能更好将电流幅值无异常概率与异常概率相结合,提供更精确的计算结果。

Description

基于可信性测度的计量回路电流测量异常评估方法及装置
技术领域
本发明涉及台区计量回路电流测量技术,具体是基于可信性测度的计量回路电流测量异常评估方法及装置。
背景技术
随着市场经济的发展,电力事业飞速发展的同时,社会对电能的依赖性越来越强。因此,保障电力系统安全、稳定的运行有着十分重要的意义。但这些年开始逐步推进电力市场化改革,电能计量部门面临一系列新的计量检测和电流电压异常评估问题,如何解决这些问题正在成为电能计量工作的当务之急。
目前,由于电力电子技术的广泛应用,电力系统中谐波源大量增加,引起局部电流波形的畸变,若其电流波形畸变与人为窃电造成电流波形畸变同时产生,会影响传统的自动计量检测系统对窃电电流异常评估的误判。若用人工对电流测量异常评估又过分依赖抄表人员和稽查人员的自身素质,且大部分情况下还是依靠经验评估电流是否异常。由于人工方法信息化程度低,不能及时对异常状况做出反应,因此采用人工评估存在较大的隐患。近些年,模糊数学作为较好的综合评估方法得到高速发展。依据模糊数学而提出的可信性测度理论,解决了模糊数学中混乱和不兼容等问题。
发明内容
本发明要解决的技术问题:针对现有技术的上述问题,提供一种利用模糊数学中的可信性测度进行计量回路电流测量异常评估方法及装置,刻画电流测量中的多重不确定性,以此评估电流测量异常发生的可能性。
为实现上述目的,本发明实施方式的一方面提供一种基于可信性测度的计量回路电流测量异常评估方法,包括以下步骤:
步骤一:计算被测台区N个不同用户群近一个月的平均电流幅值{i1、i2、i3……iN},将其平均值I作为该台区的理想电流幅值,并将平均电流幅值{i1、i2、i3……iN}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……xN};
步骤二:将被测台区不少于2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征值建立历史电流幅值偏移数据库,查询该台区电流幅值偏移特征量xn(n=1,2,3……N)在正常偏移范围内出现的偏移量总次数dsn(n=1,2,3……N),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……N),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……N)下电流无异常的概率Yn(n=1,2,3……N);
步骤三:将不同电流幅值偏移特征量xn(n=1,2,3……N)与对应电流无异常的概率Yn(n=1,2,3……N)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1,记为台区电流隶属度函数μ(x),根据该台区电流隶属度函数μ(x)计算得到电流幅值偏移正常量z的范围,由此计算第n个用户群的可信性测度Mcr(n);
步骤四:通过专家投票法设定该台区电流幅值测量预警值ap,即预警值ap由T个专家(T大于等于5)提出该台区电流幅值测量预警apt(t=1,2,3……T)求平均值得到,用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群动态负荷电流幅值越限严重度Sev(n),得到电网受扰程度;
步骤五:应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n),若出现某用户群评估结果出现异常,则判断该台区计量回路电流测量处于异常状态。
进一步的,所述步骤三中通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure BDA0002506667110000031
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量;λ为集合{x1、x2、x3……xN}中的元素;
Figure BDA0002506667110000032
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure BDA0002506667110000033
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高。
进一步的,所述步骤四中利用不同用户群的平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算不同用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure BDA0002506667110000034
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值。
进一步的,所述步骤五中第n个用户群动态负荷电流测量的异常评估结果计算公式为:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)]
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
为实现上述目的,本发明实施方式的另一方面还提供一种基于可信性测度的计量回路电流测量异常评估装置,包括:用于收集被测台区N个不同用户群近一个月的平均电流幅值{i1、i2、i3……iN}的数据收集模块,用于计算平均电流幅值{i1、i2、i3……iN}的平均值I作为该台区的理想电流幅值的理想电流计算模块,用于将平均电流幅值{i1、i2、i3……iN}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……xN}的偏移特征量计算模块,用于将被测台区不少于2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征值建立历史电流幅值偏移数据库,查询该台区电流幅值偏移特征量xn(n=1,2,3……N)在正常偏移范围内出现的偏移量总次数dsn(n=1,2,3……N),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……N),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……N)下电流无异常的概率Yn(n=1,2,3……N)的无异常概率计算模块,用于将不同电流幅值偏移特征量xn(n=1,2,3……N)与对应电流无异常的概率Yn(n=1,2,3……N)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1,记为台区电流隶属度函数μ(x)的隶属度函数计算模块,用于根据台区电流隶属度函数μ(x)计算得到电流幅值偏移正常量z的范围,由此计算第n个用户群的可信性测度Mcr(n)的可信性测度计算模块,用于通过专家投票法设定该台区电流幅值测量预警值ap,用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群动态负荷电流幅值越限严重度Sev(n)的越限严重度计算模块,用于应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n)的评估检测模块。
进一步的,通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure BDA0002506667110000051
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量;λ为集合{x1、x2、x3……xN}的元素;
Figure BDA0002506667110000052
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure BDA0002506667110000053
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高。
进一步的,利用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure BDA0002506667110000054
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值。
进一步的,第n个用户群动态负荷电流测量的异常评估结果计算公式为:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)]
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
与现有技术相比,本发明的有益效果是:
1.运用模糊数学理论,能在复杂场景下评估出台区电流测量异常的可能性,方法灵活,限制少。
2.采用改进的可信性测度方法,能更好将电流幅值无异常概率与异常概率相结合,提供更精确的计算结果。
3.与传统的电流测量异常判断不同,本发明方法计算的是电流测量异常的可能性,考虑了不确定性,在计算方法上有非常好的鲁棒性。
附图说明
图1为本发明实现基于可信性测度的计量回路电流测量异常评估方法的流程图;
图2为电流幅值偏移特征量-无异常概率图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步的详细说明。此处所描述的具体实施例仅用于解释本发明技术方案,并不限于本发明。
如图1所示,一种基于可信性测度的计量回路电流测量异常评估方法流程如下:
第一步,假设被测台区共有5(N=5)个不同的用户群,通过平均电流幅值的数据收集模块读取被测台区的用户群编号、近一个月的电流的平均幅值数据{i1、i2、i3……i5},如表1所示。
运用理想电流计算模块,将获得的5个用户群近一个月的电流的平均幅值数据{i1、i2、i3……i5}求平均值得到该台区的理想电流幅值I:
I=(80.1+70.2+89.1+72.3+60.7)/5=74.48A (1)
通过偏移特征量计算模块将不同用户群的平均电流幅值{i1、i2、i3……i5}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……x5},如表2所示。
表1不同用户群的平均电流幅值
Figure BDA0002506667110000061
表2不同用户群的电流幅值偏移特征量
Figure BDA0002506667110000071
第二步,将被测台区近2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征量建立历史电流幅值偏移数据库,如表3所示。
表3历史电流幅值偏移特征数据库
Figure BDA0002506667110000072
当x1=5.62时,通过无异常概率计算模块查询历史电流幅值偏移特征数据库中在正常偏移范围[5.339,5.901]内出现的偏移量总次数为5,并查询与之对应的没有出现台区电流异常的次数为4,得到电流幅值偏移特征量x1下电流无异常概率为80%。以此类推,查询该台区电流幅值偏移特征量xn在正常偏移范围[0.95xn,1.05xn](n=1,2,3……5)内出现的偏移量总次数dsn(n=1,2,3……5),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……5),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……5)下电流无异常的概率Yn(n=1,2,3……5);将无异常概率与电流幅值偏移特征量对应,得到图2。
第三步,通过隶属度函数计算模块将不同电流幅值偏移特征量xn(n=1,2,3……5)与对应电流无异常的概率Yn(n=1,2,3……5)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1得到的电流隶属度函数μ(x)为:
μ(x)=-0.0047x2+0.0041x+0.98 (2)
根据台区电流隶属度函数μ(x),计算得到电流无异常概率在95%以上时的电流幅值偏移特征量称为电流幅值偏移正常量z,其范围为[-2.54,2.78]。
通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure BDA0002506667110000081
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量,其范围为[-2.54,2.78];λ为集合{x1、x2、x3……x5}的元素;
Figure BDA0002506667110000082
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure BDA0002506667110000083
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高。
计算各用户群的可信性测度Mcr(n),如表4所示。
表4不同用户群的可信性测度表
Figure BDA0002506667110000084
第四步,通过专家投票法设定预警值ap,5(T=5)个专家对电流幅值测量预警值apt(t=1,2,3……5)投票结果如表5所示。
表5电流预警值aprof投票结果
Figure BDA0002506667110000085
Figure BDA0002506667110000091
计算平均值可得整个台区电流幅值测量预警值ap为:
(70+75*2+80*2)/5=76A (4)
利用不同用户群的平均电流幅值{i1、i2、i……i5}与预警值ap分别计算不同用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure BDA0002506667110000092
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值87A。
得到5个用户群动态负荷电流幅值越限严重度Sev(n),如表6所示。
表6不同用户群的动态负荷电流幅值越限严重度
Figure BDA0002506667110000093
第五步,应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n),公式如下:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)] (6)
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
得到该台区5个用户群电流异常评估结果,如表7所示。
表7不同用户群的电流异常评估结果
Figure BDA0002506667110000101
将电流发生异常的可能性分划为G个等级(G大于等于5),并设立每个指标在相应的评语集下的评判标准。通过调查问卷或者邀请多名专家将不同的Mrisk数值范围与G个等级进行对应,得到计量回路电流测量异常评估体系,具体标准举例如表8所示。
表8评判标准
Figure BDA0002506667110000102
最后结合表8评判标准可得各用户群的电流异常评估等级,如表9所示。
表9不同用户群的电流异常评估等级
Figure BDA0002506667110000103
由此可以看出该台区下的用户群3处于电流测量异常状态,故判断整个台区计量回路电流测量异常。
本发明实施方式的另一方面还提供一种基于可信性测度的计量回路电流测量异常评估装置,包括:用于收集被测台区N个不同用户群近一个月的平均电流幅值{i1、i2、i3……iN}的数据收集模块,用于计算平均电流幅值{i1、i2、i3……iN}的平均值I作为该台区的理想电流幅值的理想电流计算模块,用于将平均电流幅值{i1、i2、i3……iN}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……xN}的偏移特征量计算模块,用于将被测台区不少于2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征值建立历史电流幅值偏移数据库,查询该台区电流幅值偏移特征量xn(n=1,2,3……N)在正常偏移范围内出现的偏移量总次数dsn(n=1,2,3……N),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……N),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……N)下电流无异常的概率Yn(n=1,2,3……N)的无异常概率计算模块,用于将不同电流幅值偏移特征量xn(n=1,2,3……N)与对应电流无异常的概率Yn(n=1,2,3……N)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1,记为台区电流隶属度函数μ(x)的隶属度函数计算模块,用于根据台区电流隶属度函数μ(x)计算得到电流幅值偏移正常量z的范围,由此计算第n个用户群的可信性测度Mcr(n)的可信性测度计算模块,用于通过专家投票法设定该台区电流幅值测量预警值ap,用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群动态负荷电流幅值越限严重度Sev(n)的越限严重度计算模块,用于应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n)的评估检测模块。
进一步的,通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure BDA0002506667110000111
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量;λ为集合{x1、x2、x3……xN}的元素;
Figure BDA0002506667110000112
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure BDA0002506667110000113
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高。
进一步的,利用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure BDA0002506667110000121
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值。
进一步的,第n个用户群动态负荷电流测量的异常评估结果计算公式为:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)]
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
以上所述仅表达了本发明的优选实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形、改进及替代,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (2)

1.基于可信性测度的计量回路电流测量异常评估方法,其特征在于:包括以下步骤:
步骤一:计算被测台区N个不同用户群近一个月的平均电流幅值{i1、i2、i3……iN},将其平均值I作为该台区的理想电流幅值,并将平均电流幅值{i1、i2、i3……iN}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……xN};
步骤二:将被测台区不少于2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征值建立历史电流幅值偏移数据库,查询该台区电流幅值偏移特征量xn(n=1,2,3……N)在正常偏移范围内出现的偏移量总次数dsn(n=1,2,3……N),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……N),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……N)下电流无异常的概率Yn(n=1,2,3……N);
步骤三:将不同电流幅值偏移特征量xn(n=1,2,3……N)与对应电流无异常的概率Yn(n=1,2,3……N)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1,记为台区电流隶属度函数μ(x),根据该台区电流隶属度函数μ(x)计算得到电流幅值偏移正常量z的范围,由此计算第n个用户群的可信性测度Mcr(n);
步骤四:通过专家投票法设定该台区电流幅值测量预警值ap,用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群动态负荷电流幅值越限严重度Sev(n),得到电网受扰程度;
步骤五:应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n),若出现某用户群评估结果出现异常,则判断该台区计量回路电流测量处于异常状态;
所述步骤三中通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure FDA0004142818450000021
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量;λ为集合{x1、x2、x3……xN}中的元素;
Figure FDA0004142818450000022
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure FDA0004142818450000023
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高;
所述步骤四中利用不同用户群的平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算不同用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure FDA0004142818450000024
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值;
所述步骤五中第n个用户群动态负荷电流测量的异常评估结果计算公式为:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)]
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
2.基于可信性测度的计量回路电流测量异常评估装置,其特征在于:包括:用于收集被测台区N个不同用户群近一个月的平均电流幅值{i1、i2、i3……iN}的数据收集模块,用于计算平均电流幅值{i1、i2、i3……iN}的平均值I作为该台区的理想电流幅值的理想电流计算模块,用于将平均电流幅值{i1、i2、i3……iN}与理想电流幅值I分别相减得到电流幅值偏移特征量{x1、x2、x3……xN}的偏移特征量计算模块,用于将被测台区不少于2年的历史月度平均电流幅值与理想电流幅值I相减得到电流幅值偏移特征值建立历史电流幅值偏移数据库,查询该台区电流幅值偏移特征量xn(n=1,2,3……N)在正常偏移范围内出现的偏移量总次数dsn(n=1,2,3……N),并查询与之对应的没有出现台区电流异常的次数den(n=1,2,3……N),将没有出现台区电流异常的次数den除以总次数dsn得到整个台区在不同电流幅值偏移特征量xn(n=1,2,3……N)下电流无异常的概率Yn(n=1,2,3……N)的无异常概率计算模块,用于将不同电流幅值偏移特征量xn(n=1,2,3……N)与对应电流无异常的概率Yn(n=1,2,3……N)通过最小二乘法拟合得到多项式函数Y(x)=k1xp+k2xp-1+...+kpx+kp+1,记为台区电流隶属度函数μ(x)的隶属度函数计算模块,用于根据台区电流隶属度函数μ(x)计算得到电流幅值偏移正常量z的范围,由此计算第n个用户群的可信性测度Mcr(n)的可信性测度计算模块,用于通过专家投票法设定该台区电流幅值测量预警值ap,用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群动态负荷电流幅值越限严重度Sev(n)的越限严重度计算模块,用于应用可信性测度Mcr(n)和动态负荷电流幅值越限严重度Sev(n)计算得到各用户群动态负荷电流测量的异常评估结果Mrisk(n)的评估检测模块;
通过台区电流隶属度函数μ(x)求得各用户群的可信性测度Mcr(n)公式如下所示:
Figure FDA0004142818450000041
其中,xn为第n个用户群电流幅值偏移特征量,z为电流幅值偏移正常量;λ为集合{x1、x2、x3……xN}的元素;
Figure FDA0004142818450000042
表示所有λ小于等于xn时μ(λ)的上确界;/>
Figure FDA0004142818450000043
表示所有λ大于等于xn时μ(λ)的下确界;Mcr(n)越接近1表示电流无异常的概率越高;
利用平均电流幅值{i1、i2、i3……iN}与预警值ap分别计算第n个用户群的动态负荷电流幅值越限严重度Sev(n),公式如下:
Figure FDA0004142818450000044
其中,in为第n个用户群的平均电流幅值;amax为当前台区电流幅值测量历史结果的最高值;
第n个用户群动态负荷电流测量的异常评估结果计算公式为:
Mrisk(n)=tanh[(1-Mcr(n))*Sev(n)]
其中,Mrisk(n)为第n个用户群动态负荷电流测量的异常评估结果,其取值范围为[0,1],越接近1则该用户群电流测量异常的可能性越大;tanh表示双曲正切函数。
CN202010448059.2A 2020-05-25 2020-05-25 基于可信性测度的计量回路电流测量异常评估方法及装置 Active CN111625773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010448059.2A CN111625773B (zh) 2020-05-25 2020-05-25 基于可信性测度的计量回路电流测量异常评估方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010448059.2A CN111625773B (zh) 2020-05-25 2020-05-25 基于可信性测度的计量回路电流测量异常评估方法及装置

Publications (2)

Publication Number Publication Date
CN111625773A CN111625773A (zh) 2020-09-04
CN111625773B true CN111625773B (zh) 2023-06-16

Family

ID=72259106

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010448059.2A Active CN111625773B (zh) 2020-05-25 2020-05-25 基于可信性测度的计量回路电流测量异常评估方法及装置

Country Status (1)

Country Link
CN (1) CN111625773B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113435490B (zh) * 2021-06-18 2022-06-28 国网宁夏电力有限公司电力科学研究院 一种基于动态规划时间序列相似性算法的谐波溯源方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323469A (zh) * 2011-07-27 2012-01-18 四川大学 谐波负荷状态监测系统
US8271421B1 (en) * 2007-11-30 2012-09-18 Intellectual Assets Llc Nonparametric fuzzy inference system and method
CN102915515A (zh) * 2012-11-10 2013-02-06 四川省电力公司资阳公司 一种电网连锁性故障识别和风险评估方法
CN105976257A (zh) * 2015-12-17 2016-09-28 国家电网公司 基于隶属度函数的模糊综合评价法的电网脆弱性评估方法
CN109711663A (zh) * 2018-11-15 2019-05-03 国网山东省电力公司淄博供电公司 基于大数据分析的变电站油浸式变压器状态评估与修正方法及系统
CN110648088A (zh) * 2019-11-26 2020-01-03 国网江西省电力有限公司电力科学研究院 一种基于鸟群算法与svm的电能质量扰动源判断方法
CN110991826A (zh) * 2019-11-18 2020-04-10 国网浙江省电力有限公司电力科学研究院 一种评价低压电能表运行状态的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105205224B (zh) * 2015-08-28 2018-10-30 江南大学 基于模糊曲线分析的时间差高斯过程回归软测量建模方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8271421B1 (en) * 2007-11-30 2012-09-18 Intellectual Assets Llc Nonparametric fuzzy inference system and method
CN102323469A (zh) * 2011-07-27 2012-01-18 四川大学 谐波负荷状态监测系统
CN102915515A (zh) * 2012-11-10 2013-02-06 四川省电力公司资阳公司 一种电网连锁性故障识别和风险评估方法
CN105976257A (zh) * 2015-12-17 2016-09-28 国家电网公司 基于隶属度函数的模糊综合评价法的电网脆弱性评估方法
CN109711663A (zh) * 2018-11-15 2019-05-03 国网山东省电力公司淄博供电公司 基于大数据分析的变电站油浸式变压器状态评估与修正方法及系统
CN110991826A (zh) * 2019-11-18 2020-04-10 国网浙江省电力有限公司电力科学研究院 一种评价低压电能表运行状态的方法
CN110648088A (zh) * 2019-11-26 2020-01-03 国网江西省电力有限公司电力科学研究院 一种基于鸟群算法与svm的电能质量扰动源判断方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
辨识谐波电流监测数据中异常数据的一种方法研究;马智远等;《电力系统保护与控制》;20161101(第21期);全文 *

Also Published As

Publication number Publication date
CN111625773A (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
CN114640173B (zh) 一种基于多特征量的变压器和发电机的预警模型
Bastidas-Rodríguez et al. Model-based degradation analysis of photovoltaic modules through series resistance estimation
CN107346899B (zh) 一种光伏电站系统稳定性评估方法及系统
CN104819107B (zh) 一种风电机组功率曲线异常漂移的诊断方法及系统
CN112381476B (zh) 用于确定状态异常的电能表的方法及装置
CN111275295B (zh) 一种基于反距离权重插值的分布式光伏故障诊断方法
CN103886518A (zh) 一种基于监测点电能质量数据挖掘的电压暂降预警方法
CN113783272A (zh) 一种基于超级电容监控管理系统的安全控制方法
CN110738346A (zh) 一种基于威布尔分布的批次电能表可靠度预测方法
CN111625773B (zh) 基于可信性测度的计量回路电流测量异常评估方法及装置
CN118100446B (zh) 一种用于电力系统的节点运行状态智能检测方法
CN113985339A (zh) 智能电表的误差诊断方法及系统、设备、存储介质
CN117277435A (zh) 一种光伏逆变器的健康评估方法、系统及装置
Sarkar et al. IoT based wireless sensor network (WSN) for condition monitoring of low power rooftop PV panels
Aboshady et al. Fault detection and classification scheme for PV system using array power and cross-strings differential currents
CN108197843B (zh) 一种平坦地形风电出力特性评估方法
CN114063003A (zh) 基于小台区的电能表测量误差检测方法、系统及存储介质
Benninger et al. Anomaly detection by comparing photovoltaic systems with machine learning methods
CN112147409B (zh) 低压配电线路在线阻抗检测方法及装置
CN114039547A (zh) 一种光伏组件隐裂故障的诊断方法
CN112180194A (zh) 一种配电网多源数据质量控制方法
CN118501791B (zh) 一种计量检测数据精度评估方法及系统
Roger et al. Fault diagnostic methodology for grid-connected photovoltaic systems
CN118157325B (zh) 一种用于新能源电力的实时监测方法及系统
CN118336926B (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
TA01 Transfer of patent application right

Effective date of registration: 20200929

Address after: 330096 Jiangxi City, private science and Technology Park in the Castle Peak Park, private road, No. 88, test center, science and Technology Research Center (level 1-11)

Applicant after: STATE GRID JIANGXI ELECTRIC POWER COMPANY LIMITED Research Institute

Applicant after: STATE GRID CORPORATION OF CHINA

Address before: 330096 Jiangxi City, private science and Technology Park in the Castle Peak Park, private road, No. 88, test center, science and Technology Research Center (level 1-11)

Applicant before: STATE GRID JIANGXI ELECTRIC POWER COMPANY LIMITED Research Institute

Applicant before: STATE GRID CORPORATION OF CHINA

Applicant before: NANCHANG KECHEN ELECTRIC POWER TEST AND RESEARCH Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant