基于PLUS的佛山市高明區土地利用模擬及變化特征分析

夏義雄，康停軍，周劍鋒

(佛山市測繪地理信息研究院，廣東 佛山 528000)

1 引 言

土地利用/土地覆蓋變化(land use/land cover change，LUCC)是評估區域環境變化和可持續發展的重要組成部分[1]，模擬和預測土地利用變化機制和規律對區域合理開發利用土地資源，提高土地開發效能，促進區域可持續發展具有重要意義[2]。

目前模擬土地利用變化模型主要有系統動力學模型(SD)、CLUE-S、人工神經網絡、元胞自動機(CA)等模型。CA模型對復雜地理過程具有較好的模擬和預測能力[3～5]，近年來，許多學者通過引入多元影響因素對模型進行優化，提出了Logistics-CA[6]，FLUS[7]，PLUS[8]等多種模型，其中PLUS模型以斑塊作為演變單元分析多種土地利用類型變化的潛在驅動因素，能更準確、真實地模擬土地利用變化的情況[8]。

佛山市高明區作為粵港澳大灣區向西發展的樞紐和自然資源保護重點區域。伴隨珠三角樞紐(廣州新)機場預計近期(2035年)建設，高明區將迎來巨大的發展機遇，如何協調好城鎮的快速發展與自然資源保護的關系，成為高明區高質量發展的重要命題。本研究擬利用PLUS模型對高明區2030年土地利用情況進行模擬并分析土地利用變化特征，為城市資源利用管理、城市可持續化建設提供科學理論依據。

2 研究區概況

選取廣東省佛山市高明區作為研究區，高明區地處佛山市西部，位于東經112°22′34″～112°55′06″，北緯22°38′46″～23°01′05″。高明區占地面積 937.81 km2，森林覆蓋率高達51.17%，內部有云勇國家森林公園、凌云山森林公園、桫欏自然保護區等多處城市重要綠地林地。高明區不僅在佛山市生態環境保護中起重要作用，還是未來珠三角核心區向外輻射的重要樞紐，珠三角新干線機場、廣湛高鐵、珠江肇高鐵、肇順南城際、佛肇高速等一系列重大交通項目的建設也為高明區建設區域性綜合樞紐和提升城市能級提供了新的可能。

3 數據和方法

3.1 數據

收集了2010年、2015年和2020年三期高明區土地利用數據，包含建筑用地、耕地、林地、草地、水體等6種土地利用類型，數據來源為CLCD30M土地利用分類數據(CLCD數據下載地址：https：//zenodo.org/record/5210928)。為有效反映研究區的土地利用發展規律、社會經濟情況和地理條件，收集了地形、人口、GDP、交通、氣候、降水等12類土地利用變化驅動因素數據?？紤]到自然環境保護和生態紅線等政策影響，還收集了高明區“三區三線”數據，用以控制城市擴張邊界。

3.2 研究框架與方法

結合高明區歷史土地利用數據、土地利用變化驅動因素和Markov模型預測的土地利用目標，基于PLUS模型對高明區土地利用情況進行模擬。通過KAPPA系數和OA指標驗證PLUS模型精度，并對2030年土地利用情況進行預測，分析了2020年～2030年高明區土地利用變化特征。研究流程如圖1所示。

圖1 研究流程圖

(1)PLUS土地利用模擬模型

PLUS模型是一種基于柵格的斑塊生成土地利用模擬模型[8]。PLUS模型通過獲取土地利用類型變化與驅動因素之間的關系，挖掘土地利用類型轉化規則，使用帶有隨機種子的CA模型控制新的土地利用斑塊生成。PLUS模型相較于其他土地利用模擬模型使用了全新的土地利用轉化規則挖掘策略，重點關注土地利用轉化結果忽略其來源，降低了分析土地利用變化規則的復雜度，對真實土地利用情況具有更好的模擬和仿真能力。

(2)Markov模型

Markov模型已普遍運用于土地利用演變和模擬研究中[5，6]。根據Markov模型隨機過程原理，在土地利用的演變過程中，某時期的土地利用狀態只與其前一時期的土地利用類別和土地利用類型擴展概率有關。因此，某一時刻的土地利用狀態可以用如下方法進行預測：

LULCt+1=LULCt×Pij

(1)

其中，LULCt+1表示t+1時刻的土地利用狀態，Pij表示每類土地利用類型的擴展概率。

(2)

上式中：0≤Pij≤1；i，j為土地利用類型。

4 結果

4.1 土地利用模擬結果

利用高明區2010年和2015年兩期的數據對不同限制條件和驅動因素的2020年土地利用情況進行了模擬和精度驗證，所得結果KAPPA系數為0.85，OA為0.91，各項精度指標均表明PLUS模型能有效提取和模擬高明區土地利用變化規律。

利用PLUS模型對2030年高明區土地利用情況進行預測，結果如圖2所示。從整體上看，2030年高明區土地利用占比呈現林地>耕地>建筑用地>水體>草地>裸地，與2020年區域格局基本保持一致。

圖2 高明區2020年、2030年土地利用情況

通過預測結果與2020年的各類型土地的數據對比如表1所示，按照現有自然資源條件和社會發展政策保持不變的情況下，高明區未來土地利用變化呈現“兩減兩增”的特點。2020年-2030年增長最多的為耕地，比2020年增長了 7.48 km2，減少最多的為林地，比2020年減少了 9.19 km2，建筑用地增長了 4.88 km2，水體減少了 3.09 km2，草地和裸地的面積基本未發生變化。

表1 高明區2020年與2030年各類型土地利用面積對比(km2)

4.2 土地利用變化特征分析

通過土地利用類型轉移矩陣如表2所示，可以定量分析各土地利用類型的變化。從表中數據可以看出高明區土地利用類型轉變具有林地轉耕地，耕地轉建筑用地的特征，發生用地類型變化的林地中有94.74%轉變為耕地，發生用地類型變化的耕地中有83.78%的轉變為建筑用地，上述兩種類型用地的轉變在所有用地類型轉變中占比超過70%。這種土地利用變化特征說明高明區仍處于城市發展高速建設期，需要更多的建筑用地滿足城市發展所帶來的基礎設施、住房等需求，耕地面積的持續增加則表明高明區持續實施用地整理和耕地保護政策取得了較好成效。

表2 高明區2020年、2030年各類型土地利用面積對比(km2)

通過土地利用變化空間分布圖(圖3)可知，耕地轉變的區域主要出現在高明區東部的高明河與高明大道周邊，也是高明區城市擴張的主要區域。林地轉耕地主要出現在高明區中部的苗村、新村和橫坑村等較為偏遠鄉村，緊鄰云勇國家森林公園、鹿洞山森林公園、皂幕山森林公園。土地利用類型轉變特征與空間分布表明高明區以耕地補改措施作為橋梁，增加建筑用地指標以解決建設用地指標不足對城市發展的制約，形成了城市擴張區周邊以耕地轉變為建筑用地，偏遠鄉村周邊以林地轉變為耕地的用地變化特征。

圖3 高明區土地利用變化空間分布圖

5 結 論

通過收集佛山市高明區2010年、2015年和2020年三期土地利用數據，基于PLUS模型模擬高明區2030年土地利用情況，對高明區的土地利用變化趨勢特征進行了分析。研究發現2030年高明區的城鎮發展仍然維持在較高水平且城鎮發展與自然保護的矛盾有加劇的風險。其次，受到建設用地指標的限制，高明區采取包括土地整治、耕地補改等措施，補充建設用地指標，通過對歷史土地利用變化規律的總結和分析，發現耕地向建筑用地轉變和林地向耕地轉變是未來高明區土地利用轉變的主要方向。最后，通過分析土地利用變化的空間分布，發現城市周邊耕地的保護和偏遠鄉村周邊林地的開發是未來高明區土地利用監管需要關注的重點，云勇國家森林公園、鹿洞山森林公園、皂幕山森林公園及其周邊林地有較大可能發生用地類型轉變，應納入重點巡查范圍，減少林地退化速率，合理統籌規劃建設用地，促進區域可持續發展。