对于书中的磁异常公式

对于书中的磁异常 对于书中的磁异常公式，我所编的程序，取数据: 磁化强度M =1; 球体半径r =25; 球体埋深R =50; 球体磁倾角I =90; 球体磁偏角D =0; I地磁倾角 =57.25; 0 D地磁偏角=-6; 0 网格-100?x?100, -100?x?100 网格中y向的格数Ny=201;x向的格数Nx =101; 所得正演结果如图: 采用同样的数据，用老师给的软件得到的结果如下图: 将两图相减，所得如下: 误差基本很小，可忽略。 ,,TTcosT,以上所研究的皆是?2000nT 时的磁异常计算公式:。aa T于是在此基础上，下面开始研究不小于2000 nT的情况下的磁异常计算公式:a ,,,TTT。 0 在此，取与上面磁倾角和磁偏角相关的北京地区的H、Z值。 T即已知:的三分量——正常场水平分量H=29620 nT;垂直分量Z =46180 nT。 0 ZDHDHD (由当地磁偏角=-6可算出: =Hcos，=Hsin，=Z)00x00y00 ZTHH 的三分量——、、均可由我编的程序中算出具体值。aaxaya 已知空间内两矢量线段的位置，求其夹角就好求了，根据公式知: TT,0acos,,， TT,0a 22222THHZ,，，THZ,，其中:，， aaxaya0 TTHHHHZZ,,,，,，, 0000axaxyaya 22,,，，,,,TTTTTT2cos, aa000 可求出实际的异常值。 结果如下: 将计算值与此实际值相减得: 差值见“计算值减实际值.dat”文件。由此可见:计算值与实际值之间相差了零点几的大小。Z方向上的最大值为:0.，最小值为:-0.149438 。