基于再分析降水與融雪模塊的雅礱江上游徑流模擬研究

黃 琦，覃光華，2，向俊燕，王瑞敏

（1.四川大學(xué)水利水電學(xué)院，成都610065；2.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，成都610065）

0 引 言

中國西部部分高寒地區(qū)氣候條件惡劣，人口密度小，社會經(jīng)濟(jì)發(fā)展落后，流域內(nèi)受人類活動影響較小，但稀疏的水文站網(wǎng)對水文模型的建立極其不利。雅礱江流域新龍水文站控制區(qū)域?qū)俑吆貐^(qū)，水文資料較少［1］，且在降水、植被指數(shù)、海拔、自然環(huán)境等方面差異明顯［2，3］。而雅礱江中下游近年來開發(fā)程度較大，上游兩河口水電站將于2021年投產(chǎn)，其上游新龍水文站控制流域?yàn)槠渲匾獜搅鹘M成部分。因此，采取合適的水文模型及模型輸入數(shù)據(jù)模擬新龍流域徑流，對雅礱江中下游水資源開發(fā)利用具有重大意義。

目前，雅礱江流域上游地區(qū)在水文特性及徑流預(yù)報方面已取得一定成果。張新海等［1］利用水文氣象站點(diǎn)觀測資料，分析研究了甘孜以上區(qū)域降水、暴雨和徑流等水文要素的特性；王渤權(quán)等［4］分析了雅礱江流域徑流年內(nèi)分配特性，認(rèn)為雅礱江上游徑流年內(nèi)分配不均；部分學(xué)者采用站點(diǎn)徑流資料，基于統(tǒng)計(jì)模型預(yù)測雅礱江二灘流域中長期徑流，但數(shù)學(xué)模型無明確水文機(jī)制，模型不確定性較大［5，6］；俞烜等［7］以站點(diǎn)資料為輸入，開發(fā)了雅礱江流域水循環(huán)系統(tǒng)的分布式水文模型，對雅礱江流域月徑流進(jìn)行模擬預(yù)報；康爾泗等［8］考慮了隨著高程變化帶來降水量和蒸發(fā)量等水文要素的差異性分布，采用站點(diǎn)資料推求流域和各海拔帶降水和氣溫，并基于HBV 模型完成雅礱江流域月、年平均流量的模擬。由此可見，目前研究區(qū)徑流模擬大都以地面站點(diǎn)降水為輸入進(jìn)行模擬計(jì)算，且模型結(jié)構(gòu)改進(jìn)有限。隨著衛(wèi)星遙感數(shù)據(jù)在數(shù)據(jù)稀缺地區(qū)水文模型的日益廣泛應(yīng)用，選擇高精度的格網(wǎng)降水產(chǎn)品并應(yīng)用成為未來研究的必然趨勢；同時，采用適當(dāng)?shù)乃哪Ｐ湍M徑流對研究區(qū)徑流模擬及預(yù)報工作具有重大意義。

本文首先基于地面降水資料，評價不同時間尺度上再分析格網(wǎng)降水資料（CMFD-P）的精度；其次基于三水源新安江模型和含融雪模塊的三水源新安江模型，分別以地面降水觀測資料和CMFD-P 為輸入，模擬了新龍水文站控制流域的徑流過程并進(jìn)行評價和比較，對（CMFD-P）在研究區(qū)的水文效應(yīng)進(jìn)行評估，進(jìn)而為新龍上游流域未來徑流預(yù)報發(fā)展提供參考。

1 研究區(qū)概況

雅礱江起源于青海省玉樹縣境內(nèi)的巴顏喀拉山南麓，流向?yàn)樽晕鞅敝翓|南，于呷依寺附近進(jìn)入四川省，是金沙江左岸最大的一條支流［1］。本文研究區(qū)為新龍水文站控制的雅礱江流域上游區(qū)，位于青藏高原東南部，屬高原亞寒帶濕潤區(qū)，下墊面變化較大。流域降水觀測站點(diǎn)稀少，且降水分布的地區(qū)差異較大，年內(nèi)分配嚴(yán)重不均勻［9］。降雪為流域降水的重要組成部分，同時春季融雪徑流在多年徑流中占有較大比重［10］，區(qū)域內(nèi)徑流年內(nèi)分配不均［4］。如圖1所示，研究區(qū)地理范圍為30.77°N～34.21°N，96.82°E～101.44°E；海拔范圍為3 173～6 087 m；流域大致呈南北向條帶狀，平均海拔高程為4 430.49 m，干流長約534 km，控制面積36 685 km2。在研究時段內(nèi)（2004-2012年），研究區(qū)內(nèi)無大型工業(yè)或采礦業(yè)，人口密度小，無大型水庫或水電站，總體而言，研究區(qū)在研究時段受人類活動影響較小，植被仍然保存完好。本文研究區(qū)共四個氣象站點(diǎn)（由北至南分別為清水河站、石渠站、甘孜站、新龍站），一個水文站點(diǎn)（新龍站）。研究區(qū)范圍及氣象水文站點(diǎn)如圖1所示。

2 數(shù)據(jù)及方法

2.1 研究數(shù)據(jù)

中國區(qū)域高時空分辨率地面氣象要素驅(qū)動數(shù)據(jù)集（China Meteorological Forcing Dataset，CMFD）為一套近地面氣象與環(huán)境要素再分析數(shù)據(jù)集［11］，包含降水在內(nèi)的7 個要素，以Princeton再分析資料、GLDAS 資料、GEWEX-SRB 輻射資料，TRMM 降水資料等為背景場，融合了中國氣象局常規(guī)氣象觀測數(shù)據(jù)制作而成，時間分辨率為3 h，水平空間分辨率為0.1°。其降水部分（后稱CMFD-P）具有良好的時空分辨率，在中國區(qū)域具有較高精度［12］。其背景降水資料TRMM 在不同流域均具有較高精度［13］，且在高寒山區(qū)流域具有較強(qiáng)的適應(yīng)性［14］；因此CMFD-P 不僅在疏勒河的不同地區(qū)優(yōu)于ERA-Interim、GPCC、GPCP V2.3 和JRA-55等其他4種降水產(chǎn)品［15］；在青藏高原與TRMM 具有更高一致性且表現(xiàn)優(yōu)于CHIRPS［16］，且比IMERG，MSWEP 更適用于雅礱江流域［17］。

本文流域格網(wǎng)數(shù)據(jù)為CMFD-P 格網(wǎng)降水?dāng)?shù)據(jù)；采用中國氣象數(shù)據(jù)網(wǎng)發(fā)布的中國地面氣象站日值降水?dāng)?shù)據(jù)（下稱地面站點(diǎn)數(shù)據(jù)）作為評估CMFD-P 的基準(zhǔn)資料；水文模型輸入降水?dāng)?shù)據(jù)分別為CMFD-P 格網(wǎng)降水資料面平均值和研究區(qū)內(nèi)4個地面站點(diǎn)降水?dāng)?shù)據(jù)的平均值；潛在蒸發(fā)數(shù)據(jù)基于彭曼公式，采用CMFD再分析數(shù)據(jù)集結(jié)合中國地面氣候資料日值數(shù)據(jù)集計(jì)算而得；模型率定資料為2004-2010 新龍水文站日徑流數(shù)據(jù)，模型檢驗(yàn)資料為2011-2012 新龍水文站日徑流數(shù)據(jù)。研究數(shù)據(jù)詳情見表1。

表1 研究資料概況Tab.1 Research data

潛在蒸散發(fā)數(shù)據(jù)采用國際糧農(nóng)組織FAO 提供的彭曼公式［18］計(jì)算，其計(jì)算所需基礎(chǔ)數(shù)據(jù)中，溫度、濕度、風(fēng)速、輻射數(shù)據(jù)由CMFD 提供；日照數(shù)據(jù)由中國地面氣候資料日值數(shù)據(jù)集（V3.0）提供并結(jié)合克里金插值得到各格網(wǎng)日照數(shù)據(jù)；緯度數(shù)據(jù)由DEM格網(wǎng)資料獲取。計(jì)算公式如下：

式中：PET為可能蒸散量，mm/d；Δ 為飽和水氣壓曲線斜率，kPa/℃；Rn為地表凈輻射，MJ/（m·d）；G為土壤熱通量，MJ/（m2·d）；γ為干濕表常數(shù)，kPa/℃；Tmean為日平均溫度，℃；u2為2 m 高處風(fēng)速，m/s；es為飽和水氣壓，kPa；ea為實(shí)際水氣壓，kPa。

2.2 研究方法

2.2.1 水文模型

本文采用集總式三水源新安江模型模擬徑流和含融雪模塊的集總式三水源新安江模型，含融雪模塊的新安江原理如圖2所示。新安江模型由河海大學(xué)趙人俊教授等提出，并以近代山坡水文學(xué)和國內(nèi)外產(chǎn)匯流理論的研究成果為基礎(chǔ)，改進(jìn)成三水源新安江模型，在濕潤地區(qū)或半濕潤地區(qū)更為適用［19］。新安江模型包括蒸散發(fā)模塊、產(chǎn)流模塊、水源劃分模塊和匯流模塊；含融雪模塊的新安江模型基于三水源新安江模型，增加了融雪模塊（圖2中虛線框中部分），其融雪徑流采用度日因子法進(jìn)行計(jì)算。

圖2中參數(shù)含義如下：a為度日因子，KC為蒸散發(fā)折算系數(shù)，WUM上層土壤蓄水容量，WLM為下層土壤蓄水容量，C為深層蒸散發(fā)系數(shù)，WM為流域蓄水容量，B為流域蓄水容量曲線指數(shù)，IM為流域不透水面積比值，SM為流域自由水蓄水容量，EX為流域自由水蓄水容量曲線指數(shù)，KI與KG分別為自由水蓄水庫壤中流出流系數(shù)與自由水蓄水庫地下徑流出流系數(shù)，CS為地面徑流消退系數(shù)，CI為壤中流消退系數(shù)，CG 為地下徑流消退系數(shù)，L為滯后時間，KE為槽蓄系數(shù)，XE為流量比重因素。

其中，度日因子法基本原理如下［20］：

式中：M為日雪水當(dāng)量，cm；a為度日因子，cm/（℃·d），表示溫度每上升1 ℃所產(chǎn)生的融雪徑流深；T為度日因子數(shù)，℃-1/d，即日正積溫。

三水源新安江模型和含融雪模塊的三水源新安江模型的模型參數(shù)優(yōu)化方法均采用遺傳算法，算法的種群數(shù)量均為200，迭代次數(shù)均為100。

2.2.2 降水?dāng)?shù)據(jù)評價方法

本文基于新龍水文站控制流域4 個地面氣象站點(diǎn)2000-2017年降水資料，通過分類指標(biāo)和定量指標(biāo)評價CMFD-P 對實(shí)際降水的探測能力和模擬精度。

（1）探測能力評價。降水事件分類指標(biāo)體現(xiàn)了降水?dāng)?shù)據(jù)集對降水事件的監(jiān)測能力，包括探測概率探測率POD（Probability Of Detection）、擊中率FOH（Frequency Of Hit）和Heidke 技巧指數(shù)HSS（Heidke’s Skill Score）［21］，分別表示了降水?dāng)?shù)據(jù)的正確探測能力，擊中能力和綜合判別能力。各指標(biāo)計(jì)算公式如下：

式中：n11代表CMFD-P 和雨量站均探測到降水的頻次；n10表示CMFD-P 探測到降水而站點(diǎn)未監(jiān)測到降水的頻次；n01表示站點(diǎn)監(jiān)測到降水而CMFD-P 未監(jiān)測到降水的頻次；n00表示CMFD-P和雨量站均未監(jiān)測到降水的頻次。POD、FOH、HSS三者理想值均為1。

（2）精度評價。本文以相對偏差BIAS、絕對偏差A(yù)BIAS、相關(guān)系數(shù)R、納什效率系數(shù)（Nash-Sutcliffe Efficiency coefficient）NSE來評估降水?dāng)?shù)據(jù)集在該地區(qū)的系統(tǒng)偏差、數(shù)據(jù)相關(guān)性和數(shù)據(jù)序列模擬能力，并對不同時間尺度的精度進(jìn)行分析。各指標(biāo)計(jì)算公式如下：

式中：P代表CMFD-P，G代表雨量站點(diǎn)降水，n 為計(jì)算精度指標(biāo)所采用的數(shù)據(jù)系列長度。其中R變化范圍為-1～1；BIAS變化范圍為-∞～∞；ABIAS變化范圍為0～∞；NSE變化范圍為-∞～1。BI?AS和ABIAS越趨近0，R、NSE越趨近1，CMFD-P 精度越高。本文采用點(diǎn)對像素法（站點(diǎn)所在格網(wǎng)中心點(diǎn)數(shù)據(jù)與站點(diǎn)數(shù)據(jù)）來進(jìn)行資料的關(guān)聯(lián)與比較。

2.2.3 徑流模擬評價方法

本文將2004-2010年徑流資料作為率定期，2011-2012年徑流資料作為模型效果驗(yàn)證的基準(zhǔn)資料。本文采用納什效率系數(shù)（NSE）和合格率QR（Qualified Ratio）評價模型。納什效率系數(shù)為水文常用評價模型確定性手段，NSE約接近1，模型模擬結(jié)果越好。當(dāng)0.4＜NSE≤0.6 和NSE＞0.6 時，分別表示模型模擬結(jié)果校好和優(yōu)秀。模型評價時，以NSE為主，QR為輔，判別徑流的模擬結(jié)果。其計(jì)算公式分別為：

式中：NSE為納什效率系數(shù)（保留2 位小數(shù)），QR表示合格率，均無量綱；QSIM為預(yù)報值，m3/s；QOBS為實(shí)測值，m3/s；-QOBS為實(shí)測值的均值，m3/s；n為資料長度；k表示合格預(yù)報次數(shù)；m表示總預(yù)報次數(shù)。本文模擬日徑流過程，當(dāng)預(yù)報許可誤差（即徑流絕量對偏差A(yù)BIAS）小于20%時認(rèn)為其合格。

3 結(jié)果分析

3.1 降水產(chǎn)品評價

本文降水評價基礎(chǔ)資料為新龍水文站控制流域4個地面氣象站點(diǎn)（見圖1清水河、石渠、甘孜、新龍）2000-2017年降水資料，評價產(chǎn)品為CMFD-P，圖3為日、月、年、春季（3-5月）、夏季（6-8月）、秋季（9-11月）、冬季（12-次年2月）7 種時間尺度4 個站點(diǎn)所在格網(wǎng)CMFD-P 與站點(diǎn)降水散點(diǎn)圖。由于四季降水?dāng)?shù)據(jù)并非為連續(xù)自然時間序列，故未統(tǒng)計(jì)四季NSE。

由圖3可知，在研究區(qū)，CMFD-P在月、年、春季（3-5月）、夏季（6-8月）、秋季（9-11月）、冬季（12-次年2月）與站點(diǎn)降水均具有較高一致性，表2總體評價了各統(tǒng)計(jì)變量在4 個站點(diǎn)的平均值。在日降水探測中，POD 達(dá)到0.91，具有較高探測精度，F(xiàn)OH 為0.70，HSS 為0.63，表明CMFD-P 對降水探測能力明顯高于隨機(jī)。同時，CMFD-P 在月降水探測中也有較高的準(zhǔn)確率（POD為1.00，F(xiàn)OH為0.98，HSS為0.53）。

由圖3和表2知，除日降水外，其他時段CMFD-P 在研究區(qū)具有較低ABIAS，和較高R和NSE；其中冬季節(jié)ABIAS（0.14）較其他季節(jié)更高，說明在冬季時CMFD-P 描述降水偏差略微偏大；CMFD-P 在春夏季略微低估降水（BIAS 分別為-0.010 和-0.007），在秋冬季略微高估降水（BIAS 分別為0.007 和0.021）。造成CMFD-P 冬季精度略差的原因可能是CMFD-P 對固態(tài)降水的估計(jì)能力不足。日降水的相關(guān)系數(shù)R和NSE均較其他時間尺度低，隨著時間尺度的累積，二者均增加。相關(guān)系數(shù)R在四季尺度上分布差異不大。綜上，CMFD-P 在研究區(qū)與站點(diǎn)降水具有較高一致性。

表2 CMFD-P精度評價表Tab.2 Accuracy evaluation of the precipitation in CMFD-P

圖4為4 個站點(diǎn)不同時期降水定量評價精度指標(biāo)。其中清水河、石渠、甘孜、新龍站點(diǎn)地理位置依次為由西北至東南（圖1）。由圖4知，CMFD-P 在各時間尺度均低估清水河降水（BIAS均小于0）；而石渠除夏季外，整體輕微高估降水（BIAS均大于等于0）；CMFD-P 高估新龍站冬季降水（BIAS為0.11），低估甘孜冬季降水（BIAS為-0.03）。CMFD-P 在除日降水外的各時間尺度絕對系統(tǒng)偏差均較小（ABIAS均接近0）；但甘孜和新龍?jiān)诙鞠到y(tǒng)偏差絕對值A(chǔ)BIAS較其他兩個站點(diǎn)更大。從降水相關(guān)性來說，CMFD-P 與各站點(diǎn)降水相關(guān)性均較高，各站點(diǎn)間R值相差甚小，日降水R均小于0.75，隨著時間的累積，其他時間尺度降水相關(guān)性升高，R大于0.75。

可見，在系統(tǒng)偏差方面，4 個站點(diǎn)系統(tǒng)偏差整體較小，清水河在四季均表現(xiàn)為負(fù)偏，新龍站在冬季表現(xiàn)為較高的正偏；而其他站點(diǎn)四季偏差具有抵消作用；在系統(tǒng)絕對偏差方面，各站點(diǎn)在某一時間尺度下表現(xiàn)相似；絕對偏差在不同時間尺度具有一定差異，日尺度最高，冬季次之，其他時間尺度在低值程度上相近。在站點(diǎn)降水與CMFD-P 的相關(guān)性方面，各站點(diǎn)模擬能力在各時段有差異，但在站點(diǎn)之間差異甚小；相關(guān)性隨著時間尺度的增加而升高。CMFD-P 降水模擬能力可能與地形、氣候等因素有關(guān)，由于本文研究區(qū)地面站點(diǎn)稀疏，因此降水精度空間差異不明顯，研究降水與自然屬性關(guān)系時的代表性不足。因此，未來在獲取更多數(shù)據(jù)的基礎(chǔ)上應(yīng)當(dāng)更加深入探討影響研究區(qū)CMFD-P降水精度的因素。

3.2 徑流模擬評價

如圖5為分別以CMFD-P 和地面降水為輸入率定的新安江模型和含融雪模塊的新安江模型模擬日徑流結(jié)果圖，圖6為按時間累積后的月徑流模擬結(jié)果圖。

由圖5、圖6知，相比三水源新安江模型，含融雪模塊的新安江模型能更好地模擬日、月徑流過程，尤其對枯期（春、冬季）徑流的模擬優(yōu)于三水源新安江模型。但含融雪模塊的新安江模型對部分年份豐期（夏、秋季）具有明顯低估現(xiàn)象。

表3進(jìn)一步分析了圖5和圖6的詳細(xì)評價結(jié)果，其中月徑流評價由日徑流累加結(jié)果評價而得。由表3知，以地面站點(diǎn)降水作為三水源新安江模型的輸入率定的模型總體優(yōu)于以CMFD-P為輸入時率定的模型，前者在各時段BIAS，ABIAS均小于后者，且前者各時段R、NSE、QR均略高于后者。除此之外，兩種數(shù)據(jù)率定的三水源新安江模型在檢驗(yàn)期效果均優(yōu)于率定期。這可能是由于研究區(qū)融雪徑流占徑流較大比重，而三水源新安江模型不含融雪模塊，因此無法合理模擬研究區(qū)徑流。同時，本文研究區(qū)積雪在2004-2012年具有逐年下降的趨勢［22］，在檢驗(yàn)期時段（2011-2012年）研究區(qū)融雪徑流占的比重較率定期（2004-2010年）少，故而該模型在檢驗(yàn)期表現(xiàn)更好。

表3還展示了含融雪模塊的三水源新安江模型在各個時期的評價結(jié)果。結(jié)果表明，在日徑流模擬中，以CMFD-P 為輸入的模型優(yōu)于以站點(diǎn)數(shù)據(jù)為輸入的模型，二者均低估實(shí)際徑流（BIAS小于0），但CMFD-P 模擬結(jié)果負(fù)偏程度小。二者模擬結(jié)果的ABIAS和QR接近，但CMFD-P 模擬結(jié)果的R和NSE均高于站點(diǎn)降水模擬結(jié)果。在月徑流模擬中，CMFD-P 模擬結(jié)果相比站點(diǎn)降水模擬結(jié)果負(fù)偏程度小（BIAS小于0 且更接近0）；與實(shí)際結(jié)果總體一致性NSE更高；絕對偏差A(yù)BIAS和R相近；但合格率QR更高，體現(xiàn)了CMFD-P 模擬結(jié)果在豐期ABIAS大于20%的更高可能性。含融雪模塊的三水源新安江模型模擬結(jié)果BI?AS均為負(fù)值，表明改進(jìn)的模型對徑流具有低估現(xiàn)象，這與圖5-6中豐水期低估現(xiàn)象明顯的結(jié)論相符合。在日、月徑流模擬中，二者率定期結(jié)果均優(yōu)于檢驗(yàn)期結(jié)果，且均優(yōu)于無融雪模塊的三水源新安江模型率定期模擬結(jié)果，表明含融雪模塊的三水源新安江模型結(jié)構(gòu)合理穩(wěn)定，能夠很好地適用于研究區(qū)；CMFD-P模擬結(jié)果優(yōu)于站點(diǎn)降水模擬結(jié)果，說明由于站點(diǎn)稀疏，站點(diǎn)降水模擬徑流結(jié)果較差，而CMFD-P 模擬結(jié)果與實(shí)際徑流具有更高的一致性。

為進(jìn)一步評價不同輸入數(shù)據(jù)和模型結(jié)構(gòu)組合情況下對徑流的預(yù)報能力，分析比較表3中以站點(diǎn)降水、CMFD-P 為輸入的三水源新安江模型及含融雪模塊的三水源新安江模型在檢驗(yàn)期時的表現(xiàn)情況。從模型結(jié)構(gòu)來看，含融雪模塊的三水源新安江模型在日徑流模擬時，檢驗(yàn)期的BIAS、R均未優(yōu)于三水源新安江模型，但ABIAS、NSE和QR優(yōu)于三水源新安江模型；在月徑流模擬時，檢驗(yàn)期的BIAS未優(yōu)于三水源新安江模型，但ABIAS、NSE和QR優(yōu)于三水源新安江模型，R與之相近；表明加入融雪模塊后，模擬結(jié)果雖然整體負(fù)偏，日徑流的線性相關(guān)性降低，但模型整體模擬能力NSE得到提升。從輸入數(shù)據(jù)來看，CMFD-P作為三水源新安江模型輸入時，模擬結(jié)果未優(yōu)于站點(diǎn)降水模擬結(jié)果；但當(dāng)CMFD-P 作為含融雪模塊的三水源新安江模型的輸入時，日徑流模擬結(jié)果的負(fù)偏情況較站點(diǎn)輸入模擬結(jié)果小，ABIAS，R，NSE與之相當(dāng)，但QR減小；月徑流模擬結(jié)果的負(fù)偏情況較站點(diǎn)輸入模擬結(jié)果小，NSE增大，ABIAS，R與之相當(dāng)，QR減小。由圖CMFD-P 得到QR更小的原因可能是由于其對枯期徑流模擬結(jié)果值更大，而枯期實(shí)測徑流基數(shù)較小，因此，CMFD-P較站點(diǎn)降水模擬結(jié)果的QR下降幅度大。

表3 不同降水?dāng)?shù)據(jù)徑流模擬能力評價表Tab.3 Simulation ability evaluation of runoff calculated by different precipitation datasets

表4分析了三水源新安江模型和含融雪模塊的新安江模型在全時段（2004-2012年）春、夏、秋、冬四季模擬結(jié)果與實(shí)際徑流的水量偏差。由表4知，站點(diǎn)降水和CMFD-P 在采用三水源新安江模型模擬徑流時，對四季徑流模擬的偏差趨勢類似，即秋季高估（BIAS均大于0），春、冬季低估（BIAS均小于0），夏季略微低估。兩種數(shù)據(jù)模擬結(jié)果在夏秋季表現(xiàn)較秋冬季更好，站點(diǎn)降水模擬結(jié)果整體偏差略低于CMFD-P 模擬結(jié)果，且二者在夏秋季差異較春、冬季小。

表4 四季徑流模擬精度評價表（2004年-2012年）Tab.4 Runoff simulation accuracy at four seasons（2004-2012）

站點(diǎn)降水和CMFD-P 在采用含融雪模塊的三水源新安江模型模擬徑流時，四季偏差具有差異性。站點(diǎn)降水模擬結(jié)果在四季均未高估徑流，且在夏季低估最明顯（BIAS為-0.22），春、秋季次之（BIAS分別為-0.15 和-0.11），冬季偏差最小（BIAS約為0）。CMFD-P 模擬結(jié)果在冬季高估徑流（BIAS為0.43），在夏季低估最（BIAS為-0.16），春、秋季分別略微高估和低估徑流（BIAS分別為0.08 和-0.16）。CMFD-P 對冬季徑流模擬偏差大于站點(diǎn)降水模擬結(jié)果（BIAS和ABIAS分別為0.43 和0.45），這與3.1中CMFD-P在冬季與地面站點(diǎn)降水偏差較大的結(jié)果一致，表明CMFD-P 對固態(tài)降水的探測具有局限性；但在其他季節(jié)，尤其是春季模擬結(jié)果偏差幅度較站點(diǎn)降水模擬偏差小，表明CMFD-P能更好地描述春季徑流。

綜上，作為格網(wǎng)降水，CMFD-P 可以更全面、精確地描述區(qū)域降水，對徑流整體模擬能力較好。在采用三水源新安江模型時，模擬結(jié)果略低于站點(diǎn)降水模擬結(jié)果；在采用更適合研究區(qū)的含融雪模塊的三水源新安江模型時，優(yōu)于站點(diǎn)降水模擬結(jié)果，在冬季高估徑流，削減了春、夏秋季的徑流估計(jì)偏差。由于該地區(qū)融雪在春季為徑流的重要組成部分［10，23］，因此，融雪模塊的加入優(yōu)化了模擬結(jié)果；且CMFD-P 對流域空間尺度上的降水描述更詳細(xì)，優(yōu)于站點(diǎn)降水模擬結(jié)果。

4 結(jié)論與展望

本文首先基于地面降水資料，評價了不同時間尺度上CMFD-P 的精度；其次基于三水源新安江模型和含融雪模塊的三水源新安江模型，分別以地面降水觀測資料和CMFD-P 為輸入，通過2004-2010年雅礱江上游地區(qū)新龍水文站控制流域的日徑流來率定模型，并通過預(yù)測2011-2012年該水文站日徑流來檢驗(yàn)?zāi)Ｐ停煌瑫r將日徑流累加到月徑流對模型評價。結(jié)果表明：

（1）在降水?dāng)?shù)據(jù)精度方面：CMFD-P 在研究區(qū)4 個站點(diǎn)所在格網(wǎng)與地面站點(diǎn)降水一致性較高，在月、年、四季尺度NSE均大于0.9，因此，認(rèn)為該數(shù)據(jù)在研究區(qū)具有高精度的特點(diǎn)。

（2）在模型結(jié)構(gòu)優(yōu)化方面：含融雪模塊的三水源新安江模型在研究區(qū)比三水源新安江模型具有更高的適用性，有效優(yōu)化了模擬結(jié)果。

（3）在模型輸入方面：采用含融雪模塊的三水源新安江模型模擬徑流時，CMFD-P降水可以更全面，精確地描述研究區(qū)降水，得到較高的徑流模擬精度；優(yōu)于站點(diǎn)降水模擬結(jié)果，尤其提升了春季徑流精度。說明在描述區(qū)域降水時，CMFD-P 面平均值代表性優(yōu)于地面站點(diǎn)降水均值。

雅礱江上游徑流預(yù)報對兩河口及以下的水資源調(diào)度及梯級電站開發(fā)具有重要意義。作為集總式概念性水文模型，集總式新安江模型無法準(zhǔn)確反應(yīng)水文過程在流域的空間變異性，而以CMFD-P 計(jì)算流域面平均降水，可以減少以地面站點(diǎn)降水描述流域降水的誤差。另外，本研究區(qū)為青藏高原東南部一個狹長型流域，屬于高原亞寒帶濕潤地區(qū)，融雪占徑流比例大［10］，含融雪模塊的三水源新安江模型能夠更好地模擬融雪徑流，因此模擬新龍水文站出流過程能力較不含融雪模塊的三水源新安江模型強(qiáng)。由于研究區(qū)由南至北下墊面差異較大，在未來該流域的研究工作中，首先可以從模型結(jié)構(gòu)上完善該流域水文模型徑流形成的物理機(jī)制；其次，采用高精度的格網(wǎng)降水?dāng)?shù)據(jù)CMFD-P，并充分利用其空間優(yōu)勢，通過建立分布式水文模型，實(shí)現(xiàn)更完整精確的研究區(qū)徑流模擬及預(yù)報。 □