基于改进Runge定理的调和延拓及其相应的边值问题

利用改进的Runge定理引入了外大地水准面和外重力位等新概念 ,建立了关于外扰动位的边值问题 ,该问题同时具备了Stokes问题和Mododensky问题的优点 .作为一个整体 ,讨论了外正高的计算和地面重力的归算 ,给出了在O(T2 )量级的精度下确定外大地水准面、地面以及地球外部重力场的方法 .     

海域流动点外部扰动引力无奇异计算模型

针对海域重力场变化特征和远程飞行器机动发射保障应用需求,本文分析研究了地球外部空间扰动引力三类传统计算模型的技术特点及其适用性,指出了采用表层法作为海域流动点扰动引力计算模型的合理性及需要解决的关键问题,分析论证了空中扰动引力计算对地面观测数据的分辨率和精度要求,提出通过引入局部积分域恒等式变换、局域泰勒级数展开和非网格点内插方法,消除表层法计算模型积分奇异性固有缺陷的研究思路,进而推出了适合于海域流动点应用的扰动引力无奇异计算模型,较好地满足了全海域和全高度段对局部扰动重力场快速赋值的实际需求.以超高阶全球位模型EGM2008作为标准场,通过数值计算验证了无奇异计算模型的可行性和有效性,在重力场变化比较剧烈的海沟区,该模型的计算精度优于2×10^-5m·s^-2.     

地面扰动重力垂直梯度的确定

鉴于扰动重力垂直梯度在大地测量和物探中有重要作用,在未直接测量的情况下,对如何由地面重力异常及地形数据求取扰动垂直梯度进行了研究和分析,认为在被研点附近要求重力点分布密集,且精度不低于1×10-6ms-2;对中央区域积分奇异性问题也作了讨论;此外,扰动重力垂直梯度如何在大地测量和地球物理中进一步得到应用,以及为何用扰动重力垂直梯度代替重力异常垂直梯度也给予了简明的论述.     

顾及远区影响的向下延拓实用算法

向下延拓是Stokes边值解算中必不可少的环节.为了研究向下延拓的远区影响,本文在泊松实用公式的基础上给出了顾及远区影响的向下延拓算法的推导过程.对4°×4°实验区内5'×5'分辨率模型扰动重力的向下延拓结果表明:对于地面重力数据的向下延拓,积分半径为1°、2°时,顾及远区影响的算法可将精度分别提高0.36 m Gal、0.17 m Gal;对于4000 m航高的重力数据的向下延拓,积分半径为1°、2°时,该算法可将精度分别提高1.40 m Gal、0.67 m Gal.延拓高度越高,积分半径越小,顾及远区影响的算法对精度的改善效果越明显.     

利用海面全极化微波辐射获取毛细重力波谱的研究

通常利用双尺度模型计算海面微波辐射的原则是, 先计算大波斜率密度函数对电磁扰动权函数的影响, 再引入海面毛细重力波谱的效应. 电磁扰动权函数与海面毛细重力波谱无关, 而且海浪长波谱和毛细重力波谱可以根据分界波数分开计算; 依据这一点, 推导并获得了一个简单的计算海面微波辐射的新模型. 依据该模型, 海面微波辐射可以表示为电磁扰动权函数对海面曲率谱的积分, 其优点是电磁扰动权函数和海面曲率谱对亮温的作用可以进行分解研究. 基于该模型, 进一步讨论了利用多频全极化辐射计反演海面毛细重力波谱的可行性.     

联合不同类型重力测量数据确定地球重力场模型的迭代法

不同的重力测量数据包含了不同波段的地球重力场信息,因此要恢复更高精度的地球重力场模型,就必须对不同类型的重力测量数据进行联合处理.以地面重力异常Δg为例,推导了利用迭代法联合不同类型重力测量数据反演地球重力场模型的基本原理公式,并给出了其具体实现步骤,接着采用全球的重力异常Δg数据和扰动位T数据,基于迭代法对卫星重力梯度SGG数据解算的重力场模型进行了进一步的精化.结果表明,初始的卫星重力梯度SGG模型和经过全球重力异常Δg数据精化后的模型,其对应的累计大地水准面误差分别达到1.128cm和0.048cm、累计重力异常误差分别达到0.416mGal和0.018mGal的精度;在经过全球扰动位T数据进一步精化后的模型,其对应的累计大地水准面误差达到0.043cm、累计重力异常误差达到0.016mGal的精度.     

GPS-重力边值问题

详细研究了在O(T²)精度下如何利用地面GPS和重力观测值确定地球外部重力场的问题, 以参考椭球面为界面建立起了完整的理论和相应的精度判断标准. 具体求解方法主要分成两个步骤: 地面重力观测值到参考椭球面上的调和延拓; 参考椭球面上Neumann外问题的积分解式. 得到的积分解式在O(T²)精度下包含了重力场频谱的全部信息; 得到的公式不仅实用, 而且理论上有所突破, 在O(T²)精度下解决了椭球界面边值问题的积分解.     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解，在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式，从而为研究物理大地测量中边值问题的求解提供了新的依据     

基于轨道扰动引力谱分析的方法确定低低观测重力卫星反演重力场的空间分辨率

基于低低卫-卫跟踪重力卫星的轨道特性,从垂直和水平两个方向计算了重力卫星高空扰动引力,并根据其谱特性及星载加速度的测量噪声水平分析了重力卫星能反演重力场的阶数.利用EGM96重力场模型分别计算了400 km、450 km和500 km 轨道高度处重力卫星受到的扰动引力谱及扰动引力谱的平均量级,分析其垂直特性表明:在三个轨道高度处能分别能反演150、140和130阶的重力场模型.利用两颗同轨重力卫星相距220 km的特性,计算了400 km、450 km和500 km 轨道高度处纬度相差2°的两颗卫星纬向扰动引力差,即扰动引力水平分量,分析其谱特性,表明:重力卫星能反演至117阶的地球重力场模型.     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解，在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式，从而为研究物理大地测量中边值问题的求解提供了新的依据     

海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.     

Conditionality of inverse solutions to discretised integral equations in geoid modelling from local gravity data

The eigenvalue decomposition technique is used for analysis of conditionality of two alternative solutions for a determination of the geoid from local gravity data. The first solution is based on the standard two-step approach utilising the inverse of the Abel-Poisson integral equation (downward continuation) and consequently the Stokes/Hotine integration (gravity inversion). The second solution is based on a single integral that combines the downward continuation and the gravity inversion in one integral equation. Extreme eigenvalues and corresponding condition numbers of matrix operators are investigated to compare the stability of inverse problems of the above-mentioned computational models. To preserve a dominantly diagonal structure of the matrices for inverse solutions, the horizontal positions of the parameterised solution on the geoid and of data points are identical. The numerical experiments using real data reveal that the direct gravity inversion is numerically more stable than the downward continuation procedure in the two-step approach.     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

Mass heterogeneities and convection in the earth's mantle inferred from gravity and core-mantle boundary irregularities

Mass heterogeneities in the earth's mantle are retrieved from the gravity data and the topography of the core-mantle boundary as well as the topography of the earth's surface. A mantle circulation induced by the heterogeneities is modelled by solving the Stokes problem for incompressible Newtonian fluid. The derived models of mantle motions correlate well with the plate tectonics and point at a close relation between the surface tectonic activity and the processes in the vicinity of the core-mantle boundary.     

Non-singular time domain BEM with applications to 3D inertial soil–structure interaction

A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies.     

A RECONCILIATION OF MATHEMATICAL MODELS FOR SPONTANEOUS MINERALIZATION POTENTIALS1

Two methods for computing spontaneous mineralization potentials in the region external to the source body are reviewed. The first of these is a long-established technique in which the causation is assumed to be a distribution of simple current source on the boundary of the mineralization. The second is a more recent technique which assumes a surface distribution of current dipole moment (double layer) along the boundary of the source body. The former technique is a special case of a more general spontaneous potential (SP) model in which the source is a density of current dipole moment (current polarization) distributed throughout the mineralization. As far as the potentials in the region external to the source body are concerned this current polarization can be simply related to an equivalent double layer source function, i.e. the potential discontinuity produced over the boundary of the mineralization by an equivalent double layer model. This simple relationship suggests an integral equation technique for the exact numerical solution of boundary value problems appropriate to the polarization model for spontaneous mineralization potentials. The technique is applied to exploring the justification of interpreting mineralization self-potentials by the traditional approach.     

Least-squares modification of extended Stokes’ formula and its second-order radial derivative for validation of satellite gravity gradiometry data

The gravity anomalies at sea level can be used to validate the satellite gravity gradiometry data. Validation of such a data is important prior to downward continuation because of amplification of the data errors through this process. In this paper the second-order radial derivative of the extended Stokes’ formula is employed and the emphasis is on least-squares modification of this formula to generate the second-order radial gradient at satellite level. Two methods in this respect are proposed: (a) modifying the second-order radial derivative of extended Stokes’ formula directly, and (b) modifying extended Stokes’ formula prior to taking the second-order radial derivative. Numerical studies show that the former method works well but the latter is very sensitive to the proper choice of the cap size of integration and degree of modification.     

A theoretical study on terrestrial gravimetric data refinement by Earth gravity models

The idea of this paper is to present estimators for combining terrestrial gravity data with Earth gravity models and produce a high‐quality source of the Earth's gravity field data through all wavelengths. To do so, integral and point‐wise estimators are mathematically developed, based on the spectral combination theory, in such a way that they combine terrestrial data with one and/or two Earth gravity models. The integral estimators are developed so that they become biased or unbiased to a priori information. For testing the quality of the estimators, their global mean square errors are generated using an Earth gravity model08 model and one of the recent products of the gravity field and steady‐state ocean circulation explorer mission. Numerical results show that the integral estimators have smaller global root mean square errors than the point‐wise ones but they are not efficient practically. The integral estimator of the biased type is the most suited due to its smallest global root mean square error comparing to the rest of the estimators. Due largely to the omission errors of Earth gravity models the point‐wise estimators are not sensitive to the Earth gravity model commission error; therefore, the use of high‐degree Earth gravity models is very influential for reduction of their root mean square errors. Also it is shown that the use of the ocean circulation explorer Earth gravity model does not significantly reduce the root mean square errors of the presented estimators in the presence of Earth gravity model08. All estimators are applied in the region of Fennoscandia and a cap size of 2° for numerical integration and a maximum degree of 2500 for generation of band‐limited kernels are found suitable for the integral estimators.     

融合多源数据构建区域航空重力梯度扰动全张量

通过联合全球重力位模型(EGM2008)、航空重力扰动数据和剩余地形模型(RTM)数据,基于频谱域(二维FFT变换)和空间域(Stokes数值积分)算法对毛乌素测区GT-2A航空重力测量系统采集的空中测线后处理重力扰动数据进行解算,构建了该地区的航空重力梯度扰动全张量.(1)残余航空重力扰动延拓结果表明:残余航空重力扰动经向下延拓至大地水准面,再向上延拓至航空高度后与原数据差值的标准差为1.0078 mGal,考虑边缘效应后,内缩计算范围得到的差值标准差减小至0.1269 mGal.(2)基于残余重力扰动数据(原航空高度数据及向下延拓数据),通过不同方案解算得到的梯度扰动结果表明:两种方案得到的研究区域重力梯度扰动各分量之差的最大标准差为6.4798E(Γ_(yz)分量),最小标准差为2.6968E(Γ_(xy)分量),内缩计算范围后得到的差值标准差最大值为1.8307E(Γ_(zz)分量),最小值为0.7223E(Γ_(yz)分量).本文的思路和方法可为未来我国自主构建航空重力梯度标定场提供参考.