一种计算直流输电线路三维离子流场的特征线方法技术

本发明专利技术公开了一种计算直流输电线路三维离子流场的特征线方法，能够用于计算线路交叉跨越、地面建筑周围、以及绝缘子周围等的直流三维离子流场。首先给定一个假设的空间电场和空间电荷的分布，然后由空间电荷通过泊松方程计算新的空间电场分布，画特征线；再由新的空间电场分布通过离子流方程和电流连续性方程沿特征线求解特征线上的电荷分布；如此反复迭代，直至前后两次空间电荷的分布在允许的误差范围内，根据电荷分布进行合成电场和离子流密度的计算。本发明专利技术相比Deutsch假设法，考虑了空间电荷的影响，计算结果更准确；相比有限元法，大大地增加了电荷更新的效率，在计算效率上明显优于有限元法。

A characteristic line method for calculating three dimensional ion flow field of HVDC transmission lines

The invention discloses a DC transmission line characteristics for calculation of ion flow field method can be used in DC ion flow field calculation of line crossing, the ground around the building, and around the insulator etc.. The distribution of space electric field and space charge first given a hypothesis, and then by the space charge by the Poisson equation to calculate the electric field distribution of the new features, painting line; by the new space electric field distribution along the line on the line to solve the characteristic charge characteristics by ion flow equation and current continuity equation; so until the iteration. Before and after the two space charge distribution in the allowable error range, according to the calculation of total electric field and ion current density distribution of charge. Compared with the Deutsch hypothesis method, the invention considers the influence of space charge, and the calculation result is more accurate; compared with the finite element method, the efficiency of charge regeneration is greatly increased, and the calculation efficiency is obviously better than that of the finite element method.

【技术实现步骤摘要】
一种计算直流输电线路三维离子流场的特征线方法
本专利技术涉及高电压输电技术和电磁场计算领域，尤其涉及一种计算直流输电线路三维离子流场的特征线方法，其用于计算交叉跨越、地面存在建筑、绝缘子周围等直流输电线路由于电晕造成的三维离子流场。
技术介绍
在特高压直流输电工程的建设过程中，线路的电磁环境是关系系统设计、建设以及运行的重大问题。由于线路的电晕现象，周围的空气发生电离形成带电粒子。在双极导线之间和极导线与大地之间的整个区域内，将充满空间电荷，对由导线电荷产生的标称电场产生畸变作用，总体上的效应使合成电场增加。当空间电荷运动到地面后被地面吸收，就形成了离子流密度。这就是直流输电线路的离子流场。一般简单地考虑输电线路的物理模型，可以将离子流场问题视为一个二维问题。目前针对直流输电线路的二维离子流场计算方法已经有很多。然而考虑到线路下方跨越建筑物、植被、农业大棚等各种物体时，离子流场受其影响的畸变，准确的计算需要在三维尺度下进行。从方程形式上来说，三维离子流场的控制方程与二维时没有本质区别，因此目前三维离子流场的计算方法都从二维发展而来，分为基于Deutsch假设的解法和有限元数值解法。2010年，华北电力大学的罗兆楠等人，仍以Deutsch假设在三维下计算了输电线路下方的建筑物附近的合成电场，同时，将方法应用到了跨越森林植被时，合成电场对树木灼伤安全净空距离的确定，以及人体模型感知的合成电场和离子流密度。在有限元数值解法上，2011年，甄永赞等人以二维上流有限元法为基础，发展了三维上流有限元法，提出一种按照电荷运动方向更新节点电荷密度的方法。但计算时，仅在建筑物周围进行了三维网格剖分计算，其余区域保持二维计算。上流有限元法的稳定性较好，但当计算区域是复杂模型或大尺寸的建筑物时，计算量巨大，需要花费的代价高，因此目前的应用仅局限于单极性的离子流场，无法完全适应和推广到双极输电线路的离子流场这一实际问题。总体来说，Deutsch假设的解法由于假定实际空间电场的方向不受空间电荷产生的电荷的影响，计算结果很不精确。而三维有限元由于计算设计的区域大，边界复杂，意味着有限元法网格剖分量大，计算迭代需要占用大量的存储资源，无法针对实际中的各类环境和线路规格的开展三维离子流场计算。针对以上问题，本专利技术基于离子流场控制方程的特性，提出了计算直流输电线路三维离子流场的特征线方法。
技术实现思路
针对上述问题中存在的不足之处，本专利技术提供一种计算直流输电线路三维离子流场的特征线方法。为实现上述目的，本专利技术提供一种计算直流输电线路三维离子流场的特征线方法，包括：步骤1、空间网格剖分；步骤2、基于给定的线路参数，用模拟电荷法求解标称电场，并对节点电荷密度初始化；步骤3、根据节点电荷密度计算合成电场强度，画特征线；步骤4、沿步骤3中的特征线，求解特征线上的电荷分布，更新节点电荷密度；步骤5、重复步骤3、步骤4进行迭代，并比较前后两次电荷密度；若在允许的误差范围内，则认为达到稳定解，此时的电荷密度为离子流场的电荷分布；否则，修正节点电荷密度，返回步骤2；步骤6、根据电荷分布进行后处理，所述后处理包括合成电场和离子流密度的计算。作为本专利技术的进一步改进，该直流输电线路三维离子流场基于下述三个方程：泊松方程：式中，E是合成电场强度，ρ+和ρ-分别是正、负电荷密度，ε0是空气介电常数；离子流方程：式中，J+和J-分别是正、负离子流密度，k+和k-分别为正、负离子迁移率，W为矢量，表示风速；电流连续性方程：式中，R是正负离子的复合系数，e是电子电荷量1.602×10-19C。作为本专利技术的进一步改进，在步骤4中，根据式(12)、(13)求解特征线上的电荷分布，其具体推导方法为：以式(1)～(3)为基本方程，进行特征线法求解的方程推导；由于正、负电荷的方程类似，仅给出正电荷方程的推导过程，负电荷方程可同理得出；把式(2)代入(3)，得到：若风速W为常数，将式(1)代入式(4)中，消去得到：因此，得到关于ρ+和ρ-的梯度方程：为明确起见，三维的方程推导以负电荷为例；记式(6)展开为：等式两边都除以vx，得到：由式(8)知，三维下关于负电荷密度的偏微分方程依然可以转化为沿着特征线上关于x的常微分方程，正电荷密度也可作类似推导；其中，特征线的定义为：当不考虑风速时，W＝0，特征线即为电场线：依然以负电荷方程为例，如同二维的情况，对从导线到地面的特征线，没有正电荷，式(8)简化为：其中，M＝(ε0vx)-1为常数，若特征线起点即导线表面的电荷密度已知，微分方程式(11)的定解为：对从一极导线到另一极导线的特征线，各点都存在正、负电荷，若假设正电荷ρ+已知，式(8)存在解析解：其中：C是由特征线起点的负电荷密度，微分方程的边界条件确定的系数；因此，沿各类特征线，各点的负电荷密度可直接解析求得。作为本专利技术的进一步改进，在步骤5中，允许的误差范围设定为1％以内。与现有技术相比，本专利技术的有益效果为：1、特征线法有解析计算公式，计算速度更快、更准；2、特征线法相比Deutsch假设法，考虑了空间电荷的影响，计算结果更准确；3、特征线法相比有限元法，大大地增加了电荷更新的效率，在计算效率上明显优于有限元法；4、可以使用积分方程法计算由空间电荷产生的空间电场，使用本专利技术的特征线法计算空间电场导致的空间电荷的移动，从而在整个计算中均可以仅使用解析计算公式，不需要建立大规模的方程组，对计算机内存的要求不高。附图说明图1为本专利技术一种实施例公开的计算直流输电线路三维离子流场的特征线方法的流程图；图2为本专利技术一种实施例公开的两条垂直交叉跨越的±800kV直流输电线路图；图3为本专利技术一种实施例公开的特征线法交叉跨越线路地面电场的分布图；图4为本专利技术一种实施例公开的I回导线下电场对比图，其中：图4(a)中y取-10m，图4(b)中y取10m，图4(c)中y取-30m，图4(d)中y取30m，图4(e)中y取-100m，图4(f)中y取100m，图4(g)中y取0m；图5为本专利技术一种实施例公开的II回导线下电场对比图，其中：图5(a)中x取-10m，图5(b)中x取10m，图5(c)中x取-30m，图5(d)中x取30m，图5(e)中x取-110m，图5(f)中x取110m，图5(g)中x取0m。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。下面结合附图对本专利技术做进一步的详细描述：本专利技术提供一种计算直流输电线路三维离子流场的特征线方法，能够用于计算线路交叉跨越、地面建筑周围、以及绝缘子周围等的直流三维离子流场。空间电场由下面的三个方程描述：泊松方程：式中，E是合成电场强度，ρ+和ρ-分别是正、负电荷密度，ε0是空气介电常数；离子流方程：式中，J+和J-分别是正、负离子流密度，k+和k-分别为正、负离子迁移率，W为矢量，表示风速；电流连续性方程：式中，R是正负离子的复合系数，e是电子电荷量1.602×10-19C。本文档来自技高网...

【技术保护点】
一种计算直流输电线路三维离子流场的特征线方法，其特征在于，包括：步骤1、空间网格剖分；步骤2、基于给定的线路参数，用模拟电荷法求解标称电场，并对节点电荷密度初始化；步骤3、根据节点电荷密度计算合成电场强度，画特征线；步骤4、沿步骤3中的特征线，求解特征线上的电荷分布，更新节点电荷密度；步骤5、重复步骤3、步骤4进行迭代，并比较前后两次电荷密度；若在允许的误差范围内，则认为达到稳定解，此时的电荷密度为离子流场的电荷分布；否则，修正节点电荷密度，返回步骤2；步骤6、根据电荷分布进行后处理，所述后处理包括合成电场和离子流密度的计算。

【技术特征摘要】
1.一种计算直流输电线路三维离子流场的特征线方法，其特征在于，包括：步骤1、空间网格剖分；步骤2、基于给定的线路参数，用模拟电荷法求解标称电场，并对节点电荷密度初始化；步骤3、根据节点电荷密度计算合成电场强度，画特征线；步骤4、沿步骤3中的特征线，求解特征线上的电荷分布，更新节点电荷密度；步骤5、重复步骤3、步骤4进行迭代，并比较前后两次电荷密度；若在允许的误差范围内，则认为达到稳定解，此时的电荷密度为离子流场的电荷分布；否则，修正节点电荷密度，返回步骤2；步骤6、根据电荷分布进行后处理，所述后处理包括合成电场和离子流密度的计算。2.如权利要求1所述的计算直流输电线路三维离子流场的特征线方法，其特征在于，该直流输电线路三维离子流场基于下述三个方程：泊松方程：▽·E＝(ρ+-ρ-)/ε0(1)式中，E是合成电场强度，ρ+和ρ-分别是正、负电荷密度，ε0是空气介电常数；离子流方程：式中，J+和J-分别是正、负离子流密度，k+和k-分别为正、负离子迁移率，W为矢量，表示风速；电流连续性方程：式中，R是正负离子的复合系数，e是电子电荷量1.602×10-19C。3.如权利要求2所述的计算直流输电线路三维离子流场的特征线方法，其特征在于，在步骤4中，根据式(12)、(13)求解特征线上的电荷分布，其具体推导方法为：以式(1)～(3)为基本方程，进行特征线法求解的方程推导；由于正、负电荷的方程类似，仅给出正电荷方程的推导过程，负电荷方程可同理得出；把式(2)代入(3)，得到：若风速W为常数，将式(1)代入式(4)中，消去▽·E，得到：因此，得到关于ρ+和ρ-的梯度方程：

【专利技术属性】
技术研发人员：张波，何金良，莫江华，肖凤女，曾嵘，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京,11