Characteristics of Spatial and Temporal Distribution of Geomagnetic Field in Chinese Mainland Based on CM4 Model
-
摘要: 本文利用第四代地磁场综合模型（Comprehensive Model 4，CM4），计算了1982-2001年中国大陆地区同一经度链和同一纬度链上地磁台站的磁层源磁场及其感应场、电离层源磁场及其感应场的地磁北向分量X、东向分量Y、垂直分量Z的模型值，分析了各场源磁场随时间和空间的变化特征。结果表明：在时间上，经度链和纬度链台站的磁层源磁场及其感应场均呈现出11年和27天周期性变化。电离层源磁场及其感应场具有明显的季节变化，不同年份相同季节变化形态一致但幅度不同。在空间分布上，经度链和纬度链台站磁层源磁场及其感应场的年变化幅度呈现出不同变化特征，电离层源磁场及其感应场在经度链上变化特征不同，而纬度链台站的数值基本一致。日变化分析显示，磁静日和磁扰日期间，模型数据与台站实测数据变化一致性较好，相关性较高。Abstract: In this paper, based on the comprehensive model4 (CM4), we calculated the latitude and longitude chain stations magnetosphere and its induction magnetic field, ionosphere and its induction magnetic field of the X, Y, Z component from 1982 to 2001 in chinese mainland and also analyzed the change characteristics of the magnetic field with time and space. The results show that the magnetosphere and its induction magnetic field of the longitude chain and latitude chain show periodic variation in 11 years and 27 days. The ionosphere and its induction magnetic field have shown seasonal variation, morphological changes in the same season in which the amplitude is different in the spatial distribution. The annual variation range of the magnetosphere and its induction magnetic field in the longitude chain and latitude chain showed different characteristics. The ionosphere field and its induction magnetic field vary spatially in the longitude chain, but they are about the same of latitude chain station. Diurnal variation analysis showed that the model data is in good consistency with the measured data of the station on the magnetically quiet and disturbed date, and they are highly correlated.
-
Key words:
- CM4 /
- Magnetosphere field /
- Ionosphere field /
- Observatory data /
- Induction field /
- The mean time value /
- Correlation coefficient
-
引言
在实际地震中，地震动输入具有一定的不确定性，钢筋混凝土框架结构需要经受双向地震作用的考验。在双向受力下，节点核心区存在较大的水平剪力，产生的破坏往往是剪切破坏。通过以往的震害研究可以看出，节点的破坏是框架结构发生破坏的主要原因之一。因此，如何保证节点具有足够的抗震能力成为亟待解决的问题。
各国学者针对平面边节点展开了一系列研究，唐九如等(1985)、Ghobarah等(2005)通过大量边节点试验研究，发现普通边节点抗剪承载力低于中节点，并提出边节点存在混凝土斜压杆。针对空间节点，陈永春等(1995)、雷远德等(2014)通过研究，认为双向往复荷载作用较单向往复荷载作用更容易造成框架柱柱端或节点区破坏，承载力、延性等抗震性能退化明显，且同样出现了边柱节点在节点区附近发生脆性破坏，边柱节点承载力明显低于中柱节点的情况(框架节点专题研究组，1983；傅剑平，2002)。Hwang等(2000)基于已有研究提出了应用于平面节点的软化拉-压杆模型，这一模型是在压杆-拉杆模型的基础上考虑混凝土受压开裂后的软化特性提出的，用以计算节点区抗剪强度。但目前关于软化拉-压杆模型在空间边节点抗剪承载力计算的研究尚未完善。
因此，本文以水平双向受力相同的空间边节点为研究对象，在软化拉-压杆模型平面边节点抗剪计算模型的基础上，建立空间边节点抗剪计算方法，并通过数值模拟的方法对计算方法进行验证。
1. 空间边节点抗剪承载力计算模型
1.1 空间边节点受力分析
由于空间边节点构造的特殊性，与空间中节点相比，其一端缺少梁约束。节点区内的柱筋及梁筋造成节点上下侧的柱端、3个方向的梁端分别受拉或受压，同时受压混凝土形成相应受压区。在双向受力下，节点需要承受2个方向的力，如图 1所示，其中，Mc、Nc、Vc分别为柱端弯矩、轴力和剪力；Mbx、Nbx、Vbx分别为x向梁端弯矩、轴力和剪力；Mby、Nby、Vby分别为y向梁端弯矩、轴力和剪力。节点需要承受混凝土传递的压力、梁筋传递的力以及柱筋传递的力。节点核心区内，在2个方向分别形成节点水平剪力Vhx和Vhy，通过合成得到水平的斜向合成剪力Vh，2个方向的竖向剪力在节点核心区内合成为竖向合成剪力Vv。空间边节点抗剪机制如图 2所示，其中，hc、bs、hb分别为柱高、节点宽和梁高，${{h}_{c}}^{\prime \prime }$、${{h}_{s}}^{\prime \prime }$、${{h}_{b}}^{\prime \prime }$为相应侧柱筋及梁筋中心之间的距离。
本文主要针对空间节点中水平双向受力相同的情况进行研究，因此水平合剪力Vh与水平方向夹角为45°，竖向合剪力Vv垂直于水平面。节点核心区合剪力作用面内，由混凝土形成压杆、钢筋形成拉杆，组成拉-压杆模型承担节点合剪力。因此，在空间节点的合剪力作用面内建立拉-压杆模型。其中，θ为斜压杆倾角，α为合剪力作用面与y向主轴的夹角。
1.2 建立计算模型
在研究双向受力下节点抗剪计算时，主要针对2个方向受力相同这一情况。斜向混凝土斜压杆与水平合剪力、竖向合剪力均作用在同一平面内，在此作用面内建立节点的抗剪计算模型，该合剪力作用面与主轴夹角为45°。双向受力下的软化拉-压杆模型如图 3所示，合剪力作用面内虚线代表混凝土形成的压杆，实线代表钢筋转化成的拉杆。D为节点区斜压杆承担的斜向压力，Fh为节点内水平箍筋所形成的水平拉杆的拉力，Fv为节点内抗剪柱筋形成的竖向拉杆的拉力。
空间边节点的水平抗剪承载力Vjh表示为(Hwang等，2002)：
$$ {{V}_{jh}}=K\xi {{{f}'}_{c}}{{A}_{str}}\text{cos}\theta $$ (1) 其中，K为拉压杆系数；ξ为混凝土软化系数；${{{f}'}_{c}}$为混凝土抗压强度；Astr为斜压杆有效截面面积；θ为斜压杆与水平夹角。
节点核心区的剪力主要由混凝土形成的压杆承担，节点内的水平箍筋及竖向中间柱筋使核心区更多的混凝土形成压杆用，使得节点核心区在主斜压杆方向能够承担更多的压力，从而达到提高节点抗剪承载力的作用。因此，当节点内存在水平箍筋及竖向中间柱筋形成的拉杆时，斜向受压在斜压杆传递基础上增加了额外的传力路径，引发了更多核心区混凝土参加抗剪，从而增大了节点区抗剪强度。在计算中，将拉杆的贡献用系数K表示，计算公式为：
$$ K={{K}_{h}}+{{K}_{v}}-1 $$ (2) 其中，Kh为水平拉杆系数；Kv为竖向拉杆系数，可分别通过式(3)、(4)进行计算：
$$ {{K}_{h}}=1+({{\overline{K}}_{h}}-1)\times \frac{{{F}_{yh}}}{{{\overline{F}}_{h}}}\le {{\overline{K}}_{h}} $$ (3) $$ {{K}_{v}}=1+({{\overline{K}}_{v}}-1)\times \frac{{{F}_{yv}}}{{{\overline{F}}_{v}}}\le {{\overline{K}}_{v}} $$ (4) $$ {{\bar{K}}_{h}}=\frac{(1-{{\gamma }_{h}})+{{\gamma }_{h}}}{(1-{{\gamma }_{h}})+{{\gamma }_{h}}\left( 1-\frac{\text{si}{{\text{n}}^{2}}\theta }{2} \right)}\ge 1 $$ (5) $$ {{\bar{K}}_{v}}=\frac{(1-{{\gamma }_{v}})+{{\gamma }_{v}}}{(1-{{\gamma }_{v}})+{{\gamma }_{v}}\left( 1-\frac{\rm{co}{{\rm{s}}^{\rm{2}}}\theta\rm{ }}{\mathit{2}} \right)}\ge 1 $$ (6) $$ {{\gamma }_{h}}=\frac{2\tan \theta-1}{3}(0\le {{\gamma }_{h}}\le 1) $$ (7) $$ {{\gamma }_{v}}=\frac{2\cot \theta-1}{3}(0\le {{\gamma }_{v}}\le 1) $$ (8) 其中，${{\overline{K}}_{v}}$为水平拉杆平衡系数、${{\overline{K}}_{v}}$为竖向拉杆平衡系数；γh为水平拉杆的拉力与节点水平剪力的比值、γv为竖向拉杆的拉力与节点竖向剪力的比值。
水平及竖向拉杆的拉力可分别用式(9)、(10)进行计算：
$$ {{\overline{F}}_{h}}={{\gamma }_{h}}\times ({{\overline{K}}_{h}}\xi {{{f}'}_{c}}{{A}_{str}})\times \text{cos}\theta $$ (9) $$ {{\overline{F}}_{v}}={{\gamma }_{v}}\times ({{\overline{K}}_{v}}\xi {{{f}'}_{c}}{{A}_{str}})\times \text{sin}\theta $$ (10) 其中，${{\overline{F}}_{h}}$为水平拉杆拉力、${{\overline{F}}_{v}}$为竖向拉杆拉力。
在水平机构计算中，定义节点箍筋所提供的抗剪力作用在核心区45°方向，在计算箍筋面积时，认为空间节点箍筋的有效面积为平面节点箍筋有效面积的$\sqrt{2}$倍。因此，Fyh=$\sqrt{2}$Athfyh，Ath为节点水平箍筋各肢的总截面面积，fyh为节点水平箍筋的屈服强度。竖向机构主要由柱纵筋提供抗剪承载力，在空间模型中2个方向的剪力合力均为竖直方向。因此，Fyv=Atvfyv，Atv为节点区内抗剪柱筋截面面积，fyv为节点中间柱筋的屈服强度。
节点所能抵抗的最大水平剪力主要由主斜压杆端部所能抵抗的最大压力决定，而主斜压杆端部的抗压能力由混凝土抗压强度决定。Hwang等(2002)通过研究发现开裂后的钢筋混凝土抗压强度远小于单轴受压混凝土的抗压强度，同时混凝土的抗压强度会随柱拉应变的增大而减小，发生了软化效应。在节点抗剪过程中，将节点主斜压杆面积范围内混凝土达到极限的抗压强度，确定为节点抗剪失效的标准。因此，在计算节点抗剪承载力时需要考虑混凝土抗压强度软化系数ξ，其计算公式为：
$$ \xi \approx \frac{3.35}{\sqrt{{{{{f}'}}_{c}}}}\le 0.52 $$ (11) 其中，${{{f}'}_{c}}$为混凝土圆柱体抗压强度。
根据欧洲CEB-FIP规范(Euro-International Concrete committee，1993)，混凝土圆柱体抗压强度与混凝土立方体抗压强度标准值之间的换算关系为${{{f}'}_{c}}$=0.80fcu, k；根据中国《混凝土结构设计规范》(GB50011—2010)(中华人民共和国建设部，2010)，fck=0.88αc1αc2fcu, k, fck为混凝土轴心抗压强度标准值。
1.3 斜压杆的确定
空间边节点在受压区混凝土的作用下，节点核心区内部形成了斜向受压杆(崔建宇等，2014)。梁端及柱端传入节点核心区的压力在柱端上、下部柱端角部形成了受压区，两端受压区所形成的斜三角柱为斜压杆范围。定义下表面(ΔLON)为等效斜压杆横截面。下表面中N、L点的位置可通过2个方向梁受压区宽度确定。根据斜压杆受压特点，假设等效斜压杆横截面与节点斜对角线垂直，L、O、N点相连与水平面形成一定角度的三角形。定义ΔLON的面积为斜压杆有效抗压面积Astr，核心区立方体对角线与地面形成的夹角为斜压杆角度θ。
斜压杆角度θ及斜压杆有效抗压面积Astr可表示为：
$$ \theta =\text{arctan}\frac{{{{{h}''}}_{b}}}{\sqrt{{{{{h}''}}_{c}}^{2}+{{{{b}''}}_{s}}^{2}}} $$ (12) $$ {{A}_{str}}={{A}_{LON}}=\frac{{{A}_{LCN}}}{\text{sin}\theta } $$ (13) 其中，ALON、ALCN分别为ΔLON、ΔLCN的面积。
lc为柱下端受压区ΔLCN直角边长度，lc可按式(14)计算；ΔLON的面积可通过式(15)计算，即可得到空间边节点斜压杆有效抗压面积Astr。
$$ {{l}_{c}}=\left(\text{0}\text{.85}+\text{0}\text{.25}\frac{N}{{{A}_{\text{g}}}{{{{f}'}}_{c}}} \right){{h}_{c}} $$ (14) $$ {{A}_{str}}={{A}_{LON}}=0.707{{l}_{c}}\sqrt{0.5l_{c}^{2}+a_{b}^{2}} $$ (15) 2. 计算模型的验证
利用有限元分析软件ABAQUS进行数值模拟。模型的混凝土单元采用C3D8R实体单元；钢筋单元根据钢筋受力特点，采用桁架单元T3D2。模型边界条件柱顶端为铰接，柱底为完全铰接，梁端为加载自由端。为验证数值模拟的准确性，对王玉雷(2010)研究中的试验试件JD4进行数值分析。
表 1、图 5为模拟结果与试验结果对比及荷载位移曲线模拟值与试验值对比图，通过对比分析可知有限元模拟结果与试验结果吻合良好。
表 1 有限元模拟值与试验值对比Table 1. Comparison of simulated results and test results对比项 屈服荷载/kN 屈服位移/mm 峰值荷载/kN 峰值位移/mm 极限荷载/kN 极限位移/mm 试验值 139.9 13 163.9 45 152.8 60 模拟值 148.4 11.5 161.8 41.8 137.6 54.6 
图 5 荷载位移曲线模拟值与试验值对比Figure 5. Comparison of load-displacement curve between test results and simulated results为验证计算方法的准确度，建立了6个空间边节点模型，模型具体信息见图 6。在模型中，对柱施加轴向压力，同时在梁端双向加载，分别施加向上(向下)的力，模拟实际情况中边节点在双向受力下的受力过程。同时，将空间角节点这一特殊边节点情况考虑在内，因此加入3个角节点试验试件，并将计算结果与模拟结果或试验结果进行对比。
将计算值与模拟得到的节点抗剪承载力模拟值进行对比，模型SP1—SP3为崔建宇等(2014)给出的3个空间角节点试验试件，对其进行抗剪承载力计算，并将计算结果与试验值进行对比，具体数值见表 2。通过表 2的数据可以看出，计算值与模拟值或试验值较为接近。
表 2 空间边节点参数及计算值与模拟值对比Table 2. Comparison of spatial joints between calculated results and simulated results模型编号 节点形式 hb/mm hc/mm fck/MPa 轴压比 柱纵筋配筋 箍筋配筋 计算值Vjh, c/kN 模拟值Vjh, a/kN Vjh, c/Vjh, a B-1 边节点 450 400 27.9 0.1 12C18+4C22 A10@50 474.17 550.09 0.862 B-2 边节点 450 400 27.9 0.1 12C18+4C22 A12@50 484.77 573.71 0.845 B-3 边节点 450 400 27.9 0.15 12C18+4C22 A10@50 486.51 552.99 0.879 B-4 边节点 450 400 27.9 0.3 12C18+4C22 A10@50 524.57 589.44 0.890 B-5 边节点 450 400 27.9 0.3 12C18+4C22 A12@50 535.17 602.73 0.888 B-6 边节点 450 400 27.9 0.45 12C18+4C22 A10@50 564.20 616.74 0.915 SP1 角节点 250 220 34.3 0.2 16D13 A6@80 216.18 262* 0.825 SP2 角节点 250 220 29.1 0.6 16D13 A6@80 225.73 247.8* 0.911 SP3 角节点 250 220 29.1 0.6 16D13 A10@80 253.41 275.8* 0.919 注：fck为混凝土轴心抗压强度标准值；Vjh, a值中上标*的为试验值。 3. 结论
(1) 通过对平面软化拉—压杆模型和双向受力下钢筋混凝土边节点受力机理分析，在节点核心区双向剪力合成作用面内，建立了双向受力下钢筋混凝土框架边节点抗剪承载力计算方法。
(2) 基于软化拉-压杆模型，建立了空间边节点抗剪计算方法，确定了空间边节点核心区斜压杆截面面积及倾角的取值方法。对双向受力下软化拉-压杆模型中的拉杆进行了定义。
(3) 利用有限元分析软件ABAQUS建立不同参数下空间边节点模型，得到相应抗剪承载力模拟值。将计算值与模拟值或试验值进行对比，结果表明，双向受力下钢筋混凝土框架边节点抗剪计算方法的准确度良好。
-
图 1 经度链台站磁层源磁场及其感应场X分量时均值随时间变化
Figure 1. Changes of the mean value of magnetosphere and its induction field of X component with time in the longitude chain
图 2 经度链台站X分量年变化幅度的空间分布
Figure 2. Spatial distribution of annual variation range of X component in longitude chain
图 3 经度链台站磁层源磁场及其感应场Y分量时均值随时间变化
Figure 3. Changes of the mean value of magnetosphere and its induction field of Y component with time in the longitude chain
图 4 经度链台站Y分量年变化幅度的空间分布
Figure 4. Spatial distribution of annual variation range of Y component in longitude chain
图 5 经度链台站磁层源磁场及其感应场Z分量时均值随时间变化
Figure 5. Changes of the mean value of magnetosphere and its induction field of Z component with time in the longitude chain
图 6 经度链台站Z分量年变化幅度的空间分布
Figure 6. Spatial distribution of annual variation range of Z component in longitude chain
图 7 纬度链台站X分量年变化幅度的空间分布
Figure 7. Spatial distribution of annual variation range of X component in latitude chain
图 8 纬度链台站Y分量年变化幅度的时空分布
Figure 8. Spatial distribution of annual variation range of Y component in latitude chain
图 9 纬度链台站Z分量年变化幅度的时空分布
Figure 9. Spatial distribution of annual variation range of Z component in latitude chain
图 10 1997年经度链台站X分量的27天周期变化
Figure 10. The 27-day cycle variation of X component in the longitudinal chain in 1997
图 11 1997年经度链台站Y分量的27天周期变化
Figure 11. The 27-day cycle variation of Y component in the longitudinal chain in 1997
图 12 经度链台站1991年电离层源磁场及其感应场变化
Figure 12. Variation of ionosphere and its induction field of longitudinal chain in 1991
图 13 纬度链台站1991年电离层源磁场及其感应场变化
Figure 13. Variation of ionosphere and its induction field of latitude chain in 1991
图 14 1997年7月模型数据与实测数据的时均值变化
Figure 14. Changes of hourly mean values and Dst in Jul.1997
图 15 2001年3月模型数据与实测数据的时均值对比
Figure 15. Changes of hourly mean values and Dst in Mar., 2001
表 1 台站分布
Table 1. Location of the stations
经度链 纬度链 台站 代码 高度/m 台站 代码 高度/m 满洲里 MZL 682 格尔木 GLM 2802 红山 LYH 20 兰州 LZH 1560 武汉 WHN 45 红山 LYH 20 泉州 QZH 41 泰安 TAA 249 
下载: 导出CSV
表 2 经度链台站X分量年变化幅度
Table 2. Annual variation range of X component in longitude chain
YEAR MZL LYH WHN QZH YEAR MZL LYH WHN QZH 1982 227.2 256.0 267.4 272.9 1992 185.4 211.7 222.4 228.9 1983 133.6 153.7 162.0 167.0 1993 108.1 122.6 128.5 131.4 1984 93.6 106.1 111.2 114.2 1994 93.1 107.2 112.9 115.8 1985 92.9 107.3 113.2 116.3 1995 131.1 149.7 157.4 161.3 1986 114.2 131.2 138.2 141.7 1996 46.7 54.3 57.3 58.6 1987 85.3 94.5 97.9 99.6 1997 72.8 83.2 87.4 90.0 1988 126.1 144.1 151.7 156.1 1998 117.1 132.0 138.7 142.6 1989 213.6 245.5 258.7 266.1 1999 106.6 120.4 125.8 128.5 1990 229.6 262.4 275.9 283.3 2000 138.3 158.2 166.2 170.4 1991 172.7 194.8 203.4 208.2 2001 238.1 274.0 289.2 298.3
下载: 导出CSV
表 3 实测值与模型值差值的标准差
Table 3. The standard deviation of the difference between the measured and model values
分量 MZL LYH WHN QZH GLM LZH TAA X 11.2 10.6 8.1 12.8 5.6 11.3 32.1 Y 4.4 4.8 2.0 12.6 6.8 10.7 4.7 Z 21.4 33.9 5.9 34.8 8.9 21.1 24.7
下载: 导出CSV
-
白春华, 徐文耀, 康国发, 2008.地球主磁场模型.地球物理学进展, 23(4):1045-1057. http://manu39.magtech.com.cn/Geoprog/CN/abstract/abstract1748.shtml 丁鉴海, 索玉成, 余素荣, 2005.地磁场与电离层异常现象及其与地震的关系.空间科学学报, 25(6):536-542. doi: 10.11728/cjss2005.06.536 冯彦, 安振昌, 孙涵等, 2011.基于第四代地磁场综合模型(CM4)以改善模型边界效应的研究.地球物理学进展, 26(3):850-857. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract3446.shtml 李细顺, 高登平, 李琪等, 2015.CM4模型数据与台站实测数据的对比研究.震灾防御技术, 10(2):418-425. doi: 10.11899/zzfy20150223 李琪, 杨星, 蔡绍平, 2015.极化方法应用于地磁台阵的震例分析.震灾防御技术, 10(2):412-417. doi: 10.11899/zzfy20150222 王亶文, 2003.国际地磁参考场在中国大陆地区的误差分析.地球物理学报, 46(2):171-174. http://www.cqvip.com/QK/94718X/200302/7519358.html 徐文耀, 2002.地球主磁场的NOC模型.中国科学(D辑), 32(7):576-587. doi: 10.3321/j.issn:1006-9267.2002.07.006 徐文耀, 2009.地球电磁现象物理学.合肥:中国科学技术大学出版社, 18-29. 徐文耀, 2011.太阳风-磁层-电离层耦合过程中的能量收支.空间科学学报, 31(1):1-14. doi: 10.11728/cjss2011.01.001 杨云存, 高国明, 2014.1900-2010年地磁场水平分量梯度的全球变化.震灾防御技术, 9(S):557-571. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=2014s101&journal_id=zzfyjs 姚丽, 刘振兴, 左平兵等, 2010.行星际激波对地球磁层的压缩效应分析.空间科学学报, 30(2):113-120. doi: 10.11728/cjss2010.02.113 张素琴, 杨冬梅, 李琪等, 2008.中国部分地磁台站年均值与IGRF模型一致性分析.地震地磁观测与研究, 29(2):42-49. http://www.cqvip.com/qk/96509X/200802/27140315.html Alldredge L. R., 1987. On regional magnetic charts. Journal of Geomagnetism and Geoelectricity, 39:723-738. doi: 10.5636/jgg.39.723 Haines G. V., 1985. Spherical cap harmonic analysis. Journal of Geophysical Research:Solid Earth, 90(B3):2583-2591. doi: 10.1029/JB090iB03p02583 Hemant K., Maus S., 2005. Geological modeling of the new CHAMP magnetic anomaly maps using a geographical information system technique. Journal of Geophysical Research:Solid Earth, 110(B12):B12103. doi: 10.1029/2005JB003837 Hemant K., Thébault E., Mandea M., et al., 2007. Magnetic anomaly map of the world:merging satellite, airborne, marine and ground-based magnetic data sets. Earth and Planetary Science Letters, 260(1-2):56-71. doi: 10.1016/j.epsl.2007.05.040 Langel R. A., Sabaka T. J., Baldwin R. T., et al., 1996. The near-earth magnetic field from magnetospheric and quiet-day ionospheric sources and how it is modeled. Physics of the Earth and Planetary Interiors, 98(3-4):235-267. doi: 10.1016/S0031-9201(96)03190-1 Lesur V., Wardinski I., Rother M., et al., 2008. GRIMM:the GFZ reference internal magnetic model based on vector satellite and observatory data. Geophysical Journal International, 173(2):382-394. doi: 10.1111/gji.2008.173.issue-2 Sabaka T. J., Olsen N., Langel R. A., 2002. A comprehensive model of the quiet-time, near-Earth magnetic field:phase 3. Geophysical Journal International, 151(1):32-68. doi: 10.1046/j.1365-246X.2002.01774.x Sabaka T. J., Olsen N., Purucker M. E., 2004. Extending comprehensive models of the Earth's magnetic field with ørsted and CHAMP data. Geophysical Journal International, 159(2):521-547. doi: 10.1111/gji.2004.159.issue-2 Wardinski I., Holme R., 2006. A time-dependent model of the earth's magnetic field and its secular variation for the period 1980-2000. Journal of Geophysical Research:Solid Earth, 111(B12):B12101. doi: 10.1029/2006JB004401/full -
下载:
点击查看大图
计量
- 文章访问数: 114
- HTML全文浏览量: 19
- PDF下载量: 3
- 被引次数: 0
