何清 明輝

（四川農(nóng)業(yè)大學(xué)經(jīng)濟學(xué)院，四川 成都 611130）

一、前言

“倉廩實、天下安”，糧食是國家的根本，是經(jīng)濟社會發(fā)展的基礎(chǔ)，糧食安全關(guān)系到一個國家的命脈。作為世界上人口最多的國家，我國以不到全球7%的土地養(yǎng)活了20%的世界人口。能否保障糧食安全不僅對于我國，甚至對于世界糧食安全也有著重大影響。糧食產(chǎn)量安全是糧食安全的核心基礎(chǔ)和重要內(nèi)容。2003－2015年，我國糧食實現(xiàn)十二連增，由2003年的43069.5萬噸增長到2015年的62143.9萬噸，增長率達到 44.29%，這體現(xiàn)出我國糧食產(chǎn)量的整體提升。然而，隨著我國城鎮(zhèn)化和工業(yè)化的推進，一方面，飼料糧和工業(yè)糧需求明顯增加，糧食供給壓力不斷加大；另一方面糧食生產(chǎn)受到耕地、水資源的約束日益顯著，糧食生產(chǎn)要素及政策環(huán)境等的變化也對我國糧食安全形成潛在風(fēng)險。糧食生產(chǎn)產(chǎn)區(qū)之間由于自然環(huán)境、技術(shù)等因素對彼此產(chǎn)量存在一定的影響作用。面對新形勢，為確保我國糧食安全，應(yīng)穩(wěn)步提高糧食產(chǎn)量［1］，重新審視我國糧食安全形勢并調(diào)整和完善相關(guān)政策［2］。

二、文獻回顧

能否保障糧食產(chǎn)量安全，直接影響到人民生活和社會穩(wěn)定。由于糧食生產(chǎn)是自然再生產(chǎn)和經(jīng)濟再生產(chǎn)結(jié)合的過程，所以糧食產(chǎn)量高低是多方面因素綜合作用的結(jié)果。為了有針對性地制定相關(guān)糧食政策，有必要對糧食產(chǎn)量的主要影響因素進行分析［3］。

耕地是不可復(fù)制的農(nóng)業(yè)自然資源，是糧食生產(chǎn)最重要的物質(zhì)基礎(chǔ)，是約束糧食生產(chǎn)的本質(zhì)要素［4］?，F(xiàn)階段，全國耕地面積為 18.2574 億畝，但人均耕地面積不足1.35畝，不足世界平均水平的40%，這對于保障糧食產(chǎn)量是一個巨大的挑戰(zhàn)［5］。耕地面積是保障糧食播種面積的前提，糧食播種面積對糧食產(chǎn)量產(chǎn)生直接影響。很多學(xué)者的研究都表明糧食播種面積對糧食產(chǎn)量起到主要影響作用［6－8］。關(guān)于城鎮(zhèn)化和工業(yè)化引起的農(nóng)業(yè)勞動力轉(zhuǎn)移對糧食生產(chǎn)的影響，則存在兩種不同的觀點。一種觀點認為農(nóng)業(yè)勞動力向非農(nóng)產(chǎn)業(yè)轉(zhuǎn)移，可以促進農(nóng)業(yè)生產(chǎn)中資本、技術(shù)等的投入，進而提高生產(chǎn)效率，增加糧食產(chǎn)量［9］；而另一種觀點則認為我國農(nóng)業(yè)生產(chǎn)的“劉易斯轉(zhuǎn)折點”已經(jīng)到來［10］，農(nóng)業(yè)“內(nèi)卷化”現(xiàn)象不復(fù)存在，勞動力轉(zhuǎn)移會造成糧食生產(chǎn)的勞動力要素投入不足，影響糧食產(chǎn)量安全［11］。氣候變化對糧食產(chǎn)量的影響既有線性的，也有非線性的［12］。尹朝靜等［13］認為降水增加和氣溫上升均對糧食產(chǎn)量變動產(chǎn)生了正向影響。周曙東等［14］的研究結(jié)果顯示未來氣溫上升導(dǎo)致的糧食單產(chǎn)下降將對中國今后的糧食生產(chǎn)帶來一定程度的不利影響。

從長遠來看，技術(shù)進步對糧食產(chǎn)量的影響程度在不斷加深［15］，糧食單產(chǎn)對糧食總產(chǎn)量的貢獻率正在逐年上升。而化肥、農(nóng)藥使用量、機械化發(fā)展水平、有效灌溉面積等因素共同決定了糧食單產(chǎn)的高低［16］。化肥和農(nóng)藥在糧食增產(chǎn)中發(fā)揮了至關(guān)重要的作用，但是目前我國化肥和農(nóng)藥的利用效率很低，使用量還存在很大的節(jié)約空間［17］［18］。同時，化肥和農(nóng)藥的過度使用也讓糧食生產(chǎn)的生態(tài)環(huán)境壓力不斷增加［19］。農(nóng)業(yè)機械化水平作為技術(shù)進步的一個表現(xiàn)，農(nóng)業(yè)機械的普及和推廣，有利于解放農(nóng)業(yè)勞動力和提高糧食生產(chǎn)效率［20］。在推動農(nóng)業(yè)發(fā)展過程中，中央和地方政府通過各項農(nóng)業(yè)政策以及財政資金向農(nóng)業(yè)傾斜的方式，對糧食生產(chǎn)產(chǎn)生了顯著的正向影響［21］。有學(xué)者認為財政支農(nóng)資金應(yīng)該向糧食主產(chǎn)區(qū)傾斜，以求在最小化財政支農(nóng)支出總額的前提下，最大程度地發(fā)揮對農(nóng)業(yè)生產(chǎn)的積極作用［22］。

目前學(xué)術(shù)界已經(jīng)對糧食產(chǎn)量安全、糧食產(chǎn)量的主要影響因素都進行了比較深入和全面的研究，但是在研究方法上，已有文獻大多通過動態(tài)面板的GMM方法對糧食產(chǎn)量的影響因素進行分析，而忽略了不同地區(qū)糧食生產(chǎn)之間存在的空間互動效應(yīng)。由于糧食生產(chǎn)相鄰地區(qū)的資源環(huán)境存在相似性，經(jīng)濟活動存在互動性，所以糧食生產(chǎn)存在一定的空間溢出效應(yīng)，一個地區(qū)的糧食產(chǎn)量不僅受到自身影響因素的作用，還受到鄰近地區(qū)各影響因素的作用，所以將各研究主體當作空間獨立同分布是不合實際的，應(yīng)該將“空間數(shù)據(jù)”引入到實際分析中。將糧食生產(chǎn)的空間相關(guān)性納入研究范圍，才能更加準確地反映糧食生產(chǎn)各影響因素的具體作用大小。因此，研究我國糧食生產(chǎn)的空間效應(yīng)及糧食產(chǎn)量的影響因素，并提出針對性的政策建議具有很強的現(xiàn)實意義。本文基于2003－2014年我國糧食生產(chǎn)的省級面板數(shù)據(jù)，運用空間計量方法，從空間經(jīng)濟學(xué)視角對糧食產(chǎn)量的影響因素進行分析，并在此基礎(chǔ)上提出相關(guān)政策建議。

三、實證分析

考慮到西藏地區(qū)的數(shù)據(jù)缺失，本文基于2003－2014年我國大陸除西藏以外的30個?。ㄖ陛犑小⒆灾螀^(qū)）的面板數(shù)據(jù)進行分析。從空間經(jīng)濟學(xué)視角對糧食產(chǎn)量影響因素進行分析，其前提是各省區(qū)的糧食生產(chǎn)存在空間相關(guān)性，所以在空間計量分析之前需要通過空間相關(guān)性檢驗，確定各省之間糧食生產(chǎn)是否存在空間互動效應(yīng)，在存在空間相關(guān)性的情況下，再引入空間權(quán)重矩陣對各影響因素的作用進行研究。因此，本文的實證分析分為兩個部分，第一部分是使用全域Moran’s I指數(shù)對我國糧食生產(chǎn)進行空間相關(guān)性檢驗，以確定糧食生產(chǎn)是否存在區(qū)域聚集；第二部分構(gòu)建空間杜賓模型，通過對相關(guān)變量的估計，實證分析影響我國糧食產(chǎn)量的因素及其作用大小。

（一）糧食產(chǎn)量的空間相關(guān)性檢驗

根據(jù)地理學(xué)第一定律，包括糧食生產(chǎn)在內(nèi)的許多經(jīng)濟數(shù)據(jù)都涉及空間位置，不同省份的相關(guān)變量不是相互獨立的，而是存在一定的互動，并且距離越近的省份，經(jīng)濟數(shù)據(jù)的聯(lián)系就越密切，所以在研究糧食生產(chǎn)時要加入“空間數(shù)據(jù)”。本文采用相鄰空間權(quán)重矩陣來反映各省份之間的“空間數(shù)據(jù)”。空間權(quán)重矩陣的設(shè)置方法如下：

在確定是否使用空間計量方法前，首先考察數(shù)據(jù)是否存在空間依賴性。本文采用全域Moran’s I指數(shù)對空間相關(guān)性進行度量：

其中，為樣本方差，Wij為空間權(quán)重矩陣，為所有空間權(quán)重之和。xi和xj分別表示第i和j地區(qū)的糧食產(chǎn)量；n為地區(qū)總數(shù)。

Moran’s I指數(shù)的取值范圍為－1≤I≤1。若 I＞0表示正相關(guān)性，即高值與高值相鄰、低值與低值相鄰，糧食生產(chǎn)存在空間集聚；若I＜0表示負相關(guān)性，即高值與低值相鄰，糧食生產(chǎn)是分散的，不存在集聚效應(yīng)；若I接近于0，則表明空間分布是隨機的，不存在空間自相關(guān)。

表1是對我國大陸除西藏以外的30個?。ㄖ陛犑?、自治區(qū)）的糧食產(chǎn)量進行空間自相關(guān)檢驗的結(jié)果。2003－2014年的Moran’s I指數(shù)表明，我國不同省份之間的糧食產(chǎn)量存在一定的空間自相關(guān)性，糧食產(chǎn)量的分布存在空間聚集效應(yīng)，即存在糧食高產(chǎn)地區(qū)相互鄰近、糧食低產(chǎn)地區(qū)也相互鄰近的現(xiàn)象。在糧食生產(chǎn)具有空間相關(guān)性的情況下，通過構(gòu)建空間計量模型對影響糧食產(chǎn)量的因素進行分析。

表1 糧食產(chǎn)量空間自相關(guān)全域Moran’s I指數(shù)（2003－2014）

（二）糧食產(chǎn)量影響因素的空間計量分析

1.模型選擇。為了有效度量糧食產(chǎn)量的影響因素及空間效應(yīng)，本文選擇使用空間杜賓模型，模型的具體形式如下：

式中，被解釋變量y為糧食產(chǎn)量，X為糧食產(chǎn)量的影響因素矩陣，W代表空間權(quán)重矩陣，Wy和WX分別表示被解釋變量和解釋變量的空間依賴性。代表N×1階單位矩陣，N為省份的個數(shù)，aln代表常數(shù)項，ε代表誤差項。

2.變量選擇。結(jié)合糧食生產(chǎn)實際及相關(guān)文獻研究，可知糧食產(chǎn)量受到要素投入、資源環(huán)境、政策等多方面因素的共同作用。本文借鑒已有研究，將我國糧食產(chǎn)量的主要影響因素歸納為傳統(tǒng)投入要素、農(nóng)業(yè)生產(chǎn)條件、地理氣候及財政支持力度這四個方面：

傳統(tǒng)投入要素。主要包括糧食作物的播種面積、化肥農(nóng)藥使用量、機械總動力、從事糧食生產(chǎn)的勞動力。統(tǒng)計年鑒中提供的是針對整個農(nóng)業(yè)的總數(shù)據(jù)，本文參照蔣松穎、劉穎等［23］的研究，以各省相應(yīng)年份的糧食作物播種面積占農(nóng)作物總播種面積的比重作為權(quán)數(shù)A，用該權(quán)數(shù)A乘以相應(yīng)的農(nóng)業(yè)總物質(zhì)投入，算出糧食生產(chǎn)的物質(zhì)投入。

農(nóng)業(yè)生產(chǎn)條件。如糧食生產(chǎn)的大型機械設(shè)施等。以農(nóng)村住戶固定資產(chǎn)投資完成額進行衡量。

地理氣候。以有效灌溉面積和年平均溫度進行衡量。

糧食財政支持力度。以財政農(nóng)林水各項支出乘以權(quán)數(shù)A進行衡量。

3.數(shù)據(jù)來源。本文的變量數(shù)據(jù)主要來源于2003－2014年《全國農(nóng)產(chǎn)品成本收益資料匯編》《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》及各?。ㄊ校┑慕y(tǒng)計年鑒、農(nóng)村統(tǒng)計年鑒，個別缺失數(shù)據(jù)采用線性插值法得到，各變量的描述性統(tǒng)計如表2所示。效應(yīng)是某一解釋變量對本地區(qū)糧食產(chǎn)量的平均影響；間接效應(yīng)是某一解釋變量對其他地區(qū)糧食產(chǎn)量的平均影響；總效應(yīng)是某一解釋變量對所有地區(qū)糧食產(chǎn)量造成的平均影響。下面分別對空間杜賓模型估計結(jié)果的直接效應(yīng)、間接效應(yīng)和總效應(yīng)進行分析。

由SDM模型直接效應(yīng)的結(jié)果來看（見表4），糧食作物播種面積、糧食作物的化肥施用量對本地區(qū)糧食產(chǎn)量具有顯著的正向影響；溫度對糧食產(chǎn)量具有顯著的消極影響，這一結(jié)果的原因主要是目前我國糧食主產(chǎn)區(qū)已經(jīng)由溫度相對適宜的東部地區(qū)向西部和北部地區(qū)轉(zhuǎn)移，所以盡管東部

表2 各變量的描述性統(tǒng)計

4.估計結(jié)果分析。為了更加直觀地觀察各影響因素對糧食產(chǎn)量的影響程度，將各變量都取對數(shù)后進行參數(shù)估計，Husman檢驗拒絕固定效應(yīng)模型與隨機效應(yīng)模型之間無差別的零假設(shè)，證明固定效應(yīng)模型更加適用，空間杜賓模型的估計結(jié)果見表3。

由表3可知，空間自回歸系數(shù)ρ（rho）在5%水平上顯著為正，表明我國糧食生產(chǎn)確實存在一定的空間集聚效應(yīng)，該結(jié)果與Moran’s I指數(shù)的檢驗結(jié)果相一致。由于不同省份的糧食生產(chǎn)存在空間溢出效應(yīng)，所以空間杜賓模型的估計系數(shù)不能夠直接體現(xiàn)各解釋變量對糧食產(chǎn)量的影響。估計結(jié)果中的直接效應(yīng)、間接效應(yīng)和總效應(yīng)則能夠很好地解釋自變量對因變量的影響。其中，直接地區(qū)糧食溫度適宜，但是其糧食產(chǎn)量逐年降低，而溫度不是最適合糧食生產(chǎn)的北部和西部地區(qū)的糧食產(chǎn)量卻在逐年上升；勞動力要素對糧食生產(chǎn)具有正向影響，但是并不十分顯著，這在一定程度上說明我國農(nóng)業(yè)內(nèi)卷化問題得到了改善。

表3 空間杜賓模型的估計結(jié)果

表4 各變量的直接效應(yīng)

表5 各變量的間接效應(yīng)

由SDM模型的間接效應(yīng)結(jié)果（見表5），可以發(fā)現(xiàn)，農(nóng)業(yè)機械總動力對其他地區(qū)的糧食生產(chǎn)具有較為顯著的正向影響，即一個地區(qū)的機械化水平的提高會對其他地區(qū)的糧食生產(chǎn)產(chǎn)生積極影響，這主要是由于一個地區(qū)機械化水平的提高，會通過技術(shù)溢出效應(yīng)對鄰近地區(qū)的機械化水平產(chǎn)生一個正向促進作用，進而鄰近地區(qū)的糧食產(chǎn)量也會得到提升；糧食作物播種面積和農(nóng)藥使用量在10%的顯著水平上對其他地區(qū)的糧食生產(chǎn)具有消極影響，原因可能是本地區(qū)糧食播種面積擴大所新增的產(chǎn)量可能擠占其他地區(qū)糧食的市場份額，從而影響其他地區(qū)的糧食產(chǎn)量；而本地區(qū)使用的農(nóng)藥會通過水和土壤對其他地區(qū)產(chǎn)生負的外部性。

從SDM模型的總效應(yīng)來看（見表6），糧食作物化肥施用量、糧食作物機械總動力對糧食產(chǎn)量具有正向影響；農(nóng)藥使用量對糧食產(chǎn)量具有消極影響，其中糧食作物的化肥施用量是影響糧食產(chǎn)量的最主要因素，這與我國目前的糧食生產(chǎn)現(xiàn)實情況是吻合的。

表6 各變量的總效應(yīng)

四、結(jié)論與建議

（一）結(jié)論

全域Moran’s I指數(shù)顯示我國糧食生產(chǎn)存在正的空間相關(guān)性，鄰近省份的糧食生產(chǎn)存在空間溢出，糧食高產(chǎn)區(qū)與糧食高產(chǎn)區(qū)相鄰、糧食低產(chǎn)區(qū)與糧食低產(chǎn)區(qū)相鄰，糧食生產(chǎn)布局存在空間聚集效應(yīng)?？臻g杜賓模型的分析結(jié)果顯示，糧食作物播種面積、化肥施用量對本地區(qū)糧食產(chǎn)量具有顯著的正向影響；農(nóng)業(yè)機械總動力對其他地區(qū)的糧食生產(chǎn)具有較為顯著的正向影響；糧食作物的化肥施用量是影響其產(chǎn)量的最主要因素。

（二）政策建議

1.制定差異化的糧食生產(chǎn)扶持政策。考慮到糧食生產(chǎn)的空間相關(guān)性，應(yīng)根據(jù)不同糧食生產(chǎn)集聚區(qū)制定不同的糧食生產(chǎn)扶持政策。針對糧食高產(chǎn)集聚區(qū)，制定相關(guān)激勵政策，提高農(nóng)民的種糧積極性，鼓勵農(nóng)民擴大糧食種植規(guī)模，提高糧食單產(chǎn)，確保糧食主產(chǎn)區(qū)穩(wěn)定增產(chǎn)，為國家糧食安全做好根本保障作用；針對糧食低產(chǎn)集中區(qū)，鼓勵農(nóng)民根據(jù)當?shù)貙嶋H環(huán)境種植糧食，盡量能夠滿足本地區(qū)口糧自給。

對化肥產(chǎn)業(yè)進行重點扶持。空間杜賓模型的分析結(jié)果顯示，化肥施用量是影響我國糧食產(chǎn)量的最重要因素，所以為了保障我國糧食安全，提高糧食的綜合生產(chǎn)能力，需要國家對化肥產(chǎn)業(yè)進行大力扶持，對化肥行業(yè)進行宏觀調(diào)控。一方面，可以通過減免稅和生產(chǎn)補貼等方式減輕化肥企業(yè)的資金壓力，進而降低化肥的銷售價格，最終降低農(nóng)民種糧的成本支出，提高農(nóng)民的種糧積極性；另一方面，可以通過向化肥企業(yè)提供最新的技術(shù)和生產(chǎn)設(shè)備，在提高化肥能效的同時，盡量減少其對環(huán)境可能造成的污染和損害。

對農(nóng)藥產(chǎn)業(yè)進行嚴格把控，加大對農(nóng)民的環(huán)保意識宣傳。農(nóng)藥殘留一直以來都是我國農(nóng)產(chǎn)品出口受阻的一個重要原因，這其中既有綠色壁壘的影響，也存在我國農(nóng)產(chǎn)品自身質(zhì)量不過硬的問題。根據(jù)相關(guān)數(shù)據(jù)顯示，我國農(nóng)藥使用量排世界第一，每年使用的農(nóng)藥數(shù)量超過80萬噸［24］，單位面積的農(nóng)藥使用量是世界平均水平的2.5倍［25］。2003－2015年期間，我國農(nóng)藥使用量從132.5萬噸增加到178.3萬噸，累計增加34.57%。所以國家應(yīng)該加強對農(nóng)藥產(chǎn)業(yè)的監(jiān)控，嚴令禁止劇毒農(nóng)藥生產(chǎn)和流通。同時加強對農(nóng)民環(huán)保安全意識的宣傳，改變很多農(nóng)民根深蒂固的觀念，逐步降低糧食生產(chǎn)對農(nóng)藥的依賴性。

2.加大糧食高產(chǎn)集中區(qū)的農(nóng)業(yè)基礎(chǔ)設(shè)施和交通運輸設(shè)施的建設(shè)投入。由空間杜賓模型間接效應(yīng)可知，本地區(qū)農(nóng)業(yè)機械總動力對鄰近地區(qū)的糧食生產(chǎn)具有顯著的正向作用，所以加強糧食高產(chǎn)集中區(qū)的農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)，可以通過空間溢出效應(yīng)，顯著提高區(qū)域糧食產(chǎn)量。交通運輸設(shè)施的完善，可以使糧食高產(chǎn)區(qū)之間聯(lián)系更加緊密，增強糧食高產(chǎn)區(qū)的空間相關(guān)性，進而提高空間溢出效應(yīng)。

3.將“藏糧于地”和“藏糧于人才”戰(zhàn)略結(jié)合。土地是糧食生產(chǎn)的基本保障，要想進一步保障或提高糧食產(chǎn)量，必須繼續(xù)執(zhí)行“十三五”規(guī)劃中“藏糧于地”的戰(zhàn)略要求，確保18億畝耕地紅線，同時積極推進中低產(chǎn)田的改良，積極推進高標準農(nóng)田建設(shè)。有農(nóng)業(yè)生產(chǎn)專業(yè)知識的人才是目前推動糧食產(chǎn)業(yè)持續(xù)穩(wěn)定發(fā)展的一大動力，推進“藏糧于人才”計劃，國家應(yīng)該大力發(fā)展糧食生產(chǎn)人才培養(yǎng)計劃，盡快培養(yǎng)出一批高素質(zhì)的職業(yè)農(nóng)民和職業(yè)經(jīng)理人。

［1］陳飛，范慶泉，高鐵梅.農(nóng)業(yè)政策、糧食產(chǎn)量與糧食生產(chǎn)調(diào)整能力［J］.經(jīng)濟研究，2010，45（11）：101－114，140.

［2］高鐵生.我國糧食安全形勢與政策建議［J］.經(jīng)濟研究參考，2013，（28）：17－39.

［3］葉妍君，齊清文，姜莉莉，張岸.基于地理探測器的黑龍江墾區(qū)農(nóng)場糧食產(chǎn)量影響因素分析［J］.地理研究，2018，37（01）：171－182.

［4］聶英.中國糧食安全的耕地貢獻分析［J］.經(jīng)濟學(xué)家，2015（01）：83－93.）

［5］王曉航，鐘帆，孟凡珂.保障我國糧食安全的肥料問題探究［J］.農(nóng)村經(jīng)濟與科技，2018，29（04）：148－149

［6］孫通，封志明，楊艷昭.2003－2013 年中國縣域單元糧食增產(chǎn)格局及貢獻因素研究［J］.自然資源學(xué)報，2017，32（02）：177－185.

［7］馮超.中國谷物產(chǎn)出的“面積—質(zhì)量”導(dǎo)向因素分析［J］.干旱區(qū)資源與環(huán)境，2015，29（08）：7－13.

［8］劉守義.我國糧食主產(chǎn)區(qū)糧食產(chǎn)量波動及增長影響因素分析［J］.江西社會科學(xué)，2014，34（08）：96－100.

［9］李谷成，馮中朝，范麗霞.小農(nóng)戶真的更加具有效率嗎？——來自湖北省的經(jīng)驗證據(jù)［J］.經(jīng)濟學(xué)（季刊），2010，9（01）：95－124.

［10］蓋慶恩，朱喜，史清華.勞動力轉(zhuǎn)移對中國農(nóng)業(yè)生產(chǎn)的影響［J］.經(jīng)濟學(xué)（季刊），2014，13（03）：1147－1170.

［11］賈利軍，楊靜.農(nóng)村勞動力弱化背景下我國糧食安全保障研究［J］.中州學(xué)刊，2015，（02）：37－44.

［12］何為，劉昌義，劉杰，陳夢玫.氣候變化和適應(yīng)對中國糧食產(chǎn)量的影響——基于省級面板模型的實證研究［J］.中國人口·資源與環(huán)境，2015，25（S2）：248－253.

［13］尹朝靜，李谷成，高雪.氣候變化對中國糧食產(chǎn)量的影響——基于省級面板數(shù)據(jù)的實證［J］.干旱區(qū)資源與環(huán)境，2016，30（06）：89－94.

［14］周曙東，周文魁，林光華，喬輝.未來氣候變化對我國糧食安全的影響［J］.南京農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），2013，13（01）：56－65.

［15］曾志勇，劉穎.長江中游地區(qū)耕地資源數(shù)量與糧食生產(chǎn)關(guān)系的實證分析［J］.統(tǒng)計與決策，2016（12）：102－105.

［16］劉守義.我國糧食主產(chǎn)區(qū)糧食產(chǎn)量波動及增長影響因素分析［J］.江西社會科學(xué)，2014，34（08）：96－100.

［17］劉德偉，李強，宋孝航.中國糧食生產(chǎn)化肥施用效率分析——基于隨機前沿生產(chǎn)函數(shù)［J］.資源開發(fā)與市場，2017，33（04）：401－407.

［18］牛亮云，侯博，吳林海.基于灰關(guān)聯(lián)熵的中國農(nóng)業(yè)能源投入與糧食產(chǎn)出關(guān)系研究［J］.財貿(mào)研究，2012，23（02）：45－53.

［19］王一杰，邸菲，辛嶺.我國糧食主產(chǎn)區(qū)糧食生產(chǎn)現(xiàn)狀、存在問題及政策建議［J］.農(nóng)業(yè)現(xiàn)代化研究，2018，39（01）：37－47.

［20］王珺鑫，楊學(xué)成.山東省糧食生產(chǎn)波動及主要投入要素效應(yīng)的實證分析——基于17地市的面板數(shù)據(jù)［J］.中國農(nóng)業(yè)資源與區(qū)劃，2015，36（03）：18－23.

［21］成德寧，李燕.2003－2012年中國農(nóng)業(yè)增長的源泉——基于制度的分析［J］.福建論壇（人文社會科學(xué)版），2017（07）：68－74.

［22］韓君.糧食安全與財政支農(nóng)支出結(jié)構(gòu)調(diào)整［J］.商業(yè)研究，2013（12）：138－141.

［23］周力，周應(yīng)恒.糧食安全：氣候變化與糧食產(chǎn)地轉(zhuǎn)移［J］.中國人口·資源與環(huán)境，2011，21（07）：162－168.

［24］李玉輝.綠色壁壘對我國蔬菜出口的影響及對策［J］.經(jīng)貿(mào)實踐，2017，（23）：39－40.

［25］李翔.勞動力轉(zhuǎn)移對河南省糧食生產(chǎn)技術(shù)效率的影響［D］.華中農(nóng)業(yè)大學(xué)，2016.