张秀霞 杨冬梅 杨福喜 何宇飞 张素琴

1） 江苏省地震局，南京 210014

2 ）中国地震局地球物理研究所，北京 100081

3） 新疆地震局，乌鲁木齐 830011

地磁观测中及时了解地磁场日变化的扰动情况，可以帮助人们及时判断数据的可靠性，识别观测中的异常。研究地磁静日变化的预测方法在地磁观测数据分析处理中具有应用意义。

地磁日变化主要包含了地球外部空间电离层和磁层电流体系的变化信息，有静日变化和扰日变化之分。通常情况下人们利用每月中的5 天静日变化的平均作为当月的Sq 变化的估计。由于每月的静扰日日期通常在下一个月的中旬产出，时效性较差，不能满足日常地磁监测工作的需要。

何宇飞等曾经基于四分位法（IQR），利用研究日之前15 天内的日变化数据预测当日的日变化数据，取得了较好的效果。但其研究的时间段是11 年太阳活动周期中较平静的时段，自2012 年起，太阳活动趋于剧烈，15 天之内连续遇到磁场扰动的情况不可避免，因此需要研究一种预测地磁静日变化的新方法。

地磁静日确定的一个重要依据是Kp 指数的大小。研究表明，Vr 指数与Kp 指数有非常好的线性相关，而Vr 指数计算简单方便，因此我们考虑是否可以基于研究日之前某一时段的Vr 指数挑选日变化相对平静的日期，用于预测研究日的地磁日变化。

研究使用了新疆局乌鲁木齐地磁台2008 年的观测数据，计算了其全年的Vr 指数及其日总和，与Kp指数的日总和、每月的静日日期等进行了综合对比，发现每月最平静的5 个日期均在当月Kp 日总和最小的10 日内，但在当月Vr 日总和最小的10 日内的平均概率是82%。因此，我们选用研究日之前30 天内的数据，根据此30 天的Vr 日总和挑选出10 天相对平静的日期及其日变化数据，按每分钟将10 个数据点进行分析，去掉离均差最大的4 个点，取剩余的6 个点的平均值作为该分钟地磁场日变化的预测值。我们将此预测值与IQR 方法预测结果及传统的5 个静日平均结果进行了系统的比较，探讨了这种地磁静日变化预测方法的可行性和应用前景。