黃土丘陵區淤地壩作用下地下水位動態演化特征

白樂 殷會娟 李昭悅

摘 要：以黃土丘陵區典型淤地壩刀勞慶骨干壩為例，采用野外原位觀測及相關分析的方法，研究淤地壩作用下地下水位時空分布格局和動態演化特征，揭示上下游補排關系及驅動因素，為淤地壩壩區地下水資源保護與高效利用提供理論依據。結果表明：①淤地壩下游和壩區地下水位年最大變幅分別為1.82 m和1.27 m，壩下游地下水位整體呈下降趨勢，而壩區內地下水位呈上升趨勢;②淤地壩區和下游兩觀測點地面高差9 m，直線距離478 m，淤地壩地下水位比下游地下水位高7.51 m，淤地壩壩區內地下水補給壩下游地下水;③淤地壩下游地下水位變化具有明顯的季節性，形成高峰（非雨季）與低谷（雨季），降雨補給地下水具有時滯性，壩內夏季地下水位上升、其他季節下降;④農業灌溉用水量增加是地下水位變化的重要驅動力，其對壩下游地下水位的影響大于壩區。

關鍵詞：淤地壩;地下水位;降雨;動態演化特征;黃土丘陵區

Abstract：Taking Daolaoqing key dam of typical warping dam in loess hilly area as an example， the spatial and temporal distribution pattern and dynamic evolution characteristics of groundwater level under the action of warping dam were studied by field in-situ observation and correlation analysis. The relationship between upstream and downstream recharge and discharge and driving factors were revealed， which provided theoretical basis for the protection and efficient utilization of groundwater resources in warping dam area. The results show that a） the maximum variation of groundwater level in the downstream of the warping dam and the dam area is 1.82 m and 1.27 m respectively. The overall trend of the groundwater level in the downstream of the warping dam is downward and the groundwater level in the dam area has an upward tendency; b） the ground elevation difference between warping dam area and downstream two observation points is 9 m， straight line distance is 478 m， the groundwater level of warping dam is 7.51m higher than that of downstream water level and the groundwater of the dam area recharges the downstream groundwater; c） the variation of groundwater level in the downstream of warping dam has obvious seasonality， forming peak （non-rainy season） and trough （rainy season）， and rainfall recharge groundwater has time delay. The groundwater level in spring and autumn of the dam falls in winter and rises in summer and; d） the increase of agricultural irrigation water consumption is an important driving force of groundwater level change and its influence to the groundwater level in dam area is not as good as that in downstream dam.

Key words： warping dam; groundwater level; rainfall; dynamic evolutionary characteristics;loess hilly region

淤地壩是黃土高原水土保持的關鍵措施，大規模淤地壩建設在發揮攔沙蓄水、調節徑流[1-4]、改善生態環境等綜合功能的同時，對地下水文循環產生了深刻影響，給地下水資源保護和生態治理帶來新的挑戰。地下水位是地下水循環的“脈搏”，是衡量區域地下水資源開發合理與否的重要判據。黃土高原水資源短缺、生態環境脆弱，大規模淤地壩建設對地下水位動態特征的影響，已逐步成為水文、環境、生態等學科研究的熱點。

眾多學者從淤地壩對地下水的影響機理、補排關系、影響比例等方面進行研究，例如：高海東等[5-7]闡述了淤地壩作為障礙性節點，對徑流過程和地下水產生影響;袁水龍等[8]指出淤地壩建設明顯改變壩前壩后的水動力過程，而且改變了溝道的侵蝕動力分布，進而對地下水循環產生影響;畢慈芬等[9]研究指出淤地壩建設可形成地下水庫，調節區域水環境。Huang等[10-12]研究表明黃土高原地下水補給主要發生在7—9月的雨季;綦俊諭等[13]運用3種計算方法研究淤地壩減水作用，結果表明淤地壩等工程措施對岔巴溝流域減水作用僅為7%左右，但流域內地下水占比增加約20%;黃金柏等[14]開發了淤地壩系統水力計算模型，計算結果表明淤地壩系統對水資源再分布的影響主要體現在減少地表徑流、增加入滲補給等方面。

綜上，學者們主要通過研究淤地壩建設對地表徑流的影響，進而分析其對基流量或地下水補給及占比變化的影響[15-19]，專門研究淤地壩作用下地下水位動態變化特征的較少。因此，本文以黃土丘陵區典型淤地壩刀勞慶骨干壩為例，選取淤地壩下游和淤地壩壩區兩個觀測點，通過原位監測淤地壩作用下的地下水位變化，揭示地下水位時空分布格局和動態演化特征，以期為淤地壩壩區地下水資源保護與高效利用提供參考，也為淤地壩生態治理效益評價提供科學依據。

1 研究區概況及數據來源

1.1 研究區概況

西柳溝是內蒙古“十大孔兌”之一，發源于鄂爾多斯市東勝區柴登鎮宗兌村張家山，向北注入黃河。西柳溝河道全長106.5 km，流域面積1 356 km2。流域多年平均降水量305.9 mm，主要集中在6—9月（約占年降水量的80.4%），以7月、8月降水最為集中（分別約占年降水量的28.4%、27.8%）。流域年均水面蒸發量2 200 mm，干旱指數>7，屬典型的干旱大陸性季風氣候區。西柳溝淤地壩建設開始于2000年，至2016年有淤地壩106座，大規模的淤地壩建設對地下水位的影響已顯現。

刀勞慶骨干壩位于西柳溝流域哈他土溝右岸支溝刀勞慶，地貌地形屬典型的黃土丘陵溝壑區。刀勞慶壩址以上流域面積3.27 km2，總庫容78.3萬m3，最大壩高14 m，壩長350 m。溝道呈寬淺式，溝底比降1.2%，溝床覆蓋第四系沖積洪積砂壤土，溝底砂礫石層為地下水的滲流通道，地下水埋深5～8 m。該壩于2006年年底建成，設計淤積年限20 a。

1.2 數據來源

降雨量數據通過安置在研究區的HOBO小型自動氣象站獲得，記錄時間間隔為10 min，監測時間為2019年6月21日—2020年9月30日。兩個地下水位觀測點分別選在刀勞慶淤地壩下游和淤地壩壩區內，分別記為J1點（高程1 327 m）和J2點（高程1 336 m）。J1和J2井深分別為7 m和10 m，觀測井井口高程超出地面高程0.5 m。J1點位于J2點下游，兩觀測點直線距離478 m，高差9 m。地下水埋深數據通過安裝在監測井中的Diver自動水位記錄儀獲得，每30 min記錄一次，監測時段為2019年6月26日—2020年9月30日。安置在淤地壩壩區的J2監測井中的水位計2019年8月15—22日堵孔，其間數據通過插補獲取。

2 地下水位時空分布和動態演化特征

2.1 年尺度地下水位動態特征

淤地壩下游地下水位監測井J1水位波動幅度較大，年最大變幅1.82 m左右，地下水埋深表現為波動增大—持續減小—波動較快增大—急劇減小，春夏季水位降低，秋冬季水位逐步回升，地下水位總體呈下降趨勢。壩區地下水位監測井J2水位變化較小，年最大變幅1.27 m左右，地下水埋深呈現減小—緩慢增大—快速減小—波動減小的變化趨勢，總體變化不大（見圖1）。

2.2 季尺度地下水位動態變化特征

春季（3—5月）隨著氣溫逐漸回升，淤積凍層在垂向上沿自上而下和自下而上2個方向同步融化，凍層中部消融具時滯性[21]，壩下游地下水位緩慢下降。夏季（6—8月）降雨、蒸發和灌溉用水量增加，壩下游地下水位下降，與春季相比，地下水位下降0.88 m。秋季（9—11月）農灌用水量減少，蒸發減弱，地下水位較夏季回升了0.46 m（見圖2）。冬季（12月至次年2月）為土壤凍結期，隨著氣溫持續降低，土壤由表層向下逐漸凍結，水分向凍結鋒面遷移，發生“凍后聚墑”現象，地下水位持續回升。

春季淤地壩壩區內地下水位亦呈小幅降低趨勢;夏季隨降雨量增加，地下水位升高;秋冬季地下水埋深分別增大到7.31 m和7.37 m，這可能與其側向排泄有關。

2.3 汛期與非汛期地下水埋深變化特征

（1）淤地壩下游地下水埋深變化特征。2019年汛期（6—9月）降雨量226.6 mm，壩下游地下水埋深由6.18 m減小為5.80 m，地下水位上升了0.38 m;2020年汛期降雨量312.7 mm，較上年增加38%，汛期地下水埋深由6.25 m減小到5.44 m，減小0.81（見圖3（a））。2020年汛期地下水位升幅大于2019年同期，但淤地壩下游平均地下水埋深由2019年的6.04 m增加到6.14 m，地下水位仍呈下降趨勢。非汛期（2019年10月到2020年5月），地下水埋深由5.84 m減小到5.34 m（見圖3（b））。

（2）淤地壩壩區地下水埋深變化特征。2019年汛期（6—9月）淤地壩壩區地下水埋深由7.74 m減小為7.31 m，水位上升0.43 m。2020年汛期壩區地下水埋深由7.51 m減小到6.67 m，水位上升0.84 m（見圖3（a））。與2019年相比，2020年汛期降雨量增加38%，而平均地下水埋深由7.48 m減小到6.93 m，水位抬升0.55 m。因此，地下水位隨著降雨量增加而抬升的同時，仍受到用水增加的影響，地下水位上升幅度較小。非汛期（2019年10月到2020年5月）降水少，壩區地下水埋深由7.34 m增大到7.46 m，增加0.12 m（見圖3（b））。

2.4 月尺度地下水埋深動態特征

（1）淤地壩下游地下水位變化特征。2019年7—8月降雨量97.2 mm，地下水埋深由5.80 m增大到6.18 m（見圖4），降雨量增加，地下水埋深增大。這一時期蒸發強烈，農灌次數增加導致地下水開采量增加，一定程度上抑制了地下水位抬升。9月蒸發和開采量減少，水位小幅回升。10月到次年3月，降水減少、氣溫降低、蒸散發減弱、地下水開采量減少，地下水位回升0.50 m。4—5月隨著天氣回暖，蒸散發增強，而降雨量僅6.4 mm，地下水得不到補充，地下水位下降0.27 m。6—7月總降雨量162 mm，地下水埋深繼續增大至6.41 m，地下水開采量增加是地下水位降低的主要原因。

（2）淤地壩壩區地下水埋深特征。2019年7月壩區地下水埋深為7.74 m，之后隨著降雨量增加，地下水埋深逐步減小，到10月減小至7.29 m。2019年11月到2020年6月地下水埋深呈緩慢增大趨勢，最大埋深7.51 m，該時期降水稀少（除2020年6月降水量84.4 mm），蒸散發作用由弱變強，蒸發是影響地下水位變動的主要原因。2020年7月，隨著降雨量持續增加，地下水埋深減小到6.86 m，比6月（7.51 m）減小0.65 m，減幅8.66%。8—9月總降雨量150.7 mm，隨著氣溫逐漸降低，蒸散發作用減弱，地下水得到補充，埋深減小至6.67 m。

3 地下水補排關系及影響因素分析

3.1 地下水補排關系分析

淤地壩下游監測井位于排泄區，潛水含水層賦存于第四系沖積洪積層中。巖性以沖積砂、砂礫石為主，主要接收大氣降水的入滲補給，其次接收農業灌溉回歸水補給，地下徑流條件好[21]，故淤地壩下游監測井的地下水位可代表其附近第四系的地下水位。淤地壩壩區監測井揭露的第四系潛水含水層分布較連續，其水位可代表淤地壩壩區和附近區域地下水位變化趨勢。

監測區內淤地壩壩區與淤地壩下游監測井地面高差為9 m，井深分別為7 m和10 m。壩區和淤地壩下游地下水位分別為1 328.72 m和1 321.21 m，壩區內地下水位比下游高7.51 m。水頭差的存在導致淤地壩內地下水補給壩下游地下水，與淤地壩下游地下水位恢復較快的事實相吻合。壩區地下水位和淤地壩下游地下水位差值非汛期（1.86 m）大于汛期（1.15 m），進一步說明壩區地下水補給壩下游地下水[22]。

3.2 地下水位變化驅動因素分析

研究區第四系潛水含水層具有埋藏淺、透水性好的特點，地下水位變化受地形地貌的控制，在不同時間尺度和空間分布上表現出不盡相同的動態特征。汛期降雨量增加38%，壩下游地下水埋深增大1.27%，壩區地下水埋深減小7.35%。在其他因素保持不變的條件下，農業灌溉用水量增加是影響壩區和下游地下水位變化的重要因素，下游受影響更明顯，這與下游農民大量種植玉米等農作物密切相關。淤地壩下游地下水位變化基本上遵循大氣降水量增加地下水位抬升的規律，但農灌用水量增幅大于降水量增幅，表現為6—8月（汛期）降雨量增大，地下水位降低，汛期過后地下水位逐步回升。因此，淤地壩下游地下水位變化是降水入滲補給、蒸發排泄、農業灌溉等多種因素綜合作用的結果[21]，降水量增加是壩區地下水位抬升的主要原因，農業灌溉用水量增加是地下水位降低的重要影響因素，而蒸發量增加在一定程度上抑制了地下水位抬升。壩區地下水位受農業灌溉的影響不及下游大。

4 結 論

（1）淤地壩下游監測井地下水位年最大變幅為1.82 m左右，地下水位呈下降趨勢。淤地壩壩區監測井地下水位變化較小，年最大變幅為1.27 m左右，地下水位呈上升趨勢。

（2）汛期淤地壩下游和淤地壩壩區地下水位分別回升0.81 m和0.84 m，但淤地壩下游地下水位仍呈下降趨勢。非汛期淤地壩下游地下水埋深由5.84 m減小到5.34 m，地下水位有所回升;壩區地下水埋深由7.34 m增大到7.46 m，地下水位略有下降。

（3）淤地壩壩區地下水位比下游地下水位高7.51 m，淤地壩壩區地下水補給壩下游地下水。兩者地下水位差值非汛期大于汛期，進一步說明壩區地下水補給壩下游地下水。

（4）淤地壩下游地下水位變化基本上遵循降水量增加地下水位抬升的規律，形成高峰（非雨季）與低谷（雨季），雨季過后水位逐步回升，且降雨補充地下水具有明顯的時滯性。淤地壩下游及壩區地下水位變化是降水入滲補給、蒸發排泄、農業灌溉等多種因素綜合作用的結果，農業灌溉是地下水位變化的重要影響因素，其對淤地壩下游地下水位的影響大于對壩區的影響。

參考文獻：

[1] LI B， LIANG Z， ZHANG J， et al. Attribution Analysis of Runoff Decline in a Semiarid Region of the Loess Plateau， China[J].Theoretical and Applied Climatology，2018，131（1-2）：845-855.

[2] HASHEMI A， GHAZARIYAN S. The Effect of Sedimentation in Rock Check Dam Reservoirs on the Peak of Flood in Small Watersheds[J].European Journal of Nuclear Medicine & Molecular Imaging， 2014，39（S）：386-387.

[3] YAZDI J. Check Dam Layout Optimization on the Stream Network for Flood Mitigation： Surrogate Modelling with Uncertainty Handling[J].Hydrological Sciences Journal， 2017，62（10）：1669-1682.

[4] 高云飛，郭玉濤，劉曉燕，等.黃河潼關以上現狀淤地壩攔沙作用研究[J].人民黃河，2014，36（7）：97-99.

[5] 高海東，賈蓮蓮，李占斌，等.基于圖論的淤地壩對徑流影響的機制[J].中國水土保持科學，2015，13（4）：1-8.

[6] CALLOW J N， SMETTEM K R. The Effect of Farm Dams and Constructed Banks on Hydrologic Connectivity and Runoff Estimation in Agricultural Landscapes [J].Environmental Modelling & Software，2009，24： 959-968.

[7] 吳凱，殷會娟，何宏謀，等.基于水文連通特征的黃土高原淤地壩系[J].應用基礎與工程科學學報，2020，28（3） ：717-726.

[8] 袁水龍，李占斌，李鵬，等.基于MIKE模型的不同淤地壩型組合情景對小流域侵蝕動力和輸沙量的影響[J].水土保持學報，2019，33（4）：30-36.

[9] 畢慈芬，鄭新民，李欣，等.黃土高原淤地壩建設對水環境的調節作用[J].人民黃河，2009，31（11）：85-86.

[10] HUANG T M，PANG Z H， EDMUNDS W M. Soil Profile Evolution Following Land-Use Change： Implications for Groundwater Quantity and Quality[J].Hydrological Processes，2013，27（8） ：1238-1252.

[11] GAO Z， ZHANG L， CHENG L， et al. Groundwater Storage trends in the Loess Plateau of China Estimated from Streamflow Records[J].Journal of Hydrology，2018，53（1）：281-290.

[12] LIANG W， BAI D， WANG F， et al. Quantifying the Impacts of Climate Change and Ecological Restoration on Streamflow Changes Based on a Budyko Hydrological Model in Chinas Loess Plateau[J].Water Resources Research， 2015，51（8）：6500-6519.

[13] 綦俊諭，蔡強國，方海燕，等.岔巴溝流域水土保持減水減沙作用[J].中國水土保持科學，2010，8（1）：28-33.