|
|
The calculation method of full tensor geomagnetic gradient based on IGRF model |
Yang ZHONG, Yan-Wu GUAN, Jia-Qiang SHI, Feng XIAO |
College of Geo-Exploration Science and Technology,Jilin University,Changchun 130026,China |
|
|
Abstract The international geomagnetic reference field (IGRF) is a general international model for describing the earth’s main magnetic field. At present,this model can be used to calculate the seven elements of geomagnetic field at any point. However,with the development of aeronautical full tensor magnetic measurement technology,there is an urgent need for full tensor geomagnetic gradient data. In this paper,the calculation principle of the IGRF model is summarized and the expression of the full tensor geomagnetic gradient with spherical harmonic expansion is derived. The calculation of the seven elements of geomagnetic field and the full tensor geomagnetic gradient at any given point is realized. Comparing with the calculated data from the website of the National Oceanic and Atmospheric Administration of the United States (NOAA),the results are accurate and reliable. The contour map of the full tensor geomagnetic field in a region is drawn, and the results were in accordance with the Laplace equation. It provides the theoretical basis for the selection of learning flight working area and flight height in the aeromagnetic survey.
|
Received: 02 September 2019
Published: 24 June 2020
|
|
|
|
|
|
Geographic coordinate system and geocentric coordinate system
|
n | Pn(cosθ) | Pn(v) | 0 | 1 | 1 | 1 | cosθ | μ | 2 | (3cos2θ+1)/4 | (3v2-1)/2 | 3 | (5cos3θ+3cosθ)/8 | (5v3-3v)/2 |
|
Legendre polynomials of order 0~3
|
n | m | 伴随勒让德多项式 Pn,m(cosθ) | 高斯规格化伴随勒让德多项式 Pn,m(cosθ) | 施密特拟规格化伴随勒让德多项式 (cosθ) | 0 | 0 | 1 | 1 | 1 | 1 | 0 | cosθ | cosθ | cosθ | 1 | 1 | sinθ | sinθ | sinθ | 2 | 0 | (3cos2θ-1) | (3cos2θ-1) | (3cos2θ-1) | 2 | 1 | 3cosθsinθ | cosθsinθ | cosθsinθ | 2 | 2 | 3sin2θ | sin2θ | sin2θ | 3 | 0 | cosθ(5cos2θ-3) | cosθ(5cos2θ-3) | cosθ(5cos2θ-3) | 3 | 1 | sinθ(5cos2θ-1) | sinθ(5cos2θ-1) | sinθ(5cos2θ-1) | 3 | 2 | 15cosθsin2θ | cosθsin2θ | cosθsin2θ | 3 | 3 | 15sin3θ | sin3θ | sin3θ |
|
Order 0-3 adjoint Legendre polynomials
|
|
Contour map of the total geomagnetic field
|
|
Contour map of geomagnetic field deviation(a) and dip angle(b)
|
|
Three-component contour map of the geomagnetic field
|
|
Full-tensor geomagnetic gradient contour map
|
|
Contour map of the sum of the three tensor components on the main diagonal
|
城市 | 来源 | Bx/nT | By/nT | Bz/nT | Bh/nT | B/nT | D/° | I/° | 成都市 | 本文计算 | 33866.2 | -1213.9 | 37855.4 | 33887.9 | 50807.6 | -2.0528 | 48.1652 | NOAA | 33866.2 | -1213.9 | 37855.4 | 33887.9 | 50807.7 | -2.0529 | 48.1652 | 自贡市 | 本文计算 | 34644.5 | -1268.7 | 36051.7 | 34667.6 | 50015.7 | -2.0972 | 46.1212 | NOAA | 34644.5 | -1268.7 | 36051.7 | 34667.7 | 50015.8 | -2.0973 | 46.1212 | 泸州市 | 本文计算 | 34909.9 | -1334.7 | 35359.9 | 34935.3 | 49707.1 | -2.1895 | 45.346 | NOAA | 34909.9 | -1334.7 | 35359.9 | 34935.4 | 49707.2 | -2.1895 | 45.346 | 德阳市 | 本文计算 | 33580.9 | -1266.4 | 38443.8 | 33604.7 | 51060.8 | -2.1597 | 48.8424 | NOAA | 33580.9 | -1266.4 | 38443.8 | 33604.8 | 51060.8 | -2.1597 | 48.8424 | 最大绝对误差 | 0 | 0 | 0 | 0.1 | 0.1 | 0.0001 | 0 |
|
Table of point position magnetic field values in the test area
|
[1] |
Barraclough D R, Harwood J M, Leaton B R, et al. A definitive model of the geomagnetic field and its secular variation for 1965—I.Derivation of model and comparison with the IGRF[J]. Geophysical Journal International, 1978,55(1):111-121.
|
[2] |
安振昌, 徐元芳. 球谐分析与国际地磁参考场[J]. 地震地磁观测与研究, 1982,3(5):22-26.
|
[2] |
An Z C, Xu Y F. Spherical harmonic analysis and international geomagnetic reference[J]. Seismological and Geomagnetic Observation and Research, 1982,3(5):22-26.
|
[3] |
安振昌, 徐元芳. 利用中国地区的资料检验全球地磁场的球谐模式[J]. 地震地磁观测与研究, 1983,4(5):1-13.
|
[3] |
An Z C, Xu Y F. The spherical harmonic model of the global geomagnetic field is tested with data from China[J]. Seismological and Geomagnetic Observation and Research, 1983,4(5):1-13.
|
[4] |
徐元芳, 王月华, 安振昌. IGRF1980的检验[J]. 地震地磁观测与研究, 1986,7(2):10-17.
|
[4] |
Xu Y F, Wang Y H, An Z C. IGRF1980 inspection[J]. Seismological and Geomagnetic Observation and Research, 1986,7(2):10-17.
|
[5] |
安振昌. 第四代国际地磁参考场[J]. 地震地磁观测与研究, 1986,7(3):20-25.
|
[5] |
An Z C. Fourth generation international magnetic reference field[J]. Seismological and Geomagnetic Observation and Research, 1986,7(3):20-25.
|
[6] |
安振昌, 王月华, 徐元芳, 等. IGRF的计算与评价[J]. 物探化探计算技术, 1988,10(2):93-99.
|
[6] |
An Z C, Wang Y H, Xu Y F, et al. Calculation and evaluation of IGRF[J]. Computing Techniques for Giophysical and Geochenical Exploration, 1988,10(2):93-99.
|
[7] |
安振昌. 利用 IGRF 计算地磁台年均值的方法[J]. 地震地磁观测与研究, 1987,8(5):15-18.
|
[7] |
An Z C. Method of calculating the annual mean values of observatory by using IGRF[J]. Seismological and Geomagnetic Observation and Research, 1987,8(5) : 15-18.
|
[8] |
安振昌. 利用第七代IGRF计算1995.0年中国及邻区地磁场[J]. 地球物理学进展, 1997,12(3):22-26.
|
[8] |
An Z C. Geomagnetic field over China and adjacent at epoch 1995.0derived from the seventh generation IGRF[J]. Progress in Geophysics, 1997,12(3):22-26.
|
[9] |
安振昌, 王月华. 1900~2000年非偶极子磁场的全球变化[J]. 地球物理学报, 1999,42(2):169-177.
|
[9] |
An Z C, Wang Y H. Global changes of the non-dipole magnetic fields for 1900—2000[J]. Chinese Journal of Geophysics, 1999,42(2):169-177.
|
[10] |
林云芳, 曾小苹, 郭启华. 东亚地区地磁非偶极场长期变化的分析[J]. 地球物理学报, 1985,28(5):482-496.
|
[10] |
Lin Y F, Zeng X P, Guo Q H. Analysis of secular variations of non-dipole geomagnetic field in east asia[J]. Chinese Journal of Geophysics, 1985,28(5):482-496.
|
[11] |
Matteo N A, Morton Y T. Ionosphere geomagnetic field: Comparison of IGRF model prediction and satellite measurements 1991—2010[J]. Radio Science, 2011,46(4):1-10.
|
[12] |
Xu W Y. Revision of the high-degree Gauss coefficients in the IGRF 1945—1955 models by using natural orthogonal component analysis[J]. Earth,Planets and Space, 2002,54(7):753-761.
|
[13] |
Maus S, Fairhead D, Hemant K, et al. A near-surface geomagnetic field model to spherical harmonic degree 720[C]// AGU Fall Meeting Abstracts, 2006.
|
[14] |
张素琴, 杨冬梅, 李琪, 等. 中国部分地磁台站年均值与IGRF模型一致性分析[J]. 地震地磁观测与研究, 2008,29(2):42-49.
|
[14] |
Zhang S Q, Yang D M, Li Q, et al. The consistence analysis of IGRF model value and annual mean value of some geomagnetic observatories in China[J]. Seismological and Geomagnetic Observation and Research, 2008,29(2):42-49.
|
[15] |
Molina F, de Santis A. Considerations and proposal for a best utilization of IGRF over areas including a geomagnetic observatory[J]. Physics of the Earth and Planetary Interiors, 1987,48(3-4):379-385.
|
[16] |
任国泰. 关于东亚大陆磁场的研究[J]. 地球物理学报, 1981(4):404-414.
|
[16] |
Ren G T. A comparison between regional model and global model of the geomagnetic field[J]. Geophysical and Geochemical Exploration, 1981,?(4):404-414.
|
[17] |
安振昌. 地磁场区域模型与全球模型的比较和讨论[J]. 物探与化探, 1991,15(4):248-254.
|
[17] |
An Z C. A comparison between regional model and global model of the geomagnetic field[J]. Geophysical and Geochemical Exploration, 1991,15(4):248-254.
|
[18] |
Smart D F, Shea M A, Flückiger E O. Magnetospheric models and trajectory computations[J]. Space Science Reviews, 2000,93(2):305-333.
|
[19] |
Smart D F, Shea M A. The limitations of using vertical cutoff rigidities determined from the IGRF magnetic field models for computing aircraft radiation dose[J]. Advances in Space Research, 2003,32(1):95-102.
|
[20] |
Davis J. Mathematical modeling of earth’s magnetic field[R]. Technical Note,Virginia Tech, Blacksburg, 2004.
|
[21] |
冯春. Matlab实现IGRF国际地磁参考场模型的计算[J]. 内蒙古石油化工, 2014,40(12):43-46.
|
[21] |
Feng C. International geomagnetic reference field model(IGRF)calculated by matlab[J]. Inner Mongolia Petrochemical Industry, 2014,40(12):43-46.
|
[22] |
柴松均, 陈曙东, 张爽. 国际地磁参考场的计算与软件实现[J]. 吉林大学学报:信息科学版, 2015,33(3):280-285.
|
[22] |
Chai S J, Chen S D, Zhang S. Calculation and software realization of international geomagnetic reference field[J]. Journal of Jilin University:Information Science Edition , 2015,33(3):280-285.
|
[23] |
杨梦雨, 管雪元, 李文胜. IGRF国际地磁参考场模型的计算[J]. 电子测量技术, 2017,40(6):97-100.
|
[23] |
Yang M Y, Guan X Y, Li W S. Calculation of IGRF international geomagnetic reference field model[J]. Electronic Measurement Technology, 2017,40(6):97-100.
|
[24] |
Campbell W H. Introduction to geomagnetic fields[M]. Cambridge: Cambridge University Press, 2003: 19-33.
|
[25] |
Balmino G, Barriot J, Koop R, et al. Simulation of gravity gradients: a comparison study[J]. Bulletin Geodesique, 1991,65(4):218-229.
|
[26] |
Abramowitz M, Stegun I A, Romer R H. Handbook of mathematical functions with formulas,graphs,and mathematical tables[M]. Gaithersburg:National Bureau of Standards Department of Commerce, 1988: 331-341.
|
[27] |
Blakely R J. Potential theory in gravity and magnetic applications[M]. London: Cambridge University Press, 1996: 100-127.
|
[1] |
CHEN Xiu-Juan, LIU Zhi-Di, LIU Yu-Xi, CHAI Hui-Qiang, WANG Yong. Research into the pore structure of tight reservoirs:A review[J]. Geophysical and Geochemical Exploration, 2022, 46(1): 22-31. |
[2] |
XIAO Guan-Hua, ZHANG Wei, CHEN Heng-Chun, ZHUO Wu, WANG Yan-Jun, REN Li-Ying. Application of shallow seismic reflection surveys in the exploration of urban underground space in Jinan[J]. Geophysical and Geochemical Exploration, 2022, 46(1): 96-103. |
|
|
|
|