[发明专利]水锤冲击压力信号的自适应滤波降噪方法

说明书

本发明公开了一种水锤冲击压力信号的自适应滤波降噪方法，包括获取水锤冲击压力信号，对其进行信号两端的扩展，得到扩展后压力信号构成的时间序列为x(t)，对x(t)进行经验模式分解，将其转化为一系列水锤信号IMF分量，剔除IMF分量中的一部分高阶分量，将剩余低阶分量进行暴露，并对其扩展部分进行截取，得到截取分量，组合全部截取分量得到重构信号，评估重构信号的最大幅值损失、能量损失和R2系数，判断滤波后信号是否符合要求，通过评估后，将重构信号作为滤波后信号，否则重新剔除高阶分量，重新进行评估。本发明在滤波过程中，无需设定滑动步长、滤波器参数等，极大地减少了人为因素的影响，具有很强的泛化能力和通用性。

技术领域

本发明属于长距离输水系统的监测及安全管理技术领域，涉及一种水锤冲击压力信号的自适应滤波降噪方法。

背景技术

由于我国地域辽阔且地势复杂，导致地区间水资源分布情况差异巨大，随着大型输水工程的分阶段实施，我国已建成大批长距离有压输水管道。在密闭管路系统内，由于流体流量急剧变化而引起较大的压力波动并造成振动的现象称之为水锤效应，其瞬间压力可大大超过正常压力，并经常产生破坏性影响，水锤效应严重影响着水利设施的安全和质量，为了掌握水锤发生机理并对其开展抑制预防工作，国内外学者从产生机理、信号特征、模式识别等多个领域对其开展了研究分析，相关研究对于水利系统的安全运行有着重要意义。

水锤冲击是一种具有较高能量的瞬时冲击，具有典型的简谐波特性，其压力信号具有其自身特殊性：一方面，水锤信号主要由低频成分组成，其能量大部分集中于低频成分，高频信息主要以噪声为主，这些噪声多是由于管路、电流、电路等因素导致，这区别于机械振动等高频信号；同时，水锤冲击信号从形态上具有缓变特性，由于其能量在流体中分布呈现明显的集中特性，其主成分的周期性非常明显，这区别于其它类别的压力信号。

在实际测量中，通常通过测量管路内部压力获取水锤信号，这类信号数据是支撑各类研究和监测预警应用的基础，有效地获取水锤压力信号具有重要意义。压力信号是一类特殊的低频信号，由于压力信号有效地表征了设备内部的受力状态，在各种水力机械和输水设施中应用广泛，是水轮机、输水管道、泵站等设施的基本监测信号类型。实际中由于环境干扰、传感器误差、测试电流干扰、输水设施振动等因素，会导致水锤冲击信号中耦合大量噪声，因而需要对其进行有效地降噪滤波，但是在通过各类滤波器对信号进行降噪处理过程中，往往存在如下问题：首先，由于不同工况的压力信号存在很大差异，滤波器的参数难以确定，例如滤波步长、阶数、滤波器类型等；同时，常见滤波方法在很大程度上会造成能量的损失，导致信号幅值和能量的下降，例如滑移平均滤波、五点三次滤波、Savitzky-Golay滤波等方法，在水锤冲击工况中，这种损失可能会导致严重的偏差，埋下一定的安全隐患；此外，由于输水工程管线复杂、工况类型繁多，不同工况下相关工程参数完全不同，仅凭经验性难免耗时费力。

发明内容

本发明的目的是提供一种水锤冲击压力信号的自适应滤波降噪方法，解决了现有水锤冲击压力信号滤波降噪方法误差较大、能量损失较大的问题。

本发明所采用的技术方案是，水锤冲击压力信号的自适应滤波降噪方法，包括以下步骤：

步骤1，获取待监测水管的水锤冲击压力信号构成的时间序列为 x(t)；

步骤2，对x(t)进行经验模式分解，确定极大值序列x_max和极小值序列x_min，根据x_max和x_min确定x(t)的上下包络及两组包络的平均值序列m(t)＝(x_max+x_min)/2，将m(t)从x(t)中去掉，得到第一个筛分序列h₁(t)＝x(t)-m(t)；