Application of X-ray computed tomography for analyzing cleat spacing and cleat aperture in coal samples

New techniques to determine distributions of cleat aperture, cleat orientation and cleat spacing from CT scans have been developed. For cleat orientation and spacing distributions, two different coal blocks were scanned. The CT scans have been analyzed for the three orthogonal directions. Histograms of the cleat orientations are bimodal, expressing the typical cleat texture of face and butt cleats and bedding perpendicular relaxation fractures. Deviations up to 20° from the peak values in the cleat orientation distributions were used as input for automated image analysis of cleat spacing. Distributions of the cleat spacing measurements are related to the face and butt cleat directions. The term “relevant cleat length” is introduced as a measure to extract the amount of cleat length involved with the cleat spacing measurements. The ratio ranges from 0.03 to 0.38 and expresses the difference in cleat texture in both samples. Cleat spacing versus relevant cleat length shows sample specific patterns for face cleat, butt cleat and bedding. To describe cleat aperture quantitatively, peak height and missing attenuation have been used. The image of a cleat was seen as a convolution of a rectangular fracture profile with a Gaussian point spread function.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).     

A comparison of chironomid biostratigraphy from Lake Vuolep Njakajaure with vegetation, lake-level, and climate changes in Abisko National Park, Sweden

Chironomid remains from the sediment of Lake Vuolep Njakajaure reflect limnological conditions resulting from changing climate and vegetation throughout the Holocene, but do not strictly follow accepted climate trends or the vegetation history based on regional pollen and macrofossil analyses. Chironomid community changes appear to be influenced by organic nutrient input from the surrounding catchment vegetation and lake hydrology, both of which are indirectly responding to some combination of climate change, hypolimnetic oxygen concentration, and changes in basin morphology. The chironomid-based quantitative mean July air-temperature reconstruction differs from other regional quantitative records; this discrepancy is likely related to limnological conditions particular to Lake Vuolep Njakajaure. Comparison of a northern Swedish temperature transfer function and one from western Canada reveals differences in the mean July air-temperature optima of several common taxa, suggesting that the existing conservative estimates of Holocene climate change in northern Sweden may be underestimated due to the limited temperature gradient captured by the Swedish training set.     

Biological and Geochemical Record of Anthropogenic Impacts in Recent Sediments from Lake Pupuke, Auckland City, New Zealand

Frozen sediment cores from Lake Pupuke in Auckland City, New Zealand, contain a high resolution decadal to annual scale record of changing lake paleoenvironments and geochemistry that reflects changing landuse and hydrology in the catchment over the past c. 190 years. A reliable chronology is available from AMS ¹⁴C and ²¹⁰Pb dating of the sediments, with the timing of the older part of the record supported by the first appearance of pollen of introduced Pinus and Plantago lanceolata associated with European settlement of Auckland in the early 1840s. Diatom stratigraphy, sediment elemental and carbon isotope geochemistry reflect changes in sediment sources and lake conditions commensurate with European development of the Lake Pupuke catchment, in particular enhanced algal productivity controlled by the influx of nutrients after c. 1920 AD. Attempts to prevent nuisance algal blooms in 1933, 1934 and 1939 using CuSO₄ addition produced Cu spikes in the sediment that allowed confirmation of the accuracy of the ²¹⁰Pb chronology. Hence, the elemental and isotopic composition of the Lake Pupuke sediments reflect the timing of significant anthropogenic activities, rather than climatic variations, that have occurred within the watershed over the past c. 190 years. The comparison of records of land use change in the catchment with the multi-proxy record of changes in the sediments usually allowed unambiguous identification of the signatures of change and their causes.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

地学自由曲面径向基函数网络重建的详细算法

径向基函数(RadialBasisFunction,简称RBF)神经网络是一种理想的地学离散数据网格化工具,能够适应各种不同分布形式和边界条件的数据,收敛速度较快,可以逼近任何复杂曲面。这里详细介绍了RBF神经网络的算法。适当径向基函数的形式和偏差系数是使用RBF神经网络作地学曲面重建的关键。大量的实际数据验证结果表明,当选用Gauss型径向基函数时,一般可获得比较理想的网格化效果,同时具备"曲面平滑"和"拟合度高"的特点。