李艷偉 副教授 何任杰

（中國民航大學(xué)經(jīng)濟(jì)與管理學(xué)院 天津 300300）

引言

京津冀地區(qū)在我國社會(huì)經(jīng)濟(jì)發(fā)展過程中扮演著舉足輕重的角色，也是我國北方商貿(mào)流通業(yè)發(fā)展最活躍的地區(qū)（逄金玉等，2014），隨著我國社會(huì)的主要矛盾發(fā)生轉(zhuǎn)變，促進(jìn)商貿(mào)流通業(yè)平衡和充分發(fā)展，更有利于滿足人民美好的物質(zhì)文化需要，在京津冀協(xié)同發(fā)展上升到國家戰(zhàn)略的時(shí)代背景下，促進(jìn)京津冀地區(qū)商貿(mào)流通業(yè)的協(xié)同發(fā)展更有利于提高京津冀區(qū)域經(jīng)濟(jì)發(fā)展水平。2015年，中共中央政治局審議并通過的《京津冀協(xié)同發(fā)展規(guī)劃綱要》中指出，要促進(jìn)京津冀地區(qū)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，實(shí)現(xiàn)京津冀三地產(chǎn)業(yè)融合和協(xié)同發(fā)展（蔣敏娟等，2019），商貿(mào)流通業(yè)作為第三產(chǎn)業(yè)支柱和國民經(jīng)濟(jì)發(fā)展的引導(dǎo)性產(chǎn)業(yè)，在京津冀區(qū)域經(jīng)濟(jì)發(fā)展過程中扮演著承前啟后的作用，京津冀地區(qū)商貿(mào)流通業(yè)總產(chǎn)值占京津冀全行業(yè)總產(chǎn)值的近20%，其快速的發(fā)展趨勢和廣闊的發(fā)展前景已成為京津冀協(xié)同發(fā)展的重要產(chǎn)業(yè)動(dòng)力（張婷，2018）。雖然京津冀商貿(mào)流通業(yè)規(guī)模大，但近幾年增長速度放緩，三地的發(fā)展水平差距較大，因此，本文利用空間統(tǒng)計(jì)分析法，研究京津冀商貿(mào)流通業(yè)發(fā)展的時(shí)空差異特性以及受經(jīng)濟(jì)社會(huì)各類因素影響的顯著程度，并對其商貿(mào)流通業(yè)提出合理發(fā)展建議。

京津冀商貿(mào)流通業(yè)發(fā)展規(guī)模變化及發(fā)展水平對比

（一）京津冀商貿(mào)流通業(yè)發(fā)展規(guī)模變化

商貿(mào)流通業(yè)是以批發(fā)和零售業(yè)為主，包含飲食、倉儲(chǔ)和交通運(yùn)輸業(yè)的綜合服務(wù)型產(chǎn)業(yè)（曾蓼等，2018），本文用社會(huì)消費(fèi)品零售額作為評價(jià)指標(biāo)，衡量商貿(mào)流通業(yè)的發(fā)展水平。一般認(rèn)為，社會(huì)消費(fèi)品零售額越高，則商貿(mào)流通業(yè)發(fā)展水平也越高（任其俊，2017）。如圖1所示，從2009-2018年京津冀地區(qū)社會(huì)消費(fèi)品零售額變化情況來看，河北全省的商貿(mào)流通業(yè)規(guī)模高于北京和天津，天津的規(guī)模僅為北京的一半，最近三年京津冀三地的商貿(mào)流通業(yè)發(fā)展放緩，2018年天津出現(xiàn)負(fù)增長。

（二）京津冀商貿(mào)流通業(yè)發(fā)展水平對比

北京、天津和河北三地的商貿(mào)流通業(yè)發(fā)展存在明顯差異（溫天力，2017），北京市是我國北方的經(jīng)濟(jì)中心，密集的人口和快速發(fā)展的經(jīng)濟(jì)為北京商貿(mào)流通業(yè)發(fā)展奠定了基礎(chǔ)，商貿(mào)流通業(yè)已經(jīng)成為北京的第一大服務(wù)產(chǎn)業(yè)。而天津市近年來的經(jīng)濟(jì)增長速度有所放緩，產(chǎn)業(yè)結(jié)構(gòu)失衡現(xiàn)象日趨嚴(yán)重，天津市的發(fā)展仍以傳統(tǒng)重工業(yè)為主，而輕視了以商貿(mào)流通業(yè)為主的第三產(chǎn)業(yè)發(fā)展，加快天津市產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型，更有益于其流通產(chǎn)業(yè)的發(fā)展。如圖2所示，河北省整體的商貿(mào)流通業(yè)規(guī)模雖然最高，但各個(gè)地級市的規(guī)模發(fā)展水平較低，且差異較大，河北仍主要發(fā)展傳統(tǒng)產(chǎn)業(yè)，以交通業(yè)、服務(wù)業(yè)為主的第三產(chǎn)業(yè)發(fā)展相對落后，也間接導(dǎo)致河北各市商貿(mào)流通業(yè)發(fā)展水平存在顯著差異。

影響因素的選取與模型構(gòu)建

（一）影響因素的選取

商貿(mào)流通業(yè)是一種以批發(fā)和零售為主的綜合性服務(wù)產(chǎn)業(yè)，用社會(huì)消費(fèi)品零售總額作為衡量商貿(mào)流通業(yè)發(fā)展的指標(biāo)，社會(huì)消費(fèi)品零售額的高低反映了商貿(mào)流通業(yè)發(fā)展水平的高低。本文首先研究近十年京津冀13個(gè)城市的商貿(mào)流通業(yè)發(fā)展的時(shí)間差異特性及空間集聚效應(yīng)，同時(shí)采用2014-2018年13個(gè)城市的時(shí)空面板數(shù)據(jù)，研究商貿(mào)流通業(yè)發(fā)展受各因素影響的顯著程度，這13個(gè)城市分別為北京、天津、邯鄲、石家莊、邢臺(tái)、承德、唐山、保定、張家口、秦皇島、滄州、衡水、廊坊。商貿(mào)流通業(yè)發(fā)展受多方面因素影響，但總體來看，受社會(huì)經(jīng)濟(jì)發(fā)展水平和人文因素影響最大，所以本文選取GDP、人均GDP、人口數(shù)量、人口密度、進(jìn)出口貿(mào)易額、第一、第二和第三產(chǎn)業(yè)增加值等8個(gè)指標(biāo)，研究對商貿(mào)流通業(yè)發(fā)展影響的顯著程度（張曉倩等，2017）。

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

Global Moran's I指數(shù)模型的構(gòu)建。Global Moran's I指數(shù)（張娜，2019）是用來衡量空間自相關(guān)程度的綜合性指標(biāo)，該指標(biāo)用來反映區(qū)域?qū)傩灾档姆植际欠窦?、離散或者隨機(jī)分布，其值域?yàn)閇-1，1]，取值越接近1，則正相關(guān)性越強(qiáng)，取值越接近-1，則負(fù)相關(guān)性越強(qiáng)，取值為0表示隨機(jī)分布，同時(shí)用二元對稱空間權(quán)重距離矩陣W表示n個(gè)x區(qū)域的臨近距離關(guān)系，Global Moran's I指數(shù)及權(quán)重矩陣計(jì)算公式如下：

空間回歸模型的構(gòu)建。傳統(tǒng)OLS回歸方法基于變量空間相互獨(dú)立，當(dāng)變量之間存在空間自相關(guān)性，回歸殘差之間相互影響，傳統(tǒng)OLS回歸將失去效用，因此本文采用兩種模型研究空間相關(guān)性：一是基于因變量空間自相關(guān)的空間滯后模型（Spatial Lag Model，SLM）；二是基于回歸殘差項(xiàng)空間相關(guān)的空間誤差模型（Spatial Error Model，SEM）（郭付友等，2019）?？臻g滯后模型計(jì)算公式為：

式（3）中，y表示因變量矩陣；ρ和W分別為空間回歸系數(shù)和空間權(quán)重矩陣；Wy為空間滯后變量；X為解釋變量矩陣；β為待估參數(shù)；ε為隨機(jī)誤差項(xiàng)。對上述公式進(jìn)行變換可得：

從式（4）可以看出，每個(gè)i地區(qū)的y值不僅與這個(gè)地區(qū)的解釋變量有關(guān)，同時(shí)通過乘以一個(gè)空間因子(1-ρW)-1，也受其他地區(qū)解釋變量的影響。

空間誤差模型計(jì)算公式為：

其中空間殘差項(xiàng)μ的計(jì)算公式為：

從而公式（5）轉(zhuǎn)化為：

式中，λ表示空間殘差自回歸系數(shù)，上述公式表明每個(gè)i地區(qū)的y值受其他地區(qū)的隨機(jī)誤差ε影響，影響系數(shù)為 (1-λW)-1。

本文采用空間誤差模型和空間滯后模型中的固定效應(yīng)和隨機(jī)效應(yīng)研究各因素對商貿(mào)流通業(yè)發(fā)展影響的顯著程度，模型如下：

式中，ρ和λ分別表示空間自回歸和誤差自相關(guān)系數(shù)，wij表示i和j地區(qū)的空間權(quán)重矩陣，φ為空間自相關(guān)誤差項(xiàng)，μit表示i地區(qū)t時(shí)期的空間殘差項(xiàng)，ai和δt分別表示固定空間效應(yīng)和固定時(shí)間效應(yīng)。

模型檢驗(yàn)及結(jié)果分析

（一）商貿(mào)流通業(yè)時(shí)空分異特征分析

本文用京津冀社會(huì)消費(fèi)品零售總額來衡量京津冀地區(qū)商貿(mào)流通業(yè)的發(fā)展水平，利用Arcgis10.2軟件，分別計(jì)算2009-2018年商貿(mào)流通業(yè)發(fā)展的Global Moran's I指數(shù)，如表1所示，Global Moran's I指數(shù)均小于0，說明京津冀商貿(mào)流通業(yè)發(fā)展水平存在明顯的空間差異性，形成了發(fā)展水平高和發(fā)展水平低相互交錯(cuò)的空間格局。同時(shí)，近十年Global Moran's I指數(shù)波動(dòng)并不明顯，說明京津冀地區(qū)商貿(mào)流通業(yè)發(fā)展水平存在一定的空間集聚效應(yīng)。

圖1 2009-2018年北京、天津、河北三地商貿(mào)流通業(yè)發(fā)展規(guī)模變化

圖2 2018年京津冀13個(gè)城市商貿(mào)流通業(yè)發(fā)展水平對比

表1 京津冀地區(qū)2009-2018年商貿(mào)流通業(yè)發(fā)展Global Moran's I指數(shù)

圖3 2010年商貿(mào)流通冷熱點(diǎn)區(qū)域

圖4 2014年商貿(mào)流通冷熱點(diǎn)區(qū)域

圖5 2018年商貿(mào)流通冷熱點(diǎn)區(qū)域

為了深入了解京津冀地區(qū)商貿(mào)流通業(yè)發(fā)展的時(shí)空分異特征，利用Arcgis中的冷熱點(diǎn)區(qū)域分析（張?bào)憔甑龋?019），并結(jié)合自然斷裂點(diǎn)法，分別計(jì)算2010、2014、2018年的Getis-OrdGi*，Getis-OrdGi*可以區(qū)分出商貿(mào)流通業(yè)發(fā)展的冷熱點(diǎn)空間分布，并將商貿(mào)流通業(yè)發(fā)展分為冷熱區(qū)域4類，圖3、圖4、圖5分別為2010年、2014年和2018年冷熱點(diǎn)區(qū)域分布。

總體來看，京津冀商貿(mào)流通業(yè)發(fā)展在空間分布上呈現(xiàn)出明顯的分異特征，2010年、2014年、2018年呈現(xiàn)出中部向四周依次遞減的空間梯度特征，北京市始終是熱點(diǎn)區(qū)，張家口、承德、秦皇島等市始終是冷點(diǎn)區(qū)，區(qū)域差異顯著。此外，不同類型的城市數(shù)量隨著時(shí)間的推移也在發(fā)生改變，冷點(diǎn)區(qū)由2010年的8個(gè)減少為2018年的6個(gè)，石家莊和唐山也由2010年的次冷區(qū)，發(fā)展成為2018年的次熱區(qū)，邯鄲和滄州也由2010年的冷點(diǎn)區(qū)，演變成了2018年的次冷區(qū)，這說明京津冀地區(qū)的商貿(mào)流通業(yè)發(fā)展水平不斷提高，并呈現(xiàn)出空間相對穩(wěn)固的特征。從時(shí)空演化趨勢來看，商貿(mào)流通業(yè)發(fā)展呈現(xiàn)兩極分化態(tài)勢，以北京、天津、石家莊和唐山為主的城市群逐漸演化成了熱點(diǎn)和次熱點(diǎn)發(fā)展區(qū)域，張家口、承德和廊坊等市雖與北京相鄰，但商貿(mào)流通業(yè)發(fā)展始終處于冷點(diǎn)區(qū)域，北京對其周圍城市的商貿(mào)流通業(yè)發(fā)展產(chǎn)生了一定束縛。

為了進(jìn)一步研究京津冀地區(qū)商貿(mào)流通業(yè)發(fā)展的演化過程，利用自然斷裂點(diǎn)法計(jì)算2009-2013年和2014-2018年商貿(mào)流通業(yè)增長率的Getis-OrdGi*，并將其分為4類，如圖6和圖7所示。總體來看，增長次冷區(qū)域變化明顯，由張家口、承德、秦皇島、唐山和廊坊轉(zhuǎn)移到保定、石家莊和邯鄲，北京和天津也演變?yōu)榱嗽鲩L冷點(diǎn)區(qū)域，主要因?yàn)榫┙騼傻亟?jīng)濟(jì)發(fā)展水平和區(qū)位發(fā)展優(yōu)勢顯著，流通產(chǎn)業(yè)發(fā)展水平領(lǐng)先于京津冀其它城市。各地級市商貿(mào)流通業(yè)的快速增長，促進(jìn)了京津冀地區(qū)商貿(mào)流通業(yè)由差異化發(fā)展向區(qū)域協(xié)同發(fā)展轉(zhuǎn)變。

（二）商貿(mào)流通業(yè)影響因素分析

1.各影響因素總體相關(guān)性分析。商貿(mào)流通業(yè)發(fā)展受多種因素影響，主要分為經(jīng)濟(jì)、人文、地理三大方面因素，這些因素之間相互制約、相互作用，從而形成復(fù)雜的網(wǎng)絡(luò)系統(tǒng)，京津冀地區(qū)商貿(mào)流通業(yè)發(fā)展的各因素總體相關(guān)性（韓兆洲等，2019）如圖8所示，仍用社會(huì)消費(fèi)品零售額表示商貿(mào)流通業(yè)的發(fā)展水平，可以看出，商貿(mào)流通業(yè)發(fā)展與其他因素存在高度相關(guān)性，這表明GDP和人口數(shù)量的增長、進(jìn)出口貿(mào)易的增加、第二和第三產(chǎn)業(yè)的增長，均會(huì)促進(jìn)商貿(mào)流通業(yè)的發(fā)展，其中GDP和人口數(shù)量對商貿(mào)流通業(yè)發(fā)展的影響系數(shù)最高，其次為第三產(chǎn)業(yè)增加值。值得關(guān)注的是，第三產(chǎn)業(yè)增加值與GDP之間的影響系數(shù)最高，呈現(xiàn)出高度相關(guān)性，這表明包含交通運(yùn)輸業(yè)等服務(wù)性產(chǎn)業(yè)在內(nèi)的第三產(chǎn)業(yè)在社會(huì)經(jīng)濟(jì)水平較高的地區(qū)發(fā)展更好。

圖6 2009-2013年商貿(mào)流通增長率冷熱點(diǎn)區(qū)域

圖7 2014-2018年商貿(mào)流通增長率冷熱點(diǎn)區(qū)域

圖8 基于R語言的各因素總體相關(guān)性

2.各因素對商貿(mào)流通業(yè)影響顯著性分析。利用選取的8個(gè)影響因素，研究各因素影響商貿(mào)流通業(yè)發(fā)展的顯著程度，收集2014-2018年京津冀13個(gè)城市的影響因素?cái)?shù)值，帶入空間回歸模型進(jìn)行檢驗(yàn)，空間回歸分析結(jié)果如表2所示。空間滯后解釋變量ρ為-0.0128，空間誤差系數(shù)λ為0.0431，空間滯后模型通過了顯著性檢驗(yàn)，說明空間滯后效應(yīng)會(huì)對商貿(mào)流通業(yè)的發(fā)展產(chǎn)生影響（劉明，2018）。從空間回歸結(jié)果可知，空間滯后的系數(shù)為負(fù)值，說明一個(gè)地區(qū)的商貿(mào)流通業(yè)發(fā)展受周圍地區(qū)商貿(mào)流通業(yè)發(fā)展的負(fù)面影響，產(chǎn)業(yè)發(fā)展過程也是一個(gè)發(fā)展要素集聚與擴(kuò)散的過程，如果一個(gè)地區(qū)的經(jīng)濟(jì)、社會(huì)等產(chǎn)業(yè)發(fā)展要素在一個(gè)地區(qū)的集聚效應(yīng)大于擴(kuò)散效應(yīng)，那么勢必會(huì)影響周圍地區(qū)產(chǎn)業(yè)的發(fā)展，這一點(diǎn)也解釋了北京市商貿(mào)流通業(yè)發(fā)展遙遙領(lǐng)先于鄰近的張家口、承德和廊坊。

從各因素對商貿(mào)流通業(yè)發(fā)展的影響來看，GDP及人均GDP影響顯著，反映了人民生活水平不斷提高能夠刺激消費(fèi)，從而促進(jìn)商貿(mào)流通業(yè)的發(fā)展。第三產(chǎn)業(yè)及進(jìn)出口貿(mào)易額的快速發(fā)展同樣能促進(jìn)商貿(mào)流通業(yè)的發(fā)展，作為依賴交通運(yùn)輸業(yè)發(fā)展的商貿(mào)流通業(yè)，相比于第一產(chǎn)業(yè)和第二產(chǎn)業(yè)，受第三產(chǎn)業(yè)類型的影響較大，而京津冀地區(qū)的產(chǎn)業(yè)發(fā)展也主要依賴第三產(chǎn)業(yè)。北京產(chǎn)業(yè)高端化趨勢明顯，天津的服務(wù)業(yè)也在快速提升，而河北的制造業(yè)發(fā)展則占據(jù)較大市場份額，京津冀三地的產(chǎn)業(yè)結(jié)構(gòu)不同，也造成了其商貿(mào)流通業(yè)發(fā)展的空間差異性，而京津冀協(xié)同發(fā)展則要求京津冀各城市加速產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，實(shí)現(xiàn)京津冀區(qū)域內(nèi)產(chǎn)業(yè)協(xié)同程度不斷提升，從而同步促進(jìn)各地區(qū)商貿(mào)流通業(yè)協(xié)同發(fā)展。

結(jié)論與建議

（一）結(jié)論

本文利用空間統(tǒng)計(jì)分析方法，對京津冀商貿(mào)流通業(yè)發(fā)展的時(shí)空分異特征進(jìn)行了分析，并在此基礎(chǔ)上，研究了經(jīng)濟(jì)、人文和產(chǎn)業(yè)等因素對商貿(mào)流通業(yè)發(fā)展影響的顯著程度，得出以下結(jié)論：第一，京津冀地區(qū)商貿(mào)流通業(yè)的發(fā)展時(shí)空差異性顯著，總體呈現(xiàn)出由北京、天津和石家莊等城市為中心，依次向四周衰減的趨勢，但增長率的空間特征則呈相反的趨勢。第二，北京對臨近城市的輻射作用顯著，削弱這種發(fā)展空間上的顯著差別，是實(shí)現(xiàn)商貿(mào)流通業(yè)整體向上發(fā)展和協(xié)同發(fā)展的有效途徑。第三，GDP和產(chǎn)業(yè)結(jié)構(gòu)等因素均會(huì)影響商貿(mào)流通業(yè)的發(fā)展，特別是第三產(chǎn)業(yè)的發(fā)展，與商貿(mào)流通業(yè)發(fā)展呈正相關(guān)，調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，優(yōu)化產(chǎn)業(yè)資源配置能有效促進(jìn)京津冀地區(qū)商貿(mào)流通業(yè)快速協(xié)同發(fā)展。

表2 空間回歸分析結(jié)果

（二）建議

順應(yīng)京津冀協(xié)同發(fā)展這一國家戰(zhàn)略，努力促進(jìn)京津冀各城市商貿(mào)流通業(yè)協(xié)同發(fā)展，縮小北京對鄰近城市商貿(mào)流通業(yè)發(fā)展的束縛作用。張家口、承德和廊坊等城市的商貿(mào)流通業(yè)發(fā)展水平雖然較低，但近幾年增長率較快，發(fā)展勢頭迅猛，各市政府應(yīng)繼續(xù)支持本地商貿(mào)流通業(yè)發(fā)展，在政策扶持及資金補(bǔ)貼方面向批發(fā)和零售業(yè)傾斜，保持流通產(chǎn)業(yè)快速和平穩(wěn)發(fā)展。

統(tǒng)籌協(xié)調(diào)商貿(mào)流通業(yè)與社會(huì)經(jīng)濟(jì)發(fā)展關(guān)系，優(yōu)化京津冀商品流通網(wǎng)絡(luò)結(jié)構(gòu)，完善商品流通體制。研究發(fā)現(xiàn)GDP對商貿(mào)流通業(yè)的發(fā)展有著顯著影響，社會(huì)經(jīng)濟(jì)的快速發(fā)展也能促進(jìn)商貿(mào)流通業(yè)的快速發(fā)展，在京津冀區(qū)域經(jīng)濟(jì)快速發(fā)展的社會(huì)背景下，優(yōu)化商貿(mào)流通網(wǎng)絡(luò)，完善商貿(mào)流通環(huán)節(jié)，深化商貿(mào)流通體制改革，更能適應(yīng)經(jīng)濟(jì)發(fā)展新常態(tài)，實(shí)現(xiàn)商貿(mào)流通業(yè)與區(qū)域經(jīng)濟(jì)協(xié)調(diào)快速發(fā)展。

加速產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和資源配置，制造適合商貿(mào)流通業(yè)發(fā)展的產(chǎn)業(yè)環(huán)境。京津冀三地產(chǎn)業(yè)結(jié)構(gòu)差異明顯，北京主要發(fā)展高端服務(wù)業(yè)，商貿(mào)流通業(yè)已成為北京的第一產(chǎn)業(yè)，天津近幾年經(jīng)濟(jì)增長緩慢，產(chǎn)業(yè)結(jié)構(gòu)嚴(yán)重失衡，傳統(tǒng)重工業(yè)仍為天津主要產(chǎn)業(yè)，而河北各市的商貿(mào)流通業(yè)發(fā)展相對滯后，對商貿(mào)流通業(yè)影響較為顯著的第三產(chǎn)業(yè)發(fā)展相對落后，交通業(yè)和服務(wù)業(yè)發(fā)展緩慢，因此，加速京津冀地區(qū)的產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，重視對第三產(chǎn)業(yè)的發(fā)展，優(yōu)化產(chǎn)業(yè)資源配置，是促進(jìn)京津冀商貿(mào)流通業(yè)協(xié)調(diào)快速發(fā)展的有效手段。