劉桂東，劉紅玲，劉杰,3，趙軍產(chǎn)

(1.武漢紡織大學(xué) 數(shù)學(xué)與計(jì)算機(jī)學(xué)院應(yīng)用數(shù)學(xué)系，湖北 武漢 430073；2.湖南工商大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院應(yīng)用數(shù)學(xué)系，湖南 長沙 410205；3.武漢紡織大學(xué) 非線性科學(xué)研究中心，湖北 武漢 430073)

一、引言

紡織紗線是用各種紡織纖維加工而成的具有一定細(xì)度的紡織產(chǎn)品，可用于織布、制繩、制線、針織和刺繡等。紡織紗線可以分為短纖維紗，連續(xù)長絲等。紡織紗線是織布企業(yè)的主要原料，也是重要的民生物資。紗線貿(mào)易歷來就是古代絲綢之路和海上絲綢之路上流通的重要商品。2013年9月和10月，中國國家主席習(xí)近平在出訪哈薩克斯坦和印度尼西亞時(shí)，先后提出共建“絲綢之路經(jīng)濟(jì)帶”和“21世紀(jì)海上絲綢之路”的重大倡議(“一帶一路倡議”英文為Belt and Road Initiative，縮寫為BRI)[1]。2015年3月，中國發(fā)布《推動(dòng)共建絲綢之路經(jīng)濟(jì)帶和21世紀(jì)海上絲綢之路的愿景與行動(dòng)》。共建“一帶一路”倡議以“共商共建共享”為原則，以“和平合作、開放包容、互學(xué)互鑒、互利共贏”的絲綢之路精神為指引，以“政策溝通、設(shè)施聯(lián)通、貿(mào)易暢通、資金融通、民心相通”為重點(diǎn)，已經(jīng)從理念轉(zhuǎn)化為行動(dòng)，從愿景轉(zhuǎn)化為現(xiàn)實(shí)，從倡議轉(zhuǎn)化為全球廣受歡迎的公共產(chǎn)品[2]。2017年5月，首屆“一帶一路”國際合作高峰論壇在北京成功召開。2019年4月在北京又成功召開了第二屆“一帶一路”國際合作高峰論壇。截至2020年1月底，中國已經(jīng)同138個(gè)國家和30個(gè)國際組織簽署200份共建“一帶一路”合作文件。共建“一帶一路”國家已由亞歐延伸至非洲、拉美、南太等區(qū)域。紡織紗線作為最重要的織布原料和重要民生物資，在世界各國都有大量需求，由于地理位置和資源稟賦的差異，有些國家紡織紗線的產(chǎn)量較高，出口需求較大。由于技術(shù)和工業(yè)發(fā)展水平的差異，各國對(duì)不同種類和品質(zhì)的紗線，互相之間都有需求，這就為紡織紗線貿(mào)易提供了天然的推動(dòng)力。共建“一帶一路”倡議的提出，以及參與共建各國的共同努力應(yīng)該能夠進(jìn)一步加強(qiáng)各國在紡織紗線貿(mào)易方面的合作，但在發(fā)展的過程中也會(huì)碰到一些問題，需要認(rèn)真研究，更好地促進(jìn)各國紡織紗線貿(mào)易的發(fā)展。

Watts和Strogatz(1998)提出小世界網(wǎng)絡(luò)模型[3]，Barabasi和Albert(1999)提出無標(biāo)度網(wǎng)絡(luò)模型[4]，這兩項(xiàng)開創(chuàng)性研究為各領(lǐng)域研究者提供了一個(gè)新的分析復(fù)雜實(shí)測數(shù)據(jù)的視角：網(wǎng)絡(luò)科學(xué)。自此以來，網(wǎng)絡(luò)科學(xué)成為定量分析復(fù)雜經(jīng)濟(jì)系統(tǒng)的新工具。越來越多的經(jīng)濟(jì)學(xué)者將復(fù)雜網(wǎng)絡(luò)分析法用于商品貿(mào)易活動(dòng)的建模與分析。Serrano(2003)[5]等分析了世界貿(mào)易網(wǎng)絡(luò)的非隨機(jī)、小世界、無標(biāo)度、高集聚特征。Garlaschelli等(2005)[6]研究了世界貿(mào)易網(wǎng)絡(luò)的結(jié)構(gòu)及演化。Fagiolo等(2007)[7]研究了加權(quán)世界貿(mào)易網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)特征。Fagiolo等(2009)[8]則對(duì)世界貿(mào)易網(wǎng)絡(luò)的動(dòng)力學(xué)和演化進(jìn)行了進(jìn)一步研究。Benedictis等(2011)[9]研究發(fā)現(xiàn)世界貿(mào)易網(wǎng)絡(luò)中各成員國在貿(mào)易伙伴選擇中異質(zhì)性較強(qiáng)。

自BRI提出以來，對(duì)“一帶一路”沿線國家及共建參與國家之間商品貿(mào)易網(wǎng)絡(luò)關(guān)系特征演化的研究引起了我國學(xué)者們的廣泛關(guān)注。宋周鶯等(2017)[10]分析了“一帶一路”貿(mào)易網(wǎng)絡(luò)與全球貿(mào)易網(wǎng)絡(luò)的拓?fù)潢P(guān)系，從總體商品貿(mào)易層面對(duì)“一帶一路”貿(mào)易在全球的地位進(jìn)行了分析。許和連等(2015)[11]，運(yùn)用網(wǎng)絡(luò)中心性和模體分析法，考察“一帶一路”上65個(gè)國家的高端制造業(yè)貿(mào)易網(wǎng)絡(luò)中各國地位和貿(mào)易模式。姚星等(2018)[12]對(duì)沿線國家服務(wù)中間投入的網(wǎng)絡(luò)結(jié)構(gòu)及其影響因素運(yùn)用社會(huì)網(wǎng)絡(luò)分析法進(jìn)行探討。馬遠(yuǎn)和徐俐俐(2017)[13]則對(duì)“一帶一路”沿線國家天然氣貿(mào)易網(wǎng)絡(luò)結(jié)構(gòu)特征、微觀模式及其影響因素進(jìn)行了分析。許和連等(2015)[11]和馬遠(yuǎn)等(2017)[13]將復(fù)雜網(wǎng)絡(luò)分析的方法用到了“一帶一路”沿線國家某一大類或某一具體商品的貿(mào)易，研究了相應(yīng)網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)特征、重要國家在網(wǎng)絡(luò)中的地位及其演化規(guī)律等，對(duì)未來加強(qiáng)相關(guān)領(lǐng)域的貿(mào)易合作具有重要的參考價(jià)值。自BRI提出以來，也有越來越多的研究者對(duì)我國紡織服裝行業(yè)積極參與“一帶一路”戰(zhàn)略合作等問題進(jìn)行了理論與實(shí)證性研究。中國紡織工業(yè)聯(lián)合會(huì)會(huì)長王天凱認(rèn)為“一帶一路”戰(zhàn)略為紡織業(yè)國際布局創(chuàng)造了有利環(huán)境。國家明確了鼓勵(lì)紡織等優(yōu)勢產(chǎn)業(yè)“走出去”思路，并對(duì)管理體制進(jìn)行改革，顯著提高了對(duì)外投資的便利性。眾多的政治、外交和經(jīng)濟(jì)資源投入，將為包括紡織業(yè)在內(nèi)的中國企業(yè)在“一帶一路”沿線國家投資創(chuàng)造盡可能安全、穩(wěn)定的中長期經(jīng)濟(jì)地理環(huán)境。中國紡織行業(yè)需抓住歷史機(jī)遇，通過全球產(chǎn)能合作與優(yōu)質(zhì)資源配置，更好地實(shí)現(xiàn)建設(shè)紡織科技強(qiáng)國、品牌強(qiáng)國、可持續(xù)發(fā)展強(qiáng)國和人才強(qiáng)國的目標(biāo)[14]。中國紡織工業(yè)聯(lián)合會(huì)副會(huì)長楊紀(jì)朝(2017)[15]認(rèn)為我國紡織服裝業(yè)應(yīng)促進(jìn)與沿線國家的貿(mào)易伙伴關(guān)系逐步升級(jí)到全面伙伴關(guān)系；我國紡織服裝企業(yè)應(yīng)借助“一帶一路”發(fā)展戰(zhàn)略提高我國紡織工業(yè)國際化水平；我國紡織服裝業(yè)加入“一帶一路”對(duì)持續(xù)性體現(xiàn)雙邊互利、共享發(fā)展成果、樹立國家正面形象具有重大意義；我國紡織服裝業(yè)參與“一帶一路”為調(diào)整產(chǎn)業(yè)自身結(jié)構(gòu)、轉(zhuǎn)型升級(jí)爭取了時(shí)間和空間；“一帶一路”為地方紡織服裝產(chǎn)業(yè)集群發(fā)展帶來新機(jī)遇。丁鋒、姚新超(2018)[16]對(duì)中國與“一帶一路”沿線主要國家紡織品服裝出口現(xiàn)狀及各國競爭力進(jìn)行了分析，結(jié)果表明：中國與“一帶一路”沿線主要紡織品服裝出口國在國際市場上均具有很強(qiáng)的競爭力，中國除了紡織纖維及其廢料競爭力較弱以外，在紡紗及有關(guān)產(chǎn)品和服裝及衣著附件兩類產(chǎn)品中都具有較強(qiáng)的競爭力。因此，中國與“一帶一路”沿線紡織品服裝國家可以深化產(chǎn)業(yè)合作，不斷提升紡織品服裝出口的國際競爭力。

作為最重要的紡織原料，紡織紗線一直是“一帶一路”沿線國家間歷史悠久的傳統(tǒng)貿(mào)易產(chǎn)品，但目前對(duì)“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的結(jié)構(gòu)特征及演化趨勢還缺乏比較系統(tǒng)的研究文獻(xiàn)。探討“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的宏觀和微觀拓?fù)浣Y(jié)構(gòu)特征、網(wǎng)絡(luò)中重要國家節(jié)點(diǎn)的特征及其演化，有助于深入了解“一帶一路”沿線國家紡織紗線貿(mào)易合作的深層次特點(diǎn)，總結(jié)過去貿(mào)易合作中的經(jīng)驗(yàn)教訓(xùn)、找到現(xiàn)有貿(mào)易合作的短板和加強(qiáng)未來貿(mào)易合作的著力點(diǎn)?！耙粠б宦贰背h從提出到現(xiàn)在已有六年多的時(shí)間，這一旨在加強(qiáng)參與各國互聯(lián)互通、共同發(fā)展的倡議對(duì)“一帶一路”沿線國家紡織紗線貿(mào)易合作產(chǎn)生了何種影響？在哪些方面取得了進(jìn)步？未來應(yīng)如何進(jìn)一步加強(qiáng)合作？都是值得深入研究的問題。

基于2006—2018年間較完備的聯(lián)合國商品貿(mào)易數(shù)據(jù)，本文運(yùn)用復(fù)雜網(wǎng)絡(luò)建模分析和高階結(jié)關(guān)系分析法，擬對(duì)“一帶一路”沿線國家紡織紗線進(jìn)出口貿(mào)易網(wǎng)絡(luò)開展建模分析，研究該貿(mào)易網(wǎng)絡(luò)結(jié)構(gòu)特征、內(nèi)部高階結(jié)關(guān)系等演化發(fā)展趨勢。在使用聯(lián)合國商品貿(mào)易統(tǒng)計(jì)數(shù)據(jù)庫時(shí)，鑒于每個(gè)年度都存在部分貿(mào)易活動(dòng)參與國并未及時(shí)、主動(dòng)報(bào)告本國相關(guān)產(chǎn)品進(jìn)出口數(shù)據(jù)，本文研究過程中，對(duì)收集的原始數(shù)據(jù)均采取了雙向“取大原則”進(jìn)行修正，即將一個(gè)國家自己主動(dòng)報(bào)告的向他國進(jìn)口產(chǎn)品的數(shù)據(jù)與他國報(bào)告的向該國出口相同產(chǎn)品的數(shù)據(jù)之間的最大值作為該國向他國進(jìn)口這類產(chǎn)品數(shù)據(jù)的修正值。這一處理方法具有很好的修正效用，相關(guān)研究結(jié)果能夠比較接近真實(shí)地反映國家間實(shí)際發(fā)生的貿(mào)易行為[17，24]。

二、“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的構(gòu)建

本文所用到的紡織紗線貿(mào)易原始數(shù)據(jù)均來自聯(lián)合國商品貿(mào)易統(tǒng)計(jì)數(shù)據(jù)庫①。首先，收集了2006—2018年間聯(lián)合國商品貿(mào)易統(tǒng)計(jì)數(shù)據(jù)庫中“一帶一路”沿線65個(gè)國家報(bào)告的紡織紗線進(jìn)出口數(shù)據(jù)。如前所述，數(shù)據(jù)清洗階段發(fā)現(xiàn)，每年都有貿(mào)易活動(dòng)參與國未報(bào)告本國紡織紗線進(jìn)出口數(shù)據(jù)，但其他國家報(bào)告的數(shù)據(jù)中卻有向該國進(jìn)出口紡織紗線的逆向數(shù)據(jù)存在。為確保分析能接近真實(shí)地反映“一帶一路”沿線國家間紡織紗線貿(mào)易的真實(shí)情況，本文對(duì)原始的紡織紗線進(jìn)出口數(shù)據(jù)進(jìn)行了修正處理，該方法的可靠性分析可參見文獻(xiàn)[17，24]。

將“一帶一路”沿線國家視為擬構(gòu)建“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的65個(gè)節(jié)點(diǎn)，記之為：

VBR={v1,v2,…,v65}

它表示“一帶一路”沿線65個(gè)國家構(gòu)成的集合。記權(quán)重矩陣為：

WBR=(wij)65×65

其中，wij表示國家vi從國家vj進(jìn)口紡織紗線的修正額度。再將貿(mào)易網(wǎng)絡(luò)鄰接矩陣記為：

ABR=(aij)65×65

其中：

由此得到：

GBR=(VBR,ABR,WBR)

GBR就是“一帶一路”沿線國家紡織纖維進(jìn)口貿(mào)易的加權(quán)、有向網(wǎng)絡(luò)。

三、“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)及其閾值子網(wǎng)演化分析

(一)“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)演化分析

按照前述方法，本文構(gòu)建了歷年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)。本文研究的相關(guān)領(lǐng)域的貿(mào)易數(shù)據(jù)基本涵蓋了全部最新官方數(shù)據(jù)(因2019年度部分國家的數(shù)據(jù)未完全入庫，為確保研究數(shù)據(jù)的可靠性，構(gòu)建網(wǎng)絡(luò)覆蓋年度區(qū)間為2006—2018)。數(shù)據(jù)選擇與時(shí)間范圍設(shè)定，基本確保了相關(guān)研究結(jié)果能反映出“一帶一路”倡議提出前后涉及的兩個(gè)階段紗線貿(mào)易關(guān)系的真實(shí)演化發(fā)展特征。

下面的圖1給出了2016年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)示意圖。表1給出了2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的整體網(wǎng)逐年特征量。

圖1 2016年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)示意圖

表1 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的整體特征

從表1可以看出：自2006年起“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)一直保持為65，即“一帶一路”沿線所有國家都參與了網(wǎng)內(nèi)的紡織紗線進(jìn)出口貿(mào)易；除2008年網(wǎng)絡(luò)直徑為2，其他年份網(wǎng)絡(luò)直徑都為3，網(wǎng)絡(luò)平均路徑長度在1.379到1.508之間小幅波動(dòng)，這說明“一帶一路”沿線國家紡織紗線貿(mào)易聯(lián)系總體較為緊密，網(wǎng)絡(luò)互通性較好②。

我們注意到，在2006—2018年間，“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度、加權(quán)平均度總體呈增長態(tài)勢，但有所波動(dòng)，平均路徑長度基本呈下降態(tài)勢，2016年起稍有回升，這四項(xiàng)特征的變化曲線如圖2所示。

圖2 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)特征演化圖

從網(wǎng)絡(luò)密度和節(jié)點(diǎn)平均度這兩項(xiàng)結(jié)構(gòu)指標(biāo)來看，從2006—2007年，這兩項(xiàng)指標(biāo)大幅提高，2007—2008年增長幅度仍然較大，受2008年全球金融危機(jī)的影響，“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的密度和節(jié)點(diǎn)平均度沒有保住上升的勢頭，但隨著金融危機(jī)影響的逐步消退，2009年起繼續(xù)回升；由于國際貿(mào)易保護(hù)主義抬頭等因素的影響，2010—2012年，這兩項(xiàng)指標(biāo)回升乏力，基本持平。2013年起，BRI的提出為提升“一帶一路”沿線國家紡織紗線貿(mào)易注入了強(qiáng)大的動(dòng)力，“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的密度、節(jié)點(diǎn)平均度這兩個(gè)指標(biāo)在2013年起明顯上升，2016年達(dá)到高點(diǎn)，反映出“一帶一路”沿線國家對(duì)加強(qiáng)各國之間紡織紗線貿(mào)易有著強(qiáng)大的需求，BRI的提出促進(jìn)了“一帶一路”沿線國家之間紡織紗線貿(mào)易合作的參與度。值得注意的是，2016年以后，網(wǎng)絡(luò)密度和節(jié)點(diǎn)平均度又有所下降，即網(wǎng)絡(luò)中有一部分建立起來的貿(mào)易關(guān)系又中斷了。網(wǎng)絡(luò)加權(quán)平均度的走勢與網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度不同的地方是其最高點(diǎn)出現(xiàn)在2014年，隨后緩慢回落，但保持在較高的水平，且2017年起又恢復(fù)上升，這說明參與網(wǎng)絡(luò)內(nèi)貿(mào)易各國的參與力度維持得很好。網(wǎng)絡(luò)平均路徑長度在2016年之前是逐年下降，說明網(wǎng)絡(luò)互通性越來越好，但自2017年起，網(wǎng)絡(luò)平均路徑長度又有所上升。

綜合以上關(guān)于2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)特征演化趨勢分析可知，BRI的提出受到了“一帶一路”沿線國家的積極響應(yīng)，早期各國參與共建“一帶一路”的熱情較高，但共建“一帶一路”，實(shí)現(xiàn)“五通”合作需要建設(shè)大量基礎(chǔ)設(shè)施，加強(qiáng)各國在經(jīng)濟(jì)、社會(huì)、文化等各個(gè)領(lǐng)域的全面合作，其建設(shè)效果難以在很短的時(shí)間內(nèi)立即顯現(xiàn)。在2014年以后，“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的密度和節(jié)點(diǎn)平均度有所回落，說明一些沿線國家理性回歸，但是隨著共建“一帶一路”的持續(xù)推進(jìn)，各國之間紡織紗線貿(mào)易進(jìn)一步加強(qiáng)是大勢所趨。另外，根據(jù)以往聯(lián)合國貿(mào)易數(shù)據(jù)庫數(shù)據(jù)更新的情況，2017年、2018年的數(shù)據(jù)很可能還會(huì)有更新，網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度、加權(quán)平均度都有可能上升，平均路徑長度都有可能進(jìn)一步縮短。

(二) “一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)及其閾值子網(wǎng)演化分析

以上構(gòu)建“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的連邊規(guī)則是，連一條從節(jié)點(diǎn)i指向節(jié)點(diǎn)j的有向邊，即只要i國報(bào)告的向j國的紡織紗線進(jìn)口額，或j國報(bào)告的向i國的紡織紗線出口額至少有一個(gè)大于0，則認(rèn)為存在i國對(duì)j國的紡織紗線進(jìn)口貿(mào)易合作關(guān)系。在數(shù)據(jù)清洗階段，發(fā)現(xiàn)大量的wij取值雖然大于0，但數(shù)值很小，有的甚至只有1美元，按前面的網(wǎng)絡(luò)構(gòu)建方法，也會(huì)認(rèn)定這兩個(gè)國家之間存在紡織紗線進(jìn)口貿(mào)易合作關(guān)系。這樣的認(rèn)定，從國際貿(mào)易的角度來講難免牽強(qiáng)，同時(shí)也會(huì)對(duì)網(wǎng)絡(luò)主要貿(mào)易合作關(guān)系的分析產(chǎn)生干擾。

研究過程中，欲克服這一不利影響的一個(gè)處理方法就是：設(shè)定一個(gè)閾值λ，從原網(wǎng)絡(luò)中過濾掉權(quán)小于λ的邊，只保留權(quán)大于或等于λ的邊，并且刪去孤立節(jié)點(diǎn)，得到一個(gè)新的網(wǎng)絡(luò)：

其中：

對(duì)于紡織紗線貿(mào)易網(wǎng)絡(luò)而言，選取一個(gè)數(shù)值λ作為閾值，就是給出一個(gè)認(rèn)定兩個(gè)國家之間存在實(shí)質(zhì)性紡織紗線貿(mào)易合作關(guān)系的貿(mào)易額度臨界值λ：如果wij≥λ，則認(rèn)為存在實(shí)質(zhì)性的i國向j國的紡織紗線進(jìn)口貿(mào)易合作關(guān)系；如果wij<λ，則認(rèn)為不存在實(shí)質(zhì)性的i國向j國的紡織紗線進(jìn)口貿(mào)易合作關(guān)系。目前，在國際貿(mào)易領(lǐng)域，沒有一個(gè)公認(rèn)的認(rèn)定產(chǎn)品國際貿(mào)易實(shí)質(zhì)性合作關(guān)系存在的貿(mào)易額臨界值標(biāo)準(zhǔn)。在各種復(fù)雜網(wǎng)絡(luò)的過濾縮網(wǎng)分析中，研究者們嘗試了各種不同的規(guī)則來確定閾值λ的選取，如孫延風(fēng)，王朝勇在閾值網(wǎng)絡(luò)中提出了使用分位數(shù)來確定閾值的方法,將相關(guān)系數(shù)6等分，取第4部分的中點(diǎn)作為閾值(2018)。王芮(2018)[19]通過確定有效閾值區(qū)間和觀察最大連通子圖節(jié)點(diǎn)數(shù)的變化規(guī)律，找到一個(gè)最佳閾值。

不管采用何種規(guī)則，我們認(rèn)為：在認(rèn)定國際貿(mào)易實(shí)質(zhì)性合作關(guān)系存在的貿(mào)易額臨界值λ的選取中，都難免陷入“禿頭悖論”的境地。對(duì)于本文研究的紡織紗線國際貿(mào)易網(wǎng)絡(luò)而言，確定閾值λ，過濾縮減網(wǎng)絡(luò)的目的是為了更好地分析“一帶一路”沿線國家紡織紗線國際貿(mào)易網(wǎng)絡(luò)中主要的貿(mào)易合作關(guān)系，以及由這些主要關(guān)系交織而成的貿(mào)易合作網(wǎng)絡(luò)的拓?fù)涮卣骷捌溲莼?，進(jìn)一步可以以此為基礎(chǔ)對(duì)網(wǎng)絡(luò)進(jìn)行優(yōu)化、控制等。閾值λ的選擇應(yīng)該與這一目標(biāo)相符。

因此，針對(duì)“一帶一路”沿線國家紡織紗線國際貿(mào)易網(wǎng)絡(luò)過濾閾值的選取，本文提出如下可操作準(zhǔn)則：準(zhǔn)則(1)，閾值的選取不能導(dǎo)致網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)量的大量減少，使得網(wǎng)絡(luò)縮減為少數(shù)幾個(gè)大國的貿(mào)易網(wǎng)絡(luò)；準(zhǔn)則(2)，能有效剔除大部分符合常識(shí)的不合理連邊關(guān)系；準(zhǔn)則(3)，縮減后的子網(wǎng)絡(luò)的貿(mào)易總額相對(duì)原貿(mào)易網(wǎng)絡(luò)貿(mào)易總額的比例下降不應(yīng)太大；準(zhǔn)則(4)，閾值的選取符合國際貿(mào)易的實(shí)際操作習(xí)慣。

基于準(zhǔn)則4，本文首先選擇1萬、10萬、100萬、1000萬作為閾值，對(duì)原來的2016年度紡織紗線國際貿(mào)易網(wǎng)絡(luò)進(jìn)行過濾縮減，并將相應(yīng)的骨干網(wǎng)(子網(wǎng)絡(luò))稱為1萬俱樂部、10萬俱樂部、100萬俱樂部、1000萬俱樂部。這四個(gè)骨干網(wǎng)(子網(wǎng)絡(luò))的節(jié)點(diǎn)數(shù)、邊數(shù)、網(wǎng)絡(luò)貿(mào)易額三個(gè)特征及其與原網(wǎng)絡(luò)相應(yīng)特征的對(duì)比如表2所示：

表2 2016年“一帶一路”紡織紗線貿(mào)易閾值子網(wǎng)特征對(duì)比

從上表可知，當(dāng)閾值取為1萬、10萬、100萬、1000萬時(shí)，子網(wǎng)的節(jié)點(diǎn)數(shù)都為65，與原網(wǎng)絡(luò)完全一致；當(dāng)閾值依次取為1萬、10萬、100萬、1000萬時(shí)，子網(wǎng)的邊數(shù)依次明顯下降，子網(wǎng)邊數(shù)占原網(wǎng)絡(luò)邊數(shù)的比例依次為78%、60%、38%、17%；當(dāng)閾值依次取為1萬、10萬、100萬時(shí)，子網(wǎng)的加權(quán)平均度的下降幅度都很小，子網(wǎng)加權(quán)平均度占原網(wǎng)絡(luò)加權(quán)平均度的比例依次為99.999%、99.976%、99.732%，當(dāng)閾值變?yōu)?000萬時(shí)，子網(wǎng)加權(quán)平均度占原網(wǎng)絡(luò)加權(quán)平均度的比例為97.496%，比前面的下降幅度明顯加大。特別是當(dāng)閾值變?yōu)?000萬時(shí)，子網(wǎng)的連通分支數(shù)從此前的1變?yōu)?，網(wǎng)絡(luò)連通性發(fā)生了本質(zhì)的改變。綜合以上特點(diǎn)，本文將選取100萬為閾值，對(duì)2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)進(jìn)行過濾縮減，分析子網(wǎng)絡(luò)的主要拓?fù)涮卣骷捌溲莼?guī)律。

表3給出了2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易閾值100萬骨干網(wǎng)(子網(wǎng)絡(luò))的整體網(wǎng)逐年特征量，其中包括了網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度、加權(quán)平均度和模塊度的演化趨勢。

圖3 2016年“一帶一路”沿線國家紡織紗線貿(mào)易閾值100萬骨干網(wǎng)(子網(wǎng)絡(luò))拓?fù)涫疽鈭D

表3 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易閾值100萬骨干網(wǎng)(子網(wǎng)絡(luò))特征

圖4反映的是2006—2018年間“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)100萬閾值子網(wǎng)的網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度、節(jié)點(diǎn)加權(quán)平均度和網(wǎng)絡(luò)平均路徑長度的變化趨勢。

從表3和圖4可以看出，2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易100萬閾值子網(wǎng)在網(wǎng)絡(luò)密度、節(jié)點(diǎn)平均度、加權(quán)平均度三項(xiàng)網(wǎng)絡(luò)特征的總體演化趨勢與原網(wǎng)絡(luò)基本一致，三項(xiàng)指標(biāo)總體上是上升的，但其中的轉(zhuǎn)折點(diǎn)更明顯。受到2008年全球金融危機(jī)的影響，這三項(xiàng)指標(biāo)度出現(xiàn)了大幅下降，2009年起，隨著金融危機(jī)的消退，止跌快速回升，隨后緩慢上升，甚至小幅回落；2013年BRI提出后，三項(xiàng)指標(biāo)都又大幅提升，持續(xù)一兩年后稍有回調(diào)，然后再回升。這說明金融危機(jī)對(duì)“一帶一路”沿線國家紡織紗線實(shí)質(zhì)性的貿(mào)易產(chǎn)生了巨大的不利影響，而BRI的提出則對(duì)“一帶一路”沿線國家紡織紗線實(shí)質(zhì)性的貿(mào)易合作的廣度和深度都產(chǎn)生了積極的推動(dòng)作用，并且這種推動(dòng)作用還在持續(xù)擴(kuò)大。在網(wǎng)絡(luò)平均路徑長度這一指標(biāo)上，100萬閾值子網(wǎng)的演化趨勢與原網(wǎng)絡(luò)有比較明顯的差異，原網(wǎng)絡(luò)中，在2016年之前網(wǎng)絡(luò)平均路徑一直在逐年下降，2017年、2018年有所回升，很可能是部分?jǐn)?shù)據(jù)報(bào)告有待更新的結(jié)果；但在100萬閾值子網(wǎng)中，自2014年起，網(wǎng)絡(luò)平均路徑就呈上升態(tài)勢，“一帶一路”沿線國家之間實(shí)質(zhì)性的紡織紗線貿(mào)易鏈條拉長了。

圖4 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易100萬閾值子網(wǎng)特征演化圖

四、“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)高階結(jié)演化分析

世界貿(mào)易網(wǎng)絡(luò)(World Trade Web，簡稱WTW)的相關(guān)研究一直是貿(mào)易關(guān)系研究的熱點(diǎn)問題，近年來，Erola等(2012)[20]研究貿(mào)易網(wǎng)絡(luò)演化特征后認(rèn)為，網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的變化是研究“WTW演化→金融危機(jī)”之間關(guān)系的重點(diǎn)。Squartini和Garlaschelli(2012)[21]分析了1950年到2000年間World Trade Network (WTN)中13種三元模體的出現(xiàn)頻次及其演化規(guī)律，發(fā)現(xiàn)相對(duì)于入度、出度，三元模體分析可以提供更多的網(wǎng)絡(luò)結(jié)構(gòu)信息。Maratea等(2016)[22]研究發(fā)現(xiàn)分析WTW中的高階結(jié)關(guān)系(Higher-order ties)非常重要，而“互惠”(Reciprocal)、“三角”(Triadic)等高階結(jié)關(guān)系可以在隨機(jī)模型中提供WTW 演化的更多重要信息。

最近，魏彥杰等(2017)[23]基于高階結(jié)關(guān)系的視角，觀察貿(mào)易網(wǎng)絡(luò)演化與金融危機(jī)的關(guān)系，研究發(fā)現(xiàn)金融危機(jī)周期與WTW中的三角傳遞效應(yīng)有著高度正相關(guān)關(guān)系。他們還指出，三角傳遞效應(yīng)增長不僅預(yù)示著金融危機(jī)的早期信號(hào)，而且也是影響金融危機(jī)的重要因素。本文也將采用類似思路對(duì)相關(guān)紡織紗線貿(mào)易骨干網(wǎng)絡(luò)開展更為細(xì)致的分析。

首先介紹高階結(jié)的相關(guān)概念。眾所周知，網(wǎng)絡(luò)中節(jié)點(diǎn)之間最基本的關(guān)系是連邊，在有向網(wǎng)絡(luò)中，這種最基本關(guān)系的構(gòu)成要素包括兩個(gè)節(jié)點(diǎn)i,j，以及連接著兩個(gè)節(jié)點(diǎn)的一條有向邊eij或eji，記為R={i,j;eij}或R={i,j;eji}，這是網(wǎng)絡(luò)的一個(gè)二元連通子圖。高階結(jié)關(guān)系則是在這種基本關(guān)系的要素集合中添加節(jié)點(diǎn)或邊，形成的網(wǎng)絡(luò)子圖；高階結(jié)一般指的是基本關(guān)系R={i,j;eij}或R={i,j;eji}的要素幾何集合中添加節(jié)點(diǎn)或邊形成的網(wǎng)絡(luò)連通子圖。

按上述定義，二元高階結(jié)只有一種，即互惠關(guān)系；三元高階結(jié)就是Squartini和Garlaschelli(2012)[21]分析的13種三元模體；此外，在國際貿(mào)易網(wǎng)絡(luò)中，3個(gè)及以上國家向同一個(gè)國家進(jìn)口產(chǎn)品或3個(gè)及以上國家向同一個(gè)國家出口產(chǎn)品也是非常重要的貿(mào)易合作模式。因此，本文還將考慮3個(gè)及以上節(jié)點(diǎn)指向同一節(jié)點(diǎn)的輸入關(guān)系、同一節(jié)點(diǎn)指向3個(gè)及以上節(jié)點(diǎn)的輸出關(guān)系這兩種高階結(jié)。

綜上，我們將得到刻畫貿(mào)易網(wǎng)絡(luò)中節(jié)點(diǎn)關(guān)系的16種高階結(jié)。借鑒文獻(xiàn)[21，23]用到的模體分析法，本文將對(duì)2006—2018年間“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)的16種高階結(jié)進(jìn)行統(tǒng)計(jì)分析，以期揭示它們的變化與重要經(jīng)濟(jì)事件間的內(nèi)在關(guān)聯(lián)。本文開展分析的有向網(wǎng)絡(luò)的16種高階結(jié)關(guān)系示意圖如圖5所示：

圖5 本文關(guān)注的16種高階結(jié)關(guān)系示意圖

圖中這16種高階結(jié)關(guān)系具有較強(qiáng)的經(jīng)濟(jì)學(xué)意義：互惠關(guān)系 tie1是一種雙邊貿(mào)易關(guān)系； tie2— tie14共13種三元模體，反映的是最基礎(chǔ)的多邊貿(mào)易關(guān)系，如三角循環(huán)關(guān)系是一種對(duì)稱的三邊貿(mào)易關(guān)系；三角傳遞關(guān)系 tie9則是一種不對(duì)稱的三邊貿(mào)易關(guān)系，其中一方只有輸出，一方只有輸入，還有一方則既有輸出又有輸入；輸入關(guān)系 tie15表示某個(gè)節(jié)點(diǎn)是幾個(gè)節(jié)點(diǎn)形成的集團(tuán)的共同輸出對(duì)象；輸出關(guān)系 tie16表示某個(gè)節(jié)點(diǎn)是幾個(gè)節(jié)點(diǎn)形成的集團(tuán)的共同輸入來源。

貿(mào)易網(wǎng)絡(luò)中各類高階結(jié)數(shù)量的變化，是網(wǎng)絡(luò)內(nèi)部合作變化的直接結(jié)果，反之，通過觀察網(wǎng)絡(luò)內(nèi)各類高階結(jié)數(shù)量的變化，也可以對(duì)網(wǎng)絡(luò)合作的效果做出合理的判斷，推斷各類因素對(duì)網(wǎng)絡(luò)合作等的影響。

下面的表4給出了2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)16種高階結(jié)數(shù)量的統(tǒng)計(jì)結(jié)果。圖6反映了2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)16種高階結(jié)數(shù)量的演變情況。

表4 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易骨干網(wǎng)高階結(jié)統(tǒng)計(jì)表(閾值100萬)

從圖6可以看出，2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易100萬閾值子網(wǎng)16種高階結(jié)關(guān)系的演化趨勢可以分為以下幾類：

1.互惠關(guān)系 tie1，及含有兩個(gè)或三個(gè)互惠關(guān)系的三節(jié)點(diǎn)圈狀高階結(jié) tie13、 tie14，這三種高階結(jié)反映的貿(mào)易關(guān)系的主要特征是互惠。其數(shù)量呈現(xiàn)波浪式上升的趨勢，兩個(gè)比較明顯的轉(zhuǎn)折低點(diǎn)為2008年和2013年，其原因是受2008年世界金融危機(jī)的影響，國際雙邊貿(mào)易關(guān)系整體大幅下跌、隨著金融危機(jī)影響的消退，從2009年起，國際雙邊貿(mào)易關(guān)系開始復(fù)蘇，呈現(xiàn)上升勢頭，但由于貿(mào)易保護(hù)主義抬頭，國際雙邊貿(mào)易關(guān)系在2012年再次下降，BRI正是在這樣的形勢下提出來的，從2013年起“一帶一路”沿線國家紡織紗線貿(mào)易100萬閾值子網(wǎng)中，主要反映雙邊貿(mào)易關(guān)系的這三種高階結(jié)的數(shù)量大幅回升，其后雖然也有回調(diào)的情況，但整體呈現(xiàn)出上升的勢頭，甚至達(dá)到或超過考察期的最高水平。這說明BRI的提出對(duì)“一帶一路”沿線國家紡織紗線貿(mào)易實(shí)質(zhì)性的雙邊互惠合作有顯著的促進(jìn)作用，并讓這種雙邊互惠合作在當(dāng)前貿(mào)易保護(hù)主義抬頭的復(fù)雜環(huán)境下仍然保持了較高的水平。

圖6 2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易100萬閾值子網(wǎng)高階結(jié)關(guān)系演化圖

2. tie2、 tie5、 tie6反映的主要是單向傳遞貿(mào)易關(guān)系。這三種高階結(jié)的數(shù)量在2009年大幅降低，2009年后開始回升，隨著2013年BRI的提出，這三種高階結(jié)的數(shù)量并沒有立即上升，而是經(jīng)過一段時(shí)間的維持甚至小幅下降后又開始提升。這說明BRI的提出，對(duì)這類傳遞貿(mào)易合作關(guān)系的直接提升效果并不顯著。

3.tie3、tie4、tie7、tie8、tie9、tie10、tie11、tie12八種高階結(jié)，反映的是除互惠為主和單向傳遞貿(mào)易為主之外，其他的三邊貿(mào)易合作關(guān)系。這八種高階結(jié)數(shù)量的變化趨勢基本一致，受世界金融危機(jī)的影響，2008年這八種高階結(jié)的數(shù)量下降，但下降幅度不是很大，隨后短暫回升后又開始下降，并在2012年大幅下降，2013年又大幅反彈，反彈的力度相當(dāng)于2008年世界金融危機(jī)后的觸底反彈，此后雖然一度回調(diào)，但又恢復(fù)上升，并維持高位。這說明BRI的提出及時(shí)扭轉(zhuǎn)了“一帶一路”沿線國家紡織紗線三邊貿(mào)易關(guān)系快速下滑的頹勢，大大增強(qiáng)了“一帶一路”沿線各國在紡織紗線貿(mào)易領(lǐng)域加強(qiáng)雙邊、多邊合作的信心和意愿，對(duì)推動(dòng)“一帶一路”沿線國家紡織紗線貿(mào)易發(fā)展起到了關(guān)鍵的作用。

4.BRI提出以來輸入關(guān)系 tie15的數(shù)量大幅上升，直到2017年才有所回落，但仍保持在高位；輸出關(guān)系 tie16的數(shù)量變化則出現(xiàn)了不規(guī)則的鋸齒形，并沒有明顯的趨勢。這說明BRI的提出對(duì)“新進(jìn)口中心國家”的出現(xiàn)產(chǎn)生了推動(dòng)作用，但對(duì)“新出口中心國家”的出現(xiàn)并沒有產(chǎn)生太大的刺激作用，進(jìn)一步表明BRI對(duì)原有貿(mào)易關(guān)系產(chǎn)生的是“助力”而非“沖擊”，這一內(nèi)在機(jī)制引發(fā)的外在經(jīng)濟(jì)表現(xiàn)無疑也是各國積極制定相關(guān)經(jīng)濟(jì)政策、積極參與支持BRI的信心保障。

五、結(jié)論

本文提取聯(lián)合國貿(mào)易數(shù)據(jù)庫2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易數(shù)據(jù)，采取“雙向取大”方法對(duì)數(shù)據(jù)進(jìn)行修正，構(gòu)建了2006—2018年“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)，并以貿(mào)易額100萬美元作為閾值提取骨干子網(wǎng)絡(luò)。通過對(duì)網(wǎng)絡(luò)基本特征演化的分析，以及十六類高階結(jié)關(guān)系數(shù)量變化的分析，得到如下結(jié)論：

1.BRI的提出對(duì)網(wǎng)絡(luò)密度、網(wǎng)絡(luò)節(jié)點(diǎn)平均度、網(wǎng)絡(luò)節(jié)點(diǎn)加權(quán)平均度都有明顯提升，說明BRI顯著促進(jìn)了“一帶一路”沿線國家紡織紗線貿(mào)易合作。

2.2014年以后網(wǎng)絡(luò)密度、網(wǎng)絡(luò)節(jié)點(diǎn)平均度、網(wǎng)絡(luò)節(jié)點(diǎn)加權(quán)平均度有所下降，其中網(wǎng)絡(luò)節(jié)點(diǎn)加權(quán)平均度下降幅度較小，說明一些國家在參與網(wǎng)絡(luò)“一帶一路”沿線國家紡織紗線領(lǐng)域貿(mào)易的熱情有所消退，出現(xiàn)理性回歸，但主要國家間的貿(mào)易合作的力度維持在較高水平。

3.BRI有效提升了“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)中以“互惠”為主要特征的雙邊、三邊貿(mào)易合作關(guān)系，并讓其維持在較高的水平。

4.BRI對(duì)以單向傳遞為主要特征的三邊貿(mào)易合作關(guān)系提升幅度不大，這是因?yàn)锽RI強(qiáng)調(diào)的是互聯(lián)互通，各國之間優(yōu)先實(shí)現(xiàn)互聯(lián)互通。因此，這種單向傳遞關(guān)系數(shù)量沒有顯著提升完全符合預(yù)期。

5.BRI大幅提升了以“三角循環(huán)”“三角傳遞”為主要特征的三邊貿(mào)易合作關(guān)系，BRI推動(dòng)了“一帶一路”沿線國家在紡織紗線貿(mào)易領(lǐng)域的雙邊或多邊合作。在2014年后，這些三邊貿(mào)易合作關(guān)系又出現(xiàn)反復(fù)，一些國家在熱情過后，開始理性回歸，共建“一帶一路”這一浩大的構(gòu)建沿線國家命運(yùn)共同體的經(jīng)濟(jì)工程依舊任重道遠(yuǎn)。

6.BRI提出以后“一帶一路”沿線國家紡織紗線貿(mào)易網(wǎng)絡(luò)中“輸入”高階結(jié)的數(shù)量明顯上升并維持高位，但“輸出”高階結(jié)數(shù)量沒有趨勢性的變化，說明BRI的提出對(duì)新的進(jìn)口中心國家的出現(xiàn)起到了推動(dòng)作用，但對(duì)新的出口中心國家的出現(xiàn)并沒有產(chǎn)生太大的刺激作用，這進(jìn)一步表明BRI對(duì)原有貿(mào)易產(chǎn)生的是“助力”而非“沖擊”。

綜上所述，BRI對(duì)“一帶一路”沿線國家在紡織紗線貿(mào)易領(lǐng)域的合作起到了十分重要的推動(dòng)作用，但同時(shí)要注意到2014—2016年以后，各項(xiàng)積極指標(biāo)都出現(xiàn)了下調(diào)，有些指標(biāo)的下調(diào)幅度還較大，說明共建“一帶一路”是一項(xiàng)長期、復(fù)雜的系統(tǒng)工程，需要各參與國持續(xù)不斷地加以推進(jìn)。

注 釋：

①https://comtrade.un.org/ https://www.yidaiyilu.gov.cn/gbjg/gbgk/77073.htm

②2006年起塞爾維亞和黑山分為兩個(gè)主權(quán)國家—塞爾維亞、黑山分別報(bào)告各自的進(jìn)出口數(shù)據(jù)，此前它們是作為一個(gè)主權(quán)國家—塞爾維亞和黑山一起報(bào)告進(jìn)出口數(shù)據(jù)的。由于自2006年起塞爾維亞和黑山分為兩個(gè)國家，網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)產(chǎn)生變化，黑山是一個(gè)人口只有幾十萬的國家，它單獨(dú)作為一個(gè)國家加入網(wǎng)絡(luò)對(duì)某些指標(biāo)可能會(huì)產(chǎn)生類似離群數(shù)據(jù)的影響，在細(xì)致的網(wǎng)絡(luò)貿(mào)易行為分析中如果想保持貿(mào)易活動(dòng)參與國數(shù)目的一致性，為確保構(gòu)建的貿(mào)易網(wǎng)絡(luò)逐年節(jié)點(diǎn)數(shù)一致，本文只考慮了2006—2018年的網(wǎng)絡(luò)。