梁 策,黃元庫,劉俊飛,王萬齊,朱 軍

(1. 中國鐵道科學(xué)研究院集團(tuán)有限公司電子計算技術(shù)研究所,北京 100081; 2. 西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院,四川 成都 611756; 3. 陜西交控工程技術(shù)有限公司,陜西 西安 710117;4. 中國國家鐵路集團(tuán)有限公司科技和信息化部,北京 100844)

對在建高速鐵路工程主體進(jìn)行持續(xù)沉降觀測,用沉降變形數(shù)據(jù)分析和評估工程本體是否趨于穩(wěn)定,可輔助決策鋪軌時機(jī),以此為前提基礎(chǔ)保障高速鐵路線路平順性[1-5]。因此,沉降變形數(shù)據(jù)的客觀性和準(zhǔn)確性是沉降觀測工作的根本。然而,實際的沉降觀測過程中,由于測量錯誤、測點(diǎn)或基準(zhǔn)點(diǎn)被碰撞變形等原因,會出現(xiàn)“失真”的數(shù)據(jù),這些數(shù)據(jù)通常表現(xiàn)為沉降-時間曲線上的突變(跳動)等狀態(tài)。當(dāng)觀測數(shù)據(jù)所反映的沉降變形與工程本體正常的沉降變形規(guī)律不符,或者與設(shè)計階段對沉降變形的預(yù)期不一致時,稱這些觀測數(shù)據(jù)為異常數(shù)據(jù)[6]。如何利用已經(jīng)采集并積累的大量沉降觀測數(shù)據(jù),及時、快速判識出因測量錯誤引起、在沉降-時間曲線上呈突變特征的異常數(shù)據(jù),是在建高速鐵路沉降觀測信息化管理中亟待解決的關(guān)鍵技術(shù)問題。

在建筑物的異常沉降數(shù)據(jù)辨識過程中,對于沉降觀測頻次密集的數(shù)據(jù),小波分析法能較好地適配數(shù)據(jù)波動變化,具有應(yīng)用優(yōu)勢[7-8],但是面對高速鐵路沉降觀測時間間隔較長且隨工況多變的情況,則存在不適用的情況。設(shè)置誤差限閾值,并用最小二乘支持向量機(jī)的函數(shù)擬合預(yù)處理沉降觀測數(shù)據(jù)的方法[9-11],利用邊界值作為辨識的標(biāo)準(zhǔn),值得參考,但當(dāng)觀測數(shù)據(jù)量龐大時,擬合計算量逐級增加,且靈敏性不足。因此,有必要從適配觀測數(shù)據(jù)的波動、限值設(shè)置、計算便捷程度及靈敏性等方面進(jìn)行研究。

本文基于沉降觀測歷史大數(shù)據(jù)與數(shù)理統(tǒng)計方法,擬制定路基、橋涵、隧道各主要工況對應(yīng)的判識異常沉降數(shù)據(jù)的參考閾值,提出利用平均變形速率作為代表計算公式,在顧及計算結(jié)果靈敏性的優(yōu)先條件下,辨識高速鐵路沉降異常數(shù)據(jù)。

1 變形速率判識沉降異常數(shù)據(jù)方法

1.1 主要步驟

首先,確定鐵路路基、橋涵、隧道各主要工況判識沉降異常數(shù)據(jù)的參考閾值;其次,識別并剔除無效觀測數(shù)據(jù),僅選取最近1～5期的有效觀測數(shù)據(jù)參與計算;然后,應(yīng)用沉降變形速率代表計算公式進(jìn)行計算;最后,對照參考閾值進(jìn)行判識,對應(yīng)觸發(fā)異常數(shù)據(jù)處置選項,形成管理閉環(huán)。技術(shù)路線如圖1所示。

1.2 確定參考閾值

基于沉降觀測歷史大數(shù)據(jù),采用數(shù)理統(tǒng)計方法,分別計算路橋隧各工況對應(yīng)的沉降變形速率較大值,然后疊加測量誤差的影響,形成參考閾值[12-13]。主要公式為

(1)

以無砟軌道鐵路路基填筑完成后第1～3個月的工況為例(如圖2所示),經(jīng)過統(tǒng)計,正常情況下路基沉降變形速率的較大值處于曲線1所在位置,這一時段對應(yīng)的正常沉降變形速率較大值約為0.1～0.3 mm/d。此時觀測頻次為7 d,經(jīng)過統(tǒng)計計算,測量誤差約為1 mm,3期測量的時間跨度內(nèi)測量誤差帶來的沉降變形速率偏差約為0.1 mm/d。正常變形速率較大值與誤差疊加后為0.4 mm/d,可以作為判識異常數(shù)據(jù)的參考閾值,如圖2中的折線2。

1.3 剔除無效觀測數(shù)據(jù)

測點(diǎn)“重新布設(shè)”“接管”或“測點(diǎn)轉(zhuǎn)移”時,數(shù)據(jù)修正會消除本次測量與上次測量間實際發(fā)生的沉降,本測點(diǎn)的此次測量視為無效觀測數(shù)據(jù);無效觀測數(shù)據(jù)被識別出來后將被剔除,不參與判識計算,剔除公式為

[Si]=R(S′1,S′2,…,S′j)

(2)

式中,從S′1到S′j為由近到遠(yuǎn)各期次觀測值;[Si]為有效觀測數(shù)據(jù)序列;i最大取到5即止。

1.4 平均變形速率代表計算

相對而言,計算沉降變形速率時選擇的起止時間跨度越大,參與速率計算的有效數(shù)據(jù)越多,越能更好地平均掉隨機(jī)誤差,速率數(shù)據(jù)結(jié)果穩(wěn)定性越好。然而,時間跨度過大,反映當(dāng)前速率的及時性會變差,將會損害沉降變形速率對異常數(shù)據(jù)的靈敏性[14-15]。

為了平衡沉降變形速率計算的穩(wěn)定性和靈敏性,本文提出靈敏性優(yōu)先與計算結(jié)果輔助校驗穩(wěn)定性的策略,即選用最近1～5期有效數(shù)據(jù),優(yōu)先保障靈敏性,以各期數(shù)據(jù)的計算結(jié)果輔助校驗計算值的波動量,并應(yīng)用平均變形速率計算公式作為代表。

平均變形速率計算公式為

(3)

單次變形速率靈敏性提取公式為

Vs1=Max(V2-1,V3-2,V4-3,V5-4)

(4)

式中,V2-1、V3-2、V4-3、V5-4為相鄰兩期觀測值之間的變形速率。

平均變形速率靈敏性提取公式為

Vs2=Max(V3-1,V4-1,V5-1)

(5)

式中,V3-1、V4-1、V5-1分別為第3、4、5期觀測值的平均變形速率。

用最大平均變形速率做穩(wěn)定性輔助校驗,得到變形速率代表計算公式為

(6)

式中,k=1.5為靈敏性調(diào)整系數(shù)。

同一里程中各測點(diǎn)的計算結(jié)果不同時,取各測點(diǎn)的最大計算值,作為該里程斷面的沉降變形速率,即

Vs=Max(Vsk)

(7)

1.5 異常數(shù)據(jù)判識過程及規(guī)則

(1)用式(6)計算測點(diǎn)Pi的變形速率Vsk。

(2)用式(7)計算同一里程斷面的沉降變形速率Vs。

(3)沉降異常數(shù)據(jù)的判斷規(guī)則為Vs≥[Vs]。其中,[Vs]為判識沉降異常數(shù)據(jù)的參考閾值,按表1—表4取值。

表1 異常數(shù)據(jù)的處置選項

2 工程應(yīng)用

2.1 鐵路沉降變形觀測信息系統(tǒng)實現(xiàn)

結(jié)合實際工程應(yīng)用需求,集成本文方法,設(shè)計并研發(fā)了鐵路沉降變形觀測信息系統(tǒng)(簡稱信息系統(tǒng)),具有沉降觀測數(shù)據(jù)采集、水準(zhǔn)線路管理、測點(diǎn)沉降-時間曲線數(shù)據(jù)分析、觀測頻次提醒、數(shù)據(jù)異常識別與處置等功能。

2.2 異常沉降數(shù)據(jù)提醒

信息系統(tǒng)于每天晚上集中進(jìn)行異常數(shù)據(jù)判識,如果存在異常數(shù)據(jù),將在沉降觀測數(shù)據(jù)上傳到信息系統(tǒng)后的第2天開始提醒。以標(biāo)段為單位,按照關(guān)注、加強(qiáng)關(guān)注、重點(diǎn)關(guān)注,共3個級別進(jìn)行分類,分別統(tǒng)計橋涵、路基、隧道各專業(yè)每天出現(xiàn)的沉降異常斷面數(shù)量,界面如圖3所示。

圖3 近7日需關(guān)注的沉降異常斷面數(shù)量統(tǒng)計界面

2.3 異常數(shù)據(jù)處置選項設(shè)置

當(dāng)發(fā)現(xiàn)沉降異常數(shù)據(jù)時,在測量數(shù)據(jù)核對、測點(diǎn)或基準(zhǔn)點(diǎn)檢查、施工過程分析、周圍環(huán)境調(diào)查等方面進(jìn)行操作,并結(jié)合異常數(shù)據(jù)的沉降-時間曲線特征變化進(jìn)行綜合研判。信息系統(tǒng)設(shè)置了異常數(shù)據(jù)的處置選項(見表1),供觀測人員在完成現(xiàn)場調(diào)查、補(bǔ)測對比和初步的原因分析后選擇和填報。

2.4 部分觀測階段的參考閾值

根據(jù)1.2節(jié)的方法,對積累的數(shù)據(jù)進(jìn)行數(shù)理統(tǒng)計,制定判識路基、橋涵、隧道沉降異常數(shù)據(jù)的參考閾值,并按照嚴(yán)重程度,由低到高依次設(shè)置成3級,即:關(guān)注、加強(qiáng)關(guān)注、重點(diǎn)關(guān)注。部分觀測階段的建議參數(shù)分別見表2—表4。

表2 判識路基沉降異常數(shù)據(jù)的參考閾值 mm/d

表3 判識橋涵沉降異常數(shù)據(jù)的參考閾值 mm/d

表4 判識隧道沉降異常數(shù)據(jù)的參考閾值 mm/d

2.5 隧道沉降提醒案例

西康高速鐵路白楊坪隧道DK141+020里程,工況為“仰拱(底板)施工完成后,3個月以后”,沉降變形時-程曲線如圖4所示,通過式(4)—式(6)計算,得Vs1=0.187 mm/d,Vs2=0.065 mm/d,最終Vs=0.187 mm/d。經(jīng)過系統(tǒng)自動判識,在第1次隆起提醒的基礎(chǔ)上,觸發(fā)第2次沉降提醒,如圖5所示,對照表4,級別為重點(diǎn)關(guān)注。根據(jù)現(xiàn)場標(biāo)樁狀態(tài)檢查、數(shù)據(jù)分析和外部因素調(diào)查等,基本確定原因后,對照表1,選擇處置選項1(關(guān)注),并繼續(xù)觀察。

圖4 隧道沉降變形時-程曲線

圖5 DK141+020里程斷面自動提醒記錄(部分)

2.6 變形速率代表計算公式靈敏性對比驗證

西延高速鐵路DK189+342里程的路基在2022年12月2日的異常數(shù)據(jù)如圖6所示,工況為堆載預(yù)壓第1～3個月。每組依次取連續(xù)5期的有效觀測數(shù)據(jù),各組分別計算單次變形速率,前5、4、3期的平均變形速率。各組計算結(jié)果繪制時間-速率曲線如圖7所示。

圖6 路基DK189+342處沉降曲線

圖7 最近1—5期的平均速率取值對比

可以看出,前5、4、3期的平均變形速率(依次為V5-1、V4-1、V3-1)中,在12月2日前后,前3期的平均變形速率是最大值,為0.69 mm/d;單次變形速率(為Vs1)最大值為1.03 mm/d;此時,代入靈敏性調(diào)整系數(shù)k=1.5,得1.03

3 結(jié) 語

本文提出了利用平均變形速率判識沉降異常數(shù)據(jù)的方法。在沉降觀測積累的海量數(shù)據(jù)中,應(yīng)用此方法能極大降低人工發(fā)現(xiàn)異常數(shù)據(jù)的工作量,為如實表征高速鐵路主體工程沉降變形情況,搭建了準(zhǔn)確、客觀的數(shù)據(jù)基礎(chǔ)環(huán)境。

通過設(shè)置靈敏性調(diào)整系數(shù)k,可以抑制判識方法的誤判率,隨工程應(yīng)用和數(shù)據(jù)積累,相應(yīng)k的取值將進(jìn)一步得到優(yōu)化調(diào)整,判識路基沉降異常數(shù)據(jù)的參考閾值也同樣會得到修正,以此可進(jìn)一步提高判識的準(zhǔn)確率。

通過西康、西延等高速鐵路沉降觀測應(yīng)用驗證,面向全國大規(guī)模的鐵路建設(shè),以數(shù)據(jù)驅(qū)動沉降觀測的精細(xì)化管理,在自動、精準(zhǔn)甄別數(shù)據(jù)異常、精準(zhǔn)評估工程沉降穩(wěn)定性方面,還將發(fā)揮出巨大的數(shù)據(jù)資源優(yōu)勢。