王奉林，郭力娜，葉林

(華北理工大學(xué) 礦業(yè)工程學(xué)院，河北 唐山 063210)

休閑農(nóng)業(yè)作為一種新型產(chǎn)業(yè)，興起于20世紀(jì)的意大利、奧地利等地區(qū)，隨后迅速在歐美等一些國家發(fā)展起來，具有良好的發(fā)展前景[1]。休閑農(nóng)業(yè)是指利用農(nóng)業(yè)自然生態(tài)環(huán)境資源，同時結(jié)合農(nóng)村特有的文化、生活、生產(chǎn)活動，來給人們提供旅游度假休閑，體驗農(nóng)村生產(chǎn)生活機會等休閑體驗的一種新興農(nóng)業(yè)經(jīng)濟形態(tài)[2]。隨著我國鄉(xiāng)村轉(zhuǎn)型發(fā)展以及新農(nóng)村建設(shè)的不斷推進[3]，休閑農(nóng)業(yè)逐漸發(fā)展為提升國民福祉、提高農(nóng)業(yè)復(fù)合效益的重要產(chǎn)業(yè)。特別是在《全國休閑農(nóng)業(yè)發(fā)展第十二個五年規(guī)劃》提出之后，發(fā)展休閑農(nóng)業(yè)成為轉(zhuǎn)變農(nóng)業(yè)發(fā)展方式，促進農(nóng)民就業(yè)與增收，并成為推進新農(nóng)村建設(shè)的重要扶持發(fā)展手段，因此這對于滿足目前我國城鄉(xiāng)居民日益增長的休閑消費需求也具有重要的意義。

隨著休閑農(nóng)業(yè)的快速發(fā)展，人們對休閑旅游認(rèn)識的不斷加深，有關(guān)休閑農(nóng)業(yè)的研究也不斷得到深化與完善。目前，國內(nèi)外學(xué)者主要就休閑農(nóng)業(yè)的發(fā)展模式及戰(zhàn)略[3-5]、發(fā)展?fàn)顩r及評價[6-8]、空間布局與優(yōu)化[9-11]等開展了一些定性和定量的研究，還有部分學(xué)者對休閑農(nóng)業(yè)的影響因素[12-15]進行了研究。這些研究進一步豐富了休閑農(nóng)業(yè)形成與發(fā)展的相關(guān)理論和技術(shù)方法，然而通過查閱文獻發(fā)現(xiàn)：從研究內(nèi)容來看，目前對于休閑農(nóng)業(yè)的研究多為基于地理空間維度的靜態(tài)分析，少有學(xué)者基于時間維度從發(fā)展的角度出發(fā)對休閑農(nóng)業(yè)發(fā)展規(guī)律及分布格局[16，17]進行研究。唐山作為一個具有百年歷史的礦業(yè)城市，經(jīng)濟發(fā)展程度一直位于河北省前列，隨著資源型城市轉(zhuǎn)型升級、生態(tài)文明以及鄉(xiāng)村振興進程的不斷推進。近年來，唐山市通過科學(xué)規(guī)劃與創(chuàng)新發(fā)展，既因地制宜開發(fā)，又注重區(qū)域特色和優(yōu)勢，使得休閑農(nóng)業(yè)發(fā)展迅速已具有一定規(guī)律和規(guī)模。

基于此，該項研究選取中小尺度，同時近幾年休閑農(nóng)業(yè)發(fā)展已初具規(guī)模的唐山市為研究區(qū)，以 GIS 為平臺，并選取適合的實驗方法，從時間維度視角出發(fā)研究休閑農(nóng)業(yè)園區(qū)分布的發(fā)展情況，為唐山市休閑農(nóng)業(yè)產(chǎn)業(yè)的優(yōu)化發(fā)展提供數(shù)據(jù)支撐及參考。

1 研究區(qū)概況及數(shù)據(jù)

1.1 研究區(qū)概況

唐山位于河北省東北部，地處渤海灣中心地帶，居于燕山南麓，地勢北高南低，自西、西北向東及東南趨向平緩，直至沿海。唐山市現(xiàn)轄7個市轄區(qū)、3個縣級市、4個縣、3個開發(fā)區(qū)和1個管理區(qū)，總面積為13 472 km2，見圖1。根據(jù)統(tǒng)計數(shù)據(jù)，2017年唐山市共接待海內(nèi)外旅游者5 602.96萬人次，比2016年增長25.07%；旅游總收入587.31億元，比2016年增長34.09%；2017年全市范圍A級以上景區(qū)37個，星級飯店47家。同時全市農(nóng)業(yè)、林果、畜牧、水產(chǎn)資源豐富，能很好地為休閑農(nóng)業(yè)的發(fā)展提供有力支持。

圖1 唐山市分區(qū)概況圖

1.2 數(shù)據(jù)來源及初步處理

休閑農(nóng)業(yè)園區(qū)數(shù)據(jù)：該項研究使用的休閑農(nóng)業(yè)園區(qū)數(shù)據(jù)是通過網(wǎng)絡(luò)調(diào)研，并經(jīng)網(wǎng)絡(luò)地圖定位和人工電話核實確定的，共計收集到從2002年到2018年的340個園區(qū)數(shù)據(jù)(截止2018年4月28日)。通過進一步剔除已停業(yè)或不確定準(zhǔn)確開始營業(yè)時間的園區(qū)，最終獲得162個有效園區(qū)數(shù)據(jù)(其中包括75個省級園區(qū)[18]，87個普通園區(qū))，見圖2。結(jié)合前人研究成果與休閑農(nóng)業(yè)實際發(fā)展情況，將園區(qū)類型劃分為休閑觀光、農(nóng)事體驗、農(nóng)業(yè)科技、特色村鎮(zhèn)和鄉(xiāng)村文化5個一級類和12個二級類[15]，具體分類及其在各區(qū)縣分布情況，見表1。

圖2 唐山市休閑農(nóng)業(yè)園區(qū)分布圖

表1 唐山市休閑農(nóng)業(yè)園區(qū)類型與數(shù)量統(tǒng)計匯總表

社會經(jīng)濟數(shù)據(jù)：研究中用到的2017年接待海內(nèi)外旅游者次數(shù)、2017年旅游總收入等社會經(jīng)濟數(shù)據(jù)，通過2013～2018年《唐山市統(tǒng)計年鑒》獲取。

研究區(qū)矢量數(shù)據(jù)：唐山市各區(qū)縣、河流和道路等矢量數(shù)據(jù)下載自資源環(huán)境數(shù)據(jù)云平臺網(wǎng)站。

2 研究方法

2.1 基尼系數(shù)

基尼系數(shù)在地理學(xué)中通常被用來度量研究對象在區(qū)域內(nèi)集聚程度。該研究主要利用其對唐山區(qū)域內(nèi)休閑農(nóng)業(yè)分布差異進行分析。其計算公式如下：

(1)

C=1-G

(2)

其中，G為基尼指數(shù)；C為園區(qū)分布均勻度；pi為第i個區(qū)域內(nèi)休閑農(nóng)業(yè)園區(qū)的數(shù)量在整個唐山所占的比重；n為研究區(qū)域總數(shù)。G的值在0～1內(nèi)，G取值同均勻度成反比，當(dāng)其越接近于0時，表明區(qū)域間分布越平均；當(dāng)其越接近于1時，表明區(qū)域間分布越集中。當(dāng)0

2.2 地理集中指數(shù)

地理集中指數(shù)與洛倫茲曲線結(jié)合使用，在地理學(xué)中，通常用來分析研究對象集聚程度，該項目用其來研究休閑農(nóng)業(yè)園區(qū)在唐山區(qū)域尺度上的集聚狀況。其計算公式如下：

(3)

其中，Q為休閑農(nóng)業(yè)地理集中指數(shù)；xi為第i個區(qū)域內(nèi)休閑農(nóng)業(yè)點數(shù)量；T為休閑農(nóng)業(yè)總數(shù)；n為唐山縣區(qū)總數(shù)。Q的取值在0～100之間，其值越大，休閑農(nóng)業(yè)在唐山全域分布集聚狀況越明顯；其值越小，休閑農(nóng)業(yè)分布越分散[15]。

2.3 最鄰近指數(shù)

最鄰近指數(shù)主要通過計算平均最鄰近比率(ANN)來研究點狀要素的空間分布形態(tài)，ANN計算公式如下：

(4)

式中，R0為每個休閑農(nóng)業(yè)園區(qū)與其最近園區(qū)之間距離的平均值；Rε為假設(shè)隨機模式下休閑農(nóng)業(yè)園區(qū)的期望平均距離；N為休閑農(nóng)業(yè)園區(qū)的總數(shù)；S為唐山市全區(qū)面積或者各區(qū)縣面積。根據(jù)計算結(jié)果，ANN0.5和ANN1.5時，分別為聚集分布和均勻分布；0.5

3 休閑農(nóng)業(yè)園區(qū)演變情況

3.1 園區(qū)數(shù)量變化情況及發(fā)展階段分析

隨著經(jīng)濟的飛速發(fā)展以及持續(xù)穩(wěn)定增長，現(xiàn)代農(nóng)業(yè)的逐步深化，同時在市場多元化需求和政府政策的驅(qū)動下，唐山市休閑農(nóng)業(yè)依托其優(yōu)越區(qū)位條件和農(nóng)業(yè)產(chǎn)業(yè)資源，不斷涌現(xiàn)、發(fā)展，規(guī)模也持續(xù)壯大，如表2和圖3所示。

表2 2002～2018年唐山市休閑農(nóng)業(yè)增長情況統(tǒng)計結(jié)果

圖3 2002～2018年唐山市休閑農(nóng)業(yè)增長情況統(tǒng)計結(jié)果

由表2和圖3可知，唐山市休閑農(nóng)業(yè)在2004年以前由于園區(qū)數(shù)量基數(shù)較少導(dǎo)致其增長率過大，達到200.00%和100.00%；從2005年開始園區(qū)數(shù)量穩(wěn)步增長，增長率減小，逐漸趨于正常并且波動不大；2011年園區(qū)數(shù)量迅速增加，其增長率出現(xiàn)1個小峰值達44.74%，其后園區(qū)數(shù)量增長減緩，2014年園區(qū)增長數(shù)量最大，達到34個，同時增長率出現(xiàn)第2個小峰值達38.64%，這之后園區(qū)數(shù)量增長減少；總體來看，園區(qū)增長率從2002年到2018年該研究截止日，呈波動減小趨勢；園區(qū)數(shù)量在2005年之前呈穩(wěn)步增長趨勢，2011年開始迅速增長，到2014年園區(qū)數(shù)量增長達到最大值后，增長減少并逐漸趨于穩(wěn)定。

由于唐山市各區(qū)縣存在地理環(huán)境、自然資源及經(jīng)濟條件的差異，各類型休閑農(nóng)業(yè)園區(qū)在不同時間段內(nèi)發(fā)展情況也不同，研究參考2002～2018年唐山市休閑農(nóng)業(yè)增長情況，進一步對園區(qū)進行分類、分時間段匯總對比分析得到表3和圖4。

表3 唐山市各類型園區(qū)不同時間段統(tǒng)計匯總結(jié)果

圖4 唐山市各類型園區(qū)不同時間段統(tǒng)計結(jié)果

結(jié)合2002～2018年唐山市休閑農(nóng)業(yè)具體增長情況，并結(jié)合服務(wù)類型以及表3和圖4，將唐山市休閑農(nóng)業(yè)園區(qū)發(fā)展情況劃分為：萌芽階段(2004年及以前)，初始及起步階段(2005～2007年和2008～2010年)，品質(zhì)打造及快速發(fā)展階段(2011～2013年和2014～2016年)以及初步穩(wěn)定階段(2017年到研究截止期)等4個發(fā)展階段。

3.2 分布均衡性演變分析

由于不同的地理環(huán)境、自然資源及經(jīng)濟條件的差異，導(dǎo)致唐山市休閑農(nóng)業(yè)園區(qū)在不同時間段內(nèi)發(fā)展情況不同，園區(qū)分布的均衡性也有所差異，該研究采用基尼系數(shù)和地理集中指數(shù)2種對園區(qū)分布均衡性進行分時間段分析。

根據(jù)基尼系數(shù)計算公式，對唐山市各時間段休閑農(nóng)業(yè)園分布的基尼系數(shù)進行計算，得表4。

表4 唐山市各時間段內(nèi)休閑農(nóng)業(yè)園區(qū)的基尼系數(shù)

由表4可知，唐山市各時間段休閑農(nóng)業(yè)園區(qū)的基尼系數(shù)值最低為0.646(G>0.5)，雖然各時間段內(nèi)園區(qū)在整體上存在一定差異，但其基尼系數(shù)均大于0.6，分布均勻度均小于0.3，從整體發(fā)展來看基尼系數(shù)呈增大趨勢，分布均勻度呈減小趨勢，說明唐山市休閑農(nóng)業(yè)園區(qū)分布均衡性差異較大，園區(qū)集中程度較高，且隨時間推移園區(qū)分布的均衡性差異增大，各園區(qū)集中程度增大。

進一步利用公式(3)求得各時間段地理集中指數(shù)(實際指數(shù))，同時計算出各時間段的平均指數(shù)，詳見表5。以各區(qū)縣累積個數(shù)為橫軸，并園區(qū)累計比重為縱軸得到園區(qū)分布洛倫茲曲線圖，曲線的下凸程度，可以反映休閑農(nóng)業(yè)園區(qū)布局的集中程度，下凸程度越大，說明分布越集中，因此可以借助洛倫茲曲線來直觀地反映空間分布的不均衡程度在各時間段內(nèi)的發(fā)展情況，如圖5所示。

表5 唐山市各時間段內(nèi)休閑農(nóng)業(yè)園區(qū)的地理集中指數(shù)

圖5 唐山市各時間段內(nèi)休閑農(nóng)業(yè)園區(qū)洛倫茲曲線

由表5和圖5可知，各時間段內(nèi)得到的實際指數(shù)均大于平均分布的指數(shù)，說明各時間段內(nèi)園區(qū)在唐山全區(qū)域內(nèi)實際分布均為較集中的狀態(tài)，且從圖中6個時間段園區(qū)均衡度發(fā)展來看，地理集中指數(shù)(Q)呈減小趨勢，園區(qū)在唐山全區(qū)域內(nèi)集中分布程度減小，區(qū)域間極化程度削弱，差異有所縮小，全市園區(qū)逐漸從個別區(qū)縣集中分布發(fā)展，到呈全域發(fā)展。

3.3 聚集度演變分析

根據(jù)最鄰近分析法，對唐山市各時間段內(nèi)休閑農(nóng)業(yè)ANN值進行計算，計算結(jié)果見表6。根據(jù)表6計算結(jié)果可以看出，目前(即2017以后時間段)全市休閑農(nóng)業(yè)總體上呈“聚集—隨機型”分布模式。從各時間段對比分析來看，可以得出隨著時間推移唐山市休閑農(nóng)業(yè)園區(qū)最鄰近指數(shù)在逐漸減小，分布模式呈現(xiàn)出“隨機—離散型分布，隨機型分布，聚集—隨機型分布”逐漸向聚集型發(fā)展的趨勢，并且可以看出園區(qū)已初步實現(xiàn)局部聚集發(fā)展模式，說明唐山市休閑農(nóng)業(yè)已進入初步發(fā)展階段。

表6 唐山市休閑農(nóng)業(yè)園區(qū)各時間段內(nèi)空間分布特征

4 結(jié)論

(1)從園區(qū)發(fā)展角度來看，從2002～2018年唐山市休閑農(nóng)業(yè)園區(qū)可以劃分為，萌芽階段(2004年及以前)，初始及起步階段(2005～2007年和2008～2010年)，品質(zhì)打造及快速發(fā)展階段(2011～2013年和2014～2016年)以及初步穩(wěn)定階段(2017年到研究截止期)等4個發(fā)展階段，6個時間段。

(2)從均衡性角度來看，6個時間段內(nèi)休閑農(nóng)業(yè)園區(qū)分布均衡性呈增大趨勢，園區(qū)集中分布程度也呈增大趨勢；6個時間段內(nèi)唐山市全域內(nèi)休閑農(nóng)業(yè)園區(qū)分布均衡性呈減小趨勢，區(qū)域間極化程度削弱，差異有所縮小，說明園區(qū)從個別區(qū)縣分布情況發(fā)展到全域分布發(fā)展，園區(qū)在全域內(nèi)分布均衡性減小。

(3)從聚集度的角度來看，6個時間段內(nèi)隨著時間發(fā)展推移唐山市休閑農(nóng)業(yè)園區(qū)最鄰近指數(shù)在逐漸減小，分布模式呈現(xiàn)出“隨機—離散型分布，隨機型分布，聚集—隨機型分布”逐漸向聚集型發(fā)展的趨勢，最終表現(xiàn)為“聚集—隨機型”分布的初步“聚集”分布形態(tài)。

(4)從整體來看，唐山市全域內(nèi)園區(qū)分布呈全域發(fā)展趨勢，整個區(qū)域內(nèi)園區(qū)分布集中發(fā)展程度減小；園區(qū)間集中程度增大，園區(qū)逐漸向聚集型發(fā)展。