一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法
一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法
【技术领域】
[0001] 本发明涉及电网检测技术领域，尤其是涉及一种基于磁场逆问题求解的变电站接 地网腐蚀诊断方法。
【背景技术】
[0002] 接地网是变电站安全运行的重要保证，其工作的可靠性对电力系统的安全稳定运 行影响重大。由于国内的变电站接地网大多采用扁钢，圆钢等钢制材料，且其长期埋于地 下，不可避免地会受到土壤的电化学腐蚀，引起接地网均压导体腐蚀甚至断裂、接地引下线 截面变小、接地体之间接触不良等问题。这些故障均会使得接地网性能大幅下降，严重危害 变电站内工作人员及设备的安全。因此，找到一种能快速方便地对接地网腐蚀故障进行检 测的方法、进而采取相应的防护措施已成为电力行业现有接地网运维工作中亟待解决的问 题。
[0003] 目前，接地网腐蚀故障检测方法主要有三种：一是基于电路理论，通过建立故障诊 断方程并结合相应优化算法实现接地网的故障诊断，但该方法运算较为复杂，诊断方程病 态程度较高影响诊断精度，且当接地网仅局部存在轻微故障时难以通过该方法得到识别诊 断。二是基于电场理论，通过向接地网注入激励电流探测地表电位分布，从而对接地网进 行故障诊断，但该方法所需测点较多工作量大，且有可能出现由于现场检测不准确而导致 故障点的误判或漏判的情形。三是基于电磁场理论，通过向接地网注入激励电流探测地表 的磁感应强度，根据磁感应强度的分布特征对接地网进行故障诊断，但现有基于地表面磁 感应强度的腐蚀诊断方法测点布置需覆盖整个接地网，而实际应用中受制于变电站现场情 况，难以实现对所有测点的磁感应强度进行准确测量，且对大型接地网来说存在工作量较 大的问题。
【发明内容】
[0004] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于磁场逆问 题求解的变电站接地网腐蚀诊断方法，通过测试及计算分析变电站接地网的地表面磁感应 强度分布，能够实现对变电站接地网腐蚀故障的高效、准确判断。
[0005] 本发明的目的可以通过以下技术方案来实现：
[0006] 一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其特征在于，包括以下 步骤：
[0007] 步骤1、使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强 度；
[0008] 步骤2、应用交流信号发生器通过一可及节点向变电站接地网注入交流激励电流， 由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信 号，选取测点时，根据变电站现场情况，任意分散地选取n个便于探测的地表点为测点，其 中，n为接地网导体数；
[0009] 步骤3、根据采集到的感应电压信号计算n个测点处的磁感应强度；
[0010] 步骤4、根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段 导体的轴向电流；
[0011] 步骤5、分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将 所有导体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果；
[0012] 步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面 磁感应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应 强度曲线，并计算曲线偏差度；若某接地导体的磁感应强度曲线偏差度大于20%，则判断 该接地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。
[0013] 所述的步骤1)通过数值计算方法计算变电站接地网的地表面磁感应强度，具体 是：
[0014] la)将节点数目为m的接地网分为n段导体，计算这n段导体之间的互阻抗矩阵 R，其中，矩阵元素Ru表示i段导体和j段导体之间的互阻抗，其计算公式为：
[0016] 式中，土壤电导率；e ^为真空介电常数；e ^3土壤相对介电常数；e E= e ^ \为土壤介电常数；1 ,与h分别为第i段及第j段导体长度；1 ,，为第i段导体的镜 像长度；Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体 之间的距离；
[0017] lb)使用T型等效电路分别表示这n段导体，即1段导体对应1个T型等效电路， 所述的T型等效电路由第i段导体的自电感L、自电阻&、第i段和第j段导体之间的互感 M、第i段导体的对地电容C和对地电导G组成，其中i = 1，…，n，j = 1，…，n ;经T型电路 等效后，所述接地网共有m+n个节点及2n段导体；
[0018] lc)计算接地网经T型等效电路等效后各段导体的关联矩阵A，其中，关联矩阵A 的行对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目 2n，关联矩阵A中的任意元素的定义为：
[0020] Id)计算经T型等效电路等效后具有m+n个节点和2n条支路的接地网的阻抗矩阵 Z，其计算公式为：
[0022] 式中，My为接地网各支路经T型等效后的互感矩阵；k i为经T型等效后第i段 导体长度，其中i = 1，…，2n 』为第i段导体与第j段导体之间的距离；Z &为第i段 导体的内阻抗；h为土壤磁导率，并假设土壤和空气磁导率相同；
；y = 为导体磁导率；U ^3导体相对磁导率；P。为导体电阻率；r ^为导体半径；I Jyr。)、Mya)分别为零阶及一阶贝塞尔函数；
[0023] le)根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Y_= AZ 1T;根据阻抗 矩阵R计算导纳矩阵，为G = IT1;
[0024]If)计算接地网导体n段导体的中点电位:f和节点电位，通过求解接地网数 学模型的基本方程得到，所述的接地网数学模型的基本方程为：
[0026]式中，G为n X n矩阵；f为n个中点电位列向量；< 为m个节点列向量；
为m个节点注入电流列向量；
[0027] lg)计算接地网导体n段导体的轴向电流分布，通过方程.
求得； 其中，R"为第i段导体自阻抗； < 为第i段导体上的轴向电流； < 郝分别为第i段导 体两端点的节点电位；
[0028] lh)根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强 度，地面上任一点P处的磁感应强度计算公式为：
[0030] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y ^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0031] 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果。
[0032] 所述的步骤2)中对测点的选择方式大大减少了现场测试的工作量，对大型接地 网尤为明显。
[0033] 所述的步骤3)中的磁感应强度计算公式为：
[0035] 式中，Bini为磁感应强度某一方向分量的幅值；V。^为感应电压信号幅度；f。为激励 电流频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益。
[0036] 所述的步骤4)中计算接地网n段导体轴向电流，具体为：
[0037] 4a)建立n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式，其表达 式为：
[0038] DI1= B
[0039] 式中，D为与测点位置相关的nXn系数矩阵，n为接地网导体数；I1为n维导体轴 向电流向量；B为测点某一方向磁感应强度n维向量；
[0040] 4b)定义n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式的最小 化Tikhonov泛函的表达式为：
[0041] J^I1) = | IDI'-Bl |2+a| 1111 |2
[0042] 式中，a为正则化参数；
[0043] 4c)使用牛顿法迭代求解非线性方程| iDf-Bl | = S，其中，S为误差水平，其迭 代过程计算公式为：
[0044] f(ak) = | IDI^^-Bl |2-82= | | D (D TD+akE) _1DTB-B | |2-82
[0047] 式中，E为nXn维单位阵
[0048] 具体计算时，随机初始化设定正则化参数初值a。，计算f(ak)和f ' (ak)，此 处k = 0;进而依据公式ak+1 = a k_f (ak) /f ' (ak)迭代计算ak+1;依据公式I 1 (k+1)= (DTD+ak+iErVB 迭代计算 lYk+l);
[0049] 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体 轴向电流计算，若不满足则继续迭代。
[0050] 所述的步骤6)中的曲线偏差度计算公式X为：
[0052] 式中，N为磁感应强度曲线的点数，Ygi为第i段导体地表面磁感应强度的实际计 算结果为第i段导体地表面磁感应强度的理论计算结果。
[0053] 所述的步骤6)中绘制导体正常工况与实际 工况下的地表面磁感应强度曲线并计 算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而对接地网进行腐蚀 诊断。
[0054] 与现有技术相比，本发明可以在不开挖变电站的前提下，对接地网的腐蚀故障进 行较为准确的判断；同时，由于采用了Tikhonov正则化算法对磁场逆问题进行求解，现场 测试工作量大幅减少，也最大限度排除了由于个别测点的计算不准确或测试误差所导致的 错误诊断。
[0055] 因此本发明能在大幅减少现场测试工作量的前提下，有效地、准确地检测出变电 站接地网的腐蚀故障位置，从而可及时采取有效措施，提高变电站接地网的运行可靠性。
[0056] 本发明的进一步优点在于：步骤一中对接地网地表面磁感应强度的理论计算是基 于不等电位的接地网模型所进行的，对接地体单元采用集中参数取代分布参数，并将电路 理论与电场理论相结合，对地表面磁感应强度进行计算。
[0057] 本发明的进一步优点在于：步骤二中对测点的选择方式大大减少了现场测试的工 作量，对大型接地网尤为明显。
[0058] 本发明的进一步优点在于：步骤六中绘制导体正常工况与实际工况下的地表面磁 感应强度曲线并计算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而 对接地网进行腐蚀诊断。
【附图说明】
[0059] 图1是本发明实施例中的接地网模型示意图；
[0060] 图2是本发明接地网导体单元的T型等效图；
[0061] 图3 (a)为导体⑶的地表面磁感应强度理论计算结果；
[0062] 图3 (b)为导体EF的地表面磁感应强度理论计算结果；
[0063] 图4是本发明磁场逆问题求解流程图；
[0064] 图5是本发明接地网故障导体地表面磁感应强度的理论计算结果与根据实测得 到的实际计算结果。
【具体实施方式】
[0065] 下面结合附图和具体实施例对本发明进行详细说明。
[0066] 以图1所示的某变电站接地网模型为试验对象进行腐蚀诊断，图1中，AB段为模 拟腐蚀故障段，A点为电流注入点，X，y为坐标方向，坐标原点为0点，按照下列步骤判断该 接地网的工作状态：
[0067] (1)使用数值计算方法计算变电站接地网的地表面磁感应强度，所述的变电站接 地网地表面磁感应强度计算可用如下数值计算方法：
[0068] la.将节点数目为m的接地网分为n段导体，本实施例中m=25, n=40,计算这 n段导体之间的互阻抗矩阵R，其中，矩阵元素&表示i段导体和j段导体之间的互阻抗， 其计算公式为：
[0069]
[0070] 式中，0ES土壤电导率；e ^为真空介电常数；e 土壤相对介电常数；e E= e ^ \为土壤介电常数；1 ,与h分别为第i段及第j段导体长度；1 ,，为第i段导体的镜 像长度；Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体 之间的距离；
[0071] lb.使用T型等效电路分别表示这n段导体，即1段导体对应1个T型等效电路， 所述的T型等效电路由第i段导体（i = 1，…，n)的自电感L、自电阻4、第i段和第j段 (j = 1，…，n)导体之间的互感M、第i段导体的对地电容C和对地电导G组成，如图2所 示。经T型电路等效后，所述接地网共有m+n个节点及2n段导体，即65个节点和80段导 体；
[0072] lc.计算接地网经T型等效电路等效后各段导体的关联矩阵A，其中，关联矩阵A 的行对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目 2n，关联矩阵A中的任意元素ay的定义为：
[0074] Id.计算经T型等效电路等效后具有m+n个节点和2n条支路的接地网的阻抗矩阵 Z，其计算公式为：
[0076] 式中，My为接地网各支路经T型等效后的互感矩阵；k i为经T型等效后第i段导 体长度（i = 1，…，2n) 为第i段导体与第j段导体之间的距离；ZQi为第i段导体的内 阻抗；h为土壤磁导率，并假设土壤和空气磁导率相同；
；y = 为 导体磁导率；1^为导体相对磁导率；0。为导体电阻率；1' (|为导体半径；1(|(丫1'(|)、11(丫1' (|) 分别为零阶及一阶贝塞尔函数；
[0077] le.根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Y_= AZ11，此处，m+n =65 ;根据阻抗矩阵R计算导纳矩阵，为G = IT1;
[0078]If.计算接地网导体n段导体的中点电位f和节点电位:^，这可通过求解接地 网数学模型的基本方程得到，所述的接地网数学模型的基本方程为：
[0080] 式中，G为nXn矩阵； < 为n个中点电位列向量； < 为m个节点列向量；
为m个节点注入电流列向量；此处，m = 25, n = 40 ;
[0081] lg.计算接地网导体n段导体的轴向电流分布，这可通过方程.
求 得。其中，R"为第i段导体自阻抗；//为第i段导体上的轴向电流； < 和<1/分别为第i 段导体两端点的节点电位；此处，n = 40 ;
[0082] lh.根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强 度，地面上任一点P处的磁感应强度计算公式为：
[0084] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0085] 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果；部分导体 的地表面磁感应强度y方向分量计算结果如图3 (a)和3 (b)所示，其中图3 (a)为导体CD的 计算结果，图3 (b)为导体EF的计算结果，其中，C点坐标为（0, 3. 5)，D点坐标为（7, 3. 5)， E点坐标为（0, 1. 5)，F点坐标为（7, 1. 5);
[0086] (2)应用交流信号发生器通过一可及节点向所述变电站接地网注入交流激励电 流，由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压 信号。选取测点时，根据变电站现场情况，任意分散地选取n个便于探测的地表点为测点， 其中，n为所述接地网导体数；此处，n = 40 ;
[0087] (3)根据采集到的感应电压信号计算n个测点处的磁感应强度，所述的磁感应强 度计算公式为：
[0089] 式中，Bini为磁感应强度某一方向分量的幅值；V。_"为感应电压信号幅度；f。为激励 电流频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益；此 处，fc= 65Hz;N =350;S =6. 4Xl(T3m2;A =1 ;
[0090] (4)根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段导 体的轴向电流，所述的接地网n段导体轴向电流的计算过程为：
[0091] 4a.建立n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式，其表达 式为：
[0092] DI1= B
[0093] 式中，D为与测点位置相关的nXn系数矩阵（n为所述接地网导体数）f为n维 导体轴向电流向量；B为测点y方向磁感应强度n维向量；此处，n = 40 ;
[0094] 4b.定义n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式的最小 化Tikhonov泛函的表达式为：
[0095] J^I1) = | IDI'-Bl |2+a| 1111 |2
[0096] 式中，a为正则化参数；
[0097] 4c.使用牛顿法迭代求解非线性方程| iDf-Bl | = S，其中，S为误差水平，此处， S = 1(T2,其迭代过程计算公式为：
[0098] f(ak) = | IDI^^-Bl |2-82= | | D (D TD+akE) _1DTB-B | |2-82
[0101] 式中，E为nXn维单位阵，此处，n = 40 ;
[0102] 具体计算时，随机初始化设定正则化参数初值a。，计算f(ak)和f ' (ak)，此 处k = 0;进而依据公式ak+1 = a k_f (ak) /f ' (ak)迭代计算ak+1;依据公式I 1 (k+1)= (DTD+ak+iErVB 迭代计算 lYk+l);
[0103] 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体 轴向电流计算，若不满足则继续迭代；
[0104] (5)根据计算所得导体轴向电流，使用叠加原理计算地面上的磁感应强度，地面上 任一点P处的磁感应强度计算公式为：
[0106] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y ^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0107] 分别计算每段导体的轴向 电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果；
[0108] (6)将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感 应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应强度 曲线，并计算曲线偏差度。若某接地导体的磁感应强度曲线偏差度大于20%，则判断该接 地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。所述的曲线偏差度计算公式 为：
[0110] 式中，N为磁感应强度曲线的点数，此处，N = 21。
[0111] 图5显示了本实施例中根据上述方法得到的接地网模型AB段导体的理论地表面 磁感应强度曲线及根据实测结果计算所得的磁感应强度曲线，计算得到的曲线偏差度为 23. 3%，从而判断AB段导体为腐蚀故障段，所得结论与模拟故障设置一致，验证了本方法 的准确性。
[0112] 要注意的是，以上列举的仅为本发明的具体实施例，显然本发明不限于以上实施 例，随之有着许多的类似变化。本领域的技术人员如果从本发明公开的内容直接导出或联 想到的所有变形，均应属于本发明的保护范围。
【主权项】
1. 一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其特征在于，包括以下步 骤： 步骤1、使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强度； 步骤2、应用交流信号发生器通过一可及节点向变电站接地网注入交流激励电流，由探 测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信号，选 取测点时，根据变电站现场情况，任意分散地选取η个便于探测的地表点为测点，其中，η为 接地网导体数； 步骤3、根据采集到的感应电压信号计算η个测点处的磁感应强度； 步骤4、根据η个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网η段导体 的轴向电流； 步骤5、分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有 导体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果； 步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感 应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应强度 曲线，并计算曲线偏差度；若某接地导体的磁感应强度曲线偏差度大于20%，则判断该接 地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。2. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤1)通过数值计算方法计算变电站接地网的地表面磁感应强度，具体 是： la) 将节点数目为m的接地网分为η段导体，计算这η段导体之间的互阻抗矩阵R，其 中，矩阵元素 Rij表示i段导体和j段导体之间的互阻抗，其计算公式为：i = 1，…，n ;j = 1，…，η 式中，σ A土壤电导率；ε ^为真空介电常数；ε ^为土壤相对介电常数；ε Ε= ε ^ · ε ^ 为土壤介电常数；^与h分别为第i段及第j段导体长度；I i，为第i段导体的镜像长度； Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体之间的距 离； lb) 使用T型等效电路分别表示这η段导体，即1段导体对应1个T型等效电路，所述 的T型等效电路由第i段导体的自电感L、自电阻Ztl、第i段和第j段导体之间的互感Μ、第 i段导体的对地电容C和对地电导G组成，其中i = 1，…，n，j = 1，…，η ;经T型电路等效 后，所述接地网共有m+n个节点及2η段导体； lc) 计算接地网经T型等效电路等效后各段导体的关联矩阵Α，其中，关联矩阵A的行 对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目2η， 关联矩阵A中的任意元素的定义为：ld) 计算经T型等效电路等效后具有m+n个节点和2η条支路的接地网的阻抗矩阵Z， 其计算公式为：式中，Mu为接地网各支路经T型等效后的互感矩阵；k ,为经T型等效后第i段导体长 度，其中i = 1，…，2n 为第i段导体与第j段导体之间的距离；Ztli为第i段导体的内 阻抗；μ C1为土壤磁导率，并假设土壤和空气磁导率相同；μ = μ(ιμ^3 导体磁导率；μ 导体相对磁导率；P。为导体电阻率；r ^为导体半径；I Jyrtl)、I1(Yrtl) 分别为零阶及一阶贝塞尔函数； le) 根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Yn+m= AZ ^1At;根据阻抗矩阵R 计算导纳矩阵，为G = IT1; lf) 计算接地网导体η段导体的中点电位和节点电位，通过求解接地网数学模 型的基本方程得到，所述的接地网数学模型的基本方程为：式中，G为η X η矩阵；f为η个中点电位列向量； < 为m个节点列向量； /OT = [CC，…,/二』为m个节点注入电流列向量； lg) 计算接地网导体η段导体的轴向电流分布，通过方程尺/; 求得；其中， Rii为第i段导体自阻抗；//为第i段导体上的轴向电流； < 和分别为第i段导体两端 点的节点电位； lh) 根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强度， 地面上任一点P处的磁感应强度计算公式为：式中，Bip为第i段导体轴向电流在点P形成的磁感应强度；μ ^为真空磁导率；r为第 i段导体中点与点P之间的距离； 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导体的 计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果。3. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤2)中对测点的选择方式大大减少了现场测试的工作量，对大型接地 网尤为明显。4. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤3)中的磁感应强度计算公式为：式中，Bini为磁感应强度某一方向分量的幅值；VmiS感应电压信号幅度；f。为激励电流 频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益。5. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤4)中计算接地网η段导体轴向电流，具体为： 4a)建立η个测点处的磁感应强度与η段导体轴向电流之间的数学关系式，其表达式 为： DI1= B 式中，D为与测点位置相关的ηΧη系数矩阵，η为接地网导体数；I1为η维导体轴向电 流向量；B为测点某一方向磁感应强度η维向量； 4b)定义η个测点处的磁感应强度与η段导体轴向电流之间的数学关系式的最小化 Tikhonov泛函的表达式为： Ja(Ix) = I IDI1-Bl |2+a| Il1I I2 式中，a为正则化参数； 4c)使用牛顿法迭代求解非线性方程I IDI1-Bl I = δ，其中，δ为误差水平，其迭代过 程计算公式为：式中，E为ηΧη维单位阵 具体计算时，随机初始化设定正则化参数初值a。，计算f (ak)和f' (ak)，此处k = 0 ; 进而依据公式 ak+1=ak_f(ak)/f' (ak)迭代计算 ak+1^据公式 iHk+l) = (DTD+aj^ErWB 迭代计算iHk+l); 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体轴向 电流计算，若不满足则继续迭代。6. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤6)中的曲线偏差度计算公式λ为：式中，N为磁感应强度曲线的点数，Ygi为第i段导体地表面磁感应强度的实际计算结 果；Yi为第i段导体地表面磁感应强度的理论计算结果。7. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤6)中绘制导体正常工况与实际工况下的地表面磁感应强度曲线并 计算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而对接地网进行腐 蚀诊断。
【专利摘要】本发明涉及一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，包括以下步骤：使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强度；向变电站接地网注入交流激励电流，由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信号；根据采集到的感应电压信号计算n个测点处的磁感应强度；根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段导体的轴向电流；得到接地网的地表面磁感应强度的实际计算结果；步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感应强度进行比较。与现有技术相比，本发明具有提高变电站接地网的运行可靠性等优点。
【IPC分类】G01R31/00, G01R31/08
【公开号】CN104897996
【申请号】CN201510316304
【发明人】司文荣, 傅晨钊, 苏磊, 陈璐
【申请人】国网上海市电力公司, 华东电力试验研究院有限公司
【公开日】2015年9月9日
【申请日】2015年6月10日
【技术领域】
[0001] 本发明涉及电网检测技术领域，尤其是涉及一种基于磁场逆问题求解的变电站接 地网腐蚀诊断方法。
【背景技术】
[0002] 接地网是变电站安全运行的重要保证，其工作的可靠性对电力系统的安全稳定运 行影响重大。由于国内的变电站接地网大多采用扁钢，圆钢等钢制材料，且其长期埋于地 下，不可避免地会受到土壤的电化学腐蚀，引起接地网均压导体腐蚀甚至断裂、接地引下线 截面变小、接地体之间接触不良等问题。这些故障均会使得接地网性能大幅下降，严重危害 变电站内工作人员及设备的安全。因此，找到一种能快速方便地对接地网腐蚀故障进行检 测的方法、进而采取相应的防护措施已成为电力行业现有接地网运维工作中亟待解决的问 题。
[0003] 目前，接地网腐蚀故障检测方法主要有三种：一是基于电路理论，通过建立故障诊 断方程并结合相应优化算法实现接地网的故障诊断，但该方法运算较为复杂，诊断方程病 态程度较高影响诊断精度，且当接地网仅局部存在轻微故障时难以通过该方法得到识别诊 断。二是基于电场理论，通过向接地网注入激励电流探测地表电位分布，从而对接地网进 行故障诊断，但该方法所需测点较多工作量大，且有可能出现由于现场检测不准确而导致 故障点的误判或漏判的情形。三是基于电磁场理论，通过向接地网注入激励电流探测地表 的磁感应强度，根据磁感应强度的分布特征对接地网进行故障诊断，但现有基于地表面磁 感应强度的腐蚀诊断方法测点布置需覆盖整个接地网，而实际应用中受制于变电站现场情 况，难以实现对所有测点的磁感应强度进行准确测量，且对大型接地网来说存在工作量较 大的问题。
【发明内容】
[0004] 本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于磁场逆问 题求解的变电站接地网腐蚀诊断方法，通过测试及计算分析变电站接地网的地表面磁感应 强度分布，能够实现对变电站接地网腐蚀故障的高效、准确判断。
[0005] 本发明的目的可以通过以下技术方案来实现：
[0006] 一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其特征在于，包括以下 步骤：
[0007] 步骤1、使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强 度；
[0008] 步骤2、应用交流信号发生器通过一可及节点向变电站接地网注入交流激励电流， 由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信 号，选取测点时，根据变电站现场情况，任意分散地选取n个便于探测的地表点为测点，其 中，n为接地网导体数；
[0009] 步骤3、根据采集到的感应电压信号计算n个测点处的磁感应强度；
[0010] 步骤4、根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段 导体的轴向电流；
[0011] 步骤5、分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将 所有导体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果；
[0012] 步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面 磁感应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应 强度曲线，并计算曲线偏差度；若某接地导体的磁感应强度曲线偏差度大于20%，则判断 该接地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。
[0013] 所述的步骤1)通过数值计算方法计算变电站接地网的地表面磁感应强度，具体 是：
[0014] la)将节点数目为m的接地网分为n段导体，计算这n段导体之间的互阻抗矩阵 R，其中，矩阵元素Ru表示i段导体和j段导体之间的互阻抗，其计算公式为：
[0016] 式中，土壤电导率；e ^为真空介电常数；e ^3土壤相对介电常数；e E= e ^ \为土壤介电常数；1 ,与h分别为第i段及第j段导体长度；1 ,，为第i段导体的镜 像长度；Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体 之间的距离；
[0017] lb)使用T型等效电路分别表示这n段导体，即1段导体对应1个T型等效电路， 所述的T型等效电路由第i段导体的自电感L、自电阻&、第i段和第j段导体之间的互感 M、第i段导体的对地电容C和对地电导G组成，其中i = 1，…，n，j = 1，…，n ;经T型电路 等效后，所述接地网共有m+n个节点及2n段导体；
[0018] lc)计算接地网经T型等效电路等效后各段导体的关联矩阵A，其中，关联矩阵A 的行对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目 2n，关联矩阵A中的任意元素的定义为：
[0020] Id)计算经T型等效电路等效后具有m+n个节点和2n条支路的接地网的阻抗矩阵 Z，其计算公式为：
[0022] 式中，My为接地网各支路经T型等效后的互感矩阵；k i为经T型等效后第i段 导体长度，其中i = 1，…，2n 』为第i段导体与第j段导体之间的距离；Z &为第i段 导体的内阻抗；h为土壤磁导率，并假设土壤和空气磁导率相同；
；y = 为导体磁导率；U ^3导体相对磁导率；P。为导体电阻率；r ^为导体半径；I Jyr。)、Mya)分别为零阶及一阶贝塞尔函数；
[0023] le)根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Y_= AZ 1T;根据阻抗 矩阵R计算导纳矩阵，为G = IT1;
[0024]If)计算接地网导体n段导体的中点电位:f和节点电位，通过求解接地网数 学模型的基本方程得到，所述的接地网数学模型的基本方程为：
[0026]式中，G为n X n矩阵；f为n个中点电位列向量；< 为m个节点列向量；
为m个节点注入电流列向量；
[0027] lg)计算接地网导体n段导体的轴向电流分布，通过方程.
求得； 其中，R"为第i段导体自阻抗； < 为第i段导体上的轴向电流； < 郝分别为第i段导 体两端点的节点电位；
[0028] lh)根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强 度，地面上任一点P处的磁感应强度计算公式为：
[0030] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y ^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0031] 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果。
[0032] 所述的步骤2)中对测点的选择方式大大减少了现场测试的工作量，对大型接地 网尤为明显。
[0033] 所述的步骤3)中的磁感应强度计算公式为：
[0035] 式中，Bini为磁感应强度某一方向分量的幅值；V。^为感应电压信号幅度；f。为激励 电流频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益。
[0036] 所述的步骤4)中计算接地网n段导体轴向电流，具体为：
[0037] 4a)建立n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式，其表达 式为：
[0038] DI1= B
[0039] 式中，D为与测点位置相关的nXn系数矩阵，n为接地网导体数；I1为n维导体轴 向电流向量；B为测点某一方向磁感应强度n维向量；
[0040] 4b)定义n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式的最小 化Tikhonov泛函的表达式为：
[0041] J^I1) = | IDI'-Bl |2+a| 1111 |2
[0042] 式中，a为正则化参数；
[0043] 4c)使用牛顿法迭代求解非线性方程| iDf-Bl | = S，其中，S为误差水平，其迭 代过程计算公式为：
[0044] f(ak) = | IDI^^-Bl |2-82= | | D (D TD+akE) _1DTB-B | |2-82
[0047] 式中，E为nXn维单位阵
[0048] 具体计算时，随机初始化设定正则化参数初值a。，计算f(ak)和f ' (ak)，此 处k = 0;进而依据公式ak+1 = a k_f (ak) /f ' (ak)迭代计算ak+1;依据公式I 1 (k+1)= (DTD+ak+iErVB 迭代计算 lYk+l);
[0049] 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体 轴向电流计算，若不满足则继续迭代。
[0050] 所述的步骤6)中的曲线偏差度计算公式X为：
[0052] 式中，N为磁感应强度曲线的点数，Ygi为第i段导体地表面磁感应强度的实际计 算结果为第i段导体地表面磁感应强度的理论计算结果。
[0053] 所述的步骤6)中绘制导体正常工况与实际 工况下的地表面磁感应强度曲线并计 算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而对接地网进行腐蚀 诊断。
[0054] 与现有技术相比，本发明可以在不开挖变电站的前提下，对接地网的腐蚀故障进 行较为准确的判断；同时，由于采用了Tikhonov正则化算法对磁场逆问题进行求解，现场 测试工作量大幅减少，也最大限度排除了由于个别测点的计算不准确或测试误差所导致的 错误诊断。
[0055] 因此本发明能在大幅减少现场测试工作量的前提下，有效地、准确地检测出变电 站接地网的腐蚀故障位置，从而可及时采取有效措施，提高变电站接地网的运行可靠性。
[0056] 本发明的进一步优点在于：步骤一中对接地网地表面磁感应强度的理论计算是基 于不等电位的接地网模型所进行的，对接地体单元采用集中参数取代分布参数，并将电路 理论与电场理论相结合，对地表面磁感应强度进行计算。
[0057] 本发明的进一步优点在于：步骤二中对测点的选择方式大大减少了现场测试的工 作量，对大型接地网尤为明显。
[0058] 本发明的进一步优点在于：步骤六中绘制导体正常工况与实际工况下的地表面磁 感应强度曲线并计算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而 对接地网进行腐蚀诊断。
【附图说明】
[0059] 图1是本发明实施例中的接地网模型示意图；
[0060] 图2是本发明接地网导体单元的T型等效图；
[0061] 图3 (a)为导体⑶的地表面磁感应强度理论计算结果；
[0062] 图3 (b)为导体EF的地表面磁感应强度理论计算结果；
[0063] 图4是本发明磁场逆问题求解流程图；
[0064] 图5是本发明接地网故障导体地表面磁感应强度的理论计算结果与根据实测得 到的实际计算结果。
【具体实施方式】
[0065] 下面结合附图和具体实施例对本发明进行详细说明。
[0066] 以图1所示的某变电站接地网模型为试验对象进行腐蚀诊断，图1中，AB段为模 拟腐蚀故障段，A点为电流注入点，X，y为坐标方向，坐标原点为0点，按照下列步骤判断该 接地网的工作状态：
[0067] (1)使用数值计算方法计算变电站接地网的地表面磁感应强度，所述的变电站接 地网地表面磁感应强度计算可用如下数值计算方法：
[0068] la.将节点数目为m的接地网分为n段导体，本实施例中m=25, n=40,计算这 n段导体之间的互阻抗矩阵R，其中，矩阵元素&表示i段导体和j段导体之间的互阻抗， 其计算公式为：
[0069]
[0070] 式中，0ES土壤电导率；e ^为真空介电常数；e 土壤相对介电常数；e E= e ^ \为土壤介电常数；1 ,与h分别为第i段及第j段导体长度；1 ,，为第i段导体的镜 像长度；Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体 之间的距离；
[0071] lb.使用T型等效电路分别表示这n段导体，即1段导体对应1个T型等效电路， 所述的T型等效电路由第i段导体（i = 1，…，n)的自电感L、自电阻4、第i段和第j段 (j = 1，…，n)导体之间的互感M、第i段导体的对地电容C和对地电导G组成，如图2所 示。经T型电路等效后，所述接地网共有m+n个节点及2n段导体，即65个节点和80段导 体；
[0072] lc.计算接地网经T型等效电路等效后各段导体的关联矩阵A，其中，关联矩阵A 的行对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目 2n，关联矩阵A中的任意元素ay的定义为：
[0074] Id.计算经T型等效电路等效后具有m+n个节点和2n条支路的接地网的阻抗矩阵 Z，其计算公式为：
[0076] 式中，My为接地网各支路经T型等效后的互感矩阵；k i为经T型等效后第i段导 体长度（i = 1，…，2n) 为第i段导体与第j段导体之间的距离；ZQi为第i段导体的内 阻抗；h为土壤磁导率，并假设土壤和空气磁导率相同；
；y = 为 导体磁导率；1^为导体相对磁导率；0。为导体电阻率；1' (|为导体半径；1(|(丫1'(|)、11(丫1' (|) 分别为零阶及一阶贝塞尔函数；
[0077] le.根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Y_= AZ11，此处，m+n =65 ;根据阻抗矩阵R计算导纳矩阵，为G = IT1;
[0078]If.计算接地网导体n段导体的中点电位f和节点电位:^，这可通过求解接地 网数学模型的基本方程得到，所述的接地网数学模型的基本方程为：
[0080] 式中，G为nXn矩阵； < 为n个中点电位列向量； < 为m个节点列向量；
为m个节点注入电流列向量；此处，m = 25, n = 40 ;
[0081] lg.计算接地网导体n段导体的轴向电流分布，这可通过方程.
求 得。其中，R"为第i段导体自阻抗；//为第i段导体上的轴向电流； < 和<1/分别为第i 段导体两端点的节点电位；此处，n = 40 ;
[0082] lh.根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强 度，地面上任一点P处的磁感应强度计算公式为：
[0084] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0085] 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果；部分导体 的地表面磁感应强度y方向分量计算结果如图3 (a)和3 (b)所示，其中图3 (a)为导体CD的 计算结果，图3 (b)为导体EF的计算结果，其中，C点坐标为（0, 3. 5)，D点坐标为（7, 3. 5)， E点坐标为（0, 1. 5)，F点坐标为（7, 1. 5);
[0086] (2)应用交流信号发生器通过一可及节点向所述变电站接地网注入交流激励电 流，由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压 信号。选取测点时，根据变电站现场情况，任意分散地选取n个便于探测的地表点为测点， 其中，n为所述接地网导体数；此处，n = 40 ;
[0087] (3)根据采集到的感应电压信号计算n个测点处的磁感应强度，所述的磁感应强 度计算公式为：
[0089] 式中，Bini为磁感应强度某一方向分量的幅值；V。_"为感应电压信号幅度；f。为激励 电流频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益；此 处，fc= 65Hz;N =350;S =6. 4Xl(T3m2;A =1 ;
[0090] (4)根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段导 体的轴向电流，所述的接地网n段导体轴向电流的计算过程为：
[0091] 4a.建立n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式，其表达 式为：
[0092] DI1= B
[0093] 式中，D为与测点位置相关的nXn系数矩阵（n为所述接地网导体数）f为n维 导体轴向电流向量；B为测点y方向磁感应强度n维向量；此处，n = 40 ;
[0094] 4b.定义n个测点处的磁感应强度与n段导体轴向电流之间的数学关系式的最小 化Tikhonov泛函的表达式为：
[0095] J^I1) = | IDI'-Bl |2+a| 1111 |2
[0096] 式中，a为正则化参数；
[0097] 4c.使用牛顿法迭代求解非线性方程| iDf-Bl | = S，其中，S为误差水平，此处， S = 1(T2,其迭代过程计算公式为：
[0098] f(ak) = | IDI^^-Bl |2-82= | | D (D TD+akE) _1DTB-B | |2-82
[0101] 式中，E为nXn维单位阵，此处，n = 40 ;
[0102] 具体计算时，随机初始化设定正则化参数初值a。，计算f(ak)和f ' (ak)，此 处k = 0;进而依据公式ak+1 = a k_f (ak) /f ' (ak)迭代计算ak+1;依据公式I 1 (k+1)= (DTD+ak+iErVB 迭代计算 lYk+l);
[0103] 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体 轴向电流计算，若不满足则继续迭代；
[0104] (5)根据计算所得导体轴向电流，使用叠加原理计算地面上的磁感应强度，地面上 任一点P处的磁感应强度计算公式为：
[0106] 式中，BiP为第i段导体轴向电流在点P形成的磁感应强度；y ^为真空磁导率；r 为第i段导体中点与点P之间的距离；
[0107] 分别计算每段导体的轴向 电流在接地网各地表点形成的磁感应强度，再将所有导 体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果；
[0108] (6)将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感 应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应强度 曲线，并计算曲线偏差度。若某接地导体的磁感应强度曲线偏差度大于20%，则判断该接 地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。所述的曲线偏差度计算公式 为：
[0110] 式中，N为磁感应强度曲线的点数，此处，N = 21。
[0111] 图5显示了本实施例中根据上述方法得到的接地网模型AB段导体的理论地表面 磁感应强度曲线及根据实测结果计算所得的磁感应强度曲线，计算得到的曲线偏差度为 23. 3%，从而判断AB段导体为腐蚀故障段，所得结论与模拟故障设置一致，验证了本方法 的准确性。
[0112] 要注意的是，以上列举的仅为本发明的具体实施例，显然本发明不限于以上实施 例，随之有着许多的类似变化。本领域的技术人员如果从本发明公开的内容直接导出或联 想到的所有变形，均应属于本发明的保护范围。
【主权项】
1. 一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其特征在于，包括以下步 骤： 步骤1、使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强度； 步骤2、应用交流信号发生器通过一可及节点向变电站接地网注入交流激励电流，由探 测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信号，选 取测点时，根据变电站现场情况，任意分散地选取η个便于探测的地表点为测点，其中，η为 接地网导体数； 步骤3、根据采集到的感应电压信号计算η个测点处的磁感应强度； 步骤4、根据η个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网η段导体 的轴向电流； 步骤5、分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有 导体的计算结果进行叠加，从而得到接地网的地表面磁感应强度的实际计算结果； 步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感 应强度进行比较，对同一接地导体，分别绘制正常工况与实际工况下的地表面磁感应强度 曲线，并计算曲线偏差度；若某接地导体的磁感应强度曲线偏差度大于20%，则判断该接 地导体存在腐蚀故障，从而完成对变电站接地网的腐蚀诊断。2. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤1)通过数值计算方法计算变电站接地网的地表面磁感应强度，具体 是： la) 将节点数目为m的接地网分为η段导体，计算这η段导体之间的互阻抗矩阵R，其 中，矩阵元素 Rij表示i段导体和j段导体之间的互阻抗，其计算公式为：i = 1，…，n ;j = 1，…，η 式中，σ A土壤电导率；ε ^为真空介电常数；ε ^为土壤相对介电常数；ε Ε= ε ^ · ε ^ 为土壤介电常数；^与h分别为第i段及第j段导体长度；I i，为第i段导体的镜像长度； Du为将第i段与第j段导体之间的距离；D 为将第i段导体镜像与第j段导体之间的距 离； lb) 使用T型等效电路分别表示这η段导体，即1段导体对应1个T型等效电路，所述 的T型等效电路由第i段导体的自电感L、自电阻Ztl、第i段和第j段导体之间的互感Μ、第 i段导体的对地电容C和对地电导G组成，其中i = 1，…，n，j = 1，…，η ;经T型电路等效 后，所述接地网共有m+n个节点及2η段导体； lc) 计算接地网经T型等效电路等效后各段导体的关联矩阵Α，其中，关联矩阵A的行 对应于接地网经T型等效电路等效后的节点数目m+n，关联矩阵A的列对应于支路数目2η， 关联矩阵A中的任意元素的定义为：ld) 计算经T型等效电路等效后具有m+n个节点和2η条支路的接地网的阻抗矩阵Z， 其计算公式为：式中，Mu为接地网各支路经T型等效后的互感矩阵；k ,为经T型等效后第i段导体长 度，其中i = 1，…，2n 为第i段导体与第j段导体之间的距离；Ztli为第i段导体的内 阻抗；μ C1为土壤磁导率，并假设土壤和空气磁导率相同；μ = μ(ιμ^3 导体磁导率；μ 导体相对磁导率；P。为导体电阻率；r ^为导体半径；I Jyrtl)、I1(Yrtl) 分别为零阶及一阶贝塞尔函数； le) 根据关联矩阵A和阻抗矩阵Z计算节点导纳矩阵，为Yn+m= AZ ^1At;根据阻抗矩阵R 计算导纳矩阵，为G = IT1; lf) 计算接地网导体η段导体的中点电位和节点电位，通过求解接地网数学模 型的基本方程得到，所述的接地网数学模型的基本方程为：式中，G为η X η矩阵；f为η个中点电位列向量； < 为m个节点列向量； /OT = [CC，…,/二』为m个节点注入电流列向量； lg) 计算接地网导体η段导体的轴向电流分布，通过方程尺/; 求得；其中， Rii为第i段导体自阻抗；//为第i段导体上的轴向电流； < 和分别为第i段导体两端 点的节点电位； lh) 根据求得的接地网导体的轴向电流分布使用叠加原理计算地面上的磁感应强度， 地面上任一点P处的磁感应强度计算公式为：式中，Bip为第i段导体轴向电流在点P形成的磁感应强度；μ ^为真空磁导率；r为第 i段导体中点与点P之间的距离； 分别计算每段导体的轴向电流在接地网各地表点形成的磁感应强度，再将所有导体的 计算结果进行叠加，从而得到接地网的地表面磁感应强度的理论计算结果。3. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤2)中对测点的选择方式大大减少了现场测试的工作量，对大型接地 网尤为明显。4. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤3)中的磁感应强度计算公式为：式中，Bini为磁感应强度某一方向分量的幅值；VmiS感应电压信号幅度；f。为激励电流 频率；N为探测线圈匝数；S为探测线圈截面积；A为数据采集器对信号的放大增益。5. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤4)中计算接地网η段导体轴向电流，具体为： 4a)建立η个测点处的磁感应强度与η段导体轴向电流之间的数学关系式，其表达式 为： DI1= B 式中，D为与测点位置相关的ηΧη系数矩阵，η为接地网导体数；I1为η维导体轴向电 流向量；B为测点某一方向磁感应强度η维向量； 4b)定义η个测点处的磁感应强度与η段导体轴向电流之间的数学关系式的最小化 Tikhonov泛函的表达式为： Ja(Ix) = I IDI1-Bl |2+a| Il1I I2 式中，a为正则化参数； 4c)使用牛顿法迭代求解非线性方程I IDI1-Bl I = δ，其中，δ为误差水平，其迭代过 程计算公式为：式中，E为ηΧη维单位阵 具体计算时，随机初始化设定正则化参数初值a。，计算f (ak)和f' (ak)，此处k = 0 ; 进而依据公式 ak+1=ak_f(ak)/f' (ak)迭代计算 ak+1^据公式 iHk+l) = (DTD+aj^ErWB 迭代计算iHk+l); 4d.判断当前轴向电流向量是否满足误差水平限制，若满足则结束迭代完成导体轴向 电流计算，若不满足则继续迭代。6. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤6)中的曲线偏差度计算公式λ为：式中，N为磁感应强度曲线的点数，Ygi为第i段导体地表面磁感应强度的实际计算结 果；Yi为第i段导体地表面磁感应强度的理论计算结果。7. 根据权利要求1所述的一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，其 特征在于，所述的步骤6)中绘制导体正常工况与实际工况下的地表面磁感应强度曲线并 计算其偏差度，可直观地反映出导体地表面磁感应强度的变化情况，进而对接地网进行腐 蚀诊断。
【专利摘要】本发明涉及一种基于磁场逆问题求解的变电站接地网腐蚀诊断方法，包括以下步骤：使用数值计算方法计算正常工况下变电站接地网的理论地表面磁感应强度；向变电站接地网注入交流激励电流，由探测线圈沿接地网上方感应地表面磁场，并由数据采集器采集探测线圈的感应电压信号；根据采集到的感应电压信号计算n个测点处的磁感应强度；根据n个测点处的磁感应强度，采用Tikhonov正则化算法计算接地网n段导体的轴向电流；得到接地网的地表面磁感应强度的实际计算结果；步骤6、将实际的接地网地表面磁感应强度与正常工况下理论计算所得的地表面磁感应强度进行比较。与现有技术相比，本发明具有提高变电站接地网的运行可靠性等优点。
【IPC分类】G01R31/00, G01R31/08
【公开号】CN104897996
【申请号】CN201510316304
【发明人】司文荣, 傅晨钊, 苏磊, 陈璐
【申请人】国网上海市电力公司, 华东电力试验研究院有限公司
【公开日】2015年9月9日
【申请日】2015年6月10日
最新回复(0)