Automated Kerogen Classification in Microscope Images of Dispersed Kerogen Preparation

We develop the classification part of a system that analyses transmitted light microscope images of dispersed kerogen preparation. The system automatically extracts kerogen pieces from the image and labels each piece as either inertinite or vitrinite. The image pre-processing analysis consists of background removal, identification of kerogen material, object segmentation, object extraction (individual images of pieces of kerogen) and feature calculation for each object. An expert palynologist was asked to label the objects into categories inertinite and vitrinite, which provided the ground truth for the classification experiment. Ten state-of-the-art classifiers and classifier ensembles were compared: Naïve Bayes, decision tree, nearest neighbour, the logistic classifier, multilayered perceptron (MLP), support vector machines (SVM), AdaBoost, Bagging, LogitBoost and Random Forest. The logistic classifier was singled out as the most accurate classifier, with an accuracy greater than 90. Using a 10 times 10-fold cross-validation provided within the Weka software, we found that the logistic classifier was significantly better than five classifiers (p<0.05) and indistinguishable from the other four classifiers. The initial set of 32 features was subsequently reduced to 6 features without compromising the classification accuracy. A further evaluation of the system alerted us to the possible sensitivity of the classification to the ground truth that might vary from one human expert to another. The analysis also revealed that the logistic classifier made most of the correct classifications with a high certainty.     

Some Results on the Spatial Breakdown Point of Robust Point Estimates of the Variogram

The effect of outliers on estimates of the variogram depends on how they are distributed in space. The ‘spatial breakdown point’ is the largest proportion of observations which can be drawn from some arbitrary contaminating process without destroying a robust variogram estimator, when they are arranged in the most damaging spatial pattern. A numerical method is presented to find the spatial breakdown point for any sample array in two dimensions or more. It is shown by means of some examples that such a numerical approach is needed to determine the spatial breakdown point for two or more dimensions, even on a regular square sample grid, since previous conjectures about the spatial breakdown point in two dimensions do not hold. The ‘average spatial breakdown point’ has been used as a basis for practical guidelines on the intensity of contaminating processes that can be tolerated by robust variogram estimators. It is the largest proportion of contaminating observations in a data set such that the breakdown point of the variance estimator used to obtain point estimates of the variogram is not exceeded by the expected proportion of contaminated pairs of observations over any lag. In this paper the behaviour of the average spatial breakdown point is investigated for cases where the contaminating process is spatially dependent. It is shown that in two dimensions the average spatial breakdown point is 0.25. Finally, the ‘empirical spatial breakdown point’, a tool for the exploratory analysis of spatial data thought to contain outliers, is introduced and demonstrated using data on metal content in the soils of Sheffield, England. The empirical spatial breakdown point of a particular data set can be used to indicate whether the distribution of possible contaminants is likely to undermine a robust variogram estimator.     

Stability Analysis and Numerical Simulation of Differential Frost Heave

Differential frost heave is often implicated in the formation of patterned ground in regions subject to recurrent freezing and thawing. A linear stability analysis (LSA) indicates that a continuum model of frost heave is linearly unstable under typical natural freezing conditions of silty-clay soils. A two-dimensional non-linear numerical analysis corroborates the frozen time LSA results, and also indicates the importance of non-linear and time-dependent terms that ultimately lead to a preferred mode, which the LSA fails to predict. Instability of the one-dimensional solution occurs at shallow freezing depths and near-zero surface loads when positive perturbations in the ice content at the freezing front lead to a concomitant increase in thermomolecular pressure and upward ice velocity. Differential frost heave can then occur because of the increased heat flux from the perturbed surfaces. A three-dimensional model using random initial surface perturbations indicates that regular surface patterns will evolve with a length scale in the order of 2–4 meters, which corresponds quite closely with naturally-occurring non-sorted patterned ground.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

焦煤集团九里山矿14101工作面突水治理浅析

通过对焦煤集团九里山矿14101工作面底板突水经过的研究,认为工作面底板突水的水源是是Lg灰岩水,其直接补给水源为L2灰岩水,间接补给水源为O2灰岩水;导水通道为底板采动裂隙带。采取了在工作面打水闸墙,利用钻孔对采空区灌注骨料,增加水流阻力,切断导水通道,封堵水源,加固隔水层的地面注浆和井下工程相结合的综合治水方案。通过观察,工作面中间巷无水流出,突水前后水位变化明显,在较短时间内获得了堵水成功。取得了显著的经济效益和社会效益,为矿井突水治理积累了经验。     

Long-term seismic noise investigations on Shikotan Island: First results

In 2003–2004, long-term seismic noise observations were launched on Shikotan Island (Lesser Kuril Range) based on the “Shikotan” dormant regional seismic station. The geological and geophysical data on the registration area are reported. Information about the equipment and its technical specifications is given. The precursors to the strongest local earthquakesthat occurred in the Shikotan Island region in January 2005–March 2007 are identified.     

Detection of a new interstellar molecular ion, H2COH+ (protonated formaldehyde)

A new interstellar molecular ion, H2COH+ (protonated formaldehyde), has been detected toward Sgr B2, Orion KL, W51, and possibly in NGC 7538 and DR21(OH). Six transitions were detected in Sgr B2(M). The 1(1,0)-1(0,1) transition was detected in all sources listed above. Searches were also made toward the cold, dark clouds TMC-1 and L134N, Orion (3N, 1E), and a red giant, IRC + 10216, without success. The excitation temperatures of H2COH+ are calculated to be 60-110 K, and the column densities are on the order of 10(12)-10(14) cm-2 in Sgr B2, Orion KL, and W51. The fractional abundance of H2COH+ is on the order of 10(-11) to 10-(9), and the ratio of H2COH+ to H2CO is in the range 0.001-0.5 in these objects. The values in Orion KL seem to be consistent with the "early time" values of recent model calculations by Lee, Bettens, & Herbst, but they appear to be higher than the model values in Sgr B2 and W51 even if we take the large uncertainties of column densities of H2CO into account. We suggest production routes starting from CH3OH may play an important role in the formation of H2COH+.     

Octahedral-site Fe2+ quadrupole splitting distributions from Mössbauer spectroscopy along the (OH, F)-annite join

We have performed a detailed Mössbauer study of synthetic annites on the (OH, F)-join. Recently developed data treatment and spectral analysis methods were used to extract true intrinsic Fe²⁺ quadrupole splitting distributions (QSDs) that represent the most information that can be resolved from the spectra. The overall room temperature (RT) QSDs can be consistently interpreted in terms of four QSD contributions (or populations) centered at: QS_HH2.55 mm/s for Fe²⁺O₄(OH)₂ octahedra (cis and trans not resolved), QS_HF 2.35 mm/s for Fe²⁺O₄(OH)F octahedra (cis and trans not resolved), QS_cFF2.15 mm/s for cis-Fe²⁺O₄F₂ octahedra, and QS_tFF 1.5 mm/s for trans-Fe²⁺O₄F₂ octahedra. Each such contribution has a width ( 0.2 mm/s) caused by distortions of the octahedra. Minor contributions due to Fe²⁺O₅(OH) and Fe²⁺O₅F octahedra probably also contribute to the overall Fe²⁺ QSDs. The ferric iron spectral components were also characterized. Here, two distinct types of octahedral Fe³⁺ contributions are seen and interpreted as being due mainly to Fe³⁺O₅OH and Fe³⁺O₅F octahedra, respectively. Tetrahedral Fe³⁺ is seen only in the OH-annite end-member and the total Fe³⁺ content drops significantly on addition of F. On leave from: Department of Materials Physics, University of Science and Technology Beijing, 100083 Beijing, China     

Extreme value analysis of diamond-size distributions

Extreme value analysis provides a semiparametric method for analyzing the extreme long tails of skew distributions which may be observed when handling mining data. The estimation of important tail characteristics, such as the extreme value index, allows for a discrimination between competing distribution models. It measures the thickness of such long tailed distributions, if only a limited sample is available. This paper stresses the practical implementation of extreme value theory, which is used to discriminate a lognormal from a mixed lognormal distribution in a case study of size distributions for alluvial diamonds.