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.     

预测陕西关中地区需水量的改进GM(1,1)模型

需水量评价与预测是水资源规划与管理一项重要的研究内容。本文采用定额法计算需水量。通过对陕西关中地区需水量要素组成的分析,认为需水量要素属平稳时间序列,且一阶累加生成数据系列满足指数规律,符合灰色预测条件。为了提高GM(1,1)模型的精度,采用一种改变背景值的新方法,即:中心逼近方法。通过精度检验,中心逼近式GM(1,1)模型平均误差百分比、误差平方和分别比传统的GM(1,1)模型提高了0.147和1.579。用中心逼近式GM(1,1)模型预测需水量各要素近期和中长期值,经检验,其中非农业人口、农业人口、耕地面积、工业总产值预测模型满足一级模型,等级为良好。牲畜头数预测模型为不合格模型,分析其原因,主要是因为时间序列数据不完全符合指数规律。通过定额法计算,预计75%保证率时2005年和2010年关中地区工农业需水量分别占总需水量的92.21%和89.75%,生活需水量分别占7%和9.04%。     

复合土钉墙支护技术在某大型软土基坑围护中的应用

根据拟建建筑基坑的工程地质条件、周边环境条件、基坑开挖深度等特点综合考虑，对基坑分段采用超前微型桩支护加土钉墙、放坡加土钉墙2种复合支护方案。实践证明，该工程所采用的多种复合式支护实施方案，技术安全可靠，经济合理。     

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

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

试论洋山深水港与上海国际航运中心建设

建设国际航运中心是上海城市发展的宏伟目标之一。目前,围绕东北亚国际航运中心的竞争十分激烈,釜山、神户、高雄等港口均提出了明确的发展目标,并纷纷建设一系列深水集装箱泊位。本文以国际竞争态势为背景,着重分析了国际集装箱航运市场的最新动态、洋山深水港对建设上海国际航运中心的重要意义,并详细论证了强化上海国际航运中心功能的政策创新,探讨了优化上海国际航运中心营运水平的关键措施,为推进上海国际航运中心建设献计献策。     

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.