上宫构造蚀变岩型金矿床地质特征及化探找金经验

上宫金矿位于河南省西部熊耳山北麓,在该区1 020 km²范围内进行了1:5万的地球化学水系沉积物测量,按水系划分取样单元,取样本着"小沟多取,支流放稀,大河不管"的原则进行。本次水系沉积物测量的特点是:取样少,控制面积大,速度快,省钱、省人、省物,找矿效果好、经济效益好。本区水系沉积物测量中所发现的金异常,后经异常检查、地质普查评价及钻探验证,确定为一大型构造蚀变岩型金矿床。     

直流输电线路对磁测的影响

越来越多的远距离、大容量直流输电线路投入运行对线路附近地区的磁测工作产生了不可忽视的影响。本文阐述了直流输电线路的结构、线路磁场特征及对以地质找矿为目的的磁测工作的影响。对直流输电线路磁场在磁测图中的特点也做了讨论,列举了直流输电线路影响的实例。     

Magnetite-Fluorite Rock: A New Rock Type of Hot Water Sedimentation

The new type hot water sedimentary rock - magnetite-fluorite rock occurs as quasi-layers inflat parts of contact zones between rock body and strata in Bamiaushan of Changshan County,Zhejiang Province, China. The main mineral assemblage is fluorite magnetite cassiterite. The rockshows typical laminated structure and obvious mosaic texture. Its formation temperature is between123℃-160℃, averaging at 142℃. The major chemical composition of the rock includes CaF2, SiO2,Al2O3, FeO, and Fe2O3; the high-content microelement association includes W, Sn, Be, Rb, Sr, S, andCl; and the total content of REE is low (∑REE between 35.34x10-6-38.35x10-6), showing LREEenrichment type of distribution pattern. Diagenesis: driven by the tectonic stress, the formation waterheated in the deep strata had moved along the fissures or fractures in strata and had extractedcomponents from the strata on the way, and finally stagnated in the flat parts of contact zones betweenrock body and strata. With drop in temperature, magnetite and fluorite were separated from the hotwater and precipitated alternately, forming this hot water sedimentary rock with new typemineralogical composition, typical laminated structure, obvious mosaic texture and sub-horizontaloccurrence. The characteristics of the new type mineralogical composition, sedimentary tectonicenvironment and chemical composition are different from that of the well-known traditionalhydrothermai sedimentary rocks.     

下庄矿田湖子地区铀矿找矿方向探讨

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

祁阳山字型构造质疑

结合区域构造演化背景及构造变形特征,对祁阳弧形构造的形成机制进行研究。分析表明,传统山字型构造机制难以解释祁阳弧形构造的诸多特征,且前人认为变形动力来源于东西向挤压应力场,也与印支期区域构造应力场实为NWW向的构造背景不符,因此祁阳弧形构造并不属典型山字型构造。基于研究区构造活动的客观实际,提出祁阳弧形构造的可能形成机制:NW向基底隐伏断裂和NNE-NE向主干断裂分别于印支期和燕山期产生强烈左旋走滑活动,从而使区域NNE向构造线在中段产生左旋偏转成为NNW向,从而形成S形的祁阳弧形构造。这一机制可较好解释祁阳弧形构造的若干特征,如关帝庙穹窿呈NWW走向、北反射弧构造形迹不显著、北弧弧顶脱位及内弧曲率大于外弧、南弧弧顶脱位、紫云-中田-高峰串珠状穹窿的形成等。上述认识对华南地区弧形构造研究具有一定启示意义:除山字型构造作用和砥柱作用外,还应注意断裂走滑等其他构造活动对弧形构造的制约,以及多次构造活动叠加对构造形迹可能造成的影响。     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

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

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.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

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

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.     

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

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