[发明专利]一种电力系统潮流计算的导纳矩阵计算方法

说明书

本发明公开了一种电力系统潮流计算的导纳矩阵计算方法，包括以下步骤：读支路首节点号数组I、末节点号数组J、电阻数组R、电抗数组X、对地电纳数组B、变压器变比数组K；定义AS和AE为n×l阶稀疏矩阵；计算支路π形等值电路串联支路导纳数组；计算节点支路关联矩阵；由串联支路形成初始导纳矩阵；计算有向支路首节点和末节点的对地导纳数组；形成各节点对地导纳数组并修正导纳矩阵。本发明提出的方法在Matlab平台实现，便于科研人员使用Matlab的各种工具和函数对计算结果进行测试和分析。本发明提出的导纳矩阵计算采用矩阵运算和复数运算，减少了程序代码，简化了编程，程序更加清晰；矩阵运算也大大提高了计算速度。

技术领域

本发明涉及一种电力系统潮流计算的导纳矩阵计算方法，特别是一种适合研究目的使用的基于Matlab的潮流计算的导纳矩阵计算方法。

背景技术

电力系统潮流计算是研究电力系统稳态运行的一项基本计算，它根据给定的运行条件和网络结构确定整个网络的运行状态。潮流计算也是电力系统其他分析的基础，如安全分析、暂态稳定分析等都要用到潮流计算。潮流计算是电力系统分析的基本分析工具，科研人员经常以潮流计算为基础进行进一步地研究。实用的商业软件采用C语言等高级编程语言编写且使用稀疏矩阵技术和节点优化编号等高级技术。这些技术虽然能大幅度提高潮流计算的速度、降低内存占用量，但编程非常麻烦且难以修改和维护，不易增加新的功能，因而不适合科研人员用于研究目的使用。

Matlab软件以矩阵为最基本的数据单位，可以方便地处理各种矩阵和向量运算，也可以很方便自然地处理复数类型，其指令表达式与数学中常用的形式很接近，还有大量常见实用的函数，给编程带来很大便利。Matlab软件简单易用、代码短小易操作，易于编程和调试，计算功能强大，同时还具有非常强大的可视化图形处理和交互式功能，为科学研究以及工程应用提供了一种高效的编程工具，目前已经成为许多科学领域的基本工具和首选平台，在各种科学和工程计算领域得到了广泛的应用。为了适应越来越多的科研人员需要在Matlab平台上以潮流计算为基础进行进一步地研究的需求，迫切需要一种基于Matlab软件的易于编程、修改和调试的潮流计算方法。

目前常用的牛顿法和快速分解法潮流计算方法都是以节点电压法为基础的，都需要求取节点导纳矩阵。

节点导纳矩阵为：

式中，Y_ik为节点导纳矩阵元素，当下标i≠k时，为节点i和节点k之间的互导纳，当下标i＝k时，为节点i的自导纳；n为节点数。

可以采用追加支路法计算导纳矩阵元素，即依次扫描所有支路，每扫描一条支路在原有导纳矩阵元素基础上增加一个导纳增量。由于输电线路和变压器支路都属于支路，通常把输电线路和变压器支路都作为支路数据统一输入。作为区分，变压器支路的非标准变比侧的节点号加个负号。

输电线路采用如图2所示的π形等值电路，设第m条支路的串联支路导纳为

y_m＝1/z_m＝1/(r_m+jx_m) (2)

式中，r_m，x_m，z_m分别为输电线路等值电路的电阻、电抗和阻抗。

追加第m条输电线路时，导纳矩阵元素的计算公式如下：

式中，下标i、j分别表示支路的首节点号i_m和末节点号j_m去掉负号后的节点号，b_m为输电线路等值电路的对地电纳，符号“←”表示右端计算结果赋值给左端变量。