劉學(xué)軍，高玉峰，賀一凡，姜兆東

(1.新疆建筑科學(xué)研究院(有限責(zé)任公司)，新疆烏魯木齊830002；2.新疆大學(xué)建筑工程學(xué)院，新疆烏魯木齊830046)

0 引 言

隧道工程是交通運(yùn)輸工程的重要組成部分，受地質(zhì)災(zāi)害的影響，在施工過(guò)程中經(jīng)常發(fā)生安全問(wèn)題，甚至?xí){到工程的安全穩(wěn)定[1]。對(duì)隧道圍巖的穩(wěn)定性提前預(yù)測(cè)可以減少工程事故的發(fā)生[2]。針對(duì)圍巖穩(wěn)定性的問(wèn)題，國(guó)內(nèi)外學(xué)者已開(kāi)展了大量研究工作。許傳華等[3]采用模糊數(shù)學(xué)綜合評(píng)估法和模糊語(yǔ)言等計(jì)算方法，對(duì)圍巖穩(wěn)定性進(jìn)行了不同程度的評(píng)估。曾懷宇[4]為分析地質(zhì)水文環(huán)境及巖石自身?xiàng)l件等眾多因素對(duì)圍巖分類(lèi)的影響，在已收集數(shù)據(jù)的基礎(chǔ)上，采用粗糙集理論對(duì)數(shù)據(jù)中的隱含知識(shí)進(jìn)行挖掘，實(shí)現(xiàn)對(duì)隧道圍巖等級(jí)的預(yù)判。儲(chǔ)漢東等[5]進(jìn)行數(shù)值模擬，研究圍巖的應(yīng)力和應(yīng)變變化場(chǎng)。李健等[6]基于云模型理論探討了圍巖穩(wěn)定性分類(lèi)模型。閆天俊等[7]對(duì)弱結(jié)構(gòu)面流動(dòng)修正系數(shù)進(jìn)行修正，并提出了考慮到巖石流動(dòng)影響的圍巖分類(lèi)指標(biāo)。郝杰等[8]通過(guò)使用基于拓?fù)湓乩碚摰暮?jiǎn)單相關(guān)函數(shù)，確定影響圍巖的每個(gè)因素的權(quán)重。王佳信等[9]提出一種基于Alpha穩(wěn)定分析(SaS)-概率神經(jīng)網(wǎng)絡(luò)(PNN)模型，對(duì)圍巖穩(wěn)定性進(jìn)行了預(yù)測(cè)。貢力等[10]為對(duì)圍巖的安全穩(wěn)定性做出準(zhǔn)確預(yù)測(cè)，提出了基于PCA-IRBF的計(jì)算模型，對(duì)隧道圍巖的安全穩(wěn)定性進(jìn)行預(yù)測(cè)評(píng)價(jià)的新方法。李宏波[11]為探究TBM掘進(jìn)量圍巖等級(jí)關(guān)系，采用數(shù)據(jù)清洗、分布統(tǒng)計(jì)等多種方法，建立了以神經(jīng)網(wǎng)絡(luò)聚類(lèi)和最小二乘支持向量機(jī)相結(jié)合的圍巖等級(jí)預(yù)測(cè)方法。周浩宇等[12]以江西某隧道為例，根據(jù)已有工程經(jīng)驗(yàn)，選取合適的影響因素建立指標(biāo)體系，采用理想點(diǎn)法建立隧道坍塌風(fēng)險(xiǎn)等級(jí)預(yù)測(cè)模型。李蓬喜[13]采用支持向量回歸、卷積神經(jīng)網(wǎng)絡(luò)工智能算法，對(duì)TBM的掘進(jìn)數(shù)據(jù)、圍巖等級(jí)進(jìn)行預(yù)測(cè)分析。

在東天山隧道施工研究中，秦?fù)碥姷萚2]利用TSP物探技術(shù)與灰色關(guān)聯(lián)度等技術(shù)對(duì)東天山隧道施工前期的數(shù)據(jù)進(jìn)行了分析，并對(duì)圍巖質(zhì)量進(jìn)行預(yù)測(cè)研究，但其研究是基于東天山隧道施工前期的數(shù)據(jù)，數(shù)據(jù)量較少導(dǎo)致其研究成果仍存在不足。為此，本文在秦?fù)碥姷热搜芯康幕A(chǔ)上，采用東天山隧道施工后期的數(shù)據(jù)進(jìn)行研究分析，所得數(shù)據(jù)更加完備，監(jiān)測(cè)數(shù)據(jù)更為豐富，可使研究結(jié)果更加準(zhǔn)確可靠。此外，在實(shí)際工程中獲取的數(shù)據(jù)往往是龐大的，需要找到一個(gè)可以大規(guī)模處理工程數(shù)據(jù)的方法。隨著人工智能的不斷發(fā)展，具有大規(guī)模并行處理、高運(yùn)算速度、強(qiáng)適應(yīng)性以及強(qiáng)容錯(cuò)能力等高性能的神經(jīng)網(wǎng)絡(luò)被廣泛應(yīng)用。鑒于以上原因，本文基于優(yōu)化后的東天山隧道施工后期數(shù)據(jù)，利用高性能的神經(jīng)網(wǎng)絡(luò)等處理技術(shù)，對(duì)圍巖隧道質(zhì)量進(jìn)行預(yù)測(cè)。

1 工程概況

新疆東天山隧道全長(zhǎng)約11 km，位于高山區(qū)和中山區(qū)，融雪和降雨等因素易地表形成徑流。受多時(shí)期構(gòu)造運(yùn)動(dòng)的影響，導(dǎo)致項(xiàng)目沿線出露的基巖節(jié)理較為發(fā)育[14]。隧道地理位置見(jiàn)圖1。

圖1 隧道地理位置

根據(jù)現(xiàn)場(chǎng)踏勘與鉆探取樣分析，隧址區(qū)上部主要為粉質(zhì)黏土、卵石和塊石等，屬于第四系全新統(tǒng)～上更新統(tǒng)；隧址區(qū)下部基巖巖性主要為凝灰質(zhì)砂巖、凝灰?guī)r、熔結(jié)泥灰?guī)r和花崗巖、輝綠巖、花崗斑巖等，屬于石炭系中統(tǒng)及華力西期侵入巖。東天山隧道主要跨越了2個(gè)三級(jí)構(gòu)造單元，走向?yàn)楸蔽鳌蔽魑鞯膮^(qū)域構(gòu)造占了絕大部分，走向?yàn)楸睎|或近南北的僅有極少部分。北天山北坡陡、南坡緩，北部沖洪積扇較南部山麓發(fā)育，南部新構(gòu)造活動(dòng)明顯弱于北部。

2 基本理論分析

2.1 主成分分析關(guān)聯(lián)理論

主成分分析是用1組新的相互不相關(guān)的綜合指標(biāo)替代大量的指標(biāo)(如P指數(shù))。當(dāng)統(tǒng)計(jì)分析被用來(lái)研究1個(gè)多變量的課題時(shí)，太多的變量會(huì)增加課題的復(fù)雜性，希望以較少的變量獲得更多的信息。在存在多因素影響的條件下，多變量之間存在內(nèi)在的聯(lián)系。如果2個(gè)變量之間有內(nèi)在關(guān)聯(lián)，則認(rèn)為這2個(gè)變量反映的主體信息有一定的重疊。主成分分析是將重復(fù)的變量(密切相關(guān)的變量)從最初提出的所有變量中刪除，并創(chuàng)建盡可能少的新變量，使這些新變量是2個(gè)或更多的不相關(guān)的變量，這些新變量在反映被試的信息能盡可能多地保留原有信息。

最經(jīng)典的方法是用F1的方差來(lái)表示。若全部組合中的F1方差最高，則F1稱(chēng)為第1主成分。假如F1難以反映原始指標(biāo)P的信息，再考慮F2，使其有效地代表原始信息。同時(shí)，F(xiàn)2中的信息不再包括F1中的信息。也可以表達(dá)為要求Cov(F1，F(xiàn)2)=0，則將F2稱(chēng)為第2主成分。以此類(lèi)推可以構(gòu)造第i個(gè)直至第P個(gè)主成分[15]。

2.2 TSP物探指標(biāo)敏感性分析

考慮巖石硬度和巖石完整性2個(gè)基本因素的質(zhì)量特征，以及基本的圍巖質(zhì)量指標(biāo)BQ進(jìn)行初步分類(lèi)[16]?；举|(zhì)量指數(shù)的修正值是考慮到修正因素的影響而得到的。巖體的基本圍巖質(zhì)量指標(biāo)BQ的計(jì)算公式為

BQ=100+3Rc+250Kr

(1)

式中，RC為巖石單軸抗壓強(qiáng)度；Kr為巖石完整度系數(shù)。若RC>90Kr+30，則將RC=90Kr+30和Kr代入計(jì)算BQ值；當(dāng)Kr>0.04Rc+0.4時(shí)，應(yīng)以Kr=0.04Rc+0.4和Rc代入計(jì)算BQ值。

計(jì)算出BQ值后，結(jié)合JTG 3370.1—2018《公路隧道設(shè)計(jì)規(guī)范》給出的地下水影響修正系數(shù)K1、主要軟弱結(jié)構(gòu)面產(chǎn)狀影響修正系數(shù)K2、初始應(yīng)力狀態(tài)影響修正系數(shù)K3進(jìn)行修正計(jì)算得出圍巖質(zhì)量指標(biāo)修正值[BQ]，計(jì)算公式為

[BQ]=BQ-100(K1+K2+K3)

(2)

就東天山隧道而言，工程條件復(fù)雜，在隧道的長(zhǎng)段上，干燥和水淹的條件交替出現(xiàn)，而且結(jié)構(gòu)的表面幾何形狀也不同。僅采用3個(gè)修正條件難以全部滿足所有條件，同時(shí)修正系數(shù)的取值將受到影響，導(dǎo)致計(jì)算出的[BQ]修正值存在誤差。對(duì)地下水修正系數(shù)K1進(jìn)行加權(quán)平均計(jì)算，可以使得結(jié)果更加準(zhǔn)確[17]，計(jì)算公式為

(3)

式中，ki為第i個(gè)因子的地下水修正系數(shù)；wi為第i個(gè)因子的隧道長(zhǎng)度；n為因子個(gè)數(shù)。

例如，假設(shè)一段巷道長(zhǎng)200 m，巖體基本圍巖質(zhì)量指標(biāo)BQ為330，其中有55 m是潮濕或點(diǎn)滴狀出水，且水壓小于0.1 MPa；有145 m是淋雨?duì)罨蛴苛鳡畛鏊?，其中?05 m水壓大于0.5 MPa，有40 m水壓大于0.1 MP且小于等于0.5 MPa。如果按規(guī)范計(jì)算，200 m中有105 m為淋雨?duì)罨蛴苛鳡畛鏊宜畨捍笥?.5 MPa，則地下水修正系數(shù)K1的取值為0.8，這樣計(jì)算出來(lái)的修正值[BQ]會(huì)比實(shí)際偏低。若采用加權(quán)平均法，則

3 研究數(shù)據(jù)處理

3.1 加權(quán)修正圍巖質(zhì)量指標(biāo)

在秦?fù)碥姷萚2]研究的基礎(chǔ)上，采用東天山隧道施工后期優(yōu)化后的數(shù)據(jù)進(jìn)行研究分析，計(jì)算出[BQ]值，并列出相應(yīng)的TSP物探指標(biāo)。由于數(shù)據(jù)過(guò)多，因此只列出表1所示的部分主要數(shù)據(jù)。

3.2 關(guān)聯(lián)度計(jì)算結(jié)果

將表1的數(shù)據(jù)通過(guò)MATLAB編程進(jìn)行關(guān)聯(lián)度計(jì)算，結(jié)果見(jiàn)表2。通過(guò)TSP物探指標(biāo)敏感性分析可知，在多因素物理力學(xué)參數(shù)的影響中，波速V、縱橫波速比Vp/Vs、泊松比μ和密度ρ與圍巖質(zhì)量的相關(guān)程度較高，所以可以認(rèn)為這4個(gè)因素是影響巖石質(zhì)量及穩(wěn)定性的主控因素。而剪切模量G、拉梅系數(shù)L、體積模量K、動(dòng)態(tài)楊氏模量ED和靜態(tài)楊氏模量ES與圍巖質(zhì)量的相關(guān)程度最低，說(shuō)明這5個(gè)因素對(duì)圍巖質(zhì)量及穩(wěn)定性影響相對(duì)較小。

表1 東天山隧道TSP物探指標(biāo)及對(duì)應(yīng)[BQ]值

表2 TSP物探各指標(biāo)的圍巖穩(wěn)定性關(guān)聯(lián)系數(shù)

4 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)圍巖質(zhì)量

4.1 BP神經(jīng)網(wǎng)絡(luò)關(guān)聯(lián)理論

BP算法是神經(jīng)網(wǎng)絡(luò)眾多算法中較為經(jīng)典的算法，同時(shí)也是應(yīng)用最廣泛的算法之一。BP算法搭建的神經(jīng)網(wǎng)絡(luò)主要由輸入層、隱含層和輸出層3層神經(jīng)元構(gòu)成，3層全部存在情況下，將其稱(chēng)為全連接神經(jīng)網(wǎng)絡(luò)，見(jiàn)圖2。

圖2 全連接神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

全連接神經(jīng)網(wǎng)路中，輸入層和輸出層一般為單層神經(jīng)元，而隱含層通常為雙層甚至多層神經(jīng)元。神經(jīng)網(wǎng)絡(luò)的輸入層主要負(fù)責(zé)提取樣本數(shù)據(jù)，并將數(shù)據(jù)傳遞給隱含層做深層處理，最終將處理好的數(shù)據(jù)傳遞給輸出層，而神經(jīng)網(wǎng)絡(luò)模型的準(zhǔn)確與否，通常由隱含層來(lái)決定。神經(jīng)網(wǎng)絡(luò)中處理數(shù)據(jù)的方法采用BP算法，通過(guò)比對(duì)神經(jīng)網(wǎng)絡(luò)的實(shí)際輸出與期望輸出，降低兩者誤差。經(jīng)過(guò)對(duì)神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)及參數(shù)的大量調(diào)試，當(dāng)誤差在允許范圍時(shí)，則認(rèn)為神經(jīng)網(wǎng)絡(luò)已經(jīng)訓(xùn)練完成。

4.2 預(yù)測(cè)數(shù)據(jù)處理

由表2計(jì)算得出的結(jié)果，選取V、Vp/Vs、μ和ρ這4主控因素作為隧道圍巖質(zhì)量預(yù)測(cè)的主要依據(jù)。預(yù)測(cè)樣本采用表1中隨機(jī)選取的6組數(shù)據(jù)(見(jiàn)表3)，初始數(shù)據(jù)為剩余的512組數(shù)據(jù)(見(jiàn)表4)，所取數(shù)據(jù)是在文獻(xiàn)[2]的基礎(chǔ)上進(jìn)行優(yōu)化，數(shù)據(jù)的來(lái)源是隧道工程施工后期，且取樣的部位更加合理，因此本文所用數(shù)據(jù)更有代表性。在初始數(shù)據(jù)與預(yù)測(cè)樣本的關(guān)聯(lián)度計(jì)算中，本文使用灰色關(guān)聯(lián)算法計(jì)算。在512個(gè)初始數(shù)據(jù)計(jì)算出的關(guān)聯(lián)度中選取最大關(guān)聯(lián)度所對(duì)應(yīng)的初始數(shù)據(jù)的圍巖質(zhì)量，作為該組預(yù)測(cè)樣本的結(jié)果。

表3 預(yù)測(cè)樣本(參考序列)

表4 初始數(shù)據(jù)(比較序列)

4.3 預(yù)測(cè)結(jié)果

采用BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)計(jì)算得出的結(jié)果見(jiàn)表5。

從表5可以看出，除預(yù)測(cè)樣本6外，所有5組預(yù)測(cè)樣本都可以在表4的512組試樣中找到對(duì)應(yīng)的實(shí)際樣本數(shù)據(jù)。此外，由于工程數(shù)據(jù)的不完整性，導(dǎo)致預(yù)測(cè)樣本2的預(yù)測(cè)結(jié)果存在誤差，但誤差在允許范圍內(nèi)。從整體數(shù)據(jù)來(lái)看，[BQ]值為131.8的數(shù)據(jù)總量較少，這是誤差產(chǎn)生的主要原因之一。隨著工程數(shù)據(jù)的不斷增加和完善，BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)精度將會(huì)越來(lái)越高。因此，利用地球物探指標(biāo)的高度相關(guān)性預(yù)測(cè)圍巖質(zhì)量是可行的。

表5 樣本的BP神經(jīng)網(wǎng)絡(luò)關(guān)聯(lián)度預(yù)測(cè)結(jié)果

5 結(jié) 語(yǔ)

本文以新疆東天山隧道工程為研究對(duì)象，利用施工后期優(yōu)化后的數(shù)據(jù)，采用TSP物探指標(biāo)敏感性分析法和BP神經(jīng)網(wǎng)絡(luò)關(guān)聯(lián)分析理論，對(duì)影響隧道圍巖質(zhì)量的因素進(jìn)行分析，并對(duì)圍巖質(zhì)量進(jìn)行預(yù)測(cè)，所得結(jié)論如下：

(1)TSP物探指標(biāo)敏感性分析結(jié)果表明，V、Vp/Vs、μ和ρ與圍巖質(zhì)量的相關(guān)性最高，關(guān)聯(lián)度分別高達(dá)0.802 7，0.863 3，0.855 7和0.819 4，均是影響圍巖質(zhì)量的重要因素。

(2)根據(jù)BP神經(jīng)網(wǎng)絡(luò)對(duì)圍巖質(zhì)量主要影響因素的預(yù)測(cè)結(jié)果可知，BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)方法是可靠準(zhǔn)確的，所得結(jié)果將會(huì)隨著技術(shù)和數(shù)據(jù)的完善而有所改善。

(3)圍巖的縱橫波速比和泊松比的預(yù)測(cè)精度對(duì)TSP物探指標(biāo)預(yù)測(cè)結(jié)果影響較大，提高2個(gè)因素的預(yù)測(cè)精度可有效增強(qiáng)整體預(yù)測(cè)效果。