地磁作业(转换、延拓)-4eb5933e0912a216147929d7.pdfVIP

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