重力实习报告

第一章序言1、实习目的本次重力生产实习是重要的实践性教学环节，通过教学实习将理论联系实际，巩固校内课堂理论教学知识；熟悉、掌握重力仪器的操作和使用；培养学生树立实事求是的科学态度和不畏艰苦、勇于探索的工作作风；培养学生的动手能力和处理野外实际问题的能力，为进一步学习重力勘探专业理论和实际工作打下初步的基础。

2、实习内容本次实习主要分三个阶段进行，2011年6月26日在校园内进行了仪器的检查试验和基点网的联测；2011年7月21日在校园内教一楼前进行了测线的重力测量；2011年7月22日于教五楼319室学习了数据处理及解释的方法。

本次实习以小组为单位，运用多台仪器不同时间，不同测线，不同任务进行仪器检测及测线测量，前期工作全班分为3个大组进行仪器检查试验，后期测线测量全班分2个大组，分别对A、B共2条测线进行了重复测量。

最后利用测得的实际重力资料进行资料处理、反演解释和报告编写。

分析了二连浩特重力异常特征，并对其重力资料进行处理和解释。

第二章 重力仪检查与试验1、重力仪的基本操作步骤如下（1）将电池与仪器接上，打开温度显示，对仪器内部进行电恒温处理，直到达到仪器 指定的恒温温度（大约需要2-5小时）。

（2）将仪器小心从箱中取出，轻轻放在铝盘上。

（3）将仪器纵、横水泡调节居中。

（4）打开读数灯泡，从目镜中观察摆丝位置，这时摆丝停靠在左侧终止线上。

（5）打开摆丝锁定开关，调节刻度盘，使摆丝可以自由摆动。

（6）进一步调节刻度盘，使摆丝与仪器指定的读数线重合。

（7）锁上摆丝，记录点线号、读数时间和记数器上的读数值（xxxx.xxx ）。

（8）把仪器放回箱中。

2、仪器的检查观测主要有2个方面内容：横水准器检查、调节和灵敏度检查调节。

横水准器检查与调节方法如下：（1）利用纵、横水泡调节装置，将纵、横水泡调节居中。

（2）利用刻度盘，将摆丝与仪器的读数线对齐，读出记数器上的读数。

（3）保持纵水泡居中，调节横水泡位置，使其向某一方向偏离中心位置1-2格。

（4）利用刻度盘，将摆丝与仪器的读数线对齐，读出记数器上的读数。

（5）保持纵水泡居中，调节横水泡位置，使其继续沿某一方向再偏离1-2格。

（6）利用刻度盘，将摆丝与仪器的读数线对齐，读出记数器上的读数。

（7）向另一方向重复3-5步过程。

（8）以水泡位置为横坐标，以记数器上的读数为纵坐标，画水泡曲线。

（9）如果横水泡居中时读数最大，则横水泡位于正确的位置。

否则，将仪器调到读数最大 的位置（这时水泡不居中），保持仪器不动，利用横水泡调节螺丝将横水泡调节居中。

（10）重新检查水泡位置。

检查数据如表1所示：表1 横水准器检查数据表 根据以上数据成图1如下：水泡位置 上偏2格（-2） 上偏1格（-1） 居中（0） 下偏1格（1） 下偏2格（2） 重力读数3663.9853664.0753664.1363664.1083664.062可见，横水泡居中时读数最大，所以横水泡位于正确的位置。

只是横水泡上下偏重力读数不对称，不明确这种状况对重力仪的测量精度有无影响。

3、灵敏度检查和调节灵敏度检查和调节原理是仪器的纵水泡与摆丝平行，用来确定仪器的灵敏度。

步骤如下：（1）在仪器恒温状态下，将纵、横水泡调平，打开摆丝开关。

（2）利用刻度盘将摆丝调到左侧终止线，然后再向右调，距离左侧终止线1-2个光线格。

（3）记下刻度盘位置，再顺时针旋转一圈（1毫伽）。

（4）从目镜中观察摆丝的移动。

（5）如果摆丝的移动在9-11个光线格之间，灵敏度合适，不需要调节。

（6）如果摆丝移动小于9个光线格，灵敏度低，需用如下方法将灵敏度调高：将纵水泡向远离目镜方向移动一点（向左移动），使仪器右侧降低，检查灵敏度，如果在允许的范围内，则小心保持仪器不动，调节纵水泡螺丝使纵水泡居中。

（7）如果摆丝移动大于11个光线格，灵敏度过高，这时要抬高仪器的右侧（纵水泡向右移动一点），方法同上。

步骤（2）后摆丝读数2.8，顺时针旋转一周后后读数3.6，则摆丝移动8个光格，说明此仪器灵敏度偏低，需要进行步骤（6）的调节。

4、仪器试验观测仪器试验观测方法分为静态试验和动态试验。

静态试验是对仪器在同一地点连续观测一段时间，本次试验观测地点为教五楼下花园，测量时间段为11：30-12：35，读数间隔5min；动态试验测量地点为逸夫楼1-8层，测量时间段12：45-14：11。

静态测量原始数据见表2：表2 静态试验原始重力读数值测量时间11：3011：3511：40 11：45 11：50 11：55 12：00 重力读数3664.1723664.172 3664.172 3664.173 3664.173 3664.172 3664.172 测量时间12：05 12：10 12：15 12：20 12：25 12：30 12：35 重力读数3664.170 3664.168 3664.168 3664.175 3664.175 3664.175 3664.174 最后计算零漂值，整个计算流程由重力野外处理程序完成，得数据结果如图2所示图2 静态试验处理数据结果图运用matlab绘出零漂曲线等如图3所示重力的动态试验在逸夫楼内进行，测量过程为1-8层往复测量，8层设置一个辅助点，测量数据如表3所示表3 动态试验测量原始重力读数表点号D1 D2 D3 D4 D5 D6 D7 D8测量时间12:45 12:52 12:57 13:03 13:09 13:18 13:24 13:29 重力读数3663.750 3662.481 3661.422 3660.346 3659.417 3658.284 3657.161 3656.144 点号D8 D7 D6 D5 D4 D3 D2 D1测量时间13:36 13:42 13:47 13:51 13:55 14:02 14:07 14:11 重力读数3656.283 3657.161 3658.333 3659.282 3660.505 3661.501 3662.601 3663.715 点号辅助点F 测量时间13：33 重力读数3656.291 将以上数据经过重力野外处理程序处理后，得动态试验仪器动态误差值为0.06583652。

5、基点网联测图4 基点网联测点线图基点网联测依照上图进行，3个组各测量一个闭合回路，按照单次往返观测法进行，本组观测路线为：G3-G8-G7-G6-F-G6-G7-G8-G3。

各组基点网联测数据见下表：计算后的重力读数如表5所示计算后各点边段差如表6所示表6 各边段差及公共边段差运用matlab程序计算环路闭合差及各环路平差系数结果如表7所示第三章剖面重力测量与资料整理1、剖面重力测量及数据整理对于剖面的观测，我们小组对教一楼前A测线进行了测量，测线A上布置测点30个，点距5m。

测量开始前和结束后分别对基点进行测量，且测量有一定的精度要求。

测量结束后还要求对B测线进行检查测量。

测点测量为单次测量。

测量数据如下表：表8 测线A测点重力读数表点号测量时间重力读数点号测量时间重力读数点号测量时间重力读数基点14:203663.7231015:303663.9572216:103663.866辅点14:243663.7081115:323663.8712316:143663.840基点14:293663.7161215:363663.7922416:173663.734 114:433663.8521315:403663.7592516:203663.779 214:583663.6781415:433663.9192616:223663.851 315:023663.8931515:463663.8082716:263663.811 415:053663.8681615:483663.8192816:293663.681 515:103663.8551715:523663.8322916:323663.786 615:133663.8591815:553663.8293016:363663.827 715:183663.8241915:593663.905B2816:443663.830 815:223663.9322016:043663.804B2416:483663.7179 15:253663.9522116:073663.788基点16:523663.820经过计算得各测点的相对重力值如下表所示：表9 测线各测点相对重力值点号基点 1 2 3 4 5 6 7 8 9重力值0 0.13 -0.056 0.16 0.133 0.116 0.118 0.079 0.186 0.204 10 11 12 13 14 15 16 17 18 19 200.206 0.117 0.034 -0.002 0.158 0.043 0.053 0.063 0.057 0.131 0.025 21 22 23 24 25 26 27 28 29 30 单位0.006 0.083 0.053 -0.057 -0.014 0.058 0.014 -0.121 -0.017 0.021 g.u相对重力值曲线如下图所示：图5 相对重力异常曲线图2、布格重力异常的计算布格重力异常是对测量值进行布格校正后的重力异常，布格校正包括高度校正、地形校正、纬度校正和中间层较正，由于测区地势平坦，不做地形校正。

每项校正都有相应的计算公式，可以通过matlab软件进行相应的校正计算，同时生成相应的图件，也可以将经过计算之后的结果运用surfer软件成图。

布格重力异常计算公式：Δg B=g k+Δg h+Δg T+Δgσ-gφ，g k为相对重力异常值高度校正计算公式：{Δg h} g。

u=3.068{h}m。

中间层校正计算公式：{Δgσ}g。

u=—0.419{σ} g.cm-3{h}m。

纬度校正的计算公式：{gφ} g。

u=—8.14sin2φ{D}km。

经过matlab程序计算测线B的布格重力异常曲线如下图所示：图6 测线B布格重力异常第四章重力异常处理解释方法1、向上延拓和异常分离的原理和方法由于重力场是场源到测点距离的平方的反比，因此对于深度相差较大的两个场源体来说，进行同一个高度的延拓，他们各自的异常的减弱或增大的速度是不同的。

进行上延计算时，由浅部小的场源体引起的范围小、比较尖锐的“高频”异常，随高度增加的衰减速度比较快；而由深部大的场源体引起的范围大的宽缓的“低频”异常，随高度增加的衰减速度比较慢。

因此，向上延拓有利于突出深部异常的特征，分离深浅异常。

向上延拓空间域计算公式：由于校园内实测仅限于一条测线，所以我们运用二连浩特的实测数据来实现向上延拓的处理方法，二连浩特原始异常图如下所示：图7 二连浩特原始异常图通过matlab处理，将上述异常向上延拓500m得图8和图9所示的异常。