综合利用航磁三轴梯度计算全张量梯度.PDF

第61 卷 增刊 地地 质质 论论 评评 GEOLOGICALGEOLOGICAL REVIEWREVIEW Vol. 61 Supp. 2015 年 4 月 Apr., 2015 综合利用航磁三轴梯度计算全张量梯度* 骆遥，王平，段树岭，刘浩军，王金龙，安战锋 中国国土资源航空物探遥感中心，北京，1000831 航空磁测广泛应用于区域地球物理调查及相 1 全张量梯度计算 关领域，通常意义上仅测量地磁场强度（Total Magnetic Intensity, TMI ）。磁梯度测量基本不受地磁 航磁异常可近似为磁异常矢量在主磁场方向 场随时间变化以和区域地磁场影响，具有较高的分 的投影，对其求三个正交方向上的偏导数，可得到 辨率。随着超导量子干涉器件（Superconducting 航磁异常的三个轴向梯度。由于正常地磁场的梯度 Quantum Interference Devices，SQUID ）的发展，对 很小且近似不变，相对于磁异常的梯度则可以忽略 地磁场矢量的梯度测量成为可能。进入 21 世纪初 不计，因此航磁三轴梯度可近似为磁异常梯度即全 SQUID 磁梯度仪研究取得重大进展，德国IPHP 成 张量梯度在主磁场方向上的投影。将全张量梯度用 功利用SQUID 梯度仪进行了航磁全张量梯度测量， 磁位的偏导数表示，可利用傅里叶变换将航磁梯度 并获取了各梯度分量的实测资料（Stolz et al., 写成与磁位相关的线性表达，并通过最小二乘确定 2006 ）。磁张量梯度测量不仅对地球物理勘探意义 磁位。最终通过磁位计算全张量磁梯度，实现将航 重大，在未爆军火探测、磁反潜、地磁导航等相关 磁三轴梯度转换全张量梯度。 领域也具有广阔的应用空间。美国等西方国家也致 2 理论模型计算 力于 SQUID 梯度测量系统研制，但航磁全张量梯 度测量难度极大，除 IPHT 还未见其他航磁全张量 为了检验上述转换方法，选取长方体磁性模型 梯度实测资料。 进行数值检验。该模型为长方体，长宽高分别是 8 我国已成功研制出具有完全自主知识产权的 km 、4 km 和4 km，中心埋深为3 km，磁化强度1A/m 航磁全轴梯度勘查系统，并多次进行航磁梯度测量 方向为偏角30°倾角50° 。模型正演计算的网格距为 （骆遥等, 2011 ），具备实际应用能力。由于航磁总 100 m 高度为0 m 处的梯度数据，其中主磁场磁倾 场梯度可近似为全张量梯度主磁场方向的投影，一 角60°、磁偏角10°。将正演的ΔT 三方向梯度数据 定条件下二者存在互换关系，利用实测航磁总场梯 进行转换，图 1 给出了转换得到的全张量梯度，由 度计算磁张量梯度则成为目前条件下我国取得实 于全张量梯度具有对称性，图1 仅给出全张量梯度 际磁梯度张资料的一种途径。Nelson （1988）曾提 的6 个量。较磁异常标量梯度相比，图1 反映了磁 出利用总场的垂向梯度或水平梯度计算磁梯度张 异常矢量的空间变化率信息，其地质细节更为丰 量，但仅利用了垂向梯度，没有挖掘磁总场在三维 富。将图1 同理论模型比较，二者无论是在形态还 空间上的变化率信息。针对航磁三轴梯度测量的广 是数值上均非常一致，其中绝度误差的最大值为 1 泛应用，提出航磁总场三轴梯度转换全张量梯度的 pT/m ，绝度误差的最大标准差为0.2 pT/m，该误差 计算方法，并将该方法用于实测三轴梯度，首次获 表明转换方法的数值精度，其能满足实际资料处理 得了我国的磁张量梯度数据。 要求。 注：本文为国家高技术研究发