解林曉，段藝芳，肖超菲，李金龍

（聊城大學地理與環(huán)境學院，山東聊城 252000）

0 引言

隨著中國城鎮(zhèn)化和工業(yè)化進程的不斷加快，以及計劃生育政策的放開，人地矛盾被進一步激化，諸如土地利用結構失衡、資源配置不合理、環(huán)境污染與破壞、土地人均占有量少等問題日漸突顯[1-3]，科學分析區(qū)域土地利用生態(tài)-社會-經濟系統(tǒng)耦合協(xié)調關系及其時空演變規(guī)律，識別主要影響因素，對于協(xié)調區(qū)域人地關系、促進區(qū)域可持續(xù)發(fā)展具有重要意義。

以往對土地利用效益的研究多從城鎮(zhèn)空間尺度展開，并逐漸由單一城市向城市群及區(qū)域尺度轉變，現已涵蓋鄉(xiāng)村、縣域、城市、省域、流域和全國等諸多尺度，如王磊、李啟權、趙娜倩等分別對城鄉(xiāng)一體化進程中的鄉(xiāng)村、川中丘陵縣域、晉城市等的土地利用效益進行相關研究[4-6]，田俊峰、范樹平、Liu等則從東北三省、巢湖流域、全國尺度等更大的空間尺度開展研究[7-9]。已有研究多采用DEA數據包絡模型、綜合指數模型、耦合協(xié)調度模型等定量化數學模型開展土地利用效益及耦合關系研究[10-12]，此外還引入了GIS、遙感技術等進行綜合分析[13-14]；近年來，隨著生態(tài)效益研究的日漸成熟，研究內容已由早期的社會、經濟單項效益評價[15-16]轉向綜合效益評價[17-18]，并開始關注區(qū)域內部不同土地利用效益之間的協(xié)調關系及其空間分布特征[16-17]。

目前，對耦合協(xié)調關系的研究多以兩兩之間的靜態(tài)現狀研究為主，而針對社會、經濟、生態(tài)三者之間整體耦合協(xié)調演變趨勢以及各子系統(tǒng)間相對發(fā)展水平的分析均略顯不足，此外對影響土地利用效益提升的諸多因素的探討也相對薄弱?；诖?，本研究從生態(tài)、社會、經濟3個方面構建土地利用效益評價指標體系和生態(tài)-社會-經濟耦合協(xié)調度模型，引入相對發(fā)展度模型和障礙度模型，揭示了山東省土地利用生態(tài)-社會-經濟耦合協(xié)調度的時空演變特征及趨勢，并探討了耦合協(xié)調度的主要影響因素，以期了解山東省土地利用效益未來發(fā)展趨勢，為促進山東省區(qū)域可持續(xù)協(xié)調發(fā)展及因地制宜地優(yōu)化相關政策提供理論與數據支持。

1 研究區(qū)概況

山東省位于中國華東沿海地帶，同時北承京津冀地區(qū)，南臨長江三角洲，地理位置優(yōu)越，且境內水路、鐵路交通網發(fā)達，具備承接兩大城市群產業(yè)、技術轉移的便利（圖1）。截至2017年底，轄17個地級市，137個縣級行政區(qū)，區(qū)域總面積達1.58×105km2。山東省是東部沿海的人口大省、經濟大省，但區(qū)域間發(fā)展質量很不平衡，空間差異明顯，大致呈由東向西逐漸降低的梯度分異現象，靠近沿海的魯東地區(qū)吸引外資與社會投資的能力明顯優(yōu)于其他地區(qū)[19]。

圖1 研究區(qū)區(qū)位圖

2 研究方法

2.1 評價指標體系的構建

本研究借鑒已有研究成果[20-22]，基于系統(tǒng)性、代表性、科學性等原則，以山東省17地市為研究對象，綜合考慮土地利用現狀及相關數據的可獲取性，從土地利用的生態(tài)效益、社會效益、經濟效益3個方面共選取19項指標，構建山東省土地利用效益評價指標體系（表1）。

表1 山東省土地利用效益評價指標體系

2.2 耦合協(xié)調度模型

本研究采用極值法對指標數據進行標準化和無量綱化處理，并考慮到熵權法為客觀賦權法，能有效避免主觀賦權法所致的人為誤差，故采用熵權法確定各指標的權重[19]。在此基礎上，通過加權求和測算得到山東省17地市土地利用綜合效益指數T。同時，為了反映各系統(tǒng)之間的良性互動關系，以避免由于低生態(tài)水平與低社會經濟發(fā)展水平所引起的高耦合狀況，本研究引入耦合協(xié)調度的概念[16-18]，構建土地利用生態(tài)-社會-經濟耦合協(xié)調度模型。土地利用生態(tài)-社會-經濟耦合協(xié)調度計算見公式(1)，土地利用生態(tài)-社會-經濟系統(tǒng)的耦合協(xié)調系數計算見公式(2)，土地利用綜合效益指數計算見公式(3)，土地利用的生態(tài)、社會、經濟效益指數計算見公式(4)。

式中：D表示土地利用生態(tài)-社會-經濟耦合協(xié)調度，C表示土地利用生態(tài)-社會-經濟系統(tǒng)的耦合協(xié)調系數，T表示土地利用綜合效益指數，U1、U2、U3分別表示土地利用的生態(tài)、社會、經濟效益指數，xij’為樣本經標準化處理后的數據。其中α、β、γ均為待定系數，且α+β+γ=1，三系統(tǒng)追求協(xié)調發(fā)展，此處取值。參考已有研究成果[12,23]，本研究對土地利用生態(tài)-社會-經濟系統(tǒng)耦合協(xié)調等級及其所處階段進行相應劃分（表2），以便于分析其時空格局及演變趨勢。

表2 耦合協(xié)調度等級及階段劃分

2.3 相對發(fā)展度模型

為了更加直觀地表達土地利用生態(tài)、社會、經濟子系統(tǒng)間的相對發(fā)展程度，并據此分析各制約類型區(qū)的發(fā)展演變情況。本研究引入相對發(fā)展度模型[18,23]，對土地利用生態(tài)、社會、經濟子系統(tǒng)間的相對發(fā)展程度進行測度，計算公式見(5)、(6)、(7)。

式中，Ecm、Ecs、Ems分別表示土地利用生態(tài)-經濟、生態(tài)-社會和經濟-社會相對發(fā)展度，Um、Uc、Us分別代指經濟、生態(tài)、社會效益指數。參考已有相關研究成果[18,23]，確定Ecm、Ecs、Ems的取值區(qū)間及其對應的實際含義（表3）。

表3 相對發(fā)展度取值區(qū)間劃分

2.4 障礙度模型

為進一步細化制約山東省及各地市不同時期土地利用效益耦合協(xié)調程度的主要障礙因子，本研究通過引入因子貢獻度Fj、指標偏離度Vj和障礙度Mj等變量對指標因子的障礙度[20,24-25]進行測算。具體見公式(8)、(9)、(10)。

式中，Mj表示第j項指標的障礙度值，可根據Mj的大小判斷制約土地利用效益耦合協(xié)調程度的主要障礙因子；Fj為因子貢獻度，即各指標相對于總目標的權重，表示第j項指標對總目標的影響程度；rj表示第j項指標所屬第i個準則層的權重，rj取值均為1/3；wj表示第j項指標在各準則層中的權重；Vj表示指標偏離度，它是各項指標與土地利用效益總目標之間的差距；xij’表示單項指標的標準化值。

3 結果與分析

3.1 耦合協(xié)調度時空演變

3.1.1 時序演變趨勢 2008—2017年間，山東省土地利用生態(tài)-社會-經濟耦合協(xié)調性整體處于磨合階段，并呈現出持續(xù)上升趨勢，且各地市間耦合協(xié)調度均衡化水平逐步提高（表4）。具體來看，山東省土地利用生態(tài)-社會-經濟耦合協(xié)調度多年均值為0.63，總體處于磨合階段的初級協(xié)調狀態(tài)，且呈現出持續(xù)增長的趨勢，由2008年的勉強協(xié)調(0.54)轉變?yōu)?017年的中級協(xié)調(0.70)，可見，山東省土地利用生態(tài)、社會、經濟三子系統(tǒng)間良性互動程度逐漸增強。此外，2008—2011年耦合協(xié)調度水平增幅明顯較大，這與該時期山東省生態(tài)效益水平增長相對緩慢，而社會、經濟水平發(fā)展較為迅速，且三者總體趨于協(xié)調有關。期間，各地市耦合協(xié)調度差異性逐漸縮小，變異系數由2008年的0.126降至2017年的0.072。

表4 山東省各地市土地利用生態(tài)-社會-經濟耦合協(xié)調度

具體來看，東營市、威海市耦合協(xié)調度最高，均高于0.70，達到了協(xié)調階段的中級協(xié)調狀態(tài)，且耦合協(xié)調等級逐漸升高，但年均增速較慢，年均增長率不足2.5%，這主要與東營市、威海市生態(tài)、社會、經濟基礎較好，但受產業(yè)結構制約較為嚴重有關；棗莊市、臨沂市、聊城市和菏澤市水平較低，介于0.50～0.60之間，處于磨合階段的勉強協(xié)調狀態(tài)，亦呈持續(xù)增長態(tài)勢，其中，菏澤市耦合協(xié)調度增速最快，由2008年的0.34增至2017年的0.58，年均增長率達5.97%，這與菏澤市產業(yè)結構優(yōu)化程度及節(jié)能減排力度關系密切，但目前仍處于瀕臨協(xié)調的邊緣，應注意繼續(xù)推進新舊動能轉換，加快社會經濟發(fā)展；其余11個地級市則處于磨合階段的初級協(xié)調狀態(tài)，整體水平介于0.60～0.70之間。

3.1.2 空間格局演變 山東省土地利用生態(tài)-社會-經濟耦合協(xié)調度呈自魯東至魯中、魯西遞減的梯度分布特征，且2008—2017年間各地市耦合協(xié)調度等級明顯提高（圖2）。2008年，菏澤市、聊城市、棗莊市、臨沂市均處于輕度失調或瀕臨失調狀態(tài)，東營市和威海市均處于初級協(xié)調狀態(tài)，其余11地市均呈勉強協(xié)調狀態(tài)；較之2008年，2011年的菏澤市、聊城市、棗莊市和臨沂市均提升至勉強協(xié)調狀態(tài)，東營市率先達到中級協(xié)調狀態(tài)，萊蕪市、濟寧市和日照市變化不明顯，其余8地市則均由勉強協(xié)調提升至初級協(xié)調狀態(tài)。至2014年，菏澤市持續(xù)提升至勉強協(xié)調狀態(tài)，聊城市、濟寧市、棗莊市、臨沂市、萊蕪市和日照市提升至初級協(xié)調狀態(tài)，青島市、煙臺市和威海市提升至中級協(xié)調狀態(tài)，其余地市則無明顯變化。該階段的變化主要歸結于山東半島藍色經濟區(qū)、黃河三角洲高效生態(tài)經濟區(qū)、魯南經濟帶等重點區(qū)域發(fā)展戰(zhàn)略的有效實施，使得各地市生態(tài)、社會、經濟均有所提高；至2017年，濟南市和淄博市提升至中級協(xié)調狀態(tài)，而其余地市則無明顯變化。總體來看，山東省各地市生態(tài)、社會、經濟協(xié)同發(fā)展效果越來越明顯，但仍未出現良好協(xié)調區(qū)，且區(qū)域空間分異特征明顯，故山東省耦合協(xié)調度水平尚存在較大提升空間，區(qū)域一體化仍有待進一步加強。

圖2 山東省各地市耦合協(xié)調度空間演變

3.2 相對發(fā)展關系空間演變

從空間格局來看，山東省各地市經歷了由社會、經濟滯后向生態(tài)滯后及協(xié)調發(fā)展的轉變，其中，魯西北以生態(tài)滯后為主，魯中南部以社會滯后為主，半島地區(qū)逐漸向協(xié)調發(fā)展區(qū)轉變（圖3）。2008年，山東省主要受社會、經濟因素的制約，除聊城市、德州市、濱州市、東營市等，受生態(tài)建設步伐緩慢和社會基礎服務水平較低的影響，構成生態(tài)、社會滯后區(qū)外，其余地市生態(tài)建設力度大于社會經濟發(fā)展水平，屬于社會、經濟滯后區(qū)；至2011年，社會、經濟因素仍屬于主要制約因素，但在重點區(qū)域帶動戰(zhàn)略驅動下，社會基礎建設水平不斷提高，社會、經濟滯后區(qū)范圍明顯縮小，并出現了社會滯后區(qū)和經濟滯后區(qū)的分化，而東營市、濱州市積極加大社會基礎建設力度，但經濟建設力度有所減緩，屬于生態(tài)、經濟滯后區(qū)；2014年，受國家生態(tài)文明建設的推動，濟寧市、泰安市、煙臺市、威海市演變?yōu)閰f(xié)調發(fā)展區(qū)，而青島市及魯西北部分城市生態(tài)問題逐漸顯現，構成了明顯的生態(tài)滯后區(qū)，但因各地區(qū)受社會、經濟因素制約程度的不同，而分化為生態(tài)、社會滯后區(qū)和生態(tài)、經濟滯后區(qū)；至2017年，生態(tài)影響進一步顯現，威海市同魯中地區(qū)南部兩側及其以北構成了明顯的半圈式生態(tài)滯后區(qū)，且亦存在著生態(tài)、社會滯后區(qū)和生態(tài)、經濟滯后區(qū)的分化，而日照市得益于經濟的快速發(fā)展，該時期已演變?yōu)閰f(xié)調發(fā)展區(qū)。

圖3 2008—2017年山東省各制約類型區(qū)空間演變

3.3 障礙因子診斷分析

歷年障礙度排名前6位的障礙因子，障礙度累計貢獻率均超過50%，故選取歷年前6位障礙因子進行土地利用效益耦合協(xié)調度障礙因子分析（表5）。結果表明，2008—2017年間的障礙因子主要包括C1、C2等生態(tài)因子，C9、C11、C13等社會因子以及C15、C16、C17等經濟因子，其中，社會及經濟因子一直占據主導地位，但2013年后生態(tài)影響明顯增強。2008—2012年間，公路交通等社會基礎設施建設相對落后，社會經濟資產和生產技術雖持續(xù)積累，但受山東省雙核結構影響，邊緣城市社會經濟效益增長相對緩慢，加之耕地不合理的開發(fā)利用現象較為普遍且長期存在，使得該時期C19、C16、C15、C11、C9等社會、經濟障礙因子占據主導地位；2013—2017年間，受社會、經濟障礙因子影響的同時，C1、C2等生態(tài)障礙因子的影響程度不斷加強，該時期社會經濟發(fā)展已具備一定規(guī)模，但生態(tài)環(huán)境方面的負面影響逐步顯現，如林地面積減小，濕地生態(tài)環(huán)境惡化等，不利于區(qū)域協(xié)調可持續(xù)發(fā)展。

表5 山東省土地利用效益主要障礙因子及其障礙度

4 結論與建議

4.1 結論

本研究通過構建土地利用生態(tài)-社會-經濟耦合協(xié)調度模型、相對發(fā)展度模型和障礙度模型，深入分析了山東省17地市土地利用生態(tài)、社會和經濟3個子系統(tǒng)間的整體耦合協(xié)調性、相對發(fā)展狀態(tài)及主要障礙因子，可為山東省土地利用生態(tài)、社會、經濟間的協(xié)調可持續(xù)發(fā)展提供相應的理論支持。主要研究結論如下。

（1）2008—2017年，山東省土地利用生態(tài)-社會-經濟整體處于磨合階段的初級協(xié)調狀態(tài)，且耦合協(xié)調度呈持續(xù)上升趨勢，由勉強協(xié)調(0.54)轉變?yōu)橹屑墔f(xié)調(0.70)；各地市耦合協(xié)調度均衡化水平明顯提高，但區(qū)域間耦合協(xié)調水平仍存在一定差異?？臻g上，整體呈自魯東至魯中、魯西遞減的梯度分布特征，中部地市梯度差異逐漸縮小，并向初級或中級協(xié)調階段轉變，但尚有較大改善空間，仍未出現良好協(xié)調區(qū)，區(qū)域一體化仍有待進一步加強。

（2）各地市土地利用生態(tài)、社會和經濟子系統(tǒng)間相對發(fā)展度差異顯著。至2017年，除濟寧市、泰安市、日照市、煙臺市屬于協(xié)調發(fā)展區(qū)外，威海市同魯中地區(qū)南部兩側及北部構成了明顯的半圈式生態(tài)滯后區(qū)，且因受社會、經濟制約程度的不同而分化為生態(tài)、社會滯后區(qū)和生態(tài)、經濟滯后區(qū)。

（3）影響山東省土地利用生態(tài)、社會、經濟耦合協(xié)調發(fā)展的主要障礙因子為社會經濟因子，但2013年后生態(tài)影響有所增強。當前主要障礙因子為新增林地面積、人均公路里程、地均水資源總量、人均耕地面積、地均GDP、地均固定資產投資總額。

4.2 建議

（1）考慮到當前“生態(tài)影響有所增強”以及人均土地面積較少的現狀，未來應加強科技創(chuàng)新培養(yǎng)，大力推進高新技術產業(yè)發(fā)展，并嘗試與綠色生態(tài)產業(yè)相結合，不斷提高第三產業(yè)比重，積極倡導綠色循環(huán)經濟模式；合理規(guī)劃土地利用布局，進一步加強對林地和濕地的生態(tài)治理與保護。

（2）根據當前各地市土地利用生態(tài)、社會、經濟耦合協(xié)調性制約因素的差異化特點，可分別實施區(qū)域差別化發(fā)展對策。臨沂市、萊蕪市等經濟滯后區(qū)，應依托重點區(qū)域帶動戰(zhàn)略，不斷優(yōu)化營商環(huán)境，加快對外招商引資和新舊動能轉換，增強城市經濟實力和綜合競爭力，提高土地利用效率；棗莊市等社會滯后區(qū)，則應當加快城市產業(yè)轉型升級，加強城市基礎設施建設和提高公共服務水平；威海市同魯南地區(qū)兩側及其以北構成的半圈式生態(tài)滯后區(qū)，應當強化系統(tǒng)治理，增加對科技創(chuàng)新清潔產業(yè)的投資力度，實行“生態(tài)效益為基，社會經濟協(xié)同發(fā)展”的策略方針；濟寧市、泰安市、日照市、煙臺市等協(xié)同發(fā)展區(qū)，應當在注重保障土地集約利用政策長效性和穩(wěn)定性的同時，積極優(yōu)化區(qū)域土地利用結構，增強土地集約利用程度，推動區(qū)域可持續(xù)發(fā)展邁上新臺階；此外還應注意在區(qū)域一體化的戰(zhàn)略背景下，積極發(fā)揮省會、膠東、魯南“三圈一體化”的優(yōu)勢，加強地區(qū)間的產業(yè)協(xié)同、資源整合及技術交流，推動區(qū)域間協(xié)同發(fā)展。