欒元重 欒亨宣 馬德鵬 劉 娜 周 巖

1)山東科技大學(xué)測繪科學(xué)與工程學(xué)院，青島 266590

2)山東科技大學(xué)機械電子工程學(xué)院，青島 266590

3)山東科技大學(xué)資源與環(huán)境工程學(xué)院，青島 266590

4)重慶地礦測繪院，重慶 400042

橋梁變形數(shù)據(jù)小波去噪與混沌預(yù)測*

欒元重1)欒亨宣2)馬德鵬3)劉 娜4)周 巖1)

1)山東科技大學(xué)測繪科學(xué)與工程學(xué)院，青島 266590

2)山東科技大學(xué)機械電子工程學(xué)院，青島 266590

3)山東科技大學(xué)資源與環(huán)境工程學(xué)院，青島 266590

4)重慶地礦測繪院，重慶 400042

跨海大橋因海浪、颶風(fēng)、船體撞擊等外界因素作用，橋梁變形表現(xiàn)為動態(tài)非線性變形狀態(tài)，也存在混沌現(xiàn)象。為此，選擇db10小波和軟閾值消噪法分解與消除變形數(shù)據(jù)中的粗差與噪聲，并對小波預(yù)處理后的平滑變形時間序列通過相空間重構(gòu)，運用Lyapunov指數(shù)進行橋梁混沌狀態(tài)判定;再將混沌時間序列預(yù)測結(jié)果與指數(shù)平滑法預(yù)測值進行分析對比，證明混沌預(yù)測可靠且精度高。

軟閾值消噪法;相空間重構(gòu);Lyapunov指數(shù);混沌時間序列;指數(shù)平滑法

1 引言

近些年，我國已有多個跨海大橋建成運營。橋梁處在海洋環(huán)境中，受海浪、颶風(fēng)、船體撞擊等因素的影響，使橋梁變形噪聲顯著，產(chǎn)生撓曲變形、垂直位移及水平移動，表現(xiàn)為動態(tài)非連續(xù)變形，還有突變、混沌現(xiàn)象發(fā)生。如何將橋梁變形數(shù)據(jù)噪聲除去，顯現(xiàn)出橋梁變形的特征規(guī)律;如何判定橋梁變形混沌狀態(tài)，建立橋梁變形的預(yù)測模型;是橋梁變形數(shù)據(jù)處理的研究熱點。

橋梁變形過程隨機時間空間變化的信號，可利用小波分析對變形信息進行時頻分析、奇異性探測、提取空間變形趨勢項和粗差及變形信號多尺度分析。當(dāng)選用合適的小波函數(shù)，經(jīng)去噪處理，能消除變形數(shù)據(jù)中的高頻噪聲，提高變形數(shù)據(jù)預(yù)測精度。

經(jīng)小波處理后的變形數(shù)據(jù)，首先計算最大Lyapunov指數(shù)并判別橋梁變形系統(tǒng)的混沌特征;然后選取合適的延遲時間τ和嵌入維數(shù)m，進行相空間重構(gòu);最后對混沌時間序列進行預(yù)測，實現(xiàn)橋梁變形預(yù)測創(chuàng)新。本文將混沌時間序列預(yù)測結(jié)果與指數(shù)平滑法預(yù)測值進行了分析對比，證明了混沌時間序列預(yù)測精度值較高。

2 橋梁變形數(shù)據(jù)小波分解與重構(gòu)去噪

2.1 小波基函數(shù)選擇

對 Haar、db6、db10、Mexican Hat小波［1］分別加入高斯白噪聲的信噪比為3的blocks信號，并進行分解和低頻重構(gòu)去噪實驗，db10小波去噪后與原測試信號相比，均方根值最小，信噪比最高，平滑度指數(shù)也最低，故本文采用db10小波函數(shù)進行橋梁變形數(shù)據(jù)預(yù)處理。

2.2 小波分解與重構(gòu)去噪

將橋梁變形數(shù)據(jù)f(x)分解成不同頻率成分，信號在空間Vj上投影為:

在空間Wj上的投影為:

式中Ajf(x)是信號f(x)頻率不超過2－j成分，Djf(x)是頻率介于2－j與2－j+1之間的成分。

矩陣形式為:

式中，H為低通濾波器，H=(hk－2n);G為帶通濾波器。Mallat重構(gòu)算法為:

f(x)經(jīng)濾波后平滑表達式(x)為:

3 橋梁變形混沌特征判定

3.1 橋梁變形倍周期分叉混沌性質(zhì)

當(dāng)橋梁系統(tǒng)有擾動時，系統(tǒng)發(fā)生分岔并可能演變成混沌狀態(tài)。取Verhulst生物模型［2］擬合橋梁變形位移時間序列，其差分方程為:

3.2 Lyapunov指數(shù)

3.3 橋梁系統(tǒng)混沌特征判定

Lyapunov指數(shù)是用來衡量軌道收斂率或發(fā)散率。Lyapunov指數(shù)表示多次迭代中平均每次迭代引起指數(shù)分離中的度量。Lyapunov指數(shù)λ值表達式為:

進行橋梁變形混沌預(yù)測，首次應(yīng)判定橋梁變形時間序列中是否具有混沌特征，表(1)列出了不同運動系統(tǒng)混沌特征量的判別標(biāo)準(zhǔn)。

表1 各不同運動特征值取值Tab.1 Characteristic value of different movements

4 橋梁垂直變形觀測數(shù)據(jù)小波去噪

取青島跨海大橋2011年8月—2012年5月垂直變形觀測數(shù)據(jù)，A1點下沉速度如圖1所示。

從圖1可以看到，沉降量的變化是不均勻的。在橋梁變形觀測前段時間內(nèi)橋梁下沉速率變化較快，而在中后期逐漸變慢。下沉的幅值前期較大，中間經(jīng)過一段時間段波動到后期逐漸趨向平穩(wěn)。可見從圖中無法識別橋梁下沉的外界影響因素與突變點，需從小波分解圖上發(fā)現(xiàn)其他頻率的干擾波。

圖1 下沉速度Fig.1 Deformation speed

4.1 小波閾值去噪

A1點監(jiān)測數(shù)據(jù)降噪后信噪比和均方根［3］誤差見表2，可見改進軟閾值消噪法去噪效果最好。由圖2可見，橋梁變形在去噪后，不同時期顯現(xiàn)了不同的變形規(guī)律。

表2 不同消噪方法比較Tab.2 Comparison of different de-noising effects

圖2 去噪后與原信號對比Fig.2 Comparison between signals before and after denoising

5 消噪后橋梁變形數(shù)據(jù)混沌時間序列預(yù)測

5.1 吸引子關(guān)聯(lián)維值計算

取延遲時間τ=5，作lnC(r)和lnr關(guān)系圖，由圖分析，當(dāng)維數(shù)增加到8時，lnC(r)和lnr比值不再隨m增加而改變，故取維數(shù)8為嵌入維數(shù)。經(jīng)計算吸引子關(guān)聯(lián)維值為2.8。

5.2 混沌特性判別

降噪后A1點監(jiān)測數(shù)據(jù)，采用改進小數(shù)據(jù)量法算法，計算最大Lyapunov指數(shù)值，計算結(jié)果為2.2×10－4，大于零。

將計算吸引子關(guān)聯(lián)維值和最大Lyapunov指數(shù)值，與表1中橋梁變形混沌特征判別標(biāo)準(zhǔn)對比，可見該時間序列具有混沌性質(zhì)。

5.3 混沌時間序列預(yù)測

取時間序列中前 16 個數(shù)據(jù){x1，x2，…，x16}，采用加權(quán)零階局域法和加權(quán)一階局域法分別對橋梁一次進行短期預(yù)測［4］，預(yù)測x17的值(表3)。

分別取指數(shù)平滑系數(shù)為0.9、0.5，采用指數(shù)平滑法［5］進行了預(yù)測，并與混沌預(yù)測結(jié)果進行對比，結(jié)果見表4。

表3 混沌時間序列法預(yù)測Tab.3 Forecast of chaotic time series method

表4 指數(shù)平滑法預(yù)測Tab.4 Forecast of exponential smoothing method

通過計算結(jié)果知，加權(quán)零階局域法和加權(quán)一階局域法的預(yù)測誤差明顯小于指數(shù)平滑法，加權(quán)一階局域法計算精度優(yōu)于加權(quán)零階局域法。預(yù)測結(jié)果表明混沌預(yù)測方法可靠、精度高，能真實地反映橋梁變形狀態(tài)。

6 結(jié)論

1)對跨海大橋變形監(jiān)測數(shù)據(jù)，選用db10小波函數(shù)進行橋梁預(yù)處理后，去噪效果較好，去噪后結(jié)果顯示了橋梁變化規(guī)律。

2)根據(jù)橋梁混沌變形狀態(tài)，提出混沌特征量計算方法，確定了橋梁變形混沌判別準(zhǔn)則。

3)對小波去噪后的數(shù)據(jù)建立了混沌預(yù)測模型，并將混沌極差預(yù)測值與指數(shù)平滑預(yù)測值進行比較，結(jié)果顯示了混沌預(yù)測方法精度較高。

1 劉云成，陳振學(xué)，馬于濤.自適應(yīng)閾值的小波圖像去噪［J］.電光工程，2007，34(6):77 －81.(Liu Yunxue，Chen Zhenxue and Ma Yutao.Adaptive wavelet thresholding method for image denoising［J］.Opto-Electronic Engineering，2007，34(6):77 －81)

2 韓敏.混沌時間序列預(yù)測理論與方法［M］.北京:中國水利水電出版社，1996.(Han Min.Chaotic time series prediction theory and method［M］.Beijing:China Water Power Press，1996)

3 Donoho D L.De-noising by soft-thresholding［J］.IEEE Transactions on Information Theory，1995，41:613 －627.

4 徐銘陶，秦樹人，陳志奎.實施一維Mallat算法的幾個問題［J］.重慶大學(xué)學(xué)報，1997，20(1):7 －11.(Xu Mingtao，Qin Shuren and Chen Zhikui.Some problems about the implementat ion of one-dimensional mallat algorithm［J］.Journal of Chongqing University，1997，20(1):7 －11)

5 張忠平.指數(shù)平滑法［M］.北京:中國統(tǒng)計出版，1996.(Zhang Zhongping.Exponential smoothing［M］.Beijing:China Statistics Press，1996)

WAVELET DENOISING AND CHAOS PREDICTION OF BRIDGE DEFORMATION DATA

Luan Yuanzhong1)，Luan Hengxuan2)，Ma Depeng3)，Liu Na4)and Zhou Yan1)
1)College of Geomatics of SUST，Qingdao266590
2)College of Mechanical and Electronic Engineering，Shandong University of SUST，Qingdao266590
3)College of Resources and Environmental Engineering of SUST，Qingdao266590
4)Surveying and Mapping Institute of Chongqing，Chongqing400042

Due to external factors such as waves，hurricanes and ship collision，cross-sea bridge deformation shows the dynamic nonlinear state，there is also a chaotic phenomenon.Therefore，db10 wavelet threshold de-noising and soft law are used to break down and eliminate the gross errors in the data distortion and noise.Phase space is reconstructed by the smooth deformation of the time-series after wavelet pretreatment，the bridge chaos state is determined by Lyapunov indices，and chaotic time series prediction and exponential smoothing forecast are analyzed and compared，it is proved that chaos prediction has reliable and high precision.

soft threshold de-noising;phase space reconstruction;Lyapunov index;chaotic time series;exponential smoothing method

P207

A

1671-5942(2013)05-0133-04

2013-03-12

欒元重，男，1963年5月生，博士，教授，博士生導(dǎo)師，主要研究方向為變形數(shù)據(jù)處理理論與方法.E－mail:mdp123@163.com