廣東省土地利用效率和產業結構耦合協調特征及時空演變分析

馮維祥，廖富強，萬智巍

（1.廣東省科學院 廣州地理研究所，廣東 廣州 510070；2.江西師范大學 地理與環境學院，江西 南昌 330022；3.鄱陽湖濕地與流域研究教育部重點實驗室，江西 南昌 330022；4.贛南師范大學 地理與環境工程學院，江西 贛州 341000）

0 引言

改革開放以來，中國經濟和產業發展取得了巨大成就［1］，其中第二、三產業增加值分別從改革開放初期的0.18、0.09 萬億元增至2022 年的48.32、63.87 萬億元。作為經濟增長的關鍵要素，土地資源的合理配置和產業結構的轉型升級與現代經濟快速發展聯系緊密［2］。但在經濟高速發展的同時，也帶來諸如土地利用效率低下、產業結構不合理等諸多問題，進而影響經濟可持續發展［3］。究其原因在于長期以來粗放的發展方式［4］，當前階段，中國正處于轉變經濟發展方式的重要時期，高質量發展對土地資源利用效率和產業結構升級優化提出了更高的要求，積極探索高質量發展路徑，解決土地資源要素配置，推進產業結構合理優化，使土地利用和產業結構協調健康發展在當前經濟轉型中尤為重要。

國內外學者對于土地利用效率的測算方法，早期主要使用單一要素評價法或多要素評價法［5］計算土地利用效率，但由于其受到主觀權重的影響，逐步被其他方法所取代。近年來，當前學術界主流的方式是采用數據包絡分析（DEA）方法［6，7］和隨機前沿分析（SFA）方法［8，9］，數據包絡分析法具有非參數性的特點，該方法利用線性規劃來計算多項投入和產出指標，以此測算各項評價單元的相對效率值，但由于其前沿面固定，無法對達到前沿面的單元做進一步區分［10］，在數據包絡分析法的基礎上發展而來的含有非期望產出的超效率SBM模型［11－13］可以有效解決指標的擁擠或松弛現象以及前沿面固定的問題，并且該方法將負產出指標納入評價體系內，可以更加準確測算土地利用效率的實際表現，所以逐漸成為評價其效率的主流方法［14］。在產業結構水平測算方面，目前主要采用產業結構高級化來衡量產業結構水平，計算方法包括基于克拉克定律的非農業產值比重法［15］和以第三、第二產業增加值之比的衡量方法［16］，而空間向量夾角法因其能夠較好地反映三次產業比重沿著第一、二、三產業的順序不斷上升的規律，成為目前測算產業結構高級化水平的主要方法［17，18］。在土地利用和產業結構關系的研究中，主要集中于單向作用［19，20］、雙向作用機制［21，22］、相互影響關系［23，24］及空間的集聚溢出效應［25］等方向。綜上可見，已有文獻中對于土地利用效率與產業結構發展相互關系的成果較為豐富，但當前研究較少關注到土地利用效率和產業結構發展的耦合協調性，探究其耦合協調性的時空發展規律以及空間關聯特征在土地利用和產業結構的研究中仍有完善空間。

廣東省是全國改革開放的先行地區，具有較高土地城鎮化和產業發展水平［26］，但區域發展存在顯著的不平衡現象，而且土地利用效率和產業轉型升級仍是其面臨的重要問題，所以對土地利用效率與產業結構耦合協調水平分析，以及制定合理的土地利用和產業政策，對廣東省高質量發展具有重要意義。綜上，本文以廣東省21 個地級市作為研究對象，選取2005—2020 年共16 個年份作為研究時間段，分別采用含有非期望產出的超效率SBM模型測算土地利用效率，使用空間向量夾角法計算得出產業結構高級化指標，以此衡量產業結構水平，再運用耦合協調度模型分析土地利用效率和產業結構水平耦合協調的時空演變規律，最后基于空間自相關分析方法探索其空間相關性特征，以期為廣東省土地和產業政策制定提供科學決策依據，促進區域土地利用和產業結構協同發展。

1 研究方法、數據來源與指標體系

1.1 研究方法

1.1.1 含有非期望產出的超效率SBM模型

結合投入、期望產出和非期望產出3 個指標，本文采用由Tone 提出的含有非期望產出的超效率SBM模型測算土地利用效率，該模型將生產過程中產生的廢水、廢氣等非期望產出納入計算，不僅能有效解決投入指標過于擁擠或松弛的問題［27］，還能進一步區分處于生產前沿面的研究單元［28］。廣東省有21 個地級市即研究單元n ＝21，每個城市研究單元由4 個投入指標m、3 個期望產出指標r1和3 個非期望產出指標r2構成，投入要素的矩陣元素為x，期望產出的矩陣元素為yd，非期望產出的矩陣元素為yu，土地利用效率值為ρ，計算公式如下：

1.1.2 產業結構高級化指標

本文使用空間向量夾角法［29］測算產業結構高級化水平，該方法可以有效衡量三次產業的產值占比逐漸由低級流向高級的轉變程度。其基本原理是將國內生產總值根據第一、二、三產業劃分為3 個部分，再把各產業當年的增加值占國內生產總值的比重作為空間向量中的一個分量，隨即構成一組三維向量X0＝（x1，0，x2，0，x3，0）。再分別計算三維向量X0與第一、二、三產業空間向量X1＝（1，0，0）、X2＝（0，2，0）、X3＝（0，0，3）的夾角αj（j ＝1，2，3），計算公式如下：

再將產業結構高級化水平定義為L，L 值越大，說明產業結構的高級化水平越高，其計算公式如下：

1.1.3 耦合協調度模型

耦合協調度模型的原理是模擬兩個個體在自身和外界共同作用下產生的相互作用［30］，可以較好地表現兩者的整體性、綜合性以及協調性。本文采用該模型來測算廣東省土地利用效率和產業結構水平之間的耦合協調度，其計算公式如下：

式中：C代表耦合度；T 代表協調度；D 代表耦合協調度；U1代表土地利用效率ρ；U2代表產業結構水平，即產業結構高級化值L；a 和b 為兩個體的權重，反映個體在協調系統中的重要程度，鑒于本文認為土地利用效率和產業結構水平在經濟發展中重要程度相近，故將兩者權重均設定為0.5。

1.1.4 空間自相關分析模型

地理學中認為空間中任意事物或現象都是存在聯系的，相對鄰近的個體之間的關聯度會比非鄰近的個體間聯系程度更深［31］，所以意味著相鄰地區的土地利用效率和產業結構發展可能產生相互影響。本文通過空間自相關分析模型來研究廣東省土地利用效率和產業結構耦合協調度的空間相關性和差異性，其中全局Moran′s I指數是評價全局空間自相關常用的統計指標，計算公式如下：

式中：n 為研究單元總個數即21；xi、xj分別為研究單元i和j的耦合協調度；wij代表研究單元i和j是否相鄰，當i和j相鄰時，wij＝1，當i和j 不相鄰時，則wij＝0。全局Moran′s I 指數的取值范圍為［－1，1］，結果等于0 時說明不相關，處于［－1，0）區間說明呈空間正相關，處于（0，1］區間說明呈負相關，其絕對值越大表明空間相關性越強。

局部空間自相關可以體現某個單獨研究單元的耦合協調度同與其相鄰的研究單元的空間關聯性，通常使用Local Moran′s I指數來量化研究單元的局部空間自相關程度［32］。本文采用該指數計算廣東省土地利用效率和產業結構耦合協調度的局部空間相關性，同時反映耦合協調度在空間上熱、冷點區域的空間分布規律，計算公式如下：

式中：Zi、Zj分別為研究單元i 和j 的耦合協調度；ωij為空間權重矩陣，表現研究單元之間的連通性；Ii為研究單元i局部空間自相關指數，表示其和與它相鄰單元之間的關聯程度，∑jWijZj為相鄰研究單元偏差量的加權平均值。

1.2 數據來源

本文以廣東省21 個地級以上城市作為研究單元，研究時間跨度為中國“十五”末期到“十三五”末期，即2005—2020 年，研究中用到的產業結構、固定資產投資、二三產業從業人員、地區生產總值數據來自2006—2021 年《廣東省統計年鑒》，建成區面積數據來自2006—2021 年《中國城市建設統計年鑒》，城市公園綠地面積、地方一般公共預算收入、地方財政科技支出、工業廢水排放量、工業二氧化硫排放量、工業煙粉塵排放量來自2006—2021 年《中國城市統計年鑒》，GDP 平減指數基于統計年鑒中地區生產總值的當年價和不變價計算得到，固定資產投資價格指數、居民消費價格指數來源于國家統計局網站。部分缺失數據通過插值法進行補全，最終得到廣東省21 個地級市2005—2020 年共16 年的面板數據。本文采用的省、市行政區劃矢量數據來源于全國地理信息資源目錄服務網站（https：／ ／www.webmap.cn），數據口徑為2021 年公眾版1∶100萬基礎數據，空間坐標系統為2000 年國家大地坐標系（China Geodetic Coordinate System 2000）。

1.3 指標體系

參考王文麗、李曉靜等已有研究［11，33］，對土地利用效率投入產出的指標體系及數據統計口徑的一致性進行綜合考慮，本文選取土地、資本、勞動力、科學技術作為土地利用效率的投入指標，第二、第三產業增加值，地方一般公共預算收入，城市公園綠地面積分別代表經濟效益、社會效益、生態環境效益作為產出指標，鑒于經濟發展過程中產生廢棄物對環境造成的不利影響，對土地利用效率具有負向作用，所以本文以工業“三廢”（廢水、二氧化硫、煙粉塵排放量）作為非期望產出指標，綜上構建廣東省土地利用效率的評價指標體系（表1）。為消除時間變化帶來的價格影響，本文對部分數據進行處理：第二、第三產業增加值使用GDP 平減指數計算調整至研究初期年即2005 年，對固定資產投資額采用固定資產投資價格指數進行平減至2005 年，地方財政科技支出和地方一般公共預算收入用居民消費價格指數平減至2005 年。

表1 土地利用效率評價指標體系Table 1 Evaluation index system of land use efficiency

2 結果及分析

2.1 土地利用效率和產業結構水平

基于MaxDEA 8.0 平臺，在超效率SBM 模型、無徑向（Non－Oriented）、規模報酬可變（VRS）以及含有非期望產出（undesirable outputs）的軟件環境下，測算2005—2020 年廣東省21 個地級市土地利用效率，基于廣東省產業結構數據，采用空間向量夾角法測算得出廣東省2005—2020 年期間21 個地級市產業結構高級化指標，對研究期間內廣東省土地利用效率和產業結構高級化水平的平均值進行分析（圖1）。由圖1 可知，廣東省各地級市研究期內土地利用效率均值變化大致呈“W”型趨勢，2005 年土地利用效率最高0.99，至2007 年降至最低0.76，隨后在2007—2013 年期間回升至0.89，再次經歷短暫下降后，2015 年開始繼續回升，最終在2020 年達到0.91。廣東省各地級市產業結構高級化水平均值在研究期內總體呈現持續上升走勢，其中2005—2009 年和2013—2017 年兩個階段上升的速率較快，從2018 年開始漲勢有所停滯，并出現小幅度回落現象，從2005—2020 年廣東省產業結構高級化水平均值從6.37 增長至6.74。廣東省土地利用效率和產業結構水平在2007—2013 年和2015—2020 年兩個時間段內變化情況大體一致，說明廣東省土地利用和產業發展情況總體較為協調，且隨時間推進土地利用效率和產業結構水平均有上升趨勢。

圖1 2005—2020 年廣東省土地利用效率和產業結構水平均值Figure 1 Average level of land use efficiency and industrial structure in Guangdong Province，2005－2020

2.2 耦合協調度時間變化特征

為進一步評價廣東省土地利用效率和產業結構水平之間的協同關系，基于耦合協調度評價模型，對數據進行標準化處理，測算得到2005—2020 年廣東省21 個地級市土地利用效率和產業結構水平的耦合協調度，結合現有研究［34］建立耦合協調度的評價標準（表2）。結果表明，研究期內，廣東省各個地級市土地利用效率和產業結構之間的耦合協調度范圍處于0.156—0.862 之間，耦合協調水平跨度較大。

表2 土地利用效率和產業結構耦合協調水平評價標準Table 2 Evaluation criteria of coupling coordination of land use efficiency and industrial structure

根據廣東省土地利用效率和產業結構耦合協調度，統計2005—2020 年間各地級市的耦合協調水平類型（表3）。

表3 2005—2020 年廣東省土地利用效率和產業結構耦合協調水平類型統計Table 3 Statistics of coupling coordination between land use efficiency and industrial structure in Guangdong Province，2005－2020

由表3 可知，2005—2020 年廣東省21 個地級市耦合協調水平總體較高，到2020 年時各城市耦合協調水平均已達到勉強協調以上。具體來看，其中位于良好協調的數量較多，15 年間經歷了先減后增再減3 個階段，從15 個減少至13 個，減幅13.3%，分別為汕頭、佛山、河源、梅州、汕尾、東莞、江門、陽江、湛江、茂名、潮州、揭陽、云浮，其中江門、陽江的耦合協調水平由勉強協調達到良好協調，云浮由輕度失調達到良好協調；僅2005 年的清遠處于中度失調狀態，輕度失調的城市從2005 年的肇慶和云浮2個增至3 個，最后在2020 年降至0 個，均達到勉強協調或良好協調狀態；勉強協調的城市數量變化幅度較大，研究初期由江門和陽江2 個增加至10 個，隨后在波動中逐漸減少，最后剩下韶關、惠州、中山、肇慶、清遠5 個城市；處于優質協調狀態的城市在2005 年僅潮州1 個，經歷數次變化后2020 年達到3個，分別是廣州、深圳、珠海。總體上看，廣東省各城市土地利用效率和產業結構耦合協調度的水平隨時間推進有明顯提升，尤其體現在珠三角主要城市，可見產業結構的不斷升級優化與土地利用效率的提升具有較大關聯性，兩者的耦合協調度有顯著提升。

按照廣東省的經濟區劃將其劃分為珠三角（廣州、深圳、珠海、佛山、惠州、東莞、中山、江門、肇慶），粵東（汕頭、汕尾、潮州、揭陽），粵西（陽江、湛江、茂名），粵北（韶關、河源、梅州、清遠、云浮）四個地區，繪制各分區在研究時期內土地利用效率和產業結構之間耦合協調度均值折線圖（圖2）。

從圖2可見，2005—2020 年珠三角、粵東、粵西、粵北四大地區和廣東全省的耦合協調度均值分別為0.641、0.762、0.673、0.709、0.698，整體來看，均處在良好協調水平，說明在該時間段內土地利用投入產出效率和產業結構高級化水平協調度較好。研究期間，廣東省土地利用效率和產業結構的耦合協調度隨年份推進呈現先降后升的態勢，2005—2007 年從0.650降至0.569，2008—2014 年持續升高至0.664，經歷2015年小幅下降后，最終在2020 年達到0.698，總體增幅為7.34%。從各分區情況看，2005—2020 年廣東省土地利用效率和產業結構耦合協調度整體較高，分區差異呈現粵東＞珠三角＞粵西＞粵北，但各地區總體變化趨勢大致相當，在2005—2007 年有明顯下降，隨后從2008 年開始直至2020 年總體上均有提升。其中，珠三角地區初始耦合協調水平較高，在2006年達到0.691，經歷2007年下降隨后持久平穩升高，但上升的年均變化幅度不大，于2020 年時升高至0.709，總體增幅3.06%，但其歷年耦合協調度均值在0.6以上，總體水平良好；粵東地區變化波動較為顯著，經歷了較大幅度的升降后在2019 年達到最高值0.767，同時其歷年均值也保持在0.6 以上，研究期內增幅為5.46%；粵西地區初始協調水平較弱，達到勉強協調，2007年降至0.544，2008年起逐漸升高，期間有小幅度下降，耦合協調度在2018 年達到最高值0.685，較2000年增幅為12.48%；粵北地區的耦合協調度升降幅度較大，2000 年時為0.555，達到勉強協調水平，歷經2007—2013 年和2015—2019 年兩個時間段的顯著上升后，該地區耦合協調度在2020 年達到0.641，處于良好協調水平，研究期內整體增幅達到15.50%，為廣東省各地區最高值。總體來看，廣東省各地區土地利用效率和產業結構耦合協調水平較高，2005—2020年均有不同程度提升，珠三角和粵東地區提升幅度較小，粵西和粵北地區提升幅度較大。

2.3 耦合協調度空間變化特征

廣東省土地利用效率和產業結構的耦合協調水平存在一定地域性，為進一步探究其空間變化特征，本文選擇2005 年、2010 年、2015 年和2020 年4 個時點作為研究窗口，繪制土地利用效率和產業結構耦合協調水平狀況空間分布圖（圖3），并分析其時空變化情況。首先，從整體來看，2005年珠三角、粵東、粵西、粵北四個地區達到勉強協調以上地級市分別為8、4、2、3個，至2020年時分別為9、4、3、3 個，珠三角和粵西土地利用效率和產業結構的耦合協調水平有小幅上升趨勢。其次，土地利用效率和產業結構耦合協調空間分布的地域關聯特征顯著，珠三角、粵東、粵西地區耦合協調水平的空間格局具有較強穩定性，2005年粵東地區耦合協調度最高，達到0.723，相比珠三角、粵西、粵北分別高出0.035、0.124、0.167，到2020 年粵東地區耦合協調度仍領先其他地區，為0.762，分別高出珠三角、粵西、粵北0.053、0.089、0.121，可以發現地區間差異正在縮小，但整體具有顯著空間依賴特征。最后，廣東省各地區間耦合協調水平差異較大，但地區內部分異的程度不高，2005年珠三角、粵東、粵西、粵北地區內部各市耦合協調度的標準差分別為0.102、0.07、0.016、0.145，至2020 年時，珠三角和粵西地區耦合協調度標準差分別升高了0.035 和0.017，粵東和粵北地區標準差分別降低了0.052 和0.045，可以看出珠三角和粵西地區的內部差異在擴大，而粵東和粵北地區的內部差異在縮小，但總體變化幅度不大。

圖3 2005—2020年廣東省土地利用效率和產業結構耦合協調水平空間分布Figure 3 Spatial distribution of the coupling coordination of land use efficiency and industrial structure in Guangdong Province，2005－2020

2.4 耦合協調度空間關聯性分析

全局空間自相關可以判斷廣東省土地利用效率和產業結構是否存在空間集聚現象，通過測算得到2005—2020年土地利用效率和產業結構耦合協調度的全局Moran′s I指數（表4）。

表4 2005—2020年廣東省土地利用效率和產業結構耦合協調度全局自相關指數Table 4 Global spatial autocorrelation of coupling coordination degree between land use efficiency and industrial structure in Guangdong Province，2005－2020

由表4 可知，該指數在2005 年處于最高值0.327，至2007 年出現明顯下降，隨后在2007—2020年間始終維持在0 以下，且除2017 年以外其絕對值均未超過0.3，說明土地利用效率和產業結構耦合協調度總體呈現“集聚——隨機”分布變化趨勢，在研究初期有一定的空間依賴性，但空間聚集現象隨時間推進逐漸減弱，開始表現出隨機分布現象，2017年廣東省土地利用效率和產業結構的耦合協調度在區域分布上顯現出典型的離散性。

采用局部空間自相關分析方法，基于Arc-GIS10.2 軟件測算2005—2020 年廣東省各地級市土地利用效率和產業結構的耦合協調度的局部Moran′s I指數，在數據通過顯著性檢驗的情況下，根據計算結果的zi和∑jwijzj值將各地級市的空間相關類型劃分為4 種（表5），分別是“擴散型”“傳染型”“沉陷型”“極化型”，將未通過顯著性檢驗的地級市劃定為空間關聯特征“不顯著”，以此進一步分析各研究單元與相鄰單元的空間關聯性和依賴性特征。

表5 空間相關性類型說明Table 5 Description of spatial correlation types

根據2005—2020 年廣東省土地利用效率和產業結構耦合協調度空間關聯類型的劃分結果（表6）可以發現，土地利用效率和產業結構耦合協調性“不顯著”的地級市占比較大，研究期間平均每年數量為17.94 個，占比達到85.42%，總體來看，空間相鄰性雖對廣東省各地級市耦合協調度存在一定影響，但多數地級市受影響程度不大。16 年間擴散型、傳染型、沉陷型、極化型4 種空間相關類型總計個數分別為7、7、20、15 個，沉陷型和極化型兩種表現為負相關的地級市較多，呈現正相關類型地市較少。其中擴散型主要出現在2005—2007 年和2015—2016 年兩個時間段，均為東莞市，除此之外只有同位于珠三角地區的中山市在2011 年表現為擴散型，說明上述時間廣東省中南部有典型的空間相關性特征，耦合協調水平相比其他城市均表現為較高水平。傳染型只出現在2005 年、2006 年兩個年份，隨后直至研究末期未出現該空間相關類型的城市，2005 年陽江市、茂名、肇慶和云浮4 市表現為傳染性，主要分布在粵西及其相鄰地區，說明該年份在廣東省中西部存在一定程度的LL 集聚現象，耦合協調水平同時表現為較低，至2006 年時傳染型城市為陽江、茂名和肇慶3 市，云浮市空間關聯性轉為不顯著，粵西地區空間依賴性仍然顯著。沉陷型城市在研究時期內數量最多，平均每年1.25 個，僅出現在東莞、惠州、中山3 個城市中，均位于珠三角地區，其中惠州市在2015—2020 年持續表現為沉陷型，與之相鄰的東莞市在2010—2014 年普遍表現為沉陷型，而中山市則是在2007—2017 年大部分年份中顯示出沉陷型，說明該空間相關類型輻射范圍集中在珠三角主要城市。極化型城市數量也較多，平均每年0.94 個，分布較廣，佛山和茂名是典型的極化型城市，分別位于珠三角和粵西地區，其在研究期間與周邊城市的HL集聚現象顯著，影響范圍較大，極化效應明顯，珠海、深圳、汕尾、河源4 市也有部分年份表現為極化型，主要位于廣東省南部，上述城市與周邊城市相比耦合協調水平較高，具有極化現象。

表6 2005—2020 年廣東省土地利用效率和產業結構耦合協調度空間關聯類型統計Table 6 Spatial correlation types of coupling coordination degree of land use efficiency and industrial structure in Guangdong Province，2005－2020

將2005 年、2010 年、2015 年和2020 年4 個年份作為觀察時點，對廣東省土地利用效率和產業結構耦合協調度的空間相關類型進行可視化（圖4）。由圖4 可知，廣東省在上述年份存在3 個空間關聯性的中心區域，分別是以東莞、惠州為核心的珠三角地區，以汕尾為核心的粵東地區，以茂名、陽江為核心的粵西地區，云浮由于地理位置因素所以與珠三角和粵西地區聯系較密切。2005 年形成了以肇慶市、云浮市、陽江市和茂名市連片傳染型為主的冷點中心，與之相鄰的地級市均處于耦合協調較弱的水平，東莞市則形成了區域擴散型的熱點中心，佛山市表現為極化型，耦合協調水平普遍高于周邊地級市。隨后到2010 年，茂名市從傳染型轉變為極化型，耦合協調度有所上升，與相鄰地級市差距逐漸增大，而珠三角的東莞市、中山市則發展為沉陷型，相比與周邊城市其耦合協調的發展水平有所放緩。至2015年時，東莞市耦合協調水平逐漸提升，又成為珠三角地區的熱點中心，其周邊的中山市仍為沉陷型城市，同時惠州市發展水平也有所下滑，轉變為沉陷型，而汕尾市經歷多年發展，尤其受到珠三角地區的影響，耦合協調水平逐漸提升，成為極化型城市，與周邊的地級市差距較明顯，其他地級市則未表現出空間關聯的顯著性。最后到2020 年，廣東省土地利用效率和產業結構耦合協調水平的空間相關性逐漸減弱，大部分地級市均顯示為不顯著相關，但惠州市保持為2015 年的沉陷型，與周邊城市相比，耦合協調水平較弱，佛山市土地利用效率則在此期間發展迅速，跟上產業結構水平，帶動其耦合協調水平提升，成為極化型城市。總體來看，2005—2020 年期間廣東省土地利用效率和產業結構耦合協調水平逐漸由局部空間集聚現象有所減弱，部分城市轉變為空間性相關性不顯著，兩個典型的冷、熱點區域隨時間推進逐步消失，地理位置的鄰近性難以成為主導其耦合協調水平的關鍵因素。

圖4 廣東省土地利用效率和產業結構耦合協調度的空間相關類型Figure 4 Spatial correlation types of the coupling coordination degree between land use efficiency and industrial structure in Guangdong Province

4 結論和建議

4.1 結論

本文基于2005—2020 年廣東省21 個地級以上城市面板數據，構建土地利用效率評價指標體系，利用含有非期望產出的超效率SBM 模型以及空間向量夾角法分別計算土地利用效率和產業結構高級化程度，分析了研究區域土地利用效率和產業結構水平耦合協調度的時空分布，再結合空間自相關模型對其耦合協調度的時空演變規律和空間相關性進行探究。主要結論如下：①研究期內，廣東省土地利用效率均值變化大致呈“W”型趨勢。產業結構高級化水平均值總體持續上升，其中2005—2009 年和2013—2017 年兩個階段上升速率較快，2005—2020年產業結構水平從6.37 增長至6.74。廣東省土地利用效率和產業結構水平在研究期處于穩定上升趨勢，且兩者在大部分時間段增減趨勢類似。②從時間變化來看，土地利用效率和產業結構之間的耦合協調度范圍處于0.156—0.862 之間，耦合協調水平跨度較大，但大部分城市處于良好協調狀態，且隨時間推進有明顯提升，至2020 年時各城市耦合協調水平均已達到勉強協調以上，水平總體較高。③從空間變化來看，耦合協調度地域關聯特征顯著，珠三角、粵東、粵西地區耦合協調水平的空間格局具有較強的穩定性，差異呈現粵東＞珠三角＞粵西＞粵北，珠三角和粵東地區提升幅度較小，粵西和粵北地區提升幅度較大，地區間差異正在縮小，但珠三角和粵西地區的內部差異在擴大，粵東和粵北地區的內部差異在縮小。④從空間相關性分析來看，廣東省2005—2020 年耦合協調度Moran′s I指數逐漸降低，總體呈現“集聚——隨機”分布變化趨勢，空間聚集現象隨時間推進逐漸減弱，16 年間擴散型、傳染型、沉陷型、極化型4 種空間相關類型總計個數分別為7、7、20、15 個，呈現空間負相關的地級市較多，而呈現正相關的地級市較少。隨時間推進廣東省兩個典型的冷、熱點區域逐步消失，說明土地利用效率和產業結構水平的耦合協調度受地理位置鄰近性的影響程度越來越弱。

4.2 建議

針對本文研究結論，為優化廣東省用地結構和生產要素配置，提升土地利用效率城市和產業結構的耦合協調度，促進區域集約節約用地和高質量協調發展，提出以下針對性建議：

對土地利用效率和產業結構耦合協調度水平較低的地區，根據其產生原因，使用差異化的土地與產業調控手段。在現有研究基礎上，利用土地利用效率和產業結構高級化水平評價結果對耦合協調水平的主導因素進行分析，主要有以下兩種情況：第一種是由于土地利用效率低水平所導致的，主要以云浮市、汕尾市、茂名市為主。這類地級市應以改善土地利用模式，提高土地利用效率為重點，首先應在編制地區國土空間規劃時同時考慮其資源稟賦和承載能力，確定合理的新增建設用地規模，著力實施存量土地盤活再開發，達到建設用地集約節約化利用，其次應加大對區域內生產企業的整治力度，提高企業生產技術和內部管理方式，引導資金投入，鼓勵企業向高新產業轉型，降低無效資源投入，同時還應強調企業的環境保護意識，減少非期望產出，在提升土地利用效率的同時實現綠色發展。第二種是由于產業結構低水平所導致的，主要以東莞市、中山市、清遠市為主。這類地級市應重點推動產業結構升級，基于各地區的資源要素、產業基礎和發展需求，因地制宜提升產業水平，優化產業結構。具體而言，珠三角地區應在其經濟水平優勢的基礎上，重點發展新興化、高級化的第三產業，逐步退出高污染、高能耗的低端產業，為城市高質量發展騰挪土地空間，而粵東、粵西、粵北地區則可以在勞動力及資源要素優勢的基礎上，承接珠三角地區的產業轉移，重點發展高端制造業，淘汰低附加值產業，實現地區間協調發展，部分產業結構高級化水平相對落后的地區，可通過建設高新技術開發區和科技產業園等手段，針對性地吸引高附加值企業和高端產業入駐，以達到產業結構的有效提升和可持續發展。

目前廣東省土地利用效率和產業結構耦合協調度的空間依賴現象有所減弱，地區間差異在逐漸縮小，但地區之間依然存在一定差異，部分地級市仍瀕臨失調，所以廣東省應重視一體化發展，同步推進區域統籌協調，根據各地市情況靈活施策。對于以珠三角地區的東莞市和中山市為主的擴散型空間集聚區，應引導其繼續良性發展，基于其合理的土地利用方式和產業結構基礎，進一步優化土地要素配置，推進產業升級和發展，實現土地利用效率和產業結構水平兩者間的助力和交互提升，并在此基礎上，擴大影響范圍，帶動相鄰城市協調發展。對于以粵西地區的陽江市和茂名市為主的傳染型空間集聚區，應通過適當的政策引導，鼓勵其積極學習先進地區的土地利用和產業發展政策，整治低效閑置用地，提升土地儲備，為承接發達地區的產業和創新要素轉移打好現實基礎，提前做好產業服務。對于以珠三角地區的惠州市、中山市和佛山市為主的極化型和沉陷型空間集聚區，在保證自身良性發展的前提下，應提倡高水平帶動低水平地區，相互協助，實現落后地區土地利用效率與產業結構的協調發展的提升，進而發展成為區域熱點中心，并且繼續擴大影響范圍。

總體而言，基于上述研究結論，應統籌考慮廣東省各地級市土地利用效率與產業結構升級的相互影響關系，促進兩者協同發展。對土地利用效率和產業結構耦合協調水平較高的地市，應鼓勵其主動建立地區間的聯系和協同發展渠道，通過產業轉移和創新示范等方式，帶動周邊耦合協調度水平較低的地市發展。針對水平相對落后的地市，應重視土地利用和產業協調發展，通過完善相關政策，采取與其他高水平地區聯動和建立共享機制等方式，達到區域產業互補、合理分工，實現全省高質量協調發展。