阳耀 基于小波降噪的铁路边坡测斜监测分析 刘 伟，朱玉婷，付平勇 （华东交通大学道路与铁道工程重点实验室，江西南昌３３００１ ３） 【摘要】 铁路边坡测斜监测数据中时常含有许多锯齿状波动，为了消除或减弱波动的干扰，通常对监 测数据进行降噪处理。基于小波降噪理论，利用ＭａｔｌａｂＲ２００７ａ软件，分别选用小波变换模极大值法、非线性 小波变换软阈值法与非线性小波变换硬阈值法三种方法对峰福铁路Ｋ２９４处边坡测斜监测数据进行降噪处 理，结果表明：以上三种方法降噪效果明显，其中非线性小波变换软阈值法更有效，更适合于铁路边坡测斜监 测分析。 【关键词】 小波分析； 阈值处理；数据降噪；测斜监测；Ｍａｔｌａｂ软件 【中图分类号】Ｕ２１３．１ ３ 铁路边坡的深部水平侧移曲线时常锯齿状波动，一方面 【文献标识码】Ａ 细节；｛ｈ ｝ 与｛ｇ ｝ 是一对离散正交镜像滤波器（低通 和高通），（４）式即为Ｍａｌｌａｔ塔式分解算法。 Ｍａｌｌａｔ重构算法： ｅ ＝ ｅ川＋Ｇ Ｄ ＋ｌ， （５） 反映降雨、地震等自然环境因素或爆破、钻孔等人为施工因 素作用使边坡内部位移产生突变，另一方面，反映监测信号 由于电压不稳、电磁干扰等不良因素产生了噪声。为了提高 铁路边坡稳定性监测与评价的准确性，有必要对采集的监测 信号进行降噪处理 ，２ ７。小波分析作为一种有力的数学工 具，具有多分辨分析的特点，在时频和频域都具有表征信号 局部特征的能力，时间窗和频率窗都可以根据信号的具体形 态动态调整，基于这些特征，小波分析可以广泛应用于监测 其中， ＝Ｊ，Ｊ一１，…，１，Ｈ 和Ｇ ，分别为日和Ｇ的 共轭。 小波分解与重构降噪法就是根据Ｍａｌｌａｔ分解式（１）至式 （３），将含噪声的信号分解在不同的频带内，然后，在信号重 构时，将与噪声相应的高频细节信号有关部分Ｄ，置零，重构 信号便可达到降噪的目的。 信号降噪 Ｊ。 ｌ小波降噪基本原理 利用小波分析降噪 卜 ，即在不同尺度下作小波变 换，其实质是用不同中心频率的带通滤波器对信号进行滤 ）＝Ａ ）＝∑ｃｍ ｊ，ｋ（ ）。 ｋＥＺ （６） ）是 ）经过滤波后的平滑信号表达式，即为提取的 真实信号，这便是小波降噪的基本原理。因此信号降噪可分 三步处理：首先对原始信号进行小波分解，则噪声通常含在 波，以取得较好质量的信号。在实际的工程中，有用信号通 常表现为比较平稳的信号或低频信号，而噪声信号通常表现 为高频信号。 设 ）是观测点观测序列，采用Ｍａｌｌａｔ算法，将信号分 解成不同频率成分： ，（ ）＝Ａ，一 ）＝Ａ ）＋Ｄ ） （１） 高频系数中；然后对小波分解的高频系数以阈值等形式进行 量化处理；最后再对信号重构即可达到降噪的目的。小波降 噪的方法很多，如Ｍａｌｌａｔ算法、小波变换模极大值法、非线性 小波变换阈值法和平移不变量小波降噪法等，而小波变换模 极大值法降噪和非线性小波变换阈值法降噪最为常用。 １．１小波变换模极大值法 其中，信号在低频空间 上的投影为： Ａ ） ∑Ｃｊ， 咖 （ ） ｋＥＺ （２） 由于信号和噪声的小波变换系数在不同尺度上具有不 同的传播特性，即随着尺度的增大，噪声所对应的模极大值 迅速衰减。而信号的模极大值有三情形：若为缓变信号，则 模极大值逐渐增大；若为阶越信号，则模极大值保持不变；若 为脉冲信号，正、负极值组成的脉冲对的幅值将同时变小。 于是，连续进行若干次小波分解后，综合各个尺度上模极大 值的位置和幅值，可断定哪些模极大值是由噪声引起，哪些 在高频空间 ，上的投影为： Ｄ ） ∑Ｄｍ （ ） ￡Ｚ （３） Ａ ）是信号／（ ）的频率不超过２ 的成分，Ｄ （ ）是 之间的成分。 Ｃ ＋ｌ ＨＣ』 频率介于２ 与２ 上述小波分解式实际可以写成如下矩阵形式： ［定稿日期］２０１１—０５—１１ Ｄ川＝∞ （４） ［基金项目］江西省自然科学基金项目（２０１０ＧＺＣ００２３）； 江西省教育厅项目（ＧＪＪ１１４４１） 式中： ＝１，２，…－，；Ｈ＝（ｈ ），Ｈ为尺度函数对应的低 通滤波器；Ｇ＝（ｇ ），Ｇ为小波函数对应的带通滤波器；Ｃ， 为在２一 分辨率下的离散逼近；Ｄ，为在２ 分辨率下的离散 ［作者简介］刘伟（１９８７～），男，硕士在读，研究方向为 地下工程稳定性分析与结构优化。 １２４ 四川建筑第３１卷６期２０１ １．１２ ・ 蛰 耩龋 是由信号产生。去除噪声引起的模极大值，再重构剩余的模 极大值，便可实现降噪。 １．２非线性小波变换阈值法 ～ｊ （１）边坡测斜布置 该测区属剥蚀低丘，边坡为土质边坡，山体自然坡度２０。 ３０。，植被发育，线路沿水Ｅｌ库区北岸行走，整个病害长度 约为３００ ｍ，包括一段路堑及其两端各一段路堤。根据路段 地形地貌及地质情况，布置２排共６个测斜孔，其中一排布 置在左侧路堑坡顶以上的山坡上，且位于已有纵向裂缝的下 非线性小波变换阈值法也称小波收缩。当信号的Ｌｉｐｓ— ｃｈｉｔｚ指数为正时，其小波变换极大值的幅度将随着尺度而增 加；但噪声的Ｌｉｐｓｃｈｉｔｚ指数为负时，则噪声所对应的小波变 换模极大值将随尺度的增加而减少。利用这个特性，在不同 方，测孔间距５０ ｍ，布置有ＣＸ一２、ＣＸ一６、ＣＸ一１３与ＣＸ一７ 四个测孑Ｌ；另一排布置在路堑右侧水沟外侧，布置有ＣＸ一１ 和ＣＸ一９两个测孔，布置见图１所示。 的分解尺度上设置一定的阈值，将小于该阈值的模极大值点 认为是噪声的小波变换引起的，将其置为零；大于该阈值的 模极大值点则认为是信号小波变换引起的而予以保留。常 见的阈值规则有ｓｑｔｗｏｌｏｇ规则、ｒｉｇｒｓｕｒｅ规则、ｈｅｕｒｓｕｒｅ规则和 ｍｉｎｉｍａｘｉ规则。工程中常用的简单阈值处理法有非线性小 波变换软阈值处理法、非线性小波变换硬阈值处理法等。 １．３ 小波降噪的实现与性能评价 （１）为了选择较优的小波降噪方案，将涉及到小波函数、 阈值的选择。实际工程中可根据自相似原则选取小波基函 数，若选择的小波基函数与待处理信号有相似性，小波变换 后，利于信号能量集中和减少计算量。土木工程 中关于测 斜监测曲线，用Ｓｙｍｌｅｔ、Ｃｏｆｌｅｔ、Ｄａｕｂｅｃｈｉｅｓ小波系均可实现曲 线的降噪平滑，但以Ｄａｕｂｅｃｈｉｅｓ小波系效果最优。对于阈值 图１ 测斜孔位布置示意图 的选择，其实现十分简单，一个简单的方法就是用每个小波 （２）数据采集与分析 待测斜管安装完，使用测斜仪进行预通并采集初读数， 再按时进行其他测斜监测。下面列出峰福铁路Ｋ２９４边坡 ＣＸ一６测斜孔采集的部分水平位移监测数据，并用小波变换 模极大值法、非线性小波变换软阈值法与非线性小波变换硬 阈值法三种方法对所采集数据进行降噪处理，降噪结果见 表Ｉ。 表１ ＣＸ一６测斜监测降噪处理结果 分解系数的幅值减去阈值，若差为负值则将系数设置为０；若 差是正值则保持系数不变。在小波降噪方案实施过程中，同 时也可以采用Ｍａｔｌａｂ中的小波工具实现。Ｍａｔｌａｂ中获取阈 值的函数有：ｄｄｅｎｃｍｐ、ｔｈｓｅｌｅｃｔ、ｗｂｍｐｅｎ、ｗｄｃｂｍ；实现信号阈 值去噪的函数有：ｗｄｅｎ、ｗｄｅｎｃｍｐ、ｗｔｈｒｅｓｈ、ｗｔｈｃｏｅｆ、ｗｐｔｈｃｏｅｆ、 ｗｐｄｅｎｃｍｐ。 （２）为了对各种降噪方法进行比选，可用小波降噪的性 能评价标准来评判，工程中常用的小波降噪的性能评价标准 有信噪比（ＳＮＲ）、信噪比增益（Ａ）和均方根误差（ＲＳＭＥ），这 三种评价标准是等价的 Ｊ。其中信噪比是测量信号中噪声 量度的传统方法，常被用来作为评价降噪效果的指标，信噪 比越高则滤波效果越好，则降噪性能越好。信噪比单位是分 贝，其定义为： 对应 采集数 小波变换模 非线性小波变换阈值法 监测日期 极大值法 降噪（ｈｅｒｓｕｒｅ阈值）（ｍｍ） 周数 据（ｍｍ） 降噪（ｍｍ） 软阈值法 硬阔值法 ２０１０．１２．Ｏ５ １ —７．３１ —８．０１３４ —７．８０３８ —８．１７ｌ２ ２０１０．１２．１２ ２ —３．２４ ３．１５ Ｏ．８７ —１．７３６５ ２．０３６１ ２．０７１４ —２．６８４５ ２．７３８２ ２．ＯＯ４１ —３．１４８６ ２．８７４４ ２．７７４２ ＳＮＲ＝１０ｌｏｇ ｌ０ｌ — ————一ｆ ÷∑ ｎ） ］ ｌ２０１０．１２．１９ ３ （７） ｌ ２工程实例 ｎ）一 ｎ）］ Ｊ ２０１０．１２．２６ ４ 式中｜厂（／７，）是原始信号 ｎ）为小波降噪后的估计信号。 ２０１１．０１．０２ ５ ２０１１．Ｏ１．Ｏ９ ６ ２０１１．０１．１６ ７ ２０１１．Ｏ１．２３ ８ ６．４４ １Ｏ．Ｏ６ ｌ０．３７ ８．１４ ５．９４３７ ８．７４５２ ９．７３６２ ９．３４０２ ６．２５４４ ９．２２９５ １０．２７０３ ８．９１６４ ６．３１２９ ９．３９６８ ９．７３３４ ９．８９６７ 以峰福铁路Ｋ２９４＋６８０一Ｋ２９４＋９００段边坡ＣＸ一６测 斜孑Ｌ，从２０１０年１２月至２０１１年３月（测量间隔为１周，共间 隔１６周）的水平位移监测数据作为小波降噪处理的对象，监 测数据均通过某高科技实业有限公司ＣＸ一８０６Ｂ型数字显 示滑动式钻孔测斜仪采集获得，并使用小波变换模极大值 法、非线性小波变换软阈值法和非线性小波变换硬阈值法对 所采集数据进行降噪处理，小波函数选择Ｄｂ３，尺度选取３， ２０１１．Ｏ１．３０ ９ １１．４２ １０．５４３４ １１．５６８７ ｌ０．６５ｌ８ ９．Ｏｌ７５ １２．８５４９ １４．４９３８ ｌ５．３７９６ １０．３０９２ １３．５８９１ １０．８９６５ ９．３７２１ １３．１５６４ １３．３９４５ １４．２１３４ １１．５１３２ １２．６３５５ １０．６６３６ １Ｏ．１１９７ １４．１１５１ １３．７５８３ １５３５８５ ．２Ｏ１１．０２．Ｏ６ １０ １２．８７ ２０１１．０２．１３ １１ ２０１１．Ｏ２．２０ １２ ２０１１．Ｏ２．２７ ｌ３ ２０１１．０３．Ｏ６ １４ ２０１１．０３．１３ １５ ９．５５ ７．８６ １４．５３ ｌ３．７ １４．１３ 阈值选取常用的ｈｅｒｓｕｒｅ阈值，并用信噪比对降噪效果作对 比，提取铁路边坡体内部水平变形的信息，进行铁路边坡测 斜监测分析。 四川建筑第３１卷６期２０１ １．１２ １２５ 续表１ 非线性小波变换硬阈值法三种方法降噪效果都较明显，利于 测斜监测数据反映铁路边坡的位移变化趋势。 （２）基于ｈｅｒｓｕｒｅ阈值的非线性小波变换软阈值法，降噪 后的曲线有更高的信噪比，有较好的光滑性，且有效地保留 了测斜水平位移监测值的原有真实特性，故可以直观地获取 产生水平位移突变所在深度，有利于铁路边坡变形的判定和 稳定性的评价，更适合铁路边坡测斜监测分析。 （３）利用小波降噪理论，结合Ｍａｔｌａｂ处理工具箱，在工 程中处理深层水平位移测斜监测可节约测量人员内业处理 对应 采集数 小波变换模 非线性小波变换阈值法 监测日期 极大值法 降噪（ｈｅｒｓｕｒｅ阈值）（ｍｍ） 周数 瞎（ｍｍ） 降噪（ｍｍ） 软阈值法 硬阈值法 ２０１１．０３．２０ １６ １８．７３ １８．４５２６ １７．７６７２ １８．５ｌ８９ ２０１１．０３．２７ １７ ２０．９４ ２１．３１０２ ２１．１９１９ ２０．９７６５ 信噪比（ｄＢ） ， ７ ２０．５８２２ ２２．６２４９ ２１．２９３６ ＣＸ一６孔在第３、４、８、１０与ｌ２周采集的原始数据有明 显的波动，其中第４、１０与１２周表现较剧烈。通过降噪前后 的数据结果对比，小波变换模极大值法、非线性小波变换软 分析时间，也使降噪后的数据更具科学性与说服性，值得在 相关领域推广。 参考文献 闯值法与非线性小波变换硬阈值法三种方法降噪效果都较 明显，降噪后的数据穿行在原始数据之间。第４与ｌ２周时， 小波变换模极大值法与非线性小波变换软阈值法降噪后的 曲线，较非线性小波变换硬阈值法更好地保留了突变波动。 （３）边坡测斜降噪评价 ［１］闾振，陆阳，田斌鹏．小波分析在边坡深部水平位移监测数据 降噪中的应用［Ｊ］．中南公路工程，２００７，３２（２）：４１—４４ ［２］黄飞澜，肖红．测斜仪在高填方地基侧向水平位移监测中的应 由表１可见，非线性小波变换阈值法降噪的信噪比大于 小波变换模极大值法；而在非线性小波变换阈值法中，又以 非线性小波变换软阈值法的信噪比最高。由此可见，由于非 线性小波变换软阈值法和小波变换模极大值法都建立在奇 异性分析基础上，故降噪后的曲线对原始曲线中的尖峰和突 变进行了有效地保持。但从曲线的平滑程度上看，小波变换 模极大值法的平滑性优于非线性小波变换阈值法。同理，可 用 Ｉ｛． 轻 褪２Ｏ１０，３５（５）：１１２—１２０ ［３］冯毅，王香华．小波变换降噪处理及其Ｍａｔｌａｂ实现［Ｊ］．数据采 集与处理，２００６，２１（１２）：３７—３９ ［４］ 张鹏，李献勇，陈建平．基于小波降噪的隧道围岩监测数据分 析［Ｊ］．吉林大学学报（地球科学版），２００８，３８（６）：１０１０—１０１４ ［５］张鹏．基于小波理论的白鹤隧道围岩稳定性非线性研究［Ｄ］． 长春：吉林大学，２００９ ［６］ 田胜利，周拥军，葛修润，等．基于小波分解的建筑物变形监测 对Ｋ２９４处边坡其他测斜孔的监测数据进行降噪。通过对其 他各测孑Ｌ采集得较多数据降噪处理效果对比验证，基于ｈｅｒ． ｓｕｒｅ阈值的非线性小波变换软阈值法更合理有效，能够敏感 识别监测数据噪声和有用水平位移信息，更适合铁路边坡测 斜监测分析。 数据处理［Ｊ］．岩石力学与工程学报，２００４，２３（１５）：２６３９ ２６４２ ［７ Ｊ Ｍａｌｈ Ｓ，Ｈｗａｎｇ Ｗ Ｌ．Ｓｉｎｇｕｌａｒｉｔｙ Ｄｅｔｅｃｔｉｏｎ ａｎｄ Ｐｒｏｃｅｓｓｉｎｇ ｗｉｔｈ Ｗａｖｅｌｅｔｓ［Ｊ］．ＩＥＥＥ Ｔｒａｎｓ Ｉｎｆｏｒｍａｔｉｏｎ Ｔｈｅｏｒｙ，１９９２，３８（２）：６１７ —６４３ ３结论 ［８］Ｄｏｎｏｈｏ ＤＬ．Ｄｅ—ｎｏｉｓｉｎｇ ｂｙ Ｓｏｆｔ—ｔｈｒｅｓｈｏｌｄｉｎｇ［Ｊ］．ＩＥＥＥ Ｔｒａｎｓａｃ— ｔｉｏｎｓ ｏｎ Ｉｎｆｏｒｍａｔｉｏｎ Ｔｈｅｒｏｙ，１９９５，４１（３）：６１３—６２７ （１）小波变换模极大值法、非线性小波变换软阈值法与 １２６ 四川建筑第３１卷６期２０１ １．１２
