吳建鵬 于海洋 馬慧慧

（河南理工大學(xué)礦山空間信息技術(shù)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室 焦作 454000）

1 引言

SWAT模型是一種具有長(zhǎng)時(shí)序連續(xù)模擬能力的分布式流域水文模型，在晉江流域［1］、灤河流域［2］、渭河流域［3］等國(guó)內(nèi)水文模擬中取得了廣泛的應(yīng)用驗(yàn)證［4］。其中土地利用/覆被數(shù)據(jù)是構(gòu)建SWAT模型的基礎(chǔ)數(shù)據(jù)，已有的研究中主要采用Landsat等中低分辨率衛(wèi)星傳感器數(shù)據(jù)解譯獲取，數(shù)據(jù)精度相對(duì)較低。國(guó)內(nèi)新近完成的地理國(guó)情普查提供了大量精確的土地利用分類數(shù)據(jù)［5］，如何將其擴(kuò)展應(yīng)用到相關(guān)地學(xué)領(lǐng)域是近年來(lái)研究的熱點(diǎn)［6～7］。本文以鶴壁湯河流域?yàn)檠芯繀^(qū)，結(jié)合流域的地理國(guó)情普查數(shù)據(jù)，制定了湯河流域的土地利用分類體系，構(gòu)建SWAT水文模型，對(duì)比分析地理國(guó)情普查數(shù)據(jù)與Landsat影像解譯的土地利用/覆被數(shù)據(jù)在SWAT流域水文模型模擬中的適用性，為區(qū)域土地利用規(guī)劃、生態(tài)恢復(fù)及水資源合理配置提供理論依據(jù)［8］。

2 SWAT水文模擬模型

SWAT模型中由土壤侵蝕子模型、污染負(fù)荷子模塊和水文過(guò)程子模型三大子模型組成，本文研究涉及到水文過(guò)程子模型。水文過(guò)程模擬研究分為兩個(gè)部分，分別為水文循環(huán)陸面部分和河道匯流部分，SWAT模型計(jì)算遵守水量平衡原理，其原理方程見(jiàn)下式：

式中：SWm是第n天的土壤最終含水量（mm）；SWo是第n天的土壤前期含水量（mm）；m是時(shí)間（d）；Rd是第n天的降水量（mm）；Qsu是第n天的地表徑流量（mm）；E是第n天的蒸散發(fā)量（mm）；Wse是第n天存在于土壤剖面地層的滲漏量和側(cè)流量（mm）；Qgw是第n天的地下水含量（mm）。

3 地理國(guó)情普查數(shù)據(jù)分類系統(tǒng)的整合方法

本文主要使用湯河流域2015年的地理國(guó)情地表覆蓋數(shù)據(jù)，數(shù)據(jù)獲取主要是以湯河流域的2015年經(jīng)過(guò)預(yù)處理的高分辨率影像作為地理國(guó)情普查的基礎(chǔ)地圖數(shù)據(jù)，同時(shí)收集各項(xiàng)專題數(shù)據(jù)如二調(diào)數(shù)據(jù)等為地理國(guó)情數(shù)據(jù)作為參考；本文獲取同時(shí)期的Landsat遙感影像，對(duì)其進(jìn)行預(yù)處理并結(jié)合獲取的流域范圍內(nèi)相關(guān)資料，在圖像分析過(guò)程中利用面向?qū)ο螅?］的方法，首先采用合適的分割尺度對(duì)影像進(jìn)行分割，然后采用SVM［10～11］分類方法并結(jié)合目視解譯的方法得到湯河流域土地利用類型分布圖。根據(jù)研究區(qū)土地利用的實(shí)際情況，本文將研究區(qū)土地利用類型劃分為耕地、林地、草地、園地、工礦用地、高密度居住地、低密度居住地、水域和裸地等9類，由于我國(guó)的土地利用分類標(biāo)準(zhǔn)不符合SWAT模型下的美國(guó)制土地利用類型代碼的要求，因此需要對(duì)地理國(guó)情數(shù)據(jù)進(jìn)行重分類，轉(zhuǎn)換成SWAT模型識(shí)別的土地利用類型分類代碼（見(jiàn)表1）。

4 數(shù)據(jù)處理與結(jié)果分析

4.1 研究區(qū)概況

研究區(qū)選取鶴壁市湯河上游流域作為研究對(duì)象。流域位于河南省北部，太行山東麓和華北平原的過(guò)渡地帶。湯河屬海河流域衛(wèi)河水系，湯河上游流域主要包括鶴壁市鶴山區(qū)和山城區(qū)（圖1），面積約165 km2。研究區(qū)內(nèi)湯河水庫(kù)建于1958年，其上游水文、氣象觀測(cè)數(shù)據(jù)完整，同時(shí)為消除水庫(kù)人為調(diào)節(jié)的影響，徑流觀測(cè)選取在水庫(kù)匯水區(qū)上游區(qū)域。

表1 湯河流域土地利用重分類表

圖1 研究區(qū)位置

4.2 數(shù)據(jù)預(yù)處理

根據(jù)圖2中湯河流域的分類結(jié)果，本文以湯河流域2015年30m分辨率Landsat影像解譯結(jié)果和2015年地理國(guó)情數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)，各土地利用/覆被類型所占面積如表2所示。

由圖2和表2可知，兩期土地利用數(shù)據(jù)各土地利用類型所占面積比例相似，空間分布存在一定差異。2015年湯河流域土地利用類型以耕地、林地、草地、高密度居住地、低密度居住地和礦業(yè)用地為主，這六種土地利用類型的面積之和占流域總面積的94%以上，而裸地、果園和水域的面積相對(duì)較小。

圖3 研究區(qū)土地利用分類結(jié)果

表2 研究區(qū)不同來(lái)源土地利用類型面積統(tǒng)計(jì)結(jié)果

4.3 SWAT模型構(gòu)建

構(gòu)建SWAT模型需要獲取研究區(qū)范圍內(nèi)包括高程數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)、氣象數(shù)據(jù)等一系列數(shù)據(jù)，然后將其輸入模型，最后將在研究區(qū)范圍內(nèi)獲取的實(shí)測(cè)水文數(shù)據(jù)作為對(duì)模擬值的適用性進(jìn)行校正和驗(yàn)證的依據(jù)。

研究中采用機(jī)載LiDAR生成的DEM數(shù)據(jù)為水文模擬分析的基礎(chǔ)［12～14］。機(jī)載 LiDAR 數(shù)據(jù)是在2009年4月21日采用Leica ALS50獲取的。氣象數(shù)據(jù)是以中國(guó)氣象局大氣同化系（CLDAS）［15］技術(shù)為參考，時(shí)間周期選擇在2008年1月1日至2014年12月31日期間。降雨量數(shù)據(jù)由安陽(yáng)和鶴壁水文局提供，主要包括小河子、鶴壁、盤石頭和施家溝等雨量監(jiān)測(cè)站。土壤數(shù)據(jù)是以河南省1：20萬(wàn)分縣土壤類型圖作為參考，該數(shù)據(jù)是根據(jù)土壤發(fā)生學(xué)的觀點(diǎn)，結(jié)合河南省土壤實(shí)際情況進(jìn)行劃分。

首先基于DEM進(jìn)行子流域劃分，然后在SWAT模型中加載土地利用數(shù)據(jù)、氣象數(shù)據(jù)、坡度數(shù)據(jù)以及土壤數(shù)據(jù)，生成湯河流域水文響應(yīng)單元。在上述過(guò)程中同時(shí)加載了煤礦抽排水點(diǎn)信息，在此模擬研究區(qū)內(nèi)五個(gè)礦區(qū)煤礦抽排水過(guò)程主要是采用點(diǎn)源輸入（Point Source Input）的形式。

模型運(yùn)行模擬完成后，將利用研究區(qū)范圍內(nèi)獲取的實(shí)測(cè)水文數(shù)據(jù)對(duì)模型模擬值的適用性進(jìn)行校正與驗(yàn)證。論文利用小河子、鶴壁、盤石頭和施家溝水文站的2008-2014年間逐日、逐月徑流的數(shù)據(jù)，設(shè)置2008為預(yù)熱期，運(yùn)用2009-2011年實(shí)測(cè)徑流數(shù)據(jù)進(jìn)行模型率定，模型驗(yàn)證利用2012-2014年實(shí)測(cè)數(shù)據(jù)。本次研究中在其他輸入數(shù)據(jù)不變的情況下，分別基于2015年地理國(guó)情數(shù)據(jù)和2015年30m分辨率下Landsat解譯的土地利用數(shù)據(jù)輸入模型，并在率定期和驗(yàn)證期模擬得到相關(guān)系數(shù)（R2）和Nash-Sutcliffe效率系數(shù)（Ens），相關(guān)研究認(rèn)為R2＞0.6、Ens＞0.5時(shí)模擬取得顯著效果。

4.4 結(jié)果分析

由表3、圖3和圖4可知，SWAT模型模擬的徑流量和實(shí)測(cè)徑流量吻合較好，由此證明了模型模擬的結(jié)果符合湯河流域徑流量的變化規(guī)律。模擬結(jié)果表明，基于地理國(guó)情普查數(shù)據(jù)率定期和驗(yàn)證期相關(guān)系數(shù)模擬結(jié)果（0.74、0.82）優(yōu)于同時(shí)期Landsat影像解譯的土地利用/覆被數(shù)據(jù)模擬結(jié)果（0.71、0.80），證明基于地理國(guó)情普查數(shù)據(jù)構(gòu)建SWAT模型模擬的月徑流值與實(shí)測(cè)徑流值吻合較好；分析原因是地理國(guó)情數(shù)據(jù)更接近于真實(shí)地表覆被信息，在輸入模型所需其他數(shù)據(jù)不變的情況下，基于地理國(guó)情數(shù)據(jù)的模擬精度高于低分辨率下解譯的土地利用/覆被數(shù)據(jù)的模擬結(jié)果。

表3 湯河流域不同來(lái)源土地利用數(shù)據(jù)下SWAT模型月徑流模擬結(jié)果評(píng)價(jià)

圖4 湯河流域基于30m分辨率Landsat影像解譯數(shù)據(jù)模擬月徑流模擬值與實(shí)測(cè)值對(duì)比圖

5 結(jié)語(yǔ)

本文選取鶴壁湯河流域作為研究區(qū)，獲取該流域同一時(shí)期不同來(lái)源的土地利用/覆被數(shù)據(jù)，進(jìn)行了分布式水文模型模擬實(shí)驗(yàn)，分析了不同來(lái)源土地利用/覆被數(shù)據(jù)在水文分析和模擬中的適用性。統(tǒng)計(jì)結(jié)果顯示，流域的各土地利用類型在各自流域所占面積相似，在輸入模型所需其他數(shù)據(jù)不變的情況下，基于地理國(guó)情數(shù)據(jù)的模擬精度高于低分辨率下解譯的土地利用/覆被數(shù)據(jù)的模擬結(jié)果，因此建議在利用SWAT進(jìn)行徑流模擬時(shí)，地理國(guó)情數(shù)據(jù)較更低分辨率的Landsat影像解譯的土地利用數(shù)據(jù)有一定的優(yōu)勢(shì)，能夠在一定程度上提高徑流模擬的效果。