省域土地人口經濟多維協調格局及分區發展策略

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

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

0 引 言

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

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

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

1 材料與方法

1.1 數據來源

受2015 年開始的農業供給側結構性改革的影響，研究區糧食價格產生了巨大的波動，且這種影響至今尚未消弭，為保證數據的穩健性，故本文研究期限截止于2015年。本文城鄉社會經濟原始數據來源于2011—2016 年《黑龍江省統計年鑒》、《全國農產品成本收益資料匯編》和《中國住戶調查年鑒》等。本文人口數據來源于1976—2016 年《黑龍江統計年鑒》。本文土地利用數據來源于2016 年土地利用數據庫，提取黑龍江省各地市耕地面積并匯總得到。

1.2 研究方法

1.2.1 土地規模經營模型

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

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

1.2.2 飽和城鎮化率的計算模型

本文采用logistic 模型測度飽和城鎮化率，該方法被聯合國廣泛用于世界各地區城鎮化水平預測，故又稱“聯合國法”。即假設在預測期內城鄉人口增長率差保持不變，根據已知城鄉人口調查的數據，外推預測期末的城鎮人口占總人口比例。其式為：

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

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

1.2.3 地區經濟水平評價方法

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

對各地區經濟發展水平的評價，采用區域發展水平測度通常采用的結節性指數，函數形式如下：

式中Eco 為經濟發展水平結節性指數，Kj為第j 項指標標準化后的分值；Kj為第j 項指標因子權重；n 為指標數目，此處n=6。

1.2.4 秩相關系數

秩相關系數是描述兩要素之間相關程度的統計指標，它將兩要素樣本值位次替代實際數據而求得的統計量，是位次分析的數量化，本文用以測度系統間的總體絕對相關性。

式中R 為秩相關系數，di為要素間位次差，n 為樣本量。

表1 經濟評價指標體系 Table 1 Economic evaluation index system

1.2.5 核密度分析法

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

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

1.2.6 耦合協調度

耦合度協調度借鑒來源于物理學中容量耦合系數模型，是可以用以定量表征多要素間協調程度的指標，本文用以探究系統間各要素的局部協調性，式為

式中D 表示耦合協調度，C 為耦合度，U1表示過剩農村人口，U2表示城鎮化人口缺口，U3表示地區經濟水平。φ、ρ 和β 為權重參數。本文認為土地、人口、經濟方面因素同等重要故φ、ρ、β 取均等值0.333。在耦合協調度的計算中，為使多系統間指標可比，對數據進行極值標準化以去除量綱。

2 人口、土地、經濟系統評價

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

表2 單系統評價結果 Table 2 Single system evaluation results

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

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

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

圖1 土地、人口、經濟系統核密度圖 Fig.1 KDE of land, population and economic system

3 土地-人口-經濟多維協調性分析

秩（復）相關系數的結果表征的是系統間的絕對關系，但無法定位系統內部各要素的具體協調程度。耦合協調度能定位系統內部元素的協調程度，但表征的是相對關系。因此本文首先采用秩(復)相關系數進行總體協調水平定性，進而采用耦合協調度對內部要素間協調水平進行劃分，這種處理方式保證了結果的科學準確。

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

3.1 總體協調性驗證

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

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

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

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

3.2 局部協調性分析

通過總體協調性驗證，可知土地、人口、經濟三系統間均存在顯著的相關，并且劃分出了局部協調性的分級標準，運用耦合協調度模型對三系統間的局部協調性進行測度，結果如表3。為探索局部協調性的空間特征，將表中數據進行空間可視化表達，如圖2 所示。

表3 局部協調性評價結果 Table 3 Micro-coordination evaluation results

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

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

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

圖2 人口、土地、經濟系統耦合協調格局 Fig.2 Types of population, land and economy coupling coordination relations

4 協調性分異與空間分區

4.1 協調性多維分異

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

通常若元素位于相同象限，揭示他們的實際地情類似，面臨的矛盾趨同，應進行統一規劃和管理。由于二維協調性產生了兩種不同的分異模式，因此需進一步對三維協調關系進行探究。各地市人口-土地-經濟三維分異的類型及內涵闡釋如表4。

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

表4 土地-人口-經濟三維分異類型表 Table 4 Types of land-population-economy system relations

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

表5 發展模式類型 Table 5 Types of development model

4.2 空間分區及建議

根據協調性分異結果的發展模式聚類，第1 類以城鄉協調型大慶為代表，其城鎮人口缺口值小，農村土地壓力值小，地區經濟水平高。地區定位為城鄉協調型的主要貢獻因子是高企的地區生產總值，但其構成中石油工業占比過高，隨著地區石油資源的耗竭，地區面臨巨大的產業轉型壓力。因此，該地區應該發揮現有的資金優勢，及時尋找替代性支柱產業，本文將其劃歸為產業轉型換擋區。

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

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

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

5 結 論

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

2）黑龍江省的土地-人口、土地-經濟、人口-經濟系統均達到了很好的耦合協調水平，且耦合協調格局的高值區均出現在省內西南部地帶。揭示了地區應以人口流動促進土地權屬流動從而緩解城鄉矛盾的合理性。即促進過剩農業人口向城鎮轉移，提高城鎮發展質量的同時又可以促進土地經營規模的擴大，是促進城鄉協調發展的主要抓手。

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