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

上傳人：奔*** IP屬地：河北更新時間：2024-03-01 格式：doc 頁數(shù)：17 大?。?38.00KB 人氣指數(shù)：12 舉報(bào) 版權(quán)申訴

已閱讀1頁，還剩16頁未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

19、;p>　　本報(bào)告仍選擇了Kc=1.954作為率定值，因?yàn)檫x擇這個數(shù)值時1989年和1990年兩年的誤差最相近（都約為21%），這表明此時模型的穩(wěn)定性較好。其次，如果使Kc增加或減小，都會使其中一年的誤差減小，但另一年的誤差將變得特別大。綜上所述，本報(bào)告采用Kc=1.954是合理的。　　通過采用1991年和1992年的數(shù)據(jù)對率定的Kc進(jìn)行徑流計(jì)算并檢驗(yàn)，結(jié)果表明Kc=1.954對于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代表性較弱，甚至可能是錯誤的；2、徑流計(jì)算過程中所使用的水文模型屬于概念性的流域水文模型，其本身的精度就有限；3、在計(jì)算流域徑流過程中各個雨量站的權(quán)重因子的選擇可能沒有充分考慮流域的流域氣候和地質(zhì)的時空不均勻性等。

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

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

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

24、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ū)降雨量　　Dim E0(1 To 2200) As Single '流域的水面蒸發(fā)

25、Dim Ep(1 To 367) As Single '流域的蒸發(fā)能力　　Dim e(1 To 367) As Single '流域總的蒸發(fā)量　　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;　　Dim Q(1 To 2200) As Single '流域?qū)崪y徑流　　Dim R(1 To 367) As Single '流域計(jì)算徑流　　Dim RS(1 To 367) As Single, RG(1 To 367) As Single '地表及地下徑流

27、　　Dim PE(1 To 367) As Single　　Dim W(1 To 367) As Single '流域總的蓄水量　　Dim WU(1 To 367) As Single, WL(1 To 367) As Single, WD(1 To 367) As Single　　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　　Dim sumQo, sumRo, verQ As Single　　Const h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的權(quán)重<p

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

54、　　End If　　End If　　Next i　　JD(k, m) = sumQo - sumRo '求解絕對誤差　　XD(k,

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

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

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

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

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

60、;/b>　　Private Sub Command1_Click() '選擇資料檢驗(yàn)Kc　　Dim i As Integer　　Dim P(1 To 2400) As Single '流域的平均降雨量　　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　　Dim E0(1 To 2400) As Single '流域的水面蒸發(fā)　　Dim Ep(1 To 366) As Single '流域的蒸發(fā)能力<

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

78、　　End If　　End If　　e(i) = EU(i) + EL(i) + ED(i)　　W(i) = WU(i) + WL(i) + WD(i)　　&

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

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

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

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

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

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

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

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

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

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

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

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

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

92、lt;/p>　　Private Sub Command1_Click() '次洪水徑流預(yù)報(bào)計(jì)算　　Dim i As Integer　　Dim P(1 To 2400) As Single '流域的平均降雨量　　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　　Dim E0(1 To 2400) As Single '流域的水面蒸發(fā)　　Dim Ep(1 To 366) As Single '流域的蒸發(fā)能力

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

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

96、lt;/p>　　Dim WU(1 To 366) As Single, WL(1 To 366) As Single, WD(1 To 366) As Single　　Dim WMM As Single　　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　　Dim sumQ, sumR As Single　　Const h1 = 0.33, h2 = 0.14, h3 = 0.33, h4 = 0.2 '各雨量站的權(quán)重　　Const Wm = 140, Um = 20, Lm = 60, Dm = 60

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

119、gt;　　RS(i) = 0　　RG(i) = R(i)　　End If　　Next i　　For i = 1 To

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

121、　　Else　　m = 0　　End If　　n = m + n　　Next j&

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

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

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

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

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

眾賞文庫> 全部分類> 畢業(yè)設(shè)計(jì)

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

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

文檔簡介

溫馨提示

最新文檔

評論

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

文檔簡介

溫馨提示

最新文檔

評論

免費(fèi)下載