CN108092284B - 一种基于线性模型的三相不平衡智能配电网网络重构方法 - Google Patents

一种基于线性模型的三相不平衡智能配电网网络重构方法 Download PDF

Info

Publication number
CN108092284B
CN108092284B CN201711180647.7A CN201711180647A CN108092284B CN 108092284 B CN108092284 B CN 108092284B CN 201711180647 A CN201711180647 A CN 201711180647A CN 108092284 B CN108092284 B CN 108092284B
Authority
CN
China
Prior art keywords
power
phase
distribution network
branch
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.)
Active
Application number
CN201711180647.7A
Other languages
English (en)
Other versions
CN108092284A (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201711180647.7A priority Critical patent/CN108092284B/zh
Publication of CN108092284A publication Critical patent/CN108092284A/zh
Application granted granted Critical
Publication of CN108092284B publication Critical patent/CN108092284B/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
    • H02J3/26Arrangements for eliminating or reducing asymmetry in polyphase 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]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于线性模型的三相不平衡智能配电网网络重构方法,包括使用Distflow支路潮流模型对三相不平衡智能配电网进行建模与近似分析;建立考虑含分布式电源(Distributed Generation,DG)的智能配电网线路参数三相不对称和功率不平衡问题的网络重构数学模型;使用最佳等距分段线性逼近的方法将所建模型线性化,得到网络重构的混合整数线性模型并使用商业数学规划软件求解。本发明弥补了当前智能配电网网络重构优化过程中忽略线路参数不对称以及功率不平衡的不足,同时提高了网络重构模型的求解速度和效率,为智能配电网的经济和可靠运行提供了理论支撑。

Description

一种基于线性模型的三相不平衡智能配电网网络重构方法
技术领域
本发明涉及智能配电网运行的网络重构方法,特别是涉及一种考虑配电网三相功率不平衡和线路参数不对称的经济重构方法。
背景技术
智能配电网是一个由变电站母线、分段开关、联络开关、负荷和分布式电源等设备和元件通过一定形式连接而成的复杂网络。配网网络重构是通过合理改变系统中开关的开合状态,使配电网在满足安全运行和辐射状运行的条件下,达到既定优化目的的过程。分布式发电作为清洁能源发电技术的代表,以高效率、低能效得到日益广泛的关注。然而,接入配电网后,不仅改变潮流的流动方向,而且会影响系统的其他方面。因此,国内外的许多学者对分布式电源接入后的智能配电网重构模型进行了研究。
为实现不同的重构目标,网络重构的目标函数可以选线路网损最小、平衡线路负荷、开关操作次数最小,供电可靠性最优等。由于配网网络重构可包含一个或多个优化目标,同时结合其开关状态变量和潮流计算非线性特点,因此配网网络重构是一个含多目标的混合整数非线性优化问题,属于NP-hard难解问题。求解网络重构问题的主要方法包括启发式算法、人工智能算法和数学规划方法。虽然启发式方法的研究和开发工作实现了在实际配电网上的快速求解,但全局解的最优性仍有待确定。人工智能算法是一种随机的方法,主要由人工智能技术实现,非常适合于混合二进制优化问题的解决方案。然而,基于随机的解决方案通常是不可重复的,可能需要几个执行的算法,以获得最佳的解决方案。此外,与启发式方法相比,基于随机的方法并不要求减少计算工作量,对于单目标问题的全局最优解仍然不能被正式验证。
另一方面,确定性数学规划方法还没有被充分利用以解决网络重构问题,目前国内外的研究也比较热门,其特征在于将配电网网络重构题转化为混合整数线性和非线性模型来求解。一般认为,利用数学规划方法求解重构问题可以得到不依赖于配电网初始结构的全局最优解,并且基于确定性数学规划方法的解决方案是可重复的,通常能够保证全局最优解或者次全局最优解。但是,基于确定性数学规划方法的计算复杂度通常较高,计算效率也比较低,需要采取相应简化近似方法提高数学规划方法的实用性。
发明内容
发明目的:为弥补当前智能配电网网络重构优化过程中忽略线路参数不对称以及功率不平衡的不足,提供一种能够提高当前各相功率不平衡以及线路参数不对称的含分布式电源的智能配电网运行经济性的网络重构优化方法。
技术方案:一种基于线性模型的三相不平衡智能配电网网络重构方法,其特征在于,包括以下步骤:
(1)基于Distflow支路潮流模型对三相线路参数不对称、三相功率不平衡的含分布式电源的智能配电网进行总线路损耗的建模,并进行电压近似简化;
(2)根据步骤(1)得到的智能配电网进行总线路损耗模型,推导配电网三相总的线路损耗,并设总网损最优为配电网经济重构模型的目标函数,推导存在三相线路参数不对称和功率不平衡情况的配电网运行约束条件形式,建立起网络重构的数学模型;
(3)使用最佳等距分段线性逼近的方法将所建立的经济重构混合整数二次规划模型中的二阶函数近似为线性函数,获得近似误差可调的基于线性潮流的混合整数线性网络重构模型并使用商业数学线性规划软件求解。
进一步地,步骤(1)包括:
(a)对于辐射状运行的智能配电网,推导出其优化模型的Distflow支路潮流形式如下。
对于电网中节点j,有
Figure BDA0001479113790000021
Figure BDA0001479113790000022
式中,
Figure BDA0001479113790000023
为ABC三相模型;集合t(j)为电网中以j为末端节点的支路的首端节点集合;集合f(j)为电网中以j为首端节点的支路的末端节点集合;
Figure BDA0001479113790000024
为节点电压幅值;
Figure BDA0001479113790000025
分别为支路ij首端三相有功功率和无功功率;
Figure BDA0001479113790000026
为节点j的有功功率净注入值;
Figure BDA0001479113790000027
Figure BDA0001479113790000028
分别为节点j上所连接的DG有功功率和负荷有功功率;
Figure BDA0001479113790000029
为计及三相线路自阻抗和互阻抗的支路电阻;
(b)根据Distflow支路潮流形式,推导出智能配电网网络总损耗的表达形式如下:
Figure BDA0001479113790000031
实际配电网运行的节点电压标幺值一般在1附近,因此将节点电压近似为1,总损耗变为如下形式:
Figure BDA0001479113790000032
分母的二次项忽略后,总网损变为一个二阶函数的形式,便于进一步化简计算。
进一步地,步骤(2)包括:
(a1)根据步骤(1)得到的配电网总损耗表达式,建立配电网网络重构数学模型,目标函数设为总网损最优,其中控制变量设为网络拓扑变量和线路潮流变量,决策变量是网络拓扑变量,网络重构模型的目标函数描述如下:
Figure BDA0001479113790000033
由于智能配电网存在三相线路参数不对称以及三相线路潮流不平衡的情况,因此网络重构模型的目标函数需要考虑支路电阻的不对称和支路有功无功功率的不一致情况;其中,线路电阻
Figure BDA0001479113790000034
为一个3×3的矩阵,具体描述如下:
Figure BDA0001479113790000035
式中
Figure BDA0001479113790000036
除主对角线外的其他元素均为三相之间的互阻抗,
Figure BDA0001479113790000037
Figure BDA0001479113790000038
支路有功
Figure BDA0001479113790000039
为一个3×1的矩阵,具体描述如下:
Figure BDA00014791137900000310
根据上述的各参数矩阵形式,可以写出智能配电网网络重构模型各相总损耗之和的目标函数;
(b1)该模型的约束条件包括节点注入功率约束,网络辐射状约束以及支路潮流上下限约束。
进一步地,步骤(b1)包括的节点注入功率约束,网络辐射状约束以及支路潮流上下限约束具体如下:
配网网络重构模型的节点注入功率约束描述为:
Figure BDA0001479113790000041
Figure BDA0001479113790000042
式中,αij表示支路(i,j)的开合状态,αij=1表示该支路连通,αij=0表示该支路断开;N(i)表示网络中能够与节点i直接连接的所有节点集合;
由于智能配电网存在三相功率不平衡的情况,每个节点各相的节点注入功率都必须满足上面的功率平衡约束;其中DG的有功功率方向为正,负荷有功功率的方向为负;
配网网络重构模型的网络辐射状约束条件描述如下:
βijji=αij i=1,2,…,n j∈N(i)
Figure BDA0001479113790000043
β1j=0 j∈N(1)
βij∈{0,1} i=1,…,n j∈N(i)
0≤αij≤1 (i,j)∈E
式中,βij为生成树的父子节点关系变量,当节点j是节点i的父节点时βij=1,反之取0;其中第1式表明当支路(i,j)连接时,必有一个节点是另一个节点的父节点;第2式则表明节点i只存在一个父节点;第3式表明变电站1节点不是任何节点的父节点;每条支路的分段开关和联络开关为三相的,网络拓扑变量不用考虑各相不一致情况,因此当满足上述约束条件时,模型的拓扑优化结果即为辐射状网络;
配网网络重构模型的支路潮流上下限约束描述如下:
Figure BDA0001479113790000044
Figure BDA0001479113790000045
由于智能配电网存在三相功率不平衡的情况,因此约束条件考虑了各相的支路功率不相等的情况以及不同线路不同相的功率上下限也会存在差异的情况。
进一步地,步骤(3)包括:
(a2)将智能配电网网络重构模型目标函数中的支路功率二次项
Figure BDA0001479113790000051
表示为y=h(x)形式的二次函数,使用最佳等距分段线性逼近法对该功率二次项进行一阶线性化近似表示;首先说明最佳等距分段中对于最小三相支路功率的分段距离查找原理:假设y=h(x)为需要线性化处理的三相不平衡智能配电网有功和无功非线性支路功率函数曲线,曲线的上方和下方分别存在一条三相误差功率上限曲线和三相误差功率下限曲线,分别可用
Figure BDA0001479113790000052
yh(x)=h(x)-δ表示,其中,δ为给定的最大三相功率逼近误差常数;x1点为三相误差功率下限曲线h(x)中斜率较大的一点,以该端点(x1,h(x1))作为起点做三相误差功率上限曲线的切线,可得到切点为
Figure BDA0001479113790000058
并与三相误差功率下限曲线相交于点(x2,h(x2)),该切线斜率为
Figure BDA0001479113790000053
接着以交点(x2,h(x2))为起点做三相误差功率上限曲线的切线,可得到切点为
Figure BDA0001479113790000054
并与三相误差功率下限曲线相交于点(x3,h(x3)),并得该切线斜率
Figure BDA0001479113790000055
反复执行上述操作直至
Figure BDA0001479113790000056
的终点,并得到一系列三相不平衡配电网的三相支路功率函数的分段直线;
对上述三相不平衡配电网支路功率曲线切线的斜率进行化简可得:
Figure BDA0001479113790000057
当x1和δ为已知量时,可以类推分别求得配电网重构模型中所有三相不平衡功率函数分段直线的端点值x2,x3,…,xl;基于上述三相功率端点值,可获得相邻功率端点的最小功率间距值C,以该间距值作为等距三相支路功率的分段标准,对三相不平衡配电网支路功率函数y=h(x)进行等距分段可得到新的三相不平衡支路功率分段直线端点,分别为
Figure BDA0001479113790000061
以此可获得等距分段线性逼近三相不平衡支路功率函数以表示原模型中的三相不平衡配电网功率二次函数;
经过上述对三相不平衡配电网模型中非线性功率项的最佳等距分段线性逼近,可获得配电网三相支路功率
Figure BDA0001479113790000062
的分段线性化近似表示式:
Figure BDA0001479113790000063
式中,Ωt,Ωs分别为
Figure BDA0001479113790000064
分段后的三相功率断面数集合;
Figure BDA0001479113790000065
分别为三相不平衡支路功率函数
Figure BDA0001479113790000066
各个断面上三相线性功率函数的斜率;
Figure BDA0001479113790000067
分别为
Figure BDA0001479113790000068
各断面上三相线性功率函数的取值,且
Figure BDA0001479113790000069
通过上述对三相不平衡配电网支路非线性功率项进行线性化近似,三相不平衡配电网重构模型的三相支路功率变量
Figure BDA00014791137900000610
可转化为由三相分段线性功率变量
Figure BDA00014791137900000611
构成的线性表达式:
Figure BDA00014791137900000612
(b2)通过线性化变换,三相不平衡配电网重构模型目标函数可转化为三相支路功率的线性函数,从而使得三相不平衡智能配电网网络重构问题转化为混合整数线性规划问题;由于三相线性支路功率变量
Figure BDA00014791137900000613
被限制为非负,而实际中当系统潮流与假定潮流正方向相反时,潮流值为负值,与设定不符,因此引入以下四个附加的三相不平衡潮流功率变量
Figure BDA00014791137900000614
并增加相应约束条件如下:
建立附加三相不平衡潮流功率变量与原三相支路功率变量
Figure BDA00014791137900000615
的关系约束如下:
Figure BDA0001479113790000071
附加三相变量自身约束如下:
Figure BDA0001479113790000072
Figure BDA0001479113790000073
Figure BDA0001479113790000074
各断面上三相不平衡线性功率函数的三相线性支路变量
Figure BDA0001479113790000075
Figure BDA0001479113790000076
的上下限约束:
Figure BDA0001479113790000077
(c2)使用商业数学线性规划求解器,即可对上述混合整数线性规划模型高效求解,通过设置分段后的断面数集合大小,可以调整网络重构线性模型的计算精度和效率。
有益效果:与现有技术相比,本发明的方法具有以下优点:
(1)本发明可以弥补当前含分布式电源的智能配电网网络重构优化过程忽略三相功率不平衡和线路参数不对称情况的不足,提高了配电网网络重构优化算法的工程应用能力,为智能配电网经济运行的优化研究提供理论支撑。
(2)本发明提出的智能配电网网络重构方法相比于现有的配网重构方法而言可以更高效地实现对于数学模型的求解,并且经济重构模型的优化方案具有可重复性,能够得到全局最优解或者次全局最优解,模型具有更优的求解性能。
(3)本发明提出的智能配电网网络重构方法的目标函数计算近似误差可以根据实际工程应用需要人为自由调节,在保证求解精度的同时,可以大幅节约求解资源,提高了实际工程应用的价值。
附图说明
图1为本发明方法的流程图;
图2为最佳等距分段法分段间距查找示意图;
图3为分段线性逼近法原理示意图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
如图1所示,一种基于线性模型的三相不平衡智能配电网网络重构方法,包括以下步骤:
步骤1:对于辐射状智能配电网,基于Distflow支路潮流模型对三相线路参数不对称、三相功率不平衡的含分布式电源的智能配电网进行总线路损耗的建模,并进行电压近似简化。
(1)对于辐射状运行的智能配电网,推导出其优化模型的Distflow支路潮流形式如下。
对于电网中节点j,有
Figure BDA0001479113790000081
Figure BDA0001479113790000082
式中,
Figure BDA0001479113790000083
为ABC三相模型;集合t(j)为电网中以j为末端节点的支路的首端节点集合;集合f(j)为电网中以j为首端节点的支路的末端节点集合;
Figure BDA0001479113790000084
为节点电压幅值;
Figure BDA0001479113790000085
分别为支路ij首端三相有功功率和无功功率;
Figure BDA0001479113790000086
为节点j的有功功率净注入值;
Figure BDA0001479113790000087
Figure BDA0001479113790000088
分别为节点j上所连接的DG有功功率和负荷有功功率;
Figure BDA0001479113790000089
为计及三相线路自阻抗和互阻抗的支路电阻;
(2)根据Distflow支路潮流形式,推导出智能配电网网络总损耗的表达形式如下:
Figure BDA00014791137900000810
由于分母存在节点电压二次项,该式为典型的非凸函数,现有的求解算法容易陷入局部最优,难以获得全局最优解。而实际配电网运行的节点电压标幺值一般在1附近,因此将节点电压近似为1,总损耗变为如下形式:
Figure BDA00014791137900000811
分母的二次项忽略后,总网损变为一个二阶函数的形式,便于进一步化简计算。
步骤2:推导配电网三相总的线路损耗,并设总网损最优为配电网经济重构模型的目标函数,推导存在三相线路参数不对称和功率不平衡情况的配电网运行约束条件形式,建立起网络重构的数学模型;
基于步骤1提出的基于Distflow支路潮流模型的智能配电网总线路损耗的近似表达式,针对智能配电网的运行特点,分别列写出网络重构模型的目标函数和约束条件,主要实施步骤如下所述:
(1)根据步骤1得到的配电网总损耗表达式,建立配电网网络重构数学模型,目标函数设为总网损最优,其中控制变量设为网络拓扑变量和线路潮流变量,决策变量是网络拓扑变量,网络重构模型的目标函数描述如下:
Figure BDA0001479113790000091
由于智能配电网存在三相线路参数不对称以及三相线路潮流不平衡的情况,因此网络重构模型的目标函数需要考虑支路电阻的不对称和支路有功无功功率的不一致情况。其中,线路电阻
Figure BDA0001479113790000092
为一个3×3的矩阵,具体描述如下:
Figure BDA0001479113790000093
式中
Figure BDA0001479113790000094
除主对角线外的其他元素均为三相之间的互阻抗,
Figure BDA0001479113790000095
Figure BDA0001479113790000096
支路有功
Figure BDA0001479113790000097
为一个3×1的矩阵,具体描述如下:
Figure BDA0001479113790000098
根据上述的各参数矩阵形式,可以写出智能配电网网络重构模型各相总损耗之和的目标函数;
(2)该模型的约束条件包括节点注入功率约束,网络辐射状约束以及支路潮流上下限约束。
1)配网网络重构模型的节点注入功率约束描述如下:
Figure BDA0001479113790000101
Figure BDA0001479113790000102
式中,αij表示支路(i,j)的开合状态,αij=1表示该支路连通,αij=0表示该支路断开;N(i)表示网络中能够与节点i直接连接的所有节点集合。
由于智能配电网存在三相功率不平衡的情况,每个节点各相的节点注入功率都必须满足上面的功率平衡约束。其中DG的有功功率方向为正,负荷有功功率的方向为负;
2)配网网络重构模型的网络辐射状约束条件描述如下:
βijji=αij i=1,2,…,n j∈N(i)
Figure BDA0001479113790000103
β1j=0 j∈N(1)
βij∈{0,1} i=1,…,n j∈N(i)
0≤αij≤1 (i,j)∈E
式中,βij为生成树的父子节点关系变量,当节点j是节点i的父节点时βij=1,反之取0。其中第1式表明当支路(i,j)连接时,必有一个节点是另一个节点的父节点;第2式则表明节点i只存在一个父节点;第3式表明变电站1节点不是任何节点的父节点。每条支路的分段开关和联络开关为三相的,网络拓扑变量不用考虑各相不一致情况,因此当满足上述约束条件时,模型的拓扑优化结果即为辐射状网络;
3)配网网络重构模型的支路潮流上下限约束描述如下:
Figure BDA0001479113790000104
Figure BDA0001479113790000105
由于智能配电网存在三相功率不平衡的情况,因此约束条件考虑了各相的支路功率不相等的情况以及不同线路不同相的功率上下限也会存在差异的情况。
步骤3:使用最佳等距分段线性逼近的方法将所建立的经济重构混合整数二次规划模型中的二阶函数近似为线性函数,获得近似误差可调的基于线性潮流的混合整数线性网络重构模型并使用商业数学线性规划软件求解。
根据分段线性逼近的近似方法,将步骤2所提出的智能配电网网络重构二次规划模型线性化,对原模型的目标函数进行改写近似,同时需要增加额外的附加变量和约束条件,具体如下所示:
(a)使用最佳等距分段线性逼近法对智能配电网网络重构模型目标函数中的二次项
Figure BDA0001479113790000111
进行一阶线性化近似表示,以形如y=h(x)的二次函数为例说明最佳等距分段线性逼近法。首先说明最佳等距分段中最小分段距离查找原理:如图2所示,假设y=h(x)为需要线性化处理的三相不平衡智能配电网有功和无功非线性支路功率函数曲线,曲线的上方和下方分别存在一条三相误差功率上限曲线和三相误差功率下限曲线,分别可用
Figure BDA0001479113790000112
Figure BDA0001479113790000113
表示,其中,δ为给定的最大三相功率逼近误差常数;x1点为三相误差功率下限曲线h(x)中斜率较大的一点,以该端点(x1,h(x1))作为起点做三相误差功率上限曲线的切线,可得到切点为
Figure BDA0001479113790000114
并与三相误差功率下限曲线相交于点(x2,h(x2)),该切线斜率为
Figure BDA0001479113790000115
接着以交点(x2,h(x2))为起点做三相误差功率上限曲线的切线,可得到切点为
Figure BDA0001479113790000116
并与三相误差功率下限曲线相交于点
Figure BDA0001479113790000117
并得该切线斜率
Figure BDA0001479113790000118
反复执行上述操作直至
Figure BDA0001479113790000119
的终点,并得到一系列三相不平衡配电网的三相支路功率函数的分段直线;
对上述三相不平衡配电网支路功率曲线切线的斜率进行化简可得:
Figure BDA00014791137900001110
当x1和δ为已知量时,可以类推分别求得配电网重构模型中所有三相不平衡功率函数分段直线的端点值x2,x3,…,xl;基于上述三相功率端点值,可获得相邻功率端点的最小功率间距值C,以该间距值作为等距三相支路功率的分段标准,对三相不平衡配电网支路功率函数y=h(x)进行等距分段可得到新的三相不平衡支路功率分段直线端点,分别为
Figure BDA0001479113790000121
以此可获得等距分段线性逼近三相不平衡支路功率函数以表示原模型中的三相不平衡配电网功率二次函数,如图3所示;
经过上述对三相不平衡配电网模型中非线性功率项的最佳等距分段线性逼近,可获得配电网三相支路功率
Figure BDA0001479113790000122
的分段线性化近似表示式:
Figure BDA0001479113790000123
式中,Ωt,Ωs分别为
Figure BDA0001479113790000124
分段后的三相功率断面数集合;
Figure BDA0001479113790000125
分别为三相不平衡支路功率函数
Figure BDA0001479113790000126
各个断面上三相线性功率函数的斜率;
Figure BDA0001479113790000127
分别为
Figure BDA0001479113790000128
各断面上三相线性功率函数的取值,且
Figure BDA0001479113790000129
通过上述对三相不平衡配电网支路非线性功率项进行线性化近似,三相不平衡配电网重构模型的三相支路功率变量
Figure BDA00014791137900001210
可转化为由三相分段线性功率变量
Figure BDA00014791137900001211
构成的线性表达式:
Figure BDA00014791137900001212
(b2)通过线性化变换,三相不平衡配电网重构模型目标函数可转化为三相支路功率的线性函数,从而使得三相不平衡智能配电网网络重构问题转化为混合整数线性规划问题;由于三相线性支路功率变量
Figure BDA00014791137900001213
被限制为非负,而实际中当系统潮流与假定潮流正方向相反时,潮流值为负值,与设定不符,因此引入以下四个附加的三相不平衡潮流功率变量
Figure BDA00014791137900001214
并增加相应约束条件如下:
建立附加三相不平衡潮流功率变量与原三相支路功率变量
Figure BDA00014791137900001215
的关系约束如下:
Figure BDA0001479113790000131
附加三相变量自身约束如下:
Figure BDA0001479113790000132
Figure BDA0001479113790000133
Figure BDA0001479113790000134
各断面上三相不平衡线性功率函数的三相线性支路变量
Figure BDA0001479113790000135
Figure BDA0001479113790000136
的上下限约束:
Figure BDA0001479113790000137
(c2)使用商业数学线性规划求解器,即可对上述混合整数线性规划模型高效求解,通过设置分段后的断面数集合大小,可以调整网络重构线性模型的计算精度和效率。

Claims (2)

1.一种基于线性模型的三相不平衡智能配电网网络重构方法,其特征在于,包括以下步骤:
(1)基于Distflow支路潮流模型对三相线路参数不对称、三相功率不平衡的含分布式电源的智能配电网进行总线路损耗的建模,并进行电压近似简化;
(2)根据步骤(1)得到的智能配电网进行总线路损耗模型,推导配电网三相总的线路损耗,并设总网损最优为配电网经济重构模型的目标函数,推导存在三相线路参数不对称和功率不平衡情况的配电网运行约束条件形式,建立起网络重构数学模型;
(3)使用最佳等距分段线性逼近的方法将所建立的网络重构数学模型中的二阶函数近似为线性函数,获得近似误差可调的基于线性潮流的混合整数线性网络重构模型并使用商业数学线性规划软件求解;
步骤(2)包括:
(a1)根据步骤(1)得到的配电网总损耗表达式,建立配电网网络重构数学模型,目标函数设为总网损最优,其中控制变量设为网络拓扑变量和线路潮流变量,决策变量是网络拓扑变量,网络重构模型的目标函数描述如下:
Figure FDA0002769460450000011
式中,N为网络总的节点数;
由于智能配电网存在三相线路参数不对称以及三相线路潮流不平衡的情况,因此网络重构模型的目标函数需要考虑支路电阻的不对称和支路有功无功功率的不一致情况;其中,线路电阻
Figure FDA0002769460450000012
为一个3×3的矩阵,具体描述如下:
Figure FDA0002769460450000013
式中
Figure FDA0002769460450000014
除主对角线外的其他元素均为三相之间的互阻抗,
Figure FDA0002769460450000015
Figure FDA0002769460450000016
支路有功
Figure FDA0002769460450000017
为一个3×1的矩阵,具体描述如下:
Figure FDA0002769460450000021
根据上述的各参数矩阵形式,可以写出智能配电网网络重构模型各相总损耗之和的目标函数;
(b1)该模型的约束条件包括节点注入功率约束,网络辐射状约束以及支路潮流上下限约束;
步骤(b1)包括的节点注入功率约束,网络辐射状约束以及支路潮流上下限约束具体如下:
配网网络重构模型的节点注入功率约束描述为:
Figure FDA0002769460450000022
Figure FDA0002769460450000023
式中,αij表示支路(i,j)的开合状态,αij=1表示该支路连通,αij=0表示该支路断开;N(i)表示网络中能够与节点i直接连接的所有节点集合;
Figure FDA0002769460450000024
表示节点i上的
Figure FDA0002769460450000025
相有功功率;
Figure FDA0002769460450000026
表示节点i上的负荷
Figure FDA0002769460450000027
相有功功率需求;
Figure FDA0002769460450000028
表示节点i上实际注入的
Figure FDA0002769460450000029
相有功功率;
Figure FDA00027694604500000210
表示节点i上的
Figure FDA00027694604500000211
相无功功率;
Figure FDA00027694604500000212
表示节点i上的负荷
Figure FDA00027694604500000213
相无功功率需求;
Figure FDA00027694604500000214
表示节点i上实际注入的
Figure FDA00027694604500000215
相无功功率;
由于智能配电网存在三相功率不平衡的情况,每个节点各相的节点注入功率都必须满足上面的功率平衡约束;其中DG的有功功率方向为正,负荷有功功率的方向为负;
配网网络重构模型的网络辐射状约束条件描述如下:
βijji=αij i=1,2,…,n j∈N(i)
Figure FDA00027694604500000216
β1j=0 j∈N(1)
βij∈{0,1} i=1,…,n j∈N(i)
0≤αij≤1 (i,j)∈E
式中,βij为生成树的父子节点关系变量,当节点j是节点i的父节点时βij=1,反之取0;其中第1式表明当支路(i,j)连接时,必有一个节点是另一个节点的父节点;第2式则表明节点i只存在一个父节点;第3式表明变电站1节点不是任何节点的父节点;每条支路的分段开关和联络开关为三相的,网络拓扑变量不用考虑各相不一致情况,因此当满足上述约束条件时,模型的拓扑优化结果即为辐射状网络;
配网网络重构模型的支路潮流上下限约束描述如下:
Figure FDA0002769460450000031
Figure FDA0002769460450000032
由于智能配电网存在三相功率不平衡的情况,因此约束条件考虑了各相的支路功率不相等的情况以及不同线路不同相的功率上下限也会存在差异的情况;
步骤(3)包括:
(a2)将智能配电网网络重构模型目标函数中的支路功率二次项
Figure FDA0002769460450000033
表示为y=h(x)形式的二次函数,使用最佳等距分段线性逼近法对该功率二次项进行一阶线性化近似表示;首先说明最佳等距分段中对于最小三相支路功率的分段距离查找原理:假设y=h(x)为需要线性化处理的三相不平衡智能配电网有功和无功非线性支路功率函数曲线,曲线的上方和下方分别存在一条三相误差功率上限曲线和三相误差功率下限曲线,分别可用
Figure FDA0002769460450000034
yh(x)=h(x)-δ表示,其中,δ为给定的最大三相功率逼近误差常数;x1点为三相误差功率下限曲线h(x)中斜率较大的一点,以该端点(x1,h(x1))作为起点做三相误差功率上限曲线的切线,可得到切点为
Figure FDA0002769460450000035
并与三相误差功率下限曲线相交于点(x2,h(x2)),该切线斜率为
Figure FDA0002769460450000036
接着以交点(x2,h(x2))为起点做三相误差功率上限曲线的切线,可得到切点为
Figure FDA0002769460450000037
并与三相误差功率下限曲线相交于点(x3,h(x3)),并得该切线斜率
Figure FDA0002769460450000038
反复执行上述操作直至
Figure FDA0002769460450000039
的终点,并得到一系列三相不平衡配电网的三相支路功率函数的分段直线;
对上述三相不平衡配电网支路功率曲线切线的斜率进行化简可得:
Figure FDA0002769460450000041
当x1和δ为已知量时,可以类推分别求得配电网重构模型中所有三相不平衡功率函数分段直线的端点值x2,x3,…,xl;基于上述三相功率端点值,可获得相邻功率端点的最小功率间距值C,以该间距值作为等距三相支路功率的分段标准,对三相不平衡配电网支路功率函数y=h(x)进行等距分段可得到新的三相不平衡支路功率分段直线端点,分别为
Figure FDA0002769460450000042
以此可获得等距分段线性逼近三相不平衡支路功率函数以表示原模型中的三相不平衡配电网功率二次函数;
经过上述对三相不平衡配电网模型中非线性功率项的最佳等距分段线性逼近,可获得配电网三相支路功率
Figure FDA0002769460450000043
的分段线性化近似表示式:
Figure FDA0002769460450000044
式中,Ωt,Ωs分别为
Figure FDA0002769460450000045
分段后的三相功率断面数集合;
Figure FDA0002769460450000046
分别为三相不平衡支路功率函数
Figure FDA0002769460450000047
各个断面上三相线性功率函数的斜率;
Figure FDA0002769460450000048
分别为
Figure FDA0002769460450000049
各断面上三相线性功率函数的取值,且
Figure FDA00027694604500000410
通过上述对三相不平衡配电网支路非线性功率项进行线性化近似,三相不平衡配电网重构模型的三相支路功率变量
Figure FDA00027694604500000411
可转化为由三相分段线性功率变量
Figure FDA00027694604500000412
构成的线性表达式:
Figure FDA00027694604500000413
(b2)通过线性化变换,三相不平衡配电网重构模型目标函数可转化为三相支路功率的线性函数,从而使得三相不平衡智能配电网网络重构问题转化为混合整数线性规划问题;由于三相线性支路功率变量
Figure FDA0002769460450000051
被限制为非负,而实际中当系统潮流与假定潮流正方向相反时,潮流值为负值,与设定不符,因此引入以下四个附加的三相不平衡潮流功率变量
Figure FDA0002769460450000052
并增加相应约束条件如下:
建立附加三相不平衡潮流功率变量与原三相支路功率变量
Figure FDA0002769460450000053
的关系约束如下:
Figure FDA0002769460450000054
附加三相变量自身约束如下:
Figure FDA0002769460450000055
Figure FDA0002769460450000056
Figure FDA0002769460450000057
各断面上三相不平衡线性功率函数的三相线性支路变量
Figure FDA0002769460450000058
Figure FDA0002769460450000059
的上下限约束:
Figure FDA00027694604500000510
(c2)使用商业数学线性规划求解器,即可对上述混合整数线性规划模型高效求解,通过设置分段后的断面数集合大小,可以调整网络重构线性模型的计算精度和效率。
2.根据权利要求1所述的基于线性模型的三相不平衡智能配电网网络重构方法,其特征在于,步骤(1)包括:
(a)对于辐射状运行的智能配电网,推导出其优化模型的Distflow支路潮流形式如下:
对于电网中节点j,有
Figure FDA00027694604500000511
Figure FDA0002769460450000061
式中,
Figure FDA0002769460450000062
为ABC三相模型;集合t(j)为电网中以j为末端节点的支路的首端节点集合;集合f(j)为电网中以j为首端节点的支路的末端节点集合;
Figure FDA0002769460450000063
为节点j的电压幅值;
Figure FDA0002769460450000064
分别为支路ij首端三相有功功率和无功功率;
Figure FDA0002769460450000065
为节点j的有功功率净注入值;
Figure FDA0002769460450000066
Figure FDA0002769460450000067
分别为节点j上所连接的DG有功功率和负荷有功功率;
Figure FDA0002769460450000068
为计及三相线路自阻抗和互阻抗的支路电阻;
(b)根据Distflow支路潮流形式,推导出智能配电网网络总损耗的表达形式如下:
Figure FDA0002769460450000069
式中,N为网络总的节点数;
实际配电网运行的节点电压标幺值一般在1附近,因此将节点电压近似为1,总损耗变为如下形式:
Figure FDA00027694604500000610
分母的二次项忽略后,总网损变为一个二阶函数的形式,便于进一步化简计算。
CN201711180647.7A 2017-11-23 2017-11-23 一种基于线性模型的三相不平衡智能配电网网络重构方法 Active CN108092284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711180647.7A CN108092284B (zh) 2017-11-23 2017-11-23 一种基于线性模型的三相不平衡智能配电网网络重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711180647.7A CN108092284B (zh) 2017-11-23 2017-11-23 一种基于线性模型的三相不平衡智能配电网网络重构方法

Publications (2)

Publication Number Publication Date
CN108092284A CN108092284A (zh) 2018-05-29
CN108092284B true CN108092284B (zh) 2021-01-05

Family

ID=62172982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711180647.7A Active CN108092284B (zh) 2017-11-23 2017-11-23 一种基于线性模型的三相不平衡智能配电网网络重构方法

Country Status (1)

Country Link
CN (1) CN108092284B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109361210A (zh) * 2018-10-24 2019-02-19 广东电网有限责任公司 一种基于辐射状结构建模的配电网检修优化方法
CN109473980B (zh) * 2018-11-30 2023-11-03 中国电力科学研究院有限公司 一种配电网的重构方法及系统
CN110061494B (zh) * 2019-04-09 2023-04-28 上海电力学院 一种考虑dg出力不确定性的三相不平衡配电网重构方法
CN111082470B (zh) * 2020-01-15 2022-09-02 合肥工业大学 含低风速分散式风电的配电网多目标动态鲁棒重构方法
CN112084626A (zh) * 2020-08-06 2020-12-15 国网浙江省电力有限公司嘉兴供电公司 基于分布式光伏接入的配电网无功补偿配置容量计算方法
CN114665478B (zh) * 2022-05-23 2022-10-11 国网江西省电力有限公司电力科学研究院 一种基于多目标深度强化学习的主动配电网重构方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310065A (zh) * 2013-06-25 2013-09-18 国家电网公司 计及分布式发电和储能单元的智能配网重构方法
CN105870912A (zh) * 2016-03-03 2016-08-17 广州华工弈高科技有限公司 主动配电网孤岛划分混合整数规划模型的建模方法
CN106300344A (zh) * 2016-09-07 2017-01-04 南京工程学院 一种基于改进二进制粒子群优化算法的配电网重构方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10296988B2 (en) * 2013-08-19 2019-05-21 Board Of Trustees Of Michigan State University Linear optimal power flow system and method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310065A (zh) * 2013-06-25 2013-09-18 国家电网公司 计及分布式发电和储能单元的智能配网重构方法
CN105870912A (zh) * 2016-03-03 2016-08-17 广州华工弈高科技有限公司 主动配电网孤岛划分混合整数规划模型的建模方法
CN106300344A (zh) * 2016-09-07 2017-01-04 南京工程学院 一种基于改进二进制粒子群优化算法的配电网重构方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"采用时段动态划分和分层优化策略的配电网重构";江东林等;《电网技术》;20120228;第36卷(第2期);全文 *

Also Published As

Publication number Publication date
CN108092284A (zh) 2018-05-29

Similar Documents

Publication Publication Date Title
CN108092284B (zh) 一种基于线性模型的三相不平衡智能配电网网络重构方法
CN104410066B (zh) 一种基于关联矩阵的矿井高压供电系统自动短路计算方法
CN107104442B (zh) 计及参数模糊性的含风电场电力系统概率潮流计算方法
CN104794541B (zh) 一种基于模拟退火和锥优化的配电网运行优化方法
TWI460436B (zh) 相量計算單元位置選取方法
CN112910012B (zh) 一种配电系统弹性提升评估方法、存储介质及计算设备
CN114140022A (zh) 一种多虚拟电厂分布式动态经济调度方法及系统
CN104143161A (zh) 一种基于广度优先搜索配电网回路的搜索方法
CN111917111B (zh) 用于配电网的分布式光伏电源接纳能力在线评估的方法、系统、设备及存储介质
CN111181164B (zh) 一种改进的主从分裂输配协同潮流计算方法及系统
Xiao et al. Optimal sizing and siting of soft open point for improving the three phase unbalance of the distribution network
CN109066674A (zh) 基于互补约束全光滑牛顿法的配电网重构方法
CN108306299B (zh) 一种配电网异步迭代分布式三相潮流计算方法及系统
CN113346484B (zh) 一种考虑暂态不确定性的配电网弹性提升方法及系统
Liu et al. Review of grid stability assessment based on AI and a new concept of converter-dominated power system state of stability assessment
CN112072692A (zh) 一种新能源发电场站的阻抗等值方法及装置
Wang et al. Analysis of network loss energy measurement based on machine learning
CN111834996B (zh) 电网线损计算方法和装置
CN116613746A (zh) 一种基于os2的分布式电网建模模型
CN107425519B (zh) 含分布式电源的三相配电网最大供电能力计算方法
Huang et al. Distribution network expansion planning considering distributed generation using probabilistic voltage constraints
CN110380418B (zh) 一种基于图分解的牛顿-拉夫逊潮流计算优化方法
Huang et al. A fast modeling method of distribution system reconfiguration based on mixed-integer second-order cone programming
Liu et al. A novel acceleration strategy for nl contingency screening in distribution system
Lin et al. A nested decomposition method and its application for coordinated operation of hierarchical electrical power grids

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