鄧 銳 李 波 羅海風(fēng)

(中鐵西南科學(xué)研究院有限公司， 四川成都 610000)

凍土地區(qū)高速鐵路軟土路基沉降預(yù)測(cè)模型研究

鄧 銳 李 波 羅海風(fēng)

(中鐵西南科學(xué)研究院有限公司， 四川成都 610000)

對(duì)幾種常用沉降預(yù)測(cè)模型進(jìn)行了對(duì)比分析，介紹了組合預(yù)測(cè)模型的構(gòu)建思路及預(yù)測(cè)效果的評(píng)價(jià)標(biāo)準(zhǔn)，建立了凍土地區(qū)路基組合預(yù)測(cè)模型。

凍土地區(qū) 軟土路基 沉降預(yù)測(cè) 組合預(yù)測(cè)模型

1 常用沉降預(yù)測(cè)精度對(duì)比

沉降觀測(cè)工作是測(cè)量中變形觀測(cè)最基本的工作之一。隨著科技的不斷進(jìn)步和眾多專家學(xué)者的不斷研究與探索，目前已建立了比較完善的沉降觀測(cè)體系。常用的沉降預(yù)測(cè)模型主要有以下幾種：規(guī)范雙曲線模型、修正雙曲線模型、泊松曲線模型、三點(diǎn)法模型、Asaoka法模型和灰色系統(tǒng)GM(1，1)模型[6-7]，預(yù)測(cè)精度的評(píng)價(jià)主要有以下兩個(gè)方面：相關(guān)系數(shù)分析與預(yù)測(cè)誤差分析，根據(jù)“哈齊客運(yùn)專線”某區(qū)段實(shí)測(cè)的數(shù)據(jù)，對(duì)以上模型的預(yù)測(cè)精度進(jìn)行分析。

1.1 相關(guān)系數(shù)分析

相關(guān)系數(shù)是高速鐵路沉降觀測(cè)評(píng)估分析一個(gè)很重要的技術(shù)指標(biāo)，也是檢查沉降觀測(cè)數(shù)據(jù)質(zhì)量的重要手段[8]。本文采取若干個(gè)觀測(cè)斷面的沉降觀測(cè)數(shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù)，通過(guò)分析并計(jì)算出各模型的相關(guān)系數(shù)。規(guī)范雙曲線法的曲線回歸的相關(guān)系數(shù)為64.4%；修正雙曲線法的曲線回歸的相關(guān)系數(shù)為63.2%；泊松曲線法的曲線回歸的相關(guān)系數(shù)為95.4%；三點(diǎn)法的曲線回歸的相關(guān)系數(shù)為93.1%；Asaoka法的曲線回歸的相關(guān)系數(shù)為89.7%；灰色系統(tǒng)GM(1，1)模型曲線回歸的相關(guān)系數(shù)為40.2%。由此可見(jiàn)，泊松曲線法、三點(diǎn)法和Asaoka法曲線回歸的相關(guān)系數(shù)都比較大，規(guī)范雙曲線法和修正雙曲線法的相關(guān)系數(shù)次之，灰色系統(tǒng)GM(1，1)模型的相關(guān)系數(shù)較低。

1.2 預(yù)測(cè)誤差分析

預(yù)測(cè)誤差也是高速鐵路沉降觀測(cè)評(píng)估分析一個(gè)很重要的技術(shù)指標(biāo)。它是檢驗(yàn)沉降預(yù)測(cè)質(zhì)量的重要手段，預(yù)測(cè)誤差的大小將直接影響路基的工后沉降量的判定。表1列出了上述監(jiān)測(cè)斷面的各個(gè)預(yù)測(cè)模型曲線擬合的預(yù)測(cè)誤差范圍(△指預(yù)測(cè)誤差)。

表1 各個(gè)預(yù)測(cè)方法曲線回歸的預(yù)測(cè)誤差統(tǒng)計(jì) %

由表1可見(jiàn)，對(duì)該段測(cè)區(qū)而言，泊松曲線法和Asaoka法曲線回歸的預(yù)測(cè)誤差比較??；規(guī)范雙曲線法、修正雙曲線法和三點(diǎn)法的預(yù)測(cè)誤差次之；灰色系統(tǒng)GM(1，1)模型的預(yù)測(cè)誤差最大。

通過(guò)對(duì)以上預(yù)測(cè)模型相關(guān)系數(shù)和預(yù)測(cè)誤差進(jìn)行綜合比較，可知泊松曲線法和Asaoka法在該段區(qū)域內(nèi)的預(yù)測(cè)效果優(yōu)于其他幾種預(yù)測(cè)模型。

2 組合預(yù)測(cè)模型介紹

現(xiàn)實(shí)生活中，一個(gè)問(wèn)題常常有不同的解決方法，高速鐵路軟土路基沉降預(yù)測(cè)也一樣，有很多種預(yù)測(cè)模型可供選擇使用，但是各個(gè)預(yù)測(cè)方法將給出不相同的預(yù)測(cè)結(jié)果，預(yù)測(cè)精度也不盡相同。如果按照預(yù)測(cè)精度的好壞，從多種預(yù)測(cè)方法中選擇預(yù)測(cè)精度最好的一個(gè)，而把預(yù)測(cè)誤差較大的預(yù)測(cè)方法完全拋棄，這樣并不是提高預(yù)測(cè)精度的辦法，反而會(huì)失去一些有用的信息，造成信息的浪費(fèi)。比較科學(xué)的方法是：把各個(gè)不同的預(yù)測(cè)方法通過(guò)某種方式進(jìn)行組合搭配，也就是組合預(yù)測(cè)模型。組合預(yù)測(cè)方法就是針對(duì)同一個(gè)問(wèn)題，綜合分析各種常用預(yù)測(cè)模型的預(yù)測(cè)信息，把兩種或兩種以上的預(yù)測(cè)模型按照適當(dāng)?shù)臋?quán)重比例組合起來(lái)，從而得到組合預(yù)測(cè)模型；組合預(yù)測(cè)方法的關(guān)鍵問(wèn)題是各個(gè)常用預(yù)測(cè)模型加權(quán)系數(shù)的確定，要使得組合預(yù)測(cè)方法達(dá)到最好的預(yù)測(cè)效果；組合的主要目的是綜合各個(gè)模型的優(yōu)點(diǎn)，避其不足，盡可能地提高預(yù)測(cè)精度。一般地，組合預(yù)測(cè)方法有兩種基本形式。

等權(quán)組合：將各個(gè)預(yù)測(cè)模型的預(yù)測(cè)值按相同的權(quán)系數(shù)組合，得到新的預(yù)測(cè)值。

不等權(quán)組合：給各個(gè)預(yù)測(cè)模型的預(yù)測(cè)值賦予的權(quán)系數(shù)不相同，組合成新的預(yù)測(cè)值。

2.1 組合預(yù)測(cè)模型的構(gòu)建思路

如果一個(gè)單預(yù)測(cè)模型含有待預(yù)測(cè)系統(tǒng)的獨(dú)立信息，即使它的預(yù)測(cè)精度比較差，但是當(dāng)它與另一個(gè)較好的單預(yù)測(cè)模型組合后，就會(huì)有可能提高組合模型的預(yù)測(cè)精度。其數(shù)學(xué)表達(dá)形式如下

(t=1，2，…，n)

上式被稱為B-G模型。其中，m為單預(yù)測(cè)模型的個(gè)數(shù)；n為實(shí)測(cè)數(shù)據(jù)的個(gè)數(shù)；fit(i=1，2，…，m;t=1，2，…，n)為第i種單預(yù)測(cè)模型在t時(shí)刻的預(yù)測(cè)值，ki(i=1，2，…，m)為第i種單預(yù)測(cè)模型在線性組合預(yù)測(cè)模型中的權(quán)重，為保證組合預(yù)測(cè)的無(wú)偏性，各個(gè)權(quán)重系數(shù)需滿足：k1+k2+…+km=1，yt為線性組合預(yù)測(cè)模型在t時(shí)刻的組合預(yù)測(cè)值。

組合預(yù)測(cè)模型的一般建模步驟如圖1所示。

圖1 組合預(yù)測(cè)模型的建模步驟

2.2 組合預(yù)測(cè)模型運(yùn)用需要考慮的問(wèn)題

(1)所選各個(gè)單預(yù)測(cè)模型都能獨(dú)立對(duì)觀測(cè)數(shù)據(jù)進(jìn)行預(yù)測(cè)分析，并能達(dá)到一定的預(yù)測(cè)效果。

(2)所選單預(yù)測(cè)模型能夠滿足組合預(yù)測(cè)模型的應(yīng)用條件。

(3)所選單預(yù)測(cè)模型對(duì)觀測(cè)數(shù)據(jù)的要求或假設(shè)條件應(yīng)大致相同。

(4)所選單預(yù)測(cè)模型的復(fù)雜程度和數(shù)據(jù)處理的難易程度要適度。

(5)單預(yù)測(cè)模型個(gè)數(shù)的確定要合理。

(6)選擇合適的權(quán)系數(shù)求解方法。

2.3 組合預(yù)測(cè)模型預(yù)測(cè)效果的評(píng)價(jià)標(biāo)準(zhǔn)

根據(jù)預(yù)測(cè)效果的評(píng)價(jià)原則和慣例，常采用下面五個(gè)指標(biāo)來(lái)評(píng)價(jià)組合預(yù)測(cè)模型的預(yù)測(cè)效果：

(1)平方和誤差

(2)均方誤差

(3)平均絕對(duì)誤差

(4)均方百分比誤差

(5)平均絕對(duì)百分比誤差

一般情況下，通過(guò)滿足“預(yù)測(cè)誤差平方和最小準(zhǔn)則”和“預(yù)測(cè)誤差絕對(duì)值最小準(zhǔn)則”來(lái)衡量模型的預(yù)測(cè)精度。

預(yù)測(cè)誤差平方和最小準(zhǔn)則：假定某一預(yù)測(cè)對(duì)象的指標(biāo)序列{xt，t=1，2，…，N}，存在m種單預(yù)測(cè)模型對(duì)其進(jìn)行預(yù)測(cè)分析，設(shè)第i種單預(yù)測(cè)模型在t時(shí)刻的預(yù)測(cè)值為xit，i=1，2，…，m；t=1，2，…，N，稱eit=(xt-xit)為第i種單預(yù)測(cè)模型在t時(shí)的預(yù)測(cè)誤差。

設(shè)J1為組合預(yù)測(cè)模型預(yù)測(cè)誤差平方和，其表達(dá)式為

因此，以預(yù)測(cè)誤差平方和最小準(zhǔn)則的線性組合預(yù)測(cè)模型為下面方程式的最優(yōu)化問(wèn)題

假設(shè)m(m

上面的公式可以轉(zhuǎn)化為以下矩陣形式

模型的最優(yōu)解和目標(biāo)函數(shù)最優(yōu)解

預(yù)測(cè)誤差絕對(duì)值和最小準(zhǔn)則：

假設(shè)F表示組合預(yù)測(cè)模型的預(yù)測(cè)誤差絕對(duì)值和，則有

可以得到以誤差絕對(duì)值和最小準(zhǔn)則的最優(yōu)線性組合預(yù)測(cè)模型

則可以將預(yù)測(cè)問(wèn)題轉(zhuǎn)換為線性規(guī)劃求解最優(yōu)解的問(wèn)題，假設(shè)

；

這樣就轉(zhuǎn)換成一個(gè)線性規(guī)劃問(wèn)題，用現(xiàn)有的線性規(guī)劃軟件即可求解。

3 凍土地區(qū)路基組合預(yù)測(cè)模型的建立

結(jié)合凍土地區(qū)的實(shí)驗(yàn)數(shù)據(jù)，論證三點(diǎn)法和灰色系統(tǒng)GM(1，1)的組合預(yù)測(cè)模型在滿足預(yù)測(cè)誤差平方和最小準(zhǔn)則和以預(yù)測(cè)誤差絕對(duì)值和最小準(zhǔn)則前提下的可行性，并與其他基本預(yù)測(cè)模型進(jìn)行精度對(duì)比。

3.1 相關(guān)系數(shù)對(duì)比

這里仍然以前面所選測(cè)區(qū)的觀測(cè)數(shù)據(jù)作為考查對(duì)象，從誤差平方和最小準(zhǔn)則和誤差絕對(duì)值和最小準(zhǔn)則兩個(gè)方面考慮，將三點(diǎn)法-灰色系統(tǒng)GM(1，1)組合預(yù)測(cè)模型與泊松曲線法、Asaoka法、三點(diǎn)法、灰色系統(tǒng)GM(1，1)模型的預(yù)測(cè)結(jié)果相比較，表2給出了6種不同預(yù)測(cè)方法下的曲線回歸分析的相關(guān)系數(shù)大小。

分析表2可以得出，除了灰色系統(tǒng)GM(1，1)模型曲線回歸的相關(guān)系數(shù)不能完全滿足《評(píng)估指南》中相關(guān)系數(shù)不小于0.92的要求之外，其他幾種預(yù)測(cè)模型曲線回歸的相關(guān)系數(shù)都能滿足要求；以誤差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法比以誤差絕對(duì)值和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法的曲線回歸的相關(guān)系數(shù)大，其中以誤差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法比三點(diǎn)法和灰色系統(tǒng)GM(1，1)模型的曲線回歸的相關(guān)系數(shù)大，以誤差絕對(duì)值和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法與三點(diǎn)法的相關(guān)系數(shù)相當(dāng)，比灰色系統(tǒng)GM(1，1)模型的相關(guān)系數(shù)大。

表2 5個(gè)斷面在6種不同預(yù)測(cè)方法下的曲線回歸分析的相關(guān)系數(shù)

3.2 預(yù)測(cè)誤差對(duì)比

圖2為某一斷面在6種不同預(yù)測(cè)方法下的預(yù)測(cè)誤差。

圖2 DK12+550斷面在6種不同預(yù)測(cè)方法下的預(yù)測(cè)誤差

分析上圖可知，除了灰色系統(tǒng)GM(1，1)模型的預(yù)測(cè)誤差較大之外，其他幾種預(yù)測(cè)模型的預(yù)測(cè)誤差相對(duì)較?。灰哉`差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法、Asaoka法、泊松曲線法的預(yù)測(cè)誤差最小，以誤差絕對(duì)值和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法和三點(diǎn)法的預(yù)測(cè)誤差次之。

4 結(jié)論

從以上的介紹和對(duì)實(shí)測(cè)數(shù)據(jù)的分析，可以得出如下基本結(jié)論：

從曲線回歸的相關(guān)系數(shù)大小和預(yù)測(cè)誤差大小來(lái)看，以誤差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法的預(yù)測(cè)效果明顯優(yōu)于三點(diǎn)法與灰色系統(tǒng)GM(1，1)模型；以誤差絕對(duì)值和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法的預(yù)測(cè)效果比灰色系統(tǒng)GM(1，1)模型好很多，但與三點(diǎn)法的預(yù)測(cè)效果相當(dāng)。

從曲線回歸的相關(guān)系數(shù)大小和預(yù)測(cè)誤差大小來(lái)看，以誤差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法比以誤差絕對(duì)值和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型組合預(yù)測(cè)方法的預(yù)測(cè)效果要好，與泊松曲線法、Asaoka法的預(yù)測(cè)效果基本相同。

對(duì)于哈齊鐵路客運(yùn)專線軟土路基“沉降變形量級(jí)小、相對(duì)波動(dòng)大”的情形而言，推薦優(yōu)先考慮全過(guò)程的預(yù)測(cè)精度高的沉降預(yù)測(cè)方案，即首選以誤差平方和最小的三點(diǎn)法-灰色系統(tǒng)GM(1，1)模型的組合預(yù)測(cè)方法和泊松曲線法對(duì)沉降全過(guò)程進(jìn)行預(yù)測(cè)分析，其次選擇Asaoka法對(duì)恒載后的沉降觀測(cè)數(shù)據(jù)進(jìn)行預(yù)測(cè)分析，并和前兩種預(yù)測(cè)結(jié)果進(jìn)行對(duì)比。也可以運(yùn)用三點(diǎn)法、修正雙曲線法等預(yù)測(cè)方法進(jìn)行預(yù)測(cè)，綜合考慮各個(gè)預(yù)測(cè)結(jié)果，這樣會(huì)有更好的效果。

本文基于凍土地區(qū)高速鐵路軟土路基的沉降進(jìn)行預(yù)測(cè)模型的選取，對(duì)凍土地區(qū)其他軟土路基的沉降預(yù)測(cè)模型的選取具有一定的參考價(jià)值。

[1] 張麗麗.軟土路基沉降計(jì)算及預(yù)測(cè)方法研究[D].西安:西安建筑科技大學(xué)，2006

[2] 劉升傳，曹淵.新建鐵路軟土路基沉降規(guī)律研究[J].鐵道工程學(xué)報(bào)，2010，5(5):1-3

[3] 李濤，張儀萍，張土喬.軟土路基沉降的優(yōu)性組合預(yù)測(cè)[J].巖石力學(xué)與工程學(xué)報(bào)，2005，24(18):3282-3286

[4] 雷康.軟土路基沉降的組合預(yù)測(cè)方法研究[D].長(zhǎng)沙:長(zhǎng)沙理工大學(xué)，2010

[5] 孟慶新.軟土路基變形分析與沉降預(yù)測(cè)[D].石家莊:河北工業(yè)大學(xué)，2006

[6] 趙明華，劉煜，曹文貴.軟土路基沉降發(fā)展規(guī)律及其預(yù)測(cè)[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，35(1):157-161

[7] 許俊超.軟土地基沉降預(yù)測(cè)方法探討[J].鐵道勘測(cè)與設(shè)計(jì)，2008(4)

[8] 中華人民共和國(guó)鐵道部.客運(yùn)專線無(wú)砟軌道鐵路工程測(cè)量暫行規(guī)定[S].北京:中國(guó)鐵道出版社，2006

[9] 何飛.高速鐵路路基基礎(chǔ)工后沉降預(yù)測(cè)的相關(guān)研究[D].成都:西南交通大學(xué)，2009

[10]劉志明.客運(yùn)專線地基沉降計(jì)算、預(yù)測(cè)方法研究[D].成都:西南交通大學(xué)，2009

ResearchofPredictionModelofHighspeedRailwaySoftSoilSubgradeSettlementinFrozenEarthArea

DENG Rui LI Bo LUO Hai-feng

2014-04-24

鄧 銳(1987—)，男，2012年畢業(yè)于中國(guó)地質(zhì)大學(xué)(武漢)遙感科學(xué)與技術(shù)專業(yè)，助理工程師。

1672-7479(2014)04-0024-04

TU433

： A