解志剛 許海東

1 概述

在預(yù)測(cè)分析中，最基本的是線(xiàn)性回歸預(yù)測(cè)，但是在實(shí)際工程應(yīng)用中，我們的數(shù)據(jù)是非線(xiàn)性的、離散狀態(tài)的，隧道圍巖的變形就是一個(gè)由復(fù)雜因素所決定的非線(xiàn)性的變形，對(duì)于隧道圍巖變形的預(yù)測(cè)中線(xiàn)性預(yù)測(cè)辦法是不可取的，我們可以應(yīng)用灰色理論對(duì)圍巖的變形進(jìn)行預(yù)測(cè)?；疑碚摚址Q(chēng)灰色系統(tǒng)，是指既有已知信息，又包含未知信息的系統(tǒng);是研究分析、建模、預(yù)測(cè)、決策、控制的理論。該理論自20世紀(jì)80年代創(chuàng)立以來(lái)，已經(jīng)在社會(huì)科學(xué)和自然科學(xué)的許多領(lǐng)域得到了應(yīng)用。巖石力學(xué)中的一些問(wèn)題，例如上述圍巖變形，具有灰色系統(tǒng)的特點(diǎn)，灰色理論能很好的應(yīng)用于上述問(wèn)題。以下將討論利用灰色理論進(jìn)行圍巖變形的預(yù)測(cè)。

2 灰色理論的GM(1，1)模型

圍巖變形的測(cè)量結(jié)果，是一個(gè)時(shí)間序列數(shù)據(jù)，通過(guò)建立灰色理論模型，來(lái)預(yù)測(cè)未來(lái)圍巖變形，以便為工程的設(shè)計(jì)和施工提供有關(guān)的圍巖動(dòng)態(tài)信息。首先討論GM(1，1)模型問(wèn)題。

2.1 建立 GM(1，1)模型

假定在時(shí)間序列 t1，t2，…，tn，量測(cè)的圍巖變形為數(shù)列 x1，x2，…，xn，記為:

此為原始數(shù)據(jù)，為了弱化原始數(shù)據(jù)的隨機(jī)性，對(duì)原始數(shù)據(jù)做累加，也就是它的1-AGO，得到數(shù)列:

由原始數(shù)據(jù)(1)和累加數(shù)據(jù)(2)可以分別得到向量Y和矩陣B:

對(duì)于GM(1，1)模型，相應(yīng)的白化微分方程為:

該微分方程的解為:

其中，a，b均為待定系數(shù)，一般記為 A=［a，b］T。

對(duì)于微分方程(5)，其待定系數(shù)A=［a，b］T的最小二乘法解為:

將原始數(shù)據(jù)x(0)(i)和累加數(shù)據(jù)x(1)(i)代入式(7)中，可以得到a，b值。我們可以通過(guò)式(6)進(jìn)行預(yù)測(cè)。

2.2 GM(1，1)模型的精度

以上我們介紹了灰色理論模型的建立方法，下面我們要根據(jù)解式(6)計(jì)算殘差:

得到殘差向量和相對(duì)誤差分別為:

設(shè)原始數(shù)據(jù)的方差及殘差向量的方差分別為s21和s22，則:

方差比值Q為:

小誤差概率為:

模型的精度由Q和P來(lái)控制和表述，一般將其分成四級(jí)，見(jiàn)表1。通過(guò)計(jì)算所得P值或Q值來(lái)劃分等級(jí)，評(píng)價(jià)預(yù)測(cè)精度的好壞。

表1 預(yù)測(cè)精度等級(jí)劃分

2.3 基于灰色預(yù)測(cè)的等維灰數(shù)替補(bǔ)迭代

為了進(jìn)一步提高模型預(yù)測(cè)的精度，在灰色模型的基礎(chǔ)上還可以做進(jìn)一步的改進(jìn)，使得預(yù)測(cè)結(jié)果更加準(zhǔn)確。GM(1，1)模型具有預(yù)測(cè)意義的數(shù)據(jù)僅僅是x(0)(i)以后的前幾位數(shù)或者∧x(0)(i)數(shù)列的前幾位數(shù)，隨著時(shí)間的推移，舊的數(shù)據(jù)越來(lái)越不適應(yīng)新的情況，所以，要在原數(shù)據(jù)的基礎(chǔ)上每次增加一個(gè)新信息時(shí)，就去掉一個(gè)舊信息。這種新數(shù)據(jù)補(bǔ)充、舊數(shù)據(jù)除掉的數(shù)列，由于維數(shù)沒(méi)有變，因而叫等維灰數(shù)數(shù)列或者叫新陳代謝數(shù)列。利用新的數(shù)列重新建立GM(1，1)模型，對(duì)下一個(gè)新數(shù)據(jù)進(jìn)行預(yù)測(cè)，預(yù)測(cè)完畢后，再次生成新的數(shù)列。這樣就類(lèi)似于迭代原理，我們可以通過(guò)數(shù)值計(jì)算來(lái)實(shí)現(xiàn)GM(1，1)模型的生成和迭代，這樣的改進(jìn)可以保證精度的良好性。

3 利用GM(1，1)模型預(yù)測(cè)隧道圍巖的變形

以武廣高速鐵路衡陽(yáng)段胡子沖隧道為例，隧道的圍巖以花崗巖為主，地下水豐富，地形比較復(fù)雜。隧道初期支護(hù)的噴射混凝土厚度為25 cm，我們?cè)谒淼篱_(kāi)挖過(guò)程中對(duì)DK091+124斷面拱頂圍巖變形進(jìn)行了監(jiān)測(cè)，監(jiān)測(cè)數(shù)據(jù)見(jiàn)表2。

GM(1，1)模型的建立，是在取得一定數(shù)據(jù)之后進(jìn)行的，建模以后就可以對(duì)未來(lái)變形進(jìn)行數(shù)據(jù)預(yù)測(cè)。由于收斂變形是持續(xù)進(jìn)行的，因此，不斷有新的數(shù)據(jù)充實(shí)變形時(shí)間序列，考慮到現(xiàn)場(chǎng)的條件不斷變化，前期測(cè)量的數(shù)據(jù)對(duì)后期研究圍巖的變形有一定的作用，但是越是接近于預(yù)測(cè)時(shí)間的量測(cè)數(shù)據(jù)越有意義(見(jiàn)表3)。因此，我們采用GM(1，1)迭代模型來(lái)進(jìn)行預(yù)測(cè)，預(yù)測(cè)后我們可以對(duì)預(yù)測(cè)和實(shí)測(cè)數(shù)據(jù)進(jìn)行比較(見(jiàn)圖1)。

圖1 實(shí)測(cè)與預(yù)測(cè)數(shù)據(jù)對(duì)比圖

表2 監(jiān)測(cè)數(shù)據(jù)

表3 GM(1，1)模型的預(yù)測(cè)值

經(jīng)過(guò)對(duì)殘差的計(jì)算，本例中的GM(1，1)模型的精度達(dá)到了Ⅰ級(jí)，所以預(yù)測(cè)的數(shù)據(jù)與實(shí)際數(shù)據(jù)還是比較吻合的。

通過(guò)計(jì)算我們可以判斷在4 d以后圍巖穩(wěn)定，在50 d左右變形加大，這時(shí)胡子沖隧道進(jìn)尺達(dá)到56 m，應(yīng)該進(jìn)行二次支護(hù)或者應(yīng)停止開(kāi)挖，以保證隧道施工的安全。同時(shí)我們使用有限元軟件對(duì)胡子沖隧道進(jìn)行了施工模擬，通過(guò)模擬可以知道在隧道進(jìn)尺60 m左右時(shí)必須進(jìn)行二次支護(hù)，模擬過(guò)程這里不再贅述。從這里可以看出利用灰色理論進(jìn)行隧道圍巖變形的預(yù)測(cè)在實(shí)際工作中是有很大意義的。

4 結(jié)語(yǔ)

本文通過(guò)以胡子沖隧道的圍巖變形預(yù)測(cè)為例，對(duì)灰色理論在圍巖與初期支護(hù)相互作用變形預(yù)測(cè)的應(yīng)用進(jìn)行了闡述?；疑到y(tǒng)理論有著廣泛的應(yīng)用前景，將其用于分析和處理圍巖的變形預(yù)測(cè)是一種行之有效的方法。通過(guò)上述分析和論述，我們可以得出以下結(jié)論:

1)利用GM(1，1)模型對(duì)變形進(jìn)行預(yù)測(cè)方便之處是實(shí)測(cè)數(shù)據(jù)需要很少(灰色理論并不因?yàn)閷?shí)測(cè)數(shù)據(jù)多而變得精確度高)，并且預(yù)測(cè)數(shù)據(jù)的精度高，對(duì)于普通回歸方法則要更多的數(shù)據(jù)來(lái)預(yù)測(cè)，而且灰色理論能對(duì)一些非線(xiàn)性或者不確定的因素進(jìn)行弱化，以提高預(yù)測(cè)的準(zhǔn)確性。2)收斂變形的時(shí)效性動(dòng)態(tài)表明，預(yù)測(cè)變形必須時(shí)時(shí)更新，而灰色理論預(yù)測(cè)恰恰能解決這個(gè)問(wèn)題，這與以往回歸方法不同，這樣能更加有利于隧道施工中的變形監(jiān)控(施工要求時(shí)時(shí)監(jiān)測(cè)圍巖變形，已確定圍巖變形的穩(wěn)定情況)。3)從取得原始數(shù)據(jù)到建立GM(1，1)模型，計(jì)算量還是比較大的，但是我們可以在計(jì)算機(jī)上利用數(shù)值計(jì)算方法輕松完成建模和迭代計(jì)算。因此，用灰色理論來(lái)分析圍巖變形資料，預(yù)測(cè)圍巖收斂變形以及判斷可能失穩(wěn)的前兆及時(shí)間是有工程實(shí)用意義的。

［1］ 張忠亭，陶振宇.洞室圍巖收斂變形的時(shí)間效應(yīng)特性研究［J］.武漢水利水電學(xué)院學(xué)報(bào)，1995(12):89-92.

［2］ 鄧聚龍.灰色系統(tǒng)理論教程［M］.武漢:華中理工大學(xué)出版社，1990.

［3］ 熊合金，徐華中.灰色控制［M］.北京:國(guó)防工業(yè)出版社，2005:9.

［4］ 解 偉.灰色預(yù)測(cè)和灰色關(guān)聯(lián)度在結(jié)構(gòu)可靠度中的應(yīng)用［J］.華北水利水電學(xué)院學(xué)報(bào)，2005(6):20-25.

［5］ 劉法貴.灰色數(shù)學(xué)及其應(yīng)用［M］.開(kāi)封:河南大學(xué)出版社，2002:21-88.