CN105004926A - 一种交流电相位频率幅值跟踪重构的方法 - Google Patents

一种交流电相位频率幅值跟踪重构的方法 Download PDF

Info

Publication number
CN105004926A
CN105004926A CN201510448078.4A CN201510448078A CN105004926A CN 105004926 A CN105004926 A CN 105004926A CN 201510448078 A CN201510448078 A CN 201510448078A CN 105004926 A CN105004926 A CN 105004926A
Authority
CN
China
Prior art keywords
phase
alternating current
amplitude
frequency
phase alternating
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.)
Granted
Application number
CN201510448078.4A
Other languages
English (en)
Other versions
CN105004926B (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201510448078.4A priority Critical patent/CN105004926B/zh
Publication of CN105004926A publication Critical patent/CN105004926A/zh
Application granted granted Critical
Publication of CN105004926B publication Critical patent/CN105004926B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
  • Particle Accelerators (AREA)

Abstract

本发明公开了一种交流电相位频率幅值跟踪重构方法,按照该重构方法,对采集的交流电实际值分析出交流电的频率、相位、相序关系,其后利用二阶广义积分器、频率和相位信息得到各相实时的估计幅值,再次估计幅值代入提出的算法,减小估计的频率和相位误差。按照本发明实现的交流电相位频率幅值跟踪重构的方法,无需提供前馈参考较频率,能够有效分析任何位置角频率的交流电相位、频率和幅值,并且对于单相或三相电压、电流信号均可实现分析,而且能够实现有效、快速的实时重构。

Description

一种交流电相位频率幅值跟踪重构的方法
技术领域
本发明属于交流电相位检测与重构的技术领域,更具体地,涉及一种交流电相位频率幅值跟踪重构的方法。
背景技术
光伏、风力等新能源并网发电快速发展,分布式发电成为热点研究对象。软件锁相环是分布式发电系统中的关键技术,用以检测污染的电网电压的频率和相位,为并网发电提供支持。受污染的电网不仅电压不平衡,而且含有多次谐波分量。为越来越多的学者对其进行了研究,提出了诸多锁相环的控制与设计方案。例如基于低通滤波锁相环、基于加权最小二乘法估计的锁相环、基于扩展卡尔曼滤波器的锁相环、基于空间矢量滤波器的锁相环等等。他们都存在以下问题:锁相的精度不高,响应慢,对系统频率变化和三相电压不平衡较敏感。随着研究的深入,许多锁相环被提出,从而使锁相环的性能得到不断的改善和提升。基于单同步坐标系软件锁相环,当三相电网电压不平衡时,得到的电压正序分量的相位存在较大2次谐波。基于对称分量法的单同步坐标系软件锁相环,对频率适应性相对不足。为了提高锁相环的锁相精度和响应,提出了基于双同步坐标系的解耦软件锁相环。该方法能在电压平衡和不平衡条件下准确地获取电网电压的相位、幅值和频率信息,具有较好的电网适应性。
由于单相电网电压信号不存在构成同步旋转坐标系的两相正交电压信号,上述基于同步旋转坐标系的三相电网锁相环控制策略一般不适用于单相电网锁相环的控制。单相电网电压的锁相环方案有:基于延迟法构造虚拟两相的单相锁相环,存在周期性的延迟,动态响应慢。基于P-Q理论的单相锁相环,在电压跌落、相位突变时,响应慢。基于Park反变换的单相锁相环,滤波器设计困难,电压不平衡时响应慢。
基于双二阶广义积分器的检测算法在许多方面都用应用。基于双二阶广义积分器的三相软件锁相环能在电压平衡和不平衡条件下分析三相电网电压的相位、幅值和频率信息,响应快。而且基于双二阶广义积分器的单相软件锁相环也能准确快速的获取单相交流电的信息。但是,一方面,以上锁相环技术均需要提供被测对象的基频和三相电压的相序;另一方面,当交流电谐波分量较大时,实际的相位曲线不是标准的正弦曲线,角频率不为恒定值,而以上算法得到的角频率均为恒定值,这样相位误差就会变大,在变频调速、未知电网和锁相环相位跟踪等复杂的应用情况下,上述算法并不能实际解决问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于二阶广义积分器交流电相位频率跟踪重构的算法,其目的在于无需给定参考频率,可以对任意未知幅值、频率的交流电分析,从而得出各自的幅值、频率和相位,具备良好的跟踪和重构效果。
为实现上述目的,按照本发明的一个方面,提供了一种交流电相位频率幅值跟踪重构方法,其特征在于,该方法包括如下步骤:
(1)实时采集交流电,判断其是单相或是三相交流电,若是单相,则进入步骤(2),若是三相交流电,则进入步骤(3);
(2)实现单相交流电的重构;其包括如下子步骤:
(2-1)连续采集N组单相交流电实时数据,并将其输入二阶广义积分器求解出两个正交的输出量,利用上述两个正交的输出量求出实时交流电幅值;
(2-2)求解其中θ为单相交流电实际相位角,为单相交流电跟踪相位角,其中其中X、Y是分别与所述两个正交的输出量相关的系数;
(2-3)将所述步骤(2-2)获得的作为PI控制器的输入,进行调节,进而输出所述单相交流电的角频率,所述角频率积分之后获得所述单相交流电的相位角;
(2-4)由上述步骤获得的所述幅值、角频率和相位角,实现所述单相交流电的重构;
(3)实现三相交流电的重构;其包括如下子步骤:
(3-1)连续采集N组三相交流电实时数据,计算得出其中θ为三相交流电实际相位角,为三相交流电跟踪相位角,其中M、N是与三相电压幅值和三相电压信号相关的表达式;
(3-2)将所述步骤(3-1)获得的作为PI控制器的输入,进行调节,进而输出所述三相交流电的角频率,所述角频率积分之后获得所述三相交流电的相位角;
(3-3)将所述步骤(3-2)的输出输入二阶广义积分器,求解得出三组两个正交的输出量,计算得出各相实时幅值;
(3-4)由所述步骤(3-2)、(3-3)获得的所述幅值、角频率和相位角,实现所述单相交流电的重构。
进一步地,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为其中t和l分别为求出的所述两正交的输出量。
进一步地,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为反Park变换法,即
总体而言,通过本发明所构思的以上技术方案与现有技术相比,由于对采集的三相交流电实际值分析出三相交流电的频率、相位、相序关系,其次利用二阶广义积分器、频率和相位信息得到各相实时的估计幅值,再次估计幅值代入提出的算法,减小估计的频率和相位误差,能够取得下列有益效果:
(1)巧妙地应用了正弦差公式及二阶广义积分器,无需提供前馈参考角频率,能有效分析任意未知角频率的交流电的相位、频率和幅值;
(2)对于单相或三相电压、电流信号,均可实现信号分析;
(3)本发明能适应含有谐波和噪声的交流电信号,并分析出准确的交流电的相位、频率和幅值;
(4)本发明实用有效、计算量小,在0.2s内即计算出交流电的频率、相位和幅值,并且实现交流电的实时重构。
附图说明
图1是按照本发明实现的针对三相交流电锁相环原理框图;
图2是按照本发明实现的针对单相交流电锁相环原理框图;
图3是按照本发明实现的交流电跟踪重构方法流程图;
图4是发电机实际控制绕组三相电流的观测图;
图5是按照本发明实现的三相锁相环计算得出的相位和频率的观测图;
图6是按照本发明的方法重构的B相电流与实际滤波后的三相电流的观测图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
将三相交流电压用正弦函数表示,同理也可以用余弦函数分析求解。三相交流电压表达式如下所示:
va+vb+vc=0   (4)
A、B、C分别表示三相交流电压幅值,当三相电压幅值相等时,满足A=B=C;当三相交流电压不平衡,幅值不相等时。忽略电压幅值恒为零情况,式(1)-(3)也可表示为:
由(5)-(7)可得:
其中为a相电压实际相位角;为a相电压算法跟踪相位角和角频率。当未知三相交流电的相序时,可使用正序锁相环算法。如果计算出来的角频率为负,则三相交流电为负序;如果计算出来的角频率为正,则三相交流电为正序。
以a相电压作为输入信号为例,单相交流电的锁相算法等式为(13)-(14)。系统达到稳态时,输入信号v和输出信号v′、qv′角频率相等幅值相等|v|=|v′|=|qv′|,v′与qv′正交,v′相位角与输入信号v一致且超前qv′相位角90°。
利用改进后的二阶广义积分器,有两种方法求输入信号的幅值。以a相电压作为输入信号为例,第一种是直接求解方法式(15)所示。
第二种是反Park变换方法求解,如式(16)所示。
作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,角频率积分之后得到交流电的相位角。
综上所诉,可以得出的三相交流电锁相环原理框图,如图1所示;单相交流电锁相环原理框图,如图2所示。
基于上述分析,如图3所示的流程图,按照本发明实施例的交流电频率、相位和幅值跟踪与重构方法包括如下步骤:
(1)根据交流电的实际状况,判断交流电是否单相,是则顺序执行如下步骤(A1)~(A4);否则说明交流电是三相交流电,顺序执行如下步骤(B1)~(B4)。
步骤(A1)~(A4)分别为:
(A1)连续采集N组单相交流电实时数据D1~DN,其中,第k组电机运行数据Dk,利用二阶广义积分器求解得出两个正交的输出量。求解利用式(15)或式(16)计算得出实时交流电幅值。
(A2)利用式(13)和式(14)计算得出
(A3)作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,角频率积分之后得到交流电的相位角。
(A4)利用求解得出的交流电幅值,角频率和相位角,实现交流电的重构。
步骤(B1)~(B4)分别为:
(B1)连续采集N组三相交流电实时数据D1~DN,其中,第k组电机运行数据Dk=(va,vb,vc),利用式(9)到式(12)计算得出
(B2)作为一个PI控制器的输入,进行调节。PI控制器的输出为交流电的角频率,如果计算得出的角频率为负数,则三相交流电的相许未负;反之为正。角频率积分之后得到交流电的相位角。
(B3)利用二阶广义积分器,将三相交流电求解得出三组两个正交的输出量。再利用式(15)或式(16)计算得出交流电各相实时幅值。
(B4)利用求解得出的三相交流电各相幅值,角频率和相位角,实现三相交流电的重构。
为使本领域技术人员更好地理解本发明,下面结合具体实例对本发明的交流电频率、相位和幅值跟踪与重构方法进行详细说明。
利用发电实验平台,使用提出的方法跟踪和重构发电机控制绕组三相电流。发电机控制绕组三相电流作为研究对象,发电机的同步转速为750RPM。发电机的转速从670RPM匀加速到860RPM,控制绕组的给定电流频率从负序的5.33Hz匀减为0,再从正序的0匀加速到7.33Hz,采样周期为0.1ms,计算时间4s。
采样的三相电流如图4所示,A相(1通道所示),B相(2通道所示),C相(3通道所示)。实际控制绕组三相电流是变频器直接输出的,不仅含有测量噪声,而且含有高次谐波分量。滤波后的三相电流作为提出的算法的输入量。经过上述方法,利用式(9)到式(12)计算得出PI控制器的输入量。进行调节计算得出三相电流的角频率,角频率积分之后得到三相电流的相位角,如图5所示,三相锁相环计算得出的相位(3通道所示)、未滤波频率(2通道所示)和滤波后的频率(4通道所示)。利用二阶广义积分器和式(15)或式(16)计算得出交流电各相实时幅值。利用求解得出的三相交流电各相幅值,角频率和相位角,实现三相交流电的重构,如图6所示,三相锁相环算法重构的B相电流(4通道所示)与实际滤波后的三相电流,A相(1通道所示),B相(2通道所示),C相(3通道所示)。按照本发明实现的重构方法的三相锁相环运行0.2s后,可实时跟踪电网频率、相位和电压幅值。未滤波的频率波动在±0.3Hz,而滤波后的频率波动在±0.05Hz,相位跟踪误差很小。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种交流电相位频率幅值跟踪重构方法,其特征在于,该方法包括如下步骤:
(1)实时采集交流电,判断其是单相或是三相交流电,若是单相,则进入步骤(2),若是三相交流电,则进入步骤(3);
(2)实现单相交流电的重构;其包括如下子步骤:
(2-1)连续采集N组单相交流电实时数据,并将其输入二阶广义积分器求解出两个正交的输出量,利用上述两个正交的输出量求出实时交流电幅值;
(2-2)求解其中θ为单相交流电实际相位角,为单相交流电跟踪相位角,其中其中X、Y是分别与所述两个正交的输出量相关的系数;
(2-3)将所述步骤(2-2)获得的作为PI控制器的输入,进行调节,进而输出所述单相交流电的角频率,所述角频率积分之后获得所述单相交流电的相位角;
(2-4)由上述步骤获得的所述幅值、角频率和相位角,实现所述单相交流电的重构;
(3)实现三相交流电的重构;其包括如下子步骤:
(3-1)连续采集N组三相交流电实时数据,计算得出其中θ为三相交流电实际相位角,为三相交流电跟踪相位角,其中M、N是与三相电压幅值和三相电压信号相关的表达式;
(3-2)将所述步骤(3-1)获得的作为PI控制器的输入,进行调节,进而输出所述三相交流电的角频率,所述角频率积分之后获得所述三相交流电的相位角;
(3-3)将所述步骤(3-2)的输出输入二阶广义积分器,求解得出三组两个正交的输出量,计算得出各相实时幅值;
(3-4)由所述步骤(3-2)、(3-3)获得的所述幅值、角频率和相位角,实现所述单相交流电的重构。
2.如权利要求1所述的交流电相位频率幅值跟踪重构方法,其特征在于,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为其中t和l分别为求出的所述两正交的输出量。
3.如权利要求1所述的交流电相位频率幅值跟踪重构方法,其特征在于,所述步骤(2-1)及所述步骤(3-3)中计算所述幅值的方法为反Park变换法,即
CN201510448078.4A 2015-07-27 2015-07-27 一种交流电相位频率幅值跟踪重构的方法 Expired - Fee Related CN105004926B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510448078.4A CN105004926B (zh) 2015-07-27 2015-07-27 一种交流电相位频率幅值跟踪重构的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510448078.4A CN105004926B (zh) 2015-07-27 2015-07-27 一种交流电相位频率幅值跟踪重构的方法

Publications (2)

Publication Number Publication Date
CN105004926A true CN105004926A (zh) 2015-10-28
CN105004926B CN105004926B (zh) 2018-01-26

Family

ID=54377669

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510448078.4A Expired - Fee Related CN105004926B (zh) 2015-07-27 2015-07-27 一种交流电相位频率幅值跟踪重构的方法

Country Status (1)

Country Link
CN (1) CN105004926B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108226588A (zh) * 2017-11-01 2018-06-29 中国矿业大学(北京) 一种适用于单相和三相电力系统的软件锁相方法
CN108306569A (zh) * 2018-02-02 2018-07-20 华中科技大学 基于广义积分器的永磁同步电机无速度控制方法及系统
CN110635484A (zh) * 2019-11-13 2019-12-31 平顶山学院 一种单相有源滤波器指令电流的生成方法、装置及设备
CN110907030A (zh) * 2019-10-14 2020-03-24 中国电力科学研究院有限公司 一种用于测量三相电容器装置噪声的方法及系统
CN112325933A (zh) * 2020-10-29 2021-02-05 贵州电网有限责任公司 一种基于边缘计算的多能源系统采集及预处理方法和终端
CN113224993A (zh) * 2021-05-25 2021-08-06 中国神华能源股份有限公司哈尔乌素露天煤矿 矿用电动轮重启方法和系统
CN113447715A (zh) * 2021-07-30 2021-09-28 燕山大学 一种2次谐波式电网同步锁频环方法
CN115795253A (zh) * 2021-09-09 2023-03-14 武汉市聚芯微电子有限责任公司 一种频率跟踪方法、装置、设备以及计算机存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102185319A (zh) * 2011-05-10 2011-09-14 重庆新世杰电气股份有限公司 一种有源电力滤波器的谐波检测锁相方法和装置
JP5255930B2 (ja) * 2008-07-04 2013-08-07 新電元工業株式会社 位相検出器
CN103986458A (zh) * 2014-04-16 2014-08-13 国家电网公司 一种基于重复控制的微电网单相并网锁相环控制方法
CN104158540A (zh) * 2014-08-27 2014-11-19 河北大学 一种三相数字锁相环及锁相方法
CN104410408A (zh) * 2014-12-25 2015-03-11 国家电网公司 一种单相锁相方法及装置
CN104600696A (zh) * 2015-01-05 2015-05-06 上海电力学院 一种用于不平衡条件下的软件锁相环实现方法
CN104786873A (zh) * 2015-03-25 2015-07-22 西南交通大学 一种基于锁相环的电压重构方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5255930B2 (ja) * 2008-07-04 2013-08-07 新電元工業株式会社 位相検出器
CN102185319A (zh) * 2011-05-10 2011-09-14 重庆新世杰电气股份有限公司 一种有源电力滤波器的谐波检测锁相方法和装置
CN103986458A (zh) * 2014-04-16 2014-08-13 国家电网公司 一种基于重复控制的微电网单相并网锁相环控制方法
CN104158540A (zh) * 2014-08-27 2014-11-19 河北大学 一种三相数字锁相环及锁相方法
CN104410408A (zh) * 2014-12-25 2015-03-11 国家电网公司 一种单相锁相方法及装置
CN104600696A (zh) * 2015-01-05 2015-05-06 上海电力学院 一种用于不平衡条件下的软件锁相环实现方法
CN104786873A (zh) * 2015-03-25 2015-07-22 西南交通大学 一种基于锁相环的电压重构方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王鹿军: "分布式发电中三相三电平并网逆变器的若干关键技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108226588A (zh) * 2017-11-01 2018-06-29 中国矿业大学(北京) 一种适用于单相和三相电力系统的软件锁相方法
CN108306569A (zh) * 2018-02-02 2018-07-20 华中科技大学 基于广义积分器的永磁同步电机无速度控制方法及系统
CN110907030A (zh) * 2019-10-14 2020-03-24 中国电力科学研究院有限公司 一种用于测量三相电容器装置噪声的方法及系统
CN110907030B (zh) * 2019-10-14 2022-10-04 中国电力科学研究院有限公司 一种用于测量三相电容器装置噪声的方法及系统
CN110635484A (zh) * 2019-11-13 2019-12-31 平顶山学院 一种单相有源滤波器指令电流的生成方法、装置及设备
CN112325933A (zh) * 2020-10-29 2021-02-05 贵州电网有限责任公司 一种基于边缘计算的多能源系统采集及预处理方法和终端
CN113224993A (zh) * 2021-05-25 2021-08-06 中国神华能源股份有限公司哈尔乌素露天煤矿 矿用电动轮重启方法和系统
CN113447715A (zh) * 2021-07-30 2021-09-28 燕山大学 一种2次谐波式电网同步锁频环方法
CN113447715B (zh) * 2021-07-30 2022-04-15 燕山大学 一种2次谐波式电网同步锁频环方法
CN115795253A (zh) * 2021-09-09 2023-03-14 武汉市聚芯微电子有限责任公司 一种频率跟踪方法、装置、设备以及计算机存储介质
CN115795253B (zh) * 2021-09-09 2023-10-03 武汉市聚芯微电子有限责任公司 一种频率跟踪方法、装置、设备以及计算机存储介质

Also Published As

Publication number Publication date
CN105004926B (zh) 2018-01-26

Similar Documents

Publication Publication Date Title
CN105004926A (zh) 一种交流电相位频率幅值跟踪重构的方法
CN106501574B (zh) 一种有源电力滤波器谐波电流检测方法
Qasim et al. Optimal current harmonic extractor based on unified ADALINEs for shunt active power filters
Rani et al. A three phase PLL with a dynamic feed forward frequency estimator for synchronization of grid connected converters under wide frequency variations
CN107706929B (zh) 基于最小方差滤波的自适应锁相环方法及系统
Yazdani et al. A real-time extraction of harmonic and reactive current in a nonlinear load for grid-connected converters
CN102401858A (zh) 一种电网电压基波分量及谐波分量的检测方法
CN110557118B (zh) 一种锁相装置及锁相方法
CN105425011B (zh) 一种适用于单相畸变电网的非线性幅相检测方法
CN101771361A (zh) 无交流电压传感器并网逆变器的直接功率控制方法
CN104184148A (zh) 一种同步旋转坐标系中谐波电流分次控制方法
CN106410858A (zh) 一种基于双dq坐标变换的软件数字锁相方法
CN109066678A (zh) 一种改进的apf系统谐波电流检测方法
CN105021872A (zh) 一种电网电压畸变不对称状态下电流各分量有功成分和无功成分的检测方法
CN104143837B (zh) 具有参数自适应特性的逆变器无交流电压传感器控制方法
CN103259290B (zh) 一种无锁相环的双馈发电机网侧变流器直接电压控制方法
Karimi-Ghartemani et al. A phasor measurement algorithm based on phase-locked loop
CN106787650B (zh) 电网不平衡下有功二次脉动抑制的直接电流指令计算方法
Nwobu et al. Grid voltage synchronization for unbalanced voltages using the energy operator
Mekhilef et al. Single-phase hybrid active power filter with adaptive notch filter for harmonic current estimation
CN113804949A (zh) 一种适用于宽变频交流发电系统的幅值确定方法及装置
CN111200287B (zh) 一种指定谐波电流注入装置及其谐波电流的给定方法
Muñoz-Aguilar et al. A sensor-less sliding mode control scheme for a stand-alone wound rotor synchronous generator under unbalanced load conditions
Guo et al. Improved park inverse transform algorithm for positive and negative sequence separation of grid voltage under unbalanced grid conditions
Popa et al. Improving the Shunt Active Power Filter Control Methods under Distorted and Unbalanced Grid Voltages

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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

Granted publication date: 20180126

Termination date: 20200727

CF01 Termination of patent right due to non-payment of annual fee