姜海波，吳 鵬，張 軍

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

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

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

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

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

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

0 引 言

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

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

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

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

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

1.1 工程概況

新疆某水電站引水隧洞為城門(mén)洞形，寬4.60 m，高5.23 m，直墻高3.78 m，拱頂半徑2.55 m。引水隧洞示意見(jiàn)圖1。引水隧洞圍巖巖石較堅(jiān)硬，呈中厚層狀，圍巖類別為Ⅲ類。引水隧洞施工過(guò)程中存在高地溫問(wèn)題，圍巖開(kāi)挖最高溫度105 ℃，而運(yùn)行時(shí)的溫度低至0～5 ℃，圍巖穩(wěn)定和支護(hù)結(jié)構(gòu)的設(shè)計(jì)受到溫度場(chǎng)和應(yīng)力場(chǎng)耦合的影響。為研究圍巖與支護(hù)結(jié)構(gòu)熱力學(xué)參數(shù)的變化對(duì)應(yīng)力的敏感性，本文以引水隧洞高地溫段圍巖熱力學(xué)參數(shù)對(duì)圍巖噴層環(huán)向應(yīng)力的敏感性分析作為依據(jù)，結(jié)合現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)，重點(diǎn)分析導(dǎo)熱系數(shù)、比熱、對(duì)流系數(shù)和線膨脹系數(shù)隨溫度變化時(shí)對(duì)支護(hù)結(jié)構(gòu)應(yīng)力的敏感性，明確不同參數(shù)對(duì)支護(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ù)所要討論的具體問(wèn)題給出基準(zhǔn)參數(shù)集。如果要分析某地下洞室圍巖應(yīng)力對(duì)其巖石熱力學(xué)參數(shù)變化的敏感性，則該工程巖石熱力學(xué)參數(shù)的推薦值可取為基準(zhǔn)參數(shù)集?；鶞?zhǔn)參數(shù)集確定后，就可對(duì)各參數(shù)進(jìn)行敏感性分析。分析參數(shù)a*對(duì)系統(tǒng)特性的影響時(shí)，可令其他各參數(shù)取基準(zhǔn)值且固定不變，而令a*在其可能的范圍內(nèi)變化，則系統(tǒng)特性P表現(xiàn)為

=φk(ak)

(1)

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

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

(2)

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

(3)

(4)

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

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

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

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

2 敏感性分析

2.1 線膨脹系數(shù)

2.1.1 圍巖

根據(jù)表1方案，僅改變圍巖線膨脹系數(shù)，其他參數(shù)均采用初始值進(jìn)行計(jì)算。由于改變的是圍巖的線膨脹系數(shù)，對(duì)溫度場(chǎng)沒(méi)有影響，噴層的溫度不隨圍巖線膨脹系數(shù)的改變而改變，所以圍巖噴層的溫度根據(jù)現(xiàn)場(chǎng)實(shí)測(cè)值(內(nèi)側(cè)6.4 ℃，外側(cè)31.7 ℃)進(jìn)行分析，而噴層的應(yīng)力是隨圍巖線膨脹系數(shù)的變化而變化的，從而可分析圍巖線膨脹系數(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ù)變化時(shí)噴層拱頂、拱肩、側(cè)墻中部環(huán)向應(yīng)力的變化情況，并可以分析得出圍巖線膨脹系數(shù)變化時(shí)噴層不同部位應(yīng)力的變化及其與敏感性的關(guān)系。不同圍巖線膨脹系數(shù)時(shí)噴層應(yīng)力見(jiàn)圖2。

圖2 不同圍巖線膨脹系數(shù)時(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)出線性增加的趨勢(shì)，拱肩、側(cè)墻中部環(huán)向應(yīng)力都呈現(xiàn)出線性減小的趨勢(shì)。

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

2.1.2 噴層

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

圖3 不同噴層線膨脹系數(shù)時(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ì)。

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

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

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

2.2.1.1溫度

圍巖導(dǎo)熱系數(shù)不同時(shí)噴層內(nèi)外側(cè)溫度變化見(jiàn)圖4。從圖4可知，隨著圍巖導(dǎo)熱系數(shù)的增加，噴層的內(nèi)外側(cè)的溫度都呈現(xiàn)上升的趨勢(shì)，這是由于圍巖導(dǎo)熱系數(shù)越大，圍巖傳導(dǎo)熱系數(shù)的能力越強(qiáng)，噴層在過(guò)水時(shí)散熱能力就相對(duì)降低，所以噴層的溫度就會(huì)增大。圍巖導(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è)直接過(guò)水的結(jié)果。

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

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

2.2.1.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見(jiàn)圖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)增加的趨勢(shì)，開(kāi)始呈線性增加，幅度較大，后來(lái)慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減小；噴層拱肩和側(cè)墻中部的規(guī)律和拱頂相似。

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

2.2.2.1溫度

噴層導(dǎo)熱系數(shù)不同時(shí)噴層內(nèi)外側(cè)溫度變化見(jiàn)圖6。從圖6可知，隨著噴層導(dǎo)熱系數(shù)的增加，噴層的外側(cè)的溫度呈現(xiàn)降低的趨勢(shì)，這是由于噴層導(dǎo)熱系數(shù)越大，噴層傳導(dǎo)熱系數(shù)的能力越強(qiáng)。當(dāng)噴層的導(dǎo)熱系數(shù)很低時(shí)，就相當(dāng)于一個(gè)隔熱層，噴層外側(cè)的溫度就會(huì)很高，所以噴層外側(cè)溫度隨著噴層導(dǎo)熱系數(shù)的增大而降低；對(duì)于噴層內(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è)直接過(guò)水的結(jié)果。

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

2.2.2.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見(jiàn)圖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)減小的趨勢(shì)，開(kāi)始呈線性減小，幅度較大，后來(lái)慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減小；噴層拱肩和側(cè)墻中部和拱頂相似。

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

2.3 比熱

2.3.1 圍巖比熱敏感性

2.3.1.1溫度

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

圍巖比熱不同時(shí)噴層內(nèi)外側(cè)溫度變化見(jiàn)圖8。從圖8可知，隨著圍巖比熱的增加，噴層的內(nèi)側(cè)和外側(cè)的溫度都呈現(xiàn)上升的趨勢(shì)，這是由于圍巖比熱越大，圍巖在降低1 ℃時(shí)釋放的熱量也就越多，圍巖溫度降低的幅度就降小，傳遞給噴層的溫度也就會(huì)越高。具體來(lái)講，圍巖比熱從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)力變化見(jiàn)圖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)增加的趨勢(shì)，但增加趨勢(shì)開(kāi)始呈線性增加，幅度較大，后來(lái)慢慢趨于平緩，呈現(xiàn)出拋物線的形式，先增大后減?。粐妼庸凹缫?guī)律和拱頂相似；噴層側(cè)墻中部環(huán)向應(yīng)力隨著圍巖比熱的增大呈現(xiàn)增加的趨勢(shì)，由于噴層應(yīng)力出現(xiàn)了從壓應(yīng)力到拉應(yīng)力的轉(zhuǎn)變，其敏感度起始點(diǎn)始于突變點(diǎn)，但總的來(lái)說(shuō)，其敏感度逐漸降低。

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

2.3.2 噴層比熱敏感性

2.3.2.1溫度

噴層比熱不同時(shí)噴層內(nèi)外側(cè)溫度變化見(jiàn)表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è)的溫度幾乎沒(méi)有變，這是由于噴層特別薄，比熱對(duì)其影響不明顯。

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

2.3.2.2應(yīng)力

在上述溫度分布下，噴層的應(yīng)力變化見(jiàn)表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 對(duì)流系數(shù)

2.4.1 溫度

對(duì)流系數(shù)不同時(shí)噴層內(nèi)外側(cè)溫度變化見(jiàn)圖10。從圖10可知，隨著圍巖對(duì)流系數(shù)的增加，噴層的內(nèi)外側(cè)的溫度都呈現(xiàn)降低的趨勢(shì)。圍巖對(duì)流系數(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%。這是由于圍巖對(duì)流系數(shù)越大，圍巖與水換熱能力越強(qiáng)，噴層的溫度就會(huì)降低。

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

2.4.2 應(yīng)力

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

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

從圖11可知，對(duì)流系數(shù)的增大對(duì)噴層的環(huán)向應(yīng)力影響比較明顯。對(duì)流系數(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)力隨著對(duì)流系數(shù)的增大呈現(xiàn)出增大的趨勢(shì)。對(duì)流系數(shù)為50～1 000 W/(m2·℃)時(shí)對(duì)噴層的環(huán)向應(yīng)力影響比較大，1 000 W/(m2·℃)以上影響越來(lái)越小了。故其敏感度呈現(xiàn)減小的趨勢(shì)，且開(kāi)始減小較快，后來(lái)慢慢趨于穩(wěn)定；噴層拱肩和邊墻的環(huán)向應(yīng)力和拱頂?shù)淖兓?guī)律相似。

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

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

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

3 結(jié) 語(yǔ)

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

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

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

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

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

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

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

[1] 王輝, 陳衛(wèi)忠. 嘎隆拉隧道圍巖力學(xué)參數(shù)對(duì)變形的敏感性分析[J]. 巖土工程學(xué)報(bào), 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ù)對(duì)圍巖破損區(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] 楊林德, 馮紫良, 朱合華, 等. 巖土工程問(wèn)題的反演理論與工程實(shí)踐[M]. 北京： 科學(xué)出版社, 1996．[7] 朱維申, 何滿潮. 復(fù)雜條件下圍巖穩(wěn)定性與巖體動(dòng)態(tài)施工力學(xué)[M]. 北京： 科學(xué)出版社, 1995．

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

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

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

[11] 聶衛(wèi)平, 徐衛(wèi)亞, 周先齊. 基于三維彈塑性有限元的洞室穩(wěn)定性參數(shù)敏感性灰關(guān)聯(lián)分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 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] 黃書(shū)嶺, 馮夏庭, 張傳慶. 巖體力學(xué)參數(shù)的敏感性綜合評(píng)價(jià)分析方法研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2008, 27(z1)： 2624- 2 630．

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

[15] 楊蒙, 譚躍虎, 李二兵, 等. 基于敏感性分析的圍巖力學(xué)參數(shù)反演方法研究[J]. 地下空間與工程學(xué)報(bào), 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

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

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

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