平天湖流域土地利用空間異質性特征分析

張姣姣，胡和兵，姜志豪

（1.池州學院 地理與規劃學院，安徽 池州 247000；2.阜陽師范大學 歷史文化與旅游學院，安徽 阜陽 236037）

土地資源是人類生產和生活必不可少的自然資源，土地利用/覆被變化（LUCC）一直是相關研究的熱點[1]。早期主要對土地利用的數量變化進行研究[2]。隨著3S技術的快速發展，在更大時間和空間尺度上獲取土地利用數據成為可能，土地利用空間分異研究也日益受到關注，許多學者針對縣域、市域、省域、全國等不同尺度，以及不同流域等開展了大量研究[3-7]，成果豐碩。但針對水源地流域進行土地利用空間異質性的關注較少。

流域土地利用空間異質性是指流域土地利用類型組成特征及其在空間分布上的不均勻性和復雜性。相比于區域，流域是地球表層相對獨立的自然綜合體。以流域為單元進行綜合研究與管理是實現資源開發與環境保護相協調的最佳途徑[8]。許多研究表明，在點源得到較好控制后，人類活動影響下的土地利用結構與空間格局變化是導致流域地表水體水質惡化的主要原因[9]。其中空間格局決定了流域內的能量流動和營養物質循環過程，進而影響著流域環境的變化[10]。

鑒于此，以平天湖流域為研究區，借助遙感以及GIS空間分析方法，依據流域DEM，在將整個流域劃分為石馬河流域、白沙河流域、碧山河流域以及其他部分等4個子流域的基礎上，對平天湖流域土地利用結構、多樣性、有序性以及強度進行了空間異質性分析，可為平天湖流域的土地利用優化以及水源地水質保護提供科學依據。

1 研究方法及數據來源

1.1 研究區域概況

平天湖流域位于安徽省池州市主城區東南部（圖1），其主體部分位于齊山-平天湖國家AAAA級旅游景區內，流域面積約74 km2。屬亞熱帶季風性濕潤氣候區，年平均溫度16.1℃，平均年降雨量在1400-1700mm，流域整個地勢南高北低，最高海拔在714 m。流域南部為中低山區，土地利用類型以林地為主，中北部為丘陵區，農用地及居民建設用地是其優勢土地利用類型。平天湖現為池州城市備用飲用水水源地，水域面積12 km2，西臨清溪河，北接秋浦河故道后匯入長江，石馬河、白沙河、碧山河為其3條入湖支流。

圖1 平天湖流域地理位置及子流域劃分

1.2 研究方法

1.2.1 吉布斯—馬丁多樣化指數 研究土地利用空間差異性可以通過對各研究區土地利用類型的多樣化程度進行比較和分析，主要分析的目的在于確定研究區內土地利用類型的齊全程度和多樣化水平等[11]，其結果可基于吉布斯—馬丁多樣化指數來分析[6]。

1.2.2 信息熵 信息熵最早從熱力學中引用過來，信息熵越大，系統就越復雜，信息熵越小，系統就越簡單[12]。目前，許多對土地利用的研究中都運用到信息熵，土地利用信息熵的大小反映了某一地區土地利用系統的有序程度[8]。

1.2.3 土地利用程度綜合指數 土地利用程度綜合指數大小反映了某地區土地利用的程度，按土地利用程度分析的方法，將土地利用分為4級。第1級為未利用地；第2級為草地、林地、水域；第3級為耕地；第4級為居民建設用地[13]。

1.3 數據預處理

基于2017年遙感影像，利用ENVI5.3進行輻射校正、大氣校正、影像融合、鑲嵌與裁剪等一系列處理后，在ArcGIS10.3平臺上進行目視解譯。參考中國科學院土地分類體系，結合研究區土地利用特點，將研究區域地表類型分為6類：林地、耕地、草地、建筑用地、未利用地和水體，并采用GPS進行野外定點核查，解譯精度達87.56%，滿足研究要求。基于流域DEM（30 m，來源于ASTER GDEM），利用ArcGIS10.3的水文分析功能，進行平天湖流域邊界提取，并劃分為白沙河流域、石馬河流域、碧山河流域以及其他4個子流域。

2 結果與分析

2.1 平天湖流域土地利用結構分析

平天湖流域及其4個子流域均呈現以林地為主的土地利用結構（圖2）。在整個流域中，林地面積占比達57%，其次是耕地，占18%。各子流域占優勢的土地利用類型仍是林地。其中石馬河流域林地占比最大，達69%，其他3個子流域林地占比差別不大。其次是耕地，耕地的占比在石馬河流域、白沙河流域以及碧山河流域間差別不大，均約占各自流域的20%左右。在各子流域中，建設用地的占比均較小。建設用地占比最高的石馬河流域也僅為16.22%。平天湖流域水域的占比最小，在水域分布最多的碧山和流域也僅達2.4%。

圖2 平天湖流域及子流域土地利用結構

各子流域土地利用類型的空間分布存在較大差異。石馬河流域建筑物主要分布在北部站前區，南部村落建筑物成離散分布、小型聚集，建設用地面積相對較廣。耕地主要分布在南部，圍繞著村落分布，耕地周圍以山地為主，其土地利用類型以林地和草地為主。白沙河流域建筑物主要分布在流域的中部，耕地分布在北部和中部，其占比比建筑用地多，流域的南部以山地為主。碧山河流域建筑物主要分布在流域北部教育園區一帶。流域內耕地圍繞著建筑用地分布，集中在流域的北部和南部，其林地和草地主要分布在北部和東部，部分位于南部的山區。其他部分建筑物主要分布在東北部和南部，占地面積相對較多，耕地主要分布在北部，圍繞平天湖分布，主要是平天湖提供了充足的水資源，有利于農作物的生長。

2.2 平天湖流域土地利用多樣化分析

石馬河流域、白沙河流域、碧山河流域和其他部分土地利用多樣化指數分別為0.70、0.48、0.65和0.68，其中石馬河流域土地利用多樣化指數最高，為0.70，其主要原因是石馬河流域包含了池州站前區、包家山、區域內居民地、魚塘和田地等，土地利用類型多且占比相對較為均衡。白沙河流域土地利用多樣指數最低，為0.48，主要是由于白沙河流域林地的優勢過于明顯，占流域面積達68.74%，而其他土地利用類型分布過少導致。

通過創建平天湖流域采樣格網，計算每個格網的土地利用多樣化指數，將其作為格網中心點的土地利用多樣化指數，以每個格網中心點為基礎運用克里金插值法對平天湖流域進行插值分析，裁剪后得到平天湖各流域土地利用多樣化等級圖（圖3）。

圖3 子流域土地利用多樣化空間差異等級圖

平天湖各流域每個格網的土地利用多樣化指數都在0和0.85之間，以0.17為間隔，將平天湖流域土地利用多樣化分為5個等級。為更加清楚地分析和比較平天湖各流域土地利用多樣化，在克里金插值分析的基礎上統計了平天湖各流域土地利用多樣化每一級別的面積（表1）。

表1 子流域土地利用多樣化等級面積統計

各子流域土地利用多樣化程度總體表現為：石馬河流域、碧山河流域較高，而白沙河流域與其他部分較低。石馬河流域中等土地利用多樣化面積占比較大，為48.91%，較高土地利用多樣化面積占20.25%，二者共計達69.16%，在4個子流域中最高。高土地利用多樣化和較高土地利用多樣化集中分布在石馬河流域的北部，這主要是由于北部以池州站前區為主。隨著齊山花園、徽商江南世家等眾多樓盤的開發，致使這一區域土地利用多樣化提高。碧山河流域中等土地利用多樣化面積占比為46.43%，較高土地利用面積占比也較高，為20.58%。教育園區高和較高土地利用多樣化集中分布在碧山河流域中部的碎嶺一帶，G318國道跨牧之路（碎嶺）互通立交橋的修建以及周邊農林養殖業等，使得該地段的土地利用多樣化程度提升。

白沙河流域低土地利用多樣化面積占比最大，占比40.18%，這是由于石馬河流域南部土地類型以山地為主，土地利用類型以單一的林地為主。與石馬流域相同，高土地利用多樣化和較高土地利用多樣化集中分布在流域的北部，主要由站前區的發展、G318國道的修建所致。其他部分的中等土地利用多樣化面積占比最大，占比54.04%，但較高等級的面積僅占13.30%，土地利用多樣化較高集中分布在平天湖風景區內村莊一帶，農戶生計直接影響土地利用的類型及其變化。

2.3 平天湖流域土地利用有序性分析

石馬河、白沙河、碧山河和其他部分土地利用信息熵值分別為1.43、0.97、1.33和1.37，石馬河土地利用信息熵最高，土地利用有序性相對較低，白沙河流域信息熵最低，土地利用有序性相對較高，穩定性較強。

利用采樣格網中心點的土地利用信息熵進行克里金插值分析后，通過裁剪獲得平天湖各流域土地利用信息熵等級圖（圖4）。

圖4 子流域土地利用信息熵空間差異等級圖

平天湖各流域采樣格網的土地利用信息熵都在0和1.69之間，以0.338為間隔，將土地利用信息熵分為五個等級，再在克里金插值分析基礎上統計各等級面積（表2）。

表2 子流域土地利用信息熵等級面積統計

石馬河流域土地利用信息熵主要集中在較低信息熵和中等信息熵之間，中等信息熵面積占比最大，占比44.97%，高、較高信息熵主要分布在北部池州站前區齊山新村、清溪雅居安置房一帶，這與站前區建設、318國道池州至殷匯段一級公路改建工程等項目拆遷安置有關，促進了站前區的面積在不斷的擴張，人口流動性大，區域內土地利用結構處于相對較高的無序化狀態。白沙河流域土地利用信息熵主要集中在低信息熵和中等信息熵之間，低信息熵面積占比最大，占比42.12%，高、較高信息熵主要分布在北部。碧山河流域土地利用信息熵主要集中在低信息熵和中等信息熵之間，中等信息熵面積占比最大，占比53.04%，高信息熵和較高信息熵主要分布在中部。白沙河流域、碧山河流域受G318國道修建的影響，調節土地利用結構。其他部分土地利用信息熵集中分布低信息熵和中等信息熵之間，中等信息熵面積占比最大，占比55.99%。

2.4 平天湖流域土地利用程度分析

石馬河、白沙河、碧山河和其他部分土地利用綜合指數分別為245.01、227.78、237.96和225.21，其中石馬河流域的土地利用程度最高，其他部分土地利用程度最低。

利用采樣格網中心點的土地利用程度進行克里金插值分析后，通過裁剪獲得平天湖各流域土地利用程度等級圖（圖5）。

圖5 子流域土地利用程度空間差異等級圖

通過分析計算得知平天湖各流域每塊采樣網格的土地利用程度指數都在100和400之間，以60為間隔，將平天湖流域土地利用程度分為5個等級。在普通克里金插值分析基礎上統計平天湖各流域等級面積（表3）。

表3 子流域土地利用程度等級面積統計

石馬河流域土地利用程度低、較低的區域主要分布在平天湖旁的包家山、蓮花山，土地利用程度中等等級面積占比最大，占比69.31%，較高和高土地利用程度在流域北部的站前區、中部和南部均有分布。白沙河流域土地利用程度主要集中在較低土地利用程度和較高土地利用程度之間，較低土地利用程度等級面積占比最大，占比50.93%，主要是白沙河流域南部集中分布著廣闊的林地。白沙河流域較高和高土地利用程度分布在北部。碧山河流域土地利用程度主要集中在較低土地利用程度和較高土地利用程度之間，較低等級的土地利用程度區域為流域林地一帶，中等土地利用程度等級面積占比最大，占比60.40%，主要在教育園區、周邊村落一帶。

其他部分較低土地利用程度和中等土地利用程度之間，較低土地利用等級面積占比最大，占比66.33%，較高和高土地利用程度的區域主要分布在西部的平天湖度假區、東北部的高新技術產業開發區一帶。

3 結論與建議

平天湖整個流域及其子流域均是以林地為優勢的土地利用結構。整個流域林地面積占比達57%，其次是耕地，占比達18%，而其他土地利用類型相差不大；在4個子流域中，以白沙河流域林地的占比最大，達68.74%，石馬河占比最小，為46.12%。平天湖流域的土地利用多樣性及有序性在空間上存在一定差異。其中以白沙河流域土地利用多樣性及有序性最低，分別為0.48、0.97，而其他3個子流域間差異不大。平天湖流域土地利用多樣性、有序性高的區域主要分布在石馬河流域、白沙河流域、其他部分的北部、碧山河流域的中部。平天湖流域的土地利用程度整體較低。其中石馬河流域和碧山河流域土地利用程度相對較高，且以中等土地利用強度為主，面積占比分別為69.31%、60.40%，而高土地利用程度占比很少，均在3%以下。

土地利用程度量化了人類活動對于土地系統的影響程度，定量地描述了該地區土地利用的綜合水平和變化趨勢。隨著池州站前區、高新技術開發區建設的不斷推進，平天湖流域土地資源將遭受空前的壓力。要注重開發規劃與發展規劃相協調，注重山、水、林、田、湖、草的保護。由于區域內土地利用類型多樣，需要針對不同類型的土地實施區別開發措施。對強土地利用程度分布區，應嚴格控制建筑用地的密度，加強區域綠化和景觀建設。平天湖作為池州市備用水源地，對土地利用程度較低的待開發區域，需要堅持生態保護和修治并舉，實現人與自然和諧共生。這樣有助于更加合理的開發和利用土地資源，促進區域可持續發展。