韋敬楠，劉 偉

(1廣西科技大學(xué)經(jīng)濟(jì)與管理學(xué)院，廣西柳州 545006；2內(nèi)蒙古建筑職業(yè)技術(shù)學(xué)院經(jīng)濟(jì)管理學(xué)院，內(nèi)蒙古呼和浩特 010070；3北京林業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，北京，100083)

0 引言

甘蔗是我國(guó)最主要糖料來(lái)源，也是主要經(jīng)濟(jì)作物之一[1]。廣西作為我國(guó)最大的甘蔗生產(chǎn)基地，甘蔗年產(chǎn)量約占全國(guó)的68%，對(duì)保證國(guó)內(nèi)食糖供給安全具有不可替代的作用[2]。近年來(lái)，受甘蔗種植成本上升、機(jī)械化程度低及比較效益下降等不利影響，廣西甘蔗種植面積和產(chǎn)量呈逐漸下滑趨勢(shì)[3]。為重振蔗糖產(chǎn)業(yè)，廣西于2014年啟動(dòng)建設(shè)“雙高”甘蔗生產(chǎn)示范基地建設(shè)，旨在扭轉(zhuǎn)甘蔗生產(chǎn)嚴(yán)峻形勢(shì)，推動(dòng)廣西蔗糖產(chǎn)業(yè)高質(zhì)量發(fā)展。

目前，關(guān)于農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)相關(guān)研究成果頗為豐富。通過(guò)梳理相關(guān)文獻(xiàn)發(fā)現(xiàn)，不同學(xué)者對(duì)農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)測(cè)度方式不盡相同，主要包括：一是采用規(guī)模比較優(yōu)勢(shì)指數(shù)法、資源稟賦優(yōu)勢(shì)指數(shù)法、效率比較優(yōu)勢(shì)指數(shù)法及綜合比較優(yōu)勢(shì)指數(shù)法等對(duì)我國(guó)及省域農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)進(jìn)行測(cè)度[4-5]；二是通過(guò)采用遙感技術(shù)手段評(píng)價(jià)農(nóng)作物生產(chǎn)布局的比較優(yōu)勢(shì)度，根據(jù)高分辨率遙感影像獲取農(nóng)作物空間分布圖像[6]；三是利用探索性空間數(shù)據(jù)分析方法，結(jié)合GIS及空間計(jì)量經(jīng)濟(jì)學(xué)研究農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)時(shí)空分布特征[7]；四是基于區(qū)域比較理論，綜合考慮生產(chǎn)自然條件、成本等因素，運(yùn)用概率優(yōu)勢(shì)模型分析農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)[8]。

然而，關(guān)于甘蔗生產(chǎn)比較優(yōu)勢(shì)研究較少，吳多廣等[9]研究中國(guó)甘蔗主產(chǎn)區(qū)生產(chǎn)比較優(yōu)勢(shì)，認(rèn)為廣西具有絕對(duì)的規(guī)模比較優(yōu)勢(shì)和綜合比較優(yōu)勢(shì)。鄒曉蔓等[10]研究了我國(guó)甘蔗生產(chǎn)區(qū)域比較優(yōu)勢(shì)，發(fā)現(xiàn)廣西、云南和廣東是我國(guó)甘蔗生產(chǎn)優(yōu)勢(shì)省區(qū)。基于此，本文將采用規(guī)模比較優(yōu)勢(shì)指數(shù)法、效率比較優(yōu)勢(shì)指數(shù)法和綜合比較優(yōu)勢(shì)指數(shù)法對(duì)廣西甘蔗生產(chǎn)區(qū)域比較優(yōu)勢(shì)時(shí)空變化進(jìn)行研究，并借助空間自相關(guān)分析法，探討其空間相關(guān)性問題，通過(guò)構(gòu)建混合效應(yīng)的空間杜賓模型進(jìn)行了實(shí)證檢驗(yàn)，旨在優(yōu)化廣西甘蔗生產(chǎn)布局，提高廣西蔗糖產(chǎn)業(yè)競(jìng)爭(zhēng)力。

1 數(shù)據(jù)來(lái)源及研究方法

1.1 數(shù)據(jù)來(lái)源

根據(jù)本文研究?jī)?nèi)容需要，結(jié)合數(shù)據(jù)獲取可得性，甘蔗產(chǎn)量、種植面積等相關(guān)數(shù)據(jù)來(lái)自《廣西統(tǒng)計(jì)年鑒》(2012～2021)，由于個(gè)別年份存在數(shù)據(jù)缺失，因此，本文擬采用插值法進(jìn)行補(bǔ)充。此外，本文對(duì)廣西 14個(gè)地級(jí)市甘蔗生產(chǎn)比較優(yōu)勢(shì)進(jìn)行空間計(jì)量分析，市域矢量數(shù)據(jù)來(lái)自國(guó)家測(cè)繪地理信息局平臺(tái)網(wǎng)站。

1.2 研究方法

1.2.1 生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)法

測(cè)度農(nóng)作物生產(chǎn)比較優(yōu)勢(shì)的方法眾多，根據(jù)已有的文獻(xiàn)可知[11-12]：綜合比較優(yōu)勢(shì)指數(shù)(Aggregated Advantage Index，簡(jiǎn)稱AAI)是從生產(chǎn)層面上考察農(nóng)作物比較優(yōu)勢(shì)，能夠很好地反映地區(qū)差異。鑒于綜合比較優(yōu)勢(shì)指數(shù)(AAI)是效率比較優(yōu)勢(shì)指數(shù)(Efficiency Advantage Index，簡(jiǎn)稱EAI)和規(guī)模比較優(yōu)勢(shì)指數(shù)(Scale Advantage Index，簡(jiǎn)稱SAI)的綜合，因此，本文擬采用規(guī)模比較優(yōu)勢(shì)指數(shù)法、效率比較優(yōu)勢(shì)指數(shù)法、綜合比較優(yōu)勢(shì)指數(shù)法對(duì)廣西甘蔗生產(chǎn)區(qū)域比較優(yōu)勢(shì)進(jìn)行測(cè)度[11]。

1.2.1.1 規(guī)模比較優(yōu)勢(shì)指數(shù)

規(guī)模比較優(yōu)勢(shì)指數(shù)(SAI)是指從規(guī)模角度去衡量一個(gè)國(guó)家或地區(qū)生產(chǎn)某種產(chǎn)品的比較優(yōu)勢(shì)，本文是指某個(gè)時(shí)期某地甘蔗種植面積占該地區(qū)農(nóng)作物種植面積的比重與同一時(shí)期廣西甘蔗種植面積占廣西農(nóng)作物種植總面積比重的比值。其計(jì)算公式如下：

式(1)中，SAIij為i地區(qū)j農(nóng)作物生產(chǎn)的規(guī)模比較優(yōu)勢(shì)指數(shù)，這里的j農(nóng)作物指的是甘蔗，MSij為i地區(qū)甘蔗種植面積，單位：hm2；Si為i地區(qū)所有農(nóng)作物種植面積之和，單位：hm2；Sj為廣西甘蔗種植面積之和，單位：hm2；S為廣西農(nóng)作物種植面積之和，單位：hm2。

1.2.1.2 效率比較優(yōu)勢(shì)指數(shù)

效率比較優(yōu)勢(shì)指數(shù)(EAI)是指從土地產(chǎn)出率的角度衡量一個(gè)國(guó)家或地區(qū)內(nèi)某種農(nóng)作物的相對(duì)于同期全國(guó)水平的比較優(yōu)勢(shì)，本文是指某個(gè)時(shí)期某地甘蔗單產(chǎn)與同期廣西甘蔗平均單產(chǎn)相比的值，其計(jì)算公式如下：

式(2)中，EAIij為i地區(qū)j農(nóng)作物生產(chǎn)的效率比較優(yōu)勢(shì)指數(shù)，這里的j農(nóng)作物指甘蔗，APij為i地區(qū)甘蔗的單產(chǎn)，單位：t/hm2；APj為廣西甘蔗平均單產(chǎn)，單位：t/hm2。

1.2.1.3 綜合比較優(yōu)勢(shì)指數(shù)

綜合比較優(yōu)勢(shì)指數(shù)(AAI)是效率比較優(yōu)勢(shì)指數(shù)和規(guī)模比較優(yōu)勢(shì)指數(shù)乘積的平方根，它能夠更加全面地反映一個(gè)國(guó)家或地區(qū)在農(nóng)作物生產(chǎn)上的優(yōu)勢(shì)程度，其計(jì)算公式如下：

式(3)中，AAIij為i地區(qū)j農(nóng)作物生產(chǎn)的綜合比較優(yōu)勢(shì)指數(shù)，這里的j農(nóng)作物指甘蔗。若AAIij＞1，說(shuō)明與廣西平均水平相比，i地區(qū)甘蔗生產(chǎn)具有綜合比較優(yōu)勢(shì)；若AAIij＜1，說(shuō)明與廣西平均水平相比，i地區(qū)甘蔗生產(chǎn)不具有綜合比較優(yōu)勢(shì)。AAIij值越大，說(shuō)明i地區(qū)甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)越突出，反之說(shuō)明i地區(qū)甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)越不明顯。

1.2.2 空間自相關(guān)分析

1.2.2.1 全局空間自相關(guān)

全局空間自相關(guān)由Moran于1950年首次提出，用于識(shí)別地理單元之間是否存在相似性，最常見的統(tǒng)計(jì)量是Moran’sI指數(shù)[13]，該指數(shù)取值區(qū)間為[-1，+1]。若指數(shù)大于0時(shí)，說(shuō)明存在正相關(guān)，并且指數(shù)越接近于1其正相關(guān)性越強(qiáng)；若指數(shù)小于0時(shí)，說(shuō)明存在負(fù)相關(guān)，并且指數(shù)越接近于-1其負(fù)相關(guān)性越強(qiáng)；當(dāng)指數(shù)為0時(shí)，說(shuō)明觀測(cè)值呈獨(dú)立隨機(jī)分布，全局空間不相關(guān)。Moran’sI指數(shù)的計(jì)算公式為：

式(4)中，I是全局Moran’sI指數(shù)，n是廣西14個(gè)地級(jí)市，wij為空間權(quán)重，取值分別為1或0，表示市域之間相鄰和不相鄰，xi和xj分別是單元i和單元j的甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)；x-為廣西 14個(gè)地級(jí)市蔗生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)平均值，xi-x-為觀測(cè)值方差。

1.2.2.2 局部空間自相關(guān)

局部空間自相關(guān)用于識(shí)別屬性相似聚集區(qū)的空間分布位置，使用局部Moran’sI指數(shù)Ii表示。局部空間自相關(guān)指數(shù)能直觀顯示廣西 14個(gè)地級(jí)市域間甘蔗生產(chǎn)比較優(yōu)勢(shì)的空間差異程度與顯著性，根據(jù)變化程度劃分為第一象限(高—高)、第三象限(低—低)、第二象限(低—高)和第四象限(高—低)4個(gè)象限。局部空間自相關(guān)Moran’sI指數(shù)Ii的公式為：

式(5)中，Ii為廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)局部空間自相關(guān)指數(shù)，當(dāng)Ii取值為正數(shù)時(shí)，表示廣西市域甘蔗生產(chǎn)比較優(yōu)勢(shì)呈高集聚或低集聚；當(dāng)Ii取值為負(fù)數(shù)時(shí)，表示呈高低集聚或低高集聚。

1.3 變量選取說(shuō)明

農(nóng)業(yè)生產(chǎn)比較優(yōu)勢(shì)是由農(nóng)地資源稟賦、技術(shù)環(huán)境條件和社會(huì)經(jīng)濟(jì)條件共同決定的[15]。因此，本文分析廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響因素時(shí)，主要探討農(nóng)地設(shè)施、農(nóng)業(yè)機(jī)械化、農(nóng)作物種植結(jié)構(gòu)、城鎮(zhèn)化、農(nóng)村基礎(chǔ)設(shè)施等因素對(duì)其影響效應(yīng)。

1.3.1 農(nóng)地設(shè)施

農(nóng)地設(shè)施條件對(duì)甘蔗生產(chǎn)效益的影響表現(xiàn)在農(nóng)地具備灌溉條件，可以有效應(yīng)對(duì)持續(xù)性干旱對(duì)甘蔗生長(zhǎng)的不利影響。此外，機(jī)耕路通達(dá)性對(duì)甘蔗機(jī)械化生產(chǎn)水平產(chǎn)生十分重要的影響，而農(nóng)業(yè)機(jī)械化水平又影響甘蔗生產(chǎn)效益。因此，農(nóng)地設(shè)施條件對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)具有正向影響效應(yīng)。

1.3.2 農(nóng)業(yè)機(jī)械化

提高甘蔗生產(chǎn)效益離不開采用機(jī)械化生產(chǎn)。當(dāng)前，廣西甘蔗生產(chǎn)綜合機(jī)械化水平已經(jīng)提高到68.2%[16]，這對(duì)提高甘蔗生產(chǎn)效益具有不可替代的作用，通過(guò)減少勞動(dòng)力投入，從而降低生產(chǎn)成本，提高甘蔗生產(chǎn)效率。因此，農(nóng)業(yè)機(jī)械化水平對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)具有顯著的正向影響效應(yīng)。

1.3.3 農(nóng)作物種植結(jié)構(gòu)

在廣西農(nóng)村地區(qū)，糧食種植與經(jīng)濟(jì)作物相互爭(zhēng)奪有限的土地。為保障糧食生產(chǎn)安全就意味著需要更多的土地進(jìn)行糧食耕種，從而限制甘蔗種植面積。廣西作為我國(guó)最大的甘蔗生產(chǎn)基地，勢(shì)必要調(diào)整其他作物的比例保證甘蔗生產(chǎn)穩(wěn)定，因此，農(nóng)業(yè)種植結(jié)構(gòu)對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)有一定影響。

1.3.4 農(nóng)村城鎮(zhèn)化

農(nóng)村城鎮(zhèn)化發(fā)展有利于提高農(nóng)民非農(nóng)收入。在農(nóng)村城鎮(zhèn)化進(jìn)程中，農(nóng)民收入結(jié)構(gòu)也將發(fā)生顯著變化，具體表現(xiàn)在種糧收入占農(nóng)民純收入比重持續(xù)下降，經(jīng)濟(jì)作物與非農(nóng)收入占農(nóng)民純收入比重持續(xù)上升。甘蔗作為廣西最主要的經(jīng)濟(jì)作物，但其經(jīng)濟(jì)效益低于水果、桑蠶等特色經(jīng)濟(jì)作物，因此，農(nóng)村城鎮(zhèn)化對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)具有顯著的負(fù)向影響效應(yīng)。

1.3.5 農(nóng)村基礎(chǔ)設(shè)施

完善的基礎(chǔ)設(shè)施對(duì)降低農(nóng)業(yè)生產(chǎn)成本具有顯著影響。甘蔗作為大宗農(nóng)產(chǎn)品，其生產(chǎn)對(duì)農(nóng)村基礎(chǔ)設(shè)施的依賴程度較高，尤其是農(nóng)村道路交通設(shè)施對(duì)甘蔗生產(chǎn)的影響非常大。因此，道路交通對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)具有顯著的正向影響效應(yīng)。

根據(jù)以上分析，本文結(jié)合研究?jī)?nèi)容，考慮到數(shù)據(jù)獲取的可得性，擬選取農(nóng)地灌溉面積比重、人均農(nóng)業(yè)機(jī)械總動(dòng)力、糧食種植面積占比、城鎮(zhèn)化率、公路網(wǎng)密度等5個(gè)相關(guān)指標(biāo)作為分析廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響因素，選取的影響因素和具體變量指標(biāo)詳見表1。

表1 各變量指標(biāo)

2 測(cè)算結(jié)果與分析

2.1 生產(chǎn)比較優(yōu)勢(shì)指數(shù)時(shí)空變化分析

分別利用式(1)、式(2)、式(3)計(jì)算出廣西14個(gè)地級(jí)市 2011年和 2020年甘蔗規(guī)模比較優(yōu)勢(shì)指數(shù)(SAI)、效率比較優(yōu)勢(shì)指數(shù)(EAI)、綜合比較優(yōu)勢(shì)指數(shù)(AAI)，結(jié)果見表 2。運(yùn)用 GeoDa軟件繪制 2011和2020年廣西地級(jí)市甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)的時(shí)空變化圖（圖1、圖2）。

圖1 2011年廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)

圖2 2020年廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)

表2 廣西地級(jí)市甘蔗生產(chǎn)比較優(yōu)勢(shì)指數(shù)

由表2可知，從生產(chǎn)規(guī)模比較優(yōu)勢(shì)指數(shù)(SAI)變化來(lái)看，2011年生產(chǎn)規(guī)模比較優(yōu)勢(shì)指數(shù)大于1的城市分別為崇左、來(lái)賓、防城港、柳州和百色，說(shuō)明這些地區(qū)存在明顯的生產(chǎn)規(guī)模優(yōu)勢(shì)，其中，崇左甘蔗生產(chǎn)規(guī)模比較優(yōu)勢(shì)指數(shù)最高。此外，南寧、北海和河池這3個(gè)地區(qū)甘蔗生產(chǎn)規(guī)模比較優(yōu)勢(shì)指數(shù)(SAI)接近 1，其余地區(qū)甘蔗生產(chǎn)規(guī)模不具備比較優(yōu)勢(shì)。2020年，全區(qū)甘蔗生產(chǎn)規(guī)模比較優(yōu)勢(shì)指數(shù)(SAI)大于1的城市相比2011年明顯增多，崇左一直保持廣西最大的甘蔗產(chǎn)區(qū)[14]?？梢?，廣西甘蔗生產(chǎn)規(guī)模比較優(yōu)勢(shì)在崇左、來(lái)賓、防城港等少數(shù)地區(qū)進(jìn)一步增強(qiáng)，甘蔗種植面積下滑趨勢(shì)得到明顯扭轉(zhuǎn)。

其次，從效率比較優(yōu)勢(shì)指數(shù)(EAI)變化來(lái)看，2011年生產(chǎn)效率比較優(yōu)勢(shì)指數(shù)大于1的城市有8個(gè)，分別為玉林、崇左、貴港、防城港、北海、桂林、梧州和南寧，說(shuō)明這些地區(qū)存在明顯的生產(chǎn)效率優(yōu)勢(shì)，其中，玉林甘蔗生產(chǎn)效率比較優(yōu)勢(shì)指數(shù)最高；其余 6個(gè)地區(qū)甘蔗生產(chǎn)效率指數(shù)均小于 1，說(shuō)明這些地區(qū)不具備甘蔗生產(chǎn)效率比較優(yōu)勢(shì)。2020年，全區(qū)甘蔗生產(chǎn)效率比較優(yōu)勢(shì)指數(shù)(EAI)大于1的城市較2011年減少了2個(gè)，說(shuō)明廣西甘蔗生產(chǎn)效率比較優(yōu)勢(shì)并不突出。

最后，從綜合比較優(yōu)勢(shì)指數(shù)(AAI)變化來(lái)看，2011年甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)大于1的城市有5個(gè)，分別為崇左、來(lái)賓、防城港、柳州和北海，說(shuō)明這些地區(qū)存在明顯的生產(chǎn)綜合比較優(yōu)勢(shì)，其中，崇左甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)最高，其余9個(gè)地區(qū)不具備甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)。2020年，全區(qū)甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)大于1的城市較2011年增加了2個(gè)，說(shuō)明廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)依然保持，崇左、來(lái)賓、防城港甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)更為明顯。

2.2 廣西甘蔗生產(chǎn)規(guī)模化比較優(yōu)勢(shì)空間格局特征分析

根據(jù) 2011～2020年廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)全局自相關(guān) Moran’sI指數(shù)的計(jì)算結(jié)果（表 3），2011～2020年，Moran’sI指數(shù)都為正值，并且P值均小于 0.05，說(shuō)明廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)存在明顯的空間聚集。由于通過(guò)全局Moran’sI指數(shù)值無(wú)法識(shí)別不同地區(qū)甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)在相鄰區(qū)域空間的相關(guān)性，為了進(jìn)一步考察相鄰地市相關(guān)性，本文采用局部Moran’sI指數(shù)進(jìn)行判斷。

表3 2011～2020年廣西甘蔗生產(chǎn)全局自相關(guān)Moran’s I指數(shù)

繪制2011年和2020年LISA聚類圖，如圖3、圖4所示。廣西不同地區(qū)甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)在相鄰區(qū)域空間相關(guān)性較弱，且并無(wú)明顯變化。具體來(lái)看，2011年，廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)高高聚集的城市只有防城港和南寧，低低聚集的城市為賀州和梧州，其他地區(qū)均呈現(xiàn)不明顯的空間集聚特征。到2020年，低低聚集的城市僅有賀州，局部Moran’sI指數(shù)相比2011年基本沒有明顯變化，說(shuō)明廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)在鄰近地區(qū)不呈現(xiàn)明顯的向四周溢出效應(yīng)。

圖3 2011年局部自相關(guān)Moran’s I指數(shù)

圖4 2020年局部自相關(guān)Moran’s I指數(shù)

3 廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響因素分析

3.1 計(jì)量模型構(gòu)建

根據(jù)表 3全局自相關(guān) Moran’sI指數(shù)的P值在0.05的顯著性水平上通過(guò)檢驗(yàn)，說(shuō)明近10年廣西區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)存在明顯的正向空間自相關(guān)性，因而，考慮采用空間面板數(shù)據(jù)進(jìn)行實(shí)證分析。

3.1.1 空間杜賓模型

傳統(tǒng)面板計(jì)量模型由于忽略了空間效應(yīng)而導(dǎo)致估計(jì)結(jié)果是有偏的，面板空間計(jì)量模型則較好地克服傳統(tǒng)面板計(jì)量模型的缺陷，因而在各個(gè)領(lǐng)域得到廣泛運(yùn)用。因此，本文擬采用空間杜賓模型(SDM)來(lái)實(shí)證研究廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響因素。

3.1.2 模型構(gòu)建

根據(jù)空間杜賓模型(SDM)的一般形式，本文擬構(gòu)建如下計(jì)量模型[17]：

式(6)中，EAIit為被解釋變量，表示i地區(qū)t年的甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)，ρ、βi、δi為帶估計(jì)參數(shù)，wit為空間權(quán)重矩陣，xit為解釋變量，包括農(nóng)地灌溉面積比重(x1)、人均農(nóng)業(yè)機(jī)械總動(dòng)力(x2)、糧食種植面積占比(x3)、城鎮(zhèn)化率(x4)、公路網(wǎng)密度(x5)，εit為隨機(jī)誤差項(xiàng)。本文采用鄰接關(guān)系構(gòu)建空間權(quán)重矩陣，基本形式為：

wit=當(dāng)2個(gè)區(qū)域相鄰時(shí)，取值為1；不相鄰時(shí)取值為0。

在空間杜賓模型(SDM)中，普通最小二乘估計(jì)(OLS)方法估計(jì)得到的結(jié)果是有偏且不一致，通常采用極大似然法(QMLE)來(lái)估計(jì)[18]。

3.2 影響因素計(jì)算結(jié)果及分析

一般根據(jù)豪斯曼統(tǒng)計(jì)量(Hausman)檢驗(yàn)結(jié)果選擇固定效應(yīng)還是隨機(jī)效應(yīng)，運(yùn)行Stata 15.0軟件，結(jié)果顯示為負(fù)值，說(shuō)明應(yīng)選擇隨機(jī)效應(yīng)的空間杜賓模型進(jìn)行估計(jì)更為合理。利用Stata 15.0軟件計(jì)算空間杜賓模型系數(shù)，結(jié)果如表4所示。根據(jù)Spatial rho、lgt_theta、sigma2_e的估計(jì)結(jié)果，空間回歸系數(shù)在0.001的水平上顯著，表明存在空間溢出效應(yīng)。

表4 空間杜賓模型(SDM)估計(jì)結(jié)果

由于空間杜賓模型的估計(jì)系數(shù)不能直接反映與被解釋變量的關(guān)系，所以對(duì)空間杜賓模型的估計(jì)結(jié)果進(jìn)一步分解為直接效應(yīng)、間接效應(yīng)和總效應(yīng)?？臻g杜賓模型的直接效應(yīng)表示解釋變量的變化對(duì)本區(qū)域被解釋變量影響的均值，間接效應(yīng)表示解釋變量的變化對(duì)鄰近區(qū)域被解釋變量影響的均值，總效應(yīng)表示解釋變量的變化對(duì)所有區(qū)域被解釋變量影響的均值[19]。根據(jù)各解釋變量影響效應(yīng)估計(jì)結(jié)果見表5。

表5 影響因素效應(yīng)估計(jì)結(jié)果

從 SDM 的直接效應(yīng)來(lái)看，公路網(wǎng)密度對(duì)區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)有顯著的負(fù)向影響，本文認(rèn)為主要原因是廣西甘蔗生產(chǎn)逐漸向交通條件相對(duì)較差的區(qū)域集中，其生產(chǎn)比較優(yōu)勢(shì)進(jìn)一步凸顯。此外，區(qū)域城鎮(zhèn)化率對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)有一定的正向影響，隨著城鎮(zhèn)化水平的提高，非農(nóng)就業(yè)機(jī)會(huì)不斷增多，為甘蔗生產(chǎn)提供更多的社會(huì)化服務(wù)，促進(jìn)了甘蔗生產(chǎn)效率的提高，比較優(yōu)勢(shì)更加明顯。

從 SDM 的間接效應(yīng)來(lái)看，農(nóng)業(yè)機(jī)械化水平對(duì)區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)有顯著正向影響，這說(shuō)明某區(qū)域農(nóng)業(yè)機(jī)械化水平的提高會(huì)對(duì)鄰近區(qū)域農(nóng)業(yè)機(jī)械化水平提高存在一定的空間溢出效應(yīng)。某區(qū)域農(nóng)業(yè)機(jī)械化水平提高后，甘蔗生產(chǎn)比較優(yōu)勢(shì)逐漸凸顯，受其影響，鄰近區(qū)域甘蔗生產(chǎn)機(jī)械化水平也會(huì)逐漸提高，從而促進(jìn)甘蔗生產(chǎn)比較優(yōu)勢(shì)。糧食種植面積占比對(duì)區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)有一定的負(fù)向影響，這說(shuō)明區(qū)域農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整會(huì)影響當(dāng)?shù)馗收嵘a(chǎn)比較優(yōu)勢(shì)，主要表現(xiàn)在增加糧食種植面積后，甘蔗種植面積受到擠壓，因而降低甘蔗生產(chǎn)比較優(yōu)勢(shì)，農(nóng)業(yè)種植結(jié)構(gòu)的調(diào)整會(huì)波及鄰近區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。公路網(wǎng)密度對(duì)鄰近區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)有顯著的負(fù)向影響，甘蔗生產(chǎn)向交通條件相對(duì)較差的區(qū)域集中，相鄰區(qū)域受其影響，也會(huì)因?yàn)榻煌l件相對(duì)較差而選擇種植甘蔗，因而鄰近區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)逐漸凸顯。

從 SDM 的總效應(yīng)來(lái)看，公路網(wǎng)密度和糧食種植面積占比通過(guò)了顯著性檢驗(yàn)，并且影響方向?yàn)樨?fù)，說(shuō)明廣西區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)逐漸向交通條件相對(duì)較差的區(qū)域轉(zhuǎn)移，農(nóng)業(yè)種植結(jié)構(gòu)也在一定程度上影響區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。隨著廣西甘蔗生產(chǎn)進(jìn)一步向崇左、來(lái)賓等地區(qū)集中，其區(qū)域生產(chǎn)比較優(yōu)勢(shì)逐漸凸顯，在一定程度上也表明廣西在不斷優(yōu)化農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)，根據(jù)區(qū)域生產(chǎn)比較優(yōu)勢(shì)合理布局甘蔗種植，從而提高區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。

4 結(jié)論與建議

4.1 研究結(jié)論

本文基于生產(chǎn)綜合比較優(yōu)勢(shì)指數(shù)法、空間自相關(guān)分析法，研究2011～2020年廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)，在此基礎(chǔ)上，運(yùn)用全局Moran’sI指數(shù)檢驗(yàn)廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的空間相關(guān)性，進(jìn)一步構(gòu)建面板空間杜賓模型，實(shí)證研究廣西甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響因素，得出主要結(jié)論有：

(1)近 10年來(lái)，廣西甘蔗生產(chǎn)具備綜合比較優(yōu)勢(shì)的區(qū)域保持基本穩(wěn)定，甘蔗主產(chǎn)區(qū)生產(chǎn)綜合比較優(yōu)勢(shì)主導(dǎo)作用進(jìn)一步增強(qiáng)，區(qū)域生產(chǎn)優(yōu)勢(shì)更加凸顯，呈現(xiàn)向桂西南、桂中地區(qū)進(jìn)一步集中趨勢(shì)。從效率比較優(yōu)勢(shì)和規(guī)模比較優(yōu)勢(shì)變化來(lái)看，桂東南地區(qū)甘蔗生產(chǎn)效率比較優(yōu)勢(shì)最明顯，桂西南和桂中地區(qū)規(guī)模比較優(yōu)勢(shì)進(jìn)一步凸顯。

(2)區(qū)域交通條件對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)具有顯著的負(fù)向影響。實(shí)證結(jié)果表明，本地區(qū)甘蔗生產(chǎn)會(huì)逐漸向交通條件較差的地區(qū)轉(zhuǎn)移，受其影響，相鄰地區(qū)也會(huì)呈現(xiàn)甘蔗生產(chǎn)向交通條件相對(duì)較差的區(qū)域轉(zhuǎn)移。

(3)農(nóng)業(yè)機(jī)械化水平對(duì)甘蔗生產(chǎn)比較優(yōu)勢(shì)有一定的正向影響。機(jī)械化對(duì)甘蔗生產(chǎn)而言非常重要，甘蔗生產(chǎn)機(jī)械化水平較高的區(qū)域不僅對(duì)當(dāng)?shù)馗收嵘a(chǎn)比較優(yōu)勢(shì)具有顯著的促進(jìn)作用，通過(guò)空間溢出效應(yīng)，還對(duì)鄰近區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)有一定的正向影響。

(4)農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整對(duì)區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)在 0.1的水平上顯著，說(shuō)明調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)不僅會(huì)影響到當(dāng)?shù)馗收嵘a(chǎn)比較優(yōu)勢(shì)，鄰近區(qū)域也會(huì)通過(guò)調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)，從而影響區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。城鎮(zhèn)化對(duì)本區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)的影響效應(yīng)較弱但依然顯著，說(shuō)明區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)在一定程度上受市場(chǎng)化程度的影響。

4.2 建議

根據(jù)研究結(jié)論，廣西甘蔗生產(chǎn)綜合比較優(yōu)勢(shì)的區(qū)域差異較為明顯，應(yīng)通過(guò)合理優(yōu)化甘蔗種植布局，進(jìn)一步提高區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。

一是廣西甘蔗生產(chǎn)應(yīng)繼續(xù)向桂西南、桂中地區(qū)轉(zhuǎn)移，通過(guò)提升區(qū)域甘蔗生產(chǎn)效率優(yōu)勢(shì)來(lái)提高綜合生產(chǎn)比較優(yōu)勢(shì)。

二是加快甘蔗主產(chǎn)區(qū)農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)，形成高標(biāo)準(zhǔn)現(xiàn)代化農(nóng)業(yè)生產(chǎn)基地，推動(dòng)甘蔗全程機(jī)械化發(fā)展水平，全面提高區(qū)域甘蔗全要素生產(chǎn)率，從而提高區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)。

三是結(jié)合廣西區(qū)域甘蔗生產(chǎn)比較優(yōu)勢(shì)的分布特點(diǎn)，在崇左、來(lái)賓、南寧和柳州等甘蔗生產(chǎn)規(guī)模優(yōu)勢(shì)較明顯的區(qū)域，進(jìn)一步強(qiáng)化“雙高”甘蔗生產(chǎn)基地建設(shè)，實(shí)現(xiàn)跨區(qū)域技術(shù)共享、農(nóng)機(jī)服務(wù)，進(jìn)一步增強(qiáng)廣西區(qū)域甘蔗生產(chǎn)空間溢出效應(yīng)。