江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度的時(shí)空格局演變研究

呂 曉，史洋洋

（1.曲阜師范大學(xué)地理與旅游學(xué)院，山東 日照 276826；2.國(guó)土資源部海岸帶開發(fā)與保護(hù)重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210024）

隨著城鎮(zhèn)化、工業(yè)化進(jìn)程的快速推進(jìn)，與之相伴隨的土地利用粗放低效、城鄉(xiāng)發(fā)展不協(xié)調(diào)等問(wèn)題漸次出現(xiàn)，并進(jìn)一步導(dǎo)致區(qū)域發(fā)展差異拉大、資源環(huán)境約束加劇[1-2]。如何突破有限的土地供給和發(fā)展需求之間的瓶頸，合理規(guī)劃布局，提升土地利用效率和節(jié)約集約用地水平，已成為社會(huì)各界關(guān)注的焦點(diǎn)[3-5]。城市土地經(jīng)濟(jì)密度作為衡量土地利用效率的重要指標(biāo)已被廣泛應(yīng)用于城市經(jīng)濟(jì)發(fā)展、土地集約利用等相關(guān)研究[6-7]。城鄉(xiāng)建設(shè)用地是區(qū)域經(jīng)濟(jì)社會(huì)發(fā)展的主要空間載體，提高其節(jié)約集約利用水平、綜合效益以及開發(fā)利用有序程度是解決土地供需矛盾的重要途徑。相較于城市建成區(qū)或城鎮(zhèn)用地經(jīng)濟(jì)密度[8-10]，城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度涵蓋了鄉(xiāng)村建設(shè)用地對(duì)區(qū)域經(jīng)濟(jì)的貢獻(xiàn)，能更全面地反映土地利用效益。同時(shí)，城鄉(xiāng)建設(shè)用地在空間上具有明顯的連續(xù)性和自相關(guān)性，如空間相互作用和空間擴(kuò)散效應(yīng)，從經(jīng)濟(jì)和城鄉(xiāng)建設(shè)用地相結(jié)合的角度分析城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度的時(shí)空格局及其空間關(guān)聯(lián)[11-12]對(duì)進(jìn)一步揭示城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)效益的時(shí)空差異、推動(dòng)區(qū)域經(jīng)濟(jì)協(xié)調(diào)發(fā)展具有重要價(jià)值。

江蘇省作為中國(guó)經(jīng)濟(jì)發(fā)達(dá)區(qū)和人口密集地區(qū)，城鎮(zhèn)化、工業(yè)化發(fā)展迅速，2015年末城鎮(zhèn)化率達(dá)到66.52%，高于全國(guó)平均水平約10個(gè)百分點(diǎn)，是經(jīng)濟(jì)發(fā)達(dá)與人地矛盾并存的典型地區(qū)?；诖?，本文以江蘇省為例運(yùn)用全局和局部空間自相關(guān)分析方法，結(jié)合變異系數(shù)、泰爾指數(shù)、相對(duì)發(fā)展率等指標(biāo)，選擇2005、2010和2015年三個(gè)時(shí)點(diǎn)對(duì)江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空格局及動(dòng)態(tài)演變特征進(jìn)行分析，以期為提升城鄉(xiāng)建設(shè)用地集約利用水平、優(yōu)化用地與產(chǎn)業(yè)結(jié)構(gòu)、改進(jìn)土地利用政策等提供依據(jù)。

1 數(shù)據(jù)來(lái)源與研究方法

1.1 數(shù)據(jù)來(lái)源

本文城鄉(xiāng)建設(shè)用地?cái)?shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)——長(zhǎng)江三角洲數(shù)據(jù)共享平臺(tái)（http://nnu.geodata.cn）。該數(shù)據(jù)底圖來(lái)源于2005年以來(lái)的Landsat TM數(shù)字遙感影像、航測(cè)數(shù)據(jù)和其他輔助歷史資料等，通過(guò)對(duì)遙感圖像進(jìn)行融合增強(qiáng)、幾何矯正和圖形鑲嵌處理，并對(duì)所獲數(shù)據(jù)進(jìn)行目視解譯、空間位置配準(zhǔn)和精度檢驗(yàn)，土地利用一級(jí)類型綜合評(píng)價(jià)精度達(dá)到94.30%以上[13-15]，滿足省域?qū)用娉青l(xiāng)建設(shè)用地經(jīng)濟(jì)密度差異研究精度。江蘇省13個(gè)地級(jí)城市二、三產(chǎn)業(yè)增加值數(shù)據(jù)來(lái)源于相應(yīng)年份的《江蘇省統(tǒng)計(jì)年鑒》。

1.2 研究方法

本文所指城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度是單位面積城鄉(xiāng)建設(shè)用地上的二、三產(chǎn)業(yè)增加值。依據(jù)《土地利用現(xiàn)狀分類（GB/T 21010-2007）》和《市（地）級(jí)土地利用總體規(guī)劃編制規(guī)程（TD/T 1023-2010）》，同時(shí)考慮遙感影像解譯數(shù)據(jù)包含了與城鎮(zhèn)村緊密相連的部分交通運(yùn)輸和機(jī)場(chǎng)用地，因此本文選用的城鄉(xiāng)建設(shè)用地?cái)?shù)據(jù)為城鎮(zhèn)村及工礦用地，以及部分交通運(yùn)輸和機(jī)場(chǎng)用地的集合。

首先，通過(guò)變異系數(shù)（Cv）、泰爾指數(shù)（T）、相對(duì)發(fā)展率（Nich）等指標(biāo)來(lái)考察江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度的時(shí)空分異。其中，變異系數(shù)可以反映數(shù)據(jù)空間離散程度，泰爾指數(shù)具有可分解性，能夠通過(guò)將總體差異分解為組內(nèi)與組間差異，進(jìn)而得到組內(nèi)差異與組間差異對(duì)總體差異的貢獻(xiàn)率，具體公式詳見已有文獻(xiàn)[8，10]。

相對(duì)發(fā)展率可進(jìn)一步反映各城市在一定時(shí)間階段內(nèi)相對(duì)于江蘇省總體的發(fā)展速度，計(jì)算公式為[8]：

式（1）中，Nich是相對(duì)發(fā)展率；Yi2、Yi1分別是第i個(gè)城市末期、初期的城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度值；Y2、Y1分別是江蘇省末期、初期的城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度值。

探索性空間數(shù)據(jù)分析（ESDA，Exploratory Spatial Data Analysis）是基于空間關(guān)聯(lián)度，對(duì)事物的時(shí)空格局進(jìn)行描述和可視化分析，分為全局和局部空間自相關(guān)兩部分，現(xiàn)在已經(jīng)廣泛應(yīng)用于區(qū)域經(jīng)濟(jì)與地理空間差異研究[16-17]。全局空間自相關(guān)是描述研究對(duì)象在區(qū)域整體上的空間集聚程度，分析研究對(duì)象間相對(duì)位置的空間依賴性及其顯著性[18]。本文選用Moran’s I指數(shù)分析全局空間自相關(guān)，公式為：

式（2）中，I是全局Moran值；n是研究單元的數(shù)目；xi和xj分別為區(qū)域i和j的城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度值；為研究區(qū)域13城市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度值平均數(shù)；wij為空間單元i和j的相鄰權(quán)重矩陣，本文基于Rook’s鄰接標(biāo)準(zhǔn)構(gòu)建江蘇省13城市相鄰權(quán)重矩陣，相鄰接時(shí)值為1，否則為值0；Moran’s I的取值范圍為（-1，1），0＜I＜1表示存在正空間自相關(guān)，數(shù)值越大表示空間分布相關(guān)性越大，空間上越集聚，-1＜I＜0表示負(fù)空間自相關(guān)，數(shù)值越小代表相關(guān)性越小，I= 0表示不存在空間自相關(guān)性。其假設(shè)檢驗(yàn)采用得分檢驗(yàn)方式，即置信水平P≤0.05，|Z|＞1.96經(jīng)驗(yàn)閥值，可認(rèn)為其存在顯著的空間自相關(guān)。

局部空間自相關(guān)能夠描述區(qū)域間的空間關(guān)聯(lián)效應(yīng)，反映區(qū)域與其周圍相鄰地區(qū)的相互關(guān)系[19]。結(jié)合Moran散點(diǎn)圖和LISA集聚圖進(jìn)行可視化表達(dá)，Moran散點(diǎn)圖能反映區(qū)域局部不穩(wěn)定性和空間異常性，LISA集聚圖能反映每個(gè)研究單元的空間自相關(guān)性[20-22]，公式為：

式（3）中：Ii值是Moran’s I的分解形式，其余變量定義及檢驗(yàn)同式（2）。Ii值為正，表示i區(qū)域單元周圍存在高值或低值（HH或LL類型）空間集聚；Ii值為負(fù)，表示i區(qū)域單元周圍存在高低或低高值（HL或LH類型）空間集聚。ESDA相關(guān)計(jì)算通過(guò)Geoda軟件實(shí)現(xiàn)。

2 城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度的時(shí)空分異

2.1 時(shí)空分異的總體特征

2015年江蘇省城鄉(xiāng)建設(shè)用地平均經(jīng)濟(jì)密度為2.71×109元/km2，最高的是蘇州市，最低的是連云港市，分別為5.72×109元/km2和0.77×109元/km2，差距明顯（表1）。其中，7個(gè)城市的城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度超過(guò)平均水平，這些城市的城鄉(xiāng)建設(shè)用地面積占全省的44.97%，其二三產(chǎn)業(yè)增加值占全省的73.10%。低于江蘇省平均水平的城市有6個(gè)，這些城市的城鄉(xiāng)建設(shè)用地面積占到全省的55.03%，但其二三產(chǎn)業(yè)增加值僅占全省的26.90%。

（1）城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度總體呈增長(zhǎng)態(tài)勢(shì)。在2005—2015年的10年間，江蘇省城鄉(xiāng)建設(shè)用地平均經(jīng)濟(jì)密度由0.86×109元/km2上升到2.71×109元/km2，增長(zhǎng)了近3倍。盡管無(wú)錫市和鎮(zhèn)江市在2010—2015年間城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度有所下降，但其余城市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度都呈上漲趨勢(shì)，總體增長(zhǎng)顯著。

（2）城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度呈現(xiàn)空間集聚特征。由圖1可見，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度呈現(xiàn)明顯的空間集聚特征。2005年，無(wú)錫市最高、南京市和南通市次之，相較而言，南京市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度雖然稍低于南通市，但其輻射帶動(dòng)作用要比南通市強(qiáng)，與南京市相鄰的鎮(zhèn)江市和常州市經(jīng)濟(jì)密度要高于相鄰于南通市的蘇州市和泰州市。2010年，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度整體上升，無(wú)錫市仍明顯高于其他地區(qū)，蘇州市位居第二，蘇錫常地區(qū)成為城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度最高的區(qū)域。2015年，蘇錫常仍舊是城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度最高的區(qū)域，南京市仍是經(jīng)濟(jì)密度高值區(qū)域，南通市在沿海開發(fā)國(guó)家戰(zhàn)略實(shí)施進(jìn)程中發(fā)展迅速。而宿遷市和連云港市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度遠(yuǎn)低于省內(nèi)其他城市。

區(qū)域經(jīng)濟(jì)發(fā)展水平（GDP）與土地利用集約度呈正的相關(guān)性[23]。隨著經(jīng)濟(jì)發(fā)展水平的不斷提高，建設(shè)用地集約利用水平增加量雖呈先增加后減少的趨勢(shì)，但集約利用水平整體呈現(xiàn)上升趨勢(shì)并最終趨于穩(wěn)定[24]。從江蘇省3大區(qū)域來(lái)看，位于蘇南的蘇錫常地區(qū)緊鄰上海市，是長(zhǎng)江三角洲重要組成區(qū)域，區(qū)位優(yōu)勢(shì)明顯，外向型經(jīng)濟(jì)發(fā)達(dá)，城鄉(xiāng)建設(shè)用地集約度高，經(jīng)濟(jì)密度大；南京市作為省會(huì)城市和歷史文化名城，經(jīng)濟(jì)發(fā)展水平較高，土地集約利用水平與城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度也較大。蘇中和蘇北由于缺乏區(qū)位優(yōu)勢(shì)，加之交通、政策、觀念等因素，經(jīng)濟(jì)發(fā)展水平相對(duì)落后，城鄉(xiāng)建設(shè)用地利用相對(duì)粗放，土地集約利用水平低。

（3）城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空差異顯著，區(qū)域核心—邊緣效應(yīng)明顯。通過(guò)分析2005、2010和2015年三個(gè)時(shí)點(diǎn)的城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度空間分布，發(fā)現(xiàn)江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度存在顯著的時(shí)空差異，總體上蘇南大于蘇中，蘇中大于蘇北，呈現(xiàn)出由南向北的分層遞減態(tài)勢(shì)。2015年高于江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度平均水平的城市有7個(gè)，其中蘇南4個(gè)，蘇中3個(gè)，蘇北各城市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度均低于平均水平。這與上海市快速發(fā)展產(chǎn)生的強(qiáng)大輻射帶動(dòng)作用密切相關(guān)。同時(shí)，從全省層面看，南京市、無(wú)錫市等城市對(duì)周圍城市輻射帶動(dòng)效果也深刻體現(xiàn)了核心—邊緣效應(yīng)。

表1 江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度分布表Tab.1 The distribution of economic density of URCL in Jiangsu Province

2.2 區(qū)域差異動(dòng)態(tài)演變特征

（1）總體區(qū)域差異減小。根據(jù)變異系數(shù)和泰爾指數(shù)分析發(fā)現(xiàn)，2005—2015年江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度區(qū)域差異總體上不斷減小（表2）。其中，空間離散程度在2005—2010年階段下降幅度較大，在2010—2015年階段下降趨勢(shì)放緩。城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度空間總體差異呈現(xiàn)持續(xù)減小特征，即2005—2010年階段空間差異減小，泰爾指數(shù)由0.53下降到0.48，而2010—2015年階段泰爾指數(shù)由0.48下降到0.40，相較于上一階段空間差異減小幅度更大。

圖1 江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空格局演變Fig.1 The spatio-temporal pattern evolution of economic density of URCL in Jiangsu Province

（2）區(qū)間和區(qū)內(nèi)差異持續(xù)減小。根據(jù)江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空分異總體特征及泰爾指數(shù)的空分解特性，區(qū)域差異由蘇北、蘇南和蘇中三個(gè)地區(qū)的區(qū)間差異及各地區(qū)的區(qū)內(nèi)差異構(gòu)成[8，11]。由表2可見，在2005—2015年間，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度總體差異值不斷減小，說(shuō)明江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)效益的地區(qū)差異持續(xù)減小，總體上向均衡的方向發(fā)展。區(qū)間差異與總體差異發(fā)展態(tài)勢(shì)趨同，區(qū)間差異值由2005年的0.3543降低到2010年的0.3204，2015年又急劇下降到0.2702，區(qū)間差異縮小態(tài)勢(shì)顯著。江蘇省區(qū)間差異值與所占比重呈反向發(fā)展趨勢(shì)，2005年區(qū)間差異值最高，占總體差異的66.25%，說(shuō)明2005年之前階段，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度發(fā)展很不平衡，各地區(qū)城市間經(jīng)濟(jì)出現(xiàn)不均衡發(fā)展。2005—2015年間，江蘇省區(qū)間差異所占比重不斷上升，而區(qū)間差異值不斷減小，表明地區(qū)間城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度相對(duì)協(xié)調(diào)發(fā)展，但三大地區(qū)間發(fā)展仍不平衡。

由表2可知，江蘇省三大地區(qū)的區(qū)內(nèi)差異仍呈持續(xù)減小態(tài)勢(shì)。2005年江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度的區(qū)內(nèi)差異值為0.1804，區(qū)內(nèi)差異最為顯著的是蘇北，對(duì)全省總體差異貢獻(xiàn)率達(dá)到30.23%，其次是蘇南和蘇中。2010年全省區(qū)內(nèi)差異增長(zhǎng)到0.1591，而區(qū)內(nèi)差異對(duì)總體差異的貢獻(xiàn)率由33.74%下降到33.18%，其中蘇北的區(qū)內(nèi)差異占總體差異比重仍然最大，說(shuō)明2005—2010年階段，蘇北內(nèi)部各市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度發(fā)展極不均衡。2010—2015年階段，整體上區(qū)內(nèi)差異呈急劇下降態(tài)勢(shì)，由0.1591下降到0.1249，下降幅度較大，其中蘇北的區(qū)內(nèi)差異占總體差異的比重遠(yuǎn)高于蘇中和蘇南地區(qū)，蘇中變化最小，蘇南有所上升，但其只占總體差異的6.02%，相較于蘇北的24.98%，蘇南的區(qū)內(nèi)差異對(duì)總體差異貢獻(xiàn)很小。

表2 江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度動(dòng)態(tài)差異Tab.2 The dynamic disparity of economic density of URCL in Jiangsu Province

3 城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空演化特征

3.1 區(qū)域經(jīng)濟(jì)發(fā)展的時(shí)空分異

（1）城市間相對(duì)發(fā)展速度存在明顯區(qū)域差異。蘇州市、常州市的發(fā)展速度明顯快于其他城市，南京市、揚(yáng)州市、泰州市和南通市發(fā)展速度緊隨其后，發(fā)展速度最慢的是連云港市，其相對(duì)發(fā)展率僅為2.56%，發(fā)展最快的蘇州市相對(duì)發(fā)展率是它的近8倍。

（2）城市發(fā)展速度存在明顯的時(shí)間階段性特征。2005—2010年階段，城市發(fā)展速度快的蘇州市、無(wú)錫市在2010—2015年階段發(fā)展速度下降，其中無(wú)錫市作為2005—2010年階段發(fā)展速度較快的城市成為2010—2015年階段負(fù)增長(zhǎng)城市，發(fā)展速度較快的鎮(zhèn)江市在2010—2015年階段也出現(xiàn)了負(fù)增長(zhǎng)現(xiàn)象。其余城市發(fā)展速度都有所提高，其中常州市、南京市、揚(yáng)州市、南通市和泰州市發(fā)展速度增長(zhǎng)最為顯著。

（3）相對(duì)發(fā)展速度呈現(xiàn)明顯的空間分異特征。蘇州市、常州市和南京市仍保持較高的發(fā)展速度，而原經(jīng)濟(jì)密度較高的無(wú)錫市、鎮(zhèn)江市發(fā)展速度相對(duì)降低，原經(jīng)濟(jì)密度較低的揚(yáng)州市、泰州市相對(duì)發(fā)展速度提升很快，低經(jīng)濟(jì)密度的蘇北城市發(fā)展速度穩(wěn)步上升，從而形成了新的空間格局。城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度與區(qū)域經(jīng)濟(jì)發(fā)展和城鎮(zhèn)化水平密切相關(guān)，不同的經(jīng)濟(jì)發(fā)展階段，城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度與經(jīng)濟(jì)發(fā)展之間的關(guān)系不同。2000年以來(lái)，隨著城鎮(zhèn)化、工業(yè)化進(jìn)程的整體加快，城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度不斷提高，但各地“攤大餅”式的盲目擴(kuò)張帶來(lái)了土地粗放利用，城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度增長(zhǎng)有所減緩。依靠大量消耗土地來(lái)?yè)Q取經(jīng)濟(jì)增長(zhǎng)的粗放式擴(kuò)張模式已經(jīng)無(wú)法適應(yīng)綠色發(fā)展要求，集約高效利用城鄉(xiāng)建設(shè)用地，以轉(zhuǎn)變土地利用方式倒逼經(jīng)濟(jì)轉(zhuǎn)型成為城鄉(xiāng)轉(zhuǎn)型發(fā)展的有效途徑。

3.2 空間格局演化的全局特征

運(yùn)用Geoda軟件對(duì)2005、2010和2015年江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度值進(jìn)行空間自相關(guān)分析，計(jì)算得出全局自相關(guān)系數(shù)Moran’s I值均為正，分別為0.4097、0.6169和0.6056（圖2）。整體上，2005—2015年間，江蘇省各市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度呈現(xiàn)出顯著的正相關(guān)關(guān)系，空間上存在明顯區(qū)域集聚特征。

2005—2010年，Moran’s I值出現(xiàn)較大幅度上升，說(shuō)明各城市間的差異在不斷縮小，可見蘇南經(jīng)濟(jì)向北方輻射擴(kuò)散以及沿海開發(fā)戰(zhàn)略的實(shí)施使得沿海地區(qū)向內(nèi)陸的經(jīng)濟(jì)輻射和帶動(dòng)作用有效提升了整體經(jīng)濟(jì)水平。這也驗(yàn)證了江蘇省不斷推進(jìn)城鄉(xiāng)一體化發(fā)展、縮小蘇北、蘇中與蘇南地區(qū)差距以及促進(jìn)區(qū)域全面協(xié)調(diào)可持續(xù)發(fā)展的政策實(shí)施效用。相較于2010年，2015年Moran’s I值出現(xiàn)了小幅下降，但仍比2005年要大，表明這一階段江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度空間分布相關(guān)性減弱，但總體上空間差異呈不斷減小趨勢(shì)。

3.3 空間格局演化的局部特征

全局空間自相關(guān)并不能反映出局部空間自相關(guān)的顯著差異，為全面分析江蘇省地區(qū)間的空間關(guān)聯(lián)程度和時(shí)空演化特征，基于式（3），將不同年份的Ii值導(dǎo)入ArcGIS10.2軟件，將江蘇省各年份的Lisa集聚圖（圖3）和Moran’s散點(diǎn)圖（圖2）對(duì)照分析，探討各城市間的局部空間自相關(guān)性。

結(jié)合圖2和圖3分析可得，2005年蘇北出現(xiàn)LL類型集聚，而HH類型集聚出現(xiàn)在了常州市，表明蘇北整體經(jīng)濟(jì)密度較低，蘇中城市集聚不明顯，蘇南發(fā)展相對(duì)較好。南京市沒(méi)有出現(xiàn)明顯集聚，說(shuō)明相對(duì)而言南京市經(jīng)濟(jì)效益與蘇南仍有差距。在2005年蘇北低值集聚基礎(chǔ)上，2010年以無(wú)錫市為中心的HH類型集聚向蘇州市轉(zhuǎn)移，說(shuō)明經(jīng)濟(jì)密度高值集聚城市群進(jìn)一步擴(kuò)散。2015年出現(xiàn)了新的HH集聚中心無(wú)錫市和泰州市，說(shuō)明由于蘇錫常地區(qū)經(jīng)濟(jì)崛起產(chǎn)生的經(jīng)濟(jì)輻射擴(kuò)散效果顯著，蘇南經(jīng)濟(jì)發(fā)展成效開始向北方擴(kuò)散，高值集聚中心開始由單個(gè)城市向城市連綿區(qū)方向發(fā)展，但蘇北仍是LL類型集聚地區(qū)。

總體來(lái)看，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度空間集聚現(xiàn)象明顯。高值集聚（HH）主要在蘇南和蘇中地區(qū)，低值集聚（LL）主要在蘇北地區(qū)，南北差異顯著，這與江蘇省經(jīng)濟(jì)發(fā)展水平的時(shí)空特征有密切關(guān)系。蘇錫常地區(qū)作為江蘇省經(jīng)濟(jì)發(fā)展的增長(zhǎng)極，對(duì)區(qū)域經(jīng)濟(jì)擴(kuò)散和輻射帶動(dòng)作用十分明顯。

圖2 江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度Moran散點(diǎn)圖Fig.2 Moran scatter plots of economic density of URCL in Jiangsu Province

圖3 江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度LISA集聚圖Fig.3 LISA cluster map of economic density of URCL in Jiangsu Province

4 結(jié)論與討論

本文基于ArcGIS 10.2和Geoda軟件，綜合運(yùn)用全局和局部空間自相關(guān)分析方法，并結(jié)合變異系數(shù)、泰爾指數(shù)、相對(duì)發(fā)展率，對(duì)江蘇省2005、2010和2015年13個(gè)地級(jí)城市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空格局及動(dòng)態(tài)演變特征進(jìn)行分析，得出以下結(jié)論：

（1）2005—2015年10年間，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度總體呈增長(zhǎng)態(tài)勢(shì)，空間集聚特征顯著，高值點(diǎn)主要集中于蘇錫常地區(qū)、南京市和南通市。城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度總體上南方大于北方，呈現(xiàn)出由南向北的分層遞減態(tài)勢(shì)，區(qū)域城市之間核心—邊緣效應(yīng)明顯。

（2）2005—2015年10年間，江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度空間離散程度不斷下降，總體差異和區(qū)間差異呈持續(xù)減小態(tài)勢(shì)，說(shuō)明江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)效益的地區(qū)差異持續(xù)減小，總體上向均衡的方向發(fā)展；蘇南和蘇中地區(qū)組合的區(qū)內(nèi)差異雖有上升，但不及蘇北地區(qū)組合的區(qū)內(nèi)差異下降幅度大，整體上仍呈現(xiàn)不斷減小態(tài)勢(shì)。

（3）江蘇省城市發(fā)展速度存在明顯的時(shí)間階段性特征，城市間區(qū)域差異與空間分異特征顯著。各城市城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度在空間演變上存在明顯的區(qū)域集聚特征：高值集聚（HH）主要在蘇南和蘇中地區(qū)，其中蘇錫常地區(qū)經(jīng)濟(jì)輻射擴(kuò)散效果顯著，蘇南經(jīng)濟(jì)發(fā)展成效開始向北方遷移，城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度高值集聚中心由單個(gè)城市向城市連綿區(qū)方向發(fā)展。蘇北地區(qū)一直是城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度低值集聚（LL）地區(qū)，南北差異顯著。

（4）江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度南北差異表明，蘇北地區(qū)通過(guò)盲目擴(kuò)張和粗放利用土地來(lái)?yè)Q取經(jīng)濟(jì)發(fā)展的模式已經(jīng)無(wú)法適應(yīng)當(dāng)前的經(jīng)濟(jì)社會(huì)發(fā)展要求，以蘇南、蘇中地區(qū)為代表的追求高效集約的城鄉(xiāng)土地利用才是順應(yīng)經(jīng)濟(jì)發(fā)展階段的有效手段。目前，已有研究多側(cè)重研究城市建成區(qū)的經(jīng)濟(jì)密度，而本文將鄉(xiāng)村建設(shè)用地對(duì)于區(qū)域經(jīng)濟(jì)的貢獻(xiàn)也考慮其中，更能反映基于行政區(qū)劃單元的區(qū)域整體發(fā)展水平。此外，本文研究數(shù)據(jù)來(lái)源于遙感影像解譯的土地利用數(shù)據(jù)，在宏觀尺度上更能真實(shí)反映區(qū)域土地利用實(shí)際。當(dāng)然，不同空間尺度對(duì)于區(qū)域差異特征的反映效果不同，本文基于市域?qū)用鎸?duì)江蘇省城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度時(shí)空格局演變特征進(jìn)行了分析，下一步有待對(duì)縣域、鎮(zhèn)域等多尺度的區(qū)域差異特征、驅(qū)動(dòng)機(jī)制及其空間演化等進(jìn)行分析，以期能夠更細(xì)致地探討區(qū)域城鄉(xiāng)建設(shè)用地經(jīng)濟(jì)密度發(fā)展特征及形成機(jī)制，為高效集約用地及城鄉(xiāng)協(xié)調(diào)可持續(xù)發(fā)展提供科學(xué)依據(jù)。

參考文獻(xiàn)（References）：

［1］關(guān)興良，魏后凱，魯莎莎，等. 中國(guó)城鎮(zhèn)化進(jìn)程中的空間集聚、機(jī)理及其科學(xué)問(wèn)題［J］.地理研究，2016，35（2）：227 - 241.

［2］呂曉，黃賢金，鐘太洋，等. 土地利用規(guī)劃對(duì)建設(shè)用地?cái)U(kuò)張的管控效果分析——基于一致性與有效性的復(fù)合視角［J］.自然資源學(xué)報(bào)，2015，30（2）：177 - 187.

［3］王德起，侯圣銀. 基于引力模型的京津冀城市群土地利用強(qiáng)度研究［J］.土地經(jīng)濟(jì)研究，2016，（2）：107 - 125.

［4］夏方舟，嚴(yán)金明，徐一丹. 土地要素對(duì)中國(guó)經(jīng)濟(jì)增長(zhǎng)技術(shù)效率的隨機(jī)邊界分析［J］.土地經(jīng)濟(jì)研究，2014，（1）：91 -103.

［5］呂曉，牛善棟，張全景，等. 基于內(nèi)容分析法的集體建設(shè)用地流轉(zhuǎn)政策演進(jìn)分析［J］.中國(guó)土地科學(xué)，2015，29（4）：25 - 33.

［6］匡兵，盧新海，周敏. 中國(guó)城市土地經(jīng)濟(jì)密度的分布動(dòng)態(tài)演進(jìn)［J］.中國(guó)土地科學(xué)，2016，30（10）：47 - 54.

［7］吳一洲，吳次芳，羅文斌. 經(jīng)濟(jì)地理學(xué)視角的城市土地經(jīng)濟(jì)密度影響因素及其效應(yīng)［J］.中國(guó)土地科學(xué)，2013，27（1）：26 - 33.

［8］貝涵璐，吳次芳，馮科，等. 土地經(jīng)濟(jì)密度的區(qū)域差異特征及動(dòng)態(tài)演變格局——基于長(zhǎng)江三角洲地區(qū)的實(shí)證分析［J］.自然資源學(xué)報(bào)，2009，24（11）：1952 - 1962.

［9］呂曉，黃賢金，張全景. 城鄉(xiāng)建設(shè)用地轉(zhuǎn)型研究綜述［J］ .城市規(guī)劃，2015，39（4）：105 - 112.

［10］張改素，丁志偉，趙萌，等. 中原經(jīng)濟(jì)區(qū)縣域經(jīng)濟(jì)密度的空間分異及影響因素［J］.經(jīng)濟(jì)地理，2014，34（9）：19 -26，39.

［11］王海燕，濮勵(lì)杰，張健，等. 城鎮(zhèn)建設(shè)用地經(jīng)濟(jì)密度時(shí)空分異的實(shí)證分析［J］.中國(guó)土地科學(xué)，2012，26（4）：47 - 53.

［12］曹廣忠，白曉. 中國(guó)城鎮(zhèn)建設(shè)用地經(jīng)濟(jì)密度的區(qū)位差異及影響因素——基于273個(gè)地級(jí)及以上城市的分析［J］ .中國(guó)人口·資源與環(huán)境，2010，20（2）：12 - 18.

［13］呂曉，史洋洋，黃賢金，等. 江蘇省土地利用變化的圖譜特征［J］.應(yīng)用生態(tài)學(xué)報(bào)，2016，27（4）：1077 - 1084.

［14］劉紀(jì)遠(yuǎn)，劉文超，匡文慧，等. 基于主體功能區(qū)規(guī)劃的中國(guó)城鄉(xiāng)建設(shè)用地?cái)U(kuò)張時(shí)空特征遙感分析［J］.地理學(xué)報(bào)，2016，71（3）：355 - 369.

［15］陳昌玲，張全景，呂曉，等. 江蘇省耕地占補(bǔ)過(guò)程的時(shí)空特征及驅(qū)動(dòng)機(jī)理［J］.經(jīng)濟(jì)地理，2016，36（4）：155 - 163.

［16］馬曉東，馬榮華，徐建剛. 基于ESDA-GIS的城鎮(zhèn)群體空間結(jié)構(gòu)［J］.地理學(xué)報(bào)，2004，59（6）：1048 - 1057.

［17］謝花林. 環(huán)鄱陽(yáng)湖地區(qū)農(nóng)業(yè)經(jīng)濟(jì)空間差異分析——基于探索性空間數(shù)據(jù)分析（ESDA）方法［J］.農(nóng)業(yè)現(xiàn)代化研究，2010，31（3）：299 - 303.

［18］彭程，陳志芬，吳華瑞，等. 基于ESDA的城市可持續(xù)發(fā)展能力時(shí)空分異格局研究［J］.中國(guó)人口·資源與環(huán)境，2016，26（2）：144 - 151.

［19］楊楊，吳次芳，韋仕川，等. 中國(guó)建設(shè)用地?cái)U(kuò)展的空間動(dòng)態(tài)演變格局：基于EBI和EBIi的研究［J］.中國(guó)土地科學(xué)，2008，22（1）：23 - 31.

［20］蒲英霞，葛瑩，馬榮華，等. 基于 ESDA的區(qū)域經(jīng)濟(jì)空間差異分析——以江蘇省為例［J］.地理研究，2005，24（6）：965 - 974.

［21］鮑洪杰，劉德光．甘肅省區(qū)域經(jīng)濟(jì)空間差異分析研究：基于Geado095i. 統(tǒng)計(jì)分析［J］.工業(yè)技術(shù)經(jīng)濟(jì)，2011，30（9）：54 - 59.

［22］王雪青，陳媛，劉炳勝. 中國(guó)區(qū)域房地產(chǎn)經(jīng)濟(jì)發(fā)展水平空間統(tǒng)計(jì)分析［J］.數(shù)理統(tǒng)計(jì)與管理，2014，33（1）：59 - 71.

［23］雷廣海，劉友兆，陸效平. 江蘇省13城市土地利用集約度時(shí)空變異及驅(qū)動(dòng)因素［J］.長(zhǎng)江流域資源與環(huán)境，2009，18（1）：7 - 13.

［24］孔偉，郭杰，歐名豪. 不同經(jīng)濟(jì)發(fā)展水平下的建設(shè)用地集約利用及區(qū)域差別化管控［J］.中國(guó)人口·資源與環(huán)境，2014，24（4）：100 - 106.