于嵐嵐

（遼寧省江河流域管理局，遼寧 沈陽(yáng)110003）

0 引 言

蒸散發(fā)對(duì)于流域水循環(huán)的水量平衡方程來(lái)說(shuō)，是一個(gè)重要的輸出項(xiàng)，蒸散發(fā)大小的合理計(jì)算將會(huì)直接影響整個(gè)流域的水量計(jì)算。流域蒸散發(fā)資料主要是通過個(gè)別蒸散發(fā)觀測(cè)資料獲得，以點(diǎn)代面，通過點(diǎn)觀測(cè)蒸散發(fā)資料，得到整個(gè)流域的蒸散發(fā)資料，特別是在水文模型的輸入中，不能考慮流域的地形、土地利用對(duì)流域蒸散發(fā)的影響，造成流域蒸散發(fā)輸入的誤差。此文選用分布式雙源蒸散發(fā)模型，以淮河上游息縣水文站以上區(qū)域?yàn)檠芯苛饔?，基于流域?shù)字高程數(shù)據(jù)、土地利用數(shù)據(jù)以及流域內(nèi)外附近8個(gè)氣象站點(diǎn)的氣象要素?cái)?shù)據(jù)，計(jì)算了研究區(qū)域的蒸散發(fā)時(shí)空分布。

1 研究方法及資料準(zhǔn)備

1.1 研究區(qū)域概況

選擇淮河上游息縣水文站以上流域作為研究區(qū)域，面積為10 190 km2。流域地處我國(guó)濕潤(rùn)氣候與半干旱氣候過渡帶，淮南屬長(zhǎng)江流域氣候，一般溫暖多雨，淮北屬黃河流域型氣候，寒冷少雨，流域年平均降水量800~1 400 mm，多年平均徑流深約為371 mm，年平均水面蒸發(fā)量為800～1 000 mm，流域內(nèi)主要的土地利用方式為耕地和林地。

1.2 雙源蒸散發(fā)模型

雙源蒸散發(fā)模型，是將植被冠層和土壤分別當(dāng)做兩個(gè)相互獨(dú)立但又相互作用的2個(gè)通量源匯項(xiàng)，模型可以分別對(duì)土壤表層和植被冠層進(jìn)行能量的平衡計(jì)算,模型計(jì)算方程如下所示：

1）植被冠層截留蒸發(fā)E（i被植被葉面截留的那部分水量的蒸發(fā)）

2）植被蒸騰能力計(jì)算（土壤含水量達(dá)到田間持水量時(shí)，植被葉面氣孔所能蒸發(fā)的水分）

3）土壤蒸發(fā)能力（土壤含水量到達(dá)田間持水量時(shí)，裸土或者植物冠層以下土壤所蒸發(fā)的水分）

其中：Rnc——植物冠層所能得到的太陽(yáng)凈輻射量，W/m2；Rns——土壤表面所能獲得到的太陽(yáng)凈輻射量，W/m2；G——土壤的熱通量，Wm2；γ——空氣中的濕度常數(shù)，KPa/℃；△——飽和水汽壓梯度，KPa/℃；ρ——平均空氣的密度，kg/m3；CP——空氣的比熱，1.013×10-3kJkg-1/℃，λ——蒸發(fā)的潛熱，MJ/kg；Wfr——潮濕冠層比例；rcp——土壤含水量達(dá)到田間持水量時(shí)的冠層總氣孔阻抗，s/m；rsp——土壤含水量達(dá)田間持水量時(shí)的土壤表面阻抗，s/m；ras——土壤的表面與植物冠層源匯高度之間的空間動(dòng)力學(xué)阻抗，s/m；rac——植物冠層總邊界層的阻抗，s/m；D0——植物冠層源匯高度處的水汽壓強(qiáng)差，kPa。

1.3 模型數(shù)據(jù)處理

模型需要的數(shù)據(jù)包括以下幾項(xiàng)：

1）數(shù)字高程數(shù)據(jù)：數(shù)字高程數(shù)據(jù)選用美國(guó)國(guó)家地球物理中心提供的全球分辨率為1 km×1 km基礎(chǔ)地形數(shù)據(jù)，運(yùn)用ARCGIS軟件結(jié)合流域邊界，可得到息縣水文站以上流域的數(shù)字高程圖。

2）土地利用數(shù)據(jù)：基于中國(guó)科學(xué)院提供的2000年全國(guó)空間分辨率為1 km×1 km土地利用圖，結(jié)合流域邊界，應(yīng)用ArcGIS軟件進(jìn)行切割，得到流域2000年的土地利用圖。

3）氣象資料：收集流域附近武漢站2000—2010年凈。收集了流域內(nèi)外附近8個(gè)氣象站2000—2010年日平均氣溫、日平均水汽壓以及日平均風(fēng)速。

4）蒸發(fā)資料：收集了流域內(nèi)息縣蒸發(fā)站2000—2010年實(shí)測(cè)蒸發(fā)皿蒸發(fā)資料。

2 模型計(jì)算結(jié)果檢驗(yàn)與分析

基于已構(gòu)建的考慮植被葉面積指數(shù)動(dòng)態(tài)變化的雙源蒸散發(fā)模型，該研究計(jì)算了淮河息縣以上流域內(nèi)13 588個(gè)柵格單元2000—2010年的逐日蒸散發(fā)能力，選擇2個(gè)年份繪制了息縣站日計(jì)算蒸散發(fā)能力與實(shí)測(cè)蒸發(fā)皿蒸發(fā)日過程對(duì)比圖，見圖1。從圖1中可以看出計(jì)算的蒸散發(fā)能力和蒸發(fā)皿蒸發(fā)量波動(dòng)趨勢(shì)基本一致。

圖1 息縣站2001年和2010年蒸發(fā)皿與雙源蒸散發(fā)日過程比較圖

2.1 日尺度相關(guān)性分析

為進(jìn)一步檢驗(yàn)分布式雙源蒸散發(fā)模型計(jì)算的蒸散發(fā)能力與蒸發(fā)皿實(shí)測(cè)蒸發(fā)資料在時(shí)間上的相似性和一致性，通過建立實(shí)測(cè)蒸發(fā)皿蒸發(fā)與計(jì)算的蒸散發(fā)能力一元線性回歸模型和相關(guān)性的分析來(lái)進(jìn)行檢驗(yàn)。

在息縣流域內(nèi)共有3個(gè)蒸發(fā)站，分別為南灣站、石山口站以及息縣站，基于這3個(gè)蒸發(fā)站所在柵格單元2000—2010年的逐日蒸散發(fā)能力計(jì)算結(jié)果，建立這3個(gè)站實(shí)測(cè)蒸發(fā)皿與計(jì)算蒸散發(fā)能力之間的一元回歸方程，并采用F檢驗(yàn)］對(duì)實(shí)測(cè)蒸發(fā)皿和計(jì)算蒸散發(fā)能力之間的線性擬合程度進(jìn)行了檢驗(yàn)，結(jié)果見表1。

表1 蒸發(fā)皿日資料與所在柵格單元日蒸散發(fā)能力回歸分析

從表1中可以看出，3個(gè)站的日蒸發(fā)皿蒸發(fā)資料與其所在網(wǎng)格計(jì)算的蒸散發(fā)能力回歸方程都通過顯著性水平a=0.01的F檢驗(yàn)，其建立的回歸方程的斜率都是正數(shù)，表明計(jì)算的日蒸散發(fā)能力與各站的蒸發(fā)皿蒸發(fā)呈現(xiàn)正相關(guān)性。

2.2 蒸散發(fā)空間分布計(jì)算

基于流域內(nèi)各柵格日蒸發(fā)能力計(jì)算成果，并統(tǒng)計(jì)了各柵格年蒸發(fā)能力計(jì)算成果，并繪制息縣以上流域2000—2010年的年均蒸散發(fā)能力空間分布圖，可以發(fā)現(xiàn)，蒸散發(fā)量空間分布大體和高程分布具有一定的相似性，海拔較高的區(qū)域，氣溫較低，對(duì)蒸散發(fā)需求較低，因此蒸散發(fā)相對(duì)較??；而在海拔較低的區(qū)域，氣溫相對(duì)較高，蒸散發(fā)較為強(qiáng)烈，蒸散發(fā)量相對(duì)較大。成蒸散發(fā)量在700～1 000 mm的區(qū)域主要集中在海拔較低的區(qū)域，而蒸散發(fā)量在400～1 700 mm的區(qū)域主要集中在海拔較高的區(qū)域。

3 結(jié) 論

通過在淮河上游息縣水文站以上區(qū)域構(gòu)建分布式雙源蒸散發(fā)模型，基于研究區(qū)域數(shù)字高程、土地利用數(shù)據(jù)、以及研究區(qū)內(nèi)外附近8個(gè)氣象站點(diǎn)2000—2010年氣象要素?cái)?shù)據(jù)，運(yùn)用該模型計(jì)算了淮河息縣以上流域2000—2010年日蒸散發(fā)能力，并結(jié)合流域內(nèi)實(shí)則蒸發(fā)資料對(duì)模型計(jì)算結(jié)果進(jìn)行了檢驗(yàn)，檢驗(yàn)結(jié)果表明計(jì)算的蒸散發(fā)能力和實(shí)測(cè)的蒸發(fā)皿蒸發(fā)在時(shí)間上具有較好的相似性和一致性，模型可用來(lái)計(jì)算研究區(qū)域的蒸散發(fā)能力。

［1］郭瑾.老哈河流域蒸散發(fā)計(jì)算及其氣候變化對(duì)蒸發(fā)皿蒸發(fā)影響的研究［C］.南京：河海大學(xué)，2008.

［2］楊雨亭,尚松浩.雙源蒸散發(fā)模型估算潛在蒸散發(fā)量的對(duì)比［J］.農(nóng)業(yè)工程學(xué)報(bào),2012,28（24）：85-91.

［3］Zhou，M.C.，Ishidaira，H.P.，Magome，J.，et al.Estimating potential evaporation using the Shuttleworth-Wallace model and NOAA-AVHRR NDVIto feed a distributrd hydrological modeling over the Mekong River basin［J］.Journal of Hydrology，2006，327：151-173.

［4］王璨，周秀平，王文圣.窟野河洪水序列變異點(diǎn)綜合診斷［J］.水電能源科學(xué)，2012，30（7）：50-53.

［5］馬曉超，粟曉玲，薄永占.渭河生態(tài)水文特征變化研究［J］.水資源與水工程學(xué)報(bào)，2011，22（1）：16-21.