[发明专利]内嵌机电暂态仿真多端口等值网络元件参数快速计算方法

说明书

本发明涉及内嵌机电暂态仿真多端口等值网络元件参数快速计算方法，属于电力系统机电‑电磁暂态混合仿真分析技术领域；该方法包括：对电网某一正常运行状态或故障状态，确定电网拓扑，形成网络节点电压方程矩阵——导纳阵Y，对Y矩阵进行因子分解，得到存储因子分解的结果矩阵A；分别对矩阵A的上、下三角部分求解，得到解向量，再重构成压缩阻抗阵M并得到逆矩阵压缩导纳阵E；根据矩阵E的各元素，形成含参数的对应混合仿真分网接口实际拓扑的等值多端口网络；机电‑电磁混合仿真接口等值中实际拓扑对自阻抗、互阻抗，计算结束。本方法计算过程无简化，计算结果精确，大大提高了计算方法的效率和计算速度。

技术领域

本发明属于电力系统机电-电磁暂态混合仿真分析技术领域，特别涉及一种内嵌机电暂态仿真多端口等值网络元件参数快速计算方法。

背景技术

电力系统机电暂态仿真主要用于分析电力系统的稳定性，也用来分析电力系统某一正常运行状态下受到某种干扰后，能否经过一定的时间后回到原来的运行转态或者过渡到一个新的稳定运行状态的问题。目前面向大电网机电暂态仿真软件或程序，大多采用稀疏的导纳矩阵存储网络信息，基于此信息高效、快速、准确地计算出多端口等值网络元件参数，是电力系统机电-电磁暂态混合仿真接口等值必要的步骤和环节，也可用于电力系统仿真中网络等值简化。

利用接口处节点三相短路电流和稳态潮流求解等值阻抗也是一种常见的方法，但是在实际系统中其稳态潮流计算时采用的发电机等值电势与系统发生三相短路时采取的发电机等值暂态电势常常存在差异，大多数情况并不能严格保证等值电势的恒定不变，在求解等值阻抗时候有一定偏差。由于使用了基于潮流阻抗的对称故障信息，对于求解负序和零序等值阻抗，需要通过一系列的变换才能换算出负序和零序参数，计算速度较慢，且不能自动适应实际电网多种负载网架结构。

发明内容

本发明的目的是为克服已有技术的计算效率不足之处，提出了一种内嵌机电暂态仿真多端口等值网络元件参数快速计算方法，本方法计算过程无简化，计算结果精确，大大提高了计算方法的效率和计算速度。

本发明提出了一种内嵌机电暂态仿真多端口等值网络元件参数快速计算的方法，具体包括以下步骤：

1)记全电网节点个数为N，电力系统机电-电磁暂态混合仿真接口数量为n，在全电网N个节点中的序号记为{i₁,i₂,…,i_n}，n≤N；

2)对全电网某一正常运行状态或故障状态，确定电网拓扑，形成网络节点电压方程矩阵——导纳阵，记为Y，如式(1)所示:

式中，Y₁₁Y₂₂Y₃₃…Y_NN为自导纳、其它为互导纳；

3)对Y矩阵进行因子分解，得到存储因子分解的结果矩阵，记为A，如式(2)所示：

A＝L+D+U (2)

式中，A为矩阵Y的因子表矩阵，L为因子表矩阵A的下三角部分，D为因子表矩阵A的主对角元部分，U为因子表矩阵A的上三角部分；

4)由步骤3)得到的L、D、U分别同时并发n个线程，分别对基于矩阵L的前代、矩阵U的回代求解，如式(3)所示：

Yx_k＝b_k,k＝1,2,…,n (3)

式中，定义为第i_k个接口处元素为1.0、其它元素均为零的N维列向量，i_k∈1,2,…,n；x_k为相应的解向量；