摘要：提出了一种基于矩量法的在频域下计算超高压输电铁塔附近三维电场分布的数值方法。该方法在频域下使用复电阻率的概念，将空气和土壤联合起来视为多层介质，待求变量为各段导体的漏电流，激励源既可以为电压源又可以为电流源。该方法可以分析频域下由输电线路、铁塔以及铁塔接地部分产生的三维电场分布。对同一例题下的由该方法和电力系统电磁分析软件包CDEGS计算的结果进行了比较，两者结果是吻合的。在500 kV输电线路的实测结果也验证了该方法的有效性。

关键词：超高压输电 三维电场 频域 矩量法 电力系统

1引言

　　运行输电线路附近存在较高的电场，这些电场可能对周围物体和公用走廊的其他线路产生影响。电压等级越高，影响越大。为了人身和设备的安全，超高压输电线路产生的电场正日益受到人们的重视。为了在输电线路上进行安全带电作业，需选择操作方式并采取绝缘防护措施，这就需要对输电铁塔附近的电场分布进行研究。在电力铁塔上同塔悬挂自承式全介质光缆（ADSS）时，由于这种光缆在一定电场强度作用下会发生腐蚀（电蚀），这也需要分析铁塔周围的电场分布从而找出合适的悬挂范围。对于即将建设的西北750kV输电线路，研究输电线路及铁塔周围的电场分布有着重要的现实意义。虽然可以现场实测输电线路产生的电场，但由于工作环境的限制，具体测量难度较大，铁塔附近的测量更加复杂，且不可能对所有铁塔都进行测量，因此需要进行模拟计算。不考虑铁塔时，某些简单的输电线路的电场也可以采用近似的解析公式进行估计。但铁塔附近的电场由于实际问题比较复杂，应当使用数值方法进行计算。

　　国外出于高场强下带电作业的需要，开始研究输电线路产生的电场强度及由此引起的生物效应，并取得了一定的成就^[1]。我国目前对此问题也非常重视，但多数计算方法只考虑输电线路本身的电场分布，未充分考虑铁塔的影响^[2,3]。然而实际输电线路中铁塔附近的工频电场有较大的畸变，因此在计算输电线路的电场时必须考虑铁塔的影响。

　　通常模拟电荷法用于计算输电线路产生的电场是比较有效的^[1,4]。但将大地视为非完纯导体或铁塔没有与架空地线相连及需要将铁塔的接地部分同时考虑时，模拟电荷法计算电场比较复杂，且模拟电荷法一般用来分析直流模型。本文基于矩量法^[5]提出了一种在频域下计算输电线路及铁塔附近的三维电场分布的数值方法，将空气视为多层介质中的一层并考虑各层介质的电阻率和介电常数，将各层介质的电阻率在频域下视为复数，各段导体的泄漏电流为待求变量，从而能够将铁塔、铁塔接地系统、输电线路和避雷线同时考虑。

　　2 基本方法

　　2.1 复电阻率及多层介质中的格林函数

　　由于本文使用矩量法来计算频域下输电线路及铁塔附近的电场分布，并以各段导体的泄漏电流为待求变量，因此需要首先得到包括空气和多层土壤在内的多层介质中的格林函数。

　　空气和土壤都具有一定的电阻率和介电常数，分析交流情况下铁塔周围的电场分布时它们都会起作用。由麦克斯韦方程组的第一式

　　式中H和E 分别为介质的磁场强度和电场强度；ρ和ε分别为介质的电阻率和电容率。引入新的电流密度

　　使用上面的复电阻率，也可以将空气视为一层导电介质，这样得到的格林函数就可以用来同时分析由输电线路铁塔及其接地部分产生的电场。求此格林函数的方法与求多层土壤中的格林函数的方法类似^[6]，但要注意拉普拉斯方程的边界条件的设定。图1为多层介质示意图。

　　1）空气和土壤之间的分界面上的边界条件与土壤层之间的分界面上的边界条件相同，均为

　　2）当纵坐标z趋向于正无穷或负无穷时，电位V应当趋向于零。

　　本文采用计算速度快，计算精度高^[7]的复镜像法求解格林函数。

　　2.2 使用矩量法计算铁塔的电场分布

　　作者已根据矩量法建立了多层土壤中的大型接地网的频域分析模型^[8]，该模型以各段导体的泄漏电流为待求变量，与目前国内的计算方法和软件相比待求量比较少，其计算结果与国外的接地分析软件包的CDEGS^[9]结果一致，一个小接地网上的测试结果也证明了该方法的实用性。采用这个模型和上面介绍的复电阻率的概念来分析超高压输电线路铁塔附近的三维电场。本文与文献[8]的相同之处是使用导体表面上电位的连续性来建立方程组，即导体段表面上两点间的电位差是由各导体段上的泄漏电流决定的，而导体段内这两点间的电位差是由导体的自阻抗和流过导体段上的轴向电流产生的，这两个电位差应相等，且导体段上的轴向电流可用各导体段上的泄漏电流表示，从而建立了以各段导体的泄漏电流为待求变量的线性方程组，解之可得导体各段上的泄漏电流。通过这些泄漏电流可以求得超高压输电线路铁塔附近的三维电场。本文与文献[8]的区别在于使用的格林函数和激励源的不同：

　　（1）文献[8]使用的格林函数是多层土壤的格林函数，不考虑空气，也无须使用复电阻率的概念；而本文推导多层导电媒质中的格林函数时，将空气考虑在内，且使用2.1小节中推导得出的复电阻率。