新安江模型研究的回顧和展望

陸旻皎

（1.長岡技術(shù)科學(xué)大學(xué)，日本新潟縣長岡市 940-2188；2.重慶交通大學(xué)，重慶 400060）

1 研究背景

新安江模型是由河海大學(xué)（原華東水利學(xué)院）趙人俊教授所領(lǐng)導(dǎo)的團隊研發(fā)的降雨徑流模型，于1963年提出雛形，1973年正式發(fā)布［1-2］。該模型在濕潤、半濕潤地區(qū)得到廣泛應(yīng)用。其應(yīng)用范圍不僅僅局限于降雨徑流模擬及洪水預(yù)報，也包括水資源規(guī)劃管理、防災(zāi)、農(nóng)業(yè)等與水有關(guān)的領(lǐng)域。為國民經(jīng)濟發(fā)展發(fā)揮了重大貢獻，被列入中華人民共和國重大科技成果選集［3］。與很多世界同類知名模型，如薩克拉門托模型［4］、TOPMODEL［5］、水箱模型［6］和HBV 模型［7］等相比，在精度上也不遜色［8-9］。

新安江模型在中國及世界各國的廣泛應(yīng)用，使其結(jié)構(gòu)設(shè)計在不同氣候和地理條件下得到考驗。圖1和圖2為2018年8月初在百度學(xué)術(shù)和谷歌學(xué)術(shù)上搜索相關(guān)中英文關(guān)鍵詞得到的期刊論文數(shù)?？梢娦掳步Ｐ拖嚓P(guān)研究是相當(dāng)活躍的。學(xué)術(shù)研究如此，全國水文機構(gòu)也廣泛使用新安江模型。從1999年舉行的全國水文預(yù)報技術(shù)競賽［9］參賽隊伍可見一斑。此次競賽共有25 隊10 個模型參賽，其中14 個隊選擇了新安江模型。新安江模型的應(yīng)用流域數(shù)應(yīng)該是其他模型很難比肩的。這為其參數(shù)的區(qū)域規(guī)律研究創(chuàng)造了良好的條件。

新安江模型在國內(nèi)水文單位取得了壓倒性的利用，但在學(xué)術(shù)層面的應(yīng)用似乎有停滯或減少的趨勢。與國內(nèi)期刊論文數(shù)近十年基本保持不變或微幅增加相比（圖1），英文論文數(shù)有比較明顯的下降趨勢（圖2）。是新安江模型已經(jīng)趨于成熟，已經(jīng)變成諸多模型中的一個，無需進一步研發(fā)呢。還是隨著社會需求的多樣化，不能滿足各類用戶的需求，不能吸引新用戶。抑或是傳統(tǒng)的集總式已經(jīng)到了被分布式模型取代的時代。

圖1 百度學(xué)術(shù)新安江模型相關(guān)期刊論文數(shù)

圖2 谷歌學(xué)術(shù)新安江模型相關(guān)論文數(shù)

就分布式模型而言，在實際應(yīng)用中的優(yōu)勢并不明顯。日本所有主要流域都引進了分布式模型進行洪水預(yù)報，但精度等方面還有待改善。發(fā)揮分布式模型理論上的優(yōu)勢需要高精度分布式降水資料的穩(wěn)定獲取以及分布式水文參數(shù)的提取等方面的支撐。筆者認(rèn)為在這些方面還有很多工作要做。因此，在今后較長的一段時間內(nèi)分布式模型不可能完全取代傳統(tǒng)的集總式模型。模型用戶將會根據(jù)自己的需求做出選擇。作為實質(zhì)標(biāo)準(zhǔn)水文模型，新安江模型將繼續(xù)發(fā)揮其作用。

本文試圖回顧有關(guān)新安江模型研發(fā)的相關(guān)研究，冀望從中得到啟發(fā)，展望模型發(fā)展的方向。如有不周之處，還請指正。

2 新安江模型研究的回顧

期待新安江模型研發(fā)更注重可用性和易用性，本文主要對模型核心概念、參數(shù)空間特性、參數(shù)率定及參數(shù)可移植性進行了回顧。

2.1 蓄滿產(chǎn)流和現(xiàn)代產(chǎn)流機制研究愛爾蘭工程師Mulvaney 在1851年發(fā)表的經(jīng)驗公式［10］是第一個從雨量計算（洪峰）流量的公式。推理公式［11］雖然只是把式中的經(jīng)驗常數(shù)替換成徑流系數(shù)，卻跨出了把產(chǎn)流和匯流區(qū)分處理的一大步。而1930年代初由Sherman 提出的單位線理論［12］和Horton 提出的滲透公式［13］分別提供了匯流和產(chǎn)流計算手段，自然形成了預(yù)報流量過程線的強有力組合，對現(xiàn)代水文科學(xué)，水文模擬產(chǎn)生了深遠(yuǎn)的影響。全流域超滲產(chǎn)流的概念被廣泛接受。雨強大于下滲能力，即超滲成為產(chǎn)流的唯一機制，產(chǎn)生的地表流也被冠名為超滲地表流或Horton 地表流。Betson 于1964年證明即使是超滲產(chǎn)流也僅在流域的一部分發(fā)生，提出部分產(chǎn)流區(qū)（partial area）的概念［14］，才開始改變這種狀況。而山坡水文學(xué)研究得出的變動產(chǎn)流區(qū)（variable source area）的概念［15-17］對產(chǎn)流機制的認(rèn)知產(chǎn)生了更大的影響。在此基礎(chǔ)上，Dunne 和Black 提出了飽和產(chǎn)流機制［18-19］，由此產(chǎn)生的地表流也被命名為飽和地表流。在濕潤地區(qū)，土壤發(fā)育較好，下滲能力大，超滲產(chǎn)流難以發(fā)生。另一方面，由于雨量豐沛，地下水面一般比較淺，加上地形作用地下水向河道附近匯集，在其附近形成飽和或接近飽和的區(qū)域，一旦降雨，這個區(qū)域就會成為或迅速成為產(chǎn)流區(qū)，而且會在降雨過程中擴大或縮小。TOP?MODEL［5］正是由Beven and Kirkby 提出的基于飽和產(chǎn)流機制的水文模型。

中國的水文規(guī)律及產(chǎn)流機制的研究比國外大約要遲20年。在大規(guī)模水利建設(shè)的推動下從1950年代開始迅速發(fā)展，至1960年代中期就獲得了一批重要的研究成果［20］。由于客觀條件的限制，中國學(xué)者主要對降雨徑流關(guān)系及流量過程線進行了非常深入的研究［2，21］，并提出蓄滿產(chǎn)流機制，也就是土壤水分滿足土壤蓄水容量以后才開始產(chǎn)流。更進一步把蓄水容量的空間分布統(tǒng)計歸納后投影到一條分布曲線——即蓄水容量曲線上，得以非常簡潔地表達(dá)產(chǎn)流量和產(chǎn)流區(qū)的時變性和非線性特征。Beven［22］將新安江模型恰如其分地分類為概率分布式模型（probability distributed model）［23］。蓄水容量曲線可以說是新安江模型的核心，其統(tǒng)計投影手法也被很多水文模型，如VIC 模型［24］、ARNO 模型［25］所借鑒，為考慮流域水文特性的異質(zhì)性開創(chuàng)了一個新的途徑。

需要指出的是，蓄滿產(chǎn)流雖然與飽和產(chǎn)流非常相似，但還是有一些區(qū)別。超滲產(chǎn)流和飽和產(chǎn)流都是表述地表徑流的產(chǎn)流機制，而蓄滿產(chǎn)流則是描述總徑流量的產(chǎn)流機制。所謂蓄滿是指土層全層達(dá)到田間持水量，蓄滿以后進入土層的水分為自由水，最后在重力作用下以壤中流或地下水的形式流出。在新安江模型里，產(chǎn)流和分水源是分別通過引進張力水蓄水容量曲線和自由水蓄水容量曲線來處理的。

另一個潛在的問題是，蓄滿產(chǎn)流沒有考慮雨水進入土壤后的垂向運動，即下滲后的再分布［26］。濕潤地區(qū)一般植被良好，下滲能力很大，以日本為例，森林流域的穩(wěn)定下滲率平均在200 mm/h以上［27］，很少出現(xiàn)超滲產(chǎn)流，而且下滲的水分也在很短的時間內(nèi)完成再分布。加之當(dāng)時模型運用以日模型為主，這樣的處理有其合理性。但是隨著數(shù)據(jù)采集技術(shù)的進步，水文資料的時間分辨率會提高。同時社會也會對洪水預(yù)報提出更高的要求，時模型的利用會逐漸增加。有必要對再分布的影響進行評估。

2.2 蓄水容量曲線如前節(jié)所述，蓄水容量曲線是新安江模型的核心。迄今為止已經(jīng)有一些研究嘗試用地形、土壤和植被等信息推求蓄水容量［28-31］。這些研究的共同之處是，假設(shè)流域某一點或某一柵格的張力水蓄水容量（WM）和自由水蓄水容量（SM）為包氣帶內(nèi)的所有張力水和腐殖質(zhì)土層內(nèi)的所有自由水，即有：

式中：θs、θf、θwp分別為飽和含水量、田間持水量和凋萎含水量；Lv、Lh分別為包氣帶厚度和腐殖質(zhì)土層厚度。然后通過地形指數(shù)或土壤地形指數(shù)［5］等計算包氣帶和腐殖質(zhì)土層厚度得到蓄水容量。也有部分研究嘗試從地形指數(shù)累積曲線直接推導(dǎo)蓄水容量曲線［32-33］?？追舱芎退螘悦吞岢鲇玫乇砥露韧茖?dǎo)蓄水容量曲線的方法［34］。值得特別提出的是，向小華等提出了一種考慮毛細(xì)管邊緣（capillary fringe）的蓄水容量［35］，改善了式（1）不考慮水分毛細(xì)管上升的缺陷，可望提高對細(xì)顆粒土壤的適用性。

上述研究不但可以幫助確定有關(guān)模型參數(shù)，更重要的是可以使模型的核心概念有更明確的物理意義。但筆者認(rèn)為包氣帶厚度和蓄水容量不但取決于地形、土壤和植被，也應(yīng)該受氣候條件的影響。在此整理一下基本物理概念，考慮如何引進氣候因素。

筆者定義張力水蓄水容量為最大包氣帶張力水缺水量，為：

式中：Lv，max為包氣帶最大厚度；θ為體積含水量；z為垂向坐標(biāo)，以上方為正。包氣帶厚度是隨著時間變化的，其最大值代表該地點的干燥程度，也受氣候影響。由于最大厚度一般出現(xiàn)于很長的無降雨期的末期，含水量垂向分布比較穩(wěn)定，可以近似地假設(shè)毛管張力ψ和重力處于平衡狀態(tài)，即：

如果用Brooks-Corey 模型：

表達(dá)土壤水分特征曲線，可得含水量如下：

式中：ψb，λ為Brooks-Corey 模型參數(shù)；θr為滯留含水量。田間持水量出現(xiàn)在：

對式（3）積分可得：

通過式（6），只要有土壤參數(shù)和地下水資料就可以估算張力水蓄水容量。

2.3 模型參數(shù)和率定新安江模型參數(shù)按照模型結(jié)構(gòu)分為4 層：蒸散發(fā)計算參數(shù)；產(chǎn)流計算參數(shù)；分水源計算參數(shù)；匯流計算參數(shù)［36］。表1列出了三水源新安江模型相關(guān)參數(shù)及其物理意義，取值范圍和敏感度［36-37］。由于本文會引用筆者所建的網(wǎng)頁版新安江模型［38］，其敏感度［39］也列在表中。

筆者的網(wǎng)頁版新安江模型與原版三水源新安江模型的不同之處僅在于：引進雨量折算系數(shù)；省略子流域間匯流；省略壤中流和地下水的河網(wǎng)匯流。引進雨量折算系數(shù)是考慮多雨地區(qū)雨量誤差有可能導(dǎo)致蒸散發(fā)能力折算系數(shù)出現(xiàn)不合理數(shù)值［40］。省略子流域間匯流和省略壤中流和地下水的河網(wǎng)匯流是因為筆者的網(wǎng)頁版新安江模型主要目的是教學(xué)和研究，對象流域較小，一般不需要分子流域。而且河網(wǎng)匯流的時間尺度一般遠(yuǎn)小于壤中流和地下水的時間尺度，可以忽略。另一方面，連結(jié)各子流域的河網(wǎng)匯流無論在信息提取和模型化都已經(jīng)非常成熟，有很多水力學(xué)或水文學(xué)匯流模型可供選擇。劃分子流域在大流域的建模中帶來的不確定性相對較小。從教學(xué)和研究的觀點看，理解中小流域內(nèi)的水文規(guī)律相對比較重要。

兩個模型有14 或15 個模型參數(shù)。雖然這些參數(shù)都有明確的物理意義，參數(shù)率定即使對專業(yè)的水文工作者來說也非易事。在1999年舉行的全國水文預(yù)報技術(shù)競賽中，新安江模型的14 個隊在同一個流域用同一資料得到的參數(shù)和模擬結(jié)果還是存在較大的差異［41］。專業(yè)人員尚且如此，更何況跨專業(yè)模型用戶。怎么降低參數(shù)率定的不確定性是一個亟待解決的重要課題。

理解參數(shù)空間的結(jié)構(gòu)對參數(shù)率定，無論是手動還是自動率定，都是非常重要的。所謂結(jié)構(gòu)是指參數(shù)的取值范圍，參數(shù)的敏感度及參數(shù)間關(guān)聯(lián)性。趙人俊等［1，42］通過大量的模型率定歸納出表1所示的參數(shù)范圍，并發(fā)現(xiàn)4 個層次間的參數(shù)獨立性較好［43］。在參數(shù)敏感度方面，發(fā)現(xiàn)敏感參數(shù)只有7 個，分別為：K（蒸散發(fā)能力折算系數(shù)）、SM（表層土自由水蓄水容量）、KI（表層自由水壤中流日出流系數(shù)）、KG（表層自由水地下水日出流系數(shù)）、CG（地下水日消退系數(shù)）、L（河網(wǎng)蓄水滯后時間，表中未列出）和CS（地表水日消退系數(shù)）。

表1 三水源新安江模型參數(shù)及其取值范圍和敏感度

敏感度分析可以幫助理解參數(shù)空間的結(jié)構(gòu)。由于水文現(xiàn)象的非線性性很強，應(yīng)該盡可能使用反映參數(shù)空間整體情況的全局敏感度分析方法。Zhang 等［37］用普適似然度（GLUE）方法對新安江模型參數(shù)敏感度進行分析，指出對模型輸出最為敏感的3 個參數(shù)為K、SM和CS。K控制水量極為敏感，SM和CS影響流量過程線非常敏感。Lu 和Li［39］以實測及計算日流量計算年月日過程線的Nash 效率系數(shù)為目標(biāo)函數(shù)，在年月日這3 個不同的時間尺度上對新安江模型參數(shù)用Morris 法進行了全局敏感度分析。結(jié)果表明，在不同的時間尺度不同的參數(shù)組合顯示其敏感度。在所有時間尺度上，Cp和Cep具有壓倒性敏感度。在固定Cp和Cep后，分水源和匯流參數(shù)，SM、KG、CS和CI（壤中流消退系數(shù)）在日尺度上比較敏感。

根據(jù)對參數(shù)空間結(jié)構(gòu)的理解，就可以設(shè)計一套參數(shù)率定方案。趙人俊等的方案［42］是：以水量平衡為目標(biāo)函數(shù)率定蒸發(fā)及產(chǎn)流層的參數(shù)，以過程線絕對誤差為目標(biāo)函數(shù)率定匯流層的參數(shù)，以過程線對數(shù)絕對誤差為目標(biāo)函數(shù)率定分水源層的參數(shù)。而Li 和Lu［44］則把參數(shù)分成3 組：第一組包括兩個折算系數(shù)，第二組包括分水源及匯流的參數(shù)，第三組包括產(chǎn)流和蒸發(fā)的參數(shù)，然后依次優(yōu)化這些參數(shù)。第一組和第三組擬合年過程線，第2 組擬合日過程線。其和趙人俊等的方案有不謀而合之處，具體見下節(jié)。

2.4 模型參數(shù)自動率定自從概念性水文模型問世以來，參數(shù)自動率定一直是模型開發(fā)者和用戶的愿望。也被認(rèn)為是增加模型可用性和易用性的重要途徑。大多數(shù)自動率定以優(yōu)化算法為基礎(chǔ)。隨著計算機性能提高，涌現(xiàn)出一批出色的優(yōu)化算法，如遺傳算法、粒子群算法和模擬退火算法等等。正如Duan 等［45］指出的那樣，即使是一僅有6 個參數(shù)的簡化水文模型，其目標(biāo)函數(shù)的響應(yīng)曲面已非常復(fù)雜，有很多局部極值。這種復(fù)雜性的主要原因包括參數(shù)間相互作用， 響應(yīng)曲面的非凸性和微分的不連續(xù)性。因此，在早期的模型優(yōu)化研究中，較多使用Rosenbrock 法和單純形法等直接尋優(yōu)的方法，而不太使用基于梯度的方法。從Duan 等的結(jié)果不難知道這些局部優(yōu)化方法很難搜到全局最優(yōu)。也就是說這些方法可以或多或少改善模型參數(shù)，卻很難找到最好的參數(shù)組合。最近的參數(shù)優(yōu)化研究較多使用全局優(yōu)化算法［46-49］，也顯示了全局優(yōu)化算法的優(yōu)越性。特別是SCE-UA（或SCEM-UA）是一個考慮水文模型特性的全局優(yōu)化算法，在國內(nèi)外得到廣泛利用［50-51］。

毫無疑問，相關(guān)研究對新安江模型自動率定是非常有意義的。但必須指出，即使是全局優(yōu)化算法，也只是提高探索得到全局最優(yōu)的概率。不應(yīng)過分仰賴優(yōu)化算法，而應(yīng)該對優(yōu)化結(jié)果進行充分的論證。筆者曾經(jīng)用SCEM-UA 法對表1中15 個參數(shù)進行優(yōu)化，雖然每次優(yōu)化都能得到一組不錯的參數(shù)組合，但每次都不同。而且每次調(diào)用模型50 萬次都達(dá)不到收斂標(biāo)準(zhǔn)。Li 和Lu［44］根據(jù)敏感度分析結(jié)果，提出了多步優(yōu)化方案。把參數(shù)分成3 組：第一組包括兩個折算系數(shù)；第二組包括分水源及匯流的參數(shù)；第三組包括產(chǎn)流和蒸發(fā)的參數(shù)。從一組常用參數(shù)值出發(fā)，依次用SCEM-UA 法優(yōu)化第一組，第二組和第三組參數(shù)。前面的優(yōu)化結(jié)果在后續(xù)優(yōu)化中使用。第一組和第三組最大化年過程線效率系數(shù)，第2 組擬合最大化日過程線效率系數(shù)。為了確認(rèn)方案的收斂性，筆者也曾經(jīng)循環(huán)這三組參數(shù)的優(yōu)化，第二次、第三次結(jié)果變化不大，也就是說不用循環(huán)就可以得到較好的結(jié)果。在理想條件下，即使用人工生成流量資料時，方案基本能夠得到生成流量時使用的參數(shù)。在使用實際資料時，也能得到比較滿意的結(jié)果，包括精度和參數(shù)的合理性。收斂速度第一組最快，第二組次之，第三組最慢。模型調(diào)用次數(shù)一般在10 萬到15 萬之間，遠(yuǎn)遠(yuǎn)低于設(shè)定的優(yōu)化程序模型調(diào)用上限50 萬次。

根據(jù)Li 和Lu［44］的結(jié)果，Lu 和Li［52］提出了“降維縮容”的模型優(yōu)化戰(zhàn)略。“降維”是通過理解參數(shù)空間的結(jié)構(gòu)把高維空間分割成幾個低維子空間。在這些子空間上進行優(yōu)化，可以提高優(yōu)化的收斂性和得到全局最優(yōu)的可能性?！翱s容”分兩個部分。一個是在不失一般性的情況下盡可能縮小參數(shù)取值范圍，使參數(shù)空間容積縮小。如利用退水曲線推估地下水消退系數(shù)，然后據(jù)此設(shè)定比較小的取值范圍。另一個是考慮參數(shù)間的相互關(guān)系限制取值范圍，如設(shè)定CI ＜CG。雖然這會限制理論上可能出現(xiàn)而實際上很少出現(xiàn)的壤中流出流慢于地下水的現(xiàn)象，但卻可以縮小參數(shù)空間，并避免新安江模型結(jié)構(gòu)上潛在的異參同效［53］。從模型結(jié)構(gòu)不難發(fā)現(xiàn)，兩組參數(shù)值和互換不會對模型結(jié)果產(chǎn)生任何影響。因此在CI和CG的取值范圍有重疊時，參數(shù)空間內(nèi)會出現(xiàn)完全一模一樣的兩個子空間，使絕大多數(shù)優(yōu)化算法面臨困難。

另一個在文獻中很少被提到的是模型初始值問題。優(yōu)化過程中，每次調(diào)用模型都必須設(shè)定初始值。他們對模型輸出的影響會持續(xù)一段時間。地表水短一些，幾小時到幾天，主要因流域大小而異。壤中流較長，從幾天到幾十天。地下水最長，為數(shù)月到數(shù)年［36］。土壤水分則在數(shù)月到一年左右，氣候越干燥影響時間越長［54］。計算目標(biāo)函數(shù)時應(yīng)盡可能剔除這段時間的模型結(jié)果，以避免優(yōu)化參數(shù)值受影響。優(yōu)化次洪模型時更應(yīng)充分考慮初始值問題。

2.5 模型參數(shù)可移植性參數(shù)可移植性一方面是模型運用上的需求，也是對模型普適性的檢驗。日模型參數(shù)移植到小時進行次洪預(yù)報需要時域上的可移植性，資料短缺流域或無資料流域的洪水預(yù)報需要空域上的可移植性?？蓮臅r間，空間兩個方面討論［55］。

趙人俊等［42］指出新安江模型分水源和匯流相關(guān)參數(shù)受模型時段長影響。分水源參數(shù)受雨量時段內(nèi)均化影響比較復(fù)雜，而匯流相關(guān)參數(shù)可以進行時段轉(zhuǎn)換。Jie 等［56］對這些參數(shù)和模型時段長作了深入的分析，并建立了參數(shù)值和時段長的經(jīng)驗關(guān)系。陸旻皎［57］根據(jù)線性水庫理論得出：

式中：CΔt為時段長Δt的消退系數(shù)；T為線性水庫的蓄泄系數(shù)。并據(jù)此得到消退系數(shù)時段轉(zhuǎn)換公式：

舉例來說日消退系數(shù)是時消退系數(shù)的24 次方。因此在討論消退系數(shù)時必須明確說明時段長度。

理論上講，式（8）可以用來進行任何時段之間的轉(zhuǎn)換。但在長時段向短時段轉(zhuǎn)換時，長時段消退系數(shù)往往有很大誤差，不實用。陸旻皎［57］考慮雨量的時段內(nèi)不均勻性，簡單討論消退系數(shù)所具有的不確定性，建議率定參數(shù)或進行時段轉(zhuǎn)換時確保轉(zhuǎn)換前時段長，也就是說CΔt＞0.82 以保證退水段有大約80%洪量。趙人俊和王佩蘭［36］指出地下水日消退系數(shù)一般為0.98 ～0.998，壤中流日消退系數(shù)經(jīng)常達(dá)到0.9，可以用式（8）從日消退系數(shù)轉(zhuǎn)換得到時消退系數(shù)進行次洪計算。而河網(wǎng)蓄水消退系數(shù)決定于河網(wǎng)地貌等，匯流時間在中小流域一般為數(shù)十小時，日消退系數(shù)會有很大的不確定性。通過轉(zhuǎn)換得到時消退系數(shù)很困難。往往需要用次洪模型進行率定?？偠灾?，除了不同時段長度參數(shù)間的關(guān)系外，還必須考慮轉(zhuǎn)換前參數(shù)值所帶有的不確定性。

空域上的可移植性主要指模型參數(shù)的跨流域使用，與參數(shù)區(qū)域規(guī)律的內(nèi)涵有很大重疊。對解決缺資料或無資料地區(qū)水問題具有重要意義。

歷代水文學(xué)家都非常重視這個問題，做了大量的工作。趙人俊等［42］，陳志明［58］詳盡分析了SM與地質(zhì)條件，巖性的關(guān)系，為無資料地區(qū)確定SM推求提供了一個途徑。趙人?。?9］在1991年提出了CS隨時間、流量變化的方程。徐倩等［60］對黃山地區(qū)13 個流域進行參數(shù)率定，得出了CS與流域面積的經(jīng)驗關(guān)系。陸旻皎［57］從線性水庫理論出發(fā)得到式（7）所示消退系數(shù)和蓄泄系數(shù)T的關(guān)系，并提出蓄泄系數(shù)與流域面積的經(jīng)驗關(guān)系。在2.2 小節(jié)提及的研究也為從數(shù)值地理信息提取蓄水容量曲線提供了思路。另一方面夏自強［61］提出流域平均蓄水容量與干旱指數(shù)有關(guān)，顯示了氣候條件的影響。

Li 和Lu［62］結(jié)合年水量平衡方程式和Budyko 假設(shè)，提出了用年雨量、徑流量及蒸發(fā)能力推算Cep的關(guān)系如下：

式中：R為年徑流量；Pg為年雨量；ζ為資料干旱指數(shù)，即資料年平均蒸發(fā)量與年平均雨量之比；α為常數(shù)。式（9）在優(yōu)化時可以用來減少優(yōu)化參數(shù)。如果雨量不需要折算，即Cp =1 時，

就可以直接估算網(wǎng)頁版新安江模型的Cep，即原版新安江模型的K。這是一個很重要的參數(shù)。

利用各種水文、氣候、地理和地質(zhì)等方面的數(shù)據(jù)，信息、知識、規(guī)律和法則推求模型參數(shù)是比較理想的路子。更多是用有資料流域的參數(shù)值和各種水文，氣候，地理，地質(zhì)等方面的數(shù)據(jù)建立統(tǒng)計關(guān)系。如前所述，這方面新安江模型在應(yīng)用流域數(shù)上有絕對優(yōu)勢。每個流域的模型參數(shù)受率定者的主觀因素影響也許有很大的差異，但覆蓋全國的大數(shù)據(jù)一定會揭示中國境內(nèi)模型參數(shù)的規(guī)律性。中國幅員廣闊，覆蓋了很多氣候區(qū)，這些規(guī)律性應(yīng)該不僅僅適用于中國，對模型在全世界范圍的拓展也會有舉足輕重的影響。

3 新安江模型的未來展望

通過幾個重要方面回顧新安江模型相關(guān)研究，不難看出經(jīng)過廣泛的應(yīng)用，已經(jīng)積累了大量的實例，經(jīng)驗和知識，為新安江模型的新的發(fā)展和突破奠定了良好的基礎(chǔ)。

水文模型從根本上來說是水文規(guī)律認(rèn)知的數(shù)理表現(xiàn)。沒有正確的認(rèn)知，就沒有正確的水文模型。毫無疑問，水文規(guī)律研究是正確認(rèn)知的根本。從水文現(xiàn)象的觀察、觀測、分析到抽象化，模型化，再到模型的應(yīng)用和檢驗是一個周而復(fù)始，螺旋上升的過程。新資訊，如數(shù)字地理信息、遙測遙感資料，新技術(shù)都會提供新的視角和視野，允許我們更深入理解水文規(guī)律，最后反映到認(rèn)知和模型中去。

水文規(guī)律研究當(dāng)然也應(yīng)該包括相關(guān)領(lǐng)域。舉例來說，流域產(chǎn)沙［63-64］及化學(xué)物質(zhì)輸移［65］就和地表徑流的產(chǎn)流機制有密切關(guān)系；流域內(nèi)生態(tài)環(huán)境［66］也與濕潤狀況、植被和土壤等息息相關(guān)。在這些領(lǐng)域的應(yīng)用可以給水文模型提出新的需求，而來自相關(guān)領(lǐng)域的信息可以間接或直接地檢驗水文模型的合理性。受條件限制，水文學(xué)家經(jīng)常僅能用流量的擬合程度來判斷模型的優(yōu)劣。流量以外的檢驗手段因此具有特殊的重要意義。

從模型運用的角度來看，與其他現(xiàn)代水文模型相同，新安江模型有十多個參數(shù)而且相互影響，即使是水文工作者參數(shù)率定也非易事。依靠高速計算機和全局優(yōu)化算法盲目優(yōu)選參數(shù)也很難保證結(jié)果的最優(yōu)性。如何正確高效率定參數(shù)就成為模型應(yīng)用的一個門檻。無論是參數(shù)人工調(diào)試還是自動率定，理解參數(shù)空間的結(jié)構(gòu)都是必不可少的。趙人俊等的方案［42］和Li 和Lu 提出的多步優(yōu)化方案［44］都符合Lu 和Li［52］提出的“降維縮容”模型優(yōu)化戰(zhàn)略。通過設(shè)計不同的目標(biāo)函數(shù)進行全局敏感度分析可以從不同角度分析參數(shù)空間結(jié)構(gòu)，據(jù)此提出有效的參數(shù)率定方案應(yīng)該是有效的途徑。

另一方面，深化模型參數(shù)規(guī)律研究，減少需率定參數(shù)也是應(yīng)該努力的方向。利用各種水文、氣候、地理和地質(zhì)等方面的數(shù)據(jù)、信息、知識、規(guī)律和法則，從模型物理概念出發(fā)理論上推估模型參數(shù)值或參數(shù)間關(guān)系是一個有效的途徑。收集模型參數(shù)，形成一個覆蓋全國的大數(shù)據(jù)，分析其區(qū)域規(guī)律，也是很有價值的工作。中國幅員廣闊，覆蓋了多個氣候區(qū)，這些規(guī)律性不但適用于國內(nèi)，應(yīng)該可以適用于國外類似地區(qū)?？蔀槿辟Y料或無資料流域的模型應(yīng)用創(chuàng)造條件。

模型研發(fā)的最終目的就是效果良好且使用方便的模型。作為一個嘗試，筆者開發(fā)了一個網(wǎng)頁版新安江模型，用戶只需注冊，上傳數(shù)據(jù)，就可以開始使用模型。一般在十幾分鐘內(nèi)可以完成。不需要學(xué)習(xí)操作系統(tǒng)，不需要編程，大大降低了新安江模型使用的門檻。只要輸入?yún)?shù)值，就可以執(zhí)行模型。通過圖示模擬結(jié)果和標(biāo)示各種數(shù)值指標(biāo)，使模型參數(shù)人工調(diào)試更為順暢。作為參數(shù)率定支持功能，可估算部分參數(shù)值，推薦用戶試用。某種意義上講，我們的網(wǎng)頁版新安江模型是一個反映和測試最新研究成果的一個平臺。

在模型推廣方面，目前還看不到一個官方渠道發(fā)布模型的軟件、網(wǎng)頁和各種說明書的官方版本。要在國外推廣，英文說明書必不可少。除了說明書，設(shè)計一個好的新安江模型課件也會幫助新安江模型走進各個學(xué)校的課堂。一方面可以加強模型程序的中央控制，另一方面也可以提高模型公信力。同時減輕用戶升級模型的負(fù)擔(dān)。另外，學(xué)術(shù)界的模型應(yīng)用已經(jīng)越來越多樣化，不僅僅局限于傳統(tǒng)的洪水預(yù)報及水資源管理。迄今為止，已有很多水文研究人員做了各方面的工作。應(yīng)該有一個平臺把他們的成果整合到模型中去，形成一個能覆蓋很多與水有關(guān)的應(yīng)用領(lǐng)域的綜合性模型。