基于灰色預(yù)測的空間多點殘差修正模型研究＊

李曉蕾 劉 睿 田永瑞 謝炳科

(1)重慶科技學(xué)院建筑工程學(xué)院,重慶 401331 2)中國科學(xué)研究院地理科學(xué)與資源研究所,資源與環(huán)境信息系統(tǒng)國家重點實驗室,北京 100101 3)長安大學(xué)地質(zhì)工程與測繪學(xué)院,西安 710054)

基于灰色預(yù)測的空間多點殘差修正模型研究＊

李曉蕾1,3)劉 睿2)田永瑞3)謝炳科1)

(1)重慶科技學(xué)院建筑工程學(xué)院,重慶 401331 2)中國科學(xué)研究院地理科學(xué)與資源研究所,資源與環(huán)境信息系統(tǒng)國家重點實驗室,北京 100101 3)長安大學(xué)地質(zhì)工程與測繪學(xué)院,西安 710054)

在分析灰色單點預(yù)測模型不足的基礎(chǔ)之上,綜合考慮變形點之間的空間相關(guān)性,提出將殘差修正應(yīng)用于空間多點模型,引入空間多點殘差修正模型。對黃河小浪底水利樞紐工程 B斷面沉降監(jiān)測結(jié)果的分析表明,該模型整體模擬與預(yù)測相對誤差絕對值之和 (5個周期)為 122%,優(yōu)于空間多點模型的 179%與 G M(1,1)模型的284%。同時,短周期數(shù)據(jù)模擬(5個周期內(nèi))相對誤差能夠控制在 20%以內(nèi),具有最好的數(shù)據(jù)擬合精度。

灰色預(yù)測;變形分析;殘差修正;空間多點殘差修正模型;沉降監(jiān)測

1 引言

變形分析與預(yù)報的常規(guī)方法主要包括回歸分析、時間序列模型、Kalman濾波模型等,但以上方法適用于觀測數(shù)列較長,數(shù)據(jù)信息量大的數(shù)據(jù)序列,而現(xiàn)代大量工程中經(jīng)常存在短數(shù)據(jù)貧信息序列,信息量少、規(guī)律性不強,使得上述方法的應(yīng)用存在較大難度。所以,自從鄧聚龍教授提出灰色系統(tǒng)理論以來,由于灰色系統(tǒng)理論提供了貧信息情況下解決系統(tǒng)問題的新途徑,為變形監(jiān)測領(lǐng)域?qū)崟r、動態(tài)、準確地開展變形分析與預(yù)報提供了新的手段與方法,因此,灰色系統(tǒng)理論得到了飛速發(fā)展[1,2]。

但應(yīng)用中,灰色單點預(yù)測模型準確率不高,嚴重影響預(yù)測結(jié)果,因此,本文針對該問題,在分析灰色單點預(yù)測模型的不足基礎(chǔ)之上,綜合考慮變形點之間的空間相關(guān)性,引入空間多點殘差修正模型,并結(jié)合實際工程進行比較分析,得到了比較好的結(jié)果。

2 灰色預(yù)測模型原理

灰色預(yù)測是對本特性灰色系統(tǒng)進行的預(yù)測,利用已有的數(shù)據(jù)資料,建立灰色系統(tǒng)模型,對系統(tǒng)未來的發(fā)展作趨勢外推,其預(yù)測值為一個灰平面區(qū)間,且區(qū)間大小與預(yù)測精度成反比 (圖 1)[3,4]。基于灰色系統(tǒng)理論的灰色預(yù)測對原始數(shù)據(jù)沒有大樣本的要求,通常只要原始數(shù)列有 4個以上的數(shù)據(jù),就可以通過數(shù)據(jù)生成變換建立灰色模型 G M(1,1)[5]。

圖1 灰色模塊示意圖Fig.1 Sketch of greymodule

2.1 GM(1,1)模型

若原始數(shù)列 x(0)={x(0)(k)}的一次累加生成(1-AGO)數(shù)列 x(1)={x(1)(k)}具有明顯的指數(shù)規(guī)律,則可以建立 G M(1,1)模型：

其中:a稱為灰色系統(tǒng)的發(fā)展系數(shù);u稱作灰作用量,是一個內(nèi)生變量,可以理解為 t時刻,變量 x(t)的變化率 dx(t)/dt與變量本身 x(t)及灰作用量 u的一個線性組合,u是常量,其通解為:

基于時間點 t=0,1,2,…,k,將式 (2)離散化,有:

式(3)為 G M(1,1)模型的離散響應(yīng)函數(shù),當 k =0,1,2,…,n-1時,得到的(1)(k)為 G M(1,1)模型的擬合值;當 k=n,n+1,…,時,得到的(1)(k)為 G M(1,1)模型的預(yù)測值。對(1)(k)做一次累減生成(1-I AGO),即可還原得到原始數(shù)列的擬合值與預(yù)測值(0)(k),其建模框圖如圖 2所示。

圖 2 G M(1,1)模型建模框圖Fig.2 Block diagram ofmodel G M(1,1)

2.2 殘差 GM(1,1)模型

式(6)稱為殘差修正 G M(1,1)模型。

2.3 灰色預(yù)測模型改進思路

G M(1,1)、G M(1,1)以及本文中沒有詳細列出的 G M(1,N)、G M(2,1)[6]等灰色模型其實質(zhì)都是針對單點變形的研究,其中 G M(1,N)考慮到了相關(guān)點之間的影響,但其仍然是基于單點數(shù)據(jù)進行預(yù)測,沒有從整體考慮變形點之間的相互關(guān)系。而實際各種工程的變形中,大量的變形監(jiān)測點間相互影響、彼此關(guān)聯(lián)。拋開這些相互關(guān)聯(lián)的信息,僅僅對單點進行處理與預(yù)測,顯然不足以反映變形體的整體變形趨勢和變形規(guī)律。因此,有必要轉(zhuǎn)向多點的空間整體變形分析,充分利用各類采集到的信息,發(fā)揮整體積極作用,更接近現(xiàn)實本質(zhì)的空間動態(tài)分析,完整和全面地反映變形體的狀況。

3 空間多點殘差修正模型研究

3.1 空間多點模型

考慮 n個點相互關(guān)聯(lián)和相互影響[7],對此生成序列建立 n元一階常微分方程組:

寫成矩陣形式:

其連續(xù)時間響應(yīng)為[8]：

3.2 空間多點殘差修正模型

當 k＜m時,將式(9)離散化得到：

4)將滿足殘差修正條件的點建立 1-AGO序列代入式 (5)計算其模擬時間序列(k);

5)根據(jù)式(6)得到其修正后的時間響應(yīng),得出殘差修正值;

6)模型精度評定。

4 模型試驗與分析

選擇黃河小浪底水利樞紐工程大壩B斷面的 4個沉降監(jiān)測數(shù)據(jù)點 2001—2006年垂直方向位移年變化量作為試驗數(shù)據(jù)(表 1)。以Matlab 7為建模工作環(huán)境,綜合運用 G M(1,1)、G M(1,N)、空間多點以及空間多點殘差修正模型對實測值數(shù)據(jù)進行模擬與預(yù)測。其中前 5個周期數(shù)據(jù)用于模擬,為模擬值,最后一個周期數(shù)據(jù)為預(yù)測值 G M(1,N)無預(yù)測值,全部為模擬值)。經(jīng)計算,C13與 813點的原始序列不存在 k0滿足殘差修正模型建立的條件,因此文中僅209與B14點存在空間多點灰色模型修正值。數(shù)據(jù)結(jié)果與對比如表 1與圖 3所示。

表 1 小浪底大壩 B斷面年際沉降監(jiān)測數(shù)據(jù)(單位:mm)Tab.1 Annual subsidence mon itoring data of B section of the Xi aolangdiDam(un it:mm)

從表 1與圖 3可知：

1)空間多點殘差修正模型的數(shù)據(jù)擬合精度最高,兩個點相對誤差絕對值之和分別為 38%與15%。同時,相對誤差呈現(xiàn)出隨周期逐漸增加的趨勢,短期內(nèi)模擬精度最好。

2)G M(1,N)模型的數(shù)據(jù)擬合精度次于空間多點殘差修正模型,變化趨勢與空間多點殘差修正模型相反,短期內(nèi)擬合相對誤差較大,更適合于長期數(shù)據(jù)模擬與內(nèi)插。然其不能應(yīng)用于數(shù)據(jù)預(yù)測,局限性較大。

3)3種預(yù)測模型中,空間多點模型的預(yù)測效果最好,其平均相對誤差為 96%。空間多點殘差修正模型次之,兩個點的相對誤差平均值為 96%,不如空間多點模型的 39%。G M(1,1)模型的預(yù)測精度最差,其相對誤差平均值達到了 284%。

4)綜合考慮整體模擬與預(yù)測精度,空間多點殘差修正模型的精度最高,相對誤差絕對值之和為122%,空間多點模型為 179%,均遠遠優(yōu)于 G M(1, 1)模型的 284%。

圖3 預(yù)測結(jié)果對比Fig.3 Comparison among the prediction results

5 結(jié)論與展望

空間多點殘差修正模型的整體模擬預(yù)測精度在現(xiàn)有幾種灰色系統(tǒng)模型中最好,其整體相對誤差絕對值之和為 122%,優(yōu)于空間多點模型的 179%與G M(1,1)模型的 284%。同時,該模型能夠達到很好的數(shù)據(jù)擬合精度,尤其對于短期數(shù)據(jù)模擬(5個周期內(nèi))相對誤差能夠控制在 20%以內(nèi)。空間多點模型數(shù)據(jù)預(yù)測精度略優(yōu)于空間多點殘差修正模型,然而模擬精度遠不如空間多點殘差修正模型與 G M (1,N)模型。總體而言,多點模型將變形過程作為一個整體的系統(tǒng)進行考慮,該思想代表了變形分析與預(yù)報的發(fā)展趨勢。

空間多點殘差修正模型要求數(shù)據(jù)序列滿足殘差修正條件,在適用范圍方面不如空間多點模型。然而空間多點模型在數(shù)據(jù)擬合與內(nèi)插精度上不如 G M (1,N)模型。具體到實際應(yīng)用中需要多數(shù)據(jù)多模型綜合應(yīng)用以解決變形監(jiān)測中的實際問題。同時,考慮到文中數(shù)據(jù)量以及數(shù)據(jù)為一次殘差修正等限制,作者認為空間多點殘差修正模型能夠得到更進一步的發(fā)展應(yīng)用。

1 黃聲享,尹暉,蔣征.變形監(jiān)測數(shù)據(jù)處理[M].武漢：武漢大學(xué)出版社,2003.(Huang Shengxiang,Yin Hui and Jiang Zheng.Defor mation monitoring data processing[M]. Wuhan：Wuhan University Press,2003)

2 肖新平.灰色系統(tǒng)模型方法的研究[D].華中科技大學(xué), 2002.(Xiao Xinping.Study of grey system model[D].Wuhan：HuazhongUniversity of Scienceamp;Technology,2002)

3 鄧聚龍.灰色系統(tǒng)理論教程[M].武漢：華中理工大學(xué)出版社,1990.(Deng Julong.Tutorial of grey system theory [M].Wuhan：Publishing CompanyofHuazhongUniversityof Science and Technology,1990)

4 劉思峰,黨耀國,方志耕.灰色系統(tǒng)理論及其應(yīng)用[M].北京：科學(xué)出版社,2004.(Liu Sifeng,Dang Yaoguo and Fang Zhigeng.Gray system theory and application[M].Beijing：Science Press,2004)

5 劉思峰,鄧聚龍.G M(1,1)模型的適用范圍[J].系統(tǒng)工程理論與實踐,2000,20(5)：121-124.(Liu Sifeng and Deng Julong.The range suitable for G M(1,1)[J].Systems Engineering Theoryamp;Practice,2000,20(5)：121-124)

6 陳曉斌,張家生,安關(guān)峰.G M(1,1)與 G M(2,1)模型在基坑工程預(yù)測中的應(yīng)用 [J].巖土工程學(xué)報,2006,28 (S1)：1 401-1 405.(Chen Xiaobin,Zhang Jiasheng and An Guanfeng.Application of gray-method models G M(1,1) and G M(1,2)forecasting of pit excavation engineering[J]. Chinese Journal of Geotechnical Engineering,2006,28 (S1)：1 401-1 405)

7 潘國榮,劉大杰.顧及鄰近點變形因素項的動態(tài)模型辨識及預(yù)測[J].測繪學(xué)報,2001,30(01)：32-35.(Pan Guorong and Liu Dajie.Dynamic modeling identification and predication in consideration of the adjacent point deformation [J].Acta Geodaetica etCartographica Sinica,2001,30(1)：32-35)

8 王穗輝,潘國榮.基于MATLAB多變量灰色模型及其在變形預(yù)測中的應(yīng)用 [J].土木工程學(xué)報,2005,38(05)：24-27.(Wang Suihui and Pan Guorong.A Matlab-based multi-variate grey model and its application in deformation prediction[J].China Civil Engineering Journal,2005,38 (5)：24-27)

STUDY ON SPATIAL M ULTI-PO INT RESIDUAL MODEL BASED ON GREY PRED ICTION

Li Xiaolei1,3),Liu Rui2),Tian Yongrui3)and Xie Bingke1)

(1)School of Civil Engineering and Architecture,Chongqing University of Science and Technology,Chongqing 401331 2)Institute of Geographical Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101 3)College of Geology Engineering and Geom atics,Chang’an University,X i’an 710054)

On the basisof the analysisof the shortage of Grey single-pointpredictionmodel,by considering the spatial correlation between the deformation observation points,the residual modification is applied in the spatial multi-pointmodel and thus the spatialmulti-point residualmodel is developed.By taking the subsidence monitoring data of Xiaolangdimultipurpose dam as an example,the results shows that the total absolute values of relative error (5 cycles)of spatialmulti-point residualmodel is122%,better than 179%of spatialmulti-pointmodel and 284% of G M(1,1)model.Meanwhile,thismodel has the best fitting accuracy with the relative error less than 20%in short-cycle fitting.

grey prediction;defor mation analysis;residual modification;spatial multi-point residual model;subsidence monitoring

1671-5942(2010)05-0125-04

2010-04-19

中國科學(xué)院地理科學(xué)與資源研究所前沿領(lǐng)域創(chuàng)新項目

李曉蕾,女,1981年生,碩士,主要從事建筑工程測量理論與應(yīng)用等研究.E-mail：li.xiaolei8@gmail.com

P207

A