■阮麟

（福建省交通建設(shè)工程試驗檢測有限公司，福州350008）

新聯(lián)隧道塌方前后的監(jiān)控量測和數(shù)據(jù)分析

■阮麟

（福建省交通建設(shè)工程試驗檢測有限公司，福州350008）

監(jiān)控量測是及時掌握隧道施工中圍巖動態(tài)變化過程的重要手段，通過對量測數(shù)據(jù)的回歸分析計算可以預(yù)測可能發(fā)生危及施工安全的隱患或事故，進(jìn)而有效地指導(dǎo)隧道的設(shè)計與施工。本文結(jié)合新聯(lián)隧道施工的監(jiān)控量測，根據(jù)4種回歸函數(shù)和相關(guān)性系數(shù)對量測數(shù)據(jù)進(jìn)行對比分析，發(fā)現(xiàn)多項式回歸函數(shù)適合對隧道拱頂沉降的預(yù)測，并且成功預(yù)測了隧道的塌方。塌方經(jīng)過治理后，指數(shù)回歸函數(shù)對拱頂沉降的預(yù)測更接近實際。

隧道拱頂沉降監(jiān)控量測回歸分析拱頂沉降

1 前言

我國作為一個多山嶺的國家，山地和丘陵占到了國土總面積的三分之二。因此，今后我國的高等級公路建設(shè)將主要在山丘地區(qū)，進(jìn)而由于地形及選線的限制，必然伴隨著大量山區(qū)公路隧道的建設(shè)，隧道開挖后洞身會不可避免地產(chǎn)生變形[1]。對變形的監(jiān)控量測是隧道施工管理中不可缺少的一個重要環(huán)節(jié)，現(xiàn)場監(jiān)控量測，可以及時掌握在開挖過程中圍巖穩(wěn)定和支護(hù)受力、變形的動態(tài)信息，及時、準(zhǔn)確地預(yù)報可能發(fā)生的危及施工安全的隱患或事故,并及時反饋動態(tài)監(jiān)測信息給相關(guān)部門，使相關(guān)部門及時修正設(shè)計支護(hù)參數(shù)和調(diào)整施工措施,從而達(dá)到設(shè)計與施工安全、經(jīng)濟(jì)的目的。

本文以新聯(lián)隧道工程的監(jiān)控量測為例，根據(jù)4種回歸函數(shù)和相關(guān)性系數(shù)對量測數(shù)據(jù)進(jìn)行對比分析，得到多項式回歸函數(shù)更適合隧道拱頂沉降的預(yù)測，并且成功地預(yù)測了隧道的塌方；但塌方治理后，指數(shù)回歸函數(shù)卻更適合對拱頂沉降的預(yù)測，即對塌方前后進(jìn)行分段回歸分析計算[2]。

2 隧道概況

新聯(lián)隧道屬海西高速公路網(wǎng)廈沙線上的工程，位于三明市尤溪縣內(nèi)，走向近南北，屬分離式隧道，最大埋深165m。左線隧道起迄樁號ZK144+120～ZK145+406，全長1286m；右線隧道起迄樁號YK144+102～YK145+386,全長1295m，海拔高度為205～225m。

2.1 工程地質(zhì)條件

（1）地形地貌

隧址區(qū)屬剝蝕丘陵地貌，地形呈波狀起伏。地勢較緩，支溝發(fā)育，多呈“V”字形，山坡坡度10～35°，進(jìn)洞口自然山坡坡度約為10～15°，自然山坡穩(wěn)定；出洞口自然山坡坡度約20～35°。區(qū)內(nèi)最高點海拔約380m，溝塹較發(fā)育，見地表水，溝內(nèi)常年有水流動，地表水較豐富，水量較大。沿線植被茂盛，多為山林地。

（2）地層巖性

在隧道進(jìn)口段，巖性主要從上覆地層和下伏地層角度來描述，具體巖性分述如下：

覆蓋層。山坡為殘積砂質(zhì)粘性土夾碎石,沖溝地段為碎石、礫砂等洪積物。

下伏地層燕山早期侵入及花崗巖其風(fēng)化層，花崗巖為淺棕黃、灰黑灰褐，中—粗粒似斑狀花崗結(jié)構(gòu)，塊狀構(gòu)造，無斑晶，主要礦物為鉀長石、斜長石、正長石及石英、云母、角閃石。但節(jié)理裂隙較發(fā)育，差異風(fēng)化嚴(yán)重。風(fēng)化層的巖體破碎，呈散體狀，巖芯呈砂土狀，手易捏碎，巖體基本質(zhì)量等級屬Ⅴ級。

（3）水文地質(zhì)條件

隧道位于當(dāng)?shù)厍治g基準(zhǔn)面以下，山坡坡體起伏，有地表水經(jīng)過。隧址區(qū)地表水主要為分布于進(jìn)出口段的山間小溪，常年有水，主要接受大氣降水的補(bǔ)給，受季節(jié)性影響較大。地下水主要為風(fēng)化帶網(wǎng)狀空隙裂隙水、基巖裂隙水及構(gòu)造裂隙水。進(jìn)洞口處賦水條件非常好，勘查期間地下水高于洞頂，進(jìn)洞后地下水呈淋雨狀出水。

2.2 施工概況

隧道采用新奧法設(shè)計和施工，鉆爆法開挖。淺埋、地質(zhì)條件差的洞口Ⅴ級圍巖段采用單側(cè)壁導(dǎo)坑法（CD法）開挖[2-3]。開挖步序是左右兩個側(cè)上導(dǎo)洞交錯開挖，兩導(dǎo)洞的掌子面保持開挖距離15m以上；兩個側(cè)下導(dǎo)洞再跟進(jìn)，與兩上導(dǎo)洞保持距離15～20m；豎向的臨時鋼支撐直至仰拱封閉后拆除。

3 監(jiān)控量測的方法

3.1 斷面測點布設(shè)及測量方法

拱頂沉降監(jiān)測點每個斷面共布設(shè)3個，即在拱頂布設(shè)測點2#，在其左右兩側(cè)距離該點2m的位置布設(shè)測點1#和測點3#；布設(shè)一條水平監(jiān)測點（對測點），各斷面拱頂沉降監(jiān)測點布置示意圖如圖1所示。監(jiān)測點做法是在各測點上先安裝膨脹螺栓，再在膨脹螺栓上焊上鐵片并粘貼全站儀專用反光片，做好標(biāo)記[4]。由于隧道圍巖較差，監(jiān)測斷面布設(shè)比較密集，因此為便于量測，拱頂沉降采用全站儀而非水準(zhǔn)儀。采用全站儀量測時，先測出測點1、2、3的高程H1、H2、H3；第二天，同樣測出H1’、H2’、H3’；用△表示個監(jiān)測點的拱頂測量值，△=H1-H1’，用三點的△值求出平均值，則可認(rèn)為是該斷面的相鄰兩天的相對沉降值，以后類推則可得出每天的相對沉降值，累積后則得到該斷面的拱頂?shù)睦鄯e沉降值，其直接表示了該斷面拱頂沉降面型的大小及變化速率，同理可得到圍巖收斂的累積收斂值。

圖1 隧道監(jiān)控量測布點示意圖

3.2 測量頻率

測量頻率的選擇在隧道位移監(jiān)測中非常重要，它直觀記錄了隧道開挖施工過程中圍巖的變形情況，頻率太密集則測量強(qiáng)度大，太疏松卻又不能及時反映位移變化情況，因此參考《公路隧道施工技術(shù)規(guī)范》的規(guī)定，監(jiān)測頻率選取如下：

工作面開挖15天內(nèi)，1～2次/天；16天～1個月內(nèi)，1次/天；1～3個月內(nèi)，1～2次/周；大于3個月，1～3次/月。

4 量測數(shù)據(jù)分析

新聯(lián)隧道里程較長，監(jiān)測斷面布置較多，本文僅以右線里程樁號為YK144+159的典型斷面的拱頂沉降測量數(shù)據(jù)為例，應(yīng)用回歸方法來分析預(yù)測隧道施工過程中拱頂沉降的變形量情況，其他斷面的數(shù)據(jù)可以依此類推。

4.1 數(shù)據(jù)處理

斷面YK144+159施工完成初期支護(hù)后即開始對圍巖水平收斂及拱頂沉降，進(jìn)行監(jiān)控量測。量測所得的數(shù)據(jù)，及時整理成所需要的原始數(shù)據(jù)，并繪制累計位移—時間折線散點圖。由圖2可知拱頂最大沉降速率和圍巖最大收斂速率分別為-4.8mm/d和-5.9mm/d。

圖2 YK144+159時間-位移曲線圖

4.2 數(shù)據(jù)回歸分析

由于現(xiàn)場量測所取得的原始數(shù)據(jù)，不可避免會具有一定的離散性，其中包括測量誤差甚至測試誤差，因而不經(jīng)過數(shù)學(xué)處理的量測數(shù)據(jù)是難以直接利用的。回歸分析是目前量測數(shù)據(jù)數(shù)學(xué)處理的主要方法，主要是利用數(shù)理統(tǒng)計的原理，建立不同變量之間相關(guān)關(guān)系的數(shù)學(xué)表達(dá)式。在隧道開挖的變形分析中，可用回歸分析方法建立因變量與自變量（即變形與影響因子）之間的數(shù)學(xué)關(guān)系式，從而獲得能較準(zhǔn)確反映實際情況的u-t曲線圖，然后可以預(yù)測圍巖的最終位移值和各階段的位移速率[5-9]。

目前主要是對指數(shù)、對數(shù)和雙曲三種曲線函數(shù)進(jìn)行線性回歸和多項式函數(shù)進(jìn)行計算，三種曲線函數(shù)的原形公式與換算公式如下：

公式中A、B為回歸參數(shù)；u為變形值；t為觀測時間。

由圖2可知，隧道圍巖的凈空收斂和拱頂沉降基本呈相似的變化規(guī)律，即：加速變形—突然緩慢變形—急劇變形，故取拱頂沉降進(jìn)行回歸分析。但經(jīng)過回歸計算，指數(shù)函數(shù)和對數(shù)函數(shù)的A、B相關(guān)性系數(shù)分別為R1=0.862和R2=0.700，都小于0.9，所以不可采用，而利用雙曲線函數(shù)進(jìn)行回歸分析時相關(guān)性R=0.997，其線性回歸計算見下表1。

表1 線性回歸分析計算

但是從表中可以看出最后幾天的預(yù)測值與觀測值值相差太大，故也不可采用雙曲線函數(shù)進(jìn)行回歸分析計算，所以采用多項式函數(shù)計算。

多項式回歸函數(shù)是根據(jù)監(jiān)測數(shù)據(jù)離散點結(jié)合計算機(jī)軟件自動生成的回歸趨勢線，回歸曲線與實測曲線如下圖3，所示多項式函數(shù)表達(dá)式為：

階數(shù)n=4時，R=0.998

圖3 回歸曲線與實測曲線圖

由上述回歸曲線函數(shù)及相關(guān)系數(shù)可知，4個回歸模型中，相關(guān)系數(shù)R的排列順序為:對數(shù)＜指數(shù)＜雙曲線＜多項式＜1，而且多項式回歸曲線基本上與觀測點重合，可見多項式函數(shù)的回歸精度很高，用來預(yù)測該斷面拱頂最終沉降量最精確。

4.3 塌方預(yù)測

當(dāng)t=23時，代入函數(shù)得到預(yù)測累積沉降值u= 50.2mm，將會超過預(yù)警值，且隨著左右導(dǎo)洞的開挖掘進(jìn)，在掌子面空間效應(yīng)和應(yīng)力釋放的雙重作用下，松軟的殘坡積土，變得強(qiáng)度越來越低、越來越軟，可能發(fā)生塌方。根據(jù)現(xiàn)場的反饋，在施工的第22天中午，新聯(lián)隧道右洞的左導(dǎo)洞YK144+159～YK144+162發(fā)生突發(fā)性滑塌（見圖4）。所以運(yùn)用回歸分析計算，成功預(yù)測了隧道的塌方。

圖4 塌方形態(tài)圖

4.4 塌方治理后

塌方段YK144+159～YK144+162塌方處治后即開始對圍巖水平收斂及拱頂沉降進(jìn)行監(jiān)控量測。但是多項式計算的回歸位移與最終位移值都大于實測位移值，兩組變量的相關(guān)性系數(shù)R也不接近1。因此可進(jìn)行分段回歸分析，即塌方之前多項式函數(shù)回歸計算，塌方后運(yùn)用其他函數(shù)回歸計算。根據(jù)收斂·—時間曲線走向圖（見圖5），宜選用指數(shù)函數(shù)公式作為回歸函數(shù)，計算結(jié)果見表2。

圖5 塌方治理后時間—位移曲線

表2 指數(shù)函數(shù)回歸計算

計算最終得到回歸函數(shù)。

當(dāng)R=0.973時，可得到拱頂下沉位移的預(yù)測值的終值u=4.13mm。

綜合以上分析，可知經(jīng)過治理后，拱頂最大沉降速率為-0.73mm/d，遠(yuǎn)小于塌方前的-4.8mm/d，而且累積沉降位移趨向穩(wěn)定，最終位移預(yù)測值僅為4.13mm。

5 結(jié)論

（1）通過對拱頂沉降回歸分析計算，成功地預(yù)測了隧道的塌方。

（2）塌方經(jīng)過處治理后，最終的沉降累積位移和速率都遠(yuǎn)小于塌方前的，并且達(dá)到在一個安全的穩(wěn)定值。

（3）為了更好地指導(dǎo)施工和準(zhǔn)確預(yù)測極值，根據(jù)塌方前后得到的隧道監(jiān)控量測數(shù)據(jù)應(yīng)進(jìn)行分段回歸分析。

[1]呂蓬.內(nèi)鄧高速公路中條山隧道安全施工技術(shù)研究[D].西安:長安大學(xué).2012.

[2]張國華，陳禮彪，錢師雄，等．大斷面小凈距大帽山隧道監(jiān)控量測及分析[J].巖土力學(xué),32（2）:497-502.2011.

[3]張國華，陳禮彪，夏祥,等．大斷面隧道爆破開挖圍巖損傷試驗研究及數(shù)值計算.巖石力學(xué)與`工程學(xué)報，2011,30（9）:1610-1619.

[4]周亦濤,等.雙線隧道現(xiàn)場監(jiān)控量測數(shù)據(jù)分析[J].研究與設(shè)計，2007（3）:26-28.

[5]趙玉光.地下工程開挖與支護(hù)安全監(jiān)控量測及地質(zhì)災(zāi)害動態(tài)預(yù)測與防治[J].中國地質(zhì)災(zāi)害與防治學(xué)報，2001,12（3）:36-40.

[6]王勝濤，梁小勇，周亦濤.隧道監(jiān)控量測的數(shù)據(jù)回歸分析探討[J].隧道建設(shè)，2009,29（6）:629-632.

[7]張俊中，宋蕾，張健雄.多元回歸分析模型在變形監(jiān)測中的應(yīng)用[J].河南工程學(xué)院院報（自然科學(xué)版）,2009,21（3）:22-25..

[8]丁健，王超林，余熠，王興安.某隧道塌方前后的監(jiān)控量測及數(shù)據(jù)分析[J].采礦技術(shù)，2014,14（6）:32-36.

[9]陳秋南，張永興，劉新榮，等.隧道塌方段加固后的施工監(jiān)測與仿真分析[J].巖石力學(xué)與工程學(xué)報，2006,25（1）:158-161.