摘 要 耕地是我國(guó)最為寶貴的資源，是保障國(guó)家糧食安全的基礎(chǔ)。以廣西北部灣經(jīng)濟(jì)區(qū)為研究區(qū)，基于基礎(chǔ)地理信息數(shù)據(jù)和2014—2021年統(tǒng)計(jì)年鑒數(shù)據(jù)，利用重心模型、敏感度分析模型等分析其耕地和糧食生產(chǎn)的時(shí)空變化特征，以及糧食生產(chǎn)對(duì)耕地面積變化的敏感程度。結(jié)果表明： 1）2014—2021年北部灣經(jīng)濟(jì)區(qū)耕地面積整體呈下降趨勢(shì)，糧食產(chǎn)量表現(xiàn)為逐年遞減—緩慢增加的過程，糧食生產(chǎn)、耕地的空間格局與總量變化具有一定差異。2）2014—2021年北部灣經(jīng)濟(jì)區(qū)耕地重心和糧食生產(chǎn)重心的移動(dòng)方位和軌跡具有一定差異，耕地重心表現(xiàn)為向西偏南方向移動(dòng)，糧食生產(chǎn)重心表現(xiàn)為向西偏北方向移動(dòng)。3）對(duì)比2014—2018年、2018—2021年兩個(gè)時(shí)期，北部灣經(jīng)濟(jì)區(qū)糧食產(chǎn)量對(duì)耕地變化的敏感度逐漸減弱。

關(guān)鍵詞 糧食生產(chǎn)；耕地；時(shí)空變化；北部灣經(jīng)濟(jì)區(qū)

中圖分類號(hào)：F326.11 文獻(xiàn)標(biāo)志碼：A DOI：10.19415/j.cnki.1673-890x.2024.19.020

耕地是我國(guó)最為寶貴的資源，是保障國(guó)家糧食安全的基礎(chǔ)。近年來，我國(guó)城鎮(zhèn)化發(fā)展迅速，耕地的非農(nóng)化利用不可避免。根據(jù)第三次全國(guó)國(guó)土調(diào)查結(jié)果，2019年末全國(guó)耕地19.18億畝，自第二次全國(guó)土地調(diào)查以來的10年間減少了1.13億畝[1]，每年耕地減少約75.3萬hm2，耕地保護(hù)和糧食安全面臨嚴(yán)峻挑戰(zhàn)。2021年，國(guó)家“十四五”規(guī)劃和中央一號(hào)文件相繼提出耕地保護(hù)和糧食生產(chǎn)的重要性，開展耕地與糧食產(chǎn)量變化研究，把握各地耕地變化和糧食生產(chǎn)規(guī)律，對(duì)優(yōu)化糧食生產(chǎn)布局、因地制宜制定耕地保護(hù)和糧食發(fā)展政策具有重要意義。針對(duì)耕地空間變化及其與糧食生產(chǎn)的相互關(guān)系，國(guó)內(nèi)學(xué)者從耕地的數(shù)量、質(zhì)量、利用強(qiáng)度及其變化驅(qū)動(dòng)力因子等方面對(duì)糧食生產(chǎn)的影響，以及全國(guó)、糧食主產(chǎn)區(qū)、流域、省域、市域等不同研究尺度開展了大量研究[2-11]，隨著城鎮(zhèn)化背景下耕地不斷減少，糧食生產(chǎn)將持續(xù)受到關(guān)注。

廣西北部灣是我國(guó)大西南地區(qū)重要的出?？?，與東南亞國(guó)家聯(lián)盟海陸相連，近年來其建設(shè)用地面積顯著增長(zhǎng)，耕地面積有所減少[12-13]。隨著西部陸海新通道的加快建設(shè)，北部灣經(jīng)濟(jì)區(qū)對(duì)建設(shè)用地的需求將持續(xù)增加，給耕地保護(hù)帶來新的壓力。本文從區(qū)域經(jīng)濟(jì)合作區(qū)的尺度，分析廣西北部灣經(jīng)濟(jì)區(qū)耕地變化與糧食生產(chǎn)的時(shí)空分異，從敏感度角度分析糧食生產(chǎn)對(duì)耕地?cái)?shù)量變化的依賴程度，定量研究北部灣經(jīng)濟(jì)區(qū)耕地變化和糧食生產(chǎn)的特征規(guī)律，為區(qū)域經(jīng)濟(jì)發(fā)展過程中耕地利用和糧食安全問題解決與政策制定提供參考。

1 "數(shù)據(jù)與方法

1.1 "研究區(qū)概況

北部灣經(jīng)濟(jì)區(qū)（以下簡(jiǎn)稱經(jīng)濟(jì)區(qū)）地處我國(guó)沿海西南端，于2006年批準(zhǔn)成立，是我國(guó)第一個(gè)國(guó)際區(qū)域經(jīng)濟(jì)合作區(qū)。2014年，經(jīng)濟(jì)區(qū)由南寧、北海、欽州、防城港4市延及玉林和崇左市，形成由6市38區(qū)縣組成的區(qū)域經(jīng)濟(jì)區(qū)[14]，面積7.34萬km2?！笆濉逼谀?020年），經(jīng)濟(jì)區(qū)地區(qū)生產(chǎn)總值占全廣西的比重達(dá)到48.3%[15]。2021年經(jīng)濟(jì)區(qū)以占廣西43.2%的耕地，養(yǎng)活了廣西45.8%的人口，糧食總產(chǎn)量達(dá)到568.29萬t。

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

本文基于2021年北部灣經(jīng)濟(jì)區(qū)縣級(jí)行政區(qū)劃基礎(chǔ)地理信息數(shù)據(jù)，以組成經(jīng)濟(jì)區(qū)的38個(gè)區(qū)縣為研究單元，研究數(shù)據(jù)主要來自2014—2021年廣西統(tǒng)計(jì)年鑒（http：//tjj.gxzf.gov.cn/tjsj/tjnj/），在GIS軟件中將其作為屬性字段賦予矢量圖形并進(jìn)行空間分析。

1.3 "研究方法

本文采用耕地指數(shù)、耕地變化指數(shù)、重心模型對(duì)經(jīng)濟(jì)區(qū)的耕地和糧食產(chǎn)量時(shí)空動(dòng)態(tài)及時(shí)空耦合進(jìn)行研究，在此基礎(chǔ)上利用敏感度分析模型探討糧食產(chǎn)量對(duì)耕地變化的敏感性特征。

1.3.1 "耕地指數(shù)

耕地指數(shù)表達(dá)耕地的集聚稀疏程度，是一定范圍內(nèi)耕地面積占比，本文以縣級(jí)行政區(qū)為劃分單元。耕地指數(shù)具體計(jì)算公式如下[16]：

[Ri=SiCi×100]% （1）

式中，[Ri]表示第[i]個(gè)區(qū)縣的耕地指數(shù)，[Si]表示第[i]個(gè)區(qū)縣的耕地面積，[Ci]表示第[i]個(gè)區(qū)縣的行政區(qū)劃面積。根據(jù)各區(qū)縣的耕地指數(shù)，可將耕地集聚稀疏程度劃分為稀疏區(qū)、一般過渡區(qū)、聚集區(qū)等類型。

1.3.2 "耕地變化指數(shù)

耕地變化指數(shù)表達(dá)各區(qū)縣不同時(shí)點(diǎn)間范圍內(nèi)的變化水平，用于劃分各區(qū)縣耕地變化強(qiáng)度級(jí)別，可采用分級(jí)圖進(jìn)行表示，其計(jì)算公式如下[16]：

[δi=（Sti?St?1i）Sti×100]% （2）

式中，[δi]為第[i]個(gè)區(qū)縣的耕地變化指數(shù)，[Sti]表示第[i]個(gè)區(qū)縣在t時(shí)間點(diǎn)的耕地面積，[St?1i]表示第[i]個(gè)區(qū)縣在t-1時(shí)間點(diǎn)的耕地面積。根據(jù)各區(qū)縣的耕地變化指數(shù)，可將耕地變化強(qiáng)度劃分為減少區(qū)、基本穩(wěn)定區(qū)、增長(zhǎng)區(qū)等級(jí)別。

1.3.3 "重心模型

重心模型是分析研究區(qū)域發(fā)展中耕地、糧食等要素空間變動(dòng)的重要工具，要素的重心移動(dòng)客觀反映了區(qū)域發(fā)展各要素地理集聚特征及其位移規(guī)律。重心坐標(biāo)計(jì)算公式為[17]：

[Xt=i=1nFtixii=1nFti] （3）

[Yt=i=1nFtiyii=1nFti] "（4）

式中，[Xt]、[Yt]分別表示某要素t年的區(qū)域重心坐標(biāo)，[Fti]表示t年第[i]個(gè)區(qū)縣的某要素總量，[xi]、[yi]分別表示第[i]個(gè)區(qū)縣的地理中心坐標(biāo)。

設(shè)第t、t+m年某要素的區(qū)域重心坐標(biāo)分別為[（Xt，Yt）]、[（Xt+m，Yt+m）]，則重心由t年向t+m年移動(dòng)的距離為：

[Dm=Xt+m?Xt+（Yt+m?Yt）] "（5）

1.3.4 "敏感度分析模型

為分析糧食生產(chǎn)對(duì)耕地面積變化的響應(yīng)程度，引入敏感度評(píng)價(jià)模型[17]：

[β=（Gt+1?Gt）Gt（St+1?St）St] （6）

式中，[β]為某區(qū)縣糧食產(chǎn)量對(duì)耕地面積變化的敏感度系數(shù)，[Gt]、[Gt+1]分別為研究基期和末期的糧食產(chǎn)量，[St]、[St+1]分別為研究基期和末期的耕地面積。若[βlt;0]，表明研究期內(nèi)糧食產(chǎn)量對(duì)耕地面積變化不敏感；若[β]＞0，表明糧食產(chǎn)量變化受耕地面積變化的影響，[β]值越大，表明糧食產(chǎn)量對(duì)耕地面積變化越敏感。根據(jù)各區(qū)縣的敏感度系數(shù)，由高到低可將敏感程度劃分為高度敏感度（[β]＞10）、中度敏感度（5＜[β]≤10）、低度敏感度（0＜[β]≤5）和無敏感度（[β]≤0）4種級(jí)別。

2 "結(jié)果與分析

2.1 "耕地的時(shí)空變化

2014—2021年，經(jīng)濟(jì)區(qū)耕地總量整體上呈下降趨勢(shì)，由2014年的190.0萬hm2降至2021年的148.6萬hm2，總計(jì)減少41.4萬hm2，約占2014年耕地面積的21.8%（見表1）。隨著耕地面積減少和人口增加，2021年經(jīng)濟(jì)區(qū)的人均耕地面積下降至0.065 hm2·人-1，逐漸逼近聯(lián)合國(guó)糧食及農(nóng)業(yè)組織劃定的人均耕地面積0.053 hm2的警戒線，人地矛盾問題突出。

2014—2021年，經(jīng)濟(jì)區(qū)耕地指數(shù)平均值在20%～26%之間，根據(jù)已有耕地空間格局劃分研究[16]，結(jié)合經(jīng)濟(jì)區(qū)耕地實(shí)際情況，將耕地分布劃分為高度稀疏區(qū)（0≤Ri≤10%）、低度稀疏區(qū)（10%＜Ri≤25%）、一般過渡區(qū)（25%＜Ri≤45%）、低度聚集區(qū)（45%＜Ri≤75%）等4種類型（見圖1）。

總體上，經(jīng)濟(jì)區(qū)內(nèi)耕地分布類型主要為低度稀疏區(qū)和一般過渡區(qū)，少部分為高度稀疏區(qū)和低度聚集區(qū)。2014—2021年，經(jīng)濟(jì)區(qū)的低度稀疏區(qū)逐步擴(kuò)展范圍，一般過渡區(qū)逐步縮小范圍；低度稀疏區(qū)從19個(gè)增加到25個(gè)，一般過渡區(qū)從16個(gè)減少到10個(gè)，高度稀疏區(qū)從2個(gè)增加到3個(gè)，低度集聚區(qū)從1個(gè)減少為0，這表明經(jīng)濟(jì)區(qū)各區(qū)縣的耕地聚集程度整體上有所降低。高度稀疏區(qū)集中在經(jīng)濟(jì)區(qū)南部的防城港市，于2017年擴(kuò)展到防城區(qū)等3個(gè)區(qū)縣；一般過渡區(qū)主要集中在經(jīng)濟(jì)區(qū)北部的南寧市、崇左市和經(jīng)濟(jì)區(qū)南部的北海市。2014—2020年，崇左市扶綏縣為經(jīng)濟(jì)區(qū)唯一低度集聚區(qū)，于2021年退化為一般過渡區(qū)。

2.2 "糧食生產(chǎn)的時(shí)空變化

2014—2021年，經(jīng)濟(jì)區(qū)糧食產(chǎn)量表現(xiàn)為逐年遞減—緩慢增加的過程（見表1）。2014—2019年，雖然經(jīng)濟(jì)區(qū)耕地面積減少速度較緩慢，但由于糧食成本不斷上升和經(jīng)濟(jì)作物種植收益較高等原因調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)，導(dǎo)致糧食種植面積逐年下降，糧食總產(chǎn)量隨之下降，2019年降至548.8萬t，為2014年糧食產(chǎn)量的86.0%，是研究期內(nèi)糧食生產(chǎn)的最低谷。2019—2021年，隨著鄉(xiāng)村振興戰(zhàn)略和糧食生產(chǎn)支持等政策的實(shí)施，國(guó)家和地方種糧補(bǔ)貼等政策引導(dǎo)向種植糧食作物調(diào)整，糧食播種面積緩慢增長(zhǎng)，糧食產(chǎn)量也隨之增加，至2021年達(dá)到568.3萬t，但僅恢復(fù)到2014年糧食產(chǎn)量的89.0%，人均糧食產(chǎn)量為2014年的84.7%。

根據(jù)聯(lián)合國(guó)糧農(nóng)組織的糧食標(biāo)準(zhǔn)，人均年糧食需求量是400 kg，結(jié)合有關(guān)研究基礎(chǔ)[18]，按各區(qū)縣人均糧食占有量將糧食生產(chǎn)格局分為嚴(yán)重缺糧區(qū)（＜150 kg）、一般缺糧區(qū)（150～300 kg）、基本自給區(qū)（300～400 kg）、一般余糧區(qū)（400～600 kg）、重要余糧區(qū)（＞600 kg）（見圖2）。

如圖2所示，2014、2018、2021年經(jīng)濟(jì)區(qū)內(nèi)糧食產(chǎn)量達(dá)到基本自給及以上水平的區(qū)縣分別有17、14、16個(gè)，不足總區(qū)縣數(shù)的50%，糧食自給率較低。其中，一般余糧區(qū)和重要余糧區(qū)主要集中在經(jīng)濟(jì)區(qū)北部和玉林市部分區(qū)縣，包括了左右江沖積而成的南寧盆地及玉林盆地，盆地平原地勢(shì)平坦、土壤肥沃、水熱充足，使其成為經(jīng)濟(jì)區(qū)重要的產(chǎn)糧區(qū)。

嚴(yán)重缺糧區(qū)主要分布在經(jīng)濟(jì)區(qū)中部、南部沿海和沿邊境的部分區(qū)縣，2014—2021年嚴(yán)重缺糧區(qū)由7個(gè)擴(kuò)展到13個(gè)區(qū)縣，多為各市主城區(qū)，包括南寧市的青秀區(qū)等5個(gè)區(qū)縣，玉林市的玉州區(qū)，沿海的有北海市海城區(qū)、銀海區(qū)及防城港市港口區(qū)、東興市4個(gè)區(qū)縣，以及邊境崇左市的憑祥市等3個(gè)區(qū)縣。其中，南寧市的青秀區(qū)和玉林市的玉州區(qū)等區(qū)縣由于城鎮(zhèn)化發(fā)展，人口大量流入及耕地用途轉(zhuǎn)變等原因，糧食產(chǎn)量難以滿足當(dāng)?shù)匦枨螅谎睾：脱剡吘车牟糠謪^(qū)縣由于水土資源條件限制，糧食生產(chǎn)能力有限。

2014—2017年，武鳴區(qū)為經(jīng)濟(jì)區(qū)唯一重要余糧區(qū)，但2018年退化為一般余糧區(qū)。武鳴區(qū)于2016年撤縣設(shè)區(qū)納入南寧市城區(qū)建設(shè)，其人口吸納和產(chǎn)業(yè)集聚得到增強(qiáng)，隨著城鎮(zhèn)化建設(shè)等占用耕地和人口增加，以及經(jīng)濟(jì)類農(nóng)作物種植面積增加，使得糧食產(chǎn)量和人均糧食占有量降低。

2.3 "糧食生產(chǎn)重心和耕地重心的時(shí)空耦合

2.3.1 "耕地重心變化

如圖3所示，2014—2021年，經(jīng)濟(jì)區(qū)耕地重心在南寧市邕寧區(qū)內(nèi)移動(dòng)，其中2014—2019年向西北部緩慢移動(dòng)2.3 km，2019—2020年向西南部急速移動(dòng)4.8 km，總體上呈現(xiàn)出向西偏南方向移動(dòng)的趨勢(shì)。2014—2021年耕地重心總位移5.1 km，其中向西移動(dòng)4.2 km，向南移動(dòng)2.9 km，尤以2020年移動(dòng)幅度最大。

2.3.2 "糧食生產(chǎn)重心變化

如圖4所示，2014—2021年，經(jīng)濟(jì)區(qū)糧食生產(chǎn)重心在南寧市橫州市內(nèi)移動(dòng)，其中2014—2017年糧食生產(chǎn)重心主要向西移動(dòng)約2 km，2017—2018年重心突然向北偏西大幅度移動(dòng)3 km，整體上呈現(xiàn)出向西偏北方向移動(dòng)的趨勢(shì)。2014—2021年，糧食生產(chǎn)重心總位移4.6 km，其中向西移動(dòng)3.6 km，向北移動(dòng)2.8 km，尤以2018年移動(dòng)幅度最大。

2.3.3 "糧食生產(chǎn)重心與耕地重心的時(shí)空耦合

耕地的空間分布變化是糧食產(chǎn)量的影響因素，對(duì)糧食生產(chǎn)重心移動(dòng)有一定影響，但兩者的重心位置和移動(dòng)軌跡有很大差異。重心位置方面，2014—2021年，經(jīng)濟(jì)區(qū)耕地重心和糧食生產(chǎn)重心雖然都位于南寧市，但位于不同的區(qū)縣；重心移動(dòng)整體趨勢(shì)方面，耕地重心表現(xiàn)為向西偏南方向移動(dòng)，可能是受經(jīng)濟(jì)區(qū)北部和東部方向一般過渡區(qū)數(shù)量減少等因素影響；而糧食生產(chǎn)重心表現(xiàn)為向西偏北方向移動(dòng)，可能是受經(jīng)濟(jì)區(qū)北部和西部方向的一般余糧區(qū)數(shù)量增多等因素影響。

根據(jù)耕地重心軌跡與糧食生產(chǎn)重心軌跡對(duì)比可知，耕地空間分布變化并不是糧食生產(chǎn)重心的唯一決定性因素，糧食生產(chǎn)重心除了受耕地?cái)?shù)量變化影響外，還與耕地質(zhì)量、糧食種植占比、農(nóng)業(yè)科技投入等因素有關(guān)，因此兩者移動(dòng)方位和軌跡具有一定差異。

2.4 "糧食生產(chǎn)對(duì)耕地變化的敏感度分析

在2014—2021年經(jīng)濟(jì)區(qū)耕地面積和糧食產(chǎn)量時(shí)空變化分析的基礎(chǔ)上，將2014—2018年、2018—2021年兩個(gè)時(shí)期的糧食產(chǎn)量對(duì)耕地變化的敏感度進(jìn)行定量分析評(píng)價(jià)（見圖5）。

2014—2018年，經(jīng)濟(jì)區(qū)共計(jì)32個(gè)區(qū)縣的糧食產(chǎn)量對(duì)耕地變化具有敏感度，6個(gè)區(qū)縣為無敏感度。具有敏感度的區(qū)縣占經(jīng)濟(jì)區(qū)區(qū)縣總數(shù)的84.2%，其中高度敏感性有17個(gè)，中度敏感性有2個(gè)，低度敏感性有13個(gè)。具有敏感性的區(qū)縣在經(jīng)濟(jì)區(qū)各市均有分布，這些地區(qū)在快速城鎮(zhèn)化進(jìn)程中，占用耕地特別是城鎮(zhèn)周圍的水田等優(yōu)質(zhì)土地，而補(bǔ)充耕地的質(zhì)量需要一定耕種周期才能達(dá)到改良效果，因此耕地減少特別是優(yōu)質(zhì)耕地面積減少等原因，導(dǎo)致了糧食產(chǎn)量對(duì)耕地面積變化較為敏感。

2018—2021年，經(jīng)濟(jì)區(qū)共計(jì)15個(gè)區(qū)縣的糧食產(chǎn)量對(duì)耕地變化具有敏感度（全部為低度敏感性），占經(jīng)濟(jì)區(qū)區(qū)縣總數(shù)的39.5%，比前一時(shí)期減少了44.7%，且分布較為零散，可見經(jīng)濟(jì)區(qū)的糧食產(chǎn)量對(duì)耕地變化的敏感性明顯降低。隨著耕地面積的不斷減少和糧食生產(chǎn)能力的減弱，部分區(qū)域糧食產(chǎn)量起伏不明顯。2018—2021年，經(jīng)濟(jì)區(qū)中19個(gè)具有高度敏感性或中度敏感性的區(qū)縣轉(zhuǎn)為低度敏感性或無敏感性，耕地面積雖然較2014—2018年有所減少，但由于糧食作物播種面積占比提升和高標(biāo)準(zhǔn)農(nóng)田建設(shè)提升耕地產(chǎn)能等因素影響，糧食產(chǎn)量有小幅度增加，糧食生產(chǎn)對(duì)耕地變化的敏感性降低。

2.5 "耕地變化級(jí)別分析

2014—2021年，經(jīng)濟(jì)區(qū)耕地變化指數(shù)在[-189%，-7%]之間，結(jié)合已有研究[16，19]，以及經(jīng)濟(jì)區(qū)耕地實(shí)際變化情況，采用Quantile方法將經(jīng)濟(jì)區(qū)耕地變化劃分為高速減少區(qū)（[δt]＜-35%）、中速減少區(qū)（-35%≤[δt]＜-27%）、低速減少區(qū)（-27%≤[δt]＜-13%）和基本穩(wěn)定區(qū)（-13%≤[δt]＜0）4種類型（見圖6）。

高速減少區(qū)包括10個(gè)區(qū)縣，主要分布在南寧市，以及玉林市玉州區(qū)、防城港市東興市，該類型區(qū)縣占經(jīng)濟(jì)區(qū)國(guó)土面積的19.3%，耕地面積減少了17.0萬hm2，對(duì)耕地面積減少的貢獻(xiàn)率為41.1%。

中度減少區(qū)包括9個(gè)區(qū)縣，分布在除防城港市外的其他5市，該類型區(qū)縣占經(jīng)濟(jì)區(qū)國(guó)土面積的28.5%，耕地面積減少了13.6萬hm2，對(duì)耕地面積減少的貢獻(xiàn)率為32.8%。

低度減少區(qū)包括9個(gè)區(qū)縣，主要分布在欽州、北海和玉林3市，以及南寧市青秀區(qū)，該類型區(qū)縣占經(jīng)濟(jì)區(qū)國(guó)土面積的24.3%，耕地面積減少了5.9萬hm2，對(duì)耕地面積減少的貢獻(xiàn)率為14.2%。

基本穩(wěn)定區(qū)包括10個(gè)區(qū)縣，主要分布在崇左市，以及沿海的防城港市上思縣、欽州市浦北縣和北海市主城區(qū)，該類型區(qū)縣占經(jīng)濟(jì)區(qū)國(guó)土面積的27.9%，耕地面積減少了4.9萬hm2，對(duì)耕地面積減少的貢獻(xiàn)率為11.9%。

從總體上看，2014—2021年經(jīng)濟(jì)區(qū)各區(qū)縣耕地均有不同幅度減少，在國(guó)家及地方政策的推動(dòng)下，該區(qū)域內(nèi)耕地轉(zhuǎn)為建設(shè)用地的面積逐年增大[20]。2021年，減少區(qū)28個(gè)區(qū)縣內(nèi)的耕地占經(jīng)濟(jì)區(qū)總耕地面積的68.2%，基本穩(wěn)定區(qū)占31.8%。

3 "結(jié)論與討論

本文利用耕地指數(shù)、耕地變化指數(shù)、重心模型和敏感度分析模型等方法深入分析北部灣經(jīng)濟(jì)區(qū)耕地與糧食生產(chǎn)變化情況，主要結(jié)論如下：

1）北部灣經(jīng)濟(jì)區(qū)糧食生產(chǎn)、耕地的空間格局與總量變化具有一定差異。2014—2021年，經(jīng)濟(jì)區(qū)耕地面積總量整體呈下降趨勢(shì)，各區(qū)縣耕地均有不同幅度減少，研究期內(nèi)總共減少了41.4萬hm2。耕地分布類型主要為低度稀疏區(qū)和一般過渡區(qū)，低度稀疏區(qū)逐步擴(kuò)展范圍，一般過渡區(qū)主要集中在經(jīng)濟(jì)區(qū)北部并逐步縮小范圍，經(jīng)濟(jì)區(qū)的耕地聚集程度整體上有所降低。經(jīng)濟(jì)區(qū)糧食產(chǎn)量表現(xiàn)為逐年遞減—緩慢增加的過程，一般余糧區(qū)和重要余糧區(qū)主要集中在經(jīng)濟(jì)區(qū)北部和玉林市部分區(qū)縣，嚴(yán)重缺糧區(qū)主要分布在經(jīng)濟(jì)區(qū)中部、南部沿海和沿邊境的部分區(qū)縣。

2）2014—2021年，北部灣經(jīng)濟(jì)區(qū)耕地重心和糧食生產(chǎn)重心都位于南寧市內(nèi)，但所處區(qū)縣不同。耕地重心表現(xiàn)為向西偏南方向移動(dòng)，總位移5.1 km；糧食生產(chǎn)重心表現(xiàn)為向西偏北方向移動(dòng)，總位移4.6 km，兩者移動(dòng)方位和軌跡具有一定差異，體現(xiàn)了耕地空間分布變化對(duì)糧食生產(chǎn)重心影響的有限性。

3）對(duì)比2014—2018年、2018—2021年兩個(gè)時(shí)期，北部灣經(jīng)濟(jì)區(qū)糧食產(chǎn)量對(duì)耕地變化的敏感度逐漸減弱。2014—2018年具有敏感度的區(qū)縣占經(jīng)濟(jì)區(qū)區(qū)縣總數(shù)的84.2%，在經(jīng)濟(jì)區(qū)各市均有分布；2018—2021年具有敏感度的區(qū)縣比前一時(shí)期減少了44.7%，分布較為零散，體現(xiàn)出經(jīng)濟(jì)區(qū)糧食產(chǎn)量對(duì)耕地變化的敏感性明顯降低，表明耕地?cái)?shù)量對(duì)糧食生產(chǎn)能力的影響程度逐漸減弱。

北部灣經(jīng)濟(jì)區(qū)是廣西經(jīng)濟(jì)發(fā)展的前沿陣地，是服務(wù)國(guó)家“一帶一路”和西部陸海新通道建設(shè)的重要部分，應(yīng)理性看待該區(qū)域耕地和糧食產(chǎn)量減少的現(xiàn)實(shí)，但較低的糧食自給率會(huì)給區(qū)域耕地保護(hù)和可持續(xù)發(fā)展帶來很大壓力，應(yīng)從糧食安全、生態(tài)安全和經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展長(zhǎng)遠(yuǎn)考慮，根據(jù)主體功能定位合理分工，實(shí)施差別化耕地保護(hù)政策，控制農(nóng)業(yè)種植調(diào)整結(jié)構(gòu)方向和數(shù)量，提高就近糧食保障能力，同時(shí)通過農(nóng)業(yè)科技提高耕地生產(chǎn)力和糧食單產(chǎn)水平，統(tǒng)籌做好糧食安全、經(jīng)濟(jì)發(fā)展和生態(tài)安全工作。耕地?cái)?shù)量變化雖然是影響糧食產(chǎn)量的客觀基礎(chǔ)，但糧食生產(chǎn)還受耕地利用、政策、水熱資源、農(nóng)業(yè)科技水平、經(jīng)濟(jì)收益、勞動(dòng)力等諸多因素影響，本文僅分析了耕地空間分布和數(shù)量變化對(duì)糧食生產(chǎn)的影響，下一步需要考慮其他經(jīng)濟(jì)社會(huì)、自然條件等因素的影響，進(jìn)一步分析耕地與糧食生產(chǎn)空間格局的形成機(jī)理和變化趨勢(shì)等。

參考文獻(xiàn)：

[1] 自然資源部.第三次全國(guó)國(guó)土調(diào)查主要數(shù)據(jù)成果發(fā)布會(huì)[EB/OL]. （2021-08-26）[2023-07-14]. https：//www.mnr.gov.cn/dt/zb/2021/qggtdc/jiabin/.

[2] 傅澤強(qiáng)，蔡運(yùn)龍，楊友孝，等.中國(guó)糧食安全與耕地資源變化的相關(guān)分析[J].自然資源學(xué)報(bào)，2001（4）：313-319.

[3] 余振國(guó)，胡小平.我國(guó)糧食安全與耕地的數(shù)量和質(zhì)量關(guān)系研究[J].地理與地理信息科學(xué)，2003（3）：45-49.

[4] 徐珊，宋戈，李丹，等.東北糧食主產(chǎn)區(qū)耕地資源時(shí)空變化及其對(duì)糧食生產(chǎn)能力的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（21）：1-9.

[5] 余曉敏，李強(qiáng)，張娜.江漢平原糧食產(chǎn)量空間格局分析[J].測(cè)繪地理信息，2020，45（2）：41-44，126.

[6] 黃繪青，王偉，劉莉.長(zhǎng)株潭城市群耕地資源數(shù)量時(shí)空演變及驅(qū)動(dòng)機(jī)制分析[J].測(cè)繪地理信息，2023，48（4）：90-95.

[7] 蔡潔，馬舒悅，李曉平，等.糧食生產(chǎn)能力與耕地利用強(qiáng)度耦合協(xié)調(diào)發(fā)展研究[J].資源開發(fā)與市場(chǎng)，2023，39（4）：400-407，416.

[8] 鮑超，方創(chuàng)琳.長(zhǎng)江流域耕地-糧食-人口復(fù)合系統(tǒng)的動(dòng)態(tài)分析及調(diào)控途徑[J].中國(guó)人口·資源與環(huán)境，2007，96（2）：115-120.

[9] 潘佩佩，楊桂山，王曉萌，等.太湖流域糧食生產(chǎn)對(duì)耕地利用變化動(dòng)態(tài)響應(yīng)分析及預(yù)測(cè)[J].長(zhǎng)江流域資源與環(huán)境，2019，28（10）：2364-2375.

[10] 王國(guó)強(qiáng)，宋艷華.基于耕地質(zhì)量數(shù)量的河南省糧食生產(chǎn)能力研究[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2012，33（1）：49-55，61.

[11] 李娜，鄭新奇，董金瑋.濟(jì)南市糧食生產(chǎn)與耕地變化的分析及預(yù)測(cè)[J].安徽農(nóng)業(yè)科學(xué)，2008（5）：1983-1985.

[12] 陸耀邦，黃嫦慧，梁自力，等.廣西壯族自治區(qū)耕地利用現(xiàn)狀與糧食安全問題研究[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2014，35（5）：66-72.

[13] 潘進(jìn)彬.2000—2020年廣西北部灣沿海土地利用時(shí)空演變分析[J].中國(guó)資源綜合利用，2023，41（3）：83-88，97.

[14] 廣西壯族自治區(qū)人民政府辦公廳.廣西北部灣經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃（2014年修訂）[Z].2014：1-29.

[15] 廣西壯族自治區(qū)人民政府辦公廳.廣西北部灣經(jīng)濟(jì)區(qū)高質(zhì)量發(fā)展“十四五”規(guī)劃[Z].2021：1-28.

[16] 關(guān)興良，方創(chuàng)琳，魯莎莎.中國(guó)耕地變化的空間格局與重心曲線動(dòng)態(tài)分析[J].自然資源學(xué)報(bào)，2010，25（12）：1997-2006.

[17] 劉彥隨，王介勇，郭麗英.中國(guó)糧食生產(chǎn)與耕地變化的時(shí)空動(dòng)態(tài)[J].中國(guó)農(nóng)業(yè)科學(xué)，2009，42（12）：4269-4274.

[18] 譚言飛，濮勵(lì)杰，解雪峰，等.基于敏感度分析的江蘇省糧食生產(chǎn)與耕地?cái)?shù)量變化動(dòng)態(tài)響應(yīng)研究[J].長(zhǎng)江流域資源與環(huán)境，2019，28（5）：1102-1110.

[19] 李鑫，歐名豪.江蘇省耕地空間分布與變化研究[J].農(nóng)業(yè)現(xiàn)代化研究，2011，32（6）：730-734.

[20] 黃天能，許進(jìn)龍.近40年來廣西耕地非農(nóng)化的時(shí)空演變特征、成因與對(duì)策[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃，2023，44（10）：40-51.

（責(zé)任編輯：易 "婧）