2003年中国及邻区地磁场模型的计算与分析

根据2003年中国地区的135个测点和35个台站的地磁数据, 建立了2003年中国及邻区地磁场泰勒多项式模型和球冠谐模型. 在模型计算过程中, 细致地分析了模型的截断阶数和边界约束对模型计算的影响. 结果表明, 所使用的地磁观测资料是准确可靠的, 模型截断阶数的选取是合理的. 本文所建立的中国区域地磁场球冠谐模型与泰勒多项式模型具有良好的一致性.      

地磁正常场的选取与地磁异常场的计算

根据2003年中国地磁观测数据（包括135年地磁测点和35个地磁台）以及我国邻近地区38个IGRF计算点的地磁数据,计算中国地磁异常场的分布。选取两种地磁场模型作为地磁正常场,一是国际参考地磁场的球谐模型,二是中国地磁场泰勒多项式模型。根据各个测点的地磁异常值（观测值减去模型计算值）,用球冠谐分析方法计算地磁异常场的球冠谐模型,并绘制2003年中国地磁异常（△D,△I,△F,△X,△Y,△Z）。分析和讨论了中国地磁异常场。     

东亚地磁场模型的计算与分析

根据中国、前苏联和蒙古等地区的地磁复测点和地磁台站资料,使用泰勒多项式方法和冠谐分析方法,计算东亚地磁场(X,Y,Z)的泰勒多项式模型和冠谐模型,以及东亚剩余磁场(ΔX,ΔY,ΔZ)的冠谐模型和泰勒多项式模型,并绘制了相应的理论地磁图. 泰勒多项式模型的展开原点位于45°N和100°E,截断阶数为7.地磁场泰勒多项式模型的均方偏差为:X分量是133.0nT,Y分量是107.4nT,Z分量是14.0nT. 球冠极点位于45°N和100°E,球冠半角为42°,冠谐模型的截断阶数为10. 剩余磁场冠谐模型的均方偏差分别为131.2nT(ΔX),112.6nT(ΔY)和13.7nT(ΔZ). 对比分析了上述两种模型. 提出了确定区域模型截断阶数的判据.     

1936年中国地磁参考场的冠谐模型

根据1936年中国东部和中部地区的地磁测量资料，以及国际地磁参考场DGRF1935和DGRF1940，用冠谐分析方法计算1936年中国地磁参考场(CGRF)的冠谐模型. 冠谐模型的截断阶数为8，球冠极位于36°N和104°E，球冠半角为30°.CGRF1936比相应的国际地磁参考场能更好地表示中国地磁场的分布，本文计算的冠谐模型的均方偏差分别为104.9 nT(X分量)，84.8 nT(Y分量)和121.1 nT(Z分量). 根据地磁场的冠谐模型，绘制1936年中国地磁图(X，Y，Z，F)和异常磁场图(ΔX，ΔY，ΔZ，ΔF).     

京津冀地区地磁场球冠谐分析

2002年在京津冀地区进行了45个测点的地磁三分量测量,对测量资料进行通化处理,通化时间为2002年5月5日16～18时（世界时）. 通化后的观测均方差分别优于1.5nT（地磁场总强度F）,0.5′(磁偏角D和磁倾角I). 将国际参考地磁场（IGRF2000）作为地磁正常场,建立了京津冀地区地磁异常场的球冠谐模型（BTHASCH）. 球冠极的空间位置坐标为39.5°N和117.0°E,球冠半角为4°. 在模型计算过程中,球冠谐函数的截断阶数分别取为1～10. 经综合比较, 最终采用的截断阶数为5. 建立了京津冀地区参考地磁场的球冠谐模型（BTHGRF）.根据模型,绘制了京津冀地区地磁异常场图(ΔX、ΔY、ΔZ、ΔF、ΔD、ΔI)和京津冀地区地磁图(X、Y、Z、F、D、I).     

2000年中国地磁场及其长期变化冠谐分析

根据1998～2000年完成的118个地磁测点 和39个地磁台的三分量绝对测量资料以及IGRF2000，计算2000年中国地磁场冠谐模型(截断 阶数为8)，以及2000～2005年中国地磁长期变化冠谐模型(截断阶数为6). 球冠极位于36 °N，104°E，球冠半角为30°. 中国地磁场冠谐模型能更好地表示我国地磁场的时空变化 ，地磁场模型的均方偏差为：104.4 nT(X分量)，103.3 nT(Y分量)，123.9 nT(Z分量). 依据地磁场及其长期变化的冠谐模型，分别绘制2000年中国地磁图(F，X，Y，Z)和异常磁场图(ΔF，ΔX，ΔY，ΔZ)，以及2000～2005年地磁长期变化图(F，X，Y，Z). 指出改善地磁场模型边界效应 的途径，并对如何布设地磁复测点提出了建议.     

青藏高原地磁场模型的研究

根据青藏高原地磁三分量绝对测量资料,使用泰勒多项式方法和冠谐分析方法,计 算了青藏高原地磁场（Ｘ, Ｙ,Ｚ ）的泰勒多项式模型和青藏高原地磁剩余场（△Ｘ,△Ｙ,△Ｚ）的冠谐 模型,并绘制了相应的理论地磁图．分析了磁异常点对地磁场模型的影响,对比分析了地磁 场的多项式模型和冠谐模型,讨论了地磁场模型的边界效应问题．     

中国地区地磁场的球冠谐和分析

本文根据中国及其邻近地区的地磁三分量绝对测量资料,利用球冠谐和分析方法,计算出中国地区地磁剩余场的冠谐模型,地磁剩余场△X,△Y,△Z的模型均方根偏差分别为106.9,89.7,137.6 nT.提出由地磁场的国际参考地磁场和地磁剩余场的冠谐模型组成的联合模型作为中国参考地磁场的模型,它能较好地表示中国地磁场的分布.以地磁场的联合模型为正常背景场,计算出三分量磁异常的冠谐模型,并分析了磁异常的基本特征,它将为深入研究中国岩石层结构提供新证据.     

1950～1990年中国地磁剩余场冠谐分析

根据中国1950、1960、1970、1980年和1990年地磁三分量绝对测量资料，使用球冠谐和分析方法，分别计算1950～1990年各个年代中国地磁剩余场冠谐模型. 球冠极点位于36°N和104°E，球冠半角为30°，冠谐模型的截断阶数为8. 地磁剩余场冠谐模型的均方偏差分别为：对于X分量，1950年为93.1 nT，1960年为128.9 nT，1970年为107.2 nT，1980年为107.6 nT，1990年为95.2 nT；对于Y分量，1950～1990年依次为74.8 nT，98.1 nT，89.2 nT，89.9 nT和84.0 nT；对于Z分量，分别为122.2 nT(1950年)，135.0 nT(1960年)，137.7 nT(1970年)，110.1 nT(1980年)和107.5 nT(1990年). 根据中国地磁剩余场冠谐模型和全球地磁场DGRF模型，得到中国地磁场的冠谐模型，并对冠谐模型的边界效应进行了分析和讨论.     

欧洲及其邻区MAGSAT卫星磁异常冠谐模型

通过对欧洲及其邻近地区 MAGSAT 卫星黎明资料进行了处理,得到1°×1°卫星磁异常网格值。本文使用冠谐分析方法,计算该地区卫星矢量磁异常（ΔX,ΔY,ΔZ）冠谐模型。球冠极点位于33°N和26°E,球冠半角为40°.冠谐模型的截断指数为18。根据卫星磁异常冠谐模型和地磁场球谐模型 DGRF1980,计算卫星总强度磁异常（ΔF）的冠谐一球谐模型。根据卫星磁异常的理论模型,计算并绘制不同高度（300,400,500km）的理论卫星磁异常图。对冠谐模型和大磁异常进行了分析和讨论。     

1683-1949年中国地磁测量、地磁图和地磁模型的总考评与研究

本文全面系统回顾考察了新中国成立前中国的地磁测量、地磁图和地磁模型研究概况:主要介绍中国人独立进行的地磁测量;评述了1932年山东半岛地磁图,1915.0和1936.0年中国(部分地区)地磁图;1945.0年东亚地磁图;1946.1年四川北碚地区地磁图以及1946.1年四川北碚地区地磁场泰勒多项式模型.并根据1936年地磁数据计算研究了1936.0年中国地磁场的球冠谐模型和曲面Spline模型,以及1900.0-1940.0年5个年代的中国地磁数据集、地磁模型和地磁图.对今后我国地磁测量研究工作给出了几点建议.     

基于三维Taylor模型的福建及邻近地区地磁模型研究

基于福建及邻近地区高密度的256个地磁场矢量(D、I、H分量)实测数据,结合全球陆地1km基准海拔高度工程的网格数据以及最新的第十二代国际地磁参考场模型(IGRF12),采用三维Taylor多项式模型构建了福建地区的地磁场模型.通过比较均方偏差(RMS)、磁场实际分布以及残差等,结果显示:(1)当三维Taylor模型的截断阶数(N)为2时,可基本反映N为6时的二维Taylor模型的拟合效果;(2)三维Taylor模型计算方便,精度较高,但较易出现龙格现象,因此在确定截断阶数时既要考虑模拟精度,还需要考虑边界效应问题.     

地磁场模型和冠谐分析

地磁场是由内源场(地核场和地壳场)和外源场(电离层场和磁层场)组成的,地磁场模型是表示地磁场时空分布的数学表达式.本文简述了地磁场全球模型和区域模型的计算方法和应用;评述了冠谐分析方法及其在地磁学中的应用;介绍了地磁场模型误差的主要来源.     

基于CHAMP卫星与三维Taylor多项式模型的区域地磁建模研究

利用CHAMP卫星矢量和标量地磁测量数据,通过三维Taylor多项式模型建立了2010.0年中国及邻近地区在300km高度附近的X、Y、Z和F分量的三维磁场模型.为了比较验证所建模型,分析了其与Taylor多项式模型(二维模型)不同截断阶数所对应的均方偏差(RMSE)、残差及分布等.结果表明,三维模型每一阶的RMSE和残差绝对值的平均值均要比二维模型的小约45%.由于采用了系数完全展开的建模方式,三维模型的系数数量约为二维模型的2倍,三维模型较低的截断阶数可以反映更多的地磁信息.本研究中的5阶三维Taylor模型基本可达到8阶Taylor模型的精度.两种模型绘制的地磁场及残差分布有较好的一致性.     

基于CHAMP卫星与地面矢量数据联合建立中国地磁模型

首先对地面台站以及CHAMP卫星的矢量磁测数据进行预处理,结合12个子午工程台站最新的磁测数据,并利用第十二代国际地磁参考场模型(IGRF12)对不同数据归算至统一时间点,再结合第四代地磁场综合模型(CM4)计算并移除主磁场、磁层磁场、电离层磁场、感应磁场及较为微弱的环形场,最后通过球冠谐(SCH)模型联合建立了中国地区岩石圈磁场(地磁要素X、Y、Z和总强度F)的球冠谐模型(SCH2000),结合IGRF12得到了中国地磁模型(CGRF2000).建模同时考虑了测点数据海拔高度的影响,并讨论了截断阶数的选取以及边界效应的控制问题.通过比较均方偏差(RMS),认为8阶SCH2000模型可较好地反映中国地区岩石圈磁场,其分布与SCH1936、CM4以及Taylor多项式模型具有一定的相似性.从位置和强度来看,SCH2000模型所反映的岩石圈磁场更为准确,而与CM4存在一定差异,除了测点不同,截断阶数也是原因之一.通过验证,SCH2000模型各要素随高度的变化与CM4模型较为一致,所建立的CGRF2000与IGRF12模型的分布也较为相似.     

喀什—乌恰地区流动地磁测量数据处理及初步分析

介绍了2008年在喀什—乌恰地区32个测点的流动地磁测量数据的处理方法。利用该数据建立了喀什—乌恰地区地磁场的曲面样条模型,绘制了区域地磁基本场分布图。以2005年中国地磁场球冠谐模型为参考场剥离了岩石圈磁场,绘制了喀什—乌恰地区岩石圈磁场局部异常分布图。表明岩石圈磁场局部异常与测区内地震在空间上有一定的对应关系。     

计算震磁背景场的数学方法

如何计算地磁场 (或地磁剩余场 )及其长期变化的趋势变化 ,是提取震磁异常的关键。文章简要介绍了计算震磁背景场的数学方法 :多项式方法、球冠谐和分析方法、矩谐分析方法和曲面样条函数方法     

2010.0年中国及邻近地区地磁场

本文通过对中国地区地磁场的球冠谐和分析和曲面样条分析,建立了2010.0年代的"中国地磁参考场球冠谐和模型"和"中国地磁参考场曲面样条模型".得到如下结论:"2010.0年代中国地磁参考场球冠谐和模型"描述的各地磁要素空间分布与"2005.0年代中国地磁参考场球冠谐和模型"较相似,但局部略有变化.比如对于磁偏角D和东向分量Y大致可理解为负值区域增强;磁倾角I在我国华北、华南和东北大部大面积正异常区消失了,但在我国西北部大面积的负异常区仍然维持;总强度F和垂直分量Z整体的负异常分布仍然维持,新疆西部局部地区和东南沿海地区局部的正异常分布则略有变化;变化最大的是水平分量H和北向分量X的空间分布,在2010.0年代,东南地区出现了大范围的正异常区,而且在中国北部出现了一条东西条带状的正异常区,似乎可认为是2005.0年代存在于太平洋的正异常区向西移动扩大的结果."2010.0中国地磁参考场曲面样条模型"描述的磁偏角D和东向分量Y相对较为曲折,最显著的异常位于蒙古境内,其次为滇东北地区、东北及华北大部、新疆等地区;磁倾角I和垂直分量Z相对较为平缓,最显著的异常存在于蒙古境内,其次的异常在内蒙满洲里附近;总强度F、水平分量H和北向分量X最为显著的异常位于满洲里附近,其次的主要异常区位于新疆地区、京津附近、滇东北地区.     

利用GPS数据与Klobuchar模型估计地磁北极

基于地磁坐标和地理位置的关系,利用Newton迭代法推导了地磁北极的计算公式;根据Klobuchar模型中电离层电子含量与地磁位置的关系,给出地磁位置的表达式;进而建立了GPS观测量估计地磁北极的模型.理论分析和数据验证了方法的可行性,得出1)Klobuchar模型属于全球性的电离层估计模型,对于区域性电离层信息描述的不十分准确,建立区域化Klobuchar模型将更有利于地磁北极的估计;2)利用泰勒公式将Klobuchar模型线性化的过程中,公式仅取其一阶项,将会影响模型的精度,有待于提高模型的展开阶数;3)实际中地磁北极的估计结果受初始值得约束;4)不同计算方法和模型对结果的可靠性有所影响,试图选择其它估计过程和方法,提高估计的可靠性.