土地利用/覆被變化分析方法研究

丑慶川，胡 濤 ，何詩雨 ，史秀華

(1.深圳大學生命科學學院，廣東深圳 518060； 2.深圳大學高爾夫學院，廣東深圳 518060)

?

土地利用/覆被變化分析方法研究

丑慶川1，胡 濤1，何詩雨1，史秀華2*

(1.深圳大學生命科學學院，廣東深圳 518060； 2.深圳大學高爾夫學院，廣東深圳 518060)

土地利用/覆被變化一直以來都被認為是全球變化的重要組成部分和主要驅動因子之一，近30年來，越來越受到各國相關學者的普遍關注，目前已成為國際學術界研究的焦點和重點。國內外學者基于不同的目的，從不同的角度出發構建了大量模型，以便深入了解土地利用/覆被變化的驅動因素和機制，預測其未來發展變化趨勢。該研究從土地利用速度、土地利用程度、土地利用變化趨勢和土地利用的景觀指數等方面研究土地利用動態變化的模型，其中重點介紹當前應用最廣泛的模型，同時為土地利用/覆被變化模型的未來發展指明方向。

模型；動態變化；土地利用/覆被變化；預測

土地利用/覆被變化(Land Use and Cover Change，簡稱LUCC)包含兩層不同卻又密切相關的含義，即土地利用和土地覆被[1-2]。20世紀90年代以來，土地利用/覆被變化被認為是全球變化的重要組成部分和主要驅動因子之一[3]，越來越受到各國相關學者的普遍關注，目前已成為國際學術界研究的焦點和重點[4]。國際上自1993年國際科學聯合會與國際社會科學聯合會聯合成立了土地利用/覆被變化核心項目計劃委員會后，一些積極參與全球環境變化的國際組織和國家紛紛跟進，啟動了各自的土地利用/覆被變化研究項目[5]。

近年來國內學者也在土地利用/覆被變化方面做了大量研究[6-17]，主要集中在土地利用/覆被變化狀況的評估、土地利用/覆被變化的建模與預測、土地利用/覆被變化的集成技術研究、土地利用/覆被變化的驅動力研究、土地利用/覆被變化的景觀生態學理念等方面，從不同的角度出發構建模型，深入了解土地利用/覆被變化的成因和過程，預測未來發展變化趨勢。筆者從土地的“過去-現在-未來”出發，主要從土地利用速度、土地利用程度、土地利用變化趨勢和土地利用的景觀指數等方面重點介紹研究土地利用動態變化的模型。

1 土地利用速度

1.1 單一類型土地利用速度單一類型土地利用速度又稱土地利用類型的動態度，是指某一時間范圍之內研究區域內一種土地利用類型的數量大小變化情況[3,18-24]，其表達式為：

式中,R指某一時間范圍內某一土地利用類型的利用速度；Ua、Ub分別指研究初期、研究末期某一土地利用類型的土地利用數量；T是研究的時間跨度，當T以年記時，R值就是某種土地利用類型的年變化率。

1.2 綜合土地利用速度綜合土地利用速度又稱綜合土地利用動態度，是用來反映區域內土地利用變化的綜合速度[3,19,21-24]，其表達式為：

式中,Lc指研究時間段內研究區域的綜合土地利用速度；Ai指研究初始第i類土地利用類型的面積大??；ΔAi-j指研究時間段內第i類土地利用類型面積數量轉為非i類土地利用類型面積數量的絕對值；T指研究的時間跨度(a)。

2 土地利用程度

土地利用程度是反映人類對土地進行利用的深度與廣度的概念，它體現的是土地利用中土地的自然屬性、人為因素和環境因素的綜合效應[23]。借鑒劉紀遠[25]等建立的土地利用程度數量化綜合指標體系，并賦予不同土地利用類型土地利用程度指數。

2.1 土地利用程度綜合指數土地利用程度綜合指數反映的是土地利用/覆被變化的總體變化趨勢[3,23-26]，其表達式為：

Lj=∑(Ai×Ci)

式中,Lj指研究區域內的土地利用程度綜合指數；Ai指研究區域內第i類土地利用類型的土地利用程度指數；Ci指某研究時間跨度內研究區域內第i類土地利用類型面積占總研究區域面積的百分比。

2.2 土地利用程度變化土地利用程度變化是指在特定研究區域內，不同土地利用類型在某一時間跨度內共同變化作用的結果，它可對研究區域內土地利用的動態變化及變化方向從變化量和變化率兩個層次進行定量揭示[3,24]。土地利用程度變化量D和變化率λ表達式分別為：

式中,D為研究的時間跨度內段研究區域的土地利用程度變化量；λ指研究的時間跨度內研究區域的土地利用程度變化率；Ai指研究區域內第i類土地利用類型的土地利用程度指數；Cit1和Cit2分別指研究初始期和研究末期研究區域內第i類土地利用類型的面積百分比；t1和t2分別指研究始初期和研究末期。

2.3 土地利用類型變化劇烈程度土地利用類型的變化劇烈程度是指某一時間跨度內研究區域內不同土地利用類型之間相互轉化的速度，其可定量揭示研究區域內不同土地利用類型的動態變化及變化方向[19]，其表達式為：

式中,S為與時間T對應的究區域土地利用類型變化劇烈程度，為方便將其擴大10倍；Si為研究初始期第i類土地利用類型的面積數量；ΔSi-j為研究時間跨度內第i類土地利用類型面積數量轉化為非i類土地利用類型面積的絕對值；T為研究時間跨度的長度，通常為年。

3 土地利用類型變化趨勢模型

模型是研究土地利用/覆被變化的重要工具[27]，土地利用類型變化趨勢模型可對未來土地利用/覆被狀況進行模擬和預測[28]。目前，預測土地利用類型變化趨勢模型有：馬爾科夫模型[19,29-30]、元胞自動機模型[31-40]、系統動力學模型[29,41-44]、CLUE模型[45-53]、多元統計模型[54-59]、回歸分析模型[29-51]、類似杜能模型[60]、智能體模型[61-69]等多種模型，但其中馬爾科夫模型和元胞自動機模型的接受度最高且應用最為廣泛，故該研究重點介紹這2種模型的基本原理。

3.1 馬爾科夫模型馬爾科夫模型是利用某一運動系統的現在狀況及其發展動向預測該系統未來狀況的分析方法[29]。在馬爾科夫模型中較常用的是一階馬爾科夫過程，即一個運動系統從某一特定的t時刻狀態運動到t+1時刻狀態的概率只和t時刻狀態有關，而與其他時刻的狀態無關[19]。馬爾科夫模型的基本方程為：

馬爾科夫模型成功應用的關鍵在于對其轉移概率的確定，一階馬爾科夫模型主要應用于較小空間尺度的土地利用/覆被變化中[29]。土地利用/覆被變化中建立馬爾科夫模型的前提是土地利用/覆被轉移概率固定不變，而受到社會經濟驅動的土地利用/覆被轉移概率在長時間尺度內是不可能不變的，因此該模型只適合短期的預測[30]。

3.2 元胞自動機模型元胞自動機是一種簡單的通過局部轉換規則來模擬復雜空間結構的模型[31]，一個標準的模型由元胞、空間、鄰居和轉換規則組成，空間被分割成規則的網格，每一個網格即為一個元胞，每個元胞具有一個狀態，每個元胞各自的狀態隨時間而變化，每個元胞都是根據局部規則來更新狀態[32]。即每一個元胞在某一特定時刻的狀態只決定于上一時刻該元胞的狀態和與該元胞相鄰的所有元胞上一時刻的狀態[33-34]。常規的標準模型可以簡單地表達如下：

St+1=f(St,N)

式中,S是狀態；N是鄰居；f是轉換函數；t是時間。

在土地利用/覆被變化的模擬中，元胞一般代表地塊信息，每一個元胞的狀態由模型的轉換規則決定[35]。轉換規則被認為是元胞自動機的核心，目前已有包括多準則判斷[36]、邏輯回歸[37]、人工神經網絡模型[38]、蟻群算法[39]、人工免疫系統[40]等多種方法來確定CA模型中轉換規則的參數。

4 土地利用類型的景觀指數和數量結構多樣性模型

4.1 景觀偏離度景觀偏離度是指人為建造的人文景觀對自然景觀的偏離程度。這一指標主要從數量上揭示人類活動對景觀改變的程度，數值越大說明人文景觀占的比重越大，表明人類對自然的改造越多[20]。其表達式為：

式中，M為景觀偏離度；Sa、Sb、Sc、Sd、Se和Sf分別為研究區域內各種建設用地、人工水域、耕地、園地、人工林地和人工草地的面積；S總為研究區域的總面積。

4.2 景觀多樣性景觀多樣性是指不同的景觀類型在空間結構，功能機制和時間動態方面的多樣化和變異性，它反映的事景觀的復雜性。引入景觀多樣性的公式來計算土地利用類型的多樣性，其指數越大，說明多樣性越多[20]。表達式為：

式中，H(V)為多樣性指數;pi為i景觀所占的百分比;n為變量數。

4.3 土地數量結構多樣性指數模型區域土地利用/覆被類型數量結構多樣性指數可以用于定量評價區域內各類土地利用/覆被類型的齊全程度和多樣化狀況[16]，進而分析人類活動對土地的影響程度[23]。故采用Gibbs-Mirtin多樣化指數來度量，其模型為：

式中，GM指多樣化指數;Si指第i種土地利用類型的面積。

5 展望

近20年來，隨著土地利用/覆被變化研究的不斷深入，國內外學者在土地利用/覆被變化分析方法所代表內涵上的看法基本趨于一致；土地利用程度反映的是人類對土地進行利用的深度與廣度，可定量揭示區域內土地利用的動態變化及變化方向[23]；預測模型是對未來土地利用/覆被狀況進行模擬和預測的重要工具等[28]。

隨著GIS和遙感技術的不斷發展，不斷涌現出許多新的進行土地利用/覆被變化研究的方法，社會因素和經濟因素對土地利用/覆被變化的不斷影響在新方法的探索中越來越受到重視。傳統分析方法將不能滿足對新技術條件下獲得的數據的分析要求，因此，急需更多的理論研究來不斷提出新的分析方法以滿足大數據條件下的土地利用/覆被變化研究。

[1] ALLEN T R,WANG Y,CRAWFORD T W.Remote Sensing and GIS science in Geomorphology:Remote sensing of land cover dynamics.Academic Press,San Diego,CA,80-102.

[2] RADOSLAVA K,MIRIAM K,JOZEF N,et al.Land-use and land-cover changes in rural areas during different political systems:A case study of Slovakia from 1782 to 2006[J].Land Use Policy,2014,36:554-566.

[3] 李德一，張安定，王大鵬，等.基于遙感的10年間龍口市土地利用變化分析[J].水土保持研究，2007,14(1)：221-225.

[4] 陳懷亮，徐祥德，劉玉潔.土地利用與土地覆蓋變化的遙感監測及環境影響研究綜述[J].氣象科技，2005,33(4)：289-294.

[5] 劉新衛，陳百明，史學正.國內LUCC研究進展綜述[J].土壤，2004，36(2):132-135.

[6] 王秀蘭，包玉海．土地利用動態變化研究方法探討[J]．地理科學進展，1999，18(1)：81-87．

[7] 陳百明，劉新衛，楊紅．LUCC研究的最新進展評述[J]．地理科學進展，2003，22(1)：22-29．

[8] 張秋菊，傅伯杰，陳利頂，等．黃土丘陵溝壑區縣域耕地變化驅動要素研究：以安塞縣為例[J]．水土保持學報，2003，17(4)：146-152.

[9] 梁偉，楊勤科．基于RS的黃河中游多沙粗沙區土地利用變化分析[J]．水土保持研究，2006，13(5)：90-92．

[10] 常慶瑞，劉京，楊勤科，等．基于RS與GIS的黃土丘陵溝壑區土地覆蓋動態研究[J]．應用生態學報，2003，14(12)：2187-2190．

[11] 溫仲明，焦峰，張曉萍，等．黃土丘陵區紙坊溝流域60年來土地利用格局變化研究口]．水土保持學報，2004，18(5)：125-133.

[12] 徐勇，郭騰云．黃土丘陵區中部土地利用變化歷程及驅動力分析[J]．水土保持研究，2006，13(2)：1-4．

[13] 馬建華．現代自然地理學[M]．北京：北京師范大學出版社，2002：495-506．

[14] 蔡運龍．土地利用/土地覆被變化研究：尋求新的綜合途徑[J]．地理研究，2001，20(6)：645-652．

[15] 朱會義，李秀彬．環渤海地區土地利用的時空分析[J]．地理學報，2001，56(3)：253-259．

[16] 朱風武，彭補拙．溫州市土地利用空間格局研究[J]．經濟地理，2001，21(1)：101-104．

[17] 劉曉曼，莊大方，屈冉，等.基于遙感數據的1989-2009年拉魯濕地變化研究[J]．濕地科學，2013(11)：1-5.

[18] 張澤晗，張韜，白恒情，等.3S技術在呼和浩特市城市綠地動態變化中的應用研究[J]．內蒙古農業大學學報，2004，25(1)：31-35.

[19] 于兵.基于GIS的大慶市土地利用動態變化分析與預測[J]．東北林業大學學報，2008，36(2)：45-48.

[20] 崔振洋，李曉亮.基于RS與GIS的山西省土地利用動態變化分析[J].中國農學通報，2006,22(6)：418-420.

[21] 陳本清，徐涵秋.廈門市土地利用年際變化遙感分析[J]．地球信息科學，2004，6(3)：99-104.

[22] 趙東波，梁偉，楊勤科，等.陜北黃土丘陵區近30年來土地利用動態變化分析[J]．水土保持通報，2008，28(2)：22-26.

[23] 王安周，張貴賓，鄭杰，等.新鄉市土地利用動態變化分析[J].水土保持研究，2007，14(4)：355-357.

[24] CHOU Q C,LI R L,XU H L,et al.Study on the Land Use/Cover Changes of Mangrove Reserve of Shenzhen Bay based on Remote Sensing Date[J].Advanced Materials Research，2014,962/965:2533-2540.

[25] 劉紀遠．中國資源環境遙感宏觀調查與動態研究[M]．北京：中國科學技術出版社，1996：262-275．

[26] 莊大方，劉紀遠.中國土地利用程度的區域分異模型研究[J].自然資源學報，1997,12(2)：105-111.

[27] MOTTET A，LADET S，COQUE N，et al.Agriculture land-use change and its drivers in mountain landscapes:A case study in the Pyrenees[J].Agriculture,Ecosystems and Environment,2006,114:296-310.

[28] MATTHEWS R B,GILBERT N G,ROACH A,et al.Agent-based land-use models:a review of application[J].Landscape Ecology,2007,22(10):1447-1459.

[29] 黃秋昊，蔡運龍.國內幾種土地利用變化述評[J]．中國土地科學，2005，19(5)：25-30.

[30] 張靖，同麗嘎，張慶，等.土地利用/覆蓋變化研究方法[J]．內蒙古大學學報：自然科學版，2013，44(3)：329-336．

[31] 周成虎，孫戰利，謝一春.地理元胞自動機研究[M]．北京：科學出版社，1999.

[32] TAKEYAMA M,COUCLELI S H.Map dynamics:Integrating cellular automata and GIS through Ceo-algebra[J]．International Journal of Geographical Information Science,1997，11:73-91.

[33] 郭紅.元胞自動機模型在森林景觀生態規劃中的應用[J]．東北農業大學學報，2007,38(6)：814-817.

[34] CLARKE K C,HOPPEN S,GAYDOS L.A self-modifying cellular automaton model of historical urbanization in the San Francisco Bay area[J]．Environment and Planning B,1997,24:247-261.

[35] 王羊，高陽，趙琳，等.元胞自動機模型的尺度敏感性分析[J]．北京大學學報：自然科學版，2011,47(4)：750-758.

[36] WU F,WEBSTER C J.Simulation of land development through the integration of cellular automata and multi-criteria evaluation[J].Environment and Planning B,1998,25:103-126.

[37] WU F.Calibration of stochastic cellular automata:the application to rural-urban land conversions[J].International Journal of Geographical Information Science,2002,16(8):795-818.

[38] LI X,YEH A G O.Neural-network-based cellular automata for simulating multiple land use changes using GIS[J].International Journal of Geographical Information Science,2002,16(4):323-343.

[39] 劉小平,黎夏,葉嘉安，等.利用蟻群智能挖掘地理元胞自動機的轉換規則[J].中國科學:D 輯,2007,37(6):824-834.

[40] 劉小平,黎夏,張嘯虎，等.人工免疫系統與嵌入規劃目標的城市模擬及應用[J]．地理學報,2008,63(8):882-894.

[41] 擺萬奇.深圳市土地利用動態趨勢分析[J]．自然資源學報，2000,15(2)：112-116.

[42] 徐瑞祥.區域耕地總量動態平衡模型研究——以溫州市為例[J]．經濟地理，2002,22(4)：435-439.

[43] 何春陽，史培軍，李景剛，等.中國北方未來土地利用變化情景模擬[J]．地理學報，2004,59(4)：599-607.

[44] STEPHENNE N,LAMBIN E F.A dynamic simulation model of land-use changes in Sudano-sahelian countries of Africa (SALU)[J]．Environment,2001,85:145-161.

[45] VERBURG P H,OVERMARS K P.Combining top-down and bottom-up dynamics in land use modeling:exploring the future of abandoned farmlands in Europe with the Dyna-CLUE model[J]．Landscape Ecology,2009,24(9):1167-1181.

[46] VELDKAMP A,FRESCO L O.CLUE:a conceptual model to study the conversion of land use and its effects[J]．Ecological Modelling,1996,85(2/3):253-270.

[47] VERBURG P H,SOEPBOER W,VELDKMP A,et al.Modeling the spatial dynamics of regional land use:The CLUE-S model[J]．Environment Management,2002,30(3):391-405.

[48] 鄧祥征，劉紀遠，戰金艷，等.太仆寺旗土地利用變化時空格局的動態模擬[J]．地理研究,2004,23(2):147-156.

[49] 張永民，趙士洞，VERBURG P H.CLUE-S模型及其在奈曼旗土地利用時空動態變化模擬中的應用[J]．自然資源學報,2003(3):310-318.

[50] 蔡玉梅，劉彥隨，宇振榮，等.土地利用變化空間模擬的進展——CLUS-S模型及其應用[J]．地理科學進展,2004(4):63-71.

[51] 陳佑啟.中國土地利用變化及其影響的空間建模分析[J]．地理科學進展,2000,19(2):116-127.

[52] 張永明.科爾沁沙地及其周圍地區土地利用變化的情景分析[J]．自然資源學報， 2004， 19(1):29-37.

[53] 張永明.CLUE-S模型及其在奈曼旗土地利用時空動態變化模擬中的應用[J]．自然資源學報，2003，18(3):310-318.

[54] 甘紅，劉彥隨，王大偉.土地利用類型轉換的人文驅動因子模擬分析[J]．資源科學，2004，26(2):88-93.

[55] 史培軍，陳晉，潘耀忠.深圳市土地利用變化機制分析[J]．地理學報，2000，55(2):150-160.

[56] 王秀蘭.土地利用/土地覆蓋變化中的人口因素分析[J]．資源科學，2000，22(3):39-42.

[57] 朱會義，何書金，張明.環渤海地區土地利用變化的驅動力分析[J]．地理研究，2001，20(6):669-678.

[58] 擺萬奇.大渡河上游地區土地利用/土地覆被變化與驅動力分析[J]．地理科學進展，2004，23(1):71-78.

[59] 焦峰，秦伯強.GIS支持下的小尺度土地驅動力研究——以宜興市湖滏小流域為例[J]．長江流域資源與環境，2003，12(3):205-210.

[60] 李秀彬.土地利用變化的解釋[J]．地理科學進展，2002，21(5):195-203.

[61] PARKER D C,MANSON S M,JANSSEN M A,et al.Multi-agent systems for the simulation of land-use and land-cover change:a review[J]．Annals of the Assosiation of America Geographers,2003,93(2):314-337.

[62] CRAWFORD T W,MESSINA J P,MANSON S M,et al.Complexity science,complex systems,and land-use research[J]．Environment and planning B:Planning and Design,2005,32(6):792-798.

[63] VALBUENA D,VERBURG P H,BREGT A K.A method to define a typology for agent-based analysis in regional land-use research[J]．Agriculture Ecosystems & Environment,2008,128(1/2):27-36.

[64] JANSSEN M A,WALKER B H,LANGRIDGE J,et al.An adaptive model for analyzing co-evolution of management and policies in a complex rangeland system[J]．Ecological Modelling,2000,131(2/3):249-268.

[65] LE Q B.Multi-agent system for simulation of land-use and land-cover change:A Theoretical Framework and its First Implementation for an Upland Watershed in the Central Coast of Vietnam[M].Gottingen:ZFF,2005.

[66] GAUBE V,KAISER C,WILDENBERG M,et al.Combining agent-based and stoke-flow modelling approaches in a participative analysis of the integrated and system in Reichraming,Austria[J]．Landscape Ecology,2009,24:1149-1165.

[67] MONTICINO M,ACEVEDO M,CALLICOTT B,et al.Multi-agent model of human values and land-use change[C]//Proceedings of the fifth IASTED International Conference on Modelling,Simulation and Optimization,2005:279-284.

[68] ACEVEDO M F,CALLICOTT J B,MONTICINO M,et al.Models of natural and human dynamics in forest landscapes:Cross-site and cross-culture synthesis[J]．Geoforum,2008,39:846-866.

[69] MESSINA J P,EVANS T P,MANSON S M,et al.Complex systems models and the management of error and uncertainty[J]．Journal of Land Use Science,2008,3:11-25.

Study on the Analysis Method of Land Use / Cover Change

CHOU Qing-chuan1, HU Tao1, HE Shi-yu1, SHI Xiu-hua2*

(1. College of Life Science, Shenzhen University, Shenzhen, Guangdong, 518060; 2. Golf College of Shenzhen University, Shenzhen, Guangdong 518060)

Land use /cover change is considered to be one of the important components and the main driving factors of global change, and increasingly concern by international organizations and countries around the world, has become the forefront of the international global change research. Domestic and foreign scholars from different angle build a large number of models, and in-depth understanding the causes and process of land use/cover change, predict the future development trends. The land use dynamic change model was studied mainly from land use speed, land use degree, land use change trend and the landscape index of land use and other aspects, which highlights the current most widely used model, at the same time points out the future development direction of the model.

Model; Dynamic change; Land use / cover change; Prediction

深圳灣紅樹林結構調控及修復技術研究與示范項目(2013KJCX011-04)。

丑慶川(1988-)，男，甘肅慶陽人，碩士研究生，研究方向：景觀與濕地資源。*通訊作者，教授，博士，從事環境保護研究。

2015-03-02

S-03

A

0517-6611(2015)10-301-04