張志雁

摘要：應用EXCEL VBA編制的程序，通過讀取規(guī)劃灌區(qū)的灌溉制度數(shù)據(jù)，實現(xiàn)自動繪制灌區(qū)年度初步灌水率圖，并完成動態(tài)修正功能。該方法能夠快速地解決繪制灌水率圖中煩瑣、重復的規(guī)劃設計工作，從而提高工作效率。

關鍵詞：編程；灌水率圖；EXCEL VBA

中圖分類號：TP311 文獻標識碼：A 文章編號：1009-3044（2017）01-0098-02

Abstract： Figure of irrigation rate could be Automatic drawn and dynamic amended by getting the data of irrigation schedule using the programming based on EXCEL VBA. This method can solve problem of planning， thereby improving efficiency of design work.

Key words： programming； figure of irrigation rate； EXCEL VBA

1 概述

灌水率是指灌區(qū)單位面積上所需灌溉的凈流量，它是根據(jù)灌溉制度確定的。將所得的灌水率繪成直方圖，稱為灌水率圖。設計灌水率是設計灌區(qū)內渠首的引水流量和灌溉渠道的設計流量的重要參數(shù)，它的大小直接影響到灌溉系統(tǒng)的工程規(guī)模以及水資源利用充分與否，因此，在灌區(qū)規(guī)劃設計中受到格外的關注。

傳統(tǒng)作物灌水率圖的繪制通用性較差、效率較低，而利用EXCEL VBA編寫程序可有效提高灌區(qū)的規(guī)劃設計效率。本文主要介紹基于EXCEL VBA在灌區(qū)規(guī)劃設計中作物灌水率圖繪制程序的開發(fā)及應用。

2程序設計

2.1設計原理

在灌區(qū)規(guī)劃設計中，影響作物灌水率主要因素為作物種植比例、灌水定額、灌水延續(xù)時間等基本參數(shù)，而這些灌溉制度參數(shù)通常是可確定的。

編程過程中，通過讀取確定格式的灌溉制度參數(shù)，計算各種作物播前灌水及生育期內各次灌水的灌水率，通過編程實現(xiàn)輸出全年各灌水時間（以天為單位）內對應的作物灌水率，并應根據(jù)每次灌水延續(xù)時間繪制各種作物的灌水率過程線，將同時期各種作物灌水率相加，繪成設計灌區(qū)年度初步灌水率圖，并完成修正。

2.2程序實現(xiàn)

1）讀取灌溉制度數(shù)據(jù)

程序中灌溉制度設定為固定格式，見表1。灌溉制度中，作物數(shù)量、各作物灌水次數(shù)等參數(shù)為不確定數(shù)值。

J = 1

For I = 1 To 200

If Cells（I + 1， 4） = "" Then Exit For

If Cells（I + 1， 1） <> ""） Then

作物類別（J） = Cells（I + 1， 1）

灌溉面積（J） = Cells（I + 1， 2）

作物比例（J） = Cells（I + 1， 3）

灌溉定額（J） = Cells（I + 1， 7）

For K = 1 To 200

If （K > 1 And Cells（I + K， 1） <> "" Or Cells（I + K， 6） = ""） Then Exit For

灌水起始時間（J， K） = Cells（I + K， 4）

灌水天數(shù)（J， K） = Cells（I + K， 5） -Cells（I + K， 4） + 1

灌水定額（J， K） =Cells（I + K， 6）

Next K

作物數(shù)量 = J

J = J + 1

End If

Next I

2）計算作物灌水率

按照《灌溉與排水工程設計規(guī)范》公式計算作物灌水率。

For I = 1 To 作物數(shù)量

For J = 0 To 200

If （灌水天數(shù)（I， J） = 0） Then Exit For

灌水率（I， J） = 作物比例（I） * 灌水定額（I， J） / 灌水天數(shù)（I， J） / 864

Next J

Next I

3）繪制初步灌水率圖

在“灌水率”工作表中，按照灌水日期與全年日期對應匹配輸出各作物的灌水率，然后完成后繪制初步灌水率圖。

For I = 1 To 作物數(shù)量

Cells（2， I + 2） = 作物類別（I）

For J = 1 To 200

If （灌水天數(shù)（I， J） = ""） Then Exit For

For K = 0 To 灌水天數(shù)（I， J） - 1

Cells（灌水起始時間（I， J） + 2 + K， I + 2） = 灌水率（I， J）

Next K

Next J

Next I

S1 = Chr$（66 + 作物數(shù)量）

Range（"B2：" & S1 & "368"）.Select

ActiveSheet.Shapes.AddChart.Select

ActiveChart.ChartType = xlColumnStacked

4）修正灌水率圖

在符合規(guī)范要求的情況下，可通過提前或推遲灌水日期、延長或縮短灌水時間、調整灌水定額等方法對繪制的初步灌水率圖進行修正，并實現(xiàn)灌水率表和灌水率圖的動態(tài)調整，進而獲得設計灌水率圖。

5）計算作物月凈耗水量

同時，在灌區(qū)規(guī)劃設計中，需進行水資源供需平衡分析評價，利用該程序可快速計算出各作物月凈耗水量，可有效提高計算效率。其方法為，在步驟3中循環(huán)語句內添加如下語句：Cells（灌水起始時間（I， J） + 2 + K， I + 2） = 灌水定額（I， J） / 灌水天數(shù)（I， J），該語句用于計算、輸出各作物在不同日期的灌水定額，然后按以下程序可計算各作物月凈耗水量，并在“月凈耗水量統(tǒng)計”工作表中實現(xiàn)數(shù)據(jù)輸出功能。

For I = 1 To 作物數(shù)量

s = 0

For J = 0 To 11

月灌水量（I， J） = 0

For K = 0 To Months（J） - 1

月灌水量（I， J） = 月灌水量（I， J） + Cells（3 + K + s， I + 2）

Next K

s = s + Months （J）

Next J

月灌水量（I， J - 1） = 月均灌水量（I， J - 1）

Next I

For I = 1 To 作物數(shù)量

Cells（2 + I， 1） = 作物類別（I）

Cells（2 + I， 2） = 灌溉面積（I）

For J = 0 To 11

Cells（2 + I， 3 + J） = 月灌水量（I， J） * 灌溉面積（I）

Next J

Next I

3算例分析

以南疆某中型灌區(qū)規(guī)劃設計為例，介紹運用該程序繪制灌水率圖。該灌區(qū)種植作物主要有果園、棉花、小麥等7種，灌溉模式有節(jié)水灌和地面灌。輸入灌溉制度后，直接單擊程序設計的按鈕命令，即可運行，繪制出初步灌水率圖，并通過提前或推遲灌水日期，動態(tài)調整確定設計灌水率圖，同時獲得各作物月凈耗水量。最終的設計灌水率圖見圖1。

4結束語

本程序通過EXCEL VBA編制程序讀取EXCEL表中的灌溉制度數(shù)據(jù)，計算各種作物播前灌水及生育期內各次灌水的灌水率，輸出全年各灌水時間內對應的作物灌水率，自動繪制灌區(qū)年度初步灌水率圖，并實現(xiàn)動態(tài)修正功能，同時，計算出作物月凈耗水量。該方法能夠快速地解決灌水率圖繪制中煩瑣、重復的規(guī)劃設計工作，通用性強，可極大地提高工作效率。

參考文獻：

[1] GB 50288-99灌溉與排水工程設計規(guī)范[S] .北京：中國計劃出版社， 1999.

[2] 水利部水利水電規(guī)劃設計總院.水工設計手冊第九卷灌溉、供水[M] . 2版.北京：中國水利水電出版社，2014.

[3] 汪志農.灌溉排水工程學[ M] .北京：中國農業(yè)出版社， 2003.

[4] 龔沛曾.Visual Basic程序設計簡明教程[M].2版.北京：高等教育出版社，2006.

[5] Excel VBA程序開發(fā)自學寶典[M]. 3版.北京：電子工業(yè)出版社，2014.

[6] 邱衛(wèi)國，房寬厚. 修正灌水率圖的新方法[J]. 河海大學學報，1995（5）.

[7] 張宇亮，張禮兵，周玉良，等. 基于改進加速遺傳算法的作物灌水率圖修正研究[J]. 灌溉排水學報， 2015（11）

[8] 楊強勝，張禮兵，陳得陽，等. 基于EXCEL的作物灌水率圖自動繪制與修正[J]. 中國農村水利水電，2011（6）.