[发明专利]地磁Kp指数的反演方法、系统、电子装置及存储介质有效
申请号: | 201911051348.2 | 申请日: | 2019-10-31 |
公开(公告)号: | CN110850492B | 公开(公告)日: | 2021-09-07 |
发明(设计)人: | 薛炳森;周率;贺晓冬;李绍帅 | 申请(专利权)人: | 国家卫星气象中心(国家空间天气监测预警中心);易天气(北京)科技有限公司 |
主分类号: | G01V3/40 | 分类号: | G01V3/40 |
代理公司: | 北京汉之知识产权代理事务所(普通合伙) 11479 | 代理人: | 陈敏 |
地址: | 100081 北京市海淀*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 地磁 kp 指数 反演 方法 系统 电子 装置 存储 介质 | ||
1.一种地磁Kp指数的反演方法,其特征在于,包括以下步骤:
在地磁台记录观测日的地磁场观测值,并选取其中的水平分量观测值,形成所述水平分量观测值的时间序列Fi;
在所述水平分量观测值中去除地磁场的背景磁场Fbi,得到所述观测日的不同时刻的地磁场扰动幅度值DFi;
根据所述地磁扰动幅度值DFi计算所述观测日的每个时段的累积地磁扰动幅度值DTj;
根据所述累积地磁扰动幅度值由以下公式反演所述地磁Kp指数:
Kp=Aj[lg(DTj)+Bj],
其中,i表示选定时间点距所述观测日的00:00时刻的分钟数;一天中不同的时段对应不同的Aj和Bj,j为时间段序号。
2.根据权利要求1所述的地磁Kp指数的反演方法,其特征在于,在所述水平分量观测值中去除地磁场的背景磁场Fbi,得到所述观测日的不同时刻的地磁场扰动幅度值DFi之前还包括:
在所述观测日之前的预定天数内选择磁静日;
提取所述磁静日的地磁场水平分量观测值,将所述磁静日的地磁场水平分量观测值作为所述背景磁场。
3.根据权利要求2所述的地磁Kp指数的反演方法,其特征在于,选择自所述观测日起之前的27天为所述预定天数。
4.根据权利要求3所述的地磁Kp指数的反演方法,其特征在于,在所述观测日之前的预定天数内选择磁静日还包括以下步骤:
读取所述27天内每一天的地磁场扰动幅度指数Ap;
选择所述地磁场扰动幅度指数Ap最小的自然日作为所述磁静日。
5.根据权利要求4所述的地磁Kp指数的反演方法,其特征在于,在所述观测日之前的预定天数内选择磁静日还包括以下步骤:
如果在所述27天内出现地磁场扰动幅度指数Ap相同的多个自然日,则选择距离所述观测日最近的自然日作为所述磁静日。
6.根据权利要求1所述的地磁Kp指数的反演方法,其特征在于,在所述观测日的00:00时刻开始划分所述时段,每3个小时划分为一个时段。
7.根据权利要求6所述的地磁Kp指数的反演方法,其特征在于,由以下公式计算所述观测日的每个时段的累积地磁扰动幅度值DTj:
8.一种地磁Kp指数的反演系统,其特征在于,包括:
地磁场观测单元,用于在地磁台记录观测日的地磁场观测值,并选取其中的水平分量观测值,形成所述水平分量观测值的时间序列Fi;
地磁场扰动幅度值计算单元,与所述地磁观测单元连接,用于接收所述水平分量观测值的时间序列Fi,并在所述水平分量观测值中去除地磁场的背景磁场Fbi,得到所述观测日的不同时刻的地磁场扰动幅度值DFi;
累积地磁扰动幅度值计算单元,与所述地磁场扰动幅度值计算单元连接,用于接收所述地磁场扰动幅度值DFi,并根据所述地磁场扰动幅度值DFi计算累积地磁扰动幅度值DTj;
地磁Kp指数反演单元,与所述累积地磁扰动幅度值计算单元连接,用于接收所述累积地磁扰动幅度值DTj,并根据所述累积地磁扰动幅度值DTj根据如下公式反演所述地磁Kp指数:
Kp=Aj[lg(DTj)+Bj],
其中,i表示选定时间点距所述观测日的00:00时刻的分钟数;一天中不同的时段对应不同的Aj和Bj,j为时间段序号。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于国家卫星气象中心(国家空间天气监测预警中心);易天气(北京)科技有限公司,未经国家卫星气象中心(国家空间天气监测预警中心);易天气(北京)科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201911051348.2/1.html,转载请声明来源钻瓜专利网。