新安江流域水文模型

1、新安江流域水文模型第二章 新安江流域水文模型60年代初，河海大學(xué)(原華東水利學(xué)院)水文系趙人授等開始研究蓄滿產(chǎn)流模型，配合一定的匯流計算，將模型應(yīng)用于水文預(yù)報和水文設(shè)計。1973年，他們在對新安江水庫做人庫流量預(yù)報的工作中，把他們的經(jīng)驗歸納成一個完整的降雨徑流流域模型新安江模型。模型可用于濕潤地區(qū)和半濕潤地區(qū)的濕潤季節(jié)徑流模擬和計算。最初的新安江模型為兩水源模型，只能模擬地表徑流和地下徑流。80年代初期，模型研制者將薩克拉門托模型與水箱模型中，用線性水庫函數(shù)劃分水源的概念引入新安江模型，提出了三水源新安江模型，模型可以模擬地面徑流、壤中流、地下徑流。1984至1986年，又提出了四水源新安江模

2、型，可以模擬地面徑流、壤中流、快速地下徑流和慢速地下徑流。三水源新安江模型一般應(yīng)用效果較好，但模擬地下水豐富地區(qū)的日徑流過程精度不夠理想。在新安江三模型中增加慢速地下水結(jié)構(gòu)就成為四水源新安江模型。 當(dāng)流域面積較小時，新安江模型采用集總模型，當(dāng)面積較大時，采用分塊模型。分塊模型把流域分成許多塊單元流域，對每個單元流域做產(chǎn)、匯計算，得到單元流域的出口流量過程。再進行出口以下的河道洪水演算，求得流域出口的流量過程。把每個單元流域的出流過程相加，就求得了流域出口的總出流過程。劃分單元流域的主要目的是處理降雨分布的不均勻性，因此單元流域應(yīng)當(dāng)大小適當(dāng)，使得每塊面積上的降雨分布比較均勻并有一定數(shù)目的雨量站。

3、其次盡可能使單元流域與自然流域相一致，以便于分析與處理問題，并便于利用已有的小流域水文資料。如果流域內(nèi)有大中型水庫，則水庫以上的集水面積即應(yīng)作為一個單元流域。因為各單元流域的產(chǎn)匯、流計算方法基本相同，以下只討論一個單元流域的情況。2.1新安江兩水源模型 1模型結(jié)構(gòu)和參數(shù) 新安江兩水源模型的產(chǎn)流子模型采用蓄滿產(chǎn)流模型，蒸發(fā)計算采用三層蒸發(fā)計算模型。利用穩(wěn)定下滲率將徑流劃分為地面徑流和地下徑流兩種水源。地面徑流采用單位線匯流，地下徑流采用一次線性水庫匯流。模型把流域面積劃分為透水面積和不透水面積兩部分，不透水面積上的降水在滿足蒸發(fā)后將直接轉(zhuǎn)化為地面徑流。透水面積上將發(fā)生下滲，下滲的水量一部分存儲于

4、土壤層，后期耗于蒸發(fā)；滿足了流域土壤蓄水容量后的下滲水量才能轉(zhuǎn)化為徑流。不透水面積用參數(shù)表示，它是用流域內(nèi)不透水面積占全流域面積的百分比表示的。新安江模型的輸出是流域出流過程和流域蒸散發(fā)過程，輸入則為時段降雨量、蒸發(fā)皿觀測蒸發(fā)量。 新安江兩水源模型共有9個參數(shù)，一條單位線。 流域蒸發(fā)折算系數(shù)，是流域蒸散發(fā)能力與蒸發(fā)皿蒸發(fā)量之比； C深層蒸散發(fā)系數(shù)； 不透水面積占全流域面積的百分比，； 流域平均蓄水容量(指張力水)，； 流域平均上層蓄水容量, ； 流域平均下層蓄水容量, ；蓄水容量曲線指數(shù)；穩(wěn)定下滲率；地下水消退系數(shù)；單位線。模型結(jié)構(gòu)如下圖所示：圖0 新安江兩水源模型結(jié)構(gòu)示意圖圖中方框內(nèi)標(biāo)注為狀

5、態(tài)變量，方框外標(biāo)注的是模型參數(shù)。2模型參數(shù)的調(diào)試當(dāng)模型初定后就可上機調(diào)試，在調(diào)試的過程中，應(yīng)注意各參數(shù)的徑流響應(yīng)特征，以及參數(shù)之間相依性的影響，以便調(diào)試能有目的地順利進行。以下簡介的調(diào)試技術(shù)可供參考。流域蒸發(fā)折算系數(shù)的調(diào)試，對一個具體流域來說，此參數(shù)完全靠優(yōu)選。調(diào)試此參數(shù)主要考察年徑流模擬誤差，當(dāng)年徑流模擬誤差達最小時，此值最優(yōu)，但個別參數(shù)的最優(yōu)并不表示模型總體最優(yōu)。根據(jù)現(xiàn)有經(jīng)驗，年徑流模擬誤差控制在58左右就可以了。流域平均蓄水容量的調(diào)試，一般根據(jù)實測的降雨徑流資料分析或經(jīng)驗選取后，適當(dāng)微調(diào)即可。經(jīng)驗選定后，此值調(diào)試范圍不大，比較容易確定。通常在南方濕潤地區(qū)可選80150左右，江淮一帶約在1

6、10左右,燕山東北東部地區(qū)，選150左右。流域平均上層蓄水容量，一般取1020，流域平均下層蓄水容量一般取6090。蓄水容量曲線指數(shù)的調(diào)試，表示流域蓄水容量分布的不均勻性，當(dāng)全流域蓄滿后，這個參數(shù)就不起作用了，因此，應(yīng)該選取流域沒有達到蓄滿的那些洪水點據(jù)作為調(diào)試的依據(jù)。此值一般在0.20.5之間取值。注意值與值之間有相依性，二者對徑流模擬的結(jié)果相互有影響。穩(wěn)定下滲率的調(diào)試，對水源劃分起決定性作用，但至今對的研究還不充分，目前可對作一些簡單處理，如果取為一個常數(shù)使模擬效果不好時，可考慮作為變動參數(shù)?？疾焓欠窈线m主要觀察地下徑流的模擬精度。地下水消退系數(shù)的調(diào)試，主要觀察洪水退水段的擬合精度。不透水

7、面積參數(shù)對模型擬合精度影響很小，一般取00.05。2.2新安江三水源模型一、新安江三水源模型結(jié)構(gòu) 新安江三水源模型包括4個計算環(huán)節(jié)：流域產(chǎn)流計算；徑流的劃分；蒸散發(fā)計算；匯流計算。流域的產(chǎn)流計算和蒸散發(fā)計算與新安江二水源模型相同，水源劃分則完全不同，因此匯流計算也不相同。新安江三水源模型結(jié)構(gòu)示意圖如下 二、新安江三水源模型流域產(chǎn)流計算 新安江模型產(chǎn)流部分的計算是蓄滿產(chǎn)流模式， 蓄滿產(chǎn)流指在流域包氣帶土濕滿足田間持水量以前不產(chǎn)流，所有的降雨都被土壤吸收；而在土濕達到田間持水量之后，所有的降雨（除去同期的蒸散發(fā)）都產(chǎn)流。在產(chǎn)流后，流域包氣帶土壤的下滲能力為穩(wěn)定下滲率，下滲的水分成為地下徑流和壤中流

8、，超蓄的部分成為地面徑流?？紤]到流域內(nèi)各點的蓄水容量并不相同，實際產(chǎn)流時常常是在部分面積上產(chǎn)流，新安江模型引入流域蓄水容量曲線來刻劃流域內(nèi)各點蓄水容量的不均勻性，把流域內(nèi)各點的蓄水容量概化成如圖2.3所示的一條拋物線（也可概化成其它函數(shù)形式），其方程為： （2-1）式中： 為產(chǎn)流面積；是點蓄水容量,；是蓄水容量小于等于蓄水量的全部點的面積,；是單元流域面積, ；是蓄水容量曲線指數(shù)。是流域內(nèi)點最大蓄水容量,；其與流域平均蓄水容量的關(guān)系為: （2-2）引入蓄水容量曲線后，在降雨過程中，只有在滿足了包氣帶蓄水容量的面積上才可能產(chǎn)生徑流，其余面積上不產(chǎn)生徑流?？鄢舭l(fā)后的有效降雨量中，未轉(zhuǎn)化為凈雨產(chǎn)流

9、量的水分通過下滲進入流域包氣帶土層，補充土層蓄水量。引入流域蓄水容量曲線后的產(chǎn)流計算式為：當(dāng)時為部分面積產(chǎn)流 （2-3）當(dāng)時為全流域面積產(chǎn)流 （2-4） 式中： 為有效降雨量， ；為與流域土壤初始蓄水量相應(yīng)的前期影響雨量，由下式計算 （2-5）三、新安江三水源模型水源劃分1不考慮自由水蓄水容量分布不均勻的劃分方法求得的產(chǎn)流量包括地面徑流 、壤中流和地下徑流 三部分。新安江三水源模型用一個自由水蓄水庫解決水源劃分，自由水蓄水庫有兩個出流孔，底孔為地下徑流出流孔，邊孔為壤中流出流孔。新安江模型考慮了產(chǎn)流面積（）的變化，自由水蓄水庫實際只發(fā)生在產(chǎn)流面積上，其底寬為產(chǎn)流面積，顯然它是隨時間變化的。產(chǎn)流

10、量進入水庫即在產(chǎn)流面積上產(chǎn)生的徑流深，也就是自由水蓄水庫所增加的蓄水深，當(dāng)自由水蓄水深超過其最大值時，超過部分成為地面徑流。壤中流和地下徑流按線性水庫出流，其出流系數(shù)分別為和。底孔出流量和邊孔出流量分別進入各自的水庫，并按線性水庫的退水規(guī)律流出，分別成為地下水出流和壤中流出流。模型認為蒸散發(fā)在張力水中消耗，自由水蓄水庫的水量全部為徑流。水源劃分示意圖見圖3。水源劃分計算式為：地面徑流計算式：當(dāng) （2-6）當(dāng) （2-7）壤中流計算式： （2-8）地下徑流計算式： （2-9）注意每步計算時應(yīng)將自由水蓄水庫的蓄量折算為當(dāng)前時段產(chǎn)流面積上的深度。2.考慮自由水蓄水容量分布不均勻的劃分方法流域自由水蓄水

11、容量實際上也是不均勻的，作為圖3表示的水源劃分方法的改進，人們引入流域自由水蓄水容量分布曲線來刻化自由水蓄水量的不均勻性，其線型與流域蓄水容量曲線類似，曲線見圖4，曲線方程如下式。 （2-10）式中： 為產(chǎn)流面積；是點自由水蓄水容量；是自由水蓄水容量小于等于蓄水量的全部點的面積；是單元流域面積；是自由水蓄水容量曲線指數(shù)。是流域內(nèi)點最大自由水蓄水容量；與流域平均自由水蓄水容量的關(guān)系為： （2-11）由于認為在產(chǎn)流面積上才有自由水，因此，產(chǎn)流面積上的點最大自由水蓄水容量是隨產(chǎn)流面積變化而變化的，見圖4。對當(dāng)前時段，已知，則將其代人（2-10）得 （2-12）1.0從（2-12）解出得： （2-13

12、）同樣，在產(chǎn)流面積上的自由水平均蓄水容量與該面積上的點最大自由水蓄水容量的關(guān)系為： （2-14）在劃分水源時，產(chǎn)流面積上的只有在滿足了蓄水容量的水庫面積上才可能產(chǎn)生地面徑流；其余面積上不產(chǎn)生地面徑流。進入水庫的水分扣除地面徑流量后，剩余的補充水庫蓄水量S，這時的水源劃分計算式為：地面徑流出流當(dāng)時為部分面積產(chǎn)流 （2-15）當(dāng)時為全流域面積產(chǎn)流 （2-16）式中：為計算時段初的自由水蓄水庫的蓄水量（用產(chǎn)流面積上的平均深度表示）；為計算相應(yīng)于的自由水蓄水庫前期影響入流量，由下式計算； （2-17）進入水庫的水分扣除地面徑流量后，剩余的補充水庫蓄水量S： （2-18）自由水蓄水庫的時段入流量為時段產(chǎn)

13、流量，但由（2-3）或（2-4）計算的產(chǎn)流量是用流域平均水深表示的水量，應(yīng)將折算為產(chǎn)流面積上表示的水深?？梢宰C明折算為產(chǎn)流面積上表示的水深值恰為流域有效降雨值，證明如下：為折算為產(chǎn)流面積上表示的水深，根據(jù)降雨徑流關(guān)系的性質(zhì)，水文學(xué)原理中已經(jīng)證明相對產(chǎn)流面積與產(chǎn)流量之間有關(guān)系： （2-19）從此式可看出壤中流出流： （2-20）地下徑流出流： （2-21）由于，都是產(chǎn)流面積上的深度表示的水量，計算所得、應(yīng)折算回全流域深度表示的水量，即：、在當(dāng)前時段計算結(jié)束后應(yīng)乘。蒸散發(fā)計算四、流域蒸散發(fā)計算蒸散發(fā)模型不考慮面上分布的不均勻性，但可考慮土濕垂向分布的不均勻性，根據(jù)需要采用兩層或三層蒸散發(fā)模型計算流

14、域蒸發(fā)量。兩層蒸散發(fā)模型將土層分為上、下兩層，各層蓄水容量分別為、（）。流域土層下滲蓄水和蒸散發(fā)計算過程按下述原則進行：降雨先補充上層，上層蓄滿后再補充下層；蒸發(fā)時先蒸發(fā)上層，上層蓄水量蒸發(fā)殆盡后再蒸發(fā)下層。1.兩層蒸散發(fā)模型的計算式如下：， ： ，式中：是流域降雨量，；是流域蒸散發(fā)能力，；是上層土壤蓄水量，；是下層土壤蓄水量，； 是下層土壤蓄水容量，；、分別為上層、下層蒸發(fā)量，。2.三層蒸散發(fā)模型的計算式如下：當(dāng)：， ，時：，時：，時：，式中：是深層蒸散發(fā)系數(shù)； 五、流域匯流計算匯流計算可以采用兩種方法：1.傳統(tǒng)流域匯流這種匯流方法是：地面徑流采用單位線匯流，壤中流和地下徑流采用線性水庫匯流

15、至出口斷面后與單位線匯流結(jié)果疊加，最后得到流域出口斷面的模擬出流。即圖3中下面的兩個水庫的出流就是壤中流和地下徑流在流域出口斷面的出流流量，該兩個水庫分別稱為壤中流調(diào)蓄水庫和地下徑流調(diào)蓄水庫。壤中流匯流計算地下徑流匯流計算式中： t時段壤中流出流流量，；壤中流消退系數(shù)，可假定為；t時段地下徑流出流流量，；地下徑流消退系數(shù)，可假定為；是徑流深轉(zhuǎn)換為流量時的單位轉(zhuǎn)換系數(shù)，； 上式中假定了調(diào)蓄水庫與自由水水庫的出流系數(shù)相同。在實踐中，為延緩壤中流和地下徑流出流，有人把下面的兩個水庫的出流看作是坡地調(diào)蓄水庫，而在它們后面再分別串聯(lián)一個水庫反映壤中流和地下徑流受河槽調(diào)蓄的作用。2.河網(wǎng)匯流單位線匯流河網(wǎng)

16、匯流單位線也稱時變河網(wǎng)匯流單位線，這種匯流方法忽略地面徑流的坡面匯流階段，直接將地面徑流產(chǎn)流量轉(zhuǎn)換為流量作為河網(wǎng)地面徑流入流（看作集中在計算時段末入流）。地下徑流和壤中流匯流均采用線性水庫作坡地調(diào)蓄后得到相應(yīng)河網(wǎng)入流，即：把圖3中的轉(zhuǎn)換為流量作為河網(wǎng)地面徑流入流，下面的兩個水庫的出流看作壤中流和地下徑流經(jīng)坡地調(diào)蓄后注入河網(wǎng)的入流量。三種水源的河網(wǎng)入流流量相加就是所謂的河網(wǎng)總?cè)肓?。河網(wǎng)總?cè)肓鹘?jīng)流域河網(wǎng)調(diào)蓄后得到流域出口斷面的流量過程。河網(wǎng)匯流單位線的計算式如下： （2-22）式中：t時段河網(wǎng)的總?cè)肓?，；河網(wǎng)消退系數(shù)；河網(wǎng)匯流時間，；t時刻的河網(wǎng)出流（即t時刻的流域出口斷面流量），；t時刻的河網(wǎng)總

17、入流經(jīng)河網(wǎng)調(diào)蓄后，延遲時間的出流，。是隨時間而變化的時變參數(shù)，它取決于t時刻的河網(wǎng)的總?cè)肓髁亢土饔虻暮硬厶匦?。一般用下式計算?（2-23）式中：反映流域河槽特性的系數(shù)。2.3線性水庫不同單位時段的出流系數(shù)轉(zhuǎn)換 從(2-20)和(2-21)知道新安江模型的自由水水庫出流是按線性水庫出流處理的。在客觀世界中，流域的降雨、蒸散發(fā)、徑流等水文物理量都是連續(xù)量，流域內(nèi)的蓄水量也是連續(xù)量，觀測和計算時都只能處理為一個一個的時段平均量；從數(shù)學(xué)觀點來看，這實際上是將連續(xù)量處理為離散量，將表達物理過程的連續(xù)方程處理為表達離散量的差分方程，這樣就必然帶來差分誤差，特別是水文模型計算通常采用的時段較大，如果在模型

18、的算法設(shè)計上不注意，差分誤差可能較大。為了減小計算的差分誤差，本可以縮小計算時段，但是，在水文上常常受到限制，這是由于原始水文觀測數(shù)據(jù)時段本身較大，如果在所有的計算環(huán)節(jié)上都縮小計算時段并不能提高成果的精度。盡管如此，在模型的算法上，至少應(yīng)要求算法不會帶來更大的差分誤差。仔細考察新安江流域水文模型計算的各個環(huán)節(jié)，就會發(fā)現(xiàn)(2-20)和(2-21)的自由水水庫出流計算可能會出現(xiàn)較大的差分誤差，為此，有的模型程序設(shè)計者在算法上作了一定處理，以設(shè)法減小差分誤差?；镜南敕ㄊ窃谶@個計算環(huán)節(jié)上減小計算步長，即將模型總體的計算時段在這個環(huán)節(jié)上進一步細分；為此，需要將水庫長時段的出流系數(shù)轉(zhuǎn)換為短時段的出流系數(shù)

19、。設(shè)水庫在一個計算時段內(nèi)的入流量為，為減小差分誤差，將等分為份入流，這實際上是將原計算時段細分為個更小的計算時段，每個時段的水庫入流為?？梢宰C明水庫計算時段為的出流系數(shù)和與計算時段為的出流系數(shù)和之間存如下關(guān)系： ； （2-24）證明如下：1. 只有一個出流孔的線性水庫不同時段出流系數(shù)的轉(zhuǎn)換 先對只有一個出流孔的線性水庫進行證明。設(shè)線性水庫的出流方程和退水時的水量平衡方程如下： （2-25） （2-25）式中： 為線性水庫的瞬時出流系數(shù)；為線性水庫時刻的蓄量；為線性水庫時刻的出流流量。合并（2-25）和（2-26）得到 （2-27）在上積分上式得 （2-28）此式稱為線性水庫退水方程，稱為線性水

20、庫消退系數(shù)（又稱退水系數(shù)）。此式的另一個表達形式是： （2-29）這是（2-27）在區(qū)間上積分的結(jié)果。令，即取為一個單位計算時段，則上式記為 （2-30）在區(qū)間上積分得到水庫在時段內(nèi)流出的水量 （2-31）令，即取為一個單位計算時段，則上式記為 （2-32）（2-32）中的就是線性水庫單位時段為的出流系數(shù)。比較（2-25）和（2-32）可知和的意義是完全不同的，總是與一個單位時段相聯(lián)系的，而且單位時段不同，其值也不同。令，記單位時段為的水庫消退系數(shù)為。由（2-30）的退水方程，可推導(dǎo)得 （2-33）此式中的1代表一個時段，表示個時段。顯然。對單位時段的退水方程，有； 或 （2-34）比較（2-

21、33）最后一個式子與（2-34）式可知，線性水庫單位計算時段為的出流系數(shù)與單位計算時段為的出流系數(shù)之間有如下關(guān)系： （2-35）已知對單位計算時段，有；記單位計算時段的出流系數(shù)為，則可推得 （2-36）此式就是所需要得結(jié)果。2有二個出流孔的線性水庫不同時段出流系數(shù)的轉(zhuǎn)換記單位計算時段為的出流系數(shù)分別為和，單位計算時段為的出流系數(shù)分別為和，不同計算時段的兩孔總出流系數(shù)分別為和，代人（2-36）式得到： （2-37）要求解和還須增加一個條件，可假定不同計算時段的兩孔出流系數(shù)分配比例相同，即成立： （2-38）利用（2-38）消去得到： （2-39）利用（2-38）消去得到： （2-40）利用（2-

22、39）與（2-40）就可進行線性水庫不同時段的出流系數(shù)轉(zhuǎn)換。2.4新安江三水源模型的參數(shù)調(diào)試 采用河網(wǎng)匯流單位線的新安江三水源模型共有14個參數(shù)，根據(jù)參數(shù)的功能可分為4類 與蒸散發(fā)計算有關(guān)的參數(shù)：蒸散發(fā)折算系數(shù)，上層土壤蓄水容量下層土壤蓄水容量，深層蒸散發(fā)系數(shù)。與產(chǎn)流計算有關(guān)的參數(shù)：流域蓄水容量，流域蓄水容量曲線指數(shù)，流域不透水面積比值。與水源劃分有關(guān)的參數(shù)：流域自由水蓄水容量，流域自由水蓄水容量曲線指數(shù)，自由水蓄水庫壤中流出流系數(shù)與自由水蓄水庫地下徑流出流系數(shù)。與匯流有關(guān)的參數(shù)：壤中流消退系數(shù)和地下徑流消退系數(shù)；河網(wǎng)匯流單位線參數(shù)和。壤中流匯流調(diào)蓄水庫和地下徑流匯流調(diào)蓄水庫與自由水水庫是不同

23、的水庫，但如果假定二調(diào)蓄水庫的出流系數(shù)與自由水水庫對應(yīng)的徑流的出流系數(shù)相同，即假定，則參數(shù)可減為12個。4類參數(shù)相對獨立，同類參數(shù)有一定相依性。1參數(shù)初定根據(jù)流域自然地理資料、水文氣象資料和前人以積累的經(jīng)驗初定參數(shù)。蒸散發(fā)折算系數(shù)：可采用汛期連續(xù)大雨，流域蓄水量處于蓄滿狀態(tài)的累計降雨、累計徑流、累計蒸發(fā)觀測資料，用下式估計: （2-41）流域蓄水容量、上層土壤蓄水容量、下層土壤蓄水容量：值一般根據(jù)實測的降雨徑流資料分析或經(jīng)驗選取后，適當(dāng)微調(diào)即可。經(jīng)驗選定后，此值調(diào)試范圍不大，比較容易確定。通常在南方濕潤地區(qū)可選80150左右，北方半濕潤地區(qū)可達170，江淮一帶約在110左右,燕山東北東部地區(qū)，

24、選150左右。 中包含有植物截留，少林地510，多林地1020左右。在6090范圍。 流域蓄水容量曲線指數(shù)：此值一般在0.20.5之間取值，通常流域越小，該值越小，特小流域可能小于0.2。1000平方公里以上流域一般大于0.3。 不透水面積比值：值對模型擬合精度影響很小，一般取00.05。深層蒸散發(fā)系數(shù)： 南方多林地區(qū)左右 ；北方半干旱地區(qū)左右。 流域自由水蓄水容量： 一般由經(jīng)驗初定，土層淺薄的山區(qū)，土厚多林地區(qū)，甚至更大。 流域自由水蓄水容量曲線指數(shù)在； 根據(jù)經(jīng)驗，常在11.5左右。 自由水蓄水庫壤中流出流系數(shù)與自由水蓄水庫地下徑流出流系數(shù)：一般來說，約為0.7左右，二者分配比與具體流域有關(guān)

25、，只能調(diào)試決定。壤中流消退系數(shù)和地下徑流消退系數(shù)：一般趨近0.9；約為0.980.998，具體取值與具體流域有關(guān)，可進一步調(diào)試決定。河網(wǎng)匯流單位線參數(shù)和在0.02左右，與具體流域的河網(wǎng)特性有關(guān)，可進一步調(diào)試決定。與具體流域的河網(wǎng)特性有關(guān)，為計算方便，通常取計算單位時段的整倍數(shù)。2參數(shù)率定首先率定第一類參數(shù)，主要觀察率定期的總水量是否與實測水量相近。其次率定第二類參數(shù)，主要觀察次洪水的水量是否與實測水量相近。第三類參數(shù)的率定，主要觀察次洪水的退水段過程的水量是否與實測水量相近。第四類參數(shù)的率定，主要觀察次洪水的過程線是否擬合良好。參數(shù)率定過程中要反復(fù)對參數(shù)進行多遍調(diào)試，直至滿足精度要求。2.5新

26、安江三水源模型的VB程序介紹新安江三水源，兩層蒸發(fā)模型（時變單位線）VB源程序Function XingAnJang4() As Integer Dim Fno As Integer, I As Integer，Nian(2) As Integer，AX(20) As Single, Qo As Single, Dim Wm As Single, Wum As Single, Wlm As Single, k As Single, Sm As Single, b As Single, Bx As Single Dim Kss As Single, Kg As Single, CRss As S

27、ingle, CRg As Single, Immp As Single, F As Single, U As Single Dim P As Single, E As Single, Pe As Single, Eu As Single, El As Single, EE As Single, FR0 As Single Dim R As Single, Rs As Single, Rss As Single, Rg As Single, Q As Single, X As Single, XX As Single Dim W As Single, Wu As Single, Wl As S

28、ingle, S As Single, Wmm As Single, Smm As Single Dim QRss0 As Single, QRg0 As Single, Cr As Single, Cs As Single, Cx As Single, Qz As Single, Qa() As Single Qo = 5 'Qo是基流。 Open "e:test流域模型enshi9699parm3.txt" For Input As #1 Open "e:test流域模型enshi9699ensh9699.txt" For Input As

29、#5 Open "e:test流域模型enshi9699end96994.txt" For Output As #6 Fno = 1 ReadRecord Fno:TurnRecord Nian, AX, 11 + 1, 1: Close (Fno) Wm = AX(1): Wum = AX(2): k = AX(3): b = AX(4): Sm = AX(5): Bx = AX(6): Kss = AX(7): Kg = AX(8) CRss = AX(9): CRg = AX(10): Immp = AX(11) W = 85: Wu = 1: S = 0: FR0

30、= 0: Dt = 3: F = 2928: U = F / (3.6 * Dt): QRss0 = 0: QRg0 = 10.1 Wmm = (1 + b) * Wm: Smm = (1 + Bx) * Sm: Wlm = Wm - Wum: Wl = W - Wu: X = 1 / (1 + b): XX = 1 / (1 + Bx) Cr = 0.025: Cx = 6 Fno = 5: n = 0: Do Until EOF(Fno): Line Input #Fno, TextLien: n = n + 1: Loop: Seek #Fno, 1 ReDim Qa(n + 1000)

31、 For I = 1 To n ReadRecord Fno: TurnRecord Nian, AX, 4 + 1, 1 n1 = Nian(1): m1 = Int(AX(1): P = AX(2): E = AX(3): qb = AX(4) Pe = P - k * E If Pe > 0 Then Eu = k * E: El = 0 XuManChanLiu1 Pe, R, Wm, Wmm, Wum, Wlm, W, Wu, Wl, b, X Else R = 0 EuEl P, k, E, Wlm, Wu, Wl, Eu, El End If RsRssRg1 Pe, R,

32、 S, Sm, Smm, Rs, Rss, Rg, Bx, XX, FR0, Kss, Kg, Immp 'RsRssRg2 Pe, R, S, Sm, Smm, Rs, Rss, Rg, Bx, XX, FR0, Kss, Kg, Immp EE = Eu + El: W = Wu + Wl HeWangZongRuLiu Rs, Rss, Rg, QRss0, QRg0, Qz, CRg, CRss, Qo, U Tao = 1 Cs = 1 - Cr * Qz 0.4 'Tao = Int(Cx / (1 - Cs) / 3 Qa(I + Tao) = Cs * Qa(I

33、) + (1 - Cs) * Qz Q = Qa(I + Tao) Write #6 I, P, EE, R, Rs, Rss, Rg, W, Q, qb End If Next I Close (6) Form1.Text1.Text = CStr(666)End FunctionPublic Sp(20) As VariantFunction ReadRecord(Fno) As Integer Dim TextLien As String, Stext As String, sss As String Dim I As Integer, j As Integer, LineLen As

34、Integer j = 1 Line Input #Fno, TextLien LineLen = Len(TextLien) For I = 1 To LineLen Stext = Mid(TextLien, I, 1) Select Case Asc(Stext) Case 45, 46, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57 sss = sss + Stext If I = LineLen Then Sp(j) = sss End If Case Else If sss <> "" Then Sp(j) = sss

35、: j = j + 1: sss = "" End If End Select Next I ReadRecord = jEnd Function'蓄滿產(chǎn)流函數(shù)一Function XuManChanLiu1(Pe, R, Wm, Wmm, Wum, Wlm, W, Wu, Wl, b, X) As Integer Dim Dr As Single, A As Single A = Wmm * (1 - (1 - W / Wm) X) If Pe + A < Wmm Then R = Pe - Wm + W + Wm * (1 - (Pe + A) / Wmm)

36、 (1 + b) Else R = Pe + W - Wm End If If Wu + Pe - R < Wum Then Wu = Wu + Pe - R Else Wl = Wl + Wu + Pe - R - Wum Wu = Wum If Wl > Wlm Then Dr = Wl - Wlm Wl = Wlm R = R + Dr End If End IfEnd Function'兩層蒸發(fā)函數(shù)Function EuEl(P, k, E, Wlm, Wu, Wl, Eu, El) As Integer If Wu + P > k * E Then Eu = k * E Wu = Wu + P - Eu El = 0 Else Eu = Wu + P El = (k * E - Eu) * Wl / Wlm Wu = 0 Wl = Wl - El End IfEnd Function'三水源劃分函數(shù)一Function RsRssRg1(Pe, R, S, Sm, Smm, Rs, Rss, Rg, Bx, XX, FR0, Kss, Kg, Immp) As Integer Dim j As Integer, nn As Integer Dim FR As Single, Au As Single, Smmf As Single