一种基于Lasso理论的接地网腐蚀诊断方法与流程

一种基于lasso理论的接地网腐蚀诊断方法
技术领域
1.本发明涉及变电站接地网故障诊断技术领域，特别涉及一种基于lasso理论的接地网腐蚀诊断方法。


背景技术：

2.变电站接地网具有均压和泄流功能，对保证电力系统的安全稳定运行，保护变电站工作人员的人身安全和电气设备的运行安全起到重要作用。变电站接地网大多采用普通钢或镀锌扁钢，容易受到土壤腐蚀和电化学腐蚀的影响而导致接地体锈蚀变薄甚至断裂脱落。接地网腐蚀会降低接地网的接地性能，抬升地网的接地电位，危及设备及人身安全，给变电站的运行带来巨大的安全隐患。
3.目前接地网腐蚀诊断主要有三种方法：电磁法、电化学法及电网络法。电磁法向接地网的引出线注入交流电流，通过检测地面上磁感应强度的分布情况来进行接地网腐蚀的诊断与定位，电磁法会存在电磁干扰问题，特别是变电站这种强电磁环境会对测量结果产生较大影响；电化学法利用电化学特征来表征接地网的腐蚀情况，接地网在电化学腐蚀的作用下后接地极会发生氧化反应变成阳极，通过电化学传感器测量接地网不同位置的极化电阻来诊断接地网的腐蚀程度，电化学法需要开挖土壤埋设传感器，且周围的噪声会对检测精度产生影响；电网络法利用电网络理论和特勒根定理建立故障诊断方程组，通过端口电阻的变化情况得到接地网各支路电阻阻值的变化量，通过求解端口网络阻抗变化的方程组，得到接地网各支路的支路电阻，根据支路电阻的阻值得到接地网各支路的腐蚀情况。这种方法操作简单，便于实现，但由于接地网引出线的数目少于支路数，采用电网络法测量所得的端口电压电流方程数少于待求解的支路数，因而电网络法所建立的方程为欠定方程，目前，对这类方程一般采用优化方法进行求解，如文献[许磊,李琳.基于电网络理论的变电站接地网腐蚀及断点诊断方法[j].电工技术学报，2012,27(10):270-276]利用最小二乘法进行求解，这类方法在可及节点过少(方程数过少)的情况下求解会出现过拟合、诊断结果误差大、出现伪故障等问题。


技术实现要素：

[0004]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，利用lasso理论对支路电阻变化量进行稀疏处理，再对故障诊断方程组进行求解，以解决变电站接地网腐蚀诊断中采用电网络方法时欠定方程求解容易出现不收敛及出现伪故障的问题。
[0005]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，包括：
[0006]
步骤1：建立变电站接地网的结构关系矩阵，同时，获取变电站接地网的基本参数；
[0007]
步骤2：对变电站接地网进行检测，获得检测参数；
[0008]
步骤3：基于所述结构关系矩阵、基本参数和所述检测参数，建立变电站接地网腐蚀诊断方程组；
[0009]
步骤4：利用lasso理论迭代求解所述变电站接地网腐蚀诊断方程组，获得求解结
果；
[0010]
步骤5：基于求解结果，图形化展示变电站接地网的腐蚀程度。
[0011]
优选的，所述步骤1中，建立变电站接地网的结构关系矩阵，包括：
[0012]
分析变电站接地网的设计图纸的拓扑结构，建立描述变电站接地网的结构关系的关联矩阵a；
[0013]
其中，关联矩阵a描述变电站接地网的节点与支路之间的连接关系，a的“行”代表接地网的节点元素，a的“列”代表接地网支路元素；按照从左到右、从上到下作为参考方向，对于a的每个元素，当该元素对应的支路按照参考方向离开该元素所对应的节点时，该元素取+1；当该元素对应的支路按照参考方向进入该元素所对应的节点时，该元素取-1；当该元素对应的支路与该元素对应的节点不相关联时，该元素取0。
[0014]
优选的，所述步骤1中，获取变电站接地网的基本参数，包括：
[0015]
根据变电站接地网的支路的长度和支路采用的接地材料，计算支路的标称电阻值ri＝ρli/si，其中，i＝1,2,3...,b,b为接地网支路总数，ρ为地网支路金属的电阻率，li为地网支路的长度，si为地网支路金属的截面积，并据此得到接地网的支路标称电阻向量rk＝(r1,r2,...ri,...,rb)；
[0016]
由支路标称电阻值rk得接地网支路电导矩阵进一步得接地网节点电导矩阵gn＝ag
bat
，其中，下标n代表接地网节点总数。
[0017]
优选的，所述步骤2：对变电站接地网进行检测，获得检测参数，包括：
[0018]
将变电站地网上设有引出线的节点作为可及节点，选定变电站地网中任一节点作为参考节点，在参考节点与可及节点1之间注入直流电流i0，测量可及节点1与参考节点之间的端口电压值u
′
(1)
；
[0019]
选择下一个可及节点与参考节点之间注入直流电流i0，测量该可及节点与参考节点之间的端口电压值u
′
(2)
；
[0020]
再对其余可及节点依次重复上述操作，记录其对应的端口电压值u
′
(m)
，m为可及节点的个数，得到注入电流i0情况下的端口电压向量u’＝(u
′
(1)
,u
′
(2)
,
…
,u
′
(m)
)。
[0021]
优选的，所述步骤3：基于所述结构关系矩阵、基本参数和所述检测参数，建立变电站接地网腐蚀诊断方程组，包括：
[0022]
(31)以步骤2中选取的参考节点，根据变电站接地网的网络拓扑和步骤1中计算所得的接地网节点电导矩阵gn＝ag
bat
，计算出在步骤2中各个可及节点注入直流电流i0情况下，变电站接地网各节点电压理论值：下，变电站接地网各节点电压理论值：其中in是注入节点n电流向量，in＝(0
(1)
,0
(2)
,...,i
0(i)
,...,-i
0(j)
,...,0
(n)
)；
[0023]
(32)计算接地网支路电流理论值,i
k(1)
＝g
batun(1)
，i
k(2)
＝g
batun(2)
，...，i
k(m)
＝g
batun(m)
，并据此得到接地网支路电流向量理论值为ik＝(i
k(1)
,i
k(2)
,...,i
k(m)
)；
[0024]
(33)根据所述接地网各节点电压理论值，计算在各可及节点注入直流电流i0时，各可及节点与参考节点之间的端口电压理论值u＝(u
(1)
,u
(2)
,
……u(m)
)；
[0025]
(34)根据欧姆定律计算接地网各可及节点端口电阻的理论值：r
ij(1)
＝u
(1)
/i0，r
ij(2)
＝u
(2)
/i0，...，r
ij(m)
＝u
(m)
/i0，并据此得到接地网端口电阻理论值向量r
ij
＝(r
ij(1)
,r
ij(2)
,...,r
ij(m)
)；
[0026]
(35)根据步骤2得到接地网检测参数，计算接地网可及节点与参考节点端口电阻的实际值r
′
ij(1)
＝u
′
(1)
/i0，r
′
ij(2)
＝u
′
(2)
/i0，...，r
′
ij(m)
＝u
′
(m)
/i0，并据此得到接地网端口电阻实际值向量r
′
ij
＝(r
′
ij(1)
,r
′
ij(2)
,...,r
′
ij(m)
)；
[0027]
(36)计算接地网端口电阻实际值与理论值的变化量：
△rij
＝r
′
ij-r
ij
；
[0028]
(37)通过特勒根定理建立接地网支路电阻变化量与端口电阻变化量之间的关系：
[0029][0030]
其中，
△rij
是(36)中计算的接地网端口电阻实际值与理论值的变化量，
△rk
是待求的接地网支路电阻变化量，ik是(32)中计算所得的各支路电流理论值，i
′k是腐蚀后各支路电流的实际值，与支路电阻的阻值变化量
△rk
有关，i0是可及节点和参考节点之间注入的直流电流值。
[0031]
优选的，所述步骤4：利用lasso理论迭代求解所述变电站接地网腐蚀诊断方程组，获得求解结果，包括：
[0032]
(41)设置腐蚀后各支路电流的初始值，令于是步骤三)的接地网支路电阻方程变为
[0033]
(42)随机化初始给
△rk
随机取一个初值
[0034]
(43)从开始，到为止，依次求
[0035]
的解如下：
[0036]
其中，||||1是l1正则化，λ∈[0,1]为惩罚系数，可人为指定；
[0037]
在求解时，均为常数；
[0038]
依次求出
[0039]
(44)由上次计算所得的更新并重新计算并重新计算由此得到下一次迭代时的支路电流值
[0040]
(45)将更新后的代入方程
[0041]
(46)判断(其中ε为设置的误差常数)是否成立，若不成立，则返回(43)进行循环；若成立，则为接地网支路电阻变化量与端口电阻变化量方程的最优解；
[0042]
(47)计算接地网支路电阻实际值的解
[0043]
优选的，所述步骤5：基于求解结果，图形化展示变电站接地网的腐蚀程度，包括：
[0044]
根据步骤4求解所得接地网支路电阻的实际值，计算支路电阻的变化倍数pk＝r'k/rk，其中，k＝1,2,...,b，b代表支路数；
[0045]
基于预设的评价标准，根据支路电阻变化倍数pk＝r'k/rk对接地网各支路的腐蚀程度进行评价，获得评价结果；
[0046]
图形化展示所述评价结果，其中，轻度腐蚀的支路用绿色表示，中度腐蚀的支路用黄色表示，严重腐蚀的支路用红色表示；
[0047]
根据接地网各支路的腐蚀程度绘制接地网实际各支路腐蚀诊断结果图，并进行输出。
[0048]
优选的，基于lasso理论的接地网腐蚀诊断方法，还包括：
[0049]
步骤6：训练腐蚀应对模型，基于所述腐蚀应对模型，根据变电站接地网的腐蚀程度，制定应对策略，基于所述应对策略，进行相应应对处理。
[0050]
优选的，所述步骤6中，训练腐蚀应对模型，包括：
[0051]
获取多个第一腐蚀应对事件，同时，获取所述第一腐蚀应对事件对应的应对方；
[0052]
获取所述应对方对应的应对类型，所述应对类型包括：内部应对和外部应对；
[0053]
当所述应对方对应的应对类型为内部应对时，获取所述应对方对应的经验程度值；
[0054]
若所述经验程度值小于等于预设的第一阈值，剔除对应所述第一腐蚀应对事件；
[0055]
当所述应对方对应的应对类型为外部应对时，获取所述应对方对应的信用程度值；
[0056]
若所述信用程度值小于等于预设的第二阈值，剔除对应所述第一腐蚀应对事件；
[0057]
当需要剔除的所述第一腐蚀应对事件均剔除后，将剔除剩余的所述第一腐蚀应对事件作为第二腐蚀应对事件；
[0058]
对所述第二腐蚀应对事件进行合理性分析，获取合理值；
[0059]
若所述合理值大于等于预设的第三阈值，将对应所述第二腐蚀应对事件作为第三腐蚀应对事件；
[0060]
基于预设的模型训练算法，根据所述第三腐蚀应对事件进行模型训练，获得腐蚀应对模型。
[0061]
优选的，对所述第二腐蚀应对事件进行合理性分析，获取合理值，包括：
[0062]
获取多个专家节点，同时，获取所述专家节点对所述第二腐蚀应对事件进行合理性评价的评价值；
[0063]
获取所述专家节点对应的专家权重，赋予所述评价值对应所述专家权重，获得目标值；
[0064]
累加计算所述目标值，获得合理值。
[0065]
有益效果：
[0066]
(1)本发明利用lasso理论进行特征压缩，不仅可以解决过拟合问题，而且可以将不重要的变量直接压缩为0，删除无效变量(未被腐蚀的支路)的影响。可以在参数缩减过程中，达到提取有效的支路电阻变化量
△rk
(已被腐蚀的支路)的作用。既实现了较为准确的
接地网故障诊断，也实现了变量选择(降维)。
[0067]
(2)本发明解决了接地网腐蚀诊断的盲目性，无需大面积开挖土壤即可较为准确的进行故障诊断与定位，为接地网是否开挖进行更换以及开挖位置的选择提供参考，进而降低现场接地网的运营维护成本。
[0068]
(3)本发明通过图形化显示，使接地网腐蚀诊断结果形象直观，应用方便。
[0069]
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
[0070]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0071]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
[0072]
图1为本发明实施例中一种基于lasso理论的接地网腐蚀诊断方法的流程图；
[0073]
图2为本发明实施例中变电站接地网网络拓扑结构的示意图；
[0074]
图3为本发明实施例中变电站接地网腐蚀诊断结果的示意图；
[0075]
图4为本发明实施例中又一基于lasso理论的接地网腐蚀诊断方法的流程图。
具体实施方式
[0076]
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
[0077]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，如图1-3所示，包括：
[0078]
步骤1：建立变电站接地网的结构关系矩阵，同时，获取变电站接地网的基本参数；
[0079]
步骤2：对变电站接地网进行检测，获得检测参数；
[0080]
步骤3：基于所述结构关系矩阵、基本参数和所述检测参数，建立变电站接地网腐蚀诊断方程组；
[0081]
步骤4：利用lasso理论迭代求解所述变电站接地网腐蚀诊断方程组，获得求解结果；
[0082]
步骤5：基于求解结果，图形化展示变电站接地网的腐蚀程度；
[0083]
所述步骤1中，建立变电站接地网的结构关系矩阵，包括：
[0084]
分析变电站接地网的设计图纸的拓扑结构，建立描述变电站接地网的结构关系的关联矩阵a；
[0085]
其中，关联矩阵a描述变电站接地网的节点与支路之间的连接关系，a的“行”代表接地网的节点元素，a的“列”代表接地网支路元素；按照从左到右、从上到下作为参考方向，对于a的每个元素，当该元素对应的支路按照参考方向离开该元素所对应的节点时，该元素取+1；当该元素对应的支路按照参考方向进入该元素所对应的节点时，该元素取-1；当该元素对应的支路与该元素对应的节点不相关联时，该元素取0；
[0086]
所述步骤1中，获取变电站接地网的基本参数，包括：
[0087]
根据变电站接地网的支路的长度和支路采用的接地材料，计算支路的标称电阻值ri
＝ρli/si，其中，i＝1,2,3...,b,b为接地网支路总数，ρ为地网支路金属的电阻率，li为地网支路的长度，si为地网支路金属的截面积，并据此得到接地网的支路标称电阻向量rk＝(r1,r2,...ri,...,rb)；
[0088]
由支路标称电阻值rk得接地网支路电导矩阵进一步得接地网节点电导矩阵gn＝ag
bat
，其中，下标n代表接地网节点总数；
[0089]
所述步骤2：对变电站接地网进行检测，获得检测参数，包括：
[0090]
将变电站地网上设有引出线的节点作为可及节点，选定变电站地网中任一节点作为参考节点，在参考节点与可及节点1之间注入直流电流i0，测量可及节点1与参考节点之间的端口电压值u
′
(1)
；
[0091]
选择下一个可及节点与参考节点之间注入直流电流i0，测量该可及节点与参考节点之间的端口电压值u
′
(2)
；
[0092]
再对其余可及节点依次重复上述操作，记录其对应的端口电压值u
′
(m)
，m为可及节点的个数，得到注入电流i0情况下的端口电压向量u’＝(u
′
(1)
,u
′
(2)
,
…
,u
′
(m)
)；
[0093]
所述步骤3：基于所述结构关系矩阵、基本参数和所述检测参数，建立变电站接地网腐蚀诊断方程组，包括：
[0094]
(31)以步骤2中选取的参考节点，根据变电站接地网的网络拓扑和步骤1中计算所得的接地网节点电导矩阵gn＝ag
bat
，计算出在步骤2中各个可及节点注入直流电流i0情况下，变电站接地网各节点电压理论值：下，变电站接地网各节点电压理论值：其中in是注入节点n电流向量，in＝(0
(1)
,0
(2)
,...,i
0(i)
,...,-i
0(j)
,...,0
(n)
)；
[0095]
(32)计算接地网支路电流理论值,i
k(1)
＝g
batun(1)
，i
k(2)
＝g
batun(2)
，...，i
k(m)
＝g
batun(m)
，并据此得到接地网支路电流向量理论值为ik＝(i
k(1)
,i
k(2)
,...,i
k(m)
)；
[0096]
(33)根据所述接地网各节点电压理论值，计算在各可及节点注入直流电流i0时，各可及节点与参考节点之间的端口电压理论值u＝(u
(1)
,u
(2)
,
……u(m)
)；
[0097]
(34)根据欧姆定律计算接地网各可及节点端口电阻的理论值：r
ij(1)
＝u
(1)
/i0，r
ij(2)
＝u
(2)
/i0，...，r
ij(m)
＝u
(m)
/i0，并据此得到接地网端口电阻理论值向量r
ij
＝(r
ij(1)
,r
ij(2)
,...,r
ij(m)
)；
[0098]
(35)根据步骤2得到接地网检测参数，计算接地网可及节点与参考节点端口电阻的实际值r
′
ij(1)
＝u
′
(1)
/i0，r
′
ij(2)
＝u
′
(2)
/i0，...，r
′
ij(m)
＝u
′
(m)
/i0，并据此得到接地网端口电阻实际值向量r
′
ij
＝(r
′
ij(1)
,r
′
ij(2)
,...,r
′
ij(m)
)；
[0099]
(36)计算接地网端口电阻实际值与理论值的变化量：
△rij
＝r
′
ij-r
ij
；
[0100]
(37)通过特勒根定理建立接地网支路电阻变化量与端口电阻变化量之间的关系：
[0101][0102]
其中，
△rij
是(36)中计算的接地网端口电阻实际值与理论值的变化量，
△rk
是待求的接地网支路电阻变化量，ik是(32)中计算所得的各支路电流理论值，i'k是腐蚀后各支路电流的实际值，与支路电阻的阻值变化量
△rk
有关，i0是可及节点和参考节点之间注入的直流电流值；
[0103]
所述步骤4：利用lasso理论迭代求解所述变电站接地网腐蚀诊断方程组，获得求
解结果，包括：
[0104]
(41)设置腐蚀后各支路电流的初始值，令于是步骤三)的接地网支路电阻方程变为
[0105]
(42)随机化初始给
△rk
随机取一个初值
[0106]
(43)从开始，到为止，依次求
[0107]
的解如下：
[0108]
其中，||||1是l1正则化，λ∈[0,1]为惩罚系数，可人为指定；
[0109]
在求解时，均为常数；
[0110]
依次求出
[0111]
(44)由上次计算所得的更新并重新计算并重新计算由此得到下一次迭代时的支路电流值
[0112]
(45)将更新后的代入方程
[0113]
(46)判断(其中ε为设置的误差常数)是否成立，若不成立，则返回(43)进行循环；若成立，则为接地网支路电阻变化量与端口电阻变化量方程的最优解；
[0114]
(47)计算接地网支路电阻实际值的解
[0115]
所述步骤5：基于求解结果，图形化展示变电站接地网的腐蚀程度，包括：
[0116]
根据步骤4求解所得接地网支路电阻的实际值，计算支路电阻的变化倍数pk＝r'k/rk，其中，k＝1,2,...,b，b代表支路数；
[0117]
基于预设的评价标准，根据支路电阻变化倍数pk＝r'k/rk对接地网各支路的腐蚀程度进行评价，获得评价结果；
[0118]
图形化展示所述评价结果，其中，轻度腐蚀的支路用绿色表示，中度腐蚀的支路用黄色表示，严重腐蚀的支路用红色表示；
[0119]
根据接地网各支路的腐蚀程度绘制接地网实际各支路腐蚀诊断结果图，并进行输出。
[0120]
上述技术方案的工作原理及有益效果为：
[0121]
实施例的接地网网络拓扑及节点、支路编号如图1所示，矩形代表接地网支路，圆
形代表接地网节点，红色圆形代表接地网可及节点，实施例中假设支路4与支路8出现了腐蚀，其支路电阻为10欧姆，其他支路状况良好，支路电阻为1欧姆，可及节点为3、6、9、15节点，参考节点为16节点。
[0122]
(1)建立变电站接地网结构关系矩阵，获取接地网基本参数。根据变电站接地网的设计图纸，分析其拓扑结构，建立描述接地网结构的关联矩阵a。关联矩阵a描述接地网节点与支路之间的连接关系，a的“行”代表接地网的节点元素，a的“列”代表接地网支路元素。按照从左到右、从上到下作为参考方向，对于a的每个元素，当该元素对应的支路按照参考方向离开该元素所对应的节点时，该元素取+1；当该元素对应的支路按照参考方向进入该元素所对应的节点时，该元素取-1；当该元素对应的支路与该元素对应的节点不相关联时，该元素取0。
[0123]
以4x4节点网络数据为例，计算出该网络的关联矩阵a，该关联矩阵大小为16x24，也即该网络共有16个节点，24条支路。
[0124][0125][0126]
(2)根据接地网的基本情况，得到接地网的支路电阻向量rk＝(r1,r2,...ri,...,r
24
)＝(1,1,...,1,...,1)。
[0127]
(3)对接地网的实际情况进行仿真模拟，获取实际接地网的运行参数。利用matlab建立仿真模型，设置16号节点为参考节点，选定3、6、9、15节点作为可及节点。各可及节点与参考节点之间注入的直流电流源i0设置为10a，仿真得到可及节点3,6，9,15分别注入电流i0时的端口电压向量u’＝(u
′
(3)
,u
′
(6)
,u
′
(9)
,u
′
(15)
)＝(14.7743，14.1934，17.0429，7.0356)。
[0128]
(4)正常情况(腐蚀前)下支路电导矩阵为也即主对角线为rk的倒数，其余都为0的24x24方阵。由gn＝ag
bat
得到16x16的节点电导矩阵，再由得到16x4的节点电压矩阵。其中节点注入电流矩阵in为：
[0129][0130]
(5)得到正常情况(腐蚀前)下支路电流矩阵为ik＝g
bat
un，是24x4矩阵。接地网支路电流矩阵理论值为ik＝(i
k(3)
,i
k(6)
,i
k(9)
,i
k(15)
)。
[0131][0132]
(6)计算在各可及节点注入直流电流i0时，腐蚀前各可及节点与参考节点之间的端口电压理论值。u＝(u
(3)
,u
(6)
,u
(9)
,u
(15)
)＝(14.5982，12.9464，14.5982，7.0089)。
[0133]
(7)根据得到接地网端口电压检测参数，计算接地网可及节点与参考节点端口电阻的实际值r
′
ij(3)
＝u
′
(3)
/i0，r
′
ij(6)
＝u
′
(6)
/i0，r
′
ij(9)
＝u
′
(9)
/i0，r
′
ij(15)
＝u
′
(15)
/i0，得到接地网端口电阻实际值向量r
′
ij
＝(r
′
ij(3)
,r
′
ij(6)
,r
′
ij(9)
,r
′
ij(15)
)＝(1.47743，1.41934，1.70429，0.70356)。
[0134]
(8)根据欧姆定律计算接地网各可及节点端口电阻的理论值：r
ij(3)
＝u
(3)
/i0，r
ij(6)
＝u
(6)
/i0，r
ij(9)
＝u
(9)
/i0，r
ij(15)
＝u
(15)
/i0，得到接地网端口电阻理论值向量r
ij
＝(r
ij(3)
,rij(6)
,r
ij(9)
,r
ij(15)
)＝(1.45982，1.29464，1.45982，0.70089)。
[0135]
(9)计算接地网端口电阻实际值与理论值的变化量：
△rij
＝r
′
ij-r
ij
＝(
△rij(3)
,
△rij(6)
,
△rij(9)
,
△rij(15)
)＝(0.0176，0.1247，0.2445，0.0027)。
[0136]
(10)初始迭代计算时，将腐蚀后支路电流i'k等于腐蚀前支路电流ik，也即令再将
△rij
、ik、i0代入接地网故障诊断方程组中，得到四个欠定方程。
[0137][0138]
在matlab中调用lasso函数，使得最小，得到初次解
[0139]
得到初次lasso迭代计算的支路变化量为：(0，0，0，0.6369，0，0，0，2.1800，0，0，0，0，0，0，0，0，0，0，0，0，0，0，0，0)。
[0140]
由可以得到为：(1，1，1，1.6369，1，1，1，3.1800，1，1，1，1，1，
1，1，1，1，1，1，1，1，1，1，1)。
[0141]
(11)由上次计算所得的更新并重新计算由此得到下一次迭代时的支路电流值
[0142]
(12)将更新后的代入方程利用lasso优化方法求解使min的
[0143]
(13)判断(其中ε为设置的误差常数)是否成立，其中ε设置为0.01。若不成立，则返回(11)进行循环；若成立，则为接地网支路电阻变化量与端口电阻变化量方程的最优解。计算接地网支路电阻实际值的解
[0144]
本例中，经12次迭代后，满足收敛条件得到方程最终解为：
[0145]
△rk
＝(0，0，0，7.2231，0，0，0，7.8930，0，0，0，0，0，0，0，0，0，0，0，0，0，0，0，0)。
[0146]
根据计算实际接地网各支路电阻阻值的计算值为：
[0147][0148]
(14)根据求解所得接地网支路电阻的实际值，计算支路电阻的变化倍数pk＝r'k/rk(k＝1,2,...,b。b代表支路数)。
[0149]
pk＝r'k/rk＝(1，1，1，8.2231，1，1，1，8.8930，1，1，1，1，1，1，1，1，1，1，1，1，1，1，1，1)。
[0150]
(15)根据支路电阻变化倍数pk＝r'k/rk对接地网各支路的腐蚀程度进行评价，预设的评价标准如下表所示：
[0151][0152]
如图3所示。将接地网支路腐蚀程度图形化显示，轻度腐蚀的支路用绿色表示，中
度腐蚀的支路用黄色表示，严重腐蚀的支路用红色表示。并根据接地网各支路的腐蚀程度绘制接地网实际各支路腐蚀诊断结果图。
[0153]
本发明利用lasso理论进行特征压缩，不仅可以解决过拟合问题，而且可以将不重要的变量直接压缩为0，删除无效变量(未被腐蚀的支路)的影响。可以在参数缩减过程中，达到提取有效的支路电阻变化量
△rk
(已被腐蚀的支路)的作用。既实现了较为准确的接地网故障诊断，也实现了变量选择(降维)。解决了接地网腐蚀诊断的盲目性，无需大面积开挖土壤即可较为准确的进行故障诊断与定位，为接地网是否开挖进行更换以及开挖位置的选择提供参考，进而降低现场接地网的运营维护成本。通过图形化显示，使接地网腐蚀诊断结果形象直观，应用方便。
[0154]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，如图4所示，基于lasso理论的接地网腐蚀诊断方法，还包括：
[0155]
步骤6：训练腐蚀应对模型，基于所述腐蚀应对模型，根据变电站接地网的腐蚀程度，制定应对策略，基于所述应对策略，进行相应应对处理。
[0156]
上述技术方案的工作原理及有益效果为：
[0157]
当变电站接地网出现腐蚀情况时，需要及时进行应对处理；基于机器学习技术，训练能够学习人工分析变电站接地网腐蚀情况进行应对处理的模型，即腐蚀应对模型，基于训练的腐蚀应对模型，根据腐蚀程度，制定应对策略，进行相应应对处理。
[0158]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，所述步骤6中，训练腐蚀应对模型，包括：
[0159]
获取多个第一腐蚀应对事件，同时，获取所述第一腐蚀应对事件对应的应对方；
[0160]
获取所述应对方对应的应对类型，所述应对类型包括：内部应对和外部应对；
[0161]
当所述应对方对应的应对类型为内部应对时，获取所述应对方对应的经验程度值；
[0162]
若所述经验程度值小于等于预设的第一阈值，剔除对应所述第一腐蚀应对事件；
[0163]
当所述应对方对应的应对类型为外部应对时，获取所述应对方对应的信用程度值；
[0164]
若所述信用程度值小于等于预设的第二阈值，剔除对应所述第一腐蚀应对事件；
[0165]
当需要剔除的所述第一腐蚀应对事件均剔除后，将剔除剩余的所述第一腐蚀应对事件作为第二腐蚀应对事件；
[0166]
对所述第二腐蚀应对事件进行合理性分析，获取合理值；
[0167]
若所述合理值大于等于预设的第三阈值，将对应所述第二腐蚀应对事件作为第三腐蚀应对事件；
[0168]
基于预设的模型训练算法，根据所述第三腐蚀应对事件进行模型训练，获得腐蚀应对模型。
[0169]
上述技术方案的工作原理及有益效果为：
[0170]
训练腐蚀应对模型时，需要获取人工分析变电站接地网腐蚀情况进行应对处理的过程记录即第一腐蚀应对事件，便于基于机器学习技术，进行学习训练，但是，为了保证腐蚀应对模型的训练质量，以保证应对策略制定的适宜性，需要对第一腐蚀应对事件进行筛选；第一腐蚀应对事件对应的应对方(进行分析和应对处理指挥的人工方)的应对类型分为
内部应对(变电站公司内部的专业人员)和外部应对(其他变电站公司的专业人员)；当应对类型为内部应对时，应对方的经验程度可以精准追溯(有记录可查)，获取对应经验程度值，若经验程度值较小，剔除对应第一腐蚀应对事件；当应对类型为外部应对时，可以依据应对方历史上提供的腐蚀应对事件的总体质量进行验证即获取信用程度值，若信用程度值较小，剔除对应第一腐蚀应对事件；基于预设的模型训练算法(机器学习算法)，根据剔除剩余的第三腐蚀应对事件进行模型训练。
[0171]
本发明提供一种基于lasso理论的接地网腐蚀诊断方法，对所述第二腐蚀应对事件进行合理性分析，获取合理值，包括：
[0172]
获取多个专家节点，同时，获取所述专家节点对所述第二腐蚀应对事件进行合理性评价的评价值；
[0173]
获取所述专家节点对应的专家权重，赋予所述评价值对应所述专家权重，获得目标值；
[0174]
累加计算所述目标值，获得合理值。
[0175]
上述技术方案的工作原理及有益效果为：
[0176]
对第二腐蚀应对事件进行合理性分析时，设置多个专家节点，专家节点对应于一个接地网腐蚀应对专家，由专家节点对第二腐蚀应对事件进行合理性评价，获得评价值；专家节点的专家权重越大，说明专家节点进行合理性评价的可依据性越高，赋予评价值对应专家权重(两者相乘)，获得目标值；对目标值进行累加求和，即获得合理值，提升了合理值获取的精准性。
[0177]
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。