基于ESDA的中國(guó)城鄉(xiāng)收入差距時(shí)空分異格局分析

孫煥芳，梁 敏

（1.國(guó)家測(cè)繪地理信息局第一航測(cè)遙感院，陜西 西安 710054）

隨著我國(guó)經(jīng)濟(jì)的高速增長(zhǎng)，特別是改革開放以后，城鄉(xiāng)收入差距不斷擴(kuò)大，已成為影響我國(guó)社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展與全面建設(shè)小康社會(huì)的熱點(diǎn)問題。同時(shí)，城鄉(xiāng)收入差距也引起了學(xué)術(shù)界的廣泛關(guān)注，有關(guān)專家針對(duì)城鄉(xiāng)收入差距產(chǎn)生的根源[1-2]、影響因素[3-4]、變化趨勢(shì)[5-6]、解決對(duì)策以及與城市化的關(guān)系[7]、縮小城鄉(xiāng)收入差距的對(duì)策等進(jìn)行了廣泛討論[8]。但是，目前國(guó)內(nèi)對(duì)城鄉(xiāng)收入差距空間分布的相關(guān)性、異質(zhì)性的研究較缺乏，本文以2002～2014年中國(guó)339個(gè)地級(jí)及以上城市為研究對(duì)象，運(yùn)用空間自相關(guān)性、空間異質(zhì)性分析方法分析了城鄉(xiāng)收入的空間格局，為國(guó)家因地制宜實(shí)行差異化、區(qū)域化的發(fā)展政策提供科學(xué)參考。

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

1.1 研究方法

1.1.1 空間自相關(guān)

空間自相關(guān)是描述研究區(qū)域?qū)傩灾蹬c其相鄰位置屬性值是否存在相關(guān)性以及相關(guān)程度的一種空間數(shù)據(jù)分析方法。

1）全局空間自相關(guān)，是對(duì)屬性值在整個(gè)研究區(qū)域的空間特征的描述，主要通過Moran'sI指數(shù)來衡量，其計(jì)算公式為[9-10]：

式中，xi、xj為空間單元i、j的屬性值；n為空間單元的數(shù)目；Wij為i、j之間的空間連接矩陣。本文采用基于鄰近標(biāo)準(zhǔn)的Queen Contiguity方法建立空間權(quán)重文件，即若i與j相鄰，則Wij=1，否則Wij=0；當(dāng)i=j時(shí)，Wij=0

對(duì)Moran'sI計(jì)算結(jié)果進(jìn)行Z統(tǒng)計(jì)檢驗(yàn)，其檢驗(yàn)公式為：

2）局部空間自相關(guān)。為了更好地反映局部空間單元的空間自相關(guān)性質(zhì)，本文采用LISA指數(shù)度量各區(qū)域與周邊地區(qū)之間的局部空間關(guān)聯(lián)和空間差異程度[11-12]。其計(jì)算公式為：

式中，Zi、Zj為屬性值的均值標(biāo)準(zhǔn)化。

1.1.2 全局趨勢(shì)分析

通常，一個(gè)表面由確定的全局趨勢(shì)和隨機(jī)的短程變異組成。全局趨勢(shì)反映空間物體在空間區(qū)域上變化的主體特征，主要揭示空間物體的總體規(guī)律。本文利用ArcGIS的全局趨勢(shì)分析功能，將屬性值賦予每個(gè)空間單元，形成三維趨勢(shì)面，從而分析數(shù)據(jù)在特定方向上的趨勢(shì)。

1.1.3 空間變差函數(shù)

長(zhǎng)江委2018年工作會(huì)議提出了“由開發(fā)中保護(hù)向保護(hù)中發(fā)展轉(zhuǎn)變、實(shí)現(xiàn)以保護(hù)為前提的高質(zhì)量可持續(xù)發(fā)展”的治江新思路，確立了“四個(gè)長(zhǎng)江”的治江目標(biāo)，擬通過補(bǔ)齊補(bǔ)強(qiáng)防洪抗旱減災(zāi)體系短板，著力構(gòu)建安瀾長(zhǎng)江；強(qiáng)化優(yōu)化水資源配置和節(jié)約利用，著力構(gòu)建綠色長(zhǎng)江；建立健全流域綜合管理體制機(jī)制，著力構(gòu)建和諧長(zhǎng)江；加快推進(jìn)水生態(tài)環(huán)境保護(hù)與修復(fù)，著力構(gòu)建美麗長(zhǎng)江。

空間變差函數(shù)亦稱半變異函數(shù)，是描述區(qū)域化變量隨機(jī)性和結(jié)構(gòu)性的基本手段，設(shè)區(qū)域化變量Z(xi)和Z(xi+h)分別為Z(x)在空間位置xi和xi+h上的觀測(cè)值，則空間變差函數(shù)的計(jì)算公式為[13]：

式中，N(h)為分割距離為h的樣本量。

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

本文以2002、2014年兩個(gè)年份作為研究斷面，以城鄉(xiāng)收入比（城鎮(zhèn)居民可支配收入/農(nóng)民人均純收入）為研究指標(biāo)，以地級(jí)及以上城市（包括地級(jí)市、地區(qū)、自治州和盟，除港、澳、臺(tái)）為研究單元。數(shù)據(jù)來源于中國(guó)統(tǒng)計(jì)信息網(wǎng)、各省市統(tǒng)計(jì)年鑒和中國(guó)城市統(tǒng)計(jì)年鑒等。

2 時(shí)空差異分析

2.1 全局空間格局分析

為了更直觀地表示中國(guó)城鄉(xiāng)收入差距的時(shí)空分布，根據(jù)中國(guó)實(shí)際情況和國(guó)家經(jīng)驗(yàn)，將城鄉(xiāng)收入比重劃分為4類，分別為低收入比（<1.5），中收入比（1.5～2.5），高收入比（2.5～3.5）和極高收入比（>3.5）。基于此，制作表1說明中國(guó)2002～2014年間城鄉(xiāng)收入差距的時(shí)間序列變化情況。

表1 2002～2014年中國(guó)高、中、低城鄉(xiāng)收入比

由表1可知，13 a間中國(guó)城鄉(xiāng)收入差距先增后減；雖從2006年開始城鄉(xiāng)收入差距有縮小趨勢(shì)，但整體上中國(guó)城鄉(xiāng)收入差距還是比較大的。高收入比的地區(qū)一直占據(jù)主導(dǎo)地位，2002年高收入比的地區(qū)有215個(gè)，占地區(qū)總數(shù)的63.05%，2006年有251個(gè)，占地區(qū)總數(shù)的73.60%，2010年有209個(gè)，占地區(qū)總數(shù)的61.29%，2014年有182個(gè)，占地區(qū)總數(shù)的53.66%；低收入比的地區(qū)很少，2002年有1個(gè)（丹東市），2010年有4個(gè)，全部位于黑龍江省（丹東市、鶴崗市、伊春市、綏化市），2014年有3個(gè)也全部位于黑龍江?。档そ小⒁链菏?、雞西市）；總的來說，中收入比地區(qū)的個(gè)數(shù)有所增加，2002年中收入比的地區(qū)有125個(gè)，2014年達(dá)到156個(gè)。

由表2可以看出，中國(guó)城鄉(xiāng)居民收入比大致呈東部較小、西部最大，由東向西逐步擴(kuò)大的趨勢(shì)，其地域分布與我國(guó)經(jīng)濟(jì)發(fā)展水平三大地帶的劃分基本一致。城鄉(xiāng)收入比大于3.5的地區(qū)全部位于西部地區(qū)（除2010年的海南省外），其中西藏自治區(qū)是收入差距最大的地區(qū)，歷年城鄉(xiāng)收入比均大于3，2002年城鄉(xiāng)收入比更是高達(dá)5。這主要是因?yàn)槲鞑刈灾螀^(qū)當(dāng)時(shí)正處于經(jīng)濟(jì)相對(duì)落后的工業(yè)化起步階段，國(guó)家對(duì)其工業(yè)化發(fā)展的支持力度較大，如加大財(cái)政轉(zhuǎn)移支付的力度，為了鼓勵(lì)優(yōu)秀人才在西藏建功立業(yè)而采取的相對(duì)較高的工資制度等；但相比而言，廣大的農(nóng)牧民族地區(qū)，由于當(dāng)?shù)刈匀粭l件惡劣，市場(chǎng)經(jīng)濟(jì)落后等，導(dǎo)致農(nóng)村經(jīng)濟(jì)也很落后，農(nóng)民收入不僅低而且很難提高。

表2 2002～2014年中國(guó)各省區(qū)城鄉(xiāng)收入比重

另外，從表2還可以看出，城鄉(xiāng)收入比顯著縮小的省區(qū)包括新疆維吾爾自治區(qū)、西藏自治區(qū)、青海省、云南省、四川省和甘肅省，這些地區(qū)全部位于中國(guó)西部，說明中國(guó)城鄉(xiāng)收入差距顯著縮小的地區(qū)主要集中在西部，這與中國(guó)實(shí)施的西部大開發(fā)戰(zhàn)略密切相關(guān)。城鄉(xiāng)收入比增加的省區(qū)包括寧夏回族自治區(qū)、福建省等，其中福建省城鄉(xiāng)收入比增加最明顯，由2002年的2.4增加到2014年的2.8。東部、中部地區(qū)在2006年以前，城鄉(xiāng)收入比基本上呈增加趨勢(shì)，2006年以后逐漸呈減小趨勢(shì)，說明城鄉(xiāng)收入比在2006年達(dá)到最大值，隨后開始逐漸縮小，這與中國(guó)在2006年取消農(nóng)業(yè)稅和大部分的農(nóng)村收費(fèi)（如糧食補(bǔ)貼、農(nóng)資綜合直補(bǔ)以及農(nóng)業(yè)牲畜養(yǎng)殖補(bǔ)貼等）密切相關(guān)，這些政策在一定程度上減輕了農(nóng)民負(fù)擔(dān)，增加了農(nóng)民收入；同時(shí)社會(huì)主義市場(chǎng)經(jīng)濟(jì)體制改革也為農(nóng)村富余勞動(dòng)力提供了大量的就業(yè)崗位。

2.2 空間自相關(guān)分析

2.2.1 全局空間自相關(guān)分析

本文利用OpenGeoda軟件分別計(jì)算了中國(guó)339個(gè)地級(jí)及以上城市2002～2014年的城鄉(xiāng)居民收入比的全局Moran'sI指數(shù)和相關(guān)參數(shù)。由表3可知，4個(gè)年份的Moran'sI指數(shù)全部為正，且檢驗(yàn)結(jié)果顯著（Z值在顯著性水平P=0.001下明顯超過臨界值1.96）。

表3 2002～2014年中國(guó)城鄉(xiāng)居民收入比的Moran's I 值

根據(jù)Moran'sI指數(shù)來看，4個(gè)年份的城鄉(xiāng)收入差距在空間上的集聚程度均很明顯；由Moran'sI指數(shù)的變化情況可知，隨著時(shí)間的推移，中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距的空間依賴性逐漸減弱。

2.2.2 局部空間自相關(guān)分析

為了揭示中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距的熱點(diǎn)區(qū)和冷點(diǎn)區(qū)，本文利用Geoda軟件計(jì)算和繪制了各地區(qū)城鄉(xiāng)收入比的LISA集聚圖。由圖1可知，中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距發(fā)展格局具有明顯的動(dòng)態(tài)特點(diǎn)。

2002年城鄉(xiāng)收入差距呈“高—高”集聚格局的地區(qū)主要分布在西部地區(qū)，且連片分布，特別是西藏自治區(qū)、青海省、云南省、廣西壯族自治區(qū)和甘肅省的東南部、陜西省的北部等地區(qū)；而呈“低—低”集聚格局的地區(qū)分布較零散，主要集中在東部地區(qū)，特別是遼寧省、江蘇省和廣東省的南部、山東省的北部等地區(qū)。與2002年相比，中國(guó)2006年城鄉(xiāng)收入差距呈“高—高”，“低—低”集聚格局的地區(qū)都有所減少，呈“高—高”集聚格局的地區(qū)有45個(gè)，比2002年減少了8個(gè)，約占地區(qū)總數(shù)的13.27%，呈“低—低”集聚格局的地區(qū)有40個(gè)，比2002年減少了11個(gè)，約占地區(qū)總數(shù)的11.80%。

圖1 2002～2014年中國(guó)城鄉(xiāng)收入差距LISA集聚圖（審圖號(hào)：GS(2013)1655）

2010年，呈“高—高”、“低—低”集聚格局的地區(qū)有所增加，分別占地區(qū)總數(shù)的19.18%和14.45%，前者向西部地區(qū)、西南地區(qū)、黃河中上游地區(qū)集聚，后者主要集中在中國(guó)的東北部、新疆北部地區(qū)，另外東部和南部也有零散分布。與2010年相比，2014年呈“高—高”、“低—低”集聚格局的地區(qū)均有所減少，其中呈“低—低”集聚格局的地區(qū)減少幅度較大，由2010年的49個(gè)減少到2014年的20個(gè)?？傮w而言，4個(gè)年份中呈“高—高”、“低—低”集聚格局的地區(qū)占地區(qū)總數(shù)的比例一直大于20%，分別為32.15%、27.43%、34.51%和23.01%，說明中國(guó)城鄉(xiāng)收入差距呈不均衡的發(fā)展格局，而這種空間格局的形成與各地區(qū)的區(qū)域條件、資源稟賦以及國(guó)家政策有關(guān)。

2.3 熱點(diǎn)分析

LISA集聚圖雖能反映中國(guó)城鄉(xiāng)收入差距的局部空間相關(guān)性，但部分地區(qū)的相關(guān)性不明顯。為了彌補(bǔ)和驗(yàn)證LISA集聚圖的分析結(jié)果，本文引入G*統(tǒng)計(jì)指數(shù)對(duì)2002年、2004年、2010年和2014年的城鄉(xiāng)收入差距進(jìn)行了熱點(diǎn)分析，并按照G*指數(shù)的大小，從高到低劃分為熱點(diǎn)地區(qū)、次熱地區(qū)、次冷地區(qū)和冷點(diǎn)地區(qū)4 個(gè)類型，見圖2。

圖2 2002～2014年中國(guó)城鄉(xiāng)收入差距熱點(diǎn)分析圖（審圖號(hào)：GS(2013)1655）

由圖2可知，中國(guó)城鄉(xiāng)收入差距呈南北方向上的“帶狀”分布，熱點(diǎn)地區(qū)主要分布在西部和中部地區(qū)，并逐漸向中部地區(qū)轉(zhuǎn)移；次熱地區(qū)緊鄰熱點(diǎn)地區(qū)分布；冷點(diǎn)地區(qū)主要分布在東部和東北地區(qū)，并逐漸向東部地區(qū)轉(zhuǎn)移；次冷地區(qū)緊鄰冷點(diǎn)地區(qū)分布。空間分析結(jié)果與LISA集聚圖的分析結(jié)果是一致的。

3 結(jié) 語

本文以2002～2014年中國(guó)339個(gè)地級(jí)及以上城市的城鄉(xiāng)居民收入比為衡量城鄉(xiāng)收入差距的基礎(chǔ)數(shù)據(jù)，采用Moran'sI和空間變差函數(shù)探討了4個(gè)年份中國(guó)城鄉(xiāng)收入差距的空間格局。結(jié)果表明，中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距具有以下特點(diǎn)：①2002～2014年中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距逐漸減??；②中國(guó)城

鄉(xiāng)收入差距在空間上存在一定的正空間自相關(guān)性，即城鄉(xiāng)居民收入比相似的地區(qū)在空間上積聚分布；③根據(jù)LISA集聚圖和熱點(diǎn)圖分析可知，中國(guó)西部和中部地區(qū)的城鄉(xiāng)收入差距較大，東部地區(qū)的收入差距相對(duì)較小，呈西部>中部>東部的分布格局。本文仍存在一定的缺陷，所選的指標(biāo)比較單一，時(shí)間跨度小，不足以代表13 a以來的發(fā)展趨勢(shì)。引入其他指標(biāo)，并對(duì)其他年份進(jìn)行分析，得出的結(jié)果是否會(huì)有所差異，有待進(jìn)一步探討，多指標(biāo)、長(zhǎng)時(shí)間、小區(qū)域是今后繼續(xù)研究的方向。

[1] 高帆.中國(guó)居民收入差距變動(dòng)的因素分解:趨勢(shì)及解釋[J].經(jīng)濟(jì)科學(xué),2012(3):5-17

[2] 韓建雨.我國(guó)城鄉(xiāng)居民收入差距變動(dòng)的影響因素分析[J].統(tǒng)計(jì)與決策,2013(9):126-128

[3] 趙人偉,李實(shí).中國(guó)居民收入差距的擴(kuò)大及其原因[J].經(jīng)濟(jì)研究,1997(9):19-28

[4] 程開明,李金昌.城市偏向、城市化與城鄉(xiāng)收入差距的作用機(jī)制及動(dòng)態(tài)分析[J].數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究,2007(7):116-125

[5] 徐永發(fā).改革開放以來我國(guó)城鄉(xiāng)居民收入差距問題研究[D].武漢:武漢大學(xué),2005:25-29

[6] 陳斌開,林毅夫.發(fā)展戰(zhàn)略、城市化與中國(guó)城鄉(xiāng)收入差距[J].中國(guó)社會(huì)科學(xué),2013(4):81-102

[7] 韓其恒,李俊青.二元經(jīng)濟(jì)下的中國(guó)城鄉(xiāng)收入差距的動(dòng)態(tài)演化研究[J].金融研究,2011(8):15-30

[8] 張義博,劉文忻.人口流動(dòng)、財(cái)政支出結(jié)構(gòu)與城鄉(xiāng)收入差距[J].中國(guó)農(nóng)村經(jīng)濟(jì),2012(1):16-30

[9] 潘竟虎.中國(guó)地級(jí)及以上城市城鄉(xiāng)收入差距時(shí)空分異格局[J].經(jīng)濟(jì)地理,2014(6):60-67

[10] 王靜,張小雷,杜宏茹.新疆縣域經(jīng)濟(jì)空間格局演化特征[J].地理科學(xué)進(jìn)展,2011(4):471-478

[11] 何格,歐名豪.城鎮(zhèn)建設(shè)用地效率時(shí)空分異研究:以青島為例[J].商業(yè)研究,2007(10):123-127

[12] 方葉林,黃震方,涂瑋,等.基于地統(tǒng)計(jì)分析的安徽縣域經(jīng)濟(jì)空間差異研究[J].經(jīng)濟(jì)地理,2013(2):33-38

[13] 周揚(yáng),李寧,吳文祥,等.1982～2010年中國(guó)縣域經(jīng)濟(jì)發(fā)展時(shí)空格局演變[J].地理科學(xué)進(jìn)展,2014(1):102-113

[14] 方葉林,黃震方,陳文娣,等.2001～2010年安徽省縣域經(jīng)濟(jì)空間演化[J].地理科學(xué)進(jìn)展,2013(5):831-839

[15] 陳鵬,李錦濤,馬偉.犯罪熱點(diǎn)的分析方法研究[J].中國(guó)人民公安大學(xué)學(xué)報(bào)(自然科學(xué)版),2012(3):53-57