CN109188204B - 计及小电源的配电网理论线损计算方法 - Google Patents

计及小电源的配电网理论线损计算方法 Download PDF

Info

Publication number
CN109188204B
CN109188204B CN201811043738.0A CN201811043738A CN109188204B CN 109188204 B CN109188204 B CN 109188204B CN 201811043738 A CN201811043738 A CN 201811043738A CN 109188204 B CN109188204 B CN 109188204B
Authority
CN
China
Prior art keywords
branch
current
matrix
loss
node
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.)
Expired - Fee Related
Application number
CN201811043738.0A
Other languages
English (en)
Other versions
CN109188204A (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.)
University of Jinan
Original Assignee
University of Jinan
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 University of Jinan filed Critical University of Jinan
Priority to CN201811043738.0A priority Critical patent/CN109188204B/zh
Publication of CN109188204A publication Critical patent/CN109188204A/zh
Application granted granted Critical
Publication of CN109188204B publication Critical patent/CN109188204B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

本发明公开了一种计及小电源的配电网理论线损计算方法,获取配电网各支路的阻抗,各节点的有功功率、无功功率、电压和电流数据;统计各小电源节点的用电量,采用平均电流法计算配电网的电能损耗;根据各节点的有功功率、无功功率、电压幅值和电流幅值,计算各节点的注入电流;采用搜索法形成路径互阻矩阵;基于路径互阻矩阵和各节点的注入电流,采用回路分析法得到精确的电能损耗计算模型和基于平均电流法的损耗计算模型;利用精确的电能损耗计算模型和基于平均电流法的损耗计算模型,得到平均电流法的计算误差,即损耗的修正量;计算平均电流法得到的电能损耗与损耗的修正量之和,得到配电网的理论电能损耗。

Description

计及小电源的配电网理论线损计算方法
技术领域
本发明属于电力系统技术领域,具体涉及一种计及小电源的配电网理论线损计算方法。
背景技术
配电网线损是一项综合的经济技术指标。合理准确的线损理论计算方法对降低线损、提高经济效益具有重要的现实意义。目前,配电网理论线损计算方法主要包括均方根电流法、平均电流法、损失因数法和等值电阻法等。其中,均方根电流法计算过程简单,计算精度高,但由于配电网结构复杂,参数不完善等原因导致该方法计算工作量大,计算准确性不能保证。利用平均电流与均方根电流的等效关系,工程上多采用由均方根电流法派生的平均电流法。平均电流法需假设各负荷节点功率因数和负荷曲线与馈线首端相同,忽略了沿线电压的损失,因而计算结果存在一定的误差,且形状系数不易确定;等值电阻法在计算等值电阻时采用与平均电流法相同的条件,因而也存在同样的问题;损失因数法利用负荷曲线的最大值与均方根值等效关系进行计算,所需计算数据少,但损耗因数不易计算,计算精度较低。
针对以上各方法的局限性,有学者提出相应的改进方法。潮流及估计技术利用实时测量数据对负荷数据进行修正,提高了计算精度;等值电阻法通过潮流计算改进假设条件中不精确的参数,得到改善的等值电阻法,实际上由于此类算法数据采集困难,并且会带来计算不收敛的问题,因而应用较少;改进损失因数法中引入最小负荷率并提出了四个边界理论,使计算过程包含了更多的电流曲线信息,弱化了假定条件,其形状系数和负荷率通过统计获得,与实际情况存在偏差。
除计算精度以外,小电源的接入也是配电网线损计算中的典型问题。随着我国智能电网的建设以及电力市场的逐步推行,传统的集中式大电网供电模式已经无法满足当今社会对电力的需求,常见的小规模风力发电和光伏发电等分布式小电源在配电网中的占比越来越重,小电源接入配电网之后对线路电压及损耗的影响不容忽视。现有的等效容量法将小电源等值为具有负容量的专用配电变压器参与计算,将代表日容量分段进行等效,虽提高了计算精度,但依然没有充分利用小电源点的已知参数。以潮流计算为基础,考虑了小电源向配电网出力或者反供电力以及线路电压对损耗的影响,去除了节点电压相等的假设条件,由于配电网规模庞大,运行记录短缺,基于潮流法计算难以实现。
综上所述,现有技术中对于配电网中含分布式小电源的线损计算问题,尚缺乏有效的解决方案。
发明内容
为了克服上述现有技术的不足,本发明提供了一种计及小电源的配电网理论线损计算方法,该方法考虑了小电源接入点的电压功率以及实时发电量,以平均电流法为基础,加以修正,使线损计算更加精确。
本发明所采用的技术方案是:
一种计及小电源的配电网理论线损计算方法,该方法包括以下步骤:
获取配电网各支路的阻抗,各节点的有功功率、无功功率、电压和电流数据;
统计各小电源节点的用电量,采用平均电流法计算配电网的电能损耗;
根据各节点的有功功率、无功功率、电压幅值和电流幅值,计算各节点的注入电流;
采用搜索法形成路径互阻矩阵;
基于路径互阻矩阵和各节点的注入电流,采用回路分析法得到精确的电能损耗计算模型和基于平均电流法的损耗计算模型;
利用精确的电能损耗计算模型和基于平均电流法的损耗计算模型,得到平均电流法的计算误差,即损耗的修正量;
计算平均电流法得到的电能损耗与损耗的修正量之和,得到配电网的理论电能损耗。
作为本发明的进一步限定,所述采用平均电流法计算配电网的电能损耗的方法为:
计算各小电源节点的平均电流;
由末端小电源节点开始,逐段向首端节点方向迭代相加,得到相应的每一支路的平均电流;
利用等效系数将每一支路的平均电流转化为均方根电流;
将每一支路的均方根电流与该支路的电阻、运行时间相乘,得到该支路的电能损耗;
将所有支路的电能损耗进行累加,得到整个配电网总电能损耗。
作为本发明的进一步限定,所述小电源节点的平均电流的计算方法为:
计算日首端小电源节点的平均电流Iav0,并统计小电源节点j的日用电量Aa(j)以及所有小电源节点的总用电量;
将小电源节点j的日用电量Aa(j)与所有小电源节点的总用电量的比值,与日首端小电源节点的平均电流Iav0相乘,得到该小电源节点j的平均电流。
作为本发明的进一步限定,所述路径互阻矩阵的具体形成方法为:
将馈线首端节点作为根节点,采用深度优先原则对全网节点进行编号,并对支路进行编号;
设路径互阻矩阵内各矩阵元素初值为0,依次遍历配电网内所有支路,得到各节点的路径所经支路电阻之和,以及各节点的路径与小电源节点到每个节点的路径重叠支路电阻之和,形成路径互阻矩阵。
作为本发明的进一步限定,所述对支路进行编号原则为:支路靠近馈线首端的一端为首节点,远离馈线首端的一端为末节点,支路末节点号即为支路编号。
作为本发明的进一步限定,采用回路分析法得到精确的电能损耗计算模型的步骤包括:
利用各支路的阻抗,构建支路阻抗矩阵Zb
根据各节点的注入电流,构建回路电流矩阵is
令B为网络路径矩阵,该网络矩阵元素Bij为:
Figure GDA0002498528380000041
其中,
Figure GDA0002498528380000042
表示支路li包含于路径
Figure GDA0002498528380000043
i=1,2…b,j=1,2…N,b=N,
Figure GDA0002498528380000044
表示支路li不包含于路径
Figure GDA0002498528380000045
将网络路径矩阵的共轭转置BT与回路电流矩阵is相乘,得到支路电流矩阵ib
计算支路电流矩阵ib的共轭转置与支路阻抗矩阵Zb、支路电流矩阵ib的乘积,得到全网复功率损耗模型;
该全网复功率损耗模型的实部即为系统精确的电能损耗计算模型。
作为本发明的进一步限定,所述精确的电能损耗计算模型为:
Figure GDA0002498528380000051
其中,Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为路径互阻矩阵对角线元素和非对角线元素的实部,θij为第i,j个节点的电流相角差。
作为本发明的进一步限定,采用回路分析法得到基于平均电流法的损耗计算模型的步骤包括:
利用各支路的阻抗,构建支路阻抗矩阵Zb
依据平均电流法的假设条件,各小电源节点功率因数和负荷曲线均与馈线首端相同,则各电流相角相等;
根据各支路电流幅值,构建电流支路幅值矩阵Is,将电流支路幅值矩阵Is与电流相角的余弦值相乘,得到回路电流矩阵is
令B为网络路径矩阵,该网络矩阵元素Bij为:
Figure GDA0002498528380000052
其中,
Figure GDA0002498528380000053
表示支路li包含于路径
Figure GDA0002498528380000054
i=1,2…b,j=1,2…N,b=N,
Figure GDA0002498528380000055
表示支路li不包含于路径
Figure GDA0002498528380000056
将网络路径矩阵的共轭转置BT与回路电流矩阵is相乘,得到支路电流矩阵ib
计算支路电流矩阵ib的共轭转置与支路阻抗矩阵Zb、支路电流矩阵ib的乘积,得到全网复功率损耗模型;
该全网复功率损耗模型的实部即为系统精确的基于平均电流法的损耗计算模型。
作为本发明的进一步限定,所述基于平均电流法的损耗计算模型为:
Figure GDA0002498528380000061
其中,∑R为路径互阻矩阵的实部构成的矩阵,
Figure GDA0002498528380000062
I1,I2,…IN为支路电流幅值。
作为本发明的进一步限定,所述损耗修正量为:
Figure GDA0002498528380000063
其中,Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为路径互阻矩阵对角线元素和非对角线元素的实部,θij为第i,j个节点的电流相角差。
与现有技术相比,本发明的有益效果是:
(1)本发明提出的方法提高了含小电源的配电网理论线损计算的精度,计算结果更逼近实际值。其中小电源供电量与全网总发电量占比为41%左右时,本文算法几乎与实际值相等,修正效果最明显;
(2)本发明的该方法是在平均电流法的基础上增加一个修正量,可直接搜索得到修正部分的路径互阻矩阵Rij,在平均电流法相同的假设条件下,同时考虑了小电源接入点的功率、电压及其发电实时曲线,计算过程简单明了,实现了含分布式小电源的线损简化计算,对含分布式小电源的损耗计算具有很重要的现实意义。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1是本发明实施例一公开的计及小电源的配电网理论线损计算方法流程图;
图2是本发明实施例二公开的计及小电源的配电网理论线损计算方法流程图;
图3是一个普通配电网络图;
图4是本发明实施三的照政线路系统图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
正如背景技术所介绍的,现有的等效容量法将小电源等值为具有负容量的专用配电变压器参与计算,将代表日容量分段进行等效,虽提高了计算精度,但依然没有充分利用小电源点的已知参数。以潮流计算为基础,考虑了小电源向配电网出力或者反供电力以及线路电压对损耗的影响,去除了节点电压相等的假设条件,由于配电网规模庞大,运行记录短缺,基于潮流法计算难以实现。
为了解决如上的技术问题,本发明实施例一提出了一种计及小电源的配电网理论线损计算方法。如图1所示,该方法包括以下步骤:
S101,获取配电网各支路的阻抗,各节点的有功功率、无功功率、电压和电流数据。
S102,统计各小电源的发电量,采用平均电流法进行配电网损耗计算,得到损耗
Figure GDA0002498528380000081
S103,忽略电压损耗,根据各小电源节点的有功功率、无功功率、电压幅值和电流幅值数据,计算各小电源节点的注入电流Ik∠θk
为了计算小电源点的相角,因此对各小电源节点的注入电流进行计算,为后续修正量计算提供依据。
S104,由搜索法直接形成路径互阻矩阵。
所述路径互阻矩阵具体形成过程如下:
S1041,将馈线首端节点作为根节点,节点号为0,采用深度优先原则对全网节点进行编号。
S1042,对支路进行编号。
对支路进行编号原则为:定义支路距离馈线首端近的一端为首节点,远的一端为末节点,支路末节点号即为支路编号。
S1043,设路径互阻矩阵的矩阵元素初值:Rii=0,i=0,...,N,Rij=0,i=1,...,N-1,j=1,...,N。
S1044,依次遍历配电网内所有支路,得到路径互阻矩阵的元素Rii、Rij的值,形成路径互阻矩阵。
对于支路i,首节点为k,则Rii=Rkk+ri;若
Figure GDA0002498528380000082
即支路i包含于路径
Figure GDA0002498528380000083
中,j=i+1,...,N,则Rij=Rij+rj,其中ri为支路i的电阻。
S105,基于回路分析法推导得到精确的损耗Ploss计算模型以及基于平均电流法的损耗
Figure GDA0002498528380000091
计算模型;由精确的损耗Ploss计算模型以及基于平均电流法的损耗
Figure GDA0002498528380000092
计算模型得到的损耗的误差项,即修正量ΔP。
基于步骤S103得到的路径互阻矩阵,利用回路分析法推导得到精确的损耗Ploss计算模型为:
Figure GDA0002498528380000093
其中,is为回路电流,
Figure GDA0002498528380000094
为is的共轭转置,Ii、Ij分别为第i,j个节点的电流辐值,cosθij为第i,j个节点电流相角差的余弦值;Rii和Rij分别为路径互阻矩阵的矩阵元素,Rii表示支路li的路径
Figure GDA00024985283800000911
所经支路电阻之和;Rij表示支路li的路径
Figure GDA0002498528380000095
与路径
Figure GDA0002498528380000096
重叠支路电阻之和。
由精确的损耗Ploss计算模型以及基于平均电流法的损耗
Figure GDA0002498528380000097
计算模型的误差即为修正量ΔP,为:
Figure GDA0002498528380000098
其中,Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为路径互阻矩阵对角线元素和非对角线元素的实部,θij为第i,j个节点的电流相角差。
S106,最终损耗为
Figure GDA0002498528380000099
该最终损耗是由所述步骤S101中平均电流法计算得到的损耗
Figure GDA00024985283800000910
与所述步骤(4)中的修正量ΔP之和。
图2是本发明实施例二的计及小电源的配电网理论线损计算方法。如图2所示,该方法包括以下步骤:
S201,统计小电源的发电量,将小电源等值为负荷,采用平均电流法进行损耗计算,得到整个配电网的总电能损耗。
利用均方根电流法计算电能损耗公式如下:
Figure GDA0002498528380000101
式中,Ijf为代表日某时刻某支路的均方根电流,r为该支路的电阻,Δt为运行时间。
均方根电流可通过代表日实测电流或者通过该支路的有功功率、无功功率和电压计算得到,工作量比较大且不易获得,因而,实际工程常采用简化的平均电流法,利用等效系数将平均电流转化成均方根电流:
Figure GDA0002498528380000102
其中,Ijf0,Iav0分别为代表日首端节点的均方根电流和平均电流,K为等效系数。
已知选取代表日的配电线路各元件电阻、馈线首端代表日的负荷曲线及有功、无功电量、首端代表日的电压曲线,以及各负荷的用电量。作如下假设:各负荷节点负荷曲线的形状与首端相同;各负荷节点功率因数与首端相等;忽略沿线的电压损失对能耗的影响。则各负荷节点电流为:
Figure GDA0002498528380000103
其中,Iav(j)为第j节点的平均电流,j为负荷节点序号,n为负荷节点总数,Aa(j)为第j节点的日用电量。
由末端负荷节点开始,逐段向根节点方向迭代相加,得到相应的每一段线路的平均电流,则每一支路的电能损耗:
Figure GDA0002498528380000111
其中,r为相应支路的电阻,K为等效系数。
累加所有支路的电能损耗,得到整个配电网总电能损耗:
Figure GDA0002498528380000112
其中,ΔAL为整个配电网总电能损耗。
S202,忽略电压损耗,根据各负荷节点的有功功率P,无功功率Q,电压幅值U,电流幅值I数据,假设各负荷节点电压相角与首端相同为∠0°,根据公式
Figure GDA0002498528380000113
计算各负荷注入节点的电流Ik∠θk
S203,由搜索法直接形成路径互阻矩阵。
在小电源运行不满足假设条件下,直接将小电源等值为负荷,会带来一定的计算误差,因而有必要根据其运行参数进行修正。
假设N+1节点的配电网共有小电源m个,其负荷节点集合为Gm={g1,g2,...,gm},除小电源外的其余N-m负荷节点,各负荷曲线与首端负荷曲线相同,且功率因数也相同,则路径互阻矩阵R仅与各节点与小电源节点之间的路径互阻相关。路径互阻矩阵R中向量Rij定义如下:
当配电网中接入m个小电源时,Rij为一个m×N维的行向量,第j列元素为路径
Figure GDA0002498528380000121
遍历支路与小电源所在路径pi遍历支路重叠的支路电阻之和。
所述路径互阻矩阵的形成具体步骤如下:
S2031,将馈线首端节点作为根节点,节点号为0,采用深度优先原则实行全网节点编号;
S2032,定义支路距离馈线近的一端为首节点,远的一端为末节点,支路末节点号即为支路编号;
S2033,设路径互阻矩阵内各元素初值为0,即:
Rii=0,i=0,...,N,Rij=0,i=1,...,N-1,j=1,...,N;
S2034,依次遍历所有支路,对于支路li,首节点为k,则Rii=Rkk+ri;若
Figure GDA0002498528380000122
即支路i包含于路径
Figure GDA0002498528380000123
中,j=i+1,...,N,则Rij=Rij+rj。其中ri为支路li的电阻。
S204,计算修正量ΔP。
假定配电网共有N+1个节点,每个节点均可接入外部电路(电源、负荷或二者混合),则该网络图可用一有向图M表示,其中,图中顶点为电网中的节点,电网中支路即为有向图中的边,边的方向即为电流流动的方向。网络中各节点对应的顶点标记为1,......,N+1,以
Figure GDA0002498528380000124
表示公共接地点,该点电位
Figure GDA0002498528380000125
任意连接外部电路的节点k的电压为uk。如图4所示,该有向图可表示为M=c∪ε,其中c由配电网支路及节点构成,称之为内部子图,ε作为M的外部子图,连接配电网每个节点到公共地点,即顶点
Figure GDA0002498528380000126
对于没有连支仅由树支构成的图,共有N+1个节点,b条支路,则b=N。
令Zb=Rb+jXb为支路阻抗矩阵:
Figure GDA0002498528380000131
其中,Zb为对角阵,其中Zi为支路li的阻抗,i=1,2…b。该矩阵实部Rb为支路电阻矩阵:
Figure GDA0002498528380000132
其中,r1,......,rb为支路电阻。
图4中,M中所有支路构成配电网最大支撑树,该支撑树上,小电源节点j到每个节点的路径记为
Figure GDA0002498528380000133
路径
Figure GDA0002498528380000134
与ε的边构成一基本回路,全网共有回路数N,其回路电流为:
Figure GDA0002498528380000135
其中,i1…iN为各节点的注入电流,I1,I2,…IN为支路电流幅值,θ1…θN为电流相角。
令B为网络路径矩阵,该矩阵为N×N维的矩阵,其矩阵元素Bij定义如下:
Figure GDA0002498528380000136
其中,
Figure GDA0002498528380000137
表示支路li包含于路径
Figure GDA0002498528380000138
(i=1,2…b,j=1,2…N),
Figure GDA0002498528380000139
支路li不包含于路径
Figure GDA0002498528380000141
令支路电流:
Figure GDA0002498528380000142
其中,il1,il2…ilb为各支路电流。
可知:
ib=BTis (8)
则,全网复功率损耗:
Figure GDA0002498528380000143
其中,
Figure GDA0002498528380000144
分别为ib,is的共轭转置,∑=BZbBT为回路阻抗矩阵。
依据定义,可知∑为对称阵,其中对角元Σii为节点i的路径
Figure GDA0002498528380000149
所经支路阻抗之和,即路径自阻抗;非对角元Σij为路径
Figure GDA0002498528380000147
与路径
Figure GDA0002498528380000148
重叠支路阻抗之和,即路径互阻抗;其表达式为:
Figure GDA0002498528380000145
其中,Zi为支路li的阻抗;i,j=1,2,...,N。
式(9)的实部即为系统的有功损耗:
Figure GDA0002498528380000146
其中Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为∑矩阵对角线元素和非对角线元素的实部。ii=Ii∠θi,其中Ii为节点i电流幅值,θi为节点i电流相角。
根据平均电流法的假设条件,各负荷节点功率因数和负荷曲线均与馈线首端相同,并且忽略延线电压损耗,则有θ1=θ2....=θ0。式(6)可写为:
Figure GDA0002498528380000151
代入式(9),得到网络功率损耗:
Figure GDA0002498528380000152
式(13)实部即为有功功率损耗:
Figure GDA0002498528380000153
其中,∑R为∑的实部构成的矩阵。
比较式(11)、(14)可得平均电流法的计算误差为:
Figure GDA0002498528380000154
其中Rij为ΣR中相应元素。
式(15)即为修正量,其表明,当平均电流法的前提假设条件不能满足时,cosθij≠1,注入电流较大或者路径互阻较大时,该误差项不能忽略。
S205,最终损耗为
Figure GDA0002498528380000155
以山东省某地区照政线路为例,本发明实施例三是一种计及小电源的配电网理论线损计算方法的具体实施例,该线路结构如图4所示,共有22个节点,其中节点10和节点16接入负荷较重,其余均为一般负荷和轻负荷,节点14接入十字埠发电机。为进行比较,依据馈线首端、十字埠发电机的整点运行参数及各负荷的电能数据,利用匹配潮流法进行计算,将计算结果作为参考真值。
已知小电源的电压,功率及其发电量,分别采用平均电流法与改进算法进行线损计算,计算结果与真值进行对比。
为方便分析,将整个配电线路划分为3个片区,如图4所示。针对十字埠发电机发电量多少,可分为以下几种情形进行讨论:
场景一:小电源的供电量仅供片区1负荷,在此情形下,小电源供电量满足片区1负荷的20%,50%,80%,100%几种情况,计算结果如下表:
表1场景1计算结果
Figure GDA0002498528380000161
场景二:小电源的供电量满足片区1全部负荷需求外,同时满足片区2部分负荷需求20%,50%,80%,100%几种情况,计算结果如下表:
表2场景2计算结果
Figure GDA0002498528380000171
场景三:小电压的供电量满足1,2全部负荷需求外,同时满足片区3部分负荷需求20%,50%,80%,100%几种情况,计算结果如下表:
表3场景3计算结果
Figure GDA0002498528380000172
由上述表格数据可以看出,本发明充分利用了接入电源点的已知参数,提高了含小电源的配电网理论线损计算的精度,计算结果比传统的平均电流法更逼近实际值,验证了本发明方法的有效性。
从以上的描述中,可以看出,本申请上述的实施例实现了如下技术效果:
(1)本发明提出的方法提高了含小电源的配电网理论线损计算的精度,计算结果更逼近实际值。其中小电源供电量与全网总发电量占比为41%左右时,本文算法几乎与实际值相等,修正效果最明显;
(2)本发明的该方法是在平均电流法的基础上增加一个修正量,可直接搜索得到修正部分的路径互阻矩阵Rij,在平均电流法相同的假设条件下,同时考虑了小电源接入点的功率、电压及其发电实时曲线,计算过程简单明了,实现了含分布式小电源的线损简化计算,对含分布式小电源的损耗计算具有很重要的现实意义。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (10)

1.一种计及小电源的配电网理论线损计算方法,其特征是,包括以下步骤:
获取配电网各支路的阻抗,各节点的有功功率、无功功率、电压和电流数据;
统计各小电源节点的用电量,采用平均电流法计算配电网的电能损耗;
根据各节点的有功功率、无功功率、电压幅值和电流幅值,计算各节点的注入电流;
采用搜索法形成路径互阻矩阵;
基于路径互阻矩阵和各节点的注入电流,采用回路分析法得到精确的电能损耗计算模型和基于平均电流法的损耗计算模型;
利用精确的电能损耗计算模型和基于平均电流法的损耗计算模型,得到平均电流法的计算误差,即损耗的修正量;
计算平均电流法得到的电能损耗与损耗的修正量之和,得到配电网的理论电能损耗。
2.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,所述采用平均电流法计算配电网的电能损耗的方法为:
计算各小电源节点的平均电流;
由末端小电源节点开始,逐段向首端节点方向迭代相加,得到相应的每一支路的平均电流;
利用等效系数将每一支路的平均电流转化为均方根电流;
将每一支路的均方根电流与该支路的电阻、运行时间相乘,得到该支路的电能损耗;
将所有支路的电能损耗进行累加,得到整个配电网总电能损耗。
3.根据权利要求2所述的计及小电源的配电网理论线损计算方法,其特征是,所述小电源节点的平均电流的计算方法为:
计算日首端小电源节点的平均电流Iav0,并统计小电源节点j的日用电量Aa(j)以及所有小电源节点的总用电量;
将小电源节点j的日用电量Aa(j)与所有小电源节点的总用电量的比值,与日首端小电源节点的平均电流Iav0相乘,得到该小电源节点j的平均电流。
4.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,所述路径互阻矩阵的具体形成方法为:
将馈线首端节点作为根节点,采用深度优先原则对全网节点进行编号,并对支路进行编号;
设路径互阻矩阵内各矩阵元素初值为0,依次遍历配电网内所有支路,得到各节点的路径所经支路电阻之和,以及各节点的路径与小电源节点到每个节点的路径重叠支路电阻之和,形成路径互阻矩阵。
5.根据权利要求4所述的计及小电源的配电网理论线损计算方法,其特征是,所述对支路进行编号原则为:支路靠近馈线首端的一端为首节点,远离馈线首端的一端为末节点,支路末节点号即为支路编号。
6.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,采用回路分析法得到精确的电能损耗计算模型的步骤包括:
利用各支路的阻抗,构建支路阻抗矩阵Zb
根据各节点的注入电流,构建回路电流矩阵is
令B为网络路径矩阵,该网络矩阵元素Bij为:
Figure FDA0002498528370000031
其中,
Figure FDA0002498528370000032
表示支路li包含于路径
Figure FDA0002498528370000033
b=N,
Figure FDA0002498528370000034
表示支路li不包含于路径
Figure FDA0002498528370000035
将网络路径矩阵的共轭转置BT与回路电流矩阵is相乘,得到支路电流矩阵ib
计算支路电流矩阵ib的共轭转置与支路阻抗矩阵Zb、支路电流矩阵ib的乘积,得到全网复功率损耗模型;
该全网复功率损耗模型的实部即为系统精确的电能损耗计算模型。
7.根据权利要求6所述的计及小电源的配电网理论线损计算方法,其特征是,所述精确的电能损耗计算模型为:
Figure FDA0002498528370000036
其中,
Figure FDA0002498528370000037
为is的共轭转置,∑=BZbBT为回路阻抗矩阵,Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为路径互阻矩阵对角线元素和非对角线元素的实部,θij为第i,j个节点的电流相角差。
8.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,采用回路分析法得到基于平均电流法的损耗计算模型的步骤包括:
利用各支路的电阻和阻抗,构建支路阻抗矩阵Zb
依据平均电流法的假设条件,各小电源节点功率因数和负荷曲线均与馈线首端相同,则各电流相角相等;
根据各支路电流幅值,构建电流支路幅值矩阵Is,将电流支路幅值矩阵Is与电流相角的余弦值相乘,得到回路电流矩阵is
令B为网络路径矩阵,该网络矩阵元素Bij为:
Figure FDA0002498528370000041
其中,
Figure FDA0002498528370000042
表示支路li包含于路径
Figure FDA0002498528370000043
b=N,
Figure FDA0002498528370000044
表示支路li不包含于路径
Figure FDA0002498528370000045
将网络路径矩阵的共轭转置BT与回路电流矩阵is相乘,得到支路电流矩阵ib
计算支路电流矩阵ib的共轭转置与支路阻抗矩阵Zb、支路电流矩阵ib的乘积,得到全网复功率损耗模型;
该全网复功率损耗模型的实部即为系统精确的基于平均电流法的损耗计算模型。
9.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,所述基于平均电流法的损耗计算模型为:
Figure FDA0002498528370000046
其中,∑R为路径互阻矩阵的实部构成的矩阵,
Figure FDA0002498528370000047
I1,I2,…IN为支路电流幅值。
10.根据权利要求1所述的计及小电源的配电网理论线损计算方法,其特征是,所述损耗修正量为:
Figure FDA0002498528370000048
其中,Ii、Ij分别为第i,j个节点的电流辐值,Rii、Rij分别为路径互阻矩阵对角线元素和非对角线元素的实部,θij为第i,j个节点的电流相角差。
CN201811043738.0A 2018-09-07 2018-09-07 计及小电源的配电网理论线损计算方法 Expired - Fee Related CN109188204B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811043738.0A CN109188204B (zh) 2018-09-07 2018-09-07 计及小电源的配电网理论线损计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811043738.0A CN109188204B (zh) 2018-09-07 2018-09-07 计及小电源的配电网理论线损计算方法

Publications (2)

Publication Number Publication Date
CN109188204A CN109188204A (zh) 2019-01-11
CN109188204B true CN109188204B (zh) 2020-07-03

Family

ID=64915174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811043738.0A Expired - Fee Related CN109188204B (zh) 2018-09-07 2018-09-07 计及小电源的配电网理论线损计算方法

Country Status (1)

Country Link
CN (1) CN109188204B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110112733B (zh) * 2019-05-21 2023-03-14 华北水利水电大学 一种低压配电系统引入不平衡度的线损计算方法
CN110212521B (zh) * 2019-06-05 2022-04-22 国网河南省电力公司电力科学研究院 配电网线损计算方法及系统
CN110635471B (zh) * 2019-08-02 2023-06-16 国网四川省电力公司电力科学研究院 基于同时段功率和电量的台区电网理论线损计算方法
CN111786410B (zh) * 2020-07-13 2021-11-12 济南大学 一种基于网损误差的配电网网损分摊方法及系统
CN115201743B (zh) * 2022-09-16 2023-01-24 北京志翔科技股份有限公司 低压台区计量点误差的确定方法、装置及存储介质
CN115276004B (zh) * 2022-09-21 2023-02-03 石家庄科林电气股份有限公司 电网区块线损计算方法、装置及终端设备

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8324757B2 (en) * 2009-02-02 2012-12-04 Sony Corporation Power distributing system, power transmitting device, power receiving device, power transmitting method and power receiving method
CN103577679A (zh) * 2012-08-10 2014-02-12 深圳市龙电电气有限公司 一种低压台区理论线损实时计算方法
CN103606014A (zh) * 2013-10-18 2014-02-26 国家电网公司 一种基于多目标的海岛分布式电源优化方法
CN103904644A (zh) * 2014-03-26 2014-07-02 国家电网公司 一种基于分布式电源接入的智能变电站负荷自动分配方法
CN104200314A (zh) * 2014-08-22 2014-12-10 国家电网公司 基于数据远程同步采集技术的电量及线损同期统计方法
CN104200302A (zh) * 2014-08-04 2014-12-10 国家电网公司 用于电力企业同期线损管理的技术架构和方法
CN104821590A (zh) * 2015-04-25 2015-08-05 江阴方程电气自动化研究所有限公司 智能电网节能优化综合测评方法
CN105375477A (zh) * 2015-11-27 2016-03-02 国网重庆市电力公司电力科学研究院 一种农村配电网线损计算方法及系统
CN105427186A (zh) * 2015-11-17 2016-03-23 国网上海市电力公司 一种基于改进等值电阻法的配电网线损计算方法
CN106097156A (zh) * 2016-07-15 2016-11-09 广东电网有限责任公司电网规划研究中心 配电网线损计算方法和系统
CN107046286A (zh) * 2017-05-12 2017-08-15 国网上海市电力公司 基于前推回代法的低压配电台区理论线损计算方法
CN107846040A (zh) * 2017-11-02 2018-03-27 国网重庆市电力公司电力科学研究院 基于二阶锥松弛的分布式光伏与储能协调规划方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105117560B (zh) * 2015-09-15 2018-11-23 云南电网有限责任公司昆明供电局 一种中压配电线路理论线损计算方法
CN108254642B (zh) * 2016-01-31 2021-08-31 郑龙海 一种10kV配电网
CN107394797A (zh) * 2017-08-30 2017-11-24 华侨大学 一种考虑无功优化的分布式电源接入配电网的评估方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8324757B2 (en) * 2009-02-02 2012-12-04 Sony Corporation Power distributing system, power transmitting device, power receiving device, power transmitting method and power receiving method
CN103577679A (zh) * 2012-08-10 2014-02-12 深圳市龙电电气有限公司 一种低压台区理论线损实时计算方法
CN103606014A (zh) * 2013-10-18 2014-02-26 国家电网公司 一种基于多目标的海岛分布式电源优化方法
CN103904644A (zh) * 2014-03-26 2014-07-02 国家电网公司 一种基于分布式电源接入的智能变电站负荷自动分配方法
CN104200302A (zh) * 2014-08-04 2014-12-10 国家电网公司 用于电力企业同期线损管理的技术架构和方法
CN104200314A (zh) * 2014-08-22 2014-12-10 国家电网公司 基于数据远程同步采集技术的电量及线损同期统计方法
CN104821590A (zh) * 2015-04-25 2015-08-05 江阴方程电气自动化研究所有限公司 智能电网节能优化综合测评方法
CN105427186A (zh) * 2015-11-17 2016-03-23 国网上海市电力公司 一种基于改进等值电阻法的配电网线损计算方法
CN105375477A (zh) * 2015-11-27 2016-03-02 国网重庆市电力公司电力科学研究院 一种农村配电网线损计算方法及系统
CN106097156A (zh) * 2016-07-15 2016-11-09 广东电网有限责任公司电网规划研究中心 配电网线损计算方法和系统
CN107046286A (zh) * 2017-05-12 2017-08-15 国网上海市电力公司 基于前推回代法的低压配电台区理论线损计算方法
CN107846040A (zh) * 2017-11-02 2018-03-27 国网重庆市电力公司电力科学研究院 基于二阶锥松弛的分布式光伏与储能协调规划方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Loss allocation in radial distribution networks with various distributed;Kushal Manoharrao Jagtap;《electrical power and energy systems》;20150925;全文 *
线损精细化计算与评估方法的应用研究;林卓琼;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20161015;全文 *
配电网理论线损计算算法比较;刘静;《河北工业大学成人教育学院学报》;20081231;全文 *

Also Published As

Publication number Publication date
CN109188204A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109188204B (zh) 计及小电源的配电网理论线损计算方法
Tang et al. Real-time optimal power flow
Contreras et al. Improved assessment of the flexibility range of distribution grids using linear optimization
CN106777708B (zh) 一种电力-天然气区域综合能源系统的稳态分析方法
CN107666149B (zh) 一种中压配电网线损计算方法
US10164434B2 (en) Methods and systems for power injection or extraction in a power network
CN109474017A (zh) 一种配电网实时分布式经济调度方法
CN109149638A (zh) 基于mpc和admm算法的vsc-hvdc并网风电场分布式协调电压控制方法及系统
CN109802381A (zh) 一种基于模糊控制的直流微电网多源动态协调控制方法
CN109494735B (zh) 一种低压配电网线损计算模型及方法
CN106385035B (zh) 一种计及充电电容的直流潮流计算方法
CN100563078C (zh) 一种准确计算电力线路损耗的方法
CN113886767A (zh) 一种电力系统碳排放流实时计算方法
CN110635471A (zh) 基于同时段功率和电量的台区电网理论线损计算方法
CN112968441B (zh) 一种应用于大规模风电基地的电网规划方法
CN106410784B (zh) 变电站有功负荷对地区电网有功网损的灵敏度计算方法
CN113054688A (zh) 考虑不确定性的可再生能源集群出力分布式协调控制方法
CN107039981A (zh) 一种拟直流线性化概率最优潮流计算方法
CN109995025A (zh) 一种基于两部制电价体系的储能控制方法
CN110336322B (zh) 基于日最小负荷置信区间的光伏发电准入容量确定方法
CN107465195B (zh) 一种基于微电网结合潮流计算的最优潮流双层迭代方法
Jiang et al. An enhanced drop control method for DC microgrids with accurate current sharing and DC bus voltage restoration
CN103838962A (zh) 一种计及pmu量测的分步线性状态估计方法
CN103986158A (zh) 一种分布式电源配电网潮流计算方法
CN116111596A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200703

Termination date: 20210907