劉書暢，葉艷妹，林耀奔

基于脫鉤理論與LMDI模型的農(nóng)村居民點(diǎn)演化特征及驅(qū)動(dòng)因素分解

劉書暢，葉艷妹※，林耀奔

（1. 浙江大學(xué)公共管理學(xué)院，杭州 310058；2. 城鄉(xiāng)建設(shè)用地節(jié)約集約利用實(shí)驗(yàn)室，北京 100812）

為研究中國農(nóng)村人口與農(nóng)村居民點(diǎn)用地變化協(xié)同關(guān)系和驅(qū)動(dòng)因素，以促進(jìn)農(nóng)村居民點(diǎn)用地集約利用，該文基于農(nóng)村人地變化脫鉤和城鄉(xiāng)人口轉(zhuǎn)型視角，以全國和31個(gè)省級(jí)行政區(qū)域?yàn)檠芯繀^(qū)，將Tapio脫鉤模型和對(duì)數(shù)平均迪式指數(shù)分解法（logarithmic mean Divisia index，LMDI）應(yīng)用于農(nóng)村人地協(xié)調(diào)演化研究及驅(qū)動(dòng)效應(yīng)分解。結(jié)果表明，近20 a來，中國城鄉(xiāng)發(fā)展經(jīng)歷了快速轉(zhuǎn)型，農(nóng)村常住人口持續(xù)減少，降幅約為30.7%，而同期農(nóng)村居民點(diǎn)規(guī)模不減反增，農(nóng)村人地演化關(guān)系失衡，建設(shè)用地趨于粗放利用；中國及各地區(qū)農(nóng)村人口和居民點(diǎn)變化的脫鉤關(guān)系主要表現(xiàn)為強(qiáng)負(fù)脫鉤，其中東部和西部地區(qū)負(fù)脫鉤程度較為嚴(yán)重，且負(fù)脫鉤程度總體上呈加劇趨勢(shì)；1996—2008和2009—2016年期間，中國農(nóng)村居民點(diǎn)規(guī)模分別增加了2.020×105hm2、7.275×105hm2，其中農(nóng)村居民點(diǎn)利用集約度效應(yīng)、城鄉(xiāng)人口結(jié)構(gòu)效應(yīng)、城鎮(zhèn)化效應(yīng)、區(qū)域總?cè)丝谛?yīng)的貢獻(xiàn)量分別為33.497×105、-116.833×105、71.813×105、13.544×105hm2和36.689×105、-68.271×105、32.164×105、6.692×105hm2；不同時(shí)段各效應(yīng)貢獻(xiàn)度不同，1996—2008年城鎮(zhèn)化效應(yīng)大于農(nóng)村居民點(diǎn)利用集約度效應(yīng)，而在2009年之后，農(nóng)村居民點(diǎn)利用集約度效應(yīng)逐漸超過城鎮(zhèn)化效應(yīng)，成為農(nóng)村居民點(diǎn)增加的主導(dǎo)因素。

土地利用；農(nóng)村；模型；農(nóng)村居民點(diǎn)；演化特征；Tapio脫鉤；LMDI

0 引 言

農(nóng)村居民點(diǎn)承載了農(nóng)民的生存保障功能，鄉(xiāng)村人地關(guān)系的變化較多的映射在人口與鄉(xiāng)村聚落上。近年來，中國工業(yè)化、城鎮(zhèn)化進(jìn)程日益加快，城鄉(xiāng)建設(shè)用地?cái)U(kuò)張成為這一進(jìn)程中的顯著特征[1-3]。解釋論上，隨著城鄉(xiāng)人口轉(zhuǎn)型及經(jīng)濟(jì)快速增長，城鄉(xiāng)建設(shè)用地變化應(yīng)表現(xiàn)為城鎮(zhèn)建設(shè)用地的增加和農(nóng)村居民點(diǎn)用地的減少。但現(xiàn)實(shí)中隨著鄉(xiāng)村人口的減少和城鎮(zhèn)用地規(guī)模的擴(kuò)大，農(nóng)村居民點(diǎn)用地卻不減反增[4-5]，這一結(jié)果與一般性假設(shè)相悖。據(jù)統(tǒng)計(jì)，1996—2016年全國農(nóng)村人口凈減少2.6億，約下降30.7%，而同期農(nóng)村居民點(diǎn)用地增長了16.7%。在城鄉(xiāng)用地矛盾日漸突出的現(xiàn)狀下，深入探索農(nóng)村人口與居民點(diǎn)用地變化的協(xié)同關(guān)系及驅(qū)動(dòng)力，促進(jìn)農(nóng)村土地資源的有效配置與高效利用，將是強(qiáng)化耕地保護(hù)、平衡城鄉(xiāng)土地利用和優(yōu)化國土空間格局的關(guān)鍵突破口。

目前，有關(guān)農(nóng)村居民點(diǎn)用地演化的研究主要從以下3個(gè)方面：1）從農(nóng)村居民點(diǎn)利用的數(shù)量規(guī)模展開，研究農(nóng)村居民點(diǎn)用地?cái)U(kuò)展的動(dòng)態(tài)變化，并對(duì)其影響因素進(jìn)行分析。如王介勇等[6]研究了黃淮海平原農(nóng)區(qū)典型村莊用地?cái)U(kuò)展的主要特征，從內(nèi)生動(dòng)力和外部動(dòng)力兩方面探討了村莊用地?cái)U(kuò)張的驅(qū)動(dòng)機(jī)制；劉超等[7]運(yùn)用動(dòng)態(tài)度模型研究了關(guān)中地區(qū)農(nóng)村居民點(diǎn)用地變化的時(shí)空特征，并從非農(nóng)人口和農(nóng)業(yè)機(jī)械動(dòng)力等因素的差異對(duì)農(nóng)村居民點(diǎn)用地動(dòng)態(tài)變化進(jìn)行解釋；譚雪蘭等[8]運(yùn)用回歸分析方法研究了農(nóng)林牧漁總產(chǎn)值、糧食產(chǎn)量和固定資產(chǎn)投資等對(duì)長株潭城市群農(nóng)村居民點(diǎn)變化的作用力。2）從農(nóng)村居民點(diǎn)利用質(zhì)量的角度，開展農(nóng)村居民點(diǎn)集約利用評(píng)價(jià)，探索農(nóng)村居民點(diǎn)用地集約化管控路徑。如曲衍波等[9]基于投入-產(chǎn)出原理對(duì)北京市平谷區(qū)農(nóng)村居民點(diǎn)集約利用程度進(jìn)行評(píng)價(jià)，并分析了影響集約利用的障礙因子；周躍云等[10]從“兩型社會(huì)”建設(shè)背景出發(fā)，分析了株洲市農(nóng)村居民點(diǎn)集約狀況，利用灰色關(guān)聯(lián)分析法進(jìn)行因素探討。3）從城鄉(xiāng)建設(shè)用地協(xié)調(diào)變化的視角，王婧等[11]分析了中國城鄉(xiāng)人地變化的耦合特征，蔡芳芳等[12]以南通市為例，解析了近年來城鄉(xiāng)建設(shè)用地變化的失衡狀態(tài)及成因等。通過分析發(fā)現(xiàn)，既有研究具有視角多元、方法論多元的特征，但主要側(cè)重于某區(qū)域、省份或市（縣）等中微觀尺度和定性分析，缺乏全國宏觀層面對(duì)農(nóng)村居民點(diǎn)用地變化與驅(qū)動(dòng)因素的定量比較研究。

為此，本研究從城鄉(xiāng)人口轉(zhuǎn)型和農(nóng)村人地變化脫鉤關(guān)系的視角，采用Tapio脫鉤模型計(jì)算1996—2016年全國農(nóng)村人口與居民點(diǎn)變化的脫鉤程度，并從省域尺度分析農(nóng)村人地變化的空間差異?；跀U(kuò)展的Kaya恒等式，綜合考慮農(nóng)村居民點(diǎn)利用集約度效應(yīng)、城鄉(xiāng)人口結(jié)構(gòu)效應(yīng)、城鎮(zhèn)化效應(yīng)以及區(qū)域總?cè)丝谛?yīng)4個(gè)因素，建立農(nóng)村居民點(diǎn)規(guī)模變化的因素分解模型，并將對(duì)數(shù)平均迪式指數(shù)分解法（LMDI）引入到農(nóng)村居民點(diǎn)擴(kuò)張的驅(qū)動(dòng)因素分解中來，定量分析各驅(qū)動(dòng)效應(yīng)對(duì)農(nóng)村居民點(diǎn)規(guī)模變化的貢獻(xiàn)度，探索影響農(nóng)村居民點(diǎn)變化的主要因素，以揭示其潛在的作用機(jī)制，為有效管控農(nóng)村居民點(diǎn)規(guī)模和加強(qiáng)農(nóng)村居民點(diǎn)用地規(guī)劃管理提供依據(jù)。

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

1.1 Tapio脫鉤模型

脫鉤（decoupling）是來源于物理領(lǐng)域的概念，可用來分析2個(gè)或多個(gè)具有響應(yīng)關(guān)系的系統(tǒng)時(shí)空變化的耦合狀態(tài)[13-14]。目前，用于脫鉤分析的模型主要有OECD模型[15]、Tapio脫鉤模型[16]、IPAT方程[17]等。Tapio脫鉤模型是由OECD模型發(fā)展演變而來[18]，基于彈性變化，用增長率的比值對(duì)脫鉤狀態(tài)進(jìn)行衡量，且不受研究基期選擇的限制[19-20]，使得脫鉤程度的計(jì)算更加靈活、準(zhǔn)確。本文采用Tapio脫鉤模型測(cè)度農(nóng)村人口和居民點(diǎn)用地變化的脫鉤關(guān)系，模型定義如下

式中為農(nóng)村人口和居民點(diǎn)用地變化的脫鉤彈性值，0、L分別表示基期和第年的農(nóng)村居民點(diǎn)面積，hm2；0、P分別表示基期和第年的農(nóng)村人口，人；為農(nóng)村居民點(diǎn)面積增長率，%，為農(nóng)村人口變化率，%。

對(duì)于不同的研究領(lǐng)域或研究對(duì)象，“脫鉤”程度的定義都會(huì)有所差異，基于本文研究目的，結(jié)合Tapio[16]對(duì)脫鉤類型的劃分以及黃木易等[21]國內(nèi)學(xué)者對(duì)Tapio脫鉤模型的應(yīng)用研究，以脫鉤彈性值0.8和1.2作為劃分依據(jù)，將脫鉤程度劃分為強(qiáng)脫鉤（<0，>0且<0）、增長性脫鉤（0<<0.8，>0且>0）、衰退性脫鉤（>1.2，<0且<0）、強(qiáng)負(fù)脫鉤（<0，<0且>0）、衰退負(fù)脫鉤（0<<0.8，<0且<0）、增長負(fù)脫鉤（>1.2，>0且>0）、增長性耦合（0.8<<1.2，>0且>0）、衰退性耦合（0.8<<1.2，<0且<0）8個(gè)類型。另外，鑒于以往相關(guān)研究中以彈性值0.8和1.2來劃分脫鉤類型，而不便于比較研究單元脫鉤程度的動(dòng)態(tài)變化，同時(shí)考慮到農(nóng)村人口普遍減少的現(xiàn)實(shí)，農(nóng)村人地變化很可能多處于負(fù)脫鉤狀態(tài)，為此，本文對(duì)負(fù)脫鉤類型作了進(jìn)一步細(xì)致劃分，其中強(qiáng)負(fù)脫鉤狀態(tài)進(jìn)一步細(xì)分為“強(qiáng)負(fù)脫鉤Ⅰ、Ⅱ…Ⅵ”，增長負(fù)脫鉤狀態(tài)細(xì)分為“增長負(fù)脫鉤Ⅰ、Ⅱ…Ⅵ”，衰退負(fù)脫鉤狀態(tài)劃分為“衰退負(fù)脫鉤Ⅰ、Ⅱ、Ⅲ、Ⅳ”。負(fù)脫鉤程度由Ⅰ到Ⅵ不斷加劇，以便于比較在脫鉤狀態(tài)不變時(shí)，農(nóng)村人口與居民點(diǎn)用地脫鉤程度的動(dòng)態(tài)變化（表1）。

表1 農(nóng)村人口與農(nóng)村居民點(diǎn)用地變化脫鉤程度劃定

其中，強(qiáng)脫鉤狀態(tài)下，農(nóng)村人口增加而建設(shè)用地趨于減少，這是農(nóng)村人地關(guān)系變化的理想狀態(tài)；增長性脫鉤表示農(nóng)村人口和建設(shè)用地面積均增加，但相比建設(shè)用地?cái)U(kuò)張幅度，農(nóng)村人口增幅更大；衰退性脫鉤意味著農(nóng)村人口和建設(shè)用地面積均減少，但相比農(nóng)村人口減少幅度，建設(shè)用地降幅更明顯。強(qiáng)負(fù)脫鉤意味著隨著農(nóng)村人口的減少，農(nóng)村建設(shè)用地趨于增加，這是較不合理的狀態(tài)；衰退負(fù)脫鉤表示農(nóng)村人口與建設(shè)用地均減少，但建設(shè)用地的降幅小于人口降幅；增長負(fù)脫鉤表示農(nóng)村人口與建設(shè)用地均增加，但建設(shè)用地增幅高于人口增幅。在耦合狀態(tài)下，農(nóng)村人口與建設(shè)用地增幅或降幅相當(dāng)，是介于脫鉤和負(fù)脫鉤的中間狀態(tài)。

1.2 擴(kuò)展的Kaya恒等式

Kaya恒等式是日本的Yoichi Kaya教授提出的，當(dāng)時(shí)用于分析影響CO2排放量的驅(qū)動(dòng)因素[22-23]。本文基于擴(kuò)展的Kaya恒等式，將其運(yùn)用于農(nóng)村居民點(diǎn)變化的驅(qū)動(dòng)力測(cè)度，建立農(nóng)村居民點(diǎn)用地變化的因素分解模型，將農(nóng)村居民點(diǎn)用地變化的影響因素分解為農(nóng)村居民點(diǎn)利用集約度效應(yīng)、城鄉(xiāng)人口結(jié)構(gòu)效應(yīng)、城鎮(zhèn)化效應(yīng)以及區(qū)域總?cè)丝谛?yīng)4個(gè)因素，擴(kuò)展的Kaya恒等式形式如下

式中是農(nóng)村居民點(diǎn)總面積，0、1、2分別是農(nóng)村人口、城鎮(zhèn)人口和區(qū)域總?cè)丝?，人?/0為人均農(nóng)村居民點(diǎn)面積，用以表示農(nóng)村居民點(diǎn)利用的集約度，m2/人；=0/1為鄉(xiāng)村和城鎮(zhèn)人口規(guī)模之比，用以指示城鄉(xiāng)人口規(guī)模變動(dòng)；=1/2表示區(qū)域城鎮(zhèn)化水平。

1.3 LMDI分解模型

LMDI分解法由新加坡學(xué)者Ang[24]提出，其能夠?qū)τ绊懸蛩赝耆纸馇也划a(chǎn)生殘差[25]，主要有加法分解和乘法分解2種模式，并且這2種模式實(shí)質(zhì)上是無差異的[26]。按照擴(kuò)展的Kaya恒等式和LMDI分解模型，研究基期到第年的農(nóng)村居民點(diǎn)用地變化△為總效應(yīng)，其是由農(nóng)村居民點(diǎn)利用集約度效應(yīng)、城鄉(xiāng)人口結(jié)構(gòu)效應(yīng)、城鎮(zhèn)化效應(yīng)以及區(qū)域總?cè)丝谛?yīng)引起的，按照加法分解形式，農(nóng)村居民點(diǎn)用地變化的LMDI分解模式如下：

式中△S為農(nóng)村居民點(diǎn)用地變化量，hm2；△S表示基于加法分解模式的農(nóng)村居民點(diǎn)利用集約度效應(yīng)，用來反映農(nóng)村居民點(diǎn)利用集約度的變化對(duì)農(nóng)村居民點(diǎn)變化的影響；△S表示城鄉(xiāng)人口結(jié)構(gòu)效應(yīng)，反映城鄉(xiāng)人口規(guī)模比值的變化對(duì)農(nóng)村居民點(diǎn)變化的影響；△S表示城鎮(zhèn)化效應(yīng)，用以反映城鎮(zhèn)化水平的變化對(duì)農(nóng)村居民點(diǎn)變化的影響；△S表示區(qū)域總?cè)丝谛?yīng)，用以反映區(qū)域總?cè)丝诘淖兓瘜?duì)農(nóng)村居民點(diǎn)變化的影響。以上4個(gè)驅(qū)動(dòng)因素基于加法分解的計(jì)算公式如下：