本篇为本科课程《高电压工程基础》的笔记。
本篇为这一单元的第二篇笔记。上一篇传送门。
行波通过串联电感与旁路并联电容
由于并联电容或串联电感的存在,线路上传播的行波会发生幅值和波形的变化。
直角波通过串联电感
有一个无限长的直角波 U 1 f U_{1f} U1f,从一个波阻抗为 Z 1 Z_1 Z1的导线上传来,经过过渡电感到一个波阻抗为 Z 2 Z_2 Z2的导线上来,如下图所示。
假设 Z 2 Z_2 Z2中反行波还没有到达节点,则其等效的电路图如下所示。
可根据电路图列出关系式为:
2 u 1 f = i 2 f ( Z 1 + Z 2 ) + L d i 2 f d t 2u_{1f}=i_{2f}(Z_1+Z_2)+L\frac{\mathrm{d}i_{2f}}{\mathrm{d}t} 2u1f=i2f(Z1+Z2)+Ldtdi2f
解方程可得解为:
i 2 f = 2 u 1 f Z 1 + Z 2 ( 1 − e − t T ) i_{2f}=\frac{2u_{1f}}{Z_1+Z_2}(1-e^{-\frac{t}{T}}) i2f=Z1+Z22u1f(1−e−Tt)
其中,T是这个回路的时间常数, T = L Z 1 + Z 2 T=\frac{L}{Z_1+Z_2} T=Z1+Z2L
根据电压电流的关系可得:
u 2 f = i 2 f Z 2 = 2 Z 2 Z 1 + Z 2 u 1 f ( 1 − e − t T ) u_{2f}=i_{2f}Z_2=\frac{2Z_2}{Z_1+Z_2}u_{1f}(1-e^{-\frac{t}{T}}) u2f=i2fZ2=Z1+Z22Z2u1f(1−e−Tt)
可见,前行波都是由两部分组成的:
- 第一部分是和时间无关的强制分量。
- 第二部分是随着时间而衰减的自由分量。
还可见,无穷长直角波通过集中的电感时,波头被拉长,而电压电流的稳态值和无电感时候一样。
沿着 Z 1 Z_1 Z1返回的反射波可列出方程:
u 2 f + L d i 2 f d t = u 1 f + u 1 b u_{2f}+L\frac{\mathrm{d}i_{2f}}{\mathrm{d}t}=u_{1f}+u_{1b} u2f+Ldtdi2f=u1f+u1b
将前行电压、电带波表达式入上式可得:
u 1 b = Z 1 − Z 2 Z 1 + Z 2 u 1 f + 2 Z 1 Z 1 + Z 2 u 1 f e − t T i 1 b = − u 1 b Z 1 = − Z 1 − Z 2 Z 1 + Z 2 u 1 f Z 1 − 2 Z 1 + Z 2 u 1 f e − t T u_{1b}=\frac{Z_1-Z_2}{Z_1+Z_2}u_{1f}+\frac{2Z_1}{Z_1+Z_2}u_{1f}e^{-\frac{t}{T}}\\\\ i_{1b}=-\frac{u_{1b}}{Z_1}=-\frac{Z_1-Z_2}{Z_1+Z_2}\frac{u_{1f}}{Z_1}-\frac{2}{Z_1+Z_2}u_{1f}e^{-\frac{t}{T}} u1b=Z1+Z2Z1−Z2u1f+Z1+Z22Z1u1fe−Tti1b=−Z1u1b=−Z1+Z2Z1−Z2Z1u1f−Z1+Z22u1fe−Tt
从上面两式可见,当t=0时, i 2 f = 0 , u 2 f = 0 , u 1 b = u 1 f , i 1 b = − i 1 f i_{2f}=0,u_{2f}=0,u_{1b}=u_{1f},i_{1b}=-i_{1f} i2f=0,u2f=0,u1b=u1f,i1b=−i1f,这是因为电感中的电流不能突变。当波到达电感的时刻,电感相当于开路,全部的磁场能量转换为电场能量,使得电压增大一倍。在这之后,按照指数规律变化。
当 t → ∞ t\rightarrow\infty t→∞时, i 2 f = 2 u 1 f Z 1 + Z 2 , u 2 f = 2 Z 2 Z 1 + Z 2 u 1 f , u 1 b = Z 1 − Z 2 Z 1 + Z 2 u 1 f , i 1 b = − Z 1 − Z 2 Z 1 + Z 2 u 1 f Z 1 i_{2f}=\frac{2u_{1f}}{Z_1+Z_2},u_{2f}=\frac{2Z_2}{Z_1+Z_2}u_{1f},u_{1b}=\frac{Z_1-Z_2}{Z_1+Z_2}u_{1f},i_{1b}=-\frac{Z_1-Z_2}{Z_1+Z_2}\frac{u_{1f}}{Z_1} i2f=Z1+Z22u1f,u2f=Z1+Z22Z2u1f,u1b=Z1+Z2Z1−Z2u1f,i1b=−Z1+Z2Z1−Z2Z1u1f。当时间趋于无穷的时候,电感相当于短路,所以就是两个导线的反射和折射作用。
电压和电流反射波随着时间的变化如下图所示:
折射电压波 u 2 f u_{2f} u2f的陡度表达式如下:
d u 2 f d t = 2 u 1 f L Z 2 e − t T \frac{\mathrm{d}u_{2f}}{\mathrm{d}t}=\frac{2u_{1f}}{L}Z_2e^{-\frac{t}{T}} dtdu2f=L2u1fZ2e−Tt
陡度的最大值在t=0处取得:
( d u 2 f d t ) m a x = d u 2 f d t ∣ t = 0 = 2 u 1 f Z 2 L \left(\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\right)_{\mathrm{max}}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\Big|_{t=0}=\frac{2u_{1f}Z_2}{L} (dtdu2f)max=dtdu2f t=0=L2u1fZ2
在空间上的陡度叫做空间陡度,而最大空间陡度尾:
( d u 2 f d l ) m a x = d u 2 f d l ∣ t = 0 = d u 2 f d t ∣ t = 0 d t d l = 2 u 1 f Z 2 L v \left(\frac{\mathrm{d}u_{2f}}{\mathrm{d}l}\right)_{\mathrm{max}}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}l}\Big|_{t=0}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\Big|_{t=0}\frac{\mathrm{d}t}{\mathrm{d}l}=\frac{2u_{1f}Z_2}{Lv} (dldu2f)max=dldu2f t=0=dtdu2f t=0dldt=Lv2u1fZ2
由上可见,想要降低 Z 2 Z_2 Z2上的前行电压波 u 2 f u_{2f} u2f陡度的有效措施是增加电感L。
直角波旁过并联电容
有一个无限长的直角波 u 1 f u_{1f} u1f,旁过电容接线方式如下图所示。
假设 Z 2 Z_2 Z2中反行波还没有到达节点A,则其等效的电路图如下所示。
可根据电路图列出关系式为:
2 u 1 f = i 1 Z 1 + i 2 f Z 2 i 1 = i 2 f + C d u 2 f d t = i 2 f + C Z 2 d i 2 f d t 2u_{1f}=i_1Z_1+i_{2f}Z_2\\\\ i_1=i_{2f}+C\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}=i_{2f}+CZ_2\frac{\mathrm{d}i_{2f}}{\mathrm{d}t} 2u1f=i1Z1+i2fZ2i1=i2f+Cdtdu2f=i2f+CZ2dtdi2f
其解为:
i 2 f = 2 u 1 f Z 1 + Z 2 ( 1 − e − t T ) u 2 f = Z 2 i 2 f = 2 Z 2 u 1 f Z 1 + Z 2 ( 1 − e − t T ) = α u 1 f ( 1 − e − t T ) i_{2f}=\frac{2u_{1f}}{Z_1+Z_2}(1-e^{-\frac{t}{T}})\\\\ u_{2f}=Z_2i_{2f}=\frac{2Z_2u_{1f}}{Z_1+Z_2}(1-e^{-\frac{t}{T}})=\alpha u_{1f}(1-e^{-\frac{t}{T}}) i2f=Z1+Z22u1f(1−e−Tt)u2f=Z2i2f=Z1+Z22Z2u1f(1−e−Tt)=αu1f(1−e−Tt)
其中,T是该电路的时间常数, T = Z 1 Z 2 Z 1 + Z 2 C T=\frac{Z_1Z_2}{Z_1+Z_2}C T=Z1+Z2Z1Z2C。 α \alpha α是电压波的折射系数。
由上面两式可得, u 2 f , i 2 f u_{2f},i_{2f} u2f,i2f均从零开始,按照指数规律趋向稳定值。原来的直角波变成指数波,波首变平,且稳态值只和波阻抗 Z 1 , Z 2 Z_1,Z_2 Z1,Z2有关,与电容值C无关。
因为电容的存在,在开始的时候电容两端的电压不能突变,随着充电的进行,电容充满了电,相当于开路,对于导线中的波传播不再起到作用。
因为有 u 1 f + u 1 b = u 2 f u_{1f}+u_{1b}=u_{2f} u1f+u1b=u2f,所以在导线 Z 1 Z_1 Z1上的反射波计算式为:
u 1 b = Z 2 − Z 1 Z 1 + Z 2 u 1 f − 2 Z 2 Z 1 + Z 2 u 1 f e − t T i 1 b = − u 1 b Z 1 = − Z 2 − Z 1 u 1 f Z 1 + Z 2 Z 1 + 2 Z 2 Z 1 + Z 2 u 1 f Z 1 e − t T u_{1b}=\frac{Z_2-Z_1}{Z_1+Z_2}u_{1f}-\frac{2Z_2}{Z_1+Z_2}u_{1f}e^{-\frac{t}{T}}\\\\ i_{1b}=-\frac{u_{1b}}{Z_1}=-\frac{Z_2-Z_1u_{1f}}{Z_1+Z_2Z_1}+\frac{2Z_2}{Z_1+Z_2}\frac{u_{1f}}{Z_1}e^{-\frac{t}{T}} u1b=Z1+Z2Z2−Z1u1f−Z1+Z22Z2u1fe−Tti1b=−Z1u1b=−Z1+Z2Z1Z2−Z1u1f+Z1+Z22Z2Z1u1fe−Tt
反射波的电压电流图如下所示:
当t=0, u 1 b = − u 1 f u_{1b}=-u_{1f} u1b=−u1f,因为电容电压不能突变,波到节点A的时刻,全部电场能量转换为磁场能量,相当于末端短路的反射。
在 Z 2 Z_2 Z2线路中折射电压的最大陡度可由表达式求导可得:
( d u 2 f d t ) m a x = d u 2 f d t ∣ t = 0 = 2 u 1 f Z 1 C \left(\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\right)_{\mathrm{max}}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\Big|_{t=0}=\frac{2u_{1f}}{Z_1C} (dtdu2f)max=dtdu2f t=0=Z1C2u1f
最大空间陡度为:
( d u 2 f d l ) m a x = d u 2 f d l ∣ t = 0 = d u 2 f d t ∣ t = 0 d t d l = 2 u 1 f Z 1 C v \left(\frac{\mathrm{d}u_{2f}}{\mathrm{d}l}\right)_{\mathrm{max}}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}l}\Big|_{t=0}=\frac{\mathrm{d}u_{2f}}{\mathrm{d}t}\Big|_{t=0}\frac{\mathrm{d}t}{\mathrm{d}l}=\frac{2u_{1f}}{Z_1Cv} (dldu2f)max=dldu2f t=0=dtdu2f t=0dldt=Z1Cv2u1f
可见最大空间陡度与 Z 2 Z_2 Z2无关,仅与 Z 1 , C Z_1,C Z1,C有关。
为了限制波的陡度,可以采用串联电感或并联电容,究竟选择哪一种,需要进行经济上的核算。
行波在无损平行多导线系统中的传播
前面分析是针对一根导线的波过程,实际中的输电线路都是由多根平行的导线构成的。比如装有避雷线的三相输电线路就有4根或5根的平行导线。此时波在导线上传播,平行多导线系统中就会发生相互的电磁耦合作用。
仍然是无损导线,传播的是平面电磁波,这就只需引入波速的概念,就可以将静电场中麦克斯韦方程应用于平行多导线系统。
如下图所示就是一个多导线系统以及导线的镜像。
图中有n根相互平行的导线,他们也和大地平行,各自带有的电荷量分别为: q 1 , q 2 , q 3 , . . . , q n q_1,q_ 2,q_3,...,q_n q1,q2,q3,...,qn,各自的电压大小为: u 1 , u 2 , . . . , u n u_1,u_2,...,u_n u1,u2,...,un,那么就可以写出下面的方程组:
u 1 = α 11 q 1 + α 12 q 2 + . . . + α 1 n q 1 u 2 = α 21 q 1 + α 22 q 2 + . . . + α 2 n q 1 . . . u n = α n 1 q 1 + α n 2 q 2 + . . . + α n n q 1 u_1=\alpha_{11}q_{1}+\alpha_{12}q_{2}+...+\alpha_{1n}q_{1}\\\\ u_2=\alpha_{21}q_{1}+\alpha_{22}q_{2}+...+\alpha_{2n}q_{1}\\\\ ...\\\\ u_n=\alpha_{n1}q_{1}+\alpha_{n2}q_{2}+...+\alpha_{nn}q_{1} u1=α11q1+α12q2+...+α1nq1u2=α21q1+α22q2+...+α2nq1...un=αn1q1+αn2q2+...+αnnq1
其中, α k k \alpha_{kk} αkk是自电位系数, α k j \alpha_{kj} αkj是互电位系数,他们的值取决于导线的几何尺寸和布置。计算公式如下:
α k k = 1 2 π ε r ε 0 ln 2 h k r k α k j = 1 2 π ε r ε 0 ln D k j d k j \alpha_{kk}=\frac{1}{2\pi\varepsilon_r\varepsilon_0}\ln\frac{2h_k}{r_k}\\\\ \alpha_{kj}=\frac{1}{2\pi\varepsilon_r\varepsilon_0}\ln\frac{D_{kj}}{d_{kj}} αkk=2πεrε01lnrk2hkαkj=2πεrε01lndkjDkj
其中, h k h_k hk为导线k距离地面的高度, r k r_k rk为导线k的半径, D k j D_{kj} Dkj是导线k和导线j的镜像之间的距离, d k j d_{kj} dkj是导线k和导线j之间的距离。
如果将电压和电荷的关系式中右边的电荷 q k q_k qk乘上 v v v就得到了电路的值 i k i_k ik,即 q k v = i k q_kv=i_k qkv=ik,而电位系数除以速度就是阻抗的量纲,那么就可以写出下面的公式:
u 1 = z 11 i 1 + z 12 i 2 + . . . + z 1 n i n u 2 = z 21 i 1 + z 22 i 2 + . . . + z 2 n i n . . . u n = z n 1 i 1 + z n 2 i 2 + . . . + z n n i n u_1=z_{11}i_1+z_{12}i_2+...+z_{1n}i_n\\\\ u_2=z_{21}i_1+z_{22}i_2+...+z_{2n}i_n\\\\ ...\\\\ u_n=z_{n1}i_1+z_{n2}i_2+...+z_{nn}i_n u1=z11i1+z12i2+...+z1ninu2=z21i1+z22i2+...+z2nin...un=zn1i1+zn2i2+...+znnin
其中, z k k z_{kk} zkk是导线k的自波阻抗, z k j z_{kj} zkj是导线k和导线j间的互波阻抗。对于架空线路有:
z k k = α k k C = 60 ln 2 h k r k z k j = z j k = α k j C = 60 ln D k j d k j z_{kk}=\frac{\alpha_{kk}}{C}=60\ln\frac{2h_k}{r_k}\\\\ z_{kj}=z_{jk}=\frac{\alpha_{kj}}{C}=60\ln\frac{D_{kj}}{d_{kj}} zkk=Cαkk=60lnrk2hkzkj=zjk=Cαkj=60lndkjDkj
如果一个导线上既存在前行波,也存在反行波,各自都可以列出电压和电流的关系式。那么对于有n根平行导线的系统,每一根导线k都可以列出下面的式子:
u k = u k f + u k b i k = i k f + i k b u k f = z k 1 i 1 f + z k 2 i 2 f + . . . + z k n i n f u k b = − ( z k 1 i 1 b + z k 2 i 2 b + . . . + z k n i n b ) u_k=u_{kf}+u_{kb}\\\\ i_k=i_{kf}+i_{kb}\\\\ u_{kf}=z_{k1}i_{1f}+z_{k2}i_{2f}+...+z_{kn}i_{nf} \\\\ u_{kb}=-(z_{k1}i_{1b}+z_{k2}i_{2b}+...+z_{kn}i_{nb}) uk=ukf+ukbik=ikf+ikbukf=zk1i1f+zk2i2f+...+zkninfukb=−(zk1i1b+zk2i2b+...+zkninb)
其中, u k f , u b f u_{kf},u_{bf} ukf,ubf分别是导线k上的前行电压波和反行电压波, i k f , i b f i_{kf},i_{bf} ikf,ibf分别是导线k上的前行电流波和反行电流波。上面的四个方程,前两个说明了前行波和反行波叠加就得到了导线上的波,第三个说明了前行波的电压和电流的关系,第四个说明了反行波的电压和电流的关系,而且此时需要加上负号。
n个导线可以列出n个方程然后求解,根据边界条件可以分析无损平行多导线系统中的波过程。但是多导线之间是存在电磁耦合,所以多导线系统的波过程不能孤立的看成彼此每一联系的单根导线的波过程,他们通常需要借助相模变换来求解。
两导线系统的耦合关系
有一个两导线系统,导线1为避雷线,导线2为对地绝缘的导线,如下图所示。
有一个雷击中了塔顶,避雷线上就有一个电压波 u 1 u_1 u1传播,对地绝缘的导线2上没有电流流过,但是他又处于导线1的电磁场之内,所以也会产生电压波 u 2 u_2 u2。根据电压和电流的关系可得:
u 1 = z 11 i i + z 12 i 2 u 2 = z 21 i i + z 22 i 2 u_1=z_{11}i_i+z_{12}i_2\\\\ u_2=z_{21}i_i+z_{22}i_2 u1=z11ii+z12i2u2=z21ii+z22i2
因为导线2电流为零,所以有:
u 2 = z 21 z 11 u 1 = K c 12 u 1 u_2=\frac{z_{21}}{z_{11}}u_1=K_{c12}u_1 u2=z11z21u1=Kc12u1
其中, K c 12 K_{c12} Kc12是导线1对导线2的耦合系数,因为存在关系 z 21 < z 11 z_{21}<z_{11} z21<z11,所以耦合系数 K c 12 < 1 K_{c12}<1 Kc12<1,大约为0.2到0.3,这是防雷的重要参数。
如下图所示,导线1对地有电场,而导线2获得了和 u 1 u_1 u1相同极性的对地电压 u 2 u_2 u2,这样就可以计算出导线1和导线2之间的电压差 Δ u \Delta u Δu:
Δ u = u 1 − u 2 = ( 1 − z 21 z 11 ) u 1 = ( 1 − K c 12 ) u 1 \Delta u=u_1-u_2=\left(1-\frac{z_{21}}{z_{11}}\right)u_1=(1-K_{c12})u_1 Δu=u1−u2=(1−z11z21)u1=(1−Kc12)u1
上面这个电压也是导线2对地的电压。
分析上式可知,当忽略耦合系数的时候,绝缘子串上承受的电压即 Δ u = u 1 \Delta u=u_1 Δu=u1,如果考虑耦合系数,那么绝缘子串上承受的电压即 Δ u = ( 1 − K c 12 ) u 1 \Delta u=(1-K_{c12})u_1 Δu=(1−Kc12)u1。可见耦合系数越大,电压差越小,越有利于防雷和绝缘子串的安全运行。在有些多雷地区,有时候会在导线下面假设耦合地线,以增大耦合系数。
线路杆塔的耦合系数
如下图所示是一个220kV线路杆塔。输电线路架设两根避雷线,他们之间通过金属塔杆相连接。
根据电压和电流的关系式,可得导线1、2、3之间的关系:
u 1 = z 11 i 1 + z 12 i 2 + z 13 i 3 u 2 = z 21 i 1 + z 22 i 2 + z 23 i 3 u 3 = z 31 i 1 + z 32 i 2 + z 33 i 3 u_1=z_{11}i_1+z_{12}i_2+z_{13}i_3\\\\ u_2=z_{21}i_1+z_{22}i_2+z_{23}i_3\\\\ u_3=z_{31}i_1+z_{32}i_2+z_{33}i_3 u1=z11i1+z12i2+z13i3u2=z21i1+z22i2+z23i3u3=z31i1+z32i2+z33i3
因为两根避雷线的离地高度和半径完全相等,是完全对称的,所以有 z 11 = z 22 , z 12 = z 21 , z 13 = z 31 , z 23 = z 32 , i 1 = i 2 z_{11}=z_{22},z_{12}=z_{21},z_{13}=z_{31},z_{23}=z_{32},i_1=i_2 z11=z22,z12=z21,z13=z31,z23=z32,i1=i2。
导线3是绝缘的,其上是无电流, i 3 = 0 i_3=0 i3=0,且避雷线上的电压 u 1 = u 2 = u u_1=u_2=u u1=u2=u,所以有:
u 1 = z 11 i 1 + z 12 i 2 u 2 = z 21 i 1 + z 22 i 2 u 3 = z 31 i 1 + z 32 i 2 u_1=z_{11}i_1+z_{12}i_2\\\\ u_2=z_{21}i_1+z_{22}i_2\\\\ u_3=z_{31}i_1+z_{32}i_2 u1=z11i1+z12i2u2=z21i1+z22i2u3=z31i1+z32i2
然后得到导线3的电压表达式:
u 3 = z 13 + z 23 z 11 + z 12 u = K c 1 , 2 − 3 u u_3=\frac{z_{13}+z_{23}}{z_{11}+z_{12}}u=K_{c1,2-3}u u3=z11+z12z13+z23u=Kc1,2−3u
其中, K c 1 , 2 − 3 K_{c1,2-3} Kc1,2−3是避雷线1,2对导线3的耦合系数,且 K c 1 , 2 − 3 = z 13 + z 23 z 11 + z 12 = z 13 / z 11 + z 23 / z 11 z 11 / z 11 + z 12 / z 11 = K c 13 + K c 23 1 + K c 12 K_{c1,2-3}=\frac{z_{13}+z_{23}}{z_{11}+z_{12}}=\frac{z_{13}/z_{11}+z_{23}/z_{11}}{z_{11}/z_{11}+z_{12}/z_{11}}=\frac{K_{c13}+K_{c23}}{1+K_{c12}} Kc1,2−3=z11+z12z13+z23=z11/z11+z12/z11z13/z11+z23/z11=1+Kc12Kc13+Kc23,而 K c 13 , K c 23 , K c 12 K_{c13},K_{c23},K_{c12} Kc13,Kc23,Kc12分别是导线1-3,2-3,1-2之间的耦合系数。
上面就是求得避雷线1,2对导线3的耦合系数的方法,而对于其他导线的耦合系数求取的方法是相同的。