一种闪电三维定位方法技术

本发明专利技术公开了一种闪电三维定位方法，涉及闪电定位领域，S1获取每个测站的脉冲信号；S2读取脉冲信号的振幅，根据寻峰阈值脉冲间隔选取定位脉冲；S3每个测站定位脉冲都与其他测站定位脉冲进行脉冲时间匹配获得脉冲组合；S4将所有匹配好的脉冲组合进行解算得到闪电三维定位结果；S5优化拟合闪电三维定位结果得到优化后的定位结果；采用一种新的脉冲阈值设定和脉冲匹配方案，根据真实物理背景设计合适的剔差方法，在不影响定位精准度的情况下，降低定位计算量，提高定位效率。提高定位效率。提高定位效率。

【技术实现步骤摘要】
一种闪电三维定位方法


[0001]本专利技术涉及闪电定位领域，尤其涉及一种闪电三维定位方法。

技术介绍

[0002]闪电是严重的主要自然灾害之一，可引起森林和油库火灾、造成供电及通讯信息系统故障或损坏，对航天航空、矿山及一些重要而敏感的高技术装备等具有重大威胁。八十年代以后，闪电引起的危害显著增加，特别是与高新技术关系密切的领域，如航空航天、国防、通讯、电力、计算机、电子工业等由于广泛应用对闪电电磁干扰极为敏感的大规模及超大规模集成电路致使遭雷击的几率大大增加；据保守估计，我国每年因雷害造成的直接经济损失超过数亿元，而由此造成的间接经济损失和影响难以估计，由于其成灾迅速而给其研究、预报和防治带来了极大的困难。通过三维闪电定位算法，能够实时计算出闪电的位置与种类(“地闪”或是“云闪”)以及它们的时空演变、放电通道的发展，也可以计算出固定区域内发生闪电的数量，以及放电参数的统计特征从而实现对雷电活动监控并进行预警，从而减小闪电所导致的有害影响，因此发展新的三维闪电定位算法是当前非常迫切的任务。

技术实现思路

[0003]本专利技术的目的就在于为了解决上述问题设计了一种闪电三维定位方法。
[0004]本专利技术通过以下技术方案来实现上述目的：
[0005]一种闪电三维定位方法，包括：
[0006]S1、获取每个测站的脉冲信号；
[0007]S2、读取脉冲信号的振幅，根据寻峰阈值、脉冲间隔选取定位脉冲；
[0008]S3、每个测站定位脉冲都与其他测站定位脉冲进行脉冲时间匹配获得脉冲组合；
[0009]S4、将所有匹配好的脉冲组合进行解算得到初始闪电三维定位结果；
[0010]S5、优化拟合初始闪电三维定位结果。
[0011]本专利技术的有益效果在于：采用一种新的脉冲阈值设定和脉冲匹配方案，根据真实物理背景设计合适的剔差方法，在不影响定位精准度的情况下，降低定位计算量，提高定位效率。
附图说明
[0012]图1是本专利技术一种闪电三维定位方法的流程示意图；
[0013]图2是闪电多站探测示意图。
具体实施方式
[0014]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本专利技术实施
例的组件可以以各种不同的配置来布置和设计。
[0015]因此，以下对在附图中提供的本专利技术的实施例的详细描述并非旨在限制要求保护的本专利技术的范围，而是仅仅表示本专利技术的选定实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0016]应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
[0017]在本专利技术的描述中，需要理解的是，术语“上”、“下”、“内”、“外”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该专利技术产品使用时惯常摆放的方位或位置关系，或者是本领域技术人员惯常理解的方位或位置关系，仅是为了便于描述本专利技术和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。
[0018]此外，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
[0019]在本专利技术的描述中，还需要说明的是，除非另有明确的规定和限定，“设置”、“连接”等术语应做广义理解，例如，“连接”可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接连接，也可以通过中间媒介间接连接，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本专利技术中的具体含义。
[0020]下面结合附图，对本专利技术的具体实施方式进行详细说明。
[0021]如图1、图2所示，一种闪电三维定位方法，包括：
[0022]S1、获取每个测站的脉冲信号。
[0023]S2、读取所有脉冲信号的振幅得到累积概率分布，当累积概率达到99％时，将相应的幅值作为寻峰阈值，设置15μs的时间窗口和脉冲间隔，在时间窗口内寻一个幅值最大的脉冲信号作为定位脉冲，间隔超过15μs的脉冲，只选取幅值高的脉冲信号作为定位脉冲，设置15μs脉冲时间间隔的附加限制，超出间隔范围的2个脉冲，幅值较小的会被舍去，进一步减少计算中使用的匹配组合的数量。
[0024]S3、每个测站定位脉冲与其他测站定位脉冲进行“三角法则”时间窗口内的脉冲匹配获得脉冲组合；三角法则为辐射源到两个测站的时间差小于光在这两个测站所传播的时间。
[0025]S4、将所有匹配好的脉冲组合进行解算得到初始闪电三维定位结果；
[0026]S41、设定位系统各测站的空间位置(x
i
,y
i
,z
i
,t
i
)
T
，i＝0,1,2,3,4，其中，i＝0表示主站，i＝1,2,3,4表示辅站，脉冲信号辐射源的空间位置为(x,y,z,t)
T
；
[0027]S42、计算脉冲信号辐射源到各测站之间的距离为
[0028][0029]简化后为(x0‑
x
i
)x+(y0‑
y
i
)y+(z0‑
z
i
)z＝k
i
+r0·
Δr
i
，其中，r
i
表示辐射源与第i站
之间的距离，Δr
i
表示辐射源第i站与其到主站之间的距离差，
[0030][0031]S43、采用间接法将r0看做已知量，并将简化后的距离写作矩阵式AX＝F，
[0032][0033]S44、采用伪拟法求解得到初始闪电三维定位结果为X＝(A
T
A)
‑1A
T
F，求解时如果所有测站海拔高度为0，且解得闪电三维定位的结果高度为负数时，取绝对值，即将定位结果所在理论双曲线的位置根据水平面镜像移到对称曲线上。
[0034]S5、采用Levenberg
‑
Marquardt(LM)算法优化拟合初始闪电三维定位结果得到优化后的定位结果；
[0035]如果测站数不少于5站时，S4得到闪电三维定位结果不存在模糊解，但得到的结果仍然不是最大似然估计，因为S4中没有考虑线性化方程组的权重，当有测站探测到的数据具有较大误差的时候，算法的定位精度将急剧下降，所以采用Levenberg
‑
Marquardt(LM)算法对定位结果进行再次优化；Levenberg
‑
Marquardt(L本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种闪电三维定位方法，其特征在于，包括：S1、获取每个测站的脉冲信号；S2、读取脉冲信号的振幅，根据寻峰阈值、脉冲间隔选取定位脉冲；S3、每个测站定位脉冲都与其他测站定位脉冲进行脉冲时间匹配获得脉冲组合；S4、将所有匹配好的脉冲组合进行解算得到初始闪电三维定位结果；S5、优化拟合初始闪电三维定位结果。2.根据权利要求1所述的一种闪电三维定位方法，其特征在于，在S2中，读取脉冲信号的振幅得到累积概率分布，当累积概率达到99％时，将相应的幅值作为寻峰阈值；设置15μs的时间窗口和脉冲间隔，在时间窗口内寻一个幅值最大的脉冲信号作为定位脉冲，间隔超过15μs的脉冲，选取幅值高的脉冲信号作为定位脉冲。3.根据权利要求2所述的一种闪电三维定位方法，其特征在于，在S3中，每个测站定位脉冲与其他测站定位脉冲进行“三角法则”时间窗口内的脉冲匹配获得脉冲组合。4.根据权利要求2所述的一种闪电三维定位方法，其特征在于，在S4中包括：S41、设定位系统各测站的空间位置(x
i
,y
i
,z
i
,t
i
)
T
，i＝0,1,2,3,4，其中，i＝0表示主站，i＝1,2,3,4表示辅站，脉冲信号辐射源的空间位置为(x,y,z,t)
T
；S42、计算脉冲信号辐射源到各测站之间的距离为简化后为(x0‑
x
i
)x+(y0‑
y
i
)y+(z0‑

【专利技术属性】
技术研发人员：马子龙，蒋如斌，郄秀书，行鸿彦，刘明远，张鸿波，李宗祥，
申请(专利权)人：中国科学院大气物理研究所，
类型：发明
国别省市：