異質(zhì)性勞動(dòng)力集聚對(duì)區(qū)域碳排放的影響研究

赫永達(dá)，林伯強(qiáng)

（1.山西財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)學(xué)院，山西 太原 030006；2.廈門(mén)大學(xué) 中國(guó)能源政策研究院，福建 廈門(mén) 361005）

一、引言

在全球氣候變暖趨勢(shì)加劇和節(jié)能減排政策的驅(qū)動(dòng)下，我國(guó)正積極履行碳減排承諾，力爭(zhēng)在2030年前實(shí)現(xiàn)碳達(dá)峰，2060年?duì)幦?shí)現(xiàn)碳中和。然而，我國(guó)的能源要素投入總體上仍表現(xiàn)為高碳化特征，并且這種以煤炭和石油等化石能源消耗為主的能源使用結(jié)構(gòu)在相當(dāng)長(zhǎng)時(shí)期內(nèi)并不會(huì)發(fā)生根本性的改變。［1］尤其是伴隨著中國(guó)經(jīng)濟(jì)由高速發(fā)展向高質(zhì)量發(fā)展的轉(zhuǎn)型，勞動(dòng)力大規(guī)模跨區(qū)域流動(dòng)，在加速了地區(qū)城鎮(zhèn)化和工業(yè)化的同時(shí)也誘發(fā)了地區(qū)資源短缺、環(huán)境惡化、碳減排政策失效等問(wèn)題。如何改善勞動(dòng)力、能源等生產(chǎn)要素的市場(chǎng)配置效率已成為中國(guó)經(jīng)濟(jì)實(shí)現(xiàn)綠色低碳化發(fā)展所亟須解決的問(wèn)題。因此，研究異質(zhì)性勞動(dòng)力跨區(qū)域流動(dòng)與地區(qū)碳排放之間的影響機(jī)制及其時(shí)空演化態(tài)勢(shì)對(duì)我國(guó)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)與低碳減排協(xié)調(diào)發(fā)展具有重要意義。

目前學(xué)術(shù)界對(duì)勞動(dòng)人口和碳排放之間關(guān)系進(jìn)行研究主要從兩個(gè)角度切入。第一種研究思路以STIRPAT模型為理論基礎(chǔ)分析人口結(jié)構(gòu)對(duì)碳排放的影響。在早期研究人口結(jié)構(gòu)和碳排放的研究中，學(xué)者主要集中討論人口密度和規(guī)模對(duì)人口結(jié)構(gòu)的影響，［2-4］并進(jìn)一步延伸至人口結(jié)構(gòu)與碳排放之間的關(guān)系。該種思路需要先驗(yàn)設(shè)定人口規(guī)模對(duì)碳排放是單位彈性的。如果人口規(guī)模對(duì)碳排放的彈性不滿足單位彈性假設(shè)，則以此為基礎(chǔ)的環(huán)境庫(kù)茲涅茲曲線不再成立。遺憾的是這樣的先驗(yàn)設(shè)定并沒(méi)有獲得任何實(shí)證檢驗(yàn)，尤其在人力資本積累達(dá)到一定程度時(shí)，不同技能結(jié)構(gòu)的勞動(dòng)力集聚及其引致的地區(qū)資本深化程度是不同的。［5-6］換言之，“人口規(guī)模對(duì)碳排放的單位彈性假設(shè)”會(huì)隨著技能型勞動(dòng)力規(guī)模效應(yīng)的體現(xiàn)而不再成立。由此，Richard York 等為克服該局限提出了STIRPAT 模型，認(rèn)為人口、技術(shù)和財(cái)富對(duì)碳排放是非固定比例的隨機(jī)影響。［7］然而，基于STIRPAT 模型框架的研究，對(duì)單純的時(shí)間序列、不同國(guó)家的面板樣本、中國(guó)不同地區(qū)的省級(jí)面板樣本進(jìn)行實(shí)證檢驗(yàn)時(shí)，目前并沒(méi)有就勞動(dòng)力規(guī)模、勞動(dòng)力結(jié)構(gòu)等人口因素對(duì)碳排放的影響形成一致的結(jié)論。［8-11］

第二種研究思路則是從產(chǎn)業(yè)結(jié)構(gòu)角度出發(fā)。由于地區(qū)產(chǎn)業(yè)結(jié)構(gòu)差異決定了地區(qū)勞動(dòng)力雇傭結(jié)構(gòu)的不同，整體而言工業(yè)勞動(dòng)碳排放水平必然高于農(nóng)業(yè)和服務(wù)業(yè)碳排放水平。因此，這類文獻(xiàn)的聚焦點(diǎn)在于揭示異質(zhì)性勞動(dòng)力集聚對(duì)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)變遷產(chǎn)生的影響，認(rèn)為碳排放變化表面上來(lái)看是產(chǎn)業(yè)結(jié)構(gòu)變遷造成，實(shí)際上則由勞動(dòng)力結(jié)構(gòu)決定。［12-13］該種研究思想得到了眾多實(shí)證研究者的驗(yàn)證。如孫振清等測(cè)算了不同產(chǎn)業(yè)結(jié)構(gòu)背景下，創(chuàng)新勞動(dòng)力和技術(shù)創(chuàng)新水平對(duì)碳排放效率的影響；［14］王勇等、張同斌等、朱佩譽(yù)和凌文基于投入產(chǎn)出模型，模擬了不同產(chǎn)業(yè)結(jié)構(gòu)下中國(guó)碳排放水平和碳達(dá)峰路徑；［15-17］周亞軍和吉萍利用空間計(jì)量模型驗(yàn)證了產(chǎn)業(yè)升級(jí)過(guò)程中地區(qū)間要素配置及其空間溢出效應(yīng)對(duì)碳排放影響。［18］

通過(guò)對(duì)上述文獻(xiàn)的梳理，發(fā)現(xiàn)關(guān)于勞動(dòng)人口和碳排放之間關(guān)系的研究存在以下兩點(diǎn)仍值得進(jìn)一步探討。首先，現(xiàn)有文獻(xiàn)往往忽視了勞動(dòng)力流動(dòng)這一典型特征事實(shí)。中國(guó)改革開(kāi)放和加入WTO 以來(lái)，地區(qū)間經(jīng)濟(jì)差距逐漸顯現(xiàn)，其中一個(gè)明顯的特征是東部沿海地區(qū)依靠先發(fā)優(yōu)勢(shì)，率先完成了經(jīng)濟(jì)增長(zhǎng)由傳統(tǒng)的“要素驅(qū)動(dòng)”到“技術(shù)模仿驅(qū)動(dòng)”再到“創(chuàng)新驅(qū)動(dòng)”的逐步過(guò)渡。其造成的直接后果就是勞動(dòng)力從農(nóng)村向城市、從內(nèi)陸地區(qū)向沿海地區(qū)大規(guī)模轉(zhuǎn)移，且人口紅利逐漸由豐富轉(zhuǎn)向匱乏再到消失（見(jiàn)圖1 和圖2）。按照現(xiàn)有關(guān)于勞動(dòng)力規(guī)模和碳排放水平之間關(guān)系研究的結(jié)論推斷，地區(qū)碳排放水平的提高理應(yīng)隨著人口紅利的消失而停止。然而上述預(yù)期在實(shí)際生產(chǎn)過(guò)程中并沒(méi)有發(fā)生，勞動(dòng)力流動(dòng)規(guī)模與地區(qū)碳排放之間究竟是單調(diào)的線性關(guān)系還是復(fù)雜的非線性關(guān)系，現(xiàn)有文獻(xiàn)尚未對(duì)此給出確切的答案，尤其是結(jié)合現(xiàn)有關(guān)于產(chǎn)業(yè)結(jié)構(gòu)和碳排放之間關(guān)系的研究成果來(lái)看，勞動(dòng)力流動(dòng)常常對(duì)應(yīng)著產(chǎn)業(yè)結(jié)構(gòu)的變遷。由此，揭示異質(zhì)性勞動(dòng)力集聚及其造成的地區(qū)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整是指明地區(qū)人口結(jié)構(gòu)和碳排放關(guān)系的關(guān)鍵因素之一。

圖1 地區(qū)城鎮(zhèn)勞動(dòng)數(shù)量（萬(wàn)人）

圖2 地區(qū)就業(yè)水平（萬(wàn)人）

其二，對(duì)于地區(qū)勞動(dòng)力結(jié)構(gòu)和碳排放之間關(guān)系的研究，目前鮮有文獻(xiàn)關(guān)注異質(zhì)性勞動(dòng)力流動(dòng)這一關(guān)鍵事實(shí)。從異質(zhì)性勞動(dòng)力流動(dòng)與地區(qū)產(chǎn)業(yè)結(jié)構(gòu)和碳排放水平之間關(guān)系來(lái)看，異質(zhì)性勞動(dòng)力集聚對(duì)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整的影響不同，技能型勞動(dòng)力集聚必然促進(jìn)產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化，無(wú)論是能源利用率還是能源使用水平均有助于碳減排；從異質(zhì)性勞動(dòng)力流動(dòng)與地區(qū)技術(shù)水平和碳排放水平關(guān)系來(lái)看，技能型勞動(dòng)一方面能夠在較大程度上促進(jìn)技術(shù)進(jìn)步，進(jìn)而反向推動(dòng)能源利用效率的提高，另一方面能夠誘發(fā)技術(shù)性要素對(duì)能源、非技術(shù)性生產(chǎn)要素之間的替代，進(jìn)而減少能源的使用，控制碳排放水平；從異質(zhì)性勞動(dòng)力流動(dòng)與勞動(dòng)力結(jié)構(gòu)和碳排放水平之間關(guān)系來(lái)看，技能型勞動(dòng)力與技術(shù)進(jìn)步、要素替代、碳減排之間的內(nèi)生反饋機(jī)制并不是無(wú)限循環(huán)的。一方面技能型勞動(dòng)力集聚會(huì)提高地區(qū)生活成本，廠商必然會(huì)擴(kuò)大資本和能源的使用，平衡由技能型勞動(dòng)力集聚造成的簡(jiǎn)單生產(chǎn)要素（低成本勞動(dòng)力）的擠出，另一方面，產(chǎn)業(yè)相依性決定了技能型勞動(dòng)力和高技術(shù)產(chǎn)業(yè)并不能夠脫離低端產(chǎn)業(yè)而集聚。［5］［19］最終，異質(zhì)性要素集聚和碳排放之間的關(guān)系應(yīng)是上述幾種作用力的相對(duì)均衡。

鑒于上述分析，本文擬構(gòu)造兩地區(qū)、三部門(mén)的一般均衡模型，從理論上揭示異質(zhì)性勞動(dòng)力集聚對(duì)碳排放的作用機(jī)制。本文可能實(shí)現(xiàn)的創(chuàng)新之處在于：首先，通過(guò)區(qū)分高碳排產(chǎn)業(yè)和傳統(tǒng)服務(wù)業(yè)，將異質(zhì)性勞動(dòng)力流動(dòng)和資本流動(dòng)方程納入一般均衡框架，構(gòu)建能夠描述異質(zhì)性勞動(dòng)力集聚特征的一般均衡模型；其次，本文對(duì)理論模型的關(guān)鍵結(jié)論進(jìn)行參數(shù)校準(zhǔn)和數(shù)值模擬，結(jié)合理論框架和數(shù)值模擬結(jié)果，從理論上闡明異質(zhì)性勞動(dòng)力流動(dòng)和地區(qū)碳排放之間的作用機(jī)制和演化特征；第三，本文采用兩區(qū)制空間杜賓模型，結(jié)合中國(guó)省級(jí)宏觀數(shù)據(jù)，對(duì)異質(zhì)性勞動(dòng)力流動(dòng)過(guò)程中可能存在的異質(zhì)性空間溢出效應(yīng)及其對(duì)地區(qū)碳排放的影響進(jìn)行量化檢驗(yàn)。通過(guò)上述工作，力圖對(duì)異質(zhì)性勞動(dòng)力集聚和地區(qū)碳排放之間的內(nèi)在機(jī)制和實(shí)際運(yùn)行規(guī)律做出關(guān)鍵性闡述。

二、碳排放理論模型設(shè)定

如前文所述，解釋異質(zhì)性勞動(dòng)力流動(dòng)與地區(qū)碳排放之間內(nèi)生機(jī)制的關(guān)鍵在于厘清異質(zhì)性勞動(dòng)力流動(dòng)及其引致的資本流動(dòng)速度之間的相互關(guān)系，由此判斷地區(qū)資本深化演化趨勢(shì)，最終傳導(dǎo)到要素配置決定的要素密集型產(chǎn)業(yè)結(jié)構(gòu)調(diào)整規(guī)律。由此，本文在Dissou和Horvath等學(xué)者模型基礎(chǔ)上構(gòu)建了包含代表性家庭、低排放行業(yè)（傳統(tǒng)服務(wù)業(yè)）、高排放行業(yè)（現(xiàn)代工業(yè)）的三部門(mén)一般均衡模型。［20-21］具體而言，將勞動(dòng)力劃分為技能型勞動(dòng)力、非技能型勞動(dòng)力，其中非技能型勞動(dòng)力從事普通生產(chǎn)和服務(wù)業(yè)生產(chǎn)，技能型勞動(dòng)力更多地作為技術(shù)載體。這樣做的目的在于，不同勞動(dòng)力對(duì)應(yīng)著不同的資本——?jiǎng)趧?dòng)力替代規(guī)則，而地區(qū)資本深化則是決定碳排放的關(guān)鍵因素。由此，通過(guò)對(duì)異質(zhì)性勞動(dòng)力和資本流動(dòng)的刻畫(huà)，映射到地區(qū)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和地區(qū)碳排放結(jié)構(gòu)動(dòng)態(tài)變遷，以此揭示地區(qū)間異質(zhì)性勞動(dòng)力流動(dòng)和碳排放之間的內(nèi)生傳導(dǎo)機(jī)理。

考慮一個(gè)包含兩地區(qū)（地區(qū)s和n），且工業(yè)品（CMt）和服務(wù)產(chǎn)品（CBt）共存的經(jīng)濟(jì)體。工業(yè)品廠商雇傭技能型勞動(dòng)力（Ht）、非技能型勞動(dòng)力（Lt）、使用可變資本（Kt）和能源（Et）進(jìn)行生產(chǎn)，服務(wù)產(chǎn)品廠商雇傭簡(jiǎn)單勞動(dòng)力（Jt）、使用固定資本（F）進(jìn)行生產(chǎn)。勞動(dòng)力市場(chǎng)的價(jià)格信號(hào)（工資差距）決定了異質(zhì)性勞動(dòng)力的跨區(qū)域流動(dòng)，資本市場(chǎng)的價(jià)格信號(hào)（利率差距）決定了資本的跨區(qū)域流動(dòng)。

1.碳排放

石油、煤炭、天然氣等能源在生產(chǎn)消耗過(guò)程中會(huì)排放CO2?？紤]不同種類的能源在消耗過(guò)程中排放的CO2數(shù)量不同，因此參照張同斌等的做法，［17］設(shè)定煤炭在生產(chǎn)過(guò)程中的排放系數(shù)為θcoal，石油和天然氣在生產(chǎn)過(guò)程中的排放系數(shù)為θoil，則工業(yè)碳排放emt方程可表示為：

emt表示t期工業(yè)碳排放總量表示工業(yè)使用煤炭等能源總量表示工業(yè)使用原油、天然氣等能源總量；θcoal和θoil則表示相應(yīng)能源在生產(chǎn)過(guò)程中的碳排放轉(zhuǎn)換系數(shù)。

2.高碳排行業(yè)生產(chǎn)行為

設(shè)定每個(gè)地區(qū)只存在一個(gè)典型加總的工業(yè)廠商。本文設(shè)定工業(yè)生產(chǎn)技術(shù)滿足C-D 生產(chǎn)函數(shù)形式，并以資本和勞動(dòng)份額參數(shù)之間的差別描述地區(qū)間不同要素密集程度的差異：

為s（s∈A,B）地區(qū)工業(yè)廠商在t期的總產(chǎn)出，Kst、Hst、Lst和Est分別表示該地區(qū)工業(yè)廠商使用的總資本，雇傭的可流動(dòng)技能型勞動(dòng)力和可流動(dòng)非技能型勞動(dòng)力總量以及使用的能源總量；At表示地區(qū)的技術(shù)水平；as、bs、cs和（1-as-bs-cs）分別表示資本、技能型勞動(dòng)、能源和非技能型勞動(dòng)對(duì)工業(yè)產(chǎn)出的貢獻(xiàn)份額。

3.低碳排放廠商生產(chǎn)行為

根據(jù)傳統(tǒng)D-S框架，假設(shè)服務(wù)產(chǎn)品生產(chǎn)表現(xiàn)出顯著的規(guī)模報(bào)酬遞增特點(diǎn)，且滿足壟斷競(jìng)爭(zhēng)條件。服務(wù)產(chǎn)品廠商雇傭可流動(dòng)勞動(dòng)力Jit以及使用相應(yīng)的固定投資Fi提供不同種類的服務(wù)產(chǎn)品。除此之外，設(shè)定每個(gè)地區(qū)的市場(chǎng)存在qt個(gè)服務(wù)產(chǎn)品廠商，每個(gè)廠商生產(chǎn)存在一定差異的服務(wù)產(chǎn)品。假設(shè)每個(gè)服務(wù)產(chǎn)品廠商采用線性生產(chǎn)技術(shù)進(jìn)行生產(chǎn)，每個(gè)廠商的產(chǎn)出xit與使用的固定資本Fi和可變勞動(dòng)力Jit滿足線性生產(chǎn)技術(shù)：

Jit為廠商i雇傭的普通勞動(dòng)力數(shù)量；Fi為單個(gè)廠商所使用的固定資本投入（假設(shè)不存在資本折舊）；aT為產(chǎn)品轉(zhuǎn)化效率。

設(shè)定地區(qū)存在qt個(gè)差異性服務(wù)產(chǎn)品廠商。根據(jù)D-S條件，每個(gè)服務(wù)產(chǎn)品廠商生產(chǎn)的服務(wù)產(chǎn)品數(shù)量cit以及居民消費(fèi)的“打包”服務(wù)產(chǎn)品數(shù)量CBt之間關(guān)系滿足：

ρ為每種服務(wù)產(chǎn)品之間的替代彈性。根據(jù)居民支出的約束條件，得到居民購(gòu)買(mǎi)“打包”服務(wù)產(chǎn)品和每個(gè)服務(wù)產(chǎn)品廠商制定的出廠價(jià)格之間關(guān)系為：

假設(shè)服務(wù)產(chǎn)品廠商i提供xit單位的服務(wù)產(chǎn)品，均衡條件下地區(qū)居民對(duì)服務(wù)產(chǎn)品總量的需求等于地區(qū)服務(wù)產(chǎn)品總供給，即：

σ＞0 為相對(duì)風(fēng)險(xiǎn)規(guī)避系數(shù)①關(guān)于σ數(shù)更加詳細(xì)的說(shuō)明見(jiàn)后文居民偏好部分。數(shù)理模型推導(dǎo)可向作者索要。；Ot為地區(qū)總?cè)丝贠t=Ht+Lt+Jt。將公式（5）關(guān)于服務(wù)產(chǎn)品價(jià)格的關(guān)系式代入公式（6），求得服務(wù)產(chǎn)品廠商i供給xi,t為：

當(dāng)商品完全同質(zhì)不存在差別情況下，pi=p，xi=x。同時(shí)，由于M商品可自由貿(mào)易，則M商品價(jià)格滿足一價(jià)定律。出于計(jì)算簡(jiǎn)便考慮，假設(shè)同一市場(chǎng)內(nèi)部不存在商品運(yùn)輸成本，則PM價(jià)格可作為價(jià)格基準(zhǔn)，即PM=1。那么地區(qū)單個(gè)服務(wù)產(chǎn)品廠商的產(chǎn)出可進(jìn)一步化簡(jiǎn)為：

關(guān)于服務(wù)產(chǎn)品廠商制定的出廠價(jià)格和支付勞動(dòng)工資的求解，根據(jù)廠商0 利潤(rùn)條件，廠商的全部利潤(rùn)用于支付普通勞動(dòng)工人的工資，即Rt=ptxt=ptct=·Jt。將廠商產(chǎn)出方程（8）代入0 利潤(rùn)條件，可得廠商支付的工資水平為：

4.異質(zhì)性勞動(dòng)力流動(dòng)與碳排放

在實(shí)際要素流動(dòng)過(guò)程中，要素的流動(dòng)規(guī)模取決于地區(qū)間相應(yīng)要素收入差距以及地區(qū)間生活成本差距。假設(shè)地區(qū)間生活成本以地區(qū)間不可貿(mào)易的服務(wù)產(chǎn)品價(jià)格、可貿(mào)易的工業(yè)品價(jià)格和不可貿(mào)易的住房?jī)r(jià)格表示。定義異質(zhì)性勞動(dòng)力的實(shí)際工資差異形式為：

表示s地區(qū)技能型勞動(dòng)力的實(shí)際工資水平表示s地區(qū)非技能型勞動(dòng)力的實(shí)際工資水平；

表示s地區(qū)普通勞動(dòng)力的實(shí)際工資水平。而表示s地區(qū)技能型勞動(dòng)力的名義工資水平；表示s地區(qū)非技能型勞動(dòng)力的名義工資水平；表示s地區(qū)普通勞動(dòng)力的名義工資水平。技能型勞動(dòng)力和非技能型勞動(dòng)力的名義工資水平分別為其對(duì)產(chǎn)出的邊際貢獻(xiàn)，普通勞動(dòng)力的名義工資水平由（9）式?jīng)Q定。由此，異質(zhì)性勞動(dòng)力跨地區(qū)流動(dòng)的決策決定于地區(qū)間實(shí)際工資差異，［22］那么流入地區(qū)（不妨標(biāo)記s地區(qū)為流入地區(qū)）異質(zhì)性勞動(dòng)力的動(dòng)態(tài)積累方程分別為：

h為技能型勞動(dòng)力的自然增長(zhǎng)率，n為非技能勞動(dòng)力和簡(jiǎn)單勞動(dòng)力的自然增長(zhǎng)率；ψh、ψl、ψj分別為技能型勞動(dòng)力、非技能型勞動(dòng)力和普通勞動(dòng)力對(duì)地區(qū)間相應(yīng)勞動(dòng)報(bào)酬率缺口的反應(yīng)敏感系數(shù)；為流出地區(qū)在t時(shí)期技能型勞動(dòng)力數(shù)量，為流出地區(qū)在t時(shí)期非技能型勞動(dòng)力數(shù)量為流出地區(qū)在t時(shí)期普通勞動(dòng)力數(shù)量。上述三個(gè)勞動(dòng)轉(zhuǎn)移方程表明：地區(qū)間勞動(dòng)的實(shí)際工資差距（名義工資和價(jià)格水平）是決定相應(yīng)勞動(dòng)力流動(dòng)的根本因素。

與勞動(dòng)力流動(dòng)方程類似，資本的轉(zhuǎn)移同樣受到兩地區(qū)間利率差異的影響。與勞動(dòng)力流動(dòng)所不同的是，由于利率本身代表貨幣的價(jià)格，因此資本是否轉(zhuǎn)移取決于地區(qū)間名義利率。則資本動(dòng)態(tài)方程的化簡(jiǎn)形式為：

gk為資本積累速度；ψk為資本對(duì)地區(qū)名義利率缺口的反應(yīng)敏感系數(shù)為流入地市場(chǎng)利率，為流出地市場(chǎng)利率水平；為流出地資本數(shù)量。

流出地要素積累方程與流入地要素積累方程對(duì)應(yīng)，只需要將自然增長(zhǎng)率后面的轉(zhuǎn)移部分改為負(fù)號(hào)即可。［5］

結(jié)合碳排放方程（1）以及要素流動(dòng)方程（13）、（14）、（15），能夠得到地區(qū)碳排放總量和異質(zhì)性勞動(dòng)力流動(dòng)之間的關(guān)系為：

式（17）、式（18）即為以技能型勞動(dòng)力和非技能型勞動(dòng)力數(shù)量表示的流入地區(qū)碳排放水平的函數(shù)①流出地區(qū)只需將分子第二項(xiàng)加號(hào)改為減號(hào)即可。。從上述方程來(lái)看，發(fā)現(xiàn)技能型勞動(dòng)力和非技能型勞動(dòng)力對(duì)地區(qū)碳排放影響機(jī)制存在差異：技能型勞動(dòng)力更多通過(guò)單位勞動(dòng)力碳排放的垂直效應(yīng)來(lái)影響地區(qū)碳排放水平（技能型勞動(dòng)力工資水平和地區(qū)間工資差異較大），而非技能型勞動(dòng)力則主要通過(guò)碳排放的規(guī)模效應(yīng)影響地區(qū)（非技能型勞動(dòng)力工資水平較低但規(guī)模較大）。通過(guò)對(duì)比地區(qū)碳排放和異質(zhì)性勞動(dòng)力之間關(guān)系，進(jìn)一步能夠發(fā)現(xiàn)：有效采用技能型勞動(dòng)力和資本對(duì)非技能型勞動(dòng)力進(jìn)行替代，能夠極大地降低碳排放的規(guī)模效應(yīng)。

5.居民偏好

出于簡(jiǎn)單且不失一般性的考慮，本文假設(shè)不同群體的居民均具有不變跨期替代彈性的消費(fèi)偏好，其根據(jù)市場(chǎng)利率和工資水平提供相應(yīng)的資本和勞動(dòng)力。每類代表性家庭均追求終生效用最大化。為了便于分析，設(shè)定每個(gè)代表家庭均供給一單位相應(yīng)勞動(dòng)力獲得工資收入，并根據(jù)市場(chǎng)價(jià)格在儲(chǔ)蓄和消費(fèi)之間權(quán)衡。因此，不同群體在相應(yīng)收支約束條件下的效用函數(shù)①除了要素流動(dòng)方程外，流入地和流出地的其他所有方程形式完全一致。因此為了寫(xiě)作的便利和模型形式的簡(jiǎn)潔，本文所建立的模型，除需特殊備注說(shuō)明，所有方程只寫(xiě)出流入地相應(yīng)方程的表達(dá)形式，并略去各變量表示地區(qū)的下標(biāo)?？杀硎緸椋?/p>

和分別為d居民群體在t期對(duì)不可貿(mào)易的服務(wù)產(chǎn)品、可貿(mào)易的工業(yè)品和住房的消費(fèi)數(shù)量；為d居民投資水平為相應(yīng)的資本存量；為d居民群體獲得的名義勞動(dòng)收入。由于假定每個(gè)居民固定供給1單位勞動(dòng)，因此每個(gè)居民的勞動(dòng)收入等于市場(chǎng)工資率；rt為市場(chǎng)利率；β為主觀貼現(xiàn)因子；σ＞0為相對(duì)風(fēng)險(xiǎn)規(guī)避系數(shù)；γ為實(shí)物資產(chǎn)的效用偏好彈性。在收支約束下最大化居民終身效用，可求得居民d對(duì)相應(yīng)商品的需求為：

三、參數(shù)校準(zhǔn)與數(shù)值模擬

為揭示異質(zhì)性勞動(dòng)力流動(dòng)與地區(qū)碳排放之間的關(guān)系，針對(duì)第二部分的理論模型，根據(jù)求解出的均衡條件進(jìn)行數(shù)值模擬。揭示異質(zhì)性勞動(dòng)力和地區(qū)碳排放之間的關(guān)鍵機(jī)制在于厘清因異質(zhì)性勞動(dòng)力流動(dòng)造成的要素配置變化（資本深化）以及產(chǎn)業(yè)結(jié)構(gòu)變化（高碳排放產(chǎn)業(yè)結(jié)構(gòu)變化）。對(duì)此本節(jié)對(duì)異質(zhì)性勞動(dòng)力流動(dòng)特征、異質(zhì)性勞動(dòng)力流動(dòng)造成的產(chǎn)業(yè)結(jié)構(gòu)變遷和資本配置分別進(jìn)行數(shù)值模擬，并對(duì)異質(zhì)性勞動(dòng)力流動(dòng)造成的地區(qū)間碳排放規(guī)模演變特征進(jìn)行數(shù)值計(jì)算。數(shù)值模擬的主要目的在于揭示異質(zhì)性勞動(dòng)力流動(dòng)這一典型的內(nèi)生經(jīng)濟(jì)現(xiàn)象與產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、碳排放結(jié)構(gòu)變遷之間存在的必然聯(lián)系，對(duì)定量分析并不作過(guò)分苛刻的要求。

1.參數(shù)校準(zhǔn)

為了在極大程度上使本文設(shè)定的理論模型模擬結(jié)果與中國(guó)經(jīng)濟(jì)發(fā)展現(xiàn)實(shí)相契合，本文對(duì)理論模型中的參數(shù)和狀態(tài)變量的初始賦值按照以下規(guī)則校準(zhǔn)。

關(guān)于工業(yè)生產(chǎn)函數(shù)的參數(shù)校準(zhǔn)，要素產(chǎn)出彈性可以根據(jù)要素收入份額占總產(chǎn)出比率近似求得。本文細(xì)分了異質(zhì)性勞動(dòng)力，將本科及以上學(xué)歷作為技能型勞動(dòng)、其他學(xué)歷水平的勞動(dòng)作為非技能型勞動(dòng)力和普通勞動(dòng)力（技能型勞動(dòng)力和普通勞動(dòng)力假設(shè)為同質(zhì)，只不過(guò)雇傭?qū)ο蟛煌?。參照傅曉霞和吳利學(xué)的研究成果，［23-24］以此校準(zhǔn)的流入地資本產(chǎn)出彈性as=0.43，流出地資本產(chǎn)出彈性an=0.41；流入地技能型勞動(dòng)力產(chǎn)出彈性bs=0.24，流出地技能型勞動(dòng)力產(chǎn)出彈性bn=0.21；對(duì)于能源對(duì)工業(yè)產(chǎn)出份額參數(shù)的校準(zhǔn)，通過(guò)投入產(chǎn)出表中第一象限基本流量表的原始數(shù)據(jù)，計(jì)算能源對(duì)工業(yè)的中間投入系數(shù)c=0.17；流入地非技能型勞動(dòng)力產(chǎn)出彈性（1-as-bs-c）=0.18，流出地非技能型勞動(dòng)力產(chǎn)出彈性（1-an-bn-c）=0.15。

關(guān)于CO2排放系數(shù)，借鑒聯(lián)合國(guó)政府間氣候變化委員會(huì)（IPCC）公布的測(cè)算方法、陳詩(shī)一和張同斌等的測(cè)算方法，［17］［25］對(duì)不同能源燃燒利用率和排放系數(shù)進(jìn)行加權(quán)處理計(jì)算，校準(zhǔn)得到煤炭行業(yè)的碳排放系數(shù)θcoal為2.8①1公斤煤炭燃燒后排放2.8公斤CO2，石油、天然氣的排放系數(shù)意義同此。，石油和天然氣行業(yè)的綜合碳排放系數(shù)θoil為2.1。

關(guān)于服務(wù)業(yè)②具體參數(shù)校準(zhǔn)過(guò)程中，校準(zhǔn)的“服務(wù)業(yè)”具體包括運(yùn)輸和郵電、住宿和餐飲業(yè)、旅游業(yè)、教育、衛(wèi)生、社會(huì)服務(wù)、公共管理社會(huì)保障及其他服務(wù)業(yè)。需要指明的是，本文所指的服務(wù)業(yè)旨在表明該行業(yè)提供的服務(wù)不可跨地區(qū)貿(mào)易。因此，諸如運(yùn)輸和郵電業(yè)，盡管該行業(yè)存在跨地區(qū)業(yè)務(wù)往來(lái)，但其提供的物流服務(wù)本質(zhì)上也是本地服務(wù)的一種。生產(chǎn)函數(shù)的參數(shù)校準(zhǔn)，根據(jù)相應(yīng)行業(yè)產(chǎn)出、勞動(dòng)力數(shù)量以及固定資本形成總額，設(shè)定地區(qū)服務(wù)業(yè)固定資產(chǎn)使用總額Fs=Fn=300，產(chǎn)品轉(zhuǎn)化效率aTs=aTn=0.3。兩地區(qū)之間各個(gè)服務(wù)產(chǎn)品廠商生產(chǎn)技術(shù)和生產(chǎn)規(guī)模并不存在差異，影響服務(wù)產(chǎn)品出廠價(jià)格的主要因素源自地區(qū)對(duì)不可貿(mào)易的服務(wù)產(chǎn)品的需求。［26-27］此外，根據(jù)目前研究的一致取值方式，設(shè)定異質(zhì)性消費(fèi)品之間的替代參數(shù)ρ=4。

對(duì)于本文居民偏好等參數(shù)校準(zhǔn)，首先設(shè)定主觀貼現(xiàn)因子β=0.99；根據(jù)相關(guān)研究以及經(jīng)濟(jì)事實(shí)經(jīng)驗(yàn)，設(shè)定資本的穩(wěn)態(tài)增長(zhǎng)率按照固定資產(chǎn)投資增長(zhǎng)率的平均數(shù)取值為0.15；［28］由于居民的消費(fèi)偏好參數(shù)難以直接獲取或校準(zhǔn)，絕大多數(shù)關(guān)于發(fā)達(dá)國(guó)家的宏觀研究中風(fēng)險(xiǎn)規(guī)避系數(shù)設(shè)定為0.5，由此求得居民個(gè)體的主觀貼現(xiàn)率取值在2%～4%之間?？紤]中國(guó)居民儲(chǔ)蓄傾向較高，因此本文設(shè)定中國(guó)家庭的風(fēng)險(xiǎn)規(guī)避系數(shù)為σ=0.55。

本文要素流動(dòng)方程的參數(shù)校準(zhǔn)可分為兩部分。對(duì)于異質(zhì)性勞動(dòng)力的自然增長(zhǎng)率的校準(zhǔn)，由于本文設(shè)定的模型區(qū)分了勞動(dòng)力技能水平存在異質(zhì)性，近似地認(rèn)為普通勞動(dòng)力的自然增長(zhǎng)率約等于地區(qū)戶籍人口的自然增長(zhǎng)率，因此對(duì)于非技能型勞動(dòng)力和技能型勞動(dòng)力，近似地認(rèn)為其平均增長(zhǎng)率n=0.01；針對(duì)技能型勞動(dòng)力的測(cè)算③在針對(duì)技能型勞動(dòng)自然增長(zhǎng)率測(cè)算過(guò)程中發(fā)現(xiàn)平均數(shù)較高主要源于2009-2013 年高校擴(kuò)招所致。技能型勞動(dòng)近些年的增長(zhǎng)率趨于穩(wěn)定，流入地區(qū)的增長(zhǎng)率穩(wěn)定在2%左右，而流出地的增長(zhǎng)率則不足1%。由于本模型的目的并不在于做精確的數(shù)值模擬，而是通過(guò)構(gòu)建合理的模型，對(duì)地區(qū)間勞動(dòng)流動(dòng)所展現(xiàn)出的規(guī)律做理論上的分析，因此本文在數(shù)值模擬過(guò)程中，取2003-2018 年間技能型勞動(dòng)的平均數(shù)作為技能型勞動(dòng)的增長(zhǎng)率。，將本科及以上學(xué)歷水平的勞動(dòng)力作為技能型勞動(dòng)力，通過(guò)對(duì)地區(qū)高校招生人數(shù)的近似估算，技能型勞動(dòng)力的平均增長(zhǎng)率h=0.02。而對(duì)于異質(zhì)性勞動(dòng)力跨區(qū)域流動(dòng)系數(shù)的校準(zhǔn)，本文按照最小偏誤的原則，設(shè)定技能型勞動(dòng)力對(duì)相應(yīng)工資差距的流動(dòng)敏感系數(shù)為ψh=0.05、非技能型勞動(dòng)力對(duì)相應(yīng)工資差距的流動(dòng)敏感系數(shù)為ψl=0.03、普通勞動(dòng)力對(duì)相應(yīng)工資差距的流動(dòng)敏感系數(shù)為ψj=0.03、資本對(duì)利率差距的流動(dòng)敏感系數(shù)為ψk=0.05。

關(guān)于生產(chǎn)要素的初始值，按照各地區(qū)統(tǒng)計(jì)局發(fā)布的統(tǒng)計(jì)年鑒以及《勞動(dòng)與就業(yè)統(tǒng)計(jì)年鑒》公布的數(shù)據(jù)，取流入地區(qū)技能型勞動(dòng)力初始值=140、流出地區(qū)技能型勞動(dòng)力初始值=150；流入地區(qū)非技能型勞動(dòng)力初始值=1 400、流出地區(qū)非技能型勞動(dòng)力初始值=2 000；流入地區(qū)普通勞動(dòng)力初始值=1 400、流出地區(qū)普通勞動(dòng)力初始值=6 000（單位均為萬(wàn)人）；關(guān)于資本存量的初始值，參照張軍利用永續(xù)盤(pán)存法，［29］計(jì)算出2004 年流入地工業(yè)資本存量=48 000、流出地工業(yè)資本存量=73 000（單位均為億元）。

上述參數(shù)的校準(zhǔn)以及相應(yīng)變量初始值取值，所涉及的基礎(chǔ)數(shù)據(jù)源于中國(guó)統(tǒng)計(jì)局發(fā)布的《中國(guó)投入產(chǎn)出表-2018》《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)能源統(tǒng)計(jì)年鑒》《中國(guó)人口和就業(yè)統(tǒng)計(jì)年鑒》以及《中經(jīng)網(wǎng)統(tǒng)計(jì)數(shù)據(jù)庫(kù)》。

2.數(shù)值模擬

本文通過(guò)模擬地區(qū)產(chǎn)業(yè)結(jié)構(gòu)（地區(qū)工業(yè)產(chǎn)出占總產(chǎn)出的比重）、地區(qū)工業(yè)產(chǎn)出比（流入地和流出地工業(yè)產(chǎn)出比）、地區(qū)要素流動(dòng)規(guī)模（技能型勞動(dòng)力、非技能型勞動(dòng)力和資本流動(dòng)規(guī)模）、地區(qū)勞動(dòng)技能結(jié)構(gòu)（各地區(qū)技能型勞動(dòng)力和非技能型勞動(dòng)力比值）和地區(qū)人均資本深化程度（各地區(qū)勞均資本存量）等指標(biāo)，對(duì)地區(qū)間異質(zhì)性勞動(dòng)力流動(dòng)和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整之間規(guī)律進(jìn)行刻畫(huà)。

本文首先模擬了地區(qū)間要素流動(dòng)特征（見(jiàn)圖3）。從技能型勞動(dòng)力流動(dòng)特征來(lái)看，技能型勞動(dòng)力從技術(shù)劣勢(shì)地區(qū)向技術(shù)優(yōu)勢(shì)地區(qū)的流動(dòng)規(guī)模會(huì)逐漸減少，在可以預(yù)見(jiàn)的未來(lái)（40 期模擬期后），地區(qū)間技能型勞動(dòng)力甚至可能出現(xiàn)回流的趨勢(shì)。造成這樣的原因在于：首先，技術(shù)優(yōu)勢(shì)地區(qū)的資本回報(bào)率更高，因此在規(guī)模報(bào)酬不變的前提下，工業(yè)勞動(dòng)對(duì)產(chǎn)出的貢獻(xiàn)更小。這也就造成了逐利行為會(huì)使得資本大規(guī)模從技術(shù)劣勢(shì)地區(qū)向技術(shù)優(yōu)勢(shì)地區(qū)流動(dòng)，資本流動(dòng)一方面造成了地區(qū)工業(yè)進(jìn)程加快，在產(chǎn)業(yè)相依性的約束下，若流入地區(qū)服務(wù)業(yè)無(wú)法同規(guī)模擴(kuò)張，則工業(yè)擴(kuò)張的一個(gè)必然結(jié)果是地區(qū)服務(wù)業(yè)價(jià)格上升、生活成本提高，從而對(duì)地區(qū)工業(yè)勞動(dòng)力和服務(wù)業(yè)勞動(dòng)力產(chǎn)生擠出效應(yīng)；另一方面，資本規(guī)模擴(kuò)張以及生活成本的提高會(huì)促使資本對(duì)工業(yè)勞動(dòng)和能源使用的替代。在擠出效應(yīng)和要素的替代效應(yīng)共同作用下，最終會(huì)促使資本流向技術(shù)先進(jìn)地區(qū)，且對(duì)勞動(dòng)力和能源產(chǎn)生一定程度的擠出效應(yīng)。

圖3 地區(qū)要素流動(dòng)特征

其次，在對(duì)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)的模擬中發(fā)現(xiàn)地區(qū)間工業(yè)會(huì)向技術(shù)要素集聚的地區(qū)遷移。具體來(lái)看，對(duì)于資本和技能型勞動(dòng)力集聚的地區(qū)而言，資本的集聚會(huì)導(dǎo)致地區(qū)產(chǎn)業(yè)結(jié)構(gòu)向資本密集型調(diào)整，這種特征直觀地反映在流入地的工業(yè)結(jié)構(gòu)上，技能型勞動(dòng)力集聚和資本密集地區(qū)工業(yè)占比逐漸提高，而資本和技能型勞動(dòng)力流失的地區(qū)工業(yè)占比則會(huì)因工業(yè)要素流失的雙重疊加而產(chǎn)生萎縮。要素流失地區(qū)只能通過(guò)簡(jiǎn)單擴(kuò)大生產(chǎn)規(guī)模來(lái)阻滯實(shí)體經(jīng)濟(jì)的萎縮，并不能夠通過(guò)技術(shù)進(jìn)步和效率提升來(lái)實(shí)現(xiàn)產(chǎn)業(yè)升級(jí)，從要素投入視角上來(lái)看，技能型勞動(dòng)力流失的地區(qū)只能通過(guò)提高能源使用和簡(jiǎn)單要素投入來(lái)維持工業(yè)份額，此舉不僅不利于工業(yè)升級(jí)，同時(shí)也會(huì)造成碳排放水平提高。

圖4 地區(qū)產(chǎn)業(yè)結(jié)構(gòu)

圖5 地區(qū)勞動(dòng)技能結(jié)構(gòu)

圖6 地區(qū)勞均資本存量（萬(wàn)元/人）

圖7 地區(qū)碳排放總量（億噸）

第三，本文模擬了地區(qū)間勞動(dòng)力技能結(jié)構(gòu)和地區(qū)間勞動(dòng)力資本配置結(jié)構(gòu)。模擬結(jié)果表明：對(duì)于流入地而言，在工業(yè)勞動(dòng)力集聚的早期，因地區(qū)間產(chǎn)業(yè)相依性引致的服務(wù)業(yè)勞動(dòng)力流動(dòng)規(guī)模會(huì)快于工業(yè)勞動(dòng)力的流動(dòng)規(guī)模，并在一定程度上降低地區(qū)勞動(dòng)力平均技能水平。因此在模擬初期發(fā)現(xiàn)流入地區(qū)勞動(dòng)力技能結(jié)構(gòu)存在微弱的下降趨勢(shì)。而伴隨著簡(jiǎn)單勞動(dòng)力逐漸集聚并形成規(guī)模，則會(huì)導(dǎo)致技術(shù)集聚地區(qū)的勞動(dòng)力技能結(jié)構(gòu)持續(xù)提高。與之對(duì)應(yīng)的是地區(qū)間勞動(dòng)力的資本配置狀況。由于資本對(duì)地區(qū)間利率差異的敏感性要顯著高于勞動(dòng)力對(duì)地區(qū)間工資差異的敏感性，因此資本流動(dòng)規(guī)模的變動(dòng)速度會(huì)顯著快于勞動(dòng)力流動(dòng)規(guī)模的變動(dòng)速度，這就導(dǎo)致了地區(qū)間勞均資本存量顯著提高，資本深化促使了地區(qū)間工業(yè)集聚。

最后，本文模擬了地區(qū)碳排放總量和地區(qū)勞均碳排放量的演化路徑。從碳排放總量上來(lái)看，由于流入地在經(jīng)濟(jì)體量上和流出地存在一定差距，因此流入地碳排放總量較流出地而言存在一定差距。但從相對(duì)變化來(lái)看，伴隨著流入地技能型勞動(dòng)力集聚和技術(shù)進(jìn)步，流入地區(qū)的碳排放總量增長(zhǎng)趨勢(shì)較流出地而言更加緩慢，地區(qū)碳排放比呈現(xiàn)出下降的趨勢(shì)。這表明要素集聚，尤其是技能型勞動(dòng)力的集聚會(huì)對(duì)地區(qū)化石能源使用存在一定程度的替代，技術(shù)進(jìn)步會(huì)提高化石能源的利用效率，進(jìn)而促進(jìn)技術(shù)進(jìn)步和碳減排同時(shí)發(fā)生。

上述結(jié)論在勞均碳排放這一指標(biāo)的演化趨勢(shì)上更加顯著。從圖8看出，無(wú)論是對(duì)于勞動(dòng)力流入地，還是對(duì)于勞動(dòng)力流出地而言，地區(qū)勞均碳排放量均呈現(xiàn)下降的趨勢(shì)。對(duì)要素流出地而言，要素流失造成地區(qū)行業(yè)規(guī)模收縮，不僅減少了能源的使用，同時(shí)在碳減排政策的影響下，廠商和行業(yè)出于成本最小化原則，會(huì)利用其他要素對(duì)能源進(jìn)行替代，進(jìn)而導(dǎo)致能源使用量降低的程度要大于勞動(dòng)力流失的程度，造成流出地勞均碳排放下降。對(duì)于流入地而言，一方面勞動(dòng)力集聚會(huì)促使能源利用效率的提高，另一方面技能型勞動(dòng)力及其引致的資本會(huì)對(duì)能源產(chǎn)生更大程度的替代，綜合表現(xiàn)為流入地勞均碳排放呈現(xiàn)下降的趨勢(shì)，且較流出地而言，碳減排程度和效果更好。

圖8 地區(qū)勞均碳排放量（千噸/萬(wàn)人）

四、異質(zhì)性勞動(dòng)力流動(dòng)與地區(qū)碳排放的實(shí)證檢驗(yàn)

1.變量構(gòu)造及數(shù)據(jù)選取

上文理論方程的建立以及數(shù)值模擬的結(jié)果顯示：地區(qū)間異質(zhì)性勞動(dòng)力流動(dòng)（尤其是技能型勞動(dòng)力的流動(dòng)）對(duì)地區(qū)間要素相對(duì)配置和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整存在顯著影響。由此，揭示異質(zhì)性勞動(dòng)力流動(dòng)對(duì)地區(qū)碳排放之間的關(guān)系，本質(zhì)上在于檢驗(yàn)技能型勞動(dòng)力流動(dòng)是否能夠促進(jìn)地區(qū)間技術(shù)進(jìn)步和能源利用效率的提高以及人力資本和物質(zhì)資本對(duì)能源的替代。因此，本文利用空間計(jì)量的方法對(duì)理論模型揭示的異質(zhì)性要素配置的內(nèi)生機(jī)制以及數(shù)值模擬推演的關(guān)鍵結(jié)論進(jìn)行實(shí)證檢驗(yàn)。

本文核心被解釋變量為碳排放，選取地區(qū)CO2排放總量（千噸）作為相應(yīng)的度量指標(biāo)。關(guān)于核心解釋變量的選取，針對(duì)理論模型的關(guān)鍵結(jié)論，以“地區(qū)年末常住人口”與“地區(qū)年末戶籍人口”之間的差，比上“地區(qū)年末戶籍人口”，作為勞動(dòng)力凈遷入比率（migit）；［5］［30］選取地區(qū)本科及以上學(xué)歷勞動(dòng)力數(shù)量作為衡量技能型勞動(dòng)力的指標(biāo)（skillit）；選取城鎮(zhèn)單位就業(yè)人員總量、私營(yíng)企業(yè)就業(yè)人員和個(gè)體就業(yè)人員總量（萬(wàn)人），并扣除技能型勞動(dòng)力作為衡量地區(qū)勞動(dòng)力數(shù)量的指標(biāo)（labit）；選取工業(yè)增加值占GDP份額作為產(chǎn)業(yè)結(jié)構(gòu)（struit）的代理指標(biāo)；選取化石能源使用量占能源總使用量的份額作為能源結(jié)構(gòu)（feit）的代理指標(biāo)。

考慮數(shù)據(jù)完整性，本文所采用的指標(biāo)橫跨的時(shí)間序列為2004-2019 年。其中，碳排放數(shù)據(jù)源于Wind數(shù)據(jù)庫(kù)；化石能源使用量相應(yīng)數(shù)據(jù)源于《中國(guó)能源統(tǒng)計(jì)年鑒》；工業(yè)增加值和GDP數(shù)據(jù)在具體實(shí)證過(guò)程中均經(jīng)過(guò)價(jià)格平減，相應(yīng)數(shù)據(jù)源于《中經(jīng)網(wǎng)統(tǒng)計(jì)數(shù)據(jù)庫(kù)》；技能型勞動(dòng)力數(shù)據(jù)源于《中國(guó)人口和就業(yè)統(tǒng)計(jì)年鑒》。詳細(xì)的指標(biāo)含義、處理過(guò)程見(jiàn)表1。

表1 主要變量定義及測(cè)量方式

2.空間相關(guān)性檢驗(yàn)

盡管前文理論模型從機(jī)制上識(shí)別出了技能型勞動(dòng)力流動(dòng)導(dǎo)致的能源利用效率提高以及要素配置狀態(tài)的改變。但在異質(zhì)性勞動(dòng)力流動(dòng)過(guò)程中技術(shù)是否存在溢出效應(yīng)以及不同地區(qū)技能型勞動(dòng)力技術(shù)溢出效應(yīng)是否相同是決定各地區(qū)碳排放向特定水平收斂的關(guān)鍵。因此，在以空間計(jì)量手段檢驗(yàn)異質(zhì)性勞動(dòng)力跨區(qū)域流動(dòng)與地區(qū)碳排放之間關(guān)系的過(guò)程中，首先要解決的問(wèn)題是檢驗(yàn)各個(gè)地區(qū)是否存在空間性。對(duì)此，采用經(jīng)典的空間相關(guān)莫蘭指數(shù)對(duì)地區(qū)間可能存在的空間相關(guān)性進(jìn)行驗(yàn)證：

其中，S2=，xi為第i個(gè)地區(qū)的觀測(cè)值，在本文中xi具體代表該地區(qū)的人均CO2排放水平，n為地區(qū)總數(shù)；wij為空間矩陣W中的元素。莫蘭指數(shù)的取值范圍為［-1，1］，數(shù)值越大表示地區(qū)間的空間正相關(guān)性越強(qiáng)；數(shù)值越小則表示地區(qū)間負(fù)相關(guān)性越強(qiáng)；當(dāng)莫蘭指數(shù)接近0值時(shí)，表明地區(qū)間不存在顯著的空間相關(guān)性。本文選用經(jīng)濟(jì)距離矩陣計(jì)算莫蘭指數(shù)，在后續(xù)的實(shí)證檢驗(yàn)過(guò)程中也將經(jīng)濟(jì)距離矩陣作為主要的空間權(quán)重矩陣緊縮實(shí)證檢驗(yàn)，并通過(guò)替換不同的空間權(quán)重矩陣作為穩(wěn)健性檢驗(yàn)的方式。全局自相關(guān)莫蘭指數(shù)的具體計(jì)算結(jié)果見(jiàn)表2。

表2 莫蘭指數(shù)時(shí)間變化趨勢(shì)

根據(jù)計(jì)算的全局莫蘭指數(shù)結(jié)果，從時(shí)間趨勢(shì)來(lái)看，全國(guó)范圍內(nèi)的空間相關(guān)性逐漸減弱，顯著性也逐漸減弱。同時(shí)一個(gè)有趣的現(xiàn)象是計(jì)算所得的莫蘭指數(shù)均為正值，表明相鄰地區(qū)之間表現(xiàn)出較為一致的碳排放特征，空間鄰近的省份或地區(qū)的碳排放結(jié)構(gòu)類似。為進(jìn)一步揭示不同地區(qū)的碳排放結(jié)構(gòu)和勞動(dòng)力流動(dòng)特征，本文采取LISA分析的方式來(lái)描述不同地區(qū)勞動(dòng)力和碳排放的集聚特征①具體的結(jié)果染色圖可向作者索要。。

從各省勞動(dòng)力的地理分布來(lái)看，表現(xiàn)出顯著的High-High 集聚特征的俱樂(lè)部仍然為位于東南沿海地區(qū)的江蘇，位于珠三角地區(qū)的廣東，以及西部的中心地區(qū)四川、重慶；而表現(xiàn)出明顯虹吸效應(yīng)的High-Low 集聚特征的地區(qū)主要為北京、上海和西安，這三個(gè)地區(qū)勞動(dòng)力數(shù)量的增加并未對(duì)周?chē)貐^(qū)形成有明顯的技術(shù)溢出，也未與周?chē)貐^(qū)形成積極的知識(shí)交流；此外，大部分中部地區(qū)與東部地區(qū)之間存在顯著的負(fù)向空間自相關(guān)，表現(xiàn)出顯著的Low-High 集聚特征；最后，大部分西部地區(qū)仍然表現(xiàn)出勞動(dòng)力數(shù)量低水平的Low-Low集聚特征。而從時(shí)間變動(dòng)趨勢(shì)來(lái)看，勞動(dòng)力的流動(dòng)逐漸向發(fā)達(dá)地區(qū)集中，即勞動(dòng)力逐漸向東南沿海地區(qū)、北京地區(qū)和川渝地區(qū)集中，而中部地區(qū)并未與東部地區(qū)形成較好的知識(shí)交流，表現(xiàn)出顯著的負(fù)向空間自相關(guān)，此外西部地區(qū)的勞動(dòng)力密度則更加稀疏。這些時(shí)間特征均表明勞動(dòng)力，尤其是技能型勞動(dòng)力，逐漸在技術(shù)優(yōu)勢(shì)地區(qū)形成了集聚效應(yīng)。

從各省碳排放的地理分布來(lái)看，各省的碳排放水平仍然在珠三角、長(zhǎng)三角地區(qū)形成了較強(qiáng)的空間關(guān)聯(lián)，表現(xiàn)出明顯的High-High 集聚特征；北京地區(qū)的高技能勞動(dòng)力對(duì)周邊地區(qū)存在顯著的虹吸效應(yīng)，表現(xiàn)為顯著的High-Low 集聚特征；中部地區(qū)的碳排放水平仍然表現(xiàn)出本地水平較低、周邊水平較高的Low-High 的空間負(fù)相關(guān)集聚特征；絕大部分西部地區(qū)表現(xiàn)出顯著的低水平的俱樂(lè)部集聚現(xiàn)象。從時(shí)間趨勢(shì)上來(lái)看，全國(guó)范圍內(nèi)碳排放水平逐漸集聚于東南沿海地區(qū)和北京地區(qū)，中、西部地區(qū)的碳排放水平逐漸減少且未與周邊地區(qū)形成較強(qiáng)的空間聯(lián)動(dòng)。

將上述結(jié)果與圖1 和圖2 所顯示的規(guī)律對(duì)比，發(fā)現(xiàn)目前勞動(dòng)力和碳排放的空間集聚特征極為相似。勞動(dòng)力集聚和城市碳排放集聚之間存在較強(qiáng)的耦合關(guān)系。這初步表明地區(qū)產(chǎn)業(yè)結(jié)構(gòu)（高碳排放的工業(yè)）會(huì)隨著勞動(dòng)力的集聚而集聚。在地區(qū)工業(yè)擴(kuò)大生產(chǎn)規(guī)模過(guò)程中，勞動(dòng)力流動(dòng)和勞動(dòng)力集聚對(duì)地區(qū)能源集聚存在顯著的引致效應(yīng)。而技能型勞動(dòng)力集聚是否會(huì)反向推動(dòng)能源利用技術(shù)升級(jí)，本文進(jìn)一步采取空間計(jì)量的手段對(duì)此進(jìn)行驗(yàn)證。

3.碳排放與勞動(dòng)力流動(dòng)的實(shí)證檢驗(yàn)

結(jié)合上文的理論模型可以看出，一方面，技能型勞動(dòng)力的集聚程度會(huì)促使地區(qū)資本深化，并促進(jìn)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)向資本密集型產(chǎn)業(yè)調(diào)整，推動(dòng)地區(qū)現(xiàn)代產(chǎn)業(yè)（工業(yè)）集聚的發(fā)生，而地區(qū)工業(yè)集聚是造成地區(qū)碳排放水平提高的決定因素。另一方面，技能型勞動(dòng)力集聚對(duì)技術(shù)進(jìn)步的主要驅(qū)動(dòng)因素，不僅能夠替代資本和簡(jiǎn)單勞動(dòng)力，同時(shí)能夠反向推動(dòng)能源利用技術(shù)升級(jí)，提高能源利用效率的同時(shí)降低能源消耗總量。另外，前文理論模型發(fā)現(xiàn)產(chǎn)業(yè)相依性決定了高碳排行業(yè)的發(fā)展不能脫離低碳排行業(yè)的支持，技能型勞動(dòng)力集聚不僅對(duì)能源和簡(jiǎn)單勞動(dòng)力存在擠出效應(yīng)，同時(shí)產(chǎn)業(yè)相依性決定了技能型勞動(dòng)力集聚對(duì)能源之間存在引致效用。最終異質(zhì)性勞動(dòng)力流動(dòng)和碳排放之間關(guān)系取決于上述三種作用力的相對(duì)均衡。為驗(yàn)證理論模型的結(jié)論，本文構(gòu)建如下空間杜賓面板模型（SDM）①本文首先采用LM 檢驗(yàn)，結(jié)果表明SEM 模型、SAR 模型均適用（LM-lag 檢驗(yàn)、R-LM-lag 檢驗(yàn)、LM-err 檢驗(yàn)、R-LM-err 檢驗(yàn)均通過(guò)1%的顯著性水平檢驗(yàn)）；其次，LR 檢驗(yàn)顯著地拒絕了原假設(shè)，說(shuō)明SDM 模型不會(huì)簡(jiǎn)化為SEM 模型或SEM 模型；再次，Hausman檢驗(yàn)結(jié)果支持固定效應(yīng)；最后，聯(lián)合顯著性檢驗(yàn)拒絕了原假設(shè)，說(shuō)明選用時(shí)間和個(gè)體雙重固定的SDM模型較為合適。對(duì)此進(jìn)行檢驗(yàn)。

表3 地區(qū)碳排放與異質(zhì)性勞動(dòng)力流動(dòng)的單區(qū)制空間計(jì)量檢驗(yàn)

從上述實(shí)證結(jié)果可以看出，無(wú)論采用何種空間權(quán)重矩陣，各變量的符號(hào)均未發(fā)生系統(tǒng)性變化，表明本文采用空間杜賓模型所得到的回歸結(jié)果是穩(wěn)健可信的。就各個(gè)變量估計(jì)值所傳遞的具體含義來(lái)看，首先，簡(jiǎn)單勞動(dòng)力總量的集聚、勞動(dòng)力流動(dòng)比率與地區(qū)碳排放之間存在一定的正向關(guān)系，表明地區(qū)勞動(dòng)力總量增加會(huì)引致工業(yè)規(guī)模擴(kuò)大，從而提高碳排放水平。類似的經(jīng)濟(jì)傳導(dǎo)機(jī)制在地區(qū)產(chǎn)業(yè)結(jié)構(gòu)（stru）、化石類能源使用比率（fe）的系數(shù)上也有體現(xiàn)：當(dāng)?shù)貐^(qū)工業(yè)規(guī)模擴(kuò)張、使用化石能源比率上升，則必然增加地區(qū)碳排放水平。

地區(qū)勞動(dòng)力流動(dòng)比率和技能型勞動(dòng)力對(duì)地區(qū)碳排放存在一定程度的負(fù)向影響，這是實(shí)證結(jié)果的一個(gè)重要發(fā)現(xiàn)。從前文理論機(jī)制可以判斷出技能型勞動(dòng)力的集聚在生產(chǎn)過(guò)程中不僅存在著對(duì)其他生產(chǎn)要素的替代，同時(shí)會(huì)反向推動(dòng)能源利用技術(shù)升級(jí)，二者作用下造成技能型勞動(dòng)力集聚對(duì)碳排放產(chǎn)生負(fù)向影響。

前文理論模型一個(gè)重要的引申結(jié)論是：生產(chǎn)技術(shù)水平是決定地區(qū)間產(chǎn)業(yè)結(jié)構(gòu)、能源利用的關(guān)鍵，進(jìn)而對(duì)碳排放存在直接影響。因此，本文采用空間計(jì)量來(lái)檢驗(yàn)勞動(dòng)力流動(dòng)過(guò)程中可能存在的技術(shù)溢出效應(yīng)。然而一個(gè)尚未突破的問(wèn)題在于從不同地區(qū)流出的勞動(dòng)力其空間溢出效應(yīng)可能存在異質(zhì)性，由此造成地區(qū)異質(zhì)性勞動(dòng)力集聚和碳排放特征也會(huì)存在差異。因此，本文采用兩區(qū)制空間杜賓模型，對(duì)不同地區(qū)勞動(dòng)力流動(dòng)可能造成的異質(zhì)性空間溢出效應(yīng)進(jìn)行細(xì)致區(qū)分和討論。具體的兩區(qū)制空間杜賓模型構(gòu)造如下：

其中，d為虛擬變量。結(jié)合不同地區(qū)產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)水平和勞動(dòng)力流動(dòng)特征，本文設(shè)定北京、天津、上海、江蘇、浙江、福建和廣東作為區(qū)制1，其余地區(qū)作為區(qū)制2，以此揭示不同地區(qū)勞動(dòng)力流動(dòng)過(guò)程中可能存在的異質(zhì)性空間溢出效應(yīng)。具體計(jì)算結(jié)果見(jiàn)表4，同時(shí)表5 匯報(bào)了各變量的空間溢出效應(yīng)。

表4 地區(qū)碳排放與異質(zhì)性勞動(dòng)力流動(dòng)的雙區(qū)制空間計(jì)量檢驗(yàn)

表5 地區(qū)碳排放與異質(zhì)性勞動(dòng)力流動(dòng)的空間溢出效應(yīng)

分析兩區(qū)制空間杜賓模型的結(jié)果發(fā)現(xiàn)，整體上，兩區(qū)制空間杜賓模型和單區(qū)制空間杜賓模型所表現(xiàn)出的經(jīng)濟(jì)規(guī)律是類似的：無(wú)論是表示工業(yè)生產(chǎn)規(guī)模的勞動(dòng)力總量（lab）、勞動(dòng)力流動(dòng)比率（mig），抑或是代表地區(qū)工業(yè)化程度的產(chǎn)業(yè)結(jié)構(gòu)（stru），或者是地區(qū)化石能源消耗（fe），整體而言均會(huì)對(duì)地區(qū)碳排放水平產(chǎn)生顯著的正向影響；而技能型勞動(dòng)力集聚，整體而言則會(huì)對(duì)地區(qū)碳排放產(chǎn)生不同顯著程度的負(fù)向影響。

進(jìn)一步對(duì)比不同區(qū)制所反映出的規(guī)律，發(fā)現(xiàn)對(duì)于技術(shù)水平較為先進(jìn)、技能型勞動(dòng)力集聚程度較高的地區(qū)（d=1），上述變量的影響系數(shù)要顯著大于技術(shù)水平較為落后、技能型勞動(dòng)力集聚程度較低的地區(qū)（d=0），尤其是對(duì)于技能型勞動(dòng)力集聚程度而言。由于ρd=1·d＞ρd=0·(1-d)，當(dāng)技能型勞動(dòng)力向技術(shù)優(yōu)勢(shì)、產(chǎn)業(yè)結(jié)構(gòu)高端的地區(qū)集聚時(shí)，較高的技術(shù)水平存量更能發(fā)揮技能型勞動(dòng)力的研發(fā)效應(yīng)和替代效應(yīng)，因此對(duì)能源和其他要素的替代作用更強(qiáng)。

五、結(jié)論

本文通過(guò)構(gòu)建兩地區(qū)三部門(mén)的一般均衡模型，并對(duì)理論模型的關(guān)鍵結(jié)論進(jìn)行數(shù)值模擬，從理論層面揭示了異質(zhì)性勞動(dòng)力流動(dòng)和地區(qū)碳排放之間的內(nèi)在機(jī)制。隨后，本文利用我國(guó)省際宏觀數(shù)據(jù)對(duì)理論模型所闡釋的理論機(jī)制進(jìn)行實(shí)證檢驗(yàn)，得到如下結(jié)論：

首先，本文從理論上揭示了技能型勞動(dòng)力集聚會(huì)通過(guò)提高能源利用水平，以及對(duì)能源和其他生產(chǎn)要素的替代，從而實(shí)現(xiàn)二氧化碳排放的減少。具體而言，在理論模型建立過(guò)程中，通過(guò)識(shí)別不同勞動(dòng)力的異質(zhì)性碳排放特征，刻畫(huà)了工業(yè)勞動(dòng)力與傳統(tǒng)服務(wù)業(yè)勞動(dòng)力流動(dòng)過(guò)程中存在的聯(lián)動(dòng)機(jī)制。這種聯(lián)動(dòng)機(jī)制從總量上決定了地區(qū)碳排放的整體水平，從結(jié)構(gòu)上決定了技能型勞動(dòng)力對(duì)能源的替代程度，從而實(shí)現(xiàn)技術(shù)進(jìn)步和碳減排。最終二者之間的均衡共同決定了地區(qū)碳排放總量。

其次，本文采用數(shù)值模擬的方式，針對(duì)理論模型揭示的異質(zhì)性勞動(dòng)力流動(dòng)和碳排放之間的內(nèi)在機(jī)制進(jìn)行演化分析。從地區(qū)產(chǎn)業(yè)結(jié)構(gòu)演變這一視角來(lái)看，技能型勞動(dòng)力的集聚會(huì)帶來(lái)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化；從要素的配置角度來(lái)看，技能型勞動(dòng)力的集聚一方面對(duì)非技能型勞動(dòng)力產(chǎn)生了替代，另一方面對(duì)資本流動(dòng)存在引致作用，二者共同作用造成了地區(qū)資本深化，進(jìn)一步促使產(chǎn)業(yè)結(jié)構(gòu)向資本密集型和技能密集型演化；由于資本密集型產(chǎn)業(yè)結(jié)構(gòu)和技能密集型產(chǎn)業(yè)結(jié)構(gòu)對(duì)地區(qū)碳排放的作用機(jī)制相反，從現(xiàn)階段中國(guó)地區(qū)經(jīng)濟(jì)發(fā)展模型和碳減排的結(jié)果來(lái)看，由技能集聚反向推動(dòng)能源技術(shù)進(jìn)步，進(jìn)而促使碳減排的效應(yīng)要強(qiáng)于資本深化和規(guī)模擴(kuò)張帶來(lái)的高碳排效應(yīng)，最終促使地區(qū)平均碳排放水平降低。

最后，本文對(duì)理論模型的關(guān)鍵結(jié)論以及數(shù)值模擬對(duì)經(jīng)濟(jì)演化的基本判斷進(jìn)行實(shí)證檢驗(yàn)。本文分別采用單區(qū)制空間杜賓模型和雙區(qū)制空間杜賓模型，對(duì)異質(zhì)性勞動(dòng)力流動(dòng)過(guò)程中可能存在的異質(zhì)性溢出效應(yīng)進(jìn)行甄別和檢驗(yàn)。結(jié)果表明工業(yè)規(guī)模擴(kuò)張、生產(chǎn)要素投入以及化石能源的使用對(duì)地區(qū)碳排放水平存在不同顯著水平的正向影響，而技能型勞動(dòng)力集聚則對(duì)地區(qū)碳減排存在顯著影響。同時(shí)，在技能型勞動(dòng)力集聚和技術(shù)水平更高的地區(qū)更容易發(fā)揮技能型勞動(dòng)力對(duì)能源利用技術(shù)進(jìn)步的反向推動(dòng)作用，技能型勞動(dòng)力的碳減排效應(yīng)更強(qiáng)。

就本文的研究而言，有效疏導(dǎo)異質(zhì)性勞動(dòng)力的流動(dòng)，在吸引技能型勞動(dòng)力集聚、促進(jìn)產(chǎn)業(yè)升級(jí)的同時(shí)，加快實(shí)現(xiàn)地區(qū)經(jīng)濟(jì)增長(zhǎng)由傳統(tǒng)的“要素驅(qū)動(dòng)”向“研發(fā)驅(qū)動(dòng)”轉(zhuǎn)變，是保證經(jīng)濟(jì)高質(zhì)量發(fā)展同時(shí)實(shí)現(xiàn)碳減排的有效途徑；此外，能源使用結(jié)構(gòu)的改變和升級(jí)是一個(gè)較長(zhǎng)期的過(guò)程，在大力發(fā)展清潔能源的同時(shí)，應(yīng)強(qiáng)調(diào)優(yōu)化能源消費(fèi)結(jié)構(gòu)、促使能源利用技術(shù)創(chuàng)新、提高能源使用效率。最后，應(yīng)充分考慮經(jīng)濟(jì)環(huán)境和企業(yè)生存環(huán)境，針對(duì)行業(yè)特征實(shí)施和制定異質(zhì)性碳減排、碳達(dá)峰目標(biāo)，并給予相應(yīng)的經(jīng)濟(jì)補(bǔ)償，確保發(fā)揮企業(yè)對(duì)經(jīng)濟(jì)運(yùn)行的創(chuàng)新活力。