绝对重力仪研制中一种新的自由落体轨迹重建算法

真空系统中落体做自由下落运动的轨迹重建算法,是进行拟合求解绝对重力测量结果的关键步骤之一.在自主研制小型高精度激光干涉绝对重力仪的过程中,基于目前高速发展的数字测量技术,提出一种自由落体轨迹重建的算法.其基本原理是:首先通过高速采样,得到落体自由下落过程中产生的一系列数字化的激光干涉信号;其次通过数字处理算法提取干涉信...     

基于激光干涉的新型高精度绝对重力仪

本文介绍了一种自主研发的新型高精度绝对重力仪样机,它利用激光干涉原理测量落体自由下落时相对于参考棱镜的时间位移信息,通过振动误差补偿算法提取落体下落轨迹上的时间位移坐标,然后基于最小二乘拟合法计算测点的重力加速度.在白家疃国家地球物理观象台国家重力基准点进行的对比观测试验结果与误差分析确定,样机经误差修正后的系统偏差为(-27.46～-30.59)±3.06×10-8 ms-2,2小时测量精度优于10×10-8 ms-2.     

利用FG5绝对重力仪观测数据确定重力垂直梯度

本文建立了利用FG5实测数据求解重力垂直梯度的数据处理模型与算法.通过对多次自由落体实验的下落距离拟合残差叠加求均值,发现下落距离观测量中存在明显的有色噪声.通过对有色噪声的建模,并以剩余残差为依据选取可靠的下落时段,解算测站点的重力垂直梯度.利用本文所提出的数据处理方法分别对FG5-214绝对重力仪在两个测站上的观测数据进行处理,以相对重力仪测量的重力垂直梯度结果为参考值,本文处理得到的重力垂直梯度结果相比于未考虑有色噪声并依据经验选取下落时段的解算方法得到了显著改善.     

点源位错引起的重力、位势及其梯度变化

本文研究了点源位错引起的重力、位势及其梯度变化，完整地从理论上推证点源位错引起的位势和重力变化；对点源的处理数学简单、物理意义明确，避免了对积分的奇异性估计，使该问题理论上更趋完善．同时，分析了各型位错下空间重力及其梯度变化、地面上的重力变化、重力垂直位移梯度和视重力垂直位移梯度的特征．分析表明:① 点源位错引起的空间重力及其梯度变化随垂直位置变化很小；② 点源位错引起的重力变化比梯度变化大的多，重力梯度变化是不显著的；③ 质量再分布引起的重力变化可达总重力变化的一半左右，且走滑、倾滑垂直位错引起的总重力变化与垂直位移形变场成正负反像分布；④ 重力垂直位移梯度不是常数，对于各型位错具有不同的分布特征；⑤ 点源位错引起的重力垂直位移梯度和视重力垂直位移梯度近似相等．      

地面扰动重力垂直梯度的确定

鉴于扰动重力垂直梯度在大地测量和物探中有重要作用,在未直接测量的情况下,对如何由地面重力异常及地形数据求取扰动垂直梯度进行了研究和分析,认为在被研点附近要求重力点分布密集,且精度不低于1×10-6ms-2;对中央区域积分奇异性问题也作了讨论;此外,扰动重力垂直梯度如何在大地测量和地球物理中进一步得到应用,以及为何用扰动重力垂直梯度代替重力异常垂直梯度也给予了简明的论述.     

绝对重力仪设计中的一种新型落体控制方式

为了将绝对重力仪在真空筒中自由下落落体的运动传动控制结构在大气压条件下实现,便于检修和维护,我们提出了一种全新的落体控制方法:焊接波纹管法.从理论计算上升段和下降段电机对拖架控制过程的速度、位移曲线以及试验结果可以看出,这种焊接波纹管驱动拖架的方法实现了对落体伺服跟踪的同时,也满足了对落体自由下落时间的设计要求.     

中国几处城市高楼重力垂直梯度测定

重力垂直梯度测定是应用地球物理学和物理大地测量学的一个重要的实验性课题。研究重力垂直梯度（以下简称垂直梯度）在地球表面的分布与变化,可以增加对地球重力场和地下地质构造及其活动性的了解。1982-1984年,我们在研究垂直梯度测量原理和     

航空重力测量研究进展

近十几年来,航空重力测量技术的研究和应用日趋活跃,业已成为地球重力场研究中最为热门的领域之一．航空重力测量可分为航空标量重力测量、航空矢量重力测量和航空梯度重力测量．本文概述了航空标量重力测量的研究现状及发展动态,传统航空重力测量系统的改进,新型系统方案的研究,从GPS中导出加速度的理论和方法．滤波和估算方法及航空重力数据的应用等五个方面的重要研究成果．文中最后介绍了我国自行研制的航空重力测量系统的概况。     

绝对重力仪测量结果的振动划线处理方法初步研究

针对自主研发的高精度绝对重力仪样机测量结果中引入的振动误差, 提出了一种全新的振动抑制算法．本算法利用一簇斜率相同的平行直线与记录到的参考棱镜振动曲线相交, 通过提取这些交点的时间信息, 确定落体下落轨迹上进行重力加速度计算的时间和位移坐标． 然后, 利用最小二乘算法获取测点的绝对重力加速度值．对不同重力加速度值的适用性和不同类型振动干扰的抑制效果模拟计算发现, 对于优势频率大于0.2 Hz, 峰峰值小于500 counts的干扰, 算法可以有效提高测量结果的系统偏差和测量精度．最后, 实际实验样机的测量结果表明, 本算法使最终测量精度提高了一倍, 为绝对重力仪的小型化和流动测量提供了软件技术支持。      

利用激光雷达测量重力波三维结构

激光雷达观测得到的密度、温度等数据被广泛应用于大气重力波研究.瑞利激光雷达可以获取激光路径上的大气密度、温度数据.对于大气中的三维波动而言,单条路径上的观测参量能提取得到的波动信息有限.本文首先以单色重力波为例,分析了利用激光雷达直接观测三维波动结构的可行性.激光雷达垂直观测即可得到重力波的垂直波长,当激光雷达以一定的天顶角斜向测量时,所得到的波长包含了重力波的垂直波长以及水平波长信息.因此,利用激光雷达同时以三个方向(垂直、向南(天顶角30°)以及向西(天顶角30°))测量,可以提取得到重力波的垂直波长和水平波长.本文利用中国科学院国家空间科学中心研制的车载532nm瑞利激光雷达的经向系统和纬向系统同时以不同的指向角观测大气重力波,对利用激光雷达获取三维波动结构的方法进行了分析研究.本文给出了北京地区激光雷达观测重力波的诸多案例,分析了30~60km高度范围内北京地区大气重力波的垂直及水平波长信息.并以2017年11月7日观测的准单色重力波为例,结合再分析资料的风场数据,分析了该重力波的水平波长,垂直波长及传播方向等信息.     

RECENT DEVELOPMENTS IN UNDERGROUND GRAVITY SURVEYS1

Underground gravity observations in deep coal mines using the conventional gravity meters Worden (type Master) and LaCoste-Romberg (model D), both of which have been adapted to the fire damp regulations, can be accurate to ± 10 and ± 3 μgal, respectively. For underground determination of the vertical gradient of gravity the LaCoste-Romberg meter is used together with a specially designed measuring tower. Using this euipment an accuracy in tower gradient observations of ± 30 E was obtained. To apply the equipment to precision gravity observations in underground situations an additional correction, i.e. a gallery correction, is needed. High accuracy in correction is achieved by a new method of three-dimensional modelling. The gravity effect is computed for bodies with a surface approximated by triangular elements, which are generated from corner points of the body. The combination of gallery correction with tower gradient data leads to a new method for in situ density determination. It offers the possibility of horizontal instead of vertical density profiling. To demonstrate the effectiveness of the developments in underground observations the localization of a pump room is presented. Microgravity and tower gradient observations were carried out to detect the cavity. The horizontal gradient was also calculated to give a more reliable location.     

Three-Dimensional Gravity Modeling In All Space

We review available analytical algorithms for the gravity effect and gravity gradients especially the vertical gravity gradient due to a right rectangular prism, a right polygonal prism, and a polyhedron. The emphasis is placed on an investigation of validity, consistency, and especially singularities of different algorithms, which have been traditionally proposed for calculation of the gravity effect on ground (or outside anomalous bodies), when they are applied to all points in space. The rounding error due to the computer floating point precision is estimated. The gravity effect and vertical gradient of gravity in three dimensions caused by a cubic model are calculated by different types of algorithms. The reliability of algorithms for the calculation of gravity of a right polygonal prism and a polyhedron is further verified by using a regular polygonal prism approximating a vertical cylinder and a regular polyhedron approximating a sphere, respectively. By highlighting Haáz-Jung-Plouff and Okabe-Steiner-Zilahi-Sebess' formulae for a right rectangular prism, Plouff's algorithm for a right polygonal prism, and Gouml;tze and Lahmeyer's algorithm for a polyhedron and removing their singularities, we demonstrate that these formulae and algorithms can be used to model the gravity anomaly and its vertical gradient at all possible computation positions.     

位场垂向梯度最佳自比值的边界检测技术

位场梯度换算在地质体边界检测中有着重要的应用.但传统的梯度算法易受干扰影响,计算稳定性差,且很难在复杂的叠加异常中识别出小型地质体的边界.鉴于此,本文给出了自比值的定义,提出了能够处理高阶导数的位场垂向梯度最佳自比值的边界检测方法,阐述了方法的数学含义和物理意义.模型试验表明,垂向梯度最佳自比值算法不仅计算稳定性强,而且能清晰地检测出传统梯度算法无法检测的模型体边界.在地质条件复杂的鸭绿江盆地的重力异常实例应用中,垂向三阶导数最佳自比值计算结果识别出的构造边界与实际地质体分布有着较好的对应关系,这不但与前人的工作成果互为佐证,而且自比值圈定的负异常分布区能较好地反映出浑江煤田的工作范围.     

基于等效源的重力异常空间导数计算方法及其在川滇地区的应用

在综合分析前人研究成果的基础上, 采用基于贝叶斯准则的等效源模型滤波增强方法来识别地质体边界。 通过理论模型实验讨论了在相同的噪声水平情况下, 不同算法对于计算结果的影响, 模型实验结果表明, 与直接计算和传统高斯滤波增强技术相比, 新算法可以有效消除传统算法中普遍存在的对于噪声的过度敏感问题, 所得到的地质体边界形态失真更小。 随后, 利用此方法对川滇地区的重力数据进行了处理, 并对该地区的重力水平总梯度(THD)、 重力垂直梯度模(AVG)、 重力解析信号模量(AAS)、 重力水平梯度垂向导数(VGHD)应用效果进行了对比分析。 研究结果提供了川滇地区的构造分区信息, 对该地区构造演化、 边界定位及动力学过程分析具有重要的指导意义。     

基于水平方向解析信号的均衡重力位场边界识别方法

边界识别技术是位场数据解释中一项基本的工作,现有的边界识别方法多存在边界识别结果发散和不能均衡深浅地质体异常的缺点.目前一些均衡边界识别方法会因正负异常同时存在而引起额外的错误边界或者存在人为主观因素去除错误边界信息的缺点.本文充分利用重力位场张量梯度的多信息成分,提出利用水平方向解析信号及其垂向导数与传统的均衡边界识别方法做结合的方式定义新的探测边界的方法.通过理论模型试验证明新方法同传统方法相比,能够更加清晰、准确的圈定出深浅地质体的边界.最后将新的边界识别方法应用到实测重力异常数据解释中,取得了良好的边界识别结果并能够发现更多的构造细节.     

Allowance for the vertical gradient in an airborne gravity survey

When processing the results of an airborne gravity survey flown in the Arctic at flight heights of 2000 to 3000 m, the authors encountered conditions, when the values of the free-air anomaly of gravity exceeded the values taken from the anomaly map, with a trend towards “heavier” values. The formula for calculation of the vertical gradient as a function of the latitude of the locality has been refined. The tentative calculations reported in the paper also indicate that a free-air gravity anomaly above the mountainous areas will contain additional positive systematic components, since only the effect of the vertical gradient of the normal gravity field is taken into account here, and the flight’s height during an airborne gravity survey, above a mountainous area, should be held constant, otherwise it would be necessary to take into account additional systematic biases, obtained in different flights, which are caused by the different values of the vertical gradient at different heights.     

基于加速度计重力梯度仪分析与设计

目前,重力梯度仪成为重力仪器的主要研究方向,国内尚无此种仪器问世.本文利用加速度计在旋转载体中测量公式,推导出基于加速度计重力梯度仪测量公式.分析利用此公式进行重力梯度测量时需要注意的问题.然后以美国Bell/Aero公司生产的重力梯度仪为例,导出旋转加速度计式重力梯度仪测量结果.最后总结出基于加速度计重力梯度仪设计过程中存在的各种误差,并给出解决这些问题的方案.为国内重力梯度仪的研究工作提供理论依据.     

Experimental errors of measurements of vertical geomagnetic gradients in the stratosphere

The errors of measurements of vertical geomagnetic gradients at altitudes of 20–40 km, using a balloon magnetic gradiometer with a 6-km-long measuring base oriented along gravity, have been studied in the work. The errors related to the deviation of the measuring base position relative to the vertical have been studied during the real balloon flight with the help of the navigation GPS receivers. The deviations of the measuring base within 5°, which can sometimes reach 15°, have been obtained. This results in a decrease in the magnetic gradient measurement accuracy due to the errors introduced in the specification of the normal magnetic field used to detect magnetic anomalies. To eliminate this error, a GPS receiver was built in each magnetometer in order to observe magnetometers during synchronous measurements and to correct the measurements for the normal magnetic field. It has been indicated that the effect of deviations of the measuring base position on the results is not more than 2% of the measured value at such organization of a gradiometer.     

preliminary discussion on horizontal gravity gradient and its application to seismic gravity precursor research

In this paper, according to the synthetic gravity anomaly of a horizontally infinite cylindrical geologic body, gravity gradient in horizontal direction was calculated by potential field discrete cosine transformation in frequency domain. In the calculation, the minimum curvature method was used to extend edge lines. We found that the gravity gradient field from the potential field transformation was dependable by comparison with synthetic gravity gradient, except the data in the edges. Then, the accumulative horizontal gravity gradients before Lushan M_S7.0 earthquake were calculated for the accumulative gravity anomaly from September 2010 to October 2012. In the north-south direction, gravity gradient in Daofu-Kangding-Shimian and Markang-Lixian-Lushan exhibited a positive high value, and the strike of the high value zone was in line with the strike of Xianshuihe Faults and Markang Faults. In the east-west direction, high value zone was not as obvious as that in the north-south direction. Gravity gradients in the direction along and vertical to the strike of Longmenshan Faults were calculated by the definition of directional derivative. In the along-strike direction, high gravity gradient values appeared in Markang-Lixian areas along Markang Faults and Daofu-Kangding-Shimian areas along Xianshuihe Faults, and extremum appeared in Kangding-Shimian and the area nearby Lixian. In the direction vertical to the strike of Longmenshan fault zone, high gravity gradient values appeared in Lixian-Lushan-Kangding-Shimian areas, and the extremum appeared in the area nearby Kangding. The results indicate that gravity gradient in the direction along and vertical to the strike of faults can better show the relative gravity change on the two sides of faults. Lushan M_S7.0 earthquake is located at the transition zone between the two high value zones of gravity gradient. The total horizontal gravity gradient shows that the location and strike of the high value zone are basically consistent with regional faults, and the extremums of total horizontal gravity gradient appeared nearby Lixian, Kangding and Shimian.