新鲜煤和氧化煤自燃倾向性的FTIR对比分析

利用傅立叶变换红外光谱（FTIR）分析技术,对乌达矿区低温氧化和新鲜煤样进行了红外光谱测试,通过两类煤样的特征官能团吸光度的对比分析．揭示了煤在氧化前后特征官能团的变化规律：氧化后的煤样中脂肪族和芳烃CH减少,酸酐C=O从无到有,芳烃骨架C=C基本不变,进而研究了其对煤自燃倾向性的影响。结果表明脂肪族和羟基含量越高,煤自燃倾向性越大。     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

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.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

地震信号的相位恢复及误差分析

常见的几种地震子波有:零相位子波、最小相位子波、最大相位子波和混合相位子波。以不同类型的信号为激发子波,对不同地层模型合成地震道。仅用原始合成地震道记录的振幅谱来恢复原始记录,讨论不同方法下,由于激发子波的不同,对比分析原始记录与恢复记录的差异,及其恢复效果。     

角度格式转换过程中可靠算法研究

提出了一种可靠算法是将原始度.分秒加上一个与之同符号,且不影响计算精度的有效微小量ｅ,然后再进行格式转换,能有效减少工程返工和损失。