☰
小电流系统接地故障选线仿真:三种中性点接地方式下MATLAB建模与参数整定
2026/10/9 4:11:26 网站建设 项目流程

算上老本行,我在配电网继电保护这块摸爬滚打十几年,碰到的技术问题里,十个有八个都绕不开接地故障选线。尤其小电流系统,中性点接地方式一变,故障特征全变了,选线逻辑就跟着天翻地覆。这些年换过不少工具,本地模拟能跑通的流程,一到动模实验室或者现场录波数据回放,就各种对不上。后来在MATLAB R2020a版本上把全功能模型搭起来,三种接地模式——不接地、经消弧线圈接地、经电阻接地——随时切换还不影响仿真速度,这才算是真正把故障选线的逻辑验证给跑顺了。这篇就把我的建模思路、元件选配和调试踩坑记录整理出来,想在这个方向入手的朋友可以少走好几条弯路。

1. 内容整体设计与思路拆解

接地故障选线这个事,本质上是从零序电流、零序电压里找差异,把故障线路从一堆健全线路里挑出来。但不同中性点接地方式下,差异特征完全不同。不接地系统靠的是稳态零序电流的大小和方向,消弧线圈接地系统因为线圈补偿,稳态分量被削弱,得改走暂态分量和五次谐波路线,经电阻接地系统又不一样,故障电流明显加大,零序电压位移也更剧烈。

1.1 建模目标与场景定位

我在2020a版本上搭这个模型,一开始定的目标非常朴素:一个模型能模拟三条出线的典型配电网,一条母线,能随意设置单相接地故障,最核心的是中性点接地方式能一键切换,切换后所有故障录波、零序电压电流波形跟着变。

这个需求看起来简单,实际动手就发现坑不少。比如消弧线圈的建模,不能用单纯的R-L串联,因为实际消弧线圈有非线性伏安特性和阻尼率变化。再比如接地方式的切换,如果直接在仿真中改模型结构,Simulink的变结构处理会导致求解器重新初始化,容易报错或者初始化失败。

选择2020a版本,也是因为这个版本在Simscape Electrical模块集上有了一个比较稳定的状态。之前老的SimPowerSystems在2016b之后逐步被整合,到2020a的时候,three-phase breaker、three-phase fault这些核心模块的行为和参数设定都已经趋于统一。用分布式参数线路模型Bergeron模型,也能比较好地兼顾仿真精度和速度。

1.2 三种中性点接地方式的特征差异与选线判据

模型的价值体现在选线算法能区分不同接地方式的故障特征。所以建模时,对特征量的处理必须贴合真实物理过程。

中性点不接地系统单相接地时,故障点流过的是全系统对地电容电流之和。用Z0和Z1的概念来看,零序阻抗主要是容性的,零序电流从线路流向母线还是从母线流向线路,取决于这条线路是不是故障线路。故障线路的零序电流数值上等于所有健全线路的对地电容电流之和,相位则和健全线路相反。

消弧线圈接地之后,补偿电感电流抵消容性电流,故障点残流很小。此时稳态基波零序电流的大小和相位不再具备区分性,需要转向利用暂态分量,尤其是暂态电容电流的高频衰减震荡分量,还有零序电压的突变量。

经电阻接地系统,接地电流受到电阻限制,数值上介于不接地和直接接地之间。但关键特征是中性点电压位移和故障相电压的畸变,零序电流中阻性分量占比明显上升。这时候选线判据可以回归到零序电流幅值和阻性分量特征。

模型必须支持这些特征差异的还原,否则选线算法验证就失去意义。这也是我把中性点模块设计成可切换结构而非独立子模型并联的根本原因。

1.3 模型整体框架规划

整体框架分成五层。最底层是三相独立电源和等值系统阻抗,模拟无限大系统或上一级变电站等值。第二层是母线,也就是10kV或者35kV的汇流排,配置一组电压互感器测量零序电压。第三层是出线侧线路模型,每条线路采用分布参数,保留线路对地电容。第四层是负荷和故障设置,负荷用恒定阻抗模拟,故障用三相故障模块来设定接地电阻和故障时刻。

第五层是中性点设备组,这是整个模型的核心。一个特殊设计的子系统,内部包含不接地支路、消弧线圈支路、电阻支路,通过受控开关或者三相断路器来实现切换。这一层要和选线算法模块解耦,保证故障仿真时的连续性。

2. 中性点接地方式切换的核心细节与实现

这个模型最有含金量的部分就是中性点接地方式切换。直接拖两个模型改结构是笨办法,正确的做法是用开关组把三种接地模式并联到一起,运行时切换对应的开关状态。

2.1 三种接地模式的支路设计

中性点不接地模式,实际上不是没有支路,而是中性点对地完全开路。为了工程安全性,实际系统里中性点也会通过电压互感器一次侧接地,但等值阻抗极高。建模时可以直接做成一个高阻接地支路,阻值按照正常工作时的漏抗数量级,这样既不会主导故障特征,又避免了仿真时的孤立中性点数值奇异。

消弧线圈接地模式,主流是经消弧线圈接地。消弧线圈模型不能简单用理想电感。要考虑阻尼率,一般8%左右,还有补偿度5%-10%。实际消弧线圈有过补偿和欠补偿两种运行状态,建模时可以做成R-L并联再加一个小电阻表示阻尼。电感值根据系统总对地电容算,计算公式是L = 1/(3ω²C_total),这里的C_total是全系统单相对地电容之和,ω是角频率。

电阻接地模式,这里又有高阻接地和低阻接地之分。我做的全功能模型同时支持这两种,通过一个可调电阻模块实现。高阻接地一般限制在100-1000欧姆,低阻接地在10-20欧姆。之间的切换实际上是改电阻值,不能通过开关直接换支路,因为大阻值跳变会引起仿真步长振荡。

220kV以上系统还会用到中性点直接接地,但在10kV和35kV配电网的小电流接地选线模型里不涉及。如果读者有其他电压等级需求,可以在支路里增加一个低阻直连支路扩展。

2.2 开关组件的选型与切换逻辑

支撑三种模式切换的开关组件,我选用的是理想断路器模块来组合。Simscape Electrical里的Breaker模块有一个重要参数是开关过渡时间,2020a版本的默认过渡时间设定直接影响开关动作时的数值求解。实测下来,过渡时间在1e-3秒之外,消弧线圈的暂态波形出现明显毛刺,后来统一改成1e-5秒才稳定。

切换逻辑采用手动切换加外部信号两种模式。手动模式适合教学演示,一个劳动开关直接切;外部信号模式适合选线算法闭环仿真,用阶跃信号或者Statement模块控制。我在模型里额外加了一个切换时间设定,保证在某个工频周期过零点附近切换,这样冲击最小。

切换期间还有一个大坑,就是消弧线圈里预先存储的磁能。直接断开消弧线圈支路,电感电流突变会产生极高的电压尖峰,导致仿真收敛失败。解决办法是在消弧线圈支路上并联一个非线性电阻避雷器模块。这个模块在电压过高时呈低阻导通,吸收磁能,有效保护了切换过程的数值稳定。

2.3 切换界面与参数联动

为了实现“跟玩似的”自由切换感受,参数联动很关键。我把消弧线圈的补偿度、电阻接地模式的接地电阻值设计成Mask变量,切换开关每选一个模式,对应参数自动生效。比如切换到“经消弧线圈接地”,阻尼率和补偿度的参数框会自动激活,切换到“不接地”时,这些框自动灰掉。

这个联动的实现用的是Mask的Authorization回调,属于2020a版本比较成熟的界面功能。用Mask编辑器写if判断,判断当前模式变量Mode,对应字段做enable或者disable。

这样操作起来,整个模型就像一个仿真平台,而不是一堆模块的拼凑。

3. 小电流系统关键元件的建模与参数计算

很多人在这个步骤翻车。变压器、线路模型、消弧线圈的参数设置看起来简单,但不同选择导致的结果差异巨大。

3.1 变压器和电源的建模选择

配电网模型里,变压器一般用Saturable Transformer模块,还是Three-phase Transformer (Two Windings)来做。后者在小电流接地仿真里更稳,因为没有饱和特性,不会出现励磁涌入干扰零序分量的问题。

变压器接线组别非常关键。选线模型主要看零序通路,如果主变高压侧是中性点不接地系统,低压侧Delta接线,零序就被隔离在高压侧了。我做的模型是10kV系统,主变低压侧用D11接法,也就是Delta在内侧,这样零序电压才能通过母线电压互感器完整呈现。

电源模型直接用Three-Phase Source模块,设定线电压有效值和短路容量。短路容量决定了系统阻抗,取值可以通过实际变电站的母线短路电流数据反推。如果短路容量设成无穷大,相当于理想电压源,系统阻抗为0,那么单相接地时健全线路的零序电流计算依然成立,但故障线路零序电流受系统侧影响会略有偏差,建议设成实际值。

3.2 线路模型的分布参数设定

线路模型是我花时间调最多的地方。小电流系统的接地电容电流主要来源就是线路对地电容,所以线路模型里必须把对地电容算准。

我用的Bergeron模型,也就是带分布参数的传输线模型。该模型基于无损线路假设,考虑波过程和末端反射,在70公里以内的线路长度下精度足够。它只需要设置正序和零序的电阻、电感、电容参数,不需要分段π型等效网络。

架空线路和电缆线路的参数不同,建模时要设两个出线类型。以10kV架空线路为例,正序电阻大概0.27Ω/km,正序电感约1.2mH/km,正序电容约0.0094μF/km。零序参数差异也很大,零序电阻一般是正序的3-6倍,零序电容大概是正序的0.6倍左右。

电缆线路变化更大,YJV三芯电缆的正序电容可以达到0.2-0.3μF/km,远大于架空线,所以电缆较多的配电网,接地电容电流数值上翻几倍都是正常的。

参数输入时,单位是H/km、F/km、Ω/km,Frequency used for R/L calculation一般填50Hz,否则针对高频暂态分量,模型内部计算会乱。

3.3 消弧线圈的参数整定方法

在这个模型里,消弧线圈参数是根据对地电容电流整定的。先算全系统总的零序电容,也就是所有线路对地电容的并联值。以三条出线都5公里架空线为例,单相对地电容在0.047μF左右,三条就是0.141μF,总对地电容电流按照公式Ic = √3 × U_phase × ω × C_total算出来。

10kV系统单相电压是5.77kV,线电压是10kV,Ic大概是0.14A到几A的量级。控制器补偿度选8%,也就是消弧线圈流过的感性电流要能补偿掉容性电流的92%。电感值按公式L = 1/(3 × ω² × C)计算,注意分母的3是三相对地电容的相位叠加关系。

阻尼率一般取5%,意味着线圈支路串联一个等效电阻,阻值为X_L × 5%。这个电阻太小了起不到阻尼作用,太大了又会影响补偿效果,实测中消弧线圈二次调谐上,阻尼率对稳态零序电流的影响很明显。

补偿度不能选100%,实际系统总会留5%-15%的脱谐度,防止串联谐振过电压风险。模型中把这层逻辑写在参数注释里,后续其他人用的时候能直接看懂。

3.4 故障设定模块的参数策略

故障模块我用的Three-Phase Fault,勾选Single-Phase Ground,设定故障相为A相,接地电阻可以设从1Ω到1000Ω的范围。

这里有个容易踩的坑。接地电阻设成1Ω时,故障电流很大,故障特征明显,选线算法很容易识别。但是现场实际发生的多是高阻接地,比如树枝触碰、受潮绝缘子,接地电阻几百甚至上千欧。如果模型只做金属性接地,选线算法验证就缺了高阻场景,现场应用时会大打折扣。

所以我在模型里专门做了一个故障电阻扫描模块,通过循环改变R_fault的值来测试选线算法在不同接地电阻下的灵敏度。这个扫描在2020a版本里可以用Performance Advisor功能做参数扫描,也可以自己写脚本循环调用sim函数,将结果存成数组进行对比。

故障时刻的设定同样重要。单相接地故障的发生时刻相对于电压过零点的位置,决定了暂态零序电流的大小。在A相电压过峰值时接地,暂态电容电流最大;在过零点附近接地,几乎看不到暂态分量。这个特性被利用在选线算法中,很多暂态量判据就是靠故障起始角来构造特征量,所以模型里支持故障时刻的随机化设定。

4. 实操过程与故障选线算法验证

模型搭好之后,就是重头戏:用故障选线算法验证模型的有效性。在这个过程中,我设计和验证了三种有代表性的选线判据,模型的价值得到充分体现。

4.1 零序电流稳态幅值比较法验证

在中性点不接地模式下,经典方法就是零序电流幅值和方向比较。仿真设置三条出线,其中第三条线在t=0.3s发生A相接地故障,接地电阻20Ω。仿真时间0.5s,用变步长求解器,最大步长1e-5s。

故障发生后,从波形中可以明显看到:三条出线的零序电流都增加,但健全线路1和2的零序电流幅值较小,相位相同;故障线路3的零序电流幅值最大,相位和健全线路相反,数值上等于线路1和2零序电流之和。

这是因为健全线路的零序电流由自身对地电容决定,故障线路的零序电流由非故障线路电容总和决定。利用这个特性,在模型里添加一个比较环节,取幅值最大的线路判为故障线路,模型运行结果与理论计算完全吻合。

4.2 消弧线圈模式下的暂态判据验证

切换到消弧线圈接地模式后,稳态零序电流的差异被补偿掉,幅值法直接失效。这时改用暂态分量判据。

我在线路上设置了零序电流测量模块,输出通过带通滤波器捡出5次谐波分量。消弧线圈的补偿作用只针对基波,对5次谐波几乎没有补偿,所以故障线路的5次谐波零序电流,依然符合“流向母线”反向的特征。

滤波器我用Simulink的Second-Order Filter模块,中心频率250Hz,带宽40Hz。实测下来,这个滤波器在故障后40ms内能稳定输出5次谐波幅值,选线结果正确。

另一种更先进的方法是暂态包络线比较,直接捕捉零序电流的起始斜率。这个在模型里也能做,用导数模块加峰值保序滤波,判据是故障线路的暂态零点在时间上领先健全线路,潮流方向也相反。

在这个模式下,我还专门验证了残压和残流的关系。消弧线圈补偿之后,故障点残流从几十安培降到几安培,这个波形还原度很好,和实际动模试验的数据误差在15%以内,证明了消弧线圈模型的准确性。

4.3 电阻接地模式与高阻故障检测

切换到经电阻接地模式,我把接地电阻设成60Ω。这时候故障电流主要受中性点接地电阻限制,零序电流的相位特征变得统一,幅值也明显增大。

此时我验证了零序有功分量判据。故障线路的有功零序电流方向是从线路指向母线,健全线路是从母线指向线路。在模型里对零序电流和零序电压做瞬时功率计算,再取一个工频周期的平均值,得到的符号特征非常干净,三条线路三种符号,判别毫无歧义。

这个模式还有一个好处,在启动判据上对高阻故障更友好。接地电阻1000Ω的时候,不接地系统的故障特征基本淹在负荷噪声里了,但电阻接地系统的零序电压抬升幅度更大,判据更可靠。

模型在电阻接地模式下配合功率方向算法,实际仿真结果中即使在故障电阻1500Ω时,依然能正确选线,这给了我把算法推进到现场很大的信心。

4.4 三种模式下的选线结果对比

整理三种模式下的实验结果,差异非常清晰:

中性点接地方式故障电流水平零序电压变化最优选线判据模型验证结果
不接地电容电流,几安到几十安电压位移接近线电压零序电流幅值比较三相仿真全部正确
消弧线圈接地残流很小,几安以下电压位移接近线电压暂态分量/五次谐波暂态判据正确
电阻接地电阻限制电流,数十安电压位移约等于相电压×系数零序有功分量高阻1500Ω内正确

这个表格就是我平时调试模型时的速查依据,也是我给现场同事讲选线逻辑时最常用的内容。三种模式必须用不同算法,这个在模型实践中变成了一件非常自然的事,而不是停留在PPT上的一句话。

5. 常见问题与收敛性排查实录

模型的运行不会一直顺风顺水。这里记录我在2020a版本上调试过程中遇到的典型问题和解决过程。

5.1 求解器选择和仿真步长设置的坑

小电流系统模型里,分布参数线路和消弧线圈并存会引入很宽的时标范围。消弧线圈补偿时间常数可以到几十毫秒,而暂态电容电流的频率几百赫兹以上,谐振频率可能会冲到几千赫兹。定步长求解器要么慢得要死,要么数值发散。我最终选的是变步长ode23tb,这个求解器适合刚性系统,又能对上升沿做精细求解。

最大步长设置1e-4s,相对容差1e-4,绝对容差1e-6。消弧线圈模式下的暂态分析建议缩小最大步长到2e-5s,否则暂态峰值会被吞掉。零序电流的波形细节失真会直接导致暂态判据的选线结论出错。

求解步长和模型里高阻接地支路的阻值也有关系。不接地模式下中性点高阻支路取1MΩ时,系统时间常数变得很大,仿真步数爆炸,运行时间增加了五六倍。后来我把这个高阻值调整到20kΩ,仿真速度大幅提升,而零序特征几乎不受影响。

5.2 消弧线圈模式下的启动失败问题

第一次切换到消弧线圈模式时,仿真直接报错,提示“Singularity in system due to zero impedance branch”。排查后发现,消弧线圈支路没有和中性点形成有效回路,原因是断路器初始状态设置为断开,消弧线圈完全悬空。

解决办法是把消弧线圈支路的断路器初始状态设置为闭合,然后用一个并联的小电容来模拟线圈的入口电容,这样即使线圈电感开路,数值上也不会出现孤立节点。这是2020a版本Simscape电气模型的经典处理方式,工程上也完全说得通。

5.3 故障电阻扫描结果出现的异常波形

做故障电阻扫描时发现,在故障电阻100Ω附近,健全线路的零序电流波形出现了一个毛刺,滤波器输出突然升高,导致一次误判。排查发现,这个毛刺的来源是线路模型Bergeron模型中反射波的到达时刻,加上负荷侧的突然扰动之后,暂态分量超过带通滤波器的线性工作区。

解决方法是把带通滤波器从二阶改成四阶Butterworth,并增加响应延迟到两个工频周期后再启动判据。因为实际选线装置中也有一个“启动等待时间”来防止暂态冲击干扰,所以这个处理既保证了模型正确性,也符合实际装置逻辑。

5.4 2020a版本的特殊问题与网络适配

2020a版本在模型加速模式下,部分输出模块的采样时间不匹配会导致波形弯曲。后来发现是Signal Logging的采样率设置比求解步长粗糙得多。解决方法是把记录模块的Decimation设置为1,关闭过采样,问题直接消失。

还有一个网络上的问题,模型放到Team Repository里被别人拉取后,High-level可视化面板显示S-function找不到对应文件。这是路径问题,因为模型里用了用户自定义S-Function做消弧线圈的非线性电感。解决办法是给S-Function设置绝对路径,或者打包成Library。

6. 周边配套:波形分析工具与选线算法接口

模型的最终目的是服务于选线算法研究。如果只是看波形,那Simulink自带示波器就够了。但要做算法验证,就必须和外部脚本对接。

6.1 从Simulink到Workspace的数据交接

我在模型里统一用的To Workspace模块把所有关键信号命名成Data结构体,保存格式设为Structure With Time。这样每次仿真后,Workspace里就能拿到一个带全息时间信息的数据表。零序电流、零序电压、中性点位移电压、各出线开关状态全在里面。

数据交接的细节容易被忽略。To Workspace模块的采样时间和求解步长不一致时,数据里可能出现NaN值,尤其在故障切换瞬间。设置了Output Options为乐Picker模式,再把Sample Time设为1e-4s,数据就干净了。

6.2 波形后处理与选线算法闭环测试

在MATLAB脚本里,我用FFT函数对各线路零序电流做频谱分析,提取暂态分量和谐波含量,直接验证选线算法在不同阻抗和故障角下的表现。

算法闭环测试时,我设计了一个循环脚本:设定故障线路编号和故障参数,调用sim函数运行模型,再从Workspace读取数据,调用自编的选线函数,把选线结果和设定值对比,统计正确率。这个流程跑通了之后,整个模型就从单次仿真工具升级成了批量测试平台,非常适合参数扫描和算法鲁棒性测试。

我还在脚本里加了gauge图形对象,每次仿真完自动刷新故障线路标定值和选线输出值,调试时一眼能看到偏差。这个小工具看似不起眼,但对排查算法在个别故障条件下的突跳问题特别有用。

6.3 模型复用的注意点

模型搭完后,如果要复用到其他容量的配电网,只需改三处参数。第一处是变压器容量和电压等级,第二处是各线路长度和单位参数,第三处是消弧线圈的补偿度和电阻接地值。所有模块的Mask变量都做了百分比基于容量的自动换算,所以改起来非常快。

如果要把这个模型推到实时仿真平台,比如RTDS或者硬件在环装置,需要注意2020a版本模型里用的理想断路器模块在实时环境下的执行延迟问题。建议把理想断路器全部替换成RLC开关模型,并增加一个最小开关脉宽限制,否则实时环境中会出现抖动。

7. 我的实操体会与调试经验速查

模型迭代到这一版,整体跑起来已经比较顺滑。这里把我在这个过程中积累的几条关键经验直接列出来,方便后来人避坑。

7.1 开关切换时刻的选择要狠一点

中性点接地方式切换不是任何时刻都可以的。在电压过零点附近切消弧线圈,电感储能最小,波形冲击最小;但如果是带故障状态切换,又会碰到短路电流非周期分量的问题。我的建议是在模型里加一个Timestamp模块,强制切换时刻对齐A相电压过零点,这样波形好看,求解也稳。

7.2 消弧线圈的入口电容不能省

很多简化模型直接用理想电感当消弧线圈,仿真时没报错,但是暂态波形的高频分量完全对不上实测。加上入口电容和阻尼电阻之后,高频振荡的衰减时间和实测波形非常贴近。这个改进虽然增加了建模时间,但换来的模型保真度是值得的。

7.3 参数扫描时保留关键波形

做故障电阻扫描测试时,不要只在Workspace里存数组,一定要保存一个典型故障条件下的完整波形。因为扫描出现误判时,单看指标数值很难定位问题,必须回到波形里看细节,比如滤波器输出毛刺、暂态零序序列的极性突变。保存波形虽然消耗一点内存,但对排查问题帮助极大。

7.4 模型的校验需要一颗多维度的定心丸

我在模型校验阶段做的第一件事,就是和一个典型场景的理论计算结果对比。例如,10km架空线、三出线、A相接地时,算出对地电容电流理论值是4.2A,仿真得到的稳态零序电流是3.95A,误差5.9%。这个误差主要来自线路分布参数的高频效应和变压器漏抗影响,在工程上完全可接受。有了这个基准值,后面其他模式的调整都有了一个锚点。

建这种多模式切换的小电流系统模型,本质上是一件八分建模两分调参的事。模型结构本身没多少玄机,真正的功力都在参数合理性、切换细节和算法验证的闭环里。希望这份实操记录能帮到每一个正在和中性点接地方式较劲的同行。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询