CN106169750B - 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 - Google Patents
一种基于二阶锥松弛的主动配电网最大供电能力计算方法 Download PDFInfo
- Publication number
- CN106169750B CN106169750B CN201610618114.1A CN201610618114A CN106169750B CN 106169750 B CN106169750 B CN 106169750B CN 201610618114 A CN201610618114 A CN 201610618114A CN 106169750 B CN106169750 B CN 106169750B
- Authority
- CN
- China
- Prior art keywords
- distribution network
- node
- power distribution
- active power
- constraint
- 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
- 238000000205 computational method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000011156 evaluation Methods 0.000 claims abstract description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 abstract description 5
- 230000001131 transforming effect Effects 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009365 direct transmission Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于二阶锥松弛的主动配电网最大供电能力计算方法,属于电力系统优化评估领域。本发明方法包括:考虑所有不同主变故障的场景,引入新用于表示不同故障场景的参数,建立了一个基于二阶锥松弛的主动配电网最大供电能力计算模型,其目标函数为最大化主动配电网中总的所带负荷量减去加权后的网络功率损耗,同时需要满足电网运行的各项技术约束及设定所述主动配电网中描述潮流变量关系的二阶锥约束。通过该模型的求解,可以获得该主动配电网的最大供电能力评估结果,即满足配变N‑1约束情况下所带总负荷上限,以及取负荷上限时相应的各个节点负荷值。本发明方法建立的评估模型十分精细,符合实际情况,且易于求解,具有很强的实用性。
Description
技术领域
本发明涉及一种基于二阶锥松弛的主动配电网最大供电能力计算方法,属于电力系统优化评估领域。
背景技术
主动配电网主变压器(简称主变)故障属于最严重的故障,安全的配电网要求主变故障时所有失电负荷均能被转供。变电站中的任意变压器发生故障时,其下所带负荷如果有被同站其它变压器转供、或者通过馈线间联络线转供的能力,则称为该配网系统满足配变N-1安全约束。主动配网系统在满足配变N-1约束的情况下所带总负荷会有一个上限,即配电网最大供电能力;最大供电能力是衡量配电网供电安全程度十分重要的指标,与电网运行的安全性息息相关,需要得到合理准确计算。配电网通常呈辐射状运行,对于主变引出的每条馈线支路,主变就可视为其根节点,而配电网的最大供电能力即为所有根节点能够向馈线支路输出的最大功率总和。
然而,传统的配电网供电能力计算方法所建模型十分粗糙,如式(1-1)所示:
其中,目标函数为最大化各个主变的载荷,目标函数中的TSC表示配电网的供电能力,Ri为载荷率,Ti为主变容量;而任意一个主变所带的负荷,都能由馈线间的通道sij转移(满足通道容量上限Cij),满足了主变N-1的约束,其中和分别为同变电站和不同变电站的有直接传输通道的馈线集合(同站内转供时主变能够以系数k过载)。所以该传统模型(1-1)仅为考虑馈线间联络容量约束、主变约束,以及满足任意主变故障后失电负荷均能被转供的条件下,最大化配电网的带负荷能力,没有涉及支路潮流和节点电压的具体约束,也没有考虑馈线内具有的可重构性,这将使得供电能力评估结果过于乐观或悲观而偏离实际情况。
在数学上,二阶锥规划的标准形式如式(1-2)所示:
其中,x∈RN为决策变量;系数常量包括b∈RM、c∈RN和AM×N∈RM×N;K表示如下形式的二阶锥(1-3)或旋转二阶锥(1-4):
a)二阶锥
b)旋转二阶锥
二阶锥规划可以视作线性规划的推广,本质上属于一种凸规划,因此具有解的最优性和计算高效性。在很多情况下,一些非凸优化问题可以通过二阶锥松弛处理,转化为二阶锥规划问题以便求解。
发明内容
本发明的目的是为了克服已有技术的不足之处,提出了一种基于二阶锥松弛的主动配电网最大供电能力计算方法。该方法考虑了主动配电网中的安全运行约束和网络重构能力,能够对其最大供电能力进行准确评估,并通过二阶锥松弛技术的应用使得该问题转化成为易于求解的形式。
本发明提出的一种基于二阶锥松弛的主动配电网最大供电能力计算方法,其特征在于,该方法包括以下步骤:
1)考虑所有不同主变故障的场景,引入参数f=0,1,2,...,Ntrans用于不同故障场景的表示;其中Ntrans为已知的主变压器个数;f为场景参数,等于0代表没有任何主变故障的正常运行状态,等于1表示主变1故障的场景,等于2代表主变2故障的场景,以此类推;
2)建立所述主动配电网供电能力计算的目标函数如式(1)所示:
该目标函数为最大化主动配电网中总的所带负荷量减去加权后的网络功率损耗;其中,ΦN为主动配电网中除根节点外所有节点的集合;LP,i表示节点i处的有功负荷变量;Rij为已知的支路ij电阻值;为故障场景f下支路ij流过的电流幅值的平方变量;φ为人为设定的权重系数,取值范围为(0,10);
3)设定所述主动配电网的辐射状约束如式(2)所示:
其中,为故障场景f下描述支路ij开断状态的二元变量,等于0表示该支路处于断开状态,等于1表示该支路处于连接状态;Nnode为主动配电网中除根节点外其他所有节点的总个数,是已知参数;B表示与故障变压器直接相连的支路集合;
4)设定所述主动配电网的功率平衡约束如式(3)所示:
其中,dij为描述支路ij潮流流向的已知二元参数,取1表示从节点j流向节点i,取-1表示从节点i流向节点j;N(i)表示与节点i相连的所有节点集合;和分别为故障场景f下支路ij以节点i为终端的有功和无功潮流变量;LP,i和LQ,i分别为节点i处的有功和无功负荷变量;为节点i处已知的负荷功率系数;
5)设定所述主动配电网中分布式电源出力约束如式(4)所示:
其中,ΦDG为该主动配电网中分布式电源节点的集合;和分别为节点i处分布式电源已知的最大有功和无功出力;
6)设定所述主动配电网中每条支路传输的功率容量约束如式(5)所示:
其中,为支路ij已知的视在功率上限值;
7)设定所述主动配电网中每个节点的电压安全约束如式(6)所示:
其中,为故障场景f下节点i的电压幅值变量;为故障场景f下节点i的电压幅值的平方变量;和分别为节点i已知的电压幅值平方的下限值和上限值;
8)设定所述主动配电网中网络损耗约束如式(7)所示
其中,为故障场景f下支路ij的电流幅值变量;为故障场景f下支路ij的电流幅值平方变量;Rij和Xij分别为支路ij已知的电阻值和电抗值;
9)设定所述主动配电网中变压器容量约束如式(8)所示
其中,ΦR为主动配电网中的根节点集合;为已知的根节点电压幅值的平方;为根节点i处变压器的已知容量值;
10)设定所述主动配电网中潮流方程约束如式(9)所示
其中,M0为人为设定的很大的正数,取值范围为[1000,10000];
11)设定所述主动配电网中描述潮流变量关系的二阶锥约束如式(10)所示
式(10)是潮流变量关系约束采用二阶锥松弛后的形式;其中,||·||2表示取2-范数内积运算;
12)基于步骤3)到步骤11)建立的主动配电网运行约束条件,求解步骤2)中的目标函数式(1),从而获得该主动配电网的总供电能力评估结果,即满足配变N-1约束情况下所带总负荷上限,以及取负荷上限时相应的各个节点负荷值。
本发明提出的这种一种基于二阶锥松弛的主动配电网最大供电能力计算方法,其优点是:
1、本发明方法考虑了支路潮流和节点电压的具体约束,以及馈线内的网络重构能力,获得的主动配电网供电能力评估结果更为准确,符合实际情况。
2、本发明方法除了能获得主动配电网最大供电能力的结果外,还能得到该结果下具体各节点所带的负荷量,即最大供电能力下的负荷分布情况。
3、本发明方法所建立的优化模型为一个凸的二阶锥规划问题,易于求解,具有很强的实用性,可以应用许多商业优化求解器直接进行模型的求解。
具体实施方式
本发明提出的一种基于二阶锥松弛的主动配电网最大供电能力计算方法,进一步说明如下:
本方法包括以下步骤:
1)考虑所有主变依次发生故障的场景,引入参数f=0,1,2,...,Ntrans用于不同故障场景的表示;其中Ntrans为已知的主变压器个数;f为场景参数,等于0代表没有任何主变故障的正常运行状态,等于1表示主变1故障的场景,等于2代表主变2故障的场景,以此类推;
2)建立所述主动配电网供电能力计算的目标函数如式(1)所示:
该目标函数为最大化主动配电网中总的所带负荷量减去加权后的网络功率损耗;其中,ΦN为主动配电网中除根节点外所有节点的集合;LP,i表示节点i处的有功负荷变量;Rij为已知的支路ij电阻值;为故障场景f下支路ij流过的电流幅值的平方变量;φ为人为设定的权重系数,取值范围为(0,10),取9为典型值;
3)设定所述主动配电网的辐射状约束如式(2)所示:
该式描述了主动配电网的辐射状运行约束,以保证最终网络拓扑结构中没有环状回路;其中,为故障场景f下描述支路ij开断状态的二元变量,等于0表示该支路处于断开状态,等于1表示该支路处于连接状态;Nnode为主动配电网中除根节点外其他所有节点的总个数,是已知参数;B表示与故障变压器直接相连的支路集合;
4)设定所述主动配电网的功率平衡约束如式(3)所示:
该式描述了主动配电网中每个节点的功率平衡约束;其中,dij为描述支路ij潮流流向的已知二元参数,取1表示从节点j流向节点i,取-1表示从节点i流向节点j;N(i)表示与节点i相连的所有节点集合;和分别为故障场景f下支路ij以节点i为终端的有功和无功潮流变量;LP,i和LQ,i分别为节点i处的有功和无功负荷变量;为节点i处已知的负荷功率系数;
5)设定所述主动配电网中分布式电源出力约束如式(4)所示:
该式描述了主动配电网中分布式电源的有功和无功出力约束;其中,ΦDG为该主动配电网中分布式电源节点的集合;和分别为节点i处分布式电源已知的最大有功和无功出力;
6)设定所述主动配电网中每条支路传输的功率容量约束如式(5)所示:
该式描述了主动配电网中每条支路的功率传输容量约束,以支路视在功率容量作为其限制值;其中,为支路ij已知的视在功率上限值;
7)设定所述主动配电网中每个节点的电压安全约束如式(6)所示:
该式描述了主动配电网中每个节点的电压幅值的上下限约束;其中,为故障场景f下节点i的电压幅值变量;为故障场景f下节点i的电压幅值的平方变量;和分别为节点i已知的电压幅值平方的下限值和上限值;
8)设定所述主动配电网中网络损耗约束如式(7)所示
该式描述了主动配电网中支路的有功和无功功率损耗约束;其中,为故障场景f下支路ij的电流幅值变量;为故障场景f下支路ij的电流幅值平方变量;Rij和Xij分别为支路ij已知的电阻值和电抗值;
9)设定所述主动配电网中变压器容量约束如式(8)所示
该式描述了主动配电网中的变压器功率容量约束;其中,ΦR为主动配电网中的根节点集合;为已知的根节点电压幅值的平方;为根节点i处变压器的已知容量值;
10)设定所述主动配电网中潮流方程约束如式(9)所示
该式描述了主动配电网中的潮流方程约束,即限定了支路的有功、无功潮流与支路两端节点电压之间的关系;其中,M0为人为设定的很大的正数,取值范围为[1000,10000],取5000为典型值;
11)设定所述主动配电网中描述潮流变量关系的二阶锥约束如式(10)所示
式(10)是潮流变量关系约束采用二阶锥松弛后的形式;其中,||·||2表示取2-范数内积运算;
12)基于步骤3)到步骤11)建立的主动配电网运行约束条件,求解步骤2)中的目标函数式(1),从而获得该主动配电网的总供电能力评估结果,即满足配变N-1约束情况下所带总负荷上限,以及取负荷上限时相应的各个节点负荷值。
Claims (1)
1.一种基于二阶锥松弛的主动配电网最大供电能力计算方法,其特征在于,本方法包括以下步骤:
1)考虑所有不同主变故障的场景,引入参数f=0,1,2,...,Ntrans用于不同故障场景的表示;其中Ntrans为已知的主变压器个数;f为场景参数,等于0代表没有任何主变故障的正常运行状态,等于1表示主变1故障的场景,等于2代表主变2故障的场景,以此类推;
2)建立所述主动配电网供电能力计算的目标函数如式(1)所示:
该目标函数为最大化主动配电网中总的所带负荷量减去加权后的网络功率损耗;其中,ΦN为主动配电网中除根节点外所有节点的集合;LP,i表示节点i处的有功负荷变量;Rij为已知的支路ij电阻值;为故障场景f下支路ij流过的电流幅值的平方变量;φ为人为设定的权重系数,取值范围为(0,10);
3)设定所述主动配电网的辐射状约束如式(2)所示:
其中,为故障场景f下描述支路ij开断状态的二元变量,等于0表示该支路处于断开状态,等于1表示该支路处于连接状态;Nnode为主动配电网中除根节点外其他所有节点的总个数,是已知参数;B表示与故障变压器直接相连的支路集合;
4)设定所述主动配电网的功率平衡约束如式(3)所示:
其中,dij为描述支路ij潮流流向的已知二元参数,取1表示从节点j流向节点i,取-1表示从节点i流向节点j;N(i)表示与节点i相连的所有节点集合;和分别为故障场景f下支路ij以节点i为终端的有功和无功潮流变量;LP,i和LQ,i分别为节点i处的有功和无功负荷变量;为节点i处已知的负荷功率系数;
5)设定所述主动配电网中分布式电源出力约束如式(4)所示:
其中,ΦDG为该主动配电网中分布式电源节点的集合;和分别为节点i处分布式电源已知的最大有功和无功出力;
6)设定所述主动配电网中每条支路传输的功率容量约束如式(5)所示:
其中,为支路ij已知的视在功率上限值;
7)设定所述主动配电网中每个节点的电压安全约束如式(6)所示:
其中,为故障场景f下节点i的电压幅值变量;为故障场景f下节点i的电压幅值的平方变量;和分别为节点i已知的电压幅值平方的下限值和上限值;
8)设定所述主动配电网中网络损耗约束如式(7)所示
其中,为故障场景f下支路ij的电流幅值变量;为故障场景f下支路ij的电流幅值平方变量;Rij和Xij分别为支路ij已知的电阻值和电抗值;
9)设定所述主动配电网中变压器容量约束如式(8)所示
其中,ΦR为主动配电网中的根节点集合;为已知的根节点电压幅值的平方;为根节点i处变压器的已知容量值;
10)设定所述主动配电网中潮流方程约束如式(9)所示
其中,M0为人为设定的很大的正数,取值范围为[1000,10000];
11)设定所述主动配电网中描述潮流变量关系的二阶锥约束如式(10)所示
式(10)是潮流变量关系约束采用二阶锥松弛后的形式;其中,||·||2表示取2-范数内积运算;
12)基于步骤3)到步骤11)建立的主动配电网运行约束条件,求解步骤2)中的目标函数式(1),从而获得该主动配电网的总供电能力评估结果,即满足配变N-1约束情况下所带总负荷上限,以及取负荷上限时相应的各个节点负荷值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610618114.1A CN106169750B (zh) | 2016-07-29 | 2016-07-29 | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610618114.1A CN106169750B (zh) | 2016-07-29 | 2016-07-29 | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106169750A CN106169750A (zh) | 2016-11-30 |
CN106169750B true CN106169750B (zh) | 2018-10-19 |
Family
ID=58065050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610618114.1A Active CN106169750B (zh) | 2016-07-29 | 2016-07-29 | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106169750B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107968397B (zh) * | 2017-11-27 | 2020-09-25 | 国网江西省电力有限公司经济技术研究院 | 一种考虑运行随机性的配电网传输能力计算方法 |
CN108306342B (zh) * | 2018-01-19 | 2021-04-02 | 广东电网有限责任公司佛山供电局 | 含dg离散功率因数和设备动作限制的配网运行控制方法 |
CN108233369B (zh) * | 2018-01-30 | 2020-12-04 | 广东电网有限责任公司韶关供电局 | 一种预想事故条件下的主动配电网载荷能力安全评估方法 |
CN108493924B (zh) * | 2018-03-06 | 2020-08-04 | 深圳供电局有限公司 | 配电网供电能力评估方法和系统、计算机存储介质及设备 |
CN110391660B (zh) * | 2018-04-17 | 2023-01-24 | 中国电力科学研究院有限公司 | 一种提升配电网供电能力的网络重构方法及装置 |
CN109066654B (zh) * | 2018-08-13 | 2020-09-18 | 深圳供电局有限公司 | 基于混合整数线性规划的配电网最大供电能力评估方法 |
CN109256772B (zh) * | 2018-10-16 | 2020-09-01 | 清华大学 | 一种电力系统输电网最优开断的凸优化求解方法 |
CN109861232A (zh) * | 2019-02-22 | 2019-06-07 | 国家电网有限公司 | 一种基于二阶锥松弛方法的配电网动态重构降损方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10271695A (ja) * | 1997-03-27 | 1998-10-09 | Mitsubishi Motors Corp | バッテリ残存容量検出装置及びハイブリッド電気自動車の発電制御装置 |
CN102025153B (zh) * | 2010-12-15 | 2012-11-14 | 天津大学 | 一种中压配电网最大供电能力精确计算方法 |
CN102611106B (zh) * | 2012-03-31 | 2014-04-16 | 凯里供电局 | 一种环型供电中压配电网的最大供电能力评估方法 |
CN102769287B (zh) * | 2012-07-04 | 2014-09-17 | 天津大学 | 一种配电网最大供电能力计算方法 |
CN103279656B (zh) * | 2013-05-20 | 2016-03-09 | 天津大学 | 基于潮流的配电网最大供电能力计算方法 |
CN104485662A (zh) * | 2014-11-24 | 2015-04-01 | 国家电网公司 | 一种变电站并列运行时的配电网最大供电能力的方法 |
CN105337275B (zh) * | 2015-10-26 | 2017-12-19 | 天津大学 | 一种基于配电系统可靠性的中压配电网供电能力评估方法 |
-
2016
- 2016-07-29 CN CN201610618114.1A patent/CN106169750B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106169750A (zh) | 2016-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106169750B (zh) | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 | |
CN107464048B (zh) | 一种基于研究态的日前计划安全校核方法 | |
Han et al. | An automated impedance estimation method in low-voltage distribution network for coordinated voltage regulation | |
CN105762826B (zh) | 一种含vsc-hvdc交直流系统最大输电能力计算方法 | |
CN111680084A (zh) | 一种基于多源数据融合的线损计算方法 | |
CN105826918A (zh) | 电网理论线损的计算方法和系统 | |
CN109038551B (zh) | 交直流并联输电通道功率优化分配的混合整数凸规划方法 | |
CN109829560B (zh) | 一种配电网可再生能源发电集群接入规划方法 | |
CN104600695A (zh) | 基于在线状态估计和实时调度计划的趋势潮流计算方法 | |
CN104348158A (zh) | 地县一体化自动电压控制中协调控制区的自动生成方法 | |
CN108233383A (zh) | 一种交直流互联电网的经济调度方法 | |
CN104158199B (zh) | 对电力系统实时状态进行无功电压优化控制的系统和方法 | |
CN103746388A (zh) | 配电网无功电压三级协调控制方法 | |
CN108134394A (zh) | 一种考虑分布式电源影响的优化减载方法 | |
CN107276098A (zh) | 一种配电网多级协调控制系统及方法 | |
CN106099914A (zh) | 一种基于供电能力的主动配电网安全等级划分方法 | |
CN109377020B (zh) | 一种考虑配电网负荷转供能力的输电网规划方法 | |
CN103617569B (zh) | 基于质量评价的配网状态估计数据辨识与修正方法 | |
Neagu et al. | The optimal operation of active distribution networks with smart systems | |
CN106505576B (zh) | 一种基于Link-path的城市高压配电网负载转供方法 | |
CN107769187A (zh) | 一种适用于电压型直流潮流控制器的优化配置方法 | |
Singh et al. | A smart volt-var optimization engine for energy distribution system | |
CN105354645B (zh) | 一种220kV变电站负荷转带方法 | |
CN107221926A (zh) | 一种基于tsc的馈线接入用户容量计算方法 | |
Magham et al. | Voltage profile improvement in a microgrid including PV units using genetic algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |