何尹杰， 張占錄， 高文文

(中國人民大學(xué) 公共管理學(xué)院, 北京 100872)

規(guī)劃建設(shè)北京城市副中心，是以習(xí)近平總書記為核心的黨中央作出的重大決策部署。這不僅是優(yōu)化城市空間結(jié)構(gòu)，有序疏解非首都功能的需要，也是推動城鄉(xiāng)統(tǒng)籌協(xié)調(diào)發(fā)展、探索構(gòu)建新型城鎮(zhèn)化空間體系的需要[1]。隨著北京副中心選址確定在通州，通州的發(fā)展建設(shè)面臨著重大機(jī)遇，土地開發(fā)和環(huán)境整治等項目加快實施[2]。但近年來，隨著經(jīng)濟(jì)社會產(chǎn)業(yè)的快速發(fā)展，同時承接大量中心城區(qū)人口疏解，通州的土地利用問題也愈發(fā)顯著，人地矛盾日益突出。當(dāng)前學(xué)界關(guān)于政府搬遷的研究，可總結(jié)為3個方面，一是政府搬遷基礎(chǔ)研究[3-5]，如搬遷的意義，理論，類型，原因和后果等問題。二是政府搬遷對經(jīng)濟(jì)發(fā)展的影響[6-8]。三是政府搬遷對城市空間結(jié)構(gòu)的影響[9-11]。政府搬遷對城市空間結(jié)構(gòu)產(chǎn)生最直接的影響就是土地利用結(jié)構(gòu)的變化和建設(shè)用地擴(kuò)張[12-14]，從已有的研究看，政府搬遷前后土地利用變化進(jìn)行定量描述的研究較多，對搬遷后做出前瞻性預(yù)測的研究較少。

關(guān)于土地模擬預(yù)測的研究[15]，從土地利用變化動力機(jī)制[16]到模擬預(yù)測的方法模型[17]和演變規(guī)律總結(jié)[18]，無論是理論還是方法層面，國內(nèi)外都進(jìn)行了大量的探討。MCE-CA-Markov耦合模型不僅發(fā)揮了CA模型在各微觀因子之間相互作用機(jī)制方面的優(yōu)勢和Markov模型在時空尺度上的變化特征與規(guī)律的優(yōu)勢，而且MCE模型還可以充分考慮到政策影響因素。為使預(yù)測結(jié)果具有較高的精度與實際意義，本文結(jié)合《北京城市總體規(guī)劃(2016—2035年)》[19]與《北京城市副中心控制性詳細(xì)規(guī)劃(街區(qū)層面)(2016—2035年)》[20](以下簡稱《總規(guī)》和《詳規(guī)》)，選用MCE-CA-Markov耦合模型作為研究工具。與先前的研究相比，模擬預(yù)測方法更進(jìn)一步的地方在于：深入考慮規(guī)劃政策、自然條件、尤其是軌道交通等驅(qū)動因素對通州未來發(fā)展變化的影響，使其結(jié)果更具有可靠性。

1 研究區(qū)概況與數(shù)據(jù)處理

1.1 研究區(qū)地情概況

通州區(qū)地處北京市東南部，長安街延長線東端，京杭大運(yùn)河北端，地理坐標(biāo)在北緯39°36′—40°02′，東經(jīng)116°32′—116°52′，面積90 581.12 hm2(第二次全國土地調(diào)查數(shù)據(jù))，下轄23個街鎮(zhèn)鄉(xiāng)，南與天津、東與河北接壤。全區(qū)地處河流沖積平原，地勢平坦，平均海拔20 m。分布13條河流，總長245.3 km，主要河流有北運(yùn)河、潮白河、涼水河、鳳港減河[21]。北京行政副中心坐落通州行政區(qū)，位于長安街向東延55 km，面積約15 500 hm2，地理位置優(yōu)越。

1.2 數(shù)據(jù)源及預(yù)處理

通州區(qū)2009年，2013年、2017年三期土地利用矢量數(shù)據(jù)來源于北京市通州區(qū)自然資源局提供，該數(shù)據(jù)基于第二次全國土地利用調(diào)查及變更數(shù)據(jù)。2009—2017年通州區(qū)存在三級地類，根據(jù)研究的需要，參考2018年土地利用現(xiàn)狀(GBT21010-2018)分類系統(tǒng)以及結(jié)合研究區(qū)的土地利用特征，將部分地類進(jìn)行合并。由于通州區(qū)草地量較少，且與園地區(qū)分不明顯，故把兩類合并為一類?，F(xiàn)體系共分為以下6類：耕地、林地、園地(草地)、建設(shè)用地、水域和設(shè)施農(nóng)用地，具體的轉(zhuǎn)換見表1。地鐵數(shù)據(jù)采用論文[15]方法提取，高程和道路數(shù)據(jù)在地理空間數(shù)據(jù)云下載。計算各數(shù)據(jù)由矢量轉(zhuǎn)為柵格的面積和精度損失情況，確定最佳分辨率是30 m，所有數(shù)據(jù)進(jìn)行統(tǒng)一的投影、坐標(biāo)系、分辨率和范圍大小處理。

表1 土地利用轉(zhuǎn)換關(guān)系

2 分析方法

2.1 景觀格局變化

景觀格局變化分析同時也是對通州區(qū)2009—2017年土地利用變化情況進(jìn)行定量回顧。選取整體景觀水平(Landscape-level metrics)和類型水平(Class-level metrics)，整體景觀水平是對同一區(qū)域不同時期的整體景觀分析，而類型水平是對同一區(qū)域不同時期的不同利用類型斑塊景觀分析。前者強(qiáng)調(diào)對通州整體土地利用變化分析，而后者強(qiáng)調(diào)對通州某一種土地利用變化分析。

整體景觀水平上選擇11個指標(biāo)：斑塊數(shù)量(Number of Patches, NP)，景觀水平上NP表示圖斑總量，值越大說明破碎度越高，反之越低；斑塊密度(Patch Density, PD)，即斑塊個數(shù)與面積的比值，NP越高則PD越高；景觀分割指數(shù)(Landscape Division Index, DIVISION)表示不同類型斑塊的分離程度；邊緣密度(Edge Density，ED)反映斑塊的優(yōu)勢度和斑塊形狀的復(fù)雜程度，ED值大說明形狀越復(fù)雜；最大斑塊占景觀面積比例(Largest Patch Index，LPI)表示最大斑塊面積占景觀面積的比例；周長面積分維(Perimeter-Area Fractal Dimension，PAFRAC)數(shù)值在1～2，數(shù)值越小說明斑塊形狀越簡單，受人為干預(yù)程度大，反之越復(fù)雜；香濃多樣性(Shannon′s Diversity Index，SHDI)描述稀有圖斑的量，斑塊類型越豐富，SHDI值越大；整體性(斑塊凝聚度)指數(shù)(Patch Cohesion Index，COHESION)表示景觀黏合度，反映景觀空間均質(zhì)和連續(xù)的程度；蔓延度(Contagion，CONTAG)描述不同斑塊的集聚程度或延展趨勢，CONTAG越大說明優(yōu)勢斑塊形成了良好的連接性，反之說明多種類型并存，破碎化程度高；聚集指數(shù)(Aggregation Index，AI)描述同類型的斑塊公共邊界聚合程度；斑塊多度(景觀豐度)(Patch Richness，PR)，描述景觀的異質(zhì)性，即多少種用地類型。

類型水平上選擇9個指標(biāo)：斑塊數(shù)量(NP)、斑塊密度(PD)、最大斑塊占景觀面積比例(LPI)、邊緣密度(ED)、周長面積分維(PAFRAC)、斑塊所占景觀面積比例(PLAND)、整體性(斑塊凝聚度)指數(shù)(COHESION)、聚集指數(shù)(AI)、景觀分割指數(shù)(DIVISION)，這些指標(biāo)含義與整體景觀水平的指標(biāo)含義大體相同，區(qū)別在于類型水平上的指標(biāo)表示景觀中某種利用類型圖斑特征。

2.2 Markov模型(土地利用類型轉(zhuǎn)移矩陣)

Markov模型是根據(jù)現(xiàn)有的狀態(tài)水平來預(yù)測未來任意一個時間節(jié)點(diǎn)狀況的方法。轉(zhuǎn)移矩陣主要有兩方面特征：一是各土地利用類型轉(zhuǎn)移面積的信息；二是各土地利用類型轉(zhuǎn)移概率的信息。土地利用類型轉(zhuǎn)移矩陣可用來作為Markov鏈[22]。

2.3 MCE-CA-Markov耦合模型

MCE(multi-criteria evaluation)是一種多準(zhǔn)則評價的決策支持工具，通過一個多標(biāo)準(zhǔn)的評估，組合成適宜性圖集[23-25]。MCE，CA和Markov三者耦合建模動態(tài)模擬土地利用變化，有利于綜合考慮現(xiàn)實條件和自然要素的影響，能提高模擬的精度。對通州區(qū)6類用地類型分別制作MCE適宜性圖集，適宜性圖集分為兩類，限制性因素圖集和影響性因素圖集。在可預(yù)見的未來土地利用類型不會發(fā)生變化的區(qū)域，比如大部分河流水域和永久基本農(nóng)田，加入限制性因素圖集。在可預(yù)見的未來會影響土地利用類型變化的區(qū)域，比如交通路線附近，自然條件等，加入影響性因素圖集。

耕地圖斑：大部分河流水域設(shè)置限制因素，未來不會被開墾成為耕地。IDRISI軟件上把河流水域重分類為“0”，限制其轉(zhuǎn)變。影響因素設(shè)置，主要受地形因素影響。坡度越大高程越高，越不適宜耕種，并且坡度大于25°的地方禁止開墾為耕地。通州地勢平坦，坡度大多在2°～3°，高程普遍在10 m左右，大部分區(qū)域適宜耕種。

林地、園地(草地)圖斑：大部分河流水域和永久基本農(nóng)田設(shè)置限制因素，未來不會轉(zhuǎn)變?yōu)榱值亍@地(草地)。影響因素設(shè)置，由于通州區(qū)地形地勢平坦，而林地、園地(草地)相對于耕地和建設(shè)用地來說沒太苛刻條件，所以研究區(qū)域內(nèi)普遍適宜種植，沒有特殊地形設(shè)置。

水域圖斑：大部分河流水域被認(rèn)為未來不會有太大的變化，直接使用該圖斑的轉(zhuǎn)移矩陣作為適宜性圖。

設(shè)施農(nóng)用地圖斑：設(shè)施農(nóng)用地是指直接用于經(jīng)營性養(yǎng)殖、栽培或晾曬等以及其相關(guān)的附屬設(shè)施用地。根據(jù)2019年12月自然資源部關(guān)于該類用地管理的通知[26]，第二條提到，設(shè)施農(nóng)用地占用一般耕地屬于農(nóng)業(yè)內(nèi)部結(jié)構(gòu)調(diào)整，無需要占補(bǔ)平衡。涉及少量永久基本農(nóng)田確實難以避讓的，可以使用永久基本農(nóng)田，沒破壞耕作層的，無需補(bǔ)劃，破壞了則必須補(bǔ)劃。設(shè)施農(nóng)用地應(yīng)盡量選擇在難以利用和低效閑置的土地上。所以對設(shè)施農(nóng)用地不進(jìn)行水域和基本農(nóng)田的限制因素設(shè)置。影響因素設(shè)置，主要由地形因素影響。根據(jù)規(guī)定，設(shè)施建設(shè)應(yīng)盡量利用荒山荒坡，平整的區(qū)域盡量留給耕地，所以設(shè)置影響因素時，把不適宜耕種的區(qū)域設(shè)為設(shè)施農(nóng)用地發(fā)展區(qū)。

建設(shè)用地圖斑：限制因素的設(shè)置，把大部分河流水域、永久基本農(nóng)田、生態(tài)管控區(qū)3種類型區(qū)域重分類成“0”，限制建設(shè)用地侵蝕。影響因素設(shè)置，主要由交通因素，地形因素決定。交通因素影響：(1) 地鐵因素，目前通過通州的地鐵只有三條線路，分別是6號線，八通線，和亦莊線，根據(jù)《總規(guī)》(2016—2035年)2021年市域軌道交通規(guī)劃示意圖，預(yù)計經(jīng)過通州的地鐵還將建成17號線、7號線東延、八通線南延、22號線(平谷線)，新增加4條線路。根據(jù)論文[27]認(rèn)為國內(nèi)外地鐵線路普遍對2 000 m范圍內(nèi)產(chǎn)生影響，地鐵對建設(shè)用地的影響函數(shù)成“N”型曲線，200～800 m范圍內(nèi)是建設(shè)用地的最佳選址位置，少于200 m由于受到人流、噪音和污染的影響，不太適合建設(shè)用地發(fā)展，2 000 m后由于通勤到地鐵的距離過遠(yuǎn)，地鐵的影響則不是首位的。所以“N”曲線上的控制點(diǎn)分別設(shè)為0，200，800，2 000。(2) 鐵路因素，根據(jù)《總規(guī)》市域客運(yùn)樞紐體系規(guī)劃圖，目前通過通州區(qū)域的高速鐵路有：京津城際鐵路、城際鐵路聯(lián)絡(luò)線和小段京秦線。普通鐵路線路有：小段的京承線。到2035年內(nèi)高速鐵路和普通鐵路均無遠(yuǎn)期路線。根據(jù)我國鐵路保護(hù)條例[28]，第十條鐵路線路兩側(cè)應(yīng)當(dāng)設(shè)立鐵路線路安全保護(hù)區(qū)，安全距離為15 m左右。鐵路沿線兩側(cè)50 m內(nèi)為鐵路用地，鐵路與車站有關(guān)系，一般認(rèn)為市區(qū)開車到車站最佳時間在2 h之內(nèi)，如果按照市區(qū)車速30 km/h算，60 km將是鐵路影響最大值，鐵路對建設(shè)用地的影響函數(shù)曲線因素成“N”型曲線，曲線上的控制點(diǎn)分別設(shè)為15，100，1 000，60 000。(3) 道路因素，由于噪音和汽車尾氣的污染，一般認(rèn)為距離道路200 m左右是最適宜建設(shè)用地發(fā)展的，800 m以后則道路的影響減弱。道路對建設(shè)用地的影響函數(shù)曲線成“S”型曲線，曲線上的控制點(diǎn)分別設(shè)置為200 m和800 m。地形因素影響，由于通州的地勢平坦，地形因素對建設(shè)用地影響較小。

建設(shè)用地圖斑影響因素由于較多，利用層次分析法(AHP)來進(jìn)行因子權(quán)重賦值。AHP層次分析法是用來確定多個因子權(quán)重中最常見且有效的方法[29]，首先構(gòu)建層次結(jié)構(gòu)模型，將決策目標(biāo)分解成多個準(zhǔn)則層要素，準(zhǔn)則層要素再分解為多個備選方案3個層次；然后判斷矩陣，在準(zhǔn)則層要素的基礎(chǔ)之上對備選方案之間的重要性進(jìn)行判斷，判斷矩陣需通過一致性檢驗；最后得出計算結(jié)果，即各個備選方案對決策目標(biāo)的影響權(quán)重大小。把影響評價結(jié)果的高程、坡度、鐵路、地鐵、主要道路，5個因子看作是備選方案要素，兩要素之間在建設(shè)適宜度和可能性的準(zhǔn)則下對比其重要程度，得到相關(guān)的數(shù)據(jù)矩陣，確定建設(shè)用地影響因素指標(biāo)權(quán)重，結(jié)果是高程為12.26%，坡度為28.26%，鐵路為12.26%，地鐵為28.26%，主要道路為18.96%，檢驗一致性合格，得出建設(shè)用地的適宜性圖集。

3 結(jié)果與分析

3.1 景觀格局變化分析

3.1.1 整體景觀水平分析 由表2可知、2009—2017年、通州區(qū)景觀破碎化程度開始逐漸降低。斑塊數(shù)量(NP)減少了38.62%，17 431塊；斑塊密度(PD)也隨之降低；景觀分割指數(shù)(DIVISION)減少，表明斑塊間隔離度降低。邊緣密度(ED)2013—2017年減少較快，斑塊中心面積不斷減少。2009—2013年最大斑塊占景觀面積比例(LPI)大幅度增加，表明這4 a基礎(chǔ)設(shè)施建設(shè)用地加速推進(jìn)。周長面積分維數(shù)(PAFRAC)下降和香濃多樣性(SHDI)減少表明發(fā)展呈集聚趨勢，土地受人為干預(yù)，集中利用。蔓延度(CONTAG)指數(shù)和聚集指數(shù)(AI)上升，二者都表明了優(yōu)勢斑塊團(tuán)聚程度不斷提高，尤其在2013—2017年更為明顯，體現(xiàn)了土地受到了集中連片地開發(fā)。斑塊多度(景觀豐度)(PR)沒有發(fā)生變化、2009—2017年有6種土地利用類型。

表2 2009-2017年通州區(qū)整體景觀水平指標(biāo)變化

3.1.2 類型水平動態(tài)分析 聚散性指標(biāo)變化特征：斑塊數(shù)量(NP)，耕地和設(shè)施農(nóng)用地的斑塊數(shù)量2009—2017年增多，表明兩種用地類型被人為干擾嚴(yán)重，破碎化程度上升。建設(shè)用地斑塊數(shù)減少，但面積卻增加，表明城市化過程中，通州的建設(shè)用地集中連片開發(fā)。斑塊密度(PD)，耕地和設(shè)施農(nóng)用地的斑塊密度上升，表明破碎化程度增加。建設(shè)用地和水域的斑塊密度大幅度下降，表明兩種地類趨向均一。各用地類型的整體性(斑塊凝聚度)指數(shù)(COHESION)和景觀分割(DIVISION)指數(shù)變化不大，除設(shè)施農(nóng)用地小幅度下降以外；聚集指數(shù)(AI)普遍上升，表明類型分散程度減少，反映了人類的活動對通州景觀格局的特征影響增強(qiáng)。

面積—邊緣指標(biāo)變化特征：園地(草地)、林地和水域的邊緣密度(ED)值下降，表明單位面積的邊界長度減少，斑塊的形狀趨于單一和規(guī)則。最大斑塊占景觀面積比例(LPI)，除了建設(shè)用地，其余的地類指數(shù)都下降，表明其余地類在景觀的優(yōu)勢度下降。

形狀指標(biāo)變化特征：所有地類的周長面積分維(PAFRAC)相差不大，建設(shè)用地和耕地是景觀中占優(yōu)勢的主導(dǎo)地類，被當(dāng)成景觀基質(zhì)來處理，其他規(guī)模較小的地類斑塊鑲嵌其中(表3)。

表3 2009-2017年通州區(qū)類型水平景觀指數(shù)變化

3.2 土地利用類型轉(zhuǎn)移矩陣分析

將研究的初始時間點(diǎn)2009年的圖層和末期時間點(diǎn)2017年的圖層進(jìn)行疊加分析，得出通州的土地利用變化情況。結(jié)果顯示，2009—2017年，通州區(qū)有4 614.24 hm2土地利用類型發(fā)生變化，占總面積的5.09%(圖1)。耕地主要流向建設(shè)用地、設(shè)施農(nóng)用地，轉(zhuǎn)化的面積比例分別為4.96%和1.18%；同時，部分建設(shè)用地、林地和園地(草地)轉(zhuǎn)化為耕地。建設(shè)用地的增加主要來自于耕地、林地和水域。其中耕地的轉(zhuǎn)換面積比例為4.91%(表4)。2009年通州區(qū)各類土地利用類型轉(zhuǎn)化為2017年的各類土地利用類型的概率見表5，各類土地利用類型轉(zhuǎn)化過程中，轉(zhuǎn)化為建設(shè)用地的可能性最大；轉(zhuǎn)化為園地(草地)的可能性最小。

3.3 土地利用變化模擬預(yù)測

疊加2009年、2013年土地利用圖，得到土地利用類型轉(zhuǎn)移概率矩陣和轉(zhuǎn)移面積矩陣。采用5×5的濾波器(即5 km×5 km的矩形范圍內(nèi)的元胞對中心元胞的狀態(tài)改變有顯著影響)，模擬預(yù)測得出2017年的土地利用圖，和現(xiàn)有真實的2017年土地利用現(xiàn)狀圖比較，得出模擬精度。調(diào)整適宜性圖集，提高模擬精度(Kappa系數(shù))到90%之后，利用2009年與2017年現(xiàn)實的土地利用圖像選擇CA循環(huán)次數(shù)，模擬2035年的土地利用狀況。

圖1 2009-2017年通州區(qū)土地利用類型變化位置

預(yù)測結(jié)果如圖2所示：2035年通州區(qū)耕地面積為30 208.80 hm2，占33.35%，林地為9 040.00 hm2，占9.98%，園地(草地)為3 378.68 hm2，占3.74%，水域為9 529.13 hm2，占10.52%，建設(shè)用地為2 309.82 hm2，占39.87%，設(shè)施農(nóng)用地為2 309.82 hm2，占2.55%。其中，城市副中心的建設(shè)用地為11 234.40 hm2，占副中心區(qū)域的72.48%；亦莊新城(通州部分)建設(shè)用地為6 106.10 hm2，占該區(qū)域的93.94%。2009—2035年、通州區(qū)的各種土地類型的面積發(fā)生了較大變化。耕地、園地(草地)減少，而建設(shè)用地、林地、水域以及設(shè)施農(nóng)用地增加。其中，耕地減少了4 984.59 hm2，所占比重由38.85%下降到33.34%，下降5.51%。建設(shè)用地增加了3 773.67 hm2，所占比重由35.70%增長到39.87%，增長4.16%。由表6可知，耕地、林地、園地(草地)、水域、建設(shè)用地、設(shè)施農(nóng)用地的用地結(jié)構(gòu)由2009年的38.87%，9.36%，4.22%，9.93%，35.69%，1.90%，到2035年變?yōu)?3.35%，9.98%，3.73%，10.52%，39.87%?？臻g上，耕地被周圍的建設(shè)用地侵蝕；園地(草地)減少區(qū)域主要位于通州區(qū)東部的西集鎮(zhèn)和漷縣鎮(zhèn)；水域增加區(qū)域基本在中南部張家灣和永樂店新市鎮(zhèn)一些人工的水產(chǎn)養(yǎng)殖；設(shè)施農(nóng)用地增長區(qū)域多分布在北部宋莊靠近順義區(qū)，林地零散分布在通州各個區(qū)域，均有所增長；建設(shè)用地增長區(qū)域主要在于北京城市副中心、臺湖和亦莊新城。

表4 2009-2017年通州區(qū)土地利用類型面積轉(zhuǎn)移矩陣 hm2

表5 2009-2017年通州區(qū)土地利用類型轉(zhuǎn)移概率矩陣 %

圖2 通州區(qū)2035年土地利用類型模擬

3.4 模擬預(yù)測結(jié)果與規(guī)劃對比

按照《總規(guī)》和《詳規(guī)》(2016—2035年)城市空間布局的構(gòu)想，北京城市副中心和外圍控制區(qū)即通州全區(qū)，將承接首都核心區(qū)和主城區(qū)的人口的疏解和相關(guān)非首都功能的轉(zhuǎn)移。到2035年城市副中心人口規(guī)?？刂圃?30萬人以內(nèi)，通州區(qū)人口規(guī)?？刂圃?00～205萬人以內(nèi)，形成通州—城市副中心—亦莊新城(通州部分)—鎮(zhèn)—新型農(nóng)村社區(qū)的空間體系。提取《總規(guī)》和《詳規(guī)》(2016—2035年)中關(guān)于土地利用結(jié)構(gòu)和規(guī)模的核心指標(biāo)數(shù)據(jù)，與模擬預(yù)測結(jié)果對比。預(yù)測2035年通州區(qū)的城鄉(xiāng)建設(shè)用地規(guī)模為36 114.69 hm2，開發(fā)強(qiáng)度為39.87%，低于《總規(guī)》規(guī)定的2035年北京市平原地區(qū)開發(fā)強(qiáng)度44%，但離《詳規(guī)》明確規(guī)定的30.36%還存在一定的差距，要實現(xiàn)《詳規(guī)》目標(biāo)，則必須嚴(yán)格執(zhí)行規(guī)劃提到的1∶0.8平均拆占比，守住城鎮(zhèn)開發(fā)邊界紅線，使建設(shè)用地多減少增。按2035年城市副中心130萬人規(guī)模預(yù)估，預(yù)測中2035年人均綠地面積離預(yù)期性指標(biāo)仍有一定的差距，要達(dá)到預(yù)期性指標(biāo)則需要控制人口規(guī)模或增加綠地面積；副中心的建設(shè)用地預(yù)測規(guī)模為72.48%，減去規(guī)劃中留出的區(qū)域建設(shè)用地(特指特殊用地和對外交通用地)的6.4%，剩下城鄉(xiāng)建設(shè)用地規(guī)模66.08%，比規(guī)劃的指標(biāo)超額了1.48%，在可控制范圍之內(nèi)；生態(tài)綠地及農(nóng)林地高出預(yù)期性。亦莊新城(通州部分)的城鄉(xiāng)建設(shè)用地規(guī)模高出指標(biāo)4.71%，未來需做好減量發(fā)展，按照創(chuàng)新驅(qū)動發(fā)展的要求，不隨意擴(kuò)張建設(shè)用地規(guī)模。

表6 通州區(qū)2009—2035年土地利用類型面積/占比變化趨勢

4 結(jié) 論

景觀總體水平上，通州區(qū)總體景觀破碎化程度減少，土地按規(guī)劃集中連片使用。景觀類型水平上，建設(shè)用地和耕地是主導(dǎo)優(yōu)勢地類。耕地破碎化趨勢明顯，林地、建設(shè)用地集中連片開發(fā)，園地(草地)保持在較穩(wěn)定的水平。

2009—2017年，通州區(qū)有4 614.24 hm2土地利用類型發(fā)生變化，占比5.09%，耕地是主要的流失類型，流向建設(shè)用地和設(shè)施農(nóng)用地；建設(shè)用地的增加主要來自于耕地、林地和水域。各類土地利用類型轉(zhuǎn)化過程中，轉(zhuǎn)化為建設(shè)用地的可能性最大；轉(zhuǎn)化為園地(草地)的可能性最小。

應(yīng)用MCE-CA-Markov模型，對通州區(qū)土地利用進(jìn)行了模擬預(yù)測。在規(guī)劃期限內(nèi)非基本農(nóng)田的耕地面臨著被城市建設(shè)用地擴(kuò)張侵蝕的壓力，耕地、園地(草地)減少，而林地、建設(shè)用地、水域、以及設(shè)施農(nóng)用地增加。北京城市副中心、臺湖和亦莊新城三處建設(shè)用地有蔓延、連片發(fā)展的趨勢，開發(fā)強(qiáng)度都高達(dá)到60%以上。

模擬預(yù)測結(jié)果與規(guī)劃對比，建設(shè)用地預(yù)測的值比規(guī)劃指標(biāo)高，需要注意；但生態(tài)用地指標(biāo)超出了預(yù)期值，呈現(xiàn)良好的保護(hù)趨勢。未來需防止建設(shè)用地過快、無序地擴(kuò)張，應(yīng)嚴(yán)格執(zhí)行生態(tài)保護(hù)、永久基本農(nóng)田和城鎮(zhèn)開發(fā)邊界三條紅線，落實用途管制。