官司河流域景觀結構變化及其分形研究

陳俊華，何政偉，龔固堂，許輝熙，朱志芳，吳雪仙，慕長龍

（1.成都理工大學，四川 成都 610059；2.四川省林業科學研究院，四川 成都 610081；3.四川建筑職業技術學院 測量工程研究所，四川 德陽 618000）

0 前言

景觀結構是指斑塊、廊道和基質的時空配置形成的鑲嵌格局［1］。斑塊作為景觀結構的基本要素，其特征是由大小和形狀來表示的。其中由于斑塊的形狀復雜多樣，通常情況下難以直接測量，因此，一般采用各種指數來描述［2］。分形理論是20世紀70年代由法國數學家 Mandelbrot［3］提出的，用以定量表達自然界中傳統歐式幾何學不能描述的復雜而有規則的幾何現象，揭示了非線性系統中有序和無序的統一、確定性和隨機性的統一問題。通過分維的計算，可以對分形幾何對象的復雜性程度進行度量。景觀研究中可以在多尺度下觀察景觀，從而對復雜的邊界和斑塊形狀進行量化計算［4］。對于分形理論的研究，國內多采用從種群和群落角度計算分維或從靜態的角度來研究景觀格局與分維的關系［5～16］。

作者在本文以四川省綿陽市游仙區官司河小流域作為研究對象，在“3S”技術的支持下，應用分形理論定量地分析該區域各種景觀要素結構的復雜性與穩定性，揭示其景觀空間結構的分形特征和景觀結構變化的分維動態特點，可為區域生態環境保護以及可持續發展提供科學依據。

1 研究區概況

官司河小流域位于四川省綿陽市游仙區新橋鎮，處于龍門山前緣向盆地的過渡地帶，屬淺～深切割的丘陵地貌。氣候屬北亞熱帶濕潤季風氣候，年均溫16.1℃，年均降水量921.0mm，年均相對濕度79%，無霜期272d［17］。土壤主要有發育在第四系松散堆積物上的老沖積黃壤和發育在白堊系泥巖和砂巖風化物上的紫色土和灰白砂土。森林植被主要有針葉林、闊葉林、針闊混交林、經濟林等，森林覆蓋率達35.2%，主要森林類型有柏木林、馬尾松林、櫟類林、榿木等［17～19］。

2 研究方法

2.1 研究數據源及土地利用GIS數據庫構建

研究采用1995的森林資源分布圖和2005年的IKONOS影像數據。根據研究區實際情況，將研究區土地利用景觀類型劃分為耕地（cultivated land）、建筑用地（building land）、道路（road）、水域（water）、柏木林（Biota orientalis）、馬尾松林（Masson Pine forest）、櫟林（Quercus forest）七個類型。

在1∶10 000林相圖（2005年Ⅱ類調查成果）和1∶10 000地形圖的輔助下，利用GPS進行實地調查采樣，編制影像目視解譯標志數據庫，在地理信息系統軟件ArcGIS 9.3環境中對IKONOS影像進行人工目視解譯，得到2005年的土地利用GIS數據庫。

在GIS軟件（ArcGIS 9.3）環境中對1995年森林資源分布圖進行人機交互屏幕矢量化，得到1995年的土地利用GIS數據庫。

利用GIS軟件的制圖功能，制作兩個時期的土地利用景觀類型圖（見圖1）。

2.2 數據計算

景觀中的斑塊大都呈不規則的圖形，景觀鑲嵌體是由各種景觀要素斑塊在空間上鑲嵌分布且有機結合在一起形成的。景觀鑲嵌體是自然界中最典型的分形幾何體，并在特定時間呈二維空間分布。因此景觀的動態變化維數可采用二維歐氏空間的分數維公式來計算，從而描述景觀斑塊的復雜性。并結合景觀多樣性（Diversity index）指數、優勢度（Dominance）和景觀破碎度（Landscape fragmentation）等指標進行定量化描述［20、21］。

2.2.1 分維值

利用董連科［22］提出的公式。

式中 D表示分維數；P（r）指某一斑塊的周長；A（r）為同一圖斑的面積；C為常數。

D的理論值區間［1，2］，D值越大，表示該土地利用類型越復雜。當D＝1時，表示該類型斑塊的形狀為正方形；當D等效于“2”時，表示該類型斑塊的形狀最復雜；當D＝1.5時，則表示該類型處于一種類似于布朗運動的隨機狀態，即最不穩定狀態。D值越接近“1.5”表示該要素越不穩定。

圖1 官司河流域土地利用類型圖Fig.1 The map of land－use types of Guansihe watershed

2.2.2 穩定性指數

通常采用徐建華［23］等定義的景觀要素的穩定性指數要SI。

SI＝｜D－1.5｜

式中 SI指景觀類型在其形態復雜狀態下的一種穩定性，D為分維值。

SI值越大，表示該景觀類型越穩定。

2.2.3 景觀多樣性指數 （H）

區域內不同景觀類型分布的均勻化和復雜化程度，通常用景觀多樣性指數表示。在本文中，作者采用Shannon－Wiener指數，公式為：

式中 Pi指第i類景觀所占的面積比例；m指景觀類型的數目。

H值越大，表示景觀要素類型愈豐富，即景觀多樣性越大，其最大值為：Hmax＝logm。

2.2.4 景觀優勢度 （D）

優勢度表示少數嵌塊體在景觀中的支配程度，其公式為：

式中 D為景觀優勢度。當D值小時，反映該景觀是由多個面積大致相當的嵌塊體所組成；當D值大時，表示景觀只受一個或少數幾個嵌塊體類型支配。

由此可以看出，景觀優勢度與多樣性指數剛好相反，對于景觀類型數目相同的不同景觀，景觀多樣性指數越大，其優勢度越小。

2.2.5 景觀破碎度 （I）

破碎度表征景觀被分割的破碎程度，反映景觀空間結構的復雜性，在一定程度上反映了人類對景觀的干擾程度，其公式為：

式中 I為景觀破碎度，I值越大表示破碎化程度越高；Ni為第i類景觀斑塊數；A為景觀總面積。

3 結果與分析

3.1 不同時期景觀格局特征

根據兩期衛星遙感影像解譯出官司河流域各景觀類型屬性數據，在Excel 2003下進行統計分析面積（見下頁表1），并制作直方圖（見下頁圖2）。

分析兩個時期土地利用類型面積變化可以看出（見表1），兩個時期耕地均占絕對優勢，分別占整個流域面積的68.77%和51.47%；其次是林業用地面積，分別占22.17%和37.73%。從1995年～2005年十一年間，土地利用類型變化較大。從土地利用類型一級分類來看，林地、道路、水域的面積不斷擴大。其中，道路用地增加的幅度最大，增加了113.06%；其次是林地，增加為70.22%，水域增加了33.01%。而耕地和建筑用地面積卻不斷減少，其中耕地減少了25.16%，這是因為1995年～2005年這十一年期間，25°以上的坡耕地被退耕地還林的緣故。建筑用地面積減少是因為近些年來實施農民聚居以及部分農民外地務工在城市購買房屋后遷離該區域。有林地面積和森林植被類型也發生了較大變化。馬尾松、柏樹、櫟樹林地面積大幅度增加，其中面積增加最多的是馬尾松，十一年來增加了177.78hm2，增幅為64.36%；其次是柏木林，增加了131.5hm2，增幅為79.96%。櫟林面積增加了3.39hm2，增幅為73.38%。

從兩個時期各景觀類型結構可以看出（見圖2），兩個時期建筑用地、水域和耕地的多樣性指數最大，優勢度也較為明顯，而各類林業用地多樣性指數和優勢度均較小。從破碎度來說，建筑用地的破碎度較大，林業用地的破碎度在兩個時期也相對較小。

3.2 景觀要素的分維及變化

在地理信息系統軟件ArcGIS 9.3支持下生成各景觀要素斑塊的面積和周長，在SPSS 15軟件下進行擬合，獲取回歸模型，得出分形維數（如下頁表2所示）。從表2可以看出，官司河流域各景觀的空間分布在一定尺度范圍內均具有分形性質（0＜D＜2），大多呈顯著相關（相關系數在0.9以上）。分維能夠定量地描述景觀結構的特征參數，且不隨測量尺度的改變而變化。分維的大小反映了一定區域內景觀結構的復雜程度，當分維值越大，表示景觀要素斑塊越不規則、越復雜。不同時段的分維變化可以表征景觀斑塊在時間維上的空間變化。

（1）1995年官司河流域各景觀類型分維數按由大到小排列為：柏木林＞馬尾松林＞道路＞耕地＞櫟林＞水域＞建筑用地。柏木林的分維數最大，表明在官司河流域，柏木林的分布最為復雜和不規則。建筑用地和水域的分維值較低，接近于“1”，說明這兩種類型分布規模較小，斑塊較為規則。

（2）2005年各景觀類型分維數按由大到小排列為：馬尾松林＞柏木林＞櫟林＞耕地＞道路＞建筑用地＞水域。其中，馬尾松林在十一年間平均斑塊面積減少，破碎化程度增加，分維數上升到第一位。而水域斑塊數減少，平均斑塊面積增加，分維數下降到最后一位。

從穩定性指數來看，兩個時期建筑用地和水域的該指數均最大，表明這兩種景觀類型最為穩定，而馬尾松林的該指數最小，表明該類型最不穩定。整個林業用地的穩定性指數普遍小于其它景觀類型，表明林業用地處于不穩定狀態。

從兩個時期分維數的分布來看，林業用地和耕地的分維數均大于其它景觀類型，表明這兩種類型在官司河流域處于主導地位，是基質。從分維數的動態變化來看，除道路和水域這兩種景觀類型減少外，其它景觀類型都增加，其中建筑用地增加得最多，其次是櫟林。表明了人類對這幾種類型的利用程度更加復雜化。

表1 官司河流域景觀結構變化情況（1995年～2005年）Tab.1 The changes of landscape structure of Guansihe watershed from 1995to 2005

圖2 官司河流域景觀格局特征Fig.2 The landscape patterns of Guansihe watershed

表2 官司河流域景觀要素斑塊分維數Tab.2 The changes of fractal dimension on the landscape of Guansihe watershed from 1995to 2005

4 討論

官司河流域屬于典型的四川盆地丘陵區，人口眾多，人類活動較為頻繁，對土地的利用也較為復雜，因而該區域耕地和林業用地的分維值高于其它類型，但穩定性卻明顯低于其它類型。從1995年到2005年這十一年間，由于經濟的發展和退耕還林的實施，耕地面積明顯減少，各林業用地面積明顯增加，人類對林業用地的利用程度復雜化，使得各林業用地類型的分維值上升，但穩定性卻降低。這說明從一定程度上該區域的森林結構不合理，應該進行調整。

研究表明，景觀分維在度量景觀及其變化上具有綜合特征，分維的大小能夠反映區域內景觀斑塊的不規則程度和復雜程度。“3S技術”在快速度量區域土地利用動態變化方面具有較大優勢，利用GPS作為野外輔助調查工具，RS作為獲取土地利用詳細信息，GIS強大的空間分析功能，可獲取各景觀類型的結構特征，從而提供了一個快速的量測土地利用景觀復雜性和穩定性的工具。

［1］FORMAN RTT，GODRON M.景觀生態學［M］.肖篤寧，譯.北京：科學出版社，1990.

［2］FORMAN R.Land Mosaics：the Ecology of Landscapes and Regions［M］.Cambridge：Cambridge U－niversity Press，1995.

［3］MANDELBROT B B.The fractal geometry of nature［M］.New York：W H Freeman，1982.

［4］KRUMMEL J R，GARDNER R H，SUSIHARA G，et al.Landscape pattern in a disturbed environment［J］.Oikos，1987，48（3）：321.

［5］馬克明，祖元剛.植被格局的分形特征［J］.植物生態學報，2000，24（1）：111.

［6］鄧紅兵，周永斌，王慶禮，等.三峽庫區次生柏木林中物種分布格局的分形特征［J］.應用生態學報，1999，10（5）∶518.

［7］吳大千，張治國，王玉芳.暖溫帶石灰巖山地森林種群格局的分形特征研究［J］.山東林業科技，2006（5）：5.

［8］余世孝，郭濼.山東泰山地區景觀結構變化及其分形分析［J］.生態學報，2005，25（1）：129.

［9］王兆杰.格氏栲自然保護區生態旅游景觀評價及區劃研究［D］.福州：福建農林大學，2007.

［10］畢如田，白中科，葉寶瑩，等.安太堡露天煤礦排土場景觀結構及其分形研究［J］.中國農學通報，2007，23（4）：332.

［11］梁帥.人類活動對上海濱海地區景觀環境影響對比研究［J］.能源與環境，2010（2）：42.

［12］王水獻，董新光，寇文件.新疆焉耆盆地綠州景觀分形特征及驅動機制分析［J］.水土保持研究，2007（12）：256.

［13］朱志芳，龔固堂，許輝熙，等.基于RS／GIS的平通河流域土地利用格局變化及其分形分析［J］.四川林業科技，2009，30（1）：24.

［14］李玲慧.城鎮化進程下的澤庫縣生態環境脆弱性分析與評價［D］.西安：西安建筑科技大學，2011.

［15］晉秀龍，陸林.旅游生態學研究方法評述［J］.生態學報，2008，28（5）：2343.

［16］李義玲，喬木，楊小林，等.干旱區典型流域近30年土地利用／土地覆被變化的分形特征分析——以瑪納期河流域為例［J］.干旱區地理，2008，31（1）：75.

［17］陽小成.李旭光.葉志義.川中半自然型防護林演替過程中的物種多樣性［J］.應用與環境生物學報.2002，8（2）：127.

［18］陳俊華，龔固堂，朱志芳，等.小流域防護林體系的空間對位配置［J］.林業科學，2012，48（2）：38.

［19］駱宗詩，向成華，慕長龍，等.綿陽官司河流域主要森林類型凋落物含量及動態變化［J］.生態學報，2007，27（5）：1772.

［20］鄔建國.景觀生態學——格局、過程、尺度與等級［M］.北京：高等教育出版社，2000.

［21］李義玲，喬木，楊小林，等.干旱區典型流域近30年土地利用／土地覆蓋變化的分形特征分析－以瑪納斯河流域為例［J］.干旱區地理，2008，31（1）：75.

［22］董連科，分形理論及其應用［M］.沈陽：遼寧科學出版社，1991.

［23］徐建華，艾南山，金烔，等，北干旱區景觀要素鑲嵌結構的分形研究－以黑河流域為例［J］.干旱區研究，2001，18（1）：36.