发明名称 基于小波能量相对熵的高压直流输电线路单端保护方法
摘要 一种基于小波能量相对熵的高压直流输电线路单端保护方法,其步骤主要为:1)数据采集与预处理:电流测量装置采集电流信号,保护数据处理装置利用正交变换矩阵进行相模变换,得到线路边界的内侧1模电流x<sub>a</sub>(n)和线路边界的外侧1模电流x<sub>b</sub>(n)。2)小波能量相对熵的计算:分别对内侧1模电流向量X<sub>a</sub>和外侧1模电流向量X<sub>b</sub>进行小波变换,计算内侧1模电流向量X<sub>a</sub>与外侧1模电流向量X<sub>b</sub>之间的小波能量相对熵M。3)高压直流输电线路单端保护:判别小波能量相对熵M是否大于给定的保护阈值ε,若M大于ε,保护装置则判定故障发生在高压直流输电输电线路区内。该方法对线路故障的识别率高,保护动作更准确可靠,能更好的保证输电线路的安全稳定运行。
申请公布号 CN105591372A 申请公布日期 2016.05.18
申请号 CN201610137717.X 申请日期 2016.03.11
申请人 西南交通大学 发明人 林圣;高杉
分类号 H02H7/26(2006.01)I 主分类号 H02H7/26(2006.01)I
代理机构 成都博通专利事务所 51208 代理人 陈树明
主权项 一种基于小波能量相对熵的高压直流输电线路单端保护方法,其步骤为:A、数据采集与预处理电流测量装置在线路边界内侧采集到内侧正极线路电流i<sub>ap</sub>(n)、内侧负极线路电流i<sub>an</sub>(n),数据处理装置利用正交变换矩阵对内侧正极线路电流i<sub>ap</sub>(n)、内侧负极线路电流i<sub>an</sub>(n)进行相模变换,得到线路边界的内侧1模电流x<sub>a</sub>(n);电流测量装置在线路边界外侧采集到外侧正极线路电流i<sub>bp</sub>(n)、外侧负极线路电流i<sub>bn</sub>(n);数据处理装置利用正交变换矩阵对外侧正极线路电流i<sub>bp</sub>(n)、外侧负极线路电流i<sub>bn</sub>(n)进行相模变换,得到线路边界的外侧1模电流x<sub>b</sub>(n);其中n表示采样时刻,n=1,2,…,N,N是采样点总数;B、小波能量相对熵的计算B1、小波变换由内侧1模电流x<sub>a</sub>(n)构成内侧1模电流向量X<sub>a</sub>,X<sub>a</sub>=[x<sub>a</sub>(1),x<sub>a</sub>(2),…x<sub>a</sub>(n),…,x<sub>a</sub>(N)],对内侧1模电流向量X<sub>a</sub>进行离散小波分解,再进行单支重构,得到内侧1模电流重构系数矩阵D<sub>a</sub>,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>D</mi><mi>a</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>a</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000939146080000011.GIF" wi="1283" he="542" /></maths>其中d<sub>aj</sub>(n)是内侧1模电流向量X<sub>a</sub>在分解尺度j下、n时刻的重构系数,j为小波变换系数的分解尺度序号,j=1,2,…,J+1,J为小波分解层数;由外侧1模电流x<sub>b</sub>(n)构成外侧1模电流向量X<sub>b</sub>,X<sub>b</sub>=[x<sub>b</sub>(1),x<sub>b</sub>(2),…x<sub>b</sub>(n),…,x<sub>b</sub>(N)],对外侧1模电流向量X<sub>b</sub>进行离散小波分解,再进行单支重构后得到内侧1模电流重构系数矩阵D<sub>b</sub>,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>D</mi><mi>b</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd><mtd><mrow></mrow></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mi>K</mi></mtd><mtd><mrow><msub><mi>d</mi><mrow><mi>b</mi><mi>J</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000939146080000012.GIF" wi="1286" he="535" /></maths>其中d<sub>bj</sub>(n)是外侧1模电流向量X<sub>b</sub>在分解尺度j下、n时刻的重构系数;B2、计算小波能量相对熵计算内侧1模电流向量X<sub>a</sub>在分解尺度j下的能量E<sub>aj</sub>,<img file="FDA0000939146080000021.GIF" wi="373" he="158" />进而得到内侧1模电流向量X<sub>a</sub>的总能量E<sub>a</sub>,<img file="FDA0000939146080000022.GIF" wi="294" he="166" />进而得到分解尺度j下内侧1模电流向量X<sub>a</sub>的能量E<sub>aj</sub>与总能量E<sub>a</sub>之比p<sub>aj</sub>,<img file="FDA0000939146080000023.GIF" wi="268" he="165" />计算外侧1模电流向量X<sub>b</sub>在分解尺度j下的能量E<sub>bj</sub>,<img file="FDA0000939146080000024.GIF" wi="372" he="159" />进而得到外侧1模电流向量X<sub>b</sub>的总能量E<sub>b</sub>,<img file="FDA0000939146080000025.GIF" wi="301" he="166" />进而得到分解尺度j下外侧1模电流向量X<sub>b</sub>的能量E<sub>bj</sub>与总能量E<sub>b</sub>之比p<sub>bj</sub>,<img file="FDA0000939146080000026.GIF" wi="269" he="165" />再计算出,内侧1模电流向量X<sub>a</sub>相对于外侧1模电流向量X<sub>b</sub>的小波能量相对熵M<sub>ab</sub>,<img file="FDA0000939146080000027.GIF" wi="514" he="198" />外侧1模电流向量X<sub>b</sub>相对于内侧1模电流向量X<sub>a</sub>的小波能量相对熵M<sub>ba</sub>,<img file="FDA0000939146080000028.GIF" wi="508" he="199" />最后计算出,内侧1模电流向量X<sub>a</sub>与外侧1模电流向量X<sub>b</sub>之间的小波能量相对熵M,M=M<sub>ab</sub>+M<sub>ba</sub>C、高压直流输电线路单端保护若内侧1模电流向量X<sub>a</sub>与外侧1模电流向量X<sub>b</sub>之间的小波能量相对熵M大于给定的保护的阈值ε,数据处理装置就判断故障发生在高压直流输电线路的边界内,数据处理装置发出信号,保护控制装置切除边界内的线路。
地址 610031 四川省成都市二环路北一段111号