塔里木灌區(qū)膜下滴灌棉花水分生產(chǎn)函數(shù)及其效益

閆 映 宇

(新疆維吾爾自治區(qū)水土保持生態(tài)環(huán)境監(jiān)測總站, 新疆 烏魯木齊 830000)

?

塔里木灌區(qū)膜下滴灌棉花水分生產(chǎn)函數(shù)及其效益

閆 映 宇

(新疆維吾爾自治區(qū)水土保持生態(tài)環(huán)境監(jiān)測總站, 新疆 烏魯木齊 830000)

以2007—2012年塔里木灌區(qū)膜下滴灌田間灌溉試驗數(shù)據(jù)為基礎，采用最小二乘法原理，擬合了棉花的水分生產(chǎn)函數(shù)模型，分析并計算了棉花最高產(chǎn)量灌溉定額、最佳效益灌溉定額、高效用水灌溉定額，揭示了棉花的水分效應及需水規(guī)律。結果表明：塔里木灌區(qū)膜下滴灌棉花需水臨界期是花鈴期和蕾期；產(chǎn)量與耗水量、灌水量均呈良好的二次拋物線關系；合理灌溉定額為3 091～3 464 m3·hm-2，最高產(chǎn)量灌溉定額為3 464 m3·hm-2，高效用水灌溉定額為3 091 m3·hm-2；水資源投入的最佳效益點并非水分利用效率最高點和最高產(chǎn)量點，而是存在于3 091～3 464 m3·hm-2區(qū)間；當邊際效益等于邊際成本時，凈收益最大，為24 333.1 元·hm-2，每立方灌水量凈收益7.23 元·hm-2；現(xiàn)狀條件下最佳效益灌溉定額為3 459 m3·hm-2，產(chǎn)量為6 360.7 kg·hm-2，與最高產(chǎn)量6 360.8 kg·hm-2基本相同，但比最高產(chǎn)量節(jié)水5 m3·hm-2，每立方灌水凈收益增加0.21 元·hm-2，水分利用效率提高0.003 kg·m-3。

膜下滴灌；棉花；水分生產(chǎn)函數(shù)；水分利用效率；灌水效益

棉花是塔里木灌區(qū)主要經(jīng)濟作物之一，由于塔里木灌區(qū)屬于資源型極度缺水地區(qū)，而棉花又是高耗水作物，區(qū)域水資源不足，直接影響著棉花生產(chǎn)的發(fā)展。為解決水資源短缺的矛盾，新疆大力發(fā)展節(jié)水灌溉技術，提高水資源利用效率。膜下滴灌是新疆綠洲農(nóng)業(yè)最為有效的節(jié)水灌溉技術之一，在棉花種植中得到了大面積推廣。近年來廣大學者針對滴灌條件下土壤水鹽運動規(guī)律及影響因素、膜下滴灌條件下棉花根系生長、耗水規(guī)律等方面進行了深入研究[1-13]，為從根本上探討水資源的最合理利用方式，確立最優(yōu)灌溉制度，著力提高水資源利用效率提供了依據(jù)。

傳統(tǒng)的灌溉目標是向作物提供適宜的水分以獲得高的單位面積產(chǎn)量，隨著水資源日益緊缺，灌溉費用的增加，要求灌溉不僅要使作物高產(chǎn)，而且要獲得最優(yōu)的經(jīng)濟效益。作物水分生產(chǎn)函數(shù)是進行限水灌溉和有限水資源優(yōu)化調配的基礎。國內外有關水分生產(chǎn)函數(shù)的研究有很多[14-20]，從不同層次和不同范疇研究了作物的水分利用效率，并建立了各種水分生產(chǎn)函數(shù)模型，目的主要是通過模型找到作物需水的關鍵期，以使有限的水資源發(fā)揮更大的作用。一般資源投入與農(nóng)業(yè)產(chǎn)生都存在報酬遞減現(xiàn)象，但作物產(chǎn)量和耗水量、水分利用效率之間存在著較為復雜的關系，至今未能得到很好的解決[17]。目前作物水分生產(chǎn)函數(shù)可分為兩類：一是反映作物產(chǎn)量與全生育期總耗水量關系的水分生產(chǎn)函數(shù)稱為總水分生產(chǎn)函數(shù)；二是反映作物產(chǎn)量與各生育階段耗水量組合之間關系的水分生產(chǎn)函數(shù)稱為階段水分生產(chǎn)函數(shù)。我國最常用的模型有相加(Blank)模型和相乘(Jensen)模型，對于作物水分生產(chǎn)函數(shù)的研究主要集中在對各種模型的計算方法、建摸和水分脅迫對產(chǎn)量影響的分析上，針對棉花的研究也主要集中在常規(guī)灌溉方式上，對膜下滴灌棉花水分生產(chǎn)函數(shù)的研究較少。本文試圖通過對膜下滴灌棉花水分生產(chǎn)函數(shù)的分析，運用邊際均衡原理分析灌溉用水量與棉花產(chǎn)量之間的最佳結合點，在獲得較高產(chǎn)量的同時，又有較高的灌水效率，能使水資源投入的經(jīng)濟效益達到最佳。為優(yōu)化膜下滴灌灌溉制度和經(jīng)濟用水提供理論依據(jù)和實踐指導。

1　材料與方法

1.1研究區(qū)概況

試驗于2005—2012年在中國科學院阿克蘇農(nóng)田生態(tài)系統(tǒng)國家野外科學觀測研究站(40°37′N, 80°45′E)進行。該地氣候屬于暖溫帶干旱氣候，年平均氣溫11.2℃，與同緯度地區(qū)相比，夏季溫度偏高、冬季溫度偏低、春秋季節(jié)氣溫變化劇烈。多年平均降水量45.7 mm，年水面蒸發(fā)量2 500 mm，無霜期207 d，全年日照時數(shù)2 940 h，年太陽輻射總量6 000 MJ·m-2。2008年降水量為88 mm，其中棉花生育期4—11月份降水量為80.9 mm。試驗地土壤為粉砂壤土。

1.2試驗設計

試驗采用隨機區(qū)組設計，灌溉量共設6個水平，分別為：4 265 m3·hm-2(Ⅰ)、3 926 m3·hm-2(Ⅱ)、3 600 m3·hm-2(Ⅲ)、3 271 m3·hm-2(Ⅳ)、2 947 m3·hm-2(Ⅴ)、2 618 m3·hm-2(Ⅵ)，每處理4次重復，小區(qū)面積389.2 m2。供試棉花品種為中棉49，株行距采用寬窄行配置(10 cm+65 cm+10 cm+60 cm)，即一膜種植四行棉花，窄行距10 cm，寬行距65 cm，膜間距60 cm，株距均為10 cm。4月28日播種，播種前施磷酸二銨450 kg·hm-2作基肥。全生育期灌溉采用滴灌模式，滴灌帶鋪設在寬行中央，每7 d滴灌一次，每次灌溉持續(xù)7～11 h，6月24日(蕾期)第1次灌水，共灌水10次。在第3～6次灌溉時，隨水追施尿素450 kg·hm-2、磷酸二氫鉀150 kg·hm-2。

1.3作物水分生產(chǎn)函數(shù)理論

1.3.1膜下滴灌棉花滴灌量與產(chǎn)量的關系在灌溉制度一定的情況下，最佳滴灌量與產(chǎn)量的關系，可通過水分生產(chǎn)函數(shù)的分析來確定。研究作物水分的消耗和作物產(chǎn)量之間的關系函數(shù)，稱水分生產(chǎn)函數(shù)。為了研究灌溉對作物產(chǎn)量的影響，假定在整個生產(chǎn)過程中，只有水資源一種投入物的數(shù)量是不斷變化的，且施用技術既定，其他資源如肥料等均被視為固定資源，氣候、光照等不可控資源或供應量無法限制的資源也被對照處理所固定。

大量研究表明[17,21-26]，隨著水源條件的改善和管理水平的提高，產(chǎn)量與耗水量的關系中出現(xiàn)了一個明顯的界限值，當耗水量小于此界限值時，產(chǎn)量隨耗水量的增加而增加，開始增加的幅度較大，然后減少；當達到該界限值時，產(chǎn)量不再增加，其后Y隨ET增加而減少。因此，作物產(chǎn)量與全生育期總耗水量的函數(shù)關系用二次拋物線關系表達，即：

Y=a×ET2+b×ET+c

(1)

式中，Y為棉花的產(chǎn)量(kg·hm-2)，ET為棉花耗水量(mm)，a、b、c為經(jīng)驗系數(shù)，通過試驗資料確定。

通過膜下滴灌灌水試驗資料的統(tǒng)計分析得出：拋物線關系最符合實驗區(qū)的實際情況。用最小二乘法原理擬合的水分生產(chǎn)函數(shù)模型y=ax2+bx+c能夠較好地反映膜下滴灌棉花的連續(xù)性及生長規(guī)律，且參數(shù)a、b、c易于估計。因此，采用該模型比較合理，且擬合程度較好，精度較高，易于應用推廣。

1.3.2水分利用效率水分利用效率(Water use efficiency, WUE)即水資源的平均生產(chǎn)能力，表示一單位水資源生產(chǎn)的作物產(chǎn)量或生產(chǎn)1 kg作物產(chǎn)量的耗水量，它是節(jié)水農(nóng)業(yè)研究的最終目標[17]。作物水分利用效率表達方式有三種，一是用作物耗水量(ET)，這是普遍所指的水分利用效率；二是用灌溉水量(I)，得到的是灌溉水利用效率，它對于確定最佳灌溉定額必不可少，在節(jié)水灌溉中意義重大；三是用天然降水(P)，得到的是降水利用效率，它是旱地節(jié)水農(nóng)業(yè)中的重要指標。

膜下滴灌棉花水分利用效率可以表示為：

WUE=Y/ET

(2)

式中，WUE為棉花水分利用效率，Y為棉花產(chǎn)量，ET為棉花耗水量。

1.3.3棉花灌水邊際效益灌水的邊際效益，是指在其他情況不變的條件下，增加一單位水，給生產(chǎn)帶來產(chǎn)值的增加量，邊際值Mc是邊際曲線函數(shù)的導數(shù)。用公式來表示，則有：

Mc=(ΔY/ΔI)×PY=(dY/dI)×PY

(3)

式中，Mc為邊際效益，ΔY為產(chǎn)量的變化量，ΔI為灌水的變化量，PY為棉花單價。隨著用水量的增加，水的邊際效益是在遞減的，邊際效益曲線是一條向右下方傾斜的曲線。

灌水的邊際成本，是指灌水單價PI，當灌水量增加，邊際效益下降但邊際成本不變。

灌水的邊際利潤，是指增加單位灌水成本所增加的灌水利潤，用R表示。令灌水利潤為B，灌水成本為C，則有:

B=ΔY×PY-C

(4)

dB/dC=d(ΔY×PY-C)/dC=d(ΔY×PY)/dC-1=PY(dΔY/dI)×(dI/dC)-1

(5)

因為：dC/dI=PI

因此式(5)可轉化為：

dB/dC=PY×(dY/dI)×(1/PI)-1=(dY/dI)×(PY/PI)-1

(6)

由上式可得出：

dY/dI=(dB/dC+1)×(PI/PY)

(7)

因為：邊際利潤R=(dB/dC)

則有：

dY/dI=(R+1)×(PI/PY)

(8)

由上式可見：

當R>0時，邊際效益(dY/dI)×PY大于邊際成本PI，只要繼續(xù)增加灌水投入量，灌水的凈效益仍可繼續(xù)增加，R值越大，邊際效益越高。

當R<0時，邊際效益(dY/dI)×PY小于邊際成本PI，此時增加灌水投入量所取得的收益已經(jīng)不能抵償成本。當R=-1時，dY/dI=0，產(chǎn)量達到最高，此時的灌水量即為最高產(chǎn)量灌水量。

當R=0時，邊際效益(dY/dI)×PY等于邊際成本PI，收投抵消，凈收益最大，此時的灌水量即為最佳效益灌水量。

1.4研究方法

1.4.1土壤水分觀測在不同水分處理棉花寬行、窄行、膜間分別埋設中子管，用中子水分儀(CNC503DR)測定0～20、20～40、40～60、60～80、80～100 cm土層的土壤含水量。每隔2 d測定一次，滴灌前后加測。為了解決中子儀測定表層土壤水分不準確的問題，利用烘干法對中子儀進行標定。降雨量取自阿克蘇農(nóng)田生態(tài)系統(tǒng)國家野外科學觀測研究站。

1.4.2膜下滴灌棉花耗水量確定方法根據(jù)農(nóng)田水量平衡原理，膜下滴灌棉田水量平衡式為：

ET=I+P-(W2-W1)-Q

(9)

式中，ET為棉花耗水量(mm)；I為灌水量(mm)；P為有效降雨量(mm)；W1為時段初80cm土體儲水量(mm)；W2為時段末80cm土體儲水量(mm)；Q為根系層(80cm)下界面水分交換量(mm)，取向下為正。本研究中通過簡化非飽和土壤水運動達西定律獲得根區(qū)下界面水分通量，非飽和土壤水分運動達西定律如下式：

(10)

式中，θ為根區(qū)下界面處的土壤含水率(m3·m-3)，由實測獲得；D(θ)為非飽和土壤水分擴散率(mm2·d-1)；K(θ)為非飽和土壤導水率(mm·d-1)，利用已有的研究成果[27]；z為土壤水的高度(mm)。式(10)的中心差分方程為：

(11)

1.4.3數(shù)據(jù)處理數(shù)據(jù)采用單因子方差分析(One-way ANOVA)，如果ANOVA結果顯著，再用最小顯著差數(shù)法(LSD)對各處理進行兩兩配對多重比較。所有統(tǒng)計分析均采用Excel和 SPPS 13.0(SPSS公司，美國)完成。

2　結果與分析

2.1膜下滴灌棉花耗水量與耗水規(guī)律

由表1看出，隨滴灌量減小，棉花總耗水量減小。由于苗期土壤水分一致，各處理耗水量基本相同，從蕾期開始經(jīng)過灌水處理，耗水量隨滴灌量減小而減小。各處理棉花花鈴期耗水量最大，其次是蕾期，苗期和吐絮期較小。各處理棉花的日耗水強度，苗期差異不明顯，其值在1.2～1.3 mm·d-1；蕾期以后日耗水強度增大，并隨滴灌量的增加而增大。蕾期各處理在2.5～3.3 mm·d-1；進入花鈴期后，日耗水強度達到最大值，其值在3.7～5.5 mm·d-1變化；吐絮期又降低，在1.4～2.2 mm·d-1之間?？梢?，膜下滴灌棉花花鈴期是其生育期中耗水強度最大的階段，耗水量占全生育期耗水總量的50%以上，其次是蕾期，苗期和吐絮期最小。

表1　不同水分處理棉花各生育階段的耗水量(2008年)

注：不同小寫字母表示處理間差異達顯著水平(P<0.05)。下同。

Note: Different small letters meant significant difference among different water treatments at 0.05 level. The same below.

2.2膜下滴灌棉花產(chǎn)量與水分利用效率

2.2.1棉花產(chǎn)量、水分利用效率、邊際產(chǎn)量與耗水量的關系通過對試驗數(shù)據(jù)進行統(tǒng)計分析、數(shù)值曲線擬合，求得塔里木灌區(qū)膜下滴灌棉花籽棉產(chǎn)量Y與耗水量ET的關系如下式，相關分析得方程極顯著。

Y=-0.1976ET2+212.11ET-50506

(R2=0.9648，n=6)

(12)

式中，Y為籽棉產(chǎn)量,kg·hm-2；ET為耗水量,mm。

對此函數(shù)求一階導數(shù)，可得產(chǎn)量最大為6 415 kg·hm-2時的耗水量為537 mm和水分利用效率為11.946 kg·hm-2·mm-1。從圖1Y-ET曲線可看出，起初隨耗水量的增加棉花產(chǎn)量也逐漸增加，當耗水量達537 mm時，棉花產(chǎn)量最大(6 415 kg·hm-2)，此時灌水邊際效益趨于0，水資源在棉花生長發(fā)育過程中的短邊效應已不存在。如果沒有生產(chǎn)條件的重新組合，水資源生產(chǎn)力將表現(xiàn)為負增長，即隨灌水量增加，棉花會減產(chǎn)。

棉花水分利用效率是棉花產(chǎn)量與耗水量的比值，由式(12)可得：

WUE=212.11-0.1976ET-(50506/ET)

(13)

從圖1中WUE-ET曲線可看出，水分利用效率和耗水量的關系也近似為二次曲線。起初水分利用效率隨灌水量增加而增加，當耗水量達到506 mm時，水分利用效率最大，為12.310 kg·hm-2·mm-1；當耗水量超過506 mm后，水分利用效率隨著耗水量的增加而減小。顯然，從提高水分利用效率的角度考慮，棉花總耗水量以達到506 mm為宜。

圖1膜下滴灌棉花產(chǎn)量(Y)、邊際產(chǎn)量(M)、水分利用效率(WUE)與耗水量(ET)的關系

Fig.1Relationships between water consumption(ET) and yield (Y) with drip irrigation under mulching, marginal yield(M), and water use efficiency (WUE) of cotton

邊際產(chǎn)量M是指每增加1 mm耗水量所能增加的產(chǎn)量，說明了當水分投入增加時所引起的產(chǎn)量變動率，它是Y-ET關系的一階導數(shù)：

M=dY/dET=212.12-0.3952ET

(14)

隨著水資源投入量增加，生產(chǎn)力由遞增到遞減表現(xiàn)了生產(chǎn)函數(shù)的一般特性，這個過程體現(xiàn)了水資源生產(chǎn)力的三個不同發(fā)展階段。Y-ET曲線反映的水資源生產(chǎn)力發(fā)展的三個階段，反映了小麥產(chǎn)量與水分關系的一般特征。

第一階段，從原點至A點之間，生產(chǎn)力迅速增長階段。邊際產(chǎn)量高于水分利用效率，水資源總生產(chǎn)力迅速增長，屬于水分極缺區(qū)?？偖a(chǎn)量由最低點上升到水分利用效率最高點，水分利用效率最高為12.310 kg·hm-2·mm-1，此時的耗水量為506 mm，對應產(chǎn)量為6 229 kg·hm-2。在這個階段，平均每投入一單位水資源所取得的產(chǎn)品量達最大，單位資源投入的成本最小。說明水資源利用尚存在較大的潛力。在這一范圍內，在水量條件允許時，應盡可能增加投入，以便獲得最大產(chǎn)量。

第二階段，A點至B點之間，生產(chǎn)力緩慢增長階段。水分利用效率逐漸降低，屬于水分缺乏調節(jié)區(qū)。隨著水資源投入量的增加，資源總生產(chǎn)力繼續(xù)增長?？偤乃坑?06 mm增加至537 mm，棉花產(chǎn)量由6 229 kg·hm-2上升到6 415 kg·hm-2，進入水資源投入的高產(chǎn)區(qū)，但資源投入的平均生產(chǎn)力或資源利用效率卻在下降。從經(jīng)濟意義上看，這一階段資源得到相對充分利用，資源效益得到較充分發(fā)揮。最佳耗水量就在該段范圍內，必須嚴格控制和調節(jié)耗水量。

第三階段，最高產(chǎn)量點B以后，生產(chǎn)力回落階段。在這一階段隨著水資源用量的增加，總產(chǎn)量不增加反而下降，邊際產(chǎn)量為負值，屬于奢侈耗水區(qū)。在該段越增加水分投入，減產(chǎn)越嚴重，造成的虧損越大，屬于不合理生產(chǎn)階段。

生產(chǎn)函數(shù)三階段的劃分，說明棉花產(chǎn)量與水資源利用效率之間存在著差異，二者并非同步變化。產(chǎn)量的最高點與WUE的最高點不一致，當水資源利用效率達到最高時，棉花產(chǎn)量遠未達到最大，相應于WUE最高點的ET要低于Y最高點的ET。例如在Y-ET曲線的頂點B，棉花產(chǎn)量達最大(6 415 kg·hm-2)，其水分利用效率為11.946 kg·hm-2·mm-1。水分利用效率最高時棉花產(chǎn)量6 229 kg·hm-2，與最高產(chǎn)量相差186 kg·hm-2，說明這時水資源潛力尚未充分發(fā)揮。在生產(chǎn)中，若單純追求資源利用效率指標，無異于以資源投入的低效益換取較高的資源利用效率，必然造成水資源的隱形浪費。而若單純追求產(chǎn)量而忽視資源的利用效率，同樣會造成資源浪費和利用不經(jīng)濟。

如果以獲得單產(chǎn)最高為灌溉原則，應按B點確定灌溉水量。但在水資源有限的地區(qū)，是以有限的水使較小的耕地面積獲得最高產(chǎn)量的經(jīng)濟效益大，還是適當減少單位面積的灌水量，以獲得最高產(chǎn)量而換取較多的灌溉面積所取得的經(jīng)濟效益大，是灌溉決策者必須認真分析的一個關鍵問題。

圖2膜下滴灌棉花產(chǎn)量(Y)、平均產(chǎn)量(A)、邊際產(chǎn)量(M)與灌水量(I)的關系

Fig.2Relationships between irrigation amounts(I) and yield(Y), average yield(A), marginal yield(M) of cotton with drip irrigation under mulching

2.2.2棉花產(chǎn)量、邊際產(chǎn)量與滴灌量的關系膜下滴灌棉花籽棉產(chǎn)量(Y)與全生育期滴灌量(I)的關系，經(jīng)回歸分析，符合二次拋物線關系：

Y=-0.0026 I2+18.015 I-24845(R2=0.9592，n=6)

(15)

為了分析灌溉效應，引入灌溉平均產(chǎn)量(A)和灌溉邊際產(chǎn)量(M)兩個基本概念，且定義灌溉平均產(chǎn)量(A)是指平均每投入單位灌水量所獲得的產(chǎn)量[28]，即灌溉水分利用效率；灌溉邊際產(chǎn)量(M)是指增加單位灌水量所增加的產(chǎn)量，它是Y-I關系的一階導數(shù)，A、M的單位為kg·hm-2·m-3。

M=dY/dI=18.015-0.0052I

(16)

A=Y/I=18.015-0.0026I-(24845/I)

(17)

由圖2可見，三條曲線有兩個重要的轉折點把Y-I曲線分為三段：第一點是M-I曲線與A-I曲線的交點C，是平均產(chǎn)量曲線的最高點，該點邊際產(chǎn)量等于平均產(chǎn)量，灌溉水分利用效率最大，滴灌量為3 091 m3·hm-2；第二點是Y-I曲線的最高點D，該點邊際產(chǎn)量M等于零，產(chǎn)量最大，滴灌量為3 464 m3·hm-2，過此點邊際產(chǎn)量為負。曲線第一段是從原點到C點，邊際產(chǎn)量大于平均產(chǎn)量，產(chǎn)量增加的幅度大于灌水量增加的幅度。此段屬于棉花水分極缺區(qū)，灌溉量的增加會顯著提高產(chǎn)量，因此只要水量允許，就應予以投入；第二段是從C點到D點，雖然邊際產(chǎn)量大于零，但已小于平均產(chǎn)量，產(chǎn)量增加的幅度小于灌水量增加的幅度，屬于棉花水分缺乏調節(jié)區(qū)。第三段是D點以后，邊際產(chǎn)量小于零，產(chǎn)量隨灌水量增大而減小，屬于奢侈灌溉區(qū)?？梢?，在第一段內，灌溉的增產(chǎn)效應沒有充分發(fā)揮出來，在第三段增加灌水量反而導致減產(chǎn)，從節(jié)水和經(jīng)濟上看，在第一段和第三段均屬于不合理灌溉階段，第二段屬于合理灌溉階段。

2.3膜下滴灌棉花灌水邊際效益分析2.3.1最高產(chǎn)量灌溉定額對于膜下滴灌棉花灌溉生產(chǎn)函數(shù)而言，隨著滴灌量的增加，棉花產(chǎn)量產(chǎn)值相應增加，當?shù)喂嗔吭黾拥揭欢ㄖ禃r，棉花產(chǎn)量隨滴灌量增加的速度變慢，當產(chǎn)量達到最大時，隨滴灌量的增加，產(chǎn)量反而下降。邊際效益由大到小以至于出現(xiàn)負值。從式(8)可知，當R=-1時，dY/dI=0，此時產(chǎn)量最大(D點)，灌水投入方案為最優(yōu)方案，最高產(chǎn)量為6 360.8 kg·hm-2，相應的滴灌量為3 464 m3·hm-2，灌溉水分利用效率為1.836 kg·hm-2·m-3。

2.3.2最佳效益灌溉定額在干旱半干旱區(qū)，水資源短缺是制約農(nóng)業(yè)發(fā)展的關鍵因素，因此，dY/dI=0雖為灌溉最優(yōu)方案，但不是最經(jīng)濟方案，當水量投入所增加的產(chǎn)出大于或等于水量投入本身的價值時，即邊際收益等于邊際成本時的投入方案是最經(jīng)濟方案。設灌溉單方水的價格為PI，籽棉價格PY，從式(8)可知，當R=0時，單位面積上的灌溉利潤或收益最大，即：

dY/dI=PI/PY

(18)

由式(16)、式(18)可得到最佳效益灌溉定額為：

I=(PI/PY-18.015)/(-0.0052)

(19)

不同籽棉價格和單方水價格下的最佳效益灌溉定額見表2。

表2　不同籽棉價格和單方水價格下的最佳效益灌溉定額/(m3·hm-2)

由式(19)、表2可以看出，最佳效益灌溉定額不僅與棉花水分生產(chǎn)函數(shù)有關，而且與籽棉價格和單方水價格有關，當單方水價一定時，最佳效益灌溉定額隨籽棉價格的增加有增加趨勢，尤其是在單方水價相對較高時。當籽棉價格一定時，最佳效益灌溉定額隨單方水價的增加有減小趨勢，但變化趨勢不是很明顯。由于灌溉用水水價基本恒定，因此不管灌水量為多少，用水的增加和成本的增加都是呈等比例變化的，邊際成本不會隨著用水量的增加而改變。當水價為0.20 元·m-3，籽棉價格為7 元·kg-1，生產(chǎn)要素不變費用為19 500 元·hm-2時，邊際收益等于邊際成本時的灌水量，即最佳效益灌溉定額為3 459 m3·hm-2，此時產(chǎn)量為6 360.7 kg·hm-2，較最高產(chǎn)量時的滴灌量3 464 m3·hm-2，節(jié)水5 m3·hm-2。此時凈收益達到最大，為24 333.1 元·hm-2，每立方灌水量凈收益為7.23 元·hm-2。與最高產(chǎn)量點相比，每立方灌水產(chǎn)量由1.836 kg·hm-2增加到1.839 kg·hm-2，水分利用效率提高0.003 kg·m-3，每立方灌水凈收益增加0.21 元·hm-2，棉花總產(chǎn)量卻基本一致?？梢娮罴研б纥c灌水凈收益和灌水利用效率較最高產(chǎn)量點均有所增長，但僅比最高產(chǎn)量點節(jié)水5 m3·hm-2，且產(chǎn)量基本一致，凈收益只增加0.3 元·hm-2，差異不大。這主要是由于膜下滴灌提高了灌溉水利用效率，同時目前單方水價格較低，低水價對最佳效益灌溉定額影響較小。

2.3.3高效用水灌溉定額棉花高效用水追求的目標是單位灌水量所生產(chǎn)出的籽棉產(chǎn)量最大，即灌溉水分利用效率最大。由式(17)可得：

dY/dI=0.0026-(24845/I2)

(20)

當dY/dI=0時，灌溉水分利用效率最大。代入上式得棉花高效用水灌溉定額為3 091 m3·hm-2，相應的最大灌溉水分利用效率為1.941 kg·hm-2·m-3，此時棉花產(chǎn)量為5 998 kg·hm-2。

3　結論與討論

通過塔里木灌區(qū)膜下滴灌棉花灌水實驗，得到以下結論:

1) 塔里木灌區(qū)膜下滴灌棉花花鈴期是其生育期中耗水強度最大的階段，耗水量占全生育期耗水總量的50%之多，其次是蕾期，苗期和吐絮期最小。棉花的需水臨界期是花鈴期和蕾期，有限的水資源應首先分配在這兩個生育階段。

2) 塔里木灌區(qū)膜下滴灌棉花產(chǎn)量與耗水量呈二次拋物線關系，遵循“報酬遞減規(guī)律”。水分利用效率最高為12.310 kg·hm-2·mm-1，此時的耗水量為506 mm，產(chǎn)量為6229 kg·hm-2。棉花產(chǎn)量最高為6 415 kg·hm-2時，耗水量為537 mm，水分利用效率為11.946 kg·hm-2·mm-1。根據(jù)水分利用效率的變化規(guī)律和產(chǎn)量與耗水量的關系，將水資源生產(chǎn)力劃分為三個階段。耗水量小于506 mm為第一階段，產(chǎn)量隨耗水量的增加而增加，屬于水分極缺區(qū)，在水量條件允許時，應盡可能增加投入，以便獲得最大產(chǎn)量。耗水量從506 mm到537 mm為第二階段，產(chǎn)量仍隨耗水量的增加而增加，但增加的幅度減小，生產(chǎn)力緩慢增長，水分利用效率逐漸降低，屬于水分缺乏調節(jié)區(qū)。耗水量超過537 mm以后為第三階段，產(chǎn)量隨耗水量的增加而逐漸減小，出現(xiàn)產(chǎn)量負增長，屬于奢侈耗水區(qū)，在該段越增加水分投入，減產(chǎn)越嚴重，造成的虧損越大，屬于不合理生產(chǎn)階段。

3) 塔里木灌區(qū)膜下滴灌棉花產(chǎn)量與滴灌量呈二次拋物線關系，產(chǎn)量的最高點與灌溉水利用效率的最高點并不一致，相應于灌溉水利用效率最高點的滴灌量要低于產(chǎn)量最高點的滴灌量。在現(xiàn)有灌溉技術和地下水埋深條件下，膜下滴灌棉花合理灌溉定額為3 091～3 464 m3·hm-2。最高產(chǎn)量為6 360.8 kg·hm-2，其灌溉定額為3 464 m3·hm-2，灌溉水分利用效率為1.836 kg·hm-2·m-3；高效用水灌溉定額為3 091 m3·hm-2，此時棉花產(chǎn)量為5 998 kg·hm-2，灌溉水分利用效率最大，為1.941 kg·hm-2·m-3。

4) 膜下滴灌棉花的最佳效益點并非水分利用效率的最高點和產(chǎn)量最高點，而是存在于3 091～3 464 m3·hm-2區(qū)間內，在產(chǎn)量最高點之前?，F(xiàn)狀條件下邊際效益等于邊際成本時，凈收益達到最大(24 333.1 元·hm-2)，每立方灌水量凈收益7.23 元·hm-2，最佳效益灌溉定額為3 459 m3·hm-2，其產(chǎn)量為6 360.7 kg·hm-2，較最高產(chǎn)量滴灌量節(jié)約用水5 m3·hm-2，每立方灌水凈收益增加0.21 元·hm-2，水分利用效率提高0.003 kg·m-3，而產(chǎn)量基本相同。最佳效益點雖然比最高產(chǎn)量點節(jié)約用水5 m3·hm-2，灌水凈收益和灌水利用效率均有增長，且產(chǎn)量基本一致，但兩點之間的變化非常微小，主要是由于膜下滴灌提高了灌溉水利用效率，同時目前單方水價格較低，低水價對最佳效益灌溉定額影響較小，這也符合報酬遞減規(guī)律。

[1]呂殿青,邵明安,王全九,等.膜下滴灌土壤鹽分特性及影響因素的初步研究[J].灌溉排水學報,2001,20(1):28-31.

[2]呂殿青,王文焰,王全九.滴灌條件下土壤水鹽運移特性的研究[J].灌溉排水,2000,19(1):16-21.

[3]李毅,王文焰,王全九,等.溫度勢梯度下土壤水平一維水鹽運動特征的實驗研究[J].農(nóng)業(yè)工程學報,2002,18(6):4-8.

[4]王新平,李新榮,康爾泗.干旱沙區(qū)滴灌條件下水鹽運移過程試驗研究[J].干旱地區(qū)農(nóng)業(yè)研究,2002,20(3):44-48.

[5]王全九,王文焰,汪志榮,等.鹽堿地膜下滴灌技術參數(shù)的確定[J].農(nóng)業(yè)工程學報,2001,17(2):47-50.

[6]呂殿青,王全九,王文焰,等.膜下滴灌水鹽運移影響因素研究[J].土壤學報,2002,39(6):794-801.

[7]張瓊,李光永,柴付軍.棉花膜下滴灌條件下灌水頻率對土壤水鹽分布和棉花生長的影響[J].水利學報,2004,35(9):123-126.

[8]趙成義,閆映宇,李菊艷,等.塔里木灌區(qū)膜下滴灌的棉田土壤水鹽分布特征[J].干旱區(qū)地理,2009,32(6):892-898.

[9]Raine S R, Meyer W S, Rassam D W, et al. Soil-water and solute movement under precision irrigation: Knowledge gaps for managing sustainable root zones[J]. Irrigation Science, 2007,26(1):91-100.

[10]Coelho F E, Or D. Root distribution and water uptake patterns of corn under surface and subsurface drip irrigation[J]. Plant and soil, 1999,206(2):123-136.

[11]閆映宇,趙成義,盛鈺,等.膜下滴灌對棉花根系、地上部分生物量及產(chǎn)量的影響[J].應用生態(tài)學報,2009,20(4):970-976.

[12]王興繁,胡順軍,田長彥,等.塔里木灌區(qū)棉田蒸散與棵間蒸發(fā)變化規(guī)律[J].干旱地區(qū)農(nóng)業(yè)研究，2012，30(6)：74-78.

[13]劉新永,田長彥,馬英杰,等.南疆膜下滴灌棉花耗水規(guī)律以及灌溉制度研究[J].干旱地區(qū)農(nóng)業(yè)研究,2006,24(1):108-112.

[14]雷成霞,何新林,王振華,等.新疆地下滴灌無膜移栽棉花水分生產(chǎn)函數(shù)試驗分析[J].灌溉排水學報,2011,30(4):132-134.

[15]金建華,孫書洪,王仰仁,等.棉花水分生產(chǎn)函數(shù)及灌溉制度研究[J].節(jié)水灌溉,2011,(2):46-48.[16]胡順軍,王仰仁,康紹忠,等.棉花水分生產(chǎn)函數(shù)Jensen模型敏感指數(shù)累積函數(shù)研究[J].沈陽農(nóng)業(yè)大學學報,2004,35(5-6):423-425.

[17]劉昌明,周長青,張士鋒,等.小麥水分生產(chǎn)函數(shù)及其效益的研究[J].地理研究,2005,24(1):1-10.

[18]王會肖,劉昌明.作物水分利用效率內涵及研究進展[J].水科學進展,2000,11(1) :99-104.

[19]Kipkorir EC, Raes D, Massawe B. Seasonal water production functions and yield response factors for maize and onion in Perkerra, Kenya[J]. Agricultural Water Management, 2002,56(3):229-240.

[20]Rajputgs, Singhj. Water production function for wheat under different environmental conditions[J]. Agricultural Water Management, 1986,11:318-332.

[21]馬忠明.有限灌溉條件下作物-水分關系的研究[J].干旱地區(qū)農(nóng)業(yè)研究,1998,16(2):75-79.

[22]雷志棟,胡和平,楊詩秀,等.葉爾羌河灌區(qū)冬小麥灌溉試驗與分析[J].灌溉排水,1999,18(2):30-33.

[23]王修貴,張祖蓮,趙長友,等.作物產(chǎn)量對水分虧缺敏感性指標的初步研究[J].灌溉排水,1998,17(2):25-30.

[24]陳亞新.作物-水模型及其敏感指標的確認[J].灌溉排水,1995,14(4):1-6.

[25]Kirda C. Deficit irrigation scheduling based on plant growth stages showing water stress tolerance[C]//Deficit Irrigation Practices. Water Reports 22, Rome:FAO, 2002.

[26]Wang Huixiao, Zhang Lu, DawesW R, et al. Improving water use efficiency of irrigated crops in the North China Plain measurements and modeling[J]. Agricultural Water Management, 2001,48:151-167.

[27]閆映宇,林新慧.灌溉方式對土壤水分運動參數(shù)的影響[J].水土保持研究,2011,18(4):160-165.

[28]胡順軍,周宏飛,宋郁東,等.渭干河灌區(qū)棉花高產(chǎn)節(jié)水優(yōu)化灌溉制度試驗研究[J].灌溉排水,2000,19(4):25-29.

Water production function and efficiency of cotton under mulched drip irrigation

YAN Ying-yu

(GeneralStationofWaterandSoilConversationandEco-environmentalMonitoringofXinjiang,Urumqi,Xinjiang830000,China)

This paper discusses the relationship between the amount of irrigation and the yield of cotton. Statistical analysis of experimental data on cotton under mulched drip irrigation during 2007—2012 indicates that the relationship between yield of cotton and irrigation presents a quadratic parabola in Tarim irrigated area. The irrigation norms for the highest yield, for the highest water use efficiency, and for the highest benefits, and the critical period of water requirement in cotton were analyzed by mathematical method. Best fit model of cotton water production function was adopted by least squares theory on the basis of field experimental data of cotton. Water production function is a quadratic function and its math expression isY=0.0026I2+18.015I-24845 (R2=0.9592，n=6). Water productivity is divided into three different development stages according to the changing rule of water use efficiency and the relationship between yield and irrigation. Irrigation amounts less than 3 091 m3·hm-2are identified as the first stage when yield rises with the increase of irrigation amounts, and productivity rises gradually at this stage. Irrigation amounts at the second stage rises from 3 091 m3·hm-2to 3 464 m3·hm-2and the yield still rises with the increase of irrigation amounts, but the range is reduced. Irrigation amounts at the third stage exceed 3 464 m3·hm-2and the yield is reduced with the increase of irrigation amounts. Negative rise occurs and productivity is declined. The analysis of cotton irrigation production function and the water requirement law indicates that the critical periods of water application for cotton are at the flowering stage and budding stage. The rational irrigation norm, norm for the highest yield, and norm for the highest water use efficiency for cotton are 3 091 to 3 464 m3·hm-2, 3 464 m3·hm-2, and 3 091 m3·hm-2, respectively. The best benefits point of irrigation efficiency of cotton is not the highest point of utilization ratio of water or the highest point of yield, but within the 3 091～3 464 m3·hm-2interval. When marginal benefit equals marginal cost, net income of irrigation on cotton is the biggest. The maximal net income and net income of unit of irrigation amounts are 24 333.1 yuan·hm-2and 7.23 yuan·hm-2, respectively, and irrigation amounts is 3 459 m3·hm-2at this time. Compared with the irrigation norm for the highest yield, it saves 5 m3·hm-2of water, meanwhile water use efficiency and net income of unit of irrigation amounts are increased by 0.003 kg·m-3, 0.21 yuan·hm-2, respectively, whereas the yield is almost the same.

drip irrigation under mulching; cotton; water production function; water use efficiency; irrigation efficiency

1000-7601(2016)04-0018-07

10.7606/j.issn.1000-7601.2016.04.04

2015-11-23

973課題“綠洲化的水、土、氣過程及其相互作用機制”(2009CB421302)；中國科學院創(chuàng)新項目(KZCX2-YW-127)

S562；S274.1

A