发明名称 一种基于小波域行波信号色散校正的双端行波测距方法
摘要 本发明公开了一种基于小波域行波信号色散校正的双端行波测距方法,包括以下步骤:获取输电线M、N端的三相电流行波信号;利用凯伦贝尔变换解耦三相电流行波信号;计算故障点到所述输电线两端的距离;计算其在各离散频率下的传播系数和校正函数;在小波变换的某一尺度下,对解耦后的行波信号进行校正;最后,基于校正后的行波信号重新计算故障距离。本发明解决了行波色散对行波测距带来的影响这一行波故障测距领域的难点问题,提高了行波故障测距的精度,不仅可以减少线路故障带来的经济损失保证电网安全稳定运行,而且具有重大的社会效益和经济效益。
申请公布号 CN105866618A 申请公布日期 2016.08.17
申请号 CN201610179620.5 申请日期 2016.03.25
申请人 华北电力大学(保定) 发明人 贾惠彬;李明舒
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 一种基于小波域行波信号色散校正的双端行波测距方法,其特征在于:包括以下步骤:步骤1:在输电线M端和N端分别安装电流行波信号提取装置,获取输电线M端和N端的三相电流行波信号i<sub>M</sub>(t)和i<sub>N</sub>(t);步骤2:设定模量m为1或2,利用凯伦贝尔变换解耦输电线M、N端的三相电流行波信号i<sub>M</sub>(t)、i<sub>N</sub>(t),得到解耦后的行波信号分别为i<sub>m,M</sub>(t)、i<sub>m,N</sub>(t):i<sub>m,M</sub>(t)=T<sup>‑1</sup>i<sub>M</sub>(t)Ti<sub>m,N</sub>(t)=T<sup>‑1</sup>i<sub>N</sub>(t)T         (1)其中,T为凯伦贝尔变换:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000950076840000011.GIF" wi="1262" he="219" /></maths>步骤3:设定行波信号的采样频率为f<sub>s</sub>,小波变换的点数为2<sup>n</sup>,且n为整数,对输电线M、N端的解耦后行波信号i<sub>m,M</sub>(t)、i<sub>m,N</sub>(t)分别采样后进行快速傅立叶变换后得到输电线M、N端的解耦后FFT变换信号i<sub>m,M</sub>(ω<sub>i</sub>)、i<sub>m,N</sub>(ω<sub>i</sub>),<img file="FDA0000950076840000012.GIF" wi="194" he="118" />为离散频率,且0≤i≤2<sup>n</sup>‑1;步骤4:计算小波变换函数,设定尺度a的小波变换函数为,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msqrt><mi>a</mi></msqrt><mo>&CenterDot;</mo><mi>&Psi;</mi><mrow><mo>(</mo><msub><mi>a&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000950076840000013.GIF" wi="1461" he="87" /></maths>其中Ψ(aω<sub>i</sub>)为尺度a下的小波函数;步骤5:通过快速傅立叶反变换计算所述输电线M、N端的解耦后行波信号i<sub>m,M</sub>(t)、i<sub>m,N</sub>(t)在尺度a下的小波分解:WT<sub>a</sub>i<sub>m,M</sub>(t)=IFFT(i<sub>m,M</sub>(ω<sub>i</sub>)H(ω<sub>i</sub>))WT<sub>a</sub>i<sub>m,N</sub>(t)=IFFT(i<sub>m,N</sub>(ω<sub>i</sub>)H(ω<sub>i</sub>))       (4)步骤6:根据双端行波故障定位原理计算故障点到所述输电线M、N端的距离l<sub>M</sub>和l<sub>N</sub>;步骤7:计算输电线M、N端的解耦后行波信号i<sub>m,M</sub>(t)、i<sub>m,N</sub>(t)在各离散频率<img file="FDA0000950076840000021.GIF" wi="206" he="119" />0≤i≤2<sup>n‑1</sup>‑1下的传播系数Q<sub>M</sub>(ω′<sub>i</sub>)、Q<sub>N</sub>(ω′<sub>i</sub>):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>M</mi></msub><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&lambda;</mi><mi>m</mi></msub><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><msub><mi>l</mi><mi>M</mi></msub></mrow></msup></mrow>]]></math><img file="FDA0000950076840000022.GIF" wi="364" he="68" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>N</mi></msub><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&lambda;</mi><mi>m</mi></msub><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><msub><mi>l</mi><mi>N</mi></msub></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000950076840000023.GIF" wi="1182" he="70" /></maths>其中,λ<sub>m</sub>(ω′<sub>i</sub>)为频变参数,由Carson公式计算得到;步骤8:将所述输电线M、N端的解耦后行波信号i<sub>m,M</sub>(t)、i<sub>m,N</sub>(t)在各离散频率<img file="FDA0000950076840000024.GIF" wi="211" he="125" />0≤i≤2<sup>n‑1</sup>‑1下的传播系数Q<sub>M</sub>(ω′<sub>i</sub>)、Q<sub>N</sub>(ω′<sub>i</sub>),关于<img file="FDA0000950076840000025.GIF" wi="155" he="117" />镜像对称扩展成2<sup>n</sup>点的扩展传播系数Q<sub>M</sub>(ω<sub>i</sub>)和Q<sub>N</sub>(ω<sub>i</sub>);计算输电线M、N端的解耦后行波信号校正函数<img file="FDA0000950076840000026.GIF" wi="160" he="71" />和<img file="FDA0000950076840000027.GIF" wi="179" he="69" />步骤9:在根据下式计算输电线M、N端色散校正后的小波变换WTi<sub>m,M</sub>(t)和WTi<sub>m,N</sub>(t):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>WTi</mi><mrow><mi>m</mi><mo>,</mo><mi>M</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mi>F</mi><mi>F</mi><mi>T</mi><mo>(</mo><msub><mi>I</mi><mrow><mi>m</mi><mo>,</mo><mi>M</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>Q</mi><mi>M</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow>]]></math><img file="FDA0000950076840000028.GIF" wi="982" he="79" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>WTi</mi><mrow><mi>m</mi><mo>,</mo><mi>N</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mi>F</mi><mi>F</mi><mi>T</mi><mo>(</mo><msub><mi>I</mi><mrow><mi>m</mi><mo>,</mo><mi>N</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>Q</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>H</mi><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000950076840000029.GIF" wi="1622" he="86" /></maths>步骤10:根据双端行波定位原理,计算故障点到所述输电线M、N端的精确距离l′<sub>M</sub>和l′<sub>N</sub>。
地址 071000 河北省保定市永华北大街619号