楊紅超，胡新麗，黃凱湘，李 群，楊裕云

（1.中國地質(zhì)大學工程學院，湖北 武漢430074；2.中交第二公路勘察設(shè)計研究院有限公司，湖北 武漢430050）

水是人類生存和發(fā)展的主要物質(zhì)，其中地下水是制約國民經(jīng)濟發(fā)展的主要資源［1—2］，尤其在干旱地區(qū)或在降水和地表水資源隨季節(jié)變化大的地區(qū)，地下水資源在國民經(jīng)濟建設(shè)中的影響更大。目前，對地下水開發(fā)利用的方式較多［3—7］，諸如：直接引用泉水；打井抽取地下水；通過渠道或滲溝導(dǎo)引河水，存入下游兩岸含水層中，形成地下水庫，待枯水期備用；在巖溶山區(qū)采用人工堵塞洞穴，形成地下水庫，引水灌溉發(fā)電；在平原區(qū)修筑地下截潛工程，抬高地下水位，使之直接流入渠道中，而廣闊渠截潛工程即是這一方案的規(guī)劃。因此，在截潛工程實施前有必要對研究區(qū)枯水期地下水資源進行預(yù)測和評價。

目前在沖積平原區(qū)采用截潛自流引水工程的成功實例極少，需要在查明區(qū)域和研究區(qū)水文地質(zhì)條件的基礎(chǔ)上，對該方案的可行性開展系列問題的研究。即首先需要對研究區(qū)地下水資源進行評價，同時需要對因截潛工程施工和運行過程中可能產(chǎn)生的環(huán)境地質(zhì)問題進行預(yù)測評估，如當截潛工程實施并形成地下水庫后，地下水位必將大范圍上升而導(dǎo)致浸沒災(zāi)害；因地下水庫形成以后其徑流滯緩，導(dǎo)致水質(zhì)惡化等。為此，本文以北汝河中下游沖積平原區(qū)廣泛渠截潛自流引水工程為例，在查明和分析研究區(qū)水文地質(zhì)條件的基礎(chǔ)上，根據(jù)研究區(qū)地下水資源的類型和水文地質(zhì)條件的差異分別選擇適宜的評價方法［1，7—8］，對截潛工程實施前枯水期地下水資源量進行了評價，即對該計算區(qū)枯水期地下水資源的儲存量（采用重力水體積法計算）、地下水的主要補給量（采用過水斷面法計算）、可供截潛自流開發(fā)利用的地下水徑流量（采用水均衡法計算）進行了計算，同時根據(jù)研究區(qū)地理及地質(zhì)環(huán)境的特點，對枯水期可供自流開發(fā)利用的地下水資源不足的原因進行了分析與評價，為該區(qū)截潛自流引水工程可行性分析提供了可靠的依據(jù)。

1 工程概況及地質(zhì)環(huán)境條件

1.1 工程概況

廣闊渠引水工程渠首位于河南省郟縣城關(guān)西北20km的北汝河左岸［9］，渠首底板高程為144m，設(shè)計引河流量為10m3/s。該工程于1978年竣工，全長23km，在農(nóng)田灌溉、補充下游水資源、改善生態(tài)環(huán)境等方面發(fā)揮了很好的作用。1982年，因特大暴雨洪水使北汝河主流河槽向西改道，導(dǎo)致廣闊渠喪失自流導(dǎo)引河水的功能。2001年實施工程擴建，自渠首向西經(jīng)河床至新河槽間修建了長522m的暗涵，其底板埋深為1.5～2m，進水口底板高程為144.68m，從此恢復(fù)了渠道自流導(dǎo)引河水的功能。

近年來，由于氣候變化的影響以及勘察區(qū)及其上游城鄉(xiāng)工農(nóng)業(yè)發(fā)展用水量增大，枯水期北汝河流經(jīng)本河段的流量逐年減小，有時河水幾近干涸，廣闊渠的引水量不能滿足郟縣城鄉(xiāng)工農(nóng)業(yè)發(fā)展的需要。為了增加廣闊渠枯水期的流量，郟縣主管部門提出了廣闊渠截潛自流引水工程方案，擬于暗涵下游河床至左岸漫灘修筑地下壩，攔截河床和漫灘中的地下水，使之自流進入暗涵渠道中，規(guī)劃增加流量1.0～1.5m3/s。

1.2 研究區(qū)地質(zhì)環(huán)境條件

1.2.1 地形地貌

研究區(qū)屬北汝河中下游沖積平原區(qū)，地形平緩開闊，總地勢由西北向東南傾斜，地面高程為153～146m。按微地形特征可分為河槽、河床和漫灘三個自然微地貌單元（見圖1）。其中，河槽寬為30～110m，高程為150～146m，是枯水期排泄河水和地下水的主要通道；河床寬約為500m，高出河槽數(shù)十厘米，枯水期常干涸無水；漫灘主要在河床東北側(cè)，分布寬廣，是當?shù)剞r(nóng)耕區(qū)，與河床最大高差約3m。

1.2.2 地層巖性

研究區(qū)全新世沖積物分布廣泛，厚6～11m，按巖性可分為以下四類：

（1）卵礫石：粒徑一般5～8cm，大者達20cm，其中卵礫組含量為70%～90%，砂粒組含量為10%～30%，厚6～11m，出露于河槽、河床，遍布研究區(qū)及鄰區(qū)。據(jù)現(xiàn)場水文地質(zhì)試驗，卵礫石屬強透水層，滲透系數(shù)為480m/d，重力給水度為0.2。

（2）細砂土：青灰、灰白色，疏松，顆粒較均勻，分布于本區(qū)東北部，厚3～4m，呈薄層狀透鏡體分布于卵礫石層中，透水性中等，滲透系數(shù)為10m/d。

（3）粉砂土：灰白色，松散，夾少量薄層黏性土，厚0.5～1.7m，分布于河堤南側(cè)漫灘地表，覆蓋在卵礫石層之上，滲透系數(shù)為2.54m/d，入滲系數(shù)為0.2。

（4）粉土：顆粒較均勻，疏松，厚0.5～2.5m，分布于河堤北側(cè)漫灘表部，覆蓋于卵礫石層和細砂土層之上，滲透系數(shù)為1.16m/d，入滲系數(shù)為0.2。

第四系下伏地層為上第三系洛陽組，巖性以黏土巖為主，下部夾砂巖和礫巖，分布廣，厚度較大，區(qū)內(nèi)鉆探揭露厚11.4m，滲透系數(shù)為0.000 47m/d，屬微弱透水性。

1.2.3 氣象及水文地質(zhì)

圖1 截潛工程區(qū)水文地質(zhì)平面圖Fig.1 Hydrogeological planar graph of the intercepting underflow project area

據(jù)郟縣氣象站1951—2004年觀測資料表明，本區(qū)年降水量變化較大，年均降水量為697.9mm，最大年降水量達1 133.8mm，最小年降水量僅356.0 mm。同時，年內(nèi)降水量分配不均，多集中在6～9月份，而旱季較長，一般從當年11月至翌年4月，有時延至6月（如2006年）。本區(qū)年均蒸發(fā)量為1 861.9mm，最大蒸發(fā)量為2 537mm，最小蒸發(fā)量為1 512.9mm。

研究區(qū)位于北汝河中下游地帶，北汝河是淮河水系的重要支流，發(fā)源于西北部嵩縣山區(qū)，出山口進入沖積平原區(qū)后，具游蕩型河流特征。據(jù)汝州水文站1978—1985年觀測資料表明，年內(nèi)流量差異極大，汛期最大流量為389m3/s，枯水期最小流量為0.02m3/s，有時幾近斷流，2006年6月13日在研究區(qū)西北邊界發(fā)現(xiàn)河槽干涸無水。

據(jù)研究區(qū)大量勘探測試工作，并結(jié)合地表調(diào)查訪問和鄰近地區(qū)水文地質(zhì)資料分析，研究區(qū)水文地質(zhì)條件具有以下主要特征。

（1）研究區(qū)巖土類型不多，但滲透性差異極大，滲透系數(shù)以卵礫石層最大，第三系黏土巖最小，其余介于兩者之間。

（2）按地下水的埋藏條件和動力特點，研究區(qū)地下水類型為孔隙潛水，主要含水層為分布廣泛的卵礫石層，富水性良好，水位埋藏較淺，枯水期地下水埋深僅2～4m。

（3）地下水的補給來源較多，諸如邊界西北和東北部地下潛流的補給，降水的補給，汛期洪水的補給以及灌溉回滲水的補給，但枯水季節(jié)補給來源較少，主要為來自邊界西北側(cè)地下潛流的補給。

（4）受區(qū)域總地勢、季節(jié)變化及人類活動的影響，地下水的運動方向多變，其運動總趨勢為西北向東南，靠岸邊附近漫灘中的地下水補給河水，而汛期河水向岸邊運動補給地下水。據(jù)現(xiàn)場實測，地下水水力坡度較小，大都在0.003左右，地下水流動較緩慢。

（5）地下水的排泄方式較多，一般以地下潛流向東南下游區(qū)外排泄為主，其次為眾多民井抽水、通過暗涵向廣闊渠和河槽排泄以及蒸發(fā)消耗等。

2 研究區(qū)截潛工程實施前地下水資源量的計算

在截潛工程實施前，需在查明勘察區(qū)約4.13 km2范圍內(nèi)的水文地質(zhì)條件基礎(chǔ)上，對該區(qū)枯水期地下水資源的數(shù)量進行分析計算，以便為截潛工程方案可行性提供依據(jù)，同時為下一步對研究區(qū)截潛工程效益的數(shù)值模擬提供計算參數(shù)的初值。

本次地下水資源計算以勘察區(qū)和鄰近地區(qū)的水文地質(zhì)參數(shù)以及2006年6月13日實測編制的地下水潛水等水位線圖為依據(jù)，主要針對地下水儲存資源、補給資源和可供截潛自流開發(fā)利用的徑流資源三部分的數(shù)量進行計算。

2.1 地下水儲存量

地下水的儲存量是指地下水位以下含水層中重力水的體積，常以m3表示。它受多種因素的影響，如勘察區(qū)范圍大小、含水層的巖土類型、含水層的性質(zhì)（潛水或承壓水）、含水層的厚度及季節(jié)變化等。

研究區(qū)地下水主要潛水，在枯水季節(jié)主要儲存于卵礫石含水層中，厚5～9m，地下水位埋深較淺，一般為2～4m。為了查明枯水期潛水的儲存量QA，以2006年6月13日的地下水水位等值線圖為依據(jù)，采用以下公式計算含水層中重力水的體積：

式中：μ為卵礫石層的重力給水度，按0.2計；V為含水層體積（m3）。

在研究區(qū)4.13km2范圍內(nèi)，因含水層厚度各地不一，因此含水層的體積需分區(qū)計算，其計算公式為

式中：hi為各區(qū)含水層厚度（m）；Si為各區(qū)面積（m2）。

在4.13km2范圍內(nèi)，按式（2）計算得到研究區(qū)含水層總體積為3 149.96×104m3，將其代入式（1），可計算得到研究區(qū)枯水期地下水儲存量為629.99×104m3。

2.2 地下水補給量

一般說，補給潛水的水源較多，其中降雨入滲是主要的補給來源。據(jù)郟縣氣象站資料，2006年6月22日以前降雨很少，大部分時間干旱無雨，6月13日實測全區(qū)地下水位時發(fā)現(xiàn)勘察區(qū)西北邊界附近河槽干涸無水，因此在這種條件下降雨入滲對地下水的補給可以不予考慮。此外，勘察區(qū)東北和西北邊界大致與地下水流向平行，其兩側(cè)地下徑流的補給甚微，其補給量也可忽略不計。

由于研究區(qū)內(nèi)地下水運動的總趨勢由西北流向東南，枯水期本區(qū)地下水的主要補給量來自西北側(cè)邊界的地下徑流，因此計算地下水補給量時的計算斷面選在西北邊界附近?？紤]到該邊界附近含水層厚度和水力坡度的差異，將過水斷面分為4段（見圖2），其補給量QB按下式采用過水斷面法分段計算：

式中：QB為西北邊界地下徑流對本區(qū)潛水的總補給量（m3/d）；Ji、Fi分別為某計算段的平均水力坡度和過水斷面面積（m2）；Ki為某計算段卵礫石層的滲透系數(shù)，各段均按480m/d計。

圖2 研究區(qū)上游邊界地下水補給量計算示意圖Fig.2 Computation graph of groundwater recharge in the upstream area

按照式（3）計算研究區(qū)枯水期各段地下水的補給量，其結(jié)果見表1。

表1 各段地下水補給量計算結(jié)果Table 1 Calculation results of the groundwater increment

由表1可見，枯水期補給本區(qū)地下水的來源很少，主要來源于西北側(cè)邊界地下徑流的補給，其補給量為42 663.1m3/d（0.494m3/s）。

2.3 可供截潛自流開發(fā)利用的地下水徑流量

可供截潛自流開發(fā)利用的地下水徑流量Qk包括兩部分：一是地下水通過暗涵進入渠道的流量；一是通過暗涵和渠道底板以下及暗涵進水口西南側(cè)含水層排向下游區(qū)外的潛水流量。根據(jù)研究區(qū)水文地質(zhì)特征及大量勘察測試資料，其枯水期流量可采用水均衡法計算，當均衡區(qū)為4.13km2、均衡期以日計時，區(qū)內(nèi)潛水體積不變，則以各種方式對本區(qū)的補給量與支出量相等，可用下式表示：

式中：QB為西北邊界補給本區(qū)的潛水流量，為42 663.1m3/d；Qb為農(nóng)田潛溉民井抽排的流量，據(jù)郟縣水利局現(xiàn)場調(diào)查訪問，為8 433m3/d；Qn為地下水通過暗涵排入渠道的流量，實測為13 824 m3/d；Qs為通過暗涵和渠道底板以及暗涵進水口西南側(cè)含水層排向下游區(qū)外的潛水流量（m3/d）；Qh為灌溉回歸水入滲量，取其灌溉入滲系數(shù)為0.2，則入滲量為1 686.6m3/d；Qg、Qz分別為降雨量和蒸發(fā)量，據(jù)本均衡區(qū)和均衡期的特點，令Qg＝Qz，則（4）式可寫成：

將上述各收支項的參數(shù)代入式（5），則有

可供截潛自流開發(fā)利用的地下水徑流量Qk由下式計算：

3 截潛自流引水方案的可行性分析

通過對研究區(qū)截潛自流引水工程實施前枯水期地下水資源量的計算，可對該區(qū)廣闊渠截潛自流引水方案的可行性進行如下分析論證：

（1）在研究區(qū)4.13km2范圍內(nèi)，枯水期地下水的儲存量為629.99×104m3，相當于一座小型水庫的容積，說明本區(qū)地下水儲存資源較豐富。

（2）枯水季節(jié)補給研究區(qū)地下水的來源很少，主要來自本區(qū)西北側(cè)地下徑流，經(jīng)計算其總補給量為42 663.1m3/d（0.494m3/s）；而規(guī)劃實施的地下截滲墻抬高地下水位后，要求自流進入暗涵渠道的流量為1.0～1.5m3/s，可見所計算的補給量不能滿足規(guī)劃用水的要求。分析原因認為：研究區(qū)地形平緩開闊，地下水的天然水力坡度較小，多為0.003左右，加之西北側(cè)主要補給區(qū)農(nóng)田廣布，枯水期正值農(nóng)田灌溉期需抽取大量地下水，因而影響了對研究區(qū)的補給量。

（3）枯水期研究區(qū)地下水的補給來源和補給量有限，加上區(qū)內(nèi)農(nóng)田灌溉用水的消耗，可供自流開發(fā)利用的地下水流量更小，僅為0.416m3/s，遠不能滿足規(guī)劃用水的要求。

（4）渠首和暗涵的底板高程已先后確定，暗涵埋深較淺，凡水位在144.68m以下的地下水不能自流進入暗涵渠道中，這是截潛工程實施前在枯水期地下水水位下降與河水流量銳減共同影響下渠道流量減少的原因。

（5）當截滲墻實施后，預(yù)計此種條件下對研究區(qū)地下水的補給量將小于0.494m3/s；同時，因區(qū)內(nèi)外含水層存在良好的水力聯(lián)系，在截滲墻兩端必將產(chǎn)生繞墻流向下游，自流進入暗涵渠道的流量將會小于0.416m3/s；此外，截潛工程實施后，因地下水位上升將會導(dǎo)致大范圍沼澤化、鹽漬化和水質(zhì)污染等環(huán)境地質(zhì)問題。因此，為了滿足枯水期規(guī)劃用水的要求，在研究區(qū)采用截潛自流引水方案是不適宜的，需擴大研究思路，探討其他經(jīng)濟合理、水源可靠、環(huán)境影響小的供水方案。經(jīng)方案利弊比較［9］，考慮到研究區(qū)北汝河汛期流量較大，建議沿廣闊渠附近開挖若干大型蓄水池，通過渠道直接導(dǎo)引北汝河汛期河水存入池內(nèi)，待枯水期備用。目前該方案已被郟縣主管部門采納。

4 結(jié)論與建議

（1）經(jīng)現(xiàn)場勘探、試驗及調(diào)查訪問，查明了研究區(qū)工程地質(zhì)環(huán)境條件，為本區(qū)地下水資源評價提供了可靠依據(jù)。

（2）采用水均衡法和解析法分別對本區(qū)枯水期地下水補給量、儲存量和可供自流開發(fā)的徑流量進行計算，結(jié)果表明采用截潛自流引水的方案不能滿足枯水期規(guī)劃用水的要求。

（3）為了滿足本區(qū)枯水期用水要求，需探討其他經(jīng)濟合理的方案。考慮到本區(qū)北汝河汛期流量較大，建議沿廣闊渠附近開挖若干大型蓄水池，將汛期河水由渠道引入池內(nèi)，待枯水期備用。該方案已被郟縣主管部門采納。

［1］王增銀.供水水文地質(zhì)學［M］.武漢：中國地質(zhì)大學出版社，1994.

［2］武選民，文冬光，吳登定，等.我國主要平原盆地地下水資源可持續(xù)利用的幾點思考［J］.水文地質(zhì)工程地質(zhì)，2005，32（6）：96－99.

［3］李濤，鄧銘江，王于寶，等.臺蘭河山前地下儲水構(gòu)造及地下水庫可行性研究［J］.水文地質(zhì)工程地質(zhì)，2011，38（1）：30－34.

［4］康衛(wèi)東，仇小強，李文鵬，等.秦嶺山前截洪引滲與地下水庫調(diào)蓄及其協(xié)同效應(yīng)研究［J］.水文地質(zhì)工程地質(zhì)，2011，38（2）：8－13.

［5］李硯閣.地下水庫建設(shè)研究［M］.北京：中國環(huán)境科學出版社，2007.

［6］張兆吉，費宇紅，陳宗宇，等.華北平原地下水可持續(xù)利用調(diào)查評價［M］.北京：地質(zhì)出版社，2009.

［7］何師意，周錦忠，曾飛躍.巖溶地下河流域地下水資源評價——以湖南湘西大龍洞為例［J］.水文地質(zhì)工程地質(zhì)，2007，34（5）：33－36.

［8］楊紅超，胡新麗，李群，等.郟縣截潛自流引水工程環(huán)境負效應(yīng)及對策研究［J］.工程地質(zhì)學報，2012，20（6）：992－997.

［9］楊裕云，王亮清，李群，等.河南省郟縣廣闊渠截潛引水工程枯水期水資源數(shù)值模擬研究［R］.武漢：中國地質(zhì)大學工程學院，2006.