鄧雅文，林 耕，曾穗平

（天津城建大學(xué) 建筑學(xué)院，天津300384）

土地利用變化是生態(tài)、社會(huì)和政策等多因素綜合作用下的結(jié)果，隨著城市化和工業(yè)化的進(jìn)程不斷加速，深入研究土地利用變化特征并模擬預(yù)測未來的變化趨勢是十分重要的.因此，提升城市土地利用效率，推進(jìn)城鎮(zhèn)建設(shè)高水平發(fā)展是當(dāng)前國土空間優(yōu)化的有效途徑，且運(yùn)用科學(xué)技術(shù)模擬和預(yù)測地類變化過程是國土空間規(guī)劃編制的重要基礎(chǔ).當(dāng)前，有以下模型（見表1）可運(yùn)用于土地利用變化的預(yù)測和分析中，單一模型均有其優(yōu)缺點(diǎn)，不同模型之間的耦合建立較為常用.其中CA-Markov模型高效地集合了CA模型的空間預(yù)測技術(shù)與Markov模型的長時(shí)間序列預(yù)測模擬的優(yōu)點(diǎn)，有助于提高研究的精確性，能夠更準(zhǔn)確地預(yù)測土地利用變化，在目前的土地動(dòng)態(tài)預(yù)測中應(yīng)用較廣.相關(guān)研究有國外的泰國華欣海濱城市[1]、埃塞俄比亞的Andassa河流域[2]、印度拉敦市[3]；國內(nèi)研究聚焦于珠江三角洲地區(qū)[4]、黃土丘陵區(qū)[5]、河西走廊[6]、黑河流域[7]和新疆奎屯河流域[8]等.

表1 土地利用變化預(yù)測與分析模型

運(yùn)用CA-Markov耦合模型研究和預(yù)測天津市土地利用時(shí)空變化便于理解大城市發(fā)展的趨勢，對(duì)于了解研究區(qū)的土地利用非線性特征有理論意義.同時(shí)，研究基于該模型對(duì)把握未來大城市土地利用的變化趨勢，改進(jìn)土地利用的結(jié)構(gòu)與規(guī)劃，解決目前大城市土地利用相關(guān)問題具有實(shí)踐意義.當(dāng)前社會(huì)對(duì)于土地利用時(shí)空變化的模擬預(yù)測大多集中于兩期較大尺度的預(yù)測或從整體區(qū)域規(guī)劃的角度研究土地利用數(shù)據(jù)[9]，再進(jìn)行預(yù)測，且精度評(píng)價(jià)上多為單一kappa系數(shù)評(píng)價(jià)，使得結(jié)果不夠準(zhǔn)確.因此，文章選取三個(gè)時(shí)期的土地利用變化情況，利用CA-Markov耦合模型模擬近二十年土地利用動(dòng)態(tài)變化情況，并以擴(kuò)展的kappa系數(shù)為精度評(píng)價(jià)方法，使得研究和預(yù)測結(jié)果更加準(zhǔn)確.根據(jù)研究可在一定程度上理解天津市的土地利用變化情況，并針對(duì)預(yù)測后可能出現(xiàn)的問題提出相關(guān)策略，為類似大城市的土地利用發(fā)展提供科學(xué)的理論支撐與依據(jù)[10].

1 研究區(qū)與數(shù)據(jù)

1.1 研究區(qū)概況

天津市作為中國北方最大的沿海開放城市，2005—2019年城市化率從75.73%提高到83.48%，迅速發(fā)展的城市化對(duì)于用地資源的需求量逐漸增加，這使得農(nóng)業(yè)用地迅速減少.從2005—2018年，天津市農(nóng)業(yè)用地從7 075 km2減少到6 894.41 km2（見圖1），耕地面積在2011年到2012年變化速率達(dá)到最快，減少最多，整體呈現(xiàn)“急劇下降—趨于平緩”的走勢；建設(shè)用地在2012—2013增長速率較快，整體為“增長—下降—緩慢增長—趨于穩(wěn)定”的趨勢.與此同時(shí)，耕地土壤也受到了不同程度的污染，限制了天津的社會(huì)經(jīng)濟(jì)發(fā)展.為了提高發(fā)展速度，獲得較大的發(fā)展機(jī)遇，天津市各類用地利用類型相互轉(zhuǎn)換頻繁，人地矛盾漸漸凸顯.

圖1 2011—2018年天津市城鄉(xiāng)建設(shè)用地?cái)?shù)據(jù)

根據(jù)天津市統(tǒng)計(jì)年鑒數(shù)據(jù)資料顯示，2005—2018年，天津市農(nóng)業(yè)人口（見圖2）從376.91萬人減少到321.30萬人，城市人口從564.2萬人增加至760.33萬人，2015年前后農(nóng)業(yè)人口急劇下降，城鎮(zhèn)人口明顯上升，城市化率顯著提升，表明天津市人口聚集現(xiàn)象明顯，城市化進(jìn)程發(fā)展速度快.

圖2 2011—2018年天津市城鄉(xiāng)人口數(shù)據(jù)

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

本文所用1995年、2005年和2015年天津土地利用現(xiàn)狀數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(tái)（http：//www,gscloud.cn）.根據(jù)土地利用2015年的一級(jí)分類標(biāo)準(zhǔn)，將天津市用地類型劃分為9類用地：耕地、草地、水域、林地、未利用地、城鎮(zhèn)居住用地、農(nóng)村居民點(diǎn)、其他建設(shè)用地.研究所得到的研究區(qū)數(shù)據(jù)為矢量數(shù)據(jù)，轉(zhuǎn)為柵格數(shù)據(jù)之后將每期土地利用數(shù)據(jù)柵格像元大小、坐標(biāo)系、投影以及空間范圍調(diào)為一致，在Arcgis中轉(zhuǎn)為ascii文件導(dǎo)入IDRISI軟件.

2 研究方法

本文研究方法（見圖3）基于地理信息數(shù)據(jù)和社會(huì)經(jīng)濟(jì)數(shù)據(jù)，通過多時(shí)段分析和多種技術(shù)方法得出天津市三期土地利用文件，并運(yùn)用CA-Markov模型模擬近20年土地利用時(shí)空變化，提煉對(duì)比得出天津市三期城市擴(kuò)展方式，各類用地時(shí)空變化特征.

圖3 研究方法

2.1 土地利用動(dòng)態(tài)度

土地利用動(dòng)態(tài)度能夠定量描述某一區(qū)域的土地利用面積變化數(shù)量、速度以及差異情況[11]，其表達(dá)式如下

其中：U表示某一土地利用類型占地面積；a、b為研究開始與結(jié)束的時(shí)間節(jié)點(diǎn)；T為總時(shí)間.

2.2 元胞自動(dòng)機(jī)（CA）

CA是擁有時(shí)間空間運(yùn)算技術(shù)的網(wǎng)絡(luò)動(dòng)態(tài)系統(tǒng)模型，有時(shí)空狀態(tài)都呈離散形式并且元胞的影響都為局部的特征[12]，該技術(shù)運(yùn)用于模擬Arcgis中復(fù)雜的時(shí)空演化過程.CA模擬系統(tǒng)由狀態(tài)、領(lǐng)域、元胞和轉(zhuǎn)換規(guī)則構(gòu)成，一般運(yùn)用5×5的摩爾領(lǐng)域作為濾波器，遵循中心元胞狀態(tài)轉(zhuǎn)換影響與距離為反比的轉(zhuǎn)換規(guī)則.CA的轉(zhuǎn)換規(guī)則中每一個(gè)元胞都基于局部的規(guī)則進(jìn)行轉(zhuǎn)換，其模型轉(zhuǎn)換式為

式中：f為局部空間的元胞轉(zhuǎn)化規(guī)則；S為元胞有限、離散狀態(tài)的幾何；t、t+1為不同時(shí)間點(diǎn)；N為元胞的相鄰區(qū)域.但是往往對(duì)土地利用變化帶來影響的因素都是非線性的，只運(yùn)用CA局部的轉(zhuǎn)換規(guī)則較難對(duì)城市發(fā)展進(jìn)行精準(zhǔn)預(yù)測.

2.3 Markov模型

Markov是一個(gè)隨機(jī)系統(tǒng)，必須遵循兩個(gè)條件：一是系統(tǒng)在任意時(shí)間點(diǎn)都可以用有限的狀態(tài)之一來定義描述；二是，系統(tǒng)無后效性，在某一個(gè)階段的狀態(tài)只要確定之后，之后的演變都不再受之前的各種因素影響.土地利用變化預(yù)測符合馬爾可夫鏈的條件，具有Markov性，不同地類之間相互轉(zhuǎn)換的比例就為狀態(tài)轉(zhuǎn)移概率，利用這個(gè)概率就可建立Markov模型，公式如下

式中：Pij（i，j=1，2，…，n）代表時(shí)間循環(huán)周期中地類從i轉(zhuǎn)為j的轉(zhuǎn)移概率；P為各地類轉(zhuǎn)移概率矩陣；n為用地類型的數(shù)量；St、St+1為t時(shí)間點(diǎn)與t+1時(shí)間點(diǎn)用地狀態(tài).Markov鏈所進(jìn)行的全局約束是時(shí)間的函數(shù)，能夠有效合理地模擬長時(shí)間序列土地利用類型數(shù)量的變化.

2.4 CA-Markov模型

CA與Markov耦合模型，在基礎(chǔ)Markov模型之上加入了空間權(quán)重因素，不僅擁有CA模擬復(fù)雜時(shí)空變化的能力，同時(shí)又具備Markov模型長期進(jìn)行預(yù)測的優(yōu)點(diǎn)[13]，能夠保證土地利用轉(zhuǎn)變?yōu)槟壳艾F(xiàn)存的土地地類，而非完全碎金生成與分布，CA-Markov有助于提高該預(yù)測的精度.通過CA-Markov模塊預(yù)測研究對(duì)象未來發(fā)展趨勢與土地利用布局，主要步驟為如下.

（1）在Arcgis中將處理好的矢量文件轉(zhuǎn)為柵格文件，運(yùn)用面轉(zhuǎn)柵格功能將柵格文件轉(zhuǎn)為ASCII文本，之后再運(yùn)用IDRISI軟件中的ARCRASTER模塊將ASCII文件轉(zhuǎn)為.rst格式文件，之后對(duì)轉(zhuǎn)換之后的數(shù)據(jù)進(jìn)行RECLASS重分類處理.

（2）運(yùn)用IDRISI軟件中的Markov模塊，處理獲得1995—2005年，2005—2015年轉(zhuǎn)移面積矩陣，之中轉(zhuǎn)移面積矩陣作為轉(zhuǎn)換規(guī)則參與運(yùn)算；之后運(yùn)用Decision Wizard模塊進(jìn)行城鄉(xiāng)工礦居民用地、與耕地的適宜性圖運(yùn)算，以高程、坡度、道路通達(dá)性為影響因子形成適應(yīng)性圖，與轉(zhuǎn)移面積矩陣圖一起形成適宜性圖集，參與CA-Markov運(yùn)算過程.

（3）將CA的循環(huán)次數(shù)設(shè)置為10次，在IDRISI軟件集成的CA-Markov模塊中，將元胞自動(dòng)機(jī)迭代次數(shù)設(shè)為10，預(yù)測2015年土地利用格局；之后與2015年的實(shí)際土地利用進(jìn)行對(duì)比，進(jìn)行精度檢驗(yàn)，若符合要求，進(jìn)行研究區(qū)2025年土地利用格局預(yù)測.

2.5 模擬精度驗(yàn)證

常用的精度驗(yàn)證方法有傳統(tǒng)的kappa系數(shù)檢驗(yàn)、隨機(jī)檢驗(yàn)法、全數(shù)檢驗(yàn)法與擴(kuò)展的kappa系數(shù)檢驗(yàn)[14]等.國外著名教授提出一種更為準(zhǔn)確的拓展kappa指數(shù)，包含4種擴(kuò)展的kappa指數(shù)，為：kstandard，簡稱標(biāo)準(zhǔn)kappa指數(shù)，用來評(píng)價(jià)綜合變化信息；kno，簡稱隨機(jī)kappa指數(shù)，用來評(píng)價(jià)土地利用綜合信息變化klocation，簡稱位置kappa指數(shù)，運(yùn)用擴(kuò)展的kappa指數(shù)能夠得出2個(gè)時(shí)間點(diǎn)之間的動(dòng)態(tài)變化情況，使得綜合評(píng)價(jià)模型在土地利用格局模擬中精度更加準(zhǔn)確.

本文采用擴(kuò)展的kappa系數(shù)對(duì)研究區(qū)2015年的模擬數(shù)據(jù)與實(shí)際數(shù)據(jù)進(jìn)行精度驗(yàn)證，進(jìn)而得出CAMarkov模型在模擬天津市土地利用格局中的準(zhǔn)確性與適用性，如符合要求，能夠更好地模擬2025年天津市土地利用格局情況.

3 結(jié)果與分析

3.1 1995—2015年天津市用地變化

研究從通過ArcGis10.2分析獲取的天津市3期土地利用變化信息（見圖4）中可以得出：1995—2015年城鎮(zhèn)居住用地面積（見表2）從533.55 km2增加到1 153.35 km2，增加量為619.80 km2，20年來城鎮(zhèn)居住用地發(fā)展迅速，增長率由50.65%上升到116.16%，天津城市化發(fā)展進(jìn)程推進(jìn)速度快；農(nóng)村居民點(diǎn)用地從871.62 km2增加到1 089.05 km2，增加量為217.43 km2，2005—2015年農(nóng)村居民點(diǎn)增長率高于1995—2005年15個(gè)百分點(diǎn)，2005年之后城鄉(xiāng)居住用地增長速度逐步加快.其他建設(shè)用地為先減少后增加的變化特征，主要變化集中于濱海新區(qū)，2005年，濱海新區(qū)被寫入“十一五”規(guī)劃并納入國家發(fā)展戰(zhàn)略，成為國家重點(diǎn)支持開發(fā)開放的國家級(jí)新區(qū).2005年后濱海新區(qū)通過填海造陸等舉措，在渤海沿岸擴(kuò)建建設(shè)用地，帶動(dòng)了區(qū)域內(nèi)的經(jīng)濟(jì)發(fā)展以及與中心城區(qū)的溝通交流.其他用地均有不同程度減少，變化劇烈程度呈現(xiàn)“水域＞草地＞林地＞耕地”的特征（見表3），草地、林地等生態(tài)用地的減少，可能帶來不同程度的生態(tài)問題，如水資源短缺，空氣污染嚴(yán)重等.

圖4 天津市1995—2015年土地利用變化

表2 天津市1995—2015年土地利用面積

表3 天津市1995—2015年土地利用變化面積、變化率

3.2 2015年天津市土地利用格局模擬

3.2.1 土地利用轉(zhuǎn)移概率矩陣

利用CA-Markov可得天津市1995—2005年和2005—2015年地類轉(zhuǎn)移概率矩陣（見表4-5）.

表4 天津市1995—2005年土地利用轉(zhuǎn)移概率

通過分析得出的土地利用轉(zhuǎn)移概率等信息，1995—2015年20年內(nèi)城鎮(zhèn)用地、農(nóng)村居民點(diǎn)出現(xiàn)持續(xù)增加的趨勢，主要來源于耕地、草地和水域.從轉(zhuǎn)換頻繁程度來看，天津市各類用地轉(zhuǎn)換頻繁，其中城鄉(xiāng)建設(shè)用地面積增加明顯，其他用地都以轉(zhuǎn)出為主，土地利用轉(zhuǎn)換主要存在于城市中部以及西南沿海地區(qū).從發(fā)展區(qū)域來看，濱海新區(qū)的創(chuàng)建對(duì)于天津城市發(fā)展有一定的積極影響，在2005—2015年影響尤為明顯，隨著輕軌，津?yàn)I大道等交通設(shè)施的發(fā)展，天津中心城區(qū)與濱海新區(qū)逐漸連為一體.從整體來看，天津城市擴(kuò)展方式為中心城區(qū)著重發(fā)展，市郊區(qū)以中心城區(qū)為原點(diǎn)放射性發(fā)展，城市化進(jìn)程速度較快.其空間擴(kuò)張形式為攤餅式發(fā)展為主，放射式低密度蔓延.

1995—2015年期間，城鎮(zhèn)用地、農(nóng)村居民點(diǎn)和耕地是天津市土地利用的主要3種類型，三種用地占天津市整體用地的78.10%.20年間天津市城鎮(zhèn)用地和農(nóng)村居民點(diǎn)用地面積逐年遞增，但除了城市中心以及西南部用地以外，農(nóng)村居民點(diǎn)多呈離散式分布，郊區(qū)區(qū)域發(fā)展較快但閑置率高，整體區(qū)域的結(jié)構(gòu)與用地布局有待調(diào)整；耕地面積變化整體上是負(fù)向變化，部分轉(zhuǎn)化為建設(shè)用地，變化面積為-34.79 km2.

3.2.2 2015年天津市土地利用適宜性圖集

利用MCE方法構(gòu)建部分土地利用適宜性圖集（見圖5）作為CA的補(bǔ)充轉(zhuǎn)換規(guī)則，這里以城鎮(zhèn)用地為例，制城鎮(zhèn)用地第一至二期土地利用適宜圖.在保護(hù)生態(tài)水源的基礎(chǔ)上，選擇以水域?yàn)橄拗茥l件，選擇坡度、高程、道路通達(dá)性為約束因子，運(yùn)用模糊隸屬度函數(shù)中的J型遞減函數(shù)，運(yùn)用AHP方法得出各項(xiàng)權(quán)重，最后得出城鎮(zhèn)用地適宜性圖.從中可以看出，城鎮(zhèn)用地主要分布于中部地區(qū)與東南部濱海新區(qū)部分，天津北部高程、坡度變化較大地區(qū)較不適宜居住.適宜性圖集中的像素范圍值為0~255，0表示不合適建設(shè)該地類，255表示高度合適.

圖5 天津市第一、二期適宜性圖集

3.2.3 2015年天津市土地利用變化模擬

基于CA-Markov耦合模型，以天津市2005年二期土地利用圖為基礎(chǔ)圖（見圖6），結(jié)合適宜性圖集進(jìn)行預(yù)測，將元胞迭代次數(shù)設(shè)置為10 a，可得2015年模擬土地利用圖（見圖7）；并通過IDRISI軟件中的VADLIDATE模塊將2015年預(yù)測用地與實(shí)際用地進(jìn)行精度檢驗(yàn)，VALIDATE結(jié)果得出擴(kuò)展的4種kappa系數(shù)分別為：kno=0.838 2；kstandard=0.833 4；klocation=0.829 5，klocationstrata=0.829 5；kappa系數(shù)=0.831 4.數(shù)據(jù)表示2015年天津市土地利用模擬結(jié)果與2015年天津市實(shí)際土地利用情況一致性較高，預(yù)測結(jié)果較好，可以用于后續(xù)的預(yù)測.

圖6 天津市2015年土地利用類型

圖7 天津市2025年土地預(yù)測類型

表5 天津市2005—2015年土地利用轉(zhuǎn)移概率

3.3 2025年天津市用地格局預(yù)測

根據(jù)CA-Markov模型模擬預(yù)測結(jié)果（見表6）可以得出2025年研究區(qū)內(nèi)各類用地面積與占比，多項(xiàng)數(shù)據(jù)與2015年天津市土地利用數(shù)據(jù)比較可得：除城鎮(zhèn)用地、農(nóng)村居民點(diǎn)、其他建設(shè)用地和水域有不同程度的增長，其他用地類型出現(xiàn)了各種程度的負(fù)向變化.城鎮(zhèn)用地變化量為433.94 km2，單一土地利用動(dòng)態(tài)度為3.76%，是增加面積最多的用地類型，若不考慮其他政策因素，天津的城鎮(zhèn)化發(fā)展還將會(huì)繼續(xù)快速推進(jìn)；農(nóng)村居民點(diǎn)用地面積增長98.56 km2，鄉(xiāng)村居民點(diǎn)出現(xiàn)融合集中發(fā)展趨勢.根據(jù)研究可知，若不考慮政策影響，2015—2025年天津還將進(jìn)一步擴(kuò)大城鄉(xiāng)建設(shè)用地.預(yù)測得2015—2025年內(nèi)，耕地面積轉(zhuǎn)出程度較高，這將增加天津耕地資源壓力，可能帶來環(huán)境惡化等影響，因?yàn)榻ㄔO(shè)用地具有不可逆性，政府應(yīng)適當(dāng)進(jìn)行合理規(guī)劃防止天津市城鄉(xiāng)建設(shè)用地?zé)o序發(fā)展.

表6 2015—2025年天津市土地利用變化情況

2015—2025年天津市土地利用預(yù)測狀況大致上與1995—2015年天津市變化情況一致，說明CAMarkov模型能夠較好地預(yù)測土地利用變化趨勢.

4結(jié)論

本文基于CA-Markov耦合模型對(duì)天津市2015年土地利用進(jìn)行模擬，且拓展kappa系數(shù)達(dá)到良好標(biāo)準(zhǔn)，可以確定轉(zhuǎn)移矩陣和用地適宜性圖集約束控制和轉(zhuǎn)換規(guī)則，以預(yù)測2025年天津市土地利用情況.

（1）1995—2015年間天津市城鎮(zhèn)化發(fā)展迅速，土地利用變化顯著，城鎮(zhèn)用地面積從533.55 km2增加到1 153.35 km2，增長率由50.65%上升到116.16%；農(nóng)村居民點(diǎn)用地從871.62 km2增加到1 089.05 km2，天津城市化發(fā)展迅速；其他用地均有不同程度減少，變化劇烈程度呈現(xiàn)“水域＞草地＞林地＞耕地”的特征.天津城市擴(kuò)展方式為中心城區(qū)著重發(fā)展，市郊區(qū)以中心城區(qū)為原點(diǎn)放射性發(fā)展，城市化進(jìn)程速度較快，其空間擴(kuò)張形式為攤餅式發(fā)展為主，放射式低密度蔓延.

（2）運(yùn)用2015年實(shí)際數(shù)據(jù)與預(yù)測數(shù)據(jù)進(jìn)行對(duì)比，得出CA-Markov模型模擬程度較好，利用擴(kuò)展的kappa指數(shù)進(jìn)行精度檢驗(yàn)，可以從多個(gè)評(píng)價(jià)維度全面分析模擬結(jié)果的優(yōu)劣，降低潛在誤差對(duì)結(jié)果的影響.2015年天津市用地模擬與實(shí)際2015年解譯土地利用情況較為一致，預(yù)測結(jié)果良好，可以用于后續(xù)的2025年土地利用預(yù)測.

（3）預(yù)測2015—2025年間城鄉(xiāng)建設(shè)用地與耕地之間轉(zhuǎn)化量較大，耕地、林地、草地用地的不斷減少將會(huì)帶來水土流失，生態(tài)失調(diào).隨著國家生態(tài)文明改革進(jìn)程的不斷推進(jìn)，城市的發(fā)展應(yīng)該更加關(guān)注人為活動(dòng)對(duì)自然空間產(chǎn)生的影響[15].從天津市近二十年土地利用變化情況分析，隨著中心城區(qū)的城市化發(fā)展，以及濱海新區(qū)填海造陸等工程，城鄉(xiāng)居住以及建設(shè)用地大幅度增加，城鄉(xiāng)流動(dòng)明顯.而由土地利用變化圖和CA-Markov預(yù)測可得，2025年農(nóng)用地將明顯減少，這種變化有違基本農(nóng)田的保護(hù)政策.

天津市留存定量的農(nóng)業(yè)用地和農(nóng)村景觀不僅僅是由于生態(tài)保護(hù)紅線的控制，也應(yīng)該是研究區(qū)城鄉(xiāng)融合發(fā)展的路徑，無論從生態(tài)、社會(huì)、經(jīng)濟(jì)多個(gè)角度，農(nóng)村都應(yīng)該是需要保留合理引導(dǎo)發(fā)展的區(qū)域.科學(xué)劃定“生態(tài)保護(hù)紅線、永久基本農(nóng)田、城鎮(zhèn)開發(fā)邊界”三條控制勢在必行[16].科學(xué)系統(tǒng)地劃定三條控制線，能夠有效地調(diào)整好“三生”空間之間的格局關(guān)系.生態(tài)保護(hù)紅線應(yīng)將把生態(tài)功能極重要、生態(tài)極敏感脆弱的區(qū)域、經(jīng)過評(píng)估雖然不能確定但還是具有潛在生態(tài)價(jià)值的自然保護(hù)地區(qū)域納入生態(tài)紅線范圍內(nèi).城鎮(zhèn)開發(fā)邊界應(yīng)以城鎮(zhèn)開發(fā)現(xiàn)狀為基礎(chǔ)，進(jìn)行資源環(huán)境承載力的評(píng)價(jià)、國土開發(fā)適應(yīng)性評(píng)價(jià)等，以此為依據(jù)，進(jìn)行城鎮(zhèn)發(fā)展定位等研究，提出主要發(fā)展方向和空間布局[17].永久基本農(nóng)田的劃定應(yīng)在嚴(yán)守耕地紅線的基礎(chǔ)上，將達(dá)到質(zhì)量要求的耕地依法劃入.

未來規(guī)劃研究應(yīng)秉持環(huán)境就是生產(chǎn)力的綠色發(fā)展觀念，從始至終將保護(hù)自然環(huán)境作為經(jīng)濟(jì)發(fā)展內(nèi)在因素和內(nèi)生動(dòng)力，處理好城市化發(fā)展與生態(tài)安全之間的關(guān)系[18].當(dāng)城鎮(zhèn)化水平達(dá)到一定程度時(shí)，粗放的發(fā)展方式將會(huì)轉(zhuǎn)變?yōu)榇媪堪l(fā)展發(fā)揮資源的綜合效能的發(fā)展模式[19].但是本研究中沒有將經(jīng)濟(jì)政策等因素考慮在內(nèi)，對(duì)于天津市2015—2025年的土地利用的模擬結(jié)果是基于天津市整體的土地利用變化趨勢相對(duì)不變的前提下，進(jìn)行可能的預(yù)測，得出的結(jié)論有待于在今后的實(shí)踐中進(jìn)一步考量.