重力场与重力勘探实验指导课件

重力场与重力勘探实验指导实验一、重力仪认识操作一、实验目的了解掌握金属弹簧重力仪的测量原理；学习操作金属弹簧重力仪，了解金属弹簧重力仪的工作范围、操作步骤、性能特点等，为实际野外重力测量实习奠定良好基础。二、实验内容金属弹簧重力仪的测量原理；金属弹簧重力仪野外测量相关情况介绍；金属弹簧重力仪的认识操作。

技术指标

D型G型测量范围2000g.u. 70000g.u.测量精度0.02g.u.约0.04g.u.零点漂移约5g.u./月 约5g.u./月 (使用1年以上)(使用1年以上)约10g.u./月 约10g.u./月 (使用1年以下)(使用1年以下)重复性约0.05g.u. 约0.1g.u.电源DC12VDC12V净重3.2kg3.2kg当重力改变时，旋转测微螺旋，使杠杆上下倾斜，再带动主弹簧，让摆杆回到零点位置。主弹簧工作原理平衡方程式mglsinα=KLd因 mglsinα=Kx1y又x1=bcos(π/2-α)=bsinα mglsinα=Kbysinα 或mgl=Kby即是该仪器测量重力变化的原理表达式。美国LCR型金属弹簧重力仪1、4－水平调节螺旋(横向)；2－纵向水准器；3－读数钮鼓及计数器；5－外壳；6－检流计；7－目镜筒；8－水平调节螺旋(纵向)；9－横向水准器；10－夹固螺旋；11－温度计窗6789101112345重力仪的基本操作步骤：将仪器小心从箱中取出，轻轻放在铝盘上。将电池与仪器接上，打开温度显示，对仪器内部进行电恒温处理，直到达到仪器指定的恒温温度（大约需要2-5小时）。将仪器纵、横水泡调节居中。打开读数灯泡，从目镜中观察摆丝位置，这时摆丝停靠在左侧终止线上。打开摆丝开关，调节刻度盘，使摆丝可以自由摆动。进一步调节刻度盘，使摆丝与仪器指定的读数线重合。记录记数器上的读数。锁上摆丝，把仪器放回箱中.实验二、重力资料整理一、实验目的掌握重力基点网平差及平差系数的计算方法；掌握重力异常值的计算方法及各种改正的地质—地球物理含义；学会使用MATLAB或C语言编写相关程序。二、实验内容1、基点网平差实验下图为三环路基点网，各边段差值如图所示，顺时针方向为正。假设各边段时间差均为1，试采用线性方程组法求解平差系数，要求用MATLAB或C语言编程求解并输出各环路的平差系数。2、扇形法地形校正值计算已知测点高程为10米，在其周围进行了三环八方位地形测量（见下图），近环0~5m，中环5~10m，远环10~20m，各点高程见下表，地表岩石平均密度为2.3g/cm3。试用MATLAB或C语言编程求测点的地形校正值。1、重力基点网平差如果没有误差存在，各边段构成的每一个闭合环路内段差相加后应为零，ΣΔgi=0由于联测中的误差，上式一般是得不到满足的，往往存在一个不等于零的偏差值，称为基点网的闭合差。产生闭合差的主要原因在于重力仪混合零点校正的不完全。三、实验原理III0.6-0.8若只有一个环路，闭合差为V，每边上观测的平均时间为ti，在按各边观测时间长短来分配闭合差时，其平差系数为：第i边上的平差值为δgi=-kti这样，该闭合环满足了V+∑δgi=0的条件。当基点网是由多个环路组成，每个环上都有一个或多个公共边时，就要求用每个环的闭合差所求得的ki来进行平差，并使同一公共边上两侧的平差值大小相等而符号相反。建立线性方程组联立求解平差系数 VⅠ、VⅡ分别为第Ⅰ、Ⅱ环的闭合差， tiⅠ与tiⅡ为观测Ⅰ、Ⅱ环某边的闭合时间(包含公共边的tFC)， kⅠ、kⅡ为待求的平差系数。扇形法地形校正公式：布格重力异常是对观测值进行地形校正、布格校正(高度校正与中间层校正)和正常场（纬度）校正后获得的。布格重力异常值公式：其中，h—测点海拔高程或测点与基准点的高差，m；σ—中间层平均密度，g/cm3。φ—总基点纬度或测区平均纬度；D—测点到总基点间纬向（南北向）距离，km。Δgk—重力观测值，g.u.。

实验三、重力异常正演一、实验目的通过球体、水平圆柱体和铅垂台阶重力异常的正演计算，掌握简单规则形体重力异常正演的计算方法；通过计算认识球体、水平圆柱体和铅垂台阶重力异常的一般分布规律，了解影响简单规则形体重力异常的主要因素（如形体的几何参数和密度参数等），培养学生实际动手能力与分析问题的能力；学会使用MATLAB语言编写程序绘制异常图件。二、实验内容用MATLAB或C语言编程实现二度球体、三度球体重力异常及其导数正演，使用MATLAB语言编程绘制异常曲线或等值线图。用MATLAB或C语言编程实现二度水平圆柱体重力异常及其导数正演，使用MATLAB语言编程绘制异常曲线图。用MATLAB或C语言编程实现二度铅垂台阶重力异常及其导数正演，使用MATLAB语言编程绘制异常曲线图。三、实验原理1、二度球体、三度球体重力异常及其导数正演三度球体二度球体D=100mR=50mσ=1两个球体组合模型理论重力异常2、二度水平圆柱体重力异常及其导数正演D=60mR=22.5mσ=13、二度铅垂台阶重力异常及其导数正演∞h=20mH=60mσ=1四、实验报告内容包括实验目的、实验内容、实验原理、计算程序代码、实验结果、结果分析或小结。实验四、重力资料处理一、实验目的通过重力异常频率域滤波处理实验，掌握频率域滤波处理方法原理；通过简单规则形体重力异常特征点反演实验，掌握特征点反演方法原理。二、实验内容利用频率域滤波的MATLAB程序，对简单规则形体重力异常进行频率域向上延拓、水平导数和垂直导数处理并绘图；从第三次实验得到的重力异常中提取出任意一个有意义的异常剖面数据，试用特征点法反演密度不均匀体的几何和物性参数。三、实验原理1、频率域滤波处理方法基本原理应用傅立叶变换，计算空间域观测异常gB(x,y)的频谱GB(u,v)；将GB(u,v)与所求转换异常的滤波算子H(u,v)相乘，得到转换异常的傅立叶变换GFT(u,v)；将GFT(u,v)进行傅立叶反变换，便得到所求的转换异常gFT(x,y)。解析延拓向上延拓向下延拓

S为压制高频干扰的滤波因子。一次导数

x方向y方向z方向垂向二次导数根据异常曲线上的一些点或特征点（如极大值点、零值点、拐点）的异常值及相应的坐标求取场源体的几何或物性参数;仅适用于剩余密度为常数的几何形体。应用条件对异常作平滑处理，尽量准确确定原点的位置；对异常曲线作分离处理，获得单纯由研究对象引起的异常；对剩余（局部）异常进行分类，判明该异常的场源体接近于何种可能的几何形体，然后选用相应的反演公式。2、特征点法反演方法基本原理异常曲线形态分类第一类是单峰异常，零值点在无穷远处 如球体的Δg曲线、台阶的Vxz曲线等；第二类是具有极大值、极小值和一个零值点 如球体的Vxz曲线、台阶的Vzz、Vzzz曲线；第三类是具有一个极大值、两个极