數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理耦合及影響因素
——基于我國(guó)30個(gè)省級(jí)面板數(shù)據(jù)的實(shí)證檢驗(yàn)

李曉梅，崔 靚

（遼寧工程技術(shù)大學(xué)營(yíng)銷管理學(xué)院，遼寧葫蘆島 125105）

一、引言

數(shù)字經(jīng)濟(jì)是我國(guó)經(jīng)濟(jì)未來的發(fā)展趨勢(shì)，也是促進(jìn)經(jīng)濟(jì)高質(zhì)量發(fā)展必不可缺的一部分。隨著數(shù)字化技術(shù)不斷發(fā)展，我國(guó)物流業(yè)也呈現(xiàn)出一些新變化，一是通過融入大數(shù)據(jù)、云計(jì)算、區(qū)塊鏈等技術(shù)，構(gòu)建現(xiàn)代數(shù)字物流體系；二是在物流行業(yè)中結(jié)合物流大數(shù)據(jù)平臺(tái)、物流云、智能物流等，推動(dòng)我國(guó)生產(chǎn)、流通、消費(fèi)方式的變革。數(shù)字技術(shù)的應(yīng)用為物流業(yè)增添了新動(dòng)力，使得我國(guó)物流業(yè)快速發(fā)展。根據(jù)國(guó)家統(tǒng)計(jì)局最近發(fā)布的數(shù)據(jù)，2020年我國(guó)物流產(chǎn)業(yè)增加值為41 561.7億元，占國(guó)內(nèi)生產(chǎn)總值的4.1%。

《中華人民共和國(guó)國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展第十四個(gè)五年規(guī)劃和2035年遠(yuǎn)景目標(biāo)綱要》（以下簡(jiǎn)稱“十四五”規(guī)劃）提出，要加快數(shù)字化發(fā)展和建設(shè)數(shù)字中國(guó)。根據(jù)2021年《中國(guó)數(shù)字經(jīng)濟(jì)發(fā)展白皮書》數(shù)據(jù)，2020年我國(guó)數(shù)字經(jīng)濟(jì)規(guī)模達(dá)到了39.2 萬億元，數(shù)字經(jīng)濟(jì)占國(guó)內(nèi)生產(chǎn)總值的比重達(dá)38.6%[1]。截至2019年，我國(guó)鐵路運(yùn)營(yíng)里程、公路運(yùn)營(yíng)里程和快遞業(yè)務(wù)量均居世界首位，我國(guó)已經(jīng)成為物流強(qiáng)國(guó)[2]。這些成果的取得離不開國(guó)家政策的支持。比如，我國(guó)在“十四五”規(guī)劃中多次提到“數(shù)字”和“物流”等概念，“數(shù)字”被提及81次，“物流”出現(xiàn)20次；2020年8月，國(guó)有資產(chǎn)監(jiān)督管理委員會(huì)發(fā)布《關(guān)于加快推進(jìn)國(guó)有企業(yè)數(shù)字化轉(zhuǎn)型工作的通知》，要求加快數(shù)字化轉(zhuǎn)型工作，中國(guó)鐵物等央企明確表示，將數(shù)字化轉(zhuǎn)型工作定為重點(diǎn)任務(wù)，爭(zhēng)取實(shí)現(xiàn)數(shù)字鐵物規(guī)劃設(shè)計(jì)目標(biāo)[3]；2019年3月，國(guó)家發(fā)展和改革委員會(huì)發(fā)布《關(guān)于推動(dòng)物流高質(zhì)量發(fā)展促進(jìn)形成強(qiáng)大國(guó)內(nèi)市場(chǎng)的意見》，強(qiáng)調(diào)在今后的發(fā)展中要注重物流的高質(zhì)量發(fā)展[5]；2020年9月，中央財(cái)經(jīng)委員會(huì)提出要注重發(fā)揮流通體系在國(guó)民經(jīng)濟(jì)中的作用，構(gòu)建雙循環(huán)新發(fā)展格局必須注重流通體系的重要性等[5]。數(shù)字物流在我國(guó)經(jīng)濟(jì)發(fā)展中的重要性不言而喻，但在改善生態(tài)環(huán)境方面還有很大利用空間。二氧化碳等溫室氣體排放在很大程度上是由物流的快速發(fā)展造成的[6-7]。例如，物流業(yè)在運(yùn)輸、包裝和儲(chǔ)存過程中可能會(huì)產(chǎn)生二氧化碳和其他污染物?；茉聪M(fèi)和二氧化碳排放量對(duì)經(jīng)濟(jì)可持續(xù)發(fā)展和節(jié)能減排目標(biāo)的實(shí)現(xiàn)有重大影響[8]。根據(jù)國(guó)家郵政局2020年發(fā)布的數(shù)據(jù)，全國(guó)快遞業(yè)務(wù)完成量達(dá)833.6 億件，同比增長(zhǎng)31.2%[9]。龐大的業(yè)務(wù)量使快遞垃圾數(shù)量激增，可回收快遞包裝物的回收利用率不足20%，給我國(guó)的環(huán)境帶來了很大的壓力，影響環(huán)境的可持續(xù)發(fā)展。數(shù)字物流在促進(jìn)經(jīng)濟(jì)發(fā)展的同時(shí)，能夠解決一部分碳環(huán)境治理的問題，因此數(shù)字物流如何在促進(jìn)物流業(yè)發(fā)展的前提下實(shí)現(xiàn)碳環(huán)境治理和經(jīng)濟(jì)的可持續(xù)發(fā)展是一個(gè)值得研究的問題。

以往的研究主要關(guān)注經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境的關(guān)系，其中最著名的是環(huán)境庫茲涅茨曲線（Envi?ronmental Kuznets Curve，EKC）。此外，相當(dāng)多的學(xué)者關(guān)注經(jīng)濟(jì)發(fā)展與物流發(fā)展的關(guān)系以及物流發(fā)展與生態(tài)環(huán)境的關(guān)系。然而，從耦合協(xié)調(diào)角度研究數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理關(guān)系的卻很少。麥基利斯（Maciulis）等[10]通過構(gòu)建運(yùn)輸服務(wù)水平評(píng)價(jià)指標(biāo)來分析運(yùn)輸服務(wù)業(yè)對(duì)經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境的影響，認(rèn)為運(yùn)輸系統(tǒng)的改善可以促進(jìn)經(jīng)濟(jì)發(fā)展和改善生態(tài)環(huán)境。本文以我國(guó)30 個(gè)?。ㄗ灾螀^(qū)、直轄市）為研究對(duì)象，根據(jù)2013—2018年的統(tǒng)計(jì)數(shù)據(jù)，分析數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理的耦合協(xié)調(diào)度，通過探索性空間分析來測(cè)度其分布是否與空間相關(guān)，并利用回歸分析來探究耦合協(xié)調(diào)度影響因素的重要性，并對(duì)結(jié)論進(jìn)行穩(wěn)健性檢驗(yàn)。

二、研究綜述

（一）數(shù)字物流與區(qū)域經(jīng)濟(jì)的關(guān)系

數(shù)字物流由智慧物流發(fā)展而來，國(guó)際商業(yè)機(jī)器公司（International Business Machine，IBM）最早將智慧物流定義為：通過傳感器、射頻識(shí)別（RFID）、全球定位系統(tǒng)（GPS）等技術(shù)來生成實(shí)時(shí)信息進(jìn)行管理的活動(dòng)。隨著數(shù)字技術(shù)應(yīng)用的系統(tǒng)化，張則強(qiáng)等[11]首次提出了數(shù)字物流的概念，認(rèn)為數(shù)字物流是在互聯(lián)網(wǎng)、數(shù)據(jù)庫等技術(shù)的支持下，應(yīng)用數(shù)字技術(shù)對(duì)信息進(jìn)行處理的技術(shù)系統(tǒng)。數(shù)字物流的核心是用數(shù)字化技術(shù)為物流行業(yè)賦能，與此同時(shí)還需要先進(jìn)的物流管理技術(shù)、數(shù)據(jù)孿生和超庫存等技術(shù)[12]。在物流與區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展的關(guān)系方面，卡伊奇（Kayikci）[13]提出，物流產(chǎn)業(yè)集群會(huì)促進(jìn)經(jīng)濟(jì)的可持續(xù)性發(fā)展；楊守德[14]用經(jīng)濟(jì)總量變動(dòng)和物流經(jīng)營(yíng)績(jī)效來衡量技術(shù)創(chuàng)新對(duì)物流的影響，發(fā)現(xiàn)降低單位物流費(fèi)用能夠推動(dòng)物流業(yè)高質(zhì)量發(fā)展，最終促進(jìn)經(jīng)濟(jì)的高質(zhì)量發(fā)展；陳瑩[15]從經(jīng)濟(jì)高質(zhì)量發(fā)展視角出發(fā)，以山東省為研究對(duì)象，對(duì)其物流能力進(jìn)行評(píng)價(jià)，并提出相應(yīng)的建議；上官緒明[16]以我國(guó)30個(gè)省市區(qū)為樣本研究發(fā)現(xiàn)，物流業(yè)集聚對(duì)制造業(yè)高質(zhì)量發(fā)展具有門檻效應(yīng)。綜上所述，物流的發(fā)展可以有效地促進(jìn)區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展。本文在物流的研究基礎(chǔ)上融入數(shù)字化發(fā)展水平，旨在研究數(shù)字物流與區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的耦合關(guān)系。

（二）碳環(huán)境治理內(nèi)涵及其與區(qū)域經(jīng)濟(jì)發(fā)展的關(guān)系

2003年英國(guó)發(fā)布的能源白皮書《我們能源的未來：創(chuàng)建低碳經(jīng)濟(jì)》最早提出“低碳經(jīng)濟(jì)”一詞，但沒有界定其概念。2006年，經(jīng)濟(jì)學(xué)家斯特恩（Stern）發(fā)表了《斯特恩報(bào)告》，這份報(bào)告認(rèn)為低碳經(jīng)濟(jì)是低碳產(chǎn)業(yè)和低碳技術(shù)相結(jié)合的一種新經(jīng)濟(jì)形態(tài)。環(huán)境治理是多方合作的結(jié)果，不僅要提高中央的政策置信度和強(qiáng)化地方的問責(zé)制，還要鼓勵(lì)社會(huì)和公眾參與，以此提高環(huán)境污染的治理效率[17]。綜上，本文認(rèn)為，碳環(huán)境治理需要政府、社會(huì)和公眾通過合作的方式來做。

由于有關(guān)區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的研究較少，所以本文從與碳環(huán)境治理相關(guān)的三個(gè)方面來進(jìn)行闡述。一是區(qū)域經(jīng)濟(jì)和低碳的關(guān)系。鄔彩霞[18]從能源流、資源流兩方面重新定義低碳經(jīng)濟(jì)，并運(yùn)用2005—2019年我國(guó)低碳治理和區(qū)域經(jīng)濟(jì)等相關(guān)數(shù)據(jù)進(jìn)行實(shí)證分析，研究表明，低碳經(jīng)濟(jì)能夠促進(jìn)區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展。實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)、能源消費(fèi)和低碳之間的協(xié)調(diào)發(fā)展，能夠更好地促進(jìn)經(jīng)濟(jì)高質(zhì)量發(fā)展[19]。二是區(qū)域經(jīng)濟(jì)和碳減排的關(guān)系。孟（Meng）等[20]發(fā)現(xiàn)，碳排放能夠促進(jìn)經(jīng)濟(jì)的發(fā)展，且存在空間溢出效應(yīng)。要降低單位GDP 的二氧化碳強(qiáng)度，提高二氧化碳的產(chǎn)出效益，最終達(dá)到區(qū)域經(jīng)濟(jì)增長(zhǎng)和二氧化碳減排雙贏的目標(biāo)[21]。三是區(qū)域經(jīng)濟(jì)和生態(tài)環(huán)境的關(guān)系。1993年庫茲涅茨提出的環(huán)境庫茲涅茨假說是研究區(qū)域經(jīng)濟(jì)和生態(tài)環(huán)境的重要理論。隋建利等[22]使用非線性馬爾科夫區(qū)制轉(zhuǎn)移因果模型測(cè)度了中國(guó)經(jīng)濟(jì)增長(zhǎng)與生態(tài)環(huán)境的非線性關(guān)系，研究表明，中國(guó)的環(huán)境惡化是由經(jīng)濟(jì)的快速增長(zhǎng)造成的。因此，在做好低碳、碳減排工作的同時(shí)保護(hù)好生態(tài)環(huán)境，就能在一定程度上促進(jìn)區(qū)域經(jīng)濟(jì)的發(fā)展。

（三）物流產(chǎn)業(yè)與碳環(huán)境治理的關(guān)系

關(guān)于物流產(chǎn)業(yè)與碳效率、碳減排及生態(tài)環(huán)境的研究受到學(xué)界廣泛關(guān)注，這為物流產(chǎn)業(yè)與碳環(huán)境治理的相關(guān)研究提供了理論依據(jù)。本文將相關(guān)研究文獻(xiàn)分為三類：一是對(duì)物流與碳效率之間關(guān)系的研究。姚山季等[23]運(yùn)用三階段數(shù)據(jù)包絡(luò)分析（Data Envelopment Analysis，DEA）和曼奎斯特（Malmquist）指數(shù)對(duì)我國(guó)17 個(gè)省（自治區(qū)、直轄市）低碳物流效率進(jìn)行測(cè)度，研究表明，外部環(huán)境和隨機(jī)因素都會(huì)影響低碳物流效率。二是對(duì)物流與碳減排關(guān)系的研究。烏韋達(dá)（Ubeda）等[25]以西班牙的一家配送公司為例進(jìn)行分析，采用有約束的車輛路徑規(guī)劃（Capacitated Vehicle Routing Problem，CVRP）來達(dá)到運(yùn)輸距離和碳排放量最優(yōu)的狀態(tài)。臧新等[25]對(duì)外商直接投資與中國(guó)物流業(yè)碳排放的關(guān)系進(jìn)行了實(shí)證分析，研究表明，外商直接投資能夠顯著地減少和降低中國(guó)物流業(yè)CO2排放量和排放強(qiáng)度。三是對(duì)物流與生態(tài)環(huán)境關(guān)系的研究。劉（Liu）等[26]收集了相關(guān)的物流和環(huán)境數(shù)據(jù)，通過對(duì)42個(gè)亞洲國(guó)家進(jìn)行GMM回歸模型分析，發(fā)現(xiàn)物流運(yùn)輸過程中產(chǎn)生的二氧化碳和物流基礎(chǔ)設(shè)施建造對(duì)生態(tài)環(huán)境造成了不利影響。周泰[27]構(gòu)建了區(qū)域物流和生態(tài)環(huán)境的協(xié)調(diào)度模型，研究表明，四川省的物流和生態(tài)環(huán)境處于穩(wěn)定發(fā)展?fàn)顟B(tài)。綜上所述，要在發(fā)展好物流產(chǎn)業(yè)的同時(shí)，做好碳減排工作和生態(tài)環(huán)境保護(hù)工作，使物流產(chǎn)業(yè)和碳環(huán)境治理達(dá)到協(xié)調(diào)發(fā)展的水平。

本文研究的主要問題為：我國(guó)30 個(gè)?。ㄗ灾螀^(qū)、直轄市）數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的耦合協(xié)調(diào)水平處于哪個(gè)階段；各地區(qū)的耦合協(xié)調(diào)度是否與空間分布相關(guān)；人均可支配收入、貨運(yùn)量、互聯(lián)網(wǎng)普及率、單位GDP二氧化碳排放量和生態(tài)建設(shè)與保護(hù)投資對(duì)耦合協(xié)調(diào)度的影響程度。

本研究的貢獻(xiàn)：將數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理當(dāng)做一個(gè)整體來研究其協(xié)調(diào)水平，繼承和拓展區(qū)域均衡理論和可持續(xù)發(fā)展理論；引入數(shù)字物流研究視角，為構(gòu)建現(xiàn)代流通體系奠定了理論基礎(chǔ)，進(jìn)而促進(jìn)雙循環(huán)格局的形成；在傳統(tǒng)生態(tài)環(huán)境研究的基礎(chǔ)上提出了碳環(huán)境治理這一新研究視角，對(duì)低碳經(jīng)濟(jì)理論起到了補(bǔ)充作用，為實(shí)現(xiàn)碳達(dá)峰、碳中和目標(biāo)提供了可借鑒的分析框架；為提升各省（自治區(qū)、直轄市）協(xié)調(diào)水平提供了新路徑，剖析了各?。ㄗ灾螀^(qū)、直轄市）應(yīng)采取的做法，對(duì)提升各?。ㄗ灾螀^(qū)、直轄市）整體協(xié)調(diào)水平有著重要意義。

三、指標(biāo)選取與模型構(gòu)建

（一）指標(biāo)選取

考慮到指標(biāo)體系的系統(tǒng)性、完整性、科學(xué)性和數(shù)據(jù)的可用性，結(jié)合相關(guān)研究成果，建立數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的指標(biāo)體系。

由于數(shù)字物流方面的研究較少，所以本文關(guān)于數(shù)字物流的指標(biāo)選取借鑒了評(píng)價(jià)物流和互聯(lián)網(wǎng)發(fā)展的指標(biāo)。如蘭（Lan）等[28]和崔宏凱等[29]認(rèn)為，應(yīng)該用物流業(yè)固定資產(chǎn)投資、物流產(chǎn)業(yè)從業(yè)人員數(shù)量來衡量物流的投入能力。張毅等[30]、孫浩杰等[31]和蘭等[28]認(rèn)為，應(yīng)用貨運(yùn)周轉(zhuǎn)量、貨運(yùn)量、物流業(yè)總產(chǎn)值來衡量物流的產(chǎn)出能力。李立威等[32]、希梅內(nèi)斯（Jiménez）等[33]、劉姿均等[34]、韓先鋒等[35]、劉軍等[36]和王娟娟等[37]用互聯(lián)網(wǎng)普及率、電話普及率、光纜線路長(zhǎng)度、電子商務(wù)銷售額和電子商務(wù)采購額來測(cè)度數(shù)字化發(fā)展水平。

在選取區(qū)域經(jīng)濟(jì)的指標(biāo)時(shí)，借鑒了“十四五”時(shí)期經(jīng)濟(jì)社會(huì)發(fā)展指標(biāo)體系的指標(biāo)。唐建榮等[38]、祝捷等[39]、曼蘇爾（Manzoor）等[40]、崔宏凱等[29]、王（Wang）等[41]認(rèn)為，可用人均可支配收入、社會(huì)消費(fèi)品零售總額、地方財(cái)政收入來衡量經(jīng)濟(jì)增長(zhǎng)水平。曼蘇爾等[40]、曹翠珍等[42]、趙曉敏等[43]、張（Zhang）等[44]用居民消費(fèi)價(jià)格指數(shù)和人均GDP來衡量生活消費(fèi)水平。張毅等[30]、祝捷等[39]和王等[41]用進(jìn)出口總額、出口商品總額來衡量對(duì)外開放水平。

在構(gòu)建碳環(huán)境治理指標(biāo)時(shí)，張等[44]和李虹等[45]用生態(tài)建設(shè)與保護(hù)本年完成投資、工業(yè)污染治理投資來衡量環(huán)境污染治理，并用人均水資源量和森林覆蓋率來衡量資源環(huán)境水平。王國(guó)惠等[46]用城市綠地面積來衡量資源環(huán)境治理水平。本文借鑒相關(guān)研究構(gòu)建碳環(huán)境治理的指標(biāo)。

本研究構(gòu)建的耦合協(xié)調(diào)評(píng)價(jià)指標(biāo)體系如表1所示。

表1 耦合協(xié)調(diào)評(píng)價(jià)指標(biāo)體系

（二）模型構(gòu)建

1.面板數(shù)據(jù)熵權(quán)法

首先利用面板數(shù)據(jù)熵權(quán)法確定指標(biāo)權(quán)重，然后再對(duì)其進(jìn)行面板數(shù)據(jù)熵權(quán)法處理。

（1）指標(biāo)標(biāo)準(zhǔn)化

其中，xmax和xmin分別表示第j項(xiàng)指標(biāo)的最大值和最小值，xαij和Zαij分別表示指標(biāo)標(biāo)準(zhǔn)化處理前和標(biāo)準(zhǔn)化處理后的值。

（2）指標(biāo)歸一化

其中，m代表年份，k代表?。ㄗ灾螀^(qū)、直轄市），Zαij為標(biāo)準(zhǔn)化處理后的值，為第j 項(xiàng)指標(biāo)所有省（自治區(qū)、直轄市）及年份的觀測(cè)數(shù)據(jù)之和，Pαij為指標(biāo)歸一化后的值。

（3）熵值

（4）各項(xiàng)指標(biāo)的差異系數(shù)

（5）各項(xiàng)指標(biāo)的權(quán)重

2.數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理耦合協(xié)調(diào)模型

耦合的概念源于物理學(xué)，它描述了兩個(gè)或多個(gè)相關(guān)系統(tǒng)之間的相互作用以及作用的程度。耦合協(xié)調(diào)模型應(yīng)用廣泛。陸（Lu）等[47]以武漢市為例，建立了經(jīng)濟(jì)、資源、環(huán)境三個(gè)子系統(tǒng)之間的耦合模型，并對(duì)其協(xié)調(diào)度對(duì)城市可持續(xù)發(fā)展的重要性進(jìn)行了評(píng)價(jià)。李林紅等[48]采用耦合協(xié)調(diào)模型對(duì)云南省物流系統(tǒng)和綠色旅游系統(tǒng)進(jìn)行了實(shí)證分析，研究表明，雖然耦合協(xié)調(diào)度逐年上升，但總體為區(qū)域物流發(fā)展滯后型。由以往研究可以發(fā)現(xiàn)，耦合協(xié)調(diào)模型常用在兩者及以上的關(guān)系分析中。

為了從整體上反映區(qū)域數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理的發(fā)展水平和協(xié)調(diào)效應(yīng)，本文收集了我國(guó)30個(gè)?。ㄗ灾螀^(qū)、直轄市）的數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理相關(guān)數(shù)據(jù)，用耦合協(xié)調(diào)模型計(jì)算數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理的耦合協(xié)調(diào)度。

（1）有序度函數(shù)

設(shè)數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理演化發(fā)展過程中的序參量為ui（i=1，2，3），其中uij（j=1，2，3，...，n)為序參量ui的第j 個(gè)指標(biāo)，其值為Xij，序參量的有序度U(uij)可以表示為：

其中，αij為序參數(shù)的上限值，βij為序參數(shù)的下限值。U(uij)∈[ ]0，1 。序參量的有序度越高，對(duì)系統(tǒng)有序度的貢獻(xiàn)越大。

對(duì)所有指標(biāo)進(jìn)行加權(quán)求和：

其中，U(ui)為子系統(tǒng)的有序度，U(ui)∈[0 ，1] ，λij為各指標(biāo)權(quán)重，

（2）耦合度函數(shù)

其中，U(u1)、U(u2)、U(u3)分別表示數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理的有序度，耦合度C ∈[0 ，1] ，D 為耦合協(xié)調(diào)度，T為綜合協(xié)調(diào)指數(shù)，反映系統(tǒng)整體層面上的協(xié)同程度，a、b、c為系數(shù)。借鑒姜磊等[49]的系數(shù)選取方法，本文選取a=1/3，b=1/3，c=1/3。

3.耦合協(xié)調(diào)度的探索性空間數(shù)據(jù)分析

翁鋼民等[50]運(yùn)用探索性空間數(shù)據(jù)分析（Ex?ploratory Spatial Data Analysis，ESDA）方法來研究我國(guó)旅游和文化產(chǎn)業(yè)的協(xié)調(diào)度，研究表明，協(xié)調(diào)水平在空間上呈正向集聚。

（1）全局空間自相關(guān)

一般使用全局莫蘭指數(shù)（Global Moran's I）檢驗(yàn)相鄰空間分布對(duì)象之間的自相關(guān)性。其值分布在-1～1 之間，當(dāng)I>0 時(shí)，說明屬性值在研究區(qū)域內(nèi)呈現(xiàn)空間集聚現(xiàn)象，I 值越大，空間相關(guān)性越強(qiáng)；當(dāng)I=0時(shí)，呈隨機(jī)性分布；當(dāng)I<0時(shí)，表明屬性值內(nèi)存在差異，I值越小，差異性越大。其計(jì)算公式為：

其中，n是評(píng)價(jià)個(gè)數(shù)，xi和xj為i和j的屬性值，xˉ為平均值，Wij為空間權(quán)重矩陣。另外，用標(biāo)準(zhǔn)統(tǒng)計(jì)量Z值來檢驗(yàn)Moran's I的顯著性。

（2）局部空間自相關(guān)

全局空間相關(guān)性容易忽視內(nèi)部的不平穩(wěn)性，因此通過計(jì)算局部莫蘭指數(shù)（Local Moran'I）來反映各單元的異質(zhì)性，其計(jì)算公式為：

其中，Wij表示空間權(quán)重矩陣，xˉ為鄰近區(qū)域?qū)傩灾灯畹木怠?/p>

4.計(jì)量模型

為探尋數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理耦合協(xié)調(diào)水平的影響因素，以耦合協(xié)調(diào)度為被解釋變量，參考相關(guān)研究，選取數(shù)字物流發(fā)展水平、區(qū)域經(jīng)濟(jì)發(fā)展水平和碳環(huán)境治理發(fā)展水平三個(gè)方面的指標(biāo)作為解釋變量，并按重要性進(jìn)行排序。指標(biāo)選取情況如表2所示。

表2 耦合協(xié)調(diào)度及其影響因素指標(biāo)選取結(jié)果

基于所選取的指標(biāo)進(jìn)行實(shí)證分析，對(duì)變量進(jìn)行標(biāo)準(zhǔn)化處理后，選擇固定效應(yīng)模型進(jìn)行回歸分析，公式如下：

其中，CDit為耦合協(xié)調(diào)度，α 為常數(shù)項(xiàng)，β1、β2、β3、β4和β5為待定系數(shù)，ui為個(gè)體固定效應(yīng)，ε為隨機(jī)誤差項(xiàng)。

四、實(shí)證分析

（一）數(shù)據(jù)來源及處理

本文以物流業(yè)固定資產(chǎn)投資為例對(duì)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，利用公式（1）對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，得到的結(jié)果如表3所示。

表3 2013—2018年各省份物流固定資產(chǎn)投資標(biāo)準(zhǔn)化處理結(jié)果

（二）數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理的耦合協(xié)調(diào)度分析

采用2013—2018年我國(guó)30 個(gè)省（自治區(qū)、直轄市）的數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理面板數(shù)據(jù)。利用公式（8）、（9）、（10）計(jì)算出各?。ㄗ灾螀^(qū)、直轄市）數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的耦合協(xié)調(diào)度。計(jì)算結(jié)果如表4所示。

根據(jù)高楠等[51]對(duì)耦合協(xié)調(diào)度的等級(jí)分類標(biāo)準(zhǔn)，以2018年為例，可得到各?。ㄗ灾螀^(qū)、直轄市）的協(xié)調(diào)等級(jí)，如表5所示。

分析表4和表5數(shù)據(jù)可以得到如下結(jié)果：

表4 2013—2018年各省份耦合協(xié)調(diào)度

表5 2018年各協(xié)調(diào)等級(jí)下的省份

2018年30 個(gè)?。ㄗ灾螀^(qū)、直轄市）的耦合協(xié)調(diào)度均值在0.49 左右，總體上數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理介于瀕臨失調(diào)階段和勉強(qiáng)協(xié)調(diào)階段之間，并且其波動(dòng)幅度不大，這反映出耦合協(xié)調(diào)度格局具有一定的穩(wěn)定性。

2013—2018年，我國(guó)數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理整體發(fā)展水平逐年提高，但地區(qū)之間的發(fā)展存在差異。下面從每個(gè)等級(jí)中選取有代表性的地區(qū)進(jìn)行分析。在中度失調(diào)等級(jí)的?。ㄗ灾螀^(qū)、直轄市）中，寧夏回族自治區(qū)由于地處內(nèi)陸，數(shù)字基礎(chǔ)設(shè)施建設(shè)不足，阻礙了該地區(qū)數(shù)字物流的發(fā)展，且互聯(lián)網(wǎng)普及率、人均可支配收入和電子商務(wù)銷售額等指標(biāo)數(shù)值不高，致使該地區(qū)的數(shù)字物流和區(qū)域經(jīng)濟(jì)有序度較低；在輕度失調(diào)等級(jí)的省（自治區(qū)、直轄市）中，天津市由于物流業(yè)固定資產(chǎn)投資較少，致使其物流有序度較低，在30 個(gè)省（自治區(qū)、直轄市）中排在第27 位（見表3），與其他?。ㄗ灾螀^(qū)、直轄市）相比，天津市的森林覆蓋率較低，人均水資源量少，致使其碳環(huán)境治理有序度不高；在瀕臨失調(diào)等級(jí)的?。ㄗ灾螀^(qū)、直轄市）中，位于我國(guó)東北地區(qū)的黑龍江省近年來人才流失嚴(yán)重且資源過度開發(fā)，造成了嚴(yán)重的環(huán)境污染，與其他?。ㄗ灾螀^(qū)、直轄市）相比，東北地區(qū)的數(shù)字物流和區(qū)域經(jīng)濟(jì)較不發(fā)達(dá)；在勉強(qiáng)協(xié)調(diào)等級(jí)的七個(gè)?。ㄗ灾螀^(qū)、直轄市）中，上海市由于森林覆蓋率較低，生態(tài)建設(shè)與保護(hù)投資不足，導(dǎo)致其碳環(huán)境治理有序度較低；在初級(jí)協(xié)調(diào)等級(jí)的四個(gè)?。ㄗ灾螀^(qū)、直轄市）中，浙江省屬于經(jīng)濟(jì)發(fā)達(dá)地區(qū)，2018年6月推出的運(yùn)輸結(jié)構(gòu)調(diào)整政策促進(jìn)了物流行業(yè)貨運(yùn)量和貨物周轉(zhuǎn)量的增加，同時(shí)浙江省對(duì)二氧化碳排放量持續(xù)進(jìn)行監(jiān)測(cè)和控制，截止到2018年，浙江省關(guān)于排污的政策文件已經(jīng)出臺(tái)了181 份；廣東省處于良好協(xié)調(diào)等級(jí)，其原因可能是，該省地處東部沿海地區(qū)，不僅有全國(guó)最大的港口群，而且是我國(guó)進(jìn)出口貿(mào)易的重要基地，該省的數(shù)字基礎(chǔ)設(shè)施建設(shè)已經(jīng)相當(dāng)成熟，另外現(xiàn)代物流研究院對(duì)該省物流發(fā)展進(jìn)行的整體指導(dǎo)也是廣東省處于良好協(xié)調(diào)等級(jí)的原

因，但2018年廣東省物流業(yè)單位GDP 二氧化碳排放量在30 個(gè)?。ㄗ灾螀^(qū)、直轄市）中排在第11 位，因此其控制二氧化碳的工作還需要進(jìn)一步落實(shí)。

（三）數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理耦合協(xié)調(diào)度的探索性空間數(shù)據(jù)分析

1.耦合協(xié)調(diào)度的全局空間自相關(guān)分析

利用GEODA 軟件計(jì)算30 個(gè)?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度全局莫蘭指數(shù)I（表6），如果其通過5%的顯著性檢驗(yàn)，則表明耦合協(xié)調(diào)度在空間上有集聚性。

表6 耦合協(xié)調(diào)度的全局莫蘭指數(shù)

數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理耦合協(xié)調(diào)度的莫蘭指數(shù)2013年為0.258，2018年為0.265，均為正值，說明耦合協(xié)調(diào)度的空間聚集性呈現(xiàn)出一定的正相關(guān)性，且隨著時(shí)間的推移，耦合協(xié)調(diào)度的空間正相關(guān)性不斷增強(qiáng)。

2.耦合協(xié)調(diào)度LISA分析

高高集聚區(qū)代表本?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度高，周邊?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度也高。由表7可知，2018年，福建省為高高集聚區(qū)，這一方面說明其數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理耦合協(xié)調(diào)度較高，周邊?。ㄗ灾螀^(qū)、直轄市）的耦合協(xié)調(diào)度也較高，且具有較強(qiáng)的相關(guān)性，地區(qū)之間起到了協(xié)同作用；另一方面說明該省近幾年經(jīng)濟(jì)發(fā)展水平處于上升趨勢(shì)，相鄰省（自治區(qū)、直轄市）也更加注重三者的協(xié)調(diào)發(fā)展，從而帶動(dòng)了鄰?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度的提高。福建省處于沿海地區(qū)，地理位置優(yōu)越，地區(qū)較為開放，能夠跟周邊?。ㄗ灾螀^(qū)、直轄市）進(jìn)行更深入的交流，促進(jìn)了鄰省（自治區(qū)、直轄市）耦合協(xié)調(diào)度的提高。

表7 耦合協(xié)調(diào)度的LISA分析

低高集聚區(qū)代表本?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度低，但周邊省（自治區(qū)、直轄市）耦合協(xié)調(diào)度高。江西省為低高集聚區(qū)，且與高高集聚區(qū)相鄰，該省鄰接耦合協(xié)調(diào)度相對(duì)較高的?。ㄗ灾螀^(qū)、直轄市），其數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理耦合協(xié)調(diào)度增長(zhǎng)較慢，受耦合協(xié)調(diào)度增長(zhǎng)較快區(qū)域的影響較小。江西省耦合協(xié)調(diào)度較低，而鄰省（自治區(qū)、直轄市）耦合協(xié)調(diào)度較高，這可能與人口流失有關(guān)，大部分人口可能前往長(zhǎng)三角或者珠三角等發(fā)展較好的地區(qū)尋找工作。

低低集聚區(qū)代表本?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度低，周邊?。ㄗ灾螀^(qū)、直轄市）耦合協(xié)調(diào)度也低。新疆、甘肅、內(nèi)蒙古為低低集聚區(qū)，這些地區(qū)的數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理都欠佳，出現(xiàn)了低耦合協(xié)調(diào)度和低耦合協(xié)調(diào)度的集聚現(xiàn)象。但隨著經(jīng)濟(jì)的發(fā)展、國(guó)家政策的支持和科學(xué)技術(shù)水平的提高，這些地區(qū)數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理的耦合協(xié)調(diào)度會(huì)不斷提高。

（四）計(jì)量模型

1.描述性統(tǒng)計(jì)

描述性統(tǒng)計(jì)的具體結(jié)果如表8所示。

表8 變量描述性統(tǒng)計(jì)結(jié)果

2.相關(guān)系數(shù)矩陣

為了檢驗(yàn)核心解釋變量之間是否存在多重共線性，通過觀察相關(guān)系數(shù)矩陣（表9）和計(jì)算方差膨脹因子（表10）進(jìn)行進(jìn)一步確定。結(jié)果顯示，VIF<5，這表明核心解釋變量之間不存在多重共線性問題。

表9 解釋變量相關(guān)系數(shù)矩陣

表10 解釋變量的方差膨脹因子

3.固定效應(yīng)回歸模型分析

經(jīng)過豪斯曼（Hausman）檢驗(yàn)，最終選擇固定效應(yīng)模型對(duì)數(shù)據(jù)進(jìn)行回歸分析。結(jié)果如表11所示。

通過表11 數(shù)據(jù)可以看出，人均可支配收入對(duì)耦合協(xié)調(diào)度的影響最大，貨運(yùn)量次之，互聯(lián)網(wǎng)普及率排在第三位，生態(tài)建設(shè)與保護(hù)投資對(duì)耦合協(xié)調(diào)度的影響最小。具體而言：人均可支配收入與耦合協(xié)調(diào)度呈正相關(guān)，經(jīng)濟(jì)發(fā)展水平的提高也會(huì)促使耦合協(xié)調(diào)度提高；貨運(yùn)量的相關(guān)系數(shù)為0.221，且在1%水平下顯著，表明貨運(yùn)量與耦合協(xié)調(diào)度顯著正相關(guān)，因?yàn)樨涍\(yùn)量在很大程度上體現(xiàn)了物流業(yè)發(fā)展水平，物流業(yè)在飛速發(fā)展的同時(shí)，也在促進(jìn)經(jīng)濟(jì)的發(fā)展，因此會(huì)增加耦合協(xié)調(diào)度；互聯(lián)網(wǎng)普及率與耦合協(xié)調(diào)度呈正相關(guān)，因?yàn)榛ヂ?lián)網(wǎng)普及率代表著企業(yè)的數(shù)字化基礎(chǔ)設(shè)施水平，互聯(lián)網(wǎng)普及率的提高會(huì)促進(jìn)物流業(yè)的發(fā)展，最終達(dá)到提高耦合協(xié)調(diào)度的目的；生態(tài)建設(shè)與保護(hù)投資的相關(guān)系數(shù)為0.098，通過了顯著性檢驗(yàn)，對(duì)耦合協(xié)調(diào)度有一定的影響；單位GDP 二氧化碳排放量沒有通過顯著性檢驗(yàn)，可能是樣本容量不夠，導(dǎo)致其回歸結(jié)果不顯著。

表11 固定效應(yīng)模型回歸分析結(jié)果

4.穩(wěn)健性檢驗(yàn)

接下來通過穩(wěn)健性檢驗(yàn)進(jìn)一步檢驗(yàn)結(jié)論的可靠性。一般采用替換變量和改變時(shí)間段回歸的方法進(jìn)行檢驗(yàn)，但考慮到數(shù)據(jù)缺失等問題，本文借鑒孫黃平等[52]的方法，將人均可支配收入替換為人均GDP來進(jìn)行穩(wěn)健性檢驗(yàn)，分析結(jié)果如表12所示。

表12 計(jì)量模型回歸分析（穩(wěn)健性檢驗(yàn)）

檢驗(yàn)結(jié)果表明，結(jié)論與未替換變量之前是一致的，由此推斷本文所得結(jié)論是十分穩(wěn)健的。

五、研究結(jié)論與政策啟示

（一）研究結(jié)論

本文以2013—2018年30 個(gè)?。ㄗ灾螀^(qū)、直轄市）為樣本，利用耦合協(xié)調(diào)模型對(duì)數(shù)字物流、區(qū)域經(jīng)濟(jì)、碳環(huán)境治理之間的耦合協(xié)調(diào)度進(jìn)行分析，以耦合協(xié)調(diào)度為被解釋變量來進(jìn)行空間相關(guān)分析，并從區(qū)域經(jīng)濟(jì)發(fā)展水平、數(shù)字物流發(fā)展水平、碳環(huán)境治理水平三個(gè)層面探究三者協(xié)調(diào)發(fā)展的影響因素，對(duì)其重要性進(jìn)行排序，得出：

第一，各省（自治區(qū)、直轄市）耦合協(xié)調(diào)度多集中在0.2～0.9之間，但是地區(qū)的發(fā)展存在差異。

第二，耦合協(xié)調(diào)度在空間上呈正相關(guān)，在空間上表現(xiàn)出一定的集聚性，空間集聚呈加強(qiáng)狀態(tài)，并對(duì)鄰接省（自治區(qū)、直轄市）產(chǎn)生了影響。

第三，數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理三者回歸結(jié)果表明，各影響因素的重要性排序?yàn)椋喝司芍涫杖?貨運(yùn)量>互聯(lián)網(wǎng)普及率>生態(tài)建設(shè)與保護(hù)投資。

（二）政策建議

第一，從地區(qū)之間耦合協(xié)調(diào)度差異來看，各?。ㄗ灾螀^(qū)、直轄市）應(yīng)統(tǒng)籌規(guī)劃數(shù)字物流、區(qū)域經(jīng)濟(jì)與碳環(huán)境治理之間的關(guān)系，不能為了促進(jìn)數(shù)字物流和區(qū)域經(jīng)濟(jì)的高質(zhì)量發(fā)展，而犧牲碳環(huán)境治理。比如，遼寧省應(yīng)該加快區(qū)域物流公共信息平臺(tái)建設(shè)，發(fā)揮平臺(tái)整合資源的能力，真正實(shí)現(xiàn)物流業(yè)的資源共享和信息共用，在發(fā)展數(shù)字物流的同時(shí)，還要控制二氧化碳排放量，可以定期對(duì)城市空氣質(zhì)量進(jìn)行監(jiān)測(cè)。

第二，從耦合協(xié)調(diào)度在空間上的集聚性來看，處于低低集聚區(qū)的地區(qū)應(yīng)該借鑒東部沿海地區(qū)的先進(jìn)經(jīng)驗(yàn)和發(fā)展模式，多與耦合協(xié)調(diào)度較高的?。ㄗ灾螀^(qū)、直轄市）交流學(xué)習(xí)，打造自己獨(dú)有的數(shù)字物流發(fā)展模式。如廣東省發(fā)展較好的地區(qū)應(yīng)該在繼續(xù)維持自身發(fā)展勢(shì)頭的基礎(chǔ)上，積極向外輸送人才、技術(shù)和經(jīng)驗(yàn)，打破地域限制，帶動(dòng)全國(guó)數(shù)字物流、區(qū)域經(jīng)濟(jì)和碳環(huán)境治理的優(yōu)質(zhì)協(xié)調(diào)發(fā)展。

第三，從提高耦合協(xié)調(diào)度的方式來看，各?。ㄗ灾螀^(qū)、直轄市）應(yīng)找準(zhǔn)自身的特色產(chǎn)業(yè)，打造核心競(jìng)爭(zhēng)產(chǎn)業(yè)，提升人民的生活水平。在數(shù)字物流方面，應(yīng)該重視各?。ㄗ灾螀^(qū)、直轄市）的物流貨運(yùn)量和互聯(lián)網(wǎng)普及率，一方面運(yùn)用大數(shù)據(jù)技術(shù)使數(shù)字物流數(shù)據(jù)庫與車輛聯(lián)網(wǎng)數(shù)據(jù)庫實(shí)現(xiàn)共享，優(yōu)化運(yùn)輸結(jié)構(gòu)，提高各?。ㄗ灾螀^(qū)、直轄市）貨運(yùn)量；另一方面在國(guó)家加大對(duì)西部地區(qū)政策扶持力度的基礎(chǔ)上，進(jìn)一步加強(qiáng)數(shù)字基礎(chǔ)設(shè)施的建設(shè)，如長(zhǎng)途光纜線路的鋪設(shè)、互聯(lián)網(wǎng)的入戶工作等，提高互聯(lián)網(wǎng)普及率。在碳環(huán)境治理方面，政府應(yīng)該主動(dòng)發(fā)展節(jié)能減排技術(shù)和可再生能源技術(shù)，以此減少物流業(yè)的二氧化碳排放。