[發明專利]一種基于地磁匹配的地磁日變修正方法

有效

• 申請號: 201310028198.X
• 申請日: 2013-01-24
• 公開（公告）號: CN103115624A
• 公開（公告）日: 2013-05-22
• 發明（設計）人: 黃朝艷;鄧翠婷;趙華;田海冬;張華龍
• 申請（專利權）人: 南京航空航天大學
• 主分類號: G01C21/08
• 分類號: G01C21/08
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 朱小兵

1.一種基于地磁匹配的地磁日變修正方法，其特征在于，包括以下步驟：

步驟1，選取地磁匹配特征量，所述地磁匹配特征量為三個地磁分量，包括北向分量X、東向分量Y、垂直分量Z；

步驟2，基于FMI方法從地磁匹配特征量的實測數據中提取各地磁匹配特征量的日變量；

步驟3，通過地磁匹配特征量的實測數據減去與其對應的日變量，建立地磁實時圖；

步驟4，基于多維特征量匹配對地磁實時圖和地磁基準圖進行匹配；

步驟5，輸出匹配結果。

2.根據權利要求1所述的基于地磁匹配的地磁日變修正方法，其特征在于，所述步驟2涉及的實測數據是由磁傳感器采集的；所述基于FMI方法從實測數據中提取日變量包括以下步驟：

步驟2.1，根據地磁匹配特征量的擾動幅度a分別確定與各地磁匹配特征量相對應的K指數的初始值K_X0、K_Y0、K_Z0，取K₀=max(K_X0,K_Y0,K_Z0)，所述擾動幅度a為每3小時時段內地磁匹配特征量的實測數據的最大值和最小值間的差值；

步驟2.2，地磁匹配特征量的時均值由該小時及其前后各(m+n)分鐘的數據平均值來計算，然后利用5次諧波函數對時均值進行擬合，生成S_R曲線；所述n=K³³；m與地方時有關：地方時3點到6點和18點到21點，m取60；地方時7點到17點，m取0；其余時間，m取90；

步驟2.3，分別計算地磁匹配特征量的實測數據與S_R曲線的殘差序列，由殘差序列的3小時擾動幅度a分別確定與各地磁匹配特征量相對應的K指數的中間值K_Xm、K_Ym、K_Zm，取K_m=max(K_Xm,K_Ym,K_Zm)；

步驟2.4，重復步驟2.2和步驟2.3的計算一次，再次得到K指數的中間值

K_m=max(K_Xm,K_Ym,K_Zm)，此K_m即為K指數的最終值；

步驟2.5，由K指數的最終值，代入步驟2.2計算，生成最終S_R曲線；

步驟2.6，計算最終S_R曲線與地球主磁場和磁異常間的差值，得到地磁匹配特征量的日變量。

3.根據權利要求1所述的基于地磁匹配的地磁日變修正方法，其特征在于，所述步驟4涉及的地磁基準圖是指基于地磁模型建立的基準數據庫，或基于測量資料繪制的地磁圖。