资源描述：
《低磁纬度区_T剖面磁异常场源深度计算的切线法_郭志宏》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第27卷第5期物探与化探Vol.27,No.52003年10月GEOPHYSICAL&GEOCHEMICALEXPLORATIONOct.,2003低磁纬度区$T剖面磁异常场源深度计算的切线法1,21,21,2郭志宏,于长春,周坚鑫(1.中国地质大学,北京100083;2.中国国土资源航空物探遥感中心,北京100083)摘要:航磁$T剖面异常场源深度计算的切线法是基于中高磁纬度区建立起来的计算剖面磁异常场源深度行之有效的传统方法,以此方法为基础,明确地提出计算低磁纬度区$T剖面负磁异常场源深度的异常/反切0概念和做法
2、,并从理论上证明了其正确性。通过$T剖面负磁异常的/反切0做法,使得传统的计算中高磁纬度区$T剖面磁异常场源深度的切线法及系数表可直接用于低磁纬度区,从而圆满地解决了低磁纬度区利用切线法计算$T剖面负磁异常场源深度的问题。关键词:低磁纬度区;航磁;$T剖面负磁异常;场源深度;切线法;/反切0中图分类号:P631.2文献标识码:A文章编号:1000-8918(2003)05-0391-04传统的航磁$T剖面异常场源深度计算的切线sinI0$T=sinA##[ZcosE+HsinE](1)法主要是针对中高磁纬度区的磁异常
3、建立起来的,sini0几十年来,该方法一直在使用并获得比较好的应用式中,A为板状体倾角;I0为地磁倾角;i0为沿剖面[1]方向的有效地磁倾角;E=90b+A-i-i效果。但对于南海等低磁纬度地区,由于接近水0为异常特平的地磁场方向,使得这些地区的$T磁异常特征征角;i为沿剖面方向的总有效磁化倾角;Z、H分别发生很大变化,出现了大量的以负磁异常为主体的为垂直板状体垂直磁化情况下的磁场垂直、水平分[2]量:异常,如果在这种情况下仍然不做任何变化地套用中高磁纬度区以$T剖面正磁异常为主体的常规x+bx-bZ=2Ms#[ar
4、ctan-arctan],(2)hh切线法,由于/所认定的正磁异常0通常不是低磁纬22(x-b)+h度区实际场源引起的异常主体部分,因而由此计算H=2Ms#ln[22]。(3)(x+b)+h出的场源深度结果及标定的场源体地面投影位置将式中:Ms为沿剖面方向的有效磁化强度;b为板状是错误的。为此,笔者明确提出计算低磁纬度区体水平半宽度;h为板状体顶面深度。$T剖面负磁异常场源深度的磁异常/反切0概念和由(1)式可以看出$T剖面磁异常的曲线形态做法。所谓/反切0磁异常,就是沿剖面寻找负磁异完全由异常特征角E决定,而系数s
5、inA#sinI0/sini0常为主体的异常极值、拐点等特征点,然后将异常剖面上下/颠倒0过来使用常规切线法计算该负值主体异常的场源深度。下面笔者将通过$T剖面磁异常场源深度计算的切线法原理及公式从理论上证明低磁纬度区$T剖面负磁异常场源深度计算的/反切0法的正确性。1$T剖面磁异常切线法计算模型$T剖面磁异常切线法计算通常采用比较有代表性的板状体模型,走向无限、延深无限的倾斜厚板[2,3]状体$T剖面磁异常的解析表达式如下:图1各种异常特征角的板状体$T异常形态示意收稿日期:2003-04-14#392#物探与化探
6、27卷只影响异常幅度。随着E角度变化,$T异常曲线ym1)/b,磁化强度系数K3=($Tmax-$Tmin)/Ms,有不同的特征,归纳为8种类型(图1):Z型、-Z从而获得倾斜磁化(E=45b~90b)板状体$T剖面型、H型、-H型、ZH型、-ZH型、-HZ型、HZ型。磁异常切线法系数表(略)。2中高磁纬度区$T剖面磁异常切线法2.1$T剖面磁异常切线法正演计算在不考虑剩磁影响的情况下,中高磁纬度区板状体的$T剖面磁异常形态主要对应于图1中的垂直磁化(E=0b)Z型和倾斜磁化ZH型、-HZ型等以正磁异常为主体的异常形
7、态,而H型、-H型为2种向低磁纬度区过渡的正负等量伴生的$T剖面磁异常形态。由前面(1)~(3)式可正演计算得到板状[2]n图3倾斜磁化(E=45b~90b)$T异常曲线切线法示意体$T剖面磁异常切线法的系数表。2.2$T剖面磁异常切线法反演计算2.1.1垂直磁化(E=0b)$T剖面磁异常切线法系$T剖面磁异常切线法反演计算是根据前面正数表演计算获得的切线法系数表来进行的。首先通过如图2所示,$Tmax、$TP为异常曲线极大值和$T磁异常极大值处切线分别到两翼极小值处切线拐点处的异常值,yP、ym、y0为异常曲线拐点
8、坐标、之间垂向距离S1、S2(见图3)的比值来判断异常曲拐点处切线与极大值点处切线及零值线交点坐标。线形态,然后根据不同的异常曲线形态选择不同的令E=0b,通过B=b/h的渐变,由(1)~(3)式可计切线法系数表来进行异常场源顶面深度、宽度和磁算得到异常形态参数(y0-ym)/ym,深度系数K1=OP化强度的反演计算。(y0-ym)/h、K2=
