張 慧，蘇 航，薛 睿 ,2 ，章桂芳

（1. 東北農(nóng)業(yè)大學(xué)公共管理與法學(xué)院，哈爾濱 150030；2. 哈爾濱工業(yè)大學(xué)建筑學(xué)院，哈爾濱 150001；3.中山大學(xué)地球科學(xué)與工程學(xué)院，廣州 510275）

0 引 言

為推進(jìn)國家治理體系和治理能力的現(xiàn)代化，原國土部、國家發(fā)改委、原環(huán)保部和住建部四部委2014 年決定在若干地區(qū)開展“多規(guī)合一”試點[1]，2018 年十九屆三中全會開啟自然資源管理體制改革序幕，正式將城鄉(xiāng)規(guī)劃管理職能劃歸自然資源部，國土空間開發(fā)利用的綜合性和系統(tǒng)性正不斷加強(qiáng)。2019 年5 月23 日，中共中央、國務(wù)院正式發(fā)布《關(guān)于建立國土空間規(guī)劃體系并監(jiān)督實施的若干意見》，強(qiáng)調(diào)要“整體謀劃新時代國土空間開發(fā)保護(hù)格局，綜合考慮人口分布、經(jīng)濟(jì)布局、國土利用、生態(tài)環(huán)境保護(hù)等因素”，并指出“省級國土空間規(guī)劃是對全國國土空間規(guī)劃的落實......側(cè)重協(xié)調(diào)性”[2]。本文從土地、人口與經(jīng)濟(jì)多維協(xié)調(diào)的視角去探究省域國土空間開發(fā)策略，是破解當(dāng)前城鄉(xiāng)差距和人地矛盾問題的關(guān)鍵。只有土地、人口與經(jīng)濟(jì)都保持良好發(fā)展并且相互協(xié)調(diào)才是城鄉(xiāng)協(xié)調(diào)發(fā)展的應(yīng)有之義，新型城鎮(zhèn)化的訴求才能夠被滿足，鄉(xiāng)村振興戰(zhàn)略才得以實現(xiàn)。

早在戰(zhàn)國時期的《商君書》提出“制土分民之律”，中國就已經(jīng)認(rèn)識到了人口和土地二者關(guān)系的調(diào)控在國家管理中的重要性，認(rèn)為國家發(fā)展一直處于“民勝其地”和“地勝其民”的反復(fù)之中[3]。改革開放以后，國內(nèi)學(xué)者開始用現(xiàn)代的觀點和方法研究人口、土地、經(jīng)濟(jì)、資源、環(huán)境、社會等多要素的交互影響[4-9]。張志強(qiáng)將可持續(xù)發(fā)展問題歸結(jié)為“人口、資源、環(huán)境和發(fā)展”四種問題的復(fù)合體[10]，吳攀升指出在研究過程中要揭示人口和土地對經(jīng)濟(jì)發(fā)展的影響力和承載力，以達(dá)到一種人地關(guān)系的平衡狀態(tài)[11]。也有學(xué)者從自然、經(jīng)濟(jì)、社會和生態(tài)環(huán)境等多方視角研究土地的綜合承載能力[12-16]。還有部分學(xué)者已經(jīng)對人口、土地、生態(tài)、產(chǎn)業(yè)、科技、經(jīng)濟(jì)等要素耦合協(xié)調(diào)關(guān)系，及其對城鎮(zhèn)化的影響進(jìn)行了研究[17-21]。國外有關(guān)土地與其他系統(tǒng)協(xié)調(diào)性的研究中，Bartholomew 等對土地利用和交通方式轉(zhuǎn)變間的聯(lián)系展開了探索[22]，Deng 等 梳理了土地利用與生態(tài)系統(tǒng)關(guān)系的研究框架[23]，Paracchini 等對土地多功能性間的協(xié)調(diào)評估進(jìn)行了探討[24]。Gerundo 等分別研究了不同政策與土地利用間的協(xié)同作用[25-28]?，F(xiàn)有研究大多采用指標(biāo)組合賦權(quán)后進(jìn)行多因子綜合評價的方式進(jìn)行協(xié)調(diào)性研究，評價方法較單一且研究層次集中在空間特征、格局的探索等方面[29-31]，依據(jù)協(xié)調(diào)性的多維分異來指導(dǎo)空間規(guī)劃實踐的研究較為罕見。本文首先針對性地采用獨立模型，即土地經(jīng)營規(guī)模模型、logistic 模型、結(jié)節(jié)性指數(shù)對黑龍江各地市土地壓力、人口缺口、經(jīng)濟(jì)水平進(jìn)行了分別測度，并從總體和局部上聯(lián)合測度了三者的多維協(xié)調(diào)性。進(jìn)而，依據(jù)對協(xié)調(diào)性分異的分析來揭示各地發(fā)展模式，基于此進(jìn)行空間分區(qū)并提出分區(qū)發(fā)展策略，可以為國土空間綜合開發(fā)以及相關(guān)研究提供有利借鑒。

黑龍江省土地優(yōu)質(zhì)肥沃，適宜于大規(guī)模集約化農(nóng)業(yè)生產(chǎn)，但農(nóng)村土地零星破碎成為了開展規(guī)模經(jīng)營的桎梏。同時黑龍江省人口自然增長率低下且人口省外遷出嚴(yán)重，導(dǎo)致城鎮(zhèn)集聚效應(yīng)難以發(fā)揮，城鎮(zhèn)化動力不足。并且經(jīng)濟(jì)上嚴(yán)重依賴農(nóng)業(yè)和重工業(yè)，導(dǎo)致近年來黑龍江省經(jīng)濟(jì)增長乏力，在全國水平上排名處于末流。土地、人口、經(jīng)濟(jì)上發(fā)展不協(xié)調(diào)的矛盾已經(jīng)開始顯現(xiàn)，故本文在上述背景下選取黑龍江省進(jìn)行研究。

1 材料與方法

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

受2015 年開始的農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革的影響，研究區(qū)糧食價格產(chǎn)生了巨大的波動，且這種影響至今尚未消弭，為保證數(shù)據(jù)的穩(wěn)健性，故本文研究期限截止于2015年。本文城鄉(xiāng)社會經(jīng)濟(jì)原始數(shù)據(jù)來源于2011—2016 年《黑龍江省統(tǒng)計年鑒》、《全國農(nóng)產(chǎn)品成本收益資料匯編》和《中國住戶調(diào)查年鑒》等。本文人口數(shù)據(jù)來源于1976—2016 年《黑龍江統(tǒng)計年鑒》。本文土地利用數(shù)據(jù)來源于2016 年土地利用數(shù)據(jù)庫，提取黑龍江省各地市耕地面積并匯總得到。

1.2 研究方法

1.2.1 土地規(guī)模經(jīng)營模型

本文采用的是城鄉(xiāng)收入均等化作為農(nóng)戶理想土地經(jīng)營規(guī)模的標(biāo)尺，即農(nóng)民單純依靠務(wù)農(nóng)收入能與城鎮(zhèn)職工家庭收入持平。辛嶺等沿此思路采用農(nóng)民工打工年收入與每公頃糧食凈利潤之比來表征最優(yōu)土地經(jīng)營規(guī)模[32]（模型1）。但相同生活水平下農(nóng)村的生活成本要低于城市，為了抵消城市與農(nóng)村二者的絕對收入差異給人造成的誤導(dǎo)，因此同時參考岳立教授反映城鄉(xiāng)之間相對生活成本差異的“支出—收入”之比的城鄉(xiāng)生活成本差異模型[33]（模型2），采用（2）式中生活成本差異指數(shù)Ix 對（1）式進(jìn)行校正?？紤]數(shù)據(jù)可得性導(dǎo)致的模型擬合精度影響，本文僅考慮靜態(tài)工資水平?？紤]國情變化與地方實際，本文將模型1 中IC 進(jìn)行了改進(jìn)，改進(jìn)后的適度經(jīng)營規(guī)模模型如公式（3）

式中S 代表土地經(jīng)營規(guī)模，IC 代表打工收入（元），RP為公頃土地平均凈利潤（元）。Ix 是城鄉(xiāng)生活成本差異指數(shù)；Cv表示農(nóng)村居民支出（元），Iv表示農(nóng)村居民收入（元）；Cu為城市居民的支出（元），Iu為農(nóng)民落戶城市后的收入（元）。S’為農(nóng)戶土地經(jīng)營面積（公頃）。本文將（1）式中家庭收入IC 核算方式進(jìn)行改進(jìn)，將IC 換算為Inc×La，具體含義如下：黨中央出臺了放開落戶限制、引導(dǎo)農(nóng)民工市民化等多項政策，故本文認(rèn)為農(nóng)民向城鎮(zhèn)的轉(zhuǎn)移是長期性舉家轉(zhuǎn)移，將受供養(yǎng)人口剔除后的戶均勞動力數(shù)量視同為進(jìn)城務(wù)工人數(shù)，即式中La。由于農(nóng)村勞動力向城鎮(zhèn)轉(zhuǎn)移之后主要在私營企業(yè)從事勞動，故本文認(rèn)為工資口徑應(yīng)采用的是城鎮(zhèn)私營企業(yè)人均工資，即Inc。RP'為研究區(qū)3 種主要糧食作物（玉米、水稻、大豆）每公頃平均凈利潤。

1.2.2 飽和城鎮(zhèn)化率的計算模型

本文采用logistic 模型測度飽和城鎮(zhèn)化率，該方法被聯(lián)合國廣泛用于世界各地區(qū)城鎮(zhèn)化水平預(yù)測，故又稱“聯(lián)合國法”。即假設(shè)在預(yù)測期內(nèi)城鄉(xiāng)人口增長率差保持不變，根據(jù)已知城鄉(xiāng)人口調(diào)查的數(shù)據(jù)，外推預(yù)測期末的城鎮(zhèn)人口占總?cè)丝诒壤Ｆ涫綖椋?/p>

式中PU 為城市人口比例，t 為時間，k 為相對增長率，α為系數(shù)。通常利用式（4）預(yù)測某一年度的城鎮(zhèn)化水平，但其默認(rèn)城市化率的飽和值為100%，這顯然不符合實際，因此本文參考陳彥光教授的相關(guān)研究[34-35]，將式（4）右側(cè)分子部分飽和城鎮(zhèn)化水平因子1 換為參數(shù)β。

式中β 代表城鎮(zhèn)化率的飽和值，并將（5）式變換后取對數(shù)轉(zhuǎn)換為線性形式，利用回歸分析對β 值進(jìn)行擬合。其公式如下：

1.2.3 地區(qū)經(jīng)濟(jì)水平評價方法

由于對地區(qū)經(jīng)濟(jì)水平的單純評價不涉及多系統(tǒng)間復(fù)雜權(quán)重分配，因此本文指標(biāo)選取過程，充分考慮地區(qū)客觀實際、數(shù)據(jù)可得性和專家學(xué)者相關(guān)研究[9,13,17,21]，構(gòu)建經(jīng)濟(jì)水平評價指標(biāo)體系。各指標(biāo)權(quán)重利用層次分析軟件YAAHP賦得，CI政府投資能力=0.028＜0.1，CI農(nóng)戶投資能力=0.004＜0.1，結(jié)果均通過一致性檢驗，確定最終權(quán)重見表1。

對各地區(qū)經(jīng)濟(jì)發(fā)展水平的評價，采用區(qū)域發(fā)展水平測度通常采用的結(jié)節(jié)性指數(shù)，函數(shù)形式如下：

式中Eco 為經(jīng)濟(jì)發(fā)展水平結(jié)節(jié)性指數(shù)，Kj為第j 項指標(biāo)標(biāo)準(zhǔn)化后的分值；Kj為第j 項指標(biāo)因子權(quán)重；n 為指標(biāo)數(shù)目，此處n=6。

1.2.4 秩相關(guān)系數(shù)

秩相關(guān)系數(shù)是描述兩要素之間相關(guān)程度的統(tǒng)計指標(biāo)，它將兩要素樣本值位次替代實際數(shù)據(jù)而求得的統(tǒng)計量，是位次分析的數(shù)量化，本文用以測度系統(tǒng)間的總體絕對相關(guān)性。

式中R 為秩相關(guān)系數(shù)，di為要素間位次差，n 為樣本量。

表1 經(jīng)濟(jì)評價指標(biāo)體系 Table 1 Economic evaluation index system

1.2.5 核密度分析法

核密度估計 (kernel density estimation, KDE) 是基于數(shù)據(jù)密集度函數(shù)聚類算法的一種空間密度分析方法, 分析過程中賦予樣方中心點鄰近事件較高權(quán)重, 而距離中心點較遠(yuǎn)的事件被賦予較低權(quán)重, 其四次多項式核函數(shù)為

式中s 為待估計點的位置；sj為落在以s 為圓心，h 為半徑的圓形范圍內(nèi)的第i 個研究單元的位置，h 代表步長，即以p 為源點的曲面在空間上延展的寬帶。

1.2.6 耦合協(xié)調(diào)度

耦合度協(xié)調(diào)度借鑒來源于物理學(xué)中容量耦合系數(shù)模型，是可以用以定量表征多要素間協(xié)調(diào)程度的指標(biāo)，本文用以探究系統(tǒng)間各要素的局部協(xié)調(diào)性，式為

式中D 表示耦合協(xié)調(diào)度，C 為耦合度，U1表示過剩農(nóng)村人口，U2表示城鎮(zhèn)化人口缺口，U3表示地區(qū)經(jīng)濟(jì)水平。φ、ρ 和β 為權(quán)重參數(shù)。本文認(rèn)為土地、人口、經(jīng)濟(jì)方面因素同等重要故φ、ρ、β 取均等值0.333。在耦合協(xié)調(diào)度的計算中，為使多系統(tǒng)間指標(biāo)可比，對數(shù)據(jù)進(jìn)行極值標(biāo)準(zhǔn)化以去除量綱。

2 人口、土地、經(jīng)濟(jì)系統(tǒng)評價

模型（3）依據(jù)城鄉(xiāng)收入均等化的思路，得出黑龍江省農(nóng)村耕地適度經(jīng)營規(guī)模在戶均11.6 hm2左右，則根據(jù)耕地總量和戶均人數(shù)計算可知黑龍江省農(nóng)村人口承載規(guī)模為469 萬人，該情況下過剩農(nóng)業(yè)人口為1 102 萬人，農(nóng)村土地總壓力為51 478 940 hm2。將上述方法應(yīng)用于各地市，可分別計算出各地市農(nóng)村土地壓力；利用模型（6），當(dāng)擬合參數(shù)β 值到0.795 時回歸方程y=-0.154 9x+0.216 6 的R2達(dá)到0.950 1，可視為在5%的水平上顯著，擬合結(jié)果是可信的。得到黑龍江省總體城鎮(zhèn)化率飽和值為79.5%，以研究期末3812 萬人口為基準(zhǔn)，可得城鎮(zhèn)人口承載力為3 031 萬，該情況下城鎮(zhèn)人口缺口數(shù)值為789 萬人，同理，根據(jù)各地市現(xiàn)狀城鎮(zhèn)人口數(shù)據(jù)計算出各地市與飽和城鎮(zhèn)化率間的人口缺口；利用模型（7），以及表1 中相應(yīng)的指標(biāo)權(quán)重計算出地區(qū)經(jīng)濟(jì)水平分值。具體結(jié)果見表2。

表2 單系統(tǒng)評價結(jié)果 Table 2 Single system evaluation results

表2 顯示黑龍江省13 個地市的農(nóng)村土地壓力全部為正，各地均存在農(nóng)村人口過剩的現(xiàn)象?？臻g上表現(xiàn)出了由西向東和由南向北逐漸減弱的特征，核密度中心出現(xiàn)在哈爾濱（圖1a），西部和南部各地市農(nóng)村人口壓力明顯高于東部和北部地市，形成了包括哈爾濱-綏化-齊齊哈爾-牡丹江-大慶的高值區(qū)和包括雞西-佳木斯-雙鴨山-七臺河-伊春-黑河-鶴崗-大興安嶺的相對低值區(qū)（數(shù)值從高到低順序）。

城鎮(zhèn)人口缺口表現(xiàn)出正負(fù)共存的特征，由核密度分布圖（圖1b），在空間上同樣表現(xiàn)出了由西向東和由南向北逐漸減弱的特征，核密度中心出現(xiàn)在哈爾濱，西部和南部各地市城鎮(zhèn)化缺口明顯高于東部和北部地市，形成了包括哈爾濱-綏化-齊齊哈爾-牡丹江-大慶的高值區(qū)和包括佳木斯-雙鴨山-黑河-雞西-七臺河-伊春-鶴崗-大興安嶺的相對低值區(qū)。

經(jīng)濟(jì)水平（圖1c）空間集聚特征并不明顯。但與土地壓力和人口缺口相同，3 個最高值大慶-哈爾濱-牡丹江仍然在黑龍江省的西南部聚集，核密度中心出現(xiàn)在大慶。不同的是，人口壓力和土地壓力空間分布中的兩個傳統(tǒng)高值地區(qū)齊齊哈爾和綏化，反而成為了經(jīng)濟(jì)水平分值最低的地市。而黑龍江省北部的大興安嶺地區(qū)和黑河市是人口缺口和土地壓力的低值區(qū)，在經(jīng)濟(jì)水平中卻成為了得分僅次于大慶、哈爾濱、牡丹江、雞西、佳木斯的次高值區(qū)。

圖1 土地、人口、經(jīng)濟(jì)系統(tǒng)核密度圖 Fig.1 KDE of land, population and economic system

3 土地-人口-經(jīng)濟(jì)多維協(xié)調(diào)性分析

秩（復(fù)）相關(guān)系數(shù)的結(jié)果表征的是系統(tǒng)間的絕對關(guān)系，但無法定位系統(tǒng)內(nèi)部各要素的具體協(xié)調(diào)程度。耦合協(xié)調(diào)度能定位系統(tǒng)內(nèi)部元素的協(xié)調(diào)程度，但表征的是相對關(guān)系。因此本文首先采用秩(復(fù))相關(guān)系數(shù)進(jìn)行總體協(xié)調(diào)水平定性，進(jìn)而采用耦合協(xié)調(diào)度對內(nèi)部要素間協(xié)調(diào)水平進(jìn)行劃分，這種處理方式保證了結(jié)果的科學(xué)準(zhǔn)確。

宗教對文學(xué)的影響，屬于比較文學(xué)的跨學(xué)科研究的范疇。在文學(xué)產(chǎn)生和發(fā)展的過程中，宗教曾經(jīng)起過重要作用；宗教會對作家的世界觀發(fā)生影響，進(jìn)而影響到文學(xué)作品的思想內(nèi)容，有些重要的宗教經(jīng)典，如《圣經(jīng)》，本身就是文學(xué)作品，或者取自文學(xué)作品。[1]274-275

3.1 總體協(xié)調(diào)性驗證

利用秩相關(guān)系數(shù)驗證人口-土地兩系統(tǒng)的相關(guān)程度R=0.945＞R0.01=0.703，二者在超過0.01 的顯著水平上相關(guān)，可見黑龍江省農(nóng)村土地壓力與城鎮(zhèn)人口缺口在總體上已表現(xiàn)出了很好的協(xié)調(diào)性（圖3a），因此將人口-土地協(xié)調(diào)關(guān)系采用自然斷點法進(jìn)一步劃分為3 個等級：D 取值在0～0.323 之間為較差協(xié)調(diào)，D 取值在0.323～0.562 之間為較好協(xié)調(diào)，在0.562～1 之間為很好協(xié)調(diào)。

經(jīng)濟(jì)-土地系統(tǒng)的秩相關(guān)系數(shù)R 為0.137，協(xié)調(diào)性未通過秩相關(guān)驗證。但對比土地-經(jīng)濟(jì)協(xié)調(diào)關(guān)系，兩系統(tǒng)間大部分要素的位次變化量相似。將離群值綏化市、齊齊哈爾市、大興安嶺地區(qū)剔除后另行計算秩相關(guān)系數(shù)。得到R校正=0.618＞R0.05=0.564，二者在超過0.05 的顯著水平上相關(guān)，因此判定秩相關(guān)系數(shù)失真的現(xiàn)象是由個別離群值引起的，剔除離群值后土地與經(jīng)濟(jì)系統(tǒng)間達(dá)到很好的總體協(xié)調(diào)水平。即土地壓力越小的地區(qū)經(jīng)濟(jì)水平越高，土地的規(guī)模經(jīng)營和地區(qū)經(jīng)濟(jì)水平存在較明顯的正相關(guān)。因此本文將各地市土地-經(jīng)濟(jì)系統(tǒng)協(xié)調(diào)關(guān)系進(jìn)一步劃分為3 個等級：D 取值在0～0.460 之間的三個離群值劃定為不協(xié)調(diào)，其余采用自然斷點法，D 取值在0.460～0.693 之間的為較好協(xié)調(diào)，在0.693～1 之間的為很好協(xié)調(diào)。

經(jīng)濟(jì)-人口兩系統(tǒng)的秩相關(guān)系數(shù)R 為0.126。對比人口-經(jīng)濟(jì)位次關(guān)系將離群值綏化市、齊齊哈爾市、大興安嶺地區(qū)剔除后另行測算秩相關(guān)系數(shù)。得到 R校正=0.594＞ R0.05=0.564，其顯著性超過95%。同樣可以判定，受離群值的影響，秩相關(guān)系數(shù)對人口-經(jīng)濟(jì)實際協(xié)調(diào)關(guān)系亦存在失真的描述。剔除離群值后相關(guān)關(guān)系顯著，這表明黑龍江省人口缺口越小的地區(qū)經(jīng)濟(jì)水平越高，城鎮(zhèn)化發(fā)展水平和地區(qū)經(jīng)濟(jì)強(qiáng)度間存在比較顯著的正向關(guān)聯(lián)。因此將各地市人口-經(jīng)濟(jì)協(xié)調(diào)關(guān)系進(jìn)一步劃分為3 個等級：D 取值在0～0.405 之間的3 個離群值劃定為不協(xié)調(diào)，其余采用自然斷點法，D 取值在0.405～0.684之間的為較好協(xié)調(diào)，D 取值在0.684～1 之間的為很好協(xié)調(diào)。

在衡量人口-土地-經(jīng)濟(jì)三系統(tǒng)間協(xié)調(diào)性時，利用SPSS 計算三系統(tǒng)間復(fù)相關(guān)系數(shù)，得到R'=2.59。對復(fù)相關(guān)系數(shù)需進(jìn)行F 檢驗，經(jīng)計算R'=2.59＞F0.1=2.56，可見三系統(tǒng)在總體上具有很好的協(xié)調(diào)性。在此基礎(chǔ)上將耦合協(xié)調(diào)度采用自然斷點法劃分為三類：D 取值在0.00～0.141 之間為較差協(xié)調(diào)，D 取值在0.141～0.478 之間為較好協(xié)調(diào)，在0.478～1.00 之間為很好協(xié)調(diào)。

3.2 局部協(xié)調(diào)性分析

通過總體協(xié)調(diào)性驗證，可知土地、人口、經(jīng)濟(jì)三系統(tǒng)間均存在顯著的相關(guān)，并且劃分出了局部協(xié)調(diào)性的分級標(biāo)準(zhǔn)，運用耦合協(xié)調(diào)度模型對三系統(tǒng)間的局部協(xié)調(diào)性進(jìn)行測度，結(jié)果如表3。為探索局部協(xié)調(diào)性的空間特征，將表中數(shù)據(jù)進(jìn)行空間可視化表達(dá)，如圖2 所示。

表3 局部協(xié)調(diào)性評價結(jié)果 Table 3 Micro-coordination evaluation results

各地市土地壓力與經(jīng)濟(jì)水平耦合協(xié)調(diào)關(guān)系為很好協(xié)調(diào)的有7 個占54%；較好協(xié)調(diào)的有3 個占23%；較差協(xié)調(diào)的有3 個占23%。在空間分布上（圖2b），體現(xiàn)出一定的空間分異規(guī)律，矛盾仍然主要集中在黑龍江省西南部地帶，同樣說明了黑龍江省西部和南部個別地市是土地利用和經(jīng)濟(jì)協(xié)調(diào)發(fā)展矛盾比較尖銳的地區(qū)。綜上，可以說明黑龍江省的農(nóng)村土地壓力和地區(qū)經(jīng)濟(jì)發(fā)展水平形成了很好的耦合協(xié)調(diào)關(guān)系，土地的規(guī)模經(jīng)營程度和地區(qū)經(jīng)濟(jì)發(fā)展水平之間存在比較明顯的正相關(guān)?？梢?，通過轉(zhuǎn)移農(nóng)業(yè)過剩人口來擴(kuò)大經(jīng)營主體的種植規(guī)模和提高農(nóng)業(yè)生產(chǎn)效率是相當(dāng)必要的。

各地市經(jīng)濟(jì)水平和城鎮(zhèn)人口缺口耦合協(xié)調(diào)關(guān)系為很好協(xié)調(diào)的有7 個占54%；較好協(xié)調(diào)的有3 個占23%；不協(xié)調(diào)的的有3 個占23%?？臻g分異特征明顯（圖2c），矛盾集中點仍然處于黑龍江省西南部區(qū)域，說明該地區(qū)個別地市是城鎮(zhèn)化進(jìn)程和經(jīng)濟(jì)協(xié)調(diào)發(fā)展矛盾比較尖銳的地區(qū)。各地的城鎮(zhèn)化進(jìn)程和地區(qū)經(jīng)濟(jì)發(fā)展水平形成了很好的耦合協(xié)調(diào)關(guān)系，農(nóng)村人口的市民化程度和地區(qū)經(jīng)濟(jì)發(fā)展水平之間存在顯著的正相關(guān)。黑龍江作為人口密度稀疏的邊疆地區(qū)，城鎮(zhèn)發(fā)育水平普遍不足，只有促進(jìn)農(nóng)業(yè)過剩人口城鎮(zhèn)化，才能發(fā)揮城鎮(zhèn)的增長極核作用，從而帶動地區(qū)發(fā)展。

黑龍江省土地-人口-經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)關(guān)系為很好協(xié)調(diào)的有10 個占77%；較好協(xié)調(diào)的有1 個占8%；較差協(xié)調(diào)的有2 個占15%?？臻g上，體現(xiàn)出了明顯的分異特征，耦合協(xié)調(diào)度較差的地區(qū)仍然聚集在黑龍江省西南部地區(qū)，協(xié)調(diào)性較好的地區(qū)集中連片分布在三江平原和北部林區(qū)，最高值為大慶市。通過圖2d 可知，同樣作為協(xié)調(diào)性較差的三個地區(qū)，哈爾濱市的土地、人口、經(jīng)濟(jì)三系統(tǒng)間變化趨勢是趨同的，而齊齊哈爾和綏化的土地、經(jīng)濟(jì)、人口系統(tǒng)間變化趨勢是拮抗的。三者之間仍存在內(nèi)部差異，他們的矛盾動因以及應(yīng)采取的管控策略是不同的，下文將進(jìn)一步分析。

圖2 人口、土地、經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)格局 Fig.2 Types of population, land and economy coupling coordination relations

4 協(xié)調(diào)性分異與空間分區(qū)

4.1 協(xié)調(diào)性多維分異

通過前文分析可知，黑龍江各地市土地-人口系統(tǒng)間存在明顯的線性協(xié)調(diào)關(guān)系（如圖3a 所示）。但土地-經(jīng)濟(jì)、人口-經(jīng)濟(jì)、土地-人口-經(jīng)濟(jì)耦合協(xié)調(diào)性出現(xiàn)分異。將各地市人口-經(jīng)濟(jì)、土地-經(jīng)濟(jì)耦合協(xié)調(diào)分析的元數(shù)據(jù)（表2）制作散點分布圖，如圖5b、5c 所示?？梢姼鞯厥性谕恋?經(jīng)濟(jì)、人口-經(jīng)濟(jì)協(xié)調(diào)性的集聚結(jié)果相同。均僅有哈爾濱市位于第一象限（H-H）；均僅有大慶市位于第二象限(L-H)；大興安嶺、雞西、鶴崗、黑河、佳木斯、牡丹江、伊春、雙鴨山、七臺河均位于第三象限(L-L)；齊齊哈爾、綏化均分布于第四象限(H-L)。

通常若元素位于相同象限，揭示他們的實際地情類似，面臨的矛盾趨同，應(yīng)進(jìn)行統(tǒng)一規(guī)劃和管理。由于二維協(xié)調(diào)性產(chǎn)生了兩種不同的分異模式，因此需進(jìn)一步對三維協(xié)調(diào)關(guān)系進(jìn)行探究。各地市人口-土地-經(jīng)濟(jì)三維分異的類型及內(nèi)涵闡釋如表4。

圖3 黑龍江省土地人口經(jīng)濟(jì)二維分異圖 Fig.3 Land population economy clustering quadrant of Heilongjiang province

表4 土地-人口-經(jīng)濟(jì)三維分異類型表 Table 4 Types of land-population-economy system relations

將圖3各地市的土地-人口-經(jīng)濟(jì)水平與表4關(guān)系類型進(jìn)行匹配。可知，僅有哈爾濱市為I 型；齊齊哈爾、綏化為II型；僅有大慶市的為VII 型；大興安嶺、雞西、鶴崗、黑河、佳木斯、牡丹江、伊春、雙鴨山、七臺河的為VIII 型。排除呈線性協(xié)調(diào)關(guān)系的土地-人口系統(tǒng)外，各地區(qū)的土地-人口-經(jīng)濟(jì)的三維協(xié)調(diào)分異和經(jīng)濟(jì)-人口、土地-經(jīng)濟(jì)二維協(xié)調(diào)分異聚類結(jié)果相同，均形成以下4 種模式（表5）。

表5 發(fā)展模式類型 Table 5 Types of development model

4.2 空間分區(qū)及建議

根據(jù)協(xié)調(diào)性分異結(jié)果的發(fā)展模式聚類，第1 類以城鄉(xiāng)協(xié)調(diào)型大慶為代表，其城鎮(zhèn)人口缺口值小，農(nóng)村土地壓力值小，地區(qū)經(jīng)濟(jì)水平高。地區(qū)定位為城鄉(xiāng)協(xié)調(diào)型的主要貢獻(xiàn)因子是高企的地區(qū)生產(chǎn)總值，但其構(gòu)成中石油工業(yè)占比過高，隨著地區(qū)石油資源的耗竭，地區(qū)面臨巨大的產(chǎn)業(yè)轉(zhuǎn)型壓力。因此，該地區(qū)應(yīng)該發(fā)揮現(xiàn)有的資金優(yōu)勢，及時尋找替代性支柱產(chǎn)業(yè)，本文將其劃歸為產(chǎn)業(yè)轉(zhuǎn)型換擋區(qū)。

第2 類以農(nóng)林專業(yè)型大興安嶺、鶴崗、黑河、雞西、佳木斯、牡丹江、七臺河、雙鴨山、伊春等地市為代表。其城鎮(zhèn)人口缺口值小，農(nóng)村土地壓力值小，地區(qū)經(jīng)濟(jì)水平雖低，但與人地系統(tǒng)協(xié)調(diào)性好。因此農(nóng)林發(fā)展區(qū)的各地市，應(yīng)該順應(yīng)地區(qū)自然資源稟賦條件，把生態(tài)環(huán)境保護(hù)看作是地區(qū)經(jīng)濟(jì)發(fā)展的保證，樹立和踐行“綠水青山就是金山銀山”的理念，突出綠色農(nóng)業(yè)特色，向開放型、高效型、生態(tài)型、觀賞型現(xiàn)代化農(nóng)、林、牧、漁業(yè)發(fā)展。本文將其劃歸為生態(tài)農(nóng)業(yè)發(fā)展區(qū)。農(nóng)林專業(yè)型是黑龍江北部和東部地區(qū)的主導(dǎo)發(fā)展模式，對此，應(yīng)通過專業(yè)化分工和協(xié)作，避免不必要的內(nèi)耗，促使整體最優(yōu)化。

第3 類以城鄉(xiāng)差距型哈爾濱市為代表，其城鎮(zhèn)人口缺口值大，農(nóng)村土地壓力值大，地區(qū)經(jīng)濟(jì)水平高。造成其城鄉(xiāng)發(fā)展水平差距的主要原因在于，農(nóng)村土地難以承載如此龐大的人口數(shù)量。因此，要從根本上改變農(nóng)村發(fā)展處境，必須有效地實現(xiàn)農(nóng)村富余勞動力的轉(zhuǎn)移。同時，憑借著哈爾濱市作為省會城市的巨大資源優(yōu)勢，發(fā)揮對資金、物資、信息、技術(shù)等生產(chǎn)要素產(chǎn)生巨大的輻射效應(yīng)和帶動效應(yīng)，帶動整個區(qū)域的發(fā)展，因此本文將其劃歸為城鎮(zhèn)輻射帶動區(qū)。

第4 類以總體滯后型齊齊哈爾、綏化為代表，其城鎮(zhèn)人口缺口值大，農(nóng)村土地壓力值大，地區(qū)經(jīng)濟(jì)水平低。城鄉(xiāng)人地壓力均較大，一方面成為經(jīng)濟(jì)發(fā)展的桎梏，同時也表明地區(qū)蘊藏巨大發(fā)展?jié)摿?，因此地區(qū)必須以工業(yè)化為抓手，培育對經(jīng)濟(jì)支撐作用強(qiáng)的大企業(yè)，進(jìn)而增強(qiáng)城鎮(zhèn)經(jīng)濟(jì)實力、吸納過剩農(nóng)業(yè)人口。同時農(nóng)村地區(qū)要提高土地復(fù)合利用強(qiáng)度，提升土地單位經(jīng)濟(jì)產(chǎn)出，走農(nóng)業(yè)多功能化、農(nóng)村工業(yè)化和人口城鎮(zhèn)化多措并舉以促進(jìn)農(nóng)民增收的新路子，本文將其劃歸為城鄉(xiāng)復(fù)合優(yōu)化區(qū)。

5 結(jié) 論

1）本文得到黑龍江省戶均適宜土地經(jīng)營規(guī)模為 11.6 hm2，過剩農(nóng)業(yè)人口為1 102 萬人，農(nóng)村土地總壓力為51 478 940 hm2；黑龍江省總體城鎮(zhèn)化率飽和值為79.5%，即城鎮(zhèn)人口承載力為3 031 萬，該情況下城鎮(zhèn)人口缺口數(shù)值為789 萬人；各地市的經(jīng)濟(jì)水平呈現(xiàn)多中心分散分布，核密度中心出現(xiàn)在大慶市。土地壓力、人口缺口、經(jīng)濟(jì)水平高值區(qū)均集中在西南部地帶，核密度中心均出現(xiàn)在哈-大-齊工業(yè)走廊地帶。

2）黑龍江省的土地-人口、土地-經(jīng)濟(jì)、人口-經(jīng)濟(jì)系統(tǒng)均達(dá)到了很好的耦合協(xié)調(diào)水平，且耦合協(xié)調(diào)格局的高值區(qū)均出現(xiàn)在省內(nèi)西南部地帶。揭示了地區(qū)應(yīng)以人口流動促進(jìn)土地權(quán)屬流動從而緩解城鄉(xiāng)矛盾的合理性。即促進(jìn)過剩農(nóng)業(yè)人口向城鎮(zhèn)轉(zhuǎn)移，提高城鎮(zhèn)發(fā)展質(zhì)量的同時又可以促進(jìn)土地經(jīng)營規(guī)模的擴(kuò)大，是促進(jìn)城鄉(xiāng)協(xié)調(diào)發(fā)展的主要抓手。

3）土地-人口-經(jīng)濟(jì)系統(tǒng)多維協(xié)調(diào)類型分，顯現(xiàn)了明顯的聚類。第1 類為城鄉(xiāng)協(xié)調(diào)型的大慶，本文將其劃歸為產(chǎn)業(yè)轉(zhuǎn)型換擋區(qū)，應(yīng)抓緊尋找替代性產(chǎn)業(yè)；第2 類包含農(nóng)林發(fā)展型的大興安嶺、鶴崗、黑河、雞西、佳木斯、牡丹江、七臺河、雙鴨山、伊春，本文將其劃歸為生態(tài)農(nóng)業(yè)發(fā)展區(qū)，應(yīng)走生態(tài)綠色發(fā)展的路子，并避免同質(zhì)化競爭；第3 類為城鄉(xiāng)差距型的哈爾濱市，本文將其劃歸為城鎮(zhèn)輻射帶動區(qū)，應(yīng)充分發(fā)揮大城市的帶動作用促進(jìn)周邊地區(qū)發(fā)展；第4 類為總體滯后型的齊齊哈爾、綏化，本文將其劃歸為城鄉(xiāng)復(fù)合優(yōu)化區(qū)，應(yīng)快培育工業(yè)企業(yè)帶動城鎮(zhèn)發(fā)展，農(nóng)業(yè)向多功能化著眼以消化巨量過剩農(nóng)業(yè)人口。