馬 原

(遼寧省錦州水文局,遼寧 錦州 121000)

SWAT(Soil and Water Assessment Tool)是目前被廣泛應(yīng)用于流域污染負(fù)荷計算的模型，由美國農(nóng)業(yè)部(USDA)將SWRRB和ROTO模型綜合利用形成，目前已在50多個國家的不同流域得到應(yīng)用和推廣。水文學(xué)者Beven將水文模擬的特點主要分為尺度界限、等效分析、唯一性標(biāo)準(zhǔn)、非線性和不確定性等5個特征。近年來，利用DEM數(shù)字高程在水文學(xué)多個領(lǐng)域被廣泛的研究和應(yīng)用，其計算結(jié)果被科學(xué)、準(zhǔn)確地應(yīng)用于分布式水文模擬研究分析方面[1]。研究流域的河長、寬度、坡度以及流域邊界條件等河網(wǎng)參數(shù)的準(zhǔn)確提取受DEM空間分辨率的影響顯著[2]。一般情況下，DEM空間分辨率和精確度越高，對河流的信息提取和描述就越真實可靠，越能準(zhǔn)確反應(yīng)出研究流域下墊面層的水文結(jié)構(gòu)分布狀況，其分布式水文模型的模擬結(jié)果就越好。然而，現(xiàn)階段如何獲得高精度的DEM數(shù)據(jù)仍需要進(jìn)一步的研究和探討，且對于大尺度空間的研究流域，高分辨率的DEM數(shù)據(jù)可能會導(dǎo)致計算機(jī)運(yùn)行緩慢，甚至無法滿足計算的要求，不僅會明顯降低計算機(jī)的運(yùn)算速度，而且可能導(dǎo)致運(yùn)算結(jié)果的不準(zhǔn)確[3- 4]。

基于DEM數(shù)字高程模型的水文參數(shù)是分布式水文模型的重要輸入數(shù)據(jù)，其精度和分辨率將直接影響到研究流域的水文特征信息和徑流模擬結(jié)果的準(zhǔn)確性和可靠性，是進(jìn)行分布式水文模型模擬的重要基礎(chǔ)和參考依據(jù)。據(jù)此，文章開展相關(guān)DEM精確性對水文模型模擬結(jié)果的相關(guān)研究，為提高分布式水文模型輸入水文信息的真實性和有效性提供一定的理論支持和科學(xué)依據(jù)[5]。

1 區(qū)域概況及數(shù)據(jù)來源

1.1 區(qū)域概況

選擇大凌河流域為研究對象，該流域位于我國遼寧省西部，水土流失嚴(yán)重，泥沙含量較大，約為57kg/m3。資料顯示9月和10是該流域的降雨旺季，年降水量為450～600mm，徑流量約為16.67億m3。大凌河主要包括北、西、南三處發(fā)源區(qū)。北源區(qū)自北向南，徑流自凌源縣萬元店至凌源縣城南辛杖子匯入西源；南源自建昌縣經(jīng)喀左縣與西、北兩源匯流。三源匯流后由西南向東北方向徑流，沿線徑流包括朝陽縣、北票市、義縣、錦縣，最終流入渤海[6- 8]。

1.2 數(shù)據(jù)來源

文章參考借鑒了6個氣象站在1982～2014年的水文氣象數(shù)據(jù)，并以1982～1985年的降雨監(jiān)測數(shù)據(jù)以基準(zhǔn)，主要包括大凌河流域的降雨量、日照時長、風(fēng)速、溫度、濕度、蒸散發(fā)量、太陽輻射等水文參數(shù)；大凌河流域的數(shù)字高程DEM數(shù)據(jù)在GIS上截出，并用于水系提取和流域離散化處理，精度為1∶33萬；1∶25萬水系分布數(shù)據(jù)由水利部水文局提供；土壤類型圖精度為1∶33萬，包括土壤類型分布和土壤化學(xué)屬性；土地利用圖的精度為1∶10萬，由環(huán)境與生態(tài)科學(xué)研究中心提供。上述介紹的所有數(shù)據(jù)均采用統(tǒng)一的投影坐標(biāo)系統(tǒng)，即UTM_Zone_48N橢球投影坐標(biāo)系統(tǒng)。

2 研究方法

2.1 SWAT模型簡介

分布式SWAT水文模型被美國學(xué)者最先提出并得到研究和應(yīng)用，該模型具有參數(shù)設(shè)置少、所需資料少、操作簡便等優(yōu)點，被國內(nèi)外學(xué)者廣泛應(yīng)用于大尺度、長時段的水文模擬研究。大量的水文學(xué)者認(rèn)識到SWAT水文模型可廣泛應(yīng)用于水文資源模擬中，且模型模擬的結(jié)果精確度較高，能較好地反應(yīng)出研究流域的實際情況，是目前在水文研究領(lǐng)域應(yīng)用較為成熟的模型之一[9]。

2.2 研究方法

SWAT水文模型有不同的設(shè)置參數(shù)，可對研究流域在不同的水文影響下進(jìn)行水文過程模擬。TOPAZ軟件包是SWAT模型的主要計算元件，TOPAZ是進(jìn)行河網(wǎng)水系坡度、河長以及河網(wǎng)邊界等水文要素特征信息提取以及進(jìn)行子流域劃分的重要工具，SWAT模型是基于TOPAZ對水文要素特征信息的提取進(jìn)行徑流模擬的影響研究。故文章充分利用Arc Hydro Tools工具的計算特性和優(yōu)點，以TOPAZ模塊為基礎(chǔ)開展水文特征信息和受DEM空間分辨率大小的影響研究。地表徑流、河網(wǎng)匯流、蒸散發(fā)模型和數(shù)據(jù)差值算法均可在SWAT水文模型中應(yīng)用和計算，其已被廣泛應(yīng)用于水文模擬分析中，且模擬結(jié)果表現(xiàn)出較高的可靠性[10]。

2.3 數(shù)據(jù)處理

研究中利用NEAREST法分別對空間分辨率為30m的ASTER DEM數(shù)據(jù)進(jìn)行反復(fù)的取樣計算。為了便于計算和數(shù)據(jù)處理，文章采用25×2n指數(shù)法進(jìn)行反復(fù)取樣，當(dāng)n為0，1，2，……，7時，DEM空間分辨率分別為25m，50m，100m，200m，400m，800m，1600m，3200m。

2.4 評估標(biāo)準(zhǔn)

文章在分析比較的基礎(chǔ)上，為提高計算模型模擬效果的可靠性，選取效率系數(shù)(ENS)、相對誤差(Dv)以及徑流過程線得出確定性系數(shù)(R2)進(jìn)行模型模擬結(jié)果的評價。流域水量偏差估計一般采用徑流過程線和相對誤差(Dv)參數(shù)；模型的整體模擬效果以及準(zhǔn)確性分析通常采用效率系數(shù)(ENS)和確定性系數(shù)(R2)，計算方程如下所示。

(1)

式中，Dv—相對誤差；O—觀測的平均值；P—預(yù)測平均值。

(2)

式中，ENS—流域污染模擬效率系數(shù)；Oi—i時刻的觀測值；Pi—i時刻的模擬值；其他同上。

(3)

式中，R2—確定性系數(shù)；其他同上。

文章建立的水文模型模擬結(jié)果評價標(biāo)準(zhǔn)見表1。

表1 水文模擬結(jié)果評價標(biāo)準(zhǔn)

3 構(gòu)建SWAT模型并率定

3.1 建立模型

利用6個氣象站的降雨數(shù)據(jù)建立天氣發(fā)生器，并將15個雨量站在1982～1985年的日降雨量監(jiān)測數(shù)據(jù)輸入模型，收集DEM分辨率為30m的土地利用類型、土壤類型以及植被覆蓋率為研究流域的水文下墊面信息。河網(wǎng)集水面積閥值為120km2，以多種HRU劃分方式對土壤類型、面積等進(jìn)行設(shè)置，并設(shè)定最小面積閥值比為10%；產(chǎn)匯流計算利用SCS模型，并以馬斯京根法進(jìn)行匯流計算。

3.2 參數(shù)率定

文中選取敏感性較高的參數(shù)進(jìn)行率定，并選取典型的大城子水文站、朝陽水文站和凌海水文站進(jìn)行模型在月尺度徑流的參數(shù)率定。根據(jù)上述三個具有一定代表性的水文站實測數(shù)據(jù)，SWAT模型在不同時期的模擬結(jié)果見表2。

表2 月尺度徑流的模擬率定計算結(jié)果

由表2可知，月尺度徑流模擬結(jié)果均表現(xiàn)出良好的精度和準(zhǔn)確性，相對誤差均小于15%，三個水文站的效率系數(shù)在驗證期均大于0.80，在率定期在0.78～0.88之間，同樣表現(xiàn)出良好的準(zhǔn)確性，確定性系數(shù)R2值在0.88～0.95之間。驗證期較率定期總體計算值偏高，但變化規(guī)律基本一致，總體反映了研究區(qū)的實際情況，滿足流域面源污染模擬的精度規(guī)范要求。SWAT模型經(jīng)過率定后，可用于對其他方面的研究模擬計算，可利用SWAT模型探討水文參數(shù)及徑流量受空間分辨率的影響研究。

4 影響研究

4.1 水文特征信息提取影響

文章充分利用了Arc Hydro Tools工具的計算特性和優(yōu)點，然后以TOPAZ模塊為基礎(chǔ)建立的SWAT水文模型開展了水文特征信息提取受DEM空間分辨率大小影響的研究，8種DEM空間分辨率對水文特征信息的提取結(jié)果見表3。

表3 不同DEM空間分辨率下的大凌河流域水文特征信息匯總表

由表3可知，在不同的DEM分辨率下河流平均高度變化幅度最小，而最低高程隨著分辨率的減少逐漸增大，最高高程隨分辨率的減少而逐漸減少，研究表明河流的地勢結(jié)構(gòu)隨分辨率的降低而趨于平緩；不同分辨率條件下，模型對流域的河網(wǎng)總長計算結(jié)果存在較大差異，河網(wǎng)總長隨空間分辨率的降低而減少，DEM分辨率大于400m時河網(wǎng)總長變化幅度明顯增大；隨分辨率的降低大凌河流域面積和子流域個數(shù)呈減少發(fā)展趨勢，即分辨率越低對流域的信息判斷越不準(zhǔn)確，流域的邊界線分析趨于模糊。

4.2 徑流模擬影響研究

文章利用效率系數(shù)(ENS)、相對誤差(Dv)以及徑流過程線得出確定性系數(shù)(R2)進(jìn)行模型模擬結(jié)果的評價，計算結(jié)果見表4。

表4 各水文監(jiān)測站在不同分辨率下的模擬結(jié)果

結(jié)果表明：隨著空間分辨率的降低，文中所選的三個水文站其模擬徑流的效率系數(shù)和確定性系數(shù)均表現(xiàn)出降低的趨勢，當(dāng)分辨率在25～200m區(qū)間時，R2和ENS降低幅度較低，當(dāng)分辨率在400～3200m時，二者的降低幅度明顯提高；大凌河洪峰流量模擬結(jié)果與實測結(jié)果誤差隨DEM空間分辨率的降低而逐漸增大，模擬徑流擬合曲線吻合度逐漸朝著擴(kuò)散方向發(fā)展。采用相對誤差表征流域徑流模擬總量和實測總量之間的誤差，Dv值的大小代表模擬與實測之間的貼近程度，值越大則模擬結(jié)果偏差越大，值越小則模擬結(jié)果越好。研究表明：三個水文站的徑流量變化隨空間分辨率的降低無明顯的規(guī)律性，即DEM空間分辨率對流域的徑流量影響無明顯的相關(guān)性，而徑流總量的變化趨勢與Dv值的正負(fù)保持一致性。

5 結(jié)論

文章以大凌河流域為研究對象，在分析了SWAT分布式水文模型的基礎(chǔ)上，通過選取8中不同的DEM空間分辨率對研究流域的水文特征信息和徑流模擬影響進(jìn)行研究，結(jié)果表明：

(1)文章建立的SWAT模型，其月尺度徑流模擬結(jié)果均表現(xiàn)出良好的精度和準(zhǔn)確性，相對誤差均小于15%；三個水文站的ENS在驗證期均大于0.80，在率定期在0.78～0.88之間；R2值在0.88～0.95之間，模型適用于大凌河流域在月尺度的水文模擬研究。

(2)在不同的DEM分辨率下河流平均高度變化幅度最小，而最低、最高高程隨著分辨率的減少分別逐漸增大、減少。

(3)隨著空間分辨率的降低，文中所選的三個水文站其模擬徑流的ENS和R2均表現(xiàn)出降低的趨勢。DEM空間分辨率對流域的徑流量影響無明顯相關(guān)性，而徑流總量的變化趨勢與Dv值的正負(fù)保持一致性。

[1] 周祖昊, 仇亞琴, 賈仰文, 等. 變化環(huán)境下渭河流域水資源演變規(guī)律分析[J]. 水文, 2009(01): 21- 25.

[2] 姜廣田, 鄧程林. 遼寧省城市飲用水水源地安全評價研究[J]. 水利技術(shù)監(jiān)督, 2008(02): 24- 26.

[3] 孫龍, 臧文斌, 黃詩峰. DEM空間分辨率對流域水文特征信息提取及徑流模擬影響研究[J]. 水文, 2014(06): 21- 25.

[4] 全占東. 水足跡理論視角下遼河流域水資源評價[J]. 水土保持應(yīng)用技術(shù), 2017(01): 22- 24.

[5] 才慶欣. 南票區(qū)水資源狀況及開發(fā)利用分析[J]. 水利規(guī)劃與設(shè)計, 2014(09): 27- 29.

[6] 黃亮. 水資源可持續(xù)發(fā)展存在的問題及對策[J]. 水土保持應(yīng)用技術(shù), 2013(04): 32- 34.

[7] 李學(xué)森. 凌河流域水資源現(xiàn)狀及保護(hù)措施[J]. 水土保持應(yīng)用技術(shù), 2015(03): 36- 37.

[8] 張海軍. 水資源保護(hù)監(jiān)測存在問題及建設(shè)初探[J]. 水土保持應(yīng)用技術(shù), 2016(03): 40- 43.

[9] 曲錦艷. 遼寧省水資源存在的問題與保護(hù)對策[J]. 水土保持科技情報, 2002(06): 16- 17.

[10] 程海英. 錦州市城區(qū)水資源現(xiàn)狀及優(yōu)化配置方案[J]. 水土保持應(yīng)用技術(shù), 2015(06): 22- 23, 32.