第五章：外业重力测量.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,1,、重力测量发展简史,2,、重力测量分类,3,、相对重力测量,4,、绝对重力测量,5,、重力基准,6,、重力仪简介,第五章,重力测量,重力测量发展简史,大约,16,世纪末期，意大利物理学家,Galileo,发现了重力加速度,1673,年：惠更斯首次建立了数学摆,1733,年：拉康达明用物理摆测出精度为,1,加的重力值。,19,世纪：贝塞尔提出了可倒摆原理,二十世纪,30,年代：哈特利和格拉夫研制了金属弹簧重力仪,重力测量发展简史,20,世纪,50,年代：发展了一门新兴学科,卫星重力学,目前的主要研究方向静态重力测量、时变重力测量,（重力场潮汐变化，重力场非潮汐变化）,重力测量分类（按仪器观测位置分）,陆地重力测量（采用汽车、火车、飞机为交通工具，实际观测在地面进行）,海洋重力测量（实际观测在船上进行）,航空重力测量（实际观测在飞机上进行）,卫星重力测量（,CHAMP,，,GRACE,，,GOCE,）,重力测量分类（按观测方法）,绝对重力测量（动力法）,1,、自由落体法,2,、摆法,相对重力测量（动力法和静力法）,相对重力测量,1,相对重力测量的目的,2,相对重力测量的工作原理,3,重力观测点的选埋,4,测前仪器格值的标定,5,重力测量的外业观测和记录,6,数据处理和成果分析,二、地面位移与地球内部应力场关系,将应力与位移的关系代入得内部应力公式为：,（,5,）,式中,二、相对重力测量原理,测量方式：,绝对重力测量和相对重力测量,（,1,）绝对重力测量,：,用仪器直接测定地面上某点的绝对重力值。地球表面上的绝对重力值约在,978,983Gal,。,（,2,）相对重力测量,：,用仪器测定地面上两点之间的重力差值。地球表面上的最大重力差约为,5000mGal,。,（,3,）固定台站重力测量：,观测重力随时间的变化。,（,4,）流动站重力测量,：,观测重力随空间位置的变化,重力测量的目的,1,科学研究：（,1,）地震监测；（,2,）火山监测（,3,）地壳变形监测；（,4,）冰川监测,2资源勘察：（,1,）石油勘察；（,2,）天然气勘察，（,3,）矿产勘察等,3,工程应用,（,1,）地下溶洞、暗穴勘测（,2,）,（,2,）地下堵塞水管探测等,从上述推导过程和式,(5),看出，只要知道地面,S,上边界条件,(,如位移函数、面力函数,),和体力函数，我们可通过积分方程,(5),计算弹性地壳内任意点,P,的应力。从以上的理论推导和分析可知,二、地面位移与地球内部应力场关系,动力法：,观测物体的运动状态以测定重力，可应用于绝对重力测量和相对重力测量,。,静力法：,它是观测物体受力平衡，量测物体平衡位置受重力变化而产生的位移来测定两点的重力差，该方法只能用于相对重力测量。,重力测量原理,静力法测定相对重力,基本原理：,弹性体在重力下发生形变，而弹性体所受到的弹性力与重力平衡时，则弹性体处于某一平衡位置，当重力改变时，则弹性体的平衡位置就有所改变观测弹性体两次平衡位置的变化就可以测定两点的重力差。,垂直型弹簧重力仪：,m,：,荷重的质量,：,弹簧在无荷重作用是的长度,：,第一点上弹簧的长度,与其相应的重力,:,第二点上弹簧的长度,与其相应的重力,k,:,弹簧的弹性系数,其中,c,称为格值,振摆法测定相对重力,或,这里,，,。,重力网的布设、测网类型,测网的布设应根据科学目标和工程任务而定,测网类型,国际基准网；国家基准网；区域网；局部网,测点选择和标志设立,1,根据地质、地形交通等资料进行室内设计,2,重力点应选择在基础稳定且振动和其它干扰源小的地方,3,、相对重力点的观测平台可选用基岩、专用水泥标石、固定建筑物等,4,观测点应绘制点之记、办理委托保管书、并拍摄观测点图象资料等。,测点选择和标志设立,5,、重力点上应设置,水准标志,，并在平台表面表明,磁北方向,6,有条件的应直接测定测点的,经度、纬度和高程,，否则其平面坐标可在小于,1,：,5,万的比例尺的地形图上量取。,7,、重力测量工作结束后，应提交下列资料,（,1,）重力点点之记（包括点位交通略图）,（,2,）测量标志委托保管书,测点选择和标志设立,（,3,）全网平面位置图。,（,4,）技术总结,测点选择和标志设立,测点选择和标志设立,重力仪的标定基线及仪器标定,重力仪的标定基线分长基线、短基线和垂直基线。长基线可选用国家重力网（,85,）中的哈尔滨,-,北京,-,郑州,-,武汉,-,广州、短基线有北京灵山、江西庐山两条、垂直基线可选用武汉垂直基线,相对重力仪的标定,1,选用的标定基线应能覆盖测区的所有读数段。,2,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,重力测量的外业观测和记录,数据处理和成果分析,数据处理和成果分析,绝对重力测量,自由落体测绝对重力,(320Km),振摆测定绝对重力,绝对重力测量,自由落体测绝对重力,根据牛顿运动定理物体在下落过程中的运动方程为：,物体在下落过程中（真空）,（,1,）,（,2,）,将（,2,）代入（,1,）得：,（,3,）,对（,3,）式两边积分得：,式中 ，,和,g,为待定常数。,从（,3,）式可以看出，我们只要在,3,个不同的时刻测定,h,的大小，就可通过解三元一次线性方程求出绝对重力,g,。,自由落体三位置法：,其中,，,，,，,通过上式，我们可利用三位置观测确定绝对重力值,g.,对称自由运动的方法（上抛法）：,如图所示：将物体以初速度,v,0,上抛，测定物体在上抛和下落过程中，在同一高度时的时间差,t,1,和另一高度时,t,2,同时测定高差,h,。,的时间差,上抛过程：,上抛速度：,垂直距离：,在顶点时：,时，,下落过程：,通过比较可以看出，在上升和下落过程中通过测定的速度值相同，在上升和下落过程中经过位置之间的时间差：,当,时，物体为下落阶段,而,两边同时乘以,2g,，通过比较可以看出：,重力,和,与观测值,之间的精度关系,：,对上式两边取对数微分可得，,应用误差传播定理得：,若要求,，,则按等影响原则有：,如果物体下落距离,，下落时间,则长度测量误差应不超过,1,微米，时间量测误差不超过,振摆测定绝对重力,如图所示,将质量为,m,的物体悬挂在,o,点,在无空气阻力的情况下,物体在给一作用力后自由振荡,物体在切线方向受力为,mgsin,为物理摆的摆动角速度，,角加速度,物理摆的运动方程为：,o,线加速度,根据牛顿定理,振摆测定绝对重力,两边同乘以,振摆测定绝对重力,代入初始条件,t=0,=0,t=/4,得,振摆测定绝对重力,通过计算,最后可得,这里,，称为改化摆长。,求解上述微分方程可得,：,重力,之间的精度关系：,与观测值,和,假定对的影响相等，并要求重精度为,1mGal,，即,则摆动周期的允许观测误差为：,改化摆长的允许观测误差为：,也就是说，如果要求重力测量达到,1mGal,的精度，则当摆动周期为,1s,时，周期的观测误差不得超过,当改化摆长为,1m,时，其量测误差不得超过,1,微米。,绝对重力测量结果误差分析与改正,对于实时观测的重力值，由于它不仅包括地球质量的引力，还包含有日月引力及其产生的地球形,四、汶川强震区的,GPS,资料计算结果分析,由,GPS,资料计算得四川地区地壳,15km,深处剪应力年变化分布如图,4,所示,.,从图,3,和图,4,整体上看,整个区域的主压应力变化呈现出从西南向东北变化的格局,这与中国大陆受印度板块的挤压，方向基本一致,.,从局部看,在汶川强地震活动区的附近,沿断层构造的应力场变化明显大于其它地区,呈现沿龙门山分布的变化特征,.,这说明了局部应力,.,变的潮汐效应，环境变化以及仪器系统的误差等，这些因素影响到观测重力值。因此，需要对实时观测重力值可能产生的误差作分析及必要的改正，消除已知的各种影响因素引起的实时观测重力值变化。,光线对重力垂线的偏离,用激光干涉仪测量自由落体移动距离时，如射向垂直下落落体的光心的光线与重力垂线有一小角度,的偏离，可以,推算出重力测量的误差：,确的重力改正量。光线垂直性误差的影响是使测量值偏小，,测定,后，便可算出正,四、汶川强震区的,GPS,资料计算结果分析,通过理论推导和实际,GPS,观测资料的计算分析,本文得出如下结论和认识,:,根据格林第二内部积分公式,导出了根据地面位移场速率观测值计算地球内部应力场变化的近似计算公式,.,通过四川地区,(2001-2005,和,2004-2007,年,),间观测的,GPS,观测速率资料计算结果表明,主压应力变化较大的区域主要为龙门山断裂附近的地震活动区,这与中国大陆,应作此项改正。,2.,振动干,扰,地球表片存在着复杂频谱的微震，它对重力测量用的干涉仪的参考立体棱镜干扰很大。给测量重力值,g,带来误差。此类误差是随机性的；随测量次数增加而逐渐减小。因此,MIM-,型绝对重力仪采用的多位置符合振动检出数据处理方法防止（或消除）大地为震动对测量重力值,g,的影响。而,JILA,型和,FG5,型等国际上的绝对重力仪均采用惯性悬挂参考棱镜的方法消除微震的影响。,3.,时间频率标准,在绝对重力测量中，时间的测量是一个重要的标准量，它直接影响测定重力值,g,的精度，现在先进的绝对重力仪（,NIM-,型、,JILA,型和,FG5,型等）是用准确度优于,1,，稳定度,的铷原子频标为时间频率标准。所以完全满足,优于,3,测量重力的精度要求。,4.,长度标准,长度标准同样是直接影响测量重力值的精度的重要的数据。绝对重力仪一般都是应用迈克尔逊干涉仪测定落体自,由下落距离长度的。干涉仪的光源是碘吸收稳定的氦氖激光，,其波长准确度为,1,，稳定度优于,1,，复现性为,4,满足测量重力值,g,的要求。,5.,空气阻力,经抽高真空后，真空室内的空气浮力，压差阻力及内摩擦力完全可以忽略。,6.,静电影响,真空室内部静电场的来源有两方面，意识落体悬挂接触定位处的两部分元件采用同一种金属材料。另一是裸体传动机械的摩擦起,自转服务,（,IERS,）,提供的极坐标。,7,、重力垂直梯度转换高度改正,在经过以上各项改正后，所得到的重力值是落体质量在顶部位置的重力值。一般作为相对重力测量起始点的绝对重力点，需将重力值转换到地面重力标志上，或者为了与其他类型的绝对重力仪的观测值进行比较，将重力值转换到特定的高度如,0.8m,或,1m,处。由于测点点位收地形起伏、地下密度不均匀体以及点位周围建筑的影响，根据实测的重力垂直梯度可将,1.30m,处的重力值为归算地面。,电荷具体的结构、材料选用有关，同时与摩擦速度、接触时间、压力与张力的积累状态有关，这在仪器设计、实验过程中，应尽量采取必要的措施以消除诸因素影响。使静电对重力测量的影响控制在小于,1 ,以使该影响可以忽略不计。,8.,磁性,磁性对于测量重力,g,的影响，一个是落体剩磁与真空室内各种咋散磁场的磁力作用，另一个是涡流阻力。为减小磁性影响，真空室内的零、部件如落体、上部法兰盘等材料的选择，必须采用超级纯铁，无赐不锈钢。同时对仪器各种杂散磁场要严加控制，使磁力影响控制在很小的范围内。,9.,极移改正,极移改正是由于地球自转与台站测点之间距离随时间变化而导致的离心力变化。因此，需要根据最接近观测时间的磁极位置进行重力改正。此项改正公式为：,是极移改正,（）,，是地球旋转角速度,弧度,/,秒）,，,（,7292115,是地球长半轴,（,6378136m,）,，,和,分别是观测点的大地纬度和经度，,x,，,y,是国际地球,10.,局部气压改正,对于每次落体测量所观测的重力值均归算到台站正常大气压时的重力值，测点局部气压变化的重力改正,为：,A,为气压变化对重力的影响系数，通常在,（,0.30,0.42,）,.,是观测点的局部气压。,是根据测,采用,DIN5450,标准大气模型计算公式得：,各测点的,值进行实测后，以实测值作改正用。,11,、仪器有效高度改正,对于每一次落体观测的重力值均被归算为落体顶部处的重力值，此项改正是通过配置一下四阶多项式进行改正：,在此，,是观测重力垂直梯度；,分别是落体初始,时的位置和速度；,分别是任意观测时刻的观测,观测时刻及落体位置，仪器的有效高度改正通,常大约为,-25,.,速有限改正约为,-11,。,12,、空间重力梯度,空间存在着重力垂向梯度，以,表示。如落体悬挂在顶,在下落,t,时刻时，落到高度,h=,处的重力值为,:,因,零，故下落,h,高度所需的时间就比,=,零时所需时间,要小。使精确测定,g,时产生误差。由于,是一般平均值，与各地实际,值有差异。所以，必须对,变化而随时间变化。因此，在精确地确定地面点 的重力值,g,时，必须作固体潮汐改正以消除其影响。,13,、光源有限改正,由于光速不是无穷大而是有限的，当落体在真空中下落时，光电倍增管接收到的该落体反射回来的信号要比落体反应出来的信号延迟一个时间,.,是光原子从落体反射到光电倍,增管接收器上所需的时间。若光程距离为,S,，,则,=,会引起精确测距误差，从而使观测重力值比实际重力值,要大。因此，观测重力值需要加上光速传播时间的改正。光,激光光源、光电倍增管、放大器等器件不可避免地存在着噪声，为提高测量精度，要求有高的信噪比，为了提高信噪比，就需要选择噪声低的激光管，并在测量前仔细调节放电电流，使其工作在噪声最小值内；同时，在设计激光干涉仪时要尽可能提高光源的功率、减少光路的损失。,14.,信噪比,15.,地球潮汐改正,根据万有引力定律所有天体对地球每个质点都有引力作用，因而对地球表面形成引潮力，引潮力的垂直分量与重力方向一致，引起体面上每个点的重力值因地球与天体位置的,三、重力基准,世界重力基点：,世界公认的一个重力起始点,维也纳系统（,1900,年,-IAG,）,波茨坦系统（,1990,年,-IAG,1894-1904,年）,国际重力基准网：,1956,年,IAG,决定建立世界一等重力网（,FOWGN,）,1967,年,IAG,决定在波茨坦绝对重力值中加上,-14mGal,作为新的国际重力基准,1971,年,IUGG,决定采用,IGSN71,代替波茨坦国际重力基准,新的波茨坦国际重力基点的值为,引潮力,M,E,d,设地球质量为,E,，月亮质量为,M,，地球质心与月亮质心之间的引力为,月亮对地球上单位质量的平均引力为,P,引潮力,月球对地球质心的引力为：,月球对地球上,P,点的引力为：,因为,月球对,P,点的引潮力为月球对,P,点的引力与月求对地球质心,O,点的引力之差。,引潮位,月球对地心的引力位为,月球对,P,的引力位为,月球绕地球旋转产生的惯性离心力位差为,O,M,P,d,r,z,中国的重力基准网：,在全国范围内提供各种目的重力测量的基准和最高一级控制,中国曾在,1957,年建成第一个国家,57,重力基本网，它的平均联测精度为,:,1985,年中国又新建了国家,85,重力基本网，其平均联测精度较之,“,57,网,”,提高一个数量级，达到 的精度,该网改正了波茨坦系统的系统误差，增测了绝对重力基准点，加大了基本点的密度。,中国,2000,重力基本网,覆盖了中国的全部领土（除台湾外，包含南海海域和香港、澳门特别行政区）。全网由,133,个点组成，其中有,17,个基准点（绝对重力点）和,116,个基本点（相对重力联测点）。为便于今后联测和作为基本点的备用点，对,106,个基本点每点布设了一个引点。,重力基准点的观测精度优于,重力基本点的相对观测精度优于,，,平差后重力基本网的中误差不大于 。,四、主要重力仪简介,性能指标,：,Lacoste&Romberg G,型重力仪（美国），其测量精度为,40微伽，重复度为,10-20,微伽，测程为,7000mGal,。,特点,：精度高、观测时间短、体积小重量轻（,3.2kg,）,观测成果计算简单。,Graviton,EG,型重力仪,General Specifications,.Total shipping weight 320kg in 6 containers,.Total Volume:1.5 Cu.m,.Floor space requirement 3 sq.m,.Input voltage:110-240 VAC(switchable)50-60HZ,.Nominal power requirement 500w,Performance Specifications,.Accuracy.2,.Operating dynamic range:world-wide,.Operating temperature range:10 deg C to 30 deg C.,GEOPHYSICAL RESEARCH,.Detection of vertical crustal motion,.Complementary verification of displacements measured with GPS and VLBI,.Determination of the geoid,.Volcanic magma flow monitoring(i.e.mammoth lakes),.Postglacial rebound studies,.Uplift or subduction studies,.Earthquake research,.Long period tidal monitoring and modeling of earth anelasticity,ENVIRONMENTAL MONITORING,.Water table monitoring in deep and/or multiple aquifers,.Nuclear waste management and cleanup,.Global sea level studies for global warming,EXPLORATION AND RESOURCE MANGEMENT,.Oil exploration,.Mineral exploration,PRECISION MEASUREMENTS AND CALIBRATIONS,.Pressure transducer and load cell calibration,.Redefinition of the kilogram in the SI system of units,.Big G determinations and the equivalence principle,.Calibration of superconducting of other high precision relative gravity meters,INERTIAL VAVIGATION,.Gravity reference station determinations,.Relative gravity network control points,.Establishing geodetic tie points for gravity networks,.Defining the geoid,A-10,绝对重力仪,Instrument Features,.Quick setup,.Easy to operate,.Rugged,.Fully automatic operation,.Automated levelling,.Battery operated,.Temperature controlled sensor,.Ideal roadside operation from a vehicle,.Real-time gravity reduction,General Specifications,.Accuracy:10 Gal,.Precision:10 Gal in 10 minutes,