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

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

3、流域與自然流域相一致，以便于分析與處理問題，并便于利用已有的小流域水文資料。如果流域內有大中型水庫，則水庫以上的集水面積即應作為一個單元流域。因為各單元流域的產匯、流計算方法基本相同，以下只討論一個單元流域的情況。2.1新安江兩水源模型1 .模型結構和參數新安江兩水源模型的產流子模型采用蓄滿產流模型，蒸發計算采用三層蒸發計算模型。利用穩定下滲率FC將徑流劃分為地面徑流和地下徑流兩種水源。地面徑流采用單位線匯流，地下徑流采用一次線性水庫匯流。模型把流域面積劃分為透水面積和不透水面積兩部分，不透水面積上的降水在滿足蒸發后將直接轉化為地面徑流。透水面積上將發生下滲，下滲的水量一部分存儲于土壤層，后期

4、耗于蒸發；滿足了流域土壤蓄水容量后的下滲水量才能轉化為徑流。不透水面積用參數IMP表示，它是用流域內不透水面積占全流域面積的百分比表示的。新安江模型的輸出是流域出流過程Qt和流域蒸散發過程Et,輸入則為時段降雨量P、蒸發皿觀測蒸發量EI。新安江兩水源模型共有9個參數，一條單位線。K流域蒸發折算系數，是流域蒸散發能力與蒸發皿蒸發量之比；C深層蒸散發系數；IMP不透水面積占全流域面積的百分比，%;Wm流域平土蓄水容量（指張力水），mm；WUM流域平均上層蓄水容量,mm；WLM流域平均下層蓄水容量,mm；B蓄水容量曲線指數；FC穩定下滲率mm/h；KKG地下水消退系數；UH單位線。模型結構如下圖所示

5、：圖0新安江兩水源模型結構示意圖圖中方框內標注為狀態變量，方框外標注的是模型參數。2 .模型參數的調試當模型初定后就可上機調試，在調試的過程中，應注意各參數的徑流響應特征，以及參數之間相依性的影響，以便調試能有目的地順利進行。以下簡介的調試技術可供參考。流域蒸發折算系數K的調試，對一個具體流域來說，此參數完全靠優選。調試此參數主要考察年徑流模擬誤差，當年徑流模擬誤差達最小時，此值最優，但個別參數的最優并不表示模型總體最優。根據現有經驗，年徑流模擬誤差控制在5%8%左右就可以了。流域平均蓄水容量Wm的調試，一般根據實測的降雨徑流資料分析或經驗選取后，適當微調即可。經驗選定后，此值調試范圍不大，比

6、較容易確定。通常在南方濕潤地區Wm可選80150mm左右，江淮一帶約在110mm左右，燕山東北東部地區，選150mm左右。流域平均上層蓄水容量WUM,一般取1020mm,流域平均下層蓄水容量WLM一般取6090mm。蓄水容量曲線指數的調試，B表示流域蓄水容量分布的不均勻性，當全流域蓄滿后，這個參數就不起作用了，因此，應該選取流域沒有達到蓄滿的那些洪水點據作為調試的依據。此值一般在0.20.5之間取值。注意B值與Wm值之間有相依性，二者對徑流模擬的結果相互有影響。穩定下滲率FC的調試，FC對水源劃分起決定性作用，但至今對FC的研究還不充分，目前可對FC作一些簡單處理，如果FC取為一個常數使模擬效

7、果不好時，可考慮FC作為變動參數。考察FC是否合適主要觀察地下徑流的模擬精度。地下水消退系數KKG的調試，主要觀察洪水退水段的擬合精度。不透水面積參數IMP對模型擬合精度影響很小，一般取00.05。2.2新安江三水源模型一、新安江三水源模型結構新安江三水源模型包括4個計算環節：流域產流計算；徑流的劃分；蒸散發計算；匯流計算。流域的產流計算和蒸散發計算與新安江二水源模型相同，水源劃分則完全不同，因此匯流計算也不相同。新安江三水源模型結構示意圖如下上層蒸發EU蒸散發E深層蒸發ED不透水面積IMP降雨徑流總流量Q圖1.新安江三水源模型示意圖、新安江三水源模型流域產流計算新安江模型產流部分的計算是蓄滿

8、產流模式，蓄滿產流指在流域包氣帶土濕滿足田間持水量以前不產流，所有的降雨都被土壤吸收；而在土濕達到田間持水量之后，所有的降雨（除去同期的蒸散發）都產流。在產流后，流域包氣帶土壤的下滲能力為穩定下滲率，下滲的水分成為地下徑流和壤中流，超蓄的部分成為地面徑流。考慮到流域內各點的蓄水容量并不相同，實際產流時常常是在部分面積上產流，新安江模型引入流域蓄水容量曲線來刻劃流域內各點蓄水容量的不均勻性，把流域內各點的蓄水容量概化成如圖2.3所示的一條拋物線（也可概化成其它函數形式），其方程為：fF1一西W.mm(2-1)式中：f=fr為產流面積；F-_1t_trak-l-t_1hWm是點畜水谷量,mm;f是

9、蓄水容量小于等于蓄水量Wm的全部點的面積，Km2;F是單元流域面積，Km2;b是蓄水容量曲線指數。Wmm是流域內點最大蓄水容量，mm;其與流域平均蓄水容量Wm的關系為：Wmm=（1+b）.Wm（2-2）引入蓄水容量曲線后，在降雨過程中，只有在滿足了包氣帶蓄水容量的面積上才可能產生徑流，其余面積上不產生徑流。扣除蒸發后的有效降雨量R中，未轉化為凈雨產流量R的水分通過下滲進入流域包氣帶土層，補充土層蓄水量W。引入流域蓄水容量曲線后的產流計算式為：當巳+a<Wmm時為部分面積產流R=RWm+W0+Wm1-1Pe+a'寸WmmJ(2-3)當Pe+a之Wmm時為全流域面積產流R=Pe-Wm

10、W0(2-4)圖2:流域蓄水容量分布曲線式中：Pe為有效降雨量，Pe=P-E；a為與流域土壤初始蓄水量W0相應的前期影響雨量，由下式計,Wn、例a=Wmm.1+1(2-5)lWM)L_.三、新安江三水源模型水源劃分1 .不考慮自由水蓄水容量分布不均勻的劃分方法求得的產流量R包括地面徑流RS、壤中流RI和地下徑流RG三部分。新安江三水源模型用一個自由水蓄水庫解決水源劃分，自由水蓄水庫有兩個出流孔，底孔為地下徑流RG出流孔，邊孔為壤中流RI出流孔。新安江模型考慮了產流面積(FR=R/Pe)的變化，自由水蓄水庫實際只發生在產流面積上，其底寬為產流面積FR,顯然它是隨時間變化的。產流量R進入水庫即在產

11、流面積上產生PJ勺徑流深，也就是自由水蓄水庫所增加的蓄水深，當自由水蓄水深S超過其最大值Sm時，超過部分成為地面徑流RS。壤中流RI和地下徑流RG按線性水庫出流，其出流系數分別為KI和KGo底孔出流量RG和邊孔出流量RI分別進入各自的水庫，并按線性水庫的退水規律流出，分別成為地下水出流QG和壤中流出流QI。模型認為蒸散發在張力水中消耗，自由水蓄水庫的水量全部為徑流。水源劃分示意圖見圖3。水源劃分計算式為：圖3:自由水蓄水庫結構地面徑流計算式：當Pe+S。SmRS=Pe+So-Sm(2-6)當Pe+S0MsmRS=0(2-7)壤中流計算式：RI=KIS(2-8)地下徑流計算式：RG=KG,S(2

12、-9)注意每步計算時應將自由水蓄水庫的蓄量S折算為當前時段產流面積上的深度。2 .考慮自由水蓄水容量分布不均勻的劃分方法流域自由水蓄水容量實際上也是不均勻的，作為圖3表示的水源劃分方法的改進，人們引入流域自由水蓄水容量分布曲線來刻化自由水蓄水量的不均勻性，其線型與流域蓄水容量曲線類似，曲線見圖4,曲線方程如下式。f/Sm、似一=11-(2-10)F1smm,式中：f=FR為產流面積；FSm是點自由水蓄水容量；f是自由水蓄水容量小于等于蓄水量Sm的全部點的面積；F是單元流域面積；bx是自由水蓄水容量曲線指數。Smm是流域內點最大自由水蓄水容量；與流域平均自由水蓄水容量SM的關系為：Smm=(1十

13、bX卜Sm(2-11)由于認為在產流面積上才有自由水，因此，產流面積FR上的點最大自由水蓄水容量Smmf是隨產流面積FR變化而變化的，見圖4。對當前時段，已知FR,則將其代人(2-10)得FRbxSmmfSmm(2-12)S7mmFRFR1.0從（2-12）解出Smmf得:Smmf=Smm1-1-FRbX(2-13)同樣，在產流面積FR上的自由水平均蓄水容量Smf與該面積上的點最大自由水蓄水容量Smmf的關系為：Smmf=1bxSmf(2-14)(2-15)(2-16)在劃分水源時，產流面積上的只有在滿足了蓄水容量的水庫面積上才可能產生地面徑流RS；其余面積上不產生地面徑流。進入水庫的水分扣除

14、地面徑流量RS后，剩余的補充水庫蓄水量S,這時的水源劃分計算式為：地面徑流出流當Pe,Au：二Smmf時為部分面積產流1P+Au而RS=Pe-Smf十So十Smf1-1Smmf當EAu-Smmf時為全流域面積產流RS=Pe-SmfS0式中:So為計算時段初的自由水蓄水庫的蓄水量(用產流面積上的平均深度表小)；Au為計算相應于So的自由水蓄水庫前期影響入流量，由下式計算；1S'J4bxAu=Smmf1+1-0(2-17)SSSmfJ進入水庫的水分扣除地面徑流量RS后，剩余的補充水庫蓄水量S:S=So+PeRS(2-18)自由水蓄水庫的時段入流量為時段產流量R,但由(2-3)或(2-4)計

15、算的產流量R是用流域平均水深表示的水量，應將R折算為產流面積上表示的水深。可以證明R折算為產流面積上表示的水深值恰為流域有效降雨值P.,證明如下：RF/f為R折算為產流面積上表示的水深，根據降雨徑流關系的性質，RPe水文學原理中已經證明相對產流面積f/F與產流量R之間有關系：(2-19)從此式可看出RF/f=Pe壤中流出流:RI=KIS(2-20)地下徑流出流：RG=KG，S(2-21)由于S,R都是產流面積上的深度表示的水量，計算所得RS、RI、RG應折算回全流域深度表示的水量，即：RS、RI、RG在當前時段計算結束后應乘f/Fo蒸散發計算四、流域蒸散發計算蒸散發模型不考慮面上分布的不均勻性

16、，但可考慮土濕垂向分布的不均勻性，根據需要采用兩層或三層蒸散發模型計算流域蒸發量。兩層蒸散發模型將士層分為上、下兩層，各層蓄水容量分別為WUM、WLM(WUM=WM-WLM)。流域土層下滲蓄水和蒸散發計算過程按下述原則進行：降雨先補充上層，上層蓄滿后再補充下層；蒸發時先蒸發上層，上層蓄水量蒸發殆盡后再蒸發下層。1 .兩層蒸散發模型的計算式如下：P+WU'Ep：EU=Ep,EL=0,E=EU2 2)P+WU<Ep：EU=P+WU,EL=(Ep-EU)WL/WLM,E=EU+EL式中：P是流域降雨量，mm;Ep是流域蒸散發能力，mm;WU是上層土壤畜水量，mm;WL是下層土壤蓄水量，

17、mm;WLM是下層土壤蓄水容量，mm;EU、EL分別為上層、下層蒸發量，mm2.三層蒸散發模型的計算式如下:當P+WU之Ep：EU=Ep,EL=0,ED=0,E=EUP+WU<Ep,WL之CWLM時：EU=P+WU,EL=Ep-EUWL/WLM,ED=0,E=EUELP+WU<Ep,C(Ep-EU)<WL<CWLM時：EU=P+WU,EL=CEp-EU,ED=0,E=EUEL(4) P+WU<Ep,WL<C(Ep-EU)時：EU=P+WU,EL=CWL,ED=CEp-EU-EL,E=EUELED式中：C是深層蒸散發系數；五、流域匯流計算匯流計算可以采用兩種方

18、法：1 .傳統流域匯流這種匯流方法是：地面徑流采用單位線匯流，壤中流和地下徑流采用線性水庫匯流至出口斷面后與單位線匯流結果疊加，最后得到流域出口斷面的模擬出流。即圖3中下面的兩個水庫的出流就是壤中流和地下徑流在流域出口斷面的出流流量，該兩個水庫分別稱為壤中流調蓄水庫和地下徑流調蓄水庫。壤中流匯流計算QIt=QIt/CIRIt1-CIU地下徑流匯流計算QGt=QGCGRGt1-CGU式中：QItt時段壤中流出流流量，m3/s;CI壤中流消退系數，可假定為CI=1-KI；QGtt時段地下徑流出流流量，m3/s;CG地下徑流消退系數，可假定為CG=1-KG；U是徑流深轉換為流量時的單位轉換系數，U=

19、F/(3.64);上式中假定了調蓄水庫與自由水水庫的出流系數相同。在實踐中，為延緩壤中流和地下徑流出流，有人把下面的兩個水庫的出流看作是坡地調蓄水庫，而在它們后面再分別串聯一個水庫反映壤中流和地下徑流受河槽調蓄的作用。2 .河網匯流單位線匯流河網匯流單位線也稱時變河網匯流單位線，這種匯流方法忽略地面徑流的坡面匯流階段，直接將地面徑流產流量轉換為流量作為河網地面徑流入流QS=RSU(看作集中在計算時段末入流)。地下徑流和壤中流匯流均采用線性水庫作坡地調蓄后得到相應河網入流QI、QG,即：把圖3中的RS轉換為流量作為河網地面徑流入流，下面的兩個水庫的出流看作壤中流和地下徑流經坡地調蓄后注入河網的入

20、流量。三種水源的河網入流流量相加就是所謂的河網總入流。河網總入流經流域河網調蓄后得到流域出口斷面的流量過程。河網匯流單位線的計算式如下：Qy=CsQt+(1-Cs)Qz,t(2-22)式中：QZ,tt時段河網的總入流，QZ,t=QS+QI+QG,m3/s；Cs河網消退系數；.河網匯流時間，h；Qtt時刻的河網出流(即t時刻的流域出口斷面流量)，m3/s;Qt+t時刻的河網總入流經河網調蓄后，延遲七時間的出流，m3/s。Cs是隨時間而變化的時變參數，它取決于t時刻的河網的總入流量Qz,t和流域的河槽特性。一般用下式計算：_04，一一、Cs=1-CrQz,t(2-23)式中：Cr反映流域河槽特性的

21、系數。2.3線性水庫不同單位時段的出流系數轉換從(2-20)和(2-21)知道新安江模型的自由水水庫出流是按線性水庫出流處理的。在客觀世界中，流域的降雨、蒸散發、徑流等水文物理量都是連續量，流域內的蓄水量也是連續量，觀測和計算時都只能處理為一個一個的時段平均量；從數學觀點來看，這實際上是將連續量處理為離散量，將表達物理過程的連續方程處理為表達離散量的差分方程，這樣就必然帶來差分誤差，特別是水文模型計算通常采用的時段較大，如果在模型的算法設計上不注意，差分誤差可能較大。為了減小計算的差分誤差，本可以縮小計算時段，但是，在水文上常常受到限制，這是由于原始水文觀測數據時段本身較大，如果在所有的計算環

22、節上都縮小計算時段并不能提高成果的精度。盡管如此，在模型的算法上，至少應要求算法不會帶來更大的差分誤差。仔細考察新安江流域水文模型計算的各個環節，就會發現(2-20)和(2-21)的自由水水庫出流計算可能會出現較大的差分誤差，為此，有的模型程序設計者在算法上作了一定處理，以設法減小差分誤差。基本的想法是在這個計算環節上減小計算步長，即將模型總體的計算時段在這個環節上進一步細分；為此，需要將水庫長時段的出流系數轉換為短時段的出流系數。設水庫在一個計算時段At內的入流量為R,為減小差分誤差，將R等分為m份入流，這實際上是將原計算時段，t細分為m個更小的計算時段，每個時段的水庫入流為R/m。可以證明

23、水庫計算時段為8的出流系數KI和KG與計算時段為At/m的出流系數KI'和KG'之間存如下關系：KIKGKG11-1-KIKG】mKIKIKG(2-24)證明如下：1.只有一個出流孔的線性水庫不同時段出流系數的轉換(2-25)先對只有一個出流孔的線性水庫進行證明。設線性水庫的出流方程和退水時的水量平衡方程如下：QtStdSt(2-25)式中：為線性水庫的瞬時出流系數；St為線性水庫t時刻的蓄量；Qt為線性水庫t時刻的出流流量合并（2-25）和（2-26）得至I(2-27)(2-28)dQdtQ在b,t1上積分上式得Qt=Q0e-t此式稱為線性水庫退水方程，稱Kr=e.作為線性水

24、庫消退系數（又稱退水系數）此式的另一個表達形式是：Qt=Qte&&=QtKr（2-29）這是（2-27）在區間t,t+N上積分的結果。令&=1,即取At為一個單位計算時段，則上式記為Qt+=Qte4（2-30）在區間t,t+N】上積分Qt得到水庫在占時段內流出的水量t:同S=Qtdt=Qoe-tdt=-%e|t；%651-一1=61飛一(2-31)(2-32)比較（2-25）和令人=1,即取M為一個單位計算時段，則上式記為S=St1-e'=二St（2-32）中的a就是線性水庫單位時段為&=1的出流系數（2-32）可知口和人的意義是完全不同的，a總是與一個

25、單位時段相聯系的，而且單位時段不同，其值也不同。.1令&=記單位時段為&的水庫消退系數為Kr。由（2-30）的退水m方程，可推導得2Qt1=QtKr=Qt2-Qt1Kr=QtKr=Qt%=QtK,(2-33)此式中的1代表一個&時段，m表示m個&'時段。顯然Qt=Qt林。對單位時段&的退水方程，有Qt=QtKr；或Qt+=QtKr(2-34)比較（2-33）最后一個式子與（2-34）式可知，線性水庫單位計算時段為&'的出流系數K；與單位計算時段為&的出流系數K,之間有如下關系：,mKr=Kr(2-35)已知對單位計算時段&

26、amp;,有0（=（1一KJ；記單位計算時段t的出流系數為豆=（1一K,），則可推得ot1)1-Krm.)(2-36)此式就是所需要得結果。2.有二個出流孔的線性水庫不同時段出流系數的轉換、,一一一一一.一-、一、-一一、,一一一一一.一'記單位計算時段為&的出流系數分別為KI和KG,單位計算時段為At的出流系數分別為KI和KG,不同計算時段的兩孔總出流系數分別為KI+KG,代人（2-36）式得到：KI+KG和(KIKG)-1-KI-KG1m(2-37)要求解KI和KG還須增加一個條件，可假定不同計算時段的兩孔出流系數分配比例相同，即成立：(2-38)KI_KIKGKG(2-3

27、9)禾I用（2-38）消去KG'得至IJ:KI1-1-KIKGKG1-KI禾I用（2-38）消去KG'得至IJ:1-1-KIKG如KGKIKG(2-40)利用（2-39）與（2-40）就可進行線性水庫不同時段的出流系數轉換。2.4 新安江三水源模型的參數調試采用河網匯流單位線的新安江三水源模型共有14個參數，根據參數的功能可分為4類與蒸散發計算有關的參數：蒸散發折算系數K,上層土壤蓄水容量WUM下層土壤蓄水容量WLM,深層蒸散發系數Co與產流計算有關的參數：流域蓄水容量Wm,流域蓄水容量曲線指數b,流域不透水面積比值IMPo與水源劃分有關的參數：流域自由水蓄水容量Sm,流域自由

28、水蓄水容量曲線指數bx,自由水蓄水庫壤中流出流系數KI與自由水蓄水庫地下徑流出流系數KG0與匯流有關的參數：壤中流消退系數CI和地下徑流消退系數CG；河網匯流單位線參數Cr和e。壤中流匯流調蓄水庫和地下徑流匯流調蓄水庫與自由水水庫是不同的水庫，但如果假定二調蓄水庫的出流系數與自由水水庫對應的徑流的出流系數相同，即假定CI=1-KI,CG=1-KG,則參數可減為12個。4類參數相對獨立，同類參數有一定相依性。1 .參數初定根據流域自然地理資料、水文氣象資料和前人以積累的經驗初定參數。蒸散發折算系數K：可采用汛期連續大雨，流域蓄水量處于蓄滿狀態的累計降雨、累計徑流、累計蒸發觀測資料，用下式估計：(

29、2-41)“P-%RK二、E流域蓄水容量Wm、上層土壤蓄水容量WUM、下層土壤蓄水容量WLM:Wm值一般根據實測的降雨徑流資料分析或經驗選取后，適當微調即可。經驗選定后，此值調試范圍不大，比較容易確定。通常在南方濕潤地區Wm可選80150mm左右，北方半濕潤地區可達170mm,江淮一帶約在110mm左右，燕山東北東部地區，選150mm左右。WUM中包含有植物截留，少林地510mm,多林地1020mm左右。WLM在6090mm范圍。流域蓄水容量曲線指數b：此值一般在0.20.5之間取值，通常流域越小，該值越小，特小流域可能小于0.2。1000平方公里以上流域一般大于0.3。不透水面積比值IMP:

30、IMP值對模型擬合精度影響很小，一般取00.05。深層蒸散發系數C:南方多林地區C=0.18左右；北方半干旱地區C=0.08左右。流域自由水蓄水容量Sm:Sm一般由經驗初定，土層淺薄的山區Sm=815mm,土厚多林地區Sm=50mm,甚至更大。流域自由水蓄水容量曲線指數bx在；根據經驗，bx常在11.5左右。自由水蓄水庫壤中流出流系數與KI自由水蓄水庫地下徑流出流系數KG:一般來說，KI+KG約為0.7左右，二者分配比與具體流域有關，只能調試決定。壤中流消退系數CI和地下徑流消退系數CG:一般CI趨近0.9;CG約為0.980.998,具體取值與具體流域有關，可進一步調試決定。河網匯流單位線參

31、數Cr和ECr在0.02左右，與具體流域的河網特性有關，可進一步調試決定。工與具體流域的河網特性有關，為計算方便，通常取計算單位時段的整倍數。2 .參數率定首先率定第一類參數，主要觀察率定期的總水量是否與實測水量相近。其次率定第二類參數，主要觀察次洪水的水量是否與實測水量相近。第三類參數的率定，主要觀察次洪水的退水段過程的水量是否與實測水量相近。第四類參數的率定，主要觀察次洪水的過程線是否擬合良好。參數率定過程中要反復對參數進行多遍調試，直至滿足精度要求。2.5 新安江三水源模型的VB程序介紹新安江三水源，兩層蒸發模型(時變單位線)VB源程序FunctionXingAnJang4()AsInt

32、egerDimFnoAsInteger,IAsInteger,Nian(2)AsInteger,AX(20)AsSingle,QoAsSingle,DimWmAsSingle,WumAsSingle,WlmAsSingle,kAsSingle,SmAsSingle,bAsSingle,BxAsSingleDimKssAsSingle,KgAsSingle,CRssAsSingle,CRgAsSingle,ImmpAsSingle,FAsSingle,UAsSingleDimPAsSingle,EAsSingle,PeAsSingle,EuAsSingle,ElAsSingle,EEAsSing

33、le,FR0AsSingleDimRAsSingle,RsAsSingle,RssAsSingle,RgAsSingle,QAsSingle,XAsSingle,XXAsSingleDimWAsSingle,WuAsSingle,WlAsSingle,SAsSingle,WmmAsSingle,SmmAsSingleDimQRss0AsSingle,QRg0AsSingle,CrAsSingle,CsAsSingle,CxAsSingle,QzAsSingle,Qa()AsSingleQo=5'Qo是基流。流域模型enshi9699parm3.txt"ForInputAs#1

34、流域模型enshi9699ensh9699.txt"ForInputAs#5流域模型enshi9699end96994.txt"ForOutputAs#6Fno=1ReadRecordFno:TurnRecordNian,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=0:Dt=3:F=2928:U=F/(3.6*Dt):QRssO=0:QR

35、g0=10.1Wmm=(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=6Fno=5:n=0:DoUntilEOF(Fno):LineInput#Fno,TextLien:n=n+1:Loop:Seek#Fno,1ReDimQa(n+1000)ForI=1TonReadRecordFno:TurnRecordNian,AX,4+1,1n1=Nian(1):m1=Int(AX(1):P=AX(2):E=AX(3):qb=AX(4)Pe=P-k*EIfPe>0ThenEu=k*E:El=0Xu

36、ManChanLiu1Pe,R,Wm,Wmm,Wum,Wlm,W,Wu,Wl,b,XElseR=0EuElP,k,E,Wlm,Wu,Wl,Eu,ElEndIfRsRssRglPe,R,S,Sm,Smm,Rs,Rss,Rg,Bx,XX,FR0,Kss,Kg,Immp'RsRssRg2Pe,R,S,Sm,Smm,Rs,Rss,Rg,Bx,XX,FR0,Kss,Kg,ImmpEE=Eu+El:W=Wu+WlHeWangZongRuLiuRs,Rss,Rg,QRss0,QRg0,Qz,CRg,CRss,Qo,UTao=1Cs=1-Cr*QzA0.4'Tao=Int(Cx/(1-Cs)/

37、3Qa(I+Tao)=Cs*Qa(I)+(1-Cs)*QzQ=Qa(I+Tao)Write#6I,P,EE,R,Rs,Rss,Rg,W,Q,qbEndIfNextIClose(6)Form1.Text1.Text=CStr(666)EndFunctionPublicSp(20)AsVariantFunctionReadRecord(Fno)AsIntegerDimTextLienAsString,StextAsString,sssAsStringDimIAsInteger,jAsInteger,LineLenAsIntegerj=1LineInput#Fno,TextLienLineLen=L

38、en(TextLien)ForI=1ToLineLenStext=Mid(TextLien,I,1)SelectCaseAsc(Stext)Case45,46,48,49,50,51,52,53,54,55,56,57sss=sss+StextIfI=LineLenThenSp(j)=sssEndIfCaseElseIfsss<>""ThenSp(j)=sss:j=j+1:sss=""EndIfEndSelectNextIReadRecord=jEndFunction'蓄滿產流函數一FunctionXuManChanLiu1(Pe,R,Wm,Wmm,Wum,Wlm,W,Wu,Wl,b,X)AsIntegerDimDrAsSingle,AAsSingleA=Wmm*(1-(1-W/Wm)AX)IfPe+A<WmmThenR=Pe-Wm+W+Wm*(1-(Pe+A)/Wmm)a(1+b)ElseR=Pe+W-WmEnd