王永慶 馮仲寧 劉意立 梁軍平

(1.中電建路橋集團(tuán)有限公司，北京 100084；2.中國水電建設(shè)集團(tuán)路橋工程有限公司梁忠高速總承包部，重慶 405200；3.清華大學(xué)環(huán)境學(xué)院，北京 100084)

?

基于地質(zhì)統(tǒng)計(jì)學(xué)的巖溶隧道工作面涌水量預(yù)測研究

王永慶1馮仲寧2劉意立3梁軍平3

(1.中電建路橋集團(tuán)有限公司，北京 100084；2.中國水電建設(shè)集團(tuán)路橋工程有限公司梁忠高速總承包部，重慶 405200；3.清華大學(xué)環(huán)境學(xué)院，北京 100084)

巖溶隧道涌水量與當(dāng)?shù)氐淖匀粭l件有著密不可分的關(guān)系，在梁忠高速公路禮讓隧道地質(zhì)環(huán)境相對穩(wěn)定的條件下，對隧道施工過程中工作面涌水逐日變化和日際變化進(jìn)行統(tǒng)計(jì)學(xué)分析。結(jié)果表明，當(dāng)變異函數(shù)采用球狀模型進(jìn)行擬合并通過普通克里格插值法可在時間域或空間域?qū)ξ磥?2天或50 m內(nèi)的工作面涌水情況做出準(zhǔn)確預(yù)測。

隧道 工作面 涌水量 預(yù)測 地質(zhì)統(tǒng)計(jì)學(xué)

目前，國內(nèi)外對隧道涌水的研究主要集中在隧道涌水量預(yù)測、地質(zhì)超前預(yù)報技術(shù)和隧道涌水防治措施這三個方面[1-3]。巖溶地區(qū)隧道涌水的研究受限于隧址區(qū)地質(zhì)蓄水結(jié)構(gòu)的多樣性、復(fù)雜性以及對含水體特性探測手段的局限性，目前仍存在諸多不足，并存在難以預(yù)測的問題[4-6]。依托梁平至忠縣高速公路禮讓隧道巖溶隧道項(xiàng)目建設(shè)開展研究，通過全過程完整的觀察和記錄隧道工作面的涌水量，基于地質(zhì)統(tǒng)計(jì)學(xué)理論對各施工段工作面涌水變化進(jìn)行預(yù)測，綜合分析隧道工作面的涌水情況，提出防治措施和施工方案。

1 隧址區(qū)地質(zhì)概況

巖溶突水影響因素可概化為地質(zhì)因素和工程因素兩方面[7]。地質(zhì)因素方面：由于巖溶地區(qū)地下水含水介質(zhì)的復(fù)雜性和異質(zhì)性，地下水分布極不均勻。地下水運(yùn)動表現(xiàn)出層流與紊流共存的特性，其中較小的裂隙構(gòu)成主要貯水空間，而較大的溶蝕裂隙作為主要導(dǎo)水通道[8]。工程因素方面：隧道開挖造成巖體擾動、卸荷，同時爆破施工以及帷幕注漿堵水等都是主要的影響因素，并在很大程度上影響了原有地下水運(yùn)動路徑，導(dǎo)致涌水的不確定性增加[8]。梁忠高速公路是梁黔高速的重要組成部分，是國家高速公路網(wǎng)的重要補(bǔ)充，起于梁平碧山鎮(zhèn)，終點(diǎn)位于忠縣拔山鎮(zhèn)，線路全長約72 km，禮讓隧道總長5 517.6 m，位于四川盆地東部暖濕性亞熱帶季風(fēng)氣候區(qū)，氣候溫暖，雨量較多。隧址區(qū)屬低山丘陵地貌，各含水巖組的構(gòu)造裂隙發(fā)育，為典型的巖溶地貌。

經(jīng)勘查，該地區(qū)地下水可分為松散巖間孔隙水、碎屑巖孔隙間裂隙水、風(fēng)化后砂泥巖間裂隙水和碳酸鹽巖間巖溶水(巖溶裂隙、溶洞水)四類，其中松散巖間孔隙水分布呈零星特征，水量占比較小，對隧道涌水貢獻(xiàn)較弱，而廣泛分布的碎屑巖孔隙間裂隙水和碳酸鹽巖間巖溶水構(gòu)成隧道涌水的主要來源，可能對施工造成較大影響。

2 地質(zhì)統(tǒng)計(jì)學(xué)計(jì)算方法

上世紀(jì)60年代，G.Matheron創(chuàng)立了地質(zhì)統(tǒng)計(jì)學(xué)(以下簡稱“地統(tǒng)學(xué)”)是以區(qū)域變量為理論基礎(chǔ)，采用變異函數(shù)為研究、計(jì)算工具，對各類在時空域具有隨機(jī)性、結(jié)構(gòu)性的自然現(xiàn)象進(jìn)行統(tǒng)計(jì)分析并預(yù)測推算的學(xué)科。1977年桂林冶金地質(zhì)研究所情報室首次將地統(tǒng)學(xué)引入我國，之后孫洪泉、高歌等拓展了其在礦業(yè)、氣象方面的應(yīng)用[9，10]。同時，侯景儒、余先川等將地統(tǒng)學(xué)的理論、方法延伸到了時空域中[11，12]。

經(jīng)過多年的發(fā)展，目前在地統(tǒng)學(xué)上已形成比較認(rèn)同的共識，即凡是要研究時空分布隨機(jī)數(shù)據(jù)的隨機(jī)性、結(jié)構(gòu)性，并進(jìn)行無偏估計(jì)，或要模擬研究對象的離散性、波動性或者其他性質(zhì)多采用地統(tǒng)學(xué)理論與方法。當(dāng)某一變量表現(xiàn)在空間維度上連續(xù)分布時，可稱為“區(qū)域化變量”，這種變量可用于反映某種現(xiàn)象的空間特征，采用“區(qū)域化變量”來描述的現(xiàn)象可被稱為“區(qū)域化現(xiàn)象”。地統(tǒng)學(xué)中所指的區(qū)域化變量是指在一定的相鄰研究域內(nèi)，以時間或空間坐標(biāo)為自變量的隨機(jī)場，具有時空局限性、連續(xù)性、異向性，具有變程、套和結(jié)構(gòu)等屬性。

變異函數(shù)可用于描述區(qū)域化變量的隨機(jī)性和結(jié)構(gòu)性，并從數(shù)學(xué)理論方面給予嚴(yán)格分析，因此該函數(shù)可作為分析時空變異規(guī)律分析和時空結(jié)構(gòu)分析的有效工具。當(dāng)研究對象處于一維空間時，變異函數(shù)可被定義為“區(qū)域化變量Z(x)分別在點(diǎn)x和x+h處數(shù)值差方差的一半”，并可記為該區(qū)域化變量Z(x)在x軸方向上的變異函數(shù)(γ(h))，用數(shù)學(xué)的方法表述即為

以h為橫坐標(biāo)，以γ(h)為縱坐標(biāo)，畫出變異函數(shù)曲線圖，從而標(biāo)示出區(qū)域化變量Z(x)的空間變異特點(diǎn)。同時可直接從變異函數(shù)圖中得到反映自然現(xiàn)象空間分布結(jié)構(gòu)或相關(guān)范圍。

目前，按擬合方式的不同變異函數(shù)理論模型可分為三類：有基臺值模型、無基臺值模型和孔穴效應(yīng)模型。當(dāng)前，采用有基臺值的球狀模型在地統(tǒng)學(xué)研究中成為主流，該模型的數(shù)學(xué)表達(dá)為

式中：C0為塊金值，用于反應(yīng)表示區(qū)域化變量的非連續(xù)變異特性，該值由區(qū)域化變量的屬性或測量誤差決定。

C為拱高，C0+C為基臺值，當(dāng)探測點(diǎn)間距h逐漸增大時，變異函數(shù)γ(h)的值將從初始的塊金值逐漸增大，并隨增速下降，最終穩(wěn)定在一個常數(shù)，該常數(shù)值被稱為基臺值。實(shí)際操作中，當(dāng)探測點(diǎn)間距h拉大，使變異函數(shù)數(shù)值超過基臺值時，即變異函數(shù)數(shù)值不再隨探測點(diǎn)間距離的增加而改變，表明這兩點(diǎn)空間不存在相關(guān)性。

a為變程，即變異函數(shù)由塊金值增加到基臺值時的距離。當(dāng)探測點(diǎn)間距h≥a后，區(qū)域化變量Z(x)與Z(x+h)的空間相關(guān)性消失。

而當(dāng)0

則可以得到線性模型

y=b0+b1x1+b2x2

根據(jù)已有數(shù)據(jù)，對上式進(jìn)行最小二乘擬合，即可得到變異函數(shù)參數(shù)C0、C和a。隨后通過普通克里格插值法即可估計(jì)鄰域內(nèi)塊段或區(qū)間平均值的無偏線性估計(jì)量[13]。

3 基于地質(zhì)統(tǒng)計(jì)學(xué)的巖溶隧道工作面涌水預(yù)測

按照地統(tǒng)學(xué)基本理論，對梁平至忠縣高速公路禮讓隧道巖溶隧道工作面涌水進(jìn)行預(yù)測。

3.1 選定合理段落

根據(jù)隧道設(shè)計(jì)資料，ZK16+150～ZK16+300段隧道穿越地段為三疊系中統(tǒng)雷口坡組底層，為頁巖、薄層灰?guī)r、泥灰?guī)r互層。頁巖為塊狀結(jié)構(gòu)，頁理構(gòu)造發(fā)育，灰?guī)r、泥灰?guī)r為薄層狀。該段巖層軟硬不均，差異較大。巖體中構(gòu)造發(fā)育，巖體以完整-破碎為主。隧道出水形式以潮濕或點(diǎn)滴出水為主。據(jù)檢測，圍巖Rc=12.43MPa，完整系數(shù)Kv=0.62，K1=0.10，K2=0.20，K3=0.00，[BQ]=252.29。該段在整個隧道中具有代表性，故選用該段作為研究對象。

3.2 試驗(yàn)設(shè)計(jì)

由于事先資料缺乏，實(shí)驗(yàn)數(shù)據(jù)分為兩個層次，(1)按照每日工作面為空間變化量測出水量，從而形成時間空間、距離空間上的涌水量。(2)個別工作面測量每日出水量。其試驗(yàn)數(shù)據(jù)見圖1。

由于工作面施工直接將巖體剝離，會造成原本地下水系破壞，因此通過研究工作面涌水情況能夠更真實(shí)地反映地下水體所受到的影響。根據(jù)2014年4月16日至2014年6月9日共計(jì)63d對該隧道施工工作面涌水量的記錄(圖1)，在前38d內(nèi)分別于每天早中晚三個時段分別記錄工作面的涌水情況，結(jié)果如圖2所示。

圖1 涌水量隨位置與時間的變化關(guān)系

圖2 工作面涌水量日際變化

隧道施工工作面涌水量變化顯著，最高時可達(dá)到60 m3/h，最低時僅有8.5 m3/h，相差約7倍。因此，對于不同施工段防突涌水的要求是不同的。施工過程中工作面的平均涌水量為23 m3/h，表明隧道施工對于當(dāng)?shù)氐叵滤h(huán)境的影響較大，工作面每日早、中、晚涌水情況變化不大，圖2還表明每日涌水量較日間變化更為顯著，說明各水平向地下水體聯(lián)系較為松散，以獨(dú)立水文單元構(gòu)成為主。

3.3 時間域空間、空間域空間變異函數(shù)計(jì)算

根據(jù)統(tǒng)計(jì)的63 d內(nèi)工作面日均涌水量，按時間域和空間域分別采用變異函數(shù)理論球狀模型進(jìn)行擬合，并得到如下擬合結(jié)果(如圖3)。

時間域：塊金值C0=3.948 1，拱高C=5.703 0，變程a=12.053 8

空間域：塊金值C0=3.712 6，拱高C=653.911 4，變程a=49.088 7

圖3 時間域和空間域變異函數(shù)

由于C0/C+C0可以反映不同變量由隨機(jī)變量引起的變異，當(dāng)比值接近于1時表明隨機(jī)性變異占主導(dǎo)，而當(dāng)比值接近于0時，表明空間/時間相關(guān)性占主導(dǎo)，具體為

C0/C+C0<25% 空間/時間相關(guān)性強(qiáng)；

25%

75%

對于隧道施工工作面涌水分析，在時間域該比值為40.9%，在空間域該比值為0.56%，表明當(dāng)采用空間變量進(jìn)行預(yù)測時該函數(shù)模型具有極好的準(zhǔn)確性，單次可估算的范圍至多為49 m，而對時間域的分析表明，工作面涌水量同時間變量具有中等強(qiáng)度的相關(guān)性，單次可估算的范圍至多為12 d。與實(shí)際施工相聯(lián)系，認(rèn)為工作面涌水量更易受到相鄰地質(zhì)段的影響，而隧址區(qū)巖石裂隙發(fā)育，涌水量直接受降雨影響程度較大，因此表現(xiàn)出較強(qiáng)的隨機(jī)性特征。

3.4 模型準(zhǔn)確性評估

根據(jù)擬合出的球狀模型，采用普通克里格法進(jìn)行空間插值，并選取ZK16+170、ZK16+266兩個工作面進(jìn)行驗(yàn)證，其結(jié)果如表1。

表1 實(shí)測值與計(jì)算值對比

由表1可以看出，采用受隨機(jī)變量影響較大的時間域作為研究對象，在變程a的范圍內(nèi)，(小于12 d)計(jì)算值與實(shí)測值偏差較小，均方根誤差(RMSE)分析表明，該方法具有較高的準(zhǔn)確性。

4 結(jié)論與建議

涌水量大小對巖溶隧道施工過程具有極大影響，同時與隧址區(qū)自然條件有著密不可分的關(guān)聯(lián)。分別采用時間域和空間域作為變量建立地統(tǒng)學(xué)變異函數(shù)并通過球狀模型進(jìn)行擬合，發(fā)現(xiàn)涌水量受相鄰地質(zhì)段和隨機(jī)降雨過程的影響較大。利用普通克里格插值法可對一定時間域或空間域內(nèi)工作面涌水量進(jìn)行預(yù)測，通過選取兩個典型工作面進(jìn)行模型檢驗(yàn)表明：當(dāng)所預(yù)測點(diǎn)落于變程內(nèi)時，該模型具有較高的準(zhǔn)確性。

[1] 劉勇.北武夷山隧道嵐谷斜井涌水預(yù)報技術(shù)[J].鐵道勘察,2013,39(2):34-36

[2] 李術(shù)才,劉斌,孫懷鳳,等.隧道施工超前地質(zhì)預(yù)報研究現(xiàn)狀及發(fā)展趨勢[J].巖石力學(xué)與工程學(xué)報,2014,33(6):1090-1113

[3] 王遇國.巖溶隧道突水災(zāi)害與防治研究[D].北京：中國鐵道科學(xué)研究院,2010

[4] 李光偉.論巖溶區(qū)工程地質(zhì)勘察問題與地質(zhì)選線[J].鐵道勘察,2016(2):12-15，32

[5] 李建偉.某鐵路隧道水文地質(zhì)分析及涌水量預(yù)測[J].鐵道勘察,2011,37(6):72-76

[6] 喻成云.基于統(tǒng)計(jì)分析的西南巖溶隧道涌水量預(yù)測探析[D].成都：成都理工大學(xué),2013

[7] 楊艷娜.西南山區(qū)巖溶隧道涌突水災(zāi)害危險性評價系統(tǒng)研究[D].成都：成都理工大學(xué),2009

[8] 王國斌.滬蓉西高速公路烏池壩巖溶隧道涌水成災(zāi)機(jī)理研究[D].北京：中國地質(zhì)大學(xué),2012

[9] 孫洪泉.地質(zhì)統(tǒng)計(jì)學(xué)及其應(yīng)用[M].徐州：中國礦業(yè)大學(xué)出版社,1990

[10]高歌, 龔樂冰, 趙珊珊,等.日降水量空間插值方法研究[J].應(yīng)用氣象學(xué)報,2007,18(5):732-736

[11]侯景儒,王志民.時間-空間域中多元信息的地質(zhì)統(tǒng)計(jì)學(xué)[J].北京科技大學(xué)學(xué)報,1995(2):101-106

[12]余先川,侯景儒,程曉春,等.時空信息統(tǒng)計(jì)學(xué)的一些基本理論和方法[J].北京師范大學(xué)學(xué)報:自然科學(xué)版,2003,39(3):353-359

[13]楊勇,梅楊,張楚天,等.基于時空克里格的土壤重金屬時空建模與預(yù)測[J].農(nóng)業(yè)工程學(xué)報,2014(21):249-255

Study on the Predication of Water Inrush in Tunnel Face Based on Geostatistic

WANG Yongqing1FENG Zhongning2LIU Yili3LANG Junping3

2016-05-20

王永慶(1978—)，男，1999年畢業(yè)于黃河水利學(xué)院水利工程專業(yè)，工程師。

1672-7479(2016)05-0075-03

P332； P628+.1

B