CN104850738B - 架空电力线路杆塔处雷电感应电压计算方法 - Google Patents
架空电力线路杆塔处雷电感应电压计算方法 Download PDFInfo
- Publication number
- CN104850738B CN104850738B CN201510214027.5A CN201510214027A CN104850738B CN 104850738 B CN104850738 B CN 104850738B CN 201510214027 A CN201510214027 A CN 201510214027A CN 104850738 B CN104850738 B CN 104850738B
- Authority
- CN
- China
- Prior art keywords
- lightning
- line
- grid
- tower
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Suspension Of Electric Lines Or Cables (AREA)
Abstract
本发明请求保护一种架空电力线路杆塔处雷电感应电压的计算方法。主要内容：1.在考虑杆塔影响情况下的空间雷电电磁场计算。选择采用基于亚网格技术的时域有限差分法，对空气空间进行粗网格剖分，对杆塔进行细网格剖分。2.将雷电电磁场二维计算结果应用于架空电力线路与雷电流回击通道三维空间结构的场量分解。雷电回击电流激发形成的电磁场。3.架空电力线路雷电感应电压计算及杆塔位置雷电感应电压提取。采用Agrawal模型和FDTD法计算架空线路上的雷电感应电压；通过网格查询定位于杆塔位置，提取杆塔处的雷电感应电压。本发明能够有效且更加准确地评估架空电力线路和杆塔处绝缘子、避雷器等所承受的雷电感应电压。
Description
技术领域
本发明涉及架空电力线路防雷、绝缘配合计算领域，具体为一种架空电力线路杆塔处雷电感应电压计算方法。
背景技术
雷击大地时，雷电释放出的脉冲电流将在数千米的空间范围内产生强烈的辐射电磁场。架空输配电线路距离长，裸露于空气当中，极易与雷电电磁场相耦合，形成线路过电压，严重的会造成线路对地或相间闪络、损坏端接变压器及开关设备等。特别是对于35kV及以下的架空配电线路而言，由于绝缘水平较低，雷电感应过电压已成为线路绝缘闪络、跳闸增多的主要原因。正确评估架空电力线路，特别是位于杆塔处的绝缘子、避雷器等装置所承受的雷电感应电压的强度和时变特性，对于线路防雷、确保其安全稳定运行具有重要意义。
已有文献均在忽略杆塔影响的情况下线路雷电感应过电压的计算，所得计算结果是存在缺陷的。因为根据电磁波传播的基本原理可知，在电磁波传播的路径上出现两种媒质时，由于电磁参数发生突变，会出现反射、折射等现象。因此，当雷电电磁场触及杆塔时，会发生畸变，这样的电磁场耦合到线路上所形成的感应电压也随之发生变化。
若要更准确地计算架空电力线路以及安装在杆塔上的绝缘子、避雷器等所承受的雷电感应电压，从而更合理地进行绝缘配合和防雷设计，就有必要考虑杆塔的影响。
发明内容
针对现有技术中的问题，本发明提供了一种考虑杆塔影响的情况下更加准确地计算架空电力线路所承受的雷电感应电压，为线路导线、绝缘子选型和线路绝缘配合设计提供理论依据的架空电力线路杆塔处雷电感应电压计算方法，本发明的技术方案如下：一种架空电力线路杆塔处雷电感应电压计算方法，其包括以下步骤：
101、获取多导架空线路雷电流波形特征参数，包括电流峰值I0、波头时间τ1和波尾时间τ2、回击速度、衰减常数，建立基底雷电流模型is(0,t)和回击通道雷电流is(z',t)的数学模型，其中基底雷电流模型is(0,t)采用Heidler模型函数：
式中：I0用于限定电流峰值；为峰值电流修正系数；τ1、τ2分别表示波头、波尾时间。
回击通道雷电流is(z',t)的数学模型采用改进传输线模型：
i(z',t)＝e-λz'is(0,t-z'/vi) t≥z'/vi (2)
式中：z′为雷电通道高度；λ为衰减常数；vi为雷电流回击速度；
102、获取空气、杆塔的介质参数，具体为在空气中的介电常数ε、磁导率μ和电导率σ；在杆塔中介电常数εf、磁导率μf和电导率σf。设置二维柱坐标系，对计算区域进行亚网格剖分，即对空气区域进行粗网格剖分，对杆塔进行细网格剖分，代入步骤101中的回击通道雷电流模型is(z',t)，根据亚网格时域有限差分法即Sub-grid FDTD，迭代计算雷电电磁场；
103、获取多导架空线路空间的位置参数，即在设定的计算坐标系中，获取多导架空导线所处的坐标位置，以及线路与雷击点的垂直距离d，根据架空电力线路与雷电流回击通道三维空间关系，将二维柱坐标系中计算所得的雷电场分量Er沿架空线轴线方向进行分解转换Ex，此外，根据多导架空线路空间位置参数，计算多导架空线对应的单位长度电感系数矩阵L和单位长度电容系数矩阵C；
104、根据步骤103得到的雷电电场分量沿架空线路轴向的分解计算结果，及单位长度电感系数矩阵L和单位长度电容系数矩阵C，并输入架空线路始端和终端等值阻抗、架空线路长度，计算得出雷电电场对多导架空电力线路的电磁耦合的Agrawal模型，形成含有分布电源的传输线方程组；
105、将架空线路沿x方向以步长Δx剖分，时间离散步长取Δt，可得线路上雷电感应电压的迭代计算式
式中：和分别表示第(n+1)Δt时刻线路长度为kΔx处的感应电流列相量和散射电压向量；
获取架空线路始端和终端等值阻抗，修正边界
式中：和分别表示第(n+1)Δt时刻线路长度为kΔx处的感应电流列相量和散射电压向量；ZS和ZL分别为线路始端和终端端接阻抗向量；为导线处入射电场的水平分量向量；和分别为线路始端和终端位置入射电场的垂直分量向量；
106、求线路上的雷电感应电压，即散射电压与入射电压之和
然后根据架空线路起止位置、档距、杆塔与线路始端距离、网格剖分步长信息，确定杆塔所在位置对应的计算网格数，提取和输出杆塔处的雷电感应电压。
进一步的，步骤102中采用粗网格和细网格对二维雷电电磁场进行划分计算，的具体步骤为：
A、粗网格的空间步长和时间步长选择；Δlmax＜λmin/10,
式中：Δlmax为网格的最大空间步长；λmin为雷电电磁波的最小波长；
时间步长的选择依据Courant数值稳定条件
B、细网格的空间步长和时间步长选择；
将一个粗网格作nf等分，则细网格的空间步长为
细网格的时间步长为
C、采用中心差分离散Maxwell方程对粗网格上雷电电磁场分量的迭代计算；
D、基于亚网格FDTD方法的细网格上雷电电磁场分量的迭代计算。
进一步的，步骤C中采用中心差分离散Maxwell方程，得空间内雷电电磁场的FDTD迭代计算式
式中：i，j分别为r方向和z方向上的网格数，n为迭代时间步数；εc、μc和σc分别为空气的介电常数、磁导率和电导率。
而雷电流回击通道上的Ez需要做特别处理：
式中：表示(n+1)Δt时刻、高度为处的回击通道雷电流大小。
进一步的，步骤D中采用亚网格FDTD方法对细网格上雷电电磁场分量的迭代计算具体为：
d1、在杆塔所在位置进行细网格剖分，将原有一个粗网格作nf等分，Erf和Ezf表示细网格上的雷电电磁场在方向上的分量；
d2、对细网格从n→n+1时刻的推进分为nf个小时间步进行，即细网格内部的场分量计算进行nf次循环迭代计算，迭代计算公式为
式中：ii，jj分别为r方向和z方向上的细网格数，m为细网格内迭代时间步数；εf、μf和σf分别为杆塔介质材料的介电常数、磁导率和电导率。
进一步的，步骤d2中还包括通过插值法对粗细网格重合边界处的场分量的计算步骤，具体为利用n时刻和n+1时刻已知的场分量值以及细网格时间迭代步数进行插值，插值需满足波动方程
式中：D表示Er、Ez或其中之一；v为波速。
本发明的优点及有益效果如下：
本发明提出了基于亚网格技术的时域有限差分法(FDTD)，在考虑杆塔影响的情况下计算架空电力线路上形成的雷电感应电压的理论方法，这样既可以获得较准确的计算结果，又尽可能减小计算量、节省计算时间。
本发明还提出通过雷电场分量分解，将雷电电磁场二维计算结果应用于架空线路和雷击通道三维结构场-线耦合计算的技术方法，达到简化计算、提高效率的目的。上述要点弥补了现有理论研究的空白。
附图说明
图1是本发明优选实施例架空电力线路杆塔处雷电感应电压计算流程图；
图2是雷电电磁场计算粗网格剖分示意图；
图3是雷电电磁场计算细网格剖分示意图；
图4是亚网格FDTD计算程序实现流程图；
图5是将Er分解为Ex的坐标及网格设置；
图6是雷电电磁场与多导架空线路耦合的Agrawal模型。
具体实施方式
以下结合附图，对本发明作进一步说明：
1.基于FDTD亚网格技术的雷电电磁场计算
说明：鉴于雷电电磁场的瞬变特性，选择时域有限差分法(FDTD)进行计算；进一步地，因为在大尺度空间电磁场的计算过程中需要考虑小尺度的杆塔结构，因此采用了FDTD的亚网格技术。
1.1获取雷电流波形特征参数，包括电流峰值、波头时间和波尾时间、回击速度、衰减常数等，建立基底雷电流模型is(0,t)和回击通道雷电流is(z',t)的数学模型。
基底雷电流推荐使用Heidler模型函数：
式中：I0用于限定电流峰值；为峰值电流修正系数；τ1、τ2用于限定波头、波尾时间。
回击通道雷电流推荐使用改进传输线模型：
i(z',t)＝e-λz'is(0,t-z'/vi) t≥z'/vi (2)
式中：z′为雷电通道高度；λ为衰减常数；vi为雷电流回击速度。
1.2获取空气、杆塔的介质参数，具体为各自的介电常数ε、磁导率μ和电导率σ。
1.3粗网格的空间步长和时间步长选择
利用FDTD方法进行均匀空间雷电电磁场计算时，网格剖分的尺度在式(3)所示稳定性要求的前提下可适当取大一点，这类网格称为粗网格。网格过小会导致过大计算量，严重影响计算效率。
Δlmax＜λmin/10 (3)
式中：Δlmax为网格的最大空间步长；λmin为雷电电磁波的最小波长。
与此同时，时间步长的选择依据Courant数值稳定条件
式中：cmax为最大波速；Δx、Δy、Δz为网格在x、y、z三个方向上的步长。
雷电电磁场在空气中传播，并且计算时取Δx＝Δy＝Δz＝Δlc，则时间步长稳定性条件化简为
式中：c为光速。
1.4细网格的空间步长和时间步长选择
由于杆塔塔身直径相比粗网格的步长Δlc小得多，如要体现杆塔对雷电电磁场的影响，就需要在杆塔位置处进行精细剖分。
假设将一个粗网格作nf等分，则细网格的空间步长为
细网格的时间步长为
1.5粗网格上雷电电磁场分量的迭代计算
雷击大地或建筑物后雷电回击电流形成电磁场，雷电电磁场的空间分布以雷电流通道为轴对称分布，因此可以在二维柱坐标系内进行计算。
满足Δlmax、Δtc的要求在二维柱坐标系中进行网格剖分，如图2所示，其中Er、Ez和表示粗网格上的电磁场分量。
电场分量和磁场分量的取样分别为和采用中心差分离散Maxwell方程，可得空间内雷电电磁场的FDTD迭代计算式
式中：i，j分别为r方向和z方向上的网格数，n为迭代时间步数；εc、μc和σc分别为空气的介电常数、磁导率和电导率。
雷电流回击通道上的Ez需要做特别处理：
式中：表示(n+1)Δt时刻、高度为处的回击通道雷电流大小。
1.6细网格上雷电电磁场分量的迭代计算
在杆塔所在位置进行细网格剖分，如图3所示将原有一个粗网格作nf等分，(以nf＝3为例)。Erf和Ezf表示细网格上的场分量。
在FDTD计算中，磁场分量的计算用到相邻4个电场分量，而电场分量的计算用到相邻2个磁场分量，这不仅适用于粗网格，同样适用于细网格。
进一步地，对细网格从n→n+1时刻的推进需分为nf个小时间步进行，即细网格内部的场分量计算需进行nf次循环迭代计算，迭代计算公式为
式中：ii，jj分别为r方向和z方向上的细网格数，m为细网格内迭代时间步数；εf、μf和σf分别为杆塔介质材料的介电常数、磁导率和电导率。
直接计算粗细网格重合边界处的场分量需要并不存在的细网格场分量，解决该问题的方法是利用n时刻和n+1时刻已知的场分量值以及细网格时间迭代步数，通过插值法进行求解。插值需满足波动方程
式中：D表示Er、Ez或其中之一；v为波速。
基于FDTD亚网格技术的雷电电磁场计算程序实现流程如图4所示。
2.电场分量在架空线路-雷击通道三维结构中的分解计算
(说明：通常情况下，架空电力线路与雷击通道的空间位置是三维结构，为了避免复杂而繁琐的三维雷电电磁场计算，提高计算效率，可将二维柱坐标系内计算的场分量进行分解转换。)
2.1获取多导架空线路空间位置参数，即在设定的计算坐标系中，多导架空导线所处的坐标位置，以及线路与雷击点的垂直距离d等。
2.2设架空线路轴向为x方向，对架空导线进行网格剖分，步长取Δx。
2.3如图5所示为将二维柱坐标系内计算所得Er分解为Ex时所使用的坐标及网格设置。转换关系式满足
式中：i为r方向网格数；k为x方向网格数。
3.多导架空电力线路及绝缘子雷电感应电压计算
3.1根据多导架空线路空间位置参数，计算多导架空线对应的单位长度电感系数矩阵L和单位长度电容系数矩阵C。其中，C＝P-1
P为多导架空线路对应电位系数矩阵，为n阶方阵，n为导线数。矩阵中的元素根据镜像法可求得：
式中：ε0为真空介电常数；hi和Ri分别为i导线离地垂直距离和半径；Lij和L'ij分别为i导线及其镜像导线到j导线的距离。
架空电力线路近似可视为无损耗线，则
L＝μ0ε0P (15)
式中：μ0为真空磁导率。
3.2根据计算雷电电场场对多导架空电力线路的电磁耦合的Agrawal模型(如图6所示)，形成含有分布电源的传输线方程组：
式中：usca(x,t)为散射电压；i(x,t)为线路中的感应电流；Ex(x,h,t)为入射电场在导线处的水平分量。
3.3将架空线路沿x方向以步长Δx剖分，时间离散步长取Δt，以中心差商近似计算式(14)中的偏导，可得线路上雷电感应电压的迭代计算式
式中：和分别表示第(n+1)Δt时刻线路长度为kΔx处的感应电流列相量和散射电压向量。
获取架空线路始端和终端等值阻抗，修正边界
式中：和分别表示第(n+1)Δt时刻线路长度为kΔx处的感应电流列相量和散射电压向量；ZS和ZL分别为线路始端和终端端接阻抗向量；为导线处入射电场的水平分量向量；和分别为线路始端和终端位置入射电场的垂直分量向量。
3.4线路上的雷电感应电压为散射电压与入射电压之和
3.5根据架空线路起止位置、档距、杆塔与线路始端距离、网格剖分步长等信息，确定杆塔所在位置对应的计算网格数，提取和输出杆塔处的雷电感应电压。
以上这些实施例应理解为仅用于说明本发明而不用于限制本发明的保护范围。在阅读了本发明的记载的内容之后，技术人员可以对本发明作各种改动或修改，这些等效变化和修饰同样落入本发明权利要求所限定的范围。
Claims (5)
1.一种架空电力线路杆塔处雷电感应电压计算方法，其特征在于，包括以下步骤：
101、获取多导架空线路雷电流波形特征参数，包括电流峰值I0、波头时间τ1和波尾时间τ2、回击速度、衰减常数，建立基底雷电流模型is(0,t)和回击通道雷电流is(z',t)的数学模型，其中基底雷电流模型is(0,t)采用Heidler模型函数：
<mrow>
<msub>
<mi>i</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mi>&eta;</mi>
</mfrac>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>t</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中：I0用于限定电流峰值；为峰值电流修正系数；τ1、τ2分别表示波头、波尾时间；
回击通道雷电流模型is(z',t)采用改进传输线模型：
is(z',t)＝e-z'/λis(0,t-z'/vi) t≥z'/vi
式中：z′为雷电通道高度；λ为衰减常数；vi为雷电流回击速度；
102、获取空气、杆塔的介质参数，具体为在空气中的介电常数εc、磁导率μc和电导率σc；在杆塔中介电常数εf、磁导率μf和电导率σf；设置二维柱坐标系，对计算区域进行亚网格剖分，即对空气区域进行粗网格剖分，对杆塔进行细网格剖分，代入步骤101中的回击通道雷电流模型is(z',t)，根据亚网格时域有限差分法即Sub-grid FDTD，迭代计算雷电电磁场；
103、获取多导架空线路空间的位置参数，即在设定的计算坐标系中，获取多导架空导线所处的坐标位置，以及线路与雷击点的垂直距离d，根据架空电力线路与雷电流回击通道三维空间关系，将二维柱坐标系中计算所得的雷电场分量Er沿架空线轴线方向进行分解转换为Ex，此外，根据多导架空线路空间位置参数，计算多导架空线对应的单位长度电感系数矩阵L和单位长度电容系数矩阵C；
104、根据步骤103得到的雷电电场分量沿架空线路轴向的分解计算结果，及单位长度电感系数矩阵L和单位长度电容系数矩阵C，并输入架空线路始端和终端等值阻抗、架空线路长度，计算得出雷电电场对多导架空电力线路的电磁耦合的Agrawal模型，形成含有分布电源的传输线方程组；
105、将架空线路沿x方向以步长Δx剖分，时间离散步长取Δt，可得线路上雷电感应电压的迭代计算式
<mrow>
<msubsup>
<mi>i</mi>
<mi>k</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>i</mi>
<mi>k</mi>
<mi>n</mi>
</msubsup>
<mo>-</mo>
<msup>
<mi>&Delta;tL</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>k</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msup>
<mi>&Delta;tC</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>i</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
式中：和分别表示第(n+1)Δt时刻线路长度为kΔx处的感应电流列相量和散射电压向量；表示第nΔt时刻线路长度为kΔx处的入射电场水平方向分量；
获取架空线路始端和终端等值阻抗，修正边界
<mrow>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>C</mi>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>C</mi>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>i</mi>
<mn>1</mn>
<mi>n</mi>
</msubsup>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>hE</mi>
<mrow>
<mi>z</mi>
<mi>S</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mi>C</mi>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>L</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mi>C</mi>
<mrow>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>L</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<msubsup>
<mi>u</mi>
<mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>i</mi>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mi>n</mi>
</msubsup>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>hE</mi>
<mrow>
<mi>z</mi>
<mi>L</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mi>L</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
式中：ZS和ZL分别为线路始端和终端端接阻抗向量；和分别为线路始端和终端位置入射电场的垂直分量向量；
106、求线路上的雷电感应电压，即散射电压与入射电压之差
<mrow>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>u</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>h</mi>
</msubsup>
<msubsup>
<mi>E</mi>
<mi>z</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>d</mi>
<mi>z</mi>
</mrow>
然后根据架空线路起止位置、档距、杆塔与线路始端距离、网格剖分步长信息，确定杆塔所在位置对应的计算网格数，提取和输出杆塔处的雷电感应电压。
2.根据权利要求1所述的一种架空电力线路杆塔处雷电感应电压计算方法,其特征在于，步骤102中采用粗网格和细网格对二维雷电电磁场进行划分计算，的具体步骤为：
A、粗网格的空间步长和时间步长选择；△lmax<λmin/10,
式中：△lmax为网格的最大空间步长；λmin为雷电电磁波的最小波长；
时间步长的选择依据Courant数值稳定条件
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo><</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>/</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
式中：△x、△y、△z分别表示在x、y、z三个方向上进行剖分的步长；cmax表示最大波速；
B、细网格的空间步长和时间步长选择；
将一个粗网格作nf等分，则细网格的空间步长为
<mrow>
<msub>
<mi>&Delta;l</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;l</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>n</mi>
<mi>f</mi>
</msub>
</mfrac>
</mrow>
细网格的时间步长为
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>n</mi>
<mi>f</mi>
</msub>
</mfrac>
</mrow>
式中：△lc、△tc分别表示粗网格的空间步长和时间步长；
C、采用中心差分离散Maxwell方程对粗网格上雷电电磁场分量的迭代计算；
D、基于亚网格FDTD方法的细网格上雷电电磁场分量的迭代计算。
3.根据权利要求2所述的一种架空电力线路杆塔处雷电感应电压计算方法,其特征在于，步骤C中采用中心差分离散Maxwell方程，得空间内雷电电磁场的FDTD迭代计算式
式中：i，j分别为r方向和z方向上的网格数，n为迭代时间步数；εc、μc和σc分别为空气的介电常数、磁导率和电导率；
而雷电流回击通道上的Ez：
式中：表示(n+1)△t时刻、高度为处的回击通道雷电流大小。
4.根据权利要求2所述的一种架空电力线路杆塔处雷电感应电压计算方法,其特征在于，步骤D中采用亚网格FDTD方法对细网格上雷电电磁场分量的迭代计算具体为：
d1、在杆塔所在位置进行细网格剖分，将原有一个粗网格作nf等分，Erf和Ezf表示细网格上的雷电电磁场在r、z方向上的分量；
d2、对细网格从n→n+1时刻的推进分为nf个小时间步进行，即细网格内部的场分量计算进行nf次循环迭代计算，迭代计算公式为
式中：ii，jj分别为r方向和z方向上的细网格数，m为细网格内迭代时间步数；εf、μf和σf分别为杆塔介质材料的介电常数、磁导率和电导率。
5.根据权利要求4所述的一种架空电力线路杆塔处雷电感应电压计算方法,其特征在于：步骤d2中还包括通过插值法对粗细网格重合边界处的场分量的计算步骤，具体为利用n时刻和n+1时刻已知的场分量值以及细网格时间迭代步数进行插值，插值需满足波动方程
<mrow>
<msup>
<mo>&dtri;</mo>
<mn>2</mn>
</msup>
<mi>D</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mfrac>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>D</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
</mrow>
式中：D表示Er、Ez或其中之一；v为波速。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510214027.5A CN104850738B (zh) | 2015-04-29 | 2015-04-29 | 架空电力线路杆塔处雷电感应电压计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510214027.5A CN104850738B (zh) | 2015-04-29 | 2015-04-29 | 架空电力线路杆塔处雷电感应电压计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104850738A CN104850738A (zh) | 2015-08-19 |
| CN104850738B true CN104850738B (zh) | 2018-04-20 |
Family
ID=53850379
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510214027.5A Expired - Fee Related CN104850738B (zh) | 2015-04-29 | 2015-04-29 | 架空电力线路杆塔处雷电感应电压计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104850738B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105426671B (zh) * | 2015-11-11 | 2018-02-16 | 重庆大学 | 一种雷暴天气下架空配电线路的可靠性评测方法 |
| CN106124830A (zh) * | 2016-07-11 | 2016-11-16 | 清华大学 | 基于人工神经网络获得雷电流的方法 |
| CN107845141B (zh) * | 2017-11-27 | 2018-11-27 | 山东大学 | 一种瞬变电磁三维fdtd正演多分辨网格剖分方法 |
| CN108037373B (zh) * | 2017-12-07 | 2020-03-31 | 南方电网科学研究院有限责任公司 | 雷电回击电磁场的检测方法和系统 |
| CN108363899A (zh) * | 2018-05-17 | 2018-08-03 | 陕西省地方电力（集团）有限公司 | 一种配网感应雷过电压评估方法及系统 |
| CN109726448B (zh) * | 2018-12-11 | 2023-03-24 | 中国人民解放军63660部队 | 一种基于迭代法的分段式模式转换器设计方法 |
| CN112163054B (zh) * | 2020-09-09 | 2024-02-02 | 成都深瑞同华科技有限公司 | 一种基于路径图的变压器插入方法 |
| CN112749500B (zh) * | 2021-01-11 | 2023-09-22 | 华能新能源股份有限公司 | 高海拔山区风电场集电线路雷电感应过电压的估算方法 |
| CN113252971A (zh) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | 一种基于fdtd的配电线路感应雷过电压计算方法 |
| CN113255976B (zh) * | 2021-05-12 | 2023-01-24 | 云南电网有限责任公司丽江供电局 | 一种配电线路架设架空地线的接地优化设计方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103714239A (zh) * | 2013-11-15 | 2014-04-09 | 贵州电力试验研究院 | 雷击大地时低压线路绝缘子雷电感应电压计算方法 |
| CN103984992A (zh) * | 2014-05-12 | 2014-08-13 | 重庆大学 | 一种多回交流输电线路纵横交错区域合成磁场的预测方法 |
-
2015
- 2015-04-29 CN CN201510214027.5A patent/CN104850738B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103714239A (zh) * | 2013-11-15 | 2014-04-09 | 贵州电力试验研究院 | 雷击大地时低压线路绝缘子雷电感应电压计算方法 |
| CN103984992A (zh) * | 2014-05-12 | 2014-08-13 | 重庆大学 | 一种多回交流输电线路纵横交错区域合成磁场的预测方法 |
Non-Patent Citations (6)
| Title |
|---|
| Analysis of Lightning-Induced Voltages on Overhead Lines Using a 2-D FDTD Method and Agrawal Coupling Model;He-Ming Ren 等;《IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY》;20081231;第50卷(第3期);第1-6页 * |
| Improvements of an FDTD-Based Surge Simulation Code and Its Application to the Lightning Overvoltage Calculation of a Transmission Tower;Taku Noda 等;《Electric Power Systems Research》;20071231;第651-659页 * |
| 基于 FDTD 的雷击杆塔瞬态响应仿真研究;吴耀辉;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20140615(第6期);全文 * |
| 架空输电线路雷击感应过电压耦合机理及计算方法分析;莫付江 等;《电网技术》;20050331;第29卷(第6期);第72-77页 * |
| 架空配电线路雷电感应过电压计算研究;边凯 等;《中国电机工程学报》;20121105;第32卷(第31期);第191-199页 * |
| 配电线路感应雷过电压计算与防护的研究;徐兴发;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20131015(第10期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104850738A (zh) | 2015-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104850738B (zh) | 架空电力线路杆塔处雷电感应电压计算方法 | |
| CN109711088B (zh) | 一种场路结合的配电线路感应雷过电压计算方法 | |
| CN106940745B (zh) | 基于全波电磁暂态模型的雷电过电压计算方法 | |
| Yang et al. | Lightning impulse corona characteristic of 1000-kV UHV transmission lines and its influences on lightning overvoltage analysis results | |
| Chen et al. | Lightning-induced voltages on a distribution line with surge arresters using a hybrid FDTD–SPICE method | |
| Wang et al. | Breakdown characteristics of long air gap with negative polarity switching impulse | |
| CN103293451B (zh) | 一种高压输电线路杆塔接地装置防雷特性的评估方法 | |
| de Azevedo et al. | Porous Soil Effect on Backflashover Analysis of Power Transmission Lines | |
| Andreotti et al. | An analytical approach to calculation of lightning induced voltages on overhead lines in case of lossy ground—Part II: Comparison with other models | |
| Borghetti et al. | Calculation of lightning-induced voltages on an overhead line taking into account the presence of nearby buildings | |
| Tong et al. | High current field test of impulse transient characteristics of substation grounding grid | |
| Sun et al. | Lightning-induced Overvoltage of 10 kV Distribution Line Based on Electromagnetic Return-stroke Model Using FDTD | |
| Silveira et al. | The influence of attachment height on lightning-induced voltages | |
| Liu et al. | Magnetic field distribution inside metallic grid-like buildings struck by lightning based on finite element method | |
| Huang et al. | A distributed modelling method of the transmission tower and transient response analysis of lightning wave process | |
| De Araújo et al. | Overvoltage analysis of transmission towers considering the influence of tower-footing impedance | |
| Asadpourahmadchali et al. | Hybrid continuous circuit‐trapezoidal integration method analysis of multi‐cross structure of grounding system | |
| Zhang et al. | Validation and revision of far-field-current relationship for the lightning strike to electrically short objects | |
| CN113255975B (zh) | 一种配电线路架设架空地线的空间位置优化方法 | |
| Zhang | Computation of lightning transients in large scale multiconductor systems | |
| Silveira et al. | On the lightning-induced voltage amplitude: First versus subsequent negative strokes | |
| Zhang et al. | Measurement and analysis of impulse grounding impedance of tower grounding device subjected to artificial lightning M‐component | |
| Yang et al. | Study on the lightning surge response of the AC side of UHVDC converter station based on finite difference time domain method | |
| CN113255976B (zh) | 一种配电线路架设架空地线的接地优化设计方法 | |
| Zhang et al. | A New FDTD Model for Lightning Return Stroke Channel Above Lossy Ground and Its Validation with Rocket-triggered Lightning Data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180420 Termination date: 20190429 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |