CN108233382A - 一种提取直角坐标潮流方程雅可比矩阵的方法 - Google Patents
一种提取直角坐标潮流方程雅可比矩阵的方法 Download PDFInfo
- Publication number
- CN108233382A CN108233382A CN201711229577.XA CN201711229577A CN108233382A CN 108233382 A CN108233382 A CN 108233382A CN 201711229577 A CN201711229577 A CN 201711229577A CN 108233382 A CN108233382 A CN 108233382A
- Authority
- CN
- China
- Prior art keywords
- formula
- matrix
- represent
- power
- node
- 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
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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种提取直角坐标潮流方程雅可比矩阵的方法,包括如下步骤:从节点复功率方程出发,根据复数的节点功率方程,按节点电压的直角坐标形式分别对e和f求偏导;然后展开得到复功率不平衡量表达式;分别提取替换后矩阵的实部和虚部,并将其带入复功率不平衡量表达式;最后整理为矩阵形式,其系数矩阵即为雅可比矩阵。本发明的方法极大地降低了雅可比矩阵形成过程的复杂度,提升潮流算法的计算速度,减少算法对计算机资源占用。
Description
技术领域
本发明属于电力潮流计算方法技术领域,具体涉及一种提取直角坐标潮流方程雅可比矩阵的方法。
背景技术
电力系统潮流计算是分析电力系统稳态运行情况的一种基本计算。数学原理上,它是求解多元非线性方程组的问题,应用最为广泛的求解方法是牛顿法。使用牛顿法进行求解的计算量,很大程度在于形成雅可比矩阵这一步骤,降低这一步骤的复杂度,是提高算法的计算效率,减少算法对计算机资源占用的重要手段。
发明内容
本发明的目的在于针对现有技术中的上述问题,提供一种提取直角坐标潮流方程雅可比矩阵的方法,可降低雅可比矩阵形成过程的复杂度,提升潮流算法的计算速度。
为实现上述发明目的,本发明采用了如下技术方案:
一种提取直角坐标潮流方程雅可比矩阵的方法,包括如下步骤:
S1.根据复数的节点功率方程按节点电压的直角坐标形式分别对e和f求偏导,得:
其中,表示节点功率;表示节点功率增量;diag表示取对应矢量的对角矩阵; 表示电压向量;表示的共轭复数; 表示节点导纳矩阵;表示的共轭复数;e表示电压向量的实部;f表示电压向量的虚部;
S2.展开步骤S1中得到的式(1),得到复功率不平衡量表达式:
式中,Δe表示电压实部增量,Δf表示电压虚部增量;
S3.对步骤S2中得到的式(2)进行符号化替换,令矩阵
将式(2)变为:
S4.从步骤S3得到的式(3)中分别提取矩阵D、矩阵J的实部和虚部,得:
式中,Re表示提取对应矩阵的实部;Im表示提取对应矩阵的虚部;DR表示D的实部;JR表示J的实部;DI表示D的虚部;JI表示J的虚部;
S5.将步骤S4中的式(4)代入步骤S3中的式(3)中,得:
S6.将步骤S5得到的式(5)整理为矩阵形式,其系数矩阵即为雅可比矩阵:
式中,ΔP表示有功功率增量;ΔQ表示无功功率增量。
相比于现有技术,本发明的优势在于:
本发明所提供的一种提取直角坐标潮流方程雅可比矩阵的方法,在现有使用牛顿法求解潮流问题形成雅可比矩阵的过程中,从节点复功率方程出发,直接提取雅可比矩阵。本发明的方法极大地降低了雅可比矩阵形成过程的复杂度,提升大电网潮流计算速度,减少算法对计算机CPU、内存的占用。
附图说明
图1是本发明一种提取直角坐标潮流方程雅可比矩阵的方法的流程图。
具体实施方式
以下结合实施例及其附图对本发明技术方案作进一步非限制性的详细说明。
如图1所示,一种提取直角坐标潮流方程雅可比矩阵的方法,包括如下步骤:
S1.从节点复功率方程出发,根据复数的节点功率方程按节点电压的直角坐标形式分别对e和f求偏导,得:
其中,表示节点复功率;表示节点复功率增量;diag表示取对应矢量的对角矩阵; 表示电压向量;表示的共轭复数; 表示节点导纳矩阵;表示的共轭复数;e表示电压向量的实部;f表示电压向量的虚部;为实数向量;
S2.展开步骤S1中得到的式(1),得到复功率不平衡量表达式:
式中,Δe表示电压实部增量,Δf表示电压虚部增量;
S3.对步骤S2中得到的式(2)进行符号化替换,令矩阵
将式(2)变为:
S4.从步骤S3得到的式(3)中分别提取矩阵D、矩阵J的实部和虚部,得:
式中,Re表示提取对应矩阵的实部;Im表示提取对应矩阵的虚部;DR表示D的实部;JR表示J的实部;DI表示D的虚部;JI表示J的虚部;j表示虚部符号;
S5.将步骤S4中的式(4)代入步骤S3中的式(3)中,得:
S6.将步骤S5得到的式(5)整理为矩阵形式,其系数矩阵即为雅可比矩阵:
式中,ΔP表示有功功率增量;ΔQ表示无功功率增量。
实施例:
对IEEE4潮流分析数据进行电力系统潮流计算,包含如下步骤:
S1:给出电压向量初始值,并给出收敛精度,为:ε=10-6;
S2:形成节点导纳矩阵,得:
S3:根据复数的节点功率方程按节点电压的直角坐标形式分别对e和f求偏导,得:
其中,表示节点复功率;表示节点复功率增量;diag表示取对应矢量的对角矩阵; 表示电压向量;表示的共轭复数; 表示节点导纳矩阵;表示的共轭复数;e表示电压向量的实部;f表示电压向量的虚部;为实数向量;
S4:展开式(1),得到复功率不平衡量表达式:
式中,Δe表示电压实部增量,Δf表示电压虚部增量;
S5:对式(2)进行符号化替换,令矩阵将式(2)变为:
S6:从式(3)中分别提取矩阵D、矩阵J的实部和虚部,得:
式中,Re表示提取对应矩阵的实部;Im表示提取对应矩阵的虚部;DR表示D的实部;JR表示J的实部;DI表示D的虚部;JI表示J的虚部;
S7:将式(4)代入式(3)中,得:
S8:将式(5)整理为矩阵形式,其系数矩阵即为雅可比矩阵:
式中,P表示有功功率;ΔP表示有功功率增量;Q表示无功功率;ΔQ表示无功功率增量;
按上述步骤S3-S8,并处理PV节点电压恒定问题,形成雅可比矩阵,在第一次迭代过程中,雅可比矩阵为:
S9:进行迭代,达到收敛精度后输出计算结果,为:V1=0.9846-0.0086i,V2=0.92587-0.1084i,V3=1.0924+0.1290i,V1=1.0500-0.0000i。
运用本发明的方法进行上述电力系统潮流计算,其中的雅可比矩阵的形成过程得到了极大的简化,在实际应用中,能够极大提高大电网潮流计算效率,降低对电脑CPU、内存的占用。
综上所述,本发明在现有使用牛顿法求解潮流问题形成雅可比矩阵的过程中,从节点复功率方程出发,根据复数的节点功率方程,按节点电压的直角坐标形式分别对e和f求偏导;然后展开得到复功率不平衡量表达式;分别提取替换后矩阵的实部和虚部,并将其带入复功率不平衡量表达式;最后整理为矩阵形式,其系数矩阵即为雅可比矩阵。
以上所揭露的仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或变型,都应涵盖在本发明的保护范围之内。
Claims (1)
1.一种提取直角坐标潮流方程雅可比矩阵的方法,其特征在于,包括如下步骤:
S1.根据复数的节点功率方程按节点电压的直角坐标形式分别对e和f求偏导,得:
其中,表示节点复功率;表示节点复功率增量;diag表示取对应矢量的对角矩阵; 表示电压向量;表示的共轭复数; 表示节点导纳矩阵;表示的共轭复数;,e表示电压向量的实部;,f表示电压向量的虚部;
S2.展开步骤S1中得到的式(1),得到复功率不平衡量表达式:
式中,Δe表示电压实部增量,Δf表示电压虚部增量;
S3.对步骤S2中得到的式(2)进行符号化替换,令矩阵将式(2)变为:
S4.从步骤S3得到的式(3)中分别提取矩阵D、矩阵J的实部和虚部,得:
式中,Re表示提取对应矩阵的实部;Im表示提取对应矩阵的虚部;DR表示D的实部;JR表示J的实部;DI表示D的虚部;JI表示J的虚部;
S5.将步骤S4中的式(4)代入步骤S3中的式(3)中,得:
S6.将步骤S5得到的式(5)整理为矩阵形式,其系数矩阵即为雅可比矩阵:
式中,ΔP表示有功功率增量;ΔQ表示无功功率增量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711229577.XA CN108233382B (zh) | 2017-11-29 | 2017-11-29 | 一种提取直角坐标潮流方程雅可比矩阵的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711229577.XA CN108233382B (zh) | 2017-11-29 | 2017-11-29 | 一种提取直角坐标潮流方程雅可比矩阵的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108233382A true CN108233382A (zh) | 2018-06-29 |
CN108233382B CN108233382B (zh) | 2020-09-01 |
Family
ID=62653588
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711229577.XA Active CN108233382B (zh) | 2017-11-29 | 2017-11-29 | 一种提取直角坐标潮流方程雅可比矩阵的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108233382B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112994021A (zh) * | 2021-04-27 | 2021-06-18 | 广西大学 | 一种提取直角坐标节点复功率方程海森矩阵的计算方法 |
CN113158126A (zh) * | 2021-04-27 | 2021-07-23 | 广西大学 | 一种提取极坐标节点复功率方程海森矩阵的计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199525A (zh) * | 2013-03-27 | 2013-07-10 | 河海大学 | 一种基于等效节点注入电流变换的配电网潮流计算方法 |
US20150112498A1 (en) * | 2013-09-23 | 2015-04-23 | Sureshchandra B Patel | Methods of Patel Loadflow Computation for Electrical Power System |
CN106356859A (zh) * | 2016-09-29 | 2017-01-25 | 大连海事大学 | 一种基于Matlab的直角坐标牛顿法潮流计算方法 |
CN107181260A (zh) * | 2017-07-10 | 2017-09-19 | 大连海事大学 | 基于Matlab稀疏矩阵直角坐标牛顿法潮流计算方法 |
-
2017
- 2017-11-29 CN CN201711229577.XA patent/CN108233382B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199525A (zh) * | 2013-03-27 | 2013-07-10 | 河海大学 | 一种基于等效节点注入电流变换的配电网潮流计算方法 |
US20150112498A1 (en) * | 2013-09-23 | 2015-04-23 | Sureshchandra B Patel | Methods of Patel Loadflow Computation for Electrical Power System |
CN106356859A (zh) * | 2016-09-29 | 2017-01-25 | 大连海事大学 | 一种基于Matlab的直角坐标牛顿法潮流计算方法 |
CN107181260A (zh) * | 2017-07-10 | 2017-09-19 | 大连海事大学 | 基于Matlab稀疏矩阵直角坐标牛顿法潮流计算方法 |
Non-Patent Citations (4)
Title |
---|
YANG YUDE等: "Vectorial power flow calculation based on nonlinear programming", 《 2008 IEEE POWER AND ENERGY SOCIETY GENERAL MEETING - CONVERSION AND DELIVERY OF ELECTRICAL ENERGY IN THE 21ST CENTURY》 * |
戴承承: "基于直角坐标的电力系统潮流计算方法研究", 《电气开关》 * |
李传栋: "向量化编程技术及其在大型电网潮流计算中的应用", 《福建农林大学学报(自然科学版)》 * |
李婷婷: "小阻抗直角坐标牛顿潮流算法发散机理研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112994021A (zh) * | 2021-04-27 | 2021-06-18 | 广西大学 | 一种提取直角坐标节点复功率方程海森矩阵的计算方法 |
CN113158126A (zh) * | 2021-04-27 | 2021-07-23 | 广西大学 | 一种提取极坐标节点复功率方程海森矩阵的计算方法 |
CN112994021B (zh) * | 2021-04-27 | 2022-11-01 | 广西大学 | 一种提取直角坐标节点复功率方程海森矩阵的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108233382B (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107947192B (zh) | 一种下垂控制型孤岛微电网的无功优化配置方法 | |
CN110705762B (zh) | 基于矩阵填充的泛在电力物联网感知数据缺失修复方法 | |
Xu et al. | Trajectory sensitivity analysis on the equivalent one‐machine‐infinite‐bus of multi‐machine systems for preventive transient stability control | |
CN104156542B (zh) | 一种基于隐式投影的有源配电系统稳定性仿真方法 | |
Ke et al. | Intelligent islanding detection method for photovoltaic power system based on Adaboost algorithm | |
CN111797510A (zh) | 一种新能源场站短路比的计算方法及系统 | |
Dezelak et al. | Comparison between the particle swarm optimisation and differential evolution approaches for the optimal proportional–integral controllers design during photovoltaic power plants modelling | |
CN108233382A (zh) | 一种提取直角坐标潮流方程雅可比矩阵的方法 | |
CN103915839B (zh) | 一种用于分析含风电电力系统随机稳定性的方法 | |
CN104022507A (zh) | 一种直角坐标牛顿法潮流计算方法 | |
Abdel‐Akher et al. | Initialised load‐flow analysis based on Lagrange polynomial approximation for efficient quasi‐static time‐series simulation | |
Chang et al. | Analysis on current characteristics of PMSG under grid three‐phase fault | |
CN106991229B (zh) | 一种针对复杂拓扑的风电场等值建模方法 | |
CN105203883A (zh) | 基于双密度小波网络的无偏估计光伏逆变器参数辨识方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN110460104B (zh) | 一种光伏逆变器等效正、负序故障阻抗推导方法 | |
Yang et al. | Parallel solution of transient stability constrained optimal power flow by exact optimality condition decomposition | |
CN107957974A (zh) | 一种提取极坐标潮流方程雅可比矩阵的方法 | |
CN107425519B (zh) | 含分布式电源的三相配电网最大供电能力计算方法 | |
CN104182566B (zh) | 用于电网系统中发生部分变化后的求解元件节点电压方法 | |
Ghosh et al. | Optimal sizing and placement of DG units in radial distribution system using cuckoo search algorithm | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN111817359B (zh) | 一种基于等效电路的微电网可求解边界分析的方法 | |
CN110417020B (zh) | 一种处理非光滑约束的综合能源系统潮流计算方法及系统 | |
CN104636630B (zh) | 一种基于均值滤波和等梯度的热电厂蒸汽压力滤波方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |