朱文娟 曾 豪 鄧 彬

1.湖南工業(yè)大學商學院 湖南 株洲 412007

2. 湖南工業(yè)大學湖南省包裝經(jīng)濟研究基地 湖南 株洲 412007

3. 湖南工業(yè)大學計算機學院 湖南 株洲 412007

1 研究背景

包裝印刷是以各種包裝材料為載體的印刷，是包裝工業(yè)和現(xiàn)代產(chǎn)業(yè)體系的重要組成部分?！丁笆奈濉毖h(huán)經(jīng)濟發(fā)展規(guī)劃》和《“十四五”現(xiàn)代物流發(fā)展規(guī)劃》中明確指出，要推行重點產(chǎn)品綠色設(shè)計，推動包裝和包裝印刷減量化；提升與優(yōu)化印刷書寫紙、包裝用紙等品質(zhì)，適應(yīng)多元化、個性化消費市場的需求。由此可見，包裝印刷行業(yè)對構(gòu)建新發(fā)展格局、實施綠色發(fā)展戰(zhàn)略、促進消費市場發(fā)展都具有重要作用。

近十年來，中國包裝印刷行業(yè)發(fā)展規(guī)模持續(xù)擴大，百強企業(yè)規(guī)模不斷增長[1]（如圖1所示），但包裝印刷百強企業(yè)的分布具有較大的地區(qū)差異。目前，空間物流布局與演變已成為學者們的研究熱點，主要體現(xiàn)在以下三個方面：一是研究內(nèi)容，主要涉及物流企業(yè)空間演化特征及影響因素[2-4]、物流企業(yè)空間差異[5-6]、物流企業(yè)空間聯(lián)系及網(wǎng)絡(luò)結(jié)構(gòu)[7-9]，以及物流企業(yè)空間溢出效應(yīng)等[10-12]。二是研究對象，主要聚焦于特定區(qū)域內(nèi)物流企業(yè)[13-16]，以及某一特定功能類型的物流企業(yè)或某一產(chǎn)業(yè)與物流業(yè)融合發(fā)展等[17-18]，而針對包裝印刷企業(yè)的研究比較缺乏。三是研究方法，越來越多的學者借助探索性空間數(shù)據(jù)分析（exploring spatial data analysis，ESDA）技術(shù)[4,10,14]以及依托地理信息系統(tǒng)（geographic information system，GIS）的核密度估計方法[3,19]，開展與空間數(shù)據(jù)有關(guān)的研究。包裝印刷企業(yè)是現(xiàn)代物流產(chǎn)業(yè)體系中不可或缺的一部分，而在全國范圍內(nèi)針對包裝印刷企業(yè)發(fā)展時空演化特征的研究仍比較缺乏。包裝印刷百強企業(yè)作為行業(yè)龍頭，對產(chǎn)業(yè)鏈發(fā)展具有集聚和輻射效應(yīng)。深入研究其空間分布特征、演變趨勢及影響因素，對于進一步優(yōu)化區(qū)域布局、推動包裝產(chǎn)業(yè)實現(xiàn)更高質(zhì)量、更加廣闊的發(fā)展有一定的價值。

本研究以《印刷經(jīng)理人》界定的中國排名前100家包裝印刷企業(yè)為研究對象，以省域為基本單位，選取2012, 2017, 2022年3個時間點，運用標準差橢圓、核密度估計、熱點分析等GIS空間分析方法，探究中國包裝印刷百強企業(yè)時空演變特征，同時采用地理探測器（GeoDetector）方法對影響因素及其交互作用進行定量分析。研究結(jié)果以期能為包裝印刷產(chǎn)業(yè)的發(fā)展及合理布局提供參考。

圖1 2010—2020年百強企業(yè)平均銷售收入及變化率Fig. 1 Average sales revenue and change rate of top 100 enterprises from 2010 to 2020

2 數(shù)據(jù)來源和研究方法

2.1 數(shù)據(jù)來源和處理

以全國包裝印刷企業(yè)100強為研究對象，按照包裝印刷企業(yè)注冊地或總部所在地，利用Google Earth獲取其地理坐標數(shù)據(jù)，借助ArcGIS10.8軟件平臺對包裝印刷企業(yè)進行坐標匹配，并實現(xiàn)空間可視化。研究選取2012, 2017, 2022年數(shù)據(jù)，分析包裝印刷百強企業(yè)發(fā)展的時空演變特征及其影響因素。相關(guān)統(tǒng)計數(shù)據(jù)來源于《中國統(tǒng)計年鑒》。

2.2 研究方法與步驟

2.2.1 標準差橢圓法

標準差橢圓法以可視化的方式，從中心性、展布范圍、分布方向、分布形狀等角度描述研究對象的空間分布特征[20]，其主要參數(shù)包括：中心點經(jīng)緯度、短軸、長軸和旋轉(zhuǎn)角θ。本研究采用標準差橢圓對不同年份包裝印刷百強企業(yè)進行可視化操作。其詳細計算公式參照文獻[21]。

2.2.2 核密度估計法

核密度估計法是利用連續(xù)的密度曲線描述隨機變量空間分布形態(tài)的一種非參數(shù)方法，可用于揭示地理事物分布局部特征[22-23]。本研究將包裝印刷百強企業(yè)的地理坐標視作“點”要素，采用核密度估計對中國包裝印刷百強企業(yè)進行密度制圖與分析，其計算公式[24]為

式中：f(y)為核密度估計值；n為要素總數(shù)100；y-yi為估值點y到要素yi的距離；l為帶寬200 km。

f(y)值越大，表明包裝印刷百強企業(yè)分布越集中。

2.2.3 熱點分析法

熱點分析法是通過尋找高值或低值元素的空間集群，采用計算Getis-Ord統(tǒng)計量分析局部地區(qū)空間集聚情況，其計算公式[25]為

式中：xc為空間單元c的屬性值；為數(shù)據(jù)集中要素屬性的平均值；S為各要素標準差，wbc為空間單元b和c之間的空間權(quán)重；a為空間單元總數(shù)。

為區(qū)分樣本點空間分布的冷熱點地區(qū)，本研究將值進行標準化處理，其計算公式為

當?shù)挠^察值大于期望值并且具有統(tǒng)計學意義時，該區(qū)域為“熱點區(qū)”；當觀察值小于期望值時，該區(qū)域為“冷點區(qū)”[26]。

2.2.4 地理探測器法

地理探測器法主要用于探究不同X因子之間的相互作用，即因子X1和X2共同作用時是否會增強或減弱對因變量Y的影響力。探測因子特征的解釋力用q值度量，并用式（4）計算。兩個因子之間的關(guān)系如表1[27]所示。

表1 雙因子交互作用的類型Table 1 Types of two-factor interactions

式中：h為自變量X的分類數(shù)，h=1, 2, …,L；N、Nh分別為全區(qū)和類h的單元數(shù)；σ2、分別為全區(qū)和類h的Y值方差。

q的取值范圍為[0, 1]，表示某因子X解釋了100q%的Y。q值越大，即自變量X對Y的影響越強，反之則越弱。

3 中國包裝印刷百強企業(yè)時空演變特征分析

3.1 包裝印刷百強企業(yè)空間分布特征

采用ArcGIS10.8軟件中的方向分布工具對2012,2017和2022年中國包裝印刷百強企業(yè)進行標準差橢圓分析，并實現(xiàn)可視化。橢圓分布和相關(guān)參數(shù)分別如圖2和表2所示。本文地圖基于自然資源部地圖技術(shù)審查中心標準地圖服務(wù)網(wǎng)站下載的審圖號為GS(2023)2767號的標準地圖制作。

由圖2和表2可知，2012—2022年中國包裝印刷百強企業(yè)的空間分布與“胡煥庸線”基本保持一致，空間集聚格局均表現(xiàn)出不同規(guī)則的環(huán)狀分布特征。十年間，百強企業(yè)空間分布格局由“東北-西南”走向轉(zhuǎn)變?yōu)椤拔鞅?東南”走向，以標準橢圓的中軸線為基準，總體呈現(xiàn)向中西部地區(qū)偏移的趨勢。

圖2 百強企業(yè)空間分布及標準差橢圓Fig. 2 The spatial distribution and standard deviation ellipsoid of top 100 enterprises

彩圖

3.2 包裝印刷企業(yè)空間集聚變化

為進一步分析中國包裝印刷百強企業(yè)空間集聚變化特征，采用ArcGIS10.8軟件空間分析中集成的核密度工具，分別生成2012, 2017和2022年中國包裝印刷百強企業(yè)的核密度分布圖，如圖3所示。

由圖3可知，2012年時包裝印刷百強企業(yè)在“胡煥庸線”以東的沿海地區(qū)出現(xiàn)了明顯的空間集聚特征。2017—2022年，包裝印刷百強企業(yè)逐漸向“胡煥庸線”西側(cè)的內(nèi)陸地區(qū)擴散。從核密度分布來看，百強企業(yè)不僅在沿海地區(qū)集聚程度有所加強，在以湖北、陜西、四川等為主的內(nèi)陸省份集聚程度也均有所提高。

圖3 百強企業(yè)空間分布核密度分析圖Fig. 3 Kernel density analysis map of the spatial distribution of top 100 enterprises

彩圖

3.3 包裝印刷企業(yè)空間集聚程度

運用熱點分析工具，分別計算2012, 2017, 2022年的Getis-Ord統(tǒng)計數(shù)據(jù)，并采用“自然斷裂點法”將局部的統(tǒng)計量分類，以進一步探究中國包裝印刷百強企業(yè)的空間集聚程度情況，結(jié)果如圖4所示。

由圖4可以看出，2012—2022年，中國包裝印刷百強企業(yè)從以沿海地區(qū)為主逐漸向長江經(jīng)濟帶發(fā)展延伸，形成了海岸與長江經(jīng)濟帶兩個一級重點經(jīng)濟帶的“T字型空間戰(zhàn)略”①“T字型空間戰(zhàn)略”由中國科學院院士陸大道在1984年提出，即海岸經(jīng)濟帶和長江經(jīng)濟帶兩個一級重點經(jīng)濟帶形成“T型”，并在長三角交匯，長江經(jīng)濟帶將成渝地區(qū)、武漢地區(qū)與海岸經(jīng)濟帶聯(lián)系起來。分布。在冷熱分布上，十年間包裝印刷百強企業(yè)冷熱點分布格局較為穩(wěn)定，熱點與次熱點區(qū)中河南、安徽、湖北、浙江、廣東和山東等始終保持不變，冷點地區(qū)均處于西北部以及東北部。

圖4 百強企業(yè)空間分布熱點分析圖Fig. 4 Hot spot analysis of spatial distribution of top 100 enterprises

彩圖

4 中國包裝印刷百強企業(yè)空間分布影響因素

4.1 影響因素指標的選取及相關(guān)性分析

4.1.1 指標選取

包裝印刷百強企業(yè)的空間分布受多種因素的影響。在借鑒以往研究[4,25,28-29]的基礎(chǔ)上，基于中國包裝印刷百強企業(yè)空間布局變化特征以及數(shù)據(jù)的可獲得性，選取了區(qū)域經(jīng)濟實力（X1）、居民消費能力（X2）、區(qū)域政策（X3）、工業(yè)化水平（X4）、城鎮(zhèn)化水平（X5）、市場規(guī)模（X6）、勞動力水平（X7）、國際貿(mào)易水平（X8）8個指標，利用地理探測器模型，探究中國包裝印刷百強企業(yè)空間分異的驅(qū)動因素及其演化過程。指標的選取和計算方法如表3所示。

表3 影響因素指標及其計算方法Table 3 Indicators of impact factors and their calculation

4.1.2 相關(guān)性分析

為保證選取指標的合理性，采用SPSS軟件對各年份省域包裝印刷百強企業(yè)數(shù)量與選取指標進行Pearson相關(guān)性分析，結(jié)果如表4所示。

表4 百強企業(yè)數(shù)量與選取指標的Pearson相關(guān)系數(shù)Table 4 Pearson’s correlation coefficient between the number of top 100 companies and selected indicators

由表4可知，包裝印刷百強企業(yè)Pearson相關(guān)系數(shù)均為正數(shù)，表明選取的8個指標與包裝印刷百強企業(yè)數(shù)量之間均存在正向相關(guān)關(guān)系，從側(cè)面反映出指標選取具有一定的合理性。具體來看，區(qū)域經(jīng)濟實力、居民消費能力、區(qū)域政策、城鎮(zhèn)化水平、市場規(guī)模 、勞動力水平、國際貿(mào)易水平均在3個時期內(nèi)通過了顯著性檢驗，Pearson相關(guān)系數(shù)相對較高。

4.2 影響因素分析

4.2.1 主要影響因素及其影響力分析

采用地理探測器分析中國包裝印刷百強企業(yè)空間分異格局的影響因素。首先，借助ArcGIS10.8軟件中的轉(zhuǎn)換工具，對8個自變量（即影響因素）進行數(shù)據(jù)柵格化處理；然后利用Spatial Analyst工具進行因子重分類，采用自然斷裂點法將影響因素分為4類，再利用數(shù)據(jù)管理工具建立采樣點；最后按照采樣點提取分析，將所得數(shù)據(jù)讀入地理探測器軟件，分別計算2012、2017、2022年影響因素對包裝印刷百強企業(yè)的影響力q值，如表5所示。

表5 不同影響因素對百強企業(yè)分布的影響力q值Table 5 q value of influence of different influencing factors on the distribution of top 100 enterprises

由表5可知，國際貿(mào)易水平、區(qū)域經(jīng)濟實力、勞動力水平，是影響近十年間中國包裝印刷百強企業(yè)空間分布格局的三大主要因素（q值均大于0.5）。其中國際貿(mào)易水平對2012年和2017年百強企業(yè)空間分布的影響最高；2022年，區(qū)域經(jīng)濟實力取代國際貿(mào)易水平成為影響百強企業(yè)空間分布格局的首要因素。居民消費能力、工業(yè)化水平和城鎮(zhèn)化水平的q值相對較小，是包裝印刷百強企業(yè)空間分異的次要影響因素。2012—2022年，除市場規(guī)模的q值呈上升趨勢外，其余影響因素的q值均呈下降趨勢，且各影響因素之間的q值差異也在縮小。這說明中國包裝印刷百強企業(yè)的空間分布并非由單一因素決定，而是取決于多個因素，且各因素的影響程度趨于均衡。

4.2.2 主要影響因素的作用機制分析

從各影響因素的q值看，國際貿(mào)易水平、區(qū)域經(jīng)濟實力、勞動力水平、市場規(guī)模和區(qū)域政策是影響中國包裝印刷百強企業(yè)空間分布格局的五大核心因素。

1）國際貿(mào)易水平。2012—2022年期間國際貿(mào)易水平q值的均值為0.702 7，是影響中國包裝印刷百強企業(yè)空間分布的第一大因素。中國是世界第二包裝大國，包裝印刷作為包裝產(chǎn)業(yè)的重要組成部分，一直深度參與全球產(chǎn)業(yè)鏈、價值鏈的分工合作，進出口貿(mào)易額占包裝印刷行業(yè)總產(chǎn)值的比重較大。由于海運的經(jīng)濟性和便捷性，大量包裝印刷企業(yè)在江蘇、浙江、上海、廣東等沿海地區(qū)集聚。此外，沿海地區(qū)也擁有較多的外資企業(yè)，能夠帶來國際資源與合作機會，這進一步促進了包裝印刷產(chǎn)業(yè)鏈的集聚和發(fā)展。同時也看到，十年間國際貿(mào)易水平q值出現(xiàn)了下降趨勢，這說明隨著國內(nèi)市場規(guī)模的不斷擴大以及“雙循環(huán)”新發(fā)展格局的加快構(gòu)建，國際貿(mào)易水平對百強企業(yè)空間布局的影響力度逐漸下降。

2）區(qū)域經(jīng)濟實力。2012—2022年期間區(qū)域經(jīng)濟實力，對中國包裝印刷百強企業(yè)空間分布格局的影響力始終處于較高的水平，q值的均值為0.600 1。從包裝印刷產(chǎn)業(yè)鏈來看，上游連接紙、塑料、玻璃、金屬、竹木等包裝材料，以及油墨和印刷裝備等，下游銜接食品、藥品、電子、日用品等行業(yè)，貫通設(shè)計、生產(chǎn)、流通、消費等各環(huán)節(jié)，對國民經(jīng)濟高效運行起到重要的支撐作用。包裝印刷百強企業(yè)更加傾向于在區(qū)域經(jīng)濟實力強、產(chǎn)業(yè)覆蓋廣、產(chǎn)業(yè)鏈更齊備的長三角、珠三角等沿海、沿江地區(qū)集聚，并隨著產(chǎn)業(yè)的梯度轉(zhuǎn)移逐漸向內(nèi)陸擴散。隨著我國區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略的深入推進，內(nèi)陸地區(qū)的經(jīng)濟取得較大發(fā)展，因此十年間區(qū)域經(jīng)濟實力的q值有所下降，其對百強企業(yè)空間分布的影響力度也略有減小，但是從整體看依然是主要的影響因素之一。

3）勞動力水平。勞動力水平q值的均值為0.572 3，是影響中國包裝百強企業(yè)空間布局的第三大因素?？傮w而言，包裝印刷屬于勞動密集型產(chǎn)業(yè)，地區(qū)勞動力水平在較大程度上影響了包裝印刷百強企業(yè)的空間分布格局。這也是江蘇、浙江、上海、廣東等沿海地區(qū)作為大量外來勞動力的流入地，以及河南、安徽、湖北、四川等人口大省集聚了大量包裝印刷企業(yè)的重要原因。隨著包裝印刷行業(yè)信息化、自動化和智能化程度的提升，數(shù)字印刷、自動化倉庫、自動機器人、自動檢驗機等技術(shù)和設(shè)備正在逐步應(yīng)用，對于勞動力的需求相應(yīng)減少，這也是勞動力水平q值有所下降的重要原因之一。

4）市場規(guī)模。十年間，市場規(guī)模q值從0.436 7上升到0.533 9，不僅是影響包裝印刷百強企業(yè)空間分布的第四大因素，也是唯一具有增強趨勢的影響因素。中國包裝行業(yè)市場規(guī)模超過萬億，并將長期具有強勁增長勢頭。市場規(guī)?？梢詭懋a(chǎn)業(yè)集聚，實現(xiàn)專業(yè)化分工與合作，降低各種交易成本，形成協(xié)同效應(yīng)和規(guī)模效應(yīng)，助力包裝印刷企業(yè)高質(zhì)量發(fā)展。隨著全國統(tǒng)一大市場建設(shè)，市場規(guī)模因素對百強企業(yè)空間布局的影響力還將不斷凸顯。

5）區(qū)域政策。區(qū)域政策是中國包裝印刷百強企業(yè)空間分異的第五大影響因素。區(qū)域政策越好的地方，越有利于吸引包括包裝印刷在內(nèi)的產(chǎn)業(yè)集聚并形成規(guī)模效應(yīng)。這十年是我國深入實施區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略的關(guān)鍵時期，除東部地區(qū)領(lǐng)先發(fā)展外，中西部地區(qū)的經(jīng)濟發(fā)展也得到了大幅度提升，為相關(guān)政策紅利的釋放提供了堅實基礎(chǔ)。各地政府加大招商引資和項目落地，通過優(yōu)化營商環(huán)境、實施稅費減免、提供貸款擔保、加強配套保障等方式，大力促進相關(guān)產(chǎn)業(yè)發(fā)展，在廣東、浙江、福建、四川等地涌現(xiàn)出一批具有規(guī)?；瘍?yōu)勢的包裝印刷產(chǎn)業(yè)園和產(chǎn)業(yè)基地。

4.3 影響因素交互作用分析

進一步使用地理探測器軟件中的交互作用探測器，確定兩個影響因素在交互作用下對包裝印刷百強企業(yè)空間分布的影響程度，結(jié)果如表6所示。

表6 百強企業(yè)不同影響因素的交互探測結(jié)果Table 6 Interactive detection results of different influencing factors of top 100 enterprises

由表6知，兩個影響因素之間交互作用后，主要通過非線性增強（NE）和雙因子增強（BE）兩種方式，增強了對中國包裝印刷百強企業(yè)空間分異的解釋力，其中雙因子增強型的結(jié)果最為常見。這表明，不同影響因素共同作用對包裝印刷百強企業(yè)空間分異的解釋力大于各單一因素的解釋力，且不存在相互獨立的影響因素。

5 結(jié)語

本研究基于ArcGIS軟件平臺，采用標準差橢圓分析、核密度分析和熱點分析，探究了2012—2022年中國包裝印刷百強企業(yè)時空演變特征，并利用地理探測器軟件研究其空間格局演變的影響因素及其交互作用，可得以下主要結(jié)論：

1）十年間，包裝印刷百強企業(yè)在“胡煥庸線”以東出現(xiàn)明顯的空間集聚特征?？臻g集聚格局均表現(xiàn)出不同規(guī)則的環(huán)狀分布特征，其中長江三角洲與珠江三角洲是主要分布地。

2）中國包裝印刷百強企業(yè)空間分布格局由“東北-西南”走向轉(zhuǎn)變?yōu)椤拔鞅?東南”走向，在以湖北、江西、安徽、陜西為主的內(nèi)陸地區(qū)的集聚程度也均有加強，并呈現(xiàn)出由沿海向內(nèi)陸擴散的趨勢。

3）十年間包裝印刷百強企業(yè)冷熱點分布格局較為穩(wěn)定，熱點與次熱點區(qū)中河南、安徽、湖北、浙江、廣東和山東仍始終保持不變，冷點地區(qū)均處于西北部及東北部。

4）國際貿(mào)易水平、區(qū)域經(jīng)濟實力、勞動力水平、市場規(guī)模和區(qū)域政策是影響中國包裝印刷百強企業(yè)空間分布格局的五大核心因素。其中，市場規(guī)模因素對包裝印刷百強企業(yè)空間分布影響呈上升趨勢。

5）中國包裝印刷百強企業(yè)的空間結(jié)構(gòu)并非由單一因素決定，而是取決于多個因素，且各因素的影響程度趨于均衡。

包裝印刷是現(xiàn)代產(chǎn)業(yè)體系的重要組成部分，作為中堅力量的百強企業(yè)的合理布局與健康發(fā)展，對促進中國包裝行業(yè)的高質(zhì)量發(fā)展和國民經(jīng)濟的高效運行具有重要的戰(zhàn)略意義?？傮w而言，近十年來包裝印刷百強企業(yè)呈現(xiàn)出的空間集聚特點和時空演變特征既反映了我國經(jīng)濟、社會和產(chǎn)業(yè)發(fā)展水平的地區(qū)差異性，也在一定程度上體現(xiàn)了我國區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略與產(chǎn)業(yè)轉(zhuǎn)移政策的實施成效。按照《中國包裝工業(yè)發(fā)展規(guī)劃（2021—2025）》[30]部署，要繼續(xù)推動包裝印刷企業(yè)區(qū)域布局優(yōu)化：一方面，要利用沿海地區(qū)包裝印刷產(chǎn)業(yè)的先發(fā)優(yōu)勢，充分發(fā)揮百強企業(yè)的龍頭作用，打造集群發(fā)展示范區(qū)和創(chuàng)新發(fā)展引領(lǐng)區(qū)；另一方面，也要繼續(xù)完善內(nèi)陸地區(qū)包裝印刷產(chǎn)業(yè)鏈，加快形成產(chǎn)業(yè)新興集聚區(qū)和產(chǎn)業(yè)轉(zhuǎn)移承接區(qū)。同時，要密切關(guān)注影響百強企業(yè)的空間集聚程度的重要因素及其交互作用，不斷優(yōu)化發(fā)展環(huán)境，厚植發(fā)展優(yōu)勢，促進包裝產(chǎn)業(yè)高質(zhì)量發(fā)展和包裝強國建設(shè)。