二磁性体磁场正演

1、2222220522222222cos cos3sin3cos sin2cos sin3sin3cos cos42sin3cos cos3cos sinaxayaxyRIARxIxyIAHmHxyRIARyIxyIAxyRZxyRIRxIARyIA1、球体磁场计算公式：、球体磁场计算公式： 式中m为球体总磁矩：3043kmM VTr22222222205/2222222222(2)sin(2)coscos4(2)cossin3sin2 cos3cossin23sin2 sinmTRxyIxyRIAxyRyxRIAxRIAxyIAyRIA2、水平圆柱体磁场计算公式：、水平圆柱体磁场计算公式： 式

2、中ms为球体主剖面内的有效磁矩：ssmM V2202222202221()sin2cos 2()1()cos2sin 2()sasssassmZRxiRxixRmHRxiRxixR 220222sinsin 2902cos 290sin2ssssmITRxiRxiixR3、有效磁化强度、有效磁化强度Ms与有效磁化倾角与有效磁化倾角is计算公式：计算公式：22 1/222211()(coscossin)(sec)sxzzsxMMMMIAIMitgtgtgIAM% 测点分布范围dx=5; % X方向测点间距dy=5; % Y方向测点间距nx=81; % X方向测点数ny=81; % Y方向测点数x

3、min=-200; % X方向起点ymin=-200; % Y方向起点x=xmin:dx:(xmin+(nx-1)*dx); % X方向范围y=ymin:dy:(ymin+(ny-1)*dy); % Y方向范围X,Y=meshgrid(x,y); % 转化为排列% 球体参数i=pi/4; %有效磁化倾角isa=pi/2; %剖面磁方位角r=20; % 球体半径 mv=4*pi*R13u=4*pi*10(-7); %磁导率M=0.7 %磁化强度 A/mm=M*v; %磁矩R=30; % 球体埋深 m% 球体Za理论磁异常Za=(u*m1*(2*R.2-(X-50).2-Y.2)*sin(i)-3

4、*R*(X-50).*cos(i)*cos(a)-3*R*Y.*cos(i)*sin(a)./(4*pi*(X-50).2+Y.2+R.2).(5/2);hax=(u*m*(2*x.2-y.2-h.2)*cos(I)*cos(a)-3*h*x.*sin(I)+3*x.*y.*cos(I)*sin(a)./(4*pi*(x.2+y.2+h.2).(5/2);hay=(u*m*(2*y.2-x.2-h.2)*cos(I)*sin(a)-3*h*y.*sin(I)+3*x.*y.*cos(I)*cos(a)./(4*pi*(x.2+y.2+h.2).(5/2);za=(u*m*(2*h.2-x.2-y.2)*sin(I)-3*h*x.*cos(I)*cos(a)-3*h*y*cos(I)*sin(a)./(4*pi*(x.2+y.2+h.2).(5/2);dt=hax.*cos(I)*cos(a)+hay.*cos(I)*sin(a)+za.*sin(I);%总磁场强度异常附图1：球体磁场Za 和T特征磁化倾角I=9O(垂直磁化)磁化倾角I=O(水平磁化)附图2：球体磁场