Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite Element Modeling of Swelling Soil Pressures Using Joint/Interface Elements

There could be some discontinuities in a soil media such as layers, earth structures, cracks, and fissures so that estimation of stresses and deformations in these types of soil masses are somewhat different from continuous masses. The discontinuities in a soil mass could be considered as a special link between two blocks. Transmitted swelling pressures affected by the soil properties of the backfill especially at the interface between the backfill and retaining wall. In order to estimate the transmitted swelling pressure distribution behind a retaining wall, using with interface element, a new finite element model and a code (SWELPRES) have been developed. To evaluate the effect of backfill characteristics and interface element from the transmitted lateral swelling pressures, four cases of thickness of backfills with or without interface elements are investigated.     

钱二块铀矿床采铀注液结垢趋势的理论分析与预测

用离子系数矩阵法对钱二块铀矿床地浸采铀试验注液中存在的独立离子反应进行了确定,由溶度积规则及反应平衡原理,借助Matlab编程,估算出了钱二块铀矿床地浸采铀试验注液可能会生成的沉淀物种类及数量,为防治结垢提供了一定的理论依据.     

贵州修文小山坝铝土矿中镓等伴生元素分布规律研究

通过ICP-MS等分析方法对贵州小山坝铝土矿中伴生元素的研究,发现该矿床TiO2、Ga、Ge、Sc、Nb、Li和REE等有用伴生元素在铝土矿矿体中均有较高富集,其含量大多接近或超过了工业综合利用标准,具有较大的经济价值,特别是矿区富铁铝土矿石中Ga和Sc等含量异常高.贵州铝土矿资源丰富,其伴生元素在铝土矿冶炼过程中所产生的赤泥中更为富集,因此,贵州铝土矿综合利用前景十分巨大,这些元素潜在的经济价值,在某种程度上,甚至会超过主导产品氧化铝的价值.开展铝土矿中镓等伴生有用元素的分布及赋存状态的研究,对于镓等稀散元素在特殊地质地球化学环境中超常富集机制、矿产资源综合价值的重新认识评估和矿产资源的综合利用将十分重要,不仅可为铝工业注入新的活力,增加新的经济增长点,而且还可以解决大量赤泥堆放和利用问题,同时获得较大的社会效益与环境效益.     

冲绳海槽有孔虫壳体的微量元素Sr,Nd同位素地球化学

钙质生物壳的微量元素组合和Sr,Nd同位素组成是识别海底混合源沉积物中生物源物质相对贡献的重要参数.冲绳海槽有孔虫壳体强烈富集Sr,P,Mn和Ba,富集Li,U,Th,Sc,Co,Pb,Zn,Cr,Rb,Y,Sb和轻稀土元素,弱富集V,Ga,Zr,Nb,Cd和中稀土元素,相对贫Ge,Mo,In,Sn,Cs,Hf,Ta,W,Tl,Bi和重稀土元素,海水中微量元素的背景含量和生物活动对微量元素的选择性吸收是有孔虫壳体中微量元素发生富集和贫化的主要机制,冲绳海槽有孔虫壳体的稀土元素配分模式与海水和太平洋有孔虫的有明显差异,表现出中稀土元素相对富集,并具有微弱的负Ce异常.有孔虫壳体的Sr,Nd同位素比值也与大洋海水不同,分别为0.709769和0.512162,前者略高于大洋海水,后者略低于大洋海水,表明冲绳海槽海水明显受大陆河水影响.     

Numerical and experimental study on hydrodynamic interaction of side-by-side moored multiple vessels

This paper aims to investigate the basic interaction characteristics of side-by-side moored vessels both numerically and experimentally. A higher-order boundary element method (HOBEM) combined with generalized mode approach is applied to analysis of motion and drift force of side-by-side moored multiple vessels (LNG FPSO, LNGC and shuttle tankers). Model tests were carried out for the same floating bodies investigated in the numerical study in regular and irregular waves. Global and local motion responses and drift forces of three vessels are compared with those of calculations. Discussions is highlighted on applicability of numerical method to prediction of sophisticated multi-body interaction problem of which motion behavior is very important to analysis of mooring dynamics of deep sea floating bodies.     

Imaging of Buried Archaeological Materials: The Reflection Properties of Archaeological Wood

Effective marine archaeological site management demands detailed information on not only the spatial distribution of artefacts but also the degradation state of the materials present. Although sonar methods have frequently been used in an attempt to detect buried wooden shipwrecks they are currently unable to indicate their degradation state. To assess the sensitivity of acoustic measurements to changes in the degradation state of such material, and hence the potential for sonars to quantify degradation, laboratory measurements of compressional wave velocity, as well as bulk density for oak and pine samples, in varying states of decay, were undertaken. These data enabled the calculation of theoretical reflection coefficients for such materials buried in various marine sediments. As wood degrades, the reflection coefficients become more negative, resulting in the hypothesis that the more degraded wood becomes, the easier it should be to detect. Typical reflection coefficients of the order of −0.43 and −0.52 for the most degraded oak and pine samples in sand are predicted. Conversely, for wood exposed to seawater the predicted reflection coefficients are large and positive for undegraded material (0.35 for oak, 0.18 for pine) and decrease to zero or slightly below for the most degraded samples. This indicates that exposed timbers, when heavily degraded, can be acoustically transparent and so undetectable by acoustic methods. Corroboration of these experimental results was provided through comparison with high resolution seismic reflection data that has been acquired over two shipwrecks.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

Drag Coefficient,Dynamic Roughness and Reference Wind Speed

Surface waves are the roughness element of the ocean surface. The parameterization of the drag coefficient of the ocean surface is simplified by referencing to wind speed at an elevation proportional to the characteristic wavelength. The dynamic roughness is analytically related to the drag coefficient. Under the assumption of fetch limited wave growth condition, various empirical functions of the dynamic roughness can be converted to equivalent expressions for comparison. For datasets covering a wide range of the dimensionless frequency (inverse wave age), it is important to account for the variable rate of wave development at different wave ages. As a result, the dependence of the Charnock parameter on wave age is nonmonotonic. Finally, the analysis presented here suggests that the significant wave steepness is a sensitive property of the ocean surface and a single variable normalization of the dynamic roughness using a wavelength or wave height parameter actually produces more robust functions than bi-variable normalizations using wave height and wave slope.