一种下垂控制型孤岛微电网的三相潮流分析方法技术

本发明专利技术公开了一种下垂控制型孤岛微电网的三相潮流分析方法，首先建立孤岛微电网中各下垂控制型DG、负荷及线路元件的三相模型并计算系统的节点导纳矩阵Y，列写下垂控制型孤岛微电网的节点功率平衡方程，得到三相潮流方程组，然后推导三相潮流方程组在采用牛顿法求解时的修正方程式及雅可比矩阵J的计算公式，最后计算下垂节点功率及支路功率，输出潮流计算结果，本发明专利技术解决了现有技术中存在的下垂控制型孤岛微电网的三相潮流方法中信赖域法潮流计算量大而常规牛顿法难收敛的问题。

【技术实现步骤摘要】
一种下垂控制型孤岛微电网的三相潮流分析方法
本专利技术属于电力系统分析
，具体涉及一种下垂控制型孤岛微电网的三相潮流分析方法。
技术介绍
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算，是电网技术研究的一个重要领域，是进行电网稳定性分析、规划和运行研究的基础。微电网作为一种灵活高效利用分布式能源的形式，有并网和孤岛运行两种形式。孤岛微电网有主从控制和对等控制两种结构。主从控制结构的微电网内有主电源提供电压频率支撑，其运行机理与传统电网类似，在潮流计算中，由主电源作为平衡节点，其他DG一般均采用恒功率控制，可处理为PQ或PV节点。而对等控制结构中，由多个可控型分布式电源(DistributedGeneration,DG)共同参与电压、频率调节和控制，这些DG通常采用下垂控制方法，系统的运行机理与传统电网有明显区别。因此，传统的潮流计算方法对其不再适用，有必要结合下垂控制的特点对此类孤岛微电网的潮流计算问题进行研究。目前，对微电网潮流计算的研究，大多未考虑对等控制策略或分散下垂控制策略下的孤岛微电网的潮流计算。在考虑分散下垂控制策略时，目前文献中有两类计算方法：信赖域法及牛顿法。采用信域法求解下垂控制的孤岛微电网的潮流时，在迭代过程必须反复重新计算海森矩阵及其逆矩阵，导致在大规模潮流中计算量过大，计算速度慢于牛顿拉夫逊法；而牛顿法用在计算下垂控制型孤岛微电网潮流尤其是三相不对称的孤岛微电网的潮流计算时，收敛性较差，常出现难收敛的情况。
技术实现思路
本专利技术的目的是提供一种下垂控制型孤岛微电网的三相潮流分析方法，解决了现有技术中存在的下垂控制型孤岛微电网的三相潮流方法中信赖域法潮流计算量大而常规牛顿法难收敛的问题。本专利技术所采用的技术方案是，一种下垂控制型孤岛微电网的三相潮流分析方法，具体按照以下步骤实施：步骤1、建立孤岛微电网中各下垂控制型DG、负荷及线路元件的三相模型并计算系统的节点导纳矩阵Y；步骤2、列写下垂控制型孤岛微电网的节点功率平衡方程，得到三相潮流方程组f(x)＝0；步骤3、推导步骤2的三相潮流方程组在采用牛顿法求解时的修正方程式及雅可比矩阵J的计算公式；步骤4、设定迭代次数k＝0，设定PV节点电压幅值的初值为其给定值，其它各节点电压幅值的初值为标幺值1，节点电压相角初值为0，系统频率初值为标幺值1，并得到对应待求未知量的对应初值x(0)；步骤5、将x(k)代入f(x)的表达式，计算得到f(x(k))，并判断其是否满足收敛判据max|f(x(k))|＜ε，ε为要求的计算精度；若满足收敛判据，则跳转至步骤10；否则进行步骤6；步骤6、计算雅可比矩阵J(k)，计算修正量Δx(k)；步骤7、利用潮流方程的泰勒级数的高阶展开信息确定最优乘子步骤8、在每步迭代时以作为修正量进行修正；步骤9、迭代次数k＝k+1；转至步骤5；步骤10、计算下垂节点功率及支路功率，输出潮流计算结果。本专利技术的特点还在于，步骤1具体按照以下步骤实施：步骤1.1、首先建立下垂控制型微电源的三相模型：设微电网中所有节点构成的集合为B，而BPQ、BPV及BDroop分别表示微电网中PQ节点、PV节点及Droop下垂节点的集合，当Droop节点按有功/无功限制运行，按PQ节点计算，三相不对称孤岛微电网中，对含有下垂控制型DG的节点有：其中：分别为下垂控制型DG输出的三相总有功功率和无功功率；分别为其单相有功功率和无功功率值，其中，m＝{a,b,c}；mPi、nQi为下垂节点有功和无功的等效调节系数；ω为系统角频率；和ω0为下垂控制DG的电压和角频率设定值；为下垂节点i的正序电压幅值，设和分别为为节点i的m相电压幅值和相角，其中，m＝{a,b,c}，则有：同时，下垂节点的三相电压幅值相等，波形对称，即：步骤1.2、建立三相负荷的静态特性模型如下：式中，分别为m相负荷的有功功率、无功功率；分别为m相负荷的有功功率额定值、无功功率额定值；ULi0为节点i的电压的额定值；ω为系统的角频率，ωL0为角频率的额定值；和分别是m相负荷有功功率静态电压特性系数，和分别是m相负荷无功功率静态电压特性系数，和分别是m相负荷有功功率和无功功率的静态频率特性系数；步骤1.3、建立三相线路的等值电路：三相不对称孤岛微电网中，节点i和j之间支路的线路用三相阻抗矩阵表示：其中，为节点i和j之间支路的三相阻抗矩阵，为节点i和j之间支路的三相之间的阻抗，其中，m＝{a,b,c},n＝{a,b,c}；步骤1.4、建立节点导纳矩阵Y。步骤2具体按照以下步骤实施：列写下垂控制型孤岛微电网的节点功率平衡方程，对任意节点i，设注入节点i的p相有功功率和无功功率分别记为Pip和其中，p＝{a,b,c}，则其功率平衡方程表示为：其中：ΔPip、分别为节点i的p相有功和无功功率不平衡量；Pip、分别为节点i的p相有功和无功注入功率；和分别为节点导纳矩阵中节点i的p相和节点k的m(m＝{a,b,c})相的导纳阵元素的实部和虚部；为节点i的p相和节点k的m相的电压的相角差。对于每个PQ节点，其功率不平衡量表达式对应有6个方程：其中，和分别为节点i接入的除下垂控制型DG之外如风力发电机组或光伏单元等DG的p相有功功率和无功功率，若无风电机组或光伏单元，则取值为零；对于每个PV节点，由于节点电压恒定，且不考虑无功功率不平衡量方程，其功率平衡方程组为3个方程：每个下垂节点的功率平衡方程组将包含12个等式约束：对于每个下垂节点都有12个未知量，用表示下垂节点的各相电压幅值，用表示a相相角，并将b和c相相角分别用和代入，同时将和用下垂特性关系式表示并消去和则下垂节点的未知量化简为6个其节点功率平衡方程为：假设取节点1的相角为参考相角，即为方便表述，设表示除下垂节点外其它节点的三相有功功率不平衡量列向量，表示PQ节点的三相有功功率不平衡量列向量，和分别表示下垂节点三相有功和无功功率不平衡量列向量，并记设表示除参考节点及下垂节点外其它节点的三相电压相角列向量，表示PQ节点的电压幅值列向量，和分别表示各下垂控制节点的a相电压相角和幅值列向量，及分别表示下垂节点的a相和b相有功及无功功率列向量，并记则式(6)～(10)的三相潮流方程组写成以下形式：将式(11)孤岛微电网潮流方程进一步统一表示为如下形式：f(x)＝0(12)式中，步骤3具体按照以下步骤实施：采用牛顿法求解上述三相潮流方程组，推导其修正方程式，如下所示：其中，Δx1Δx2Δx3为变量x1、x2及x3的修正量；雅可比矩阵J的各个子块计算公式如下：式(14)～式(18)中各雅可比矩阵子块J11、J12、J21及J22中的H、N、M、L分块的各元素定义与传统三相潮流计算中对应分块元素的定义本质相同，只是每个子块对应的相别和节点集合不同，具体来说，H、M分块与传统三相潮流计算雅可比矩阵Jorig元素表达式没有区别，而N、L分块需在Jorig元素表达式的基础上，对其主对角元素需进行修改调整，设和分别为Jorig中N和L分块的任一对角元素，和分别为孤岛微电网雅可比矩阵J中N和L分块的对应对角元素，用和分别表示其修改叠加量，即有：J11中的E11和F11分块和J21中的E21和F21分块均本文档来自技高网...

【技术保护点】
1.一种下垂控制型孤岛微电网的三相潮流分析方法，其特征在于，具体按照以下步骤实施：步骤1、建立孤岛微电网中各下垂控制型DG、负荷及线路元件的三相模型并计算系统的节点导纳矩阵Y；步骤2、列写下垂控制型孤岛微电网的节点功率平衡方程，得到三相潮流方程组f(x)＝0；步骤3、推导步骤2的三相潮流方程组在采用牛顿法求解时的修正方程式及雅可比矩阵J的计算公式；步骤4、设定迭代次数k＝0，设定PV节点电压幅值的初值为其给定值，其它各节点电压幅值的初值为标幺值1，节点电压相角初值为0，系统频率初值为标幺值1，并得到对应待求未知量的对应初值x(0)；步骤5、将x(k)代入f(x)的表达式，计算得到f(x(k))，并判断其是否满足收敛判据max|f(x(k))|＜ε，ε为要求的计算精度；若满足收敛判据，则跳转至步骤10；否则进行步骤6；步骤6、计算雅可比矩阵J(k)，计算修正量Δx(k)；步骤7、利用潮流方程的泰勒级数的高阶展开信息确定最优乘子

【技术特征摘要】
1.一种下垂控制型孤岛微电网的三相潮流分析方法，其特征在于，具体按照以下步骤实施：步骤1、建立孤岛微电网中各下垂控制型DG、负荷及线路元件的三相模型并计算系统的节点导纳矩阵Y；步骤2、列写下垂控制型孤岛微电网的节点功率平衡方程，得到三相潮流方程组f(x)＝0；步骤3、推导步骤2的三相潮流方程组在采用牛顿法求解时的修正方程式及雅可比矩阵J的计算公式；步骤4、设定迭代次数k＝0，设定PV节点电压幅值的初值为其给定值，其它各节点电压幅值的初值为标幺值1，节点电压相角初值为0，系统频率初值为标幺值1，并得到对应待求未知量的对应初值x(0)；步骤5、将x(k)代入f(x)的表达式，计算得到f(x(k))，并判断其是否满足收敛判据max|f(x(k))|＜ε，ε为要求的计算精度；若满足收敛判据，则跳转至步骤10；否则进行步骤6；步骤6、计算雅可比矩阵J(k)，计算修正量Δx(k)；步骤7、利用潮流方程的泰勒级数的高阶展开信息确定最优乘子步骤8、在每步迭代时以作为修正量进行修正；步骤9、迭代次数k＝k+1；转至步骤5；步骤10、计算下垂节点功率及支路功率，输出潮流计算结果。2.根据权利要求1所述的一种下垂控制型孤岛微电网的三相潮流分析方法，其特征在于，所述步骤1具体按照以下步骤实施：步骤1.1、首先建立下垂控制型微电源的三相模型：设微电网中所有节点构成的集合为B，而BPQ、BPV及BDroop分别表示微电网中PQ节点、PV节点及Droop下垂节点的集合，当Droop节点按有功/无功限制运行，按PQ节点计算，三相不对称孤岛微电网中，对含有下垂控制型DG的节点有：其中：分别为下垂控制型DG输出的三相总有功功率和无功功率；分别为其单相有功功率和无功功率值，其中，m＝{a,b,c}；mPi、nQi为下垂节点有功和无功的等效调节系数；ω为系统角频率；和ω0为下垂控制DG的电压和角频率设定值；为下垂节点i的正序电压幅值，设和分别为为节点i的m相电压幅值和相角，其中，m＝{a,b,c}，则有：同时，下垂节点的三相电压幅值相等，波形对称，即：步骤1.2、建立三相负荷的静态特性模型如下：式中，分别为m相负荷的有功功率、无功功率；分别为m相负荷的有功功率额定值、无功功率额定值；ULi0为节点i的电压的额定值；ω为系统的角频率，ωL0为角频率的额定值；和分别是m相负荷有功功率静态电压特性系数，和分别是m相负荷无功功率静态电压特性系数，和分别是m相负荷有功功率和无功功率的静态频率特性系数；步骤1.3、建立三相线路的等值电路：三相不对称孤岛微电网中，节点i和j之间支路的线路用三相阻抗矩阵表示：其中，为节点i和j之间支路的三相阻抗矩阵，为节点i和j之间支路的三相之间的阻抗，其中，m＝{a,b,c},n＝{a,b,c}。3.根据权利要求2所述的一种下垂控制型孤岛微电网的三相潮流分析方法，其特征在于，所述步骤2具体按照以下步骤实施：列写下垂控制型孤岛微电网的节点功率平衡方程，对任意节点i，设注入节点i的p相有功功率和无功功率分别记为Pip和其中，p＝{a,b,c}，则其功率平衡方程表示为：其中：ΔPip、分别为节点i的p相有功和无功功率不平衡量；Pip、分别为节点i的p相有功和无功注入功率；和分别为节点导纳矩阵中节点i的p相和节点k的m(m＝{a,b,c})相的导纳阵元素的实部和虚部；为节点i的p相和节点k的m相的电压的相角差；对于每个PQ节点，其功率不平衡量表达式对应有6个方程：其中，和分别为节点i接入的除下垂控制型DG之外如风力发电机组或光伏单元等DG的p相有功功率和无功功率，若无风电机组或光伏单元，则取值为零；对于每个PV节点，由于节点电压恒定，且不考虑无功功率不平衡量方程，其功率平衡方程组为3个方程：每个下垂节点的功率平衡方程组将包含12个等式约束：对于每个下垂节点都有12个未知量，用表示下垂节点的各相...

【专利技术属性】
技术研发人员：潘忠美，刘健，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：陕西,61