煤炭资源勘查误差理论研究

对煤炭资源勘查资源量与客观物质量之间差值的误差理论探讨,提出误差由真值误差、系统误差和离差3个基本部分组成,其中真值误差是不可消除的观测误差、系统误差是控制网度所决定的准确度标准、而离差是由地质构造及矿床形态的多解性产生的资源量波动。在此基础上,构筑了一个关于误差的数学模型,最后应用误差理论对现行勘查工作中常见的错误进行分析,并就如何避免错误及减少误差提出了见解。     

随机Dirichlet问题解的调和展开

针对重力学随机Dirichlet问题,通过适当地对边界检验函数的分解,并在随机边界样本空间中提取确定性部分的对偶基,本文将随机Dirichlet问题的一般解展开为一随机系数的调和级数形式。     

辽宁铁岭市土壤侵蚀时空演变研究

近年来,人为活动的加剧造成铁岭地区水土流失日益严重。本文采用遥感和GIS技术,以部颁标准SL190-96对铁岭地区1991年和2001年的土壤侵蚀情况进行研究,并结合1986年和1996年的土壤侵蚀调查数据,分析了铁岭地区近15年来土壤侵蚀强度的时空变化,结果表明铁岭地区在此期间土壤侵蚀面积缓慢增长,其中强度侵蚀面积在1986~1996年间持续增长,2001年明显下降;中度侵蚀面积在1986~1996年间持续下降,2001年则明显增强;空间分析显示明显恶化的区域有西丰、开原和清河区,具明显改善的为铁岭县,而昌图和调兵山地区基本保持稳定。今后该区水保的重点应放在强度侵蚀的治理,以及土壤侵蚀等级增强的区域。     

Tectonic relation between northeastern China and the Korean peninsula revealed by interpretation of GRACE satellite gravity data

The major continental blocks in northeastern Asia are the North China block and the South China block, which have collided starting from the Korean peninsula. Geologic and geophysical interpretations reveal a well defined suture zone in northeastern China from Qinling through Dabie to Jiaodong. The discovery of high-pressure metamorphic rocks in the Hongseong area of the Korean peninsula, prominent evidence for the collision zone, indicates extension of the collision zone in northeastern China into the Korean peninsula. Interpretation of the GRACE satellite gravity dataset shows two prominent structural boundaries in the Yellow Sea. One extends from the Jiaodong Belt in eastern China to the Imjingang Belt in the Korean peninsula. The other extends from near Nanjing, eastern China, to Hongseong. Tectonic movement in or near the suture zone may be responsible for seismic activity in the western Korean peninsula and the development of the Yellow Sea sedimentary basin.     

The concentration and distribution of different mercury species in the water columns and sediment of Aha Lake

In order to find out whether Aha Lake was polluted by the acid mining waste water or not, the concentration and distribution of different mercuryspecies in the water columns and sediment profile collected from Aha Lake were investigated. It was found that discernible seasonal variation of different mercury species in water body were obtained in the Aha Reservoir. With regards to the whole sampling periods, the concentrations of HgP in the Aha Reservoir water body were evidently correlated to the concentrations of total mercury, showing that total mercury was mostly associated with particle mercury. The concentrations of methylmercury in water body were also evidently correlated to the concentrations of dissolved mercury. The dissolved mercury evidently affects the distribution and transportation of methylmercury. However, there is no correlation between methylmercury and total mercury. The dissolved mercury, reactive mercury, dissolved methylmercury levels in the water body of high flow period were much higher than those in low flow period. The distribution, speciation and levels of mercury within the Aha Reservoir water body were governed by several factors, such as the output of river, the release of sediment . Discernible seasonal variation of total mercury and methylmercury in porewater was described during the sampling periods, with the concentrations in high flow period generally higher than those in low flow period. The methylmercury in pore water column was evidently correlated to that of the sediment. The results indicated that highly elevated MeHgD concentrations in the porewater were produced at the depths from 2 to 5 cm in the sediment profile, and decreased sharply with depth. A positive correlation has been found between MeHgD formation and sulfate reducing bacterial activity. These highly elevated concentrations of MeHgD at the intersurface between waters and sediments suggest a favorable methylation condition. Moreover,     

A comparative error analysis of manual versus automated methods of data acquisition for algebraic strain estimation

This paper presents a statistical analysis of the algebraic strain estimation algorithm of Shimamoto and Ikeda [Shimamoto, T., Ikeda, Y., 1976. A simple algebraic method for strain estimation from deformed eillipsoidal objects: 1. Basic theory. Tectonophysics 36, 315–337]. It is argued that the error in their strain estimation procedure can be quantified using an expected discrepancy measure. An analysis of this measure demonstrates that the error is inversely proportional to the number of clasts used. The paper also examines the role of measurement error, in particular that incurred under (i) a moment based and (ii) manual data acquisition methods. Detailed analysis of these two acquisition methods shows that in both cases, the effect of measurement error on the expected discrepancy is small relative to the effect of the sample size (number of objects). Given their relative speed advantage, this result favours the use of automated measurement methods even if they incur more measurement error on individual objects. Validation of these results is carried out by means of a simulation study, as well as by reference to studies appearing in previous literature. The results are also applied to obtain an upper bound on the error of strain estimation for various studies published in the literature on strain analysis.     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.