分析变电站空载母线波经过数值.docx

《分析变电站空载母线波经过数值.docx》由会员分享，可在线阅读，更多相关《分析变电站空载母线波经过数值.docx（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、分析变电站空载母线波经过数值hesp导语：我国一直把继电保护等二次设备安装在变电站主控楼内。将这些二次设备下放到变电站开关场内，不仅可以把电流、电压信号就地信息化，而且具有重要的经济价值。1引言长期以来，我国一直把继电保护等二次设备安装在变电站主控楼内。将这些二次设备下放到变电站开关场内，不仅可以把电流、电压信号就地信息化，而且具有重要的经济价值。在三峡输变电建立工程中将有一批500kV变电站的二次设备下放到开关场内。制约二次设备下放的一个重要技术问题是电磁干扰EMI问题，一些兴旺国家十分是美国电力科学研究院EPRI，近二十年来对此问题进展过较为深化的研究1。而国内的研究只处于起步阶段。根据国

2、内外经历，变电站二次设备较严重的干扰源主要有：变电站开关操纵、短路故障、辐射电磁场以及自然界的雷电作用。其中尤以用断路器和隔分开关切、合短空载母线产生的EMI最为严重。此时，母线上产生的电磁波经过不仅通过传导耦合，而且通过空间电磁场辐射耦合到二次设备。所以为了研究空载母线波经过产生的辐射电磁场，首先应计算母线上的波经过。国内外输电线路波经过的分析、计算方法已近成熟2，但对于变电站空载母线的波经过计算方法的研究报道极少。本文采用多导体传输线MTL理论建立了变电站空载母线波经过的计算模型，讨论了求解MTL的时域有限差分FDTD法，对无损耗、无负载的变电站空载母线波经过进展了数值分析，计算结果可用于

3、变电站开关场内二次设备EMI问题的研究。b2多导体传输线的根本方程/b均匀MTL模型如图1所示3。对于非均匀MTL，一般可以分段均匀化，然后按照均匀MTL处理4或者按照卷积特征线方法处理5。MTL模型可用时域电报方程描绘如下4img=355,116zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/40-1.jpg/img式中z轴为传输方向；V、I分别为线上z点时刻t的电压和电流列矢量；VF、IF分别为z处时刻t的鼓励电压源和电流源列矢量；L、C、R、G分别为MTL的单位长电感、电容、电阻和电导矩阵。这些分布参数矩阵可以通过电磁场方法根据MTL

4、几何尺寸和媒质参数计算。img=347,99zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/40-2.jpg/imgb3MTL的FDTD法/bFDTD是一种应用非常广泛的电磁场数值计算方法6，所需内存少、易处理复杂物体、算法简单，十分适于电磁场的时域分析。Paul教授首先应用FDTD法对无分支的MTL进展了计算，并与实际测试进展了比拟，获得较好的计算结果7。将MTL按图2进展离散，沿线电压离散为NDZ1个点，沿线电流离散为NDZ个点。时域电报方程1和2通过空间时域差分近似离散为：img=389,200zszl.cepee/cepee_kjl

5、w_pic/files/wx/zgdjgcxb/2000-6/40-3.jpg/img4式中z为沿线空间步长，t为时间步长；下标表示空间位置点，式3中k取1到NDZ的整数，式4中k取1到NDZ1的整数；上标为时间序列，n为自然数。电压量上标为零、电流量上标为1/2时表示初始值，始端边界条件I0=0，INDZ+1=0。显然，式3和4为一组空间-时间差分方程，可以采用迭代法求解。迭代经过采用如图3所示的跳蛙式方法7，依次对电流、电压进展计算。为了保证算法稳定，要求tzv。其中v为电磁波在MTL中传播的最大形式速度，可用MTL形式分析法获得3。img=343,365zszl.cepee/cepee_

6、kjlw_pic/files/wx/zgdjgcxb/2000-6/40-4.jpg/img图3FDTD迭代经过b4FDTD算法验证/b在输电线路波经过分析中，多采用Bergeron法。由于该法将传输线等效为诺顿电路或者戴维南电路，按电路理论进展计算2，4，所以只能求解节点处的电压和电流波经过。对于传输线上的某一离散点，应人为设置一个假想的节点，因此难以求解传输线上所有离散点的电压和电流波经过。为验证FDTD法的有效性，分别应用FDTD法和Bergeron法计算图4所示三导体无损耗传输线的串音问题，其中，VS2为单位阶跃电压源，RS1、RS2、RL1、RL2分别为50，线长05m，初始条件为零

7、，其单位长电感和电容矩阵可计算得：img=344,212zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/40-5.jpg/img图4两线串音问题图5分别为FDTD法和Bergeron法的计算结果，其中曲线A、B、C和D分别为图4中相应的四点的电压波形。从图5可以看出，FDTD法和Bergeron法所得结果根本吻合，只是FDTD法产生了一些毛刺，但不会导致大的误差。img=378,147zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/41-1.jpg/img图5FDTD法Bergeron法

8、计算结果比拟5变电站空载母线的波经过计算图6为某变电站500kV空载母线简化模型，其中只考虑一组母线和引线，且母线上不带负载，忽略线路和构架的影响。L2、L3为母线，L1为引线，L150m，L290m，L350m，h16m。母线间隔65m，引线间隔8m。母线采用LGJQT1400型特种轻型钢芯铝绞线。在模型中，忽略母线间的邻近效应和两端的边沿效应，且设母线无损耗，大地为无限大完纯导体平面，其母线单位长电感和电容矩阵分别为img=388,339zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/41-2.jpg/img显然，与图1不同的是，图6中

9、的MTL具有分支构造。为此将三段MTL分别离散，使分支点2为电压离散点，并利用分支点2的电压和电流边界条件，可推导获得如下2点的FDTD法迭代公式img=354,32zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/41-3.jpg/imgimg=390,278zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/41-4.jpg/img图6变电站母线简化模型将分支点2的方程56和各段MTL的迭代方程34联立，并利用初始条件和边界条件即可进展具有分支构造MTL的波经过数值分析。设A、B和C相电源

10、为单位余弦对称三相电压源，电源内阻设为零，空载母线初始状态为零。时刻t0时对空载母线加电。采用FDTD法计算线上各离散点的电压、电流波经过。图7、图8和图9分别给出了A、B和C相母线和引线上1、2、3和4四个点的电压波形，图10给出A相母线和引线上四个点的电流波形。img=360,434zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/41-5.jpg/img图9C相四点电压波形从图可以看出，由于屡次反射、透射的影响，使得电压、电流波形极其复杂。其中B、C相电压波形在所研究的时间内几乎一致。由于此时线间耦合影响相差不大，在所研究的较短时间内，

11、B相和C相的电压源根本一样。在图7中，由于A相电压源峰值为1V，所以在其内阻为0时，1点的电压在所研究的较短时间内保持不变。当向母线充电时，随着时间的变化，L1上的各点按照光速，依次带上1V的电压，当到达分支点2时，由于阻抗不匹配，此时发生波的反射和透射，反射系数为033，到L2和L3的透射系数均为066，此时B、C相的电压也传播到分支点2处，对A相也有影响，但影响较小。于是L2、L3上的电压变为066V，时间为01668s。然后，在L2、L3上分别以066V的电压传播下去，而在L1上出现的反射波电压为033V。当03337s时，L3上的电压波到了终端3，于是发生反射，反射系数为1，所以电压值

12、变为原来的2倍。同时，L1上的波到1点，由于电压源的限制其电压应为1V，所以反射系数为1。此后到了0467s，L2的终端4发生反射，反射系数为1，电压值变为原来的2倍。到了050s时，L1、L2上的波又到了2点，由于其反射和透射，使得电压发生变化，为原来的2倍。此后，循环往复进展下去。在所研究的时间内线路上最大的电压可到原来电压源幅值的29倍，即到29V，且作用时间极短，约为0066s。img=359,116zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/42-1.jpg/img图10A相四点电流波形b、c、d分别表示L1、L2、L3侧的点

13、2应该指出，以上计算结果是在假设MTL无损耗、大地为理想导体面的情况下获得的，所得结果虽不能代表实际情况，但其波经过的规律应和实际情况相似。对电压波经过作快速傅立叶变换，可以进一步得到线上各离散点的电压和电流波经过频谱。图11为图7中线4的频谱，可以看出变电站500kV空载母线加电瞬间，电磁干扰的频谱主要集中在10MHz的频率范围内，这与现场实际测量结果是一致的8。img=333,152zszl.cepee/cepee_kjlw_pic/files/wx/zgdjgcxb/2000-6/42-2.jpg/img图11A相点4的电压频谱b6结论/b采用多导体传输线对变电站空载母线建模是可行的。应

14、用文中提出的具有分支构造多导体传输线的时域有限差分公式，可对变电站空载母线波经过进展有效的数值计算。其计算结果可进一步应用于变电站二次设备下放后的电磁干扰问题研究。b参考文献：/b1EPRITR102006：ElectromagnetictransientsinsubstationsRVolI、，PaloAlto，CA，19932吴维韩，张芳榴，等电力系统过电压数值计算M北京：科学出版社，19893TescheFM，IanozMV，KarlssonTEMCanalysismethodsandcomputationmodelsMNewYork：JohnWileySonsPress，19964Pa

15、ulCRAnalysisofmulticondutortransmissionlinesMNewYork：JohnWileySonsPress，19945MaoJunma，LiZhengFanAnalysisofthetimeresponseofMTLwithfrequencydependentlossesbythemethodofconvolutioncharacteristicsJIEEEonMTT，1992，404：6376446KunzKS，LuebbersRJTheFinitedifferencetimedomainmethodinelectromagneticsMBocaRaton，FL：CRCPress，19937PaulCRFiniteDifferenceTimeDomainanalysisoflossytransmissionlinesJIEEEonEMC，1996，381：15238何彬电力系统二次设备的电磁兼容问题J中国电力，0