各地的磁倾角_各地磁偏角

用飞机携带磁力仪在空中进行的地磁测量,它比陆地磁测和海洋磁测速度快,费用省。

航空磁测可分为两种类型,一种是用磁通门磁力仪,或质子旋进磁力仪,或光泵磁力仪测量地磁场的总强度(标量);另一种是用磁通门分量磁力仪,或质子旋进分量磁力仪测量地磁场的各个分量,有的测量磁偏角、水平强度和磁倾角,有的测量地磁场的北向强度、东向强度和垂直强度。测量地磁总强度时,飞行高度较低,通常是几十米或几百米,测线也较密,线距为几百米或几千米。为了减少飞机本身产生的磁场对测量结果的影响,要把探头放在机舱外面(见彩图),用一定长度的电缆同飞机联接。测量地磁场分量时,飞行高度为几公里,线距为几十公里。测量地磁场分量的难度比测量其总强度的难度大得多,这是因为不仅要测定探头相对于参考系统的方向,而且要补偿飞机磁场对测量结果的影响。

在进行航空磁测时,除了磁力仪之外,定位和导航是很重要的辅助技术。根据测量的目的和测区的条件,可以使用不同的定位和导航方法,例如可以用地形图和照相机、惯性导航仪和多普勒导航仪,以及高精度的圆系统定位导航.航空磁测数据用模拟记录器或数字磁带来收集,并用微处理机进行处理。

尽管三分量航空磁测的精度低于地面磁测的精度,但它可以在交通不便或不可能进行地面磁测的地区进行测量,为研究这些地区的地磁场及其长期变化提供资料。另外,航空磁测还广泛地应用于金属矿普查、石油普查和地质构造研究等方面(见航空地球物理勘探)。

卫星磁测

把磁力仪放在人造卫星上进行的地磁测量。卫星磁测技术发展迅速,最初只是当卫星飞过地面接收站上空时,卫星才发射信息,目前已使用记忆装置,能获得整个卫星轨道上磁场数据。

1958年,苏联发射了世界上第一颗测量地磁场的卫星(“人造地球卫星”3号),上面装有磁通门矢量磁力仪,由于不能准确地确定仪器的方向,所以只能得到总强度的资料。以后,苏联和美国又先后发射了几颗飞行不高的测量地磁场的卫星,如美国的“先锋”3号、“宇宙”26号、“宇宙”49号、“宇宙”321号、“奥戈”2号、“奥戈”4号和“奥戈”6号,这些卫星都只携带测量地磁场总强度的磁力仪(质子旋进磁力仪或光泵磁力仪),飞行高度通常是几百公里,能够准确、迅速地测量地磁场总强度。1979年10月30日美国发射了一颗“地磁卫星”,它的轨道通过两极上空,能够覆盖整个地球表面。卫星上除装有光泵磁力仪和磁通门矢量磁力仪外,还装有星象照相机,能较准确地确定卫星的飞行姿态,因而有可能作出较准确的地磁三分量的全球测量。

通过卫星磁测,人们在很短的时间里,就能取得整个地球磁场的资料。根据卫星磁测资料,可以建立全球范围的地磁场模型,研究全球范围的磁异常,并可以研究地磁场的空间结构(见磁层)。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
圆柱体的磁异常可以通过计算圆柱体内部和外部的磁场差异得到。在Matlab中,可以使用数值计算工具箱中的函数来计算圆柱体的磁场分布和磁异常。以下是一个简单的示例代码: ```matlab % 圆柱体参数 R = 10; % 圆柱体半径 L = 50; % 圆柱体长度 chi = 0.01; % 圆柱体磁化率 inc = 30; % 磁倾角 dec = 0; % 磁偏角 % 地球磁场参数 B0 = 50000; % 地球磁场强度 inc0 = 60; % 地球磁场倾角 dec0 = 0; % 地球磁场偏角 % 计算磁场分布 [x, y, z] = meshgrid(-100:5:100, -100:5:100, -100:5:100); [Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0); % 计算磁异常 [Bx0, By0, Bz0] = cylinder_mag(R, L, 0, inc, dec, x, y, z, B0, inc0, dec0); dBx = Bx - Bx0; dBy = By - By0; dBz = Bz - Bz0; dB = sqrt(dBx.^2 + dBy.^2 + dBz.^2); % 绘制磁异常图像 figure; slice(x, y, z, dB, [], [], [-50, 0, 50]); colorbar; xlabel('X'); ylabel('Y'); zlabel('Z'); title('Cylinder Magnetic Anomaly'); ``` 该代码使用了一个名为`cylinder_mag`的自定义函数,用于计算圆柱体的磁场分布和磁异常。以下是该函数的代码: ```matlab function [Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0) % 计算圆柱体的磁场分布和磁异常 % R: 圆柱体半径 % L: 圆柱体长度 % chi: 圆柱体磁化率 % inc: 磁倾角 % dec: 磁偏角 % x, y, z: 空间坐标 % B0: 地球磁场强度 % inc0: 地球磁场倾角 % dec0: 地球磁场偏角 mu0 = 4*pi*1e-7; % 真空磁导率 % 圆柱体内部的磁场分布 r = sqrt(x.^2 + y.^2); theta = atan2(y, x); Bz1 = zeros(size(x)); Br1 = zeros(size(x)); Bx1 = Br1.*cos(theta) - Bz1.*sin(theta); By1 = Br1.*sin(theta) + Bz1.*cos(theta); Bz2 = chi*B0*L/2/mu0*(r.^2-R^2)./((r.^2+L^2/4-2*L/4*r.*cos(theta)).^(3/2)); Br2 = chi*B0*L/4/mu0*L./((r.^2+L^2/4-2*L/4*r.*cos(theta)).^(3/2)).*(L/2-r.*cos(theta)); Bx2 = Br2.*cos(theta) - Bz2.*sin(theta); By2 = Br2.*sin(theta) + Bz2.*cos(theta); Bz = Bz1 + Bz2; Bx = Bx1 + Bx2; By = By1 + By2; % 圆柱体外部的磁场分布 r = sqrt(x.^2 + y.^2); theta = atan2(y, x); phi = atan2(z, r); Bz3 = zeros(size(x)); Br3 = zeros(size(x)); Bx3 = Br3.*cos(theta) - Bz3.*sin(theta); By3 = Br3.*sin(theta) + Bz3.*cos(theta); for n = 0:20 k = n/R; Jnkr = besselj(n, k*r); Jnkp = besselj(n, k*R); Ynkp = bessely(n, k*R); Hnkr = besselh(n, 1, k*r); Hnkp = besselh(n, 1, k*R); Hnkm = besselh(n, 1, k*L/2); An = (Jnkp*Ynkp - Jnkr*Hnkp)./(Jnkp*Hnkm - Hnkp*Jnkm); Bz3 = Bz3 + k*(n+1)*An.*(chi*B0/mu0).*Jnkr.*cos(n*phi); Br3 = Br3 - k*(n+1)*An.*(chi*B0/mu0).*Jnkr.*sin(n*phi); Bx3 = Br3.*cos(theta) - Bz3.*sin(theta); By3 = Br3.*sin(theta) + Bz3.*cos(theta); end % 地球磁场 Bx0 = B0*cosd(inc0)*cosd(dec0); By0 = B0*cosd(inc0)*sind(dec0); Bz0 = B0*sind(inc0); % 总磁场 Bx = Bx + Bx3 + Bx0; By = By + By3 + By0; Bz = Bz + Bz3 + Bz0; end ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值