劉利霞 王威 王小煥 王芳 楊英 黨永峰

摘要：基于（GIS，運用熵權(quán)—TOPSIS模型，充分考慮京西林場特殊的自然條件、社會經(jīng)濟狀況，并遵循科學性、綜合性、可行性、動態(tài)性、區(qū)域性等原則，選取10個指標因子，以小流域為單位對京西林場資源環(huán)境承載能力進行了評價，通過對綜合得分的排序，共分為了3個級別，分別為好（0.32

關(guān)鍵詞：京西林場;資源環(huán)境;承載力;熵權(quán)TOPSIS模型

中圖分類號：S7

文獻標識碼：A

文章編號：1674-9944（2018）6-0129-04

1引言

區(qū)域承載力是區(qū)域條件在維持發(fā)展的情況下，培育健康穩(wěn)定、優(yōu)質(zhì)高效的森林，滿足生態(tài)、社會、經(jīng)濟發(fā)展需要，意義重大，是區(qū)域規(guī)劃的基礎(chǔ)。從承載力的角度研究區(qū)域資源、環(huán)境、生態(tài)問題，一直是地理學、生態(tài)學研究的熱點問題之一。目前關(guān)于承載力的研究在區(qū)域尺度上多集中于省級、市級、區(qū)縣級，在研究方法的選擇與應(yīng)用上主要是采用層次分析法、主成分分析等方法。過于主觀，盡管現(xiàn)有關(guān)于資源、環(huán)境與經(jīng)濟作用機制及人類活動相互作用規(guī)律這方面取得相應(yīng)研究成果，但現(xiàn)有文獻對縣級以下尤其是林場尺度的研究涉及甚少。因此對林場區(qū)域資源環(huán)境評價研究具有重要意義，為此本研究在總結(jié)以往研究的基礎(chǔ)上，結(jié)合京西林場特有的自然環(huán)境、區(qū)域經(jīng)濟條件，對其環(huán)境承載力進行評價，以期為轉(zhuǎn)型期山區(qū)林場建設(shè)提供參考。

2區(qū)域概況

京西林場位于北京西部太行山深處，地跨門頭溝區(qū)和房山區(qū)，林場總部位于門頭溝區(qū)人民政府附近，共有大臺、大安山、長溝峪、珠窩林、雁翅、河南臺、二斜井等7個林區(qū)，下設(shè)石房溝、八二零、玉皇廟、北港溝、曹家鋪、大安山、長溝峪、齋堂山8個分場，總面積17.44萬畝，地理位置介于東經(jīng)115°40'-116°38'，北緯39°40'～40°10'。林區(qū)山體屬太行山系小五臺山余脈，海拔200～1610m。氣候?qū)贉貛Т箨懶约撅L氣候，四季分明，無霜期150～180d，年平均氣溫7～10℃，全年平均氣溫比北京市區(qū)低5～8℃，尤其夏季（7～9月）的平均氣溫比市區(qū)低6℃，年平均降雨量500--600mm;溝深坡陡，林區(qū)屬于海河水系上游的永定河流域和大清河流域，河流以季節(jié)性河流為主;土壤類型主要為褐土、棕壤土;場內(nèi)森林資源豐富，其中，林地16.74萬畝，非林地0.70萬畝，森林覆蓋率為31.3%（圖1）。

3數(shù)據(jù)來源與研究方法

3.1數(shù)據(jù)來源

文中所用數(shù)據(jù)主要包括：1：1萬地形圖、1：1萬土地利用圖、2016年2.5m分辨率SPOT-5衛(wèi)星影像、2017年京西林場二類調(diào)查數(shù)據(jù)庫、1：1萬道路分布圖、1：1萬DEM數(shù)據(jù)、1：1萬行政邊界數(shù)據(jù)等相關(guān)數(shù)據(jù)，均來自國家林業(yè)局調(diào)查規(guī)劃設(shè)計院（圖2）。

3.2熵權(quán)—TOPSIS模型

熵權(quán)—TOPSIS模型的優(yōu)勢在于其客觀性和可操作性，對樣本需求不大，且結(jié)果合理。計算步驟如下：

（1）構(gòu)建評價指標體系矩陣（S）。假設(shè)被評價對象有，m個，每個被評價對象有個，則矩陣為：

（7）計算目標對象與最優(yōu)解最劣解間的歐氏距離，計算公式如下：

（8）計算綜合評價指數(shù)，計算公式如下：

公式（9）中，Wi越大表示評價對象越優(yōu)。

4評價指標體系構(gòu)建

京西林場原為北京礦務(wù)局（2001年更名北京京煤集團有限責任公司）所屬林場，于1964年成立，旅游資源、礦產(chǎn)資源較為豐富，但場內(nèi)水資源匱乏，經(jīng)濟基礎(chǔ)薄弱，在選取評價指標時，充分考慮了其特殊的自然條件、社會經(jīng)濟狀況，并遵循科學性、綜合性、可行性、動態(tài)性、區(qū)域性等原則。本文共選取10個指標因子（表1），基于熵權(quán)TOPSIS模型，以小流域為單位對京西林場資源環(huán)境承載能力進行評價。

指標獲取方法如下：

人口密度=場內(nèi)人口/區(qū)域面積;森林覆蓋率=（有林地面積十國家特別規(guī)定的灌木林地）/區(qū)域土地面積×100%;公路網(wǎng)密度一區(qū)域公路長度/區(qū)域土地面積.

土地利用狀況表達式如下：

公式（10）中：n表示土地利用狀況指數(shù)（表1），LU表示土地利用狀況指數(shù)，Si表示為第i級土地利用狀況分級指數(shù)（表2），M，表示第i級土地利用狀況面積百分比。

海拔和高程根據(jù)研究區(qū)1：1萬DEM影像及京西林場二類調(diào)查數(shù)據(jù)，運用Arcgis軟件，得出研究區(qū)海拔和高程數(shù)據(jù)。

植被蓋度一場內(nèi)植被面積/土地總面積×100%

蓄積保有量以植被生長的Landsat ETM+衛(wèi)星數(shù)據(jù)及二類調(diào)查數(shù)據(jù)中的地類字段數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)，運用Arcgis、ENVI軟件對遙感影像進行幾何校正、配準及相關(guān)數(shù)據(jù)處理。

樹種多樣性指數(shù)和傳統(tǒng)的林分變量相比，樹種多樣性反映了森林結(jié)構(gòu)更多的細節(jié)信息，揭示了森林客觀存在的樹種和結(jié)構(gòu)的多樣性，計算表達式采用Shannon-Wiener指數(shù)（SHI）：

公式（11）中，N表示樹種的個體數(shù);Ai表示第i個樹種的個體占所有個體數(shù)的比例。

生態(tài)環(huán)境敏感性指數(shù)=N×（0.05×一般敏感性土地面積+0.25×較敏感性土地面積）×100/區(qū)域總面積 （12）

公式（12）中，N表示敏感性指數(shù)歸一化后的系數(shù)，取1。

6結(jié)果與分析

資源環(huán)境承載力評價指標權(quán)重最優(yōu)劣值排序。通過熵權(quán)-TOPSIS模型，計算得出林場每個流域的綜合分值（如表3），通過對得分進行降序排序，將區(qū)域資源環(huán)境承載力分為3個等級，分別為好（0.32

通過表4可看出：京西林場17個小流域中，二斜井流域及齋堂山流域得分最高，河南臺流域得分最低，分別為0.59、0.14。

二斜井流域、位于齋堂山分場的齋堂山流域、位于石房溝分場及代管林地分場和北港溝分場及八二零分場的清水澗流域、大安山分場的花園溝流域、曹家鋪分場的曹家鋪流域、長溝峪分場的賈島峪溝流域、石房溝分場的櫻桃溝流域，資源環(huán)境承載力好，開發(fā)潛力相對較大;北港溝分場的北港溝流域和南港溝流域、大安山分場的大安山流域、石房溝分場的石房溝流域、長溝峪分場的長溝峪流域、八二零分場的五二五溝流域、大臺分場的大臺流域，資源環(huán)境承載力等級為較好，開發(fā)潛力屬中等水平;石房溝分場的雁翅流域及珠窩流域、石房溝分場的河南臺流域，資源環(huán)境承載力一般，開發(fā)力度較大，應(yīng)加大生態(tài)保護和植被恢復(fù)力度（表3、表4）。