王 東

（遼寧省水文局，遼寧沈陽110003）

?

柵格蒸散發(fā)模型在區(qū)域蒸散發(fā)時空分布模擬中的應(yīng)用研究

王 東

（遼寧省水文局，遼寧沈陽110003）

摘要：文章應(yīng)用地理信息技術(shù)結(jié)合區(qū)域數(shù)字高程數(shù)據(jù)，將區(qū)域離散成1km×1km的單元網(wǎng)格，應(yīng)用雙源蒸散發(fā)模型構(gòu)建基于柵格尺度的蒸散發(fā)模型，進行區(qū)域蒸散發(fā)時空模擬研究，并和區(qū)域內(nèi)的實測蒸發(fā)皿蒸發(fā)進行對比分析。研究結(jié)果表明：柵格蒸散發(fā)模型和實測蒸發(fā)皿蒸發(fā)年相關(guān)系數(shù)為0.6575，在春季相關(guān)性最高，在夏季相關(guān)性最低。柵格蒸散發(fā)模型可應(yīng)用于區(qū)域的蒸散發(fā)計算。蒸散發(fā)空間分布主要受到區(qū)域高程的影響，和高程分布較為一致。研究成果對于區(qū)域蒸散發(fā)時空模擬以及無資料地區(qū)的蒸散發(fā)計算提供參考價值。

關(guān)鍵詞：柵格蒸散發(fā)模型；雙源蒸散發(fā)模型；蒸散發(fā)空間模擬；蒸散發(fā)計算

0 引言

流域蒸散發(fā)空間分布的研究一直是國內(nèi)外學者研究的熱點問題，早期，由于技術(shù)條件的限制，常常以點蒸發(fā)皿蒸發(fā)來代替區(qū)域面蒸發(fā)，但區(qū)域內(nèi)蒸發(fā)站點較少，很難以點代面，且這種以點代面的計算方法，沒有考慮流域自然地理條件（土地利用、數(shù)字高程）空間上的差異性對流域潛在蒸散發(fā)的影響。其次，在無資料地區(qū)暴雨洪水計算時，由于缺少實測蒸發(fā)皿蒸發(fā)資料，需要運用蒸散發(fā)模型來計算區(qū)域的蒸散發(fā)。為此，當前構(gòu)建基于柵格的蒸散發(fā)模型來模擬計算區(qū)域面蒸散發(fā)，成為許多學者研究和計算區(qū)域蒸散發(fā)的有效途徑，并取得一定的研究成果［1-5］。這些研究大都基于P-M公式進行的柵格蒸散發(fā)計算研究，而雙源蒸散發(fā)模型由于可獨立考慮植被和土壤的蒸散發(fā)，被許多學者應(yīng)用于區(qū)域的蒸散發(fā)模擬計算中，也取得一定的研究成果［6-8］，但構(gòu)建基于柵格尺度的雙源蒸散發(fā)模型的研究相對較少，為此，本文基于柵格尺度的雙源蒸散發(fā)模型，進行區(qū)域蒸散發(fā)時空模擬研究，并和區(qū)域內(nèi)的實測蒸發(fā)皿蒸發(fā)進行對比分析，研究成果對于區(qū)域蒸散發(fā)時空模擬以及無資料地區(qū)的蒸散發(fā)計算提供參考價值。

1 柵格蒸散發(fā)模型原理

基于地理信息技術(shù)，結(jié)合研究區(qū)域數(shù)字高程數(shù)據(jù)，將區(qū)域空間離散成1km×1km的單元網(wǎng)格，結(jié)合雙源蒸散發(fā)模型，構(gòu)建基于柵格的蒸散發(fā)模型。

1.1 雙源蒸散發(fā)模型單獨計算植被和土壤的潛在蒸散發(fā)計算

（1）植物葉面氣孔蒸散發(fā)計算。

式中：Rnc為植物葉面冠層能夠吸收到的太陽的輻射量，Wm-2；ρ為空氣的平均密度，kgm-3；Cp為空氣的比熱常數(shù)值，1.013×10-3KJkg-1℃-1；D0為植物葉面冠層與蒸發(fā)源之間的水汽壓的壓強差，kPa1；rac為植物葉面冠層邊界總的氣孔阻抗，sm-1；λ為區(qū)域蒸發(fā)的潛熱值，MJkg-1；Δ為水汽壓的飽和梯度值，kPa℃-1；γ為空氣濕度的一個常數(shù)，kPa℃-1；Wft為植物冠層潮濕的比例。

（2）植被根系含水經(jīng)過葉面氣孔蒸發(fā)的計算。

式中：Rns為土壤表面上能吸收到的太陽的輻射量，Wm-2；rcp為土壤含水量達到田間持水量后土壤表層所產(chǎn)生的阻抗，sm-1。

（3）區(qū)域土壤表層的水分蒸發(fā)量計算。

式中：G為通過土壤表層的熱能通量，Wm-2；ras為土壤表層與植物葉面冠層之間高度所產(chǎn)生的空間動力阻抗值，sm-1；

1.2 流域水面的蒸發(fā)計算

在雙源蒸散發(fā)計算流域蒸散發(fā)能力的基礎(chǔ)上考慮流域水面的蒸發(fā)計算，蒸發(fā)計算公式如下：

式中：Rn為水面所獲得的凈輻射量，Wm-2；G為水體的熱通量，Wm-2；U2為2m處的風速，m/s；es為飽和水氣壓強，kPa；ea為實際水汽壓強，kPa。

3 模型數(shù)據(jù)

模型需要輸入的數(shù)據(jù)有土地利用數(shù)據(jù)、數(shù)字高程數(shù)據(jù)以及蒸發(fā)、氣象數(shù)據(jù)，其中土地利用數(shù)據(jù)和土地利用數(shù)據(jù)下載來源見參考文獻［9-10］，土地利用數(shù)據(jù)、及數(shù)字高程數(shù)據(jù)下載后結(jié)合流域邊界獲取研究區(qū)域的土地利用數(shù)據(jù)、數(shù)字高程數(shù)據(jù)，數(shù)據(jù)見圖1。蒸發(fā)皿數(shù)據(jù)采用區(qū)域內(nèi)3處蒸發(fā)站（1#蒸發(fā)站、2#蒸發(fā)站、3#蒸發(fā)站），其中1#蒸發(fā)站位于上游、2#蒸發(fā)站位于下游，1#蒸發(fā)站位于中游區(qū)域。氣象數(shù)據(jù)采用附近區(qū)域氣象站2001～2010年氣象數(shù)據(jù)。

4 模型應(yīng)用

4.1 年尺度模擬

為對比柵格蒸散發(fā)模擬的適用性，結(jié)合研究區(qū)內(nèi)3個實測蒸發(fā)站2000～2010年實測蒸發(fā)皿蒸發(fā)數(shù)據(jù)，對比柵格蒸散發(fā)模擬的適用性，并結(jié)合假設(shè)檢驗的方法對柵格蒸散發(fā)模型計算的蒸散發(fā)和柵格蒸發(fā)皿蒸發(fā)進行了相關(guān)性分析。分析結(jié)果見圖1、圖2及表1。

表1 柵格蒸散發(fā)計算與實測蒸發(fā)皿線性回歸分析結(jié)果

圖1 柵格蒸散發(fā)模型和實測蒸發(fā)皿過程對比圖

圖2 柵格蒸散發(fā)模型和實測蒸發(fā)皿過程相關(guān)圖

表1為區(qū)域柵格蒸散發(fā)模型計算的蒸散發(fā)和3個蒸散發(fā)站實測蒸發(fā)皿蒸發(fā)相關(guān)分析結(jié)果，從表中可以看出，3個蒸發(fā)站的蒸發(fā)皿數(shù)據(jù)和計算的蒸發(fā)均通過了顯著性a=0.01的假設(shè)檢驗分析。從表中可以看出相關(guān)系數(shù)為正數(shù)，則表示計算的蒸散發(fā)和實測蒸散發(fā)在年尺度是正相關(guān)性。從表中的回歸方程斜率和截距可看出，3個站點的回歸方程具有一定的相似性，但也有一定的差異性，這主要是由于各站點所在的不同的植被類型所造成的，1#蒸發(fā)站的植被覆蓋為林地，2#蒸發(fā)站的植被覆蓋為草地，而3#蒸發(fā)站的植被類型為耕地，從表中可以看出植被為林地的1#蒸發(fā)站站斜率和截距都大于其他2站，其次是植被類型為草地的2#蒸發(fā)站，植被類型為耕地的2#蒸發(fā)站的斜率和截距都小于其他2站，由此也可以看出，就蒸發(fā)能力而言，林地的年蒸發(fā)能力最大，其次是草地，耕地年蒸發(fā)能力最小。以3#蒸發(fā)站為列，本文繪制了3#蒸發(fā)站2000～2010年年計算蒸發(fā)能力與年實測蒸發(fā)皿蒸發(fā)散點的相關(guān)圖及其回歸分析。從圖2中可以看出，3#蒸發(fā)站2000～2010年年計算蒸發(fā)能力與年實測蒸發(fā)皿蒸發(fā)散點的相關(guān)達到了0.6575，相關(guān)系數(shù)較高。

4.2 季節(jié)尺度模擬

基于柵格蒸散發(fā)模型，計算了月尺度的蒸散發(fā)，并和3個蒸發(fā)皿所在柵格計算的蒸散發(fā)進行月尺度回歸分析，分析結(jié)果見表2。

表2 柵格蒸散發(fā)模型和實測蒸發(fā)皿季相關(guān)性分析

表2為區(qū)域3個蒸發(fā)站實測蒸發(fā)皿數(shù)據(jù)和柵格蒸散發(fā)模型在季節(jié)尺度的相關(guān)性分析結(jié)果，從表中可以看出，不同季節(jié)實測蒸發(fā)皿數(shù)據(jù)和柵格蒸散發(fā)模型在季節(jié)尺度的相關(guān)性分析有所差異，其中，在春季，由于植被開始生長，區(qū)域的潛在蒸散發(fā)受植被蒸散發(fā)影響較小，主要受到氣象條件的影響。因此和蒸發(fā)皿蒸發(fā)的相關(guān)性都好于其他3個季節(jié)。進入夏季后，植物處于生長最為旺盛的時期，植物蒸散發(fā)占據(jù)區(qū)域蒸散發(fā)的比重較大。使得在夏季，柵格蒸散發(fā)模型計算的潛在蒸散發(fā)和實測蒸發(fā)皿的蒸散相關(guān)系數(shù)較低。進入秋季以后，植被開始凋落，區(qū)域潛在蒸散發(fā)量受到植被蒸散發(fā)影響減弱，和蒸發(fā)皿蒸發(fā)的相關(guān)性加大。進入冬季后，植物開始凋零，柵格蒸散發(fā)模型計算的潛在蒸散發(fā)和實測蒸發(fā)皿的相關(guān)系數(shù)大于秋季的相關(guān)系數(shù)。

4.3 蒸散發(fā)空間模擬

利用構(gòu)建的柵格蒸散發(fā)模型，計算了區(qū)域2000～2010年蒸散發(fā)的時空分布。從計算結(jié)果可以看出，區(qū)域多年平均蒸散發(fā)整體呈現(xiàn)北大南小的趨勢，這主要是因為區(qū)域蒸散發(fā)很大程度上受到區(qū)域高程的影響，區(qū)域高程高值主要分布在北邊，南北區(qū)域?qū)儆诤０屋^低的區(qū)域。對于高程較高的區(qū)域，蒸發(fā)量受高程影響一般較小，而在南邊海拔較低的區(qū)域，受到高程影響，蒸發(fā)量較低。

5 結(jié)論

本文建立柵格蒸散發(fā)模型，實現(xiàn)區(qū)域蒸散發(fā)時空模擬，并結(jié)合實測蒸發(fā)皿蒸發(fā)進行相關(guān)性分析，研究取得以下結(jié)論：

（1）柵格蒸散發(fā)模型可考慮下墊面條件差異對蒸散發(fā)的影響，和實測蒸散發(fā)相關(guān)系數(shù)在年尺度的相關(guān)系數(shù)好于月尺度上的相關(guān)系數(shù)；

（2）蒸散發(fā)空間分布主要受高程影響，海拔較低的區(qū)域，蒸散發(fā)量一般加大，反之，海拔較高的區(qū)域，蒸散發(fā)量較小。

參考文獻

［1］俞健.于什蓋水庫蒸發(fā)與滲漏量分析［J］.水利規(guī)劃與設(shè)計，2014（07）：45-47.

［2］王積強，陸旭，劉巽民.中國水面蒸發(fā)器的發(fā)展簡史與相關(guān)技術(shù)問題探討［J］.水利技術(shù)監(jiān)督，2011（03）：9-11+39.

［3］王積強.推薦E_（601）型水面蒸發(fā)器的檢測方法［J］.水利技術(shù)監(jiān)督，2010（01）：14-15.

［4］王海波，馬明國.基于遙感和Penman-Monteith模型的內(nèi)陸河流域不同生態(tài)系統(tǒng)蒸散發(fā)估算［J］.生態(tài)學報，2014（19）：5617-5626.

［5］李紅霞，張永強，張新華，等.遙感Penman-Monteith模型對區(qū)域蒸散發(fā)的估算［J］.武漢大學學報（工學版），2011（04）：457-461.

［6］徐靜，任立良，劉曉帆，等.基于雙源蒸散與混合產(chǎn)流的Pa1mer旱度模式構(gòu)建及應(yīng)用［J］.水利學報，2012（05）：545-553.

［7］楊雨亭，尚松浩.雙源蒸散發(fā)模型估算潛在蒸散發(fā)量的對比［J］.農(nóng)業(yè)工程學報，2012，24：85-91.

［8］劉曉帆，任立良，袁飛，等.雙源蒸散發(fā)能力計算模型在半干旱區(qū)的適用性［J］.河海大學學報（自然科學版），2009（02）：138-142.

［9］土地利用數(shù)據(jù)：http：//westdc.westGIS.ac.cn/data/1cad1a63-ca8d-431a-b2b2-45d9916d860d/.

［10］DEM數(shù)據(jù)：http：//srtm.csi.cgiar.org/SELECTION/inputCoord.asp.

中圖分類號：P333

文獻標識碼：B

文章編號：1008-1305（2016）02-0053-03

DOI：10.3969/j.issn.1008-1305.2016.02.023

收稿日期：2015-12-28

作者簡介：王 東（1983年—），男，工程師。