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.     

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.     

Dissolved trace elements in river water: spatial distribution and the influencing factor, a study for the Pearl River Delta Economic Zone, China

Twenty-nine water samples were collected from different river channels of the Pearl River Delta Economic Zone, China. An inductively coupled plasma-mass spectromonitor (ICP-MS) was used to measure concentrations of the trace elements in these samples. The results suggest that the average concentrations of rare earth elements in river water show an increasing trend from the West River, the North River, the rivers of the Pearl River Delta, and the Shenzhen River to the East River. Relatively high concentrations of heavy metals appear in the East River, the rivers of the Pearl River Delta and the Shenzhen River, while the West River and the North River have relatively low heavy metal concentrations. Trace element concentrations in samples collected near urban or industrial areas are much higher than those of samples collected from distant areas, away from urban and industrial areas. After natural conditions, human activities have significant influence on the trace element concentrations in river water. This trace element concentration’s spatial distribution in the river water from the Pearl River Delta Economic Zone is actually an integrated effect of natural conditions and human activity.     

要十分重视另一类型地质灾害的防治

文章提出另一类型地质灾害,即化学型地质灾害的定义及其分类。并认为化学型地质灾害不像物理型地质灾害那样发生瞬时灾害,而有一个延缓过程,但是一旦发生,难于挽救。化学型地质灾害体的主体是化学元素,开展地球化学调查,研究化学元素的赋存规律、土壤及沉积物承受各类元素的极限、各类元素活化的条件是揭示诱发化学型地质灾害机理的关键。因此要十分重视与加强化学型地质灾害的调查、研究、监测与治理工作。     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

Petrogenesis and geotectonic setting of the PanAfrican basement rocks in Bamenda Massif, Obudu Plateau, southeastern Nigeria： Evidence from trace element geochemistry

1 IntroductionThe study area is part of the Obudu Plateau,which forms part of the Precambrian Pan-African tec-tonothermal belt, lying between the West African Cra-ton and the Gabon-Congo Craton ( Fig. 1). It isbounded by the Benue Trough, in the NE-SW axi…     

REE/trace element characteristics of sandstone-type uranium deposits in the Ordos Basin

1 Introduction The Ordos Basin is the second largest sedimentary basin in China. During the last 10 years, a great progress has been achieved in the aspects of tectonic evolution, dynamics process, inner and outer geological processes during Mesozoic-Cen…     

中国大陆科学钻探工程主孔100～2000 m 放射性产热元素的垂向分布特征及其成因

地球内部放射性产热元素U、Th及K(HPE)含量如何随深度的变化而变化是固体地球科学中的一个重要参数,在限定地壳的热和流变学结构、地球化学、岩石和构造模型中起关键性的作用。对HPE深部分布的认识主要来自于对大型花岗岩岩基的研究及对地表热流值和产热率之间关系的模拟,对高压-超高压变质地体的HPE随深度的分布缺乏认识。在苏鲁超高压变质带中实施的中国大陆科学钻探工程届时将提供超过5km的岩心,为确定苏鲁超高压带的HPE结构提供了最好的机会。对CCSD的100～2000m岩心的732块样品详细的放射性产热元素含量的测试及现今产热率计算的初步结果表明:(1)花岗质片麻岩具有最高的产热率,平均为1665×10-11W/kg;(2)副片麻岩(变沉积岩)具有中等的产热率,为994×10-11W/kg;(3)金红石榴辉岩及石榴石橄榄岩具有最低的产热率,为17×10-11～20×10-11W/kg;(4)放射性产热元素的浓度及相应的产热率随岩性的变化而变化,呈现阶梯状的分布特征。产热率随深度的变化特征表明CCSD主孔中的HPE三明治结构。该结构可能代表着高压-超高压变质地体中的典型HPE结构,比大型花岗岩岩基的HPE结构更复杂,与传统的HPE指数衰减分布模式不吻合。CCSD主孔中所揭示的三明治式HPE结构是大陆被动陆缘中的基性及超基性岩、沉积岩及花岗岩在深