2-D Simulation of Non-isothermal Fate and Transport of a Drip-applied Fumigant in Plastic-mulched Soil Beds. I. Model Development and Performance Investigation

Pre-plant application of toxic fumigants to soil beds covered by plastic film is commonly used in agriculture to control soil-borne pathogens. Plastic mulch covers tend to physically suppress the emissive loss of gaseous fumigant to the atmosphere. When liquid fumigant metham sodium (MS) is applied in irrigation water to field soil, it is rapidly transformed to the gaseous methyl isothiocyanate (MITC). The gaseous MITC is a potential atmospheric contaminant, and any untransformed MS is a potential contaminant of underlying groundwater due to the high water solubility of MS. A finite element numerical model was developed to investigate two-dimensional MITC fate/transport under non-isothermal soil conditions. Directional solar heating on soil beds, coupled heat and water flow in the soil, and non-isothermal chemical transport were included in the model. Field soil data for MITC distribution, soil water content, meteorological data, and laboratory data were used to verify the model for soil beds covered with plastic mulch. Four possible scenarios were considered: low and high drip-irrigation rates and low and high water contents. The movement of the center of MITC mass in the soil profile was effectively simulated. The lower drip-irrigation rate of MS yielded more extensive coverage of MITC in the plastic-covered soil bed. The lower soil air contents due to higher soil water contents for the higher irrigation rate resulted in high concentrations of soil MITC. NRMSE (normalized root mean square error) calculations further verified that the model predicted fumigant fate/transport well under these non-isothermal field conditions.     

建筑物内火灾烟羽流和顶蓬射流诱发的空气运动盐水模拟实验研究

盐水在清水中的扩散与火灾烟气在空气中的蔓延相似，因而可用盐水的扩散来模拟火灾烟气的蔓延；同样地，由于盐水的扩散而导致的环境清水的运动也可用来模拟烟气蔓延诱发的空气运动。基于这一原理，本文用实验的方法，研究了建筑物内初起火灾产生的烟羽流和顶蓬射流所诱发的空气运动，并对走廊空气层生长规律和运动特性进行了初步的测量和分析。所得结果基本合理，对建筑物通风排烟及疏散通道设计有帮助，从而为建筑火灾环境空气的运动的实验研究开辟了新的思路，也为火灾区域模拟提供了实验依据     

Systematic drainage of irrigated soils underlain by an impermeable layer

Problems of flow through porous media in which irrigation and drainage are simultaneously taken into account are very often encountered in connection with reclamation projects in arid zones. This article investigates the percolating flow in a soil of finite thickness under conditions of seepage and horizontal closed drainage. The soil is drained by a system of horizontal parallel equidistant circular pipes buried at the same depth. It is assumed that the drains are completely full of water; the pressure head in the soil along the line of contact with the drain is equal to the elevation of its upper point, at which there is no excess pressure. At the surface of the soil the moisture content is constant. This may be the result of a nonuniform irrigation rate distribution. If the irrigation is uniform, the surface layer of soil is nonuniformly moistened: maximum saturation, often supersaturation, occurs midway between the drains, while over the drains the soil is practically dry. The corresponding problem of steady state percolating flow reduces to a boundary-value problem of the theory of analytic functions, by solving which it is possible to determine the irrigation rate distribution needed to ensure uniform humidity of the soil surface. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 66–71, November–December, 1986.     

PS溶液渗透对非饱和土的增强效应研究

特定模数硅酸钾溶液(PS)已广泛的用于我国西北干旱地区土遗址加固保护工程中,并且取得了显著的效果。但是PS对潮湿地区遗址土加固效果,目前尚缺乏成熟的理论和实验验证。本文从室内试验的资料出发,研究PS对潮湿地区土遗址非饱和土的增强效应。本次研究土样取自位于潮湿环境中的浙江杭州良渚土遗址(3000 B C.),在实验室内对高含水量(13％ ~25％)的重塑土样进行了PS渗透加固。通过研究对比加固前后重塑土样的固结特性和抗剪强度,发现 PS对潮湿地区遗址土也有一定的加固效果。加固效果随初始饱和度的不同而不同,但初始含水量达25%的土在PS作用下仍可能得到增强。经PS渗透的土样不仅提高了初始固结性能和抗剪强度,而且加固土样而且呈现出一定的弹脆性特征,呈现出塑脆性转化,表明PS渗透可能产生对土初始结构的稳固作用。这一探索研究将为扩大PS在土遗址加固保护工程的应用提供一定的理论支持。     

Role of NAPL Thermal Properties in the Effectiveness of Hot Water Flooding

Hot water flooding is a thermal nonaqueous phase liquid (NAPL) recovery technology originally developed in the petroleum industry that has recently been proposed for enhanced recovery of NAPLs in the contaminated subsurface. This technology, however, has received relatively little laboratory or numerical modeling investigation in the contaminant hydrology community. In this study the utility of flooding NAPL contaminated source zones at elevated water temperatures was investigated. Simulations were conducted using 16 different geostatistical representations of an actual field site. Two NAPLs were selected for this study—a light NAPL with hydraulic properties that have moderate temperature dependencies and a dense NAPL with significant viscosity temperature dependency. For these two NAPLs, flooding the source zone with water at elevated temperatures resulted in enhanced NAPL recovery. However, injection of hot water also resulted in accelerated downward movement of coal tar DNAPL due to the reduced viscosity at elevated temperatures. NAPL recovery was also dependent on the source zone architecture with greater NAPL mass recovery when the NAPL was localized in a small volume at high saturations. These results suggest that hot water flooding can significantly speed up the recovery of viscous NAPLs and, as such, is a powerful technique for the remediation of viscous NAPLs.     

盐渍土工程与力学性质研究进展

介绍我国西部地区大开发高等级公路、铁路网建设中盐渍土病害的力学相关问题. 综述了对盐胀、溶陷和腐蚀3大工程危害的机理研究以及对盐渍土的工程及力学特性的研究进展, 回顾了盐渍土水、盐、热的数学运移模型的发展, 盐渍土水、盐、热、力4场耦合模型的数值与实验探索, 对盐渍土研究领域中亟待进一步研究的问题进行了展望.     

Electrokinetic Salt Removal from Porous Building Materials Using Ion Exchange Membranes

The removal of salt from porous building materials under the influence of an applied voltage gradient normally results in high pH gradients due to the formation of protons and hydroxyl ions at the electrodes. The formed acidic and alkaline regions not only lead to disintegration of the porous material, but also affect the salt transport. In this work we use ion exchange membranes between the electrodes and the porous material to prevent the protons and hydroxyl ions from intruding into the material. The porous material used in this study is fired clay brick, which has been saturated with a 4?mol/l sodium chloride solution prior to the desalination treatment. In order to experimentally determine the salt removal, we monitored the sodium ion concentration profiles across the material with nuclear magnetic resonance (NMR). In addition, we present theoretical predictions for the salt removal according to a model based on the Poisson?CNernst?CPlanck theory for ion transport. From the work reported here, we can conclude that the use of ion exchange membranes to desalinate porous building materials is not useful since it reduces the salt removal rate to such an extent that desalination with poultices, which is driven by diffusion only, is more efficient. The reason behind this is twofold. First, the ion exchange membranes provide a penalty for the ions to leave the material. Second, in the absence of acidic and alkaline regions, the salt concentration at the edges of the porous material will reduce to almost zero, which leads to a locally increased electrical resistance, and thus a reduction of the electrical field in the bulk of the material. Due to this reduction the effect of the applied voltage gradient across the material vanishes, and the salt removal is limited by diffusion.     

粉煤灰作为排水垫层筑坝的动力分析

淮南某电厂灰坝设计拟采用粉煤灰作为坝基排水垫层筑坝以节约资金。为研究其可行性,在已有的大量静力和动力特性试验研究基础上,对该灰坝在地震荷载作用下的动力反应进行研究。研究中以Biot波动方程为基础,考虑水土动力耦合作用和土体和非线性性质,引入地震孔隙水压力和残余变形项,把动力分析与动力渗流、液化发展及骨架变形加以综合考虑。研究表明,7°地震烈度条件下用本文配制的粉煤灰代替砂作为排水垫层筑坝是可行的。     

膨胀土胀缩性评价中有关问题的研究

膨胀土吸水膨胀、失水收缩，产生地表变形，进而导致建筑物破坏。本文就现行膨胀土胀缩性指标试验及地基变形量评价计算中有关问题进行了研究，认为以天然含水量为起始含水量测定土样胀缩率是不合适的，建议以缩限为起点进行膨胀试验，以胀限为起点进行收缩试验，采用加荷膨胀量（而非卸荷膨胀量）参加有荷载下的膨胀率计算。膨胀土地基勘察必须结合场地条件分析可能影响土体含水量变化的因素，对同一工程采用几种方法评价，互相对比、参照。     

The numerical analysis of water and solute movement in scale heterogeneous profiles

A computer based numerical method is presented for the analysis of water and solute movement in unsaturated heterogeneous porous materials. Such a method is necessary since, for those field studies where solute movement is of concern, the soil profiles under consideration are invariably heterogeneous. The numerical analysis is based on a general one-dimensional finite difference soil water flow model which includes a numerical technique combining the concepts of scale heterogeneity with an interpolative soil water hysteresis model. An explicit finite difference solute movement subroutine is incorporated into the unsaturated flow model to describe the transport of nonreactive solutes. A velocity dependent longitudinal dispersion coefficient is used in the solution of the hydrodynamic dispersion equation. The resulting hysteretic scale heterogeneous solute movement model permits the study of solute dynamics during infiltrating and redistribution in realistically complex spatially varying soil profiles. Results are presented for the leaching of both coarse grading to fine and fine grading to coarse sand profiles. Both vertical and horizontal profiles are studied using either a constant flux or a constant concentration input boundary condition. The four cases studied demonstrate the versatility of the numerical method and emphasise the substantial differences in transport behavior that can arise between heterogeneous and homogeneous profiles.Now with BHP Petroleum Pty. Ltd., GPO Box 1911R, Melbourne, Vic. 3001, Australia.     

Time-Dependent Measurement of Strongly Density-Dependent Flow in a Porous Medium via Nuclear Magnetic Resonance Imaging

Solute transport in saturated artificial porous media was observed in a series of laboratory experiments using magnetic resonance imaging. The objective was to study a situation of density-dependent flow in three dimensions both qualitatively and quantitatively. The time-dependent measurements visualised inflow from below of dense salt water into a freshwater reservoir, internal density-driven flow and the behaviour of a salt water layer below freshwater flow including plume development by dispersion. The main feature of the flow experiment was the strong tendency for the salt water to remain stagnant and to resist being swept out by the freshwater. Additional measurements were performed to gain information about reproducibility, flow field and breakthrough curves.     

高压驱动研究方法的改进及其对质点运动的研究

从分析二维运动的测试原理入手 ,利用高速狭缝扫描摄影技术建立了测试系统 ,该系统能获得二维运动的轨迹和速度 ,从而把高速狭缝扫描摄影的研究领域由一维拓展到了二维 ,并对不同种测试系统进行了比较。利用新的测试系统获得了物体运动的照片 ,研究了质点高压下的运动 ,同时 ,获得了变形速度和变形角。     

真空预压法加固吹填土的孔隙水压力试验研究

利用真空预压法处理吹填土时,孔隙水压力变化常常反映土体固结程度的好坏。通过6个模型箱试验,监测不同排水系统下孔隙水压力变化,确定有效排水体间距。研究发现0.4m间距的土内孔压下降效果比0.8m间距的土内孔压下降效果好; 排水体内的孔隙水压力与排水体类型有关,且距离排水体10cm处土体内的孔隙水压力仅为排水体内孔隙水压力的1/2弱; 滤膜排水系统中的吹填土孔隙水压力下降幅度最快,B型排水板系统次之,而砂井系统最慢。另外,对于吹填土而言,排水体有效间距介于0.4m与0.8m之间,其中滤膜的有效间距最大,B型排水板次之,砂井远小于前两者。     

Effects of Fine Particles on the Suffusion of Cohesionless Soils. Experiments and Modeling

Internal erosion has long been a major problem associated with earthen structures, and its field investigation has been limited because of its complexity. Laboratory experiments provide a potential insight into the induced processes. Soils suitable to suffusion involve an easy movement of fine particles between the coarse ones. In this study, suffusion experiments were performed on a laboratory column packed with different sand–fines mixtures and subjected to controlled flow conditions. The mixtures are made of Fontainebleau sand (NE34) and fines (kaolinite or illite or silt). The initiation and the progression of soil suffusion were investigated for soil mixtures, and the hydraulic conductivity variation was derived from pressure measurements. Attempts are made to assess the influence of the type and content of fine particles on mixtures suffusion. The results show that among the tested fine particles, illite exhibits a great resistance against suffusion. Fine content increase involves less suitability to suffusion till a threshold value. Recorded pressures indicate that detachment of fines may be followed by clogging in the soil matrix, leading to a decrease of hydraulic conductivity. Furthermore, analytical solutions of simplified mathematical model based on the mass conservation of eroded particles (Govindaraju et al., J Hydrol 172:331–350, 1995) were used to simulate the recorded soil suffusion curves, and the model performance was evaluated. Comparisons between numerical and experimental results indicated a quite good agreement, and allowed the estimation of the model parameters.     

大庆地区湿陷性黑土研究

大庆南部地区分布着一种颜色发黑、结构松散、成分复杂、对水极度敏感的工程性质极差的地质材料 ,它兼有湿陷性土、有机质土、碱土和分散土等诸多特殊土的劣质特性。通过野外调查、原位测试和室内各种试验研究 ,对这种黑土的野外特性、成分结构、物理力学性质、承载力和湿陷特性、地层年代等进行了较系统的阐述和分析 ,并通过综合比较 ,将其定名为湿陷性黑土.     

Estimability Analysis for Optimization of Hysteretic Soil Hydraulic Parameters Using Data of a Field Irrigation Experiment

To improve the quality of parameter optimization, estimability analysis has been proposed as the first step before inverse modeling. When using field data of irrigation experiments for the determination of soil hydraulic parameters, wetting and drying processes may complicate optimization. The objectives of this study were to compare estimability analysis and inverse optimization of the soil hydraulic parameters in the models with and without considering hysteresis of the soil water retention function. Soil water pressure head data of a field irrigation experiment were used. The one-dimensional vertical water movement in variably saturated soil was described with the Richards equation using the HYDRUS-1D code. Estimability of the unimodal van Genuchten–Mualem hydraulic model parameters as well as of the hysteretic parameter model of Parker and Lenhard was classified according to a sensitivity coefficient matrix. The matrix was obtained by sequentially calculating effects of initial parameter variations on changes in the simulated pressure head values. Optimization was carried out by means of the Levenberg-Marquardt method implemented in the HYDRUS-1D code. The parameters \(\alpha , K_{s}, \theta _{s}\) , and \(n\) in the nonhysteretic model were found sensitive and parameter \(\theta _{s}\) strongly correlated with parameter \(n\) . When assuming hysteresis, the estimability was decreased with soil depth for \(K_{s}\) and \(\alpha ^{d}\) , and increased for \(\theta _{s}\) and n. Among the shape parameters, \(\alpha ^{w}\) was the most estimable. The hysteretic model could approximate the pressure heads in the soil by considering parameters from wetting and drying periods separately as initial estimates. The inverse optimization could be carried out more efficiently with most estimable parameters. Despite the remaining weaknesses of the local optimization algorithm and the inflexibility of the unimodal van Genuchten model, the results suggested that estimability analysis could be considered as a guidance to better define the optimization scenarios and then improved the determination of soil hydraulic parameters.     

饱和非饱和土壤中溶质运移的数值模拟

本文结合拉格朗日观点的运动坐标与欧拉观点的固定坐标,分别求解对流—弥敌方程的对流效应和弥散效应,并在从固定坐标向运动坐标投影时利用了三次样条插值方法,使模型更为精确而简单。通过与特定情况下的理论解比较,说明模型是准确可靠的。通过与实验室室内试验结果的比较,模型的可靠性得到了进一步的验证。利用模型对淡水洗盐结合浅群井抽排咸水治理盐碱地过程的模拟分析,认为结合浅群井抽排可较大地提高淡水洗盐的效果。     

矿坑涌水量的地质分析与模拟──以三山岛金矿为例

本文充分利用地表及地下井巷揭露资料，配合物探成果推测深部水文地质条件。利用同位素 ¹⁸O分析结果量化各种水体补给量，然后利用数值模型检验、校正上述成果，优选水文地质参数并进一步预测深部中段涌水量。考虑到模拟过程中某些未知因素，或未来开拓过程中工程措施的影响，本文对预测结果进行了必要的讨论和修订。     

动力排水固结法加固吹填黏性土的模型试验研究

为了研究动力排水固结法加固吹填土地基的效果,研发了大比例尺的室内模拟试验装置并开展了试验。重点对夯沉量、土体内孔隙水压力的瞬时和长期变化进行了监测和分析,得到如下结论:夯击产生的能量以波的形式向土体内部传播,由于波的作用造成土体拉伸形成裂隙,这些裂隙成为初期的排水通道; 孔隙水压力的动态变化具有明显的脉动特性,夯击瞬间出现负增长,然后再累积增长; 不同深度土体的孔隙水压力增幅不一样; 土体表面铺设的砂垫层不仅可以用做排水,还可以作为静力荷载,有助于土体的排水固结。     

渗透固结试验装置的改进

针对传统固结仪不能进行渗透试验的不足,对固结仪进行了改装。改装后的固结仪不仅可以交叉进行固结、渗透试验,以研究不同荷载作用下土样的渗透特性,而且可以通过三通阀控制排水条件,进行双面排水固结试验,或在单面排水条件下量测固结过程中试样不排水面处的孔隙水压力,也可通过串联固结仪分别研究试样不同截面处孔隙水压力的消散规律。新仪器具有构造简单,密封性能好的特点。通过与常规固结、渗透试验的结果对比,验证了采用改装固结仪测定固结模型参数的可靠性。另外,通过单面排水固结试验发现了孔压滞后现象。