渠系水利用系數測算方法研究

謝亨旺，羅云英，靳偉榮

(1.江西省灌溉試驗中心站，南昌 330201；2. 暨南大學環境學院，廣州 511443；3.武漢大學水利水電學院，武漢 430072)

0 引 言

灌溉水有效利用系數是最嚴格水資源管理的指標之一，也是農業灌溉用水管理、灌溉技術水平評價、區域水資源配置和制定節水灌溉發展規劃的基礎參數[1]。渠系水利用系數是灌溉水有效利用系數的重要組成部分，其測算方法的研究有助于正確認識灌溉水利用系數指標及內涵。目前測算渠系水利用系數國外常用的經驗公式有Davis-Wilson公式、USBR公式、Ingham公式、Molesworth公式、前蘇聯公式和Kositiakov經驗公式等[2,3]。我國廣泛采用Kositiakov經驗公式的改進公式[4]，對土壤透水性能經驗常數做了細化規定，但由于各地區普遍缺少渠道狀況的實測參數資料，使該公式在實際應用中存在較多問題。榮豐濤、江崇安等分析了現行渠道滲漏量計算方法的不足，如滲水量減少系數、防滲材料滲透系數等難以確定[5,6]；門寶輝、鐘玲等利用積分方法推導了新公式，對傳統渠道損失流量計算公式的不足進行了修正，提高了渠道損失流量和渠道水利用系數的計算精度[7,8]；白美健等采用回歸分析方法建立了渠道滲漏損失水量與流量間的相關關系，總結出了依據渠道流量估算渠道水利用系數的經驗公式[9]。本文應用渠道連續滲漏理論，對傳統計算渠道單位長度入滲損失系數的Kostiakov經驗公式進行改進，并在江西省贛撫平原灌區進行驗證。改進的計算方法將渠道流量損失過程看成水在沿渠道流動方向的連續動態過程，符合水流損失規律，但相關計算精度仍然受到防滲渠道損失修正系數、地下水頂托修正系數等參數的影響。這些數據需要大量的試驗資料進行率定。

1 國內常用的Kostiakov經驗公式及其改進

1.1 國內常用的Kostiakov經驗公式

我國在計算渠道滲漏損失的時候多以Kostiakov經驗公式為基礎，并對反映土壤透水性能的經驗常數A、m進行了細化，也考慮了渠道襯砌方式及地下水位對渠道滲水損失的影響。常用的的經驗公式形式為：

σ=0.01AQ-m

(1)

S=0.01γβAQ1-m

(2)

式中：σ為每公里渠道輸水損失率，%；Q為渠道流量，m3/s；A、m為渠床土壤透水參數；S為每公里渠道的渠道輸水損失流量，m3/s；σ為地下水頂托修正系數；β為采取防滲措施后渠道滲水損失修正系數。

在《灌溉與排水工程設計規范》(GB 50288-99)和《農田水利學》[4]中，式(2)中的A、m、γ、β的取值分別見表1、表2和表3。

表1 Kostiakov公式土壤透水參數A、m

表2 地下水頂托修正系數γ

表3 采取防滲措施渠道滲水損失修正系數β經驗取值

已知某類型渠道的長度 (單位：km)、渠道毛流量為Qg、渠道凈流量為Qn。當灌區地下水位較高時，地下水壅阻會減少渠道滲漏，用系數γ修正；當渠道采取了防滲措施時，不同防滲措施會對渠道滲水產生不同程度的降低折減，用系數β進行修正。綜合考慮地下水和防滲措施的影響，渠道滲漏水量可按下式計算：

ΔQ=βγσLQn

(3)

該段渠道的渠道水利用系數為：

(4)

1.2 Kostiakov經驗公式的改進

Kostiakov經驗公式計算簡便，在我國應用廣泛，但它還存在一些明顯的不足和局限性。作為估算公式，其估算精度取決于參數A、m、γ、β的取值是否合理。單位渠道輸水損失率σ是渠道流量的函數，而渠道流量沿程變化，一般按式(3)～(4)計算時，假定了σ不隨沿程渠道流量變化而變化，這使得計算渠道越長誤差越大。李淑榮通過論述多種求解渠系水利用系數思路，探索推導出利用系數法、損失系數法及定義法的誤差公式，在此基礎上揭示了誤差產生的原因，損失系數法在計算過程中忽略了沿程渠道流量變化[10]。

榮豐濤分析了Kostiakov經驗公式未考慮渠道長度對計算結果的影響，用現行公式計算了渠道滲漏損失相對誤差小于5%時的最大允許計算渠長，結果表明，在A、m、β一定的情況下，最大允許渠長隨著設計流量Q的增加而增加；在其他參數不變的情況下，渠道防滲效果越好(β越小)，最大允許渠長越長[5]。

根據公式(2)可以發現，Kostiakov經驗公式求得的每公里渠長的輸水損失是一樣的，而實際上每公里渠長上的損失流量隨流量變化，越往渠尾處流量越小，相應的滲水損失也越小，渠首段單位渠長的流量損失一般大于渠尾段。針對這種情況，門登輝等學者提出采用積分方法對渠道設計中的流量損失計算公式進行改進[7-9]。設渠首毛流量為Qg，經過流程L后凈流量為Qn。設渠道流量Q在經過渠段dl后的損失流量為dQ。在dl渠段內，單位流量在單位流程上的損失σ=dQ/(Qdl)，則：

(5)

將上式(5)在渠長L范圍內進行積分求解得：

(6)

(7)

(8)

該計算公式將渠道流量損失過程看成水在沿渠道流動方向的連續動態過程，符合水流損失規律，會顯著提高估算精度。

2 改進公式的應用效果

將改進的公式與一般使用的Kostiakov經驗公式進行比較分析，求解不同情況下的相對誤差，相對誤差 用下式進行表示：

(9)

式中：ΔQ改進為積分方法改進公式推求的渠道流量損失值，m3/s；ΔQ為Kostiakov公式推求的渠道流量損失的近似值，m3/s。

取A=1.9、m=0.4，不考慮頂托和渠道防滲襯砌措施。分別采用Kostiakov經驗公式和積分方法改進的公式計算不同流量和渠長組合情況下的渠道輸水損失量，計算結果見表4。

表4 改進前后渠道損失流量對比

由表4可知，2種公式損失流量的計算結果均與渠道凈流量和渠道長度有關，無論是Kostiakov經驗公式計算還是積分改進公式計算，渠道流量和長度均與渠道損失流量成正相關關系。渠道凈流量一定時，2種方法的計算結果在較小的渠道長度范圍內，相對誤差很小，計算結果都可以為工程采用。但是隨著渠道長度的增加，Kostiakov經驗公式計算結果與積分改進公式計算結果的相對誤差逐漸增大。

3 改進經驗公式在灌區進行驗證

3.1 改進經驗公式與實測數據比較分析

本文采用動水測定法的實測數據對經驗公式的適用性進行驗證。

實際測定工作在位于江西省贛撫平原灌區二干渠上進行。選擇二干二支渠、五支渠為典型支渠，二支三斗渠、五支五斗渠為典型斗渠，典型支、斗渠所處位置見圖1、圖2。對于二支三斗渠，根據《基本農田建設設計規范》布置20條農渠，受地形限制，每條農渠的灌溉面積不盡相等；根據農渠控制灌溉面積大致相等分2組輪灌(三斗渠上左、右前6條農渠為第1輪灌組，7～10條為第2輪灌組)，渠系概化圖見圖1。在五支五斗渠上布置了20條農渠，分2組輪灌(斗渠上左、右1～5條農渠為一組，6～10為一組)，渠系分布情況見概化圖圖2。

圖1 二干渠二支渠概化圖

圖2 二干渠五支渠概化圖

根據實地調查，研究區干渠、支渠全部為混凝土襯砌，斗渠、農渠一半為襯砌渠道一半為土渠，斗渠防滲形式為混凝土護面，農渠為U形槽。計算襯砌渠道時，假設襯砌優先從渠首位置開始。采取防滲措施的各級渠道都采用混凝土材料，參考表3可知，混凝土護面的襯砌折減系數β為0.05～0.15。本文在計算過程中，對采用防滲措施的干渠和支渠的β取中間值0.1；由于斗渠和農渠襯砌厚度較薄且日常維護較少，所以取相應的襯砌渠道β值為0.15。二干渠灌溉區域的渠道工作制度為干渠、支渠、斗渠續灌，農渠輪灌。

3.1.1 農渠渠道水利用系數

采用田間水利用系數實測結果0.917 6推算農渠田間凈流量：

(10)

式中：A農為農渠灌溉面積，萬hm2；Q農凈為農渠供給田間的凈流量，m3/s；q設為設計凈灌水率，m3/(s·hm2)。

根據公式(6)有：

(11)

式中：Q農凈為農渠的凈流量，m3/s；Q農毛為農渠的毛流量，m3/s；研究區域土質為輕壤土，參考表1，取A=2.65，m=0.45。

由Cropwat 8.0軟件計算得到研究區域1951-2013年騰發量ET，進而推求作物灌溉制度和灌溉定額，根據種植比例求得研究區年灌溉用水定額，由1963年灌溉用水定額排頻繪制P-Ⅲ曲線，得灌溉用水保證率為90%時設計典型年為1966年，繪制灌水率圖并修正，得設計凈灌水率為0.000 705 m3/(s·hm2)。

根據動水測定法計算得到的典型渠段上下游流量Q典1和Q典2，就可以求得典型渠段的水利用系數η典段：

η典段=Q典1/Q典2

(12)

動水測定法渠道輸水損失系數改進公式為：

(13)

由公式(12)、(13)得每條農渠毛流量和渠道水利用系數，由控制灌溉面積加權平均得典型農渠平均渠道水利用系數為0.933，將其視為研究區域農渠渠道水利用系數，二支渠計算結果見表5。

3.1.2 斗渠渠道水利用系數

由于農渠輪灌，以農渠輪灌組最大流量值作為斗渠凈流量，二支渠三斗渠結果見表5。

表5 二支渠農渠渠道水利用系數和三斗渠凈流量計算結果

典型斗渠長度為1.68 km，自渠首起有一半長度采用混凝土材料襯砌，根據公式(6)得斗渠毛流量為0.134 3 m3/s，最終得到斗渠渠道水利用系數為0.964 6。以斗渠控制面積的灌溉水利用系數0.794為擴大指標來推求其他斗渠毛流量，進而得到二支渠凈流量為0.509 2 m3/s，計算結果見表6。

表6 斗渠毛流量和支渠凈流量計算結果

3.1.3 支渠渠道水利用系數。

典型支渠全部為混凝土材料襯砌，根據公式(12)求得支渠毛流量為0.516 8 m3/s，最終得到支渠渠道水利用系數為0.985 2，計算結果見表7。

表7 支渠毛流量計算結果

3.1.4 干渠渠道水利用系數

干渠全部為混凝土襯砌，以支渠控制面積的灌溉水利用系數0.782為擴大指標來推求其他支渠毛流量，進而得到二干渠凈流量為8.683 6 m3/s，結果見表8。

求得二干渠的渠道水利用系數為0.984 3：

對于五支渠，同理推導得出各級渠道水利用系數，由面積加權平均得研究區域的各級渠道水利用系數和相應的渠系水利用系數，計算結果見表9。

此外，也可以由推算得到的凈流量與毛流量的比值來計算渠系水利用系數(見表10)，二支渠和五支渠結果按面積加權平均后為0.839。

表8 二干渠流量推算結果

表9 改進經驗公式法的渠道水利用系數計算結果

表10 由渠道凈流量與毛流量比值推求的渠系水利用系數結果

3.2 小 結

改進經驗公式法的計算中，由公式計算的各級渠道水利用系數連乘求得的渠系水利用系數為0.873，較由級渠道凈流量與逐級推算得到的干渠渠首流量(毛流量)的比值得到的渠系水利用系數結果0.839有所不同。由于經驗公式法計算過程中采用某級渠道的灌溉水利用系數作為擴大指標來計算同級其他渠道流量不能真實反映整個渠系情況，較實際情況存在偏差，因此從渠系水利用系數定義出發結果0.839更合理些。同時在用改進經驗公式計算時，考慮了襯砌和地下水影響時，混凝土襯砌材料對滲漏水量的折減比例很大(根據表1，折減了約90%)。根據實地調研，灌區渠道防滲材料遭到不同程度的破壞，渠道大多存在不同程度的漏水情況，但經驗公式中的相關參數未及時修正，因此改進經驗公式求得的結果反映了渠道在某種襯砌狀況下的渠系水利用系數。

4 結 論

(1)為解決現行渠道利用系數測算工作量大、影響因素復雜、測算精度難以控制的問題，根據渠道連續滲漏理論分別對傳統計算渠道單位長度入滲損失系數的Kostiakov經驗公式進行了改進，提出了改進單位長度渠道損失流量的計算公式，公式如下：

(2)以贛撫平原灌區二干渠土壤、襯砌材料等為背景，對渠道輸水凈流量為0.05～20.00 m3/s，渠道長度為0.05～40.00 km的測量和分析表明，改進Kostiakov經驗公式與原公式相比，2者的相對誤差為0.01%～46.00%，渠道的流量越小、渠道越長，Kostiakov經驗公式計算結果與積分改進公式計算結果的相對誤差逐漸增大。

(3)以贛撫平原灌區二干渠的實測資料為背景，改進經驗公式計算得出渠系水利用系數為0.873，較由級渠道凈流量與逐級推算得到的干渠渠首流量(毛流量)的比值得到的渠系水利用系數結果0.839有所不同，由于改進經驗公式法計算過程中采用某級渠道的灌溉水利用系數作為擴大指標來計算同級其他渠道流量不能真實反映整個渠系情況，較實際情況存在偏差，因此從渠系水利用系數定義出發結果0.839更合理些。根據實地調研，灌區渠道防滲材料遭到不同程度的破壞，渠道大多存在不同程度的漏水情況，但經驗公式中的相關參數未及時修正，因此改進經驗公式求得的結果反映了渠道在某種襯砌狀況下的渠系水利用系數。