湯同旭，岳東杰

(河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 210098)

特長(zhǎng)鋼箱梁橋GPS監(jiān)測(cè)數(shù)據(jù)的奇異譜分析

湯同旭，岳東杰

(河海大學(xué) 地球科學(xué)與工程學(xué)院，江蘇 南京 210098)

實(shí)時(shí)準(zhǔn)確地監(jiān)測(cè)和分析大型建筑物的變形對(duì)于大型建筑物的施工與運(yùn)行是十分重要的。奇異譜分析(SSA)是一種與經(jīng)驗(yàn)正交函數(shù)相關(guān)聯(lián)的一種統(tǒng)計(jì)技術(shù)。文中采用奇異譜分析(SSA)的方法，結(jié)合功率譜峰值與原序列和RC1-2重建序列均方誤差最小的方法確定最優(yōu)潛入維數(shù)M，研究了特長(zhǎng)鋼箱梁橋索塔的變形趨勢(shì)和震蕩周期。小波分析技術(shù)也可以很好地描述時(shí)間序列的時(shí)頻分布，但SSA重構(gòu)后的結(jié)果比小波要平滑些，更能反映出時(shí)間序列的特征。結(jié)合蘇通大橋北索塔監(jiān)測(cè)序列表明，北索塔站的N，E,U方向上均存在明顯的趨勢(shì)和顯著的變動(dòng)周期，而且也含有較多的噪聲信息。

奇異譜；潛入維數(shù)；變形監(jiān)測(cè)；特征信息；小波分析

近年來，隨著橋梁監(jiān)測(cè)理論和信號(hào)處理的不斷發(fā)展和深入研究，應(yīng)用橋梁采集的數(shù)據(jù)進(jìn)行橋梁的損傷識(shí)別，荷載識(shí)別成為學(xué)術(shù)界的熱點(diǎn)。通過從這些原始數(shù)據(jù)中提取有用的信息，發(fā)揮變形監(jiān)測(cè)的作用尤為重要。由于受到日照、風(fēng)速、荷載、溫度和地脈動(dòng)的影響，鋼箱梁橋索塔的監(jiān)測(cè)數(shù)據(jù)一般具有非線性、隨機(jī)性的特征。這些影響可以分為周期項(xiàng)、趨勢(shì)向與噪聲3個(gè)部分，表現(xiàn)為多尺度信號(hào)的混合疊加。趨勢(shì)項(xiàng)可以反映出鋼箱梁橋面長(zhǎng)期變形的一種趨勢(shì)，是鋼箱梁橋面變形的主要部分。周期項(xiàng)反映出鋼箱梁橋面的變化規(guī)律，噪聲項(xiàng)一方面可以評(píng)估檢測(cè)技術(shù)的整體質(zhì)量，另一方面也可以反演分析檢測(cè)環(huán)境的變化。鋼箱梁索塔信號(hào)的處理就是去除噪聲項(xiàng)，提取周期項(xiàng)和趨勢(shì)向。

1 奇異譜分析的原理

Vautard，Yiou，Ghil 3人于1992年提出奇異譜分析(SSA)，在一個(gè)標(biāo)準(zhǔn)時(shí)間序列xi中，樣本指數(shù)i是從1到N變化的，嵌入維數(shù)(窗口長(zhǎng)度)為M，按照給定的潛入維數(shù)M建立時(shí)滯矩陣

(1)

該時(shí)滯矩陣(每個(gè)對(duì)角線都有統(tǒng)一固定值)共有N-M+1個(gè)狀態(tài)，對(duì)該時(shí)滯矩陣進(jìn)行經(jīng)驗(yàn)正交展開，就可以得到一個(gè)滯后協(xié)方差陣

Tx=

(2)

C(j)為時(shí)間序列時(shí)遲為j的自協(xié)方差，用公式表示為

(3)

這個(gè)時(shí)滯矩陣中的特征值是降序排列的，即

λ1≥λ2≥…≥λk≥0.

(4)

M.

(5)

它反映的是特征向量在原序列Xi+1，Xi+2，Xi+3，…，Xi+M時(shí)的權(quán)重。

每個(gè)由奇異譜分析識(shí)別的原始時(shí)間序列的主成分都可以被重建，奇異譜分析的重要功能是由成分重建RC(reconstructed component)實(shí)現(xiàn)的，第k個(gè)重建成分序列為

(6)

重建成分RC 還具有疊加性質(zhì)，所有重建成分疊加之和與原序列相同，即

(7)

包含在第k個(gè)重建成分的原始時(shí)間序列(相當(dāng)于標(biāo)準(zhǔn)化的時(shí)間序列)總體方差的部分就是特征值λk，因此，重建成分用過降低原始時(shí)間序列的信息來排序。大部分的方差都包含在前面的一些重建成分里面，大部分的或者說全部的剩余的重建成分包含著噪聲。SSA通常將一個(gè)周期小于M的周期性的時(shí)間序列分解為重建成分，在周期大于M的時(shí)間序列里面一或兩個(gè)重建成分含有方差。

2 蘇通大橋北索塔GPS變形監(jiān)測(cè)信號(hào)分析

2.1 蘇通大橋GPS變形監(jiān)測(cè)系統(tǒng)

為了便于蘇通大橋運(yùn)行期間安全監(jiān)測(cè)，建立一整套自動(dòng)化橋梁監(jiān)測(cè)系統(tǒng)，如圖1所示。該系統(tǒng)采用徠卡公司的GRX1200Pro系統(tǒng)，該系統(tǒng)包括1臺(tái)GPS基準(zhǔn)站和8臺(tái)GPS移動(dòng)站?；鶞?zhǔn)站架設(shè)在長(zhǎng)江岸邊一開闊地帶的帶有強(qiáng)制對(duì)中裝置的觀測(cè)墩上，2臺(tái)GPS移動(dòng)站架設(shè)在南北索塔頂上，其余6臺(tái)架設(shè)在主橋跨中以及跨中兩邊位置。該系統(tǒng)的坐標(biāo)系為橋軸線坐標(biāo)系，X方向平行于橋軸線方向，長(zhǎng)江北方向(N方向)為X軸的正方向，Y方向垂直于橋軸線的方向，往長(zhǎng)江下游去(E方向)為正方向。

圖1 蘇通大橋GPS變形監(jiān)測(cè)系統(tǒng)

大跨度橋梁的主頻一般都是0～0.5 Hz，根據(jù)奈奎斯斯(Nyqust)抽樣定理，接收機(jī)采樣頻率需要滿足

(8)

式中：fs為接收機(jī)的采樣頻率，Δt為接收機(jī)的采樣間隔，fc為信號(hào)能夠接收到的頻率。保證接收機(jī)的采樣頻率大于信號(hào)能接收頻率的2倍，因此,需要將接收機(jī)的采樣頻率設(shè)置為1 Hz，衛(wèi)星的截止高度角設(shè)限為15°。GPS雙頻接收機(jī)兩頻率對(duì)電離層的延遲有所不同，該系統(tǒng)利用雙頻接收機(jī)可以很大程度上消除電離層對(duì)電磁波延遲的影響。對(duì)于電離層折射對(duì)信號(hào)接收信號(hào)的影響，雙頻機(jī)能以L2觀測(cè)值加以改正。所以雙拼接收機(jī)適宜于一些中、長(zhǎng)基線(大于20 km)的測(cè)量，因此擁有快速靜態(tài)測(cè)量的功能，而且還可以升級(jí)為RTK測(cè)量的功能。利用RTK定位可以快速解算出載波相位的模糊度，在幾分鐘之內(nèi)就能完成初始化。

2.2 數(shù)據(jù)預(yù)處理

圖2 監(jiān)測(cè)點(diǎn)索塔三維E，N，U方向的時(shí)程曲線圖

2.3 潛入維數(shù)M的選取

對(duì)于SSA，如果嵌入維數(shù)M太大，譜分辨率越靈敏，但會(huì)導(dǎo)致奇異值的分解得到不同的部分產(chǎn)生混疊的情況。如果M太小，則無法使信號(hào)從弱到強(qiáng)進(jìn)行漸進(jìn)劃分，從而使得一些信號(hào)不能獲取。另外，外界環(huán)境條件對(duì)于嵌入維數(shù)M的選擇也是有影響的。對(duì)于給定的時(shí)間序列長(zhǎng)度N，為了使得最大滯后M-1的自協(xié)方差的統(tǒng)計(jì)誤差小于估計(jì)量本身，M不應(yīng)超過N/3。潛入維數(shù)M的選取也要考慮變形信號(hào)自身的周期性，最好的辦法就是選擇變形信號(hào)序列自身的周期作為潛入維數(shù)M。

為了得到功率譜估值，首先將信號(hào)進(jìn)行離散傅立葉變換，然后再用它的平方除以數(shù)據(jù)序列長(zhǎng)度N，即

(9)

(10)

由于數(shù)據(jù)序列x(n)的離散傅立葉變換是有周期的，因此這種功率譜也是具有周期的，通常稱作周期圖(PSD)，式中：f為時(shí)間序列的頻率，n為數(shù)據(jù)序列長(zhǎng)度。利用蘇通大橋2013-09-01至2012-09-27北索塔坐標(biāo)時(shí)間序列。對(duì)每天的時(shí)間序列取54個(gè)平均值，得到由這27每天總共1 458個(gè)時(shí)間序列組成的時(shí)間序列。則對(duì)此測(cè)點(diǎn)的E，N，U方向分別作PSD見圖3。

圖3 周期圖得到的功率譜曲線

由圖3可看出，E方向?qū)?yīng)的功率譜峰值為292，512，682；N方向?qū)?yīng)的功率譜峰值為227，292，512，682；U方向?qū)?yīng)的功率譜峰值為186，341，682。

理論上窗口長(zhǎng)度不應(yīng)大于N/3(N=1 461)，并且窗口長(zhǎng)度太小不能有效地反映頻分辨率，因此對(duì)于N，E方向取窗口長(zhǎng)度為2～450，U方向取窗口長(zhǎng)度為300～600，分別用第一第二時(shí)間主分量進(jìn)行重建，計(jì)算RC1-2重建序列跟實(shí)際序列的均方根誤差

(11)

式中:Xi為重建序列向量，Yi為實(shí)際序列向量。由此可以作出RC1-2重構(gòu)序列與實(shí)際序列的均方根誤差隨窗口的變化圖，如圖4所示。

圖4 RC1-2重建序列與實(shí)際序列均方根誤差隨窗口M的變化

由圖4可看出，N方向在M=292處存在拐點(diǎn)，E方向在M=292處存在拐點(diǎn)，U方向在M=341處存在拐點(diǎn)，結(jié)合圖3，所以可設(shè)N，E，U方向的窗口長(zhǎng)度為292，292，341，此時(shí)的均方根誤差RMSE分別為0.097 13 mm，0.116 8 mm，0.003 258 mm。

2.4 特征信息分析

蘇通大橋索塔坐標(biāo)時(shí)間序列是GNSS基準(zhǔn)站點(diǎn)長(zhǎng)時(shí)間觀測(cè)分析的結(jié)果，研究分析坐標(biāo)時(shí)間序列的特征，可以分析該索塔基準(zhǔn)點(diǎn)的位置變化規(guī)律，實(shí)現(xiàn)對(duì)基準(zhǔn)站點(diǎn)的穩(wěn)定性評(píng)價(jià)，對(duì)于GNSS觀測(cè)數(shù)據(jù)的應(yīng)用、索塔變形監(jiān)測(cè)以及如何更加合理地對(duì)索塔結(jié)構(gòu)健康狀況進(jìn)行安全評(píng)價(jià)都有重要意義。

國(guó)內(nèi)外的許多學(xué)者對(duì)GNSS坐標(biāo)時(shí)間序列進(jìn)行了分析研究，發(fā)現(xiàn)GNSS坐標(biāo)時(shí)間序列具有以下3個(gè)特征：①周期項(xiàng)成分；②趨勢(shì)項(xiàng)成分；③噪聲特性。本文主要對(duì)蘇通大橋北索塔坐標(biāo)時(shí)間序列的趨勢(shì)項(xiàng)成分和周期項(xiàng)成分進(jìn)行提取和分析。

2.4.1 周期項(xiàng)成分分析

周期項(xiàng)成分是指觀測(cè)序列中形態(tài)較為一致、時(shí)段較為同步、周期較為固定的一種變化。這種變化主要是由地球表面質(zhì)量負(fù)荷的季節(jié)性變化造成的，引發(fā)此類變化的原因有重力激發(fā)、熱效應(yīng)以及水文動(dòng)力學(xué)等。

分別對(duì)N，E，U3個(gè)方向進(jìn)行周期成分提取，將3個(gè)方向SSA重構(gòu)后的時(shí)間序列與原始時(shí)間序列進(jìn)行對(duì)比，如圖5所示，可以看出，經(jīng)過SSA重構(gòu)后的時(shí)間序列比原始時(shí)間序列更加光滑，因此SSA具有明顯的降噪效果。

2.4.2 趨勢(shì)項(xiàng)成分分析

趨勢(shì)項(xiàng)成分是指觀測(cè)序列以較為穩(wěn)定的速率保持平穩(wěn)、上升或者下降的變化。站點(diǎn)所在區(qū)域的地質(zhì)構(gòu)造運(yùn)動(dòng)引起了這種變化，該變化受此地區(qū)構(gòu)造應(yīng)力場(chǎng)的控制，且在一定的時(shí)間內(nèi)保持穩(wěn)定。

分別對(duì)N，E，U3個(gè)方向信號(hào)的有效主分量進(jìn)行重建序列，如圖6(a)、圖6(b)分別為N，E方向有效主分量信號(hào)重建序列圖，其中藍(lán)線代表處理后的原始時(shí)間序列，紅線代表變形信號(hào)趨勢(shì)成分。圖6(c)為U方向信號(hào)前6個(gè)有效主分量進(jìn)行重建的序列圖，橫坐標(biāo)為Time series/days，縱坐標(biāo)為dU/mm，藍(lán)線代表處理后的原始時(shí)間序列，紅線代表變形信號(hào)趨勢(shì)成分(RC1—RC6)。

圖6 有效主分量信號(hào)重建序列

2.5 與小波分析方法比較

小波分析也是分析時(shí)間序列的有效方法，它可以獲取時(shí)間序列的特征信息。該方法使用一種窗函數(shù)(小波函數(shù))，時(shí)頻窗面積不變，但是形狀可以改變。它可以根據(jù)需要調(diào)整時(shí)間和頻率分辨率，具有多分辨率分析的特點(diǎn)，在低頻部分具有較低的時(shí)間分辨率和較高的頻率分辨率，在高頻部分具有較低的頻率分辨率和較高的時(shí)間分辨率，克服了傅立葉變換分析非平穩(wěn)信號(hào)單一分辨率的困難，使得小波變換具有對(duì)信號(hào)的自適應(yīng)性。

將原始時(shí)間序列、SSA重構(gòu)序列、小波重構(gòu)序列進(jìn)行對(duì)比，如圖7所示。

圖7 原始時(shí)間序列、SSA重構(gòu)序列、小波重構(gòu)序列對(duì)比

由圖7可以看出，SSA和小波重構(gòu)都保留了時(shí)間序列的主要特征，且不同程度地剔除了部分噪聲，但是SSA重構(gòu)后的結(jié)果比小波要平滑些，更能反映時(shí)間序列的特征。

3 結(jié) 論

1)SSA可以有效地分離并且提取出時(shí)間序列中的趨勢(shì)項(xiàng)成分，進(jìn)而得到比較純粹的趨勢(shì)結(jié)果，減少了周期項(xiàng)對(duì)趨勢(shì)項(xiàng)的影響。

2)SSA可以較好地識(shí)別和提取時(shí)間序列中的周期項(xiàng)成分，北塔在N和E方向上存在著明顯的日周期，U方向上也具有日周期，變形大小不太具有規(guī)律性，可能與每日的溫度、日照不同有關(guān)。實(shí)驗(yàn)結(jié)果分析表明，SSA在周期選取方面具有較好的應(yīng)用價(jià)值。小波分析也能較好地識(shí)別和提取趨勢(shì)向。

3)SSA和小波重構(gòu)都保留了時(shí)間序列的主要特征，且都能剔除了部分噪聲，但是SSA重構(gòu)后的結(jié)果比小波要平滑些，更能反映出時(shí)間序列的特征。

4)通過對(duì)蘇通大橋北索塔原始序列進(jìn)行重構(gòu)，可以發(fā)現(xiàn)不論是E，N還是U方向上變形都具有可恢復(fù)性，不存在不和諧的趨勢(shì)性震蕩，表明索塔運(yùn)行安全可靠。

5)特長(zhǎng)鋼箱梁大橋的變形受內(nèi)部變形因子和外部環(huán)境因素的影響很多，本文是基于蘇通大橋外部GPS監(jiān)測(cè)數(shù)據(jù)進(jìn)行橋梁變形分析，內(nèi)部變形信號(hào)的提取與分析還需進(jìn)一步研究。

[1] 楊柏寧，徐愛林.GPS在蘇通大橋索塔動(dòng)態(tài)特性測(cè)試中的應(yīng)用[J].工程勘察，2009，37(3):56-59.

[2] 徐克飛，程鵬飛，文漢江.太陽黑子數(shù)時(shí)間序列的奇異譜分析和小波分析[J].測(cè)繪科學(xué)，2007，32(6):35.

[3] 王解先，連麗珍，沈云中.奇異譜分析在GPS站坐標(biāo)監(jiān)測(cè)序列分析中的應(yīng)用[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2013,41 (2):282-288.

[4] 許昌，岳東杰.基于RTK-GPS技術(shù)的高索塔振動(dòng)試驗(yàn)與分析[J].振動(dòng)與沖擊，2010，29(3):134-136.

[5] 劉成龍，張德強(qiáng)，黃澤純.大跨懸索橋施工索塔變形成因分析與監(jiān)測(cè)[J].西南交通大學(xué)學(xué)報(bào)，2000，35(5)：501-504.

[6] 潘國(guó)榮，谷川.GPS監(jiān)測(cè)橋梁擺動(dòng)及其數(shù)據(jù)的小波分析[J].工程勘察，2008(10):55-58.

[7] 姚連璧，姚平，王人鵬，等.南浦大橋形變GPS動(dòng)態(tài)監(jiān)測(cè)試驗(yàn)及結(jié)果分析[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2008,36(12):1633-1636.

[8] 許福友，陳艾榮.蘇通大橋三維顫振分析[J].工程力學(xué)，2008，25(8):139-144.

[9] 許映梅，岳東杰，袁豹.基于小波去噪的蘇通大橋索塔GPS監(jiān)測(cè)數(shù)據(jù)分析[J].水利與建筑工程學(xué)報(bào)，2013，11(2):1-6，22.

[責(zé)任編輯：劉文霞]

Singular spectrum analysis of GPS monitoring data on specialty steel box girder bridge

TANG Tongxu,YUE Dongjie

(School of Earth Science and Engineering，Hohai University，Nanjing 210098，China)

Real-time accurate monitoring and analysis on deformation of large buildings for construction and operation are very important.Singular spectrum analysis (SSA) is a kind of statistical technique associated with the empirical orthogonal function.This paper adopts the method of singular spectrum analysis (SSA),combined with a method of power spectrum peak rate with the original sequence and RC1-2 reconstruction sequences of minimum mean square error to determine the optimal dimension M. The deformation trend and the oscillation cycle of cable tower are studied for specialty steel box girder bridge.Wavelet analysis technology can well describe the time sequence of the time-frequency distribution.But SSA reconstruction result is more smoothing than wavelet analysis,which can reflect the characteristics of time series better.North cable tower of the Sutong bridge monitoring sequence shows that N,E and U direction of north cable tower stay with obvious trends and significant changes in the cycle,and also contain more information noise.

singular spectrum;embedding dimensions;deformation monitoring;feature information;wavelet analysis

2016-07-01

湯同旭(1992-)，男，碩士研究生.

著錄：湯同旭，岳東杰.特長(zhǎng)鋼箱梁橋GPS監(jiān)測(cè)數(shù)據(jù)的奇異譜分析[J].測(cè)繪工程，2017，26(8)：54-61.

10.19349/j.cnki.issn1006-7949.2017.08.012

P228

A

1006-7949(2017)08-0054-08