CN112667962A - 一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法 - Google Patents

一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法 Download PDF

Info

Publication number
CN112667962A
CN112667962A CN202011257225.7A CN202011257225A CN112667962A CN 112667962 A CN112667962 A CN 112667962A CN 202011257225 A CN202011257225 A CN 202011257225A CN 112667962 A CN112667962 A CN 112667962A
Authority
CN
China
Prior art keywords
harmonic
matrix
admittance
transformer
power
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
Application number
CN202011257225.7A
Other languages
English (en)
Other versions
CN112667962B (zh
Inventor
张逸
林金榕
李传栋
邵振国
刘雄飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN202011257225.7A priority Critical patent/CN112667962B/zh
Publication of CN112667962A publication Critical patent/CN112667962A/zh
Application granted granted Critical
Publication of CN112667962B publication Critical patent/CN112667962B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，包括步骤S1：设置谐波次数，各种谐波源统计模型初始值，迭代误差限；步骤S2：由输入的网络参数生成基波和谐波导纳矩阵；步骤S3：计算基波潮流；步骤S4：检测存在谐波源的节点，并生成索引矩阵，利用基波潮流结果，由谐波源统计模型预测谐波数据，进行谐波潮流计算；步骤S5：利用谐波网损功率更新基波负荷功率；步骤S6：判断前后节点注入功率误差是否小于迭代误差限，若是则计算收敛，输出谐波潮流计算结果，否则返回步骤S3。本发明基于支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法以及基于相量矩阵的思维进行程序架构，实现了大规模电网谐波潮流仿真的快速开发和高效计算。

Description

一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算 方法
技术领域
本发明涉及电力系统安全领域，特别是一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法。
背景技术
电力系统的谐波问题在20世纪30年代就引起了人们的关注，随后由于现代工业技术的发展，电网中的非线性负荷大量增加，如工业中广泛应用的电弧和接触焊设备、矿热炉、硅铁炉、高频炉等，同时，电力电子技术的迅速发展，使晶闸管整流和换流技术得到广泛应用，这都使得电网中的谐波大量增加。谐波对电力系统的影响和危害是十分严重的，其主要表现为降低供电系统的功率因数、使电器设备过热，绝缘老化，产生振动和噪声，缩短用电设备的寿命、引起电力系统的并联和串联谐振损坏电容器等。为了解决上述谐波危害问题，70年代以来，各工业国家纷纷投入力量对谐波问题进行研究，作为谐波问题研究的重要内容，谐波潮流计算也得到了相应发展。
谐波潮流计算是研究谐波问题的重要内容，通过谐波潮流计算可以描绘电网络的谐波潮流分布，得出电网络各节点的谐波指标，这是评估电力系统安全运行的重要依据，根据谐波潮流计算结果还可分析谐波产生的原因，并进一步研究谐波治理措施。
当前的谐波潮流计算不适合运用于大规模网络，对谐波源的拓展不方便，因此程序的效率不高，不适合于工程实际情况。
发明内容
有鉴于此，本发明的目的是提供一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，基于相量矩阵的思维进行程序架构，实现了大规模电网谐波潮流仿真的快速开发和高效计算。
本发明采用以下方案实现：一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，包括以下步骤：
步骤S1：设置谐波次数，设置谐波源统计模型初始值，算法迭代误差限；
步骤S2：由包括变压器参数、线路参数、节点电压、发电机参数和负荷参数生成基波和谐波导纳矩阵；
步骤S3：计算基波潮流；
步骤S4：检测存在谐波源的节点，利用基波潮流结果和谐波源统计模型预测谐波数据得到谐波电流，然后再进行谐波潮流计算，得到谐波损耗；
步骤S5：利用谐波网损功率更新基波负荷功率；即利用步骤S4 得到的谐波损耗更新谐波源对应节点的功率，重复执行步骤S3，重新进行基波潮流计算；
步骤S6：判断前后节点注入功率误差是否小于迭代误差限，若是则计算收敛，输出谐波潮流计算结果，否则返回步骤S3。
进一步地，所述步骤S1的具体内容为：
设置要计算的谐波次数h为谐波次数，2≤h≤25，输入谐波源的统计模型的谐波电流参数，以及基波、谐波算法的迭代收敛误差值，小于10-3。
进一步地，步骤S2中所述建立基波导纳矩阵的具体内容为：
步骤S2中所述建立基波导纳矩阵的具体内容为：
基波线路：
Figure BDA0002773376590000031
Figure BDA0002773376590000032
表示第i条线路的电阻、
Figure BDA0002773376590000033
表示第i条线路的电抗,
Figure BDA0002773376590000034
表示第i 条线路t端电纳，
Figure BDA0002773376590000035
表示第i条线路f端电纳；
Figure BDA0002773376590000036
表示第i条线路t端的自导纳,
Figure BDA0002773376590000037
表示第i条线路f端的自导纳，
Figure BDA0002773376590000038
表示第i条线路f端对t端的互导纳，
Figure BDA0002773376590000039
表示第i条线路t端对f端的互导纳；i表示第i条线路，f表示首端，t表示末端；
关联矩阵：Cfl表示f端索引矩阵，Ctl表示t端索引矩阵，为依据节点编号生成的矩阵；
Figure BDA00027733765900000310
m表示总线路数,[*]表示矩阵；
基波变压器：
Figure BDA00027733765900000311
Figure BDA00027733765900000312
Figure BDA0002773376590000041
表示第q台变压器的电阻，
Figure BDA0002773376590000042
表示第q台变压器的电抗，tapq表示第q台变压器变比，
Figure BDA0002773376590000043
表示第q台变压器t端自导纳，
Figure BDA0002773376590000044
表示第q台变压器f端自导纳；
Figure BDA0002773376590000045
表示第q台变压器f端对t端互导纳，
Figure BDA0002773376590000046
表示第q台变压器t端对f端互导纳，f表示首端，t 表示末端，q表示变压器编号；
关联矩阵：CTf表示f端索引矩阵，CTt表示t端索引矩阵，根据变压器节点编号顺序生成的矩阵；
Figure BDA0002773376590000047
n表示总变压器台数，[*]表示矩阵；
基波总导纳：
Figure BDA0002773376590000048
进一步地，步骤S2中所述建立谐波导纳矩阵的具体内容为：
发电机：
Figure BDA0002773376590000049
Figure BDA00027733765900000410
Figure BDA00027733765900000411
为发电机谐波导纳矩阵，
Figure BDA00027733765900000412
为发电机第h次谐波导纳子相量矩阵，
Figure BDA00027733765900000413
为第k台发电机的基波电阻,
Figure BDA00027733765900000414
为第k台发电机的基波电抗，h为谐波次数，由步骤S1输入的谐波次数参数范围决定，
Figure BDA00027733765900000415
为第k台发电机第h次谐波的导纳，b为总发电机数；
负荷：
Figure BDA0002773376590000051
Figure BDA0002773376590000052
Figure BDA0002773376590000053
Figure BDA0002773376590000054
为负荷等效谐波导纳矩阵，
Figure BDA0002773376590000055
为负荷第h次谐波等效谐波导纳子相量矩阵，
Figure BDA0002773376590000056
为第z个负荷等效基波电阻，
Figure BDA0002773376590000057
为第z个负荷等效基波电抗，h为谐波次数，2≤h≤25，h为步骤S1输入的谐波次数参数范围决定，p为负荷数；
输电线路：
Figure BDA0002773376590000058
Figure BDA0002773376590000059
Figure BDA00027733765900000510
Figure BDA00027733765900000511
Figure BDA00027733765900000512
Figure BDA00027733765900000513
Figure BDA00027733765900000514
m为总支路数；
关联矩阵：CLfh线路f端谐波索引矩阵，CLth线路t端谐波索引矩阵，由线路节点编号生成：
Figure 1
Figure BDA0002773376590000062
Figure BDA0002773376590000063
表示第i条线路的电阻,
Figure BDA0002773376590000064
表示第i条线路的电抗，
Figure BDA0002773376590000065
表示第 i条线路t端电纳，
Figure BDA0002773376590000066
表示第i条线路f端电纳,
Figure BDA0002773376590000067
表示第i 条线路第h次谐波t端自导纳，
Figure BDA0002773376590000068
表示第i条线路第h次谐波f 端自导纳，
Figure BDA0002773376590000069
表示第i条线路第h次谐波t端对f端互导纳，
Figure BDA00027733765900000610
表示第i条线路第h次谐波f端对t端互导纳，
Figure BDA00027733765900000611
为第h 次谐波f端自导纳相量矩阵，
Figure BDA00027733765900000612
为第h次谐波t端自导纳相量矩阵，
Figure BDA00027733765900000613
为第h次谐波f端对t端互导纳相量矩阵,
Figure BDA00027733765900000614
为第h次谐波t端对f端互导纳相量矩阵，f表示首端，t表示末端；
变压器：
Figure BDA00027733765900000615
Figure BDA00027733765900000616
关联矩阵：CTfh表示变压器f端谐波索引矩阵，CTth表示变压器t端谐波索引矩阵，由节点相对编号顺序生成；
Figure 100002_9
Figure 100002_10
Figure BDA0002773376590000072
Figure BDA0002773376590000073
表示第q台变压器的基波电阻，
Figure BDA0002773376590000074
表示第q台变压器的基波电抗,
Figure BDA0002773376590000075
表示第q台变压器第h次谐波t端自导纳，
Figure BDA0002773376590000076
表示第q台变压器第h次谐波f端自导纳，
Figure BDA0002773376590000077
表示第q台变压器第 h次谐波t端对f端互导纳，
Figure BDA0002773376590000078
表示第q台变压器第h次谐波f 端对t端互导纳
Figure BDA0002773376590000079
为变压器第h次谐波f端自导纳相量矩阵，
Figure BDA00027733765900000710
为变压器第h次谐波t端自导纳相量矩阵，
Figure BDA00027733765900000711
为变压器第h次谐波f端对t端互导纳相量矩阵；
Figure BDA00027733765900000712
为变压器第h次谐波t端对f端互导纳相量矩阵f表示首端，t表示末端，h为谐波次数；
谐波总导纳矩阵：
Figure BDA00027733765900000713
进一步地，步骤S3中所述计算基波潮流的具体内容为：
Figure BDA00027733765900000714
Figure BDA0002773376590000081
Figure BDA0002773376590000082
Figure BDA0002773376590000083
Figure BDA0002773376590000084
Figure BDA0002773376590000085
Figure BDA0002773376590000086
Figure BDA0002773376590000087
Figure BDA0002773376590000088
Figure BDA0002773376590000089
Figure BDA00027733765900000810
[Ix]为电流实部矩阵，[Iy]为电流虚部矩阵，[Ux]为电压实部矩阵， [Uy]为电压虚部矩阵，[Pin]为节点注入有功功率矩阵，[Qin]为节点注入无功功率矩阵，[B]为电纳矩阵，[G]为电导矩阵电网运行中发电机母线电压均是受控的，电压受控的母线处无功功率是不可控，因此对于基础雅可比方程修订如下：
Figure BDA0002773376590000091
PVBUS、PQBUS为向量矩阵，由初始输入节点命名顺序决定；
不平衡量生成：
Figure BDA0002773376590000092
雅可比
Figure RE-GDA0002953358030000103
误差修正量：
[Δx]＝[J]\[F]
利用上述误差修正量修正电压：
[Ux](PQBUS)＝[Ux](PQBUS)+[Δx](i1)
[Uy](PQBUS)＝[Uy](PQBUS)+[Δx](i2)
[Ux](PVBUS)＝[Ux](PVBUS)+[Δx](i3)
[Uy](PVBUS)＝[Ux](PVBUS)+[Δx](i4)
i1为PQBUS电压实部的索引向量，i2为PQBUS电压虚部的索引向量， i3为PVBUS电压实部的索引向量,i4为PVBUS电压虚部的索引向量，由初始输入节点命名顺序决定；
判断F是否大于步骤S1给定的基波迭代误差值，若是，则再次执行上述修正潮流方程；若否，则结束基波潮流计算，得到基波潮流电压值，进而得到电流和功率值。
进一步地，所述步骤S4的具体内容为：
利用参数标记检测含有谐波源标记的节点，并判断是否有多个谐波源在同一节点；若有，则考虑谐波源之间的相互影响；若无，则不考虑谐波源之间的相互影响；
各次谐波电流：
Figure BDA0002773376590000101
若同一节点有多个谐波源时考虑它们的相互影响，由国标得到考虑谐波源之间相互影响的公式：
Figure BDA0002773376590000102
harmonich为第h次谐波电流含有率，Ix、Iy由基波潮流得到，
Figure BDA0002773376590000103
分别为谐波源1与谐波源2的谐波电流，Kh为给定系数；
由
Figure BDA0002773376590000104
得到各节点各次谐波电压，由谐波电压得到谐波损耗,利用谐波源节点的谐波损耗更新此节点的功率。
进一步地，步骤S6中判断前后节点注入功率误差是否小于迭代误差限；具体内容为：
Figure BDA0002773376590000111
Figure BDA0002773376590000112
为第k次循环得到的功率值矩阵，
Figure BDA0002773376590000113
为第k-1次循环得到的功率值矩阵；判断[ΔS(k)]是否小于步骤S1给定的谐波迭代误差值，若否，则返回步骤S3，若是，则结束谐波潮流，输出谐波电流，谐波电压，谐波损耗。
与现有技术相比，本发明具有以下有益效果：
本发明可以将谐波潮流计算中复杂高度稀疏的谐波阻抗矩阵生成转化和矩阵分解计算与计算机底层的数据结构组织，稀疏矩阵寻址等过程解耦开。这样一方面可以持续运用最先进的通用矩阵运算算法和商业程序来提高计算效率，同时也大幅度降低了谐波潮流算法扩展新模型和根据配网实际应用场景快速迭代增加新功能的效率。因此，本基于相量矩阵运算的谐波潮流运算具有在工程实际中广阔的应用空间。
附图说明
图1为本发明实施例的方法流程图。
图2为IEEE 3机9节点接线图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示，本实施例提供一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，包括以下步骤：
步骤S1：设置谐波次数，设置谐波源统计模型初始值，算法迭代误差限；
步骤S2：由包括变压器参数、线路参数、节点电压、发电机参数和负荷参数生成基波和谐波导纳矩阵；
步骤S3：计算基波潮流；
步骤S4：检测存在谐波源的节点，利用基波潮流结果和谐波源统计模型预测谐波数据得到谐波电流，然后再进行谐波潮流计算，得到谐波损耗；
步骤S5：利用谐波网损功率更新基波负荷功率；即利用步骤S4 得到的谐波损耗更新谐波源对应节点的功率，重复执行步骤S3，重新进行基波潮流计算；
步骤S6：判断前后节点注入功率误差是否小于迭代误差限，若是则计算收敛，输出谐波潮流计算结果，否则返回步骤S3。
如图2所示，以IEEE3机9节点电力系统为例进一步说明本实施例，表1给出接入电弧炉谐波源的模型信息，表2给出电铁谐波源的模型信息，表3给出光伏谐波源的模型信息，将电弧炉与电铁接在同一节点，系统其他信息与IEEE标准例程相同，不在赘述。
本实施例采用标幺值进行计算，输出结果以标幺值形式显示。
表1：电弧炉谐波电流理论值(Ii表示谐波电流含有率)
h(谐波次数) 2 3 4 5 6 7 8 9 10 11 12 13
I<sub>i</sub>(％) 5 7 4 5 4 2 2 1 1.5 1.5 0.5 0.5
表2：电铁谐波电流理论值
Figure BDA0002773376590000131
表3：光伏谐波电流理论值
Figure BDA0002773376590000132
Figure BDA0002773376590000141
表4.输出结果
Figure BDA0002773376590000142
在本实施例中，所述步骤S1的具体内容为：
设置要计算的谐波次数h为谐波次数，2≤h≤25，输入谐波源的统计模型的谐波电流参数，以及基波、谐波算法的迭代收敛误差值，小于10-3。
在本实施例中，步骤S2中所述建立基波导纳矩阵的具体内容为：
步骤S2中所述建立基波导纳矩阵的具体内容为：
Figure BDA0002773376590000151
Figure BDA0002773376590000152
Figure BDA0002773376590000153
表示第i条线路的电阻、
Figure BDA0002773376590000154
表示第i条线路的电抗,
Figure BDA0002773376590000155
表示第i 条线路的t端电纳，
Figure BDA0002773376590000156
表示第i条线路的f端电纳；
Figure BDA0002773376590000157
表示第 i条线路t端的自导纳,
Figure BDA0002773376590000158
表示第i条线路f端的自导纳，
Figure BDA0002773376590000159
表示第i条线路f端对t端的互导纳，
Figure BDA00027733765900001510
表示第i条线路t端对f端的互导纳；i表示第i条线路，f表示首端，t表示末端；
关联矩阵：Cfl表示f端索引矩阵，Ctl表示t端索引矩阵，为依据节点编号生成的矩阵；
Figure RE-GDA00029533580300001611
m表示总线路数,[*]表示矩阵；
基波变压器：
Figure BDA00027733765900001512
Figure BDA00027733765900001513
Figure BDA00027733765900001514
表示第q台变压器的电阻，
Figure BDA00027733765900001515
表示第q台变压器的电抗，tapq表示第q台变压器变比，
Figure BDA00027733765900001516
表示第q台变压器t端自导纳，
Figure BDA00027733765900001517
表示第q台变压器f端自导纳；
Figure BDA00027733765900001518
表示第q台变压器f端对t端互导纳，
Figure BDA00027733765900001519
表示第q台变压器t端对f端互导纳，f表示首端，t 表示末端，q表示变压器编号；
关联矩阵：CTf表示f端索引矩阵，CTt表示t端索引矩阵，根据变压器节点编号顺序生成的矩阵；
Figure BDA0002773376590000161
n表示总变压器台数，[*]表示矩阵；
基波总导纳：
Figure BDA0002773376590000162
在本实施例中，步骤S2中所述建立谐波导纳矩阵的具体内容为：
发电机：
Figure BDA0002773376590000163
Figure BDA0002773376590000164
Figure BDA0002773376590000165
为发电机谐波导纳矩阵，
Figure BDA0002773376590000166
为发电机第h次谐波导纳子相量矩阵，
Figure BDA0002773376590000167
为第k台发电机的基波电阻,
Figure BDA0002773376590000168
为第k台发电机的基波电抗，h为谐波次数，由步骤S1输入的谐波次数参数范围决定，
Figure BDA0002773376590000169
为第k台发电机第h次谐波的导纳，b为总发电机数；
负荷：
Figure BDA00027733765900001610
Figure BDA00027733765900001611
Figure BDA00027733765900001612
Figure BDA0002773376590000171
为负荷等效谐波导纳矩阵，
Figure BDA0002773376590000172
为负荷第h次谐波等效谐波导纳子相量矩阵，
Figure BDA0002773376590000173
为第z个负荷等效基波电阻，
Figure BDA0002773376590000174
为第z个负荷等效基波电抗，h为谐波次数，2≤h≤25，h为步骤S1输入的谐波次数参数范围决定，p为负荷数；
输电线路：
Figure BDA0002773376590000175
Figure BDA0002773376590000176
Figure BDA0002773376590000177
Figure BDA0002773376590000178
Figure BDA0002773376590000179
Figure BDA00027733765900001710
Figure BDA00027733765900001711
m为总支路数；
关联矩阵：CLfh线路f端谐波索引矩阵，CLthCLfh线路t端谐波索引矩阵，由线路节点编号生成：
Figure 2
Figure BDA0002773376590000181
Figure BDA0002773376590000182
表示第i条线路的电阻,
Figure BDA0002773376590000183
表示第i条线路的电抗，
Figure BDA0002773376590000184
表示第 i条线路t端电纳，
Figure BDA0002773376590000185
表示第i条线路f端电纳,
Figure BDA0002773376590000186
表示第i 条线路第h次谐波t端自导纳，
Figure BDA0002773376590000187
表示第i条线路第h次谐波f 端自导纳，
Figure BDA0002773376590000188
表示第i条线路第h次谐波t端对f端互导纳，
Figure BDA0002773376590000189
表示第i条线路第h次谐波f端对t端互导纳，
Figure BDA00027733765900001810
为第h 次谐波f端自导纳相量矩阵，
Figure BDA00027733765900001811
为第h次谐波t端自导纳相量矩阵，
Figure BDA00027733765900001812
为第h次谐波f端对t端互导纳相量矩阵,
Figure BDA00027733765900001813
为第h次谐波t端对f端互导纳相量矩阵，f表示首端，t表示末端；
变压器：
Figure BDA00027733765900001814
Figure BDA00027733765900001815
关联矩阵：CTfh表示变压器f端谐波索引矩阵，CTth表示变压器t端谐波索引矩阵，由节点相对编号顺序生成；
Figure 6
Figure 7
Figure BDA0002773376590000191
Figure BDA0002773376590000192
表示第q台变压器的基波电阻，
Figure BDA0002773376590000193
表示第q台变压器的基波电抗,
Figure BDA0002773376590000194
表示第q台变压器第h次谐波t端自导纳，
Figure BDA0002773376590000195
表示第q台变压器第h次谐波f端自导纳，
Figure BDA0002773376590000196
表示第q台变压器第 h次谐波t端对f端互导纳，
Figure BDA0002773376590000197
表示第q台变压器第h次谐波f 端对t端互导纳，
Figure BDA0002773376590000198
为变压器第h次谐波f端自导纳相量矩阵，
Figure BDA0002773376590000199
为变压器第h次谐波t端自导纳相量矩阵，
Figure BDA00027733765900001910
为变压器第h次谐波f端对t端互导纳相量矩阵；
Figure BDA00027733765900001911
为变压器第h 次谐波t端对f端互导纳相量矩阵f表示首端，t表示末端，h 为谐波次数；
谐波总导纳矩阵：
Figure BDA00027733765900001912
在本实施例中，步骤S3中所述计算基波潮流的具体内容为：
Figure BDA00027733765900001913
Figure BDA0002773376590000201
Figure BDA0002773376590000202
Figure BDA0002773376590000203
Figure BDA0002773376590000204
Figure BDA0002773376590000205
Figure BDA0002773376590000206
Figure BDA0002773376590000207
Figure BDA0002773376590000208
Figure BDA0002773376590000209
Figure BDA00027733765900002010
[Ix]为电流实部矩阵，[Iy]为电流虚部矩阵，[Ux]为电压实部矩阵， [Uy]为电压虚部矩阵，[Pin]为节点注入有功功率矩阵，[Qin]为节点注入无功功率矩阵，[B]为电纳矩阵，[G]为电导矩阵电网运行中发电机母线电压均是受控的，电压受控的母线处无功功率是不可控，因此对于基础雅可比方程修订如下：
Figure BDA0002773376590000211
PVBUS、PQBUS为向量矩阵，由初始输入节点命名顺序决定；
不平衡量生成：
Figure BDA0002773376590000212
雅可比
Figure RE-RE-GDA0002953358030000103
误差修正量：
[Δx]＝[J]\[F]
利用上述误差修正量修正电压
[Ux](PQBUS)＝[Ux](PQBUS)+[Δx](i1)
[Uy](PQBUS)＝[Uy](PQBUS)+[Δx](i2)
[Ux](PVBUS)＝[Ux](PVBUS)+[Δx](i3)
[Uy](PVBUS)＝[Ux](PVBUS)+[Δx](i4)
i1为PQBUS电压实部的索引向量，i2为PQBUS电压虚部的索引向量， i3为PVBUS电压实部的索引向量,i4为PVBUS电压虚部的索引向量，由初始输入节点命名顺序决定；
判断F是否大于步骤S1给定的基波迭代误差值，若是，则再次执行上述修正潮流方程；若否，则结束基波潮流计算，得到基波潮流电压值，进而得到电流和功率值。
在本实施例中，所述步骤S4的具体内容为：
利用参数标记检测含有谐波源标记的节点，并判断是否有多个谐波源在同一节点；若有，则考虑谐波源之间的相互影响；若无，则不考虑谐波源之间的相互影响；
各次谐波电流：
Figure BDA0002773376590000223
若同一节点有多个谐波源时考虑它们的相互影响，由国标得到考虑谐波源之间相互影响的公式：
Figure BDA0002773376590000221
harmonich为第h次谐波电流含有率，为步骤S1谐波源统计模型参数给定，Ix、Iy由基波潮流得到，
Figure BDA0002773376590000222
分别为谐波源1与谐波源2 的谐波电流，Kh为给定系数；
由
Figure BDA0002773376590000231
得到各节点各次谐波电压，由谐波电压得到谐波损耗,利用谐波源节点的谐波损耗更新此节点的功率。
在本实施例中，步骤S6中判断前后节点注入功率误差是否小于迭代误差限；具体内容为：
Figure BDA0002773376590000232
Figure BDA0002773376590000233
为第k次循环得到的功率值矩阵，
Figure BDA0002773376590000234
为第k-1次循环得到的功率值矩阵；判断[ΔS(k)]是否小于步骤S1给定的谐波迭代误差值，若否，则返回步骤S3，若是，则结束谐波潮流，输出谐波电流，谐波电压，谐波损耗。
较佳的，本实施例基于相量矩阵的思维进行程序架构，实现了大规模电网谐波潮流仿真的快速开发和高效计算，并充分运用多重索引矩阵等索引技术极大的简化了相应算法的复杂性，同时可以充分运用高度优化的成熟矩阵计算程序来实现算法的速度优化；基于高度便捷的模型扩展，面向矩阵和相量以及算法框架中谐波源的接口设计，实现了统计类谐波源模型的高度可扩展性，实现了模型扩展与算法主体维护的解耦，大幅度提高了算法程序不断迭代提升的效率；基于兼容树形网络和环形网络的潮流算法，适用电网末端接入各类分布式电源的情况也能很好的适应，能快速适应目前日益复杂的配电网设备接入后的仿真。
较佳的，本实施例基于相量矩阵的思维进行程序架构，实现了大规模电网谐波潮流仿真的快速开发和高效计算。通过矩阵相量的方式来描述，并充分运用多重索引矩阵等索引技术极大的简化了相应算法的复杂性，同时可以充分运用高度优化的成熟矩阵计算程序来实现算法的速度优化。
较佳的，本实施例实现了高度便捷的模型扩展，基于算法面向矩阵和相量以及算法框架中谐波源的接口设计，实现了统计类谐波源模型的高度可扩展性，实现了模型扩展与算法主体维护的解耦，大幅度提高了算法程序不断迭代提升的效率。
较佳的，本实施例本算法兼容树形网络和环形网络，适用电网末端接入各类分布式电源的情况也能很好的适应，能快速适应目前日益复杂的配电网设备接入后的仿真。
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。

Claims (7)

1.一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：包括以下步骤：
步骤S1：设置谐波次数，设置谐波源统计模型初始值，算法迭代误差限；
步骤S2：由包括变压器参数、线路参数、节点电压、发电机参数和负荷参数生成基波和谐波导纳矩阵；
步骤S3：计算基波潮流；
步骤S4：检测存在谐波源的节点，利用基波潮流结果和谐波源统计模型预测谐波数据得到谐波电流，然后再进行谐波潮流计算，得到谐波损耗；
步骤S5：利用谐波网损功率更新基波负荷功率；即利用步骤S4得到的谐波损耗更新谐波源对应节点的功率，重复执行步骤S3，重新进行基波潮流计算；
步骤S6：判断前后节点注入功率误差是否小于迭代误差限，若是则计算收敛，输出谐波潮流计算结果，否则返回步骤S3。
2.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：所述步骤S1的具体内容为：
设置要计算的谐波次数h为谐波次数，2≤h≤25，输入谐波源的统计模型的谐波电流参数，以及基波、谐波算法的迭代收敛误差值，小于10-3。
3.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：步骤S2中所述建立基波导纳矩阵的具体内容为：
基波线路：
Figure RE-FDA0002953358020000021
rl i表示第i条线路的电阻、
Figure RE-FDA0002953358020000022
表示第i条线路的电抗；
Figure RE-FDA0002953358020000023
表示第i条线路的t端电纳，
Figure RE-FDA0002953358020000024
表示第i条线路的f端电纳，
Figure RE-FDA0002953358020000025
表示第i条线路t端的自导纳,
Figure RE-FDA0002953358020000026
表示第i条线路f端的自导纳，
Figure RE-FDA0002953358020000027
表示第i条线路f端对t端的互导纳，
Figure RE-FDA0002953358020000028
表示第i条线路t端对f端的互导纳；i表示第i条线路，f表示首端，t表示末端；
关联矩阵：Cfl表示f端索引矩阵，Ctl表示t端索引矩阵，为依据节点编号生成的矩阵；
Figure RE-FDA0002953358020000029
m表示总线路数,[*]表示矩阵；
基波变压器：
Figure RE-FDA0002953358020000031
Figure RE-FDA0002953358020000032
Figure RE-FDA0002953358020000033
Figure RE-FDA0002953358020000034
表示第q台变压器的电阻，
Figure RE-FDA0002953358020000035
表示第q台变压器的电抗，tapq表示第q台变压器变比，
Figure RE-FDA0002953358020000036
表示第q台变压器t端自导纳，
Figure RE-FDA0002953358020000037
表示第q台变压器f端自导纳；
Figure RE-FDA0002953358020000038
表示第q台变压器f端对t端互导纳，
Figure RE-FDA0002953358020000039
表示第q台变压器t端对f端互导纳，f表示首端，t表示末端，q表示变压器编号；
关联矩阵：CTf表示f端索引矩阵，CTt表示t端索引矩阵，根据变压器节点编号顺序生成的矩阵；
Figure RE-FDA00029533580200000310
n表示总变压器台数，[*]表示矩阵；
基波总导纳：
Figure RE-FDA0002953358020000041
4.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：步骤S2中所述建立谐波导纳矩阵的具体内容为：
发电机：
Figure RE-FDA0002953358020000042
Figure RE-FDA0002953358020000043
Figure RE-FDA0002953358020000044
为发电机谐波导纳矩阵，
Figure RE-FDA0002953358020000045
为发电机第h次谐波导纳子矩阵，
Figure RE-FDA0002953358020000046
为第k台发电机的基波电阻,
Figure RE-FDA0002953358020000047
为第k台发电机的基波电抗，h为谐波次数，由步骤S1输入的谐波次数参数范围决定，
Figure RE-FDA0002953358020000048
为第k台发电机第h次谐波的导纳，b为总发电机数；
负荷：
Figure RE-FDA0002953358020000049
Figure RE-FDA00029533580200000410
Figure RE-FDA00029533580200000411
Figure RE-FDA00029533580200000412
为负荷等效谐波导纳矩阵，
Figure RE-FDA00029533580200000413
为负荷第h次谐波等效谐波导纳子相量矩阵，
Figure RE-FDA00029533580200000414
为第z个负荷等效基波电阻，
Figure RE-FDA00029533580200000415
为第z个负荷等效基波电抗，h为谐波次数，2≤h≤25，h为步骤S1输入的谐波次数参数范围决定，p为负荷数；
输电线路：
Figure RE-FDA0002953358020000051
Figure RE-FDA0002953358020000052
Figure RE-FDA0002953358020000053
Figure RE-FDA0002953358020000054
Figure RE-FDA0002953358020000055
Figure RE-FDA0002953358020000056
Figure RE-FDA0002953358020000057
m为总支路数；
关联矩阵：CLfh线路f端谐波索引矩阵，CLth线路t端谐波索引矩阵，由线路节点编号生成：
Figure 9
Figure 10
Figure RE-FDA0002953358020000063
rl i表示第i条线路的电阻,
Figure RE-FDA0002953358020000064
表示第i条线路的电抗，
Figure RE-FDA0002953358020000065
表示第i条线路t端电纳，
Figure RE-FDA0002953358020000066
表示第i条线路f端电纳,
Figure RE-FDA0002953358020000067
表示第i条线路第h次谐波t端自导纳，
Figure RE-FDA0002953358020000068
表示第i条线路第h次谐波f端自导纳，
Figure RE-FDA0002953358020000069
表示第i条线路第h次谐波t端对f端互导纳，
Figure RE-FDA00029533580200000610
表示第i条线路第h次谐波f端对t端互导纳，
Figure RE-FDA00029533580200000611
为第h次谐波f端自导纳相量矩阵，
Figure RE-FDA00029533580200000612
为第h次谐波t端自导纳相量矩阵，
Figure RE-FDA00029533580200000613
为第h次谐波f端对t端互导纳相量矩阵,
Figure RE-FDA00029533580200000614
为第h次谐波t端对f端互导纳相量矩阵，f表示首端，t表示末端；
变压器：
Figure RE-FDA0002953358020000071
Figure RE-FDA0002953358020000072
Figure RE-FDA0002953358020000073
关联矩阵：CTfh表示变压器f端谐波索引矩阵，CTth表示变压器t端谐波索引矩阵，由节点相对编号顺序生成；
Figure 11
Figure RE-FDA0002953358020000075
Figure RE-FDA0002953358020000076
表示第q台变压器的基波电阻，
Figure RE-FDA0002953358020000077
表示第q台变压器的基波电抗,
Figure RE-FDA0002953358020000078
表示第q台变压器第h次谐波t端自导纳，
Figure RE-FDA0002953358020000079
表示第q台变压器第h次谐波f端自导纳，
Figure RE-FDA00029533580200000710
表示第q台变压器第h次谐波t端对f端互导纳，
Figure RE-FDA00029533580200000711
表示第q台变压器第h次谐波f端对t端互导纳,
Figure RE-FDA00029533580200000712
为变压器第h次谐波f端自导纳相量矩阵，
Figure RE-FDA00029533580200000713
为变压器第h次谐波t端自导纳相量矩阵，
Figure RE-FDA0002953358020000081
为变压器第h次谐波f端对t端互导纳相量矩阵；
Figure RE-FDA0002953358020000082
为变压器第h次谐波t端对f端互导纳相量矩阵f表示首端，t表示末端，h为谐波次数；
谐波总导纳矩阵：
Figure RE-FDA0002953358020000083
5.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：步骤S3中所述计算基波潮流的具体内容为：
Figure RE-FDA0002953358020000084
Figure RE-FDA0002953358020000091
[ΔP]＝[Pin]-[Ux].*[Ix]-[Uy].*[Iy]
[ΔQ]＝[Qin]+[Ux].*[Iy]-[Uy].*[Ix]
Figure RE-FDA0002953358020000092
[ΔP]＝[Pin]-[Ux].*([G][Ux]-[B][Uy])-[Uy].*([G][Uy]+[B][Ux])
＝[Pin]-[Ux].*[G][Ux]+[Ux].*[B][Uy]-[Uy].*[G][Uy]-[Uy].*[B][Ux][ΔQ]＝[Qin]+[Ux].*([G][Uy]+[B][Ux])-[Uy].*([G][Ux]-[B][Uy])
＝[Qin]+[Ux].*[G][Uy]+[Ux].*[B][Ux]-[Uy].*[G][Ux]+[Uy].*[B][Uy]
Figure RE-FDA0002953358020000093
Figure RE-FDA0002953358020000094
Figure RE-FDA0002953358020000095
Figure RE-FDA0002953358020000096
[Ix]为电流实部矩阵，[Iy]为电流虚部矩阵，[Ux]为电压实部矩阵，[Uy]为电压虚部矩阵，[Pin]为节点注入有功功率矩阵，[Qin]为节点注入无功功率矩阵，[B]为电纳矩阵，[G]为电导矩阵；
电网运行中发电机母线电压均是受控的，电压受控的母线处无功功率是不可控，因此对于基础雅可比方程修订如下：
Figure RE-FDA0002953358020000101
PVBUS、PQBUS为向量矩阵，由初始输入节点命名顺序决定；
不平衡量生成：
Figure RE-FDA0002953358020000102
雅可比
Figure 12
误差修正量：
[Δx]＝[J]\[F]
利用上述误差修正量修正电压
[Ux](PQBUS)＝[Ux](PQBUS)+[Δx](i1)
[Uy](PQBUS)＝[Uy](PQBUS)+[Δx](i2)
[Ux](PVBUS)＝[Ux](PVBUS)+[Δx](i3)
[Uy](PVBUS)＝[Ux](PVBUS)+[Δx](i4)
i1为PQBUS电压实部的索引向量，i2为PQBUS电压虚部的索引向量，i3为PVBUS电压实部的索引向量,i4为PVBUS电压虚部的索引向量，由初始输入节点命名顺序决定；
判断F是否大于步骤S1给定的基波迭代误差值，若是，则再次执行上述修正潮流方程；若否，则结束基波潮流计算，得到基波潮流电压值，进而得到电流和功率值。
6.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：所述步骤S4的具体内容为：
利用参数标记检测含有谐波源标记的节点，并判断是否有多个谐波源在同一节点；若有，则考虑谐波源之间的相互影响；若无，则不考虑谐波源之间的相互影响；
各次谐波电流：
Figure RE-FDA0002953358020000111
若同一节点有多个谐波源时考虑它们的相互影响，由国标得到考虑谐波源之间相互影响的公式：
Figure RE-FDA0002953358020000112
harmonich为第h次谐波电流含有率，Ix、Iy由基波潮流得到，
Figure RE-FDA0002953358020000113
分别为谐波源1与谐波源2的谐波电流，Kh为给定系数；
由
Figure RE-FDA0002953358020000121
得到各节点各次谐波电压，由谐波电压得到谐波损耗,利用谐波源节点的谐波损耗更新此节点的功率。
7.根据权利要求1所述的一种支持快速开发和模型扩展的电网谐波潮流相量矩阵计算方法，其特征在于：步骤S6中判断前后节点注入功率误差是否小于迭代误差限；具体内容为：
Figure RE-FDA0002953358020000122
Figure RE-FDA0002953358020000123
为第k次循环得到的功率值矩阵，
Figure RE-FDA0002953358020000124
为第k-1次循环得到的功率值矩阵；判断[ΔS(k)]是否小于步骤S1给定的谐波迭代误差值，若否，则返回步骤S3，若是，则结束谐波潮流，输出谐波电流，谐波电压，谐波损耗。
CN202011257225.7A 2020-11-11 2020-11-11 一种电网谐波潮流相量矩阵计算方法 Active CN112667962B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011257225.7A CN112667962B (zh) 2020-11-11 2020-11-11 一种电网谐波潮流相量矩阵计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011257225.7A CN112667962B (zh) 2020-11-11 2020-11-11 一种电网谐波潮流相量矩阵计算方法

Publications (2)

Publication Number Publication Date
CN112667962A true CN112667962A (zh) 2021-04-16
CN112667962B CN112667962B (zh) 2023-02-14

Family

ID=75404043

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011257225.7A Active CN112667962B (zh) 2020-11-11 2020-11-11 一种电网谐波潮流相量矩阵计算方法

Country Status (1)

Country Link
CN (1) CN112667962B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114267191A (zh) * 2021-12-10 2022-04-01 北京理工大学 缓解交通拥堵驾驶员控制系统、方法、介质、设备及应用

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002040067A (ja) * 2000-07-27 2002-02-06 Nissin Electric Co Ltd 次数間高調波検出方法
CN104659816A (zh) * 2015-03-13 2015-05-27 贵州电力试验研究院 一种基于改进粒子群算法的分布式电源接入配电系统优化配置方法
CN106849093A (zh) * 2016-11-29 2017-06-13 国网冀北电力有限公司承德供电公司 用于家庭能源的大型光伏电站集电系统谐波潮流计算方法
CN107017637A (zh) * 2017-04-21 2017-08-04 重庆大学 一种三相整流机组谐波功率的简化建模评测方法
CN110601198A (zh) * 2019-10-30 2019-12-20 国网浙江省电力有限公司宁波供电公司 计及谐波和电压不平衡约束的混合微电网优化运行方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002040067A (ja) * 2000-07-27 2002-02-06 Nissin Electric Co Ltd 次数間高調波検出方法
CN104659816A (zh) * 2015-03-13 2015-05-27 贵州电力试验研究院 一种基于改进粒子群算法的分布式电源接入配电系统优化配置方法
CN106849093A (zh) * 2016-11-29 2017-06-13 国网冀北电力有限公司承德供电公司 用于家庭能源的大型光伏电站集电系统谐波潮流计算方法
CN107017637A (zh) * 2017-04-21 2017-08-04 重庆大学 一种三相整流机组谐波功率的简化建模评测方法
CN110601198A (zh) * 2019-10-30 2019-12-20 国网浙江省电力有限公司宁波供电公司 计及谐波和电压不平衡约束的混合微电网优化运行方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵德勇: "两相静止坐标系下并网变流器若干关键技术的研究", 《万方数据学位论文库》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114267191A (zh) * 2021-12-10 2022-04-01 北京理工大学 缓解交通拥堵驾驶员控制系统、方法、介质、设备及应用

Also Published As

Publication number Publication date
CN112667962B (zh) 2023-02-14

Similar Documents

Publication Publication Date Title
CN112670987B (zh) 一种电网三相谐波潮流相量矩阵计算方法
Jiang et al. A reduced-space interior point method for transient stability constrained optimal power flow
CN107329044B (zh) 一种基于电弧暂态分量的配电网单相接地故障选线方法
CN102751720B (zh) 一种Flexible HVDC潮流计算模型及其计算方法
CN105226668A (zh) 一种用于upfc的选址和容量配置方法
CN103259472B (zh) 大型发电机组次同步振荡的反振荡自适应抑制系统
Seifi et al. Power flow study and comparison of FACTS: Series (SSSC), Shunt (STATCOM), and Shunt-Series (UPFC)
CN111797510A (zh) 一种新能源场站短路比的计算方法及系统
CN112667962B (zh) 一种电网谐波潮流相量矩阵计算方法
CN104899396A (zh) 一种修正系数矩阵的快速分解法潮流计算方法
CN106655195B (zh) 一种主动式配电网高频谐波潮流的计算方法
CN108808681A (zh) 基于混合注入模型的光伏并网潮流计算方法
CN114282383A (zh) 基于传输线路解耦的有源配电网电磁暂态并行仿真方法
CN111900718B (zh) 基于多级优化追赶变分迭代法的有源配电网动态仿真方法
CN110266045B (zh) 光伏发电系统lcl-t模块的降阶方法
CN108281967B (zh) 一种改善电力系统潮流计算收敛性的算法启动方法
CN111339624B (zh) 基于psasp和emtp/atp短路电流直流分量计算方法
Xu et al. System‐level dynamic phasor models of hybrid AC/DC microgrids suitable for real‐time simulation and small signal analysis
CN114696329A (zh) 一种中压不接地配电网的三相潮流统一计算方法
CN110046450B (zh) 一种适合含超导电缆电网的牛顿法潮流计算的初值选取方法
Selmi et al. A simple method for the steady state performances of self-excited induction generators
CN111898282A (zh) 一种改进的模块化多电平换流器戴维南等效建模方法
Mashhour et al. Three-phase backward/forward power flow solution considering three-phase distribution transformers
CN103928941A (zh) 一种发电机聚合的解析求解方法
Liu et al. Power flow algorithm and practical contingency analysis for distribution systems with distributed generation

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