于元赫 吳 健* 李子君 胡 蕾

(1.中國人民大學(xué) 環(huán)境學(xué)院，北京 100872； 2.山東師范大學(xué) 地理與環(huán)境學(xué)院，濟南 250014)

糧食安全是人類共同關(guān)注的全球性重大問題之一，是一個國家和地區(qū)經(jīng)濟發(fā)展和社會穩(wěn)定的重要基礎(chǔ)[1]。黨的十九大報告指出，“確保國家糧食安全，把中國人的飯碗牢牢端在自己手中。”[2]國家糧食安全的重要前提是保障區(qū)域糧食供需平衡，但由于不同區(qū)域農(nóng)業(yè)資源稟賦和農(nóng)業(yè)生產(chǎn)條件存在明顯差異，其糧食生產(chǎn)水平和糧食安全狀況也不同[3]。因此，深入探討區(qū)域糧食生產(chǎn)變化及時空分布格局，對于優(yōu)化糧食生產(chǎn)布局、保障國家糧食安全具有重要的實踐意義。

我國耕地資源短缺，人口眾多，糧食供需矛盾一直存在[4]。在這種背景下，大量研究關(guān)注了糧食生產(chǎn)的時空格局[5-8]、區(qū)域差異[9-12]、影響因素[13-16]、糧食安全[17-21]等相關(guān)領(lǐng)域。在全國尺度上，Liu等[22]采用GAEZ模型分析了中國耕地糧食生產(chǎn)潛力空間格局，指出中國耕地糧食生產(chǎn)潛力呈現(xiàn)南減北增、總量減少的特征；陳秧分等[23]對中國糧食產(chǎn)量變化的時空格局進行了探究，發(fā)現(xiàn)中國糧食產(chǎn)量變化可分為3個階段，東中西部省份糧食生產(chǎn)存在較大差異。在區(qū)域尺度上，劉玉等[24]運用LMDI模型定量分析了黃淮海地區(qū)糧食生產(chǎn)格局及影響因素，指出糧食增產(chǎn)主要來源于糧食單產(chǎn)的提高；潘佩佩等[25]揭示了太湖流域糧食生產(chǎn)時空格局及糧食安全水平，發(fā)現(xiàn)該流域存在較大的糧食安全風(fēng)險。這些成果較好地揭示了國家、地區(qū)和省、市等層面糧食生產(chǎn)格局的演化規(guī)律，為客觀認識我國及區(qū)域糧食生產(chǎn)在宏觀格局上的變化提供了依據(jù)，但對于區(qū)域內(nèi)部糧食生產(chǎn)變化的異質(zhì)性及影響機制未及深入討論。

山東省是我國13個傳統(tǒng)糧食生產(chǎn)和消費大省之一，是保障國家糧食安全的重要基礎(chǔ)。然而，隨著城市化和工業(yè)化的推進，大量耕地被侵占，糧食生產(chǎn)面積萎縮；加之人口逐年增多，家庭經(jīng)濟水平不斷提升，糧食需求量大幅提高，糧食生產(chǎn)壓力不斷加大；再加上居民膳食結(jié)構(gòu)和消費格局的轉(zhuǎn)變，導(dǎo)致糧食生產(chǎn)結(jié)構(gòu)也發(fā)生了很大變化。鑒于此，本研究以山東省134個縣(市、區(qū))作為基本分析單元，運用重心遷移模型和探索性空間數(shù)據(jù)分析等方法，探討1995—2015年山東省糧食生產(chǎn)時空格局演變特征，并引入人均糧食占有量、糧食產(chǎn)量波動系數(shù)和變異系數(shù)等指標評價山東省糧食安全水平，以期為山東省糧食可持續(xù)生產(chǎn)及政策制定提供學(xué)術(shù)參考。

1 研究方法與數(shù)據(jù)來源

1.1 研究方法

1.1.1探索性空間數(shù)據(jù)分析(ESDA)

探索性空間數(shù)據(jù)分析(ESDA)是一系列空間數(shù)據(jù)分析方法的集合，它通過將數(shù)據(jù)進行可視化表達來展現(xiàn)事物的空間集聚特征以及與周邊事物的異質(zhì)性。該方法分為全局空間自相關(guān)和局部空間自相關(guān)。本研究采用Global Moran’sI統(tǒng)計量來分析山東省糧食生產(chǎn)空間相關(guān)性的總體趨勢。計算公式為：

(1)

(2)

式中：Z為檢驗統(tǒng)計量；E(I)為期望；Var(I)為方差。

熱點分析G*指數(shù)常用于分析事物在局部空間上的集聚特征，可進一步測度糧食生產(chǎn)熱點區(qū)(高值區(qū))和冷點區(qū)(低值區(qū))的空間分布。計算公式為：

(3)

本研究結(jié)合全局Moran指數(shù)和熱點分析G*指數(shù)來測度山東省糧食生產(chǎn)的空間集聚格局。其中全局Moran指數(shù)反映了糧食生產(chǎn)在整個省域內(nèi)的空間特征，以此來分析城市間空間關(guān)聯(lián)程度；在此基礎(chǔ)上，局部G*指數(shù)可反映糧食生產(chǎn)局部空間異質(zhì)特征，進而判斷局部空間分異規(guī)律。

1.1.2重心遷移模型

糧食生產(chǎn)重心遷移模型能夠反映糧食生產(chǎn)空間分布及遷移規(guī)律，若一個區(qū)域由多個次級區(qū)域組成，第i個次級區(qū)域的重心坐標為(Mi，Ni)，P為該區(qū)域某種屬性值的權(quán)重。計算公式為：

(4)

式中：M、N分別為糧食生產(chǎn)重心的經(jīng)緯度坐標；Pi為第i個縣(市、區(qū))的糧食產(chǎn)量；Mi、Ni分別為第i個縣(市、區(qū))的中心經(jīng)緯度坐標，利用ArcGIS軟件計算得出。

1.1.3糧食產(chǎn)量變異系數(shù)

糧食產(chǎn)量變異系數(shù)為糧食產(chǎn)量標準差與平均值之比，可以反映糧食真實增長率偏離平均值的程度，變異系數(shù)越高，代表糧食生產(chǎn)越不穩(wěn)定。其表達式為：

(5)

1.1.4糧食產(chǎn)量波動指數(shù)

糧食產(chǎn)量波動指數(shù)為糧食總產(chǎn)量年際間的波動幅度，波動指數(shù)越大，代表糧食實際產(chǎn)量對趨勢產(chǎn)量的偏離程度越大，穩(wěn)定性越差。其表達式為：

αt=[|Tt-Ft|/Ft]×100%

(6)

式中：αt為t年的波動指數(shù)；Tt為t年的糧食產(chǎn)量；Ft表示t年的趨勢產(chǎn)量，用糧食產(chǎn)量時間序列計算得到。

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

數(shù)據(jù)來源主要為山東省及各縣(市、區(qū))統(tǒng)計年鑒，包括《山東省統(tǒng)計年鑒》[26]、《中國縣(市) 社會經(jīng)濟統(tǒng)計年鑒》[27]。主要數(shù)據(jù)包括全省及各縣(市)糧食產(chǎn)量、糧食播種面積、主要糧食作物產(chǎn)量與播種面積以及年末常住人口數(shù)。其中各縣(市、區(qū))的糧食產(chǎn)量選取了1995、2005、2010和2015年的4期數(shù)據(jù)。1995年以來山東省部分縣(市、區(qū))行政區(qū)劃發(fā)生調(diào)整，因此以2015年山東省縣域行政區(qū)劃為基準對1995、2005和2010年的縣級單元及相關(guān)統(tǒng)計數(shù)據(jù)進行修正，最終共獲取134個縣(市、區(qū))單元。

2 糧食生產(chǎn)的時序變化

2.1 糧食總產(chǎn)的時序變化

1995—2015年，山東省糧食總產(chǎn)的變化趨勢大體可分為2個階段(圖1)：1995—2002年為波動下降期，糧食總產(chǎn)從4 245.0萬t下降至3 292.7萬t。由于該階段山東省城市化進程加快，建設(shè)用地蠶食耕地，導(dǎo)致復(fù)種指數(shù)明顯下降，加之農(nóng)業(yè)結(jié)構(gòu)調(diào)整，農(nóng)民種糧積極性受挫，導(dǎo)致糧食產(chǎn)量急劇下滑。另外，1997和2002年下降最為嚴重主要是由于這兩年全省發(fā)生了百年不遇的特大干旱，導(dǎo)致糧食大幅減產(chǎn)；2003—2015年為穩(wěn)步上升期，由3 435.5萬t上升至4 712.7萬t，尤其是2004年之后，中央頒布了農(nóng)業(yè)稅減免、直糧補貼、糧食最低收購價等惠農(nóng)舉措，同時，加強了耕地保護，積極推動農(nóng)業(yè)科技進步，一系列重大利好政策使糧食生產(chǎn)得到極大保障和扶持，農(nóng)民種糧積極性顯著提高，糧食連年豐收。

圖1 山東省糧食總產(chǎn)和單產(chǎn)的變化趨勢Fig.1 Changes of grain yield and the per unit area yield of grain in Shandong Province

2.2 糧食單產(chǎn)的時序變化

1995—2015年山東省糧食單產(chǎn)總體呈現(xiàn)出上升趨勢(圖1)，大體可分為2個階段：1995—2002年為波動性變化，期間出現(xiàn)2次明顯的階段性下降(1997年和2002年)；2003—2015年為持續(xù)提高階段，至2015年，糧食單產(chǎn)達到近20多年來最高值，為6 290.0 kg/hm2，相比2002年提升約32個百分點。2003年以后，在農(nóng)作物品種改良、農(nóng)業(yè)生產(chǎn)條件改善和管理技術(shù)進步等因素的綜合作用下，糧食單產(chǎn)能力得到了極大地提高，為近10多年來山東省糧食增產(chǎn)做出了重要貢獻。

2.3 糧食生產(chǎn)結(jié)構(gòu)的時序變化

2.3.1糧食作物產(chǎn)量結(jié)構(gòu)的變化

1995—2015年，在山東省糧食作物產(chǎn)量結(jié)構(gòu)中，小麥和玉米是占糧食作物總產(chǎn)量比重最高的兩種作物(圖2)。其中，小麥的比重高于45%，為山東省“第一主糧”；玉米的比重逐年增加，由1995年的33.9%上升到2015年的43.5%；甘薯作為山東省第三大主糧，比重由1995年的9.6%下降至2015年的3.7%；大豆產(chǎn)量也呈現(xiàn)下降趨勢；而稻谷、谷子、高粱和其他糧食作物產(chǎn)量較低，比重變化不明顯。

圖2 山東省主要糧食作物產(chǎn)量結(jié)構(gòu)Fig.2 The structure of main grain crop yield in Shandong Province

2.3.2糧食作物播種面積結(jié)構(gòu)的變化

1995—2015年，在山東省糧食作物播種面積結(jié)構(gòu)中，小麥和玉米是占糧食作物播種面積比重最高的2種作物(圖3)。其中，小麥的比重最大，接近50%；玉米比重不斷上升，由1995年的29.7%上升到2015年的42.4%，提高了12.7個百分點；甘薯和大豆的播種面積均呈不斷下降的趨勢；而稻谷、谷子和高粱等糧食作物播種面積較小，比重變化不明顯。

圖3 山東省主要糧食作物播種面積結(jié)構(gòu)Fig.3 The structure of sown area of main grain crop in Shandong Province

通過以上分析可以看出，山東省糧食生產(chǎn)結(jié)構(gòu)以小麥和玉米為主，小麥的總產(chǎn)量和播種面積占比最高，玉米次之，而甘薯和大豆占比呈持續(xù)下降的趨勢。究其原因，一方面隨著勞動力務(wù)農(nóng)機會成本上升，進城務(wù)工的農(nóng)民傾向于選擇省時省力的小麥和玉米；另一方面，小麥和玉米優(yōu)良品種的普及以及管理技術(shù)的進步，促進其單產(chǎn)水平持續(xù)提高。然而，甘薯作為山東省第三大主糧，所占比重卻呈現(xiàn)持續(xù)下降趨勢。甘薯與小麥和玉米相比，種植成本低、產(chǎn)量高，但由于市場需求小，往往出現(xiàn)供大于求的情況，種植效益低，導(dǎo)致農(nóng)民種植甘薯的積極性不高，再加之種植技術(shù)落后和粗放管理，因此其總產(chǎn)量和播種面積持續(xù)下滑。

3 糧食生產(chǎn)的空間格局演變

3.1 糧食生產(chǎn)空間格局變化特征

參照相關(guān)文獻[28]，本研究按照山東省縣域糧食產(chǎn)量平均值的50%、100%和150%(即18.81、37.62和56.43萬t為分節(jié)點)，將山東省縣域糧食產(chǎn)量分為低產(chǎn)區(qū)、一般區(qū)、中產(chǎn)區(qū)和高產(chǎn)區(qū)(圖4)。山東省糧食產(chǎn)量存在明顯的區(qū)域差異。在數(shù)量變化上，高產(chǎn)區(qū)的縣(市、區(qū))由1995年的12個增加至2015年的46個；中產(chǎn)區(qū)和一般區(qū)的縣(市、區(qū))呈持續(xù)減少的趨勢；低產(chǎn)區(qū)的縣(市、區(qū))則經(jīng)歷了“增加—減少—增加”的變化。在空間變化上，高產(chǎn)區(qū)的空間范圍明顯擴大，呈現(xiàn)出向魯西地區(qū)集中的趨勢，魯西地區(qū)大部分縣(市、區(qū))，如鄄城縣、商河縣等，已發(fā)展為全省主要糧食高產(chǎn)區(qū)；中產(chǎn)區(qū)和一般區(qū)的空間分布比較分散，波動頻繁；低產(chǎn)區(qū)則主要分布于魯北、魯中山區(qū)和魯東沿海地區(qū)，如墾利縣、沂源縣、牟平區(qū)等，空間分布相對穩(wěn)定。

圖4 山東省糧食生產(chǎn)時空格局變化Fig.4 The spatial-temporal variation of grain production in Shandong Province

3.2 糧食生產(chǎn)空間集聚格局變化特征

基于各縣(市、區(qū))的糧食總產(chǎn)量，利用GeoDa軟件計算得到了1995、2005、2010和2015年的全局Moran 指數(shù)分別為0.205、0.267、0.281、0.337，呈上升趨勢，且Z值均通過1%的置信度檢驗，表明1995—2015年山東省糧食生產(chǎn)的空間分布呈正相關(guān)，整體空間集聚態(tài)勢增強。

為了進一步分析糧食生產(chǎn)的局部空間集聚格局變化特征，通過測度各縣(市、區(qū))糧食產(chǎn)量的局部G*指數(shù)，利用自然斷點分級法將其分成4類，按數(shù)值從高到低依次記為：熱點區(qū)、次熱區(qū)、次冷區(qū)和冷點區(qū)。通過圖5可知，有相當數(shù)量的縣(市、區(qū))糧食產(chǎn)量高值(低值)區(qū)域在空間上集聚，形成高高集聚或低低集聚的空間分異格局。分階段來看：1995—2005年，糧食生產(chǎn)的空間集聚格局波動較大，魯東地區(qū)的熱點區(qū)明顯收縮，而次熱區(qū)與次冷區(qū)有所擴張；2005—2010年，魯西地區(qū)的熱點區(qū)明顯增加，次熱區(qū)、次冷區(qū)和冷點區(qū)變化微弱；2010—2015年，熱點區(qū)繼續(xù)擴大且在山東省西部形成集中連片的分布格局，而冷點區(qū)基本穩(wěn)定，分散分布于魯北、魯中山區(qū)以及魯東沿海地區(qū)?？傮w而言，1995—2015年山東省糧食生產(chǎn)逐漸形成了一定的時間穩(wěn)定性空間分布，魯西地區(qū)形成大規(guī)模穩(wěn)定性熱點，魯北、魯中山區(qū)和魯東沿海地區(qū)形成分散的穩(wěn)定性冷點。

3.3 糧食生產(chǎn)分布重心變化特征

將4個時期的糧食產(chǎn)量重心坐標依次連接，得到糧食生產(chǎn)重心遷移軌跡見圖6?？梢钥闯?，糧食生產(chǎn)重心的地理坐標為117.50°～117.92° E，36.25°～36.28° N，均位于山東省幾何中心(118.15° E，36.36° N)的西南方向，表明山東省糧食生產(chǎn)空間分布不均衡，魯西地區(qū)逐漸成為糧食生產(chǎn)重點區(qū)域。這主要是由于魯西平原地勢平坦、土壤肥沃，適合大規(guī)模機械化作業(yè)，農(nóng)業(yè)優(yōu)勢明顯，僅菏澤市糧食總產(chǎn)量就約占全省的1/8。在當?shù)卣r(nóng)業(yè)政策的扶持下，逐漸發(fā)展為山東省糧食生產(chǎn)的相對優(yōu)勢區(qū)。1995—2015年，山東省糧食生產(chǎn)重心逐漸向西遷移，呈現(xiàn)出糧食生產(chǎn)和經(jīng)濟發(fā)展區(qū)域分工的空間格局。隨著魯東地區(qū)經(jīng)濟和人口不斷增長，糧食生產(chǎn)重心的移動方向與經(jīng)濟和人口重心移動方向相反，將導(dǎo)致山東省糧食供需關(guān)系愈發(fā)緊張。因此，探討山東省糧食生產(chǎn)時空格局演變過程中的糧食安全水平對于協(xié)調(diào)經(jīng)濟發(fā)展與耕地保護之間的關(guān)系，促進糧食可持續(xù)生產(chǎn)及政策制定具有重要意義。

圖5 山東省糧食生產(chǎn)空間集聚格局變化Fig.5 The spatial agglomeration pattern variation of grain production in Shandong Province

圖6 山東省糧食生產(chǎn)重心遷移軌跡Fig.6 The trajectories of the center of gravity for grain production in Shandong Province

4 糧食安全評價

目前，國內(nèi)外學(xué)術(shù)界對于糧食安全評價并沒有統(tǒng)一的標準，研究內(nèi)容比較寬泛。我國對糧食安全的思路正面臨重構(gòu)，其中，保障和提高糧食生產(chǎn)能力是國家糧食安全觀的一個重要體現(xiàn)[28-29]。姚成勝[19,21]、宋小青等[30]和張鳳榮等[31]研究表明，糧食安全的核心就在于確保足夠的、穩(wěn)定的糧食生產(chǎn)能力。因此，保障國家糧食安全的最基本要求就是要穩(wěn)定和提高區(qū)域糧食產(chǎn)量，特別是糧食主產(chǎn)區(qū)的糧食產(chǎn)量。由于不同省(市、自治區(qū))農(nóng)業(yè)資源稟賦存在較大地域差異，其糧食安全的定義和要求也是不同的。對于農(nóng)業(yè)生產(chǎn)條件差，必須依靠調(diào)入糧食滿足需求的省(市、自治區(qū))目前的糧食安全核心是滿足基本的糧食需求，提高糧食生產(chǎn)能力；山東省作為重要的糧食主產(chǎn)省，其對糧食安全的要求更加嚴格，不僅需要具備較強的糧食生產(chǎn)能力，更需要具備較穩(wěn)定的糧食供給能力，能夠維持省內(nèi)外糧食供需動態(tài)平衡并穩(wěn)定糧食市場價格。因此，本文著重從糧食生產(chǎn)能力和糧食供給穩(wěn)定性兩方面評價山東省糧食安全水平。

4.1 糧食生產(chǎn)能力評價

結(jié)合上述分析可知，1995—2015年山東省糧食總產(chǎn)量和單位面積產(chǎn)量得到明顯提升。雖然在一定程度上反映出山東省糧食生產(chǎn)能力較強，但由于每年的糧食播種面積、人口數(shù)量等不斷發(fā)生變化，僅單純從總產(chǎn)量和單位面積產(chǎn)量方面考慮糧食生產(chǎn)能力缺乏一定的科學(xué)性和完整性[32]。人均糧食占有量既是反映糧食生產(chǎn)能力的綜合指標，也是衡量糧食安全程度的評價指標[19,24]，人均糧食占有量越高，說明糧食生產(chǎn)能力越高、糧食安全程度也越好；反之亦然。

不同統(tǒng)計數(shù)據(jù)口徑下的人口數(shù)據(jù)會對人均糧食占有量的計算結(jié)果產(chǎn)生重要影響，這種影響在人口數(shù)量大、人口流動活躍的山東省更加突出。參考相關(guān)文獻[33]，本研究在計算人均糧食占有量時采用年末常住人口數(shù)據(jù)。綜合張利國[20]、姚成勝等[21]的觀點，將人均糧食占有量表征的糧食安全類型大致分成以下6種：<300 kg為嚴重不安全型，≥300～350 kg為較不安全型，≥350～400 kg為溫飽型，≥400～450 kg為較安全型，≥450～500 kg為小康型，≥500 kg為富裕型，且定義人均糧食占有量小于350 kg為糧食不安全。

1995—2015年山東省人均糧食占有量呈先下降后上升的變化趨勢(圖7)，糧食安全狀況整體趨于好轉(zhuǎn)，由溫飽型過渡到較安全型，又重新邁入小康型。1995年山東省人均糧食占有量為488 kg，處于小康型。此后人均糧食占有量呈波動下降趨勢，2002—2004年均跌破400 kg，糧食安全狀況惡化，下降為溫飽型。2005年開始，糧食安全狀況出現(xiàn)持續(xù)好轉(zhuǎn)，人均糧食占有量再次突破400 kg，從溫飽型過渡到較安全型，并持續(xù)到2007年。這期間山東省人均糧食占有量雖然有所波動，但一直未能突破450 kg，始終處于較安全型。2008年人均糧食占有量再次超過450 kg，重新回到小康型，并一直保持上升趨勢。這期間山東省人均糧食占有量基本穩(wěn)定在450—500 kg之間，糧食安全狀況明顯改善。

圖7 山東省人均糧食占有量的變化Fig.7 The variation of the per capita share of grain in Shandong Province

4.2 糧食供給穩(wěn)定性評價

一個區(qū)域不同時期糧食總產(chǎn)量波動幅度的大小以及糧食真實增長率偏離平均值的程度反映了糧食供給的平穩(wěn)程度，也體現(xiàn)了其糧食安全風(fēng)險水平。為此，本研究結(jié)合糧食產(chǎn)量波動系數(shù)和糧食產(chǎn)量變異系數(shù)兩個指標來表征山東省糧食供給穩(wěn)定性。

通過計算可知(圖8)，1995—2015年，山東省糧食產(chǎn)量波動系數(shù)的范圍為0.35%～18.24%，其中有近一半以上年份的波動系數(shù)高于同時期全國平均水平，存在糧食安全風(fēng)險。尤其是2002—2004年的波動系數(shù)均超過14%，這一階段受農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和自然災(zāi)害的雙重影響，糧、棉、蔬果等農(nóng)產(chǎn)品產(chǎn)量均大幅減產(chǎn)，從而導(dǎo)致糧食作物產(chǎn)量降至近20年最低值(圖1)。可見，波動系數(shù)較大的年份與減產(chǎn)幅度較大的年份(2002—2004年)在時間上具有很好的對應(yīng)性，表明糧食減產(chǎn)越多，糧食安全風(fēng)險越大。但2004年以后波動系數(shù)明顯下降，均未超過6%，僅有2010年和2011年略高于全國平均水平，表明該時期糧食安全風(fēng)險逐漸降低。

圖8 山東省和全國不同時期糧食產(chǎn)量波動系數(shù)比較Fig.8 Comparison of fluctuation coefficients of grain yield in different periods in Shandong Province and China

通過計算可知，1995—2015年山東省糧食產(chǎn)量變異系數(shù)為9.28%，高于全國平均水平(5%)，表明山東省糧食供給穩(wěn)定性較差。目前山東省二、三產(chǎn)業(yè)發(fā)展迅速，經(jīng)濟的快速發(fā)展也帶來耕地非農(nóng)化速度加快、耕地質(zhì)量下降等問題。雖然2004年以來糧食總產(chǎn)量有所增長，但同期糧食總產(chǎn)量占全國百分比由8.98%下降至7.58%。利用Pearson相關(guān)系數(shù)進行相關(guān)性分析發(fā)現(xiàn)，山東省糧食總產(chǎn)量占全國百分比與同期GDP占全國百分比在5%的顯著性水平下呈負相關(guān)，相關(guān)系數(shù)為-0.387。由此可見，在未來相當長的一段時間內(nèi)耕地面積減少的趨勢難以逆轉(zhuǎn)的態(tài)勢下，糧食安全問題值得高度關(guān)注。

5 結(jié)論與討論

5.1 結(jié)論

本研究基于重心遷移模型、探索性空間數(shù)據(jù)分析等方法，探討了山東省糧食生產(chǎn)變化及時空分布格局，并引入人均糧食占有量、糧食產(chǎn)量波動系數(shù)和變異系數(shù)等指標評價了山東省糧食安全狀況，主要結(jié)論如下：

1)1995—2015年，山東省糧食總產(chǎn)經(jīng)歷了從波動下降到穩(wěn)步上升2個階段，而糧食單產(chǎn)整體呈上升趨勢，糧食單產(chǎn)水平的提高為糧食增產(chǎn)做出了重要貢獻；在主要糧食作物中，小麥的總產(chǎn)量和播種面積占比最高，玉米次之，而甘薯和大豆占比呈持續(xù)下降的趨勢。

2)山東省糧食生產(chǎn)存在明顯的空間差異，魯西地區(qū)大部分縣(市、區(qū))已發(fā)展為全省主要糧食高產(chǎn)區(qū)，糧食生產(chǎn)重心向西遷移；糧食生產(chǎn)空間集聚態(tài)勢增強，呈現(xiàn)出高高集聚或低低集聚特征，魯西地區(qū)成為大規(guī)模穩(wěn)定性熱點，魯北、魯中山區(qū)以及魯東沿海地區(qū)成為穩(wěn)定性冷點。

3)1995—2015年，山東省糧食生產(chǎn)能力逐漸提高，糧食安全狀況整體趨于好轉(zhuǎn)，由溫飽型過渡到較安全型，又重新邁入小康型；然而糧食產(chǎn)量變異系數(shù)和近一半以上年份的糧食產(chǎn)量波動系數(shù)均超全國平均水平，表明山東省糧食供給穩(wěn)定性較差，糧食安全問題依然嚴峻。

5.2 討論

糧食生產(chǎn)空間分布格局和集聚格局表明，山東省糧食生產(chǎn)存在明顯的空間差異，這可歸因于不同區(qū)域的自然資源狀況、城鎮(zhèn)化水平和經(jīng)濟發(fā)展水平不同。魯西地區(qū)農(nóng)業(yè)生產(chǎn)條件優(yōu)越，種糧優(yōu)勢較大，而魯東地區(qū)經(jīng)濟較發(fā)達，城鎮(zhèn)化水平較高。李明杰等[34]研究發(fā)現(xiàn)，魯西北地區(qū)光熱水土資源較好，農(nóng)業(yè)基礎(chǔ)設(shè)施較為完善，糧食生產(chǎn)能力較強。王兆華等[35]也進一步指出山東省糧食生產(chǎn)應(yīng)向種糧自然條件相對較好的地區(qū)進行轉(zhuǎn)移，可有效降低糧食生產(chǎn)成本。因此，未來應(yīng)盡快建立區(qū)域糧食生產(chǎn)補償制度，加大對糧食生產(chǎn)大縣(市、區(qū))的補償力度(如定陶縣、鄄城縣、商河縣等)，建立區(qū)域糧食生產(chǎn)利益補償?shù)拈L效機制(如鼓勵魯東地區(qū)經(jīng)濟較發(fā)達的城市到魯西地區(qū)投資建設(shè)糧食生產(chǎn)基地等)，實現(xiàn)糧食主產(chǎn)地區(qū)與主銷地區(qū)之間的利益平衡；因地制宜、分類指導(dǎo)省內(nèi)不同地區(qū)糧食生產(chǎn)，對于糧食生產(chǎn)優(yōu)勢區(qū)(如聊城、菏澤等地區(qū))，應(yīng)加大種糧補貼和農(nóng)機具購置補貼力度，快速推動糧食生產(chǎn)向機械化、規(guī)模化、專業(yè)化方向發(fā)展，同時，加快推進高標準基本農(nóng)田建設(shè)，力爭將魯西地區(qū)建成“山東糧倉”；對于糧食生產(chǎn)非優(yōu)勢區(qū)(如煙臺、威海等地區(qū))，應(yīng)協(xié)調(diào)好城市化、經(jīng)濟發(fā)展與耕地保護的關(guān)系，嚴守耕地紅線，穩(wěn)定耕地面積，同時，促進耕地集約利用，保證糧食作物播種面積。

糧食安全狀況的變化表明，山東省糧食生產(chǎn)能力逐漸提高，但糧食供給穩(wěn)定性較差，依然存在糧食安全風(fēng)險。一方面可能與農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整有關(guān)，該階段，山東省農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)經(jīng)歷著由糧食—經(jīng)濟作物二元結(jié)構(gòu)向糧食—經(jīng)濟—飼料作物三元結(jié)構(gòu)轉(zhuǎn)變的“陣痛期”，糧食作物播種面積大幅減少；另一方面可能與自然災(zāi)害頻發(fā)有關(guān)，欒健等[36]研究指出自然災(zāi)害對山東省糧食總產(chǎn)量具有顯著的負向影響；葉明華等[37]研究也表明山東省農(nóng)業(yè)氣象災(zāi)害會對糧食產(chǎn)量短期波動產(chǎn)生沖擊作用。因此，為增強糧食供給穩(wěn)定性，應(yīng)進一步完善農(nóng)田水利基礎(chǔ)設(shè)施建設(shè)，保證農(nóng)業(yè)用水有效供給；建立多層次的農(nóng)業(yè)氣象災(zāi)害防范和管理體系，提高氣象災(zāi)害的監(jiān)測與預(yù)報能力；積極研發(fā)和培育抗災(zāi)新品種，增強糧食作物自身的耐旱、耐澇特性，提高糧食作物的抗災(zāi)穩(wěn)產(chǎn)能力。

本研究僅從數(shù)據(jù)上分析了1995—2015年山東省糧食生產(chǎn)時空格局演變規(guī)律，但對于這種時空格局變化帶來的更深層次的影響未能展開討論。山東省糧食生產(chǎn)重心的遷移和糧食生產(chǎn)的空間聚集，清晰地預(yù)示著區(qū)域功能分異的特征，由此產(chǎn)生的經(jīng)濟和環(huán)境影響值得進一步研究。