有入滲時地下水位變動下的分水嶺移動規(guī)律

胡 強(qiáng)， 肖先煊， 袁文杰， 蔡國軍， 梅 杰

(成都理工大學(xué) 環(huán)境與土木工程學(xué)院； 國家級地質(zhì)工程實驗教學(xué)示范中心, 成都 610059)

·儀器設(shè)備研制與開發(fā)·

有入滲時地下水位變動下的分水嶺移動規(guī)律

胡 強(qiáng)， 肖先煊， 袁文杰， 蔡國軍， 梅 杰

(成都理工大學(xué) 環(huán)境與土木工程學(xué)院； 國家級地質(zhì)工程實驗教學(xué)示范中心, 成都 610059)

根據(jù)相似原理，在室內(nèi)建立能模擬地下水分水嶺移動特征的模型。模擬一定降雨條件下，保持一側(cè)河渠水位不變，改變另一側(cè)水位，并加入示蹤劑，觀察不同水位差時地下水運動特征并記錄測壓管水位，繪制潛水位與分水嶺形態(tài)，實測分水嶺位置，得到試驗數(shù)據(jù)，并對比不同水位差下地下水分水嶺的移動特點。同時，根據(jù)已知的模型參數(shù)，由河間地塊潛水流浸潤曲線方程與分水嶺位置計算式得到理論潛水位和分水嶺位置，并與實測值相比。結(jié)果表明：通過公式計算得到的理論潛水位和分水嶺位置與實際試驗測得的數(shù)據(jù)相吻合，但存在臨界值Cr，當(dāng)兩側(cè)水位差Δh>Cr時不再適用，臨界值大小與流體密度、動力黏性系數(shù)、滲流斷面平均流速、砂樣平均粒徑等有關(guān)。

潛水; 入滲補(bǔ)給; 河間地塊; 水位差; 分水嶺移動

0 引 言

近年來，我國水資源形勢日趨緊張，政府批準(zhǔn)了一系列引水調(diào)水工程緩解城市水資源匱乏的現(xiàn)狀，如滇中引水工程[1]、南水北調(diào)工程。我國一些糧食主產(chǎn)區(qū),如河南、山東近年來多次遭遇干旱天氣，農(nóng)業(yè)部門積極興修水利以消除旱情對農(nóng)業(yè)生產(chǎn)的影響。水庫選址也要充分考慮蓄水后側(cè)向滲漏問題。在一系列水利工程建設(shè)中有很大可能出現(xiàn)地下水分水嶺移動引起水文地質(zhì)災(zāi)害或環(huán)境問題[2-4]。研究地下水潛水位及分水嶺的移動規(guī)律，進(jìn)而通過人工控制水位[5]，防止水庫發(fā)生側(cè)向泄漏[6-8]和河間地塊鹽漬化[9-11]或沼澤化[11]。城市垃圾填埋場或核廢料埋置點選址，也要充分研究地下水分水嶺運動規(guī)律防止發(fā)生泄漏造成污染[12]。因此，研究有入滲補(bǔ)給時地下水水位及分水嶺移動規(guī)律具有一定意義。

前人對河間地塊潛水位與分水嶺位置的研究多是通過理論推導(dǎo)而來，通過相應(yīng)理論公式計算得到的潛水位和分水嶺位置與實際情況是否相近需要實踐證明。本文試圖在前人研究的基礎(chǔ)上，在實驗室建立河間地塊砂槽模型，一定入滲條件下，改變兩側(cè)河渠水位高度，通過測壓管和示蹤劑，測得河間地塊潛水位和分水嶺位置，比較實測潛水位和分水嶺位置與公式理論計算值，驗證公式實用性，為實際生產(chǎn)中確定分水嶺位置和控制潛水位提供理論依據(jù)。

1 試驗?zāi)Ｐ?/h2>

薛禹群等[11]提出了潛水流浸潤曲線方程和河渠間有入滲時，河渠間分水嶺位置計算式:

(1)

(2)

式中：x為距左端起始斷面的距離；h為河間地塊距左端起始斷面x處的潛水流厚度；h1為左側(cè)河渠邊河流高度；h2為右側(cè)河渠邊河流高度；L為兩河渠之間的距離；K為河間地塊含水層的滲透系數(shù)；W為入滲強(qiáng)度；a為分水嶺位置。

上述公式推導(dǎo)條件為：① 含水層均質(zhì)各向同性，底部相對隔水層水平，上部有均勻入滲，可用入滲強(qiáng)度W表示；② 河渠間彼此平行，應(yīng)用裘布依假設(shè)，通過垂向平均，潛水流可視為一維流，符合達(dá)西定律；③ 潛水流是漸變流并趨于穩(wěn)定[11]。此外兩側(cè)河流為定水頭邊界。因為在河間地塊研究剖面內(nèi)河渠的水位是由剖面以上大面積的、多種形式的水量補(bǔ)給所決定的,而受局部范圍潛水排泄補(bǔ)給的影響很小，所以，這種假定是允許的[13]。

以均質(zhì)各向同性河間地塊為模擬對象，建立砂槽模型(見圖1)。建模過程中，遵循幾何相似原則，即① 幾何相似；② 運動相似；③ 動力相似；④ 邊界條件相似[14-15]。

圖1 河間地塊潛水砂槽模型

河間地塊潛水模擬箱(長×寬×高)1.9 m×0.4 m×0.8 m。單側(cè)設(shè)置9排測壓管，每排1支，管徑φ=10 mm,L=1 m，測壓管距左側(cè)河流的距離分別為0，20,40.7, 61, 75.7, 90.5, 111, 131 ,148 cm。同側(cè)有8支25 mL的示蹤劑注入管, 為增大試驗測量分水嶺位置的精確度，在模擬箱頂部增加5支示蹤劑注入管。蓄水箱容量(長×寬×高)=1.8 m×0.5 m×0.35 m。模擬箱頂部安裝雨篩，使給水能模擬區(qū)域降雨，實驗中要確保模擬降雨全部入滲，不能形成地表徑流。為控制左、右側(cè)河渠水位，在滲流試驗箱體的左、右端附有可通過升降系統(tǒng)控制高程的排水溢流箱[14]。模擬箱內(nèi)接受降雨入滲的區(qū)域面積為2 910.48 cm2。試驗填充為均質(zhì)石英砂，粒徑0.1～1.0 mm，通過篩分實驗得到粒徑級配曲線(見圖2)。得到d50=0.6 mm,Cu(不均勻系數(shù))=5.14，Cc(曲率系數(shù))=2.1。裝置所用石英砂為級配良好的不均勻土[16]。

圖2 粒徑級配曲線

表1 滲流試驗數(shù)據(jù)

圖3 毛細(xì)水上升高度隨時間變化曲線

2 結(jié)果處理

河間地塊為石英砂填充，試驗中不予變化，因此在本次試驗中K值不變，K=0.238 cm/s。打開抽水泵，在模擬箱頂部形成降雨，試驗中保證降雨全部均勻入滲，不形成地表徑流。通過升降裝置調(diào)節(jié)兩側(cè)水位高度h1、h2,實驗中h2水頭保持不變。待溢流箱流出的水量穩(wěn)定后，測量一定時間t內(nèi)左、右兩側(cè)溢流箱水量V1、V2，采用水量時間法計算裝置流量Q=(V1+V2)/t，得到入滲強(qiáng)度W=Q/S，S為入滲區(qū)面積。加入示蹤劑，觀測示蹤劑運動特征和分水嶺位置，測量分水嶺到左側(cè)河渠距離a1，觀察測壓管記錄潛水水位，將試驗參數(shù)代入式(1)、(2)，得到理論潛水位與理論分水嶺位置a2，用計算機(jī)繪制潛水位。重復(fù)試驗3次，觀測分水嶺位中間、偏移、消失的現(xiàn)象。將實際測量結(jié)果與式(1)、(2)理論計算結(jié)果比較，驗證其實用性。3次試驗數(shù)據(jù)記錄見表2。從表2可知，3次試驗計算得到W分別為0.033、0.034、0.034 cm/s，可以認(rèn)為試驗中入滲強(qiáng)度W保持不變。通過實測和公式計算分別得到實測分水嶺a1與理論分水嶺位置a2，見表2。

2.1 兩側(cè)水位相等

調(diào)節(jié)溢流箱高度，h1=37.7 cm、h2=37.6 cm,因讀數(shù)誤差和儀器精度，可以認(rèn)為河流水頭相等。試驗現(xiàn)象見圖4。由圖5可見，實測水位h與理論水位H存在差值,但均不超過5%，可以認(rèn)為結(jié)果一致。分析誤差由讀數(shù)誤差和測壓管水頭損失以及儀器精度等造成。試驗表明，當(dāng)兩側(cè)水位相等時，通過式(1)、(2)計算分水嶺位置、潛水位高度與實際相吻合。

表2 3次試驗數(shù)據(jù)記錄

圖4 兩側(cè)河水位相等時地下水運動軌跡圖

圖5 兩側(cè)水位相等時地下水位曲線

2.2 兩側(cè)有水位差且分水嶺存在

待第1次試驗所加入示蹤劑完全消失后，升高左側(cè)溢流箱高度，待流量穩(wěn)定后，測得h1= 48.8 cm、h2=38 cm,由于試驗裝置排水孔較小，h2由37.6 cm升高至38 cm。試驗現(xiàn)象見圖6，可以觀測到，分水嶺向水位較高的左側(cè)發(fā)生偏移。

圖6 兩側(cè)有水位差且分水嶺存在地下水運動軌跡

由圖7可知，實測水位h與理論水位H相對誤差不超過5%，考慮到儀器精度、讀數(shù)誤差以及測壓管水頭損失，可以認(rèn)為理論結(jié)果與實測數(shù)據(jù)一致。實測分水嶺位置a1=53.80 cm,通過式(2)計算得到理論分水嶺。位置a2=51.56 cm，兩者相差為2.24 cm, 因儀器精度和讀數(shù)誤差，可以認(rèn)為a1=a2。試驗表明，當(dāng)兩側(cè)河渠有水位差且分水嶺存在時，通過式(1)、(2)計算分水嶺位置、潛水位高度與實際相吻合。

圖7 兩側(cè)有水位差且分水嶺存在地下水位曲線

2.3 兩側(cè)有水位差且分水嶺消失

待示蹤劑完全消失后，繼續(xù)升高左側(cè)溢流箱高度，待流量穩(wěn)定后，測得h1= 54.7 cm、h2=38 cm，試驗現(xiàn)象見圖8。由圖9可知，實測結(jié)果與理論計算結(jié)果差別較大。隨著兩側(cè)水位差Δh繼續(xù)增大，潛水流浸潤曲線方程和河渠間分水嶺位置計算式與實際結(jié)果相差較大，不再適用。

圖8 兩側(cè)有水位差且分水嶺消失地下水運動軌跡

圖9 兩側(cè)有水位差且分水嶺消失地下水位曲線

2.4 臨界水位差Cr

由3次試驗可知，兩側(cè)河渠水位差Δh增大，分水嶺向水位較高一側(cè)偏移。同時，隨著Δh增大，地下水流速v增大，在試驗2.3中出現(xiàn)公式便不再適用，由達(dá)西定律計算得到試驗2.3中的流速v=270 μm/s與通過水量時間法計算地下水流速v=500 μm/s相差較大。根據(jù)雷諾數(shù)Re計算式[19]，計算試驗2.3中Re，其中ρ=1 g/cm3，v=500 μm/s，d=0.6 mm，10℃下水的黏性系數(shù)[19]μ=1.308 Pa/s。計算得到Re=22.9>10。因此，隨著Δh的增大,試驗條件下Re超過達(dá)西定律適用范圍雷諾數(shù)1Cr時，式(1)、(2)不再適用。在試驗條件下，計算得到Re=10時，地下水流速v=190 μm/s，根據(jù)達(dá)西定律得到Δh=VL/K=13.55 cm，即Cr=13.56 cm。

從表6可知，試驗3.1、試驗3.2中ΔhCr時，不能在實際中利用潛水浸潤曲線方程和分水嶺計算公式來確定分水嶺位置與潛水位。

表6 水位差Δh與理論計算吻合度

3 結(jié) 論

通過本研究，可獲得以下結(jié)論與認(rèn)識：

(1) 在模型試驗中，河間地塊水位受到兩側(cè)河渠水位高度影響，在均質(zhì)各向同性、底部隔水層水平河間地塊中，入滲強(qiáng)度不變，地下水運動滿足一維穩(wěn)定流條件下，通過式(1)計算得到潛水水位與實際潛水位吻合，因此，可用式(1)來計算潛水位。

(2) 在均質(zhì)各向同性、底部隔水層水平的河間地塊中，入滲強(qiáng)度不變，地下水運動滿足一維穩(wěn)定流條件下，分水嶺位置受到兩側(cè)水位高度的影響，偏向水位較高一側(cè)。在已知滲透系數(shù)、入滲強(qiáng)度時即可通過式(2)來確定分水嶺位置。

(3) 一定入滲強(qiáng)度下，兩側(cè)水位差影響公式的實用范圍。存在一個臨界值,當(dāng)水位差大于臨界值時，不再利用潛水浸潤曲線方程和分水嶺計算式來確定潛水位和分水嶺位置。臨界值受流體密度ρ和動力黏性系數(shù)μ、滲流斷面平均流速v、砂樣的平均粒徑d影響，由式Cr=ReμL/(Kdρ)確定，其中Re=10。

本研究可為均質(zhì)各向同性河間地塊中兩側(cè)河渠間距確定、河間地塊潛水位控制等問題提供理論依據(jù)。

[1] 瞿霜菊，黃 輝，曹正浩. 云南省滇中引水工程規(guī)劃研究[J].人民長江，2013,44(10)：80-83.

[2] 焦致嫻. 滇中引水工程對水源區(qū)及水源工程下游段陸生生態(tài)評價區(qū)植被影響初步研究[D].武漢：華中師范大學(xué),2012.

[3] 馮 健，周懷東，彭文啟,等.滇中引水工程對洱海水環(huán)境影響[J].中國水利水電科學(xué)研究院學(xué)報，2012, 10(4):241-246.

[4] 廖曉超，許 模，蔣 莉.中壩隧道地下分水嶺演化的數(shù)值模擬分析[J].地下水，2014(1):48-49.

[5] 陳文芳.中國典型地區(qū)地下水位控制管理研究[D].北京：中國地質(zhì)大學(xué)(北京),2010.

[6] 夏金梧，吳玉華，羅小杰,等.江坪河水利樞紐庫首右岸河間地塊巖溶滲漏分析[J].人民長江, 1994(8):41-45.

[7] 周志芳，李 艷.清江隔河巖水利工程河間地塊滲漏估算[J].河海大學(xué)學(xué)報, 1996(5):101-104.

[8] 陳運東.重慶市藻渡水庫左岸河間地塊滲漏評價[J].工程地質(zhì)學(xué)報，2014(22)：244-248.

[9] 張 紅，章光新，楊建鋒,等.人工控制潛水位下蘇打鹽漬土水鹽動態(tài)分析[J].干旱區(qū)資源與環(huán)境,2008,22(5):149-154.

[10] 趙海卿.吉林西部平原區(qū)地下水生態(tài)水位及水量調(diào)控研究[D].北京：中國地質(zhì)大學(xué)(北京),2012.

[11] 薛禹群，吳吉春.地下水動力學(xué)[M].3版.北京：地質(zhì)出版社，2011：53-56.

[12] 鄭 佳.北京西郊垃圾填埋場對地下水污染的預(yù)測與控制研究[D].北京：中國地質(zhì)大學(xué)(北京),2009.

[13] 張明泉.河間地塊潛水分水嶺移動規(guī)律的初步探討[J].蘭州大學(xué)學(xué)報,1986(4):112-117.

[14] 肖先煊，許 模，蔡國軍,等.抽水試驗?zāi)Ｐ团c地下水動力學(xué)實驗教學(xué)試驗[J].實驗室研究與探索，2014,33(1):164-168.

[15] 虞修競，蔡國軍，付小敏,等.水文地質(zhì)實驗裝置的研制及應(yīng)用[J]. 實驗室研究與探索，2011,30(3)：209-212.

[16] 陳仲頤，周景星，王洪瑾.土力學(xué)[M].北京：清華大學(xué)出版社,1994：3-7.

[17] 虞修競，付小敏，蘇道剛.水文地質(zhì)學(xué)實驗教程[M].成都：西南交通大學(xué)出版社,2008：30-34.

[18] 張人權(quán)，梁 杏，靳孟貴,等.水文地質(zhì)學(xué)基礎(chǔ)[M].6版.北京：地質(zhì)出版社,2011：20.

[19] 柯 葵，朱立明.流體力學(xué)與流體機(jī)械[M].上海：同濟(jì)大學(xué)出版社,2009：237.

The Test of Watershed Movement Regulation Induced by Water Table Fluctuation When Infiltrated

HUQiang,XIAOXianxuan,YUANWenjie,CAIGuojun,MEIJie

(Collage of Environment and Civil Engineering; Laboratorial Center for Geo-engineering,Chengdu University of Technology, Chengdu 610059, China)

On the condition of constant infiltration strength, how to acquire the level of groundwater and the watershed position between two parallel ditches are worth studying. Previous scholars issued the theoretical formulas of water table curve and watershed position between two parallel ditches. However, the practice of formulas still need verifying. As a result, we built a model which could imitate the watershed movement indoor according to the similar principle. On the precondition of stable recharge, we kept one side water level constant and changed the other side level. By observing the path of groundwater after adding to tracer, we noted the level of piezometer tube, drew the form of watershed, and compared the watershed characteristic in the different water level. Meanwhile, we obtained a group of academic data by the formulas of water table curve and watershed position and compared with test data. The result showed we could acquire the level of groundwater and the watershed position by formulas. But there is a critical value if water level difference beyond a threshold the formulas will not work. The critical value is related to fluid density, dynamic viscosity coefficient, average flow velocity and average particle size.

phreatic water; infiltrated recharge; parallel ditch; water level difference; watershed moving

2016-05-16

國家自然科學(xué)青年基金項目(41202209)；成都理工大學(xué)中青年骨干教師培養(yǎng)計劃(HG0092)；省級專業(yè)綜合改革與實踐教學(xué)項目(SZH1103ZX08、16Z0547)

胡 強(qiáng)(1995-)，男，山西交城人，本科生，地下水科學(xué)與工程專業(yè)。

Tel.:18428381167； E-mail.：2540746843@qq.com

肖先煊(1985-)，男，江西吉安人，講師，主要從事地下水實驗教學(xué)工作。

Tel.:028-84077667； E-mail.：sassoon1123@126.com

P 641.2

A

1006-7167(2017)02-0014-04