Spherical Harmonic Analysis of Gravitational Curvatures and Its Implications for Future Satellite Missions

In this study we assume that a gravitational curvature tensor, i.e. a tensor of third-order directional derivatives of the Earth’s gravitational potential, is observable at satellite altitudes. Such a tensor is composed of ten different components, i.e. gravitational curvatures, which may be combined into vertical–vertical–vertical, vertical–vertical–horizontal, vertical–horizontal–horizontal and horizontal–horizontal-horizontal gravitational curvatures. Firstly, we study spectral properties of the gravitational curvatures. Secondly, we derive new quadrature formulas for the spherical harmonic analysis of the four gravitational curvatures and provide their corresponding analytical error models. Thirdly, requirements for an instrument that would eventually observe gravitational curvatures by differential accelerometry are investigated. The results reveal that measuring third-order directional derivatives of the gravitational potential imposes very high requirements on the accuracy of deployed accelerometers which are beyond the limits of currently available sensors. For example, for orbital parameters and performance similar to those of the GOCE mission, observing third-order directional derivatives requires accelerometers with the noise level of \({\sim}10^{-17}\,\hbox {m}\,\hbox {s}^{-2}\) Hz\(^{-1/2}\).     

引力梯度归算的模拟计算

引入引力梯度不变量后可将GOCE卫星观测数据简化成扰动位的径向二阶导数边界条件,由于卫星轨道不规则性给解算带来了很大的困难,因此解算过程中需要将在轨道上建立的边界条件延拓到如平均球面这样的规则曲面上来,由此便需要对引力梯度的延拓或归算等问题展开研究.本文依据Taylor展开讨论了引力梯度的归算方法,并针对EGM2008...     

卫星重力梯度测量与地球引力场的精度研究

本文根据地球引力位的球谐函数展开式,利用重力梯度张量各分量导出了位系数模型的精度估计公式.从三方面进行了研究:假定卫星重力梯度仪测量精度,探讨用重力梯度数据确定地球重力场模型的精度;求出位系数模型和大气阻力引起的重力梯度卫星的轨道误差;最后,反求轨道误差和位系数误差对重力梯度测量值的影响.数值计算表明,与地面技术和常规卫星方法相比,卫星梯度测量可使重力场模型的精度至少提高3-5倍;利用重力梯度张量全分量求得的重力值精度比单用径向分量V_rr的结果提高40％以上;若仅顾及位系数模型和大气阻力误差,则轨道误差对梯度测量值的影响△V_i3（i=3,2,1）至少可分别在1/4和1/3弧圈内达到△V_i3≤σ（仪器精度）.     

磁变仪布局的研究

本文提出了改进磁变仪布局的必要性和布局标准,对磁变仪定向及其稳定性、标度值稳定性和分量仪间允许的最小间距进行了研究,提出了可以提高磁变仪记录质量的改进布局和满足布局全部标准的新布局,改进布局特别是新布局可以显著提高磁变仪保持定向、抗干扰能力,获得记录标度值稳定的较理想的磁照图。     

重力全张量数据联合欧拉反褶积法研究及应用

全张量测量技术是在空中或海上用加载了多个加速度计的移动平台技术测量位场的五个独立分量.各张量分量包含不同方向的地下地质体信息,水平张量分量T_(xx)、T_(yy)、T_(xy)、T_(xz)、T_(yz)通常用于识别和映射与地质构造或地层变化有关的测量区域中的目标,垂直张量分量Tzz用于估计深度.然而,这些分量传统上是彼此分开解释,经常遇到错失关键信息的风险.本文所用全张量欧拉反褶积是在单独z方向的欧拉反演基础上发展而来的,它融合了重力异常垂直分量以及其三个方向导数、水平分量以及其三个方向导数.全张量数据信息得以有效应用的同时,欧拉反褶积结果也比常规欧拉反褶积结果更加收敛.最后,结合美国墨西哥湾地区实测航空FTG数据,用重力梯度张量数据进行联合欧拉三维反演研究,有效的识别岩盖的边界信息,划分岩盖范围,为进一步研究盖层底下深部复杂地质情况提供可靠的解释结果.     

Non-singular expression for the gradient of the Earth gravitational potential in the local north-oriented ellipsoidal reference frame

In this contribution we continue our earlier research, concerning the ellipsoidal harmonic expansions of the Earth disturbing gravitational potential and its derivatives on an external reference ellipsoid confocal with respect to the normal ellipsoid and close to it. One of the results of the previous investigation is represented by a new expression for the derivative of the Jekeli’s Legendre function of the second kind, entering the ellipsoidal harmonics in the potential derivative. The derived expression depends on two Gauss hypergeometric functions which converge better than the hypergeometric functions of other authors. In the present paper we construct another expression for the derivative of the Jekeli’s Legendre function, depending on two alternative hypergeometric functions. While our earlier hypergeometric series in the expression for the derivative of this function converge better when the orders of the terms do not exceed a half of their degrees, the series constructed in the present paper converge more rapidly when the orders surpass a half of the degrees. We deduce an improved expression for the derivative of the Jekeli’s Legendre function by combining these results and then construct a corresponding new expression for the derivative of the disturbing potential. This expression is applied for constructing non-singular expressions for the components of the gradient of the potential in the local north-oriented ellipsoidal reference frame. The new expressions for these components have no these deficiencies and the expression for the potential gradient depends on very quickly convergent hypergeometric series.     

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.     

Computing three-dimensional gravitational fields with equivalent sources

Existing techniques for computing the gravitational field due to a homogeneous polyhedron all transform the required volume integral, expressing the field due to a volume distribution of mass, into a surface integral, expressing the potential due to a surface mass distribution over the boundary of the source body. An alternative representation is also possible and results in a surface integral expressing the potential due to a variable-strength double layer located on the polyhedral source boundary. Manipulation of this integral ultimately allows the gravitational field component in an arbitrary direction to be expressed as a weighted sum of the potentials due to two basic source distributions. These are a uniform-strength double layer located on all faces and a uniform-strength line source located along all edges. The derivatives of the gravitational field components can also be expressed in a similar form as can the magnetic field components due to a homogeneous magnetic polyhedron. It follows that the present approach can be used to generate a universal program capable of modelling all the commonly used potential field responses due to 3D bodies of arbitrary shape.     

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.     

Spectral expressions for modelling the gravitational field of the Earth’s crust density structure

We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling of major known density structures within the Earth’s crust (excluding the ocean density contrast) beneath the geoid surface. The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples, are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined relative to the adopted value of the reference crustal density of 2670 kgm⁻³. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface. The maxima of the gravitational signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits. The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental crustal thickness with variable geological structure.     

MD－2型磁变仪数据采集系统的研制

介绍了为俄罗斯石英光电磁变仪配接的ＭＤ－２型磁变仪数据采集系统的设计、性能特点，及其与俄磁变仪组成的数字磁变仪系统在北京地磁中心台完成数字化记录的应用结果．     

MD－2型磁变仪数据采集系统的研制

介绍了为俄罗斯石英光电磁变仪配接的ＭＤ－２型磁变仪数据采集系统的设计、性能特点，及其与俄磁变仪组成的数字磁变仪系统在北京地磁中心台完成数字化记录的应用结果．     

Experimental and numerical study of the validity of Hele–Shaw cell as analogue model for variable-density flow in homogeneous porous media

Several laboratory experiments were conducted to identify the validity domain under which a Hele–Shaw cell may serve as a suitable analogue for variable-density flow in homogeneous porous media. These experiments are concerned with the injection into a Hele–Shaw cell of a salt solute at different concentrations and flow rates. The experimental data analysis highlighted two types of mixing zone shape: with and without ‘fingers’. A semi-empirical criterion based on the ratio between gravitational and injected velocities was used to forecast the change from one shape to another. The experimental data were then analysed using numerical solutions of the classical Hele–Shaw equations by taking into account an anisotropic dispersion tensor whose components depend on fluid density gradients. The good agreement between experimental and numerical results clearly shows that the validity of the concentration-dependent dispersion tensor strongly depends on the local Péclet number variation. For Péclet numbers lower 50, the Hele–Shaw cell can be considered as an analogous model of a homogeneous and isotropic 2D porous medium. It can be successfully used to study, at the laboratory scale, the gravitational instability effects induced by flow and transport phenomena into a porous medium.     

数字磁变仪的设计与实现

本文介绍了石英光电磁变仪及其数字化采集，记录的设计和实现，并讨论了数字磁变仪系统在地磁台站的应用结果。     

多功能照相记录磁变仪电源

介绍了多功能照相记录磁变仪电源的主要特性、各部分电路的工作原理、制作方法和使用注意事项，为仿制提供了可靠的依据。该电源还适用于倾斜等照相记录的仪器。     

Noise filtering of full-gravity gradient tensor data

In oil and mineral exploration, gravity gradient tensor data include higher-frequency signals than gravity data, which can be used to delineate small-scale anomalies. However, full-tensor gradiometry (FTG) data are contaminated by high-frequency random noise. The separation of noise from high-frequency signals is one of the most challenging tasks in processing of gravity gradient tensor data. We first derive the Cartesian equations of gravity gradient tensors under the constraint of the Laplace equation and the expression for the gravitational potential, and then we use the Cartesian equations to fit the measured gradient tensor data by using optimal linear inversion and remove the noise from the measured data. Based on model tests, we confirm that not only this method removes the high-frequency random noise but also enhances the weak anomaly signals masked by the noise. Compared with traditional low-pass filtering methods, this method avoids removing noise by sacrificing resolution. Finally, we apply our method to real gravity gradient tensor data acquired by Bell Geospace for the Vinton Dome at the Texas-Louisiana border.     

Time-space multiscale analysis by use of tensor product wavelets and its application to hydrology and GRACE data

This paper presents a wavelet analysis of temporal and spatial variations of the Earth’s gravitational potential based on tensor product wavelets. The time-space wavelet concept is realized by combining Legendre wavelets for the time domain and spherical wavelets for the space domain. In consequence, a multiresolution analysis for both temporal and spatial resolution is formulated within a unified concept. The method is then numerically realized by using first synthetically generated data and finally two real data sets.     

CHFR测井原理及影响因素研究

应用电磁场理论及传输线模型可建立CHFR测井响应方程,应用传输线模型可研究CHFR测井的测量原理及其影响因素.研究表明中场区电场电位二阶导数稳定,电位二阶导数对地层电阻率变化敏感,在套管外仅存在径向分量,该区延伸范围大,通过测量测点中场区的电位二阶导数和电位可实现过套管地层电阻率的测量.测量频率、电极距、水泥环、套管非均质性等因素都对测井响应产生影响.     

On the possibility of detecting complex sources of earthquakes from data of the centroid moment tensor catalog

In their analysis of the crustal stress-strain state and the study of strong earthquake sources, a number of authors use the parameter ? measuring the deviation of the focal mechanism from the double-couple model. Can one obtain reliable estimates of ? from regular determinations of the seismic moment tensor components made at world seismological centers and use these estimates to gain constraints on the rupture at the source? Comparing centroid moment tensor solutions of Harvard University with determinations made at the National Earthquake Information Center (United States) and analyzing some strong earthquakes, we arrive at a negative answer to this question. The values of ? and the consideration of determination uncertainties in seismic moment tensor components cannot provide a reliable idea of the source complexity, so that independent geological and geophysical information should be invoked.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.