基于威胁度评价的巴基斯坦境内恐怖袭击事件空间分布特征研究

徐顺聪，杨永斌，张世伟，张 磊

中国人民警察大学，河北 廊坊 065000

一、引言

近年来，全球恐怖主义风险仍处于高发态势，恐怖主义问题引发全世界的广泛关注。巴基斯坦是“一带一路”沿线的重要国家，“中巴经济走廊”是“六廊六路”的重要组成部分，我国在巴项目众多，经济及战略意义重大，巴基斯坦境内恐怖袭击事件也同时威胁我国在巴项目安全。2021 年7 月14 日，由中方企业承建的开普省达苏水电站项目出勤班车在赴施工现场途中遭遇袭击爆炸，造成包括9 名中方人员和3名巴方人员在内的12人死亡。2021年8月20日，一名自杀式袭击者在俾路支省瓜达尔地区袭击载有中国公民的车队，造成2 名儿童死亡、3 人受伤。巴基斯坦境内恐怖主义风险较为突出，根据《2023年全球恐怖主义指数》，巴基斯坦的全球恐怖主义指数排名已经上升到第六名，是恐怖主义较为严重和活跃的地区［1］。研究巴基斯坦境内恐怖袭击事件对保护我国在巴项目安全，增强我国海外利益保护具有重要意义。

当前，对于恐怖袭击事件的预防预警工作已有一定的研究成果。hao 等为揭示东南亚恐怖袭击的分布规律提出基于反距离加权和基于多标签K-最近邻的评估模型［2］。师原兵在城市轨道交通领域中应用大数据技术来加强对数据的分析与研判，并提出从多个方面提升反恐防范和应急处置能力，有效打击恐怖主义犯罪并维护国家安全［3］。邓世成等采用投影寻踪评价模型和灰色新陈代谢模型对全球恐怖袭击事件的危害程度进行量化分级，预测了未来5年全球恐怖活动的蔓延趋势［4］。罗澜峻等利用长短时记忆神经网络模型对恐怖袭击事件的发生时间进行预测，结果表明该模型能够较准确预测出短时内恐怖袭击事件的发生时间［5］。这些文献大多研究恐怖袭击事件的动机和成因，很少与地理空间信息相结合，通过可视化的方法分析恐怖袭击事件。

在恐怖袭击时空演变分析方面，国内外已经有许多学者对此进行了研究。王诚聪等采用热点分析和层次聚类分析的方法对全球恐怖袭击造成的伤亡人数和事件频次进行时空演变特征分析［6］。姜春艳等对全球恐怖袭击事件进行探索，根据恐怖袭击事件的时间序列以及空间特征指出恐怖袭击活动主要集中在南亚、中东和北非、撒哈拉以南的非洲及东南亚地区［7］。Hao 等采用地理空间统计方法并结合随机森林算法，对中南半岛恐怖袭击的时空演变规律进行分析［8］。同样，陈伟等将核密度分析用于恐怖袭击聚集度的可视化以及标准差椭圆用于分析恐怖袭击的空间演变特征［9］。郝蒙蒙等利用空间统计方法和核密度分析的方法分析中南半岛恐怖袭击事件的时空演变特征［10］。李继萌等运用统计分析法和核密度估计法指出“一带一路”沿线国家遭受的恐怖袭击事件在空间上呈现出高度聚集特征并且向四周扩散，国家交界处是恐怖袭击事件多发地带［11］。

当前对巴基斯坦境内恐怖袭击的研究除集中于定性研究外，也开始出现利用时空统计技术对其分布特征的研究［12］。但在前期研究中将恐怖袭击事件等同对待，更多的是以恐怖主义袭击的次数作为分析对象，对不同恐怖袭击事件之间危害性差异考虑较少。本文选用2011—2020 年巴基斯坦境内恐怖袭击数据，通过构建衡量不同恐怖袭击事件危害性的评价指数——威胁度，建立威胁度评价指标体系及量化模型（如图1 所示），对恐怖袭击事件进行威胁度分级。随后，采用核密度分析、空间自相关分析以及新兴时空热点分析等方法研究巴基斯坦境内恐怖袭击事件的危害程度及空间演变特征，以期对未来预防恐怖袭击提供重要的参考价值。

图1 威胁度评价指标体系及量化模型

二、研究区域与研究方法

（一）研究区域概况

巴基斯坦，全称巴基斯坦伊斯兰共和国，位于南亚次大陆的西北部，地理范围是北纬24°37′至37°06′，东经61°33′至75°31′。西部与伊朗接壤，南部和东南部与印度接壤，东部与阿富汗毗邻，北部与中国相邻，国土面积（不含巴控克什米尔地区）为79.6 万平方公里。巴基斯坦的行政区划包括：四个省份，分别是旁遮普省、开伯尔-普赫图赫瓦省（以下简称开普省）、俾路支省和信德省；两个行政区域，分别为自由克什米尔和吉尔吉特-巴尔蒂斯坦、伊斯兰堡首都地区。

（二）研究方法

1. 威胁度分级

对2011—2020 年巴基斯坦境内恐怖袭击事件的数据进行综合分析，以地理空间统计学的视角，建立一个完善的指标体系和分析模型，以期获得更准确的结论。参照《生产安全事故报告和调查处理条例》将恐怖袭击事件划分为四个等级，对数据进行离散化处理，将恐怖袭击事件威胁度评价指标按1～4进行赋值，通过K-Means 聚类分析的方法对事件进行分级处理。

2. 核密度分析

核密度估计法（Kernel Density Estimation，KDE）认为地理事件可能发生在空间的任何位置，但不同位置上事件发生的概率不同。相对于样方计数法，KDE更适合使用可视化方法来表示分布模式［13］。研究中使用高斯核密度函数，设X1，X2，…，Xn是从分布密度函数f中抽取的独立同分布样本，则估计f在某点x处的密度函数为：

式中：n为样点个数，k(·)为核函数，h为带宽，(x-Xi)为点x到事件Xi处的距离。

3. 空间自相关

在对空间进行计量时，需考虑空间与空间之间是否存在相关性。空间相关性是指相邻位置上区域事物发展水平相似。若是高值与高值、低值与低值相邻，则称为“正空间自相关”。反之，若是高值与低值相邻，则称为“负空间自相关”。如果高值与低值完全随机分布，则不存在空间自相关。全局Moran’s I指数法计算公式如下［14］：

式中：S0为所有空间权重之和，S0可以表示为S0=为要素i的属性与其平均值(x-)的偏差；n为样本总数；wij为要素i和j之间的空间权重。

全局性Moran’s I 指数法检验的标准化统计量为：

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

4. 热点识别法

可通过Gi*指数法对巴基斯坦境内恐怖袭击局部情况进行分析，Gi*指数可以比较准确地探测出局部空间聚集区域［15］。其公式如下［16］：

式中：xj为要素j的属性值，wij为要素i和j之间的空间权重，n为要素总数，期望值

热点分析各模式如表1所示：

表1 热点分析各模式定义

5.Mann-Kendall趋势分析

Mann-Kendall 趋势分析方法是一种非参数统计检验方法，它可用于对独立的时空条柱序列进行测试。对于由n个样本组成的时间序列，趋势检验统计量如下［17］：

式中：sgn(·)是符号函数，若xj＞xi取1，若xj=xi取0，若xj＜xi取-1。

当S服从正态分布时，其方差公式如下：

标准化统计量Z计算公式如下：

根据Z得分对恐怖袭击发生率的冷热点趋势的显著性进行归类，分为新增热点、连续热点、加强热点等冷热点模式（如表2所示）。

表2 新兴时空热点分析各模式定义（部分）

三、巴基斯坦境内恐怖袭击事件威胁度评价

采用威胁度评价对恐怖袭击事件进行量化分析，有助于准确地预测恐怖袭击事件的发生概率。通过量化评估，还可以确定哪些地区、哪些群体、哪些场所面临更大的威胁，从而有针对性地配置防恐资源，提高反恐效率。

（一）数据来源及处理

1. 数据来源

本文的数据来源为全球恐怖主义数据库（Global Terrorism Database，GTD），GTD 是由马里兰大学开发维护的一个开源数据库，下载地址为https：//www. start. umd. edu/gtd/，记录了1970—2020年全球各地恐怖袭击事件的详细信息［18-19］。其中包括恐怖袭击事件发生的时间、国家、经纬度、袭击组织以及造成的伤亡人数等内容。本文选取巴基斯坦2011—2020年的恐怖袭击数据。

2. 数据预处理

（1）缺失值处理。在排除经纬度坐标缺失和伤亡人数有缺失的数据后，共得到10 639 条有效恐怖袭击数据。已有研究表明，当缺失值小于5%时，使用随机森林填补法效果较优，并且随机森林填补法不仅在软件操作上简便，还对数据结构和分布的要求较为宽松，在反映真实情况方面也较为准确［20］。故选择使用sklearn库中的随机森林算法来对有缺失值的变量进行填充。

（2）无量纲处理。在该研究中使用极值处理法消除各指标单位和数值间的系统误差，将原始数据处理成0～1范围内的数据。公式如下［21］：

式中：xij(i= 1，2…n；j= 1，2…m)表示第i个被评价单元在第j项指标上的观测值是经过无量纲化处理后的标准观测值。

（二）评价指标体系的构建

本研究根据GTD 数据库信息，选取7 个指标（如图2 所示），包括死亡人数、受伤人数、财产损失、目标类型、袭击地点、攻击类型和武器类型。其中，定量指标有死亡人数、受伤人数和财产损失。定性指标有目标类型、袭击地点、攻击类型和武器类型。

图2 巴基斯坦境内恐怖袭击事件威胁度指标体系

（三）指标体系权重的确定

在对恐怖袭击事件威胁度进行计算时，通过等级赋分将定性指标量化，研究中分为4 个等级（如表3所示）［22］。

确定指标体系权重的方法可分为两类：一类是主观赋权方法，它的缺点在于过于依赖专家自身的经验与知识，主观性较强；另一类是客观赋权方法，缺点在于过分依赖样本数据，不能体现决策者对各个指标的重视程度并且计算较复杂。本文结合层次分析法和熵权法并基于加法合成法的原理来确定各个指标的权重。公式如下［23］：

式中：wi(i= 1，2，…，n)表示第i评价指标的综合权重，αi表示第i评价指标通过层次分析法得到的权重值，βi表示第i评价指标通过熵权法得到的权重值，k1表示综合权重中层次分析法所占的比重，k2表示综合权重中熵权法所占的比重。

利用上文构建评价指标体系中的7 个指标并结合GTD 中2011—2020 年的恐怖袭击事件数据，计算出各个指标的综合权重值，如表4所示。

表4 指标权重得分表

（四）威胁度评价

根据上述模型计算巴基斯坦境内恐怖袭击事件的综合得分，其得分区间在［1.073 2， 3.758 6］，综合得分越高则表明该恐怖袭击事件的威胁程度越大。

将7 个二级指标作为K-Means 聚类方法的输入项，定义K为［3，8］，分别计算其轮廓系数，轮廓系数的取值范围在［0，1］，系数值越高，聚类效果越佳。轮廓系数公式如下［24］：

式中：xi为样本点，a(xi)为样本点xi在同类别内到其他点的平均距离，b(xi)为样本点xi到最近不同类别样本的平均距离。

根据K-Means 聚类方法，分别计算不同聚类数的轮廓系数（如图3 所示），结果表明，当K为5 时，其对应的轮廓系数值最大。故将2011—2020 年巴基斯坦境内恐怖袭击事件按威胁程度分为五级，并利用可视化方法将巴基斯坦境内恐怖袭击事件各个等级数据进行统计。其中：一级事件的综合得分为［1.073 2，1.656 2］，共计1 985 起；二级事件的综合得分为［1.656 3，1.905 3］，共计4 155 起；三级事件的综合得分为［1.912 2，2.159 1］，共计3 061 起；四级事件的综合得分为［2.159 2，2.554 0］，共计1 201起；五级事件的综合得分为［2.554 1，3.758 6］。

图3 不同聚类数对应轮廓系数图

由图4 可知，2011—2020 年巴基斯坦境内恐怖袭击事件数量呈现先升后降的趋势，其中在2013 年达到最高值。每年发生的恐怖袭击事件中，不同等级事件的数量呈现与恐怖袭击事件总数大致相同的趋势，二级事件数量最多，五级事件分数最高，事件数相对应的较少。

图4 2011—2020年巴基斯坦境内恐怖袭击事件数量变化图

从空间分布上来看（如图5 所示），各等级恐怖袭击事件在空间上趋同，即在白沙瓦及其附近区域最为集中，在其他省份的恐怖袭击事件分布相对较为分散。从时间上来看，2011—2014年，巴基斯坦境内恐怖袭击呈上升态势，其中2011—2012 年呈爆发式增长，增长率达到61.45%。2014年恐怖袭击有所下降，但总体事件数仍较高，同比降幅为2.63%。2014—2020 年，巴基斯坦恐怖袭击影响区域逐年下降，其中2014—2015年间同比降幅高达42.18%。此后5年，恐怖袭击影响区域下降趋势较平缓。

图5 2011—2020年巴基斯坦境内各个等级恐怖袭击事件分布图

四、巴基斯坦境内恐怖袭击事件空间特征分析

（一）恐怖袭击事件空间聚集特征分析

利用核密度分析方式探讨巴基斯坦境内恐怖袭击事件在各个威胁度等级的聚集地区（如图6 所示）。全部事件以及各等级事件的核心活跃聚集地区大致相同，主要形成三大核心活跃区以及不同程度的次级活跃区。次级活跃区主要是以白沙瓦为中心的核心活跃区蔓延至周边产生的许多次级活跃区。向巴基斯坦东北方向蔓延至开普省的马拉坎，向东蔓延至巴基斯坦首都伊斯兰堡，向西南蔓延至北瓦济里斯坦。以奎达为中心的核心活跃区向东南蔓延产生2 个次级活跃区，分别是马希拉巴德和锡比地区。将单元区内恐怖袭击事件的威胁度作为权重进行三维核密度分析（如图7 所示），对比图6（a）中全部事件二维核密度图，发现2011—2020 年巴基斯坦境内恐怖袭击的三个核心活跃区中威胁程度最高的地区是卡拉奇，其次是开普省的白沙瓦及其附近区域，再者是俾路支省的奎达及其附近区域。

图7 2011—2020年巴基斯坦境内恐怖袭击事件威胁度评估图

（二）恐怖袭击事件空间热点分析

1. 空间自相关分析

依据巴基斯坦境内恐怖袭击各等级事件的威胁度进行统计，利用全局性Moran’s I指数和Gi*指数法对恐怖袭击事件相关数据进行分析。

先分别对巴基斯坦境内恐怖袭击各等级事件进行全局性空间自相关检验。Z得分表示标准差的倍数，标准差反映的是数据的离散程度。一般情况下，当Z得分大于2.56 时，表明数据呈现出明显的聚类现象。P值反映的是显著性检验的结果，P值越小则说明显著性越高。一级事件至五级事件的全局性空间自相关检验结果如表5 所示。通过Moran’s I指数的结果可知一级事件至五级事件的Z值均大于2.56，表明其检验具有显著性，巴基斯坦境内恐怖袭击事件各个等级的威胁度在空间分布上存在相关性并且呈聚集状态。

表5 全局性Moran’s I指数

2. 热点分析

利用Gi*指数法对巴基斯坦境内全部和各等级恐怖袭击事件的局部情况进行分析（如图8 所示），探讨巴基斯坦各个地区恐怖袭击事件的频数是形成高值的簇或是形成低值的簇，从而可以进一步在地图上刻画出空间聚集较为明显的位置和区域相关性较强的程度。（1）全部和各等级恐怖袭击事件的热点地区主要是白沙瓦及其附近区域、奎达及其附近区域、卡拉奇地区，并且这些热点地区较为稳定。此外，全部和各等级恐怖袭击事件中均无冷点地区。（2）从图8 可以看出，在低等级恐怖袭击事件中锡比为热点地区，但随着等级的升高其热点范围在不断缩小。（3）威胁度最为严重的五级事件中，其热点地区范围相较于其他四个等级最小，主要分布在3 个核心活跃区。由此可知，巴基斯坦境内各个级别的恐怖袭击事件在空间分布上呈现着动态平衡。

图8 2011—2020年巴基斯坦境内恐怖袭击事件热点分析图

3. 新兴时空热点分析

通过Gi*指数法分析NetCDF 输入时空立方体的变量，利用Mann-Kendall 趋势测试评估恐怖袭击事件的热点和冷点趋势。利用新兴时空热点分析工具，构建巴基斯坦境内恐怖袭击事件时空立方体（如图9 所示）。由图9（a）可知，2011—2020 年巴基斯坦境内恐怖袭击的热点数较为集中，冷点数相对分散。恐怖袭击的热点区域主要集中在开普省的白沙瓦及其附近区域、奎达及其附近区域以及卡拉奇。少数加强的热点分布于莫赫曼德地区。持续热点分布于白沙瓦附近区域和卡拉奇。少数连续热点分布于莫赫曼德、巴焦尔和奎达及其附近地区。连续冷点主要出现在卡拉特、海得拉巴、拉尔卡纳、苏库尔以及旁遮普省的大部分地区。这些地区已经维持了较长时间的恐怖袭击低发率，同时，在这些地区出现聚集现象的强度也没有随时间的推移发生变化。

图9 2011—2020年巴基斯坦境内恐怖袭击事件新兴时空热点分析图

由图9（b）可以看出，巴基斯坦境内恐怖袭击一级事件的连续热点主要分布于白沙瓦、奎达、锡比、卡拉奇和联邦直辖部落地区，少数分散热点出现在联邦直辖部落地区。由图9（c）可以看出，二级事件大量持续热点聚集在白沙瓦及其附近区域，还有一些持续热点分散于海得拉巴地区。少量持续热点分布于奎达及其附近地区，持续冷点地区主要分布于苏库尔、木尔坦、马拉坎和旁遮普省的大部分地区。由图9（d）可以看出，持续热点地区与二级事件中大致相同，但其冷点地区主要分布于巴基斯坦的中部地区。由图9（e）可以看出，热点地区减少，有少量持续热点地区聚集在卡拉奇和白沙瓦附近地区，持续冷点地区主要分布于俾路支省和伊斯兰堡附近地区。由图9（f）可以看出，五级事件中已经没有热点地区，振荡的冷点地区分布于白沙瓦及其附近区域，一些持续冷点地区数量少但较为分散，在巴基斯坦各个省份都有涉及。

五、结论与讨论

（一）结论

本文对巴基斯坦境内2011—2020 年恐怖袭击事件进行威胁度评价，并运用核密度、空间自相关和新兴时空热点等分析方法，对巴基斯坦境内恐怖袭击事件的空间特征进行分析，研究其恐怖袭击事件的空间分布特征。

一是巴基斯坦境内恐怖袭击事件的危害程度从低到高分为五个等级，级别越高危害程度越大。其中2011—2020 年恐怖袭击事件中，二级事件的数量最多，达到4 155 起，五级事件相对最少，达到228起。在巴基斯坦境内恐怖袭击事件中，恐怖分子主要选择平民、军队和警察作为攻击目标。这种选择涉及政治动机、军事战略、维持社会秩序和意识形态冲突等多个因素。爆炸和武装袭击成为巴基斯坦境内恐怖分子的主要攻击方式。由于巴基斯坦内部安全形势复杂，社会面武器存在量较多，恐怖分子可以利用这种环境获取武器和实施袭击。此外，宗教和政治动机也推动了恐怖组织选择攻击巴基斯坦政府和与其价值观和体制相悖的目标。综合这些因素，爆炸和武装冲突成为巴基斯坦境内恐怖分子的主要攻击手段。

二是巴基斯坦境内恐怖袭击聚集区域大体上可以分为三个核心活跃区，分别是白沙瓦、奎达、卡拉奇。白沙瓦及其附近地区是巴基斯坦境内恐怖分子活动最为频繁的地方之一，位于与阿富汗接壤的地区，恐怖分子利用边境的薄弱点在两国之间进行跨境袭击。此外，白沙瓦周边地区地形复杂，包括山区和丘陵地带，为恐怖分子提供了藏匿的便利条件。而奎达作为俾路支省的省会，也是恐怖分子活动的热点地区之一。俾路支省位于巴基斯坦西南部，该省存在着政治、社会和经济的不稳定因素，为恐怖组织提供了渗透和扩张的机会。奎达作为该省的经济中心和交通枢纽，吸引了大量的流动人口和商业活动，同时也成为恐怖分子的袭击目标。卡拉奇作为信德省的省会和巴基斯坦的最大城市，是巴基斯坦的商业、金融和工业中心，这使其成为恐怖组织实施经济和政治攻击的首选目标，同时，卡拉奇自身的治安问题，包括犯罪、极端主义和分裂主义等，也为恐怖分子发动袭击提供了空间。

三是在对恐怖袭击事件进行全局自相关分析时，各等级事件的Z值均大于2.56，表明其检验具有显著性，巴基斯坦境内恐怖袭击事件的发生在空间分布上存在相关性并且呈聚集状态。巴基斯坦境内恐怖主义活动在一定程度上与恐怖组织的策略和行动有关，恐怖组织往往倾向于选择边境地区、交通要道、人口密集区等目标区域进行袭击，以追求最大的影响力和威胁效果。因此，当这些区域发生恐怖袭击事件时，周边地区也会受到波及，导致事件数量在空间上呈现聚集状态。巴基斯坦境内恐怖主义活动的空间分布也受到政治和社会环境的影响。此外，巴基斯坦的地理环境对恐怖主义的聚集现象也具有一定影响。例如，接壤阿富汗的边境地区，由于地形复杂和边境管控困难，恐怖分子更容易利用边境进行跨境袭击，导致该地区和周边地区恐怖袭击事件数量呈现出聚集状态。

四是通过Gi*指数法对巴基斯坦境内各级别恐怖袭击事件频次热点分析，对比发现：巴基斯坦的高热点地区基本不变，只有在北瓦济里斯坦其热点地区不断缩小。通过新兴时空热点分析发现，巴基斯坦境内恐怖袭击事件的新兴热点模式主要以连续热点和持续热点为主且相对稳定。冷点区域大部分以连续冷点的形式出现，新增冷点则主要出现在马拉坎附近区域。巴基斯坦高热点地区变化较小，说明这些地区长期以来一直是恐怖袭击事件的高发区。这与这些地区的政治和社会环境有关，比如与恐怖组织有关的权力斗争、民族冲突、薄弱的治安控制等导致这些地区一直处于高风险状态。北瓦济里斯坦地区的热点地区在不断缩小，主要原因是政府在该区域采取了更加积极的反恐措施。例如，政府加强了在这一地区的治安力量部署，加大了对恐怖分子的打击力度，以及推行了一些提高当地居民生活水平和社会稳定的措施，使恐怖袭击事件的频次和热点范围减少。此外，连续热点和持续热点是巴基斯坦境内恐怖袭击事件的主要模式，表明在某些地区，恐怖袭击事件会连续在一段时间内持续发生，这与地区治安控制的困难、恐怖组织的庇护与藏身地的存在有关。新增冷点主要出现在马拉坎附近区域，这是由于政府在该地区加大了反恐力度，使恐怖组织在这一地区的庇护或活动受到打击。

（二）讨论

目前，“中巴经济走廊”在建和拟建的项目主要集中在巴基斯坦的俾路支省、旁遮普省和信德省。从巴基斯坦国内各地区项目分布和恐怖袭击威胁程度的角度来看，这些地区都属于高风险地带。巴基斯坦国内的恐怖组织错综复杂，尤其是俾路支省当地的恐怖组织对“中巴经济走廊”建设的威胁最大，近年来绝大多数针对中国的恐怖袭击事件主要发生在俾路支省［25］。

我国应积极采取相应的政策措施来应对这一挑战。第一，要加强安全合作与情报交流。中巴两国应加强情报共享，加强边境管理，共同打击恐怖组织。同时，与国际社会加强联合反恐行动，分享情报和经验，共同应对跨国恐怖威胁。第二，要加大基础设施安全投入，建立风险管控机制。两国应加强对“中巴经济走廊”建设沿线地区的基础设施安全投入，包括加大巡逻力度、建立监控系统、加强边境管控以及加强恐怖主义防范和处理能力的培训。第三，要提升经济合作与社会发展水平。通过加大对俾路支省等地区的经济援助和发展投入，改善当地民众的生活水平和社会稳定，减少恐怖组织渗透的空间。同时，加强职业培训和提供就业机会，降低恐怖组织招募成员的可能性。第四，要增强民众反恐意识和宣传教育。通过加大对民众的安全意识和反恐意识的宣传教育，提高民众对恐怖主义的认识和对“中巴经济走廊”建设的支持。