多孔介质中非水相流体运移的数值模拟

针对多孔介质中水、气和非水相流体(NAPLs)的多相流动特点,建立了非水相流体(NAPLs)污染物迁移模型,分析了非水相流体在土壤非饱和区和地下水系统中的运移规律。通过有限元数值解对轻非水相流体和重非水相流体在土壤系统中的迁移过程进行模拟,得到了污染物的时空分布特征和污染范围。计算结果表明,数值模拟方法能够合理地描述非水相流体的运移过程和污染特征。土体渗透性和污染物残余饱和度是其重要影响因素。     

土壤中可挥发性污染物清除的离心试验研究

土工离心模拟试验技术是研究环境岩土工程问题的有效手段。本文研究了非水相流体污染物在非饱和土中的迁移以及随后的抽气清除过程。 当离心机运行到要求的加速度时,汽油污染物从地下油罐中释放并在非饱和土中迁移一年,之后采用土壤通气法对污染土壤进行修复。对土壤取样分析,得到污染物在土体中的迁移规律和分布特征。试验结果表明,土壤通气法可以清除非饱和土体中的挥发性有机污染物,是一种有效的原位土壤修复技术。     

基于竖管法的轻非水相液体毛细上升特性研究

被称为“工业血液”的轻非水相液体(LNAPL)及其衍生物在开采、生产、运输过程中所产生的污染已经成为常见的污染物,目前诸多学者对LNAPL污染进行了广泛研究,但对于LNAPL在土壤中的迁移及毛细作用研究尚不充分。本次试验的目的是通过室内模拟试验分析不同竖管直径条件下LNAPL在不同介质中的毛细上升规律,为研究LNAPL对地下水污染提供一定的理论依据。结果表明：影响毛细上升高度的因素大小依次为：溶液>介质>竖管直径;竖管直径与最大毛细上升高度并不是完全成比例关系,并且对毛细上升高度的影响相对较小;水与柴油在不同介质中毛细上升高度、毛细上升速率变化趋势基本一致,但是具体数值上存在差异,柴油的最大毛细上升高度与水相比降低了40%～50%,柴油的最大毛细上升速率与水相比降低30%～50%。这些特征都能够较好地体现LNAPL在不同介质中的毛细上升规律,在认识LNAPL对地下水的污染以及污染土地修复方面具有重要意义。     

非饱和带中非均质条件下LNAPL运移与分布特性实验研究

本文通过室内实验建立物理模型，研究非饱和带中层状非均质条件下轻非水相液体LNAPL（Light Non-Aqueous Phase Liquid）运移机制和分布特性。实验结果发现，介质的结构突变面改变了LNAPL的迁移模式。在LNAPL到达干湿界面以后，由于毛细水压力的作用，其入渗速度明显加快。当LNAPL迁移到夹层上界面时，纵向迁移受到阻碍。随着入渗过程的继续，污染物沿着界面产生向右横向迁移；在夹层上界面下部出现“优势流”现象，污染物由集中走向分散，随着分散块的体积不断增大，最终污染物又连接到一起，由分散走向集中；在夹层下界面附近锋面线集中在一起，说明污染物向下的纵向迁移基本停止。另外，介质的局部密实“透镜体”对LNAPL的局部迁移起到“阻滞”作用，并迫使其改变迁移路径。地下水位的变动将使污染物的产生重新分布。当地下水位降低以后，污染物产生明显的向下迁移。     

饱和度-压力曲线法预测柴油在毛细带中形成的透镜体厚度

轻非水相液体（LNAPL）进入到地下环境中,首先在非饱和带进行垂向迁移,到达毛细带后在地下水面以上形成饱和透镜体。准确预测LNAPL渗漏在毛细带形成的透镜体厚度,对于LNAPL的去除及污染含水层的修复具有重要的指导意义。通过模拟实验测定了水-油两相吸湿及脱湿曲线,采用水-油饱和度-压力曲线拟合出了吸湿及脱湿进入压力。根据透镜体达到稳定状态时脱湿、吸湿压力水头之间的关系,建立了透镜体厚度的预测计算方法,并以中砂和粗砂为例,预测出实验条件下柴油在毛细带形成透镜体稳定时的厚度分别为4.61 cm和1.29 cm。通过室内模拟槽实验,测得柴油在中砂和粗砂中迁移时形成的透镜体厚度分别为5.30 cm和1.50 cm。预测与实验结果的相对误差分别为13.0%和14.0%。误差产生的主要原因为脱湿与吸湿曲线是在非干燥情况下测得,从而导致预测的透镜体厚度偏小。     

冻土离心模型试验相似准则分析

在研究土壤冻融问题时,土工离心模型试验在时间效应和应力全等型模拟方面更具优势。然而,目前关于冻土离心模型试验的研究成果较少,其中的相似准则也不够完善。针对这一现状,在充分考虑土体的水-热-力耦合作用过程的基础上,采用Butterfield量纲分析法确定了控制饱和土冻融变形性状的无量纲项,建立了冻土离心模型试验中孔隙压力、热扩散、未冻水迁移、融土固结以及冻融变形的相似准则。分析结果表明,冻土离心模型试验中未冻水迁移、融土固结及热扩散效应的非稳态时间具有统一的比尺,即离心模型是原型的1/N2倍,而未冻水迁移流速的比尺为离心模型是原型的N倍。此外,利用冻土离心模拟装置完成了一则对渠道基土冻融作用的离心模拟。     

重金属迁移转化模型研究

大量重金属污染物迁移转化的现象和研究成果都表明重金属以泥沙颗粒为载体迁移转化,描述重金属在天然水体中的迁移转化必须紧紧抓住泥沙颗粒运动及重金属与泥沙之间的转化关系进行.为此在对现有的重金属迁移转化数学模型进行概括分类的基础上,根据水沙运动与污染物相互作用关系,分析了泥沙颗粒运动及重金属吸附解吸不平衡过程,并结合水沙数学模型,建立了重金属迁移转化的耦合模型.同时在模型合理性分析的基础上,对模型进行了计算分析,分析表明模型能够合理地反映重金属污染物在水体中的迁移转化过程.     

湿度对盐溶液在壁画地仗中的毛细迁移影响研究

盐溶液在地仗中毛细迁移引起的盐分迁移、富集已成为壁画病害发生的重要原因。研究不同湿度条件下盐溶液在模拟地仗中的毛细迁移过程对可溶盐再分布的影响。分析不同湿度条件下毛细上升高度随时间的变化关系,含水率和电导率随高度变化关系,毛细饱和后试样表面变化情况及可溶盐含量随试样高度的分布。试验结果表明,空气相对湿度越低,毛细上升过程中水分向空气中迁移越多,水分难以向上迁移。试样的含水率随高度呈递减趋势,而电导率随试样高度增加而增大。可溶盐在毛细水上升过程中发生结晶分异,毛细前锋以NaCl结晶富集为主,亚前锋以Na2SO4结晶富集为主。其研究结果可以为壁画盐害的防治提供基础科学依据。     

离心实验在污染物迁移研究中的应用

离心实验模拟以其可以获得与原型一致的应力水平并且能够大大缩短原型历时而被应用于污染物迁移实验研究。在离心相似理论及离心实验模拟的相似基础上,重点回顾了离心机在饱和带水分、非饱和带水分、保守性溶质、NAPLs、重金属、核素迁移以及污染场地修复方面的实验应用。最后讨论了离心模拟中土壤预制和加速度选择的问题,并简单介绍了离心监测方法。可以得出结论:离心机能够成功用于各类物质迁移的实验研究中,离心实验模拟能够为理论和数值等分析方法快速提供真实可靠的参数依据,但离心实验的理论基础和监测方法需要进一步完善,此外也应积极开展更接近实际情况下的离心实验模拟研究。     

多孔介质界面对重非水相液体迁移过程影响的图像法研究

重非水相液体(dense nonaqueous phase liquid,DNAPL)污染土壤和地下水的问题已引起广泛关注,研究其在不同粒径多孔介质及其界面的运移特征形态是确定污染区域、修复治理土壤和地下水环境的前提。文章通过室内试验研究多孔介质界面对DNAPL运移与分布特性的影响。首先在二维砂槽上进行DNAPL污染物的入渗试验,试验过程中用数码相机拍照,将DNAPL扩散过程以图像的形式记录下来;然后用AutoCAD对图片进行处理,绘制出DNAPL迁移过程的锋面变化图。结果表明:DNAPL入渗过程中,迁移主要受到重力作用与毛细作用的控制,毛细作用力随着介质粒径的增大逐层减小,重力作用逐渐起主导作用使污染物入渗速度逐层增大;介质结构影响DNAPL的迁移形态,介质粒径逐层增大,DNAPL污染物的渗流面与指进扩散宽度逐层减小,扩散方式由面状变为指状;在不同粒径介质界面介质结构发生突变时,DNAPL迁移锋面线曲率也相应变大,此时DNAPL的迁移呈现“凸”型特征,另外,不同的界面横向扩散的滞留宽度不同,随着介质粒径的增大,界面的横向扩散宽度相对变短。     

An elastoplastic strainhardening model for soil allowing for hydraulic bonding–debonding effects

The paper presents a strainhardening constitutive model for unsaturated soil behaviour based on energy conjugated stress variables in the framework of superposed continua. The proposed constitutive law deals with hydro‐mechanical coupling phenomena. The main purpose is to develop within a consistent framework a model that can deal with possible mechanical instabilities occurring in partially saturated materials. The loss of capillary effects during wetting processes can, in fact, play a central role in unstable processes. Therefore, it will be shown that the bonding effects due to surface tensions can be described in a mathematical framework similar to that employed for bonded geomaterials to model weathering or diagenesis effects, either mechanically or chemically induced. The results of several simulations of common laboratory tests on partially saturated soil specimens are shown. The calculated behaviour appears to be in good qualitative agreement with that observed in the laboratory. In particular it is shown that volumetric collapse phenomena due to hydraulic debonding effects can be successfully described by the model. Finally, it will be highlighted the ability of the model to naturally capture the transition to a fully saturated condition and to deal with possible mechanical instabilities in the unsaturated regime. Copyright © 2008 John Wiley & Sons, Ltd.     

Morphological transitions for pore water and pore air during drying and wetting processes in partially saturated sand

Acta Geotechnica - Water retention characteristics are important for modeling the mechanical and hydraulic behavior of partially saturated sand. It is well known that the soil water characteristic...     

Molecular dynamics modeling of a partially saturated clay-water system at finite temperature

The mechanical and hydraulic properties of unsaturated clay under nonisothermal conditions have practical implications in geotechnical engineering applications such as geothermal energy harvest, landfill cover design, and nuclear waste disposal facilities. The water menisci among clay particles impact the mechanical and hydraulic properties of unsaturated clay. Molecular dynamics (MD) modeling has been proven to be an effective method in investigating clay structures and their hydromechanical behavior at the atomic scale. In this study, we examine the impact of temperature increase on the capillary force and capillary pressure of the partially saturated clay-water system through high-performance computing. The water meniscus formed between two parallel clay particles is studied via a full-scale MD modeling at different elevated temperatures. The numerical results have shown that the temperature increase impacts the capillary force, capillary pressure, and contact angle at the atomic scale. The capillary force on the clay particle obtained from MD simulations is also compared with the results from the macroscopic theory. The full-scale MD simulation of the partially saturated clay-water system can not only provide a fundamental understanding of the impact of temperature on the interface physics of such system at the atomic scale, but also has practical implication in formulating physics-based multiscale models for unsaturated soils by providing interface physical properties of such materials directly through high-performance computing.     

Selection of Multifractal Scaling Breaks andSeparation of Geochemical andGeophysical Anomaly

1Multifractalities have been observed and estimated formany quantities in physics, chemistry, as well as earth andspace sciences (see a review in Agterberg, 20011 Cheng,1999al Agterberg et al., 1993 1 Schertzer and Love joy, 1991).The multifractal theory has provided a superior distributiontype for describing the complex measures in comparison withthe ordinary statistical distributions. The former usually involves the multiple moments and parameters (Cheng, 1999alAgterberg et al., 1993 1 Sc…     

Centrifuge modelling of capillary rise

This paper reports results from centrifuge tests designed to investigate capillary rise in soils subjected to different gravitational fields. The experimental programme is part of the EU-funded NECER project (Network of European Centrifuges for Environmental Geotechnic Research), whose objective is to investigate the appropriateness of geotechnical centrifuge modelling for the investigation of geoenvironmental problems, particularly with reference to partially saturated soils. The tests were performed at the geotechnical centrifuge laboratories of Cardiff, Bochum, Manchester, and LCPC in Nantes. The aim was to determine the scaling laws of capillary rise under both equilibrium and transient conditions.

In all laboratories, column wetting tests in fine poorly graded sands (Congleton Sand, Bochum Normsand, HPF5 Sand, and Fontaineblau Sand) were performed. Capillary rise above the phreatic surface of the sand model was distinguished in a continuous capillary zone (completely saturated) and a discontinuous capillary zone (partially saturated).

The Cardiff Geotechnical Centrifuge Laboratory used matrix potential probes to follow the capillary rise of the continuous zone and, therefore, determine the suction above the phreatic zone during centrifuge testing. At Bochum, two cameras were used for optical and volumetric measurements, in order to follow the rise of the visible wetting front (upper limit of discontinuous zone) in the sand within the sample column. At Manchester, the movement of the wetting front was observed by video cameras over periods up to 8 h, whereas in LCPC pore pressure transducers recorded the changes in pressure caused by capillarity.

A simple centrifuge similitude law for capillary rise in these sands has been established and the kinetic phenomena have been measured as a function of the gravitational field. The results from these experiments verify that both the continuous and discontinuous capillary zones are scaled at a factor 1/N whereas the time for rise seems to be scaled at a factor 1/N². This research suggests that capillary phenomena can be modelled using a geotechnical centrifuge. Therefore, centrifuge testing can be a useful tool for future modelling of boundary value problems involving complex transport phenomena.     

An elasto‐plastic three phase model for partially saturated soil for the finite element simulation of compressed air support in tunnelling

This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd.     

A linear constitutive model for unsaturated poroelasticity by micromechanical analysis

This paper extends the Biot theory of poroelasticity from the saturated to unsaturated case. The Biot phenomenological model uses parameters that are easily observable, such as the deformation of porous frame, total stress, pore pressure, and fluid specific discharge. Such model is preferred for engineering applications. At this macroscopic level, the extension of Biot theory from saturated to unsaturated is straightforward. The constitutive constants, however, are combined properties of solid, pore space, and fluids. In the unsaturated case, the constants are functions of the degree of saturation. Their measurements and tabulation over a range of saturation is generally not feasible for practical applications. In this work, a Biot-Willis type analysis is performed for the unsaturated case to provide a theory that the bulk material constants can be evaluated using laboratory measurable micromechanical constants under saturated condition, plus a capillary pressure curve (saturation versus suction pressure) typically available for unsaturated porous medium, without the need of measurement at each state of saturation. In particular, it is demonstrated that the surface energy contained in the meniscus interface manifests as a “capillary modulus,” given by the negative inverse slope of the capillary pressure curve. A rigorous analysis based on the thermodynamic variational energy approach is also conducted to lend theoretical support to the phenomenological approach. The presented model can bring a closure to the practical engineering modeling of the deformation of partially saturated porous medium that lacks the information of material constants over the range of saturation.     

Application of In-situ Lysimetric Studies for Determining Soil Hydraulic Conductivity

Several techniques such as laboratory column studies, under normal and accelerated gravity environments, numerical modeling and pedo-transfer functions have been employed by previous researchers to determine hydraulic conductivity of soils. However, these methods suffer with several limitations, particularly, as far as regeneration of in-situ soil conditions and boundary conditions, in an exceptionally small model, are concerned. Under these circumstances, in-situ lysimeters are found to be quite useful for conducting the investigations to obtain hydraulic conductivity of fully saturated or partially saturated soils. These lysimeters can easily be installed in the soil mass, without disturbing its state, and are large enough to exhibit representative soil mass that exist at the site. With this in view, an attempt was made in this study to determine hydraulic conductivity of the soil mass beneath a waste disposal site in India by installing a lysimeter. Details of the instrumentation are presented in this paper along with the methodology to determine hydraulic conductivity of the soil mass.