一种含大量风机的电力系统短路电流的计算方法
【专利摘要】本发明公开一种含大量风机的电力系统短路电流的计算方法，所述计算方法包括：根据馈线Lij中接入N个风机后的电路确定等效电路，所述等效电路分为N+1段线路，且每一段线路中连接一个等效电阻；根据所述等效电路确定所述馈线Lij中支路电压的通解；根据所述支路电压的通解确定所述馈线Lij对节点导纳矩阵的对角元素及非对角元素的贡献；遍历电力系统中的所有馈线，根据各所述馈线对节点导纳矩阵的对角元素和非对角元素的贡献确定节点导纳矩阵Y；根据节点导纳矩阵Y确定节点阻抗矩阵Z；根据节点阻抗矩阵Z确定馈线Lij的短路电流ii?j。本发明含大量风机的电力系统短路电流的计算方法可避免节点导纳矩阵的维数随着风机接入数量的增加而增加。
【专利说明】
一种含大量风机的电力系统短路电流的计算方法
技术领域
[0001] 本发明涉及电力系统技术领域，特别是涉及一种含大量风机的电力系统短路电流 的计算方法。
【背景技术】
[0002] 风电的大规模开发与利用将会给电力系统运行、控制、保护、稳定分析等各方面带 来显著影响。数量众多的风机接入电力系统后，可能显著改变电力系统的短路电流分布。为 了分析风机接入对电力系统短路电流带来的影响，需要对含风机的电力系统进行短路电流 计算。
[0003] 其中，短路电流的计算方法包括(假设故障发生在第f个节点）：
[0004] 故障之后从第f个节点流进
[0005] 故障之后第i个节点的电压
[0006] 其中，Zlf表示节点阻抗矩阵Z中第i行第f列的元素，Zff表示节点阻抗矩阵Z中第f 行第f列的元素，?表示过渡电阻。
[0007] 在计算出故障后各节点的电压后，进而可确定馈线L1J馈线Llj的首段节点号为i， 末端节点号为j)的短路电流i i-j:
[0008
_____公式（丨3 );其中，Zi-j表示馈线Lij的阻抗。
[0009] 短路电流的计算关键在于获得节点阻抗矩阵Z，通常情况下，往往是先获得节点导 纳矩阵Y，然后在通过求逆运算得到节点阻抗矩阵Z。而目前节点导纳矩阵Y的实现方法为：
[0010] 阻抗为ZM的任意馈线L1P对节点导纳矩阵Y的对角元素及非对应元素的贡献分别
(14 >根据公式(14)遍历所有馈线，以确定节点导纳矩阵Y。
[0011]但是，在一个含大量风机的电力系统，设外电网的节点数目为nac。假设系统中接入 一座风电场，风电场中含有M条馈线，每条馈线上含有N个台风机(如图1所示），风电场总计 有风机F台。根据目前的电路电流技术方法：
[0012] (I)N台风机的接入点将馈线Llj分成了 (N+1)段，每一小段均被视作一条新的馈线。 N台风机的接入点均被视作新的节点。节点导纳矩阵Y的节点数增加 N。
[0013] (2)遍历这N+1段新的馈线，并根据公式（14)式，将这每一条馈线对节点导纳矩阵Y 的对角元素及非对角元素的贡献计入到系统的节点导纳矩阵Y当中。
[0014] (3)当F个风机接入后，系统节点的维数从Nac增加到了（NAC+F)维。相应地，节点导 纳矩阵Y的存储空间从·2增加到了（NAC+F)2。而节点导纳矩阵Y的求逆运算F 1的复杂度从〇 (Nac3)增加到o((NAG+F)3)，耗时也随之增加。
[0015]如果进一步考虑，系统中接入了多座风电场，那么系统节点维数、导纳矩阵存储空 间数、矩阵求逆复杂度将进一步扩大，从而大大增大计算机存储与计算开销。
【发明内容】
[0016] 本发明的目的是提供一种含大量风机的电力系统短路电流的计算方法，可避免节 点导纳矩阵的维数随着风机接入数量的增加而增加。
[0017] 为实现上述目的，本发明提供了如下方案：
[0018] -种含大量风机的电力系统短路电流的计算方法，所述计算方法包括：
[0019] 步骤一:根据馈线Llj中接入N个风机后的电路确定等效电路，所述等效电路分为N+ 1段线路，且每一段线路中连接一个等效电阻；其中，首段节点号为i，末端节点号为j，从首 段节点i到末端节点号j的每一段线路中等效电阻的阻抗依次为 Z1，...，zn，...，ZN+1，对应段 线路中接入风机的阻抗依次为X 1, . . .X，n，XN，，对应段线路的支路电压分别为
[0020] 步骤二:根据所述等效电路确定所述馈线Llj中支路电压g的通解；
[0021] 步骤三:根据所述支路电压&的通解确定所述馈线Llj对节点导纳矩阵Y的对角元 素及非对角元素的贡献；
[0022] 步骤四：遍历电力系统中的所有馈线，根据各所述馈线对节点导纳矩阵Y的对角元 素和非对角元素的贡献确定所述节点导纳矩阵Y;
[0023] 步骤五:根据所述节点导纳矩阵Y确定节点阻抗矩阵Z;
[0024] 步骤六:根据所述节点阻抗矩阵Z确定所述馈线Llj的短路电流i^。
[0025] 可选的，所述确定支路电压的通解的方法包括：
[0026] 根据以下公式确定支路电压& :
[0027]
____公式q ),其中，Ifj为第η段线路中流过的电流； [0028]将公式(1)变形得到公式(2):
[0029]
------公式（2 ):
[0030] 根据以下公式确定待定系数αη+1和βη+1:
[0031]
[0032]
[0033]
[0034] 通解。
[0035] 可选的，在步骤三中，根据以下公式计算所述馈线Llj对节点导纳矩阵的对角元素 及非对角元素的贡献：
[0036]
其中，yii表示对角元素，yij表示非对角元素。
[0037] 可选的，在步骤四中，所述确定节点导纳矩阵Y的方法包括：
[0038] 将各所述馈线对节点导纳矩阵的对角元素和非对角元素的贡献叠加求和。
[0039] 可选的，在步骤五中，根据以下公式所述确定节点阻抗矩阵Z:
[0040]
------公式(7)。
[0041] 可选的，，在步骤六中，根据以下公式确定所述馈线Llj的短路电流
[0042：
-----公式（8 )，其中，ZH为所述馈线Llj的阻抗，η为过渡阻 抗，Zlf为节点阻抗矩阵Z中第i行第f列的元素，Zjf为节点阻抗矩阵Z中第j行第f列的元素， Zff为节点阻抗矩阵Z中第f行第f列的元素。
[0043] 根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0044] 本发明含大量风机的电力系统短路电流的计算方法根据馈线Llj中接入N个风机后 的电路确定等效电路，进而确定所述馈线Llj中支路电压$的通解，根据所述支路电压&的 通解确定所述馈线Llj对节点导纳矩阵的对角元素及非对角元素的贡献，进而可确定节点导 纳矩阵Y，可避免节点导纳矩阵的维数随着风机接入数量的增加而增加。
【附图说明】
[0045] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所 需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施 例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图 获得其他的附图。
[0046]图1为馈线Lij中接入N台风机的电路图；
[0047]图2为图1的等效电路图；
[0048]图3为本发明含大量风机的电力系统短路电流的计算方法的流程图。
【具体实施方式】
[0049]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0050] 本发明的目的是提供一种含大量风机的电力系统短路电流的计算方法根据馈线 Llj中接入N个风机后的电路确定等效电路，进而确定所述馈线Llj中支路电压ξ的通解，根 据所述支路电压&的通解确定所述馈线Llj对节点导纳矩阵的对角元素及非对角元素的贡 献，进而可确定节点导纳矩阵Υ，可避免节点导纳矩阵的维数随着风机接入数量的增加而增 加。
[0051] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实 施方式对本发明作进一步详细的说明。
[0052]如图3所示，本发明含大量风机的电力系统短路电流的计算方法包括：
[0053]步骤110:根据馈线Llj中接入N个风机后的电路确定等效电路；
[0054] 步骤:120:根据所述等效电路确定所述馈线Llj中支路电压I；的通解；
[0055] 步骤130:根据所述支路电压g的通解确定所述馈线Llj对节点导纳矩阵Y的对角元 素及非对角元素的贡献；
[0056] 步骤140:遍历电力系统中的所有馈线，根据各所述馈线对节点导纳矩阵Y的对角 元素和非对角元素的贡献确定所述节点导纳矩阵Y;
[0057]步骤150:根据所述节点导纳矩阵Y确定节点阻抗矩阵Z;
[0058]步骤160:根据所述节点阻抗矩阵Z确定所述馈线Llj的短路电流i^。
[0059] 如图2所示，所述等效电路分为N+1段线路，且每一段线路中连接一个等效电阻;其 中，首段节点号为i，末端节点号为j，从首段节点i到末端节点号j的每一段线路中等效电阻 的阻抗依次为Z1, ...，zn，...，zN+1，对应段线路中接入风机的阻抗依次为Jf1, . .? ^ , 对应段线路的支路电压分别为，n = l，. . .，N-1。
[0060] 在步骤120中，所述确定支路电压的通解的方法包括：
[0061] 步骤121:根据以下公式确定支路电压#· η ·
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085]
[0086]
[0087]
[0088]
[0089]
[0090] 确定&的通解。
[0091] 其中，可根据公式(4)和(5)得到。
[0092]在步骤130中，根据以下公式计算所述馈线Llj对节点导纳矩阵的对角元素及非对 角元素的贡献：
[0093:
Μ >其中，yii表示对角元素，yij表示非对角元素。
[0094] 在步骤140中，所述确定节点导纳矩阵Y的方法包括:将各所述馈线对节点导纳矩 阵的对角元素和非对角元素的贡献叠加求和。
[0095] 在步骤150中，根据以下公式所述确定节点阻抗矩阵Z:
[0096] Z = Y-1--------公式（7)。
[0097]在步骤160中，根据以下公式确定所述馈线Llj的短路电流
[0098]
--公式（8 )，其中，ZH为所述馈线Llj的阻抗，rf为过渡阻 抗，Zlf为节点阻抗矩阵Z中第i行第f列的元素，Zjf为节点阻抗矩阵Z中第j行第f列的元素， Zff为节点阻抗矩阵Z中第f行第f列的元素。
[0099] 本发明方法在计算含风机的馈线Llj对应节点导纳矩阵Y的对角元素及非对角元素 时，没有将N台风机的接入点视作新的节点，因此不会导致节点导纳矩阵Y的节点数随风机 数量的增加而增加，进而节点导纳矩阵Y的存储空间数和矩阵求逆复杂度也不会随风机数 量的增加而增加，从而节约了计算机的存储空间及成本。
[0100] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他 实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
[0101] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说 明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据 本发明的思想，在【具体实施方式】及应用范围上均会有改变之处。综上所述，本说明书内容不 应理解为对本发明的限制。
【主权项】
1. 一种含大量风机的电力系统短路电流的计算方法，其特征在于，所述计算方法包括： 步骤一:根据馈线以^中接入N个风机后的电路确定等效电路，所述等效电路分为N+1段 线路，且每一段线路中连接一个等效电阻;其中，首段节点号为i，末端节点号为j，从首段节 点i到末端节点号j的每一段线路中等效电阻的阻抗依次为Z1，. . .，Zn，. . .，ZN+1，对应段线路 中接入风机的阻抗依次为义1，..，《V.，对应段线路的支路电压分别为巧，...义，...，巧、W， η 二 1,. . . , N-1; 步骤二:根据所述等效电路确定所述馈线中支路电压?ζ的通解； 步骤Ξ:根据所述支路电压的通解确定所述馈线以^对节点导纳矩阵Υ的对角元素及 非对角元素的贡献； 步骤四：遍历电力系统中的所有馈线，根据各所述馈线对节点导纳矩阵Υ的对角元素和 非对角元素的贡献确定所述节点导纳矩阵Υ; 步骤五:根据所述节点导纳矩阵Υ确定节点阻抗矩阵Ζ; 步骤六:根据所述节点阻抗矩阵Ζ确定所述馈线k北勺短路电流il-J。2. 根据权利要求1所述的大量风机的电力系统短路电流的计算方法，其特征在于，在步 骤二中，所述确定支路电压A的通解的方法包括： 根据W下公式确定支路电压：根据所述等效电路中的边界条件为g =1, 二0---------公或(5 )，确定ξ的通解。3. 根据权利要求1所述的大量风机的电力系统短路电流的计算方法，其特征在于，在步 骤Ξ中，根据W下公式计算所述馈线节点导纳矩阵的对角元素及非对角元素的贡献：其中，yi康示对角元素 ，yi康示非对角元素。4. 根据权利要求1所述的含大量风机的电力系统短路电流的计算方法，其特征在于，在 步骤四中，所述确定节点导纳矩阵Y的方法包括： 将各所述馈线对节点导纳矩阵的对角元素和非对角元素的贡献叠加求和。5. 根据权利要求1所述的含大量风机的电力系统短路电流的计算方法，其特征在于，在 步骤五中，根据W下公式所述确定节点阻抗矩阵Z: Ζ=γ-ι--------公式(7)。6. 根据权利要求1所述的含大量风机的电力系统短路电流的计算方法，其特征在于，在 步骤六中，根据W下公式确定所述馈线k北勺短路电流ii-j:-…--公乂（8)，其中，ZW为所述馈线Lu的阻抗，rf为过渡阻抗， zif为节点阻抗矩阵Z中第i行第f列的元素，ZW为节点阻抗矩阵Z中第j行第f列的元素，zff为 节点阻抗矩阵Z中第f行第f列的元素。
【文档编号】G01R19/00GK105842516SQ201610160725
【公开日】2016年8月10日
【申请日】2016年3月21日
【发明人】孙毅卫, 汪雨辰, 李珂, 陈西平, 陈来军, 黄少伟, 王志文, 王召健, 唐可翾
【申请人】陕西省地方电力（集团）有限公司, 清华大学