流域水文模型屬于數(shù)學模型,它在水旱災害防治和水資源開發(fā)利用中有廣泛應用,也是研究流域產沙過程和污染物質在水中輸移過程的必要前提。水文模型大致有以下四種分類方法;是按研究的范圍分類;二是按研究的對象分類;三是按研究的內容分類;四是按建模的技術分類。
    流域水文模型的研究大約始于上世紀50年代,70年代至80年代中期是其蓬勃發(fā)展時期,一些比較著名的模型,如美國的斯坦福Stanford模型和薩克拉門托Sscrsment模型,日本的水箱Tank模型以及中國的新安江模型等都在這一時期提出并得到應用。80年代后期至今,全世界范圍內流域水文模型的發(fā)展處于緩慢階段,除了對原有模型有一些修改外,幾乎沒有什么突破性的進展,也沒有有影響的新模型出現(xiàn)。
    1傳統(tǒng)模型的局限
    1.1預報效果不理想
    現(xiàn)有流域水文模型都是概念性模型,水文現(xiàn)象十分復雜,使得人們至今還不能用數(shù)學物理方程嚴格地描述其中每一個子過程。因此,現(xiàn)有流域水文模型在許多結構環(huán)節(jié)上仍主要借助于概念性元素模擬或經驗函數(shù)關系描述。這樣的模擬往往只涉及現(xiàn)象的表面而不涉及現(xiàn)象的本質或物理機制,因此,現(xiàn)有流域水文模型包含的許多參數(shù)都缺乏明確的物理意義,只能主要依據實測降雨和徑流資料來反求,而這樣求得的模型參數(shù)必然帶有經驗統(tǒng)計性,只能反映有關影響因素對流域徑流形成過程的平均作用。這是現(xiàn)有流域水文模型擬合一組資料中的大部分雖可達到令人滿意的程度,但對該組資料中的個別特殊情況,或者該組資料以外的另一些資料卻不一定能獲得令人滿意的擬合結果的癥結所在。
    1.2用最優(yōu)化方法確定模型參數(shù)對實測降雨徑流資料的依賴性很大
    現(xiàn)有流域水文模型一般都包含有兩個以上的參數(shù)需要由實測降雨和徑流資料來反求。在數(shù)學上處理這類“逆問題”的方法通常是最優(yōu)化方法,它首先要求設計一個目標函數(shù),然后在一定約束條件或無約束條件下尋求一組模型參數(shù)使該目標函數(shù)達到極值,F(xiàn)有流域水文模型目標函數(shù)通常根據模型的狀態(tài)變量或輸出變量的模擬值與實測值之間的誤差來構造,一般與模型本身的結構沒有什么聯(lián)系,所考慮的約束條件一般也與模型本身的結構沒有什么聯(lián)系,因此,按此方法求得的模型參數(shù)只能反映模擬值與實測值擬合程度,而不能揭示參數(shù)的物理意義。此外,由于模型中各參數(shù)之間可能存在相依性,以及所構造的目標函數(shù)的非單峰谷性,按最優(yōu)化方法求得的最優(yōu)參數(shù)組可能不是唯一的。
    1.3模型輸入的空間分散性和不均勻性
    流域水文模型的輸入是流域上各點的降雨過程,輸出是流域出口斷面的流量或水位過程,因此,它是一種輸入具有分散性和輸出具有集中性的模型,F(xiàn)有流域水文模型在結構上一般與此并不匹配,在實際應用中考慮這一問題時,幾乎無一例外地采用將全流域按雨量站劃分成若干個單元面積的方法,認為當單元面積的尺度小到一定程度時,即可作集中輸入和集中輸出的流域水文模型來模擬該單元面積的徑流形成,最后將各單元面積對全流域出口斷面輸出的貢獻迭加起來作為其出口斷面的輸出。顯然,這種實用的處理方法仍是不完善的。
    1.4有效的計算方法問題
    與計算機技術發(fā)展不可分割的流域水文模型,總是以離散方式來模擬流域徑流形成過程,因此采用有效的計算方法是十分重要的。流域產、匯流過程是一個連續(xù)過程,當用數(shù)學模型模擬時,必然因離散帶來離散誤差。計算機的數(shù)字有效位數(shù)總量有限,這又必然帶來誤差。如果離散化方法不適當,計算方法不夠好,則非物理的數(shù)值效應可能會導致背離物理圖景的計算結果。這就表明,流域水文模型是否成功,除了合理的結構和參數(shù)外,還有賴于采用合適的離散化格式和有效的計算方法,換言之,流域水文模擬存在數(shù)值計算是否滿足重系數(shù):
    Wj=〔WW∧W 〕 j表示預報時段數(shù),
    通過加權后可得到總輸出Q
    Qj=WX+§=WX+§ 式中θ為閥值。
    2.2黑箱模型預報原理
    在實際應用中,影響一場洪水大小的因素,可歸納為降雨量、當前流域斷面出口流量、蒸散發(fā)能力三個方面。在黑箱模型中,可將三個影響因素作為模型的輸入,計算流域出口的產流量Q。
    現(xiàn)假設有某流域,需要預報該流域t時刻后v個時段(一般采用一小時作為一個時段)的出流過程,且流域面用泰森多邊形劃分為f?zhèn)單元面積,同時具有t時刻的降雨量m、蒸發(fā)量e斷面出流量r。
    設:M=〔MM∧M〕其中m為各單元降雨總量
    E=ee∧e其中e為各單元蒸發(fā)總量R=rk其中r為流域斷面總出流量根據黑箱數(shù)學模型可得到輸入量:
    X=MER=mmKmee∧er
    假設權重系數(shù):
    Wv=
    其中θ為計算閥值
    其中每一行系數(shù),表示一個預報時段的權重系數(shù)。
    則流域出口斷面t時刻出流過程為:Qv==WvX
    2.2 權重系數(shù)的求解方法
    權重系數(shù)的求解過程,實際上是記錄洪水漲落規(guī)律的過程,這些規(guī)律的表現(xiàn)形式就是權重系數(shù)。對于一個流域的權重系數(shù)求解,必須使用該流域的多次洪水進行求解,一般預報時段數(shù)就是需要使用的洪水場數(shù)。
    為清晰分析權重系數(shù)的求解過程,現(xiàn)假設流域有n場洪水資料,其中包括降雨M、蒸發(fā)E、流域斷面總出流量Q。且已知輸入向量
    X=MER=mmKmee∧er
    Qv=
    求解權重系數(shù)
    Wv=其中θ為計算閥值
    根據黑箱數(shù)學模型
    Qv=X Wv 可得到j個線性方程組
    其中
    j=1,2,3, ┄,v;
    k=經驗常數(shù)(一般k=3);
    f=流域單元數(shù);
    v=預報時段數(shù);
    n=洪水資料場次數(shù)其中v+1<n;
    4 結語
    洪水預報的黑箱模型還不成熟,但對大量無實測資料的中小水庫進行洪水預報有作較大的意義,可以一勞永逸的解決洪水預報模型問題,為標準化的編制預報軟件打下理論基礎,是很有價值的一種思路。
    參考文獻
    趙人俊,流域水文模擬,水利電力出版社,1984
    苑希民 李鴻雁 劉樹坤 崔廣濤,神經網絡和遺傳算法在水科學領域的應用,2002
    朱華,水情自動測報系統(tǒng),1992