中国地质大学（北京） 程名称课：应用地磁学 班号：1 0 1 0 5 2 ———————————————————————————————————————————指导老师： 李淑玲 日期：2 0 0 8 年 5 月 2 7 日 应用地磁学实验报告 实验2 ——磁异常转换计算 学号： 10105218 姓名： 朱占升 一、实验目的 1、掌握水平圆柱体磁场异常分布； 2、用Matlab 实现水平圆柱体的磁异常场正演计算； 3、利用正演结果进行磁异常分量之间的换算； 4、通过程序换算认知测点间距即采样点数对换算效果的影响； 5、加深对磁法勘探的理解认识； 二、程序代码 %磁法异常换算 %剖面为北向 A’=0 度，则有 I=is。 %所测数据均在同一水平面，柱体深 30m，半径 8m,测点数为 100, %测点间距依次选取 2 4 6 8m clc clear fprintf('\n 柱体深 30m，半径 8m,测点数为 100,点距依次选取 2 4 6 8m\n'); for b=2:2:8 fprintf('\n 点距取%dm\n',b); figure('color','w','NumberTitle','off','name','za-→ha'); x2=1:b:100*b; %点距为 bm x1=x2-100*b/2; h=30; %柱深 R=8; %柱体半径 m s=pi*(R^2);%柱体截面积 k=0.2; %磁铁矿磁化率 u=4*pi*10^(-7); %磁导率 B=50000; %nT 磁感应强度 H=B/u ; %磁化场强度 M=k*H; %磁化强度 m=M*s; %磁矩 a=0; %剖面为北向 A’=0 度. I=90/180*pi; %倾斜角 is=atan(tan(I)*csc(pi/2-a)); hold on za=u*m*((h.^2-x1.^2)*sin(is)-2*h*x1.*cos(is))./(2*pi*(x1.^2+h.^2).^2); hax=-1*u*m*((h.^2-x1.^2)*cos(is)+2*h*x1.*sin(is))./(2*pi*(x1.^2+h.^2).^2); plot(x2, za,'.-m'); plot(x2, hax,'.-g'); title('za 转换为 hax') xlabel(' X 剖面走向/m'); ylabel('磁异常 nT'); c=[0.4268 0.1749 0.1103 0.0813 0.0645 0.0536 0.0458 0.0400 0.0355 0.1759]; n=length(c); %转换系数个数 m=length(za); for i=(n+1):(m-n) haxz=0; for j=1:n haxz=haxz+c(j)*(za(i+j)-za(i-j)); end haxzh(i)=haxz; end for i=1:80 haxzh1(i)=haxzh(i+10); end x=((n+1)*b):b:(100*b-n*b); plot(x,haxzh1,'.-') legend('za','hax','za→hax'); end 中国地质大学（北京） 程名称课：应用地磁学 班号：1 0 1 0 5 2 ———————————————————————————————————————————指导老...
