Distribution and environmental impact of coal-mining wastes in Upper Silesia, Poland

 About 50 million tonnes/year of waste rock from coal-mining is generated in the limited area of the thickly populated Upper Silesian Coal Basin (USCB) in Poland. There are 380 coal-mining waste dumps, including 76 active dump sites covering over 2,000 ha. About 15-16 million tonnes/year of waste rock is being reused for civil engineering purposes in the same area. This brings about a problem of ground water deterioration by constituents leached from waste rock exposed to atmospheric conditions. The major factors determining the ground water contamination potential from waste rock are chloride salinity, sulfur content and acid generation potential. The concept behind the presented studies was to provide data for correct evaluation and prediction of contaminant release from the waste rock, based on the characterization of coal-mining waste properties, as well as on long-term laboratory, lysimetric and field studies. The results show that coal-mining waste dumps can be a long-term source of ground water contamination, lasting for decades and increasing with time. Ground water down-gradient from the disused 15–30-years-old part of the studied dump displays high and increasing acidification, high TDS, SO₄, and the highest, still increasing concentrations of Mn, Fe and Zn. Cost-effective and efficient pollution control measures, similar to the presented design and construction elements of the dump site, can mitigate the negative environmental impacts. Received: 3 July 1997 · Accepted: 9 September 1998     

CALCULATING THE DYNAMIC STIFFNESS MATRIX OF 2-D FOUNDATIONS BY DISCRETE WAVE NUMBER INDIRECT BOUNDARY ELEMENT METHODS

Apart from some special cases, calculating the dynamic stiffness matrix of foundations on a layered half-space, especially for embedded foundations, is computationally expensive. An efficient method for two-dimensional foundations in a horizontally layered soil media is presented in this paper. This method is based on indirect boundary element methods and uses discrete wave number solution methods for calculating Green's functions for displacements and analytical methods for the integrations over the boundary. For surface foundations, the present method applies at all frequencies. For embedded foundations or for constructing energy transmitting boundaries, because the free-field part is modelled by boundary elements and the excavated part is modelled by finite elements, the present method applies only at low frequencies for the spring coefficients (the real parts of the dynamic stiffness matrix) but applies at all frequencies for the damping coefficients (the imaginary part of the dynamic stiffness matrix) for undamped sites. The novelty of the method can be used for three-dimensional foundations. © 1997 by John Wiley & Sons, Ltd.     

Third-Order Stokes Wave Solutions of the Free Surface Capillary-Gravity Wave and the Interfacial Internal Wave

Based on the Stokes wave theory, the capillary-gravity wave and the interfacial internal wave in two-layer constant depth''s fluid system are investigated. The fluids are assumed to be incompressible, inviscid and irrotational. The third-order Stokes wave solutions are given by using a perturbation method. The results indicate that the third-order solutions depend on the surface tension, the density and the depth of each layer. As expected, the first-order solutions are the linear theoretical results (the small amplitude wave theoretical results). The second-order and the third-order solutions describe the nonlinear modification and the nonlinear interactions. The nonlinear impact appears not only in the n (n≥2) times'' high frequency components, but also in the low frequency components. It is also noted that the wave velocity depends on the wave number, depth, wave amplitude and surface tension.     

2010年6月广西凌云-凤山交界3级震群应力水平分析

2010年6月28日至8月15日,在广西凌云与凤山交界持续发生了高密集、高集中、高频率的震群活动,共记录地震3 071次,其中ML2~2.9级地震41次,3级以上地震3次,最大为7月1日10时27分3.2级地震。对此次震群的空间分布范围、2级以上地震震源机制解P轴优势方向、b值特征、应力降等进行了分析研究。结果表明,此次震群震中集中,深度较浅,震源机制解P轴没有形成优势分布,序列参数b值较高,约为1.21,应力降主要分布在0.19~2.0 Mpa之间,认为广西凌云-凤山交界3级震群是属于低应力背景下的地震活动,不具有前兆震群和前震序列的性质。分析结果对正确了解此次震群的特征、正确判断地震趋势具有一定的帮助作用。     

非线性四阶四点边值问题的正解

利用锥不动点定理,并结合Green函数性质,证明一类四阶四点边值问题的正解的存在性.     

Numerical modelling of rise and fall of a dense layer in salt diapirs

Numerical models are used to study the entrainment of a dense anhydrite layer by a diapir. The anhydrite layer is initially horizontally embedded within a viscous salt layer. The diapir is down-built by aggradation of non-Newtonian sediments ( n = 4, constant temperature) placed on the top of the salt layer. Several parameters (sedimentation rate, salt viscosity, perturbation width and stratigraphic position of the anhydrite layer) are studied systematically to understand their role in governing the entrainment of the anhydrite layer. High sedimentation rates during the early stages of the diapir evolution bury the initial perturbation and, thus, no diapir forms. The anhydrite layer sinks within the buried salt layer. For the same sedimentation rate, increasing viscosity of the salt layer decreases the rise rate of the diapir and reduces the amount (volume) of the anhydrite layer transported into the diapir. Model results show that viscous salt is capable of carrying separate blocks of the anhydrite layer to relatively higher stratigraphic levels. Varying the width of the initial perturbation (in our calculations 400–800 m), from which a diapir triggers, shows that wider diapirs can more easily entrain an embedded anhydrite layer than the narrower diapirs. The anhydrite layer is entrained as long as rise rate of the diapir exceeds the descent rate of the denser anhydrite layer. We conclude that the four parameters mentioned above govern the ability of a salt diapir to entrain an embedded dense layer. However, the model results show that the entrained blocks inevitably sink back if the rise rate of the diapir is less than the rate of descent of the anhydrite layer or the diapir is permanently covered by a stiff overburden in case of high sedimentation rates.     

Interpolation using a generalized Green's function for a spherical surface spline in tension

A variety of methods exist for interpolating Cartesian or spherical surface data onto an equidistant lattice in a procedure known as gridding. Methods based on Green's functions are particularly simple to implement. In such methods, the Green's function for the gridding operator is determined and the resulting gridding solution is composed of the superposition of contributions from each data constraint, weighted by the Green's function evaluated for all output–input point separations. The Green's function method allows for considerable flexibility, such as complete freedom in specifying where the solution will be evaluated (it does not have to be on a lattice) and the ability to include both surface heights and surface gradients as data constraints. Green's function solutions for Cartesian data in 1-, 2- and 3-D spaces are well known, as is the dilogarithm solution for minimum curvature spline on a spherical surface. Here, the spherical surface case is extended to include tension and the new generalized Green's function is derived. It is shown that the new function reduces to the dilogarithm solution in the limit of zero tension. Properties of the new function are examined and the new gridding method is implemented in Matlab® and demonstrated on three geophysical data sets.