CN114152836A - 一种双端行波测距模量选择同步方法及装置 - Google Patents
一种双端行波测距模量选择同步方法及装置 Download PDFInfo
- Publication number
- CN114152836A CN114152836A CN202010937655.7A CN202010937655A CN114152836A CN 114152836 A CN114152836 A CN 114152836A CN 202010937655 A CN202010937655 A CN 202010937655A CN 114152836 A CN114152836 A CN 114152836A
- Authority
- CN
- China
- Prior art keywords
- modulus
- current
- maximum value
- mode
- traveling wave
- 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.)
- Pending
Links
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
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种双端行波测距模量选择同步方法,包括:获取三相电流行波数据,进行相模变换得到各个模量电流,对各模量电流进行小波分析求取模极大值集合,根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值,将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值,比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流,实现两侧模量选择的严格同步。本申请同时公开了相应的同步装置。本申请的方案能够有效解决发生多相不同时故障时,两侧因模量选择不一致而无法正确测距的问题。
Description
技术领域
本发明涉及输电线路故障定位技术领域,具体涉及一种双端行波测距模量选择同步方法及装置。
背景技术
输电线路双端行波测距是利用线路两端的初始行波波头时刻,借助光纤通道和卫星精确授时进行测距的故障定位技术,双端行波测距具有不需要识别故障点反射波、测距精度高的优点。理想情况下,各个模量电流的初始行波波头时刻相同,双端行波测距可以任意选择模量电流进行计算。但实际故障工况往往非常复杂,比如发展性故障,某一相先发生接地故障,随后发生两相接地故障,造成各个模量电流的初始行波波头时刻出现较大差异。这种情况下线路两侧如果仍然随意选择模量电流,可能造成选择的模量不一致,两侧的初始行波波头时刻差不能正确反映故障距离,导致测距结果的错误。
此外,双端行波测距的精度与初始行波波头时刻的识别精度密切相关,一般来说,行波波头越陡,小波变换模极大值幅值越大,波头时刻的提取精度越高,所以实际工程中应尽量选择行波波头上升速度最快的模量电流进行测距。
发明内容
本申请的目的,在于提供一种双端行波测距模量选择同步方法及装置,通过对两侧各模量电流模极大值的综合判断,实现两侧模量选择的严格同步,同时确保所选择模量电流的模极大值为最优,从而有效提高双端行波测距的精度。
为了达成上述目的,本申请的解决方案是:
一方面,本申请提供了一种双端行波测距模量选择同步方法,包括如下步骤:
获取三相电流行波数据,进行相模变换得到各个模量电流;
对各模量电流进行小波分析求取模极大值集合;
根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
优选的方案中,所述三相电流行波数据采样周期为1μs。
优选的方案中,所述各个模量电流包括0模电流i0、α模电流iα、β模电流iβ。
优选的方案中,所述相模变换采用Clarke变换、Karenbauer变换或者Wedpohl变换。
优选的方案中,所述对各模量电流进行小波分析求取模极大值集合计算步骤如下:
对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β)。
所述根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值计算步骤包括:
遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。系数K1的取值范围为[0.2,0.5]。优选的,系数K1的取值为0.3。
优选的方案中,两侧行波测距装置通过专用光纤或复用光纤交换各模量电流初始行波波头的模极大值。
优选的方案中,所述模极大值综合最优的模量选取方法为:
①求取两侧各模量电流初始行波波头模极大值的最小值lk:
②求取l0、lα、lβ三者中的最大值,该最大值对应的模量即为模极大值综合最优的模量。
优选的方案中,所述双端行波测距中初始行波波头时刻的计算采用所述模极大值综合最优的模量电流。
另一方面,本申请提供了一种双端行波测距模量选择同步装置,包括依次连接的采集单元、模极大值集合求取单元、初始行波波头模极大值计算单元、通讯单元和模量选取单元;
所述采集单元,用于获取三相电流行波数据,进行相模变换得到各个模量电流;
所述模极大值集合求取单元,用于对各模量电流进行小波分析求取模极大值集合;
所述初始行波波头模极大值计算单元,用于根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
所述通讯单元,用于将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
所述模量选取单元,用于比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
优选的,所述三相电流行波数据采样周期为1μs。
优选的,所述各个模量电流包括0模电流i0、α模电流iα、β模电流iβ。
优选的,所述相模变换采用Clarke变换、Karenbauer变换或者Wedpohl变换。
优选的,所述对各模量电流进行小波分析求取模极大值集合计算步骤如下:
对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β)。
所述根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值计算步骤包括:
遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。
优选的,所述模极大值综合最优的模量选取方法为:
①求取两侧各模量电流初始行波波头模极大值的最小值lk:
②求取l0、lα、lβ三者中的最大值,该最大值对应的模量即为模极大值综合最优的模量。
优选的,所述双端行波测距中初始行波波头时刻的计算采用所述模极大值综合最优的模量电流。
采用上述方案后,双端行波测距装置通过获取两侧各模量电流初始行波波头的模极大值,对各模量电流的模极大值进行综合比较,选出两侧模极大值均较大的模量电流进行测距,提高了波头时刻识别的精确度;同时由于两侧装置处理的数据和程序是完全相同的,保证了两侧模量选择的严格一致,有效解决发生多相不同时故障时,两侧可能因模量选择不一致而无法正确测距的问题。
附图说明
图1是本申请一种双端行波测距模量选择同步方法的实施例流程图之一。
图2是本申请一种双端行波测距模量选择同步方法的实施例流程图之二。
图3(a)是本申请一种双端行波测距模量选择同步方法的实施例分析结果之一。
图3(b)是本申请一种双端行波测距模量选择同步方法的实施例分析结果之二。
图4是本申请一种双端行波测距模量选择同步装置实施例示意图。
具体实施方式
下面结合附图和具体的实施例对本申请技术方案作进一步的详细描述,以使本领域的技术人员可以更好的理解本申请并能予以实施,但所举实施例不作为对本申请的限定。
如图1所示,一种双端行波测距模量选择同步方法的实施例,包括如下步骤:
S11:获取三相电流行波数据,进行相模变换得到各个模量电流;
S12:对各模量电流进行小波分析求取模极大值集合;
S13:根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
S14:将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
S15:比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
如图2所示的另一种双端行波测距模量选择同步方法的实施例,包括如下步骤:
步骤一:获取三相电流行波数据,进行相模变换得到各个模量电流。
优选地,所述三相电流行波数据采样周期为1μs,所述模量电流包括0模电流i0、α模电流iα、β模电流iβ,所述相模变换可以采用Clarke变换、Karenbauer变换或者Wedpohl变换。
步骤二:对各模量电流进行小波分析求取模极大值集合。
具体计算步骤如下:
①对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
②遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β)。
步骤三:根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
具体计算步骤如下:
①遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
②再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。系数K1的取值范围为[0.2,0.5]。
本实施例中系数K1取值为0.3,上述步骤则为再次遍历模极大值集合Mk,取第一个数值大于0.3mk的模极大值作为该模量电流初始行波波头的模极大值,记为hk。
步骤四:将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值。
优选地,两侧行波测距装置通过专用光纤或复用光纤交换各模量电流初始行波波头的模极大值。
步骤五:比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
模极大值综合最优的模量选取方法为:
①求取两侧各模量电流初始行波波头模极大值的最小值lk:
②求取l0、lα、lβ三者中的最大值,该最大值对应的模量即为模极大值综合最优的模量。
双端行波测距中初始行波波头时刻的计算采用上述模极大值综合最优的模量电流。
图3(a)、图3(b)分别为本申请一种双端行波测距模量选择同步方法的实施例某故障工况下线路两侧的三相行波电流及各模量小波分析结果,可以看到该故障工况为三相电流不同时故障,各模量电流的初始行波波头时刻不相同,下面假设图(a)为M侧、图(b)为N侧。由图可知,按照本申请方法,M侧各模量电流的模最大值m0、mα、mβ分别为:80.38A、123A、108.2A,各模量电流的初始行波波头模极大值h0、hα、hβ分别为62.12A、67.39A、103.7A;N侧各模量电流的模最大值m0、mα、mβ分别为:20.31A、15.9A、22.76A,各模量电流的初始行波波头模极大值h0、hα、hβ分别为13.82A、11.82A、19.12A。根据上述数据可求得两侧各模量电流初始行波波头模极大值的l0、lα、lβ最小值分别为13.82A、11.82A、19.12A,其中l0、lα、lβ的最大值为19.12A,对应的模量为β模,因此线路两侧都将选择β模分量进行双端行波测距。
如图4所示为本申请提供的一种双端行波测距模量选择同步装置实施例,包括依次连接的采集单元、模极大值集合求取单元、初始行波波头模极大值计算单元、通讯单元和模量选取单元。其中:
采集单元,用于获取三相电流行波数据,进行相模变换得到各个模量电流;
模极大值集合求取单元,用于对各模量电流进行小波分析求取模极大值集合;
初始行波波头模极大值计算单元,用于根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
通讯单元,用于将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
模量选取单元,用于比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
优选的装置实施例中,所述三相电流行波数据采样周期为1μs。
优选的装置实施例中,所述各个模量电流包括0模电流i0、α模电流iα、β模电流iβ。
优选的装置实施例中,所述相模变换采用Clarke变换、Karenbauer变换或者Wedpohl变换。
优选的装置实施例中,所述对各模量电流进行小波分析求取模极大值集合计算步骤如下:
对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β)。
所述根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值计算步骤包括:
遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。
优选的装置实施例中,所述模极大值综合最优的模量选取方法为:
①求取两侧各模量电流初始行波波头模极大值的最小值lk:
②求取l0、lα、lβ三者中的最大值,该最大值对应的模量即为模极大值综合最优的模量。
优选的装置实施例中,所述双端行波测距中初始行波波头时刻的计算采用所述模极大值综合最优的模量电流。
本申请通过对线路两侧的各模量电流的模极大值进行综合比较,实现了两侧模量选择的严格同步,能够有效解决发生多相不同时故障时,两侧可能因模量选择不一致而无法正确测距的问题;同时由于采用了模极大值综合最优的模量,提高了双端行波测距的精度。
本申请不受上述实施例子的限制,上述实施例中的描述只用于帮忙理解本申请的核心思想,凡是依据本申请的思想,对本申请进行修改或等同替换,在具体实施方式及应用范围上所做的任何改动,都应该属于本申请的保护范围。
Claims (17)
1.一种双端行波测距模量选择同步方法,其特征在于,包括如下步骤:
获取三相电流行波数据,进行相模变换得到各个模量电流;
对各模量电流进行小波分析求取模极大值集合;
根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
2.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:所述三相电流行波数据采样周期为1μs。
3.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:所述各个模量电流包括0模电流i0、α模电流iα、β模电流iβ。
4.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:所述相模变换采用Clarke变换、Karenbauer变换或者Wedpohl变换。
5.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:所述对各模量电流进行小波分析求取模极大值集合计算步骤如下:
对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β);
所述根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值计算步骤包括:
遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。
6.如权利要求5所述的双端行波测距模量选择同步方法,其特征在于:所述系数K1的取值范围为[0.2,0.5]。
7.如权利要求6所述的双端行波测距模量选择同步方法,其特征在于:所述系数K1的取值为0.3。
8.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:两侧行波测距装置通过专用光纤或复用光纤交换各模量电流初始行波波头的模极大值。
10.如权利要求1所述的双端行波测距模量选择同步方法,其特征在于:所述双端行波测距中初始行波波头时刻的计算采用所述模极大值综合最优的模量电流。
11.一种双端行波测距模量选择同步装置,其特征在于,包括依次连接的采集单元、模极大值集合求取单元、初始行波波头模极大值计算单元、通讯单元和模量选取单元;
所述采集单元,用于获取三相电流行波数据,进行相模变换得到各个模量电流;
所述模极大值集合求取单元,用于对各模量电流进行小波分析求取模极大值集合;
所述初始行波波头模极大值计算单元,用于根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值;
所述通讯单元,用于将各模量电流初始行波波头的模极大值传送给对侧,同时接收对侧各模量电流初始行波波头的模极大值;
所述模量选取单元,用于比较两侧的模极大值,选取模极大值综合最优的模量电流为双端行波测距模量电流。
12.如权利要求11所述的双端行波测距模量选择同步装置,其特征在于:所述三相电流行波数据采样周期为1μs。
13.如权利要求11所述的双端行波测距模量选择同步装置,其特征在于:所述各个模量电流包括0模电流i0、α模电流iα、β模电流iβ。
14.如权利要求11所述的双端行波测距模量选择同步装置,其特征在于:所述相模变换采用Clarke变换、Karenbauer变换或者Wedpohl变换。
15.如权利要求11所述的双端行波测距模量选择同步装置,其特征在于:所述对各模量电流进行小波分析求取模极大值集合计算步骤如下:
对各模量电流ik(k=0,α,β)进行小波分析,其中k代表模量,取值包括0模,α模,β模;
遍历各模量电流小波分析结果,求得各模量电流的模极大值集合Mk(k=0,α,β);
所述根据各模量电流的模极大值集合计算各模量电流初始行波波头的模极大值计算步骤包括:
遍历各模量电流的模极大值集合Mk(k=0,α,β),求得各模量电流的模最大值mk(k=0,α,β);
再次遍历各模量电流的模极大值集合Mk(k=0,α,β),取第一个数值大于K1*mk(k=0,α,β)的模极大值作为该模量电流初始行波波头的模极大值,记为hk(k=0,α,β),其中K1为系数,0<K1<1。
17.如权利要求11所述的双端行波测距模量选择同步装置,其特征在于:所述双端行波测距中初始行波波头时刻的计算采用所述模极大值综合最优的模量电流。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010937655.7A CN114152836A (zh) | 2020-09-08 | 2020-09-08 | 一种双端行波测距模量选择同步方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010937655.7A CN114152836A (zh) | 2020-09-08 | 2020-09-08 | 一种双端行波测距模量选择同步方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114152836A true CN114152836A (zh) | 2022-03-08 |
Family
ID=80462118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010937655.7A Pending CN114152836A (zh) | 2020-09-08 | 2020-09-08 | 一种双端行波测距模量选择同步方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114152836A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115754588A (zh) * | 2022-11-07 | 2023-03-07 | 国网江苏省电力有限公司南通市海门区供电分公司 | 一种节约通道资源的输电线路行波测距方法 |
-
2020
- 2020-09-08 CN CN202010937655.7A patent/CN114152836A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115754588A (zh) * | 2022-11-07 | 2023-03-07 | 国网江苏省电力有限公司南通市海门区供电分公司 | 一种节约通道资源的输电线路行波测距方法 |
CN115754588B (zh) * | 2022-11-07 | 2023-12-29 | 国网江苏省电力有限公司南通市海门区供电分公司 | 一种节约通道资源的输电线路行波测距方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8909069B2 (en) | Signal framing in a coherent optical receiver | |
CN106406174B (zh) | 一种多模块多通道采集同步系统及工作方法 | |
EP2899570B1 (en) | Method and device for digitalizing scintillation pulse | |
CN105259410B (zh) | 一种强噪声干扰下的欠采样波形的频率估计方法及其装置 | |
JP2010019625A (ja) | 事故点標定方法および事故点標定装置 | |
CN110750484B (zh) | 一种转速和多路振动通道数据同步采集系统及采集方法 | |
CN114152836A (zh) | 一种双端行波测距模量选择同步方法及装置 | |
CN110780110B (zh) | 一种基波电压过零点自动检测方法、系统及采样装置 | |
CN109870283B (zh) | 叶端定时传感器信号序列转位移序列的方法及系统 | |
WO2016131374A1 (zh) | 干扰信号记录装置及有线电视网损伤点定位的系统和方法 | |
CN108390752A (zh) | 信号接收方法 | |
CN107431614A (zh) | 用于自动偏移补偿的方法和装置 | |
CN106841918A (zh) | 一种利用单‑双端结合的电缆低阻故障定位方法 | |
US11336554B2 (en) | Universal semiconductor-based automatic highspeed serial signal testing method | |
WO2020207066A1 (zh) | 故障指示器装置、系统及配电网系统 | |
EP0717525A3 (en) | Method and device for serial communication | |
CN114740329A (zh) | 数据比对装置、数据比对方法及芯片测试仪 | |
CN114487695A (zh) | 基于单双端测距的线缆混合直流线路故障测距方法及系统 | |
CN113630296A (zh) | 一种lvds传输延迟窗口自动测试方法及系统 | |
CN116520095B (zh) | 故障测距方法、系统以及计算机可读存储介质 | |
CN207718361U (zh) | Pcm信号采集装置及系统 | |
CN112014669A (zh) | 一种合并单元角差、比差的在线检测方法及系统 | |
CN103199879A (zh) | 数字接收机信号的检测方法 | |
CN111585742B (zh) | 一种多路信号峰值同步检测系统 | |
CN109492185B (zh) | 一种采样数据处理方法及装置 |
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 |