一种基于数字孪生的避雷器运行状态检测方法与流程

本发明属于避雷器运行状态检测，尤其涉及一种基于数字孪生的避雷器运行状态检测方法。

背景技术：

1、避雷器是电力系统中抑制过电压、释放冲击大电流和保护电力设备不受雷击损坏的重要保护装置，其运行状态极大影响到电网系统的稳定性。然而由于长期受外界环境、供电电压、电网谐波及设备内部受潮、老化现象等因素的影响，避雷器内部阀片的绝缘性能和灵敏度会出现不同程度的偏差，使得传统在线检测方法对避雷器误差较大，如何建立考虑老化、受潮等因素影响的避雷器模型，成为避雷器在线检测的难题。同时现电网站内设备往往排列较为密集，带电导体之间存在的耦合电容容易产生相间及空间干扰电流，会对避雷器泄漏电流等检测带来较大干扰，需要建立考虑相间和空间干扰的相关模型和算法。目前最常用的避雷器模型为仅考虑了正常状态的电阻和电容特性的小电流区模型，该模型对于泄漏电流的检测虽然简单有效，但随着老化、受潮等影响，误差会逐渐变大，使得在线检测有效性不足。

技术实现思路

1、本发明要解决的技术问题是：提供一种基于数字孪生的避雷器运行状态检测方法，以解决常用的避雷器模型为仅考虑了正常状态的电阻和电容特性的小电流区模型，该模型对于泄漏电流的检测虽然简单有效，但随着老化、受潮等影响，误差会逐渐变大，使得在线检测有效性不足等技术问题。

2、本发明技术方案是：

3、一种基于数字孪生的避雷器运行状态检测方法，所述方法包括：

4、步骤1：获取待建模避雷器的设备参数、运行数据以及外部信息，获取避雷器在电网拓扑中的安装信息；

5、步骤2：从电学角度对避雷器实体进行特征映射，根据改进的阻性电流算法建立考虑避雷器老化因素的避雷器数字孪生模型；

6、步骤3：考虑相间杂散电容干扰，对避雷器数字孪生模型进行修正；

7、步骤4：基于数字孪生模型，考虑谐波影响的避雷器在线检测；

8、步骤5：基于数字孪生模型信息，利用谐波分析法提出避雷器老化状态评估，根据计算结果将避雷器分为未老化、老化中期和完全老化三个等级。

9、获取待建模避雷器设备参数包括电气参数和结构参数；电气参数包括伏安特性曲线、额定电压、持续运行电压、标称放电电流下残压和最大放电电流，结构参数包括避雷器型号、阀片材质、质量和比热容参数。

10、电网拓扑中的安装信息的获取方法为：获取电网电压等级，明确避雷器的接入方式与保护范围，得到电网拓扑图和避雷器在电网安装分布图。

11、步骤2具体包括：

12、步骤2.1、基于获取的避雷器设备参数，建立避雷器基本数字孪生模型；

13、根据避雷器伏安特性曲线将避雷器等效为一个电容c支路与一个非线性电阻r支路的并联电路；避雷器的泄漏电流ix分为阻性电流ir和容性电流ic两部分，计算公式如下：

14、其中u为施加于避雷器两端电压，k为避雷器自身老化程度决定的常数：

15、

16、步骤2.2、基于容性电流补偿法消去容性电流，从而得到阻性电流；考虑谐波对容性电流的影响，容性电流分量根据u计算得到，如下式：

17、

18、ir＝ix-ic

19、其中ω为交流电频率，n为谐波次数；

20、步骤2.3、根据改进的阻性电流算法，对非线性等效电阻进行修正，建立考虑避雷器老化对电阻影响的避雷器数字孪生模型；

21、将非线性电阻r产生基波阻性电流的电阻等效为线性电阻r1，其余部分等效为非线性电阻r2，从而使非线性电阻r等效为r1与r2并联组成；根据基本孪生电路模型，对阻性电流ir进行fft分解，如下式所示：

22、

23、

24、

25、其中θx为x次谐波电流的相位；

26、计算x次谐波在基波处的阻性电流幅值imx1，如下式所示：

27、

28、计算得到r1线性等效电阻，然后根据r1得到r2，如下式所示，其中θ为基波电流的相位：

29、

30、r2＝tanu/k-r1。

31、步骤3对避雷器数字孪生模型进行修正具体包括：

32、步骤3.1、考虑相间杂散电容干扰，对避雷器数字孪生模型进行修正，设相间杂散电容为cab、cba、cbc、cbc，当避雷器以一字形方式接入三相电路中时，四个杂散电容参数相等，则通过谐波分析法得到各项阻性电流的真实值修正公式，如下式所示，其中为当前时刻电流与电压夹角，为各相电流与电压的夹角：

33、

34、

35、

36、步骤3.2、采集避雷器的三相泄漏电流ia、ib、ic以及三相电压信号ua、ub、uc，并计算对应的幅值和相位信息；

37、步骤3.3、计算三相电流的零序和正序电流分量，计算b相补偿后的全电流幅值和夹角信息；

38、步骤3.4、计算b相初始运行时刻补偿后的幅值和夹角信息，并计算对应的阻性电流，进而按步骤2计算修正避雷器电容与电阻参数；

39、步骤3.5、各相均以b相计算的初始值为基准，加上各相幅值及夹角的变化量作为计算结果，并计算对应的阻性电流，进而按步骤2计算修正避雷器电容与电阻参数。

40、步骤4基于数字孪生模型，考虑谐波影响的避雷器在线检测的具体方法包括：

41、步骤4.1、泄漏电流包括阻性分量和容性分量，如下式所示，

42、

43、其中αk为各次谐波电流的相角；

44、步骤4.2、两边同时乘sin(nωt+αn)，并对等式两边做在一周期内的定积分计算，化简得到下式：

45、irk＝ikmcos(βk-αk)

46、步骤4.3、两边同时乘cos(nωt+αn)，并对等式两边做在一周期内的定积分计算，化简得到下式：

47、ick＝ikmsin(βk-αk)

48、实时计算各相全电流、阻性电流和容性电流，并根据正常运行和各类故障运行时计算的各电流大小，设定各相全电流、阻性电流和容性电流的告警阈值，从而实现避雷器本体的故障预警。

49、步骤5具体包括：

50、步骤5.1、对步骤4中泄漏电流的阻性电流分量计算进行简化，由下式快速计算得到，其中m为中间常数变量，可来利用平均功率定义求解，其中u为避雷器测得的端电压，i为测得的泄漏电流：

51、ir＝u/m

52、

53、步骤5.2、泄露电流的基波阻性电流ir1可由下式计算得到，

54、ir1＝u1/m

55、其中基波电压u1可通过对电压fft计算得到；

56、(3)根据步骤2数字孪生模型仿真计算得到k与泄露电流的基波阻性电流ir1关系，并分别按照k＝40,25,12对应的基波阻性电流ir1值作为避雷器未老化、老化中期和完全老化三个等级的判断值。

57、当测得避雷器ir1小于k＝40对应ir1值时，为未老化状态；当测得避雷器ir1大于k＝40对应ir1值且小于k＝25对应ir1值时，为老化中期状态；当测得避雷器ir1大于k＝12对应ir1值时，为完全老化状态。

58、本发明的有益效果是：

59、本发明建立的数字孪生模型考虑了避雷器老化、相间杂散电容等因素，相比传统仿真模型对避雷器运行状态计算更精确，减小了避雷器内部阀片老化导致绝缘性能降低带来的误差，有利于更准确在线检测避雷器运行状态。同时提出基于数字孪生模型的在线检测方法计算速度快，更适用于实际工程应用中，而提出的避雷器老化状态评估方法能够清楚评估避雷器老化状况，有利于安排检修等工作。

60、解决了常用的避雷器模型为仅考虑了正常状态的电阻和电容特性的小电流区模型，该模型对于泄漏电流的检测虽然简单有效，但随着老化、受潮等影响，误差会逐渐变大，使得在线检测有效性不足等技术问题。