李銀泉 陳長(zhǎng)生 張海平

摘 要 ?工程隧洞施工不僅可能存在因突涌水導(dǎo)致的安全隱患問(wèn)題，而且還會(huì)在一定程度上破壞工程區(qū)地下水滲流場(chǎng)平衡，影響地下水的正常循環(huán)，存在影響或疏干地下水而導(dǎo)致的環(huán)境影響問(wèn)題，處置不當(dāng)還可能造成嚴(yán)重的環(huán)境影響及社會(huì)穩(wěn)定問(wèn)題。本文以滇中引水工程某深埋長(zhǎng)隧洞為依托，采用GMS數(shù)值模型對(duì)不同施工工況下隧洞涌水量、地下水水位變化、主要泉點(diǎn)影響程度、運(yùn)行期水位恢復(fù)情況等進(jìn)行預(yù)測(cè)研究，為合理評(píng)價(jià)隧洞工程建設(shè)可能導(dǎo)致的地下水環(huán)境影響程度、研究制定滲控標(biāo)準(zhǔn)等提供合理、準(zhǔn)確、可行的依據(jù)。

關(guān)鍵詞 ?地下水環(huán)境影響;GMS;隧洞施工;預(yù)測(cè) ;滇中引水工程

中圖分類(lèi)號(hào)：X523 ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

Abstract： Tunnel construction is not only the safe potential risk that may be caused by sudden flood water problems，but also to a certain extent，destroy the existing groundwater seepage field balance，affect the normal cycle of the groundwater，influence or drainage of groundwater environmental impact caused by the problem，improper disposal can also cause serious problems of environmental impact and social stability. This paper，taking a deep water diversion project in Yunnan long tunnel as the backing，the GMS numerical model is adopted to different construction conditions influence radius of tunnel water inflow，water level changes of main stream points influence degree，operation period were used to predict the water level recovery，etc，for reasonable evaluation of tunnel engineering construction may lead to groundwater seepage control standards in the environmental impact，the research provides the basis of a reasonable， precise and feasible.

Keywords： groundwater environmental impact; GMS; tunnel construction; prediction; water diversion project in Yunnan

隨著國(guó)家實(shí)施西部大開(kāi)發(fā)戰(zhàn)略進(jìn)程的加快，西南地區(qū)的水利水電工程建設(shè)項(xiàng)目越來(lái)越多。其中水利工程中為了引調(diào)水，常以水工隧洞的型式穿越山體，水工隧洞總體呈現(xiàn)出“長(zhǎng)、大、深”的特點(diǎn)，代表性的有引大入秦工程、引黃入晉工程、新疆ABH隧洞等引調(diào)水工程[1-3]。工程隧洞施工，特別是巖溶富水隧洞施工極易引起周?chē)鷧^(qū)域的地下水環(huán)境改變。隧洞開(kāi)挖形成地表水及地下水匯流新的排泄通道，不但造成施工涌突水隱患，而且擾動(dòng)天然地表水及地下水環(huán)境[4] [5]。因此，國(guó)內(nèi)外專(zhuān)家學(xué)者一直在對(duì)工程建設(shè)引發(fā)周?chē)叵滤h(huán)境影響進(jìn)行預(yù)測(cè)研究。一方面從宏觀上以現(xiàn)狀和長(zhǎng)期預(yù)測(cè)的思路橫向指出隧洞施工可能引發(fā)的地下水環(huán)境改變問(wèn)題，另一方面以具體隧洞實(shí)例進(jìn)行分析，縱向系統(tǒng)的研究其引發(fā)的地下水環(huán)境影響[6]。地下水環(huán)境影響預(yù)測(cè)方法國(guó)內(nèi)外專(zhuān)家學(xué)者也提出了一系列的評(píng)價(jià)和研究方法，主要包括數(shù)學(xué)模型法和類(lèi)比預(yù)測(cè)法。數(shù)學(xué)模型法中的數(shù)值法是目前地下水環(huán)境預(yù)測(cè)主要的發(fā)展趨勢(shì)。數(shù)值法以其直觀、全方位、精度高、靈活性高的特點(diǎn)被廣泛運(yùn)用于地下水相關(guān)領(lǐng)域。

本文以滇中引水工程某深埋長(zhǎng)隧洞為依托，以GMS為平臺(tái)，結(jié)合其地下水含水介質(zhì)特性、巖體水文地質(zhì)結(jié)構(gòu)、地下水流動(dòng)系統(tǒng)特征、地下水化學(xué)及其環(huán)境特點(diǎn)建立三維滲流場(chǎng)數(shù)值模型，對(duì)隧洞涌水量及對(duì)地下水環(huán)境影響程度進(jìn)行預(yù)測(cè)研究，為合理評(píng)價(jià)隧洞工程建設(shè)可能導(dǎo)致的地下水環(huán)境影響程度、研究制定滲控標(biāo)準(zhǔn)等提供合理、準(zhǔn)確、可行的依據(jù)。

1 ?工程及地質(zhì)背景條件概況

滇中引水工程是云南省經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的戰(zhàn)略性基礎(chǔ)工程，擬從金沙江干流石鼓河段取水，工程多年平均引水量34.03億m3。本文所涉及的深埋隧洞全長(zhǎng)約62km，最大埋深約1450m，埋深大于600m洞段占隧洞總長(zhǎng)的67.38%，屬于典型的深埋長(zhǎng)隧洞。

該隧洞地處橫斷山脈與滇中高原交接部位，跨越金沙江與瀾滄江分水嶺，地形起伏較大，隧洞區(qū)構(gòu)造背景及巖溶水文地質(zhì)條件均極其復(fù)雜，穿越多條活動(dòng)斷裂、多段巖溶地層和多個(gè)巖溶水系統(tǒng)、以及泥頁(yè)巖等軟弱地層。涉及的工程地質(zhì)、水文地質(zhì)及環(huán)境地質(zhì)問(wèn)題眾多。其中線路地表涉及可溶巖地層分布長(zhǎng)度29.703km，占比47.45%，隧洞穿越可溶巖地層長(zhǎng)度17.866km，占比28.54%，沿線可溶巖地層分布有一系列巖溶泉，這些泉水是當(dāng)?shù)厝嗣裆詈凸まr(nóng)業(yè)生產(chǎn)、以及生態(tài)補(bǔ)水的主要補(bǔ)給水源。

2 ?模型的建立

2.1 水文地質(zhì)概念模型

水文地質(zhì)模型評(píng)價(jià)區(qū)面積約1575km2，包含了5個(gè)巖溶系統(tǒng)：白漢場(chǎng)、拉什海、文筆海、清水江-劍川、鶴慶西山巖溶系統(tǒng)。地下水主要賦存于灰?guī)r和局部的第四系中，導(dǎo)水通道主要為導(dǎo)水性斷層和巖溶管道。淺層水的補(bǔ)給來(lái)源主要為降水，開(kāi)采、蒸發(fā)、泉為主要排泄途徑。深層水補(bǔ)給來(lái)源主要來(lái)自上層下滲、區(qū)域地下水的循環(huán)，開(kāi)采、泉和側(cè)向排泄為主要排泄途徑。將評(píng)價(jià)區(qū)地下水水流系統(tǒng)概化為非均值各項(xiàng)異性三維流動(dòng)系統(tǒng)[7]。根據(jù)區(qū)內(nèi)發(fā)育的主要泉水點(diǎn)、河流及主要斷層，對(duì)研究區(qū)進(jìn)行了三維滲流模型概化，邊界區(qū)域與邊界類(lèi)型設(shè)置如圖1。

2.2 數(shù)學(xué)模型

根據(jù)水文地質(zhì)概念模型，研究區(qū)巖溶含水層含水介質(zhì)以裂隙為主，地下水流符合達(dá)西定律。將研究區(qū)地下水概化為非均質(zhì)、各向異性的穩(wěn)定三維地下水流系統(tǒng)?？捎萌缦挛⒎址匠痰亩ń鈫?wèn)題來(lái)描述：

式中：H為地下水水頭，m;Kx，Ky，Kz，為x，y，z方向滲透系數(shù)，m/d;W為單位時(shí)間從單位體積含水層流入或流出的水量，1/d;ε為單位時(shí)間在垂向上從單位面積含水層自由面流入或流出的水量，m/ d;H1為各層邊界水位，m;q為含水層第二類(lèi)邊界單位面積過(guò)水?dāng)嗝嫜a(bǔ)給流量，m/d;Ω為滲流區(qū)域;B1為水頭已知邊界，第一類(lèi)邊界;B2為流量已知邊界，第二類(lèi)邊界;B3為自由面邊界;n為滲流區(qū)邊界的單位外發(fā)線方向。地下水流場(chǎng)的初值設(shè)置為地表標(biāo)高。

2.3 數(shù)值模型

此基礎(chǔ)上，地下水?dāng)?shù)值模型采用GMS10.2軟件進(jìn)行數(shù)值離散（圖2）。巖體滲透系數(shù)根據(jù)壓水試驗(yàn)和振蕩試驗(yàn)結(jié)果進(jìn)行取值。對(duì)于非飽和帶巖體的滲透系數(shù)以及鉆孔未揭示的巖體的滲透系數(shù)，需要采用類(lèi)比和反演的方法來(lái)確定。利用隧洞區(qū)水位長(zhǎng)觀孔的長(zhǎng)期觀測(cè)資料、主要泉流量值和已經(jīng)開(kāi)工建設(shè)的多條施工支洞涌水量監(jiān)測(cè)值等資料，對(duì)模型及參數(shù)進(jìn)行調(diào)整與修改。通過(guò)擬合發(fā)現(xiàn)，總體上觀測(cè)值與計(jì)算值相差較小[8] [9]。

3 ?隧洞灌漿防滲不同工況擬定

模型擬結(jié)合以下三種工況預(yù)測(cè)不同工況下的隧洞涌水量、地下水水位變化和主要泉點(diǎn)影響程度等。

工況一：天然裸洞;

工況二：部分灌漿（對(duì)裸洞涌水量q>3m3/m*d 洞段實(shí)施普通灌漿，灌漿標(biāo)準(zhǔn)1×10-5cm/s）;

工況三：二次灌漿（對(duì)普通灌漿后涌水量q>3～5m3/m*d洞段修改為磨細(xì)水泥灌漿，灌漿標(biāo)準(zhǔn)5×10-6cm/s;對(duì)普通灌漿后涌水量q>5m3/m*d洞段修改為化學(xué)灌漿，灌漿標(biāo)準(zhǔn)1×10-6cm/s）。

4 ?隧洞施工對(duì)地下水環(huán)境影響預(yù)測(cè)

4.1 隧洞涌水量預(yù)測(cè)

隧洞施工貫通后全隧洞排水的情況，根據(jù)數(shù)值模擬計(jì)算結(jié)果，統(tǒng)計(jì)了隧洞各巖性分段在不同工況條件下的涌水量，隧洞在各工況條件下的總涌水量見(jiàn)表1。

裸洞工況下隧洞涌水量較大，總涌水量達(dá)到4.496m3/s（枯水期）、4.962m3/s（豐水期）。涌水量較大的洞段為巖溶水系統(tǒng)段、區(qū)域性斷層影響帶以及向斜（背斜）影響段;在二次灌漿后，涌水量明顯減小，總涌水量減小為1.398m3/s（枯水期）、1.426m3/s（豐水期），降幅達(dá)到60%以上。

為進(jìn)一步分析隧洞開(kāi)挖對(duì)隧洞區(qū)地下水總徑流量的影響，分別統(tǒng)計(jì)了隧洞開(kāi)挖涌水量與地下水總徑流量占比（表2）。

在天然裸洞施工工況下，地下水徑流量變幅占比約18%。采用防滲灌漿措施后，徑流量變幅變小。二次灌漿（工況三） 后，徑流量變幅占比降到約5%。

4.2 地下水水位變化預(yù)測(cè)

為定量計(jì)算隧洞排水疏干可能引起的地下水位的變化，通過(guò)所建的地下水滲流模型，結(jié)合MapGis的高程疊加功能，分別計(jì)算了不同工況下的水位下降和影響范圍?？菟诤拓S水期不同工況下隧洞開(kāi)挖地下水水位下降影響可見(jiàn)表3和圖3、4天然狀態(tài)（裸洞）下隧洞地下水位變化及影響寬度較大，地下水位降深最大范圍300～500m，最大影響寬度5～7km。隧洞在采取工況三防滲措施后，其地下水位降深范圍及影響寬度減小明顯，地下水位降深最大范圍減小為150～200m，減幅約50%，最大影響寬度減小為3～5km，減幅約40%。

4.3 主要巖溶大泉影響程度預(yù)測(cè)

隧洞區(qū)域可溶巖地層分布廣，地下水環(huán)境影響敏感，隧洞穿越的馬耳山兩側(cè)盆地（鶴慶盆地和劍川盆地）山腳巖溶大泉分布密集。根據(jù)隧洞穿越區(qū)水文地質(zhì)結(jié)構(gòu)，主要泉點(diǎn)的分布位置，以及隧洞與泉點(diǎn)在平面和垂向上的位置關(guān)系，通過(guò)數(shù)值模擬計(jì)算，預(yù)測(cè)了在不同工況下枯水期和豐水期主要巖溶泉點(diǎn)的影響程度。不同工況下主要泉點(diǎn)影響程度可見(jiàn)表4。

由預(yù)測(cè)結(jié)果可知，裸洞施工工況下對(duì)距離隧洞較近的泉點(diǎn)影響較大，會(huì)疏干部分泉點(diǎn)，如清水江村泉。其余影響較大的泉點(diǎn)有黃龍?zhí)叮ㄓ绊懗潭却笥?0%），東山寺泉、水鼓樓泉、蝙蝠洞等（影響程度大于10%）;二次灌漿工況下對(duì)各泉點(diǎn)影響顯著下降，影響較大的泉點(diǎn)有清水江村泉（影響程度大于20%），對(duì)于其它泉點(diǎn)影響均減小至5%以下。

4.4 運(yùn)行期水位恢復(fù)情況預(yù)測(cè)

在開(kāi)挖完成后，地下水位下降和地下水量的疏干或外排過(guò)程將結(jié)束，且會(huì)因?yàn)榻涤耆霛B和側(cè)向補(bǔ)給等補(bǔ)給通道的輸入，進(jìn)入地下水恢復(fù)過(guò)程。為計(jì)算區(qū)域地下水位恢復(fù)時(shí)間，建立相應(yīng)區(qū)域地下水滲流場(chǎng)的非穩(wěn)定模型，以開(kāi)挖后的滲流場(chǎng)為初始條件，計(jì)算地下水滲流場(chǎng)恢復(fù)的程度和時(shí)間，以枯水期為例展開(kāi)說(shuō)明。

圖5為隧洞地下水位分別恢復(fù)2年和5年后的地下水滲流場(chǎng)情況。 由計(jì)算預(yù)測(cè)結(jié)果可知，隧洞完全封堵后，在恢復(fù)5年后地下水滲流場(chǎng)可基本恢復(fù)至開(kāi)挖前的狀態(tài)。需要說(shuō)明的是，以上計(jì)算是基于最為理想的邊界條件和補(bǔ)給條件下的恢復(fù)時(shí)間，實(shí)際的地下水水位線很難100%恢復(fù)至初始狀態(tài)。

5 ? 結(jié)論

本文以滇中引水工程某深埋長(zhǎng)隧洞為依托，通過(guò)建立GMS數(shù)值模型對(duì)不同施工工況下隧洞施工的涌水量、地下水水位變化、主要巖溶泉點(diǎn)影響程度、運(yùn)行期水位恢復(fù)情況等進(jìn)行了預(yù)測(cè)。

通過(guò)預(yù)測(cè)，在裸洞施工工況下對(duì)地下水環(huán)境影響較大，隧洞涌水量占比地下水總徑流量的18%左右，地下水位降深最大范圍300～500m，最大影響寬度5～7km，對(duì)距離隧洞較近的泉點(diǎn)影響較大，會(huì)疏干部分泉點(diǎn)。

在采用防滲灌漿措施后，效果明顯。二次灌漿（工況三） 條件下，預(yù)測(cè)隧洞涌水量占比地下水總徑流量減小到5%左右，地下水位降深最大范圍減幅約50%，最大影響寬度減幅約40%，對(duì)主要巖溶大泉的影響均降到了5%以下。

隧洞施工對(duì)地下水環(huán)境影響的研究貫穿于從勘察設(shè)計(jì)到施工、運(yùn)行的整個(gè)過(guò)程，本次預(yù)測(cè)研究為制定滲控標(biāo)準(zhǔn)等提供了合理、準(zhǔn)確、可行的依據(jù)。在施工過(guò)程中需要對(duì)設(shè)計(jì)階段的計(jì)算成果不斷地進(jìn)行反饋修正，以完善預(yù)測(cè)的準(zhǔn)確率，更好地服務(wù)于施工，為地下水環(huán)境保護(hù)和施工決策提供科學(xué)依據(jù)。

參考文獻(xiàn)/References

[1]傅冰駿.國(guó)際隧道及地下工程發(fā)展動(dòng)向[J].探礦工程， 2002（5）：54--57.

[2]戴文亭，白寶玉.我國(guó)隧道及地下工程發(fā)展現(xiàn)狀和前景展望，東北公路，2000 （4）：90—92.

[3]潘家錚，何憬.中國(guó)大壩五十年[M] .北京：中國(guó)水利水電出版社， 2000.

[4]唐運(yùn)剛.基于解析法的隧洞施工對(duì)地下水環(huán)境影響預(yù)測(cè)[J].人民長(zhǎng)江， 2018， 49（8）：67-71.

[5]蒙彥，雷明堂.巖溶區(qū)隧洞涌水研究現(xiàn)狀及建議[J].中國(guó)巖溶，2003，22（4）：287-292.

[6]龔星，陳植華，孫璐.地下水環(huán)境影響評(píng)價(jià)若干關(guān)鍵問(wèn)題探討[J].安全與環(huán)境工程，2013，20（12）：95—99.

[7]唐輝明，晏幫川，胡新麗.工程地質(zhì)數(shù)值模擬的理論與方法[M].武漢：中國(guó)地質(zhì)大學(xué)出版社，2001.

[8]陳長(zhǎng)生，周云，王家祥，張海平.復(fù)雜巖溶區(qū)深埋長(zhǎng)隧洞涌水量預(yù)測(cè)分析[J].人民長(zhǎng)江.2015， 46（10） ：29-32.

[9]徐棟棟，鄔愛(ài)清，孫玉杰.某水電站引水隧洞突水?dāng)?shù)值模擬[J].長(zhǎng)江科學(xué)院院報(bào).2010（08）.