连锁故障的概率性负载依赖模型和停电的可能影响

1、连锁故障的概率性负载依赖模型和停电的可能影响摘要大型、互联的基础设施系统的灾难性的破坏往往是由于连锁故障。例如，电力系统的大型停电事故通常是由负载过重造成的系统连锁故障。我们引入一个连锁故障的连锁模型，该模型系统是一个具有许多随机加载的相同组成部件的系统。初始扰动会导致某些组成部件超出其承载极限，故障部件会导致其他部件的固定负载增加；而随着故障部件增加，系统更进一步发生连锁故障成为可能。故障部件数量的概率分布是一个扩展的准二项分布。使用递归推导扩展准二项分布的显式公式。在一个受限制的参数范围内的连锁模型提出了服从准二项分布的新模型。一些性质特点由扩展准二项分布表明，包括电源的末端，总指数末端和

2、整个系统故障的显著概率。1、连锁模型连锁故障是大，互联基础设施系统中的一个普遍的灾难性故障原因。例如，负荷依赖性连锁故障发生在电网大停电时。这些基础设施对社会的重要性，促使建设和研究模型来捕捉连锁故障的显着特点。我们定义概率连锁故障的连锁模型具有以下的一般特征：1、有多个相同的组成部件，每一个都具有一个随机的初始负载和初始扰动。2、当一个部件过载，它不能将负载转移到其他部件。特征2会导致连锁故障：其他部件产生另外一个的故障造成一些其他部件也可能故障，导致一个连锁故障。部件因为连锁效应逐步地过载。连锁模型的初始版本是用来研究电力传输系统的关键的负荷变化和功率末端停电大小的概率分布 8 。1.1、

3、对模型的描述连锁模型具有N个相同随机初始荷载部件。对于每一个部件的最低初始负载是L，最大初始负载是Lmax。J = 1,2，N，J部件具有初始载荷L J，是一个随机变量（L min，L max 。L 1，L 2，···，L n是独立均匀分布。当负荷超过L时，部件就会发生故障。当一个组件发生故障时，一定量的负荷P传递到其他每个组件开始连锁反应，我们假设的初始扰动加载到其他每个组件的额外值为D.其他组件是否发生故障取决于它们的初始载荷L J，并且这些组件任意一个发生故障时将分发的附加载荷P0，可能会导致进一步的连锁故障。连锁模型可以更精确地定义在算法的形式：连锁模型算

4、法：0、所有n个部件都开始有允许初始荷载L 1，L 2，L N···，确定为独立随机变量均匀分布在Lmin，Lmax。1、将初始扰动D添加到j=1,2，n组件的负载。2、测试正常组件到发生故障：如果组件j是正常的，当负载L发生故障，假定组件在Mi这一水平发生故障。3、如果Mi=0，停止测试；连锁过程结束。4、如果Mi> 0，然后根据故障的数量Mi增加的组件负载：增加Mi功率到组件j，j=1，5、增量迭代计数i，转到步骤2。一个简单的连锁模型的例子有5组分产生连锁列于表1。在这个连锁中，2个组件在迭代1次故障，一个组件在迭代2次故障和一个组件在迭代3次故障，然

5、后连锁结束，共有4个组件故障。也就是说，m1 = 2，m2 = 1，m3 = 1，然后连锁结束，共有m1+m2+m3= 4个组成部分故障。图1显示了这一系列负载在这一连锁反应中参照迭代数的增加量。表1、一个5组分连锁反应中载荷的增加实例：由于连锁故障是由的组件荷载与L的关系唯一确定，一个等效的方式来可视化连锁显示了组件负载好像固定在它们的初始值并且L的水平线向下变化，如图2所示。在这一点上，该行在第一次向下移动是因为初始扰动D，然后取决于随后的迭代Mi P。在图2中，Mi 是一个组件故障的前一个数，很方便地显示为0，该线路略过。图1、5个组件构成的级联实例中的组件负载增加。0表示初始荷载和1,

6、2,3,4表示迭代1,2,3,4分别加载。当其负载超过一个水平线L故障时，组件发生故障。图2。另一种可视化的连锁例子与负载（零），好像固定在其初始值并有L故障的水平线向下移动。该行第一次移动下降了D=3，随后按比例转移到以前的故障数（0的数目的零线刚刚通过）。考虑连锁模型的算法中的连锁控制中的组件故障的事件。为了方便描述，我们按故障发生顺序重新编号。对事件的准确描述相当于是，R分量不存在正整数k和M 1，M 2，M···，K，4页上的方程（16）。我们举一些例子来解释（2 - 6）。由于组件1在第一次迭代后故障，此时其负载L1因为初始扰动D增加， L 1

7、 + D > L 故障。这是例子解释（2）。第一次迭代后，所有的负载都增加了M 1 P.因为M 1 + 1组件第二次迭代时故障，L M 1 + 1  +M1 P + D > L故障。此外，由于M 1 + 1组件不能在第一次迭代故障，L M 1 + 1 + DL故障。这是例子解释（3）。组件r+ 1永不故障，并且组件r+1的负载因D+rP增加在连锁反应的末期。因此L r + 1 + D + rPL故障。这是例子解释（6）。在表2中总结了连锁模型。1.2标准化模型假设初始扰动D改变为D + K，故障载荷L改变为L + K 。K添加到初始扰动增加了所有荷载的在开始时的连锁效应，

8、但由于故障负载也因为K增加了，组件的故障是不变的，并将在整个连锁过程中持续。也就是说，组件故障的连锁是不变的，所有的负载和故障负载只是因为K增加数值，此结论在K的正反两面都成立，因此我们有：定理一、如果初始扰动和组件的故障载荷L增加相同数量的K，连锁过程不变。原则1也遵循从表达式（2 - 6），仅仅取决于不同LD。模型规范在两个步骤是很方便的：第一步是改变初始扰动D为D + L max L fail并改变故障载荷L fail 为L fail + L max L fail = L max。根据定理1 和K = L max L fail，这些变化对连锁过程中的组件故障没有影响。其次定义了标准化的初

9、始负载：其中lj是一个随机变量在0,1均匀分布。此外，故障负载lj= 1。令P是任意组件的负载增长量，当另一组件发生故障后它表现为一个负载，其值在Lmax到Lmin之间。类似地，D是初始扰动表示一小部分的负载波动。在表3中总结了连锁模型及其标准化参数。此外，标准化的条件下的组件的故障显示在（1）和（7 - 11）中。2 故障组件数分布本节导出了故障组件数概率分布的递归和显式公式。定义：F（r、d、p、n）是r组件在在连锁模型中发生故障的概率，模型是标准化初始扰动d、标准化的负载转移量p、以及n个组件。2.1 当d 0和 d 1时的情况当初始扰动 d 0时，没有组件发生故障，并且当 d 1时,

10、所有的n个组件立刻发生故障，并且2.2 从0 < d < 1导出递归本段假定在0d1。初始扰动会d导致初始负载的组件立即失效，初始负载在（1d，1 。因此，任何组件立即故障的概率是d，所以任何组件没有立即失败的概率是1d。由于初始荷载独立，R = 0部件故障的概率是：此外，在一个组件的情况下（n1）其他部分，我们假设，1rn；如果r1组件故障，那么一定数量组件k，k有 1kr必须因初始扰动d立即故障。Ek表示该事件，k为立即故障的组件。初始组件的负载是独立的。由于E 1，E2，.，E r为互斥和完备集合的事件，全概率的公式给出：我们主张：建立公式（20），考虑到nK组件在E K条件

11、下没有立即发生故障。由于n-k组件没有立即故障，他们的负荷l必须在0，d）并均匀分布在0，d （即，条件Ek的均匀分布在0，d 。），对n-k组件而言，L min = 0和l max = 1d。在Ek，每个nk组件有一个负荷从初始扰动增加到d，一个额外的负载从k组件的立即故障后增加Kp。因此，Ek后，每个n-k组件总的初始扰动是d = kp + d。总结一下，Ek后，nk最初未故障的组件是由连续模型的初始扰动d = kp + d决定，负载转移P = p，L min = 0，L max = 1d，L fail= 1，和n-k组件。使用（13）标准化的收益率n-k最初未故障的组件是由具有归一化初始

12、扰动Kp 1d的连续模型，归一化负载转移1d和n-k组件决定。因此，r组件故障率由Ek表示，n-k个组件中r-k的故障率由式（20）给出。结合（18），（19），（20）产生递归：2.3 扩展的准二项分布公式f本节建立扩展准二项分布公式为故障的组件数量分布f。以下对准二项分布的特殊情况进行了讨论。我们总结的公式（15 - 17）和递归（21）：方程（2225）和（26）的递归定义f（R，D，P，N）对所有n1和d > 0成立。这是一个简单的证明，通过归纳证明f（r、d、p、n）是一个对于所有的n1的概率分布，在附录A中详细介绍考虑表达式：附录B证明（27-29）满足（2226），因此对于

13、扩展准二项分布概率分布的显式公式。所以由np + d1给出的概率分布（27-29）化简为（27）：这是准二项分布引进康瑟尔 5 模型anurn问题，模型由玩家的战略决策。因此（27-29）是一个扩展准二项分布，参数范围的扩展到让np + d > 1。我们下面将看到这个扩展的参数范围所有组件故障的概率很高。康瑟尔 5 模型的准二项分布的平均值（30）为：3 连锁分布实例本节展示了组件n = 1000的梯级分布的定性行为建模实例和如何选择标准化的连锁模型的参数。3.1故障平均次数一种方法来概括的连锁分布的参数变化的行是检查的故障组件的平均数。图3所示为参数d和p的函数曲线图。在n = 100

14、0情况下考虑图3，曲线np + d = 1相连于（p，d）=（1 / n，0）=（0.001,0）（p，d）=（0,1）。图3的区域在曲线np +d= 1右侧并且曲线d = 0以上为黑色部分是由于高概率造成的，此时所有组件故障的平均数接近1000（但注意的角落附近（P，D）= 0.001,0的黑色区域是圆形，使附近的（P，D）=（0.001,0）是白色的）。该地区的曲线np + d = 1的左边和曲线d = 0上边，准二项分布（30）适用，故障平均数< r >是由（31）决定。特别是，< r >对d是成比例的.下面的小节检查扩展准二项分布更详细地在图3中显示了常数p和d

15、的不同。3.2 p=0的二项分布连锁模型的在p=0的限制的情况下，分布成为二项分布：众所周知，对于大n和小d，（32）是很好的近似高斯分布的平均nd和方差nd（1d）。3.3 常数小d的性质我们假设的归一化初始扰动是固定小的值d= 0.0001，并检查如何分布变化随着p的增加。对于p= 0，分布是二项分布像以上提到的。如图4所示，从零开始的分布。p=0.0001的分布末端不同于二项分布。但仍然近似指数分布。随着p的增加尾部变得较重，并且p的分布在某范围内的是一个近似的幂尾的分布，p=0.001近似满足条件的np + d1= 1。对p=0.002的分布有一个小r的约指数尾部，中间r的概率为零，并

16、且1000个组成部分故障的概率为0.08。（如果一个中间数的组件故障，则连锁总的1000个组件故障。）以前，限制版本2的连锁定性获得了类似的结果，通过假设P =d，增加p.对这些结果进行了类似定性仿真得到的结果通过增加负荷停电模型引起了电力传输系统的连锁故障 8 。提高标准化功率传输标准被认为是可以加强交互连锁故障的组件。在pdf的进展，指数尾随着p增加。当电力系统发生故障时，它们对连锁故障发生风险有重要的影响。幂尾指数可以近似如下。假设np + d1。通过写作= np和=nd，让n0 ，p和d0，在这样一种方式下= np和= nd固定，它可以在 6 利用Stirling公式，准二项分布分布f

17、（r，d，p，n）可以被近似为：这是广义泊松分布 6 。当r >> 1, 使用Stirling公式得到：在极限条件下，np + d = 1成为=1，设置了=1得到：当此时有。3.4 最大限度减少负荷的影响在级联模型中，试图减轻连锁故障的一种方法是减少最大的组件负载量，使Lmax低于Lfail。让K = LfailL max是最大负荷量的减少量。此时，该系统在任何初始扰动D小于K时都不会故障，我们将观察k从零开始增加造成的的影响。标准化初始扰动和归一化功率转移（k）是k的函数，特别是根据（13），因此并且组件r故障的概率是：这说明了减少最大负载的影响。如果kd（0）或d（k）0，将不

18、存在故障情况（14）。对于小k，d（k）d（0）k（1d（0），使得归一化的初始扰动减小。然而，标准化的负载转移量。一个例子如图5所示。没有故障的概率（r= 0，而不是如图5所示）是0.37为0和0.90，k= 0.0009。在这种情况下，无故障的概率缓慢增加，并减少故障的概率在在r1的范围内的。然而，增加在P（k）使末端分布较重。3.5 平均负荷变化的影响我们假设初始扰动的初始值为0，并通过改变初始载荷的平均值，考察初始载荷变化的影响A。初始载荷W =LmaxL min范围。然后，初始载荷范围可以被参数化为L和W：我们选择让W是恒定的，检查不同L的效果。于是，根据（13），归一化初始扰动D（

19、L）成为一个仿射函数L,并且传递功率P是一个常数：组件r故障的概率是：L<LfailW / 2，如果d（L）0，将不会发生故障（14）。如果L>Lfail+ W / 2，则d > 1并且所有组件失败（15）。现在考虑LfailW / 2LLfail+ W / 2，0D（L）1。特别是，我们选择Lfail= 1和W = 0.2，0.9L1.1，选择P = 0.00075和n = 1000。图6显示为d = 0.0005，0.05，0.2，0.25，分别的PDF，L = 0.9001，0.91，0.94，0.95。在这种情况下，np + d = 1发生在d = 0.25，对于d

20、> 0.25概率下pdf迅速造成所有组件故障。在范围0d0.25内的准二项分布公式（30）和（31）适用，因此故障的组件随d和L线性增加。4停电应用连锁模型是非常简单的代表了现实的一个大型电力系统的各个详细方面。因此，连锁模型的一般特征可以存在于大停电的传输线和发电机是合理的。电力系统组件故障的交互方式有很多，其中包括通过保护系统、功率流的重新分配、动态操作或规划错误。然而，按照连锁模型，这些相互作用是负载的依赖性，一个组件的故障倾向于强调其他组件，使更可能发生故障。如果我们专注于由于线路超载和和中断通过再分配的功率流（直流潮流与LP调度）的相互作用产生的连锁故障，然后结果的OPA模型1

21、，2，7 表明，各种运行机制是可能的 1，2，4 。这些机制下大面积停电的分布的例子在连锁模型中观察到的类似于3节。3节中的例子也说明了不同的模型和参数化可以在连锁模型定义并进行归一化处理通过梯级分布公式的应用。例如在3.4节的连锁故障大致模型中运用n-k安全准的整体特征使得最大负荷降低。（如果我们把初始组件故障认为是给n个剩余的组件提供的初始扰动，然后在最大负荷的减少大致接近的连锁效应的n-k准则中，系统可能存在一定数k的初始组件故障比例K.分析第3.4认为只有的n剩余的部分，此处n=n-k。该连锁模型可用于想法的测试，如连锁故障在复杂情况的大型电力系统大停电。事实上，在第3.3节中，一个限制版本的连锁模型已经被用来显示电源的末端行为定性，并与在大停电模型 8 中获得的相似。利用连锁故障模型更好地了解连锁故障在停电的OPA模型 1，2，4 的进一步应用，电力系统本身是有前途的工作。我们已经开始探索的连锁故障模型与光纤束材料破坏模型的的关系。而不是关于初始负载Lj作