一种卫星重力梯度数据的ARMA最优滤波模型构建方法

由于卫星重力梯度观测的有色噪声特性和海量观测特征,在利用直接法进行重力场模型的最小二乘求解时,观测值的协方差阵为超大型的非对角阵,这给数值求解带来了极大困难.本文提出了一种基于先验误差功率谱密度的最优ARMA滤波模型构建方法,结合法方程的分块求解策略,可实现对卫星重力梯度观测值的高效滤波处理.数值仿真结果表明,利用最优ARMA滤波器进行时域滤波后,法方程的态性得到了明显改善,重力梯度观测值中的有色噪声得到了有效的"白化"处理,大地水准面精度得到了显著提升.     

Smoothing impact of isostatic crustal thickness models on local integral inversion of satellite gravity gradiometry data

The effects of topographic masses on satellite gradiometric data are large and in order to reduce the magnitude of these effects some compensation mechanisms should be considered. Here we use the isostatic hypotheses of Airy-Heiskanen and the recent Vening Meinesz-Moritz for compensating these effects and to smooth the data prior to their downward continuation to gravity anomaly. The second-order partial derivatives of extended Stokes’ formula are used for the continuations over a topographically rough territory like Persia. The inversions are performed and compared based on two schemes of the remove-compute-restore technique and direct downward continuation. Numerical results show that the topographic-isostatic effect based on Vening Meinesz-Mortiz’s hypothesis smoothes the data better than that based on Airy-Heiskanen’s hypothesis. Also the quality of inversions of the smoothed data by this mechanism is twice better than that of the nonsmoothed ones.     

重力和重力梯度数据三维相关成像

本文提出了重力异常三维相关成像方法和重力梯度数据三维相关成像方法,并提出了基于异常分离的三维相关成像方法来提高成像分辨率.通过合成Y型岩脉模型和合成多个直立长方体组合模型的重力异常和重力梯度数据试验分析,验证了本文三维相关成像方法可显示出异常地质体的空间赋存状态和等效剩余质量分布,具有良好的纵向和横向分辨率.     

卫星重力梯度观测数据L1级构建方法

高精度重力梯度观测数据L1级构建的系统方法是推进我国自主重力卫星任务重要的基础数据处理技术.本文以GOCE卫星L1级数据预处理技术和关键载荷原始数据为参考,面向我国发展的梯度测量卫星的任务需要,系统研究并初步实现了卫星重力梯度观测数据L1级构建方法,主要包括加速度计电压数据转换、多星敏感器联合姿态数据的角速度重建、卫星重力梯度分量构建等技术内容.计算结果表明,加速度计超灵敏轴精度为10^-10~10^-11m·s^-2·Hz^-1/2,达到重力梯度仪设计精度要求;多星敏感器联合解算最佳姿态角速度w_y、w_z在10~100 mHz内精度约提升1个量级,其精度约达到10^-5 rad·s^-1·Hz^-1/2量级,能够有效抑制低精度角速度分量在坐标系转换中导致的噪声传播;基于维纳滤波方法恢复的角速度在5~100 mHz频段内的平方根功率谱密度提升了（5.21~6.56）×10^-11 rad·s^-1·Hz^-1/2,显示了基于高精度角速度解算重力梯度分量的必要性;构建重力梯度各分量计算值与全球重力场和海洋环流探测器（GOCE）官方公布的重力梯度分量精度相当,其梯度张量的迹在20~100 mHz频段范围内约为10 mE·Hz^-1/2,验证了本文构建方法的有效性.研究工作可为下一步我国推进实施民用重力梯度测量卫星任务提供自主的原始数据处理技术支撑与储备.

     

基于地面重力的卫星重力梯度检校方法

卫星重力梯度仪在轨检校是提高梯度模式重力卫星观测质量的关键.本文面向中国未来梯度模式重力卫星规划任务,研究提出一种基于地面重力的卫星在轨检校方法,该方法顾及卫星设计指标,从地面先验重力数据的精度、空间尺度以及卫星重力梯度仪的观测噪声等卫星检校要素开展分析研究,成功实现地面数据与卫星观测数据的比对检校.研究结果表明,在1...     

Heli-borne gravity gradiometry in rugged terrain

For airborne gravity gradiometry in rugged terrain, helicopters offer a significant advantage over fixed-wing aircraft: their ability to maintain much lower ground clearances. Crucially, this provides both better signal-to-noise and better spatial resolution than is possible with a fixed-wing survey in the same terrain. Comparing surveys over gentle terrain at Margaret Lake, Canada, and over rugged terrain at Mount Aso, Japan, demonstrates that there is some loss of spatial resolution in the more rugged terrain. The slightly higher altitudes forced by rugged terrain make the requirements for terrain correction easier than for gentle terrain. Transforming the curvature gradients measured by the Falcon gravity gradiometer into gravity and the complete set of tensor components is done by a Fourier method over gentle terrain and an equivalent source method for rugged terrain. The Fourier method is perfectly stable and uses iterative padding to improve the accuracy of the longer wavelengths. The equivalent source method relies on a smooth model inversion, and the source distribution must be designed to suit the survey design.     

Comparison of satellite and terrestrial gravity data

Summary The integral mean values of gravity on the surface W=W ₀ , obtained from satellite observations with the use of harmonic coefficients[3, 7] and from terrestrial gravity measurements[12], are compared. The squares and products of the harmonic coefficients were neglected, with the exception of [J ₂ ⁽⁰⁾ ] ² , which was taken into account. The Potsdam correction and the geocentric constant are being discussed. The paper ties up with[13–15] and the symbols used are the same. The given problem was treated, e.g., in[2, 4, 6, 8–10]; in the present paper the values of gravity are compared directly.     

Satellite gravity gradiometry: Secular gravity field change over polar regions

The ESA Gravity and steady state Ocean and Circulation Explorer, GOCE, mission will utilise the principle of satellite gravity gradiometry to measure the long to medium wavelengths in the static gravity field. Previous studies have demonstrated the low sensitivity of GOCE to ocean tides and to temporal gravity field variations at the seasonal scale. In this study we investigate the sensitivity of satellite gradiometry missions such as GOCE to secular signals due to ice-mass change observed in Greenland and Antarctica. We show that unaccounted ice-mass change signal is likely to increase GOCE-related noise but that the expected present-day polar ice-mass change is below the GOCE sensitivity for an 18-month mission. Furthermore, 2–3 orders of magnitude improvement in the gradiometry in future gradiometer missions is necessary to detect ice-mass change with sufficient accuracy at the spatial resolution of interest.     

Focusing iterative migration of gravity gradiometry data acquired in the Nordkapp Basin,Barents Sea

Geological interpretation based on gravity gradiometry data constitutes a very challenging problem. Rigorous 3D inversion is the main technique used in quantitative interpretation of the gravity gradiometry data. An alternative approach to the quantitative analysis of the gravity gradiometry data is based on 3D smooth potential field migration. This rapid imaging approach, however, has the shortcomings of providing smooth images since it is based on direct integral transformation of the observed gravity tensor data. Another limitation of migration transformation is related to the fact that, in a general case, the gravity data generated by the migration image do not fit the observed data well. In this paper, we describe a new approach to rapid imaging that allows us to produce the density distribution which adequately describes the observed data and, at the same time, images the structures with anomalous densities having sharp boundaries. This approach is based on the basic theory of potential field migration with a focusing stabilizer in the framework of regularized scheme, which iteratively transfers the observed gravity tensor field into an image of the density distribution in the subsurface formations. The results of gravity migration can also be considered as an a priori model for conventional inversion subsequently. We demonstrate the practical application of migration imaging using both synthetic and real gravity gradiometry data sets acquired for the Nordkapp Basin in the Barents Sea.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料,包括K波段、星载加速度和卫星轨道数据,反演了80阶地球重力场模型IGGGRACE01S,该模型在半波长为500km的空间分辨率上,确定大地水准面的精度约为0012m,中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比,该模型系数最接近于EIGEN_GRACE02S,与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏,结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近,与EGM96模型结果差异较大,差别较大地区主要在南极等地区. 对于中国大陆,比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格）,两者之间的标准偏差为48mGal.     

Determination of gravity anomaly at sea level from inversion of satellite gravity gradiometric data

Gravity gradients can be used to determine the local gravity field of the Earth. This paper investigates downward continuation of all elements of the disturbing gravitational tensor at satellite level using the second-order partial derivatives of the extended Stokes formula in the local-north oriented frame to determine the gravity anomaly at sea level. It considers the inversion of each gradient separately as well as their joint inversion. Numerical studies show that the gradients T_zz, T_xx, T_yy and T_xz have similar capability of being continued downward to sea level in the presence of white noise, while the gradient T_yz is considerably worse than the others. The bias-corrected joint inversion process shows the possibility of recovering the gravity anomaly with 1 mGal accuracy. Variance component estimation is also tested to update the observation weights in the joint inversion.     

Atmospheric effects in the derivation of geoid-generated gravity anomalies

Parameters of the gravity field harmonics outside the geoid are sought in solving the Stokes boundary-value problem while harmonics outside the Earth in solving the Molodensky boundary-value problem. The gravitational field generated by the atmosphere is subtracted from the Earth’s gravity field in solving either the Stokes or Molodensky problem. The computation of the atmospheric effect on the ground gravity anomaly is of a particular interest in this study. In this paper in particular the effect of atmospheric masses is discussed for the Stokes problem. In this case the effect comprises two components, specifically the direct and secondary indirect atmospheric effects. The numerical investigation is conducted at the territory of Canada. Numerical results reveal that the complete effect of atmosphere on the ground gravity anomaly varies between 1.75 and 1.81 mGal. The error propagation indicates that precise determination of the atmospheric effect on the gravity anomaly depends mainly on the accuracy of the atmospheric mass density distribution model used for the computation.     

On integral approach to regional gravity field modelling from satellite gradiometric data

Solution of the gradiometric boundary value problems leads to three integral formulas. If we are satisfied with obtaining a smooth solution for the Earth’s gravity field, we can use the formulas in regional gravity field modelling. In such a case, satellite gradiometric data are integrated on a sphere at satellite level and continued downward to the disturbing potential (geoid) at sea level simultaneously. This paper investigates the gravity field modelling from a full tensor of gravity at satellite level. It studies the truncation bias of the integrals as well as the filtering of noise of data. Numerical studies show that by integrating T _zz with 1 mE noise and in a cap size of 7°, the geoid can be recovered with an error of 12 cm after the filtering process. Similarly, the errors of the recovered geoids from T _xz,yz and T _{xx-yy, 2xy} are 13 and 21 cm, respectively.     

Low-degree gravity change from GPS data of COSMIC and GRACE satellite missions

This paper demonstrates estimation of time-varying gravity harmonic coefficients from GPS data of COSMIC and GRACE satellite missions. The kinematic orbits of COSMIC and GRACE are determined to the cm-level accuracy. The NASA Goddard's GEODYN II software is used to model the orbit dynamics of COSMIC and GRACE, including the effect of a static gravity field. The surface forces are estimated per one orbital period. Residual orbits generated from kinematic and reference orbits serve as observables to determine the harmonic coefficients in the weighted-constraint least-squares. The monthly COSMIC and GRACE GPS data from September 2006 to December 2007 (16 months) are processed to estimate harmonic coefficients to degree 5. The geoid variations from the GPS and CSR RL04 (GRACE) solutions show consistent patterns over space and time, especially in regions of active hydrological changes. The monthly GPS-derived second zonal coefficient closely resembles the SLR-derived and CSR RL04 values, and third and fourth zonal coefficients resemble the CSR RL04 values.     

Some aspects of interpretation of tensor gradiometry data

Aspects of the interpretation of measured data on the gravity gradient tensor (GGT) are examined. The problem is posed in relation to the great progress achieved in recent years in the development of instrumentation and the method of GGT measurements on mobile carriers. In our opinion, the new methods of measurement and the new data obtained with their help require the development of new methods of interpretation of potential fields. The paper addresses two methods taking advantage of simultaneous measurements of all components of the GGT and anomalies of the gravitaty field V _z. It is shown that the joint analysis of all GGT components can provide independent constraints on the noise level in various components. The method of tensor deconvolution proposed in the paper is a tensor analogue of the Euler method. The method is based on the calculation of invariants and is, therefore, stable with respect to the orientation uncertainties of the measuring system. The method provides means for estimating the structural index and, therefore, is particularly effective in the treatment of fields that contain isometric and/or elongated anomalies. The calculation of invariants and the tensor ratio can also be used for the development of procedures enabling automatic estimation of the axis strike azimuths of elongated anomaly-forming bodies.     

航空重力梯度测量之地形改正方法及讨论

本文简要介绍了航空重力梯度测量中地形改正的基本原理及其地改公式,并采用基于质量线和质量棱柱两种地形模型的数值积分公式进行了一些简单模拟试算.在研究航空重力梯度测量的精密地形改正方面,作者首次将考虑地球曲率影响的质量棱柱体模型用于梯度数据的地改.文章最后对所有结果做了一些必要统计分析.     

Locating anomalies by gravity gradiometry using the matched filter with a probabilistic assessment

The matched filter can be used to identify the location of a specific signal embedded in background correlated noise. The maximum filter output indicates the likely location and, with appropriate statistical assumptions on the noise, it also serves as a test statistic in the probabilistic evaluation of the filter’s performance. Different setups of the Neyman-Pearson statistical hypothesis test yield predictions of either the probability of a miss or that of a false alarm. The needed statistics of the maximum filter output are properly obtained using the distribution of order statistics. Through Monte Carlo simulations, we analyzed the ability of the matched filter to identify a sub-surface anomaly in typically correlated gravity fields using observations of elements of the gravity gradient tensor. We also evaluated the reliability of the hypothesis testing and the associated predicted probabilities of misses and false alarms. Our simulations and statistical analyses confirm that the power of the tests increases as the signal strength increases and as more gradient tensor components per observation point are included. We also found that the hypothesis test that is designed to predict the probability of a miss is more robust and powerful than the one for predicting a false alarm. Moreover, the probability of a miss is somewhat smaller than the probability of a false alarm under otherwise equal circumstances.     

基于卫星重磁资料的南海深部构造研究

南海深部构造对研究南海构造演化和油气勘探具有重要意义.本文对南海地区的自由空气重力异常进行布格校正、海水层校正和沉积层校正,得到布格重力异常,再对布格重力异常进行区域异常和局部异常分离,利用位场界面反演方法对区域布格异常进行反演计算得到研究区域的莫霍面深度分布;采用全变倾角化极方法对研究区域的卫星磁异常数据进行化极处理,并进一步对化极磁异常作向上延拓,得到延拓后化极磁异常结果.分析布格重力异常、莫霍面深度及化极磁异常特征,结合天然地震层析成像的证据,得到以下结论：推测南海北部陆缘的古俯冲带位置是从118.5°E,24°N沿北东向延伸至109°E,15°N;红河断裂入海后经过莺歌海盆地在海南岛南部转为南北向与越东断裂相接并延伸至万安盆地;推测中特提斯洋的部分闭合位置是从110°E,2°N到101°E,21°N.

     

利用GRACE卫星重力数据监测关中地区地下水储量变化

关中地区作为一带一路重要的工农业发达地区之一,开展针对该地区地下水储量变化的监测和分析工作对揭示地下水储量变化特征与经济社会发展具有重要现实意义.本文基于2003—2014年GRACE卫星重力场模型数据,采用组合滤波及单一尺度因子方法反演了关中地区陆地水储量变化,扣除GLDAS地表水平均结果,对关中地区地下水储量变化进行了监测分析.将陆地水储量变化与GLDAS进行相关性分析,将地下水储量变化与WGHM地下水模型及实测地下水位结果进行对比分析.研究结果表明：①关中地区陆地水变化与GLDAS模型结果具有较强的相关性,相关系数多数大于0.7,其中与模型平均结果的相关系数可达0.8.② 2003—2008年关中地区地下水呈正增长趋势,增加速率为0.25 cm·a^-1,与同期实测数据变化趋势一致;但2003—2013年地下水存在长期亏损,亏损速率为-0.37 cm·a^-1等效水高,这与同时期WGHM估算结果-0.35 cm·a^-1十分吻合.③关中地区地下水存在明显的年变化特征,在2003—2014年期间地下水减少速率为-0.44 cm·a^-1,与该地区降雨量有较好的一致性,在降雨偏少的2008、2012和2013年,地下水也显著减少.

     

星间加速度解算模式的模拟与分析

低低卫星跟踪卫星的观测量是两低轨卫星的星间距离或星间速度,星间加速度由星间速度通过数值微分导出,用星间加速度作为观测量可以避免解算卫星运动的变分方程,简化观测方程的建立,但数值微分会使观测噪声放大,从而影响重力位的解算精度.为了定量给出星间加速度观测模式的精度,本文分析并模拟验证了数值微分公式计算星间加速度的精度,导出了基于星间加速度的一般形式的观测方程,模拟计算了基于星间加速度的重力位模型.结果表明,采用星间加速度观测模式的解算精度要明显低于星间速度观测模式的解算精度.