CN103715684A - 一种评估励磁涌流引发电压严重畸变风险的方法 - Google Patents

一种评估励磁涌流引发电压严重畸变风险的方法 Download PDF

Info

Publication number
CN103715684A
CN103715684A CN201310660744.1A CN201310660744A CN103715684A CN 103715684 A CN103715684 A CN 103715684A CN 201310660744 A CN201310660744 A CN 201310660744A CN 103715684 A CN103715684 A CN 103715684A
Authority
CN
China
Prior art keywords
harmonic
risk
voltage
inrush current
impedance
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
CN201310660744.1A
Other languages
English (en)
Other versions
CN103715684B (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Sichuan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Sichuan Electric Power Co Ltd
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Sichuan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201310660744.1A priority Critical patent/CN103715684B/zh
Publication of CN103715684A publication Critical patent/CN103715684A/zh
Application granted granted Critical
Publication of CN103715684B publication Critical patent/CN103715684B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种评估励磁涌流引发电压严重畸变风险的方法，包括以下步骤：（a）生成待求系统各次的谐波导纳阵；（b）求取电力系统各次谐波阻抗阵；（c）求出重点考察变压器的容量标幺值；（d）求取指标，包括3、5次谐波畸变率风险指标、电压总畸变率风险指标、电压瞬时值风险指标；（e）根据指标的风险判定。本发明采用上述方法，可以较直观地反映系统中因励磁涌流导致的谐波电压畸变风险，并寻找出系统中风险最大的厂站，该指标使用方便，可以与电磁暂态仿真方法相配合，为仿真提供前期的方案指导，并对仿真结果提供相应的理论参考。

Description

一种评估励磁涌流引发电压严重畸变风险的方法
技术领域
本发明涉及电网安全评估领域，具体是一种评估励磁涌流引发电压严重畸变风险的方法。 
背景技术
变压器正常运行时，励磁电流很小。但是当空载合闸到电网的瞬间，励磁电流可能急剧增加为正常励磁电流值得几十倍，甚至上百倍。这种空载合闸时出现的过电流现象称为励磁涌流。谐波是一个周期量的正弦波分量，其频率为基波频率的整数倍。由于非线性负荷的接入以及供电系统中本身存在的非线性元件，使得电力系统中会引入频率不等于50Hz的谐波分量，使得电力网中的电压不再是光滑的正弦波，从而发生畸变。 
    谐波互阻抗：电力系统的节点方程可以满足以下的形式： 
       
Figure 167094DEST_PATH_IMAGE002
其中，
Figure 867459DEST_PATH_IMAGE004
是节点电压的列向量，而
Figure 162436DEST_PATH_IMAGE006
为节点注入电流的列向量。
Figure 117836DEST_PATH_IMAGE008
则为系统的阻抗阵。同理，上式也可用来表征电力系统中n次谐波的电压电流关系，即：
       
Figure 498264DEST_PATH_IMAGE010
    
其中，是节点n次电压的列向量，而
Figure 886225DEST_PATH_IMAGE014
为节点n次注入电流的列向量。则为系统的n次阻抗阵，它可以通过下式表示：
       
Figure 2013106607441100002DEST_PATH_IMAGE017
  
式中则为节点i与节点j之间的n次谐波互阻抗。
    电力系统中当出现空载变压器合闸操作时，变压器会向系统中注入励磁涌流。分析表明，正常情况下变压器注入的励磁涌流可达额定电流的6~8倍。长期以来，励磁涌流对继电保护的影响受到了广泛的关注，许多学者研究了多种励磁涌流的识别方法，并将其用于变压器差动保护中实现制动功能，以防止励磁涌流引发变压器的保护误动。 
    然而，励磁涌流对于系统的危害并不单纯地表现在对继电保护的影响上。励磁涌流本身含有丰富的谐波分量，对电力系统的电能质量也会造成严重的污染。特别是在建设初期或者故障恢复期，电网本身就具有长线路、轻负荷的弱联系特性，系统中的电感电容参数接近低次谐波的谐振区间。此时，励磁涌流的注入很可能会在电网的末端产生较大的谐波电压，从而导致系统末端电压的畸变率明显恶化，而电压的瞬时值也可能因谐波电压的叠加而超出允许的范围。在这种工况下，设备以及负荷存在着很大的损毁风险。这一风险在偏远地区的运行电网中已多次发生。 
    在励磁涌流导致谐波电压畸变风险的分析方面，大部分学者主要采用以下三种方法： 
1）对系统中长线路的分布特性进行分析。分析认为，当线路长度接近谐波电磁波的1/4波长时，在线路末端谐波电压可能会被严重放大，从而造成严重的电压畸变，导致设备损毁。该方法能够较为直观地说明具有长线路、轻负荷特性的弱联系系统中励磁涌流引发谐波电压严重畸变的机理，然而针对一个较为复杂的网络，由于负荷、机组、变压器等元件与线路的耦合，长线路长度与末端谐波电压放大风险的关系变得不再明显，该方法则显得有所前欠缺。
2）对空充主变母线的谐振频率进行扫描。如果扫描结果发现电网内存在着某次谐波的谐振点，那么根据经验就可以判定出系统内某些节点可能会产生严重的谐波电压畸变，甚至出现较高的谐波过电压。该方法综合考虑了负荷、发电机等多种因素，然而却无法明确的给出系统内谐波电压畸变风险最高的母线。 
    3）利用PSCAD/EMTDC等电磁暂态分析软件，对系统进行详细的仿真研究。该方法可以很直观的发现系统中谐波电压畸变风险，并较为准确的计算出畸变率以及谐波过电压数值。然而，该方法缺少较好的理论指导，导致其针对性与可信性不强。 
发明内容
本发明提供了一种评估励磁涌流引发电压严重畸变风险的方法，目的是在基于网络的谐波互阻抗，并综合考虑空充变压器容量以及系统中线路、负荷、无功补偿等因素的基础上，可以较直观地反应系统中因励磁涌流导致的谐波电压畸变风险，并寻找出系统中风险最大的厂站问题。 
    本发明的目的通过下述技术方案实现：一种评估励磁涌流引发电压严重畸变风险的方法，包括以下步骤： 
    （a）生成待求系统各次的谐波导纳阵；
    （b）求取电力系统各次谐波阻抗阵；
    （c）求出重点考察变压器的容量标幺值：选取系统中容量较大、电压等级较高的变压器，将其作为重点考察的谐波注入源，查阅技术资料，获取这些重点考察变压器的容量
Figure 2013106607441100002DEST_PATH_IMAGE021
，并利用下式计算出容量的标幺值：
                                               
Figure 2013106607441100002DEST_PATH_IMAGE023
    （d）指标的求取
求出阻抗阵中的谐波互阻抗，计算出三类电压畸变风险指标：
(b)步骤中求得的谐波互阻抗Z nij 反应了电网中n次谐波的网络分布特性。利用该参数，可以对系统中励磁涌流引发的谐波电压畸变进行初步估算。
设定合空变操作时其变压器注入的励磁涌流峰值为额定电流的
Figure 2013106607441100002DEST_PATH_IMAGE025
倍。根据文献，变压器励磁涌流可达变压器额定电流的6~8倍。为考虑最恶劣的情况，取励磁涌流峰值最大的情况，即令，则在进行合空变操作时，其产生的励磁涌流的峰值可以通过下式进行计算。 
式中，为变压器的额定容量，
Figure 2013106607441100002DEST_PATH_IMAGE035
为变压器在操作侧的额定电压。 
由此在最恶劣的情况下，主变产生的励磁涌流中n次谐波的大小可通过下式计算得出： 
Figure 2013106607441100002DEST_PATH_IMAGE037
式中：
Figure 2013106607441100002DEST_PATH_IMAGE039
为励磁涌流中n次谐波有效值与涌流峰值的比值。
假设系统中没有其他谐波源，则该工况下母线j的n次谐波电压可由下式加以计算： 
Figure 2013106607441100002DEST_PATH_IMAGE041
利用求得的谐波电压，可得对励磁涌流引发的谐波电压畸变风险进行分析，主要体现在以下三个方面：
（1）单次谐波畸变率超标风险；
（2）电压总畸变率超标风险；
（3）电压最大瞬时值超标风险。
以下将针对这三个依次进行说明。 
（1）单次谐波畸变率超标风险及单次谐波畸变率风险指标 
由于励磁涌流引发的谐波电流持续时间相对较短，且谐波含量处于衰减状态，因此一般可以不按照《公共电网谐波》标准对各次谐波含有率进行限制。但是根据文献可知，如果系统中3、5次谐波电压含有率超过了20%，则系统中感应电动机可能在短时间内被损坏。因此，根据这一描述，需要对励磁涌流在系统中造成的3、5次谐波畸变进行分析。
根据前文可以算得，当母线i上注入励磁涌流时，母线j上可能出现的最大3、5次电压的含量D3、D5分别为： 
Figure 2013106607441100002DEST_PATH_IMAGE043
假设系统中各条母线的基波电压标幺值为1.0pu，因此可以算得，上式的标幺制形式为：
Figure 2013106607441100002DEST_PATH_IMAGE045
由此可以得到，在最恶劣的条件下，励磁涌流造成的3、5次电压畸变不会造成设备损坏风险的条件为：
Figure 2013106607441100002DEST_PATH_IMAGE047
即：
由此可以定义单次谐波畸变率风险指标
Figure 2013106607441100002DEST_PATH_IMAGE051
，该指标主要体现了系统中3、5次谐波畸变过大造成电动机损坏的风险，定义该指标如下：
Figure 2013106607441100002DEST_PATH_IMAGE053
（2） 电压总畸变率超标风险及电压总畸变率风险指标
由于我国对正弦波形畸变率极限值的定义都是基于稳态值的，对于持续时间短暂的谐波都没有规定。不过，文献给出了一个参考值，即对于谐波持续时间少于2s，间隔时间大于30s的谐波畸变率允许值可放宽至稳态值的2倍。同时，考虑到励磁涌流中谐波含量呈现逐渐衰减趋势，因此本文中对于35kV电网的短时谐波畸变率允许值设定为10%。
根据前文计算结果，可以得到当母线i上注入励磁涌流时，母线j上可能出现的最大电压畸变率THD为 
Figure 955033DEST_PATH_IMAGE054
同样，假设系统中各条母线的基波电压标幺值为1.0pu，因此可以算得，上式的标幺制形式为：
Figure 856256DEST_PATH_IMAGE056
若希望母线j的电压畸变率不超过标准，则需满足：
Figure 125825DEST_PATH_IMAGE058
即：
由此可以定义电压总畸变率风险指标，该指标主要体现了励磁涌流注入后，系统电压畸变率超过短时允许值的风险，定义该指标如下：
Figure 147023DEST_PATH_IMAGE064
（3）电压最大瞬时值超标风险及谐波电压畸变风险指标
此外，在励磁涌流向电网注入谐波后，还存在着因各次谐波叠加后，造成电压瞬时值偏高而击穿绝缘的风险。
为计算该风险，设定一种极端的情况，即各次谐波的峰值均在同一时刻出现，此时电压的瞬时值VP达到最大，为： 
Figure 271100DEST_PATH_IMAGE066
根据文献，本文选择电压瞬时值最大允许值为额定相电压峰值的1.5倍，因此可以得到瞬时值不超标的关系式为：
Figure 668845DEST_PATH_IMAGE068
假设系统中各条母线的基波电压标幺值为1.0pu，因此可以算得，上式的标幺制形式为：
Figure 222145DEST_PATH_IMAGE070
即：
Figure 855382DEST_PATH_IMAGE072
由此可以定义电压瞬时值超标风险指标，该指标主要体现了励磁涌流注入后，系统中各次谐波电压叠加后瞬时值超过允许值的风险，定义该指标如下：
Figure 402612DEST_PATH_IMAGE076
根据推导，上述三种谐波电压畸变风险指标
Figure 466645DEST_PATH_IMAGE078
反应了在系统母线i处对一台容量为
Figure 513361DEST_PATH_IMAGE080
的主变进行空载合闸操作时，在母线j处产生谐波电压畸变的风险。当系统中某条母线的
Figure DEST_PATH_IMAGE081
指标计算结果大于1时，说明该电网在某些恶劣的工况下，可能存在着对应的风险，
Figure 238128DEST_PATH_IMAGE078
值越大，则风险的程度就越高。
这些指标综合考虑了被操作主变的容量以及系统中发电机、负荷、线路以及变压器等元件对谐波电压畸变的影响。因此，它比单纯的谐波互阻抗Znij更能体现出不同变压器在空载合闸操作时对谐波电压畸变风险的影响。 
最后对总结出来的指标归纳如下： 
3、5次谐波畸变率风险指标：
Figure DEST_PATH_IMAGE083
电压总畸变率风险指标：
Figure DEST_PATH_IMAGE085
电压瞬时值风险指标：
Figure DEST_PATH_IMAGE087
    （e）根据指标的风险判定
    当系统中某条母线的
Figure DEST_PATH_IMAGE089
指标计算结果大于1时，说明该电网在某些恶劣的工况下，可能存在着对应的风险，
Figure 630142DEST_PATH_IMAGE089
值越大，则风险的程度就越大。
    进一步地，所述步骤（a）的具体过程为： 
（a1）将输电线路等效为π型等值电路，得出谐波阻抗：
Figure DEST_PATH_IMAGE091
，式中，
Figure DEST_PATH_IMAGE093
是n次谐波下线路的特征阻抗，
Figure DEST_PATH_IMAGE095
是n次谐波下线路的传播系数；
（a2）求出变压器的谐波等值阻抗，
Figure DEST_PATH_IMAGE097
，式中，
Figure DEST_PATH_IMAGE099
、
Figure DEST_PATH_IMAGE101
为变压器n次谐波的谐波电阻与电抗；
Figure DEST_PATH_IMAGE103
、
Figure DEST_PATH_IMAGE105
为变压器基次谐波的谐波电阻与电抗；
（a3）求出发动机的谐波等效阻抗：
Figure DEST_PATH_IMAGE107
，式中，
Figure DEST_PATH_IMAGE109
为电动机在n次谐波下的电抗，
Figure DEST_PATH_IMAGE111
为基波正序阻抗，n为谐波次数；
（a4）求出负荷的谐波阻抗：将负荷等值为一个综合等值电动机处理，其谐波阻抗
Figure DEST_PATH_IMAGE113
为
Figure DEST_PATH_IMAGE115
                        
式中，、
Figure DEST_PATH_IMAGE119
为等值电动机的基波负序电阻与电抗。
（a5）求出无功补偿元件的谐波阻抗：将无功补偿元件近似的等值为单一的电容、电感元件，因此感性补偿元件与容性补偿元件的谐波阻抗
Figure DEST_PATH_IMAGE121
、分别为： 
Figure DEST_PATH_IMAGE125
式中，
Figure DEST_PATH_IMAGE127
为基波角频率，
Figure DEST_PATH_IMAGE129
、
Figure DEST_PATH_IMAGE131
分别为感性补偿元件与容性补偿元件的电感电容值。
    进一步地，本方案设定励磁涌流的峰值为变压器额定电流的8.0倍，同时选取电流中各次谐波电流的有效值与励磁涌流峰值之间的比例为典型值。 
进一步地，所述步骤（d）中，由于励磁涌流引入的谐波属于暂态谐波，其持续时间最多为分钟级，因此上述风险的判断依据分别为：1）系统中3、5次谐波电压含有率不得超过20%；2）电网中短时谐波畸变率允许值不得超过10%；3）设定每次谐波的峰值均在同一时刻出现，瞬时值最大允许值不得超过额定峰值的1.5倍。 
本发明具有以下有益效果：本发明提出了励磁涌流引发谐波电压畸变风险指标
Figure 578810DEST_PATH_IMAGE132
，该指标基于网络的谐波互阻抗，并综合考虑了空充变压器容量，系统中线路、负荷、无功补偿等因素，可以较直观地反映系统中因励磁涌流导致的谐波电压畸变风险，并寻找出系统中风险最大的厂站。该指标使用方便，可以与电磁暂态仿真方法相配合，为仿真提供前期的方案指导，并对仿真结果提供相应的理论参考。 
附图说明
图1为本发明的实现步骤框图； 
图2为我国西南某藏区联网工程图；
图3为 2号母线主变合闸时10号母线电压瞬时值；
图4为2号母线主变合闸时18号母线电压瞬时值；
图5为4号母线主变合闸时10号母线电压瞬时值；
图6为4号母线主变合闸时18号母线电压瞬时值。
具体实施方式
下面结合实施例及附图，对本发明作进一步的详细说明，但本发明的实施方式不仅限于此。 
实施例1： 
    如图1所示，本实施例所述的一种评估励磁涌流引发电压严重畸变风险的方法，包括以下步骤：
    （a）生成待求系统各次的谐波导纳阵：
    （b）求取电力系统各次谐波阻抗阵；
    （c）求出重点考察变压器的容量标幺值：选取系统中容量较大、电压等级较高的变压器，将其作为重点考察的谐波注入源。查阅相关资料，获取这些重点考察变压器的容量
Figure DEST_PATH_IMAGE133
，并利用下式计算出容量的标幺值：
                                                  
Figure DEST_PATH_IMAGE135
    （d）指标的求取
求出阻抗阵中的谐波互阻抗，计算出三类电压畸变风险指标：
3、5次谐波畸变率风险指标：
Figure 806047DEST_PATH_IMAGE083
电压总畸变率风险指标：
Figure 493643DEST_PATH_IMAGE085
电压瞬时值风险指标：
Figure 404092DEST_PATH_IMAGE136
    （e）根据指标的风险判定
    当系统中某条母线的
Figure DEST_PATH_IMAGE137
指标计算结果大于1时，说明该电网在某些恶劣的工况下，可能存在着对应的风险，
Figure 68816DEST_PATH_IMAGE137
值越大，则风险的程度就越大。
    进一步地，所述步骤（a）的具体过程为： 
（a1）将输电线路等效为π型等值电路，得出谐波阻抗：
Figure 722914DEST_PATH_IMAGE138
，式中，
Figure 530595DEST_PATH_IMAGE093
是n次谐波下线路的特征阻抗，
Figure 346367DEST_PATH_IMAGE095
是n次谐波下线路的传播系数；
（a2）求出变压器的谐波等值阻抗，
Figure 848851DEST_PATH_IMAGE097
，式中，
Figure 41060DEST_PATH_IMAGE099
、
Figure 93461DEST_PATH_IMAGE101
为变压器n次谐波的谐波电阻与电抗；、
Figure 93351DEST_PATH_IMAGE105
为变压器基次谐波的谐波电阻与电抗；
（a3）求出发动机的谐波等效阻抗：
Figure 89251DEST_PATH_IMAGE107
，式中，
Figure 605946DEST_PATH_IMAGE109
为电动机在n次谐波下的电抗，
Figure 29099DEST_PATH_IMAGE111
为基波正序阻抗，n为谐波次数；
（a4）求出负荷的谐波阻抗：将负荷等值为一个综合等值电动机处理，其谐波阻抗
Figure 264033DEST_PATH_IMAGE113
为
Figure 63624DEST_PATH_IMAGE115
                        
式中，、
Figure 130396DEST_PATH_IMAGE119
为等值电动机的基波负序电阻与电抗。
（a5）求出无功补偿元件的谐波阻抗：将无功补偿元件近似的等值为单一的电容、电感元件，因此感性补偿元件与容性补偿元件的谐波阻抗
Figure 603360DEST_PATH_IMAGE121
、
Figure 442527DEST_PATH_IMAGE123
分别为： 
                           
式中，
Figure DEST_PATH_IMAGE139
为基波角频率，
Figure 449502DEST_PATH_IMAGE129
、
Figure 659029DEST_PATH_IMAGE131
分别为感性补偿元件与容性补偿元件的电感电容值。
    本实施例设定励磁涌流的峰值为变压器额定电流的8.0倍，同时选取电流中各次谐波电流的有效值与励磁涌流峰值之间的比例为典型值。 
本实施例的步骤（d）中，由于励磁涌流引入的谐波属于暂态谐波，其持续时间最多为分钟级，因此上述风险的判断依据分别为：1）系统中3、5次谐波电压含有率不得超过20%；2）电网中短时谐波畸变率允许值不得超过10%；3）设定每次谐波的峰值均在同一时刻出现，瞬时值最大允许值不得超过额定峰值的1.5倍。 
为了验证本文所提出指标的有效性，本节以我国西南某藏区联网工程为算例进行指标计算与风险分析。算例系统示意图如下图所示。图1中方框内所写的数字为母线的编号。 
图2所示的系统包含了500kV、220kV以及110kV三个电压等级，系统各220kV变电站之间采用双回线串供的输电方式。整个系统输电线路距离较长，其中500kV线路总长度为344km，220kV输电线路总长度为1414km。系统中110kV输电线路长度也很长，母线7与母线10之间输电线路距离也达到了320km以上。在负荷方面，由于藏区经济相对落后，负荷水平较低，在小方式下220kV5号母线下网负荷功率大约为10MW，而6号母线下网功率大约为20MW。由此可见，图中所示的系统是一个典型的长线路、轻负荷的弱联系系统，谐波过电压风险较高。 
在图2所示的系统中，分别以下述两个操作为例，对系统的谐波过电压风险指标进行计算： 
1）在500kV2号母线对第二台1000MVA空载变压器进行合闸操作；
2）在220kV4号母线对第二台120MVA空载变压器进行合闸操作。
本次计算中，忽略6次及以上的谐波分量，同时设定励磁涌流中各次谐波的含量k n 为如下表所示的典型值。 
  
表1 各次谐波含量
 
Figure 410209DEST_PATH_IMAGE140
在上述假设下，根据指标的定义，可得系统中指标值较大的几条母线对应的计算结果如下表所示。
表2 指标计算结果
 
Figure 2013106607441100002DEST_PATH_DEST_PATH_IMAGE075
通过表中风险指标的计算，我们可以轻易的做出到以下几点判断：
1）当在母线2、4号分别进行空载变压器合闸操作时，系统中计算出的风险指标均大于1，这说明系统中可能存在着因励磁涌流导致谐波电压严重畸变的风险。
2）同时由于母线2上并入的主变容量超过4号母线上并入主变容量的8倍以上，因此在2号母线上进行空载主变合闸操作其造成的风险比4号母线上造成的风险要大。 
3）通过对比
Figure DEST_PATH_IMAGE143
与两个指标的计算结果，可以发现当2号母线进行空载主变合闸操作时，图1所示的系统中5次谐波的放大水平比3次谐波要大很多。而在4号母线进行空载主变合闸操作时，两者的差距会有所减小。 
4）图2所示的系统中，谐波电压畸变风险较大的母线包括10号母线与9号母线，而8号母线与18号母线的风险则偏小一些。 
为验证计算得准确性，利用PSCAD搭建了相同的模型。为说明特点，选取表1中风险最大的10号母线以及相对风险较小的18号母线进行分析。分别设置在2号、4号母线A相电压瞬时值为0时对各自的主变进行合闸操作，忽略主变的剩磁，可以得到10、18号两条母线35kV侧所仿真得到的电压瞬时值如图3~图6所示。 
对仿真结果进行分析，可以得到当在2号、4号母线进行空载主变合闸操作时，10号与18号母线35kV侧出现的最高电压瞬时值、最大电压畸变率以及最大的3、5次谐波畸变率如下表所示。 
表3  仿真结果分析 
 
Figure DEST_PATH_IMAGE146
通过图3~图6及表3的对比，我们可以得到以下判断：
1）当2号、4号分别空充主变时，10号、18号母线出现了明显的电压畸变。其中2号主变空充时其最高电压畸变率以及最高电压瞬时值均明显高于4号主变空充的时候。
2）对比10号母线以及18号母线的结果，10号母线的电压峰值以及畸变率也超过了18号母线。 
3）在2号母线空充主变时，母线10与18上的5次谐波畸变明显大于3次谐波，而在4号母线空充时，3次谐波与5次谐波的畸变率没有明显的区别。 
这些分析结果与表2计算结果相符，从而说明了本文所提出的指标确实对电磁暂态仿真有重要的指导意义。 
由此可见，图2所示的联网工程存在着较大的因励磁涌流导致谐波电压畸变的风险，包括电压瞬时值超标、单次谐波畸变率超标以及谐波电压总畸变超标等，建议采取相应的措施进行抑制。 
值得说明的是，为简化计算，谐波过电压风险指标计算时假定励磁涌流与额定电流的倍数为最严重的8.0倍。同时各次谐波的含量
Figure DEST_PATH_IMAGE149
也只能取典型值。但是由于实际操作中，励磁电流与额定电流的倍数会随着系统阻抗等因素而变化，且每次谐波的含量也会有所不同。因此该指标无法精确地获得谐波过电压的数值。比如，算例计算中出现了风险指标超过1，但是仿真得出的结果却没有超标的情况。这是由于风险计算时考虑的是最恶劣的情况，而仿真时励磁涌流的峰值很可能没有达到额定电流的8倍。由此可见，指标并不能完全的取代电磁仿真分析。该指标的作用主要是对谐波过电压风险进行定性的描述，利用简单的方法寻找出系统中过电压风险最高的母线，进而与电磁暂态仿真方法相配合，为后续的仿真分析提供参考和理论依据。 
    本发明已经通过上述实施例进行了说明，但应当理解的是，上述实施例只是用于举例和说明的目的，而非意在将本发明限制于所描述的实施例范围内。此外，本领域技术人员可以理解的是，本发明并不局限于上述实施例，根据本发明的教导还可以做出更多种的变型和修改，这些变型和修改均落在本发明所要求保护的范围以内。 

Claims (4)

1.一种评估励磁涌流引发电压严重畸变风险的方法，其特征在于：包括以下步骤：
    （a）生成待求系统各次的谐波导纳阵；
    （b）求取电力系统各次谐波阻抗阵；
    （c）求出重点考察变压器的容量标幺值：选取系统中容量较大、电压等级较高的变压器，将其作为重点考察的谐波注入源，查阅技术资料，获取这些重点考察变压器的容量                                                
Figure 257031DEST_PATH_IMAGE002
，并利用下式计算出容量的标幺值：
                                             
Figure 49537DEST_PATH_IMAGE004
    （d）指标的求取
求出阻抗阵中的谐波互阻抗，计算出三类电压畸变风险指标：
3、5次谐波畸变率风险指标：
Figure 206106DEST_PATH_IMAGE006
电压总畸变率风险指标：
Figure 476681DEST_PATH_IMAGE008
电压瞬时值风险指标：
Figure 815956DEST_PATH_IMAGE010
（e）根据指标判定风险
当系统中某条母线的
Figure 15249DEST_PATH_IMAGE012
指标计算结果大于1时，说明该电网在某些恶劣的工况下，可能存在着对应的风险，
Figure DEST_PATH_IMAGE013
值越大，则风险的程度就越大。
2.根据权利要求1所述的一种评估励磁涌流引发电压严重畸变风险的方法，其特征在于：所述步骤（a）的具体过程为：
（a1）将输电线路等效为π型等值电路，得出谐波阻抗：，式中，是n次谐波下线路的特征阻抗，
Figure DEST_PATH_IMAGE019
是n次谐波下线路的传播系数；
（a2）求出变压器的谐波等值阻抗，
Figure DEST_PATH_IMAGE021
，式中，
Figure DEST_PATH_IMAGE023
、
Figure DEST_PATH_IMAGE025
为变压器n次谐波的谐波电阻与电抗；
Figure DEST_PATH_IMAGE027
、
Figure DEST_PATH_IMAGE029
为变压器基次谐波的谐波电阻与电抗；
（a3）求出发动机的谐波等效阻抗：
Figure DEST_PATH_IMAGE031
，式中，
Figure DEST_PATH_IMAGE033
为电动机在n次谐波下的电抗，
Figure DEST_PATH_IMAGE035
为基波正序阻抗，n为谐波次数；
（a4）求出负荷的谐波阻抗：将负荷等值为一个综合等值电动机处理，其谐波阻抗
Figure DEST_PATH_IMAGE037
为
                      
式中，
Figure DEST_PATH_IMAGE041
、
Figure DEST_PATH_IMAGE043
为等值电动机的基波负序电阻与电抗；
（a5）求出无功补偿元件的谐波阻抗：将无功补偿元件近似的等值为单一的电容、电感元件，因此感性补偿元件与容性补偿元件的谐波阻抗
Figure DEST_PATH_IMAGE045
、
Figure DEST_PATH_IMAGE047
分别为：
Figure DEST_PATH_IMAGE049
                         
式中，
Figure DEST_PATH_IMAGE051
为基波角频率，、
Figure DEST_PATH_IMAGE055
分别为感性补偿元件与容性补偿元件的电感电容值。
3.根据权利要求1所述的一种评估励磁涌流引发电压严重畸变风险的方法，其特征在于：设定励磁涌流的峰值为变压器额定电流的8.0倍，同时选取电流中各次谐波电流的有效值与励磁涌流峰值之间的比例为典型值。
4.根据权利要求1所述的一种评估励磁涌流引发电压严重畸变风险的方法，其特征在于：所述步骤（d）中，由于励磁涌流引入的谐波属于暂态谐波，其持续时间最多为分钟级，因此上述风险的判断依据分别为：1）系统中3、5次谐波电压含有率不得超过20%；2）电网中短时谐波畸变率允许值不得超过10%；3）设定每次谐波的峰值均在同一时刻出现，瞬时值最大允许值不得超过额定峰值的1.5倍。
CN201310660744.1A 2013-12-10 2013-12-10 一种评估励磁涌流引发电压严重畸变风险的方法 Active CN103715684B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310660744.1A CN103715684B (zh) 2013-12-10 2013-12-10 一种评估励磁涌流引发电压严重畸变风险的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310660744.1A CN103715684B (zh) 2013-12-10 2013-12-10 一种评估励磁涌流引发电压严重畸变风险的方法

Publications (2)

Publication Number Publication Date
CN103715684A true CN103715684A (zh) 2014-04-09
CN103715684B CN103715684B (zh) 2016-03-02

Family

ID=50408420

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310660744.1A Active CN103715684B (zh) 2013-12-10 2013-12-10 一种评估励磁涌流引发电压严重畸变风险的方法

Country Status (1)

Country Link
CN (1) CN103715684B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035030A (zh) * 2014-06-09 2014-09-10 中广核工程有限公司 导纳型失磁保护装置的校验方法及系统
CN105093153A (zh) * 2015-08-18 2015-11-25 国网辽宁省电力有限公司大连供电公司 一种基波频率无功功率测量非线性的评价方法
CN107332267A (zh) * 2017-07-25 2017-11-07 中国电力科学研究院 一种避免特高压直流功率持续偏低的控制方法及系统
CN107765076A (zh) * 2016-08-19 2018-03-06 中国电力科学研究院 一种励磁涌流识别方法及识别装置
CN107908906A (zh) * 2017-12-12 2018-04-13 国网四川省电力公司电力科学研究院 一种风电场升压变压器励磁涌流现象仿真方法
CN108199368A (zh) * 2017-12-28 2018-06-22 国网河南省电力公司电力科学研究院 一种确定空载变压器最严重合闸角的方法
CN108268999A (zh) * 2017-12-12 2018-07-10 国网河南省电力公司电力科学研究院 用于电力网络节点谐波特性分析的多主元模态分析法
CN109613324A (zh) * 2018-11-07 2019-04-12 全球能源互联网研究院有限公司 一种谐波放大的检测方法及装置
CN110165670A (zh) * 2019-05-07 2019-08-23 南京师范大学 一种引起换流母线电压畸变的非工频电流快速定位方法
CN112858784A (zh) * 2021-04-03 2021-05-28 国网四川省电力公司电力科学研究院 一种牵引供电系统-区域电网并联谐波谐振频率识别方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6440988A (en) * 1987-08-07 1989-02-13 Toshiba Corp Drive controller for fluorescent indicator tube
CN102790405A (zh) * 2012-07-30 2012-11-21 能科节能技术股份有限公司 一种斩波调压励磁涌流抑制装置及方法
WO2013036855A2 (en) * 2011-09-07 2013-03-14 Quanta Associates, L.P. Real-time monitoring of electric power system voltage stability margins

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6440988A (en) * 1987-08-07 1989-02-13 Toshiba Corp Drive controller for fluorescent indicator tube
WO2013036855A2 (en) * 2011-09-07 2013-03-14 Quanta Associates, L.P. Real-time monitoring of electric power system voltage stability margins
CN102790405A (zh) * 2012-07-30 2012-11-21 能科节能技术股份有限公司 一种斩波调压励磁涌流抑制装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
索南加乐等: "基于波形系数的变压器励磁涌流快速识别算法", 《电网技术》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104035030A (zh) * 2014-06-09 2014-09-10 中广核工程有限公司 导纳型失磁保护装置的校验方法及系统
CN105093153A (zh) * 2015-08-18 2015-11-25 国网辽宁省电力有限公司大连供电公司 一种基波频率无功功率测量非线性的评价方法
CN107765076B (zh) * 2016-08-19 2021-01-15 中国电力科学研究院 一种励磁涌流识别方法及识别装置
CN107765076A (zh) * 2016-08-19 2018-03-06 中国电力科学研究院 一种励磁涌流识别方法及识别装置
CN107332267A (zh) * 2017-07-25 2017-11-07 中国电力科学研究院 一种避免特高压直流功率持续偏低的控制方法及系统
CN107332267B (zh) * 2017-07-25 2021-11-02 中国电力科学研究院 一种避免特高压直流功率持续偏低的控制方法及系统
CN107908906A (zh) * 2017-12-12 2018-04-13 国网四川省电力公司电力科学研究院 一种风电场升压变压器励磁涌流现象仿真方法
CN108268999A (zh) * 2017-12-12 2018-07-10 国网河南省电力公司电力科学研究院 用于电力网络节点谐波特性分析的多主元模态分析法
CN108268999B (zh) * 2017-12-12 2021-05-25 国网河南省电力公司电力科学研究院 用于电力网络节点谐波特性分析的多主元模态分析法
CN108199368A (zh) * 2017-12-28 2018-06-22 国网河南省电力公司电力科学研究院 一种确定空载变压器最严重合闸角的方法
CN108199368B (zh) * 2017-12-28 2021-05-04 国网河南省电力公司电力科学研究院 一种确定空载变压器最严重合闸角的方法
CN109613324A (zh) * 2018-11-07 2019-04-12 全球能源互联网研究院有限公司 一种谐波放大的检测方法及装置
CN109613324B (zh) * 2018-11-07 2021-10-08 全球能源互联网研究院有限公司 一种谐波放大的检测方法及装置
CN110165670A (zh) * 2019-05-07 2019-08-23 南京师范大学 一种引起换流母线电压畸变的非工频电流快速定位方法
CN112858784A (zh) * 2021-04-03 2021-05-28 国网四川省电力公司电力科学研究院 一种牵引供电系统-区域电网并联谐波谐振频率识别方法
CN112858784B (zh) * 2021-04-03 2023-05-26 国网四川省电力公司电力科学研究院 一种牵引供电系统-区域电网并联谐波谐振频率识别方法

Also Published As

Publication number Publication date
CN103715684B (zh) 2016-03-02

Similar Documents

Publication Publication Date Title
CN103715684A (zh) 一种评估励磁涌流引发电压严重畸变风险的方法
CN103944507B (zh) 基于逆变器型式试验的光伏电站低电压穿越性能评价方法
CN103472325B (zh) 一种波浪能和潮流能独立发电系统的检测平台
CN106610456A (zh) 一种风电场无功补偿装置并网性能测试方法
CN103745267A (zh) 一种分布式光伏系统并网影响的评价方法
CN105184419A (zh) 一种配网电压暂降及短时中断治理方案的优选方法
CN105548741A (zh) 一种用于风电场无功补偿装置低电压运行能力的检测方法
CN103050933A (zh) 基于单端电流突变量的大型电池储能电站接口保护方法
CN203243056U (zh) 中高压无功补偿复合结构
CN109870615B (zh) 适用于工矿企业的谐波责任评定方法
CN105305421A (zh) 一种统一电能质量控制器的实验系统与方法
CN102542355A (zh) 一种运行规划分负荷方法
Fan et al. A comparative investigation of flexible arc suppression devices in distribution network
CN201548607U (zh) 变电站谐振过电压分析系统
CN107390014A (zh) 波动性负荷闪变发射水平的测量方法
Zhang et al. Electromagnetic Transient Simulation Research on Operation Characteristics of Power Grid with Large-Scale New Energy
Geng et al. Study on the Harmonic Amplification Influence of DC Inverter Substation With Cable Outgoing Line on the Receiving-end Power Grid
Li et al. Power-quality of distribution networks with high penetrated intermittent distributed generation: A survey
Wang et al. Harmonic line loss suppression and rectification measures
Alghamdi et al. Analysis of harmonic propagations in albaha power network due to the implementation of an mvdc converter
CN108512228A (zh) 一种电气化铁路对沿途风电场电能质量影响快速评估方法
Lei et al. Effect of load on zero-sequence current in low resistance grounding system with complex grounding fault
Kang et al. Power quality problem and key improvement technology for regional power grids
CN104092197B (zh) 基于差动系数矩阵的特高压交流输电线路继电保护方法
CN108667000B (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
C14 Grant of patent or utility model
GR01 Patent grant