范 昕, 李江風(fēng), 陳萬旭, 朱麗君, 韋啟杰, 林正松

(1.中國地質(zhì)大學(xué)(武漢) 公共管理學(xué)院, 武漢 430074; 2.中國地質(zhì)大學(xué)(武漢)國土資源部法律評價工程重點實驗室, 武漢 430074； 3.中國地質(zhì)大學(xué)(武漢) 環(huán)境學(xué)院, 武漢 430074)

1978年中國城鎮(zhèn)化率低于20%，到2014年已超過50%，快速的城鎮(zhèn)化有力支撐了中國經(jīng)濟(jì)社會的高速增長和快速轉(zhuǎn)型[1]。與此同時隨著城鎮(zhèn)化的推進(jìn)，我國出現(xiàn)了城鄉(xiāng)發(fā)展失衡、分配失衡、社會權(quán)利失衡、國土空間開發(fā)利用效率低下、經(jīng)濟(jì)發(fā)展模式粗放、“空城”或“鬼城”等一系列問題和矛盾[3-4]，進(jìn)而導(dǎo)致不同類型的城鎮(zhèn)化失衡[5]、城鎮(zhèn)化協(xié)調(diào)性發(fā)展不均衡[6-7]、城鎮(zhèn)化內(nèi)部不協(xié)調(diào)等現(xiàn)象[3,8-9]，限制了城鎮(zhèn)化健康發(fā)展。這些問題的出現(xiàn)，說明我國城鎮(zhèn)化已步入城鎮(zhèn)化轉(zhuǎn)型關(guān)鍵時期[2]，并且未來5年仍是我國城鎮(zhèn)化快速發(fā)展時期[10]，為此解決城鎮(zhèn)化協(xié)調(diào)發(fā)展問題，促進(jìn)城鎮(zhèn)化順利轉(zhuǎn)型，確保我國城鎮(zhèn)化健康發(fā)展顯得尤為重要。

近年來國內(nèi)外學(xué)者對城鎮(zhèn)化進(jìn)行了大量的研究，取得了豐富的成果。部分學(xué)者對城鎮(zhèn)化的概念、內(nèi)涵、機(jī)理、演進(jìn)軌跡進(jìn)行了完善擴(kuò)充[8,11-12]，提出城鎮(zhèn)化是一個綜合過程[13]，包括人口城鎮(zhèn)化、社會城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化、土地城鎮(zhèn)化等多方面的內(nèi)涵[14]；一些學(xué)者構(gòu)建了評價指標(biāo)體系對我國城鎮(zhèn)化發(fā)展水平進(jìn)行測度[15-17]，對我國城鎮(zhèn)化發(fā)展時空特征及驅(qū)動因素進(jìn)行研究；同時一些學(xué)者還對城鎮(zhèn)化與新型工業(yè)化、信息化、農(nóng)業(yè)現(xiàn)代化的耦合協(xié)調(diào)關(guān)系進(jìn)行研究，也不乏“四化”系統(tǒng)中某兩個或者三個子系統(tǒng)之間融合協(xié)調(diào)發(fā)展、互動關(guān)系及作用機(jī)理的研究[18-22]。較多的研究集中在我國城鎮(zhèn)化系統(tǒng)中各子系統(tǒng)的相互協(xié)調(diào)性[23]，其中包括兩個[24-25]、三個[26]或者多個系統(tǒng)的城鎮(zhèn)化協(xié)調(diào)性的研究[27]。研究尺度涉及全國、省域、市域、縣域、城市群和都市圈等各個層次；在研究方法上，指標(biāo)權(quán)重的確定大多采用專家打分法、層次分析法等主觀賦權(quán)法和熵值法、主成分分析法等客觀賦權(quán)法，城鎮(zhèn)化發(fā)展水平的測度主要采用TOPSIS法、變異系數(shù)法、灰色關(guān)聯(lián)度法、投影尋蹤模型、隸屬函數(shù)法、地理探測器、空間距離測度模型、DEA、GWR和ESDA等方法進(jìn)行定量分析，對于城鎮(zhèn)化各子系統(tǒng)協(xié)調(diào)性的測度主要采用耦合協(xié)調(diào)度模型。

本文在相關(guān)研究成果基礎(chǔ)上，對城鎮(zhèn)化協(xié)調(diào)發(fā)展理論進(jìn)一步創(chuàng)新，構(gòu)建城鎮(zhèn)化協(xié)調(diào)發(fā)展的三維模型，并利用該模型對2000年、2005年、2010年、2014年河南省108個縣市城鎮(zhèn)化系統(tǒng)中的人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化以及土地城鎮(zhèn)化3個子系統(tǒng)的協(xié)調(diào)發(fā)展水平進(jìn)行測度，采用ESDA和GWR模型對河南省縣域城鎮(zhèn)化協(xié)調(diào)發(fā)展時空分異及驅(qū)動因素進(jìn)行分析，以期為河南省制定縣域城鎮(zhèn)化健康發(fā)展路徑提供借鑒。

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

1.1 研究區(qū)概況

河南省位于中國中東部，黃河中下游，110°21′—116°39′E，31°23′—36°22′N，包含17個地級市，1個省直管縣級市，是我國的人口大省、產(chǎn)糧大省和經(jīng)濟(jì)大省，人均資源量排名靠后。2014年河南省以占全國1.74%的土地面積，承載了全國6.90%的人口，生產(chǎn)了全國糧食產(chǎn)量的9.50%，經(jīng)濟(jì)總量居于全國第五位，中西部第一位；經(jīng)濟(jì)發(fā)展主要依托農(nóng)業(yè)，一、二、三產(chǎn)業(yè)增加值比重分別為5.3%，62.4%，32.3%，產(chǎn)業(yè)結(jié)構(gòu)不合理；城鎮(zhèn)化水平為45.20%，與全國城鎮(zhèn)化水平54.77%相去甚遠(yuǎn)；2000年河南省建成區(qū)面積為785 km2，到2014年建成區(qū)面積為1 877 km2，增加了1.39倍。河南省經(jīng)濟(jì)區(qū)位和交通區(qū)位優(yōu)勢明顯，擁有水、陸、空綜合立體交通運輸體系，起著承東啟西、連南貫北、吸引四面、輻射八方的作用，河南省城鎮(zhèn)化質(zhì)量對于中部地區(qū)城鎮(zhèn)化乃至全國城鎮(zhèn)化健康發(fā)展有著重要意義。

1.2 數(shù)據(jù)來源與數(shù)據(jù)處理

本研究涉及河南省108個縣域，研究時間選取2000年、2005年、2010年、2014年4個節(jié)點，跨越3個五年計劃，在所選的研究期中包含了中國城鎮(zhèn)化發(fā)展最快速的10年。本研究涉及的基礎(chǔ)數(shù)據(jù)主要來源于2001年、2006年、2011年、2015年《河南省統(tǒng)計年鑒》、《中華人民共和國分縣人口統(tǒng)計資料》及河南省18個省轄市相關(guān)年份統(tǒng)計年鑒及權(quán)威網(wǎng)站，其中一些指標(biāo)是經(jīng)過相關(guān)指標(biāo)計算得出，數(shù)據(jù)缺失年份取相鄰年份均值或采用插值法獲取。研究中總?cè)丝诓捎玫氖强h域常住人口，城鎮(zhèn)化率采用非農(nóng)人口比重來衡量，非農(nóng)人口是指從事第二、第三產(chǎn)業(yè)的人口及其所撫養(yǎng)的部分人口。把非農(nóng)人口比重作為衡量城鎮(zhèn)化水平的指標(biāo)，是從農(nóng)民真正實現(xiàn)進(jìn)城落戶、安居樂業(yè)的角度出發(fā)，來研究河南省縣域城鎮(zhèn)化水平的時空差異。由于各指標(biāo)統(tǒng)計口徑的量綱不一致，本研究采用標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化法對原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。各行政單元的行政界線、城市駐地、DEM數(shù)據(jù)等位置信息數(shù)據(jù)來自國家基礎(chǔ)地理信息中心1∶400萬數(shù)據(jù)庫。

2 指標(biāo)構(gòu)建與研究方法

2.1 指標(biāo)構(gòu)建

本文在參照有關(guān)專家學(xué)者的研究成果基礎(chǔ)上[7,24,27-28]，遵循科學(xué)、可靠、層次性、可操作性的原則，圍繞可以集中體現(xiàn)河南省城鎮(zhèn)化水平的人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化、土地城鎮(zhèn)化3個子系統(tǒng)的內(nèi)涵來構(gòu)建指標(biāo)體系，通過對城鎮(zhèn)化的這3大子系統(tǒng)進(jìn)行測度分析得到河南省城鎮(zhèn)化發(fā)展水平以及協(xié)調(diào)發(fā)展度狀況。具體的指標(biāo)構(gòu)建見表1。

表1 城鎮(zhèn)化子系統(tǒng)評價指標(biāo)體系

2.2 研究方法

2.2.1 城鎮(zhèn)化發(fā)展水平和協(xié)調(diào)發(fā)展度 人口城鎮(zhèn)化子系統(tǒng)、經(jīng)濟(jì)城鎮(zhèn)化子系統(tǒng)和土地城鎮(zhèn)化子系統(tǒng)是城鎮(zhèn)化系統(tǒng)中3個重要的組成部分，其中人口城鎮(zhèn)化是城鎮(zhèn)化的核心，經(jīng)濟(jì)城鎮(zhèn)化是動力，土地城鎮(zhèn)化是保障[28]，三者相互制約、相互影響、相互作用、相互促進(jìn)，形成互為動力、互為支撐、良性諧振、良性循環(huán)的城鎮(zhèn)化發(fā)展綜合體，任何一方面的偏頗都將降低城鎮(zhèn)化綜合效益。

首先構(gòu)建人口城鎮(zhèn)化發(fā)展指數(shù)Pp、經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)Ee、土地城鎮(zhèn)化發(fā)展指數(shù)Ll，用上述3個指數(shù)測度城鎮(zhèn)化子系統(tǒng)的發(fā)展水平：

(1)

圖1城鎮(zhèn)化協(xié)調(diào)發(fā)展理論的基本模式

Ci=cosαi

(2)

式中：Ci為協(xié)調(diào)發(fā)展系數(shù)；Di為綜合發(fā)展度；Xi，Yi和Zi分別表示各縣市的人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化和土地城鎮(zhèn)化發(fā)展指數(shù)；αi為各區(qū)域的實際發(fā)展行為向量與最佳發(fā)展行為向量之間的夾角。

為進(jìn)一步衡量區(qū)域人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化與土地城鎮(zhèn)化之間的協(xié)調(diào)發(fā)展?fàn)顩r，本文采用協(xié)調(diào)發(fā)展度(Ti)作為衡量尺度[7]，具體公式如下：

(3)

圖2改進(jìn)后的城鎮(zhèn)化協(xié)調(diào)發(fā)展理論模型

2.2.2 城鎮(zhèn)化區(qū)域差異 采用泰爾指數(shù)(Theil Index)對河南省地級市內(nèi)部以及地級市的城鎮(zhèn)化水平差異進(jìn)行測度，將城鎮(zhèn)化水平的空間差異(T)分解為市內(nèi)差異(Tw)和市際差異(Tb)。泰爾指數(shù)越大，表示區(qū)域之間差異越大，反之則越小[31]。

(4)

(5)

(6)

2.2.3 探索性空間分析 為了研究河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度是否在空間上表現(xiàn)出分布的聚集性和空間異常等特征，引入空間探索性分析手段ESDA進(jìn)行分析，利用全局Moran′sI指數(shù)和局域Moran′sI指數(shù)(LISA)對河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度空間關(guān)聯(lián)模式(趨同或異質(zhì))進(jìn)行度量和檢驗，揭示其空間分布特征。其中全局空間自相關(guān)用來驗證區(qū)域的某種屬性值在整個區(qū)域空間分布的集聚態(tài)勢；局部空間自相關(guān)(Local Moran′sI)是度量研究單元某一屬性值與臨近地區(qū)相同屬性值的空間差異程度，具體公式參照相關(guān)文獻(xiàn)[17]。根據(jù)Local Moran′sI指數(shù)，計算結(jié)果出現(xiàn)High—High，Low—Low，Low—High，High—Low共4種類型，分別表示相鄰空間的關(guān)系為高高、低低、低高、高低區(qū)。

2.2.4 城鎮(zhèn)化協(xié)調(diào)發(fā)展影響因素 引入地理加權(quán)回歸模型對河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展區(qū)域差異及其驅(qū)動因子進(jìn)行研究，該模型可以很好地納入數(shù)據(jù)的空間特性，客觀實際地探測空間數(shù)據(jù)的空間非平穩(wěn)性。應(yīng)用加權(quán)最小二乘法(WLS)對每個觀測點的參數(shù)進(jìn)行估計，加權(quán)時一般采用高斯函數(shù)構(gòu)建加權(quán)函數(shù)，利用(赤池信息準(zhǔn)則)[32]來確定最優(yōu)帶寬，然后構(gòu)建河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展的GWR模型進(jìn)行模擬，并根據(jù)模擬和模型檢驗結(jié)果進(jìn)行分析討論。其模型公式如下：

(7)

式中:Yi為被解釋變量值；(Si，Ti)是第i個觀測點空間地理位置坐標(biāo)；αj(Si，Ti)為觀測點回歸參數(shù)；xij為觀解釋變量值；εi是隨機(jī)誤差。

3 結(jié)果與分析

在構(gòu)建人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化和土地城鎮(zhèn)化評價指標(biāo)體系的基礎(chǔ)上，利用公式(1)、(2)、(3)測算出河南省各縣市城鎮(zhèn)化子系統(tǒng)發(fā)展水平、協(xié)調(diào)發(fā)展度以及綜合發(fā)展度(圖3)。

圖3 河南省城鎮(zhèn)化平均發(fā)展指數(shù)、協(xié)調(diào)發(fā)展度與綜合發(fā)展度

3.1 河南省城鎮(zhèn)化發(fā)展水平與協(xié)調(diào)性的時空演變特征分析

3.1.1 河南省城鎮(zhèn)化發(fā)展水平分析 研究期間城鎮(zhèn)化子系統(tǒng)平均發(fā)展水平逐年提高，人口城鎮(zhèn)化平均發(fā)展水平由2000年的0.178 0增加到2014年的0.468 1，經(jīng)濟(jì)城鎮(zhèn)化平均發(fā)展水平由2000年的0.202 6增加到2014年的0.475 8，土地城鎮(zhèn)化平均發(fā)展水平由2000年的0.207 6增加到2014年的0.286 6。通過對城鎮(zhèn)化各子系統(tǒng)發(fā)展水平進(jìn)行對比發(fā)現(xiàn)，2000年河南省土地城鎮(zhèn)化平均發(fā)展水平優(yōu)于經(jīng)濟(jì)城鎮(zhèn)化和人口城鎮(zhèn)化平均發(fā)展水平，2005年、2010年、2014年經(jīng)濟(jì)城鎮(zhèn)化平均發(fā)展水平高于人口城鎮(zhèn)化平均發(fā)展水平，土地城鎮(zhèn)化處于最低水平,說明2000—2014年河南省城鎮(zhèn)化歷經(jīng)了經(jīng)濟(jì)人口城鎮(zhèn)化滯后型到土地城鎮(zhèn)化滯后型的過程，即先經(jīng)歷建設(shè)用地快速擴(kuò)張后再有經(jīng)濟(jì)的高速發(fā)展以及城鎮(zhèn)人口的迅速增長，從側(cè)面揭示出了河南省的城鎮(zhèn)土地利用從粗放利用逐步向集約利用方向發(fā)展的趨勢。河南省東部地區(qū)及中原城市群內(nèi)部地區(qū)的人口城鎮(zhèn)化發(fā)展水平相對較高，另外河南省省轄市周邊地區(qū)人口城鎮(zhèn)化水平顯著高于其他縣市；經(jīng)濟(jì)城鎮(zhèn)化發(fā)展水平較低的地區(qū)主要集中于周口市、商丘市以及河南省黃河灘區(qū)，鄭州市周邊縣市經(jīng)濟(jì)城鎮(zhèn)化水平明顯高于其他地區(qū)；土地城鎮(zhèn)化發(fā)展水平較低的地區(qū)主要集中于豫西太行山深石山區(qū)、西南部伏牛山區(qū)以及東南部大別山區(qū)，較高的地區(qū)是義馬市與長葛市，并且平原地區(qū)土地城鎮(zhèn)化發(fā)展指數(shù)明顯高于山地丘陵地區(qū)，說明河南省城鎮(zhèn)化3個子系統(tǒng)發(fā)展水平空間分異明顯，地形地貌和地理區(qū)位條件對城鎮(zhèn)化各子系統(tǒng)的發(fā)展水平有著顯著的影響。

河南省城鎮(zhèn)化綜合發(fā)展度平均值從2000年的0.196 1增加到2014年的0.410 2，增長趨勢明顯，2000年大部分縣市綜合發(fā)展度處于較低水平，其中47個縣市發(fā)展水平高于河南省平均發(fā)展水平，2014年有45個縣市發(fā)展水平高于河南省平均發(fā)展水平，從中可以看出2000—2014年河南省城鎮(zhèn)化整體發(fā)展水平較快。從河南省內(nèi)城鎮(zhèn)化綜合水平地區(qū)差異來看，三門峽的義馬市城鎮(zhèn)化綜合發(fā)展度一直處于最高水平，可能是由于義馬市工業(yè)基礎(chǔ)較好，礦產(chǎn)資源豐富，經(jīng)濟(jì)發(fā)展速度較快，產(chǎn)業(yè)結(jié)構(gòu)相對合理，城鎮(zhèn)化水平較高。另外城鎮(zhèn)化綜合發(fā)展度較高的地區(qū)主要集中于中原城市群內(nèi)部一些縣市，以及一些地級市周邊的縣市，而河南省的西南部山地及丘陵地區(qū)，受地形地貌地理區(qū)位的限制，城鎮(zhèn)化綜合發(fā)展度普遍較低。

3.1.2 河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展分析 河南省城鎮(zhèn)化平均協(xié)調(diào)發(fā)展度如同城鎮(zhèn)化發(fā)展綜合水平以及城鎮(zhèn)化3個子系統(tǒng)發(fā)展水平一樣處于逐年提升的狀態(tài)，從2000年的0.427 6增加到2014年的0.625 1。從各協(xié)調(diào)類型(表2)的縣市數(shù)量來看，2000年只有義馬市1個市的城鎮(zhèn)化發(fā)展水平達(dá)到初級協(xié)調(diào)水平，其他的縣市都處于初級協(xié)調(diào)水平以下，其中有15個縣市處于勉強(qiáng)協(xié)調(diào)水平，53個縣市處于瀕臨失調(diào)水平，39個縣市處于輕度失調(diào)水平；到2014年河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度有很大的提高，原來代表河南省最高協(xié)調(diào)發(fā)展度水平的義馬市由2000年的初級協(xié)調(diào)水平提升到了優(yōu)質(zhì)協(xié)調(diào)水平，此外有9個縣市處于中級協(xié)調(diào)水平，59個縣市處于初級協(xié)調(diào)水平，39個縣市處于勉強(qiáng)協(xié)調(diào)水平?？臻g分布上河南省的城鎮(zhèn)化協(xié)調(diào)發(fā)展度也受地形地貌和地理區(qū)位條件的影響表現(xiàn)出河南省西南部、南部山區(qū)協(xié)調(diào)度低于平原地區(qū)，市區(qū)周邊縣市城鎮(zhèn)化協(xié)調(diào)度高于其他縣市的規(guī)律。

表2 河南省城鎮(zhèn)化協(xié)調(diào)度等級及協(xié)調(diào)發(fā)展類型[4]

3.2 河南省城鎮(zhèn)化發(fā)展水平空間差異分析

研究期間，河南省城鎮(zhèn)化區(qū)域間差異總體上不斷減小，其中區(qū)域間人口城鎮(zhèn)化發(fā)展指數(shù)、經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)、城鎮(zhèn)化綜合發(fā)展度與城鎮(zhèn)化協(xié)調(diào)發(fā)展度總體差異呈下降趨勢，而土地城鎮(zhèn)化發(fā)展指數(shù)在區(qū)域之間差異有逐漸變大的趨勢。經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)的泰爾指數(shù)變化最為明顯，由表3可以看出，由2000年的0.117 6減少到2014年的0.022 6，城鎮(zhèn)化協(xié)調(diào)發(fā)展度的泰爾指數(shù)變化相對較小，從2000年的0.010 0減少到2014年的0.004 3，土地城鎮(zhèn)化發(fā)展指數(shù)從2000年的0.115 8增加到2014年的0.129 8。對比市際泰爾指數(shù)(Tb)和市內(nèi)泰爾指數(shù)(Tw)可以發(fā)現(xiàn)，城鎮(zhèn)化協(xié)調(diào)發(fā)展度、城鎮(zhèn)化綜合發(fā)展度與土地城鎮(zhèn)化發(fā)展指數(shù)的市內(nèi)泰爾指數(shù)大于市際泰爾指數(shù)，經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)的市際泰爾指數(shù)大于市內(nèi)泰爾指數(shù)，2000年、2010年人口城鎮(zhèn)化發(fā)展指數(shù)的市際泰爾指數(shù)高于市內(nèi)泰爾指數(shù)，2005年、2014年則相反。

表3 河南省城鎮(zhèn)化空間總體差異及分維度的泰爾分解

3.3 河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展空間關(guān)聯(lián)性分析

由表4可知，2000年、2005年、2010年、2014年河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展的全局Moran′sI指數(shù)都通過了1%的顯著性檢驗，說明河南省城鎮(zhèn)化整體有較強(qiáng)的空間正相關(guān)，各縣市城鎮(zhèn)化協(xié)調(diào)發(fā)展存在相互關(guān)聯(lián)，在空間上存在集群趨勢。

表4 河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展Moran′s I統(tǒng)計值

注：(1) 采用了K值最臨近空間權(quán)值矩陣，K=5；隨機(jī)性檢驗中，采用9999 permutation；(2)***表示顯著性水平為1%。

河南省各縣市城鎮(zhèn)化協(xié)調(diào)發(fā)展度存在明顯的空間分異格局，主要表現(xiàn)為：高—高區(qū)主要集中于鄭州市、洛陽市、焦作市、平頂山市等地區(qū)，說明中原城市圈內(nèi)部縣市城鎮(zhèn)化發(fā)展水平較高；低—低區(qū)主要分布于伏牛山區(qū)和大別山區(qū)，這些區(qū)域是河南省城鎮(zhèn)化發(fā)展的低洼區(qū)，發(fā)展水平相對較低。對比2000年、2014年的Moran′sI散點圖可以發(fā)現(xiàn)，2000年有75個地區(qū)表現(xiàn)出正向的空間關(guān)系，有29個地區(qū)表現(xiàn)出不相似的空間關(guān)系，其中13個位于第二象限，16個位于第四象限；2014年有79個地級市表現(xiàn)出正向的空間關(guān)系，26個地區(qū)表現(xiàn)出不相似的空間關(guān)系，其中15個位于第二象限，11個位于第四象限，3個跨了第二、三象限。

3.4 河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展影響因素的空間計量分析

基于以上分析，本文選取人口分布密度、人均教育經(jīng)費支出、人均地方財政一般預(yù)算收入、地方財政人均農(nóng)林水利支出、農(nóng)民人均純收入、人均社會消費品零售總額、建筑業(yè)總產(chǎn)值、平均海拔、人均糧食產(chǎn)量9個指標(biāo)進(jìn)行城鎮(zhèn)化協(xié)調(diào)發(fā)展影響因素的空間計量分析。利用SPSS 23.0的描述性統(tǒng)計功能將所選9個指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，然后利用回歸分析方法對標(biāo)準(zhǔn)化后的所有變量進(jìn)行共線性檢驗，避免各指標(biāo)的多重共線性。由檢驗結(jié)果(表5)可知所選變量方差膨脹因子(VIF)都小于10，這9個指標(biāo)都作為GWR模型的解釋變量。本文取2000年、2005年、2010年、2014年的因變量和解釋變量的均值來避免某個時間截面上數(shù)據(jù)誤差影響分析結(jié)果，以所選9個指標(biāo)作為解釋變量，河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度為被解釋變量。

表5 回歸方程自變量的多重共線性檢驗

分別利用ArcGIS軟件提供的最小二乘法(OLS)模型和地理加權(quán)回歸模型(GWR)對影響河南省城鎮(zhèn)化的因素進(jìn)行全局回歸分析，得到的結(jié)果如表6所示:

表6 OLS和GWR回歸分析結(jié)果對比

對比R2可知，GWR的模擬結(jié)果優(yōu)于OLS的結(jié)果。若OLS模型和GWR模型的AICc的差值大于3，表明GWR模型擬合性能優(yōu)于OLS模型[32]，而本文兩者之間的差距為24左右，再次證明GWR的擬合結(jié)果更合理。

從河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度回歸系數(shù)來看，9個解釋變量對于河南省各個縣級單元的參數(shù)估計的結(jié)果均各不相同，有正有負(fù)，存在明顯的空間差異。影響程度由大到小依次為：人口分布密度、人均教育經(jīng)費支出、人均地方財政一般預(yù)算收入、平均海拔、農(nóng)民人均純收入、人均社會消費品零售總額、人均糧食產(chǎn)量、建筑業(yè)總產(chǎn)值、地方財政人均農(nóng)林水利支出。(1) 人口分布密度、人均教育經(jīng)費支出、人均地方財政一般預(yù)算收入對河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響均為正向，三者在一定程度上可以促進(jìn)城鎮(zhèn)化的協(xié)調(diào)發(fā)展。其中人口分布密度的回歸系數(shù)由南向北逐漸增加，人均教育經(jīng)費支出對京廣線以東地區(qū)的影響力明顯高于以西地區(qū)，加大教育經(jīng)費投資力度可以有效地促進(jìn)河南省商丘市、周口市、駐馬店市和信陽等地方城鎮(zhèn)化的進(jìn)程，人均地方財政一般預(yù)算收入對河南省東南地區(qū)的影響強(qiáng)于西南地區(qū)。(2) 農(nóng)民人均純收入對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響有正有負(fù)，沿河南省中部南北一線向東西兩側(cè)逐漸減弱，其中商丘市、濮陽市等部分縣市回歸系數(shù)為負(fù)，西南部伏牛山區(qū)回歸系數(shù)同樣較低，其中表現(xiàn)為負(fù)影響的地區(qū)主要集中在商丘市、濮陽市等部分縣市以及西南部伏牛山區(qū)，主要由于這些縣市地區(qū)氣溫降雨、地形地貌等自然條件較差，農(nóng)民人均純收入較低，城鎮(zhèn)化協(xié)調(diào)發(fā)展受到負(fù)面影響。而與此相反的是在地理區(qū)位和氣候條件較好的地區(qū)農(nóng)民人均純收入對城鎮(zhèn)化的影響大體上表現(xiàn)為正。(3) 人均社會消費品零售總額對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響有正有負(fù)，從濮陽到開封市、鄭州市、許昌市、漯河市一線的影響系數(shù)為負(fù)，河南省東部、東南部影響系數(shù)稍微較高，河南省的西南部影響系數(shù)最高。(4) 建筑業(yè)總產(chǎn)值對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響系數(shù)從河南省南部向北部顯著增加，由負(fù)變?yōu)檎?。從三門峽市、洛陽市、鄭州市、開封市、商丘市以北的縣市，建筑業(yè)總產(chǎn)值對于城鎮(zhèn)化協(xié)調(diào)度的影響為正向，以南的縣市影響系數(shù)基本為負(fù)。(5) 平均海拔對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響有正有負(fù)，西南山區(qū)平均海拔較高抑制了城鎮(zhèn)化的協(xié)調(diào)發(fā)展，對于河南省東部地區(qū)平原地區(qū)，海拔高度在一定程度上促進(jìn)了城鎮(zhèn)化的協(xié)調(diào)發(fā)展。(6) 以京廣線為分界線，以東地區(qū)人均糧食產(chǎn)量對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度有正向影響；以西地區(qū)人均糧食產(chǎn)量對于城鎮(zhèn)化協(xié)調(diào)發(fā)展度有負(fù)向影響。河南省東部地區(qū)地勢平坦，是河南省糧食主產(chǎn)地區(qū)，有利地支撐了城鎮(zhèn)化的快速發(fā)展，西南部山區(qū)自然條件相對較差，人均糧食產(chǎn)量較低影響了城鎮(zhèn)化的健康發(fā)展。(7) 地方財政人均農(nóng)林水利支出對城鎮(zhèn)化協(xié)調(diào)發(fā)展度的影響基本為負(fù)向，在河南省東北部的少部分地區(qū)影響系數(shù)為正，影響系數(shù)最小的在河南的東南地區(qū)和西北地區(qū)。安陽市、濮陽市、開封市和商丘市的部分地區(qū)需要加大地方財政人均農(nóng)林水利支出以促進(jìn)城鎮(zhèn)化的協(xié)調(diào)發(fā)展。河南省東北地區(qū)濮陽市、新鄉(xiāng)市、焦作市的部分縣市處于黃河灘區(qū)，洪澇災(zāi)害的治理投入對于城鎮(zhèn)化協(xié)調(diào)發(fā)展有重要影響。

4 結(jié)論與討論

4.1 結(jié) 論

(1) 在研究期內(nèi)河南省城鎮(zhèn)化的整體水平、各子系統(tǒng)發(fā)展水平及其協(xié)調(diào)發(fā)展度都處于逐年提升的狀態(tài)；2000—2014年河南省城鎮(zhèn)化歷經(jīng)了經(jīng)濟(jì)人口城鎮(zhèn)化滯后到土地城鎮(zhèn)化滯后的轉(zhuǎn)變；地區(qū)間人口城鎮(zhèn)化發(fā)展指數(shù)、經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)、城鎮(zhèn)化綜合發(fā)展度與城鎮(zhèn)化協(xié)調(diào)發(fā)展度總體差異呈下降趨勢，而土地城鎮(zhèn)化發(fā)展指數(shù)在區(qū)域之間差異有逐漸變大的趨勢。

(2) 從全局空間自相關(guān)來看，河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展整體有較強(qiáng)的空間正相關(guān)；從局部空間自相關(guān)來看，高—高區(qū)主要分布在鄭州市及其周邊地區(qū)，而低—低區(qū)主要分布在河南省南部地區(qū)。

(3) 河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展的9個影響因素的影響力由大到小依次為人口分布密度、人均教育經(jīng)費支出、人均地方財政一般預(yù)算收入、平均海拔、農(nóng)民人均純收入、人均社會消費品零售總額、人均糧食產(chǎn)量、建筑業(yè)總產(chǎn)值、地方財政人均農(nóng)林水利支出，建筑業(yè)總產(chǎn)值與地方財政人均農(nóng)林水利支出整體上表現(xiàn)為負(fù)相關(guān)，其余指標(biāo)表現(xiàn)為正相關(guān)。

4.2 討 論

(1) 2000—2014年河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展快速推進(jìn)，但是到2014年仍然有很多縣市處于低水平的協(xié)調(diào)狀態(tài)，也有部分地區(qū)由于城鎮(zhèn)化子系統(tǒng)發(fā)展水平差異較大，導(dǎo)致協(xié)調(diào)發(fā)展度較低。河南省政府在推動城鎮(zhèn)化協(xié)調(diào)發(fā)展的過程中應(yīng)結(jié)合區(qū)域?qū)嶋H情況，對于人口城鎮(zhèn)化發(fā)展水平較低的縣市應(yīng)該不斷完善城鎮(zhèn)體系、產(chǎn)業(yè)發(fā)展、交通網(wǎng)絡(luò)、基礎(chǔ)設(shè)施和公共服務(wù)布局，保證公共服務(wù)均衡發(fā)展，同時注重保障和改善民生，堅持民生優(yōu)先、共建共享，系統(tǒng)推進(jìn)民生幸福工程，促進(jìn)社會公平正義，保證就業(yè)機(jī)會充分，增強(qiáng)廣大市民的幸福感；對于經(jīng)濟(jì)城鎮(zhèn)化發(fā)展指數(shù)較低的地區(qū)應(yīng)該加強(qiáng)投資力度，推動城鄉(xiāng)生產(chǎn)要素優(yōu)化組合、促進(jìn)城鄉(xiāng)協(xié)調(diào)發(fā)展，引進(jìn)新科學(xué)、新技術(shù)、新方法，推動經(jīng)濟(jì)結(jié)構(gòu)調(diào)整和產(chǎn)業(yè)結(jié)構(gòu)不斷升級，加快產(chǎn)業(yè)轉(zhuǎn)型步伐，推動結(jié)構(gòu)優(yōu)化、集群發(fā)展、技術(shù)進(jìn)步，提高經(jīng)濟(jì)效益；對于土地城鎮(zhèn)化發(fā)展水平較低的地區(qū)，需要統(tǒng)籌考慮各區(qū)域的資源、生態(tài)、環(huán)境和經(jīng)濟(jì)社會發(fā)展需求，根據(jù)生態(tài)健康、生活宜居、生產(chǎn)便利的要求，促進(jìn)布局優(yōu)化，統(tǒng)籌協(xié)調(diào)布局城鎮(zhèn)村生活空間、各產(chǎn)業(yè)生產(chǎn)空間和生態(tài)環(huán)境保護(hù)空間，形成功能清晰、協(xié)調(diào)有序的國土空間開發(fā)利用格局。在實現(xiàn)經(jīng)濟(jì)質(zhì)量增長和效益同步提高的同時，注重資源節(jié)約和環(huán)境保護(hù)，堅持速度和效益相結(jié)合，走高技術(shù)、低污染、全面、協(xié)調(diào)、可持續(xù)發(fā)展的新型工業(yè)化道路，達(dá)到更高級、更有組織、更高效率的經(jīng)濟(jì)發(fā)展水平。河南省“三山一灘”地區(qū)多為貧困縣，經(jīng)濟(jì)發(fā)展水平較低，自然條件較為惡劣，成為河南省發(fā)展的洼地和連片貧困區(qū)，這些地區(qū)城鎮(zhèn)化發(fā)展質(zhì)量較低。對于黃河灘區(qū)，應(yīng)該結(jié)合實際情況，合理轉(zhuǎn)移居民，加大農(nóng)林水利支出，解決防洪安全；對于“三山”地區(qū)應(yīng)該加快基礎(chǔ)設(shè)施建設(shè)，改造提升、轉(zhuǎn)移就業(yè)、實現(xiàn)產(chǎn)業(yè)發(fā)展和生態(tài)保護(hù)相結(jié)合，發(fā)展特色旅游，探索扶貧新路徑。

(2) 本文以河南省108個縣市為研究對象，分析了河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展的時空格局變化、區(qū)域差異、空間相關(guān)性和影響因素。在指標(biāo)選取上受到數(shù)據(jù)獲取的原因略顯簡單，在影響因子的選取上較多地考慮了社會經(jīng)濟(jì)因素，對自然因素的考慮有所欠缺，在以后的研究中有待進(jìn)一步深化。

參考文獻(xiàn)：

[1] 國務(wù)院發(fā)展研究中心和世界銀行聯(lián)合課題組,李偉, Sri Mulyani,等.中國推進(jìn)高效、包容、可持續(xù)的城鎮(zhèn)化[J].管理世界,2014(4):5-41.

[2] 樊杰,劉毅,陳田,等.優(yōu)化我國城鎮(zhèn)化空間布局的戰(zhàn)略重點與創(chuàng)新思路[J].中國科學(xué)院院刊,2013,28(1):20-27.

[3] 冷智花,付暢儉.城鎮(zhèn)化失衡發(fā)展對糧食安全的影響[J].經(jīng)濟(jì)學(xué)家,2014(11):58-65.

[4] 呂添貴,吳次芳,李洪義,等.人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)性測度及優(yōu)化:以南昌市為例[J].地理科學(xué),2016,36(2):239-246.

[5] Lin G C S. Chinese urbanism in question: State, society, and the reproduction of urban spaces[J]. Urban Geography, 2007,28(1):7-29.

[6] 李鑫,李興校,歐名豪.江蘇省城鎮(zhèn)化發(fā)展協(xié)調(diào)度評價與地區(qū)差異分析[J].人文地理,2012(3):50-54.

[7] 曹文莉,張小林,潘義勇,等.發(fā)達(dá)地區(qū)人口、土地與經(jīng)濟(jì)城鎮(zhèn)化協(xié)調(diào)發(fā)展度研究[J].中國人口·資源與環(huán)境,2012,22(2):141-146.

[8] 姚士謀,張平宇,余成,等.中國新型城鎮(zhèn)化理論與實踐問題[J].地理科學(xué),2014,34(6):641-647.

[9] 李杰,崔許鋒.山地民族地區(qū)城鎮(zhèn)化的非均衡性與時空分異研究[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2016(4):84-93.

[10] 國務(wù)院.國家新型城鎮(zhèn)化規(guī)劃(2014—2020年)[EB/OL]. http:∥www. gov. cn/xinwen/2014-03/16/content_ 2639841. html.2014.

[11] 美克拉衣,艾克拜爾.中國城鎮(zhèn)化推進(jìn)模式研究[J].城市地理,2016(6):97.

[12] 王新越,秦素貞,吳寧寧.新型城鎮(zhèn)化的內(nèi)涵、測度及其區(qū)域差異研究[J].地域研究與開發(fā),2014,33(4):69-75.

[13] Friedmann J. Four theses in the study of China′s urbanization[J]. International Journal of Urban and Regional Research, 2006,30(2):440-451.

[14] 陳明星,陸大道,劉慧.中國城市化與經(jīng)濟(jì)發(fā)展水平關(guān)系的省際格局[J].地理學(xué)報,2010,65(12):1443-1453.

[15] 凌筱舒,王立,薛德升.江西省縣域城鎮(zhèn)化水平測度及其分異研究[J].人文地理,2014(3):89-94.

[16] 王洋,方創(chuàng)琳,王振波.中國縣域城鎮(zhèn)化水平的綜合評價及類型區(qū)劃分[J].地理研究,2012,31(7):1305-1316.

[17] 張榮天,焦華富.長江三角洲地區(qū)城鎮(zhèn)化效率測度及空間關(guān)聯(lián)格局分析[J].地理科學(xué),2015,35(4):433-439.

[18] 丁志偉,張改素,王發(fā)曾.中原經(jīng)濟(jì)區(qū)“三化”協(xié)調(diào)的內(nèi)在機(jī)理與定量分析[J].地理科學(xué),2013,33(4):402-409.

[19] 王瑜煒,秦輝.中國信息化與新型工業(yè)化耦合格局及其變化機(jī)制分析[J].經(jīng)濟(jì)地理,2014,34(2):93-100.

[20] 肖海平,肖地楚,朱佩芬.湖南郴州市新型城鎮(zhèn)化與新型工業(yè)化協(xié)調(diào)發(fā)展研究[J].地域研究與開發(fā),2015,34(3):39-44.

[21] 祖佳嬉,葉長盛,雷慧敏.江西省縣域“三化”協(xié)調(diào)發(fā)展度空間格局演變及其障礙因子診斷[J].地域研究與開發(fā),2015,34(2):34-39.

[22] 徐維祥,舒季君,唐根年.中國工業(yè)化、信息化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化協(xié)調(diào)發(fā)展的時空格局與動態(tài)演進(jìn)[J].經(jīng)濟(jì)學(xué)動態(tài),2015(1):76-85.

[23] 李小建,羅慶.新型城鎮(zhèn)化中的協(xié)調(diào)思想分析[J].中國人口·資源與環(huán)境,2014,24(2):47-53.

[24] 郭付友,李誠固,陳才,等.2003年以來東北地區(qū)人口城鎮(zhèn)化與土地城鎮(zhèn)化時空耦合特征[J].經(jīng)濟(jì)地理,2015,35(9):49-56.

[25] 王亞力,彭保發(fā),熊建新,等.2001年以來環(huán)洞庭湖區(qū)經(jīng)濟(jì)城鎮(zhèn)化與人口城鎮(zhèn)化進(jìn)程的對比研究[J].地理科學(xué),2014,34(1):67-75.

[26] 劉法威,許恒周,王姝.人口—土地—經(jīng)濟(jì)城鎮(zhèn)化的時空耦合協(xié)調(diào)性分析：基于中國省際面板數(shù)據(jù)的實證研究[J].城市發(fā)展研究,2014,21(8):7-11.

[27] 雒海潮,李國梁.河南省城鎮(zhèn)化協(xié)調(diào)發(fā)展評價與空間差異分析[J].地理科學(xué),2015,35(6):749-755.

[28] 胡燕燕,曹衛(wèi)東.近三十年來我國城鎮(zhèn)化協(xié)調(diào)性演化研究[J].城市規(guī)劃,2016,40(2):9-17.

[29] 王富喜,毛愛華,李赫龍,等.基于熵值法的山東省城鎮(zhèn)化質(zhì)量測度及空間差異分析[J].地理科學(xué),2013,33(11):1323-1329.

[30] 鹿晨昱,張琳,薛冰,等.遼寧省區(qū)域可持續(xù)發(fā)展時空綜合測度研究[J].經(jīng)濟(jì)地理,2015,35(8):32-39.

[31] 徐維祥,陳國亮,舒季君,等.基于空間連續(xù)性的“四化同步區(qū)”形成與演化機(jī)理研究[J].中國工業(yè)經(jīng)濟(jì),2015(5):18-31.

[32] Mikelbank B A. Quantitative geography: Perspectives on spatial data analysis, by A. S. Fotheringham, C. Brunsdon, and M. Charlton[J]. Geographical Analysis, 2001,33(4):370-372.

[33] 徐建華.地理建模方法[M].北京:科學(xué)出版社,2010.