atp报告-500kV变电站中雷电暂态过程的仿真计算.docx

电力系统过电压计算实验报告500kV变电站中雷电暂态过程的仿真计算1研究背景与意义在雷电活动频繁的地区，由于雷击造成变电站内电气设备损坏的情况时常发生。变电站内设备因雷击造成的绝缘损坏大多由输电线侵入的雷电过电压波引起。为保证电网的安全稳定运行，减少变电站电气设备雷击损坏，应对变电站的雷电侵入波特性进行研究分析，以确保在各种运行方式下各电气设备因雷电引起的过电压值低于其额定雷电冲击耐受水平（BIL）。研究沿线路传来的雷电波侵入变电站的波过程时，由于输电线的分布参数复杂，且变电站内设备和回路分支众多，雷电侵入波在变电站内将发生复杂的折反射。变电站的雷害可能来自两个方面：一是雷直击于变电站；二是雷击输电线路产生的雷电过电压波沿线路侵入变电站。对直击雷的防护一般采用避雷针或避雷线。我国运行经验表明，凡按规程标准要求正确安装避雷针、避雷线和接地装置的变电站，绕击和反击的事故率都很低，防雷效果很可靠。另一方面，线路落雷比较频繁，虽然线路入侵变电站的雷电波电压受到线路绝缘的限制，但变电站电气设备的绝缘水平较线路低。同时，变电站虽设有进线段保护，然其与非进线段并无本质差异，同样可能遭受雷击而形成侵入波。由于变电站内有许多联络短线（如变压器到母线和避雷器的连线等），它们和输电线路一样，在持续时间极短的雷电波作用下，表现为各个联络线段间行波的传播和快速折、反射过程，并且通常在此过程中产生瞬时幅值极高的过电压，容易对设备造成危害。因此若防护措施不可靠，势必造成变电站电气设备的损坏事故。由此可见，沿线路侵入变电站的雷电过电压波是对变电站电气设备构成威胁的主要原因。因此，对雷电侵入波在变电站电气设备上所产生的过电压进行仿真计算，找出过电压的分布变化规律，能够为限制雷电过电压（以保护变电站电气设备的绝缘）提供有价值的参考依据，从而合理配置避雷器保护，以进一步优化变电站的工程设计。为此，对某500kV变电站进行仿真计算，分析雷电侵入变电站时，在不同雷击点与变电站的距离、不同的运行方式、不同的杆塔冲击接地电阻、不同的等效入口电容、避雷器的配置以及其他因素下，雷电过电压对设备的影响。图1-1为500kV变电站电气主接线图。图1 500kV变电站单线图2 系统建模及参数设置2.1 雷电流模型雷电的物理过程虽然复杂，但是从地面感受的实际效果和防雷保护的工程实用角度，可以把它简化为一个沿着一条固定波阻抗的雷电通道向地面传播的电磁波过程，据此建立计算模型如图2.7。在雷电放电过程中，人们能够测知的电气参数主要是雷击地面时流过被击物体的电流i，然后再根据计算模型反推出雷电波的参数以供工程应用。由图2.7（a）电流源等值电路，有公式如下：其中：Z为被击中物体的阻抗，Z0为雷电通道波阻抗，I0为雷电流波幅值，A显然i与Z有关。当Z=Z0时，恰好i=i0；当Z=0时，i=2i0，而Z不可能为零；但若ZZ0，仍可测得i2i0。在雷电流的实际测量中一般都能满足条件ZZ0）物体时，流过该物体的电流定义为雷电流。定义中的雷电流i恰好等于沿雷电通道传播而来的雷电流波的两倍。因而在防雷保护计算的彼得逊等值电路中，如图2.7（a）所示，等值雷电流源通常就用雷电流来表示。根据国内外的实测统计，75%90%的雷电流是负极性的。因此电力系统的防雷保护和绝缘配合通常都取负极性的雷电冲击波进行研究分析。电力系统的防雷计算、保护设计都要求将雷电流波形等值为典型化的可以解析表达的波形。常用等值波形有三种：标准冲击波、等值斜角波和等值余弦波。本文的仿真中，雷电流波形选用2.6/50ms的等值斜角波。前苏联科学家通过观测和计算得出雷电通道波阻抗Z0在3003000之间，如图2.7（b）所示。I100kA时，Z0稳定于300附近，在反击计算中雷电流幅值很高，Z0较小可以取为300W。（a）雷电流源等值电路 （b）Z0和I关系曲线图图2.7 雷电流源等值电路和雷电通道波阻抗现有的雷电流幅值概率分布是由一般高度的杆塔上测量数据（基本上是下行雷）得出的，不包括上行雷。上行雷是由高建筑物引发的特殊放电形式。防雷电反击计算应该把上行雷和下行雷分开计算，否则会带来较大的误差。超过100m的高杆塔，随着杆塔高度的增加，落雷次数增加，上行雷发生的概率也越来越大。3条线路杆塔高度不超过70m。根据上行雷发生概率和建筑物高度的关系图（图2.8），上行雷的概率不到10。因此在本研究中，计算反击闪络率时只考虑下行雷的落雷次数，而把上行雷部分除外。图2.8 上行雷发生概率和建筑物高度关系图对于本次计算的500kV系统线路，其反击耐雷水平均大于100kA，所以从图2.7(b)中看出，当I100kA时，Z0稳定于300附近，在反击计算中雷电流幅值很高，Z0较小可以取为300W，而雷电流波形选用2.6/50ms的等值斜角波。图2.9 雷电流源的仿真模型与数据仿真计算使用的雷电流在2.6us处到达峰值，波尾时间为50us，在雷电流源侧设置检测模块观察雷电流波形，以峰值为50kA的雷电流为例，雷电流波形如图2.10所示。图2.10 雷电流波形图（50kA为例）2.2 杆塔模型近年来，国外一些专家尝试通过试验的方法来建立新的杆塔模型，并且建立了多种多波阻抗模型。多波阻抗模型是基于垂直导体不同高度处的波阻抗是不同的原理，使用于尺寸较大的杆塔。如图2.4所示，在t=t1时，行进波到达半径x1=vct1覆盖的区域内，此时该波不能到达其他区域，因此距起点为x1的这段垂直导体的波阻抗Z(x1)是可以通过x1区域内的几何和物理参数得出的，由于几何参数x1不同于t=t2时的x2，所以Z(x1)不同于Z(x2)。因此，距起点不同的地方，垂直导体的波阻抗是不同的，所以就可以将垂直导体分割成几段，每一段都可以计算出一个波阻抗。图2.4 垂直导体上的波行进由波阻抗的特性可知，垂直圆柱体的波阻抗值仅依赖于该圆柱体的半径和对地高度，由此我们可以用下面的经验公式来描述单根垂直导体的波阻抗。式中：r 是垂直圆柱体的半径/m；h 是垂直圆柱体的高度/m。以上式为基础，我们继续研究多导体系统的波阻抗。现有一个多导体系统，包含数根平行的圆柱体，则这n根平行圆柱体组成的系统平均波阻抗为：式中：n 为圆柱体的数目；ZT,kk 为第k根圆柱体的自波阻抗/W；ZT,kl 为第k根和第l根之间的互波阻抗/W。假定上式对于该系统是有效的，那么ZT,kk和ZT,kl可表示为：式中：Rkl 是第k根和第l根间的距离/m。由于杆塔最多是4导体系统，所以给出了2、3、4多导体系统的波阻抗计算公式。式中的re为多导体系统的等效半径，可以由下式算出： 式中：R 是临近导体间的距离，如图2.5(a)所示/m。 对于非平衡多导体系统，如图2.5(b)等效半径公式是基于测量结果的经验公式： ，(a)平行导体组成的多导体系统 (b)包含不平行导体的多导体系统图2.5 多导体系统Hara T等学者认为，杆塔可以等效为图2.6所示的等值分布多波阻模型。杆塔的主体可以分为四段，每段杆塔可以用波阻抗分别为ZT1-ZT4，波速都为光速的分布参数表示。图2.6杆塔的多波阻抗模型主支架部分每部分的波阻抗ZTk计算公式为： （k=1 4）式中rek为多导体系统的等效半径，可以由式算出：如上式所示，由于本文采用的杆塔塔基半径比通常的要大，所以在计算re1re3时用来替代，而在计算re4时才采用。实验表明，测得在有无支架的情况下的波阻抗，从测量结果中可看出：导体系统的波阻抗在增加了支架之后减小了10%左右，从而每段杆塔对应支架部分的波阻抗ZLk计算公式为： （k=14）同时，电磁波通过含有支架的多导体系统时需要更长时间，在仿真中，模型支架部分的长度可取为对应主体部分的1.5倍。相关文献表明，可以用分布参数来模拟杆塔臂。取其波阻抗为ZAk，波速为光速，长度可以取它自身的长度。波阻抗ZAk可以由下式得到： （k=14）式中rAk为第k部分横担的等值半径，可取为杆塔臂与主杆塔连接处长度的1/4。多波阻抗计算模型结果如下：2.3 线路模型对变电站内的线路采用连续换位的Clarke模型。具体参数如下：对变电站外输电线采用8线的JMarti LCC 模型(相导线加架空地线)，即6条相导线和2条地线（含一条避雷线）的模型。因这种架空线的模型直接计算了地线与导线之间的耦合关系，所以在计算过程中就不必再考虑被击避雷线与已闪络导线对未闪络导线的耦合电压，这样简化了计算，提高了计算精度。2.4 绝缘子串模型在EMTP-ATP仿真计算软件中，将绝缘子串等效为压控开关，绝缘子串的冲击闪络电压U50%则为开关的导通电压，计算中考虑到感应电压Ui的影响，故压控开关导通电压设置为（U50%-Ui），其中感应电压Ui的计算公式参见GB/T 50064“交流电气装置的过电压保护和绝缘配合设计规范”。2.5 冲击接地电阻杆塔冲击接地电阻的范围一般在7-5，仿真建模时冲击接地电阻选择为10。2.6 避雷器的模型由于氧化锌电阻片具有非常优异的非线性伏安特性，可以实现避雷器的无间隙无续流，且造价低廉，因此氧化锌避雷器（MOA）已经得到越来越广泛的应用。本次仿真变电站内避雷器用金属氧化物避雷器，其参数设置如表：表 某1000KV避雷器伏安特性I(KA)0.11245102030U(KV)650760800834850935108212002.7 互感器模型在雷电陡波作用下，电磁电压互感器用入口等值电容模型，此500kV站内电压互感器等效入口电容为3000pF。2.8 变压器模型在雷电陡波作用下，变压器用入口等值电容模型，此500kV站内变压器等效入口电容为5000pF。2.9 断路器模型对于断路器，在分闸状态下，有断口电容和两端触头的对地电容，合闸状态下，只有两端触头的对地电容，其等效模型如图 2-4 所示，图 b 中开关的开闭分别代表断路器的开关状态。取C1 =140pF，C2 =830pF。 a 电容等效模型 b 断路器不同状态等效模型图 2-4 断路器雷电冲击下等效模型3. ATP-EMTP建模根据上述分析,在ATP中搭建的仿真模型如图所示：A 整体电路图B 杆塔模型电路图C 500KV变电站电路图4. 计算结果与分析4.1 反击过电压4.1.1 雷击点的影响在正常运行方式（1、2馈线两回进线，3馈线接电源）下进行雷电流的仿真实验，即上面四个断路器闭合，其他打开，単母线的运行，寻找合适的雷击位置，对各种雷击位置下的各设备电压进行测量，包括断路器（CB1、CB2、CB3、CB4）、主变（TR400）、馈线（LINE1、LINE2、LINE3、LINE4、LINE5）、电压互感器（PT1、PT2、PT3）、母线1末端（BUS1、BUS2），受雷击输电线路杆塔顶端（TWR）上的过电压，实验计算结果，各设备的电压如表3所示（本算例中并未体现0塔）：表1不同受雷位置下，各设备的最高电压，及避雷器中的最大电流大小电压(MV)#1塔#2塔#3塔#4塔#5塔#6塔#7塔TWR3.3483.3483.3483.3483.3482.992TR4001.0930.9540.8550.7670.8220.8130.848LINE12.7961.9651.5841.3201.1361.0021.010PT11.6751.1111.0731.0360.9990.9630.857LINE22.6150.8400.8450.8020.8300.8481.055PT21.4340.8100.8260.7780.7920.7890.849LINE30.8690.8370.8300.8060.8320.7780.837PT30.8710.8220.8030.7860.7970.7790.782DB11.8521.3331.1821.0991.2341.0261.160DB21.3060.9530.8710.8120.80407580.826DB31.1970.9330.8490.8350.7720.8620.756DB41.4231.2431.1621.1770.9601.0701.112BUS11.3391.1170.9400.9020.9000.7860.849BUS21.2061.0390.9820.9200.7700.9460.837电流(KA)MOA127.438.8857.9046.9525.9675.0422.889MOA221.172.0822.3041.2401.6191.5673.016MOA34.3532.3121.8701.5021.6821.0501.202MOATR7.5683.4901.1260.7661.6171.3891.938由结果看出，雷击#1塔塔顶会在变电站形成最严重的侵入波过电压，这种想法在某些情况下可能是正确的，但在我国，大多数情况下不合适。大量研究表明，#1塔和变电站的终端门型构架（也称#0 塔）（本算例中并未体现）距离一般较近，再加上门型构架的冲击接地电阻比较小，雷击#1塔塔顶时，经地线由#0塔返回的负反射波很快返回#1塔，降低了#1塔电位，使侵入波过电压减小。而#2塔、#3塔离#0塔较远，受负反射波的影响较小，过电压较高。所以选择计算雷击#1塔侵入波过电压不合适。另外，进线段各杆塔的塔型、高度、绝缘子串的伏秒特性、杆塔接地电阻不同，也影响着雷击进线段各杆塔时的侵入波过电压。在雷击#6塔时出现了电压反升的现象，且雷击点越远，最高过电压出现的时刻越大（从5us一直到50us不等），推断是电压折反射叠加的结果，根据经验一般选择近区雷，由于雷击#2塔或#3塔时的过电压较高。本算例中后续仿真中选择#2塔作为雷击点。4.1.2 避雷器对过电压大小的影响配置一：将加在变压器出口处避雷器隐藏，仅配置进线处避雷器。表2 配有3个进线避雷器，各设备的最高电压，及避雷器的最大电流大小最高过电压大小(MV)TWTR400LINE1PT1LINE2PT2LINE3PT33.3481.1311.9651.1110.8480.8550.8450.833最高过电压大小(MV)最高过电流大小(KA)DB1DB2DB3DB4MOA1MOA2MOA3MOATR1.3351.0370.9401.2948.8853.0663.311配置二：将加在进线处避雷器隐藏，仅配置变压器出口处避雷器。表3 配有1个变压器避雷器，各设备的最高电压，及避雷器的最大电流大小最高过电压大小(MV)TWTR400LINE1PT1LINE2PT2LINE3PT33.3481.1201.9641.4401.1731.1611.0801.080最高过电压大小(MV)最高过电流大小(KA)DB1DB2DB3DB4MOA1MOA2MOA3MOATR1.5611.1721.0781.6096.124配置三：将加在其余进线处与变压器的避雷器隐藏，仅配置进线1出口处避雷器。表4 配有1个进线避雷器，各设备的最高电压，及避雷器的最大电流大小最高过电压大小(MV)TWTR400LINE1PT1LINE2PT2LINE3PT33.3481.2621.9651.1111.1031.1031.1651.165最高过电压大小(MV)最高过电流大小(KA)DB1DB2DB3DB4MOA1MOA2MOA3MOATR1.3541.1911.1431.6778.885配置四：不配置避雷器表5不配有避雷器，各设备的最高电压，及避雷器的最大电流大小最高过电压大小(MV)TWTR400LINE1PT1LINE2PT2LINE3PT33.3481.5601.9641.4391.3791.3621.4741.474最高过电压大小(MV)最高过电流大小(KA)DB1DB2DB3DB4MOA1MOA2MOA3MOATR1.6011.4231.4052.063二四对比可以看出，变压器出口处避雷器对变压器过电压作用明显；三四对比，一四可以看出，进线出口处避雷器对全网均有影响；综上可以看出，加入避雷器可以有效降低过电压，且接在线路入口处可以同时降低变电站各处电压，一进线出口处避雷器起最大作用。一三对比，可以看出避雷器对雷电流的幅值有一定的衰减，但是衰减程度较小，只有过电压幅值较高时才有明显的衰减，这是因为MOA的伏安特性决定的，在电压较低时，避雷器几乎不流过电流，因此起不到抑制设备过电压的作用。同时由于采用的是间隙式的避雷器，阈值电压的设定也很重要，在不同时刻开始泄放电流对过电压的值也有影响，算例中阈值为1100KV，阈值越低，放电时刻越早，对过电压的抑制越有利，但阈值有一定要求标准，否则正常运行或其他工况时也会放电，本算例中不细究。4.1.3不同电源初始值对过电压大小的影响仅电源的初始相角，在不带避雷器方式下进行雷电流的仿真实验，实验计算结果，各设备的电压如表8所示：表8不同入口电容下，雷电流波形的变化电压(MV)最大过电压TWTR400LINE1PT1LINE2PT2PT303.3481.5601.9651.4401.3931.3801.468903.1261.4732.3901.3972.1161.3511.4271803.4311.4591.4141.4462.1691.5131.378DB1DB2DB3DB4BUS1绝缘子闪络01.6011.4271.3702.0261.490/1.702一A901.6861.3811.4311.8531.482/1.539一A、二C1801.5611.5691.5921.8071.490/1.706二C在0时，#2塔，一回线A相绝缘子击穿；在90时，#2塔，一回线A相，二回线C相绝缘子击穿；在180时，#2塔，二回线C相绝缘子击穿；由于不同的初始相位下，各绝缘子串到达闪络的时刻不同，将影响各设备的过电压值，上表数据0时，最大过电压均来自于一回A相，而180时，均来自二回C相，故绝缘子串的闪络将很大程度上影响设备过电压，变压器的绝缘水平通常低于互感器，并且由于本算例的初衷是为了检验变压器出口处过电压，故选用0出相位。4.1.4 接地电阻对过电压大小的影响在母线及线路上都不加装避雷器（效果较明显），在正常运行方式下进行雷电流的仿真实验，实验计算结果，各设备的电压如表9所示：表9不同接地电阻情况下，各设备的最高电压大小电压(MV)最大过电压TWTR400LINE1PT1LINE2PT2PT3全部取103.3481.5591.9641.4401.3931.3801.468#2塔取73.2371.2481.9371.2701.1081.0941.179全部取73.2371.2271.9361.2681.0991.0881.157#2塔取153.5402.0182.0631.8552.0481.8821.908全部取153.5402.0462.0641.9042.0511.9441.950从表中可以看出接地电阻的大小对雷电压的大小有较大影响，且靠近雷击点的杆塔起主要作用，接地电阻越大，系统中各处的雷电压均越大。其次土壤电阻率对雷电流大小也有影响，在接地电阻全部取10情况下，改变土壤电阻率300*m ，500*m，1000*m，如表17所示。表10不同接地电阻情况下，各设备的最高电压大小电压(MV)最大过电压TWTR400LINE1PT1LINE2PT2PT3300*m3.3121.5922.0361.4811.3951.4081.499500*m3.3481.5591.9641.4401.3931.3801.4681000*m3.4521.5231.8711.2681.3351.3151448可以看出土壤电阻率小将使塔顶电压减小，但会提升变电站内的电压。查阅资料了解到，高土壤电阻率的地方易于形成连续的放电。而土壤电阻率低的地方，如粘土、泥炭和耕地等，落雷时雷电流比较大，因而感应过电压也比较高，易于造成事故。4.1.5不同运行方式下对设备过电压大小的影响方式一：单母线，1、2进线，3接电源方式二：单母线，1进线，2停运，3接电源方式三：单母线，1停运，2进线，3接电源方式四：双母线，1、2进线，3进线，2连接1母线断路器检修表11. 不同运行方式下，各设备的最高电压，及避雷器中的最大电流大小电压(MV)方式1方式2方式3方式4TWR3.3483.0273.4783.348TR4000.9540.9210.8370.837LINE11.9652.0432.3221.965PT11.1110.9841.0461.111LINE20.8401.9212.0150.826PT20.8101.0280.9790.795LINE30.8370.9020.8190.863PT30.8220.8620.8210.832DB11.3331.2531.891/1.7781.627DB20.9531.494/1.6400.9650.943DB30.9331.0450.9640.847/0.928DB41.2431.0361.026BUS1_N11.1171.0711.1071.062BUS1_N21.0391.1120.9760.908BUS2_N10.880BUS2_N20.957电流(KA)MOA18.8855.64111.568.870MOA22.0828.5175.9102.805MOA32.3123.5652.3751.809MOATR3.4904.5242.3062.294表11 不同运行方式下，各设备的最高电压，及避雷器中的最大电流大小（不带避雷器）电压(MV)最大过电压TWTR400LINE1PT1LINE2PT2PT3方式一3.3481.5601.9651.4401.3931.3801.468方式二3.0272.0692.0421.6952.5672.1781.916方式三3.4781.7342.7602.4502.0141.8601.715方式四3.3481.2631.9641.4371.2931.3021.313DB1DB2DB3DB4BUS1BUS2方式一1.6011.4271.3702.0261.490/1.702方式二2.2852.7611.9132.5091.980/2.185方式三3.1601.6111.6341.9240.821/0.619方式四2.1411.3121.1901.338/1.3111.284/1.342由结果对比可知，变电站内部设备接线方式的大小对过电压的计算结果有较大的影响，接线方式越大，主变上所出现的过电压幅值就越低。有两方面的原因，一是因为雷击杆塔放电的总能量是一定的，如果变电站接线方式越复杂，设备越多，则网络的总电容量就越大，各电气设备上出现的过电压数值就越低。比如，与运行方式一比较，运行方式四由于使用了2条母线，对雷电流分流将增大，因此致使侵入波过电压降低；又如，运行方式二与运行方式一比较，由于只有一回进线，从避雷器泄放的电流与杆塔顶端电压可以看出，过电压明显升高。另一方面，电气设备上所出现的过电压幅值，往往是多个来自不同支路的行波叠加的结果，这些电压波有正极性，也有负极性的。接线方式越大，支路越多，行波的数量也越多，但波峰和波峰或波谷和波谷相遇的机会却很少，因而使得叠加的电压幅值并不高。在研究变电站的保护方案时，宜选择过电压最严重的运行方式，既一线一母一变来进行计算，这样得到的过电压水平较高，而由此确定的防雷保护方案也比较可靠。从结果还可以看出，由于有避雷器的保护作用，变压器、电压互感器的过电压都在允许范围内，而断路器之间的电压却始终呈现偏高，特别是最左侧一回进线处的上下两个断路器，在运行方式二与三时，过电压水平已经超出了500KV一般绝缘水平。在运行方式二中，避雷器泄放的电流超过10KA也需引起注意。4.2 绕击对过电压的影响变电站的雷电侵入波除了反击侵入波外，还有雷电流绕过避雷线击中导线而形成的绕击侵入波。计算绕击侵入波过电压时，最重要的参数是最大绕击电流Im。根据电气几何模型，当雷电流幅值增大到Im时，雷电流或击中避雷线，或击中大地，而不会发生绕击，Im称为最大绕击电流。结合实际进线段杆塔的结构尺寸，按下式可计算出最大击距rsm： 式中：hs:避雷线高度， m；hc:导线高度， m；: 避雷线保护角求出rsm后，根据雷电流幅值与击距的关系，求出最大绕击电流Im：式中 k6.72； p0.8。 根据假设的进线段杆塔计算出其最大绕击雷电流幅值为一般不超过40kA，雷电通道的波阻抗取800。由于绕击雷电流幅值远小于反击雷电流幅值，因此绕击时，变电站各电气设备上的雷电过电压水平要低于反击时的过电压水平，所以说绕击侵入波过电压对变电站的威胁不大。绕击和反击的过电压对比如图所示。表12 绕击情况下，各设备的最高电压大小电压(MV)最大过电压TWTR400LINE1PT1BUS1DB1反击3.3481.5591.9641.4401.4901.601绕击0.7620.6231.954