专利名称：直流输电线下小型接地物体表面合成电场的计算方法
技术领域：
本发明属于超高压和特高压直流输电线路电磁环境分析技术领域，尤其涉及一种直流输电线下物体表面合成电场的计算方法。
背景技术：
我国是世界上能源分布和经济发展极不平衡的国家之一。国家电网公司提出了通过建设以特高压电网为骨干网架的坚强智能电网，以实现全国能源资源的优化配置，保障国家能源安全的发展战略。直流输电系统在这一过程中扮演着重要角色，我国目前已经建成了大量直流工程，包括葛南、天广、三广、三常、贵广、三沪、灵宝，向上等多项大型工程，而且未来十年内，还会有一批直流工程将会兴建。高压直流输电线路运行中不可避免的产生电晕现象，进而由其产生的定向运动的空间电荷(称为离子流)将大大加强线路导线电荷产生的静电场(称为标称电场)。这种由离子流和线路导线电荷共同作用产生的电场称为合成电场。直流输电线路线下的地面最大合成电场可能达到地面标称电场的3 3. 5倍，是影响直流输电线路电磁环境的一个重要因素。线下人员活动时的合成电场分布问题是直流输电电磁环境的核心问题，当人体位于直流输电线下合成电场最大位置附近，人体接地电阻较小时，人体电位近似视为零，人体畸变电场最大，工程中主要关注此时的合成电场分布，该问题可归结为直流输电线下小型接地物体的合成电场计算问题。随着社会发展和环保意识的增强，人们对其关注度也逐渐上升。但由于问题过于复杂，一直未能得到很好的分析和计算。应用三维合成电场的计算方法，解决该问题代价仍然巨大，不能满足大量工程分析的需要，因此迫切需求一种简单， 准确的计算方法。
发明内容
本发明的目的在于，提供一种直流输电线下小型接地物体表面合成电场的计算方法，用于解决目前使用的计算方法过于复杂的问题。为实现上述目的，本发明提供的技术方案是，一种直流输电线下小型接地物体表面合成电场的计算方法，其特征是所述方法包括步骤1 在直流输电线下不存在小型接地物体时，利用二维离子流场计算方法计算空间各点的电荷密度；步骤2 在直流输电线下存在小型接地物体时，利用所述二维离子流场计算方法计算的空间各点的电荷密度，使用三维静电场计算方法获得物体表面合成电场。所述步骤1具体是步骤11:设定边界条件；步骤12 设定空间各点的电荷密度值；步骤13 由空间各点的电荷密度值，计算空间各点的电位和场强；
步骤14 :由空间各点的电位和场强，计算空间各点的电荷密度；步骤15 判断边界条件和空间各点的电荷密度是否满足收敛条件，如果边界条件和空间各点的电荷密度满足收敛条件，则所述空间各点的电荷密度即为所求解；否则，执行步骤16 ；步骤16 修正空间各点的电荷密度，并将修正后的空间各点的电荷密度作为空间各点的电荷密度值，返回步骤13。所述步骤2具体是步骤21 设定三维区域，使物体位于三维区域底部中心位置；步骤22 以所述二维离子流场计算方法计算的空间各点的电荷密度作为所述三维区域对应空间位置的点的电荷密度；步骤23 获取所述三维区域的各边界电位、物体表面电位；步骤M 利用计算物体表面电场，即为所求合成电场。本发明利用二维离子流场计算方法结合简单的三维静电场计算方法确定直流输电线路下方的电场强度，在确保计算结果精确性的前提下，降低了计算复杂度。
图1是利用上流有限元法计算空间各点的电荷密度示意图；图2是某接地物体的示意图；图3是本发明提供的方法和三维上流有限元离子流场算法沿图2中A线计算合成电场的结果对比图；图4是本发明提供的方法和三维上流有限元离子流场算法沿图2中B线计算合成电场的结果对比图；图5是本发明提供的方法和三维上流有限元离子流场算法沿图2中C线计算合成电场的结果对比图；图6是本发明提供的方法和三维上流有限元离子流场算法沿图2中D线计算合成电场的结果对比图；图7是本发明提供的方法和三维上流有限元离子流场算法沿图2中E线计算合成电场的结果对比图；图8是本发明提供的方法和三维上流有限元离子流场算法沿图2中F线计算合成电场的结果对比图；图3-图8中，方法1为本发明提供的方法，方法2为三维上流有限元离子流场算法。
具体实施例方式下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。实施例1给定直流线路为双极六分裂士SOOkV的直流输电线路，线路高度18米，极间距22 米，导线分裂间距45cm，子导线型号为LGJ-630。
本发明提供的直流输电线下小型接地物体表面合成电场的计算方法包括步骤1 在直流输电线下不存在物体吋，利用ニ维离子流场计算方法计算空间各 点的电荷密度。根据描述直流离子流场的基本方程，正负离子在离子流场中运动的速度有三方面 的因素影响在电场力作用下离子迁移、离子扩散和风速。忽略电荷扩散的影响后，其正负 离子的运动速度可表示为v+=K+Es+w(1)v_ = -K_Es+w(2)描述双极输电线路离子流场的基本方程可表达为
V2cp = -(p+-p_)/s0(3)
Es = -V^(4)j+ = p+v+(5)j = p V(6)V 又=-Rp+p— Ie(7)=Rp+p_le(8)对于单极直流线路(以正极性为例)的离子流场，由于只存在ー种空间电荷，不存 在正负电荷的复合过程，所以其描述方程可简化为
V2CP = -Pjs0(9)
Es = -V (p(10)V./ =0(11)]+ = p+v+(12)式中，尾为场強，又、又分别为正负离子流密度，P+、P-分别为正负空间电荷密 度，K+、K_分別是正负离子迁移率，R为离子的复合系数，#为风速矢量，e C1为真空介电常数， e是基本电子电量，タ是空间任意一点的电位。上述式(3)-(8)中各点的空间电荷密度P +和P—以及电位タ相互耦合，但又都是 未知的，因此本文下面将采用迭代法进行求解，其具体过程是步骤11 设定边界条件。由微分方程理论可知，上述方程要有足够的定解条件才能确定方程的解。由于本 发明将根据泊松方程的有限元法计算空间各点的电位和场强，而有限元法的计算需要将输 电线路的无限域转化为有限域考虑，因此具体到本发明，需要设定ー个计算场域。具体方法 是在离导线较远处认为划出一个边界，称作人工边界。这样，导线表面、地面和人工边界所 围成的封闭平面区域即成为计算场域，而泊松方程的边界条件也就是导线表面、地面和人 エ边界三种边界所需满足的条件。通常，人工边界的半径取大于等于4倍的导线高度为宜。 据此，三种边界的条件为
i)导线表面滑=士￠/ii)地面<3 = 0iii)假定人工边界离导线足够远，以至于空间电荷对该处电场的影响与导线电荷相比可基本忽略，则在人工边界上有-.(P = UmOiv)再将Kapazov假设(即导体起晕后表面场强认为保持在起晕场强不变)作为
边界条件之一，则有·β^ = Ε +
on其中，U为直流输电线路运行电压，Uffl为标称场电压，Eoi分别为正负极起晕场强。步骤12 设定空间各点的电荷密度值。步骤13 由空间各点的电荷密度值，计算空间各点的电位和场强。由空间各点的电荷密度值，可以根据泊松方程的有限元法计算空间各点的电位， 进而求出场强。步骤14 :由空间各点的电位和场强，计算空间各点的电荷密度。本实施例采用上流有限元法计算空间各点的电荷密度。首先寻找和判断上流元。 如图1所示，若要根据已知电荷密度的点j和点m来求点i的电荷密度，就必须判断与i相关的三角形ijm是否构成上流元。逆时针旋转，节点i的下一个点定义为j，第二个点定义为m。向量^和5分别为向量-和;/i逆时针旋转90°所得到，若j点的空间电荷迁移速度f 与 -石和的夹角都小于90°，则该三角元为上流元。符合这个条件的的方向限制在图1的虚线区域，即^和，延长线所围成的区域。相应的数学描述为bjVx+CjVy ^ 0(13)bmVx+cmVy 彡 0(14)这样计算出来的结果满足0< PiCmadPp Pm)，即速度下方点的电荷密度总是小于速度上方的点，保证了解的收敛，也与事实相符。在三角形i jm内对电荷密度进行插值
权利要求
1.一种直流输电线下小型接地物体表面合成电场的计算方法，其特征是所述方法包括步骤1 在直流输电线下不存在小型接地物体时，利用二维离子流场计算方法计算空间各点的电荷密度；步骤2 在直流输电线下存在小型接地物体时，利用所述二维离子流场计算方法计算的空间各点的电荷密度，使用三维静电场计算方法获得物体表面合成电场。
2.根据权利要求1所述的直流输电线下物体表面合成电场的计算方法，其特征是所述步骤1具体是步骤11 设定边界条件；步骤12 设定空间各点的电荷密度初值；步骤13 由空间各点的电荷密度值，计算空间各点的电位和场强； 步骤14 由空间各点的电位和场强，计算空间各点的电荷密度； 步骤15 判断边界条件和空间各点的电荷密度是否满足收敛条件，如果边界条件和空间各点的电荷密度满足收敛条件，则所述空间各点的电荷密度即为所求解；否则，执行步骤 16 ；步骤16 修正空间各点的电荷密度，并将修正后的空间各点的电荷密度作为空间各点的电荷密度值，返回步骤13。
3.根据权利要求1所述的直流输电线下物体表面合成电场的计算方法，其特征是所述步骤2具体是步骤21 设定三维区域，使物体位于三维区域底部中心位置； 步骤22:以所述二维离子流场计算方法计算的空间各点的电荷密度作为所述三维区域对应空间位置的点的电荷密度；步骤23 获取所述三维区域的各边界电位和物体表面电位； 步骤M 计算物体表面电场，即为所求合成电场。
全文摘要
本发明公开了超高压和特高压直流输电线路电磁环境分析技术领域中的一种直流输电线下小型接地物体表面合成电场的计算方法。包括在直流输电线下不存在小型接地物体时，利用二维离子流场计算方法计算空间各点的电荷密度；在直流输电线下存在小型接地物体时，利用所述二维离子流场计算结果中空间各点的电荷密度，使用三维静电场计算方法获得物体表面合成电场。本发明在确保计算结果精确性的前提下，降低了计算复杂度。
文档编号G06F19/00GK102368286SQ201110324649
公开日2012年3月7日 申请日期2011年10月21日 优先权日2011年10月21日
发明者卢铁兵, 崔翔, 甄永赞 申请人:华北电力大学