(12)发明专利申请
(10)申请公布号 CN 111950155 A(43)申请公布日 2020.11.17
(21)申请号 202010819720.6(22)申请日 2020.08.14
(71)申请人 江苏师范大学
地址 221000 江苏省徐州市铜山区上海路
101号(72)发明人 陈辰 王淑良 顾喜丰 (51)Int.Cl.
G06F 30/20(2020.01)G06F 16/951(2019.01)G06F 16/9537(2019.01)G06Q 10/06(2012.01)G06Q 50/26(2012.01)
权利要求书2页 说明书5页 附图1页
(54)发明名称
一种城市公共交通系统脆弱性评估方法(57)摘要
一种城市公共交通系统的脆弱性评估方法，包括以下步骤：1)首先利用公交网站Web端JavaScript API接口获取公交站点和线路信息，并对数据进行处理。2)利用Space L方法对公共交通系统建模，并将站点、线路信息作为属性赋予给节点。3)根据灾害事件的不同特点，对灾害事件进行建模。4)提出一种新的公共交通网络脆弱性指标，表示网络性能的变化。5)模拟不同灾害事件下对公共交通网络脆弱性进行仿真验证，观察记录仿真过程中网络模型的脆弱性指标变化，得出脆弱性评估。本发明的一种城市公共交通系统的脆弱性评估方法对公共交通系统的日常维护和管理、灾后资源分配提供支持。
CN 111950155 ACN 111950155 A
权 利 要 求 书
1/2页
1.一种城市公共交通网络脆弱性评估方法,其特征在于，包括以下步骤：S1：首先利用公交网站Web端JavaScript API接口获取数据进行处理，获取城市公共交通站点线路信息。
S2：利用space L的方法对公共交通系统建模，并将站点线路信息作为属性赋予给节点。
S3：根据灾害事件的不同特点，对灾害事件进行分析建模。S4：提出一种新的公共交通网络脆弱性指标，表示网络灾后的性能变化。S5：模拟不同灾害事件下公共交通网络脆弱性，进行仿真验证，观察记录仿真过程中网络模型的脆弱性指标变化，得出脆弱性评估。
2.根据权力要求1所述的一种城市公共交通网络脆弱性评估方法，其特征在于，所述步骤S1的用公交网站Web端JavaScript API接口获取公共交通站点线路信息具体为：
公共交通信息包括站点信息和线路信息，本文利用Python语言编写数据爬取算法从公交网上得到公共交通信息；站点经纬度信息获取则是基于公共交通信息，结合地图软件API；利用Web服务API中的Place API得到公共交通站点的经纬度信息；主要信息获取流程为：首先读取已获取的站点信息数据；调用Place API，采用顺序遍历法，遍历站点名称得到对应站点的经纬度信息；接着修正数据，对由于站点名称差异性存在的遗漏站点信息进行处理，并对自环及独立线路进行删除；最后存储数据，对站点及对应线路信息、经纬度信息进行存储。
3.根据权力要求1所述的一种城市公共交通网络脆弱性评估方法，其特征在于，所述步骤S2中Space L的方法对于公共交通网络建模规则具体为：
Space L建模方法是以自然网络形态呈现出的网络结构，Space L方法将公共交通网络中的站点抽象为节点，将相邻两站点之间的区间；抽象为边，当两个车站间存在直接相连的轨道线路时，那么在拓扑结构中这两个节点间具有一条连边，反应两者之间的直接相连关系。采用Space L方法构建的城市轨道交通网络拓扑结构，反映了网络的自然连接状态，这种自然结构便于分析网络最基本的特征及其形态；
轨道交通网络用图G＝{V,E}来表示，V＝{vi|i＝1,2,L,N}表示节点的集合，E＝{lij|lij＝(vi,vj)}表示为边的集合，N是节点的数量，当节点i和节点j直接相连的时候lij＝1，相反当节点i和j不直接相连lij＝0。
4.根据权力要求1所述的一种城市公共交通网络脆弱性评估方法，其特征在于，所述步骤S4中的网络脆弱性指标具体为：
假设网络中节点i的度为ki，则其重要度可以定义为：
对于ki＝0的节点不做考虑，可以定义网络结构熵为：
为消除节点数目N对SE的影响，对网络结构熵进行归一化，定义为网络的标准结构熵。
2
CN 111950155 A
权 利 要 求 书
2/2页
5.根据权力要求1所述的一种城市公共交通系统脆弱性评估方法，其特征在于，所述步骤S5中的脆弱性模拟仿真具体为：
从四种不同的灾害情形来分析城市公共交通网络的脆弱性，分别为对最大节点度、最高节点介数、最大局部中心值和最高k-shell值的站点蓄意攻击来模拟城市公共交通遭受攻击的情景，来评估和对比在不同灾害情形下网络标准结构熵的变化，找出网络中的关键节点，从而进行重点保护。
3
CN 111950155 A
说 明 书
一种城市公共交通系统脆弱性评估方法
1/5页
技术领域
[0001]本发明涉及公共交通系统的脆弱性和应急管理领域，具体涉及一种城市公共交通系统脆弱性评估方法。
背景技术
[0002]城市化进程的加快，使得我国的城市公共交通系统逐渐趋向规模化和网络化。城市公共交通相比于其他的交通方式存在明显的优势，运输能力大并且具有确定的线路等特点使其很大程度上优化了城市空间布局，缓解了城市的交通压力。城市公共交通系统是居民出行和城市职能发挥的基础工程设施，正常有序的交通网络对国家或地区的社会经济发展有着至关重要的作用。因此，应重视对公共交通系统方面的研究，降低事故发生的强度和频率。
发明内容
[0003]本发明的内容是提供一种城市公共交通系统脆弱性评估方法，用于评估网络中的重要的站点，对运营商和后期的公共交通建设维护提供理论支持。[0004]为了实现上述目的，本发明采取了如下技术方案；[0005]S1：首先利用公交网站Web端JavaScript API接口获取数据进行处理，获取城市公共交通站点线路信息。[0006]S2：根据space L的方法对公共交通系统建模，并将站点线路信息作为属性赋予给节点。
[0007]S3：根据灾害事件的不同特点，对灾害事件进行分析建模。[0008]S4：提出一种新的公共交通网络脆弱性指标，表示灾后网络的性能变化。[0009]S5：模拟不同灾害事件下公共交通网络脆弱性，进行仿真验证，观察记录仿真过程中网络模型的脆弱性指标变化，得出脆弱性评估。[0010]进一步地，所述步骤S1的用公交网站Web端JavaScript API接口获取公共交通站点线路信息具体为：
[0011]公共交通信息包括站点信息和线路信息，本文利用Python语言编写数据爬取算法从8684公交网(https://www.8684.cn/)上得到公共交通信息。站点经纬度信息获取则是基于公共交通信息，结合高德地图API。利用Web服务API中的Place API得到公共交通站点的经纬度信息。主要信息获取流程为：
[0012]首先读取以获取的站点信息数据；调用Place API，采用顺序遍历法，遍历站点名称得到对应站点的经纬度信息；接着修正数据，对由于站点名称差异性存在的遗漏站点信息进行处理，并对自环及独立线路进行删除；最后存储数据，对站点及对应线路信息、经纬度信息进行存储。[0013]进一步地，所述步骤S2中的Space L的方法对于公共交通网络建模规则具体为：[0014]Space L建模方法是以自然网络形态呈现出的网络结构，Space L方法将公共交通
4
CN 111950155 A
说 明 书
2/5页
网络中的站点抽象为节点，将相邻两站点之间的区间(路段)抽象为边，当两个车站间存在直接相连的轨道线路时，那么在拓扑结构中这两个节点间具有一条连边，反应两者之间的直接相连关系。采用Space L方法构建的城市轨道交通网络拓扑结构，反映了网络的自然连接状态，这种自然结构便于分析网络最基本的特征及其形态。[0015]对于公共交通系统的建模，通常我们分别将站点和线路抽象为节点和边，一般来说，轨道交通网络可以用图G＝{V,E}来表示，这里V＝{vi|i＝1,2,L,N}表示节点的集合，E＝{lij|lij＝(vi,vj)}表示为边的集合，N是节点的数量，当节点i和节点j直接相连的时候lij＝1，相反当节点i和j不直接相连lij＝0。[0016]进一步地，所述步骤S3中灾害事件建模具体为：[0017]脆弱性是系统固有属性，不同的扰动策略会暴露出系统不同程度的脆弱性。通过模拟不同的灾害模式，在随机扰动和蓄意扰动下分析网络脆弱性。随机扰动从研究的网络节点中等概率选择节点进行扰动失效。蓄意扰动根据节点的重要性对节点进行排序，按照重要性顺序对节点进行攻击。根据灾害事件的特点，扰动分为以下四类：随机攻击，比如因设备故障和操作失误引起的扰动，这种扰动类型繁杂，并且具有高度的不确定性。[0018]基于局部属性的攻击。这种扰动通常是指系统中关键元件受到破坏，主要分为基于度值的攻击和基于局部中心性的攻击。基于度值的扰动即是根据度值度量节点重要性，根据节点重要性的排序列表进行攻击。节点i的度ki为与节点i连接的边的数目，网络中与该节点相连的边越多，该节点在网络中发挥的作用越显著。基于局部中心性的攻击是根据节点最近邻节点的度值和次近邻节点的度值来评估节点重要性的方式，其计算方式如下：
[0019]
其中Ui为节点i最近邻居集合,Uj为节点j最近邻居集合,N(u)为节点u最近邻居数和次近邻居数之和。
[0021]基于全局属性的攻击(也称基于介数的攻击)。例如地震、洪水等自然灾害引起的大范围元件同时故障。节点i的介数为经过节点i的最短路径数与网络中所有最短路径数的比值，其计算方式如下：
[0022][0023]
[0020]
其中，gst为从节点s到t的最短路径数目，为节点s到节点t的最短路径gst中经过
节点i的最短路径数目。
[0024]基于距离属性的攻击(也称基于k-shell值的攻击)。通过k壳分解法，对节点进行分层，通过初始网络节点的壳值来度量节点重要性，根据节点重要性的排序对网络进行扰动。
[0025]进一步地，所述步骤S4中的网络脆弱性指标具体为：[0026]熵是系统的一种无序的度量。网络结构熵定义的引出可以更简洁地度量复杂网络有序状态，它是定义在度分布上的。假设网络中节点i的度为ki，则其重要度可以定义为：
[0027]
5
CN 111950155 A[0028][0029][0030]
说 明 书
3/5页
对于ki＝0的节点不做考虑，可以定义网络结构熵为：
为消除节点数目N对SE的影响，可以对网络结构熵进行归一化，定义为网络的标准
结构熵。
[0031]
进一步地，所述步骤S5中的脆弱性模拟仿真具体为：
[0033]本发明将从四种不同的灾害情形来分析城市公共交通网络的脆弱性。分别为对最大节点度、最高节点介数、最大局部中心值和最高k-shell值的站点蓄意攻击来模拟城市公共交通遭受攻击的情景，来评估和对比在不同灾害情形下网络标准结构熵的变化，找出网络中的关键节点，从而对其进行重点保护。[0034]本发明的有益效果是：本发明通过获取公共交通站点、线路的地理信息，结合Space L模型建立合乎实际的城市公共交通网络模型。模拟不同场景的灾害模式，对比在不同灾害情形下网络脆弱性指标的变化，并且可以找出对网络中影响较大站点。本发明目的在于提供一种城市公共交通系统脆弱性评估方法，脆弱性研究是为了保证城市公共交通系统在各种灾害作用下，找出其脆弱性部分，加强日常维护，从而保持系统的稳定性，履行公共交通系统的基本运输能力。
附图说明
[0035]图1为本发明一种城市公共交通系统脆弱性评估方法的总理流程。
具体实施方式
[0036]为了便于理解下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅是本发明的一部分实施例。[0037]如图1所示，一种城市公共交通系统脆弱性评估方法，主要包括以下步骤：[0038]S1：首先利用公交网站Web端JavaScript API接口获取数据进行处理，获取城市公共交通站点线路信息。
[0039]公共交通信息包括站点信息和线路信息，本文利用Python语言编写数据爬取算法从8684公交网(https://www.8684.cn/)上得到公共交通信息。站点经纬度信息获取则是基于公共交通信息，结合高德地图API。利用Web服务API中的Place API得到公共交通站点的经纬度信息。主要信息获取流程为：
[0040]首先读取已获取的站点信息数据；调用Place API，采用顺序遍历法，遍历站点名称得到对应站点的经纬度信息；接着修正数据，对由于站点名称差异性存在的遗漏站点信息进行处理，并对自环及独立线路进行删除；最后存储数据，对站点及对应线路信息、经纬度信息进行存储。[0041]S2：利用space L的方法对公共交通系统建模，并将站点线路信息作为属性赋予给节点。
6
[0032]
CN 111950155 A[0042]
说 明 书
4/5页
Space L建模方法是以自然网络形态呈现出的网络结构，Space L方法将公共交通
网络中的站点抽象为节点，将相邻两站点之间的区间(路段)抽象为边，当两个车站间存在直接相连的轨道线路时，那么在拓扑结构中这两个节点间具有一条连边，反应两者之间的直接相连关系。采用Space L方法构建的城市轨道交通网络拓扑结构，反映了网络的自然连接状态，这种自然结构便于分析网络最基本的特征及其形态。[0043]对于公共交通系统的建模，通常我们分别将站点和线路抽象为节点和边，一般来说，轨道交通网络可以用图G＝{V,E}来表示，这里V＝{vi|i＝1,2,L,N}表示节点的集合，E＝{lij|lij＝(vi,vj)}表示为边的集合，N是节点的数量，当节点i和节点j直接相连的时候lij＝1，相反当节点i和j不直接相连lij＝0。[0044]S3：根据灾害事件的不同特点，对灾害事件进行分析建模。[0045]脆弱性是系统固有属性，不同的扰动策略会暴露出系统不同程度的脆弱性。通过模拟不同的灾害模式，在随机扰动和蓄意扰动下分析网络脆弱性。随机扰动从研究的网络节点中等概率选择节点进行扰动失效。蓄意扰动根据节点的重要性对节点进行排序，按照重要性顺序对节点进行攻击。根据灾害事件的特点，扰动分为以下四类：随机攻击，比如因设备故障和操作失误引起的扰动，这种扰动类型繁杂，并且具有高度的不确定性。[0046]基于局部属性的攻击。这种扰动通常是指系统中关键元件受到破坏，主要分为基于度值的攻击和基于局部中心性的攻击。基于度值的扰动即是根据度值度量节点重要性，根据节点重要性的排序列表进行攻击。节点i的度ki为与节点i连接的边的数目，网络中与该节点相连的边越多，该节点在网络中发挥的作用越显著。基于局部中心性的攻击是根据节点最近邻节点的度值和次近邻节点的度值来评估节点重要性的方式，其计算方式如下：
[0047]
其中Ui为节点i最近邻居集合,Uj为节点j最近邻居集合,N(u)为节点u最近邻居数和次近邻居数之和。
[0049]基于全局属性的攻击(也称基于介数的攻击)。例如地震、洪水等自然灾害引起的大范围元件同时故障。节点i的介数为经过节点i的最短路径数与网络中所有最短路径数的比值，其计算方式如下：
[0050][0051]
[0048]
其中，gst为从节点s到t的最短路径数目，为节点s到节点t的最短路径gst中经过
节点i的最短路径数目。
[0052]基于距离属性的攻击(也称基于k-shell值的攻击)。通过k壳分解法，对节点进行分层，通过初始网络节点的壳值来度量节点重要性，根据节点重要性的排序对网络进行扰动。
[0053]S4：提出一种新的公共交通网络脆弱性指标，表示网络灾后的性能变化。[0054]熵是系统的一种无序的度量。网络结构熵定义的引出可以更简洁地度量复杂网络有序状态，它是定义在度分布上的。假设网络中节点i的度为ki，则其重要度可以定义为：
7
CN 111950155 A
说 明 书
5/5页
[0055]
对于ki＝0的节点不做考虑，可以定义网络结构熵为：
[0059]为消除节点数目N对SE的影响，可以对网络结构熵进行归一化，定义为网络的标准结构熵。
[0060]
[0058]
以上结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述
实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。
[0061]
8
CN 111950155 A
说 明 书 附 图
1/1页
图1
9
因篇幅问题不能全部显示，请点此查看更多更全内容