CN104749497B - 对超声波放电检测信号处理后的音频数据可视化的方法 - Google Patents
对超声波放电检测信号处理后的音频数据可视化的方法 Download PDFInfo
- Publication number
- CN104749497B CN104749497B CN201410722389.0A CN201410722389A CN104749497B CN 104749497 B CN104749497 B CN 104749497B CN 201410722389 A CN201410722389 A CN 201410722389A CN 104749497 B CN104749497 B CN 104749497B
- Authority
- CN
- China
- Prior art keywords
- dft
- voice data
- fft
- domain
- real
- 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.)
- Active
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种对超声波放电检测信号处理后的音频数据可视化的方法。包含有如下步骤:首先对通过检测器存储的音频数据进行处理得到其离散时域图,一是对离散时域图经过FFT变换和频率平移得到原超声波的频谱图,二是对离散时域图进行取绝对值处理和低频率重采样得到原超声波的幅值相位图。本发明能够得到超声波的数据统计图谱,便于数据直观显示,更好的分析电力设备绝缘介质发生局部放电原因。
Description
技术领域
本发明涉及电力行业输变电技术领域,具体涉及到一种用于电力设备检测的对超声波放电检测信号处理后的音频数据可视化的方法。
背景技术
电力设备绝缘介质发生局部放电时,不但产生高频脉冲电信号,同时会产生超声波。超声波以球面波的方式向四周传播。超声波检测器可以把接收到的超声波信号经过外差处理然后转化成人耳可听的声音信号。这种可听的信号可以用耳机收听,也可以储存下来供其它人员使用。超声波检测器还可以显示超声波信号的强度,用于保存或比较。一般都把信号的强度作为历史资料保存。为了更好的分析电力设备绝缘介质发生局部放电,就要得到原始的超声波数据统计图谱,这就需要研发一种对超声波局部放电检测信号外差处理后存储的音频信号进行还原的可视化方法。
发明内容
本发明所要解决的技术问题是针对上述存在的不足,提供一种能够得到超声波的数据统计图谱,便于数据直观显示,更好的分析电力设备绝缘介质发生局部放电原因的对超声波放电检测信号处理后的音频数据可视化的方法。
本发明解决其技术问题所采用的技术方案是:
对超声波放电检测信号处理后的音频数据可视化的方法,包含有如下步骤:
首先利用传感器收集电力设备绝缘介质放电时产生的超声波信号,并通过放大器进行处理,传入调制器;将调制器连入本地振荡器,在调制器内对信号进行外差法处理,将局部放电的超声波信号经过超声波检测器的放大和外差处理后按一定的采样率和分辨率采集得到音频数据;其特征在于:
将得到的音频数据按采样率计算采样间隔,设置横轴的时间,取线性图,得到音频数据的离散时域图。
对所得到的离散时域图谱进行FFT变换,得到音频数据的频谱图,然后变换为柱状图,其中,FFT含义如下:f(t)是t的周期函数,如果t满足狄里赫莱条件:在一个周期内具有有限个间断点,且在这些间断点上,函数是有限值;在一个周期内具有有限个极值点;绝对可积。则有下式成立,称为积分运算f(t)的傅立叶变换;
为了在科学计算和数字信号处理等领域使用计算机进行傅里叶变换,必须将函数定义在离散点上而非连续域内,且须满足有限性或周期性条件;这种情况下,序列的离散傅立叶变换(discreteFouriertransform,DFT)为
FFT是一种DFT高效算法,称为快速傅立叶变换(fastFouriertransform)。由上式计算DFT对于X[k]的每个K值,需要进行4N次实数相乘和(4N-2)次相加,对于N个k值,共需N*N乘和N(4N-2)次实数相加。改进DFT算法,减小它的运算量,利用DFT中e-j2π/N的周期性和对称性,使整个DFT的计算变成一系列迭代运算,可大幅度提高运算过程和运算量,这就是FFT的基本思想。信号是实数序列,任何实数都可看成虚部为零的复数,利用复数据FFT对实数据进行有效计算。一个N点FFT同时计算两个N点实序列的DFT,设X1(n),X2(n)是彼此独立的两个N点实序列,且X1(k)=DFT[X1(n)],X2(k)=DFT[X2(n)],可通过一次FFT运算同时获得X1(k),X2(k)。算法如下:
首先将X1(n),X2(n)分别当作一复序列的实部及虚部,
令X(n)=X1(n)+jX2(n)
通过FFT运算可获得x(n)的DFT值,X(k)=DFT[X1(n)]+jDFT[X2(n)]=X1(k)+jX2(k)利用离散付里叶变换的共轭对称性
有了X(n)的FFT运算结果X(k),由上式即可得到X1(k),X2(k)的值。
将音频数据频谱图的柱状图平移本振频率fo个单位即得到原超声波频谱图。
截取音频数据离散时域图的单周期图;对得到的单周期离散时域图纵轴取绝对值;
使采样点数不变,根据重采样频率,重新计算采样时间和采样间隔,根据采样间隔重新设置时间轴,得到低频率重采样图;
把得到的低频率重采样图的时间轴变换成相位轴,关系如下:
p=ft×360°
f为被检测电力设备的工频交流电频率,然后变换为柱状图,即得到幅值相位图。
本发明的有益效果在于:
本发明方法能够得到超声波的数据统计图谱,便于数据直观显示,更好的分析电力设备绝缘介质发生局部放电原因。
附图说明
图1是本发明实施例超声波检测示意图;
图2是本发明实施例音频数据的离散时域图;
图3是本发明实施例音频数据的频谱图;
图4是本发明实施例音频数据频谱图的柱状形式;
图5是本发明实施例超声波数据的频谱图;
图6是本发明实施例音频数据离散时域的单周期图;
图7是本发明实施例音频数据单周期离散时域的绝对值图;
图8是本发明实施例音频数据的低频率重采样图;
图9是本发明实施例超声波数据的幅值相位图。
具体实施方式
下面结合具体实施方式,对本发明作进一步的说明:
如图1至图9所示的对超声波放电检测信号处理后的音频数据可视化的方法,包含有如下步骤:
如图1所示,首先利用传感器收集电力设备绝缘介质放电时产生的超声波信号,并通过放大器进行处理,传入调制器;将调制器连入本地振荡器,在调制器内对信号进行外差法处理,将局部放电的超声波信号经过超声波检测器的放大和外差处理后按一定的采样率和分辨率采集得到音频数据。
将得到的音频数据按采样率计算采样间隔,设置横轴的时间,取线性图,得到音频数据的离散时域图,如图2所示。
对所得到的离散时域图谱进行FFT变换,得到音频数据的频谱图,如图3所示,然后变换为柱状图,如图4所示。
其中,FFT含义如下:f(t)是t的周期函数,如果t满足狄里赫莱条件:在一个周期内具有有限个间断点,且在这些间断点上,函数是有限值;在一个周期内具有有限个极值点;绝对可积。则有下式成立,称为积分运算f(t)的傅立叶变换。
为了在科学计算和数字信号处理等领域使用计算机进行傅里叶变换,必须将函数定义在离散点上而非连续域内,且须满足有限性或周期性条件。这种情况下,序列的离散傅立叶变换(discreteFouriertransform,DFT)为
FFT是一种DFT高效算法,称为快速傅立叶变换(fastFouriertransform)。由上式计算DFT对于X[k]的每个K值,需要进行4N次实数相乘和(4N-2)次相加,对于N个k值,共需N*N乘和N(4N-2)次实数相加。改进DFT算法,减小它的运算量,利用DFT中e-j2π/N的周期性和对称性,使整个DFT的计算变成一系列迭代运算,可大幅度提高运算过程和运算量,这就是FFT的基本思想。信号是实数序列,任何实数都可看成虚部为零的复数,利用复数据FFT对实数据进行有效计算。一个N点FFT同时计算两个N点实序列的DFT,设X1(n),X2(n)是彼此独立的两个N点实序列,且X1(k)=DFT[X1(n)],X2(k)=DFT[X2(n)],可通过一次FFT运算同时获得X1(k),X2(k)。算法如下:
首先将X1(n),X2(n)分别当作一复序列的实部及虚部,
令X(n)=X1(n)+jX2(n)
通过FFT运算可获得x(n)的DFT值,X(k)=DFT[X1(n)]+jDFT[X2(n)]=X1(k)+jX2(k)利用离散付里叶变换的共轭对称性
有了X(n)的FFT运算结果X(k),由上式即可得到X1(k),X2(k)的值。
将音频数据频谱图的柱状图向右平移图1中的本振频率fo(38KHz)个单位即得到原超声波频谱图,如图5所示。
截取音频数据离散时域图的单周期图,如图6所示。
对得到的单周期离散时域图纵轴取绝对值,如图7所示。
使采样点数不变,根据重采样频率,重新计算采样时间和采样间隔,根据采样间隔重新设置时间轴,得到低频率重采样图,如图8所示。
把得到的低频率重采样图的时间轴变换成相位轴,关系如下:
p=ft×360°
f为被检测电力设备的工频交流电频率,然后变换为柱状图,即得到幅值相位图,如图9所示。
以上说明仅为本发明的应用实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明申请专利范围所作的等效变化,仍属本发明的保护范围。
Claims (1)
1.对超声波放电检测信号处理后的音频数据可视化的方法,包含有如下步骤:
首先利用传感器收集电力设备绝缘介质放电时产生的超声波信号,并通过放大器进行处理,传入调制器;将调制器连入本地振荡器,在调制器内对信号进行外差法处理,将局部放电的超声波信号经过超声波检测器的放大和外差处理后按一定的采样率和分辨率采集得到音频数据;其特征在于:
将得到的音频数据按采样率计算采样间隔,设置横轴的时间,取线性图,得到音频数据的离散时域图;
对所得到的离散时域图谱进行FFT变换,得到音频数据的频谱图,然后变换为柱状图,其中,FFT含义如下:f(t)是t的周期函数,如果t满足狄里赫莱条件:在一个周期内具有有限个间断点,且在这些间断点上,函数是有限值;在一个周期内具有有限个极值点;绝对可积;则有下式成立,称为积分运算f(t)的傅立叶变换;
将函数定义在离散点上而非连续域内,且须满足有限性或周期性条件;这种情况下,序列的离散傅立叶变换为
由上式计算DFT对于X[k]的每个K值,需要进行4N次实数相乘和(4N-2)次相加,对于N个k值,共需N*N乘和N(4N-2)次实数相加;改进DFT算法,减小它的运算量,利用DFT中e-j2π/N的周期性和对称性,使整个DFT的计算变成一系列迭代运算;信号是实数序列,任何实数都可看成虚部为零的复数,利用复数据FFT对实数据进行有效计算;一个N点FFT同时计算两个N点实序列的DFT,设X1(n),X2(n)是彼此独立的两个N点实序列,且X1(k)=DFT[X1(n)],X2(k)=DFT[X2(n)],可通过一次FFT运算同时获得X1(k),X2(k);算法如下:
首先将X1(n),X2(n)分别当作一复序列的实部及虚部,
令X(n)=X1(n)+jX2(n)
通过FFT运算可获得x(n)的DFT值,X(k)=DFT[X1(n)]+jDFT[X2(n)]=X1(k)+jX2(k)
利用离散付里叶变换的共轭对称性
有了X(n)的FFT运算结果X(k),由上式即可得到X1(k),X2(k)的值;
将音频数据频谱图的柱状图平移本振频率fo个单位即得到原超声波频谱图;
截取音频数据离散时域图的单周期图;对得到的单周期离散时域图纵轴取绝对值;
使采样点数不变,根据重采样频率,重新计算采样时间和采样间隔,根据采样间隔重新设置时间轴,得到低频率重采样图;
把得到的低频率重采样图的时间轴变换成相位轴,关系如下:
p=ft×360°
f为被检测电力设备的工频交流电频率;然后变换为柱状图,即得到幅值相位图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410722389.0A CN104749497B (zh) | 2014-12-02 | 2014-12-02 | 对超声波放电检测信号处理后的音频数据可视化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410722389.0A CN104749497B (zh) | 2014-12-02 | 2014-12-02 | 对超声波放电检测信号处理后的音频数据可视化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104749497A CN104749497A (zh) | 2015-07-01 |
CN104749497B true CN104749497B (zh) | 2016-05-04 |
Family
ID=53589478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410722389.0A Active CN104749497B (zh) | 2014-12-02 | 2014-12-02 | 对超声波放电检测信号处理后的音频数据可视化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104749497B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI582420B (zh) * | 2015-08-26 | 2017-05-11 | 行政院原子能委員會核能研究所 | 超音波探頭性能檢測裝置及方法 |
CN106199612A (zh) * | 2016-06-22 | 2016-12-07 | 杨越 | 无人船监控区域相控阵超声波低频成像方法 |
CN112652037A (zh) * | 2017-03-05 | 2021-04-13 | 杭州小影创新科技股份有限公司 | 一种实时音乐频谱矢量图形的绘制方法 |
CN107990974B (zh) * | 2017-05-16 | 2021-11-09 | 萨姆株式会社 | 辐射超声波可视化方法及执行辐射超声波可视化方法的电子装置 |
CN109490723A (zh) * | 2018-11-06 | 2019-03-19 | 红相股份有限公司 | 一种超声波可视化检测方法及超声波可视化检测仪 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103038599A (zh) * | 2010-04-16 | 2013-04-10 | U·E·系统公司 | 机载超声波频谱和图像生成 |
CN203773013U (zh) * | 2014-03-24 | 2014-08-13 | 三泰电力技术(南京)股份有限公司 | 一种局部放电检测器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060209632A1 (en) * | 2002-11-12 | 2006-09-21 | U-E Systems, Inc. | General purpose signal converter |
-
2014
- 2014-12-02 CN CN201410722389.0A patent/CN104749497B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103038599A (zh) * | 2010-04-16 | 2013-04-10 | U·E·系统公司 | 机载超声波频谱和图像生成 |
CN203773013U (zh) * | 2014-03-24 | 2014-08-13 | 三泰电力技术(南京)股份有限公司 | 一种局部放电检测器 |
Also Published As
Publication number | Publication date |
---|---|
CN104749497A (zh) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104749497B (zh) | 对超声波放电检测信号处理后的音频数据可视化的方法 | |
Chen et al. | Compound faults detection of rotating machinery using improved adaptive redundant lifting multiwavelet | |
Li et al. | Research on test bench bearing fault diagnosis of improved EEMD based on improved adaptive resonance technology | |
Hong et al. | Compound faults diagnosis based on customized balanced multiwavelets and adaptive maximum correlated kurtosis deconvolution | |
CN104188689B (zh) | 基于超声回波射频信号的组织位移估算方法和系统 | |
CN108051189B (zh) | 一种旋转机械故障特征提取方法及装置 | |
CN103267571B (zh) | 采用单面声压和振速测量的实时声场分离方法 | |
CN102519725B (zh) | 通过非线性冗余提升小波包处理轴承设备振动信号的方法 | |
CN103197141A (zh) | 一种测量电力系统信号频率及谐波参数的方法 | |
CN103393435A (zh) | 一种胎心音信号包络的取得方法及装置 | |
CN106908663A (zh) | 一种基于小波变换的电动汽车充电谐波辨识方法 | |
CN103646011A (zh) | 一种基于线性调频z变换的信号频谱细化方法 | |
CN105445541A (zh) | 一种任意频率下自适应功率计算方法 | |
CN104089699A (zh) | 一种变电站设备声音重建算法 | |
Ma et al. | Analysis and design of modified window shapes for S-transform to improve time–frequency localization | |
CN102008302A (zh) | 一种多频率同步信号源的合成方法 | |
CN102998528A (zh) | 动态正弦畸变信号中频率分量的提取方法及其实现装置 | |
CN104778342A (zh) | 一种基于小波奇异熵的心音特征提取方法 | |
CN103412188B (zh) | 基于贝塞尔函数与Toeplitz算法的SFM信号参数估计方法 | |
CN112328956A (zh) | 一种强频变信号时频分析方法 | |
CN105606893B (zh) | 基于空间平滑修正music的电力间谐波检测方法 | |
CN104459318A (zh) | 电力谐波的测量方法和系统 | |
CN104483563A (zh) | 电力信号的同步采样方法及系统 | |
Liu et al. | Power system harmonic analysis based on windowed FFT and wavelet transform | |
Ukte et al. | Two empirical methods for improving the performance of statistical multirate high-resolution signal reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |