播期和密度對白銀沿黃灌區冬小麥產量的影響

俞華林 楊繼忠 鄭彩霞 趙寶勰 師學豪

摘要：采用飽和D最優設計，在白銀沿黃灌區試驗分析了播期與密度對冬小麥產量的影響。結果表明，隨著播期的推遲，小麥產量降低;在一定范圍內，通過增加密度可以提高小麥成穗數，進而達到提高產量的目的。小麥產量在6 300 kg/hm2以上時，95%的置信區間播期為-0.855 9～0.054 7，密度為-0.275 4～0.672 8，相應的播期在9月17日至9月29日、密度為231.52～338.19 kg/hm2時增產效果為佳。

關鍵詞：冬小麥;播期;密度;產量;飽和D最優設計

中圖分類號：S512.1? ? ? ?文獻標志碼：A? ? ? ? 文章編號：1001-1463（2021）06-0018-06

doi：10.3969/j.issn.1001-1463.2021.06.006

Effect of Sowing Date and Density on Winter Wheat Yield in

Baiyin Irrigation Area Along Yellow River

YU Hualin， YANG Jizhong， ZHENG Caixia， ZHAO Baoxie， SHI Xuehao

（Baiyin Institdue of Agricultural Sciences， Baiyin Gansu 730900， China）

Abstract：The effects of sowing date and density on winter wheat yield were analyzed by saturated D optimal design in Baiyin Irrigation Area Along Yellow River. The results showed that wheat yield decreased with the delay of sowing date. In a certain range， the number of panicles can be increased by increasing the density， and then the yield can be increased. When wheat yield was above 6 300 kg/hm2， the sowing date of 95% confidence interval was -0.855 9～0.054 7， and the density was -0.275 4～0.672 8. The corresponding sowing date was September 17 solstice on September 29， and the density was 231.52～338.19 kg/hm2， which yield increase effect was best.

Key words：Winter wheat;Sowing date;Density;Yield;Saturated D-optimal design

小麥作為我國的重要糧食作物之一，對保證國家糧食安全戰略目標起著極為重要的作用[1 ]，隨著國家提出黃河流域生態保護與高質量發展的戰略，白銀沿黃灌區開始引種冬小麥，主要是由于冬小麥收獲后復種蔬菜可以提高經濟效益，且可增加冬春季節的覆蓋度，防止水土流失和沙塵暴危害，具有良好的生態效益。播期和密度是影響小麥群體性狀、農藝性狀、產量和品質的重要因素[2 ]。同期播種的不同品種其適宜種植密度不同[3 ]，但不同地區間播期和密度對產量構成因素的影響不同[4 - 5 ]，種植密度過高或過低均會影響小麥的產量和品質[6 ]。近年來隨著全球氣候變暖，不同地區的生態環境均有不同程度的改變，有關白銀沿黃灌區冬小麥適宜播種時間和種植密度尚不確定。我們采用兩因素二項飽和D最優試驗設計，研究不同播期和種植密度對冬小麥產量極其構成的影響，旨在為確定白銀沿黃灌區冬小麥適宜的播期和密度，加快冬小麥在白銀沿黃灌區的推廣及冬小麥高產優質栽培提供參考。

1? ?材料與方法

1.1? ?供試材料

指示冬小麥品種為輪選169，由白銀市農業科學研究所提供。

1.2? ?試驗區基本概況

試驗于2017 — 2018年指在靖遠縣劉川鎮來窯村進行。海拔1 763 m，前茬蘿卜，土質沙壤，耕層土壤含有機質19.84 g/kg、全氮0.89 g/kg、全磷1.37 g/kg、氧化鉀16.13 g/kg、堿解氮49.45 mg/kg、速效磷35.37 mg/kg、速效鉀74.14 mg/kg。

1.3? ?試驗方法

試驗采用兩因素二項飽和D最優試驗設計（表1）[7 ]。對各因素的水平進行編碼，把各因素有量綱的自然變量zj變成無量綱的規范變量xj，播期編碼值用x1表示，密度時間編碼值用x2表示，進行線性變換，其公式為：

編碼值xj=（zj-z0 j）/？駐j? ? ?故：zj=xj×？駐j+z0 j

其中z0 j=該因素的（上限水平+下限水平）/2

？駐j=該因素的（上限水平-z0 j）/1

根據以上已知數據，對播期和密度分別進行編碼，求出各自的編碼值xj與自然變量Zj的對應值。

根據水平編碼值，試驗設6個處理（表2），其中播期設9月15日、9月25日、10月2日、10月10日4個水平，密度設150、248、307、375 kg/hm2 4個水平。播種行距為15 cm，小區面積為20 m2（5 m×4 m），3次重復，周圍設1 m寬保護行。播前施磷酸二銨300 kg/hm2、普通過磷酸鈣600 kg/hm2、尿素300 kg/hm2，結合灌二水追施尿素425 kg/hm2。生育期灌水4次，中耕除草3次，防治蚜蟲1次。其他管理同當地大田。

1.4? ?測定項目及方法

每小區去掉邊行收獲計產。成熟期在田間采用1 m雙行測定有效穗數，每小區取60穗室內考種統計穗粒數，籽粒風干后用電子天平稱重取平均值，統計數據均為3個重復的平均值。

1.5? ?數據分析

試驗數據采用Excel軟件進行統計分析。

2? ?結果與分析

2.1? ?播期和密度對小麥產量及其構成要素影響的反應方程式建立與檢驗

以y為目標性狀，采用二次多項式建立播期與密度的數學模型。

y=b0+b1 x1+b2 x2+b3 x12+b4 x22+b5 x1x2

式中，x1為播種期（日/月），x2為播種密度（kg/hm2），bi（i=0、1 、2 ……5）為回歸系數，y為目標性狀[8 ]。

2.1.1? ? 播期和密度對小麥產量影響的反應方程式建立與檢驗? ? 由表3中的小麥產量可求得輪選169的產量與播期和密度之間的多元回歸方程：

對所得方程進行回歸式顯著性檢驗，F（回歸）=（Q回/f 回）/（Q誤/f誤），Q回=∑∑（- Y）2，Q 誤=∑∑（Yij-i）2，f 回=6-1=5，f 誤= 10，其中Yij為第i個處理第j個重復數值，Yi為第i個處理各重復的平均值，為回歸值，Y為總平均值。計算F值：F（回歸）=（Q回/f 回）/（Q誤/f 誤）=（1479002.5/5）/（45895/10）= 64.45>F0.01 （5， 10）=5.64，差異達極顯著水平，則方程成立，所建立方程可靠，可用于生產中產量及優化栽培方案選擇。

2.1.2? ? 播期和密度對成穗數影響的反應方程式建立與檢驗? ? 由表3中的小麥成穗數，可求得輪選169的成穗數與播期和密度之間的多元回歸方程為：

對所得方程進行回歸式顯著性檢驗，F（回歸）=（Q回/f回）/（Q誤/f誤），Q回= ∑∑（-Y）2，Q誤=∑∑（Yij-i）2，f回= 6-1=5，f誤=10，其中Yij為第i個處理第j個重復數值，Yi為第i個處理各重復的平均值，為回歸值，Y為總平均值。計算F值：F（回歸）=（Q回/ f回）/（Q誤/f誤）=（83696.5/5）/（834.5/10）= 200.59>F0.01 （5，10）=5.64，差異達極顯著水平，則方程成立。

2.1.3? ?播期和密度對穗粒數影響的反應方程式建立與檢驗? ?由表3中的小麥穗粒數，可求得輪選169的穗粒數與播期和密度之間的多元回歸方程為：

對所得方程進行回歸式顯著性檢驗，F（回歸）=（Q回/f回）/（Q誤/f誤），Q回=∑∑（- Y）2，Q誤=∑∑（Yij-Yi）2，f回= 6-1=5，f誤= 10，其中Yij為第i個處理第j個重復數值，Yi為第i個處理各重復的平均值，為回歸值，Y為總平均值。計算F 值：F（回歸）=（Q回/f 回）/（Q誤 / f 誤）=（39.985/5）/（1.9/10）= 42.09>F0.01（5，10）=5.64，差異達極顯著水平，則方程成立。

2.1.4? ?播期和密度對千粒重影響的反應方程式建立與檢驗? ? 由表3中小麥的千粒重，可求得輪選169的千粒重與播期和密度之間的多元回歸方程為：

對所得方程進行回歸式顯著性檢驗， F（回歸）=（Q回/f回）/（Q誤/f誤），Q回=∑∑（- Y）2，Q誤=∑∑（Yij-Yi）2，f回=6-1= 5，f誤=10，其中Yij為第i個處理第j個重復數值，Yi為第i個處理各重復的平均值，為回歸值，Y為總平均值。計算F值：F（回歸）=（Q回/f回）/（Q誤/f誤）=（5.28/5）/（1.66/10）= 6.36>F0.01（5，10）=5.64，差異達極顯著水平，則方程成立。

2.2? ?播期對小麥產量及產量構成因素的影響

2.2.1? ? 播期對小麥產量的影響? ? 在方程（1）中，將密度編碼設為零，則可得到關于播期與小麥產量的一元二次方程為：

從一元二次方程（5）可以得出，函數以x=-0.776為對稱軸，函數開口向下。當x1在[-1，-0.776]范圍內，函數y值是隨x1的增大而增大，即在此范圍內小麥產量隨播期推遲而上升;當x1在[-0.776，1]范圍內，函數y值是隨x1的增大而減小，即在此范圍內小麥產量隨播期的推遲而下降。

2.2.2? ? 播期對小麥成穗數的影響? ? 在方程（2）中，將密度編碼設為零，則可得到關于播期與小麥成穗數的一元二次方程為：

從一元二次方程（6）式可以得出，函數以x=-0.903為對稱軸，函數開口向下。當x1在[-1，-0.903]范圍內，函數y值是隨x1的增大而增大，即在此范圍內小麥成穗數隨播期推遲而上升;當x1在[-0.903，1]范圍內，函數y值是隨x1的增大而減小，即在此范圍內小麥成穗數隨播期的推遲而下降。

2.2.3? ? 播期對小麥穗粒數的影響? ? 在方程（3）中，將密度編碼設為零，則可得到關于播期與小麥穗粒數的一元二次方程為：

從一元二次方程（7）式可以得出，函數以x=3.979為對稱軸，函數開口向上。當x1在[-1，1]范圍內，函數y值是隨x1的增大而減小，即在此范圍內小麥穗粒數隨播期推遲而下降。

2.2.4? ? 播期對小麥千粒重的影響? ? 在方程（4）中，將密度編碼設為零，則可得到關于播期與小麥千粒重的一元二次方程為：

從一元二次方程（8）式可以得出，函數以x=-0.302為對稱軸，函數開口向下。當x1在[-1，-0.302]范圍內，函數y值是隨x1的增大而增大，即在此范圍內小麥千粒重隨播期推遲而上升;當x1在[-0.302，1] 范圍內，函數y值是隨x1的增大而減小，即在此范圍內小麥千粒重隨播期的推遲而下降。

2.3? ?密度對小麥產量及產量構成要素的影響

2.3.1? ? 密度對小麥產量的影響? ? 在方程（1）式中，將播期編碼設為零，則可得到關于密度與小麥產量的一元二次方程為：

從一元二次方程（9）式可以得出，函數以x=0.087為對稱軸，函數開口向下。當x2在[-1，0.087]范圍內，函數y值是隨x2的增大而增大，即在此范圍內小麥產量隨密度增加而上升;當x2在[0.087，1]范圍內，函數y值是隨x2的增大而減小，即在此范圍內小麥產量隨密度的增加而下降。

2.3.2? ? 密度對小麥成穗數的影響? ? 在方程（2）中，將播期編碼設為零，則可得到關于密度與小麥成穗數的一元二次方程為：

從一元二次方程（10）式可以得出，函數以x=0.502為對稱軸，函數開口向下。當x2在[-1，0.502]范圍內，函數y值是隨x2的增大而增大，即在此范圍內小麥成穗數隨密度增加而上升;當x2在[0.502，1]范圍內，函數y值是隨x2的增大而減小，即在此范圍內小麥成穗數隨密度的增加而下降。

2.3.3? ? 密度對小麥穗粒數的影響? ? 在方程（3）中，將播期編碼設為零，則可得到關于密度與小麥穗粒數的一元二次方程為：

從一元二次方程（11）式可以得出，函數以x=-0.172為對稱軸，函數開口向下。當x2在[-1，-0.172]范圍內，函數y值是隨x2的增大而增大，即在此范圍內小麥穗粒數隨密度增加而上升;當x2在[-0.172，1]范圍內，函數y值是隨x2的增大而減小，即在此范圍內小麥穗粒數隨密度的增加而下降。

2.3.4? ? 密度對小麥千粒重的影響? ? 在方程（4）中，將播期編碼設為零，則可得到關于密度與小麥千粒重的一元二次方程為：

y=41.158-0.098x2-0.046x22 （12）

從一元二次方程（12）式可以得出，函數以x=-1.065為對稱軸，函數開口向下。當x2在[-1，1]范圍內，函數y值是隨x2的增大而減小，即在此范圍內小麥千粒重隨密度增加而減小。

2.4? ?播期與密度最優選擇

將4個編碼值兩兩組合，形成16個全因子試驗，分別代入方程（1）式中的x1和x2中，即可模擬出16 個試驗結果（表4），其中包括已經實施的6個試驗在內。根據生產實踐，本試驗產量在6 300 kg/hm2以上的組合均為最優組合，從16個組合中選出滿足約束條件的組合為9個。對入選的9個最優組合進行平均數、標準差、標準誤、95%的置信區間以及播期和密度的計算。

經計算可得出，小麥產量在6 300 kg/hm2以上時，兩因素95%的置信區間x1為-0.855 9～0.054 7，x2為-0.275 4～0.672 8，相應的播期在9月17日至9月29日、密度為231.52～338.19 kg/hm2時效果為佳。

3? ?小結與討論

通過采用播期與密度的飽和D最優設計，分析播期與密度對冬小麥產量的影響。播期對產量有顯著的影響，增加密度能夠在一定范圍內增加穗數，進而增加產量。小麥產量在6 300 kg/hm2以上時，播期為9月17日至9月29日，密度為231.52～338.19 kg/hm2的效果為佳。

小麥產量的高低主要決定于單位面積有效穗數、穗粒數和千粒重，播期主要通過影響單株分蘗來影響穗數，最終影響產量，密度是決定穗數的主要因素[9 ]。研究表明，在一定范圍內，隨播期推遲，小麥產量下降。在播種密度同為150 kg/hm2時，9月15日播種的產量顯著高于10月10日的產量，這是因為小麥晚播使生長發育所需要的積溫不足，幼苗發育不良，造成穗數減少，進而影響小麥產量，這與張耀輝等[10? ]、呂麗華等[1 ]研究結果基本一致。一定范圍內，隨種植密度的增大，小麥產量增加。在播種日期同為9月15日時，播種密度為375 kg/hm2的產量高于150 kg/hm2的產量，這是因為，增加種植密度能夠在一定范圍內增加穗數，進而增加產量[12 - 14 ]。

本研究還表明，成穗數隨著播期的推遲而下降，在一定范圍內增加密度可以提高小麥成穗數，這與胡煥煥等[15 ]研究結果一致;穗粒數隨播期推遲而下降，在一定范圍內隨密度減小而上升，這與楊健等[16 ]研究結果一 致;千粒重在一定范圍內隨密度的增加而增加，播期對千粒重影響不大，這與楊健[16 ]等研究結果不同，這可能與輪選169品種特性和白銀沿黃灌區的氣候有關，有待進一步研究。

參考文獻：

[1] 楊文雄.? 中國西北春小麥[M].? 北京：中國農業出版社，2016.

[2] 王? ?煒，陳? ?琛，劉? ?風，等.? 播期播量對隴春31號籽粒性狀的影響[J].? 中國種業，2020（4）：53-56.

[3] 劉? ?萍，郭文善，徐月明，等.? 種植密度對中、弱筋小麥籽粒產量與品質的影響[J].? 麥類作物學報，2006，26（5）：116-121.

[4] 劉芳亮，張保軍，張正茂，等.? 播期和密度對普冰151小麥籽粒產量與品質的影響[J].? 干旱地區農業研究，2016，34（6）：1-7;14.

[5] 劉萬代，陳現勇，尹? ?鈞，等.? 播期和密度對冬小麥豫麥49-198群體性狀和產量的影響[J].? 麥類作物學報，2009，29（3）：464-469.

[6] 王? ?東，于振文，賈效成，等.? 播期對優質強筋冬小麥籽粒產量和品質的影響[J].? 山東農業科學，2004（2）：25-26.

[7] 張永成.? 飽和D最優設計方法在農業試驗中的應用[J].? 中國馬鈴薯，1997（3）：44-49.

[8] 高洪澤.? 運用Excel建立二次飽和D-最優設計回歸模型[J].? 現代農業科技，2010（22）：20-21.

[9] 趙廣才，常旭虹，王德梅，等.? 小麥生產概況及其發展[J].? 作物雜志，2018（4）：1-7.

[10] 張耀輝，宋建榮，岳維云，等.? 隴南雨養旱區播期與密度對冬小麥產量與品質的影響[J].? 干旱地區農業研究，2011，29（6）：74-78.

[11] 呂麗華，梁雙波，張麗華，等.? 不同小麥品種產量對冬前積溫變化的響應[J].? 作物學報，2016，42（1）：149-156.

[12] 柴芳梅，高甜甜，柴守璽，等.? 種植密度對甘肅不同生態區小麥產量形成的影響[J].? 作物雜志，2020（3）：154-160.

[13] 馬小黎，劉建華，鐘新榕.? 施肥量和種植密度對冬小麥葉面積及產量的調控[J].? 甘肅農業科技，2019（7）：19-25.

[14] 王希恩，韓? ?瑜，王德賢，等.? 播期與密度對旱地冬小麥天選55號產量及品質的影響[J].? 甘肅農業科技，2018（2）：51-55.

[15] 胡煥煥，劉麗萍，李瑞奇，等.? 播種期和密度對冬小麥品種河822產量形成的影響[J].? 麥類作物學報，2008，28（3）：490-495.

[16] 楊? ? 健，張保軍，毛建昌，等.? 播期與密度對冬小麥西農9871籽粒產量的影響[J].? 麥類作物學報，2011，31（3）：529-534.

（本文責編：陳? ? 偉）