一种油浸式变压器热点温度监测方法
技术领域
本发明涉及电气绝缘在线检测与故障诊断领域,特别是一种油浸式变压器热点温度监测方法。
背景技术
油浸式变压器作为电力系统的核心设备,其高效安全运行对电力行业的运营发展起着至关重要的作用。变压器的工作寿命取决于内部绝缘性能,而热老化是其绝缘性能劣化的主要影响因素,故变压器的工作寿命与其内部温升息息相关。目前对变压器内部温升包括顶层油温升及热点温升的确定的一般手段是根据IEEE及IEC标准中“两段式”负载导则模型进行的,但随着电力负荷的增长,变压器时常面临着过载运行的情况,这将导致标准所推荐的绕组指数值用于热点温度计算时具有一定偏差。
发明内容
本发明的目的是提出一种油浸式变压器热点温度监测方法,该方法考虑了变压器绕组指数在多倍过负载下的计算修正值,使得其更适用于具有宽泛负载区间的油浸式变压器热点温度的计算。
实现本发明目的的技术方案如下:
一种油浸式变压器热点温度监测方法,包括
步骤1:确定修正后的绕组指数m,其方法为:
(1)令所述油浸式变压器绕组的平均温度梯度g计算模型为:
式中,gr为额定负载下绕组平均温度梯度;m为绕组指数,用来描述绕组平均温度梯度随负载的变化趋势;Ipu为负载系数,为实际负载电流I与额定电流IR的比值;Tw为绕组平均温度,Ttop为顶层油温,Tbom为底层油温;
(2)利用双对数线性回归模型对绕组指数m进行回归估计,如下:
ln(g/gr)=2mln(Ipu)
式中,令Y=ln(g/gr),X=ln(Ipu);
利用最小二乘法辨识参数m,即:
式中,j为样本个数,i为标号;
上述方法中,所述绕组平均温度Tw,顶层油温Ttop和底层油温Tbom的求取方法,包括(1)获取所述油浸式变压器结构及物性参数,包括冷热芯高位差Δh、绕组高度hw、散热器高度hr、绕组竖直油道热工水力直径Dw、散热器油道热工水力直径Dr、油比热容coil、空气比热容cair,油密度ρoil、空气密度ρair、油热膨胀系数βoil、绕组沿程阻力系数fw、散热器沿程阻力系数fr、散热器总传热系数U、绕组表面传热系数uw、绕组区域流通面积Aw、散热器流通面积Ar、散热器有效散热面积AR,绕组与油流圆周方向接触表面积As、温差指数λ、环境温度Tamb;
(2)联立以下公式及条件,利用Newton-Raphson法迭代求解绕组区域油流体积流量Gw、散热器油流体积流量Gr、顶层油温Ttop、底层油温Tbom和绕组平均温度Tw;所述油浸式变压器内部为单循环,Gw=Gr;
1)稳态条件下绕组竖直油道内的热浮升力与流体阻力达到平衡,如下:
式中,ga为重力加速度,S为冷却循环面积,
S=Δh(Ttop-Tbom)+hr[Ttop-Tbom-ΔTlm-0.5(Ttop-Tamb)],
其中,ΔTlm为散热器中油流与外界空气的对数平均温差,
2)稳态时绕组产生负载损耗Qw等于其周围油流吸收的热量,即
Qw=ρoilcoilGw(Ttop-Tbom);
3)绕组周围油流吸收的热量等于线圈沿其表面向外对流扩散的热量,即
4)油循环达到最终稳态,散热器中油相对空气温升为
Toil-Tair=Ce-λh,
式中,Toil、Tair分别为沿散热器高度方向h的油流温度及空气温度;C为底部油温和环境温度的差值;
绕组负载损耗Qw产生的热量将通过油流全部传输至外界空气,即
Qw=UAR(Toil-Tair)=UARCe-λh;
其中,Qw,R为所述油浸式变压器的额定负载损耗;
5)选定标准大气压下油浸式变压器运行时的环境温度作为参考温度;
步骤2:监测实时环境温度θamb和绕组实际负载电流I;
步骤3:实时监测修正绕组指数后的热点温度θhs,如下:
式中,n为油指数;R为额定负载损耗与空载损耗之比;Ipu为负载系数,为实际负载电流I与额定电流IR的比值;Δθoil,R为额定负载下的顶层油温升,由所述油浸式变压器出厂时确定;τoil,R为油时间常数;m为步骤1确定的修正后的绕组指数;Δθhs,R额定负载下的热点温升,由所述油浸式变压器出厂时确定;θoil为修正绕组指数m后的顶层油温;τw,R为绕组时间常数。
进一步地,
所述油时间常数τoil,R根据所述油浸式变压器的类型取值:配电变压器,油时间常数τoil,R取值为210;中型和大型电力变压器,还根据其散热类型为ONAN、ONAF或OF,油时间常数τoil,R分别取值为210、150或90;
所述绕组时间常数τw,R根据所述油浸式变压器的类型取值:配电变压器,绕组时间常数τw,R取值为4;中型和大型电力变压器,还根据其散热类型为ONAN、ONAF或OF,绕组时间常数τw,R分别取值为10、10或7;
所述油指数n根据所述油浸式变压器的散热类型为ONAN、ONAF或OF,分别取值为0.8、0.8或1。
本发明的有益效果在于,借助基于动量守恒及能量守恒的绕组区域热模型、双对数线性回归模型对反映负载率与绕组温升关系的绕组指数进行解析计算,得到其多倍负荷下的绕组指数修正值,并将其引入“两段式”温升模型进行热点温度的计算。相比于标准所推荐的参考值,该方法提升了多倍负荷下热点温度的计算精度。具有以下优点:
1)获取油浸式变压器绕组平均温升数据时,主要基于传热学及动量守恒下的绕组区域热模型进行计算,无负载系数的限制,适用于时常面临过载运行的油浸式变压器热点温度监测时绕组指数的修正计算;
2)该方法中绕组指数的计算可以选择多倍负载区间,考虑了时常发生过载运行的油浸式变压器的特殊工况,对热点温度的监测计算更准确且全面;
3)该方法可以用于不同结构的油浸式变压器的热点温度监测计算,具有普适性。
附图说明
图1为本发明的热点温度监测流程图。
具体实施方式
以下对本发明进行进一步说明。
第一步、获取已知的油浸式变压器相关结构及物性参数:
冷热芯高位差Δh、绕组高度hw、散热器高度hr、绕组竖直油道热工水力直径Dw、散热器油道热工水力直径Dr、油比热容coil、空气比热容cair,油密度ρoil、空气密度ρair、油热膨胀系数βoil、绕组沿程阻力系数fw、散热器沿程阻力系数fr、散热器总传热系数U、绕组表面传热系数uw、绕组区域流通面积Aw、散热器流通面积Ar、散热器有效散热面积AR,绕组与油流圆周方向接触表面积As、温差指数λ、环境温度Tamb。
第二步、求解待定未知量,包括:绕组区域油流体积流量Gw、散热器油流体积流量Gr、顶层油温Ttop、底层油温Tbom、绕组平均温度Tw,如下:
1)稳态条件下绕组竖直油道内的热浮升力与流体阻力达到平衡,如式所示:
式中,ga为重力加速度,S为冷却循环面积,定义如下:
S=Δh(Ttop-Tbom)+hr[Ttop-Tbom-ΔTlm-0.5(Ttop-Tamb)] (2)
式中,ΔTlm为散热器中油流与外界空气的对数平均温差,定义为:
当变压器内部为单循环时,关系式Gw=Gr成立;
2)稳态时绕组产生负载损耗Qw等于其周围油流油流吸收的热量,即:
Qw=ρoilcoilGw(Ttop-Tbom) (4)
3)绕组周围油流吸收的热量等于线圈沿其表面向外对流扩散的热量:
4)油循环达到最终稳态,散热器中油相对空气温升:
Toil-Tair=Ce-λh (6)
式中,Toil、Tair分别为沿散热器高度方向h的油流温度及空气温度;C为底部油温和环境温度差值;
绕组损耗产生的热量将通过油流全部传输至外界空气:
Qw=UAR(Toil-Tair)=UARCe-λh (7)
5)选定标准大气压下油浸式变压器运行时的环境温度作为参考温度;
综上,结合第一步已知参数,联立第1)至第5)中的公式及条件,利用Newton-Raphson法(牛顿-拉夫逊方法)迭代求解未知量:Gw、Gr、Ttop、Tbom、Tw;同时,由于绕组损耗与额定负载损耗Qw,R及负载系数Ipu存在以下关系:
因此,可根据负载系数Ipu计算出对应的负载损耗,进而求解出该负载下的上述待求解参量。
第三步、列出绕组平均温度梯度g计算模型:
式中,gr为额定负载下绕组平均温度梯度;m为绕组指数,用来描述绕组平均温度梯度随负载的变化趋势;
利用双对数线性回归模型对绕组指数m进行回归估计,如下式所示:
ln(g/gr)=2mln(Ipu) (10)
式中,令Y=ln(g/gr),X=ln(Ipu)。利用最小二乘法辨识参数m,即:
式中,j为样本个数,i为标号;
第四步、将第二步所计算的绕组平均温度数据与负载系数代入第三步中的回归模型,求出绕组指数m。
第五步、计算热点温度:
(1)、由安装在油浸式变压器外部的温度传感器(距离变压器5m及以上)将环境温度θamb及由电流传感器测得的绕组负载电流I上传至上位机;
(2)、将第二步中修正的绕组指数m代入“两段式”热点温度计算模型,如下式:
式中,θoil为顶层油温,θhs为经修正绕组指数后的热点温度;θamb是环境温度;其中Ipu为上位机中实际负载电流I与额定电流IR的比值;τoil,R及τw,R分别为油时间常数和绕组时间常数;n为油指数;
(3)将步骤(1)中监测所得的随时间变化的环境温度θamb、负载系数Ipu带入步骤(2)中的计算模型,即可求得热点温度θhs。
上述方法中,
油时间常数τoil,R根据所述油浸式变压器的类型取值:配电变压器,油时间常数τoil,R取值为210;中型和大型电力变压器,还根据其散热类型为ONAN、ONAF或OF,油时间常数τoil,R分别取值为210、150或90;
绕组时间常数τw,R根据所述油浸式变压器的类型取值:配电变压器,绕组时间常数τw,R取值为4;中型和大型电力变压器,还根据其散热类型为ONAN、ONAF或OF,绕组时间常数τw,R分别取值为10、10或7;
油指数n根据所述油浸式变压器的散热类型为ONAN、ONAF或OF,分别取值为0.8、0.8或1。