滇南開遠(yuǎn)市中心城鎮(zhèn)土地利用演變及其影響因素

李琛，趙筱青

(云南大學(xué) 地球科學(xué)學(xué)院，昆明 650500)

0 引言

城鎮(zhèn)是經(jīng)濟(jì)發(fā)展的重要支撐點(diǎn)，隨著“一帶一路”戰(zhàn)略的不斷推進(jìn)，沿線國(guó)家的城鎮(zhèn)化發(fā)展加快。特別在面向東南亞開發(fā)中的地區(qū)，由于其區(qū)位優(yōu)勢(shì)和“一帶一路”的戰(zhàn)略支撐，逐漸改變中小城鎮(zhèn)發(fā)展相對(duì)滯后的情況[1]。然而急速城鎮(zhèn)化的典型特征就在于劇烈的土地利用變化，出現(xiàn)了區(qū)域生態(tài)環(huán)境嚴(yán)重破壞[2]、城鎮(zhèn)化與耕地保護(hù)矛盾突出等問題[3]，土地利用演變特征是人類塑造活動(dòng)空間的綜合反映，為人類合理、可持續(xù)地使用土地提供依據(jù)。因此探討土地利用時(shí)空演變特征及影響因素分析具有重要意義。國(guó)內(nèi)外學(xué)者在土地利用研究方面進(jìn)行了諸多的探討，主要集中于土地變化[4]和優(yōu)化方法[5-6]，探討影響因素[7-8]，并建立起了一系列土地利用變化的研究方法和模型，例如土地利用動(dòng)態(tài)度、信息熵、轉(zhuǎn)移矩陣及空間動(dòng)態(tài)模型、主成分分析、回歸分析、地理探測(cè)器法及Logistic回歸模型等方法[9-13]。然而在空間尺度方面，研究主要集中于東部經(jīng)濟(jì)發(fā)達(dá)地區(qū)以及生態(tài)脆弱等大尺度區(qū)域，分析其土地利用演變特征，而對(duì)于西部中小城鎮(zhèn)尺度的研究較少[14]，特別是面向東南亞國(guó)家“一帶一路”沿線的中小城鎮(zhèn)研究。

云南省紅河州的個(gè)舊、開遠(yuǎn)、蒙自為核心組成的滇南中心城市處于與越南“兩廊一圈”中西廊-昆河經(jīng)濟(jì)走廊，不僅是面向東南亞的門戶城市，也是國(guó)土空間開發(fā)利用與保護(hù)的重點(diǎn)區(qū)域[15]。其中，開遠(yuǎn)市地處中國(guó)-東盟兩大市場(chǎng)結(jié)合部和中國(guó)昆明、越南河內(nèi)兩大城市輻射交匯點(diǎn)，在“一帶一路”戰(zhàn)略背景下，作為滇南地區(qū)的商貿(mào)物流中心和面向東盟開放的交通樞紐的開遠(yuǎn)市得到了迅速的發(fā)展，其中心城鎮(zhèn)土地利用變化大、速度快。然而，快速城鎮(zhèn)化中建設(shè)用地的迅速擴(kuò)張給區(qū)域土地資源承載力和生態(tài)環(huán)境帶來了巨大壓力。為促進(jìn)區(qū)域可持續(xù)發(fā)展，緩解用地矛盾，以開遠(yuǎn)市中心城鎮(zhèn)為研究區(qū)，選取2009年、2014年和2018年三期土地利用數(shù)據(jù)，在系統(tǒng)分析土地利用時(shí)空演變特征的基礎(chǔ)上，運(yùn)用二元Logistic回歸模型探究土地利用演變的影響因素,以期為開遠(yuǎn)市中心城鎮(zhèn)土地資源的合理利用、深入推進(jìn)滇南中心城市建設(shè)及“一帶一路”沿線中小城鎮(zhèn)的國(guó)土空間規(guī)劃提供參考。

1 研究區(qū)概況

開遠(yuǎn)市位于云南省東南部(103°04′～103°43′E，23°30′～23°58′N)，隸屬于云南省紅河州，地處中國(guó)-東盟兩大市場(chǎng)結(jié)合部和中國(guó)昆明、越南河內(nèi)兩大城市輻射交匯點(diǎn)，是滇南中心城市群的重要組成部分，下轄靈泉街道、樂白道街道、小龍?zhí)舵?zhèn)等7個(gè)鄉(xiāng)鎮(zhèn)。中心城鎮(zhèn)位于中部，由樂白道街道和靈泉街道兩個(gè)鄉(xiāng)鎮(zhèn)組成。其中樂白道街道是開遠(yuǎn)市的政府駐地，與靈泉街道一起構(gòu)成了開遠(yuǎn)市的政治、經(jīng)濟(jì)、文化和交通中心(圖1)。中心城鎮(zhèn)的人口約占開遠(yuǎn)市總?cè)丝诘?5%、財(cái)政收入占全市的73.5%[16]，在開遠(yuǎn)市處于重要地位，因此中心城鎮(zhèn)城市用地、建制鎮(zhèn)、村莊以及交通運(yùn)輸用地面積擴(kuò)大成為了必然趨勢(shì)。

圖 1 開遠(yuǎn)市中心城鎮(zhèn)區(qū)位示意Figure 1 Location of the Kaiyuan Central Town

2 數(shù)據(jù)來源及處理

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

土地利用數(shù)據(jù)采用2009年開遠(yuǎn)市第二次全國(guó)土地調(diào)查數(shù)據(jù)、2014年和2018年土地變更調(diào)查成果；社會(huì)經(jīng)濟(jì)數(shù)據(jù)主要來源于《2009年開遠(yuǎn)市統(tǒng)計(jì)年鑒》《開遠(yuǎn)經(jīng)濟(jì)工作手冊(cè)(2016年)》《開遠(yuǎn)經(jīng)濟(jì)工作手冊(cè)(2018年)》；部分自然和經(jīng)濟(jì)數(shù)據(jù)主要參考《開遠(yuǎn)市樂百街道辦事處土地利用總體規(guī)劃(2010—2020)》《開遠(yuǎn)市靈泉道街道辦事處土地利用總體規(guī)劃(2010—2020)》。

2.2 數(shù)據(jù)處理

開遠(yuǎn)市土地利用分類以第二次全國(guó)土地調(diào)查土地利用現(xiàn)狀分類標(biāo)準(zhǔn)(GB/T 21010-2007)為基礎(chǔ)，結(jié)合實(shí)際情況，為進(jìn)一步體現(xiàn)開遠(yuǎn)市中心城鎮(zhèn)土地利用演變情況，對(duì)土地利用類型進(jìn)行調(diào)整，分為耕地、園地、林地、草地、交通運(yùn)輸用地、水域及水利設(shè)施用地、設(shè)施農(nóng)用地、裸地、城市、建制鎮(zhèn)、村莊和工礦用地等12個(gè)土地利用類型。

2.3 研究方法

2.3.1 土地利用時(shí)間演變特征分析方法

引入土地利用面積變化量和單一土地利用類型動(dòng)態(tài)度[17]，反映開遠(yuǎn)市中心城鎮(zhèn)土地利用時(shí)間變化，公式如下：

S=Ub-Ua

(1)

(2)

式(1～2)中：S表示土地利用面積變化量，K表示單一土地利用類型動(dòng)態(tài)度，Ua、Ub分別為研究初末期土地利用類型的面積；T為研究時(shí)段。

2.3.2 土地利用空間演變特征分析方法

(1)土地利用程度空間分布

借鑒劉紀(jì)遠(yuǎn)的土地利用程度綜合指數(shù)模型和土地利用程度變化模型[18-20]，參考相關(guān)研究[21-22]，確定開遠(yuǎn)市中心城鎮(zhèn)土地利用程度分級(jí)指數(shù)，具體公式如下：

(3)

(4)

式(3～4)中:Lj為j時(shí)間研究區(qū)土地利用程度綜合指數(shù);△Lb-a為土地利用程度變化量；Ai為研究區(qū)內(nèi)第i級(jí)土地利用程度分級(jí)指數(shù)。裸地為1級(jí)；林地、草地、水域?yàn)?級(jí)；耕地、園地為3級(jí)；城鎮(zhèn)村及工礦用地、交通運(yùn)輸用地為4級(jí)。

(2)土地利用空間動(dòng)態(tài)變化模型土地利用空間動(dòng)態(tài)模型通過土地轉(zhuǎn)移矩陣反映土地利用變化的類型和轉(zhuǎn)移狀態(tài)，為土地利用空間類型變化過程定量分析提供依據(jù)[23-24]。土地利用動(dòng)態(tài)變化空間分析模型為[12]：

(5)

(6)

CCLi=TRLi+IRLi

(7)

式(5～7)中:IRLi表示第i種土地利用類型轉(zhuǎn)入速率；TRLi表示第i種土地利用類型轉(zhuǎn)出速率；CCLi表示第i種土地利用類型的空間變化速率；LA(i,i1)-ULA,表示第i種土地利用類型轉(zhuǎn)移部分的面積；LA(i,i2)-ULA,表示新增部分面積；t1、t2分別表示研究期初、期末的時(shí)間。

2.3.3 土地利用影響因素分析方法

(1)二元Logistic回歸模型

采用二元Logistic回歸模型方法對(duì)影響因素進(jìn)行分析。當(dāng)假設(shè)某一柵格出現(xiàn)某種土地利用的概率為P時(shí)，則二元Logistic回歸模型可以得到表達(dá)式[25]：

(8)

式(8)中：Pi是每個(gè)柵格可以出現(xiàn)第i類型土地的概率，X是影響因子，βn是Xn的系數(shù)，βn值越大，則表示第n種影響因子X與土地利用類型的相關(guān)程度越高。

通過函數(shù)變化，Logistic回歸模型表達(dá)式可以表示為：

(9)

式(9)中：該比值稱為OR值為發(fā)生比，即發(fā)生與不發(fā)生之間的比例，EXP(β)>1發(fā)生概率增加，EXP(β)<1發(fā)生概率減少，EXP(β)=1發(fā)生概率不變。

(2)檢驗(yàn)方法

利用ROC曲線下面積值來驗(yàn)證二元Logistic回歸模型的準(zhǔn)確性[26]。當(dāng)ROC值越接近1，則效果越好，ROC>0.7時(shí)，說明影響因子對(duì)土地利用類型變化具有較強(qiáng)的解釋性[27]。

(3)因子選擇及處理

土地利用變化不僅深受自然因素的影響，也受社會(huì)經(jīng)濟(jì)因素的影響。依據(jù)參考文獻(xiàn)[28-30]，結(jié)合開遠(yuǎn)市中心城鎮(zhèn)實(shí)際情況，其地勢(shì)起伏較大，山區(qū)地帶地形交錯(cuò)復(fù)雜，地形地貌條件對(duì)研究區(qū)有重要影響；中心城鎮(zhèn)在整個(gè)開遠(yuǎn)市中人口社會(huì)經(jīng)濟(jì)發(fā)展相對(duì)較快，直接影響研究區(qū)的土地利用變化；作為面向東盟開放的交通樞紐，區(qū)域可達(dá)性對(duì)土地利用變化有著明顯的影響作用?？紤]所能獲取的相關(guān)數(shù)據(jù)，從自然因素、社會(huì)經(jīng)濟(jì)因素和可達(dá)性3個(gè)方面選擇土地利用時(shí)空演變影響因子(表1)。

表 1 土地利用時(shí)空演變影響因素 Table 1 Impact factors of spatial and temporal evolution of the land use

將數(shù)據(jù)進(jìn)行柵格化，并把所有柵格數(shù)據(jù)統(tǒng)一為投影坐標(biāo)體系(CGCS2000)和空間分辨率(30 m×30 m)，將各個(gè)因子?xùn)鸥駭?shù)據(jù)轉(zhuǎn)化為TIF格式，由于城市用地和建制鎮(zhèn)用地在建設(shè)過程中的主要影響因素相似且建制鎮(zhèn)用地在中心城鎮(zhèn)中所占面積較小，因此因子分析時(shí)將城市用地和建制鎮(zhèn)合并為城鎮(zhèn)用地進(jìn)行分析。

將柵格數(shù)據(jù)轉(zhuǎn)為ASCII數(shù)據(jù)，再利用CLUE-S模型中的convert工具將ASCII數(shù)據(jù)轉(zhuǎn)化為SPSS中可識(shí)別的數(shù)據(jù)并對(duì)其進(jìn)行Z標(biāo)準(zhǔn)化。

3 結(jié)果分析

3.1 土地利用時(shí)空演變特征

3.1.1 土地利用時(shí)間演變特征

(1)土地利用面積變化

利用公式(1～2)得到土地利用面積變化量和單一土地利用類型動(dòng)態(tài)度結(jié)果(表2)。開遠(yuǎn)市中心城鎮(zhèn)以林地、耕地和草地為主，2009—2018年林地是最主要的土地利用類型，均占總面積的40%以上；其次是耕地，均占總面積的20%以上；城鎮(zhèn)村及工礦用地僅占整個(gè)中心城鎮(zhèn)的5%～6%左右。耕地、園地、林地所占百分比分別下降0.53%、0.38%、0.20%，而城鎮(zhèn)村及工礦用地占比則上升0.63%。

表 2 土地利用面積變化量及單一土地利用動(dòng)態(tài)度Table 2 Area change of land use and dynamics of single land use

2009—2018年土地?cái)?shù)量變化主要以耕地、園地、林地、城市用地和交通運(yùn)輸用地為主，耕地、園地和林地減少，而城市用地和交通運(yùn)輸用地增加。其中，在2009—2014年中心城鎮(zhèn)土地利用數(shù)量變化最多的是城市用地，增加了283.31 hm2，占總變化量的23.81%，其次為耕地，耕地減少了250.33 hm2，占總變化量的21.04%。2014—2018年中心市鎮(zhèn)增加的土地利用數(shù)量變化最多的是耕地，減少了75.44 hm2，占總變化量的26.67%，其次是水域及水利設(shè)施用地，增加了53.48 hm2，占總變化量的18.90%。

(2)土地利用動(dòng)態(tài)變化

2009—2018年土地利用動(dòng)態(tài)度在10%以上的是交通運(yùn)輸用地、設(shè)施農(nóng)用地以及村莊，分別增速57.36%、35.14%、14.09%，土地利用動(dòng)態(tài)度在1%以上的是城市和水域及水利設(shè)施用地，分別增速2.18%和1.18%，其他土地利用類型的土地利用動(dòng)態(tài)度則不明顯。其中2009—2014年交通運(yùn)輸用地和設(shè)施農(nóng)用地增速高達(dá)92.25%和45.83%，建制鎮(zhèn)以及城市增速則為21.1%和3.69%； 2014—2018年設(shè)施農(nóng)業(yè)地增速為6.62%、建制鎮(zhèn)增速為2.59%，交通運(yùn)輸用地增速則為2.45%，而城市增速僅為0.25%。

總之，前5年為開遠(yuǎn)市中心城鎮(zhèn)土地利用動(dòng)態(tài)度最為活躍的時(shí)期，城鎮(zhèn)村、交通和設(shè)施農(nóng)用地均迅速增長(zhǎng)，其中建制鎮(zhèn)土地利用動(dòng)態(tài)度明顯高于城市用地，交通運(yùn)輸用地增速最為顯著；后5年土地利用動(dòng)態(tài)度明顯緩慢，僅有水利及水域設(shè)施用地和村莊的動(dòng)態(tài)度超過了上一階段，其中水利及水域設(shè)施則是由于溝渠建設(shè)而引起的。

3.1.2 土地利用空間演變特征

(1)土地利用程度空間變化

通過構(gòu)建100 m×100 m、200 m×200 m、500 m×500 m、1 000 m×1 000 m尺度計(jì)算土地利用程度綜合指數(shù)，認(rèn)為在200 m×200 m尺度上選擇采樣點(diǎn)能較好地反映中心城鎮(zhèn)城市用地、建制鎮(zhèn)用地等高利用程度地區(qū)形態(tài)和擴(kuò)展方向等空間變化情況，根據(jù)公式(3～4)，通過插值得到土地利用程度綜合指數(shù)和土地利用程度變化量的空間分布(圖2～3)。

2009—2018年土地利用程度較高的區(qū)域主要集中在中部城區(qū)，并且面積在增加，主要是由于沿交通運(yùn)輸用地的增多，城市用地、建制鎮(zhèn)和村莊的擴(kuò)展，導(dǎo)致土地利用程度提高。

當(dāng)ΔL>0時(shí)，說明土地利用程度得到提高；當(dāng)ΔL≤0時(shí)，說明土地利用程度不變或者下降。由圖3可知，2009—2014年土地利用程度變化呈現(xiàn)出隨城鎮(zhèn)用地內(nèi)部填充和沿交通路線擴(kuò)散而提高的趨勢(shì)，利用程度提高區(qū)域主要集中在中部和東南部，沿交通路線和圍繞城鎮(zhèn)用地呈團(tuán)塊狀分布。2014—2018年土地利用程度有升有降，但總體是在提高，主要集中在南部地區(qū)，城鎮(zhèn)用地?cái)U(kuò)展形態(tài)逐漸由面狀擴(kuò)展轉(zhuǎn)向點(diǎn)狀擴(kuò)展，分散在各個(gè)建制鎮(zhèn)和村莊。

圖 2 土地利用程度綜合指數(shù)空間變化Figure 2 Spatial change of land use degree composite index

圖 3 綜合土地利用程度空間變化示意Figure 3 Spatial change of comprehensive land use degree

(2)土地利用空間動(dòng)態(tài)變化分析

依據(jù)土地利用空間動(dòng)態(tài)變化模型并結(jié)合土地利用轉(zhuǎn)移矩陣，利用公式(5～7)得到土地利用類型轉(zhuǎn)入速率(IRL)、轉(zhuǎn)出速率(TRL)和空間變化速率(CCL)(表3)。

表 3 各年土地利用動(dòng)態(tài)變化空間分析/%Table 3 Area change of land use and dynamic of single land use

2009—2018年開遠(yuǎn)市中心城鎮(zhèn)的CCL變化較為明顯，其中CCL交通運(yùn)輸用地為57.37%、CCL設(shè)施農(nóng)用地為35.47%和CCL建制鎮(zhèn)為15.45%，大量耕地、園地以及林地轉(zhuǎn)變?yōu)槌鞘杏玫亍⒔ㄖ奇?zhèn)、交通運(yùn)輸用地、水域及水利設(shè)施用地，表明城鎮(zhèn)擴(kuò)張和耕地短缺之間的矛盾及城鎮(zhèn)發(fā)展與生態(tài)環(huán)境保護(hù)的矛盾越來越激烈。

2009—2014年開遠(yuǎn)市中心城鎮(zhèn)的CCL變化顯著。變化明顯的土地利用類型為交通運(yùn)輸用地、設(shè)施農(nóng)用地和建制鎮(zhèn)，其CCL分別為92.25%、45.87%和22.32%，以耕地、園地向城市用地、交通運(yùn)輸用地和建制鎮(zhèn)轉(zhuǎn)移為主。這一時(shí)期，城鎮(zhèn)的快速發(fā)展侵占了大量的耕地，城鎮(zhèn)擴(kuò)張和耕地短缺之間的矛盾愈發(fā)突出。

2014—2018年開遠(yuǎn)市中心城鎮(zhèn)的CCL變化相對(duì)放緩，變化明顯的土地利用類型為：CCL設(shè)施農(nóng)用地為6.88%、CCL水域及設(shè)施為3.07%、CCL建制鎮(zhèn)為2.60%，呈現(xiàn)出多種土地利用類型轉(zhuǎn)變，土地利用結(jié)構(gòu)正在逐漸優(yōu)化。與上一階段相比，IRL與TRL之間的差值進(jìn)一步縮小，區(qū)域多種土地利用類型轉(zhuǎn)入轉(zhuǎn)出；耕地、園地保持一定的動(dòng)態(tài)平衡，耕地、園地的IRL迅速提高，耕地和園地被侵占的趨勢(shì)有所減緩，但I(xiàn)RL

表 4 ROC值檢驗(yàn)結(jié)果 Table 4 ROC value test result

表 5 耕地時(shí)空演變影響因子二元Logistic回歸結(jié)果 Table 5 The logistic regression results of the impact factors of the spatial-temporal evolution of cultivated land

3.2 土地利用變化影響因素分析

3.2.1 模型檢驗(yàn)

采用Pearson相關(guān)系數(shù)進(jìn)行多重共線性分析，將高度相關(guān)的因子剔除(|r|≥0.8)消除多重共線性[31]。按照剔除后的影響因子對(duì)不同地類進(jìn)行二元Logistic回歸，利用ROC對(duì)結(jié)果進(jìn)行檢驗(yàn)(表4)。

開遠(yuǎn)市中心城鎮(zhèn)各土地利用類型除了林地和草地的ROC值小于0.7以外，其他土地利用類型的ROC值均大于0.7(表4)，說明所構(gòu)建的模型對(duì)于耕地、林地和草地的擬合效果一般，但整體來說模型較為可靠。

3.2.2 土地利用時(shí)空演變影響因子分析

重點(diǎn)對(duì)開遠(yuǎn)市中心城鎮(zhèn)土地利用變化明顯的耕地和城鎮(zhèn)用地進(jìn)行分析。

(1)耕地時(shí)空演變影響因子分析

由表5可知，2009—2018年耕地的影響因子主要為距農(nóng)村道路距離(X9)、坡度(X2)以及海拔(X1),其中距農(nóng)村道路距離(X9)是最主要的影響因子，2009年系數(shù)β為-0.621，2014年為-0.635、2018年為-0.631，發(fā)生比EXP(β)2009為 0.573 7、EXP(β)2014為 0.530、EXP(β)2018為0.532，距農(nóng)村道路距離與耕地呈負(fù)相關(guān)，說明耕地多分布于距農(nóng)村道路較近的區(qū)域，農(nóng)村道路的建設(shè)，有利于農(nóng)業(yè)生產(chǎn)資料的運(yùn)輸和農(nóng)業(yè)的快速發(fā)展。距農(nóng)村道路距離發(fā)生變化對(duì)耕地的時(shí)空演變產(chǎn)生的影響最大。次要影響因子是坡度，耕地與坡度有明顯的負(fù)相關(guān)關(guān)系，耕地多分布于坡度較為平緩的地區(qū)，平緩的土地有利于耕作業(yè)的進(jìn)行。海拔與耕地為正相關(guān)，表明耕地多分布于海拔相對(duì)較高的區(qū)域。 除此之外，農(nóng)民人均所得(X7)、距公路道路距離(X8)以及距鄉(xiāng)村點(diǎn)距離(X11)也對(duì)耕地產(chǎn)生影響，只是相對(duì)于主要影響因子而言影響較弱。

2009—2018年距農(nóng)村道路距離、距公路道路距離的系數(shù)β絕對(duì)值總體上在增加，但坡度、海拔系數(shù)及人口密度的系數(shù)β絕對(duì)值在減少，表明可達(dá)性對(duì)耕地農(nóng)業(yè)生產(chǎn)的影響因素在逐漸增大，耕地雖然在靠近鄉(xiāng)村點(diǎn)的區(qū)域，但多分布在人口較少的地區(qū)；區(qū)域坡度和海拔等自然因素雖然對(duì)耕地有著限制作用，但這種作用在減弱。

(2)城鎮(zhèn)用地時(shí)空演變影響因子分析

2009—2018年城鎮(zhèn)用地影響因子主要為海拔(X1)、距公路道路距離(X8)、人口密度(X4)和人均財(cái)政收入(X6)(表6)。其中，海拔在3年中均是最主要的影響因子，自然因素對(duì)城鎮(zhèn)用地的發(fā)展起著基礎(chǔ)作用，城鎮(zhèn)用地多分布于海拔較低、坡度較緩、人口密度較高和公路密集的區(qū)域?？傮w來說，城鎮(zhèn)用地與海拔的系數(shù)β絕對(duì)值有所下降，而與人口密度、距公路道路距離的系數(shù)β絕對(duì)值有所上升，城鎮(zhèn)用地?cái)U(kuò)展過程中，自然因素的影響作用在不斷減弱，社會(huì)經(jīng)濟(jì)因素與可達(dá)性因素的影響作用在不斷增強(qiáng)。但自然因素依舊是最主要的影響因子，影響建設(shè)的適宜性程度，是城鎮(zhèn)用地布局的基礎(chǔ)，而人口密度越高，更需要城鎮(zhèn)用地的擴(kuò)展，交通運(yùn)輸條件有利于促進(jìn)城鎮(zhèn)用地的發(fā)展。

4 結(jié)論與討論

4.1 結(jié)論

(1)中心城鎮(zhèn)時(shí)間演變特征

土地利用類型以林地和耕地為主。2009—2018年城市、交通運(yùn)輸、建制鎮(zhèn)用地面積分別增長(zhǎng)301.65、234.99、53.05 hm2，增長(zhǎng)比例各為0.49%、0.38%、0.09%；而耕地、園地和林地分別減少325.78、231.96、123.73 hm2，減少比例各為0.53%、0.38%、0.20%；耕地、園地、林地面積減少帶來的區(qū)域可持續(xù)問題不容忽視。

土地利用時(shí)間演變有明顯的階段性，2009—2014年為迅速擴(kuò)張階段，2014—2018年為緩慢發(fā)展和內(nèi)部調(diào)整階段。2009—2014年變化最為明顯的是城市用地。2014—2018年變化最為明顯的是耕地。

(2)中心城鎮(zhèn)空間演變特征

2009—2014年土地利用空間變化數(shù)量大、速度快，土地利用程度提高的區(qū)域集中在中部，沿交通路線和圍繞城鎮(zhèn)用地周圍呈塊狀和帶狀分布，耕地、園地面積的89.7%向城市用地、交通運(yùn)輸用地和建制鎮(zhèn)轉(zhuǎn)移。

2014—2018年土地利用空間變化速度減緩，土地利用程度變化由面狀擴(kuò)展轉(zhuǎn)向點(diǎn)狀擴(kuò)展。多種土地利用類型發(fā)生轉(zhuǎn)入和轉(zhuǎn)出。城鎮(zhèn)擴(kuò)張和耕地短缺之間的矛盾越來越突出，減少的耕地面積中有78.8%轉(zhuǎn)移為城鎮(zhèn)用地及交通運(yùn)輸用地，解決城鎮(zhèn)與生態(tài)環(huán)境協(xié)調(diào)發(fā)展的問題迫在眉睫。

(3)中心城鎮(zhèn)影響因素分析

土地利用時(shí)空演變的主要影響因素為海拔、坡度、人口密度、距農(nóng)村道路距離、距公路道路距離等。隨著區(qū)域經(jīng)濟(jì)發(fā)展，自然因素的影響在不斷減弱，社會(huì)經(jīng)濟(jì)以及交通運(yùn)輸條件的影響作用在不斷增強(qiáng)。

4.2 討論

在土地利用變化方面，綜合多種研究方法可以較好地反映土地利用變化特征。本研究從地理學(xué)的時(shí)空觀出發(fā)，采用土地利用的數(shù)量、動(dòng)態(tài)度等方法綜合分析其時(shí)間演變特征，借助土地利用程度空間分布、空間動(dòng)態(tài)變化綜合分析其空間演變特征，結(jié)合時(shí)間和空間視角，在一定程度上完善了對(duì)該區(qū)域的研究。

在土地利用影響因素方面，周銳等研究表明二元Logistic模型不僅可以在短時(shí)間序列中，反映土地利用變化在時(shí)間和空間上與各個(gè)因子之間的關(guān)系，并在中小城鎮(zhèn)尺度中效果較好[14]。本研究從自然因素、社會(huì)經(jīng)濟(jì)因素、可達(dá)性因素3個(gè)方面選擇影響因子，運(yùn)用二元Logistic模型揭示開遠(yuǎn)市中心城鎮(zhèn)土地利用變化的影響因素。然而，土地利用變化必然是多個(gè)因素共同影響的，對(duì)于文化、政策等難以量化的因素仍難以定量探討。針對(duì)這一情況，王興友將基本農(nóng)田和自然保護(hù)區(qū)作為政策限制因子進(jìn)行討論[32]，然而不能將其作為快速發(fā)展的中心城鎮(zhèn)土地變化的影響因素。因此，中心城鎮(zhèn)土地利用變化雖然受到文化環(huán)境和相關(guān)政策的影響，但本研究并未考慮。今后應(yīng)加強(qiáng)影響因子的可量化和空間化，使選取的因子更加廣泛，結(jié)論更加準(zhǔn)確。