ARTICLE DETAIL

资讯详情

深耕编程入门与网站建设的一线实战洞察。

小电流接地系统单相接地故障Simulink仿真建模与波形分析

小电流接地系统单相接地故障Simulink仿真建模与波形分析 写这篇博文之前我先说一句大实话小电流接地系统在10kV、35kV中压配电网里太常见了但很多刚接触继电保护的人对“单相接地之后到底发生了什么”其实是糊涂的。我早期做故障仿真的时候也踩过不少坑后来用Simulink把中性点不接地和经消弧线圈接地这两种系统完整搭了出来才把零序电压、零序电流、故障相电压恢复这些概念真正串起来。这篇文章就是围绕这套小电流接地系统故障仿真模型写的从整体设计思路到模块参数再到波形对比和调试经验一次性说清楚。适合正在做继电保护课程设计、毕业设计或者想搞懂配电网单相接地故障机理的研究生和现场工程师参考。1. 仿真模型搭建前的整体设计思路1.1 仿真到底想解决什么问题小电流接地系统单相接地故障最大的特点是故障电流很小线路可以带故障继续运行一段时间这对供电可靠性是有好处的。但副作用也很大非故障相电压会升高到线电压长时间运行可能击穿绝缘发展成两点接地短路同时接地故障不容易被检测出来选线困难是现场出了名的头疼事。仿真的目标就是要把这个“看不见摸不着”的过程变得可视化。我们希望通过模型看到三个核心现象第一单相接地后故障相电压怎么变化非故障相电压怎么升高第二零序电压和零序电流的幅值、相位、分布有什么特点第三中性点不接地系统里靠电容电流维持的接地电流换成消弧线圈补偿之后残流能压到多少。这三个问题搞清楚了后面做保护判据、选线算法才有基础。1.2 模型结构选型不搞花架子仿真要贴近真实工况很多初学者一上来就喜欢用理想电压源直接带负载这样跑出来的波形“干净”到不像话但和真实系统差太远。我的做法是模拟一座简化的10kV配电网络包含电源、变压器、两条出线、负荷和接地故障点。模型选型的几个关键考虑点我列一下电源用三相可编程电压源带内阻避免无穷大电源让故障特征失真。变压器用10kV/0.4kV Dyn11还是10kV中性点不接地取决于你要仿真的系统。这里我们仿真的是10kV中压侧所以变压器高压侧中性点要么悬空要么通过消弧线圈接地。线路用分布参数模型Three-Phase PI Section Line而不是简单的RLC串联。因为单相接地故障的暂态过程包含电容充放电集中阻抗根本模拟不出零序电容电流的分布。故障模块直接使用Simscape Electrical库里的Three-Phase Fault可以设置A相接地、B相接地或者两相接地也能调过渡电阻。这样的结构既不复杂又能真实反映单相接地故障时各序网络之间的关系。我给自己的标准是仿真波形能跟理论计算对上误差不超过工程允许范围。1.3 两种接地方式的本质区别一个靠系统电容一个靠电感补偿中性点不接地系统单相接地后故障点流过的电流是线路对地电容产生的电容电流这个电流大小由系统对地电容决定。系统越大、线路越长电容电流越大故障点电弧越难熄灭。当电容电流超过一定值比如10kV系统超过30A就需要把中性点改成经消弧线圈接地。消弧线圈的核心作用就是用一个可调电感产生的感性电流去抵消接地点的电容电流让通过故障点的残流变得很小。这个补偿关系是理解整套仿真的钥匙。我在模型里把消弧线圈单独做成一个可调电感支路就是为了方便观察不同补偿度下残流的变化。后面我会专门讲补偿度的计算。2. Simulink核心模块搭建与关键参数设置2.1 元件库盘点用到的模块和它们在哪儿在Simulink里做电力系统仿真核心库是Simscape Electrical老版本叫SimPowerSystems。我用的是MATLAB R2020a以上版本库路径叫Simscape Electrical Specialized Power Systems下面这些模块是这次仿真必须用到的模块名称库路径作用Three-Phase SourceElectrical Sources提供三相电压源可设内阻Three-Phase Transformer (Two Windings)Elements用于连接不同电压等级Three-Phase PI Section LineElements分布参数线路模型含正序零序参数Three-Phase Series RLC BranchElements搭消弧线圈回路Three-Phase FaultElements设置接地故障Three-Phase V-I MeasurementMeasurements测量三相电压电流Three-Phase BreakerElements模拟断路器分合闸powerguiSpecialized Power Systems求解器设置必须拖入模型Scope / To WorkspaceSimulink波形显示和数据导出powergui这个模块经常被忽略实际上是整套仿真的“地基”。它决定了用连续还是离散求解器也决定了你怎样观察稳态波形。如果模型里忘了放powergui仿真根本跑不起来。我一般用连续模式配合自动步长观察暂态过程比较方便如果模型过大跑太慢再切到离散模式。2.2 主电源与变压器参数给系统一个真实的内阻抗三相电压源的线电压设置为10kV频率50Hz。这里有一个很容易忽略的细节默认的三相电压源是理想电压源内阻为零这在接地故障仿真里会造成短路电流偏大、电压跌落过于“硬”。我习惯给电源叠加一个星形连接的RL阻抗等效短路容量法比如假设系统短路容量为200MVA那么系统阻抗Z U²/S 10.5² / 200 ≈ 0.55Ω再按X/R7把电阻和电感拆出来。这样波形会平滑很多也更接近现场。变压器在模型里承担两个职责一是把不同的电压等级连接起来二是提供一个中性点引出位置。仿真10kV系统时变压器高压侧要么不接地要么通过消弧线圈接地。如果你的仿真从10kV直接开始变压器可以省掉但要保证电源中性点悬空这个很容易被忽略导致系统变成中性点直接接地。2.3 线路模型参数正序和零序参数是仿真准确的关键线路参数一定要区分正序和零序这是继电保护仿真的基本功。我用的典型架空线路参数如下参数正序零序电阻 R (Ω/km)0.170.23电感 L (H/km)1.2e-33.6e-3电容 C (F/km)9.7e-96.0e-9电缆的参数差异很大比如YJV22-3×95型10kV电缆正序电容一般在0.2μF/km以上零序电容约为正序的2到3倍。为什么零序电容重要因为单相接地故障的电容电流主要就是零序电容电流它决定了故障特征幅度。对于中性点不接地系统接地电容电流估算公式是( I_C 3 \times \omega \times C_0 \times U_\phi )其中Uφ是相电压约5.77kVω100π。假设单条10km架空线每相零序电容C₀6nF/km总计60nF那么单相接地时对地电容电流约0.6A左右。如果换成10km电缆电容电流能到几十安培。这就是为什么电缆为主的配电网必须考虑消弧线圈或小电阻接地。在PI Section Line模块里正序和零序参数是分开填的千万别填错。2.4 消弧线圈回路建模电感、电阻和补偿度计算消弧线圈实际装置是有损耗的不能只填一个纯电感。我在模型里用Three-Phase Series RLC Branch搭了一条中性点对地支路串联一个约60mH到300mH的可调电感和一个约10Ω到20Ω的电阻。电阻用来模拟电抗器绕组损耗和阻尼电阻这个阻尼电阻非常重要它能把谐振过电压峰值压下来。补偿度的计算过程我举个例子。假设系统对地总电容电流为15A相电压为5.77kV那么系统的对地容抗是( X_C U_\phi / I_C 5770 / 15 \approx 384.7Ω )消弧线圈的感抗要接近这个值考虑到过补偿实际取值比容抗小10%左右( X_L 384.7 / 1.1 \approx 349.7Ω )对应电感( L X_L / (2πf) 349.7 / 314 \approx 1.11H )注意这个电感是折算到10kV系统的一相等效电感实际消弧线圈接在中性点三相零序回路等效为三个并联支路所以中性点处消弧线圈电感在零序网络中表现为3L。很多初学者直接把L填进模块里结果发现补偿效果对不上其实就是没有做这个折算。Zero-Order和正序的电感参数填法不一样调模型的时候一定要先算清楚。3. 两种接地方式仿真实现与结果对比3.1 中性点不接地系统单相接地后能带病运行但不该一直带病运行先跑最简单的中性点不接地系统。设置故障时间为0.1sA相在0.1s时通过10Ω过渡电阻接地持续0.2s后切除。仿真结果和理论完全对得上。A相电压在故障瞬间跌到接近零B相和C相电压升到线电压水平大约10kV。零序电压从接近零突变到约5.77kV也就是一倍的相电压。这说明一个重要的工程判据中性点不接地系统发生单相接地时零序电压会达到相电压级别这是最典型的启动判据。零序电流的分布更有意思。故障线路的零序电流幅值等于所有非故障线路的电容电流之和方向是从线路流向母线非故障线路的零序电流只等于自身电容电流方向是从母线流向线路。也就是说零序电流方向相反是选线的关键特征。这个“反相”现象在波形图上一眼就能看出来我在仿真里专门把两条线路的零序电流放到同一个Scope里对比非常直观。不过有一条要注意如果故障点过渡电阻太大比如经过1000Ω以上的高阻接地故障相电压不是跌到零非故障相电压升高也不明显整个特征会变得很钝。现场很多单相接地故障是高阻接地仿真时多扫几个过渡电阻值你会对保护的灵敏度有更直观的认识。3.2 经消弧线圈接地系统电感电流把电容电流“吃”掉了换成经消弧线圈接地后同样的故障条件重新仿真。这时消弧线圈提供的感性电流滞后于电容电流90度左右在故障点互相抵消。我用1.11H左右的中性点电感残流从原来的十几安培降到了几安培。波形上的另一个明显变化是熄弧后的故障相电压恢复过程。消弧线圈补偿电流会维持故障相电压接近零的特性消失取而代之的是一种缓慢恢复过程电弧过零熄灭后不容易重燃。这就是消弧线圈能有效降低接地故障危害的根本原因。但也带来一个麻烦故障线路和非故障线路的零序电流方向不再像不接地系统那样清晰反相特别是过补偿状态下故障线路零序电流方向变得和非故障线路一致。现场选线难本质就是这个原因。仿真里如果你只盯着稳态零序电流的幅值大小很容易判错。3.3 波形对比一个表格看懂两种系统的差异我在同一个模型里做了两组仿真把关键数据直接摘出来对比如下对比项中性点不接地经消弧线圈接地故障相电压降至近零降至近零但略高非故障相电压升至线电压约10kV同样升高到线电压附近零序电压幅值约5.77kV相电压约5.77kV略有差异故障点电流电容电流十几安补偿后残流几安故障线路零序电流方向与非故障线路相反可能相同过补偿时电弧重燃风险较高显著降低为什么零序电压两种情况下都接近相电压因为零序电压是由系统对地参数不对称决定的中性点接地方式改变的是故障电流通路并不会改变故障相电压跌落这个基本事实。所以零序电压通常被用作“系统有没有接地故障”的启动判据而零序电流的方向、大小才用来做选线和保护。3.4 为什么现场更愿意用消弧线圈但又不全靠消弧线圈工程上消弧线圈把接地电流压小故障电弧容易自行熄灭线路无需跳闸供电可靠性大大提高。所以在以架空线为主、电容电流不大的农村配电网里消弧线圈是主流。但如果系统电容电流特别大消弧线圈需要做得很大响应速度也可能跟不上暂态过程这时就会考虑改用中性点经小电阻接地让保护直接跳闸。到底选哪种方案要算清楚电容电流还要看供电可靠性要求。仿真模型里两种接地方式随意切换正好可以拿来对比分析给方案决策提供直观的数据支持。4. 常见问题与调试经验实录4.1 三相断路器当故障模块用是新手最常见的坑我见过太多人把Three-Phase Breaker当成故障设置模块结果怎么仿真都不出现接地效果。断路器断开的是主回路不是对地故障。要模拟单相接地故障必须用Three-Phase Fault模块在这个模块里勾选A相、B相或者C相并且把接地选项打开设置好故障电阻和故障时间。如果发现仿真时波形没有变化第一件事就检查故障模块的接地复选框是不是勾上了。另一个细节故障模块的时间参数用的是绝对时间还是相对于仿真启动时间默认是绝对时间。如果你的仿真起始时间不是0注意调整。我一般把故障起始设在0.1s让系统先稳定0.1s再故障这样能清楚看到故障前后的对比。4.2 波形毛刺多、仿真跑不动多半是求解器和阻尼参数的问题电力电子和电力系统仿真常见的毛病就是数值振荡。我在调消弧线圈模型时遇到过零序电压波形高频振荡特别严重的情况原因是故障瞬间电感电流突变连续求解器自动步长太小导致仿真速度极慢。我的处理方法是先切到powergui里的离散模式采样时间设为50μs或20μs这样仿真速度快很多波形也能接受。如果一定要用连续模式给消弧线圈并联一个很小的电阻比如1kΩ以上或者给故障模块并联一个snubber电路都能抑制高频振荡。snubber电阻默认一般是1MΩ对于10kV系统可能不够尝试调到10kΩ到100kΩ之间。4.3 零序分量测不准波形乱跳问题出在测量点很多初学者直接在三相电压测量模块后面接一个减法器来算零序电压结果得到一堆乱七八糟的波形。标准做法是用Three-Phase V-I Measurement测出三相电压电流后接一个零序计算模块。最简单的方式是用一个Gain模块把三相电压求和后乘以1/3得到的就是零序电压。零序电流同理对三相电流求和取1/3。不过要注意测量模块的电压测量参考点必须和系统接地参考一致。如果你把测量模块的地和系统接地分开测出来的零序量会莫名其妙。我在第一次搭模型的时候就是没注意这个接地参考零序电压一直不对排查了半天才发现是测量模块的接线问题。为了更容易观察也可以把三相电压直接接到一个Demux再用数学运算模块求零序。这个手动搭法虽然笨一点但对理解零序定义特别有帮助。4.4 消弧线圈参数明明算了补偿后残流还是大检查你的等效关系还有一种常见情况是按理论计算填入电感值后故障点电流完全没有减小。排查思路是这样的先确认你填的是不是中性点对地支路的电感值而不是折算前的每相电感值。前面说过零序网络中消弧线圈等效为3L如果你把计算得到的1.11H直接填进去实际等效感抗偏小补偿度就不对。其次检查消弧线圈电阻是不是太大。阻尼电阻用来抑制谐振过电压但太大会削弱补偿效果。我一般取感抗的3%到5%。比如感抗350Ω阻尼电阻取10Ω到15Ω就够。电阻太大残流会明显上升。4.5 高阻接地故障特征不明显怎么办如果你想模拟经过1000Ω甚至3000Ω高阻接地会发现非故障相电压升高不到线电压零序电压可能只有几百伏波形特征是“软绵绵”的。这其实是真实情况不是模型错误。此时需要通过提高测量精度、加长仿真时长、或者采用小波变换等暂态分析工具来提取特征。Simulink里可以把零序电压电流导出到工作区用MATLAB脚本做FFT分析提取五次谐波分量。实际故障选线装置很多时候就是靠这些暂态特征来识别高阻接地单纯看稳态量容易漏判。5. 从仿真到工程应用这些结论还能怎么用5.1 仿真结果如何与现场故障选线装置对标搭建这套模型不光是为了交作业更大的价值是验证选线算法。现场应用的故障选线装置常见判据有零序电流幅值法、零序电流方向法、五次谐波法、暂态分量法等。仿真模型里每条线路的零序电流都可以单独导出方便直接测试不同选线算法在各种过渡电阻、各种故障角下的表现。我在MATLAB里写了一个简单的脚本把故障后两个周波的零序电流数据导出来计算每条线路零序电流的幅值和相位然后根据“幅值最大且方向相反”的规则选线。在中性点不接地系统里这个规则很奏效换到消弧线圈接地系统后就得改用暂态分量或谐波分量。这个过程比单纯看波形有意思得多也更能理解现场选线为什么难。5.2 可以做哪些扩展从稳态仿真到暂态算法验证如果时间允许建议在这个模型基础上做三个方向的扩展。第一把固定电阻故障改成可控过渡电阻故障模拟电弧接地电阻的非线性变化市面上很多论文里的电弧接地模型就是这么加的。第二增加多条出线把两条改到五条、八条观察零序电流分布随线路数量的变化规律。第三把零序电压电流接入自己写的保护逻辑模块Simulink里可以用Stateflow或者MATLAB Function写启动判据搭一个简易的选线保护装置原型。还有一个实用技巧Simulink模型里所有需要对比波形的地方尽量用To Workspace把数据导出然后用MATLAB统一画图。Scope看波形方便但截图不够清晰而且不利于批量分析。数据导出一个矩阵用plot或者FFT工具统一后处理效率高得多。5.3 关于消弧线圈补偿度调整的一次实测心得最后分享一个我个人调试的小经验。我在模型里搭建消弧线圈时一开始贪心想同时对比欠补偿、全补偿、过补偿三种情况就把电感值设成了三个可切换的变量。但跑仿真时发现全补偿状态下系统出现了明显的谐振过电压零序电压飙得比故障相电压还高波形一直在振荡。这个现象在理论上是预料之中的但真正在波形上看到感觉还是很不一样。后来我把阻尼电阻加大谐振峰值就压下来了。这个经验告诉我消弧线圈补偿度的整定必须考虑阻尼电阻。实际装置里的自动调谐本质上就是在各种运行方式下找到一个既能补偿电容电流、又不至于引发谐振过电压的平衡点。仿真模型的好处就是可以大胆试错多跑几组参数你对这个“平衡点”的理解会比只看书本要深得多。这套模型我后来一直保留着做馈线自动化测试、做故障指示器算法验证的时候都会从这套模型里取数据。只要参数设置正确仿真波形和现场录波的特征是一致的用来做算法前期验证完全够用。新手如果用我这套方法搭一遍踩过的坑也基本能绕过去。
返回列表