基于MCR模型的城鎮(zhèn)生態(tài)安全格局構建和建設用地開發(fā)模式

邱 碩，王宇欣，王平智，林 聰

?

基于MCR模型的城鎮(zhèn)生態(tài)安全格局構建和建設用地開發(fā)模式

邱 碩1,2，王宇欣1※，王平智1，林 聰1

（1.中國農業(yè)大學水利與土木工程學院，北京 100083；2. 河北省香河縣住房和城鄉(xiāng)建設局，廊坊 065400）

近年來，在中國城市土地開發(fā)的頂層設計中越發(fā)強調生態(tài)空間的重要性，充分識別城市生態(tài)安全格局并探討建設用地空間擴展格局成為劃定生態(tài)保護紅線、確定城市開發(fā)邊界和“三規(guī)合一”空間管制等政策實施的迫切需求。該文以唐山市為例，基于“生態(tài)過程—格局”理論和GIS的空間相關性分析、距離分析和疊加分析等方法，從綜合水安全、地質災害預警、生物生境保護、水土保持、游憩安全5個方面，整合構建了底線安全格局、緩沖安全格局和最優(yōu)安全格局3種不同等級的綜合生態(tài)安全格局。該文在MCR模型算法引入不同等級“源”的權重系數，以綜合生態(tài)安全格局為阻力因子，并結合城鎮(zhèn)中心的吸引力、主要道路的吸引力和政策調控（區(qū)域開發(fā)重點）等城鎮(zhèn)用地擴展影響因素，得到在不同約束條件下的區(qū)域空間擴展阻力面，分別確定了唐山市在“生態(tài)安全約束型”“經濟增長主導型”和“生態(tài)與經濟并重”3種模式下城鎮(zhèn)空間可能的發(fā)展范圍，經過多方案比較，“生態(tài)與經濟并重”城鎮(zhèn)建設用地空間發(fā)展模式是未來城鎮(zhèn)健康有序發(fā)展的最優(yōu)模式，這種發(fā)展模式既有利于減少大規(guī)模、高強度的城鎮(zhèn)化和工業(yè)化過程所造成的城市發(fā)展和環(huán)境保護之間的用地沖突，也可為城市用地布局規(guī)劃和生態(tài)紅線劃定提供科學依據。

土地利用；遙感；生態(tài)安全格局；建設用地；擴展；最小累積阻力模型；唐山市

0 引 言

中國正處于城鎮(zhèn)化快速發(fā)展時期，社會經濟高速發(fā)展與產業(yè)結構調整步伐加快，建設用地空間擴展作為城鎮(zhèn)化土地景觀轉變最明顯的標志，城鎮(zhèn)間對區(qū)域空間的競爭與控制所引發(fā)的矛盾日益突出。蔓延式的城鎮(zhèn)發(fā)展對生態(tài)環(huán)境造成了一系列不良影響，包括動植物棲息地的破壞、生物多樣性的減少以及流域生態(tài)系統的衰退[1-3]。同時，由于圍繞城市建成區(qū)的周邊多為耕地、林地等農田景觀，城市擴張必然導致生態(tài)環(huán)境質量的下降，如何控制城市蔓延已經成為許多國家面臨的挑戰(zhàn)，必須從生態(tài)安全的高度予以重視[4]。

生態(tài)安全是對生態(tài)系統的完整性和健康水平的整體反映，也是指所處的生態(tài)環(huán)境條件與具備生態(tài)系統服務可以有效支撐經濟發(fā)展，當生態(tài)遭到破壞時，也能保障人民生活、健康不受其損害的狀態(tài)與能力。如何進行精準地、實際地優(yōu)化和改善區(qū)域生態(tài)安全狀況，區(qū)域生態(tài)安全格局優(yōu)化成為目前國內外研究的熱點[5-6]。在表征和量化城市擴張的驅動力方面，Couch等[7]針對城市混合景觀和可持續(xù)性問題，探討了歐洲城市、郊區(qū)、后郊區(qū)之間的空間城市二元論與擴張的原因，并對混合驅動力進行了比較與探索。Fallah等[8]借鑒城市群理論，對美國都市圈的經濟績效與各自的擴張程度之間的關系進行了實證研究。近年來，國內部分學者開始嘗試通過城市空間的擴展進程，研究社會發(fā)展、區(qū)域經濟和自然環(huán)境等對城市擴張的影響機制。張俊鳳等[9]通過長三角城市土地利用效率近20年的時空變化特征，探討了長三角地區(qū)不同級別的城市建成區(qū)擴張與經濟增長之間關系以及土地利用戰(zhàn)略。一些學者借助約束性ACO-CA模型[10]、RF-CA模型[11]、神經網絡分析[12]等，在城市空間形態(tài)的情景分析、城市空間擴展模擬及城市不同擴展模式的多方案比較方面開展了研究探索[13]。然而，這些分析大多局限于在歷史數據的基礎上，構建城鎮(zhèn)空間增長模型[14-15]，對城市生態(tài)空間需求、城市內生以及城市規(guī)劃政策調控等因素在建設用地空間擴展上的影響等問題卻較少涉及，缺乏對城市擴展機制的理解和相關理論的探討。以生態(tài)安全格局為技術方法進行中國城鎮(zhèn)建設用地空間擴展研究，將對指導中國城鎮(zhèn)化健康發(fā)展起到重要作用，為城鎮(zhèn)建設用地的空間發(fā)展提供決策支持。

本文利用高分辨率多源遙感影像提取土地利用信息，結合GIS空間相關性、距離分析等方法，以唐山市為例進行實證研究。通過綜合考慮城鎮(zhèn)發(fā)展的影響因素，利用改進的最小累積阻力模型（MCR）分析3種模式下城鎮(zhèn)建設用地空間擴展格局，探討城鎮(zhèn)發(fā)展的適宜生態(tài)安全格局，從而為城鎮(zhèn)建設用地的規(guī)劃布局和生態(tài)紅線劃定提供科學依據。

1 研究方法與數據處理

1.1 阻力面模型

最小累積阻力模型（minimal cumulative resistance，MCR）廣泛應用于物種保護和景觀格局分析等生態(tài)領域[16]，國內學者將其在城市空間增長、土地利用、旅游規(guī)劃等社會經濟領域開始進行一些嘗試，如模擬城市擴展的空間演變過程[17]、城市土地生態(tài)適宜性評價[18]、生態(tài)連通性[19]等。本文借助MCR 模型引入影響區(qū)域城鎮(zhèn)擴展的阻力和動力因素，模擬源向周圍空間擴張過程。根據不同擴展源的向外擴展能力大小[20]，引入相應的阻力權重系數S。阻力面模型可表達為

式中″為城鎮(zhèn)擴展源穿越所有景觀單元最小累積阻力值；D為從源到評價單元的空間距離；R表示景觀單元對城鎮(zhèn)空間擴展的阻力系數；是一個未知的正函數，反映空間中任一點的最小阻力與其到所有源的距離和阻力面特性的正相關關系；S表示源所屬等級的相對阻力因子，城鎮(zhèn)用地的等級越高，擴展能力越強，其相對阻力因子就越?。籱in表示某景觀單元對不同的源取累積阻力最小值。

1.2 數據處理

唐山市位于河北省東部（117°31′-119°19′E、38°55′-40°28′N之間），土地總面積約為14 286 km2。本文采用2014年唐山市土地利用現狀圖斑矢量數據以及縣級行政區(qū)劃圖、土壤圖、地質地貌圖等矢量數據；對2014年的遙感影像圖和DEM數據進行拼接、裁剪等，利用ENVI4.7軟件計算其植被覆蓋指數；基于研究區(qū)DEM數字高程模型，提取研究區(qū)高程、地形起伏度、河流積水量，并對其進行重分類，生成所需矢量圖件。

2 研究區(qū)生態(tài)安全格局構建

2.1 單一生態(tài)安全格局構建

2.1.1 綜合水安全格局

水資源安全格局：根據唐山河流等級、水庫等級以及水源地保護級別（一級保護、二級保護），確定源等級系數，分別為1=0.80，2=0.90，3=1.00，越小反映水源的重要程度越高，利用ARCGIS10.0的平臺距離分析工具，利用公式（1）構建各等級源阻力面，綜合疊加利用幾何間斷分級法進行重分類。

洪水淹沒分析：利用水文分析模塊模擬流向、流量和填洼分析，找出徑流停滯的匯水盆域；通過Snap Pour Point和Stream link 的計算，找得到出水點或匯合點，即控制水流運動的戰(zhàn)略點。洪水淹沒區(qū)[21]是利用研究區(qū)的降水數據和數字高程模型，利用GIS水文分析工具以“無源淹沒”方法分析低于海拔低于指定水位的易澇區(qū)域。根據歷年洪水資料及城市防洪標準，建立10、20、50 a一遇3種風險頻率的洪水淹沒范圍。表1為不同安全格局下唐山市10、20、50 a一遇洪水淹沒情況。將地表水源安全格局范圍與開放式的防洪和雨洪管理范圍進行等權疊加，最終得到不同安全格局下的唐山市綜合水安全格局見圖1a。

表1 10、20、50 a一遇洪水淹沒分析和淹沒面積計算

2.1.2 地質災害預防安全格局

本文對研究區(qū)的地震液化、活動斷裂引起的地震災害及崩塌-滑坡和泥石流等突發(fā)性地質災害問題做主要研究。在唐山地區(qū)，活動斷裂是引起地震根本原因，區(qū)域性活動斷裂帶多為歷史或現代的地震帶，工程活動場地根據活動斷裂的空間分布及地震峰值加速度值等因素采取避讓措施。1976 年唐山大地震時，唐山市平原區(qū)及遵化遷安盆地均出現了地震液化地質災害，考慮到地震液化災害敏感性會受到多種因子綜合影響，災害發(fā)育程度參考1976年唐山市的地震液化分區(qū)[22]；根據地質液化危害嚴重程度，分為5個等級。

在崩塌-滑坡、泥石流空間熱點分析中，運用ARCGIS10.0中的Getis-Ord General*功能可以較準確地預測出局域空間自相關[23]，計算Getis-Ord General*公式如下

式中是點的總數；X和X分別表示第個和第個點的觀測值；W()是空間權重矩陣。

期望值的計算公式如下

正態(tài)分布條件下的*統計值的計算公式如下：

式中Var(G*)是方差。<0.05顯著呈現集聚，在的絕對值小于1.96的情況下，(G*)值較高表明相似且較高屬性值的熱點聚集；較低表明相似且較低的屬性值的冷點聚集；接近0表明不存在顯著的空間關聯模式[24]。

本文選擇1 km×1 km格網重采樣模擬崩塌-滑坡、泥石流空間分布特征，值為14.87，似然值小于1%。結合崩塌－滑坡災害點和泥石流災害點分布情況根據敏感性等級劃分成5級，參考蘇泳嫻等[25]研究成果，確定地質災害的影響因子及其敏感性劃分標準（表2），采用層次分析法，確定各評價指標權重，多因子綜合評價法，構建唐山市地質災害綜合敏感性評價模型，計算每個柵格單元的敏感性值，得到地質災害預防安全格局分布圖1b。

圖1 唐山市各單一及綜合生態(tài)安全格局

表2 地質災害的影響因子及其敏感性劃分標準

注：*是隨機一致性比率，< 0.1表明判斷矩陣具有良好的一致性，判斷是合理的。唐山市地質災害綜合敏感性評價模型= 0.0074<0.1。

Note:*is the ratio of random consistency,< 0.1, it shows that the comparison matrix has good consistency and the result is reasonable. Comprehensive sensitivity evaluation model of geological hazards in Tangshan City= 0.007 4<0.1.

2.1.3 生物生境保護安全格局

本研究綜合考慮生境質量（自身特征屬性）和生境適宜性進行評價來識別源地。以生態(tài)系統類型作為評價指標，依據研究區(qū)的土地覆被類型和NDVI指數劃分生態(tài)系統類型（見表3），并對指示物種生境適宜性進行評價。通過以上指示物種生境適宜性評價確定高度適宜、中等適宜和比較適宜的生物生境源地等級分別為1=0.80，2=0.90，3=1.00。

表3 生物生境保護適宜性評價指標及分級標準

生物生境保護最小累計阻力可通過ARCGIS10.0中柵格分析模塊（費用距離分析cost-distance）建立的實現公式（1），以不同阻力系數的土地覆被類型為阻力因子（見表4）。隨著阻力值不均勻的增大，在阻力值與（一定柵格數量）面積的關系曲線上出現跳躍性的門檻和突變點，其可以確定阻力閾值，根據阻力面來判別安全格局圖1c。

表4 關鍵物種空間運動阻力因子（土地覆被類型）與阻力系數

2.1.4 水土保持安全格局

本文通過水土保持重要性來構建區(qū)域水土保持安全格局，綜合考慮土壤侵蝕敏感性程度及其可能造成的對下游河流和水資源的危害程度[26]，評價標準參照生態(tài)功能區(qū)劃技術規(guī)程，得出水土保持安全格局圖1d。

2.1.5 游憩安全格局

研究區(qū)內自然保護區(qū)、森林公園、風景名勝區(qū)和濕地公園等自然景觀及文化和自然遺產保護區(qū)等人文景觀要素游憩價值和適宜性極高，因此，本文將其一一并選作游憩的景觀源，不同土地覆被類型對人類游憩可達性阻力系數（見表5），基于最小累積阻力模型在ARCGIS10.0平臺中運用費用距離分析（cost-distance）功能建立游憩安全格局阻力面，找到突變點的阻力閾值，構建游憩安全格局見圖1e。

2.2 宏觀生態(tài)安全格局整合

本文將綜合水安全格局、地質災害預防安全格局、生物生境保護安全格局、水土保持安全格局和游憩安全格局5個單一的生態(tài)安全格局進行等權疊加，根據安全水平的不同，形成底線、緩沖、最優(yōu)3種質量生態(tài)安全格局見圖1f，以此指導城市空間擴展并成為城鎮(zhèn)空間形態(tài)的基底條件。

底線安全格局是區(qū)域生態(tài)安全的核心區(qū)，此格局應納入到城市的基本生態(tài)用地，劃入城市的禁止建設區(qū)和生態(tài)紅線區(qū)內，并實行最嚴格的保護。此區(qū)域面積為2 187.53km2，占總用地面積的15.31%，人均生態(tài)用地288.69m2，主要分布在遷西縣北部山區(qū)、南部濱海濕地和各大中型水庫周圍，該區(qū)域景觀主要以林地、草地及灘涂濕地為主，存在水土流失嚴重，水資源匱乏等生態(tài)環(huán)境問題。緩沖安全格局面積為7 154.93km2，占總用地面積的50.08%，是底線安全格局外圍的緩沖地帶，位于生態(tài)環(huán)境較好的城市周邊地區(qū)，該區(qū)用于維護最優(yōu)景觀格局的區(qū)域生態(tài)功能服務，在進行開發(fā)時這個區(qū)域被用來保護和恢復當地的生態(tài)系統。最優(yōu)安全格局面積為1 814.05km2，占總用地面積的12.70%，它達到了維護區(qū)域生態(tài)安全的較為優(yōu)化狀態(tài)，最大限度上保護關鍵生態(tài)過程和生態(tài)基礎設施的完整性，此范圍可根據具體情況進行有條件的開發(fā)。

表5 游憩過程阻力因子（土地覆被類型）與阻力系數

3 基于改進的MCR模型的城鎮(zhèn)建設用地空間擴展

3.1 城鎮(zhèn)空間擴展“源”的設定

以市域城鎮(zhèn)建設用地作為空間擴展“源”[20]，同時考慮各“源”的擴展能力如區(qū)位條件、經濟實力、人口狀況和開發(fā)程度等影響城市空間增長的內在因素以及外在的政策調控，并結合市域的城鎮(zhèn)等級結構劃分對“源”進行等級劃分，源等級系數越小反映城鎮(zhèn)空間擴張能力越強（表6）。

3.2 城鎮(zhèn)空間擴展阻力層構建

在市域城鎮(zhèn)空間擴展的影響因素選取方面，近年來，規(guī)劃政策調控和生態(tài)空間質量等因素對于空間擴展的影響作用逐步增大，為了比較準確地探測出空間發(fā)展重點聚集區(qū)域，主要結合研究區(qū)土地利用總體規(guī)劃、城市總體規(guī)劃等空間規(guī)劃的重點發(fā)展戰(zhàn)略，將規(guī)劃政策調控因素（區(qū)域開發(fā)重點）引入空間擴展影響因素中。根據各阻力因子的影響程度進行阻力賦值，采用層次分析法，確定出不同阻力因子相應權重（表7），其中唐山市城鎮(zhèn)建設用地擴展影響因子=0.0755<0.1，表明判斷矩陣具有良好的一致性。本文以生態(tài)安全格局為生態(tài)阻力面，作為城鎮(zhèn)空間發(fā)展控制的基本保障要素，又考慮到建設用地擴展主要受城鎮(zhèn)中心的吸引力、主要道路的吸引力等發(fā)展條件影響一同作為成本距離柵格（圖2）。利用ARCGIS10.0空間熱點分析理論選擇1km×1km網格重新采樣模擬重點開發(fā)區(qū)域空間分布特征見圖3，值為16.75，似然值小于1%。根據阻力系數分級標準，建立不同等級源的城鎮(zhèn)空間擴展綜合阻力面（圖4）。

表6 城鎮(zhèn)的用地源的分級及影響因素

表7 城鎮(zhèn)建設用地擴展的影響因子與阻力系數

a. 綜合生態(tài)安全格局

a.Comprehensive ecological security pattern

b. 與主要交通干道距離

b. Distance from main traffic road

c. 重點開發(fā)區(qū)域

c. Major development zone

d. 與城鎮(zhèn)中心距離

圖3 唐山重點開發(fā)區(qū)域空間熱點分析

圖4 唐山市城鎮(zhèn)空間擴展綜合最小累積阻力面

3.3 城鎮(zhèn)建設用地空間擴展模式

3.3.1 生態(tài)安全約束型的城鎮(zhèn)空間擴展模式

本文以綜合生態(tài)安全格局為最小累積阻力面確定城鎮(zhèn)空間擴展范圍，判別最小累積阻力值（像元值）與用地面積（柵格數量）的直方圖中阻力閾值，像元值分別在1.1、3.7、8.9和131.5處確定為突變點，劃分出不同生態(tài)安全水平約束下的城鎮(zhèn)建設用地增長范圍（圖5a）。減去該區(qū)已建成區(qū)、基本農田面積和水域面積，最優(yōu)生態(tài)安全水平下的城鎮(zhèn)擴展模式能提供的實際建設用地面積為208.21km2，占研究范圍總面積的1.46%，以人均城鎮(zhèn)建設用地110 m2/人粗略計算，可以滿足134萬人的城鎮(zhèn)建設用地規(guī)模；緩沖生態(tài)安全水平下的城鎮(zhèn)空間擴展模式可提供的建設用地和農業(yè)用地規(guī)模為3 382.23km2，占研究范圍總面積的23.68%。底線生態(tài)安全水平下的城鎮(zhèn)空間擴展可提供的建設用地和農業(yè)用地規(guī)模為4 626.78 km2，占研究范圍總面積的32.39%。從圖上可以看出，在遷西縣、遷安市等城鎮(zhèn)大部分地段位于北部丘陵區(qū)，大部分屬于生態(tài)嚴格控制區(qū)，生態(tài)服務功能顯著，阻力值較高，城鎮(zhèn)擴展生態(tài)安全性越低，越不適宜城鎮(zhèn)開發(fā)性建設。

圖5 唐山市城鎮(zhèn)空間擴展模式

3.3.2 經濟增長主導型的城鎮(zhèn)空間擴展模式

非飽和供應和適度供應的下土地供應總量能得到有效控制[27]，本文以城鎮(zhèn)建設用地的最飽和擴展量來探討經濟增長主導型建設用地擴展模式，根據《唐山市城市總體規(guī)劃（2010－2020年）》做出的人口發(fā)展預測，研究范圍常住城鎮(zhèn)人口在2020年和2030年將分別達到720萬人和900萬人，“經濟增長主導”型的發(fā)展模式下，主要考慮影響建設用地擴展的經濟因素，將以底線生態(tài)安全格局邊界作為生態(tài)景觀安全臨界狀態(tài)，同時要在維持規(guī)劃期內基本農田保護面積保持在4 822km2以上的前提下，剩余范圍為其可擴展范圍（見圖5b）。分別模擬2020年、2030年的建設用地的擴展趨勢，外圍可擴展實際面積上限為1 038.20km2，占研究范圍總面積的7.27%。由圖6分析發(fā)現，中部平原地區(qū)城鎮(zhèn)建設用地沿主要公路空間擴展范圍明顯，灤縣、灤南縣和玉田縣與中心城區(qū)之間的綠帶幾乎消失并聯通成為一個整體；北部各城鎮(zhèn)建設用地組團雖沒有完全連成片，空間特征表現為沿112國道至唐南公路城鎮(zhèn)發(fā)展軸的縱向發(fā)展較為明顯。阻力從低值到高值的過渡過程中有更多的生態(tài)用地被侵占，生態(tài)要素間的連通度和生態(tài)服務功能降低，區(qū)域生態(tài)安全受到嚴重威脅。

3.3.3 “生態(tài)與經濟”并重型的城鎮(zhèn)空間擴展模式

城鎮(zhèn)用地的空間擴展受到多種因素的制約，不僅僅要考慮區(qū)域生態(tài)環(huán)境因素，而且要特別關注區(qū)域社會經濟條件和政策導向因素，在“生態(tài)控制”條件約束和經濟活動驅動下的空間擴展格局（圖5c），低阻力區(qū)主要分布在建設用地外圍，減去該區(qū)的水域面積、基本農田面積和已建成的城鎮(zhèn)建設用地面積，可提供的建設用地規(guī)模為249.75km2，占區(qū)域總面積的1.75%；以人均城鎮(zhèn)建設用地110m2/人粗略計算，可以滿足227萬人的城鎮(zhèn)建設用地規(guī)模。該區(qū)城市擴展的累積阻力很小，未來城市化概率極大，京沈高速軸線低密度橫向蔓延嚴重，集中在中心城區(qū)、南部沿海地區(qū)和一些重點發(fā)展城鎮(zhèn)，南部沿海公路軸線和運河沿岸連綿式擴展趨勢明顯。中度阻力區(qū)指中等累積阻力擴展范圍，其可擴展實際面積為308.53km2，占區(qū)域總面積的2.16%，可以滿足280萬人的城鎮(zhèn)建設用地規(guī)模，此區(qū)域短期內城市化概率較小，主要分布在低阻力區(qū)域的外圍城郊地區(qū)鄉(xiāng)鎮(zhèn)，以城市和鄉(xiāng)村景觀為主，低密度擴展；其中，中心城區(qū)沿G205國道呈條帶狀發(fā)展、遠郊呈小組團狀分布，灤河外圍沿岸帶狀縱向擴張，但擴展壓力較大、擴展空間有限。高阻力區(qū)為擴展累積阻力較大，近期不可能進行城市擴展的地區(qū)，該區(qū)面積為4 227.02km2，占區(qū)域總面積29.59%，其可擴展實際面積為121.20km2，可以滿足110萬人的城鎮(zhèn)建設用地規(guī)模。極高阻力區(qū)是生態(tài)脆弱區(qū)和自然環(huán)境較好且對城鎮(zhèn)空間擴展具有很強約束性的地區(qū)，該區(qū)面積為3 135.43 km2，占區(qū)域面積的21.95%，該區(qū)域基本農田和生態(tài)用地范圍內應禁止發(fā)展建設用地。

4 結 論

1）本研究融合了地統計分析、空間相關性分析、距離分析和疊加分析等空間分析方法，通過綜合水安全、地質災害預防、生物生境安全、水土保持、游憩安全等5個單一生態(tài)安全格局構建了底線安全格局、緩沖安全格局和最優(yōu)安全格局等3種不同等級的綜合生態(tài)安全格局，面積分別占研究范圍總面積的15.31%、50.08%和12.70%。

2）城鎮(zhèn)用地的空間擴展受到多種因素的制約，不僅僅要考慮區(qū)域生態(tài)環(huán)境因素，而且要特別關注區(qū)域社會經濟條件和政策導向因素。通過在MCR模型算法引入不同等級“源”的權重系數，分別確定在“生態(tài)安全約束型”、“經濟增長主導型”和“生態(tài)與經濟并重”3種模式下城鎮(zhèn)空間可能的發(fā)展范圍。研究表明，“生態(tài)與經濟并重”下的城鎮(zhèn)建設用地空間發(fā)展模式能夠達到“生態(tài)—經濟—社會”效益的平衡，更符合城鎮(zhèn)建設用地的空間發(fā)展規(guī)律，能有效預測城市空間發(fā)展范圍，是未來城鎮(zhèn)健康有序發(fā)展的適宜模式，對于優(yōu)化土地利用布局和生態(tài)環(huán)境保護均具有重要意義。

3）在城鎮(zhèn)用地擴展源的設定方面，工礦用地的等級應考慮生態(tài)和經濟發(fā)展的特殊性，根據土地利用規(guī)劃的導向，具體判斷建設用地擴展源的等級。但在城鎮(zhèn)空間擴展影響因子選擇方面，可以增加科技水平、區(qū)域流通等指標進行綜合分析，同時，在城市擴展的內生機制方面還需進行深入研究。

[1] Jian Peng, Liu Yanxu, Li Tianyi, et al. Regional ecosystem health response to rural land use change: A case study in Lijang City, China[J]. Ecological Indicators, 2017, 72: 399－410.

[2] Wang Wangjing, Guo Huacheng, Chuai Xiaowei, et al. The impact of land use change on the temporospatial variations of ecosystems services value in China and an optimized land use solution[J]. Environmental Science & Policy, 2014, 44: 62－72.

[3] 孫萍，唐瑩，Robert J Mason，等. 國外城市蔓延控制及對我國的啟示[J]. 經濟地理，2011，31(5)：748－753. Sun Ping, Tang Ying, Robert J Mason, et al. Inspiration of urban sprawl control studies in foreign countries[J]. Economic Geography, 2011, 31(5): 748－753.（in Chinese with English abstract）

[4] 俞孔堅，王思思，李迪華，等. 北京市生態(tài)安全格局及城市增長預景[J]. 生態(tài)學報，2009，29(3)：1189－1204. Yu Kongjian, Wang Sisi, Li Dihua, et al. The function of ecological security patterns as an urban growth framework in Beijing[J]. Acta Ecologica Sinica, 2009, 29(3): 1189－1204.（in Chinese with English abstract）

[5] 劉國華. 西南生態(tài)安全格局形成機制及演變機理[J]. 生態(tài)學報，2016，36(22)：7088－7091. Liu Guohua. Formation and evolution mechanism of ecological security pattern in Southwest China[J]. Acta Ecologica Sinica, 2016, 36(22): 7088－7091.（in Chinese with English abstract）

[6] 馬克明，傅伯杰，黎曉亞，等. 區(qū)域生態(tài)安全格局：概念與理論基礎[J]. 生態(tài)學報，2004，24(4)：761－768. Ma Keming, Fu Bojie, Li Xiaoya, et al. The regional pattern for ecological security （RPES）: The concept and theoretical basis[J]. Acta Ecologica Sinica, 2004, 24(4): 761―768.（in Chinese with English abstract）

[7] Couch C, Leontidou L, Petschel H G. Urban Sprawl in Europe:Landscapes, Land-Use Change & Policy[M]. Oxford:

Blackwell Publishing Ltd , 2007.

[8] Fallah B N, Partridge M D, Olfert M R. Urban sprawl and productivity, Evidence from US metropolitan areas[J]. Papers in Regional Science, 2011, 90(3): 451－472.

[9] 張俊鳳，劉友兆. 城市建成區(qū)擴張與經濟增長間的關系一以長三角地區(qū)為例[J]. 城市問題，2013，211(2)：11－16. Zhang Junfeng, Liu Youzhao . Relation between built-up area expansion and economic growth: Taking the Yangtze Delta Area for example[J]. Urban Problems, 2013, 211(2): 11－16. (in Chinese with English abstract)

[10] Liu Ming-hao, An Guangwen, Li Chao, et al. Urban land use simulation based on ACO-CA model embedded dynamic neighborhood theory: A case study of Shapingba District in Chongqing[J].Geography and Geo-Information Science, 2016, 32(3): 74－80.

[11] Kamusoko C, Gamba J. Simulating urban growth using a random forest-cellular automata (RF-CA) Model[J]. ISPRS International Journal of Geo-Information, 2015, 4(2): 447－470.

[12] 穆飛翔，蒲春玲，劉祥鑫. 烏魯木齊市轄區(qū)建成區(qū)規(guī)模擴張及預測分析[J]. 上海國土資源，2017，38(3)：30－35. Mu Feixiang, Pu Chunling, Liu Xiangxin. Scale expansion and forecast analysis of built-up area in Urumqi Municipal district[J]. Shanghai Land & Resources, 2017, 38(3): 30－35. (in Chinese with English abstract)

[13] 張榮天，張小林. 國內外城市空間擴展的研究進展及其述評[J]. 中國科技論壇，2012(8)：151－155. Zhang Rongtian，Zhang Xiaolin. Research progress and review of urban spatial expansion[J]. Forum on Science and Technology in China, 2012(8): 151－155.（in Chinese with English abstract）

[14] 崔福全，徐新良，孫希華. 上海城市空間擴展過程模擬預測的多模型對比[J]. 生態(tài)學雜志，2012，31(10)：2703－2708. Cui Fuquan, Xu Xinliang, Sun Xihua. Simulation and prediction of urban spatial expansion in Shanghai: A comparison of multiple models[J]. Acta Ecologica Sinica, 2012, 31(10): 2703－2708. (in Chinese with English abstract)

[15] 徐昔保，楊桂山，張建明. 近50年蘭州城市空間擴展模擬及其未來預測[J]. 干旱區(qū)研究，2009，26(5)：763－769. Xu Xibao, Yang Guishan, Zhang Jianming. Simulation and prediction of urban spatial expansion of Lanzhou city[J]. Arid Zone Research, 2009, 26(5): 763－769.（in Chinese with English abstract）

[16] 羅曉，李雙江. 石家莊市生態(tài)安全格局的識別與優(yōu)化[J]. 安徽農業(yè)科學，2011，39(27)：16730－16733. Luo Xiao, Li Shuangjiang. Recognition and optimization of ecological security pattern in Shijiazhuang[J]. Journal of Anhui Agricultural Sciences, 2011, 39(27): 16730－16733.（in Chinese with English abstract）

[17] 普鵑鵬，劉雪峰. 基于最小累計阻力模型的城市擴張空間布局優(yōu)化及其生態(tài)效應一以武漢市為例[J]. 現代商業(yè)，2017(6)：152－155.

[18] 劉孝富，舒儉民，張林波. 最小累積阻力模型在城市土地生態(tài)適宜性評價中的應用一以廈門為例[J]. 生態(tài)學報，2010，20(2)：421－428. Liu Xiaofu, Shu Jianmin, Zhang Linbo. Research on applying minimal cumulative resistance model in urban land ecological suitability assessment: As an example of Xiamen[J]. Acta Ecologica Sinica, 2010, 20(2): 421－428. (in Chinese with English abstract)

[19] 李謙，戴靚，朱青，等. 基于最小阻力模型的土地整治中生態(tài)連通性變化及其優(yōu)化研究[J]. 地理科學，2014，34(6)：733－739. Li Qian, Dai Jing, Zhu Qing, et al. Ecological connectivity changes and its pattern optimization during land consolidation based Minimal Accumulative Resistance model[J]. Scientia Geographica Sinica, 2014, 34(6): 733－739. (in Chinese with English abstract)

[20] 葉玉瑤，蘇泳嫻，張虹鷗，等. 生態(tài)阻力面模型構建及其在城市擴展模擬中的應用[J]. 地理學報，2014，69(4)：485－496. Ye Yuyao, Su Yongxian, Zhang Hong'ou, et al. Ecological resistance surface model and its application in urban expansion simulations[J]. Acta Geographica Sinica, 2014, 69(4): 485－496. (in Chinese with English abstract)

[21] 周銳，王新軍，蘇海龍，等. 平頂山新區(qū)生態(tài)用地的識別與安全格局構建[J]. 生態(tài)學報，2015，35(6)：2003－2012. Zhou Rui, Wang Xinjun, Su Hailong, et al. Identification and security pattern of ecological land in Pingdingshan newly developed area[J]. Acta Ecologica Sinica, 2015, 35(6): 2003－2012. (in Chinese with English abstract)

[22] 魏風華. 河北省唐山市地質災害風險區(qū)劃研究[D]. 北京：中國地質大學（北京），2006. Wei Fenghua. Researches on Geological Hazard and Risk Zonation in Tangshan[D]. Beijing: China University of Geosciences (Beijing), 2006. (in Chinese with English abstract)

[23] 羅文，羅畏. 基于空間統計的水質相關性分析[J].水電能源科學，2011，29(3)：27－30. Luo Wen, Luo Wei. Correlation analysis of water quality based on spatial statistics[J]. Water Resources and Power, 2011, 29(3): 27－30. (in Chinese with English abstract)

[24] 陳菁，吳麗芳，陳鑫彬. 城市景觀安全格局及其預案研究一以廈門市為例[J]. 中國農學通報，2014，30(2)：199－207. Chen Jing, Wu Lifang, Chen Xinbin. Research of the urban landscape security pattern and plan: As an example of Xiamen City[J]. Chinese Agricultural Science Bulletin, 2014, 30(2): 199－207. (in Chinese with English abstract)

[25] 蘇泳嫻，張虹鷗，陳修治，等. 佛山市高明區(qū)生態(tài)安全格局和建設用地擴展預案[J]. 生態(tài)學報，2013，33(5)：1524－1534. Su Yongxian, Zhang Hong’ou, Chen Xiuzhi, et al. The ecological security patterns and construction land expansion simulation in Gaoming[J]. Acta Ecologica Sinica, 2013, 33(5): 1524－1534.（in Chinese with English abstract）

[26] 楊偉州，邱碩，付喜廳，等.河北省生態(tài)功能區(qū)劃研究[J].水土保持研究，2016，23(4)：269－276. Yang Weizhou, Qiu Shuo, Fu Xiting, et al. Research for the ecological function regionalization of Hebei Province[J]. Research of Soil and Water Conservation, 2016, 23(4): 269－276. (in Chinese with English abstract)

[27] 劉志玲，李江風. 土地供應市場化對城市用地擴展的影響：以武漢市為例[J]. 國土資源科技管理，2009，26(2)：11－14. Liu Zhiling, Li Jiangfeng. On effect of land marketization on urban expansion of land use-with Wuhan City as an example[J]. Scientific and Technological Management of Land and Resources, 2009, 26(2): 11－14. (in Chinese with English abstract)

Construction of urban ecological security pattern and construction land development based on MCR Model

Qiu Shuo1,2,Wang Yuxin1※,Wang Pingzhi1, Lin Cong1

(1,,100083,; 2.,065400,)

In recent years, the importance of ecological space has been emphasized in the top-level design of urban land development in China. China is now in the stage of fast urbanization and the pace of industrial structure adjustment, the expansion of the construction land space is the most obvious sign of the transformation of the urban land landscape, and contradiction between cities and towns' competition and control of regional space is becoming more and more prominent. Identifying the urban ecological security pattern and discussing the spatial expansion pattern of construction land have become an urgent need for the implementation of the policy for delineating the red line of ecological protection, determining the boundaries of urban development and the “Three-planning united” spatial control. Using Tangshan City as a case, this paper was based on the theory of “Ecological process—Pattern” and the functions of spatial correlation, distance and superposition analysis of GIS. In this study, we integrated geostatistics analysis, spatial correlation analysis, distance analysis and overlay analysis. It obtained comprehensive ecological security pattern of Tangshan City by 5 single ecological security patterns: comprehensive hydrology security pattern, geological disaster pre-warning security pattern, biological habitat protection security pattern, soil and water conservation security pattern, and recreational security pattern. The integrated ecological security pattern under 3 different levels accounted for 15.31%, 50.08% and 12.70% of the total area respectively. In this paper, we built and improved MCR (minimal cumulative resistance) model, from which the weight coefficient of different levels of “source” was introduced. This paper took full account of the expansion capability of different extended sources. By using comprehensive ecological security pattern, attractiveness of the center of town and the main road and policy adjustment factor as resistance, and then the resistance surface of the regional space under different constraint conditions was obtained to determine the urban development scope under the condition of “Ecological security constraint pattern”, “Economic dominance pattern” and “Ecology & Economy win-win pattern of urban spatial expansion”, respectively. The spatial expansion of urban land was restricted by many factors, not only to consider the regional ecological and environmental factors, but also to pay special attention to the regional social and economic conditions and policy oriented factors. Through multi-scheme comparison, we found that the spatial expansion pattern driven by the conditions of “Ecology & Economy win-win pattern of urban spatial expansion” development model was the best model for the healthy and orderly development of cities and towns in the future, and was of great significance for optimizing the layout of land use and the protection of the ecological environment, moreover, which can effectively predict the scope of urban space development .On the basis of all above, considering the influence factors of urban development, the improved MCR was used to analyze the spatial expansion pattern of urban construction land use under different patterns, and to explore the appropriate ecological security pattern of urban development, which reduced the land conflict between urban development and environmental protection caused by large-scale and high-intensity urbanization and industrialization process, and put forward a scientific basis and reference for urban land layout planning and ecological red line delineation.

land use; remote sensing; ecological security pattern; construction land; expansion; minimal cumulative resistance model; Tangshan City

2018-04-07

2018-06-19

農業(yè)部“現代生態(tài)農業(yè)基地技術試驗示范項目”（2110402）

邱 碩，主要從事城鄉(xiāng)規(guī)劃與建設領域管理研究。Email：w6899366@126.com.

王宇欣，男，副教授，主要從事農業(yè)生物環(huán)境及能源工程研究。Emall：Wangyuxin@cau.edu.cn

10.11975/j.issn.1002-6819.2018.17.034

TU984

A

1002-6819(2018)-17-0257-09

邱 碩，王宇欣，王平智，林 聰. 基于MCR模型的城鎮(zhèn)生態(tài)安全格局構建和建設用地開發(fā)模式[J]. 農業(yè)工程學報，2018，34(17)：257－265. doi：10.11975/j.issn.1002-6819.2018.17.034 http://www.tcsae.org

Qiu Shuo, Wang Yuxin, Wang Pingzhi, Lin Cong. Construction of urban ecological security pattern and construction land development based on MCR Model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(17): 257－265. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.17.034 http://www.tcsae.org