CN108957131B - 基于剩余电流的低压配电网谐波责任划分方法 - Google Patents
基于剩余电流的低压配电网谐波责任划分方法 Download PDFInfo
- Publication number
- CN108957131B CN108957131B CN201810880436.2A CN201810880436A CN108957131B CN 108957131 B CN108957131 B CN 108957131B CN 201810880436 A CN201810880436 A CN 201810880436A CN 108957131 B CN108957131 B CN 108957131B
- Authority
- CN
- China
- Prior art keywords
- distribution network
- low
- current
- voltage distribution
- harmonic
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000001939 inductive effect Effects 0.000 claims abstract description 58
- 238000005070 sampling Methods 0.000 claims description 78
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 12
- 238000002474 experimental method Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于剩余电流的低压配电网谐波责任划分方法，包括步骤：一、测量低压配电网PCC点的相电压和相电流；二、计算相电压的有效值和相电流的有效值；三、计算低压配电网PCC点的有功功率；四、计算低压配电网PCC点的等效电导；五、计算低压配电网PCC点的有功电流和有功电流的有效值；六、计算低压配电网PCC点的感性无功电流、感性无功电流的有效值、容性无功电流、容性无功电流的有效值、感性剩余电流、感性剩余电流的有效值、容性剩余电流和容性剩余电流的有效值；七、进行低压配电网谐波责任的划分。本发明方法计算量小，能够准确划分系统侧和用户侧谐波责任，为降低电能损耗、谐波治理及经济处罚提供了科学依据。
Description
技术领域
本发明属于配电网电能质量，具体涉及一种基于剩余电流的低压配电网谐波责任划分方法。
背景技术
随着电力电子技术的不断发展，大量的非线性设备接入到低压配电网中，使得电力负载侧谐波含量增多。谐波使得电压电流波形发生畸变，进而增加电力损耗，缩短用电设备寿命，对低压配电系统的安全、经济运行产生了严重威胁。另外，各种光伏、风电等新能源大量接入电网，使得低压配电网中的背景谐波电压增大。为了有效地控制电网中的谐波，避免出现电能质量纠纷问题，国际上提出了一种“奖惩性”方案，该方案能够执行的前提条件就是在公共连接点(point of common coupl ing，PCC)处能够正确划分系统以及用户各自的谐波责任。因此，正确判别谐波来源，明确电力用户负载及系统侧的谐波发射电流，对进一步改善电能质量、提高低压配电网的安全运行有着非常重要的意义。
目前，谐波责任划分方法主要是要区分两端的系统侧和用户侧对PCC点谐波畸变的影响，只能判断出系统侧和电力用户谁负主要谐波源责任，而无法明确计算双方各自的谐波发射量，而且由于谐波阻抗求取困难，进而造成谐波责任评估不准确。这类方法主要用于中压配电网的谐波责任划分，由于低压配电网背景谐波电压含量较高，使用该方法造成的误差更大，谐波责任划分结果准确性比较低。而且，传统的谐波责任划分方法大多是频域方法，要对电流电压进行频谱分解，计算工作量大，也无法准确划分系统侧和用户侧谐波责任。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于剩余电流的低压配电网谐波责任划分方法，其方法步骤简单方便，易于实现，不用对信号进行频谱分解，计算量小，能够准确划分系统侧和用户侧谐波责任，为降低电能损耗、谐波治理及经济处罚提供了科学依据。
为解决上述技术问题，本发明采用的技术方案是：一种基于剩余电流的低压配电网谐波责任划分方法，其特征在于，该方法包括以下步骤：
步骤一、测量低压配电网PCC点的相电压和相电流；
步骤二、计算低压配电网PCC点的相电压的有效值U和相电流的有效值I；
步骤三、计算低压配电网PCC点的有功功率P；
步骤四、计算低压配电网PCC点的等效电导G；
步骤五、计算低压配电网PCC点的有功电流和有功电流的有效值Ip；其中，将t时刻低压配电网PCC点的有功电流瞬时值记为ip(t)；
步骤六、计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqlr(t)和容性剩余电流的有效值Iqcr；
步骤七、进行低压配电网谐波责任的划分，具体方法为：
步骤701、根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，当判断为低压配电网系统中有谐波源时，再判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧；
步骤702、当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish；当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich；当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算系统侧的谐波电流大小Ish和用户侧的谐波电流大小Ich，并根据谐波电流大小进行谐波责任划分，谐波电流越大责任越大，谐波电流越小责任越小。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤一中所述测量低压配电网PCC点的相电压和相电流时，将一个采样周期内第n个相电压的瞬时采样值记为u(n)，将t时刻低压配电网PCC点的相电压的瞬时采样值记为u(t)；将一个采样周期内第n个相电流的瞬时采样值记为i(n)，将t时刻低压配电网PCC点的相电流的瞬时采样值记为i(t)；其中，n的取值为1～N的自然数，N为一个采样周期内相电压或相电流的采样点总数。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤二中计算低压配电网PCC点的相电压的有效值U采用的计算公式为
步骤二中计算低压配电网PCC点的相电流的有效值I采用的计算公式为
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤三中计算低压配电网PCC点的有功功率P采用的计算公式为
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤四中计算低压配电网PCC点的等效电导G采用的计算公式为
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤五中所述计算低压配电网PCC点的有功电流时，计算t时刻低压配电网PCC点的有功电流瞬时值ip(t)采用的计算公式为ip(t)＝G·u(t)；步骤五中所述计算低压配电网PCC点的有功电流的有效值Ip采用的计算公式为
ip(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的有功电流瞬时值。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤五之后还包括计算低压配电网PCC点的无功电流的步骤，将t时刻低压配电网PCC点的无功电流瞬时值记为iq(t)；iq(t)的计算公式为iq(t)＝i(t)-ip(t)。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤六中所述计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqcr(t)和容性剩余电流的有效值Iqcr的具体过程为：
步骤601、对低压配电网PCC点的相电压采用积分算子使相电压相位滞后半个采样周期，得到低压配电网PCC点的相电压积分瞬时值m(t)＝∫u(t)dt，将一个采样周期内第n个相电压的瞬时采样值对应的积分瞬时值记为m(n)；对低压配电网PCC点的相电压采用微分算子使电压相位超前半个采样周期，得到低压配电网PCC点的相电压微分瞬时值v(t)＝u′(t)，将一个采样周期内第n个相电压的瞬时采样值对应的微分瞬时值记为v(n)；
步骤602、根据公式
计算低压配电网PCC点的相电压积分的有效值M；根据公式
计算低压配电网PCC点的相电压微分的有效值V；
步骤603、根据公式
计算低压配电网PCC点的感性无功电流iql(t)，并根据公式
计算低压配电网PCC点的感性无功电流的有效值Iql；
步骤604、根据公式
计算低压配电网PCC点的容性无功电流iqc(t)，并根据公式
计算低压配电网PCC点的容性无功电流的有效值Iqc；
步骤605、根据公式iqlr(t)＝i(t)-ip(t)-iql(t)计算低压配电网PCC点的感性剩余电流iqlr(t)，并根据公式
计算低压配电网PCC点的感性剩余电流的有效值Iqlr，iqlr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的感性剩余电流瞬时值；
步骤606、根据公式iqcr(t)＝i(t)-ip(t)-iqc(t)计算低压配电网PCC点的容性剩余电流iqcr(t)，并根据公式
计算低压配电网PCC点的容性剩余电流的有效值Iqcr，iqcr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的容性剩余电流瞬时值。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤701中所述根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，采用的判断方法为：当Iqlr＝0或Iqcr＝0时，判断为低压配电网系统中无谐波源；当Iqlr≠0或Iqcr≠0时，判断为低压配电网系统中有谐波源；
步骤701中所述判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧，采用的方法为：当Iql＝-Iqc时，判断为低压配电网系统的谐波源来自用户侧；当Iql≠-Iqc且Iqlr＜＜Iql时，或当Iql≠-Iqc且Iqcr＜＜Iql时，判断为低压配电网系统的谐波源来自系统侧；当Iql≠-Iqc且Iqlr≈Iql时，或当Iql≠-Iqc且Iqcr≈Iql时，或当Iql≠-Iqc且Iqlr>Iql时，或当Iql≠-Iqc且Iqcr>Iql时，判断为低压配电网系统的谐波源同时来自系统侧和用户侧。
上述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤702中所述当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish，采用的计算公式为Ish＝Ih；其中，Ih为低压配电网的总谐波电流的有效值且
ih(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的低压配电网的总谐波电流的瞬时值，将t时刻低压配电网PCC点的相电流的瞬时采样值i(t)对应的低压配电网的总谐波电流的瞬时值记为ih(t)，
C32为A、B、C三相坐标轴转化到α、β坐标轴的矩阵，且
C为瞬时无功功率理论中求解瞬时有功电流ip′(t)和瞬时无功电流iq′(t)的系数矩阵且
ω为低压配电网中与A相电网电压同相位的正弦信号的角速度，
为瞬时有功电流ip′(t)的直流分量，
为瞬时无功电流iq′(t)的直流分量，
为t时刻低压配电网PCC点的A相电流的瞬时采样值，iB(t)为t时刻低压配电网PCC点的B相电流的瞬时采样值，iC(t)为t时刻低压配电网PCC点的C相电流的瞬时采样值；
步骤702中所述当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich，采用的计算公式为Ich＝Iqlr；
步骤702中所述当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算用户侧的谐波电流大小Ich用的计算公式为Ich＝Iqlr，计算系统侧的谐波电流大小Ish采用的计算公式为Ish＝Ih-Ich。
本发明与现有技术相比具有以下优点：
1、本发明利用低压配电网PCC点处测量到的相电压和相电流划分谐波责任，避免了网络参数和谐波阻抗的计算，简单方便，易于实现。
2、本发明针对传统的方法在含有背景谐波电压条件下谐波责任评估不准确的情况提出了一种解决方案，能够将低压配电网系统的背景谐波电压转化为相应的谐波剩余电流，实现系统侧谐波责任量化计算。
3、本发明提出了一种时域的谐波责任划分方法，不用对信号进行频谱分解，计算量小，且能够准确划分系统侧和用户侧谐波责任。
4、本发明针对现有方法只能判断出主要谐波源责任的情况，提出了一种基于剩余电流的谐波责任划分方法，能够实现低压配电网PCC点处系统侧和用户侧的谐波责任判别，并能准确量化双方谐波电流大小，为降低电能损耗、谐波治理及经济处罚提供了科学依据，能够对低压配电系统的安全、经济运行提供帮助。
综上所述，本发明的方法步骤简单方便，易于实现，不用对信号进行频谱分解，计算量小，能够准确划分系统侧和用户侧谐波责任，为降低电能损耗、谐波治理及经济处罚提供了科学依据。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明基于剩余电流的低压配电网谐波责任划分方法的方法流程框图。
图2A为本发明试验一中低压配电网PCC点电压波形图。
图2B为本发明试验一中低压配电网PCC点电流波形图。
图3A为本发明试验一中低压配电网PCC点容性无功电流波形图。
图3B为本发明试验一中低压配电网PCC点容性剩余电流波形图。
图3C为本发明试验一中低压配电网PCC点感性无功电流波形图。
图3D为本发明试验一中低压配电网PCC点感性剩余电流波形图。
图4A为本发明试验二中低压配电网PCC点电压波形图。
图4B为本发明试验二中低压配电网PCC点电流波形图。
图5A为本发明试验二中低压配电网PCC点容性无功电流波形图。
图5B为本发明试验二中低压配电网PCC点容性剩余电流波形图。
图5C为本发明试验二中低压配电网PCC点感性无功电流波形图。
图5D为本发明试验二中低压配电网PCC点感性剩余电流波形图。
图6A为本发明试验二中低压配电网PCC点电压波形图。
图6B为本发明试验二中低压配电网PCC点电流波形图。
图7A为本发明试验二中低压配电网PCC点容性无功电流波形图。
图7B为本发明试验二中低压配电网PCC点容性剩余电流波形图。
图7C为本发明试验二中低压配电网PCC点感性无功电流波形图。
图7D为本发明试验二中低压配电网PCC点感性剩余电流波形图。
具体实施方式
如图1所示，本发明的基于剩余电流的低压配电网谐波责任划分方法，包括以下步骤：
步骤一、测量低压配电网PCC点的相电压和相电流；
本实施例中，步骤一中所述测量低压配电网PCC点的相电压和相电流时，将一个采样周期内第n个相电压的瞬时采样值记为u(n)，将t时刻低压配电网PCC点的相电压的瞬时采样值记为u(t)；将一个采样周期内第n个相电流的瞬时采样值记为i(n)，将t时刻低压配电网PCC点的相电流的瞬时采样值记为i(t)；其中，n的取值为1～N的自然数，N为一个采样周期内相电压或相电流的采样点总数。
步骤二、计算低压配电网PCC点的相电压的有效值U和相电流的有效值I；
本实施例中，步骤二中计算低压配电网PCC点的相电压的有效值U采用的计算公式为
步骤二中计算低压配电网PCC点的相电流的有效值I采用的计算公式为
步骤三、计算低压配电网PCC点的有功功率P；
本实施例中，步骤三中计算低压配电网PCC点的有功功率P采用的计算公式为
步骤四、计算低压配电网PCC点的等效电导G；
本实施例中，步骤四中计算低压配电网PCC点的等效电导G采用的计算公式为
步骤五、计算低压配电网PCC点的有功电流和有功电流的有效值Ip；其中，将t时刻低压配电网PCC点的有功电流瞬时值记为ip(t)；
本实施例中，步骤五中所述计算低压配电网PCC点的有功电流时，计算t时刻低压配电网PCC点的有功电流瞬时值ip(t)采用的计算公式为ip(t)＝G·u(t)；步骤五中所述计算低压配电网PCC点的有功电流的有效值Ip采用的计算公式为
ip(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的有功电流瞬时值。
本实施例中，步骤五之后还包括计算低压配电网PCC点的无功电流的步骤，将t时刻低压配电网PCC点的无功电流瞬时值记为iq(t)；iq(t)的计算公式为iq(t)＝i(t)-ip(t)。
步骤六、计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqlr(t)和容性剩余电流的有效值Iqcr；
本实施例中，步骤六中所述计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqcr(t)和容性剩余电流的有效值Iqcr的具体过程为：
步骤601、对低压配电网PCC点的相电压采用积分算子使相电压相位滞后半个采样周期，得到低压配电网PCC点的相电压积分瞬时值m(t)＝∫u(t)dt，将一个采样周期内第n个相电压的瞬时采样值对应的积分瞬时值记为m(n)；对低压配电网PCC点的相电压采用微分算子使电压相位超前半个采样周期，得到低压配电网PCC点的相电压微分瞬时值v(t)＝u′(t)，将一个采样周期内第n个相电压的瞬时采样值对应的微分瞬时值记为v(n)；
步骤602、根据公式
计算低压配电网PCC点的相电压积分的有效值M；根据公式
计算低压配电网PCC点的相电压微分的有效值V；
步骤603、根据公式
计算低压配电网PCC点的感性无功电流iql(t)，并根据公式
计算低压配电网PCC点的感性无功电流的有效值Iql；
步骤604、根据公式
计算低压配电网PCC点的容性无功电流iqc(t)，并根据公式
计算低压配电网PCC点的容性无功电流的有效值Iqc；
步骤605、根据公式iqlr(t)＝i(t)-ip(t)-iql(t)计算低压配电网PCC点的感性剩余电流iqlr(t)，并根据公式
计算低压配电网PCC点的感性剩余电流的有效值Iqlr，iqlr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的感性剩余电流瞬时值；
步骤606、根据公式iqcr(t)＝i(t)-ip(t)-iqc(t)计算低压配电网PCC点的容性剩余电流iqcr(t)，并根据公式
计算低压配电网PCC点的容性剩余电流的有效值Iqcr，iqcr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的容性剩余电流瞬时值。
步骤七、进行低压配电网谐波责任的划分，具体方法为：
步骤701、根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，当判断为低压配电网系统中有谐波源时，再判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧；
本实施例中，步骤701中所述根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，采用的判断方法为：当Iqlr＝0或Iqcr＝0时，判断为低压配电网系统中无谐波源；当Iqlr≠0或Iqcr≠0时，判断为低压配电网系统中有谐波源；
步骤701中所述判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧，采用的方法为：当Iql＝-Iqc时，判断为低压配电网系统的谐波源来自用户侧；当Iql≠-Iqc且Iqlr＜＜Iql时，或当Iql≠-Iqc且Iqcr＜＜Iql时，判断为低压配电网系统的谐波源来自系统侧；当Iql≠-Iqc且Iqlr≈Iql时，或当Iql≠-Iqc且Iqcr≈Iql时，或当Iql≠-Iqc且Iqlr>Iql时，或当Iql≠-Iqc且Iqcr>Iql时，判断为低压配电网系统的谐波源同时来自系统侧和用户侧。
步骤702、当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish；当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich；当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算系统侧的谐波电流大小Ish和用户侧的谐波电流大小Ich，并根据谐波电流大小进行谐波责任划分，谐波电流越大责任越大，谐波电流越小责任越小。
本实施例中，步骤702中所述当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish，采用的计算公式为Ish＝Ih；其中，Ih为低压配电网的总谐波电流的有效值且
ih(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的低压配电网的总谐波电流的瞬时值，将t时刻低压配电网PCC点的相电流的瞬时采样值i(t)对应的低压配电网的总谐波电流的瞬时值记为ih(t)，
C32为A、B、C三相坐标轴转化到α、β坐标轴的矩阵，且
C为瞬时无功功率理论中求解瞬时有功电流i′p(t)和瞬时无功电流iq′(t)的系数矩阵且
ω为低压配电网中与A相电网电压同相位的正弦信号的角速度，
为瞬时有功电流i′p(t)的直流分量，
为瞬时无功电流iq′(t)的直流分量，
iA(t)为t时刻低压配电网PCC点的A相电流的瞬时采样值，iB(t)为t时刻低压配电网PCC点的B相电流的瞬时采样值，iC(t)为t时刻低压配电网PCC点的C相电流的瞬时采样值；
步骤702中所述当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich，采用的计算公式为Ich＝Iqlr；
步骤702中所述当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算用户侧的谐波电流大小Ich用的计算公式为Ich＝Iqlr，计算系统侧的谐波电流大小Ish采用的计算公式为Ish＝Ih-Ich。
为了验证本发明能够产生的效果，进行了以下试验：
试验一、在低压配电网PCC点处，系统侧为谐波源，有谐波电压；电力用户侧为非谐波源，不发射谐波电流，低压配电网PCC点电压波形如图2A所示，低压配电网PCC点电流波形如图2B所示，可见电压波形、电流波形有畸变，低压配电网PCC点容性无功电流波形如图3A所示，低压配电网PCC点容性剩余电流波形如图3B所示，低压配电网PCC点感性无功电流波形如图3C所示，低压配电网PCC点感性剩余电流波形如图3D所示，根据本发明方法计算的计算结果如表1所示。
表1低压配电网PCC点的各电流有效值及总谐波电流有效值
根据表1列出的计算结果可知，系统的剩余电流不为零，感性和容性无功电流分量不相等，剩余电流远远小于无功电流时，低压配电网的谐波来自于系统侧，系统侧发射的谐波电流大小为Ish＝Ih＝0.13024(A)。
可见，该谐波电流是由背景谐波电压产生，本发明实现了将背景谐波电压转为谐波电流的目的。
试验二、在低压配电网PCC点处，系统侧不发射谐波，电力用户侧为谐波源，发射谐波电流。低压配电网PCC点处的电压波形如图4A所示，低压配电网PCC点电流波形如图4B所示，可见电流波形有畸变，低压配电网PCC点容性无功电流波形如图5A所示，低压配电网PCC点容性剩余电流波形如图5B所示，低压配电网PCC点感性无功电流波形如图5C所示，低压配电网PCC点感性剩余电流波形如图5D所示，根据本发明方法计算的计算结果如表2所示。
表2低压配电网PCC点的各电流有效值及总谐波电流有效值
根据表2列出的计算结果可知，系统的剩余电流不为零，感性和容性无功电流分量大小相等，剩余无功电流和无功电流大小为同一数量级，低压配电网的谐波来自于电力用户侧，用户侧发射的谐波电流大小为Ich＝Iqlr＝13.35073(A)。
可见，由于用户侧负载的非线性引起谐波电流。
试验三、在低压配电网PCC点处，系统侧有谐波电压；用户为谐波源，发射谐波电流，低压配电网PCC点处的电压波形如图6A所示，低压配电网PCC点电流波形如图6B所示，可见电压、电流波形有畸变，可知在PCC点测量到的电压和电流信号是用户及背景谐波电压造成的谐波耦合，无法区分谐波源的责任。低压配电网PCC点容性无功电流波形如图7A所示，低压配电网PCC点容性剩余电流波形如图7B所示，低压配电网PCC点感性无功电流波形如图7C所示，低压配电网PCC点感性剩余电流波形如图7D所示，根据本发明方法计算的计算结果如表3所示。
表3低压配电网PCC点的各电流有效值及总谐波电流有效值
根据表3列出的计算结果可知，系统的剩余电流不为零，感性和容性无功电流分量不相等，剩余无功电流大于无功电流时，系统的剩余电流主要来自于负载侧(非正弦电压引起的剩余电流太小可以忽略)，因此，电力用户侧发射的谐波电流大小为Ich＝Iqlr＝14.20683(A)，系统侧发射的谐波电流大小为Ish＝Ih-Ich＝15.07558-14.20683＝0.86875(A)。
综上所述，通过本发明的方法能够很好地判别谐波责任是来自于系统侧还是电力用户侧，即背景谐波电压还是非线性负载引起的谐波，本发明实现了将背景谐波电压转为谐波电流的目的，并且能量化各谐波电流的大小。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。
Claims (7)
1.一种基于剩余电流的低压配电网谐波责任划分方法，其特征在于，该方法包括以下步骤：
步骤一、测量低压配电网PCC点的相电压和相电流；
步骤二、计算低压配电网PCC点的相电压的有效值U和相电流的有效值I；
步骤三、计算低压配电网PCC点的有功功率P；
步骤四、计算低压配电网PCC点的等效电导G；
步骤五、计算低压配电网PCC点的有功电流和有功电流的有效值Ip；其中，将t时刻低压配电网PCC点的有功电流瞬时值记为ip(t)；
步骤六、计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqcr(t)和容性剩余电流的有效值Iqcr；
步骤七、进行低压配电网谐波责任的划分，具体方法为：
步骤701、根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，当判断为低压配电网系统中有谐波源时，再判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧；
步骤702、当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish；当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich；当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算系统侧的谐波电流大小Ish和用户侧的谐波电流大小Ich，并根据谐波电流大小进行谐波责任划分，谐波电流越大责任越大，谐波电流越小责任越小；
步骤一中所述测量低压配电网PCC点的相电压和相电流时，将一个采样周期内第n个相电压的瞬时采样值记为u(n)，将t时刻低压配电网PCC点的相电压的瞬时采样值记为u(t)；将一个采样周期内第n个相电流的瞬时采样值记为i(n)，将t时刻低压配电网PCC点的相电流的瞬时采样值记为i(t)；其中，n的取值为1～N的自然数，N为一个采样周期内相电压或相电流的采样点总数；
步骤六中所述计算低压配电网PCC点的感性无功电流iql(t)、感性无功电流的有效值Iql、容性无功电流iqc(t)、容性无功电流的有效值Iqc、感性剩余电流iqlr(t)、感性剩余电流的有效值Iqlr、容性剩余电流iqlr(t)和容性剩余电流的有效值Iqcr的具体过程为：
步骤601、对低压配电网PCC点的相电压采用积分算子使相电压相位滞后半个采样周期，得到低压配电网PCC点的相电压积分瞬时值m(t)＝∫u(t)dt，将一个采样周期内第n个相电压的瞬时采样值对应的积分瞬时值记为m(n)；对低压配电网PCC点的相电压采用微分算子使电压相位超前半个采样周期，得到低压配电网PCC点的相电压微分瞬时值v(t)＝u′(t)，将一个采样周期内第n个相电压的瞬时采样值对应的微分瞬时值记为v(n)；
步骤602、根据公式
计算低压配电网PCC点的相电压积分的有效值M；根据公式
计算低压配电网PCC点的相电压微分的有效值V；
步骤603、根据公式
计算低压配电网PCC点的感性无功电流iql(t)，并根据公式
计算低压配电网PCC点的感性无功电流的有效值Iql；
步骤604、根据公式
计算低压配电网PCC点的容性无功电流iqc(t)，并根据公式
计算低压配电网PCC点的容性无功电流的有效值Iqc；
步骤605、根据公式iqlr(t)＝i(t)-ip(t)-iql(t)计算低压配电网PCC点的感性剩余电流iqlr(t)，并根据公式
计算低压配电网PCC点的感性剩余电流的有效值Iqlr，iqlr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的感性剩余电流瞬时值；
步骤606、根据公式iqcr(t)＝i(t)-ip(t)-iqc(t)计算低压配电网PCC点的容性剩余电流iqcr(t)，并根据公式
计算低压配电网PCC点的容性剩余电流的有效值Iqcr，iqcr(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的容性剩余电流瞬时值；
步骤702中所述当判断为低压配电网系统的谐波源来自系统侧时，将低压配电网谐波责任划分给系统侧，并计算系统侧的谐波电流大小Ish，采用的计算公式为Ish＝Ih；其中，Ih为低压配电网的总谐波电流的有效值且
ih(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的低压配电网的总谐波电流的瞬时值，将t时刻低压配电网PCC点的相电流的瞬时采样值i(t)对应的低压配电网的总谐波电流的瞬时值记为ih(t)，
C32为A、B、C三相坐标轴转化到α、β坐标轴的矩阵，且
C为瞬时无功功率理论中求解瞬时有功电流i′p(t)和瞬时无功电流iq′(t)的系数矩阵且
ω为低压配电网中与A相电网电压同相位的正弦信号的角速度，
为瞬时有功电流i′p(t)的直流分量，
为瞬时无功电流iq′(t)的直流分量，
iA(t)为t时刻低压配电网PCC点的A相电流的瞬时采样值，iB(t)为t时刻低压配电网PCC点的B相电流的瞬时采样值，iC(t)为t时刻低压配电网PCC点的C相电流的瞬时采样值；
步骤702中所述当判断为低压配电网系统的谐波源来自用户侧时，将低压配电网谐波责任划分给用户侧，并计算用户侧的谐波电流大小Ich，采用的计算公式为Ich＝Iqlr；
步骤702中所述当判断为低压配电网系统的谐波源同时来自系统侧和用户侧时，计算用户侧的谐波电流大小Ich用的计算公式为Ich＝Iqlr，计算系统侧的谐波电流大小Ish采用的计算公式为Ish＝Ih-Ich。
2.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤二中计算低压配电网PCC点的相电压的有效值U采用的计算公式为
步骤二中计算低压配电网PCC点的相电流的有效值I采用的计算公式为
3.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤三中计算低压配电网PCC点的有功功率P采用的计算公式为
4.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤四中计算低压配电网PCC点的等效电导G采用的计算公式为
5.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤五中所述计算低压配电网PCC点的有功电流时，计算t时刻低压配电网PCC点的有功电流瞬时值ip(t)采用的计算公式为ip(t)＝G·u(t)；步骤五中所述计算低压配电网PCC点的有功电流的有效值Ip采用的计算公式为
ip(n)为一个采样周期内第n个相电流的瞬时采样值i(n)对应的有功电流瞬时值。
6.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤五之后还包括计算低压配电网PCC点的无功电流的步骤，将t时刻低压配电网PCC点的无功电流瞬时值记为iq(t)；iq(t)的计算公式为iq(t)＝i(t)-ip(t)。
7.按照权利要求1所述的基于剩余电流的低压配电网谐波责任划分方法，其特征在于：步骤701中所述根据低压配电网PCC点的感性剩余电流的有效值Iqlr和容性剩余电流的有效值Iqcr，判断低压配电网系统中是否有谐波源，采用的判断方法为：当Iqlr＝0或Iqcr＝0时，判断为低压配电网系统中无谐波源；当Iqlr≠0或Iqcr≠0时，判断为低压配电网系统中有谐波源；
步骤701中所述判断低压配电网系统的谐波源来自系统侧还是用户侧，还是同时来自系统侧和用户侧，采用的方法为：当Iql＝-Iqc时，判断为低压配电网系统的谐波源来自用户侧；当Iql≠-Iqc且Iqlr＜＜Iql时，或当Iql≠-Iqc且Iqcr＜＜Iql时，判断为低压配电网系统的谐波源来自系统侧；当Iql≠-Iqc且Iqlr≈Iql时，或当Iql≠-Iqc且Iqcr≈Iql时，或当Iql≠-Iqc且Iqlr＞Iql时，或当Iql≠-Iqc且Iqcr＞Iql时，判断为低压配电网系统的谐波源同时来自系统侧和用户侧。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810880436.2A CN108957131B (zh) | 2018-08-03 | 2018-08-03 | 基于剩余电流的低压配电网谐波责任划分方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810880436.2A CN108957131B (zh) | 2018-08-03 | 2018-08-03 | 基于剩余电流的低压配电网谐波责任划分方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108957131A CN108957131A (zh) | 2018-12-07 |
| CN108957131B true CN108957131B (zh) | 2020-07-28 |
Family
ID=64467560
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810880436.2A Expired - Fee Related CN108957131B (zh) | 2018-08-03 | 2018-08-03 | 基于剩余电流的低压配电网谐波责任划分方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108957131B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111781536B (zh) * | 2020-08-04 | 2023-03-28 | 蚌埠依爱消防电子有限责任公司 | 一种监测单相交流系统对地绝缘电阻及泄漏电容的方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6215316B1 (en) * | 1998-08-11 | 2001-04-10 | The Governor Of The University Of Alberta | Method and apparatus for measuring harmonic current sources in electric power distribution systems |
| CN103163373B (zh) * | 2013-04-01 | 2015-08-19 | 华北电力大学（保定） | 适用于电能质量数据交换格式的谐波污染责任的计算方法 |
| CN104407214A (zh) * | 2014-11-27 | 2015-03-11 | 云南电网公司电力科学研究院 | 一种谐波源辨识方法 |
| CN105759117B (zh) * | 2015-10-15 | 2018-08-24 | 华北电力大学 | 一种基于谐波分析综合等值电路的谐波责任量化方法 |
| CN105548736A (zh) * | 2015-12-02 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | 一种谐波责任分摊量化评估方法 |
| CN106093570B (zh) * | 2016-06-08 | 2018-10-09 | 南瑞（武汉）电气设备与工程能效测评中心 | 一种判断系统侧为主要谐波源的谐波溯源方法 |
| CN107064744A (zh) * | 2017-04-14 | 2017-08-18 | 国网上海市电力公司 | 一种谐波源定位方法 |
| CN107144732B (zh) * | 2017-06-29 | 2019-08-30 | 西安科技大学 | 基于集总功率的用户侧谐波源定位方法 |
| CN107167659B (zh) * | 2017-06-29 | 2019-08-13 | 西安科技大学 | 含背景谐波电压的谐波责任区分方法 |
| CN108279363B (zh) * | 2018-01-26 | 2019-01-29 | 四川大学 | 基于改进有功功率方向法的主谐波源确定方法及装置 |
-
2018
- 2018-08-03 CN CN201810880436.2A patent/CN108957131B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN108957131A (zh) | 2018-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108318852B (zh) | 一种智能电能表方波影响试验方法 | |
| Cataliotti et al. | Current transformers effects on the measurement of harmonic active power in LV and MV networks | |
| CN109633262A (zh) | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 | |
| CN204536528U (zh) | 一种充电桩直流电能表检定装置 | |
| CN104316894A (zh) | 一种电能表实际运行环境的仿真校验方法 | |
| WO2018018467A1 (zh) | 一种电器短路试验测量系统的校准方法 | |
| Abdullah et al. | A new vector draft method for harmonic source detection at point of common coupling | |
| Müller et al. | Characterization of small photovoltaic inverters for harmonic modeling | |
| CN106501755A (zh) | 一种基于动态负荷模型的智能电能表动态误差测量方法 | |
| Artale et al. | Measurement of simplified single-and three-phase parameters for harmonic emission assessment based on IEEE 1459-2010 | |
| Peterson et al. | Evaluation of harmonic distortion from multiple renewable sources at a distribution substation | |
| CN108957131B (zh) | 基于剩余电流的低压配电网谐波责任划分方法 | |
| CN106093811A (zh) | 一种基于最小二乘法的变压器铁心剩磁检测方法 | |
| Lei et al. | Behavior of voltage transformers under distorted conditions | |
| CN109188080A (zh) | 一种仿真试验平台及相量测量单元的误差测算方法 | |
| Lavers et al. | Real-time measurement of electric arc-furnace disturbances and parameter variations | |
| Demerdziev et al. | Analysis of errors in active power and energy measurements under random harmonic distortion conditions | |
| Artale et al. | Measurement uncertainty of harmonic emission indicators based on IEEE Std. 1459-2010 | |
| Yamada et al. | Uncertainty estimation and performance evaluation of a non-sinusoidal power measurement standard | |
| Tentzerakis et al. | Measurement of wind farm harmonic emissions | |
| CN208224473U (zh) | 一种电容型设备相对介损测试带电检测装置 | |
| CN110007129A (zh) | 一种应用于动态电能计量的三相电压实时估计方法 | |
| Arseneau et al. | A test system for calibrating flickermeters | |
| Xiaojuan et al. | Evaluation of digital energy meter error by Monte Carlo method | |
| Cipriani et al. | Uncertainty evaluation in the measurements for the electric power quality analysis |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200728 |