山東蔬菜生產(chǎn)格局演變及空間集聚效應(yīng)分析

苗曉穎，胡繼連，2，王秀鵑

（1.山東農(nóng)業(yè)大學(xué)經(jīng)濟管理學(xué)院，山東 泰安 271018；2.山東農(nóng)業(yè)大學(xué)馬克思主義學(xué)院，山東 泰安 271018；3.山東農(nóng)業(yè)大學(xué)公共管理學(xué)院，山東 泰安 271018）

蔬菜產(chǎn)業(yè)是我國傳統(tǒng)的農(nóng)業(yè)產(chǎn)業(yè)，是推動農(nóng)村經(jīng)濟發(fā)展的重要支柱產(chǎn)業(yè)。黨的“十九大”提出實施“鄉(xiāng)村振興”戰(zhàn)略，指出農(nóng)業(yè)是國民經(jīng)濟的基礎(chǔ)，農(nóng)村經(jīng)濟是現(xiàn)代化經(jīng)濟體系的重要組成部分。蔬菜作為經(jīng)濟作物，在農(nóng)村區(qū)域經(jīng)濟發(fā)展中占有舉足輕重的地位，發(fā)展蔬菜產(chǎn)業(yè)對推動農(nóng)村經(jīng)濟發(fā)展，促進“鄉(xiāng)村振興”戰(zhàn)略落實具有重要意義。我國蔬菜加工程度較低，主要仍以直接消費為主，蔬菜空間布局在一定程度上影響著區(qū)域?qū)用娴氖卟斯┙o。

近年來，我國學(xué)者對農(nóng)業(yè)生產(chǎn)布局的研究主要集中在農(nóng)業(yè)生產(chǎn)重心演變及其驅(qū)動因素、生產(chǎn)區(qū)域布局變動特征三個方面。丁存振等［1］運用重心模型，對羊肉生產(chǎn)布局演變進行分析，指出1980—2014年我國羊肉生產(chǎn)重心的移動方向、距離和速度。林正雨等［2］運用重心模型和ESDA方法對四川茶葉生產(chǎn)重心演變及空間格局進行分析。楊宗輝等［3］基于空間計量理論，通過對我國玉米空間布局的實證研究，發(fā)現(xiàn)農(nóng)業(yè)技術(shù)和政策對玉米空間溢出效應(yīng)顯著，各影響因素的局部性效應(yīng)使玉米空間格局發(fā)生變化。

對于蔬菜產(chǎn)業(yè)布局聚集程度及影響因素的研究，已形成成熟的研究方法。吳建寨等［4］對1995—2012年我國蔬菜生產(chǎn)空間集聚特征、效應(yīng)與機制進行研究，通過集中度指數(shù)、基尼系數(shù)及柯布－道格拉斯生產(chǎn)函數(shù)等方法，揭示了蔬菜生產(chǎn)聚集對產(chǎn)業(yè)發(fā)展的影響效應(yīng)。彭暉等［5］運用空間面板計量模型對中國蔬菜生產(chǎn)集聚度影響因素和空間相關(guān)性進行實證研究，認為技術(shù)條件、市場需求規(guī)模、自然資源及交通運輸成本是影響蔬菜集聚的主要因素。張哲晰等［6］利用空間Durbin模型、兩部制模型、地理加權(quán)回歸模型等方法，以蔬菜產(chǎn)業(yè)為例，研究了產(chǎn)業(yè)聚集對農(nóng)民收入的影響。張倩等［7］運用莫蘭指數(shù)和莫蘭散點圖從縣域維度對京津冀地區(qū)蔬菜生產(chǎn)相關(guān)性進行了分析。許偉［8］通過構(gòu)建蔬菜生產(chǎn)空間優(yōu)勢度評價指標體系，采用多因素綜合評價、空間相關(guān)性分析等方法對上海市崇明區(qū)蔬菜生產(chǎn)空間布局進行了研究。

目前，針對蔬菜生產(chǎn)格局演變及集聚的研究多從國家層面進行，以省為研究對象的較少。山東是我國傳統(tǒng)蔬菜產(chǎn)區(qū)，蔬菜播種面積和產(chǎn)量一直處于我國前列，但近些年有下降趨勢。開展山東蔬菜生產(chǎn)格局演變及集聚效應(yīng)研究，探索山東蔬菜空間布局變動的驅(qū)動力，有利于根據(jù)山東各地區(qū)比較優(yōu)勢，優(yōu)化山東蔬菜生產(chǎn)空間布局，給缺水地區(qū)農(nóng)業(yè)結(jié)構(gòu)調(diào)整提供樣板，在全國范圍內(nèi)推廣以農(nóng)業(yè)結(jié)構(gòu)調(diào)整為手段的農(nóng)業(yè)節(jié)水模式。

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

1.1 研究方法

1.1.1 重心模型 農(nóng)業(yè)產(chǎn)業(yè)重心模型源于力學(xué)重心模型，最早由美國學(xué)者Walker［9］運用到人口空間分布的研究中，隨后被用于經(jīng)濟重心分析、環(huán)境問題研究及農(nóng)業(yè)產(chǎn)業(yè)重心移動等多個方面，由于產(chǎn)業(yè)經(jīng)濟發(fā)展水平和速度不同，產(chǎn)業(yè)重心也處在不斷變化當中。本文將重心模型引入到蔬菜產(chǎn)業(yè)生產(chǎn)布局的研究中，選取蔬菜產(chǎn)量、蔬菜種植面積兩個指標，運用重心模型對山東省蔬菜生產(chǎn)重心進行研究，對山東省蔬菜產(chǎn)業(yè)生產(chǎn)格局空間演變特征進行分析。計算公式為:

式中，CXj、CYj分別表示第j年待評價指標重心坐標的經(jīng)度值和緯度值，cxi、cyi分別表示山東省i市（區(qū)）幾何中心的經(jīng)度值和緯度值，zij表示第j年i市（區(qū)）某一指標的量值。在本文中，評價指標為蔬菜產(chǎn)量和蔬菜種植面積，用各市（區(qū)）的幾何中心作為行政中心的經(jīng)緯度。

通過不同時期蔬菜生產(chǎn)的重心坐標可以計算出蔬菜生產(chǎn)重心空間移動距離，具體公式為:

式中，d(t＋k)－t表示t到t＋k期蔬菜生產(chǎn)重心移動距離，λ為常數(shù)，是地球表面坐標單位（°）轉(zhuǎn)化為平面距離（km）的系數(shù)，取值為111.11。

1.1.2 區(qū)位商和產(chǎn)業(yè)集中度指數(shù)CRn產(chǎn)業(yè)集中度常用的評價指標有區(qū)位商和產(chǎn)業(yè)集中度指數(shù)CRn，區(qū)位商可以測算一個產(chǎn)業(yè)的相對集中度，產(chǎn)業(yè)集中度指數(shù)CRn可以測算產(chǎn)業(yè)的絕對集中度。在本文中，使用修正后的區(qū)位商和產(chǎn)業(yè)集中度指數(shù)CRn分析中國蔬菜生產(chǎn)的集聚特征。

修正后的區(qū)位商（FLQ*）是傳統(tǒng)區(qū)位商（LQ）的修正和改進，其表達式為:

式中，λ*為修正系數(shù)；根據(jù)敏感性檢驗，本文中δ取0.3；FLQ*的值越高，表明該地區(qū)蔬菜生產(chǎn)的專業(yè)化程度越高。eij、Ej分別為i市和全省蔬菜產(chǎn)值，ei為i市農(nóng)業(yè)產(chǎn)值，∑Ej為全省農(nóng)業(yè)產(chǎn)值。

產(chǎn)業(yè)集中度指數(shù)CRn指某產(chǎn)業(yè)相關(guān)市場內(nèi)規(guī)模最大的前n個地區(qū)或企業(yè)所占的市場份額，可以反映產(chǎn)業(yè)的集中程度，在本文中使用產(chǎn)業(yè)集中度指數(shù)來測度蔬菜生產(chǎn)的絕對集中度。

1.1.3 探索性空間數(shù)據(jù)分析方法 探索性空間數(shù)據(jù)分析方法（ESDA）是對某一地理變量空間分布中相鄰位置間相關(guān)性進行檢驗的一種統(tǒng)計方法，利用該方法可以揭示區(qū)域變化量的空間分布特征，包括全局自相關(guān)分析和局部自相關(guān)分析。

全局自相關(guān)指標Moran’s I主要探索屬性數(shù)值在整個區(qū)域分布上空間趨同特征強弱。其計算公式為:

局部空間自相關(guān)指標LISA用于揭示一個區(qū)域單元上的某種地理現(xiàn)象或某一屬性值與鄰近單元上同一地理現(xiàn)象或?qū)傩灾档南嚓P(guān)程度。其計算公式為:

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

研究數(shù)據(jù)包括屬性數(shù)據(jù)和空間數(shù)據(jù)，空間數(shù)據(jù)來源于國家基礎(chǔ)地理信息數(shù)據(jù)庫的山東省縣級行政區(qū)劃圖，屬性數(shù)據(jù)來源于《山東統(tǒng)計年鑒》及各縣（市、區(qū)）主要經(jīng)濟指標。

2 山東省蔬菜生產(chǎn)概況

山東省位于我國東部，黃河下游，大體分為魯中（濟南、淄博、濰坊、泰安、濱州鄒平）、魯南（菏澤、濟寧、臨沂、日照、棗莊）、魯西（聊城、德州的南部）、魯北（東營、濱州、濰坊的西北部、德州的北部）、魯東（青島、煙臺、威海）地區(qū)五個板塊。全省土地總面積1 571.26萬hm2，其中農(nóng)用地1 156.60萬hm2，占土地總面積的73.61%，為傳統(tǒng)農(nóng)業(yè)大省。全省年均氣溫11～14℃，年均降水量550～950 mm，由東南向西北遞減，降水季節(jié)分配不均，多集中在6～8月份，年均日照2 290～2 890 h，適宜發(fā)展優(yōu)質(zhì)蔬菜生產(chǎn)。

山東蔬菜種植歷史悠久，在產(chǎn)量、播種面積、產(chǎn)值、出口方面均居全國前列。2019年全省蔬菜播種面積達146.42萬hm2，僅低于河南省，位居全國第二，產(chǎn)量為8 181.15萬t，雖然自2016年有所下降，但仍是全國蔬菜產(chǎn)量最高的省區(qū)。蔬菜出口量連續(xù)15年居全國首位，2019年蔬菜出口額288.5億元。

如表1所示，1995—2015年山東蔬菜產(chǎn)量逐年遞增，但增幅有所放緩，2016年后有所下降，與全國蔬菜產(chǎn)量變動趨勢基本一致，可以看出，山東蔬菜產(chǎn)量的規(guī)模變化在一定程度上影響著我國蔬菜產(chǎn)量。除1996、2000、2002、2007年外，山東蔬菜產(chǎn)量增長率均低于全國蔬菜產(chǎn)量增長率，近年來山東蔬菜產(chǎn)量占全國蔬菜產(chǎn)量的比重有所下降。自1995年，山東蔬菜產(chǎn)量變動大致可以分為四個階段，第一階段為1995—2004年的快速增長期，蔬菜產(chǎn)量一直保持穩(wěn)定增長，且增長率較大；第二階段為2005—2007年的緩慢回落期，蔬菜產(chǎn)量有所下降；第三階段為2008—2015年的穩(wěn)定上升期，2015年蔬菜產(chǎn)量達到頂峰，但增長率低于全國；第四階段為2016—2019年的回落穩(wěn)定期，蔬菜產(chǎn)量于2016年出現(xiàn)了斷崖式下降，隨后產(chǎn)量基本保持不變。

表1 1995—2019年山東蔬菜生產(chǎn)變化

從播種面積來看，山東蔬菜播種面積經(jīng)歷了先上升后下降的過程，這一過程與農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整有關(guān)。2019年山東產(chǎn)出全國11.35%的蔬菜，但播種面積僅占全國7.02%，可見，山東仍是我國主要蔬菜產(chǎn)區(qū)且蔬菜生產(chǎn)效率高。

3 山東蔬菜生產(chǎn)重心演變路徑分析

3.1 山東蔬菜播種面積重心演變路徑分析

測算結(jié)果顯示，1995—2019年山東蔬菜播種面積重心在東經(jīng)117°28′04″～117°50′28″、北緯36°03′29″～36°15′42″之間波動。

運用蔬菜播種面積生產(chǎn)重心坐標可以得到山東蔬菜播種面積歷年重心演變路徑（圖1）。從總體上來看，1995—2019年山東蔬菜播種面積重心演變波動性較大，移動過程中出現(xiàn)了1997、2000、2008、2014年等重要拐點。從重心的空間分布集中情況來看，蔬菜播種面積重心在1995—1999年階段波動較大，除1997年外，重心位移距離均在萬米以上，其中1996年重心位移距離最大，達32 269.65 m；2000—2004年為蔬菜播種面積的調(diào)整過度階段，山東蔬菜播種面積重心穩(wěn)定向東南方向移動；在2005—2019年階段重心分布較為集中，波動較為平穩(wěn)。從演變的方向來看，1995—1997年山東蔬菜播種面積重心主要向西南方向移動，位移變化大，1998—2000年山東蔬菜播種面積重心主要向東北方向移動，位移變化較大，2001—2005年山東蔬菜播種面積重心主要向西南方向移動，位移變化不大且比較平穩(wěn)，2006—2019年山東蔬菜種植面積重心呈小幅度不規(guī)律變動，位移變化較小，其中2006—2008、2013—2019年間主要向南移動，2009—2012年間主要向北移動。

圖1 山東蔬菜播種面積重心歷年演變路徑（1995—2019）

3.2 山東蔬菜產(chǎn)量重心演變路徑分析

測算結(jié)果顯示，1995—2019年山東蔬菜產(chǎn)量重心 在 東 經(jīng)117°36′42″～118°03′54″、北 緯36°08′12″～36°15′56″之間波動，波動多集中于東經(jīng)117°36′42″～117°45′45″之間。

運用蔬菜產(chǎn)量生產(chǎn)重心坐標可以得到山東蔬菜產(chǎn)量歷年重心演變路徑（圖2）。從總體上來看，1995—2019年山東蔬菜產(chǎn)量重心演變波動性較大，移動過程中出現(xiàn)了1997、1999、2008、2014年等重要拐點。從重心的空間分布集中情況來看，在1995—1998年階段波動較大，重心位移距離均在萬米以上，其中1996年重心位移距離最大，達32 578.20 m；1999—2006年為山東蔬菜產(chǎn)量的調(diào)整過度階段，重心隨著蔬菜種植面積重心向西南方向移動；在2007—2019年階段重心分布較為集中，移動距離較小，波動較為平穩(wěn)。從演變的方向來看，1995—1997年山東蔬菜產(chǎn)量重心主要向西南方向移動，位移變化大，1998—2000年山東蔬菜產(chǎn)量重心主要向東北方向移動，位移變化較大，2001—2008年山東蔬菜產(chǎn)量重心主要向西南方向移動，位移變化不大且比較平穩(wěn)，2009—2019年山東蔬菜產(chǎn)量重心呈小幅度不規(guī)律變動，位移變化較小，其中2009—2014年主要向西北方向移動，2015—2019年主要向東南方向移動。

圖2 山東蔬菜產(chǎn)量重心歷年演變路徑（1995—2019）

4 山東蔬菜生產(chǎn)的空間集中度

4.1 山東蔬菜生產(chǎn)的地區(qū)專業(yè)化分析

蔬菜生產(chǎn)的專業(yè)化程度可以反映各地區(qū)的農(nóng)業(yè)分工情況和農(nóng)業(yè)結(jié)構(gòu)差異，在一定程度上也可以反映各地區(qū)農(nóng)業(yè)產(chǎn)業(yè)中的主導(dǎo)產(chǎn)業(yè)和發(fā)展方向。一般來說，區(qū)位商越高的地區(qū)，其蔬菜生產(chǎn)專業(yè)化程度越高，自然資源、勞動資源或社會經(jīng)濟基礎(chǔ)的優(yōu)勢越強，蔬菜產(chǎn)業(yè)規(guī)模較大，生產(chǎn)集中度較高。

2000—2015年，山東蔬菜生產(chǎn)專業(yè)化水平總體穩(wěn)步上升，2019年農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)化調(diào)整后，專業(yè)化水平出現(xiàn)下降趨勢。自1995年以來，蔬菜生產(chǎn)專業(yè)化水平一直處在山東市區(qū)前六名的為棗莊和濰坊（表2），泰安在2015年前位于前六名，2019年下滑至第九名。就變化趨勢來說，專業(yè)化水平靠前的市區(qū)由分散分布轉(zhuǎn)變?yōu)橄蛭鞑亢湍喜康貐^(qū)聚集的趨勢。1995年蔬菜專業(yè)化水平前六名的城市在魯中、魯東、魯南地區(qū)均有分布；2019年演變?yōu)橹饕性趦蓚€區(qū)域，即魯南地區(qū)（棗莊、菏澤、臨沂）和魯西地區(qū)（聊城、德州），魯中地區(qū)泰安、濟南、萊蕪（現(xiàn)濟南市萊蕪區(qū)和鋼城區(qū)）專業(yè)化水平大幅下降，僅濰坊市處在排名前六的隊列中。魯南地區(qū)專業(yè)化生產(chǎn)水平提升較快，菏澤專業(yè)化水平穩(wěn)定提升，2019年僅位于濰坊、棗莊之后，臨沂專業(yè)化水平也有上升趨勢，2019年位于第六名。魯西地區(qū)的聊城市、德州市專業(yè)化水平在2000年后一直處于較高的狀態(tài)。

表2 山東蔬菜生產(chǎn)專業(yè)化水平前6個市（區(qū)）的區(qū)位商

4.2 山東蔬菜生產(chǎn)的地理聚集分析

山東蔬菜的主產(chǎn)區(qū)（本文將蔬菜生產(chǎn)占全省蔬菜生產(chǎn)比重高于7%的市（區(qū)）定義為山東蔬菜主產(chǎn)區(qū)）主要位于魯西地區(qū)和魯中地區(qū)，這兩大區(qū)域蔬菜種植的面積和產(chǎn)量占全省面積和產(chǎn)量的比重均高于60%，所以山東蔬菜生產(chǎn)具有地理區(qū)位上的聚集性。從生產(chǎn)的絕對集中度來看（表3），山東蔬菜生產(chǎn)的集中度一直處于較高的水平，2019年，蔬菜種植面積和產(chǎn)量的CR8分別為76.89%和76.44%，比1995年增加了10.90%和11.59%，可以看出大部分蔬菜主產(chǎn)區(qū)保持著較快的增長速度，且播種面積的增長速度與產(chǎn)量的增長速度基本保持一致；蔬菜種植面積和產(chǎn)量的CR6分別為63.01%和61.02%，比1995年增加了13.10%和9.61%，播種面積的增長速度要大于產(chǎn)量的增長速度；蔬菜種植面積和產(chǎn)量的CR4分別為46.38%和44.36%，比1995年增加了16.30%和5.17%，播種面積的增長速度遠大于產(chǎn)量的增長速度。可見，山東蔬菜生產(chǎn)有集中分布的趨勢，但產(chǎn)量集中度要小于播種面積集中度，這也反映了山東蔬菜在空間布局演變過程中，傳統(tǒng)優(yōu)勢產(chǎn)區(qū)由于配套設(shè)施健全、技術(shù)成熟，蔬菜生產(chǎn)效率較高，產(chǎn)量增長速度較快，帶動其周邊地區(qū)蔬菜種植規(guī)模不斷擴大，但由于勞動力、物質(zhì)投資等資源條件的限制，產(chǎn)量增長速度慢于播種面積增長速度，造成了一定的自然資源浪費。但無論是在播種面積上還是產(chǎn)量上，山東蔬菜生產(chǎn)呈現(xiàn)空間上的連片分布特征，1995年以濰坊為中心的魯中地區(qū)是山東蔬菜生產(chǎn)重心，但隨著生產(chǎn)布局演變，2019年，濟寧市、聊城市、菏澤市及其周邊地區(qū)成為魯西地區(qū)和魯南地區(qū)的蔬菜集聚中心。

表3 山東蔬菜生產(chǎn)的空間集中度（%）

5 山東蔬菜生產(chǎn)空間相關(guān)性及驅(qū)動力分析

5.1 山東蔬菜生產(chǎn)全局自相關(guān)分析

本文用莫蘭指數(shù)來分析山東省蔬菜生產(chǎn)的全局空間相關(guān)性。對山東地區(qū)的137個縣（區(qū)）構(gòu)建鄰接矩陣，計算2000—2019年的Moran’s I和Z值，結(jié)果如表4所示。全局Moran’s I均為正值，Z值均大于1.96且P值均小于0.05，因此可以拒絕原假設(shè)，山東省蔬菜生產(chǎn)的空間分布不是隨機的，表現(xiàn)出顯著的空間相關(guān)性，即山東省蔬菜生產(chǎn)具有空間聚集性。

表4 2000—2019年山東省縣域蔬菜產(chǎn)量的 全局空間自相關(guān)Moran’s I和Z值

在2000—2008年間，Moran’s I呈增長態(tài)勢，從0.198259增長至0.253316，表明在該期間，山東蔬菜生產(chǎn)聚集性逐漸增強，原因可能是隨著社會的不斷發(fā)展，科技水平的提高，蔬菜生產(chǎn)不再局限于自然條件優(yōu)越、有蔬菜種植歷史的傳統(tǒng)蔬菜產(chǎn)區(qū)，具體表現(xiàn)為，蔬菜產(chǎn)量較高的縣帶動周邊地區(qū)蔬菜生產(chǎn)的發(fā)展，即蔬菜產(chǎn)量較高的縣趨于和蔬菜產(chǎn)量較高的縣相鄰，蔬菜產(chǎn)量較低的縣趨于和蔬菜產(chǎn)量較低的縣相鄰。

在2008—2019年間，Moran’s I出現(xiàn)下降趨勢，從0.253316下降至0.184385，表明在2008年后，山東蔬菜生產(chǎn)聚集性在不斷減弱，原因可能是一方面當科技提高帶來的邊際收益逐漸減少時，理智的農(nóng)戶會放棄盲目追尋因周邊蔬菜生產(chǎn)的收益高而種植蔬菜，政府也會降低對蔬菜生產(chǎn)的財政投入，導(dǎo)致蔬菜的輻射帶動效應(yīng)較弱；另一方面，由于蔬菜屬于勞動力聚集產(chǎn)業(yè)和高耗水作物，勞動力資源與自然資源對蔬菜生產(chǎn)的影響力逐漸增強，蔬菜生產(chǎn)逐漸向勞動力和自然資源豐富的地區(qū)轉(zhuǎn)移，區(qū)位優(yōu)勢的影響相對較弱，聚集程度有所下降。

5.2 山東蔬菜生產(chǎn)局部自相關(guān)分析

全局空間相關(guān)性分析可以對整個研究區(qū)域的空間模式進行分析，但不能揭示不同地區(qū)空間關(guān)聯(lián)模式，不能確定聚集發(fā)生的具體位置，也不能確定聚集是屬于高－高聚集還是低－低聚集等聚集形式，所以對于具體空間區(qū)位和聚集關(guān)系需要進行局部自相關(guān)分析。

為了進一步分析各個區(qū)域與周邊區(qū)域的局部空間關(guān)聯(lián)、空間差異程度和空間格局分布，利用局部莫蘭指數(shù)和ArcGIS軟件生成2000、2005、2010、2015、2019年5個不同年份縣域蔬菜產(chǎn)量局部自相關(guān)（LISA）集聚圖（圖3），可以更直觀地觀察到山東省縣域蔬菜產(chǎn)量的空間聚集情況。

圖3 2000—2019年5個不同年份縣域蔬菜產(chǎn)量局部自相關(guān)LISA集聚

2000年以來山東省蔬菜生產(chǎn)縣域呈現(xiàn)出顯著的局部自相關(guān)關(guān)系，不同類型蔬菜產(chǎn)量縣域空間格局變動規(guī)律如下:

（1）蔬菜高產(chǎn)縣域聚集于魯中和魯西地區(qū)，魯南地區(qū)和魯東地區(qū)有少量分布。2000年山東省蔬菜高產(chǎn)縣域的高－高集聚區(qū)主要集中在魯中地區(qū)和魯東地區(qū)，分布在魯中地區(qū)以濰坊市壽光市為中心的濰坊市、青州市、廣饒縣以及章丘區(qū)、萊蕪市、岱岳區(qū)、肥城市以及魯東地區(qū)的平度市、即墨市、膠州市。相比2000年，2005年高－高集聚區(qū)逐步向魯西和魯南集中，在莘縣、嘉祥縣、金鄉(xiāng)縣形成新的集聚，魯東地區(qū)的高－高集聚區(qū)消失。2010年高－高集聚中心進一步向魯西地區(qū)集中，東昌府區(qū)、陽谷縣形成新的聚集，魯中和魯南地區(qū)的集聚區(qū)減少。2015年，除新增昌樂縣外，高－高集聚區(qū)沒有發(fā)生變化。2019年高－高集聚區(qū)進一步分散，形成魯西地區(qū)的莘縣、陽谷縣、冠縣集聚，魯南地區(qū)的金鄉(xiāng)縣、單縣集聚，魯中地區(qū)的岱岳區(qū)集聚和青州市、壽光市、昌樂縣集聚，魯東地區(qū)的平度市集聚五大蔬菜生產(chǎn)高－高集聚區(qū)。

從整體上看，高－高集聚區(qū)空間分布變化較大，分布越來越分散，但魯中和魯西地區(qū)始終是高－高集聚的主要分布區(qū)，主要原因一是魯中地區(qū)降水較多，與蔬菜高耗水的屬性吻合，魯西地區(qū)土壤有機質(zhì)高、酸堿度適宜，有大量的耕地，適合蔬菜的種植。二是魯中濰坊市和魯西聊城市地區(qū)蔬菜栽培歷史長，尤其是濰坊市，是全國著名的菜籃子。三是由于蔬菜生產(chǎn)的特殊性，由科技進步帶來的收益已接近頂峰，擴散效應(yīng)降低，聚集區(qū)逐漸分散。

（2）蔬菜生產(chǎn)高－低聚集類型和低－高聚集類型的縣域很少。2000、2005年，高－低類型集聚只有滕州市，其蔬菜產(chǎn)量遠高于周邊縣域，2010年后，高－低類型集聚新增蘭陵縣，均分布在魯南地區(qū)。低－高類型集聚僅在2000年出現(xiàn)在泰山區(qū)。其原因是滕州市、蘭陵縣作為山東的富水區(qū)，同時也有悠久的蔬菜種植歷史，蔬菜產(chǎn)量高于周邊縣域。而泰山區(qū)位于省會濟南市周邊，隨著濟南市經(jīng)濟發(fā)展，農(nóng)作物種植類型逐步多樣化，蔬菜種植的優(yōu)勢逐漸下降。

另外，山東省大部分縣域都屬于局部自相關(guān)關(guān)系不顯著的類型，說明山東省蔬菜生產(chǎn)水平仍有提升空間，目前僅有小部分縣域形成規(guī)?；N植。同時，山東省應(yīng)將蔬菜生產(chǎn)向精細化生產(chǎn)方向轉(zhuǎn)變，按照趨勢向資源豐富地區(qū)移動，充分突出各個地區(qū)的優(yōu)勢。

5.3 驅(qū)動力分析

5.3.1 自然條件 蔬菜是高耗水作物，蔬菜生產(chǎn)發(fā)展依賴于一定的自然地理環(huán)境，自然區(qū)位優(yōu)劣直接影響蔬菜產(chǎn)業(yè)的生產(chǎn)發(fā)展水平，進而造成了各縣域間蔬菜生產(chǎn)水平的空間差異。山東省蔬菜種植高產(chǎn)區(qū)在二十一世紀初具有明顯的集聚趨勢，主要集中在山東省中部和東部，但近些年，蔬菜種植高產(chǎn)聚集片區(qū)數(shù)目增多，分布越來越分散，在魯西、魯南、魯中、魯東地區(qū)均有分布。一是這些地區(qū)的土壤質(zhì)地多為潮土、棕壤、褐土，能提供蔬菜生長所需營養(yǎng)物質(zhì)，適合蔬菜生長。二是山東省主要農(nóng)業(yè)用水來源于大氣降水、地下水，這些地區(qū)大多有黃河客水資源流入或區(qū)域內(nèi)河流較多，光照充足，滿足蔬菜生長發(fā)育需要，而魯東地區(qū)由于地下水過度使用，海水滲透，水土條件惡化，蔬菜生產(chǎn)受限，其蔬菜生產(chǎn)逐漸向山東省西部和南部轉(zhuǎn)移。

5.3.2 社會經(jīng)濟發(fā)展水平 隨著農(nóng)業(yè)科技的發(fā)展，設(shè)施蔬菜逐漸成為山東蔬菜生產(chǎn)的主要模式，社會經(jīng)濟水平的提升有助于蔬菜產(chǎn)區(qū)農(nóng)業(yè)基礎(chǔ)設(shè)施的完善與蔬菜生產(chǎn)效率的提高，從而影響蔬菜產(chǎn)業(yè)的空間布局。尤其對山東蔬菜高產(chǎn)區(qū)的集聚與分散有明顯的影響。一方面，2010年前，科技創(chuàng)新使得一些自然條件較差地區(qū)的蔬菜生產(chǎn)條件變得相對容易，降低了蔬菜對自然條件的依賴性，如大棚技術(shù)、自動化技術(shù)是壽光蔬菜產(chǎn)業(yè)集群發(fā)展與擴張的重要條件。同時，交通、物流條件的改善，使蔬菜打破本地生產(chǎn)，跨地域銷售成為可能。另一方面，在經(jīng)濟報酬的拉動下，蔬菜生產(chǎn)者收益等經(jīng)濟因素對蔬菜空間布局起到了更大的作用。消費者經(jīng)濟水平的提升，引起消費習(xí)慣的改變，在經(jīng)濟發(fā)達地區(qū)，對高價蔬菜的接受能力變強，再加上政府補貼，即使經(jīng)濟發(fā)達地區(qū)種植蔬菜也變得有利可圖，再加上完全依靠外地蔬菜供應(yīng)體系脆弱性日益體現(xiàn)，使山東東部地區(qū)蔬菜種植出現(xiàn)“回遷”，這也是山東蔬菜生產(chǎn)由集聚變?yōu)榉稚⒌囊粋€重要原因。

5.3.3 政策因素 區(qū)域政策是影響山東省蔬菜產(chǎn)業(yè)發(fā)展最直接也最有效的因素。國家政策和地方調(diào)控政策都對山東蔬菜空間布局有重要影響。《山東省蔬菜產(chǎn)業(yè)提質(zhì)增效轉(zhuǎn)型升級實施方案（2016—2020年）》就是山東2015年后山東蔬菜空間布局變動的重要原因。在自然資源尤其是水資源在蔬菜生產(chǎn)中地位越來越重要的背景下，《山東省農(nóng)業(yè)水價綜合改革實施方案》及《山東省發(fā)展和改革委員會山東省水利廳關(guān)于建立完善農(nóng)業(yè)用水超定額累進加價制度的指導(dǎo)意見》中超定額累進加價制度使水資源成為影響蔬菜空間布局的重要因素。與其他因素相比，政策因素具有強制性、計劃性、可預(yù)見性的特點，能根據(jù)城市發(fā)展需要科學(xué)全面地引導(dǎo)產(chǎn)業(yè)發(fā)展。

6 結(jié)論與建議

6.1 主要結(jié)論

本文以山東省為例，選用1995—2019年山東省市域和縣域蔬菜相關(guān)數(shù)據(jù)，采用重心模型、ESDA等方法，對山東省蔬菜生產(chǎn)格局演變進行分析，并進一步分析了影響蔬菜生產(chǎn)空間格局演變的因素。得出以下主要結(jié)論:

（1）山東省蔬菜生產(chǎn)重心呈現(xiàn)向南部移動的趨勢。通過對比蔬菜播種面積重心與蔬菜產(chǎn)量重心兩個指標，可以得到其演變路徑的變化趨勢基本一致，即蔬菜播種面積重心與蔬菜產(chǎn)量重心在總體上均向山東省南部移動，相對于蔬菜播種面積重心演變，蔬菜產(chǎn)量重心移動方向更偏向西南一些。二者移動趨勢也基本一致，在1995—1997年大幅向西部移動，1998—2008年穩(wěn)定向南部移動，2009年后出現(xiàn)不規(guī)律小幅波動。這說明2008年之前山東蔬菜處于不斷調(diào)整狀態(tài)，2009年后形成穩(wěn)定的集聚群。

（2）山東省蔬菜生產(chǎn)地理集中度較高，專業(yè)化水平較低。山東省蔬菜主要分布在魯中地區(qū)和魯西地區(qū)，從絕對集中度測算結(jié)果來看，山東蔬菜生產(chǎn)表現(xiàn)為較高的集中程度，總體上排名前六的主產(chǎn)區(qū)蔬菜生產(chǎn)規(guī)模呈不斷擴大的態(tài)勢。但從修正的區(qū)位商測算結(jié)果來看，山東省蔬菜生產(chǎn)專業(yè)化程度仍處于一個比較低的水平。

（3）山東省蔬菜生產(chǎn)分布呈顯著空間集聚的特征。從空間自相關(guān)分析來看，全局自相關(guān)表明，山東蔬菜空間集聚性經(jīng)歷了先增強后減弱的過程。部分地區(qū)蔬菜生產(chǎn)分布呈現(xiàn)出局部空間集聚現(xiàn)象，高－高集聚區(qū)由集中變?yōu)榉稚⒎植迹?019年分布在山東省中部、西部和南部地區(qū)，高－低集聚和低－高集聚區(qū)較少。

6.2 建議

（1）發(fā)展特色蔬菜產(chǎn)業(yè)。從研究結(jié)果來看，山東蔬菜生產(chǎn)水平在各縣區(qū)之間存在顯著差異，一些地區(qū)并不適合種植蔬菜。建議充分發(fā)揮各個地區(qū)的自然資源優(yōu)勢和社會經(jīng)濟優(yōu)勢，根據(jù)自然條件規(guī)劃不同品種的蔬菜種植，在經(jīng)濟發(fā)達區(qū)域種植特色、優(yōu)質(zhì)蔬菜，從而提升市場競爭力，建設(shè)一批符合區(qū)域特色的蔬菜生產(chǎn)基地。

（2）完善基礎(chǔ)設(shè)施，加大科技投入。山東蔬菜生產(chǎn)空間相關(guān)性處于減弱態(tài)勢且高－高集聚區(qū)逐漸分散，要實現(xiàn)山東蔬菜產(chǎn)業(yè)持續(xù)發(fā)展，就要完善各項配套設(shè)施，加大科技投入。一方面，要完善灌溉工程、施肥工程等相關(guān)基礎(chǔ)設(shè)施的建設(shè)，節(jié)約水資源，解放勞動力，提升蔬菜生產(chǎn)效益；另一方面，加大科研投入，根據(jù)各地不同的自然環(huán)境條件，培育新型蔬菜品種，增強蔬菜抵御自然災(zāi)害和病蟲害的能力，不斷提升蔬菜品質(zhì)。

（3）優(yōu)化蔬菜生產(chǎn)布局。目前山東省蔬菜生產(chǎn)仍主要分布在魯中和魯西地區(qū)，與水資源稟賦存在矛盾，要進一步推動蔬菜生產(chǎn)向豐水市域移動，建設(shè)一批與自然資源稟賦匹配的優(yōu)勢蔬菜生產(chǎn)基地，例如可將魯西缺水地區(qū)的蔬菜逐步向魯南地區(qū)轉(zhuǎn)移，逐步建立以臨沂為中心的魯南蔬菜生產(chǎn)基地。同時，要控制魯西缺水市域的蔬菜種植規(guī)模，推廣應(yīng)用農(nóng)業(yè)節(jié)水技術(shù)，使蔬菜布局更合理，促進山東蔬菜產(chǎn)業(yè)可持續(xù)發(fā)展。