基于蒙特卡罗方法线路雷击跳闸仿真.docVIP

基于蒙特卡罗方法线路雷击跳闸仿真 　　【摘 要】随机抽样雷电流幅值和陡度、导线工频瞬时电压，用较符合雷击实际过程的蒙特卡罗法计算线路雷击跳闸率。模拟结果显示考虑工频电压叠加对绕击闪络次数影响不大，但对反击闪络次数有一定影响。小幅值大陡度的雷电流也可以引起绝缘闪络，而规程法计算结果不会引起绝缘闪络，规程法计算雷击跳闸率的结果偏向于保守。对雷电参数进行大范围、充分的普查，才能提高雷电参数概率分布函数的可信性和合理性，也才能更精确的进行雷击跳闸率计算。 　　【关键词】雷击 跳闸 　　规程法耐雷水平计算雷电流波形的波头时间统一为2.6μs，这样的计算没有准确反应雷电活动的客观规律。实际上每次雷击参数都不一样，有些小幅值、大陡度的雷电流也能引起线路跳闸。雷击输电线路瞬间，导线上工频电压瞬时值是随机量，虽然瞬时值对220kV及以下线路的影响不大，但对较高电压等级的线路却有着较大影响，直接关系到了绝缘子串两端的闪络电压取值。传统计算方法一般不考虑工频电压的影响，或者采用半周工作电压的平均值作为修正参数。计算输电线路雷击跳闸率应该通过一个统计的方法，才更为科学合理，更符合实际运行情况。蒙特卡罗算法正是满足这些要求的一种有效的统计模拟算法，通过随机抽样产生满足各种分布的抽样序列，在计算机上进行统计模拟试验，然后给出问题数值解的统计估计值。 　　1 雷电流分布规律 　　根据浙江220kV新安江―杭州输电线路Ⅰ回路在1962～1988年27年间用磁钢棒实测到的104个雷击输电线路塔顶数据，我国电力行业规程推荐雷电流幅值概率分布表达式为 　　（1） 　　2005年IEEE工作组对全球雷电研究工作进行总结和回顾时，推荐Anderson提出的雷电流幅值概率分布的近似对数正态分布表达式[6]。假定雷击参数x的统计变化遵从自然对数正态分布，则其概率密度函数为 　　（2） 　　式中：σlnx为lnx的标准偏差；xm为x的中值。 　　当两个雷击参数符合二维正态分布时，其联合概率密度函数为 　　（3） 　　式中 　　， ， 　　，ρc为相关系数。 　　雷电流I遵从自然对数正态分布，即X=lnI（-∞X∞）遵从分布参数为μX=lnIm、σX=σlnI的正态分布，Im为雷电流的中值，则I的概率密度函数、分布函数分别为 　　（4） 　　（5） 　　2 算法流程 　　蒙特卡罗法计算线路雷击跳闸率步骤如下： 　　（1）输入输电线路的几何数据和雷电参数的初始数据。计算避雷线对外侧导线的耦合系数，避雷线、导线的平均高度、杆塔电感、线路绝缘水平、线路绕击率、建弧率、雷电流幅值和陡度的分布等，相关参数可以参照规程法的计算。 　　（2）控制模拟线路遭受雷击的次数。生成的雷击次数太少容易导致结果精确度降低，因此在计算时可以放大样本的容量。 　　（3）随机抽样雷电流幅值和陡度、导线工频瞬时电压。规程法实际上是把雷电流幅值和陡度看成相互关联的雷电参数，如果认为雷电流幅值和陡度不相关，而是互相独立，则可以分别生成两组独立的随机抽样序列。 　　（4）生成雷击输电线路部位的随机数。当随机数小于导线绕击率时则看做发生绕击，当随机数大于导线绕击率而小于击杆率时则看做雷击塔顶，当随机数大于击杆率时则看做雷击档距中央，本文不考虑该种情况。 　　（5）分别计算绕击和雷击塔顶的闪络次数。绝缘闪络的判据以绝缘子串的50%冲击放电电压为标准，当计算出的线路绝缘子串两端电压大于50%冲击放电电压时则定义为发生一次雷击闪络。据试验统计分析，绝缘子串的50%冲击放电电压是遵从正态分布的随机量，本文暂时以常数处理，这会带来一定的误差。 　　（6）计算输电线路的雷击跳闸率。按每年40个雷暴日的雷电活动频率，线路长度为100km，折算雷击跳闸次数。 　　3 算例分析 　　算例采用220kV线路，全长100km，架设在平原，酒杯型铁塔，绝缘子链由13片LXP-70型绝缘子组成，绝缘子串长度为2.2m，避雷线弧垂7m、半径5.5mm。避雷线间距11.6m，与外侧导线间距1.7m，距横担3.5m。导线弧垂12m，距离地面23.4m。保护角α=16.5o，杆塔冲击接地电阻以7Ω计。 　　避雷线、导线的平均高度为：hb=（23.4+2.2+3.5）-2×7/3≈24.5（m），hd=23.4-2×12/3=15.4（m）。计算冲击放电电压：U50%=100+13×84.5=1199（kV），取击杆率：g=1/6，绕击率： 　　。 　　雷击塔顶时的耐雷水平： 　　绕击导线时的耐雷水平：I2=1199/100=11.99（kA），由此计算出跳闸率：n=0.249 （次/100km?年） 　　用蒙特卡罗法计算跳闸率，假定雷电流幅值和陡度服从自然对数正态分布：Im=34kA，