城市發展模式變化對地下水補給的影響

崔一嬌 杜旋 孫趙爽 尚子琦 王樹芳

摘 要：近年來，城市發展模式由資源消耗型轉變為智慧生態型。人類活動對地下水環境不再僅帶來“負效應”，更通過管網滲漏、人工補給、節水壓采等途徑逐漸突顯出對地下水資源的“正效應”，因此，需要從“正負效應”綜合分析城市化對地下水資源的影響。北京市大興區具有城市發展模式轉變的典型特征，以1996年、2006年、2016年為時間節點，分析該地區淺層水在城市發展過程中的時空變化特征。分析表明大興區城市化擴張明顯，建設用地的面積占比由1996年的19.80%增長至2016年的53.61%，而水位變幅小于0.3 m的區域面積從1996年的5.17 km2增長至2016年507.99 km2，與建設用地的范圍高度重合。不透水面空間的連續性和集中性，更易阻礙雨季降水對地下水的補給，致使水位變幅減弱。隨降水量的增加和再生水的推廣利用，淺層水位得到控制并回升，局部范圍甚至改變了原有的地下水流場狀態。2006—2016年淺層地下水儲量增加約0.72億m3，灌溉區和永定河沿岸的增加量占其中的83.33%。人工補給已經是除自然降水外不可忽視的地下水補給方式，并使人為活動產生的“正效應”在對地下水資源的影響中發揮主導作用。

關鍵詞：城市發展模式;地下水補給;土地利用;再生水灌溉;地下水動態特征

Abstract： In the recent 20 years， the urban development model has changed from that of resources consumption to intelligent ecology. Human activities no longer only bring “negative effects” to the groundwater environment， but also “positive effects” through pipeline systems leakage， artificial recharge， water saving measures， groundwater exploitation control and so on. Therefore， it is necessary to comprehensively analyze the impact of urbanization on groundwater resources balancing the “positive and negative effects”. Daxing District in Beijing has experienced the typical transformations of urban development models. This paper analyzes the temporal and spatial variation characteristics of shallow water in the process of urban development in Daxing District， based on the data in 1996， 2006 and 2016. The analysis shows that the urbanization of Daxing District is obviously expanding， the area proportion of construction land increased from 19.80% in 1996 to 53.61% in 2016， while the area with water level variation of less than 0.3 m increased from 5.17 km2 in 1996 to 507.99 km2 in 2016， which is highly coincident with the scope of construction land. The continuity and concentration of the impervious space are more likely to hinder the recharge of groundwater by precipitation in the rainy season， resulting in the weakening of water level variation. With the increase of precipitation and the popularization of reclaimed water in more towns， the shallow water level is under-control and rises again and in some areas，? and changes even occur to the original state of groundwater flow fields. Shallow water storage increased about 72 million m3 from 2006 to 2016， in which 83.33% was contributed by the increase in irrigation zones and along the Yongding River. In addition to natural precipitation， the groundwater recharge through human activities is a way that cannot be ignored， and the “positive effect” produced by human activities plays a leading role in the impact on groundwater resources.

Keyword： urban development model; groundwater recharge; land use/land cover; reclaimed water irrigation; dynamic characteristics of groundwater

在缺水性城市，地下水資源對城市發展具有戰略性意義，“極速城市化”和“空間蔓延式”的大擴張與資源承載能力的矛盾也更為突出（陳明星，2015）。在城市發展的過程中，人類活動通過地下水開采、土地利用方式的改變、管網系統的滲漏、灌溉等多種途徑，在地下水補徑排方式上已經極大地改變了自然的水循環模式，造成區域地下水失衡（王浩等，2016;王超等，2018）。

有研究表明，土地利用及其變化能夠高度敏感的影響水文過程，甚至在短期尺度上，比氣候變化的作用更明顯（張成鳳等，2019）。城市的擴張過程中，增加的不透水面會對降雨入滲補給能力產生干擾（宋全香等，2014;Chithra et al.，2015;Jacobson et al.，2011;夏軍等，2017a;夏軍等，2017b）。朱琳等（2013）利用WetSpass模型估算出北京平原區1982—2007年城鎮用地擴張了517 km2，而降水入滲補給量減少約3000萬m3。王新娟等（2017）通過入滲系數估算北京平原區降水補給量，2000—2013年多年平均值比1961—1980 年多年平均值減少3.2億m3。實際上人類活動對地下水資源量的影響并不只存在“負效應”。由于管網的自然損失以及主觀的人為回補也會增加地下水的補給量，即城市化進程帶來的“正效應”（黃婉彬等，2020）。但由于這些非自然補給在一些地區補給量較小同時又較難準確進行估算，所以在研究中往往被省略不計（于開寧等，2004;Hibbs et al.，2012）。于開寧（2001）在針對石家莊市的研究中發現，供排水系統的管（渠）道滲漏、綠地的過量灌溉、地下水回灌、再生水回灌等已經成為城市地下水補給的重要組成部分，在城市中心地區單位面積的地下水補給量甚至大于市區周邊。

如今，城市地下水的開采逐步得到有效的控制管理，越來越多的城市重視自然生態的可持續發展，這種城市發展模式的轉變給地下水補給帶來的“正效應”越來越不應被忽視。因此，考量城市化與地下水資源相互作用的影響亟需從人類活動“正負效應”的雙向作用綜合分析，才能更準確的評估和管理地下水資源。北京大興區在研究城市化與地下水補給關系上具有典型示范性，本次研究以該地區為研究對象，對應大興區的發展階段和區域規劃，以1996年、2006年、2016年為分析節點，分析大興區淺層地下水的動態時空演變，以及城市發展模式變化過程中對地下水資源“正負效應”的影響作用。研究成果以期為首都的規劃發展以及北方缺水城市的生態建設提供一定的參考，助力于處理好城市建設與水資源可持續發展之間的關系。

1 研究區及數據概況

1.1 地理環境條件

北京大興區位于北京南郊，總面積1030 km2，是距離城市中心區最近的郊區。北臨豐臺區、朝陽區，東臨通州區，西與房山區隔永定河相望，南與河北省接壤。該地區由永定河沖、洪積作用而成，地勢平坦，自西北向東南緩傾。區內有6條主要河道，分別為永定河、涼水河、天堂河、大龍河、小龍河和新鳳河。由西北向東南，地下水含水層由單一的潛水到多層的承壓水，巖性由以砂卵石為主逐漸過渡為中砂、細砂，富水性由強到弱（圖1）。垂向100 m左右深度普遍分布一層30～40 m的粉質黏土，100 m以淺的第四系含水層基本連續。以此為界，100 m以淺的潛水和微承壓水可劃為淺層地下水，其補給來源主要依靠大氣降水、地下水側向徑流、灌溉水回歸和河渠入滲。多年平均降水量為518 mm（1996—2016年），主要集中在夏季。

1.2 土地利用類型的變化

大興區因其地理位置的優勢成為了北京城市空間結構戰略轉型的支撐地。北京的城市化現象表現為中心區舊城改造，近郊區通過農用耕地轉化到建設用地來不斷擴張城市空間 （劉芳，2010）。北部建設發展了黃村、亦莊兩個新興衛星城，成為高技術制造業和戰略性新興產業的聚集區，南部建有北京大興國際機場，成為了京津冀的空間樞紐。在臨近北京市中心的北部呈現完全城市化的趨勢，農村地區小城鎮建設用地面積也大幅度增加。

2009—2016年期間，建設用地面積由35809.79 hm2上漲到39743.85 hm2，農用地面積由58635.74 hm2縮減至54793.42 hm2，未利用地面積基本持平（圖2）。這證實了大興區的城市化擴張主要是依靠減少農用地面積來實現的。

1.3 用水結構的調整

地下水是大興區的主要供水來源之一，在北京城市化發展加速時期，水資源剛性需求的增加使該地區地下水一直處于超采狀態，開發利用率一度可達146.7%（楊廣慶等，2014），并導致地面沉降、地下水漏斗問題突出（周超凡等，2017）。面對水資源短缺的嚴峻形勢，在“十一五”規劃后，大興區用水政策從一味的過度開采消耗轉變為重視地下水的涵養恢復。通過高效節水措施，提高灌溉水利用系數，實現用水定額管理，精確計量，有效控制了用水總量，并著力構建了“一廊、兩線、三網、四帶、五園”的水系格局（趙雙，2015）。

其中，挖掘再生水利用潛能、代替地下水開采是大興區重要的用水結構調整。2002年南紅門灌區開始投入使用再生水，灌區地下水位下降速度明顯緩解（張鵬卿等，2013）。2006年再生水灌溉經驗在多個鄉鎮得到進一步推廣（圖3a），利用量得到大幅提高，從0.2億m3使用量增長維持在1.1億m3左右（王勇等，2012;溫先高，2012），使地下水的開采量從3億m3縮減并穩定在了2億m3左右。地下水的減采有利于地下水的涵養恢復，而再生水的灌溉能夠補給影響淺層地下水。2005—2014年大興區的用水量統計結果顯示，再生水已經成為大興區主要的供水來源之一，主要用于環境，其余用于工業、農業和市政用水，在全區總用水量的占比達31%～33%（圖3b）。

1.4 數據資料

由于地表補給環境的變化直接作用于淺層地下水，故搜集了54組80 m以淺的潛水或微承壓水地下水監測井歷史數據（大興區內38組），范圍覆蓋大興區及周邊，用于分析區域地下水位變化（圖1）。

在此基礎上，結合地理位置、水文地質條件、城市發展特點選取了3處典型點進行深入的對比分析，3處典型點分別處于不同富水性分區內。海子角位于大興北部，屬于城區中心，建設用地集中，降水入滲受地表不透水面影響較多。青云店位于大興東部，屬于再生水灌溉重點地區。小店位于大興南部，近些年受大興機場建設和永定河治理影響，受城市化擴張影響最晚。時間節點對應大興區發展，以10年為間隔便于對比城市建設的變化，選擇1996年、2006年和2016年，年降水量分別為695.4 mm、447.1 mm、722.4 mm。1996年地下水開采現象仍十分嚴重，隨后逐漸得到控制。2006年降水量偏少，但再生水灌溉量開始大幅提升，補給方式產生明顯變化。2016年降水量與1996年相近，便于對比。

2 結果與分析

2.1 地下水位變化趨勢

圖4為海子角、青云店、小店3處典型點1992—2018年的地下水位變化趨勢，2006年是再生水大范圍投入使用的時間，1999年后降水量出現由枯轉豐的趨勢。結合該情況，擬合了再生水使用前、再生水使用后以及降水枯豐轉化期（2000—2018年）3個時間段的地下水位趨勢線。

1999年后，北京進入連續枯水年，在此期間大興區地下水開采量維持在3億m3～3.5億m3之間，2000年后降水量逐漸增加，2007年后降水量基本達到或超過多年平均水平。2006再生水投入后地下水開采量降至2.1億m3～2.5億m3。一般規律下，淺層水位主要受降水和農業開采的影響而變化。開采量的穩定更突顯出降雨補給的影響作用。

由圖4可以看到隨降水量的減少，地下水位均呈現持續的下降趨勢。但隨2000年降水量的逐步增加，僅有青云店一處水位和降水變化有較為一致的趨勢（圖4b）。海子角和小店仍持續下降至2012年，較1998年下降分別為10.77 m和11.95 m（圖4a，圖4c）。此外，海子角年內變幅（最小值與最大值范圍）多年處于2 m左右，明顯小于其他兩處。小店的地下水位年內變幅也出現明顯的縮減，從20世紀90年代的7 m左右減弱至3 m左右。地下水位的持續下降和變幅減小的特征說明降水對地下水的補給作用在逐漸減弱，在開采量較為穩定甚至減少的情況下，降水補給作用的減弱很可能與地表覆蓋變化有一定關系。

此外，從典型點的水位變化趨勢中可以看到，再生水的投入使用對地下水位有明顯的影響。小店監測井在2013年所在鎮投入使用再生水后，地下水位扭轉回升，相比2012年水位回升2.44 m。而青云店由于處在再生水南紅門灌溉區內，在再生水灌溉的作用下，枯水年期間水位仍能保持平穩。并且需要注意的是大興地區地下水流向自西北向東南，但2006年以后在再生水灌溉和降水量回升的影響下，青云店監測點的水位17～21 m，海子角和小店的水位在17 m以下，青云店水位相對其他兩處已成為高值，在局部區域內已改變了原有地下水流場方向。可見，在降水不足的條件下，再生水灌溉已經是該地區地下水的一種重要的補給來源，并已影響了局部的地下水流場。

2.2 地表入滲條件的影響

為分析大興區地下水位變化趨勢和年內變幅與地表覆蓋的關系，采用監督分類方法提取了建設用地、農用地、未利用地（包括休耕地）共3類土地利用類型。選用的支持向量機（SVM）方法可以自動尋找對分類有較大區分能力的支持向量，將類與類之間的間隔最大化，具有較好的推廣性和較高的分類準確率（祁亨年，2004）。 本次的樣本分類度達到1.94以上，樣本之間可分離性好，為合格樣本。水位變幅是利用區內54組淺水監測數據計算了6—9月的水位月變幅的平均值。該時段的選擇考慮了降水的時間段、補給的滯后性，并且開采量穩定，能夠反應降水入滲的補給程度。圖5為1996年、2006年和2016年的土地利用類型和地下水位變幅的空間分布結果，表1為土地分類的面積統計結果。

從圖5可以看出，月變幅與土地利用結果在空間分布上有明顯的關聯性。3個年份的水位變幅大于0.9 m的地區多為農用地，水位變幅小于0.3 m的區域多為建設用地。1996年農用地占比最大，占70.06%，而水位變幅大于0.9 m的面積為849.55 km2，水位變幅小于0.3 m的分布僅占5.17 km2。2006年，建設用地面積有所增長，農用地面積略有下降，城市擴張更多的是利用未利用地。而到2016年建設用地面積大幅增長，占比達到53.61%，對應水位變幅大于0.9 m的面積僅占8.56 km2，而小于0.3 m的區域面積達到507.99 km2，與建設用地的范圍高度重合。

盡管1996年農業開采的影響較大，2006年的降水量相對偏小，但從整體的趨勢變化上仍可以看出地表條件對降水入滲有明顯影響，建設用地表現出最小的水位變幅，并且城市發展的不透水面的增加具有面積的連續性和空間集中性，比起零散碎斑式的分布會更加不利于降水的入滲，人類活動對地下水資源有明顯的“負效應”影響。

2.3 再生水灌溉的影響

大興區典型點的地下水位趨勢顯示出再生水使用對地下水位變化的扭轉性影響。為分析再生水利用對地下水資源的影響，對比了1996年、2006年、2016年，3年12月份的同期水位（圖6），12月份的地下水位動態受降水和開采影響較小，能夠真實的反映地下水位穩定后的情況。

與1996年相比，至2006年，地下水位以下降為主（圖6a）。下降幅度集中在5～10 m，面積703.32 km2，北部城中心區和中部魏善莊開采區地下水位下降更明顯，達10 m以上，面積99.76 km2。南部永定河北岸和東部南紅門灌溉區地下水位下降幅度稍小，在5 m以內，面積213.12 km2。地下水位回升的地區面積僅有17.87 km2。按照大興含水層介質的條件，給水度取值0.117，經計算，2006年相比1996年的淺層地下水儲量累計虧損7.92億m3。

至2016年，相比1996年，20年來整體地下水位仍處于下降狀態，淺層地下水儲量累計虧損0.71億m3（圖6b）。在南紅門灌溉區和永定河沿岸地下水位的降幅減弱，下降趨勢已得到明顯控制。下降幅度集中在5～10 m，面積達585.42 km2，下降0～5 m的地區面積有187.13 km2，下降大于10 m的地區仍集中在北部和中部的城鎮地區，面積94.13 km2。而榆垡上升面積擴大至46.74 km2。

而2016年相比2006年，地下水位以回升為主（圖6c），大面積地區上升0～5 m，面積達536.68 km2。特別在南紅門灌溉區和永定河南部出境處地下水位上升可達到5 m以上，面積68.09 km2，而下降地區的面積429.29 km2。相比2006年，2016年的淺層地下水儲量增加約0.72億m3。其中南紅門灌溉區內地下水儲量增加約0.49億m3，回升面積占灌溉區的72.46%。永定河2 km范圍內的地下水儲量增加量約0.11億m3，上升面積占計算范圍的88.99%。

水位同期對比的結果可以看到，1996—2006年期間，由于降水量低于多年平均水平，加之開采等因素的影響，地下水儲量虧損明顯。至2016年，在壓采、降水量增加等因素的影響下，地下水位已開始明顯回升。其中，人工補給起到了不可忽視的作用，灌區和永定河沿岸2006—2016年的地下水增加量占到了全區總量的83.33%。人工補給已經是地下水得到涵養恢復的重要影響因素。

3 結論及建議

（1）快速的城市化進程下，城市地表環境的變化能夠影響降水對地下水的自然補給。水位變幅大于0.9 m的地區多為農用地，水位變幅小于0.3 m的地區多為建設用地。并且隨著城市化發展，建設用地的面積占比由1996年19.80%增長至2016年的53.61%，而水位變幅小于0.3 m的面積從1996年的5.17 km2增長至2016年507.99 km2，與建設用地的范圍高度重合。不透水面的連續性和空間集中性，更易阻礙雨季降水對地下水的補給，致使水位變幅減弱。

（2）積極的生態恢復工作對地下水位的連續下降起到了明顯的遏制作用，人為補給已經是除自然降水外不可忽視的地下水補給方式，并極大程度的抵消了城市化發展所帶來的負面影響。在再生水利用大范圍推廣前，2006年與1996年的潛水儲量相比虧損約7.92億m3。隨降水量的增加及2006年再生水的推廣利用，至2016年，淺層地下水儲量相比2006年增加約0.72億m3。再生水灌溉區以及永定河沿岸地區的地下水位回升明顯高于其他地區。

（3）城市化發展對大興區淺層地下水資源量的影響，已從“負效應”為主表現為“正效應”為主。由于連年降水量偏少、多年的過量開采、連續不透水面的增加等因素，造成地下水位受降水量變化的影響減弱，水位持續下降，年內變幅減小，這一階段人類活動以“負效應”影響占主。在城市發展模式轉變后，特別是再生水的推廣使用使地下水位轉降為升，并且在局部地區由于再生水的使用，改變了原有的地下水流場狀態，人類活動對地下水資源的影響開始以“正效應”為主導。人為補給已逐漸成為地下水補給的重要途徑之一，彌補了地下水的自然補給缺失。

（4）城市在發展建設時應綜合規劃，考慮自然生態系統的聯系空間，多利用透水性材料，避免不透水性地表的連續和集中。再生水的使用和人工灌溉等人為措施，從減少開采和增加補給兩方面，能夠有效涵養地下水資源，因此，城市化發展中應優化配置水資源利用，增強再生水的利用潛能和優化管理。同時，由于人為補給方式已成為淺層地下水重要的補給來源之一，應格外加強關注人為方式補給對地下水水質的影響。

參考文獻：

陳明星，2015. 城市化領域的研究進展和科學問題[J]. 地理研究，34（4）： 614-630.

黃婉彬，鄢春華，張曉楠，等，2020. 城市化對地下水水量、水質與水熱變化的影響及其對策分析[J]. 地球科學進展，35（5）：497-512.

劉芳，2010. 北京城市蔓延的特征及成因分析[D]. 北京交通大學.

祁亨年，楊建剛，方陸明，2004. 基于多類支持向量機的遙感圖像分類及其半監督式改進策略[J]. 復旦學報： 自然科學版，43（5）： 781-784.

宋全香，左其亭，楊峰，2004. 城市化建設帶來的水問題及解決措施[J]. 水資源與水工程學報，15（1）：56-58.

溫先高，2012. 北京南紅門灌區再生水灌溉對地下水環境影響研究[D]. 北京：清華大學.

王超，董少剛，劉曉波，等，2018. 城市化對呼和浩特市潛水補給影響研究[J]. 現代地質，32（3）： 574-583.

王浩，游進軍，2016. 中國水資源配置30年[J]. 水利學報，47（3）：265-271+282.

王新娟，張院，孫穎，等，2017. 人類活動對北京平原區地下水的影響[J]. 人民黃河，39（2）：77-81.

王勇，陳杰，張旭宗，2012. 南紅門再生水灌區建設管理中的經驗與對策探討[J]. 北京水務（2）： 54-57.

夏軍，石衛，王強，等，2017a. 海綿城市建設中若干水文學問題的研討[J]. 水資源保護，33（1）：1-8.

夏軍，張永勇，張印，等，2017b. 中國海綿城市建設的水問題研究與展望[J]. 人民長江，48（20）：1-5.

楊廣慶，孫麗芳，鄂喜梅，等，2014. 北京市大興區水資源及開發利用調查評價[J]. 北京水務（6）： 34-37.

于開寧，2001. 城市化對地下水補給的影響：以石家莊市為例[J]. 地球學報，22（2）： 175-178.

于開寧，婁華君，郭振中，等，2004. 城市化誘發地下水補給增量的機理分析[J]. 資源科學（2）：70-75.

周超凡，宮輝力，陳蓓蓓，等，2017. 北京地面沉降時空分布特征研究[J]. 地球信息科學學報，19（2）：205-215.

朱琳，劉暢，李小娟，等，2013. 城市擴張下的北京平原區降雨入滲補給量變化[J]. 地球科學（中國地質大學學報），38（5）：1065-1072.

張鵬卿，張平，楊再俊，等，2013. 大興區再生水利用模式與政策研究[J]. 北京水務（1）：42-45.

張成鳳，楊曉甜，劉酌希，等，2019. 氣候變化和土地利用變化對水文過程影響研究進展[J]. 華北水利水電大學學報（自然科學版），40（4）：46-50.

趙雙，2015，大興區再生水利用現狀[J]. 北京水務（5）：42-45.

CHITHRA S V，NAIR M V H，AMARNATH A，et al.，2015. Impacts of impervious surfaces on the environment[J]. International Journal of Engineering Science Invention，4（5）： 2319-6726.

HIBBS B J， SHARP J M，2012. Hydrogeological impacts of urbanization[J]. Environmental & Engineering Geoscience，18（1）： 3-24.

JACOBSON C R，2011. Identification and quantification of the hydrological impacts of imperviousness in urban catchments： a review[J]. Journal of Environmental Management，92（6）：1438-1448.