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

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

【技术实现步骤摘要】

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

【技术保护点】
一种基于核函数的谐波电流概率密度计算方法，其特征在于，具体包括以下步骤：步骤1：采集用户接入系统的电流瞬时值，通过傅里叶变换得到谐波电流相量数据序列；步骤2：计算窗宽；步骤3：计算每一个分量的核函数；步骤4：利用多元随机变量的联合概率密度函数计算谐波电流概率密度。

【技术特征摘要】

【专利技术属性】
技术研发人员：陆兵能，李瑞，车静文，李小华，袁愿，
申请(专利权)人：云南电网有限责任公司曲靖供电局，厦门红相电力设备股份有限公司，
类型：发明
国别省市：云南;53