劉凱華，繳錫云，劉 懿

(1. 河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，南京 210098；2. 河海大學(xué)水利水電學(xué)院，南京 210098；3. 海南省三亞市水利水電勘測(cè)設(shè)計(jì)院，海南 三亞 572000

在地面灌溉過(guò)程中，由于實(shí)際灌水流量不易控制，所以相較于設(shè)計(jì)灌水流量有偏差，以致不能得到理想的灌水質(zhì)量，甚至可能出現(xiàn)灌水質(zhì)量不合格的情況[1-4]。為了克服流量及其他因素的時(shí)空變異性，王維漢進(jìn)行溝畦灌的穩(wěn)健設(shè)計(jì)[5]。進(jìn)行穩(wěn)健設(shè)計(jì)工作的基礎(chǔ)之一就是對(duì)灌水流量等灌水技術(shù)要素的控制誤差進(jìn)行調(diào)查。目前，流量的精確測(cè)量方法有很多，包括流量計(jì)法、流速儀法、量水堰法等[6,7]，但對(duì)于以統(tǒng)計(jì)流量控制誤差為目標(biāo)的調(diào)查來(lái)說(shuō)，這些方法均具有難以實(shí)現(xiàn)無(wú)干擾測(cè)量、成本過(guò)高等缺點(diǎn)。示蹤劑稀釋法根據(jù)注入水體里示蹤劑的濃度與混合均勻斷面上示蹤劑濃度的比值推算水體流量[8-10]，可低成本、無(wú)干擾測(cè)量流量，但可能造成局部水質(zhì)污染，且對(duì)施測(cè)水體有較高要求，并不適用于存在明顯入滲現(xiàn)象的灌溉過(guò)程的流量測(cè)量。楊靜等根據(jù)畦首水深和水流推進(jìn)的數(shù)據(jù)，引入地表儲(chǔ)水形狀系數(shù)和地下儲(chǔ)水形狀系數(shù)，進(jìn)而估算畦灌的單寬流量。但這種反演方法的精度不高，相對(duì)誤差在10%左右[11]。為此，有必要提出具有較高精度且方便可行的新的入畦流量反演方法。本文擬利用2013年冬小麥春灌試驗(yàn)得到的諸如水流推進(jìn)消退過(guò)程、灌水前后土壤含水率等數(shù)據(jù)，提出通過(guò)對(duì)控制點(diǎn)的入滲量進(jìn)行積分求得單寬入滲量，再由水量平衡原理得到灌水量的反演方法，從而使大范圍地、準(zhǔn)確地調(diào)查灌水流量控制誤差成為可能，為溝畦灌的穩(wěn)健設(shè)計(jì)奠定基礎(chǔ)。

1 材料與方法

1.1 實(shí)驗(yàn)區(qū)概況

田間試驗(yàn)于2013年4月中旬在中科院南皮生態(tài)農(nóng)業(yè)試驗(yàn)站結(jié)合冬小麥春灌進(jìn)行。試驗(yàn)區(qū)距南皮縣城約5 km，地處北緯38°06′，東經(jīng)116°40′，區(qū)內(nèi)主要種植作物為冬小麥和夏玉米，屬華北平原井灌區(qū)，灌溉方式為畦灌。土壤屬粉壤土，0～100 cm土層的土壤干密度見表1。

表1 0～100 cm深度土層的土壤干密度Tab.1 Dry density of soil at 0～100 cm

1.2 試驗(yàn)方案與觀測(cè)項(xiàng)目

選擇3條畦田進(jìn)行灌水試驗(yàn)，畦田規(guī)格均為長(zhǎng)80 m，寬3.7 m，水源為機(jī)井抽取的地下水。具體灌水試驗(yàn)設(shè)計(jì)見表2。

表2 灌水試驗(yàn)設(shè)計(jì)Tab.2 Design of irrigation experiment

田間試驗(yàn)主要觀測(cè)項(xiàng)目包括：土壤初始含水率、水流推進(jìn)消退過(guò)程、畦首水深、灌水流量及灌后土壤含水率等。土壤含水率的測(cè)量采用取土烘干法。從畦首開始，沿畦長(zhǎng)方向每隔5 m布置1個(gè)觀測(cè)斷面，每個(gè)斷面3個(gè)取土點(diǎn)，每條畦田共布置17個(gè)取土斷面，共計(jì)51個(gè)取土點(diǎn)。每個(gè)取土點(diǎn)觀測(cè)0～5、5～10、10～20、20～40、40～60、60～80、80～100 cm等7個(gè)土層的土壤含水率。

2 入畦流量反演

本文提出了2種基于水量平衡原理的入畦流量反演方法：利用土壤含水率計(jì)算入畦流量的方法(以下簡(jiǎn)稱含水率法)和基于土壤入滲模型計(jì)算入畦流量的方法(以下簡(jiǎn)稱入滲模型法)，并與傳統(tǒng)的基于畦首水深及地表水流運(yùn)動(dòng)過(guò)程計(jì)算入畦流量的方法(以下簡(jiǎn)稱畦首水深法)進(jìn)行了反演精度和工作量等方面的比較。

2.1 含水率法

含水率法首先利用控制點(diǎn)灌水前后土壤含水率的變化求得該點(diǎn)的入滲水深，然后利用辛普森公式積分得到畦田單寬入滲水量，根據(jù)水量平衡方程即可求得入畦單寬流量。運(yùn)用此方法時(shí)，每條畦田上17個(gè)斷面的土壤灌前灌后含水率均需觀測(cè)，對(duì)每個(gè)斷面3點(diǎn)含水率取均值，另需記錄灌水時(shí)間。

儲(chǔ)存于各點(diǎn)土壤計(jì)劃濕潤(rùn)層內(nèi)的水量可表示為：

(1)

式中：V(x)為觀測(cè)點(diǎn)入滲量，cm；θ(z,t)為灌水后的含水率分布；θ(z,0)為灌水前的含水率分布；D為土層厚度，cm，這里取100 cm。

對(duì)各觀測(cè)斷面含水率均值在畦長(zhǎng)方向上采用復(fù)合辛普森公式積分，得到整條畦田的單寬入滲水量。畦田長(zhǎng)80 m，每隔5 m設(shè)立一個(gè)觀測(cè)斷面，從畦首開始觀測(cè)點(diǎn)依次編號(hào)為1、2、3、…、16、17。根據(jù)辛普森公式原理將畦田沿長(zhǎng)度方向按步長(zhǎng)h=10 m分為8等份，在每個(gè)子區(qū)間[xk,xk+2]上采用辛普森公式，k為奇數(shù)觀測(cè)點(diǎn)，即1、3、5、…、13、15；xk為第k個(gè)觀測(cè)點(diǎn)距畦首的距離。若記xk+1=xk+h/2，則得：

(2)

式中：I為畦田單寬入滲量，cm2；V(xk)為第k個(gè)觀測(cè)點(diǎn)入滲量，cm；其余符號(hào)意義同前。

根據(jù)水量平衡原理，總的入滲量即為灌水量，忽略畦寬的空間變異性，則單寬入滲量即為單寬灌水量。入畦單寬流量為：

q=I/t

(3)

式中：q為入畦單寬流量，cm2/s；I為畦田單寬灌水量，cm2；t為灌水時(shí)間，s。

2.2 入滲模型法

入滲模型法首先利用隨機(jī)取的畦田上2個(gè)觀測(cè)斷面灌水前后土壤含水率的變化求得Kostiakov模型參數(shù)，利用Kostiakov模型計(jì)算得到點(diǎn)入滲量，積分得畦田單寬入滲量，進(jìn)而得到入畦單寬流量。此方法需要觀測(cè)的項(xiàng)目包括水流推進(jìn)消退資料，沿畦長(zhǎng)方向2個(gè)觀測(cè)斷面的灌前、灌后土壤含水率值。

Kostiakov入滲模型屬于經(jīng)驗(yàn)?zāi)Ｐ?，被廣泛應(yīng)用。該模型描述的土壤累積入滲量和入滲時(shí)間的關(guān)系如下：

I=ktα

(4)

式中：I為土壤累計(jì)入滲量，mm；k為土壤入滲系數(shù)，mm/minα；α為入滲指數(shù)，無(wú)量綱；α和k統(tǒng)稱為入滲參數(shù)。

土壤入滲參數(shù)估算方法較多，但大多數(shù)方法都需要用到入畦流量，本文采用的是繳錫云介紹的灌溉前后利用土壤含水率分布資料估算畦田入滲參數(shù)[12]，此方法僅需灌溉前后觀測(cè)點(diǎn)的土壤含水率以及推進(jìn)消退資料，計(jì)算結(jié)果能真實(shí)反應(yīng)田面水流運(yùn)動(dòng)對(duì)入滲的影響。

采用Kostiakov模型描述土壤入滲過(guò)程，則觀測(cè)點(diǎn)處入滲的水量可表示為：

Vi=k(tr-ta)α

(5)

式中：tr為該點(diǎn)處的消退時(shí)間，min；ta為水流前鋒推進(jìn)至該點(diǎn)處的時(shí)間，min；其余符號(hào)意義同前。

對(duì)田間任意2個(gè)觀測(cè)斷面，建立方程入滲參數(shù)k和α計(jì)算公式如下：

(7)

式中：觀測(cè)斷面推進(jìn)時(shí)間(ta1，ta2)、消退時(shí)間(tr1，tr2)和灌水前后含水率[θ1(z,0)、θ1(z,t)、θ2(z,0)、θ2(z,t)]均通過(guò)試驗(yàn)測(cè)得；D為取土層深度，cm。

根據(jù)求得的土壤入滲參數(shù)和觀測(cè)得到的各點(diǎn)的入滲時(shí)間，可得到沿畦長(zhǎng)方向各觀測(cè)點(diǎn)的入滲量。同樣利用辛普森公式積分求出整條畦田的入滲量，根據(jù)水量平衡原理求得入畦單寬流量：

(8)

式中：q為入畦單寬流量，cm2/s；I為畦田單寬灌水量(即單寬入滲量)，cm2；t為灌水時(shí)間，s；L為畦長(zhǎng)，m。

2.3 反演方法的比較

在大田試驗(yàn)過(guò)程中，利用水表與秒表觀測(cè)了流量，以此為基準(zhǔn)分析含水率法、入滲模型法以及傳統(tǒng)的畦首水深法的反演精度，具體見表3。各反演入畦流量方法觀測(cè)的內(nèi)容不盡相同，具體每種方法需要觀測(cè)的內(nèi)容見表4。

表3 反演精度比較Tab.3 Comparison of inversion accuracy

由表3可以看出，入畦流量的反演方法中，畦首水深法、含水率法和入滲模型法的平均相對(duì)誤差分別為10.84%、5.81%和5.83%。含水率法和入滲模型法明顯具有較高的反演精度，且兩者之間的差異不大。

表4 工作量比較Tab.4 Comparison of workload

入畦流量的反演方法中，含水率法和入滲模型法均具有較高精度，但是含水率法需要觀測(cè)714個(gè)土樣，遠(yuǎn)遠(yuǎn)大于入滲模型法和畦首水深法的84個(gè)土樣，工作量偏大。

3 結(jié) 語(yǔ)

針對(duì)入畦流量的調(diào)查，本文提出了含水率法和入滲模型法2種反演入畦流量的方法。含水率法利用控制點(diǎn)灌水前后土壤含水率的變化求得該點(diǎn)的入滲水深，進(jìn)而利用辛普森公式積分得到畦田單寬灌水量。入滲模型法利用Kostiakov模型計(jì)算得到點(diǎn)入滲量，積分得畦田單寬入滲水量，進(jìn)而得到入畦單寬流量。與傳統(tǒng)的反演方法相比，含水率法和入滲模型法均具有更高的精度，且入滲模型法需要更少的觀測(cè)量。因此，一般宜選用入滲模型法反演入畦流量。

□

[1] Trout T J, Mackey B E. Furrow inflow and infiltration variability[J]. Transactions of the ASAE, 1988,31(2):531-537.

[2] Mailhol J C, Priol M, Benali M. A furrow irrigation model to improve irrigation practices in the Gharb valley of Morocco[J]. Agricultural Water Management, 1999,42(1):65-80.

[3] Schwankl L J, Raghuwanshi N S, Wallender W W. Furrow irrigation performance under spatially varying conditions.[J]. Journal of Irrigation & Drainage Engineering, 2000,126(6):355-361.

[4] 陳云武，朱春光. 基于Monte Carlo模擬流量變異性對(duì)灌水質(zhì)量影響的評(píng)價(jià)[J]. 安徽農(nóng)業(yè)科學(xué), 2011,39(36):22 849-22 851.

[5] 王維漢. 畦灌影響因素變異規(guī)律及灌水技術(shù)要素穩(wěn)健設(shè)計(jì)[D]. 南京:河海大學(xué), 2009.

[6] Samani Z, Magallanez H. Simple Flume for Flow Measurement in Open Channel[J]. Journal of Irrigation & Drainage Engineering, 2000,126(2):127-129.

[7] 王成李. 明渠流量測(cè)量方法及裝置研究[D]. 杭州:中國(guó)計(jì)量學(xué)院, 2014.

[8] 韓惠蘭，孫新熙，孟慶元，等. 食鹽電導(dǎo)稀釋法測(cè)試明渠水流量技術(shù)[J]. 環(huán)境工程，1990,8(1):51-55.

[9] Lei T, Xia W, Zhao J, et al. Method for measuring velocity of shallow water flow for soil erosion with an electrolyte tracer[J]. Journal of Hydrology, 2005,301(1-4):139-145.

[10] 董月群. 溶質(zhì)示蹤測(cè)量徑流流速方法及其在融水侵蝕中的應(yīng)用研究[D]. 北京：中國(guó)農(nóng)業(yè)大學(xué), 2014.

[11] 楊 靜，繳錫云，劉 懿. 基于畦首水深和田面水流推進(jìn)過(guò)程的畦灌流量反演方法[J]. 水資源與水工程學(xué)報(bào)，2016,27(1):238-240.

[12] 繳錫云，王維漢. 溝畦灌溉穩(wěn)健設(shè)計(jì)[M]. 南京: 河海大學(xué)出版社, 2012.