向茂紅 車月英 包國舉

摘要：在不同肥力條件下，進行了玉米施肥“3414”試驗，采用Excel軟件回歸分析，依據(jù)《測土配方施肥技術(shù)規(guī)范》，以玉米相對產(chǎn)量50%以下、50%～70%、70%～90%、90%以上為標(biāo)準(zhǔn)，土壤豐缺養(yǎng)分級別分為高、中、低、極低4個級別，建立了試驗區(qū)土壤豐缺指標(biāo)體系，提出了高臺縣飼糧兼用玉米高、中、低肥力N、P、K配合使用及經(jīng)濟最佳施肥量。

關(guān)鍵詞：高臺縣;測土施肥; 飼糧兼用玉米;養(yǎng)分豐缺指標(biāo); 施肥量;最佳施肥量

中圖分類號：S513;S147.2 ? ? 文獻標(biāo)志碼：A ? ? 文章編號：1001-1463（2021）07-0004-06

doi：10.3969/j.issn.1001-1463.2021.07.002

Abstract：The effects of different planting methods of Lycium barbarum on its seedling and soil salinity were studied in Qaidam Basin. The results showed that it is all increased leaf area of Lycium barbarum which mulching with plastic film and straw at 30 cm underground， mulching with plastic film at 30 cm underground， mulching with straw at 30 cm underground， planting in furrows and planting on ridges， but had little effect on soil pH. Among them， the soil electrical conductivity of planting in furrows and planting on ridges was low， which greatly reduces the effect of salt on the plant， and can improve the survival rate of Lycium barbarum seedlings and the ability of leaf to absorb water. In conclusion， planting in furrows and ridges are the suitable planting methods of Lycium barbarum in saline-alkali land.

Key words：Saline-alkali land;Planting method;Lycium barbarum;Survival rate;Qaidam Basin

高臺縣是全國商品糧基地縣之一，2019年全縣農(nóng)作物播種面積3.65萬hm2，其中糧食作物播種1.86 萬hm2，全年糧食產(chǎn)量12 740.6萬kg，平均6 840 kg/hm2。近年來，高臺縣積極發(fā)展飼糧兼用型玉米，推行玉米整株青貯，發(fā)展草畜產(chǎn)業(yè)，不僅農(nóng)業(yè)結(jié)構(gòu)調(diào)整成效顯著，同時也讓農(nóng)民實現(xiàn)了增產(chǎn)增收。我們針對高臺縣不同區(qū)域土壤狀況、地貌類型等，在不同區(qū)域分高、中、低肥力條件開展了飼糧兼用玉米“3414”試驗，提出了高臺縣不同肥力條件土壤豐缺指標(biāo)體系，制定了土壤養(yǎng)分豐缺指標(biāo)下施肥量。

1 ? 試驗區(qū)概況

高臺縣位于河西走廊中部，黑河中游下段。南抵祁連，北依合黎，黑河縱貫東西。轄9鎮(zhèn)，136個行政村，總?cè)丝?5.83萬人。全縣總面積459.7 km2，耕地面積3.96萬hm2。全境海拔為1 260～3 140 m，大陸沙漠干旱型氣候，全年無霜期150 d左右，多年平均降水量103 mm，蒸發(fā)量2 000 mm左右，年日照時數(shù)3 088 h，熱量豐富，晝夜溫差大。

2 ? 材料與方法

2.1 ? 供試作物

指示作物為玉米品種正德305。

2.2 ? 供試肥料

供試肥料為尿素（新疆心連心能源化工生產(chǎn)，含N 46%），普通過磷酸鈣（云南玉溪星海肥業(yè)生產(chǎn)，含P2O5 12%），硫酸鉀（青海格爾木鉀肥廠生產(chǎn)，含K2O 50%）。

2.3 ? 試驗設(shè)計及方法

試驗在新壩鎮(zhèn)、南華鎮(zhèn)、駱駝城鎮(zhèn)、羅城鎮(zhèn)按照高肥力3個、中肥力4個、低肥力3個點開展。試驗采用“3414”完全實施方案，因素選擇氮、磷、鉀三個大量元素，水平設(shè)4個，0水平不施肥，2水平指當(dāng)?shù)赝扑]施肥水平，1水平=2水平×0.5，3水平=2水平×1.5（過量施肥水平），隨機區(qū)組排列，不設(shè)重復(fù)，小區(qū)間及試驗四周設(shè)有土埂，單排單灌，試驗四周設(shè)有保護行，小區(qū)面積保證在30 m2以上，試驗各處理及施肥量見表1。

2.4 ? 數(shù)據(jù)處理

采用Excel數(shù)據(jù)分析軟件，分析多元回歸、主效應(yīng)分析、交互效應(yīng)分析、邊際產(chǎn)量分析、最佳施肥量、最大施肥量。

3 ? 結(jié)果與分析

3.1 ? 土壤養(yǎng)分豐缺指標(biāo)

根據(jù)10個飼糧兼用玉米“3414”的肥效試驗，飼糧兼用玉米施肥體系豐缺指標(biāo)用缺肥區(qū)產(chǎn)量占全肥區(qū)產(chǎn)量百分數(shù)，即相對產(chǎn)量的高低來表征土壤養(yǎng)分的豐缺情況[1 - 2 ]。各試驗區(qū)處理6為全肥區(qū)（NPK），處理2、4、8為缺肥區(qū)（PK、NK、NP），收獲后計算產(chǎn)量相對數(shù)。以對數(shù)類型建立玉米相對產(chǎn)量和土壤養(yǎng)分之間的回歸曲線，得出兩者之間的回歸方程。

3.1.1 ? 相對產(chǎn)量與耕層堿解氮含量的關(guān)系

土壤堿解氮和玉米相對含量的回歸方程為Y=27.513 1 ln（x）-25.069 0（R2=0.400 3），相關(guān)性達到顯著水平（圖1）。

3.1.2 ? 相對產(chǎn)量與耕層有效磷含量的關(guān)系

土壤有效磷和玉米相對含量的回歸方程為Y=37.470 0 ln（x）-28.320 0（R2=0.510 0），相關(guān)性達到顯著水平（圖2）。

3.1.3 ? 相對產(chǎn)量與耕層速效鉀含量的關(guān)系

土壤速效鉀和玉米相對含量的回歸方程為Y=43.628 0 ln（x） - 53.239 0（R2 = 0.806 6），相關(guān)性達到極顯著水平（圖3）。

依據(jù)《甘肅省測土配方施肥技術(shù)規(guī)范》，用相對產(chǎn)量50%、70%、90%分別計算出對應(yīng)土壤養(yǎng)分含量，即得到玉米土壤養(yǎng)分豐缺指標(biāo)體系（表3），分極低（相對產(chǎn)量<50%）、低（相對產(chǎn)量50%～70%）、中（相對產(chǎn)量70%～90%）、高（相對產(chǎn)量≥90%）4個水平。同時確定對應(yīng)的肥料施用量。

3.2 ? 模型的建立

肥力效應(yīng)方程式模擬施肥量和產(chǎn)量之間的數(shù)量關(guān)系模擬三元二次效應(yīng)方程：Y=b0+

b1 x1 ?+ b2 ?x2 + b3 ?x3 +b4 x1+ b5 ?x2 + b6 ?x3 + b7 ?x1 x2 + b8 ?x1 x3 + b9 x2 x3，建立數(shù)學(xué)模型對玉米施肥進行擬合。將試驗方案設(shè)計與得到產(chǎn)量結(jié)果在Excel中進行回歸分析，得到肥料效應(yīng)方程。

經(jīng)分析，玉米肥料效應(yīng)回歸模型全部為x1、x2、x3三元二次多項式回歸模型，F(xiàn)值 > F0.05，高肥力F=154.05 > F0.05=0.001，中肥力F=23.34 > F0.05=0.004，低肥力F=19.58 > F0.05= 0.006。經(jīng)顯著性檢驗達到極顯著水平，說明肥料效應(yīng)回歸模型具有實際應(yīng)用價值，可以用于高臺縣飼糧兼用型玉米施肥量的預(yù)測。

3.3 ? 單因素分析

x1、x2、x3各施肥因子對產(chǎn)量影響的獨立效應(yīng)，所得回歸方程為：

Y=b0+b1x1+b4x12

Y=b0+b2 x2+b5 x22

Y=b0+b3 x3+b6 x32

經(jīng)分析，各產(chǎn)量等級F值均大于F0.05，反映了玉米產(chǎn)量與主控因子x1（氮）、x2（磷）、x3（鉀）、P、K的相關(guān)因子達到顯著水平，說明玉米產(chǎn)量氮、磷、鉀肥之間有顯著回歸關(guān)系。經(jīng)回歸，R2值最低的低肥力鉀回歸復(fù)相關(guān)系數(shù)為0.887 1，其他較高，表明擬合性好。

從表5看出，增產(chǎn)效果居首位的為氮肥，中肥力下施N增產(chǎn)1 260.0 kg/hm2，增產(chǎn)率29.7%;低肥力下磷肥的增產(chǎn)幅度次之，施P2O5增產(chǎn)557.4 kg/hm2，增產(chǎn)率15.8%;高肥力下施K2O最高增產(chǎn)489.4 kg/hm2，增產(chǎn)率8.1%，增產(chǎn)效果順序為氮>磷>鉀。從單位養(yǎng)分的增產(chǎn)效果看， 鉀（K2O）素的增產(chǎn)效果最高，為7.36～8.16 kg/kg;磷（P2O5）素的增產(chǎn)效果次之，為4.64～6.56 kg/kg;氮（N）素的增產(chǎn)效果最低，為3.23～5.57 kg/kg。

3.4 ? 交互作用分析

營養(yǎng)元素間的相互作用必然影響作物營養(yǎng)與產(chǎn)量。對三元二次肥料效應(yīng)方程， 分別采用“降維法”固定一個因素在零水平， 這相當(dāng)于在特定條件下所做的交互效應(yīng)試驗， 所得回歸方程為：

Y=b0+b1 x1+b2 x2+b4 x1x2+b7 x1 2+b8 x2 2

Y=b0+b1 x1+b3 x3+b5 x1x3+b7 x1 2+b9 x3 2

Y=b0+b2 x2+b3 x3+b6 x2 x3+b8 x2 2+b9 x32

以處理2、處理3、處理4、處理5、處理6、處理7、處理11、處理12為NP交互效應(yīng)研究，處理2、處理3、處理6、處理8、處理9、處理10、處理11、處理13為NK交互效應(yīng)研究，處理4、處理5、處理6、處理7、處理8、處理9、處理10、處理14為PK交互研究。在不同肥料的交互作用中，NP的交互作用在低肥力條件下影響最大，NK的交互作用在高肥力條件下影響最大，PK的交互作用在中肥力條件下影響最大。

3.5 ? 模型優(yōu)化

作物最佳施肥量是指邊際成本等于邊際利潤、肥料的邊際利潤為零時的施肥量[3 - 6 ]。運用Excel分析軟件，對高臺縣10個玉米“3414”肥效試驗進行分析，肥料效應(yīng)方程式擬合結(jié)果呈顯著水平。用三元二次肥料效應(yīng)函數(shù)分別求產(chǎn)量對N、P、K三個因素求偏導(dǎo)， 即得其邊際產(chǎn)量。邊際產(chǎn)量等于肥料的單價除以玉米的單價施用的施肥量為最優(yōu)施肥量。本試驗采用的N、P2O5、K2O和玉米的價格分別為4.89元/kg、5.83元/kg、9.42元/kg、2.00元/kg。經(jīng)過Excel數(shù)據(jù)分析軟件對數(shù)據(jù)進行統(tǒng)計和分析，建立高臺縣不同肥力條件下最佳施肥料量及最大施肥量（表6）。高臺縣飼糧兼用玉米平均經(jīng)濟最佳施肥量及配比為：高肥力下N 251.55 kg/hm2、P2O5 105.77 kg/hm2、K2O 51.94 kg/hm2，氮、磷、鉀質(zhì)量比為1∶0.42∶0.21，產(chǎn)量為10 083.9 kg/hm2。中肥力下N 208.38 kg/hm2、P2O5 118.68 kg/hm2、K2O 49.22 kg/hm2，氮、磷、鉀質(zhì)量比為1∶0.57∶0.24，產(chǎn)量為 ? ? 7 547.28 kg/hm2。低肥力下N 182.21 kg/hm2、P2O5 124 kg/hm2、K2O 53.71 kg/hm2，氮、磷、鉀質(zhì)量比為1∶0.68∶0.29，產(chǎn)量為5 831.82 kg/hm2。

4 ? 小結(jié)

飼糧兼用玉米相對產(chǎn)量與土壤堿解氮、速效磷、速效鉀含量存在對數(shù)函數(shù)關(guān)系，且所得對數(shù)相關(guān)方程達顯著水平。增產(chǎn)效果居首位的為氮肥，對產(chǎn)量影響的大小順序為：N >P >K。以玉米相對產(chǎn)量50%、50%～70%、70%～90%、90%以上為臨界值標(biāo)準(zhǔn)，土壤養(yǎng)分豐缺級別分為極低、低、中、高4個級別，建立了該地區(qū)土壤豐缺指標(biāo)體系。

分析得出高臺縣飼糧兼用玉米高中低肥力N、P、K配合施用質(zhì)量比，高肥力為1∶0.42∶0.21，中肥力為1∶0.57∶0.24，高肥力為1∶0.68∶0.29。

參考文獻：

[1] 張明學(xué)，馬志超，周倉軍，等. ?關(guān)中玉米氮磷鉀養(yǎng)分豐缺指標(biāo)及最佳施肥量研究[J]. ?西北農(nóng)林科技大學(xué)學(xué)報，2015（11）：146-151.

[2] 趙宗海，李志明. ?定西市安定區(qū)測土配方施肥土壤養(yǎng)分豐缺指標(biāo)及施肥比例研究[J]. ?現(xiàn)代農(nóng)業(yè)科技，2014（18）：154-155.

[3] 孫建好，李偉綺，趙建華. ?高臺縣小麥及玉米施肥現(xiàn)狀調(diào)查與評價[J]. ?甘肅農(nóng)業(yè)科技，2019（6）：51-56.

[4] 李 ? 旭，杜曉明，洛俊峰，等. ?基于“3414”試驗設(shè)計的玉米氮、磷、鉀效應(yīng)研究[J]. ?吉林農(nóng)業(yè)，2019（2）：62-63.

[5] 潘進紅，王艷群，劉建平. ?基于“3414+1 ”模型的春玉米產(chǎn)量效應(yīng)及推薦施肥量研究[J]. ?河北北方學(xué)院學(xué)報，2015（6）：24-29.

[6] 胡建軍，張月琴，溫學(xué)飛，等. ?基于“3414”模型研究寧夏鹽池縣玉米氮磷鉀施肥效應(yīng) [J]. ?水土保持研究， 2011（8）：268-272.

（本文責(zé)編：陳 ? 珩）