CN105703359A - 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 - Google Patents
一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 Download PDFInfo
- Publication number
- CN105703359A CN105703359A CN201610130735.5A CN201610130735A CN105703359A CN 105703359 A CN105703359 A CN 105703359A CN 201610130735 A CN201610130735 A CN 201610130735A CN 105703359 A CN105703359 A CN 105703359A
- Authority
- CN
- China
- Prior art keywords
- factor table
- nonzero element
- row
- array
- sparse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 107
- 238000004364 calculation method Methods 0.000 title claims abstract description 46
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 26
- 238000006467 substitution reaction Methods 0.000 claims abstract description 20
- 238000005516 engineering process Methods 0.000 claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 230000003068 static effect Effects 0.000 claims description 7
- 238000012804 iterative process Methods 0.000 claims description 6
- 238000003491 array Methods 0.000 claims description 4
- 238000012856 packing Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000003860 storage Methods 0.000 description 20
- 239000011159 matrix material Substances 0.000 description 6
- 230000003993 interaction Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007935 neutral effect Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Complex Calculations (AREA)
Abstract
一种对称稀疏因子表法在直角坐标PQ分解法潮流计算中的应用，属于电力系统分析计算领域。本发明从给定结构的数据文件中读取数据，大大提高了数据读取速度和Ipi、Iqi或ΔPi、ΔQi的计算速度；利用对称稀疏技术快速形成因子表；记录其上三角非零元素的列号及其个数；将对称稀疏技术用于后续的快速前代及回代计算；在有功、无功迭代中引入塞德尔方式加快潮流计算速度。如对IEEE-118系统，本发明与传统方法相比，读取数据文件、形成因子表、有功及无功迭代的时间及潮流计算总时间分别为后者的6.12％、6.15％、11.86％和6.16％。且系统节点数越多，本发明的优势越大。
Description
技术领域
本发明属于电力系统分析计算领域，涉及直角坐标PQ分解法潮流计算。
背景技术
牛顿-拉夫逊法(牛顿法)是电力系统潮流计算中最常用的方法，按计算过程可分直角坐标和极坐标牛顿法，简化后可分别派生出直角坐标和极坐标PQ分解法。从计算过程上来看，极坐标牛顿法比直角坐标牛顿法修正方程式的方程数要少，雅可比矩阵元素的种类和数量也少，有时潮流迭代次数也少一次，因此极坐标比直角坐标牛顿法应用似乎更为广泛。同理，极坐标PQ分解法得到了广泛应用，而几乎无文献介绍直角坐标PQ分解法。实际上，由于直角坐标潮流计算中不涉及三角函数，因此在考虑稀疏性和直角坐标潮流计算系数矩阵结构特点的情况下其存贮单元数、潮流迭代速度等指标并不亚于极坐标法潮流计算。
电力系统中PQ分解法潮流计算中涉及到多个对称稀疏矩阵应用。若不考虑稀疏性，大量零元素的存贮、读取及计算过程中会导致存贮空间的极大浪费、读写数据文件及计算过程时间过长。因此，结合直角坐标PQ分解法系数矩阵的特点、考虑稀疏技术对元素的存贮方式及形成因子表的过程不但可大量节省存贮单元，还可大大减少数据文件的存贮、读取及计算过程的时间。
传统直角坐标PQ分解法潮流计算需3个数组，导纳矩阵Y用于计算节点电流(Ipi、Iqi)和节点功率(ΔPi、ΔQi)，系数矩阵B′、B″用于求取节点电压的虚部增量Δfi和实部增量Δei。传统方法中对3个数组数据的存贮、读取、应用等存在以下不足：
(1)B′、B″阵元素的获取方式不合理。
B′、B″阵元素的构成对潮流收敛性或收敛速度影响较大，因此其构成应各不相同。如果B′、B″阵元素直接取自Y阵元素的虚部，或如果仅B′、B″阵元素相同、而不同于Y阵元素的虚部，都可能影响PQ分解法潮流的迭代次数或收敛性。
(2)Y、B′、B″阵元素的存贮和读取方式不合理。
设系统节点数为n，PQ节点数为m，传统直角坐标PQ分解法中对应的数组分别为Y(n,2n)、B′(n-1,n-1)、B″(n-1,n-1)，这种存贮方式会造成存贮空间的极大浪费和读写数据文件及计算时间过长，且用Y(n,2n)数组计算Ipi、Iqi或ΔPi、ΔQi效率极为低下。如果用坐标存贮、顺序存贮、链表存贮等方法，尽管可节省存贮单元，但由于对角元与非对角元分开存贮使得数据的读取过程繁琐，不利于其计算及处理。此外，坐标存贮、顺序存贮、链表存贮等方法比较适用于实数元素的存贮，对复数元素的存贮则略为复杂。此外，在对Y、B′、B″阵元素的存贮和读取时未利用3个矩阵中绝大部分非零元素的位置相同的特点。
(3)Y、B′、B″阵数据文件的存贮个数较多、读写时间较长。
如不考虑稀疏性按传统方法分别存贮Y、B′、B″阵的3个数据文件，则在潮流计算程序中要分别打开3个数据文件并读取数据；如考虑稀疏性按链表存贮方法分别对Y、B′、B″阵进行存贮，则存贮文件的个数可达到9个，不利于PQ分解法的实时计算。
(4)Y、B′、B″阵非零元素位置之间的关系未充分利用。
形成Y、B′阵的过程中，在1～(n-1)阶范围内，其静态的非零元素的位置和个数完全一样；而在对Y、B′阵的消元过程中，在1～(n-1)阶范围内，其动态的非零元素的位置和个数也完全一样。在第1～m行，B′、B″阵各行非零元素的位置和个数在静态和动态情况下完全相同；而在第(m+1)～(n-1)行，B″阵的中除对角元为-2外，其余元素均为零，与B′阵第(m+1)～(n-1)行中在静态和动态情况下非零元素的位置和个数不同，且在消元过程中B″阵除对角元由于规格化发生变化外，其余元素均不变化。传统方法中由于未利用该特性而导致存贮单元的浪费和因子表形成过程中一些不必要的计算。
(5)对B′、B″阵形成因子阵的方式不合适。
1)形成因子阵的方法不合适。因Pi-Δfi和Qi-Δei迭代均为常系数线性方程，可用因子表法或LR、LDU、CU三角分解法等求解，一般多用LDU三角分解法。但LDU三角分解法需求解L、D、U三个因子阵，而因子表法计算过程简单，且比LDU三角分解法的计算速度更快。
2)形成因子表时稀疏技术应用不足。B′、B″阵和Y阵元素结构几乎相同，且极度稀疏，若不使用稀疏技术，则计算效率极为低下。三角分解法中稀疏技术应用并不多见，而因子表法中稀疏技术的应用极为广泛，但并不完全到位。
3)对B′、B″阵因子表分开形成。在PQ节点的范围内，B′、B″阵的元素结构完全相同，如果不能根据对B′阵形成因子表的过程快速形成B″阵的因子表，则形成B″阵因子表时要浪费时间。
(6)用因子表进行前代和回代过程中稀疏技术应用不到位。
形成因子表后，有时传统方法会利用因子表上三角的非零元素快速回代求取Δfi、Δei，但各行非零元素的个数和列号分开存贮未能使其计算方式达到最优，而对ΔPi/ei、ΔQi/ei的前代计算则未利用稀疏技术。
(7)有功、无功迭代过程中未引入塞德尔迭代方式。
引入塞德尔迭代方式不但对PQ分解法的计算速度有较大影响，有时还会影响其收敛性，而传统PQ分解法中较少引入塞德尔迭代方式。
由于综上所述原因，传统直角坐标PQ分解法的计算速度远未达到最优，从而几乎不被应用。
发明内容
为了克服上述现有技术的不足，本发明提供了一种对称稀疏因子表法在直角坐标PQ分解法潮流计算中的应用。
本发明读取的数据文件为仅由非零元素构成的三个虚拟数组Y(n,3d1)、B′(n-1,d2)、B″(m,d3)所组成的A(n,d)数据文件，加速数据文件的读取和Ipi、Iqi或ΔPi、ΔQi的计算；利用对称稀疏技术加速因子表的形成；记录B′(n-1,n-1)数组因子表B′(n-2)的形成过程中上三角非零元素的个数和列号，并将其用于加速B″(m,n-1)数组因子表B″(m-1)的形成；将对称稀疏技术继续用于因子表B′(n-2)、B″(m-1)中加速后续的前代及回代计算；在迭代过程中引入塞德尔迭代方式，加速有功、无功迭代。
本发明是通过以下技术方案实现的。
本发明所述的一种对称稀疏因子表法在直角坐标PQ分解法潮流计算中的应用，包括以下步骤：
步骤1：定义数组Y(n,3d1+1)、B′(n-1,n-1)、B″(m,n-1)；
步骤2：从A(n,d)数据文件中将虚拟数组Y(n,3d1)、B′(n-1,d2)及B″(m,d2)的数据分别读入至Y(n,3d1+1)、B′(n-1,n-1)、B″(m,n-1)数组；
(1)在应用PQ分解法潮流程序之前先建立无非零元素的A(n,d)数据文件，将虚拟数组Y(n,3d1)、B′(n-1,d2)、B″(m,d2)对应的数据全部存于A(n,d)数组中，其中d＝3d1+2d2+3，Y(n,3d1)对应的列数为3d1，B′(n-1,d2)、B″(m,d2)对应的列数均为d2，“+3”分别为行号列和S1max、S2max所在的计数列，即表1中的第1～3列；且根据数组结构有d1＝S1max，d2＝S2max，S1max、S2max为数组各行中最大的非零元素节点之和。
(2)传统方式中，Y阵的阶数为n，而B′、B″阵阶数均为(n-1)。在1～(n-1)阶范围内，Y、B′阵的阶数相同时其非零元素的个数和列号完全一样；在第1～m行范围内，B′、B″阵各行非零元素的个数和列号完全相同，而在第(m+1)～(n-1)行，B″阵除对角元为-2外，其余元素均为零，与B′阵不同。因此，只要给定Y(n,3d1)各行中对应的父节点和非零元素子节点之和为S1max，B′(n-1,d2)和B″(m,d2)对应的为S2max。所以Y(n,3d1)中只需存贮主节点及非零元素子节点的列号j和参数，而B′(n-1,d2)、B″(m,d2)中只需存贮所有对应的非零元素参数，而非零元素的列号j可根据S2max由Y(n,3d1)中得到。
(3)利用Si1快速将A(n,d)中Y(n,3d1)的数据读入Y(n,3d1+1)数组，“+1”表示将Si1读入到Y(n,3d1+1)数组的第1列；利用Si2快速将B′(n-1,d2)、B″(m,d2)的数据分别读入B′(n-1,n-1)、B″(m,n-1)数组，但不用将Si2读入到B′(n-1,n-1)、B″(m,n-1)数组，且B″阵第(m+1)～(n-1)行的元素除对角元为-2外，其余元素均为零，因此也无需读入。
(4)A(n,d)数组分为“行号组”、“节点数组”、“Y阵组”、“B′阵组”、“B″阵组”，其元素存贮方式如表1所示。
表1A(n,d)数组元素的存贮方式
注：表1中除含最大非零元素节点对应的行外，并非所有的存贮单元都有数据，而Si1、Si2可保证实际存贮单元数远远小于最大存贮单元数，进一步提高数据的存贮和读写效率。
(5)从A(n,d)数据文件中利用Si1值和Si2值快速将虚拟数组Y(n,3d1)、B′(n-1,d2)、B″(m,d3)对应的数据分别读入Y(n,3d1+1)、B′(n-1,n-1)及B″(m,n-1)数组，利用Si1值和Y(n,3d1+1)可快速计算Ipi、Iqi或ΔPi、ΔQi，用B′(n-1,n-1)、B″(m,n-1)求解Δfi、Δei。
步骤3：根据对称稀疏技术快速形成B′(n-1,n-1)、B″(m,n-1)的因子表B′(n-2)、B″(m-1)，并记录形成B′(n-2)因子表中上三角每行非零元素的个数和列号，而B′(n-2)因子表中第1～m行非零元素的个数及列号可直接用于B″(m-1)的形成过程中；
(1)用对称稀疏技术快速形成B′(n-1,n-1)的因子表B′(n-2)，自动统计记录因子表中每行上三角非零元素u′ij的个数S′i和列号，以便在后续计算中根据S′i和相应的列号直接取用B′(n-2)的非零元素反复对Pi-Δfi方程进行快速前代和回代，计算要点如下。
1)快速确定非零元素和计算元素。根据对角元以右非零元素的值和位置按对称性获得对角元以下非零元素的值和位置，仅计算各非零元素行与列交互点上的元素。
在“逐行规格化、按列消元”过程中，B′(n-2)的每行元素在规格化前，其对角元以右的非零元素与对角元以下的非零元素数值相等、位置对称；在规格化后，对角元以右的非零元素与对角元以下的非零元素数值不等、但位置仍然对称。因此可根据对角元以右的非零元素的值和位置按对称性快速获得对角元以下的非零元素的值和位置，并利用相应的非零元素所在行与所在列的交互点来确定需计算的元素，大大提高对非零元素的判断和计算效率。
2)省去对下三角非零元素的计算。根据1)，对对角元以右的元素进行规格化前，先将其赋给相应位置的对角元以下的素。因此在形成因子表的过程中只需计算对角元及其以右的、非零元素行列交互点上的元素，可减少约50％非零的非对角元素的计算。
3)自动统计并记录上三角非零元素个数和列号。自动统计记录B′(n-2)中每行上三角非零元素u′ij的个数S′i和列角标，并存入相应数组，以便在后续计算中根据S′i和相应的列号直接取用B′(n-2)的非零元素反复快速进行前代和回代，免去该过程对零元素的无效计算。
(2)根据对B′(n-2)中记录的S′i和列号快速形成B″(m,n-1)的因子表B″(m-1)，以便在后续计算中根据S′i和相应的列号直接取用B″(m-1)的非零元素反复对Qi-Δei方程进行快速前代和回代，计算要点如下；
1)B′、B″阵的实际阶数均为(n-1)。但在形成B′、B″阵的过程中，第1～m行B′、B″阵各行静态非零元素b′ij、b″ij的个数Si2和列号完全相同，而在第(m+1)～(n-1)行，B″阵除对角元为-2外，其余元素均为零，与B′阵各行静态非零元素的个数Si2和列号完全不同。
2)在形成B′、B″阵因子表的过程中，第1～m行B′、B″阵因子表各行动态非零元素u′ij、u″ij的个数S′i和列号完全相同；而在第(m+1)～(n-1)行，B″阵因子表除对角元外，其余元素仍全部为零，与B′阵因子表各行动态非零元素的个数S′i和列号完全不同。
3)在形成因子表的消元过程中，对B″阵的前m行的m列，在对对角元以右的元素进行规格化前，可根据对角元以右非零元素的值和位置获得对角元以下非零元素的值和位置，省去对下三角非零元素的计算；对B″阵的前m行的(m+1)～(n-1)列，仅需计算其非零元素交互点上的元素；对B″阵(m+1)～(n-1)行的所有元素，由于只有对角元-2，其它元素均为零，且在消元过程中，除对角元-2在规格化时变化外，所有零元素均不变化，因此无需在B″阵的因子表中考虑(m+1)～(n-1)行非对角元的动态元素，即对B″阵的B″(m,n-1)数组只要形成含PQ节点部分的因子表B″(m-1)。
4)由于在形成B″阵因子表B″(m-1)的过程中，其前m行非零元素u″ij的个数和列号与B′阵因子表B′(n-2)中非零元素u′ij的个数S′i和列号完全一样，因此无需对其再进行非零判断和自动统计记录，而在后续计算中可根据对应B″(m-1)的S′i和相应的列号直接取用B″(m-1)的非零元素反复对Qi-Δei方程进行快速前代和回代。
步骤4：根据对称稀疏技术利用因子表B′(n-2)中下三角非零元素l′ji和S′i快速前代计算ΔPi/ei，利用其上三角非零元素u′ij和S′i快速回代求取Δfi。利用因子表B″(m-1)中下三角非零元素l″ji和S′i快速前代计算ΔQi/ei，利用其上三角非零元素u″ij和S′i快速回代求取Δei，并在有功、无功迭代过程中引入塞德尔方式；
(1)通过因子表形成过程中记录的B′(n-2)中上三角非零元素u′ij的位置，按对称性得到下三角非零元素l′ji的位置，用l′ji元素和S′i快速前代计算ΔPi/ei；用u′ij元素和S′i快速回代求取Δfi；通过B″(m-1)中上三角非零元素u″ij的位置，按对称性得到下三角非零元素l″ji的位置，用l″ji元素和S′i快速前代计算ΔQi/ei；用u″ij元素和S′i快速回代求取Δei，大大加快前代、回代过程。
(2)引入塞德尔迭代方式后，PQ分解法的迭代过程为： 即先计算有功，然后迭代有功求出fi的修正量和新值，再用fi的新值计算无功，继而迭代无功求出ei的修正量和新值，再用ei的新值计算有功，依次循环加快收敛。
步骤5：判断是否满足收敛条件；
如果不满足收敛条件，则利用本次迭代得到的Δfi、Δei继续进行下一次的前代和回代计算；如果满足收敛条件，则执行步骤6。
步骤6：计算平衡节点的功率及支路功率并输出计算结果。
本发明方法与传统的PQ分解法相比，具有以下优点。
(1)仅读取A(n,d)数组1个数据文件，而A(n,d)数组中仅存贮相应非零元素的数据，存贮单元数及数据文件的读取时间大大减少。
(2)可直接根据Y(n,3d1+1)数组高效计算Ipi、Iqi或ΔPi、ΔQi。
(3)形成B″阵和B″阵的因子表时考虑其PV节点在形成和消元过程的特点，对(n-1)阶的B″阵仅需形成B″(m,n-1)数组，而对B″(m,n-1)数组只要形成含PQ节点部分的B″(m-1)因子表。
(4)根据对称稀疏性快速获得B′、B″阵的因子表，并仅记录B′阵因子表中上三角非零元素的个数S′i及列号，以便后续计算中直接取用B′(n-2)、B″(m-1)阵的非零元素反复快速前代和回代。
(5)利用下三角非零元素l′ji、l″ji和S′i可快速前代计算ΔPi/ei、ΔQi/ei。
(6)利用上三角非零元素u′ij、u″ij和S′i可快速回代求取Δfi、Δei。
(7)有功、无功迭代过程中引入塞德尔迭代方式，加快潮流收敛。
附图说明
图1为传统直角坐标PQ分解法潮流的计算流程图。
图2为本发明方法用于直角坐标PQ分解法潮流的计算流程图。
具体实施方式
本发明将通过以下实施例作进一步说明。
实施例1。分别用传统的直角坐标PQ分解法和本发明法用于直角坐标PQ分解法对IEEE-30、-57、-118节点系统进行潮流计算，比较其读取数据文件、形成因子表、有功无功迭代和潮流计算(总)的平均计算时间。计算结果如表2所示。
表2传统方法和本发明对IEEE系统进行PQ分解潮流计算过程时间的比较
tr.c、tf.c、ti.c、tp.c：传统PQ分解法中不考虑稀疏性时，读取数据文件、形成因子表、有功无功迭代的平均计算时间及潮流计算(总)时间；
tf.s.c、tp.s.c：传统方法PQ分解法中仅判断对角元以下的非零元素形成因子表的平均计算时间及潮流计算(总)时间(除形成因子表外，其他同传统方法中不考虑元素稀疏性)；
tr.new、tf.new、ti.new、tp.new：本发明中读取数据文件时间、形成因子表、有功无功迭代的平均计算时间及PQ分解法潮流计算(总)时间；
tr.new/tr.c、tf.new/tf.c、ti.new/ti.c、tp.new/tp.c：本发明与传统方法不考虑稀疏性时读取数据文件时间、形成因子表时间、有功无功迭代时间、潮流计算(总)时间的百分比。
ttfs.c/tx.c、tp.s.c/tp.c：传统方法中仅判断对角元以右的非零元素形成因子表的时间和PQ分解法潮流计算(总)时间与传统方法中不考虑任何稀疏性时相应计算结果的百分比。
以IEEE-118节点系统为例，计算结果比较分析如下：
(1)读取数据文件数及时间的比较
传统因子表法分别将Y、B′、B″阵元素分别存放在Y(n,2n)、B′(n-1,n-1)和B″(n-1,n-1)数组中，需读取3个数据文件，而本发明方法存放在A(n,d)数组中，仅需读取1个数据文件。本发明读取数据文件的时间仅为传统方法的6.12％。
(2)形成因子表时间的比较
本发明形成因子表的时间仅为传统方法的6.15％；如果在传统方法中判断对角元以下的非零元素形成因子表，则其为传统方法中不考虑稀疏性时的7.84％。
(3)本发明方法有功、无功迭代时间仅为传统方法的11.86％。
(4)潮流计算(总)时间的比较
本发明潮流计算(总)的时间仅为传统方法的6.16％。传统方法中判断对角元以下的非零元素形成因子表与传统方法中不考虑稀疏性时形成因子表的速度相比虽然大大提高，但前者潮流计算的时间与后者相比几乎没有变化。
上述计算结果表明，本发明提出的方法与传统的直角坐标PQ分解法相比，无论在读取数据文件、形成因子表、有功无功迭代还是潮流计算方面，计算速度均大大加快，说明本发明方法的计算效率极高。且电力系统节点数越多，本发明方法优势越大。
本发明可以采用任何一种编程语言和编程环境实现，这里采用VisualC++语言，开发环境是VisualC++，运行平台是Intel(R)Corei7-4790CPU3.60GHZ，内存8.00GB。
Claims (1)
1.一种对称稀疏因子表法在直角坐标PQ分解法潮流计算中的应用，其特征包括以下步骤：
步骤1：定义数组Y(n,3d1+1)、B′(n-1,n-1)、B″(m,n-1)；
步骤2：从A(n,d)数据文件中将虚拟数组Y(n,3d1)、B′(n-1,d2)及B″(m,d2)的数据分别读入至Y(n,3d1+1)、B′(n-1,n-1)、B″(m,n-1)数组，其中Y(n,3d1+1)用于快速计算Ipi、Iqi或ΔPi、ΔQi，B′(n-1,n-1)、B″(m,n-1)数组用于求解Δfi、Δei；
步骤3：根据对称稀疏技术快速形成B′(n-1,n-1)、B″(m,n-1)的因子表B′(n-2)、B″(m-1)，并记录形成B′(n-2)因子表中每行上三角非零元素的个数和列号，而B′(n-2)因子表中第1～m行非零元素的个数及列号可直接用于B″(m-1)的形成过程中；
(1)用对称稀疏技术快速形成B′(n-1,n-1)的因子表B′(n-2)，自动统计记录因子表中上三角每行非零元素u′ij的个数S′i和列号，在后续计算中根据S′i和相应的列号直接取用B′(n-2)的非零元素反复对Pi-Δfi方程进行快速前代和回代；
(2)根据对B′(n-2)中记录的S′i和列号用对称稀疏技术快速形成B″(m,n-1)的因子表B″(m-1)，在后续计算中根据S′i和相应的列号直接取用B″(m-1)的非零元素反复对Qi-Δei方程进行快速前代和回代；其中，在形成B′、B″阵和形成B′、B″阵因子表的过程中，第1～m行B′、B″阵静态和动态非零元素的个数Si2和S′i以及列号完全相同，而在第(m+1)～(n-1)行，B″阵除对角元外，其余元素均为零，与B′阵静态和动态非零元素的个数Si2和S′i以及列号完全不同；
步骤4：根据对称稀疏技术利用因子表B′(n-2)中下三角非零元素l′ji和S′i快速前代计算ΔPi/ei，利用其上三角非零元素u′ij和S′i快速回代求取Δfi；利用因子表B″(m-1)中下三角非零元素l″ji和S′i快速前代计算ΔQi/ei，利用其上三角非零元素u″ij和S′i快速回代求取Δei，并在有功、无功迭代过程中引入塞德尔方式；
步骤5：判断是否满足收敛条件；
如果不满足收敛条件，则利用本次迭代得到的Δfi、Δei继续进行下一次的前代和回代计算；如果满足收敛条件，则执行步骤6；
步骤6：计算平衡节点的功率及支路功率并输出计算结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610130735.5A CN105703359B (zh) | 2016-03-08 | 2016-03-08 | 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610130735.5A CN105703359B (zh) | 2016-03-08 | 2016-03-08 | 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105703359A true CN105703359A (zh) | 2016-06-22 |
| CN105703359B CN105703359B (zh) | 2021-01-01 |
Family
ID=56221081
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610130735.5A Active CN105703359B (zh) | 2016-03-08 | 2016-03-08 | 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105703359B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109446551A (zh) * | 2018-09-13 | 2019-03-08 | 南昌大学 | 电力系统pq法潮流中随机非零元素因子表的形成及应用 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292520A1 (en) * | 2006-07-27 | 2009-11-26 | Drexel University | Solver for hardware based computing |
| CN102013680A (zh) * | 2010-12-13 | 2011-04-13 | 大连海事大学 | 一种电力系统快速分解法潮流计算方法 |
| CN103065059A (zh) * | 2013-01-29 | 2013-04-24 | 河海大学 | 一种基于变量代换的辐射型配电网潮流计算方法 |
| US20130304266A1 (en) * | 2012-04-13 | 2013-11-14 | Regents Of The University Of Minnesota | State estimation of electrical power networks using semidefinite relaxation |
| CN104158182A (zh) * | 2014-08-18 | 2014-11-19 | 国家电网公司 | 一种大规模电网潮流修正方程并行求解方法 |
| CN104218577A (zh) * | 2014-09-19 | 2014-12-17 | 国家电网公司 | 一种基于节点电压的主动配电网三相潮流计算方法 |
| CN104317553A (zh) * | 2014-10-13 | 2015-01-28 | 南昌大学 | 一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法 |
| CN104933528A (zh) * | 2015-06-24 | 2015-09-23 | 南昌大学 | 一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法 |
| CN105045767A (zh) * | 2015-06-24 | 2015-11-11 | 南昌大学 | 一种快速存贮及读取电力系统稀疏矩阵数据的方法 |
| CN105354422A (zh) * | 2015-11-12 | 2016-02-24 | 南昌大学 | 一种基于对称稀疏矩阵技术快速求取极坐标牛顿-拉夫逊法潮流的方法 |
| CN105375468A (zh) * | 2015-11-12 | 2016-03-02 | 南昌大学 | 一种基于对称稀疏矩阵技术快速求取直角坐标牛顿-拉夫逊法潮流的方法 |
| CN105786769A (zh) * | 2016-02-15 | 2016-07-20 | 南昌大学 | 一种基于快速数据读取及对称稀疏因子表法在极坐标pq分解法潮流中的应用 |
-
2016
- 2016-03-08 CN CN201610130735.5A patent/CN105703359B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090292520A1 (en) * | 2006-07-27 | 2009-11-26 | Drexel University | Solver for hardware based computing |
| CN102013680A (zh) * | 2010-12-13 | 2011-04-13 | 大连海事大学 | 一种电力系统快速分解法潮流计算方法 |
| US20130304266A1 (en) * | 2012-04-13 | 2013-11-14 | Regents Of The University Of Minnesota | State estimation of electrical power networks using semidefinite relaxation |
| CN103065059A (zh) * | 2013-01-29 | 2013-04-24 | 河海大学 | 一种基于变量代换的辐射型配电网潮流计算方法 |
| CN104158182A (zh) * | 2014-08-18 | 2014-11-19 | 国家电网公司 | 一种大规模电网潮流修正方程并行求解方法 |
| CN104218577A (zh) * | 2014-09-19 | 2014-12-17 | 国家电网公司 | 一种基于节点电压的主动配电网三相潮流计算方法 |
| CN104317553A (zh) * | 2014-10-13 | 2015-01-28 | 南昌大学 | 一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法 |
| CN104933528A (zh) * | 2015-06-24 | 2015-09-23 | 南昌大学 | 一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法 |
| CN105045767A (zh) * | 2015-06-24 | 2015-11-11 | 南昌大学 | 一种快速存贮及读取电力系统稀疏矩阵数据的方法 |
| CN105354422A (zh) * | 2015-11-12 | 2016-02-24 | 南昌大学 | 一种基于对称稀疏矩阵技术快速求取极坐标牛顿-拉夫逊法潮流的方法 |
| CN105375468A (zh) * | 2015-11-12 | 2016-03-02 | 南昌大学 | 一种基于对称稀疏矩阵技术快速求取直角坐标牛顿-拉夫逊法潮流的方法 |
| CN105786769A (zh) * | 2016-02-15 | 2016-07-20 | 南昌大学 | 一种基于快速数据读取及对称稀疏因子表法在极坐标pq分解法潮流中的应用 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109446551A (zh) * | 2018-09-13 | 2019-03-08 | 南昌大学 | 电力系统pq法潮流中随机非零元素因子表的形成及应用 |
| CN109446551B (zh) * | 2018-09-13 | 2023-03-14 | 南昌大学 | 电力系统pq法潮流中随机非零元素因子表的形成及应用 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105703359B (zh) | 2021-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104317553B (zh) | 一种基于稀疏矩阵技术快速形成及读写电力系统节点导纳矩阵数据的方法 | |
| CN104933528B (zh) | 一种基于稀疏矩阵技术快速形成电力系统潮流计算中雅可比矩阵的方法 | |
| WO2018027706A1 (zh) | Fft处理器及运算方法 | |
| US20220129523A1 (en) | Method, circuit, and soc for performing matrix multiplication operation | |
| CN104714928B (zh) | 一种基于对称稀疏矩阵技术的高斯消元法求取电力系统节点阻抗矩阵的方法 | |
| CN115270633A (zh) | 燃料电池三维物理场的预测方法、系统、设备及介质 | |
| CN111327048A (zh) | 一种含三端snop的配电网鲁棒运行优化方法 | |
| CN105703359A (zh) | 一种对称稀疏因子表法在直角坐标pq分解法潮流计算中的应用 | |
| CN105354422A (zh) | 一种基于对称稀疏矩阵技术快速求取极坐标牛顿-拉夫逊法潮流的方法 | |
| CN105786769A (zh) | 一种基于快速数据读取及对称稀疏因子表法在极坐标pq分解法潮流中的应用 | |
| CN105591388B (zh) | 一种基于稀疏技术的电力系统直角坐标pq分解法潮流数据存贮方法 | |
| CN105045767B (zh) | 一种快速存贮及读取电力系统稀疏矩阵数据的方法 | |
| CN115222870B (zh) | 非结构网格隐式lusgs均衡着色方法、设备及介质 | |
| CN115718986B (zh) | 一种基于分布式内存架构的多核并行时域仿真方法 | |
| CN116231657A (zh) | 一种输配电网全局碳流分布式确定方法及装置 | |
| CN116227276A (zh) | 多重材料分布情况下电场电位多尺度有限元计算方法、设备及存储介质 | |
| CN113379327B (zh) | 一种电力系统多区域并行状态估计方法及系统 | |
| CN109447839B (zh) | 非零元素的随机生成及形成电力系统雅可比矩阵的方法 | |
| Gao et al. | Revisiting thread configuration of SpMV kernels on GPU: A machine learning based approach | |
| CN111614078A (zh) | 电力系统小干扰稳定性分析方法、装置、设备及存储介质 | |
| CN110968953B (zh) | 一种基于嵌套对角加边形式的电力系统暂态稳定仿真并行计算方法 | |
| CN105375468B (zh) | 一种基于对称稀疏矩阵技术快速求取直角坐标牛顿‑拉夫逊法潮流的方法 | |
| CN110991034B (zh) | 基于全并行嵌套bbdf的电力系统暂态稳定仿真并行计算方法 | |
| CN113312312B (zh) | 一种基于lsm高效查询流数据的分布式索引方法及系统 | |
| CN109886398A (zh) | 神经网络矩阵乘法运算方法及相关产品 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20160622 Assignee: Jiangxi Jimu Technology Co.,Ltd. Assignor: Nanchang University Contract record no.: X2023980053717 Denomination of invention: Application of a Symmetric Sparse Factor Table Method in Power Flow Calculation of Cartesian Coordinate PQ Decomposition Method Granted publication date: 20210101 License type: Common License Record date: 20231222 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| OL01 | Intention to license declared | ||
| OL01 | Intention to license declared |