摘要：基于GIS平臺(tái)，利用射洪縣的土地利用調(diào)查數(shù)據(jù)，并結(jié)合該縣的DEM、聚落、水域和道路數(shù)據(jù)，建立一、二級(jí)退耕還林決策模型。結(jié)果表明，射洪縣可以退耕還林的總面積為8 907.86 hm2，其中一級(jí)退耕還林的對(duì)象為坡度大于25°的旱地，退耕面積為6 314.23 hm2，占退耕總面積的70.9%，占旱地總面積的10.5%，占耕地總面積的7.8%；二級(jí)退耕還林的對(duì)象為坡度在15°～25°，且在聚落500 m和道路1 000 m范圍外與水域1 000 m范圍內(nèi)的耕地，退耕面積為2 593.63 hm2，占退耕總面積的29.1%，占耕地總面積的3.2%。本研究通過構(gòu)建退耕還林決策模型，可以快速準(zhǔn)確地獲取退耕還林?jǐn)?shù)據(jù)，為新一輪的退耕還林還草工程提供決策依據(jù)。

關(guān)鍵詞：退耕還林；GIS；DEM；空間疊加分析；決策模型

中圖分類號(hào)：P208 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2016）15-4048-04

DOI：10.14088/j.cnki.issn0439-8114.2016.15.062

Abstract： Based on GIS platform，this paper set up primary and secondary decision models of conversion from farmland to forests by using the data of land use survey and the data combination of DEM，settlements，water area and roads in Shehong county. The results show that the total area of giving farmland back to the forest is 8 907.86 hm2 in Shehong county，in which 6 314.23 hectare area can be regarded as the object of the primary model and 2 593.63 hectare area can be regarded as the object of the second model. The former area of which the slope is greater than 25 degree accounts for 70.9 percent of the reforesting area，10.5 percent of dry land area and 7.8 percent of total cultivated area. The latter area of which slope is between 15 degree and 25 degree accounts for 29.1 percent of the reforesting area，3.2 percent of the arable land area. The reforesting land area of the second model is 500 meter away from the villages and towns， and is beyond 1000 meter of the roads，and is within 1 000 meter of the rivers. In this study，through the construction of decision model of returning farmland into forest，people can quickly and accurately obtain the data in demand. Thus， it would provide some references for the implement of the new project of returning cultivated land into forest and grass.

Key words： conversion from farmland to forests； GIS； DEM； spatial overlay analysis； decision model

生態(tài)文明建設(shè)是社會(huì)、經(jīng)濟(jì)、資源可持續(xù)發(fā)展的必然選擇[1]， 良好的生態(tài)環(huán)境是人與社會(huì)走可持續(xù)發(fā)展道路的根基。四川省屬于多山地區(qū)，地形復(fù)雜、起伏大，耕地的坡度一般較大，進(jìn)而導(dǎo)致耕作難度偏大和耕作效率低下。又由于人們長期的盲目開墾，水土流失較為嚴(yán)重，這不僅給局部地區(qū)造成了生態(tài)安全威脅，也給生態(tài)文明建設(shè)提出了巨大挑戰(zhàn)[2-4]。鑒于坡耕地和陡坡耕地的水土保持能力較差，生產(chǎn)力較低，生產(chǎn)成本較高，中國存在大量可用于退耕還林還草的坡耕地和沙化耕地[5]。

傳統(tǒng)的地面調(diào)查和統(tǒng)計(jì)方式確定退耕還林耕地的面積及空間范圍分布比較耗費(fèi)人力、物力、財(cái)力和時(shí)間，而且受地面地形條件的限制和人為不確定性因素的影響，使得退耕耕地定位的精準(zhǔn)性不高。基于3S技術(shù)的研究應(yīng)用，科學(xué)退耕還林的生態(tài)建設(shè)方法已經(jīng)在研究和實(shí)踐中逐步應(yīng)用，如楊存建等[6，7]、楊鑫等[8]利用GIS和RS技術(shù)對(duì)四川省和云南省構(gòu)建了退耕還林還草決策模型，并定量地分析出適合退耕還林、還草的具體類型和范圍；湯國安等[9]、周丹等[10]、郭建英等[11]利用GIS技術(shù)對(duì)研究區(qū)的土地利用和坡度等級(jí)狀況進(jìn)行快速獲?。煌趵杳鞯萚12]利用GIS技術(shù)，設(shè)計(jì)了基于GRID的坡耕地退耕人口壓力指數(shù)模型，揭示了坡耕地與周圍耕地的關(guān)系。這些研究采用了不同的研究尺度和技術(shù)方法，但受數(shù)據(jù)來源和遙感影像分類精度的影響，對(duì)于退耕還林還草的精確性還有待于進(jìn)一步研究提高。本文利用GIS技術(shù)結(jié)合DEM和土地利用現(xiàn)狀數(shù)據(jù)，能夠快速、準(zhǔn)確地分析出精確到各級(jí)行政區(qū)的退耕耕地圖斑，實(shí)現(xiàn)退耕還林還草工程的科學(xué)規(guī)劃、高效實(shí)施和自動(dòng)化管理。

1 研究區(qū)概況

本研究區(qū)是地處四川盆地中部，涪江中游的遂寧市射洪縣。該區(qū)介于30°38′-31°10′N，105°10′-105°39′E之間，海拔299.0-674.4 m，西北高，東南低，全縣地貌類型以丘陵為主，丘陵占地區(qū)總面積的89.1%。射洪縣位于成渝經(jīng)濟(jì)區(qū)腹心，屬于全省“一極一軸一區(qū)塊”規(guī)劃中的成都都市圈增長極，其耕地面積有80 580.24 hm2，耕地包括旱地、水田和水澆地，以旱地為主，占耕地總面積的75%。

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

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

研究數(shù)據(jù)主要有射洪縣2013年土地利用變更調(diào)查數(shù)據(jù)，覆蓋全縣的1∶1萬地形圖數(shù)據(jù)，利用GIS軟件派生出格網(wǎng)大小為5 m×5 m的DEM和坡度數(shù)據(jù)，以及射洪縣的行政區(qū)劃界線矢量數(shù)據(jù)。

2.2 數(shù)據(jù)處理

將坡度分級(jí)數(shù)據(jù)與耕地?cái)?shù)據(jù)疊加分析得到耕地坡度分級(jí)數(shù)據(jù)，并在此基礎(chǔ)上結(jié)合射洪縣的耕地類型和聚落、水域與道路緩沖區(qū)數(shù)據(jù)，建立分級(jí)退耕決策模型，具體的研究技術(shù)路線如圖1所示。在ARCGIS 10.1軟件平臺(tái)中，采用DEM數(shù)據(jù)逐格網(wǎng)計(jì)算坡度生成坡度柵格數(shù)據(jù)，依據(jù)第二次全國土地調(diào)查耕地坡度分級(jí)標(biāo)準(zhǔn)（表1），利用工具箱3D Analyst Tools-Raster Reclass中的Reclassify工具，對(duì)坡度數(shù)據(jù)進(jìn)行重新分類，生成坡度分級(jí)圖，其中：一級(jí)平坡；二級(jí)緩坡；三級(jí)斜坡；四級(jí)陡坡；五級(jí)急、險(xiǎn)坡。

將坡度分級(jí)數(shù)據(jù)與耕地圖斑進(jìn)行疊加，為保證耕地圖斑的幾何和屬性的完整性、一致性，利用工具箱Spatial Analyst Tools-Zonal中的Zonal Statistics as Table工具統(tǒng)計(jì)每個(gè)耕地圖斑中各坡度級(jí)的權(quán)重，以權(quán)重值最高的坡度級(jí)作為該圖斑的坡度級(jí)。

2.3 不同坡度級(jí)耕地面積分析

統(tǒng)計(jì)各級(jí)坡度的耕地面積（表2）。射洪縣1～4級(jí)坡度耕地面積為74 192.18 hm2，占耕地總面積的92.08%，其中以3級(jí)坡度耕地為主，占耕地總面積的50.64%，5級(jí)坡度耕地占耕地總面積的7.92%。

對(duì)耕地中的旱地進(jìn)行分類統(tǒng)計(jì)分析（表3）。由表3可知，1～4級(jí)坡度旱地面積為54 058.07 hm2，占旱地總面積的89.54%，其中以3級(jí)坡度旱地為主，占旱地總面積的63.71%，5級(jí)坡度旱地占旱地總面積的10.46%。

3 退耕還林決策分析

本研究根據(jù)射洪縣退耕還林實(shí)際情況，采用分級(jí)退耕還林決策的方法，分為一級(jí)退耕還林的空間決策分析和二級(jí)退耕還林的空間決策分析。

3.1 一級(jí)退耕還林的空間決策分析

根據(jù)國家退耕還林還草相關(guān)政策規(guī)定，確定一級(jí)退耕還林的對(duì)象為坡度大于25°的耕地，由于25°以上的耕地基本為旱地，水田和水澆地較少，又考慮到旱地的耕地質(zhì)量、水土保持性、土壤侵蝕程度和生態(tài)環(huán)境效益等因素，本研究確定射洪縣一級(jí)退耕還林對(duì)象為坡度大于25°的旱地。在ARCGIS軟件中將坡度分級(jí)數(shù)據(jù)與耕地?cái)?shù)據(jù)進(jìn)行疊加分析，根據(jù)耕地坡度大于25°且耕地類型為旱地的退耕要求，建立一級(jí)退耕決策模型，得出退耕旱地面積為 6 314.23 hm2，占旱地總面積的10.5%。

3.2 二級(jí)退耕還林的空間決策分析

二級(jí)退耕還林是指在一級(jí)退耕還林完成退耕的基礎(chǔ)上，根據(jù)改善生態(tài)環(huán)境的需要可選擇退耕的對(duì)象[6]。二級(jí)退耕還林對(duì)象主要考慮15°<坡度≤25°的耕地?？紤]到勞動(dòng)力成本和耕地澆灌難度隨著耕地與居民地間距離的增加而增加，且距離河流較近的耕地土壤侵蝕嚴(yán)重，容易造成水土流失，河道泥沙堆積、堵塞，因此本研究結(jié)合射洪縣耕地與地形的實(shí)際情況，確定二級(jí)退耕還林對(duì)象為4級(jí)坡度耕地與該縣村莊和建制鎮(zhèn)、河流水庫面數(shù)據(jù)和道路線數(shù)據(jù)的緩沖區(qū)疊加分析得到的耕地。

根據(jù)該縣的二級(jí)退耕還林對(duì)象，將其二級(jí)退耕還林分為3種情況：①利用村莊、建制鎮(zhèn)數(shù)據(jù)在ARCGIS軟件中生成500 m緩沖區(qū)數(shù)據(jù)，與4級(jí)坡度耕地?cái)?shù)據(jù)進(jìn)行疊加分析，根據(jù)耕地坡度在15°～25°且距離村莊、建制鎮(zhèn)500 m范圍外的退耕要求，建立聚落退耕模型，得出退耕耕地面積為255.96 hm2；②利用河流、水庫面數(shù)據(jù)生成1 000 m緩沖區(qū)數(shù)據(jù)，與4級(jí)坡度耕地?cái)?shù)據(jù)進(jìn)行疊加分析，根據(jù)耕地坡度在15°～25°且距離河流、水庫1 000 m范圍內(nèi)的退耕要求，建立水域退耕模型，得出退耕耕地面積為3 413.38 hm2；③利用道路數(shù)據(jù)生成1 000 m緩沖區(qū)數(shù)據(jù)，與4級(jí)坡度耕地?cái)?shù)據(jù)進(jìn)行疊加分析，根據(jù)耕地坡度在15°～25°且距離道路1 000 m范圍外的退耕要求，建立道路退耕模型，得出退耕耕地面積為6 463.21 hm2。

基于聚落、水域和道路退耕模型得到的耕地地塊存在重疊的現(xiàn)象，需要對(duì)以上3種類型的疊加結(jié)果進(jìn)行復(fù)合疊加處理，在此基礎(chǔ)上，根據(jù)耕地坡度在15°～25°的條件下，且在村莊、建制鎮(zhèn)500 m和道路1 000 m范圍外與河流、水庫1 000 m范圍內(nèi)的退耕條件，建立二級(jí)退耕決策模型，得出退耕耕地面積為2 593.63 hm2，占耕地總面積的3.2%。再將一、二級(jí)退耕數(shù)據(jù)與射洪縣的行政區(qū)劃界線數(shù)據(jù)進(jìn)行疊加統(tǒng)計(jì)分析，即可得到各鄉(xiāng)鎮(zhèn)的一、二級(jí)退耕還林地塊和面積數(shù)據(jù)，統(tǒng)計(jì)結(jié)果見表4，空間分布見圖2。

4 結(jié)論與討論

1）射洪縣退耕還林總面積為8 907.86 hm2，其中一級(jí)退耕還林的對(duì)象為坡度大于25°的旱地，退耕面積為6 314.23 hm2，占旱地總面積的10.5%，占耕地總面積的7.8%，占退耕總面積的70.9%，主要分布于射洪縣的西北角；二級(jí)退耕還林的對(duì)象為坡度在15°～25°，且在村莊、建制鎮(zhèn)500 m和道路1 000 m范圍外與水域1 000 m范圍內(nèi)的耕地，退耕面積為2 593.63 hm2，占耕地總面積的3.2%，占退耕總面積的29.1%。

2）在射洪縣各鄉(xiāng)鎮(zhèn)中，金華鎮(zhèn)的一級(jí)退耕面積和退耕總面積最多，一級(jí)退耕面積為784.59 hm2，占全縣一級(jí)退耕總面積的12.4%，退耕總面積為872.56 hm2，占全縣退耕總面積的9.8%；伏河鄉(xiāng)的一級(jí)退耕面積和退耕總面積最少，一級(jí)退耕面積為35.19 hm2，占全縣一級(jí)退耕總面積的0.6%，退耕總面積為50.23 hm2，占全縣退耕總面積的0.6%；柳樹鎮(zhèn)的二級(jí)退耕面積最多，二級(jí)退耕面積為205.99 hm2，占全縣二級(jí)退耕總面積的7.9%；青堤鄉(xiāng)的二級(jí)退耕面積最少，二級(jí)退耕面積為11.09 hm2，占全縣二級(jí)退耕總面積的0.4%。

本研究利用已有的土地利用數(shù)據(jù)成果在GIS平臺(tái)可以快速、準(zhǔn)確地構(gòu)建科學(xué)、合理的退耕還林決策模型，得到精確到地形、地類、各行政級(jí)別以及地塊圖斑的退耕還林?jǐn)?shù)據(jù)，從而極大地提高退耕還林工程的工作、經(jīng)濟(jì)效率。由于本研究采用的分析數(shù)據(jù)較少，無法進(jìn)行具體有效的還林還草決策分析，可以結(jié)合研究區(qū)的土壤侵蝕數(shù)據(jù)、土壤類型數(shù)據(jù)、氣候數(shù)據(jù)、植被類型數(shù)據(jù)、林業(yè)經(jīng)濟(jì)數(shù)據(jù)等定量數(shù)據(jù)以及相關(guān)的林業(yè)規(guī)劃政策，建立精確到分析出各柵格單元的退耕還林還草決策模型，從而提出更加科學(xué)的退耕還林還草決策，也可以為其他地區(qū)的退耕還林還草工程提供參考依據(jù)。

參考文獻(xiàn)：

[1] 白 楊，黃宇馳，王 敏，等.我國生態(tài)文明建設(shè)及其評(píng)估體系研究進(jìn)展[J].生態(tài)學(xué)報(bào)，2011，31（20）：6295-6304.

[2] 江忠善，劉 志，賈志偉.地形因素與坡地水土流失關(guān)系的研究[J].中國科學(xué)院水利部西北水土保持研究所集刊（黃土高原試驗(yàn)區(qū)土壤侵蝕和綜合治理減沙效益研究專集），1990（12）：1-8.

[3] 李文華.長江洪水與生態(tài)建設(shè)[J].自然資源學(xué)報(bào)，1999，14（1）：1-8.

[4] 王效科，歐陽志云，肖 寒，等.中國水土流失敏感性分布規(guī)律及其區(qū)劃研究[J].生態(tài)學(xué)報(bào)，2001，21（1）：14-19.

[5] 李世東.中外退耕還林還草之比較及其啟示[J].世界林業(yè)研究，2002（2）：22-27.

[6] 楊存建，劉紀(jì)遠(yuǎn)，張?jiān)鱿?，?遙感和GIS支持下的云南省退耕還林還草決策分析[J].地理學(xué)報(bào)，2001，56（2）：181-188.

[7] 楊存建，馮 涼，徐育建，等.基于ARCGIS的四川省分級(jí)退耕還林還草空間決策分析[J].自然資源學(xué)報(bào)，2007，22（6）：986-993.

[8] 楊 鑫，楊武年，李雅輝，等.“3S”技術(shù)在達(dá)縣退耕還林（草）動(dòng)態(tài)監(jiān)測(cè)中的應(yīng)用[J].中國水土保持，2010（12）：68-70.

[9] 湯國安，楊瑋瑩，秦鴻儒，等.GIS技術(shù)在黃土高原退耕還林草工程中的應(yīng)用[J].水土保持通報(bào)，2002，22（5）：46-50.

[10] 周 丹，楊武年.基于3S技術(shù)的隆昌縣土地利用/覆被動(dòng)態(tài)監(jiān)測(cè)[J].地理空間信息，2014（5）：72-74.

[11] 郭建英，李錦榮，劉鐵軍，等.退耕還林工程建設(shè)對(duì)吳起縣土地利用/覆被變化及其土壤侵蝕的影響[J].水土保持研究，2013， 20（5）：1-6.

[12] 王黎明，楊燕風(fēng)，歐陽宏斌.基于GRID的坡耕地退耕壓力模型研究[J].地理研究，2001，20（4）：389-396.