雅安4·20震后寶興河流域農(nóng)村居住區(qū)適宜性評價

米琳迪, 周 英, 王 青

(西南科技大學 環(huán)境與資源學院, 四川 綿陽 621010)

?

雅安4·20震后寶興河流域農(nóng)村居住區(qū)適宜性評價

米琳迪, 周 英, 王 青

(西南科技大學 環(huán)境與資源學院, 四川 綿陽 621010)

在汶川—蘆山兩次強震作用下，先天脆弱與高風險并存的武都—瀘定一線山區(qū)成為受災最為嚴重的地理單元。災區(qū)生境重建是4·20災后重建的重點、難點，也是維護災區(qū)社會穩(wěn)定的關鍵之一。以寶興河流域農(nóng)村居住地為研究對象，選取了地形、氣象、植被、地質(zhì)災害、大熊貓棲息地世界遺產(chǎn)等影響適宜性的關鍵因子，開展農(nóng)村居住地環(huán)境適宜性評價。結果表明：(1) 山區(qū)農(nóng)村居住地適宜性分布呈現(xiàn)出明顯的海拔、坡度、坡向指向，且集中分布在世界遺產(chǎn)保護區(qū)及其外圍。海拔2 900 m以下、坡度小于30°且陽坡、半陽坡地帶是適宜性及次適宜性的主要分布區(qū)，具有沿河流及兩側(cè)集聚分布的特征，說明地形條件是影響適宜性分布的主要因子；(2) 海拔2 900 m以上、坡度大于30°且半陰坡、陰坡地帶是較不適宜及不適宜的分布區(qū)域，這與山高坡陡、熱量不足等局地植被生境限制條件有關；(3) 農(nóng)村居住地適宜性由東南部向西北部逐漸降低，適宜區(qū)占總面積的16.79%，次適宜區(qū)占26.93%，較不適宜區(qū)占32.47%，不適宜區(qū)占23.81%。研究結果可為地震災區(qū)安全選址與生境修復提供科技支撐。

雅安4·20地震; 農(nóng)村居住地; GIS; 適宜性評價

農(nóng)村居住地環(huán)境的形成和發(fā)展，是長期以來人類利用自然資源、適應自然環(huán)境的產(chǎn)物。山區(qū)地形、氣候及突發(fā)的地質(zhì)災害活動嚴重影響到居民的生存環(huán)境，易造成居住地適宜性發(fā)生變化。雅安“4·20”地震后，寶興河流域居民住房倒損約215萬m2，農(nóng)村房屋倒損約14.4萬間，各類基礎設施和公共服務系統(tǒng)受到嚴重損毀[1]，地震引發(fā)的環(huán)境破壞已被列為是導致農(nóng)村居住區(qū)環(huán)境變化的主要原因，受損生境的修復已成為當前地理學和生態(tài)學研究的熱點課題[2-4]。伴隨著3S技術、計算機技術的出現(xiàn)和發(fā)展，特別是GIS技術的應用，在大尺度范圍上開展適宜性研究已成為新的研究熱點[5-11]。然而傳統(tǒng)的適宜性評價往往依靠專家的野外觀察選取幾種生境因子，通過簡單的地圖疊加功能產(chǎn)生適宜性分布圖，數(shù)據(jù)來源及質(zhì)量難以保證、專業(yè)模型缺乏通用性、系統(tǒng)集成性差[12-14]。

雅安地震災區(qū)的農(nóng)村居住地環(huán)境適宜性評價是保證居民安居樂業(yè)和實現(xiàn)災后山區(qū)安全選址的基本前提。對災區(qū)居住地環(huán)境進行適宜性評價，不僅在于反映其地質(zhì)條件及自然環(huán)境狀況對生物的生存、繁衍等生命活動的影響及適宜狀況，更重要的是按適宜性級別制定出相應的災后生態(tài)恢復措施，因地制宜地進行生態(tài)重建和恢復。當前，災后重建研究多集中于地震引發(fā)的次生災害的評估與預測，對災區(qū)地質(zhì)條件及自然本底狀況的分析研究涉及較少。因此，依據(jù)國內(nèi)外適宜性評價的最新研究成果，綜合采用層次分析法、熵值法并結合GIS空間分析技術，選取寶興河流域農(nóng)村居住地為研究對象，針對農(nóng)村居住地環(huán)境特征開展定量評價，明確適宜性等級，研究結果可為地震災區(qū)重建規(guī)劃與場地選址提供科學依據(jù)。

1　研究區(qū)概況

寶興河流域位于雅安市北部的青衣江上游，寶興河系其主源，發(fā)源于夾金山東段巴朗山南麓螞蟥溝，上游分為東、西兩河，東河為主流，兩河在寶興縣城上游2 km處的河口相匯。本文的研究區(qū)域為雅安市寶興縣，全縣位于30°09′—30°56′N，102°29′—103°01′E，土地面積3 114 km2。以高山峽谷為主要地貌特征，地勢西北高東南低，海拔為700～5 300 m，受山地海拔影響，具有從亞熱帶到永凍帶的立體氣候，年均降水量939.15 mm，年均氣溫14.6℃[2]。同時該區(qū)植被類型多樣，垂直分布明顯：(1) 海拔1 800 m以下的基帶為偏濕性亞熱帶常綠闊葉林；(2) 海拔1 800～2 200 m為常綠和落葉闊葉混交林；(3) 海拔2 200～2 900 m為針闊混交林；(4) 海拔2 900～3 200 m為亞高山針葉林；(5) 海拔3 200 m以上為高山灌叢、草甸帶；(6) 海拔4 500 m以上為高山稀疏植被[15]。

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

2.1數(shù)據(jù)來源及預處理

2.1.1數(shù)據(jù)來源研究采用的數(shù)據(jù)主要源于：(1) 遙感數(shù)據(jù)來源于中科院對地觀測與數(shù)字地球科學中心的SPOT-5影像和Google公司的Google Earth在線數(shù)據(jù)；(2) 聚落高程、坡度、坡向數(shù)據(jù)來源于30 m精度的ASTER GDEM數(shù)據(jù)和90 m精度的ARTM數(shù)據(jù)；(3) 統(tǒng)計資料包括縣統(tǒng)計年鑒、縣志等；(4) 基礎地理信息數(shù)據(jù)包括寶興縣行政區(qū)劃圖、大熊貓棲息地世界遺產(chǎn)分布圖(1∶1 000 000)、寶興縣土地利用現(xiàn)狀圖、寶興縣居住地分布圖和中國植被類型圖(1∶1 000 000)；(5) 氣象數(shù)據(jù)來源于四川省氣象信息中心觀測站1999—2013年的統(tǒng)計資料；(6) 地質(zhì)災害數(shù)據(jù)來自于雅安市地質(zhì)災害分布圖(1∶100 000)。

2.1.2數(shù)據(jù)預處理(1) 基于ArcGIS 9.3軟件平臺，利用Arc Toolbox對DEM數(shù)據(jù)進行幾何校正及數(shù)據(jù)中雜點的消除，選擇WGS_1984_UTM_Zone_48N坐標系統(tǒng)，以實地采取的控制點為地理參考，獲取海拔、坡度、坡向等數(shù)據(jù)；(2) 利用植被類型圖、大熊貓棲息地世界遺產(chǎn)分布圖、地質(zhì)災害分布圖、居住地分布圖等，在ArcGIS軟件中利用Editor工具進行地圖數(shù)字化，完成植被、地質(zhì)、農(nóng)村居住地及大熊貓棲息地世界遺產(chǎn)等專題圖的制作；(3) 搜集、整理寶興河流域降水、溫度、日照等氣象資料，考慮自然環(huán)境如海拔、坡度等對氣溫、降水、日照等氣候因素的影響，利用空間插值法生成。

2.2研究方法

居住地環(huán)境適宜性反映了地質(zhì)條件及自然環(huán)境狀況對人類生存活動的適宜程度。對農(nóng)村居住地環(huán)境的適宜性評價采用多目標線性加權函數(shù)評價模型，充分考慮多種因素對農(nóng)村居住地環(huán)境適宜性評價的影響，計算公式如下：

(1)

式中：S——適宜性綜合評價指數(shù)；Wi——第i個評價因子的相對權重；Si——第i個評價因子的分值；n——評價因子數(shù)。

3　適宜性評價

3.1評價指標及權重

適宜性評價指標體系的選擇與量化是進行居住地環(huán)境評價的基礎，其合理與否直接關系到評價結果的科學性和完整性。因此，在指標選取過程中，以地形、植被、地質(zhì)災害及大熊貓棲息地世界遺產(chǎn)分布狀況為主要評定條件，同時考慮氣象因素的影響。遵循主導性、穩(wěn)定性及可操作性的原則，選擇了海拔、坡度、坡向3個地形因子；年均降水量、干濕季日照時數(shù)、干濕季均溫、干濕季相對濕度、干濕季蒸發(fā)量等9個氣象因子；植被類型、地層巖性、地震烈度、大熊貓棲息地世界遺產(chǎn)等共16個評定指標，考慮區(qū)域不同生態(tài)因子對適宜性影響程度不同，研究中采用問卷調(diào)查和專家打分法[16-19]，并運用熵值法修正的層次分析法確定不同生態(tài)因子對適宜性影響的權重及分級(表1)，其中，分值越高，表示適宜性程度越大。

表1　寶興河流域農(nóng)村居住地環(huán)境適宜性評價因子權重及分級

注：A1海拔(m)，A2坡度(°)，A3坡向，A4植被，A5地層巖性，A6大熊貓棲息地世界遺產(chǎn)，A7地震烈度，A8年均降水量(mm)，A9干季日照時數(shù)(h)，A10濕季日照時數(shù)(h)，A11干季平均溫度(℃)，A12濕季平均溫度(℃)，A13干季相對濕度(%)，A14濕季相對濕度(%)，A15干季蒸發(fā)量(mm)，A16濕季蒸發(fā)量(mm)。

3.2評價因子分級

在ArcGIS軟件中，以Spatial Analyst為分析工具，結合寶興河流域的地形和通過空間插值法生成的氣溫、降水、日照等氣象資料，完成地形和氣象因子分級圖的制作；以Editor為分析工具，完成植被、地層巖性、地震烈度、大熊貓棲息地世界遺產(chǎn)等評價因子分級圖的制作(附圖14)。

3.3結果與分析

根據(jù)地形、氣象、植被、地層巖性、地震烈度及大熊貓棲息地世界遺產(chǎn)等因子的分級標準(表1)，借助ArcGIS 9.3軟件平臺，應用多目標線性加權函數(shù)評價模型，在統(tǒng)一的空間框架中，對寶興河流域各評價因子的柵格數(shù)據(jù)進行空間疊加，計算出寶興河流域農(nóng)村居住地環(huán)境適宜性指數(shù)(0.740～3.804)，完成對研究區(qū)適宜性的綜合評價(圖1)。為了更加直觀地表示評價結果，根據(jù)適宜性指數(shù)值的大小，采用自然斷點法、按照數(shù)值由大到小的順序?qū)⒕C合評價指數(shù)劃分為4級，其代表的適宜性等級分別為適宜、次適宜、較不適宜、不適宜(表2)。

農(nóng)村居住地環(huán)境適宜性空間分布特征反映了地質(zhì)條件和自然因子對人類生產(chǎn)、生活產(chǎn)生的影響。通過對寶興河流域農(nóng)村居住地環(huán)境適宜數(shù)據(jù)的量化處理，得出該區(qū)適宜性面積及適宜性空間分布特征。將寶興河流域農(nóng)村居住地環(huán)境適宜性指數(shù)圖與人口分布圖疊加，兩者輪廓基本吻合，為研究結果的科學性提供了證明(圖1—2)。

3.3.1適宜性面積通過對寶興河流域16個評價因子的定量研究，結果表明：寶興河流域農(nóng)村居住地適宜性面積522.85 km2，比例最小，僅為16.79%；次適宜面積838.50 km2，所占比例較大，達26.93%；較不適宜面積最大，達1 011.13 km2，占總面積的32.47%；不適宜面積741.52 km2，占總面積的23.81%(圖3)。

3.3.2適宜性空間分布特征(1) 農(nóng)村居住地環(huán)境適宜性空間分布與地質(zhì)條件密切相關。農(nóng)村居住地環(huán)境適宜性分布呈現(xiàn)出明顯的海拔、坡度、坡向指向，其中具有較大影響作用的因素依次為海拔、坡度、坡向、地層巖性及地震烈度，氣象因子間差異并不顯著。適宜性區(qū)域主要集中在南部、東南部，且集聚特征明顯。海拔2 900 m以下、坡度小于30°、陽坡及半陽坡是適宜性和次適宜性的主要分布區(qū)。該區(qū)域海拔較低、坡度較小、朝南或盡可能朝南的坡向日照充足，耕地和林地資源豐富，更適宜居住。

(2) 農(nóng)村居住地環(huán)境適宜性空間分布具有沿河流及兩側(cè)集聚分布的特征。寶興河流域適宜性及次適宜性區(qū)域大部分沿河流及兩側(cè)分布，該區(qū)主要位于世界遺產(chǎn)保護區(qū)及外圍，生態(tài)服務功能明顯；海拔適宜，土壤穩(wěn)定性較高；氣溫、降水、日照等氣候條件適宜，有利于植物生長，植被類型多為喬木，生境質(zhì)量較好。

(3) 不適宜區(qū)域主要分布在高山區(qū)。在寶興河流域的中山區(qū)和高山區(qū)，適宜性指數(shù)偏低，適宜性等級表現(xiàn)為較不適宜和不適宜，生境修復難度較大，主要與該區(qū)的地形、地質(zhì)災害及氣候條件有關。該區(qū)海拔較高，地質(zhì)構造復雜，災害頻發(fā)，人類活動較少。地震、滑坡、泥石流等地質(zhì)災害及氣象因子成為該區(qū)環(huán)境質(zhì)量的主要限制因子。植被類型以高山灌叢、草甸為主，生境質(zhì)量差，生態(tài)風險較高。

表2　寶興河流域農(nóng)村居住地環(huán)境適宜性等級劃分標準

圖1　寶興河流域農(nóng)村居住地環(huán)境適宜性等級

圖2　寶興河流域農(nóng)村居住地分布

圖3　寶興河流域農(nóng)村居住地環(huán)境適宜性面積

4　結論與建議

地震災區(qū)生境修復是實現(xiàn)災區(qū)社會穩(wěn)定、居民安居樂業(yè)，加快山區(qū)新型城鎮(zhèn)化建設的重要內(nèi)容。針對寶興河流域農(nóng)村居住區(qū)的適宜性特征選取了16個評價因子，應用多目標線性加權函數(shù)評價模型，得出寶興河流域農(nóng)村居住地的環(huán)境適宜性空間分布特征，結論如下：寶興河流域農(nóng)村居住地環(huán)境適宜性程度在空間上呈現(xiàn)出東南部適宜程度較高，西北部山區(qū)適宜程度偏低的分布特征。適宜性及次適宜性區(qū)域集中分布于海拔2 900 m以下、坡度小于30°、陽坡及半陽坡地帶且大部分沿河流及兩側(cè)分布，具有沿河流集聚分布的特征；而中山區(qū)和高山區(qū)，適宜性指數(shù)偏低，適宜性等級表現(xiàn)為較不適宜和不適宜，生境修復難度較大，反映了地形和地質(zhì)災害等是決定寶興河流域農(nóng)村居住地環(huán)境適宜程度差異的主導因素。較不適宜、不適宜面積總和占總面積的56.28%，說明該區(qū)環(huán)境保護的急迫性。因此，在地震災區(qū)重建規(guī)劃過程中，寶興河流域災區(qū)生境修復可分區(qū)進行。

(1) 適宜性區(qū)域。該區(qū)域主要分布在寶興河一級支流兩側(cè)，地勢平坦，自然生態(tài)環(huán)境適宜，人口密集，是區(qū)域水電開發(fā)、生態(tài)旅游的集中地，同時也是世界遺產(chǎn)外圍保護區(qū)。該區(qū)可在原基礎上進行修復，充分利用現(xiàn)有資源，合理規(guī)劃用地類型。

(2) 次適宜性區(qū)域。該區(qū)域多分布在寶興河二級支流流域，大部分位于世界遺產(chǎn)保護區(qū)內(nèi)，地理位置和經(jīng)濟發(fā)展較好，具有改造和發(fā)展的潛力。應根據(jù)發(fā)展現(xiàn)狀改造交通條件，恢復植被，從而引導和促進該區(qū)域向適宜性區(qū)域的轉(zhuǎn)變。

(3) 較不適宜性區(qū)域。該區(qū)域主要分布在海拔較高的中山區(qū)，生態(tài)環(huán)境的敏感性強，經(jīng)濟發(fā)展水平較低，地形起伏大，難以滿足用地要求，且生境修復難度較大。該區(qū)域應加強生態(tài)環(huán)境保護，整治水土流失，提高土地的生產(chǎn)效益，進一步優(yōu)化空間格局。

(4) 不適宜性區(qū)域。該區(qū)域主要集中在高山區(qū)的地震斷裂活動帶和泄洪通道，易引發(fā)崩塌、滑坡、泥石流等次生地質(zhì)災害，不宜建設居民住房和永久性設施。該類區(qū)域的居民應遷至交通條件好、有發(fā)展空間的區(qū)位。

[1]王海茳,郭亞琳,王青.雅安地震災區(qū)寶興縣聚落空間特征分析[J].寧夏大學學報:自然科學版,2015,36(2):191-196.

[2]張建強,蘇鳳環(huán),范建容.“4·20”蘆山地震崩塌滑坡與公路危險性評價:以寶興縣省道S210沿線為例[J].山地學報,2013(5):616-623.

[3]方炫,劉德林,李壁成,等.黃土高原鄉(xiāng)級尺度土地利用適宜性評價:以寧南山區(qū)河川鄉(xiāng)為例[J].水土保持研究,2014,21(5):174-178.

[4]王強,袁興中,劉紅,等.山地河流生境快速評價模型與應用[J].水利學報,2011,42(8):928-933.

[5]Store R, Jokim?ki J. A GIS-based multi-scale approach to habitat suitability modeling[J]. Ecological Modelling,2003,169(1):1-15.

[6]Store R, Kangas J. Integrating spatial multi-criteria evaluation and expert knowledge for GIS-based habitat suitability modelling[J]. Landscape and Urban Planning,2001,55(2):79-93.

[7]蘇雷,朱京海,傅立群,等.GIS支持下的錦州葫蘆島沿海土地生態(tài)適宜性評價[J].水土保持研究,2013,20(4):207-212.

[8]張艷紅,何春光.基于GIS扎龍自然保護區(qū)丹頂鶴適宜生境動態(tài)變化[J].東北林業(yè)大學學報,2009,37(4):43-45.

[9]王筱明.基于生態(tài)適宜性的濟南市山地丘陵區(qū)荒地資源開發(fā)研究[J].水土保持研究,2013,20(4):264-267.

[10]王金亮,陳姚.3S技術在野生動物生境研究中的應用[J].地理與地理學,2004,20(6):44-47.

[11]任平,蘭亭超,周介銘,等.城鄉(xiāng)建設用地增減掛鉤區(qū)域適宜性評價與空間布局規(guī)劃研究:以成都龍泉驛區(qū)為例[J].水土保持研究,2014,21(1):272-282.

[12]蔡靖疆,辜寄蓉,王成善,等.岷江上游土地資源適宜性評價研究[J].礦物巖石,2001,21(1):91-97.

[13]肖燚,歐陽志云,朱春全,等.岷山地區(qū)大熊貓生境評價與保護對策研究[J].生態(tài)學報,2004,24(7):1373-1379.

[14]高彥凈,謝余初,錢大文,等.甘肅白龍江流域植被覆蓋度及景觀格局變化[J].水土保持研究,2015,22(1):181-187.

[15]汪麗娟,王青.寶興綜合生態(tài)系統(tǒng)管理功能分區(qū)與生態(tài)農(nóng)業(yè)建設[J].安徽農(nóng)業(yè)科學,2009,37(1):296-298.

[16]王介勇,劉彥隨,張富剛.海南島土地生態(tài)適宜性評價[J].山地學報,2007,25(3):290-294.

[17]胡彩婷,李巧云,關欣,等.永興縣發(fā)展冰糖橙產(chǎn)業(yè)的土地適宜性評價[J].水土保持研究,2012,19(6):228-236.

[18]王世波,崔金鑒,王世良.基于層次分析法的風險投資項目評價[J].工業(yè)技術經(jīng)濟,2003,22(6):87-89.

[19]王青,石敏球,郭亞琳,等.岷江上游山區(qū)聚落生態(tài)位垂直分異研究[J].地理學報,2013,68(11):1559-1567.

Suitability Evaluation of Rural Residence in Baoxing River Basin After Ya′an 4·20 Earthquake

MI Lindi, ZHOU Ying, WANG Qing

(SchoolofEnvironmentandResources,SouthwestUniversityofScienceandTechnology,Mianyang,Sichuan621010,China)

The mountainous area of Wudu-Luding line with fragile environment has become the most serious geographical unit after strong earthquake of Wenchuan and Lushan. It is intractable to reconstruct the damaged habitat after earthquake and inevitable to maintain the social stability of disaster areas. In order to evaluate residence environmental suitability, we selected topography, weather, vegetation, geological disaster, and world natural heritage of giant panda habitat as key factors in the rural residence of Baoxing River Basin. The results are as follows. (1) The distribution of rural residence suitability is characterized by altitude, slope, and aspect, and especially concentrated in the world heritage area and its periphery. This distribution characterization indicates that the terrain is the main influence factor of habitat suitability. The suitability and sub-suitability habitats mainly distribute in the sunny slope and half sunny slope positions where the altitudes are less than 2 900 m and slope gradients are less than 30°； (2) The half shady slope and shady slope places with altitude of greater than 2 900 m and slope greater than 30° are the main distribution areas of sub-unsuitability and unsuitable habitat. This distribution characterization is related to limited factors of vegetation growth, such as high altitude, steep slope and shortage of heat； (3) The suitability of rural residence reduced gradually from the southeast to the northwest. The suitable region accounts for 16.79% of the study area, and sub-suitable region, sub-unsuitable region and unsuitable region account for 26.93%, 32.47% and 23.81%, respectively. According to the above main evaluation results of residence environmental suitability, this research could provide scientific support for the selection of safe site and habitat restoration of earthquake areas.

Ya′an 4·20 earthquake; rural residence; GIS; suitability evaluation

2015-05-02

2015-06-01

科技部“十二五”國家科技支撐計劃(2015BAC05B04-01);四川省科技廳項目(2014SZ0063,2013SZ0101,2013SZ0159)

米琳迪(1989—),女,河北邢臺人,碩士研究生,研究方向為環(huán)境規(guī)劃與管理。E-mail:milindi1223@163.com

王青(1967—),男,山西陽高人,博士,教授,主要從事環(huán)境規(guī)劃與管理研究。E-mail:qingw@imde.ac.cn

F301; X826

A

1005-3409(2016)03-0199-05