[发明专利]一种抑制管道堵塞声信号中宽带背景噪声的方法

权利要求书

1.一种抑制管道堵塞声信号中宽带背景噪声的方法，包括对输入的含平稳宽带背景噪声的油气管道堵塞低频回波信号提取其频域谱、对背景噪声的频域谱进行估计、利用维纳滤波器抑制背景噪声的步骤，其特征在于，

油气管道堵塞低频回波信号的频域谱提取方法如下：

1)首先定义油气管道堵塞回波信号模型：

y(n)＝s(n)+b(n)            (1)

y(n)是数据采集系统采集到的信号，s(n)是管道真实堵塞回波信号，b(n)是背景噪声信号；

2)将输入信号y(n)分成长度为含有L个互不重叠样点的时间窗，对每个时间窗中的L个样点，用傅里叶变换提取频域谱：

Y(k)=Σn=0N-1y(i)·e-j2πNk·n---(2)]]>

其中k＝1，2，...，N-1，y(i)表示主传感器接收到的时域信号，Y(k)就是当前输入信号的频域谱，也称为傅里叶变换谱；

背景噪声的频域谱估计方法如下：

1)把式(1)变换到频域谱：

Y(k)＝S(k)+B(k)    (3)

Y(k)是采集到信号的频域谱，S(k)是目标信号频域谱，B(k)是背景噪声频域谱，

2)系统背景噪声频域谱初始化：

B(k0)=Σk=1Nα(k)B(k)---(4)]]>

3)背景噪声频域谱估计：

B(kn)=λB(kn-1)+(1-λ)Y(kn)ifY(kn)<βB(kn-1)B(kn-1)others---(5)]]>

其中λ是背景噪声频域谱的更新因子，β是更新控制因子，用来控制堵塞回波信号频域谱和背景噪声频域谱的强弱比率；

利用维纳滤波器抑制背景噪声的方法如下：

1)维纳滤波器参数估计：

在完成背景噪声频域谱估计的情况下，用下式获取维纳滤波系数：

H(k)=Y(k)-B(k)Y(k)---(6)]]>

2)获取维纳滤波：

S^(k)=Y(k)H(k)---(7)]]>

3)获取消除背景噪声后的堵塞反射回波时域信号：

将式(6)中做傅里叶逆变换，得

S^(n)=Σk=0N-1S^(k)·ej2πNk·n---(8)]]>

即是堵塞低频反射回波信号的时域表示。