郝愛民，任 禛

（鄭州航空工業(yè)管理學(xué)院經(jīng)濟學(xué)院，河南鄭州 450046）

受國際市場環(huán)境和新冠肺炎疫情的影響，全球產(chǎn)業(yè)鏈體系的脆弱性日益凸顯，我國產(chǎn)業(yè)鏈安全穩(wěn)定也面臨著重大考驗，亟須依托國內(nèi)大循環(huán)推進產(chǎn)業(yè)間協(xié)調(diào)聯(lián)動、優(yōu)化市場資源配置，加快技術(shù)補鏈、市場強鏈。流通業(yè)和制造業(yè)作為產(chǎn)業(yè)鏈分工體系內(nèi)的重要主體，其協(xié)同融合的日漸深化引起社會各界的廣泛關(guān)注，并逐漸上升至國家戰(zhàn)略高度?！丁笆奈濉眹鴥?nèi)貿(mào)易發(fā)展規(guī)劃》明確指出要促進流通與制造協(xié)同聯(lián)動，深化現(xiàn)代流通業(yè)與制造業(yè)融合創(chuàng)新發(fā)展，為構(gòu)建新發(fā)展格局提供有力支撐。流通業(yè)以其促進供需精準(zhǔn)匹配、有效反饋市場信息、引領(lǐng)產(chǎn)能資源配置等先導(dǎo)性作用，推動制造業(yè)績效改善［1］；制造業(yè)規(guī)模持續(xù)擴張、效率日益提升，展現(xiàn)出強大的韌性，又為加速流通基礎(chǔ)設(shè)施建設(shè)、提升流通業(yè)服務(wù)質(zhì)量提供有力保障，尤其是在當(dāng)前“雙循環(huán)”新發(fā)展格局下，制造業(yè)韌性的不斷增強對保障我國產(chǎn)業(yè)鏈供應(yīng)鏈安全、推動產(chǎn)業(yè)高質(zhì)量發(fā)展以及加快制造強國建設(shè)具有重要影響。但是，我國制造業(yè)的發(fā)展仍面臨著價值鏈低端鎖定、高端要素集聚不足、自主創(chuàng)新能力薄弱、區(qū)域發(fā)展不平衡等諸多困境。而伴隨新一代信息技術(shù)的蓬勃發(fā)展，流通業(yè)的數(shù)字化改造能夠有效降低搜尋成本、拓寬供銷渠道、提高運營效率，有助于彌補制造業(yè)產(chǎn)業(yè)鏈供應(yīng)鏈的薄弱環(huán)節(jié)，增強我國制造業(yè)在國際競爭中的韌性和活力。因此，深入探究流通數(shù)字化與制造業(yè)韌性的交互作用機理和空間溢出效應(yīng)，形成流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展并輻射帶動周邊地區(qū)共同進步的整體局面，對于加快產(chǎn)業(yè)間良性互動、區(qū)域間創(chuàng)新聯(lián)動，實現(xiàn)經(jīng)濟高質(zhì)量發(fā)展具有重要意義。

1 文獻綜述

首先是流通業(yè)數(shù)字化轉(zhuǎn)型的相關(guān)研究。俞彤暉等［2］、汪旭暉等［3］認為，數(shù)字技術(shù)在流通領(lǐng)域的廣泛應(yīng)用是提升流通效率、強大國內(nèi)市場的必然趨勢。隨著以“互聯(lián)網(wǎng)+流通”為特征的新一輪流通革命的興起，商貿(mào)流通活動逐漸向數(shù)字化轉(zhuǎn)變［4］。Piotrowicz 等［5］、Grewal 等［6］、Choi 等［7］認為，由線上與線下流通渠道邊界日益模糊衍生而來的全渠道零售是具備可視化、智能化等特點的新型業(yè)態(tài)模式，能夠有效降低流通成本、縮短渠道結(jié)構(gòu)。徐鋒等［8］、丁俊發(fā)［9］、謝莉娟等［10］認為，“互聯(lián)網(wǎng)+流通”作為推進流通業(yè)從粗放到集約的突破口，不僅能使信息資源在產(chǎn)業(yè)鏈中高效傳遞，還能實現(xiàn)產(chǎn)銷關(guān)系由“推”式向“拉”式轉(zhuǎn)變的媒介機制升級。張建軍等［11］、唐任伍等［12］的研究表明，新零售、智慧物流等新模式在傳統(tǒng)流通產(chǎn)業(yè)鏈的重構(gòu)與升級中扮演著重要角色。

其次有關(guān)韌性的測度以及影響制造業(yè)韌性的因素分析。國外對韌性的研究起步較早，如Rose 等［13］、Martin 等［14］、Sensier 等［15］對經(jīng)濟韌性的內(nèi)涵、測度以及發(fā)展模式等展開了較為詳細的論述。此后，Doran 等［16］利用靜態(tài)凡登定律衡量美國都市區(qū)韌性指數(shù)；借鑒該方法，徐圓等［17］對我國城市經(jīng)濟韌性進行測度；李蘭冰等［18］、劉瑞等［19］對我國制造業(yè)韌性進行衡量，從產(chǎn)業(yè)適應(yīng)性調(diào)整視角分析了生產(chǎn)性服務(wù)業(yè)專業(yè)化集聚、制造業(yè)空間集聚對制造業(yè)韌性的影響。此外，Rajesh［20］、Hao 等［21］、Agarwal等［22］學(xué)者通過實證研究發(fā)現(xiàn)，柔性經(jīng)營戰(zhàn)略、精益生產(chǎn)、技術(shù)創(chuàng)新等有助于提升制造業(yè)韌性。

最后是關(guān)于數(shù)字技術(shù)影響制造業(yè)轉(zhuǎn)型升級的相關(guān)探討。部分學(xué)者從產(chǎn)業(yè)結(jié)構(gòu)升級的角度展開研究，如Caputo 等［23］構(gòu)建了理論分析框架，闡述了物聯(lián)網(wǎng)技術(shù)推動制造業(yè)結(jié)構(gòu)升級的重要作用；嚴(yán)北戰(zhàn)等［24］運用省級細分行業(yè)數(shù)據(jù)，從供需雙向驅(qū)動的視角分析了“互聯(lián)網(wǎng)+”對制造業(yè)升級的影響。還有部分學(xué)者從價值鏈升級的角度進行分析，如Marilungo等［25］認為信息通信技術(shù)（ICT）賦能制造業(yè)將帶來可持續(xù)的創(chuàng)新模式，推動制造業(yè)升級、完善基礎(chǔ)設(shè)施建設(shè)，進而為客戶提供更多服務(wù)；劉德學(xué)等［26］基于跨國面板數(shù)據(jù)，從數(shù)量、質(zhì)量雙重視角實證檢驗了ICT 有助于制造業(yè)全球價值鏈攀升；何文彬［27］對制造業(yè)細分行業(yè)的貿(mào)易增加值、垂直專業(yè)化分工等進行多維解析發(fā)現(xiàn)，數(shù)字化投入對制造業(yè)知識密集型部門價值鏈高端化轉(zhuǎn)型具有顯著的正向影響。

總體來看，現(xiàn)有相關(guān)研究已經(jīng)關(guān)注到數(shù)字經(jīng)濟背景下流通業(yè)變革、制造業(yè)升級的重要作用，但缺乏對流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展水平的系統(tǒng)測度與異質(zhì)性分析，以及對兩者之間交互影響與空間溢出效應(yīng)的系統(tǒng)闡述與機制檢驗?；诖耍狙芯拷⒘魍〝?shù)字化與制造業(yè)韌性（以下簡稱“兩系統(tǒng)”）的綜合評價指標(biāo)體系，運用耦合協(xié)調(diào)模型分析我國各地區(qū)兩系統(tǒng)協(xié)同發(fā)展現(xiàn)狀及時空演化特征，同時建立空間聯(lián)立方程模型對兩系統(tǒng)的交互影響及空間溢出效應(yīng)進行實證檢驗，以期豐富產(chǎn)業(yè)間協(xié)調(diào)發(fā)展的相關(guān)實證研究，并為進一步挖掘流通數(shù)字化潛力、增強制造業(yè)產(chǎn)業(yè)鏈韌性提供可資借鑒的政策啟示。

2 理論分析

2.1 流通數(shù)字化與制造業(yè)韌性的雙向影響機理

2.1.1 流通數(shù)字化對制造業(yè)韌性的影響

伴隨數(shù)字經(jīng)濟的快速發(fā)展，流通業(yè)的數(shù)字化轉(zhuǎn)型逐漸成為擊破制造業(yè)粗放式發(fā)展痛點的關(guān)鍵力量，驅(qū)動制造業(yè)產(chǎn)業(yè)鏈、供應(yīng)鏈韌性不斷增強。理由如下：

第一，流通業(yè)的數(shù)字化轉(zhuǎn)型能夠在更大程度上發(fā)揮由市場反向引領(lǐng)制造業(yè)柔性生產(chǎn)的先導(dǎo)作用。數(shù)字技術(shù)的廣泛應(yīng)用推動流通業(yè)產(chǎn)業(yè)鏈上下游的柔性化發(fā)展，去中心化的生產(chǎn)流通機制使得生產(chǎn)邊界日益模糊，生產(chǎn)制造企業(yè)的運營體系也隨之轉(zhuǎn)型［28］。其中，數(shù)字金融、金融科技融入流通領(lǐng)域有助于搭建數(shù)字供應(yīng)鏈金融平臺，以其規(guī)模效應(yīng)、長尾效應(yīng)、普惠效應(yīng)等優(yōu)勢拓寬制造業(yè)投融資渠道。流通業(yè)逐漸由去中心化階段過渡至融合化發(fā)展階段，以其更加完善的線上線下全渠道營銷平臺服務(wù)制造業(yè)生產(chǎn)、流通、消費各環(huán)節(jié)的全產(chǎn)業(yè)鏈良性互動，引領(lǐng)高端制造業(yè)動力變革，提升制造業(yè)企業(yè)抵御風(fēng)險的能力。

第二，流通基礎(chǔ)設(shè)施數(shù)字轉(zhuǎn)型的集約效應(yīng)和擴展效應(yīng)有助于推進制造業(yè)提質(zhì)增效，可以顯著提升制造業(yè)的創(chuàng)新與再造韌性。具體來說，新型流通設(shè)施不僅能夠加快制造業(yè)由勞動密集型向技術(shù)集約型轉(zhuǎn)變，促進裝備制造國產(chǎn)化，還能夠有效降低能源密集型制造業(yè)的能耗消耗和污染排放，建設(shè)更加完備的綠色制造體系，在可持續(xù)與精益化變革中發(fā)揮乘數(shù)效應(yīng)助力制造業(yè)全方位全鏈條改造［29］。此外，由于各地區(qū)資源分布的顯著差異，我國制造業(yè)供應(yīng)鏈的空間布局復(fù)雜冗長，而數(shù)字網(wǎng)絡(luò)平臺的應(yīng)用可以極大節(jié)約時間成本，提升制造業(yè)供應(yīng)效率，并為增強其產(chǎn)業(yè)鏈韌性創(chuàng)造條件。

第三，流通數(shù)字化的虛擬集聚效應(yīng)有助于推進制造業(yè)服務(wù)化進程。數(shù)字信息技術(shù)融入流通業(yè)能夠有效打破市場壁壘，形成更為靈活的產(chǎn)銷供應(yīng)模式，構(gòu)建供需雙方均為關(guān)鍵節(jié)點的價值網(wǎng)絡(luò)［28］。伴隨新零售、智慧物流等流通模式發(fā)展而產(chǎn)生的虛擬集聚效應(yīng)，一方面可以打破地理空間對商品運輸、信息傳遞的限制，節(jié)約搜尋成本和合約成本，優(yōu)化制造業(yè)分工水平，提升制造業(yè)生產(chǎn)效率［30］；另一方面可以極大程度提高生產(chǎn)制造過程中吸收消費端衍射的可能性，加快制造業(yè)產(chǎn)業(yè)鏈向“微笑曲線”兩端的高附加值領(lǐng)域邁進［31］。

2.1.2 制造業(yè)韌性對流通數(shù)字化的影響

制造業(yè)韌性的提升意味著產(chǎn)業(yè)鏈供應(yīng)鏈在遭受沖擊后的恢復(fù)能力增強，并且能探尋新的路徑以實現(xiàn)經(jīng)濟效益的穩(wěn)增長。制造業(yè)韌性的增強助力流通數(shù)字化轉(zhuǎn)型主要表現(xiàn)在以下方面：

第一，制造業(yè)韌性提升的重要表現(xiàn)之一就是供應(yīng)鏈韌性的提升，其對流通數(shù)字化改造有著極強的推動作用。制造業(yè)韌性的提升能夠通過自動扶梯效應(yīng)與其他產(chǎn)業(yè)融合發(fā)展，以自身為核心實現(xiàn)收益在經(jīng)濟體中的擴散［30］。制造業(yè)韌性增強往往意味著制造業(yè)產(chǎn)業(yè)結(jié)構(gòu)得以優(yōu)化且創(chuàng)新擴散能力逐漸增強，一方面能夠整合更多創(chuàng)新資源，在各地區(qū)強化其自身發(fā)展的比較優(yōu)勢，挖掘當(dāng)?shù)厣藤Q(mào)流通業(yè)數(shù)字化轉(zhuǎn)型潛能，為打通流通供銷創(chuàng)造條件；另一方面能夠提升對數(shù)字化、智能化平臺資源的吸納能力，直接擴大當(dāng)?shù)匦滦土魍ㄆ脚_的輻射范圍，以融合發(fā)展提升市場資源的配置效率。

第二，制造業(yè)韌性的增強既依賴于制造業(yè)企業(yè)在研發(fā)、管理、生產(chǎn)等各環(huán)節(jié)優(yōu)化升級，又可以為流通企業(yè)的數(shù)字化轉(zhuǎn)型提供更多場景支持。流通企業(yè)通過吸收制造業(yè)柔性生產(chǎn)、精益生產(chǎn)的精髓，構(gòu)建具有虛擬集成特點的大規(guī)模通信網(wǎng)絡(luò)系統(tǒng)，轉(zhuǎn)變傳統(tǒng)流通企業(yè)的經(jīng)營理念和運營模式，著力發(fā)展依托互聯(lián)網(wǎng)的私人訂制、網(wǎng)絡(luò)直銷以及云制造等新型流通模式，不僅可以提升供給體系和流通體系的適配性，為流通價值的創(chuàng)造提供必要支持，還能夠?qū)崿F(xiàn)“互聯(lián)網(wǎng)+流通+生產(chǎn)”的全面升級，形成商流、物流與信息流融合發(fā)展的供應(yīng)鏈服務(wù)體系。

第三，制造業(yè)韌性的增強勢必會瞄準(zhǔn)具有高技術(shù)含量和高附加值的戰(zhàn)略性新興產(chǎn)業(yè)、裝備制造業(yè)等高端制造業(yè)繼續(xù)發(fā)力，為流通基礎(chǔ)設(shè)施的智能化轉(zhuǎn)型提供技術(shù)支持。高端制造業(yè)一方面能夠以更高的產(chǎn)業(yè)關(guān)聯(lián)度和控制力發(fā)揮技術(shù)外溢效應(yīng)，為信息化、智能化發(fā)展?fàn)I造必要的創(chuàng)新氛圍［32］；另一方面能夠以更低的污染排放踐行制造業(yè)向低碳環(huán)保方向發(fā)展，提升智能制造的綠色化水平，在加快綠色物流、智慧商流服務(wù)體系的建設(shè)中實現(xiàn)供應(yīng)鏈數(shù)字化轉(zhuǎn)型，為打造可持續(xù)發(fā)展的流通環(huán)境助力。

基于上述分析，提出假設(shè)1：流通數(shù)字化與制造業(yè)韌性之間存在交互影響。

2.2 流通數(shù)字化與制造業(yè)韌性的空間溢出效應(yīng)

2.2.1 流通數(shù)字化的空間溢出

數(shù)字化技術(shù)的空間溢出與擴散效應(yīng)能夠幫助落后地區(qū)追趕先進生產(chǎn)技術(shù)［33］。那么，流通數(shù)字化發(fā)展的空間溢出至少存在以下原因：一是隨互聯(lián)網(wǎng)發(fā)展而來的數(shù)字化流通供銷平臺的應(yīng)用使得地區(qū)間傳統(tǒng)邊界限制逐漸弱化，商品貿(mào)易往來的便利化程度得以提升，流通供銷紅利共享推動各地區(qū)共建數(shù)字化流通網(wǎng)絡(luò)。二是數(shù)字金融與商貿(mào)流通業(yè)的融合能夠提升當(dāng)?shù)亓魍I(yè)的融資效率，伴隨而來的是金融資源的跨地區(qū)流動，幫助鄰近地區(qū)獲得更多資金支持，讓智慧物流、智能倉儲等更好地服務(wù)于實體經(jīng)濟的高質(zhì)量發(fā)展，融資外部性釋放的空間溢出得以顯現(xiàn)。三是在數(shù)字基礎(chǔ)設(shè)施的運營過程中，數(shù)據(jù)要素的廣泛傳播調(diào)控知識要素的跨區(qū)域流動，基于知識溢出效應(yīng)，各地區(qū)技術(shù)人才的相互學(xué)習(xí)、科研機構(gòu)的深度合作、各類資源的頻繁交換有助于縮小地區(qū)間數(shù)據(jù)要素的稟賦差異［34］。

2.2.2 制造業(yè)韌性的空間溢出

隨著地區(qū)間經(jīng)濟社會關(guān)聯(lián)日益加深，制造業(yè)韌性的增強受益于鄰近地區(qū)的集聚效應(yīng)［18］。其空間溢出至少來源于以下途徑；一是制造業(yè)產(chǎn)業(yè)鏈條的延伸拓展以及產(chǎn)業(yè)分工逐漸向?qū)I(yè)化、精細化轉(zhuǎn)型，各地區(qū)在產(chǎn)業(yè)鏈上下游合作密切，形成前向后向關(guān)聯(lián)，推動產(chǎn)業(yè)一體化發(fā)展。二是隨著制造業(yè)內(nèi)競爭日益激烈，為適應(yīng)全球競爭的需要，技術(shù)密集型制造業(yè)應(yīng)運而生，其利用信息外溢效應(yīng)形成的技術(shù)生產(chǎn)網(wǎng)絡(luò)推動技術(shù)、資本等要素流動，使得地區(qū)間交往日益密切［35］。三是基于各地區(qū)產(chǎn)業(yè)招商引資的競爭效應(yīng)和示范效應(yīng)，本地高技術(shù)制造業(yè)的快速發(fā)展吸引周邊地區(qū)競相模仿，引進同類產(chǎn)業(yè)以保證地方政府具備競爭優(yōu)勢［36］。

2.2.3 流通數(shù)字化與制造業(yè)韌性的交互空間溢出

由于流通數(shù)字化轉(zhuǎn)型和制造業(yè)韌性增強均需要高技術(shù)人才的引進，因此極易形成本地稀缺型高技能人才向具有更優(yōu)人才引進政策的鄰近省份流動的局面，人才競爭加劇了各地區(qū)市場競爭壓力，這種回波效應(yīng)將對本地就業(yè)水平和產(chǎn)業(yè)發(fā)展帶來消極影響［37］。具體而言，若鄰近地區(qū)制造業(yè)企業(yè)較本地率先獲得技術(shù)突破或?qū)崿F(xiàn)產(chǎn)業(yè)價值鏈大幅升級，會產(chǎn)生鄰近地區(qū)制造業(yè)企業(yè)人員工資增加、吸引就業(yè)等直接效應(yīng)，同時間接引發(fā)其對流通服務(wù)智慧化升級的迫切需要［38］。為此，基于優(yōu)質(zhì)的人才引進政策或自身價值的實現(xiàn)，本地高技術(shù)人才涌入該地區(qū)為其流通供銷數(shù)字化轉(zhuǎn)型貢獻技術(shù)支持，直接導(dǎo)致本地技能人才流失，從而產(chǎn)生負向的空間溢出效應(yīng)；同樣，鄰近地區(qū)流通數(shù)字化發(fā)展也會因人才、技術(shù)、信息等資源的擠占，對本地制造業(yè)韌性提升產(chǎn)生負向的空間溢出效應(yīng)。長此以往，各地區(qū)人力資本質(zhì)量、產(chǎn)業(yè)創(chuàng)新績效、經(jīng)濟發(fā)展水平的差異化將逐漸擴大，為區(qū)域間協(xié)調(diào)聯(lián)動發(fā)展埋下隱患。

基于以上分析，提出假設(shè)2：流通數(shù)字化與制造業(yè)韌性均存在空間溢出效應(yīng)且具有地區(qū)交互影響。即鄰近地區(qū)流通數(shù)字化改造能夠促進本地流通數(shù)字化發(fā)展，但對本地制造業(yè)韌性提升產(chǎn)生阻礙作用；鄰近地區(qū)制造業(yè)韌性增強能夠促進本地制造業(yè)韌性提升，但對本地流通數(shù)字化轉(zhuǎn)型具有抑制作用。

3 研究設(shè)計

3.1 指標(biāo)體系與數(shù)據(jù)來源

遵循系統(tǒng)性、科學(xué)性、可得性等原則，本研究以我國31 個省份（未含港澳臺地區(qū)）2014—2019 年的有關(guān)數(shù)據(jù)為樣本，參考郝愛民［39］、魯飛宇等［40］的方法，分別構(gòu)建流通數(shù)字化和制造業(yè)韌性的綜合評價指標(biāo)體系，如表1 所示。具體而言，分別從流通供銷數(shù)字化、數(shù)字供應(yīng)鏈金融、流通設(shè)施數(shù)字化3個維度衡量流通數(shù)字化發(fā)展水平；從抵抗力、恢復(fù)力、再造力3 個方面對制造業(yè)韌性進行測度。二級指標(biāo)金融科技公司數(shù)量借鑒宋敏等［41］的方法進行檢索，其余各指標(biāo)數(shù)據(jù)均來源于歷年的《中國統(tǒng)計年鑒》《中國第三產(chǎn)業(yè)統(tǒng)計年鑒》《中國工業(yè)統(tǒng)計年鑒》《中國稅務(wù)年鑒》《中國高技術(shù)產(chǎn)業(yè)統(tǒng)計年鑒》《中國電子商務(wù)發(fā)展指數(shù)報告》《北京大學(xué)數(shù)字普惠金融指數(shù)》。本研究中商貿(mào)流通業(yè)主要包括批發(fā)和零售業(yè)、交通運輸、倉儲和郵政業(yè)以及住宿和餐飲業(yè)。部分缺失數(shù)據(jù)運用插值法進行補充。同時為消除價格因素的影響，將涉及固定資產(chǎn)投資的指標(biāo)和與產(chǎn)值相關(guān)的經(jīng)濟變量分別以2013 年為基期的固定資產(chǎn)投資價格指數(shù)和價格指數(shù)進行平減處理。

表1 流通數(shù)字化與制造業(yè)韌性系統(tǒng)評價指標(biāo)體系

3.2 研究方法與模型構(gòu)建

3.2.1 流通數(shù)字化與制造業(yè)韌性綜合評價模型

在構(gòu)建評價指標(biāo)體系的基礎(chǔ)上，為了盡可能避免主觀因素造成的選擇偏誤，運用熵值法分別測算兩系統(tǒng)的綜合發(fā)展水平。首先，利用正向標(biāo)準(zhǔn)化法對數(shù)據(jù)進行極差標(biāo)準(zhǔn)化處理，以消除各指標(biāo)的量綱和數(shù)量級影響。具體公式如下：

式（1）中：i、j、t分別為地區(qū)、指標(biāo)和年份；Xijt和Sijt分別為原始數(shù)據(jù)和標(biāo)準(zhǔn)化后的數(shù)據(jù)；Xmax和Xmin分別為樣本期間第j項指標(biāo)的最大值和最小值。

其次，測算信息熵冗余度與指標(biāo)權(quán)重。計算公式如下：

式（2）～（4）中：k、m、n分別為樣本的時間范圍、地區(qū)總數(shù)和指標(biāo)總數(shù)；djt為信息熵冗余度；wjt為各指標(biāo)的最終權(quán)重。

最后，利用線性加權(quán)法得到各年度流通數(shù)字化（DCit）與制造業(yè)韌性（MRit）的綜合評價值。計算公式如下：

3.2.2 耦合協(xié)調(diào)模型

基于物理學(xué)的耦合概念衡量兩系統(tǒng)的耦合協(xié)調(diào)度，以反映兩系統(tǒng)之間的相互作用與依賴程度。模型設(shè)定如下：

式（7）～（9）中：Cit為i地區(qū)第t年兩系統(tǒng)的耦合度，其值介于0～1 之間，越接近1 表明兩系統(tǒng)總體發(fā)展越有序；Tit為兩系統(tǒng)綜合協(xié)調(diào)指數(shù)；α和β為待定系數(shù)，且滿足α+β=1，考慮到流通數(shù)字化與制造業(yè)韌性的互促互補關(guān)系，故將兩者視為同等重要地位，即α、β均取值0.5；Dit為兩系統(tǒng)的耦合協(xié)調(diào)度，參考葛鵬飛等［42］的做法，將耦合協(xié)調(diào)度進行級別劃分，如表2 所示。

表2 流通數(shù)字化與制造業(yè)韌性二元系統(tǒng)耦合協(xié)調(diào)度級別劃分

3.2.3 空間自相關(guān)檢驗

為了識別兩系統(tǒng)耦合協(xié)調(diào)度的空間相關(guān)性與集聚特征，運用全局自相關(guān)莫蘭指數(shù)（Moran'sI）和局部自相關(guān)莫蘭指數(shù)（Local Moran's I）進行分析。指數(shù)計算如下：

式（10）（11）中：Di和Dh為i地區(qū)和h地區(qū)的兩系統(tǒng)耦合協(xié)調(diào)度；為m個地區(qū)兩系統(tǒng)的耦合協(xié)調(diào)度均值；wih為地理二值的空間權(quán)重矩陣，當(dāng)?shù)貐^(qū)i與h不相鄰或i=h時，wih=0，否則wih=1；全局莫蘭指數(shù)取值范圍為［-1,1］，值為正時表明兩地區(qū)空間正相關(guān)，反之則呈現(xiàn)空間負相關(guān)，且指數(shù)的絕對值越大表明空間關(guān)聯(lián)性和集聚性越強；局部莫蘭指數(shù)取值范圍也為［-1,1］，值大于0 表明地區(qū)間呈現(xiàn)相似的空間集聚特征，即H-H 型（高-高集聚）或L-L 型（低-低集聚）；反之則呈現(xiàn)非相似的空間集聚特征，即H-L 型（高-低集聚）或L-H型（低-高集聚）。

3.2.4 空間聯(lián)立方程模型

為了檢驗流通數(shù)字化與制造業(yè)韌性之間的交互影響與空間溢出效應(yīng)，采用空間聯(lián)立方程模型，運用廣義三階段空間最小二乘估計法（GS3SLS）進行估計。此方法的優(yōu)點在于可以在一定程度上緩解方程聯(lián)立性偏誤、遺漏變量等引起的內(nèi)生性問題。模型形式如下：

式（12）（13）中：μ為空間相關(guān)系數(shù)，反映流通數(shù)字化和制造業(yè)韌性的空間溢出效應(yīng)；wih為地區(qū)間的空間關(guān)聯(lián)，當(dāng)i=h或時wih=0，當(dāng)i≠h或時wih=1/dih；dih為地區(qū)i與h之間的距離；為地理閾值；控制變量Zit包括經(jīng)濟發(fā)展水平，以平減處理后的地區(qū)生產(chǎn)總值衡量，人口密度以地區(qū)人口總數(shù)與行政區(qū)面積的比值衡量，交通基礎(chǔ)設(shè)施以公路運輸路線長度衡量；貿(mào)易規(guī)模以進出口總額與地區(qū)生產(chǎn)總值的比值衡量；控制變量Xit在Zit的基礎(chǔ)上增加城市失業(yè)率變量，以城鎮(zhèn)登記失業(yè)率衡量。

4 實證結(jié)果分析

4.1 流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展分析

如表3 所示，2014—2019 年我國各地區(qū)流通數(shù)字化與制造業(yè)韌性的耦合協(xié)調(diào)度逐年遞增，由2014年的0.379 增長到2019 年的0.455，年均增長率為3.72%?？傮w上，兩系統(tǒng)處于中度耦合協(xié)調(diào)階段，在自身的平穩(wěn)發(fā)展中深化良性互動。

具體來看，由于不同地區(qū)經(jīng)濟基礎(chǔ)、產(chǎn)業(yè)發(fā)展存在一定差異，因此兩系統(tǒng)的耦合協(xié)調(diào)程度具有區(qū)域異質(zhì)性特征。為探究兩系統(tǒng)協(xié)同發(fā)展的空間演進規(guī)律，借鑒唐曉華等［43］的分類方法，將我國31 個省份劃分為東北、東部沿海、南部沿海、北部沿海、黃河中游、長江中游、大西南和大西北八大經(jīng)濟區(qū)，直觀清晰地分析各區(qū)域的兩系統(tǒng)耦合協(xié)調(diào)差異。由表3 可知，東部沿海、南部沿海、北部沿海、長江中游經(jīng)濟區(qū)的兩系統(tǒng)耦合協(xié)調(diào)發(fā)展水平始終位列前四。其中，東部沿海經(jīng)濟區(qū)兩系統(tǒng)整體處于高度耦合協(xié)調(diào)階段，協(xié)調(diào)發(fā)展水平遠超其他經(jīng)濟區(qū)，主要源于其以上海為核心打造產(chǎn)業(yè)數(shù)字服務(wù)平臺、布局?jǐn)?shù)字商經(jīng)濟集聚區(qū)，帶動流通業(yè)數(shù)字化轉(zhuǎn)型，依托于浙西南、蘇西北制造業(yè)基地，在提升制造業(yè)產(chǎn)業(yè)鏈、供應(yīng)韌性的同時大力建設(shè)中國特色現(xiàn)代流通體系，促進兩系統(tǒng)協(xié)調(diào)發(fā)展；長江中游經(jīng)濟區(qū)以其增速優(yōu)勢于2018 年趕超北部沿海經(jīng)濟區(qū)躋身第3 位，得益于長江中游地區(qū)不斷推進市場一體化、完善智慧流通工程，輔助支持制造業(yè)服務(wù)化轉(zhuǎn)型，加速推進兩系統(tǒng)良性互動；而黃河中游、大西南、東北、大西北經(jīng)濟區(qū)的兩系統(tǒng)耦合協(xié)調(diào)水平均低于31 個省份平均水平，特別是大西北經(jīng)濟區(qū)制造業(yè)水平整體偏低、應(yīng)對風(fēng)險能力整體偏弱、智慧物流建設(shè)起步較晚，兩系統(tǒng)始終處于低度耦合協(xié)調(diào)階段，雖整體呈現(xiàn)出逐年上升的趨勢，但距離實現(xiàn)高度甚至極度耦合協(xié)調(diào)仍有較大差距。

進一步觀察各經(jīng)濟區(qū)兩系統(tǒng)耦合協(xié)調(diào)的年均增速發(fā)現(xiàn)（見表3），八大經(jīng)濟區(qū)中年均增長率最高的為長江中游經(jīng)濟區(qū)，增速達5.53%；大西南和大西北經(jīng)濟區(qū)緊隨其后，增速分別為5.21%和4.69%；排名末尾的為東北經(jīng)濟區(qū)，增速為1.48%。可見，多數(shù)處于低度和中度耦合協(xié)調(diào)階段的經(jīng)濟區(qū)的兩系統(tǒng)耦合協(xié)調(diào)度年均增長率相對較高，表明我國各經(jīng)濟區(qū)兩系統(tǒng)耦合協(xié)調(diào)發(fā)展空間分布不均的局面逐漸緩解，但仍具有較大的提升空間。

表3 我國區(qū)域流通數(shù)字化與制造業(yè)韌性耦合協(xié)調(diào)水平

表3（續(xù)）

4.2 流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性協(xié)同發(fā)展分析

如圖1 所示。就31 個省份整體而言，流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性的耦合協(xié)調(diào)度呈現(xiàn)出穩(wěn)步發(fā)展的態(tài)勢，其值均介于0.3～0.5 之間，處于中度耦合協(xié)調(diào)階段，仍有待進一步提升。從圖1 中均值線的高度來看，流通設(shè)施數(shù)字化與制造業(yè)韌性的耦合協(xié)調(diào)度均值相對最低，說明我國流通領(lǐng)域數(shù)字基礎(chǔ)設(shè)施建設(shè)與制造業(yè)韌性的協(xié)調(diào)適配有待提高；數(shù)字供應(yīng)鏈金融與制造業(yè)韌性的耦合協(xié)調(diào)度增勢尤為明顯，2014—2019 年的年均增長率達11.6%，表明數(shù)字供應(yīng)鏈金融子系統(tǒng)在推動流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展中最具潛力，反映出大力發(fā)展數(shù)字普惠金融、助推金融科技進步在實現(xiàn)產(chǎn)業(yè)融合創(chuàng)新中意義重大。

圖1 我國31 個省份流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性耦合協(xié)調(diào)水平對比

如表4 所示，大部分經(jīng)濟區(qū)的數(shù)字供應(yīng)鏈金融與制造業(yè)韌性的耦合協(xié)調(diào)均值最高；而東部沿海與南部沿海經(jīng)濟區(qū)不同于其他區(qū)域，其流通供銷數(shù)字化與制造業(yè)韌性的耦合協(xié)調(diào)均值高于數(shù)字供應(yīng)鏈金融與制造業(yè)韌性的耦合協(xié)調(diào)均值，表明東南沿海地區(qū)依托其經(jīng)濟優(yōu)勢加速電商、物流在供銷流通環(huán)節(jié)的廣泛應(yīng)用，相對更加適應(yīng)制造業(yè)發(fā)展韌性的提升。從耦合協(xié)調(diào)度年均增長率來看，大西北、大西南經(jīng)濟區(qū)數(shù)字供應(yīng)鏈金融與制造業(yè)韌性的耦合協(xié)調(diào)增長率位列前二，分別為15.95%和13.68%，可見在經(jīng)濟相對欠發(fā)達的西部地區(qū)，以供應(yīng)鏈金融為核心激發(fā)流通數(shù)字化轉(zhuǎn)型并與制造業(yè)發(fā)展韌性提升形成良性互動頗具潛力。

表4 2014—2019 年我國區(qū)域流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性耦合協(xié)調(diào)水平對比

4.3 流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展的空間自相關(guān)檢驗

由表5 可知，考察期間31 個省份兩系統(tǒng)協(xié)同發(fā)展的全局莫蘭指數(shù)均大于0.4，且在1%的水平上顯著，表明兩系統(tǒng)協(xié)同發(fā)展整體呈現(xiàn)出顯著的空間自相關(guān)性，空間集聚態(tài)勢不斷增強。為了進一步探究各地區(qū)兩系統(tǒng)協(xié)同發(fā)展的空間相關(guān)性差異，將局部莫蘭指數(shù)分析得出的LISA 集聚情況予以呈現(xiàn)，可以發(fā)現(xiàn)山東、江蘇、安徽、上海、浙江、福建始終位于H-H 型集聚區(qū)；江西自2016 年起從L-H 型集聚區(qū)轉(zhuǎn)變至H-H 型集聚區(qū)；新疆、青海、甘肅始終位于L-L 型集聚區(qū)，而內(nèi)蒙古自2015 年起加入L-L型集聚區(qū)，西藏自2017 年起退出L-L 型集聚區(qū)；四川始終位于H-L 型集聚區(qū)；其余省份LISA 集聚結(jié)果不顯著，表明其與相鄰省份間的聯(lián)系相對較弱。由此可見，31 個省份兩系統(tǒng)的協(xié)同發(fā)展主要聚集于H-H 型集聚區(qū)和L-L 型集聚區(qū)，具有一定的兩極分化趨勢：東部沿海經(jīng)濟區(qū)及其部分鄰近省份兩系統(tǒng)的自身耦合協(xié)調(diào)度較高，且能夠通過技術(shù)擴散效應(yīng)、知識溢出效應(yīng)、產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化等對周邊地區(qū)形成較強的輻射帶動作用；大西北經(jīng)濟區(qū)則由于產(chǎn)業(yè)發(fā)展相對滯后、交通運輸設(shè)施較匱乏、科技創(chuàng)新人才短缺等問題，區(qū)域內(nèi)各省份間兩系統(tǒng)合作基礎(chǔ)欠佳、協(xié)調(diào)發(fā)展能力不足，因此分布在L-L 型集聚區(qū)。

表5 我國區(qū)域流通數(shù)字化與制造業(yè)韌性協(xié)同發(fā)展的空間自相關(guān)檢驗結(jié)果

4.4 流通數(shù)字化與制造業(yè)韌性的交互影響與空間溢出效應(yīng)檢驗

由表6 可知，在地理二值和地理閾值兩種空間矩陣下，各變量的回歸系數(shù)均未出現(xiàn)符號變化，表明顯著性水平基本穩(wěn)健。以地理二值空間權(quán)重矩陣的GS3SLS 回歸結(jié)果為例，制造業(yè)韌性對流通數(shù)字化的估計系數(shù)為正且在1%的水平上顯著，制造業(yè)韌性程度每增加1%，可以促進流通數(shù)字化水平提升1.55%；反之，流通數(shù)字化對制造業(yè)韌性的估計系數(shù)也在1%的水平下顯著為正，流通數(shù)字化水平每增加1%，可以使得制造業(yè)韌性程度提升0.49%。從空間溢出效應(yīng)來看，流通數(shù)字化與制造業(yè)韌性的空間滯后性的溢出效應(yīng)分別為0.019、0.020，其空間交互溢出效應(yīng)分別為-0.012、-0.031，且均在5%或1%的水平上顯著。由此可見，實證結(jié)果驗證了本研究提出的假設(shè)1 和假設(shè)2。

表6 2014—2019 年31 個省份流通數(shù)字化與制造業(yè)韌性交互影響與空間溢出效應(yīng)檢驗結(jié)果

表6（續(xù)）

為了進一步探究不同的流通數(shù)字化轉(zhuǎn)型方式與制造業(yè)韌性增強的交互影響是否存在異質(zhì)性，將流通數(shù)字化各子系統(tǒng)的綜合評價值進行再次回歸。如表7 的地理二值空間權(quán)重矩陣下的實證結(jié)果所示，制造業(yè)韌性對流通供銷數(shù)字化和流通設(shè)施數(shù)字化的估計系數(shù)分別為1.529、1.755，且均在1%的水平上顯著；流通數(shù)字化各子系統(tǒng)對制造業(yè)韌性的估計系數(shù)為正，且均在1%或10%的水平上顯著，流通供銷數(shù)字化、數(shù)字供應(yīng)鏈金融和流通設(shè)施數(shù)字化每提升1 個百分點，制造業(yè)韌性分別增強0.615、0.144、0.563 個百分點。由此可見，制造業(yè)韌性的增強對流通基礎(chǔ)設(shè)施數(shù)字化改造的促進作用最強，流通供銷環(huán)節(jié)的數(shù)字化變革對提升制造業(yè)韌性的效果最優(yōu)。從空間互動來看，流通設(shè)施數(shù)字化與制造業(yè)韌性對鄰近地區(qū)的正向空間溢出效應(yīng)較為明顯，鄰近地區(qū)流通設(shè)施數(shù)字化程度每提升1 個百分點，本地流通設(shè)施數(shù)字化水平提高0.029 個百分點；鄰近地區(qū)制造業(yè)韌性每增強1 個百分點，本地制造業(yè)韌性增強0.021 個百分點。數(shù)字供應(yīng)鏈金融與制造業(yè)韌性的交互空間溢出效應(yīng)最弱，即鄰近地區(qū)數(shù)字供應(yīng)鏈金融發(fā)展對本地打造制造業(yè)韌性的資源擠占較??；反之亦是如此。因此，就31 個省份整體而言，各地區(qū)推動流通基礎(chǔ)設(shè)施智慧化轉(zhuǎn)型能夠有效借助信息技術(shù)的外溢效應(yīng)輻射帶動周邊地區(qū)流通設(shè)施改造，并形成區(qū)域間產(chǎn)業(yè)聯(lián)動發(fā)展的良好局面；就部分經(jīng)濟欠發(fā)達地區(qū)而言，各省份的相互輻射帶動能力有限，在加快新型基礎(chǔ)設(shè)施建設(shè)的同時可以考慮加大金融科技、數(shù)字普惠金融在流通領(lǐng)域的應(yīng)用深度和廣度，更好地推動實體產(chǎn)業(yè)優(yōu)化升級并向產(chǎn)業(yè)鏈高端延伸。

表7 2014—2019 年31 個省份流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性交互影響與空間溢出效應(yīng)檢驗結(jié)果

5 結(jié)論與政策建議

5.1 研究結(jié)論

（1）2014—2019 年間，我國31 個省份流通數(shù)字化與制造業(yè)韌性的協(xié)同發(fā)展度逐年遞增，呈現(xiàn)出由低度耦合協(xié)調(diào)向中度耦合協(xié)調(diào)過渡的趨勢，且存在明顯的區(qū)域異質(zhì)性特征，即東部高于西部、沿海優(yōu)于內(nèi)陸的分布狀態(tài)。具體地，從均值來看，東部沿海和南部沿海經(jīng)濟區(qū)始終位居前二，達到高度耦合協(xié)調(diào)水平，而黃河中游、大西南、東北、大西北經(jīng)濟區(qū)均低于31 個省份平均水平；從增速來看，長江中游經(jīng)濟區(qū)年均增長率最高，大西南和大西北經(jīng)濟區(qū)緊隨其后，可見區(qū)域間差異正在逐漸縮小。

（2）從流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性耦合協(xié)調(diào)情況來看，我國流通設(shè)施數(shù)字化與制造業(yè)韌性耦合協(xié)調(diào)度的均值最低，但在考察期內(nèi)具有穩(wěn)步發(fā)展的態(tài)勢；數(shù)字供應(yīng)鏈金融與制造業(yè)韌性耦合協(xié)調(diào)度的年均增長率最高，尤其是在經(jīng)濟相對欠發(fā)達的西部地區(qū)增勢尤為明顯。

（3）由空間效應(yīng)分析可知，31 個省份流通數(shù)字化與制造業(yè)韌性的協(xié)同發(fā)展具有顯著的空間自相關(guān)性，且兩者之間存在雙向空間溢出效應(yīng)，即本地與鄰近地區(qū)流通數(shù)字化水平顯著正相關(guān)，鄰近地區(qū)流通數(shù)字化發(fā)展抑制本地制造業(yè)韌性增強；本地與鄰近地區(qū)制造業(yè)韌性顯著正相關(guān)，鄰近地區(qū)制造業(yè)韌性提升抑制本地流通數(shù)字化發(fā)展。

（4）通過分組回歸發(fā)現(xiàn)，31 個省份流通數(shù)字化各子系統(tǒng)與制造業(yè)韌性的交互影響和空間效應(yīng)存在一定的顯著性差異，制造業(yè)韌性提升對流通基礎(chǔ)設(shè)施數(shù)字化改造的促進作用最強，流通供銷環(huán)節(jié)的數(shù)字化變革對增強制造業(yè)韌性的效果最優(yōu)，數(shù)字設(shè)施數(shù)字化與制造業(yè)韌性的空間溢出效應(yīng)較為明顯。

5.2 政策建議

第一，統(tǒng)籌流通數(shù)字化與制造業(yè)韌性耦合發(fā)展的協(xié)同建設(shè)，以第二產(chǎn)業(yè)和第三產(chǎn)業(yè)良性互動引領(lǐng)經(jīng)濟高質(zhì)量發(fā)展。一方面，發(fā)揮數(shù)字技術(shù)在流通領(lǐng)域的長尾效應(yīng)，為制造業(yè)創(chuàng)新發(fā)展提供信息、技術(shù)、渠道支持，增強制造業(yè)產(chǎn)業(yè)鏈供應(yīng)鏈韌性；另一方面，發(fā)揮制造業(yè)效率提升的資源整合效應(yīng)，在推進制造業(yè)服務(wù)化發(fā)展中打破數(shù)字化資源流動壁壘，進一步推進流通數(shù)字化轉(zhuǎn)型。此外，由于當(dāng)前我國流通數(shù)字化與制造業(yè)韌性耦合協(xié)調(diào)發(fā)展存在區(qū)域差異，有必要強化地區(qū)間協(xié)同發(fā)展意識，優(yōu)化兩系統(tǒng)耦合協(xié)調(diào)機制，以減小地區(qū)間發(fā)展的不平衡。通過打破各經(jīng)濟區(qū)之間的制度壁壘，建立以東南沿海為代表的高水平經(jīng)濟區(qū)和以大西北、東北為代表的低水平經(jīng)濟區(qū)有效合作體制機制，發(fā)揮高水平經(jīng)濟區(qū)的示范帶動作用，以“數(shù)字流通+智能制造”的融合發(fā)展模式，為各經(jīng)濟區(qū)營造良好的科技創(chuàng)新氛圍，在智慧城市建設(shè)中帶動城市經(jīng)濟可持續(xù)發(fā)展。

第二，基于比較優(yōu)勢，因地制宜發(fā)揮流通數(shù)字化各子系統(tǒng)的區(qū)域優(yōu)勢，加快構(gòu)建“雙循環(huán)”新發(fā)展格局。實證結(jié)果表明我國數(shù)字供應(yīng)鏈金融子系統(tǒng)與制造業(yè)韌性耦合協(xié)調(diào)增長率最高，因此，應(yīng)進一步激發(fā)金融科技、數(shù)字普惠金融等在流通領(lǐng)域的發(fā)展?jié)摿?，深化智慧供銷渠道建設(shè)，緩解經(jīng)濟欠發(fā)達地區(qū)制造業(yè)融資約束，在打通國內(nèi)生產(chǎn)、分配、流通、消費各環(huán)節(jié)的同時增強制造業(yè)產(chǎn)業(yè)鏈韌性，以重塑后疫情時代我國產(chǎn)業(yè)的國際競爭優(yōu)勢，推進國內(nèi)國際雙循環(huán)相互促進的新發(fā)展格局建設(shè)。當(dāng)前，東南沿海經(jīng)濟區(qū)擁有相對較多的高新技術(shù)資源，可率先運用大數(shù)據(jù)、區(qū)塊鏈、云計算等數(shù)字技術(shù)實現(xiàn)傳統(tǒng)產(chǎn)業(yè)變革，以數(shù)字化轉(zhuǎn)型引領(lǐng)高質(zhì)量發(fā)展；對西北、西南和東北經(jīng)濟區(qū)而言，不能過度投資以追求短期內(nèi)經(jīng)濟快速發(fā)展而忽視勞動密集型產(chǎn)業(yè)的重要作用，應(yīng)依托“一帶一路”建設(shè)，充分挖掘當(dāng)?shù)禺a(chǎn)業(yè)發(fā)展的比較優(yōu)勢，制定適合當(dāng)?shù)匕l(fā)展的產(chǎn)業(yè)政策，為經(jīng)濟相對欠發(fā)達地區(qū)探索“彎道超車”的可行方案。

第三，立足空間思維，促進資源的跨區(qū)流動，加快城市群一體化建設(shè)，發(fā)揮大國城市的規(guī)模經(jīng)濟效應(yīng)。實證結(jié)果表明我國流通數(shù)字化和制造業(yè)韌性的空間滯后項均表現(xiàn)為正向空間溢出效應(yīng)，為此，應(yīng)依托具備地緣優(yōu)勢的核心城市發(fā)展城市群經(jīng)濟，通過頂層設(shè)計推動區(qū)域一體化發(fā)展和市場化改革，構(gòu)建互利共生的區(qū)域多產(chǎn)業(yè)協(xié)同、多鏈條融合發(fā)展格局。一方面，加速各地區(qū)技術(shù)人才交流，利用知識和信息的外溢效應(yīng)提升數(shù)字化要素的空間配置效率；另一方面，運用鄰近區(qū)域資源優(yōu)勢互補，建立跨區(qū)協(xié)調(diào)發(fā)展的組織協(xié)同機構(gòu)與運行體系，將流通數(shù)字化與制造業(yè)韌性耦合發(fā)展推向更高水平。而針對我國流通數(shù)字化與制造業(yè)韌性之間存在負向交互空間溢出效應(yīng)的現(xiàn)狀，各級地方政府在制定相關(guān)政策時，更應(yīng)考慮將各城市置于同一空間體系中，以整體空間經(jīng)濟結(jié)構(gòu)發(fā)展為前提推動資本、技術(shù)、信息的跨區(qū)自由流動，加強與周邊地區(qū)在產(chǎn)業(yè)鏈供應(yīng)鏈上下游的協(xié)同合作。