[摘要]以廣西百色市為例，在 GIS及 RS技術(shù)的支持下，結(jié)合回歸分析方法，對研究年份的TM遙感圖像進(jìn)行解譯分類，利用土地利用動態(tài)度、轉(zhuǎn)移矩陣分析模型、回歸擬合分析等方法對土地利用變化進(jìn)行研究，揭示百色市新型城鎮(zhèn)化水平與土地利用變化的關(guān)系。結(jié)果表明，①2005 ～2015百色市各類土地的面積都在變動;②百色市新型城鎮(zhèn)化水平的高低與各類土地面積存在著相關(guān)性;③土地城鎮(zhèn)化、產(chǎn)業(yè)城鎮(zhèn)化及人口城鎮(zhèn)化之間存在著正反饋關(guān)系。

[關(guān)鍵詞]土地利用變化;新型城鎮(zhèn)化;廣西百色市

[中圖分類號]F301.24[文獻(xiàn)標(biāo)識碼]A

土地作為人類活動的空間載體，在城鎮(zhèn)化進(jìn)程中發(fā)揮著極其重要的作用。隨著我國經(jīng)濟(jì)的逐步發(fā)展， “空間城鎮(zhèn)化”并沒有帶來“人口城市化”。與此同時，土地撂荒、人地矛盾、優(yōu)質(zhì)耕地漏失等問題愈演愈烈，我國城鎮(zhèn)化的模式亟待轉(zhuǎn)型。隨著相關(guān)規(guī)定的出臺以及有關(guān)新型城鎮(zhèn)化試點方案的推行，使得新型城鎮(zhèn)化的概念深入人心。在新型城鎮(zhèn)化發(fā)展的進(jìn)程中國內(nèi)外的學(xué)者開始關(guān)注新型城鎮(zhèn)化與土地利用變化之間的關(guān)系，主要集中在新型城鎮(zhèn)化中土地利用中的結(jié)構(gòu)與布局、效益耦合發(fā)展、集約利用程度等方面的研究。也有一些學(xué)者利用3S技術(shù)結(jié)合回歸分析法，對究土地利用變化與城市化及土壤侵蝕程度進(jìn)行研究。

縱觀國內(nèi)外，學(xué)者們對僅是對一些發(fā)達(dá)城市的土地利用變化與城鎮(zhèn)化關(guān)系進(jìn)行研究，對邊境地區(qū)研究較少;且多數(shù)從土地利用效益、土地利用結(jié)構(gòu)的角度對城鎮(zhèn)化進(jìn)行剖析，很少從土地利用變化這個大的角度研究其與新型城鎮(zhèn)化發(fā)展的關(guān)系。且這些地區(qū)基本上都屬于廣西區(qū)外的研究對廣西區(qū)內(nèi)城市的研究極少，特別是對對革命老區(qū)的研究更少。因此，以廣西百色市為例研究2005～2015百色市土地利用的變化趨勢，并探討新型城鎮(zhèn)化前后的其變化情況，揭示新型城鎮(zhèn)化與土地利用變化的內(nèi)在聯(lián)系。

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

1.1 研究區(qū)域概況

百色市位于廣西西部及右江流域，旅游與礦產(chǎn)資源豐富，隨著新型城鎮(zhèn)化的推進(jìn)，百色市的經(jīng)濟(jì)與城鎮(zhèn)化發(fā)展顯著加快。2015年百色市總?cè)丝谶_(dá)413.19萬人，其中常住人口占總數(shù)的87%，整個市的GDP達(dá)980.4億元。由于山區(qū)面積占比較大，城鎮(zhèn)化的快速發(fā)展對土地資源稀缺的百色市提出來更高的要求，迫切需要在新型城鎮(zhèn)化過程中提高土地的利用率，形成更為合理的發(fā)展模式。

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

通過互聯(lián)網(wǎng)資源共享下載2005、2010、2015百色市遙感影像，在ArcGis10.2及ERDAS IMAGINE_9.2軟件的技術(shù)支持下，對百色市遙感衛(wèi)星影像進(jìn)行處理，從解譯的影像中提取變化信息，最后利用土地利用轉(zhuǎn)移矩陣算出各地類變化。新型城鎮(zhèn)化評價指標(biāo)體系的數(shù)據(jù)主要來源于《百色市統(tǒng)計年鑒》以及政府門戶網(wǎng)站。

1.3 研究方法

1.3.1 土地利用動態(tài)度（LUDI）。用此方法研究百色市三個時間段的土地在數(shù)量、質(zhì)量、類型等方面的面積變化，。其表達(dá)式為：

（1）

式中Ua、Ub分別為研究初、末期某土地利用類型的數(shù)量;T為研究時段。

1.3.2 轉(zhuǎn)移矩陣分析模型。轉(zhuǎn)移矩陣分析模型可全面反映區(qū)域土地利用類型變化的方向、特征與結(jié)構(gòu)。其數(shù)學(xué)形式為

（2）

公式中：n為地類的數(shù)量，S為圖斑面積，i j分別為研究初、期末地類，該模型可以清楚的反映研究時間段內(nèi)初、末期的土地利用類型結(jié)構(gòu)以及各時段內(nèi)各地類的組成、來源。

（3）回歸擬合模型。其公式如下：

Y = aX2 + bX + c （3）

式中：Y為某土地利用類型占總面積的比例; X為新型城鎮(zhèn)化發(fā)展水平。

2 2005 ～2015百色市土地利用時空變化特征分析

2.1 土地利用數(shù)量變化分析

通過分析2005、2010及2015百色市土地利用數(shù)量結(jié)構(gòu)表（表1）和2005、2010及2015百色市土地利用現(xiàn)狀圖（圖1），結(jié)果表明：從時間上看，建設(shè)用地和水域面積不斷增加;耕地和草地面積先減后增;林地面積先增后減;未利用地在三個時期的變化不大。從空間上看，林地、耕地數(shù)量的增減在百色市南部及東部地區(qū)較為明顯;建設(shè)用地的增加則主要是體現(xiàn)在百色市中部及東部地區(qū);草地、水域及未利用地在空間上的變化表現(xiàn)不明顯。

2.2 土地利用變化速度分析

本研究采用土地利用動態(tài)度來分析百色市土地利用變化速度，依據(jù)土地利用動態(tài)度計算公式，結(jié)合不同年份的土地利用類型總面積可得到 3個時期的土地利用動態(tài)度，如表2所示。

由表2可知，2005 ～2015間，耕地面積先減后增，土地利用動態(tài)度為-0.10%;林地面積先增后減，土地利用動態(tài)度為-0.01%;草地面積變化量較小土地利用動態(tài)度約接近于0;建設(shè)用地面積、水域面積持續(xù)增加，土地利用動態(tài)度分別是1.91%和5.60%，其中建設(shè)用地面積的變化程度最大，水域次之;未利用地先增后減，土地利用動態(tài)度為0.23%。

2.3 土地利用類型變化分析

本研究運用ArcGIS10.2對百色市土地利用分類矢量數(shù)據(jù)進(jìn)行空間疊置分析，得出2005～2010百色市土地利用格局轉(zhuǎn)移矩陣（見表3）、2010 ～2015百色市土地利用格局轉(zhuǎn)移矩陣（見表4）。

2005 ～2010百色市各地類的變化，由表3可知：①草地轉(zhuǎn)為建設(shè)用地的量為590.97hm2;②耕地轉(zhuǎn)為建設(shè)用地的量為5032.21hm2;③林地轉(zhuǎn)為建設(shè)用地的量為2828.98;④水域轉(zhuǎn)出為建設(shè)用地的量為518.2hm2。轉(zhuǎn)出水域的量變成林地的面積占最多，草地和耕地次之，未利用地與建設(shè)用地最小;⑤轉(zhuǎn)出的建設(shè)用地主要變成耕地和林地。

2010 ～2015各地類的變化，由表4可知：①耕地轉(zhuǎn)為建設(shè)用地的面積是4982.7hm2;②林地轉(zhuǎn)出為建設(shè)用地的面積是2407.21hm2;③草地轉(zhuǎn)出為建設(shè)用地的面積是601.05hm2;④水域轉(zhuǎn)出為建設(shè)用地的面積是831.89hm2;⑤轉(zhuǎn)入的建設(shè)用地面積小于轉(zhuǎn)出的。

結(jié)合表3、表4可知：在這幾個時期中，建設(shè)用地變化最大，未利用地的變化率次之，說明建設(shè)用地與未利用地的面積還處于波動狀態(tài);從相互轉(zhuǎn)移的數(shù)量來看，建設(shè)用地多是由耕地轉(zhuǎn)入而來的，表明百色市城鎮(zhèn)化發(fā)展建設(shè)多以占用耕地為主。此外，占建設(shè)用地轉(zhuǎn)入面積較大的還有林地。

3土地利用變化與新型城鎮(zhèn)化水平的相關(guān)分析

3.1 百色市新型城鎮(zhèn)化發(fā)展水平測度結(jié)果

基于已有的新型城鎮(zhèn)化評價指標(biāo)體系相關(guān)研究成果，對備選指標(biāo)進(jìn)行了仔細(xì)篩選，指標(biāo)的選取涵蓋生產(chǎn)、生活、生態(tài)等方面的32個具體指標(biāo)體系。通過計算，最后得出百色市新型城鎮(zhèn)化發(fā)展水平。采用改進(jìn)的熵權(quán)法確定權(quán)重，對新型城鎮(zhèn)化的綜合水平進(jìn)行分析，以消除在權(quán)重賦值過程中主觀因素的影響。

通過計算，可得到百色市2005～2015年新型城鎮(zhèn)化發(fā)展水平測度結(jié)果如表5所示。

由表5可知，從2005～2015年，百色市新型城鎮(zhèn)化發(fā)展水平總體上呈現(xiàn)逐年上升趨勢。不斷上升的新型城鎮(zhèn)化發(fā)展水平反映出百色市經(jīng)濟(jì)穩(wěn)步增長的同時，人民無論在生活上還是精神上都得到質(zhì)量上的提高，生產(chǎn)、生態(tài)、生活配置結(jié)構(gòu)日趨合理，使得城市的聚集效益不斷增強，影響力也不在斷擴大。

3.2 土地利用變化與新型城鎮(zhèn)化的比較分析

借鑒已有研究，本文將各地類變化與百色市新型城鎮(zhèn)化進(jìn)行定量比較分析。然后把新型城鎮(zhèn)化發(fā)展水平與各土地利用類型占比值制成散點圖（見圖2、圖3），得出土地利用變化與百色市新型城鎮(zhèn)化發(fā)展水平的定量研究關(guān)系。

由圖2可知，農(nóng)用地占總面積比隨著新型城鎮(zhèn)化水平的上升不斷的下降，兩者在數(shù)量上呈現(xiàn)出顯著的線性關(guān)系，關(guān)聯(lián)系數(shù) R2=0.991，回歸方程為：

Y1=-26.686X2+10.62X+97.943 （11）

式中：Y1為農(nóng)用地占總面積比;X為新型城鎮(zhèn)化發(fā)展水平。從曲線圖可以看出，新型城鎮(zhèn)化水平與農(nóng)用地占總面積比有很密切的相關(guān)性，表明了新型城鎮(zhèn)化的不斷發(fā)展導(dǎo)致城市周邊農(nóng)用地減少。

由圖3可知，建設(shè)用地占總面積比隨著新型城鎮(zhèn)化水平的上升而上升，兩者在數(shù)量上同樣呈現(xiàn)出典型的線性相關(guān)，關(guān)聯(lián)系數(shù) R2=0.9891，回歸方程為：

Y2=20.219X2-8.2605X+1.2603 （12）

式中：Y2為建設(shè)用地占總面積比;X為新型城鎮(zhèn)化發(fā)展水平。從回歸擬合曲線來看，新型城鎮(zhèn)化水平的高低與建設(shè)用地占總面積比之間在數(shù)量上相關(guān)性很好，且新型城鎮(zhèn)化的快速發(fā)展表現(xiàn)為建設(shè)用地規(guī)模的不斷擴大。

3.3 土地利用變化與新型城鎮(zhèn)化的互饋機制分析

本文將土地利用變化與城市化發(fā)展作用機制應(yīng)用于百色市新型城鎮(zhèn)化發(fā)展與土地利用變化的關(guān)聯(lián)分析中，結(jié)合圖4分析，得出新型城鎮(zhèn)化過程土地、產(chǎn)業(yè)、人口城鎮(zhèn)化之間存在正反饋關(guān)系 。

3.3.1 三個時間段百色市的新型城鎮(zhèn)化上水平都呈現(xiàn)上升的趨勢，新型城鎮(zhèn)化水平由2005年的0.1850提高到2015年的0.3138，隨著城鎮(zhèn)化進(jìn)程的推進(jìn)，農(nóng)村人口不斷向城鎮(zhèn)積聚，城鎮(zhèn)人口占比由25%上升至34%。在新型城鎮(zhèn)化初期，農(nóng)村人口向城市遷移，人口的大量轉(zhuǎn)移導(dǎo)致農(nóng)村土地撂荒，城鎮(zhèn)發(fā)展空間不斷擴大以滿足人們的生產(chǎn)生活需要，這需要占用城市周邊的土地以此來擴大城鎮(zhèn)建設(shè)用地規(guī)模，最終導(dǎo)致土地利用結(jié)構(gòu)發(fā)生改變，2005～2015建設(shè)用地比重由0.42%增加至0.65%;農(nóng)業(yè)用地比重由99.00%減少至98.66%。

但是，在新型城鎮(zhèn)化水平不斷提升的過程中，其經(jīng)濟(jì)基礎(chǔ)和社會功能等基本能滿足遷入人口的需求，為了提升并優(yōu)化新型城鎮(zhèn)化發(fā)展水平的目標(biāo)，2010 ～2015百色市城鎮(zhèn)化進(jìn)程的發(fā)展，不再是一味的以建設(shè)用地擴張為主，而是注重以人為核心，注重轉(zhuǎn)變發(fā)展方式，提升集約水平，故2010 ～2015百色市建設(shè)用地增加的幅度在減小，但是新型城鎮(zhèn)化水平仍在穩(wěn)步發(fā)展。

3.3.2 土地利用結(jié)構(gòu)的重組也使得產(chǎn)業(yè)結(jié)構(gòu)重新配置，2005～2015第一產(chǎn)業(yè)所占比重從26.7%變成17.3%，二、三產(chǎn)業(yè)占比上升，由73.3%增至82.7%。城鎮(zhèn)邊緣區(qū)建設(shè)用地擴張多是占用了其邊周邊的農(nóng)用地。

3.3.3 因二、三產(chǎn)業(yè)的收入遠(yuǎn)高于第一產(chǎn)業(yè)，大量農(nóng)村剩余勞動力會向二、三產(chǎn)業(yè)靠攏，城鎮(zhèn)聚集效應(yīng)增強，農(nóng)村人口從事非農(nóng)生產(chǎn)占比將持續(xù)上升，勢必進(jìn)一步加快了新型城鎮(zhèn)化的發(fā)展。但是，這種正向作用機制并非長久存在，當(dāng)新型城鎮(zhèn)化水平達(dá)到其應(yīng)有水平后，正向影響機制就會停止，土地應(yīng)該更多的向空間發(fā)展而不是攤“大餅式”。

4 結(jié)論

以廣西百色市為研究區(qū)，選取2005年、2010年和2015年為研究節(jié)點，在 GIS及 RS技術(shù)的支持下，結(jié)合回歸分析方法，研究百色市新型城鎮(zhèn)化水平與土地利用變化的關(guān)系，結(jié)果表明，2005～2010年土地利用大體表現(xiàn)為建設(shè)用地、林地、水域面積的增加，耕地、草地、未利用地的減少，2010～2015年土地利用大體表現(xiàn)為耕地、草地、建設(shè)用地、水域面積的增加，林地、未利用地的減少。

新型城鎮(zhèn)化水平的高低與土地利用類型比存在著明顯的相關(guān)性，低水平的新型城鎮(zhèn)化導(dǎo)致了農(nóng)用地的萎縮和建設(shè)用地的擴大;當(dāng)新型城鎮(zhèn)化水平達(dá)到一定高度時，粗濫的發(fā)展方式轉(zhuǎn)變成更注重發(fā)揮資源的綜合效益的發(fā)展模式，農(nóng)業(yè)人口合理向非農(nóng)產(chǎn)業(yè)轉(zhuǎn)移，建設(shè)用地多向空間發(fā)展，土地利用類型的變化比一般城鎮(zhèn)化進(jìn)程的結(jié)構(gòu)優(yōu)化，這一變化會在城鎮(zhèn)周邊地區(qū)表現(xiàn)得非常明顯。

新型城鎮(zhèn)化水平較低時，僅表現(xiàn)為建設(shè)用地的擴張;隨著各方面的不斷完善，土地、產(chǎn)業(yè)及人口城鎮(zhèn)化之間存在著正向顯著相關(guān)，當(dāng)發(fā)展理念轉(zhuǎn)變、水平不斷提高時，土地城鎮(zhèn)化將帶動產(chǎn)業(yè)及人口城鎮(zhèn)化，產(chǎn)業(yè)和人口城鎮(zhèn)化也將推動土地城鎮(zhèn)化。

[參考文獻(xiàn)]

[1] 余霜，李光，楊筱珊，等.新型城鎮(zhèn)化建設(shè)對安順市土地利用結(jié)構(gòu)的影響研究[J].南方農(nóng)業(yè)，2016，10（22）：40-42.

[2] 張超，李丁，魏秀梅，等.西北河谷型城市新型城鎮(zhèn)化與土地利用效益耦合協(xié)調(diào)發(fā)展研究——以蘭州市為例[J].蘭州大學(xué)學(xué)報（自然科學(xué)版），2015，51（02）：173-179.

[3] 彭沖，陳樂一，韓峰.新型城鎮(zhèn)化與土地集約利用的時空演變及關(guān)系[J].地理研究，2014，33（11）：2005-2020.

[4] 羅琛，胡寶清，韋俊敏，等.縣域土地變化效應(yīng)與新型城鎮(zhèn)化耦合分析——以廣西上林縣為例[J].廣西師范學(xué)院學(xué)報（自然科學(xué)版），2015，32（02）：78-85.

[5] Kenneth F.Wieand.An extension of the monocentric urban spatial equilibrium model to a multicenter setting：The case of the two-center city. Journal of Urban Economics，1987.

[6] Downs A.What does smart growth really mean？. Planning . 2001.

[7] 陳永林，謝炳庚，李曉青，等.2003—2013年長沙市土地利用變化與城市化的關(guān)系[J].經(jīng)濟(jì)地理，2015，35（01）：149-154.

[8] 邸利，張仁陟，張富，等.基于RS與GIS的定西市安定區(qū)土地利用變化與土壤侵蝕研究[J].干旱區(qū)資源與環(huán)境，2011，25（02）：40-45.

[9] 趙永平.新型城鎮(zhèn)化發(fā)展水平測度及其時空差異分析[J].西安電子科技大學(xué)學(xué)報（社會科學(xué)版），2016，26（5）：60-68.

[10] 魯春陽，楊慶媛，文楓.城市化與城市土地利用結(jié)構(gòu)關(guān)系的協(xié)整檢驗與因果分析——以重慶市為例[J].地理科學(xué)，2010，30（04）：551-557.

[11] 李永樂，吳群，舒幫榮.城市化與城市土地利用結(jié)構(gòu)的相關(guān)研究[J].中國人口·資源與環(huán)境，2013，23（04）：104-110.

[12] 王福海，陳丹.基于RS和GIS的土地利用變化時空特征研究——以重慶主城區(qū)為例[J].重慶理工大學(xué)學(xué)報（自然科學(xué)），2016，30（07）：78-85.

[13] 張會英，劉萍萍，李旭祥.陜北紅堿淖流域土地利用動態(tài)變化研究[J].環(huán)境科學(xué)與技術(shù)，2011，34（03）：180-183

[14] 王福海.特大山地城市建設(shè)用地擴展及其遙感地表特征參量響應(yīng)研究[D]. 重慶：重慶工商大學(xué)，2015.

[15] 趙永平.新型城鎮(zhèn)化發(fā)展水平測度及其時空差異分析[J].西安電子科技大學(xué)學(xué)報（社會科學(xué)版），2016，26（05）：60-68.

[收稿日期]2018-12-03

[作者簡介]張宗刻（1990—），男，廣西百色人，助理工程師，從事國有土地管理工作。