山東省縣域土地利用多功能時(shí)空演變及影響因素

孫進(jìn)，張全景，王世清

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

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

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

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

1.1 研究區(qū)概況

山東省地處我國(guó)東部沿海地區(qū)，位于黃河下游、京杭大運(yùn)河的中北段處，介于114°47′—112°42.3′E、34°22′—38°24′N(xiāo) 之間，全省國(guó)土總面積約15.71 萬(wàn)km2（圖1）。境內(nèi)地形地貌差異明顯，東部的山東半島主要為起伏和緩的波狀丘陵區(qū)；西部和北部多為黃河沖積平原，地勢(shì)相對(duì)低洼平坦；中部則主要為魯中南山地丘陵區(qū)。山東省地處北溫帶半濕潤(rùn)季風(fēng)氣候區(qū)，四季分明、雨熱同期、水系發(fā)達(dá)，境內(nèi)河流分屬淮河、黃河、海河三大流域，山東半島獨(dú)流入海水系。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ù)來(lái)源

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

2 研究方法

2.1 土地利用多功能測(cè)度方法

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

表1 山東省土地利用多功能性評(píng)價(jià)指標(biāo)體系

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

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

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

第一步，計(jì)算指標(biāo)標(biāo)準(zhǔn)化值Zij的均值Ej：

第二步，計(jì)算指標(biāo)標(biāo)準(zhǔn)化值Zij的均方差值σj：

第三步，計(jì)算各指標(biāo)權(quán)重Wj：

土地利用功能測(cè)度：本文運(yùn)用線(xiàn)性加權(quán)法計(jì)算了各準(zhǔn)則層和目標(biāo)層的土地利用功能值，計(jì)算公式為：

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

2.2 地理探測(cè)器

地理探測(cè)器是探測(cè)空間分異性，揭示其背后驅(qū)動(dòng)因子的一種統(tǒng)計(jì)學(xué)方法［21］。本文運(yùn)用地理探測(cè)器中的因子探測(cè)模塊分別探測(cè)了各功能指標(biāo)對(duì)土地利用多功能時(shí)空分異的解釋程度，用q值表征，表達(dá)式為：

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

3 結(jié)果及分析

3.1 土地利用總功能時(shí)空演變特征

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

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

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

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

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

3.2 土地利用單項(xiàng)功能時(shí)空演變特征

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

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

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

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

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

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

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

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

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

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

4 結(jié)論與建議

4.1 結(jié)論

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

4.2 建議

本文基于以上研究結(jié)論，綜合考慮山東省新舊動(dòng)能轉(zhuǎn)換發(fā)展背景和土地利用現(xiàn)狀，提出以下發(fā)展建議：①魯西南和魯西北平原地區(qū)應(yīng)借助現(xiàn)代科技手段和經(jīng)營(yíng)方式提升農(nóng)業(yè)生產(chǎn)功能，積極轉(zhuǎn)變經(jīng)濟(jì)發(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)勢(shì)打造優(yōu)勢(shì)產(chǎn)業(yè)，為經(jīng)濟(jì)發(fā)展增速。同時(shí)，完善林地分級(jí)管護(hù)措施，增強(qiáng)森林資源的環(huán)境效益。③魯東地區(qū)應(yīng)立足區(qū)位和海洋資源優(yōu)勢(shì)，打造海洋生產(chǎn)和生態(tài)保護(hù)模式，加快新舊動(dòng)能轉(zhuǎn)換，實(shí)現(xiàn)經(jīng)濟(jì)增長(zhǎng)與生態(tài)保護(hù)并舉；優(yōu)化基礎(chǔ)設(shè)施建設(shè)，進(jìn)一步提升居民生活質(zhì)量，打造沿海宜居城市帶。