Statistical analysis of sand grain/bed collision process recorded by high-speed digital camera

A high‐speed digital camera was employed to record the sand grain/bed collision process. With image processing and a statistical method, a series of parameters of the collision process were obtained. The results show that the collision process of a grain with rebounding can be represented by two parameters: the kinetic energy restitution coefficient and the collision angle. Both parameters satisfy a normal distribution, and they are dependent on one another. With an increase of the collision angle, the distribution of the kinetic energy restitution gradually reduces from a broad to a narrow range with low values. The percentage of vertical velocity restitution coefficients greater than 1 can reach 70% or more, which ensures that the settling time of the sand grains in the air increases and that they receive more energy from the air to progress the saltation movement.     

Techniques to measure grain-size distributions of loamy sediments: a comparative study of ten instruments for wet analysis

Ten instruments for measuring the grain‐size distribution of loamy sediments in water were tested in a comparative study. Instruments tested were: the Malvern Mastersizer S, the Coulter LS 200, the Fritsch Analysette 22 (version C), the Horiba Partica LA‐950, the Sedigraph 5100, the Atterberg cylinder, the Coulter Multisizer 3, the CIS‐100, the EyeTech and the image analysis software Histolab. The first four instruments use laser diffraction. The Sedigraph and Atterberg cylinder are based on sedimentation. The Coulter Multisizer 3 uses impedance measurements and the last three instruments use optical techniques. Grain‐size analyses were carried out on four sediments with median grain diameters of approximately 35, 30, 12 and 9 μm, respectively, and with no particles >90 μm. The laser diffraction instruments produced the best results for the various criteria considered in this study, followed by the Coulter Multisizer 3 (electrical sensing zone technique), the sedimentation techniques and the optical methods. However, the study also showed that there is no such thing as an ‘ideal’ method because techniques may score differently according to the criterion tested. An appropriate evaluation thus requires that various criteria be considered. Issues to be taken into account are: type of sediment, quantity of sediment available, speed of measurement, complexity of the measurement protocol, processing of the data, reproducibility of the results and the specific aims of the study undertaken.     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

EPR investigation of the methyl radical,the hydrogen atom and carbon oxide radicals in Maxixe-type beryl

The EPR spectra of Maxixe-type beryl contain a large number of overlapping signals. The angular dependence of the 1:3:3:1 signal typical for the CH₃ radical shows that this radical is located at the center of the channel cavity with its symmetry axis parallel to the crystal c-axis and is rotating around this axis. Its EPR spectrum is axially symmetric with g _// = 2.00263, g _⊥ = 2.00249 and A_// = 2.288 mT, A_⊥ = 2.256 mT. These anisotropies have the opposite signs of those found for surface-adsorbed methyl radicals. Hydrogen atoms are located at position 2a at the center of the beryl cavity and the EPR parameters of the narrow doublet signal are A ₀ = 1,407 MHz and g = 2.00230. Another doublet signal, which is broader and has axial symmetry with g _// = 2.00265, g _⊥ = 2.00625 and A_// = 0.895 mT, A_⊥ = 0.885 mT, could come from a HCO₃ radical. One narrow and easily saturated signal with g _// = 2.00227 and g _⊥ = 2.00386 is interpreted to arise from a carbon monoxide radical in the beryl channel, oriented with its axis parallel to the crystal c-axis. Additional weak doublet lines, which have similar g values as the carbon monoxide radical, are created by nearby hydrogens. A powder spectrum with g _// = 2.0017 and g _⊥ = 2.0004 appears upon UV irradiation of the single crystal and is easily saturated. This spectrum is interpreted to arise from a carbon dioxide radical, which rotates around its symmetry axis.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

滑坡体三维地质建模与可视化分析

针对滑坡地质研究的自身特点,提出了面向滑坡地质体三维建模的NURBS-BRep混合数据结构和地质结构单元实体构造技术.通过对滑坡基础地质数据的预处理、滑坡数字地形和滑动面的三维建模、地质结构面的系统构造、地质结构体的生成和显示,形成了一套完整的滑坡三维地质模型的构建方法.将该方法应用于某水库滑坡,建立了相应的三维滑坡地质模型,并基于此模型进行了三维剖切分析、数字钻孔、等值线生成、滑块自动剖分、滑坡失稳可视化动态模拟和滑坡体方量精确计算等一系列实用的可视化分析,为滑坡稳定性的准确计算和客观评价提供了有力的支持.     

一种新的波浪变形三维数值模式──0-1混合型边界元

为提高边界元法模拟三维波场波浪变形的数值计算精度,借鉴常数元和线性元剖分方式、波势函数及波势函数法向导数对单元节点设定的各自适应性,提出了一种新的单元剖分模式──0-1混合型边界元,以控制和减缓由于计算误差累计而造成的波浪数值计算上的“横向振动”,借此结合边界元法的分区模式可实现较大范围的波场线性波浪变形计算,并为时域内的波浪非线性变形计算提供时间步长的数值保证.     

渤海三维海洋温度和海流数值预报

"十·五"期间,我国开展了三维海洋温度和海流数值预报的业务化研究工作.经过3年的努力,渤海三维海洋温度和海流数值预报系统研制完成,并于2003年10月,开始试预报.本文对该预报系统以及运行情况进行了介绍,并分析了所存在的问题和发展方向.