黃淮海平原地下水儲量與埋深變化特征分析

田益民，鮑振鑫，宋曉猛，莫昱晨，王國慶，劉翠善

（1.中國礦業大學資源與地球科學學院，江蘇徐州 221008；2.南京水利科學研究院水文水資源與水利工程科學國家重點實驗室，南京 210029；3.水利部應對氣候變化研究中心，南京 210029）

0 引 言

黃淮海平原是中國第二大平原，面積廣袤，地勢低平，是中國政治、經濟、文化、交通中心，也是我國主要的糧食生產區。糧食作物以小麥、玉米為主，其中小麥的面積及總產分別占全國的45%和51%以上，麥田面積為糧食作物種植面積的45%以上；玉米的面積及總產占30%和50%以上，玉米播種面積為73 333 km2。然而黃淮海平原也是我國北方缺水較嚴重的地區，人均水資源量僅有460 m3，低于國際公認的極度缺水的標準（500 m3）［1］。黃淮海平原旱災平均受災面積占全國受災面積的比例最高，達28%［2］。近年來的氣候變暖又加劇了水資源短缺等問題更加制約了該地區社會和經濟的發展［3］。正是因為其農業生產帶來的用水問題，使得該地區的地下水開采情況日益加重。據有關統計，目前黃淮海平原的地下水開采井約為200 萬眼，其多年平均的開采量超過了100 億m3［4］。由于多年的地下水超采，黃淮海平原已經成為世界上最大的“漏斗區”：包括淺層漏斗和深層漏斗在內的復合地下水漏斗，面積約為73 288 km2［5］。所以研究該地區地下水儲量與地下水埋深的時空變化特征，對合理和有效管理區域水資源，以及水資源可持續利用具有重要意義。

傳統檢測地下水的方法一般為觀測井監測，局限性比較大，觀察井監測方法受人為的影響比較大，并且這種方法難以充分反映陸地水儲量在空間上和時間上的連續性變化特征。根據陸地水儲量的分布情況，除去冰川水外，地下水約占陸地水儲量的97.5%，占據了陸地水儲量很大的比重，所以在平原區陸地水儲量的變化可以近似表征地下水儲量的變化。隨著衛星技術的發展，GRACE（Gravity Recovery and Climate Experiment）重力衛星作為一種監測陸地水儲量變化的新方法，打破了傳統地面觀測在時間和空間上的局限性，并且在國內外分析地下水與水儲量變化的應用中得到較好的效果［6-8］。

GRACE 衛星因為其精確的觀測精度以及大范圍的時空監測被許多研究者廣泛運用到各個領域。在國外的研究中，GRACE 衛星被運用到美國、尼羅河流域和北冰洋地區［9-11］。不僅如此，國內的研究者們也基于GRACE 衛星開展了非常廣泛的研究工作［12-17］，其中鐘敏等［16］利用2003年至2007年的GRACE月重力位資料對華北區域的水儲量進行了研究，發現華北地區存在一個陸地水的減少信號。蘇曉莉［18］利用2002-2010年的數據得到華北地區陸地水儲量以1.1 cm/a的速率減少。

許多國內的研究中更側重于結合水儲量和地下水的關系去研究問題［18，19］，在研究區內很少考慮到降水與蒸發的因素。在地下水數據獲取方面，之前的研究中大都采用的是地下水位觀測井的資料，缺少大范圍、高精度的地下水位數據，導致一部分地下水位結果可能與實際情況耦合不準確。而本文的研究中在數據獲取方面采用了水利部信息中心的黃淮海平原地下水埋深實測數據，空間分辨率為0.043°×0.043°，取得的結果可能更貼近于黃淮海平原的地下水變化的真實情況。在研究區劃分上，不僅將黃淮海平原從整體去分析它的概況，還將黃淮海平原分割成海河與淮河兩個平原去分析比較兩個區域的不同結果，同時運用了M-K檢驗在時間上去衡量水儲量和地下水埋深的顯著性。

1 研究區、數據來源和方法

1.1 研究區

黃淮海平原位于我國東部，分布范圍為E113.0°～121.5°，N30.0°～40.5°。北靠燕山，西傍太行山，南抵黃河下游，東臨渤海，包括京津2市與河北省的全部平原區，河南省黃河以北的平原區和江蘇省與山東省的部分地區，面積約為31 萬km2，總人口約1.33 億。黃淮海平原地勢較低，最高海拔為100 m 左右。氣候屬中緯度大陸性季風氣候，具有四季分明的特點。研究區位置分布如圖1所示。

1.2 數據來源

研究分析了黃淮海平原14年以來地下水儲量的趨勢變化，采用的是基于GRACE 衛星2003年1月至2016年12月的月重力位資料，GRACE 衛星是美國國家航空航天局（NASA）跟德國航空中心的合作項目，是觀測地球重力場變化的衛星。通過其重力場的變化，推測出水儲量的變化。GRACE 衛星于2002年3月升空后開始將重力數據傳回地球。數據由GSOC 收集并有三家機構共同處理、分發，分別是：JPL，UT/CSR（University of Texas/ Center for Space Research，德州大學空間研究中心）和GFZ（GeoForschungsZentrum，德國地學中心，2008年之后改名為HelmholtzCentre Potsdam—GFZ German Research Centre for Geosciences）。研究所選的資料來源（JPL：https：//grace.jpl.nasa.gov/data/get-data/jpl_global_mascons/）精度為1°×1°的逐月數據，對重力衛星中缺失的數據采用相鄰平均法進行插值處理。

降水數據采用了中國氣象數據網（http：//data.cma.cn）的地面氣候日值數據集站點資料，選取了黃淮海平原內的45個氣象站點資料，時間范圍為2003年1月1日-2016年12月31日。

地下水埋深數據采用的是水利部信息中心根據站點觀測數據插值成的空間分辨率為0.043°×0.043°的逐月網格資料，共計200×200個網格，本研究選取了研究區域內的21 000個網格。地下水供水量、地下水資源量、農業用水量與總用水量采用了中國水資源公報、淮河水資源公報與海河水資源公報數據，時間范圍為2003-2016年。

蒸發采用Gleam v3.3 的遙感數據產品，數據一共分為a、b兩套，該產品由英國布里斯托大學地理科學學院水文氣象系發布，利用來自不同的遙感衛星觀測遙感數據得到空間分辨率為0.25°的逐月蒸發數據，選取數據時間范圍為2003年1月-2016年12月。

1.3 方法介紹

Mann-Kendall突變分析方法［20-21］是一種非參數統計檢驗方法，常用于氣候變化影響下的降水、干旱頻次趨勢檢測。應用Mann-Kendall 檢驗法可以判斷氣候序列中是否存在氣候突變，如果存在，可確定出突變發生的時間。原理如下：

對具有n個樣本量的時間序列x，構造一秩序列：

秩序列Sk是第i時刻大于j時刻數值個數的累計數。在時間序列隨機獨立的假定下，定義統計量：

式中：UF1=0，E（SK），Var（SK）是累計數SK的均值和方差，在x1，x2，…，xn相互獨立，且有相同連續分布時，它們可由下式算出：

UFi為標準正態分布，它是按時間序列x順序x1，x2，…，xn計算出的統計量序列，給定顯著性水平α，查正態分布表，若|UFi|＞Uα，則表明序列存在明顯的趨勢變化。按時間序列x逆序xn，xn-1，…，x1，再重復上述過程，同時使UBk=-UFk，k=n，n-1，…，1），UB1=0。

2 結果分析

2.1 水儲量與地下水埋深時間演變趨勢

研究將黃淮海平原劃分成海河平原與淮河平原兩個區域，運用GRACE月重力位資料估計了黃淮海地區平均地下水儲量的按年變化趨勢，得出水儲量與地下水埋深在時間上分布的規律性如圖2所示。研究顯示：

（1）在海河平原和淮河平原，地下水儲量存在明顯的下降趨勢，地下水埋深呈現明顯的增大趨勢。分析得到黃淮海平原地下水儲量以1.2 cm/a 的速率減小，其中海河平原的下降速率為1.6 cm/a，淮河平原的下降速率為0.9 cm/a；地下水埋深平均每年增大0.09 m，其中海河平原平均每年增大0.16 m，淮河平原平均每年增大0.05 m。本研究與如下研究結果具有相似的結論：蘇曉莉［18］在華北地區得出水儲量以-1.1 cm/a 速率減少，鐘敏等［16］研究得到京津冀地區的陸地水儲量下降的速率大于1 cm/a。

（2）從時間上看地下水儲量的變化與地下水埋深的變化存在高度的負相關性且海河平原的地下水儲量下降程度明顯大于淮河平原，其地下水埋深的增大程度也明顯大于淮河平原，海河平原的地下水開采的嚴重程度大于淮河平原。

分別對2003年至2016年整個黃淮海平原、海河平原與淮河平原的水儲量以及地下水埋深做M-K 檢驗，所得結果（表1）。

表1 M-K檢驗結果匯總Tab.1 Summary of mann-Kendall test results

檢驗結果顯示無論在淮河平原、海河平原還是在黃淮海平原，其地下水儲量都超過了0.01 的顯著性界限，表現出顯著下降的趨勢。對于地下水埋深，在海河平原以及整個黃淮海平原，地下水埋深處于顯著增大的趨勢，而在淮河平原，地下水埋深的增大趨勢卻只超過了0.05 的顯著性界限，可以看出海河平原的地下水開采情況較淮河平原嚴重。

對黃淮海平原、海河平原以及淮河平原在年尺度［圖3中（a）、（b）、（c）］與按月的年內尺度［圖3中（d）、（e）、（f）］上對水儲量（W）和地下水埋深（D）做相關性分析如圖3所示。

分析結果表明：黃淮海平原、海河平原以及淮河平原地下水埋深和地下水儲量都呈現出強負相關的趨勢，說明了地下水儲量的變化與地下水埋深變化在時間尺度上具有一致性。

2.2 水儲量與地下水埋深空間演變趨勢

2.3 水儲量與地下水埋深變化成因分析

2.3.1 降水的時空變化趨勢

對黃淮海平原、海河平原以及淮河平原的降水進行時間與空間上的分析如圖5所示，研究表明：在時間上，2003-2016年期間黃淮海平原降雨量總體呈現出下降的趨勢，且海河平原的降雨量明顯小于淮河平原的降雨量；在空間上，黃淮海平原中部和南部的降水年變化率都小于0呈現出減少的趨勢。

2.3.2 降水變化的影響

對黃淮海平原、海河平原以及淮河平原在年尺度（圖6）上對降雨量（P）與地下水儲量（W）和地下水埋深（D）做相關性分析如圖6所示。

研究顯示：①海河平原內，降雨量與地下水儲量的相關性很低甚至呈現出不顯著的負相關性，降水對于水儲量的促進作用幾乎完全喪失。在年尺度上，淮河平原降雨量與地下水儲量呈現出較低的正相關性，降水對水儲量的增加有微弱促進作用；②無論在黃淮海平原、海河平原和淮河平原，降水與地下水埋深都呈現出非常弱的相關性，降水對于地下水埋深的影響較弱。

降水、水儲量與地下水埋深的特征分析對降水（P）、地下水儲量（W）和地下水埋深（D）在年內尺度上觀察三者之間的關系，如（圖7）所示。

研究結果顯示，地下水儲量的低谷大概集中在6月份左右，分析原因：①人為因素：在農業需求上，因為黃淮海平原作為九大農業區，農作物以冬小麥和夏玉米為主，夏玉米的主要灌溉期為5-8月［23］。所以5-6月灌溉的需求量開始增大，水儲量和地下水位開始下降；②氣候因素：處于3、4、5月份時，由降水減去蒸發的年內分配［如圖8（b）］所示，此時降水減去蒸發為負值，水儲量呈現出下降的趨勢同時蒸散發能力E比較大，而降雨P不明顯，所以這段時間內水儲量在6月份出現低谷。

地下水儲量的峰值要滯后降雨量的峰值大概2～3 個月左右，這是因為在年際尺度上蓄水變量可近似等于0；當時間處于4、5、6月份時，蒸散發能力比較大，而降雨不明顯，所以這段時間，蒸發和降水都占據主要因素，水儲量和降雨量都呈現出增加的趨勢；當時間處于7、8月份時，降水比較多，蒸發也較多［如圖8（a）］所示，所以水儲量會呈現增加的趨勢，降雨量開始逐漸下降，同時降水減去蒸發出現較大值［如圖8（b）］所示，所以此時水儲量剛好處于峰值；當時間處于1、2、10、11、12月份時，此時蒸發開始減小，而同樣降水也開始減小，降水減去蒸發的變化趨于平緩［如圖8（b）］，此時水儲量開始逐步趨于平緩。地下水埋深與水儲量的規律剛好相反呈現高度負相關，與地下水埋深具有“記憶性”的結論一致［24］。

2.4 人類活動的影響

2.4.1 地下水超采影響

不過，前述論證并不全面。畢竟，修改后的民事訴訟法擴展了法律監督的范圍，也豐富了法律監督的形式，除去啟動再審程序之外，還有權以檢察建議形式對審判監督程序以外的其他審判程序中審判人員的違法行為進行監督。再審程序之外的法律監督會影響審判機關在民事訴訟中的獨立地位嗎？這個問題在之前以抗訴制度存廢為焦點的討論中一直被忽略，現在已經躍然其上，亟待思考解答。

由水資源公報2003-2016年農業用水量、總用水量、地下水供水量與地下水資源量數據分析得到，淮河區年平均總用水量為600.4 億m3，農業用水量為405.9 億m3占總用水量的67.6%；海河區多年平均總用水量為373.1 億m3，農業用水量為250.6 億m3占總用水量的67.2%，農業生產用水量占據了很大的比重。淮河區多年平均地下水供水量為167.6 億m3，地下水資源量為412.9 億m3，海河區多年平均地下水供水量為235 億m3大于地下水資源量的232.1 億m3，海河平原明顯出現地下水超采現象，地下水位下降更明顯。李玲［25］通過建立華北平原地下水流區域數值模型方法反演估算得出華北平原2002-2008年的地下水年平均開采量為249.2億m3，黃淮海平原面積約為31 萬km2，得到地下水位下降幅度約為0.08 m/a，由于本本研究的時間尺度為2003-2016年得到的結果為黃淮海平原地下水位下降幅度約為0.09 m/a，考慮到之后農業生產的增加帶來的地下水水位下降程度的增加，故本研究具有可靠性。

2.4.2 小麥月灌溉分布影響

魏童彤［26］分析了華北地區1959-2018年小麥在月尺度灌溉需水量的變化趨勢，得出小麥4月灌溉需水量多年平均值為84.73 mm，5月灌溉需水量多年平均值為144.56 mm，6月灌溉需水量多年平均值為30.83 mm，4、5、6月為需水關鍵期，灌溉需水量約占全生育期灌溉需水量的79%，所以黃淮海平原4、5、6月農業灌溉用水量較大，地下水埋深開始增大，地下水儲量開始趨于低谷。

3 結 論

本文利用GRECE 重力衛星資料揭示出黃淮海平原地區地下水儲量的時空演變規律，同時利用降水量、地下水埋深等相關數據對整個黃淮海平原進行地下水開采情況的成因分析，得到的結論如下。

（1）黃淮海平原的地下水埋深在2003-2016年期間呈現增大趨勢，在淮河平原地下水埋深變化率分布不均勻，變化趨勢不明顯。

（2）在黃淮海平原中部以及西北部，即保定市、石家莊市與邢臺市地區，地下水水位和水儲量的下降趨勢最顯著。

（3）黃淮海平原的地下水水位和地下水儲量的下降程度大致從海河平原—淮河平原呈現遞減的趨勢，降雨量從海河平原—淮河平原大致呈現出遞增的趨勢，降雨量較少的區域地下水埋深的增大程度相對較高，地下水開采程度越顯著，水儲量下降的趨勢越明顯。

研究結果表明，黃淮海平原的地下水儲量與地下水埋深形勢不容樂觀，特別是在海河平原，地下水開采程度日益劇烈，與其農業灌溉用水巨大的需求量有著密切的關系。如何在保持其農業生產的同時修復這個巨型“漏斗區”已成為黃淮海平原地區未來必須面對的難題。□