孫進，張全景，王世清

（1.曲阜師范大學地理與旅游學院，山東日照 276826；2.日照市國土空間規(guī)劃與生態(tài)建設(shè)重點實驗室，山東日照 276826）

近年來，隨著城市化與工業(yè)化進程的持續(xù)推進，我國土地利用方式和土地利用格局變化顯著，土地資源的稀缺性日益顯現(xiàn)［1］，各種土地利用類型之間的爭奪與沖突日趨激烈［2］，不合理的用地結(jié)構(gòu)導致資源環(huán)境沖突形勢也變得越來越嚴峻［3，4］。隨著黃河流域生態(tài)保護和高質(zhì)量發(fā)展上升為國家重大戰(zhàn)略，以及我國首個新舊動能轉(zhuǎn)換綜合試驗區(qū)的建設(shè)推進，山東省的創(chuàng)新發(fā)展、綠色發(fā)展和生態(tài)保護得到高度重視，土地利用的功能多樣化迎來良好發(fā)展契機?；谕恋刭Y源狀況、國情需要、國家戰(zhàn)略等多方形勢需要的背景下，土地利用功能也呈現(xiàn)出多樣化趨勢。因此，開展土地利用多功能研究能夠為當前形勢下科學指導土地管理、協(xié)調(diào)各類用地的關(guān)系提供理論支持和經(jīng)驗借鑒。

土地利用多功能性的概念最早源自國際上對農(nóng)業(yè)、生態(tài)系統(tǒng)服務(wù)和景觀功能的研究。2004 年啟動的全球土地計劃（GLP）支持的歐盟第六框架計劃下的SENSOR 項目提出了土地利用多功能性的概念，并應(yīng)用于區(qū)域可持續(xù)性評價［5］。此后，國外學者聚焦于土地利用多功能性綜合評價［6］，探究不同土地利用政策對土地利用多功能的影響［7，8］。國內(nèi)學者從土地利用多功能性的概念界定與分析框架［9］、識別與分類［10－12］、可持續(xù)性評價［13］、時空演變特征分析［14，15］、權(quán)衡與協(xié)同［16－18］、影響因素探究［19］等方面展開。多數(shù)研究在評價尺度上以國家、省、市為主，而以縣域為評價單元進行區(qū)域時空演變特征的研究則相對較少。本文以山東省為例，綜合運用GIS 空間分析、地理探測器等方法模型分析了2005—2018年縣域尺度上土地利用多功能性時空演變特征及影響因素，為土地資源優(yōu)化配置和可持續(xù)利用提供科學參考。

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

1.1 研究區(qū)概況

山東省地處我國東部沿海地區(qū)，位于黃河下游、京杭大運河的中北段處，介于114°47′—112°42.3′E、34°22′—38°24′N 之間，全省國土總面積約15.71 萬km2（圖1）。境內(nèi)地形地貌差異明顯，東部的山東半島主要為起伏和緩的波狀丘陵區(qū)；西部和北部多為黃河沖積平原，地勢相對低洼平坦；中部則主要為魯中南山地丘陵區(qū)。山東省地處北溫帶半濕潤季風氣候區(qū)，四季分明、雨熱同期、水系發(fā)達，境內(nèi)河流分屬淮河、黃河、海河三大流域，山東半島獨流入海水系。2020年，山東省地區(qū)生產(chǎn)總值73129.0 億元，三次產(chǎn)業(yè)結(jié)構(gòu)調(diào)整為7.3 ∶39.1 ∶53.6。

圖1 山東省研究區(qū)行政區(qū)劃

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

數(shù)據(jù)來源：①社會經(jīng)濟數(shù)據(jù)主要來源于《中國縣域統(tǒng)計年鑒》（2006—2019 年）、《中國區(qū)域經(jīng)濟統(tǒng)計年鑒》（2011 年）、《山東統(tǒng)計年鑒》（2006—2019 年）、各縣域統(tǒng)計年鑒與統(tǒng)計公報，個別年份存在統(tǒng)計數(shù)據(jù)缺失，則按照年均增長率推算得出。②水域面積、林地面積、區(qū)域總面積等數(shù)據(jù)來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心（http：//www.resdc.cn/）對應(yīng)年份的山東省土地利用現(xiàn)狀圖（2005 年、2010 年、2018年）。

2 研究方法

2.1 土地利用多功能測度方法

指標體系構(gòu)建：本文借鑒國土空間“三生”功能理論［20］，從山東省農(nóng)業(yè)種植、社會經(jīng)濟發(fā)展、生態(tài)環(huán)境保護現(xiàn)狀等實際情況出發(fā)，將土地利用功能區(qū)分為生產(chǎn)功能、生活功能、生態(tài)功能3 類，并基于以上3 個維度，選取相應(yīng)指標構(gòu)建山東省土地利用多功能性評價指標體系（表1）。

表1 山東省土地利用多功能性評價指標體系

指標權(quán)重確定：首先，運用極差標準化方法對原始數(shù)據(jù)進行去量綱化處理。計算公式為：

式中，Xij為第i年的第j 項指標值；maxXj、minXj分別為研究期內(nèi)第j 項指標的最大值和最小值；Zij為第i年第j項指標的標準化值。

其次，運用均方差決策法確定指標權(quán)重。

第一步，計算指標標準化值Zij的均值Ej：

第二步，計算指標標準化值Zij的均方差值σj：

第三步，計算各指標權(quán)重Wj：

土地利用功能測度：本文運用線性加權(quán)法計算了各準則層和目標層的土地利用功能值，計算公式為：

式中，F(xiàn)i表示評價單元i的土地利用功能值。

2.2 地理探測器

地理探測器是探測空間分異性，揭示其背后驅(qū)動因子的一種統(tǒng)計學方法［21］。本文運用地理探測器中的因子探測模塊分別探測了各功能指標對土地利用多功能時空分異的解釋程度，用q值表征，表達式為：

式中，h為對因變量或自變量的分級數(shù)；L 為各分級區(qū)域的個數(shù)；N和Nh分別為全區(qū)和第h 級的單元樣本數(shù)；和σ2分別為第h 級和全區(qū)的土地利用總功能值的方差；q的分布區(qū)間為［0，1］。

3 結(jié)果及分析

3.1 土地利用總功能時空演變特征

采用地統(tǒng)計分析（Geostatistical Analyst）模塊中的趨勢分析工具（Trend Analysis Tool）分析山東省縣域土地利用多功能值時空分異的總體格局特征及趨勢，結(jié)果見圖2。從圖2 可見，2005—2018 年山東省縣域東西方向總體上地域差異不顯著，呈由西至東土地利用多功能水平緩慢上升的分布格局，空間形態(tài)呈現(xiàn)為一條坡度變化平緩的曲線，且差距逐漸縮小。研究區(qū)域南北方向總體呈現(xiàn)中部高、南北低的倒“U”型空間格局，地域差異較東西方向顯著。

圖2 2005—2018 年山東省土地利用多功能水平趨勢研究

從圖3 可見，2005—2010 年山東省土地利用總功值明顯上升，高水平和較高水平分布區(qū)明顯增多，低水平分布區(qū)則明顯減少。其中，高水平分布區(qū)由2005 年的泰安中心城區(qū)、張店區(qū)、青島中心城區(qū)3個增至2010 年的濟南中心城區(qū)、煙臺中心城區(qū)、青島中心城區(qū)等7 個區(qū)縣；處于較高水平的分布區(qū)由2005年的濟南中心城區(qū)、滕州市、博山區(qū)等12 個區(qū)縣增至2010 年的龍口市、榮成市、平度市等25 個區(qū)縣。與2005 年相比，魯西南的單縣、菏澤中心城區(qū)、鄆城縣等，魯西北的陽谷縣、冠縣、臨清市等和魯中地區(qū)的昌樂縣、安丘市、沂水縣等眾多區(qū)縣由低水平功能區(qū)轉(zhuǎn)變?yōu)檩^低水平及中等水平分布區(qū)，魯西北的武城縣、平原縣、陵城區(qū)等、魯東的蓬萊區(qū)、文登區(qū)、黃島市、日照中心城區(qū)等區(qū)縣由較低水平轉(zhuǎn)為中等水平和較高水平。

圖3 2005—2018 年山東省土地利用總功能時空分異

從圖3 可以發(fā)現(xiàn)，2018 年山東省土地利用總功值與2010 年相比有了明顯提高，高水平與較高水平分布區(qū)明顯增多，較低水平與低水平分布區(qū)明顯減少。其中，2010—2018 年山東省土地利用多功能性高水平分布區(qū)由點狀分布轉(zhuǎn)變?yōu)槊鏍罘植?，形成了四大面狀高水平集聚區(qū)，分別包括以泰安中心城區(qū)與濟南中心城區(qū)為核心的岱岳區(qū)、章丘區(qū)、萊蕪區(qū)等9個區(qū)縣，以張店為核心的臨淄區(qū)、周村區(qū)、桓臺縣等7 個區(qū)縣，榮成市、威海中心城區(qū)、文登區(qū)3 個區(qū)縣，龍口市、蓬萊區(qū)、煙臺中心城區(qū)3 個區(qū)縣。另外，還包括滕州市、棗莊中心城區(qū)、城陽區(qū)、青島中心城區(qū)、黃島區(qū)等較為零散的高水平點狀分布區(qū)。2018年山東省土地利用多功能性較高水平分布區(qū)與2010年相比，片狀集聚分布明顯，魯東地區(qū)形成了由東北向西南延伸呈面狀分布的較高水平集聚區(qū)。此外，在魯西南的菏澤中心城區(qū)、鄆城縣、鄄城縣等區(qū)縣，魯西北的武城縣、寧津縣、陵城區(qū)、平原縣等區(qū)縣、魯中的寧陽縣、曲阜市、兗州區(qū)等區(qū)縣也形成了較小范圍的較高水平集聚區(qū)。2018 年較低水平與低水平分布區(qū)與2010 年相比則明顯減少。其中，2018年低水平分布區(qū)僅剩余微山縣，較低水平分布區(qū)僅剩余無棣縣、沾化區(qū)、夏津縣、山亭縣等15 個區(qū)縣，其他均轉(zhuǎn)變?yōu)檩^低水平和中等水平分布區(qū)。中等水平分布區(qū)呈現(xiàn)出減少的趨勢，其中新增中等水平分布區(qū)主要來自于冠縣、曹縣、郯城區(qū)、樂陵市等低水平、較低水平分布區(qū)，減少的中等水平分布區(qū)主要提升至蓬萊區(qū)、乳山市、牟平區(qū)、定陶區(qū)等高水平及較高水平分布區(qū)。

3.2 土地利用單項功能時空演變特征

生產(chǎn)功能時空演變特征：從圖4 可見，2005—2010年較高水平分布區(qū)及低水平分布區(qū)變化較大，與2005 年相比，2010 年山東省土地生產(chǎn)功能低水平分布區(qū)僅剩余淄川區(qū)和博山區(qū)兩個區(qū)縣。2018 年高水平分布區(qū)數(shù)量與2010 年相比未發(fā)生變化，其中新泰由較高水平轉(zhuǎn)為高水平分布區(qū)，臨淄區(qū)由高水平分布區(qū)轉(zhuǎn)為較高水平分布區(qū)。2005—2010 年山東省土地生產(chǎn)功能中等水平分布區(qū)轉(zhuǎn)入與轉(zhuǎn)出并存，新增中等水平分布區(qū)主要來源于慶云縣、沾化區(qū)、墾利區(qū)等較低水平、低水平分布區(qū)，減少的中等分布區(qū)主要轉(zhuǎn)化為郯城區(qū)、乳山市、單縣等較高水平與高水平分布區(qū)；2010—2018 年高水平分布區(qū)由7 個增加到9 個，其中，棲霞市、青島中心城區(qū)、蘭陵縣、陽谷縣等區(qū)縣由較高水平分布區(qū)轉(zhuǎn)變?yōu)楦咚椒植紖^(qū)，肥城市由高水平分布區(qū)轉(zhuǎn)變?yōu)檩^高水平分布區(qū)。2018年低水平分布區(qū)的區(qū)縣數(shù)目由2010 年的2 個增加至3 個，沂源縣由中等水平轉(zhuǎn)化為低水平。2010—2018 年，較高水平、中等水平與較低水平分布區(qū)轉(zhuǎn)入與轉(zhuǎn)出并存，與高水平和低水平分布區(qū)相比變化較大。

圖4 2005—2018 年山東省土地生產(chǎn)功能時空分異

生活功能時空演變特征：從圖5 可見，2005—2010年山東省土地生活功能時空變化差異顯著，以低水平與較低水平、低水平與中等水平之間的變化為主。高水平和較高水平呈零星分布，其中，2010年高水平分布區(qū)由2005 年的5 個轉(zhuǎn)化為2010 年的15 個，主要來自于濟南中心城區(qū)、濱州中心城區(qū)、臨沂中心城區(qū)等較高水平、中等水平分布區(qū)。2010 年低水平分布區(qū)與2005 年相比大面積減少，主要轉(zhuǎn)化為日照中心城區(qū)、莒南縣、冠縣等較低水平、中等水平分布區(qū)。2005—2010 年中等水平分布區(qū)顯著增加，這一時期中等水平分布區(qū)的轉(zhuǎn)入與轉(zhuǎn)出并存，新增中等水平分布區(qū)主要來自于低水平、較低水平分布區(qū)，減少的中等水平分布區(qū)主要轉(zhuǎn)為了較高水平、高水平分布區(qū)；2010—2018 年山東省土地生活功能較高水平分布區(qū)顯著增加，呈片狀分布，形成了魯西南定陶區(qū)、曹縣、巨野縣等區(qū)縣，魯中鋼城區(qū)、萊蕪區(qū)、廣饒縣、壽光市等區(qū)縣，魯東日照中心城區(qū)、諸城市、莒縣等區(qū)縣，魯西北武城縣、寧津縣、平原縣等區(qū)縣四大面狀分布區(qū)。相較于2000 年2018 年低水平分布區(qū)明顯減少，僅剩余沂源縣一個區(qū)縣，與2010年相比，2018 年較低水平、中等水平分布區(qū)均呈現(xiàn)減少的趨勢，向高水平、較高水平轉(zhuǎn)化明顯。2018年，山東省生活功能高水平分布區(qū)明顯增多且形成了小范圍集聚的面狀分布。

圖5 2005—2018 年山東省土地生活功能時空分異

生態(tài)功能時空演變特征：從圖6 可見，2005—2010年山東省土地生態(tài)功能與其他功能區(qū)相比變化較小，高水平分布區(qū)由2005 年的2 個增加至2010年的4 個，淄川區(qū)和威海中心城區(qū)由較高水平轉(zhuǎn)為高水平分布區(qū)。2010 年較高水平分布區(qū)由2005 年的12 個增加至14 個，其中，淄川區(qū)、威海中心城區(qū)由較高水平轉(zhuǎn)為高水平分布區(qū)，東營中心城區(qū)、昌邑市、文登區(qū)、臨朐縣等區(qū)縣由中等水平轉(zhuǎn)化為較高水平分布區(qū)。2010 年低水平、較低水平和中等水平分布區(qū)轉(zhuǎn)入與轉(zhuǎn)出并存，與2005 年相比數(shù)目變化較?。?010—2018 年山東省土地生態(tài)功能與2005—2010 年相比時空分異變化較為顯著。2018 年高水平分布區(qū)顯著增多，由2010 年的4 個增加至9 個，且在魯東地區(qū)的煙臺中心城區(qū)、牟平區(qū)、威海中心城區(qū)出現(xiàn)了小范圍的面狀分布。2010—2018 年低水平分布區(qū)顯著減少，僅剩余慶云縣、金鄉(xiāng)縣、滕州市、嶧城區(qū)等4 個區(qū)縣。2018 年較高水平分布區(qū)在魯中地區(qū)出現(xiàn)面狀分布，費縣、平邑縣、岱岳區(qū)等較低水平和中等水平分布區(qū)轉(zhuǎn)化為較高水平分布區(qū)。2018年較低水平和中等水平分布區(qū)與2010 年相比均有所增加，主要來源于單縣、陽谷縣、萊陽市等低水平分布區(qū)。

圖6 2005—2018 年山東省土地生態(tài)功能時空分異

3.3 山東省土地利用多功能時空演變的驅(qū)動因素

本文根據(jù)山東省縣區(qū)尺度土地利用單項功能及指標，并基于ArcGIS10.2 平臺，將該省2005 年、2010年、2018 年指標值利用自然斷裂法劃分為7 個等級，運用地理探測器測算各指標的q 值，具體如表2所示。

表2 山東省土地利用多功能性影響因子探測

由表2 可知，研究期內(nèi)經(jīng)濟發(fā)展水平、空間承載力、社會保障能力等對山東省土地利用多功能性空間分異的解釋力度較強，且呈逐步增強的趨勢，而農(nóng)業(yè)生產(chǎn)水平、環(huán)境凈化能力、資源供給能力對該省土地利用多功能性空間分異的解釋力度較弱，且在2005—2018 年間呈現(xiàn)減弱的趨勢。主要表現(xiàn)在：①經(jīng)濟發(fā)展水平與土地利用多功能的空間分布關(guān)系密切。經(jīng)濟密度（X4）的解釋力在2005 年、2010 年、2018年分別為0.558、0.594、0.640。經(jīng)濟密度代表單位面積上經(jīng)濟活動的效率和土地利用的密集程度，表征著一個地區(qū)高質(zhì)量發(fā)展水平的高低。經(jīng)濟密度越大，說明能夠以更低的土地資源消耗支撐更高效益的經(jīng)濟發(fā)展。經(jīng)濟增長是經(jīng)濟密度演化的宏觀動力，經(jīng)濟的發(fā)展使得居民收入水平和消費需求擴大刺激產(chǎn)業(yè)加速發(fā)展，從而推動區(qū)域內(nèi)土地經(jīng)濟功能提升，促進土地利用多功能水平的提高。第二、三產(chǎn)業(yè)占比（X5）在3 個時間截面的解釋力分別為0.515、0.512、0.406，產(chǎn)業(yè)結(jié)構(gòu)演化是一個長期的歷史過程，在其不同發(fā)展歷史階段，經(jīng)濟增長的動力也隨之演變，促使土地利用由單一的農(nóng)業(yè)生產(chǎn)轉(zhuǎn)向多產(chǎn)業(yè)融合發(fā)展模式，通過調(diào)整優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，刺激土地利用多功能潛力的發(fā)揮，有效提升土地利用效率。②社會保障能力對于土地利用多功能性的空間分異也具有顯著影響。萬人醫(yī)療衛(wèi)生床位數(shù)（X8）在3 個時間截面上的解釋力分別為0.555、0.627、0.699，隨著社會經(jīng)濟發(fā)展、人民生活水平改善，基本醫(yī)療衛(wèi)生服務(wù)和居民健康保障水平不斷上升，醫(yī)療設(shè)施建設(shè)成為城鄉(xiāng)社會服務(wù)體系的重要組成，然而優(yōu)質(zhì)醫(yī)療資源在山東省各區(qū)域的分布不夠均衡。完善欠發(fā)達地區(qū)醫(yī)療基礎(chǔ)設(shè)施建設(shè)，對于提升區(qū)域公共衛(wèi)生服務(wù)水平，促使土地更有效地發(fā)揮生活功能，提高人類福祉具有重要意義。人均城鄉(xiāng)居民儲蓄存款余額（X9）的解釋力分別為0.489、0.601、0.595，其增加一定程度上受傳統(tǒng)習慣和居民對于未來收入的信心的影響，但關(guān)鍵在于社會基礎(chǔ)保障制度的完善和居民收入水平的提升，人們面對突發(fā)事件時的抗風險能力增強。應(yīng)指出的是，人均城鄉(xiāng)居民儲蓄存款余額的持續(xù)增長從側(cè)面反映了人們投資意向薄弱、投資渠道閉塞，引導人們選擇較低風險的投資理財積累財富，這也是提高社會保障功能的關(guān)鍵要義。③就業(yè)支撐能力對于土地利用多功能性的空間分異解釋力度較強。從業(yè)人口密度（X10）在3 個時間截面上的解釋力分別為0.647、0.650、0.699，其大小能夠反映所在區(qū)域經(jīng)濟社會發(fā)展程度和就業(yè)服務(wù)水平。近年來，山東省經(jīng)濟運行穩(wěn)中有進，經(jīng)濟結(jié)構(gòu)不斷優(yōu)化，對穩(wěn)定基礎(chǔ)性崗位、拓寬就業(yè)創(chuàng)業(yè)渠道、提高就業(yè)質(zhì)量形成了有力拉動，形成了保障民生和增強經(jīng)濟發(fā)展活力的良性循環(huán)。

4 結(jié)論與建議

4.1 結(jié)論

本文綜合運用GIS 空間分析、地理探測器模型分析了2005—2018 年山東省縣域土地利用多功能時空演變及影響因素，主要結(jié)論如下：①運用趨勢分析、GIS空間分析方法發(fā)現(xiàn)，2005—2018 年山東省縣域土地利用多功能水平總體呈現(xiàn)“東高西低、南北低中部高”的空間格局。2005—2018 年山東省縣域土地利用多功能差距逐漸縮小，空間分布格局差異明顯，具體表現(xiàn)為高水平和較高水平分布區(qū)顯著增加，低水平和較低水平分布區(qū)顯著減少，中等水平分布區(qū)在這一時期轉(zhuǎn)入與轉(zhuǎn)出并存，總體變化相對較小。2005—2018 年山東省縣域土地生產(chǎn)功能時空變化較小，生活功能與生態(tài)功能變化顯著。②運用地理探測器分析土地利用多功能時空演變的影響因素，發(fā)現(xiàn)2005—2018 年山東省經(jīng)濟發(fā)展水平、社會保障能力、就業(yè)支撐能力對山東省縣域土地利用多功能的時空分異影響最大，經(jīng)濟密度、人口密度、建設(shè)用地比例、萬人醫(yī)療衛(wèi)生床位數(shù)的解釋力較強且對山東省土地利用多功能時空分異的影響逐漸增強；而水果產(chǎn)量，第二、三產(chǎn)業(yè)占比，人均水域面積，森林覆蓋率等探測因素對研究區(qū)土地利用多功能空間分異的影響相對較弱。

4.2 建議

本文基于以上研究結(jié)論，綜合考慮山東省新舊動能轉(zhuǎn)換發(fā)展背景和土地利用現(xiàn)狀，提出以下發(fā)展建議：①魯西南和魯西北平原地區(qū)應(yīng)借助現(xiàn)代科技手段和經(jīng)營方式提升農(nóng)業(yè)生產(chǎn)功能，積極轉(zhuǎn)變經(jīng)濟發(fā)展方式，創(chuàng)新型發(fā)展生態(tài)農(nóng)業(yè)，并加大鹽堿地改造、退耕還林、水污染治理等生態(tài)建設(shè)力度，提高生態(tài)功能。②魯中南山地丘陵區(qū)應(yīng)優(yōu)化城鄉(xiāng)居民生產(chǎn)生活布局，保障就業(yè)、醫(yī)療、住房、交通等公共服務(wù)水平，完善相關(guān)基礎(chǔ)設(shè)施建設(shè)，充分利用資源優(yōu)勢打造優(yōu)勢產(chǎn)業(yè)，為經(jīng)濟發(fā)展增速。同時，完善林地分級管護措施，增強森林資源的環(huán)境效益。③魯東地區(qū)應(yīng)立足區(qū)位和海洋資源優(yōu)勢，打造海洋生產(chǎn)和生態(tài)保護模式，加快新舊動能轉(zhuǎn)換，實現(xiàn)經(jīng)濟增長與生態(tài)保護并舉；優(yōu)化基礎(chǔ)設(shè)施建設(shè)，進一步提升居民生活質(zhì)量，打造沿海宜居城市帶。