川西新场气田地气测量试验   总被引：2，自引：0，他引：2  

鲁人齐  王多义  刘亚伟  左银辉  邓美洲  刘应汉 《物探与化探》2008,32(6)

通过对新场气田地气测量分析,发现地气测量的纳米微粒金属元素在新场气田上方有明显的异常显示。综合土壤地球化探结果对比分析,认为Ba、Cd、Cu、Mo、Zn、P、V、Ni等20多种纳米微粒金属元素在新场气田与相同的土壤微量元素都有着共同的响应。同时分析了有机地球化学指标在气田地表上方的表现模式,其特征和地气、土壤元素表现基本一致。地气测量方法在一定条件下对寻找天然气藏、确定气田范围具有一定的应用前景。  相似文献   

下庄矿田湖子地区铀矿找矿方向探讨   总被引：2，自引：0，他引：2       下载免费PDF全文

邱列  毛孟才  王建雄 《地质找矿论丛》2008,23(1):11-15,26

湖子地区位于NNE向诸广-新兴铀成矿带与近EW向大东山-漳州大断裂复合部位,是华南早、晚两期铀成矿热液活动叠加区、放射性高场区,既有"硅化带大脉"型铀矿产出,又有"交点"型铀矿存在,找矿前景良好.文章在论述下庄矿田铀成矿地质环境、铀成矿特征及铀矿定位条件基础上,分析了湖子地区铀成矿条件与找矿前景,指出该区今后铀矿找矿方向是:① 6009号硅化断裂带,在其北段找硅化带大脉型铀矿,在其南段找"交点"型铀矿;②新桥-下庄硅化断裂带和6009号带之间成矿部位,寻找硅化带型和"交点"型铀矿;③NW向、近EW向辉绿岩脉与NNE向、NE向构造交汇部位(交点),寻找"交点"型铀矿.  相似文献   

柴达木盆地三湖地区磁力异常形成机理研究与异常分析   总被引：1，自引：0，他引：1       下载免费PDF全文

王志强  王财富  徐子远  李瑞 《地质找矿论丛》2008,23(4):339-342,347

柴达木盆地的油气资源丰富,近年来在该区发现了大量异常显示,但是各种异常形成机理各不相同,存在较大多解性,使得该区油气储量探明率很低.根据柴达木盆地以往开展的常规非地震勘探的现状及效果,以磁测找油机理为基础,采用磁力资料化极、剩余磁力异常提取和磁力平剖图等处理手段,对三湖地区生物气藏磁异常进行分析研究,取得明显的异常显示效果.  相似文献   

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C   总被引：2，自引：1，他引：1  

Christoph Förste  Roland Schmidt  Richard Stubenvoll  Frank Flechtner  Ulrich Meyer  Rolf König  Hans Neumayer  Richard Biancale  Jean-Michel Lemoine  Sean Bruinsma  Sylvain Loyer  Franz Barthelmes  Saskia Esselborn 《Journal of Geodesy》2008,82(6):331-346

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

Methodology and use of tensor invariants for satellite gravity gradiometry   总被引：2，自引：1，他引：1  

Oliver Baur  Nico Sneeuw  Erik W. Grafarend 《Journal of Geodesy》2008,82(4-5):279-293

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.  相似文献   

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem   总被引：4，自引：1，他引：3  

Maria Clara de Lacy  Mirko Reguzzoni  Fernando Sansò  Giovanna Venuti 《Journal of Geodesy》2008,82(9):527-542

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.  相似文献   

GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals   总被引：19，自引：8，他引：11  

Yanming Feng 《Journal of Geodesy》2008,82(12):847-862

This paper presents a general modeling strategy for ambiguity resolution (AR) and position estimation (PE) using three or more phase-based ranging signals from a global navigation satellite system (GNSS). The proposed strategy will identify three best “virtual” signals to allow for more reliable AR under certain observational conditions characterized by ionospheric and tropospheric delay variability, level of phase noise and orbit accuracy. The selected virtual signals suffer from minimal or relatively low ionospheric effects, and thus are known as ionosphere-reduced virtual signals. As a result, the ionospheric parameters in the geometry-based observational models can be eliminated for long baselines, typically those of length tens to hundreds of kilometres. The proposed modeling comprises three major steps. Step 1 is the geometry-free determination of the extra-widelane (EWL) formed between the two closest L-band carrier measurements, directly from the two corresponding code measurements. Step 2 forms the second EWL signal and resolves the integer ambiguity with a geometry-based estimator alone or together with the first EWL. This is followed by a procedure to correct for the first-order ionospheric delay using the two ambiguity-fixed widelane (WL) signals derived from the integer-fixed EWL signals. Step 3 finds an independent narrow-lane (NL) signal, which is used together with a refined WL to resolve NL ambiguity with geometry-based integer estimation and search algorithms. As a result, the above two AR processes performed with WL/NL and EWL/WL signals respectively, either in sequence or in parallel, can support real time kinematic (RTK) positioning over baselines of tens to hundreds of kilometres, thus enabling centimetre-to-decimentre positioning at the local, regional and even global scales in the future.  相似文献   

A data-driven approach to local gravity field modelling using spherical radial basis functions   总被引：3，自引：0，他引：3  

R. Klees  R. Tenzer  I. Prutkin  T. Wittwer 《Journal of Geodesy》2008,82(8):457-471

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.  相似文献   

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition   总被引：3，自引：0，他引：3  

Oliver Baur  Gerrit Austen  Jürgen Kusche 《Journal of Geodesy》2008,82(4-5):207-221

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.  相似文献   

中国西部现今地壳形变与地震活动特征     

李金平 《测绘工程》2008,17(1):6-9

以1999-2004年间中国西部地区近650个GPS站点的观测资料为基础,采用二维“高张力样条”函数内插算法获得了连续地壳形变场。结果表明,相对于稳定欧亚参考框架,中国西部现今地壳运动西强东弱,北向运动从西向东逐渐减弱,东向运动逐渐增强,青藏高原东缘及附近地区是东向运动的消减区带;主应变从西向东由压缩应变转变为拉伸应变;最大剪应变主要位于喜马拉雅地块及其东部地区;中国西部地区地震活动在空间分布上似乎集中发生在印度板块北东向挤压欧亚板块的两条共扼带内,与地壳形变场有一定对应关系,地壳形变的高应变率区为中强地震多发地带。  相似文献