伍文鋒，稅思梅

(1.長江勘測設(shè)計研究院 長江空間信息技術(shù)工程有限公司，湖北 武漢 430010；2.中國水電顧問集團成都勘測設(shè)計研究院 監(jiān)測中心，四川 成都 610072)

1 前 言

溪洛渡水電站地下廠房洞室群規(guī)模巨大，左、右岸地下廠房基本對稱布置于壩軸線上游庫區(qū)，廠房縱軸線方位分別為N24°W和N70°W。左、右岸地下廠區(qū)由主機間、安裝間、副廠房、主變室、9條壓力管道、9條母線洞、9條尾水管及尾水連接洞、3條尾水洞、2條出線井以及通排風(fēng)系統(tǒng)、防滲排水系統(tǒng)等組成，構(gòu)成以三大洞室為主體、縱橫交錯、上下分層的大規(guī)模的地下洞室群。三大洞室的設(shè)計開挖尺寸：廠房439.74m×75.6m×28.4m，主變室設(shè)計開挖尺寸349.289m×33.32m×19.8m，尾調(diào)室設(shè)計開挖尺寸317m×95.0m×25.0m(下部尺寸)。

地下廠區(qū)位于壩肩上游山體內(nèi)，水平埋深300～450m，垂直埋深340～480m。圍巖主要由P2β4、P2β5、P2β6層的斑狀玄武巖、含斑玄武巖、致密狀玄武巖及各層上部的角礫集塊熔巖組成，巖石堅硬，新鮮完整。

2 開挖支護方式

2.1 開挖方式

洞室開挖采用“新澳法”分層施工，開挖一層，支護一層。主廠房分十層開挖，主變室分四層開挖，尾調(diào)室分十五層開挖。主廠房和尾水調(diào)壓室先行施工，主變室在主廠房和尾水調(diào)壓室頂拱層開挖支護完成及主變室排風(fēng)洞開挖完成后滯后跟進。左岸主廠房7號機組開挖進度見圖1。

圖1 左岸主廠房7號機組開挖進度曲線

2.2 支護方式

洞室圍巖支護包括永久支護和臨時支護，針對不同的地質(zhì)情況，設(shè)計采用掛鋼筋網(wǎng)、噴射混凝土、安裝錨桿及預(yù)應(yīng)力錨索等多種支護方式。此外，對于局部不良地質(zhì)地段，采用中空式、自進式錨桿或鋼支撐等進行臨時支護。預(yù)應(yīng)力錨索均采用1 860MPa高強低松弛無粘結(jié)預(yù)應(yīng)力鋼絞線，設(shè)計安裝噸位分別為1 500kN、1 750kN。

3 回歸分析

溪洛渡地下廠房圍巖以Ⅱ類和Ⅲ1類為主，局部出現(xiàn)Ⅲ2類圍巖。地質(zhì)條件較好。從多點位移計成果分析，圍巖變形曲線大致呈臺階狀，存在突變點，與周圍巖體開挖爆破關(guān)系密切，累計變形絕大部分由開挖爆破引起，開挖間歇時變形緩慢增長，巖體卸荷量值很小。下面以主廠房7號機組斷面為例闡述回歸分析在溪洛渡地下廠房施工過程中的應(yīng)用。

3.1 回歸因子的初選

在溪洛渡地下廠房開挖過程中，影響邊墻變形的主要因素是施工開挖進度、開挖爆破方式、支護的方式及時間、時間效應(yīng)、溫度及相應(yīng)部位的地質(zhì)條件等。主廠房7號機組斷面實測位移量、溫度、開挖面高程與時間對應(yīng)關(guān)系見表1。

表1 主廠房7號機組斷面實測位移量、溫度、開挖面高程與時間對應(yīng)關(guān)系

由于地質(zhì)條件對測點的影響難以具體的量化，且認(rèn)為該部位的地質(zhì)條件屬于恒定常數(shù)，故將地質(zhì)條件的影響歸為常系數(shù)；由于開挖后支護的具體時間難以把握，且支護的進度基本與開挖進度存在一定的相關(guān)性；本次對地下廠房作的回歸分析主要考慮開挖進度、時間和溫度三個因子，可表達(dá)為：

y=yH+yT+yθ

(1)

3.1.1 開挖進度分量yH

開挖進度分量主要是因開挖造成應(yīng)力釋放而引起的邊墻變形，可表達(dá)為:

式中H——開挖面與儀器高程差值(隨時間變化)，i初選1，k初選2。

3.1.2 溫度分量yT

圍巖具有熱脹冷縮的特性，洞室溫度變化影響淺表層的變形，所以選擇洞室溫度作為溫度因子。安裝于孔口的多點位移計表筒內(nèi)裝有溫度計，可測量洞室溫度，所以選擇溫度計測值作溫度因子。圍巖變形與溫度呈線性關(guān)系，溫度分量表達(dá)為：

yT=bT

3.1.3 時效分量yθ

地下廠房邊墻變形體現(xiàn)出的時效分量的原因較復(fù)雜，它主要反映了開挖后應(yīng)力釋放造成的圍巖的徐變、塑性變形，結(jié)構(gòu)支護后造成的巖體間壓縮變形等。從位移變化的規(guī)律考慮，選取如下函數(shù):

k初選為4，l初選為2。

對于地下廠房開挖施工階段的統(tǒng)計模型，可選擇如下：

3.2 方差分析與逐步回歸

由于多種原因，回歸分析中初選的因子不一定對模型的建立都是重要和顯著的。如果把那些作用不顯著的因子剔除出模型后對模型影響不大，則應(yīng)該剔除不顯著因子而使模型達(dá)到優(yōu)化及方便使用。因此，統(tǒng)計模型建立中應(yīng)對各因子進行回歸效果顯著性的檢驗，原理如下：若設(shè)初選m個因子，其回歸方程為:

y=a0+a1x1+a2x2+…+am-1xm-1+amxm

經(jīng)最小二乘解算選擇其中k個因子，得統(tǒng)計模型為:

上式反映了統(tǒng)計模型中減少了一個因子xk后殘差平方和的增加量，它表明xk因子對回歸平方和貢獻的大小。

(2)

由此兩個x2變量，構(gòu)成統(tǒng)計檢驗量F：

(3)

式中n——觀測次數(shù)；

k——回歸方程因子的個數(shù)。

根據(jù)F檢驗以自由度(1,n-k)和所選的置信水平α，在F分布表中查取Fα之值。若由式(3)計算得F>Fα，則可知在置信水平α下，應(yīng)對原假設(shè)拒絕，表明xk因子對模型作用顯著，不能剔除；反之，可剔除出統(tǒng)計模型，使模型得到優(yōu)化。

若對已經(jīng)初步建立的統(tǒng)計模型各因子都按上述方法逐個地進行檢驗，那么各因子的顯著性就可以得到檢驗。把作用甚微的因子剔除而保留效果顯著的因子，使建立的最終模型達(dá)到最佳，這就是逐步回歸分析。

7號機組斷面下游374.5m高程測點是主廠房變形最大的測點，為了較好地了解各影響因素對該部位變形的影響，對該處多點位移計測點M4091-CL進行回歸分析。

取2008年4月23日～2008年11月20日時間段內(nèi)測點位移數(shù)據(jù)按上述方法取置信水平α=0.05進行逐步回歸分析，建立如下回歸方程式：

y(θ)= -16.08lnθ+4.60θ^0.5-45.60θ^(-0.5) +

0.56T+3.01H-0.08H^2+26.83

(4)

3.3 回歸分析精度與回歸效果

利用下式進行回歸分析的精度估計及剩余標(biāo)準(zhǔn)差：

回歸效果可根據(jù)復(fù)相關(guān)系數(shù)R值來判定:

式中y——直接觀測的變形值;

由未知數(shù)的個數(shù)(k+1)及自由度(n-k-1)、置信水平α，查復(fù)相關(guān)系數(shù)表Rα(k+1,n-k-1)之值。若計算得的R>Rα，則表明在置信水平α下方程的效果顯著。

如果沒有復(fù)相關(guān)系數(shù)表，也可用F分布表查取對應(yīng)值，再經(jīng)計算而求得Rα。其原理如下：由F統(tǒng)計量的定義可知:

由上式可以解得:

(5)

以置信水平α查F分布表得Fα之值，代入上式，求得的值即為Rα。

對建立的回歸方程進行精度估計，得m=0.64mm；復(fù)相關(guān)系數(shù)R=0.999，取置信水平α=0.05，用公式(5)計算Rα，得Rα=0.76，表明回歸方程在置信水平α=0.05下效果顯著。

對2008年4月23日～2008年11月20日時間段內(nèi)位移量數(shù)據(jù)進行擬合，擬合結(jié)果與實測數(shù)據(jù)之差的標(biāo)準(zhǔn)差為0.57。2008年11月20日實測位移值為45.60mm，擬合值為45.17mm。實測數(shù)據(jù)序列與擬合數(shù)據(jù)序列過程線見圖2。

3.4 測點位移的短期預(yù)測

3.4.1 預(yù)報取值范圍確定

圖2 7號機組下游邊墻374.5m高程M4091-CL測點實測值與擬合值比較

回歸分析所確定的剩余標(biāo)準(zhǔn)差m反映了所有隨機因素及方程外的有關(guān)因子對y(θ)的一次測值影響的平均變差的大小。它的單位與y(θ)相同，可以作為精度的標(biāo)志。

3.4.2 短期預(yù)報

截至2008年11月，地下廠房Ⅶ層開挖已全部結(jié)束，Ⅶ層以上支護也全部完成，準(zhǔn)備進行Ⅷ層基窩開挖。在基窩開挖之前，邊墻的結(jié)構(gòu)性態(tài)不會發(fā)生較大的變化，所以此模型對一定時間內(nèi)的短期預(yù)測仍能發(fā)揮一定的作用。 根據(jù)上述回歸方程進行短期預(yù)測，結(jié)果見表2和圖3。

表2 7號機組下游邊墻374.5m高程實測值與預(yù)測值對比

圖3 M4091-CL測點位移預(yù)測值與實測值比較

利用回歸方程對2008年11月21日～12月25日的位移進行預(yù)測，位移量為0.45mm。實測值與預(yù)測值最大差值為0.70mm，在2倍剩余標(biāo)準(zhǔn)差范圍內(nèi)，預(yù)測值有效。表明在未進行第Ⅷ層開挖的情況下，邊墻變形位移量不大，變形趨于穩(wěn)定。

3.5 各分量對效應(yīng)量影響大小分析

逐步回歸方程中的每一個因子都經(jīng)過統(tǒng)計檢驗，對回歸方程方差貢獻顯著。因此，在被考察的k-1個環(huán)境因子中，在因子間相關(guān)系數(shù)很小時，被選入方程的l個環(huán)境因子對效應(yīng)量來說都是主要因子，而未被選入方程的k-l-1個環(huán)境因子都不是主要因子。各因子變量影響的程度可用偏回歸平方和來衡量，也可用標(biāo)準(zhǔn)回歸系數(shù)或偏相關(guān)系數(shù)來反映。

按上述方法對回歸方程(4)各分量分別計算其Δi值，得出ΔH=25.45，Δθ=18.80，ΔT=4.35。由此可見，開挖進度分量對位移的影響最大，溫度因子對位移的影響最小。

4 結(jié) 論

(1)從7號機組下游邊墻374.5m高程回歸分析成果分析，在溪洛渡地下廠房施工過程中，開挖進度分量對圍巖位移影響最大，溫度因子對位移影響最小。

(2)溪洛渡水電站地下廠房施工期間建立回歸分析對其圍巖的變形進行短期預(yù)測是可行的，且其精度能基本滿足穩(wěn)定性評判的要求。

(3)施工開挖造成地下廠房結(jié)構(gòu)性態(tài)的不斷變化，需要對回歸方程不斷進行調(diào)整，使其適用于新的結(jié)構(gòu)性態(tài)。

(4)在Ⅷ層開挖尚沒有進行的情況下，邊墻變形主要受時間效應(yīng)的影響，從對7號機下游374.5m高程多點位移計的短期預(yù)測來看，其2008年11月21日～12月25日一月內(nèi)的位移量為0.45mm，變形量不大，表明圍巖變形在沒有新的施工進展的情況下已趨于穩(wěn)定。