地球化学方法在福安堡钨钼矿床的应用

吉林省中部大面积森林覆盖区,分布着大面积华力西期、印支期花岗岩类,受传统成矿理论的束缚,这一地区成为有色金属矿产勘查的"处女地"。1:20万化探扫面之后,相继开展了不同方法的地球化学异常查证工作,发现了很有前景的福安堡钨钼矿床。结果表明,在东北森林覆盖区化探找矿方法十分有效。     

江西金山金矿田东蒋金矿床地质特征及找矿远景

阐述了东蒋金矿床的成矿地质背景,通过对该金矿区地质特征的分析和总结,认为东蒋韧性剪切带为一卓具规模的控矿、容矿构造,其具有构造岩分带、蚀变分带及金矿化强度分带“三位一体”的规律性分布特点。指出该区业已发现韧性剪切带中的似层状金矿体和陡倾斜断裂带(F1)中的脉状金矿体,而该带目前无矿的东段分布有物探激电、航磁、重力异常及化探异常,这些都显示该区尚具有良好的找矿前景。     

Ellipsoidal geoid computation

Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e². They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.     

A technique for modeling the Earth’s gravity field on the basis of satellite accelerations

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

A very large scale GIS-based groundwater flow model for the Nubian sandstone aquifer in Eastern Sahara (Egypt,northern Sudan and eastern Libya)

A three-dimensional GIS-based groundwater flow model for the Nubian Sandstone Aquifer in the eastern Sahara was developed and calibrated under steady-state and transient conditions. The model was used to simulate the response of the aquifer to climatic changes that occurred during the last 25,000 years. The simulation results indicated that the groundwater in this aquifer was formed by infiltration during the wet periods 20,000 and 5,000 years b.p. The recharge of groundwater due to regional groundwater flow from more humid areas in the south was excluded. It also indicates that the Nubian Aquifer System is a fossil aquifer, which had been in an unsteady state condition for the last 3,000 years.

---

Resumen Fue desarrollado un modelo de flujo de agua subterránea en tres dimensiones, basado en un SIG, para el Acuífero Arenisca Nubian en el Sahara Oriental, el cual fue calibrado para condiciones de estado estacionario y transitorio. El modelo se usó para simular la respuesta del acuífero a los cambios climáticos que ocurrieron durante los últimos 25000 años. Los resultados de esta simulación indicaron que el agua subterránea en este acuífero, se formó por infiltración, durante los períodos húmedos que hubo hace 20000 y 5000 años, antes del presente. Fue excluida la recarga del acuífero debida a un flujo regional de agua subterránea proveniente de áreas con un clima más húmedo en el sur. El modelo también muestra, que el Sistema Acuífero Nubian es un acuífero fósil, el cual ha permanecido en una condición de estado no estacionario, durante los últimos tres mil años.

---

Résumé Pour laquifère gréseux Nubien de Sahara -Est on a mis au points un modèle tridimensionnel, basé sur GIS. Le modèle a été calibré tant pour lécoulement stationnaire que pour lécoulement transitoire. On a simulé après la réponse de laquifère aux changements climatiques des derniers 25000 ans. Les résultats des simulations indiquent que la nappe a été rechargée par des infiltrations pendant une période humide qui sétend 5000 et 20000 ans, dès temps actuel. On na pas pris en compte la recharge de laquifère par la zone plus humide située dans sa partie sud. Le modèle indique aussi que leau de laquifère Nubien est une eau fossile qui a eu un écoulement transitoire pendant les derniers 3000 ans.

     

Non-Continuation of Integrals of the Rotating Two-body Problem in Hill's Lunar Problem

We consider Hill's lunar problem as a perturbation of the integrable two-body problem. For this we avoid the usual normalization in which the angular velocity of the rotating frame of reference is put equal to unity and consider as the perturbation parameter. We first express the Hamiltonian H of Hill's lunar problem in the Delaunay variables. More precisely we deduce the expressions of H along the orbits of the two-body problem. Afterwards with the help of the conserved quantities of the planar two-body problem (energy, angular momentum and Laplace–Runge–Lenz vector) we prove that Hill's lunar problem does not possess a second integral of motion, independent of H, in the sense that there exist no analytic continuation of integrals, which are linear functions of in the rotating two-body problem. In connection with the proof of this main result we give a further restrictive statement to the nonintegrability of Hill's lunar problem.     

中国-日本绝对重力仪器测量比对结果

本文给出了中科院测地所FG5-112和日本京都大学FG5-201及日本国土地理院FG5-104等几台绝对重力仪在武汉九峰和香港天文台绝对重力测量比对结果。分析说明不同仪器间具有良好的一致性,其中FG5-112与FG5-104结果间的互差达8×10-8m/s2;FG5-112与FG5-201和FG5-104等的观测互差小于2×10-8m/s2。这种比对观测,可为深入了解仪器性能和精确测定不同台站重力非潮汐变化提供可靠参考。     

航空重力测量数据向下延拓中空间协方差函数特性研究

本文通过对航空重力测量数据的分析,建立起具有空间相关特性的空间协方差函数模型。利用所建立的模型,将航空重力测量试验数据采用最小二乘配置的方法向下延拓到地面,与地面已知重力数据进行比较。由比较结果,分析了空间协方差函数模型对航空重力测量数据向下延拓结果的影响,并表明所建立的空间协方差函数模型,满足协方差函数特性,且在局部区域使用能够达到满意的精度要求。     

卫星重力梯度测量数据的粗差探测

分别讨论了阈值法、Dixon检验法和小波算法应用于卫星重力梯度数据粗差探测的效果。为了克服这3种单一粗差探测方法的不足,提出了卫星重力梯度数据的粗差探测组合方案,即联合Dixon检验和小波算法,以及联合Dixon检验、阈值法和小波算法的组合方案。模拟结果表明,两种组合方案均能有效改善粗差探测的效果。     

绝对重力仪对日本海域6.0级地震的精确纪录

FG5绝对重力仪在上海佘山基准站观测前后发生了日本海域6.0级地震。数据处理分析表明,地震在该基准站引起的绝对重力值变化约为-16微伽,这对地震进行预报以及对地震造成的危害性进行评估有着积极意义。