CN114002559A - 一种柔性直流输电线路行波双端测距方法及*** - Google Patents
一种柔性直流输电线路行波双端测距方法及*** Download PDFInfo
- Publication number
- CN114002559A CN114002559A CN202111437331.8A CN202111437331A CN114002559A CN 114002559 A CN114002559 A CN 114002559A CN 202111437331 A CN202111437331 A CN 202111437331A CN 114002559 A CN114002559 A CN 114002559A
- Authority
- CN
- China
- Prior art keywords
- frequency
- traveling wave
- double
- signal
- calculating
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 55
- 238000000691 measurement method Methods 0.000 title claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 53
- 230000010355 oscillation Effects 0.000 claims abstract description 43
- 230000008878 coupling Effects 0.000 claims abstract description 21
- 238000010168 coupling process Methods 0.000 claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims description 27
- 238000005070 sampling Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000009466 transformation Effects 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 8
- 230000001902 propagating effect Effects 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 13
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 3
- 239000004020 conductor Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
本发明涉及一种柔性直流输电线路行波双端测距方法及***,属于电力***继电保护技术领域。本发明首先从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号;再分别提取双端行波耦合箱容性电流行波的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱;最后计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。本发明利用柔性直流输电线路两端的故障行波信号在线路两端平波电抗器边界和故障点来回折反射的等间隔分布频谱规律实现故障定位,无需识别行波波头性质、无需标定波到时刻、无需双端时钟同步,抗干扰能力、抗过渡电阻能力强,测距结果准确、可靠,具有良好的稳定性、鲁棒性。
Description
技术领域
本发明涉及一种柔性直流输电线路行波双端测距方法及***,属于电力***继电保护技术领域。
背景技术
柔性直流输电技术在新能源并网、电网异步区域互联、无源网络供电等多方面都具有突出的技术优势。随着新型电力电子装置、多电平模块化换流器的研发、有故障自清除功能的子模块拓扑、新型直流断路器拓扑的出现,进一步促进了柔性直流输电技术应用的推广。与传统直流输电相比,柔性直流换流站不存在换相失败问题,相比于传统直流输电更为可靠。与直流电缆相比,架空线路更适用于远距离输电,架设及维修较为方便,成本较低,但容易受到气象和环境影响而引起故障,且直流线路故障具有发展速度快,故障电流大的特点,因此,研究适用于柔性直流输电线路故障测距方法具有重要的理论意义和工程价值。
柔性直流输电线路故障测距方法主要借鉴交流输电线路和传统直流输电线路故障测距。主流测距方法按原理上的不同主要包括故障分析法、行波法和固有频率法。故障分析法根据线路参数和电气量之间的关系列写方程,通过优化求解实现故障测距,该方法对采样率要求较低,但依靠固定线路参数实现算法,其精度受线路参数品变特性影响较大。行波法测距主要利用故障行波的传播特性,其关键在于正确找出反应故障位置的波头,现有行波法测距技术研究主要从波头奇异性检测可靠性出发,提出了多种波头检测方法和手段,在波头奇异性可靠提取的前提下,行波法测距具有精度高,速度快的优势,但其波头奇异性检测受干扰影响严重,波头极性辨识不易、需要较高采样率。相比之下,固有频率法测距对采样率要求不高,其测距原理是根据故障行波在故障点与线路边界来回折反射的周期性规律,提取故障行波的固有频率主频,根据波速和固有频率主频、故障距离之间的数学关系计算故障距离。固有频率法测距的优势在于其计算结果比较稳定,无需进行波头的标定、抗干扰、抗噪声能力较强等,但是其测距精度受固有频率主频值大小影响,而该频率与***边界条件有关,与计算时窗长短有关,与波速有关,固有频率法测距结果会因此出现偏差。
发明内容
本发明要解决的技术问题是提供一种柔性直流输电线路行波双端测距方法及***,用以解决上述故障测距的问题。
本发明的技术方案是:一种柔性直流输电线路行波双端测距方法,具体步骤为:
Step1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号,构造容性电流测距信号。
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱。
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
所述Step1具体为:
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号。
Step1.2:通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
所述Step2具体为:
Step2.1:分别提取双端容性电流测距信号的自由振荡分量。
Step2.2:计算Step2.1所得信号的高阶奇次方后的信号。
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱。
Step2.4:判断两端获得的频谱是否为等间隔,若否,则改变数据时窗长度,返回执行Step2.3,若是,则分别计算两端频域信号的频差。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
所述Step3具体为:
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差。
Step3.2:按照式(1)求出双端频差的比值ΩMN,M表示双端柔性直流装置的一端,N表示双端柔性直流装置的另一端。
式(1)中,ΔωM、ΔωN分别表示M端和N端计算所得的角频差、ΔfM、ΔfN分别表示M端和N端计算所得的频差,双端频差和角频差分别满足式(2)、(3)的关系式。
Step3.3:利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%:
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l表示直流线路全长。
Step3.4:按线路全长计算校验频差值Δωl。
Step3.5:比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度,返回Step2.3。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
所述Step3.2还可以为,按照式(5)求取双端频差的比值ΩMN。
所述Step3.3还可以为,利用双端频差的比值ΩNM所构造的测距公式(6),得到测距结果(1-k)%。
所述Step3.5还可以为,比较校验频差值Δωl与ΔωN的比值与所得到的测距估计结果(1-k)%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度,返回Step2.3。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
一种柔性直流输电线路行波双端测距***,包括:
采集模块,用于从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号。
计算模块,用于构造容性电流测距信号,提取双端容性电流测距信号的自由振荡分量。
信号分析模块,分别计算两端测距信号的频域信号,获得自由振荡分量对应的等间隔频谱。
信号处理模块,用于计算频域信号的等间隔频差值,并计算两端频差值的比值。
测距模块,用于计算和记录故障测距结果。
所述采集模块包括:
变送单元,用于变换互感器二次侧的电流信号为行波装置A/D采集的信号。
采集单元,用于将电流模拟量信号转换成数字信号。
存储单元,用于按时标命名录波数据文件,并存储于本地存储器。
启动单元,用于判断波形突变是否大于设定的启动阈值,若是,则将电流信号存储成录波数据文件。
所述计算模块包括:
计算单元1,通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
计算单元2,分别计算双端容性电流测距信号自由振荡分量的高阶奇次方。
所述信号分析模块包括:
时频变换单元,对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱。
频谱判断单元,用于判断两端频域信号的频谱是否为等间隔。
所述信号处理模块包括:
频差计算单元,用于提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差。
比值计算单元,用于求出双端频差的比值ΩMN、ΩNM,M表示双端直流装置的一端,N表示双端直流装置的另一端。
校验计算单元,用于按线路全长计算校验频差值Δωl。
校验比较单元,比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度。
本发明的有益效果是:
本发明不受装置边界条件影响,规避了行波波头识别的不可靠、波到时刻的标定的不精确所带来的误差。
本发明利用双端频差比构造测距方程,测距精度不受波头畸变、波形缺损、首波头丢失的影响,数学上减少了波速衰减的影响,不依赖双端同步对时,对高阻故障等弱故障模态具有较好的鲁棒性。
本发明可利用行波耦合箱所获得的行波信号,作为现有行波测距算法的补充和辅助,测距结果以故障位置占全线长的百分比表示,线路长度不参与测距计算,其精度不受线路长度实际变化影响。
附图说明
图1是实施例1所对应的流程图;
图2是实施例1中所对应的Step1的具体步骤流程图;
图3是实施例1中所对应的Step2的具体步骤流程图;
图4是实施例1中所对应的Step3的具体步骤流程图;
图5是实施例1中所涉及的***功能框图;
图6是实施例中的柔性直流输电***拓扑示意图;
图7是实施例中的M端和N端行波耦合箱所获得的电流行波示意图;
图8是实施例2中的M端和N端所对应的高阶奇次方后信号的等间隔频谱图。
图9是实施例3中柔性直流输电线路两端行波耦合箱分别采集的容性电流行波信号iMc、iNc的图;
图10是实施例3中极空间模量容性电流行波信号图;
图11是实施例3中双端容性电流测距信号的自由振荡分量图;
图12是实施例3中M端高阶奇次方信号的频谱图;
图13是实施例3中N端高阶奇次方信号的频谱图;
图14是实施例4中柔性直流输电线路两端行波耦合箱分别采集的容性电流行波信号iMc、iNc的图;
图15是实施例4中极空间模量容性电流行波信号图;
图16是实施例4中双端容性电流测距信号的自由振荡分量图;
图17是实施例4中M端高阶奇次方信号的频谱图;
图18是实施例4中N端高阶奇次方信号的频谱图;
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
实施例1:如图1所示,一种柔性直流输电线路行波双端测距方法,具体步骤为:
Step1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号,构造容性电流测距信号。
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱。
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
如图2所示,Step1具体为:
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号。
Step1.2:通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
如图3所示,Step2具体为:
Step2.1:分别提取双端容性电流测距信号的自由振荡分量。
Step2.2:计算Step2.1所得信号的高阶奇次方后的信号。
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱。
Step2.4:判断两端获得的频谱是否为等间隔,若否,则改变数据时窗长度,返回执行Step2.3,若是,则分别计算两端频域信号的频差。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
如图4所示,Step3具体为:
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差。
Step3.2:按照式(1)求出双端频差的比值ΩMN,M表示双端柔性直流装置的一端,N表示双端柔性直流装置的另一端。
式(1)中,ΔωM、ΔωN分别表示M端和N端计算所得的角频差、ΔfM、ΔfN分别表示M端和N端计算所得的频差,双端频差和角频差分别满足式(2)、(3)的关系式。
Step3.3:利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%。
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l表示直流线路全长。
Step3.4:按线路全长计算校验频差值Δωl。
Step3.5:比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度,返回Step2.3。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
如图5所示,一种柔性直流输电线路行波双端测距***,包括:
采集模块101,用于从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号。
计算模块102,用于构造容性电流测距信号。提取双端容性电流测距信号的自由振荡分量。
信号分析模块103,分别计算两端测距信号的频域信号,获得自由振荡分量对应的等间隔频谱。
信号处理模块104,用于计算频域信号的等间隔频差值,并计算两端频差值的比值。
测距模块105,用于计算和记录故障测距结果。
所述采集模块101包括:
变送单元1011,用于变换互感器二次侧的电流信号为行波装置A/D采集的信号。
采集单元1012,用于将电流模拟量信号转换成数字信号。
存储单元1013,用于按时标命名录波数据文件,并存储于本地存储器。
启动单元1014,用于判断波形突变是否大于设定的启动阈值,若是,则将电流信号存储成录波数据文件。
所述计算模块102包括:
计算单元1021A,通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
计算单元1021B,分别计算双端容性电流测距信号自由振荡分量的高阶奇次方。
所述信号分析模块103包括:
时频变换单元1031,对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱。
频谱判断单元1032,用于判断两端频域信号的频谱是否为等间隔。
所述信号处理模块104包括:
频差计算单元1041,用于提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差。
比值计算单元1042,用于求出双端频差的比值ΩMN、ΩNM,M表示双端直流装置的一端,N表示双端直流装置的另一端。
校验计算单元1043,用于按线路全长计算校验频差值Δωl。
校验比较单元1044,比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度。
实施例2:如图6所示,在PSCAD/EMTDC环境下搭建的±500kV真双极双端柔性直流输电***(MMC-HVDC)。整流侧为M端,逆变侧为N端。双端MMC换流器单个桥臂各有200个半桥子模块,桥臂电抗器Larm=100mH。线路全长l=500km,架空线采用频变参数模型,保留避雷线,线路双端限流电抗器L=150mH,采样频率取100kHz。假设在线路距离M端100km处发生单极金属性接地故障情况下,利用柔性直流输电线路行波双端测距方法进行测距,其具体步骤为:
Step1:从柔性直流输电线路M、N两端行波耦合箱分别采集容性电流行波信号构造容性电流测距信号。
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱。
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
所述Step1具体包括:
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号iMc、iNc。
Step1.2:通过式(5)所示的Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号
式(5)中,i1、i0分别表示空间模量和地模分量,I+、I_分别表示正极电流和负极电流,根据式(4)计算可分别得到M端和N端极空间模量容性电流行波信号iMc1、iNc1,对M端和N端分别用下角标M和N表示,容性电流用下角标c表示。
按式(6)求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
Δic1(k)=ic1(k+1)-ic1(k) (6)
分别求得M端和N端的容性电流测距信号用△iMc1、△iNc1表示,如图7所示。
所述Step2具体包括:
Step2.1:分别提取双端容性电流测距信号的自由振荡分量。按行波在全线长范围传播所需要的时间τ为基准,分别取首故障行波到达前后时间窗为τ,2τ,4τ长度的测距信号,其包含故障行波在故障点到量测端来回多次折反射的振荡特性,作为自由振荡分量。
Step2.2:计算Step2.1所得信号的α次方信号Signalc1,如式(7)所示。
Signalc1(k)=(Δic1(k))α,α=2n+1,n=1,2,3... (7)
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱,如图8所示。
Step2.4:判断两端获得的频谱是否为等间隔,若是,则分别计算两端频域信号的频差。
所述改变数据时窗长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
所述Step3具体包括:
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差△fM=372.64Hz,△fN=248.38Hz。
Step3.2:按照式(1)求出双端频差的比值ΩMN=1.5003,M表示双端柔性直流装置的整流侧,N表示双端柔性直流装置的逆变侧。
式(1)中,ΔfM、ΔfN分别表示M端和N端计算所得的频差,
Step3.3利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%=39.99%。
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l=500km表示直流线路全长。
Step3.4:按线路全长计算校验频差值△fl。
Step3.5:比较校验频差值△fl与△fM的比值为39.98%。
39.98%与所得到的测距估计结果k%之差为-0.0005,小于所设定的校验阈值±0.01,输出测距结果k%=39.99%。
实施例3:在PSCAD/EMTDC环境下搭建的±500kV真双极双端柔性直流输电***(MMC-HVDC)。整流侧为M端,逆变侧为N端。双端MMC换流器单个桥臂各有200个半桥子模块,桥臂电抗器Larm=100mH。线路全长l=500km,架空线采用频变参数模型,保留避雷线,线路双端限流电抗器L=150mH,采样频率取200kHz。假设在线路中点处发生单极接地故障,过渡电阻为500欧姆的情况下,利用柔性直流输电线路行波双端测距方法进行测距,其具体步骤为:
Step1:从柔性直流输电线路M、N两端行波耦合箱分别采集容性电流行波信号构造容性电流测距信号。
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号iMc、iNc,如图9所示。
Step1.2:通过式(5)所示的Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,如图10所示:
式(5)中,i1、i0分别表示空间模量和地模分量,I+、I_分别表示正极电流和负极电流,根据式(5)计算可分别得到M端和N端极空间模量容性电流行波信号iMc1、iNc1,对M端和N端分别用下角标M和N表示,容性电流用下角标c表示。
按式(6)求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
Δic1(k)=ic1(k+1)-ic1(k) (6)
分别求得M端和N端的容性电流测距信号用△iMc1、△iNc1表示,如图11所示。
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱。
Step2.1:分别提取双端容性电流测距信号的自由振荡分量。按行波在全线长范围传播所需要的时间τ为基准,分别取首故障行波到达前后时间窗为τ,2τ,4τ长度的测距信号,其包含故障行波在故障点到量测端来回多次折反射的振荡特性,作为自由振荡分量,如图11所示。Step2.2:计算Step2.1所得信号的α次方信号Signalc1,如式(7)所示。
Signalc1(k)=(Δic1(k))α,α=2n+1,n=1,2,3... (7)
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱,如图12、图13所示。
Step2.4:判断两端获得的频谱是否为等间隔,若是,则分别计算两端频域信号的频差。
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差△fM=296.86Hz,△fN=296.87Hz。
Step3.2:按照式(1)求出双端频差的比值ΩMN=1,M表示双端柔性直流装置的整流侧,N表示双端柔性直流装置的逆变侧。
式(1)中,ΔfM、ΔfN分别表示M端和N端计算所得的频差。
Step3.3:利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%=50%。
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l=500km表示直流线路全长。
Step3.4:按线路全长计算校验频差值△fl。
Step3.5:比较校验频差值△fl与△fM的比值为50%。
50%与所得到的测距估计结果k%之差为0,小于所设定的校验阈值±0.01,输出测距结果k%=50%。
实施例4:在PSCAD/EMTDC环境下搭建的±500kV真双极双端柔性直流输电***(MMC-HVDC)。整流侧为M端,逆变侧为N端。双端MMC换流器单个桥臂各有200个半桥子模块,桥臂电抗器Larm=100mH。线路全长l=500km,架空线采用频变参数模型,保留避雷线,线路双端限流电抗器L=150mH,采样频率取1MHz。假设在线路距离M近端(50km)处发生双极短路故障情况下,利用柔性直流输电线路行波双端测距方法进行测距,其具体步骤为:
Step1:从柔性直流输电线路M、N两端行波耦合箱分别采集容性电流行波信号构造容性电流测距信号。
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号iMc、iNc,如图14所示。
Step1.2:通过式(5)所示的Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,如图15所示。
式(5)中,i1、i0分别表示空间模量和地模分量,I+、I_分别表示正极电流和负极电流,根据式(5)计算可分别得到M端和N端极空间模量容性电流行波信号iMc1、iNc1,对M端和N端分别用下角标M和N表示,容性电流用下角标c表示。
按式(6)求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
Δic1(k)=ic1(k+1)-ic1(k) (6)
分别求得M端和N端的容性电流测距信号用△iMc1、△iNc1表示,如图16所示。
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱。
Step2.1:分别提取双端容性电流测距信号的自由振荡分量。按行波在全线长范围传播所需要的时间τ为基准,分别取首故障行波到达前后时间窗为τ,2τ,4τ长度的测距信号,其包含故障行波在故障点到量测端来回多次折反射的振荡特性,作为自由振荡分量,如图16所示。
Step2.2:计算Step2.1所得信号的α次方信号Signalc1。如式(7)所示。
Signalc1(k)=(Δic1(k))α,α=2n+1,n=1,2,3... (7)
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱,如图17、图18所示。
Step2.4:判断两端获得的频谱是否为等间隔,若是,则分别计算两端频域信号的频差。
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差△fM=1489.6Hz,△fN=165.87Hz。
Step3.2:按照式(1)求出双端频差的比值ΩMN=8.98,M表示双端柔性直流装置的整流侧,N表示双端柔性直流装置的逆变侧。
式(1)中,ΔfM、ΔfN分别表示M端和N端计算所得的频差。
Step3.3:利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%=10.02%。
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l=500km表示直流线路全长。
Step3.4:按线路全长计算校验频差值△fl。
Step3.5:比较校验频差值△fl与△fM的比值为10%。
10%与所得到的测距估计结果k%之差为-0.00017,小于所设定的校验阈值±0.01,输出测距结果k%=10.02%。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (12)
1.一种柔性直流输电线路行波双端测距方法,其特征在于:
Step1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号,构造容性电流测距信号;
Step2:分别提取双端容性电流测距信号的自由振荡分量,计算其频域信号,获得自由振荡分量对应的等间隔频谱;
Step3:计算频域信号的等间隔频差值,计算两端频差值的比值,利用该比值进行故障测距。
2.根据权利要求1所述的柔性直流输电线路行波双端测距方法,其特征在于,所述Step1具体为:
Step1.1:从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号;
Step1.2:通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号。
3.根据权利要求1所述的柔性直流输电线路行波双端测距方法,其特征在于,所述Step2具体为:
Step2.1:分别提取双端容性电流测距信号的自由振荡分量;
Step2.2:计算Step2.1所得信号的高阶奇次方后的信号;
Step2.3:对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱;
Step2.4:判断两端获得的频谱是否为等间隔,若否,则改变数据时窗长度,返回执行Step2.3,若是,则分别计算两端频域信号的频差。
4.根据权利要求1所述的柔性直流输电线路行波双端测距方法,其特征在于,所述Step3具体为:
Step3.1:分别提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差;
Step3.2:按照式(1)求出双端频差的比值ΩMN,M表示双端柔性直流装置的一端,N表示双端柔性直流装置的另一端;
式(1)中,ΔωM、ΔωN分别表示M端和N端计算所得的角频差、ΔfM、ΔfN分别表示M端和N端计算所得的频差,双端频差和角频差分别满足式(2)、(3)的关系式;
Step3.3:利用双端频差的比值ΩMN所构造的测距公式(4),得到测距估计结果k%:
式(4)中,xf表示故障点相对于M端的距离,k%表示故障点相对于M端的距离占线路总长的百分比,l表示直流线路全长;
Step3.4:按线路全长计算校验频差值Δωl;
Step3.5:比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度,返回Step2.3。
6.根据权利要求5所述的柔性直流输电线路行波双端测距方法,其特征在于:所述Step3.5还可以为,比较校验频差值Δωl与ΔωN的比值与所得到的测距估计结果(1-k)%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度,返回Step2.3。
7.根据权利要求3、4或6所述的柔性直流输电线路行波双端测距方法,其特征在于:所述改变数据时窗,其长度以行波在线路全长范围传播的时间τ为单位,每次改变Kτ,K=2n,n=1,2,3…。
8.一种柔性直流输电线路行波双端测距***,其特征在于,包括:
采集模块,用于从柔性直流输电线路两端行波耦合箱分别采集容性电流行波信号;
计算模块,用于构造容性电流测距信号,提取双端容性电流测距信号的自由振荡分量;
信号分析模块,分别计算两端测距信号的频域信号,获得自由振荡分量对应的等间隔频谱;
信号处理模块,用于计算频域信号的等间隔频差值,并计算两端频差值的比值;
测距模块,用于计算和记录故障测距结果。
9.根据权利要求8所述的柔性直流输电线路行波双端测距方法,其特征在于,所述采集模块包括:
变送单元,用于变换互感器二次侧的电流信号为行波装置A/D采集的信号;
采集单元,用于将电流模拟量信号转换成数字信号;
存储单元,用于按时标命名录波数据文件,并存储于本地存储器;
启动单元,用于判断波形突变是否大于设定的启动阈值,若是,则将电流信号存储成录波数据文件。
10.根据权利要求8所述的柔性直流输电线路行波双端测距方法,其特征在于,所述计算模块包括:
计算单元1,通过Karenbauer变换矩阵解耦直流线路,求取极空间模量容性电流行波信号,并求取该信号的相邻采样点的值之差,作为前一个采样点的新值,以此构造容性电流测距信号;
计算单元2,分别计算双端容性电流测距信号自由振荡分量的高阶奇次方。
11.根据权利要求8所述的柔性直流输电线路行波双端测距方法,其特征在于,所述信号分析模块包括:
时频变换单元,对高阶奇次方信号进行离散傅里叶变换分别得到两端高阶奇次方信号的频谱;
频谱判断单元,用于判断两端频域信号的频谱是否为等间隔。
12.根据权利要求8所述的柔性直流输电线路行波双端测距方法,其特征在于,所述信号处理模块包括:
频差计算单元,用于提取双端等间隔频谱的谱峰位置的频率值,对任意相邻两个频率值作差,得到两端各自频谱所对应相等的频差;
比值计算单元,用于求出双端频差的比值ΩMN、ΩNM,M表示双端直流装置的一端,N表示双端直流装置的另一端。
校验计算单元,用于按线路全长计算校验频差值Δωl;
校验比较单元,比较校验频差值Δωl与ΔωM的比值与所得到的测距估计结果k%之差是否小于所设定的校验阈值,若是,则输出测距结果,若否,则改变数据时窗长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111437331.8A CN114002559B (zh) | 2021-11-29 | 2021-11-29 | 一种柔性直流输电线路行波双端测距方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111437331.8A CN114002559B (zh) | 2021-11-29 | 2021-11-29 | 一种柔性直流输电线路行波双端测距方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114002559A true CN114002559A (zh) | 2022-02-01 |
CN114002559B CN114002559B (zh) | 2022-10-14 |
Family
ID=79930851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111437331.8A Active CN114002559B (zh) | 2021-11-29 | 2021-11-29 | 一种柔性直流输电线路行波双端测距方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114002559B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115963358A (zh) * | 2023-03-13 | 2023-04-14 | 昆明理工大学 | 一种混合三端柔性直流输电线路故障测距方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3812432A1 (de) * | 1988-04-14 | 1989-10-26 | Siemens Ag | Verfahren und vorrichtung zur bestimmung des fehlerortes auf einer elektrischen leitung |
CN101860000A (zh) * | 2010-05-14 | 2010-10-13 | 河南电力试验研究院 | 一种输电线路单相重合闸前的永久性故障快速识别方法 |
CN105738759A (zh) * | 2014-12-12 | 2016-07-06 | 国家电网公司 | 一种基于暂态录波数据的直流输电线路故障定位方法 |
CN107861024A (zh) * | 2017-10-29 | 2018-03-30 | 天津大学 | 柔性直流输电线路的电压行波差动保护时域计算方法 |
CN112083272A (zh) * | 2020-08-18 | 2020-12-15 | 昆明理工大学 | 一种基于护层电流行波固有频率频差法的10kV电缆故障测距方法 |
-
2021
- 2021-11-29 CN CN202111437331.8A patent/CN114002559B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3812432A1 (de) * | 1988-04-14 | 1989-10-26 | Siemens Ag | Verfahren und vorrichtung zur bestimmung des fehlerortes auf einer elektrischen leitung |
CN101860000A (zh) * | 2010-05-14 | 2010-10-13 | 河南电力试验研究院 | 一种输电线路单相重合闸前的永久性故障快速识别方法 |
CN105738759A (zh) * | 2014-12-12 | 2016-07-06 | 国家电网公司 | 一种基于暂态录波数据的直流输电线路故障定位方法 |
CN107861024A (zh) * | 2017-10-29 | 2018-03-30 | 天津大学 | 柔性直流输电线路的电压行波差动保护时域计算方法 |
CN112083272A (zh) * | 2020-08-18 | 2020-12-15 | 昆明理工大学 | 一种基于护层电流行波固有频率频差法的10kV电缆故障测距方法 |
Non-Patent Citations (2)
Title |
---|
ZHAO ZHAO 等: "VSC-HVDC transmission line protection based on current natural frequency", 《2016 IEEE PES ASIA-PACIFIC POWER AND ENERGY ENGINEERING CONFERENCE (APPEEC)》 * |
刘兵 等: "基于小波分析频差性的电力电缆故障测距", 《电力自动化设备》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115963358A (zh) * | 2023-03-13 | 2023-04-14 | 昆明理工大学 | 一种混合三端柔性直流输电线路故障测距方法及*** |
CN115963358B (zh) * | 2023-03-13 | 2023-08-04 | 昆明理工大学 | 一种混合三端柔性直流输电线路故障测距方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN114002559B (zh) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114152840B (zh) | 一种lcc-mmc混合直流输电线路故障测距方法及*** | |
CN103513159A (zh) | 一种直流接地极线路上的故障测距方法及装置 | |
CN103901324A (zh) | 一种基于单端故障信息的配电网混合线路组合式测距方法 | |
CN109541392A (zh) | 一种适用于柔性直流输电***的单端故障测距方法 | |
CN116338525B (zh) | 一种风电交流送出线路故障测距方法及*** | |
CN111766470B (zh) | 高压直流输电线路的故障定位方法、***及直流输电线路 | |
CN107621591B (zh) | 一种基于零模行波波速变化特性的输电线路迭代测距方法 | |
CN114966301A (zh) | 一种三端混合高压直流输电线路故障测距方法及*** | |
CN115015686A (zh) | 一种lcc-vsc混合高压直流输电线路故障测距方法及*** | |
CN112083272A (zh) | 一种基于护层电流行波固有频率频差法的10kV电缆故障测距方法 | |
CN114002559B (zh) | 一种柔性直流输电线路行波双端测距方法及*** | |
CN113945806B (zh) | 一种柔性直流输电线路单端故障测距方法及*** | |
CN115097253A (zh) | 一种mmc-hvdc直流输电线路故障测距方法及*** | |
CN115015687A (zh) | 一种四端环状柔性直流电网故障测距方法及*** | |
CN113311287B (zh) | 一种应用于双端vsc控制的直流配电网的故障测距方法 | |
CN114487695A (zh) | 基于单双端测距的线缆混合直流线路故障测距方法及*** | |
CN114578184B (zh) | 一种直流输电线路双端行波频差比值故障测距方法及*** | |
CN110927539B (zh) | 一种单端暂态能量谱相似性的配电网故障区段定位方法 | |
CN114200246B (zh) | 一种直流输电线路故障测距方法及*** | |
CN114460412B (zh) | 一种直流输电线路单端故障测距方法及*** | |
CN114137356A (zh) | 一种直流输电线路测距方法及*** | |
CN114325226A (zh) | 一种输电线路单端自适应校正的高频故障定位方法及*** | |
CN114047406A (zh) | 一种基于派克变换的半波长线路双端故障测距方法 | |
CN111308267A (zh) | 基于c-evt的混合多端直流输电线路行波测距方法 | |
CN111665419A (zh) | 基于波头对称性的t接线路故障单端测距方法及*** |
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 |