CN108075469B - 一种配电变压器零序阻抗计算方法 - Google Patents
一种配电变压器零序阻抗计算方法 Download PDFInfo
- Publication number
- CN108075469B CN108075469B CN201711382991.4A CN201711382991A CN108075469B CN 108075469 B CN108075469 B CN 108075469B CN 201711382991 A CN201711382991 A CN 201711382991A CN 108075469 B CN108075469 B CN 108075469B
- Authority
- CN
- China
- Prior art keywords
- phase
- voltage
- zero sequence
- distribution transformer
- impedance
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 230000002068 genetic effect Effects 0.000 claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 4
- 230000007935 neutral effect Effects 0.000 claims description 9
- 238000013139 quantization Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 10
- 238000012360 testing method Methods 0.000 abstract description 6
- 238000002847 impedance measurement Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
一种配电变压器零序阻抗计算方法,所述方法在三相负荷不平衡时,通过台区采集装置获取配电变压器运行电气量,建立零序阻抗待求量的电压回路方程组,将电压回路方程组求值问题转化成目标函数优化问题,应用遗传算法计算出配电变压器零序阻抗。本发明针对配电变压器出厂试验一般不做零序阻抗试验的情况,在利用台区配电变压器已有运行数据的基础上,建立含零序阻抗待求量非线性方程组,考虑非线性方程组难以用一般数值方法计算,将非线性方程组求解问题转化成目标优化问题,应用遗传算法对目标优化函数进行求解,最终求出不同零序电流下配电变压器零序阻抗。
Description
技术领域
本发明涉及一种配电变压器零序阻抗计算方法,属变压器技术领域。
背景技术
配电变压器零序阻抗会对0.4kV低压配电网运行带来影响。低压配电网三相负荷均衡运行时,三相电压不会产生偏移,然而实际情况中三相负荷往往是不平衡的,加上配电变压器零序阻抗的存在,使得中性点电压出现偏移,至少一相负荷端电压升高,用户家用电器在非正常电压下运行容易烧毁,零序阻抗值大小直接影响到中性点电压偏移程度。
目前,配电变压器短路阻抗参数通过出厂试验获得,零序阻抗一般不做出厂试验,导致实际运行中无法根据零序阻护有效掌握中性点电压偏移程度,相关测量零序阻抗方法需在投运前进行测量,无法根据配变运行工况进行在线计算,另外根据仿真软件等方法来计算零序阻抗大小,计算过程复杂,需考虑变压器内部结构,同样无法满足实际测量需求,实用性不强。
台区配电变压器实际运行工况复杂,柱上配电变压器常年在户外运行,配电变压器零序阻抗有可能发生变化,如何随时掌握配电变压器零序阻抗情况对低压配电网运行具有重要的意义。本文提出了一种配电变压器零序阻抗计算方法,在三相负荷不平衡时,通过台区采集装置获取配电变压器运行电气量,建立含零序阻抗待求量的电压回路方程组,将电压回路方程组求值问题转化成目标函数优化问题,应用遗传算法计算出配电变压器零序阻抗。
发明内容
本发明的目的是,为了解决配电变压器零序阻抗一般不做出厂试验,相关测量零序阻抗方法需在投运前进行测量,又无法在线计算的问题,本发明公开一种配电变压器零序阻抗计算方法。
实现本发明的技术方案是,在三相负荷不平衡时,通过台区采集装置获取配电变压器运行电气量,建立含零序阻抗待求量的电压回路方程组,将电压回路方程组求值问题转化成目标函数优化问题,应用遗传算法计算出配电变压器零序阻抗。
所述零序阻抗待求量的电压回路方程组即零序阻抗复数方程组,包括以下表达式:
其中,Ua、Ub、Uc分别为低压综合配电柜采集侧各相电压;Ia、Ib、Ic分别为低压综合配电柜采集侧各相电流;Ra、Rb、Rc分别为配电变压器各相短路阻抗的线圈直流电阻;Xa、Xb、Xc分别为配电变压器各相短路阻抗的线圈漏抗;Rn为零序阻抗中的电阻分量;Xn为零序阻抗中的电抗分量。
所述零序阻抗复数方程组标量化处理为:
In(re)=-(Ia*cos(αa-θa)+Ib*cos(αb-θb)+Ic*cos(αc-θc))
In(im)=-(Ia*sin(αa-θa)+Ib*sin(αb-θb)+Ic*sin(αc-θc))
相应的中性点电压实部和虚部分量分别为:
Un(re)=Rn*In(re)-Xn*In(im)
Un(im)=Rn*In(im)+Xn*In(re)
A相阻抗电压实部和虚部分量分别为:
A相电压回路方程实部和虚部分别为:
Ea*cos(0)=Ua*cos(αa)+Uat(re)+Un(re);
Ea*sind(0)=Ua*sin(αa)+Uat(im)+Un(im);
同理求出B、C相电压路方程实部和虚部等式:
Eb*cos(-120)=Ub*cos(αb)+Ubt(re)+Un(re);
Eb*sin(-120)=Ub*sin(αb)+Ubt(im)+Un(im);
Ec*cos(120)=Uc*cos(αc)+Uct(re)+Un(re);
Ec*sin(120)=Uc*sin(αc)+Uct(im)+Un(im)。
所述应用遗传算法计算出配电变压器零序阻抗为:
由A、B、C相电压路方程实部和虚部等式,得到适应度函数,将求解零序阻抗问题转换成适应度函数最优问题:
本发明的有益效果是,本发明提出的一种配电变压器零序阻抗计算方法,针对配电变压器出厂试验一般不做零序阻抗试验的情况,在利用台区配电变压器已有运行数据的基础上,建立含零序阻抗待求量非线性方程组,考虑非线性方程组难以用一般数值方法计算,将非线性方程组求解问题转化成目标优化问题,应用遗传算法对目标优化函数进行求解,最终求出不同零序电流下配电变压器零序阻抗。本发明计算方法,对两种不同联结组别配电变压器零序阻抗计算结果与实测值相比,误差在±10%之内,计算结果可满足运行实际需求。
本发明计算方法可实现配电变压器零序阻抗在线计算,而不需要通过试验方法获取,实用性强。
附图说明
图1为本发明实施例低压台区三相四线等值电路图;
图2为本发明遗传算法流程图。
具体实施方式
本发明的具体实施方式如图1和图2所示。
(1)建立零序阻抗复数方程组
配电变压器短路阻抗由线圈直流电阻和漏抗组成,分别计为Ra、Rb、Rc、Xa、Xb、Xc,零序阻抗中的电阻分量和电抗分量分别计为Rn、Xn,低压综合配电柜采集侧电压、电流分别计为Ua、Ub、Uc、Ia、Ib、Ic。
根据图1建立含零序阻抗复数方程:
(2)三相测量点电压与电压源夹角计算
设三相测量点电压与电压源的夹角分别αa、αb、αc,夹角关系为:
αa+αb+αc=0……………………………………………………6
(3)三相测量点电压与电流之间夹角计算
三相测量点电压与电流之间的夹角分别为θa、θb、θc,可通过低压综合配电柜采集到的三相有功、无功值计算出三相测量点电压与电流之间的夹角:
θa=artan(Qa/Pa)
θb=artan(Qb/Pb)
θc=artan(Qc/Pc)…………………………………………………7
(4)短路电阻与短路电抗计算
配电变压器短路电阻Rk和短路电抗Xk均可通过阻抗电压Uk%、短路损耗Pf、额定电流Ie、额定电压Ue4个出厂参数计算得出,并假定三相短路电阻和短路电阻值相等,计算公式如下:
Zk=Uk%*Ue/Ie…………………………………………………8
Rk=Pf/(3*Ie*Ie)…………………………………………………9
Ra=Rb=Rc=Rk…………………………………………………11
Xa=Xb=Xc=Xk…………………………………………………12
(5)零序阻抗方程组标量化
根据公式(4),中性点电流分解成实部和虚部分量:
In(re)=-(Ia*cos(αa-θa)+Ib*cos(αb-θb)+Ic*cos(αc-θc))………………13
In(im)=-(Ia*sin(αa-θa)+Ib*sin(αb-θb)+Ic*sin(αc-θc))………………14
相应的中性点电压实部和虚部分量分别为:
Un(re)=Rn*In(re)-Xn*In(im)………………………………………15
Un(im)=Rn*In(im)+Xn*In(re)………………………………………16
A相阻抗电压实部和虚部分量分别为:
A相电压回路方程实部和虚部分别为:
Ea*cos(0)=Ua*cos(αa)+Uat(re)+Un(re)………………………………19
Ea*sind(0)=Ua*sin(αa)+Uat(im)+Un(im)………………………………20
同理求出B、C相电压路方程实部和虚部等式:
Eb*cos(-120)=Ub*cos(αb)+Ubt(re)+Un(re)……………………………21
Eb*sin(-120)=Ub*sin(αb)+Ubt(im)+Un(im)……………………………22
Ec*cos(120)=Uc*cos(αc)+Uct(re)+Un(re)………………………………23
Ec*sin(120)=Uc*sin(αc)+Uct(im)+Un(im)………………………………24
(6)遗传算法求解非线性方程组
将公式19、20、21、22、23、24组成的方程组求解问题转化成如下优化问题:
将公式19、20、21、22、23、24右侧全部移至左侧,分别取绝对值,再将6个公式相加,便得适应度函数,将求解零序阻抗问题转换成适应度函数最优问题,通过遗传优化算法便可求解。
Claims (1)
1.一种配电变压器零序阻抗计算方法,其特征在于,所述方法在三相负荷不平衡时,通过台区采集装置获取配电变压器运行电气量,建立零序阻抗待求量的电压回路方程组,将电压回路方程组求值问题转化成目标函数优化问题,应用遗传算法计算出配电变压器零序阻抗;
所述零序阻抗待求量的电压回路方程组即零序阻抗复数方程组,包括以下表达式:
其中,分别为配电变压器低压侧A相、B相、C相三相电势; 分别为低压综合配电柜采集侧各相电压;分别为低压综合配电柜采集侧各相电流;Ra、Rb、Rc分别为配电变压器各相短路阻抗的线圈直流电阻;Xa、Xb、Xc分别为配电变压器各相短路阻抗的线圈漏抗;Rn为零序阻抗中的电阻分量;Xn为零序阻抗中的电抗分量;为配电变压器零序电流,为配电变压器零序电压;分别为配电变压器低压侧A相、B相、C相三相电势;Un为零序电压;
所述零序阻抗复数方程组标量化处理为:
In(re)=-(Ia*cos(αa-θa)+Ib*cos(αb-θb)+Ic*cos(αc-θc))
In(im)=-(Ia*sin(αa-θa)+Ib*sin(αb-θb)+Ic*sin(αc-θc))
其中,αa、αb、αc分别为A相、B相、C相三相电势与相应地A相、B相、C相测量点电压之间的夹角;θa、θb、θc分别为A相、B相、C相三相测量点电压与相应地A相、B相、C相电流之间的夹角;Ia、Ib、Ic分别表示三相电流的模;In(re)表示零序电流的实部,In(im)表示零序电流的虚部;
相应的中性点电压实部和虚部分量分别为:
Un(re)=Rn*In(re)-Xn*In(im)
Un(im)=Rn*In(im)+Xn*In(re)
其中,Un(re)表示零序电压的实部,Un(im)表示零序电压的虚部;
A相阻抗电压实部和虚部分量分别为:
其中,Uat(re)表示A相阻抗电压的实部;Uat(im)表示A相阻抗电压的虚部;
A相电压回路方程实部和虚部分别为:
Ea*cos(0°)=Ua*cos(αa)+Uat(re)+Un(re);
Ea*sin(0°)=Ua*sin(αa)+Uat(im)+Un(im);
Eb*cos(-120°)=Ub*cos(αb)+Ubt(re)+Un(re);
Eb*sin(-120°)=Ub*sin(αb)+Ubt(im)+Un(im);
Ec*cos(120°)=Uc*cos(αc)+Uct(re)+Un(re);
Ec*sin(120°)=Uc*sin(αc)+Uct(im)+Un(im);
其中,Ubt(re)、Uct(re)表示B相、C相阻抗电压的实部,Ubt(im)、Uct(im)表示B相、C相阻抗电压的虚部;Ub、Uc表示B相、C相采集点电压的模,Eb、Ec表示B相、C相电势的模;
所述应用遗传算法计算出配电变压器零序阻抗为:
由A、B、C相电压回路方程实部和虚部等式,得到适应度函数,假定Ea=Eb=Ec=E,已知Ua、Ub、Uc、Ia、Ib、Ic、θa、θb、θc、Ra、Rb、Rc、Xa、Xb、Xc,求解六个未知数E、αa、αb、αc、Rn、Xn,将求解零序阻抗问题转换成适应度函数最优问题:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711382991.4A CN108075469B (zh) | 2017-12-20 | 2017-12-20 | 一种配电变压器零序阻抗计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711382991.4A CN108075469B (zh) | 2017-12-20 | 2017-12-20 | 一种配电变压器零序阻抗计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108075469A CN108075469A (zh) | 2018-05-25 |
CN108075469B true CN108075469B (zh) | 2021-07-02 |
Family
ID=62158584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711382991.4A Active CN108075469B (zh) | 2017-12-20 | 2017-12-20 | 一种配电变压器零序阻抗计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108075469B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108776869B (zh) * | 2018-06-07 | 2022-03-04 | 国网河南省电力公司电力科学研究院 | 基于智能电表大数据的台区变压器三相不平衡治理方法 |
CN109192472A (zh) * | 2018-09-12 | 2019-01-11 | 合肥多元节电科技有限公司 | 一种新式调整变压器阻抗的方法 |
CN111242459A (zh) * | 2020-01-07 | 2020-06-05 | 中国南方电网有限责任公司 | 全网设备参数异常值辨识方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101419254A (zh) * | 2008-12-17 | 2009-04-29 | 河南电力试验研究院 | 一种特高压输电线路参数测量***及方法 |
CN102361321A (zh) * | 2011-09-05 | 2012-02-22 | 四川省电力公司资阳公司 | 110/220kV变压器中性点经小电抗接地 |
CN104166767A (zh) * | 2014-08-18 | 2014-11-26 | 国家电网公司 | 一种将变压器零序试验数据转换为计算用零序参数的方法 |
US8964423B2 (en) * | 2011-08-17 | 2015-02-24 | Honeywell International Inc. | Low weight 3-phase 5-output wire power conversion system for micro-grid |
CN206657058U (zh) * | 2017-04-11 | 2017-11-21 | 广东电网有限责任公司惠州供电局 | 一种变压器短路阻抗在线检测*** |
-
2017
- 2017-12-20 CN CN201711382991.4A patent/CN108075469B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101419254A (zh) * | 2008-12-17 | 2009-04-29 | 河南电力试验研究院 | 一种特高压输电线路参数测量***及方法 |
US8964423B2 (en) * | 2011-08-17 | 2015-02-24 | Honeywell International Inc. | Low weight 3-phase 5-output wire power conversion system for micro-grid |
CN102361321A (zh) * | 2011-09-05 | 2012-02-22 | 四川省电力公司资阳公司 | 110/220kV变压器中性点经小电抗接地 |
CN104166767A (zh) * | 2014-08-18 | 2014-11-26 | 国家电网公司 | 一种将变压器零序试验数据转换为计算用零序参数的方法 |
CN206657058U (zh) * | 2017-04-11 | 2017-11-21 | 广东电网有限责任公司惠州供电局 | 一种变压器短路阻抗在线检测*** |
Non-Patent Citations (2)
Title |
---|
变压器的零序阻抗及测量;赵旺初;《华北电力技术》;19951231(第4期);第53-56页 * |
基于阻抗模型故障特征匹配技术的配电网单相接地故障测距研究;康忠健 等;《电网与清洁能源》;20150831;第31卷(第8期);第40-46页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108075469A (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108075469B (zh) | 一种配电变压器零序阻抗计算方法 | |
CN102269787B (zh) | 电子式电压互感器谐波特性检测方法及其检测装置 | |
CN104316894B (zh) | 一种电能表实际运行环境的仿真校验方法 | |
CN107248736B (zh) | 一种配电网线路正序参数的在线辨识方法 | |
CN105720594B (zh) | 纯有功三相不平衡负荷补偿容量的补偿方法 | |
CN103293504B (zh) | 现场电流互感器综合检测装置的测试方法 | |
CN107359612A (zh) | 一种电能质量对配电网能耗影响的综合评估方法 | |
CN102435858A (zh) | 变压器短路损耗和开路损耗的在线测量方法及*** | |
CN108303606A (zh) | 一种配电变压器容量在线评估方法 | |
Islam et al. | Time-frequency-based instantaneous power components for transient disturbances according to IEEE standard 1459 | |
CN103870703A (zh) | 一种基于戴维南等值参数跟踪的动态短路比计算方法 | |
CN103592528A (zh) | 一种基于动态轨迹灵敏度的光伏逆变器模型参数辨识方法 | |
Karmakar et al. | Machine learning approach to estimation of internal parameters of a single phase transformer | |
CN111141995A (zh) | 一种基于比幅原理的线路双端稳态量测距方法和*** | |
Orts-Grau et al. | Discussion on useless active and reactive powers contained in the IEEE standard 1459 | |
WO2015192273A1 (zh) | 一种三相交流电弧炉自适应建模装置及其仿真算法 | |
CN107167707A (zh) | 一种基于参数未知的同杆双回线故障测距方法及装置 | |
Han et al. | A new assessment for the total harmonic contributions at the point of common coupling | |
CN108181600A (zh) | 电容式电压互感器测试装置 | |
CN102004191A (zh) | 消弧线圈并列运行***中计算线路电容的方法 | |
CN203250017U (zh) | 现场电流互感器综合检测装置 | |
CN203981735U (zh) | 多变比电流互感器自动调节装置 | |
CN113054747B (zh) | 一种配电台区非线性负荷的非侵入式检测方法 | |
CN204855583U (zh) | 一种精密大电流的暂态升流发生装置 | |
CN104934969B (zh) | 一种电力线路参数的计算方法 |
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 |