-
In real complex systems, the overall function is maintained through the connections among nodes. Failures of some nodes may destroy the connectivity of the system and thus damage the function of the system. In some complex systems, some nodes can form “interdependency groups” through hidden interdependency. The failure of one node may damage the rest of the nodes in the interdependency group. In this paper, we investigate the effects of the interdependency strength of the nodes, the size distribution, and the size of the interdependency groups on the cascading dynamics and the robustness of complex networks. Through numerical simulation and theoretical analysis, it is found that the cascading failures of the networks can be divided into two processes at a scale level: “intra-group cascading” and “inter-group cascading”. In the intra-group cascading process, the failure of one node will result in damage to the other nodes in the group through the interdependence among nodes, thus inducing more nodes to be unworkable and resulting in greater destructive force. In the inter-group cascading process, the failed nodes will cause the networks to be fragmented, which leads some nodes outside the interdependency group to isolate from the giant component and go to failure. Under the synergistic effects of these two processes, it is found that there are continuous and discontinuous phase transition phenomena in the cascade dynamics of the network. The occurrence of these two kinds of phase transition phenomena is related to the interdependency strength of nodes, the network degree distribution and the size distribution of the interdependency group. This means that by controlling the characteristics of interdependency groups, such as the interdependence strength of the nodes in the interdependency group or the size distribution of interdependency groups, the system can avoid collapsing suddenly and thus the robustness of the network can be improved.
-
Keywords:
- interdependency group /
- weak interdependence /
- cascading failure /
- giant component /
- robustness
1. 引 言
真实世界中的许多复杂系统都可以抽象成复杂网络, 如常见的交通网络[1-4]、电力网络[5]、社交网络[6]、疾病传播网络[7,8]和生物网络[9,10] 等. 在复杂系统鲁棒性的研究中, 人们常常借助于复杂网络级联失效模型对复杂系统的崩溃过程进行建模和研究, 节点的失效不仅会对网络整体的连通性造成影响, 也可能会因节点之间的依赖关系而引起失效在网络中蔓延. 最初复杂网络鲁棒性的研究是将两种效应分开进行的, 即针对网络节点之间的连接关系来探讨网络在失去部分连接或节点之后的连通性[11-13], 以及节点间的耦合或依赖性所引起的失效节点的雪崩或级联效应[14]. 对于前者, 人们借助于渗流模型发现, 随机网络发生渗流相变的临界点与网络平均度呈反比关系[15]; 无标度网络具备“鲁棒性和脆弱性共存”的特征[16-18], 即对随机故障的鲁棒性较强, 但面对蓄意攻击时则显得十分脆弱; 对于小世界网络, “捷径”使得网络的鲁棒性增强[19], 但也更容易在蓄意攻击下破碎[20]. 对于复杂网络级联失效的研究, 如何对节点之间的耦合和依赖进行建模是一个核心问题. 研究人员考虑了真实复杂系统的特征, 主要考虑了动态负载分配、局部依赖和依赖组的引入等几个方案.
1) 动态负载分配[21-27]. 当网络中部分节点失效后, 网络节点之间最短路径的分布发生较大的变化, 网络中信息或能量的流动路径就会重新分配, 从而引起部分节点的运行效率下降[21], 甚至过载失效[27]. 在这类模型中, 一个负载较小节点的失效可能不会引起负载失衡并诱发较大范围的过载故障; 而失效节点的负载较大时, 这些节点的失效可能会显著影响其他节点的负载, 并引发一系列大范围的过载故障[23,27]. 为了探讨级联失效的产生条件并为级联失效的防御和控制提供思路和策略, 节点失效后荷载的分配方案对网络级联失效动力学的影响也受到广泛关注, 如随机荷载再分配机制[28]、基于边权重的再分配策略[29,30]、供需网络中的供应策略[31]等. 除了流量分配机制外, 网络结构调整策略[32]、攻击策略[33]、网络中流量的涨落[34]、网络的结构特征[35]等因素对于级联失效的影响也受到广泛的关注.
2)局部依赖. Watts模型研究了邻居中相反状态节点比例对节点状态的影响, 发现了复杂网络发生级联失效的临界判据[36,37]. 此外, 考虑到网络高度值互连部分节点(k 核)所组成结构的重要性, Dorogovtsev等[38]提出了k核渗流模型, 该模型将通过一个反复迭代的剪枝过程获得的最小度为k的节点所形成的连通子图视为k核. 度分布的二阶矩有限时, 随着初始删除节点的增加, k 核巨分支在某一临界点会突然消失, 表现为不连续的混合相变现象. 如果度分布的第二阶矩发散, 则网络包含一个无限序列的k核, 这些k核对随机攻击具有很强的鲁棒性. 类似地, 靴攀渗流模型考虑了失效邻居的数量对节点状态的影响, 当节点失效邻居数达到某个阈值时该节点就会失效, 研究发现, 在不同的度分布情形下, 网络中也会出现一阶不连续相变和连续相变的现象[39].
3)依赖组的引入. 上述两种模型中节点都是通过相互邻接关系进行状态传递的, 但现实网络中不相邻的节点也可能存在隐含的相互影响. Parshani等[40]考虑到不相邻节点之间也可能存在依赖性, 于是提出了依赖组概念来描述复杂网络中节点的隐含依赖性, 这种依赖组的存在导致网络能够发生级联失效并导致网络更加脆弱, 如在依赖组的密度较高或依赖组规模较大时网络会在遭受攻击时以一阶不连续相变的方式发生突然性的崩溃[41], 同时依赖组的规模和分布[41]以及依赖规则[42-44] 等特征的变化对网络鲁棒性和级联失效有着不同程度的影响.
在关于依赖组的研究中, 大多数模型基于一个简单的假设: 一个节点的失效可以直接导致其依赖节点的全部失效, 这种“节点对节点”的“强依赖”的假设导致了网络在遭受攻击时极为脆弱, 这一结果不但与现实中某些观测不符, 也不能充分描述部分复杂系统节点之间的依赖关系[45,46]. 在现实网络中, 一个节点的失效可能对其依赖节点的一部分功能造成影响, 而不是导致其完全失效[47,48]. 如金融系统中, 所属同一个投资者的公司之间相互依赖, 由于风险分散, 一家公司的倒闭不太可能直接导致投资者的其他公司全部破产, 而现实的情况是, 其余公司可能会因现金流的减少而收缩自己的业务. 因此与以往模型所描述的“节点对节点”依赖模式不同, 在一些普遍场景下系统节点之间的依赖关系体现在“节点对边”的相互作用上, 即一个节点失效后, 依赖组内的其余节点会因受到冲击而损失部分连接. 这种描述“节点对边”的相互作用机制可以描述更为一般化的场景, 如属于同一个依赖组的若干节点, 当其中一个节点失效的时候, 不会引起其余全部或部分节点立即失效, 而是失效节点给剩余依赖节点带来一定程度的冲击和影响.
本文引入了“弱依赖组”概念, 通过一个可调参数来控制同一依赖组内节点的相互依赖强度并研究依赖组的大小、分布和网络度分布对网络的级联失效的影响. 通过理论研究和数值模拟的结果, 发现改变依赖组的参数时, 网络能够发生连续和不连续的两种渗流相变, 且无标度网络在渗流过程中可能会出现双重相变现象.
2. 具有弱依赖组的网络级联失效模型
考虑一个由N个节点通过一些连接边所组成的网络, 节点的连接度服从概率分布函数P(k)[7], 同时网络中的每一个节点都随机地与其他节点形成一个包含M个节点的依赖组, 组内的节点互相依赖. 当依赖组中一个节点失效时, 其余节点会受到一定程度的冲击, 影响强度由解耦系数α 控制, 即其余每个节点的每条边都以α 的概率保留下来, 而以
1−α 的概率被删除. 参数α在描述失效节点对其余依赖节点所造成破坏性的同时, 也描述了一种“节点对边”的作用机理, 在描述依赖节点之间的相互作用上更具有一般性的意义, 当α→1 时, 依赖组内节点的耦合强度最弱, 一个节点的失效不能对组内其余节点造成任何影响, 此时模型与普通的座渗流等价; 而当α→0 时, 依赖组内节点的耦合强度最强, 一个节点的失效就能造成组内所有节点全部失效, 此时模型与文献[41]中强依赖组的模型等价. 通过调整模型参数α 可以研究不同耦合强度对网络级联失效动力学和鲁棒性的影响[49,50].网络级联失效通过随机攻击网络中比例为1–p的节点来触发, 这些被攻击节点及其连接被从网络中全部删除, 进而导致网络发生破碎. 在网络中, 能够连接到巨分支的节点视为功能节点, 其他节点视为失效节点. 由于弱依赖组内节点之间存在一定程度的依赖性, 网络中一个节点的失效会造成两种影响: 1)组外影响, 节点失效会引起网络破碎, 进而导致组外部分节点不能通过该节点连接到网络巨分支而失效; 2)组内影响, 失效节点所在的依赖组内的其余节点会受到损害, 其每一条连接边都以1 – α的概率被删除, 即以α 的概率被保留. 当一个节点失效后, 组外影响导致失效能够跨依赖组传播, 从而使失效扩展到更大的范围, 而组内影响导致组内的其余节点的边受到破坏, 从而引起组内更多节点失效. 在这两种影响的交替作用下, 网络上的级联失效就会发生. 我们将这两种扩散过程称为“组间级联失效”和“组内级联失效”. 当这两个过程停止时, 网络不会产生新的失效节点而发生进一步的破碎, 所有依赖组内也不会有新的连接边被删除, 网络即达到稳态. 本文用网络稳态巨分支的相对大小S 来衡量网络的完整性, 随着初始保留节点比例p的降低, 网络巨分支的大小S会以二阶连续相变或一阶不连续相变的方式在一个临界点
pIIc 或pIc 处降低至零, 二阶相变的临界点pIIc 或一阶相变的临界点pIc 通常用于度量网络的鲁棒性[15,51,52].图1为具有弱依赖组的网络的级联失效示意图. 图中虚线椭圆内的节点形成一个依赖组, 实线表示网络内功能正常的连接边, 点状虚线表示被删除的连接边; 功能正常的节点用绿色标记, 失效节点用红色标记, 黄色节点表示功能正常但受到损坏的节点. 图1中(a)表示的是一个未受到攻击的网络; (b)网络中存在一个最初被删除的节点2; (c)由于节点2被删除, 其所有连接边失效; (d)由于节点1与网络巨分支之间的连接边被移除, 导致节点1脱离网络巨分支而失效; (e)由于节点7和节点1属于同一个弱依赖组, 相互依赖, 故节点1 的失效会导致节点7受损; (f)节点7 受到了一定的冲击, 节点7与节点6之间的连接边失效; (g)因节点7与节点6之间的连接边失效, 导致节点6脱离网络巨分支, 发生故障; (h)由于处在同一弱依赖组, 节点5, 4受节点6失效的影响而受损; (i)节点5受损导致节点3和节点5之间的连接边失效, 并且网络达到了最终的稳定状态.
3. 理 论
用概率生成函数理论来对模型求解. 在求解过程中, 采用度分布的生成函数
G0(x)=∑kP(k)xk 来生成网络的度分布, 即一个随机节点度值为k的概率是P(k); 类似地, 采用余度分布生成函数G1(x)=∑kkP(k)⟨k⟩xk−1 来生成一条随机边到达一个度值为k 的节点的概率, 即排除这条随机边之外到达节点的余度为k – 1的概率[53]. 在这两个生成函数中,⟨k⟩≡∑kP(k)k 表示网络的平均度. 为了求解稳态时网络巨分支的规模, 定义R为网络中任意一条随机边能够连接到巨分支的概率. 假定在网络中随机挑选一条边后, 沿着这条边的任意一端可到达一个度值为k的节点i, 若该节点所属依赖组内有t个失效节点, 他的每条边能够保留下来的概率为αt . 由于任意一个节点度值k服从概率分布P(k), 节点i的任意一条随机边能够连到巨分支的概率分布遵循αt[1−P(k)k⟨k⟩(1−αtR)k−1] . 根据生成函数G1(x) 并考虑k的概率分布, 可将一条随机边连接到巨分支的概率简写为αt[1−G1(1−αtR)] [54,55]. 如果将任意一个节点所在的依赖组内失效节点数t所服从的概率分布记为f(t), 考虑所有可能的t, 可以得到R 满足自洽方程:R=pM−1∑t=0αt[1−G1(1−αtR)]f(t)≡h(R). (1) 同理, 可以得到一个随机节点属于网络稳态巨分支的概率S的方程:
S=pM−1∑t=0[1−G0(1−αtR)]f(t). (2) 求解方程(1)和方程(2)需要用到概率分布函数f(t), 该函数可以通过概率R来获得. 当网络处在稳定状态时, 假如随机选择一个节点具有t个失效的依赖节点, 那么该节点还能保留下来的概率为
p[1−G0(1−αtR)] , 与其同在一个依赖组中的其他M–t–1个依赖节点全部存活的概率为{p[1−G0(1−αtR)]}M−t−1 . 此外, 依赖组内的失效节点分为两种, 一种是因遭受初始攻击或因组外节点失效导致其脱离巨分支而失效, 另一种是组内其他节点的失效导致其连接的减少而失效. 假设有s个节点的失效是由前一种原因造成的, 其发生概率为{1−p[1−G0(1−R)]}s , 那么t−s 个节点的失效是由后一种原因导致的, 概率为{p[G0(1−αsR)−G0(1−R)]}t−s . 因此, t的概率分布函数f(t) 为\begin{split} f(t) =\;&\binom{M-1}{t}\{p[1-G_0(1-\alpha^tR)]\}^{M-t-1}\\ &\times \sum\limits_{s=0}^t\binom{t}{s} \{\{1-p[1-G_0 (1- R)]\}^s\\ &\times\{p[G_0 (1-\alpha^s R)-G_0(1-R)]\}^{t-s}\}. \end{split} (3) 对于给定的度分布P(k)和弱依赖组大小M, 结合方程(3)并通过联立方程(1)和方程(2)可以得到网络中巨分支的最终大小S.
当
\alpha\rightarrow1 时, 依赖组内节点依赖强度趋向0, 节点之间无相互依赖性而不会发生级联失效, 随着保留节点比例p的降低, 网络以二阶连续相变的方式发生破碎[15]. 当\alpha\rightarrow0 时, 弱依赖组内节点之间的依赖强度最强, 一旦组内一个节点被删除, 组内其余节点将全部失效, 网络以一阶不连续相变的方式发生破碎. 因此, 弱依赖组内节点之间的依赖强度α的值决定着网络在遭受攻击时的破碎形式, 当α 超过临界点\alpha_{\rm{c}} 时, 网络的破碎方式可以从二阶渗流相变转变为一阶渗流相变. 在网络以二阶相变方式发生破碎的情况下, 当p 接近二阶渗流相变点p_{\rm{c}}^{{\rm{II}}} 时, 随机边属于稳态巨分支的概率R趋于零. 在R\equiv \varepsilon \rightarrow 0 和p\rightarrow p_{\rm{c}}^{{\rm{II}}} 时, 可以将(1)式进行泰勒展开表示为\begin{equation} h(\varepsilon)=h'(0)\varepsilon+\frac{1}{2} h''(0)\varepsilon^2+O(\varepsilon^3 )=\varepsilon. \end{equation} (4) 忽略高阶项, 当
p\rightarrow p_{\rm{c}}^{{\rm{II}}} 时,h'(0)=1 . 因此可以得到网络发生二阶渗流相变的条件\begin{equation} p_{\rm{c}}^{{\rm{II}}}\alpha^{2(M-1)} G_1'(1)=1. \end{equation} (5) 进而可以求出二阶渗流相变点:
\begin{equation} p_{\rm{c}}^{{\rm{II}}}= {1}/\big[{\alpha^{2(M-1)}G_1'(1)}\big]. \end{equation} (6) 当α = 1或M = 1时, 网络不存在弱依赖组, 其结果和普通网络的渗流结果相同, 二阶渗流相变点
p_{\rm{c}}^{{\rm{II}}}= {1}/{G_1'(1)} [17]. 同时还可以发现, 二阶渗流相变点p_{\rm{c}}^{{\rm{II}}} 随着弱依赖组大小M的增大而减小, 这意味着网络中弱依赖组规模越大, 网络就越脆弱. 这一结果与弱依赖多层网络的结果类似, 即二阶相变点与网络的层数和参数α有关[56].在一阶相变点
p_{\rm{c}}^{\rm{I}} 处, R从零跳跃到一个非平凡值R_{\rm{c }} , 由图2(a)可知此时y=h(R)-R 的曲线与直线y=0 相切:图 2 在不同初始保留节点比例p和解耦系数α时(1)式的图解 (a) 的结果; (b)\alpha=0.5 的结果; (c)\alpha=0.63 的结果. 本图中所用网络的度分布服从幂率分布\alpha=0.8 , 其中P(k)\sim k^{-\gamma} (k_{{\rm{min}}}\leqslant k \leqslant k_{{\rm{max}}}) ,k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 \gamma=2.3 Fig. 2. The graphical solutions of Eq. (1) for different values of p and α: (a) Result for ; (b) result for\alpha=0.5 ; (c) result for\alpha= 0.63 . For each panel, the degree distribution of networks follows a truncated power-law distribution\alpha= 0.8 P(k)\sim k^{-\gamma} with(k_{{\rm{min}}}\leqslant k \leqslant k_{{\rm{max}}}) ,k_{{\rm{min}}}=2 , andk_{{\rm{max}}}=141 \gamma=2.3 \begin{equation} \frac{{\rm{d}}h(R)}{{\rm{d}}R}\bigg|_{(R=R_{\rm{c}},p=p_{\rm{c}}^{\rm{I}} )}=1. \end{equation} (7) 通过联立方程(7)和方程(1), 可求得一阶渗流相变点
p_{\rm{c}}^{\rm{I }} , 本文采用数值方法求解.当
\alpha=\alpha_{\rm{c}} 时, 一阶和二阶渗流相变的条件同时满足, 即p_{\rm{c}}^{\rm{I}}=p_{\rm{c}}^{{\rm{II}}} . 在这种情况下, 当p\rightarrow p_{\rm{c}} (p_{\rm{c}}^{\rm{I}} 或p_{\rm{c}}^{{\rm{II}}}) 时, 则R_{\rm{c}}\rightarrow0 . 此时, (4)式简化为\begin{equation} \frac{1}{2} h''(0) \varepsilon^2 + O(\varepsilon^3)=0, \end{equation} (8) 即
h''(0)=0 , 结合p_{\rm{c}}={1}/\big[{\alpha_{\rm{c}}^{2(M-1)}G_1'(1)}\big] 可进一步推出\begin{equation} \alpha_{\rm{c}}^{2M+1} G_1''(1)+2(\alpha_{\rm{c}}^2-1)(M-1)G_0'(1)=0. \end{equation} (9) 通过对(9)式求解, 可以得到一阶相变和二阶相变发生转变的临界点
\alpha_{\rm{c}} . 在随机网络中, 节点度值服从泊松分布, 方程(9)可以简化为\begin{equation} \alpha_{\rm{c}}^{2M+1}\langle k \rangle+2(\alpha_{\rm{c}}^2-1)(M-1)=0. \end{equation} (10) 从(10)可以看出,
\alpha_{\rm{c}} 依赖于网络平均度和依赖组的规模M. 对于无标度网络, (9)式可以写成\begin{equation} \frac{\langle k(k-1)(k-2)\rangle}{\langle k\rangle^2}\alpha_{\rm{c}}^{2M+1} +2(\alpha_{\rm{c}}^2-1)(M-1)=0, \end{equation} (11) 其中括号
\langle\cdot\rangle 表示依据度分布对括号内的项取平均. 因此,\alpha_{\rm{c}} 不仅取决于依赖组的规模M, 而且依赖于无标度网络的度分布P(k).图2给出了度值服从幂率分布的无标度网络在参数p和α取不同值时
y=h(R)-R 的函数曲线, 该图给出了系统发生一阶相变和二阶相变的条件, 并证实了模型中存在一阶相变和二阶相变的转变点\alpha_{{\rm{c}}} . 图2(a)显示了当\alpha=0.5 (\alpha < \alpha_{\rm{c}} )时, 网络巨分支的解随节点保留概率p的变化, 当p < p_{\rm{c}}^{\rm{I}} 时, R 的解由原点给出. 当p=p_c^{\rm{I}} 时, R的解由曲线与横坐标轴的切点给出, 这说明网络巨分支在p_{\rm{c}}^{\rm{I }} 处存在跳跃. 图2(c)给出了当\alpha=0.8 (\alpha > \alpha_{\rm{c }} ) 时, (1)式的图解, 即网络巨分支的解随节点保留概率p的变化. 当p=p_{\rm{c}}^{{\rm{II}}} 时,y=h(R)-R 在R=0 处与直线y=0 相切, 证实了网络在此处发生了二阶连续渗流相变. 图2(b)展示了一个特别的现象, 当\alpha=0.63 时, 网络巨分支在p\approx0.8357 处先以二阶相变的形式涌现, 随着p的增加, 而后在p\approx 0.8773 处又以一阶相变的形式再次发生相变. 这说明当网络度分布服从幂率分布P(k)\sim k^{-\gamma} (k_{{\rm{min}}}\leqslant k \leqslant k_{{\rm{max}}}) 时, 网络存在双重相变的现象. 由于不论网络以什么形式破碎, 在相变点p_{\rm{c}}^{{\rm{II}}} 和p_{\rm{c}}^{\rm{I}} 处都依次满足条件(5)式和(7)式, 可以利用这些条件来确定双重相变和单次一阶相变之间的边界, 以及双重相变和单次二阶相变之间的边界.4. 数值模拟结果与分析
本文选择了两种网络模型, 即度分布满足泊松分布的ER随机网络和度分布满足近似幂律的无标度网络, 来研究弱依赖组内节点的依赖强度、依赖组的规模和分布等不同特征对网络级联失效动力学和鲁棒性的影响, 并通过理论计算对数值模拟的结果予以验证. 对于ER随机网络, 采用随机向N个节点随机加边的方式生成. 而对于无标度网络, 首先生成一个近似满足幂率分布的度序列, 然后采用配置模型生成其对应的无标度网络. 在两种模型的生成过程中, 自环和重复连接被禁止.
4.1 弱依赖组规模均匀分布
为了描述具有弱依赖组的网络上的“组间级联”和“组内级联”两种失效过程对网络鲁棒性的影响, 本文以依赖组大小
M=3 和M=4 的两种情况为例, 研究当弱依赖组大小均匀分布时网络级联失效动力学的特点. 图3(a)和图3(b)分别给出了当M=3 和M=4 时, 平均度\langle k\rangle=4 的随机网络稳态巨分支的大小S随初始节点保留概率p的增加而变化情况. 图3(c)和图3(d)给出了平均度\langle k\rangle= 5 的随机网络, 当弱依赖组大小M分别为3和4时, 网络稳态巨分支规模S随p的增加而变化的情况. 从图3(a)和图3(b)可以发现, 无论是M= 3 还是M=4 的情况, 当α 较小时, 随着初始保留节点p的降低, 网络都是以一阶不连续相变的方式发生崩溃; 而对于较大的α值, 随着初始保留节点p的降低, 网络以二阶连续相变的形式逐步瓦解. 另外, 可以发现M=3 的网络的临界点小于M=4 的网络的临界点, 这意味着前者比后者具有更强的鲁棒性. 对于\langle k\rangle=5 的随机网络, 虽然定量的结果与\langle k\rangle=4 时存在差异, 但定性的结果是相似的. 这些结果证实了依赖组的规模和组内节点的依赖强度对网络的破碎形式有着重要影响, 依赖组规模越大, 网络的鲁棒性就越差, 同时依赖强度对网络的鲁棒性也有重要影响, 当依赖强度较大时网络以一阶不连续相变的形式发生破碎, 而依赖强度较小时网络以二阶连续相变的形式发生破碎.图 3 依赖组规模 和M=3 的随机网络上的渗流相变 (a), (b)平均度为M=4 的随机网络, 弱依赖组大小分别为\langle k\rangle=4 和M=3 的结果; (c), (d)在M=4 时, 弱依赖组平均大小分别为\langle k\rangle=5 和M=3 的随机网络的结果. 图中网络规模M=4 , 实线表示由(1)式和(2)式得到的理论预测, 而垂直虚线分别表示由(5)式和(7)式所预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均N=10^5 Fig. 3. Simulation results for the percolation transitions on random networks with group sizes andM=3 : (a), (b) Results for random networks with the interdependency group sizesM=4 andM=3 respectively, where the average degreeM=4 is set to 4 for both panels; (c), (d) results for random networks with the interdependency group sizes\langle k\rangle andM=3 respectively, where average degreeM=4 is set to 5 for both panels. The network size is\langle k\rangle . The solid lines behind the symbols denote the theoretical predictions that were obtained by Eqs.(1) and (2). The vertical dashed lines denote the first-order and second-order percolation transition points predicted by Eqs.(7) and (5), respectively. The simulation results were averaged fromN=10^5 realizations of different initial conditions.100 图3还提供了理论预测的结果, 理论预测的网络巨分支规模S与模拟结果非常符合, 另外, 数值模拟所得到的网络巨分支的规模S在理论预测的相变点处消失为0. 这些结果证实了本文理论的正确性. 此外, 在不同初始条件下, 初始节点保留概率p值对于极大簇的规模分布有着非常重要的影响. 在p值远离临界点时, 极大簇的规模分布是非常均匀的, 而在p值处于临界点附近时, 极大簇规模的涨落会非常大. 因此除了理论预测的方法之外, 还可以通过数值模拟来获得极大簇方差的大小来定位网络发生破碎的临界点[47].
图4给出了弱依赖组大小分布均匀时, 无标度网络的巨分支规模S随初始保留节点比例p的变化. 此时, 网络节点度服从
P(k)\sim k^ {-\gamma} (k_{{\rm{min}}}\leqslant k\leqslant k_{{\rm{max}}}) , 其中k_{{\rm{max}}} 和k_{{\rm{min}}} 分别是节点度的上下限,\gamma 是幂律指数, 其平均度\langle k\rangle = \displaystyle\sum\nolimits^{k_{{\rm{max}}}}_{k_{{\rm{min}}}}P(k)k . 与ER随机网络的结果类似, 对于较小的α 值, 网络以一阶不连续的方式发生崩溃, 对于较大的α值, 网络以二阶连续的方式逐步瓦解. 两种渗流相变形式由临界点\alpha_{\rm{c}} 隔开. 且当M=3 时, 网络的渗流相变点小于M=4 的网络的渗流相变点, 这意味着M=3 的网络比M=4 的网络具有更强的鲁棒性. 有趣的是, 我们还发现当α增大到一定程度时无标度网络会发生双重相变, 如\alpha=0.63 . 对于双重相变, 网络首先以二阶连续相变的方式渗流, 然后随着p的增加又会发生一阶不连续相变. 随着α 的增加, 网络破碎形式转变为单个二阶连续相变; 随着α的减小, 二阶连续相变消失, 网络的破碎形式转变为单个一阶不连续相变.图 4 依赖组规模 和M=3 的无标度网络渗流相变 (a), (b)在M=4 (\langle k\rangle=4 ,k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 )时, 弱依赖组平均大小分别为\gamma=2.5 和M=3 时的结果; (c), (d)在M=4 (\langle k\rangle=5 ,k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 )时, 弱依赖组规模分别为\gamma=2.3 和M=3 时的结果. 网络规模M=4 . 图中实线表示由(1)式和(2)式得到的理论预测. 垂直虚线分别表示(7)式和(5)式预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均N=10^5 Fig. 4. Simulation results for the percolation transitions on scale-free networks with the interdependency group sizes andM=3 : (a), (b) Results for scale-free networks with the interdependency group sizesM=4 andM=3 respectively, where the average degreeM=4 is\langle k\rangle with4 ,k_{{\rm{min}}}=2 andk_{{\rm{max}}}=63 ; (c), (d) results for scale-free networks with the interdependency group sizes\gamma=2.5 andM=3 respectively, where the average degreeM=4 is\langle k\rangle with5 ,k_{{\rm{min}}}=2 andk_{{\rm{max}}}=141 . The network size is\gamma=2.3 . The solid lines behind the symbols denote the theoretical predictions by Eqs.(1) and (2), and the vertical dashed lines denote the first-order or second-order percolation transition points predicted by Eqs.(7) and (5), respectively. The simulation results are averaged fromN=10^5 realizations of different initial conditions100 此外, 图4的结果还表明, 双重相变发生于度分布异质性较强的无标度网络且依赖组内节点依赖强度适中时. 这是由于无标度网络节点的度分布异质性较强, 当网络受到攻击而发生级联失效时, 较多度值较小的节点会率先脱离网络的巨分支而失效, 当失效节点的规模发散时网络就会发生一阶不连续的相变, 而此时网络中度值较大的节点仍然能够形成稳定的簇结构, 需要更大的攻击强度才能够将其彻底破坏. 当依赖组中节点的依赖强度较大时, 度值较大的中心节点很容易受到度值较小节点的影响而失效, 无法形成规模较大的簇结构来对抗级联失效, 此时网络只能发生一阶相变; 在当依赖组中节点的依赖强度较弱时, 节点会因相互之间的依赖太弱无法形成较大规模的级联失效而发生二阶连续相变.
图5(a)和图5(b)分别给出了在
M=3 和M= 4 时随机网络的渗流相变点p_{\rm{c}} 随α变化的情况, 其中不同颜色的曲线表示不同网络平均度的结果. 图中\alpha_{\rm{c}} 表示一阶相变和二阶相变的转变点, 其值由弱依赖组的大小M和网络平均度\langle k\rangle 决定. 图5(c)和图5(d)分别给出了在M=3 和M=4 时无标度网络的渗流相变点p^{{\rm{I}}}_{\rm{c}} 或p^{{\rm{II}}}_{\rm{c}} 随α的变化.图 5 随机网络和无标度网络在弱依赖组规模服从平均分布时渗流相变点 与α的关系 (a)依赖组规模p_c 的随机网络的渗流相变点M=3 或p_{\rm{c}}^{\rm{I}} 与α的关系, 其中平均度p_{\rm{c}}^{{\rm{II}}} 为4, 5和6; (b)依赖组规模\langle k\rangle 的随机网络的结果; (c)平均度M=4 为4, 5和6 (分别对应度分布的幂律指数\langle k\rangle 为2.5, 2.3和2.1)且依赖组规模\gamma 的无标度网络的结果; (d)依赖组规模M=3 的无标度网络的结果. 图中实线为一阶渗流相变点M=4 , 虚线为二阶渗流相变点p_{\rm{c}}^{\rm{I}} . 模拟结果来自于100次不同初始条件下的平均p_{\rm{c}}^{{\rm{II}}} Fig. 5. The percolation transition point as functions of α on random networks and scale-free networks with uniform interdependency group size: (a) Results for random networks with the interdependency group sizep_c , where the average degreeM=3 is 4, 5 and 6 from high to low; (b) results for random networks with the group size\langle k\rangle , where the average degreeM=4 is 4, 5 and 6 from high to low; (c) results for scale-free networks with the group size\langle k\rangle , where the average degreeM=3 is 4, 5 and 6 (corresponding to a power-law exponent of degree distribution\langle k\rangle is 2.5, 2.3 and 2.1, respectively); (d) results for scale-free networks with the interdependency group size\gamma , where the average degreeM=4 is 4, 5 and 6. The solid line is the first-order percolation transition point\langle k\rangle , and the dashed line is the second-order percolation transition pointp_{\rm{c}}^{\rm{I}} . The simulation results were averaged from 100 realizations of different initial conditionsp_{\rm{c}}^{{\rm{II}}} 从中可以发现, 相变的方式以
\alpha_{\rm{c }} 为界分为不连续和连续两种, 但\alpha_{\rm{c}} 的值取决于弱依赖组的大小M 和网络度分布P(k). 另外, 当无标度网络出现双重相变情况时, 网络先以二阶连续相变的方式渗流, 因此网络渗流相变点为二阶相变点p_{\rm{c}}^{{\rm{II}}} . 同时也可以发现, 在大多数情况下, 当依赖组内节点依赖强度较弱时, 随机网络鲁棒性比无标度网络差, 而当依赖组内节点依赖强度较强时, 无标度网络鲁棒性比随机网络差. 这是由于在随机攻击时, 无标度网络中度值较大的中心节点不容易受到攻击, 在依赖组内节点依赖强度较弱时也不容易受到其他节点的影响而受到损害, 然而在依赖强度较大时, 无标度网络中度值较大的节点容易受到组内其他节点失效的影响而造成大量连接被删除, 从而对网络整体的连通性造成较大的影响. 这一结果表明, 无标度网络的鲁棒性与依赖组内节点的依赖强度有着密切的关联.4.2 弱依赖组大小泊松分布
考虑到现实中一些复杂系统节点所形成的依赖组的规模并不一定是均匀的, 本节以弱依赖组规模服从泊松分布的情况为例, 用ER随机网络和无标度网络来研究弱依赖组规模的非均匀分布对网络级联失效动力学的影响. 考虑弱依赖组规模服从泊松分布函数
Q(M=i)=\dfrac{\langle M \rangle^i}{i!}{\rm{e}}^{-\langle M \rangle}, i=1, 2, \cdots, M_{{\rm{max}}} , 其中\langle M \rangle 表示依赖组的平均规模. 此时方程(1) 转换成\begin{equation} R = p\sum\limits_{M=1}^{M_{{\rm{max}}}}Q(M)\sum\limits_{t=0}^{M-1}\alpha^t[1-G_1(1-\alpha^tR)]f(t) \equiv h(R). \end{equation} (12) 同理, 方程(2)要转换成
\begin{equation} S=p\sum\limits_{M=1}^{M_{{\rm{max}}}}Q(M)\sum\limits_{t=0}^{M-1}[1-G_0(1-\alpha^tR)]f(t). \end{equation} (13) 关于渗流相变点的求值与上文所述的弱依赖组平均分布的求值方法相同, 可得到二阶渗流相变点:
\begin{equation} p_{\rm{c}}^{{\rm{II}}}=\frac{1}{\sum_{M=1}^{M_{{\rm{max}}}}Q(M)\alpha^{2M-2}G_1'(1)}. \end{equation} (14) 根据方程(14)可以发现, 当网络中弱依赖组规模分布服从泊松分布时, 网络的二阶渗流相变点与依赖组平均规模
\langle M \rangle 相关, 二阶渗流相变点随着弱依赖组平均规模\langle M \rangle 的增大而减小, 这意味着网络中弱依赖组平均规模越大, 网络鲁棒性就越差.当网络发生一阶渗流相变时, R在相变点
p_{\rm{c}}^{{\rm{I}}} 处从零跳跃到一个非零值R_{\rm{c}} , 此时方程(7)同样成立, 通过联立方程(7) 和方程(12)求得弱依赖组规模服从泊松分布时的一阶渗流相变点p_{\rm{c}}^{{\rm{I}}} .与弱依赖组规模平均分布时一阶与二阶渗流相变的临界点
\alpha_{\rm{c}} 满足的条件相同, 当p\rightarrow p_{\rm{c}} 和\alpha= \alpha_{\rm{c}} 时,h''(0)=0 , 此时得到\begin{split} &\sum\limits_{M=1}^{M_{{\rm{max}}}}Q(M)\{\alpha_c^{3M-3}G_1'' (1)+2p_{\rm{c}} (M-1)\\ & \times (\alpha_{\rm{c}}^{3M-4}-\alpha_{\rm{c}}^{3M-6})G_0'(1) G_1'(1)\}=0. \end{split} (15) 在度分布服从泊松分布的随机网络中, 方程(15)可简化为
\begin{split} &\sum\limits_{M=1}^{M_{{\rm{max}}}}Q(M) [\alpha_{\rm{c}}^{3M-3}+2p_{\rm{c}} \\ & \times (M-1)(\alpha_{\rm{c}}^{3M-4}-\alpha_{\rm{c}}^{3M-6})]=0. \end{split} (16) 在度分布近似幂率分布的无标度网络中, 方程(15)可简化为
\begin{split} &\sum\limits_{M=1}^{M_{{\rm{max}}}}Q(M)[\frac{\langle k(k-1)(k-2)\rangle}{\langle k\rangle \langle k(k-1)\rangle}\alpha_{\rm{c}}^{3M-3}\\ &+2p_{\rm{c}}(M-1)(\alpha_{\rm{c}}^{3M-4}-\alpha_{\rm{c}}^{3M-6})]=0. \end{split} (17) 由于此时
p_{\rm{c}}^{{\rm{II}}}=p_{\rm{c}}^{\rm{I}} , 所以渗流相变点p_{\rm{c}} 的表达式等同于方程(14), 再根据方程(16)和方程(17), 得到了和弱依赖组规模均匀分布相似的结论: 在随机网络中, 临界点\alpha_{\rm{c}} 与依赖组平均规模\langle M \rangle 和随机网络平均度\langle k\rangle 相关, 在无标度网络中, 临界点\alpha_{\rm{c}} 不仅受到依赖组平均规模\langle M \rangle 影响, 还受到度分布P(k)的影响.图6(a)和图6(b)分别给出了弱依赖组规模服从泊松分布且平均值
\langle M \rangle 分别为3和4时, 平均度\langle k\rangle= 4 随机ER网络稳态巨分支的大小S随p增加而变化的情况. 图6(c)和图6(d)分别给出了弱依赖组规模服从泊松分布且平均值\langle M \rangle 分别为3和4时, 平均度\langle k\rangle=5 的随机网络稳态巨分支的大小S随p增加而变化情况. 从图中可以发现在随机网络中, 对于较小的α值, 网络以一阶不连续相变方式渗流, 对于较大的α值, 网络以二阶连续相变方式渗流. 两种相变的形式在临界点\alpha_{\rm{c }} 处发生改变. 另外, 可以发现\langle M \rangle=3 时, 网络渗流相变点小于\langle M \rangle= 4 的网络, 这意味着依赖组平均规模较小的网络鲁棒性更强. 当\langle k\rangle=5 时, 弱依赖组平均规模\langle M \rangle=3 和\langle M \rangle=4 的随机网络也具备类似的现象. 图6还提供了理论预测, 结果与模拟结果基本符合.图 6 弱依赖组平均规模 和\langle M \rangle=3 的随机网络渗流相变的模拟结果 (a), (b) 平均度\langle M \rangle=4 的ER随机网络在依赖组平均规模分别为\langle k\rangle=4 和\langle M \rangle=3 时的结果; (c), (d)平均度\langle M \rangle=4 的ER 随机网络在依赖组平均规模分别为\langle k\rangle=5 和\langle M \rangle=3 时的结果. 网络规模设为\langle M \rangle=4 . 图中实线表示由(12)式和(13)式得到的理论预测. 垂直虚线分别表示理论预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均N=10^5 Fig. 6. Simulation results for percolation transitions on random networks with average group sizes and\langle M\rangle=3 : (a), (b) Results for random networks with the average group sizes\langle M\rangle=4 and\langle M\rangle=3 , respectively, where the average degree\langle M\rangle=4 is\langle k\rangle ; (c), (d) results for random networks with the average group sizes4 and\langle M\rangle=3 respectively, where\langle M\rangle=4 is 5. The network size is\langle k\rangle . The solid lines behind the symbols denote the theoretical predictions that were obtained by Eqs. (12) and (13), and the vertical dashed lines denote the first-order and second-order percolation transition points predicted by theory. The simulation results were averaged fromN=10^5 realizations of different initial conditions100 图7给出了弱依赖组大小分布为泊松分布
Q(M) , 且规模平均值\langle M\rangle=5 和\langle M\rangle=6 时, 无标度网络上的渗流相变. 从图7可以发现类似于随机网络的结果, 对于较小的α 值, 网络以一阶不连续相变方式渗流, 对于较大的α值, 网络以二阶连续相变方式渗流. 两种渗流相变的形式的转变点发生于\alpha_{\rm{c}} 处, 且\langle M \rangle=5 的网络的渗流相变点小于\langle M \rangle=6 的网络, 这意味着前者比后者具有更强的鲁棒性. 同时发现, 当弱依赖组大小服从泊松分布时, 无标度网络的双重相变现象依然存在. 图7还提供了理论预测, 理论结果与模拟结果的误差相对于随机网络较大, 主要存在两个原因. 首先对于无标度网络在给定度序列的情况下, 由于度值较大的节点的数量较小, 但这些节点对网络的鲁棒性产生了非常关键的作用, 不同网络连接结构的差异性会对模拟结果产生较大的影响, 因此模拟结果的误差较大. 第二, 极大簇在相变点附近有较大的涨落, 也会导致模拟结果存在一定的误差.图 7 弱依赖组平均规模 和\langle M\rangle=5 的无标度网络渗流相变的模拟结果 (a), (b)在\langle M\rangle=6 , 即\langle k\rangle=4 ,k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 时, 弱依赖组平均规模分别为\gamma=2.5 和\langle M\rangle=5 的结果; (c), (d) 在\langle M\rangle=6 , 即\langle k\rangle=5 ,k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 时, 弱依赖组平均规模分别为\gamma=2.3 和\langle M\rangle=5 的结果. 网络规模为\langle M\rangle=6 . 图中实线表示由(12)式和(13)式得到的理论预测. 垂直虚线分别表示预测的一阶和二阶渗流相变点N=10^5 Fig. 7. Simulation results for percolation transitions on scale-free network with average group sizes and\langle M\rangle=5 : (a), (b) Results for scale-free networks with average group sizes\langle M\rangle=6 and\langle M\rangle=5 respectively, where the average degree\langle M\rangle=6 is 4 with\langle k\rangle ,k_{{\rm{min}}}=2 andk_{{\rm{max}}}=63 ; (c), (d) results for scale-free networks with\gamma=2.5 and\langle M\rangle=5 respectively, where\langle M\rangle=6 is 5 with\langle k\rangle ,k_{{\rm{min}}}=2 andk_{{\rm{max}}}=141 . The network size is\gamma=2.3 . The solid lines behind the symbols denote the theoretical predictions by Eqs.(12) and (13), and the vertical dashed lines denote the first-order and second-order percolation transition points predicted by theory. The simulation results were averaged fromN=10^5 realizations of different initial conditions.100 图8(a)和图8(b)分别给出了当网络平均度
\langle k\rangle 不同时, 弱依赖组平均规模\langle M\rangle=3 和\langle M\rangle=4 的随机网络的渗流相变点p_{\rm{c}} 随α的变化. 图中网络渗流相变的形式在\alpha_{\rm{c}} 处发生改变, 且\alpha_{\rm{c}} 的大小由弱依赖组平均规模\langle M\rangle 和随机网络的平均度\langle k\rangle 决定. 图8(c)和图8(d) 分别给出了不同平均度\langle k\rangle 时, 弱依赖组平均大小\langle M\rangle=3 和\langle M\rangle=4 无标度网络的渗流相变点p_{\rm{c}} 随α的变化. 从中也可以发现类似依赖组大小平均分布的相变情况, 相变的方式以\alpha_{\rm{c}} 为临界点分为不连续和连续两种, 但临界点\alpha_{\rm{c}} 还取决于度分布的参数设置. 另外, 由于无标度网络出现双重相变情况时, 网络先以连续相变的形式渗流, 所以将p_{\rm{c}}^{{\rm{II}}} 视为此时渗流相变点p_{\rm{c}} . 通过比较图8和图5, 可以发现依赖组规模不均匀分布的网络比依赖组规模均匀分布的网络具有更强的鲁棒性.图 8 随机网络和无标度网络在弱依赖组规模服从泊松分布时渗流相变点 与α的关系 (a)p_{\rm{c}} 的随机网络的渗流相变点\langle M\rangle=3 (p_{\rm{c}} 或p_c^{\rm{I}} )与α的关系, 其中平均度p_{\rm{c}}^{{\rm{II}}} 分别为4, 5和6; (b)\langle k\rangle 的随机网络的结果; (c)平均度\langle M\rangle=4 分别为4, 5和6 (度分布的幂指数\langle k\rangle 分别为2.5, 2.3和2.1)且\gamma 的无标度网络的结果; (d)\langle M\rangle=3 的无标度网络的结果.图中实线为一阶渗流相变点\langle M\rangle=4 , 虚线为二阶渗流相变点p_c^{\rm{I}} p_c^{{\rm{II}}} Fig. 8. The percolation transition point as functions of α on random networks and scale-free networks when the interdependency groups size follows poisson distribution: (a) Results for random networks with the average interdependency group sizep_{\rm{c}} , where the average degree\langle M\rangle=3 is\langle k\rangle ,4 and5 from high to low; (b) the same results to panel (a) but for6 ; (c) results for scale-free networks with the average interdependency group size\langle M\rangle=4 , where the average degree\langle M\rangle=3 is 4, 5 and 6 from high to low (corresponding to a power-law exponent of degree distribution 2.5, 2.3 and 2.1, respectively);(d) the same results to panel (c) but for\langle k\rangle . The solid line is the first-order percolation transition point\langle M\rangle=4 , and the dashed line is the second-order percolation transition pointp_{\rm{c}}^{\rm{I}} p_{\rm{c}}^{{\rm{II}}} 5. 总结与展望
在一些真实的复杂系统中, 网络节点可能通过隐含的依赖关系而形成一些依赖组, 其中一个节点失效后会对依赖组内其余节点造成一定程度的损害, 进而会引发更多节点的失效而形成级联过程. 基于此, 本文提出了一种具有“弱依赖组”的复杂网络级联失效模型, 依赖组中节点的依赖方式由一种“节点对边”的相互作用机制决定, 与过去“节点到节点”的相互作用机制不同, 当依赖组内一个节点失效时, 组中其余节点都会受到一定程度的损害, 而非完全失效. 本文通过引入一个可调参数α来描述复杂网络中弱依赖组内节点的依赖强度, 并研究了网络节点所形成的弱依赖组的规模, 规模分布和依赖强度等几个特征对于复杂网络级联失效动力学和网络鲁棒性的影响.
由于“节点对边”的弱依赖机制的引入, 本文发现网络在级联失效过程中存在两种效应: 组间级联失效和组内级联失效. 具体而言, 组间级联失效表示一个节点的失效使依赖组外的部分节点不能通过失效节点连接到网络的巨分支而失效, 这些新增失效节点进而会造成更多节点的失效, 组间级联使得失效节点能够跨依赖组传播; 组内级联失效表示一个节点的失效会给依赖组内其余功能节点造成损害而失去部分连接, 从而使这些节点脱离网络巨分支而失效, 这些新增失效节点进而给组内其余节点造成更大的损害而诱发更多节点失效, 导致故障在依赖组内传播, 这一微观效应在之前强依赖机制的模型中是不存在的, 因在强依赖机制下, 一旦依赖组中一个节点被删除, 其余节点将会被完全删除. 由于组外级联失效的存在, 某一个依赖组内产生的失效会在扩张到其他依赖组后又传播到该依赖组内, 从而进一步导致组内更多节点失效. 而组内级联失效的存在打破了网络在部分节点在脱离网络巨分支后的平衡. 因此可以看出, 二者存在协同作用.
在两种级联动力学的共同作用下, 网络以一阶或二阶相变的形式瓦解, 且网络瓦解的形式由节点之间的依赖强度决定. 这说明依赖组内节点的依赖强度对网络发生级联失效的形式产生了至关重要的影响. 当节点之间的依赖性较弱时, 节点的失效不能对组内依赖节点造成较大的破坏性, 从而抑制了故障的传播, 使系统通过二阶渗流相变的形式发生破碎. 当节点之间的依赖性较强时, 节点的故障会以组内级联的方式传播到较为广泛的范围, 从而导致系统的突然崩溃. 这些结果表明, 组内级联失效过程影响着组间级联失效的过程, 并对网络级联失效的形式和网络鲁棒性起到了至关重要的作用.
数值模拟和理论分析的结果表明, 随着依赖组的规模和依赖组内节点的依赖强度的增大, 网络的鲁棒性会逐渐下降, 这是由于弱依赖组规模更小的网络可将组内级联限制在一个较小的范围内, 避免了失效在较大范围内扩散. 同时, 依赖组规模不均匀分布的网络比依赖组规模均匀分布的网络具有更强的鲁棒性, 这是因为依赖组规模分布不均匀的网络虽允许更大规模的依赖组出现, 但大规模依赖组出现的概率较低, 不容易受到攻击, 且攻击数量较多且规模较小的依赖组也不会造成太大范围的组内级联, 从而限制了组内级联的范围.
此外, 在大多数情况下, 当依赖组内节点依赖强度较弱时, 随机网络鲁棒性比无标度网络差, 而当依赖组内节点依赖强度较强时, 无标度网络鲁棒性比随机网络差. 这一结果表明, 存在依赖组的情况下, 度分布的异质性和依赖组规模分布的异质性都能够对系统的鲁棒性造成显著的影响.
最后, 本文还发现了无标度网络中在依赖组中节点依赖强度适中时会发生二重相变的现象, 这是由于无标度网络度分布的异质性较强, 当网络受到一定程度的攻击时所发生的级联失效只能破坏为数众多的度值较小的节点, 而对于度值较大的节点还能形成稳定的簇结构, 如将网络彻底破坏还需更强的攻击强度.
本文将“节点对边”作用机理引入到了具有依赖组的的级联失效模型中, 描述了节点之间更为一般的依赖机制, 不但对于理解依赖组的规模、规模分布和节点依赖强度等因素对网络级联失效的作用和级联失效的发生机理具有重要的科学意义, 而且对于通过依赖组的特征来对复杂网络级联失效的控制、预防干预策略的设计和鲁棒性的优化具有重要的启发意义. 本文所研究的依赖组基于组内节点的相互依赖性完全相同的假设, 然而在现实中依赖组内节点的相互影响可能不是完全等同的. 因此后续的研究仍然存在较大的拓展空间, 如依赖强度的非均匀性或非对等性对于网络的级联失效动力学是否具有显著的影响?此外, 依赖组内的节点是否能够形成某种结构, 依赖组的内节点的依赖结构是否对网络级联失效有较大的影响?这些问题都值得我们进一步进行实证或理论研究.
[1] Du W B, Zhou X L, Lordan O, Wang Z, Zhao C, Zhu Y B 2016 Transp. Res. Pt. e-Logist. Transp. Rev. 89 108
Google Scholar
[2] Tan F, Wu J, Xia Y, Chi K T 2014 Phys. Rev. E 89 062813
Google Scholar
[3] Liu R R, Jia C X, Lai Y C 2019 New J. Phys. 21 045002
Google Scholar
[4] Zhang L, Zeng G, Li D, Huang H-J, Stanley H E, Havlin S 2019 Proc. Natl. Acad. Sci. U.S.A. 116 8673
Google Scholar
[5] Ouyang M 2014 Reliab. Eng. Syst. Saf. 121 43
Google Scholar
[6] Newman M E, Watts D J, Strogatz S H 2002 Proc. Natl. Acad. Sci. U.S.A. 99 2566
Google Scholar
[7] Pastor-Satorras R, Castellano C, Van Mieghem P, Vespignani A 2015 Rev. Mod. Phys. 87 925
Google Scholar
[8] Zhao D, Wang L, Li S, Wang Z, Wang L, Gao B 2014 PLoS One 9 e112018
Google Scholar
[9] Zhao H 2004 Neural Networks 17 47
Google Scholar
[10] Liu X, Maiorino E, Halu A, Glass K, Prasad R B, Loscalzo J, Gao J, Sharma A 2020 Nat. Commun. 11 1
Google Scholar
[11] Shim P-S, Lee H K, Noh J D 2012 Phys. Rev. E 86 031113
Google Scholar
[12] Cho Y, Lee J, Herrmann H, Kahng B 2016 Phys. Rev. Lett. 116 025701
Google Scholar
[13] Hoppe K, Rodgers G J 2014 Phys. Rev. E 90 012815
Google Scholar
[14] Hackett A, Gleeson J P 2013 Phys. Rev. E 87 062801
Google Scholar
[15] Callaway D S, Newman M E, Strogatz S H, Watts D J 2000 Phys. Rev. Lett. 85 5468
Google Scholar
[16] Albert R, Jeong H, Barabási A L 2000 Nature 406 378
Google Scholar
[17] Cohen R, Erez K, Ben-Avraham D, Havlin S 2000 Phys. Rev. Lett. 85 4626
[18] Liu J G, Wang Z T, Dang Y Z 2005 Mod. Phys. Lett. B 19 785
Google Scholar
[19] Watts D J, Strogatz S H 1998 Nature 393 440
Google Scholar
[20] Albert R, Barabási A L 2002 Rev. Mod. Phys. 74 47
Google Scholar
[21] Crucitti P, Latora V, Marchiori M 2004 Phys. Rev. E 69 045104
[22] Moreno Y, Pastor-Satorras R, Vázquez A, Vespignani A 2003 EPL 62 292
Google Scholar
[23] Zhao L, Park K, Lai Y C, Ye N 2005 Phys. Rev. E 72 025104
Google Scholar
[24] Wang W X, Lai Y C 2009 Phys. Rev. E 80 036109
Google Scholar
[25] Wang J, Rong L, Zhang L 2009 Mod. Phys. Lett. B 23 1323
Google Scholar
[26] Wu X, Gu R, Ji Y, Stanley H E 2019 Phys. Rev. E 100 022309
[27] Motter A E, Lai Y C 2002 Phys. Rev. E 66 065102
Google Scholar
[28] Lehmann J, Bernasconi J 2010 Phys. Rev. E 81 031129
Google Scholar
[29] Wang W X, Chen G 2008 Phys. Rev. E 77 026101
[30] Yang R, Wang W X, Lai Y C, Chen G 2009 Phys. Rev. E 79 026112
Google Scholar
[31] Zhang S P, Huang Z G, Dong J Q, Eisenberg D, Seager T P, Lai Y C 2015 New J. Phys. 17 063029
Google Scholar
[32] Wang, J W 2012 Nonlinear Dyn. 70 1959
Google Scholar
[33] Motter A E 2004 Phys. Rev. Lett. 93 098701
Google Scholar
[34] Simonsen I, Buzna L, Peters K, Bornholdt S, Helbing D 2008 Phys. Rev. Lett. 100 218701
Google Scholar
[35] Xia Y, Fan J, Hill D 2010 Physica A 389 1281
Google Scholar
[36] Watts D J 2002 Proc. Natl. Acad. Sci. U.S.A. 99 5766
Google Scholar
[37] Gleeson J P, Cahalane D J 2007 Phys. Rev. E 75 056103
Google Scholar
[38] Dorogovtsev S N, Goltsev A V, Mendes J F F 2006 Phys. Rev. Lett. 96 040601
Google Scholar
[39] Baxter G J, Dorogovtsev S N, Goltsev A V, Mendes J F 2010 Phys. Rev. E 82 011103
Google Scholar
[40] Parshani R, Buldyrev S V, Havlin S 2011 Proc. Natl. Acad. Sci. U.S.A. 108 1007
Google Scholar
[41] Bashan A, Parshani R, Havlin S 2011 Phys. Rev. E 83 051127
Google Scholar
[42] Wang Z, Zhou D, Hu Y 2018 Phys. Rev. E 97 032306
[43] Niu D, Yuan X, Du M, Stanley H Eugene, Hu Y 2016 Phys. Rev. E 93 042312
[44] Wang H, Li M, Deng L, Wang B H 2018 Physica A 502 195
Google Scholar
[45] Wang H, Li M, Deng L, Wang B H 2015 PLoS One 10 e0126674
Google Scholar
[46] Li M, Wang B H 2014 Chin. Phys. B 23 076402
Google Scholar
[47] Chen M, Song M, Zhang M, Jin L, Gong X 2019 Int. J. Mod. Phys. C 30 1950043
Google Scholar
[48] Li M, Liu R R, Lyu L, Hu M B, Xu S, Zhang Y C 2021 Phys. Rep. 907 1
Google Scholar
[49] Cao Y Y, Liu R R, Jia C X, Wang B H 2021 Commun. Nonlinear Sci. Numer. Simul. 92 105492
Google Scholar
[50] Liu R R, Eisenberg D A, Seager T P, Lai Y C 2018 Sci. Rep. 8 2111
Google Scholar
[51] Kong L W, Li M, Liu R-R, Wang B H 2017 Phys. Rev. E 95 032301
[52] Wu J, Deng H Z, Tan Y J, Zhu D Z 2007 J. Phys. A-Math. Theor. 40 2665
Google Scholar
[53] Cohen R, Havlin S 2010 Complex Networks: Dtructure, Tobustness and Function (Cambridge: Cambridge University Press) p31
[54] Son S W, Bizhani G, Christensen C, Grassberger P, Paczuski M 2012 EPL 97 16006
Google Scholar
[55] Feng L, Monterola C P, Hu Y 2015 New J. Phys. 17 063025
Google Scholar
[56] Jiang W J, Liu R R, Jia C X 2020 Complexity 2020 3578736
-
图 2 在不同初始保留节点比例p和解耦系数α时(1)式的图解 (a)
\alpha=0.5 的结果; (b)\alpha=0.63 的结果; (c)\alpha=0.8 的结果. 本图中所用网络的度分布服从幂率分布P(k)\sim k^{-\gamma} (k_{{\rm{min}}}\leqslant k \leqslant k_{{\rm{max}}}) , 其中k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 ,\gamma=2.3 Figure 2. The graphical solutions of Eq. (1) for different values of p and α: (a) Result for
\alpha=0.5 ; (b) result for\alpha= 0.63 ; (c) result for\alpha= 0.8 . For each panel, the degree distribution of networks follows a truncated power-law distributionP(k)\sim k^{-\gamma} (k_{{\rm{min}}}\leqslant k \leqslant k_{{\rm{max}}}) withk_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 , and\gamma=2.3 图 3 依赖组规模
M=3 和M=4 的随机网络上的渗流相变 (a), (b)平均度为\langle k\rangle=4 的随机网络, 弱依赖组大小分别为M=3 和M=4 的结果; (c), (d)在\langle k\rangle=5 时, 弱依赖组平均大小分别为M=3 和M=4 的随机网络的结果. 图中网络规模N=10^5 , 实线表示由(1)式和(2)式得到的理论预测, 而垂直虚线分别表示由(5)式和(7)式所预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均Figure 3. Simulation results for the percolation transitions on random networks with group sizes
M=3 andM=4 : (a), (b) Results for random networks with the interdependency group sizesM=3 andM=4 respectively, where the average degree\langle k\rangle is set to 4 for both panels; (c), (d) results for random networks with the interdependency group sizesM=3 andM=4 respectively, where average degree\langle k\rangle is set to 5 for both panels. The network size isN=10^5 . The solid lines behind the symbols denote the theoretical predictions that were obtained by Eqs.(1) and (2). The vertical dashed lines denote the first-order and second-order percolation transition points predicted by Eqs.(7) and (5), respectively. The simulation results were averaged from100 realizations of different initial conditions.图 4 依赖组规模
M=3 和M=4 的无标度网络渗流相变 (a), (b)在\langle k\rangle=4 (k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 ,\gamma=2.5 )时, 弱依赖组平均大小分别为M=3 和M=4 时的结果; (c), (d)在\langle k\rangle=5 (k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 ,\gamma=2.3 )时, 弱依赖组规模分别为M=3 和M=4 时的结果. 网络规模N=10^5 . 图中实线表示由(1)式和(2)式得到的理论预测. 垂直虚线分别表示(7)式和(5)式预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均Figure 4. Simulation results for the percolation transitions on scale-free networks with the interdependency group sizes
M=3 andM=4 : (a), (b) Results for scale-free networks with the interdependency group sizesM=3 andM=4 respectively, where the average degree\langle k\rangle is4 withk_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 and\gamma=2.5 ; (c), (d) results for scale-free networks with the interdependency group sizesM=3 andM=4 respectively, where the average degree\langle k\rangle is5 withk_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 and\gamma=2.3 . The network size isN=10^5 . The solid lines behind the symbols denote the theoretical predictions by Eqs.(1) and (2), and the vertical dashed lines denote the first-order or second-order percolation transition points predicted by Eqs.(7) and (5), respectively. The simulation results are averaged from100 realizations of different initial conditions图 5 随机网络和无标度网络在弱依赖组规模服从平均分布时渗流相变点
p_c 与α的关系 (a)依赖组规模M=3 的随机网络的渗流相变点p_{\rm{c}}^{\rm{I}} 或p_{\rm{c}}^{{\rm{II}}} 与α的关系, 其中平均度\langle k\rangle 为4, 5和6; (b)依赖组规模M=4 的随机网络的结果; (c)平均度\langle k\rangle 为4, 5和6 (分别对应度分布的幂律指数\gamma 为2.5, 2.3和2.1)且依赖组规模M=3 的无标度网络的结果; (d)依赖组规模M=4 的无标度网络的结果. 图中实线为一阶渗流相变点p_{\rm{c}}^{\rm{I}} , 虚线为二阶渗流相变点p_{\rm{c}}^{{\rm{II}}} . 模拟结果来自于100次不同初始条件下的平均Figure 5. The percolation transition point
p_c as functions of α on random networks and scale-free networks with uniform interdependency group size: (a) Results for random networks with the interdependency group sizeM=3 , where the average degree\langle k\rangle is 4, 5 and 6 from high to low; (b) results for random networks with the group sizeM=4 , where the average degree\langle k\rangle is 4, 5 and 6 from high to low; (c) results for scale-free networks with the group sizeM=3 , where the average degree\langle k\rangle is 4, 5 and 6 (corresponding to a power-law exponent of degree distribution\gamma is 2.5, 2.3 and 2.1, respectively); (d) results for scale-free networks with the interdependency group sizeM=4 , where the average degree\langle k\rangle is 4, 5 and 6. The solid line is the first-order percolation transition pointp_{\rm{c}}^{\rm{I}} , and the dashed line is the second-order percolation transition pointp_{\rm{c}}^{{\rm{II}}} . The simulation results were averaged from 100 realizations of different initial conditions图 6 弱依赖组平均规模
\langle M \rangle=3 和\langle M \rangle=4 的随机网络渗流相变的模拟结果 (a), (b) 平均度\langle k\rangle=4 的ER随机网络在依赖组平均规模分别为\langle M \rangle=3 和\langle M \rangle=4 时的结果; (c), (d)平均度\langle k\rangle=5 的ER 随机网络在依赖组平均规模分别为\langle M \rangle=3 和\langle M \rangle=4 时的结果. 网络规模设为N=10^5 . 图中实线表示由(12)式和(13)式得到的理论预测. 垂直虚线分别表示理论预测的一阶和二阶渗流相变点. 模拟结果来自于100次不同初始条件下的平均Figure 6. Simulation results for percolation transitions on random networks with average group sizes
\langle M\rangle=3 and\langle M\rangle=4 : (a), (b) Results for random networks with the average group sizes\langle M\rangle=3 and\langle M\rangle=4 , respectively, where the average degree\langle k\rangle is4 ; (c), (d) results for random networks with the average group sizes\langle M\rangle=3 and\langle M\rangle=4 respectively, where\langle k\rangle is 5. The network size isN=10^5 . The solid lines behind the symbols denote the theoretical predictions that were obtained by Eqs. (12) and (13), and the vertical dashed lines denote the first-order and second-order percolation transition points predicted by theory. The simulation results were averaged from100 realizations of different initial conditions图 7 弱依赖组平均规模
\langle M\rangle=5 和\langle M\rangle=6 的无标度网络渗流相变的模拟结果 (a), (b)在\langle k\rangle=4 , 即k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 ,\gamma=2.5 时, 弱依赖组平均规模分别为\langle M\rangle=5 和\langle M\rangle=6 的结果; (c), (d) 在\langle k\rangle=5 , 即k_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 ,\gamma=2.3 时, 弱依赖组平均规模分别为\langle M\rangle=5 和\langle M\rangle=6 的结果. 网络规模为N=10^5 . 图中实线表示由(12)式和(13)式得到的理论预测. 垂直虚线分别表示预测的一阶和二阶渗流相变点Figure 7. Simulation results for percolation transitions on scale-free network with average group sizes
\langle M\rangle=5 and\langle M\rangle=6 : (a), (b) Results for scale-free networks with average group sizes\langle M\rangle=5 and\langle M\rangle=6 respectively, where the average degree\langle k\rangle is 4 withk_{{\rm{min}}}=2 ,k_{{\rm{max}}}=63 and\gamma=2.5 ; (c), (d) results for scale-free networks with\langle M\rangle=5 and\langle M\rangle=6 respectively, where\langle k\rangle is 5 withk_{{\rm{min}}}=2 ,k_{{\rm{max}}}=141 and\gamma=2.3 . The network size isN=10^5 . The solid lines behind the symbols denote the theoretical predictions by Eqs.(12) and (13), and the vertical dashed lines denote the first-order and second-order percolation transition points predicted by theory. The simulation results were averaged from100 realizations of different initial conditions.图 8 随机网络和无标度网络在弱依赖组规模服从泊松分布时渗流相变点
p_{\rm{c}} 与α的关系 (a)\langle M\rangle=3 的随机网络的渗流相变点p_{\rm{c}} (p_c^{\rm{I}} 或p_{\rm{c}}^{{\rm{II}}} )与α的关系, 其中平均度\langle k\rangle 分别为4, 5和6; (b)\langle M\rangle=4 的随机网络的结果; (c)平均度\langle k\rangle 分别为4, 5和6 (度分布的幂指数\gamma 分别为2.5, 2.3和2.1)且\langle M\rangle=3 的无标度网络的结果; (d)\langle M\rangle=4 的无标度网络的结果.图中实线为一阶渗流相变点p_c^{\rm{I}} , 虚线为二阶渗流相变点p_c^{{\rm{II}}} Figure 8. The percolation transition point
p_{\rm{c}} as functions of α on random networks and scale-free networks when the interdependency groups size follows poisson distribution: (a) Results for random networks with the average interdependency group size\langle M\rangle=3 , where the average degree\langle k\rangle is4 ,5 and6 from high to low; (b) the same results to panel (a) but for\langle M\rangle=4 ; (c) results for scale-free networks with the average interdependency group size\langle M\rangle=3 , where the average degree\langle k\rangle is 4, 5 and 6 from high to low (corresponding to a power-law exponent of degree distribution 2.5, 2.3 and 2.1, respectively);(d) the same results to panel (c) but for\langle M\rangle=4 . The solid line is the first-order percolation transition pointp_{\rm{c}}^{\rm{I}} , and the dashed line is the second-order percolation transition pointp_{\rm{c}}^{{\rm{II}}} -
[1] Du W B, Zhou X L, Lordan O, Wang Z, Zhao C, Zhu Y B 2016 Transp. Res. Pt. e-Logist. Transp. Rev. 89 108
Google Scholar
[2] Tan F, Wu J, Xia Y, Chi K T 2014 Phys. Rev. E 89 062813
Google Scholar
[3] Liu R R, Jia C X, Lai Y C 2019 New J. Phys. 21 045002
Google Scholar
[4] Zhang L, Zeng G, Li D, Huang H-J, Stanley H E, Havlin S 2019 Proc. Natl. Acad. Sci. U.S.A. 116 8673
Google Scholar
[5] Ouyang M 2014 Reliab. Eng. Syst. Saf. 121 43
Google Scholar
[6] Newman M E, Watts D J, Strogatz S H 2002 Proc. Natl. Acad. Sci. U.S.A. 99 2566
Google Scholar
[7] Pastor-Satorras R, Castellano C, Van Mieghem P, Vespignani A 2015 Rev. Mod. Phys. 87 925
Google Scholar
[8] Zhao D, Wang L, Li S, Wang Z, Wang L, Gao B 2014 PLoS One 9 e112018
Google Scholar
[9] Zhao H 2004 Neural Networks 17 47
Google Scholar
[10] Liu X, Maiorino E, Halu A, Glass K, Prasad R B, Loscalzo J, Gao J, Sharma A 2020 Nat. Commun. 11 1
Google Scholar
[11] Shim P-S, Lee H K, Noh J D 2012 Phys. Rev. E 86 031113
Google Scholar
[12] Cho Y, Lee J, Herrmann H, Kahng B 2016 Phys. Rev. Lett. 116 025701
Google Scholar
[13] Hoppe K, Rodgers G J 2014 Phys. Rev. E 90 012815
Google Scholar
[14] Hackett A, Gleeson J P 2013 Phys. Rev. E 87 062801
Google Scholar
[15] Callaway D S, Newman M E, Strogatz S H, Watts D J 2000 Phys. Rev. Lett. 85 5468
Google Scholar
[16] Albert R, Jeong H, Barabási A L 2000 Nature 406 378
Google Scholar
[17] Cohen R, Erez K, Ben-Avraham D, Havlin S 2000 Phys. Rev. Lett. 85 4626
[18] Liu J G, Wang Z T, Dang Y Z 2005 Mod. Phys. Lett. B 19 785
Google Scholar
[19] Watts D J, Strogatz S H 1998 Nature 393 440
Google Scholar
[20] Albert R, Barabási A L 2002 Rev. Mod. Phys. 74 47
Google Scholar
[21] Crucitti P, Latora V, Marchiori M 2004 Phys. Rev. E 69 045104
[22] Moreno Y, Pastor-Satorras R, Vázquez A, Vespignani A 2003 EPL 62 292
Google Scholar
[23] Zhao L, Park K, Lai Y C, Ye N 2005 Phys. Rev. E 72 025104
Google Scholar
[24] Wang W X, Lai Y C 2009 Phys. Rev. E 80 036109
Google Scholar
[25] Wang J, Rong L, Zhang L 2009 Mod. Phys. Lett. B 23 1323
Google Scholar
[26] Wu X, Gu R, Ji Y, Stanley H E 2019 Phys. Rev. E 100 022309
[27] Motter A E, Lai Y C 2002 Phys. Rev. E 66 065102
Google Scholar
[28] Lehmann J, Bernasconi J 2010 Phys. Rev. E 81 031129
Google Scholar
[29] Wang W X, Chen G 2008 Phys. Rev. E 77 026101
[30] Yang R, Wang W X, Lai Y C, Chen G 2009 Phys. Rev. E 79 026112
Google Scholar
[31] Zhang S P, Huang Z G, Dong J Q, Eisenberg D, Seager T P, Lai Y C 2015 New J. Phys. 17 063029
Google Scholar
[32] Wang, J W 2012 Nonlinear Dyn. 70 1959
Google Scholar
[33] Motter A E 2004 Phys. Rev. Lett. 93 098701
Google Scholar
[34] Simonsen I, Buzna L, Peters K, Bornholdt S, Helbing D 2008 Phys. Rev. Lett. 100 218701
Google Scholar
[35] Xia Y, Fan J, Hill D 2010 Physica A 389 1281
Google Scholar
[36] Watts D J 2002 Proc. Natl. Acad. Sci. U.S.A. 99 5766
Google Scholar
[37] Gleeson J P, Cahalane D J 2007 Phys. Rev. E 75 056103
Google Scholar
[38] Dorogovtsev S N, Goltsev A V, Mendes J F F 2006 Phys. Rev. Lett. 96 040601
Google Scholar
[39] Baxter G J, Dorogovtsev S N, Goltsev A V, Mendes J F 2010 Phys. Rev. E 82 011103
Google Scholar
[40] Parshani R, Buldyrev S V, Havlin S 2011 Proc. Natl. Acad. Sci. U.S.A. 108 1007
Google Scholar
[41] Bashan A, Parshani R, Havlin S 2011 Phys. Rev. E 83 051127
Google Scholar
[42] Wang Z, Zhou D, Hu Y 2018 Phys. Rev. E 97 032306
[43] Niu D, Yuan X, Du M, Stanley H Eugene, Hu Y 2016 Phys. Rev. E 93 042312
[44] Wang H, Li M, Deng L, Wang B H 2018 Physica A 502 195
Google Scholar
[45] Wang H, Li M, Deng L, Wang B H 2015 PLoS One 10 e0126674
Google Scholar
[46] Li M, Wang B H 2014 Chin. Phys. B 23 076402
Google Scholar
[47] Chen M, Song M, Zhang M, Jin L, Gong X 2019 Int. J. Mod. Phys. C 30 1950043
Google Scholar
[48] Li M, Liu R R, Lyu L, Hu M B, Xu S, Zhang Y C 2021 Phys. Rep. 907 1
Google Scholar
[49] Cao Y Y, Liu R R, Jia C X, Wang B H 2021 Commun. Nonlinear Sci. Numer. Simul. 92 105492
Google Scholar
[50] Liu R R, Eisenberg D A, Seager T P, Lai Y C 2018 Sci. Rep. 8 2111
Google Scholar
[51] Kong L W, Li M, Liu R-R, Wang B H 2017 Phys. Rev. E 95 032301
[52] Wu J, Deng H Z, Tan Y J, Zhu D Z 2007 J. Phys. A-Math. Theor. 40 2665
Google Scholar
[53] Cohen R, Havlin S 2010 Complex Networks: Dtructure, Tobustness and Function (Cambridge: Cambridge University Press) p31
[54] Son S W, Bizhani G, Christensen C, Grassberger P, Paczuski M 2012 EPL 97 16006
Google Scholar
[55] Feng L, Monterola C P, Hu Y 2015 New J. Phys. 17 063025
Google Scholar
[56] Jiang W J, Liu R R, Jia C X 2020 Complexity 2020 3578736
期刊类型引用(2)
1. 李耀华,刘肖迪,路向阳,王洁. 基于相依失效的“空-基-地”编队可靠性分析. 兵器装备工程学报. 2024(07): 19-27 . 百度学术
2. 李静,路庆昌,徐鹏程,王琴,王世鑫. 考虑站点耦合关系的大型城市轨道交通网络级联失效分析. 浙江大学学报(工学版). 2024(09): 1945-1955 . 百度学术
其他类型引用(10)
Catalog
Metrics
- Abstract views: 5856
- PDF Downloads: 119
- Cited By: 12