肥密組合對寒地半干旱區膜下滴灌玉米產量的影響

張發明，楊克軍，李佐同，趙長江，王玉鳳，徐晶宇，張翼飛，王智慧，付　健，

谷英楠，楊系玲，劉慧迪，劉　承，劉，方永江，唐春雙，王　聰

(寒地作物種質改良與栽培重點實驗室/黑龍江八一農墾大學農學院， 黑龍江 大慶 163319)

?

肥密組合對寒地半干旱區膜下滴灌玉米產量的影響

張發明，楊克軍，李佐同，趙長江，王玉鳳，徐晶宇，張翼飛，王智慧，付健，

(寒地作物種質改良與栽培重點實驗室/黑龍江八一農墾大學農學院， 黑龍江 大慶 163319)

摘要：為揭示黑龍江西部半干旱區膜下滴灌栽培方式下玉米高產措施，以鄭單958為材料，選取氮肥、磷肥、鉀肥和密度四因素作為試驗因素，每個試驗因素設計5個水平，按照二次通用旋轉組合設計(1/2實施)統計分析方法建立回歸模型，分析不同肥料和密度組合對玉米產量的影響。結果表明：四因素各單因子對籽粒產量有較明顯的影響，均呈開口向下的拋物線趨勢變化，影響順序為施氮>密度>施鉀>施磷，且單位水平施入量引起邊際產量的減少量為施磷>密度>施鉀>施氮；氮肥與磷肥之間的配合對產量的增加具有相互促進作用；要獲得≥12231.97 kg·hm(-2)的產量，氮肥、磷肥、鉀肥及種植密度的最優組合取值范圍為：氮肥246.24～279.59 kg·hm(-2)，磷肥133.77～151.52 kg·hm(-2)，鉀肥82.85～97.16 kg·hm(-2)，種植密度79 676～85 324 株·hm(-2)。

關鍵詞：玉米；膜下滴灌；施肥種類；種植密度；產量；寒地半干旱區

中國是世界第二大玉米生產國，產量約占全球玉米產量的1/5左右。玉米是中國種植面積最大的糧食作物，作為東北的第一大作物，占總播種面積的42.9%。目前，黑龍江省的玉米播種面積、總產量和商品化率均躍居全國第一，成為全國玉米第一大省，因其得天獨厚的自然條件而被世界公認為“黃金玉米帶”[1-4]。然而，干旱缺水、土壤肥力低下是限制黑龍江省西部半干旱地區玉米產量提高的兩個主要因素[5]。作物高產的形成囿于多方面的因素，提高種植密度是玉米高產的重要措施之一[6]。密植必然導致作物對礦質養分需求的增加，而肥料投入運籌也是實現高產和提高效益最有效的方法之一[7]。膜下滴灌技術具有提高土壤耕層溫度，并使作物近根區保持適宜濕度和養分濃度，同時具有節水、節肥等優點[8]。在國外，膜下滴灌多應用于有較高經濟價值園藝作物。美國水管理研究所對番茄、棉花和甜玉米等作物15年地下滴灌的研究結果表明，利用高頻度的滴灌可以減少深層滲漏量和提高水分利用率，并顯著提高作物產量[9]。滴灌番茄的產量和水分利用率均比溝灌高，株高和果實飽滿度也有相同的趨勢[10]。前人的研究多見于膜下滴灌條件下，水肥耦合對作物產量效應[11-15]及水肥施用最優配比研究，而對肥密這兩個因子交互作用增產效應的研究相對較少。因此，本試驗采用二次通用旋轉組合設計和膜下滴灌的方法，研究不同肥密組合對玉米產量的影響，在建立該類型地區玉米產量數學模型的基礎上，尋求高產高效優化方案，以期為玉米節水節肥高產高效栽培和產業發展提供理論和實踐的指導。

1材料與方法

1.1試驗地概況及品種

試驗于2014年5月至10月在黑龍江八一農墾大學農學院試驗基地(46°62′N，125°19′E)進行，平均海拔146 m，該區屬于典型的北溫帶亞干旱季風氣候區，全年降水較少，平均氣溫在5℃左右，年平均無霜期在143 d左右，2014年試驗區氣象數據如表1所示。供試土壤為草甸土。試驗田0～20 cm耕作層有機質含量26.62 g·kg-1,堿解氮130.42 mg·kg-1,速效磷55.99 mg·kg-1,速效鉀31.37 mg·kg-1，pH值8.13，試驗田地力均勻，地勢平坦，田間持水量為28.5%。試驗選用當地高產且主栽玉米雜交種鄭單958。

表1　2014年試驗區氣象數據

1.2試驗設計

1.3數據分析

采用DPS軟件進行統計分析；用Excel 2003和 Sigmaplot 10.0進行數據整理和作圖。

2結果與分析

2.1產量模型的建立

表2　試驗組合設計及試驗方案

表3　各時期土壤持水量指標

20種處理下玉米的產量如表4。利用產量結果，以二次通用旋轉組合設計進行回歸模擬，得到產量(Y)與氮(X1)、磷(X2) 鉀(X3)、密度(X4)4因子的回歸模型：

(1)

2.2模型分析

2.2.1試驗因子的產量效應分析主因子效應分析：試驗設計中各因素均經無量綱線性編碼處理，且各一次項系數與交互項、平方項的回歸系數間都是不相關的，所得偏回歸系數已標準化。故其回歸系數絕對值的大小可直接反映變量X對產量Y的影響程度。分析產量模型，其主效應表現為：

企業是社會經濟的主體，必須改變企業的經營模式，在企業內部實施綠色生產，追求經濟利益的同時，重視環境及社會效益。一方面，企業要在生產的各個環節減少污染物的排放，重視對環境及生態的保護；另一方面，生產過程中應該不斷地提高資源利用率，降低對各種資源的消耗。

1) 一次項X1、X2、X4的系數均為正值，X3的一次項系數是負值，說明在試驗設計范圍內，氮、鉀和密度單因子都有增產效應，且由系數絕對值大小判斷出兩因子對產量影響順序為施氮量>密度>施鉀量；鉀肥單因子對提高籽粒產量有負作用，即隨著施鉀肥增加，玉米籽粒產量有降低趨勢；

2) 交互項X1X2系數為正值，其余交互項系數為負值，說明施氮與施磷之間的配合對產量的增加具有相互促進作用；

3) 二次項系數均為負值，說明產量隨各因素增加均呈開口向下的拋物線趨勢變化，即在最佳水平以下時，產量隨該因素的增大而增加，當產量水平超過最佳水平時，產量出現下降趨勢。

2.2.2單因子效應分析在玉米籽粒產量的回歸模型中，通過降維分析，分析各因素對產量的影響。將其它三個因素固定在接近適宜的“0水平”，得到各因素的一元二次回歸效應模型如下：

(2)

(3)

(4)

(5)

在試驗設計的因素水平值1.682～-1.682范圍內，根據這些回歸子模型分別作出各因子對玉米籽粒產量影響的效應圖，如圖1所示。

從圖1可以看出，施氮、施磷、施鉀以及種植密度4個因素都能影響玉米籽粒產量，并且產量隨各因素水平的提高均呈開口向下的拋物線趨勢變化，存在產量最高點。符合報酬遞減定律。四因素都有明顯的增產效應，各拋物線的頂點就是各單因素的最高產量值，與其相對應的便是各因素的最適投入量。在本試驗中，氮肥的最佳投入量為1(碼值)，實際用量則為285 kg·hm-2，此時產量可達15 645.38 kg·hm-2；磷肥和鉀肥的最佳投入量為0(碼值)，實際用量則為130 kg·hm-2和90 kg·hm-2，此時產量都達15 634.15 kg·hm-2；最適種植密度為0(碼值)即種植密度為82 500 株·hm-2，此時產量可達15 634.15 kg·hm-2。到達最適投入量時，產量最高；投入量繼續加大，產量則隨之減小。由圖1中還可以看出，在較低投入量時，密度的增產效果略高于其他措施。

圖1試驗因子的產量效應

Fig.1Effects of experimental factors on yield

2.2.3單因素邊際效應分析邊際產量可反映各因素的最適投入量和單位水平投入量變化對產量增減速率的影響，各因素在不同水平下的邊際產量可通過對回歸子模型(2)、(3)、(4)、(5)求一階偏導，則分別得到氮、磷、鉀和密度的各因素的邊際效應方程。

施氮：Y=996.03-1969.6X1

(6)

施磷：Y=74.22-3855.92X2

(7)

施鉀：Y=-177.11-2953.14X3

(8)

種植密度：Y=126.47-1186.08X4

(9)

水肥單因子效應如圖2所示：當另一因素取編碼值為零水平時，隨著氮、磷、鉀和密度投入量的增加，單位氮、磷、鉀和密度投入量的增產作用下降，說明四因素邊際效益均呈遞減趨勢；且邊際效益遞減率為施磷>施鉀>施氮>種植密度。各因子單位水平施入量引起邊際產量的減少量為施磷>密度>施鉀>施氮。

圖2單因子邊際效應

Fig.2Marginal effects of single factor

2.2.4因素間交互作用分析各因素對產量的影響不是孤立的，它們之間存在交互作用。根據回歸模型用Sigmaplot 10.0作出因素間的交互作用圖，以氮和磷，氮和密度以及氮和鉀的交互作用回歸效應子模型作圖(圖3)。曲面圖上各點的高度代表兩因子一定配比水平時的玉米的產量，曲面的高度越高，說明玉米的產量越高。從圖中還可以看出，當一個因子固定在某一水平時，玉米產量隨另一因子水平變化的規律。

從圖3(A)可見，當施氮量一定時，施磷量在-1.682～0區間水平的范圍內，玉米產量隨著施磷量的增加而增加；當施磷量在0～1.682區間水平的高施磷量范圍內，玉米產量隨施磷量的增加而降低，這說明在施氮量一定的情況下，施磷量太多或者過少，氮肥效果都得不到最大程度發揮，產量都不能達到最大值。與此同時，在較低的施氮條件下，施磷的增產效應不明顯；在較高施氮水平下磷肥的增產效果顯著提升。但當施磷量水平較低時，隨著施氮量的大量增加可能會造成減產，這種效應符合報酬遞減函數。此時，如果加大施用氮肥，則肥料利用率降低。從圖3(A)也可以得出玉米產量的最高值并不產生在施氮量和施磷量最大時，施氮量的高產臨界值在0.5水平左右，施磷量在0水平左右。

圖3氮和磷的交互作用(A)、氮和密度的交互作用(B)及氮和鉀的交互作用(C)對玉米籽粒產量的影響

Fig.3The coupling effects of nitrogen and phosphorus (A), nitrogen and density (B), nitrogen and potassium (C) on corn grain yield

圖3(B)表明，在對玉米籽粒產量的影響方面，施氮量和種植密度具有明顯的交互作用，當固定施氮量在一個水平時，玉米產量隨著密度的增加呈現先增高后降低的趨勢。而當種植密度一定時，施氮量<0編碼水平時，玉米籽粒產量隨著施氮量的增加，產量明顯增加，并且種植密度中等水平(即0編碼水平時)，增產效應的更顯著。所以，玉米籽粒產量與施氮量和種植密度的協調配合有密切關系。

圖3(C)表明，施氮量和施鉀量表現出了顯著的協同增產效應。施氮量處在0編碼水平，施鉀量在0.5水平時，此時玉米籽粒的產量最高，為15 885.96 kg·hm-2；施氮量處在-1.682編碼水平，施鉀量在-1.682低水平時，此時玉米產量最低，為3 454.77 kg·hm-2。當固定施氮量在一個水平時，玉米產量隨著施鉀水平的提高呈先升高后降低的動態變化趨勢，而當施鉀水平固定時，且施氮量<0編碼水平時，玉米籽粒產量隨著施氮量的增加，增產效果顯著且增幅較大，但當施氮量>0編碼水平時，隨著施氮量的增加可能會出現籽粒減產。所以，氮肥和鉀肥的施入量過高或過低都會限制產量的提高，不僅降低了肥料利用率，而且發揮不出作物最大的增產潛力。

2.3肥密組合優化方案

在本試驗條件下獲得的最高產量對大面積生產應用來說并不一定能代表實際的最佳水平，為了取得四因素在生產中應用的可靠性，采用頻數法進一步解析，在-1.682～1.682約束區間，所得方案中，有72套方案玉米產量≥12 231.97 kg·hm-2。其優化組合的置信區間見表5。通過模擬尋優分析，大慶地區玉米要獲得≥12 231.97 kg·hm-2的產量，在膜下滴灌條件下，氮肥與磷肥、鉀肥、種植密度最優組合取值范圍為：氮肥246.24～279.59 kg·hm-2，磷肥133.77～151.52 kg·hm-2，鉀肥82.85～97.16 kg·hm-2，種植密度79 676～85 324 株·hm-2，玉米的經濟效益和生態效益達到最佳。

3結論與討論

不同處理單因子對寒地半干旱區玉米籽粒產量有較明顯的影響，且施氮的影響達到了顯著水平。施氮與施磷之間的配合，具有明顯的相互促進、協同增產作用。產量隨各因素增加均呈開口向下的拋物線趨勢變化，符合報酬遞減定律。大慶地區玉米要獲得≥12 231.97 kg·hm-2的產量，氮肥與磷肥、鉀肥、種植密度配合最優組合取值范圍為：氮肥246.24～279.59 kg·hm-2，磷肥133.77～151.52 kg·hm-2，鉀肥82.85～97.16 kg·hm-2，種植密度79 676～85 324 株·hm-2，玉米的經濟效益、社會效益、生態環境效益達到最佳。

表5　優化提取方案中Xi取值頻率分布表

近年來，膜下滴灌種植模式在政府、水利部門的扶持下，已經得到大面積推廣應用。膜下滴灌與傳統耕種模式相比，可增產78.5%，收益增加11 425.5元·hm-2[17]，同時可明顯增溫保墑、節水抑鹽、促進作物提早成熟、提質增效、節時省工。玉米生產中，高密度及高肥不一定是最高產的，合理的肥密組合能夠擴庫增源、提高肥料利用率，進而獲得高產。將肥密組合與膜下滴灌進行創新集成應用到生產實踐中，可減小當前肥密管理的盲目性，能進一步提高肥料的利用率、減輕污染、節本增收。所以，在膜下滴灌技術的基礎上，調整合理的肥密組合尤其對水資源匱乏、肥料利用率低的寒地半干旱地區可持續農業發展具有重要的現實意義。本試驗主要分析了2014年秋季玉米籽粒產量的數據，具有一定的局限性，尚存在一些不確定性，需要結合一些單因子試驗進一步深入研究。

參 考 文 獻：

[1]姚麗鳳.借鑒美國經驗發展中國的玉米產業[J].農業經濟,2008,(7):46-47,76.

[2]孫德全,李綏艷,林紅,等.黑龍江省玉米主要病害發生原因分析及抗病育種對策[J].作物雜志,2009,(2):90-93.

[3]楊國航,孫世賢,張春原,等.東北早熟春玉米區玉米生產現狀和發展趨勢[J].玉米科學,2007,(4):143-145.

[4]劉海燕.黑龍江省玉米持續增產的主要限制因素與對策[J].作物雜志,2014,(2):9-13.

[5]李楠楠,張忠學.黑龍江半干旱區玉米膜下滴灌水肥耦合效應試驗研究[J].中國農村水利水電,2010,(6):88-90,94.

[6]陳傳永,侯玉虹,趙明,等.密植對不同玉米品種產量性能的影響及其耐密性分析[J].作物學報,2010,(7):1153-1160.

[7]王國棟,曾勝和,陳云,等.新疆滴灌春玉米密植高產栽培施肥效應研究[J].農業現代化研究,2014,(3):376-380.

[8]孫文濤,孫占祥,王聰翔,等.滴灌施肥條件下玉米水肥耦合效應的研究[J].中國農業科學,2006,(3):563-568.

[9]Ayars J E，Phene C J，Hutmacher R B，et al. Subsurface drip imiagtion of row crops: areview of 15 years of research at the Wooter Mangament Research Laboratory[J]. Agricultural water Managerment，1999，42(1):1-27.

[10]Ahanners F， Tadees T. Effect of drip furrow irrigation and plant spacing on yield of tomato at dire Dawa，Ethiopia[J]. Agricultural Water Management.1998，35(3):201-207.

[11]張秋英, 劉曉冰,金劍,等.水肥耦合對玉米光合特性及產量的影響[J].玉米科學,2001,(2):64-67.

[12]馮鵬,王曉娜,王清酈,等.水肥耦合效應對玉米產量及青貯品質的影響[J].中國農業科學,2012,(2):376-384.

[13]胡安焱,董新光,魏光輝,等.滴灌條件下水肥耦合對干旱區紅棗產量的影響[J].灌溉排水學報,2010,(6):60-63.

[14]牛新湘,許詠梅,馬興旺,等.地面灌溉條件下高產棉田水肥耦合的產量效應[J].干旱地區農業研究,2009,(1):53-55,61.

[15]王榮蓮,于健,趙永來,等.滴灌施肥水肥耦合對溫室無土栽培水果黃瓜產量的影響[J].節水灌溉,2009,(3):15-17,22.

[16]姬景紅,李玉影,劉雙全,等.覆膜滴灌對玉米光合特性、物質積累及水分利用效率的影響[J].玉米科學,2015,(1):128-133.

[17]李慧萍,趙樂軍.玉米膜下滴灌高產增收效益分析[J].現代農業科技,2011,(1):74,82.

Effects of density and fertilization combination on yield of maize under mulched drip irrigation in semi-arid and cold areas

ZHANG Fa-ming, YANG Ke-jun, LI Zuo-tong, ZHAO Chang-jiang, WANG Yu-feng, XU Jing-yu,ZHANG Yi-fei, WANG Zhi-hui, FU Jian, GU Ying-nan, YANG Xi-ling, LIU Hui-di, LIU Cheng, LIU Yu,FANG Yong-jiang, TANG Chun-shuang, WANG Cong

(KeyLaboratoryofCropGermplasmimprovementandcultivationinColdArea/CollegeofAgronomy,HeilingjiangBayiAgriculturalUniversity,Daqing,Heilongjiang163319,China)

Abstract：To investigate the cultivation practices for high yield of maize under mulched drip irrigation in the semi-arid and cold areas in western Heilongjiang province, maize variety Zhengdan 958 was used as the study material. In addition, nitrogen, phosphate, potash and planting density were selected as the four experimental factors, each of which was tested at five levels. A quadratic regression model was thereby employed to establish in accordance with rotation design statistical analysis methods (1/2 implement) to analyze the effects of different fertilizers and planting densities on grain yield. The results showed that each single factor had obvious influence on grain yield, exhibiting a trend of downward opening parabola. The effects of these four factors were in an order of nitrogen (N) > density > potash (K) > phosphate (P). Also, the marginal yield reduction caused by unit volume level applied was in a order of phosphate (P) > density > potash (K) > nitrogen (N). Coordination between nitrogen and phosphate had a mutual promoting effect on the increase of grain yield. To achieve a yield higher than 12231.97 kg·hm(-2), the best combination of nitrogen and phosphate fertilizer, potash fertilizer, and planting density were as follows: 246.24～279.59 kg·hm(-2) nitrogen; 133.77～151.52 kg·hm(-2) phosphate; 82.85～97.16 kg·hm(-2) potash and 79 676～85 324 plants·hm(-2) planting density.

Keywords:maize; semi-arid and cold areas; mulched drip irrigation; fertilizer and density; yield

中圖分類號：S513

文獻標志碼：A

作者簡介：張發明(1988—)，男，黑龍江延壽人，在讀碩士，研究方向為玉米高產理論與技術。E-mail：15164568108@126.com。通信作者：楊克軍(1968—)，男，山東莒縣人，博士，教授，主要從事作物高產栽培與耕作方面的研究。 E-mail:byndykj@163.com。

基金項目：糧食豐產科技工程(2011BAD16B11-03)；國家科技計劃課題(2013BAD07B01-07)；黑龍江省農墾總局科技攻關項目(HNK11A-03-02)；大慶市科技計劃項目(SJH-2013-35)

收稿日期：2015-02-18

doi:10.7606/j.issn.1000-7601.2016.02.05

文章編號：1000-7601(2016)02-0032-06