徐順聰，楊永斌，張世偉，張 磊

中國(guó)人民警察大學(xué)，河北 廊坊 065000

一、引言

近年來，全球恐怖主義風(fēng)險(xiǎn)仍處于高發(fā)態(tài)勢(shì)，恐怖主義問題引發(fā)全世界的廣泛關(guān)注。巴基斯坦是“一帶一路”沿線的重要國(guó)家，“中巴經(jīng)濟(jì)走廊”是“六廊六路”的重要組成部分，我國(guó)在巴項(xiàng)目眾多，經(jīng)濟(jì)及戰(zhàn)略意義重大，巴基斯坦境內(nèi)恐怖襲擊事件也同時(shí)威脅我國(guó)在巴項(xiàng)目安全。2021 年7 月14 日，由中方企業(yè)承建的開普省達(dá)蘇水電站項(xiàng)目出勤班車在赴施工現(xiàn)場(chǎng)途中遭遇襲擊爆炸，造成包括9 名中方人員和3名巴方人員在內(nèi)的12人死亡。2021年8月20日，一名自殺式襲擊者在俾路支省瓜達(dá)爾地區(qū)襲擊載有中國(guó)公民的車隊(duì)，造成2 名兒童死亡、3 人受傷。巴基斯坦境內(nèi)恐怖主義風(fēng)險(xiǎn)較為突出，根據(jù)《2023年全球恐怖主義指數(shù)》，巴基斯坦的全球恐怖主義指數(shù)排名已經(jīng)上升到第六名，是恐怖主義較為嚴(yán)重和活躍的地區(qū)［1］。研究巴基斯坦境內(nèi)恐怖襲擊事件對(duì)保護(hù)我國(guó)在巴項(xiàng)目安全，增強(qiáng)我國(guó)海外利益保護(hù)具有重要意義。

當(dāng)前，對(duì)于恐怖襲擊事件的預(yù)防預(yù)警工作已有一定的研究成果。hao 等為揭示東南亞恐怖襲擊的分布規(guī)律提出基于反距離加權(quán)和基于多標(biāo)簽K-最近鄰的評(píng)估模型［2］。師原兵在城市軌道交通領(lǐng)域中應(yīng)用大數(shù)據(jù)技術(shù)來加強(qiáng)對(duì)數(shù)據(jù)的分析與研判，并提出從多個(gè)方面提升反恐防范和應(yīng)急處置能力，有效打擊恐怖主義犯罪并維護(hù)國(guó)家安全［3］。鄧世成等采用投影尋蹤評(píng)價(jià)模型和灰色新陳代謝模型對(duì)全球恐怖襲擊事件的危害程度進(jìn)行量化分級(jí)，預(yù)測(cè)了未來5年全球恐怖活動(dòng)的蔓延趨勢(shì)［4］。羅瀾峻等利用長(zhǎng)短時(shí)記憶神經(jīng)網(wǎng)絡(luò)模型對(duì)恐怖襲擊事件的發(fā)生時(shí)間進(jìn)行預(yù)測(cè)，結(jié)果表明該模型能夠較準(zhǔn)確預(yù)測(cè)出短時(shí)內(nèi)恐怖襲擊事件的發(fā)生時(shí)間［5］。這些文獻(xiàn)大多研究恐怖襲擊事件的動(dòng)機(jī)和成因，很少與地理空間信息相結(jié)合，通過可視化的方法分析恐怖襲擊事件。

在恐怖襲擊時(shí)空演變分析方面，國(guó)內(nèi)外已經(jīng)有許多學(xué)者對(duì)此進(jìn)行了研究。王誠(chéng)聰?shù)炔捎脽狳c(diǎn)分析和層次聚類分析的方法對(duì)全球恐怖襲擊造成的傷亡人數(shù)和事件頻次進(jìn)行時(shí)空演變特征分析［6］。姜春艷等對(duì)全球恐怖襲擊事件進(jìn)行探索，根據(jù)恐怖襲擊事件的時(shí)間序列以及空間特征指出恐怖襲擊活動(dòng)主要集中在南亞、中東和北非、撒哈拉以南的非洲及東南亞地區(qū)［7］。Hao 等采用地理空間統(tǒng)計(jì)方法并結(jié)合隨機(jī)森林算法，對(duì)中南半島恐怖襲擊的時(shí)空演變規(guī)律進(jìn)行分析［8］。同樣，陳偉等將核密度分析用于恐怖襲擊聚集度的可視化以及標(biāo)準(zhǔn)差橢圓用于分析恐怖襲擊的空間演變特征［9］。郝蒙蒙等利用空間統(tǒng)計(jì)方法和核密度分析的方法分析中南半島恐怖襲擊事件的時(shí)空演變特征［10］。李繼萌等運(yùn)用統(tǒng)計(jì)分析法和核密度估計(jì)法指出“一帶一路”沿線國(guó)家遭受的恐怖襲擊事件在空間上呈現(xiàn)出高度聚集特征并且向四周擴(kuò)散，國(guó)家交界處是恐怖襲擊事件多發(fā)地帶［11］。

當(dāng)前對(duì)巴基斯坦境內(nèi)恐怖襲擊的研究除集中于定性研究外，也開始出現(xiàn)利用時(shí)空統(tǒng)計(jì)技術(shù)對(duì)其分布特征的研究［12］。但在前期研究中將恐怖襲擊事件等同對(duì)待，更多的是以恐怖主義襲擊的次數(shù)作為分析對(duì)象，對(duì)不同恐怖襲擊事件之間危害性差異考慮較少。本文選用2011—2020 年巴基斯坦境內(nèi)恐怖襲擊數(shù)據(jù)，通過構(gòu)建衡量不同恐怖襲擊事件危害性的評(píng)價(jià)指數(shù)——威脅度，建立威脅度評(píng)價(jià)指標(biāo)體系及量化模型（如圖1 所示），對(duì)恐怖襲擊事件進(jìn)行威脅度分級(jí)。隨后，采用核密度分析、空間自相關(guān)分析以及新興時(shí)空熱點(diǎn)分析等方法研究巴基斯坦境內(nèi)恐怖襲擊事件的危害程度及空間演變特征，以期對(duì)未來預(yù)防恐怖襲擊提供重要的參考價(jià)值。

圖1 威脅度評(píng)價(jià)指標(biāo)體系及量化模型

二、研究區(qū)域與研究方法

（一）研究區(qū)域概況

巴基斯坦，全稱巴基斯坦伊斯蘭共和國(guó)，位于南亞次大陸的西北部，地理范圍是北緯24°37′至37°06′，東經(jīng)61°33′至75°31′。西部與伊朗接壤，南部和東南部與印度接壤，東部與阿富汗毗鄰，北部與中國(guó)相鄰，國(guó)土面積（不含巴控克什米爾地區(qū)）為79.6 萬平方公里。巴基斯坦的行政區(qū)劃包括：四個(gè)省份，分別是旁遮普省、開伯爾-普赫?qǐng)D赫瓦?。ㄒ韵潞?jiǎn)稱開普?。①侣分『托诺率?；兩個(gè)行政區(qū)域，分別為自由克什米爾和吉爾吉特-巴爾蒂斯坦、伊斯蘭堡首都地區(qū)。

（二）研究方法

1. 威脅度分級(jí)

對(duì)2011—2020 年巴基斯坦境內(nèi)恐怖襲擊事件的數(shù)據(jù)進(jìn)行綜合分析，以地理空間統(tǒng)計(jì)學(xué)的視角，建立一個(gè)完善的指標(biāo)體系和分析模型，以期獲得更準(zhǔn)確的結(jié)論。參照《生產(chǎn)安全事故報(bào)告和調(diào)查處理?xiàng)l例》將恐怖襲擊事件劃分為四個(gè)等級(jí)，對(duì)數(shù)據(jù)進(jìn)行離散化處理，將恐怖襲擊事件威脅度評(píng)價(jià)指標(biāo)按1～4進(jìn)行賦值，通過K-Means 聚類分析的方法對(duì)事件進(jìn)行分級(jí)處理。

2. 核密度分析

核密度估計(jì)法（Kernel Density Estimation，KDE）認(rèn)為地理事件可能發(fā)生在空間的任何位置，但不同位置上事件發(fā)生的概率不同。相對(duì)于樣方計(jì)數(shù)法，KDE更適合使用可視化方法來表示分布模式［13］。研究中使用高斯核密度函數(shù)，設(shè)X1，X2，…，Xn是從分布密度函數(shù)f中抽取的獨(dú)立同分布樣本，則估計(jì)f在某點(diǎn)x處的密度函數(shù)為：

式中：n為樣點(diǎn)個(gè)數(shù)，k(·)為核函數(shù)，h為帶寬，(x-Xi)為點(diǎn)x到事件Xi處的距離。

3. 空間自相關(guān)

在對(duì)空間進(jìn)行計(jì)量時(shí)，需考慮空間與空間之間是否存在相關(guān)性。空間相關(guān)性是指相鄰位置上區(qū)域事物發(fā)展水平相似。若是高值與高值、低值與低值相鄰，則稱為“正空間自相關(guān)”。反之，若是高值與低值相鄰，則稱為“負(fù)空間自相關(guān)”。如果高值與低值完全隨機(jī)分布，則不存在空間自相關(guān)。全局Moran’s I指數(shù)法計(jì)算公式如下［14］：

式中：S0為所有空間權(quán)重之和，S0可以表示為S0=為要素i的屬性與其平均值(x-)的偏差；n為樣本總數(shù)；wij為要素i和j之間的空間權(quán)重。

全局性Moran’s I 指數(shù)法檢驗(yàn)的標(biāo)準(zhǔn)化統(tǒng)計(jì)量為：

式中：E(I)為其期望值，V(I)為其方差值。

4. 熱點(diǎn)識(shí)別法

可通過Gi*指數(shù)法對(duì)巴基斯坦境內(nèi)恐怖襲擊局部情況進(jìn)行分析，Gi*指數(shù)可以比較準(zhǔn)確地探測(cè)出局部空間聚集區(qū)域［15］。其公式如下［16］：

式中：xj為要素j的屬性值，wij為要素i和j之間的空間權(quán)重，n為要素總數(shù)，期望值

熱點(diǎn)分析各模式如表1所示：

表1 熱點(diǎn)分析各模式定義

5.Mann-Kendall趨勢(shì)分析

Mann-Kendall 趨勢(shì)分析方法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，它可用于對(duì)獨(dú)立的時(shí)空條柱序列進(jìn)行測(cè)試。對(duì)于由n個(gè)樣本組成的時(shí)間序列，趨勢(shì)檢驗(yàn)統(tǒng)計(jì)量如下［17］：

式中：sgn(·)是符號(hào)函數(shù)，若xj＞xi取1，若xj=xi取0，若xj＜xi取-1。

當(dāng)S服從正態(tài)分布時(shí)，其方差公式如下：

標(biāo)準(zhǔn)化統(tǒng)計(jì)量Z計(jì)算公式如下：

根據(jù)Z得分對(duì)恐怖襲擊發(fā)生率的冷熱點(diǎn)趨勢(shì)的顯著性進(jìn)行歸類，分為新增熱點(diǎn)、連續(xù)熱點(diǎn)、加強(qiáng)熱點(diǎn)等冷熱點(diǎn)模式（如表2所示）。

表2 新興時(shí)空熱點(diǎn)分析各模式定義（部分）

三、巴基斯坦境內(nèi)恐怖襲擊事件威脅度評(píng)價(jià)

采用威脅度評(píng)價(jià)對(duì)恐怖襲擊事件進(jìn)行量化分析，有助于準(zhǔn)確地預(yù)測(cè)恐怖襲擊事件的發(fā)生概率。通過量化評(píng)估，還可以確定哪些地區(qū)、哪些群體、哪些場(chǎng)所面臨更大的威脅，從而有針對(duì)性地配置防恐資源，提高反恐效率。

（一）數(shù)據(jù)來源及處理

1. 數(shù)據(jù)來源

本文的數(shù)據(jù)來源為全球恐怖主義數(shù)據(jù)庫(kù)（Global Terrorism Database，GTD），GTD 是由馬里蘭大學(xué)開發(fā)維護(hù)的一個(gè)開源數(shù)據(jù)庫(kù)，下載地址為https：//www. start. umd. edu/gtd/，記錄了1970—2020年全球各地恐怖襲擊事件的詳細(xì)信息［18-19］。其中包括恐怖襲擊事件發(fā)生的時(shí)間、國(guó)家、經(jīng)緯度、襲擊組織以及造成的傷亡人數(shù)等內(nèi)容。本文選取巴基斯坦2011—2020年的恐怖襲擊數(shù)據(jù)。

2. 數(shù)據(jù)預(yù)處理

（1）缺失值處理。在排除經(jīng)緯度坐標(biāo)缺失和傷亡人數(shù)有缺失的數(shù)據(jù)后，共得到10 639 條有效恐怖襲擊數(shù)據(jù)。已有研究表明，當(dāng)缺失值小于5%時(shí)，使用隨機(jī)森林填補(bǔ)法效果較優(yōu)，并且隨機(jī)森林填補(bǔ)法不僅在軟件操作上簡(jiǎn)便，還對(duì)數(shù)據(jù)結(jié)構(gòu)和分布的要求較為寬松，在反映真實(shí)情況方面也較為準(zhǔn)確［20］。故選擇使用sklearn庫(kù)中的隨機(jī)森林算法來對(duì)有缺失值的變量進(jìn)行填充。

（2）無量綱處理。在該研究中使用極值處理法消除各指標(biāo)單位和數(shù)值間的系統(tǒng)誤差，將原始數(shù)據(jù)處理成0～1范圍內(nèi)的數(shù)據(jù)。公式如下［21］：

式中：xij(i= 1，2…n；j= 1，2…m)表示第i個(gè)被評(píng)價(jià)單元在第j項(xiàng)指標(biāo)上的觀測(cè)值是經(jīng)過無量綱化處理后的標(biāo)準(zhǔn)觀測(cè)值。

（二）評(píng)價(jià)指標(biāo)體系的構(gòu)建

本研究根據(jù)GTD 數(shù)據(jù)庫(kù)信息，選取7 個(gè)指標(biāo)（如圖2 所示），包括死亡人數(shù)、受傷人數(shù)、財(cái)產(chǎn)損失、目標(biāo)類型、襲擊地點(diǎn)、攻擊類型和武器類型。其中，定量指標(biāo)有死亡人數(shù)、受傷人數(shù)和財(cái)產(chǎn)損失。定性指標(biāo)有目標(biāo)類型、襲擊地點(diǎn)、攻擊類型和武器類型。

圖2 巴基斯坦境內(nèi)恐怖襲擊事件威脅度指標(biāo)體系

（三）指標(biāo)體系權(quán)重的確定

在對(duì)恐怖襲擊事件威脅度進(jìn)行計(jì)算時(shí)，通過等級(jí)賦分將定性指標(biāo)量化，研究中分為4 個(gè)等級(jí)（如表3所示）［22］。

確定指標(biāo)體系權(quán)重的方法可分為兩類：一類是主觀賦權(quán)方法，它的缺點(diǎn)在于過于依賴專家自身的經(jīng)驗(yàn)與知識(shí)，主觀性較強(qiáng)；另一類是客觀賦權(quán)方法，缺點(diǎn)在于過分依賴樣本數(shù)據(jù)，不能體現(xiàn)決策者對(duì)各個(gè)指標(biāo)的重視程度并且計(jì)算較復(fù)雜。本文結(jié)合層次分析法和熵權(quán)法并基于加法合成法的原理來確定各個(gè)指標(biāo)的權(quán)重。公式如下［23］：

式中：wi(i= 1，2，…，n)表示第i評(píng)價(jià)指標(biāo)的綜合權(quán)重，αi表示第i評(píng)價(jià)指標(biāo)通過層次分析法得到的權(quán)重值，βi表示第i評(píng)價(jià)指標(biāo)通過熵權(quán)法得到的權(quán)重值，k1表示綜合權(quán)重中層次分析法所占的比重，k2表示綜合權(quán)重中熵權(quán)法所占的比重。

利用上文構(gòu)建評(píng)價(jià)指標(biāo)體系中的7 個(gè)指標(biāo)并結(jié)合GTD 中2011—2020 年的恐怖襲擊事件數(shù)據(jù)，計(jì)算出各個(gè)指標(biāo)的綜合權(quán)重值，如表4所示。

表4 指標(biāo)權(quán)重得分表

（四）威脅度評(píng)價(jià)

根據(jù)上述模型計(jì)算巴基斯坦境內(nèi)恐怖襲擊事件的綜合得分，其得分區(qū)間在［1.073 2， 3.758 6］，綜合得分越高則表明該恐怖襲擊事件的威脅程度越大。

將7 個(gè)二級(jí)指標(biāo)作為K-Means 聚類方法的輸入項(xiàng)，定義K為［3，8］，分別計(jì)算其輪廓系數(shù)，輪廓系數(shù)的取值范圍在［0，1］，系數(shù)值越高，聚類效果越佳。輪廓系數(shù)公式如下［24］：

式中：xi為樣本點(diǎn)，a(xi)為樣本點(diǎn)xi在同類別內(nèi)到其他點(diǎn)的平均距離，b(xi)為樣本點(diǎn)xi到最近不同類別樣本的平均距離。

根據(jù)K-Means 聚類方法，分別計(jì)算不同聚類數(shù)的輪廓系數(shù)（如圖3 所示），結(jié)果表明，當(dāng)K為5 時(shí)，其對(duì)應(yīng)的輪廓系數(shù)值最大。故將2011—2020 年巴基斯坦境內(nèi)恐怖襲擊事件按威脅程度分為五級(jí)，并利用可視化方法將巴基斯坦境內(nèi)恐怖襲擊事件各個(gè)等級(jí)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)。其中：一級(jí)事件的綜合得分為［1.073 2，1.656 2］，共計(jì)1 985 起；二級(jí)事件的綜合得分為［1.656 3，1.905 3］，共計(jì)4 155 起；三級(jí)事件的綜合得分為［1.912 2，2.159 1］，共計(jì)3 061 起；四級(jí)事件的綜合得分為［2.159 2，2.554 0］，共計(jì)1 201起；五級(jí)事件的綜合得分為［2.554 1，3.758 6］。

圖3 不同聚類數(shù)對(duì)應(yīng)輪廓系數(shù)圖

由圖4 可知，2011—2020 年巴基斯坦境內(nèi)恐怖襲擊事件數(shù)量呈現(xiàn)先升后降的趨勢(shì)，其中在2013 年達(dá)到最高值。每年發(fā)生的恐怖襲擊事件中，不同等級(jí)事件的數(shù)量呈現(xiàn)與恐怖襲擊事件總數(shù)大致相同的趨勢(shì)，二級(jí)事件數(shù)量最多，五級(jí)事件分?jǐn)?shù)最高，事件數(shù)相對(duì)應(yīng)的較少。

圖4 2011—2020年巴基斯坦境內(nèi)恐怖襲擊事件數(shù)量變化圖

從空間分布上來看（如圖5 所示），各等級(jí)恐怖襲擊事件在空間上趨同，即在白沙瓦及其附近區(qū)域最為集中，在其他省份的恐怖襲擊事件分布相對(duì)較為分散。從時(shí)間上來看，2011—2014年，巴基斯坦境內(nèi)恐怖襲擊呈上升態(tài)勢(shì)，其中2011—2012 年呈爆發(fā)式增長(zhǎng)，增長(zhǎng)率達(dá)到61.45%。2014年恐怖襲擊有所下降，但總體事件數(shù)仍較高，同比降幅為2.63%。2014—2020 年，巴基斯坦恐怖襲擊影響區(qū)域逐年下降，其中2014—2015年間同比降幅高達(dá)42.18%。此后5年，恐怖襲擊影響區(qū)域下降趨勢(shì)較平緩。

圖5 2011—2020年巴基斯坦境內(nèi)各個(gè)等級(jí)恐怖襲擊事件分布圖

四、巴基斯坦境內(nèi)恐怖襲擊事件空間特征分析

（一）恐怖襲擊事件空間聚集特征分析

利用核密度分析方式探討巴基斯坦境內(nèi)恐怖襲擊事件在各個(gè)威脅度等級(jí)的聚集地區(qū)（如圖6 所示）。全部事件以及各等級(jí)事件的核心活躍聚集地區(qū)大致相同，主要形成三大核心活躍區(qū)以及不同程度的次級(jí)活躍區(qū)。次級(jí)活躍區(qū)主要是以白沙瓦為中心的核心活躍區(qū)蔓延至周邊產(chǎn)生的許多次級(jí)活躍區(qū)。向巴基斯坦東北方向蔓延至開普省的馬拉坎，向東蔓延至巴基斯坦首都伊斯蘭堡，向西南蔓延至北瓦濟(jì)里斯坦。以奎達(dá)為中心的核心活躍區(qū)向東南蔓延產(chǎn)生2 個(gè)次級(jí)活躍區(qū)，分別是馬希拉巴德和錫比地區(qū)。將單元區(qū)內(nèi)恐怖襲擊事件的威脅度作為權(quán)重進(jìn)行三維核密度分析（如圖7 所示），對(duì)比圖6（a）中全部事件二維核密度圖，發(fā)現(xiàn)2011—2020 年巴基斯坦境內(nèi)恐怖襲擊的三個(gè)核心活躍區(qū)中威脅程度最高的地區(qū)是卡拉奇，其次是開普省的白沙瓦及其附近區(qū)域，再者是俾路支省的奎達(dá)及其附近區(qū)域。

圖7 2011—2020年巴基斯坦境內(nèi)恐怖襲擊事件威脅度評(píng)估圖

（二）恐怖襲擊事件空間熱點(diǎn)分析

1. 空間自相關(guān)分析

依據(jù)巴基斯坦境內(nèi)恐怖襲擊各等級(jí)事件的威脅度進(jìn)行統(tǒng)計(jì)，利用全局性Moran’s I指數(shù)和Gi*指數(shù)法對(duì)恐怖襲擊事件相關(guān)數(shù)據(jù)進(jìn)行分析。

先分別對(duì)巴基斯坦境內(nèi)恐怖襲擊各等級(jí)事件進(jìn)行全局性空間自相關(guān)檢驗(yàn)。Z得分表示標(biāo)準(zhǔn)差的倍數(shù)，標(biāo)準(zhǔn)差反映的是數(shù)據(jù)的離散程度。一般情況下，當(dāng)Z得分大于2.56 時(shí)，表明數(shù)據(jù)呈現(xiàn)出明顯的聚類現(xiàn)象。P值反映的是顯著性檢驗(yàn)的結(jié)果，P值越小則說明顯著性越高。一級(jí)事件至五級(jí)事件的全局性空間自相關(guān)檢驗(yàn)結(jié)果如表5 所示。通過Moran’s I指數(shù)的結(jié)果可知一級(jí)事件至五級(jí)事件的Z值均大于2.56，表明其檢驗(yàn)具有顯著性，巴基斯坦境內(nèi)恐怖襲擊事件各個(gè)等級(jí)的威脅度在空間分布上存在相關(guān)性并且呈聚集狀態(tài)。

表5 全局性Moran’s I指數(shù)

2. 熱點(diǎn)分析

利用Gi*指數(shù)法對(duì)巴基斯坦境內(nèi)全部和各等級(jí)恐怖襲擊事件的局部情況進(jìn)行分析（如圖8 所示），探討巴基斯坦各個(gè)地區(qū)恐怖襲擊事件的頻數(shù)是形成高值的簇或是形成低值的簇，從而可以進(jìn)一步在地圖上刻畫出空間聚集較為明顯的位置和區(qū)域相關(guān)性較強(qiáng)的程度。（1）全部和各等級(jí)恐怖襲擊事件的熱點(diǎn)地區(qū)主要是白沙瓦及其附近區(qū)域、奎達(dá)及其附近區(qū)域、卡拉奇地區(qū)，并且這些熱點(diǎn)地區(qū)較為穩(wěn)定。此外，全部和各等級(jí)恐怖襲擊事件中均無冷點(diǎn)地區(qū)。（2）從圖8 可以看出，在低等級(jí)恐怖襲擊事件中錫比為熱點(diǎn)地區(qū)，但隨著等級(jí)的升高其熱點(diǎn)范圍在不斷縮小。（3）威脅度最為嚴(yán)重的五級(jí)事件中，其熱點(diǎn)地區(qū)范圍相較于其他四個(gè)等級(jí)最小，主要分布在3 個(gè)核心活躍區(qū)。由此可知，巴基斯坦境內(nèi)各個(gè)級(jí)別的恐怖襲擊事件在空間分布上呈現(xiàn)著動(dòng)態(tài)平衡。

圖8 2011—2020年巴基斯坦境內(nèi)恐怖襲擊事件熱點(diǎn)分析圖

3. 新興時(shí)空熱點(diǎn)分析

通過Gi*指數(shù)法分析NetCDF 輸入時(shí)空立方體的變量，利用Mann-Kendall 趨勢(shì)測(cè)試評(píng)估恐怖襲擊事件的熱點(diǎn)和冷點(diǎn)趨勢(shì)。利用新興時(shí)空熱點(diǎn)分析工具，構(gòu)建巴基斯坦境內(nèi)恐怖襲擊事件時(shí)空立方體（如圖9 所示）。由圖9（a）可知，2011—2020 年巴基斯坦境內(nèi)恐怖襲擊的熱點(diǎn)數(shù)較為集中，冷點(diǎn)數(shù)相對(duì)分散?？植酪u擊的熱點(diǎn)區(qū)域主要集中在開普省的白沙瓦及其附近區(qū)域、奎達(dá)及其附近區(qū)域以及卡拉奇。少數(shù)加強(qiáng)的熱點(diǎn)分布于莫赫曼德地區(qū)。持續(xù)熱點(diǎn)分布于白沙瓦附近區(qū)域和卡拉奇。少數(shù)連續(xù)熱點(diǎn)分布于莫赫曼德、巴焦?fàn)柡涂_(dá)及其附近地區(qū)。連續(xù)冷點(diǎn)主要出現(xiàn)在卡拉特、海得拉巴、拉爾卡納、蘇庫(kù)爾以及旁遮普省的大部分地區(qū)。這些地區(qū)已經(jīng)維持了較長(zhǎng)時(shí)間的恐怖襲擊低發(fā)率，同時(shí)，在這些地區(qū)出現(xiàn)聚集現(xiàn)象的強(qiáng)度也沒有隨時(shí)間的推移發(fā)生變化。

圖9 2011—2020年巴基斯坦境內(nèi)恐怖襲擊事件新興時(shí)空熱點(diǎn)分析圖

由圖9（b）可以看出，巴基斯坦境內(nèi)恐怖襲擊一級(jí)事件的連續(xù)熱點(diǎn)主要分布于白沙瓦、奎達(dá)、錫比、卡拉奇和聯(lián)邦直轄部落地區(qū)，少數(shù)分散熱點(diǎn)出現(xiàn)在聯(lián)邦直轄部落地區(qū)。由圖9（c）可以看出，二級(jí)事件大量持續(xù)熱點(diǎn)聚集在白沙瓦及其附近區(qū)域，還有一些持續(xù)熱點(diǎn)分散于海得拉巴地區(qū)。少量持續(xù)熱點(diǎn)分布于奎達(dá)及其附近地區(qū)，持續(xù)冷點(diǎn)地區(qū)主要分布于蘇庫(kù)爾、木爾坦、馬拉坎和旁遮普省的大部分地區(qū)。由圖9（d）可以看出，持續(xù)熱點(diǎn)地區(qū)與二級(jí)事件中大致相同，但其冷點(diǎn)地區(qū)主要分布于巴基斯坦的中部地區(qū)。由圖9（e）可以看出，熱點(diǎn)地區(qū)減少，有少量持續(xù)熱點(diǎn)地區(qū)聚集在卡拉奇和白沙瓦附近地區(qū)，持續(xù)冷點(diǎn)地區(qū)主要分布于俾路支省和伊斯蘭堡附近地區(qū)。由圖9（f）可以看出，五級(jí)事件中已經(jīng)沒有熱點(diǎn)地區(qū)，振蕩的冷點(diǎn)地區(qū)分布于白沙瓦及其附近區(qū)域，一些持續(xù)冷點(diǎn)地區(qū)數(shù)量少但較為分散，在巴基斯坦各個(gè)省份都有涉及。

五、結(jié)論與討論

（一）結(jié)論

本文對(duì)巴基斯坦境內(nèi)2011—2020 年恐怖襲擊事件進(jìn)行威脅度評(píng)價(jià)，并運(yùn)用核密度、空間自相關(guān)和新興時(shí)空熱點(diǎn)等分析方法，對(duì)巴基斯坦境內(nèi)恐怖襲擊事件的空間特征進(jìn)行分析，研究其恐怖襲擊事件的空間分布特征。

一是巴基斯坦境內(nèi)恐怖襲擊事件的危害程度從低到高分為五個(gè)等級(jí)，級(jí)別越高危害程度越大。其中2011—2020 年恐怖襲擊事件中，二級(jí)事件的數(shù)量最多，達(dá)到4 155 起，五級(jí)事件相對(duì)最少，達(dá)到228起。在巴基斯坦境內(nèi)恐怖襲擊事件中，恐怖分子主要選擇平民、軍隊(duì)和警察作為攻擊目標(biāo)。這種選擇涉及政治動(dòng)機(jī)、軍事戰(zhàn)略、維持社會(huì)秩序和意識(shí)形態(tài)沖突等多個(gè)因素。爆炸和武裝襲擊成為巴基斯坦境內(nèi)恐怖分子的主要攻擊方式。由于巴基斯坦內(nèi)部安全形勢(shì)復(fù)雜，社會(huì)面武器存在量較多，恐怖分子可以利用這種環(huán)境獲取武器和實(shí)施襲擊。此外，宗教和政治動(dòng)機(jī)也推動(dòng)了恐怖組織選擇攻擊巴基斯坦政府和與其價(jià)值觀和體制相悖的目標(biāo)。綜合這些因素，爆炸和武裝沖突成為巴基斯坦境內(nèi)恐怖分子的主要攻擊手段。

二是巴基斯坦境內(nèi)恐怖襲擊聚集區(qū)域大體上可以分為三個(gè)核心活躍區(qū)，分別是白沙瓦、奎達(dá)、卡拉奇。白沙瓦及其附近地區(qū)是巴基斯坦境內(nèi)恐怖分子活動(dòng)最為頻繁的地方之一，位于與阿富汗接壤的地區(qū)，恐怖分子利用邊境的薄弱點(diǎn)在兩國(guó)之間進(jìn)行跨境襲擊。此外，白沙瓦周邊地區(qū)地形復(fù)雜，包括山區(qū)和丘陵地帶，為恐怖分子提供了藏匿的便利條件。而奎達(dá)作為俾路支省的省會(huì)，也是恐怖分子活動(dòng)的熱點(diǎn)地區(qū)之一。俾路支省位于巴基斯坦西南部，該省存在著政治、社會(huì)和經(jīng)濟(jì)的不穩(wěn)定因素，為恐怖組織提供了滲透和擴(kuò)張的機(jī)會(huì)。奎達(dá)作為該省的經(jīng)濟(jì)中心和交通樞紐，吸引了大量的流動(dòng)人口和商業(yè)活動(dòng)，同時(shí)也成為恐怖分子的襲擊目標(biāo)?？ɡ孀鳛樾诺率〉氖?huì)和巴基斯坦的最大城市，是巴基斯坦的商業(yè)、金融和工業(yè)中心，這使其成為恐怖組織實(shí)施經(jīng)濟(jì)和政治攻擊的首選目標(biāo)，同時(shí)，卡拉奇自身的治安問題，包括犯罪、極端主義和分裂主義等，也為恐怖分子發(fā)動(dòng)襲擊提供了空間。

三是在對(duì)恐怖襲擊事件進(jìn)行全局自相關(guān)分析時(shí)，各等級(jí)事件的Z值均大于2.56，表明其檢驗(yàn)具有顯著性，巴基斯坦境內(nèi)恐怖襲擊事件的發(fā)生在空間分布上存在相關(guān)性并且呈聚集狀態(tài)。巴基斯坦境內(nèi)恐怖主義活動(dòng)在一定程度上與恐怖組織的策略和行動(dòng)有關(guān)，恐怖組織往往傾向于選擇邊境地區(qū)、交通要道、人口密集區(qū)等目標(biāo)區(qū)域進(jìn)行襲擊，以追求最大的影響力和威脅效果。因此，當(dāng)這些區(qū)域發(fā)生恐怖襲擊事件時(shí)，周邊地區(qū)也會(huì)受到波及，導(dǎo)致事件數(shù)量在空間上呈現(xiàn)聚集狀態(tài)。巴基斯坦境內(nèi)恐怖主義活動(dòng)的空間分布也受到政治和社會(huì)環(huán)境的影響。此外，巴基斯坦的地理環(huán)境對(duì)恐怖主義的聚集現(xiàn)象也具有一定影響。例如，接壤阿富汗的邊境地區(qū)，由于地形復(fù)雜和邊境管控困難，恐怖分子更容易利用邊境進(jìn)行跨境襲擊，導(dǎo)致該地區(qū)和周邊地區(qū)恐怖襲擊事件數(shù)量呈現(xiàn)出聚集狀態(tài)。

四是通過Gi*指數(shù)法對(duì)巴基斯坦境內(nèi)各級(jí)別恐怖襲擊事件頻次熱點(diǎn)分析，對(duì)比發(fā)現(xiàn)：巴基斯坦的高熱點(diǎn)地區(qū)基本不變，只有在北瓦濟(jì)里斯坦其熱點(diǎn)地區(qū)不斷縮小。通過新興時(shí)空熱點(diǎn)分析發(fā)現(xiàn)，巴基斯坦境內(nèi)恐怖襲擊事件的新興熱點(diǎn)模式主要以連續(xù)熱點(diǎn)和持續(xù)熱點(diǎn)為主且相對(duì)穩(wěn)定。冷點(diǎn)區(qū)域大部分以連續(xù)冷點(diǎn)的形式出現(xiàn)，新增冷點(diǎn)則主要出現(xiàn)在馬拉坎附近區(qū)域。巴基斯坦高熱點(diǎn)地區(qū)變化較小，說明這些地區(qū)長(zhǎng)期以來一直是恐怖襲擊事件的高發(fā)區(qū)。這與這些地區(qū)的政治和社會(huì)環(huán)境有關(guān)，比如與恐怖組織有關(guān)的權(quán)力斗爭(zhēng)、民族沖突、薄弱的治安控制等導(dǎo)致這些地區(qū)一直處于高風(fēng)險(xiǎn)狀態(tài)。北瓦濟(jì)里斯坦地區(qū)的熱點(diǎn)地區(qū)在不斷縮小，主要原因是政府在該區(qū)域采取了更加積極的反恐措施。例如，政府加強(qiáng)了在這一地區(qū)的治安力量部署，加大了對(duì)恐怖分子的打擊力度，以及推行了一些提高當(dāng)?shù)鼐用裆钏胶蜕鐣?huì)穩(wěn)定的措施，使恐怖襲擊事件的頻次和熱點(diǎn)范圍減少。此外，連續(xù)熱點(diǎn)和持續(xù)熱點(diǎn)是巴基斯坦境內(nèi)恐怖襲擊事件的主要模式，表明在某些地區(qū)，恐怖襲擊事件會(huì)連續(xù)在一段時(shí)間內(nèi)持續(xù)發(fā)生，這與地區(qū)治安控制的困難、恐怖組織的庇護(hù)與藏身地的存在有關(guān)。新增冷點(diǎn)主要出現(xiàn)在馬拉坎附近區(qū)域，這是由于政府在該地區(qū)加大了反恐力度，使恐怖組織在這一地區(qū)的庇護(hù)或活動(dòng)受到打擊。

（二）討論

目前，“中巴經(jīng)濟(jì)走廊”在建和擬建的項(xiàng)目主要集中在巴基斯坦的俾路支省、旁遮普省和信德省。從巴基斯坦國(guó)內(nèi)各地區(qū)項(xiàng)目分布和恐怖襲擊威脅程度的角度來看，這些地區(qū)都屬于高風(fēng)險(xiǎn)地帶。巴基斯坦國(guó)內(nèi)的恐怖組織錯(cuò)綜復(fù)雜，尤其是俾路支省當(dāng)?shù)氐目植澜M織對(duì)“中巴經(jīng)濟(jì)走廊”建設(shè)的威脅最大，近年來絕大多數(shù)針對(duì)中國(guó)的恐怖襲擊事件主要發(fā)生在俾路支?。?5］。

我國(guó)應(yīng)積極采取相應(yīng)的政策措施來應(yīng)對(duì)這一挑戰(zhàn)。第一，要加強(qiáng)安全合作與情報(bào)交流。中巴兩國(guó)應(yīng)加強(qiáng)情報(bào)共享，加強(qiáng)邊境管理，共同打擊恐怖組織。同時(shí)，與國(guó)際社會(huì)加強(qiáng)聯(lián)合反恐行動(dòng)，分享情報(bào)和經(jīng)驗(yàn)，共同應(yīng)對(duì)跨國(guó)恐怖威脅。第二，要加大基礎(chǔ)設(shè)施安全投入，建立風(fēng)險(xiǎn)管控機(jī)制。兩國(guó)應(yīng)加強(qiáng)對(duì)“中巴經(jīng)濟(jì)走廊”建設(shè)沿線地區(qū)的基礎(chǔ)設(shè)施安全投入，包括加大巡邏力度、建立監(jiān)控系統(tǒng)、加強(qiáng)邊境管控以及加強(qiáng)恐怖主義防范和處理能力的培訓(xùn)。第三，要提升經(jīng)濟(jì)合作與社會(huì)發(fā)展水平。通過加大對(duì)俾路支省等地區(qū)的經(jīng)濟(jì)援助和發(fā)展投入，改善當(dāng)?shù)孛癖姷纳钏胶蜕鐣?huì)穩(wěn)定，減少恐怖組織滲透的空間。同時(shí)，加強(qiáng)職業(yè)培訓(xùn)和提供就業(yè)機(jī)會(huì)，降低恐怖組織招募成員的可能性。第四，要增強(qiáng)民眾反恐意識(shí)和宣傳教育。通過加大對(duì)民眾的安全意識(shí)和反恐意識(shí)的宣傳教育，提高民眾對(duì)恐怖主義的認(rèn)識(shí)和對(duì)“中巴經(jīng)濟(jì)走廊”建設(shè)的支持。