基于链路随机打击的后勤保障网络毁伤过程研究
邓 群
(陆军勤务学院油料系， 重庆 401331)
摘 要： 以ER(Erdös-Rényi)网络与WS(Watts-Strogtz)网络为例，对链路随机打击下后勤保障网络的毁伤过程进行了仿真与论述。通过分析仿真结果，提出后勤保障网络在毁伤过程中具有明显的阶段性特征，不同修复速率下的毁伤过程存在明显的阶跃现象，并且阶跃点的分布呈现出显著的幂次规律。尝试以新的视角探索后勤保障网络毁伤过程蕴含的动力学特征，对于提高战时后勤保障网络的战场生存能力具有一定的现实意义。
关键词： 后勤保障； 复杂网络； 系统动力学； 毁伤过程； 抗毁性
与计算机网络、社交网络、交通网络等其他复杂网络不同，后勤保障网络隐含的战争背景，意味着其可能面临更多频繁的主动攻击，而非偶尔的自然灾害或不可预料的系统风险。当后勤保障网络遭受外部攻击时，网络内部的各个要素会主动地相互协同，进而完成系统修复，在这一过程中包含了大量的非线性相互作用，甚至还有可能在宏观层面发生涌现现象。
随着复杂网络理论的兴起，后勤保障领域相关研究逐渐从以往对保障体系的定性论述深入到对保障网络的定量分析。其中，作为评判复杂网络的关键指标，抗毁性始终是后勤保障网络研究的重要课题之一。近年来，学术界关于复杂网络抗毁性的理论成果日益丰富完善，相关研究的焦点主要集中在2个方面：一是关于网络节点[1-3]与链路[4-6]的针对性研究，其本质上是对复杂网络组成要素的度、介数或其他改进指标的考察[7-10]，通过对比不同节点或链路失效对网络的影响，实现对其重要性的评判；二是对复杂网络级联失效[11-13]的探讨，其核心思想是从一个动态的系统性视角，以复杂网络的连通度为落脚点[14-15]，更加注重分析由节点或链路失效导致的连锁反应，进而从不同攻击策略[16-17]、不同负载方案[18]以及不同信息条件[19]等角度对复杂网络抗毁性研究加以细化和深入。
总体来看，复杂网络抗毁性研究从理论上进一步揭示了现实世界网络系统的规律和特性[20-22]，但仍存在几点不足：一是大多数研究单纯考虑网络损毁，忽略了网络节点与链路的修复重建，导致研究结果相对单调、片面，具有较强的主观倾向；二是对于网络节点与链路的刻画过于简单，往往认为节点“非有即无”、链路“非通即断”，很少考虑网络要素不完全损毁的情况，相关模型描述过于苛刻，与实际情况存在一定偏差；三是部分研究必须以假定网络无限大为前提，否则很难凸显复杂网络的特性，而现实世界中很多网络是有限规模的，这就使得部分理论成果转化应用效果不尽理想。此外，复杂网络组成要素间关系复杂，难以从宏观角度加以描述，以往相关研究大多数采用具体的实(算)例网络对提出的技术方法进行验证，缺乏对复杂网络共性规律的认识。
黍子对土壤和茬口要求不严，除严重盐碱地及低洼易涝地以外均可种植。在粘性土壤种黍子，由于土壤结构紧密，扎根困难，加之黍粒小，内含营养物质少，造成缺苗或死苗现象，但这种土壤保水保肥力强，对黍子后期生长有利。沙性土壤结构疏松，土温上升快，利于黍苗顶土，较适宜黍子生长，但后期易脱肥旱哀，要加强肥料供应。
综上所述，不论是比表面积的测试还是吸光度的测试，硅胶与硅酸盐复合吸附剂材料的性能是最好的，其次便是硅胶。经过一系列的测试可以得出硅酸盐及其复合材料吸附剂具有良好的吸附性能。
综上所述，笔者以链路随机打击下的后勤保障网络毁伤过程为研究对象，运用系统动力学方法，结合复杂网络特点，重新审视后勤保障网络抗毁性，旨在探索简单拓扑规则下后勤保障网络毁伤过程的规律特点，从而为复杂网络抗毁性研究提供新的思路。
1 后勤保障网络毁伤过程的动力学描述
1.1 后勤保障网络模型
利用无向图D n =(V n ,E n )完整描述某时刻的后勤保障网络。其中：
为该时刻全部节点的集合；
为该时刻全部边的集合。需要特别说明的是，e ij 与e ji 表示同一条边，在集合E n 中均记为
且仅记1次。
令
表示某时刻该节点被摧毁；
表示某时刻该节点完好；
表示某时刻该节点被摧毁后未被完全修复。当
时，意味着该节点暂时无法提供后勤保障功能。
对网络链路的持续摧毁行为从形式上可分为2类：1)当敌方军事优势明显时，可从始至终对后勤保障网络保持一定强度的打击。此种摧毁行为可抽象为持续按数量m 去除链路(式(2)、(4))，即每一轮打击均从剩余未被完全摧毁的链路中去除m 条链路；2)当敌方军事实力并不明显占优时，则难以长期保持较高的打击强度。此时，可将摧毁行为抽象为持续按比例α 去除链路(式(3)、(4))，即未被完全摧毁的链路在每一轮打击中被去除的概率恒为α 。
令Y 表示后勤保障网络的状态(0≤Y ≤M )，其中M 为后勤保障网络中链路总数。Y =M ，表示网络状态完好；Y =0，表示网络被彻底摧毁。限于篇幅，本文假定网络中所有节点始终完好，即不考虑对网络节点的打击。由此可知，若单纯考虑针对网络链路的毁伤，则某时刻后勤保障网络的状态
(1)
1.2 保障链路毁伤模型
令
表示某时刻节点v i 与v j 不直接连通或两节点之间的链路被摧毁；
表示某时刻节点v i 与v j 直接连通且链路完好；
表示某时刻节点v i 与v j 直接连通但链路被摧毁后未被完全修复。假定当
时，该链路暂时无法用于后勤保障。
假定后勤保障网络中未被完全摧毁的链路
在某一时刻是否被摧毁为独立事件
且
服从1重伯努利分布，记作
对于持续按数量m 去除链路的摧毁行为，每条链路被摧毁的概率
(2)
对于持续按比例α 去除链路的摧毁行为，链路被摧毁的概率
P n =α 。
(3)
综上，可将每轮打击对链路的摧毁行为描述为
(4)
式中：
和
分别为某时刻链路e ij 经历摧毁行为前后的状态。
1.2.4 加强安全管理。精神分裂症病患可能存在敏感、多疑、妄想、幻听等症状，为自己或他人造成伤害，故而在无陪护护理中，需对病患住院安全性作出评估，并增设对应的保障，如监控、门禁、巡视、呼叫等，尽量避免风险事故的发生。同时，加强病房巡视力度，预防一些睡眠少的病患睡醒后出现躁动；一旦出现躁动情绪，需作镇静处理，维持病房相对安静。
1.3 保障链路修复模型
假定被摧毁链路
必须依赖链路两端的完好节点才能完成修复，故可将被摧毁链路的修复行为描述为
(5)
(3)变异操作 采用单点变异策略对机器选择链进行变异操作，随机选择一道工序，改变该工序的当前加工机器。采用互换变异策略对工序顺序链进行变异操作，如图7所示，随机选择不同位置两道工序，互换其基因值。
(6)
其中，
为某时刻节点v k 正在修复的链路数量，故
意味着节点v k 在该时刻未对任何链路实施修复，可以被指派链路修复任务。
由式(6)可知：能够对被摧毁链路
实施修复的节点只可能是该链路的端点v i 或v j ，节点对被毁伤链路实施修复的前提是节点自身必须完好，且允许被指派链路修复任务。
为使模型更加贴近实际，本文特别提出2点假设：
建立由状态空间Y ∞和参数空间ψ 构成的坐标系，绘制出后勤保障网络毁伤过程终态链路变化的轨线图，从而考察不同修复速率对后勤保障网络毁伤过程的影响。
2.1 MAS 治疗4周时，3组患者的肌张力较治疗前明显改善(均P<0.05)，C组患者的肌张力改善较其余2组更显著(均P<0.05)，A组和B组差异无统计学意义。第8周随访时，A组患者肌张力较第4周及同时间点B、C组明显上升(均P<0.05)，B组与C组患者的肌张力较治疗4周时差异无统计学意义。见表2。
1) 每个完好节点同一时间只能对一条受损链路实施修复；
本文将每组仿真实验的后1 000次迭代结果的集合视为后勤保障网络毁伤过程的终态，即Y ∞={Y i |i >1 000}。若毁伤过程收敛，则Y ∞为包含一个或少数几个邻近元素的点集；若毁伤过程相对复杂，则Y ∞中将存在更多的元素，并且点集Y ∞往往具有较大的直径。不同打击强度下后勤保障网络毁伤过程终态轨线图如图2所示。
式中：
为链路e ij 下一时刻的状态；ψ 为被摧毁链路的修复速率
为节点集，包含某时刻后勤保障网络遭敌摧毁后，能够对被摧毁链路实施修复的全部节点，即
2 后勤保障网络毁伤过程的动力学仿真
针对图1所示4个网络算例，令ψ ∈{0.001x |0≤x ≤1 000,x ∈Z }，对于参数ψ 的每个取值，均以Y 0=M 、m =0.15M 为初值(不同网络初始链路数M 的取值详见表1)。根据式(2)-(6),将模型迭代2 000次，舍弃前1 000个迭代值，并将Y n 后续的迭代值全部绘制到参数ψ 对应取值的纵方向上，即可得到不同修复速率下后勤保障网络毁伤过程终态轨线图，如图3所示。
记录手术时间、术中出血量、并发症和置钉失败率，采用腰背痛视觉模拟评分（visual analogue scale,VAS）、Oswestry功能障碍指数（Oswestry disability in⁃dex,ODI）和日本骨科协会（Japanese Orthopaedic As⁃sociation,JOA）评分评定治疗效果。定期复查腰椎正侧位X线片、CT平扫并三维重建，评估椎间融合情况，观察螺钉有无松动、断裂，以及椎间融合器有无移位。
2.1 构造网络
后勤保障网络相对稀疏[23]，其中大部分节点的度相对较小，极少出现度很大的节点。由此，本文主要考察ER随机网络(Erdös-Rényi Network,ERN)[24]与WS小世界网络(Watts-Strogtz Network,WSN)[25]形式的后勤保障网络在毁伤过程中的动力学行为。为便于分析比较，本文针对ERN和WSN小世界网络分别构造2个平均度不同的全连通网络作为仿真算例。算例网络及度分布如图1所示，其基本参数如表1所示。
2.2 不同打击强度下的毁伤过程
建立由状态空间Y ∞和参数空间m 构成的坐标系，绘制出保障网络毁伤过程终态链路变化的轨线图，从而考察不同打击强度对油料保障网络毁伤过程的影响。
针对图1所示的4个网络算例，令m ∈{0.001x |0≤x ≤1 000,x ∈Z }，对于参数m 的每个取值，均以Y 0=M 、ψ =0.5为初值。根据式(2)-(6)将模型迭代2000次，舍弃前1000个迭代值，
图1 算例网络及其度分布
表1 算例网络基本参数
注：N 为网络规模；M 为初始时刻链路总数；
为平均度；
为平均路径长度；C 为聚类系数。
并将Y n 后续的迭代值全部绘制到参数m 对应取值的纵方向上，即可得到不同打击强度下后勤保障网络毁伤过程终态的变化情况。
2) 每条受损链路同一时间只能由一个完好节点实施修复。
2.3 不同修复速率下的毁伤过程
甲洛洛找话：怎么生这么多孩子啊？登子苦笑：我们也不想要啊，可我这个老婆，像头母猪，怎么都挡不住。甲洛洛：医院里不是有药买的吗？登子任然苦笑：填饱肚子都难，还买什么药啊。甲洛洛更加担心：你这么生下去，估计塔公村都装不下了？
限于篇幅，本文仅讨论基于持续按数量m 去除链路的毁伤过程。如前文所述，后勤保障网络毁伤过程主要取决于2个参数，即打击强度m 与修复速率ψ 。
图2 不同打击强度下后勤保障网络毁伤过程终态轨线图
图3 不同修复速率下后勤保障网络毁伤过程终态轨线图
需要说明的是，由式(2)可知针对后勤保障网络毁伤过程的模拟仿真存在一定的或然性，为尽可能减小极端情形对计算结果的影响，将上述实验重复10次，并将结果叠加绘制在其各自对应的轨线图上。
3 后勤保障网络毁伤过程的动力学分析
如图2所示，不同打击强度下后勤保障网络毁伤过程的轨线图相对简单。当打击强度较小时，后勤保障网络终态能够保持较高的完整度；随着打击强度的逐步增大，能够清晰地观察到轨线图中出现拐点；此后，轨线图纵向波动范围扩大，轨线总体呈现快速下降的趋势；当打击强度继续增大，轨线图开始趋于平缓，轨线纵向波动逐渐减小，最终收敛于下限Y min。根据1.3节提出的2点假设，后勤保障网络状态轨线图下限值取决于网络规模与相邻节点对受损链路的修复能力，即Y min=N ·ψ 。由表1可知：各算例网络规模均为100，链路修复速率均为0.5。因此图2中各后勤保障网络状态轨线图下限均为Y min=50。
4.3.1 观念僵化。在访谈到大部分毕业生的时候，对于创业的问题，他们的态度基本上是犹豫的，认为缺少资金、经验、人脉，并表示读书阶段创业不现实，也不知道从哪里着手，其实这也很大程度上反映出了大学生的创业意识薄弱。在另一方面来说，传统教育下，无论是学校还是家庭，对于在校期间创业都大体上视为是“不务正业”，甚至还会招到讽刺和批评，所以传统观念下的创业观念僵化死板。
将无线传感器网络应用于穿戴式枪声定位系统中，有效地解决了系统中单点麦克风阵列口径小，检测精度低的问题。多节点数据融合有效地改善了系统的检测结果，帮助边缘节点获取枪声声源位置。多次实验中，不同节点从感知枪声信号到接收融合结果之间时间不同，其最大值在2 s以内，完全符合枪声定位系统的实时性要求。
可见，当针对链路进行随机打击时，后勤保障网络在不同打击强度下的毁伤过程轨线图比较符合主观预期。这是因为链路的修复与链路自身无关，而是主要取决于受损链路的相邻接点。由于本文在模拟网络链路毁伤的过程中始终假定网络节点全部完好，故此时后勤保障网络毁伤过程的动力学特征并不复杂。
然而，不同修复速率下的后勤保障网络毁伤过程则呈现出完全不同的动力学特征，如图4所示。由图4可知：不同修复速率下毁伤过程的轨线不再平滑，轨线图中存在明显的阶跃行为。对于定义域内的某些特定取值(即阶跃点)，修复速率的微小变化将会对后勤保障网络毁伤过程产生较大影响。然而，当修复速率在相邻阶跃点之间移动时，后勤保障网络毁伤过程的总体变化却并不显著。
图4 不同修复速率下后勤保障网络毁伤过程的阶段性特征
进一步考察图4不难发现:各阶跃点的横坐标恰为0.5，0.33，0.25，0.2,…，即1/n (n >1)。随着n 的逐渐增大，阶跃行为的间隔按照1/[n ·(n -1)]的规律迅速收窄，毁伤过程的阶段性特征愈发模糊，在轨线图上则体现为修复速率越小，毁伤过程终态的轨线越平滑。
阶跃点横坐标表现出显著的幂次规律，意味着后勤保障网络毁伤过程的阶段性特征并非偶然。尽管从微观层面上看，针对链路的模拟打击具有较大的随机性，并且对链路的修复主要是基于节点主观的个体行为而非全局性的统筹协调。但是，包含了大量复杂非线性行为的毁伤过程却在宏观层面上产生了极为规律有序的动力学特征，即非线性动力学系统的涌现现象。
4 结论
本文通过大量模拟仿真，从理论上刻画了基于链路随机打击的后勤保障网络毁伤过程的完整轨迹，揭示和论述了毁伤过程中可能存在的动力学行为及其特征。提出了在不同修复速率下后勤保障网络的毁伤过程存在阶跃行为，并且阶跃点的分布具有显著的幂次规律，这对进一步加强我军后勤保障网络建设，提升战时后勤保障能力提供了理论依据和启示。不过，相关结论只是基于持续按数量去除链路的完全随机攻击模式得到的，并且未考虑网络节点同时遭受毁伤的情形。事实上，信息条件、打击策略、协同规则、节点与链路重要度差异等都将对后勤保障网络毁伤过程带来较大影响，这些也有待于进一步的深入研究。
参考文献：
[1] 尹峻松，王晓明，张新强，等.基于谱分析的复杂网络脆弱节点发现算法[J].无线电通信技术,2016,42(5):48-52.
[2] 周治彪，裘杭萍，杨志才，等.军事通信网节点重要度评估方法[J].指挥控制与仿真,2016,38(1):101-103,111.
[3] 吴辉，彭敏放，张海艳，等.基于复杂网络理论的配电网节点脆弱度评估[J].复杂系统与复杂性科学,2017,14(1):38-45.
[4] 崔文岩，孟相如，康巧燕，等.基于复合边权重的加权复杂网络级联抗毁性优化[J].系统工程与电子技术,2017,39(2):355-361.
[5] 张豫翔，吴明功，温祥西，等.边攻击下的航线网络抗毁性测度[J].电光与控制,2017,24(9):64-68.
[6] 张肖，戴明强，邵帅.基于边失效的舰艇装备保障网络抗毁性分析[J].舰船电子工程,2016,36(8):38-41.
[7] ARIANOS S,BOMPARD E,CARBONE A,et al.Power grid vulnerability:a complex network approach[J].Chaos,2009,19(1):13119.
[8] 程光权，陆永中，张明星，等.复杂网络节点重要度评估及网络脆弱性分析[J].国防科技大学学报,2017,39(1):120-127.
[9] 李晶，李宝德，王爽.基于突变理论的海上运输关键节点脆弱性度量[J].系统管理学报,2018,27(1):32-39.
[10] 李际超，吴俊，谭跃进，等.基于有向自然连通度的作战网络抗毁性研究[J].复杂系统与复杂性科学,2015,12(4):25-31.
[11] 李勇，邓宏钟，吴俊，等.不同流量的复杂保障网络抗毁性仿真分析[J].火力与指挥控制,2010,35(3):9-13.
[12] 谭跃进，吴俊，邓宏钟.复杂网络抗毁性研究进展[J].上海理工大学学报,2011,33(6):653-668.
[13] 李勇，邓宏钟，吴俊，等.基于级联失效的复杂保障网络抗毁性仿真分析[J].计算机应用研究,2008(11):3451-3454.
[14] 田田，吴俊，谭跃进.基于自然连通度的复杂网络抗毁性仿真优化研究[J].复杂系统与复杂性科学,2013,10(2):88-94.
[15] 吴俊，谭索怡，谭跃进，等.基于自然连通度的复杂网络抗毁性分析[J].复杂系统与复杂性科学,2014,11(1):77-86.
[16] 黄仁全，李为民，董雯，等.不同攻击策略下作战体系网络抗毁性研究[J].复杂系统与复杂性科学,2012,9(3):62-69.
[17] 彭兴钊，姚宏，张志浩，等.基于节点蓄意攻击的无标度网络级联抗毁性研究[J].系统工程与电子技术,2013,35(9):1974-1978.
[18] 段东立，吴俊，邓宏钟，等.基于可调负载重分配的复杂网络级联失效模型[J].系统工程理论与实践,2013,33(1):203-208.
[19] 田旭光，张成名.不完全信息条件下的装备保障网络抗毁性模型[J].系统工程理论与实践,2017,37(3):790-798.
[20] GAO Z,SMALL M,KURTHS J.Complex network analysis of time series[J].EPL (Europhysics Letters),2016,116(5):50001.
[21] SOL
CONDE L E,USED J,ROMANCE M.Optimal distributions for multiplex logistic networks[J].Chaos:an interdisciplinary journal of nonlinear science,2016,26(6):65312.
[22] DONNER R V,HEITZIG J,DONGES J F,et al.The geometry of chaotic dynamics:a complex network perspective[J].The European physical journal B,2011,84(4):653-672.
[23] 李勇.物流保障网络级联失效抗毁性研究[M].武汉:武汉大学出版社,2014:24-30.
[24] ERDÖS P,RÉNYI A.On random graphs.I.[J].Publicationes mathematicae,1959,4:3286-3291.
[25] WATTS D J,STROGATZ S H.Collective dynamics of ‘small-world’ networks[J].Nature,1998,393:440-442.
Research on Damage Process of Logistics Support Networks Based on Link Random Attacks
DENG Qun
(POL Department, Army Logistical Academy, Chongqing 401331, China)
Abstract ：Taking Erdös-Rényi(ER) network and Watts-Strogtz(WS) network as examples, the damage process of logistic support networks is simulated and discussed under link random attacks. By analyzing the simulation results, it is found that the damage process of logistic support networks has obvious phased characteristics, the damage process under different repair rates exists step phenomenon, and the distribution of step points shows a remarkable power law. This paper attempts to explore the dynamic characteristics contained in the damage process of logistic support networks from a new perspective, which has certain practical significance for improving survivability of logistic support networks on the battlefield.
Keywords : logistic support; complex network; system dynamics; damage process; invulnerability
文章编号： 1672-1497(2019)03-0078-06
中图分类号： E92
文献标志码： A
DOI： 10.3969/j.issn.1672-1497.2019.03.014
收稿日期： 2019-04-16
基金项目： 全军军事类研究生课题
作者简介： 邓 群(1991-)，男，博士研究生。
(责任编辑：牛燕平)