CN105137278B - 基于sogi的单相变压器短路参数在线实时辨识方法 - Google Patents
基于sogi的单相变压器短路参数在线实时辨识方法 Download PDFInfo
- Publication number
- CN105137278B CN105137278B CN201510617218.6A CN201510617218A CN105137278B CN 105137278 B CN105137278 B CN 105137278B CN 201510617218 A CN201510617218 A CN 201510617218A CN 105137278 B CN105137278 B CN 105137278B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- transformer
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004804 winding Methods 0.000 claims abstract description 35
- 238000012544 monitoring process Methods 0.000 claims abstract description 9
- 238000000605 extraction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
本发明公开一种基于SOGI的单相变压器短路参数在线实时辨识方法,建立单相变压器的T型等效电路,列出变压器短路阻抗参数方程和T型等效电路方程;分别得到T型等效电路在dq坐标系下的d轴方程和q轴方程;利用SOGI‑FLL模型提取各输入信号基波分量和基波正交分量,并在dq坐标系下获得各输入信号在d轴分量和q轴分量,计算得到单相变压器的短路阻抗值Rk和Xk,实现单相变压器绕组运行状态的在线监测。当变压器绕组发生形变和故障时,变压器的短路阻抗也会发生变化,通过实时测量变压器原边、副边的电压、电流来辨识变压器绕组的短路参数,实现对变压器绕组运行状态的在线监测,及早发现变压器隐患,及时切除故障点,避免运行事故的发生,延长变压器的使用寿命。
Description
技术领域
本发明涉及一种对单向变压器在线辨识方法,具体的说,涉及一种基于SOGI的单相变压器短路参数在线实时辨识方法。
背景技术
变压器是电力系统中非常重要的电气设备,它的安全运行直接关系到电力系统的供电可靠性。在变压器运行过程中,不可避免会发生各种形式的故障,变压器诸多故障中发生数量最多的一直是绕组故障,而诊断变压器绕组的方法中,频率响应法、阻抗分析法、低压脉冲法虽然有可取之处,但是都属于离线方法,不能及时发现变压器的绕组故障,不适合在线测量。
中国专利号CN200810011051公开一种电力变压器绕组参数在线实时辨识装置以及方法,具体公开如下特征:步骤4,阻抗在线辨识,变电站现场上位机通过对接收到的数据进行变压器绕组阻抗计算,利用最小二乘法实现绕组阻抗的在线识别,对变压器进行预警和故障报告。但当最小二乘法辨识的输入矩阵|AAΤ|=0时,参数具有不可辨识性,导致辨识结果值发散或者收敛到不正确的值上,且计算量大、对参数跟踪存在数据饱和问题。
发明内容
本发明的目的是提供一种基于SOGI的单相变压器短路参数在线实时辨识方法,实现对变压器绕组运行状态的在线监测,及早发现变压器隐患,及时切除故障点,避免运行事故的发生,延长变压器的使用寿命。
本发明的技术方案是:一种基于SOGI的单相变压器短路参数在线实时辨识方法,该辨识方法包括以下辨识步骤:
建立单相变压器的T型等效电路,根据T型等效电路列出变压器短路阻抗参数方程,同时建立T型等效电路方程;
在dq坐标系下,根据T型等效电路方程,将变压器的原边电压相量定相到d轴,分别得到T型等效电路在dq坐标系下的d轴方程和q轴方程;
利用SOGI-FLL模型提取各输入信号基波分量和基波正交分量,并在dq坐标系下,获得各输入信号在d轴分量和q轴分量,并带入至d轴方程和q轴方程中,计算得到单相变压器的短路阻抗值Rk和Xk,实现单相变压器绕组运行状态的在线监测。
上述T型等效电路的变压器短路阻抗参数方程为:
Rk=R1+R′2=R1+k2R2 (1)
Xk=X1σ+X′2σ=X1σ+k2X2σ
变压器的T型等效电路方程为:
公式(1)和公式(2)中,k为变压器的原副边匝数比;R1为变压器的原边绕组电阻;R′2为变压器的副边绕组电阻;X1σ为变压器原边的漏电抗;X′2σ为折算到原边以后的变压器副边的漏电抗;Rk为变压器的短路电阻;Xk为变压器的短路电抗;分别为变压器的原边电压相量和原边电流相量,为折算到一次侧以后的变压器的副边电压相量和副边电流相量。
上述T型等效电路在dq坐标系下的d轴方程和q轴方程分别为:
d轴方程:
-U2d-I2dR′2-I2qX′2σ+I1dR1+I1qX1σ=U1d (3)
q轴方程:
-U2q-I2qR′2-I2dX′2σ+I1qR1+I1dX1σ=0 (4)
公式(3)和公式(4)中,U1为变压器的原边电压;U1d为U1的d轴分量;U2为变压器的副边电压;U2d、U2q分别为U2的d轴分量、q轴分量;I1为变压器的原边电流;I2为变压器的副边电流;I1d、I1q分别为I1的d轴分量、q轴分量;I2d、I2q分别为I2的d轴分量、q轴分量。
上述单相变压器副边绕组向原边折算时,副边匝数N2折算为原边匝数N1,且R′2=R1且X′2σ=X1σ,带入至公式(3)和公式(4)得到:
采用SOGI-FLL模型提取各输入信号的基波分量及基波正交分量,SOGI-FLL模型包括SOGI-QSG模块和FLL模块两部分,SOGI-QSG模块中的估计误差εv和正交输出qv′的乘积εf作为FLL模块的输入,FLL模块中具有负增益为-γ的积分器,积分器根据εf的变化,逐渐调整输出角频率ω′,最终使ω′=ω,SOGI-FLL模型的状态方程为:
稳态时,ω=ω′,输出相量为:
公式(7)~公式(10)中,x1为输入信号基波分量;x2为输入信号基波正交分量;A和B为状态矩阵;输入信号v=Vsin(ωt+φ);v′为输入信号的估计值;qv′为v′的正交值;εv为估计误差;ω′为FLL模块估计出的输入信号基波角频率。
通过上述公式(10),获得任意输入信号在dq坐标系下的d轴分量和q轴分量,具体为:在dq坐标系下,将d轴定向于输入信号输入信号在d轴和q轴的分量分别为:
其中,和为输入信号;v1、qv1分别为的基波分量及的基波正交分量;v2、qv2分别为的基波及的基波正交分量。
本发明与现有技术相比的有益效果为:
当变压器绕组发生形变和故障时,变压器的短路阻抗也会发生变化,可以通过实时测量变压器原边、副边的电压、电流来辨识变压器绕组的短路参数,实现对变压器绕组运行状态的在线监测,从而及早发现变压器隐患,及时切除故障点,避免运行事故的发生,延长变压器的使用寿命。
附图说明
图1为本发明的T型等效电路;
图2为T型等效电路的相量图;
图3为SOGI-FLL模型的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
实施例1
本发明公开一种基于SOGI的单相变压器短路参数在线实时辨识方法,包括以下辨识步骤:
(一)建立单相变压器的T型等效电路,根据T型等效电路列出变压器短路阻抗参数方程,同时建立T型等效电路方程。
T型等效电路参见图1。图1中k为变压器变比,分别表示变压器的原边电压相量和原边电流相量,表示折算到一次侧以后的变压器的副边电压相量和副边电流相量,R1、L1σ、R′2、L′2σ分别表示变压器的原边绕组电阻和原边绕组漏感、副边绕组电阻和副边绕组漏感,Rm、Lm分别为变压器的激磁电阻和激磁电感。
图1中,变压器短路阻抗表示为:
Rk=R1+R′2=R1+k2R2 (1)
Xk=X1σ+X′2σ=X1σ+k2X2σ
其中:公式(1)中,
R′2=k2R2;
X′2σ=k2X2σ;
X1σ表示变压器原边的漏电抗;
X′2σ表示折算到原边以后的变压器副边的漏电抗;
Rk表示变压器的短路电阻;
Xk表示变压器的短路电抗。
图1中,变压器的原边电压相量表示为:
(二)在dq坐标系下,为了计算方便,根据T型等效电路方程,将变压器的原边电压相量定相到d轴,画出变压器T型等效电路的相量图如图2所示,并相量图得到T型等效电路在dq坐标系下的d轴方程和q轴方程。
根据图2列出公式(2)分别在d轴和q轴的方程,公式(2)分解如下:
d轴方程:
-U2d-I2dR′2-I2qX′2σ+I1dR1+I1qX1σ=U1d (3)
q轴方程:
-U2q-I2qR′2-I2dX′2σ+I1qR1+I1dX1σ=0 (4)
其中,公式(3)和公式(4)中,
U1表示变压器的原边电压;
U1d表示U1的d轴分量;
U2表示变压器的副边电压;
U2d、U2q分别表示U2的d轴分量、q轴分量;
I1表示变压器的原边电流;
I2表示变压器的副边电流;
I1d、I1q分别表示I1的d轴分量、q轴分量;
I2d、I2q分别表示I2的d轴分量、q轴分量。
在公式(3)和公式(4)中:
令:R′2=R1;X′2σ=X1σ;
联立公式(3)和公式(4)得到:
变压器短路阻抗表示为:
由公式(5)和公式(6)可知,测得变压器绕组原边电压、副边电压、原边电流和副边电流的信号幅值和相位,就能够实时地计算出变压器的短路阻抗参数Xk和Rk,从而实现变压器绕组运行状态的在线监测。
变压器绕组原边电压、副边电压、原边电流和副边电流的信号幅值和相位可通过下述方法得到。采用现有SOGI-FLL模型提取各输入信号基波分量和基波正交分量,并在dq坐标系下,获得各输入信号在d轴分量和q轴分量,并带入至d轴方程和q轴方程中,计算得到单相变压器的短路阻抗值Rk和Xk,实现单相变压器绕组运行状态的在线监测。
SOGI-FLL模型为二阶广义积分器-锁相环(Second Order GeneralizedIntegrator-Frequency Locked Loop)。
现有SOGI-FLL模型大都用于估算异步电机同步角频率,本发明利用SOGI-FLL模型提取各输入信号的基波分量和基波正交分量。图3为SOGI-FLL模型的结构框图,其主要由SOGI-QSG模块和FLL模块两部分组成。对于给定的正弦输入信号v=Vsin(ωt+φ),图3中v′为输入信号的估计值,qv′为v′的正交值,εv为估计误差,ω′为FLL模块估计出的输入信号基波角频率。
SOGI-QSG模块中的估计误差εv和正交输出qv′的乘积εf作为FLL模块的输入,当正弦输入信号v的角频率ω小于FLL模块的输入信号基波角频率ω′时,即ω<ω′时,qv′和εv同相位,εf>0;当正弦输入信号v的角频率ω大于FLL模块的输入信号基波角频率ω′时,即ω>ω′时,qv′和εv相位相反,εf<0;当ω=ω′时,εf=0。FLL模块中具有负增益为-γ的积分器,积分器根据εf的变化,逐渐调整输出角频率ω′,最终使ω′=ω。
图3所示SOGI-FLL模型的状态方程为:
公式(7)为基波提取算法的状态方程,其中和y分别表示SOGI-FLL模型的状态变量和输出变量,A和B为状态矩阵;当输入信号v为正弦信号时,通过对公式(7)求解可以得到公式(8)。
其中:
x1为输入信号基波分量;
x2为输入信号基波正交分量;
ω′为输入信号基波角频率。
对于给定的输入信号v=Vsin(ωt+φ),稳态时,ω=ω′,输出相量为:
由此可知,SOGI-FLL结构可以实现对输入信号基波分量x1及其输入信号基波正交分量x2的准确提取。
当电力系统存在谐波导致SOGI-FLL的输入信号含有谐波分量时,可以多个SOGI-QSG模块并联实现对各次谐波分量的准确提取。所以谐波条件下同样能够准确获得电压信号基波分量、电流信号基波分量及电压信号基波正交分量和电流信号基波正交分量。
对于输入信号和输入信号v1、qv1分别为的基波分量及的基波正交分量,v2、qv2分别为的基波及的基波正交分量。在dq坐标系下,将d轴定向于输入信号那么,输入信号在d轴和q轴的分量分别为:
因此本技术方案的原理为:利用SOGI-FLL模型提取各输入信号的基波分量及基波正交分量,然后根据公式(11)、(12)分别计算各输入信号的d轴分量和q轴分量,再带人式(5)、(6)就可以实时计算出单相变压器的短路参数Xk和Rk,从而实现单相变压器绕组运行状态的在线监测。
Claims (6)
1.一种基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于该辨识方法包括以下辨识步骤:
建立单相变压器的T型等效电路,根据T型等效电路列出变压器短路阻抗参数方程,同时建立T型等效电路方程;
在dq坐标系下,根据T型等效电路方程,将变压器的原边电压相量定相到d轴,分别得到T型等效电路在dq坐标系下的d轴方程和q轴方程;
利用SOGI-FLL模型提取各输入信号基波分量和基波正交分量,并在dq坐标系下,获得各输入信号在d轴分量和q轴分量,并带入至d轴方程和q轴方程中,计算得到单相变压器的短路阻抗值Rk和Xk,实现单相变压器绕组运行状态的在线监测。
2.根据权利要求1所述的基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于:所述T型等效电路的变压器短路阻抗参数方程为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>X</mi>
<mrow>
<mn>1</mn>
<mi>&sigma;</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mrow>
<mn>2</mn>
<mi>&sigma;</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>X</mi>
<mrow>
<mn>1</mn>
<mi>&sigma;</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<msub>
<mi>X</mi>
<mrow>
<mn>2</mn>
<mi>&sigma;</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
变压器的T型等效电路方程为:
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>jX</mi>
<mrow>
<mn>1</mn>
<mi>&sigma;</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>+</mo>
<msubsup>
<mi>jX</mi>
<mrow>
<mn>2</mn>
<mi>&sigma;</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(1)和公式(2)中,k为变压器的原副边匝数比;R1为变压器的原边绕组电阻;R′2为变压器的副边绕组电阻;X1σ为变压器原边的漏电抗;X′2σ为折算到原边以后的变压器副边的漏电抗;Rk为变压器的短路电阻;Xk为变压器的短路电抗;分别为变压器的原边电压相量和原边电流相量,为折算到一次侧以后的变压器的副边电压相量和副边电流相量。
3.根据权利要求2所述的基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于:所述T型等效电路在dq坐标系下的d轴方程和q轴方程分别为:
d轴方程:
-U2d-I2dR′2-I2qX′2σ+I1dR1+I1qX1σ=U1d (3)
q轴方程:
-U2q-I2qR′2-I2dX′2σ+I1qR1+I1dX1σ=0 (4)
公式(3)和公式(4)中,U1为变压器的原边电压;U1d为U1的d轴分量;U2为变压器的副边电压;U2d、U2q分别为U2的d轴分量、q轴分量;I1为变压器的原边电流;I2为变压器的副边电流;I1d、I1q分别为I1的d轴分量、q轴分量;I2d、I2q分别为I2的d轴分量、q轴分量。
4.根据权利要求3所述的基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于:所述单相变压器副边绕组向原边折算时,副边匝数N2折算为原边匝数N1,且R′2=R1、且X′2σ=X1σ,带入至公式(3)和公式(4)得到:
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mn>1</mn>
<mi>&sigma;</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>1</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
5.根据权利要求1所述的基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于:所述SOGI-FLL模型包括SOGI-QSG模块和FLL模块两部分,SOGI-QSG模块中的估计误差εv和正交输出qv′的乘积εf作为FLL模块的输入,FLL模块中具有负增益为-γ的积分器,积分器根据εf的变化,逐渐调整输出角频率ω′,最终使ω′=ω,SOGI-FLL模型的状态方程为:
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mo>+</mo>
<mi>B</mi>
<mi>v</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>k&omega;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>k&omega;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>v</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>v</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>qv</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msup>
<mi>&omega;</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&gamma;x</mi>
<mn>2</mn>
</msub>
<msup>
<mi>&omega;</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
稳态时,ω=ω′,输出相量为:
<mrow>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>v</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>qv</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>V</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>t</mi>
<mo>+</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mi>t</mi>
<mo>+</mo>
<mi>&phi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(7)~公式(10)中,x1为输入信号基波分量;x2为输入信号基波正交分量;A和B为状态矩阵;输入信号v=Vsin(ωt+φ);v′为输入信号的估计值;qv′为v′的正交值;εv为估计误差;输出角频率ω′为FLL模块估计出的输入信号基波角频率。
6.根据权利要求5所述的基于SOGI的单相变压器短路参数在线实时辨识方法,其特征在于:在dq坐标系下,将d轴定向于输入信号输入信号在d轴和q轴的分量分别为:
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<msub>
<mi>v</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>qv</mi>
<mn>1</mn>
</msub>
<msub>
<mi>qv</mi>
<mn>2</mn>
</msub>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>qv</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mn>2</mn>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>qv</mi>
<mn>2</mn>
</msub>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>v</mi>
<mn>2</mn>
</msub>
<msub>
<mi>qv</mi>
<mn>1</mn>
</msub>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>qv</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,和为输入信号;v1、qv1分别为的基波分量及的基波正交分量;v2、qv2分别为的基波及的基波正交分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510617218.6A CN105137278B (zh) | 2015-09-24 | 2015-09-24 | 基于sogi的单相变压器短路参数在线实时辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510617218.6A CN105137278B (zh) | 2015-09-24 | 2015-09-24 | 基于sogi的单相变压器短路参数在线实时辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137278A CN105137278A (zh) | 2015-12-09 |
CN105137278B true CN105137278B (zh) | 2018-01-09 |
Family
ID=54722693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510617218.6A Active CN105137278B (zh) | 2015-09-24 | 2015-09-24 | 基于sogi的单相变压器短路参数在线实时辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137278B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106324354B (zh) * | 2016-08-29 | 2017-05-31 | 张正 | 一种变压器漏阻抗检测方法及装置、故障检测方法及装置 |
CN106357183B (zh) * | 2016-09-22 | 2018-09-07 | 东南大学 | 一种直线振荡电机的谐振频率跟踪控制方法 |
CN106707028A (zh) * | 2016-12-30 | 2017-05-24 | 广东电网有限责任公司惠州供电局 | 变压器短路阻抗在线监测方法及装置 |
CN111077384B (zh) * | 2019-10-17 | 2021-12-10 | 华北电力大学(保定) | 一种变压器内部故障在线监测方法 |
CN110619192B (zh) * | 2019-10-17 | 2023-04-28 | 华北电力大学(保定) | 一种变压器参数在线计算方法 |
CN111123162A (zh) * | 2019-12-30 | 2020-05-08 | 国网山东省电力公司淄博供电公司 | 一种基于基波正序分量三相变压器短路参数在线监测方法 |
CN111044828B (zh) * | 2019-12-30 | 2022-04-22 | 国网山东省电力公司淄博供电公司 | 基于正、负序方程组的三相变压器绕组参数在线监测方法 |
CN111141956B (zh) * | 2019-12-31 | 2022-03-11 | 国网山东省电力公司淄博供电公司 | 一种基于微积分方程组单相变压器短路参数在线监测方法 |
CN111722152B (zh) * | 2020-06-29 | 2023-04-28 | 成都工百利自动化设备有限公司 | 一种变压器绕组变形监测方法及监测系统 |
CN112327217B (zh) * | 2020-10-13 | 2023-02-10 | 上海电力大学 | 基于变压器稳态和暂态运行的绕组漏电感在线辨识方法 |
CN113419196A (zh) * | 2021-05-07 | 2021-09-21 | 国网山东省电力公司淄博供电公司 | 一种基于负荷变化的单相变压器绕组参数在线监测方法 |
CN113625201B (zh) * | 2021-08-31 | 2024-03-29 | 山东理工大学 | 一种单相变压器绕组变形故障的在线检测方法 |
CN114156949B (zh) * | 2021-11-30 | 2024-07-09 | 中冶南方都市环保工程技术股份有限公司 | 一种单相光伏同步的方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101261297A (zh) * | 2008-04-17 | 2008-09-10 | 沈阳工业大学 | 电力变压器绕组参数在线实时辨识装置及方法 |
CN101930047A (zh) * | 2010-08-27 | 2010-12-29 | 西安交通大学 | 一种变压器绕组状态在线监测装置及其监测方法 |
CN102721898A (zh) * | 2012-02-27 | 2012-10-10 | 衢州电力局 | 变压器绕组变形在线测量方法及系统 |
CN103954865A (zh) * | 2014-05-06 | 2014-07-30 | 国家电网公司 | 一种变压器绕组机械状态在线监测装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008055499A2 (en) * | 2006-11-06 | 2008-05-15 | Gamesa Innovation & Technology, S.L. | Advanced real-time grid monitoring system and method |
-
2015
- 2015-09-24 CN CN201510617218.6A patent/CN105137278B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101261297A (zh) * | 2008-04-17 | 2008-09-10 | 沈阳工业大学 | 电力变压器绕组参数在线实时辨识装置及方法 |
CN101930047A (zh) * | 2010-08-27 | 2010-12-29 | 西安交通大学 | 一种变压器绕组状态在线监测装置及其监测方法 |
CN102721898A (zh) * | 2012-02-27 | 2012-10-10 | 衢州电力局 | 变压器绕组变形在线测量方法及系统 |
CN103954865A (zh) * | 2014-05-06 | 2014-07-30 | 国家电网公司 | 一种变压器绕组机械状态在线监测装置 |
Non-Patent Citations (4)
Title |
---|
基于二阶广义积分器-锁频环的异步电机同步角频率估计方法;辛振 等;《电工技术学报》;20140131;第29卷(第1期);第116-122页 * |
基于参数辨识的变压器绕组变形在线监测方法;邓祥力 等;《中国电机工程学报》;20141005;第34卷(第28期);第4950-4957页 * |
漏电抗的参数辨识技术在线监测变压器绕组变形的研究;欧小波 等;《高压电器》;20101231;第46卷(第12期);第41-44页 * |
电力变压器绕组状态实时监测算法;范竞敏 等;《电力自动化设备》;20100331;第30卷(第3期);第81-85页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105137278A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105137278B (zh) | 基于sogi的单相变压器短路参数在线实时辨识方法 | |
CN101981774B (zh) | 用于产生故障信号的方法、装置和现场设备 | |
CN103675445B (zh) | 一种利用电容式电压互感器进行高压谐波测量的方法 | |
CN102818973B (zh) | 一种基于相序电压幅值比较的故障选相方法及其装置 | |
CN102592003A (zh) | 一种用于机电与电磁暂态混合仿真的数据交换方法 | |
CN104007321B (zh) | 一种无功补偿电容器组的参数在线辨识方法 | |
CN111141956B (zh) | 一种基于微积分方程组单相变压器短路参数在线监测方法 | |
CN111123162A (zh) | 一种基于基波正序分量三相变压器短路参数在线监测方法 | |
CN104375025A (zh) | 中性点不接地的10kV系统中铁磁谐振的诊断方法 | |
CN105242135B (zh) | 三相电力变压器参数在线辨识装置及其实现算法 | |
CN103558481A (zh) | 同塔并架多回直流线路耦合作用引发换相失败的判别方法 | |
CN103954885A (zh) | 基于分布参数的双回线故障单端定位系统及定位方法 | |
CN104699996A (zh) | 基于双馈风电机组等值的电力系统非对称短路工频电气量计算方法 | |
CN102508098A (zh) | 用于交直流混联电网的快速零序方向元件判别方法 | |
EP3519837B1 (en) | Method to estimate an equivalent grid impedance for power generators | |
CN102621446B (zh) | 变压器三角形绕组中环流的测量方法 | |
CN105095590A (zh) | 一种基于三序等值阻抗的机电暂态仿真系统的建模方法 | |
CN111044828A (zh) | 基于正、负序方程组的三相变压器绕组参数在线监测方法 | |
Yue et al. | Coupling coefficient and load estimation for wireless power transfer systems with transmitter side input current | |
CN101866376A (zh) | Pt铁磁谐振的判别方法 | |
CN103560494A (zh) | 配电网保护整定的短路电流获取方法 | |
CN103308821A (zh) | 基于改进型锁相环的间歇性电弧接地故障选线方法 | |
CN103217585A (zh) | 一种t型输电线路正序电容测量方法 | |
Liu et al. | A selection method of mutual inductance identification models based on sensitivity analysis for wireless electric vehicles charging | |
CN103779875B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |