发明名称 基于分形和小波变换的电力电缆故障识别方法
摘要 本发明公开了一种基于分形和小波变换的电力电缆故障识别系统及方法,其系统包括依次相接的电流检测电路模块、信号调理电路模块、数据采集卡和主控计算机,电流检测电路模块包括A、B、C相霍尔电流传感器,信号调理电路模块包括A、B、C相I/V转换电路模块,A、B、C相信号放大电路模块和A、B、C相滤波电路模块;其方法包括步骤:一、信号实时检测及同步上传,二、信号采集及存储,三、电缆故障识别。本发明设计新颖合理,采用霍尔电流传感器对电缆电流进行检测,检测精度高、稳定性和可靠性高,将分形法与小波分析法相结合,能够简单、迅速、准确、有效地识别出电缆短路故障的类型,实用性强,使用效果好,便于推广使用。
申请公布号 CN103323748B 申请公布日期 2015.06.10
申请号 CN201310269904.X 申请日期 2013.06.29
申请人 西安科技大学 发明人 汪梅;杨晨;付燕
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 西安创知专利事务所 61213 代理人 谭文琰
主权项 一种基于分形和小波变换的电力电缆故障识别方法,该方法采用的识别系统包括依次相接的电流检测电路模块(1)、信号调理电路模块(2)、数据采集卡(3)和主控计算机(4),所述电流检测电路模块(1)包括用于对三相电缆的A相电流进行实时检测的A相霍尔电流传感器(1‑1)、用于对三相电缆的B相电流进行实时检测的B相霍尔电流传感器(1‑2)和用于对三相电缆的C相电流进行实时检测的C相霍尔电流传感器(1‑3),所述信号调理电路模块(2)包括用于对A相霍尔电流传感器(1‑1)输出的信号进行调理的A相电流信号调理电路模块(2‑1)、用于对B相霍尔电流传感器(1‑2)输出的信号进行调理的B相电流信号调理电路模块(2‑2)和用于对C相霍尔电流传感器(1‑3)输出的信号进行调理的C相电流信号调理电路模块(2‑3),所述A相电流信号调理电路模块(2‑1)由依次相接的A相I/V转换电路模块(2‑11)、A相信号放大电路模块(2‑12)和A相滤波电路模块(2‑13)构成,所述B相电流信号调理电路模块(2‑2)由依次相接的B相I/V转换电路模块(2‑21)、B相信号放大电路模块(2‑22)和B相滤波电路模块(2‑23)构成,所述C相电流信号调理电路模块(2‑3)由依次相接的C相I/V转换电路模块(2‑31)、C相信号放大电路模块(2‑32)和C相滤波电路模块(2‑33)构成,其特征在于:该方法包括以下步骤:步骤一、信号实时检测及同步上传:通过A相霍尔电流传感器(1‑1)对被检测电缆的A相电流进行实时检测,并将实时所检测的A相电流信号i<sub>a</sub>依次通过A相I/V转换电路模块(2‑11)、A相信号放大电路模块(2‑12)和A相滤波电路模块(2‑13)进行I/V转换、放大和滤波处理后同步上传至数据采集卡(3);通过B相霍尔电流传感器(1‑2)对被检测电缆的B相电流进行实时检测,并将实时所检测的B相电流信号i<sub>b</sub>依次通过B相I/V转换电路模块(2‑21)、B相信号放大电路模块(2‑22)和B相滤波电路模块(2‑23)进行I/V转换、放大和滤波处理后同步上传至数据采集卡(3);通过C相霍尔电流传感器(1‑3)对被检测电缆的C相电流进行实时检测,并将实时所检测的C相电流信号i<sub>c</sub>依次通过C相I/V转换电路模块(2‑31)、C相信号放大电路模块(2‑32)和C相滤波电路模块(2‑33)进行I/V转换、放大和滤波处理后同步上传至数据采集卡(3);步骤二、信号采集及存储:数据采集卡(3)对A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>进行采集并相应进行A/D转换后,同步传送至主控计算机(4);同时,通过主控计算机(4)对A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>的采样点数量、所有采样点和各采样点所对应的采样时刻进行存储;步骤三、电缆故障识别:当被检测电缆发生故障时,数据采集卡(3)传送来的A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>为故障信号,此时通过主控计算机(4)对由数据采集卡(3)传送来的A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>进行分析处理,并相应得出被检测电缆的故障类型属于A相接地短路故障、B相接地短路故障、C相接地短路故障、AB两相短路接地故障、AB相间短路故障、BC两相短路接地故障、BC相间短路故障、AC两相短路接地故障、AC相间短路故障和ABC三相短路故障中的哪一种,其分析处理过程如下:步骤301、采用分形法进行电缆故障初识别,具体过程为:步骤3011、电流波形图绘制:所述主控计算机(4)调用电流波形绘制模块绘制出A相电流信号i<sub>a</sub>随采样时间t变化的A相电流波形图、B相电流信号i<sub>b</sub>随采样时间t变化的B相电流波形图和C相电流信号i<sub>c</sub>随采样时间t变化的C相电流波形图;步骤3012、分形维数计算:首先,假设A相电流波形是欧氏空间的有界点集F<sub>A</sub>,用边长为ε<sub>A</sub>的小方格覆盖该点集F<sub>A</sub>时,假设N<sub>A</sub>(ε<sub>A</sub>)为边长为ε<sub>A</sub>时其中包含点集F<sub>A</sub>中的点的非空小方格数,所述主控计算机(4)调用分形维数计算模块且根据公式<img file="FDA0000685816190000021.GIF" wi="408" he="154" />计算出A相分形维数D<sub>A</sub>;接着,假设B相电流波形是欧氏空间的有界点集F<sub>B</sub>,用边长为ε<sub>B</sub>的小方格覆盖该点集F<sub>B</sub>时,假设N<sub>B</sub>(ε<sub>B</sub>)为边长为ε<sub>B</sub>时其中包含点集F<sub>B</sub>中的点的非空小方格数,所述主控计算机(4)调用分形维数计算模块且根据公式<img file="FDA0000685816190000031.GIF" wi="403" he="154" />计算出B相分形维数D<sub>B</sub>;然后,假设C相电流波形是欧氏空间的有界点集F<sub>C</sub>,用边长为ε<sub>C</sub>的小方格覆盖该点集F<sub>C</sub>时,假设N<sub>C</sub>(ε<sub>C</sub>)为边长为ε<sub>C</sub>时其中包含点集F<sub>C</sub>中的点的非空小方格数,所述主控计算机(4)调用分形维数计算模块且根据公式<img file="FDA0000685816190000032.GIF" wi="408" he="157" />计算出C相分形维数D<sub>C</sub>;步骤3013、三相分形维数比对及电缆故障初识别:所述主控计算机(4)将A相分形维数D<sub>A</sub>、B相分形维数D<sub>B</sub>和C相分形维数D<sub>C</sub>进行比对,并根据比对结果分五种情形对电缆故障进行初识别:当D<sub>B</sub>&lt;D<sub>C</sub>&lt;D<sub>A</sub>时,判定为A相接地短路故障,识别结束并输出识别结果;当D<sub>A</sub>&lt;D<sub>B</sub>&lt;D<sub>C</sub>时,判定为C相接地短路故障,识别结束并输出识别结果;当D<sub>C</sub>&lt;D<sub>A</sub>&lt;D<sub>B</sub>时,判定为B相接地短路故障、AB两相短路接地故障或AB相间短路故障,执行步骤302;当D<sub>A</sub>&lt;D<sub>C</sub>&lt;D<sub>B</sub>时,判定为ABC三相短路故障、BC两相短路接地故障或BC相间短路故障,执行步骤302;当D<sub>B</sub>&lt;D<sub>A</sub>&lt;D<sub>C</sub>时,判定为AC两相短路接地故障或AC相间短路故障,执行步骤302;步骤302、所述主控计算机(4)调用凯伦布尔变换模块将A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>从相空间变换到模空间,得到0模电流分量i<sub>0</sub>、α模电流分量i<sub>α</sub>和β模电流分量i<sub>β</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>i</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><msub><mi>i</mi><mi>a</mi></msub><mo>+</mo><msub><mi>i</mi><mi>b</mi></msub><mo>+</mo><msub><mi>i</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000685816190000033.GIF" wi="368" he="124" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>i</mi><mi>&alpha;</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><msub><mi>i</mi><mi>a</mi></msub><mo>-</mo><msub><mi>i</mi><mi>b</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000685816190000034.GIF" wi="289" he="122" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>i</mi><mi>&beta;</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><msub><mi>i</mi><mi>a</mi></msub><mo>-</mo><msub><mi>i</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000685816190000035.GIF" wi="296" he="125" /></maths>步骤303、采用离散小波变换进行电缆故障再识别,具体过程为:步骤3031、离散小波变换及小波系数求解:首先,所述主控计算机(4)选取0模电流分量i<sub>0</sub>作为一维信号f(n<sub>0</sub>)并调用离散小波变换模块对一维信号f(n<sub>0</sub>)进行离散小波变换,相应求得离散小波变换后的各层小波系数,各层所述小波系数包括各层近似系数和各层细节系数,其中各层所述细节系数记作d<sub>j,k</sub>,其中,n<sub>0</sub>=0,1,2,…N‑1且N为一维信号f(n<sub>0</sub>)的采样序列中的采样点数量,一维信号f(n<sub>0</sub>)的采样序列中的采样点数量、所有采样序列点和各采样序列点所对应的采样时刻均与步骤二中主控计算机(4)所存储的A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>的采样点数量、所有采样点和各采样点所对应的采样时刻一一对应;j=1,2,…,J且J为离散小波变换的层数,k=0,1,2…,N‑1且k为一维信号f(n<sub>0</sub>)的采样序列中的N个采样序列点的序号;接着,所述主控计算机(4)选取α模电流分量i<sub>α</sub>作为一维信号f(n<sub>α</sub>)并调用离散小波变换模块对一维信号f(n<sub>α</sub>)进行离散小波变换,相应求得离散小波变换后的各层小波系数,各层所述小波系数包括各层近似系数和各层细节系数,其中各层所述细节系数记作d′<sub>j,k</sub>,其中,n<sub>α</sub>=0,1,2,…N‑1且N为一维信号f(n<sub>α</sub>)的采样序列中的采样点数量,一维信号f(n<sub>α</sub>)的采样序列中的采样点数量、所有采样序列点和各采样序列点所对应的采样时刻均与步骤二中主控计算机(4)所存储的A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>的采样点数量、所有采样点和各采样点所对应的采样时刻一一对应;j=1,2,…,J且J为离散小波变换的层数,k=0,1,2…,N‑1且k为一维信号f(n<sub>α</sub>)的采样序列中的N个采样序列点的序号;然后,所述主控计算机(4)选取β模电流分量i<sub>β</sub>作为一维信号f(n<sub>β</sub>)并调用离散小波变换模块对一维信号f(n<sub>β</sub>)进行离散小波变换,相应求得离散小波变换后的各层小波系数,各层所述小波系数包括各层近似系数和各层细节系数,其中各层所述细节系数记作d″<sub>j,k</sub>,其中,n<sub>β</sub>=0,1,2,…N‑1且N为一维信号f(n<sub>β</sub>)的采样序列中的采样点数量,一维信号f(n<sub>β</sub>)的采样序列中的采样点数量、所有采样序列点和各采样序列点所对应的采样时刻均与步骤二中主控计算机(4)所存储的A相电流信号i<sub>a</sub>、B相电流信号i<sub>b</sub>和C相电流信号i<sub>c</sub>的采样点数量、所有采样点和各采样点所对应的采样时刻一一对应;j=1,2,…,J且J为离散小波变换的层数,k=0,1,2…,N‑1且k为一维信号f(n<sub>β</sub>)的采样序列中的N个采样序列点的序号;步骤3032、模极大值点检测:首先,按照一维信号f(n<sub>0</sub>)的采样序列的采样顺序,由前至后对第J层细节系数d<sub>J,k</sub>的模极大值点进行检测与记录,且检测得出当k=i时,d<sub>J,k</sub>取得模极大值点,即0模电流分量i<sub>0</sub>的初始行波模极大值I<sub>0</sub>=d<sub>J,i</sub>;接着,按照一维信号f(n<sub>α</sub>)的采样序列的采样顺序,由前至后对第J层细节系数d′<sub>J,k</sub>的模极大值点进行检测与记录,且检测得出当k=i时,d′<sub>J,k</sub>取得模极大值点,即α模电流分量i<sub>α</sub>的初始行波模极大值I<sub>α</sub>=d′<sub>J,i</sub>;然后,按照一维信号f(n<sub>β</sub>)的采样序列的采样顺序,由前至后对第J层细节系数d″<sub>J,k</sub>的模极大值点进行检测与记录,且检测得出当k=i时,d″<sub>J</sub><sub>,k</sub>取得模极大值点,即β模电流分量i<sub>β</sub>的初始行波模极大值I<sub>β</sub>=d″<sub>J,i</sub>;其中,i=0,1,2…,N‑1;步骤3033、三个初始行波模极大值比对及电缆故障再识别:所述主控计算机(4)分D<sub>C</sub>&lt;D<sub>A</sub>&lt;D<sub>B</sub>、D<sub>A</sub>&lt;D<sub>C</sub>&lt;D<sub>B</sub>和D<sub>B</sub>&lt;D<sub>A</sub>&lt;D<sub>C</sub>三种情况将0模电流分量i<sub>0</sub>的初始行波模极大值I<sub>0</sub>、α模电流分量i<sub>α</sub>的初始行波模极大值I<sub>α</sub>和β模电流分量i<sub>β</sub>的初始行波模极大值I<sub>β</sub>进行比对,并根据比对结果对电缆故障进行再识别,具体为:当D<sub>C</sub>&lt;D<sub>A</sub>&lt;D<sub>B</sub>时,首先判断是否有I<sub>0</sub>=0,当I<sub>0</sub>=0时,判断是否有I<sub>α</sub>=2I<sub>β</sub>,当I<sub>0</sub>=0且I<sub>α</sub>=2I<sub>β</sub>时,判断为AB相间短路故障;当I<sub>0</sub>≠0时,首先判断是否有I<sub>β</sub>=0且|I<sub>0</sub>|=|I<sub>α</sub>|,当I<sub>0</sub>≠0、I<sub>β</sub>=0且|I<sub>0</sub>|=|I<sub>α</sub>|时,判断为B相接地短路故障;然后判断是否有I<sub>0</sub>=2I<sub>β</sub>‑I<sub>α</sub>,当I<sub>0</sub>≠0且I<sub>0</sub>=2I<sub>β</sub>‑I<sub>α</sub>时,判断为AB两相短路接地故障;识别结束并输出识别结果;当D<sub>A</sub>&lt;D<sub>C</sub>&lt;D<sub>B</sub>时,首先判断是否有I<sub>0</sub>=0,当I<sub>0</sub>=0时,判断是否有I<sub>α</sub>=‑I<sub>β</sub>,当I<sub>0</sub>=0且I<sub>α</sub>=‑I<sub>β</sub>时,判断为BC相间短路故障;当I<sub>0</sub>≠0时,判断是否有I<sub>0</sub>=‑I<sub>α</sub>‑I<sub>β</sub>,当I<sub>0</sub>≠0且I<sub>0</sub>=‑I<sub>α</sub>‑I<sub>β</sub>时,判断为BC两相短路接地故障;当I<sub>0</sub>≠0且I<sub>0</sub>≠‑I<sub>α</sub>‑I<sub>β</sub>时,判断为ABC三相短路故障;识别结束并输出识别结果;当D<sub>B</sub>&lt;D<sub>A</sub>&lt;D<sub>C</sub>时,首先判断是否有I<sub>0</sub>=0,当I<sub>0</sub>=0时,判断是否有<img file="FDA0000685816190000062.GIF" wi="221" he="129" />当I<sub>0</sub>=0且<img file="FDA0000685816190000061.GIF" wi="194" he="134" />时,判断为AC相间短路故障;当I<sub>0</sub>≠0时,判断是否有I<sub>0</sub>=2I<sub>α</sub>‑I<sub>β</sub>,当I<sub>0</sub>≠0且I<sub>0</sub>=2I<sub>α</sub>‑I<sub>β</sub>时,判断为AC两相短路接地故障;识别结束并输出识别结果。
地址 710054 陕西省西安市雁塔中路58号