An analysis of coal bump liability in a bump prone mine

Summary This paper presents an analysis of causes of bump in a deep coal mine. The analysis consists of in-mine measurement determination, analysis of physical and mechanical properties of mine materials and simulation of mine layout and mining sequence. Major factors contributing to the causes of bump have been reviewed and those that are related to this case have been identified and simulated analytically using a three-dimensional finite element method.     

Complex behaviour of a coal sample in a large scalein situ compression test

Summary A large scalein situ compression test of a laminated coal cube is described. The cube was uniformly loaded by a system of hydraulic jacks supplied simultaneously from one pump. At the same time both load and deformations were controlled. Deformations were measured by 60 longitudinal gauges located at the side surface of the cube in the regular net. Due to the heterogeneity of the rock the readings showed scatter. As a first approximation, the scattered data were fitted to the plane using the least square method. Orientations of the planes were analysed and show permanent relative movement between adjacent lamina indicating complex behaviour of the cube under loading. Detailed analysis of this behaviour indicated that the cube was subjected to both bending and twisting.     

桩基低应变完整性检测的分析研究

本文简述了现场浇注混凝土桩的低应变完整性试验的基本理论和方法。文中设定了几个简单的分析模型，利用三维显式有限元法，对有局部材料缺陷桩完成了冲击荷载作用下的动力反应分析。同时，根据分析结果，对低应变完整性检测试验的实施方案，各种因素的影响效果以及试验结果的分析方法等若干问题进行了初步探讨。     

Suction-controlled triaxial testing: Laboratory procedures in relation to resistance envelope methods

In certain tropical slopes it appears possible that soil suction may play a significant role in maintaining stability. This communication outlines a laboratory suction-controlled triaxial test that can be undertaken to validate the threshold suction predicted by resistance envelope methods. It is shown that for two sites examined in St. Lucia, West Indies, the results from such a test confirm results obtained by use of resistance envelopes.     

Hypotheses for earthquake occurrences

Very little work has been done in generating alternatives to the Poisson process model. The work reported here deals with alternatives to the Poisson process model for the earthquakes and checks them using empirical data and the statistical hypothesis testing apparatus. The strategy used here for generating hypotheses is to compound the Poisson process. The parameter of the Poisson process is replaced by a random variable having prescribed density function. The density functions used are gamma, chi and extended (gamma/chi). The original distribution is then averaged out with respect to these density functions. For the compound Poisson processes the waiting time distributions for the future events are derived. As the parameters for the various statistical models for earthquake occurrences are not known, the problem is basically of composite hypothesis testing. One way of designing a test is to estimate these parameters and use them as true values. Momentmatching is used here to estimate the parameters. The results of hypothesis testing using data from Hindukush and North East India are presented.     

Modern Numerical Ephemerides of Planets and the Importance of Ranging Observations for Their Creation

The JPL planetary and lunar ephemerides – DE200/LE200, DE403/LE403, DE405/LE405 and the planetary and lunar ephemerides, EPM87, EPM98, and EPM2000, constructed in the Institute of Applied Astronomy of RAS are described. Common properties and differences of the various ephemerides are given. Graphical comparisons of the DE ephemerides with each other and with the EPM ephemerides are presented. A fairly good agreement of planetary orbits is between DE403, DE405 and EPM98, EPM2000, respectively, over the interval of 120 years (1886–2006) covered by EPM98 and EPM2000. Some differences are explained by a slight disagreement in representing the orbits of Ceres, Pallas, and Vesta as they affect the planets. The accurate radar observations of planets and spacecraft make it possible not only to improve the orbital elements of planets but to determine a broad set of astronomical constants as well: km/AU, parameters of Mars rotation including its precessional rate, the masses of Jupiter, Ceres, Pallas, and Vesta, relativistic parameters of the PPN formalism, the variability of the gravitational constant G. These have been obtained in the fitting process of the DE405 and EPM2000 ephemerides to observational data, including nearly 80000 American and Russian radar observations of planets (1961–1997), ranging and doppler to the Viking and Pathfinder landers, and other miscellaneous measurements from various sources and spacecraft.     

Assessing spatial and attribute errors in large national datasets for population distribution models: a case study of Philadelphia county schools

Geospatial technologies and digital data have developed and disseminated rapidly in conjunction with increasing computing efficiency and Internet availability. The ability to store and transmit large datasets has encouraged the development of national infrastructure datasets in geospatial formats. National datasets are used by numerous agencies for analysis and modeling purposes because these datasets are standardized and considered to be of acceptable accuracy for national scale applications. At Oak Ridge National Laboratory a population model has been developed that incorporates national schools data as one of the model inputs. This paper evaluates spatial and attribute inaccuracies present within two national school datasets, Tele Atlas North America and National Center of Education Statistics (NCES). Schools are an important component of the population model, because they are spatially dense clusters of vulnerable populations. It is therefore essential to validate the quality of school input data. Schools were also chosen since a validated schools dataset was produced in geospatial format for Philadelphia County; thereby enabling a comparison between a local dataset and the national datasets. Analyses found the national datasets are not standardized and incomplete, containing 76 to 90 percent of existing schools. The temporal accuracy of updating annual enrollment values resulted in 89 percent inaccuracy for 2003. Spatial rectification was required for 87 percent of NCES points, of which 58 percent of the errors were attributed to the geocoding process. Lastly, it was found that by combining the two national datasets, the resultant dataset provided a more useful and accurate solution.     

A case study to verify methods for estimating transmissivity distributions along boreholes

Two methods, both using data from fixed-interval-length hydraulic tests and geological mapping on different scales, have been utilised to estimate transmissivity distributions of individual fractures. Individual fracture properties are of importance in estimating transport of water and solutes; they influence the spreading of grout when sealing tunnels. One of the methods is non-parametric, while the other is based on combinatorics and the multiplication principle. The study verified the usefulness of both methods for estimating the probability of conductive and non-conductive fractures, and the transmissivity distributions. Data for estimating transmissivity distributions originate from flow logging of a borehole in Sweden. Fracture frequency data are from core logging and images using a borehole image processing system (BIPS). Estimated transmissivity distributions, based on flow log data from 3-m sections (test scale length, L: 3 m, and step length, dL: 3 m), were compared with transmissivities obtained directly from detailed, overlapping flow logging which identify individual conductive fractures (test scale, L: 0.5 m, step length, dL: 0.1 m). The latter data are considered to be an acceptable approximation of the individual fracture transmissivities; the agreement, compared with the calculated transmissivity distributions based on 3-m section data, is good for both methods.

---

Resumen Se han utilizado dos métodos, ambos procesando datos de pruebas hidráulicas de longitud e intervalo fijo y mapas geológicos en distintas escalas, para estimar las distribuciones de transmisividad de fracturas individuales. Las propiedades de fracturas individuales son de importancia para estimar el transporte de agua y solutos; las propiedades influyen en la propagación de lechada al sellar túneles. Uno de los métodos es no paramétrico mientras que el otro se basa en el principio de multiplicación y combinatorio. El estudio verifica la utilidad de ambos métodos para estimar la probabilidad de fracturas conductivas y no conductivas y las distribuciones de transmisividad. Los datos para estimar las distribuciones de transmisividad se originan del registro de flujo de un pozo en Suecia. Los datos de frecuencia de fracturas provienen de registros de núcleo e imágenes usando un Sistema de Procesamiento de Imágenes de Pozo (BIPS). Se compararon las distribuciones estimadas de transmisividad basadas en datos del registro de flujo en secciones de 3 m (longitud de escala de la prueba L: 3 m, y escalón de longitud dL: 3 m), con transmisividades obtenidas directamente de registros detallados de flujo superpuesto el cual identifica fracturas conductivas individuales (escala de la prueba, L: 0.5 m, escalón de longitud, dL: 0.1 m). Los últimos datos se consideran una aproximación aceptable de las transmisividades de fracturas individuales; la consistencia, comparada con las distribuciones de transmisividad calculada basadas en datos de sección de 3 m, es buena para ambos métodos.

---

Résumé Deux méthodes ont été testées afin d’estimer la distribution des transmissivités de fractures unitaires; toutes deux utilisent des données issues de tests hydrauliques effectués à intervalles réguliers et de cartographies géologiques à plusieurs échelles. La connaissance des propriétés des fractures unitaires est importante pour estimer le transport d’eau et de solutés; elles influencent en outre l’étalement du coulis lors du scellement des tunnels. La première méthode est non-paramétrique, et la seconde repose sur l’analyse combinatoire et le principe de multiplication. Cette étude a vérifié l’utilité des deux méthodes pour estimer la probabilité d’occurrence de fractures conductrices et non-conductrices, et la distribution des transmissivités. Les distributions des transmissivités ont été estimées à partir de diagraphies effectuées sur un forage situé en Suède. Les densités de fracturation proviennent de carottages et d’imagerie BIPS (Borehole Image Processing System). Les distributions estimées des transmissivités, basées sur des diagraphies de flux effectuées sur des sections de 3 m (échelle du test L:3 m; pas dL:3 m), ont été comparées aux transmissivités obtenues sur des diagraphies de flux détaillées et se chevauchant, identifiant les fractures conductrices unitaires (échelle du test L: 0.5 m ; pas dL: 0.1 m). Il est admis que ces dernières données représentent une approximation satisfaisante des transmissivités des fractures unitaires; les deux méthodes sont en accord avec les distributions des transmissivités calculées à partir des données acquises sur des sections de 3 m.

     

Hydraulic properties of the crystalline basement

Hydraulic tests in boreholes, up to 4.5 km deep, drilled into continental crystalline basement revealed hydraulic conductivity (K) values that range over nine log-units from 10⁻¹³−10⁻⁴ m s⁻¹. However, K values for fractured basement to about 1 km depth are typically restricted to the range from 10⁻⁸ to 10⁻⁶ m s⁻¹. New data from an extended injection test at the KTB research site (part of the Continental Deep Drilling Program in Germany) at 4 km depth provide K=5 10⁻⁸ m s⁻¹. The summarized K-data show a very strong dependence on lithology and on the local deformation history of a particular area. In highly fractured regions, granite tends to be more pervious than gneiss. The fracture porosity is generally saturated with Na–Cl or Ca–Na–Cl type waters with salinities ranging from <1 to >100 g L⁻¹. The basement permeability is well within the conditions for advective fluid and heat transport. Consequently, fluid pressure is hydrostatic and a Darcy flow mechanism is possible to a great depth. Topography-related hydraulic gradients in moderately conductive basement may result in characteristic advective flow rates of up to 100 L a⁻¹ m⁻² and lead to significant advective heat and solute transfer in the upper brittle crust. An erratum to this article can be found at     

Parameter identification of framed structures using an improved finite element model‐updating method—Part I: formulation and verification

In this study, we formulate an improved finite element model‐updating method to address the numerical difficulties associated with ill conditioning and rank deficiency. These complications are frequently encountered model‐updating problems, and occur when the identification of a larger number of physical parameters is attempted than that warranted by the information content of the experimental data. Based on the standard bounded variables least‐squares (BVLS) method, which incorporates the usual upper/lower‐bound constraints, the proposed method (henceforth referred to as BVLSrc) is equipped with novel sensitivity‐based relative constraints. The relative constraints are automatically constructed using the correlation coefficients between the sensitivity vectors of updating parameters. The veracity and effectiveness of BVLSrc is investigated through the simulated, yet realistic, forced‐vibration testing of a simple framed structure using its frequency response function as input data. By comparing the results of BVLSrc with those obtained via (the competing) pure BVLS and regularization methods, we show that BVLSrc and regularization methods yield approximate solutions with similar and sufficiently high accuracy, while pure BVLS method yields physically inadmissible solutions. We further demonstrate that BVLSrc is computationally more efficient, because, unlike regularization methods, it does not require the laborious a priori calculations to determine an optimal penalty parameter, and its results are far less sensitive to the initial estimates of the updating parameters. Copyright © 2006 John Wiley & Sons, Ltd.