姜海波，吳 鵬，張 軍

(1.石河子大學(xué)水利建筑工程學(xué)院，新疆石河子832000；2.新疆維吾爾自治區(qū)水利水電勘測設(shè)計研究院，新疆烏魯木齊830000)

高地溫引水隧洞圍巖與噴層結(jié)構(gòu)熱力學(xué)參數(shù)敏感性分析

姜海波1，吳 鵬2，張 軍2

(1.石河子大學(xué)水利建筑工程學(xué)院，新疆石河子832000；2.新疆維吾爾自治區(qū)水利水電勘測設(shè)計研究院，新疆烏魯木齊830000)

以新疆某水電站引水隧洞高地溫段為研究對象，以現(xiàn)場監(jiān)測獲得的環(huán)向應(yīng)力、溫度數(shù)據(jù)為基礎(chǔ)，對圍巖和噴層結(jié)構(gòu)的熱力學(xué)參數(shù)(導(dǎo)熱系數(shù)、比熱、對流系數(shù)和線膨脹系數(shù))與噴層環(huán)向應(yīng)力的關(guān)系進(jìn)行了分析，并采用參數(shù)敏感性分析法，研究圍巖與噴層熱力學(xué)參數(shù)對噴層應(yīng)力的敏感性。結(jié)果表明，圍巖和噴層的熱力學(xué)參數(shù)對噴層拱頂環(huán)向應(yīng)力的敏感度排序為：噴層線膨脹系數(shù)>噴層導(dǎo)熱系數(shù)>圍巖的導(dǎo)熱系數(shù)>對流系數(shù)>圍巖的線膨脹系數(shù)>圍巖的比熱>噴層的比熱，拱頂環(huán)向應(yīng)力對噴層線膨脹系數(shù)、噴層導(dǎo)熱系數(shù)、圍巖導(dǎo)熱系數(shù)、對流系數(shù)比較敏感，而對圍巖線膨脹系數(shù)、圍巖比熱和噴層比熱不敏感。

引水隧洞；圍巖熱力學(xué)參數(shù)；圍巖溫度；環(huán)向應(yīng)力；敏感性分析

0 引 言

在圍巖開挖、支護(hù)結(jié)構(gòu)的設(shè)計中，參數(shù)的敏感性分析作為系統(tǒng)研究圍巖穩(wěn)定、支護(hù)結(jié)構(gòu)力學(xué)特性的重要方法，逐步引入到巖土工程領(lǐng)域并得到了廣泛的應(yīng)用[1]。在高地溫條件下，圍巖及其支護(hù)結(jié)構(gòu)的力學(xué)特性具有很大的時空變異性，究其原因，其本質(zhì)是圍巖及其支護(hù)結(jié)構(gòu)熱力學(xué)參數(shù)的敏感性問題[2-3]。熱力學(xué)參數(shù)隨溫度、局部地質(zhì)構(gòu)造等因素的變化，直接影響著地下工程圍巖的穩(wěn)定性。在實(shí)際工程中，若能分析出各參數(shù)的敏感性和影響程度，有選擇地對高敏感性熱力學(xué)參數(shù)予以精確測定，抓住主要參數(shù)，使熱力學(xué)參數(shù)的確定更為有的放矢，同時也為洞室整體布局和支護(hù)結(jié)構(gòu)設(shè)計提供依據(jù)。

目前，獲取巖體、支護(hù)結(jié)構(gòu)有效熱力學(xué)參數(shù)的主要途徑是采用等效概化的方法，結(jié)合現(xiàn)場參數(shù)信息的反分析進(jìn)行的[4- 6]。朱維申和何滿潮[7]提出了敏感性分析方法，對影響圍巖穩(wěn)定的一系列參數(shù)進(jìn)行了單因素敏感性分析。黃書嶺等[8]提出基于敏感度熵權(quán)的屬性識別綜合評價模型，為參數(shù)敏感性分析提供了一種新的思路。李曉靜等[9]以瑯琊山抽水蓄能電站地下廠房為工程背景，選取對地下洞室穩(wěn)定性影響較為重要的4個參數(shù)(變形模量、洞室埋深、主廠房高度、側(cè)壓力系數(shù))進(jìn)行大量塑性數(shù)值模擬分析并對位移進(jìn)行了參數(shù)敏感性分析。侯哲生等[10]利用非線性彈塑性有限元法，研究了金川二礦區(qū)某巷道圍巖力學(xué)參數(shù)對變形的敏感性，得到了不同參數(shù)對變形的敏感性。聶衛(wèi)平等[11]采用基于彈塑性有限元的洞室穩(wěn)定性參數(shù)敏感性灰關(guān)聯(lián)分析法，對地下洞室穩(wěn)定影響參數(shù)進(jìn)行了敏感性分析。

前人的研究成果主要集中在圍巖力學(xué)參數(shù)的敏感性分析及其可能出現(xiàn)的變化對圍巖穩(wěn)定、變形的影響程度上。工程實(shí)踐表明，地下工程的深埋地段存在如高地溫、滲透水流等，影響巖體穩(wěn)定、支護(hù)結(jié)構(gòu)力學(xué)特性的參數(shù)較多，由有限的信息求解眾多參數(shù)的敏感性問題仍存在一定的困難。以往的研究成果都是通過類比和經(jīng)驗等方法確定力學(xué)參數(shù)的取值范圍[12-14]，并對其進(jìn)行敏感性分析。地下洞室圍巖及其支護(hù)結(jié)構(gòu)的熱力學(xué)參數(shù)是數(shù)值模擬分析的基礎(chǔ)，常溫時材料的熱力學(xué)參數(shù)變化很小，一般可以認(rèn)為是常數(shù)，但溫度較高時會發(fā)生明顯的變化，從而影響圍巖和支護(hù)結(jié)構(gòu)的應(yīng)力及其穩(wěn)定性。因此，開展圍巖和支護(hù)結(jié)構(gòu)熱力學(xué)參數(shù)的變化對其應(yīng)力的敏感性研究具有極其重要的現(xiàn)實(shí)意義。

1 熱力學(xué)參數(shù)敏感性分析

熱力學(xué)參數(shù)的敏感性分析是為了得到參數(shù)隨溫度變化時對噴層應(yīng)力和溫度的敏感性的大小。因應(yīng)力場和溫度場耦合的復(fù)雜性，難以求得噴層和圍巖應(yīng)力的顯示表達(dá)式，故本文采用ABAQUS有限元程序計算噴層的溫度場和應(yīng)力場，結(jié)合新疆某水電站引水隧洞工程，選取拱頂、拱肩和側(cè)墻中部的環(huán)向應(yīng)力，對熱力學(xué)參數(shù)的敏感性進(jìn)行分析。

1.1 工程概況

新疆某水電站引水隧洞為城門洞形，寬4.60 m，高5.23 m，直墻高3.78 m，拱頂半徑2.55 m。引水隧洞示意見圖1。引水隧洞圍巖巖石較堅硬，呈中厚層狀，圍巖類別為Ⅲ類。引水隧洞施工過程中存在高地溫問題，圍巖開挖最高溫度105 ℃，而運(yùn)行時的溫度低至0～5 ℃，圍巖穩(wěn)定和支護(hù)結(jié)構(gòu)的設(shè)計受到溫度場和應(yīng)力場耦合的影響。為研究圍巖與支護(hù)結(jié)構(gòu)熱力學(xué)參數(shù)的變化對應(yīng)力的敏感性，本文以引水隧洞高地溫段圍巖熱力學(xué)參數(shù)對圍巖噴層環(huán)向應(yīng)力的敏感性分析作為依據(jù)，結(jié)合現(xiàn)場監(jiān)測數(shù)據(jù)，重點(diǎn)分析導(dǎo)熱系數(shù)、比熱、對流系數(shù)和線膨脹系數(shù)隨溫度變化時對支護(hù)結(jié)構(gòu)應(yīng)力的敏感性，明確不同參數(shù)對支護(hù)結(jié)構(gòu)應(yīng)力的影響程度。

圖1 引水隧洞示意(單位：mm)

1.2 敏感性分析方法

敏感性分析首先是建立敏感分析的系統(tǒng)模型，即系統(tǒng)特性P與各因素(a1，a2，…，an)之間的函數(shù)關(guān)系P=f(a1，a2，…，an)。這種函數(shù)關(guān)系，盡可能用解析式表示。

建立系統(tǒng)模型后，根據(jù)所要討論的具體問題給出基準(zhǔn)參數(shù)集。如果要分析某地下洞室圍巖應(yīng)力對其巖石熱力學(xué)參數(shù)變化的敏感性，則該工程巖石熱力學(xué)參數(shù)的推薦值可取為基準(zhǔn)參數(shù)集?；鶞?zhǔn)參數(shù)集確定后，就可對各參數(shù)進(jìn)行敏感性分析。分析參數(shù)a*對系統(tǒng)特性的影響時，可令其他各參數(shù)取基準(zhǔn)值且固定不變，而令a*在其可能的范圍內(nèi)變化，則系統(tǒng)特性P表現(xiàn)為

=φk(ak)

(1)

利用式(1)繪制出系統(tǒng)特性P與參數(shù)ak關(guān)系曲線，通過P-ak曲線可大致了解系統(tǒng)特性P對參數(shù)ak變化的敏感性。

在實(shí)際系統(tǒng)中，決定系統(tǒng)特性的各參數(shù)往往是不同的物理量，憑借以上的分析，無法對各參數(shù)之間的敏感程度進(jìn)行比較。因此，有必要對各參數(shù)進(jìn)行無量綱化的處理[15]，即將系統(tǒng)特性P*的相對誤差δPk=|ΔPk|/P*與參數(shù)a*的相對誤差δak=|Δak|/a*的比值定義為參數(shù)ak的敏感函數(shù)Sk(ak)，即

(2)

在|Δak|/a*較小的情況下，Sk(ak)可近似地表示為

(3)

(4)

表1 敏感性分析方案及其參數(shù)

注：參數(shù)計算中，圍巖線膨脹系數(shù)的第1次取值為0.5×10-6/℃；第2次取值為1.0×10/℃；其他均為2×10-6/℃。

1.3 參與敏感性分析的參數(shù)及其分析方案

由于高地溫引水隧洞圍巖和支護(hù)結(jié)構(gòu)材料及受力的復(fù)雜性，在計算中所涉及到的參數(shù)比較多，不僅包括力學(xué)參數(shù)(彈性模量、泊松比、線膨脹系數(shù)等)，還包括熱學(xué)參數(shù)(導(dǎo)熱系數(shù)、比熱、對流系數(shù)等)。從隧洞工程的設(shè)計與支護(hù)結(jié)構(gòu)的設(shè)計出發(fā)，分析重要的敏感參數(shù)是非常重要的。通過對熱力學(xué)參數(shù)如線膨脹系數(shù)、導(dǎo)熱系數(shù)、比熱、對流系數(shù)等進(jìn)行敏感性分析，從而確定出重要的敏感性的參數(shù)，為洞室的穩(wěn)定性分析和支護(hù)結(jié)構(gòu)設(shè)計提供依據(jù)。表1給出了敏感性分析方案及其參數(shù)，計算步數(shù)表示對某參數(shù)進(jìn)行敏感性分析時，其他參數(shù)保持初值不變，在此參數(shù)的變化范圍內(nèi)，按其所設(shè)計的步長分析所需的計算次數(shù)。

2 敏感性分析

2.1 線膨脹系數(shù)

2.1.1 圍巖

根據(jù)表1方案，僅改變圍巖線膨脹系數(shù)，其他參數(shù)均采用初始值進(jìn)行計算。由于改變的是圍巖的線膨脹系數(shù)，對溫度場沒有影響，噴層的溫度不隨圍巖線膨脹系數(shù)的改變而改變，所以圍巖噴層的溫度根據(jù)現(xiàn)場實(shí)測值(內(nèi)側(cè)6.4 ℃，外側(cè)31.7 ℃)進(jìn)行分析，而噴層的應(yīng)力是隨圍巖線膨脹系數(shù)的變化而變化的，從而可分析圍巖線膨脹系數(shù)變化時噴層不同部位應(yīng)力的敏感性。

首先采用上述敏感性分析方法得出敏感度函數(shù)，由圍巖線膨脹系數(shù)和噴層環(huán)向應(yīng)力的關(guān)系曲線采用曲線擬合的方法，建立噴層拱頂環(huán)向應(yīng)力σ與線膨脹系數(shù)a的函數(shù)關(guān)系為：σ=0.083 6α+5.24，可得噴層拱頂環(huán)向應(yīng)力敏感度函數(shù)Sα為

Sα=(0.083 6α)/(0.083 6α+5.24)

(5)

綜上分析，可得圍巖線膨脹系數(shù)變化時噴層拱頂、拱肩、側(cè)墻中部環(huán)向應(yīng)力的變化情況，并可以分析得出圍巖線膨脹系數(shù)變化時噴層不同部位應(yīng)力的變化及其與敏感性的關(guān)系。不同圍巖線膨脹系數(shù)時噴層應(yīng)力見圖2。

圖2 不同圍巖線膨脹系數(shù)時噴層應(yīng)力變化

從圖2可知，圍巖線膨脹系數(shù)從0.5×10-6/℃提高到10.0×10-6/℃，增大了20倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從5.37 MPa增大到5.85 MPa，增大了8.94%；噴層拱肩的環(huán)向應(yīng)力從3.1 MPa減小到2.96 MPa，減小了4.5%；噴層側(cè)墻中部的環(huán)向應(yīng)力從1.42 MPa減小到0.82 MPa，減小了40.14%。隨著圍巖線膨脹系數(shù)的增大，噴層拱頂環(huán)向應(yīng)力呈現(xiàn)出線性增加的趨勢，拱肩、側(cè)墻中部環(huán)向應(yīng)力都呈現(xiàn)出線性減小的趨勢。

從敏感性分析可以得出，隨著圍巖線膨脹系數(shù)的增大，噴層拱頂、拱肩、側(cè)墻中部環(huán)向應(yīng)力敏感度都呈現(xiàn)線性增加的趨勢，這說明隨著圍巖線膨脹系數(shù)在0.5×10-6～10.0×10-6/℃范圍內(nèi)增大時，噴層環(huán)向應(yīng)力越來越敏感。

2.1.2 噴層

根據(jù)表1敏感性分析方案設(shè)計，僅改變噴層線膨脹系數(shù)，其他參數(shù)均采用初值進(jìn)行計算。由于改變的是噴層的線膨脹系數(shù)，對溫度場沒有影響，噴層的溫度不隨噴層線膨脹系數(shù)的改變而改變，所以噴層的溫度根據(jù)現(xiàn)場實(shí)測值(內(nèi)側(cè)6.4 ℃，外側(cè)31.7 ℃)進(jìn)行分析，而噴層的應(yīng)力是隨噴層線膨脹系數(shù)的變化而變化的，從而可分析噴層線膨脹系數(shù)變化時噴層不同部位應(yīng)力的敏感性。不同噴層線膨脹系數(shù)時噴層應(yīng)力見圖3。

圖3 不同噴層線膨脹系數(shù)時噴層應(yīng)力變化

從圖3可知，噴層線膨脹系數(shù)從2.0×10-6/℃提高到12.0×10-6/℃，增大了6倍，拱頂?shù)沫h(huán)向應(yīng)力從0.72 MPa增大到6.87 MPa，增大了8.54倍；拱肩的環(huán)向應(yīng)力從-1.89 MPa增大到4.23 MPa，由壓應(yīng)力變?yōu)槔瓚?yīng)力；側(cè)墻中部的環(huán)向應(yīng)力從-3.82 MPa 增大到2.29 MPa，由壓應(yīng)力變?yōu)槔瓚?yīng)力。隨著噴層線膨脹系數(shù)的增大，噴層拱頂、拱肩、側(cè)墻中部環(huán)向應(yīng)力都呈現(xiàn)出線性增加的趨勢。

從敏感性分析可以得出，隨著噴層線膨脹系數(shù)的增大，噴層拱頂環(huán)向應(yīng)力敏感性呈現(xiàn)減小的趨勢，且減小幅度越來越小，逐步趨于穩(wěn)定；噴層拱肩、側(cè)墻中部噴層環(huán)向應(yīng)力敏感性在壓應(yīng)力變?yōu)槔瓚?yīng)力時出現(xiàn)突變。

2.2 導(dǎo)熱系數(shù)

2.2.1 圍巖導(dǎo)熱系數(shù)

2.2.1.1溫度

圍巖導(dǎo)熱系數(shù)不同時噴層內(nèi)外側(cè)溫度變化見圖4。從圖4可知，隨著圍巖導(dǎo)熱系數(shù)的增加，噴層的內(nèi)外側(cè)的溫度都呈現(xiàn)上升的趨勢，這是由于圍巖導(dǎo)熱系數(shù)越大，圍巖傳導(dǎo)熱系數(shù)的能力越強(qiáng)，噴層在過水時散熱能力就相對降低，所以噴層的溫度就會增大。圍巖導(dǎo)熱系數(shù)從5 W/(m·℃)升為30 W/(m·℃)，升高了6倍，噴層內(nèi)側(cè)溫度從6.3 ℃升為7.6 ℃，升高了20.6%；噴層外側(cè)溫度從30.1 ℃升高到55.3 ℃，升高了83.72%?？梢钥闯?，噴層外側(cè)的溫度變化遠(yuǎn)遠(yuǎn)大于內(nèi)側(cè)，主要是噴層內(nèi)側(cè)直接過水的結(jié)果。

圖4 不同圍巖導(dǎo)熱系數(shù)噴層溫度變化

圖5 不同圍巖導(dǎo)熱系數(shù)時噴層應(yīng)力變化

2.2.1.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見圖5。圍巖導(dǎo)熱系數(shù)從5 W/(m·℃)升到30 W/(m·℃)，增大了6倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從5.48 MPa增大到6.63 MPa，增大了20.98%；噴層拱肩的環(huán)向應(yīng)力從2.95 MPa增大到3.78 MPa，增大了28.14%；噴層側(cè)墻中部的環(huán)向應(yīng)力從1 MPa增大到1.69 MPa，增大了69%。噴層拱頂環(huán)向應(yīng)力隨著圍巖導(dǎo)熱系數(shù)的增大呈現(xiàn)增加的趨勢，開始呈線性增加，幅度較大，后來慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減??；噴層拱肩和側(cè)墻中部的規(guī)律和拱頂相似。

2.2.2 噴層導(dǎo)熱系數(shù)

2.2.2.1溫度

噴層導(dǎo)熱系數(shù)不同時噴層內(nèi)外側(cè)溫度變化見圖6。從圖6可知，隨著噴層導(dǎo)熱系數(shù)的增加，噴層的外側(cè)的溫度呈現(xiàn)降低的趨勢，這是由于噴層導(dǎo)熱系數(shù)越大，噴層傳導(dǎo)熱系數(shù)的能力越強(qiáng)。當(dāng)噴層的導(dǎo)熱系數(shù)很低時，就相當(dāng)于一個隔熱層，噴層外側(cè)的溫度就會很高，所以噴層外側(cè)溫度隨著噴層導(dǎo)熱系數(shù)的增大而降低；對于噴層內(nèi)側(cè)，由于噴層導(dǎo)熱系數(shù)的增加，能更好地把圍巖的熱量傳遞給噴層，所以噴層內(nèi)側(cè)溫度升高。噴層導(dǎo)熱系數(shù)從0.5 W/(m·℃)升為3 W/(m·℃)，升高了6倍；噴層內(nèi)側(cè)溫度從6.1℃升為7.1℃，升高了16.39%；噴層外側(cè)溫度從60.5 ℃降低到25.9 ℃，降低了57.19%?？梢钥闯?，噴層外側(cè)的溫度變化遠(yuǎn)遠(yuǎn)大于內(nèi)側(cè)，主要是噴層內(nèi)側(cè)直接過水的結(jié)果。

圖6 不同噴層導(dǎo)熱系數(shù)時噴層外側(cè)溫度變化

2.2.2.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見圖7。從圖7可知，噴層導(dǎo)熱系數(shù)從0.5 W/(m·℃)升到3.0 W/(m·℃)，增大6倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從7.06 MPa減小到5.31 MPa，減小了22.78%；噴層拱肩的環(huán)向應(yīng)力從4.58 MPa減小到2.47 MPa，減小了46.1%；噴層側(cè)墻中部的環(huán)向應(yīng)力從2.75 MPa減小到0.38 MPa，減小了86.18%。噴層拱頂環(huán)向應(yīng)力隨著噴層導(dǎo)熱系數(shù)的增大呈現(xiàn)減小的趨勢，開始呈線性減小，幅度較大，后來慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減小；噴層拱肩和側(cè)墻中部和拱頂相似。

圖7 不同導(dǎo)熱系數(shù)時噴層應(yīng)力變化

2.3 比熱

2.3.1 圍巖比熱敏感性

2.3.1.1溫度

圖8 不同圍巖比熱噴層溫度變化

圍巖比熱不同時噴層內(nèi)外側(cè)溫度變化見圖8。從圖8可知，隨著圍巖比熱的增加，噴層的內(nèi)側(cè)和外側(cè)的溫度都呈現(xiàn)上升的趨勢，這是由于圍巖比熱越大，圍巖在降低1 ℃時釋放的熱量也就越多，圍巖溫度降低的幅度就降小，傳遞給噴層的溫度也就會越高。具體來講，圍巖比熱從200 J/(kg·℃)升為1 400 J/(kg·℃)，圍巖比熱升高了7倍，噴層內(nèi)側(cè)溫度從6.3 ℃升為6.9 ℃，升高了9.52%；噴層外側(cè)溫度從29.3 ℃升高到40.4 ℃，升高了37.88%。

2.3.1.2應(yīng)力在上述溫度分布情況下，噴層的應(yīng)力變化見圖9。從圖9可知，圍巖比熱從200 J/(kg·℃)升為1 400 J/(kg·℃)，圍巖比熱升高了7倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從5.9 MPa增大到6.48 MPa，增大了9.83%；噴層拱肩的環(huán)向應(yīng)力從2.73 MPa增大到3.31 MPa，增大了21.24%；噴層側(cè)墻中部的環(huán)向應(yīng)力從-0.01 MPa增大到1.40 MPa，從壓應(yīng)力轉(zhuǎn)化為拉應(yīng)力。噴層拱頂環(huán)向應(yīng)力隨著圍巖比熱的增大呈現(xiàn)增加的趨勢，但增加趨勢開始呈線性增加，幅度較大，后來慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減?。粐妼庸凹缫?guī)律和拱頂相似；噴層側(cè)墻中部環(huán)向應(yīng)力隨著圍巖比熱的增大呈現(xiàn)增加的趨勢，由于噴層應(yīng)力出現(xiàn)了從壓應(yīng)力到拉應(yīng)力的轉(zhuǎn)變，其敏感度起始點(diǎn)始于突變點(diǎn)，但總的來說，其敏感度逐漸降低。

圖9 不同圍巖比熱時噴層應(yīng)力變化

2.3.2 噴層比熱敏感性

2.3.2.1溫度

噴層比熱不同時噴層內(nèi)外側(cè)溫度變化見表2。從表2可知，隨著噴層比熱的增加，噴層比熱從200 J/(kg·℃)升到1 400 J/(kg·℃)，升高了7倍，噴層內(nèi)側(cè)溫度一直為6.7 ℃；噴層外側(cè)溫度從38.0 ℃升高到38.1 ℃，升高了0.26%。噴層的內(nèi)側(cè)和外側(cè)的溫度幾乎沒有變，這是由于噴層特別薄，比熱對其影響不明顯。

表3 不同噴層比熱時噴層關(guān)鍵部位環(huán)向應(yīng)力 MPa

2.3.2.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見表3。從表3可知，噴層比熱從200 J/(kg·℃)升到1 400 J/(kg·℃)，升高了7倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從5.93 MPa減小到5.91 MPa，減小了0.34%；噴層拱肩的環(huán)向應(yīng)力從3.22 MPa減小到3.21 MPa，減小了0.31%；噴層側(cè)中的環(huán)向應(yīng)力從1.22 MPa減小到1.21 MPa，減小了0.82%。

2.4 對流系數(shù)

2.4.1 溫度

對流系數(shù)不同時噴層內(nèi)外側(cè)溫度變化見圖10。從圖10可知，隨著圍巖對流系數(shù)的增加，噴層的內(nèi)外側(cè)的溫度都呈現(xiàn)降低的趨勢。圍巖對流系數(shù)從50 W/(m2·℃)升到5 000 W/(m2·℃)，升高了100倍，噴層內(nèi)側(cè)溫度從11.6 ℃降為5.1 ℃，降低了56.03%；噴層外側(cè)溫度從42.2 ℃降為36.4 ℃，降低了13.74%。這是由于圍巖對流系數(shù)越大，圍巖與水換熱能力越強(qiáng)，噴層的溫度就會降低。

圖10 不同對流系數(shù)噴層溫度變化

2.4.2 應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見圖11。

圖11 不同對流系數(shù)時噴層應(yīng)力變化

從圖11可知，對流系數(shù)的增大對噴層的環(huán)向應(yīng)力影響比較明顯。對流系數(shù)從50 W/(m2·℃)升為5 000 W/(m2·℃)，升高了100倍，噴層拱頂?shù)沫h(huán)向應(yīng)力從4.57 MPa增大到6.39 MPa，增大了39.82%；噴層拱肩的環(huán)向應(yīng)力從1.93 MPa增大到3.65 MPa，增大了89.12%；噴層側(cè)墻中部的環(huán)向應(yīng)力從-0.01 MPa增大到1.62 MPa，由壓應(yīng)力變化為拉應(yīng)力。噴層拱頂環(huán)向應(yīng)力隨著對流系數(shù)的增大呈現(xiàn)出增大的趨勢。對流系數(shù)為50～1 000 W/(m2·℃)時對噴層的環(huán)向應(yīng)力影響比較大，1 000 W/(m2·℃)以上影響越來越小了。故其敏感度呈現(xiàn)減小的趨勢，且開始減小較快，后來慢慢趨于穩(wěn)定；噴層拱肩和邊墻的環(huán)向應(yīng)力和拱頂?shù)淖兓?guī)律相似。

表4 各參數(shù)的敏感度因子

2.5 各參數(shù)的敏感度比較

為了定量計算各參數(shù)的敏感度因子，將噴層拱頂環(huán)向應(yīng)力設(shè)為計算分析依據(jù)，對噴層拱頂環(huán)向應(yīng)力的敏感度因子進(jìn)行匯總，比較分析各參數(shù)的敏感程度。各參數(shù)的敏感度因子見表4。由表4可知，圍巖和噴層的熱力學(xué)參數(shù)對噴層拱頂環(huán)向應(yīng)力的敏感度排序為：噴層線膨脹系數(shù)>噴層導(dǎo)熱系數(shù)>圍巖導(dǎo)熱系數(shù)>對流系數(shù)>圍巖的線膨脹系數(shù)>圍巖的比熱>噴層的比熱。敏感度為0.06時為不敏感參數(shù)，本文認(rèn)為，敏感度因子小于0.1時參數(shù)不敏感。因此，圍巖線膨脹系數(shù)、圍巖比熱和噴層比熱為不敏感參數(shù)。

3 結(jié) 語

本文以新疆某水電站引水隧洞高地溫段為研究背景，對圍巖和噴層結(jié)構(gòu)的熱力學(xué)參數(shù)與噴層環(huán)向應(yīng)力的關(guān)系進(jìn)行了分析，并對影響噴層溫度場和應(yīng)力場的圍巖和噴層熱學(xué)參數(shù)進(jìn)行了敏感性分析，得出以下結(jié)論：

(1)隨著圍巖線膨脹系數(shù)的增加，噴層拱頂?shù)睦瓚?yīng)力線性增大；拱肩和邊墻的拉應(yīng)力呈線性減小的趨勢；各部位的敏感度隨著圍巖線膨脹系數(shù)的增大都呈現(xiàn)線性增加的趨勢；隨著噴層線膨脹系數(shù)的增加，噴層拱頂拉應(yīng)力線性增大，而拱肩和側(cè)墻中部出現(xiàn)壓應(yīng)力轉(zhuǎn)化為拉應(yīng)力的狀態(tài)；噴層拱頂環(huán)向應(yīng)力敏感度隨著噴層線膨脹系數(shù)的增大呈現(xiàn)減小趨勢，且減小趨勢越來越緩慢，拱肩和側(cè)墻中部噴層應(yīng)力從壓應(yīng)力轉(zhuǎn)化為拉應(yīng)力，其敏感度在變化時出現(xiàn)突變。

(2)隨著圍巖導(dǎo)熱系數(shù)的增加，噴層拱頂拉應(yīng)力呈現(xiàn)增大的趨勢，開始呈線性增加，幅度較大，后來慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減?。粐妼庸凹绾蛡?cè)墻中部的規(guī)律和拱頂相似。隨著噴層導(dǎo)熱系數(shù)增加，噴層拱頂拉應(yīng)力呈現(xiàn)減小的趨勢，開始呈線性減小，幅度較大，后來慢慢趨于平緩，呈現(xiàn)拋物線的形式，先增大后減小；噴層拱肩和側(cè)墻中部和拱頂相似。

(3)隨著圍巖比熱的增加，噴層拱頂拉應(yīng)力呈現(xiàn)增大的趨勢，開始呈線性增加，幅度較大，后來慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減小；噴層拱肩規(guī)律和拱頂相似；側(cè)墻中部由于出現(xiàn)壓應(yīng)力到拉應(yīng)力的轉(zhuǎn)變，其敏感度起始點(diǎn)始于突變點(diǎn)，但總的來說，其敏感度逐漸降低。由于噴層特別薄，噴層比熱增加對噴層的溫度幾乎沒有影響，噴層比熱的變化不影響噴層的環(huán)向應(yīng)力。

(4)隨著對流系數(shù)的增加，噴層拱頂拉應(yīng)力呈現(xiàn)增大的趨勢。對流系數(shù)為50～1 000 W/(m2·℃)時，對噴層的環(huán)向應(yīng)力影響比較大，1 000 W/(m2·℃)以上影響越來越小了，其敏感度呈現(xiàn)出減小的趨勢，且開始減小較快，后來慢慢趨于穩(wěn)定；噴層拱肩和邊墻的環(huán)向應(yīng)力和拱頂?shù)南嗨啤?/p>

(5)通過敏感度因子的計算，圍巖和噴層的熱力學(xué)參數(shù)對噴層拱頂環(huán)向應(yīng)力的敏感度排序為：噴層線膨脹系數(shù)>噴層導(dǎo)熱系數(shù)>圍巖的導(dǎo)熱系數(shù)>對流系數(shù)>圍巖的線膨脹系數(shù)>圍巖的比熱>噴層的比熱。圍巖線膨脹系數(shù)、圍巖比熱和噴層比熱為不敏感參數(shù)。

需要指出的是，本文所得到的結(jié)論綜合考慮了多種因素對應(yīng)力和溫度的影響。現(xiàn)場監(jiān)測數(shù)據(jù)表明，支護(hù)結(jié)構(gòu)的應(yīng)力具有一定的變異特性，本文分析未考慮支護(hù)結(jié)構(gòu)應(yīng)力的變異特性，需要做進(jìn)一步的研究。

[1] 王輝, 陳衛(wèi)忠. 嘎隆拉隧道圍巖力學(xué)參數(shù)對變形的敏感性分析[J]. 巖土工程學(xué)報, 2012, 34(8)： 1548- 1553．

[2] BEIKI M, BASHARI A, MAJDI A. Genetic programming approach for estimating the deformation modulus of rock mass using sensitivity analysis by neural network[J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(1)： 1091- 1103．

[3] 朱維申, 章光. 節(jié)理巖體參數(shù)對圍巖破損區(qū)影響的敏感性分析[J]. 地下空間, 1994, 14(1)： 10- 15．

[4] 馮夏庭. 智能巖石力學(xué)導(dǎo)論[M]. 北京： 科學(xué)出版社, 2000．

[5] SAKURAI S, TAKEUCHI K. Back analysis of measured displacement of tunnel[J]. Rock Mechanics and Rock Engineering, 1983, 16(3)： 173- 180．

[6] 楊林德, 馮紫良, 朱合華, 等. 巖土工程問題的反演理論與工程實(shí)踐[M]. 北京： 科學(xué)出版社, 1996．[7] 朱維申, 何滿潮. 復(fù)雜條件下圍巖穩(wěn)定性與巖體動態(tài)施工力學(xué)[M]. 北京： 科學(xué)出版社, 1995．

[8] 黃書嶺, 馮夏庭, 張傳慶. 巖體力學(xué)參數(shù)的敏感性綜合評價分析方法研究[J]. 巖石力學(xué)與工程學(xué)報, 2008, 27(z1)： 2624- 2630.

[9] 李曉靜, 朱維申, 向建, 等. 考慮參數(shù)影響的系統(tǒng)分析方法及其應(yīng)用[J]. 巖土工程學(xué)報, 2005, 27(10)： 1207- 1210．

[10] 侯哲生, 李曉, 王思敬, 等. 金川二礦某巷道圍巖力學(xué)參數(shù)對變形的敏感性分析[J]. 巖石力學(xué)與工程學(xué)報, 2005, 24(3)： 406- 410．

[11] 聶衛(wèi)平, 徐衛(wèi)亞, 周先齊. 基于三維彈塑性有限元的洞室穩(wěn)定性參數(shù)敏感性灰關(guān)聯(lián)分析[J]. 巖石力學(xué)與工程學(xué)報, 2009, 28(z2)： 3885- 3893．

[12] LI Shihui, YANG Jie, HAO Weidong, et al. Intelligent back-analysis of displacements monitored in tunneling[J]. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(2)： 1118- 1127．

[13] 黃書嶺, 馮夏庭, 張傳慶. 巖體力學(xué)參數(shù)的敏感性綜合評價分析方法研究[J]. 巖石力學(xué)與工程學(xué)報, 2008, 27(z1)： 2624- 2 630．

[14] 章光, 朱維申. 參數(shù)敏感性分析與試驗方案優(yōu)化[J]. 巖土力學(xué), 1993, 14 (1)： 51- 57．

[15] 楊蒙, 譚躍虎, 李二兵, 等. 基于敏感性分析的圍巖力學(xué)參數(shù)反演方法研究[J]. 地下空間與工程學(xué)報, 1993, 14 (1)： 51- 57．

SensitivityAnalysisofThermodynamicsParametersofSurroundingRockandSupportingStructureinHighGroundTemperatureDiversionTunnel

JIANG Haibo1, WU Peng2, ZHANG Jun2

(1. School of Hydraulic Engineering, Shihezi University, Shihezi 832000, Xinjiang, China;2. Xinjiang Survey and Design Institute for Water Resources and Hydropower, Urumqi 830000, Xinjiang, China)

Taking the section in high temperature field of the diversion tunnel of a hydropower station in Xinjiang as study object, the relationships between the mechanical parameters (thermal conductivity, specific heat, convection coefficient and linear expansion coefficient) of surrounding rock and shotcrete layer structure and the circumferential stress of shotcrete layer are studied based on on-site stress and temperature monitoring data, and the sensitivity of thermodynamic parameters of surrounding rock and shotcrete layer to shotcrete stress is also studied by using parameter sensitivity analysis method. The results show that the main influence factors on the circumferential stress of tunnel vault are ordered as linear expansion coefficient of shotcrete layer, thermal conductivity of shotcrete layer, thermal conductivity of surrounding rock, convection coefficient, linear expansion coefficient of surrounding rock, specific heat of surrounding rock and specific heat of shotcrete layer, in which, the linear expansion coefficient of shotcrete layer, thermal conductivity of shotcrete layer, thermal conductivity of surrounding rock and convection coefficient are more sensitive and the linear expansion coefficient of surrounding rock, specific heat of surrounding rock and specific heat of shotcrete layer are not sensitive to the circumferential stress of tunnel vault．

diversion tunnel; mechanical parameter of surrounding rock; temperature of surrounding rock; circumferential stress; sensitivity analysis

TV223.1(245)

A

0559- 9342(2017)09- 0031- 08

2017- 05- 30

國家自然科學(xué)基金項目(51408377)；石河子大學(xué)杰出青年科技人才培育計劃(2015ZRKXJQ06)

姜海波(1982—)，男，湖南長沙人，副教授，博士，主要從事地下洞室抗凍及穩(wěn)定性研究工作．

(責(zé)任編輯楊 健)