Study on numerical simulation to water lowering in an open coal mine digging

In view of the situation of excavation of open coal mine for the underground water disaster,we should carry out simulation studies for the numerical value of the water lowering project and improve the accuracy and the level of the water lowering project.On the basis of the hydrological geological conditions of certain open mine digging,a more reasonable seepage numerical model was built according to MODFLOW.It was simulated in advance that the process of the confined water level descending with the time,and combining with the actual observations to test the correctness of the model.The calculation showed that the results coincided well with the results of actual measurement.Based on this,different water lowering numerical simulations were built for the open coal mine digging.It could be simulated and forecast that the changes of the groundwater level in drainage process within and outside the mine pit,and it was quantitatively assessed that the possible water lowering result of the opencast water drainage process,which provide an important basis for the actual water lowering project and the possible project disposal.     

On salinity regimes and the vertical structure of residual flows in narrow tidal estuaries

An examination is made of the circulation in narrow estuaries subject to a predominant tidal forcing. Velocity structures are derived separately for residual flow components associated with (a) river flow, (b) wind stress, (c) a well-mixed longitudinal density gradient and (d) a fully stratified saline wedge. Dimensionless parameters are introduced to indicate the magnitude of each component and these parameters are evaluated for 9 major estuaries, thereby revealing their sensitivity to each component.For a channel of constant breadth and depth, formulae are deduced for the length of saline intrusion, L. Comparisons with observed data show that such formulae may be used with confidence to predict changes in L arising from variations in river flow, tidal range or channel depths.The level of stratification is shown to be related to a product of two parameters, one associated with velocity structure and a second involving the square of the ‘flow ratio’ $\text{u}\text{u}\text{?}$ (i.e. residual velocity/amplitude of the tidal velocity). This relationship provides a simple classification system for estuarine stratification which can be used to indicate the sensitivity of any particular estuary to changing conditions.     

广西芒场锡多金属矿田稳定同位素组成对矿床成因探讨

笔者曾对广西芒场锡多金属矿田的稳定同位素进行研究。本文根据硫、铅、氢、氧、碳等稳定同位素组成和锶的初始值提供的信息.探讨了矿床成因。并结合矿田矿床地质特征、控矿条件及有关统计参数,参考前人对矿床认识的基础上,修正提出了该矿田混合热液成矿模式,可供类似矿床研究和找矿的参考。     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of double cold cores of the 35°N section in the Yellow Sea with a wave-tide-circulation coupled model

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter. This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

Recharge and salinization in the Madrid Basin

The purpose of this work was to reinvestigate the existing hydrogeological conceptual model of the basin of Madrid, Spain. A cumulative chemical isotopic diagram which enabled the distinction between different groups of water as well as calculation of the mode of their blending was used for this investigation. It was found that the groups of discharge were lighter in their isotopic composition than that of recharge. The previous explanation of this fact, backed by carbon-14 dating, was the long residence time due to flow lines going down to depths of more than 1000 m. This flow model assumes homogenous conditions to these depths. This assumption can not be supported by evidence from deep wells. Thus a modified model is suggested which maintains homogenous conditions only to about 300 m and a deep confined aquifer below containing paleowater. The higher degree of depletion of this water has been explained by a colder climate on top of an altitude effect. Another interesting observation was the correlation between the isotopic composition of the rains, the month of the rain event and the composition of the recharge group groundwater. It could be seen that the winter rains resemble the groundwater composition, which shows that practically all the spring and summer rains were evapotranspirated.     

A laboratory test of the soil chemical submodels of two models of catchment acidification

Parameters for ion exchange selectivity and aluminium hydroxide dissolution in the soil chemical submodels used in applications of the Birkenes model and of MAGIC are compared and several discrepancies identified for organic soils. A laboratory column simulation of the soil chemical submodels is proposed and applied to soils from the Loch Dee area in Galloway. Experimental results were well predicted by a simplified version of MAGIC, with ion exchange selectivity parameters similar to those used in a previous simulation of one subcatchment of Loch Dee. The aluminium hydroxide dissolution parameter used previously was found to be too low for the organic soil materials, where a value of 106 predicted the experimental results more closely. The model developed also included a simple silicate weathering reaction to release base cations into the system. It is concluded that such simple laboratory simulations are useful for independent calibration of the soil chemical submodel of catchment models.     

Tidal Breakup of Comets

In the context of the survival of periodic comets of different origins, rotational breakup and tidal disruption could be important, especially of the short period comets injected from the Kuiper belt. This is because long-period comets from the distant Oort cloud tend to be subject to thermal stress and volatile 'explosion' far more severely. A simple calculation using the Öpik method of random planetary close encounters was performed to estimate the probability of tidal disruption of comets and scattered Kuiper belt objects (SKBOs) during their orbital migration. It was found that a large fraction of the short period comets and SKBOs might have been internally fragmented by single or multiple close encounters with the outer planets.     

Cosmic structure growth and dark energy

Dark energy has a dramatic effect on the dynamics of the Universe, causing the recently discovered acceleration of the expansion. The dynamics are also central to the behaviour of the growth of large-scale structure, offering the possibility that observations of structure formation provide a sensitive probe of the cosmology and dark energy characteristics. In particular, dark energy with a time-varying equation of state can have an influence on structure formation stretching back well into the matter-dominated epoch. We analyse this impact, first calculating the linear perturbation results, including those for weak gravitational lensing. These dynamical models possess definite observable differences from constant equation of state models. Then we present a large-scale numerical simulation of structure formation, including the largest volume to date involving a time-varying equation of state. We find the halo mass function is well described by the Jenkins et al. mass function formula. We also show how to interpret modifications of the Friedmann equation in terms of a time-variable equation of state. The results presented here provide steps toward realistic computation of the effect of dark energy in cosmological probes involving large-scale structure, such as cluster counts, the Sunyaev–Zel'dovich effect or weak gravitational lensing.