CN111181160B - 一种10kV中压配电网理论线损计算方法及系统 - Google Patents

一种10kV中压配电网理论线损计算方法及系统 Download PDF

Info

Publication number
CN111181160B
CN111181160B CN202010085895.9A CN202010085895A CN111181160B CN 111181160 B CN111181160 B CN 111181160B CN 202010085895 A CN202010085895 A CN 202010085895A CN 111181160 B CN111181160 B CN 111181160B
Authority
CN
China
Prior art keywords
equivalent
voltage
medium
line loss
voltage distribution
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
Application number
CN202010085895.9A
Other languages
English (en)
Other versions
CN111181160A (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010085895.9A priority Critical patent/CN111181160B/zh
Publication of CN111181160A publication Critical patent/CN111181160A/zh
Application granted granted Critical
Publication of CN111181160B publication Critical patent/CN111181160B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Abstract

一种10kV中压配电网理论线损计算方法及系统,所述方法包括中压线路变压器出口端电气数据采集,计算每个采集点的等效电压值和等效电流值,对等效电压序列和等效电流序列进行多项式拟合,并对拟合的多项式求导,得到不同等效电流值下的等效电阻;最后根据理论线损计算公式得到每个采集时刻的理论线损功率,进而求得不同统计周期下的中压配电理论线损值及理论线损率。本发明方法,相比于常规的等效电阻法计算方法,无需考虑中压配电线路的拓扑关系、线路参数、配电变压器参数,具有较强的实用性。本发明所涉及的等效电阻计算公式,充分考虑单个采集点的随机性,通过多项式拟合得到实用合理的等效电阻,提高了中压配电线路理论线损值的准确率。

Description

一种10kV中压配电网理论线损计算方法及系统
技术领域
本发明涉及一种10kV中压配电网理论线损计算方法及系统,属配用电技术领域。
背景技术
在中压配电网电能传输过程中,流经中压配电网一次各类设备将产生相应的电能损耗,如何对中压配电网理论线损进行快速准确的计算,及时发现真实高线损配电线路,对供电企业科学精准施策降损具有重要的意义。
目前,中压配电网理论线损计算方法,主要有基于配电变压器容量等值电阻法、基于配电变压器电量等值电阻法、前推回代潮流法。三种中压配电理论线损计算方法均要考虑到配电网拓扑关系、线路参数等因素,但中压线路涉及配电变压器数量多、运行状态多变、设备异动频繁,采用前三种方法无法大规模地、准确有效地对中压线路理论线损计算。
发明内容
本发明的目的是,为了有效快速计算中压配电网理论线损值,无需考虑中压配电拓扑关系、线路参数等因素,本发明提出一种10kV中压配电网理论线损计算方法及系统。
实现本发明的技术方案如下:一种10kV中压配电网理论线损计算方法,包括中压线路变压器出口端电气数据采集,并对其进行预处理;所述方法计算每个采集点的等效电压值和等效电流值,对等效电压序列和等效电流序列进行多项式拟合,并对拟合的多项式求导,得到不同等效电流值下的等效电阻;最后根据理论线损计算公式得到每个采集时刻的理论线损功率,进而求得不同统计周期下的中压配电理论线损值及理论线损率。
所述等效电压值和等效电流值计算式如下:
等效电流值:ieqt=(ict-iat)/ict (1)
其中:iat、ict分别为t时刻的中压配电线路a、c两相出口电流值,ieqt为t时该中压配电线路的等效电流值;
等效电压值:ueqt=uabt/iat-ubct×iat/(ict×ict) (2)
其中:uabt、ubct分别为t时该中压配电线路的ab、bc两相出口线电压值,ueqt为t时该中压配电线路的等效电压值;
经式(1)和式(2)计算,分别得到相应等效电流序列ieq和等效电压序列ueq
所述多项式拟合,包括对等效电流序列ieq和等效电压序列ueq经采用3次多项式拟合,拟合公式下式:
y=a0+a1x+a2x2+a3x3 (3)
其中:x、y分别为自变量、因变量;a0、a1、a2、a3分别拟合参数。
所述拟合参数计算方法如下:
给定等效电流序列ieq和等效电压序列ueq,根据下式得到误差函数Q:
Figure GDA0002446005210000021
/>
求拟合参数使得误差函数最小,拟合参数计算公式如下:
Figure GDA0002446005210000022
其中,
Figure GDA0002446005210000031
AT为A的转置矩阵;
进而便可得到3次多项式拟合公式。
所述等效电阻计算公式如下:
对3次多项式拟合公式进行求导,便可得到等效电阻Reqt
Figure GDA0002446005210000032
其中,a1、a2、a3为根据式(5)得到的拟合参数;Reqt为t时该中压配电线路的等效电阻值。
所述理论线损计算如下:
根据每个采集点的Reqt和出口电流值,计算得到该采集时刻的理论线损值Pxst,如下式:
Pxst=3Reqt((iat+ict)/2)2 (7)
其中:Pxst为t时刻理论线损值。
所述理论线损率Pxst%计算如下:
Figure GDA0002446005210000033
其中,Pt为采集装置得到t时刻中压配电压出口有功功率;n为序列采集个数。
一种10kV中压配电网理论线损计算系统,包括变电站出线电气采集装置、计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时,实现以下步骤:
(1)对变电站出线电气采集装置采集的电压、电流数据进行预处理;
(2)计算每个采集点的等效电压值和等效电流值;
(3)对等效电压序列和等效电流序列进行多项式拟合;
(4)对多项式公式进行求导,得到不同等效电流值下的等效电阻;
(5)计算理论线损;
(6)计算不同时间尺度下的理论线损率。
所述采集装置采集中压配电线路变电站出口侧AB、BC电压值和A、C相电流值;采样间隔为15分钟或30分钟;电压互感器采用电压互感器V/V接法;电流互感器采用不完全星型接法。
本发明的有益效果在于,本发明所涉及的等效电阻计算公式,相比于常规的等效电阻法计算方法,无需考虑中压配电线路的拓扑关系、线路参数、配电变压器参数,具有较强的实用性。本发明所涉及的等效电阻计算公式,充分考虑单个采集点的随机性,通过多项式拟合得到实用合理的等效电阻,提高了中压配电线路理论线损值的准确率。
本发明提出的一种中压配电网理论线损计算方法及系统,充分利用变电站已有电气采集装置,无需新增硬件,通过计算机可执行程序代码便可实现中压配电网理论线损计算。
附图说明
图1为多时间尺度中压配电网理论线损计算方法流程图;
图2为电压互感器V/V接法。
具体实施方式
如图1所示,为本发明多时间尺度中压配电网理论线损计算方法流程图。
本发明的具体实施步骤如下:
(1)数据采集
本实施例利用变电站出线电气采集装置采集中压配电线路变电站出口侧AB、BC电压和A、C相电流值,采样间隔一般为15分钟或30分钟,电压互感器采用电压互感器V/V接法,电流互感器采用不完全星型接法,电压互感器标准接线如图2所示。
(2)计算等效电流和等效电压
将采集到的中压配电线路电压、电流数据序列按时间对齐,若某一个采集时刻至少有一项存在缺失值、零值或异常值(取电流值大于1500A或电压值大于20000A以上判断为异常值),则删除该采集时刻所有采集值,每个采集点等效电流和等效电压值计算公式如式(1)、式(2)所示:
1)等效电流值
ieqt=(ict-iat)/ict (1)
其中:iat、ict分别为t时刻的中压配电线路a、c两相出口电流值,ieqt为t时该中压配电线路的等效电压值。
2)等效电压值:
ueqt=uabt/iat-ubct×iat/(ict×ict) (2)
其中:uabt、ubct分别为t时该中压配电线路的ab、bc两相出口线电压值,ueqt为t时该中压配电线路的等效电压值。
经式(1)、式(2)计算分别得到相应等效电流序列ieq和等效电压序列ueq
(3)多项公式拟合
对等效电流序列ieq和等效电压序列ueq经采用3次多项式拟合,拟合公式如式(3)所示:
y=a0+a1x+a2x2+a3x3 (3)
其中:x、y分别为自变量、因变量,a0、a1、a2、a3分别拟合参数。
给定等效电流序列ieq和等效电压序列ueq,根据式(4)得到误差函数:
Figure GDA0002446005210000061
求拟合参数使得误差函数最小,拟合参数计算公式如式(5)所示:
Figure GDA0002446005210000062
其中:
Figure GDA0002446005210000063
AT为A的转置矩阵。
进而便可得到3次多项式拟合公式,对3次多项式拟合公式进行求导,便可得到等效电阻计算公式,如式(6)所示:
Figure GDA0002446005210000064
其中,a1、a2、a3根据式(5)得到的拟合参数,Reqt为t时该中压配电线路的等效电阻值。
(4)理论线损计算
根据每个采集点的Reqt和出口电流值,便可计算得到该采集时刻的理论线损值Pxst,如式(7)所示:
Pxst=3Reqt((iat+ict)/2)2 (7)
进而可求得统计周期下相应地理论线损率Pxst%,如式(8)所示:
Figure GDA0002446005210000065
其中:Pxst为公式(7)计算得到t时刻理论线损值,Pt为采集装置得到t时刻中压配电压出口有功功率,n为序列采集个数。
由于n视情况取值,进而可以求得统计同期下的中压配电网理论线损值及理论线损率。
(5)实施例验证
为验证本实施例的一种多时间尺度中压配电网理论线损计算方法的有效性,以江西某公司一条中压配电线路2019年6月理论线损计算为例,计算结果如下表所示:
Figure GDA0002446005210000071
Figure GDA0002446005210000081
实际理论线损率在拓扑关系、线路参数、配电变压器参数准确的基础上计算得到,由上表可知,本发明计算得到的理论线损率与实际理论线损率最大偏差不超过±10%,可以满足工程实用要求。

Claims (6)

1.一种10kV中压配电网理论线损计算方法,包括中压线路变压器出口端电气数据采集,并对其进行预处理,其特征在于,所述方法计算每个采集点的等效电压值和等效电流值,对等效电压序列和等效电流序列进行多项式拟合,并对拟合的多项式求导,得到不同等效电流值下的等效电阻;最后根据理论线损计算公式得到每个采集时刻的理论线损功率,进而求得不同统计周期下的中压配电理论线损值及理论线损率;
所述等效电压值和等效电流值计算式如下:
等效电流值:ieqt=(ict-iat)/ict(1)
其中:iat、ict分别为t时刻的中压配电线路a、c两相出口电流值,ieqt为t时该中压配电线路的等效电流值;
等效电压值:ueqt=uabt/iat-ubct×iat/(ict×ict)(2)
其中:uabt、ubct分别为t时该中压配电线路的ab、bc两相出口线电压值,ueqt为t时该中压配电线路的等效电压值;
经式(1)和式(2)计算,分别得到相应等效电流序列ieq和等效电压序列ueq
所述多项式拟合,包括对等效电流序列ieq和等效电压序列ueq经采用3次多项式拟合,拟合公式下式:
y=a0+a1x+a2x2+a3x3 (3)
其中:x、y分别为自变量、因变量;a0、a1、a2、a3分别为拟合参数;
所述拟合参数计算方法如下:
给定等效电流序列ieq和等效电压序列ueq,根据下式得到误差函数Q:
Figure QLYQS_1
求拟合参数使得误差函数最小,拟合参数计算公式如下:
Figure QLYQS_2
其中,
Figure QLYQS_3
AT为A的转置矩阵,其中,ieq1、ieq2、...、ieqn分别为t1、t2、...、tn时刻等效电流值,ueq1、ueq2、...、ueqn分别为t1、t2、...、tn时刻等效电压值;
进而便可得到3次多项式拟合公式,得到a1、a2、a3
2.根据权利要求1所述的一种10kV中压配电网理论线损计算方法,其特征在于,所述等效电阻计算公式如下:
对3次多项式拟合公式进行求导,便可得到等效电阻Reqt
Figure QLYQS_4
其中,a1、a2、a3为拟合参数计算公式得到的拟合参数;Reqt为t时该中压配电线路的等效电阻值。
3.根据权利要求1所述的一种10kV中压配电网理论线损计算方法,其特征在于,所述理论线损计算如下:
根据每个采集点的Reqt和出口电流值,计算得到该采集时刻的理论线损值Pxst,如下式:
Pxst=3Reqt((iat+ict)/2)2 (7)
其中:Pxst为t时刻理论线损值,Reqt为t时该中压配电线路的等效电阻值,iat、ict分别为t时刻的中压配电线路a、c两相出口电流值。
4.根据权利要求1所述的一种10kV中压配电网理论线损计算方法,其特征在于,所述理论线损率Pxst%计算如下:
Figure QLYQS_5
其中,Pt为采集装置得到t时刻中压配电压出口有功功率;n为序列采集个数。
5.实现如权利要求1~4任一项所述10kV中压配电网理论线损计算方法的系统,包括变电站出线电气采集装置、计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时,实现以下步骤:
(1)对变电站出线电气采集装置采集的电压、电流数据进行预处理;
(2)计算每个采集点的等效电压值和等效电流值;
(3)对等效电压序列和等效电流序列进行多项式拟合;
(4)对多项式公式进行求导,得到不同等效电流值下的等效电阻;
(5)计算理论线损;
(6)计算不同时间尺度下的理论线损率。
6.根据权利要求5所述10kV中压配电网理论线损计算方法的系统,其特征在于,所述采集装置采集中压配电线路变电站出口侧AB、BC电压值和A、C相电流值;采样间隔为15分钟或30分钟;电压互感器采用电压互感器V/V接法;电流互感器采用不完全星型接法。
CN202010085895.9A 2020-02-11 2020-02-11 一种10kV中压配电网理论线损计算方法及系统 Active CN111181160B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010085895.9A CN111181160B (zh) 2020-02-11 2020-02-11 一种10kV中压配电网理论线损计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010085895.9A CN111181160B (zh) 2020-02-11 2020-02-11 一种10kV中压配电网理论线损计算方法及系统

Publications (2)

Publication Number Publication Date
CN111181160A CN111181160A (zh) 2020-05-19
CN111181160B true CN111181160B (zh) 2023-03-28

Family

ID=70658254

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010085895.9A Active CN111181160B (zh) 2020-02-11 2020-02-11 一种10kV中压配电网理论线损计算方法及系统

Country Status (1)

Country Link
CN (1) CN111181160B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5036364B2 (ja) * 2007-03-29 2012-09-26 中部電力株式会社 電力系統の線路損失低減装置、電力系統及び電力系統構築方法
CN103439608B (zh) * 2013-08-29 2016-06-29 国家电网公司 降低10kV配电网线损的方法
CN109507536A (zh) * 2018-12-12 2019-03-22 国家电网有限公司 一种低压台区线损电量智能测量分析设备
CN109494735B (zh) * 2018-12-25 2020-07-07 中国电力科学研究院有限公司 一种低压配电网线损计算模型及方法

Also Published As

Publication number Publication date
CN111181160A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN102567603B (zh) 基于实测拓扑和量测数据自动生成bpa计算文件的方法
CN102495288B (zh) 基于scada量测数据的输电线路参数的辨识与估计方法
CN107257130B (zh) 基于区域量测解耦的低压配电网损耗计算方法
Ashraf et al. Performance analysis of static network reduction methods commonly used in power systems
CN107037322A (zh) 基于稳态特征的配电网小电流接地故障定位方法
CN107069741A (zh) 一种新型线性化潮流计算方法
CN107944680A (zh) 一种基于智能电能表的变电站电能量平衡监测方法
CN106054023A (zh) 一种输电线路单端测距中估计两侧系统阻抗的方法
CN104794531A (zh) 基于网损等值负荷模型的改进直流最优潮流方法
CN112039108B (zh) 高压直流输电系统单端暂态响应快速估计模型及计算方法
CN106203889A (zh) 非健全信息下基于最小二乘拟合法的配电网损耗评估方法
CN113777440A (zh) 一种光伏电站集电线路单相接地故障的定位方法
CN107462810B (zh) 一种适用于有源配电网的故障区段定位方法
CN112098888A (zh) 一种基于小电流接地选线装置的配电网系统电容电流计算方法
CN105320810B (zh) 短路电流扫描优化方法
CN113381408B (zh) 一种基于配电自动化数据的用电异常配变定位方法及装置
CN111181160B (zh) 一种10kV中压配电网理论线损计算方法及系统
CN107167707A (zh) 一种基于参数未知的同杆双回线故障测距方法及装置
CN108448609B (zh) 一种考虑dg等值容量利用小时系数的配网网损分摊方法
CN110133420A (zh) 一种用于不平衡配电网的单相接地故障快速线性计算方法
CN104573384B (zh) 一种适用于配电网在线分析的改进直流潮流算法
CN110048428B (zh) 基于概率守恒原理的电力系统概率潮流计算方法
CN111257607A (zh) 一种稽查单相电能表并联电阻窃电的方法
CN106953331A (zh) 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法
CN111130109B (zh) 一种Yyn0低压配电网理论线损计算方法及系统

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