[发明专利]一种K指数现报方法有效
申请号: | 201810982711.1 | 申请日: | 2018-08-27 |
公开(公告)号: | CN108983301B | 公开(公告)日: | 2019-11-29 |
发明(设计)人: | 付梦印;康嘉鹏;刘彤;王美玲;王凯;陈钊苇 | 申请(专利权)人: | 北京理工大学 |
主分类号: | G01V3/40 | 分类号: | G01V3/40;G06F17/50 |
代理公司: | 11120 北京理工大学专利中心 | 代理人: | 郭德忠;李爱英<国际申请>=<国际公布> |
地址: | 100081 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 地磁场 预处理 变化曲线 统计学 频谱 | ||
1.一种K指数现报方法,其特征在于,包括以下步骤:
S1:获取求K指数日的前N日的规则日变化曲线SR,具体为:
将所述前N日中每一日平均划分为8个时段;
分别获取每一日各时段的地磁场H分量的每3小时平均值,其中,所述每3小时平均值的计算方法为:
将每3小时以及每3小时的前后一段时间的地磁场H分量的平均值,作为该时段对应的每3小时平均值;
分别对各日的8个时段对应的每3小时平均值进行谐波拟合,获取各日的规则日变化曲线SR;
S2:将步骤S1获取的N条规则日变化曲线SR的平均值,作为求K指数日的规则日变化曲线SR*;
S3:将求K指数日的地磁场H分量观测曲线与规则日变化曲线SR*的差值曲线,作为求K指数日的K变化曲线;
S4:根据所述K变化曲线,获取求K指数日的K指数。
2.如权利要求1所述的一种K指数现报方法,其特征在于,所述获取求K指数日的前N日的规则日变化曲线SR前,确定求取规则日变化曲线SR所需的天数N,具体为:
按照步骤S1的方法,分别获取不同天数对应的规则日变化曲线SR#;
分别获取各规则日变化曲线SR#对应的备选求K指数日的磁扰幅度;
分别获取备选求K指数日的各个磁扰幅度对数与K指数之间相关系数,其中,所述K指数为备选求K指数日的K指数;
剔除不满足预设条件的相关系数,其中,所述预设条件为:相关系数大于设定经验值,且此时大于设定经验值的相关系数对应的天数小于设定经验天数;
将满足预设条件的相关系数的最大值对应的天数,作为求取规则日变化曲线SR所需的天数N。
3.如权利要求1所述的一种K指数现报方法,其特征在于,当计算凌晨3~6点和黄昏18~21点的每3小时平均值时,将每3小时以及每3小时的前60分钟与后n分钟的地磁场H分量的平均值,作为该时段对应的每3小时平均值;
当计算夜间21~3点的每3小时平均值时,将每3小时以及每3小时的前90分钟与后n分钟的地磁场H分量的平均值,作为该时段对应的每3小时平均值;
当计算其余时间的每3小时平均值时,将每3小时以及每3小时的后n分钟的地磁场H分量的平均值,作为该时段对应的每3小时平均值;
其中,n=K3.3,K根据每3小时时段中,地磁场H分量与地磁场H分量的平均值的差值最大值与差值最小值确定,具体的,所述差值最大值与差值最小值之间的差值,与所述K的取值满足预设的一一对应关系。
4.如权利要求1所述的一种K指数现报方法,其特征在于,步骤S4中所述根据所述K变化曲线,获取求K指数日的K指数具体为:
将所述K变化曲线平均划分为8个时段,获取每个时段的K变化的最大值和最小值,并根据各时段最大值和最小值的差值与预设的K指数/H变幅对应表的对应关系,获取每个时段的K指数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京理工大学,未经北京理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810982711.1/1.html,转载请声明来源钻瓜专利网。