磁方位角与真方位角的关系.docxVIP

磁方位角与真方位角的关系 1 磁偏角随经度、自由度及高程的变化规律 直线方向是确定地球直线与标准方向之间的角度。就不同的标准方向（真子普通话方向和磁子ub短路方向），它们之间的角度不同，直线方向也不同。基于磁北方向的磁方位角与真北方向的真方位角在数值大小上存在差异, 即磁北方向与真北方向之间存在一定的夹角, 这一夹角称磁偏角。由于地磁会随着时间、地点的变化而变化, 从而决定了磁偏角的变化性。首先, 讨论了几种方位角的关系;然后基于地磁理论, 推导了在已知大地坐标时磁偏角的计算公式, 并对我国经度、纬度范围内对应的磁偏角和磁倾角进行了求解, 分析了磁偏角和磁倾角随经度、纬度及高程的变化规律。这对于磁偏角的估算和实际应用具有一定的现实意义。 2 磁偏角面为轴子古代角 直线的定向就是确定地面直线与标准方向的夹角, 其标准方向可分为: 1) 真子午线方向, 即地面任一点在真子午线处的切线方向; 2) 磁子午线方向, 即地面任一点在磁子午线处的切线方向; 3) 轴子午线 (坐标纵轴) 方向, 即地面任一点与高斯平面直角坐标系或假定坐标系的坐标纵轴平行的方向。磁偏角就是真、磁子午线的夹角δ。由子午线北端顺时针方向量到测线上的夹角, 由于标准方向的不同, 产生了不同的方位角, 有:真方位角α, 磁方位角β, 坐标方位角γ。各方位角的关系表示为: α=β+δ=γ+θ(1) 式 (1) 中,θ为子午线收敛角, 由于真子午线与轴子午线的不重合而产生的夹角。其中真方位角与磁方位角之间的关系可用图1表示。 3 磁极角法的求解 3.1 地磁场强度 在地磁研究中, 对于地磁场强度以矢量B表示, 同时引入参考坐标系来描述这一矢量。对应坐标系定义为:原点O位于研究点, 直角坐标系中OX指向地理北,OY指向地理东,OZ垂直向下。B为地磁场总强度;BH为B在水平面内的投影, 称为水平强度或水平分量;X为H在OX轴上的投影, 称为北向强度或北向分量;Y为H在OY轴上的投影, 称为东向强度或东向分量;Z为F在OZ轴上的投影, 称为垂直强度或垂直分量;D为H偏离OX轴即偏离地理北的角度, 称为磁偏角,H向东偏为正;I为B同水平面的夹角, 称为磁倾角,B向下倾为正。各地磁要素之间的关系如图2所示。 3.2 dgrf模型 国际参考地磁场 (IGRF) 是描述地球主磁场的标准全球模型, 在IGRF模型中, 主磁场的标量磁位可用球谐函数表达: V(r,θ,φ)=a∑k=0k∑m=0k(ar)n+1(gmkcosmφ+hmksinmφ)Pmk(θ)(2)V(r,θ,φ)=a∑k=0k∑m=0k(ar)n+1(gkmcosmφ+hkmsinmφ)Ρkm(θ)(2) 式中,a为地球半径,r为地球径向距离,gmkkm\,hmkkm为地磁场球谐系数 (或高斯系数) ,Pmkkm(θ) 为伴随勒让德函数,θ\,φ分别为经度和余纬,φ=90-λ(λ是纬度) ,k为模型的截断水平。相应的地磁分量X\,Y\,Z为: X=∑k=0kmax∑m=0k(ar)nk+1(gmkcosmφ+hmksinmφ)dPmnk(cosθ)dθY=∑k=0kmax∑m=1kmsinθ(ar)nk+1(gmksinmφ?hmkcosmφ)Pmnk(cosθ)Z=?∑k=0kmax∑m=1k(nk+1)(ar)nk+1(gmkcosmφ+hmksinmφ)Pmnk(cosθ)???????????????????????????????????????????????(3)X=∑k=0kmax∑m=0k(ar)nk+1(gkmcosmφ+hkmsinmφ)dΡnkm(cosθ)dθY=∑k=0kmax∑m=1kmsinθ(ar)nk+1(gkmsinmφ-hkmcosmφ)Ρnkm(cosθ)Ζ=-∑k=0kmax∑m=1k(nk+1)(ar)nk+1(gkmcosmφ+hkmsinmφ)Ρnkm(cosθ)}(3) 地磁学和高空物理学国际协会 (IAGA) 有一专门小组进行以5年为间隔的国际参考地磁场 (IGRF) 研究, 目前已有8代IGRF模型。对于确定了的地磁场模型称为DGRF(Definite Geomagnetic Re-ference Field) , 而且每5年改变一次, 从1945年起, 至今一共有13个确定的地磁参考场, 地磁高斯系数今后不再修改。在地磁各要素计算中, 可采用相应时间对应的DGRF模型。根据磁偏角的定义 (图2) , 其表达式可以写为: D=arctan (Y/X) (4) 根据磁倾角的定义, 其计算表达式可以写为: I=arctan(Z/X2+y2???????√)Ι=arctan(Ζ/X2+y2)(5) 在地磁分量计算时, 基于DGRF, 分别采用内插和外推的方法求