姜慶飛

(江西明澤水利建設(shè)有限公司，江西 上饒 334000)

我國(guó)水資源豐富，技術(shù)可開(kāi)發(fā)總量約為4億kW，主要集中在西南高山峽谷地帶[1]。隨著我國(guó)水利水電工程的蓬勃發(fā)展，越來(lái)越多的大長(zhǎng)水工隧洞正在修建或擬建中，如40km的雅魯藏布江大拐彎處引水隧洞[2]、98.3km的秦嶺特長(zhǎng)輸水隧洞[3]、7.2km的花溪輸水隧洞等[4]。由于隧洞開(kāi)挖打破了巖層原始應(yīng)力狀態(tài)，引起應(yīng)力重新分布，使隧洞圍巖發(fā)生不同程度的沉降和收斂變形，造成失穩(wěn)破壞。因此，在隧洞施工過(guò)程中對(duì)圍巖變形進(jìn)行監(jiān)測(cè)，掌握隧洞的沉降與收斂變形規(guī)律，對(duì)隧洞支護(hù)時(shí)間和支護(hù)措施的選取提供決策依據(jù)，進(jìn)而保障隧洞施工安全。

目前，國(guó)內(nèi)外學(xué)者對(duì)隧道/隧洞的變形規(guī)律及預(yù)測(cè)做了較多研究。其中，回歸分析方法在變形監(jiān)測(cè)分析中使用較為廣泛，包括指數(shù)型函數(shù)、對(duì)數(shù)型函數(shù)和雙曲線函數(shù)等，如張俊中等采用多元回歸分析對(duì)某建筑物沉降進(jìn)行研究[5]，朱興禮采用回歸分析方法對(duì)富水巖溶隧道的不同變形段建立了動(dòng)態(tài)分析模型[6]，夏才初等采用雙曲線函數(shù)對(duì)茶林頂隧道變形數(shù)據(jù)進(jìn)行回歸擬合分析，并取得了較好擬合效果[7]。此外，隨著灰色系統(tǒng)理論的不斷發(fā)展，學(xué)者們逐漸將灰色理論引入隧道等變形預(yù)測(cè)中，何偉將灰色理論應(yīng)用在紅濟(jì)石調(diào)水工程隧洞變形分析中，對(duì)比了GM、DGM和Verhulst 3種模型的預(yù)測(cè)效果[8]；夏才初采用GM(1，1)和GM(2，1)灰色模型對(duì)某隧洞水平收斂值進(jìn)行分析，得出GM(1，1)模型的適用性更好[7]。然而，現(xiàn)有關(guān)于隧道/隧洞的變形規(guī)律及預(yù)測(cè)研究主要集中在單一地質(zhì)工況，且多為隧道變形研究，對(duì)隧洞不同級(jí)別圍巖段的變形規(guī)律差異缺乏探討。因此，本文依托貴州某引水隧洞工程，對(duì)其Ⅳ級(jí)和Ⅴ級(jí)圍巖典型里程段進(jìn)行變形監(jiān)測(cè)，分析其沉降和收斂變形規(guī)律，綜合采用回歸分析和灰色理論對(duì)隧洞沉降和收斂變形數(shù)據(jù)進(jìn)行分析，評(píng)價(jià)2種方法的適用情況。研究結(jié)果將為隧洞的沉降和收斂變形分析及預(yù)測(cè)提供借鑒。

1 工程概況

該引水隧洞位于貴州省境內(nèi)，隧洞全長(zhǎng)8.7km，最大埋深518m，采用鉆爆法施工開(kāi)挖。隧洞地表屬中高山地貌單元，山坡較陡，植被較發(fā)育；該區(qū)域地層結(jié)構(gòu)較為復(fù)雜，揭露的地層從上至下主要有第四系全新統(tǒng)(Q4el+dl)粉質(zhì)黏土，三疊系上統(tǒng)(T2)、下統(tǒng)(T1)灰?guī)r、頁(yè)巖，二疊系上統(tǒng)(P2)、下統(tǒng)(P1)灰?guī)r、頁(yè)巖。隧洞線路所經(jīng)地區(qū)為川黔南北向構(gòu)造帶與北東向構(gòu)造帶交接的復(fù)合部位，隧洞強(qiáng)風(fēng)化灰?guī)r中節(jié)理裂隙較發(fā)育，呈張開(kāi)-微張狀。

為掌握該隧洞的變形收斂規(guī)律，對(duì)該隧洞主要里程段圍巖沉降變形和收斂變形進(jìn)行了監(jiān)測(cè)。本文主要目的是分析不同圍巖級(jí)別段的變形收斂規(guī)律，故選取了Ⅴ級(jí)和Ⅳ級(jí)典型里程段的監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析。其中，隧洞里程段108+250～108+455圍巖級(jí)別為Ⅴ級(jí)，巖體完整性：破碎-較破碎，穩(wěn)定性：差-較差；里程段109+010～109+415圍巖級(jí)別為Ⅳ級(jí)，巖體較完整，穩(wěn)定性較好。

2 隧洞不同圍巖級(jí)別段沉降收斂規(guī)律

2.1 隧洞變形監(jiān)測(cè)方法

在隧洞沿線上間隔一定距離布置1個(gè)監(jiān)測(cè)斷面，其中Ⅴ級(jí)圍巖間隔15m布置1個(gè)監(jiān)測(cè)斷面，Ⅳ級(jí)圍巖間隔25m布置1個(gè)監(jiān)測(cè)斷面。隧洞沉降測(cè)點(diǎn)和收斂測(cè)線布置示意如圖1所示，每一監(jiān)測(cè)斷面中布置3個(gè)測(cè)點(diǎn)(拱頂B點(diǎn)、拱肩A點(diǎn)和C點(diǎn))監(jiān)測(cè)隧洞沉降變形；在隧洞周邊布置2條測(cè)線，即測(cè)試1-2和測(cè)線3-4，目的是監(jiān)測(cè)隧洞的水平收斂變形。隧洞變形監(jiān)測(cè)頻率為每天監(jiān)測(cè)1次。

圖1 隧洞沉降測(cè)點(diǎn)和收斂測(cè)線布置示意圖

2.2 隧洞沉降和收斂變形監(jiān)測(cè)結(jié)果分析

本文分別選取隧洞Ⅴ級(jí)圍巖和Ⅳ級(jí)圍巖某一斷面的監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析。Ⅴ級(jí)和Ⅳ級(jí)圍巖的拱頂(拱肩)沉降和周邊水平收斂監(jiān)測(cè)結(jié)果如圖2—3所示。由圖可知，對(duì)于2種級(jí)別的隧洞圍巖，隧洞累積沉降值均隨時(shí)間呈先快速增加后緩慢增長(zhǎng)趨勢(shì)；在監(jiān)測(cè)初期隧洞沉降速率快，在前5d內(nèi)幾乎完成了累積沉降值的50%以上，沉降速率達(dá)到了1.34～2.06mm/d(Ⅴ級(jí)圍巖)和0.74～1.10mm/d(Ⅳ級(jí)圍巖)；在監(jiān)測(cè)時(shí)間為1個(gè)月時(shí)，隧洞沉降趨于穩(wěn)定。此外，拱頂?shù)某两得黠@高于拱肩，測(cè)點(diǎn)A和測(cè)點(diǎn)C的沉降值較為接近，主要是因?yàn)?測(cè)點(diǎn)為對(duì)稱分布。隧洞累積水平收斂值隨時(shí)間的增加而增加，但收斂變形速率逐漸降低，整體上比沉降值的變化更為緩慢。測(cè)線1-2的收斂值略高于測(cè)線3-4，差別不明顯，Ⅴ級(jí)和Ⅳ級(jí)圍巖2條測(cè)線的收斂值差異均小于1mm。對(duì)隧洞不同級(jí)別圍巖而言，Ⅴ級(jí)圍巖的沉降量和收斂量均大于Ⅳ級(jí)圍巖，隧洞拱頂和拱肩趨于穩(wěn)定時(shí)的沉降值分別為18.2mm和14.6mm、8.4mm和6.2mm。此外，Ⅴ級(jí)圍巖拱頂和拱肩的沉降變形差異、以及2條測(cè)試的收斂變形差異均更大，主要原因是由于隧洞Ⅳ級(jí)圍巖段巖石的完整性和穩(wěn)定性均更好，在開(kāi)挖之后，隧洞整體上變形更小，故體現(xiàn)出不同位置的沉降和收斂變形差異小。

圖2 隧洞Ⅴ級(jí)圍巖沉降和收斂值隨時(shí)間的變化曲線

3 隧洞不同圍巖級(jí)別段變形預(yù)測(cè)

基于隧洞變形現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)，對(duì)隧洞沉降及收斂變形規(guī)律進(jìn)行理論分析，找出適用該工程地質(zhì)條件下隧洞沉降和收斂變形規(guī)律的描述理論，進(jìn)而對(duì)隧洞該斷面或臨近斷面的沉降值和水平收斂值進(jìn)行預(yù)測(cè)，從而能準(zhǔn)確掌握隧洞圍巖的變形趨勢(shì)和穩(wěn)定狀態(tài)，為后續(xù)支護(hù)時(shí)間和支護(hù)措施的選取提供理論保障。

圖3 隧洞Ⅳ級(jí)圍巖沉降和收斂值隨時(shí)間的變化曲線

3.1 基于回歸分析理論的隧洞變形預(yù)測(cè)

回歸分析方法是現(xiàn)有隧道沉降變形分析及預(yù)測(cè)的常用方法，通過(guò)對(duì)大量的變形監(jiān)測(cè)數(shù)據(jù)進(jìn)行回歸分析，基于統(tǒng)計(jì)學(xué)理論建立自變量(多為時(shí)間)與因變量(如沉降值等)的擬合關(guān)系。根據(jù)建立的回歸擬合函數(shù)，可對(duì)不同時(shí)刻的隧道變形值進(jìn)行計(jì)算，預(yù)測(cè)隧道的最終沉降值。根據(jù)現(xiàn)有研究表明，隧道變形值與監(jiān)測(cè)時(shí)間的擬合關(guān)系一般為非線性關(guān)系。因此，本文借鑒隧道變形預(yù)測(cè)中所采用的回歸分析理論，對(duì)該引水隧洞的沉降變形和收斂變形進(jìn)行分析。根據(jù)對(duì)隧洞變形數(shù)據(jù)進(jìn)行初步分析，變形值與時(shí)間為非線性關(guān)系，故分別采用對(duì)數(shù)函數(shù)和指數(shù)函數(shù)對(duì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行擬合分析，發(fā)現(xiàn)該隧洞的沉降變形和收斂變形均更符合指數(shù)函數(shù)關(guān)系，其具體關(guān)系式如下：

S=A·(1-e-B·t)

(1)

式中，S—隧洞累積沉降或累積水平收斂變形值；t—監(jiān)測(cè)時(shí)間；A和B—擬合回歸系數(shù)，與隧洞圍巖性質(zhì)相關(guān)。

第2節(jié)分析發(fā)現(xiàn)，拱肩測(cè)點(diǎn)A和測(cè)點(diǎn)C的沉降值較為接近，且測(cè)線1-2和測(cè)線3-4的水平收斂值差異也較小，因此本節(jié)只對(duì)測(cè)點(diǎn)A和測(cè)線1-2的變形數(shù)據(jù)進(jìn)行了分析。隧洞Ⅴ級(jí)圍巖和Ⅳ級(jí)圍巖變形數(shù)據(jù)的回歸擬合結(jié)果見(jiàn)表1。對(duì)于隧洞不同圍巖級(jí)別，其拱頂和拱肩的沉降量、水平收斂量的擬合相關(guān)系數(shù)均大于0.9，表明擬合程度好，可采用指數(shù)函數(shù)對(duì)該隧洞的沉降和收斂變形進(jìn)行預(yù)測(cè)分析。此外，根據(jù)擬合結(jié)果可以看出，擬合系數(shù)A和B分別隨隧洞圍巖級(jí)別的增高而增大和降低，對(duì)拱頂沉降變形和水平收斂變形而言，與Ⅳ級(jí)圍巖相比，Ⅴ級(jí)圍巖的擬合系數(shù)A提高了131%和128%，系數(shù)B降低了34%和35%。結(jié)合監(jiān)測(cè)數(shù)據(jù)變化趨勢(shì)可知，系數(shù)A代表著隧洞的最終變形值，A值越大，隧洞最終變形也越大；系數(shù)B代表曲線趨于穩(wěn)定時(shí)的變化速率，B值越大，隧洞變形值越早趨于穩(wěn)定。

表1 隧洞不同級(jí)別圍巖段變形值的回歸擬合參數(shù)

3.2 基于灰色理論的隧洞變形預(yù)測(cè)

隨著灰色系統(tǒng)理論的不斷發(fā)展，學(xué)者們將其引入到隧道、路基等沉降變形分析及預(yù)測(cè)中[13]?；疑碚摰闹饕枷胧且浴安糠中畔⒁阎糠中畔⑽粗钡男颖?、“貧信息”的不確定系統(tǒng)為研究對(duì)象，通過(guò)對(duì)已有信息進(jìn)行分析，從而把握系統(tǒng)的規(guī)律?；诂F(xiàn)有學(xué)者的研究結(jié)果發(fā)現(xiàn)，灰色理論中GM(1，1)模型的適用性最好[12]，因此本文采用GM(1，1)模型對(duì)該隧洞變形數(shù)據(jù)進(jìn)行分析。

隧洞在第i天的變形值及其累積變形值分別為x(0)(i)和x(1)(i)，因此可以得出隧道在不同監(jiān)測(cè)時(shí)間的變形值序列：

X(0)={x(0)(1)，x(0)(2)，x(0)(3)，…，x(0)(n)}

(2)

X(1)={x(1)(1)，x(1)(2)，x(1)(3)，…，x(1)(n)}

(3)

(4)

式中，a和b—模型參數(shù)，通過(guò)根據(jù)最小二乘法方法獲得，具體計(jì)算如下：

A=[ab]T=(BTB)-1BTYn

(5)

式中，A—參數(shù)序列。矩陣B和Yn分別定義為：

(6)

通過(guò)求解微分方程(4)，結(jié)合初始條件x(1)(1)=x(0)(1)，可獲得GM(1，1)模型的解x(1)(t)：

(7)

基于上式，可獲得隧洞變形數(shù)據(jù)的響應(yīng)序列x(1)(k+1)：

(8)

對(duì)于上述預(yù)測(cè)模型，采用平均相對(duì)誤差Δ進(jìn)行評(píng)價(jià)其預(yù)測(cè)精度，其計(jì)算公式為：

(9)

根據(jù)隧洞變形監(jiān)測(cè)數(shù)據(jù)，帶入式(6)可計(jì)算出矩陣B和Yn，在結(jié)合矩陣運(yùn)算公式(5)可獲得GM(1，1)模型參數(shù)a和b，最后帶入式(8)可獲得不同時(shí)刻隧洞的變形預(yù)測(cè)值。根據(jù)上述計(jì)算原理，帶入該隧洞拱頂B點(diǎn)和拱肩A點(diǎn)的沉降監(jiān)測(cè)數(shù)據(jù)，以及測(cè)線1-2的收斂數(shù)據(jù)，可計(jì)算出該隧洞Ⅴ級(jí)圍巖和Ⅳ級(jí)圍巖的沉降和收斂變形值的GM(1，1)預(yù)測(cè)模型，模型參數(shù)見(jiàn)表2。模型預(yù)測(cè)的相對(duì)誤差大部分均小于0.05，預(yù)測(cè)精度為二級(jí)，說(shuō)明預(yù)測(cè)效果較好。此外，從GM(1，1)預(yù)測(cè)模型參數(shù)可以看出，系數(shù)a和b分別隨隧洞圍巖級(jí)別的增高而降低和增大，致使b/a值隨圍巖級(jí)別的增高而增大；對(duì)拱頂沉降變形和水平收斂變形而言，Ⅴ級(jí)圍巖的b/a值相比Ⅳ級(jí)圍巖分別提高了133%和124%。根據(jù)隧洞變形監(jiān)測(cè)數(shù)據(jù)的響應(yīng)序列(式(8))，b/a值代表著隧洞的最終變形值，b/a值越大，隧洞最終變形也越大。

表2 隧洞不同級(jí)別圍巖段變形值的灰色理論有關(guān)參數(shù)

3.3 隧洞變形預(yù)測(cè)模型評(píng)價(jià)

根據(jù)3.1節(jié)和3.2節(jié)中對(duì)隧洞變形監(jiān)測(cè)數(shù)據(jù)進(jìn)行回歸擬合和GM(1，1)灰色模型分析，獲得了隧洞不同級(jí)別圍巖段的沉降和收斂變形值的預(yù)測(cè)模型，將指數(shù)函數(shù)擬合曲線和GM(1，1)預(yù)測(cè)模型曲線繪制在同一圖中，并和實(shí)際變形監(jiān)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析。圖4—7給出了隧洞Ⅴ級(jí)圍巖和Ⅳ級(jí)圍巖拱頂/拱肩的沉降變形值和收斂變形值的預(yù)測(cè)曲線?？傮w上而言，指數(shù)型擬合曲線和GM(1，1)預(yù)測(cè)曲線較為接近，且靠近隧洞實(shí)際變形監(jiān)測(cè)值，表明2種方法均能較好的對(duì)該隧洞的沉降和收斂變形進(jìn)行預(yù)測(cè)分析。然而，2種分析預(yù)測(cè)方法間也存在一定的差異，相比指數(shù)型回歸預(yù)測(cè)曲線，GM(1，1)預(yù)測(cè)曲線更為靠近實(shí)際監(jiān)測(cè)點(diǎn)，尤其是在隧洞沉降變形預(yù)測(cè)分析中(如圖4和圖6)，表明在現(xiàn)有監(jiān)測(cè)數(shù)據(jù)體量下，采用灰色理論對(duì)隧洞沉降變形分析及預(yù)測(cè)的精度更高，更能反映該隧洞的沉降變形；而在隧洞水平收斂預(yù)測(cè)方面，2種方法預(yù)測(cè)精度差別不大。此外，從2種方法的理論曲線的變化趨勢(shì)可以看出，Ⅳ級(jí)圍巖的沉降變形和收斂變形相比Ⅴ級(jí)圍巖更早趨于穩(wěn)定，Ⅴ級(jí)圍巖發(fā)生較大沉降變形和收斂變形的時(shí)間持續(xù)更長(zhǎng)。

圖4 隧洞Ⅴ級(jí)圍巖沉降變形的預(yù)測(cè)曲線

圖5 隧洞Ⅴ級(jí)圍巖收斂變形的預(yù)測(cè)曲線

圖6 隧洞Ⅳ級(jí)圍巖沉降變形的預(yù)測(cè)曲線

圖7 隧洞Ⅳ級(jí)圍巖收斂變形的預(yù)測(cè)曲線

基于3.1節(jié)和3.2節(jié)中所建立的隧洞變形數(shù)據(jù)的指數(shù)回歸模型和GM(1，1)灰色模型，可以對(duì)后續(xù)不同時(shí)刻的沉降及收斂變形值進(jìn)行預(yù)測(cè)。基于指數(shù)函數(shù)模型和GM(1，1)模型的隧洞最終變形預(yù)測(cè)值見(jiàn)表3，首先利用2種理論對(duì)30d的累積沉降值和水平收斂值進(jìn)行預(yù)測(cè)，并與30d的實(shí)際變形監(jiān)測(cè)值進(jìn)行對(duì)比分析，發(fā)現(xiàn)2種理論的相對(duì)誤差均小于8%，且絕大部分的預(yù)測(cè)值誤差小于3%，再次表明2種理論對(duì)該隧洞變形值預(yù)測(cè)的可行性和準(zhǔn)確性。另外，觀察到GM(1，1)灰色模型的預(yù)測(cè)誤差整體上相比指數(shù)回歸模型更低。利用2種理論對(duì)隧洞的最終變形值進(jìn)行預(yù)測(cè)，對(duì)拱頂B點(diǎn)沉降而言，指數(shù)回歸模型和GM(1，1)灰色模型的預(yù)測(cè)值分別為7.71mm和7.89mm(Ⅴ級(jí)圍巖)、17.84mm和18.44mm(Ⅳ級(jí)圍巖)；測(cè)線1-2最終水平收斂值分別為4.36mm和4.43mm(Ⅴ級(jí)圍巖)、9.97mm和9.91mm(Ⅳ級(jí)圍巖)。總體上，2種理論預(yù)測(cè)值差別較小，但GM(1，1)灰色模型的預(yù)測(cè)值稍大，更加保守。綜上所述，對(duì)于該隧洞而言，采用GM(1，1)灰色模型對(duì)隧洞沉降及收斂變形預(yù)測(cè)更為精確。

表3 基于指數(shù)函數(shù)模型和GM(1，1)模型的隧洞最終變形預(yù)測(cè)值

4 結(jié)論

本文對(duì)比分析了某隧洞不同級(jí)別圍巖段拱頂/拱肩沉降和水平收斂的變化規(guī)律，并引入回歸分析理論和灰色系統(tǒng)理論對(duì)變形數(shù)據(jù)進(jìn)行擬合分析和評(píng)價(jià)，預(yù)測(cè)了隧洞最終變形值。主要結(jié)論如下：

(1)隧洞Ⅴ級(jí)和Ⅳ級(jí)圍巖的拱頂-拱肩沉降值均隨時(shí)間呈先快速增加后緩慢增長(zhǎng)趨勢(shì)，且拱頂沉降大于拱肩、V級(jí)圍巖沉降大于Ⅳ級(jí)圍巖。

(2)指數(shù)型函數(shù)和GM(1，1)模型均能較好地描述該隧洞沉降及收斂變形，但較小體量的監(jiān)測(cè)數(shù)據(jù)使GM(1，1)模型預(yù)測(cè)精度更高，Ⅳ級(jí)和Ⅴ級(jí)圍巖的最終沉降值和收斂值分別為7.89mm和4.43mm、18.44mm和9.91mm。

研究成果對(duì)類似隧洞沉降收斂變形分析、預(yù)測(cè)和控制等具有重要意義，但沒(méi)有綜合考慮開(kāi)挖等施工因素的影響。