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.     

Depth from extreme points method for gravity gradient tensor data

The ‘depth from extreme points’ method is an important tool to estimate the depth of sources of gravity and magnetic data. In order to interpret gravity gradient tensor data conveniently, formulas for the tensor data form regarding depth from the extreme points method were calculated in this paper. Then, all of the gradient tensor components were directly used to interpret the causative source. Beyond the g_zz component, also the g_xx and g_yy components can be used to obtain depth information. In addition, the total horizontal derivative of the depth from extreme points of the gradient tensor can be used to describe the edge information of geologic sources. In this paper, we investigated the consistency of the homogeneity degree calculated by using the different components, which leads to the calculated depth being confirmed. Therefore, a more integrated interpretation can be obtained by using the gradient tensor components. Different synthetic models were used with and without noise to test the new approach, showing stability, accuracy and speed. The proposed method proved to be a useful tool for gradient tensor data interpretation. Finally, the proposed method was applied to full tensor gradient data acquired over the Vinton Salt Dome, Louisiana, USA, and the results are in agreement with those obtained in previous research studies.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.

     

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

传统动力学法的观测方程以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模型.     

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

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

地形/均衡重力场模型构制及其拟稳分析

重力空白区数据填补的一个主要方法是基于地壳均衡理论进行的,该方法亦用于EGM系列模型的构建中.本文研究了地形数据在构制地形/均衡重力场模型中的应用,分析了补偿深度对Airy位模型和面凝聚位模型的影响,给出二者的最佳补偿深度分别为50 km和40 km.以纯卫星重力模型为参考,后者在前120阶的精度要高于前者,但在121~250阶的精度较低,组合模型精度高于单一模型精度.对地形/均衡地球重力场模型进行了EGM2008拟稳分析,研究了不同分辨率基准的拟稳效果,分析表明：30'分辨率的拟稳基准所得拟稳模型对应的阶方差与参考阶方差曲线直到360阶都有较好的一致性,以EGM2008为基准,其相对累计大地水准面高误差在140阶时为6.83cm,相对累计重力异常误差在220阶时为1.10 mGal.     

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.     

Main crustal discontinuities of Morocco derived from gravity data

Sharp linear gradients in maps of potential field data are generally assumed to result from sharp discontinuities or boundaries between rocks having different densities or magnetic susceptibilities and are usually associated with faults or other geological contacts. The computation of the horizontal gradients of the gravity field permits us to localize the limits of such blocks and then the fault locations. The horizontal derivative maxima of the Bouguer anomaly and its upward continuation at several heights show lineaments that could reflect the layout of faults and/or contacts and their dip directions. The application of this method to the Bouguer anomaly map of Morocco (with 19,571 points, using an average crustal density ρ = 2.67 g/cm³) allowed us to perform a multiscale analysis of the gravimetric lineaments of the country. The obtained structural map is consistent with several faults already identified in previous studies, and highlights five new major subsurface faults systems with location and dip: the Saghro fault system; Bou-Arfa Midelt fault system; Sidi Slimane Mezquitem fault; Ksar El Kebir–Chefchaouen fault and the Rifan West Mediterranean fault. In addition, this study suggests a new shape and localization for the Agadir-Oujda trans-Moroccan major fault with a NE-SW direction and 900 km length, subdividing Morocco into two main domains. The results of this study contribute to the improvement of the regional structural map of the north western part of Africa, which is situated within the convergence zone between Africa and Eurasia.     

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.     

三维重力、重力梯度和大地电磁数据的平滑聚焦结构约束联合反演

多种地球物理数据联合反演解释能够显著提高勘探分辨率,近年来,基于交叉梯度结构约束的联合反演得到了广泛研究.该方法有效克服了单一方法反演易受多解性影响和不同方法反演结果结构不一致等问题,已成为目前最流行的联合反演技术之一.传统的重力和大地电磁联合反演常因重力纵向分辨能力弱,导致密度模型可靠性较低,进而使联合反演结果有失准...