搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于正交传播算子的闪电宽带甚高频辐射源定位方法研究

李书磊 邱实 石立华 李云 段艳涛

基于正交传播算子的闪电宽带甚高频辐射源定位方法研究

李书磊, 邱实, 石立华, 李云, 段艳涛
PDF
HTML
导出引用
导出核心图
  • 闪电甚高频辐射源定位技术为闪电放电特征及其物理机制的研究提供了重要手段. 基于空间谱估计理论, 可将正交传播算子方法应用于闪电放电通道时空演变过程的成像. 该方法将阵列数据协方差矩阵进行线性分解形成正交传播算子, 然后以子空间的正交性构造空间谱, 通过空间谱搜索实现辐射源定位. 针对闪电宽带甚高频信号, 采用非相干子空间处理方法将带宽内的有效频点进行平均, 减小噪声干扰. 利用数值仿真分析了该方法的定位性能, 验证该方法定位弱辐射源的有效性, 并与时间反转技术进行了对比. 针对人工触发闪电过程的定位结果表明该方法可以较高的时空分辨率清晰地描绘出闪电通道的基本结构及放电通道的发展过程, 并且其对双源同窗事件的定位能力优于时间反转方法. 该方法对提高宽带甚高频阵列在闪电弱辐射源定位、闪电起始机制的研究中的应用价值具有重要意义.
      通信作者: 邱实, zeustone@yeah.net
    • 基金项目: 国家自然科学基金(批准号: 41675002) 资助的课题.
    [1]

    Shao X M, Holden D N, Rhodes C T 1996 Geophys. Res. Lett. 23 1917

    [2]

    Dong W, Liu X, Yu Y, Zhang Y 2001 Chin. Sci. Bull. 46 1561

    [3]

    Shao X M, Blaine W, Dingus B, Smith D, Ho C, Caffrey M, Graham P, Haynes B, Bowers G, Rassoul H 2018 J. Geophys. Res.: Atmospheres 123 1

    [4]

    Lyu F, Cummer S A, Qin Z, Chen M 2019 J. Geophys. Res.: Atmospheres 124 2994

    [5]

    Shao X M, Krehbiel P R, Thomas R J, Rison W 1995 J. Geophys. Res. Atmospheres 100 2749

    [6]

    Shao X M, Krehbiel P R 1996 J. Geophys. Res.: Atmospheres 101 26641

    [7]

    Ushio T O, Kawasaki Z I, Ohta Y, Matsuura K 1997 Geophys. Res. Lett. 24 2769

    [8]

    董万胜, 刘欣生, 张义军 2001 科学通报 46 427

    Dong W S, Liu X S, Zhang Y J 2001 Chin. Sci. Bull. 46 427

    [9]

    邱实, 杨波, 董万胜, 高太长 2009 气象科学 29 92

    Qiu S, Yang B, Dong W S, Gao T C 2009 J. Meteorol. Sci. 29 92

    [10]

    Sun Z, Qie X, Liu M, Cao D, Wang D 2013 Atmosph. Res. 129-130 58

    [11]

    Stock M G, Akita M, Krehbiel P R, Rison W, Edens H E, Kawasaki Z, Stanley M A 2014 J. Geophys. Res.: Atmospheres 11 9

    [12]

    Kasemir H W 1960 J. Geophys. Res. 65 1873

    [13]

    Mazur V 2002 C. R. Phys. 3 1393

    [14]

    苟学强, 张义军, 李亚珺, 陈明理 2018 物理学报 67 205201

    Gou X Q, Zhang Y J, Li Y J, Chen M L 2018 Acta Phys. Sin. 67 205201

    [15]

    Mazur V, Ruhnke L H, Warner T A, Orville R E 2013 J. Electrostat. 71 763

    [16]

    Wang T, Shi Q, Shi L H, Li Y 2017 IEEE Trans. Electromag. Compatibility 59 1949

    [17]

    Wang T, Shi L, Qiu S, Sun Z, Zhang Q, Duan Y, Liu B 2018 IEEE Access 6 26558

    [18]

    Marcos S, Benidir M 1990 Acoustics, Speech, & Signal Processing, on IEEE International Conference France, Albuquerque, New Mexico,USA, 3-6 April 1990

    [19]

    Marcos S, Marsal A, Benidir M 1994 Proceedings of ICASSP '94. IEEEInternational Conference on Acoustics, Speech and Signal Processing, Adelaide, Australia, 19-22 April,1994.

    [20]

    王永良, 陈辉, 彭应宁, 万群 2004 空间谱估计理论与算法 (北京: 清华大学出版社)

    Wang Y L, Chen H, Peng Y N, Wan Q 2004 Spatial Spectrum Estimation Theory and Algorithms (Beijing: Peking University Press) pp18–30 (in Chinese)

    [21]

    张小飞, 汪飞, 徐大专 2010 阵列信号处理的理论和应用 (北京: 国防工业出版社)

    Zhang X F, Wang F, Xu D Z 2010 Theory and Application of Array Signal Processing (Beijing: National Defense Industry Press) pp118–119 (in Chinese)

    [22]

    Ester M, Kriegel H P, Sander J, Xu X Proceedings of 2nd International Conference on Knowledge Discovery and Data Mining Portland, USA, August 02-04,1996 pp226–231

  • 图 1  定位误差与信噪比的关系

    Fig. 1.  Relationship between locating error and SNR.

    图 2  方位角和仰角的HPW与SNR的关系 (a) HPW示意; (b)方位角; (c)仰角

    Fig. 2.  Relationship between HPW and SNR of azimuth and elevation: (a) Sketch map of HPW; (b) azimuth; (c) elevation.

    图 3  MARCOS系统天线布置示意图

    Fig. 3.  Schematic diagram of MARCOS system antennas layout.

    图 4  Trig230553的快电场波形

    Fig. 4.  Fast electric field waveform of Trig230553.

    图 5  Trig230553放电通道的定位结果 (a) VHF阵列定位结果; (b)高速摄像结果

    Fig. 5.  Locating results of discharge channels of Trig230553: (a) Mapping results of VHF arrays; (b) HSV results

    图 6  第二次直窜先导过程

    Fig. 6.  The second dart-leader process.

    图 7  Trig230553上行正先导成像结果比较 (a)宽带OPM方法定位结果; (b)宽带FDTR方法定位结果

    Fig. 7.  Comparison of mapping results for upward positive leader of Trig230553: (a) Mapping results of OPM method; (b) mapping results of FDTR method

    图 8  宽带OPM方法双源定位结果

    Fig. 8.  Mapping results of broad band OPM method for two sources of Trig230553.

    表 1  低信噪比条件下两种算法的识别概率的比较

    Table 1.  Comparison of recognition probabilities of two algorithms under low SNR.

    SNR
    –15 dB–14 dB
    宽带OPM定位算法65.22%98.76%
    宽带FDTR定位算法65.23%98.68%
    下载: 导出CSV

    表 2  不同入射角度时两种方法的最小可分辨角度

    Table 2.  Minimum distinguishable angles of two methods at different incident angles.

    入射角度(θ, φ)
    (30°,30°)(120°,60°)(220°,45°)(320°,80°)
    宽带OPM
    定位方法
    12°
    宽带FDTR
    定位方法
    13°16°
    下载: 导出CSV
  • [1]

    Shao X M, Holden D N, Rhodes C T 1996 Geophys. Res. Lett. 23 1917

    [2]

    Dong W, Liu X, Yu Y, Zhang Y 2001 Chin. Sci. Bull. 46 1561

    [3]

    Shao X M, Blaine W, Dingus B, Smith D, Ho C, Caffrey M, Graham P, Haynes B, Bowers G, Rassoul H 2018 J. Geophys. Res.: Atmospheres 123 1

    [4]

    Lyu F, Cummer S A, Qin Z, Chen M 2019 J. Geophys. Res.: Atmospheres 124 2994

    [5]

    Shao X M, Krehbiel P R, Thomas R J, Rison W 1995 J. Geophys. Res. Atmospheres 100 2749

    [6]

    Shao X M, Krehbiel P R 1996 J. Geophys. Res.: Atmospheres 101 26641

    [7]

    Ushio T O, Kawasaki Z I, Ohta Y, Matsuura K 1997 Geophys. Res. Lett. 24 2769

    [8]

    董万胜, 刘欣生, 张义军 2001 科学通报 46 427

    Dong W S, Liu X S, Zhang Y J 2001 Chin. Sci. Bull. 46 427

    [9]

    邱实, 杨波, 董万胜, 高太长 2009 气象科学 29 92

    Qiu S, Yang B, Dong W S, Gao T C 2009 J. Meteorol. Sci. 29 92

    [10]

    Sun Z, Qie X, Liu M, Cao D, Wang D 2013 Atmosph. Res. 129-130 58

    [11]

    Stock M G, Akita M, Krehbiel P R, Rison W, Edens H E, Kawasaki Z, Stanley M A 2014 J. Geophys. Res.: Atmospheres 11 9

    [12]

    Kasemir H W 1960 J. Geophys. Res. 65 1873

    [13]

    Mazur V 2002 C. R. Phys. 3 1393

    [14]

    苟学强, 张义军, 李亚珺, 陈明理 2018 物理学报 67 205201

    Gou X Q, Zhang Y J, Li Y J, Chen M L 2018 Acta Phys. Sin. 67 205201

    [15]

    Mazur V, Ruhnke L H, Warner T A, Orville R E 2013 J. Electrostat. 71 763

    [16]

    Wang T, Shi Q, Shi L H, Li Y 2017 IEEE Trans. Electromag. Compatibility 59 1949

    [17]

    Wang T, Shi L, Qiu S, Sun Z, Zhang Q, Duan Y, Liu B 2018 IEEE Access 6 26558

    [18]

    Marcos S, Benidir M 1990 Acoustics, Speech, & Signal Processing, on IEEE International Conference France, Albuquerque, New Mexico,USA, 3-6 April 1990

    [19]

    Marcos S, Marsal A, Benidir M 1994 Proceedings of ICASSP '94. IEEEInternational Conference on Acoustics, Speech and Signal Processing, Adelaide, Australia, 19-22 April,1994.

    [20]

    王永良, 陈辉, 彭应宁, 万群 2004 空间谱估计理论与算法 (北京: 清华大学出版社)

    Wang Y L, Chen H, Peng Y N, Wan Q 2004 Spatial Spectrum Estimation Theory and Algorithms (Beijing: Peking University Press) pp18–30 (in Chinese)

    [21]

    张小飞, 汪飞, 徐大专 2010 阵列信号处理的理论和应用 (北京: 国防工业出版社)

    Zhang X F, Wang F, Xu D Z 2010 Theory and Application of Array Signal Processing (Beijing: National Defense Industry Press) pp118–119 (in Chinese)

    [22]

    Ester M, Kriegel H P, Sander J, Xu X Proceedings of 2nd International Conference on Knowledge Discovery and Data Mining Portland, USA, August 02-04,1996 pp226–231

  • [1] 冯桂力, 郄秀书, 赵阳, 孔祥贞, 张广庶, 张彤, 杨静, 张其林, 王东方. 人工触发闪电电流波形特征参数分析. 物理学报, 2009, 58(9): 6616-6626. doi: 10.7498/aps.58.6616
    [2] 李艳阳, 杨仕娥, 陈永生, 周建朋, 李新利, 卢景霄. 甚高频电容耦合氢等离子体特性研究. 物理学报, 2012, 61(16): 165203. doi: 10.7498/aps.61.165203
    [3] 王彩霞, 郄秀书, 蒋如斌, 杨静. 一次人工触发闪电上行正先导的传输特征. 物理学报, 2012, 61(3): 039203. doi: 10.7498/aps.61.039203
    [4] 张晓丹, 郑新霞, 王光红, 许盛之, 岳强, 林泉, 魏长春, 孙建, 张德坤, 熊绍珍, 耿新华, 赵颖. 单室沉积高效非晶硅/微晶硅叠层太阳电池的研究. 物理学报, 2010, 59(11): 8231-8236. doi: 10.7498/aps.59.8231
    [5] 曹冬杰, 郄秀书, 段树, 宣越建, 王东方. 基于VHF辐射源短基线定位系统对闪电放电过程的研究. 物理学报, 2012, 61(6): 069202. doi: 10.7498/aps.61.069202
    [6] 刘冬霞, 郄秀书, 王志超, 吴学珂, 潘伦湘. 飑线系统中的闪电辐射源分布特征及云内电荷结构讨论. 物理学报, 2013, 62(21): 219201. doi: 10.7498/aps.62.219201
    [7] 郭学军, 卢景霄, 陈永生, 张庆丰, 文书堂, 郑 文, 申陈海, 陈庆东. 甚高频高速沉积微晶硅薄膜的研究. 物理学报, 2008, 57(9): 6002-6006. doi: 10.7498/aps.57.6002
    [8] 王俊, 王涛, 唐成双, 辛煜. 甚高频激发容性耦合氩等离子体的电子能量分布函数的演变. 物理学报, 2016, 65(5): 055203. doi: 10.7498/aps.65.055203
    [9] 杨郁, 唐成双, 赵一帆, 虞一青, 辛煜. 甚高频激发的容性耦合Ar+O2等离子体电负特性研究. 物理学报, 2017, 66(18): 185202. doi: 10.7498/aps.66.185202
    [10] 张晓丹, 赵 颖, 高艳涛, 朱 锋, 魏长春, 孙 建, 王 岩, 耿新华, 熊绍珍. 甚高频等离子体增强化学气相沉积制备微晶硅太阳电池的研究. 物理学报, 2005, 54(4): 1899-1903. doi: 10.7498/aps.54.1899
    [11] 葛 洪, 张晓丹, 岳 强, 赵 颖, 赵 静. 甚高频等离子体增强化学气相沉积大面积平行板电极间真空电势差分布研究. 物理学报, 2008, 57(8): 5105-5110. doi: 10.7498/aps.57.5105
    [12] 丁艳丽, 朱志立, 谷锦华, 史新伟, 杨仕娥, 郜小勇, 陈永生, 卢景霄. 沉积速率对甚高频等离子体增强化学气相沉积制备微晶硅薄膜生长标度行为的影响. 物理学报, 2010, 59(2): 1190-1195. doi: 10.7498/aps.59.1190
    [13] 吴春亚, 赵 颖, 薛俊明, 耿新华, 熊绍珍, 杨恢东. 甚高频等离子体增强化学气相沉积法沉积μc-Si∶H薄膜中氧污染的初步研究. 物理学报, 2003, 52(11): 2865-2869. doi: 10.7498/aps.52.2865
    [14] 高喜, 杨梓强, 侯钧, 亓丽梅, 兰峰, 史宗君, 梁正, 李大治. 具有变态光子带隙结构的相对论Cherenkov辐射源的研究. 物理学报, 2009, 58(2): 1105-1109. doi: 10.7498/aps.58.1105
    [15] 李爽, 王建国, 童长江, 王光强, 陆希成, 王雪锋. 大功率0.34 THz辐射源中慢波结构的优化设计. 物理学报, 2013, 62(12): 120703. doi: 10.7498/aps.62.120703
    [16] 李树, 邓力, 田东风, 李刚. 基于能量密度分布的辐射源粒子空间抽样方法研究. 物理学报, 2014, 63(23): 239501. doi: 10.7498/aps.63.239501
    [17] 施卫, 闫志巾. 雪崩倍增GaAs光电导太赫兹辐射源研究进展. 物理学报, 2015, 64(22): 228702. doi: 10.7498/aps.64.228702
    [18] 傅涛, 欧阳征标. 等离子体填充金属光子晶体Cherenkov辐射源模拟研究. 物理学报, 2016, 65(7): 074208. doi: 10.7498/aps.65.074208
    [19] 李娜, 白亚, 刘鹏. 激光等离子体太赫兹辐射源的频率控制. 物理学报, 2016, 65(11): 110701. doi: 10.7498/aps.65.110701
    [20] 胡耀垓, 赵正予, 项薇, 张援农. 人工电离层洞形态调制及其对短波传播的影响. 物理学报, 2011, 60(9): 099402. doi: 10.7498/aps.60.099402
  • 引用本文:
    Citation:
计量
  • 文章访问数:  194
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-10
  • 修回日期:  2019-06-03
  • 上网日期:  2019-09-11
  • 刊出日期:  2019-08-01

基于正交传播算子的闪电宽带甚高频辐射源定位方法研究

  • 陆军工程大学, 电磁环境效应与光电工程国家级重点实验室, 南京 210007
  • 通信作者: 邱实, zeustone@yeah.net
    基金项目: 国家自然科学基金(批准号: 41675002) 资助的课题.

摘要: 闪电甚高频辐射源定位技术为闪电放电特征及其物理机制的研究提供了重要手段. 基于空间谱估计理论, 可将正交传播算子方法应用于闪电放电通道时空演变过程的成像. 该方法将阵列数据协方差矩阵进行线性分解形成正交传播算子, 然后以子空间的正交性构造空间谱, 通过空间谱搜索实现辐射源定位. 针对闪电宽带甚高频信号, 采用非相干子空间处理方法将带宽内的有效频点进行平均, 减小噪声干扰. 利用数值仿真分析了该方法的定位性能, 验证该方法定位弱辐射源的有效性, 并与时间反转技术进行了对比. 针对人工触发闪电过程的定位结果表明该方法可以较高的时空分辨率清晰地描绘出闪电通道的基本结构及放电通道的发展过程, 并且其对双源同窗事件的定位能力优于时间反转方法. 该方法对提高宽带甚高频阵列在闪电弱辐射源定位、闪电起始机制的研究中的应用价值具有重要意义.

English Abstract

    • 闪电甚高频(VHF)辐射源定位技术可以实现闪电放电通道时空演变过程的高分辨率成像分析, 对闪电放电机理的认识及其电磁环境效应防护具有重要的意义[14].

      自Shao等[5,6]提出宽带干涉仪后, 基于短基线(基线长度通常几米到几十米)宽带VHF定位的各类系统在闪电物理研究中发挥了重要作用. Ushio等[7]采用分段触发的方式对一次空中触发闪电的下行负先导产生的VHF辐射进行观测和特征分析, 但没有探测到上行正先导产生的辐射. 董万胜等[8]利用三阵元正交基线宽带干涉仪阵列基于干涉法对人工触发闪电进行了观测研究, 首次探测到人工触发闪电上行正先导的辐射. 针对阵列天线间的信号时延估计, 邱实等[9]将时域互相关技术应用于宽带VHF辐射源定位. Sun等[10]提出基于小波变换的广义互相关时延估计, 利用基小波在不同尺度下的频谱对互功率谱密度进行加权, 降低了噪声的影响. Stock等[11]设计了连续采集宽带干涉仪, 并利用广义互相关算法提高了弱辐射源的探测能力.

      闪电起始的双向先导理论认为强电场环境下的闪电通道将分别以正、负先导的方式同时双向发展[12,13]. 然而, 正击穿的辐射功率比负击穿至少低一个量级, 传统VHF辐射源定位方法往往导致负击穿掩盖掉同时发展的正击穿. 正先导的发展过程观测和特征分析是澄清闪电发展过程的一个亟待解决的问题[14]. 此外, 反冲先导是否同样存在未被VHF探测到的正先导部分仍尚待观测[15]. 这对定位算法的弱辐射源定位和多源定位能力提出了较高要求. 针对弱辐射源定位, 通过多阵元对辐射信号进行聚焦, 具有较高的空间分辨能力、较高的信号增益及较强的抗干扰能力. Wang等[16,17]基于七阵元“L”型天线阵, 首次采用时间反转技术对人工触发闪电放电通道进行定位, 并提出了频域时间反转(FDTR)成像方法, 证明了该方法在弱辐射源探测中的有效性. 该方法实质上是谱估计理论在空域的扩展形式, 但阵列的角度分辨力受到空域“傅里叶限”的限制, 在一定程度上降低了定位精度, 导致位于一个波束宽度内的多个空间目标不可分辨. 另外, Wang等[17]对该方法的双源定位能力进行了实验验证, 但并未在实际闪电定位中应用. Marcos等[18,19]提出传播算子算法进行波达方向(direct of arrival, DOA)估计, 通过线性运算构造传播算子, 利用子空间的正交性获得空间谱峰, 提高了算法的空间分辨力, 同时避免矩阵特征分解所带来的计算量. 该方法因其计算复杂度低、易于工程实现等优点广泛应用于雷达目标定位、移动通信领域.

      针对时间反转波束宽度较宽、角度分辨力低的缺点, 本文将正交传播算子(OPM)方法应用于闪电宽带VHF辐射源定位, 针对宽带VHF信号, 采用非相干子空间处理方法(ISM)将带宽内的有效频点进行平均, 以减小噪声干扰. 并将定位结果与FDTR算法进行了对比, 表明该算法定位弱辐射源的有效性和同窗双源事件定位的优势.

    • 正交传播算子方法是一种基于空间谱估计理论的经典DOA估计方法. 利用宽带VHF阵列接收数据线性构造特征子空间, 突破了瑞利限的限制, 可实现高分辨率的DOA估计[20,21]. 该方法最初针对窄带DOA估计提出, 而闪电辐射的VHF信号是一种典型的宽带信号, 频率的不同导致了信号子空间的差异, 因此需要通过ISM方法将有效频点加权, 然后进行闪电辐射源DOA估计.

      假设K个闪电辐射信号${s_1}\left( t \right)$, ${s_2}\left( t \right)$, …, ${s_k}\left( t \right)$(k = 1, 2, …, K)分别来自$\left\langle {{\theta _1},{\varphi _1}} \right\rangle $, $\left\langle {{\theta _2},{\varphi _2}} \right\rangle $, …, $\left\langle {{\theta _k},{\varphi _k}} \right\rangle $方向, 其中$\theta $$\varphi $分别表示辐射源的方位角和仰角, 阵元噪声假设为不相关的高斯噪声. 则对于M(要求M > K)元VHF阵列, 输出矢量可表示为

      ${{X}}\left( {{f_j}} \right) = A\left( {{f_j},\left\langle {\theta,\varphi } \right\rangle } \right)S\left( {{f_j}} \right) + N\left( {{f_j}} \right),$

      其中, ${f_j}$表示对VHF信号进行J点离散傅里叶变换后的第$j$个频率点(j = 0, 1, …, J–1). ${{X}}\left( {{f_j}} \right) = {\left[ {{X_1}\left( {{f_j}} \right),{X_2}\left( {{f_j}} \right), \cdots ,{X_M}\left( {{f_j}} \right)} \right]^{\rm{T}}}$表示频率点${f_j}$处的阵列接收数据, $A\left( {{f_j},\left\langle {\theta,\varphi } \right\rangle } \right)$为VHF天线阵列流型, 其形式为

      $A\left( {{f_j},\left\langle {\theta,\varphi } \right\rangle } \right) = \left[\!\!\begin{array}{l} {\rm{e} ^{ - {\rm{j}}2{\rm{\pi }}{f_j}{\tau _1}({\theta _1},{\varphi _1})}}\;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _1}({\theta _2},{\varphi _2})}}\; \cdots \;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _1}({\theta _K},{\varphi _K})}} \\ {\rm{e} ^{ - {\rm{j}}2\rm{p} {f_j}{\tau _2}({\theta _1},{\varphi _1})}}\;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _2}({\theta _2},{\varphi _2})}}\; \cdots \;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _2}({\theta _K},{\varphi _K})}} \\ \quad\quad\quad\vdots \quad\quad\quad\quad\quad\quad\vdots \quad\quad\quad\quad\quad\quad\quad\vdots \\ {\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _M}({\theta _1},{\varphi _1})}}\;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _M}({\theta _2},{\varphi _2})}}\; \cdots \;{\rm{e} ^{ - {\rm{j2\pi }}{f_j}{\tau _M}({\theta _K},{\varphi _K})}} \end{array}\!\!\!\right],$

      其中, ${\tau _m}(\theta,\varphi ) = \dfrac{{{x_m}\cos \theta \cos \varphi + {y_m}\sin \theta \cos \varphi }}{c}$表示第m个VHF天线相对参考点的时延, $({x_m},{y_m})$表示其二维坐标, c为光速.

      假设闪电辐射源信号相互独立, 则天线阵列流型矩阵A是列满秩的, 则矩阵AK行线性无关, 其他MK行可由这K行线性表示. 因此可将天线阵列流形矩阵分解为

      ${{A}} = \left[ \begin{array}{l} {{{A}}_1}^{K \times K} \\ {{{A}}_2}^{\left( {M - K} \right) \times K} \end{array} \right].$

      由于矩阵${{{A}}_1}$为非奇异矩阵, 因此矩阵${{{A}}_2}$可以用矩阵${{{A}}_1}$的线性表达, 即

      ${{{A}}_2} = {{{P}}^{\rm{H}}}{{{A}}_1},$

      其中, 矩阵${{P}} \in {C^{K \times \left( {M - K} \right)}}$称为传播算子. 令${{{Q}}^{\rm{H}}} = \left[ {{{{P}}^{\rm{H}}}, - {{{I}}_{M - K}}} \right]$, 则有

      ${{{Q}}^{\rm{H}}}{{A}} = 0,$

      其中, ${{{I}}_{M - K}}$0分别为$M - K$阶单位矩阵和$\left( {M - K} \right) \times K$阶零矩阵. 上式表明矩阵Q的列向量和VHF天线阵列方向矢量正交, 则矩阵Q的列向量所张成的子空间即为噪声子空间. 将矩阵Q正交化即为正交传播算子${Q_0}$:

      ${Q_0} = {{Q}}{\left( {{{{Q}}^{\rm{H}}}{{Q}}} \right)^{ - 1/2}}.$

      定义窄带OPM定位算法的空间谱如下:

      $P\left( {{f_j},\left\langle {\theta ,\varphi } \right\rangle } \right) = \frac{1}{{{a^{\rm{H}}}\left( {{f_j},\left\langle {\theta ,\varphi } \right\rangle } \right){Q_0}{Q_0}^{\rm{H}}a\left( {{f_j},\left\langle {\theta ,\varphi } \right\rangle } \right)}}.$

      最后, 针对闪电宽带VHF信号, 根据ISM算法的原理, 将频带内有效频点进行平均, 从而得到宽带VHF信号的OPM算法空间谱:

      $P\left( {\left\langle {\theta ,\varphi } \right\rangle } \right) = \frac{1}{J}\sum\limits_{j = 1}^J {\frac{1}{{{a^{\rm{H}}}\left( {{f_j},\left\langle {\theta ,\varphi } \right\rangle } \right){Q_0}{Q_0}^{\rm{H}}a\left( {{f_j},\left\langle {\theta ,\varphi } \right\rangle } \right)}}} .$

      在实际应用中, 由于VHF天线阵列的流型矩阵是未知的, 因此一般从阵列快拍数据协方差矩阵来估计${Q_0}$, 基本步骤为:

      1)利用VHF天线阵列接收数据计算协方差矩阵R, 表达式为${{R}}\left( {{f_j}} \right) = {{X}}\left( {{f_j}} \right){{{X}}^{\rm{H}}}\left( {{f_j}} \right)$;

      2)对矩阵R进行分块, ${{R}} = \left[ {array}{l} {G^{M \times K}} \\ {H^{M \times \left( {M - K} \right)}} \\ {array} \right]$;

      3)由于VHF天线接收数据中噪声的存在, 传播算子的估计由最小化问题得到$J\left( P \right) \!=\! \left\| {H - GP} \right\|_F^2$, 此问题的解为$\widehat P = {\left( {{G^{\rm{H}}}G} \right)^{ - 1}}{G^{\rm{H}}}H$, 并进一步得到${\widehat Q^{\rm{H}}} = \left[ {{{\widehat P}^{\rm{H}}}, - {I_{M - K}}} \right]$;

      4)根据(8)式构造宽带VHF空间谱, 搜索空间谱获得闪电辐射源DOA估计结果; 实际上, 通过上述步骤构造的“倒谱”是利用子空间正交性使得原来在辐射源位置零值的谱变成趋于无限大, 而非辐射源点处的值为相对有限值, 这就将辐射源在区域中与非辐射源位置上的谱幅度区别开来; 而实际由于噪声、阵列互耦、辐射源相干等因素, 使得空间谱中出现非源点的伪峰; 为了实现辐射源同窗双源定位, 需要首先判断是否存在双源; 实测数据分析表明, 搜索次峰值与主峰值的比大于0.5时, 可以认为空间中同时存在两个辐射源, 并分别进行定位;

      5)质量控制方法, 首先, 采用密度聚类分析[22]的方法将离散的伪点去除. 设置一组邻域参数(簇的最小点个数和临近半径, 此处分别设置为20和0.03)判断辐射源点是否紧密相连, 从而去除离散伪点; 然后, 使用能量比(ER)作为辅助判据, 进一步滤除聚成一簇的伪点, 其表达式如(9)式所示, 即信号方向的能量与所有可能方向的总能量的比值,

      $ER = \frac{{P\left( {\theta ,\varphi } \right)}}{{\sum {P\left( {{\theta _i},{\varphi _i}} \right)} }}.$

    • 为了分析算法性能, 通过数值仿真研究了不同信噪比时算法的定位误差、以及空间谱半峰值宽度、角度分辨力等表征双源定位性能的参数, 并与FDTR算法进行了对比.

      参考实际宽带VHF信号波形特征, 仿真时采用正弦高斯调制信号作为辐射源信号, 其数学表达式如下:

      $s\left( t \right) = {A_0}\sin \left( {2{\rm{\pi }}{f_0}t} \right){{\rm{e}}^{ - 4{\rm{\pi }}{{\left( {\frac{{t - {\tau _1}}}{{{\tau _2}}}} \right)}^2}}},$

      式中, ${A_0}$表示信号幅度, 实际处理时进行归一化处理; ${f_0}$表示信号的中心频率, 此处设置为100 MHz; ${\tau _1}$${\tau _2}$用于调节信号的带宽, 参考VHF天线实际采集信号带宽, 此处均设置为20 ns.

      天线阵列的设置型式, 参考实验室野外观测试验的设置方式, 采用7天线“L”型平面阵列, 阵元间距设置为9 m, 详细情况参见4.1节. 仿真时辐射源的方位角和仰角事先指定, 然后根据辐射源与各天线之间的空间坐标关系获得各天线接收信号的延迟, 最后利用定位算法计算辐射源的位置.

    • 设置单一辐射源S的位置为$\theta = {120.25^ \circ }$, $\varphi = {30.50^ \circ }$. 为了定量评估定位结果的准确性, 定义定位误差为DOA估计结果与辐射源实际位置的球面距离, 以角度表示. 此处规定定位误差小于1°时, 即视为算法正确识别目标. 针对不同信噪比(SNR)分别进行10000次独立重复试验, 平均定位误差与信噪比的关系如图1所示.

      图  1  定位误差与信噪比的关系

      Figure 1.  Relationship between locating error and SNR.

      图1可得, 仅对单一辐射源而言, 宽带OPM定位算法和宽带FDTR算法的定位性能基本一致, 二者均能在信噪比大于–12 dB时对辐射源进行有效定位, 且其定位误差均在0.3°以下, 在一定程度上反映了算法的鲁棒性. 当信噪比小于–12 dB时, 两种算法仅能以一定的概率正确识别目标, 识别概率如表1所列. 随信噪比下降, 识别概率相应降低, 同一信噪比下两种算法的识别概率相近.

      SNR
      –15 dB–14 dB
      宽带OPM定位算法65.22%98.76%
      宽带FDTR定位算法65.23%98.68%

      表 1  低信噪比条件下两种算法的识别概率的比较

      Table 1.  Comparison of recognition probabilities of two algorithms under low SNR.

    • 闪电发展过程中多分支通道同时发展以及闪电起始等双向先导事件, 增加了双源同窗事件发生的概率.

      为了定量说明宽带OPM辐射源定位方法对于双辐射源的定位性能, 分别定义空间谱半峰值宽度(HPW)和角度分辨力. 其中, 空间谱HPW为空间谱峰值半功率点的宽度, 主要反映空间谱的尖锐程度, 如图2(a)所示. 图2(b)图2(c)分别为方位角和仰角所对应的HPW与SNR的关系.

      图  2  方位角和仰角的HPW与SNR的关系 (a) HPW示意; (b)方位角; (c)仰角

      Figure 2.  Relationship between HPW and SNR of azimuth and elevation: (a) Sketch map of HPW; (b) azimuth; (c) elevation.

      图2可得, 在信噪比大于–12 dB时, 宽带OPM定位方法的半峰值宽度比FDTR方法小, 表明其谱峰更加尖锐, 旁瓣相对平缓, 能量较为集中. 并且随信噪比的增大, 其HPW呈减小趋势, 而FDTR方法的HPW几乎不随SNR而变化. 当信噪比小于–12 dB时, 两种方法的HPW都随着信噪比的减小而增大, 表明算法的分辨力下降, 此时两种方法的HPW几乎相同.

      角度分辨力则是指定位算法分辨两个相近辐射源的能力. 在进行空间谱搜索时, 若相近的两个辐射源可分辨, 那么其对应的主峰与次峰要大于最小可分辨变角度. 设置两个等幅不相干辐射源, SNR设为0dB. 分别利用两种算法对不同入射方向的辐射源进行DOA估计, 其最小可分辨角度如表2所列. 其中, 最小可分辨变角度通过球面距离表示.

      入射角度(θ, φ)
      (30°,30°)(120°,60°)(220°,45°)(320°,80°)
      宽带OPM
      定位方法
      12°
      宽带FDTR
      定位方法
      13°16°

      表 2  不同入射角度时两种方法的最小可分辨角度

      Table 2.  Minimum distinguishable angles of two methods at different incident angles.

      表2可得, 宽带OPM方法的最小可分辨角度小于宽带FDTR方法, 表明基于子空间正交构造空间谱的定位方法具有更好的角度分辨能力, 从而具有更好的双源定位能力和定位精度, 表明OPM方法在双源定位方面的优势.

    • 自2016年起, SLOT(Jiangsu Lightning Observation Team)团队持续在江苏苏北地区开展人工触发闪电试验, 该地区属江苏省内雷电多发区, 且视野开阔, 具备闪电的长期观测条件. 试验场包括两个闪电观测站点和试验站点, 其中观测站点A处布置了一套多天线辐射源连续观测系统(multiple-antennas radiation continuous obser vation system, MARCOS), 采用7个天线呈“L”型平面阵列, 阵元间距设置为9 m, 阵列设置示意图如图3所示. 同时利用8通道LeCroy高速数字存储示波器采集阵元的VHF信号, 采样率为500 MHz, 存储深度为250 Mpts/ch, 可连续记录500 ms闪电数据. 该系统通过增加阵元数量、增大动态接收范围、采用连续记录模式, 提高了对弱辐射源的探测能力. 此外, 该站点还布置有同步记录电磁场的快、慢天线、磁场天线等, 以及进行光学观测的高速摄像等.

      图  3  MARCOS系统天线布置示意图

      Figure 3.  Schematic diagram of MARCOS system antennas layout.

    • 2017年8月1日23时左右, 雷暴云自东南向西北运动, 于23时05分53秒采用经典触发方式成功触发闪电(记为Trig230553). 触发时地面大气电场为–13.5 kV/m(大气电学符号规定), 据此判断此次触发为负极性触发. 触发闪电的快电场波形记录如图4所示, 0时刻代表记录信号的触发时刻. 此次触发闪电共有5次直窜先导-回击过程, 整个放电过程持续时间约600 ms.

      图  4  Trig230553的快电场波形

      Figure 4.  Fast electric field waveform of Trig230553.

      利用宽带OPM闪电VHF辐射源定位算法, 对此次触发闪电的放电通道进行了定位, 其中滑动窗的窗长设置为512个采样点(1.024 μs), 重叠率为75%. 整个记录时段的闪电成像结果如图5(a)所示, 并将其与同站设置的高速摄像拍摄的结果(图5(b))进行了对比. 两者在闪电通道形状上具有较高的一致性, 但是高速摄像受限于视野范围, 没有获得该次闪电其他的发展过程, 而基于VHF系统的成像可以对高速摄像视野以外的通道有较好的描述. 此次触发闪电的上行正先导起始于A点, 且放电通道分支较少, 当上行正先导发展到仰角为45°左右时出现明显分支通道, 之后仅在正先导发展后期可以观察到两条细小分支通道.

      图  5  Trig230553放电通道的定位结果 (a) VHF阵列定位结果; (b)高速摄像结果

      Figure 5.  Locating results of discharge channels of Trig230553: (a) Mapping results of VHF arrays; (b) HSV results

      在上行正先导发展末期, 在正先导主通道和分支通道末端分别引发了数次反冲先导过程, 该阶段持续时间大约为100 ms. 大约210 ms时, 引发了三次沿正先导主通道反向传输的直窜先导, 其持续时间分别为3, 1, 0.7 ms, 到达地面后均引发回击, 其中第二次直窜先导过程如图6所示. 这三次直窜先导过程均起始于相同的放电区域B处, 沿之前上行正先导开辟的通道向地面方向发展, 定位通道连续、清晰, 且没有其他分支通道.

      图  6  第二次直窜先导过程

      Figure 6.  The second dart-leader process.

    • 为了比较宽带OPM闪电定位方法和FDTR方法对弱辐射源的定位性能, 选取该触发闪电上行正先导阶段的数据, 分别利用两种方法进行通道成像, 结果如图7所示.

      图  7  Trig230553上行正先导成像结果比较 (a)宽带OPM方法定位结果; (b)宽带FDTR方法定位结果

      Figure 7.  Comparison of mapping results for upward positive leader of Trig230553: (a) Mapping results of OPM method; (b) mapping results of FDTR method

      图7可得, 两种方法均可以清晰地分辨出闪电发展通道的结构, 且具有较高的一致性, 对通道发展过程的描述基本类似. 需要注意的是, 由于两种方法空间谱能量集中的程度不同, 为了便于比较, 所以在滤除伪点时只采用CR阈值法, 即两种方法设定相同阈值(此处设置为0.3). 针对上行正先导阶段, 对定位点数进行统计, OPM方法和FDTR方法的定位辐射源数分别为27015和24446. 宽带OPM算法成功解算的定位点数量更多, 这主要是由于宽带OPM方法的角度分辨力高于FDTR方法, 对于多通道同时发展的情形, 其成像结果的细节表现更优.

    • 闪电发展过程中, 多分支通道同时发展时使得VHF天线接收信号发生交错, 因此滑动窗口时可能存在双源同窗的情况. 传统干涉法或互相关方法通过设置高重叠率滑动窗可在一定程度上提高时间分辨率, 但对于时间上相互交错的多辐射源信号无法处理. 宽带OPM定位方法在空域对信号进行处理, 其空间谱表示辐射源在空间各个方向上的能量分布. 在搜索空间谱最优解时, 按照次峰和主峰能量比大于0.5时判定双源同时存在, 否则为单一辐射源. Trig230553发展到–5 ms时, 上行正先导出现明显分支(图7中以P1和P2标识), 采用OPM方法对5 ms时间内的两分支进行了双源自动识别定位, 结果如图8所示.

      图  8  宽带OPM方法双源定位结果

      Figure 8.  Mapping results of broad band OPM method for two sources of Trig230553.

      图8反映了该时间内仰角随时间的变化, P1和P2分别对应于图7中的上行正先导的分支. 红色菱形表示每个窗口空间谱主峰所对应的定位结果, 黄色正方形表示次峰所对应的结果, 两者为同一时间窗对应的辐射源, 如图8中局部放大图所示. 在5 ms时间段内, 宽带OPM方法可定位到有效双源的窗口数为37. 由于估计偏保守, 实际同时存在双源的窗口个数可能更多. 自VHF定位技术应用闪电探测以来, 均只能实现同一时间窗口单个辐射源的定位, 图8给出的定位结果首次实现了对同一时间窗的双源同时定位. 双辐射源的定位结果证实了闪电各个分支上的一些击穿事件存在同时发展的现象, 这对于揭示闪电放电发生发展的物理本质具有重要的参考价值和意义, 尤其对于闪电起始等双向先导发展模式的研究以及模型的建立提供了新的工具和思路.

      利用双源定位算法, 对此次触发闪电的上行正先导阶段进行了双源定位, 定位到有效点数为32943. 但通道形状基本与单源定位结果基本类似, 这主要是由于此次触发闪电双源定位的通道单源基本也可以交替发展的方式进行定位. 另外, 此处双源定位只能给出幅度比小于2的两个辐射源的位置, 需要增加天线阵列数目等方法进一步提高物理孔径, 增强弱辐射源的定位能力.

      由于OPM方法的角度分辨力优于FDTR方法, 具有较小的最小分辨角度, 旁瓣抑制能力更强, 信源方向的空间谱更为突出, 从而更适于空间相距较近的分支通道同时发展时的定位.

    • 本文将基于空间谱估计理论的OPM方法应用于闪电宽带VHF辐射源定位, 实现了闪电放电通道时空演变过程的成像. 该方法首先通过VHF阵列接收数据协方差矩阵的线性分解形成正交化传播算子, 并利用子空间正交构造空间谱, 然后进行空间谱搜索确定辐射源的位置. 针对宽带VHF信号, 采用ISM方法将带宽内的有效频点进行平均, 减小噪声干扰.

      为了验证算法的有效性, 对OPM定位方法和FDTR方法的定位误差、空间谱半峰值宽度、角度分辨力进行了数值仿真. 结果表明, 针对单一辐射源, 宽带OPM定位算法的定位误差与FDTR算法相近, 但OPM算法的角度分辨力优于FDTR算法, 对于分辨双辐射源性能更优. 利用宽带OPM定位算法分析了一次人工触发闪电放电通道的时空发展过程, 结果表明该算法可以较高的时空分辨率清晰地描绘出闪电通道的基本结构及放电通道的发展过程, 并且OPM算法对双辐射源的定位能力有较大优势. 下一步将利用该算法双源定位的优势, 进一步开展对闪电起始先导发展过程的定位研究, 揭示正、负击穿同时存在时的发展规律.

参考文献 (22)

目录

    /

    返回文章
    返回