Spectral analysis for recovering the Earth's gravity potential by satellite gravity gradients

根据仪器功率谱密度和重力位系数阶方差的定义,本文建立了卫星重力梯度测量噪声功率谱密度与重力场模型的误差阶方差的直接对应关系,并基于此讨论了重力梯度测量精度、卫星轨道高度以及运行时间对地球重力场恢复精度的影响.相比于传统的基于最小二乘法评估卫星载荷噪声对地球重力场恢复精度的影响而言,本文提出的方法简单、直接,有助于快速设计和确定卫星重力测量计划的有关参数.     

A study of the upward continuation of gravity data from a plane surface

Summary A summation method of upward continuation of gravity data has been considered under the assumption that observations are available at regular intervals. The upward continued value has been obtained as the sum of products of individual gravity values and corresponding theoretical coefficients. Besides the usual parameter involving horizontal and vertical distances, the theoretical coefficients have been generalized to be dependent also on i) the order of a low order polynomial assumed to represent the gravity variation around a grid point and ii) the weights assigned to the gravity values at the nearest four grid points used for least-squares determination of the polynomial. Since the observations in practical cases are available over a finite area only, the effect of truncation of the area of the integration has also been discussed separately. The method has been programmed and tested on a three-dimensional model, whose true gravity effects were computed at various levels over a finite area. Upward continued values have been computed under various assumptions about the gravity field in the outside region. Comparisons of these results with the true values indicate that the truncation effect becomes increasingly important for larger values of the ratio of elevation to grid separation and/or when the gravity field is not a local one. It has also been found that the greater is the above ratio, the less important is the effect of weights on the theoretical coefficients and practically vanishes(<10 ^–4 ) when the ratio is greater than5.0.     

基于重力卫星几何轨道线性化的地球重力场反演方法

传统动力学法的观测方程以6个初始轨道参数和先验力模型为初值进行线性化,其线性化误差随积分弧长拉长而增大.本文直接以重力卫星的几何观测轨道为初值进行线性化,其线性化误差与轨道弧长无关,且不需要初始重力场模型和初始轨道参数.导出了基于卫星轨道观测值反演重力场模型的相关公式,利用JPL公布的RL02版本2008年全年的GRACE双星轨道数据和加速度计数据解算了90阶次的地球重力场模型TJGRACE01S,并以EGM2008模型为基准与其他模型进行了比较分析,结果表明:TJGRACE01S模型直到90阶次的大地水准面累积误差为17.6 cm,优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型,前27阶位系数整体精度优于EIGEN-GRACE01S,前15阶位系数整体精度与EIGEN-GRACE02S模型精度大致相当.利用美国8221个GPS水准点数据的分析结果也表明,本文模型也优于同阶次的EIGEN-CHAMP03S和EIGEN-CHAMP05S模型.     

不同高程基准位差计算的严密理论研究

高程基准是大地测量参考框架的重要组成部分,也是数字地球基础框架的重要内容.但目前各个国家或地区使用的大多是局部高程基准,这严重影响了地球信息数据的全球共享机制.本文利用高程基准与重力位之间的关系,推导了不同高程系统之间的严密位差计算公式和近似计算公式.该理论将为我国和全球高程基准的统一提供理论支持.     

A generalized method for various data processing techniques in gravity interpretation

Summary The various data processing techniques (downward continuation, first and second derivatives and their downward continuation) used in gravity interpretation, are analogous to different types of linear filtering operations whose theoretical filter (amplitude) responses can be derived from by suitably choosingN andd, whereu andv are angular frequencies in two perpendicular directions,d the height or depth of continuation in unit of grid interval; andN denotes the order of the vertical derivative. By incorporating a mathematical smoothing function, ( being the smoothing parameter) in the theoretical filter response function, it has been possible, by selecting a suitable value of smoothing parameter, to establish an approximate equivalence of the effect of the mathematical smoothing with the inherent smoothing introduced, because of the numerical approximation (approximation error) for practically all data-processing techniques. This approximate equivalence leads to a generalized method of computing sets of weight coefficients for various data-processing techniques from filter response matching method. Several sets of weight coefficients thus have been computed with different smoothing parameter. The amplitude response curves of the various existing sets of weight coefficients have also been calculated for assessing the quality of the approximation in achieving the desired filtering operation.     

On computation of second derivatives from gravity data

Summary Second derivative formulae ofElkins andRosenbach are developed from a unified approach and their merits and demerits are discussed. Some new formulae are suggested, which, from theoretical considerations it is felt, will prove more efficient.     

重力位场谱分析方法研究综述

利用Fourier变换将空间域位场数据变换到频率域进行谱分析是重力位场数据处理中的一种重要方法．这种方法在位场的延拓及转换、场分离、物性界面的正反演和位场的曲化平处理等方面已得到广泛的应用．本文详细介绍了位场谱分析的基本理论和方法,总结了国内外最近几十年来在这方面取得的主要成果,并提出了未来值得深入研究的方向．     

A geoid solution for airborne gravity data

Airborne gravity data is usually attached with satellite positioning of data points, which allow for the direct determination of the gravity disturbance at flight level. Assuming a suitable gridding of such data, Hotine’s modified integral formula can be combined with an Earth Gravity Model for the computation of the disturbing potential (T) at flight level. Based on T and the gravity disturbance data, we directly downward continue T to the geoid, and we present the final solution for the geoid height, including topographic corrections. It can be proved that the Taylor expansion of T converges if the flight level is at least twice the height of the topography, and the terrain potential will not contribute to the topographic correction. Hence, the simple topographic bias of the Bouguer shell yields the only topographic correction. Some numerical results demonstrate the technique used for downward continuation and topographic correction.     

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.     

Interpretation of tensor gravity data using an adaptive tilt angle method

Full Tensor Gravity Gradiometry (FTG) data are routinely used in exploration programmes to evaluate and explore geological complexities hosting hydrocarbon and mineral resources. FTG data are typically used to map a host structure and locate target responses of interest using a myriad of imaging techniques. Identified anomalies of interest are then examined using 2D and 3D forward and inverse modelling methods for depth estimation. However, such methods tend to be time consuming and reliant on an independent constraint for clarification. This paper presents a semi‐automatic method to interpret FTG data using an adaptive tilt angle approach. The present method uses only the three vertical tensor components of the FTG data (T_zx, T_zy and T_zz) with a scale value that is related to the nature of the source (point anomaly or linear anomaly). With this adaptation, it is possible to estimate the location and depth of simple buried gravity sources such as point masses, line masses and vertical and horizontal thin sheets, provided that these sources exist in isolation and that the FTG data have been sufficiently filtered to minimize the influence of noise. Computation times are fast, producing plausible results of single solution depth estimates t hat relate directly to anomalies. For thick sheets, the method can resolve the thickness of these layers assuming the depth to the top is known from drilling or other independent geophysical data. We demonstrate the practical utility of the method using examples of FTG data acquired over the Vinton Salt Dome, Louisiana, USA and basalt flows in the Faeroe‐Shetland Basin, UK. A major benefit of the method is the ability to quickly construct depth maps. Such results are used to produce best estimate initial depth to source maps that can act as initial models for any detailed quantitative modelling exercises using 2D/3D forward/inverse modelling techniques.     

Dip angle determination of fault planes from gravity data

This paper illustrates an approach for determining uniquely the dip angle of a faulted structure, approximated by two semi-infinite horizontal slabs, displaced vertically from each other, by making use of the maximum positive and the maximum negative gravity amplitudes. The method is simple and rapid, and it does not necessarily require computing facilities.Two field examples from Aberdeenshire, Scotland and from the Northern Eastern Desert of Egypt are presented in which the angles of dip of two geological contacts are estimated at 45 and 40 degrees, respectively, and verified from surface geology, gravity modeling, and seismic data.     

川滇地区重力场与深部结构特征

川滇地区地处青藏高原东缘强烈变形区,区内的龙门山断裂带、安宁河断裂带与鲜水河断裂带构成了华南块体、川滇块体与松潘-甘孜块体的重要分界.为研究川滇地区深部地壳结构特征,本文利用布格重力数据,并应用小波多尺度分析方法对重力场进行了分离,进而对该区地壳介质密度结构进行了分析.研究结果表明:川滇块体与松潘-甘孜块体中上地壳密度...     

Application of relaxation method to upward continuation of gravity and magnetic data

Summary The method of relaxation has been applied to the upward continuation of gravity and magnetic data. The method yields results that are more accurate than those given byPeters' scheme, and has certain added advantages. In an alternative approach, the finite Fourier sine transform has been used to reduce an area relaxation to a series of line relaxations. The results, again, are very satisfactory. Only the two-dimensional problem is considered, but, with the help of a digital computor, the method is directly extensible to three-dimensional problems.Presented at the symposium on Geophysical Prospecting held by the Central Board of Geophysics at Baroda (India) between Aug. 15–17, 1959. Published by kind permission of Secretary, Central Board of Geophysics.     

一种改进的重力梯度全张量数据的边界识别方法

边界识别是重力资料解释中的一项重要任务.随着重力梯度测量技术的迅速发展,重力梯度张量数据在边界识别中的应用越来越广泛.本文重点研究了随着深度的增加,边界识别能力下降,正负异常中出现假边缘的问题.另外,有些边缘检测方法对走向不同的地质体识别能力有所差异.本文对基于重力梯度张量的水平方向Theta法进行改进,通过选择合适的...     

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.     

重力数据3D密度成像中EXTR方法的各参数变化对反演结果的影响(英文)

Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的Extrapolation Tikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。     

Quantitative interpretation of three dimensional finite mass distribution from gravity data

Summary It is well known that the problem of determination of subsurface masses from gravity observations does not offer any unique solution. Yet, as has been shown in this paper, under some simplified but realistic assumptions some parameters of the subterrainean body are uniquely determinable in terms of some weighted double summations of gravity values (residuals). Physical significances of these parameters are also discussed. They themselves suggest about the possible mass distribution, more over they can be utilised in «cut and try method» of the determining masses.Published with the kind permission of the Director, Geological Survey of India, Calcutta.     

GOCE引力梯度的频谱分析及滤波

GOCE卫星提供的梯度数据含有非常大的低频误差,如何处理这种误差是GOCE数据处理中最为关键的工作之一.本文根据GOCE卫星的运行情况,首先分析了梯度数据的频率特性,推导了频率与阶次的对应关系;并在此之上,介绍了针对低频误差的滤波方法,即移去恢复和向前向后滤波方法,前者可解决滤波中的低频信号损失问题,后者则主要解决了滤波中的相位漂移问题.最终结果表明:引力梯度的时间频谱与球谐展开中的阶次虽不是一一对应的,但各阶所对应的最大截止频率与阶次却有一定的显式表达.同时也表明,本文所采用的滤波方法是有效的,达到了消除低频误差但保留观测频段信号的目的.     

A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.