水文預(yù)報(bào)課程設(shè)計(jì)

1、<p><b>　　基本任務(wù)</b></p><p><b>　　任務(wù)一：</b></p><p>　　根據(jù)已給數(shù)據(jù)資料及參數(shù)（本報(bào)告采用87-92年的歷史數(shù)據(jù)），將流域作為整體：</p><p>　　（1）進(jìn)行日模型產(chǎn)流量計(jì)算；</p><p>　?。?）比較計(jì)算年徑流與實(shí)測(cè)年徑流；&l

2、t;/p><p>　?。?）通過誤差分析，優(yōu)選蒸發(fā)折算系數(shù)Kc；</p><p>　?。?）87~90年的歷時(shí)數(shù)據(jù)作為率定參數(shù)，91~92年的數(shù)據(jù)作為模型檢驗(yàn)。</p><p><b>　　任務(wù)二：</b></p><p>　　根據(jù)已給的設(shè)計(jì)暴雨資料和任務(wù)一率定的Kc，將流域作為整體進(jìn)行如下計(jì)算：</p>&l

3、t;p>　?。?）次洪產(chǎn)流量計(jì)算，劃分水源；</p><p>　?。?）直接徑流匯流，地下徑流匯流的計(jì)算。</p><p>　　（3）采用2004年暴雨數(shù)據(jù)進(jìn)行預(yù)報(bào)。 </p><p>　　根據(jù)已給的資料、參數(shù)及做過的習(xí)題，自己編寫程序，將流域作為整體進(jìn)行產(chǎn)流量計(jì)算；將計(jì)算年徑流與實(shí)測(cè)年徑流進(jìn)行比較。</p><p><b>

4、　　第二章 基本資料</b></p><p><b>　　1、流域概況</b></p><p>　　某流域集雨面積856km2。流域地處南方，海洋性氣候顯著，氣候溫和，雨量豐沛。暴雨成因主要是鋒面雨和臺(tái)風(fēng)雨，常受熱帶風(fēng)暴影響。降雨年際間變化大，年內(nèi)分配不均，多年平均降雨量為1800mm，實(shí)測(cè)年最大降雨量為3417mm，汛期4—9月降雨量占年降雨量的81％左

5、右：徑流系數(shù)0.5~0.7。流域內(nèi)地勢(shì)平緩，土壤主要有黃壤和砂壤，具有明顯的腐殖層，淀積層和母質(zhì)土等層次結(jié)構(gòu)，透水性好。臺(tái)地、丘陵多生長(zhǎng)松、杉、樟等高大喬木；平原則以種植農(nóng)作物和經(jīng)濟(jì)作物為主，植被良好。</p><p>　　流域上游有一水文站，控制流域面積553km2，占流域集雨面積的64.6%。該水文站以上有4個(gè)雨量站。雨量站分布較均勻，有10年逐日降水資料和時(shí)段降水資料；該水文站具有10年以上水位、流量資料；

6、流域?qū)偕絽^(qū)性小流域且受到地形、地貌等下墊面條件影響，洪水陡漲緩落，匯流時(shí)間一般2—3h，有時(shí)更短；一次洪水總歷時(shí)2~5d。</p><p><b>　　2、基本資料：</b></p><p>　　(1) 計(jì)算流域面積為553km2。</p><p>　　(2) 流域內(nèi)有四個(gè)雨量站，權(quán)重系數(shù)分別為0.33、0.14、0.33、0.20。</

7、p><p>　　(3) 資料文件名：87-92data.xls，數(shù)據(jù)格式為：</p><p>　　T(i) Q(i) E(i)P1(i) P2(i) P3(i) P4(i)</p><p><b>　　3.計(jì)算參數(shù)：</b></p><p>　　計(jì)算參數(shù)見表2-1。</p><p>　　表

8、2-1 計(jì)算參數(shù)表</p><p>　　Qg=55.3m3/s</p><p><b>　　第三章 計(jì)算公式</b></p><p><b>　　1、蒸散發(fā)計(jì)算</b></p><p>　　根據(jù)流域特點(diǎn)，蒸散發(fā)計(jì)算采用的是三層蒸散發(fā)計(jì)算模式。</p><p>　　三層蒸發(fā)模式

9、的具體計(jì)算如下：</p><p>　　1)當(dāng)WU+P>=EP, </p><p>　　EU=Ep,EL=0,ED=0;</p><p>　　2)當(dāng)WU+P<EP, WL>=C*WLM, </p><p>　　EU=WU+P,EL=(EP-EU)*WL/WLM,ED=0;</p><p>　　3)當(dāng)W

10、U+P<EP, C* (EP-EU)<=WL<C*WLM, </p><p>　　EU=WU+P,EL=C*(EP-EU),ED=0;</p><p>　　4)當(dāng)WU+P<EP，WL<C*(EP-EU)，</p><p>　　EU=WU+P,EL=WL,ED=C*(EP-EU)-EL</p><p><b&

11、gt;　　2、產(chǎn)流量計(jì)算：</b></p><p>　　根據(jù)流域特點(diǎn)，產(chǎn)流量計(jì)算系根據(jù)蓄滿產(chǎn)流理論得出的。蓄滿產(chǎn)流，即任一地點(diǎn)上，土壤含水量達(dá)田間持水量前，降雨量全部補(bǔ)充土壤含水量，不產(chǎn)流；當(dāng)土壤蓄滿后，其后續(xù)降雨量全部產(chǎn)生徑流。 </p><p>　　a=WMM(1-(1-W0/WM)^(1/(b+1))</p><p>　　PE>0，則

12、產(chǎn)流；否則不產(chǎn)流。產(chǎn)流時(shí)，</p><p>　　當(dāng)PE+a<WMM：</p><p>　　R=PE+W-WM+WM*(1-(PE+a)/WMM)^(b+1)</p><p>　　當(dāng)PE+a>=WMM：</p><p><b>　　R=PE+W-WM</b></p><p><b&

13、gt;　　3、水源劃分：</b></p><p>　　流域坡地上的降雨產(chǎn)流量因產(chǎn)流過程的條件和運(yùn)動(dòng)路徑不同，受流域的調(diào)蓄作用不同，各徑流成分在流量過程線上的反應(yīng)是不一樣的。</p><p><b>　　各水源分量如下：</b></p><p>　　1）當(dāng)PE〈=FC時(shí)：RD=0.0</p><p><b

14、>　　RG=R</b></p><p>　　2） 當(dāng)PE〉FC時(shí)：RG=FC*R/PE</p><p><b>　　RD=R-RG</b></p><p><b>　　4、匯流計(jì)算：</b></p><p>　　根據(jù)流域凈雨和流域徑流單位線，采用卷積的差分形式算出流域出口的流量過

15、程。</p><p><b>　　計(jì)算公式：</b></p><p>　　QS(I)=RS(I)*U</p><p>　　QG(I)=CG*RG(I)*U+(1-CG)*QG(I-1)</p><p>　　Q(I)=QS(I)+QG(I)</p><p><b>　　第四章 基本數(shù)據(jù)&l

16、t;/b></p><p><b>　　1、暴雨過程</b></p><p>　　1987~1992年數(shù)據(jù)見文件87-92daa.xls，2004年暴雨過程數(shù)據(jù)見表4-1。</p><p>　　表4-1 2004年暴雨過程數(shù)據(jù)表</p><p>　　用給定權(quán)重計(jì)算流域面平均雨量。</p><p&

17、gt;<b>　　2、單位線</b></p><p>　　單位線過程（m3/s）為：0，40，80，130，100，80，48，20，10，5，0</p><p>　　第五章 結(jié)果評(píng)定及誤差分析</p><p>　　1、率定Kc值并驗(yàn)證</p><p>　　Kc優(yōu)選原則：計(jì)算的2年資料的Kc應(yīng)相同并使得2年內(nèi)每年的年徑流

18、相對(duì)誤差盡可能不超過5％，編程中解釋為1989年和1990年兩年的年徑流量相對(duì)誤差的絕對(duì)值之差不超過5％，即：Abs(Abs(XD(k, 3)) - Abs(XD(k, 4)))<= 5％。通過vb編程優(yōu)選結(jié)果顯示如下圖5-1，統(tǒng)計(jì)結(jié)果如下表5-1。</p><p>　　圖5-1 程序運(yùn)行顯示結(jié)果</p><p>　　表5-1 Kc率定結(jié)果和檢驗(yàn)結(jié)果統(tǒng)計(jì)</p><

19、;p>　　本報(bào)告仍選擇了Kc=1.954作為率定值，因?yàn)檫x擇這個(gè)數(shù)值時(shí)1989年和1990年兩年的誤差最相近（都約為21%），這表明此時(shí)模型的穩(wěn)定性較好。其次，如果使Kc增加或減小，都會(huì)使其中一年的誤差減小，但另一年的誤差將變得特別大。綜上所述，本報(bào)告采用Kc=1.954是合理的。</p><p>　　通過采用1991年和1992年的數(shù)據(jù)對(duì)率定的Kc進(jìn)行徑流計(jì)算并檢驗(yàn)，結(jié)果表明Kc=1.954對(duì)于1991年

20、和1992年的徑流計(jì)算誤差較大，且其相應(yīng)的確定系系數(shù)較小不能滿足項(xiàng)目預(yù)報(bào)的要求。另一方面，通過檢驗(yàn)發(fā)現(xiàn)檢驗(yàn)?zāi)攴莸恼`差和確定性系數(shù)都不如率定年份，究其原因分析如下：1、用來率定參數(shù)的資料序列太短，資料的代表性不強(qiáng)，則率定出的Kc代表性較弱，甚至可能是錯(cuò)誤的；2、徑流計(jì)算過程中所使用的水文模型屬于概念性的流域水文模型，其本身的精度就有限；3、在計(jì)算流域徑流過程中各個(gè)雨量站的權(quán)重因子的選擇可能沒有充分考慮流域的流域氣候和地質(zhì)的時(shí)空不均勻性等。

21、由于多種因素的作用使得水文模型整體上計(jì)算誤差都較大。</p><p><b>　　2、次洪徑流計(jì)算</b></p><p>　　在任務(wù)一率定Kc=1.954的基礎(chǔ)上，通過vb編程計(jì)算的到次洪徑流過程如下表5-2，次洪徑流過程如圖5-2。通過計(jì)算模擬的次洪徑流過程發(fā)現(xiàn)，該次洪水過程為單峰，洪峰流量為2745 m3/s出現(xiàn)時(shí)間為2004-9-25 6:00，該次洪水直接徑

22、流的持續(xù)時(shí)間為2004-9-23 15:00—— 2004-9-26 12:00。另一方面，通過對(duì)比此次洪水的直接徑流和地下徑流不難發(fā)現(xiàn)，該次洪水的直接徑流所占比例較大，在洪峰形成過程中起決定作用。</p><p>　　表5-2 次洪徑流計(jì)算結(jié)果</p><p>　　圖5-2 洪水流量過程線圖</p><p>　　第六章 計(jì)算程序及說明</p><

23、;p>　　報(bào)告vb編程計(jì)算程序：</p><p>　　1、確定參數(shù)kc ；</p><p>　　Private Sub Command1_Click() '確定KC</p><p>　　Dim i As Integer</p><p>　　Dim P(1 To 2200) As Single '流域的平均降雨量</

24、p><p>　　Dim P1(1 To 2200) As Single, P2(1 To 2200) As Single, P3(1 To 2200) As Single, P4(1 To 2200) As Single '流域各地區(qū)降雨量</p><p>　　Dim E0(1 To 2200) As Single '流域的水面蒸發(fā)</p><p>　　

25、Dim Ep(1 To 367) As Single '流域的蒸發(fā)能力</p><p>　　Dim e(1 To 367) As Single '流域總的蒸發(fā)量</p><p>　　Dim EU(1 To 367) As Single, EL(1 To 366) As Single, ED(1 To 367) As Single '三層蒸發(fā)模式參數(shù)</p&g

26、t;<p>　　Dim Q(1 To 2200) As Single '流域?qū)崪y(cè)徑流</p><p>　　Dim R(1 To 367) As Single '流域計(jì)算徑流</p><p>　　Dim RS(1 To 367) As Single, RG(1 To 367) As Single '地表及地下徑流</p><p>

27、　　Dim PE(1 To 367) As Single</p><p>　　Dim W(1 To 367) As Single '流域總的蓄水量</p><p>　　Dim WU(1 To 367) As Single, WL(1 To 367) As Single, WD(1 To 367) As Single</p><p>　　Dim WMM, su

28、mQ(1 To 4), sumR(1 To 4), a(1 To 367), kaishi, tianshu, ymdhm(1 To 2200), XD(1 To 4000, 1 To 4), XDH(4 To 4000), JD(1 To 4000, 1 To 4), qo(1 To 2200), E0o(1 To 2200), P1o(1 To 2200), P2o(1 To 2200), P3o(1 To 2200), P4o(1

29、 To 2200), DC(1 To 4), qidian, minXDH, mink, s, n, f, d As Single</p><p>　　Dim sumQo, sumRo, verQ As Single</p><p>　　Const h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的權(quán)重</p><p

30、>　　Const Wm = 140, Um = 20, Lm = 60, Dm = 60, B = 0.3, C = 0.16, IM = 0.002, Fc = 24 '模型參數(shù)</p><p>　　Open "f:課程設(shè)計(jì)資料數(shù)據(jù)1.txt" For Input As #1 '讀取數(shù)據(jù)庫(kù)數(shù)據(jù)</p><p>　　For i = 1 To 21

31、92</p><p>　　Input #1, ymdhm(i), qo(i), E0o(i), P1o(i), P2o(i), P3o(i), P4o(i)</p><p><b>　　Next i</b></p><p><b>　　Close #1</b></p><p>　　minXDH =

32、4</p><p>　　For k = 50 To 3000 'Kc=k/1000</p><p>　　XD(k, 1) = 0</p><p>　　XD(k, 2) = 0</p><p>　　XD(k, 3) = 0</p><p>　　XD(k, 4) = 0</p><p>&l

33、t;b>　　n = 0</b></p><p><b>　　s = 0</b></p><p>　　For m = 3 To 4</p><p>　　If m = 3 Then</p><p>　　kaishi = 732</p><p>　　tianshu = 365</

34、p><p>　　ElseIf m = 4 Then</p><p>　　kaishi = 1097</p><p>　　tianshu = 365</p><p><b>　　End If</b></p><p><b>　　sumQo = 0</b></p>&l

35、t;p><b>　　sumRo = 0</b></p><p>　　For j = 1 To tianshu '使帶入循環(huán)的各年降雨和蒸發(fā)量對(duì)應(yīng)</p><p>　　Q(j) = qo(kaishi + j - 1)</p><p>　　E0(j) = E0o(kaishi + j - 1)</p><p>

36、;　　P1(j) = P1o(kaishi + j - 1)</p><p>　　P2(j) = P2o(kaishi + j - 1)</p><p>　　P3(j) = P3o(kaishi + j - 1)</p><p>　　P4(j) = P4o(kaishi + j - 1)</p><p>　　P(j) = h1 * P1(j)

37、 + h2 * P2(j) + h3 * P3(j) + h4 * P4(j) '流域平均降雨量計(jì)算</p><p>　　sumQo = sumQo + Q(j) * 24 * 3600 / 553000</p><p><b>　　Next j</b></p><p>　　verQ = sumQo / tianshu</p&g

38、t;<p><b>　　'流域三層蒸發(fā)計(jì)算</b></p><p>　　W(1) = 110</p><p>　　WU(1) = 10</p><p>　　WL(1) = 40</p><p>　　WD(1) = 60</p><p>　　W(1) = WU(1) + WL(

39、1) + WD(1)</p><p>　　For i = 1 To tianshu</p><p>　　Ep(i) = E0(i) * k / 1000</p><p>　　If WU(i) + P(i) >= Ep(i) Then</p><p>　　EU(i) = Ep(i)</p><p><b>

40、;　　EL(i) = 0</b></p><p><b>　　ED(i) = 0</b></p><p><b>　　End If</b></p><p>　　If WU(i) + P(i) < Ep(i) Then</p><p>　　If WL(i) >= C * Lm

41、Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * WL(i) / Lm</p><p><b>　　ED(i) = 0</b></p><p>　　ElseIf WL(i) < C * Lm And WL(i) >=

42、 C * (Ep(i) - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * C</p><p><b>　　ED(i) = 0</b></p><p>　　ElseIf WL(i) < C * (Ep(i)

43、 - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = WL(i)</p><p>　　ED(i) = (Ep(i) - EU(i)) * C - EL(i)</p><p><b>　　End If</b></p><p>

44、<b>　　End If</b></p><p>　　e(i) = EU(i) + EL(i) + ED(i)</p><p>　　W(i) = WU(i) + WL(i) + WD(i)</p><p><b>　　'流域產(chǎn)流計(jì)算</b></p><p>　　WMM = Wm * (1 +

45、 B)</p><p>　　a(i) = WMM * (1 - (1 - (W(i) / Wm)) ^ (1 / (1 + B)))</p><p>　　PE(i) = P(i) - e(i)</p><p><b>　　'產(chǎn)流時(shí)</b></p><p>　　If PE(i) > 0 Then</p&

46、gt;<p>　　If PE(i) + a(i) < WMM Then</p><p>　　R(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)</p><p>　　W(i + 1) = W(i) + PE(i) - R(i)</p><p><b>　

47、　Else</b></p><p>　　R(i) = PE(i) + W(i) - Wm</p><p>　　W(i + 1) = Wm</p><p><b>　　End If</b></p><p><b>　　End If</b></p><p><b

48、>　　'不產(chǎn)流時(shí)</b></p><p>　　If PE(i) <= 0 Then</p><p><b>　　R(i) = 0</b></p><p>　　W(i + 1) = W(i) + PE(i)</p><p><b>　　End If</b></p&

49、gt;<p>　　sumRo = sumRo + R(i)</p><p>　　If WU(i) + P(i) - EU(i) - R(i) <= Um Then '求下一時(shí)刻地下各層含水量</p><p>　　WU(i + 1) = WU(i) + P(i) - EU(i) - R(i)</p><p>　　WL(i + 1) = W

50、L(i) - EL(i)</p><p>　　WD(i + 1) = WD(i) - ED(i)</p><p>　　ElseIf WU(i) + P(i) - EU(i) - R(i) > Um Then</p><p>　　WU(i + 1) = Um</p><p>　　If WL(i) - EL(i) + (WU(i) + P(

51、i) - EU(i) - R(i) - Um) <= Lm Then</p><p>　　WL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)</p><p>　　WD(i + 1) = WD(i) - ED(i)</p><p>　　ElseIf WL(i) - EL(i) + (WU(

52、i) + P(i) - EU(i) - R(i) - Um) > Lm Then</p><p>　　WL(i + 1) = Lm</p><p>　　If WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then</p><p>　　WD(i +

53、 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm</p><p><b>　　Else</b></p><p>　　WD(i + 1) = Dm</p><p><b>　　End If</b></p>

54、<p><b>　　End If</b></p><p><b>　　End If</b></p><p><b>　　Next i</b></p><p>　　JD(k, m) = sumQo - sumRo '求解絕對(duì)誤差</p><p>　　XD(k,

55、m) = JD(k, m) / sumQo '求解相對(duì)誤差</p><p>　　sumQ(m) = sumQo</p><p>　　sumR(m) = sumRo</p><p><b>　　Next m</b></p><p>　　XDH(k) = Abs(Abs(XD(k, 3)) - Abs(XD(k, 4

56、))) '對(duì)兩年相對(duì)誤差絕對(duì)值之差取絕對(duì)值</p><p>　　If minXDH > XDH(k) Then '尋找相對(duì)誤差之差最小的點(diǎn)</p><p>　　minXDH = XDH(k) '輸出最優(yōu)KC下各值</p><p><b>　　mink = k</b></p><p&

57、gt;　　Text1.Text = Round(sumQ(3), 5)</p><p>　　Text2.Text = Round(sumR(3), 5)</p><p>　　Text3.Text = Round(JD(k, 3), 5)</p><p>　　Text4.Text = Round(XD(k, 3), 5)</p><p>　　T

58、ext5.Text = Round(sumQ(4), 5)</p><p>　　Text6.Text = Round(sumR(4), 5)</p><p>　　Text7.Text = Round(JD(k, 4), 5)</p><p>　　Text8.Text = Round(XD(k, 4), 5)</p><p><b>

59、　　End If</b></p><p><b>　　Next k</b></p><p>　　Text9.Text = mink / 1000 '輸出最優(yōu)的KC值</p><p><b>　　End Sub</b></p><p><b>　　2、檢驗(yàn)參數(shù)Kc；<

60、;/b></p><p>　　Private Sub Command1_Click() '選擇資料檢驗(yàn)Kc</p><p>　　Dim i As Integer</p><p>　　Dim P(1 To 2400) As Single '流域的平均降雨量</p><p>　　Dim P1(1 To 2400) As Si

61、ngle, P2(1 To 2400) As Single, P3(1 To 2400) As Single, P4(1 To 2400) As Single</p><p>　　Dim E0(1 To 2400) As Single '流域的水面蒸發(fā)</p><p>　　Dim Ep(1 To 366) As Single '流域的蒸發(fā)能力</p><

62、p>　　Dim e(1 To 366) As Single '流域總的蒸發(fā)量</p><p>　　Dim EU(1 To 366) As Single, EL(1 To 366) As Single, ED(1 To 366) As Single</p><p>　　Dim Q(1 To 2400) As Single '流域?qū)崪y(cè)徑流</p><p

63、>　　Dim R(1 To 366) As Single '流域計(jì)算徑流</p><p>　　Dim RS(1 To 366) As Single, RG(1 To 366) As Single</p><p>　　Dim PE(1 To 366) As Single</p><p>　　Dim W(1 To 367) As Single '流

64、域總的蓄水量</p><p>　　Dim WU(1 To 367) As Single, WL(1 To 367) As Single, WD(1 To 367) As Single</p><p>　　Dim WMM As Single</p><p>　　Dim a(1 To 366), kaishi, tianshu, ymdhm(1 To 2200), qi

65、dian, f, d6 As Single</p><p>　　Dim sumQ, sumR, verQ As Single</p><p>　　Const h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的權(quán)重</p><p>　　Const Wm = 140, Um = 20, Lm = 60, Dm = 6

66、0, B = 0.3, C = 0.16, Fc = 24 '模型參數(shù)</p><p><b>　　sumQ = 0</b></p><p><b>　　sumR = 0</b></p><p>　　Kc = 1.954</p><p>　　Open "f:\課程設(shè)計(jì)資料數(shù)據(jù)1.

67、txt" For Input As #1</p><p>　　For i = 1 To 2192</p><p>　　Input #1, ymdhm(i), Q(i), E0(i), P1(i), P2(i), P3(i), P4(i)</p><p><b>　　Next i</b></p><p><

68、b>　　Close #1</b></p><p>　　If Option1.Value = True Then '判別所選擇的資料</p><p>　　kaishi = 1462</p><p>　　tianshu = 365</p><p>　　ElseIf Option2.Value = True Then&l

69、t;/p><p>　　kaishi = 1827</p><p>　　tianshu = 366</p><p><b>　　End If</b></p><p>　　For j = 1 To tianshu</p><p>　　ymdhm(j) = ymdhm(kaishi + j - 1)<

70、/p><p>　　Q(j) = Q(kaishi + j - 1)</p><p>　　E0(j) = E0(kaishi + j - 1)</p><p>　　P1(j) = P1(kaishi + j - 1)</p><p>　　P2(j) = P2(kaishi + j - 1)</p><p>　　P3(j) =

71、 P3(kaishi + j - 1)</p><p>　　P4(j) = P4(kaishi + j - 1)</p><p>　　P(j) = h1 * P1(j) + h2 * P2(j) + h3 * P3(j) + h4 * P4(j) '流域平均降雨量計(jì)算</p><p>　　sumQ = sumQ + Q(j) * 24 * 3600 / 5

72、53000</p><p><b>　　Next j</b></p><p>　　verQ = sumQ / tianshu '流域三層蒸發(fā)計(jì)算</p><p>　　W(1) = 110</p><p>　　WU(1) = 10</p><p>　　WL(1) = 40</p>

73、<p>　　WD(1) = 60</p><p>　　W(1) = WU(1) + WL(1) + WD(1)</p><p>　　For i = 1 To tianshu</p><p>　　Ep(i) = E0(i) * Kc</p><p>　　If WU(i) + P(i) >= Ep(i) Then</p&

74、gt;<p>　　EU(i) = Ep(i)</p><p><b>　　EL(i) = 0</b></p><p><b>　　ED(i) = 0</b></p><p><b>　　End If</b></p><p>　　If WU(i) + P(i) &l

75、t; Ep(i) Then</p><p>　　If WL(i) >= C * Lm Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * WL(i) / Lm</p><p><b>　　ED(i) = 0</b></

76、p><p>　　ElseIf WL(i) < C * Lm And WL(i) >= C * (Ep(i) - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * C</p><p><b>　　ED(i) = 0<

77、;/b></p><p>　　ElseIf WL(i) < C * (Ep(i) - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = WL(i)</p><p>　　ED(i) = (Ep(i) - EU(i)) * C - EL(i)</p>

78、<p><b>　　End If</b></p><p><b>　　End If</b></p><p>　　e(i) = EU(i) + EL(i) + ED(i)</p><p>　　W(i) = WU(i) + WL(i) + WD(i)</p><p><b>　　&

79、#39;流域產(chǎn)流計(jì)算</b></p><p>　　WMM = Wm * (1 + B)</p><p>　　a(i) = WMM * (1 - (1 - (W(i) / Wm)) ^ (1 / (1 + B)))</p><p>　　PE(i) = P(i) - e(i)</p><p><b>　　'產(chǎn)流時(shí)&l

80、t;/b></p><p>　　If PE(i) > 0 Then</p><p>　　If PE(i) + a(i) < WMM Then</p><p>　　R(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)</p><p>　　W(i

81、 + 1) = W(i) + PE(i) - R(i)</p><p><b>　　Else</b></p><p>　　R(i) = PE(i) + W(i) - Wm</p><p>　　W(i + 1) = Wm</p><p><b>　　End If</b></p><

82、p><b>　　End If</b></p><p><b>　　'不產(chǎn)流時(shí)</b></p><p>　　If PE(i) <= 0 Then</p><p><b>　　R(i) = 0</b></p><p>　　W(i + 1) = W(i) + PE

83、(i)</p><p><b>　　End If</b></p><p>　　sumR = sumR + R(i)</p><p>　　If WU(i) + P(i) - EU(i) - R(i) <= Um Then</p><p>　　WU(i + 1) = WU(i) + P(i) - EU(i) - R(i

84、)</p><p>　　WL(i + 1) = WL(i) - EL(i)</p><p>　　WD(i + 1) = WD(i) - ED(i)</p><p>　　ElseIf WU(i) + P(i) - EU(i) - R(i) > Um Then</p><p>　　WU(i + 1) = Um</p><

85、p>　　If WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm Then</p><p>　　WL(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)</p><p>　　WD(i + 1) = WD(i) - ED(i)</p>&

86、lt;p>　　ElseIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) > Lm Then</p><p>　　WL(i + 1) = Lm</p><p>　　If WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <=

87、 Dm Then</p><p>　　WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm</p><p><b>　　Else</b></p><p>　　WD(i + 1) = Dm</p><p>&l

88、t;b>　　End If</b></p><p><b>　　End If</b></p><p><b>　　End If</b></p><p><b>　　Next i</b></p><p>　　JD = sumQ - sumR</p>

89、<p>　　XD = JD / sumQ</p><p>　　If Option1.Value = True Then</p><p>　　Text1.Text = Round(sumQ, 5)</p><p>　　Text2.Text = Round(sumR, 5)</p><p>　　Text3.Text = Round(JD

90、, 5)</p><p>　　Text4.Text = Round(XD, 5)</p><p>　　ElseIf Option2.Value = True Then</p><p>　　Text5.Text = Round(sumQ, 5)</p><p>　　Text6.Text = Round(sumR, 5)</p>&

91、lt;p>　　Text7.Text = Round(JD, 5)</p><p>　　Text8.Text = Round(XD, 5)</p><p><b>　　End If</b></p><p><b>　　End Sub</b></p><p>　　3、2004年洪水徑流預(yù)報(bào)計(jì)算；&

92、lt;/p><p>　　Private Sub Command1_Click() '次洪水徑流預(yù)報(bào)計(jì)算</p><p>　　Dim i As Integer</p><p>　　Dim P(1 To 2400) As Single '流域的平均降雨量</p><p>　　Dim P1(1 To 2400) As Single, P

93、2(1 To 2400) As Single, P3(1 To 2400) As Single, P4(1 To 2400) As Single</p><p>　　Dim E0(1 To 2400) As Single '流域的水面蒸發(fā)</p><p>　　Dim Ep(1 To 366) As Single '流域的蒸發(fā)能力</p><p>　　

94、Dim e(1 To 366) As Single '流域總的蒸發(fā)量</p><p>　　Dim EU(1 To 366) As Single, EL(1 To 366) As Single, ED(1 To 366) As Single</p><p>　　Dim Q(1 To 2400) As Single '流域?qū)崪y(cè)徑流</p><p>　　D

95、im R(1 To 366) As Single '流域計(jì)算徑流</p><p>　　Dim RS(1 To 366) As Single, RG(1 To 366) As Single</p><p>　　Dim PE(1 To 366) As Single</p><p>　　Dim W(1 To 366) As Single '流域總的蓄水量&

96、lt;/p><p>　　Dim WU(1 To 366) As Single, WL(1 To 366) As Single, WD(1 To 366) As Single</p><p>　　Dim WMM As Single</p><p>　　Dim a(1 To 366), QS(0 To 38), h(1 To 38), QG(0 To 38), ymdhm(

97、1 To 2400), m, n As Single</p><p>　　Dim sumQ, sumR As Single</p><p>　　Const h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的權(quán)重</p><p>　　Const Wm = 140, Um = 20, Lm = 60, Dm = 60

98、, B = 0.3, C = 0.16, Fc = 11, Cg = 0.978 '模型參數(shù)</p><p>　　QG(0) = 55.3</p><p><b>　　QS(0) = 0</b></p><p>　　Kc = 1.954</p><p><b>　　Dim Exl</b>&l

99、t;/p><p>　　Set Exl = CreateObject("Excel.Application")</p><p>　　Set myworkbook = Exl.Workbooks.Open("f:\暴雨過程") '打開excel文件等待輸入數(shù)據(jù)</p><p>　　Exl.Visible = True</

100、p><p>　　Set xlsheet = myworkbook.Worksheets("次洪流量結(jié)果")</p><p>　　Open "f:\課程設(shè)計(jì)資料數(shù)據(jù)2.txt" For Input As #2 '讀入次洪徑流數(shù)據(jù)</p><p>　　For i = 1 To 28</p><p>　

101、　Input #2, E0(i), P1(i), P2(i), P3(i), P4(i), h(i)</p><p>　　P(i) = h1 * P1(i) + h2 * P2(i) + h3 * P3(i) + h4 * P4(i) '流域平均降雨量計(jì)算</p><p><b>　　Next i</b></p><p><b&

102、gt;　　Close #2</b></p><p><b>　　'流域三層蒸發(fā)計(jì)算</b></p><p>　　W(1) = 120</p><p>　　WU(1) = 20</p><p>　　WL(1) = 60</p><p>　　WD(1) = 60</p>

103、<p>　　W(1) = WU(1) + WL(1) + WD(1)</p><p>　　For i = 1 To 28</p><p>　　Ep(i) = E0(i) * Kc</p><p>　　If WU(i) + P(i) >= Ep(i) Then</p><p>　　EU(i) = Ep(i)</p>

104、;<p><b>　　EL(i) = 0</b></p><p><b>　　ED(i) = 0</b></p><p><b>　　End If</b></p><p>　　If WU(i) + P(i) < Ep(i) Then</p><p>　　If

105、 WL(i) >= C * Lm Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * WL(i) / Lm</p><p><b>　　ED(i) = 0</b></p><p>　　ElseIf WL(i) < C

106、* Lm And WL(i) >= C * (Ep(i) - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = (Ep(i) - EU(i)) * C</p><p><b>　　ED(i) = 0</b></p><p>　　ElseIf W

107、L(i) < C * (Ep(i) - EU(i)) Then</p><p>　　EU(i) = WU(i) + P(i)</p><p>　　EL(i) = WL(i)</p><p>　　ED(i) = (Ep(i) - EU(i)) * C - EL(i)</p><p><b>　　End If</b>&

108、lt;/p><p><b>　　End If</b></p><p>　　e(i) = EU(i) + EL(i) + ED(i)</p><p>　　W(i) = WU(i) + WL(i) + WD(i)</p><p><b>　　'流域產(chǎn)流計(jì)算</b></p><p&

109、gt;　　WMM = Wm * (1 + B)</p><p>　　a(i) = WMM * (1 - (1 - (W(i) / Wm)) ^ (1 / (1 + B)))</p><p>　　PE(i) = P(i) - e(i)</p><p><b>　　'產(chǎn)流時(shí)</b></p><p>　　If PE(i

110、) > 0 Then</p><p>　　If PE(i) + a(i) < WMM Then</p><p>　　R(i) = PE(i) + W(i) - Wm + Wm * (1 - (PE(i) + a(i)) / WMM) ^ (B + 1)</p><p>　　W(i + 1) = W(i) + PE(i) - R(i)</p>

111、<p><b>　　Else</b></p><p>　　R(i) = PE(i) + W(i) - Wm</p><p>　　W(i + 1) = Wm</p><p><b>　　End If</b></p><p><b>　　End If</b></p

112、><p>　　If WU(i) + P(i) - EU(i) - R(i) <= Um Then</p><p>　　WU(i + 1) = WU(i) + P(i) - EU(i) - R(i)</p><p>　　WL(i + 1) = WL(i) - EL(i)</p><p>　　WD(i + 1) = WD(i) - ED(i)&

113、lt;/p><p>　　ElseIf WU(i) + P(i) - EU(i) - R(i) > Um Then</p><p>　　WU(i + 1) = Um</p><p>　　If WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) <= Lm Then</p><p>　　W

114、L(i + 1) = WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um)</p><p>　　WD(i + 1) = WD(i) - ED(i)</p><p>　　ElseIf WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) > Lm Then</p><p&

115、gt;　　WL(i + 1) = Lm</p><p>　　If WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) - R(i) - Um) - Lm <= Dm Then</p><p>　　WD(i + 1) = WD(i) - ED(i) + WL(i) - EL(i) + (WU(i) + P(i) - EU(i) -

116、 R(i) - Um) - Lm</p><p><b>　　Else</b></p><p>　　WD(i + 1) = Dm</p><p><b>　　End If</b></p><p><b>　　End If</b></p><p><

117、b>　　End If</b></p><p><b>　　'不產(chǎn)流時(shí)</b></p><p>　　If PE(i) <= 0 Then</p><p><b>　　R(i) = 0</b></p><p>　　W(i + 1) = W(i) + PE(i)</p

118、><p><b>　　End If</b></p><p>　　If PE(i) > Fc Then</p><p>　　RG(i) = Fc * R(i) / PE(i)</p><p>　　RS(i) = R(i) - RG(i)</p><p><b>　　Else</b&

119、gt;</p><p><b>　　RS(i) = 0</b></p><p>　　RG(i) = R(i)</p><p><b>　　End If</b></p><p><b>　　Next i</b></p><p>　　For i = 1 To

120、 38 '求解地表徑流</p><p><b>　　n = 0</b></p><p>　　For j = 1 To 38</p><p>　　If 1 + i - j > 0 Then</p><p>　　m = h(j) * RS(1 + i - j) / 10</p><p>

121、<b>　　Else</b></p><p><b>　　m = 0</b></p><p><b>　　End If</b></p><p><b>　　n = m + n</b></p><p><b>　　Next j</b>&

122、lt;/p><p><b>　　QS(i) = n</b></p><p>　　U = 553 / (3.6 * 3)</p><p>　　QG(i) = Cg * QG(i - 1) + (1 - Cg) * RG(i) * U '求解地下徑流</p><p>　　Q(i) = QS(i) + QG(i) '

123、;求解總徑流量</p><p>　　xlsheet.cells(i + 1, 8) = QS(i) '在打開的excel中輸出地表、地下和總徑流過程</p><p>　　xlsheet.cells(i + 1, 9) = QG(i)</p><p>　　xlsheet.cells(i + 1, 10) = Q(i)</p><p>&

124、lt;b>　　Next i</b></p><p><b>　　End Sub</b></p><p><b>　　第七章 總結(jié)和心得</b></p><p>　　水文預(yù)報(bào)在防洪、抗旱、水資源開發(fā)利用、國(guó)民經(jīng)濟(jì)建設(shè)和國(guó)防等領(lǐng)域都有廣泛的應(yīng)用，經(jīng)濟(jì)效益巨大，應(yīng)用單位眾多。水文預(yù)報(bào)是與國(guó)家社會(huì)經(jīng)濟(jì)發(fā)展與安全生

125、活生產(chǎn)等密切相關(guān)的基礎(chǔ)業(yè)務(wù)，它看似單一但卻與我們的生產(chǎn)生活息息相關(guān)，十分重要。</p><p>　　通過本次課程設(shè)計(jì)，有效的溫習(xí)了產(chǎn)、匯流和蓄滿產(chǎn)流模型（三層模式）相關(guān)的知識(shí)，經(jīng)過相應(yīng)程序的編寫和計(jì)算使自己對(duì)產(chǎn)、匯流的過程和蓄滿產(chǎn)流模型的機(jī)理更加的熟練。通過對(duì)程序?qū)c的率定和檢驗(yàn)，結(jié)果顯示各年的計(jì)算徑流和實(shí)際徑流之間的相對(duì)誤差都較大。這與率定參數(shù)時(shí)資料序列較短有關(guān)，另一方面也顯示出了蓄滿產(chǎn)流模型本身的不確定，以

126、及穩(wěn)定性和適應(yīng)性差的缺點(diǎn)。當(dāng)然流域氣候和地質(zhì)的時(shí)空不均勻性，由于流域本身的地形、土壤、植被、氣候等條件各不相同，且同一流域不同地區(qū)、不同時(shí)間各要素也并不相同，都處于一個(gè)時(shí)變的狀態(tài)，這是導(dǎo)致水文模型的計(jì)算結(jié)果誤差較大的本質(zhì)原因。</p><p>　　雖然當(dāng)前的水文模型仍然存在很多的弊端和不足，但隨著計(jì)算機(jī)科學(xué)和信息技術(shù)的發(fā)展和在水文與水資源中的應(yīng)用不斷地推廣，使得流域水文模型的研究在近年取得了迅速的發(fā)展，尤其是3S