一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法，属于电力系统分析计算领域。本发明专利技术是建立在一种仅存贮节点导纳矩阵Y上下三角非零元素构成的数据文件的基础之上。这种存贮方式大大减少了Y阵数据文件的存贮单元以及数据文件的存贮时间，在形成雅可比矩阵J过程中大大减少了对Y阵数据文件的读取时间，并且利用Y阵与J阵元素结构相似、非零元素对应的特点，在计算J阵元素和节点电流及节点功率时免去了大量零元素的无效计算以及非零元素的判断时间，从而大大加快了形成J阵的速度。用本发明专利技术方法对IEEE-30、-57、-118节点系统进行验算，与传统方法相比，形成J阵的计算速度上可提高约61％～92％。
【技术实现步骤摘要】
本专利技术属于电力系统分析计算领域。
技术介绍
在大型电力系统的潮流计算中，经常需根据节点导纳矩阵Y形成雅可比矩阵J。但 由于Y阵中具有大量零元素，如不考虑利用Y阵元素的稀疏性、J阵与Y阵元素结构的相似 性以及J阵元素自身结构的特点，会导致J阵形成过程中由于Y阵数据文件较大而严重影 响J阵的形成效率以及大量无效J阵元素计算，使J阵形成时间过长。 传统方法形成J阵的过程中，有以下几点不足： (1) 一般形成J阵所用Y阵数组的形式是Y(n，2n)，而Y(n，2n)数组由于大量零元 素的存贮不但需要较大的存贮单元，且数据文件的读取速度也不高。一般情况下，读取Y 阵数据文件所需的时间远远高于用Y阵元素形成J阵的时间，且读取Y(n，2n)数组数据文 件的时间在整个J阵形成过程中所占比例较大。所以快速形成J阵的关键在于Y阵数组的 存贮和读取方式。 (2)由于Y阵和J阵元素的稀疏性结构极为相似，如果不考虑Y阵元素稀疏性直接 根据Y阵形成J阵，则对J阵有大量零元素的无效计算，导致形成J阵的计算效率极低。 (3)如果根据Y阵元素的稀疏性形成J阵，在利用Y(n，2n)数组形成J阵时可对 Y(n，2n)数组的每个元素进行非零判断（可同时判断实部和虚部元素或仅判断虚部元素， 后者速度快，本文采用后者），仅计算Yu辛0时的J ^子阵（H ^ I、Mp 或H ^ Si：j)兀素，且同时利用同组兀素之间的关系，即利用Hi』=L i：j，_M i：j和RS u = 0的 关系。虽然计算效率比（2)略有提高，如对IEEE-118节点系统可提高约27%，大量的判断 语句仍然使得形成J阵时的计算效率不高。 (4)如果考虑Y阵元素的稀疏性进行按坐标存贮、按顺序存贮或按链表存贮，尽管 可以省去不少贮存单元，但在这些存贮方式中，Y阵元素与J阵元素没有明显的对应关系， 不但数据文件的读取所需时间较长，也不便于直接用来计算J阵元素。 (5)如果利用Y阵元素的对称稀疏性只存贮Y阵上三角的非零元素，则在计算形成 J阵元素的过程中，根据上三角的非零元素按对称性获得下三角的非零元素时很不方便，且 计算程序中过多的下标转换及赋值语句也使得计算J阵元素速度等受到影响。此外，在节 点电流Ii、节点功率Pi、Qi等连续求和的计算过程中，下三角非零元素的快速获取也会受到 影响。 (6)如果不能同时两行/次计算J阵元素、利用同组元素之间的关系以及0 时有Jij# 〇和子Jji# 〇阵的关系，形成J阵的计算效率仍然较低。
技术实现思路
为了克服上述现有技术的不足，本专利技术在一种快速存贮及读取电力系统稀疏矩阵 数据方法的基础上，提出一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵 的方法。 本专利技术方法是建立在一种快速存贮及读取电力系统稀疏矩阵数据方法的基础之 上。该数据存贮方法是将Y阵中的节点分为主节点和子节点。仅存贮主节点的行号、节点 数Si (包括主节点和子节点）、主节点以及相连的具有非零元素子节点的列号和参数（自导 纳和互导纳）。 这种存贮方式简单明了，易于数据的检索、修改和计算应用。与传统的不考虑元素 稀疏性的Y阵元素存贮形式相比，大大减少了 Y阵数据的存贮单元，以及Y阵数据文件的存 贮和读取时间。与按坐标存贮、按顺序存贮、按链表存贮等等方案相比，也较多地减少了存 贮空间，并提高了数据的读取速度，且后续数据的计算处理等更为方便。与考虑对称性仅存 贮上三角非零元素的方式相比，几乎不需增加任何存贮单元，且获得下三角的非零元素更 方便、快捷。 这种Y阵元素存贮方式完全反映了网络结构的特征和数组元素的稀疏性，并且利 用了 Y阵与J阵元素结构相似、非零元素对应的特点，在形成J阵和计算1111、〇1时，可快速 获取Y阵上下三角的所有非零元素，免去对非零元素的判断以及大量零元素的无效计算， 从而大大的减少了形成J阵的时间。 由于程序的模块化，在建立J阵之前，首先建立快速存贮及读取电力系统稀疏矩 阵数据的数据文件。该数据文件将系统中Y阵的节点分为主节点和子节点。假设主节点 最多连接7个子节点，则行号和节点数以及一行中主节点和子节点的列号及参数（实部、 虚部）共需26个存贮单元。因此，可定义Y阵数据文件的数组为Y(n，26)，将其每行分为3 组，分别存贮该行的行号、与行号对应的节点数、主节点和所有非零元素子节点的列号及参 数，具体方式如下。 第1组："主节点行号i"，位于第1列，为主节点对应的行号，为检验数据之用； 第2组："节点数S/'，位于第2列，为主节点及与该主节点连接的子节点数之和，S 值由程序自动累加以保证快速读取对应主节点和子节点的参数； 第3组："主节点、子节点列号j及参数"，位于第3~26列；按递增顺序存贮主节 点及与该主节点连接的子节点的列号j，分别位于第3*k列；相应主节点或子节点的自导纳 和互导纳（均分实部、虚部）分别位于3*k+l、3*k+2列，k的取值为：1~S i; 其元素的存贮方式如下表所示： 如果主节点最多连接的子节点数超过或少于7个，则可根据超过的个数将 Y(n，26)数组改为Y(n，26±3m)数组，m为超过或少于7个子节点后的子节点个数。对 Y(n，26)形式数据文件的读取时间远远低于对Y(n，2n)数据文件的读取时间。一般情况下， 出现子节点数超过7个的情况极少。即使超过7个，对Y(n，26+3m)形式数据文件的读取时 间也仍然远远低于对Y(n，2n)数据文件的读取时间。 在形成Y(n，26)数据文件的基础上，再应用本专利技术提出的基于稀疏矩阵技术快速 形成电力系统潮流计算中J阵的方法。本专利技术方法利用了 Y阵与J阵元素结构相似、非零 元素对应、以及&子阵同组元素之间的关系等特点，在形成J阵时可以免去对大量Y阵元 素的非零判断以及对J阵元素的无效计算，从而大大加快形成J阵的速度。 本专利技术是通过以下技术方案实现的，主要包括以下步骤： 步骤1 :打开数据文件，读取Y阵数据文件到Y(n，26)数组； 与读取Y(n，2n)数组形式的数据文件相比，读取Y(n，26)数组形式的数据文件可 节省了大量的时间。另外根据Y(n，26)数组中的节点数&，还能进一步改善对Y(n，26)数 组元素的读取效率。 步骤2 :根据Y(n，26)数组计算各个节点的Ipi、、或AP^AQi; (1)由于Ipi、Iqi、APp AQi的计算过程中也涉及到大量零元素，而Y(n，26)数组 仅存贮上下三角的非零元素，因此可直接根据Y (n，26)数组计算各个节点的Ipi、、或A P p A Qp不但可免去对所有元素的非零判断，还可省去对所有Y阵零元素的无效计算。 (2)Y(n，26)数组中主节点的行号便于对数据的检索和检验，节点数51可保证更有 效地减少对Y(n，26)数组不必要的循环和计算。 步骤3 :根据Y(n，26)数组和&子阵中元素的关系初步形成J阵； (1)根据Y(n，26)数组中的每行元素对应两行/次形成J阵中的元素、并利用Ju 子阵元素化=L u，-M ^和R S 0的关系仅计算J阵中基数行的非零元素，即 仅当前第1页1 2 本文档来自技高网...
【技术保护点】
一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法，其特征是包括以下步骤：步骤1：打开数据文件，读取Y阵数据文件到Y(n,26)数组；步骤2：根据Y(n,26)数组计算各个节点的有功和无功电流Ipi、Iqi或功率ΔPi、ΔQi；步骤3：根据Y(n,26)数组和Jij子阵中元素的关系初步形成J阵；(1)根据Y(n,26)数组中的每行元素对应两行/次形成J阵中的元素、并利用Jij子阵元素Hij＝Lij，Nij＝‑Mij和Rij＝Sij＝0的关系计算J阵中基数行的非零元素，而根据Jij子阵元素的关系得到J阵偶数行的非零元素；(2)通过Y(n,26)数组中的节点数Si确定J阵中每2行非零元素的计算个数；步骤4：根据步骤2计算出的Ipi、Iqi，修正J阵中所有的对角元素，形成完整的J阵；所述的Y(n,26)数组是：其每行分为3组，分别存贮该行的行号、与行号对应的节点数、主节点和所有非零元素子节点的列号及参数，具体方式如下：第1组：“主节点行号i”，位于第1列，为主节点对应的行号，为检验数据之用；第2组：“节点数Si”，位于第2列，为主节点及与该主节点连接的子节点数之和，S值由程序自动累加以保证快速读取对应主节点和子节点的参数；第3组：“主节点、子节点列号j及参数”，位于第3～26列；按递增顺序存贮主节点及与该主节点连接的子节点的列号j，分别位于第3*k列；相应主节点或子节点的自导纳和互导纳分别位于3*k+1、3*k+2列，k的取值为：1～Si。...
【技术特征摘要】
【专利技术属性】
技术研发人员：陈恳，邵尉哲，万新儒，彭丽君，
申请(专利权)人：南昌大学，
类型：发明
国别省市：江西;36
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。