王 濤，田林亞，侯建梅

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

小波與奇異譜分析在地鐵保護(hù)區(qū)監(jiān)測(cè)數(shù)據(jù)處理中的應(yīng)用

王 濤，田林亞，侯建梅

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

以南京市地鐵沉降監(jiān)測(cè)數(shù)據(jù)為例，將小波分析和奇異譜分析方法應(yīng)用在地鐵保護(hù)區(qū)監(jiān)測(cè)數(shù)據(jù)的處理中。實(shí)驗(yàn)證明本方法在探測(cè)異常值、獲取沉降變形趨勢(shì)、變形周期等方面具有可行性。

地鐵；沉降監(jiān)測(cè)；變形趨勢(shì)；異常值探測(cè)；小波分析；奇異譜分析

全國(guó)范圍內(nèi)興建了大量的地鐵，為保證地鐵結(jié)構(gòu)的安全，應(yīng)該對(duì)其進(jìn)行全方位監(jiān)測(cè)[1-3]。掌握既有地鐵工程結(jié)構(gòu)在新建項(xiàng)目施工過(guò)程中發(fā)生的變化，為建設(shè)單位及地鐵單位提供及時(shí)、可靠的信息，準(zhǔn)確判斷工程施工對(duì)既有地鐵工程結(jié)構(gòu)的影響，避免惡性事故的發(fā)生。其中，自動(dòng)化、智能化監(jiān)測(cè)是今后地鐵安全監(jiān)測(cè)的主要發(fā)展方向[4]，如何實(shí)時(shí)準(zhǔn)確處理大量監(jiān)測(cè)數(shù)據(jù)，剔除監(jiān)測(cè)異常值、獲取地鐵變形趨勢(shì)和變形周期是一個(gè)值得研究的問(wèn)題。相關(guān)學(xué)者進(jìn)行大量的研究，并采用自適應(yīng)卡爾曼濾波、時(shí)間序列、小波時(shí)間序列等方法進(jìn)行地鐵數(shù)據(jù)處理[5-7]，但是沒(méi)有將奇異譜和小波分析應(yīng)用于地鐵沉降監(jiān)測(cè)的數(shù)據(jù)處理中。

小波分析和奇異譜分析是當(dāng)前處理時(shí)間序列數(shù)據(jù)十分有效的方法[8-10]，通過(guò)小波分析和奇異譜分析可以有效去除觀測(cè)數(shù)據(jù)中的噪聲和異常值，獲取變形趨勢(shì)和周期。本文將小波分析和奇異譜分析方法應(yīng)用于地鐵監(jiān)測(cè)數(shù)據(jù)處理，并以正在實(shí)施監(jiān)測(cè)的南京地鐵2號(hào)線某地鐵區(qū)間的監(jiān)測(cè)數(shù)據(jù)進(jìn)行實(shí)驗(yàn)分析，獲得較好的效果，本文的研究為解讀和處理地鐵監(jiān)測(cè)數(shù)據(jù)提供參考。

1 小波與奇異譜分析理論

1.1 小波分析

小波分析(Wavelet Analysis)是調(diào)和分析領(lǐng)域半個(gè)多世紀(jì)以來(lái)工作的結(jié)晶，近年來(lái)成為國(guó)際上眾多學(xué)科共同關(guān)注的熱點(diǎn)[11-12]。小波分析的核心是小波變換，小波變換可以視為被變換函數(shù)與小波基函數(shù)作卷積的結(jié)果。設(shè)滿足一定條件的小波基函數(shù)為φ(t)，經(jīng)過(guò)平移τ和伸縮a后得到小波基函數(shù)系：

(1)

然后將φa,τ(t)與被變換函數(shù)f(t)作卷積，得到連續(xù)小波變換公式：

WTf(a,τ)=[f(t),φa,τ(t)]=

(2)

Z.

(3)

相應(yīng)的離散小波變換公式為

WTf(a,τ)=〈f(t),φa,τ(t)〉=

(4)

小波分析具有優(yōu)異的局部時(shí)-頻特性，通過(guò)伸縮和平移等運(yùn)算功能對(duì)信號(hào)進(jìn)行多尺度細(xì)化分析，可以聚焦信號(hào)的任意細(xì)節(jié)，有效地從信號(hào)中提取信息。將監(jiān)測(cè)數(shù)據(jù)序列視為不同頻率成分組成的數(shù)字信號(hào)，采用小波分析方法進(jìn)行去噪和變形趨勢(shì)提取，被認(rèn)為是有效的方法[13]。

1.2 奇異譜分析

奇異譜分析(Singular Spectrum Analysis,SSA)是近年來(lái)興起的用于研究非線性時(shí)間序列十分有效的數(shù)據(jù)分析方法[14]。核心思想：將一維時(shí)間序列轉(zhuǎn)換為多維序列進(jìn)行主成分分析，分解為多個(gè)有效成分，并根據(jù)需要選擇若干分量進(jìn)行重構(gòu)，提取出信號(hào)的變化趨勢(shì)、周期和噪聲等不同成分的信號(hào)。

奇異譜分析時(shí)，首先將一維時(shí)間序列{xi}=x1,x2,x3，…，xN.轉(zhuǎn)化為時(shí)滯矩陣X:

(5)

其中，M為嵌入維數(shù)。

然后，對(duì)時(shí)滯矩陣X進(jìn)行時(shí)間經(jīng)驗(yàn)正交函數(shù)展開(kāi)(TEOFs)，得到滯后協(xié)方差陣D

(6)

求解D的特征向量Ek和特征值：λ1≥λ2≥…≥λM，這些特征值被稱(chēng)為時(shí)間序列{xi}的奇異譜，對(duì)奇異譜進(jìn)行重建的運(yùn)算稱(chēng)為奇異譜分析。λk對(duì)應(yīng)的特征向量Ek為時(shí)間經(jīng)驗(yàn)正交函數(shù)(TEOF)，定義第k個(gè)時(shí)間主分量(TPC)為序列{xi}在Ek上的正交投影系數(shù)，aik代表權(quán)重，

M.

(7)

SSA中最重要的步驟是重建(Reconstruction component；RC)[15]，RC的過(guò)程就是由TEOF和TPC重建一個(gè)與原始序列等長(zhǎng)的序列，其式為

(8)

所有重建成分之和等于原始序列。

2 工程背景

2.1 工程概況

本文以正在實(shí)施監(jiān)測(cè)的南京市地鐵2號(hào)線保護(hù)區(qū)監(jiān)測(cè)項(xiàng)目為例進(jìn)行分析，由于2號(hào)線地鐵附近新建項(xiàng)目的實(shí)施，必須要實(shí)時(shí)監(jiān)測(cè)地鐵受新建項(xiàng)目開(kāi)挖及建設(shè)過(guò)程中的影響，以保障地鐵的安全運(yùn)營(yíng)。擬建酒店及酒店式公寓，設(shè)置一層聯(lián)體地下室。其中，基坑面積約18 540 m2，周長(zhǎng)約560 m，開(kāi)挖深度約為9.2 m。地形較平坦，根據(jù)野外勘探，原位測(cè)試和室內(nèi)巖土試驗(yàn)成果綜合分析，場(chǎng)地巖土層自上而下分別是雜填土、粉質(zhì)粘土混粗砂、卵礫石、殘積土、淤泥質(zhì)粉質(zhì)粘土、中風(fēng)化粉砂質(zhì)泥巖、淤泥質(zhì)粉質(zhì)粘土。建設(shè)場(chǎng)地東部以及西側(cè)坡角處各有一條泄洪溝，走向南北，地表水主要接受大氣降水以及上游生活污水的補(bǔ)給，以流向下游、蒸發(fā)和逐漸下滲的方式排泄，和地下水之間存有互補(bǔ)關(guān)系。

基坑與地鐵的位置關(guān)系：①平面位置關(guān)系?；游挥诘罔F二號(hào)線孝陵衛(wèi)站至鐘靈街站礦山法隧道北側(cè)，其邊線對(duì)應(yīng)地鐵里程大概為K21+017～K21+145，長(zhǎng)度約為128 m。基坑支護(hù)結(jié)構(gòu)邊線距地鐵隧道結(jié)構(gòu)邊線最近約為25.4 m。②豎向位置關(guān)系?；舆吘€對(duì)應(yīng)地鐵礦山法隧道頂標(biāo)高為2.96～5.92 m(相對(duì)±0.0為-10.88～-13.84 m)，基坑墊層底標(biāo)高為7.8 m，比隧道頂部高1.88～4.84 m。其位置示意圖如圖1所示。

圖1 基坑基監(jiān)測(cè)點(diǎn)位置示意圖

2.2 監(jiān)測(cè)方案

如圖1所示，沉降監(jiān)測(cè)基準(zhǔn)點(diǎn)布設(shè)在遠(yuǎn)離地鐵變形區(qū)120 m的軌道底板上，在左、右線各布設(shè)4個(gè)基準(zhǔn)點(diǎn)(JZ1,JZ2,JZ3,JZ4；JY1,JY2,JY3,JY4)，在基坑邊線對(duì)應(yīng)區(qū)間隧道范圍內(nèi)，左線約每10 m布設(shè)一個(gè)監(jiān)測(cè)點(diǎn)(Z3～Z16)，右線約每20 m布設(shè)一個(gè)監(jiān)測(cè)點(diǎn)(Y3～Y16)。監(jiān)測(cè)采用Trimble DINI03電子水準(zhǔn)儀，配合銦瓦條碼尺和夜間照明燈具，根據(jù)相關(guān)規(guī)范采用科傻平差軟件以測(cè)站數(shù)為權(quán)進(jìn)行嚴(yán)密平差?；鶞?zhǔn)網(wǎng)初期進(jìn)行兩次獨(dú)立觀測(cè)，誤差范圍內(nèi)取均值確定各基準(zhǔn)點(diǎn)的初始高程。每月定期對(duì)基準(zhǔn)網(wǎng)進(jìn)行復(fù)測(cè)，檢測(cè)實(shí)測(cè)高程值與原采用值進(jìn)行比較，若檢測(cè)的高程值與原值的較差大于2倍高程中誤差時(shí)采用新值，否則采用原值。平時(shí)采用附和水準(zhǔn)路線，從基準(zhǔn)點(diǎn)JZ2(JY2)測(cè)至JZ3(JY3)。

3 實(shí)驗(yàn)分析

從2015年11月15開(kāi)始，到目前為止監(jiān)測(cè)了57期。實(shí)驗(yàn)選擇左線的Z4、Z10、Z15點(diǎn)的觀測(cè)值進(jìn)行分析，Z4和Z15均位于基坑邊線，Z10位于基坑中間位置，原始觀測(cè)序列如圖2所示。從圖2中無(wú)法直接獲取觀測(cè)點(diǎn)變化趨勢(shì)，無(wú)法得到觀測(cè)中存在的異常值，這使得事后分析和解讀地鐵變形情況變得困難，實(shí)際上，隨著自動(dòng)化監(jiān)測(cè)工作的開(kāi)展，必須采用新的方法處理大量的監(jiān)測(cè)數(shù)據(jù)。本文以沉降監(jiān)測(cè)為例，研究小波分析和奇異譜分析在地鐵保護(hù)區(qū)監(jiān)測(cè)數(shù)據(jù)處理中的應(yīng)用。

圖2 原始觀測(cè)序列

圖3 小波分析結(jié)果

首先對(duì)觀測(cè)序列采用小波分析方法處理，實(shí)驗(yàn)時(shí)選擇去噪效果好的Sym3小波基進(jìn)行3層小波分解，并采用全局閾值的方法進(jìn)行小波閾值降噪[16]，然后進(jìn)行小波重構(gòu)，得到去噪后的觀測(cè)序列如圖3所示。觀測(cè)序列變得光滑，并且得到了趨勢(shì)明顯的變化，說(shuō)明采用小波分析取得了很好的效果。整體上看，Z4,Z10,Z15三個(gè)監(jiān)測(cè)點(diǎn)的變化趨勢(shì)基本相同，認(rèn)為基坑開(kāi)挖對(duì)地鐵線路中心區(qū)段Z3—Z16的沉降影響基本一致。容易看出在Z15觀測(cè)序列的第55期觀測(cè)值為異常值，經(jīng)對(duì)比原始數(shù)據(jù)檢驗(yàn)該點(diǎn)在觀測(cè)時(shí)前后兩次觀測(cè)值較差大于10 mm，可以判斷其為粗差，應(yīng)給予剔除。

進(jìn)一步觀察觀測(cè)序列變化趨勢(shì)為在前20期表現(xiàn)為在0值附近浮動(dòng)，可以認(rèn)為在觀測(cè)初期新建建筑對(duì)地鐵的影響不明顯；在20～50期表現(xiàn)為上升的趨勢(shì)，監(jiān)測(cè)到該時(shí)期基坑進(jìn)行大量的開(kāi)挖，并且地鐵內(nèi)部出現(xiàn)滲漏的情況，圖4為其中幾處明顯滲漏處，說(shuō)明這個(gè)時(shí)期基坑開(kāi)挖對(duì)地鐵的影響明顯，需要特別注意此階段的沉降量，并加強(qiáng)監(jiān)測(cè)的頻率；51期以后表現(xiàn)為逐漸下降到0值以下的趨勢(shì)，說(shuō)明地鐵線路在受到基坑影響上升后又隨著建筑物基礎(chǔ)的建設(shè)又逐步下沉。以上表現(xiàn)均符合相關(guān)工程的一般變化規(guī)律。

現(xiàn)在對(duì)原始觀測(cè)序列進(jìn)行奇異譜分析，根據(jù)選擇窗口大小不大于觀測(cè)序列1/3的原則，選擇窗口M=18，進(jìn)行奇異譜分析獲取觀測(cè)序列的趨勢(shì)項(xiàng)和周期項(xiàng)，結(jié)果如圖5和圖6所示。從圖5中可以看出，奇異譜分析較好的擬合了觀測(cè)序列的變化趨勢(shì)，呈現(xiàn)出了地鐵沉降變化特性，相對(duì)于小波分析，奇異譜分析結(jié)果更光滑，擬合效果更好，但是無(wú)法探測(cè)出可能出現(xiàn)的異常值。在變化趨勢(shì)上，兩者的表現(xiàn)基本一致，但是仔細(xì)觀察不難發(fā)現(xiàn)奇異譜分析對(duì)變形趨勢(shì)的具體特征的表達(dá)不如小波分析。圖6為提取出了周期項(xiàng)，在地鐵監(jiān)測(cè)時(shí)可以據(jù)此了解監(jiān)測(cè)周期和頻率，制定合理的監(jiān)測(cè)計(jì)劃。在圖6中，以Z15觀測(cè)點(diǎn)為例進(jìn)行分析，根據(jù)變形大小不同將其變形周期分為三段，第一段為前20期，第二段為21～52期，第三段為53～57期。在第一段中，由于沉降量值較小且周期穩(wěn)定，可以據(jù)此減小監(jiān)測(cè)頻率，以減少工程監(jiān)測(cè)成本；而從第三段明顯可以看出沉降值變大且周期也變大，說(shuō)明當(dāng)下的監(jiān)測(cè)頻率已經(jīng)不能完全表達(dá)變形趨勢(shì)，需要增加監(jiān)測(cè)頻率。

圖4 地鐵隧道內(nèi)部滲漏情況

圖5 奇異譜分析趨勢(shì)項(xiàng)

圖6 奇異譜分析周期項(xiàng)

面對(duì)海量的地鐵觀測(cè)數(shù)據(jù)，傳統(tǒng)的數(shù)據(jù)處理方法已經(jīng)無(wú)法滿足要求，借助小波分析和奇異譜分析及時(shí)發(fā)現(xiàn)地鐵軌道的變形現(xiàn)狀及發(fā)展趨勢(shì)，采取有效地處理措施和預(yù)案來(lái)保障地鐵的運(yùn)營(yíng)安全是行之有效的方法。小波分析和奇異譜分析都能夠去除原始觀測(cè)序列噪聲，獲取觀測(cè)體的變化趨勢(shì)，其中奇異譜分析擬合出的趨勢(shì)更光滑，并且可以獲得變形周期，小波分析表達(dá)了更多的變化細(xì)節(jié)，而且具有探測(cè)異常值的能力。兩者對(duì)地鐵監(jiān)測(cè)數(shù)據(jù)的處理都是有效的，在工程應(yīng)用時(shí)要根據(jù)不同目的選擇使用哪一種方法。

4 結(jié)束語(yǔ)

地鐵監(jiān)測(cè)數(shù)據(jù)是掌握地鐵結(jié)構(gòu)工作狀態(tài)的重要指標(biāo)，在對(duì)地鐵的沉降、水平位移、收斂等方面進(jìn)行監(jiān)測(cè)時(shí)，準(zhǔn)確實(shí)時(shí)的處理大量監(jiān)測(cè)數(shù)據(jù)，獲取變形趨勢(shì)和異常值，進(jìn)一步預(yù)測(cè)變化趨勢(shì)和周期對(duì)地鐵監(jiān)測(cè)工作的實(shí)施和安排是十分有益的，本文將當(dāng)前處理時(shí)間序列十分熱門(mén)的小波分析和奇異譜分析理論應(yīng)用于地鐵安全監(jiān)測(cè)中，并以南京地鐵的沉降監(jiān)測(cè)數(shù)據(jù)為例進(jìn)行實(shí)驗(yàn)分析，獲得較好的處理結(jié)果，本文的研究為地鐵監(jiān)測(cè)中的數(shù)據(jù)處理提供了有益的參考。

[1] 王超,夏穎. 軌道交通安全保護(hù)區(qū)信息化監(jiān)測(cè)技術(shù)研究[J]. 湖南文理學(xué)院學(xué)報(bào)(自然科學(xué)版),2013(2):57-60，87.

[2] 鄒積亭,江恒彪. 北京地鐵沉降監(jiān)測(cè)方法及數(shù)據(jù)處理[J]. 工程勘察,2006(1):43-45.

[3] 賀躍光,吳盛才,徐鵬. 城市地鐵安全監(jiān)測(cè)數(shù)據(jù)分析與處理[J]. 工程勘察,2011(12):46-49，54.

[4] 王超,夏穎. 軌道交通安全保護(hù)區(qū)信息化監(jiān)測(cè)技術(shù)研究[J]. 湖南文理學(xué)院學(xué)報(bào)(自然科學(xué)版),2013(2):57-60，87.

[5] 文小勇. 基于小波分析的自適應(yīng)卡爾曼濾波在地鐵變形監(jiān)測(cè)中的應(yīng)用[D].西安：長(zhǎng)安大學(xué),2015.

[6] 郭超. 地鐵地表沉降監(jiān)測(cè)數(shù)據(jù)的預(yù)處理及預(yù)測(cè)方法的探討[D].南昌：東華理工大學(xué),2015.

[7] 周紅斌. 基于時(shí)間序列分析的地鐵基坑變形監(jiān)測(cè)數(shù)據(jù)處理方法研究[D].沈陽(yáng)：東北大學(xué),2008.

[8] 章浙濤. 小波分析理論及其在變形監(jiān)測(cè)中的應(yīng)用研究[D].長(zhǎng)沙：中南大學(xué),2014：12-14.

[9] 曹奇,岳東杰,高永攀,等. 基于非平穩(wěn)時(shí)間序列的不同趨勢(shì)項(xiàng)提取方法對(duì)比研究[J]. 大地測(cè)量與地球動(dòng)力學(xué),2013(6):150-154.

[10] 周天清. 基于奇異譜分析的金融時(shí)間序列自適應(yīng)分解預(yù)測(cè)研究[D].南昌：華東交通大學(xué),2012：9-13.

[11] 黨星海,趙麗潔,孔令杰,等. 小波分析在GPS振動(dòng)監(jiān)測(cè)數(shù)據(jù)中的應(yīng)用[J]. 大地測(cè)量與地球動(dòng)力學(xué),2013(2):147-150.

[12] VEER K, AGARWAL R. Wavelet and short-time Fourier transform comparison-based analysis of myoelectric signals[J].Journal of Applied Statistics, 2015,42(7)：1591-1601.

[13] 程龍，張曉梅.小波神經(jīng)網(wǎng)絡(luò)模型在高鐵路基沉降預(yù)測(cè)中的應(yīng)用研究[J].測(cè)繪與空間地理信息，2016,39(4)：218-221.

[14] 張一,惠曉峰. 基于奇異譜分析的匯率預(yù)測(cè)研究[J]. 統(tǒng)計(jì)與決策,2012(6):29-31.

[15] 翟長(zhǎng)治,岳東杰,岳順. 大型橋梁GPS監(jiān)測(cè)數(shù)據(jù)的奇異譜分析[J]. 勘察科學(xué)技術(shù),2014(4):20-23，61.

[16] 吳繼忠. GPS觀測(cè)數(shù)據(jù)的小波閾值法消噪[J]. 大地測(cè)量與地球動(dòng)力學(xué),2009(4):79-82，87.

[責(zé)任編輯：李銘娜]

Application of wavelet analysis and singular spectrum analysis to subway reserve monitoring data processing

WANG Tao，TIAN Linya，HOU Jianmei

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

Monitoring the urban subway safety is very important. This paper, taking Nanjing subway subsidence monitoring data as an example, applies the wavelet analysis and singular spectrum analysis method to the subway reserves monitoring data processing and analysis. The experiment proves the effectiveness of this method on detecting the deformation trend and deformation cycle.

subway; subsidence monitoring; deformation trend; outlier detection; wavelet analysis; singular spectrum analysis

2016-08-28

江蘇省研究生科研與實(shí)踐創(chuàng)新計(jì)劃項(xiàng)目(KYCX17-0506)

王 濤(1992-)，男，碩士研究生.

著錄：王濤，田林亞，侯建梅，等.小波與奇異譜分析在地鐵保護(hù)區(qū)監(jiān)測(cè)數(shù)據(jù)處理中的應(yīng)用[J].測(cè)繪工程，2017，26(9)：60-64.

10.19349/j.cnki.issn1006-7949.2017.09.013

P208

A

1006-7949(2017)09-0060-05