磁法勘探_06_磁异常的处理与换算分解
磁异常处理与转换

a0 ∞ T ( x) = + ∑ (a k cos ku 0 x + bk sin ku 0 x ) 2 K =1
T(x)
2L
三、频率域转换的思路
2、利用欧(尤)拉公式将富氏级数的三角形式 、利用欧( 变换成富氏级数的复数形式,并用富氏积分表示, 变换成富氏级数的复数形式,并用富氏积分表谱
(u,v)﹒ Φ(u,v)﹒ ΔTa﹋(u,v,0) → ΔTb﹋(u,v,0) 其中Φ为转换因子(响应函数、滤波因子),变换实际上就是 其中Φ为转换因子(响应函数、滤波因子),变换实际上就是 ), Φ的累积
第三步
通过反变换由转换后的频谱变为空间域, 通过反变换由转换后的频谱变为空间域,获得转换后 的磁异常, 的磁异常,从而达到转换的目的
S T (u ) = ∫ T ( x)e −iux dx
−∞ ∞
正变换 逆变换
1 T ( x) = 2π
∫
∞
−∞
S T (u )e −iux du
三、频率域转换的思路
3、ST(u)称为 (x)的频谱。通过正变换可计 、 )称为T( )的频谱。 算ST(u),反之可求T(x) ),反之可求 ( ) ),反之可求
sincos2l2lttxx三频率域转换的思路2利用欧尤拉公式将富氏级数的三角形式变换成富氏级数的复数形式并用富氏积分表示即付氏变换的正变换和逆变换正变换正变换逆变换逆变换三频率域转换的思路正变换正变换逆变换逆变换四实施步骤第一步通过正变换计算原始异常的频谱taxy0或tax0第二步计算转换后的频谱其中为转换因子响应函数滤波因子变换实际上就是的累积第三步通过反变换由转换后的频谱变为空间域获得转换后的磁异常从而达到转换的目的1褶积变换为乘积使计算简单2转换的实质就是变成了一种响应函数的累积因此可以连续实现几种转换3在频率域中实现转换因根据换算的目的及要求设计或选定相应的频率响应函数滤波因子频率域计算速度快简单各种转换仅是相应因子的乘积过程无边界损失延拓高度无局限存在一些不良效应空间域计算速度慢复杂由于边界损失会使转换后的面积减小转换受点距控制无边界损失
磁异常的处理、解释及应用

上面,简单概述了磁异常各类正问题及其正演方法。其中,均匀磁化规则形体正(演) 问题、正演方法及场的解析表达式,是磁法的基础,具有重要的理论意义和实际意义将重点 讨论。
第二节 有效磁化强度矢量与总磁场异常 Δ T 的一般表达
一、有效磁化强度矢量
已知总磁化强度矢量由感应磁化强度与剩余磁化强度两矢量组成。设总磁化强度矢量 M
K 这种磁性体的参数 k 和 M r 需用张量来描述,其正演问题是磁法中最复杂的正问题。从
70 年代后期,国内外学者相继研究出一些数值正演方法。我国学者把有限元和边界元等数 值计算方法引用到这一复杂正演问题中来,取得了一系列有理论和实际价值的成果。
(五)磁场的模拟测定
前述各类正问题的求解还可以通过实验室模拟测定来解决。模拟测定方法分为静磁场 模拟方法与低频交变场模拟方法。实践已经证明,两类模拟测定方法是可行的。
3、组合体近似法:把磁化强度均匀或分区均匀的任意形态磁性体,用多个均匀磁化规 则形体的组合形体近似代替;各个均匀磁化规则形体的磁化强度可以相同或不同。该磁性体 磁场的近似值,等于各规则体解析场值之和。作为组件的规则形体有正方体、直立长方体、 倾斜长方体、有限长水平n棱柱体等。因为直立长方体的多个ln项可以合并成一项计算,而 且在一定条件下多个tg-1项亦可合并计算,使计算速度大大加快,又因其组合任意形体的能 力较强,故直立长方体组合法得到了普遍应用。
磁法勘探-磁异常的数据处理

第九章磁异常的数据处理前面正演计算的假设条件:形状规则、均匀磁化、观测面水平、单个异常体………在此条件下建立磁体与异常特征之间的关系作为解释理论。
这与实际情况有很大的偏差,需要对观测数据进行处理实际:①剩余磁化强度;②地形起伏不平;③测量偶然误差;④地表干扰磁场;⑤多个磁性体。
一、磁异常处理和转换的目的1、使实际异常满足或接近解释理论所要求的假设条件例:曲面的观测数据→水平面上的数据分解叠加异常→独立异常2、使实际异常满足解释方法的要求①某分量→其它分量②磁场值→频谱提供多方面异常信息3、突出异常某一方面的特点如:上延压制浅部磁性体的异常方向导数突出某一走向方向的磁异常特征磁异常处理的主要内容:① 磁异常的圆滑滤波和插值;区域与局部场,深部场与浅源场 ② 异常空间换算:场的空间解析延拓;③ 磁异常参数之间的换算;由实测异常进行,T ,a Z ,a H ,a T 之间互算④ 磁异常的导数换算;计算水平、垂向导数⑤ 不同磁化方向的磁异常换算;⑥ 频率中磁异常的各种换算和数字滤波;磁异常的地形影响校正等内容。
方法分为:空间域频率域:频率域速度快,方法简单,现已成为主要方法。
各种处理方法尤其不同的物理原理和数学方法,处理的目的也不同。
对某一地区而言,并非一定要进行所有的数据处理方法,而应根据具体情况和异常特点,合理的选择,进行恰当的处理,这跟磁异常的解释效果有很大关系。
二、磁异常的圆滑滤波和插值计算1、主要作用这种数据处理的主要作用是消除磁测过程中的随机误差,地表附近的随机干扰以及磁化不均匀的影响。
这些影响在磁异常曲线上表现为无规律的高频跳动,影响了主体异常。
所谓高频干扰,是把磁异常曲线类比为电学或波动学上的震动曲线,随机干扰的频率比较高,起伏不规则,这些服从正态分布规律,起伏平均值为零。
特别是这些干扰在进行场的相似解析延拓和导数换算时,还会得到放大,使磁异常发生更大畸变。
2、目前常用的方法①最小二乘圆滑方法是一个函数的拟合问题,用一个拟和函数(一般常用多项式)去拟合离散的实测异常值,是多项式与实测异常的偏差平方和最小,以达到光滑异常曲线的目的。
第二章 磁法勘探 第一节(一) 地球的磁场及磁异常解读

磁法勘探可用于地质调查的各个阶段。在地
质填图时﹐磁法勘探可以划分沉积岩﹑喷出岩﹑ 基性岩﹑超基性岩及变质岩的分布范围﹔可以研 究沉积岩下面的基底构造﹔查明各种控制成矿的 构造﹐如深大断裂和火山口等。在普查找矿时﹐ 磁法勘探可用来直接寻找磁铁矿床﹐并可与其他 物探方法配合﹐间接寻找或预测石油﹑天然气 ﹑
向上。
由地磁场的基本特征,如地球有两个磁极,
磁极处的地磁场约等于磁赤道上的地磁场的两倍 及地磁场的等强度线,等倾线大致与纬度线平行 等,说明地磁场与一个磁偶极子的磁场相近。确 切地说,现代地磁场与一个磁心位于地心、磁轴 与地理轴夹角为11.5°、磁矩约等于7.9×1022 A.m2的磁偶极子的磁场拟合的最佳。通常称这个 磁偶极子为地心偶极子。
2.地磁要素
地磁场总强度 T 是矢量,为描述地磁场总强
度T 在地表某一点的状态,我们定义若干个地磁
要素。将空间直角坐标系的原点置于考察点,x轴
指地理北(或真北)N,z轴铅直向下。
正,图南中半,球I为T 地上磁倾倾,角规,定北I为半负球;ZT 为下地倾磁,场规垂定直I为
分水理分量平北量,分向,北 量 东 全半 , 偏 球D球 全 皆为Z球 指为正皆向正,指真,西磁北南偏北,半DY;为为球D负地Z为为;磁地负X场磁;为东偏H地向角为磁分,地场量H磁北,自场H向地 东偏Y为正,H 西偏Y为负。 以上七个量称为地磁 要素,它们的关系如下:
煤﹑铜﹑铝﹑镍和其他金属﹑金刚石等。在勘探 磁铁矿床时﹐结合钻探资料﹐可以推定矿体的形 状﹐指导正确布置钻孔和寻找钻孔旁侧及深部的 盲矿体。此外﹐磁法勘探还可用于研究深部地质 构造和解决其他地质问题﹐以及应用于考古学等 方面。
钒钛磁铁矿
轻型飞机航空磁力/磁梯度测量
磁法勘探06磁异常的处理与换算资料

第一节 磁异常的处理与换算的目的意义
应当指出,磁异常处理和转换时,有两个问题必须要明确: 1.应当合理的选择处理和转换的方法。由于转换、处理方法 较多,具有各自的特点、作用、适用条件,不应盲目的对各 种方法都使用一遍。应当认真分析磁异常特征、测区内地质、 物性情况及所要解决的地质问题,根据各个方法的功能和适 用条件,合理的选择若干种处理方法; 2.磁异常的处理和转换只是一种数学加工处理,它能使资料 中某些信息更加突出和明显。但不能获得在观测数据中不包 含的信息。数学变换只能改变异常的信噪比,而不能提供新 信息;因此,在应用各个方法时必须要注意到实际资料的精
15
第二节 磁异常的处理
1.剖面网格化
16
17
第三节 磁异常的空间转换
延拓是把原观测面的磁异常通过一定的数学方法换算到高 于或低于原观测面上,分为向上延拓与向下延拓。向上延拓 是一种常用的处理方法,它的主要用途是削弱局部干扰异 常,反映深部异常。我们知道,磁场随距离的衰减速度与 具磁性的地质体体积有关。体积大,磁场衰减慢;体积小, 磁场衰减快。对于同样大小的地质体,磁场随距离衰减的 速度与地质体埋深有关。埋深大,磁场衰减慢;埋深小, 磁场衰减快。因此小而浅的地质体磁场比大而深的地质体 磁场随距离衰减要快得多。这样就可以通过向上延拓来压 制局部异常的干扰,反映出深部大的地质体。
是很重要的。随着磁测量精度的不断提高,实测异常中所包含 的可靠信息也不断增加。如何有效地提取和利用这些信息,就 成为磁异常解释理论研究的重要课题。早在20世纪40、50年代, 诸如导数异常的计算,磁场解析延拓,化磁极等处理方法已相 继问世。到60、70年代,由于电子计算机的广泛应用,使磁异 常的处理和转换容易实现,从而其理论和方法得到了迅速的发 展,并不断得到完善。由于在实践中磁异常的转换和处理对提 高磁方法解决问题的能力和改善地质效果起到了应有的作用, 因此它已成为当今磁异常推断解释中不可缺少的重要环节。
磁法数据处理、异常反演与解释的常用方法及常见问题探讨

异 常推 断解释变 得较复杂 , 在实 际应 用中 , 产生不少使 人 困惑 的 问题 。结合 笔 者 实 际 工 作 经验 , 其进 行 初 步 探 讨 , 同行 参 对 供
考。
模 型, 每个模型的几何形状要尽可能简单 。在计算 时, 关键 问题
是 给 定 模 型 的初 值 , 为 并 非 任 意 给 一 组 初 值 , 算 机 都 会 自 因 计 动 选 择 出接 近 实 际地 质 体 的模 型 。所 以要 根 据 实 际 的 地 质 、 钻 井 及 其 它 地 球 物 理 资料 给 出接 近 实 际 地 质 体模 型参 数 的初 值 。 () 用 合 适 的 计算 方法 3应
反演计 算过程 能否快速 、 收敛并得 到精确 的结果 , 了给 除
出 适 当 的模 型体 参 数 初 值 外 , 要 取 决 于 所 用 的计 算 方 法 。应 主 针 对 不 同 的 问题 选 择 适 当 的计 算 方 法 , 且 在 计 算 过 程 的 不 同 而 阶段 不 断调 整 、 改 所 用 的计 算方 法 。 修 () 加 压 制 多 解 性 ( 唯 一 性) 资 料 和 技术 4增 非 的
1 磁 法数 据 预 处 理 常用 方 法
对 实 测 数 据 进 行 日变 、 点 、 常 场 等 改 正 后 , 注 意 消 除 基 正 应 异 常 数 据 的 误 差 与 干扰 。误 差主 要 源 于 仪 器 的状 态和 操 作 及 点 位 误 差 ; 扰 主 要 是 指 人 文 或 地 质 因 素 的干 扰 。 在严 格 执 行 技 干 术 规 范 下 , 有 人 文 干 扰 的 数 据 一 般 不 作 为 成 果 ; 质 干 扰 通 含 地 常 指 与 勘 探 目标 无 关 的地 质 因素 引起 的异 常 。浅 表局 部 的地 质
磁各向异性条件下的磁法勘探正问题及其解法

磁各向异性条件下的磁法勘探正问题及其解法
磁法勘探在今天依然是不可替代的地质勘探方法,由于重力、地磁场的影响,磁法在磁各
向异性条件下也常常会出现求解正确结果的问题。
但当今已非常成熟的计算机技术为解决该问题提供了可能。
首先,磁法在磁各向异性条件下勘探所要解决的问题是可以推求出磁异常的极大值及其极值点。
由于地磁场的大小和方向都随时间而变化,所以磁法勘探的结果也随之变化。
而地磁磁性型介质的磁异常特征是各向异性的,所以需要考虑各向异性的因素去估算出极大值
和极值位置,这样也就需要地质师工作量大。
为了解决磁各向异性条件下勘探的问题,一般采用能有效解决磁各向异性条件下磁勘探问题的奇异矩阵反演算法。
该算法是将被测区域的磁异常测量数据映射到一个正规地质格架中,再在这个相互参照的网格中进行扩散算子的处理过程,能够根据测量数据求取地质倾
角和应力参数等数据,得到更准确的结果。
于此同时,在实施这种反演算法的一般步骤中,还会考虑到地理变量的影响因素,如真实地质起伏,岩石组合等,以便获得更准确的测量结果。
从上述可以看出,磁法在磁各向异性条件下的勘探问题可以通过运用奇异矩阵反演算法来
解决,该算法能够有效反映磁各向异性条件,解决磁勘探正确对抗磁各向异性条件而导致
的磁异常测量偏移,有效得到磁异常及其极值的精确位置,有效辅助精确测量地质倾角和应力参数等。
总之,磁法在磁各向异性条件下的勘探正确求解问题可以采用奇异矩阵反演算法,它能有效考虑到磁各向异性条件,从而获得更准确的测量数据。
磁异常的处理、解释及应用

K
K
K
M s 引起,故称 M s 为有效磁化强度。 M H 与ox轴夹角为A′,磁性体走向与 M H 的夹角为A
KK
K
K
(y方向为走向), M H 与 M 夹角为I,设 M 的方向余弦为(α ,β, γ ); M s 的方向余弦为
K (α s , γ s ), M s 与ox轴夹角为is;则有:
M = (M α = Mx
2、三角形面多面体近似法:这是与多边形面多面体近似法类似的近似正演方法。该方 法正演时,首先,把任意形态磁性体外表面用多个三角形平面构成的封闭面代替;其次,由 已知的磁化强度算出每个三角形面的磁荷面密度;然后,采用高斯求积公式,对每一个三角 形面的磁场作数值计算,再迭加起来。为了应用三角形的高斯求积公式对每个三角面进行数 值积分,采取了一种将任意三角面上的坐标变为二维坐标的方法。
1、直接积分法:由磁性体磁位和磁场的积分表达式出发,在确定了积分上、下限之后, 直接通过积分运算,求得磁位和磁场的解析表达式。为了减少重复性的积分运算,根据位场 迭加原理,人们找到了这样一种积分求解途径:先求得简单形体磁位、磁场表达式,再由其
71
表达式出发进一步积分,求解更复杂形体磁位、磁场表达式。例如:由球体→水平圆柱体→ 薄板状体→厚板状体,等等。
Ta=T-T0 而△T是T与T0的模量差,即:
(一)均匀磁化规则磁性体磁场的正演方法
研究磁力勘探正问题的初期,人们首先致力于求解最为简单的均匀磁化规则形体的正 问题。这些规则形体有球体、水平圆柱体、板状体、长方体、断层、对称背斜等。正演求解 时,假定磁化强度为常向量,即体内各点磁化强度大小相等,方向相同。
磁化均匀和形态规则的假设,使磁性体的正问题大为简化,并给出了解析表达式。求 得它的解析表达式的方法有如下几种:
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
18
图1是内蒙某地用磁力勘探普查超基性岩的实例。该地区 浅部盖有一层不厚的玄武岩,使磁场表现为强烈的跳动。 为压制玄武岩的干扰,将磁场向上延拓了500m。由图可 知,向上延拓的磁场压制了玄武岩的干扰。同时右侧部分 反映了深部的超基性岩磁场。
图1 用向上延拓压制浅部玄武岩异常的影响 1.玄武岩 2.沉积岩
19
图2用向上延拓压制了浅部矿体的异常,突出了深部盲矿体产生 的低缓异常。
20
图2 用向上延拓压制浅部矿体的异常
通过向上延拓来研究深部磁性基底构造也是其应用的一个重要方 面。如把一个地区航磁资料先化极向上延拓20km,消除了浅部磁 性体影响后再作磁场分区。
第一节 磁异常的处理与换算的目的意义
应当指出,磁异常处理和转换时,有两个问题必须要明确: 1.应当合理的选择处理和转换的方法。由于转换、处理方法 较多,具有各自的特点、作用、适用条件,不应盲目的对各 种方法都使用一遍。应当认真分析磁异常特征、测区内地质、 物性情况及所要解决的地质问题,根据各个方法的功能和适 用条件,合理的选择若干种处理方法; 2.磁异常的处理和转换只是一种数学加工处理,它能使资料 中某些信息更加突出和明显。但不能获得在观测数据中不包 含的信息。数学变换只能改变异常的信噪比,而不能提供新 信息;因此,在应用各个方法时必须要注意到实际资料的精 度和处理方法本滑法
是依据异常剖面上的变化应具有一定的连续、渐变的规律,徒手修改(平 滑)掉某些明显的突变点。这种做法的要求是: ①平滑前后各相应点的异常值的偏差不应超过实测异常的均方误差; ②尽可能使平滑前后剖面曲线所围成的面积相等,重心不变。 (2)最小二乘平滑法 尽管偶然误差会使异常曲线不光滑而成锯齿状,但并不会改变异常曲 线变化的基本趋势。可以用线性公式或者多项式来拟合这种变化趋势
7
① 线性平滑公式
8
9
10
11
2.平面异常平滑法
12
平滑中包括线性平滑和二次曲线平滑
13
第二节 磁异常的处理
14
15
◆区域异常 分布较广的中深部地质因素引 起的磁力异常,其特征是异常幅值较 大,异常范围也较大,但异常梯度小。 ◆局部异常 相对区域因素而言,范围有限 的研究对象(如构造矿产)引起的范 围和幅值较小的异常,但异常梯度相 对较大。局部异常也称剩余异常。
2
第一节 磁异常的处理与换算的目的意义
3
磁异常处理和转换意义有: 实践表明,磁异常的处理和转换对于提高推断解释的效果 是很重要的。随着磁测量精度的不断提高,实测异常中所包含 的可靠信息也不断增加。如何有效地提取和利用这些信息,就 成为磁异常解释理论研究的重要课题。早在20世纪40、50年代, 诸如导数异常的计算,磁场解析延拓,化磁极等处理方法已相 继问世。到60、70年代,由于电子计算机的广泛应用,使磁异 常的处理和转换容易实现,从而其理论和方法得到了迅速的发 展,并不断得到完善。由于在实践中磁异常的转换和处理对提 高磁方法解决问题的能力和改善地质效果起到了应有的作用, 因此它已成为当今磁异常推断解释中不可缺少的重要环节。
5
第二节 磁异常的处理
磁异常的平滑 (消除干扰、分离叠加异常) 磁异常的插值(补缺、数据网格化、分离叠加 异常) 磁异常的数据网格化(将不规则测网上测的磁 异常化为规则测网磁异常)。
6
第二节 磁异常的处理
野外实测异常中总包含有测量的偶然误差和近地表不均匀体产生的干扰, 使实测异常表现出不规则的起伏。所以对异常处理时要先进行平滑,消除 干扰,突出主要异常体。下面分剖面和平面异常两种情况介绍圆滑方法。
第二节 磁异常的处理
1.剖面网格化
16
17
第三节 磁异常的空间转换
延拓是把原观测面的磁异常通过一定的数学方法换算到高 于或低于原观测面上,分为向上延拓与向下延拓。向上延拓 是一种常用的处理方法,它的主要用途是削弱局部干扰异 常,反映深部异常。我们知道,磁场随距离的衰减速度与 具磁性的地质体体积有关。体积大,磁场衰减慢;体积小, 磁场衰减快。对于同样大小的地质体,磁场随距离衰减的 速度与地质体埋深有关。埋深大,磁场衰减慢;埋深小, 磁场衰减快。因此小而浅的地质体磁场比大而深的地质体 磁场随距离衰减要快得多。这样就可以通过向上延拓来压 制局部异常的干扰,反映出深部大的地质体。
第一节 磁异常的处理与换算的目的意义
目前磁异常的转换处理主要有: 1.圆滑、划分异常(如区域场和局部场的分离,深源场与浅源场 的分离等); 2.磁异常的空间转换(由实测异常换算其它无源空间部分的磁场, 也称解析延拓); 3.磁异常换算(由实测异常计算垂向导数、水平方向导数等)、 不同磁化方向之间的换算(如化磁极等)以及曲面上磁异常转换 等。分量换算(由实测异常进行ΔT、Za、Ha及Ta之间的分量换 算)、 4.波数域处理和换算 磁异常转换处理的方法包括空间域和频率域 两类。频率域方法由于速度快,方法简单等优点,已成为磁异常 4 转换处理的主要方法。
第七章 磁异常的处理与换算
概述:
目的
磁异常的处理与转换是磁解释理论的一个 重要组成部分。尤其在磁勘探已从孤立异常解 释转向大面积资料(叠加异常)解释的今天,为了 突出某些信息,压制另一些信息,磁异常处理 与转换的意义变得更为重要。
1
第一节 磁异常的处理与换算的目的意义
磁异常处理和转换的目的有: ① 使实际异常满足或接近解释理论所要求的假设条件。例如把 分布在曲面上的实测异常换算成分布在同一平面上的异常;把叠 加异常分解为孤立异常;或把似二度异常转换为二度异常等。即 把复杂异常处理成简单异常,以便于解释。 ② 使实际异常满足解释方法的要求。例如由磁场某单分量测量 结果换算其它分量的值;或者由磁场值转换成为频谱值等。从而 可以提供多方面的异常信息来满足一些解释方法本身的要求。 ③ 突出重、磁异常某一方面的特点。例如通过向上延拓等方法 来压制浅部密度体、磁性体的异常,相对突出深部密度体、磁性 体的异常,通过走向滤波或换算方向导数来相对突出某一走向方 向的重、磁异常特征等。