呂晶日，張 強(qiáng)，韓曉磊，張志強(qiáng)，竹 娜，胡元進(jìn)

(成都理工大學(xué)，四川 成都610059)

近年來(lái)我國(guó)經(jīng)濟(jì)建設(shè)迅速發(fā)展，基礎(chǔ)設(shè)施建設(shè)也在日益加快。隧洞作為基礎(chǔ)設(shè)施建設(shè)的重要組成部分大大縮短了施工建設(shè)的里程，但同時(shí)也帶來(lái)了不少環(huán)境問(wèn)題。隧洞施工排水改變了地下水的徑流和排泄條件，導(dǎo)致地下水位降低，甚至井泉干涸。在干旱少雨的地區(qū)或季節(jié)影響當(dāng)?shù)鼐用竦纳a(chǎn)生活用水。因此，預(yù)測(cè)隧洞施工工程中地下水位變化規(guī)律，為工程設(shè)計(jì)和建設(shè)提供依據(jù)就顯得意義重大。本文利用Visual Modflowd地下水三維數(shù)值模擬軟件對(duì)海東隧洞施工過(guò)程中地下水滲流場(chǎng)的動(dòng)態(tài)變化進(jìn)行了模擬預(yù)測(cè)。

1 海東隧洞萂村段三維數(shù)值模擬

1.1 隧址區(qū)概況

賓川縣大營(yíng)鎮(zhèn)萂村地處大理白族自治州賓川縣大營(yíng)鎮(zhèn)西北邊，距縣城28 km，轄東莊、黑家營(yíng)、茅草坪等21個(gè)村民小組。全村國(guó)土面積81.99 km2，海拔1 860 m，年平均氣溫20.4℃，年降水量600 mm，現(xiàn)有人口5 784人。全村供水來(lái)源于井水、泉水和地表水。

海東隧道萂村段位于大理白族自治州賓川縣大營(yíng)鎮(zhèn)境內(nèi)，隧址區(qū)出露多種地層，包括第四系全新統(tǒng)(Qh)、二迭系下統(tǒng)(P1)、玄武巖系(Pβ1)、石炭系中上統(tǒng)(C2+3)、下統(tǒng)(C1)、中統(tǒng)長(zhǎng)育村組(D2ch)、下統(tǒng)蓮花曲組(D1l)、青山組(D1q)、康廊組(D1k)(見(jiàn)圖1)。

圖1 海東隧洞萂村段模型范圍

靠近黑家營(yíng)馬頭山至大黑山地表巖溶發(fā)育，溶蝕現(xiàn)象強(qiáng)烈，發(fā)育大型巖溶洼地，降雨后，地表水通過(guò)洼地的匯聚，再由落水洞對(duì)地下水巖溶管道進(jìn)行落入式補(bǔ)給，在斷裂帶及地層分界附近以串珠狀泉點(diǎn)形式出露地表。

表1 研究區(qū)監(jiān)測(cè)點(diǎn)情況

泉水和地表水是隧址區(qū)內(nèi)村民的重要生產(chǎn)生活用水。而隧洞施工影響到地下水滲流場(chǎng)，地下水位降低，形成地下水降落漏斗，導(dǎo)致局部區(qū)域地表水疏干等，會(huì)給部分村民用水帶來(lái)困難，所以預(yù)測(cè)隧洞施工工程中，地下水位的動(dòng)態(tài)變化預(yù)測(cè)就顯得非常重要。

1.2 三維地下水?dāng)?shù)值模型建立

通過(guò)Visual Modflow對(duì)隧址區(qū)地質(zhì)模型的模擬，實(shí)現(xiàn)海東隧道萂村段概念模型的概化。建立一個(gè)既符合巖溶水滲流特點(diǎn)又符合實(shí)際水文地質(zhì)條件的滲流模型模擬，實(shí)現(xiàn)對(duì)天然狀態(tài)下的地下水水位的擬合，進(jìn)而對(duì)隧道開(kāi)挖進(jìn)程不同時(shí)期的地下水滲流場(chǎng)的變化進(jìn)行預(yù)測(cè)。

1.2.1 計(jì)算單元與邊界條件概化

根據(jù)研究區(qū)域大小及計(jì)算精度要求，采用矩形單元對(duì)模型進(jìn)行剖分。模擬隧址區(qū)X方向總長(zhǎng)度為7 900 m，Y方向總長(zhǎng)度8 900 m，按照100×100 m網(wǎng)格將其剖分成89列，79行，層單元格數(shù)目為7 031個(gè);垂向Z方向取760 m，劃分為16層，模型總共112 496個(gè)單元格;模擬選取模型底板為1 750 m高程，模型頂板高程由Sufer軟件插值萂村地形等高線(xiàn)生成GRD文件導(dǎo)入。

模型邊界:海東隧道萂村段地下水補(bǔ)給主要來(lái)源于大氣降水，研究區(qū)內(nèi)隧道北東方向的溝底河流為排泄基準(zhǔn)面;隧址區(qū)地表水系發(fā)育，萂頭至萂尾三個(gè)地表定水頭邊界，具體的模型邊界概化示意圖見(jiàn)圖2。

圖2 模型邊界概化示意圖

1.2.2 模型的參數(shù)取值與模擬校驗(yàn)

模型參數(shù)的選取和確定，是模型建立過(guò)程中的一個(gè)重要步驟，主要包括滲透系數(shù)、巖土層給水度、貯水率、有效孔隙度和總孔隙度等指標(biāo)。

根據(jù)已有的鉆孔抽水、壓水試驗(yàn)，經(jīng)過(guò)統(tǒng)計(jì)分析、篩選及整理。結(jié)合以往過(guò)程計(jì)算經(jīng)驗(yàn)，在此范圍內(nèi)反復(fù)試算，初步確定模型參數(shù)，對(duì)模型進(jìn)行多次反演、校驗(yàn)最終確定模型參數(shù)如表2所示。

表2 模型水文地質(zhì)參數(shù)選取表

從表3可以看出模擬計(jì)算的觀測(cè)水位與實(shí)測(cè)水位存在一定的誤差，但誤差均在一個(gè)合理的比例范圍內(nèi)，差值比例范圍在0.09%～1.1%。創(chuàng)建物理模型成功，模擬結(jié)果所得滲流場(chǎng)可作為研究區(qū)的天然地下水滲流場(chǎng)，如圖3所示，為隧址區(qū)地下水三維滲流場(chǎng)。從圖3中可以看出:上新村、栽秧菁、冷水箐、萂頭村、萂中村、萂尾村均處于高山峽谷盆地，地下水位比較低;此外，萂村的三個(gè)水庫(kù)地下水揭示顏色都為藍(lán)色，地下水位最低。

模型中6個(gè)監(jiān)測(cè)點(diǎn)也都出露在地下水位較低的地方。其中，萂頭村的 DL-S04、萂尾村的 DL-S05和黑家營(yíng) DLS06的水位最低。

本次模型的模擬運(yùn)行運(yùn)用了非穩(wěn)定流對(duì)地下水的初始水頭進(jìn)行擬合，模擬時(shí)間為十年。選取6個(gè)泉點(diǎn)，分別為DL-S01、DL-S02、DL-S03、DL-S04、DL-S05、DL-S06，目的在于初始水頭的校驗(yàn)，以確定擬合準(zhǔn)確度，通過(guò)模型計(jì)算得出更貼近實(shí)際水位，來(lái)完成模型天然狀態(tài)滲流場(chǎng)的模擬。

表3 模型計(jì)算泉點(diǎn)水位與實(shí)際水位擬合情況

圖3 海東隧道萂村段天然地下水滲流場(chǎng)

2 施工條件下萂村段地下水動(dòng)態(tài)變化

運(yùn)行模型后得到萂村段天然地下水滲流場(chǎng)(如圖4所示)。從圖中可以看出:從圖4中可以看出局部地區(qū)滲流場(chǎng)發(fā)生變化，藍(lán)色區(qū)域明顯擴(kuò)大，該區(qū)域的監(jiān)測(cè)點(diǎn)DL-S01、DL-S04和DL-S05的地下水位都有明顯降深;監(jiān)測(cè)點(diǎn)DLS01、DL-S02和DL-S06此時(shí)在三維地下水滲流場(chǎng)中還沒(méi)有發(fā)生明顯變化。

運(yùn)行Visual MODFLOW模型后，將6個(gè)監(jiān)測(cè)點(diǎn)的每一年地下水位做一次記錄，得到6個(gè)監(jiān)測(cè)點(diǎn)十年內(nèi)地下水位的動(dòng)態(tài)變化并整理如表4所示。

圖4 海東隧道萂村段地下水滲流場(chǎng)

圖5 地下水降落漏斗影響范圍

表4說(shuō)明了:6個(gè)觀測(cè)點(diǎn)的地下水位均發(fā)生了明顯的降深，并呈現(xiàn)地下水位逐年遞減的趨勢(shì)。其中，監(jiān)測(cè)點(diǎn)DLS01和DL-S04的降深最為明顯，分別降低了51.52 m和59.56 m。此外監(jiān)測(cè)點(diǎn) DL-S02、DL-S03、DL-S05、DLS06的地下水位也都有不同程度的降深，分別為40.25 m、36.14 m、50 m、41.8 m。

表4 萂村段地下水動(dòng)態(tài)變化

隨著隧洞施工的進(jìn)行，降落漏斗的影響范圍也逐漸擴(kuò)大。從圖5可得知，隧洞引水工程施工后，沿引水線(xiàn)路3處出現(xiàn)降落漏斗，分別位于上新村、萂頭箐、黑家營(yíng)一帶。地下水位有明顯降深，上新村受地下水形成影響面積為4.45 km的降落漏斗，降深也最為明顯，最大中心水位降達(dá)240 km;萂頭箐地下水形成影響面積為7.52 km的橢圓形降落漏斗，最大中心水位降達(dá)210 m;黑家營(yíng)地下水形成影響面積為0.748 km的圓形降落漏斗，較其他兩個(gè)降落漏斗的影響范圍小，最大中心水位降達(dá)180 m。而隧址區(qū)的6個(gè)監(jiān)測(cè)點(diǎn)中，編號(hào)為DL-S01和DL-S04的監(jiān)測(cè)點(diǎn)均在降落漏斗的影響范圍內(nèi)，而編號(hào)為 DL-S02、DL-S05、DLS06的監(jiān)測(cè)點(diǎn)距離降落漏斗的最大距離也不過(guò)266.43 m。上新村、萂頭箐、黑家營(yíng)地表水體滲流場(chǎng)已經(jīng)改變，可能導(dǎo)致地表水位下降，部分村落泉點(diǎn)局部疏干，影響到居民的生產(chǎn)生活用水。

3 結(jié)語(yǔ)

利用地下水?dāng)?shù)值計(jì)算模型建立巖溶地區(qū)萂村段三維地下水?dāng)?shù)值模型，模擬隧址區(qū)地下水水位變化，確定水文地質(zhì)參數(shù)，對(duì)模擬結(jié)果進(jìn)行預(yù)測(cè)分析。分析表明:隨著隧洞施工的推進(jìn)，隧址區(qū)滲流場(chǎng)逐漸受到干擾，水位逐漸降低。沿著隧洞線(xiàn)形成地下水降落漏斗，漏斗范圍逐年擴(kuò)大，并且影響范圍已至上新村、栽秧箐、萂頭箐和黑家營(yíng)村。其最大影響范圍達(dá)到7.52 km。地下水位的降低，地下水降落漏斗的形成，將導(dǎo)致部分村落的井和地表水干涸，部分泉點(diǎn)衰竭、枯竭或斷流，從而影響到當(dāng)?shù)鼐用竦纳a(chǎn)生活用水。隧洞開(kāi)挖時(shí)采取合理的恢復(fù)水源措施，以保證當(dāng)?shù)鼐用竦纳a(chǎn)生活用水是非常必要的。

[1]陳喜，劉傳杰，胡忠明，等.泉域地下水?dāng)?shù)值模擬及泉流量動(dòng)態(tài)變化預(yù)測(cè)[J].水文地質(zhì)工程地質(zhì).2006(2).

[2]張繼勛，任旭華，嚴(yán)冬青，等.巖溶富水區(qū)地下工程施工期地下水位變化規(guī)律分析[J].水文水資源.2013 39(12).

[3]孫晉玉，張強(qiáng)，許模，等.隧道施工涌水后地下水特征分析[J].工程地質(zhì)學(xué)報(bào).2002(10).

[4]白細(xì)民，賀小榮，彭衛(wèi)星.南昌市地下水位降落漏斗發(fā)展趨勢(shì)R/S分析[J].南水北調(diào)與水利科技.2013 11(01).

[5]陳亞峰.施工條件下海東隧道巖溶水文地質(zhì)條件演化預(yù)測(cè)[D].成都理工大學(xué)碩士學(xué)位論文.2013.