一种基于上流无网格方法的换流变空间电荷的计算方法技术

本发明专利技术请求保护一种基于上流无网格方法的换流变压器空间电荷计算方法，包括确定特高压换流变绝缘油纸材料数学非线性模型，以及二维轴对称换流变空间电荷计算模型；针对换流变内部模型尺寸差异较大，采用改进无网格方法进行空间电荷分布的计算，该计算模型可根据尺寸自适应调整布点方式，对换流变空间电荷分布进行计算，在临近计算边界区域自适应调整局部子域大小，避免其与全局边界交叉，在全局边界采用配点法进行计算。本发明专利技术可实现特高压换流变空间电荷计算，避免网格剖分和非规则边界局部子域与全局域交界部分判识，可有效减少计算时间、提高计算精度。

A Computation Method of Converter Space Charge Based on Upstream Meshless Method

The invention requests to protect a space charge calculation method of converter transformer based on upstream meshless method, including determining the mathematical non-linear model of UHV converter insulation oil-paper material and two-dimensional axisymmetric converter space charge calculation model; aiming at the large difference in the size of the internal model of converter, the improved meshless method is adopted to calculate the space charge distribution, and the calculation is carried out. The model can adaptively adjust the distribution of space charge according to size, adjust the size of local sub-domain adaptively in the adjacent computing boundary area, avoid intersecting with the global boundary, and use the collocation method to calculate the global boundary. The invention can realize space charge calculation of UHV converter flow, avoid mesh generation and identification of local and global boundary parts of irregular boundary, effectively reduce calculation time and improve calculation accuracy.

【技术实现步骤摘要】
一种基于上流无网格方法的换流变空间电荷的计算方法
本专利技术属于(特)高压直流换流变压器内绝缘电场计算和优化设计领域，尤其涉及一种换流变空间电荷分布的计算方法。
技术介绍
特高压输电工程是我国实现大容量、远距离输电的重要发展方向。特高压换流变压器作为变换与传递电能、隔离交流与直流线路的核心装备.其质量与可靠性对于保障我国输电系统的安全与稳定起着至关重要的作用。绝缘问题是导致变压器故障的主要原因。较之交流变压器，换流变压器承受更为复杂的电场作用形式，运行工况更加严酷，因而更容易产生绝缘破坏与失效。同时，特高压工程所需电压等级更高，对换流变压器的绝缘性能、可靠性与稳定性提出了更高要求。加之国内外对此缺乏实际的运行经验和试验数据，因而，迫切需要展开有关换流变压器绝缘失效规律与机理的研究。空间电荷是引起绝缘故障与失效的重要因素。由于换流变压器工作环境的特殊性，其承受的电压类型也有交流电压、直流电压和极性反转电压等多种类型，所以对换流变内油纸绝缘结构的设计提出了更高的要求。在直流电场的作用下，油纸绝缘材料内部存在着载流子注入、复合、入陷、脱陷等复杂过程，介质内部空间电荷积聚、运动和消散会直接影响到电场的分布，削弱或加强油纸绝缘内部的局部场强。因此，计算空间电荷对电场的影响对于换流变油-纸绝缘特性的研究至关重要。目前，关于油纸绝缘系统空间电荷的研究主要集中于3个方面：变压器油中流注放电发展过程的数值仿真.基于电声脉冲法(pulsedelectro—acoustic，PEA)的油纸绝缘介质内空间电荷的测量，以及复合电场下换流变压器中的空间电荷浓度与电场强度分布数值仿真。其中，流注过程仿真与PEA测量属于机理研究，复合电场仿真属于工程应用研究。
技术实现思路
本专利技术旨在提供一种针对特高压换流变压器空间电荷分布规律的上流无网格方法，分析交直流复合作用下换流变绝缘部件的电场和空间电荷分布规律，从而得到换流变空间电荷分布的特点与规律，从而得以指导生产设计。提出了一种基于上流无网格方法的换流变空间电荷的计算方法。本专利技术的技术方案如下：一种基于上流无网格方法的换流变空间电荷的计算方法，其包括以下步骤：步骤1、通过三电极法确定特高压换流变计算模型结构参数及非线性的材料属性，并根据运行工况，确定其电场特性，建立二维轴对称空间电荷计算模型；步骤2、根据空间电荷计算模型的结构形状，对计算模型进行布点离散，并设置相应的材料参数和边界条件；步骤3、设置计算初始值，对极板采用肖特基电子发射模型，根据双极性空间电荷运动模型和空间电荷陷阱理论，基于上流无网格方法求解载流子的输运方程，求解对流方程得到换流变内部的电荷密度，并将电荷密度代入泊松方程计算空间电场强度，在边界处调整节点局部子域的大小，避免与全局边界相交，对边界节点采用配点法进行计算；步骤4、通过不断迭代泊松方程求解电场强度和上流元法求解电荷密度这两个过程，得到不同时刻电荷密度和场强的瞬态结果。进一步的，所述步骤1)二维轴对称空间电荷计算模型具体包括：肖特基发射模型(如公式8-9)，正负载流子的复合消散模型(如公式10-13)。进一步的，所述步骤2)中的电场和空间电荷瞬态空间电场的求解控制方程描述如下：1)控制方程s(x,t)＝seu+set+shu+sht(4)n0＝neu+net+nhu+nht(5)γ＝qnμ(6)其中，ε为介电常数；ρ为空间电荷密度；j为空间电荷运动运动形成的电流密度，E为电场强度；t是时间；n为电荷密度；s为载流子源项；set，seu，sht，shu为入陷电子，自由电子，入陷空穴，自由空穴的源项；net，neu，nht，nhu为入陷电子，自由电子，入陷空穴，自由空穴的电荷密度；r为电导率；q为电荷电量，μ为迁移率。模型内边界为正极板，施加电压V0；外边界为负极板，电位为零电位，用数学公式表示为：进一步的，所述步骤3)边界节点采用基于RPIMp(附加多项式基的径向基函数)形函数的配点法进行计算。进一步的，所述肖特基极板电荷发射模型方程及载流子的输运方程描述如下：1)肖特基极板电荷发射模型为：极板表面载流子的发射受电场强度、温度以及势垒的影响，其中，Je(t)和Jh(t)分别为阴极电子和阳极空穴注入的电流密度；T是开氏温度，K；A为Richardson常数，计算中取120A/(cm2·K2)；ωei和ωhi分别为电子和空穴的注入势垒；Ec(t)和Ea(t)分别为阴极和阳极表面的电场强度，ε表示介电常数；k为普朗克常数；2)载流子的输运方程其中，eμ，et，hμ，ht分别代表自由电子、入陷电子、自由空穴、入陷空穴；neμ、nht、net、nhu分别为自由电子、入陷空穴、入陷电子、自由空穴浓度；Seμ,ht、Seμ,hu、Set,hu、Set,ht表示为复合系数；Be、Bh分别是电子和空穴的入陷系数；Net，Net0分别是电子陷阱浓度和空穴陷阱浓度。进一步的，所述上流无网格方法计算换流变内部的电场强度和电荷密度方程描述如下：1)泊松方程见2)正负空间电荷的迁移速度：γ＝qnμ(6)式。本专利技术的优点及有益效果如下：本专利技术的创新点具体在于确定特高压换流变绝缘油纸材料数学非线性模型，以及二维轴对称换流变空间电荷计算模型；针对换流变内部模型尺寸差异较大，采用改进无网格方法进行空间电荷分布的计算，该计算模型可根据尺寸自适应调整布点方式，对换流变空间电荷分布进行计算，在临近计算边界区域自适应调整局部子域大小，避免其与全局边界交叉，在全局边界采用配点法进行计算。本专利技术可实现特高压换流变空间电荷计算，避免网格剖分和非规则边界局部子域与全局域交界部分判识，可有效减少计算时间、提高计算精度。附图说明图1是本专利技术提供优选实施例提供的基于上流无网格方法的换流变空间电荷计算流程图。图2为本专利技术实施例提供的油纸绝缘材料中四种载流子的复合入陷示意图。图3为本专利技术实施例提供的上流无网格法的上流节点分析的图。图4为本专利技术实施例提供的弯曲油纸绝缘结构的节点剖分图。图5为本专利技术实施例提供的一维截线下油纸界面的电荷密度分布示意图。图6为本专利技术实施例提供的考虑空间电荷下换流变角环直流试验电场分布云图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本专利技术的一部分实施例。本专利技术解决上述技术问题的技术方案是：基于上流无网格方法的换流变空间电荷计算。此方案包含以下步骤：步骤1、确定特高压换流变计算模型结构参数及材料属性，并根据运行工况，确定其电场特性，进而建立其二维轴对称空间电荷计算模型；步骤2、根据计算模型的结构形状，对计算模型进行布点离散，同时为了加强计算效率和精度，在尺寸较小区域采用密集布点的方式，在边界也同样设计较多的节点，以满足电场和空间电荷计算精度对离散方式的需求，并设置相应的材料参数和边界条件。步骤3、针对换流变承受交直流复合作用情况，考虑其中的双极电荷输运和空间电荷陷阱理论。采用肖特基极板电荷发射模型，计算极板表面载流子的发射电流密度。步骤4、利用上流无网格方法分别计算换流变内部的电场强度和电荷密度。首先，利用泊松方程求解电场，其次，利用上流无网格法根据电场计算某一时刻各点的电荷密度。在求解空间电荷分布时，由于不同时刻得到的电场强度本文档来自技高网...

【技术保护点】
1.一种基于上流无网格方法的换流变空间电荷的计算方法，其特征在于，包括以下步骤：步骤1、通过三电极法确定特高压换流变计算模型结构参数及非线性的材料属性，并根据运行工况，确定其电场特性，建立二维轴对称空间电荷计算模型；步骤2、根据空间电荷计算模型的结构形状，对计算模型进行布点离散，并设置油纸绝缘结构的材料参数和边界条件；步骤3、设置计算绝缘材料的温度、场强初始值，对极板采用肖特基电子发射模型，根据双极性空间电荷运动模型和空间电荷陷阱理论，基于上流无网格方法求解载流子的输运方程，求解对流方程得到换流变内部的电荷密度，并将电荷密度代入泊松方程计算空间电场强度，在边界处调整节点局部子域的大小，避免与全局边界相交，对边界节点采用配点法进行计算；步骤4、通过不断迭代泊松方程求解电场强度和上流元法求解电荷密度这两个过程，得到不同时刻电荷密度和场强的瞬态结果。

【技术特征摘要】
1.一种基于上流无网格方法的换流变空间电荷的计算方法，其特征在于，包括以下步骤：步骤1、通过三电极法确定特高压换流变计算模型结构参数及非线性的材料属性，并根据运行工况，确定其电场特性，建立二维轴对称空间电荷计算模型；步骤2、根据空间电荷计算模型的结构形状，对计算模型进行布点离散，并设置油纸绝缘结构的材料参数和边界条件；步骤3、设置计算绝缘材料的温度、场强初始值，对极板采用肖特基电子发射模型，根据双极性空间电荷运动模型和空间电荷陷阱理论，基于上流无网格方法求解载流子的输运方程，求解对流方程得到换流变内部的电荷密度，并将电荷密度代入泊松方程计算空间电场强度，在边界处调整节点局部子域的大小，避免与全局边界相交，对边界节点采用配点法进行计算；步骤4、通过不断迭代泊松方程求解电场强度和上流元法求解电荷密度这两个过程，得到不同时刻电荷密度和场强的瞬态结果。2.根据权利要求1所述的基于上流无网格方法的换流变空间电荷的计算方法，其特征在于，所述步骤1)二维轴对称空间电荷计算模型具体包括：肖特基发射模型，正负载流子的复合消散模型。3.根据权利要求1所述的基于上流无网格方法的换流变空间电荷的计算方法，其特征在于，所述步骤2)中的电场和空间电荷瞬态空间电场的求解控制方程描述如下：1)控制方程s(x,t)＝seu+set+shu+sht(4)n0＝neu+net+nhu+nht(5)γ＝qnμ(6)其中，ε为介电常数；ρ为空间电荷密度；j为空间电荷运动运动形成的电流密度，E为电场强度；t是时间；n为电荷密度；s为载流子源项；set，seu，sht，shu为入陷电子，自由电子，入陷空穴，自由空穴的源项；ne...

【专利技术属性】
技术研发人员：池骋，高兵，杨帆，廖瑞金，屈凯峰，司战超，
申请(专利权)人：重庆大学，北京云道智造科技有限公司，
类型：发明
国别省市：重庆,50