变电站噪声分离方法技术

一种变电站噪声分离方法，涉及音频处理技术领域，所解决的是单独测量高压变电站三种噪声的技术问题。该方法首先对高压变电站的噪声实施多层小波变换得到六组小波系数；其次将低频小波系数经过通带梳状滤波器之后做小波反变换分离出本体噪声；同时将低频小波系数经过阻带梳状滤波器之后与另五组小波系数做小波反变换，得到包含电晕噪声及风扇噪声的混合信号；再利用高频小波系数检测电晕噪声，在没有电晕噪声的时间段估计风扇噪声功率谱，在有电晕噪声的时间段估计电晕噪声功率谱，最后将估计的风扇噪声以及电晕噪声功率谱重构时域信号，从而分离出风扇噪声及电晕噪声。本发明专利技术提供的方法，适用处理大型高压交流、直流变电站噪声信号。

【技术实现步骤摘要】
变电站噪声分离方法
本专利技术涉及音频处理技术，特别是涉及一种变电站噪声分离方法的技术。
技术介绍
在电力传输系统中，在相同额定功率条件下，输电电压越高，相应的输电电流越小，从而线路的能量损耗越小。因此，高压、超高压输电线路是目前以及未来电力建设的重点发展方向。在高压输电系统中，导线附近会产生超高电场引起附近空气电离，发生放电，空气放电会产生“嗞嗞”声音，形成电晕噪声。同时由于电压周期性变化，导线附近的带电粒子在电场作用下周期性运动，运动周期与交流电周期相同，在我国，交流电工频频率为50Hz。带电粒子与空气分子发生碰撞，形成声波向外辐射，这种周期性声压变化构成本体噪声，由于非线性效应，本体噪声还包括工频频率的各次谐波信号。高压输电变压器会产生大量热量，使得变电器温度升高，影响变压器使用寿命，一般高压变电器都带有冷却设施，主要是降温风扇，降温风扇会产生噪声，称为风扇噪声。在大型高压变电站内，由于变压器、高压传输线相对集中，变电站内的噪声会对周围居民生产生活产生不良影响，根据以上分析可知高压变电站内的噪声主要有电晕噪声、本体噪声以及风扇噪声三种。三种噪声分别具有不同的噪声来源，并且是同时存在与变电站内，目前还没有能实现对三种噪声进行单独测量的方法。
技术实现思路
针对上述现有技术中存在的缺陷，本专利技术所要解决的技术问题是提供一种能有效分离高压变电站的电晕噪声、本体噪声、风扇噪声，从而能实现对该三种噪声进行单独测量的变电站噪声分离方法。为了解决上述技术问题，本专利技术所提供的一种变电站噪声分离方法，其特征在于，具体步骤如下：1)利用传声器拾取高压变电站的噪声，高压变电站的噪声包括电晕噪声、本体噪声、风扇噪声，其信号模型为：s(n)＝c(n)+b(n)+e(n)w0＝2πf0/fs其中，s(n)为高压变电站噪声，c(n)为电晕噪声，b(n)为本体噪声，e(n)为风扇噪声，w0为数字域频率，f0为本体振动基频频率，f0＝50Hz，fs为信号采样频率，P为谐波个数；Bk为谐波幅值，φk为谐波相位，n为离散采样时间索引值；2)对s(n)实施多层小波变换，得到六组小波系数分别为：si(n),i＝1,2,...6；其中，s1(n)为低频小波系数，s1(n)对应的频带范围为[0,fs/64]，另五组小波系数si(n),i＝2,...6对应的频率范围分别为[fs/28-i,fs/27-i]；其中，小波变换中所采用的小波变换滤波器组是Daubechies滤波器组，小波系数采用Mallat快速算法计算；3)将s1(n)经过通带梳状滤波器之后做小波反变换，得到本体噪声的估计信号并将作为分离出来的本体噪声，通带梳状滤波器的通带频率是本体噪声的基频频率f0以及各次谐频；4)将s1(n)经过阻带梳状滤波器之后与另五组小波系数做小波反变换，得到包含电晕噪声及风扇噪声的混合信号y(n)；5)将s5(n)和s6(n)做小波反变换，得到电晕检测信号sd(n)；6)利用sd(n)检测有否有电晕噪声，并根据检测结果，在没有电晕噪声的时间段对风扇噪声功率谱进行估计，在有电晕噪声的时间段将y(n)的功率谱与估计的风扇噪声功率谱相减，得出估计的电晕噪声功率谱，功率谱计算采用快速傅里叶变换实现；6.1)利用sd(n)检测电晕噪声的步骤如下：6.1.1)对sd(n)、y(n)实施分帧，得到sd(n)、y(n)中各帧的帧信号为：Sq＝[sd(qL-L+1),sd(qL-L+2),...sd(qL)]Yq＝[y(qL-L+1),y(qL-L+2),...y(qL)]其中，Sq为sd(n)中第q帧的帧信号，Yq为y(n)中第q帧的帧信号，L为帧长；6.1.2)设定一个判决门限，将sd(n)、y(n)中各帧的最大幅值标准差与判决门限进行比较；如果有则判定y(n)中的第q帧有电晕噪声，反之则判定y(n)中的第q帧没有电晕噪声；其中，Th为判决门限，Th＝10，max{Sq}是指取Sq的最大值，E{sd(n)}是指取sd(n)的方差；6.2)风扇噪声功率谱的估计方法如下：如果则令反之则令如果则令反之则令如果y(n)中的第q帧有电晕噪声，且则反之则其中，y(n)中第q帧的功率谱为|Yq(k)|2，为y(n)中第q帧的最小谱值，中第q帧的最大谱值，为估计的风扇噪声第q帧的功率谱，α、β、γ、κ为可调参数；6.3)电晕噪声功率谱的估计方法如下：如果y(n)中的第q帧有电晕噪声，则令其中的|Cq(k)|2为y(n)中第q帧的电晕噪声功率谱；7)对估计的风扇噪声功率谱采用傅里叶反变换重构时域信号，并将重构后得到的信号作为分离出来的风扇噪声；对估计的电晕噪声功率谱采用傅里叶反变换重构时域信号，并将重构后得到的信号作为分离出来的电晕噪声；由于功率谱计算将y(n)分帧处理，信号重构也需要分帧重构，重构方法如下：在无电晕噪声发生时，风扇噪声帧信号为：在有电晕噪声发生时，风扇噪声及电晕噪声的帧信号分别为：其中，为风扇噪声重构信号中第q帧的帧信号，为电晕噪声重构信号中第q帧的帧信号，L为帧长，IFFT为傅里叶反变换。进一步的，步骤6.2)中，α∈[0.01,0.03]，β∈[0.4,0.65]，γ∈[4,6]，κ∈[0.1,0.3]。本专利技术提供的变电站噪声分离方法，利用电晕噪声、本体噪声、风扇噪声的时、频域特性，根据本体噪声的线谱特性，设计相应通带梳状滤波器实现本体噪声的分离，利用电晕噪声与风扇噪声在时—频域上的可分性，采用小波变换，实现电晕噪声与风扇噪声的分离，从而来达到分离这三种噪声的目的，能有效分离高压变电站的电晕噪声、本体噪声、风扇噪声，从而能实现对该三种噪声进行单独测量。附图说明图1是本专利技术实施例的变电站噪声分离方法的分离流程图；图2是直流电噪声的时域波形测量图；图3是直流电噪声在时域上放大后的时域波形测量图；图4是交流电噪声的时域波形测量图；图5是交流电噪声在时域上放大后的时域波形测量图。图6是采用本专利技术实施例的方法对直流电噪声进行分离后的分离效果图；图7是采用本专利技术实施例的方法对交流电噪声进行分离后的分离效果图。具体实施方式以下结合附图说明对本专利技术的实施例作进一步详细描述，但本实施例并不用于限制本专利技术，凡是采用本专利技术的相似结构及其相似变化，均应列入本专利技术的保护范围。如图1所示，本专利技术实施例所提供的一种变电站噪声分离方法，其特征在于，具体步骤如下：1)利用传声器，以65536Hz的采样率拾取高压变电站的噪声，高压变电站的噪声包括电晕噪声、本体噪声、风扇噪声，其信号模型为：s(n)＝c(n)+b(n)+e(n)w0＝2πf0/fs其中，s(n)为高压变电站噪声，c(n)为电晕噪声，b(n)为本体噪声，e(n)为风扇噪声，w0为数字域频率，f0为本体振动基频频率，f0＝50Hz，fs为信号采样频率，P为谐波个数，Bk为谐波幅值，φk为谐波相位，n为离散采样数据时间索引值，一般在100Hz、200Hz、300Hz等谐波频率上本体噪声能量较大，电晕噪声c(n)表现为时间上随机出现的脉冲信号，传输电压越高，电晕噪声越密集，同时脉冲时间间隔要远大于脉冲持续时间，风扇噪声e(n)为宽带噪声；2)对s(n)实施多层小波变换(本例中为五层)，得到六组小波系数分别为：si(n),i＝1,2,...本文档来自技高网...

【技术保护点】
一种变电站噪声分离方法，其特征在于，具体步骤如下：1)利用传声器拾取高压变电站的噪声，高压变电站的噪声包括电晕噪声、本体噪声、风扇噪声，其信号模型为：s(n)＝c(n)+b(n)+e(n)b(n)=Σk=1PBkcos(kw0n+φk)]]>w0＝2πf0/fs其中，s(n)为高压变电站噪声，c(n)为电晕噪声，b(n)为本体噪声，e(n)为风扇噪声，w0为数字域频率，f0为本体振动基频频率，f0＝50Hz，fs为信号采样频率，P为谐波个数；Bk为谐波幅值，φk为谐波相位，n为离散采样时间索引值；2)对s(n)实施多层小波变换，得到六组小波系数分别为：si(n),i＝1,2,...6；其中，s1(n)为低频小波系数，s1(n)对应的频带范围为[0,fs/64]，另五组小波系数si(n),i＝2,...6对应的频率范围分别为[fs/28‑i,fs/27‑i]；其中，小波变换中所采用的小波变换滤波器组是Daubechies滤波器组，小波系数采用Mallat快速算法计算；3)将s1(n)经过通带梳状滤波器之后做小波反变换，得到本体噪声的估计信号并将作为分离出来的本体噪声，通带梳状滤波器的通带频率是本体噪声的基频频率f0以及各次谐频；4)将s1(n)经过阻带梳状滤波器之后与另五组小波系数做小波反变换，得到包含电晕噪声及风扇噪声的混合信号y(n)；5)将s5(n)和s6(n)做小波反变换，得到电晕检测信号sd(n)；6)利用sd(n)检测有否有电晕噪声，并根据检测结果，在没有电晕噪声的时间段对风扇噪声功率谱进行估计，在有电晕噪声的时间段将y(n)的功率谱与估计的风扇噪声功率谱相减，得出估计的电晕噪声功率谱，功率谱计算采用快速傅里叶变换实现；6.1)利用sd(n)检测电晕噪声的步骤如下：6.1.1)对sd(n)、y(n)实施分帧，得到sd(n)、y(n)中各帧的帧信号为：Sq＝[sd(qL‑L+1),sd(qL‑L+2),...sd(qL)]Yq＝[y(qL‑L+1),y(qL‑L+2),...y(qL)]其中，Sq为sd(n)中第q帧的帧信号，Yq为y(n)中第q帧的帧信号，L为帧长；6.1.2)设定一个判决门限，将sd(n)、y(n)中各帧的最大幅值标准差与判决门限进行比较；如果有则判定y(n)中的第q帧有电晕噪声，反之则判定y(n)中的第q帧没有电晕噪声；其中，Th为判决门限，Th＝10，max{Sq}是指取Sq的最大值，σsd2=E{(sd(n)-E{sd(n)})2},]]>E{sd(n)}是指取sd(n)的方差；6.2)风扇噪声功率谱的估计方法如下：如果则令反之则令如果则令反之则令如果y(n)中的第q帧有电晕噪声，且则E^q(k)=E^q-1(k),]]>反之则E^q(k)=E^q-1(k)+κ·(|Yq(k)|-E^q-1(k));]]>其中，y(n)中第q帧的功率谱为|Yq(k)|2，为y(n)中第q帧的最小谱值，为y(n)中第q帧的最大谱值，为估计的风扇噪声第q帧的功率谱，α、β、γ、κ为可调参数；6.3)电晕噪声功率谱的估计方法如下：如果y(n)中的第q帧有电晕噪声，则令其中的|Cq(k)|2为y(n)中第q帧的电晕噪声功率谱；7)对估计的风扇噪声功率谱采用傅里叶反变换重构时域信号，并将重构后得到的信号作为分离出来的风扇噪声；对估计的电晕噪声功率谱采用傅里叶反变换重构时域信号，并将重构后得到的信号作为分离出来的电晕噪声；由于功率谱计算将y(n)分帧处理，信号重构也需要分帧重构，重构方法如下：在无电晕噪声发生时，风扇噪声帧信号为：E^q=[e^(qL-L+1),e^(qL-L+2),...e^(qL)]=Yq]]>在有电晕噪声发生时，风扇噪声及电晕噪声的帧信号分别为：E^q=IFFT{|E^q(k)|ej∠Yq(k)}C^q=IFFT{|C^q(k)|ej∠Yq(k)}]]>其中，为风扇噪声重构信号中第q帧的帧信号，为电晕噪声重构信号中第q帧的帧信号，L为帧长，IFFT为傅里叶反变换。...

【技术特征摘要】
1.一种变电站噪声分离方法，其特征在于，具体步骤如下：1)利用传声器拾取高压变电站的噪声，高压变电站的噪声包括电晕噪声、本体噪声、风扇噪声，其信号模型为：s(n)＝c(n)+b(n)+e(n)w0＝2πf0/fs其中，s(n)为高压变电站噪声，c(n)为电晕噪声，b(n)为本体噪声，e(n)为风扇噪声，w0为数字域频率，f0为本体振动基频频率，f0＝50Hz，fs为信号采样频率，P为谐波个数；Bk为谐波幅值，φk为谐波相位，n为离散采样时间索引值；2)对s(n)实施多层小波变换，得到六组小波系数分别为：si(n),i＝1,2,...6；其中，s1(n)为低频小波系数，s1(n)对应的频带范围为[0,fs/64]，另五组小波系数si(n),i＝2,...6对应的频率范围分别为[fs/28-i,fs/27-i]；其中，小波变换中所采用的小波变换滤波器组是Daubechies滤波器组，小波系数采用Mallat快速算法计算；3)将s1(n)经过通带梳状滤波器之后做小波反变换，得到本体噪声的估计信号并将作为分离出来的本体噪声，通带梳状滤波器的通带频率是本体噪声的基频频率f0以及各次谐频；4)将s1(n)经过阻带梳状滤波器之后与另五组小波系数做小波反变换，得到包含电晕噪声及风扇噪声的混合信号y(n)；5)将s5(n)和s6(n)做小波反变换，得到电晕检测信号sd(n)；6)利用sd(n)检测有否有电晕噪声，并根据检测结果，在没有电晕噪声的时间段对风扇噪声功率谱进行估计，在有电晕噪声的时间段将y(n)的功率谱与估计的风扇噪声功率谱相减，得出估计的电晕噪声功率谱，功率谱计算采用快速傅里叶变换实现；6.1)利用sd(n)检测电晕噪声的步骤如下：6.1.1)对sd(n)、y(n)实施分帧，得到sd(n)、y(n)中各帧的帧信号为：Sq＝[sd(qL-L+1),sd(qL-L+2),...sd(qL)]Yq＝[y(qL-L+1),y(qL-L+2),...y(qL)]其中，Sq为sd(n)中第q帧的帧信号，Yq为y(n)中第q帧的帧信号，L为帧长；6.1.2...

【专利技术属性】
技术研发人员：彭任华，郑成诗，杨鹤飞，李晓东，
申请(专利权)人：中国科学院上海高等研究院，
类型：发明
国别省市：上海;31