雷电冲击波在变压器绕组中的振荡过电压计算

雷电冲击波（如雷击）侵入变压器时，会在绕组内部引发复杂的电磁振荡过程，产生可能危及绝缘的过电压。准确计算这一振荡过电压，是变压器绝缘配合与设计的关键。本指南将手把手解析其核心原理与计算步骤。

第一步：理解物理过程与关键参数

计算前，必须先弄清现象背后的物理本质和决定性的参数。

1. 核心物理过程：雷电波（一个陡峭上升的电压波）从线路传入变压器绕组时，不能瞬间建立均匀的电压分布。绕组各匝、各线段间的电感和对地电容、匝间电容构成了一个复杂的分布参数网络。波前会在其中激发振荡，导致绕组某些点（如端部、中部）的电压远高于侵入波幅值。

2. 关键参数与简化模型：为便于工程计算，常将连续分布的绕组等效为集中参数的梯形网络或更简单的单绕组等效电路。你需要获取或计算以下参数：

   • U0：侵入绕组的雷电冲击电压幅值（单位：kV）。
   • t1：侵入波的波前时间（单位：µs），通常取标准雷电波 1.2µs。
   • L：绕组的总等值电感（单位：H）。它与绕组结构、匝数有关。
   • C：绕组的总等值对地电容（单位：F）。
   • K：绕组的等值纵向电容（或称为串联电容，单位：F）。它反映了匝间、饼间电容的影响，是决定初始电压分布均匀性和振荡频率的关键。
   • α：绕组系数（或称为分布系数）。这是一个无量纲数，计算公式为：
     $$ \alpha = \sqrt{\frac{C}{K}} $$
     它直接决定了初始电压分布的均匀性。α 值越大，初始电压分布越不均匀（首端梯度大），振荡越剧烈。

第二步：建立计算模型与公式

基于第一步的参数，我们采用经典的“单绕组等效振荡电路”模型进行计算。此模型将复杂的分布参数振荡，简化为由 L、C 构成的串联谐振回路来求解最大过电压。

1. 计算初始电压分布：雷电波侵入瞬间（t=0），由于电感电流不能突变，电压分布仅由电容链决定。距绕组首端 x 处（x 为归一化长度，首端为0，末端为1）的初始电压 U_i(x) 为：
   $$ U_i(x) = U_0 \cdot \frac{\sinh[\alpha(1-x)]}{\sinh(\alpha)} $$

   • 重点关注：绕组首端（x=0）的初始电压梯度最大，其值为：
     $$ \left. \frac{dU_i}{dx} \right|_{x=0} = -\alpha U_0 \coth(\alpha) \approx -\alpha U_0 \quad (\text{当 } \alpha > 3 \text{时}) $$
     这表明 α 越大，首端匝间承受的初始电压越高。

2. 计算稳态电压分布：振荡结束后（t→∞），电压按绕组电阻（直流）均匀分布。距首端 x 处的稳态电压 U_s(x) 为：
   $$ U_s(x) = U_0 \cdot (1 - x) $$

3. 计算自由振荡分量与最大过电压：初始分布与稳态分布之差，即为激发自由振荡的电压源。整个绕组可以看作一个 LC 振荡回路。其自然振荡角频率 ω0 为：
   $$ \omega_0 = \frac{1}{\sqrt{LC}} $$
   对应的振荡频率 f0 为：
   $$ f_0 = \frac{\omega_0}{2\pi} = \frac{1}{2\pi\sqrt{LC}} $$
   在振荡过程中，绕组任意点 x 的最大对地过电压 U_max(x) 可近似由下式估算：
   $$ U_{max}(x) \approx U_s(x) + [U_i(x) - U_s(x)] \cdot e^{-\frac{\delta}{2} \cdot \frac{\pi}{\omega_0}} $$
   其中 δ 为回路的衰减系数（与绕组电阻、铁芯损耗等有关）。更保守的工程估算常忽略衰减，取振荡幅值为初始差值的绝对值，则：
   $$ U_{max}(x) \approx U_s(x) + | U_i(x) - U_s(x) | $$
   最危险点通常出现在绕组中部或首端附近，需要将 x=0.5 等值代入上述公式计算。

4. 计算梯度放大系数：为评估匝间绝缘应力，需计算最大梯度放大系数 β。它定义为振荡过程中最大匝间电压与侵入波平均梯度 (U0 / N，N为总匝数) 的比值。简化公式为：
   $$ \beta \approx 1 + \frac{\alpha}{\pi} \cdot \frac{\omega_0 t_1}{1 + \omega_0 t_1} $$
   其中 t1 为波前时间。由此可得最大匝间电压 U_turn_max：
   $$ U_{turn\_max} = \beta \cdot \frac{U_0}{N} $$

第三步：执行计算流程（案例演示）

假设一台变压器绕组参数如下：

• U0 = 1000 kV
• t1 = 1.2 µs
• L = 0.5 H
• C = 5000 pF = 5 × 10^{-9} F
• K = 200 pF = 2 × 10^{-10} F
• 绕组总匝数 N = 1000

按步骤计算：

1. 计算绕组系数 α：
   $$ \alpha = \sqrt{\frac{C}{K}} = \sqrt{\frac{5 \times 10^{-9}}{2 \times 10^{-10}}} = \sqrt{25} = 5 $$
   α=5 > 3，属于初始分布极不均匀的绕组。

2. 计算自然振荡频率：
   $$ \omega_0 = \frac{1}{\sqrt{LC}} = \frac{1}{\sqrt{0.5 \times 5 \times 10^{-9}}} \approx \frac{1}{\sqrt{2.5 \times 10^{-9}}} \approx \frac{1}{5 \times 10^{-5}} \approx 20000 \text{ rad/s} $$
   $$ f_0 = \frac{\omega_0}{2\pi} \approx \frac{20000}{6.283} \approx 3183 \text{ Hz} $$

3. 计算绕组中部 (x=0.5) 的最大对地过电压：

   • 初始电压：U_i(0.5) = 1000 * sinh[5*(1-0.5)] / sinh(5) = 1000 * sinh(2.5) / sinh(5)。
     查表或计算得 sinh(2.5)≈6.05，sinh(5)≈74.2，故 U_i(0.5) ≈ 1000 * 6.05 / 74.2 ≈ 81.5 kV。
   • 稳态电压：U_s(0.5) = 1000 * (1 - 0.5) = 500 kV。
   • 初始与稳态差值：81.5 - 500 = -418.5 kV（绝对值 418.5 kV）。
   • 最大过电压（忽略衰减）：U_max(0.5) ≈ 500 + 418.5 = 918.5 kV。
     结论：绕组中部可能出现的最高对地电压约为 918.5 kV，低于侵入波幅值 1000 kV，但振荡过程电压可能瞬时更高。

4. 计算最大匝间电压（梯度放大系数法）：
   $$ \beta \approx 1 + \frac{\alpha}{\pi} \cdot \frac{\omega_0 t_1}{1 + \omega_0 t_1} = 1 + \frac{5}{3.1416} \cdot \frac{20000 \times 1.2 \times 10^{-6}}{1 + 20000 \times 1.2 \times 10^{-6}} $$
   先计算 ω0 t1 = 20000 * 1.2e-6 = 0.024。
   则 β ≈ 1 + (5/3.1416) * (0.024 / 1.024) ≈ 1 + 1.592 * 0.02344 ≈ 1.037。

   • 平均梯度：U0 / N = 1000 kV / 1000 = 1 kV/匝。
   • 最大匝间电压：U_turn_max = β * (U0/N) = 1.037 * 1 ≈ 1.037 kV。
     结论：最大匝间电压被放大至约 1.037 kV，需确保匝间绝缘能承受此电压。

第四步：影响因素分析与抑制措施

计算完成后，必须理解如何影响结果并采取对策。

1. 关键影响因素：

   • 波前时间 t1：波越陡（t1 越小），初始分布越不均匀，振荡越强，过电压越高。
   • 绕组系数 α：降低 α 是核心优化方向。α = sqrt(C/K)，因此需减小对地电容 C 或增大纵向电容 K。
   • 损耗：上述计算忽略了损耗（衰减系数 δ），实际铁芯损耗、电阻损耗会阻尼振荡，使实际过电压低于计算值。

2. 常用抑制措施：

   • 采用静电屏（电容环）：在绕组端部设置金属屏，人为增大首端区域的纵向电容 K，使 α 减小，从而均匀初始电压分布。
   • 使用纠结式或内屏蔽式绕组：通过复杂换位，大幅增加绕组自身的纵向电容 K，是降低 α、改善冲击分布最有效的方法。
   • 接入冲击吸收装置：在变压器入口安装金属氧化物避雷器 (MOA)，直接钳制侵入波幅值 U0。
   • 优化进波方式：通过三相进波比单相进波更能降低振荡过电压。

第五步：高级方法与流程总结

对于更精确的分析或特殊结构，需采用以下方法：

1. 多单元梯形网络模型：将绕组分割为数十个 LC 单元，利用电路仿真软件（如 EMTP、ATP-EMTP、Simulink）建立模型，输入标准雷电波形，直接获取各节点的时域电压波形。这是目前工程设计和研究中最主流的方法。

   graph TD A["输入: 标准雷电波 U(t)"] --> B[建立多段LC梯形网络模型] B --> C["设置模型参数\nL1...Ln, C1...Cn, K1...Kn"] C --> D[运行瞬态时域仿真] D --> E["输出: 各节点电压波形 U1(t)...Un(t)"] E --> F["分析波形\n提取峰值过电压与梯度"]

2. 完整工作流程总结：

   • 步骤1：获取 绕组结构参数，计算或测量得到 L, C, K。
   • 步骤2：计算 绕组系数 α 和自然频率 f0。
   • 步骤3：选择 计算方法（简化公式估算或建立仿真模型）。
   • 步骤4：执行 计算，得出关键点的最大对地过电压和匝间电压。
   • 步骤5：评估 结果是否在绝缘耐受范围内，若否，则设计 抑制措施（如增大 K），并返回步骤1进行迭代。