長三角城市群貿(mào)易韌性綜合評(píng)價(jià)及時(shí)空演化特征研究

孫玉婷，常志朋，2

（1.安徽工業(yè)大學(xué)商學(xué)院，安徽馬鞍山，243002；2.安徽工業(yè)大學(xué)復(fù)雜系統(tǒng)多學(xué)科管理與控制安徽普通高校重點(diǎn)實(shí)驗(yàn)室，安徽馬鞍山，243002）

如何減少外部不確定性風(fēng)險(xiǎn)對(duì)我國貿(mào)易增長的消極影響，實(shí)現(xiàn)區(qū)域外向型經(jīng)濟(jì)的高質(zhì)量發(fā)展已逐步成為學(xué)界關(guān)注的重點(diǎn)問題。近年來，區(qū)域韌性概念逐步成為不同空間尺度下評(píng)價(jià)區(qū)域發(fā)展質(zhì)量和可持續(xù)性的重要指標(biāo)，為解釋區(qū)域經(jīng)濟(jì)抵御外部風(fēng)險(xiǎn)沖擊和經(jīng)濟(jì)復(fù)蘇能力差異性提供了新視角。

20 世紀(jì)70 年代加拿大生態(tài)學(xué)家Holling 最早將韌性概念引入生態(tài)學(xué)領(lǐng)域，具體表現(xiàn)為系統(tǒng)在受創(chuàng)后能否迅速自我修復(fù)的能力［1］。隨著時(shí)代的發(fā)展，韌性的概念被應(yīng)用到更加廣闊的領(lǐng)域，涉及工程學(xué)、社會(huì)學(xué)以及經(jīng)濟(jì)學(xué)等。賀燦飛等（2019）提出區(qū)域韌性是指在外部擾動(dòng)下，區(qū)域經(jīng)濟(jì)通過各種要素再組織進(jìn)行自我恢復(fù)、更新、轉(zhuǎn)型的能力［2］；除了經(jīng)濟(jì)韌性的研究，部分研究會(huì)測量城市工業(yè)韌性并分析其空間特征［3］；出口貿(mào)易韌性也逐漸成為學(xué)者研究的熱點(diǎn)，賀燦飛等人（2019）指出相關(guān)多樣化的產(chǎn)業(yè)結(jié)構(gòu)可能降低城市-產(chǎn)業(yè)出口韌性［4］。

長三角地區(qū)對(duì)中國出口貿(mào)易的發(fā)展做出了巨大的貢獻(xiàn)，疫情發(fā)生后，長三角地區(qū)在2021 上半年實(shí)現(xiàn)貨物貿(mào)易進(jìn)出口總額65 229 億元，同比增長23.0%，占我國貨物貿(mào)易進(jìn)出口總額的比重為36.1%。長三角地區(qū)貿(mào)易能夠從外部沖擊中迅速恢復(fù)，作為抵御外部風(fēng)險(xiǎn)沖擊的出口貿(mào)易韌性在其中起到了至關(guān)重要的作用。因此，準(zhǔn)確、科學(xué)地測度域內(nèi)城市的出口貿(mào)易韌性并提高長三角地區(qū)貿(mào)易韌性可以為國內(nèi)其他地區(qū)帶來經(jīng)驗(yàn)。目前，關(guān)于區(qū)域貿(mào)易韌性研究成果較少且存在不足：第一，對(duì)于測算貿(mào)易韌性時(shí)，構(gòu)建的指標(biāo)較為單一可能導(dǎo)致測算不精準(zhǔn)；第二，研究對(duì)象多從宏觀角度出發(fā)，較少從城市群的角度出發(fā)分析貿(mào)易韌性；第三，較少文獻(xiàn)分析疫情沖擊對(duì)貿(mào)易韌性的影響。為此，本文將以長三角城市群為研究對(duì)象，構(gòu)建城市貿(mào)易韌性評(píng)價(jià)指標(biāo)體系，采用馬田系統(tǒng)、突變級(jí)數(shù)模型、探索性空間分析等方法來深入分析長三角城市群的貿(mào)易韌性水平及空間特征。

一、研究對(duì)象與指標(biāo)體系構(gòu)建

性；第二階段為突發(fā)風(fēng)險(xiǎn)后的恢復(fù)能力，表現(xiàn)出區(qū)域經(jīng)濟(jì)從沖擊中恢復(fù)的速度和程度；第三階段為區(qū)域經(jīng)濟(jì)完成恢復(fù)后的自我調(diào)整能力，即一個(gè)區(qū)域的經(jīng)濟(jì)在遭受沖擊后會(huì)調(diào)整自身的經(jīng)濟(jì)結(jié)構(gòu)來適應(yīng)新的外部環(huán)境；第四階段為區(qū)域經(jīng)濟(jì)的重新定位能力，即一個(gè)區(qū)域的經(jīng)濟(jì)改變?cè)械慕?jīng)濟(jì)增長模式，創(chuàng)造新的發(fā)展路徑的長期適應(yīng)能力。區(qū)域貿(mào)易韌性是區(qū)域韌性的一個(gè)組成部分，根據(jù)現(xiàn)有研究，本文將長三角城市群貿(mào)易韌性分為抵抗能力、恢復(fù)能力、再組織能力和創(chuàng)新能力四個(gè)維度，構(gòu)建的指標(biāo)體系具體見圖2。

（一）研究對(duì)象與數(shù)據(jù)來源

本文以長三角城市群為研究對(duì)象，根據(jù)《長江三角洲區(qū)域一體化發(fā)展規(guī)劃綱要》，選取27 個(gè)城市，具體如圖1 所示。本文數(shù)據(jù)選取的時(shí)間范圍為2008—2020 年，主要來源于《中國城市統(tǒng)計(jì)年鑒》、各省市統(tǒng)計(jì)年鑒以及各市統(tǒng)計(jì)公報(bào)。

（二）指標(biāo)體系構(gòu)建

Martin［5］將區(qū)域韌性劃分為四個(gè)不同的階段，分別賦予四種不同的能力：第一階段為遭遇風(fēng)險(xiǎn)時(shí)的抵御能力，反映出區(qū)域經(jīng)濟(jì)系統(tǒng)的脆弱性和敏感

根據(jù)貿(mào)易韌性的四維能力，參考丁建軍等人［6］對(duì)經(jīng)濟(jì)韌性作用于內(nèi)外沖擊的過程描述，繪制長三角地區(qū)貿(mào)易系統(tǒng)在遭遇內(nèi)外部沖擊時(shí)的反應(yīng)過程，如圖3 所示，e1、e2、e3分別代表三種類型的風(fēng)險(xiǎn)，T0時(shí)刻是貿(mào)易經(jīng)濟(jì)遭遇風(fēng)險(xiǎn)時(shí)，T0-T1時(shí)段是貿(mào)易經(jīng)濟(jì)遭遇風(fēng)險(xiǎn)后，抵抗與再組織能力發(fā)揮作用，減緩風(fēng)險(xiǎn)帶來的負(fù)面影響，T1-T2時(shí)段是貿(mào)易韌性的恢復(fù)與創(chuàng)新能力幫助長三角城市貿(mào)易經(jīng)濟(jì)恢復(fù)原有水平或者實(shí)現(xiàn)新的發(fā)展路徑。對(duì)于不同貿(mào)易韌性水平的城市在受到?jīng)_擊時(shí)，導(dǎo)致對(duì)內(nèi)外部沖擊的反應(yīng)程度存在差異，主要有三種結(jié)果：實(shí)現(xiàn)新的發(fā)展路徑、恢復(fù)原有水平但未實(shí)現(xiàn)新的“路徑創(chuàng)造”、貿(mào)易經(jīng)濟(jì)處于持續(xù)衰退中。

二、長三角城市群貿(mào)易韌性綜合評(píng)價(jià)

（一）基于馬田系統(tǒng)的指標(biāo)權(quán)重測度

本文采用馬田系統(tǒng)計(jì)算指標(biāo)權(quán)重，馬田系統(tǒng)是20 世紀(jì)90 年代初日本著名質(zhì)量工程學(xué)家田口玄一博士提出的一種多元數(shù)據(jù)模式識(shí)別方法，它將馬氏距離和田口方法相結(jié)合［7］。該方法的一個(gè)重要功能是可以利用基于馬氏距離的信噪比測度屬性集在分類過程中的重要程度，即屬性集對(duì)MTS 能夠正確判斷類別的貢獻(xiàn)［8］。具體步驟如下：

1.構(gòu)建異常樣本矩陣

設(shè)Yi(i= 1,2,…,m) 為城市貿(mào)易韌性水平，xp(p= 1,2,…,n)為城市韌性指標(biāo)，則城市貿(mào)易韌性綜合評(píng)價(jià)問題可用矩陣形式A［9］表示為：

將城市貿(mào)易韌性評(píng)估矩陣A設(shè)定為異常樣本，根據(jù)A構(gòu)建正常矩陣B如下：

2.計(jì)算馬氏距離

對(duì)矩陣A進(jìn)行標(biāo)準(zhǔn)化處理，(xp)和(xp)分別是正常樣本xp的均值和標(biāo)準(zhǔn)差，zi(xp)是異常樣本標(biāo)準(zhǔn)化后的值，具體公式如下：

接下來，利用標(biāo)準(zhǔn)化后的數(shù)值計(jì)算異常樣本的馬氏距離，其中，p表示指標(biāo)總數(shù)量，Zi= (zi(x1)，zi(x2)，…，zi(xn))T，i= 1，2，…，m，ziT為zi的轉(zhuǎn)置矩陣，S為正常樣本的協(xié)方差矩陣，S+是S的廣義逆矩陣，具體公式如下：

3.計(jì)算貿(mào)易韌性指標(biāo)權(quán)重

使用正交表和信噪比來確定貿(mào)易韌性指標(biāo)的權(quán)重。本文在實(shí)驗(yàn)中采用二水平的正交表Lq(2p)測度貿(mào)易韌性指標(biāo)的影響對(duì)區(qū)域貿(mào)易的影響程度，其中q表示實(shí)驗(yàn)次數(shù)，p表示參與試驗(yàn)的特征變量個(gè)數(shù)，將p個(gè)變量分配到前p列中，正交表的每一行都會(huì)生成一個(gè)基準(zhǔn)空間，即對(duì)應(yīng)于每一行中水平為“1”的特征變量作為多維變量的一維參與基準(zhǔn)空間的生成，其中“1”表示使用該變量，“2”表示不使用該變量。一般常用的二水平正交表有L4( 23)、L8( 27)、L16( 215)、L20( 219)、L32(231)，本文共有16 個(gè)指標(biāo)，所以選用L20( 219)，將16 個(gè)指標(biāo)分布在前16列，安排正交實(shí)驗(yàn)。

信噪比是有用信息（信號(hào)）同有害信息（干擾）的比值，主要分為望目特性信噪比、望小特性信噪比以及望大特性信噪比，本文采用望大特性信噪比，公式如下：

其中，ηq表示第q次試驗(yàn)的信噪比，表示異常樣本的馬氏距離，m為長三角城市群數(shù)量。以2020年為例，將正交表和信噪比列為表1 所示。

信息增益（Gain）反映了信噪比中變量的重要程度，通常利用信噪比計(jì)算得出，計(jì)算公式如下：

其中，+表示貿(mào)易韌性指標(biāo)參與所有試驗(yàn)的信噪比均值，-表示貿(mào)易韌性指標(biāo)未參與試驗(yàn)的信噪比均值，信息增益圖如圖4 表示：

圖4 信息增益圖

接下來，計(jì)算城市貿(mào)易韌性指標(biāo)的權(quán)重，因?yàn)闄?quán)重w(xp) ∈(0,1)，而信息增益Gainp∈(-∞, +∞)，所以使用Sigmoid 函數(shù)進(jìn)行轉(zhuǎn)換，σ∈(0, +∞)為函數(shù)的權(quán)值，當(dāng)σ→+∞時(shí)，Sigmoid 函數(shù)圖像越接近跳躍式函數(shù)圖像，σ越小Sigmoid 函數(shù)曲線越平坦［10-11］。公式如下：

最后，對(duì)區(qū)域貿(mào)易韌性評(píng)價(jià)指標(biāo)的轉(zhuǎn)化結(jié)果進(jìn)行歸一化處理φ(xp)歸一化處理，得到指標(biāo)權(quán)重為：

由于計(jì)算出的是每個(gè)年份的指標(biāo)權(quán)重，需要得出統(tǒng)一的權(quán)重值，引入時(shí)間權(quán)重向量，設(shè)H=(h1,h2,…,hk)為時(shí)間權(quán)重向量，，ht表示第t個(gè)時(shí)間段的權(quán)重，時(shí)間權(quán)重體現(xiàn)了決策過程中對(duì)不同時(shí)間段樣本的重視程度［12］。根據(jù)“厚古薄今”思想，越靠近當(dāng)今的信息越重要，越有利于評(píng)價(jià)結(jié)果［12-13］，所以賦予最近時(shí)間段的權(quán)重較大，賦予較遠(yuǎn)時(shí)間段的權(quán)重較小，得出最終的權(quán)重值如圖5所示。

圖5 貿(mào)易韌性指標(biāo)體系權(quán)重圖

（二）基于突變級(jí)數(shù)貿(mào)易韌性水平測度

本文選用突變級(jí)數(shù)模型來計(jì)算貿(mào)易韌性指標(biāo)具體值。突變理論是法國數(shù)學(xué)家Rene Thom 所創(chuàng)立，用于研究不連續(xù)變化的現(xiàn)象，突變級(jí)數(shù)模型是在突變理論的基礎(chǔ)上發(fā)展起來的一種綜合評(píng)價(jià)方法［14］。根據(jù)突變級(jí)數(shù)法的建模思想，首先以突變理論和模糊數(shù)學(xué)相結(jié)合生成突變模糊隸屬函數(shù)，再由歸一公式進(jìn)行綜合量化運(yùn)算，最終求出總的隸屬函數(shù)［15］。具體步驟如下：

1.確定突變系統(tǒng)函數(shù)類型

由于本文貿(mào)易韌性指標(biāo)體系的每一層級(jí)的指標(biāo)有4 個(gè)，根據(jù)突變系統(tǒng)選擇規(guī)則，確定選擇蝴蝶突變系統(tǒng)。對(duì)蝴蝶突變系統(tǒng)的勢函數(shù)進(jìn)行一次和二次求導(dǎo)得到f'(x)和f″(x)，通過聯(lián)立f'(x) = 0 和f″(x) = 0 消去x，求得突變模型的分歧集方程，進(jìn)而推導(dǎo)出歸一分歧方程，具體公式如表2 所示。

表2 常見突變函數(shù)類型及相關(guān)內(nèi)容

2.取值準(zhǔn)則

計(jì)算各層指標(biāo)的相關(guān)系數(shù)矩陣，若同一指標(biāo)層變量存在強(qiáng)相關(guān)關(guān)系，那么取變量的平均值作為指標(biāo)值，若同一指標(biāo)層的變量存在弱相關(guān)關(guān)系，則按照“大中取小”的原則取變量的最

小值作為指標(biāo)值，具體韌性指標(biāo)值如圖6 所示。從表中可以看出，上海、江蘇及浙江部分城市貿(mào)易韌性值的變化呈現(xiàn)“動(dòng)態(tài)均衡”的狀態(tài)，貿(mào)易韌性有微小波動(dòng)，但總體持平且處于高位狀態(tài)；安徽和江浙部分城市的貿(mào)易韌性總體上處于上升的狀態(tài)；2018 年中美貿(mào)易戰(zhàn)以及2019 年底新冠疫情暴發(fā)對(duì)城市群貿(mào)易韌性產(chǎn)生負(fù)面影響，在連續(xù)兩年危機(jī)沖擊下，總體來看，大部分城市貿(mào)易韌性先降后升，小部分城市貿(mào)易韌性值先升后降或者持續(xù)下降，其中，宣城和馬鞍山兩市韌性值保持持續(xù)上升。

圖6 長三角城市群貿(mào)易韌性水平

除測算韌性總值外，本文從不同維度分析長三角地區(qū)貿(mào)易韌性總體變化情況（圖7）。長三角城市群總體韌性值呈現(xiàn)緩慢平穩(wěn)上升，到2020 年已處于中高韌性水平；貿(mào)易韌性不同維度的能力呈現(xiàn)不同演化特征，再組織能力數(shù)值最高，說明長三角城市群貿(mào)易調(diào)整功能和轉(zhuǎn)型的能力最強(qiáng)，能夠較好地適應(yīng)經(jīng)濟(jì)環(huán)境的變化；恢復(fù)能力數(shù)值最低，表明長三角城市群從危機(jī)中恢復(fù)原有水平的能力較弱；抵抗能力和創(chuàng)新能力差距不大且數(shù)值較高，表明長三角城市群抵抗風(fēng)險(xiǎn)和創(chuàng)造新的發(fā)展路徑的能力較強(qiáng)。

圖7 長三角城市群四個(gè)維度貿(mào)易韌性水平分布圖

三、長三角城市群貿(mào)易韌性時(shí)空演化分析

（一）時(shí)序演化特征

為分析長三角城市群貿(mào)易韌性的時(shí)序演化特征，選用核密度函數(shù)進(jìn)行估計(jì)。核密度估計(jì)是在概率論中用來估計(jì)隨機(jī)變量的密度函數(shù)，屬于非參數(shù)檢驗(yàn)方法之一。假設(shè)一組連續(xù)型隨機(jī)變量X的密度函數(shù)為f(x)，則公式如下：

其中，n為城市個(gè)數(shù)，hn為影響核密度曲線平滑程度的帶寬，K(?)為核函數(shù)，ui為i市的貿(mào)易韌性值。

由于2008 年、2018 年和2019 年分別發(fā)生了金融危機(jī)、中美貿(mào)易戰(zhàn)和新冠疫情，所以選取了2008年、2009 年、2014 年、2018 年、2019 年和2020 年六個(gè)年份，以用于更好地分析危機(jī)發(fā)生前后貿(mào)易韌性的變化（圖8）。結(jié)合圖1，總體上，從2008—2020年，核密度函數(shù)中心由2008 年集中于中低韌性到2020 年集中于中高韌性，說明自2008 年以來各城市的貿(mào)易韌性水平提升較為顯著；從曲線形態(tài)上來看，開口寬度逐漸縮小，表明長三角城市群各城市貿(mào)易韌性水平差距不斷縮小。從左圖得出，金融危機(jī)發(fā)生后，核密度函數(shù)由2008 年的尖峰型向2009年寬峰型發(fā)展變化，貿(mào)易韌性高峰值有所下降，且曲線開口寬度變大，城市貿(mào)易韌性水平差距變大；觀察右圖可得，2018 年中美貿(mào)易戰(zhàn)的發(fā)生使得2019 年核密度函數(shù)中心向左移動(dòng)，略傾向于中低韌性水平，而2019 年新冠疫情的沖擊，核密度函數(shù)中心向中高韌性方向移動(dòng)且轉(zhuǎn)變?yōu)榧夥逍?，說明對(duì)大部分城市而言，貿(mào)易韌性增強(qiáng)，高峰值有所上升。

圖8 長三角城市群貿(mào)易韌性核密度分布圖

（二）空間分異格局

根據(jù)Arcgis 自然斷點(diǎn)法，將長三角城市群的貿(mào)易韌性評(píng)價(jià)值劃分為低韌性（0.284 321～0.395 403）、中低韌性（0.395 404～0.478 704）、中高韌性（0.478 705～0.621 219）、高韌性地區(qū)（0.621 220～0.844 737）四個(gè)等級(jí)（圖9）。從選取的部分年份來看，2008 年城市群貿(mào)易韌性水平集中于低韌性和中低韌性地區(qū)，2020 年貿(mào)易韌性水平已轉(zhuǎn)變?yōu)榧杏谥懈唔g性地區(qū)。其中，2008 年金融危機(jī)發(fā)生后，貿(mào)易韌性增速放緩，2009 年整體貿(mào)易韌性值與2008 年持平，部分年份出現(xiàn)輕微下降；2018 年受中美貿(mào)易戰(zhàn)疊加新冠疫情影響，貿(mào)易韌性逆勢增長，其中南京由中高韌性區(qū)上升為高韌性地區(qū)，揚(yáng)州、泰州、鹽城三市貿(mào)易韌性提高，成為中高韌性城市；從區(qū)域整體來看，貿(mào)易韌性高值區(qū)集中于上海、蘇南、浙北等地區(qū)，安徽省大部分城市處于低韌性和中低韌性水平。

圖9 長三角城市群貿(mào)易韌性總體變化圖

（三）空間關(guān)聯(lián)性分析

1.全局空間自相關(guān)

空間自相關(guān)是用來度量空間單元屬性值的平均集聚程度，分為全局空間自相關(guān)和局部空間自相關(guān)，全局空間自相關(guān)是檢驗(yàn)長三角城市群整體空間的關(guān)聯(lián)性，主要通過Moran' I 指數(shù)進(jìn)行分析，公式如下：

其中，n為城市個(gè)數(shù)，wij為空間權(quán)重矩陣，xi和xj為第i和第j個(gè)城市的觀察值。

計(jì)算2003—2019 年長三角城市群貿(mào)易韌性水平的全局Moran' I 指數(shù)，從表3 可以看出，Moran' I指數(shù)均為正值，且通過了p＜0.01 水平的顯著性檢驗(yàn)，表明長三角城市群貿(mào)易韌性存在空間自相關(guān)性。Moran' I 指數(shù)呈現(xiàn)先升后降再升的趨勢，其中在2010 年達(dá)到峰值，說明長三角城市群貿(mào)易韌性的空間自相關(guān)性表現(xiàn)為“先增強(qiáng)后減弱再增強(qiáng)”的趨勢，2010 年后，空間異質(zhì)性逐漸增強(qiáng)，而2020 年空間異質(zhì)性開始減弱。

表3 長三角城市群貿(mào)易韌性的全局Moran'I 指數(shù)

2.局部空間自相關(guān)

局部空間自相關(guān)通常用局部Moran' I 指數(shù)即LISA 指數(shù)來測度，公式如下：

LISA 指數(shù)反映了長三角城市群貿(mào)易韌性局部地區(qū)空間關(guān)聯(lián)特征，從圖10 可以看出，高-高型地區(qū)由2008 年上海、蘇州、南通三市轉(zhuǎn)變?yōu)?020 年的上海、蘇州、南通、嘉興、湖州五市，說明上海都市圈具有較高貿(mào)易韌性水平；從2008 年到2020 年，低-低型地區(qū)總體在增加，且集中于安徽省的安慶、池州、銅陵、蕪湖等城市。低—高型區(qū)域主要在嘉興、湖州南通等地；高—低型區(qū)域在所示年份中集中合肥、南京兩市，表明作為省會(huì)城市的貿(mào)易韌性水平高于周邊城市，對(duì)周邊城市有輻射作用。

圖10 長三角城市群貿(mào)易韌性局部自相關(guān)分布圖

3.空間重心分析

重心遷移可以反映長三角城市群的整體性演變特征，重心的移動(dòng)方向和距離可以揭示其貿(mào)易韌性的變化趨勢。計(jì)算公式如下，式中（ci，di）表示各城市的坐標(biāo)，ui為貿(mào)易韌性值。

如圖11，分別計(jì)算部分年份的長三角地區(qū)重心距離，其中心變動(dòng)范圍位于東經(jīng)119.585°～119.663°和北緯30.844°～30.893°，以上年份均位于湖州市的西北部，長三角的中東部地區(qū)，這與該地區(qū)聚集貿(mào)易韌性較強(qiáng)的中心城市和省會(huì)城市有關(guān)。從重心的移動(dòng)軌跡來看，2008—2020 年間，貿(mào)易韌性整體轉(zhuǎn)移方向?yàn)闁|北-西南-西北，其原因在于江蘇省重視蘇北地區(qū)的經(jīng)濟(jì)發(fā)展，不斷給予政策支持，使得蘇北貿(mào)易韌性顯著提升，從而重心向東北轉(zhuǎn)移，2009 年以后，隨著安徽省積極融入長三角經(jīng)濟(jì)區(qū)，使得省內(nèi)合肥、蕪湖等城市的貿(mào)易韌性不斷增強(qiáng)，使得貿(mào)易韌性重心由東北轉(zhuǎn)向西南，2019 年長三角一體化全面推進(jìn)，帶動(dòng)安徽北部城市發(fā)展，重心向西北地區(qū)遷移。同時(shí)，從重心的移動(dòng)幅度來看，2008—2014 年移動(dòng)范圍最大，2014—2019 年移動(dòng)范圍最小，2019—2020 年移動(dòng)范圍又開始增大，呈現(xiàn)“初始較大—中間較小—末端增大”的變化表征。

圖11 長三角城市群貿(mào)易韌性重心遷移圖

四、結(jié)論與啟示

基于2008—2020 年長三角城市群數(shù)據(jù)，運(yùn)用馬田系統(tǒng)和突變級(jí)數(shù)模型測度其貿(mào)易韌性水平，并基于核密度函數(shù)、探索性空間分析方法及重心分析方法描繪貿(mào)易韌性水平的時(shí)空演化特征，得到的結(jié)論如下：

第一，從不同維度看，貿(mào)易韌性不同維度的能力呈現(xiàn)不同演化特征，總體呈現(xiàn)緩慢上升的趨勢，其中再組織能力較強(qiáng)。第二，從時(shí)間變化看，2008—2020 年，長三角城市群貿(mào)易韌性水平總體呈現(xiàn)波動(dòng)上升趨勢，核密度函數(shù)中心向右移動(dòng)，說明2008 年以來各城市的貿(mào)易韌性水平提升較為顯著，其中2008 年和2018 年及2019 年外部風(fēng)險(xiǎn)的沖擊對(duì)貿(mào)易造成一定的危害，部分城市的貿(mào)易韌性水平出現(xiàn)輕微下降。第三，從空間格局看，區(qū)域貿(mào)易韌性整體存在顯著的空間相關(guān)性，且空間分布不平衡，其中，以上海為中心的蘇南、浙北等周邊城市的整體貿(mào)易韌性水平較高，合肥、南京等市對(duì)周邊城市的貿(mào)易韌性拉動(dòng)作用有待進(jìn)一步增強(qiáng)，安慶、池州、銅陵、蕪湖等城市由于自身空間溢出效應(yīng)偏弱，仍處于相對(duì)低水平的聚集過程。但是，隨著安徽省不斷積極融入長三角經(jīng)濟(jì)區(qū)，整體貿(mào)易韌性重心由東北向西南再向西北轉(zhuǎn)移。

根據(jù)以上結(jié)論，建議：從提升貿(mào)易韌性的恢復(fù)能力和創(chuàng)新能力出發(fā)，需著力提升長三角地區(qū)經(jīng)濟(jì)實(shí)力，加快上海、南京、蘇州等經(jīng)濟(jì)發(fā)達(dá)城市對(duì)周邊城市的經(jīng)濟(jì)輻射作用，同時(shí)扶持安徽省經(jīng)濟(jì)較弱的城市，以調(diào)整區(qū)域非均衡的空間格局；加強(qiáng)科研經(jīng)費(fèi)投入，強(qiáng)化科技人才支撐，增強(qiáng)貿(mào)易創(chuàng)新能力，促進(jìn)長三角城市高質(zhì)量開放水平提升；不斷優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，充分發(fā)揮第二產(chǎn)業(yè)和第三產(chǎn)業(yè)對(duì)貿(mào)易的推動(dòng)作用。要以市場為導(dǎo)向，加快轉(zhuǎn)變外貿(mào)發(fā)展方式，著力優(yōu)化出口主體結(jié)構(gòu)、商品結(jié)構(gòu)、市場結(jié)構(gòu)和貿(mào)易結(jié)構(gòu)。