磁法反演matlab程序,磁法反演解释.ppt

磁法反演解释

磁法勘探 磁法勘探 磁法勘探 磁法勘探 三度异常：异常平面分布总体特征无明显延伸方向，呈现圆形或椭圆异常。 异常周围常被负异常包围——称等轴状异常，球体或似球形多层叠加磁性体 异常北侧出现较大负值——称似二度异常，场源走向长度有限或其长度远小于其埋深 二度、三度异常判断标准： 与磁性体走向长度(L)、水平宽度(2b)、埋藏深度(h)等因素有关。 L/h>5时，可近似作为二度体异常 (实际应用时常以异常平面等值线的半极大值围线的长、短轴之比来判断。长轴/短轴>3以上为二度异常，反之为似二度或三度异常。能较明显判别出长短轴之比划分磁性体中心剖面异常，近似当作二度体异常。) 2.双分量参量图 利用不同分量的组合就可以绘制成双分量参量图，如Za-Ha、△T-Za及△Th- △Tz等。以便进一步判断磁性体的产状等参数。下面例举Za-Ha参量图。 Za-Ha参量图的方程式，可以由磁性体磁场的直角坐标表达式Za、Ha中消去坐标变量x得到，例如斜磁化水平圆柱体Za-Ha参量图的方程式为 3.利用不同高度上磁场特征 利用磁异常断面等值线的特征 由不同高度上的剖面磁场值可绘制出XOZ断面内等值线图，又称空间等值线图。 利用不同高度上的Zamax值 利用下式迭代计算可以求得磁性体形状参数(n)，并通过查表来判断磁性体的形状： 在确定了磁性体形状之后，进一步利用其不同高度上的Zamax值及其坐标位置还可以计算其埋深和磁化特征角等参数，现以无限延深薄板和水平圆柱体为例，计算公式如下： 无限延深薄板有： 利用不同高度上Za特征值横坐标连线 由于各种磁性体Za异常的特征值的横坐标均随着剖面高度变化而发生规律性变化，因此，可利用它们连线之间的夹角和交点来判断磁性体的形状并确定其位置。 磁法勘探 由以上两式可知，Ta异常曲线均为正值，不受磁化方向和磁化特征角的影响，仅与磁性体的几何形状有关。且与板状体顶底板中心到观测点矢径，rA的rB乘积成反比关系。因此，能有效的利用极值两翼梯度变化来判断倾向。当Ta曲线对称时，则为直立有限延深板状体，当不对称时，则梯度缓的一侧即为倾向方向。并且不难证明其Ta曲线的极值横坐标位于板状体上顶中心在剖面上投影位置附近。 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 磁法勘探 最优化选择法反演的原理是将实测资料与已知形状产状的模型体所产生的理论异常对比，反复修改模型体的大小、产状等参量，使理论曲线和实测曲线较满意的符合为止。我们就取该理论曲线对应的模型体作为实际地质体的解释。 最优化选择法进行重磁异常反问题计算的步骤为： 1、给出实测异常和测点位置信息(包括地形高程数据)； 2、选择进行解释计算用的地质体模型； 3、对地质体模型各个参数(包括位置、大小、产状与物性参数)给初值； 4、计算模型体的理论异常，将该异常与实测资料进行对比； 5、评定理论异常和实测资料的符合程度，即使得它们之间的差异最小，判断是否需要修改地质体参量重新计算对比； 6、若差异不符合标准，自动修改地质体参量，使理论异常和实测资料之间的差异减小； 7、如此反复迭代下去，直到满足条件为止； 8、记下最后模型体的参数，作为解释结果。 第六节 可视化交互反演方法 磁异常反演的人机联作(又称人机对话)选择法，就是根据选定数学地质模型及磁场正演计算公式，由给定的初始模型参数计算的理论场值与实际观测场值进行对比，利用输出曲线图形对比不一致性或残差，依靠解释人员的经验和技能来修改模型及参数，重新计算理论场值进行下一次对比，反复多次直至达到满意为止。 二度半剖面人机联作的数学地质模型的基本形态选择为多边形截面水平棱柱体。 重磁人机交互反演系统(中国地质大学，武汉) 二度半重磁人机交互反演系统(中国地质大学，武汉) 三度体重磁人机交互反演系统(中国地质大学，武汉) 高精度磁测资料人机交互进精细解释 第七节 成像反演方法 磁异常是由地下磁性体在地磁场的作用下引起的，而地下介质是连续分布的，地下空间中的磁性体可视为由一系列较高或者较低磁化率块体组成的磁异常区域。根据地面观测的磁异常反演地下磁性体分布时，可将观测区域地下空间介质划分为一系列块体单元(2D反演将划分为矩形截面的无限长棱柱体)，当块体单元足够小时，假设每个块体单元内部磁化率均匀分布，允许不同块体单元具有不同的磁化率。根据位场理论，当不考虑剩磁和退磁影响时，磁化率均匀的块体单元所产生的磁异常与磁化率成正比，比例系数