河北省地下水壓采工程區坑塘水體動態變化分析

趙 浩，張瑞芳，齊永青，周大邁

(1. 河北農業大學資源與環境科學學院，河北 保定 071001；2. 國家北方山區農業工程技術研究中心，河北 保定 071000；3. 河北省山區農業工程技術研究中心，河北 保定 071000；4. 中國科學院遺傳與發育生物學研究所農業資源研究中心，石家莊 050021)

在北方的華北平原一帶，淡水資源嚴重匱乏，人均水資源量僅為300 m3。人們長期以來靠開采地下水維持日常的運轉，而且水資源利用率低，存在很多浪費現象，導致深層地下水的開采量已經遠遠超出了自然所能補給的能力，造成華北地區地下水位大幅下降，并形成巨大的地下水漏斗[1,2]。而坑塘作為重要的內陸水體，不僅對社會經濟發展至關重要，而且也是維持區域生態平衡與安全的必要條件[3]。因此地表的坑塘水域資源的功能性日益顯著，尤其是其對農業灌溉的積極影響，地表坑塘可以在豐水期存儲部分水資源，到了枯水期，人們可以利用其中儲存的水資源進行灌溉。所以如何提高坑塘水資源的利用效率，完善坑塘周圍水域配套水利工程措施成為當下熱門話題。

一般情況下，將坑塘定義為人工開挖或天然形成的蓄水量小于10 萬m3的土地。劉乙玄等[4]歸納總結了坑塘在生態系統中的功能以及研究現狀并提出了對坑塘生態功能的展望。王艷穎等[5]就3次降雨徑流后坑塘濕地對生源要素氮的凈化效果進行了研究，結果表明坑塘還能對湖泊水體的富營養化緩解發揮明顯的作用。由此可見，坑塘對生態環境的調節有著巨大作用。宋求明等[6]基于MODIS遙感影像用歸一化植被指數(NDVI)閾值法提取了洞庭湖水面的面積，但是MODIS數據在提取較小水體時其250 m的分辨率略顯不足。Charles Verpoorter[7,8]基于Landsat 7 ETM+數據開發了一種GWEM的水體提取算法，并建立了全球尺度的湖泊數據庫。羅玲等[9]基于Landsat衛星數據通過監督分類和歸一化水體指數(NDWI)方法提取了湖泊和水庫/坑塘的信息并系統分析了松嫩平原西部湖泊和水庫/坑塘蓄水面積的時空動態變化?？梢奓andsat衛星數據在提取水體上具有優勢。在嚴重依賴開采地下水的華北平原，坑塘作為一種非常規水體，有著很重要的調蓄作用。目前華北地區地下水限制開采的背景下，有必要較系統的分析和統計坑塘的分布和動態變化情況，使坑塘充分發揮調蓄作用和相應的灌溉能力。

本文基于Landsat系列衛星的TM/ETM+/OLI遙感影像，對2000-2014年間夏季的河北省地下水壓采工程區的坑塘水體進行了提取，并查明了地表坑塘有關數據，估算了其蓄水能力，分析了坑塘作為灌溉水源的保障能力波動情況，對此地區的水資源保護和農業灌溉提供有力的科學支持。

1 研究區概況及數據來源

1.1 研究概況

本文研究區為河北省地下水壓采工程區項目分區，所在位置位于河北省東部環渤海低平原區，即渤海西部海拔低于60 m的低平原區，屬于黃淮海平原的一部分，主要包括滄州、衡水、邢臺、邯鄲、保定、廊坊等市的44個縣市，經度范圍大致為114°35′E ～117°52′E，緯度范圍大致為36°8′N ～39°24′N，全區面積約1.964 萬km2。本地區屬于溫帶大陸性季風氣候區，冬春寒冷干燥，夏季炎熱多雨。多年平均氣溫 10.0～14.2 ℃，平均年降水量400～500 mm，降水時空分布不均，主要集中在夏季 6-8月，春季和秋冬季較為干旱，年均蒸發量900～1 400 mm，干燥度1.5左右[10,11]。研究區基本處在海河流域內，區內主要河流有大清河、子牙河、滏陽河，主要水庫、湖泊有白洋淀、衡水湖、大浪淀等(見圖1)。

圖1 研究區范圍Fig.1 The scope of study area

1.2 數據資料的獲取

在地理空間數據云網站(http:∥www.gscloud.cn/)上獲取2000年landsat-5 TM衛星影像、 2007年landsat-7 ETM+衛星影像(首先進行去條帶處理)和2014年的Landsat-8 OLI衛星影像。為了不影響水體提取結果，研究區內所有的影像的云量均小于5%。

氣象資料由中國氣象數據網(http://data.cma.cn/site/index.html)提供。為2000-2013年間七八月份研究區范圍內20個氣象臺站的月值氣象資料，包括平均氣溫和累月降水量。

2 研究方法

和其他水體，如湖泊或者河流的不同之處在于，坑塘水面較小，而且形狀不規則，其與周圍植被和潮濕的土壤緊密相連容易混淆，所以提取水體時需要極為精細。Landsat衛星具有長期數據積累的優勢,方便進行多年的動態監測，同時Landsat TM/ETM+/OLI傳感器都包含植被和土壤等背景因素且能夠很好地區分出來，如反映水體信息的短波紅外波段(0.85～1.7 μm)。因此，本研究中選用Landsat衛星影像作為數據源，并對下載的多幅圖像進行篩選對比，選取了圖像質量相對較好的幾景圖像。

利用ENVI 5.1軟件對下載的遙感影像進行輻射定標、大氣校正，再進行拼接、裁剪等預處理。然后通過公式推算和4、3、2波段標準假彩色合成圖像、歸一化水體指數(NDWI)與改進的歸一化差異水體指數(MNDWI)提取水體的效果進行比較(見圖2)，確定最適合在本研究中提取坑塘水體的方法，即福州大學徐涵秋的改進的歸一化差異水體指數方法[12-15]。

歸一化差異水體指數(NDWI)、改進的歸一化差異水體指數(MNDWI)的計算公式如下：

(1)

圖2 MNDWI、NDWI與RGB(近紅、中紅和藍波段)組合對比(從左到右依次排列)0值代表水體Fig.2 Comparison of MNDWI, NDWI and combination of red, green and blue (Near infrared, mid infrared and blue band) 0 value represents water

(2)

式中：BMIR為中紅外波段；BNIR為近紅外波段；BGreen為綠光波段。

MNDWI原理是利用土壤的反射率從近紅外波段到中紅外波段(0.76～1.75 μm)驟然轉強，而水體在該波長范圍內的反射率降低的規律，采用中紅外波段替換近紅外波段。相較于NDWI得出的水體指數將會增大，可以較精確地區分出細小水體和土壤地物。從而得出結論：MNDWI對細小水體或淺水體的反應程度比NDWI強。

然后在ArcGIS平臺中進行二值化處理，把背景和水體區分出來。最后去除單個像元，即面積在900 m2以下的，對坑塘蓄水的研究沒有太大意義的小坑塘。此外，在遙感圖像中也提取出一些線狀水體如河流、溝渠等。由于河道中的水是流動的，沒有蓄水能力，所以本研究不將其作為研究對象。研究區西部，黃驊、海興存在部分近海及海岸濕地，由于土地表面鹽漬化嚴重，且土壤水溶性鹽含量較高，因此其地表坑塘水資源并不能適于農業灌溉[16,17]，因此提取出的水體不計算在坑塘面積之內。此外，面積較大的湖泊、水庫，如白洋淀、衡水湖等水體也不再研究范圍之內。提取坑塘水體的過程中為避免出現較大誤差，在ArcGIS 平臺下結合遙感影像和Google Earth采用目視解譯的方法，對水體信息進行佐證。最后將研究區范圍裁剪出來，得出研究區內的坑塘水體圖(見圖3)。

圖3 2000、2007、2014年夏季河北省地下水壓采工程區坑塘分布情況Fig.3 Distribution of pit-ponds in 2000, 2007, 2014 in groundwater restricted exploitation area in Hebei province

3 結果與分析

(1)坑塘分布情況。從圖3可以看出從坑塘的提取結果圖可以看出，研究區范圍內坑塘分布呈北多南少、東多西少的格局。滄州市、廊坊市、衡水市境內坑塘分布最密集，這與區域氣候、地理條件有著很大關系。滄州環渤海地區屬于暖溫帶半干旱季風氣候區，較之研究區西部的大陸性季風氣候，氣候比較濕潤，降水條件較好。研究區北部小白河、大清河、中亭河、子牙河、等流經此區域，從而地表水資源較為豐富，使得大部分坑塘蓄水充足。

(2)坑塘數量、面積變化情況。河北省地下水壓采工程區內2000年的坑塘數量最多，為13937個；到了2007年坑塘數量驟降，減少到9 313個，萎縮現象明顯；2014年坑塘數量減少到8 423個，變化不顯著?？犹列钏娣e在2000年最大，達到12 644.1 hm2，在2000-2007年間坑塘蓄水總面積呈現下降趨勢，減少到9 349.0 hm2。而到2014年坑塘面積縮減到4 835.2 hm2，萎縮現象明顯。2000、2007和2014年坑塘數量和面積呈明顯下降趨勢，就坑塘數量來說2000-2007年減少明顯，減少了4 624個，而從2007-2014年的7 a間坑塘個數減少了890個。就坑塘的水域面積來說，2014年的坑塘總面積不足2000年的40%。盡管2007-2014年間坑塘數量減少不明顯，但是坑塘蓄水總面積明顯減少。這與區域氣候條件及人類開發與利用活動密切相關[18](見表1)。

表1 坑塘個數及面積統計Tab.1 Statistics of pit-ponds number and area

(3)不同面積等級坑塘變化分析。根據研究區坑塘特點，將坑塘按蓄水面積分為3個等級：一級(<0.5 hm2)、二級(0.5～1 hm2)、三級(>1 hm2)。根據圖4可以看出河北省地下水壓采工程區的坑塘數量以面積小于0.5 hm2的居多，面積在0.5～1 hm2之間的坑塘和面積大于1 hm2的坑塘數量相差不大。總體來說，從2000-2014年不同等級的坑塘數量均呈下降趨勢。其中2000-2014年間小于0.5 hm2的一級坑塘數量分別為8 847、4 525、6 604個呈先減少后增加的趨勢。2000-2007年坑塘數量減少4 525個，減少約50%，而從2007-2014年間坑塘數量增加了2 079個。面積在0.5～1.0 hm2的坑塘在2000-2007年間數量小幅增加了200個左右，在2007-2014年間面積減少一半。面積大于1 hm2的坑塘個數呈現減少趨勢，在2000-2007年間數量小幅減少了500個左右，在2007-2014年間則銳減了近一半。對于較小的一級坑塘，它們受降水等氣候條件影響較大，當某一年份夏季降水量少的情況下，小坑塘內的水體伴隨著較大的蒸發量其水域面積會大幅縮減甚至干涸。

圖4 不同面積等級的坑塘數量變化統計Fig.4 Capacity of pit-ponds water storage during 2000 to 2014

(4)估算坑塘的蓄水量。圖5是經過對南皮縣坑塘實測面積和蓄水量得到的二者之間的線性回歸關系。

圖5 坑塘面積與蓄水量的關系Fig.5 Relationship between area and water storage

根據這個關系可以提取出公式(3)：

Y=4.03x+0.18

(3)

估算出河北省地下水壓采工程區內坑塘的蓄水量見表2。

表2 2000-2014年坑塘蓄水能力Tab.2 Capacity of pit-ponds water storage during 2000 to 2014

由于坑塘蓄水量是由坑塘蓄水面積根據蓄水公式估算出來的，所以坑塘蓄水量與坑塘蓄水面積呈正相關。在2000-2014年的15 a間也是坑塘的蓄水能力也呈逐年減少的趨勢，尤其以2007-2014年間蓄水能力減少的較為迅速。

4 結 語

河北省地下水壓采工程區內的坑塘作為一種非常規水體，其季節性差異非常明顯，它可以通過收集利用汛期的水資源用來排澇和灌溉實現其調蓄作用。研究區在8月份正處于一年之中降水最多的時期，并且空氣濕度大，坑塘基本可以儲滿水。到了秋季，研究區內的玉米作物收獲，轉而種植冬小麥的時期，可以盡可能多的利用坑塘夏季儲存的水來灌溉，這樣不僅充分利用坑塘的蓄水能力，同時還可以減少區內地下水的開采，如果對查明的坑塘周圍建立配套的工程措施，對緩解河北省的淡水資源短缺有重大意義。本文在河北平原地下水壓采工程區，以Landsat數據為遙感數據源，利用改進的歸一化差異水體指數(MNDWI)方法，提取了2000、2007和2014年的坑塘水體信息。

(1)通過對比NDWI、MNDWI和標準假彩色合成3種方法在提取水體時的精度，認為改進的歸一化差異水體指數(MNDWI)在提取坑塘這種面積較小、與周圍差異不明顯的細小水體時較為準確。

(2)河北省地下水壓采工程區內坑塘主要分布在西部、北部；坑塘的面積和個數呈現明顯減少的趨勢；坑塘個數由2000年的13 937個減少到2014年的8 423個；坑塘面積由2000年的12 644.1 hm2減少到2014年的4 835.2 hm2；較之坑塘個數，坑塘面積萎縮更為嚴重。

(3)通過實測估算出的坑塘蓄水量對冬小麥的澆灌和減少地下水的開采都有意義。隨著坑塘灌溉保障作用進一步增強，它將能成為華北平原農村糧食增產和人居環境改善的重要手段。

[1] Chen J Y, Tang C Y, Shen Y J, et al. Use of water balance calculation and tritium to examine the dropdown of groundwater table in the piedmont of the North China Plain (NCP)[J]. Environmental Geology, 2003, 44:564-571.

[2] 胡春勝, 程一松. 中國生態系統定位觀測與研究數據集:農田生態系統卷, 河北欒城站 (1998-2008)[M]. 北京: 中國農業出版社, 2011:137-149.

[3] Piao S L, Ciais P, Huang Y, et al. The impacts of climate change on water resources and agriculture in China[J]. Nature, 2010,467(7311):43-51.

[4] 劉乙玄, 關 欣. 坑塘在土地生態系統中的功能研究[J]. 安徽農業科學, 2011,39(29):18 110-18 111.

[5] 王艷穎，王沛芳，徐海波. 坑塘濕地對生源要素氮的凈化效果[J].河海大學學報(自然科學版)，2006,34(6):630.

[6] 宋求明, 熊立華, 肖 義, 等. 基于 MODIS 遙感影像的洞庭湖面積與水位關系研究[J]. 節水灌溉, 2011,(6):20-26.

[7] Verpoorter C, Kutser T, Tranvik L. Automated mapping of water bodies using Landsat multispectral data[J]. Limnol. Oceanogr. Methods, 2012,(10):1 037-1 050.

[8] Verpoorter C, Kutser T, Seekell D A, et al. A global inventory of lakes based on high‐resolution satellite imagery[J]. Geophysical Research Letters, 2014,41(18):6 396-6 402.

[9] 羅 玲, 毛德華, 王宗明, 等. 松嫩平原西部湖泊與水庫/坑塘動態變化及驅動因素分析[J]. 農業工程學報, 2015,(22):285-291.

[10] 吳喜芳, 沈彥俊, 張 叢, 等. 基于植被遙感信息的作物蒸散量估算模型——以華北平原冬小麥為例[J]. 中國生態農業學報, 2014,22(8)：920-927.

[11] 高明秀，李 冉，鞏騰飛，等.“渤海糧倉”建設推進策略分析[J].農業現代化研究，2015，36(2)：245-251．

[12] McFeeters S K. The use of Normalized Difference Water Index(NDWI) in the delineation of open water features [J].International Journal of Remote Sensing, 1966,17(7):1 425-1 432．

[13] 徐涵湫.利用改進的歸一化差異水體指數(MNDWI)提取水體信息的研究[J].遙感學報，2005,9(5):589-595．

[14] Frazier PS, Page KJ. Water body detection and delineation with Landsat TM data[J].Photogrammetric Engineering and Remote Sensing,2000,66:1 467-2 467．

[15] Xiaofeng Gu, Jinxin Shi, Jianping Li. Application of wavelets analysis in image denoising[C]∥ IEEE International Conference on Apperceiving Computing and Intelligence Analysis,2008:49-52．

[16] 周在明, 張光輝, 王金哲, 等. 環渤海低平原微咸水區土壤鹽漬化與鹽分剖面特征[J]. 地理科學, 2011,31(8):929-934.

[17] 王秀萍, 張國新, 魯雪林, 等. 河北沿海區耕地土壤質量綜合評價[J]. 中國農學通報, 2013,29(30):136-142.

[18] Yang X, Lu X. Drastic change in China's lakes and reservoirs over the past decades[J]. Scientific reports, 2014,4:6 041.