巖溶區某高速公路隧道地下河最大涌水量的確定

覃蘭麗，張云良

（廣西壯族自治區桂林水文工程地質勘察院，廣西桂林541002）

巖溶區某高速公路隧道地下河最大涌水量的確定

覃蘭麗*，張云良

（廣西壯族自治區桂林水文工程地質勘察院，廣西桂林541002）

廣西大部分地區屬典型的喀斯特地貌，地下河發育，以西南巖溶地區某高速公路隧道穿越地下河為例，分析巖溶地下水動態特征，采用3種方法計算地下河最大涌水量，為設計、施工提供必要的水文參數，為采取科學的方法解決隧道水患問題提供基礎數據。

巖溶；隧道；地下河；最大涌水量

1　隧道及巖溶地下河概況

隧道為分離式長隧道，全長1070m，左、右線縱坡均為1.9%。位于峰叢洼地地貌，區域地下水類型屬裸露型巖溶水，富水性中等。2013年11月，開挖左線隧道時揭露巖溶地下河，水面標高408.4m，水深約25m；水清、無色無味；完全揭露后，巖溶地下河水面寬13～15m，走向與左線隧道走向呈40°斜交。抽水5h，水位未發生變化。為不影響工期，施工方采用洞渣直接回填巖溶地下河，在大、暴雨時節，地下水上涌，上涌最大高度為開挖路面標高（409.31m）以上約0.5m，沿隧道開挖面流出時間約1～5d。

2　巖溶地下水動態特征

區域地下水動態主要受大氣降雨的影響。區內巖溶地下水動態以劇變型為主。在降大、暴雨后，地表水、地下水由北向南匯入洼地；地下水位迅速上升，流量急劇增大。明流在北部洼地巖溶地下河伏流入口（S2）進入地下，穿過山體后，在東南部洼地巖溶地下河伏流出口（S1）流出，轉為明流。伏流段地下水排泄比明流排泄暢通性差，導致洼地形成內澇洪淹，洼地積水，水深1～2.0m。開挖隧道揭露的巖溶地下河流量特征與巖溶峰叢洼地、谷地巖溶地下河流量分布特征相似。表現在豐水期水量大，枯水期流量小，流量變化大，水位變幅大、時空分布不均等特點。

3　地下河最大涌水量計算

經調查分析，隧道開挖揭露的巖溶地下河是多條匯水于S2-S1巖溶地下河的其中一條。對其水位標高和流量影響較大的是隧道西、西北側巖溶山區，尤以降雨滲入、集中經洼地匯聚巖溶地下河的影響較大。本次做巖溶地下河連通試驗，在地下河流動天窗（S4）投入大量木屑和稻草（1～2cm），分別在第9d和26d在隧道內地下河（S3點）及S1點發現有木屑、稻草流出，從而確定S4-S3-S1為連通關系。

同一天實地測量S1水面標高為395.78m，S2水面標高為406.12m，S3水面標高為409.31m（路面最終標高409.80m），S4水面標高為410.42m；圖上測量S1至S2的直線距離為650m，S4至S1距離為890m，見圖1。自然狀態下S2至S1地下水水力坡度i=（406.12-395.78）÷650=1.6%。地下河天窗S4至伏流出口S1水力坡度i=（410.42-395.78）÷890=1.64%。長期的地質地貌發展和長期的巖溶管道發育，使地下水遵循阻力最小，最暢通的排泄通道運移，這符合自然規律，也符合該地下水流場。

下面用斷面法、降雨滲入匯流法及地下水徑向輻合和徑向輻射法分別計算地下河S3處100年一遇最大涌水量。

3.1斷面法

（1）地下河可通過的最大流量計算。巖溶地下河伏流出口S1處，已修建有人工擋水壩，壩上有2個矩形狀決口，橫斷面矩形寬B1=1.8m，高h1=0.45m，橫斷面B2=1.20m，高h2=0.45m，水深0.40m；伏流入口S2以上有多道漿砌石攔水壩，現已沖垮，形成近似矩形決口，寬1.1m，高0.8～1.20m，水深0.22m。利用以上數據計算隧道建成后地下河S3處可通過的最大流量Q′（測量時水面標高為409.31m，該處路面標高409.80m）。

圖1　水文計算簡圖（注：圖中標高均為水面標高）

采用矩形堰堰測流量法公式計算：

Q=0.018Bh3/2

式中：Q——過堰流量，L/s；

B——堰檻寬，cm；

h——過堰水深，cm；

0.018——流量系數。

①伏流出口S1流量：

Q1=0.018 B1h13/2+0.018 B2h23/2

=0.018×180×403/2+0.018×120×403/2

=1366.10（L/s）（即4917.97m3/h）

②伏流入口S2地表水流量：

Q2=0.018 B2h23/2

=0.018×110×223/2=204.31（L/s）（即735.53m3/h）

③伏流出口S1地下水流量：

Q=Q1-Q2=1366.10-204.31=1161.79（L/s）（即4182.44m3/h）

④經隧道地下河S3處，地下水過水斷面L1= 1600m（垂直地下水流向，至分水嶺），單位寬度（m）內過水流量q=Q÷L1=1161.79÷1600=0.726（L/s）（即2.614m3/h）。

⑤隧道巖溶地下河處，地下河水S3最大流量（Q′）。地下水上游匯水至隧道地下河S3處斷面長L2= 760m，由 Q′=L2×q=760×0.726=551.76（L/s）（即1986.34m3/h）。

（2）計算地下河水（S3）的最大涌水量。經實地調查并訪問幾位60～70歲村民，封閉巖溶洼地內（S4所在洼地）有村民種植的玉米，幾十年從未被水淹過，根據洼地標高及天窗S4內地衣和苔蘚生長特點，并實地測量，確定S4最高水位標高約為Hmax=413.70m（相當于50年一遇），因隧道設計年限為100年，將50年一遇的洪水標高再抬高1.0m，即H100=Hmax+1.0=414.70m（主要參考隧道出口北側村民老房子的標高，并含有一定的安全系數）。由此按地下水水力坡度1.64%計算，得出S1水位標高為H=414.70-（890×0.0164）=400.10（m），根據地形標高、攔水壩及近似矩形明流溝渠特點，可假定成B=2.0m，H=2.27m的矩形堰口。S2水位標高為H= 400.10+（650×0.016）=410.50（m），可假定成B=1.5m，H=1.25m的矩形堰口。再用矩形堰堰測流量法公式計算得：

①伏流出口S1流量：

Q1=12312.36L/s（即44324.49m3/h）

②伏流入口水S2流量：

Q2=3773.36L/s（即13584.09m3/h）

③伏流出口S1地下水流量：

Q=Q1-Q2=8539L/s（即30740.4m3/h）

④經隧道地下河S3處，地下水過水斷面L1= 1600m，單位寬度（m）內過水流量q=Q÷L1=8539÷ 1600=5.34（L/s）（即19.22m3/h）。

⑤隧道巖溶地下河處，地下河水S3最大流量（Q′）。地下水上游匯水至隧道地下河S3處斷面長L2= 760m，由 Q′=L2×q=760×5.34=4058.4（L/s）（即4.058m3/s）。

隧道內最大涌水量：

Q100=4058.4-551.76=3506.64（L/s）（3.506m3/s）。

3.2降雨滲入匯流法

Q100=F100·A·α·t-1

式中：Q100——100年一遇最大涌水量，m3/s；

F100——100年一遇最大日降水量，10-3m，按天等縣歷史實測記錄資料，最大日降水量172mm（1955年）代替；

A——巖溶地下河匯水面積，106m2；共2.535km2，其中巖溶區（A1）約2.0km2，碎屑巖和粘性土區（A2）約0.535 km2；

α——降水入滲系數，巖溶區取α1=0.45，碎屑巖區和粘性土區取α2=0.15；

t——實際調查降大、暴雨后，地下河水點出現最高水位（m）和最大涌水量時間，s；訪問和調查［在歷年的大、暴雨過程中，地下河伏流出口（S1）和大口井（S5）等主要地下水點水動態特征］獲得時間為12h。

代入上式計算得：

Q100=F100×（A1×α1+A2×α2）t-1=172×10-3×（2×106× 0.45+0.535×106×0.15）×（12×3600）-1=3.9028（m3/s）

3.3地下水徑向輻合和徑向輻射法

Q100=2πLHK/ln（r2/r1）

式中：Q100——100年一遇最大涌水量，m3/s；

π——圓周率，取3.14；

L——地下水滲流距離，m，巖溶地下河流動天窗處認為由滲流轉管道流補給邊界條件，隧道S3至S4流動天窗處取值300m；

H——滲流水頭，m，取值由地下河流動天窗S4地下水位標高414.70m，隧道巖溶地下河處地下水位水頭標高降至路面底板約407.00m（均質和網狀巖溶裂隙水均勻）為7.70m；

K——地下水所處塊段的滲透系數，按經驗值，類比為清潔的礫石層，取經驗值7×10-2cm/s（即7×10-4m/s）；

r2——徑向輻合和徑向輻射流地下水的影響半徑，m，取至上游補給邊界，為3250m；

r1——徑向輻合和徑向輻射流的半徑，m，取由隧道地下河S3至S4地下河流動天窗距離，即300m。

將數據代入上式，計算隧道巖溶地下河最大涌水量：

Q100=（2×3.14×300×7.7×7×10-4）÷ln（3250÷300）= 4.2621（m3/s）。

綜合以上，根據過水斷面法計算得到100年一遇最大涌水量為Q100=3.506m3/s，降雨滲入匯流法計算得Q100=3.9028m3/s，地下水徑向輻合和徑向輻射流法計算得Q100=4.2621m3/s，3種方法計算的結果吻合得較好，計算結果可信。

4　結語

隧道揭露的巖溶地下河屬于S2至S1地下河的一條支流。除其上游有個別地下水點露頭外，其余隱蔽于巖溶山體內，按照常規工作過程，需進行一個水文年的詳細水文地質專項勘察后，方能提供完整的水文地質參數。限于工程建設時間緊迫（沒有一個完整水文年），我們采用以地下河伏流出口S1以上整個地下河系統的流量及其特征，按匯水面積、流量、水位變化等綜合因素進行科學的預測，計算隧道巖溶地下河S3處最大涌水量，為設計、施工提供了必要的水文參數，為采取科學的方法解決隧道水患問題提供了基礎數據。

［1］《水文地質手冊》編纂委員會.水文地質手冊［M］.2版.北京：地質出版社，2012.

［2］《工程地質手冊》編員會.工程地質手冊［M］.4版.北京：中國建筑工業出版社，2007.

［3］張人權，梁杏，靳孟貴，等.水文地質學基礎［M］.6版.北京：地質出版社，2011.

［4］程俊賢，熊炳坤，孔繁業，等.1∶20萬靖西幅區域水文地質普查報告［R］.廣西壯族自治區水文工程地質隊，1977.

［5］陳崇希，林敏.地下水動力學［M］.武漢：中國地質大學出版社，1999.

［6］袁道先，蔡桂鴻.巖溶環境學［M］.重慶出版社，1988.

U45

A

1004-5716（2016）07-0202-03

2016-03-25

2016-03-29

覃蘭麗（1980-），女（壯族），廣西來賓人，工程師，現從事礦山水文地質詳查、地下水資源論證、地質災害評估等相關工作。