共查询到20条相似文献，搜索用时 183 毫秒
1.
提出了利用地面重力异常数据计算地面扰动位径向二阶梯度时，远距离带区影响的误差公式，结果以级数形式给出，以其数学模型为基础，利用全球重力位模型GPM98，计算确定了不同平均异常网格大小情况下截断半径取一系列值时的远区影响，在将卫星重力梯度和航空重力梯度观测数据与地面重力数据进行联合处理时，文中提出了远区影响的数量关系，对于确定合理的截断半径大小具有一定的参考价值。 相似文献
2.
由地面重力数据确定扰动位径向二阶梯度 总被引：5，自引：0，他引：5
提出利用地面重力异常数据计算地面扰动位径向二阶梯度，将该梯度的积分表达式转换为卷形式的谱表达式，便于应用ＦＦＴ／ＦＨＴ技术进行快速计算。 相似文献
3.
提出了利用地面重力异常数据计算地面扰动位径向二阶梯度时,远距离带区影响的误差公式,结果以级数形式给出;以其数学模型为基础,利用全球重力位模型GPM98,计算确定了不同平均异常网格大小情况下截断半径取一系列值时的远区影响;在将卫星重力梯度和航空重力梯度观测数据与地面重力数据进行联合处理时,文中提出的远区影响的数量关系,对于确定合理的截断半径大小具有一定的参考价值. 相似文献
4.
为实现大范围、高精度基准重力梯度数据库的构建,考虑到重力梯度场对地形质量的敏感效应,一般利用恒密度数字高程模型来求取重力梯度值,从而忽略了地形密度变化以及水准面以下密度异常对重力梯度的影响。根据重力位理论中求解边值问题的数值应用方法,直接利用重力异常数据求取重力梯度场,弥补了密度变化和密度异常在重力梯度上的反映。根据模型算例和实测重力异常数据求取了剖面重力梯度值,结果表明,限于重力数据空间分辨率的影响,利用重力异常数据可恢复中长波段重力梯度场。该方法与地形数据求取重力梯度和卫星重力梯度测量等方法技术相结合,对重力梯度数据库的建设具有实际应用价值。 相似文献
5.
利用重力异常信息获取梯度全张量 总被引：1，自引：0，他引：1
重力位二阶梯度张量反应了重力异常的空间变化率[1],特别适合于探测或研究局部的小地质体及其细节.本文在地球外部重力场基本理论的基础上,围绕如何利用重力异常信息计算获取重力梯度张量的理论与方法展开了研究.通过实验对比,验证了重力梯度较之于重力异常更能反应细致的局部特征. 相似文献
6.
相对于传统的重力测量手段，重力梯度测量能够以更高的灵敏度和分辨能力反映出地下密度异常体的结构特征。由于拉格朗日经验参数在实测数据反演中存在不确定性，对预条件共轭梯度反演算法加以改进，利用L曲线的拐点值代替原反演算法中的拉格朗日经验参数作为正则化参数；为改善反演中存在的病态性问题并减弱核函数的快速衰减，将地下模型改进为不等间隔模型；为改善反演中解的非唯一性，利用重力梯度的5个独立分量进行联合反演；通过对澳大利亚Kauring试验场航空重力梯度张量进行联合反演，得到该地区异常体的三维密度分布，将重力梯度联合反演结果与之前的重力反演结果对比分析，发现在中心异常体附近沿线还分布着多个异常块体。结果表明，改进后的算法能够有效地利用实测重力梯度数据反演出密度异常体的分布信息。 相似文献
7.
重力梯度的空间传播特性 总被引：1，自引：0，他引：1
基于球谐函数谱分析理论，导出了重力梯度张量住全球平均意义下的功率谱表达式，从理论上研究和揭示了重力梯度的空间传播特性。利用重力异常阶方差模型和EGM96重力位模型，研究了苇力梯度张量和扰动重力梯度张量随高度的衰减特性，据此分析了卫星重力梯度测量所能恢复重力场的最高阶数，并讨论了对卫星重力梯度仪和卫星定轨的精度要求。 相似文献
8.
给出了由地面重力数据计算外空扰动重力矢量的公式，并根据Wong and Gore截断理论给出了外区重力异常对计算点影响的公式。根据外空扰动重力与地面数据分辨率及其覆盖范围之间的关系，将重力异常分成不同频段、不同分辨率，分别计算了截断误差。利用全球重力位模型，计算出不同频段的截断误差，并给出了各频段相应的积分半径。对于模拟和检验航空矢量重力数据有一定的参考价值。 相似文献
9.
断裂构造研究是重力解释的一项重要工作,与构造单元划分密切相关。全张量重力梯度数据以其信息量大、含有更高频的信号成分,能更好地描述小的异常特征等优点在地球物理领域中得到广泛应用。基于全张量重力梯度组合研究中国南海断裂识别及提取方法。首先,比较多种重力梯度边界识别方法,包括直接利用重力梯度三分量法和全张量梯度组合法,分析它们的优缺点。通过对比分析,传统重力梯度三分量方法不能有效地均衡深浅异常的振幅,当异常中同时出现正负异常可能产生假的边界结果。全张量重力梯度组合法不仅可以有效地避免传统方法的缺陷,而且获得的边界还具有良好的连续性和收敛性。其次,利用改进的边缘检测计算理论边界提取法确定断裂的精确平面位置,得到了与全张量梯度组合法一致的结果。由此推断,南海断裂以北东走向和北西走向为主,北东东、北西、东西和近南北走向为辅。 相似文献
10.
重力测量数据存在地形数据产生的高频分量的影响,高精度地形数据正演重力梯度也能较好地反映重力局部高频特征。为获得高精度重力梯度数据,实现基准梯度数据库精确快速构建,研究了利用数字高程模型正演重力梯度的频率域快速计算方法,推导出基于余弦变换的Parker正演重力梯度理论公式。数值实验结果表明,余弦变换频率域正演方法平均绝对误差可达到0.5E左右精度要求,与傅里叶变换正演方法相比误差可减小3dB左右,与棱柱法等空间域正演方法相比,该方法计算规模小,速度优势明显。 相似文献
11.
根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。 相似文献
12.
利用最小二乘直接法反演卫星重力场模型的MPI并行算法 总被引：2，自引：0，他引：2
针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和GOCE的解算模型能够实现二者独立信息的频谱互补,表明本文方法可高效稳定地恢复高阶次地球重力场模型。 相似文献
13.
GOCE采用的高低卫-卫跟踪和卫星重力梯度测量技术在恢复重力场方面各有所长并互为补充,如何有效利用这两类观测数据最优确定地球重力场是GOCE重力场反演的关键问题。本文研究了联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的最小二乘谱组合法,基于球谐分析方法推导并建立了卫星轨道面扰动位T和径向重力梯度Tzz、以及扰动位T和重力梯度分量组合{Tzz-Txx-Tyy}的谱组合计算模型与误差估计公式。数值模拟结果表明,谱组合计算模型可以有效顾及各类数据的精度和频谱特性进行最优联合求解。采用61天GOCE实测数据反演的两个180阶次地球重力场模型WHU_GOCE_SC01S（扰动位和径向重力梯度数据求解）和WHU_GOCE_SC02S（扰动位和重力梯度分量组合数据求解）,结果显示后者精度优于前者,并且它们的整体精度优于GOCE时域解,而与GOCE空域解的精度接近,验证了谱组合法的可行性与有效性。 相似文献
14.
多种类型高分辨率重力场数据的不断增加,使得在局部范围内精化重力场模型成为了可能。本文采用Abel-Poisson核将重力场量表示成有限个径向基函数线性求和的形式,对局部区域的多种重力场数据进行联合建模。为了提高运算速度,运用了基于自适应精化格网算法的最小均方根误差准则(RMS)来求解径向基函数平均带宽。以南海核心地区为例,联合两种不同类型、不同分辨率的重力场资料(大地水准面起伏6'×6'、重力异常2'×2'),构建了局部区域高分辨率的重力场模型。所建模型表示的重力场参量达到了2'×2'的分辨率,对原始的重力异常数据(2'×2')拟合的符合程度达到±0.8×10-5m/s2。结果表明,利用径向基函数方法进行局部重力场建模,避免了球谐函数建模收敛慢的问题,有效提高了模型表示重力场的分辨率。 相似文献
15.
导航用重力梯度基准图构建方法的比较研究 总被引：1，自引：0，他引：1
基准图是重力梯度辅助导航实用化研究的关键技术,由于依靠实测构造大范围基准图还需要较长时间,理论模型逼近重力梯度开始成为可行的构图方法。系统研究了一种基于自由空间异常或布格异常数据的构图方法,避免了相关研究构建梯度图时较少考虑标准测量数据,构图结果物理意义不明确和缺乏实用性的缺点。以重力归算理论推导了构图方法原理,然后用频域算法实现构图过程中两个关键计算过程,通过算例分析了构图效果。所构基准图数值、轮廓等关键特征均与实测重力梯度吻合较好,证明了该方法在重力梯度辅助导航基准图构建中的实用价值。 相似文献
16.
Impact of accounting for coloured noise in radar altimetry data on a regional quasi-geoid model 总被引：1，自引：1，他引：0
We study the impact of an accurate computation and incorporation of coloured noise in radar altimeter data when computing a regional quasi-geoid model using least-squares techniques. Our test area comprises the Southern North Sea including the Netherlands, Belgium, and parts of France, Germany, and the UK. We perform the study by modelling the disturbing potential with spherical radial base functions. To that end, we use the traditional remove-compute-restore procedure with a recent GRACE/GOCE static gravity field model. Apart from radar altimeter data, we use terrestrial, airborne, and shipboard gravity data. Radar altimeter sea surface heights are corrected for the instantaneous dynamic topography and used in the form of along-track quasi-geoid height differences. Noise in these data are estimated using repeat-track and post-fit residual analysis techniques and then modelled as an auto regressive moving average process. Quasi-geoid models are computed with and without taking the modelled coloured noise into account. The difference between them is used as a measure of the impact of coloured noise in radar altimeter along-track quasi-geoid height differences on the estimated quasi-geoid model. The impact strongly depends on the availability of shipboard gravity data. If no such data are available, the impact may attain values exceeding 10 centimetres in particular areas. In case shipboard gravity data are used, the impact is reduced, though it still attains values of several centimetres. We use geometric quasi-geoid heights from GPS/levelling data at height markers as control data to analyse the quality of the quasi-geoid models. The quasi-geoid model computed using a model of the coloured noise in radar altimeter along-track quasi-geoid height differences shows in some areas a significant improvement over a model that assumes white noise in these data. However, the interpretation in other areas remains a challenge due to the limited quality of the control data. 相似文献
17.
Gravity gradient modeling using gravity and DEM 总被引：2，自引：0，他引：2
A model of the gravity gradient tensor at aircraft altitude is developed from the combination of ground gravity anomaly data
and a digital elevation model. The gravity data are processed according to various operational solutions to the boundary-value
problem (numerical integration of Stokes’ integral, radial-basis splines, and least-squares collocation). The terrain elevation
data are used to reduce free-air anomalies to the geoid and to compute a corresponding indirect effect on the gradients at
altitude. We compare the various modeled gradients to airborne gradiometric data and find differences of the order of 10–20 E
(SD) for all gradient tensor elements. Our analysis of these differences leads to a conclusion that their source may be primarily
measurement error in these particular gradient data. We have thus demonstrated the procedures and the utility of combining
ground gravity and elevation data to validate airborne gradiometer systems. 相似文献
18.
Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition 总被引：3，自引：0，他引：3
We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach
is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm
for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and
block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure.
Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the
algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation
of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime
results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming
satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month
of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations
and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion
strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development
stage, the LSQR method appears ready to deal with real-data applications. 相似文献
19.
20.
在重力聚焦反演基础上提出多尺度源网聚焦反演算法.首先,对源网进行粗网格剖分,用共轭梯度法求解粗网格源网模型的聚焦解,直到拟合差下降至设定的数值;然后,将粗网格得到的密度映射到细网格;最后,以细网格模型为初始模型,进一步迭代直到拟合差下降至符合反演要求.模型试验结果显示,相比于固定源网反演,多尺度源网聚焦反演迭代的总耗时... 相似文献