變化環境對華北平原地下水可持續利用的影響研究

賴冬蓉， 陳益平， 秦歡歡,， 高 柏， 孫占學,

(1.東華理工大學 核資源與環境國家重點實驗室， 江西 南昌 330013； 2.東華理工大學 水資源與環境工程學院， 江西 南昌 330013)

1 研究背景

以氣候變化和人類活動為代表的變化環境對大多數國家的陸面水文循環過程和水資源利用產生了不可忽視的影響[1-5]，針對變化環境下的水資源可持續利用問題也成為水科學領域的研究熱點之一[2]。氣候變化突出表現為全球氣溫升高，據聯合國政府間氣候變化專門委員會(Intergovernmental Panel on Climate Change, IPCC)的預測，全球平均氣溫到2081-2100年將上升0.3～4.8 ℃[6]，將會對全球陸面水文循環過程產生嚴重影響，導致區域水資源總量和時空分布的劇烈變化，從而加劇區域水資源短缺矛盾[7]，深刻影響區域的社會經濟和生態環境[1]。隨著社會經濟的發展和科學技術的進步，人類對水資源的需求和利用程度越來越大，導致人類活動(農業灌溉、引水工程、土地利用、修建水庫大壩等)[8]對陸面水文循環過程及水量和水質的時空分布等產生了顯著影響[9-11]，這成為加劇區域水資源短缺矛盾的重要驅動因素。因此，針對變化環境下區域水資源利用進行研究，探討氣候變化和人類活動的影響，對于區域水資源短缺矛盾的解決、社會經濟的發展及資源環境的保護都具有重要意義。

大尺度分布式水文模型可以模擬陸面水文循環全過程，在區域水資源可持續利用和管理中發揮了重要作用，許多學者采用這種模型研究水資源利用和管理[12-15]。MIKE SHE模型是丹麥DHI公司開發的基于物理過程的、確定性的大尺度分布式水文模型軟件，它涵蓋了水文循環主要過程及它們之間的相互作用，在水文、生態、環境等領域得到了廣泛應用[16-17]。2011年，Sultana等[18]以加拿大安大略南部Spencer Creek流域為研究區，建立了MIKE SHE/MIKE 11耦合的水文模型，模擬氣候變化對該流域水文過程的影響，結果表明，氣候變化會導致年均河流流量增加約10%～25%，耦合的MIKE SHE/MIKE 11模型能有效地體現氣候變化對復雜流域尺度水文過程的綜合影響。因此，采用MIKE SHE模型對變化環境下區域水資源利用與管理進行研究是適用的。

華北平原是我國水資源壓力最大的區域之一，人均水資源占有量僅有全國平均值的23%[19]，70%的用水量來自地下水[20-21]，淺層和深層含水層地下水開發利用率分別為112%和139%[22]，地下水超采給華北平原帶來了諸多生態環境問題，如地面沉降、土壤鹽漬化、地下水降落漏斗擴大、含水層疏干等。氣候變化和人類活動是造成華北平原水資源供需矛盾最主要的兩大因素，而且隨著全球變暖的持續發展及地下水開采、引水工程、農業灌溉等人類活動的持續進行，華北平原地下水可持續利用形勢將更加嚴峻。

本文在文獻[19]和[23]建立、校準的MIKE SHE模型的基礎上，采用情景分析和數值模擬方法，考慮氣候變化(A1B CO2排放情景下代表濕潤、正常和干旱的3種大氣環流模型)和人類活動(農業節水和南水北調工程)對華北平原地下水利用的影響，從地下水水位、水均衡和含水層儲量等方面進行定量分析與比較，由此探討變化環境下適合于華北平原社會經濟和水資源環境可持續發展的策略和方案，為華北平原地下水可持續利用提供科學依據和建議。

2 資料來源與研究方法

2.1 研究區概況

位于我國黃河以北、燕山以南、太行山以東的華北平原(34°48′～40°30′N，113°10′～119°26′E)是我國政治、經濟和文化中心，該平原地區人口密集，大中城市眾多，總面積為13.73×104km2，人口密度為1 456人/km2。華北平原屬于大陸性半干旱氣候，年平均降水量為554 mm，年平均氣溫為12.5 ℃，年平均蒸發量約為1 550 mm，降水量年內分布不均(汛期降水量占全年的75%以上)、年際變化大，是我國水資源供需壓力最大的區域之一[19,23-24]。華北平原總供水量的70%來自地下水，農業用水量占總用水量的85%。該地區典型農作物輪作方式為冬小麥和夏玉米，這兩種作物約占總耕地面積的80%和糧食產量的90%，冬小麥生長季降水量僅占全年降水量的20%～30%。南水北調工程(South-to-North Water Diversion Project, SNWDP)通過跨流域的水資源合理配置，可以大大緩解我國北方水資源嚴重短缺狀況，促進南北方經濟、社會與人口、資源、環境等多方面的協調發展，該工程分為東、中、西3條調水線路。

調查表明，華北平原天然地下水量為227.4×108m3/a，淺層和深層地下水可開采量分別為168.3×108和24.2×108m3/a。由于地下水開采布局不合理，華北平原深層地下水水位低于海平面的范圍已達76 732 km2，占平原區總面積的56%。地下水超采嚴重、水資源綜合利用率低、浪費和污染嚴重、生態環境惡化及干旱頻率增加等生態環境問題已成為華北平原水資源利用過程中比較突出的問題。華北平原地理位置及地下水水位降深分布見圖1，其數據來源于1993-2006年該地區的165口觀測井，圖1中還標出了本研究選用的6口觀測井，各觀測井信息見表1。

圖1 華北平原地理位置及地下水水位降深分布

表1 本文選用的研究區6口觀測井信息

2.2 華北平原MIKE SHE模型概述

MIKE SHE是一個確定性、分布式、基于物理過程的水文建模系統，涵蓋了水文循環的主要過程[25-27]，采用有限差分法求解描述地表水流(二維圣維南(Saint-Venant)方程)、河道水流(運動學路徑求解器，kinematic routing solver)、非飽和水流(兩層水平衡法)和飽和水流(三維布西涅斯克(Boussinesq)方程)過程的偏微分方程，同時利用解析解對蒸散量和截留量進行模擬。

本文的MIKE SHE模型是基于文獻[19]和[23]所建立和校準的模型而進行地表水與地下水耦合模擬的模型。氣候數據為熱帶降雨測量任務(tropical rainfall measuring mission, TRMM)的3B42日產品，分辨率為0.25°；參考ET數據使用19個氣象站數據生成，并使用Penman-Monteith方程進行計算。由于自2000年以來，華北平原的其他河流基本都已經干涸，因此模型中考慮了5條主要河流(灤河、永定河、滹沱河、漳河和衛河)。華北平原MIKE SHE模型中土壤分為3種不同參數類型(壤土、粉砂壤土和砂土)，飽和含水層被概念化為3層(一層淺層含水層和兩層深層含水層)，用于描述多層、非均質和各向異性的含水層系統[19,24]。模型中采用水文地質單元進行參數賦值，共劃分為27個單元。經過校準的水平滲透系數取值范圍為4.91×10-5～5.59×10-5m/s，給水度取值范圍為0.06～0.19。利用226口觀測井的地下水水位資料，對模型進行了校準(2000-2005年)和驗證(2006-2008年)，模型中地下水流模擬的時間步長為1 d，河流模擬的時間步長為5 min，文獻[19]中模擬期末結果作為本文模型模擬的初始條件。

2.3 華北平原地下水利用情景設計

本文考慮氣候變化和人類活動兩類變化環境的影響。氣候變化情景設計是在A1B CO2排放情景下基于IPCC第四次評估報告的預測，選取了代表濕潤(UKMO_HADCM3)、正常(CSIRO_MK3)和干旱(CNRM_ CM3)氣候條件的3種大氣環流模型(global circulation model, GCM)對華北平原未來氣候變化情況進行預測。人類活動情景設計考慮了農業節水措施和南水北調工程供水兩類人類活動。在此基礎上，設計了4大類情景，包括現狀保持型情景SQM(statusquomaintenance)、農業節水型情景AWS(agricultural water saving)、南水北調工程供水型情景SNWDP和綜合利用型情景CU(comprehensive utilization)。這4大類地下水利用情景中每個情景均增加了濕潤(humid climate, HC)、正常(normal climate, NC)和干旱(arid climate, AC)3個氣候狀況的子情景，即情景SQM、AWS、SNWDP和CU均考慮了氣候變化影響。這4大類情景中考慮氣候變化情景名稱分別為SQM-HC、SQM-NC和SQM-AC，AWS-HC、AWS-NC和AWS-AC，SNWDP-HC、SNWDP-NC和SNWDP-AC及CU-HC、CU-NC和CU-AC，具體情景設計見表2。模擬年份為2019-2028年，以經過校準的MIKE SHE模型2018年的模擬結果作為本研究的初始值，其他設置保持不變。

表2 華北平原地下水利用情景設計

3 結果與分析

3.1 地下水水位變化分析

地下水水位是華北平原MIKE SHE模型中非常重要的變量，是地下水含水層儲量變化的直接表現。本研究選取了6口觀測井(相關信息見圖1和表1)來展示不同地下水利用情景下華北平原地下水水位的動態變化，該6口觀測井的監測深度不同，代表了不同深度的含水層類型。模擬期2019-2028年4種情景下華北平原6口觀測井地下水水位變化曲線及年均變幅分別見圖2和表3。

由圖2和表3可知，除情景CU外，其余3個情景下華北平原6口觀測井地下水水位整體上均呈下降趨勢。情景SQM、AWS和SNWDP中觀測井地下水水位的年均下降幅度分別為0.307～1.841、0.239～1.717和0.083～1.292 m/a，下降幅度以情景SNWDP最小，情景SQM最大，情景AWS居中。相對于情景SQM，情景AWS和SNWDP中觀測井地下水水位有一定的恢復，說明農業節水措施和南水北調工程均能對地下水水位的恢復起到一定的作用，且南水北調工程的作用要大于農業節水措施的作用。情景CU中觀測井地下水水位年均上升幅度為0.043～0.437 m/a，相比于其他3個情景，地下水水位從下降扭轉為上升，說明綜合各種節水和引水措施對研究區地下水水位的恢復有較大作用。

圖2 2019-2028年4種情景下華北平原6口觀測井地下水水位變化曲線

表3 2019-2028年4種情景下華北平原6口觀測井地下水水位年均變幅 m/a

3.2 地下水均衡分析

華北平原地下水系統補給量與排泄量之間的大小關系決定了地下水利用狀況，補給量小于排泄量(負均衡)意味著含水層儲量的消耗，導致地下水水位下降；補給量大于排泄量(正均衡)意味著含水層儲量的增加，導致地下水水位回升。圖3為模擬期2019-2028年4種情景下華北平原地下水系統年均補給量、排泄量和含水層儲量，由圖3可以看出，總體上情景SQM、AWS和SNWDP下華北平原地下水系統存在負均衡，地下水含水層儲量消耗量分別為119、108和81 mm/a，情景SNWDP消耗量最小，情景AWS次之，情景SQM最大；情景CU下華北平原地下水系統存在正均衡，地下水含水層儲量增加9 mm/a。

圖3 2019-2028年4種情景下華北平原地下水系統水均衡情況圖

表4為模擬期2019-2028年4種情景下華北平原主要水均衡項的平均值，其中補給項包括降水量、總灌溉量、邊界流入量和南水北調引水量，排泄項包括實際蒸散發量、總抽水量和邊界流出量。

由表4可看出，降水量為華北平原的主要補給來源，4種情景下占總補給量的比例分別為73.41%、77.27%、68.69%和74.42%；實際蒸散發量是華北平原的主要排泄項，4種情景下占總排泄量的比例分別為66.16%、69.09%、66.44%和76.45%；總抽水量是華北平原的次要排泄項，4種情景下占總排泄量的比例分別為33.37%、30.30%、32.65%和22.30%。由于農業用水量占華北平原總用水量的比例在60%以上，而實際蒸散發量是農業用水量的主要消耗部分，因此，實際蒸散發量是華北平原主要的排泄項，控制實際蒸散發量就可以控制華北平原地下水系統的排泄項，也就可以對華北平原地下水利用進行有效管理。作為中國主要的糧食產量基地之一，農業是華北平原非常重要的用水部門，而作物生長過程中的蒸散發是華北平原的主要蒸散發項，因此，對于農業及農作物進行管理是控制華北平原實際蒸散發量的重要方法。

表4 2019-2028年4種情景下華北平原主要水均衡項的平均值 mm/a

3.3 含水層儲量變化分析

含水層儲量會隨著地下水利用而不斷發生變化，負均衡會導致含水層儲量下降，正均衡則會使含水層儲量上升。圖4為考慮與未考慮氣候變化和人類活動影響的4類情景下華北平原飽和含水層儲量變化曲線。

由圖4可以看出，對于4類情景來說，人類活動(農業節水措施和南水北調工程)影響下華北平原含水層儲量有不同程度的恢復，其中南水北調工程的作用要大于農業節水措施的作用，而將它們結合起來的作用將進一步加大。具體來說，相對于情景SQM，模擬期末(2028年12月31日)情景AWS、SNWDP和CU下華北平原飽和含水層儲量分別增加0.168、0.558和1.433 m。

圖4 考慮與未考慮氣候變化和人類活動影響的4類情景下華北平原飽和含水層儲量變化曲線

氣候變化對于4類情景的影響相類似，代表濕潤氣候的情景(情景SQM-HC、AWS-HC、SNWDP-HC和CU-HC)下華北平原含水層儲量有較大恢復；代表干旱氣候的情景(情景SQM-AC、AWS-AC、SNWDP-AC和CU-AC)下華北平原含水層儲量則出現較大下降；代表正常氣候的情景(情景SQM-NC、AWS-NC、SNWDP-NC和CU-NC)所起的作用則居前兩者之中。具體來說，到模擬期末(2028年12月31日)，相對于情景SQM，情景SQM-HC、AWS-HC、SNWDP-HC和CU-HC下華北平原飽和含水層儲量分別增加0.972、1.129、1.505和2.239 m，情景SQM-AC、AWS-AC、SNWDP-AC和CU-AC下華北平原飽和含水層儲量分別減少0.940、0.736、0.372和0.607 m，情景SQM-NC、AWS-NC、SNWDP-NC和CU-NC下華北平原飽和含水層儲量分別增加0.119、0.285、0.677和1.540 m。

圖5為模擬期2019-2028年氣候變化對4種情景下華北平原地下水飽和含水層儲量造成影響的不確定性范圍。

圖5 2019-2028年氣候變化對4種情景下華北平原飽和含水層儲量造成影響的不確定性范圍

由圖5可以看出，4種情景下氣候變化均給華北平原飽和含水層儲量帶來顯著影響，由此造成飽和含水層儲量變化的不確定性也不可忽視。結合圖4和5可知，氣候變化對華北平原地下水利用的影響與人類活動(農業節水措施和南水北調工程)的影響具有相同的量級，在重視人類活動影響的同時，也必須考慮氣候變化產生的影響。

3.4 地下水影響因素分析

人類活動和氣候變化是華北平原地下水利用過程中非常重要的兩類因素。一方面，地下水開采即為人類活動對地下含水層系統的影響，過度開采地下水會導致地下水水位下降、含水層儲量消耗，影響地下水可持續利用。同時，南水北調工程和農業節水是從引水和節水角度影響地下水利用。本文研究表明，南水北調工程對華北平原地下水利用的影響大于農業節水的影響。要徹底解決華北平原地下水含水層儲量消耗問題，保證地下水可持續利用，必須將農業節水措施與南水北調工程相結合，采取綜合措施，才能起到較好效果。

另一方面，氣候變化對華北平原地下水利用會產生顯著影響，其對飽和含水層儲量影響的量級與人類活動的影響相當。華北平原地下水水位、水均衡和飽和含水層儲量變化方向和大小與氣候的干濕程度正相關。氣候越濕潤，地下水利用越容易維持可持續發展狀態；反之，氣候越干旱，地下水利用越容易出現不可持續狀態。同時，華北平原面臨著人類活動和氣候變化雙重因素影響，對于地下水利用將起到“放大器”的作用，可以放大2～3倍[28]，即氣候變化導致降水量減少，人類活動的加劇將加大地下水開采力度，這將對地下水造成非常大的間接影響。

4 結 論

針對人類活動和氣候變化等變化環境下華北平原地下水可持續利用問題，采用MIKE SHE模型，設定4大類情景(每類情景中均設定了3類氣候變化情景)，模擬農業節水措施、南水北調工程和氣候變化對2019-2028年華北平原地下水利用的影響，得到如下結論：

(1)農業節水措施和南水北調工程等人類活動對華北平原地下水水位、水均衡和含水層儲量有較大影響，它們均能使地下水水位和含水層儲量有一定程度的恢復，但南水北調工程的影響大于農業節水措施的影響。相對于情景SQM，模擬期末(2028年12月31日)農業節水型情景AWS、南水北調工程供水型情景SNWDP和綜合利用型情景CU下華北平原飽和含水層儲量分別增加0.168、0.558和1.433 m。

(2)氣候變化對華北平原地下水水位、水均衡和含水層儲量影響的量級與人類活動的影響相當，且這三者變化方向和大小均與未來氣候條件的干濕程度正相關。到模擬期末(2028年12月31日)，相對于情景SQM，濕潤氣候情景下華北平原飽和含水層儲量增加0.972～2.239 m，干旱氣候情景下華北平原飽和含水層儲量減少0.372～0.940 m，而正常氣候情景下華北平原飽和含水層儲量增加0.119～1.540 m。

(3)僅靠單一方法無法徹底解決華北平原地下水不可持續利用問題，必須將農業節水措施和南水北調工程相結合，采取綜合措施，才能保證華北平原地下水的可持續利用。人類活動和氣候變化往往共同起著“放大器”作用，因此在未來制定相關政策時，必須考慮氣候變化的影響。