含煤层地质环境下地震波场的数值模拟

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。     

虚拟样机技术在地质工程机械领域的应用研究

基于对国内地质工程机械的研发现状的考虑，根据虚拟样机技术的特点，讨论了在地质工程机械领域实施虚拟样机技术的必要件，并提出了建立地质工程机械虚拟样机的策略。联合PRO／E、ANSYS及ADAMS建立了气动潜孔锤的虚拟样机模型，对气动潜孔锤虚拟样机进行了仿真，阐明了虚拟样机技术在地质工程机械领域的应用前景。     

Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite

The physical properties of the hydrous phyllosilicate lizardite have been investigated by atomistic simulation using the GULP code based on transferable semi-empirical interatomic potentials. Lizardite behavior was first investigated during structure relaxation at room temperature. The Helmholtz free energy is minimum for an equilibrium structure that is in agreement with experiment. The bulk, shear, and Young modulii for lizardite were calculated along with the Poisson ratio. From the shear and bulk modulii, we also calculated translational and longitudinal acoustic wave velocities that are important quantities for tectonophysics models. As expected, lizardite is stiffer in the a direction parallel to the layers than in the c perpendicular direction; the variation of the unit cell parameters with pressure is in good agreement with experiment. The cohesive energy between two successive layers along c direction was calculated at 0.33 eV (i.e., 0.11 eV per OH bond) in good agreement with recent ab initio calculations. Upon pressure and temperature variations, we evidenced that structural changes are mainly pressure induced; pressure being accommodated by a decrease of the c parameter up to 10 GPa. We also found that the change of slope in the derivative of the c cell parameter with respect to pressure occurring around 2 GPa originates from the bending of the interlayer hydroxyl groups with respect to the layer normal direction.     

Atomistic simulation of the structure and elastic properties of pyrite (FeS2) as a function of pressure

Interatomic potential parameters have been derived at simulated temperatures of 0 K and 300 K to model pyrite FeS₂. The predicted pyrite structures are within 1% of those determined experimentally, while the calculated bulk modulus is within 7%. The model is also able to simulate the properties of marcasite, even though no data for this phase were included in the fitting procedure. There is almost no difference in results obtained for pyrite using the two potential sets; however, when used to model FeS₂ marcasite, the potential fitted at 0 K performs better. The potentials have also been used to study the high-pressure behaviour of pyrite up to 44 GPa. The calculated equation of state gives good agreement with experiment and shows that the Fe–S bonds shorten more rapidly that the S–S dimer bonds. The behaviour of marcasite at high pressure is found to be similar to that of pyrite.     

Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

Simulation of oasis breeze circulation in the arid region of the Northwestern China

In this paper, the process of oasis-desert circulation (ODC) is simulated by MM5V3.5 model through designing an ideal oasis-desert scheme and assuming that initial atmosphere is at rest (V = 0). The findings showed that the key of forming special oasis boundary structure is the difference of energy and water between oasis and desert. The evaporation of oasis surface consumes heat energy, and the low temperature of oasis causes an oasis breeze circulation (OBC), which drives an ODC with a downdraft over the oasis and an updraft over the desert. Later, the cold, dry and stable boundary over oasis is gradually formed, on the contrary,the atmospheric boundary over desert on the edge of oasis is hot, humid and unstable and its height is about 600 hPa. The updraft over the desert forms a wet ring that acts as a vertical wall weakening the low-level moisture exchange between the oasis and desert. The downdraft of OBC increases the atmospheric stability that reduces the oasis evaporation. The low-level outflow from the oasis (into the desert) prevents the dry, hot air flowing from the desert into the oasis.Thus an oasis self-preservation mechanism may be formed due to OBC. The horizontal area influenced by oasis is twice as oasis area and the vertical range is four times as oasis. The ODC is strong in the daytime and reaches the strongest at 17:00, and the influenced area is the largest at 20:00.     

遥感影像三维可视化研究

利用DEM高程资料实现遥感影像的三维显示可大大提高遥感监测产品的使用价值，而商用软件成本高操作繁琐普及困难。利用计算机图形学知识和Delphi编程语言，开发了一套遥感影像三维可视化软件。对在开发过程中遇到的有关DEM数据的裁剪与拼接、图像配准、地表模拟、透视变换、消隐处理、光照和阴影效果以及颜色变换等问题的解决方案和实现方法作了简要的介绍。     

Indicator Simulation Accounting for Multiple-Point Statistics

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.