Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

BACKGROUND STUDY ON NON-RESIDUAL PHASE OF ELEMENTS IN CHANGJIANG AND HUANGHE ESTUARINE SEDIMENTS

Analysis of Zn, Cu, Pb, Co, Cr, Li, Ni, K, Al, Fe extracted by 1 mol/L HCl or 0.5 mol/LHCl/H_2O_2, showed concentrations of Zn, Cu, Pb, Co, Cr, Fe, Ni were significantly correlated with Li, Al,K, and clay. Two methods are used to indicate the background value of the non-residual phase of elementsin sediments, and are the same as the methods used to indicate the background value of totalconcentrations in sediments. The first method uses correlograms and regression equations,the second usesthe mean element concentrations normalized with grain size. Li, Al, K can be used as reference elements to determine the background value of Zn, Cu, Pb, Co,Cr, Ni, Fe, while the clay concentration's correlation with some extractable concentrations can be used tocalculate the background value of the non-residual phase of elements as a percentage of clay concentrationin the sediments. Based on this study, the concept of using the background value of the non-residualphase of elements to compare the pollution level in differ     

广西沙田柚农业地质

沙田柚生长好坏、产量品质高低 ,与所在地质环境有密切关系 ,种植在碎屑岩和花岗岩环境区及河流阶地、江心洲、山麓地貌环境区的沙田柚 ,长势好、产量品质高。种植在碳酸盐岩环境区及山坡型、岩溶洼地地貌区的沙田柚 ,长势差、产量品质低。影响沙田柚长势、产量品质的关键元素是Na,Na是沙田柚的特需元素。     

塔中4号油藏渗透率分布的非均质性及渗流方向研究

研究油藏内部油水渗流方向时,主要考虑渗透率影响,引入两个水平方向渗透率比、水平与垂向渗透率比参数来分析。前者反映样品所在储层平面上两个水平方向渗流能力的差异,后者反映样品所在储层水平与垂向渗流能力的差异,并经小柱样测试和全岩心测试两种方法确定两个参数。通过对塔中4号油藏CⅠ油组9口井2306个样品两种渗透率比值分布直方图分析,研究了CⅠ油组渗透率的非均质性及其渗流方向:CI油组储层在水平方向的渗流能力有很大差异,具有明显的方向性。储层垂向渗流能力明显小于水平方向的渗流能力。     

成都市冬季重污染日气象条件成因分析及预报

本文通过对成都10个重污染日进行天气学分析,将污染浓度与气象要素进行聚类、研究了重污染日形成原因及污染浓度与气象要素的关系。在此基础上,建立了SO_2、TSP日平均浓度分级预报方程。     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.