物探资料处理方法与MORPAS.pptVIP

基于MORPAS软件物化探资料处理方法 （中国地质大学数学地质遥感地质研究所） Wuhan 430074 关于重磁资料 物探资料不仅反映浅部的物性分布特征，而且也同时反映深部地质体空间展布的信息。对推测隐伏地质异常如隐伏基底（可能的矿源）、隐伏构造、隐伏岩体和隐伏矿床等都具有重要作用。 物探信息特点： 叠加场：地下所有重磁物性异常体的综合反映 多解性：不同的地质体可以有相近的地球物理场，也就是说地质体的大小、形状、深度与产状等参数的不同组合能引起相同或相近的地球物理异常。 处理目的: 实现不同性质地质体的场分解，从多角度将场源与地质体相联系，从各个侧面显示地质构造特征，使多解趋向于单一，此外还要结合地质、地球化学等其他地学信息的研究，来达到不同地质体识别的目的。 重磁资料的共性/区别 共性: 天然场 位场原理 反演的多解性 区别 重力异常:密度差引起 磁力异常:磁性差引起 处理方法 滤波与圆滑: 圆滑处理、窗口滑动平均、高低通滤波、递归滤波 位场转换： 空间域／频率域 解析延拓：上延、下延 空间导数：一阶方向导数、水平梯度；一、二阶垂向导数 磁场化极 参量转换：ΔT和ΔZ的三分量转换，伪重力异常换算 分离场： 滑动窗口法、趋势分析法、平均场法、匹配滤波法 界面反演 单界面二维、三维反演(包括Parker方法三维重力界面反演) 多界面二维、三维反演 变物性反演 局部体重磁二、三维人机联作反演 似磁化率填图 圆滑 圆滑处理方法是采用二维多项式函数来拟合局部范围内数据这一原理，利用最小二乘法确定多项式的系数，从而用多项式的函数值作为数据的圆滑值。 1）五点线圆滑 2）九点线圆滑 3）九点二次圆滑 4）25点二次圆滑 5）49点二次圆滑 滑动平均 窗口滑动平均是压制浅部和局部干扰异常最有效的方法之一，它是把窗口内数据的平均值作为其中心处的圆滑值，逐点计算。不同尺度窗口的平均结果，与一定深度上的异常体有着密切的联系。 补偿圆滑 补偿圆滑是一种傅立叶频率域低通滤波方法，其特点是能产生“圆滑”的滤波器，可压制或避免振荡效应的出现。其原理是利用一系列频率的指数函数的递归来构造滤波器进行滤波。 位场转换的方法原理 延拓 延拓的主要目的 导数 水平导数 求取水平方向导数的目的是为了突出线性构造在重力场中的反映，确定区内线性构造，结合剩余图划分出区内构造骨架。通常做四个方向的水平导数，00、450、900、1350，其结果突出于之垂直走向线性构造。 此外,还可以编制水平导数模量图，它可了解编图区内各方向构造线特征。 垂向二阶导数 不同阶次的重力位垂向导数对不同埋深物质反映不一样，右图所示为两个大小不一，中心埋深也不相同的球体。浅部小球的异常在实测中并不明显，但对该异常求二次垂向导数Wzzz后，深部大球引起的“区域异常”受压制，而浅部小球引起的“局部异常”得到放大，因此，通常把高次导数异常当作局部异常。 1、哈克公式?? 2、艾勒金斯公式（Ⅰ） 3、艾勒金斯公式（Ⅱ） 4、艾勒金斯公式（Ⅲ） 5、罗森巴赫公式 空间域位场换算 用空间域差分方法求水平导数和垂直导数。 频率域位场换算 用FFT方法进行延拓、水平导数、垂直导数、高低通滤波 磁场化极 磁化方向对磁异常特征有很大影响。斜磁化影响使异常变得十分复杂。把斜磁化的异常化为垂直磁化或顺层磁化异常，异常形态变得较为简单，便于定性分析和定量计算解释。 通常，垂直磁化下的异常垂直分量Za⊥还用来与重力异常的导数对比，判别分析引起异常的地质原因。 参量转换/磁源重力异常 计算水平分量 重磁场的分离往往是解释的一个前提，对迭加异常直接进行解释比较困难。重磁解释工作中，通常要进行局部异常与区域场分离、不同规模异常迭加的分离，以及从整体异常中提取目标地质体异常等工作。 在分离场的计算中，常用的方法有滑动窗口法、趋势分析法、平均场法、匹配滤波法。 但从本质上讲，他们的效果大同小异。适当调整他们的参数因子，可以得到几乎完全一样的结果。 趋势分析 利用多项式对数据进行拟合，从而得到数据的趋势值或原始值与趋势值之差——剩余值 平均场法 又称网格平均法，是用于区分区域与局部特征的有效方法。它是将平面数据划分成若干单元，每个单元内的数据取平均并作为单元中心处的区域场值，最后用这些区域场值内插出全区的区域场值。 匹配滤波 匹配滤波是基于不同深度的磁性体异常（傅立叶频率域）对数功率谱与其埋深之间近似关系来实现深源场与浅源场分离的方法。这种方法也常被用来区分区域场和局部场。 匹配滤波：多次分解 密度界面反演 基于Parker的傅立叶频率域正演公式设计的一种迭代反演方法。 视密度／视磁化率填图 视磁化率反演是一种物性填图方法。它是利用磁等