我國茶銷售市場的空間布局與影響因素研究

宋 娜，何亞麗，毛婭琪

燕山大學(xué) 經(jīng)濟(jì)管理學(xué)院，河北 秦皇島 066004

我國是世界上茶生產(chǎn)和消費(fèi)大國，茶葉的種植、生產(chǎn)與消費(fèi)文化歷史悠久。隨著現(xiàn)代經(jīng)濟(jì)和休閑產(chǎn)業(yè)的發(fā)展，人們對茶產(chǎn)品的需求已經(jīng)成為生活慣常，其需求偏好和消費(fèi)行為不斷呈現(xiàn)出新的變化和趨勢。從地理空間范圍研究我國各省和主要城市的茶銷售的分布格局及影響因素，有助于全面了解國內(nèi)茶市場特征，有效銜接我國茶產(chǎn)業(yè)的供給與需求市場，促進(jìn)茶產(chǎn)業(yè)的健康發(fā)展。

我國茶產(chǎn)業(yè)是涵蓋種植、生產(chǎn)、銷售的完整產(chǎn)業(yè)鏈條[1]，學(xué)者們分別從產(chǎn)業(yè)發(fā)展水平[2]、產(chǎn)業(yè)競爭力[3]、產(chǎn)業(yè)融合[4]、茶產(chǎn)業(yè)的空間分布格局[5]等不同角度對茶產(chǎn)業(yè)展開了研究。馬永強(qiáng) 等認(rèn)為生產(chǎn)者、消費(fèi)需求、產(chǎn)業(yè)金融平臺(tái)服務(wù)等構(gòu)成茶產(chǎn)業(yè)生態(tài)圈，通過三者的協(xié)同促進(jìn)茶產(chǎn)業(yè)生態(tài)圈的發(fā)展[6]。可見，茶產(chǎn)業(yè)的健康發(fā)展需要整個(gè)產(chǎn)業(yè)鏈的協(xié)同互促。從當(dāng)前的研究成果來看，學(xué)者們多側(cè)重于對茶產(chǎn)業(yè)鏈前端的研究，尤其是在空間布局的相關(guān)研究方面幾乎全部集中于茶產(chǎn)業(yè)資源種植和茶產(chǎn)業(yè)的生產(chǎn)空間分布格局。如趙芳利用GIS 軟件對銅陵市白茶的種植區(qū)域進(jìn)行評(píng)級(jí)和規(guī)劃[7]；張莉等以貴州中部地區(qū)茶園為研究范圍研究了茶棍薊馬的空間分布[8]；杜廣祖等以滇西南的高山茶園為例分析了茶黃薊馬的空間格局[9]；黃修杰等總結(jié)了25 年來廣東省茶葉的生產(chǎn)格局特征與集聚效應(yīng)[10]。此外，吳芹瑤等對近70 年的茉莉花茶的生產(chǎn)空間進(jìn)行了梳理，并分析了驅(qū)動(dòng)因素[11]；羅富民證明了18 個(gè)產(chǎn)茶省的茶產(chǎn)業(yè)生產(chǎn)空間集聚有助于農(nóng)民增收[12]。關(guān)于茶銷售市場的研究，一部分學(xué)者關(guān)注于出口市場，如張成鳳以親茶度和競爭態(tài)兩個(gè)指標(biāo)將茶葉出口市場分為四種類型[13]；王雨認(rèn)為可以把跨境電商作為出口平臺(tái)提高茶產(chǎn)品出口率。國內(nèi)消費(fèi)市場的研究側(cè)重關(guān)注于消費(fèi)者的個(gè)體行為[14]，石琳用調(diào)查問卷的形式對新式茶飲的消費(fèi)者進(jìn)行了人群畫 像[15]，張曉琳等對都勻毛尖茶的消費(fèi)群體特征進(jìn)行了深入研究[16]，張?jiān)吕驅(qū)Π布撞鑵^(qū)域品牌的消費(fèi)認(rèn)知結(jié)構(gòu)進(jìn)行了研究分析[17]；另外王東輝等認(rèn)識(shí)到電商市場的重要性，分析了電商市場中茶葉消費(fèi)特征，提出影響茶葉消費(fèi)的因素有可支配收入、禮品與社交、茶文化以及飲食習(xí)慣等[18]。

綜上可見，我國目前對茶產(chǎn)業(yè)的研究主要側(cè)重于產(chǎn)業(yè)鏈的種植和生產(chǎn)環(huán)節(jié)，與消費(fèi)需求市場銜接的銷售環(huán)節(jié)的研究相對較弱。而對銷售市場的研究一般以人員訪談、市場問卷調(diào)查及官方統(tǒng)計(jì)數(shù)據(jù)等為主，以地理空間屬性的大數(shù)據(jù)為基礎(chǔ)的研究還較為少見。本文以全國茶銷售市場作為對象，從宏觀視角把握當(dāng)前茶銷售的空間分布類型、分布均衡性和分布密度等，并對其空間特征形成的影響因素進(jìn)行深入分析，從作用機(jī)制的視角進(jìn)行市場調(diào)控，有助于我國茶產(chǎn)業(yè)的供需平衡發(fā)展。

1 數(shù)據(jù)來源及研究方法

1.1 數(shù)據(jù)來源

茶銷售市場非常龐大，由于電子商務(wù)平臺(tái)的興起，人們購茶的途徑也逐漸多樣化，但是目前線下實(shí)體消費(fèi)仍然是茶銷售市場的主體。為了數(shù)據(jù)統(tǒng)計(jì)的規(guī)范性和科學(xué)性，采用爬蟲技術(shù)在高德地圖餐飲服務(wù)類和購物服務(wù)類POI 數(shù)據(jù)中獲取了除港澳臺(tái)外全國31 個(gè)省份346 個(gè)城市的茶類經(jīng)營點(diǎn)信息，涵蓋茶館、茶莊、茶樓、茶店等不同稱呼的經(jīng)營店鋪。31 個(gè)省份的茶類經(jīng)營點(diǎn)分布如圖1 所示。

以全國346 個(gè)地級(jí)市為單位進(jìn)行茶類經(jīng)營點(diǎn)統(tǒng)計(jì)，并將結(jié)果按自然間斷點(diǎn)分析法進(jìn)行5級(jí)分類，在地圖中標(biāo)出，如圖2 所示。

對獲取的數(shù)據(jù)進(jìn)行選擇和整理，生成表格，再利用Arc GIS 10.2 軟件對所得茶類經(jīng)營點(diǎn)數(shù) 據(jù)進(jìn)行數(shù)字化處理，得到茶類經(jīng)營點(diǎn)的空間分布圖，并通過最鄰近點(diǎn)指數(shù)、地理集中指數(shù)和核密度分析方法，從空間分布類型、均衡性和分布密度3 個(gè)方面來分析茶銷售市場空間分布特征。

1.2 研究方法

1.2.1 空間特征分析

1.2.1.1 茶類經(jīng)營點(diǎn)分布類型分析

最鄰近指數(shù)是能夠體現(xiàn)點(diǎn)狀要素在空間中的相互鄰近程度的地理指標(biāo)，測算集聚情況，呈現(xiàn)研究對象隨機(jī)、均勻或聚集分布狀態(tài)。其中點(diǎn)狀要素最鄰近指數(shù)最大的為均勻分布，隨機(jī)分布次之，聚集分布最小。計(jì)算公式為：

圖1 茶類經(jīng)營點(diǎn)省份空間分布圖Figure 1 Spatial distribution map of tea business sites in provinces

圖2 茶類經(jīng)營點(diǎn)城市空間分布圖Figure 2 Urban spatial distribution map of tea business sites

R 為最鄰近指數(shù)，r1是任意一點(diǎn)與其最鄰近點(diǎn)之間的歐式距離，為距離的平均值， 是理論最鄰近距離，n 是茶類經(jīng)營點(diǎn)數(shù)量，A 是所在區(qū)域的面積，D 為點(diǎn)密度。

1.2.1.2 茶類經(jīng)營點(diǎn)分布均衡性分析

首位度、赫芬達(dá)爾系數(shù)和地理集中指數(shù)這3 個(gè)指標(biāo)可以反映全國茶類經(jīng)營點(diǎn)空間分布的集中程度和均衡性。首位度的公式為：

P1、P2分別代表指數(shù)最大省份和第二大省份的茶類經(jīng)營點(diǎn)數(shù)量。P 值越大，茶類經(jīng)營點(diǎn)數(shù)量就越集中，不同省份之間的研究對象就越不均衡。

赫芬達(dá)爾系數(shù)的公式為：

赫芬達(dá)爾系數(shù)H 可以反映研究對象在區(qū)域上的集聚程度，P1表示某個(gè)省份的茶類經(jīng)營點(diǎn)數(shù)量占總量的比值。赫芬達(dá)爾系數(shù)在1/n～1 之間變動(dòng)。如果茶類經(jīng)營點(diǎn)數(shù)量集中于某個(gè)省份，則指數(shù)值等于1；如果所有省份的數(shù)量均勻分布，指數(shù)值等于1/n。

地理集中指數(shù)的公式為：

其中，G 為全國茶類經(jīng)營點(diǎn)的地理集中指數(shù)，Xi為第i 個(gè)省份茶類經(jīng)營點(diǎn)的數(shù)量，T 為全國茶類經(jīng)營點(diǎn)總數(shù)，n 為全國省份總數(shù)。G 取值范圍為0 ～ 100，G 值越大，表明茶類經(jīng)營點(diǎn)在各地分布越集中；反之，越小則越分散。假設(shè)G0為茶類經(jīng)營點(diǎn)平均分布在各省份情況下的地理集中指數(shù)，若G ＞ G0，則表示茶類經(jīng)營點(diǎn)分布較為集中。

1.2.1.3 茶類經(jīng)營點(diǎn)分布密度分析

利用Arc GIS 10.2 軟件對全國茶類經(jīng)營點(diǎn)進(jìn)行核密度分析，直觀反映離散對象在連續(xù)區(qū)域的空間分布狀況。核密度分析使用的公式為：

其中，f(x)是第x 個(gè)茶類經(jīng)營點(diǎn)的核密度，h 為寬帶(h ＞0)，(x-Xi)表示估值點(diǎn)x 到事件Xi處的距離，n 為所有茶類經(jīng)營點(diǎn)數(shù)量，K 為核函數(shù)。f(x)越大表示該區(qū)域內(nèi)茶類經(jīng)營點(diǎn)分布越密集。

1.2.2 影響因素分析

1.2.2.1 研究方法——地理探測器

王勁峰等學(xué)者開發(fā)研究的地理探測器空間統(tǒng)計(jì)學(xué)方法，可以很好地反映同一區(qū)域的相似性以及不同區(qū)域間的差異性，進(jìn)而得出研究對象的空間異質(zhì)性，主要用于探測空間分異性以及地理要素對研究對象空間分布的影響因素[19]。

（1）因子探測：探測影響因子X 對茶銷售市場Y 的解釋力，其表達(dá)式為：

（2）交互作用探測：交互探測可以識(shí)別茶銷售市場不同的影響因子兩兩共同作用的結(jié)果，是否會(huì)增加或減弱對茶銷售市場的解釋力，或?qū)Σ桎N售市場的影響是否相互獨(dú)立。主要會(huì)產(chǎn)生五種結(jié)果類型（表1）。

1.2.2.2 影響因素與代表指標(biāo)的選擇

茶類經(jīng)營點(diǎn)的空間分布特征是由多種因素影響形成的，通過茶產(chǎn)業(yè)的自身特點(diǎn)以及對其他零售業(yè)空間分布影響因素的借鑒，選取以下影響因素作為研究內(nèi)容（表2）。

氣候：氣候主要對茶種植有直接的影響，同樣對人們飲茶習(xí)慣的養(yǎng)成也具有一定的作用。因此從中國氣象數(shù)據(jù)網(wǎng)中獲取各地年平均氣溫值作為氣候代表指標(biāo)。

茶供給：茶銷售作為茶產(chǎn)業(yè)鏈的重要節(jié)點(diǎn)，茶種植與茶生產(chǎn)能力作為產(chǎn)業(yè)鏈前端必然會(huì)影響茶銷售區(qū)域范圍與銷售能力，分別選取各地的茶播種面積和茶產(chǎn)量作為代表指標(biāo)。

表1 交互作用的五種類型Table 1 Five result types of factor interaction

表2 茶類經(jīng)營點(diǎn)空間分布影響因素Table 2 Influencing factors of spatial distribution of tea business sites

人口：人口作為零售市場的主要影響因素在指標(biāo)體系中必不可少，選取各地區(qū)總?cè)丝谧鳛楣潭ň幼⌒匀丝诘拇碇笜?biāo)，旅客周轉(zhuǎn)量作為流動(dòng)性人口的代表指標(biāo)。

經(jīng)濟(jì)水平：一個(gè)地區(qū)的經(jīng)濟(jì)發(fā)展水平最常見的表征指標(biāo)是地區(qū)GDP，考慮到茶葉消費(fèi)屬于非生活必需品，在禮品交往、休閑社交等方面作用較大，如果一個(gè)地區(qū)國內(nèi)貿(mào)易比較活躍，相對能夠帶動(dòng)茶消費(fèi)水平，因此把社會(huì)消費(fèi)品零售總額作為代表指標(biāo)進(jìn)行分析。

居民生活水平：國內(nèi)茶消費(fèi)市場的主體是居民，居民購買能力的大小會(huì)受到居民生活水平的限制，用居民人均可支配收入、居民最終消費(fèi)支出、城鄉(xiāng)居民養(yǎng)老保險(xiǎn)基本參保人數(shù)三個(gè)指標(biāo)分別代表居民的收入水平、消費(fèi)水平和生活保障水平。

2 結(jié)果與分析

2.1 茶銷售市場空間分布特征

2.1.1 茶類經(jīng)營點(diǎn)分布類型

運(yùn)用Arc GIS 10.2 軟件中的平均最鄰近分析方法測算得出：茶銷售市場空間分布的最鄰近比率R = 0.829854 ＜ 1，標(biāo)準(zhǔn)差Z = -5.921981，概率P = 0，通過了1%顯著性水平檢驗(yàn)，茶銷售市場各茶類經(jīng)營點(diǎn)平均觀測距離為90.44 km，預(yù)期平均距離為108.99 km。整體呈集聚分布，但集聚性稍弱。

2.1.2 茶類經(jīng)營點(diǎn)分布均衡性

四川省茶類經(jīng)營點(diǎn)數(shù)量最大P1= 150 40，廣東省茶類經(jīng)營點(diǎn)數(shù)量第二P2= 8580，通過式（2）計(jì)算得出首位度P = 1.7529 ＜ 2，表明茶類經(jīng)營點(diǎn)數(shù)量最大的省份和數(shù)量第二的省份不存在很大差距，結(jié)合圖1 可以看出茶類經(jīng)營點(diǎn)并未聚集在某一省份，在全國各省份都有分布，差距雖有但跨度不大，集中程度適當(dāng)。

根據(jù)式（3），全國茶類經(jīng)營點(diǎn)的赫芬達(dá)爾系數(shù)H 變動(dòng)范圍處于0.0323 ～ 1 之間，計(jì)算得出H = 0.2920，表明各省份茶類經(jīng)營點(diǎn)在全國分布有一定的集中性，但集聚程度不是很強(qiáng)。

從地理集中指數(shù)來看，全國茶類經(jīng)營點(diǎn)總數(shù)量T = 91112，全國省份n = 31，通過式（4）計(jì)算得出G = 25.5562288。假設(shè)全國91112 個(gè)茶類經(jīng)營點(diǎn)平均分布在31 個(gè)省份，即每個(gè)省份平均茶類經(jīng)營點(diǎn)數(shù)量為91112/31 個(gè)，通過公式得出G0= 17.96055302。可發(fā)現(xiàn)G ＞G0，茶類經(jīng)營點(diǎn)在全國分布相對集中。

綜合3 個(gè)指數(shù)值來看，全國各省份間茶類經(jīng)營點(diǎn)呈現(xiàn)一定程度的集中分布，但程度不突出，結(jié)構(gòu)較為正常。

2.1.3 茶類經(jīng)營點(diǎn)分布密度

通過Arc GIS 10.2 軟件對全國各茶類經(jīng)營點(diǎn)分布核密度實(shí)現(xiàn)了可視化表達(dá)，繪制得出核密度分布圖，直觀展現(xiàn)茶類經(jīng)營點(diǎn)散布狀況。按照核密度值將茶類經(jīng)營點(diǎn)分布劃分為五個(gè)等級(jí)層次。如圖可以發(fā)現(xiàn)大部分省份大多是以省會(huì)為主形成高密度核心并向外擴(kuò)散，全國范圍內(nèi)最高密度核心地區(qū)包括以廣東廣州、江蘇南京和河南鄭州為核心的周邊部分城市；次一級(jí)密度核心地區(qū)如上海市、四川成都、重慶、湖南長沙、湖北武漢、遼寧沈陽等及其周邊區(qū)域，主要涵蓋了3 個(gè)高密度核心的周邊地區(qū)以及部分省會(huì)城市；另外，我國中、東部其他大部分地區(qū)以及西部部分地區(qū)如廣西、云南、貴陽、四川、陜西、甘肅、寧夏、新疆烏魯木齊周邊等形成再次一級(jí)密度核心；西部大部分地區(qū)、中部內(nèi)蒙古和東部黑龍江等大多未形成相對規(guī)模密集區(qū)。

圖3 茶類經(jīng)營點(diǎn)核密度分布圖Figure 3 Core density distribution map of tea business sites

以胡煥庸線為界，茶類經(jīng)營點(diǎn)在我國版圖上“東密西疏”的特征非常明顯，東南半壁茶類經(jīng)營點(diǎn)分布核密度值較高，反之西北半壁的分布較為稀疏，說明大部分茶類經(jīng)營點(diǎn)分布聚集狀況與我國整體經(jīng)濟(jì)發(fā)展水平與人口分布狀況密切相關(guān)，與茶種植和生產(chǎn)的空間布局并未完全重合。

2.2 茶銷售市場影響因素

2.2.1 茶銷售市場的單因子探測

利用地理探測器中的因子探測器，對10 個(gè)代表指標(biāo)探測對茶銷售市場空間分布的影響。結(jié)果（表3）顯示：9 個(gè)變量基本上通過5%的顯著性水平檢驗(yàn)，上述影響因素對茶銷售市場空間分布具有較為顯著的影響。

表3 茶銷售市場影響因子探測結(jié)果Table 3 Detection results of influencing factors of tea sales market

由表3 結(jié)果可知，固定人口中的指標(biāo)因子總?cè)丝冢╔4）和國內(nèi)貿(mào)易水平中的指標(biāo)因子社會(huì)消費(fèi)品零售總額（X7）的q 值分別為0.6081 和0.5295，對我國茶消費(fèi)市場的空間分布解釋力最強(qiáng)，為主導(dǎo)因子。經(jīng)濟(jì)水平中的指標(biāo)因子地區(qū)生產(chǎn)總值（X6）、居民生活保障中的指標(biāo)因子城鄉(xiāng)居民基本參保人數(shù)（X10）、居民消費(fèi)水平中的居民最終消費(fèi)支出（X9）、茶生產(chǎn)規(guī)模中的茶產(chǎn)量（X3）、茶種植規(guī)模中的茶種植面積（X2）以及流動(dòng)人口中的旅客周轉(zhuǎn)量（X5）也對我國茶消費(fèi)市場有較大影響，地理探測器的q 值均大于35%，對應(yīng)的q 值分別為0.4767、0.4206、0.3990、0.3727、0.3684 和0.3532。此外，氣候中的指標(biāo)因子平均氣溫（X1）通過了5%的p 值檢驗(yàn)，q 值為0.1549，居民收入水平中的居民人均可支配收入（X8）沒有通過p 值檢驗(yàn)，而且驅(qū)動(dòng)力相對較小，q 值為0.1055，表明平均氣溫和居民人均可支配收入對我國茶銷售市場的作用偏弱。

2.2.2 茶銷售市場影響因素的交互探測

由表4 可見，我國茶銷售市場的10 個(gè)探測因子在交互后得到的55 個(gè)雙因子探測q 值，交互作用有兩種結(jié)果：雙因子增強(qiáng)和非線性增強(qiáng)，不存在獨(dú)立作用或作用減弱的因子。即任何兩個(gè)因子的交互作用對茶銷售市場的影響都高于單因子作用。

與單因子作用相比，任意雙因子交互探測q 值對茶銷售市場的解釋力增強(qiáng)，表明各個(gè)驅(qū)動(dòng)因子之間并非相互獨(dú)立，而是具有協(xié)同增強(qiáng)的作用。如氣候中的指標(biāo)因子平均氣溫（X1）和居民收入水平中的居民人均可支配收入（X8）在所有驅(qū)動(dòng)因子中的單個(gè)影響力是最弱的，但二者交互后產(chǎn)生的影響力不僅加強(qiáng)，還大于二者的q值之和，產(chǎn)生了明顯的非線性增強(qiáng)的效果。此外，總?cè)丝冢╔4）、社會(huì)消費(fèi)品零售總額（X7）、地區(qū)生產(chǎn)總值（X6）是對茶銷售市場的解釋作用最強(qiáng)的單因子，兩兩交互后結(jié)果為q（X4∩X7）= 0.9244，q （X4∩X6） = 0.6762，q （X7∩X6） = 0.5531， 三組交互因子的結(jié)果中X4∩X7的q 值最大。進(jìn)一步分析表4 的交互結(jié)果可知，對我國茶銷售市場解釋作用最強(qiáng)的兩個(gè)交互因子為X3∩X4，對應(yīng)的q 值為0.9423，大于X4∩X7交互得出的q 值，可見對我國茶銷售市場驅(qū)動(dòng)力最強(qiáng)的單因子再與其他影響因子交互之后解釋力不一定最強(qiáng)，而部分影響因子如平均氣溫（X1）在單獨(dú)發(fā)揮作用時(shí)產(chǎn)生的驅(qū)動(dòng)力弱一些，但是同其他影響因子通過動(dòng)態(tài)交互得出的效果卻非常突出，因此不能輕易忽視影響因素之間的關(guān)系。我國茶銷售市場的空間分布是多種因素非線性耦合的結(jié)果，需要用綜合系統(tǒng)的眼光來協(xié)調(diào)看待所有影響因子對我國茶銷售市場空間分布的作用。

表4 茶銷售市場空間分異的交互探測結(jié)果Table 4 Interactive detection results of spatial differentiation of tea sales market

3 結(jié)論

通過以高德地圖POI 數(shù)據(jù)獲取的除港澳臺(tái)外全國31 個(gè)省份346 個(gè)城市的茶類經(jīng)營點(diǎn)信息作為研究數(shù)據(jù)，在不考慮茶類經(jīng)營點(diǎn)的銷售能力和經(jīng)營水平的前提下，以空間結(jié)構(gòu)特征的視角進(jìn)行銷售市場的布局分析，并探析了影響因素的作用機(jī)制，得出以下結(jié)論：

我國茶銷售市場在各省份城市間的分布呈現(xiàn)一定程度的聚集特征，但程度不突出，結(jié)構(gòu)較為均衡。通過核密度分析，茶類經(jīng)營點(diǎn)在我國版圖上“東密西疏”的特征明顯，廣州、南京、鄭州、上海、成都、重慶、長沙、武漢、沈陽等城市及周邊是核密度較集中的區(qū)域。

通過地理探測器對茶銷售市場影響因素的空間異質(zhì)性探測，得出影響因素作用從大到小的排序依次是：總?cè)丝凇⑸鐣?huì)消費(fèi)品零售總額、地區(qū)生產(chǎn)總值、城鄉(xiāng)居民基本參保人數(shù)、居民最終消費(fèi)支出、茶葉產(chǎn)量、茶種植面積、旅客周轉(zhuǎn)量、平均氣溫、居民人均可支配收入。人口基數(shù)和經(jīng)濟(jì)發(fā)達(dá)程度是茶銷售市場空間格局的主導(dǎo)因素，居民生活的保障與消費(fèi)水平的高度是重要影響因素，茶產(chǎn)業(yè)的種植與生產(chǎn)空間雖然也會(huì)產(chǎn)生影響，但是作用程度較弱，氣候與居民收入的作用力最低，以上影響因素在兩兩交互后作用能力全部增強(qiáng)，構(gòu)成一個(gè)復(fù)雜化非線性的茶銷售市場影響系統(tǒng)。