重磁异常的处理与转换

1、 目的：复杂简化（曲面平面；迭加孤立； 等） 另一分量； 突出重磁异常某一方面特点压制浅部、突出深度可突出浅或深，或某一方面）满足解释方法某分量 重磁场值频谱值。 （上延 滤波 主要内容1、圆滑、划分异常、插值、数据网格化2、低、高、带通滤波3、分量换算（以磁为主）4、上、下延拓5、导数换算6、（对磁）不同磁化方向换算7、求磁源重力异常8 、视磁化率填图02004006008001000020040060080010000200400600800100002004006008001000020040060080010000200400600800100002004006008001000020

2、04006008001000上延压制浅部矿体异常0200400600800100002004006008001000020040060080010000200400600800100002004006008001000020040060080010000200400600800100002004006008001000地磁场02004006008001000I=30度0200400600800100002004006008001000I=60度0200400600800100002004006008001000I= - 30度0200400600800100002004006008001000

3、I= - 60度0200400600800100002004006008001000I=90度02004006008001000原则 选择合适的处理转换方法（要求使 用者必须掌握各种处理与转换方法的原理）。 处理转换只是一种数学加工方法，能突出某些信息，而不能获得数据中不包含的信息。 方法1、空间域2、频率域（波数域）一、磁异常频谱(一)、点极rqUm04100tUT. . . 2202242vueqUvuhm2222000022)(24),(vuevunvmuliqvuTvuhm如果点极不在原点下方，乘上 )(vuie(二)、磁偶极子磁异常谱(三)、直立矩形柱体000( , ), ,( ,

4、; , ; ,)( , ; , )( , ; , )T u vA H u v hM u v l m n l m nS u v a bD u v 22);,(vuhehvuH二、重磁异常频谱特点)()(),(2200022vunvmulivunmvluihvuM)2sin(2)2sin(2),(bvvauuhvuS)(),(vuiehvuD深度因子磁化因子水平尺度因子位移因子二、重磁异常频谱特点 特点：重磁异常频谱都是一些因子的乘积，即，深度因子、方向因子、水平尺度因子、位移因子连乘积000( , ), ,( , ; , ; ,)( , ; , )( , ; , )T u vA H u v hM

5、 u v l m n l m nS u v a bD u v 空间域的褶积运算频率域的乘积运算具有快速和简单的特点。寻求)()0 ,()()0 ,(),(112xxTdxTyxT)()0 ,(),(12uuTzuTFTFTFT三、频率谱处理的核心运算降级)(u向上延拓： (, 0)P xx(, 0)Uxz( , )Q rv向上延拓： (, 0)P xx(, 0)Uxz( , )Q rv向上延拓：20(0)( , )0UzUU x oz (, 0)P xx(, 0)Uxz( , )Q rv)(22hxhaaHZ )(22hxx)(x)(u2222(0)ifxfhxedxi eihxhfhfxie

6、dxehxh2222fhe2i已知：);,(,111hvuHvuT）（ 未知：);,(,222hvuHvuT）（其中：221);,(11vuhehvuH);,(22222vuhehvuH二维延拓：将h1高度异常延拓至h2高度2212)(122121vuhheTTHHTT22),(),(12vuhevuTvuT延拓过程：(1) (2) (3) ),(,),(11vuTTFyxT),(),(2122vuTevuTvuh),(,),(212yxTTFvuT两个异常在频率域的关系 二维延拓：将h1高度异常延拓至h2高度垂向导数：1222TvnT12TiuT12TivTx方向n阶水平导数：112)(TiuxTTnnn二维导数x方向水平导数：y方向水平导数：112)(TivyTTnnny方向n阶水平导数： 任意方向的导数：STxTSxyTSyzTSzxTcosyTcoszTcos2TST1xT1coscoscosyT1zT12221( coscos)cosi uvuvTsxyz0、45、90、135方向的水平导数彩色图0、45、90、135的方向导数灰度图22220002200022000()()()aaxayuvZTi u lv mnuviuHTi u lv mnuvivHTi u lv mnuv 其中：00cos;l0000cos; cosmn分量转换：由T计算Za,H