元件失效模型0609217分析.docVIP

0609217.pdf 1. 随机线路故障模型 线路l上发生的随机故障由泊松过程描述。即在时间段t中发生的故障的次数服从泊松分布，均值为。因此线路l上发生的事故次数具有下述概率分布： （式1） 期望值为 （式2） 仿真时通过产生随机的故障时间t，从而确定式（2）中的t。 两个相邻时间t1, t2之间的时间间隔符合指数分布： （式3） 即：t0是起始点，，通过对式（3）抽样获得。 2. 线路过负荷故障模型 对输电线由于输送的潮流超过其传输容量极限而故障这一过程建模。即描述的是在不同线路潮流大小及不同气候条件下，经过多长时间，导体温度达到临界温度，发生故障。 主要考虑导体过负荷发热问题。 假设，线路表面的热量正比于导体与周围空气温度之差。，T是导体的温度，T0是周围环境的温度。H是表面的热传导系数。这样，热传导问题就变成了一个由时间和和位置（即距离导体端点的距离）决定的线性的热流问题。 假设线路的初始温度是T0，通过线路的潮流是常数P，线路温度的演化可以由下式描述： 其中，（式1），V是电压，，p为半径，。。。 当线路上的潮流大于其上限一定时间，导线就会由于发热而下垂，发生故障。 Tcl为导线通过最大允许潮流的温度（由式1计算），则故障时间（也就是当导线过负荷时，多长时间导线故障）则由下式计算： 其中，tcl为故障时间，T(0)为导线的初始温度。 tcl即为导体温度达到临界温度的时间。 3．线路修复模型 其实就是如何对线路经过多长时间之后恢复运行进行建模，也就是线路修复时间 其中，c为常量，为最小修复时间，w为随机变量，具有指数分布。 w由下述指数分布抽样获得： 假设所有线路的相同。 4. 调度员的响应 （1）以概率p1断开线路，防止其故障及被损害； （2）以概率p2，不进行任何操作，任由线路由于过负荷达到其临界温度； （3）以概率p3，调整发电机出力或部分甩负荷，消除线路过负荷。 P1+p2+p3=1； 上述不同策略会对线路再次投入运行的时间产生影响。响应（1），线路完好，可以立刻投入运行；响应（2），线路被损害，则需要经过一段时间后，修复后，投入运行。 线路再次投入运行时，进行潮流计算，根据计算结果是否存在过负荷，调度员可以采取下述两种策略： （4）以概率p4，将线路再次投入运行，即使存在过负荷； （5）以概率p5，再次调整发电机出力及甩负荷，消除过负荷。 P4+p5=1； 5. 仿真过程： 问题：一共有多少种不同的事件？ Type0，线路恢复； Type1，调度员在线路存在过负荷后断开线路； Type2，调度员不作任何调控，任由线路过负荷存在； Type3，线路过负荷一段时间后，线路温度达到临界温度，故障。 另外：存在负荷波动，每一次负荷波动之后，进行潮流计算，计算结果会显示是否有过负荷。 （1）抽样，决定每一条线路发生随机故障的时间；ti （2）当任何一条线路随机故障后，则会经过一个随机时间tr恢复运行（type0）；ti+tr后，则会再次为这条线路产生一个随机故障的时间。 （3）如果有过负荷，则，调度员采取可能策略之一（根据产生的随机数及概率判断），例如type2事件，则经过一段时间tc之后，线路故障断开，然后又是经过一段随机修复时间后恢复运行。。。。如何调度员断开线路，则可以在随后的任何时间使线路恢复运行，恢复运行后，进行潮流计算，则调度员可能采取（4）或（5）中之一。如果（4），则再次经过一段时间后，线路过负荷温度过高故障。。。 clear global genchsh %定义全局标量发电机的出力 global xdpyuan %系统中的负荷节点容量 [baseMVA, bus, gen, branch, areas, gencost] = case30;% in put data gon=find(gen(:,8)0); % 找出运行的发电机 bon=find(branch(:,11)0); % 找出运行的支路 gen=gen(gon,:); % 形成运行的发电机数据 剔除未运行的发电机 branch=branch(bon,:); % 形成运行支路的数据 剔除未运行的支路 [i2e, bus, gen, branch, areas]=ext2int(bus, gen, branch, areas); % 重新形成 包括运行发电机和运行支路的算例数据 nb=size(bus,1); % jie dian shu nd1=find(bus(:,3)0);