留言板

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

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

基于FIR滤波的磁静期电离层-磁层磁场信号分离

郭斐 邱耀东 王正涛

郭斐, 邱耀东, 王正涛. 基于FIR滤波的磁静期电离层-磁层磁场信号分离[J]. 武汉大学学报 ● 信息科学版, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
引用本文: 郭斐, 邱耀东, 王正涛. 基于FIR滤波的磁静期电离层-磁层磁场信号分离[J]. 武汉大学学报 ● 信息科学版, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
GUO Fei, QIU Yaodong, WANG Zhengtao. Signal Separation of Ionospheric and Magnetospheric Magnetic Field During the Magnetic Quiet Period Based on FIR Filter[J]. Geomatics and Information Science of Wuhan University, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
Citation: GUO Fei, QIU Yaodong, WANG Zhengtao. Signal Separation of Ionospheric and Magnetospheric Magnetic Field During the Magnetic Quiet Period Based on FIR Filter[J]. Geomatics and Information Science of Wuhan University, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157

基于FIR滤波的磁静期电离层-磁层磁场信号分离

doi: 10.13203/j.whugis20180157
基金项目: 

国家杰出青-科学基金 41525014

国家自然科学基金 41774019

国家自然科学基金 41274032

国家自然科学基金 41474018

国家自然科学基金 41429401

详细信息
    作者简介:

    郭斐, 硕士, 主要从事卫星地磁数据处理、卫星摄动力与地应力研究。fguo_sgg@whu.edu.cn

    通讯作者: 王正涛, 博士, 教授。ztwang@sgg.whu.edu.cn
  • 中图分类号: P223

Signal Separation of Ionospheric and Magnetospheric Magnetic Field During the Magnetic Quiet Period Based on FIR Filter

Funds: 

The National Science Fund for Distinguished Youth Scholars 41525014

the National Natural Science Foundation of China 41774019

the National Natural Science Foundation of China 41274032

the National Natural Science Foundation of China 41474018

the National Natural Science Foundation of China 41429401

More Information
    Author Bio:

    GUO Fei, master, specializes in satellite magnetic data processing, satellite perturbation forces and geostress. E-mail:fguo_sgg@whu.edu.cn

    Corresponding author: WANG Zhengtao, PhD, professor. E-mail:ztwang@sgg.whu.edu.cn
  • 摘要: Swarm是欧空局第4个"地球探测者"任务,其主要目的是确定地球磁场及其时空演化特征。利用Swarm构建地球各圈层磁场模型的关键问题之一是对不同场源的地磁信号实施准确分离。以Swarm磁静期观测数据为例,采用FIR沿轨高通滤波对外源场信号进行处理,分析并给出了滤波器窗口长度的确定方法,利用欧空局2016年发布的MIO_SHA_2D和MMA_SHA_2C电离层-磁层磁场模型验证了该方法的可行性。结果表明:FIR对外源场长波信号的滤波结果与模型改正结果的一致性吻合较好,在给定通带波长1 200 km、阻带波长3 000 km的情况下,将滤波器长度设置为87~129阶,可同时顾及水平方向和径向方向的滤波精度。
  • 图  1  基于卫星磁测数据提取磁静期剩余磁场的预处理流程

    Figure  1.  Preprocessing Process of Extracting Magnetic Remaining Magnetic Field Based on Satellite Magnetic Data

    图  2  剩余磁场XYZ分量强度

    Figure  2.  Remaining Magnetic Field X, Y, Z Components Intensity

    图  3  剩余磁场Z分量的幅频响应特性

    Figure  3.  The Amplitude-Frequency Response of the Z Component of the Residual Magnetic Field

    图  4  径向分量Z的外源场改正强度

    Figure  4.  The Corrected Intensity of External Field of Radial Component Z

    图  5  改正后的径向残差ΔZ大小

    Figure  5.  Corrected Radial Residual ΔZ

    图  6  径向残差ΔZ的幅频响应特性

    Figure  6.  The Amplitude-Frequency Response of Radial Residual ΔZ

    表  1  基于均方根误差(RMSE)的滤波器长度对剩余磁场长波信号改正的精度影响/nT

    Table  1.   Effect of Filter Length Based on Root Mean Square Error (RMSE) on the Correction Precision of Long Wave Signal of Remaining Magnetic Field/nT

    N RMSE最小值 RMSE最大值 RMSE平均值
    X Y Z X Y Z X Y Z
    55 0.461 0.562 0.430 4.555 3.550 4.001 1.808 1.715 1.537
    87 0.438 0.504 0.498 4.498 3.944 3.543 1.825 1.767 1.397
    97 0.442 0.501 0.500 4.500 3.920 3.472 1.825 1.761 1.391
    129 0.452 0.562 0.497 4.503 3.899 3.214 1.827 1.780 1.347
    161 0.461 0.675 0.498 4.500 4.004 3.108 1.831 1.859 1.315
    193 0.472 0.779 0.511 4.512 4.246 2.738 1.846 1.939 1.246
    225 0.497 0.462 0.525 4.594 4.463 2.033 1.862 1.855 1.192
    下载: 导出CSV
  • [1] 徐文耀.地磁学[M].北京:地震出版社, 2003

    Xu Wenyao. Geomagnetism[M]. Beijing:Seismological Press, 2003
    [2] 常宜峰.卫星磁测数据处理与地磁场模型反演理论与方法研究[D].郑州: 信息工程大学, 2015 http://cdmd.cnki.com.cn/Article/CDMD-90005-1016058210.htm

    Chang Yifeng. Research on Satellite Geomagnetic Data Process and Geomagnetic Model Recovery Theory and Method[D]. Zhengzhou: Information Engineering University, 2015 http://cdmd.cnki.com.cn/Article/CDMD-90005-1016058210.htm
    [3] 安振昌, 谭东海, 王月华, 等.利用MAGSAT卫星数据研究青藏高原及其邻近地区磁异常[J].测绘学报, 1996, 25(3):221-225 doi:  10.3321/j.issn:1001-1595.1996.03.011

    An Zhenchang, Tan Donghai, Wang Yuehua, et al.Mangnetic Anomalies for Qinghai-Xizang (Tibetan) Plateau and Adjacent Region Derived from MAGSAT Data[J]. Acta Geodaetica et Cartographica Sinica, 1996, 25(3):221-225 doi:  10.3321/j.issn:1001-1595.1996.03.011
    [4] Ou J, Gillet N, Du A. Geomagnetic Observatory Monthly Series, 1930 to 2010:Empirical Analysis and Unmodeled Signal Estimation[J]. Earth, Planets and Space, 2015, 67(1):1-20
    [5] 杜劲松, 陈超.基于卫星磁测数据的全球岩石圈磁场建模进展与展望[J].地球物理学进展, 2015, 30(3):1017-1033 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxjz201503005

    Du Jinsong, Chen Chao. Progress and Outlook in Global Lithospheric Magnetic Field Modelling by Satellite Magnetic Measurements[J]. Progress in Geophysics, 2015, 30(3):1017-1033 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxjz201503005
    [6] 曲春凯, 李斐, 杨轩, 等.利用零相位Kaiser窗滤波器改善MEX多普勒数据定轨精度[J].武汉大学学报·信息科学版, 2018, 43(7):1071-1077 http://ch.whu.edu.cn/CN/abstract/abstract6155.shtml

    Qu Chunkai, Li Fei, Yang Xuan, et al.Improving MEX Orbit Determination Accuracy by Using Zero-Phase Kaiser Window Filter[J]. Geomatics and Information Science of Wuhan University, 2018, 43(7):1071-1077 http://ch.whu.edu.cn/CN/abstract/abstract6155.shtml
    [7] Thébault E, Vigneron P, Maus S, et al. Swarm SCARF Dedicated Lithospheric Field Inversion Chain[J]. Earth, Planets and Space, 2013, 65(11):1257-1270 doi:  10.5047/eps.2013.07.008
    [8] Sabaka T J, Olsen N, Langel R A. A Comprehensive Model of the Quiet-time, Near-Earth Magnetic Field:Phase 3[J]. Geophysical Journal International, 2002, 151(1):32-68 doi:  10.1046/j.1365-246X.2002.01774.x
    [9] Sabaka T J, Olsen N, Purucker M E. Extending Comprehensive Models of the Éarth's Magnetic Field with Ørsted and CHAMP Data[J]. Geophysical Journal International, 2004, 159(2):521-547 doi:  10.1111/gji.2004.159.issue-2
    [10] Sabaka T J, Olsen N. Enhancing Comprehensive Inversions Using the Swarm Constellation[J]. Earth, Planets and Space, 2006, 58(4):371-395 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=J-STAGE_595992
    [11] Lühr H, Maus S. Solar Cycle Dependence of Quiet-time Magnetospheric Currents and a Model of Their Near-Earth Magnetic Fields[J]. Earth, Planets and Space, 2010, 62(10):843-848 doi:  10.5047/eps.2010.07.012
    [12] Olsen N, Sabaka T J, Lowes F. New Parameterization of External and Induced Fields in Geomagnetic field Modeling, and a Candidate Model for IGRF 2005[J]. Earth, Planets and Space, 2005, 57(12):1141-1149 doi:  10.1186/BF03351897
    [13] 冯彦, 安振昌, 毛飞, 等.磁静日期间中国大陆地区外源弱磁场研究[J].空间科学学报, 2012, 32(3):321-335 http://www.cnki.com.cn/Article/CJFDTOTAL-KJKB201203004.htm

    Feng yan, An Zhenchang, Mao Fei, et al. Study on External Weak Magnetic Field During Magnetic Quiet Periods over Chinese Mainland[J]. Chinese Journal of Space Science, 2012, 32(3):321-335 http://www.cnki.com.cn/Article/CJFDTOTAL-KJKB201203004.htm
    [14] Olsen N, Friischristensen E, Floberghagen R, et al. The Swarm Satellite Constellation Application and Research Facility (SCARF) and Swarm Data Produ-cts[J]. Earth, Planets and Space, 2013, 65(11):1189-1200 doi:  10.5047/eps.2013.07.001
    [15] Chulliat A, Vigneron P, Thébault E, et al.Swarm SCARF Dedicated Ionospheric Field Inversion Chain[J]. Earth, Planets and Space, 2013, 65(11):1257-1270 doi:  10.5047/eps.2013.07.008
    [16] Thébault E. Global Lithospheric Magnetic Field Modelling by Successive Regional Analysis[J]. Earth, Planets and Space, 2006, 58(4):485-495 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=J-STAGE_595999
    [17] 周波阳, 罗志才, 宁津生, 等.航空矢量重力测量中有限冲激响应低通数字滤波器的设计[J].武汉大学学报·信息科学报, 2015, 40(6):772-778 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201506012

    Zhou Boyang, Luo Zhicai, Ning Jinsheng, et al. Design of FIR Lowpass Digital filters for Airborne Vector Gravimetry[J]. Geomatics and Information Science of Wuhan University, 2015, 40(6):772-778 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201506012
    [18] 孙中苗, 夏哲仁, 李迎春.L&R航空重力仪摆杆尺度因子的确定与分析[J].武汉大学学报·信息科学版, 2002, 27(4):367-371 http://ch.whu.edu.cn/CN/abstract/abstract4974.shtml

    Sun Zhongmiao, Xia Zheren, Li Yingchun. Determination and Analysis of the Beam Scale Factor for the LaCoste&Romberg Airborne Gravimeter[J]. Geomatics and Information Science of Wuhan University, 2002, 27(4):367-371 http://ch.whu.edu.cn/CN/abstract/abstract4974.shtml
    [19] Maus S, Yin F, Lühr H, et al. Resolution of Direction of Oceanic Magnetic Lineations by the Sixth-generation Lithospheric Magnetic Field Model from CHAMP Satellite Magnetic Measurements[J]. Geo-chemistry Geophysics Geosystems, 2008, 9(7):1335-1346
    [20] Kother L, Hammer M D, Finlay C C, et al. An Equivalent Source Method for Modelling the Global Lithospheric Magnetic Field[J]. Geophysical Journal International, 2015, 203(1):553-566 doi:  10.1093/gji/ggv317
    [21] 黄珹.利用LAGEOS激光测距资料精确测定地球自转参数[D].上海: 中国科学院上海天文台, 1985

    Huang Cheng. Using LAGEOS Laser Ranging Data to Accurately Measure the Earth Rotation Parameters[D]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Sciences, 1985
    [22] 严豪健.天基GPS气象学与反演技术[M].北京:中国科学技术出版社, 2007

    Yan Haojian. Space-based GPS Meteorology and Inversion Technique[M]. Beijing:China Science and Technology Press, 2007
    [23] 陈怀琛. MATLAB及其在理工课程中的应用指南[M].西安:西安电子科技大学出版社, 2007

    Chen Huaichen. MATLAB and Its Application in the Course of Science and Engineering[M]. Xi'an:Xidian University Press, 2007
    [24] 孙中苗, 夏哲仁.航空重力测量中低通滤波器参数的优选[J].测绘科学技术学报, 2002(1):11-14, 18 doi:  10.3969/j.issn.1673-6338.2002.01.004

    Sun Zhongmiao, Xiazheren. Optimal Determination of the Lowpass Filter Parameters in Airborne Gravimetry[J]. Journal of the PLA Institute of Surveying & Mapping, 2002(1):11-14, 18 doi:  10.3969/j.issn.1673-6338.2002.01.004
  • [1] 朱永兴, 谭述森, 任夏, 贾小林.  GNSS全球广播电离层模型精度分析 . 武汉大学学报 ● 信息科学版, 2020, 45(5): 768-775. doi: 10.13203/j.whugis20180439
    [2] 李夫鹏, 王正涛, 超能芳, 冯建迪, 张兵兵, 田坤俊, 韩亚坤.  利用Swarm星群探测亚马逊流域2015-2016年干旱事件 . 武汉大学学报 ● 信息科学版, 2020, 45(4): 595-603. doi: 10.13203/j.whugis20180273
    [3] 汤俊, 毛文飞.  多尺度ARMA残差修正模型震前电离层TEC异常探测 . 武汉大学学报 ● 信息科学版, 2019, 44(6): 791-798. doi: 10.13203/j.whugis20170321
    [4] 周仁宇, 胡志刚, 苏牡丹, 李军正, 李鹏博, 赵齐乐.  北斗全球系统广播电离层模型性能初步评估 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1457-1464. doi: 10.13203/j.whugis20170428
    [5] 赵齐乐, 杨超, 耿涛, 马福建, 张强.  NTCM广播电离层模型性能评估 . 武汉大学学报 ● 信息科学版, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
    [6] 王成, 毛大智, 施闯, 章红平.  全球电离层模型的分布式并行解算 . 武汉大学学报 ● 信息科学版, 2018, 43(8): 1207-1213. doi: 10.13203/j.whugis20160274
    [7] 施闯, 王成, 张涛.  基于超算的全球电离层模型快速并行解算 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2093-2098. doi: 10.13203/j.whugis20180114
    [8] 王成, 施闯, 章红平.  全球电离层模型的OpenMP多线程并行解算 . 武汉大学学报 ● 信息科学版, 2018, 43(6): 813-819. doi: 10.13203/j.whugis20150223
    [9] 宋福成, 时爽爽, 冯建迪.  基于Chapman函数的电离层TEC同化模型构建 . 武汉大学学报 ● 信息科学版, 2016, 41(6): 784-790. doi: 10.13203/j.whugis20150101
    [10] 宁津生, 王正涛, 超能芳.  国际新一代卫星重力探测计划研究现状与进展 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 1-8. doi: 10.13203/j.whugis20150732
    [11] 王成, 王解先, 段兵兵.  附有国际参考电离层约束的全球电离层模型 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1340-1346.
    [12] 阮仁桂, 吴显兵, 冯来平.  单频精密单点定位观测模型和电离层处理方法比较 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1023-1028.
    [13] 汤俊, 姚宜斌, 陈鹏, 张顺.  利用EMD方法改进电离层TEC预报模型 . 武汉大学学报 ● 信息科学版, 2013, 38(4): 408-411.
    [14] 施闯, 耿长江, 章红平, 唐卫明.  基于EOF的实时三维电离层模型精度分析 . 武汉大学学报 ● 信息科学版, 2010, 35(10): 1143-1146.
    [15] 章红平, 施闯, 唐卫明.  地基GPS区域电离层多项式模型与硬件延迟统一解算分析 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 805-809.
    [16] 柳景斌, 王泽民, 王海军, 章红平.  利用球冠谐分析方法和GPS数据建立中国区域电离层TEC模型 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 792-795.
    [17] 李征航, 陈锴, 刘万科, 王文丽.  GNSS电离层延迟模型的数学统一与方法扩展 . 武汉大学学报 ● 信息科学版, 2007, 32(8): 699-703.
    [18] 李征航, 陈锴, 刘万科, 黄欢.  顾及f~3项的电离层延迟模型 . 武汉大学学报 ● 信息科学版, 2007, 32(2): 139-143.
    [19] 程鹏飞, 李夕银, 马建平.  两种GPS测定电离层电子密度模型的探讨 . 武汉大学学报 ● 信息科学版, 2001, 26(6): 524-528.
    [20] 张小红, 李征航, 蔡昌盛.  用双频GPS观测值建立小区域电离层延迟模型研究 . 武汉大学学报 ● 信息科学版, 2001, 26(2): 140-143,159.
  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  637
  • HTML全文浏览量:  67
  • PDF下载量:  61
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-26
  • 刊出日期:  2019-06-05

基于FIR滤波的磁静期电离层-磁层磁场信号分离

doi: 10.13203/j.whugis20180157
    基金项目:

    国家杰出青-科学基金 41525014

    国家自然科学基金 41774019

    国家自然科学基金 41274032

    国家自然科学基金 41474018

    国家自然科学基金 41429401

    作者简介:

    郭斐, 硕士, 主要从事卫星地磁数据处理、卫星摄动力与地应力研究。fguo_sgg@whu.edu.cn

    通讯作者: 王正涛, 博士, 教授。ztwang@sgg.whu.edu.cn
  • 中图分类号: P223

摘要: Swarm是欧空局第4个"地球探测者"任务,其主要目的是确定地球磁场及其时空演化特征。利用Swarm构建地球各圈层磁场模型的关键问题之一是对不同场源的地磁信号实施准确分离。以Swarm磁静期观测数据为例,采用FIR沿轨高通滤波对外源场信号进行处理,分析并给出了滤波器窗口长度的确定方法,利用欧空局2016年发布的MIO_SHA_2D和MMA_SHA_2C电离层-磁层磁场模型验证了该方法的可行性。结果表明:FIR对外源场长波信号的滤波结果与模型改正结果的一致性吻合较好,在给定通带波长1 200 km、阻带波长3 000 km的情况下,将滤波器长度设置为87~129阶,可同时顾及水平方向和径向方向的滤波精度。

English Abstract

郭斐, 邱耀东, 王正涛. 基于FIR滤波的磁静期电离层-磁层磁场信号分离[J]. 武汉大学学报 ● 信息科学版, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
引用本文: 郭斐, 邱耀东, 王正涛. 基于FIR滤波的磁静期电离层-磁层磁场信号分离[J]. 武汉大学学报 ● 信息科学版, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
GUO Fei, QIU Yaodong, WANG Zhengtao. Signal Separation of Ionospheric and Magnetospheric Magnetic Field During the Magnetic Quiet Period Based on FIR Filter[J]. Geomatics and Information Science of Wuhan University, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
Citation: GUO Fei, QIU Yaodong, WANG Zhengtao. Signal Separation of Ionospheric and Magnetospheric Magnetic Field During the Magnetic Quiet Period Based on FIR Filter[J]. Geomatics and Information Science of Wuhan University, 2019, 44(6): 844-850. doi: 10.13203/j.whugis20180157
  • 地球外源场约占地球总磁场的1%[1],它起源于近地空间电流体系,这些电流体系主要分布在电离层-磁层,行星际空间电流变化虽有影响但贡献甚微[2]。外源场最明显的特征是具有快速的时空变化分布,且空间分布具有全球性[1]。由于外源场携带着丰富而复杂的地球空间电磁环境信息,因此,研究如何从地球各圈层磁场中提取外源场信号,分析外源场的时空变化特征,对于认知空间电磁环境的变化规律,理解日地能量耦合机制,实时准确预报空间天气变化具有重要的科学意义和潜在的应用价值。

    地磁场模型参数化的关键问题之一是对地球各圈层磁场的贡献实施有效分离。对于外源磁场的分离,早期主要根据电离层-磁层磁场的时空特征,借助于各类筛选条件(如地方时条件、纬度条件、Dst/ Kp指数条件)对外源场扰动变化较大的数据进行剔除,然后采用多项式拟合[3]、球谐分析[4]或沿轨滤波技术对外源场的贡献进行分离[5],如基于沿轨快速傅里叶变换的Kaiser滤波[6]或Butterworth滤波[7],但滤波难点在于截止波长的确定[5, 7],容易引入信号泄露误差,应用时需分析未模型化的残差信号强度。尽管早期国内外学者基于卫星磁测数据,依据外源场的时空特征对内外场源分离和模型参数化进行了相关理论分析[8-12],但由于单颗卫星完成覆盖全球的观测周期通常为几天,难以实现全球范围内电离层-磁层磁场的准同步观测[13],因此,很难建立实时高精度的外源场模型。2013-11-22,Swarm卫星成功发射,不仅为地磁场提供了迄今为止的最高精度测量,而且首次实现了对地磁场进行1 h到数天时间尺度上变化磁场的全球观测,Swarm星群的高时空覆盖采样、优化的数据处理流程[5, 14],有利于分析空间电磁环境的时空变化规律,为外源场的完整分离和模型化提供了强有力的数据和技术支撑。2016年,欧空局Swarm星群应用与研究中心分别发布了最新的电离层磁场MIO_SHA_2D[15]和磁层磁场MMA_SHA_2C[10]系列模型,为外源场的改正提供了重要的参考和依据。

    内源场和外源场的信号频率存在一定重叠[1, 16],很难通过滤波方法对其进行准确分离。若对全球或局部岩石圈磁场进行建模,在实施主磁场模型改正之后,剩余磁场主要包括岩石圈磁场、外源场及其感应磁场,主磁场残留的长波信号比较微弱。相对于外源场的全球尺度,起源于磁性岩石圈的岩石圈磁场,在卫星高度上的空间尺度通常为几百公里,甚至更小,因此在建立全球或局部岩石圈磁场模型过程中,可以尝试采用沿轨高通滤波对外源场进行分离。

    高通滤波器从实现方法上可分为无限冲激响应(infinite impulse response,IIR)和有限冲激响应(finite impulse response,FIR)[17-18]。由于外源场时变特征明显,为使不同频率分量的信号经滤波后具有精确统一的时标,本文采用具有精确线性相位特征的FIR滤波器。FIR滤波器的常用设计方法主要包括窗函数法和频率采样法[13]。窗函数法相对简单,物理意义清晰,通过选择合适的窗函数可以有效抑制信号泄露。

    为了分离剩余磁场中外源场的长波信号,本文研究基于FIR滤波器的Hamming窗函数法,通过设计合理的滤波参数,对Swarm卫星磁静日中低磁纬度(±50°)的观测数据进行沿轨高通滤波,并引入最新的MIO_SHA_2D、MMA_SHA_2C电离层-磁层磁场改正模型,对滤波方法的可行性和精度进行评估。

    • 本文的主要目的是采用沿轨高通滤波对外源场信号进行改正,并通过最新的国际高精度外源场改正模型验证该方法的可行性。为此,本文选取了2015年1月Swarm A星磁静日中低磁纬度(±50°)的低频观测数据(采样频率为1 Hz),在进行数据筛选之前,首先进行半圆轨道分离和10 s间隔的二次重采样,轨道分离的目的是为了对观测数据进行沿轨滤波。在卫星轨道高度,10 s间隔的二次重采样大约对应于75 km的空间分辨率,可减弱小尺度磁异常信号的混频效应。

      考虑到MIO_SHA_2D电离层磁场模型主要描述了中低纬度地区磁静期Sq和赤道电集流磁场及其感应磁场[15],本文参考国际上大多数地磁场建模中所用到的地磁指数条件[16, 19-20],利用式(1)选取地磁扰动相对平静的观测时段数据:

      $$ \left\{ {\begin{array}{*{20}{l}} {|{\rm{Dst}}| \le 20\;\;{\rm{nT}}\;\;{\rm{ and }}\;\;{\rm{d}}|{\rm{Dst}}| \le 10\;\;{\rm{nT}}}\\ {Kp \le 2\;\;{\rm{ and }}\;\;{\rm{d}}|Kp| \le 2} \end{array}} \right. $$ (1)

      式中,d| Dst|和d| Kp |分别表示Dst、Kp在前3 h内的变化值。为了充分保证数据的观测质量,通过flags值剔除了Star-camera运行状态欠佳的观测数据,并对XYZ矢量和标量F互差大于10 nT的观测数据予以剔除,最后,剔除了连续超过240个历元没有数据的短弧轨道。图 1给出了基于卫星观测值提取磁静期剩余磁场的预处理流程。

    • 由于MIO_SHA_2D和MMA_SHA_2C是基于地磁坐标系所建立的电离层磁场模型,因此在进行模型计算之前,需先将各测点的位置坐标从球坐标系(λ, θ)转换到地磁坐标系(Λ, Θ),然后再将模型所计算的地磁分量从地磁坐标系(N, E, U)dp转换到球坐标系(N, E, U)geo。此外,MIO_SHA_2D模型还涉及到协调世界时(coordinated universal time, UTC)与地磁世界时tm之间的转换。

    • 1) 球坐标系与地磁坐标系之间转换,可通过球面三角公式实现[1, 3]

      $$ \left\{\begin{array}{l}{\cos \mathit{\Theta}=\cos \theta \cos \theta_{0}+\sin \theta \sin \theta_{0} \cos \left(\lambda-\lambda_{0}\right)} \\ {\sin \mathit{\mathit{\Lambda}}=\sin \theta \sin \left(\lambda-\lambda_{0}\right) / \sin \mathit{\Theta}}\end{array}\right. $$ (2)

      式中,(Λ, Θ)表示点P的地磁经度和地磁余纬;(λ, θ)表示点P的地心经度和地心余纬;(λ0, θ0)表示地磁北极Nm的地心经度和地心余纬。

      2) J2000惯性坐标系与地固空间直角坐标系之间的转换。

      计算地磁世界时tm需要确定太阳在地磁坐标系中的日下点经度。对于给定的UTC时间,通常先根据JPL提供的星历文件计算太阳在J2000惯性坐标系下的位置和速度,然后按照式(3)将其转换到地心地固直角坐标系下。

      $$ {(x, y, z)^{\rm{T}}} = (\mathit{\boldsymbol{PR}})(\mathit{\boldsymbol{NR}})(\mathit{\boldsymbol{B}}1)(\mathit{\boldsymbol{B}}2){(X, Y, Z)^{\rm{T}}} $$ (3)

      式中,(PR)是岁差矩阵;(NR)是章动矩阵;(B1)是极移矩阵;(B2)是地球自转矩阵[20-22],其定义可参考文献[23]。

      3) 地磁观测分量从地磁坐标系转换到球坐标系下,可通过如下步骤实现:

      ① 在地磁坐标系参考框架下,利用式(4)将站心地磁计算分量(N, E, U)dp转换到地磁空间直角坐标系(X, Y, Z)dp

      $$ (X, Y, Z)_{\mathrm{dp}}=\boldsymbol{R}_{1} \cdot(N, E, U)_{\mathrm{dp}} $$ (4)
      $$ \boldsymbol{R}_{1}=\left[\begin{array}{ccc}{-\sin \mathit{\Phi} \cos \mathit{\Lambda}} & {-\sin \mathit{\Lambda}} & {-\cos \mathit{\Phi} \cos \mathit{\Lambda}} \\ {-\sin \mathit{\Phi} \sin \mathit{\Lambda}} & {\cos \mathit{\Lambda}} & {-\cos \mathit{\Phi} \sin \mathit{\Lambda}} \\ {\cos \mathit{\Phi}} & {0} & {-\sin \mathit{\Phi}}\end{array}\right] $$ (5)

      式中,R1为旋转矩阵;Φ为点P的地磁纬度。

      ② 利用式(6)将(X, Y, Z)dp转换到地心球空间直角坐标系(X, Y, Z)geo

      $$ (X, Y, Z)_{\mathrm{geo}}=\boldsymbol{R}_{2}^{-1} \cdot(X, Y, Z)_{\mathrm{dp}} $$ (6)
      $$ \boldsymbol{R}_{2}=\left[\begin{array}{ccc}{\cos \theta_{0}} & {0} & {-\sin \theta_{0}} \\ {0} & {1} & {0} \\ {\sin \theta_{0}} & {0} & {\cos \theta_{0}}\end{array}\right]\left[\begin{array}{ccc}{\cos \lambda_{0}} & {\sin \lambda_{0}} & {0} \\ {-\sin \lambda_{0}} & {\cos \lambda_{0}} & {0} \\ {0} & {0} & {1}\end{array}\right] $$ (7)

      式中,R2为旋转矩阵。

      ③ 在地心球坐标参考框架下,利用式(8)将(X, Y, Z)geo转换为站心地磁观测分量(N, E, U)geo

      $$ (N, E, U)_{\mathrm{geo}}=\boldsymbol{R}_{3} \cdot(X, Y, Z)_{\mathrm{geo}} $$ (8)
      $$ \boldsymbol{R}_{3}=\left[\begin{array}{ccc}{-\sin \varphi \cos \lambda} & {-\sin \varphi \sin \lambda} & {\cos \varphi} \\ {-\sin \lambda} & {\cos \lambda} & {0} \\ {-\cos \varphi \cos \lambda} & {-\cos \varphi \sin \lambda} & {-\sin \varphi}\end{array}\right] $$ (9)

      式中,R3为旋转矩阵;φ为点P的地心纬度。

    • 1) 地磁世界时tm的计算[15]

      $$ t_{m}=\left(180^{\circ}-\varphi_{d, s}\right) / 15 $$ (10)

      式中,φd, s为地磁坐标系日下点的地磁经度。

      2) 格里高利历与儒略日JD之间的转换:

      $$ \begin{array}{l} {\rm{JD}} = {\mathop{\rm int}} (365.25y) + {\mathop{\rm int}} (30.6001(m + 1)) + \\ D + h{\rm{ / 24 + 1720981}}{\rm{.5}} \end{array} $$ (11)

      式中,当M≤2时,则y=Y-1,m=M+12;当M>2时,则y=Ym=Mh=H+Min/60+S/3 600,其中Y为年,M为月,D为日,H、Min、S分别为时、分、秒;JD为儒略日;int表示向下取整。

    • 相对于岩石圈磁场,外源场的空间分布具有全球性,其信号波长较长,频率较低,因此可借助高通滤波器对外源场信号进行分离。高通滤波器从实现方法上可分为无限冲激响应IIR和有限冲激响应FIR。考虑到外源场时变特征明显,不同频率分量的信号经滤波后需要精确统一的时标,本文选用具有精确线性相位特性的FIR滤波器。FIR滤波器的常用设计方法主要包括窗函数法和频率采样法[24]。窗函数法是从时域进行设计,频率采样法是从频域进行设计。窗函数法相对简单、物理意义清晰,通过选择合适的窗函数可以有效抑制信号泄露,本节将重点研究利用Hamming窗函数法分离外源场信号的FIR高通滤波器参数设置问题。

    • Hamming窗函数的基本参数包括通带、阻带截止频率、过渡带宽及滤波器长度。卫星磁测数据中,信号频率f与卫星速度v、波长分辨率λ之间的关系为f=v/λ[25]。当卫星速度一定时,通带截止频率fp决定了输出信号的波长分辨率。若已知岩石圈磁场和外源场的临界频率,则可以明确确定滤波器的通带截止频率,事实上岩石圈磁场和外源场信号频率存在一定的重叠,因此需要设计合理的滤波参数,进而有效抑制吉布斯效应。

      阻带截止频率fs取决于过渡带宽Δω,Δω=ωpωs,对于Hamming窗,过渡带宽与滤波器的窗口长度N之间的关系为:

      $$ \Delta \omega=\omega_{p}-\omega_{s} \geqslant 6.6 \pi / N $$ (12)

      式中,ωpωs分别为滤波器的通带和阻带数字截止频率,

      $$ \left\{\begin{array}{l}{\omega_{p}=f_{p} \times 2 \pi / f_{s t}} \\ {\omega_{s}=f_{s} \times 2 \pi / f_{s t}}\end{array}\right. $$ (13)

      相应的归一化的模拟截止频率分别为:

      $$ \left\{\begin{array}{l}{f_{p}^{{N}}=f_{p} \times 2 / f_{s t}} \\ {f_{s}^{{N}}=f_{s} \times 2 / f_{s t}}\end{array}\right. $$ (14)
    • 在给定通带、阻带截止频率的情况下,可根据式(12)近似计算滤波器窗口长度N。滤波器窗口长度并不具有唯一性,在确定窗口长度时,应着重考虑两个原则:①滤波器长度N的奇偶性,对于长度相当的滤波器,N为偶数时性能比N为奇数时好,但显著的缺点是滤波时产生的时延(N-1)/2不具有整数特性,需进行内插[24];②窗口长度越短,滤波器群时延越小;窗口长度越长,过渡带减小,幅频响应更加精准,但当信号频率与采样频率之比很小时,可能引起类似码间串扰效应,无法得到稳定输出。

    • 本文分析了磁静日期间2015年1月Q1~Q10(平静日)10 d的Swarm A星观测数据,采用图 1的数据筛选方法选取Q1~Q10观测时段相对平静的半圆轨道,然后采用EMM2015模型(1~720阶)对内源场进行改正,提取剩余磁场。

      图  1  基于卫星磁测数据提取磁静期剩余磁场的预处理流程

      Figure 1.  Preprocessing Process of Extracting Magnetic Remaining Magnetic Field Based on Satellite Magnetic Data

      以2015-01-20某条半圆轨道为例,图 2图 3分别给出了地磁纬度为±50°范围内沿轨方向上XYZ分量的剩余磁场强度及径向剩余磁场Z分量的幅频响应特性。从图 2可以看出,在中低纬度地区,剩余磁场强度为±10~±20 nT;随着纬度的升高,XY分量的剩余磁场表现较为平稳,而Z分量则表现为增强趋势,但其强度无明显突变,说明在中低纬度地区的磁静期,相对于水平方向,外源场对径向分量Z的影响相对较大。从图 3可以看出,剩余磁场主要集中在低频段,表现为长波信号,这些长波信号主要源于大尺度磁层磁场,其次源于夜间电离层磁场及EMM2015模型的残差项。

      图  2  剩余磁场XYZ分量强度

      Figure 2.  Remaining Magnetic Field X, Y, Z Components Intensity

      图  3  剩余磁场Z分量的幅频响应特性

      Figure 3.  The Amplitude-Frequency Response of the Z Component of the Residual Magnetic Field

      下面分别采用Hamming窗函数法和MMA_SHA_2C、MIO_SHA_2D模型对剩余磁场的外源场成分进行改正和精度分析。根据式(12),设定滤波窗口长度为87阶,通带和阻带波长分别为1 200 km和3 000 km。图 4给出了两种方法计算的径向分量(Z)强度,图 5图 6分别给出了两种方法改正后的径向残差(ΔZ)大小及相应的幅频响应特性,其中黑色曲线代表模型改正结果,红色曲线代表滤波结果。

      图  4  径向分量Z的外源场改正强度

      Figure 4.  The Corrected Intensity of External Field of Radial Component Z

      图  5  改正后的径向残差ΔZ大小

      Figure 5.  Corrected Radial Residual ΔZ

      图 4可以看出,滤波结果与模型改正结果的一致性吻合较好。图 5表明:①采用两种方法改正后的径向残差均优于±3 nT,对比图 2可以看出,两种方法对大尺度的剩余磁场均能作出有效改正;②两种方法改正后的残差互差约为2 nT,相对于模型改正,滤波方法对大尺度的磁场信号改正更为彻底,但滤去的信号并不能完全界定为外源场,因为这些大尺度的磁场含有未模型化的主磁场信号,甚至可能含有大尺度的磁异常信号;③在沿轨方向上,经滤波和模型改正后的残差分布具有相似的变化特征,表明两种方法对大尺度剩余磁场的改正具有一致性,图 5所描述的径向残差Z分量具有相似的幅频响应特性,验证了该滤波方法的可行性。由图 6可知,滤波后的低频信号强度大大减弱,这表明基于Hamming窗函数法能有效滤去大尺度的剩余磁场信号。

      图  6  径向残差ΔZ的幅频响应特性

      Figure 6.  The Amplitude-Frequency Response of Radial Residual ΔZ

      为了进一步分析滤波器长度对剩余磁场长波信号改正的精度影响,在给定通带波长1 200 km、阻带波长3 000 km的情况下,以外源场模型改正后的残差信号强度为参考,采用Hamming窗函数法,对2015年1月Q1~Q10共10 d经筛选后的14条半圆轨道的剩余磁场进行沿轨高通滤波,通过引入均方根误差(root mean square error,RMSE),对比分析不同窗口长度滤波后的残差与模型改正后的残差,表 1给出了统计分析的结果。

      表 1  基于均方根误差(RMSE)的滤波器长度对剩余磁场长波信号改正的精度影响/nT

      Table 1.  Effect of Filter Length Based on Root Mean Square Error (RMSE) on the Correction Precision of Long Wave Signal of Remaining Magnetic Field/nT

      N RMSE最小值 RMSE最大值 RMSE平均值
      X Y Z X Y Z X Y Z
      55 0.461 0.562 0.430 4.555 3.550 4.001 1.808 1.715 1.537
      87 0.438 0.504 0.498 4.498 3.944 3.543 1.825 1.767 1.397
      97 0.442 0.501 0.500 4.500 3.920 3.472 1.825 1.761 1.391
      129 0.452 0.562 0.497 4.503 3.899 3.214 1.827 1.780 1.347
      161 0.461 0.675 0.498 4.500 4.004 3.108 1.831 1.859 1.315
      193 0.472 0.779 0.511 4.512 4.246 2.738 1.846 1.939 1.246
      225 0.497 0.462 0.525 4.594 4.463 2.033 1.862 1.855 1.192

      表 1可以看出,随着滤波器长度的增加,在水平方向上,RMSE呈增大趋势,在径向上,RMSE的最小值在增大,最大值在减小,其均值呈现减小趋势,表明剩余磁场水平分量和径向分量对滤波器长度具有不同的响应特性。因此,为了均衡滤波器在水平方向和径向上的滤波精度,滤波器的长度不宜选得过小或过大。

    • 本文利用Swarm A星2015年1月的磁静日观测数据,采用Hamming窗函数法沿轨高通滤波对大尺度的外源场信号进行分离,并借助于最新的MIO_SHA_2D和MMA_SHA_2C电离层-磁层磁场模型验证了该方法的可行性,同时给出了滤波器窗口长度的确定方法,为滤波器的参数设计提供参考。研究表明:

      1) 经EMM2015(1~720阶)模型改正后的剩余磁场主要以外源场信号为主,其次源于EMM2015模型的残差项,利用Hamming窗函数法的沿轨高通滤波可对中低纬度地区大尺度的外源场信号实施有效分离,其结果与外源场模型的计算值一致性吻合较好;

      2) 采用Hamming窗函数法和外源场模型改正后的径向残差均优于3 nT,相对于模型改正,滤波方法对大尺度的磁场信号分离更为彻底,若以外源场模型为参考,则滤波后的残差约为2 nT,表明该方法对外源场的分离具有较高的精度,但滤波存在信号泄漏问题,因此,该方法只能为外源场信号的分离提供借鉴和参考,其可靠性尚需作进一步的论证;

      3) 剩余磁场水平分量和径向分量对滤波器的窗口长度具有不同的响应特性,在给定通带波长1 200 km、阻带波长3 000 km的情况下,将滤波器长度设置为87~129阶,可均衡滤波器在水平方向和径向上的滤波精度。

参考文献 (24)

目录

    /

    返回文章
    返回