首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We compute globally the consolidated crust-stripped gravity disturbances/anomalies. These refined gravity field quantities are obtained from the EGM2008 gravity data after applying the topographic and crust density contrasts stripping corrections computed using the global topography/bathymetry model DTM2006.0, the global continental ice-thickness data ICE-5G, and the global crustal model CRUST2.0. All crust components density contrasts are defined relative to the reference crustal density of 2,670 kg/m3. We demonstrate that the consolidated crust-stripped gravity data have the strongest correlation with the crustal thickness. Therefore, they are the most suitable gravity data type for the recovery of the Moho density interface by means of the gravimetric modelling or inversion. The consolidated crust-stripped gravity data and the CRUST2.0 crust-thickness data are used to estimate the global average value of the crust-mantle density contrast. This is done by minimising the correlation between these refined gravity and crust-thickness data by adding the crust-mantle density contrast to the original reference crustal density of 2,670?kg/m3. The estimated values of 485 kg/m3 (for the refined gravity disturbances) and 481?kg/m3 (for the refined gravity anomalies) very closely agree with the value of the crust-mantle density contrast of 480?kg/m3, which is adopted in the definition of the Preliminary Reference Earth Model (PREM). This agreement is more likely due to the fact that our results of the gravimetric forward modelling are significantly constrained by the CRUST2.0 model density structure and crust-thickness data derived purely based on methods of seismic refraction.  相似文献   

2.
物理大地测量学与地球物理学结合中的有关问题评注   总被引:2,自引:0,他引:2  
本文依据大地重力学反问题的有关研究,对重力场边界扰动量的本征含义进行分析,由此揭示出它们在逼近地球实际重力场源过程中所具有的物理意义及存在的不足。然后进一步论及正常场源的非物理性上,并对其解决途径进行了探讨,最后对大地得力学反问题研究中的物理核--正交分解定理描述的场源量所具有的实际意义及其在地球物理反演中的具有的作用进行了分析和评注。  相似文献   

3.
Compilation of the bathymetrically and topographically corrected gravity disturbance, the so called BT disturbance, for the purpose of gravity interpretation/inversion, is investigated from the numerical point of view, with special emphasis on regions of negative heights. In regions of negative ellipsoidal (geodetic) heights, such as the Dead Sea region onshore or offshore areas of negative geoidal heights, two issues complicate the compilation and subsequently the inversion of the BT disturbance. The first is associated with the evaluation of normal gravity below the surface of the reference ellipsoid (RE). The latter is tied to the legitimacy of the harmonic continuation of the BT disturbance in these regions. These two issues are proposed to be resolved by the so called reference quasi-ellipsoid (RQE) approach. New bathymetric and topographic corrections are derived based on the RQE and the inverse problem is formulated based on the RQE. The RQE approach enables the computation of normal gravity by means of the international gravity formula, and makes the harmonic continuation in the regions of negative heights of gravity stations legitimate. The gravimetric inversion is then transformed from the RQE approach back to the RE approach, following the now legitimate harmonic upward continuation of the gravity data to stations on or above the RE. Stripping, the removal of an effect of a known density contrast, is considered in the context of the RQE approach. A numerical case study is presented for the RQE approach in a region of NW Canada.  相似文献   

4.
Discussion of Mean Gravity Along the Plumbline   总被引:2,自引:1,他引:1  
According to the definition of the orthometric height, the mean value of gravity along the plumbline between the Earth's surface and the geoid is defined in an integral sense. In Helmert's (1890) definition of the orthometric height, a linear change of the gravity with depth is assumed. The mean gravity is determined so that the observed gravity at the Earth's surface is reduced to the approximate mid-point of the plumbline using Poincaré-Prey's gravity gradient. Niethammer (1932) and later Mader (1954) took into account the mean value of the gravimetric terrain correction within the topography considering the constant topographical density distribution along the plumbline (for more details see Heiskanen and Moritz, 1967). Vaníek et al. (1995) included the effect of the lateral variation of the topographical density into the definition of Helmert's orthometric height. Recently, Hwang and Hsiao (2003) discussed the influence of the vertical gradient of disturbing gravity on the orthometric heights. In this paper, the mean integral value of gravity along the plumbline within the topography is defined so that the actual topographical density distribution and the change of the disturbing gravity with depth are taken into account. Based on the definition of the mean gravity, the relation between the orthometric and normal heights is discussed.  相似文献   

5.
储层重力密度反演后验约束正则化方法   总被引:2,自引:1,他引:1       下载免费PDF全文
本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果.  相似文献   

6.
于磊  张健  陈石  董淼  徐长仪 《地震学报》2015,37(4):565-574
珍贝—黄岩海山链作为我国南海的残留扩张中心, 对其研究具有重要的科学意义. 本文运用均衡学方法, 通过重力异常数据反演了过珍贝—黄岩海山链剖面的地壳界面变化, 同时计算了岩石圈热结构状态, 在此基础上建立了珍贝—黄岩海山链的岩石圈地温结构模型. 通过均衡分析方法, 对剖面上测点的海底地形数据进行了热均衡和重力均衡分析, 得到了热均衡和重力均衡形变量. 结果表明, 在珍贝—黄岩海山链高热流区域, 热均衡作用可以产生最大约0.55 km的形变, 其重力均衡形变范围为0.77—1.89 km. 热均衡通过改变海底地形和地壳物质密度不断作用于重力均衡, 重力又反过来作用于热均衡, 形成了热均衡-重力均衡动态调节机制.   相似文献   

7.
Based on the line integral (LI) and maximum difference reduction (MDR) methods, an automated iterative forward modelling scheme (LI‐MDR algorithm) is developed for the inversion of 2D bedrock topography from a gravity anomaly profile for heterogeneous sedimentary basins. The unknown basin topography can be smooth as for intracratonic basins or discontinuous as for rift and strike‐slip basins. In case studies using synthetic data, the new algorithm can invert the sedimentary basins bedrock depth within a mean accuracy better than 5% when the gravity anomaly data have an accuracy of better than 0.5 mGal. The main characteristics of the inversion algorithm include: (1) the density contrast of sedimentary basins can be constant or vary horizontally and/or vertically in a very broad but a priori known manner; (2) three inputs are required: the measured gravity anomaly, accuracy level and the density contrast function, (3) the simplification that each gravity station has only one bedrock depth leads to an approach to perform rapid inversions using the forward modelling calculated by LI. The inversion process stops when the residual anomalies (the observed minus the calculated) falls within an ‘error envelope’ whose amplitude is the input accuracy level. The inversion algorithm offers in many cases the possibility of performing an agile 2D gravity inversion on basins with heterogeneous sediments. Both smooth and discontinuous bedrock topography with steep spatial gradients can be well recovered. Limitations include: (1) for each station position, there is only one corresponding point vertically down at the basement; and (2) the largest error in inverting bedrock topography occurs at the deepest points.  相似文献   

8.
Explicit formula for the geoid-quasigeoid separation   总被引:1,自引:0,他引:1  
The explicit formula for the geoid-to-quasigeoid correction is derived in this paper. On comparing the geoidal height and height anomaly, this correction is found to be a function of the mean value of gravity disturbance along the plumbline within the topography. To evaluate the mean gravity disturbance, the gravity field of the Earth is decomposed into components generated by masses within the geoid, topography and atmosphere. Newton’s integration is then used for the computation of topography-and atmosphere-generated components of the mean gravity, while the combined solution for the downward continuation of gravity anomalies and Stokes’ boundary-value problem is utilized in computing the component of mean gravity disturbance generated by mass irregularities within the geoid. On application of this explicit formulism a theoretical accuracy of a few millimetres can be achieved in evaluation of the geoid-to-quasigeoid correction. However, the real accuracy could be lower due to deficiencies within the numerical methods and to errors within the input data (digital terrain and density models and gravity observations).  相似文献   

9.
Introduction Since the middle of the century, gravitational isostasy has been a fundamental hypothesis for inverting the gravity data to find the crust thickness. Geophysicists have done a lot of researches on using gravity data to investigate the depth of Moho discontinuity. Since 1980, the International Lithosphere Program emphasized the importance of investigating the Moho depth variation. Thereafter a lot of results have been published in the world (Braitenberg et al, 2000; Kaban et al,…  相似文献   

10.
The reduction of gravity-field related quantities (e.g., gravity anomalies, geoid heights) due to the topography plays a crucial role in both geodetic and geophysical applications, since in the former it is an intermediate step towards geoid prediction and in the latter it reveals lateral as well as radial density contrasts and infers the geology of the area under study. The computations are usually carried out by employing a DTM and/or a DBM, which describe the topography and bathymetry, respectively. Errors in these DTMs/DBMs will introduce errors in the computed topographic effects, while poor spatial resolution of the topography and bathymetry models will result in aliasing effects to both gravity anomalies and geoid heights, both influencing the accuracy of the estimated solutions. The scope of this work is twofold. First, a validation and accuracy assessment of the SRTM 3″ (90 m) DTM over Greece is performed through comparisons with existing global models as well as with the Greek 450 m national DTMs. Whenever a misrepresentation of the topography is identified in the SRTM data, it is “corrected” using the local 450 m DTM. This process resulted in an improved SRTM DTM called SRTMGr, which was then used to determine terrain effects to gravity field quantities. From the fine-resolution SRTMGr DTMs, coarser models of 15″, 30″, 1′, 2′ and 5′ have been determined in order to investigate aliasing effects on both gravity anomalies and geoid heights by computing terrain effects at variable spatial resolutions. From the results acquired in two test areas, it was concluded that SRTMGr provides similar results to the local DTM making the use of other older global DTMs obsolete. The study for terrain aliasing effects proved that when high-resolution and accuracy gravity and geoid models are needed, then the highest possible resolution DTM should be employed to compute the respective terrain effects. Based on the results acquired from two the test areas a corrected SRTMGr DTM has been compiled for the entire Greek territory towards the development of a new gravimetric geoid model. Results from that analysis are presented based on the well-known remove-compute-restore method, employing land and marine gravity data, EGM08 as a reference geopotential model and the SRTMGr DTM for the computation of the RTM effects.  相似文献   

11.
Gravity data inversion can provide valuable information on the structure of the underlying distribution of mass. The solution of the inversion of gravity data is an ill-posed problem, and many methods have been proposed for solving it using various systematic techniques. The method proposed here is a new approach based on the collocation principle, derived from the Wiener filtering and prediction theory. The natural multiplicity of the solution of the inverse gravimetric problem can be overcome only by assuming a substantially simplified model, in this case a two-layer model, i.e. with one separation surface and one density contrast only. The presence of gravity disturbance and/or outliers in the upper layer is also taken into account. The basic idea of the method is to propagate the covariance structure of the depth function of the separation surface to the covariance structure of the gravity field measured on a reference plane. This can be done since the gravity field produced by the layers is a functional (linearized) of the depth. Furthermore, in this approach, it is possible to obtain the variance of the estimation error which indicates the precision of the computed solution. The method has proved to be effective on simulated data, fulfilling the a priori hypotheses. In real cases which display the required statistical homogeneity, good preliminary solutions, useful for a further quantitative interpretation, have also been derived. A case study is discussed.  相似文献   

12.
常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好.  相似文献   

13.
The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.  相似文献   

14.
Power spectrum analysis of the Bouger gravity values in the Eastern Alps suggests that the gravity field may be separated into long and short wavelength components. The long wavelength component is assumed to be caused by Alpine crustal thickening. This long wavelength component was subjected to a gravimetric single density-interface inversion procedure, giving a gravimetric Mohorovi?i? model which is generally in good agreement with Moho-depths derived by refraction and reflection seismology. The residual high-frequency gravity component correlates well with the main surface geological units in the Eastern Alps. Apparent density mapping by applying an inverse density deconvolution filter to the short wavelength gravity component gives density values for the upper crust which correspond well with averaged density values from rock samples.  相似文献   

15.
高精度金星重力场的获取,是金星探测的重要内容.本文利用最新的金星地形和重力模型,通过高通滤波后的残差地形(RTM)并在考虑均衡改正的情况下改进了重力的短波成分,最终提出了一个新的金星重力模型VGM2013,该模型赤道分辨率达10 km量级,大大高于现有的金星重力场模型,最终结果是金星表面重力加速度和重力扰动.研究中同时发现金星在Airy-Heiskanen均衡模型下的全球最优补偿深度为30 km,金星地壳的密度可能小于当前认为的2700~2900 kg·m-3.VGM2013模型的结果可为将来的金星探测器定轨和着陆导航提供参考,作为重力计算的先验模型.但由于该模型没有包含短波重力观测信息,不建议直接用于更小尺度的地质和地球物理解释.  相似文献   

16.
In mapping the topography of the basement of deep sedimentary basins by gravity modelling, the accuracy can be improved by incorporating an exponential increase in density with depth. For calculating the gravity effect of a three-dimensional (3D) structure with such an exponential density-depth relation a frequency-domain forward algorithm based on series expansion is presented, the numerical evaluation of which can be performed efficiently by fast Fourier transform. The algorithm can be applied in a recursive procedure to give the inverse solution in terms of basement relief. The inversion procedure is satisfactorily tested on a 2D synthetic example and a 3D field example of gravity data from the western margin of the Pannonian Basin in eastern Austria, where up to 2.2 km of Tertiary sediments overlie an igneous or metamorphic basement. The results are confirmed by basement intersections in several wells.  相似文献   

17.
This paper deals with the analysis of gravity anomaly and precise levelling in conjunction with GPS-Levelling data for the computation of a gravimetric geoid and an estimate of the height system bias parameter No for the vertical datum in Pakistan by means of least squares collocation technique. The long term objective is to obtain a regional geoid (or quasi-geoid) modeling using a combination of local data with a high degree and order Earth gravity model (EGM) and to determine a bias (if there is one) with respect to a global mean sea surface. An application of collocation with the optimal covariance parameters has facilitated to achieve gravimetric height anomalies in a global geocentric datum. Residual terrain modeling (RTM) technique has been used in combination with the EGM96 for the reduction and smoothing of the gravity data. A value for the bias parameter No has been estimated with reference to the local GPS-Levelling datum that appears to be 0.705 m with 0.07 m mean square error. The gravimetric height anomalies were compared with height anomalies obtained from GPS-Levelling stations using least square collocation with and without bias adjustment. The bias adjustment minimizes the difference between the gravimetric height anomalies with respect to residual GPS-Levelling data and the standard deviation of the differences drops from 35 cm to 2.6 cm. The results of this study suggest that No adjustment may be a good alternative for the fitting of the final gravimetric geoid as is generally done when using FFT methods.  相似文献   

18.
融合多源数据的高精度、高分辨率的局部重力场建模是物理大地测量学的前沿和热点问题.本文研究了基于径向基函数融合多源数据的局部重力场建模方法,利用Monte-Carlo方差分量估计实现了不同类型的观测数据的合理定权,引入了最小标准差法确定基函数的适宜网络,分析了地形因素对于基函数网络确定及局部重力场建模精度的影响.以泊松小波基函数为构造基函数,结合残差地形模型,融合实测的陆地重力异常、船载重力异常及航空重力扰动数据构建了局部区域陆海统一的似大地水准面模型.研究结果表明:引入残差地形模型平滑了地形质量引入的高频扰动信号,简化了基函数的网络设计;并提高了重力似大地水准面的精度,平原地区其精度提高了4mm,地形起伏较大的山区其精度提高了约5cm.总体而言,基于"三步法"构建的局部重力似大地水准面在荷兰、比利时及德国相关区域,其精度分别达到1.12cm、2.80cm以及2.92cm.  相似文献   

19.
密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.  相似文献   

20.
为评估全球潮汐模型在我国潮汐改正中的适用性,本文首先对10个重力站2016—2018年的观测数据进行了精度评定,而后基于均方根、和方根、纬度依赖关系以及重力残差等指标对7个全球潮汐模型进行了精度评定.结果表明:10个重力站的一些评价指标达到甚至超越了早期超导重力仪,例如M2波潮汐因子的中误差普遍小于0.00070,其中...  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号