1.一种电离层VTEC扰动与雷电相关性分析方法，其特征在于，包括以下步骤：S1、收集GNSS测站双频观测数据，进行电离层VTEC反演计算；

S2、根据各卫星穿刺点经纬度及相应的时间，画图分析测站各卫星不同时刻的轨迹路线；

S3、分析各卫星轨迹路线区域内的雷电数据，选择轨迹内雷电次数多的卫星作为目标卫星，选择雷电与目标卫星轨迹、时间重合的区域作为目标区域；

S4、根据目标卫星的编号提取该卫星的电离层VTEC值；

S5、依次对目标卫星的电离层VTEC值进行求差计算，得到电离层VTEC三次差值；计算电离层VTEC三次差值的公式为：ΔΔΔvteci＝ΔΔvteci-(ΔΔvteci-1+ΔΔvteci+1)/2其中，i表示历元数，i＝1,2…n-1，n表示目标卫星电离层VTEC的历元总数，ΔΔΔvteci表示历元数i时刻的电离层VTEC三次差值，ΔΔvteci历元数i时刻的电离层VTEC二次差值，ΔΔvteci-1是历元数i-1时刻的电离层VTEC二次差值，ΔΔvteci+1是历元数i+1时刻的电离层VTEC二次差值；

S6、计算目标卫星电离层VTEC三次差值序列的均方差，根据误差理论计算目标卫星电离层VTEC三次差值序列的上下限；

S7、分时段统计目标卫星轨迹内的雷电相关参数；

S8、获取东京地磁台站发布的磁情指数Dst和地磁活动指数Kp，分析太阳和地磁活动对电离层的影响；

S9、分时段计算目标卫星电离层VTEC序列与雷电参数序列的相关系数；

S10、根据GNSS测站目标卫星轨迹的穿刺点位置、区域雷电的分布状况、电离层VTEC与雷电参数的相关系数，分析电离层VTEC扰动与雷电的相关性。

2.根据权利要求1所述的一种电离层VTEC扰动与雷电相关性分析方法，其特征在于，步骤S1采用相位平滑伪距方法或者精密单点定位方法，解算测站各卫星电离层VTEC值和各卫星的穿刺点经纬度。

3.根据权利要求1所述的一种电离层VTEC扰动与雷电相关性分析方法，其特征在于，步骤S6中目标卫星电离层VTEC三次差值序列的均方差为std(ΔΔΔ)，设定2std(ΔΔΔ)为目标卫星电离层VTEC三次差值序列的上限，设定-2std(ΔΔΔ)为目标卫星电离层VTEC三次差值序列的下限。

4.根据权利要求1所述的一种电离层VTEC扰动与雷电相关性分析方法，其特征在于，步骤S7中的雷电相关参数包括雷电的发生次数、雷电最大强度和雷电平均强度。

5.根据权利要求1所述的一种电离层VTEC扰动与雷电相关性分析方法，其特征在于，步骤S9目标卫星电离层VTEC序列与雷电参数序列的相关系数满足公式：其中，corr(x,y)表示序列x、y的相关系数，σx表示序列x的方差，σy表示序列y的方差，cov(x,y)表示序列x、y的协方差，x是目标卫星电离层VTEC序列，y是雷电参数序列。