人工神经网络在白云岩类型识别中的应用

利用地球化学数据,运用人工神经网络方法对美国密苏里州东南Ｂｏｎｎｅｔｅｒｒｅ组（寒武纪）滨海相的白云岩进行了分类、识别,判别率达１００％,结果表明,该方法性能良好,可望成为岩石分类、判别的有效手段。     

Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

区域地下水信息管理系统研究

本文通过《河南省地下水动态信息管理系统》的建立,讨论了区域地下水信息管理系统的设计目标,数据结构和软件实现,着重介绍了系统的软件设计特点,最后指出了这种管理系统在水文地质部门中具有广泛的应用前景。     

A realistic model for point-sources imaged on array detectors: The model and initial results

We have constructed a computer model for simulation of point-sources imaged on two-dimensional detectors. An attempt has been made to ensure that the model produces data that mimic real data taken with 2-D detectors. To be realistic, such simulations must include randomly generated noise of the appropriate type from all sources (e.g. source, background, and detector). The model is generic and accepts input values for parameters such as pixel size, read noise, source magnitude, and sky brightness. Point-source profiles are then generated with noise and detector characteristics added via our model. The synthetic data are output as simple integrations (onedimensional), as radial slices (two-dimensional), and as intensity-contour plots (three-dimensional). Each noise source can be turned on or off so that they can be studied separately as well as in combination to yield a realistic view of an image. This paper presents the basic properties of the model and some examples of how it can be used to simulate the effects of changing image position, image scale, signal strength, noise characteristics, and data reduction procedures.Use of the model has allowed us to confirm and quantify three points: 1) The use of traditionalsize apertures for photometry of faint point-sources adds substantial noise to the measurement which can significantly degrade the quality of the observation; 2) The number of pixels used to estimate the background is important and must be considered when estimating errors; and 3) The CCD equation normally used by the astronomical community consistently overestimates the signal-to-noise obtainable by a measurement while a revised equation, discussed here, provides a better estimator.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

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.     

用于湿地气候效应模拟的三江平原下垫面数据获取

三江平原是我国沼泽湿地的重要分布区。在20世纪的后50 a中,三江平原在经过几次大规模的开荒后,其下垫面发生了明显的阶段性变化,湿地大面积消失,农田面积不断增加。为了认识20世纪下半叶不同时期三江平原湿地的气候效应,拟采用第三代区域气候模式(RegCM3)来开展研究。而目前被广泛用于RegCM3的GL-CC数据,所反映的研究区下垫面几乎不存在沼泽湿地,这与事实严重不符,所以必须获取三江平原有沼泽湿地存在的真实的下垫面数据,才能得到可靠的湿地气候效应的模拟结果。分别选择了4期(1954年、1985年、1996年和2000年)有阶段代表性的航片、TM影像数据为主要数据源,利用ERDAS8.7和Arcview3.2等软件,首先得到4期研究区的土地利用数据,然后通过土地利用数据类型向模式数据类型的转换,从4期研究区土地利用数据中提取出了三江平原下垫面的模式数据。经过误差分析,认为获取的4期三江平原下垫面模式数据准确的反映出研究区几个时期的湿地下垫面实况,可以作为三江平原湿地气候效应模拟的基础数据。     

Improved reliability in the interpretation of geochemical measurements by the quantification of uncertainty from sampling

All geochemical measurements require the taking of field samples, but the uncertainty that this process causes is often ignored when assessing the reliability of the interpretation, of the geochemistry or the health implications. Recently devised methods for the estimation, optimisation and reduction of this uncertainty have been evaluated by their application to the investigation of contaminated land. Uncertainty of measurement caused by primary sampling has been estimated for a range of six different contaminated land site investigations, using an increasingly recognized procedure. These site investigations were selected to reflect a wide range of different sizes, contaminants （organic and metals）, previous land uses （e.g. tin mining, railway sidings and gas works）, intended future use （housing to nature reserves） and routinely applied sampling methods. The results showed that the uncertainty on measurements was substantial, ranging from 25% to 186% of the concentration values at the different sites. Sampling was identified as the dominant source of the uncertainty （〉70% of measurement uncertainty） in most cases. The fitness-for-purpose of the measurements was judged using the optimized contaminated land investigation （OCLI） method. This identifies the optimal level of uncertainty that reduces to overall financial loss caused by the measurement procedures and the misclassification of the contamination, caused by the uncertainty. Generally the uncertainty of the actual measurements made in these different site investigations was found to be sub-optimal, and too large by a factor of approximately two. The uncertainty is usually limited by the sampling, but this can be reduced by increasing the sample mass by a factor of 4 （predicted by sampling theory）. It is concluded that knowing the value of the uncertainty enables the interpretation to be made more reliable, and that sampling is the main factor limiting most investigations. This new approach quantifies this problem for the first time, and allows sampling procedures to be critically evaluated, and modified, to improve the reliability of the geochemical assessment.     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.