CN106405203A - 基于最小二乘法的单相电压跌落检测方法 - Google Patents
基于最小二乘法的单相电压跌落检测方法 Download PDFInfo
- Publication number
- CN106405203A CN106405203A CN201610788799.4A CN201610788799A CN106405203A CN 106405203 A CN106405203 A CN 106405203A CN 201610788799 A CN201610788799 A CN 201610788799A CN 106405203 A CN106405203 A CN 106405203A
- Authority
- CN
- China
- Prior art keywords
- voltage
- square
- model
- error
- voltage drop
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16533—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application
- G01R19/16538—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application in AC or DC supplies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16566—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533
- G01R19/16576—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533 comparing DC or AC voltage with one threshold
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
本发明涉及一种基于最小二乘法的单相电压跌落检测方法,首先将含有各次谐波的周期电压信号表示成三角函数形式的傅里叶级数,建立电压模型;其次将得到的电压模型离散化,并用矩阵形式进行表示;然后根据最小二乘原理,通过不断缩小实际电压采样值与电压模型相减得到的误差的平方和,修正电压模型中各分量参数,得到符合实际情况的电压模型;最后再次利用最小二乘法原理对电压信号的实时采样值与建立模型的差值进行检测,利用协方差矩阵得到相应的基波及各次谐波的参数信息,将得到的电压基波幅值与设置的电压跌落补偿阀值进行比较,当基波幅值低于电压跌落补偿阀值时视为发生跌落。本发明稳态误差较小,敏感度高,可以快速检测到电压跌落。
Description
技术领域
本发明属于电力***状态监测领域,尤其涉及一种基于最小二乘法的单相电压跌落检测方法。
背景技术
动态电压恢复器(Dynamic Voltage Restorer,DVR)是一种用于保证敏感负荷供电电压稳定的有效串联补偿装置,可在毫秒级时间内向***与负荷之间串联注入幅值和相位可调的电压,补偿电压暂降,从而保证敏感负荷的电压在受到***电压扰动时仍处于可接受的范围内。
DVR为了满足检测的实时性要求,常见的三相暂降检测方法有dq变换、αβ变换等,上述算法虽然简单,但是检测的速度受谐波影响大。而单相检测方法有基于傅立叶变换的检测法、考虑相位跳变的检测法等,这些方法也都无法很好地解决由谐波造成的带宽增加、检测时间延长等问题。
发明内容
本发明的目的在于提供一种解决谐波存在时的快速检测电网电压跌落的检测方法,涉及到***中的谐波含量检测和跌落快速检测。
实现本发明目的的技术方案为:一种基于最小二乘法的单相电压跌落检测方法,包括以下步骤:
步骤1,将含有各次谐波的周期电压信号表示成三角函数形式的傅里叶级数,建立电压模型;
步骤2,将步骤1中得到的电压模型离散化,并用矩阵形式进行表示;
步骤3,根据最小二乘原理,通过不断缩小实际电压采样值与步骤2建立的电压模型相减得到的误差的平方和,修正电压模型中各分量参数,得到符合实际情况的电压模型;
步骤4,再次利用最小二乘法原理对电压信号的实时采样值与建立模型的差值进行检测,利用协方差矩阵得到相应的基波及各次谐波的参数信息,将得到的电压基波幅值与设置的电压跌落补偿阀值进行比较,当基波幅值低于电压跌落补偿阀值时视为发生跌落。
与现有技术相比,本发明的显著优点为:
(1)本发明减小了谐波模型与实际电网电压波形的差异,通过基于FFT的谐波分析,将模型参数进行修正,得到最优的模型参数;
(2)本发明利用优化后的模型参数,进行跌落电压检测,稳态误差较小,敏感度高,可以快速检测到电压跌落。
附图说明
图1为本发明的基于最小二乘法的单相电压跌落检测方法流程图。
图2(a)为电压跌落20%时基波检测波形图。
图2(b)为电压跌落20%时三次谐波检测波形图。
图2(c)为电压跌落20%时五次谐波检测波形图。
图2(d)为电压跌落20%时七次谐波检测波形图。
图3(a)为电压跌落60%时基波检测波形图。
图3(b)为电压跌落60%时三次谐波检测波形图。
图3(c)为电压跌落60%时五次谐波检测波形图。
图3(d)为电压跌落60%时七次谐波检测波形图。
具体实施方式
结合图1,本发明的一种基于最小二乘算法的电压跌落检测方法,包括以下步骤:
步骤1,将含有各次谐波的周期电压信号表示成三角函数形式的傅里叶级数,建立电压模型:
其中,u(t)为周期电压信号,周期为T1,基波角频率Un,ωn,分别为单相电网电压中对应的基波和各次谐波分量的幅值、角频率及起始相角,n=1,2,...N,N为最高谐波次数;
步骤2,将步骤1中得到的电压模型离散化,并用矩阵形式进行表示:
u(ti)=H(ti)x(ti)
H(ti)=[cosω1ti-sinω1ti…cosωNti-sinωNti] (2)
其中,u(ti)表示离散后的电压信号,H(ti)和x(ti)为(1)式中等式右边的展开矩阵表达形式;
步骤3,根据最小二乘原理,通过不断缩小实际电压采样值与步骤2建立的电压模型相减得到的误差的平方和,修正电压模型中各分量参数,得到符合实际情况的电压模型;具体为:
步骤3-1、将采样数据分成两组A(n)和B(n),长度分别为j1和j-j1,其中j是采样数据总长度,且j1>j-j1;
步骤3-2、计算A(n)的误差的平方和:
其中,f是遗忘因子,f∈(0,1),X(ti)为x(ti)的估计值,对EA进行最小化,得到相应的估计值X(ti);
步骤3-3、将步骤3-2得到的X(ti)带入组B(n)进行验证,相应的误差平方和为:
定义误差阈值ε1,若EB≥ε1则估计值误差大,重新调整谐波次数N和j1,然后重复步骤3-2和步骤3-3,若EB<ε1,则估计值作为***电压模型的参数;
步骤4,再次利用最小二乘法原理对电压信号的实时采样值与建立模型的差值进行检测,利用协方差矩阵得到相应的基波及各次谐波的参数信息,将得到的电压基波幅值与设置的电压跌落补偿阀值进行比较,当基波幅值低于电压跌落补偿阀值时视为发生跌落;具体为:
根据步骤3得到的估计值X(ti),最小二乘法算法表达式为:
其中,Hn(ti)=[cosωnti -sinωnti],Xn(ti)是状态向量xn(ti)的估计值;Pn(ti-1)是第i-1次迭代的二维协方差矩阵,初始值为p1...nI,p1...n=(p1,p1,p2,p2,…,pn,pn)为初始变量,I为单位矩阵;rn(ti)与kn(ti)为最小二乘法算式的中间变量;
单相电网电压的基波幅值为:
假设电压跌落补偿阀值为Usag,若Uamp_1>Usag则表明未发生跌落,若Uamp_1≤Usag则表明发生跌落;
最小二乘法的误差e(ti)为:
e(ti)=u(ti)-H(ti)X(ti) (7)
设置阈值ε2,当误差e(ti)大于阈值ε2时,判定当前电网电压出现暂降过程,将进入待补偿状态,同时将协方差矩阵Pn(ti-1)重置为初始值;当误差e(ti)小于等于阈值ε2时,则判定电网电压在正常范围内,不进行补偿,无需对协方差矩阵进行重置;
设置半个基波周期延时,当检测到误差e(ti)大于阈值ε2时重置协方差,同时进入延时程序,屏蔽协方差重置直到延时结束才能再次重置协方差。
下面结合具体实施例对本发明作进一步说明。
实施例
工业电网中的谐波通常是各种整流设备造成,这样的谐波次数为奇数,而且随着次数越高,谐波幅值越小。因此假如通过基于最小二乘法的谐波检测,得到实际***中的谐波分量,测得***中的谐波为3、5、7次,有效值分别为20V,10V,5V。建立的含谐波电压方程如下所示:
为了验证在不同跌落情况下检测的准度和快速性,分别设计了跌落20%及60%的实验。
图2为用测试仪模拟电网电压20%的跌落,即跌落至原来的80%,实线为采样波形,虚线为检测波形,图2(a)为基波量,图2(b)为三次谐波量,图2(c)为五次谐波量,图2(d)为七次谐波量;图中,标签值代表了横纵坐标值,横坐标代表采样点,纵坐标代表幅值。
图3为用测试仪模拟电网电压60%的跌落,即跌落至原来的40%,实线为采样波形,虚线为检测波形,图3(a)为基波量,图3(b)为三次谐波量,图3(c)为五次谐波量,图3(d)为七次谐波量;图中,标签值代表了横纵坐标值,横坐标代表采样点,纵坐标代表幅值。
以上波形图采用Matlab仿真得到,采样点的采样间隔为50us。从图2和图3中可以看出,检测精度在0.1V,检测速度为5-10ms,跌落深度对检测的速度影响不大。
Claims (5)
1.一种基于最小二乘法的单相电压跌落检测方法,其特征在于,包括以下步骤:
步骤1,将含有各次谐波的周期电压信号表示成三角函数形式的傅里叶级数,建立电压模型;
步骤2,将步骤1中得到的电压模型离散化,并用矩阵形式进行表示;
步骤3,根据最小二乘原理,通过不断缩小实际电压采样值与步骤2建立的电压模型相减得到的误差的平方和,修正电压模型中各分量参数,得到符合实际情况的电压模型;
步骤4,再次利用最小二乘法原理对电压信号的实时采样值与建立模型的差值进行检测,利用协方差矩阵得到相应的基波及各次谐波的参数信息,将得到的电压基波幅值与设置的电压跌落补偿阀值进行比较,当基波幅值低于电压跌落补偿阀值时视为发生跌落。
2.根据权利要求1所述的基于最小二乘法的单相电压跌落检测方法,其特征在于,步骤1中含有各次谐波的周期电压信号表示为三角函数形式的傅里叶级数:
其中,u(t)为周期电压信号,周期为T1,基波角频率Un,ωn,分别为单相电网电压中对应的基波和各次谐波分量的幅值、角频率及起始相角,n=1,2,...N,N为最高谐波次数。
3.根据权利要求2所述的基于最小二乘法的单相电压跌落检测方法,其特征在于,步骤2具体为:将式(1)中函数离散化变为离散形式,用矩阵形式表示:
其中,u(ti)表示离散后的电压信号,H(ti)和x(ti)为(1)式中等式右边的展开矩阵表达形式。
4.根据权利要求3所述的基于最小二乘法的单相电压跌落检测方法,其特征在于,步骤3具体包括以下步骤:
步骤3-1、将采样数据分成两组A(n)和B(n),长度分别为j1和j-j1,其中j是采样数据总长度,且j1>j-j1;
步骤3-2、计算A(n)的误差的平方和:
其中,f是遗忘因子,f∈(0,1),X(ti)为x(ti)的估计值,对EA进行最小化,得到相应的估计值X(ti);
步骤3-3、将步骤3-2得到的X(ti)带入组B(n)进行验证,相应的误差平方和为:
定义误差阈值ε1,若EB≥ε1则估计值误差大,重新调整谐波次数N和j1,然后重复步骤3-2和步骤3-3,若EB<ε1,则估计值作为***电压模型的参数。
5.根据权利要求4所述的基于最小二乘法的单相电压跌落检测方法,其特征在于,步骤4具体为:
根据步骤3得到的估计值X(ti),最小二乘法算法表达式为:
其中,Hn(ti)=[cosωnti-sinωnti],Xn(ti)是状态向量xn(ti)的估计值;Pn(ti-1)是第i-1次迭代的二维协方差矩阵,初始值为p1...nI,p1...n=(p1,p1,p2,p2,…,pn,pn)为初始变量,I为单位矩阵;rn(ti)与kn(ti)为最小二乘法算式的中间变量;
单相电网电压的基波幅值为:
假设电压跌落补偿阀值为Usag,若Uamp_1>Usag则表明未发生跌落,若Uamp_1≤Usag则表明发生跌落;
最小二乘法的误差e(ti)为:
e(ti)=u(ti)-H(ti)X(ti) (7)
设置阈值ε2,当误差大于阈值ε2时,判定当前电网电压出现暂降过程,将进入待补偿状态,同时将协方差矩阵Pn(ti-1)重置为初始值;当误差小于等于阈值ε2时,则判定电网电压在正常范围内,不进行补偿,无需对协方差矩阵进行重置;
设置半个基波周期延时,当检测到误差e(ti)大于阈值ε2时重置协方差,同时进入延时程序,屏蔽协方差重置直到延时结束才能再次重置协方差。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610222580 | 2016-04-11 | ||
CN2016102225808 | 2016-04-11 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106405203A true CN106405203A (zh) | 2017-02-15 |
Family
ID=58000523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610788799.4A Pending CN106405203A (zh) | 2016-04-11 | 2016-08-31 | 基于最小二乘法的单相电压跌落检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106405203A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507542A (zh) * | 2018-12-27 | 2019-03-22 | 国网北京市电力公司 | 一种单相电压暂降快速检测方法 |
CN109900975A (zh) * | 2019-02-27 | 2019-06-18 | 山东汇科工控技术有限公司 | 一种变频器输入电源缺相检测方法 |
CN110672974A (zh) * | 2019-09-27 | 2020-01-10 | 长沙理工大学 | 一种新型配电网电压跌落检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508008A (zh) * | 2011-10-25 | 2012-06-20 | 甘肃电力科学研究院 | 风力发电***并网母线电压跌落幅相检测***及方法 |
CN103838959A (zh) * | 2013-12-18 | 2014-06-04 | 国网上海市电力公司 | 偏最小二乘回归应用于配电网谐波源定位与检测中的方法 |
-
2016
- 2016-08-31 CN CN201610788799.4A patent/CN106405203A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508008A (zh) * | 2011-10-25 | 2012-06-20 | 甘肃电力科学研究院 | 风力发电***并网母线电压跌落幅相检测***及方法 |
CN103838959A (zh) * | 2013-12-18 | 2014-06-04 | 国网上海市电力公司 | 偏最小二乘回归应用于配电网谐波源定位与检测中的方法 |
Non-Patent Citations (1)
Title |
---|
雷何等: "基于包含谐波模型的加权最小二乘估计算法的电压暂降检测方法", 《电力自动化设备》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507542A (zh) * | 2018-12-27 | 2019-03-22 | 国网北京市电力公司 | 一种单相电压暂降快速检测方法 |
CN109900975A (zh) * | 2019-02-27 | 2019-06-18 | 山东汇科工控技术有限公司 | 一种变频器输入电源缺相检测方法 |
CN109900975B (zh) * | 2019-02-27 | 2022-01-11 | 山东汇科工控技术有限公司 | 一种变频器输入电源缺相检测方法 |
CN110672974A (zh) * | 2019-09-27 | 2020-01-10 | 长沙理工大学 | 一种新型配电网电压跌落检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101806832B (zh) | 一种低频率信号的频率测量方法 | |
Lopes et al. | A traveling-wave detection method based on Park's transformation for fault locators | |
CN101833036B (zh) | 一种交流电的瞬时相位测量方法 | |
US20140330443A1 (en) | Method for measuring frequency of phasor of power system | |
CN105092997B (zh) | 一种特高压输电线路雷电绕击与反击的识别方法 | |
CN1200177A (zh) | 测定配电网中单相接地故障位置的方法 | |
CN106018958B (zh) | 孤网变频***中二次侧电流电压频率跟踪方法 | |
CN101566663B (zh) | 一种配电***电压跌落源定位方法 | |
CN102879639A (zh) | 一种电力***中频率的实时测量方法 | |
CN102087332A (zh) | 一种基于波速寻优的直流行波故障测距方法 | |
CN106405203A (zh) | 基于最小二乘法的单相电压跌落检测方法 | |
CN101813725A (zh) | 一种低频率信号的相位差测量方法 | |
CN106324340B (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN102095929B (zh) | 一种快速测量交流电信号频率的方法 | |
CN104808035A (zh) | 一种基于广义双曲s变换的电压暂降检测方法 | |
CN103675544A (zh) | 基于优化算法的电力***故障信号检测与波形识别方法 | |
JP7133442B2 (ja) | 落雷電流波形推定方法および落雷電流波形推定装置 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN109284933A (zh) | 一种基于数理统计的电子式互感器状态评估***及方法 | |
CN102508026A (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN105183938A (zh) | 电网不良数据辨识与估计方法 | |
Zhao et al. | Recent development of frequency estimation methods for future smart grid | |
CN114016041B (zh) | 一种用于阴极保护的油气管道电位检测方法及其装置 | |
CN102508022A (zh) | 采用最优乘子牛顿算法检测电网频率的方法 | |
CN103543331A (zh) | 一种计算电信号谐波和间谐波的方法 |
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: 20170215 |