关于布格重力异常计算及资料处理与反演和解释的报告.docVIP

关于 布格重力异常计算及资料处理与反演和解释 姓名：林俊 班级：061084-27 学号：2081003195 指导老师：陈超 日期：2011.4.14 目录 前言…………………………………………2 目的……………………………………………2 任务要求………………………………………2 工作过程………………………………………2 成果……………………………………………2 工作内容及步骤………………………………3 §1-布格重力异常计算…………………………3 §2-布格重力异常处理…………………………3 1.绘制平面等值线图………………………………3 2.异常处理（分离区域异常和局部异常）…………6 §3-布格重力异常反演——特征点法反演…………11 §4-布格重力异常的解释…………………………13 评述与结论………………………………………13 评述………………………………………………13 结论………………………………………………14 关于布格重力异常计算及资料处理与反演和解释 前言 目的：熟悉并掌握布格重力异常计算及资料处理与反演和解释 任务要求： 根据在一个地区重力测量的结果，计算出布格重力异常，并根据异常进行资料处理和解释，并完成一份工作报告。 工作过程： (1)利用实测的相对重力值、相对高程值和X,Y坐标值，计算各种校正（地形校正除外），纬度校正用 计算，自由空间（或高度）校正用 计算，中间层校正用 计算，已知地表物质密度为2.50g/cm3，起算点纬度为45°； (2)获得各点处的布格重力异常值后，绘出平面等值线图，等值线距为0.5mGal； (3)根据异常（平面或剖面）特征，选用适当的方法进行处理（如压制干扰、消除区域场等）进行处理，并对处理效果进行描述； (4)将处理后的异常进行反演； (5)写出全部过程和所采用的处理与反演方法之应用理由。 成果：根据布格重力异常数据计算及资料处理与反演初步结果判断，该异常应由地区下一球体引起，球体埋深98.8m，剩余质量，球体中心在地面的投影点坐标为（248.8,248.8）m。 工作内容及步骤 §1-布格重力异常计算 用excel先算出各项校正（除地形校正外），需要注意的是在纬度校正中 为测点到总基点间纬向距离，由于测点都位于总基点以北，故取正值； 为总基点纬度即45°; 单位要划为km。中间层校正中 取值为2.50 。 §2-布格重力异常处理 1.绘制平面等值线图 获得各点处的布格重力异常值后，用surfer绘图软件中的自然邻点网格法绘出平面等值线图，等值线距为0.5mGal，旁边为经过九点平滑法处理后的等值线图如下： 图例： 等值线 图例： 等值线 可以看到经过平滑后的等值线光滑了很多，相当于虑去了高频成分。平滑处理的目的是消除异常数据中由观测引起的偶然误差以及由地表附近密度分布不均匀引起的杂乱无章的重力效应，以得到有意义的地质体引起的异常。 为了更好的了解该区域重力异常特性，再用surfer做出其一阶导平面等值线图和二阶导平面等值线图： 一阶导在（y=250）线上的剖面图 二阶导在（y=250）线上的剖面图 由上两图跟以下球形正演得出的平面等值线图和剖面图相比较，可以基本判定地下引起重力异常的物体形态似球，且根据二阶导平面等值线图还可以推断下面可能是由一个大球和四个小球叠加而产生的重力异常。这为后面的反演过程选用什么模型提供了依据。 下面是球体重力异常一阶及二阶导正演剖面图： 2.异常处理（分离区域异常和局部异常） 根据异常特征，我们知道该地区异常比较稳定，高频成分较少，在前面已摆出经过九点二次平滑法压制干扰后的等值线图，前后基本上变化不大，只是等值线更加光滑滑而已；而我们看到异常的导数等值线图把局部异常突出的比较好，因此我们可以用高阶导数法