中國糧食生產的時空格局及分布動態演進

摘要：從縣域尺度全面考察糧食生產的差異及時空演化特征對提高我國糧食產量和保障糧食安全具有重要參考價值，基于2000—2020年2 545個縣域尺度數據，運用標準差橢圓、三階段嵌套泰爾指數和隨機核密度估計等方法對我國糧食生產的時空格局及分布動態演進進行研究。結果發現，我國糧食生產整體水平以緩慢的速度上升，糧食生產增長率在2000—2020年呈現明顯的波動趨勢，在未來很長一段時期內我國糧食產量將保持穩定但糧食生產增長率會略有下滑。空間分布格局上，我國糧食生產的重心向西北方向遷移，呈現出東北-西南的方向性、東南-西北的發展分散性特征。空間差異上，我國糧食生產的總體差異越來越高，雖然我國糧食生產的縣域差異貢獻率從49.94%下降到29.51%，但縣域尺度依然是總體差異最重要的來源。所以，空間差異主要以縣域差異為主，其次分別是省間差異、市間差異和區域間差異。空間效應上，我國糧食生產水平呈現增長趨勢并存在局部收斂的可能，其空間效應具有連續性，并且在不同時間、縣域空間下具有明顯的正相關空間效應。協同我國區域經濟結構轉型發展趨勢與現實需求，“雙循環”發展格局下細密且深入捕捉區域糧食生產韌性差異，利于各區域因地制宜制定合理有效的糧食生產政策，對于突破糧食生產困境及協調糧食生產安全具有重要意義。

關鍵詞：糧食生產；時空差異；動態演進；隨機Kernel密度估計

doi：10.13304/j.nykjdb.2023.0142

中圖分類號：F326.11 文獻標志碼：A 文章編號：1008‐0864（2024）09‐0001‐11

改革開放以來，我國高度重視糧食生產，相繼頒布了一系列糧食生產支持政策以促進糧食增產、農民增收。在確保國家糧食穩產保供上，取得了舉世矚目的成就，實現了自2004年以來糧食生產的“十七連豐”，牢牢地將糧食安全的主動權掌握在自己手中［1］。1949—1978年，我國糧食生產的政策目標是“糧食增產、保障國內供給、以農促工”，從1978年至今逐漸演變為“保障農民增收、國家糧食安全和糧食生產可持續、高質量發展”，這些變化體現了政策重心從重數量向重結構、重質量發展的趨勢。進入21世紀，我國糧食生產出現了價格“天花板”封頂、生產成本“地板”抬升、資源環境“硬約束”等突出矛盾，糧食生產的長期壓力客觀存在。為了保證糧食供應，我國各省（自治區、直轄市）仍在努力提高糧食產量［2］，國家統計局數據顯示，2018年我國北方人均糧食占有量和南方人均糧食占有量具有較大差異，北方糧食主產區逐步成為南方重要糧倉［3］，致使我國糧食生產的空間結構發生了顯著變化。與此同時，我國農業耕地面積在工業化和城鎮化中逐漸減少，表明現有糧食生產支持政策在應對農業生態環境約束、轉變農業發展方式和促進糧食生產可持續、高質量發展上面臨諸多現實困境，需進一步完善糧食生產的空間格局[4]。因此，新時期保障中國糧食安全、實現糧食供求平衡的目標較以往任何時期都更復雜和艱巨。實踐證明，科學的糧食生產布局是實現中國糧食生產安全目標的重要途徑[5]。在此背景下，以2000—2020年縣域尺度數據作為“十四五”開局之年糧食生產新的開始，從縣域尺度考察我國近20年糧食生產的數量差異及時空演化特征，對于準確把握并合理構建我國糧食生產的新格局，進一步優化配置農業資源、穩步提高我國糧食產量、保障糧食安全具有重要參考價值。

從20世紀90年代開始，我國糧食生產格局的變化就引起廣泛關注。已有研究主要通過糧食產量數據和數學模型對我國糧食生產的時空格局及分布動態演進進行分析[6-11]。隨著我國糧食產量呈現上升的趨勢，糧食生產在地理空間上表現出明顯差異，具體表現為糧食供給水平在東西部地區呈現出顯著差異，糧食生產重心表現出向中部地區集聚的現象[12-14]，呈“ 北上西進”的變動格局［15‐16］，我國糧食生產銷售的空間格局也逐步轉變為“北糧南運”，整體格局發生了巨大變化［17‐18］。從生產潛力上看，我國糧食生產在地理空間上的顯著變化是糧食的總產量呈下降趨勢，這主要是因為耕地保護政策以及城市化發展過程中，糧食生產力強的城市所能利用的耕地越來越少，導致我國糧食生產的空間規模呈現逐漸縮減的變化趨勢[19]。在以糧食生產時空格局為主要內容的研究中發現，低水平協調省區呈現出從西部集聚向全國分散的格局，拮抗型省區主要向沿海地區分散呈帶狀集聚，磨合型省區主要向中西部集聚[20]。綜合分析1978—2018年我國糧食種植面積、總產量和單產水平發現，我國的糧食產量逐年增高，縣域糧食生產的空間效應具有連續性，并且在不同省域空間下具有明顯的正相關空間效應，糧食主產區的空間溢出效應最明顯[21]。在這種現象下，主產區、主銷區、平衡區的糧食生產空間格局呈現“核心-邊緣”式的分布狀態，應在空間效應多重疊加溢出的特征下促進省間糧食生產在波動中穩定提高[22]。通過對我國各類糧食生產能力的時空差異特征研究發現，我國的糧食生產能力逐年增強，但是區域間的差異也在逐年增大[23]。在糧食生產格局顯著變化的情況下，區域間的差異會導致傳統南北糧食生產格局隨時空的演變發生變化，具體表現為由南向北轉移、從東向西擴張的特征[24]。基于對糧食生產的空間差異相關研究發現，糧食生產在不同尺度上均表現出不同的空間特征。以省級單位為研究尺度發現，我國糧食生產的空間分布表現出由南向北、從東到西的遷移特征［25‐26］，導致省級單位之間的糧食生產差距越來越明顯［27］。更深入地從地級行政單元和縣域層次深入揭示我國糧食生產時空分布差異的內在結構特征發現，地級行政單元糧食生產同樣存在空間差異［28］，縣級也如此［29‐30］。總體上看，我國糧食生產在地區、省際、地級行政單元、縣域等不同層次呈現出的空間集聚特征，在縣域尺度上表現得更明顯［31］，說明對縣域尺度糧食生產的研究更能揭示我國糧食生產的時空差異。

從不同尺度、內容上對我國糧食生產的空間格局和動態演進已開展了廣泛研究，但針對全國縣域尺度研究糧食生產的相關報道較少［32］。經濟發展的正確定義是縮小區域發展差距，基于縣級糧食生產的數據，進一步研究區域異質性差異有利于解決社會糧食安全問題[33]。隨著我國經濟的發展轉型，越發重視各地協調發展，同時也應強調糧食生產的韌性，利用糧食生產韌性的杠桿效應突破人多地少、結構老齡化、成本上升等糧食生產困境。因此，在當前“雙循環”新發展格局下[34]，應在時空格局上深入捕獲區域糧食生產差異，制定合理有效的糧食生產政策，推動各地區穩步推進協調糧食安全生產。本研究運用標準差橢圓、三階段嵌套泰爾（Theil）指數和隨機核（Kernel）密度估計等方法，基于2000—2020 年2 545 個縣域尺度的數據，對我國糧食生產的時空格局及分布動態演進進行研究。深入考察我國糧食生產的時空分布差異，并從空間地理層面探究糧食生產在不同階段及區域差異的來源有助于促進區域協調發展，進一步保障我國糧食生產安全。

1 材料與方法

1.1 數據來源

1.1.1 區域劃分

本研究基于2000—2020 年2 545個縣域尺度數據，對我國糧食生產的時空格局及分布動態演進進行研究。根據國家統計局區域劃分辦法及王宏宇等[35]的區域劃分方式，將我國31個省（自治區、直轄市，不包括港澳臺地區）的2 545 個縣域劃分為東部、中部、西部和東北4大區域；其中，東部地區包含625個縣域，中部地區包含639個縣域，西部地區包含1 056個縣域，東北部地區包含225個縣域。

1.1.2 數據來源

本研究所用的糧食產量、年末總人口等數據主要來源于2000—2020年《中國縣域統計年鑒》《中國區域經濟統計年鑒》以及各省（自治區、直轄市）、地級市、縣域的統計年鑒數據。首先，通過搜集與整理相關年鑒形成縣域面板數據；然后，將部分數據的缺失值采用插值法補齊；最后，在剔除無數據的縣域后形成本研究所使用的2 545個縣域的平衡面板數據。

1.1.3 數據選取與預處理

本研究的主體為中國糧食生產的變化情況，因研究視角不同，不同研究者采用了不同的替代變量作為糧食生產衡量指標。有研究者將糧食總產量作為糧食生產的產出變量[36-40]，王躍梅等[40]認為，糧食生產的變化主要得益于糧食產量的提高，糧食產量出現停滯不前現象是自然災害與糧食生產政策共同作用的結果。雖然自然資源稟賦因素（土地、水、糧食播種面積、灌溉面積等）是影響糧食生產格局演變的重要因素[41]，但本研究主要運用糧食總產量作為替代變量研究糧食生產的發展水平以及生產格局演變，同時將數據擴展到縣域層次，避免了自然資源稟賦因素對研究對象糧食生產所產生的數據誤差。所以，鑒于全國縣域數據可得性與完整性，本研究采用糧食總產量作為替代變量來衡量糧食生產的發展水平，能夠比較真實、客觀地反映糧食生產真實情況，以上數據的選取自2000 年開始至2020年為止。

1.2 研究方法

1.2.1 標準差橢圓

標準差橢圓是利用空間地理分布來揭示各種經濟現象空間分布的方法，通常是根據橢圓的重心、方位角、長軸和短軸等參數判斷糧食生產時空發展的重心、方向及演變規律，反映糧食生產在空間上的離散程度，進而確定其時空發展的分布方向以及重心遷移軌跡[42-44]。

式中，SDEe、SDEf分別為糧食重心的經度、緯度坐標；es、fs分別為第s 個縣域的經度、緯度坐標，本研究中經緯度坐標取各縣域多邊形的中心；n為區域樣本量；eˉ、fˉ分別為方位角的經度、緯度坐標。

1.2.2 三階段嵌套Theil指數

2002年首次提出Theil指數可進行2階段分解后，Gustafsson等[45]在同一時間提出了用三階段嵌套Theil指數模型估計不同地理層的空間差異程度。本研究以縣域為基本空間單位拓展糧食生產在空間差異上的來源，從地區、省際、地級行政單元、縣域等不同層次深入揭示我國糧食生產時空分布差異的內在結構特征[46‐47]。

i、j、k、m 分別代表地區、省際、地級市（含省直轄縣）、區縣4 個層級，根據Theil 指數的分解原理，得到第ijk 組內的泰爾指數分解式如下。

式中，Tijk為i 地區j 省k 市的組內差距；Yijkm為i地區j 省k 市m 縣的糧食產量；Pijkm為i 地區j 省k 市m 縣的年末總人口；Yijk為i 地區j 省k 市的糧食產量；Pijk為i 地區j 省k 市的年末總人口。

根據公式（3）可以將Tijk定義為衡量縣域糧食生產的差異程度。然后進一步得到第ij 組內的Theil指數分解式如下。

式中，Yij為i 地區j 省的糧食產量；Pij為i 地區j 省的年末總人口。將Tij分解為地級市的組內差距和組間差距，根據公式（4）將其定義為衡量地級市的糧食生產差異程度，再進一步得到第i 組內的Theil指數分解式如下。

式中，Yi為i 地區的糧食產量；Pi為i 地區的年末總人口。將Ti分解為省內和省間的Theil指數，根據公式（4）將其定義為衡量省際的糧食生產差異程度。最后，得到第i 組內的Theil 指數分解式如下。

式中，Y 為區域糧食產量；P 為地區年末總人口；T 為總的泰爾指數。

將式（3）～（5）逐級嵌套入式（6）中，得到三階段嵌套Theil指數分解如下。

1.2.3 隨機Kernel密度估計

隨機Kernel密度估計是傳統馬爾可夫鏈方法的改進版本，其通過Kernel密度函數得到糧食生產的概率密度估計和三維動態分布圖，描述不同條件下糧食生產的分布形態[48]。此分析基于雙變量核估計得出的結果導入高斯核函數進行計算，其表達式如下。

g （y|x） = f （x，y ） ／f （x） （9）

式中，x 表示獨立分布觀測值位置的經度；y表示獨立分布觀測值位置的緯度；g （y|x）表示位置（y，x）的核密度估計；f（x）表示為x 的邊際核密度函數；f（x，y）表示x 與y 的聯合核密度函數，其表達式如下。

式中，N 是觀察次數；hx表示獨立樣本觀測值經度的帶寬；hy表示獨立樣本觀測值緯度的帶寬；K 表示正態核函數；Xi-x 表示不同年份該縣糧食生產的觀察值；Yi-x 表示不同年份該縣糧食生產的變化值。

2 結果與分析

2.1 我國糧食生產的空間分布特征

由圖1可知，在考察期內，我國糧食生產水平呈現出波動中平穩增長的趨勢，具體表現為全國糧食產量從2003 的4.30 億t 增長至2011 年的6.30 億t，但在2012—2020 年從6.50 億t 下降至6.35 億t，然后又回升到6.50 億t，表明全國糧食增長的趨勢較為穩定。從區域糧食生產水平來看，截至2020年，中部地區的糧食生產水平達到1.90億t，糧食產量約占全國總產量的30%，是四大區域中的最高水平；其次分別是我國西部和東部地區，2000—2020 年分別從1.38 億t 增長到1.74億t、從1.44億t增長到1.56億t，西部地區約占全國糧食產量的26.5%，東部地區約占全國糧食產量的23.8%；四大區域中的東北部地區糧食生產水平最低，2000—2020年從0.54億t增長到1.31億t。從糧食生產增長率來看，四大區域糧食生產增長率的變化趨勢較相似，從圖1可以看出，2000—2020年總體變化趨勢不大，但是從糧食生產增長率來看，呈現出略有下滑的趨勢。特別是近年來，我國多項文件以“穩定糧食產量、耕地面積”為總基調，成為穩定糧食產量的前提。

2.2 我國糧食生產的重心轉移及標準差橢圓分析

基于2000—2020年我國縣域糧食生產數據，運用重心?標準差橢圓模型，以2000、2005、2010、2015 及2020 年的數據生成不同時期區域的標準差橢圓，并得到重心遷移軌跡以及分布趨勢等參數分析我國糧食生產重心的轉移路徑。從重心分布看，2000—2020 年我國糧食生產的各特征時點重心在33°32′56.0″—35°48′56.4″N、114°00′51.6″—115°31′11.9″E范圍內變動，我國糧食生產的所有重心位于我國幾何中心（34°32′27.0″N、108°55′25.0″E）東南方向，表明東部地區縣域的糧食生產平均水平高于西部地區和北部地區。從重心遷移軌跡來看，2000—2005年由河南省駐馬店市遂平縣向西北方向移動至開封市通許縣，2005—2010年由河南省開封市通許縣向西北方向移動至山東省菏澤市東明縣，2010—2015年逐漸向北轉移至河南省長恒市與安陽市滑縣交界處附近，2015—2020年向西北方向移動至河南省濮陽市范縣。總體來看，我國糧食生產分布的重心向西北方向移動的距離和跨度最大，表明2000—2020年時空動態演變中我國西部地區縣域的糧食產量增長速度高于全國平均水平。

從標準差橢圓分布來看，2000—2020年我國糧食生產基本呈現“東北?西南”的空間分布格局。縱觀橢圓的面積移動趨勢，我國糧食生產的標準差橢圓面積呈現逐漸增大趨勢，2000—2015年從285.47 km2增長到351.25 km2，隨后又下降至2020年的327.41 km2，說明我國糧食生產存在空間發散的趨勢。從轉角的變化趨勢看，2000—2020年轉角在58°07′04.9″—60°41′41.1″范圍內呈現下降趨勢，其中2015 年的轉角最大，達到60°41′41.1″，整體呈現波動式下降的態勢，說明我國糧食生產具有向西北方向發散的可能性。從半軸長度來看，2000—2020年長半軸和短半軸的長度分別由12 km 增加至14 km、7 km 增加至8 km，與橢圓面積變化趨勢相似，也呈現增長趨勢，到2015年達到最高，說明我國縣域的糧食生產分布演進趨勢呈現出以東北到西南方向主導，從東南向西北方向發展的趨勢。我國糧食生產重心向西北轉移、向主產區和功能區聚集是各地經濟發展條件、自然稟賦等不盡相同以及糧食生產適應市場調整的結果，是經濟社會發展到一定階段的必然結果，但也使得北方地區在糧食生產中面臨生態壓力進一步加大、產銷距離進一步拉大、倉儲交通和物流成本升高、不利于分散風險等問題。

2.3 我國糧食生產的空間差異及其分解

根據三階段嵌套Theil 指數分別從地區、省際、地級行政單元、縣域等不同層次深入揭示我國糧食生產的時空分布差異和內在結構特征。由圖2可知，我國糧食生產的總體差異呈一定范圍內波動式上升趨勢，2000—2020 年由0.069 上升至0.199。我國糧食生產的縣域差異在0.034～0.058內呈現上升趨勢，是我國糧食生產總體差異的最主要來源且對總體差異的貢獻率最高，在2000—2020年內其貢獻率從49.94%下降到29.51%。省間差異和市間差異變化趨于一致，呈現緩慢上升趨勢，分別由2000 年的0.012、0.019 上升到2020年的0.051、0.507，但是省間差異和市間差異對總體差異貢獻率變化趨勢卻截然相反，省間差異的貢獻率由17.41% 上升到25.37%，但市間差異的貢獻率由27.99% 下降到25.41%，即省間差異和市間差異分布是我國糧食生產總體差異的第2、第3大來源。地區間差異也呈現出增長的變化趨勢，2000—2020 年其Theil 指數由0.003 上升到0.039，其對總體差異的貢獻率也由4.66%上升到19.71%，成為我國糧食生產總體差異的第4 大來源。

2.4 我國糧食生產的分布動態演進

2.4.1 我國糧食生產的無條件Kernel 密度估計

由圖3可知，在無條件Kernel密度估計下，我國糧食生產主要表現為2個趨勢。當糧食生產水平低于200萬t時，Kernel密度等高線集中分布于45°對角線上，說明第r+3年的我國糧食生產水平和第r 年（某一年記作第r 年）相一致，呈現逐漸增長的趨勢。當糧食生產水平高于200萬t時，Kernel密度等高線主要集中分布于45°對角線下方，較少的部分分布于45°對角線上方，說明從第r 年到第r+3年我國糧食生產呈現出高水平縣域向低水平轉移的傾向。進一步分析等高線分布的方向可以看出，我國糧食產量在200萬～300萬t間存在局部收斂的可能性，說明在無條件下我國高水平縣域的糧食生產水平向低水平縣域轉移，而低水平縣域相對穩定。此結論與對我國各地區縣域糧食生產水平對比分析的結果相符，在樣本考察期內，我國糧食生產整體水平以緩慢的速度上升，在未來很長一段時期我國糧食產量將保持穩定但糧食生產增長率會略有下滑。

2.4.2 我國糧食生產空間條件的靜態Kernel密度估計

由圖3可知，Kernel密度等高線呈現出3個主要集群。若將糧食生產水平以100萬和200萬t作為臨界點分別定義為低、中和高水平，可以發現，當地理鄰接縣域在r 年的糧食生產水平小于100萬t時，Kernel密度等高線集中分布于45°對角線上，呈現低?低集聚的現象，在此階段我國縣域的糧食生產存在顯著的空間正相關關系；當地理鄰接縣域在r 年的糧食生產水平在100萬～200萬t時，Kernel密度等高線集中分布于45°對角線下方位置，呈現中?低集聚的現象，在此階段r 年地理鄰接縣域與r 年本縣域的糧食生產整體處于低水平的空間效應；當地理鄰接縣域在r 年的糧食生產水平大于200萬t時，呈現高?中集聚的現象，密度等高線和呈現出的空間效應與糧食生產水平在100萬～200萬t時相似。這與前述結果相符，我國縣域糧食生產的水平普遍處于低水平階段并表現出空間收斂模式，r 年地理鄰接縣域與r 年本縣域的糧食生產均有明顯的正相關空間效應，使得本縣域的糧食生產與地理鄰接縣域會互相影響。

2.4.3 我國糧食生產空間條件的動態Kernel密度估計

由圖3可知，我國縣域糧食生產空間條件下的動態Kernel估計結果與靜態估計Kernel估計結果基本一致，我國縣域糧食生產的空間效應具有連續性，r 年地理鄰接縣域與r 年以及空間滯后條件下r+3年本縣域的糧食生產均有明顯的正相關空間效應。

從糧食生產水平及發展趨勢可知，改革開放以來中國縣域糧食生產動態分布呈現空間收斂特征，我國東部和中部地區縣域糧食生產的收斂模式呈現由南到北的帶狀聚集模式，西部地區縣域糧食生產呈現西北和西南的塊狀分布。同時，我國各省區糧食生產水平及增長速度差異較大，2000—2020年，北部“糧倉基地”縣域糧食生產持續高速增長，成為我國縣域糧食生產高水平聚集區，東部沿海及多數中部省區的縣域相對糧食生產產量持續下滑，與東部相鄰的河南、河北等縣域糧食生產水平持續增長；緊貼中部省區的西藏和青海的糧食生產長期在低水平范圍內波動；其他西北省區糧食生產水平總體持續上升，但西部省區仍然是我國糧食生產的低水平聚集區。這是由于西部省區在水資源利用方面存在總量少、空間分布不均、利用效率低等問題；在耕地資源方面存在耕地數量少、耕地質量差等問題；在人力資源方面存在農村勞動力轉移、勞動力資源浪費等問題；在糧食資源配置方面存在糧食生產資源配置效率偏低、糧食浪費嚴重等問題。

3 討論

本研究通過對我國糧食生產空間格局的研究發現，我國糧食生產整體水平以緩慢的速度上升，但糧食生產增長率從2000—2020年呈現明顯的波動趨勢，在未來很長一段時期我國糧食產量將保持穩定趨勢但糧食生產增長率會略微下滑。盧奕帆等[49]和褚曉彤等[50]研究發現，基于水足跡視角，2000—2020年中國糧食生產的水足跡和灰水足跡具有明顯的空間差異，且不同糧食種類的生產水足跡增長速度也有所不同。這主要是因為技術進步水平對不同區域糧食生產效率的提升作用不一致，東部地區擁有先進技術、管理和政策支持等條件，其糧食生產效率明顯高于中部及西部地區。但在中后期，在激勵政策支持引導下，我國糧食生產的重心向西北方向移動，呈現出以東北到西南方向主導并從東南向西北方向發展的趨勢。所以，為進一步提升各地區糧食生產效率，應通過“科技”基建，消除地區糧食生產不平衡。雖然我國糧食生產整體水平以緩慢的速度上升，但大部分區域省市縣的糧食生產水平低于全國平均水平。首先，需要加大技術推廣力度，因地制宜發展科技創新驅動發展戰略，以科技進步引領帶動糧食高質量生產；其次，需促進以科技創新服務體系為核心的糧食生產戰略，促進先進科創企業滿足不同地區、各類農業經營主體的先進科技需求，打造專業水平高、覆蓋全產業鏈的糧食生產科技支撐體系，消除地區糧食生產發展的不平衡；最后，需形成以現代農業合作社、現代化農場等為主體的經營類單位，促進糧食生產專業化和規模化發展，同時，發揮我國糧食生產縣域之間的正向影響作用，有效促進周邊地區的糧食生產水平穩步增長，以合理縮小地區糧食生產不平衡。

本研究發現，我國糧食生產的空間差異在一定范圍內呈波動式上升趨勢，縣域差異是總體差異的最主要來源且對總體差異的貢獻率最高，然后分別是省間差異、市間差異和區域間差異，這與江立君等[51]、張啟楠等[52]研究結果一致。主要是因為縣域作為較小的單元，糧食生產的集中分布趨勢和區域差異趨勢呈現雙增強的趨勢。與此同時，我國縣域糧食生產的空間效應具有連續性并且會相互影響，各縣域間糧食生產差距較明顯。所以，為進一步以縣域單元為基準推進糧食生產區域協作機制，以提升糧食生產韌性，應根據區域糧食生產空間分布格局，劃分糧食生產水平較集中的縣域群，通過構建糧食生產區域協作機制，進一步提升糧食生產韌性。此外，需以縣域生產單元為主體構建合理的農業區域協作機制，降低縣域糧食生產對我國糧食生產總體差異形成的影響。

本研究通過對我國糧食生產的分布動態研究發現，在無條件Kernel密度估計下，我國糧食總體生產水平呈現增長趨勢并存在局部收斂的可能性；在空間條件Kernel密度估計下，我國縣域糧食生產的空間效應具有連續性，并且在不同時間、縣域空間下具有明顯的正相關空間效應。在未來一段時間和縣域空間格局下我國的糧食生產會相互影響，伴隨著我國糧食生產區間差異的不斷提升，這勢必會成為今后我國保障糧食安全生產的重要羈絆，這與趙丹丹等[53]的研究結論一致，但與李興鋒等[54]的研究結果有所不同，其以縣域水平非均衡發展為切入點，基于糧食生產要素配置視角，深入分析縣域金融集聚對糧食生產供給的影響，實證檢驗了縣域金融集聚對區域糧食生產水平具有顯著的非線性倒“U”形影響，不同要素投入并非均呈明顯的正相關空間效應。本研究僅考慮其他因素不變情況下全國糧食生產的分布動態演進，并未將金融集聚及投入因素考慮在內，導致結果不完全一致。所以，為進一步縮小區域生產差異，應科學劃分耕地管制區，建立彈性調控管理機制。根據糧食生產空間分布格局及動態演進趨勢，結合南北方現存縣域耕地空間分布差異、數量、利用方式的不同，將潛力耕地依次劃分為核心區、管制區、調控區、儲備調控區等多個分級區域，不同管控區域建立差異化管理機制。整體管控中嚴守永久基本農田剛性約束機制，依據“占補平衡”規則科學統籌規劃基本農田利用方式，利用科技手段預留耕地地力恢復空間，每個縣域單元均應因地制宜建立梯次耕地空間調控約束管理機制，以保障剛性及彈性機制靈活運用，從而科學構建國土空間緩沖調控區域，以保障我國糧食生產安全的合理性、永久性及可持續性。

參 考 文 獻

［1］ 趙聰佳，董曉光，王海帆，等.河南省糧食生產時空格局變化及其驅動因素[J].河南農業大學學報，2022，56（2）：312-322.

ZHAO C J， DONG X G， WANG H F， et al .. Spatial-temporalpattern changes and driving factors of grain yield in Henanprovince [J]. J. Henan Agric. Univ.， 2022，56（2）：312-322.

［2］ WANG J， ZHANG Z， LIU Y. Spatial shifts in grain productionincreases in China and implications for food security [J]. LandUse Policy， 2018， 74： 204-213.

［3］ 何可，宋洪遠.資源環境約束下的中國糧食安全： 內涵、挑戰與政策取向[J]. 南京農業大學學報（社會科學版）， 2021，21（3）： 45-57.

HE K， SONG H Y. China’s food security under the constraintsof resources and environment： connotation， challenges andpolicy orientation [J]. J. Nanjing Agric. Univ. （Soc. Sci.）， 2021，21（3）： 45-57.

［4］ ZHANG Q， ZHANG F， WU G， et al .. Spatial spillover effects ofgrain production efficiency in China： measurement and scope[J/OL]. J. Clean. Prod.， 2021， 278： 121062 [2023-02-01]. https：//doi.org/10.1016/j.jclepro.2020.121062.

［5］ HUANG J， RIDOUTT B G， THORP K R， et al .. Water-scarcityfootprints and water productivities indicate unsustainablewheat production in China [J/OL]. Agric. Water. Manage.，2019， 224： 105744 [2023-02-01]. https：//doi.org/10.1016/j.agwat.2019.105744.

［6］ SHANG E P， XU E Q， ZHANG H Q， et al .. Temporal-spatialtrends in potentially toxic trace element pollution in farmlandsoil in the major grain-producing regions of China [J/OL]. Sci.Rep.， 2019， 9（1）： 19463 [2023-02-01]. https：//doi.org/10.1038/s41598-019-55278-5.

［7］ PROSEKOV A Y， IVANOVA S A. Food security： thechallenge of the present [J]. Geoforum， 2018， 91： 73-77.

［8］ 蘆千文.中國農業生產性服務業：70年發展回顧、演變邏輯與未來展望[J].經濟學家，2019（11）：5-13.

LU Q W. China’s agricultural productive service industry：review of 70 years of development， evolutionary logic andfuture prospects [J]. Economist， 2019（11）：5-13.

［9］ 李亞婷，潘少奇，苗長虹.中國縣域人均糧食占有量的時空格局——基于戶籍人口和常住人口的對比分析[J].地理學報，2014，69（12）：1753-1766.

LI Y T， PAN S Q， MIAO C H. The spatial-temporal patterns of per capita share of grain at the county level in China： acomparation between registered population and residentpopulation [J]. Acta Geogr. Sin.， 2014， 69（12）：1753-1766.

［10］ 盧新海，柯楠，匡兵.中國糧食生產能力的區域差異和影響因素[J].中國土地科學，2020，34（8）：53-62.

LU X H， KE N， KUANG B. Regional differences andinfluencing factors of China’s grain production capacity [J].China Land Sci.， 2020， 34（8）：53-62.

［11］ 于智媛，梁書民.我國不同區域在新階段糧食連增中的貢獻因素分析[J].中國農業資源與區劃， 2017，38（8）：145-150.

YU Z Y， LIANG S M. The contribution factor of the cultivated landuse efficiency in the new stage of grain in different regions in China[J]. Chin. J. Agric. Resour. Region. Plan.， 2017， 38（8）：145-150.

［12］ 杜鷹.中國的糧食安全戰略（下）[J].農村工作通訊，2020（22）：17-21.

［13］ 戈大專，龍花樓，屠爽爽，等.黃淮海地區土地利用轉型與糧食產量耦合關系研究[J].農業資源與環境學報，2017，34（4）：319-327.

GE D Z， LONG H L， TU S S， et al .. Coupling relationshipbetween land use transition and grain yield in Huang-Huai-Haiplain， China [J]. J. Agric. Resour. Environ.， 2017， 34 （4）：319-327.

［14］ 陳喜，陸建珍，汪翔，等. 中國甘薯生產布局變遷及動因分析[J].中國農業資源與區劃，2022，43（2）：1-12.

CHEN X， LU J Z， WANG X， et al .. Analysis on layoutevolution of sweet potato production in China and itsinfluencing factors [J]. Chin. J. Agric. Resour. Region. Plan.，2022， 43（2）：1-12.

［15］ 翟榮新，劉彥隨.20世紀90年代以來中國糧食生產的區域格局變動[J].中國農業資源與區劃，2008，29（6）：49-54.

ZHAI R X， LIU Y S. Changes in regional pattern of grainproduction in China since 1990s [J]. Chin. J. Agric. Resour.Region. Plan.， 2008， 29（6）：49-54.

［16］ 何蒲明，魏君英.主銷區與產銷平衡區糧食安全責任分擔問題研究[J].統計與決策， 2020，36（10）：73-76.

［17］ 徐海亞，朱會義.基于自然地理分區的1990—2010年中國糧食生產格局變化[J].地理學報， 2015，70（4）：582-590.

XU H Y， ZHU H Y. Spatial change of China’s grainproduction based on geographical division of natural factorsduring 1990—2010 [J]. Acta Geogr. Sin.， 2015， 70（4）：582-590.

［18］ 鄧宗兵，封永剛，張俊亮，等.中國糧食生產區域格局演變研究[J].農業技術經濟，2013（9）：108-114

［19］ 張秋夢，楊方社，李飛.改革開放以來中國糧食生產空間重構[J].自然資源學報，2021，36（6）：1426-1438.

ZHANG Q M， YANG F S， LI F. The grain production spacereconstruction in China since the reform and opening up [J]. J.Nat. Resour.， 2021， 36（6）：1426-1438.

［20］ 殷偉，姚成勝，黃琳.我國糧食生產與經濟發展耦合協調性的時空演變[J].中國農業資源與區劃，2020，41（11）：110-121.

YIN W， YAO C S， HUANG L，et al .. Spatio-temporal evolutionof the coupling coordination between grain production andeconomic development in China [J]. Chin. J. Agric. Resour.Region. Plan.， 2020， 41（11）：110-121.

［21］ 許紅.中國糧食生產的變化趨勢及空間分異研究[J].中國農業資源與區劃，2020，41（9）：146-154.

XU H. Study on the change trend and spatial differentiation ofgrain production in China [J]. Chin. J. Agric. Resour. Region.Plan.， 2020， 41（9）：146-154.

［22］ 馮穎，侯孟陽，姚順波.中國糧食生產空間關聯網絡的結構特征及其形成機制[J].地理學報，2020，75（11）：2380-2395.

FENG Y， HOU M Y， YAO S B， et al .. Structural characteristicsand formation mechanism of spatial correlation network ofgrain production in China [J]. Acta Geogr. Sin.， 2020， 75（11）：2380-2395.

［23］ 閆琰，宋莉莉，王秀東.我國糧食“十一連增”主要因素貢獻分析及政策思考[J].中國農業科技導報， 2016，18（6）：1-8.

YAN Y， SONG L L， WNAG X D. Analysis and rethinking onmain factor contribution for increasing grain production inChina during the eleventh increasing [J]. J. Agric. Sci.Technol.， 2016， 18（6）：1-8.

［24］ LU C X， LIU A M， XIAO Y， et al .. Changes in China’s grainproduction pattern and the effects of urbanization and dietarystructure [J]. J. Resour. Ecol.， 2020， 11： 358-365.

［25］ 鄭有貴，鄺嬋娟，焦紅坡.南糧北調向北糧南運演變成因的探討——兼南北方兩個區域糧食生產發展優勢和消費比較[J].中國經濟史研究，1999（1）：99-106.

［26］ 崔奇峰，周寧，孫翠清.主產區糧食生產貢獻率及地區比較優勢分析——以水稻、小麥和玉米為例[J].農業經濟與管理，2013（2）：35-42.

CUI Q F， ZHOU N， SUN C Q. Contribution and regionalcomparative advantage of grain production in major grainproducing area： taking rice， wheat and corn as example [J].Agric. Econ. Manage.， 2013（2）：35-42.

［27］ 熊友云，張明軍，劉園園，等.中國糧食產量省區差距變化及其成因分析[J].自然資源學報，2009，24（6）：965-974.

XIONG Y Y， ZHANG M J， LIU Y Y， et al .. Analysis of theshift of provinces gaps of China’s grain production andcontributing factors [J]. J. Nat. Resour.， 2009，24（6）： 965-974.

［28］ 張杰，盧李朋，姜朋輝，等.湖北省糧食生產的多尺度及關聯動態分析[J].經濟地理，2013，33（7）：140-147.

ZHANG J， LU L P， JIANG P H， et al .. Multi-scale analysis andgrey analysis of grain production in Hubei province [J]. Econ.Geogr.， 2013， 33 （7）： 140-147.

［29］ 華樹春，鐘鈺.我國糧食區域供需平衡以及引發的政策啟示[J].經濟問題，2021（3）：100-107.

HUA S C， ZHONG Y. Regional balance of grain supply anddemand in China and the policy implication [J]. Econ.Problems， 2021（3）：100-107.

［30］ 劉玉，郭麗英，劉彥隨.1980—2008年環渤海地區縣域糧食的時空動態變化及分析[J].農業工程學報，2012，28（12）：230-236.

LIU Y， GUO L Y， LIU Y S. Spatial-temporal dynamic changeand analysis of grain production in county scale in Bohai Rimfrom 1980 to 2008 [J]. Trans. Chin. Soc. Agric. Eng.， 2012，28 （12）： 230-236.

［31］ 顏曉飛.云南糧食主產地區演變的特征及原因探析[J].中國農業資源與區劃，2015，36（3）：48-55.

YAN X F. Characteristics and causes of the evolution of maingrain producing areas in Yunnan [J]. Chin. J. Agric. Resour.Region. Plan.， 2015， 36（3）：48-55.

［32］ 謝坤，丁明軍，辛良杰，等.中國縣域糧食產量的時空格局及其與經濟發展的空間關系[J]. 經濟地理， 2021， 41（11）：167-175.

XIE K， DING M J， XIN L J， et al .. Spatial-temporal patterns ofgrain output and its spatial relationship with economicdevelopment at the county level in China [J]. Econ. Geogr.，2021，41 （11）： 167-175.

［33］ TORRES-PRECIADO V H， POLANCO-GAYTáN M，TINOCO-ZERME?O M A， et al .. Crime and regional economicgrowth in Mexico： a spatial perspective [J]. Pap. Region. Sci.，2017， 96（3）： 477-494.

［34］ 蔡雪雄，李倩.中國糧食主銷區的糧食生產安全問題研究[J].亞太經濟，2018（5）：130-136.

［35］ 王宏宇，溫紅梅，杜磊.化肥投入量對糧食產量的空間溢出效應分析——基于1978—2016年省級面板數據[J].農業經濟與管理，2020（4）：52-64.

WANG H Y， WEN H M， DU L. Analysis of spatial spillovereffect of fertilizer input on grain yield—based on provincialpanel data from 1978 to 2016 [J]. Agric. Econ. Manage.， 2020（4）：52-64.

［36］ 欒軍強，王榮成，朱子媛，等.山東省糧食供需變動下農業生產格局研究[J].中國農業資源與區劃，2021，42（8）：201-20.

LUAN J Q， WANG R C， ZHU Z Y， et al .. Study on agriculturalproduction pattern under grain supply and demand change inShandong province [J]. Chin. J. Agric. Resour. Region. Plan.，2021， 42（8）：201-209.

［37］ 王佳月，辛良杰.基于GlobeLand30數據的中國耕地與糧食生產的時空變化分析[J].農業工程學報，2017，33（22）：1-8.

WANG J Y， XIN L J. Spatial-temporal variations of cultivatedland and grain production in China based on GlobeLand30 [J].Trans. Chin. Soc. Agric. Eng.， 2017， 33（22）：1-8.

［38］ 王松南，穆月英.我國糧食生產格局演變下的糧地匹配[J].中國農業大學學報，2022，27（9）：1-11.

WANG S N， MU Y Y. Spatial pattern evolution of grainproduction and its impact on grain-land matching in chain [J].J. China Agric. Univ.， 2022， 27 （9）： 1-11.

［39］ 楊騫，司祥慧，王玨.減排增匯目標下中國糧食生產效率的測度及分布動態演進[J].自然資源學報，2022，37（3）：600-615.

YANG Q， SI X H， WANG J， et al .. Analysis of the shift ofprovinces gaps of China’s grain production and contributingfactors [J]. J. Nat. Resour.， 2009， 24 （6）： 965-974.

［40］ 王躍梅，姚先國，周明海.農村勞動力外流、區域差異與糧食生產[J].管理世界，2013（11）：67-76.

［41］ 劉超，朱滿德，陳其蘭.農業機械化對我國糧食生產的影響：產出效應、結構效應和外溢效應[J].農業現代化研究，2018，39（4）：591-600.

LIU C， ZHU M D， CHEN Q L. The effect of agriculturalmechanization on grain production in China： output effect，structure effect and spillover effect [J]. Res. Agric. Mod.， 2018，39（4）：591-600.

［42］ LEFEVER D W. Measuring geographic concentration bymeans of the standard deviational ellipse [J]. Am. J. Soc.，1926， 32（1）： 88-94.

［43］ YAO F G， XUE L Q， LIANG J Y. Spatial and temporalevolution of urban economic development efficiency in China’sYangtze River economic belt from the perspective of sustainabledevelopment [J/OL]. PLoS One， 2022， 17（9）： e0273559 [2023-02-01]. https：//doi.org/10.1371/journal.pone.0273559.

［44］ 孔祥智，張琛，張效榕.要素稟賦變化與農業資本有機構成提高——對1978年以來中國農 業發展路徑的解釋[J].管理世界，2018，34（10）：147-160.

［45］ GUSTAFSSON B， SHI L. Income inequality within and acrosscounties in rural China 1988 and 1995 [J]. J. Dev. Econ.， 2002，69（1）： 179-204.

［46］ HE Q， WEI F， DENG X， et al .. Spatiotemporal pattern ofcarbon productivity and carbon offset potential in Chinesecounties [J/OL]. Sci. Total Environ.， 2022， 846： 157153 [2023-02-01]. https：//doi.org/10.1016/j.scitotenv.2022.157153.

［47］ WU D， YUAN L， LI R， et al .. Decomposing inequality inresearch funding by university-institute sub-group： a threestagenested Theil index [J]. J. Informetr.， 2018， 12（4）： 1312-1326.

［48］ CHEONG T S， LI V J， SHI X P， et al .. Regional disparity andconvergence of electricity consumption in China： a distributiondynamics approach [J/OL]. China Econ. Rev.， 2019， 58： 101154[2023-02-01]. https：//doi.org/10.1016/j.chieco.2018.02.003.

［49］ 盧奕帆，林錦耀.城市協同發展對糧食生產區域的影響預測——以珠三角城市群為例[J].自然資源學報，2023，38（6）：1532-1549.

LU Y F， LIN J Y. Predicting the impact of urban agglomerationdevelopment on grain production regions： a case study of PearlRiver Delta [J]. J. Nat. Resour.， 2023， 38（6）：1532-1549.

［50］ 褚曉彤，王躍朋，周建.關中地區耕地顯性與隱性非糧化的時空變化及對糧食產量的影響[J].農業工程學報，2023，39（19）：241-249.

CHU X T， WANG Y P， ZHOU J. Effects of spatiotemporalvariations in visible and invisible non-grain productioncultivated land on grain yield in the Guanzhong Region ofShaanxi in China [J]. Trans. Chin. Soc. Agric. Eng.， 2023，39（19）：241-249.

［51］ 江立君，張宗毅，鐘鈺，等.動力要素成本結構與糧食生產格局變遷：機理與實證[J].經濟地理，2023，43（8）：154-164，176.

JIANG L J， ZHANG Z Y， ZHONG Y， et al .. The location layoutof grain production and its driving factor cost structure：mechanism and empirical evidence [J]. Econ. Geogr.， 2023，43（8）：154-164，176.

［52］ 張啟楠，張凡凡，麥強，等.中國糧食生產效率空間溢出網絡及提升路徑[J].地理學報，2022，77（4）：996-1008.

ZHANG Q N， ZHANG F F， MAI Q， et al .. Spatial spillovernetworks and enhancement paths of grain production efficiencyin China [J]. Acta Geogr. Sin.， 2022， 77（4）：996-1008.

［53］ 趙丹丹，周宏.農業生產集聚：如何提高糧食生產效率——基于不同發展路徑的再考察[J].農業技術經濟，2020（8）：13-28.

ZHAO D D， ZHOU H. Agricultural production agglomeration：how to improve food production efficiency？ —re-examinationbased on different development paths [J]. J. Agrotechnic. Econ.，2020（8）：13-28.

［54］ 李興鋒，王力.中國糧食可持續安全水平評價、時空演變及驅動機制研究[J].統計與決策，2023，39（19）：92-96.

基金項目：旱區現代灌區生態農業發展范式及效果評價項目（2023A02002-6）。