CN110045186B - 电力系统公共连接点处的系统谐波阻抗估计方法和系统 - Google Patents

电力系统公共连接点处的系统谐波阻抗估计方法和系统 Download PDF

Info

Publication number
CN110045186B
CN110045186B CN201910341895.8A CN201910341895A CN110045186B CN 110045186 B CN110045186 B CN 110045186B CN 201910341895 A CN201910341895 A CN 201910341895A CN 110045186 B CN110045186 B CN 110045186B
Authority
CN
China
Prior art keywords
harmonic
harmonic voltage
system side
value
fluctuation
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
CN201910341895.8A
Other languages
English (en)
Other versions
CN110045186A (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201910341895.8A priority Critical patent/CN110045186B/zh
Publication of CN110045186A publication Critical patent/CN110045186A/zh
Application granted granted Critical
Publication of CN110045186B publication Critical patent/CN110045186B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供了一种电力系统公共连接点处的系统谐波阻抗估计方法和系统，其中方法包括：建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据；计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值；对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，基于系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。本发明提供的方法和系统能够准确而快速得到公共连接点处的系统侧谐波阻抗估计值。

Description

电力系统公共连接点处的系统谐波阻抗估计方法和系统
技术领域
本发明涉及电力系统技术领域，尤其是涉及一种电力系统公共连接点处的系统谐波阻抗估计方法和系统。
背景技术
近年来，人们对电能的质量提出越来越高的要求。然而电网中大量电力电子装置和非线性负荷的接入，导致谐波问题愈发严重，不仅对电力系统的稳定运行造成干扰，还给各方带来巨大的经济损失。为监测、治理谐波问题，国内外均出台相关政策，依据谐波发射水平进行相应奖惩，合理量化谐波污染责任至关重要，而评估谐波发射水平的关键在于系统侧谐波阻抗的准确计算。然而，现有对系统侧谐波阻抗进行计算的方法包括“干预式”和“非干预式”，都存在计算结果不准确和计算复杂的技术问题。
发明内容
有鉴于此，本发明的目的在于提供一种电力系统公共连接点处的系统谐波阻抗估计方法和系统，本发明能够准确而快速地得到电力系统系统侧谐波阻抗估计值。
第一方面，本发明实施例提供了一种电力系统公共连接点处的系统谐波阻抗估计方法，包括：建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；所述戴维南等效电路模型用于表征所述系统侧电路的谐波电压源和系统侧谐波阻抗之间的连接关系；在所述电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组；基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值；基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；所述第一系统侧谐波电压源波动量函数用于描述所述电力系统公共连接点处的系统侧谐波电压源波动量的大小；将所述电学关系方程代入到所述第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；所述第二系统侧谐波电压源波动量函数用于描述所述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系；基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值。
进一步地，基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值包括：计算所述第二系统侧谐波电压源波动量函数的最小值；确定所述第二系统侧谐波电压源波动量函数的数值为所述最小值时，所述系统侧谐波阻抗的取值，并将所述系统侧谐波阻抗的取值作为所述电力系统公共连接点处的系统侧谐波阻抗估计值。
进一步地，基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值包括：通过以下公式计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值：
Figure BDA0002040750800000021
其中，
Figure BDA0002040750800000022
为第n个谐波电压波动值，
Figure BDA0002040750800000023
为所述谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure BDA0002040750800000031
为所述谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数；
Figure BDA0002040750800000032
其中，
Figure BDA0002040750800000033
为第n个谐波电流波动值，
Figure BDA0002040750800000034
为所述谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure BDA0002040750800000035
为所述谐波电流观测数据组中第n个采样数据的谐波电流值。
进一步地，基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程包括：基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程为
Figure BDA0002040750800000036
对所述戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure BDA0002040750800000037
其中，Zu为系统侧谐波阻抗，
Figure BDA0002040750800000038
为系统侧谐波电压源波动量，
Figure BDA0002040750800000039
为系统侧谐波电压源电压值。
进一步地，根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数包括：根据2范数构建数学表达式为
Figure BDA00020407508000000310
的第一系统侧谐波电压源波动量函数，其中，N为所述多个谐波电压数据或所述多个谐波电流数据的采集个数。
第二方面，本发明实施例还提供了一种电力系统公共连接点处的系统谐波阻抗估计系统，包括：等效电路模块，数据采集模块，数据处理模块，第一计算模块，第二计算模块和第三计算模块，其中，所述等效电路模块，用于建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；所述戴维南等效电路模型用于表征所述系统侧电路的谐波电压源和系统侧谐波阻抗之间的连接关系；所述数据采集模块，用于在所述电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组；所述数据处理模块，用于基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值；所述第一计算模块，用于基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；所述第二计算模块，用于根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；所述第一系统侧谐波电压源波动量函数用于描述所述电力系统公共连接点处的系统侧谐波电压源波动量的大小；以及将所述电学关系方程代入到所述第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；所述第二系统侧谐波电压源波动量函数用于描述所述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系；所述第三计算模块，用于基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值。
进一步地，所述第三计算模块还用于：计算所述第二系统侧谐波电压源波动量函数的最小值；确定所述第二系统侧谐波电压源波动量函数的数值为所述最小值时，所述系统侧谐波阻抗的取值，并将所述系统侧谐波阻抗的取值作为所述电力系统公共连接点处的系统侧谐波阻抗估计值。
进一步地，所述数据处理模块包括：谐波电压处理单元和谐波电流处理单元，其中，所述谐波电压处理单元，用于通过公式
Figure BDA0002040750800000041
Figure BDA0002040750800000042
计算所述电力系统公共连接点处的谐波电压波动值；其中，
Figure BDA0002040750800000043
为第n个谐波电压波动值，
Figure BDA0002040750800000044
为所述谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure BDA0002040750800000045
为所述谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数；所述谐波电流处理单元，用于通过公式
Figure BDA0002040750800000051
计算所述电力系统公共连接点处的谐波电流波动值；其中，
Figure BDA0002040750800000052
为第n个谐波电流波动值，
Figure BDA0002040750800000053
为所述谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure BDA0002040750800000054
为所述谐波电流观测数据组中第n个采样数据的谐波电流值。
进一步地，所述第一计算模块还用于：基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程为
Figure BDA0002040750800000055
对所述戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure BDA0002040750800000056
其中，Zu为系统侧谐波阻抗，
Figure BDA0002040750800000057
为系统侧谐波电压源波动量，
Figure BDA0002040750800000058
为系统侧谐波电压源电压值。
进一步地，所述第二计算模块还用于：根据2范数构建数学表达式为
Figure BDA0002040750800000059
的第一系统侧谐波电压源波动量函数，其中，N为所述多个谐波电压数据或所述多个谐波电流数据的采集个数。
本发明实施例带来了以下有益效果：本发明实施例通过在电力系统中，建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据；然后计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值；对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，最后基于系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。本发明通过对系统侧谐波电压源波动量的估计来计算系统侧谐波阻抗，这种方式能够使得到的系统侧谐波阻抗的估计值更为准确；同时通过以上方式能够以较少的计算步骤快速地得到电力系统的系统侧谐波阻抗估计值。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种电力系统公共连接点处的系统谐波阻抗估计方法的流程图；
图2为本发明实施例提供的一种戴维南等效电路模型的示意图；
图3为本发明实施例提供的一种电力系统公共连接点处的系统谐波阻抗估计系统的示意图；
图4为本发明实施例提供的另一种电力系统公共连接点处的系统谐波阻抗估计系统的示意图；
图5a为本发明实施例提供的一种应用实施例中谐波阻抗的幅值计算结果示意图；
图5b为本发明实施例提供的一种应用实施例中谐波阻抗的相角计算结果示意图；
图6a为本发明实施例提供的另一种应用实施例中谐波阻抗的幅值计算结果示意图；
图6b为本发明实施例提供的另一种应用实施例中谐波阻抗的相角计算结果示意图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一：
图1是根据本发明实施例的一种电力系统公共连接点处的系统谐波阻抗估计方法的流程图。如图1所示，该方法包括如下步骤：
步骤S102，建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；其中，戴维南等效电路模型用于表征电力系统的谐波电压源和系统侧谐波阻抗之间的连接关系。
可选地，图2为本发明实施例提供的一种戴维南等效电路模型，即，将公共连接点(Point of Common Coupling，简称PCC)系统侧复杂的谐波网络等效为一个谐波电压源
Figure BDA0002040750800000071
和谐波阻抗(Zu)的串联。
步骤S104，在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组；
具体地，假定一个观测周期采集的点数为N，N为大于1的整数，则采集的谐波电压和谐波电流，分别构成了N点的谐波电压观测数据组和N点的谐波电流观测数据组。
可选地，在本发明实施例中，PCC处采集的谐波电压和谐波电流数据分别记为
Figure BDA0002040750800000072
和
Figure BDA0002040750800000073
其中，n表示采样点序号，为大于零的正整数。
步骤S106，基于谐波电压观测数据组和谐波电流观测数据组，计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值。
步骤S108，基于戴维南等效电路模型得到系统侧电路的戴维南等效电路方程，并对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程。
步骤S110，根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；第一系统侧谐波电压源波动量函数用于描述电力系统公共连接点处的系统侧谐波电压源波动量的大小。
步骤S112，将电学关系方程代入到第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；第二系统侧谐波电压源波动量函数用于描述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系。
步骤S114，基于第二系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。
本发明实施例通过在电力系统中，建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据；然后计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值；对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，最后基于系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。本发明通过对系统侧谐波电压源波动量的估计来计算系统侧谐波阻抗，这种方式能够使得到的系统侧谐波阻抗的估计值更为准确；同时通过以上方式能够以较少的计算步骤快速地得到电力系统的系统侧谐波阻抗估计值。
具体地，步骤S106中，通过如下公式计算电力系统的谐波电压波动值和谐波电流波动值：
(1)谐波电压波动值计算：
Figure BDA0002040750800000081
其中，
Figure BDA0002040750800000082
为第n个谐波电压波动值，
Figure BDA0002040750800000083
为谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure BDA0002040750800000084
为谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数。
(2)谐波电流波动值计算：
Figure BDA0002040750800000091
其中，
Figure BDA0002040750800000092
为第n个谐波电流波动值，
Figure BDA0002040750800000093
为谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure BDA0002040750800000094
为谐波电流观测数据组中第n个采样数据的谐波电流值。
具体地，步骤S108包括如下步骤：
如图2所示，基于戴维南等效电路模型得到系统侧电路的戴维南等效电路方程为
Figure BDA0002040750800000095
对戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure BDA0002040750800000096
其中，Zu为系统侧谐波阻抗，
Figure BDA0002040750800000097
为系统侧谐波电压源波动量，
Figure BDA0002040750800000098
为系统侧谐波电压源电压值。
具体地，步骤S110中，根据2范数构建数学表达式为
Figure BDA0002040750800000099
的第一系统侧谐波电压源波动量函数，其中，N为多个谐波电压数据或多个谐波电流数据的采集个数。
具体地，步骤S112的数学形式如下：
Figure BDA0002040750800000101
其中，*为共轭向量。
具体地，步骤S114还包括如下步骤：
步骤S1141，计算第二系统侧谐波电压源波动量函数的最小值。
步骤S1142，确定第二系统侧谐波电压源波动量函数的数值为最小值时，系统侧谐波阻抗的取值，并将系统侧谐波阻抗的取值作为电力系统公共连接点处的系统侧谐波阻抗估计值。
具体地，由步骤S112得到的第二系统侧谐波电压源波动量函数的数学形式可知，当
Figure BDA0002040750800000102
时，J有最小值：
Figure BDA0002040750800000103
即，当第二系统侧谐波电压源波动量函数J的数值为最小值时，系统侧谐波阻抗
Figure BDA0002040750800000111
将从步骤S106中计算得到的谐波电压波动值和谐波电流波动值代入上式可以得到此时系统侧谐波阻抗的取值，并将此时系统侧谐波阻抗的取值作为系统侧谐波阻抗的估计值。
本发明与现有技术相比，具有如下的优点和有益效果：
(1)本发明求取使
Figure BDA0002040750800000112
最小时的Zu的值，因为
Figure BDA0002040750800000113
且Zu在一段时间内变化很小，所以认为J取最小值时，
Figure BDA0002040750800000114
取最小值，因此Zu是在背景谐波最稳定时计算得到的。故该发明在背景谐波剧烈波动时，利用背景谐波最稳定时的数据计算得到Zu，求得的Zu也最为准确。
(2)本发明计算步骤少且计算简单，易于操作实现，计算速度更快。
实施例二：
图3是根据本发明实施例提供的一种电力系统公共连接点处的系统谐波阻抗估计系统的示意图，如图3所示，该系统具体包括：等效电路模块10，数据采集模块20，数据处理模块30，第一计算模块40，第二计算模块50和第三计算模块60。
具体地，等效电路模块10，用于建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；戴维南等效电路模型用于表征系统侧电路的谐波电压源和系统侧谐波阻抗之间的连接关系。
数据采集模块20，用于在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组。
数据处理模块30，用于基于谐波电压观测数据组和谐波电流观测数据组，计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值。
第一计算模块40，用于基于戴维南等效电路模型得到系统侧电路的戴维南等效电路方程，并对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程。
第二计算模块50，用于根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；第一系统侧谐波电压源波动量函数用于描述电力系统公共连接点处的系统侧谐波电压源波动量的大小；以及
将电学关系方程代入到第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；第二系统侧谐波电压源波动量函数用于描述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系。
第三计算模块60，用于基于第二系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。
在本发明实施例中，通过在电力系统中，通过等效电路模块建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；通过数据采集模块在电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据；然后通过数据处理模块计算电力系统公共连接点处的谐波电压波动值和谐波电流波动值；通过第一计算模块对戴维南等效电路方程两端作差分运算，得到关于谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；第二计算模块根据2范数构建电力系统公共连接点处的系统侧谐波电压源波动量函数，最后通过第三计算模块基于系统侧谐波电压源波动量函数的最小值确定电力系统公共连接点处的系统侧谐波阻抗估计值。本发明通过对系统侧谐波电压源波动量的估计来计算系统侧谐波阻抗，这种方式能够使得到的系统侧谐波阻抗的估计值更为准确；同时通过以上方式能够以较少的计算步骤快速地得到电力系统的系统侧谐波阻抗估计值。
可选地，第三计算模块60还用于：
计算第二系统侧谐波电压源波动量函数的最小值；
确定第二系统侧谐波电压源波动量函数的数值为最小值时，系统侧谐波阻抗的取值，并将系统侧谐波阻抗的取值作为电力系统公共连接点处的系统侧谐波阻抗估计值。
可选地，图4为根据本发明实施例提供的另一种电力系统公共连接点处的系统侧谐波阻抗估计系统的示意图，如图4所示，该系统中的数据处理模块30还包括：谐波电压处理单元31和谐波电流处理单元32。
具体地，谐波电压处理单元31，用于通过公式
Figure BDA0002040750800000131
Figure BDA0002040750800000132
计算电力系统公共连接点处的谐波电压波动值；其中，
Figure BDA0002040750800000133
为第n个谐波电压波动值，
Figure BDA0002040750800000134
为谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure BDA0002040750800000135
为谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数；
谐波电流处理单元32，用于通过公式
Figure BDA0002040750800000136
计算电力系统公共连接点处的谐波电流波动值；其中，
Figure BDA0002040750800000137
为第n个谐波电流波动值，
Figure BDA0002040750800000138
为谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure BDA0002040750800000139
为谐波电流观测数据组中第n个采样数据的谐波电流值。
可选地，第一计算模块40还用于：
基于戴维南等效电路模型得到系统侧电路的戴维南等效电路方程为
Figure BDA00020407508000001310
对戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure BDA00020407508000001311
其中，Zu为系统侧谐波阻抗，
Figure BDA00020407508000001312
为系统侧谐波电压源波动量，
Figure BDA00020407508000001313
为系统侧谐波电压源电压值。
可选地，第二计算模块50还用于：
根据2范数构建数学表达式为
Figure BDA00020407508000001314
的第一系统侧谐波电压源波动量函数，其中，N为多个谐波电压数据或多个谐波电流数据的采集个数。
本发明实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。
本发明实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质，程序代码使处理器执行上述实施例一提供的方法。
本发明实施例所提供的进行一种电力系统公共连接点处的系统侧谐波阻抗估计方法的计算机程序产品，包括存储了处理器可执行的非易失的程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
实施例三：
本实施例是利用上述实施例一提供的方法计算系统侧谐波阻抗估计值的应用实施例。
采集到某100MVA直流电弧炉150kV母线(PCC)的数据，采样频率是6400Hz，每分钟对采样的电压电流数据进行快速傅里叶变换，得到各次谐波的测量值。将三次谐波电压电流值代入该方法计算程序中，得到计算结果如图5a和图5b所示，计算结果取均值得系统侧三次谐波阻抗幅值为10.0374，相角为68.0109。
实施例四：
本实施例是利用上述实施例一提供的方法计算系统侧谐波阻抗估计值的另一种应用实施例。
已知某城市电网负荷接入点处的信息，其电压等级为35kV，系统侧谐波阻抗为感性，用户侧谐波阻抗为容性。采集其PCC处某日11次谐波电压和电流数据，带入该方法的程序中，得到系统侧谐波阻抗在24小时内的变化结果，如图6a和图6b所示。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (10)

1.一种电力系统公共连接点处的系统谐波阻抗估计方法，包括：建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；所述戴维南等效电路模型用于表征所述系统侧电路的谐波电压源和系统侧谐波阻抗之间的连接关系；其特征在于，所述方法还包括：
在所述电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组；
基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值；
基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；
根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；所述第一系统侧谐波电压源波动量函数用于描述所述电力系统公共连接点处的系统侧谐波电压源波动量的大小；
将所述电学关系方程代入到所述第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；所述第二系统侧谐波电压源波动量函数用于描述所述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系；
基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值。
2.根据权利要求1所述的方法，其特征在于，基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值包括：
计算所述第二系统侧谐波电压源波动量函数的最小值；
确定所述第二系统侧谐波电压源波动量函数的数值为所述最小值时，所述系统侧谐波阻抗的取值，并将所述系统侧谐波阻抗的取值作为所述电力系统公共连接点处的系统侧谐波阻抗估计值。
3.根据权利要求1所述的方法，其特征在于，基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值包括：
通过以下公式计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值：
Figure FDA0002467615990000021
其中，
Figure FDA0002467615990000022
为第n个谐波电压波动值，
Figure FDA0002467615990000023
为所述谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure FDA0002467615990000024
为所述谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数；
Figure FDA0002467615990000025
其中，
Figure FDA0002467615990000026
为第n个谐波电流波动值，
Figure FDA0002467615990000027
为所述谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure FDA0002467615990000028
为所述谐波电流观测数据组中第n个采样数据的谐波电流值。
4.根据权利要求3所述的方法，其特征在于，基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程包括：
基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程为
Figure FDA0002467615990000029
对所述戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure FDA00024676159900000210
其中，Zu为系统侧谐波阻抗，
Figure FDA00024676159900000211
为系统侧谐波电压源波动量，
Figure FDA00024676159900000212
为系统侧谐波电压源电压值。
5.根据权利要求4所述的方法，其特征在于，根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数包括：
根据2范数构建数学表达式为
Figure FDA0002467615990000031
的第一系统侧谐波电压源波动量函数，其中，N为所述多个谐波电压数据或所述多个谐波电流数据的采集个数。
6.一种电力系统公共连接点处的系统谐波阻抗估计系统，包括：等效电路模块，所述等效电路模块，用于建立电力系统公共连接点处的系统侧电路的戴维南等效电路模型；所述戴维南等效电路模型用于表征所述系统侧电路的谐波电压源和系统侧谐波阻抗之间的连接关系；其特征在于，所述系统还包括：数据采集模块，数据处理模块，第一计算模块，第二计算模块和第三计算模块，其中，
所述数据采集模块，用于在所述电力系统公共连接点处，采集多个谐波电压数据和多个谐波电流数据，得到谐波电压观测数据组和谐波电流观测数据组；
所述数据处理模块，用于基于所述谐波电压观测数据组和所述谐波电流观测数据组，计算所述电力系统公共连接点处的谐波电压波动值和谐波电流波动值；
所述第一计算模块，用于基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程，并对所述戴维南等效电路方程两端作差分运算，得到关于所述谐波电压波动值、谐波电流波动值和系统侧谐波阻抗的电学关系方程；
所述第二计算模块，用于根据2范数构建所述电力系统公共连接点处的系统侧谐波电压源波动量函数，得到第一系统侧谐波电压源波动量函数；所述第一系统侧谐波电压源波动量函数用于描述所述电力系统公共连接点处的系统侧谐波电压源波动量的大小；以及
将所述电学关系方程代入到所述第一系统侧谐波电压源波动量函数，得到第二系统侧谐波电压源波动量函数；所述第二系统侧谐波电压源波动量函数用于描述所述系统侧电路的谐波电压源波动量的大小随系统侧谐波阻抗的变化关系；
所述第三计算模块，用于基于所述第二系统侧谐波电压源波动量函数的最小值确定所述电力系统公共连接点处的系统侧谐波阻抗估计值。
7.根据权利要求6所述的系统，其特征在于，所述第三计算模块还用于：
计算所述第二系统侧谐波电压源波动量函数的最小值；
确定所述第二系统侧谐波电压源波动量函数的数值为所述最小值时，所述系统侧谐波阻抗的取值，并将所述系统侧谐波阻抗的取值作为所述电力系统公共连接点处的系统侧谐波阻抗估计值。
8.根据权利要求6所述的系统，其特征在于，所述数据处理模块包括：谐波电压处理单元和谐波电流处理单元，其中，
所述谐波电压处理单元，用于通过公式
Figure FDA0002467615990000041
Figure FDA0002467615990000042
计算所述电力系统公共连接点处的谐波电压波动值；其中，
Figure FDA0002467615990000043
为第n个谐波电压波动值，
Figure FDA0002467615990000044
为所述谐波电压观测数据组中第n+1个采样数据的谐波电压值，
Figure FDA0002467615990000045
为所述谐波电压观测数据组中第n个采样数据的谐波电压值；n为大于零的正整数；
所述谐波电流处理单元，用于通过公式
Figure FDA0002467615990000046
计算所述电力系统公共连接点处的谐波电流波动值；其中，
Figure FDA0002467615990000047
为第n个谐波电流波动值，
Figure FDA0002467615990000048
为所述谐波电流观测数据组中第n+1个采样数据的谐波电流值，
Figure FDA0002467615990000051
为所述谐波电流观测数据组中第n个采样数据的谐波电流值。
9.根据权利要求8所述的系统，其特征在于，所述第一计算模块还用于：
基于所述戴维南等效电路模型得到所述系统侧电路的戴维南等效电路方程为
Figure FDA0002467615990000052
对所述戴维南等效电路方程两端作差分运算，得到电学关系方程为
Figure FDA0002467615990000053
其中，Zu为系统侧谐波阻抗，
Figure FDA0002467615990000054
为系统侧谐波电压源波动量，
Figure FDA0002467615990000055
为系统侧谐波电压源电压值。
10.根据权利要求9所述的系统，其特征在于，所述第二计算模块还用于：
根据2范数构建数学表达式为
Figure FDA0002467615990000056
的第一系统侧谐波电压源波动量函数，其中，N为所述多个谐波电压数据或所述多个谐波电流数据的采集个数。
CN201910341895.8A 2019-04-25 2019-04-25 电力系统公共连接点处的系统谐波阻抗估计方法和系统 Active CN110045186B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910341895.8A CN110045186B (zh) 2019-04-25 2019-04-25 电力系统公共连接点处的系统谐波阻抗估计方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910341895.8A CN110045186B (zh) 2019-04-25 2019-04-25 电力系统公共连接点处的系统谐波阻抗估计方法和系统

Publications (2)

Publication Number Publication Date
CN110045186A CN110045186A (zh) 2019-07-23
CN110045186B true CN110045186B (zh) 2020-06-23

Family

ID=67279439

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910341895.8A Active CN110045186B (zh) 2019-04-25 2019-04-25 电力系统公共连接点处的系统谐波阻抗估计方法和系统

Country Status (1)

Country Link
CN (1) CN110045186B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110261682A (zh) * 2019-08-14 2019-09-20 四川大学 实测数据相位缺失下的系统侧谐波阻抗获取方法和系统
CN110763920B (zh) * 2019-11-21 2021-11-09 国网四川省电力公司电力科学研究院 基于阻抗及电压范数最小的系统谐波阻抗估计方法及系统
EP3828557B1 (en) * 2019-11-28 2024-01-24 ABB Schweiz AG Determining thevenin equivalent model for a converter
CN114660362B (zh) * 2022-05-19 2022-08-16 四川大学 基于互信息数据优选的系统侧谐波阻抗估计方法及装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001004685A (ja) * 1999-06-22 2001-01-12 Chubu Electric Power Co Inc 電力系統の高調波特性測定方法
CN101706532B (zh) * 2009-11-25 2011-12-14 国网电力科学研究院武汉南瑞有限责任公司 一种谐波阻抗测量方法及测量装置
FR3042075B1 (fr) * 2015-10-06 2018-10-12 Universite De Poitiers Determination d'une pollution harmonique sur un reseau de distribution d'energie electrique
CN108957126A (zh) * 2018-06-20 2018-12-07 广西电网有限责任公司电力科学研究院 一种以谐波注入有功功率认定谐波污染责任的计算方法
CN108802491B (zh) * 2018-07-06 2019-07-19 四川大学 基于谐波源独立性的谐波发射水平估计求解新方法
CN109061300B (zh) * 2018-09-19 2020-11-10 江苏中凌高科技股份有限公司 一种pcc点的特征谐波源定位方法

Also Published As

Publication number Publication date
CN110045186A (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN110045186B (zh) 电力系统公共连接点处的系统谐波阻抗估计方法和系统
Melo et al. Harmonic state estimation for distribution networks using phasor measurement units
CN110031680B (zh) 一种系统侧谐波阻抗估计方法和系统
CN107342586B (zh) 一种用于确定配电网静态等值的方法及装置
CN108268999B (zh) 用于电力网络节点谐波特性分析的多主元模态分析法
CN105098763A (zh) 一种基于wams和ems的暂态电压稳定在线评估方法
CN108196127A (zh) 高压直流输电系统谐波阻抗扫描方法和装置
CN112305485B (zh) 电容式电压互感器的谐波电压测量误差修正方法及装置
Tripathy et al. Sparse S‐transform for location of faults on transmission lines operating with unified power flow controller
CN109541305A (zh) 一种谐波责任划分模型及谐波责任计算方法
CN112485598A (zh) 结合波形相关性与幅值特征的小电流接地故障定位方法
CN110333404A (zh) 一种非侵入式的负荷监测方法、装置、设备及存储介质
Artale et al. Measurement of Simplified Single-and Three-Phase Parameters for Harmonic Emission Assessment Based on IEEE 1459-2010
Samal et al. New signal subspace approach to estimate the inter‐area oscillatory modes in power system using TLS‐ESPRIT algorithm
Goh et al. Dynamic estimation of power system stability in different Kalman filter implementations
Ge et al. A dynamic parameter model of harmonic source networks
Zhou et al. Fault locating and severity assessment for power distribution systems based on elasticity network topology mapping
Zhurkin et al. Application of the Short Time Fourier Transform for Calculating Electric Energy Losses under the Conditions of a Non-Sinusoidal Voltage
Morsi A wavelet-based approach for reactive power metering in modern three-phase grids considering time-varying power quality disturbances
CN116613751A (zh) 一种新能源并网系统的小干扰稳定性分析方法及系统
CN112444675B (zh) 一种输电网节点的谐波超标预警方法、装置、设备及介质
CN111400649B (zh) 谐波责任量化方法、装置、计算机设备和存储介质
Bucci et al. Development of a low cost power meter based on a digital signal controller
CN102937676B (zh) 负荷谐波注入污染预警的实现方法及其系统
CN112698260B (zh) Cvt电压暂降测量误差校正方法及系统

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