摘 要：隨著經(jīng)濟(jì)的快速發(fā)展和城鎮(zhèn)化進(jìn)程的不斷推進(jìn)，肥城市的土地利用格局發(fā)生了較大變化。精準(zhǔn)把握肥城市土地利用變化特征，對優(yōu)化土地利用空間格局具有重要意義。基于2018—2020年肥城市的Landsat 8 OLI_TIRS遙感影像數(shù)據(jù)，采用土地利用動(dòng)態(tài)度、土地利用類型轉(zhuǎn)移矩陣、重心遷移模型等方法對肥城市土地利用變化狀況進(jìn)行分析。結(jié)果表明：2018—2020年，肥城市土地利用類型以耕地、城鄉(xiāng)建設(shè)用地、林地為主；2018—2020年，肥城市土地利用類型變化、時(shí)空變化相對較快，綜合土地利用動(dòng)態(tài)度為8.73%，其中土地利用動(dòng)態(tài)度最大的地類是城鄉(xiāng)建設(shè)用地，最小的是草地；2018—2020年，肥城市城鄉(xiāng)建設(shè)用地、林地、水域、未利用地面積擴(kuò)大，而耕地、草地面積縮小。

關(guān)鍵詞：GIS；土地利用變化；肥城市

中圖分類號：F321.1；P237 文獻(xiàn)標(biāo)志碼：A 文章編號：1674-7909（2024）5-122-5

DOI：10.19345/j.cnki.1674-7909.2024.05.028

0 引言

土地利用是指人類開展的所有有目的的開發(fā)利用土地資源的行為［1-2］。分析土地利用的時(shí)空變化是各領(lǐng)域的重要研究內(nèi)容，能為土地資源合理利用提供重要的科學(xué)依據(jù)［3］。目前，諸多學(xué)者已對土地利用變化進(jìn)行了研究。陳伊多等［4］基于各期Landsat TM/ETM遙感影像，采用土地利用動(dòng)態(tài)度、轉(zhuǎn)化趨勢等方法分析了西藏自治區(qū)土地利用/覆被特征。馮應(yīng)斌等［5］綜合運(yùn)用空間自相關(guān)模型、地理探測器等，分析了貴州省耕地顯性邊際化的總體態(tài)勢、空間集聚效應(yīng)及相關(guān)影響因素。廖美玉等［6］采用Mann-Kendall、滑動(dòng)t檢驗(yàn)、標(biāo)準(zhǔn)差橢圓和重心遷移模型等方法，探究了近40年黃河流域土地利用/覆被時(shí)空變化特征。張晨基等［7］在采用土地利用動(dòng)態(tài)度及土地利用轉(zhuǎn)移矩陣方法的基礎(chǔ)上，疊合DEM高程數(shù)據(jù)，對新疆土地利用時(shí)空變化進(jìn)行研究。安彬［8］采用土地利用程度綜合指數(shù)、土地利用信息圖譜等方法，研究了1985—2020年漢江生態(tài)經(jīng)濟(jì)帶土地利用/覆蓋演變特征。以上研究多從大的時(shí)間和空間尺度上對土地利用進(jìn)行分析，不利于捕捉小的時(shí)間和空間尺度上的土地利用時(shí)空變化?；诖耍P者以肥城市為研究區(qū)域，基于2018—2020年的土地利用狀況，探究該區(qū)域的時(shí)空變化特征，以期為肥城市土地利用空間規(guī)劃提供理論依據(jù)。

1 研究區(qū)概況

肥城市隸屬于山東省泰安市，地處山東中部、泰山西麓，位于北緯35°53′～36°19′，東經(jīng)116°28′～116°59′；地勢北高南低，由東北向西南傾斜；中部隆起地帶為丘陵山地，北部是康匯平原，南部是汶陽平原，自然形成了平原、山地、丘陵等多種地貌形態(tài)［9］；屬溫帶季風(fēng)氣候區(qū)，四季分明，雨熱同期，光照充足，氣候溫暖，年平均降水量為659 mm。肥城市礦產(chǎn)資源豐富，礦業(yè)勘探開發(fā)歷史悠久，探明資源量大。目前，肥城市已發(fā)現(xiàn)礦產(chǎn)29種，占山東省已發(fā)現(xiàn)礦產(chǎn)總數(shù)的47.6%。

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

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

2018—2020年肥城市土地利用遙感監(jiān)測數(shù)據(jù)來自地理空間數(shù)據(jù)云。各年份土地利用類型數(shù)據(jù)通過對各期Landsat 8 OLI_TIRS 遙感影像（分辨率為30 m）人工目視解譯生成，利用ENVI 5.6裁剪、輻射定標(biāo)、大氣校正等預(yù)處理，在ArcGIS 10.2和Excel中進(jìn)行處理統(tǒng)計(jì)。

2.2 研究方法

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

土地利用動(dòng)態(tài)度反映一定時(shí)期內(nèi)某一區(qū)域不同土地利用類型的數(shù)量變化，一般分為單一土地利用動(dòng)態(tài)度和綜合土地利用動(dòng)態(tài)度［10］。

2.2.1.1 單一土地利用動(dòng)態(tài)度

單一土地利用動(dòng)態(tài)度表示一定時(shí)期內(nèi)某一區(qū)域某一種土地利用類型的數(shù)量變化情況［11］，其表達(dá)式見式（1）。

式（1）中：Ki為研究區(qū)域內(nèi)第i類土地利用類型在a～b時(shí)段內(nèi)的土地利用動(dòng)態(tài)度；Uia、Uib分別為a、b時(shí)段第i類土地類型的面積；Ta、Tb分別為研究初期時(shí)間與研究末期時(shí)間。

Ki值越大，表示該土地利用類型轉(zhuǎn)化為其他類型土地的數(shù)量越多，即在研究時(shí)段內(nèi)相對變化程度較大；反之，轉(zhuǎn)化數(shù)量越少，在研究時(shí)段內(nèi)相對變化程度較小。

2.2.1.2 綜合土地利用動(dòng)態(tài)度

綜合土地利用動(dòng)態(tài)度表示一定時(shí)期內(nèi)某一區(qū)域不同土地利用類型之間的數(shù)量變化情況［12］，其表達(dá)式見式（2）。

式（2）中：LC為一定時(shí)期內(nèi)綜合土地利用動(dòng)態(tài)度；LUi為研究初期第i類土地利用類型的面積；ΔLUi為研究期內(nèi)第i類土地利用類型轉(zhuǎn)化為非i類土地利用類型面積的絕對值；T為研究時(shí)段。

LC值越大，表示該區(qū)域土地利用類型的變化程度越大；反之，表示該區(qū)域土地利用類型的變化程度越小。

2.2.2 土地利用類型轉(zhuǎn)移矩陣

土地利用類型轉(zhuǎn)移矩陣以矩陣的形式將2個(gè)不同時(shí)期土地利用類型之間互相轉(zhuǎn)換的數(shù)量關(guān)系展現(xiàn)出來，是具體反映土地利用變化的結(jié)構(gòu)特征和各類型之間轉(zhuǎn)移方向的重要方法［13］。其表達(dá)式見式（3）。

式（3）中：Sij表示從研究初期到研究末期土地類型由i轉(zhuǎn)化為j的土地面積；n表示土地利用類型的數(shù)量；i、j分別為研究初期、研究末期的土地利用類型。

2.2.3 重心遷移模型

重心遷移模型用于反映研究區(qū)域內(nèi)各土地利用類型在時(shí)空分布上的動(dòng)態(tài)變化，可以精確把握土地利用類型的空間遷移趨勢［14］。其表達(dá)式見式（4）。

式（4）中：Xt和Yt分別為t時(shí)期某種土地利用類型質(zhì)心的經(jīng)度和緯度坐標(biāo)；Cti為t時(shí)期某種土地利用類型第i個(gè)斑塊的面積；Xti和Yti分別為t時(shí)期土地利用類型第i個(gè)斑塊的幾何重心經(jīng)度和緯度坐標(biāo)。

3 結(jié)果分析

3.1 土地利用分類結(jié)果

研究利用ENVI 5.6與ArcGIS 10.2，根據(jù)土地利用分級系統(tǒng)，對2018年、2019年、2020年肥城市3期土地利用遙感影像進(jìn)行目視解譯，得到3個(gè)時(shí)間節(jié)點(diǎn)的土地利用類型分類格局。

肥城市土地利用類型統(tǒng)計(jì)情況見表1。由表1可知，在肥城市的土地利用構(gòu)成中，耕地、城鄉(xiāng)建設(shè)用地和林地是三大主導(dǎo)類型，這3種類型的土地總面積占據(jù)了全市土地總面積的90%以上。在所有土地利用類型中，耕地面積所占比重一直最大但呈逐年遞減的趨勢，城鄉(xiāng)建設(shè)用地面積逐年增加，林地面積整體呈上升趨勢，草地、水域及未利用地面積占比較小且波動(dòng)較小。

具體到空間分布上，肥城市的耕地面積最大，主要分布在肥城市北部的康匯平原和肥城市南部的汶陽平原，而在肥城市中部的丘陵地帶則有相對零星分布；城鄉(xiāng)建設(shè)用地在肥城市北部地區(qū)更為密集，而在肥城市南部較為分散；林地主要集中分布在肥城市北部山區(qū)和肥城市中部丘陵地帶；草地多分布在肥城市北部山麓和肥城市中部的丘陵地帶；水域則主要分布在肥城市南部的大汶河、肥城市中部的康王河及肥城市北部的礦區(qū)塌陷區(qū)；未利用地在肥城市極為少見，只在肥城市北部山區(qū)的某些區(qū)域有少量分布。

3.2 土地利用時(shí)間變化

2018—2020年肥城市土地利用動(dòng)態(tài)度見表2。由表2可知，2018—2020年，土地利用動(dòng)態(tài)度最大的是城鄉(xiāng)建設(shè)用地，最小的是草地，說明在此期間城鄉(xiāng)建設(shè)用地增加的幅度和草地減少的幅度最為明顯；耕地的土地利用動(dòng)態(tài)度小于０，表明耕地面積逐年減少，這種現(xiàn)象可能是由城市化進(jìn)程加快造成的；林地、水域、草地、未利用地的土地利用動(dòng)態(tài)度均從小于0到大于0，表明其面積均為先減少后增加，2019年后這4類土地利用類型面積增大，可能與人類工農(nóng)業(yè)活動(dòng)、出行減少有關(guān)；城鄉(xiāng)建設(shè)用地的土地利用動(dòng)態(tài)度均大于0，表明城鄉(xiāng)建設(shè)用地的土地利用面積持續(xù)增加，說明2018—2020年肥城市經(jīng)濟(jì)發(fā)展速度較快、人類活動(dòng)強(qiáng)度加劇。2018—2020年肥城市的綜合土地利用動(dòng)態(tài)度為8.73%，綜合土地利用動(dòng)態(tài)度大于2.50%，說明土地利用變化為快速變動(dòng)型［15］，整體土地利用類型數(shù)量變化相對較快。

3.3 土地利用空間變化

3.3.1 土地利用轉(zhuǎn)移矩陣

2018—2019年、2019—2020年肥城市土地利用類型轉(zhuǎn)移矩陣分別見表3和表4。由表3和表4可知，2018—2020年，耕地的流出面積最多，流出面積為184.81 km2，主要流向了城鄉(xiāng)建設(shè)用地和林地；林地的流出面積次之，流出面積為47.13 km2，主要流向了耕地。流入面積最多的是城鄉(xiāng)建設(shè)用地，流入面積為110.18 km2，主要來源于耕地和林地；林地的流入面積次之，流入面積為85.13 km2，主要來源于耕地和草地。

總體而言，肥城市土地利用變化表現(xiàn)為城鄉(xiāng)建設(shè)用地、林地和水域的面積增加，而耕地、草地和未利用地的面積減少。這一變化主要是因?yàn)槌鞘谢M(jìn)程加快，城鄉(xiāng)邊界不斷擴(kuò)展，城市化水平逐漸提高。同時(shí)，退耕還林政策的實(shí)施和礦區(qū)塌陷區(qū)的植被恢復(fù)，尤其是速生楊的種植，也是林地面積增加的重要原因；礦區(qū)塌陷形成的水體導(dǎo)致水域面積增加。這些變化反映了肥城市土地利用的動(dòng)態(tài)調(diào)整既是城市化發(fā)展的必然結(jié)果，也是生態(tài)恢復(fù)和環(huán)境治理政策實(shí)施的具體體現(xiàn)。

3.3.2 重心遷移

2018—2020年肥城市土地利用類型的重心變化情況見表5。由表5可知，肥城市耕地重心向肥城市東南方向遷移，城鄉(xiāng)建設(shè)用地重心向肥城市西北方向遷移，二者遷移方向相反，主要是因?yàn)槌鞘袛U(kuò)張占用大量耕地資源、消耗大量水資源，西北部地形平坦且有河流經(jīng)過，適宜城市建設(shè)，因此城市向西北方向遷移。

4 結(jié)束語

筆者根據(jù)地理空間數(shù)據(jù)云的遙感監(jiān)測數(shù)據(jù)，采用土地利用動(dòng)態(tài)度、土地利用類型轉(zhuǎn)移矩陣、重心遷移模型等方法，并借助ArcGIS10.2、ENVI5.6和Excel，對2018—2020年肥城市土地利用類型時(shí)空變化特征進(jìn)行研究，得出以下結(jié)論。

①2018—2020年，肥城市土地利用類型以耕地、城鄉(xiāng)建設(shè)用地、林地為主；肥城市的綜合土地利用動(dòng)態(tài)度為8.73%，整體土地利用類型數(shù)量變化相對較快，多數(shù)土地利用類型的單一動(dòng)態(tài)度數(shù)量變化相對劇烈，與城市快速擴(kuò)張有關(guān)。因此，政府部門應(yīng)重視土地利用空間規(guī)劃，合理規(guī)劃利用土地，確保城市可持續(xù)發(fā)展。

②2018—2020年，肥城市城鄉(xiāng)建設(shè)用地、林地、水域、未利用地面積擴(kuò)大，耕地、草地面積縮小；耕地重心向肥城市東南方向遷移，城鄉(xiāng)建設(shè)用地重心向肥城市西北方向遷移。因此，政府部門應(yīng)加強(qiáng)對非農(nóng)業(yè)用地侵占基本農(nóng)田的管理，減少耕地非農(nóng)化、非糧化，確保耕地面積合理穩(wěn)定。

此研究以縣級市為例，從多種土地利用類型的綜合利用到單一土地利用類型的綜合利用的變化，綜合研究了肥城市土地利用類型特征，彌補(bǔ)了從大的時(shí)間和空間尺度上研究的局限性與片面性。但由于數(shù)據(jù)的可獲取性有限，僅研究了部分現(xiàn)象，沒有具體分析肥城市土地利用類型變化的驅(qū)動(dòng)因子，并且直接利用遙感影像數(shù)據(jù)進(jìn)行監(jiān)督分類，土地利用分類精度可能會存在一定的誤差。因此，在未來的研究中，研究者應(yīng)完善對肥城市土地利用類型驅(qū)動(dòng)因子的分析，提高分類精度，以得到更精確的研究結(jié)果。

參考文獻(xiàn)：

［1］何英彬，陳佑啟.土地利用/覆蓋變化研究綜述［J］.中國農(nóng)業(yè)資源與區(qū)劃，2004（2）：61-65.

［2］陳佑啟，楊鵬.國際上土地利用/土地覆蓋變化研究的新進(jìn)展［J］.經(jīng)濟(jì)地理，2001（1）：95-100.

［3］陳小蘭，胡娉婷.基于GIS的東莞市土地利用時(shí)空變化規(guī)律和特征分析［J］.廣西城鎮(zhèn)建設(shè)，2021（7）：19-21.

［4］陳伊多，楊慶媛.西藏自治區(qū)土地利用/覆被變化時(shí)空演變特征及驅(qū)動(dòng)因素［J］.水土保持學(xué)報(bào)，2022，36（5）：173-180.

［5］馮應(yīng)斌，楊溪源.基于土地利用變化的貴州省耕地顯性邊際化時(shí)空演變及其影響因素［J］.水土保持通報(bào)，2023，43（6）：263-271.

［6］廖美玉，方秀琴，蔣心遠(yuǎn)，等.黃河流域近40年土地利用/覆被時(shí)空變化特征［J/OL］.水土保持學(xué)報(bào)，1-14［2024-04-19］.https：//doi.org/10.13870/j.cnki.stbcxb.2024.02.030.

［7］張晨基，張永福.基于DEM的新疆土地利用變化時(shí)空格局［J］.中南農(nóng)業(yè)科技，2024，45（1）：116-121，126.

［8］安彬.1985—2020年漢江生態(tài)經(jīng)濟(jì)帶土地利用/覆蓋演變特征分析［J］.安康學(xué)院學(xué)報(bào)，2023，35（6）：16-23.

［9］周攀清.魯西南地區(qū)村落喪葬習(xí)俗探析：以肥城市村落調(diào)查為例［J］.學(xué)理論，2014（15）：164-165.

［10］郭健斌，張英，張志偉，等.基于地理探測器的藏東南高山峽谷區(qū)土地利用變化及其驅(qū)動(dòng)機(jī)制：以西藏林芝市為例［J］.中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2023，28（4）：210-226.

［11］王秀蘭，包玉海.土地利用動(dòng)態(tài)變化研究方法探討［J］.地理科學(xué)進(jìn)展，1999（1）：83-89.

［12］陳伊多，楊慶媛.西藏自治區(qū)土地利用/覆被變化時(shí)空演變特征及驅(qū)動(dòng)因素［J］.水土保持學(xué)報(bào)，2022，36（5）：173-180.

［13］任斐鵬，江源，熊興，等.東江流域近20年土地利用變化的時(shí)空差異特征分析［J］.資源科學(xué)，2011，33（1）：143-152.

［14］張瑞，杜國明，張樹文.1986—2020年東北典型黑土區(qū)耕地資源時(shí)空變化及其驅(qū)動(dòng)因素［J］.資源科學(xué)，2023，45（5）：939-950.

［15］薛晉聰，陶志梅.天津市各區(qū)土地利用程度差異及動(dòng)態(tài)變化分析［J］.合作經(jīng)濟(jì)與科技，2019（8）：8-11.

作者簡介：石舒心（2003—），女，本科生，研究方向：GIS。