离散时变重力数据的可视化、指标量定义与解释

针对流动重力观测获得的数据成果特点,提出一种基于流动重力段差变化的时变重力数据可视化方法,并定义了两个指标量G值和C值,用以评价区域性重力场变化的显著性程度.在此基础上,应用该方法处理和分析了首都圈地区的流动重力数据.结果表明:与传统采用等值线方式来刻画时变重力场特征不同,该方法更能突出发生变化的重力测点位置、测量误差及其变化的显著性程度,可为研究与地震孕育、发生过程有关的重力场变化提供更多的定量依据.      

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

Indirect evaluation of Mars Gravity Model 2011 using a replication experiment on Earth

Curtin University??s Mars Gravity Model 2011 (MGM2011) is a high-resolution composite set of gravity field functionals that uses topography-implied gravity effects at medium- and short-scales (??125 km to ??3 km) to augment the space-collected MRO110B2 gravity model. Ground-truth gravity observations that could be used for direct validation of MGM2011 are not available on Mars??s surface. To indirectly evaluate MGM2011 and its modelling principles, an as-close-as-possible replication of the MGM2011 modelling approach was performed on Earth as the planetary body with most detailed gravity field knowledge available. Comparisons among six ground-truth data sets (gravity disturbances, quasigeoid undulations and vertical deflections) and the MGM2011-replication over Europe and North America show unanimously that topography-implied gravity information improves upon space-collected gravity models over areas with rugged terrain. The improvements are ??55% and ??67% for gravity disturbances, ??12% and ??47% for quasigeoid undulations, and ??30% to ??50% for vertical deflections. Given that the correlation between space-collected gravity and topography is higher for Mars than Earth at spatial scales of a few 100 km, topography-implied gravity effects are more dominant on Mars. It is therefore reasonable to infer that the MGM2011 modelling approach is suitable, offering an improvement over space-collected Martian gravity field models.     

Determination of gravimetric parameters of the gravity pole tide using observations recorded with superconducting gravimeters

In this study, the loading gravity effect of air mass changes calculated with the three-dimension (3D) meteorological data from the European Centre for Medium-range Weather Forecasts (ECMWF) are removed from superconducting gravimeter (SG) observations. The global hydrological gravity effect is computed and removed with hydrological data from the Global Land Data Assimilation System (GLDAS). Otherwise, the gravity influences induced by a theoretical self-consistent ocean pole tide and variations in length of day (LOD) are considered in the calculation. After removing the influences mentioned previously and also considering the long term trend in the data, a very nice linear relationship between the theoretical gravity pole tide and observed gravity residual (containing the observed gravity pole tide) for each of the selected 9 GGP stations we considered can be obtained. Therefore, the gravimetric factor of the gravity pole tide can be estimated with a simple linear regression. The results show that no clear phase lag is found between the theoretical gravity pole tide and observed gravity residuals from the nine SG stations.     

基于流动重力资料的平原地区深部构造解释

以河北省流动重力测点覆盖区域内邢台地区4条断裂为例,将卫星布格重力异常与实测布格重力异常进行对比,发现流动重力资料可更为详细地反映测点周围构造的空间分布特征,为精细刻画华北平原地区深部构造背景提供可靠资料。采用向上延拓方法,对实测布格重力异常进行处理,获得不同上延高度的布格重力异常图像,结合当地地震地质与构造资料进行分析,并对流动重力资料在构造解释上的应用前景进行展望,认为该方法可提高现有资料的利用价值,适合在流动重力测量地区推广。     

基于变差函数的插值方法计算华北地区重力场变化

定点重复重力测量是获取区域重力场变化的主要手段之一, 重力场特征与地形起伏、 构造走向等因素相关。 以华北地区为例, 考虑区域内地形和构造的北东向分布规律, 从EGM2008重力模型中拟合各向异性变差函数参数, 利用变差函数网格化插值, 对华北地区2009—2013年期间重力场观测数据进行网格化重建, 获取华北地区重力场时空变化结果。 研究结果表明基于区域地形、 构造特征的各向异性变差函数插值方法, 获得的空间重力场变化在重力异常梯级带上更加明显, 重力变化与活动构造分布具有更好的一致性。 本文研究方法对于恢复区域时空重力场异常具有重要意义, 有助于提高应用重力资料划分潜在地震危险区空间位置的精度, 为华北地区震情研判和构建地震预报定量指标体系提供可靠的地球物理场数据。     

地震重力观测技术的新发展方向:绝对重力与重力梯度一体化观测

高精度绝对重力与重力梯度一体化和阵列式观测可有效提高不同深度重力场源体的识别精度。本文基于中国地震局发布的《中国地球物理站网（重力）规划（2020—2030）》中绝对重力观测技术的不足和发展方向,提出了绝对重力与重力梯度一体化的观测技术思路。在监测中发现,强震孕育过程中壳—幔热物质运移引起的微重力变化信号,有利于提升破坏性地震的中、短期和临震预测能力。有必要进一步拓展重力站网绝对重力控制能力等方面的应用,为我国专业地震台站探索重力观测技术新发展方向。      

山西地震带不同基准的重力场变化

对2011—2017年山西及周边的流动重力观测资料,在绝对点控制下重新做平差计算,从不同尺度分析重力场的动态变化和跨断裂重力段差变化,重点对不同基准的重力场累积变化特征进行比较,结果显示:①文中的平均值基准可以明显削弱单期数据或个别测点对重力场造成的非正常影响;②在山西地震带持续出现了区域性的重力异常,形成了与构造断裂较为一致的重力梯度带;③山西地震带北部测段重力段差的变化幅值最大,南部次之但强于中部,反映了北部的构造活动比较活跃。     

Roughness of three types of gravity disturbances and their correlation with topography in rugged mountains and flat regions

We investigate the roughness of and the correlation with topography of the observed, topographically corrected (T), and bathymetrically and topographically corrected (BT) gravity disturbances. The numerical investigation is carried out for the gravity disturbances at the Earth’s surface and for the upward continued gravity disturbances at different altitudes. The area of study comprises a rough part of the Canadian Rockies surrounded by flat regions. The smoothest at the Earth’s surface are the BT gravity disturbances. The evolution of roughness with altitude shows an interesting phenomenon, diverse for the three types of gravity disturbances. The correlation with topography over the study area of the observed gravity disturbances is bellow 0.6, and of the BT gravity disturbances approximately −0.6. The largest absolute value, of about −0.75, is found between the topography and the T gravity disturbances. This large negative correlation indicates a presence of the isostatic compensation in mountainous regions of the Canadian west coast.     

云南省定点重力异常短临指标体系研究

本对云南内昆明,弥渡,下关等重力台站的重力固体潮观测资料进行了系统的研究。阐述了重力异常的提取判据,统计下重力异常与５级以上地震的历史对应情况,总结了重力异常变化特征,给出了重力异常短临指标体系,从而为今后的地震预报在重力学科提供了一些异常指标的定量依据。     

用重力测量技术观测城市地表下沉的实验研究

本文用重力测量技术对城市地表下沉进行了实验研究,从2016年3月到2017年5月在武汉市内地表下沉较大的部分城区进行了7期流动重力观测实验,并用D-InSAR观测的垂向位移进行了验证.数值结果表明重力观测每期整网平差后点值平均精度都小于10×10~(-8)m·s~(-2),说明采用重力观测能在城市内获得高精度的区域重力变化.第7期相对于第1期的结果与D-InSAR在大致相同时间段内地表垂直位移结果比较表明,重力增加的大部分区域与D-InSAR观测到的地表下沉区域相一致,说明这些区域的重力增加主要是由地表下沉引起的.从第2到7期相对于第1期的重力变化说明在近12个月的时间内测区最大重力变化约40×10~(-8)m·s~(-2),且局部区域的重力值是逐渐增加的,说明地表下沉是持续进行的.本实验结果说明重力观测技术能为城市地表下沉提供重力观测约束和机制解释.     

Regional correction of gravity data

Summary Methods for the elimination of the regional effect from gravity data are discussed in this paper. A historical review of the average gradient method used in torsion balance data is given. The connection between the averaging, or grid method, used on the gravity meter data and the average gradient method is pointed out and a criticism of these methods is given. Due to the ambiguity of the interpretation of gravity data no unique determination of the regional effect is possible and no mechanical method for its removal can be found. The determination of a regional gravity effect is an interpretation problem and for this reason it is unavoidably the subject of geological consideration and of the personal judgement of the interpreter. The method of smoothing the isogams and gravity curves is also discussed. A few suggestions are made for breaking the gravity data into two components (regional and residual) by the analysis of the gravity data. Finally the application of the second derivative methods for the elimination of the regional effect is discussed.     

2016年射阳M_S 4.4地震前后重力场变化

利用江苏重力测网2014—2017年重力场观测资料,采用绝对重力控制与相对重力联测相结合的平差方法,获取2016年射阳M_S 4.4地震前后重力场变化图。根据射阳M_S 4.4地震前后射阳地区各测线重力段差变化特点,绘制重力场等值线并进行对比分析,结合相关机理,探讨重力场变化与该地震的内在联系。同震观测数据显示:射阳M_S 4.4地震发生在重力异常值高梯度带附近,发震时震中地区位于NS挤压正异常、EW张拉负异常状态,震后区域重力梯度变化量开始减小,是一种典型的重力异常调整现象。     

华南地区时变重力场建模实验和异常分析

陆面时变重力测量是监测地壳内部密度变化和物质运移的重要手段.为确定华南时变重力观测网络的场源监测能力和重力场变化特征,本文采用球面六面体单元构建重力场模型,开展重力场建模实验,对比不同建模方法与噪声条件下的局部重力场恢复效果,并对2015-2017年来5期实测流动重力观测数据进行计算和分析.结果表明,最小二乘配置方法的...     

孕震过程中重力位、重力、重力梯度变化的数值模拟研究

简要介绍了为探索和解释地震前后重力变化的各种孕震模式,推导了由孕震引起的密度变化和位移与地面重力位、重力、重力梯度之间的变化关系。摸拟计算了孕震位移和密度变化引起的重力位、重力、重力梯度变化的空间分布并分析了重力位、重力、重力梯度变化的空间分布特征。与此同时,采用广义司托克斯云积分和有限差分方法对云南丽江7. 0级地震前重力位、重力、重力水平梯度进行了计算。结果表明,强震前重力位、重力、重力梯度有其自身的变化特征,这对预测强震有实际意义。     

Error sources and data limitations for the prediction of surface gravity: a case study using benchmarks

Gravity-based heights require gravity values at levelled benchmarks (BMs), which sometimes have to be predicted from surrounding observations. We use the Earth Gravitational Model 2008 (EGM2008) and the Australian National Gravity Database (ANGD) as examples of model and terrestrial observed data respectively to predict gravity at Australian National Levelling Network (ANLN) BMs. The aim is to quantify errors that may propagate into the predicted BM gravity values and then into gravimetric height corrections (HCs). Our results indicate that an approximate ±1 arc-min horizontal position error of the BMs causes maximum errors in EGM2008 BM gravity of ~22 mGal (~55 mm in the HC at ~2200 m elevation) and ~18 mGal for ANGD BM gravity because the values are not computed at the true location of the BM. We use RTM (residual terrain modelling) techniques to show that ~50% of EGM2008 BM gravity error in a moderately mountainous region can be accounted for by signal omission. Non-representative sampling of ANGD gravity in this region may cause errors of up to 50 mGals (~120 mm for the Helmert orthometric correction at ~2200 m elevation). For modelled gravity at BMs to be viable, levelling networks need horizontal BM positions accurate to a few metres, while RTM techniques can be used to reduce signal omission error. Unrepresentative gravity sampling in mountains can be remedied by denser and more representative re-surveys, and/or gravity can be forward modelled into regions of sparser gravity.     

Spectral Combination of Spherical Gradiometric Boundary-Value Problems: A Theoretical Study

The Earth’s gravity potential can be determined from its second-order partial derivatives using the spherical gradiometric boundary-value problems which have three integral solutions. The problem of merging these solutions by spectral combination is the main subject of this paper. Integral estimators of biased- and unbiased-types are presented for recovering the disturbing gravity potential from gravity gradients. It is shown that only kernels of the biased-type integral estimators are suitable for simultaneous downward continuation and combination of gravity gradients. Numerical results show insignificant practical difference between the biased and unbiased estimators at sea level and the contribution of far-zone gravity gradients remains significant for integration. These contributions depend on the noise level of the gravity gradients at higher levels than sea. In the cases of combining the gravity gradients, contaminated with Gaussian noise, at sea and 250?km levels the errors of the estimated geoid heights are about 10 and 3 times smaller than those obtained by each integral.     

Variations of the earth's gravity field due to the free nutation

Summary Expressions for the time variation of gravity and for deformations of equipotential surfaces due to the free nutation of the Earth have been derived. The possibility has been shown of determining the polar motion by means of accurate gravity observations under the assumption that the other effects, particularly tidal, may be expressed with sufficient accuracy theoretically. The necessity of introducing gravity corrections due to the polar motion has been pointed out in investigating the secular variations of gravity and in accurate gravity measurements (standard error of the order of 0.01 mgal and less), in general.     

2016年射阳MS 4.4地震前后重力场变化

利用江苏重力测网2014-2017年重力场观测资料,采用绝对重力控制与相对重力联测相结合的平差方法,获取2016年射阳M_S 4.4地震前后重力场变化图。根据射阳M_S 4.4地震前后射阳地区各测线重力段差变化特点,绘制重力场等值线并进行对比分析,结合相关机理,探讨重力场变化与该地震的内在联系。同震观测数据显示：射阳M_S 4.4地震发生在重力异常值高梯度带附近,发震时震中地区位于NS挤压正异常、EW张拉负异常状态,震后区域重力梯度变化量开始减小,是一种典型的重力异常调整现象。     

连续重力观测潮汐与非潮汐信号提取与处理平台设计与实现

介绍了连续重力观测潮汐与非潮汐信号提取与处理平台的平台架构、数据库、主要功能模块、技术特点与功能特性等。该平台面向中国地震重力站网的大网运维业务需求和海量数据处理需求,以连续重力数据处理中的潮汐信号处理、非潮汐信号处理、时频变换等科学计算为核心,基于大数据技术实现了中国地震重力站网海量数据的自动清理、数据质量自动评估、数据产品自动产出,产出重力潮汐因子、时频图等24种数据产品,提升了连续重力观测数据的应用效能和地震重力站网的运行、管理和服务水平。