謝銀昌

(云南省電力設(shè)計(jì)院，云南 昆明 650011)

Excel函數(shù)在小流域洪水計(jì)算中的應(yīng)用

謝銀昌

(云南省電力設(shè)計(jì)院，云南 昆明 650011)

在電力工程設(shè)計(jì)的水文計(jì)算工作中，通常需要計(jì)算設(shè)計(jì)洪水，常規(guī)采用手工方法計(jì)算，計(jì)算繁冗，工作量很大。應(yīng)用Excel函數(shù)功能，計(jì)算小流域洪峰流量、繪制洪水過(guò)程線，只要經(jīng)過(guò)一些簡(jiǎn)單設(shè)置，即可完成。此方法簡(jiǎn)單、實(shí)用。

Excel函數(shù)；林平一法；小流域洪水。

在電力工程設(shè)計(jì)的水文計(jì)算工作中，經(jīng)常推求百年一遇洪水或各種保證率的流量及洪量，經(jīng)全國(guó)幾大設(shè)計(jì)院多年設(shè)計(jì)生產(chǎn)實(shí)踐表明，小流域洪水計(jì)算方法—林平一法，對(duì)電力工程洪水計(jì)算實(shí)用性較強(qiáng)，但林平一法的實(shí)現(xiàn)需進(jìn)行手工重復(fù)試算，工作量較大，Excel的“INDEX”和“MATCH”函數(shù)基本上可以解決林平一法小流域洪水計(jì)算中手工重復(fù)試算問(wèn)題，既簡(jiǎn)單又經(jīng)濟(jì)實(shí)用。實(shí)現(xiàn)雖有的專業(yè)軟件，但價(jià)格昂貴，且使用時(shí)需要專門(mén)的培訓(xùn)，經(jīng)濟(jì)和實(shí)用性上都不盡如人意。

1 計(jì)算原理

根據(jù)《電力工程水文勘測(cè)計(jì)算手冊(cè)》，林平一法適用于匯水面積在150km2以下的小流域。

1.1 全面匯流的洪峰流量公式

式中：Qm為洪峰流量(m3/s)；F為流域面積(km2)；n 為暴雨衰減指數(shù)；S為雨力(mm/min)；τ 為全面匯流時(shí)間(min)；C 為徑流系數(shù)；K 為全面匯流時(shí)間參數(shù)；μ 為穩(wěn)定下滲率(mm/min)；L為河長(zhǎng)(km)；Jc為主河道的坡度(‰)；Nc為河道匯流糙率。

1.2 部分匯流的洪峰流量公式

按全面匯流公式計(jì)算出洪峰流量和匯流時(shí)間后，若τ > tc，則屬于部分匯流情況，按下式改正流量和匯流時(shí)間：

1.3 坡地全面匯流的洪峰流量公式

式中：Q為洪峰流量(m3/s)；τ0為全面匯流時(shí)間(min)；K0為全面匯流時(shí)間參數(shù)；L0為坡面長(zhǎng)度(km)；J0為坡面坡度(‰)；N0為坡地匯流糙率。

2 Excel應(yīng)用函數(shù)

“林平一法”小流域洪水計(jì)算中手工重復(fù)試算涉及到一些特殊函數(shù)， Excel電子表格軟件中的INDEX和MATCH函數(shù)可以很好地解決上述水文計(jì)算中的問(wèn)題。

2.1 INDEX函數(shù)

定義：返回表或區(qū)域中的值或?qū)χ档囊?。INDEX函數(shù)有兩種形式：數(shù)組形式和引用形式。

語(yǔ)法：INDEX(array，row_num，column_num)，數(shù)組形式。Array是一個(gè)單元格區(qū)域或數(shù)組常量；Row_num用于選擇要從中返回值的數(shù)組中的行。如果省略 row_num，則需要使用column_num；Column_num用于選擇要從中返回值的數(shù)組中的列。如果省略 column_num，則需要使用 row_num。

INDEX(reference，row_num，column_num,area_num)，引用形式。Reference是對(duì)一個(gè)或多個(gè)單元格區(qū)域的引用；Row_num是要從中返回引用的引用中的行編號(hào)；Column_num是要從中返回引用的引用中的列編號(hào)；Area_num用于選擇要從中返回 row_num 和column_num 的交叉點(diǎn)的引用區(qū)域。選擇或輸入的第一個(gè)區(qū)域的編號(hào)是1，第二個(gè)區(qū)域的編號(hào)是2，依此類推。如果省略 area_num，則INDEX 將使用區(qū)域1。

2.2 MATCH函數(shù)

定義：返回在指定方式下與指定數(shù)值匹配的數(shù)組中元素的相應(yīng)位置。如果需要找出匹配元素的位置而不是匹配元素本身，則應(yīng)該使用MATCH 函數(shù)而不是 LOOKUP 函數(shù)。

語(yǔ)法：MATCH(lookup_value,lookup_array,match_type)。Lookup_value為需要在數(shù)據(jù)表中查找的數(shù)值；Lookup_array可能包含所要查找的數(shù)值的連續(xù)單元格區(qū)域。Lookup_array 應(yīng)為數(shù)組或數(shù)組引用；Match_type為數(shù)字 —1、0或1。Match—type指明 Microsoft Excel 如何在lookup_array 中查找 lookup_value。

3 “林平一法小流域洪水計(jì)算”Excel計(jì)算程序的編制

根據(jù)“林平一法小流域洪水計(jì)算”公式1～公式12編制如圖1所示(F6：G28)的公式計(jì)算表格。其中F6、F7、F9、F10、F11、F12、F14、F15；G6、G7、G9、G10、G11、G12、G14、G15為林平一法洪水計(jì)算相應(yīng)參數(shù)，具體為流域面積F(km2)、河長(zhǎng)L(km)、坡度Jc(‰)、糙率N、下滲率u(mm/min)、日降水量H日(mm)和暴雨衰減指數(shù)n1、n2，為計(jì)算程序中的可變單元格，其單元格中的數(shù)值可變。而其余單元格均為相應(yīng)洪水計(jì)算公式利用Excel內(nèi)部函數(shù)計(jì)算得出，均已被鎖定。

圖1 “林平一法”小流域洪水計(jì)算程序

“林平一法”洪水計(jì)算公式4和12中，徑流系數(shù)C需手工多次(從0.001～0.999中，每次以0.001為步長(zhǎng)取值)重復(fù)試算后才能計(jì)算出結(jié)果，若洪水計(jì)算中改變了任意一個(gè)計(jì)算參數(shù)，又得重復(fù)試算，直到得到正確的計(jì)算結(jié)果。手工試算周而復(fù)始，工作量大、且繁重。而Excel電子表格軟件中的INDEX和MATCH函數(shù)可以很好地解決上述問(wèn)題。

在圖1中A1單元格公式欄輸入公式：(1—B1)*B1^($F$14/(4—$F$14))，拖動(dòng)單元格右下角的填充柄，將單元格欄中的公式復(fù)制到其他單元格中即可。B1單元格直接賦值0.01，B2單元格公式欄輸入公式：B2+0.001，拖動(dòng)單元格右下角的填充柄，將單元格欄中的公式復(fù)制到其他單元格中即可。C1，……、D1，……單元格同理可得。

F20單元格公式欄輸入公式：INDEX(B1:B900,MATCH(F18,A1:A900,))，F(xiàn)21單元格公式欄輸入公式：INDEX(D1:D900,MATCH(G18,C1:C900,))。

完成以上步驟后，徑流系數(shù)C不需人工多次(從0.001～0.999中，每次以0.001為步長(zhǎng)取值)重復(fù)試算后才能計(jì)算出結(jié)果，而計(jì)算過(guò)程將自動(dòng)完成，真正達(dá)到事半功倍的目的。

4 “林平一法小流域洪水計(jì)算”Excel程序在工程水文計(jì)算中的應(yīng)用

某“220kV變電站”站址西北側(cè)存在一匯水面，其流域面積為0.5km2，主河道長(zhǎng)865.4m，坡度為116.8‰，24h降雨量為146.63mm。采用本Excel“林平一法小流域洪水計(jì)算程序”計(jì)算，得該變電站站址西北側(cè)頻率1%的山洪流量為14.3m3/s，洪水總量2.27×104m3，計(jì)算結(jié)果詳見(jiàn)圖2。

取同樣的參數(shù)，采用《工程水文氣象分析計(jì)算制圖系統(tǒng)》(EHP v2.0)之小面積計(jì)算程序計(jì)算出的結(jié)果基本一致。

綜上所述，該Excel“林平一法”小流域洪水計(jì)算程序，以Microsoft Excel 軟件為基礎(chǔ)編制，窗口界面簡(jiǎn)潔、大方，且簡(jiǎn)單、實(shí)用。用此Excel計(jì)算程序進(jìn)行小流域暴雨洪水計(jì)算，方法簡(jiǎn)便，易于掌握，在實(shí)際查算中，一般5min以內(nèi)可以查算完一個(gè)工程，工作效率大大提高，查算成果以圖、表輸出，整潔美觀，特別是對(duì)分塊流域較多的工程進(jìn)行計(jì)算時(shí)，應(yīng)用本查算軟件優(yōu)勢(shì)更為明顯。

圖2 “林平一法小流域洪水計(jì)算”Excel程序計(jì)算結(jié)果

5 結(jié)語(yǔ)

本文給出的Excel程序通過(guò)大量試算，考慮到了各種不同類型的數(shù)據(jù)，沒(méi)有發(fā)現(xiàn)明顯的問(wèn)題。當(dāng)然，本人由于水平有限，而且只是筆者的一點(diǎn)經(jīng)驗(yàn)之談，旨在拋磚引玉，但愿能讓廣大水文專業(yè)技術(shù)人員從中得到啟發(fā)，編寫(xiě)山更多更好的應(yīng)用程序來(lái)為水文信息化建設(shè)服務(wù)，歡迎大家提出寶貴意。

[1]電力工程水文勘測(cè)計(jì)算手冊(cè)編寫(xiě)組．電力工程水文勘測(cè)計(jì)算手冊(cè)[K]．1981．

[2]云南省水文總站革命委員會(huì)．云南省水文手冊(cè)[K]．1971．

[3]云南省水利水電廳．云南省暴雨洪水查算實(shí)用手冊(cè)[K]．1992．

[4]劉炳文．Visual Basic程序設(shè)計(jì)教程[M]．清華大學(xué)出版社，2006．

Application of Excel Function to Small Basin Flood Calculation

XIE Yin-chang
(Yunnan Provincial Electric Power Design Institute, Kunming 650011, China)

In power engineering design of hydrological calculation work, we often need to calculate the design fl ood, the conventional manual method, calculation cumbersome, heavy workload. The application of Excel function calculation of small basin flood peak flow draw the flood hydrograph, only through some simple settings, can be completed. This method is simple and practical.

Excel function; lin ping method; small Basin— fl ood.

TV12

B

1671-9913(2012)02-0030-04

2012-02-29

謝銀昌(1981- )，男，白族，云南人，主要從事電力工程水文氣象設(shè)計(jì)工作。