Inverse gravity problems for models composed of bodies of Sretenskii’s class

An inverse gravity problem is solved for a geological model consisting of bodies of Sretenskii’s class. The position of the middle plane is fixed for each body. It is required to determine the upper and lower boundaries of a body, which are described by analytical functions and are parameterized. The solution of the problem is illustrated by an example.     

Solving the structural inverse gravity problem by the modified gradient methods

New methods for solving the three-dimensional inverse gravity problem in the class of contact surfaces are described. Based on the approach previously suggested by the authors, new algorithms are developed. Application of these algorithms significantly reduces the number of the iterations and computing time compared to the previous ones. The algorithms have been numerically implemented on the multicore processor. The example of solving the structural inverse gravity problem for a model of four-layer medium (with the use of gravity field measurements) is constructed.     

利用重力资料反演三维密度界面的质面系数法

文中给出一深度为h的水平矩形质面在地面产生的重力公式,对比了此公式算出的重力异常与三维长方体公式算出的重力异常,得出:当三维地质体的厚度与其埋深之比小于0.2-0.4时,二者的重力效应几乎相等,它们之间最大偏差在1-3％以内。用给出的重力公式进一步导出了重力反演三维密度界面的质面系数法。用理想模式检验,结果表明,本文给出的方法收敛快,反演误差较小,计算效果较好。最后用一实例,验证了方法的适用性。     

On a method of computing the gravitational fields of inhomogeneous bodies

Summary A solution of the direct gravity problem for a finite body with variable density is given. The method is based on Green's formula and is applicable when a particular solution of Poisson's equation is known. The attraction due to the body is expressed by integrals over its surface The exact solution of the direct gravity problem, as known from the theory of two-dimensional fields [1–3], is closely connected with the problem of the analytic continuation of the exterior field of the attracting mass system into its interior. In the first place, this is a problem of determining the singularities of the exterior field, their distribution within the system and their nature. This approach to the solution of the direct problem is also meaningful from the point of view of determining the characteristics of the attracting system and, therefore, also of solving the inverse problem. In the case of two-dimensional fields the methods of analytical continuation were widely developed in a series of well-known papers by V. N. Strakhov, and they are mainly based on the methods of the theory of the functions of the complex variable. These methods were also successfully applied by Tsirulskii and Golizdra [1, 2] in treating the homogeneous and inhomogeneous, two-dimensional direct problem by means of Cauchy's integrals. However, as regards three-dimensional fields a number of fundamental problems has not been solved in this respect.Dedicated to 90th Birthday of Professor Frantiek Fiala     

Hidden equivalence and efficiency of gravity data interpretation

Analytical properties of the solutions of the inverse problem of gravimetry are studied in the context of the approximative approach and method of linear integral representations. A new effect of the inheritance of specific analytical properties in the solutions, which is named the effect of hidden equivalence, is revealed and investigated. This effect significantly influences the interpretation informativity in the complex models of a medium, and it should be taken into account in the geological interpretation of gravity data. Hidden equivalence is studied for both linear and nonlinear inverse problems. As an example of a nonlinear problem, the inverse problem of structural gravimetry is analyzed. The correlation between the boundary equations and boundary values of the harmonic functions is demonstrated. Methods are suggested to allow for the effects that occur by expanding the approximative approach for complex conditions (the models of the media with spatially distributed parameters) by the dynamical and criterion principles.     

一种二度体单一界面的空间域重力迭代反演方法

对二度体单一密度界面的空间域求解提出了一种改进的迭代反演方法。通过详细的改正求得深部重力异常。     

重力和重力梯度数据联合聚焦反演方法

重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合.     

Determination of the excess density in a horizontal layered model

The inversion gravity problem formulated as follows is solved. The excess density in each layer of a fixed horizontal stratified model is a function of horizontal coordinates (σ(ξ, η)) approximated by a specially constructed function. The problem is to reconstruct the function σ = σ(ξ, η) from the external gravity field. If the geological model includes more than one layer, the problem is solved with the use of a set of reference points at which the sought function is given. Variations in a gravity anomaly with respect to the field at a fixed point are used in solving this problem.     

On inversion of the second- and third-order gravitational tensors by Stokes’ integral formula for a regional gravity recovery

A regional recovery of the Earth’s gravity field from satellite observables has become particularly important in various geoscience studies in order to better localize stochastic properties of observed data, while allowing the inversion of a large amount of data, collected with a high spatial resolution only over the area of interest. One way of doing this is to use observables, which have a more localized support. As acquired in recent studies related to a regional inversion of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) data, the satellite gravity-gradient observables have a more localized support than the gravity observations. Following this principle, we compare here the performance of the second- and third-order derivatives of the gravitational potential in context of a regional gravity modeling, namely estimating the gravity anomalies. A functional relation between these two types of observables and the gravity anomalies is formulated by means of the extended Stokes’ integral formula (or more explicitly its second- and third-order derivatives) while the inverse solution is carried out by applying a least-squares technique and the ill-posed inverse problem is stabilized by applying Tikhonov’s regularization. Our results reveal that the third-order radial derivatives of the gravitational potential are the most suitable among investigated input data types for a regional gravity recovery, because these observables preserve more information on a higher-frequency part of the gravitational spectrum compared to the vertical gravitational gradients. We also demonstrate that the higher-order horizontal derivatives of the gravitational potential do not necessary improve the results. We explain this by the fact that most of the gravity signal is comprised in its radial component, while the horizontal components are considerably less sensitive to spatial variations of the gravity field.     

先验地质信息约束下的三维重磁反演建模研究 ——以安徽泥河铁矿为例

安徽泥河铁(硫)矿床是近年发现的大深度隐伏矿床,它的发现再次引发了庐枞地区深部找矿的热潮.建立矿区的三维精细地质模型,对寻找深、边部隐伏矿体,深入认识深部成矿、控矿规律意义重大.本文以泥河矿区为实例,开展先验地质信息约束的三维重磁建模研究,获得了矿区面积5.6 km²(2.8 km×2.0 km)、深度1.2 km内的三维地质模型.在三维可视化平台上对该模型进行了地质解释,全面分析了矿体、地层与次火山岩之间的空间分布及对应关系,发现铁矿主要赋存于闪长玢岩与砖桥组火山岩之间,且在玢岩穹窿地段矿体厚大,这对认识"玢岩型"铁矿的成矿模式具有重要的实际意义;三维模型的重磁正演响应基本拟合了实际重磁异常,说明在建模区域内已没有新的矿体;本文提出的三维地质建模流程可为其它地区开展类似工作提供借鉴,同时研究表明地质信息约束下的三维重磁建模研究在深、边部找矿和重磁异常的精细解剖等方面具有潜在的价值和广阔的应用前景.     

大庆外围盆地地球物理场与盆地基底特征

根据收集到的大庆外围盆地的布格重力异常和航磁异常数据,对外围盆地中的8个盆地进行了地球物理场特征的分析.对各盆地的重力异常进行了由下地壳与上地幔的密度差引起的重力效应的剥离,对磁力异常进行了由居里面起伏引起的磁力效应的剥离.在此基础上应用调和级数法和遗传算法分别反演计算了这8个盆地的重力基底和磁性顶界面,并分析了其基底特征. 各盆地重力基底在0.2~9.0 km之间变化,磁性顶界面在1.8~9.8 km之间变化,其特征反映了各盆地的基本现状.     

Global ellipsoid-referenced topographic,bathymetric and stripping corrections to gravity disturbance

This paper revisits several aspects of defining and computing the anomalous gravity data for purposes of gravimetric inversion/interpretation. Attention is paid to evaluation of a refined global topographic correction to the gravity disturbance based on the reference ellipsoid (RE) and constant reference density for solid topography onshore and sea water density for liquid topography offshore. The global bathymetric correction is discussed. Two issues associated with compilation and inversion of bathymetrically and topographically corrected gravity disturbances in regions of negative ellipsoidal (geodetic) heights are pointed out: the evaluation of normal gravity and the harmonic continuation of the gravity data. Stripping, the removal of an effect of a known density contrast, is considered also for additional geological elements such as lakes, glaciers, sedimentary basins, isostatic mountain roots, etc. The stripping corrections are discussed in the context of the gravimetric inverse problem.     

An approximation method for solving the inverse MTS problem with the use of neural networks

The work develops the approximation approach to solving the inverse MTS problem with the use of neural networks. The inverse problem is considered in model classes of parametrized geoelectric structures, whose electric conductivity is controlled by a few hundreds of macroparameters (N ∼ 300). An approximate inverse operator of the problem is constructed for each model class as a neural network, whose coefficients are determined in the process of training on a representative sample of standard examples of forward problem solutions. The problem of determination of the model class of geolectric structures corresponding to the presented input MT data is solved with the use of the neural network classifier constructed for the available set of model classes of structures. Regularizing factors and errors of the neural network method are analyzed. The operation of the algorithm is illustrated by examples of the 2-D inversion of synthetic MT data.     

Inverse problems of gravity prospecting as a decision-making problem under uncertainty and risk

A new class of algorithms for solving the inverse problems of gravity prospecting is considered. The best interpretation is selected from the set Q of the admissible versions by the optimality criteria that are borrowed from the solution-making theory and adapted for the geophysical problems. The concept of retrieving the information about the sources of gravity anomalies, which treats the result of the interpretation as a set of locally optimal solutions of the inverse problem but not as a single globally optimal solution is discussed. The locally optimal solutions of the inverse problem are sort of singularity points of set Q. They are preferable to the other admissible solutions by a certain criterion formulated in terms of the geologically important information about the anomalous bodies. The admissible versions of the interpretation of the gravimetry data that meet the criteria of the decision-making theory are the primary candidates for the singularity points. The results of the numerical calculations are presented. The set of the admissible solutions from which the locally optimal versions of interpretation are selected is formed by the modifications of the assembly method developed by V.N. Strakhov.     

Non‐linear stochastic inversion of gravity data via quantum‐behaved particle swarm optimisation: application to Eurasia–Arabia collision zone (Zagros,Iran)

Potential field data such as geoid and gravity anomalies are globally available and offer valuable information about the Earth's lithosphere especially in areas where seismic data coverage is sparse. For instance, non‐linear inversion of Bouguer anomalies could be used to estimate the crustal structures including variations of the crustal density and of the depth of the crust–mantle boundary, that is, Moho. However, due to non‐linearity of this inverse problem, classical inversion methods would fail whenever there is no reliable initial model. Swarm intelligence algorithms, such as particle swarm optimisation, are a promising alternative to classical inversion methods because the quality of their solutions does not depend on the initial model; they do not use the derivatives of the objective function, hence allowing the use of L₁ norm; and finally, they are global search methods, meaning, the problem could be non‐convex. In this paper, quantum‐behaved particle swarm, a probabilistic swarm intelligence‐like algorithm, is used to solve the non‐linear gravity inverse problem. The method is first successfully tested on a realistic synthetic crustal model with a linear vertical density gradient and lateral density and depth variations at the base of crust in the presence of white Gaussian noise. Then, it is applied to the EIGEN 6c4, a combined global gravity model, to estimate the depth to the base of the crust and the mean density contrast between the crust and the upper‐mantle lithosphere in the Eurasia–Arabia continental collision zone along a 400 km profile crossing the Zagros Mountains (Iran). The results agree well with previously published works including both seismic and potential field studies.     

Analytical Solution of Direct and Inverse Problems in the Internal Gravity Waves Studies by the Doppler Frequency Shift Method

An analytical solution of direct and inverse problems arising in the study of the internal gravity waves (IGWs) dynamic via recording of the Doppler frequency shift, is presented. The direct problem is to determine the response of the Doppler shift to IGWs in the region of the radio wave reflection point; the inverse problem is the determination of IGW parameters from data on the Doppler frequency shift. Solutions were obtained in an approximation of the isothermal ionosphere for the heights of the F-region. They are presented in a form convenient for their practical use and can have a wide range of applications, including the detection of soliton-like wave structures in the F-region of the ionosphere.     

Non-linear inversion of controlled source multi-receiver electromagnetic induction data for unexploded ordnance using a continuation method

This work adopts a continuation approach, based on path tracking in model space, to solve the non-linear least-squares problem for discrimination of unexploded ordnance (UXO) using multi-receiver electromagnetic induction (EMI) data. The forward model corresponds to a stretched-exponential decay of eddy currents induced in a magnetic spheroid. We formulate an over-determined, or under-parameterized, inverse problem. An example using synthetic multi-receiver EMI responses illustrates the efficiency of the method. The fast inversion of actual field multi-receiver EMI responses of inert, buried ordnances is also shown. Software based on the continuation method could be installed within a multi-receiver EMI sensor and used for near-real-time UXO decision-making purposes without the need for a highly-trained operator.     

General inverse of Stokes, Vening-Meinesz and Molodensky formulae

The undulation of the geoid, the gravity anomaly and the deflection of the vertical are the three basic observations describing the shape and the gravity field of the earth. The Stokes’ formula that computes the undulation of the geoid using the gravity anomaly on the geoid under spherical approximate conditions was first put forward by Stokes[1]. According to Stokes’ theory, The Vening-Meinesz formula that computes the meridian and the prime vertical components of the deflection of the ve…     

Method of finding point (dipole) solutions of the potential field inverse problem

Summary Instead of solving the potential field inverse problem by means of optimal point masses a set of point masses is considered to represent one of the most concentrated (pressed) bodies of the family. This body is divided into layers, each of them represented by some quantity of point masses.     

Geoid determination by spectral combination of an Earth gravitational model with airborne and terrestrial gravimetry — a theoretical study

Today air-gravimetry is a versatile technique to quickly collect gravity data over large regions, where terrestrial gravity data are sparse and/or of poor quality. The method requires the data to be downward continued to sea level for use in geoid determination, an inverse problem operation that calls for smoothing of the data and/or the kernel function involved (in either spectral or space domain). In this purely theoretical study we avoid this separate computational step by performing direct geoid estimation by so-called spectral combination/filtering of the data, which includes terrestrial gravimetry, airgravimetry, an Earth Gravitational Model (EGM) as well as their signal and error degree variances. Each derived geoid estimator is presented as the sum of one or two integral formulas and the harmonic series of the EGM together with the expected mean square error of the estimator. The article is limited to a theoretical study, leaving its practical tests for future investigation.