配电网电缆线损工程估算简化迭代法及python实现

肖永威

456人浏览 · 2026-03-17 14:51:15

肖永威 · 2026-03-17 14:51:15 发布

1. 目的与适用范围

本方法用于快速估算地沟敷设电缆在给定负载电流、环境温度及谐波条件下的电能损耗（线损），并评估温度、负载和谐波对线损的影响。

适用场景：综合能源站充电桩供电回路，谐波含量较低（THD < 5%）。
方法特点：基于 IEC 60287 简化，采用迭代热平衡，参数取值结合工程经验，适合设计阶段和运营估算。

2. 理论基础

2.1. 导体电阻的温度效应

$Rdc(T)=R20×[1+α20(T−20)]R_{dc}(T) = R_{20} \times [1 + \alpha_{20} (T - 20)]$

$R_{20}$ ：20℃时导体直流电阻（Ω/km，查电缆手册）
$α20\alpha_{20}$ ：电阻温度系数，铜取 0.00393 /℃，铝取 0.00403 /℃
$T$ ：导体实际工作温度（℃）

2.2. 集肤效应与邻近效应（交流电阻）

交流电阻 $R_{ac}$ 为直流电阻乘以交流系数：
$R_{ac} = R_{dc} (1 + Y_s + Y_p)$
集肤效应因子 $Y_s$ 和邻近效应因子 $Y_p$ 按 IEC 60287-1-1 简化公式计算：

集肤效应因子
$xs2=8πfRdc×10−7ksx_s^2 = \frac{8\pi f}{R_{dc}} \times 10^{-7} k_s$
$Ys=xs4192+0.8xs4Y_s = \frac{x_s^4}{192 + 0.8 x_s^4}$

$f$ ：频率（Hz），若考虑谐波可用等效频率 $f_{eq}$ （见 2.4 节）
$R_{dc}$ ：直流电阻（Ω/m，注意单位）
$k_s$ ：导体结构系数，圆形绞合导体取 1

邻近效应因子（适用于三根单芯电缆呈等边三角形排列）
$xp2=8πfRdc×10−7kpx_p^2 = \frac{8\pi f}{R_{dc}} \times 10^{-7} k_p$
$Yp=xp4192+0.8xp4(ds)2×2.9Y_p = \frac{x_p^4}{192 + 0.8 x_p^4} \left( \frac{d}{s} \right)^2 \times 2.9$

$k_p$ ：系数，圆形绞合导体取 1
$d$ ：导体直径（mm）
$s$ ：导体中心间距（mm）

若电缆为多芯或排列不同，公式需调整，本简化版本采用三角形排列典型公式。

2.3. 热平衡与导体温度迭代

导体稳态温度由热平衡决定：
$Ploss=I2×Rac×nP_{loss} = I^2 \times R_{ac} \times n$
$Tcond=Tamb+Ploss×RthT_{cond} = T_{amb} + P_{loss} \times R_{th}$

$I$ ：负载电流有效值（A）
$n$ ：相数（三相取 3）
$T_{amb}$ ：地沟内环境温度（℃）
$R_{th}$ ：从导体到环境的总热阻（℃·m/W），包含绝缘、护层及外部环境热阻，需根据电缆结构和敷设条件估算。

迭代步骤：

假设初始导体温度 $T_{cond}^{(0)}$ （如 $T_{amb} + 10$ ）。
计算当前温度下的 $R_{dc}$ 、 $Y_s$ 、 $Y_p$ 、 $R_{ac}$ 。
计算损耗 $P_{loss}$ 及新温度 $T_{cond}^{(1)} = T_{amb} + P_{loss} R_{th}$ 。
重复 2~3 步直至 $∣Tcond(k)−Tcond(k−1)∣<ε|T_{cond}^{(k)} - T_{cond}^{(k-1)}| < \varepsilon$ （如 0.5℃）。

2.4. 谐波影响的简化处理

当谐波含量低于 5% 时，可近似忽略其对交流电阻的增量，仅按基波 50 Hz 计算。
若需粗略估计，可用等效频率法：
$feq=∑hhIh2∑Ih2⋅f1f_{eq} = \frac{\sum_{h} h I_h^2}{\sum I_h^2} \cdot f_1$
其中 $h$ 为谐波次数， $I_h$ 为各次谐波电流有效值。将 $f_{eq}$ 代入 $Y_s$ 、 $Y_p$ 公式。

3. 参数取值依据与建议

3.1. 参数取值依据表

参数	符号	典型值 / 依据	来源/备注
20℃直流电阻	$R_{20}$	查电缆样本（如 YJV-0.6/1kV，截面 95mm² 铜芯 $R_{20}=0.193 Ω/km$ ）	GB/T 3956
温度系数	$α20\alpha_{20}$	铜 0.00393，铝 0.00403	IEC 60287
导体直径	$d$	根据截面计算或查表（如 95mm² 铜导体直径约 11.0 mm）	经验公式 $\approx 1.13\sqrt{S}$
导体间距	$s$	对于品字形排列， $\approx d + 2\times$ 绝缘厚度；例如 YJV-0.6/1kV 绝缘厚度约 1.5mm，则 $\approx 11+3=14 mm$	电缆结构
频率	$f$	基波 50 Hz；若有谐波，用等效频率
集肤/邻近系数	$k_s, k_p$	圆形绞合导体取 1	IEC 60287
总热阻	$R_{th}$	需根据电缆结构和敷设环境估算，典型范围 0.8~3.5 ℃·m/W。	具体根据敷设条件查表3.2
		内部热阻 $T_1$ （绝缘）： $T1=ρT2πln⁡(1+2t1/d)T_1 = \frac{\rho_T}{2\pi} \ln(1+2t_1/d)$ ，XLPE 热阻系数 $ρT≈3.5\rho_T \approx 3.5$ K·m/W，t1 为绝缘厚度。	IEC 60287
		外部热阻 $T_4$ （地沟空气）：可简化为 $T4=1πDehT_4 = \frac{1}{\pi D_e h}$ ， $D_e$ 电缆外径，h 对流换热系数，静止空气取 5~10 W/(m²·K)。	工程手册
环境温度	$T_{amb}$	地沟内夏季最高可取 40℃，冬季最低 -20℃，依据当地气候。	设计规范
负载电流	$I$	充电桩实际运行电流（有效值）	运营数据

总热阻估算示例：
对于 YJV-0.6/1kV 3×95 电缆，导体直径 11.0 mm，绝缘厚度 1.5 mm，外护层厚度 2.0 mm，外径约 35 mm。

内部热阻 $T1=3.52πln⁡(1+2×1.5/11.0)≈0.14T_1 = \frac{3.5}{2\pi} \ln(1+2\times1.5/11.0) \approx 0.14$ ℃·m/W
外部热阻（对流）：设 h=8 W/(m²·K)， $T4=1π×0.035×8≈1.14T_4 = \frac{1}{\pi\times0.035\times8} \approx 1.14$ ℃·m/W
总热阻 $Rth=T1+T4≈1.28R_{th} = T_1 + T_4 \approx 1.28$ ℃·m/W

实际应用中，建议根据现场条件（地沟通风、电缆密集程度）调整 $T_4$ ，或参考类似项目实测数据。

3.2. 热阻取值方法建议

在电缆数量多、型号杂的工程现场，逐一获取每根电缆的绝缘厚度、外径、材料热阻系数确实不现实。针对地沟敷设场景，确实有更简洁的工程估算方法，核心思路是将复杂的热阻计算简化为基于敷设条件和经验的“热阻区间值”。

3.2.1. 热路合并法（半简化）

对于多根电缆，不需要每根单独计算内部热阻（T1、T2、T3），而是采用“等效电缆”概念：

将电缆群视为一个组合热源
散热路径简化为：电缆表面→地沟空气→沟壁→土壤
主要关注外部热阻T4，因为它占主导地位（通常占总热阻60%~80%）

3.2.2. 经验取值表法（最简洁）

根据电缆数量、敷设密度、地沟通风条件，直接查表取总热阻Rth的工程推荐值。这是目前工程设计中常用的方法。

地沟敷设电缆总热阻Rth经验取值表，下表基于IEC 60287和多个工程项目的实测数据综合给出：

敷设条件	电缆数量	通风情况	推荐Rth范围 (℃·m/W)	典型取值 (℃·m/W)	说明
稀疏敷设	≤3根	良好自然通风	0.8 ~ 1.2	1.0	电缆间距大，热量易散
中等密集	4~9根	一般通风	1.2 ~ 1.8	1.5	常见于综合能源站地沟
密集敷设	≥10根	通风较差	1.8 ~ 2.5	2.0	热量积聚明显
密闭地沟	任意	无通风	2.5 ~ 3.5	3.0	需特别注意温升

取值依据：

土壤热阻系数：一般取1.2 K·m/W（中等湿度土壤），干燥砂土可达2.5~3.0
电缆表面散热：静止空气中表面散热系数约5~10 W/(m²·K)
地沟尺寸影响：沟宽>0.8m、沟深>1.0m时，散热条件优于窄浅地沟

修正系数（可选）：
若已知土壤类型，可对上述取值进行修正：

潮湿粘土/沙土：×0.8
中等湿度土壤：×1.0
干燥砂土/碎石：×1.3~1.5

4. Python 实现代码

基于简化思路，让总热阻Rth由用户直接输入（根据上表选择），不再分解计算内部热阻。同时增加敏感性分析功能，便于评估温度、负载、热阻变化的影响。

"""
电缆线损工程估算简化迭代法（适用于地沟通风条件）
直接输入总热阻Rth，无需分解电缆参数
"""

import math
import matplotlib.pyplot as plt
import numpy as np

def cable_loss_simple(
    R20,                # 20℃直流电阻，Ω/m
    I,                  # 负载电流有效值，A
    Rth,                # 总热阻，℃·m/W（直接查表输入）
    Tamb=30.0,          # 环境温度，℃
    alpha=0.00393,      # 电阻温度系数，铜取0.00393
    harmonics=None,     # 等效频率
    f=50.0,             # 频率，Hz
    n_phases=3,         # 相数
    ks=1.0,             # 集肤系数
    kp=1.0,             # 邻近系数
    d_s_ratio=None,     # 可选：若需要估算集肤邻近效应，输入直径/间距比(d/s)
    max_iter=20,
    tol=0.5,            # 收敛判断中的 tol=0.1 可能过严，工程上 0.5℃ 已足够，但保留也无妨
    T_max=None          # 可选，绝缘允许温度，如 XLPE 取 90
):
    """
    简化迭代法计算线损

    参数:
        R20: 20℃直流电阻 (Ω/m)
        I: 负载电流有效值 (A) —— 若提供 harmonics，此 I 应为总电流有效值
        Rth: 总热阻 (℃·m/W) —— 根据敷设条件查表输入
        Tamb: 环境温度 (℃)
        alpha: 电阻温度系数 (/℃)
        harmonics: 可选，谐波列表，格式 [(次数, 电流有效值), ...]
        f: 基波频率 (Hz)，默认 50
        n_phases: 相数
        ks: 集肤效应系数，默认 1
        kp: 邻近效应系数，默认 1
        d_s_ratio: 可选，导体直径与间距比 (d/s)，用于估算邻近效应；若不提供则忽略邻近效应
        max_iter: 最大迭代次数
        tol: 收敛容差 (℃)

    返回:
        dict: 包含导体温度、交流电阻、线损等
    """
    # 若提供了谐波，计算等效频率
    if harmonics:
        I_total_sq = I**2
        f_eq_num = sum(h * Ih**2 for h, Ih in harmonics)
        f_eff = f_eq_num / I_total_sq * f
    else:
        f_eff = f    
    # 初始温度假设
    T_cond = Tamb + 10.0
    converged = False
        
    for _ in range(max_iter):
        # 1. 直流电阻温度修正
        R_dc = R20 * (1 + alpha * (T_cond - 20))
        
        # 2. 集肤效应因子
        xs2 = (8 * math.pi * f_eff / R_dc) * 1e-7 * ks
        Ys = xs2**2 / (192 + 0.8 * xs2**2) if xs2 > 0 else 0.0
        
        # 3. 邻近效应因子（如果提供了d/s比）
        Yp = 0.0
        if d_s_ratio:
            xp2 = (8 * math.pi * f_eff / R_dc) * 1e-7 * kp
            Yp = (xp2**2 / (192 + 0.8 * xp2**2)) * (d_s_ratio**2) * 2.9
        
        # 4. 交流电阻
        R_ac = R_dc * (1 + Ys + Yp)
        
        # 5. 总损耗
        P_loss = n_phases * I**2 * R_ac
        
        # 6. 新温度
        T_new = Tamb + P_loss * Rth
        
        if abs(T_new - T_cond) < tol:
            T_cond = T_new
            converged = True
            break
        T_cond = T_new
    
    # 重新计算最终值（基于收敛温度）
    R_dc_final = R20 * (1 + alpha * (T_cond - 20))
    xs2_final = (8 * math.pi * f_eff / R_dc_final) * 1e-7 * ks
    Ys_final = xs2_final**2 / (192 + 0.8 * xs2_final**2) if xs2_final > 0 else 0.0
    Yp_final = 0.0
    if d_s_ratio:
        xp2_final = (8 * math.pi * f_eff / R_dc_final) * 1e-7 * kp
        Yp_final = (xp2_final**2 / (192 + 0.8 * xp2_final**2)) * (d_s_ratio**2) * 2.9
    R_ac_final = R_dc_final * (1 + Ys_final + Yp_final)
    P_loss_final = n_phases * I**2 * R_ac_final

    # 可选温度警告
    warning = None
    if T_max and T_cond > T_max:
        warning = f"导体温度 {T_cond:.1f}℃ 超过允许值 {T_max}℃"

    return {
        'T_cond': round(T_cond, 1),
        'R_ac': round(R_ac_final, 8),
        'P_loss': round(P_loss_final, 3),      # W/m
        'loss_per_km': round(P_loss_final, 3),    # kW/km（数值相等）
        'converged': converged,
        'warning': warning
    }

def sensitivity_analysis(
    R20, I_base, Rth_base, Tamb_base,
    I_range=None, Tamb_range=None, Rth_range=None
):
    """
    敏感性分析：分别考察负载、环境温度、热阻对线损的影响
    """
    results = {}
    
    # 1. 负载影响
    if I_range:
        I_vals = I_range
        losses = []
        temps = []
        for I in I_vals:
            res = cable_loss_simple(R20, I, Rth_base, Tamb_base)
            losses.append(res['loss_per_km'])
            temps.append(res['T_cond'])
        results['I'] = {'vals': I_vals, 'loss': losses, 'temp': temps}
    
    # 2. 环境温度影响
    if Tamb_range:
        Tamb_vals = Tamb_range
        losses = []
        temps = []
        for Tamb in Tamb_vals:
            res = cable_loss_simple(R20, I_base, Rth_base, Tamb)
            losses.append(res['loss_per_km'])
            temps.append(res['T_cond'])
        results['Tamb'] = {'vals': Tamb_vals, 'loss': losses, 'temp': temps}
    
    # 3. 热阻影响
    if Rth_range:
        Rth_vals = Rth_range
        losses = []
        temps = []
        for Rth in Rth_vals:
            res = cable_loss_simple(R20, I_base, Rth, Tamb_base)
            losses.append(res['loss_per_km'])
            temps.append(res['T_cond'])
        results['Rth'] = {'vals': Rth_vals, 'loss': losses, 'temp': temps}
    
    return results


def plot_sensitivity(results):
    """绘制敏感性分析图"""
    plt.rcParams['font.sans-serif']=['SimHei'] #显示中文标签
    plt.rcParams['axes.unicode_minus']=False          
    fig, axes = plt.subplots(1, 3, figsize=(15, 4))
    
    # 负载影响
    if 'I' in results:
        ax = axes[0]
        ax2 = ax.twinx()
        ax.plot(results['I']['vals'], results['I']['loss'], 'b-o', label='线损')
        ax2.plot(results['I']['vals'], results['I']['temp'], 'r--s', label='导体温度')
        ax.set_xlabel('负载电流 (A)')
        ax.set_ylabel('线损 (kW/km)', color='b')
        ax2.set_ylabel('导体温度 (℃)', color='r')
        ax.set_title('负载电流影响')
    
    # 环境温度影响
    if 'Tamb' in results:
        ax = axes[1]
        ax2 = ax.twinx()
        ax.plot(results['Tamb']['vals'], results['Tamb']['loss'], 'b-o', label='线损')
        ax2.plot(results['Tamb']['vals'], results['Tamb']['temp'], 'r--s', label='导体温度')
        ax.set_xlabel('环境温度 (℃)')
        ax.set_ylabel('线损 (kW/km)', color='b')
        ax2.set_ylabel('导体温度 (℃)', color='r')
        ax.set_title('环境温度影响')
    
    # 热阻影响
    if 'Rth' in results:
        ax = axes[2]
        ax2 = ax.twinx()
        ax.plot(results['Rth']['vals'], results['Rth']['loss'], 'b-o', label='线损')
        ax2.plot(results['Rth']['vals'], results['Rth']['temp'], 'r--s', label='导体温度')
        ax.set_xlabel('总热阻 Rth (℃·m/W)')
        ax.set_ylabel('线损 (kW/km)', color='b')
        ax2.set_ylabel('导体温度 (℃)', color='r')
        ax.set_title('热阻影响')
    
    plt.tight_layout()
    plt.show()


# ==================== 使用示例 ====================

if __name__ == '__main__':
    # 电缆参数：YJV-0.6/1kV 3×95 铜芯
    R20 = 0.193e-3  # 0.193 Ω/km → Ω/m
    
    print("="*60)
    print("电缆线损估算（简化迭代法）")
    print("="*60)
    
    # 场景1：稀疏敷设，通风良好
    print("\n【场景1】稀疏敷设，地沟通风良好")
    print("-"*40)
    res1 = cable_loss_simple(
        R20=R20,
        I=200.0,
        Rth=1.0,        # 查表：稀疏敷设典型值
        Tamb=30.0,
        d_s_ratio=11/14  # 可选：估算邻近效应
    )
    print(f"导体温度: {res1['T_cond']} ℃")
    print(f"交流电阻: {res1['R_ac']:.8f} Ω/m")
    print(f"线损: {res1['P_loss']:.3f} W/m ({res1['loss_per_km']:.3f} kW/km)")
    
    # 场景2：中等密集，常规地沟
    print("\n【场景2】中等密集，常规地沟")
    print("-"*40)
    res2 = cable_loss_simple(
        R20=R20,
        I=200.0,
        Rth=1.5,        # 查表：中等密集典型值
        Tamb=30.0,
        d_s_ratio=11/14
    )
    print(f"导体温度: {res2['T_cond']} ℃")
    print(f"交流电阻: {res2['R_ac']:.8f} Ω/m")
    print(f"线损: {res2['P_loss']:.3f} W/m ({res2['loss_per_km']:.3f} kW/km)")
    
    # 场景3：密集敷设，通风差
    print("\n【场景3】密集敷设，通风较差")
    print("-"*40)
    res3 = cable_loss_simple(
        R20=R20,
        I=200.0,
        Rth=2.0,        # 查表：密集敷设典型值
        Tamb=30.0,
        d_s_ratio=11/14
    )
    print(f"导体温度: {res3['T_cond']} ℃")
    print(f"交流电阻: {res3['R_ac']:.8f} Ω/m")
    print(f"线损: {res3['P_loss']:.3f} W/m ({res3['loss_per_km']:.3f} kW/km)")
    
    # 敏感性分析：评估温度、负载、热阻变化的影响
    print("\n" + "="*60)
    print("敏感性分析（用于设计评估）")
    print("="*60)
    
    # 基准工况：中等密集，负载200A
    results = sensitivity_analysis(
        R20=R20,
        I_base=200.0,
        Rth_base=1.5,
        Tamb_base=30.0,
        I_range=[100, 150, 200, 250, 300],
        Tamb_range=[-20, -10, 0, 10, 20, 30, 40],
        Rth_range=[0.8, 1.0, 1.2, 1.5, 1.8, 2.0, 2.5]
    )
    
    # 输出关键结论
    print("\n【负载影响】")
    for i, loss in zip(results['I']['vals'], results['I']['loss']):
        print(f"  {i:3d} A → 线损 {loss:.3f} kW/km")
    
    print("\n【环境温度影响】")
    for t, loss in zip(results['Tamb']['vals'], results['Tamb']['loss']):
        print(f"  {t:3d} ℃ → 线损 {loss:.3f} kW/km")
    
    print("\n【热阻影响】")
    for r, loss in zip(results['Rth']['vals'], results['Rth']['loss']):
        print(f"  Rth={r:.1f} → 线损 {loss:.3f} kW/km")
    
    # 绘制敏感性分析图（如果需要）
    plot_sensitivity(results)

在这里插入图片描述

5. 使用说明

参数准备：

从电缆样本获取 $R_{20}$ （注意单位换算为 Ω/m）、导体直径 $d$ 、导体间距 $s$ （可估算）。
根据地沟通风条件估算总热阻 $R_{th}$ ：无强制通风且电缆密集时取较大值（如 1.5），通风良好取较小值（如 0.8）。
若考虑谐波，提供 harmonics 列表，格式 [(次数, 电流有效值)]，函数将自动计算等效频率。注意：此时 I 应为总电流有效值（包含基波和谐波）。

结果解读：
- T_cond：导体稳态温度，可用于判断是否超过绝缘允许温度（如 XLPE 允许 90℃）。
- R_ac：单位长度交流电阻，用于计算线损。
- P_loss：三相总损耗（W/m），乘以电缆长度即为总损耗功率。
敏感性分析：可改变 Tamb、I、Rth 等参数，观察对线损的影响。