本发明涉及一种基于光学集成电场传感器的雷电能量频域分布统计方法，属于波形分析技术领域。
背景技术：
高压架空输电线路是电力系统中重要组成部分，其分布很广、纵横交错、绵延数百里乃至上千公里，所以极易受雷击。统计资料显示，我国高压架空输电线路由雷击引起的跳闸次数占运行总跳闸次数的40％-70％。2012年，南方电网110kv及以上电压等级输电线路雷击跳闸1894次，占线路总跳闸次数的65.4％。在多雷、山区以及土壤电阻率较高的地区，雷击架空线引起输电的事故率则更高。当雷击架空输电线路时可能引起线路开关跳闸造成停电事故，还会造成电气设备损坏，甚至导致电力系统瓦解等恶性事故，造成重大的经济损失，严重影响电力系统运行的稳定性、可靠性、安全性和经济性。
我国电网近年来的雷击闪络呈上升趋势,其原因是：1.雷电活动加剧；2.新增电网提速，输电线路增长，杆塔增高，易受雷击；3.因缺乏对雷电活动及其分布规律的了解，防雷措施事半功倍。因此在气候变化异常和电网迅猛发展的今天,研究雷电过电压能量的频域分布是十分迫切的任务。
技术实现要素：
本发明的目的是提出一种基于光学集成电场传感器的雷电能量频域分布统计方法，该方法基于光学集成电场传感器的测量数据，通过滤波、过电压提取、解耦等方法，得到三相独立的雷电过电压波形，经归一化和fft后即可得到雷电过电压能量的频域分布。
本发明提出的基于光学集成电场传感器的雷电能量频域分布统计方法，包括以下步骤：
(1)将三台光学集成电场传感器分别安装在变电站中待测量线路的正下方，使三台光学集成电场传感器的离地高度相同，并使三台光学集成电场传感器的极性朝向一致，分别对三台光学集成电场传感器通光，以在线监测待测量线路上雷电过电压的波形；
(2)当待测量线路有直击雷或者感应雷发生时，三台光学集成电场传感器分别测量得到电场波形，该电场信号由雷电过电压电场波形和三相工频电场波形叠加而成。利用小波变化或平均滤波方法，对测量得到的电场波形进行滤波处理，得到滤波后的电场波形，并对测量得到的电场波形分别乘以相应光学集成电场传感器的转化因数，得到测量电压波形。其中，光学集成电场传感器的转化因数为测量电场是1v/m时所对应的实际电压值；
(3)将上述步骤(2)的测量电压波形，与由变电站中三台光学集成电场传感器，在未发生雷击时对待测量线路测量得到同相位下的工频电压波形相减，得到待测量线路的互有耦合三相雷电过电压波形；
(4)对上述互有耦合三相雷电过电压波形进行解耦，设互有耦合三相雷电过电压波形为ua、ub、uc，解耦后的三相雷电过电压波形为ua、ub、uc，建立如下解耦方程：
其中，qij为解耦系数，i＝1，2，3，j＝1，2，3，qij∈[0,100]；
(5)对上述解耦后的三相雷电过电压波形ua、ub、uc进行归一化处理，然后对归一化处理结果进行快速傅里叶变换(fft)，得到此雷电过电压能量的频域分布；
(6)当三台光学集成电场传感器采集到大量的雷电数据后，重复步骤(2)-步骤(5)，得到大量雷电过电压能量的频域分布，对大量雷电过电压能量的频域分布进行叠加，得到雷电能量的频域分布统计结果。
本发明提出的基于光学集成电场传感器的雷电能量频域分布统计方法，其优点是：
1、本发明方法使用的光学集成电场传感器，具有安装方便，便于测量的特点。
2、本发明方法较为清楚的阐述了基于光学集成电场传感器的雷电过电压能量的频域分布统计方法，步骤简单，容易得到准确的雷电波形参数在频域内的分布情况。
具体实施方式
本发明提出的基于光学集成电场传感器的雷电能量频域分布统计方法，包括以下步骤：
(1)将三台光学集成电场传感器(本发明的一个实施例中，使用的三相光学集成电场传感器由清华大学生产，可参见公开号为cn102854403b的中国专利申请)分别安装在变电站中待测量线路的正下方，使三台光学集成电场传感器的离地高度相同，并使三台光学集成电场传感器的极性朝向一致，分别对三台光学集成电场传感器通光，以在线监测待测量线路上雷电过电压的波形；
(2)当待测量线路有直击雷或者感应雷发生时，三台光学集成电场传感器分别测量得到电场波形，该电场信号由雷电过电压电场波形和三相工频电场波形叠加而成。利用小波变化或平均滤波方法，对测量得到的电场波形进行滤波处理，得到滤波后的电场波形，并对测量得到的电场波形分别乘以相应光学集成电场传感器的转化因数，得到测量电压波形。其中，光学集成电场传感器的转化因数为测量电场是1v/m时所对应的实际电压值；
(3)将上述步骤(2)的测量电压波形，与由变电站中三台光学集成电场传感器，在未发生雷击时对待测量线路测量得到同相位下的工频电压波形相减，得到待测量线路的互有耦合三相雷电过电压波形；
(4)对上述互有耦合三相雷电过电压波形进行解耦，设互有耦合三相雷电过电压波形为ua、ub、uc，解耦后的三相雷电过电压波形为ua、ub、uc，建立如下解耦方程：
其中，qij为解耦系数，i＝1，2，3，j＝1，2，3，qij∈[0,100]；
(5)对上述解耦后的三相雷电过电压波形ua、ub、uc进行归一化处理，然后对归一化处理结果使用快速傅里叶变换(fft)，得到此雷电过电压能量的频域分布；
(6)当三台光学集成电场传感器采集到大量的雷电数据后，重复步骤(2)-步骤(5)，得到大量雷电过电压能量的频域分布，对大量雷电过电压能量的频域分布进行叠加，得到雷电能量的频域分布统计结果。