摘要 以Z市2013—2023年糧食生產(chǎn)相關(guān)數(shù)據(jù)為基礎(chǔ)，選擇糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、糧食生產(chǎn)化肥施用量、平均氣溫等14項(xiàng)指標(biāo)，運(yùn)用灰色關(guān)聯(lián)度分析和相關(guān)性分析法，分析影響研究區(qū)糧食產(chǎn)量的關(guān)鍵因子，分析糧食單產(chǎn)與部分糧食作物面積占糧食種植面積比例的相關(guān)性。結(jié)果表明，研究區(qū)糧食產(chǎn)量與糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、薯類面積占糧食面積比例、糧食播種面積關(guān)聯(lián)度較大；與糧食生產(chǎn)化肥施用量、糧食播種面積相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.01），與玉米面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05），與薯類面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05）；糧食單產(chǎn)與稻谷面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.01），與豆類面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05）。綜合考慮化肥施用零增長(zhǎng)、耕地面積有限等因素，增加稻谷和玉米種植面積占糧食面積的比例，同時(shí)適當(dāng)降低薯類面積占糧食面積的比例利于糧食產(chǎn)量提高，有利于推動(dòng)研究區(qū)糧食增產(chǎn)。

關(guān)鍵詞 相關(guān)性分析；耕地面積；糧食產(chǎn)量；化肥施用；地膜使用量

中圖分類號(hào) F326.11" "文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 1007-7731（2024）11-0001-05

Analysis of factors influencing grain production in Z City

MU Shaodong1" " WANG Fang2" " WU Quanhua2" " LI Zhengli1" " ZHANG Xuejiao3" " SUN Jing1

（1Zunyi Planting Industry Development Service Center， Zunyi 563000， China;2Zunyi Agriculture and Rural Bureau，" Zunyi 563000， China;3Guizhou Agriculture and Rural Department， Guiyang 550000， China）

Abstract Based on the data related to grain production in Z City from 2013 to 2023， 14 factors including the total power of agricultural machinery in grain production， the amount of chemical fertilizers used in grain production， and the average temperature were selected. The grey correlation method and correlation analysis method were used to analyze the key factors affecting the grain yield in the study area， and the correlation between grain yield and the proportion of some grain crops area to grain planting area was analyzed. The results indicated that there was a strong correlation between the grain yield in the study area and the total power of agricultural machinery in grain production， the proportion of potato area to grain area， and the sowing area of grain. There was a statistically significant correlation （Plt;0.01） between the amount of chemical fertilizer applied in grain production and the sowing area of grain， a statistically significant correlation （Plt;0.05） with the proportion of corn area to grain area， and a statistically significant correlation （Plt;0.05） with the proportion of potato area to grain area. The correlation between grain yield and the proportion of rice area to grain area was statistically significant （Plt;0.01）， and the correlation with the proportion of legume area to grain area was statistically significant （Plt;0.05）.Taking into account factors such as zero growth in fertilizer application and limited arable land area， increasing the proportion of rice and corn planting area to grain area， while reducing the proportion of potato planting area to grain area， was beneficial for improving grain yield and promoting grain yield increase in the research area.

Keywords correlation analysis; cultivated land area; grain yield; fertilization; usage of plastic film

保障糧食安全的關(guān)鍵是糧食生產(chǎn)的穩(wěn)定發(fā)展[1]。糧食生產(chǎn)受經(jīng)濟(jì)、社會(huì)、自然、人文、技術(shù)和人力資源等多因素的綜合影響，分析糧食產(chǎn)量影響因素及關(guān)系對(duì)糧食生產(chǎn)有重要意義[2-3]。近年來，在糧食產(chǎn)量影響因素及其作用方面，有關(guān)學(xué)者圍繞氣候條件[4-6]、自然災(zāi)害[7-10]、耕地質(zhì)量[11-14]和生產(chǎn)要素投入[1，3，15-20]等方面分析了影響糧食產(chǎn)量的因素，但由于研究角度和方法不同，且各地糧食生產(chǎn)條件不同，研究結(jié)果存在較大差異。Z市糧食生產(chǎn)面積較大，農(nóng)業(yè)生產(chǎn)水平較高，基礎(chǔ)條件較好，分析影響糧食產(chǎn)量的因素，探明糧食產(chǎn)量和各因素之間的關(guān)系，對(duì)推動(dòng)該市糧食安全穩(wěn)定供給有重要意義。本文以該地區(qū)2013—2023年糧食生產(chǎn)相關(guān)數(shù)據(jù)為基礎(chǔ)，通過灰色關(guān)聯(lián)度分析和相關(guān)性分析等統(tǒng)計(jì)方法，探析糧食產(chǎn)量與影響因素之間的關(guān)系，明確研究區(qū)域糧食增產(chǎn)關(guān)鍵因子，為糧食生產(chǎn)提供參考。

1 材料與方法

1.1 數(shù)據(jù)來源

2013—2022年部分農(nóng)業(yè)生產(chǎn)條件指標(biāo)、農(nóng)作物播種面積和糧食產(chǎn)量等數(shù)據(jù)來源于研究區(qū)統(tǒng)計(jì)年鑒，2023年耕地面積、糧食播種面積和糧食產(chǎn)量數(shù)據(jù)由研究區(qū)統(tǒng)計(jì)調(diào)查隊(duì)提供。

1.2 研究方法

基于SPSS 26.0和Excel 2019軟件，運(yùn)用灰色關(guān)聯(lián)度分析和相關(guān)性分析法對(duì)糧食產(chǎn)量與影響因素的關(guān)聯(lián)度，糧食產(chǎn)量與影響因素的相關(guān)性進(jìn)行分析。數(shù)據(jù)計(jì)算主要涉及以下公式。

年某糧食作物面積占糧食總播種面積比例（%）=年該糧食作物播種面積（萬hm2）/年糧食總播種面積（萬hm2）×100 （1）

年某糧食作物單產(chǎn)（t/hm2）=年該糧食作物產(chǎn)量（萬t）/年該糧食作物播種面積（萬hm2） （2）

年糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力（萬kW）=年農(nóng)機(jī)總動(dòng)力（萬kW）×年糧食總播種面積（萬hm2）/該年農(nóng)作物總播種面積（萬hm2） （3）

年糧食生產(chǎn)化肥施用量（萬t）=年農(nóng)用化肥施用量（萬t）×年糧食總播種面積（萬hm2）/年農(nóng)作物總播種面積（萬hm2） （4）

年糧食生產(chǎn)地膜使用量（萬t）=年地膜使用量（萬t）×年糧食總播種面積（萬hm2）/年農(nóng)作物總播種面積（萬hm2） （5）

2 結(jié)果與分析

2.1 糧食產(chǎn)量與影響因素的關(guān)聯(lián)度分析

根據(jù)研究區(qū)2013—2022年統(tǒng)計(jì)年鑒，對(duì)糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、糧食生產(chǎn)化肥施用量、平均氣溫、年平均降水量、糧食生產(chǎn)地膜使用量、農(nóng)村平均每戶勞動(dòng)力、農(nóng)村居民家庭人均年種植業(yè)生產(chǎn)支出、稻谷面積占糧食面積比例、玉米面積占糧食面積比例、高粱面積占糧食面積比例、小麥面積占糧食面積比例、豆類面積占糧食面積比例、薯類面積占糧食面積比例和糧食播種面積14項(xiàng)指標(biāo)和糧食產(chǎn)量分別進(jìn)行統(tǒng)計(jì)（表1）。

將14項(xiàng)指標(biāo)與糧食產(chǎn)量的關(guān)聯(lián)度進(jìn)行分析（表2），結(jié)果表明，糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、糧食播種面積、薯類面積占糧食面積比例與糧食產(chǎn)量的關(guān)聯(lián)度較大，分別為0.977 0、0.808 8和0.783 7，糧食生產(chǎn)地膜使用量、農(nóng)村居民家庭人均年種植業(yè)生產(chǎn)支出、年平均降水量與糧食產(chǎn)量的關(guān)聯(lián)度較小，分別為0.746 7、0.746 5和0.448 2。

2.2 糧食產(chǎn)量與影響因素的相關(guān)性分析

對(duì)研究區(qū)2013—2022年糧食產(chǎn)量與糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、糧食生產(chǎn)化肥施用量、平均氣溫、年平均降水量、糧食生產(chǎn)地膜使用量、農(nóng)村平均每戶勞動(dòng)力、農(nóng)村居民家庭人均年種植業(yè)生產(chǎn)支出、稻谷面積占糧食面積比例、玉米面積占糧食面積比例、高粱面積占糧食面積比例、小麥面積占糧食面積比例、豆類面積占糧食面積比例、薯類面積占糧食面積比例和糧食播種面積14項(xiàng)指標(biāo)進(jìn)行相關(guān)性分析（表3）。結(jié)果表明，糧食產(chǎn)量與糧食生產(chǎn)化肥施用量、糧食播種面積相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.01），與玉米面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05），與薯類面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05），與其他因素相關(guān)性無統(tǒng)計(jì)學(xué)意義（Pgt;0.05）。這說明提高糧食生產(chǎn)化肥施用量、糧食播種面積和玉米面積占糧食面積比例，一定程度降低薯類面積占糧食面積比例，有利于糧食產(chǎn)量的增加。

依靠擴(kuò)大種植面積來增加糧食產(chǎn)量、保障糧食安全的空間有限，在當(dāng)前農(nóng)業(yè)生產(chǎn)條件下，增加糧食單產(chǎn)是提高糧食產(chǎn)量的重要途徑之一[23]，為進(jìn)一步分析糧食產(chǎn)量的提高方法，本研究對(duì)糧食單產(chǎn)與部分糧食作物占糧食種植面積比例的相關(guān)性進(jìn)行分析。對(duì)研究區(qū)域各地2022和2023年糧食單產(chǎn)與稻谷、玉米、高粱、豆類和薯類面積占糧食面積比例具體數(shù)據(jù)及相關(guān)性分析表明（表4—5），糧食單產(chǎn)與稻谷面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.01），與豆類面積占糧食面積比例相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05）。這說明提高稻谷面積占糧食面積比例和一定程度降低豆類面積占糧食面積比例，有利于糧食單產(chǎn)水平提高，進(jìn)一步增加糧食產(chǎn)量。

3 結(jié)論與討論

2013—2022年，Z市糧食生產(chǎn)農(nóng)機(jī)總動(dòng)力、薯類面積占糧食面積比例、糧食播種面積與糧食產(chǎn)量的關(guān)聯(lián)度較大；相關(guān)性分析結(jié)果表明，糧食產(chǎn)量與糧食生產(chǎn)化肥施用量、糧食播種面積的相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.01），與玉米面積占糧食面積比例的相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05），與薯類面積占糧食面積比例的相關(guān)性具有統(tǒng)計(jì)學(xué)意義（Plt;0.05），增加糧食生產(chǎn)化肥施用量、糧食播種面積和玉米面積占糧食面積比例，適當(dāng)減少薯類面積占糧食面積比例利于糧食產(chǎn)量的增加。同時(shí)，增加稻谷面積占糧食面積比例和適當(dāng)降低豆類面積占糧食面積比例利于糧食單產(chǎn)水平提高，進(jìn)一步增加糧食產(chǎn)量。由于依靠擴(kuò)大糧食面積來增加糧食產(chǎn)量、保障糧食安全的空間有限，擴(kuò)種大豆是保障糧食安全的重要途徑之一，降低大豆種植面積可能性不大，且考慮到化肥使用邊際遞減效應(yīng)[24]及當(dāng)前倡導(dǎo)施用有機(jī)肥料及生物肥料推進(jìn)“化肥零增長(zhǎng)”行動(dòng)[25]等因素，較難大規(guī)模增加糧食播種面積，研究區(qū)域糧食產(chǎn)量與化肥施用量的具體關(guān)系還須進(jìn)一步研究。綜上，在農(nóng)業(yè)生產(chǎn)條件一定的情況下，適當(dāng)增加稻谷和玉米種植面積占糧食面積的比例，同時(shí)一定程度降低薯類面積占糧食面積的比例有利于糧食產(chǎn)量提高，同時(shí)對(duì)推動(dòng)研究區(qū)糧食增產(chǎn)具有積極作用。

參考文獻(xiàn)

[1] 劉曉敏，石喆. 山東省糧食產(chǎn)量影響因素分析及增產(chǎn)對(duì)策研究[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2016，32（12）：171-174.

[2] 劉浩然，吳克寧，宋文，等. 黑龍江糧食產(chǎn)能及其影響因素研究[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃，2019，40（7）：164-170.

[3] 張莉，張賢，梁根宏，等. 貴州糧食產(chǎn)能影響因素分析及其趨勢(shì)預(yù)測(cè)[J]. 作物研究，2021，35（4）：368-375.

[4] 李志紅，謝東輝. 影響我國(guó)氣候旱區(qū)糧食生產(chǎn)能力的灰色關(guān)聯(lián)分析[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃，2016，37（8）：124-129.

[5] 王紅麗，張緒成，魏勝文. 氣候變化對(duì)西北半干旱區(qū)旱作農(nóng)業(yè)的影響及解決途徑[J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2015，32（6）：517-524.

[6] 雷秋良，徐建文，姜帥，等. 氣候變化對(duì)中國(guó)主要作物生育期的影響研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2014，30（11）：205-209.

[7] 龍方，楊重玉，彭澧麗. 自然災(zāi)害對(duì)中國(guó)糧食產(chǎn)量影響的實(shí)證分析：以稻谷為例[J]. 中國(guó)農(nóng)村經(jīng)濟(jì)，2011（5）：33-44.

[8] 莊道元，陳超，趙建東. 不同階段自然災(zāi)害對(duì)我國(guó)糧食產(chǎn)量影響的分析：基于31個(gè)省市的面板數(shù)據(jù)[J]. 軟科學(xué)，2010，24（9）：39-42.

[9] 劉佩，劉峰貴，周強(qiáng)，等. 河南省水旱災(zāi)害時(shí)空分布特征及與糧食產(chǎn)量關(guān)系[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2011，27（29）：290-295.

[10] 洪名勇，周歡，劉洪. 自然災(zāi)害對(duì)貴州省糧食波動(dòng)的影響研究[J]. 農(nóng)業(yè)現(xiàn)代化研究，2016，37（1）：35-42.

[11] 方彥杰，張緒成，侯慧芝，等. 近20年來甘肅省黃土高原區(qū)糧食產(chǎn)量的影響因素及未來糧食需求分析[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)，2015，17（4）：165-175.

[12] 楊建波，王朝暉，邱士可，等. 耕地分等因素對(duì)糧食產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26（4）：306-311.

[13] 李玉浩，王紅葉，張駿達(dá)，等. 華南區(qū)稻田耕地質(zhì)量空間分布與產(chǎn)能提升潛力[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)（中英文），2023，31（10）：1613-1625.

[14] 孫宏勇，張雪佳，田柳，等. 咸水灌溉影響耕地質(zhì)量和作物生產(chǎn)的研究進(jìn)展[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)（中英文），2023，31（3）：354-363.

[15] 李青松，鄧素君，徐國(guó)勁，等. 河南省糧食產(chǎn)量波動(dòng)特征及影響因素分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2015，31（18）：226-230.

[16] 劉忠，黃峰，李保國(guó). 2003—2011年中國(guó)糧食增產(chǎn)的貢獻(xiàn)因素分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（23）：1-8.

[17] 朱晶，李天祥，林大燕，等. “九連增” 后的思考：糧食內(nèi)部結(jié)構(gòu)調(diào)整的貢獻(xiàn)及未來潛力分析[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題，2013，34（11）：36-43，110-111.

[18] 李東升，邱宇婷，潘春燕. 基于群組結(jié)構(gòu)算法的貴州省糧食產(chǎn)量影響因素分析[J]. 河南科技，2023，42（18）：19-22.

[19] 李靜，李晶瑜. 中國(guó)糧食生產(chǎn)的化肥利用效率及決定因素研究[J]. 農(nóng)業(yè)現(xiàn)代化研究，2011，32（5）：565-568.

[20] 張利庠，彭輝，靳興初. 不同階段化肥施用量對(duì)我國(guó)糧食產(chǎn)量的影響分析：基于1952—2006年30個(gè)省份的面板數(shù)據(jù)[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2008（4）：85-94.

[21] 金京淑，劉妍. 吉林省糧食單產(chǎn)影響因素分析[J]. 吉林農(nóng)業(yè)科學(xué)，2010，35（3）：57-59，64.

[22] 王美青，衛(wèi)新，胡豹，等. 浙江省糧食單產(chǎn)影響因素分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2006，22（8）：617-620.

[23] 聶振邦，劉韌，王正友，等. 世界糧食供求現(xiàn)狀、趨勢(shì)和對(duì)策研究[J]. 中國(guó)稻米，2004，10（5）：1-5.

[24] 趙海英，張明旭. 化肥施用量對(duì)中國(guó)糧食產(chǎn)量的影響分析[J]. 經(jīng)濟(jì)研究導(dǎo)刊，2013（35）：71-72.

[25] 劉芳，劉帥. 化肥施用量對(duì)中國(guó)糧食產(chǎn)量的影響分析[J]. 中國(guó)農(nóng)機(jī)化學(xué)報(bào)，2021，42（8）：92-100.

（責(zé)編：張 蓓）