雷 雪, 王倩娜, 鄧靖琳, 高黃根, 沈 一*

(1. 四川大學(xué) 建筑與環(huán)境學(xué)院, 四川 成都 610065; 2. 四川省國土空間規(guī)劃研究院, 四川 成都 610081)

土地資源是國土空間資源的核心要素,土地利用的結(jié)構(gòu)和組成隨著城市發(fā)展發(fā)生動態(tài)變化,對區(qū)域生態(tài)環(huán)境帶來重要影響.不合理土地利用模式易導(dǎo)致土壤退化、水土流失、生物多樣性下降等生態(tài)問題[1-2].在國土空間規(guī)劃中,對土地資源可持續(xù)利用,是實(shí)現(xiàn)山水林田湖草保護(hù)與城鎮(zhèn)發(fā)展平衡的關(guān)鍵[3].

目前,在土地利用時空變化研究領(lǐng)域,常用的方法包括轉(zhuǎn)移矩陣模型[4]、CLUE-S模型[5]、CA-Markov模型[6]、FLUS模型[7]、重力模型[8]等,其中FLUS模型的模擬精度高于CLUE-S等其他常用模型[9].譚遠(yuǎn)模等[10]利用轉(zhuǎn)移矩陣、土地利用重心等方法分析廣州市土地利用發(fā)展的時空變化,發(fā)現(xiàn)城市發(fā)展有不平衡趨勢;郭麗英等[11]利用重力模型研究環(huán)渤海區(qū)域耕地重心變化情況,反映了城市化背景下耕地銳減和空間萎縮等實(shí)際問題.鹿瑤[12]利用FLUS模型對徐州市2020年土地利用進(jìn)行模擬,計(jì)算該市生態(tài)用地變化量,提出生態(tài)優(yōu)化方法,并對比模擬結(jié)果證明了優(yōu)化方法的有效性;趙林峰等[13]通過FLUS模型預(yù)測珠三角地區(qū)2025—2045年土地利用情況,對未來城市擴(kuò)張作出了預(yù)警.總體而言,FLUS模型相對成熟,能夠更高精度地挖掘出城市化背景下生態(tài)用地的變化特征.

眉山市作為四川省“成德眉資一體化”戰(zhàn)略部署城市之一,位于成渝地區(qū)雙城經(jīng)濟(jì)圈發(fā)展主軸,是成都極核的重要組成部分,具有“千湖之城,百園之市”[14]的美譽(yù).近年來,隨著城市高速擴(kuò)張,眉山市土地利用呈現(xiàn)出保護(hù)開發(fā)矛盾明顯,空間格局協(xié)調(diào)不足的問題.眉山市國土空間開發(fā)適宜性評價結(jié)果顯示,洪雅縣為生態(tài)保護(hù)極重要區(qū),青神縣、丹棱縣為農(nóng)業(yè)發(fā)展適宜區(qū),而在省級國土空間功能區(qū)規(guī)劃中,東坡區(qū)、彭山區(qū)、青神縣、仁壽縣、丹棱縣主體功能定位為城鎮(zhèn)發(fā)展區(qū),洪雅縣主體功能定位為農(nóng)產(chǎn)品主產(chǎn)區(qū),與眉山市各區(qū)縣的土地利用情況匹配度不高.與此同時,生態(tài)空間、城鎮(zhèn)空間、農(nóng)業(yè)空間存在沖突,城鎮(zhèn)用地和生態(tài)紅線、永久基本農(nóng)田紅線范圍有所重疊,導(dǎo)致生態(tài)退化、基本農(nóng)田保有量減少.因此,需重視對眉山市土地利用的研究,權(quán)衡各土地利用類型的關(guān)系,以實(shí)現(xiàn)生態(tài)、農(nóng)業(yè)、城鎮(zhèn)空間的平衡發(fā)展.封煒[15]利用轉(zhuǎn)移矩陣模型分析了1990—2016年眉山市土地利用變化情況及其驅(qū)動因子,但是目前眉山市相關(guān)研究仍相對缺乏.本研究以眉山市為例,結(jié)合“十四五”規(guī)劃中于2035年基本實(shí)現(xiàn)“建設(shè)人與自然和諧共生的現(xiàn)代化”遠(yuǎn)景目標(biāo),選定2035年為時間節(jié)點(diǎn),利用轉(zhuǎn)移矩陣模型、FLUS模型、重心模型,分析2005—2018年土地利用變化規(guī)律,預(yù)測2035年土地利用模式及重心變化趨勢,定量與定性分析相結(jié)合,判斷眉山市未來土地利用情況,為眉山市國土空間規(guī)劃與城市可持續(xù)發(fā)展提供理論依據(jù).

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

眉山市位于成都盆地西南部,東經(jīng)(102.49°～104.30°)、北緯(29.30°～30.16°)[16],全市土地總面積為7 138 km2,下設(shè)東坡區(qū)、彭山區(qū)、仁壽縣、洪雅縣、丹棱縣、青神縣6個行政區(qū)縣(圖1).眉山市屬長江流域,岷江、沱江水系,過境水資源充沛.南部以山區(qū)為主,地勢陡峭,海拔均在1 000 m以上;北部低山、丘陵、地勢緩平;中部為開闊的岷江河谷平原.

圖1 眉山市行政區(qū)劃分及自然環(huán)境概況

本研究采用2005、2010、2015、2018年共4期研究區(qū)域Landsat衛(wèi)星遙感影像數(shù)據(jù)[17]、眉山市DEM高程數(shù)據(jù)[18]和2018年眉山市基礎(chǔ)地理數(shù)據(jù)(由眉山市自然資源局提供).基于ArcGIS 10.2軟件對研究區(qū)域進(jìn)行圖像鑲嵌、行政區(qū)域裁剪等預(yù)處理,生成4期眉山市土地利用基礎(chǔ)柵格數(shù)據(jù)(30 m×30 m).參考我國《城市用地分類與規(guī)劃建設(shè)用地標(biāo)準(zhǔn)(GB50137—2011)》[19],將眉山市土地利用數(shù)據(jù)分為耕地、林地、草地、水域、建設(shè)用地、未利用地6個大類.

2 研究方法

2.1 土地利用轉(zhuǎn)移矩陣土地利用轉(zhuǎn)移矩陣能夠定量反映不同時期各地類的結(jié)構(gòu)以及轉(zhuǎn)移變化情況[4],其表達(dá)式為

S

(1)

其中,Sij為i類土地轉(zhuǎn)換成j類土地的面積,n為土地利用類型的數(shù)量.基于ArcGIS 10.2軟件,利用Dissolve、Intersect工具,結(jié)合Excel軟件,生成并統(tǒng)計(jì)了2005—2010、2010—2015、2015—2018共3期土地利用轉(zhuǎn)移數(shù)據(jù),以分析眉山市土地利用數(shù)量及方向特征.

2.2 單一土地利用動態(tài)度單一土地利用動態(tài)度表示一定時間內(nèi)某種特定土地利用類型的數(shù)量變化情況,能夠明確該地類的動態(tài)變化速度,其計(jì)算公式[20]為

(2)

其中,K為該地類的單一土地利用動態(tài)度,U為該地類的起始面積,U為該地類變化后的面積,T為研究時段.本次研究以2005年為起始值,2018年為末值,計(jì)算各地類面積的變化情況.

2.3 FLUS時空模擬模型FLUS模型是一種新型土地利用模擬模型[12],能夠預(yù)測并反映研究區(qū)域未來土地利用情景.模型的優(yōu)勢在于通過人工神經(jīng)網(wǎng)絡(luò)(ANN)訓(xùn)練模塊,導(dǎo)入社會、自然等空間自相關(guān)因子,在CA模塊中引入自適應(yīng)慣性競爭機(jī)制,有效處理各地類相互轉(zhuǎn)化的不確定性,主要包括以下2個模塊:

1) 基于神經(jīng)網(wǎng)絡(luò)的適宜性概率計(jì)算.神經(jīng)網(wǎng)絡(luò)算法(ANN)主要由輸入層、隱藏層、輸出層組成[12],其表達(dá)式為

(3)

其中,p(p,k,t)為第k種土地在柵格p、時間t上的適宜性概率,wj,k是隱含層與輸出層的權(quán)值,sigmoid是隱含層到輸出層的激勵函數(shù),netj(p,t)是第j個隱含層?xùn)鸥駊在時間t上所接收的信號.

2) 自適應(yīng)慣性競爭機(jī)制.FLUS模型運(yùn)用自適應(yīng)慣性機(jī)制模擬各地類的競爭問題,其核心是自適應(yīng)慣性系數(shù),其表達(dá)式[12]為

Interiatk=

TPtp,k=Pp,k×Ωtp,k×

Interiatk×(1-SCc→k),

(5)

土地利用變化是自然和社會因子長期作用的結(jié)果,自然環(huán)境本質(zhì)上決定了土地利用的類型、社會因子對土地利用變化產(chǎn)生重要影響[21].研究基于FLUS V2.2軟件,以2015年土地利用數(shù)據(jù)為基礎(chǔ),利用眉山市DEM高程數(shù)據(jù)、眉山市基礎(chǔ)地理數(shù)據(jù),生成自然驅(qū)動因子(高程、坡度、離水體距離)以及社會驅(qū)動因子(離農(nóng)村居民點(diǎn)距離、離城鎮(zhèn)距離、離公路距離、離鐵路距離),計(jì)算土地變化適宜性概率.借鑒已有鄰域權(quán)重參數(shù)成果[22],結(jié)合眉山市土地利用現(xiàn)狀,設(shè)置鄰域因子參數(shù)(表1),生成2018年眉山市土地利用模擬結(jié)果,并使用“Precision Validation”模塊對土地利用模擬結(jié)果進(jìn)行精度檢驗(yàn).同理,生成2035年土地利用模擬結(jié)果.

表1 土地利用鄰域因子參數(shù)

2.4 重心模型土地利用重心模型能夠反映各地類重心的空間轉(zhuǎn)移情況,體現(xiàn)土地利用的空間聚集性和位移趨勢,其表達(dá)式[11]為

x

(6)

y

(7)

(6)和(7)式中,Tij(i=1,2,3,…,n)表示第i評價單元第j年某要素的總量,Pi(xi,yi)為該評價單元的地理中心坐標(biāo),Pj(xj,yj)為該要素第j年的重心坐標(biāo).

基于2005、2018、2035年土地利用數(shù)據(jù),研究眉山市“過去-現(xiàn)在-未來”的土地利用重心轉(zhuǎn)移情況.利用ArcGIS 10.2中Mean center工具,基于CGCS 2000坐標(biāo)系統(tǒng),計(jì)算不同時期各地類的重心位置,為重心變化軌跡分析提供數(shù)據(jù)基礎(chǔ).

3 結(jié)果與分析

3.1 土地利用變化基于眉山市土地利用數(shù)據(jù),生成了2005—2018年眉山市土地利用圖(圖2),其中耕地為主要用地類型,林地主要集中在西部瓦屋山、總崗山區(qū)域以及中部彭祖山、龍泉山區(qū)域,水域面積占比較小,草地主要分布于西南區(qū)域,建設(shè)用地集中分布于各區(qū)縣的核心區(qū)域.

圖2 2005、2010、2015、2018年土地利用圖

3.1.1土地利用變化數(shù)量特征分析 基于轉(zhuǎn)移矩陣模型,生成2005—2018年3期土地利用轉(zhuǎn)移矩陣數(shù)據(jù)(表2).結(jié)果表明:2005—2010年,水域面積逐漸增加,后期基本穩(wěn)定在130 km2;林地面積呈現(xiàn)波動變化,2005—2010年林地退化29 km2,隨著近年來生態(tài)建設(shè)的加強(qiáng),林地面積恢復(fù)23.08 km2;2010—2015年草地面積增長28.97 km2,其余時段均在減少,說明草地生態(tài)系統(tǒng)不穩(wěn)定,易遭到侵蝕;耕地面積逐年遞減,其中2010—2015退化情況變緩;建設(shè)用地面積逐年遞增,2015—2018年增長率最大,高達(dá)27.78%,說明隨著天府新區(qū)的設(shè)立,眉山市城鎮(zhèn)化速度達(dá)到新的高點(diǎn).

表2 2005—2018年眉山市土地利用轉(zhuǎn)移矩陣數(shù)據(jù)

表2(續(xù))

2005—2018年眉山市土地利用總體變化情況表明(表3):草地減少幅度大、動態(tài)度高(-0.68%),為生態(tài)用地中退化最嚴(yán)重的用地類型;林地、水域動態(tài)度分別為0.13%、0.76%,呈增長趨勢,說明近年來眉山市生態(tài)建設(shè)已取得一定成效;建設(shè)用地增長幅度大,動態(tài)度最高(3.3%),為眉山市主要發(fā)展的土地類型.耕地減少幅度大但年動態(tài)度低(-0.21%),說明耕地面積萎縮嚴(yán)重,但由于總量較大,整體保留較為完整;未利用地變化幅度小但動態(tài)度高(-1.44%),可能成為眉山市最早消失的土地類型.

表3 2005—2018年眉山市土地利用總體變化情況

3.1.2土地利用變化方向特征分析 基于土地利用轉(zhuǎn)移矩陣,計(jì)算各地類相互轉(zhuǎn)化的方向及面積占比,提取出占比大于1%的轉(zhuǎn)化類型進(jìn)一步分析.圖3和表4表明,2005—2018年,林地、耕地相互轉(zhuǎn)化面積占比均在25%左右,說明退耕還林(主要在龍泉山中部、北部區(qū)域)和開墾占林(青衣江流域)的情況同時存在.耕地轉(zhuǎn)化為建設(shè)用地的面積占比約為15%,說明城鎮(zhèn)擴(kuò)張主要以侵蝕耕地為主,東坡區(qū)中部、彭山區(qū)北部、仁壽縣中北部以及與天府新區(qū)接壤的視高區(qū)域耕地退化明顯.草地轉(zhuǎn)化為林地的面積占比從10.2%下降到3.56%,林地轉(zhuǎn)化為草地面積占比從3.64%上升到9.24%,說明草地、林地邊界不穩(wěn)定,瓦屋山南部邊緣地區(qū)林草置換一直處于動態(tài)變化中,且草地受侵蝕程度更高.耕地、水域相互轉(zhuǎn)化面積占比均在2%以上,說明水域面積增加以向耕地擴(kuò)張為主,主要集中在瓦屋山水庫區(qū)域,農(nóng)田占用水資源的情況也同時存在,主要分布在青衣江下流流域.林地轉(zhuǎn)化為建設(shè)用地的面積逐年遞增,2015—2018年占比增加到1.09%,說明城市向林地擴(kuò)張的趨勢日益明顯,需加強(qiáng)對林地資源的保護(hù)與管控.

圖3 2005—2018年主要土地利用轉(zhuǎn)移類型

表4 2005—2018年眉山市主要土地利用轉(zhuǎn)化類型面積統(tǒng)計(jì)情況

3.2 基于FLUS模型的2035年土地利用模擬將基于FLUS模型的2018年眉山市土地利用模擬數(shù)據(jù)(圖4)與實(shí)際數(shù)據(jù)對比(表5),結(jié)果顯示各地類模擬面積相對誤差均小于1%,Kappa系數(shù)為0.896 5(大于0.75),精度較高,具有可信度,對比結(jié)果表明FLUS模型能夠用于模擬眉山市2035年土地利用結(jié)構(gòu).

圖4 2018年土地利用模擬結(jié)果

表5 2018年土地利用面積模擬誤差檢驗(yàn)

基于FLUS模型,導(dǎo)入7類土地利用自相關(guān)因子(圖5),生成了眉山市2035年土地利用模擬結(jié)果(表6和圖6).表6表明,2018—2035年眉山市建設(shè)用地、水域面積呈增長趨勢,草地面積將趨于穩(wěn)定,耕地、林地、未利用地均呈減少趨勢.其中,林地將減少31.4 km2,變化率為-1.9%,為唯一減少的生態(tài)用地類型,需增強(qiáng)對林地的保護(hù);建設(shè)用地將增加198.43 km2,增長率高達(dá)69.09%,其中彭山區(qū)、東坡區(qū)、青神縣城鎮(zhèn)擴(kuò)張明顯,其建設(shè)用地沿岷江河岸在南北方向上形成連接,若不加以人工干預(yù),2035年岷江沿岸建設(shè)用地將呈現(xiàn)出一體化趨勢,或?qū)︶航鷳B(tài)帶造成一定破壞,需在保護(hù)重要生態(tài)節(jié)點(diǎn)的基礎(chǔ)上,科學(xué)引導(dǎo)岷江兩岸城鎮(zhèn)發(fā)展方向;耕地將減少178.86 km2,變化率為-3.83%,耕地退化進(jìn)一步嚴(yán)重,需加強(qiáng)對耕地的管控力度.

(a) DEM

表6 2035年眉山市土地利用模擬結(jié)果

圖6 眉山市2035年土地利用模擬結(jié)果

3.3 土地利用重心轉(zhuǎn)移分析基于重心模型計(jì)算主要用地類型的重心位置如表7和表8所示,均采用CGCS 2000坐標(biāo)系統(tǒng),生成了2005—2035年眉山市土地利用重心轉(zhuǎn)移軌跡圖(圖7).

圖7 眉山市主要土地類型重心變化情況

表7 2005—2035年眉山市主要用地類型重心位置變化情況

表8 2005—2035年眉山市主要用地類型重心位置移動情況

結(jié)果表明,2005—2018年,草地重心向東北方向移動6.49 km,水域重心向西移動1.4 km,林地重心向東北移動2.13 km,建設(shè)用地、耕地重心無明顯變化.結(jié)合2005—2018年眉山市土地變化情況,發(fā)現(xiàn)土地重心變化與土地利用面積變化具有一定耦合性,其中瓦屋山南部大量草地退化、瓦屋山西部水庫面積大幅度增加、東北部龍泉山林地面積顯著提升分別與草地向北部丹棱縣蔓延、水域向西部延伸、林地向東北部發(fā)展的情況相對應(yīng),說明不均衡的土地面積變化是土地重心移動的重要原因之一.

2018—2035年,土地利用重心轉(zhuǎn)移速度較上一階段有所加快.其中,草地重心將向西南方向移動27.83 km,說明在面積無明顯變化的情況下,草地重心移動受到氣候條件、人類活動、地形地貌的影響,草地更趨向于往高海拔、遠(yuǎn)離城鎮(zhèn)、溫度更低的區(qū)域發(fā)展.林地重心將向西移動4.51 km,說明受到岷江限制、中東部城市建設(shè)影響,林地整體傾向于往人為干擾較低的西部山地區(qū)域延伸.水域重心將向西南方向移動7.18 km,說明水域趨于向西部生態(tài)環(huán)境良好的地區(qū)漫延,西部水域面積(以青衣江、瓦屋山水庫、槽魚灘水庫為主)有增加趨勢,應(yīng)注重彈性濕地空間建設(shè).建設(shè)用地重心將向東南方向移動1.39 km,說明受相關(guān)政策、地形地貌影響,建設(shè)用地整體傾向于往有利于經(jīng)濟(jì)發(fā)展、水源充足的岷江東部平原地區(qū)延伸.耕地重心將向西移動16.81 km,從仁壽縣西部轉(zhuǎn)移到東坡區(qū)中部,說明耕地受岷江流域建設(shè)用地擴(kuò)張影響,整體向西部淺山地區(qū)發(fā)展.

4 結(jié)論及討論

4.1 結(jié)論以眉山市為研究對象,采用“土地利用轉(zhuǎn)移矩陣-FLUS模擬-重心轉(zhuǎn)移”綜合方法,在分析2005—2018年土地利用變化規(guī)律的基礎(chǔ)上,基于FLUS模型,對2035年的土地利用情況進(jìn)行模擬預(yù)測,最后明確2005—2035年各地類的重心轉(zhuǎn)移特征.研究結(jié)果表明:

1) 2005—2018年,眉山市土地利用變化動態(tài)度表明:建設(shè)用地(3.3%)>水域(0.76%)>林地(0.13%)>耕地(-0.21%)>草地(-0.68%)>未利用地(1.44%),其中建設(shè)用地增長最為迅速,主要侵蝕中北部地區(qū)耕地以實(shí)現(xiàn)城鎮(zhèn)規(guī)模擴(kuò)張;水域增長平穩(wěn),主要由耕地轉(zhuǎn)入,其中瓦屋山水庫為近年來增加面積最大的水源地;草地、林地生態(tài)穩(wěn)定性低,西南部邊界區(qū)域草地階段性衰退趨勢明顯,林地總體呈增長趨勢,主要由龍泉山退耕還林以及瓦屋山草地轉(zhuǎn)化而來;耕地退化最為嚴(yán)重,未利用地逐年減少.

2) FLUS模型模擬檢驗(yàn)結(jié)果顯示,Kappa系數(shù)為0.896 5,相對誤差小于1%,證明FLUS模型能夠較準(zhǔn)確地模擬出未來土地利用情況.基于FLUS模型的2035年土地利用模擬結(jié)果顯示,建設(shè)用地將增加69.09%,其中岷江沿岸建設(shè)用地擴(kuò)張明顯;林地將減少31.4 km2,成為唯一減少的生態(tài)用地類型;耕地將退化178.76 km2,成為減少最多的土地利用類型;97.31%的未利用地轉(zhuǎn)化為其他用地.

3) 眉山市地處四川盆地邊緣地帶,擁有低山、丘陵、平原、山地等多種地形地貌,土地利用空間移動方向受地形地貌因子影響,在市域空間分布上具備移動的規(guī)律性:2005—2035年,林地、水域、草地總體傾向于向眉山市西部、西南部山地區(qū)域轉(zhuǎn)移,耕地呈現(xiàn)出向西北淺丘地區(qū)轉(zhuǎn)移的趨勢,建設(shè)用地是唯一向岷江東側(cè)平原地區(qū)偏移的土地利用類型.不難發(fā)現(xiàn),城鎮(zhèn)用地趨于向水源充足、地勢平坦、開發(fā)難度小的河谷平原擴(kuò)張;生態(tài)用地(林地、草地、水域)傾向于向高海拔、自然生境完整、人為干擾小的山區(qū)延展;農(nóng)業(yè)用地主要受城鎮(zhèn)擴(kuò)張影響,西北部向低山、丘陵地區(qū)發(fā)展.

4.2 討論1) 眉山市土地利用變化總體呈現(xiàn)出建設(shè)用地高速擴(kuò)張、農(nóng)業(yè)用地大量萎縮、生態(tài)用地不穩(wěn)定波動的規(guī)律.城鎮(zhèn)過度擴(kuò)張將導(dǎo)致生態(tài)空間退化、生產(chǎn)空間縮減,引發(fā)的自然災(zāi)害、糧食安全等問題將反向制約城市發(fā)展,因此需權(quán)衡各土地利用類型的關(guān)系,合理分配土地資源,促進(jìn)優(yōu)勢地區(qū)重點(diǎn)發(fā)展,加強(qiáng)生態(tài)功能區(qū)重點(diǎn)保護(hù).

根據(jù)眉山市土地利用分析及模擬結(jié)果,可將眉山市劃分為生態(tài)保護(hù)區(qū)、農(nóng)田保護(hù)區(qū)以及城鎮(zhèn)發(fā)展區(qū).其中,生態(tài)保護(hù)區(qū)主要包含洪雅縣南部瓦屋山區(qū)域、丹棱縣西部總崗山區(qū)域以及仁壽縣西部龍泉山區(qū)域;農(nóng)業(yè)保護(hù)區(qū)主要包括洪雅縣北部、丹棱縣中東部、青神縣、仁壽縣西部的丘陵、低山區(qū)域;城鎮(zhèn)發(fā)展區(qū)主要包括東坡區(qū)、彭山區(qū)、視高區(qū)域、各縣中心城區(qū)及周邊的平原地區(qū).具體邊界范圍需結(jié)合國土空間規(guī)劃中的“雙評估”“雙評價”結(jié)果,遵循全域“全覆蓋、不交叉、不重疊”的原則進(jìn)行進(jìn)一步劃定.各區(qū)應(yīng)突出土地利用的主導(dǎo)功能,其中生態(tài)保護(hù)區(qū)內(nèi)需減少不必要的人為活動,強(qiáng)化生態(tài)保育和生態(tài)建設(shè),限制開發(fā)性、生產(chǎn)性建設(shè)活動;農(nóng)業(yè)保護(hù)區(qū)內(nèi)需加強(qiáng)基本農(nóng)田保護(hù),因地制宜發(fā)展特色農(nóng)業(yè),防止建設(shè)用地對耕地進(jìn)一步侵蝕;城鎮(zhèn)發(fā)展區(qū)需科學(xué)引導(dǎo)城鎮(zhèn)用地擴(kuò)張方向,重點(diǎn)優(yōu)化岷江沿岸城市用地布局,落實(shí)最嚴(yán)格的節(jié)約用地制度,保護(hù)區(qū)域內(nèi)部生態(tài)用地,最大程度減少對農(nóng)業(yè)用地的侵蝕.土地利用的可持續(xù)性是實(shí)現(xiàn)生態(tài)環(huán)境與城市發(fā)展平衡的關(guān)鍵,眉山市作為“成德眉資一體化”戰(zhàn)略部署城市之一,具有一定典型意義,研究成果可以為其他西部中小型城市的土地利用時空變遷分析及國土空間規(guī)劃等提供參考.

2) FLUS模型在傳統(tǒng)CA模型上引入人工神經(jīng)網(wǎng)絡(luò)(ANN)訓(xùn)練模塊,使模擬精度進(jìn)一步提高.由于土地利用變化驅(qū)動因子錯綜復(fù)雜,本次研究只選取了最直接的7個自然社會因子進(jìn)行適宜性概率計(jì)算,對經(jīng)濟(jì)、氣候、文化等因子的考慮不夠充分,對模擬結(jié)果準(zhǔn)確性造成一定影響.與此同時,本次研究基于2005—2018年土地利用變化規(guī)律對鄰域參數(shù)進(jìn)行提取,具有一定主觀性.因此,在之后的研究中將對FLUS模型的因子選取與參數(shù)設(shè)置領(lǐng)域進(jìn)行更深入的探索.

致謝四川大學(xué)橫向科技項(xiàng)目(19H0949)對本文給予了資助，謹(jǐn)致謝意.