粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)評價(jià)、區(qū)域差異及驅(qū)動(dòng)因素研究

田璞玉，張丹婷，周燦芳

（廣東省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)經(jīng)濟(jì)與信息研究所/農(nóng)業(yè)農(nóng)村部華南都市農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室/廣東農(nóng)村研究院，廣東 廣州 510640）

【研究意義】據(jù)統(tǒng)計(jì)，2022 年我國常住人口城鎮(zhèn)化率達(dá)65.22%，開始邁進(jìn)以城市型社會(huì)為主體的時(shí)代。但耕地面積不斷減少、農(nóng)業(yè)生產(chǎn)效率低下、農(nóng)產(chǎn)品供給和需求不平衡加劇等問題仍然制約著我國農(nóng)業(yè)農(nóng)村的發(fā)展［1-3］，如何將現(xiàn)代高效農(nóng)業(yè)與快速城鎮(zhèn)化有機(jī)結(jié)合是實(shí)現(xiàn)我國農(nóng)業(yè)農(nóng)村現(xiàn)代化亟需解決的重要問題。都市農(nóng)業(yè)作為推動(dòng)一二三產(chǎn)業(yè)融合、促進(jìn)城鄉(xiāng)協(xié)調(diào)發(fā)展、保障城市高質(zhì)量“菜籃子”供應(yīng)的重要手段備受關(guān)注。目前全國36 個(gè)大中城市耕地面積接近全國耕地總面積的1/9，蔬菜產(chǎn)量占全國的1/6。新冠肺炎疫情防控期間，都市農(nóng)業(yè)在確保城市地區(qū)“米袋子”“菜籃子”量足價(jià)穩(wěn)的過程中發(fā)揮了重要作用［4］。隨著鄉(xiāng)村振興、農(nóng)業(yè)農(nóng)村現(xiàn)代化等戰(zhàn)略的實(shí)施，現(xiàn)代生物、信息等高新技術(shù)和“互聯(lián)網(wǎng)+”等經(jīng)營模式正逐漸改變傳統(tǒng)農(nóng)業(yè)發(fā)展模式，都市農(nóng)業(yè)也呈現(xiàn)出“現(xiàn)代化”特征。因此，在農(nóng)業(yè)農(nóng)村發(fā)展新階段背景下，對現(xiàn)代都市農(nóng)業(yè)發(fā)展進(jìn)行評價(jià)，分析其空間差異和動(dòng)態(tài)演進(jìn)具有重要意義。

【前人研究進(jìn)展】目前，國內(nèi)外關(guān)于都市農(nóng)業(yè)發(fā)展評價(jià)的研究中，指標(biāo)體系的選取不僅強(qiáng)調(diào)經(jīng)濟(jì)產(chǎn)出方面，也考慮生態(tài)效益、社會(huì)效益、農(nóng)業(yè)可持續(xù)發(fā)展等［5-8］，且各個(gè)維度存在權(quán)衡-協(xié)同關(guān)系［9］?？傮w來看，都市農(nóng)業(yè)發(fā)展水平評價(jià)指標(biāo)體系選取的標(biāo)準(zhǔn)因作者對都市農(nóng)業(yè)理解的不同而各有側(cè)重［10］。在評價(jià)方法上主要有層次分析法、信息熵法、綜合指數(shù)法、冪函數(shù)法等［11-13］。針對粵港澳大灣區(qū)的都市農(nóng)業(yè)，周燦芳［14］定量分析珠三角九市的農(nóng)業(yè)發(fā)展優(yōu)勢與劣勢，認(rèn)為其優(yōu)勢在于農(nóng)業(yè)現(xiàn)代化水平和農(nóng)業(yè)社會(huì)化服務(wù)水平相對較高，短板在于農(nóng)業(yè)勞動(dòng)生產(chǎn)效率偏低和財(cái)政對農(nóng)業(yè)的支持保護(hù)力度不足。佟宇競［15］對比國內(nèi)先進(jìn)城市都市農(nóng)業(yè)發(fā)展經(jīng)驗(yàn)，認(rèn)為廣州都市農(nóng)業(yè)發(fā)展存在綜合城市賦能農(nóng)業(yè)優(yōu)勢尚未充分發(fā)揮、農(nóng)業(yè)產(chǎn)業(yè)規(guī)?；徒M織化程度偏低、農(nóng)業(yè)全產(chǎn)業(yè)鏈開發(fā)能力較弱等短板。田璞玉等［16］總結(jié)了粵港澳大灣區(qū)都市農(nóng)業(yè)發(fā)展模式，包括市場流通型、特色產(chǎn)品型、休閑觀光型、科普教育型、生態(tài)保障型等。都市農(nóng)業(yè)影響因素的研究顯示，城鎮(zhèn)化引致都市農(nóng)業(yè)高質(zhì)量發(fā)展，兩者之間耦合關(guān)系顯著［17］。區(qū)域資源環(huán)境的變化決定都市農(nóng)業(yè)發(fā)展的規(guī)模、結(jié)構(gòu)以及發(fā)展水平［18］，資源環(huán)境對都市農(nóng)業(yè)發(fā)展具有約束作用［19］，兩者之間存在協(xié)同發(fā)展關(guān)系，其中資源環(huán)境包括自然和社會(huì)兩方面。城市經(jīng)濟(jì)發(fā)展拉力、資源環(huán)境壓力、制度政策的支持和城市規(guī)劃的空間保障是都市農(nóng)業(yè)發(fā)展的主要驅(qū)動(dòng)力［20］。此外，政策支持也是都市農(nóng)業(yè)發(fā)展的重要因素之一［21］。

【本研究切入點(diǎn)】粵港澳大灣區(qū)是我國開放程度最高、經(jīng)濟(jì)活力最強(qiáng)的都市圈之一，其城鎮(zhèn)化率最高（超過85%），已達(dá)到世界發(fā)達(dá)國家水平。2021 年，大灣區(qū)中的珠三角地區(qū)農(nóng)林牧漁業(yè)總產(chǎn)值為2 926.68 億元、占廣東全省比重35.2%。統(tǒng)計(jì)結(jié)果顯示，大灣區(qū)有900 多萬人口生活在農(nóng)村，傳統(tǒng)的農(nóng)業(yè)經(jīng)營模式已與粵港澳大灣區(qū)的發(fā)展不相匹配。隨著粵港澳大灣區(qū)建設(shè)不斷深入，都市農(nóng)業(yè)成為灣區(qū)內(nèi)部發(fā)達(dá)城市的應(yīng)然選擇，也成為解決灣區(qū)內(nèi)部城鄉(xiāng)發(fā)展不平衡的必然選擇。這為本文研究現(xiàn)代都市農(nóng)業(yè)發(fā)展提供了最佳觀測樣本。前人對都市農(nóng)業(yè)的評價(jià)研究缺乏對農(nóng)業(yè)“現(xiàn)代化”特征的體現(xiàn)，且對粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展的評價(jià)及其驅(qū)動(dòng)因素研究較少?！緮M解決的關(guān)鍵問題】本研究首先構(gòu)建粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)評價(jià)指標(biāo)體系，利用熵值法進(jìn)行綜合評價(jià)，使用Dagum 基尼系數(shù)分解分析區(qū)域差異和動(dòng)態(tài)演進(jìn)，采用地理探測器方法對粵港澳大灣區(qū)都市農(nóng)業(yè)發(fā)展的驅(qū)動(dòng)影響因素進(jìn)行研究。建立對粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平的評價(jià)體系，分析區(qū)域差異與驅(qū)動(dòng)因素，研究驅(qū)動(dòng)因素的作用程度與方式。

1 研究區(qū)域與數(shù)據(jù)來源

1.1 研究區(qū)域

粵港澳大灣區(qū)，包括香港特別行政區(qū)、澳門特別行政區(qū)和廣東省廣州市、深圳市、珠海市、佛山市、惠州市、東莞市、中山市、江門市、肇慶市（以下簡稱“珠三角九市”），總面積5.6 萬km2，第七次全國人口普查數(shù)據(jù)顯示，2020年大灣區(qū)人口為8 336 萬。其中，香港農(nóng)業(yè)生產(chǎn)穩(wěn)定但占比極低，澳門幾乎沒有農(nóng)業(yè)用地和農(nóng)產(chǎn)品生產(chǎn)，其農(nóng)產(chǎn)品主要來自進(jìn)口。因此，本研究重點(diǎn)考察珠三角九市都市農(nóng)業(yè)的發(fā)展。

1994 年編制的《珠三角經(jīng)濟(jì)區(qū)城市群規(guī)劃》和2003 年編制的《珠三角城鎮(zhèn)群協(xié)調(diào)發(fā)展規(guī)劃》相繼提出探索構(gòu)建以廣佛肇（廣州、佛山、肇慶）、深莞惠（深圳、東莞、惠州）和珠中江（珠海、中山、江門）“三大都市區(qū)”為核心的“多中心”區(qū)域空間結(jié)構(gòu)。2008 年《珠三角改革發(fā)展規(guī)劃綱要》出臺(tái)，并在延續(xù)“都市區(qū)”格局的基礎(chǔ)上首次提出“三大都市圈”的多中心空間結(jié)構(gòu)。據(jù)此，本研究將粵港澳大灣區(qū)分為廣佛肇、深莞惠和珠中江3 個(gè)區(qū)域。

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

省級(jí)休閑農(nóng)業(yè)與鄉(xiāng)村旅游示范點(diǎn)數(shù)據(jù)來源于廣東省農(nóng)業(yè)農(nóng)村廳（http://dara.gd.gov.cn/），其余指標(biāo)數(shù)據(jù)來源于歷年《廣東統(tǒng)計(jì)年鑒》《廣東農(nóng)村統(tǒng)計(jì)年鑒》和各地市統(tǒng)計(jì)年鑒，缺失值使用插值法補(bǔ)齊。深圳農(nóng)村、農(nóng)民和農(nóng)民可支配收入數(shù)據(jù)使用珠三角地區(qū)均值替代。現(xiàn)代都市農(nóng)業(yè)發(fā)展驅(qū)動(dòng)因素及計(jì)算公式見表1。

表1 現(xiàn)代都市農(nóng)業(yè)發(fā)展驅(qū)動(dòng)因素及計(jì)算公式Table 1 Driving factors and calculation of modern urban agriculture development

2 研究方法

2.1 指標(biāo)體系構(gòu)建

農(nóng)業(yè)農(nóng)村部都市農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室（上海交通大學(xué)）提出的《中國都市現(xiàn)代農(nóng)業(yè)發(fā)展評價(jià)指標(biāo)體系（UASJTU）》和《2019 年中國現(xiàn)代都市農(nóng)業(yè)競爭力綜合指數(shù)》［22］中的指標(biāo)體系體現(xiàn)了都市農(nóng)業(yè)“現(xiàn)代化”特征?；诖耍狙芯繕?gòu)建了粵港澳大灣區(qū)評價(jià)指標(biāo)體系（表2）。其中，一級(jí)指標(biāo)包含菜籃子產(chǎn)品保障能力、農(nóng)業(yè)生態(tài)與可持續(xù)發(fā)展水平、三產(chǎn)融合發(fā)展水平、農(nóng)業(yè)先進(jìn)生產(chǎn)要素集聚水平、現(xiàn)代農(nóng)業(yè)經(jīng)營水平。二級(jí)指標(biāo)需說明的有以下方面：

主要菜籃子產(chǎn)品保障水平=7 d 平均產(chǎn)量/7 d最低需求量

7 d 平均產(chǎn)量=主要菜籃子產(chǎn)品全年總產(chǎn)量×7/365

7 d 最低需求量=常住人口數(shù)量×（蔬菜人均最低每周需求量+鮮奶人均最低每周需求量+水產(chǎn)品人均最低每周需求量+肉類人均最低每周需求量+禽蛋人均最低每周需求量）

其中，主要菜籃子產(chǎn)品人均最低需求量參照《中國居民膳食指南》（2016）的最低標(biāo)準(zhǔn)，如蔬菜和鮮奶為300 g/d、水產(chǎn)品為250 g/周、肉類和禽蛋為280 g/周。

依據(jù)國家現(xiàn)代農(nóng)業(yè)示范區(qū)指標(biāo)體系中的農(nóng)業(yè)生產(chǎn)全過程能量消耗計(jì)算方式，農(nóng)業(yè)生產(chǎn)全過程的能量消耗總量=機(jī)電排灌能源消耗系數(shù)×機(jī)電排灌面積+農(nóng)用柴油能源消耗系數(shù)×農(nóng)用柴油量+化肥能源消耗系數(shù)×化肥施用量（折純量），其中機(jī)電排灌、農(nóng)用柴油、化肥能源消耗系數(shù)（用氮肥能源消耗系數(shù)代替）分別為0.3869 t/hm2、1.4571 t/t、4.6308 t/t。

2.2 驅(qū)動(dòng)因素分析

選取8 個(gè)因素作為現(xiàn)代都市農(nóng)業(yè)發(fā)展驅(qū)動(dòng)因素進(jìn)行分析，各因素的統(tǒng)計(jì)描述詳見表3。

表3 現(xiàn)代都市農(nóng)業(yè)發(fā)展驅(qū)動(dòng)因素及統(tǒng)計(jì)描述Table 3 Driving factors and statistical description of modern urban agriculture development

2.3 計(jì)算模型

2.3.1 熵值法 熵值法是一種客觀賦權(quán)法，其根據(jù)各項(xiàng)指標(biāo)觀測值所提供的信息大小確定指標(biāo)權(quán)重。本研究將時(shí)間因素考慮在內(nèi)，詳細(xì)參照文獻(xiàn)［23］計(jì)算方式。

2.3.2 Dagum 基尼系數(shù)計(jì)算及分解方法 采用Dagum 基尼系數(shù)計(jì)算及分解方法，參照文獻(xiàn)［24］計(jì)算區(qū)域相對差異及研究分布動(dòng)態(tài)。總基尼系數(shù)G的計(jì)算公式如下：

式中，k為劃分的區(qū)域數(shù)量，本研究中k=3，即廣佛肇、深莞惠、珠中江3 個(gè)區(qū)域，j、h分別代表k個(gè)區(qū)域不同地市，n為每個(gè)區(qū)域內(nèi)的地市數(shù)量，nj、nh為區(qū)域j、h內(nèi)地市的數(shù)量，yji、yhr分別為j、h區(qū)域內(nèi)地市i、r的現(xiàn)代都市農(nóng)業(yè)發(fā)展水平，y為粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展平均值。

Dagum 基尼系數(shù)能夠分解為區(qū)域內(nèi)差距貢獻(xiàn)（Gw）、區(qū)域間凈值差距貢獻(xiàn)（Gnb）和超變密度貢獻(xiàn)（Gt），三者之和等于總基尼系數(shù)。

式中，Gjj、Gjh表示區(qū)域j的基尼系數(shù)和區(qū)域間基尼系數(shù)，計(jì)算公式如下：

2.3.3 地理探測器 地理探測器是探測某一事物的空間分異性并揭示其驅(qū)動(dòng)因素的一種統(tǒng)計(jì)學(xué)方法［25-27］。地理探測器包括因子探測、交互探測、風(fēng)險(xiǎn)探測和生態(tài)探測4 大部分，其中因子探測主要用于檢驗(yàn)?zāi)撤N影響因子是否是導(dǎo)致特定變量空間分異的原因，而交互探測可以識(shí)別影響因子是獨(dú)立起作用還是交互起作用。本研究主要采用因子探測和交互探測對現(xiàn)代都市農(nóng)業(yè)的驅(qū)動(dòng)因素進(jìn)行分析。采用交互探測，假設(shè)U1和U2是影響W空間分異的兩個(gè)因素，通過空間疊加U1和U2形成新圖層U1∩U2，新圖層屬性由U1和U2共同決定。首先計(jì)算U1、U2和疊加圖層U1∩U2對W空間分異的影響程度，即q(U1)、q(U2)、q(U1∩U2)；其次通過比較q(U1)、q(U2)、q(U1∩U2)揭示不同因素Ui之間對W空間分異的交互作用，并評估這種交互作用是否會(huì)增強(qiáng)或減弱單個(gè)因素Ui對W空間分異的影響程度。采用因子探測，假設(shè)存在變量W、Ui（i=1,…,I），比較W的空間分異與Ui的空間分異是否具有顯著一致性，如具有一致性，表明Ui對W的空間分異有影響，其程度用q值度量，如下式所示：

式中，q表示影響空間分異的因素作用，取值范圍是［0,1］，q值越大表明Ui對W空間分異的影響越大，反之越弱；q=0 表明W的空間分異不受Ui的影響，q=1 表明Ui完全控制W的空間分布。N表示全部區(qū)域樣本數(shù)，h=1,…H表示W(wǎng)和Ui的分區(qū)，Nh表示次級(jí)區(qū)域h樣本數(shù)。需要強(qiáng)調(diào)的是，q值是自變量Ui的空間分異解釋因變量W的空間分異，不對q值進(jìn)行顯著性檢驗(yàn)，q仍有明確物理意義，不服從正態(tài)分布假設(shè)［25］。此外，在進(jìn)行空間分異驅(qū)動(dòng)因素的分析時(shí)，首先需對驅(qū)動(dòng)因素進(jìn)行分區(qū)，本研究采用自然斷點(diǎn)法將各驅(qū)動(dòng)因素分為5 個(gè)等級(jí)區(qū)，具體方法詳見http://www.geodetector.cn/。

3 結(jié)果與分析

3.1 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展整體趨勢

從整體來看，粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平呈現(xiàn)逐年增加趨勢（表4），發(fā)展水平由2013 年的0.294 增加到2020 年的0.419，年平均增長率為5.3%。從各地市來看，廣州市由2013年的0.365 增加到2020 年的0.495，年平均增長率為4.4%，是粵港澳大灣區(qū)都市農(nóng)業(yè)發(fā)展水平最高的城市；深圳現(xiàn)代都市農(nóng)業(yè)發(fā)展水平從2013 年的0.351 增加到2020 年的0.361，增長較為緩慢；東莞市從2013 年的0.312 增加到2020 年的0.477，年平均增長率為6.6%，都市農(nóng)業(yè)發(fā)展水平僅次于廣州市；佛山市從2013 年的0.261 增加到2020年的0.410，平均增長率達(dá)到7.2%，在所有城市中增長速度最快；惠州市從2013 年的0.289 增加到2020 年的0.428，平均增長率為6.0%；江門市從2013 年的0.242 增加到2020 年的0.343，年平均增長率為5.2%；肇慶市從2013 年的0.301 增加到2020 年的0.471，年平均增長率為7.1%；中山市從2013 年的0.255 增加到2020 年的0.381，年平均增長率為6.1%；珠海市從2013 年的0.270增加到2020 年的0.409，年平均增長率為6.4%。從不同區(qū)域情況（圖1）來看，各區(qū)域現(xiàn)代都市農(nóng)業(yè)發(fā)展水平均呈上升趨勢。從2015 年開始，廣佛肇地區(qū)超過深莞惠地區(qū)，成為都市農(nóng)業(yè)發(fā)展水平最高的區(qū)域，而珠中江區(qū)域一直處于平均水平以下，與深莞惠、廣佛肇地區(qū)有較大差距。

表4 2013—2020 年粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平評價(jià)結(jié)果Table 4 Evaluation results of modern urban agriculture development level in Guangdong-Hong Kong-Macao Greater Bay Area from 2013 to 2020

圖1 粵港澳大灣區(qū)與各區(qū)域現(xiàn)代都市農(nóng)業(yè)發(fā)展水平Fig.1 Development level of modern urban agriculture in Guangdong-Hong Kong-Macao Greater Bay Area and various regions

3.2 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)各維度評價(jià)結(jié)果

由圖2 可知，2013—2020 年粵港澳大灣區(qū)菜籃子產(chǎn)品保障能力處于平穩(wěn)發(fā)展?fàn)顟B(tài)，發(fā)展水平維持在0.078～0.080 之間。從各城市現(xiàn)代都市農(nóng)業(yè)各維度得分均值（表5）可以看出，菜籃子保障能力較高的城市分別為肇慶（0.114）、惠州（0.098）和廣州（0.097），較低為東莞（0.053）、珠海（0.055）和深圳（0.056）。農(nóng)業(yè)生態(tài)與可持續(xù)發(fā)展水平呈現(xiàn)波動(dòng)上升趨勢，年平均增長率為3.3%，其中惠州和珠海發(fā)展水平較高、年均水平分別為0.121 和0.105，中山較低、平均水平為0.045。三產(chǎn)融合發(fā)展水平年平均增長率為12.1%，其中，東莞和廣州三產(chǎn)融合水平較高、均值分別為0.113和0.090，江門和中山水平較低、均值分別為0.021和0.031。農(nóng)業(yè)先進(jìn)生產(chǎn)要素集聚水平變化不大，呈現(xiàn)上升-下降—上升的波動(dòng)變化趨勢，年平均增長率僅為0.45%，其中肇慶、江門先進(jìn)生產(chǎn)要素集聚程度最高，均值分別為0.063 和0.050，佛山較低（0.032）?，F(xiàn)代農(nóng)業(yè)經(jīng)營水平增長較快，年均增速達(dá)17.2%，其中，中山和佛山相對水平較高、均值分別為0.117 和0.105，江門和肇慶相對較低、均值分別為0.046 和0.051。

圖2 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)各維度發(fā)展情況Fig.2 Development of modern urban agriculture in various dimensions in Guangdong-Hong Kong-Macao Greater Bay Area

表5 2013—2020 年粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)不同維度發(fā)展均值Table 5 Average development of modern urban agriculture in different dimensions in Guangdong-Hong Kong-Macao Greater Bay Area from 2013 to 2020

3.3 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展的區(qū)域差異與驅(qū)動(dòng)因素

3.3.1 區(qū)域差異及其來源 為分析粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平的區(qū)域差異，使用Dagum基尼系數(shù)法計(jì)算2013—2020 年發(fā)展水平的總體基尼系數(shù)，并進(jìn)一步分解、測算出廣佛肇、深莞惠、珠中江三大都市群的基尼系數(shù)，結(jié)果見表6。

表6 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展區(qū)域差異的Dagum 基尼系數(shù)計(jì)算結(jié)果Table 6 Dagum Gini coefficient calculation results of regional differences in modern urban agriculture development in Guangdong-Hong Kong-Macao Greater Bay Area

總體差異與區(qū)域內(nèi)差異方面，粵港澳大灣區(qū)現(xiàn)代都市發(fā)展水平的基尼系數(shù)均呈波動(dòng)下降趨勢，即“下降—上升—下降”，總體從2013 年的0.077 下降到2020 年的0.068；地區(qū)內(nèi)差異由大到小依次為深莞惠、廣佛肇、珠中江。三大城市群中，廣佛肇地區(qū)呈現(xiàn)“下降—短暫上升—逐漸下降”的發(fā)展趨勢，觀察期內(nèi)下降幅度為4.65%；深莞惠地區(qū)呈現(xiàn)“下降—上升—下降”的趨勢，觀察期內(nèi)下降幅度為18.67%；珠中江地區(qū)呈現(xiàn)“上升—下降—上升”趨勢，觀察期內(nèi)區(qū)域內(nèi)基尼系數(shù)上升62.5%，表明珠海、中山和江門3 個(gè)城市都市農(nóng)業(yè)發(fā)展程度差異有逐漸擴(kuò)大的趨勢。

區(qū)域間差異方面，廣佛肇—珠中江間差異最大，其次是深莞惠—珠中江，廣佛肇—深莞惠間差異最小。廣佛肇—深莞惠區(qū)域間差異先增加后縮小，從2013 年的0.068 下降到2020 年的0.065，下降幅度為4.41%；廣佛肇—珠中江區(qū)域間基尼系數(shù)先表現(xiàn)為下降后上升趨勢，從2013 年的0.108下降到2017 年的0.056，再上漲到2020 年的0.097，觀測期內(nèi)總體下降10.19%；深莞惠—珠中江區(qū)域間差距也呈現(xiàn)波動(dòng)下降趨勢，從2013 年的0.098下降到2020 年的0.074，下降幅度為24.49%。

區(qū)域間現(xiàn)代都市農(nóng)業(yè)發(fā)展水平差距是粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平總體差異的主要來源，且呈先下降后上升趨勢，貢獻(xiàn)率從2013 年的60.789%下降到2017 年的36.361%，再上升到2020 年的63.062%；地區(qū)內(nèi)差距貢獻(xiàn)呈“下降—上升—下降”趨勢。超變密度貢獻(xiàn)率總體波動(dòng)變化從2013 年的17.963%上漲到2017 年的38.358%，再下降到2020 年的13.777%。由分解結(jié)果可以看出，城市群間差異是總差異的主要來源，需要進(jìn)一步縮小三大城市圈間的都市農(nóng)業(yè)發(fā)展差距。

3.3.2 驅(qū)動(dòng)因素（1）單因素驅(qū)動(dòng)力?；诘乩硖綔y器方法分析各區(qū)域都市農(nóng)業(yè)發(fā)展差異的驅(qū)動(dòng)因素，表7 顯示了粵港澳大灣區(qū)樣本中各驅(qū)動(dòng)因素對都市農(nóng)業(yè)發(fā)展影響的q值。耕地面積占行政區(qū)域面積的比重（LR）作用強(qiáng)度為0.294，在全部驅(qū)動(dòng)因素中作用程度最大，表明耕地面積越大，則現(xiàn)代都市農(nóng)業(yè)發(fā)展程度越高，耕地資源作為農(nóng)業(yè)發(fā)展的必備要素，對都市農(nóng)業(yè)發(fā)展具有較大的約束作用；其次為三產(chǎn)占比（IS）、人均地區(qū)生產(chǎn)總值（GP）、城鎮(zhèn)化水平（UL）、人均水資源量（WP）、城市人口密度（PD）、路網(wǎng)密度（ND）、都市農(nóng)業(yè)政策（AP），作用強(qiáng)度為0.046，對區(qū)域都市農(nóng)業(yè)發(fā)展也有影響但相對影響有限，表明其對都市農(nóng)業(yè)發(fā)展均有較大的正向影響。

表7 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)各維度單獨(dú)驅(qū)動(dòng)因素影響程度Table 7 Influence degree of all individual driving factors of modern urban agriculture in various dimension in Guangdong-Hong Kong-Macao Greater Bay Area

從各區(qū)域情況來看，廣佛肇地區(qū)都市農(nóng)業(yè)主要驅(qū)動(dòng)因素是三產(chǎn)占比（0.686）、人均地區(qū)生產(chǎn)總值（0.495）和城市人口密度（0.468），都市農(nóng)業(yè)政策的作用程度相對其他區(qū)域較大、達(dá)到0.165。深莞惠地區(qū)都市農(nóng)業(yè)主要驅(qū)動(dòng)因素為人口密度（0.541）和路網(wǎng)密度（0.505），其他因素作用程度相對較低。珠中江地區(qū)都市農(nóng)業(yè)發(fā)展主要驅(qū)動(dòng)因素為三產(chǎn)占比（0.577）、人均地區(qū)生產(chǎn)總值（0.493）和耕地面積占比（0.453）。此外，人均水資源量、城鎮(zhèn)化水平和路網(wǎng)密度也對珠中江區(qū)域都市農(nóng)業(yè)發(fā)展有重要影響。

（2）各因素交互驅(qū)動(dòng)作用。由圖3 可以看出，都市農(nóng)業(yè)驅(qū)動(dòng)因素中任意兩個(gè)因素之間的交互影響均大于其各自的單獨(dú)作用強(qiáng)度，這一結(jié)論在粵港澳大灣區(qū)和各區(qū)域樣本中均成立。以各因素的交互作用為“雙因子增強(qiáng)型”，任意兩個(gè)驅(qū)動(dòng)因素共同作用時(shí)都會(huì)增強(qiáng)其對都市農(nóng)業(yè)發(fā)展的作用強(qiáng)度。但各雙因子組合的驅(qū)動(dòng)力強(qiáng)化幅度有所不同?；浉郯拇鬄硡^(qū)樣本下，交互作用最強(qiáng)為三產(chǎn)占比（IS）與路網(wǎng)密度（ND）、耕地面積占比（LR）與三產(chǎn)占比（IS）、耕地面積占比（LR）與人均地區(qū)生產(chǎn)總值（GP）、路網(wǎng)密度（ND）與人口密度（PD）因素，交互強(qiáng)度達(dá)0.56～0.60；此外，人均地區(qū)生產(chǎn)總值（GP）與人口密度（PD）（0.54）、均地區(qū)生產(chǎn)總值（GP）與城鎮(zhèn)化水平（UL）（0.54）交互作用也較強(qiáng)，表明這些因素產(chǎn)生的驅(qū)動(dòng)力量遠(yuǎn)大于其各自獨(dú)立作用。綜上所述，各個(gè)因素需要彼此配合才能較好地推動(dòng)都市農(nóng)業(yè)發(fā)展。

圖3 粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)驅(qū)動(dòng)因素交互作用Fig.3 Interaction of driving factors of modern urban agriculture in Guangdong-Hong Kong-Macao Greater Bay Area

從各區(qū)域情況（圖4～圖6）來看，各因素的交互作用也均屬于“雙因子增強(qiáng)型”。其中廣佛肇地區(qū)三產(chǎn)占比（IS）與城鎮(zhèn)化水平（UL）（0.81）、三產(chǎn)占比（IS）與耕地面積占比（LR）（0.79）、三產(chǎn)占比（IS）與都市農(nóng)業(yè)政策（IS）（0.78）交互作用最強(qiáng)；三產(chǎn)占比（IS）與人均地區(qū)生產(chǎn)總值（GP）（0.73）、三產(chǎn)占比（IS）與路網(wǎng)密度（ND）（0.70）也顯示出較高交互作用，可見，三產(chǎn)占比（IS）對廣佛肇都市農(nóng)業(yè)發(fā)展具有關(guān)鍵作用。深莞惠地區(qū)交互作用最強(qiáng)的分別為人均地區(qū)生產(chǎn)總值（GP）與路網(wǎng)密度（ND）（0.68）、人均地區(qū)生產(chǎn)總值（GP）與人口密度（PD）（0.67）、三產(chǎn)占比（IS）與人口密度（PD）（0.64）；人均水資源量（WP）與人口密度（PD）（0.63）、城鎮(zhèn)化水平（UL）與人口密度（PD）（0.62）、三產(chǎn)占比（IS）與路網(wǎng)密度（ND）（0.61）也有較強(qiáng)的交互作用。珠中江地區(qū)交互作用最強(qiáng)的分別為三產(chǎn)占比（IS）與城鎮(zhèn)化水平（UL）（0.81）、耕地面積占比（LR）與三產(chǎn)占比（IS）(0.79)、都市農(nóng)業(yè)政策（AP）與三產(chǎn)占比（IS）（0.78）；三產(chǎn)占比（IS）與人均地區(qū)生產(chǎn)總值（GP）（0.73）、三產(chǎn)占比（IS）與路網(wǎng)密度（ND）（0.70）、三產(chǎn)占比（IS）與人口密度（PD）（0.73）等也有較強(qiáng)的交互作用。綜合來看，都市農(nóng)業(yè)發(fā)展的分異是由多個(gè)因素發(fā)揮合理共同作用的結(jié)果，其交互作用均為雙因子增強(qiáng)型。

圖4 廣佛肇地區(qū)現(xiàn)代都市農(nóng)業(yè)驅(qū)動(dòng)因素交互作用Fig.4 Interaction of driving factors of modern urban agriculture in Guangzhou,Foshan and Zhaoqing area

圖5 深莞惠地區(qū)都市農(nóng)業(yè)驅(qū)動(dòng)因素交互作用Fig.5 Interaction of driving factors of modern urban agriculture in Shenzhen,Dongguan and Huizhou area

圖6 珠中江地區(qū)現(xiàn)代都市農(nóng)業(yè)驅(qū)動(dòng)因素交互作用Fig.6 Interaction of driving factors of modern urban agriculture in Zhuhai,Zhongshan and Jiangmen area

4 結(jié)論與建議

4.1 結(jié)論

本研究構(gòu)建了粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)評價(jià)指標(biāo)體系，采用2013—2020 年粵港澳大灣區(qū)九市數(shù)據(jù)測算了現(xiàn)代都市農(nóng)業(yè)發(fā)展水平，并分析了區(qū)域差異和驅(qū)動(dòng)因素，得到主要結(jié)論如下：

整體來看，粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平呈逐年增加趨勢，各區(qū)域發(fā)展水平有差異。大灣區(qū)整體都市農(nóng)業(yè)發(fā)展水平年平均增長率為5.3%，廣州是大灣區(qū)中發(fā)展水平最高的城市。分區(qū)域來看，各區(qū)域現(xiàn)代都市農(nóng)業(yè)發(fā)展水平均呈上漲趨勢。從2015 年開始，廣佛肇地區(qū)超過深莞惠地區(qū)，成為都市農(nóng)業(yè)發(fā)展水平最高的區(qū)域，珠中江地區(qū)處于平均水平以下。分維度來看，粵港澳大灣區(qū)菜籃子保障能力發(fā)展平穩(wěn)，肇慶菜籃子供應(yīng)能力較高，深圳、珠海、東莞等地較低；農(nóng)業(yè)生態(tài)與可持續(xù)發(fā)展水平呈現(xiàn)波動(dòng)上升趨勢，惠州和珠海發(fā)展水平較高；三產(chǎn)融合水平發(fā)展較快、年平均增長率為12.1%，東莞和廣州三產(chǎn)融合水平較高，江門和中山水平較低；農(nóng)業(yè)先進(jìn)生產(chǎn)要素集聚水平增長幅度小、年平均增長率僅為0.45%，肇慶、江門農(nóng)業(yè)先進(jìn)生產(chǎn)要素集聚程度最高，佛山較低；現(xiàn)代農(nóng)業(yè)經(jīng)營水平增長較快、年平均增長率達(dá)17.2%，中山和佛山相對水平較高，江門和肇慶相對較低。

從相對差異來看，粵港澳大灣區(qū)現(xiàn)代都市發(fā)展水平的基尼系數(shù)呈波動(dòng)下降趨勢。其整體差異從2013 年的0.070 下降到2020 年的0.068。地區(qū)內(nèi)差異由大到小依次為深莞惠、廣佛肇和珠中江，觀察期內(nèi)前兩者基尼系數(shù)下降幅度分別為18.67%、4.65%。后者上升62.5%，表明珠海、中山和江門3 個(gè)城市都市農(nóng)業(yè)發(fā)展程度差異有逐漸擴(kuò)大的趨勢。從基尼系數(shù)分解結(jié)果來看，城市群間差異是總差異的主要來源，需要進(jìn)一步縮小三大城市圈間的都市農(nóng)業(yè)發(fā)展差距。

從驅(qū)動(dòng)因素來看，對都市農(nóng)業(yè)作用最大的因素是耕地面積占城市行政區(qū)域比重，作用強(qiáng)度為0.294，三產(chǎn)占比（0.280）、人均地區(qū)生產(chǎn)總值（0.238）、城鎮(zhèn)化水平（0.234）等也有一定程度的影響。各區(qū)域驅(qū)動(dòng)因素作用程度有區(qū)別，廣佛肇地區(qū)都市農(nóng)業(yè)主要驅(qū)動(dòng)因素是三產(chǎn)占比（0.686）、人均地區(qū)生產(chǎn)總值（0.495）、城市人口密度（0.468）；深莞惠地區(qū)主要驅(qū)動(dòng)因素為人口密度（0.541）、路網(wǎng)密度（0.505）；珠中江地區(qū)主要驅(qū)動(dòng)因素是三產(chǎn)占比（0.577）、人均地區(qū)生產(chǎn)總值（0.493）、耕地面積比例（0.453）。各因素的交互作用均屬于“雙因子增強(qiáng)型”，任意兩個(gè)驅(qū)動(dòng)因素共同作用時(shí)都會(huì)增強(qiáng)其對都市農(nóng)業(yè)發(fā)展的作用強(qiáng)度，但各雙因子組合的驅(qū)動(dòng)力強(qiáng)化幅度有所不同。

4.2 建議

粵港澳大灣區(qū)現(xiàn)代都市農(nóng)業(yè)發(fā)展水平存在一定的空間非均衡特征，區(qū)域間發(fā)展差距有擴(kuò)大趨勢，統(tǒng)籌區(qū)域都市農(nóng)業(yè)發(fā)展并避免極化程度加劇是亟需解決的問題。因此，提出以下建議：

一是加強(qiáng)區(qū)域規(guī)劃，促使粵港澳大灣區(qū)各地重視現(xiàn)代都市農(nóng)業(yè)的發(fā)展。各市按照自然稟賦和社會(huì)經(jīng)濟(jì)條件發(fā)展都市農(nóng)業(yè)，推動(dòng)都市農(nóng)業(yè)功能模式推廣與功能擴(kuò)展；在有條件的村鎮(zhèn)探索建設(shè)現(xiàn)代都市農(nóng)業(yè)試驗(yàn)示范區(qū)，重點(diǎn)圍繞都市農(nóng)業(yè)中的菜籃子保障、生態(tài)與可持續(xù)、產(chǎn)業(yè)融合、先進(jìn)要素集聚和現(xiàn)代農(nóng)業(yè)經(jīng)營各維度進(jìn)行試驗(yàn)，探索和推廣與粵港澳大灣區(qū)發(fā)展相適應(yīng)的現(xiàn)代都市農(nóng)業(yè)發(fā)展模式。

二是加快補(bǔ)齊粵港澳大灣區(qū)各城市現(xiàn)代都市農(nóng)業(yè)發(fā)展短板。廣佛肇區(qū)域和深莞惠區(qū)域要強(qiáng)化菜籃子產(chǎn)品供給的主體責(zé)任，減緩區(qū)域間保供能力不平衡、差距大的情況。培育新型農(nóng)業(yè)經(jīng)營主體，鼓勵(lì)企業(yè)化管理機(jī)制，提高農(nóng)民收入，增加肇慶等地的現(xiàn)代農(nóng)業(yè)經(jīng)營水平。深莞惠區(qū)域還要注重倡導(dǎo)發(fā)展綠色農(nóng)業(yè)，鼓勵(lì)支持化肥減量增效，提高農(nóng)業(yè)可持續(xù)發(fā)展能力。珠中江區(qū)域應(yīng)確保耕地保有率，加強(qiáng)農(nóng)產(chǎn)品質(zhì)量安全監(jiān)管，確保區(qū)域內(nèi)菜籃子產(chǎn)品穩(wěn)定供應(yīng)。加強(qiáng)三產(chǎn)融合水平，積極發(fā)展農(nóng)業(yè)加工業(yè)和服務(wù)業(yè)，充分利用地方自然資源稟賦和文化特色發(fā)展鄉(xiāng)村旅游和休閑農(nóng)業(yè)。加強(qiáng)農(nóng)業(yè)技術(shù)推廣服務(wù)，提高農(nóng)業(yè)金融服務(wù)和機(jī)械化水平，增強(qiáng)農(nóng)業(yè)先進(jìn)生產(chǎn)要素集聚水平。

三是聚焦可控變量，保障現(xiàn)代都市農(nóng)業(yè)發(fā)展要素和政策供給。相比三產(chǎn)占比、城鎮(zhèn)化水平等短時(shí)間內(nèi)無法改變的變量，耕地面積、都市農(nóng)業(yè)政策等要素是可控變量，需要給予保障，促使都市農(nóng)業(yè)高質(zhì)量發(fā)展。統(tǒng)籌土地資源，保障都市農(nóng)業(yè)產(chǎn)業(yè)用地。在嚴(yán)防都市農(nóng)地“非糧化”“非農(nóng)化”風(fēng)險(xiǎn)的前提下，應(yīng)適當(dāng)拓寬設(shè)施農(nóng)業(yè)用地范圍，完善建設(shè)用地配給機(jī)制。推進(jìn)農(nóng)村土地經(jīng)營權(quán)集中連片流轉(zhuǎn)，適當(dāng)延長土地租期，穩(wěn)定經(jīng)營主體投入預(yù)期；重點(diǎn)圍繞都市農(nóng)業(yè)的業(yè)態(tài)方向、產(chǎn)業(yè)融合的形態(tài)模式、制度與標(biāo)準(zhǔn)、產(chǎn)業(yè)用地創(chuàng)新、新一代信息技術(shù)與新型經(jīng)營模式（電商、直播）的應(yīng)用、社會(huì)資本與農(nóng)村農(nóng)民的利益聯(lián)結(jié)機(jī)制等方面給予政策支持。強(qiáng)化資金與人才保障，建立和完善公共財(cái)政對都市農(nóng)業(yè)投入持續(xù)增長的長效機(jī)制，鼓勵(lì)各地出臺(tái)完善都市農(nóng)業(yè)財(cái)政支持政策。培育都市農(nóng)業(yè)人才，以培養(yǎng)新型經(jīng)營主體為重點(diǎn)，制定都市農(nóng)業(yè)人才建設(shè)規(guī)劃，通過實(shí)施農(nóng)民素質(zhì)提升戰(zhàn)略和城市人才服務(wù)農(nóng)村戰(zhàn)略。加強(qiáng)都市農(nóng)業(yè)科技研發(fā)，加大對農(nóng)業(yè)科技創(chuàng)新的扶持，充分發(fā)揮科研機(jī)構(gòu)功能，提升服務(wù)都市農(nóng)業(yè)經(jīng)營主體的范圍和力度。通過成果轉(zhuǎn)化建立一批高技術(shù)、高投入、高產(chǎn)出的農(nóng)業(yè)科技創(chuàng)新示范項(xiàng)目。