CN105403754A - 一种基于核函数的谐波电流概率密度计算方法 - Google Patents

一种基于核函数的谐波电流概率密度计算方法 Download PDF

Info

Publication number
CN105403754A
CN105403754A CN201510666441.XA CN201510666441A CN105403754A CN 105403754 A CN105403754 A CN 105403754A CN 201510666441 A CN201510666441 A CN 201510666441A CN 105403754 A CN105403754 A CN 105403754A
Authority
CN
China
Prior art keywords
harmonic current
probability density
kernel function
harmonic wave
current probability
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.)
Pending
Application number
CN201510666441.XA
Other languages
English (en)
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.)
XIAMEN RED PHASE INSTRUMENTS Inc
Qujing Power Supply Bureau Yunnan Power Grid Co Ltd
Original Assignee
XIAMEN RED PHASE INSTRUMENTS Inc
Qujing Power Supply Bureau 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 XIAMEN RED PHASE INSTRUMENTS Inc, Qujing Power Supply Bureau Yunnan Power Grid Co Ltd filed Critical XIAMEN RED PHASE INSTRUMENTS Inc
Priority to CN201510666441.XA priority Critical patent/CN105403754A/zh
Publication of CN105403754A publication Critical patent/CN105403754A/zh
Pending legal-status Critical Current

Links

Abstract

本发明公开了电力技术领域的一种基于核函数的谐波电流概率密度计算方法。其技术方案是,采集用户接入系统的电流瞬时值,通过傅里叶变换得到谐波电流相量数据序列;根据谐波电流相量数据序列计算窗宽;计算每一个分量的核函数;根据窗宽及各分量的核函数计算谐波电流概率密度。本发明的有益效果是,基于核函数的谐波电流概率密度计算方法可以比较精确地估计谐波电流概率密度,准确分析非线性负荷的谐波电流特点及其分布规律,这对进一步解决谐波污染问题、提高电能质量管理水平具有重要意义。

Description

一种基于核函数的谐波电流概率密度计算方法
技术领域
本发明属于谐波电流密度计算方法设计领域,尤其涉及一种基于核函数的谐波电流概率密度计算方法。
背景技术
随着现代工业的快速发展,电力系统不断接入大量的非线性负荷(如:电弧炉、电力机车、电动汽车充电桩等)。电力系统中不同类型的非线性负荷会呈现不同的谐波分布规律,甚至在电网中不同接入点处,相同类型的非线性负荷呈现的谐波分布规律也不一样。研究谐波的分布规律有助于理解非线性负荷的谐波特点,从而有利于提高非线性负荷的谐波管理水平。
以往的文献针对不同的非线性负荷(如:电气化铁路、新型直流输电、煤矿电网、油田电网、小型船舶电力推进系统以及分布式电源)的谐波分布计算进行了讨论。对于谐波分布,其研究方法主要分为两类,一类是从物理特性的角度出发,构建谐波的等值模型,探讨非线性负荷的谐波特点;另一类是从数学方法的角度出发,
根据非线性负荷的谐波测量数据,探索不同的计算方法,研究非线性负荷的谐波分布规律。常见的数学方法有随机理论方法、最佳平方逼近法。这两种方法所获得的概率分布模型都具有简明的解析式,因此在实际应用中具有较好的适用性。但是这些方法的共同特点是默认谐波电流相量的幅值和相位是相互独立的,在这个假设前提下将相量的幅值和相角分开考虑然后计算各自的概率密度函数,将幅值和相位各自的概率密度函数作为谐波电流的分布规律。而在实际情况中,谐波电流的幅值与相角是不独立的,谐波电流的分布规律应该是同时包含有幅值和相位的联合概率密度函数。因此在求取概率密度函数等谐波分布特征参数的计算中,传统方法所用理论是不合适的。只有将实部和虚部作为整体计算谐波数据的联合概率分布,才能在理论上保证谐波分布计算的正确性。
本文首先采用多元核密度函数估计法计算样本数据的概率密度函数,然后利用正态分布抽样法实现了概率密度函数的抽样,在此基础上得到了谐波的统计特征。与传统方法相比,多元核密度函数估计法的应用将谐波矢量数据的实部和虚部作为一个整体计算其分布规律更符合实际,有效地避免了实部与虚部相互独立为假设条件所引起的误差。而正态分布抽样法的应用,可以在保持样本期望不变的条件下,改变原始数据的抽样重心,使抽取到的样本点更多的来自于概率密度大的区域,提高了联合概率密度函数计算结果的可信度。准确的计算结果能更好的反映谐波的分布特征,而基于这些分布特征构建的谐波负载概率模型,对后续相关分析起到关键作用。
发明内容
本发明提供一种基于核函数的谐波电流概率密度计算方法,其主要目的在于克服上述背景技术中提到的传统理论方法将相量的幅值和相角分开考虑然后各自计算概率密度函数的不足。
为解决上述技术问题,本发明采用如下技术方案:
一种基于核函数的谐波电流概率密度计算方法,其特征在于,具体包括以下步骤:
步骤1:采集用户接入系统的电流瞬时值,通过傅里叶变换得到谐波电流相量数据序列;
步骤2:计算窗宽;
步骤3:计算每一个分量的核函数;
步骤4:利用多元随机变量的联合概率密度函数计算谐波电流概率密度。
步骤1中,谐波电流相量数据序列
其中,表示公共连接点处接入的用户接入系统的谐波电流相量的幅值,表示公共连接点处接入的用户接入系统的谐波电流相量的相角;n为谐波电流相量数据序列的个数。
步骤2中,计算窗宽的公式为:
其中为样本中不同样本的个数。
步骤3中,计算每一个分量的核函数,第个分量的核函数公式如下:
其中为随机变量的样本协方差,表示公共连接点处接入的用户接入系统的谐波电流相量的幅值,表示公共连接点处接入的用户接入系统的谐波电流相量的相角,n为谐波电流相量数据序列的个数。
步骤4中,可以求出谐波电流概率密度。
公式为:
其中,表示公共连接点处接入的用户接入系统的谐波电流的幅值,表示公共连接点处接入的用户接入系统的谐波电流的相角,n为谐波电流相量数据序列的个数。
本发明的有益效果是,基于核函数的谐波电流概率密度计算方法可以比较精确地估计谐波电流概率密度,这对进一步解决谐波污染问题、提高电能质量管理水平具有重要意义。
附图说明
图1是本发明提供的一种基于核函数的谐波电流概率密度计算方法的实测系统主接线图。
图2是本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流幅值有效值曲线。
图3是本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流相位有效值曲线。
图4是本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流概率密度。
图5是本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流概率密度的等高线图。
具体实施方式
下面结合附图,对优选的实施例作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。
图1是本发明提供的一种基于核函数的谐波电流概率密度计算方法的实测系统主接线图。
实测数据来自某变电站35kV用户进线侧,测量点位于35kV母线1钢厂所在支路,利用电能监测仪采集测量点处电能质量数据。
步骤1:采集测量点处电流瞬时值,通过傅里叶变换得到5次谐波电流向量数据序列:,
图2显示了本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流幅值有效值曲线。
图3显示了本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流相位有效值曲线。
步骤2:利用步骤1得到的谐波电流数据序列,计算窗宽:
其中为样本中不同样本的个数。
步骤3:利用步骤1得到的谐波电流相量数据序列,计算每一个分量的核函数,其中第l个分量的核函数计算公式如下:
其中为随机变量的样本协方差,表示测量点处的5次谐波电流幅值,表示测量点处的5次谐波电流相位,5次谐波电流相量数据序列的样本数为56866。
步骤4:计算5次谐波电流概率密度:
其中,表示第l个分量核函数,表示测量点处的5次谐波电流幅值,表示测量点处的5次谐波电流相位。
图4显示了本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流概率密度。
图5显示了本发明提供的一种基于核函数的谐波电流概率密度计算方法实测系统测量点处的5次谐波电流概率密度的等高线图。在图5中,有两个等高线峰,最左边相角100°、幅值0.05A左右的等高线峰的谐波电流概率密度最大,处于0.04以上,然后随着等高线布置,向外围扩张时依次变小,其右侧和下侧的谐波电流概率密度均降低到约0.01附近,右下侧的谐波电流概率密度降低到约0.02附近,靠右边的相角35°、幅值0.2A左右的等高线峰的谐波电流概率密度也达到最大,处于0.04以上,然后随着等高线布置,向外围扩张时依次变小,其上侧、右侧和左侧的谐波电流概率密度均降低到约0.01附近,左上方的谐波电流概率密度降低到约0.02附近。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应该涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (5)

1.一种基于核函数的谐波电流概率密度计算方法,其特征在于,具体包括以下步骤:
步骤1:采集用户接入系统的电流瞬时值,通过傅里叶变换得到谐波电流相量数据序列;
步骤2:计算窗宽;
步骤3:计算每一个分量的核函数;
步骤4:利用多元随机变量的联合概率密度函数计算谐波电流概率密度。
2.根据权利要求1所述的一种基于核函数的谐波电流概率密度计算方法,其特征在于,所述步骤1中,谐波电流相量数据序列为:
其中,表示公共连接点处接入的用户接入系统的谐波电流相量的幅值,表示公共连接点处接入的用户接入系统的谐波电流相量的相角;n为谐波电流相量数据序列的个数。
3.根据权利要求1所述的一种基于核函数的谐波电流概率密度计算方法,其特征在于,步骤2中,计算窗宽的公式为:
其中为样本中不同样本的个数。
4.根据权利要求1所述的一种基于核函数的谐波电流概率密度计算方法,其特征在于,步骤3中,计算每一个分量的核函数,第个分量的核函数公式如下:
其中为随机变量的样本协方差,表示公共连接点处接入的用户接入系统的谐波电流相量的幅值,表示公共连接点处接入的用户接入系统的谐波电流相量的相角,n为谐波电流相量数据序列的个数。
5.根据权利要求1所述的一种基于核函数的谐波电流概率密度计算方法,其特征在于,步骤4中,可以求出谐波电流概率密度,
公式为:
其中,表示公共连接点处接入的用户接入系统的谐波电流相量的幅值,表示公共连接点处接入的用户接入系统的谐波电流相量的相角,n为谐波电流相量数据序列的个数,表示窗宽。
CN201510666441.XA 2015-10-16 2015-10-16 一种基于核函数的谐波电流概率密度计算方法 Pending CN105403754A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510666441.XA CN105403754A (zh) 2015-10-16 2015-10-16 一种基于核函数的谐波电流概率密度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510666441.XA CN105403754A (zh) 2015-10-16 2015-10-16 一种基于核函数的谐波电流概率密度计算方法

Publications (1)

Publication Number Publication Date
CN105403754A true CN105403754A (zh) 2016-03-16

Family

ID=55469361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510666441.XA Pending CN105403754A (zh) 2015-10-16 2015-10-16 一种基于核函数的谐波电流概率密度计算方法

Country Status (1)

Country Link
CN (1) CN105403754A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106651150A (zh) * 2016-12-01 2017-05-10 全球能源互联网研究院 一种评估电网谐波影响的方法及装置
CN107505519A (zh) * 2017-07-26 2017-12-22 中国电力科学研究院 一种分布式电源接入电网电能质量分析方法及装置
CN114217132A (zh) * 2021-11-11 2022-03-22 江苏弈赫能源科技有限公司 一种基于wgan的电力系统谐波规律计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008064684A (ja) * 2006-09-08 2008-03-21 Nissin Electric Co Ltd 電力系統のインピーダンス測定装置
CN202018504U (zh) * 2011-01-25 2011-10-26 华北电力科学研究院有限责任公司 一种电能质量分析仪检测系统
CN103842922A (zh) * 2011-07-19 2014-06-04 智能信号公司 通过模式序列使用内核回归建模的监控系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008064684A (ja) * 2006-09-08 2008-03-21 Nissin Electric Co Ltd 電力系統のインピーダンス測定装置
CN202018504U (zh) * 2011-01-25 2011-10-26 华北电力科学研究院有限责任公司 一种电能质量分析仪检测系统
CN103842922A (zh) * 2011-07-19 2014-06-04 智能信号公司 通过模式序列使用内核回归建模的监控系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张宇辉等: "基于KDICA的间谐波分层提取方法", 《电测与仪表》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106651150A (zh) * 2016-12-01 2017-05-10 全球能源互联网研究院 一种评估电网谐波影响的方法及装置
CN106651150B (zh) * 2016-12-01 2021-02-09 全球能源互联网研究院有限公司 一种评估电网谐波影响的方法及装置
CN107505519A (zh) * 2017-07-26 2017-12-22 中国电力科学研究院 一种分布式电源接入电网电能质量分析方法及装置
CN107505519B (zh) * 2017-07-26 2021-03-02 中国电力科学研究院 一种分布式电源接入电网电能质量分析方法及装置
CN114217132A (zh) * 2021-11-11 2022-03-22 江苏弈赫能源科技有限公司 一种基于wgan的电力系统谐波规律计算方法

Similar Documents

Publication Publication Date Title
CN103163373B (zh) 适用于电能质量数据交换格式的谐波污染责任的计算方法
Wei et al. Waveform matching approach for fault diagnosis of a high-voltage transmission line employing harmony search algorithm
CN102025152B (zh) 电气化铁路接入电力系统电能质量评估方法
CN103901273B (zh) 电网谐波检测方法和电网谐波检测装置
CN105403754A (zh) 一种基于核函数的谐波电流概率密度计算方法
CN104793082A (zh) 基于谐波相关分析的电力系统谐波源识别装置
CN104008512A (zh) 一种电力系统在线稳定评估指标系统
Liu et al. Frequency-scanning harmonic generator for (inter) harmonic impedance tests and its implementation in actual 2× 25 kV railway systems
CN104320157B (zh) 一种电力线双向工频通信上行信号检测方法
Mohan et al. A weighted multi-objective index based optimal distributed generation planning in distribution system
CN104281736A (zh) 一种城市轨道交通牵引供配电系统供电电缆的选型方法
CN104392397A (zh) 一种基于熵权法的短路电流抑制方案评价方法
CN105048473B (zh) 变电站自动电压控制中低压无功设备电压灵敏度统计方法
CN105891612A (zh) 一种三相不平衡治理效果量化评估方法
CN103605856B (zh) 基于分数阶线路模型的输电线路参数估计方法
CN109870615B (zh) 适用于工矿企业的谐波责任评定方法
CN103560497B (zh) 一种基于电网拓扑调整的短路电流限制方法
Wang et al. Novel travelling wave fault location principle for VSC-HVDC transmission line
CN105140927A (zh) 一种含谐波的配电网无功补偿电容器容量的计算方法
Liao et al. A fault location method based on traveling wave natural frequency used on±800kV UHVDC transmission lines
CN102185314A (zh) 电能质量综合评判方法
CN109599885A (zh) 一种特高压直流闭锁换流母线暂态过电压计算方法
Gupta et al. Power system network equivalents: Key issues and challenges
CN104732113B (zh) 一种基于谐波源注入的并联型apf装机容量预估方法
CN103384070A (zh) 一种基于psasp的电能质量分析计算方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160316