基于生態(tài)安全的青龍滿族自治縣土地利用沖突識別

蔡天雨,徐 磊,陳亞恒,馬立軍,郭鳳玉,陳召亞

(1.河北農(nóng)業(yè)大學(xué) 國土資源學(xué)院,河北 保定 071000;2.河北農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院,河北 保定 071000;3.河北農(nóng)業(yè)大學(xué) 城鄉(xiāng)建設(shè)學(xué)院,河北 保定 071000)

在快速城市化發(fā)展背景下，土地資源超強度開發(fā)、空間資源分配不均、土地利用方式與結(jié)構(gòu)不盡合理等問題日益嚴峻，土地利用沖突矛盾突出[1]。隨著十九大報告中“堅持人與自然和諧共生”理念的提出，生態(tài)安全問題逐漸引起各界的高度重視[2]。土地利用沖突是國內(nèi)外學(xué)者的研究熱點，研究內(nèi)容集中于概念辨析、產(chǎn)生機制、分類、診斷與管理等方面。國內(nèi)相關(guān)研究仍處于探索階段，基于不同視角，阮松濤[3]、周德[4]等對土地利用沖突概念進行了闡述;陳威[5]、賀艷華[6]等分別針對土地利用沖突發(fā)生領(lǐng)域、研究內(nèi)容的差異對其類型進行了劃定;土地利用沖突是眾多因素綜合作用的結(jié)果，李永清[7]，劉嘉[8]等認為社會政策、權(quán)益爭奪是引發(fā)沖突的主要原因;土地利用沖突診斷主要有定性與定量分析兩類方法，常用定性分析方法有參與式調(diào)查法[9]、邏輯框架法[10]、博弈論分析法[11]，常用定量分析方法有PSR模型[12]、統(tǒng)計回歸分析[13]、適宜性評價[14]，通過對沖突的診斷評價以提出相應(yīng)的管理對策?？v觀相關(guān)研究，雖已取得一定的研究成果，但僅有少量研究針對土地利用沖突識別體系展開，基于生態(tài)安全視角構(gòu)建整個區(qū)域土地利用沖突識別體系的研究更是少之又少，因此開展相關(guān)研究十分必要。

本研究以京津冀生態(tài)屏障區(qū)——青龍滿族自治縣為研究區(qū)，基于生態(tài)系統(tǒng)服務(wù)功能重要性評價結(jié)果選取生態(tài)源地，利用生態(tài)敏感性評價結(jié)果與土地利用類型等因子構(gòu)建綜合阻力面，利用最小累積阻力(MCR)模型及原理識別生態(tài)廊道與緩沖區(qū)，在考慮廊道200 m范圍內(nèi)區(qū)域保護重要性的基礎(chǔ)上，構(gòu)建生態(tài)安全格局。通過分析耕地、建設(shè)用地對不同等級生態(tài)安全用地的占用，識別2016年青龍縣土地利用沖突空間分布特征，有效提出預(yù)防對策，以促進土地資源的可持續(xù)利用以及生態(tài)保護與城鎮(zhèn)經(jīng)濟的協(xié)調(diào)發(fā)展，為研究區(qū)土地利用規(guī)劃與管理提供參考依據(jù)。

1 研究區(qū)概況

青龍滿族自治縣(以下簡稱青龍縣)隸屬河北省秦皇島市，地處東經(jīng)118°33′31″—119°36′30″，北緯40°04′40″—40°36′52″。全縣設(shè)11個鎮(zhèn)、14個鄉(xiāng)，縣域總面積約3 506.10 km2。縣域內(nèi)土壤類型豐富，地形地貌特征鮮明，地貌類型以山地、丘陵為主，地形整體呈現(xiàn)東部、南部偏低，西部、北部偏高的特點。

青龍縣資源豐富，2016年被納入國家重點生態(tài)功能區(qū)?？h域內(nèi)物種繁多，是國家野生動物資源研究開發(fā)和保護基地，設(shè)有都山、老嶺自然保護區(qū);水資源豐富，水庫眾多，是青龍縣及周邊地區(qū)的重要水源地;森林覆蓋率高，青龍縣于90年代后先后實施多項造林工程，是京津冀地區(qū)的重要生態(tài)屏障。但是青龍縣地處北方土石山區(qū)，且屬于河北省暴雨中心地區(qū)之一，晝夜溫差大，降水多，年降雨量高達914.6 mm，產(chǎn)生水土流失、地質(zhì)災(zāi)害等自然災(zāi)害的潛在威脅較大，生態(tài)敏感性較高。

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

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

研究所用數(shù)據(jù)主要包括:DEM數(shù)據(jù)，來源于地理空間數(shù)據(jù)云平臺的GDEM產(chǎn)品，空間分辨率為30 m，通過GIS空間分析工具提取坡度、地形起伏度等信息;NDVI數(shù)據(jù)，來源于美國NASA網(wǎng)站16 d合成的MODIS數(shù)據(jù)，空間分辨率為250 m，通過最大合值合成法(MVC)計算NDVI值與植被覆蓋度;氣象數(shù)據(jù)，來源于中國氣象數(shù)據(jù)網(wǎng)的中國地面氣候標準值數(shù)據(jù)集，選取2016年1—12月逐日觀測數(shù)據(jù)，利用GIS空間插值工具得到降雨侵蝕力因子空間分布數(shù)據(jù);土壤數(shù)據(jù)、地質(zhì)災(zāi)害數(shù)據(jù)，分別來源于1∶10萬青龍縣土壤類型圖、《青龍滿族自治縣地質(zhì)災(zāi)害防治規(guī)劃》，經(jīng)配準后矢量化，得到研究區(qū)土壤類型、地質(zhì)災(zāi)害點與地質(zhì)災(zāi)害易發(fā)等級空間分布數(shù)據(jù);土地利用數(shù)據(jù)，來源于青龍縣2016年土地利用變更數(shù)據(jù)庫，參考《土地利用現(xiàn)狀分類標準(GB/T 21010—2017)》，將地類重新歸并分為耕地、園地、林地、建設(shè)用地、水域、未利用地六大地類。本研究中所有數(shù)據(jù)均采用統(tǒng)一的地理坐標和投影坐標，空間分辨率為30 m×30 m。

2.2 研究方法

2.2.1 生態(tài)系統(tǒng)服務(wù)功能重要性評價與生態(tài)敏感性評價 生態(tài)系統(tǒng)服務(wù)功能重要性與生態(tài)敏感性評價是生態(tài)安全格局成功構(gòu)建的基礎(chǔ)。生態(tài)系統(tǒng)服務(wù)功能指生態(tài)系統(tǒng)為滿足人類生存所提供生態(tài)空間的能力，本研究從水源涵養(yǎng)、生物多樣性、水土保持三方面展開研究;生態(tài)敏感性指生態(tài)系統(tǒng)因人類活動所引發(fā)生態(tài)變化的敏感性大小，本研究從水土流失與地質(zhì)災(zāi)害兩方面展開研究[15]。

(1) 生態(tài)系統(tǒng)服務(wù)功能重要性評價。水源涵養(yǎng)功能指生態(tài)系統(tǒng)用于協(xié)調(diào)區(qū)域水資源平衡的重要生態(tài)系統(tǒng)服務(wù)功能，本研究選取植被覆蓋度、土地利用類型、坡度三項指標展開水源涵養(yǎng)功能重要性評價。首先，將三項指標參考白楊、初東等[16]研究的分類標準進行五級賦值;其次，利用幾何平均數(shù)模型計算水源涵養(yǎng)重要性指數(shù)，利用自然斷點法將評價結(jié)果分為5個等級，具體分級標準見表1。

表1 生態(tài)系統(tǒng)服務(wù)功能重要性評價指標體系及分級標準

生物多樣性保護功能指生態(tài)系統(tǒng)為物種生存與繁衍所提供場所與保護作用的能力。生物多樣性與生境質(zhì)量存在顯著的正相關(guān)性[17]，因此本研究以Invest模型計算的生境質(zhì)量大小代表生物多樣性保護功能的重要性，模型的相關(guān)參數(shù)參照Invest模型(3.2.0)使用說明書[18]與相關(guān)研究[19-20]綜合確定，具體參數(shù)見表2—3。

表2 威脅因子及其最大影響距離、權(quán)重及衰退類型

根據(jù)公式(1)計算得出研究區(qū)生境質(zhì)量大小，其范圍在0～1之間，將其等間隔分為5個等級。計算公式如下[12]:

(1)

式中:Qxj為生境類型j中柵格單元x的生境質(zhì)量;Hj表示生境類型j的生境適宜度;Dxj表示生境類型j中x柵格單元的生境退化度;z為歸一化指數(shù)，默認設(shè)置2.5;k為半飽和常數(shù)，通常為生境最大退化值的一半，取值0.5。

紹興城投在高速發(fā)展的同時，挑戰(zhàn)和困難也一直隨行.就目前的發(fā)展現(xiàn)狀來看，紹興城投在發(fā)展過程中面臨的最主要困難有以下幾點：第一，工程任務(wù)重，但缺少相應(yīng)的資金平衡方案，建設(shè)資金很大程度上需要城投融資解決.第二，資產(chǎn)規(guī)模較小，缺少優(yōu)質(zhì)資產(chǎn)，融資成本較高.第三，自身造血功能偏弱，經(jīng)營能力、盈利水平不高.第四，城市資源相對分散.城投承擔城市基礎(chǔ)設(shè)施的主要建設(shè)任務(wù)，應(yīng)該同時享有市區(qū)各類城市經(jīng)營性資源.而這也是國內(nèi)相當一部分城投公司或多或少存在的通病.

表3 生境適宜度及其對不同威脅因子的相對敏感程度

水土保持功能指生態(tài)系統(tǒng)保持區(qū)域水土平衡、減少因水蝕導(dǎo)致土壤侵蝕的重要生態(tài)系統(tǒng)服務(wù)功能。本研究利用修正通用土壤流失方程(RUSLE)計算土壤保持量的大小，根據(jù)土壤保持量的大小判斷水土保持功能的強弱，按照表1中的規(guī)則將水土保持量分為5個等級。計算公式如下:

Ac=R×K×LS×(1-C×P)

(2)

式中:Ac為土壤保持量[t/(hm2·a)];R為降雨侵蝕力因子[MJ·mm/(hm2·h·a)]，采用查良松等[21]計算巢湖流域降雨侵蝕力因子的經(jīng)驗公式計算;K為土壤可蝕性因子[t·hm2·h/(MJ·mm·hm2)]，根據(jù)研究區(qū)土壤類型圖與門明新等[22]學(xué)者的研究結(jié)果，進行對比賦值；R,K兩個因子需在其計算結(jié)果的基礎(chǔ)上分別乘以17.02，0.131 7以轉(zhuǎn)為國際制單位;LS為坡度坡長因子(無量綱)，坡度因子采用劉寶元[23]的坡度公式計算，坡長因子采用Wischmeier[24]提出的公式計算;C為植被與經(jīng)營管理因子(無量綱)，采用蔡崇法等[25]學(xué)者建立的模型計算;P為水土保持措施因子(無量綱)，根據(jù)相關(guān)研究結(jié)果[26]與研究區(qū)實際情況對六大地類進行賦值，未利用地:1，林地:0.8，園地:0.8，耕地:0.5，建設(shè)用地:0，水域:0。

(2) 生態(tài)敏感性評價。水土流失敏感性評價可以有效評估因人類活動引發(fā)水土流失的可能性大小，本研究以RUSLE模型為基礎(chǔ)，展開水土流失敏感性評價(各指標含義同上);通過評價地質(zhì)災(zāi)害敏感性以識別地質(zhì)災(zāi)害易發(fā)區(qū)，達到加強預(yù)防的目的，本研究選取距災(zāi)害點距離、災(zāi)害易發(fā)程度、地形起伏度、人類活動強度四項指標對青龍縣地質(zhì)災(zāi)害敏感性展開評價。首先，將各指標參照周銳、王新軍等[27]學(xué)者的研究及《生態(tài)保護紅線劃定指南》[28]進行分級賦值，評價因子及其分級標準見表4;其次，利用幾何平均數(shù)模型計算水土流失、地質(zhì)災(zāi)害敏感性指數(shù)，依據(jù)自然斷點法將敏感性指數(shù)分為5個等級。

表4 生態(tài)敏感性評價指標體系及分級標準

(3) 幾何平均數(shù)模型。針對單項生態(tài)系統(tǒng)服務(wù)功能重要性與生態(tài)敏感性，本研究將評價指標進行五級賦值，參照楊遠琴等[29]研究，利用幾何平均數(shù)模型通過計算其單項重要性與敏感性指數(shù)來判斷其生態(tài)安全水平，計算公式如下:

(3)

式中:Sj為評價單元生態(tài)重要性/敏感性指數(shù);j=1,2,3，即分別對應(yīng)上文中水源涵養(yǎng)、水土流失、地質(zhì)災(zāi)害三項評價因子;Cij為第j項評價因子第i項指標的生態(tài)重要性/敏感性分值;n為第j類評價因子對應(yīng)的評價指標個數(shù)。

2.2.2 生態(tài)安全格局構(gòu)建

(1) 生態(tài)源地的識別。生態(tài)源地是生態(tài)過程向外擴散的起點，利于維護區(qū)域生態(tài)安全[30]。本研究選取水源涵養(yǎng)、生物多樣性保護、水土保持三項生態(tài)系統(tǒng)服務(wù)功能均屬于極重要等級的地區(qū)作為研究區(qū)生態(tài)源地。具體利用GIS空間分析工具，依據(jù)“木桶原理”提取三項生態(tài)系統(tǒng)服務(wù)功能極重要區(qū)，通過合并臨近斑塊、刪除零星碎小斑塊等處理，取并集得到研究區(qū)生態(tài)源地。

(2) 綜合阻力面的構(gòu)建。生態(tài)系統(tǒng)中物質(zhì)和能量的流動以及物種的遷徙都需要通過克服阻力來實現(xiàn)[17]。結(jié)合研究區(qū)生態(tài)基底，本研究利用綜合生態(tài)敏感性評價結(jié)果構(gòu)建綜合阻力面，同時考慮自然與社會因素，最終選取7個阻力因子，通過加權(quán)求和的方法構(gòu)建綜合阻力面。其中，各阻力因子的權(quán)重依據(jù)特爾斐法確定，綜合生態(tài)敏感性通過對單項敏感性結(jié)果取最大值鑲嵌的方法得到。各阻力因子及其分級標準與權(quán)重見表5。

表5 生態(tài)阻力評價指標體系及分級標準、權(quán)重

(4)

式中:MCR為最小累積阻力值;Dij為從生態(tài)源點j到空間單元i的距離;Ri為空間單元i的阻力系數(shù)。

根據(jù)公式(4)，將綜合阻力面作為MCR模型的成本數(shù)據(jù)，利用GIS空間分析模塊中的費用距離工具，計算每個像元到成本面上最近源地的最小累積阻力成本距離。接著利用成本路徑工具生成潛在生態(tài)廊道，由于該工具生成的成本路徑是雙向的，為了去除冗余，需將各廊道的累積阻力值量化，即以數(shù)值的形式將各廊道的累積阻力值提取出來，去除相同、相似廊道，并保留低阻力廊道，最終得到研究區(qū)優(yōu)化后的生態(tài)廊道。

(4) 緩沖區(qū)的識別。緩沖區(qū)是生態(tài)源地外圍為物種生存提供保護與恢復(fù)功能的生態(tài)潛力區(qū)。依據(jù)最小累積阻力面，利用自然斷點法劃定除去生態(tài)源地的4個不同等級生態(tài)安全水平的緩沖區(qū)。緩沖區(qū)作為物種修復(fù)的調(diào)節(jié)地帶，需要執(zhí)行嚴格的保護政策。隨著源地向外擴散，最小累積阻力值逐漸增加，代表生態(tài)保護用地擴張的阻力值越來越大，生態(tài)敏感性逐漸變?nèi)?，可進行合理的開發(fā)建設(shè)。

(5) 綜合生態(tài)安全格局構(gòu)建。本研究構(gòu)建的生態(tài)安全格局組分包括:源地、廊道、緩沖區(qū)。緩沖區(qū)的等級代表了生態(tài)安全格局水平的高低情況，因此依據(jù)緩沖區(qū)的劃定結(jié)果綜合劃定研究區(qū)不同安全水平的生態(tài)安全格局。同時，考慮到土地利用沖突識別體系中生態(tài)廊道與生態(tài)源地的同等重要性，在構(gòu)建綜合生態(tài)安全格局時，創(chuàng)新性的將生態(tài)廊道一定寬度內(nèi)的區(qū)域納入保護范圍。參考相關(guān)研究成果[31-32]，并考慮青龍縣生態(tài)廊道的生物多樣性保護、大型生物遷移通道、隔離、綠帶等功能，最終選取200 m作為生態(tài)廊道的保護寬度，將生態(tài)廊道保護區(qū)列為低生態(tài)安全水平行列。綜合考慮緩沖區(qū)與生態(tài)廊道劃定結(jié)果的基礎(chǔ)上，最終劃定青龍縣低、較低、中等、較高、高5個安全水平的生態(tài)安全格局。

2.2.3 土地利用沖突識別與分類 本研究分析生態(tài)安全視角下土地利用沖突的空間分布特征。基于以上研究，定義生態(tài)安全視角下土地利用沖突為:多利益主體因農(nóng)業(yè)耕作、城鎮(zhèn)擴張等經(jīng)濟建設(shè)活動對生態(tài)安全空間的不合理占用，由此引發(fā)的土地資源競爭結(jié)果。具體體現(xiàn)為耕地、建設(shè)用地對生態(tài)安全用地的占用。構(gòu)建土地利用沖突識別體系:首先，根據(jù)引發(fā)土地利用沖突的人類經(jīng)濟活動主導(dǎo)用地類型，將土地利用沖突定義為耕地、建設(shè)用地、綜合土地利用沖突三大類型，耕地、建設(shè)用地沖突分別指農(nóng)業(yè)耕作、開發(fā)建設(shè)活動對生態(tài)安全用地的占用，綜合土地利用沖突是耕地與建設(shè)用地沖突疊加效應(yīng)的綜合體現(xiàn);其次，確定土地利用沖突類型，即利用GIS的空間分析功能，根據(jù)耕地、建設(shè)用地分別對低、較低、中等、較高、高安全水平空間格局的占用依次將土地利用沖突劃定為極嚴重、嚴重、較嚴重、一般、不沖突5個類型，具體劃分標準見表6。生態(tài)安全水平的等級越低，其生態(tài)重要性越高，越應(yīng)執(zhí)行嚴格的保護政策，禁止人類建設(shè)開發(fā)活動，因此產(chǎn)生土地利用沖突時造成的后果越嚴重。

表6 土地利用沖突分類體系

3 結(jié)果與分析

3.1 生態(tài)系統(tǒng)服務(wù)功能重要性評價與生態(tài)敏感性評價

研究區(qū)單項生態(tài)系統(tǒng)服務(wù)功能重要性、生態(tài)敏感性評價結(jié)果空間差異較大。由附圖4可知，研究區(qū)水源涵養(yǎng)、生物多樣性、水土保持功能極重要區(qū)占比分別為28.34%,18.49%,3.14%。三項生態(tài)系統(tǒng)服務(wù)功能極重要等級區(qū)都分布于植被覆蓋度高值區(qū)，由于植物的保水固土能力強，生境質(zhì)量較高。由附圖5可知，研究區(qū)生態(tài)敏感性較高，其中水土流失、地質(zhì)災(zāi)害敏感性中極敏感與高度敏感區(qū)總和占比依次為20.56%,12.91%。生態(tài)敏感性高值區(qū)多分布在河流水系附近與生態(tài)環(huán)境脆弱的丘陵地區(qū)。

3.2 生態(tài)安全格局構(gòu)建

基于生態(tài)系統(tǒng)服務(wù)功能重要性評價結(jié)果提取生態(tài)源地，見附圖6，共得到生態(tài)源地16處，面積達698.92 km2，占研究區(qū)總面積的19.93%。綜合生態(tài)敏感性評價結(jié)果構(gòu)建綜合阻力面，見圖1A，總體來看研究區(qū)整體阻力值偏高，高阻力值主要集中于人類活動強度大的經(jīng)濟發(fā)展密集區(qū)。圖1B為最小累積阻力成本距離表面，隨著源地向外擴張受到的最小累積阻力值逐漸變大，代表生物進行遷徙所需克服的阻力越來越大，生態(tài)重要性逐漸降低。

圖1 綜合阻力面與最小累積阻力成本距離表面

圖2A,B分別為潛在與優(yōu)化后生態(tài)廊道，由圖2A可知，兩兩源地間均存在遷徙廊道，過于繁雜。通過量化各廊道的累積阻力值，可以發(fā)現(xiàn)各源地間廊道的累積阻力值大小各異，值越大就說明該廊道連接的源地之間進行物種遷徙、能量流動等活動越困難。經(jīng)處理，優(yōu)化后的生態(tài)廊道總長度495.52 km，可以充分滿足物種在源地之間的遷徙與能量的流動，應(yīng)避免人類經(jīng)濟活動對生態(tài)廊道的占用。

圖2 潛在生態(tài)廊道與優(yōu)化后生態(tài)廊道

如圖3所示，源地外圍為生態(tài)緩沖區(qū)，隨著源地向外擴張，最小累積阻力值逐漸增大，其生態(tài)敏感性逐漸降低，可以進行合理的開發(fā)建設(shè)。依據(jù)緩沖區(qū)的劃定結(jié)果與生態(tài)廊道200 m范圍內(nèi)區(qū)域的保護重要性綜合劃定研究區(qū)生態(tài)安全格局。如圖4所示，低、較低、中等、較高、高生態(tài)安全水平區(qū)依次占研究區(qū)總面積22.46%,26.96%,22.20%,18.57%,9.81%。其中低生態(tài)安全水平區(qū)包括生態(tài)源地與生態(tài)廊道200 m范圍內(nèi)的保護區(qū)，是需要進行重點保護的地區(qū)，該類地區(qū)生態(tài)敏感性較高，人類活動會對其造成較大的破壞，因此需要執(zhí)行嚴格的保護政策;較低生態(tài)安全水平區(qū)即為生態(tài)緩沖區(qū)，可起到保護生態(tài)源地的作用;中等生態(tài)安全水平區(qū)作為源地向外擴張的低阻力區(qū)，是自然與社會系統(tǒng)的過渡帶，仍需控制人類的社會開發(fā)活動;較高生態(tài)安全水平區(qū)作為中阻力區(qū)，生態(tài)敏感性較低，可適當進行人類經(jīng)濟活動;高生態(tài)安全水平區(qū)作為高阻力區(qū)，生態(tài)水平穩(wěn)定，適宜開展人類經(jīng)濟建設(shè)活動。

圖3 緩沖區(qū)

圖4 生態(tài)安全格局

3.3 土地利用沖突識別與分類

由圖5A可知，2016年青龍縣耕地沖突空間分布較為分散，主要由于青龍縣境內(nèi)以山地、丘陵為主，耕地面積小且分布分散。不、一般、較嚴重、嚴重、極嚴重沖突類型的面積分別為38.35 km2,74.83 km2,90.65 km2,102.93 km2,23.34 km2。一般、不沖突區(qū)主要分布在研究區(qū)西部的肖營子鎮(zhèn)、婁丈子鎮(zhèn)、朱杖子鄉(xiāng);較嚴重和嚴重沖突區(qū)主要分布在研究區(qū)東部的祖山鎮(zhèn)、龍王廟鄉(xiāng)、木頭凳鎮(zhèn);極嚴重沖突區(qū)主要分布在研究區(qū)中東部的祖山鎮(zhèn)、龍王廟鄉(xiāng)、青龍鎮(zhèn)，主要由于該類地區(qū)為人口密集區(qū)，隨著城鎮(zhèn)的快速擴張，加大了對農(nóng)產(chǎn)品的需求，引發(fā)耕地對生態(tài)安全用地的侵占，但是因該類地區(qū)生態(tài)環(huán)境現(xiàn)狀更適于發(fā)展生態(tài)安全用地，所以耕地沖突矛盾突出。

由圖5B可知，2016年青龍縣建設(shè)用地沖突相對集中。不、一般、較嚴重、嚴重、極嚴重沖突類型的面積分別為16.03 km2,37.18 km2,39.07 km2,31.81 km2,4.65 km2。一般、不沖突區(qū)主要分布在研究區(qū)西部的青龍鎮(zhèn)、肖營子鎮(zhèn)、八道河鎮(zhèn);較嚴重與嚴重區(qū)沖突主要分布在研究區(qū)北部與東南部的青龍鎮(zhèn)、祖山鎮(zhèn)、馬圈子鎮(zhèn);極嚴重沖突區(qū)分布相對分散，主要分布在研究區(qū)北部和東部的青龍鎮(zhèn)、祖山鎮(zhèn)、大巫嵐鎮(zhèn)。由以上格局分布特征可以發(fā)現(xiàn)，建設(shè)用地沖突高等級與低等級沖突格局在諸如青龍鎮(zhèn)、祖山鎮(zhèn)等經(jīng)濟發(fā)達區(qū)出現(xiàn)了同時分布的特征。經(jīng)分析可知，出現(xiàn)以上特征的原因為該類地區(qū)以經(jīng)濟發(fā)展為主，建設(shè)用地需求高，生態(tài)敏感性小，適宜開展人類經(jīng)濟活動，因此建設(shè)用地沖突矛盾不突出。但是，隨著人口的增長與城鎮(zhèn)擴張速度的加快，建設(shè)用地對生態(tài)安全用地的不合理占用問題突出，人類經(jīng)濟建設(shè)活動超出了土地資源的正常負荷，生態(tài)安全用地與建設(shè)用地間的土地利用沖突矛盾日益嚴重。

由圖5C可知，2016年青龍縣綜合土地利用沖突問題較為嚴重，且耕地沖突問題更為突出。極嚴重、嚴重、較嚴重綜合土地利用沖突面積達292.45 km2，占研究區(qū)總面積的8.34%。其中，一般、不沖突區(qū)主要分布在研究區(qū)西部的肖營子鎮(zhèn)、婁杖子鎮(zhèn)、八道河鎮(zhèn)。嚴重與較嚴重沖突區(qū)主要分布在研究區(qū)北部與東南部的祖山鎮(zhèn)、木頭凳鎮(zhèn)、青龍鎮(zhèn);極嚴重沖突區(qū)主要分布在研究區(qū)東部與西北部的祖山鎮(zhèn)、龍王廟鄉(xiāng)、青龍鎮(zhèn)。綜合來看，青龍縣土地利用沖突高等級發(fā)生區(qū)主要分布于土地利用類型復(fù)雜、經(jīng)濟與人口發(fā)展密集的東、北部生態(tài)脆弱區(qū);土地利用沖突低等級發(fā)生區(qū)主要分布于土地利用類型簡單、經(jīng)濟欠發(fā)達的西部偏遠地區(qū)。

圖5 土地利用沖突類型

4 結(jié)論與討論

(1) 青龍縣作為京津冀生態(tài)屏障區(qū)生態(tài)地位十分重要。水源涵養(yǎng)功能與生物多樣性保護功能優(yōu)勢明顯，水土保持功能相對較弱，三項生態(tài)系統(tǒng)服務(wù)功能極重要區(qū)都分布于植被覆蓋度高值區(qū)，該類地區(qū)物種豐富，生境質(zhì)量高，需重點保護;生態(tài)敏感性整體偏高，水土流失、地質(zhì)災(zāi)害敏感性高敏區(qū)多分布于河流水系附近與生態(tài)環(huán)境脆弱的丘陵地區(qū)，需執(zhí)行嚴格的保護政策。

(2) 本研究識別青龍縣生態(tài)安全格局組分有:源地、廊道、緩沖區(qū)，并依據(jù)識別結(jié)果綜合劃定研究區(qū)生態(tài)安全格局。生態(tài)源地分布廣泛，占研究區(qū)總面積的19.93%;緩沖區(qū)作為生態(tài)源地的保護屏障，對維持區(qū)域生態(tài)穩(wěn)定有重要作用;優(yōu)化后生態(tài)廊道總長495.52 km，應(yīng)避免人類活動對其占用;青龍縣大部分區(qū)域?qū)儆诘?、較低生態(tài)安全水平，二者占比高達49.42%，屬于重要生態(tài)保護區(qū);較高、高生態(tài)安全水平占比28.38%，生態(tài)水平穩(wěn)定，可以進行合理的建設(shè)活動。

(3) 青龍縣土地利用沖突分布廣泛。耕地沖突矛盾相對突出，且分布比較分散，極嚴重沖突區(qū)主要分布于研究區(qū)中東部，針對耕地沖突高發(fā)區(qū)，應(yīng)通過推進土地整治、復(fù)墾、退耕還林還草等工程的實施提升耕地質(zhì)量，逐步形成合理的耕地空間布局;建設(shè)用地沖突分布相對集中，極嚴重沖突區(qū)主要分布在研究區(qū)北部和東部，針對建設(shè)用地沖突高發(fā)區(qū)，應(yīng)通過提高土地節(jié)約集約利用程度緩解用地矛盾;面對青龍縣土地利用沖突現(xiàn)狀，可以從調(diào)整土地利用結(jié)構(gòu)、促進產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級、控制人口數(shù)量等方面入手展開防治措施的實施，降低土地利用沖突發(fā)生等級，以加快青龍縣生態(tài)城鎮(zhèn)的建設(shè)。

本研究基于生態(tài)安全視角識別青龍縣土地利用沖突空間分布特征，識別體系不僅考慮了常規(guī)生態(tài)安全格局的劃定，同時創(chuàng)新性的考慮了生態(tài)廊道的重要性，拓寬了土地利用沖突識別方法。研究中生態(tài)安全格局的構(gòu)建方法國內(nèi)外學(xué)者均廣泛使用，具有科學(xué)性。但是，本研究基于生態(tài)安全視角分析土地利用沖突分布情況，未考慮非生態(tài)安全用地之間的土地利用沖突情況。在今后的研究中，應(yīng)逐步完善生態(tài)安全視角下的土地利用沖突識別體系。