一种电力系统潮流计算方法与流程
本发明涉及电力系统分析计算技术领域，具体涉及一种电力系统潮流计算方法。
背景技术：
电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网络结构确定整个系统的运行状态，如各母线上的电压(幅值及相角)、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础，通过电力系统潮流计算已比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。
基于牛顿-拉夫逊的潮流计算是电力系统分析计算的一种常用方法，该方法通常将负荷节点设置为pq节点，即负荷的有功功率和无功功率为恒定值。而在实际系统下，负荷的有功功率和无功功率随负荷节点电压值的变化而变化，因此采用恒定功率负荷模型进行潮流计算得到的结果并不符合实际。
传统的牛顿-拉夫逊在处理负荷节点时主要的不足在于：在计算功率的不平衡量时，将负荷功率与电压初值计算功率做差值，没有考虑到负荷的功率-电压特性。
技术实现要素：
针对传统的牛顿-拉夫逊在处理负荷节点时在计算功率的不平衡量时，将负荷功率与电压初值计算功率做差值，没有考虑到负荷的功率-电压特性，分析结果不准确的问题，本发明提供一种电力系统潮流计算方法。
本发明的技术方案是：
一种电力系统潮流计算方法，包括如下步骤：根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定zp负荷模型参数，在形成节点导纳矩阵时，将各负荷节点zp负荷模型恒阻抗分量的导纳值加入节点导纳矩阵对应对角元，在设定各节点电压初值后，采用zp负荷模型中恒功率分量功率值计算负荷节点恒功率分量功率的不平衡量和电压不平衡量，进行潮流计算。
优选地，所述的根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定zp负荷模型参数的步骤包括：
根据有载调压分接开关变化引起变电站低压侧的特征量计算zp负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值bp、bq：
其中，式中ua、ub为主变最近一次有载调压分接开关动作前后主变低压侧电压值，δp1为有载调压分接开关动作后低压侧有功功率瞬时变化值，δp2为功率瞬时变化后恢复到功率稳定的变化量，δq1为有载调压分接开关变化后低压侧无功功率瞬时变化值，δq2为功率瞬时变化后恢复到功率稳定的变化量，pi、qi为i节点负荷有功、无功功率。
优选地，根据zp负荷模型参数计算节点导纳矩阵的具体步骤包括：
设定n个节点的电力网络，其网络方程可表示为：
非负荷节点的对角元yii等于各个节点与该节点所连导纳的总和：
负荷节点i的对角元yii等于各个节点与该节点所连导纳的总和再加上zp负荷模型中恒阻抗分量的导纳值：
式中，i≠j，当j＝0时，zij为节点对地导纳，yli可由下式计算得出：
式中，ui为负荷节点电压值，pi、qi为负荷功率，bp、bq该节点zp负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值。
优选地，设定各节点电压初始值的步骤包括：
设定各节点电压初值实部
优选地，根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量的步骤包括：
将各节点电压初值或上次修正计算值代入下式ei、fi，求负荷节点恒功率分量功率的不平衡量δpi、δqi和电压不平衡量δui：
式中，(1-bp)pi、(1-bq)qi分别为去除zp负荷模型中恒阻抗分量负荷的负荷有功、无功功率，ei、fi分别为节点电压实部、虚部，gij为节点导纳矩阵第i行第j列元素实部，bij为节点导纳矩阵第i行第j列元素虚部。
优选地，计算各节点电压修正量及电压修正值的步骤具体包括：
将各节点电压初值代入，计算修正方程式系数矩阵元素，系数矩阵元素表达式如下：
当j＝i
当j≠i
rij＝sij＝0
计算各节点电压的修正量：
计算各节点电压修正值公式如下：
优选地，判断是否满足收敛条件具体为：
若不满足上式的收敛条件，则以各节点电压修正值作为初值，返回步骤：根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量进行下一轮迭代，ε为收敛精度，可根据实际要求设定。
从以上技术方案可以看出，本发明具有以下优点：实际工作过程中，由于负荷存在功率-电压特性，当系统电压变化时负荷功率随之变化，本方法解决了传统的方法当电压变化时，负荷功率保持不变的问题，本发明计算潮流结果负荷节点功率随电压降低而降低，较为符合实际情况，为后续的电力系统的分析提供可靠的参考。
此外，本发明设计原理可靠，结构简单，具有非常广泛的应用前景。
由此可见，本发明与现有技术相比，具有突出的实质性特点和显著地进步，其实施的有益效果也是显而易见的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种电力系统稳定性分析方法的示意性流程图。
图2是本发明实施例的算例模型图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
如图1所示，一种电力系统潮流计算方法，包括如下步骤：
s1：确定zp负荷模型参数；
本步骤中，利用各负荷节点变电站与当前时间最近一次有载调压分接开关变化确定zp负荷模型参数，步骤具体包括：
根据有载调压分接开关变化引起变电站低压侧的特征量计算zp负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值bp、bq：
其中，式中ua、ub为主变最近一次有载调压分接开关动作前后主变低压侧电压值，δp1为有载调压分接开关动作后低压侧有功功率瞬时变化值，δp2为功率瞬时变化后恢复到功率稳定的变化量，δq1为有载调压分接开关变化后低压侧无功功率瞬时变化值，δq2为功率瞬时变化后恢复到功率稳定的变化量，pi、qi为i节点负荷有功、无功功率。
s2：根据zp负荷模型参数计算节点导纳矩阵；
具体步骤包括：
设定n个节点的电力网络，其网络方程可表示为：
非负荷节点的对角元yii等于各个节点与该节点所连导纳的总和：
负荷节点i的对角元yii等于各个节点与该节点所连导纳的总和再加上zp负荷模型中恒阻抗分量的导纳值：
式中，i≠j，当j＝0时，zij为节点对地导纳，yli可由下式计算得出：
式中，ui为负荷节点电压值，pi、qi为负荷功率，bp、bq该节点zp负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值。
s3：设定各节点电压初始值；设定各节点电压初值实部
s4：根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量；
将各节点电压初值或上次修正计算值代入下式ei、fi，求负荷节点恒功率分量功率的不平衡量δpi、δqi和电压不平衡量δui：
式中，(1-bp)pi、(1-bq)qi分别为去除zp负荷模型中恒阻抗分量负荷的负荷有功、无功功率，ei、fi分别为节点电压实部、虚部，gij为节点导纳矩阵第i行第j列元素实部，bij为节点导纳矩阵第i行第j列元素虚部。
s5：将各节点电压初值代入，求计算修正方程式系数矩阵元素；
将各节点电压初值代入，计算修正方程式系数矩阵元素，系数矩阵元素表达式如下：
当j＝i
当j≠i
rij＝sij＝0
s6：计算各节点电压修正量；
计算各节点电压的修正量：
s7：计算电压修正值修正各节点电压初值；计算各节点电压修正值公式如下：
s8：判断是否满足收敛条件
若不满足上式的收敛条件，则各节点电压修正值作为初值，返回步骤s4：根据电压修正值计算负荷节点恒功率分量功率不平衡量和电压不平衡量进行下一轮迭代，ε为收敛精度，可根据实际要求设定。
s9：输出计算结果。
实施例1
本发明实施例如图2所示选用ieee9节点系统来进行测试，对比分析现有技术和本发明潮流计算结果。
表1给出了各负荷节点zp负荷模型参数；表2给出了bus1发电机电压为100％额定值时，现有技术潮流计算功率、电压值；表3给出了bus1发电机电压为90％额定值时，现有技术和本发明潮流计算功率、电压值比较。
表1各负荷节点zp负荷模型参数
表2bus1发电机电压为额定值100％时潮流计算电压值比较
表3bus1发电机电压为额定值90％时潮流计算电压值比较
由表2和表3可看出，降低bus1发电机电压后，系统各节点电压均下降。由于负荷存在功率-电压特性，当电压降低时，负荷功率响应降低。比较表2、表3现有技术计算潮流结果，负荷节点bus5、bus6、bus8有功、无功功率并没有随电压降低而降低，本发明计算潮流结果负荷节点功率随电压降低而降低，较为符合实际情况。
尽管通过参考附图并结合优选实施例的方式对本发明进行了详细描述，但本发明并不限于此。在不脱离本发明的精神和实质的前提下，本领域普通技术人员可以对本发明的实施例进行各种等效的修改或替换，而这些修改或替换都应在本发明的涵盖范围内/任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
- 还没有人留言评论。精彩留言会获得点赞！