CN114611278A - 一种多阶段电压暂降状态估计方法 - Google Patents
一种多阶段电压暂降状态估计方法 Download PDFInfo
- Publication number
- CN114611278A CN114611278A CN202210201077.XA CN202210201077A CN114611278A CN 114611278 A CN114611278 A CN 114611278A CN 202210201077 A CN202210201077 A CN 202210201077A CN 114611278 A CN114611278 A CN 114611278A
- Authority
- CN
- China
- Prior art keywords
- fault
- voltage sag
- sag
- line
- matrix
- 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 43
- 230000004224 protection Effects 0.000 claims abstract description 131
- 239000011159 matrix material Substances 0.000 claims abstract description 99
- 230000009471 action Effects 0.000 claims abstract description 31
- 230000035772 mutation Effects 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000005457 optimization Methods 0.000 claims description 20
- 238000012806 monitoring device Methods 0.000 claims description 10
- 238000010248 power generation Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 4
- 230000001681 protective effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/25—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
- G01R19/2513—Arrangements for monitoring electric power systems, e.g. power lines or loads; Logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16528—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values using digital techniques or performing arithmetic operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16533—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application
- G01R19/16538—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application in AC or DC supplies
- G01R19/16547—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application in AC or DC supplies voltage or current in AC supplies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16566—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533
- G01R19/16576—Circuits and arrangements for comparing voltage or current with one or several thresholds and for indicating the result not covered by subgroups G01R19/16504, G01R19/16528, G01R19/16533 comparing DC or AC voltage with one threshold
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Computing Systems (AREA)
- Strategic Management (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Water Supply & Treatment (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Emergency Protection Circuit Devices (AREA)
Abstract
本发明公开了一种多阶段电压暂降状态估计方法,首先分析不同原因导致的多阶段电压暂降特征,提出利用暂降幅值跳变特征确定多阶段暂降幅值突变时刻成因的方法;然后提出继电保护动作矩阵和故障线路集的计算方法,实现利用继电保护装置的故障清除时间等特征初步划分故障线路集的目的,有效减少暂降状态估计的计算量;接着利用不同事件导致系统阻抗变化这一本质特征,提出多阶段电压暂降各阶段事件成因的推断方法,提高状态估计的准确性;最后基于事件成因推断结果,提出了多阶段电压暂降状态估计方法,解决了现有状态估计方法难以适用于多阶段暂降的问题。
Description
技术领域
本发明涉及电压暂降状态估计技术领域,具体为一种多阶段电压暂降状态估计方法。
背景技术
现有电压暂降状态估计方法主要针对由短路引起的单阶段矩形暂降,缺乏适用于多阶段电压暂降的状态估计方法。然而现有研究表明,多阶段电压暂降事件的发生频次逐渐上升,多阶段暂降问题不容忽视。例如:随着新型电力系统的发展,系统中接入的DG比例也越来越高。在短路故障发生期间,暂降可能导致单个DG跳闸或多个DG连锁跳闸,从而形成复杂的多阶段暂降。此外,电力系统中常用的阶段式保护,比如距离保护也可能因为其动作特性导致多阶段暂降。
电压暂降状态估计需要利用短路计算公式来估计未监测母线的暂降幅值。短路计算方法为公知性方法,不同故障类型下短路计算方法如下。
设m为待求母线,母线i、j之间的线路lij故障,f为故障点,p为点f到母线i的归一化距离。p取值为[0,1]表示f从i到j,如下式所示:
式中,Lif为母线i到故障点f的距离;Lij为线路lij的长度。
①三相短路故障
母线m发生三相故障导致的三相电压幅值为:
②单相接地故障
设A相发生单相接地故障,母线m各相电压幅值为:
式中,α=ej120°为旋转因子。
③两相故障
设B、C相发生两相故障,母线m各相电压幅值为:
④两相接地故障
设B、C相发生两相接地故障,母线m各相电压幅值为:
暂降是对敏感设备影响最严重的扰动,且无法避免。要解决电压暂降问题,首先应该能够准确反映系统中不同位置的电压暂降水平,因此需要安装必要的监测装置用于监测暂降。由于电网庞大而复杂,出于监测成本和信息处理能力考虑,系统中通常只能在部分母线安装监测装置。如何利用有限监测装置的监测数据估计未监测母线的暂降水平(主要指电压暂降幅值与持续时间)是一个重要的研究任务,这一研究又称之为电压暂降状态估计。现有电压暂降状态估计方法,主要针对短路故障引起的单阶段暂降,难以适用于多阶段电压暂降的状态估计。随着多阶段电压暂降发生频率的逐渐增加,提出适用于多阶段电压暂降的状态估计方法具有重要意义。
多阶段状态估计的难点主要有如下几个。1.形成多阶段暂降的原因主要是短路期间DG脱网或者线路两端的保护装置不同时跳闸,如何判断不同阶段的暂降是何种原因导致是第一个挑战。2.短路发生的位置是影响电压幅值的重要因素,如何初步划分可能的故障线路集,减小MVSSE的计算量是第二个挑战。3.短路期间,同一故障位置处的电压幅值还受DG的脱网个数、顺序以及故障线路上保护装置的跳闸情况等事件影响,如何推断出暂降事件不同阶段的具体成因,提高MVSSE的准确性是第三个挑战。
术语解释:
电压暂降:国际电气与电子工程师协会(IEEE)将电压暂降(voltage sag)定义为供电电压有效值快速下降到0.1~0.9p.u.,持续时间为0.5个周波至1min的电能质量现象。
电压暂降状态估计:一种通过有限监测装置的监测数据估计非监测母线电压暂降状态(指一定时间内的暂降频次、母线的电压暂降幅值或暂降持续时间)的方法。本发明中的电压暂降状态指单一暂降事件中的电压暂降幅值与电压暂降持续时间。
电压暂降突变点:电压暂降有效值波形中电压幅值突然大幅变化(超过正常电压波动范围)的点。短路故障、短路后DG(Distributed Generation分布式发电装置)脱网、短路后故障线路单侧继电保护装置动作都会形成电压暂降突变点。电压暂降突变点可以通过波形点检测方法确定。
发明内容
针对上述问题,本发明的目的在于提供一种多阶段电压暂降状态估计方法,通过分析不同原因造成多阶段暂降的特点,考虑多阶段暂降监测数据中的各阶段突变时间与继电保护装置的故障清除时间,推断出多阶段暂降各个阶段可能发生的具体事件,最终得到符合工程需求的多阶段状态估计方法。技术方案如下:
一种多阶段电压暂降状态估计方法,包括以下步骤:
步骤1:分析多阶段电压暂降的成因,确定不同成因下多阶段电压暂降特点,据此对不同突变点处电压暂降成因进行识别;
步骤2:根据电网不同位置处整定的故障清除时间,构建继电保护装置动作矩阵;
步骤3:求解出一条或多条可能的故障线路,构成故障线路集,从而减少多阶段电压暂降状态估计的计算量;
步骤4:计算故障前后系统阻抗矩阵,以及故障线路单侧继电保护动作时的阻抗矩阵和分布式发电装置脱网时的阻抗矩阵;
步骤5:利用优化模型从步骤3中求解得到的故障线路集中推断具体的故障线路,基于电压暂降类型推断出多阶段暂降各阶段具体事件;
步骤6:基于电压暂降类型,通过各相邻电压暂降突变点之间的时间长度对暂降各阶段的持续时间进行估计,基于步骤5中的事件推断结果对暂降幅值进行估计。
进一步的,所述步骤1中,多阶段电压暂降的成因和特点为:
成因Ⅰ:短路故障期间,电力系统中的一个或多个分布式发电装置跳闸,导致电网丢失电源;其特点为:电压暂降发生时刻对应的电压暂降突变点右侧电压幅值小于左侧;
成因Ⅱ:短路故障期间,故障线路两侧的继电保护装置不同时跳闸,导致电网拓扑发生变化;其特点为:电压暂降发生时刻对应的电压暂降突变点右侧电压幅值大于左侧。
更进一步的,所述对不同突变点处电压暂降成因进行识别的具体方法为:
假设对一段包含电压暂降的波形使用波形点检测方法一共检测出s+1个电压暂降突变点,分别表示为MP0,MP1,...,MPx,…,MPs,且对应s+1个突变时刻t0,t1,…,tx,...,ts,这些电压暂降突变点将电压暂降部分的波形分为s段;不同段波形的电压分别表示为u1,u2,...,ux,...,us,且第m个电压um计算如下:
其中,ff表示采样率,单位为个/秒;f表示工频电流频率,单位为Hz;xj是有效值序列中的第j个值;m∈[1,s];k为……;
根据上式形成二进制向量W:
其中,向量W中的元素wx有0和1两种取值;当wx取值为0时表示tx-1时刻对应的MPx-1的右侧电压小于左侧,这意味着该时刻因为成因Ⅰ而形成了多阶段电压暂降,其中,x>1;反之,当wx取值为1时,意味着tx-1时刻因为成因Ⅱ而形成了多阶段电压暂降,其中,x>1。
更进一步的,所述步骤2中构建继电保护装置动作矩阵具体包括:
步骤2.1:构建基本动作矩阵,即描述继电保护装置排除故障行为的通用矩阵:
其中,n表示系统中的母线条数,γij和γji分别表示线路lij上靠近i和j侧母线的保护的故障清除时间,i,j∈[1,n];i≠j;γij,γji=0表示母线i和j之间物理上没有连接;
步骤2.2:将所述通用矩阵改进并解耦为两个矩阵:
上式中的上下三角矩阵分别表示线路两个端子母线上同一类型保护装置的故障清除时间;下三角矩阵中的下三角元素代表线路lij的i侧母线的保护装置故障清除时间,且γij=0,i>j,表示母线i和j之间物理上没有连接;上三角矩阵中的上三角元素代表线路另一侧的参数;
为表示二段式保护,构造另外两个类似矩阵表示后备保护参数:
步骤2.3:确定保护装置的配合关系
对于二段式保护,线路中的主保护和后备保护的配合关系包括四种,将改进的保护动作矩阵表示为以下几种配合关系:
其中,ΘⅠ表示故障线路两侧的主保护相互配合切除故障;ΘⅡ表示故障线路两侧的后备保护相互配合切除故障;ΘⅢ和ΘⅣ表示故障线路某侧的主保护与另一侧的后备保护相互配合切除故障;
步骤2.4:修正动作矩阵中的故障清除时间
将故障清除时间修正如下:
其中,γij,set和λij,set分别表示主保护、后备保护的整定值;δ1,ij和δ2,ij分别表示实际的故障清除时间与整定值之间的偏差;偏差是[0,δ]内的随机数,其中δ表示同型号继电保护装置试验或历史运行过程中的最大误差值。
更进一步的,所述步骤3具体包括:
步骤3.1:基于步骤1的暂降原因识别,将电压暂降分为如下四种类型:
类型Ⅰ:单阶段矩形暂降:
类型Ⅱ:由成因Ⅰ导致的多阶段电压暂降;
类型Ⅲ:由成因Ⅱ导致的多阶段电压暂降;
类型Ⅳ:由成因Ⅰ和成因Ⅱ共同导致的多阶段电压暂降;
步骤3.2:不同类型暂降下故障线路集求解
类型Ⅰ和类型Ⅱ:
假设步骤1中检测出电压暂降波形的第一个和最后一个电压暂降突变点时刻分别是t0和ts,则系统中故障发生到被切除的时间长度为ts-t0;若步骤2.3四个矩阵中某个矩阵主对角线对称位置的两个元素在误差阈值内等于ts-t0,则与这两个元素相对应的线路为可能的故障线路;因此故障线路集LF的求解模型如下:
其中,θij是矩阵ΘⅠ~ΘⅣ中第i行中第j列的元素,LN表示监测装置母线处暂降域交集中的线路组成的集合;
类型Ⅲ和类型Ⅳ:
假设步骤1中检测出电压暂降波形的第一个电压暂降突变点为t0,两个对应继电保护动作的电压暂降突变点时刻分别是tx-1和ts;则故障线路两侧继电保护装置从故障发生到保护动作的时间长度分别为tx-1-t0和ts-t0;那么故障线路集的求解模型如下:
更进一步的,所述步骤4具体包括:
步骤4.1:计算故障前系统阻抗矩阵
式中,se=1,2,0表示正序、负序和零序,αij表示节点i、j的互导纳,αii表示节点i的自导纳,i≠j;i,j∈[1,n];
式中,βii=0表示母线处没有发电机;
则计算系统的阻抗矩阵如下:
步骤4.2:计算短路之后系统阻抗矩阵
式中,和是矩阵中的元素,是线路lij阻抗,进一步将阻抗表示为函数的形式gmf1(i,j,p)和gff1(i,j,p),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;
步骤4.3:计算故障线路单侧继电保护动作时的阻抗矩阵
假设线路lij发生短路故障后,j侧保护装置动作将单侧线路切除,系统的阻抗计算如下:
其中,ΔZ表示过程量,使用下式计算:
与相比增加的一行一列代表故障位置fl与各目标母线m的互阻抗或故障位置的自阻抗进一步将阻抗表示为函数的形式gmf2(i,j,p,d)和gff2(i,j,p,d),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障后,由线路d侧保护装置动作切除部分线路,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;如下两式所示:
步骤4.4:计算分布式发电装置脱网时的阻抗矩阵
假设母线集h中所有母线的分布式发电装置脱网,此时首先需要将对角线矩阵修正为下式:
根据上式计算故障前系统的阻抗矩阵,最后计算短路情况下分布式发电装置脱网时的系统阻抗,并进一步将阻抗表示为函数的形式gmf3(i,j,p,h)和gff3(i,j,p,h),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障且母线集h内所有母线处的DG脱网时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;如下两式所示:
更进一步的,所述步骤5中优化模型根据暂降类型分为如下四种:
(1)类型Ⅰ使用如下优化模型推断故障线路及其短路情况,即i,j,p;
其中,f(·)为短路计算的函数;
(2)类型Ⅱ使用如下优化模型推断故障线路及其短路情况、DG脱网时刻及脱网顺序,即i,j,p,hq;
其中,s为多阶段电压暂降的阶段数;hq为第q阶段电压暂降期间,脱网DG组成的集合;
(3)类型Ⅲ使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序,即i,j,p,d;
(4)类型Ⅳ使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序、DG脱网时刻及其顺序,即i,j,p,d,hq;
其中,o为故障线路某侧继电保护装置首先动作跳闸的时刻;uo+1为第o+1段电压暂降幅值。
更进一步的,所述步骤6中电压暂降的状态估计根据暂降类型分为如下四种:
1)对于类型Ⅰ,任何未监测母线m的暂降幅值根据下式进行估算:
f(gmf1(i,j,p),gff1(i,j,p))
2)对于类型Ⅱ,任何未监测母线m的暂降幅值根据下式进行估算:
3)对于类型Ⅲ,任何未监测母线m的暂降幅值根据下式进行估算:
4)对于类型Ⅳ,任何未监测母线m的暂降幅值根据下式进行估算。
其中,y和z分别为多阶段电压暂降的第y和z阶段。
本发明的有益效果是:
1)本发明分析不同原因导致的多阶段电压暂降特征,提出利用暂降幅值跳变特征确定多阶段暂降幅值突变时刻成因的方法。
2)本发明提出继电保护动作矩阵和故障线路集的计算方法,实现利用继电保护装置的故障清除时间等特征初步划分故障线路集的目的,有效减少暂降状态估计的计算量。
3)本发明利用不同事件导致系统阻抗变化这一本质特征,提出多阶段电压暂降各阶段事件成因的推断方法,提高状态估计的准确性。
4)本发明基于事件成因推断结果,提出了多阶段电压暂降状态估计方法,解决了现有状态估计方法难以适用于多阶段暂降的问题。
附图说明
图1为电力系统的短路计算模型。
图2为本发明故障线路被单侧切除后的电力系统短路计算模型。
图3为本发明多阶段电压暂降状态估计模型构建流程。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。
本发明技术方案主要分为六个大步骤:多阶段电压暂降成因识别、继电保护动作矩阵构建、故障线路集求解、阻抗矩阵求解、推断多阶段暂降各阶段具体事件和多阶段电压暂降状态估计方法。其中每个大步骤及其小步骤的详细阐述如下:
一、多阶段电压暂降成因识别
①多阶段电压暂降的成因有两个:
成因1是:短路故障期间,电力系统中的一个或多个DG跳闸,导致电网丢失电源。成因2是:短路故障期间,故障线路两侧的继电保护装置不同时跳闸,导致电网拓扑发生变化。
②不同成因下多阶段电压暂降特点
成因1意味着故障之后,系统进一步丢失电源,因此监测到的电压暂降幅值进一步降低。其特点是,成因1发生时刻对应的电压暂降突变点右侧电压幅值小于左侧。
相反,成因2意味着故障发生之后,线路的单侧保护装置将故障与系统部分隔离,因此监测到的电压暂降幅值会有所上升。其特点是,成因2发生时刻对应的电压暂降突变点右侧电压幅值大于左侧。
③不同突变点处电压暂降成因识别
假设对一段包含电压暂降的波形使用波形点检测方法一共检测出s+1个电压暂降突变点,分别表示为MP0,MP1,…,MPx,…,MPs,且对应s+1个突变时刻t0,t1,…,tx,…,ts,这些电压暂降突变点将电压暂降部分的波形分为s段。不同段波形的电压可以分别表示为u1,u2,…,ux,…,us,且第m(m∈[1,s])个电压um可以计算如下:
其中ff表示采样率,单位为个/秒,f表示工频电流频率,单位为Hz,xj是有效值序列中的第j个值。根据式(7)可以形成二进制向量W如式(8):
向量W中的元素wx有0和1两种取值。当取值为0时表示tx-1(x>1)时刻对应的MPx-1的右侧电压小于左侧,这意味着该时刻因为成因1而形成了多阶段电压暂降。反之,当取值为1时,意味着tx-1(x>1)时刻因为成因2而形成了多阶段电压暂降.
二、继电保护装置动作矩阵
基于不同电压暂降突变点对应时刻的多阶段暂降成因,可以大致获取暂降事件中继电保护装置的实际跳闸时间。为了利用该信息,第二大步骤根据电网不同位置处整定的故障清除时间提出继电保护装置动作矩阵的构建方法。
①基本动作矩阵
继电保护有不同的类型,包括低压保护、过电流保护和距离保护等,每种类型都有其独有的触发规则与故障清除时间来清除故障。本发明提出了一种描述继电保护装置排除故障行为的通用矩阵,如下式所示。
其中n表示系统中的母线条数,γij和γji(i,j∈[1,n];i≠j)分别表示线路lij上靠近i和j侧母线的保护的故障清除时间,单位是ms。γij,γji=0表示母线i和j之间物理上没有连接。故障清除时间是配电自动化系统等电力企业的软件中记录的继电保护装置整定值。
②改进动作矩阵
在实际应用过程中,阶段式保护广泛用于配电系统,比如二段式保护。二段式保护是指在一条线路上配备两种类型的保护,一种作为主保护,另一种作为后备保护。两种保护的故障清除时间不同,主保护的故障清除时间更短。当发生短路故障时,主保护优先触发,后备保护则在主保护无法触发的情况下工作。一般情况下故障可以由主保护独立切除,有时也可以由后备保护独立切除或者由主保护与后备保护配合切除。本发明考虑到式(9)中定义的基本动作矩阵难以呈现阶段式保护的配合关系,因此对其改进并解耦为两个矩阵,分别如式(10)和(11)所示。
上式中的上下三角矩阵分别表示线路两个端子母线上同一类型保护装置的故障清除时间。下三角矩阵中的下三角元素代表线路lij的i侧母线的保护装置故障清除时间,且γij=0(i>j)表示母线i和j之间物理上没有连接。上三角矩阵中的上三角元素代表线路另一侧的参数,含义类似。设式(10)和(11)代表主保护参数,为了表示二段式保护,构造另外两个类似矩阵表示后备保护参数,分别如式(12)和(13)所示。
③保护装置的配合关系
对于广泛使用的二段式保护,线路中的主保护和后备保护可能有四种配合关系。改进的保护动作矩阵可以方便的表示这几种配合关系如下:
关系1:如式(14)所示,ΘⅠ表示故障线路两侧的主保护相互配合切除故障。
关系2:如式(15)所示,ΘⅡ表示故障线路两侧的后备保护相互配合切除故障。
关系3/4:如式(16)和(17)所示,ΘⅢ和ΘⅣ表示故障线路某侧的主保护与另一侧的后备保护相互配合切除故障。
④修正动作矩阵中的故障清除时间
式(9)~(17)中使用到的故障清除时间是从电力公司获取到的整定值。理论上线路两端的主保护或后备保护各自的故障清除时间应该相同,且等于整定值。然而实际运行中,由于误差的存在,这些矩阵中的元素可能偏离整定值,因此需要做出必要修正如式(18)所示。
其中,γij,set和λij,set分别表示主保护、后备保护的整定值。δ1,ij和δ2,ij分别表示实际的故障清除时间与整定值之间的偏差。偏差是[0,δ]内的随机数,其中δ表示同型号继电保护装置试验或历史运行过程中的最大误差值。
三、故障线路集求解
基于第一步暂降原因识别和第二步保护动作矩阵的内容,可以求解出一条或多条可能的故障线路,从而减少多阶段电压暂降状态估计的计算量。
①暂降类型
基于第一步暂降原因识别的内容,将电压暂降分为如下四种类型:
类型1:单阶段矩形暂降。
类型2:由原因1导致的多阶段电压暂降。
类型3:由原因2导致的多阶段电压暂降。
类型4:由原因1和2共同导致的多阶段电压暂降。
②不同类型暂降下故障线路集求解
类型1/2:这两类电压暂降事件中,故障线路两侧继电保护装置在最后一个电压暂降突变时刻同时动作将故障切除。假设步骤一中检测出电压暂降波形的第一个和最后一个电压暂降突变点时刻分别是t0和ts,则系统中故障发生到被切除的时间长度为:ts-t0。如果式(14)~(17)中某个矩阵主对角线对称位置的两个元素在误差阈值内等于ts-t0,则与这两个元素相对应的线路就是可能的故障线路。因此故障线路集LF的求解模型如下:
其中,θij是矩阵ΘⅠ~ΘⅣ中第i行中第j列的元素,LN表示监测装置母线处暂降域交集中的线路组成的集合。
类型3/4:这两类电压暂降事件中,故障线路两侧继电保护装置不能同时动作将故障切除。假设步骤一中检测出电压暂降波形的第一个电压暂降突变点及两个对应继电保护动作的电压暂降突变点时刻分别是t0、tx-1和ts。则故障线路两侧继电保护装置从故障发生到保护动作的时间长度分别为tx-1-t0和ts-t0。那么故障线路集的求解模型如下:
四、阻抗矩阵求解
无论是成因1还是成因2造成的多阶段暂降,其本质都是故障过程中由事件引起的系统阻抗矩阵的变化。
①故障前系统阻抗矩阵
式中,se=1,2,0表示正序、负序和零序,αij(i≠j;i,j∈[1,n])表示节点i、j的互导纳,αii(i∈[1,n])表示节点i的自导纳。
式中βii=0(i∈[1,n])表示母线处没有发电机。
利用式(21)可以计算系统的阻抗矩阵如下。
②短路之后系统阻抗矩阵
式中和是矩阵中的元素,是线路lij阻抗,进一步将阻抗表示为函数的形式gmf1(i,j,p)和gff1(i,j,p),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗。
③故障线路单侧继电保护动作时的阻抗矩阵
假设线路lij发生短路故障后,j侧保护装置动作将单侧线路切除,系统的阻抗可以通过图2所示模型计算。
其中ΔZ可以使用下式计算:
与相比增加的一行一列代表故障位置f与各目标母线m的互阻抗或故障位置的自阻抗进一步将阻抗表示为函数的形式gmf2(i,j,p,d)和gff2(i,j,p,d),两个函数分别表示线路lij上距离i端p的位置f发生短路故障后,由线路d侧保护装置动作切除部分线路,故障位置与目标节点m的互阻抗以及故障位置的自阻抗。如下两式所示:
④DG脱网时的阻抗矩阵
当短路发生时,某些节点处的DG可能因为低压穿越能力不够而脱网。假设母线集h中所有母线的DG脱网,此时首先需要将式(23)所示矩阵修正为下式。
然后将上式带入式(24)计算故障前系统的阻抗矩阵,最后利用式(1)~(26)计算短路情况下DG脱网时的系统阻抗,并进一步将阻抗表示为函数的形式gmf3(i,j,p,h)和gff3(i,j,p,h),两个函数分别表示线路lij上距离i端p的位置f发生短路故障且母线集h内所有母线处的DG脱网时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗。如下两式所示:
五、推断多阶段暂降各阶段具体事件
第四步求解出不同情况下系统阻抗矩阵之后,可以利用第五步中的优化模型从第三步中求解得到的故障线路集中推断具体的故障线路,并推断出多阶段暂降各阶段具体事件。推断模型同样根据暂降类型分为如下四种:
①类型1
对于这一类暂降,可以使用如下优化模型推断故障线路及其短路情况(即i,j,p)
其中f()为短路计算的函数,不同短路类型的计算公式见式(2)~(6)。
②类型2
对于这一类暂降,可以使用如下优化模型推断故障线路及其短路情况、DG脱网时刻及脱网顺序(即i,j,p,hq)
③类型3
对于这一类暂降,可以使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序(即i,j,p,d)
④类型4
假设一段电压暂降波形中检测到的第o个电压暂降突变点对应于单侧保护跳闸。对于这一类暂降,可以使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序、DG脱网时刻及其顺序(即i,j,p,d,hq)
需要注意的是,式(35)~(38)是单台监测装置的优化模型,这些优化模型可以通过智能算法比如量子遗传算法求解。当有多个监测装置时,优化模型中的目标函数是监测装置各自目标函数的总和。
六、多阶段电压暂降的状态估计方法
电压暂降的状态估计包括未监测母线的暂降幅值和持续时间的估计。暂降各阶段的持续时间可以通过各相邻电压暂降突变点之间的时间长度进行估计。暂降幅值的估计则基于第五步中的事件推断结果。同样,根据电压暂降类型可以分为如下四种。
①类型1
对于这类暂降,任何未监测母线m的暂降幅值可以根据下式进行估算。
f(gmf1(i,j,p),gff1(i,j,p)) (39)
②类型2
对于这类暂降,任何未监测母线m的暂降幅值可以根据下式进行估算。
③类型3
对于这类暂降,任何未监测母线m的暂降幅值可以根据下式进行估算。
④类型4
对于这类暂降,任何未监测母线m的暂降幅值可以根据下式进行估算。
本发明的方法分析了不同原因造成多阶段暂降的特点,考虑了多阶段暂降监测数据中的各阶段突变时间与继电保护装置的故障清除时间,推断出多阶段暂降各个阶段可能发生的具体事件,最终得到符合工程需求的多阶段状态估计方法。
Claims (8)
1.一种多阶段电压暂降状态估计方法,其特征在于,包括以下步骤:
步骤1:分析多阶段电压暂降的成因,确定不同成因下多阶段电压暂降特点,据此对不同突变点处电压暂降成因进行识别;
步骤2:根据电网不同位置处整定的故障清除时间,构建继电保护装置动作矩阵;
步骤3:求解出一条或多条可能的故障线路,构成故障线路集,从而减少多阶段电压暂降状态估计的计算量;
步骤4:计算故障前后系统阻抗矩阵,以及故障线路单侧继电保护动作时的阻抗矩阵和分布式发电装置脱网时的阻抗矩阵;
步骤5:利用优化模型从步骤3中求解得到的故障线路集中推断具体的故障线路,基于电压暂降类型推断出多阶段暂降各阶段具体事件;
步骤6:基于电压暂降类型,通过各相邻电压暂降突变点之间的时间长度对暂降各阶段的持续时间进行估计,基于步骤5中的事件推断结果对暂降幅值进行估计。
2.根据权利要求1所述的多阶段电压暂降状态估计方法,其特征在于,所述步骤1中,多阶段电压暂降的成因和特点为:
成因Ⅰ:短路故障期间,电力系统中的一个或多个分布式发电装置跳闸,导致电网丢失电源;其特点为:电压暂降发生时刻对应的电压暂降突变点右侧电压幅值小于左侧;
成因Ⅱ:短路故障期间,故障线路两侧的继电保护装置不同时跳闸,导致电网拓扑发生变化;其特点为:电压暂降发生时刻对应的电压暂降突变点右侧电压幅值大于左侧。
3.根据权利要求2所述的多阶段电压暂降状态估计方法,其特征在于,所述对不同突变点处电压暂降成因进行识别的具体方法为:
假设对一段包含电压暂降的波形使用波形点检测方法一共检测出s+1个电压暂降突变点,分别表示为MP0,MP1,...,MPx,...,MPs,且对应s+1个突变时刻t0,t1,...,tx,...,ts,这些电压暂降突变点将电压暂降部分的波形分为s段;不同段波形的电压分别表示为u1,u2,...,ux,...,us,且第m个电压um计算如下:
其中,ff表示采样率,单位为个/秒;f表示工频电流频率,单位为Hz;xj是有效值序列中的第j个值;m∈[1,s];k为采样点总个数;
根据上式形成二进制向量W:
其中,向量W中的元素wx有0和1两种取值;当wx取值为0时表示tx-1时刻对应的MPx-1的右侧电压小于左侧,这意味着该时刻因为成因Ⅰ而形成了多阶段电压暂降,其中,x>1;反之,当wx取值为1时,意味着tx-1时刻因为成因Ⅱ而形成了多阶段电压暂降,其中,x>1。
4.根据权利要求3所述的多阶段电压暂降状态估计方法,其特征在于,所述步骤2中构建继电保护装置动作矩阵具体包括:
步骤2.1:构建基本动作矩阵,即描述继电保护装置排除故障行为的通用矩阵:
其中,n表示系统中的母线条数,γij和γji分别表示线路lij上靠近i和j侧母线的保护的故障清除时间,i,j∈[1,n];i≠j;γij,γji=0表示母线i和j之间物理上没有连接;
步骤2.2:将所述通用矩阵改进并解耦为两个矩阵:
上式中的上下三角矩阵分别表示线路两个端子母线上同一类型保护装置的故障清除时间;下三角矩阵中的下三角元素代表线路lij的i侧母线的保护装置故障清除时间,且γij=0,i>j,表示母线i和j之间物理上没有连接;上三角矩阵中的上三角元素代表线路另一侧的参数;
为表示二段式保护,构造另外两个类似矩阵表示后备保护参数:
步骤2.3:确定保护装置的配合关系
对于二段式保护,线路中的主保护和后备保护的配合关系包括四种,将改进的保护动作矩阵表示为以下几种配合关系:
其中,ΘⅠ表示故障线路两侧的主保护相互配合切除故障;ΘⅡ表示故障线路两侧的后备保护相互配合切除故障;ΘⅢ和ΘⅣ表示故障线路某侧的主保护与另一侧的后备保护相互配合切除故障;
步骤2.4:修正动作矩阵中的故障清除时间
将故障清除时间修正如下:
其中,γij,set和λij,set分别表示主保护、后备保护的整定值;δ1,ij和δ2,ij分别表示实际的故障清除时间与整定值之间的偏差;偏差是[0,δ]内的随机数,其中δ表示同型号继电保护装置试验或历史运行过程中的最大误差值。
5.根据权利要求4所述的多阶段电压暂降状态估计方法,其特征在于,所述步骤3具体包括:
步骤3.1:基于步骤1的暂降原因识别,将电压暂降分为如下四种类型:
类型Ⅰ:单阶段矩形暂降:
类型Ⅱ:由成因Ⅰ导致的多阶段电压暂降;
类型Ⅲ:由成因Ⅱ导致的多阶段电压暂降;
类型Ⅳ:由成因Ⅰ和成因Ⅱ共同导致的多阶段电压暂降;
步骤3.2:不同类型暂降下故障线路集求解
类型Ⅰ和类型Ⅱ:
假设步骤1中检测出电压暂降波形的第一个和最后一个电压暂降突变点时刻分别是t0和ts,则系统中故障发生到被切除的时间长度为ts-t0;若步骤2.3四个矩阵中某个矩阵主对角线对称位置的两个元素在误差阈值内等于ts-t0,则与这两个元素相对应的线路为可能的故障线路;因此故障线路集LF的求解模型如下:
其中,θij是矩阵ΘⅠ~ΘⅣ中第i行中第j列的元素,LN表示监测装置母线处暂降域交集中的线路组成的集合;
类型Ⅲ和类型Ⅳ:
假设步骤1中检测出电压暂降波形的第一个电压暂降突变点为t0,两个对应继电保护动作的电压暂降突变点时刻分别是tx-1和ts;则故障线路两侧继电保护装置从故障发生到保护动作的时间长度分别为tx-1-t0和ts-t0;那么故障线路集的求解模型如下:
6.根据权利要求5所述的多阶段电压暂降状态估计方法,其特征在于,所述步骤4具体包括:
步骤4.1:计算故障前系统阻抗矩阵
式中,se=1,2,0表示正序、负序和零序,αij表示节点i、j的互导纳,αii表示节点i的自导纳,i≠j;i,j∈[1,n];
式中,βii=0表示母线处没有发电机;
则计算系统的阻抗矩阵如下:
步骤4.2:计算短路之后系统阻抗矩阵
式中,和是矩阵中的元素,是线路lij阻抗,进一步将阻抗表示为函数的形式gmf1(i,j,p)和gff1(i,j,p),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;
步骤4.3:计算故障线路单侧继电保护动作时的阻抗矩阵
其中,ΔZ为过程量,使用下式计算:
与相比增加的一行一列代表故障位置fl与各目标母线m的互阻抗或故障位置的自阻抗进一步将阻抗表示为函数的形式gmf2(i,j,p,d)和gff2(i,j,p,d),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障后,由线路d侧保护装置动作切除部分线路,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;如下两式所示:
步骤4.4:计算分布式发电装置脱网时的阻抗矩阵
假设母线集h中所有母线的分布式发电装置脱网,此时首先需要将对角线矩阵修正为下式:
根据上式计算故障前系统的阻抗矩阵,最后计算短路情况下分布式发电装置脱网时的系统阻抗,并进一步将阻抗表示为函数的形式gmf3(i,j,p,h)和gff3(i,j,p,h),两个函数分别表示线路lij上距离i端p的位置fl发生短路故障且母线集h内所有母线处的DG脱网时,故障位置与目标节点m的互阻抗以及故障位置的自阻抗;如下两式所示:
7.根据权利要求6所述的多阶段电压暂降状态估计方法,其特征在于,所述步骤5中优化模型根据暂降类型分为如下四种:
(1)类型Ⅰ使用如下优化模型推断故障线路及其短路情况,即i,j,p;
其中,f(·)为短路计算的函数;
(2)类型Ⅱ使用如下优化模型推断故障线路及其短路情况、DG脱网时刻及脱网顺序,即i,j,p,hq;
其中,s为多阶段电压暂降的阶段数;hq为第q阶段电压暂降期间,脱网DG组成的集合;
(3)类型Ⅲ使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序,即i,j,p,d;
(4)类型Ⅳ使用如下优化模型推断故障线路及其短路情况、继电保护装置跳闸时刻及其跳闸顺序、DG脱网时刻及其顺序,即i,j,p,d,hq;
其中,o为故障线路某侧继电保护装置首先动作跳闸的时刻;uo+1为第o+1段电压暂降幅值。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210201077.XA CN114611278B (zh) | 2022-03-03 | 2022-03-03 | 一种多阶段电压暂降状态估计方法 |
US18/082,611 US20230305050A1 (en) | 2022-03-03 | 2022-12-16 | State estimation method for multi-stage voltage sag |
US18/622,847 US20240250528A1 (en) | 2022-03-03 | 2024-03-29 | Methods and systems for controlling multi-stage voltage sags |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210201077.XA CN114611278B (zh) | 2022-03-03 | 2022-03-03 | 一种多阶段电压暂降状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114611278A true CN114611278A (zh) | 2022-06-10 |
CN114611278B CN114611278B (zh) | 2022-09-13 |
Family
ID=81860938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210201077.XA Active CN114611278B (zh) | 2022-03-03 | 2022-03-03 | 一种多阶段电压暂降状态估计方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20230305050A1 (zh) |
CN (1) | CN114611278B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115411729A (zh) * | 2022-09-30 | 2022-11-29 | 南方电网科学研究院有限责任公司 | 一种电力系统的电压稳定判别方法和装置 |
CN117743807A (zh) * | 2024-02-19 | 2024-03-22 | 四川大学 | 一种配电网多阶段电压暂降成因辨识方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060241881A1 (en) * | 2005-03-30 | 2006-10-26 | Gasperi Michael L | Networked power line parameter analysis method and system |
CN105842580A (zh) * | 2015-11-20 | 2016-08-10 | 武汉大学 | 一种适用于智能配电网的电压暂降源精确定位优化方法 |
CN109188138A (zh) * | 2018-09-11 | 2019-01-11 | 广东电网有限责任公司 | 一种电压暂降频次评估方法及装置 |
CN109523165A (zh) * | 2018-11-19 | 2019-03-26 | 东南大学 | 一种基于有限信息的电压暂降状态估计方法 |
CN110082639A (zh) * | 2019-05-15 | 2019-08-02 | 广东电网有限责任公司 | 一种基于最小幅值搜索的电压暂降监测点优化配置方法 |
CN112529432A (zh) * | 2020-12-17 | 2021-03-19 | 无锡泓大恒能科技发展有限公司 | 一种电压暂降严重程度评估方法、装置及电子设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120136598A1 (en) * | 2010-08-04 | 2012-05-31 | Vladimir Dmitriev-Zdorov | Optimization of Decoupling Device Choice for Electronic Design |
US9874593B2 (en) * | 2014-07-11 | 2018-01-23 | Abb Inc. | Decision support system for outage management and automated crew dispatch |
-
2022
- 2022-03-03 CN CN202210201077.XA patent/CN114611278B/zh active Active
- 2022-12-16 US US18/082,611 patent/US20230305050A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060241881A1 (en) * | 2005-03-30 | 2006-10-26 | Gasperi Michael L | Networked power line parameter analysis method and system |
CN105842580A (zh) * | 2015-11-20 | 2016-08-10 | 武汉大学 | 一种适用于智能配电网的电压暂降源精确定位优化方法 |
CN109188138A (zh) * | 2018-09-11 | 2019-01-11 | 广东电网有限责任公司 | 一种电压暂降频次评估方法及装置 |
CN109523165A (zh) * | 2018-11-19 | 2019-03-26 | 东南大学 | 一种基于有限信息的电压暂降状态估计方法 |
CN110082639A (zh) * | 2019-05-15 | 2019-08-02 | 广东电网有限责任公司 | 一种基于最小幅值搜索的电压暂降监测点优化配置方法 |
CN112529432A (zh) * | 2020-12-17 | 2021-03-19 | 无锡泓大恒能科技发展有限公司 | 一种电压暂降严重程度评估方法、装置及电子设备 |
Non-Patent Citations (3)
Title |
---|
YING WANGDENG: ""Optimized Placement of Voltage Sag Monitors Considering Distributed Generation Dominated Grids and Customers Demand"", 《FRONTIERS IN ENERGY RESEARCH》 * |
刘旭娜 等: ""电压暂降严重程度及其测度、不确定性评估方法"", 《中国电机工程学报》 * |
胡文曦 等: ""考虑多维特征刻画的电压暂降严重程度评估方法"", 《电网技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115411729A (zh) * | 2022-09-30 | 2022-11-29 | 南方电网科学研究院有限责任公司 | 一种电力系统的电压稳定判别方法和装置 |
CN115411729B (zh) * | 2022-09-30 | 2024-04-16 | 南方电网科学研究院有限责任公司 | 一种电力系统的电压稳定判别方法和装置 |
CN117743807A (zh) * | 2024-02-19 | 2024-03-22 | 四川大学 | 一种配电网多阶段电压暂降成因辨识方法 |
CN117743807B (zh) * | 2024-02-19 | 2024-04-26 | 四川大学 | 一种配电网多阶段电压暂降成因辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
US20230305050A1 (en) | 2023-09-28 |
CN114611278B (zh) | 2022-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114611278B (zh) | 一种多阶段电压暂降状态估计方法 | |
US20030040846A1 (en) | Stability prediction for an electric power network | |
CN110783946B (zh) | 用于定位微电网中相故障的方法 | |
Emami et al. | Tracking changes in the external network model | |
KR102099568B1 (ko) | 분산 전원이 연계된 배전 계통의 자동 고장 처리 장치 및 그 방법 | |
CN110635463B (zh) | 一种基于改进搜索保护与差动保护的微电网综合保护方法 | |
JP5147739B2 (ja) | 配電系統の地絡保護システム、配電系統の地絡保護方法、プログラム | |
CN113169558A (zh) | 用于嵌套式微电网的分散式错误数据减轻 | |
CN110176755B (zh) | 一种故障限流器配置方法及装置 | |
CN109861178B (zh) | 基于SCADA稳态信息的220kV电网复杂保护动作判断方法 | |
CN109546631B (zh) | 适用于同塔四回线路跨电压等级跨线故障的距离保护方法 | |
CN113541111B (zh) | 一种基于断线故障判据的输电线路断线故障判别方法 | |
CN114461982B (zh) | 一种输电线路保护特性识别及电压暂降持续时间估计方法 | |
CN113468758B (zh) | 考虑用户经济损失风险的暂降监测装置多目标优化配置方法 | |
CN112838610B (zh) | 一种基于二端口网络考虑中性点电流的mmc换流站接地故障分析方法 | |
CN107976612B (zh) | 一种基于单相接地故障线路跳闸信息的极性校验方法 | |
CN110672973A (zh) | 一种基于实时状态数据的广域单线接地故障检测方法 | |
CN113725855B (zh) | 一种基于电气DebtRank算法的电力系统脆弱线路辨识方法 | |
CN117743807B (zh) | 一种配电网多阶段电压暂降成因辨识方法 | |
CN113673085B (zh) | 一种考虑开关控制失效的配电网可靠性评估方法 | |
CN110429572B (zh) | 一种直流配电网极间故障快速保护方法 | |
CN112531641B (zh) | 基于直轴电流的有源配电网差动保护方法及系统 | |
Guo et al. | Single Phase Short Circuit Fault Location Method Based on Innovation Graph | |
Hashim et al. | Weak area analysis based on the apparent impedance and voltage indices | |
Hatefi Torshizi et al. | Improving Coordination and Operating Speed of Overcurrent Relay against Contingency of Presence of Distributed Generators |
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 |