贵州中东部剩余重力异常与区域矿产分布关系

重点研究1:20万剩余重力异常,结合航磁、物性、地质等资料,初步圈定区内隐伏、半隐伏岩浆岩体,对研究与重力异常、岩体分布关系密切的矿产的时间、空间分布及其成因,有一定参考价值和指示意义。     

一个新的沉积变质型钛矿床的成矿模式及找矿方法

内蒙古羊蹄子山沉积变质型钛矿床是一个新的钛矿床类型。通过对地质、地球物理前提条件的研究和找矿方法试验,认为高精度磁法是找沉积变质型钛矿床的基本方法,重力为有效的找矿方法。在矿区北带磨石山以金红石为主的富钛矿体上有600 μGal的重力异常,无磁异常;南带的羊蹄子山钛铁矿体上重力异常为500 μGal,有60 nT的弱磁异常。     

相控建模技术在羊二庄油田的应用

羊二庄油田的明馆储层非均质性严重,高含水、高开发后期的剩余油分布比较复杂,第一次油藏描述中应用传统的确定性建模方法已经不能满足剩余油挖潜的需要.采用RMS相控建模软件,通过构造建模、沉积微相建模、储层物性建模,建立羊二庄油田三维可视化相控地质模型,并实现了与数模软件的无缝接口.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.     

Using Topological Relationships to Inform a Data Integration Process

When spatial datasets are overlaid, corresponding features do not always coincide. This may be a result of the datasets having differing quality characteristics, being captured at different scales or perhaps being in different projections or datums. Data integration methods have been developed to bring such datasets into alignment. Although these methods attempt to maintain topological relationships within each dataset, spatial relationships between features in different datasets are generally not considered. The preservation of inter‐dataset topology is a research area of considerable current interest. This research addresses the preservation of topology within a data integration process. It describes the functional models established to represent a number of spatial relationships as observation equations. These are used to provide additional information concerning the relative positions of features. Since many topological relationships are best modelled as inequalities, an algorithm is developed to accommodate such relationships. The method, based on least squares with inequalities (LSI), is tested on simulated and real datasets. Results are presented to illustrate the optimal positioning solutions determined using all of the available information. In addition, updated quality parameters are provided at the level of the individual coordinate, enabling communication of local variation in the resultant quality of the integrated datasets.     

利用参考重力场模型基于能量法确定GRACE加速度计校准参数

利用多个参考重力场模型分别对GRACE一个月的实测加速度计观测数据进行检校.数值计算结果的比较分析表明了利用参考重力场模型确定加速度计校准参数是有效的.     

附加约束条件的序贯平差

从观测量分组的参数平差出发,详细推导了附加约束条件的序贯平差及其精度估计,讨论了附加约束条件的抗差序贯平差的计算过程.算例表明,推导的公式正确有效,简单实用,附加约束条件的序贯平差抗差估计能有效地抵制粗差的影响.     

邹平火山岩盆地铜矿地质成因地球物理特征探讨

邹平铜矿处于齐河-广饶深大断裂带南部的邹平火山岩盆地中,形成于破火山口火山通道充填的石英正长闪长岩岩颈中央上部,包括伟晶岩型铜矿和细脉浸染状斑铜矿床两种类型。前者矿体较小,但品位高;后者品位较低,但规模中等。含矿石英正长闪长岩等密度小、磁性弱,故在火山岩系中呈现高背景重力场上的重力低和杂乱高磁场背景中的低负异常,即“重磁同低”,且高极化。重磁同低异常区和高极化率异常带,是本区寻找铜矿的有利部位。