一种输电线路雷击定位和雷电流反演方法
【专利摘要】本发明涉及一种输电线路的雷击定位反演方法，属于高电压试验设备及测量【技术领域】。本方法通过查询雷电流在线监测系统监测数据，获取击中杆塔的雷电流经杆塔和输电线路分流后的波形参数；通过查询雷电定位系统监测数据，获取雷击时间、雷击区域、雷电流幅值、雷电回击次数数据；结合理论关系表和理论曲线，考虑雷电定位系统的误差，得出雷击点的位置区间，并通过仿真计算反演雷电流波形的特征参数。本方法定位精度高，有利于线路雷击损坏的快速查找和修复，节省了输电线路雷击故障查找和线路修复时间，能够还原真实的雷电流波形，为雷电成因和雷电流特征的科学研究提供依据，为避雷措施和雷电过电压防护的研究提供数据。
【专利说明】一种输电线路雷击定位和雷电流反演方法
【技术领域】
[0001 ] 本发明涉及一种输电线路雷击定位和雷电流反演方法，属于高电压试验设备及测量【技术领域】。
【背景技术】
[0002]据不完全统计，我国每年因雷击造成人员伤亡达3000?4000人，财产损失更是达到惊人的200亿元，雷电灾害已成为破坏性日趋严重的气象灾害之一。雷电带来的电磁暂态现象属于电磁污染范围，而电磁污染是世界公认的继水质污染、大气污染、噪声污染之后的第四大污染。随着工农业生产和人民生活的现代化，对铁路、航空、金融、电力、电信、电视、网络等服务系统和设施可靠性及服务质量要求也越来越高，人类活动对这些公用事业的依赖性也越来越大，对雷电防护提出了更高的要求，因此对雷电流准确测量显得越来越重要。
[0003]从电力系统的角度出发，雷电已经成为了电力系统故障的最主要原因，无论是线路、变电站还是用电设备，如何使它们能够更好地抵御雷电和雷电感应产生的电磁暂态过程的冲击是电力系统雷电防护中最为关心的问题，这需要建立在对雷电的放电过程和雷电流的各项参数有充分认识和了解的基础上。
[0004]目前对雷电流特征参数的研究大都依靠模拟和仿真的方法，而对自然界的雷电流特征参数，包括波形、幅值、上升时间、持续时间等都缺少实际数据。如何能够获得真实的雷电流数据，这对于雷电放电过程的认识和雷电流及其空间场的研究都十分有意义。而获得自然界雷电参数最直接的方法就是测量，以现在最常见的人工火箭引雷或者自然雷直击高大物体为例，都是将测量装置的探头置于雷电流入地通道上，而将主机置于通道附近或者室内进行测量和记录的。但由于雷电流幅值往往很高，基本为千安培量级，且上升沿往往只有几百纳秒到几微秒，所以选择什么样的测量探头能够满足大的动态测量范围和宽的频率响应是最为关心的问题，罗氏线圈是目前比较通用的选择，它可以满足上述两点的要求，而且它的性能已经为很多实验所证实。
[0005]因为雷电是一种随机程度很高的自然现象，所以要大量的测量雷电流参数，最好的方法就是大量的安装测量装置，但是受到各种实际条件的制约，这样的做法显然是不合适的，所以在测量装置安装地点的选择上首先要考虑以往的雷电活动规律，特别是雷击点的分布规律和雷击次数的统计规律。根据不同的雷电活动强度可以将测量地点划分为多雷击区域和少雷击区域，比如通信基站的信号塔，广播电视发射塔这些独立、高大的建筑物，以及容易遭受雷击的输电线路的局部和某些变电站，这些都属于高雷击区域。所以在这些地点安装相对较复杂，数据记录能力较强的测量装置。而在其它的少雷击区域可以选择记录功能较单一(只记录幅值)，结构较为简单，造价较低的测量装置。这样可以充分利用可用的测量资源，最大限度地测量和收集自然界雷电流的参数信息。
[0006]当雷电击中杆塔或输电线路时，雷击点所处的线路中，会有幅值很大的雷电流通过，可能造成线路损坏。同时，雷击点临近的绝缘子处可能会发生闪络，发生闪络后，雷电流流过闪络区域，会发出热量烧蚀绝缘子，影响绝缘子的强度和绝缘特性。如能通过技术手段对雷击杆塔或输电线路进行定位，可以实现对由雷击损坏的线路进行快速查找和修复 。
【发明内容】
[0007]本发明目的是提出一种输电线路雷击定位和雷电流反演方法，依据上述雷电定位系统和雷电流在线监测系统给出的测量数据，计算雷击线路故障发生的雷击位置和该雷击的雷电流波形参数，为电力系统由雷击造成的线路损坏的快速查找和修复提供位置依据，并通过雷电流的反演，得出雷电流波形参数，为雷电成因和雷电流特征的科学研究提供依据，为避雷措施和雷电过电压防护的研究提供数据。
[0008]本发明提出的输电线路雷击定位和雷电流反演方法，包括以下步骤:
[0009](I)从雷电流在线监测系统中获取雷电流在线监测系统所安装的杆塔的位置Pp雷击时间^、以及击中杆塔的雷电流经杆塔和输电线路分流后通过在线监测系统中雷电流传感器的雷电流特征参数，雷电流特征参数包括雷电流的电流幅值I1、雷电流波形、上升时间tupl和半波时间thl，并从雷电流波形中，得到波峰次数N1 ；
[0010](2)从雷电定位系统中获取监测数据，监测数据包括同时刻同区域的雷击时间t2、雷电流击中区域P2、雷电流幅值I2和雷电回击次数N2 ；
[0011](3)将上述雷电流在线监测系统中获取的雷击时间、雷电流在线监测系统所安装的杆塔的P1、雷电流特征参数中的电流幅值I1和波形中的波峰次数N1,与雷电定位系统监测数据中的雷击时间&、雷电流击中区域P2、雷电流幅值I2、雷电回击次数N2进行比对，若同时满足h与t2相差I毫秒以内、P1位于区域P2内、N1与N2相同，且I1小于I2，则判定雷电流在线监测系统记录的数据和雷电定位系统记录的数据由同一雷电引起；
[0012](4)根据雷电击中输电线路时，输电线路中，地线上的雷电流幅值1、上升时间tup和半波时间th与雷击点和观测点之间的杆塔数量及杆塔间线路的段数的理论关系表和理论关系曲线，将上述雷电流幅值I1与雷电流幅值I2的比值I1A2、上升时间tupl和半波时间thl代入上述理论关系表和理论关系曲线中，检索到雷电流在线监测系统所在位置与雷电流击中输电线路的位置之间的杆塔数量G和线路的段数D ；
[0013](5)通过对雷击线路的仿真计算，得到雷电流沿线路传播时，经过G个杆塔和D段线路广生幅值裳减后的雷电流幅值占裳减如的雷电流幅值的百分比Q，将上述雷电流在线监测系统测量的雷电流幅值I1除以Q，得出衰减前距雷电流击中线路位置最近处杆塔的雷电流幅值I。:10=Ii/Q ；
[0014](6)根据雷电击中输电线路杆塔之间的地线时，雷电流沿地线向两侧分流比例与雷击点距两侧杆塔的距离之间的理论关系表，将上述雷电流幅值IjP I2的比例IcZI2代入该理论关系表中，检索到雷电流击中位置与最近杆塔的距离L，若L为0，则判定雷电流直接击中杆塔，若L不为0，则判定雷电流击中杆塔之间的线路；
[0015](7)利用下式，计算得到雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S:S=GXA+L，其中，A为从电力系统的输电线路数据中得到雷电流在线监测系统所在位置处的杆塔之间的平均距离；
[0016](8)对上述雷电流幅值的测量值I2分别增加50%和减少50%后，分别重复步骤(4)至步骤(7)，从得到的雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S中，将最大值Smax和最小值Smin分别定义为雷电流击中输电线路位置与雷电流在线监测系统之间距离的上限和下限；[0017](9)根据上述步骤(8)的最小值Smin和最大值Smax，得到雷电流击中输电线路位置所在的线段[SMIN，Smax],记该线段的中点为Sav，再根据上述雷电流在线监测系统中获取的雷电流在线监测系统所安装的杆塔的位置P1，得到雷电流击中输电线路的位置Ptl =P0=P^Sav,将Ptl与雷电定位系统监测的雷击区间P2进行比较，若Ptl处于区域P2内，则将Ptl作为雷击线路的位置，用于雷击故障修复和线路维护，若Ptl处于区域P2外，则结束定位计算；
[0018](10)根据上述雷电流击中输电线路位置，对击中杆塔的雷电流的电流幅值、电流波形、上升时间、半波时间等雷电流特征参数的反演，其过程为:
[0019](10-1)根据从电力系统的数据手册中获取的杆塔类型、杆塔参数、线路类型和线路参数，建立电路仿真模型；
[0020]( 10-2)将上述雷电流击中输电线路位置Ptl代入上述电路仿真模型中，并在电路仿真模型中的Ptl处设置标准双指数雷电流发生模型，进行仿真计算；
[0021](10-3)在仿真计算过程中，分别调节双指数雷电流发生模型中的雷电流幅值、上升时间和半波时间，使得在电路仿真模型中的模拟雷击发生时，电路仿真模型中的雷电流在线监测系统处，通过地线的雷电流幅值、上升时间和半波时间分别与上述步骤(1)中的雷电流的电流幅值I1、上升时间tupl和半波时间thl之间的相差在5%以内，仿真计算得到的模拟雷电流的幅值、上升时间和半波时间，即为反演出的雷电流的幅值、上升时间和半波时间。
[0022]本发明提出的输电线路雷击定位和雷电流反演方法，其优点是:
[0023](I)本发明综合了雷电定位系统和雷电流在线监测系统给出的测量数据，计算雷击线路的位置和该雷击的雷电流波形参数，为电力系统由雷击造成的线路损坏的快速查找和修复提供位置依据，其定位精度高于单独由雷电定位系统或雷电流在线监测系统给出的雷电定位结果，定位误差小于一至两个杆塔的距离，对线路雷击损坏的快速查询，具有参考意义；
[0024](2)本发明方法中的雷击定位过程，考虑了雷电定位系统的测量误差，通过雷电定位系统和雷电流在线监测系统互相比对，减小误差，得出雷电定位的距离区间，保证了数据的有效性和可靠性；
[0025](3)通过本发明方法中的雷电流反演，可以还原真实雷电流的波形、电流幅值、上升时间和半波时间，为雷电成因和雷电流特征的科学研究提供依据，为避雷措施和雷电过电压的防护的研究提供数据。
【专利附图】
【附图说明】
[0026]图1为雷电流击中杆塔顶部或杆塔之间示意图。
[0027]图2为雷电流击中杆塔之间线路不同位置示意图。
[0028]图3为输电线路中地线上的雷电流幅值I与雷击点和观测点间的杆塔数量的关系O
[0029]图4为输电线路中地线上的雷电流上升时间tup与雷击点和观测点间的杆塔数量的关系。[0030]图5为输电线路中地线上的雷电流半波时间th与雷击点和观测点间的杆塔数量的关系。
[0031]图6为雷电流沿地线向一侧分流比例与雷击点相对该侧杆塔距离的关系。
[0032]图7为考虑雷电定位系统误差后根据雷电流幅值判断雷击点的位置示意图。
【具体实施方式】
[0033]本发明提出的输电线路雷击定位和雷电流反演方法，包括以下步骤:
[0034](I)从雷电流在线监测系统中获取雷电流在线监测系统所安装的杆塔的位置Pp雷击时间^、以及击中杆塔的雷电流经杆塔和输电线路分流后通过在线监测系统中雷电流传感器的雷电流特征参数，雷电流特征参数包括雷电流的电流幅值I1、雷电流波形、上升时间tupl和半波时间thl，并从雷电流波形中，得到波峰次数N1 ；
[0035](2)从雷电定位系统中获取监测数据，监测数据包括同时刻同区域的雷击时间t2、雷电流击中区域P2、雷电流幅值I2和雷电回击次数N2 ；
[0036](3)将上述雷电流在线监测系统中获取的雷击时间h、雷电流在线监测系统所安装的杆塔的P1、雷电流特征参数中的电流幅值I1和波形中的波峰次数N1,与雷电定位系统监测数据中的雷击时间&、雷电流击中区域P2、雷电流幅值I2、雷电回击次数N2进行比对，若同时满足h与t2相差I毫秒以内、P1位于区域P2内、N1与N2相同，且I1小于I2，则判定雷电流在线监测系统记录的数据和雷电定位系统记录的数据由同一雷电引起；
[0037](4)根据雷电击中输电线路时，输电线路中，地线上的雷电流幅值1、上升时间tup和半波时间th与雷击点和观测点之间的杆塔数量及杆塔间线路的段数的理论关系表和理论关系曲线，将上述雷电流幅值I1与雷电流幅值I2的比值I1A2、上升时间tupl和半波时间thl代入上述理论关系表和理论关系曲线中，检索到雷电流在线监测系统所在位置与雷电流击中输电线路的位置之间的杆塔数量G和线路的段数D ；
[0038](5)通过对雷击线路的仿真计算，得到雷电流沿线路传播时，经过G个杆塔和D段线路广生幅值裳减后的雷电流幅值占裳减如的雷电流幅值的百分比Q，将上述雷电流在线监测系统测量的雷电流幅值I1除以Q，得出衰减前距雷电流击中线路位置最近处杆塔的雷电流幅值I。:10=Ii/Q ；
[0039](6)根据雷电击中输电线路杆塔之间的地线时，雷电流沿地线向两侧分流比例与雷击点距两侧杆塔的距离之间的理论关系表，将上述雷电流幅值IjP I2的比例IcZI2代入该理论关系表中，检索到雷电流击中位置与最近杆塔的距离L，若L为0，则判定雷电流直接击中杆塔，若L不为0，则判定雷电流击中杆塔之间的线路；
[0040](7)利用下式，计算得到雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S:S=GXA+L，其中，A为从电力系统的输电线路数据中得到雷电流在线监测系统所在位置处的杆塔之间的平均距离；
[0041](8)对上述雷电流幅值的测量值I2分别增加50%和减少50%后，分别重复步骤(4)至步骤(7)，从得到的雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S中，将最大值Smax和最小值Smin分别定义为雷电流击中输电线路位置与雷电流在线监测系统之间距离的上限和下限；
[0042](9)根据上述步骤(8)的最小值Smin和最大值Smax，得到雷电流击中输电线路位置所在的线段[SMIN，Smax],记该线段的中点为Sav，再根据上述雷电流在线监测系统中获取的雷电流在线监测系统所安装的杆塔的位置P1，得到雷电流击中输电线路的位置Ptl =P0=P^Sav,将Ptl与雷电定位系统监测的雷击区间P2进行比较，若Ptl处于区域P2内，则将Ptl作为雷击线路的位置，用于雷击故障修复和线路维护，若Ptl处于区域P2外，则结束定位计算；
[0043](10)根据上述雷电流击中输电线路位置，对击中杆塔的雷电流的电流幅值、电流波形、上升时间、半波时间等雷电流特征参数的反演，其过程为:
[0044](10-1)根据从电力系统的数据手册中获取的杆塔类型、杆塔参数、线路类型和线路参数，建立电路仿真模型；
[0045](10-2)将上述雷电流击中输电线路位置P。代入上述电路仿真模型中，并在电路仿真模型中的Ptl处设置标准双指数雷电流发生模型，进行仿真计算；
[0046](10-3)在仿真计算过程中，分别调节双指数雷电流发生模型中的雷电流幅值、上升时间和半波时间，使得在电路仿真模型中的模拟雷击发生时，电路仿真模型中的雷电流在线监测系统处，通过地线的雷电流幅值、上升时间和半波时间分别与上述步骤(1)中的雷电流的电流幅值I1、上升时间tupl和半波时间thl之间的相差在5%以内，仿真计算得到的模拟雷电流的幅值、上升时间和半波时间，即为反演出的雷电流的幅值、上升时间和半波时间。
[0047]本发明定位方法中，输电线路中地线上的雷电流幅值、上升时间和半波时间与雷击点与观测点的距离的关系为，随着雷电流经过杆塔数量的增加，雷电波波形上升时间逐渐变短，幅值逐渐变小，半波时间变长，具体计算结果如下:
[0048]使用国家电网公司出版的《输电线路设计标准手册》作为仿真的设计标准，其中杆塔采用二分裂式酒杯塔，整个 输电线路包含11个典型杆塔，线路末端各设置为1000km (即线路为无穷长)。档距均为400m。三相导线采用LGJ-50，水平排列，均距7.5m，两塔头同三相导线的垂直距离为5.7m，三相导线距离地面20m，导线为二分裂，分裂间距0.4572m，导线半径0.02m,电阻率0.03984 Ω /km,垂驰12m ;接地线半径0.0057m,电阻率为2.8645 Ω /km,垂驰7m。
[0049]雷电流击中杆塔顶部或杆塔之间的示意图和雷击杆塔之间线路不同位置的示意图分别如图1和图2所示。
[0050]雷电流采用标准雷电流波形进行计算，得出当雷击塔顶或杆塔之间的线路时，如图1所示，通过仿真计算，得到通过临近几个杆塔的雷电流仿真波形，其总体规律如表1所
/Jn ο
[0051]表1雷电波形经过杆塔后衰减数据
[0052]
【权利要求】
1.一种输电线路雷击定位和雷电流反演方法，其特征在于该方法包括以下步骤: (1)从雷电流在线监测系统中获取雷电流在线监测系统所安装的杆塔的位置P1、雷击时间、以及击中杆塔的雷电流经杆塔和输电线路分流后通过在线监测系统中雷电流传感器的雷电流特征参数，雷电流特征参数包括雷电流的电流幅值I1、雷电流波形、上升时间tupl和半波时间thl，并从雷电流波形中，得到波峰次数N1 ； (2)从雷电定位系统中获取监测数据，监测数据包括同时刻同区域的雷击时间t2、雷电流击中区域P2、雷电流幅值I2和雷电回击次数N2 ； (3)将上述雷电流在线监测系统中获取的雷击时间h、雷电流在线监测系统所安装的杆塔的位置P1、雷电流特征参数中的电流幅值I1和波形中的波峰次数N1，与雷电定位系统监测数据中的雷击时间&、雷电流击中区域P2、雷电流幅值I2、雷电回击次数N2进行比对，若同时满足h与t2相差I毫秒以内、P1位于区域P2内、N1与N2相同，且I1小于I2，则判定雷电流在线监测系统记录的数据和雷电定位系统记录的数据由同一雷电引起； (4)根据雷电击中输电线路时，输电线路中，地线上的雷电流幅值1、上升时间tup和半波时间th与雷击点和观测点之间的杆塔数量及杆塔间线路的段数的理论关系表和理论关系曲线，将上述雷电流幅值I1与雷电流幅值I2的比值I1A2、上升时间tupl和半波时间thlR入上述理论关系表和理论关系曲线中，检索到雷电流在线监测系统所在位置与雷电流击中输电线路的位置之间的杆塔数量G和线路的段数D ； (5)通过对雷击线路的仿真计算，得到雷电流沿线路传播时，经过G个杆塔和D段线路产生幅值衰减后的雷电流幅值占衰减前的雷电流幅值的百分比Q，将上述雷电流在线监测系统测量的雷电流幅值 I1除以Q，得出衰减前距雷电流击中线路位置最近处杆塔的雷电流幅值 10 =10=VQ ； (6)根据雷电击中输电线路杆塔之间的地线时，雷电流沿地线向两侧分流比例与雷击点距两侧杆塔的距离之间的理论关系表，将上述雷电流幅值Itl和I2的比例IcZI2代入该理论关系表中，检索到雷电流击中位置与最近杆塔的距离L，若L为O，则判定雷电流直接击中杆塔，若L不为O，则判定雷电流击中杆塔之间的线路； (7)利用下式，计算得到雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S:S=GXA+L，其中，A为从电力系统的输电线路数据中得到雷电流在线监测系统所在位置处的杆塔之间的平均距离； (8)对上述雷电流幅值的测量值I2分别增加50%和减少50%后，分别重复步骤(4)至步骤(7)，从得到的雷电流击中输电线路位置与雷电流在线监测系统位置之间的距离S中，将最大值Smax和最小值Smin分别定义为雷电流击中输电线路位置与雷电流在线监测系统之间距离的上限和下限； (9)根据上述步骤(8)的最小值Smin和最大值Smax，得到雷电流击中输电线路位置所在的线段[SMIN，Smax]，记该线段的中点为SAV，再根据上述雷电流在线监测系统所安装的杆塔的位置P1，得到雷电流击中输电线路的位置Ptl =Pci=PJSav，将Ptl与雷电定位系统监测的雷击区间P2进行比较，若Ptl处于区域P2内，则将Ptl作为雷击线路的位置，用于雷击故障修复和线路维护，若Ptl处于区域P2外，则结束定位计算； (10)根据上述雷电流击中输电线路位置，对击中杆塔的雷电流的电流幅值、电流波形、上升时间、半波时间等雷电流特征参数的反演，其过程为:(10-1)根据从电力系统的数据手册中获取的杆塔类型、杆塔参数、线路类型和线路参数，建立电路仿真模型； (10-2)将上述雷电流击中输电线路位置Ptl代入上述电路仿真模型中，并在电路仿真模型中的Ptl处设置标准双指数雷电流发生模型，进行仿真计算； (10-3)在仿真计算过程中，分别调节双指数雷电流发生模型中的雷电流幅值、上升时间和半波时间，使得在电路仿真模型中的模拟雷击发生时，电路仿真模型中的雷电流在线监测系统处，通过地线的雷电流幅值、上升时间和半波时间分别与上述步骤(1)中的雷电流的电流幅值I1、上升时间tupl和半波时间thl之间的相差在5%以内，仿真计算得到的模拟雷电流的幅值、上升时间和半波时间，即为反演出的雷电流的幅值、上升时间和半波时间。
【文档编号】G01R19/00GK103543387SQ201310483424
【公开日】2014年1月29日 申请日期:2013年10月16日 优先权日:2013年10月16日
【发明者】余占清, 王运超, 王泽众, 李 雨, 彭珑, 赵媛, 陈水明, 何金良, 曾嵘, 端木林楠, 牛晓民, 王萍 申请人:清华大学, 国家电网公司, 华北电力科学研究院有限责任公司, 华北电网有限公司