楊錫永++段必挺++王朝軍++楊海濤++楊堃++張平++戴曙偉++劉洪華

摘要 本文利用因子分析法和系統(tǒng)聚類分析法對(duì)施甸縣781個(gè)土壤樣品的養(yǎng)分進(jìn)行了分析與綜合評(píng)價(jià)，旨在為施甸縣煙區(qū)特色優(yōu)質(zhì)煙葉生產(chǎn)提供科學(xué)合理的施肥方案。結(jié)果表明，土壤肥力最高的樣本17個(gè)，占2.18%；土壤肥力較高的樣本119個(gè)，占15.24%；土壤肥力較低的樣本423個(gè)，占54.16%；土壤肥力最低的樣本222個(gè)，占28.42%。對(duì)于土壤肥力處于較高水平以上的地塊，建議農(nóng)家肥施用量為4 500 kg/hm2左右、復(fù)合肥施用量控制在純N 60～75 kg/hm2；對(duì)于土壤肥力處于較低水平的地塊，建議農(nóng)家肥施用量為7 500～15 000 kg/hm2、復(fù)合肥施用量控制在純N 90～120 kg/hm2。

關(guān)鍵詞 土壤養(yǎng)分；因子分析；系統(tǒng)聚類分析；綜合評(píng)價(jià)；云南施甸；煙區(qū)

中圖分類號(hào) S158.2 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-5739（2017）19-0172-03

Abstract The nutrient of 781 soil samples were analyzed and evaluated by the Factor Analysis and Hierachical Cluster Analysis in this paper，in order to provide scientific and reasonable fertilization scheme for characteristic and high-quality tobacco production in tobacco-growing area of Shidian County. The results showed that the soil fertility of 17 samples were the highest，accounting for 2.18% of the total；119 samples were the second，accounting for 15.24% of the total；423 samples were the third，accounting for 54.16% of the total；222 samples were the lowest，accounting for 28.42% of the total. Therefore，the 4 500 kg/hm2 manure and 60-75 kg/hm2 pure nitrogen of compound-fertilizer were recommended for the higher soil fertility of fields，and 7 500-15 000 kg/hm2 manure and 90-120 kg/hm2 pure nitrogen of compound-fertilizer were recommended for the lower soil fertility of fields.

Key words soil nutrition；Factor Analysis；Hierachical Cluster Analysis；integrated evaluation；Shidian Yunnan；tobacco-growing area

土壤對(duì)煙葉質(zhì)量的影響最大，優(yōu)質(zhì)煙葉生產(chǎn)的基礎(chǔ)之一就是肥力適宜的土壤[1]，土壤肥力可以綜合體現(xiàn)出土壤多種性質(zhì)[2]。在同一區(qū)域內(nèi)、統(tǒng)一的栽培和調(diào)制技術(shù)條件下，相同煙草品種產(chǎn)出的煙葉會(huì)因土壤不同而存在顯著的質(zhì)量差異[3]。煙草的生長(zhǎng)發(fā)育、煙葉產(chǎn)量和質(zhì)量都受到土壤養(yǎng)分供應(yīng)的影響[4]。許多研究表明，煙草根系的生長(zhǎng)發(fā)育、煙株的抗逆性及煙葉品質(zhì)與土壤鉀含量有密切關(guān)系；煙葉根系生長(zhǎng)條件受pH值影響，pH值還對(duì)煙葉刺激性、雜氣、香氣等指標(biāo)有一定影響；土壤氮和有機(jī)質(zhì)含量適宜，有利于提高煙葉品質(zhì)；土壤中鎂、氯等元素含量對(duì)煙株生長(zhǎng)和煙葉品質(zhì)有重要的影響[5-7]。因此，如何利用土壤養(yǎng)分狀況科學(xué)、合理地評(píng)價(jià)土壤肥力，為指導(dǎo)烤煙種植和合理施肥提供可靠依據(jù)，顯得尤為重要。本文采用了因子分析法和系統(tǒng)聚類分析法對(duì)施甸縣煙區(qū)近10年土壤養(yǎng)分狀況進(jìn)行綜合評(píng)價(jià)，以期促進(jìn)科學(xué)平衡施肥和充分發(fā)揮區(qū)域土壤潛力，為提高肥料利用效率和煙草品質(zhì)提供科學(xué)依據(jù)。

1 材料與方法

1.1 土樣采集

于2008—2017年在施甸縣13個(gè)烤煙種植鄉(xiāng)鎮(zhèn)采集土樣781個(gè)。鄉(xiāng)鎮(zhèn)具體包括姚關(guān)、擺榔、萬興、酒房、太平、何元、水長(zhǎng)、由旺、老麥、仁和、甸陽、木老元、舊城，樣本數(shù)分別為96、20、56、80、118、56、42、36、93、75、75、18、16個(gè)。采用隨機(jī)布點(diǎn)法采集植煙土壤耕層0～20 cm土壤，采集時(shí)間為煙草移栽以及施底肥以前，以充分反映采樣地塊的真實(shí)養(yǎng)分狀況和供肥能力。

1.2 測(cè)定項(xiàng)目及方法

測(cè)定項(xiàng)目包括土壤有機(jī)質(zhì)、堿解氮、速效磷、速效鉀、有效硼、有效鋅、有效鎂、水溶性氯含量及pH值，測(cè)定方法均參照文獻(xiàn)[8]。

1.3 分析軟件

試驗(yàn)數(shù)據(jù)統(tǒng)計(jì)分析采用軟件Excel 2003和SPSS 19.0。

2 結(jié)果與分析

2.1 標(biāo)準(zhǔn)化土壤養(yǎng)分指標(biāo)相關(guān)系數(shù)矩陣

由表1可知，大部分養(yǎng)分指標(biāo)之間存在相關(guān)性，并且達(dá)到顯著或者極顯著水平。因此，如果土壤肥力綜合評(píng)價(jià)直接采用土壤養(yǎng)分指標(biāo)進(jìn)行，將因?yàn)樾畔⒅丿B而產(chǎn)生較大偏差，影響評(píng)價(jià)結(jié)果的準(zhǔn)確性。

2.2 土壤養(yǎng)分指標(biāo)的因子分析

2.2.1 Bartlett球度檢驗(yàn)。為了判斷土壤養(yǎng)分?jǐn)?shù)據(jù)是否適宜進(jìn)行因子分析，采用Bartlett球度檢驗(yàn)法對(duì)養(yǎng)分指標(biāo)原始數(shù)據(jù)進(jìn)行檢驗(yàn)，結(jié)果表明，其概率為0<0.05（顯著性水平），即數(shù)據(jù)適宜進(jìn)行因子分析。

2.2.2 主成分分析。對(duì)表1中的9個(gè)土壤養(yǎng)分指標(biāo)進(jìn)行因子分析，提取出5個(gè)主成分，結(jié)果見表2?？梢钥闯觯@5個(gè)主成分方差貢獻(xiàn)率分別為35.61%、13.97%、12.08%、9.97%、8.52%，累積貢獻(xiàn)率達(dá)80.15%，可以體現(xiàn)各個(gè)原變量的差異、特性及關(guān)系。

進(jìn)行方差極大正交旋轉(zhuǎn)，得到因子載荷矩陣見表3。可以看出，第1主成分主要反映了有機(jī)質(zhì)、堿解氮的作用；第2主成分主要反映了速效磷、速效鉀、有效硼的作用；第3主成分主要反映了有效鋅、有效鎂的作用；第4主成分主要反映了土壤pH值的作用；第5主成分主要反映了水溶性氯的作用。因?yàn)橐陨?個(gè)主因子的累積貢獻(xiàn)率達(dá)到80.15%，所以利用它們對(duì)土壤肥力進(jìn)行不同方面的綜合評(píng)價(jià)是合理、可行的。

2.2.3 因子得分函數(shù)模型構(gòu)建。主成分得分系數(shù)矩陣見表4，由此對(duì)5個(gè)主成分得分進(jìn)行計(jì)算，即得分系數(shù)與相對(duì)應(yīng)的標(biāo)準(zhǔn)化變量（Zi）之積，其數(shù)學(xué)模型如下：

F1=-0.031Z1+0.534Z2+0.473Z3+0.006Z4-0.323Z5+0.034Z6+0.187Z7-0.109Z8-0.052Z9；

F2=0.022Z1-0.218Z2-0.107Z3+0.369Z4+0.641Z5+0.426Z6-0.089Z7-0.106Z8+0.002Z9；

F3=-0.085Z1-0.004Z2-0.049Z3+0.040Z4-0.056Z5-0.198Z6+0.462Z7+0.813Z8+0.013Z9；

F4=1.000Z1+0.049Z2-0.081Z3+0.004Z4+0.062Z5-0.048Z6-0.032Z7-0.081Z8-0.017Z9；

F5=-0.017Z1-0.042Z2+0.043Z3-0.053Z4+0.122Z5-0.097Z6-0.158Z7+0.099Z8+0.978Z9。

將各養(yǎng)分標(biāo)準(zhǔn)化處理后的數(shù)據(jù)帶入上述5個(gè)公式，可以求出各樣品的因子得分F1、F2、F3、F4、F5的值。

2.2.4 主成分因子得分。依據(jù)主成分因子的權(quán)重計(jì)算方法，以5個(gè)主因子各自的方差貢獻(xiàn)率與入選因子累計(jì)貢獻(xiàn)率的比例為權(quán)數(shù)進(jìn)行加權(quán)綜合，可以得出施甸縣煙區(qū)土壤肥力綜合得分評(píng)價(jià)數(shù)學(xué)模型如下：

F=0.444×F1+0.174×F2+0.151×F3+0.124×F4+0.106×F5

將數(shù)據(jù)帶入公式，得出土壤肥力綜合得分F值，F(xiàn)值的高、低相應(yīng)地表征土壤肥狀況的好、差。

2.3 土壤肥力的聚類分析

采用重心法對(duì)土壤樣本主成分綜合得分進(jìn)行系統(tǒng)聚類分析，將781個(gè)樣本的綜合得分劃為4類（表5）。Ⅰ、Ⅱ、Ⅲ、Ⅳ類綜合得分F分別為[1.260 6，2.050 5]、[0.470 7，1.260 6）、[-0.319 2，0.470 7）、（-1.109 1，-0.319 2]，樣本數(shù)分別為17、119、423、222個(gè)，分別占總樣本數(shù)的2.18%、15.24%、54.16%、28.42%。根據(jù)F值可以看出，Ⅰ、Ⅱ、Ⅲ、Ⅳ類分別表示土壤肥力最高、較高、較差、差。

在全縣范圍內(nèi)，17.42%的土壤肥力較高（Ⅰ類、Ⅱ類），還有82.58%的土壤肥力不足（Ⅲ類、Ⅳ類）。因此，全縣大部分植煙土壤的肥力較低。

3 結(jié)論與討論

因子分析和系統(tǒng)聚類分析結(jié)果表明，施甸縣土壤肥力較高的樣本有17個(gè)，占2.18%；土壤肥力較高的有119個(gè)，占15.24%；土壤肥力較低的有423個(gè)，占54.16%；土壤肥力最低的有222個(gè)，占28.42%。由此可見，在全縣范圍內(nèi)，全縣很大部分植煙土壤的肥力仍較低，占82.58%，僅有17.42%的土壤肥力處于較高水平以上。因此，要依據(jù)不同地塊的土壤分析結(jié)果，制訂相應(yīng)的施肥方案，做到科學(xué)合理施肥，對(duì)于土壤肥力處于較高水平以上的地塊，要施用適量的有機(jī)肥和復(fù)合肥，建議農(nóng)家肥施用量為4 500 kg/hm2左右、復(fù)合肥施用量控制在純N 60～75 kg/hm2；對(duì)于土壤肥力處于較低水平的，要確保有機(jī)肥和復(fù)合肥的足量、配合施用，建議農(nóng)家肥施用量為7 500～15 000 kg/hm2、復(fù)合肥施用量控制在純N 90～120 kg/hm2。

4 參考文獻(xiàn)

[1] 王彥亭，謝劍平，李志宏.中國(guó)煙草種植區(qū)劃[M].北京：科學(xué)出版社，2010：28-35.

[2] 趙文田，劉海輪，呂家瓏，等.陜西隴縣烤煙種植區(qū)土壤養(yǎng)分的測(cè)定與肥力綜合評(píng)價(jià)[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)，2009，37（9）：123-128.

[3] 中國(guó)農(nóng)業(yè)科學(xué)院煙草研究所.中國(guó)煙草栽培學(xué)[M].上海：科學(xué)技術(shù)出版社，1987：1-133.

[4] 曹志洪.優(yōu)質(zhì)烤煙生產(chǎn)的土壤與施肥[M].南京：江蘇科學(xué)技術(shù)出版社，1991.

[5] 王林，盧秀萍，許自成，等.瀏陽植煙土壤肥力狀況的綜合評(píng)價(jià)[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2006，40（6）：597-601.

[6] 胡國(guó)松，鄭偉，王震東，等.烤煙營(yíng)養(yǎng)原理[M].北京：科學(xué)出版社，2000：57-89.

[7] 黃燕翔，劉淑欣，熊德中，等.福建煙區(qū)土壤條件與烤煙品質(zhì)的關(guān)系[J].福建農(nóng)業(yè)大學(xué)學(xué)報(bào)，1995，24（2）：201-204.

[8] 鮑士旦.土壤農(nóng)化分析[M].北京：中國(guó)農(nóng)業(yè)出版社，2000：64-237.