陳 青，謝 釗，唐南奇

（福建農(nóng)林大學(xué)資源與環(huán)境學(xué)院，福建 福州350002）

土地利用和覆被變化（Land Use and Land Cover Change，LUCC）作為“國際地圈生物圈計劃”（IGBP）和“國際全球環(huán)境變化人文因素計劃”（IHDP）等研究的核心領(lǐng)域之一，近年來受到越來越多專家學(xué)者的關(guān)注[1]。土地利用和覆被變化是一個動態(tài)演變過程，建立模型是對其進行動態(tài)研究的有效手段，亦是深入了解其變化機制和動力的重要途徑[2]。常用于LUCC 研究的模型有：灰色模型、回歸分析模型、Markov 模型、系統(tǒng)動力學(xué)模型、人工神經(jīng)網(wǎng)絡(luò)模型和元胞自動機模型等[3-4]。Markov 模型作為一種經(jīng)典的預(yù)測模型，在LUCC 研究中應(yīng)用最廣泛。目前，關(guān)于土地利用和覆被變化的研究主要集中在城市及區(qū)域尺度上，如戴靚等[5]通過解譯江蘇省2000年和2008年的遙感影像，得到南京市轄區(qū)在這兩個時期的土地利用數(shù)據(jù)，然后采用Markov 模型模擬預(yù)測南京市2016年的土地利用趨勢。近年來也有一些專家將Markov 模型應(yīng)用于縣域尺度的土地利用變化研究，如朱萌等[6]在對武功縣土地利用和覆被動態(tài)研究中，通過Markov 模型預(yù)測了武功縣未來12 a 土地利用和覆被格局的動態(tài)演化趨勢，且預(yù)測結(jié)果通過了模型效應(yīng)系數(shù)的檢驗，說明將Markov 模型用于縣域尺度的土地利用變化模擬是可行的。此外，還有部分學(xué)者將Markov 模型引入到其他模型之中，利用Markov 模型的優(yōu)勢彌補其他模型的不足，以提高相關(guān)模型模擬的可靠性及精確度。如吳季秋等[7]將Markov 模型與CA 元胞自動機模型耦合成CA-Markov 模型，利用Markov 長期預(yù)測的優(yōu)勢以及元胞自動機模型模擬復(fù)雜空間演變的優(yōu)勢，從土地利用適宜性原則與生態(tài)保護的角度對海南省八門灣的土地利用變化進行動態(tài)預(yù)測，結(jié)果表明人類活動尤其是相關(guān)土地政策在很大程度上影響著八門灣的土地利用變化趨勢。

福清市地處海峽西岸經(jīng)濟區(qū)中部沿海，位于閩江口金三角經(jīng)濟圈南翼，東臨臺灣海峽，是環(huán)南、北太平洋兩大經(jīng)濟繁榮圈的交接點。近年來隨著海峽西岸經(jīng)濟區(qū)的全面建設(shè)，福清市城市建設(shè)和經(jīng)濟發(fā)展飛速，人地矛盾不斷加劇，土地利用和覆被變化顯著。研究利用2009年和2012年福清市土地利用現(xiàn)狀矢量化數(shù)據(jù)，在ArcGIS、Matlab 等軟件的支持下，建立Markov 模型對該地區(qū)的土地利用和覆被變化進行模擬預(yù)測，分析福清市土地利用和覆被變化的演化特點，以期對福清市合理利用和保護土地資源具有現(xiàn)實意義和參考價值。

1 Markov 模型預(yù)測的基本原理

Markov 模型預(yù)測是一類特殊的隨機過程，它是由俄國數(shù)學(xué)家馬爾科夫（A.A Markov）于1907年提出來的，該隨機過程最大特點是具有無后效性[8]。無后效性是指某動態(tài)系統(tǒng)在T＋1 時刻的狀態(tài)只與T 時刻的狀態(tài)有關(guān)，而與T 時刻之前的狀態(tài)無關(guān)，不受T 時刻之前的狀態(tài)所影響[9]。將Markov 模型應(yīng)用于預(yù)測土地利用和覆被類型的動態(tài)變化趨勢是非常合適的，因為土地利用和覆被類型變化是受人為和自然因素共同影響的，而這些影響因素又具有隨機性、不確定性，使得土地利用和覆被類型的變化也成為一種隨機運動過程。Markov 模型預(yù)測土地利用和覆被變化的基本原理，具體而言是在利用不同歷史時期土地利用類型轉(zhuǎn)移矩陣生成區(qū)域土地利用和覆被類型變化轉(zhuǎn)移概率矩陣的基礎(chǔ)上，預(yù)測某一特定情景下未來土地利用變化的趨勢[10-11]。

研究表明，轉(zhuǎn)移概率矩陣的確定是成功應(yīng)用Markov 模型的關(guān)鍵[12]。土地利用和覆被類型轉(zhuǎn)移概率是指在土地利用變化的過程中，某一土地利用類型在下一時刻轉(zhuǎn)化為其他類型的可能性。由各土地利用狀態(tài)的轉(zhuǎn)移概率組成的矩陣P，即為轉(zhuǎn)移概率矩陣。該矩陣可以定量地說明不同土地利用類型間的相互轉(zhuǎn)化情況及轉(zhuǎn)移速率，其數(shù)學(xué)表達式為：

式中：n 表示研究區(qū)土地利用和覆被類型的數(shù)目；Pij表示土地利用/覆被類型由i 轉(zhuǎn)化為j 的概率。該矩陣的每一個元素需要同時滿足以下兩個條件：

如果已知土地利用和覆被類型系統(tǒng)中n=0 時的初始狀態(tài)矩陣那么經(jīng)過n 次轉(zhuǎn)移后的狀態(tài)矩陣為：

式中：P（n）表示初始狀態(tài)土地利用/覆被類型概率矩陣經(jīng)n 步轉(zhuǎn)移后的概率矩陣，經(jīng)過n 步轉(zhuǎn)移后的系統(tǒng)狀態(tài)概率矩陣為A（n）=A（n-1）P=A（n-2）PP=…=A（0）Pn，所以P（n）= Pn。由此可見，由初始狀態(tài)矩陣和轉(zhuǎn)移概率矩陣可以確定任何時刻的系統(tǒng)狀態(tài)概率矩陣。

2 建模過程及趨勢預(yù)測

2.1 基礎(chǔ)數(shù)據(jù)的獲取及處理

研究所用基礎(chǔ)數(shù)據(jù)來源于福清市2009年和2012年土地利用變更調(diào)查數(shù)據(jù)。參照我國《土地利用分類》國家標(biāo)準(zhǔn)（CB/T 21010-2007），并結(jié)合土地修編現(xiàn)狀分類標(biāo)準(zhǔn)以及福清市土地資源的具體情況和研究目的，將福清市土地利用和覆被類型劃分為9 大類，分別為：耕地、園地、林地、草地、建設(shè)用地（包括商服用地、工礦倉儲用地、住宅用地、公共管理與公共服務(wù)用地及水利設(shè)施用地）、交通用地、水域（水庫水面及河流水面）、其他農(nóng)用地（農(nóng)村道路、坑塘水面及溝渠）和特殊及未利用地（未利用地包括內(nèi)陸灘涂、沿海灘涂及裸地）。借助ArcGIS軟件，對福清市2009年和2012年的土地利用變更調(diào)查數(shù)據(jù)進行矢量化處理并進行數(shù)據(jù)融合，獲取福清市2009年和2012年各土地利用和覆被類型的面積，見表1。

表1 福清市2009/2012 年各土地利用和覆被類型面積（hm2）

表2 福清市2009年各土地利用和覆被類型向2012年的轉(zhuǎn)化情況（hm2）

在Arcgis 軟件中對兩個時期的數(shù)據(jù)進行空間疊置分析，得到轉(zhuǎn)移的各用地類型及屬性數(shù)據(jù)，獲取2009年各土地利用和覆被類型向2012年的轉(zhuǎn)化情況，見表2。

2.2 轉(zhuǎn)移概率矩陣的建立

2.2.1 初始狀態(tài)矩陣的確定 將土地利用和覆被類型劃分為一系列相互轉(zhuǎn)化的狀態(tài)，各狀態(tài)在系統(tǒng)中所占的比例即每種土地利用類型面積占全部土地利用類型面積的百分比作為各狀態(tài)的初始概率，從而形成初始狀態(tài)矩陣A（0）。統(tǒng)計福清市2009年各土地利用類型的面積，進一步計算得到各覆被類型的面積百分比，從而形成初始狀態(tài)矩陣A（0）（表3）。

表3 初始狀態(tài)矩陣A（0）

2.2.2 轉(zhuǎn)移概率矩陣的確定 根據(jù)福清市2009年到2012年各用地類型的面積轉(zhuǎn)移情況，以年為基本時間單位，求出時間段內(nèi)某類土地利用和覆被類型的年平均轉(zhuǎn)移概率。如研究區(qū)2009年草地中有部分到2012年變?yōu)楦?、建設(shè)用地、交通用地、其他農(nóng)用地、特殊及未利用地和園地等，后者的面積占2009年草地總面積的比例再除以年數(shù)，即為其年平均轉(zhuǎn)移概率。將草地轉(zhuǎn)化為其他土地利用和覆被類型的概率作為第一行，耕地轉(zhuǎn)化為其他土地利用和覆被類型的概率作為第二行，依次類推，便構(gòu)成了初始狀態(tài)轉(zhuǎn)移概率矩陣P，見表4。

2.3 動態(tài)模擬及檢驗

借助MATLAB R2013a 平臺，以1 a 為一個步長，根據(jù)公式①、②建立Markov 預(yù)測模型。由此可模擬預(yù)測出初始年份（2009年）后各時期的土地利用和覆被類型的轉(zhuǎn)移概率和對應(yīng)的轉(zhuǎn)移概率矩陣。如通過計算n=3 的轉(zhuǎn)移概率矩陣，可以模擬得到福清市2012年各土地利用/覆被類型的轉(zhuǎn)移概率矩陣，見表5。

由表5 的模擬結(jié)果，結(jié)合表3 中的數(shù)據(jù)可獲得福清市2012年各土地利用和覆被類型面積的預(yù)測值，將計算得到的預(yù)測值和2012年的實際值進行對比（表6），在SPSS18.0 軟件中對預(yù)測值和實際值進行配對樣本T 檢驗，設(shè)置置信水平為95%，結(jié)果顯示：t=0.386，df=8，雙側(cè)P=0.709＞0.05，即預(yù)測值與實際值之間并無顯著性差異，表明通過Markov 模型模擬和預(yù)測福清市土地利用和覆被類型的動態(tài)演變趨勢是可行的。

表5 2012年福清市土地利用和覆被類型的轉(zhuǎn)移概率矩陣

表6 2012年福清市土地利用和覆被類型面積的實際值與預(yù)測值 （hm2）

2.4 變化趨勢預(yù)測

對模型檢驗后，在保持當(dāng)前干擾因素不變的前提下，在MATLAB R2013a 平臺上編程依次輸入n=6、11、16、21，運用Markov 模型進行預(yù)測，計算得到2015年、2020年、2025年和2030年福清市各土地利用和覆被類型的面積比例，見表7。

表7 2015～2030年福清市各土地利用和覆被類型的面積比例預(yù)測結(jié)果 （%）

3 結(jié)果與分析

由表7 可知，福清市未來18 a 內(nèi)，耕地、林地、其他農(nóng)用地、特殊及未利用地、園地等呈逐年減少趨勢，而建設(shè)用地、交通用地則呈逐年上升趨勢，這與該市前期土地利用/覆被類型變化趨勢相似。在減少的用地類型中，耕地減少幅度較大，由2015年的21.30%下降到2030年的20.46%，減少的面積主要轉(zhuǎn)化為了建設(shè)用地和交通用地；特殊及未利用地減少幅度次之，由2015年的16.95%下降至2030年的16.38%，主要是其中的未利用地轉(zhuǎn)化為了建設(shè)用地和交通用地；林地由2015年的30.67%降至2030年的30.45%，減少的部分主要轉(zhuǎn)化為了耕地和建設(shè)用地；其他農(nóng)用地由2015年的8.76%下降至2030年的8.48%，主要轉(zhuǎn)化為了耕地、建設(shè)用地和交通用地等；園地由2015年的5.73%下降至2030年的5.59%，減少面積主要轉(zhuǎn)化為了草地和建設(shè)用地。在增加的土地利用類型中增加幅度最大的是建設(shè)用地，由2015年的11.38%上升至2030年的13.17%，主要由耕地與特殊及未利用地轉(zhuǎn)化而來；其次是交通用地，雖然在各覆被類型中所占的面積較少，但一直處于增長狀態(tài)，由2015年的1.56%上升至2030年的1.84%；此外，水域和草地所占的比例雖有變化但變幅較小可以忽略，可認(rèn)為其沒有變動，表明福清市的河流水域及草地得到了有效的保護。綜上分析可知，按照當(dāng)前的發(fā)展趨勢，在未來的18 a 中，福清市城市建設(shè)用地將會不可避免的占用大量農(nóng)林用地及未利用地。福清市建設(shè)用地和交通用地的持續(xù)增加說明該市的經(jīng)濟發(fā)展態(tài)勢良好，但耕地與特殊及未利用地的減少也說明這種良好態(tài)勢是以糧食安全和生態(tài)保護為代價的。

4 結(jié)論與建議

研究以2009年和2012年福清市不同土地利用類型之間的相互轉(zhuǎn)化數(shù)據(jù)為基礎(chǔ)，利用Matlab 編程Markov 模型對福清市2012年后的土地利用和覆被類型的變化趨勢進行預(yù)測。模型預(yù)測結(jié)果與實際情況基本吻合，這說明將Markov 模型用于縣域尺度土地利用和覆被格局變化趨勢的預(yù)測是可行的。

根據(jù)Markov 模型預(yù)測福清市土地利用和覆被類型動態(tài)變化趨勢，結(jié)果顯示該市耕地、林地、其他農(nóng)用地、特殊及未利用地比例呈逐年減少趨勢，而建設(shè)用地和交通用地則呈逐年增加趨勢。這表明隨著福清市城市的建設(shè)和經(jīng)濟的發(fā)展，建設(shè)用地對土地的需求與農(nóng)林用地保護之間的矛盾將日漸突出。該結(jié)果應(yīng)引起相關(guān)政府部門的高度重視，并采取切實有效的措施，如耕地占卜平衡、城鄉(xiāng)建設(shè)用地增減掛鉤、節(jié)約集約利用建設(shè)用地及嚴(yán)厲打擊違法違規(guī)用地等，既保證經(jīng)濟發(fā)展又保障生態(tài)糧食安全。

Markov 預(yù)測模型建立的前提是假定外界干擾因素不變，由已知時期的土地利用變化情況推導(dǎo)出其未來的發(fā)展演化趨勢，但實際情況中人為干擾仍然是影響福清市土地利用和覆被類型變化的主要因素，所以該預(yù)測模型存在一定的局限性。在今后的研究中應(yīng)盡力把對土地利用影響較大的相關(guān)土地政策考慮進去，使研究結(jié)果更加準(zhǔn)確可靠；此外，在建立Markov 預(yù)測模型時，時間尺度是建立該模型的一個重要環(huán)節(jié)，不同的時間尺度可能會導(dǎo)致不同的模擬結(jié)果。研究僅采取2009年和2012年兩個時期的數(shù)據(jù)，時間尺度有點小，在今后的研究中應(yīng)盡量選取較長及較多的時間尺度來進預(yù)測研究。

[1]楊維鴿.基于CA-Markov模型和多層次模型的土地利用變化模擬和影響因素研究——以陜西省米脂縣高西溝村為例[D].西安：西北大學(xué)碩士學(xué)位論文，2011.

[2]韓春建，吳克寧，劉德元，等.基于馬爾科夫模型的鄭州市郊區(qū)多方案耕地保有量預(yù)測[J].土壤，2011，43（3）：453-458.

[3]夏小梅，陳柯林.基于遙感與GIS的LUCC及馬爾科夫過程趨勢預(yù)測——以重慶市南岸區(qū)為例[J].水土保持研究，2009，16（5）：19-24.

[4]劉慧璋，郭青霞，王曰鑫，等.基于Markov的山西岔口小流域土地利用變化預(yù)測[J].山西農(nóng)業(yè)大學(xué)學(xué)報，2012，32（1）：53-57.

[5]戴 靚，吳紹華.基于Markov模型的南京土地利用時空變化研究[J].資源開發(fā)與市場，2011，27（5）：417-420.

[6]朱 萌，馬孝義，劉雪嬌.基于馬爾科夫模型的武功縣土地利用/覆被動態(tài)變化研究[J].水土保持研究，2013，20（5）：64-68.

[7]吳季秋，俞花美，肖 明，等.基于CA-Markov模型的海灣土地利用動態(tài)監(jiān)測[J].環(huán)境工程技術(shù)學(xué)報，2012，（6）：531-539.

[8]Hulst R.On the dynamics of vefetation：Markov chains asmodels of successi on[J].Vegetation，1979，40（1）：3-14.

[9]GengX M，LiL.Markov process functional-sinfinan ceand in surance[J].Applied Mathematics，2009，24（1）：21-26.

[10]喻 峰.基于Markov-CA的土地利用變化預(yù)測研究[J].國土資源情報，2009，（4）：38-46.

[11]羅 平，姜仁榮，李紅旮，等.基于空間Logistic和Markov模型集成的區(qū)域土地利用演化方法研究[J].中國土地科學(xué)，2010，（1）：31-36.

[12]郜魯豪，徐 旌.基于Markov模型的安寧市土地利用預(yù)測[J].云南地理環(huán)境研究，2010，22（2）：87-89.