CN107230982A - 一种微电网线性化潮流计算方法 - Google Patents
一种微电网线性化潮流计算方法 Download PDFInfo
- Publication number
- CN107230982A CN107230982A CN201710549492.3A CN201710549492A CN107230982A CN 107230982 A CN107230982 A CN 107230982A CN 201710549492 A CN201710549492 A CN 201710549492A CN 107230982 A CN107230982 A CN 107230982A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 23
- 239000013598 vector Substances 0.000 claims abstract description 44
- 238000011217 control strategy Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000005611 electricity Effects 0.000 claims description 3
- PXUQTDZNOHRWLI-OXUVVOBNSA-O malvidin 3-O-beta-D-glucoside Chemical group COC1=C(O)C(OC)=CC(C=2C(=CC=3C(O)=CC(O)=CC=3[O+]=2)O[C@H]2[C@@H]([C@@H](O)[C@H](O)[C@@H](CO)O2)O)=C1 PXUQTDZNOHRWLI-OXUVVOBNSA-O 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims 1
- 239000007924 injection Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect 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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H02J3/382—
-
- 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]
Abstract
一种微电网线性化潮流计算方法,包括如下步骤:获取微电网的网络参数;为各分布式电源DG成比例分配负荷功率,计算DG节点注入功率;依据微电网在不同控制策略下所要达到的平衡点,建立如下微电网线性化潮流方程:其中,ΔV为节点电压偏移相量所组成的矢量,S为各节点注入复功率矢量,VN为微电网无负荷接入时的各节点电压矢量;依据所述线性化潮流方程计算接入负荷后节点电压的偏移量、节点电压幅值、节点电压相角以及包含线损的DG输出功率。本发明对微电网控制系统所要达到的平衡点进行建模,大幅降低方程维度,并将非线性潮流方程线性化,提升计算结果的精确度。
Description
技术领域
本发明涉及微电网领域,尤其是一种微电网线性化潮流计算方法。
背景技术
微电网是有多种分布式电源(Distributed Generation,DG)、配电线路以及负荷构成的中低压电力系统,它的稳态模型是潮流方程。而潮流方程的计算是微电网的系统规划、经济调度、状态估计、电压控制等应用技术的基础。因此,对微电网潮流方程的建模、分析与计算,是微电网最重要的科研工作之一。
为大量消纳可再生能源,同时提高用户电能质量,大量微电网控制策略被提出,如今形成了以分层控制为标准,多控制器之间相互协调的体系。微电网内部的控制系统分为两层,一次控制是微电源控制器采用本地信息来迅速调节DG出力,保持系统电压和频率稳定,一般采用各种下垂控制策略;二次控制用来调节一次控制所产生的静差,根据通讯条件不同分为集中式、分散式和分布式三大类。现有的方法根据控制策略的不同将微电网潮流计算进行分类,并将控制策略的具体模型带入非线性潮流方程中考虑。然而,这一做法导致:
1)当逆变器数量较多时,将大幅增加方程的维数,而且需要对不同的控制策略进行单独建模,增加了问题的复杂程度。
2)非线性交流潮流方程的各状态变量高度耦合,需要依赖数值迭代方法求解。然而,迭代方法计算效率较低,而且收敛性难以保证,无法微电网的在线应用,如实时调度、静态安全分析等对计算时间的要求。
由于线性方程组具有成熟的求解方法,不存在收敛性问题,且可以求得未知量的解析表达式。因此,将非线性潮流方程线性化是一种常用方法,典型的是广泛应用于输电系统的直流潮流。然而,直流潮流假设了所有母线电压为常数、各相邻节点间的相位差很小、忽略线路电阻以及无功功率。这些假设严重不符合较高的中低压微电网工况,会导致较差。另外一个常用的潮流方程线性近似是配电网线性化支路潮流方程,即“LinDistFlow”,它假设系统中网络为辐射状拓扑,消去了电压相角,且忽略了系统线路损耗。因此,该方法对微电网的拓扑结构做出了限制,求解各节点电压相角,且在线路重载的情况下会产生较大的误差。
综上所述,现有的微电网非线性潮流计算方法难以保证计算的高效性和收敛的稳定性;而现有的线性化潮流计算方法无法保证结果的高精确性。
发明内容
本发明提供一种微电网线性化潮流计算方法,对微电网控制系统所要达到的平衡点进行建模,大幅降低方程维度,并将非线性潮流方程线性化,提升计算结果的精确度。
一种微电网线性化潮流计算方法,包括如下步骤:
获取微电网的网络参数;
为各分布式电源DG成比例分配负荷功率,计算DG节点注入功率,具体为
其中,分别为DG注入节点i的不包含线损的有功功率估计值和无功功率估计值,PLj、QLj分别为节点j的负荷功率,mi和ni分别为由DG的控制参数所决定的功率分配因子,NG、NL分别是微电网节点集合N={0,1,…,n}的电源节点和负荷节点子集合;
依据微电网在不同控制策略下所要达到的平衡点,建立如下微电网线性化潮流方程:
其中,ΔV为节点电压偏移相量所组成的矢量,S为各节点注入复功率矢量,VN为微电网无负荷接入时的各节点电压矢量;
依据所述线性化潮流方程计算接入负荷后节点电压的偏移量、节点电压幅值、节点电压相角以及包含线损的DG输出功率。
优选的,所述节点电压的偏移量的计算过程具体为:
解析所述线性化潮流方程,依据下列计算式计算负荷接入后的节点电压偏移相量:
其中,Δe和Δf分别为ΔV的实部和虚部矢量,eN和fN分别为VN的实部和虚部矢量,I为n维单位矩阵,H、N、J和K为子矩阵。
优选的,所述子矩阵H、N、J和K的计算过程具体为:
建立如下矩阵:
DV=diag(eN)2-diag(fN)2
其中,为无负荷节点电压幅值;
子矩阵H、N、J和K依据如下计算式获得:
其中,*表示矩阵的Hadamard积。
优选的,所述节点电压幅值和节点电压相角的计算过程具体为:
计算直角坐标系下的节点电压相量:V=VN+ΔV,式中,为各节点电压相量所组成的矢量;
依据节点电压相量计算节点电压幅值:E=EN+Δe,式中,为节点电压幅值矢量,EN=[|VN1|,...,|VNn|]T为无负荷电压幅值矢量;
计算节点电压相角:θ=θN+Δf,式中,θ节点电压相角矢量,θN=[θN1,...,θNn]T为无负荷电压相角矢量,其中各节点无负荷电压相角的计算公式为:
优选的,所述获取微电网的网络参数的步骤具体包括:
读取微电网的各网络参数,生成节点导纳矩阵Y,将节点导纳矩阵分解为
定义节点阻抗矩阵为
优选的,所述DG输出功率的计算过程如下:
取出矩阵YM中所有DG节点所在的行,并按照DG节点在前,负荷节点在后的排列顺序生成矩阵YG=[YGG,YGL],其中
DG输出功率的计算式为:SG=diag(VG)YGGVG+diag(VG)YGLVL;
其中,VG为DG节点电压矢量,VL为负荷节点电压矢量,SG为包含线损的DG输出功率矢量。
优选的,所述无负荷接入时各节点电压的计算式为:
VN=Zy0V0,式中,V0为PCC节点的电压值。
本发明的有益效果是,本发明根据微电网不同控制策略所要达到的共同目的,依据控制系统的稳态平衡点进行建模,无需考虑微电源控制器的内部动态,大幅降低方程维数,简化了建模。同时,将非线性潮流方程线性化,从而得到各未知量的解析表达式,仅需要求解线性方程组,即可得到,计算结果的精确度更高。
附图说明
图1为本发明一种实施例的微电网线性化潮流计算方法的流程图;
图2为本发明一种实施例的微电网的结构示意图;
图3为基础负荷下的电压幅值结果比较图;
图4为基础负荷下的电压相角结果比较图;
图5为2.5倍负荷下的电压幅值结果比较图;
图6为2.5倍负荷下的电压幅值结果比较图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。
本发明实施例提供一种微电网线性化潮流计算方法,包括如下步骤:
S101:获取微电网的网络参数;
微电网的网络参数是指为建立微电网线性化潮流方程所需的网络参数,包括微电网的节点数量、节点电压等,具体可以通过微电网内设置的各种传感器获得。
在获取到各网络参数后,可依据微电网的各网络参数,生成节点导纳矩阵Y,将节点导纳矩阵分解为
定义节点阻抗矩阵为
S102:为各分布式电源DG成比例分配负荷功率,计算DG节点注入功率;
将微电网不同控制策略所要达到的平衡点,即各DG成比例分配负荷功率,进行建模。微电网的一次、二次控制的目标在于将有功负荷、无功负荷在各DG之间成比例分配,有:
当不计线路损耗时,可以得出DG节点注入功率的估计值,具体为:
其中,分别为DG注入节点i的不包含线损的有功功率估计值和无功功率估计值,PLj、QLj分别为节点j的负荷功率,mi和ni分别为由DG的控制参数所决定的功率分配因子,NG、NL分别是微电网节点集合N={0,1,…,n}的电源节点和负荷节点子集合。
S103:依据微电网在不同控制策略下所要达到的平衡点,建立微电网线性化潮流方程;
在传统的高斯-赛德尔潮流计算方法中,非线性潮流方程可以写成如下隐式Z-bus形式:V=Zdiag(V)-1S+Zy0V0,此时,负荷接入后的节点电压偏移相量为ΔV=V-VN=Zdiag(VN+ΔV)-1S,式中,ΔV为节点电压偏移相量所组成的矢量;S为各节点注入复功率矢量。
当|ΔVi|<1时,该式可以展开为如下诺依曼级数形式:
忽略式中高阶项,得到如下形式的微电网线性化潮流方程:
其中,ΔV为节点电压偏移相量所组成的矢量,S为各节点注入复功率矢量,VN为微电网无负荷接入时的各节点电压矢量。
S104:依据所述线性化潮流方程计算接入负荷后节点电压的偏移量、节点电压幅值、节点电压相角以及包含线损的DG输出功率。
利用上述待求解参数与网络参数间的线性关系,可通过解析表达式方式求解,所述节点电压的偏移量的计算过程具体为:
解析所述线性化潮流方程,依据下列计算式计算负荷接入后的节点电压偏移相量:
其中,Δe和Δf分别为ΔV的实部和虚部矢量,eN和fN分别为VN的实部和虚部矢量,I为n维单位矩阵,H、N、J和K为子矩阵。
具体的,所述子矩阵H、N、J和K的计算过程具体为:
建立如下矩阵:
DV=diag(eN)2-diag(fN)2
其中,为无负荷节点电压幅值;
子矩阵H、N、J和K依据如下计算式获得:
其中,*表示矩阵的Hadamard积。
节点电压幅值和节点电压相角的计算过程具体为:
计算直角坐标系下的节点电压相量:V=VN+ΔV,式中,为各节点电压相量所组成的矢量;
依据节点电压相量计算节点电压幅值:E=EN+Δe,式中,为节点电压幅值矢量,EN=[|VN1|,...,|VNn|]T为无负荷电压幅值矢量;
计算节点电压相角:θ=θN+Δf,式中,θ节点电压相角矢量,θN=[θN1,...,θNn]T为无负荷电压相角矢量,其中各节点无负荷电压相角的计算公式为:
DG输出功率的计算过程如下:
取出矩阵YM中所有DG节点所在的行,并按照DG节点在前,负荷节点在后的排列顺序生成矩阵YG=[YGG,YGL],其中
DG输出功率的计算式为:SG=diag(VG)YGGVG+diag(VG)YGLVL;
其中,VG为DG节点电压矢量,VL为负荷节点电压矢量,SG为包含线损的DG输出功率矢量。
其中,上述实施例中的无负荷接入时各节点电压的计算方式可采用现有方法,在传统的高斯-赛德尔潮流计算方法中,非线性潮流方程可以写成如下隐式Z-bus形式
V=Zdiag(V)-1S+Zy0V0
此时,无负荷节点电压可由如下公式计算:VN=Zy0V0,式中,V0为PCC节点的电压值。
至此,所有潮流计算中的待求量均可以由上述解析表达式经简单的矩阵运算得出,且无需迭代,不存在收敛性问题。
本发明实施例可采用IEEE-123节点三相配电网算例,选取其中一相并加入分布式电源进行计算,如图2所示,网络与负荷参数来自IEEE下属配电系统分析委员会提供的标准算例系统集。
本发明提出方法得到的节点电压、相角与交流潮流方程、传统直流潮流方程、以及配电系统线性支路潮流方程在给定负荷以及重载条件下的对比如图3、图4、图5和图6所示。从图中的对比可知,本发明所提出的方法的计算精度均高于现有线性化潮流计算方法,而且在重载条件下依然可以给出高精度的计算结果。相对于精确模型依赖于数值迭代方法,本方法仅需要求解线性方程组,计算效率高且不存在收敛性问题。该方法对微电网在线监测、优化规划、状态估计以及控制算法设计等方面都有重要意义,具有广泛的推广应用前景。
以上内容是结合具体的实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换。
Claims (7)
1.一种微电网线性化潮流计算方法,其特征在于,包括如下步骤:
获取微电网的网络参数;
为各分布式电源DG成比例分配负荷功率,计算DG节点注入功率,具体为
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>L</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>G</mi>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>L</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>G</mi>
</msub>
<mo>;</mo>
</mrow>
其中,分别为DG注入节点i的不包含线损的有功功率估计值和无功功率估计值,PLj、QLj分别为节点j的负荷功率,mi和ni分别为由DG的控制参数所决定的功率分配因子,NG、NL分别是微电网节点集合N={0,1,…,n}的电源节点和负荷节点子集合;
依据微电网在不同控制策略下所要达到的平衡点,建立如下微电网线性化潮流方程:
<mrow>
<mi>&Delta;</mi>
<mi>V</mi>
<mo>=</mo>
<msub>
<mi>V</mi>
<mi>N</mi>
</msub>
<mo>+</mo>
<mi>Z</mi>
<mover>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>S</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<mi>Z</mi>
<mover>
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>V</mi>
<mo>)</mo>
</mrow>
<mi>S</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>;</mo>
</mrow>
其中,ΔV为节点电压偏移相量所组成的矢量,S为各节点注入复功率矢量,VN为微电网无负荷接入时的各节点电压矢量;
依据所述线性化潮流方程计算接入负荷后节点电压的偏移量、节点电压幅值、节点电压相角以及包含线损的DG输出功率。
2.根据权利要求1所述的方法,其特征在于,所述节点电压的偏移量的计算过程具体为:
解析所述线性化潮流方程,依据下列计算式计算负荷接入后的节点电压偏移相量:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>e</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>f</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>I</mi>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>H</mi>
</mtd>
<mtd>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>N</mi>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>H</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>J</mi>
</mtd>
<mtd>
<mi>K</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</mtd>
<mtd>
<mi>J</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Δe和Δf分别为ΔV的实部和虚部矢量,eN和fN分别为VN的实部和虚部矢量,I为n维单位矩阵,H、N、J和K为子矩阵。
3.根据权利要求2所述的方法,其特征在于,所述子矩阵H、N、J和K的计算过程具体为:
建立如下矩阵:
DV=diag(eN)2-diag(fN)2
其中,为无负荷节点电压幅值;
子矩阵H、N、J和K依据如下计算式获得:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>H</mi>
<mo>=</mo>
<msubsup>
<mi>RM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>V</mi>
</msub>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>R</mi>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>Q</mi>
<mo>*</mo>
<msub>
<mi>e</mi>
<mi>N</mi>
</msub>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>XM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>V</mi>
</msub>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<mi>X</mi>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mo>*</mo>
<msub>
<mi>e</mi>
<mi>N</mi>
</msub>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
1
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<msubsup>
<mi>XM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>V</mi>
</msub>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>X</mi>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>*</mo>
<msub>
<mi>e</mi>
<mi>N</mi>
</msub>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>RM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>V</mi>
</msub>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>R</mi>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mo>*</mo>
<msub>
<mi>e</mi>
<mi>N</mi>
</msub>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>J</mi>
<mo>=</mo>
<msubsup>
<mi>RM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>XM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>K</mi>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>RM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>Q</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>XM</mi>
<mi>V</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,*表示矩阵的Hadamard积。
4.根据权利要求1所述的方法,其特征在于,所述节点电压幅值和节点电压相角的计算过程具体为:
计算直角坐标系下的节点电压相量:V=VN+ΔV,式中,为各节点电压相量所组成的矢量;
依据节点电压相量计算节点电压幅值:E=EN+Δe,式中,为节点电压幅值矢量,EN=[|VN1|,...,|VNn|]T为无负荷电压幅值矢量;
计算节点电压相角:θ=θN+Δf,式中,θ节点电压相角矢量,θN=[θN1,...,θNn]T为无负荷电压相角矢量,其中各节点无负荷电压相角的计算公式为:
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>e</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mi>N</mi>
<mo>.</mo>
</mrow>
5.根据权利要求1所述的方法,其特征在于,所述获取微电网的网络参数的步骤具体包括:
读取微电网的各网络参数,生成节点导纳矩阵Y,将节点导纳矩阵分解为
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mn>00</mn>
</msub>
</mtd>
<mtd>
<msubsup>
<mi>y</mi>
<mn>0</mn>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mi>M</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
定义节点阻抗矩阵为
6.根据权利要求5所述的方法,其特征在于,所述DG输出功率的计算过程如下:
取出矩阵YM中所有DG节点所在的行,并按照DG节点在前,负荷节点在后的排列顺序生成矩阵YG=[YGG,YGL],其中
DG输出功率的计算式为:SG=diag(VG)YGGVG+diag(VG)YGLVL;
其中,VG为DG节点电压矢量,VL为负荷节点电压矢量,SG为包含线损的DG输出功率矢量。
7.根据权利要求1所述的方法,其特征在于,所述无负荷接入时各节点电压的计算式为:
VN=Zy0V0,式中,V0为PCC节点的电压值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549492.3A CN107230982B (zh) | 2017-07-07 | 2017-07-07 | 一种微电网线性化潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710549492.3A CN107230982B (zh) | 2017-07-07 | 2017-07-07 | 一种微电网线性化潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107230982A true CN107230982A (zh) | 2017-10-03 |
CN107230982B CN107230982B (zh) | 2021-02-02 |
Family
ID=59957016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710549492.3A Active CN107230982B (zh) | 2017-07-07 | 2017-07-07 | 一种微电网线性化潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107230982B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769201A (zh) * | 2017-11-23 | 2018-03-06 | 武汉大学 | 一种计及单一tcsc等效电抗参数的有功静态安全域构建方法 |
CN108183487A (zh) * | 2018-01-23 | 2018-06-19 | 中国石油大学(华东) | 一种基于线性数学模型的配电网潮流快速分析方法 |
CN109193664A (zh) * | 2018-09-12 | 2019-01-11 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1034593A1 (en) * | 1997-11-28 | 2000-09-13 | Abb Ab | Power flow control |
CN101034809A (zh) * | 2007-04-06 | 2007-09-12 | 清华大学 | 基于动态连续潮流的电力系统静态电压稳定的评估方法 |
CN101944742A (zh) * | 2010-08-30 | 2011-01-12 | 天津大学 | 改进的潮流可行解恢复方法 |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
CN103729537A (zh) * | 2013-09-12 | 2014-04-16 | 广西电网公司电力科学研究院 | 一种感知配电网实时态势的扩展锥方法 |
-
2017
- 2017-07-07 CN CN201710549492.3A patent/CN107230982B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1034593A1 (en) * | 1997-11-28 | 2000-09-13 | Abb Ab | Power flow control |
CN101034809A (zh) * | 2007-04-06 | 2007-09-12 | 清华大学 | 基于动态连续潮流的电力系统静态电压稳定的评估方法 |
CN101944742A (zh) * | 2010-08-30 | 2011-01-12 | 天津大学 | 改进的潮流可行解恢复方法 |
CN102570477A (zh) * | 2012-01-04 | 2012-07-11 | 河南省电力公司 | 基于线性化潮流算法模型的svc及tcsc混合布点规划方法 |
CN103729537A (zh) * | 2013-09-12 | 2014-04-16 | 广西电网公司电力科学研究院 | 一种感知配电网实时态势的扩展锥方法 |
Non-Patent Citations (1)
Title |
---|
SAIRAJ V DHOPLE: "Linear Approximations to AC Power Flow in Rectangular Coordinates", 《2015 53RD ANNUAL ALLERTON CONFERENCE》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769201A (zh) * | 2017-11-23 | 2018-03-06 | 武汉大学 | 一种计及单一tcsc等效电抗参数的有功静态安全域构建方法 |
CN108183487A (zh) * | 2018-01-23 | 2018-06-19 | 中国石油大学(华东) | 一种基于线性数学模型的配电网潮流快速分析方法 |
CN108183487B (zh) * | 2018-01-23 | 2022-03-15 | 中国石油大学(华东) | 一种基于线性数学模型的配电网潮流快速分析方法 |
CN109193664A (zh) * | 2018-09-12 | 2019-01-11 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流计算方法 |
WO2020052104A1 (zh) * | 2018-09-12 | 2020-03-19 | 国电南瑞科技股份有限公司 | 含dg并网的配电网潮流确定方法及计算机存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107230982B (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Happ | Optimal power dispatch | |
CN105514971A (zh) | 一种适用于各种运行模式微电网的潮流计算方法 | |
CN105048468B (zh) | 基于分布式计算的输配电网一体化电压稳定评估方法 | |
CN102931663B (zh) | 一种大规模交直流电力系统的动态无功优化方法 | |
CN105322546A (zh) | 交直流解耦混合潮流算法 | |
CN107230982A (zh) | 一种微电网线性化潮流计算方法 | |
CN104578045A (zh) | 独立直流微网智能功率分配方法 | |
CN108683191A (zh) | 一种下垂控制型孤岛微电网的三相潮流分析方法 | |
Gao et al. | A distributed coordinated economic droop control scheme for islanded AC microgrid considering communication system | |
Liu et al. | Remote voltage control using the holomorphic embedding load flow method | |
Liu et al. | The averaged-value model of a flexible power electronics based substation in hybrid AC/DC distribution systems | |
Ghaedi et al. | Transient stability improvement of power system with UPFC control by using transient energy function and sliding mode observer based on locally measurable information | |
Mok et al. | Tuning the fuzzy damping controller for UPFC through genetic algorithm with comparison to the gradient descent training | |
CN104734148A (zh) | 计及分布式电源的三相配电网连续潮流算法 | |
Feng et al. | Distributed hierarchical control for fast frequency restoration in VSG-controlled islanded microgrids | |
CN108718091B (zh) | 一种应用于主动配电网的三相极坐标系线性潮流计算方法 | |
Benedito et al. | Port-Hamiltonian based Optimal Power Flow algorithm for multi-terminal DC networks | |
CN104979840B (zh) | 一种主动配电网三相无功优化方法 | |
Dubyago | Mathematical description of interaction of the elements in the power network system | |
Xia et al. | Optimal secondary control of islanded AC microgrids with communication time-delay based on multi-agent deep reinforcement learning | |
CN106953328A (zh) | 一种柔性直流电网潮流线性化确定方法及系统 | |
Mahmud et al. | ANFISPID-based voltage regulation strategy for grid-tied renewable DG system with ESS | |
CN114329980A (zh) | 模块化多电平矩阵变换器的机电暂态建模方法及拓扑结构 | |
Fetzer | Development of a MATLAB/Simulink Framework for Phasor-Based Power System Simulation and Component Modeling Based on State Machines | |
Bian et al. | Identification and improvement of probabilistic voltage instability modes of power system with wind power integration |
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 |