【电力系统】基于Matlab的中压电缆的局部放电传输模型
✅作者简介：热爱科研的Matlab仿真开发者，擅长毕业设计辅导、数学建模、数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。
🍎 往期回顾关注个人主页：Matlab科研工作室
👇 关注我领取海量matlab电子书和数学建模资料
🍊个人信条：做科研，博学之、审问之、慎思之、明辨之、笃行之，是为：博学慎思，明辨笃行。
🔥 内容介绍
一、引言
中压电缆作为电力传输系统的关键组成部分，其运行可靠性直接关系到供电的稳定性和连续性。局部放电是中压电缆绝缘劣化的重要征兆，深入研究中压电缆的局部放电传输模型，对于准确检测和定位局部放电，及时评估电缆绝缘状态，预防电缆故障具有至关重要的意义。
二、中压电缆结构与局部放电概述
（一）中压电缆结构
中压电缆通常由导体、绝缘层、屏蔽层和护套等部分组成。导体作为电流传输的通道，一般采用铜或铝材质。绝缘层是电缆的核心部分，常用的绝缘材料有交联聚乙烯（XLPE）等，其作用是将导体与周围环境隔离，防止电流泄漏。屏蔽层分为内屏蔽层和外屏蔽层，内屏蔽层紧贴导体，外屏蔽层位于绝缘层外，它们能够均匀电场分布，减少电场集中现象。护套则起到保护电缆内部结构，防止机械损伤、水分侵入和化学腐蚀的作用。
（二）局部放电现象
局部放电是指在电缆绝缘内部，由于电场集中、绝缘缺陷等原因，在局部区域发生的放电现象，但放电未贯穿整个绝缘层。局部放电会产生电磁辐射、超声波、光、热等物理效应，长期的局部放电会逐渐侵蚀绝缘材料，导致绝缘性能下降，最终可能引发电缆绝缘击穿故障。例如，当绝缘层中存在气隙时，气隙内的电场强度会高于周围绝缘介质，从而使气隙内的气体发生电离，形成局部放电。
（二）电磁理论
局部放电会产生电磁辐射，根据电磁理论，电磁波在电缆中的传播受到电缆结构、绝缘材料电磁参数等因素的影响。例如，绝缘材料的相对介电常数和电导率会影响电磁波的传播速度和衰减。在分析局部放电传输模型时，需要考虑电磁理论来准确描述局部放电信号的电磁特性，如电场强度、磁场强度的分布以及它们随时间和空间的变化规律。
四、常见的局部放电传输模型
（一）频域传输模型
原理：频域传输模型基于傅里叶变换，将时域的局部放电信号转换到频域进行分析。通过研究局部放电信号在不同频率下的传输特性，如幅值衰减、相位变化等，来建立模型。在频域中，电缆的传输特性可以用传输函数来表示，传输函数描述了输入信号经过电缆传输后在不同频率下的输出响应。
应用：该模型适用于分析局部放电信号的频谱特性，对于研究局部放电信号在电缆中的高频分量传输和衰减情况具有重要意义。例如，通过频域传输模型可以确定局部放电信号中哪些频率成分在电缆传输过程中更容易衰减，从而为选择合适的检测频率提供依据。
（二）时域传输模型
原理：时域传输模型直接在时域内研究局部放电信号的传播过程。它通常基于传输线理论，通过数值方法求解传输线方程，得到局部放电信号在电缆中随时间和位置的变化情况。常用的数值方法有有限差分法（FDTD）、有限元法（FEM）等。这些方法将电缆结构离散化，对每个离散单元进行分析，从而模拟局部放电信号在整个电缆中的传播。
应用：时域传输模型能够直观地展示局部放电信号在电缆中的传播波形、传播速度以及反射、折射等现象。通过时域传输模型，可以准确地模拟局部放电信号在电缆中的传播过程，为局部放电的定位提供理论支持。例如，利用时域传输模型可以分析局部放电信号在电缆接头、终端等不连续点处的反射情况，从而实现对局部放电源的精确定位。
（三）混合传输模型
原理：混合传输模型结合了频域和时域传输模型的优点，在不同的分析阶段或针对不同的问题采用不同的方法。例如，在分析局部放电信号的初始传播和整体特性时，可以采用频域传输模型，快速获取信号的频谱特性和大致的传输衰减情况；而在对局部放电信号进行精确的定位和详细的波形分析时，则切换到时域传输模型，利用其直观、准确的特点进行深入研究。
应用：这种模型适用于复杂的中压电缆系统，能够综合考虑电缆的多种特性和局部放电信号的复杂传播情况。在实际工程中，中压电缆可能存在多种绝缘缺陷和不连续结构，混合传输模型可以更好地应对这些复杂情况，提高局部放电分析的准确性和可靠性。
五、模型的建立与验证
（一）模型建立步骤
电缆参数测量：首先需要测量中压电缆的各种参数，包括导体电阻、绝缘层的介电常数、电导率、电缆的几何尺寸等。这些参数是建立局部放电传输模型的基础，其准确性直接影响模型的精度。例如，通过高精度的测量仪器测量绝缘层的介电常数，为模型提供准确的电磁参数。
模型选择与参数设定：根据电缆的结构特点、研究目的以及局部放电信号的特性，选择合适的局部放电传输模型（频域、时域或混合模型）。然后，根据测量得到的电缆参数，设定模型中的相应参数，如传输线方程中的 R、L、G、C 等。
模型构建与求解：利用选定的模型和设定的参数，通过数学方法或数值计算软件构建局部放电传输模型，并进行求解。例如，采用 FDTD 方法求解时域传输模型时，需要将电缆结构离散化为网格，对每个网格单元进行电场和磁场的更新计算，从而模拟局部放电信号的传播过程。
⛳️ 运行结果
📣 部分代码
function out = optGaussPulseParam(in,t,w,t0,measPulse)
% OPTGAUSSPULSEPARAM - function for optimizing the gaussian pulse model parameters A and sigma
% based on measured pulse data -> objective function
% function call as function handle within a fminsearch simplex optimization
%
% Syntax: out = optGaussPulseParam(in,t,w,t0,measPulse)
%
% Inputs:
% in - double - vector of initial guess values of the gaussian model parameters
% t - double - sample time vector of corresponding measured pulse data in s
% w - double - angular frequency vector of corresponding measured pulse data in rad/s
% t0 - double - timeshift of the pulse event relative to time 0 in s
% measPul - double - measured data samples of one pulse
%
% Outputs:
% out - double - vector of optimized gaussian model parameters
%
% Other m-files required: none
% Subfunctions: gaussPulse
% MAT-files required: none
%
% Author: Martin Fritsch
% email: martin.fritsch@ovgu.de
% Mar 2020; Last revision: 19.02.2021
%------------- BEGIN CODE --------------
%% initial guess parameters from input - changed by the fminsearch in each iteration
A = in(1); % amplitude factor - determines the pulse height
sigma = in(2); % stretching factor - determines the pulse width in s
%% generate model pulse data with the initial guess
[GaussT,~] = gaussPulse(t,w,A,sigma,t0);
%% objective function - minimizing the difference between measured and model pulse data
out = sum((measPulse-GaussT).^2);
%------------- END OF CODE -------------
end