论文部分内容阅读
近年来,雷电定位系统所提供的数据已成为雷电科学研究和防灾减灾的重要依据。针对不同的使用目的,世界各地建立了多种雷电探测网,探测网由分布于不同位置的多个探测器组成。探测器能接收并记录雷电辐射的电磁波的到达方向和到达时间,定位系统可利用这些观测数据计算出雷电的发生位置和时间等参数。快速地解算出雷电未知参数并进行精度估计是雷电预警、防灾减灾及雷电科学研究中的重要问题,本文在综述了各种雷电探测网特点的基础上,归纳了雷电定位的数学原理,深入研究了雷电目标的无源定位算法和误差分析理论与方法,并评估了武汉雷电探测示范网的性能,主要研究内容和成果包括:(1)对各种雷电定位模型进行了一次全面总结,分别讨论了地闪、云闪定位的数学原理。列出了地闪二维平面定位、地闪椭球面定位、云闪三维时差定位、云闪三维方向定位的基本方程。给出了求算未知数的最大似然估计量的准则和点位误差椭圆的融合方法。建立了雷电定位问题的数学框架,是后续研究的基础。(2)给出了基于泰勒级数展开的迭代算法的计算公式,首次研究了雷电定位迭代算法收敛区和发散区的分布,结果表明雷电定位对初始值的精度要求较宽松。若采用探测网内任意一点作为初始位置,对发生于探测网内外绝大部分区域的目标进行定位时,迭代算法都能收敛,但迭代计算发散的区域会一直存在,即“万能”的初始值是不存在的。此外,把网格搜索算法和粒子群搜索算法引入到雷电定位领域。搜索算法具有较好的稳定性,但计算量较大。进一步提出了采用粒子群算法和迭代算法进行联合定位计算的策略。该联合算法能保证迭代算法的收敛,并大幅减少算法的计算量。(3)采用假设检验的统计方法判断并剔除雷电观测数据中的含粗差数据,以适应函数值作为评价指标,当观测数据不含粗差时适应函数值服从卡方分布。分析了判断粗差是否存在及进一步能否找出含粗差数据的条件,表明只有在雷电探测数据较多时才能进行粗差剔除。利用算例验证了粗差剔除的可行性,通过不断地组合探测数据,对应适应函数值最小的组合通常不含粗差,结果表明剔除了粗差后定位精度大幅提高。(4)研究了不同定位模型对地闪定位误差造成的影响。由于地球为扁椭球体,平面和球面定位模型只适用于覆盖范围较小的探测网。因而地闪定位只有采用椭球面定位模型方可保证精确定位。比较了两套不同大地线长度和方位角计算公式的计算结果,证明公式计算结果对地闪定位应用足够精确。(5)给出了利用蒙特卡罗仿真结果、协方差矩阵或CRB下限公式计算定位误差参数的方法,用数值算例表明三种计算方法所得结果是等价的。利用上述误差分析方法研究了地闪方向定位、到达时间定位及三维云闪到达时间定位的误差分布特征。针对到达时间法的特点,对二维和三维目标分别提出了利用单位圆内接三角形的面积和单位球内接四面体的体积估计点位误差大小的几何模型,这两种误差模型可用于评估目标定位误差、解释误差分布特征并指导探测网站址的选择。(6)评估了武汉示范网的性能,性能指标包括定位误差分布及有效探测范围。研究了武汉示范网地闪定位误差分布特征,误差分布计算结果说明了方向观测数据对该小型探测网的重要性。分析了不同高度云闪目标的定位误差分布,结果表明目标越高,定位误差越小;同时目标离探测网中心的平面距离越远,误差越大。讨论了地球曲率及山脉遮挡对云闪探测率的影响,并计算了不同平面位置处可被定位出的云闪的高度下限。建议未来建立云闪探测网时应提高探测站高度,可使探测范围增大。