Characterizing the hydraulic properties of unsaturated sand considering various porosities and drying‐wetting paths of the retention curve

Previous studies have shown that water retention curve (WRC) and the hydraulic conductivity vary because of changes of the void ratio or porosity of soil. However, limited documents pointed out the change of hydraulic properties of soil when compacted to different porosities while considering both of the drying and wetting processes of the WRC. This information is sometimes necessary for research like finger flow analysis or the occurrence of wetting and drying cycles as what would be seen in the field. Therefore, this study aims to examine the change of WRC characteristics with varied porosity considering both of the drying and wetting path in WRC by conducting a sand box experiment. Results show that the same type of sand compacted to various porosities have different hydraulic parameters. Hydraulic conductivities generally decrease with reduced porosities; shape parameter α of the van Genuchten equation (1980) linearly decreases with declining porosity and shape parameter n in a reversal manner for the sands of interest whether in the drying process or wetting process. The unsaturated properties of sand are further characterized by inspecting the variations of moisture content, matric suction and vertical displacement of soil body subject to periodic changes of the water level by another sand box experiment. The outcomes suggest that the saturated water content and residual water content are changing during the wetting–drying process, which can be an implication of the changed properties of WRC. The characteristics of volumetric deformation might be varied as well because of the observation of the dissimilar patterns of the changing vertical displacements among each wetting–drying process. Infiltration patterns of the sands also are identified through numerical modelling by introducing a constant infiltration flux from the surface followed by a no‐influx condition. Results indicate that less water accumulates in the sand near the surface for the sand compacted to higher porosity, but water can move deeper. Hydraulic conductivity is found as the prime factor dominating the evolvement of wetting fronts. However, shape parameters of water retention curves also affect the infiltration pattern to some extent. In addition, different sands with similar porosities can have quite different infiltrating characteristics. Copyright © 2012 John Wiley & Sons, Ltd.     

VARIABILITY OF MELTWATER AND SOLUTE FLUXES FROM HOMOGENEOUS MELTING SNOW AT THE LABORATORY SCALE

Four experiments were performed to examine the relationship between the meltwater flow field and ion release from melting snow. A 0.4 m³ volume of snow was placed in a Plexiglass box and melted from above using a heating plate. The meltwater and solute fluxes issuing from the bottom of the snow were monitored. In experiments with NaCl tracer added to the snow, the solute concentrations were generally lower in the flow fingers than in the background wetting front. Dye tracer experiments revealed contemporaneous areas of concentrated dye and dilute meltwater in flow fingers. This suggests that the meltwater in flow fingers is diluted by low concentration water from the top of the snowpack. Flow fingers contribute more meltwater flux primarily because the flow is maintained for a longer period of time than in the non-finger areas; however, the relative contribution of flow fingers to solute flux was apparently not as great as that of the background wetting front because of dilution of solute in the flow finger areas.     

Fingered flow: The creator of sand columns in dune and beach sands

Sand columns, sand cones, sand mushrooms and other striking sand forms are frequently observed in the Dutch and German beach and dune sands. This paper aims to clarify the mechanism of sand column formation. Recently it has become evident that homogeneous beach and dune sands often become irregularly wetted by infiltrating rainwater. In otherwise dry sandy soils, wet preferential flow paths (‘fingers’) may develop. At two test sites the volumetric soil moisture content varied between 0·2 and 12·0 per cent. The wet fingers represent the premature state of sand columns. When the dry sand in between these fingers is blown away by the wind, the more resistant wet sand of the fingers will remain in its place and appear as sand columns at the surface. As a result of wind and erosive sand drifts, striking sand forms may be formed.     

部分饱和砂岩储层地震物理模拟及含气饱和度预测分析

含气饱和度预测是天然气储层地震解释工作的重要目标.本文将岩石物理分析与地震物理模拟技术相结合,构建了部分;饱和砂岩储层物理模型并进行含气饱和度预测分析.物理模型中设置了高孔渗常规砂岩和低孑孔渗致密砂岩两种模拟储层,每种储层都是由具有不同含水饱和度的气-水双相饱和砂体组成.岩石物理分析结果显示在低孔渗致密砂岩中气-水混合流体更加倾向于非均匀的斑块分布,而结合了Brie等效流体公式的Gassmann流体替换理论可以更准确地描述纵波速度随含水饱和度的变化趋势.对物理模型进行地震资料采集处理后,对比了AVO特征和叠前同步反演结果对两种砂岩储层含气饱和度预测能力的差异.AVO特征结果显示,对于混合流体均匀分布的高孔渗砂岩储层,AVO响应曲线和属性变化很难对含气饱和度进行估算;对于混合流体斑块分布的致密砂岩储层,AVO特征可以定性地分辨出储层是否为高、中、低含气情况.反演结果显示,密度及纵横波速度比分别对高孔渗及致密砂岩储层的含气饱和度有着较好的指示能力.     

Using quartz grain size and shape analysis to distinguish between aeolian and fluvial deposits in the Dallol Bosso of Niger (West Africa)

Size and Fourier-shape characteristics of quartz sand grains were determined by computerized image analysis in order to distinguish between aeolian and fluvial soil parent materials in the Dallol Bosso in Niger. Factor analysis of grain-size distributions gave four sand end-members that can be related to fluvial transport dynamics operating when the sediments were initially deposited. The medium to fine (and more angular shaped) sand fractions are being reworked by wind. Aeolian deposits were well sorted whereas fluvial deposits were poorly sorted in both size and shape. Although gross-shape characteristics (lower harmonics of Fourier series expansion) indicated a common source rock for all sands, the aeolian sands were well rounded whereas the fluvial sands tended to be more angular (upper harmonics of Fourier series).     

Air and water entrapment in the vicinity of the water table

A laboratory tank was used to study entrapment of water in coarse sand lenses above the water table and of air in coarse sand lenses below the water table. Monitoring of these experiments involved a combination of visual inspection, measurement of moisture content, and measurement of air/water pressure. The medium consisted of coarse sand lenses with various degrees of vertical connectivity embedded within a fine sand matrix. Experiments were performed under conditions of both drainage (from a fully saturated medium) and imbibition. Observations during drainage included: (1) water was trapped in the coarse sand zones above the water table at heights significantly greater than anticipated from consideration of capillary rise in the coarse sand; (2) rapid drainage of these same coarse zones occurred when air penetrated into these zones through the surrounding fine sands; and (3) prior to the time of penetration of the coarse sand by air, water pressure in the coarse zone dropped significantly below atmospheric pressure. Observations during imbibition included: (1) entrapment of air within coarse sands below the water table, (2) the pore fluids in these zones varied spatially from predominantly air to predominantly water, and (3) pressure in the trapped air phase was significantly greater than pressure in the water phase in the surrounding fine sand. Overall, these results demonstrated significant sensitivity to the geometry of the coarse sand inclusions, particularly the vertical connectivity of the coarse sand lens.     

不同泥质分布形式泥质砂岩导电规律实验研究

本文利用人工制作的不同含量分散泥质和层状泥质砂岩岩心样品,测量不同矿化度和不同含油饱和度的岩心电阻率,从实验角度研究了不同泥质分布形式和含量的岩心导电规律,结果表明,泥质分布形式或含量不同,则泥质砂岩导电规律不同.基于层状泥质与分散泥质砂岩的并联导电实验规律,以及分散粘土和地层水混合物的导电规律可用HB电阻率方程描述,建立了考虑泥质分布形式影响的泥质砂岩电阻率模型.通过1组不同泥质分布形式泥质砂岩人造岩心实验数据的测试,表明该模型可以描述分散泥质砂岩、层状泥质砂岩和混合泥质砂岩地层的导电规律.分散泥质,层状泥质,人造岩样,实验测量,并联导电,HB方程,电阻率模型     

Electrical Resistance Heating of Clay Layers in Water‐Saturated Sand

Electrical resistance heating (ERH) experiments were performed in a two‐dimensional water‐saturated porous medium comprising an electrically conductive, low‐permeability clay lens embedded within a less electrically conductive, higher permeability silica sand. These were compared to experiments performed in homogeneous silica sand. All experiments were performed in the absence of a non‐aqueous phase liquid (NAPL) or dissolved volatile organic compound (VOC). Temperature monitoring showed preferential heating in the clay lens and higher overall heating rates throughout the test cell compared to the homogeneous case. Gas production was localized around the sand–clay interface due to high temperature and low capillary displacement pressure. Above the clay lens, unexpected temperature plateaus were observed, similar to those observed in previous experiments during NAPL–water co‐boiling. A conceptual model based on the consumption of thermal energy as latent heat of vaporization in the highly localized heating and gas production region adjacent to the clay lens is proposed to explain the temperature plateaus. Supporting data is drawn from images of the gas phase and electric current measurements. These results show that the use of co‐boiling plateaus as an indicator of NAPL–water co‐boiling could be misleading during applications of ERH at sites containing electrically conductive, low‐permeability clay lenses embedded within less electrically conductive, higher‐permeability sands.     

Measurement of small strain shear modulus of clean and natural sands in saturated condition using bender element test

Bender element (BE) tests of saturated sand have increased interest to researchers currently. However, the measurement of small strain modulus from BE tests shows large difference between saturated and dry conditions. In this study, BE tests of a type of clean sand (Fujian sand) and two types of natural sands (Hangzhou sand and Nanjing sand) were performed. For the purposes of comparison, resonant column (RC) test and torsional shear (TS) test were also carried out on the same specimen. The factors that influence the determination of the travel time of shear wave in BE tests are discussed and a reliable method for the determination of the shear-wave velocity is obtained. It is found that the shear-wave velocities V_s of saturated Fujian sand (clean sand) and Hangzhou sand (natural sand) obtained from BE tests are 5–10% greater than those obtained from RC and TS tests. However, the V_s of saturated Nanjing sand (natural sand) obtained from BE, RC and TS tests show good agreement with a maximum difference of about 3%. Sands with various fines contents were also tested in an attempt to explain the differences between the two saturated natural sands. Biot׳s theory accounting for the dispersion of shear wave was employed to interpret the results of BE tests. The results indicate that the fines content of natural sand plays an important effect on the hydraulic conductivity, which affects the relative motion between soil particles and fluid when a high frequency shear wave propagates in the specimen. Based on this, a method for the determination of small strain shear modulus in BE test was proposed for both saturated clean sands and natural sands.     

热传导测井方法检测油气

本文通过对孔隙介质中热源传热的数学物理方程分析，并结合沙桶中热源传热物理试验，发现在对储层加热过程中，油气层中点热源处的温度与水层（或干层）中点热源处的温度相比增温的速度明显地快，而在加热过程后，油气层中点热源处的温度与水层（或干层）中点热源处的温度相比降低的慢。据此提出用人工热源（炸药热源、电热源等）对地层进行加热，然后测量井中温度的动态变化，据此思想，提出了一种从孔隙介质传热出发解释油气层、预测油气产能的传热测井新方法。     

Test of rock physics models for prediction of seismic velocities in shallow unconsolidated sands: a well log data case‡

This paper tests the ability of various rock physics models to predict seismic velocities in shallow unconsolidated sands by comparing the estimates to P and S sonic logs collected in a shallow sand layer and ultrasonic laboratory data of an unconsolidated sand sample. The model fits are also evaluated with respect to the conventional model for unconsolidated sand. Our main approach is to use Hertz‐Mindlin and Walton contact theories, assuming different weight fractions of smooth and rough contact behaviours, to predict the elastic properties of the high porosity point. Using either the Hertz‐Mindlin or Walton theories with rough contact behaviour to define the high porosity endpoint gives an over‐prediction of the velocities. The P‐velocity is overpredicted by a factor of ～1.5 and the S‐velocity by a factor of ～1.8 for highly porous gas‐sand. The degree of misprediction decreases with increasing water saturation and porosity.Using the Hertz‐Mindlin theory with smooth contact behaviour or weighted Walton models gives a better fit to the data, although the data are best described using the Walton smooth model. To predict the properties at the lower porosities, the choice of bounding model attached to the Walton Smooth model controls the degree of fit to the data, where the Reuss bound best captures the porosity variations of dry and wet sands in this case since they are caused by depositional differences. The empirical models based on lab experiments on unconsolidated sand also fit the velocity data measured by sonic logs in situ, which gives improved confidence in using lab‐derived results.     

Velocity and attenuation of compressional waves in nearly saturated soils

Based on the two-phase theory of Biot, we present exact and approximate expressions for the velocity and attenuation of compressional waves within nearly saturated poroelastic media. We use the approximate solutions to model the low-frequency compressional waves within nearly saturated soils. The model accounts for the effective stress, degree of saturation, and void ratio, and is capable of describing experimental results on Ottawa sand. The three-phase theory of Vardoulakis and Beskos and the two-phase theory of Biot similarly describe the velocity and attenuation of compressional waves in most soils. However, the former theory breaks down for nearly saturated gravels and dense sands.     

深水钻探面临的挑战:浅水流灾害问题

浅水流(shallow water flow,简称SWF)灾害被看作深水钻探中面临的重要挑战之一,它是指深水钻探中,钻头钻过一超压砂层时,孔隙压力太大以致砂和水激烈流进到井眼里,导致了井和钻井平台损坏的事件.SWF问题仅发生在深水区,位于海底泥线下几百米以下,主要出现在一种原位超压的未固结的斜状砂体中,这种砂体一般被低渗透的泥覆盖.超压的形成常认为是快速沉积压实所致,也可能与构造因素或水合物分解有关.SWF砂体具有低密度、低速度和高Vp/Vs的性质,是可能被地震勘探方法所检测的.我国即将进入深水钻探领域,为了避免或减轻深水钻探中可能面对的SWF灾害,我们建议开展有关SWF灾害问题的前期研究.     

High strain dynamic modulus and damping of chemically grouted sand

Laboratory resonant column and cyclic triaxial tests are performed to determine the dynamic response (i.e. shear modulus and damping) of chemically grouted sand. The effect of chemical grouting is evaluated as a function of shearing-strain amplitude, confining stress, cycling prestrain, number of cycles, grout type, concentration, and curing time. The test results show that the shearing-strain amplitude, grout type and grout concentration have significant effects on the shear modulus and damping ratio of the test specimens. The increased addition of sodium silicate grout, which produces stiff gels, improved the shear modulus of the test sand. The acrylate (AC-400) and polyurethane (CG5610) grout, which produces flexible (rubber-like) gels, improved the damping capacity of the sand with increasing grout concentration. The addition of chemical grout greatly reduces the effect of cyclic prestraining over untreated sands. In the case of dense sands, the reduction of cyclic prestraining is less pronounced than in loose sands, which have a higher potential for particle movement and reorientation.     

地震作用下不同厚度饱和砂土中直群桩结构动力响应试验研究

基于振动台试验,设计制作2×2直群桩结构相似模型,通过输入一定峰值加速度的迁安波,对非液化砂土、300 mm厚饱和砂土和380 mm厚饱和砂土中群桩承台横向动力响应特性展开研究。研究结果表明:在迁安波输入下,非液化砂土中群桩承台加速度和位移时程与台面输入时程相比,其波形变化规律与峰值大小没有明显差异;而对于两种不同厚度饱和砂土中承台加速度放大较多,承台位移峰值较台面位移峰值相差不大。在300 mm饱和砂土中群桩承台加速度峰值较台面输入放大了约1.35倍,在380 mm饱和砂土中承台加速度峰值放大了1.42倍,说明在相同输入条件下,较厚的饱和砂土层在发生液化后群桩承台的动力响应更加显著。     

Liquefaction characteristics of gap-graded gravelly soils in K0 condition

A series of undrained cyclic direct simple shear tests, which used a soil container with a membrane reinforced with stack rings to maintain the K₀ condition and integrated bender elements for shear wave velocity measurement, were performed to study the liquefaction characteristics of gap-graded gravelly soils with no fines content. The intergrain state concept was employed to categorize gap-graded sand–gravel mixtures as sand-like, gravel-like, and in-transition soils, which show different liquefaction characteristics. The testing results reveal that a linear relationship exists between the shear wave velocity and the minor fraction content for sand–gravel mixtures at a given skeleton void ratio of the major fraction particles. For gap-graded gravelly sand, the gravel content has a small effect on the liquefaction resistance, and the cyclic resistance ratio (CRR) of gap-graded gravelly sands can be evaluated using current techniques for sands with gravel content corrections. In addition, the results indicate that the current shear wave velocity (V_s) based correlation underestimates the liquefaction resistance for V_s values less than 160 m/s, and different correlations should be proposed for sand-like and gravel-like gravelly soils. Preliminary modifications to the correlations used in current evaluations of liquefaction resistance have thus been proposed.     

Directions of preferential flow in a hillslope soil, 1. Quasi‐steady flow

Preferred infiltration is mainly perceived as vertically down whereas subsurface storm flow is thought to occur parallel to slopes. The transition from vertical to lateral flow in a layered hillslope soil is the focus of the contribution. Transient flow is assumed to move as a wetting front. Three time‐domain reflectometry (TDR) wave‐guides, each 0·15 m long, were mounted in the shape of a truncated tetrahedron with its peak pointing down. Each wave‐guide focuses the front velocity along its axis. The three front‐velocity vectors are decomposed into their x, y and z components, which are then assembled to the resultant velocity vector. The volume density flux of preferred flow is the product of the front velocity and the mobile water content. The latter is the amplitude of transient soil moisture measured with each wave‐guide. The resultant vector of the volume flux density is computed similarly to the velocity vector. The experimental approach allows for the rapid assessment of transient flows without relying on the variation of water potentials. The experiments indicate that the directions of the resultant vectors of velocity and volume flux density can be estimated if the moisture variations of the three TDR wave‐guides are strongly correlated during the passing of the wetting front. Copyright © 2004 John Wiley & Sons, Ltd.     

Volcanic sands of Iceland ‐ Diverse origins of aeolian sand deposits revealed at Dyngjusandur and Lambahraun

The origin, formation and evolution of volcanic sands are less well known than the formation of the much more common quartz‐rich sand sheets. Combining active volcanism and a cold climate, Iceland is covered for about 21% of its surface by sandy areas. The sands were analyzed in detail at two sites and results reveal their diverse origins. The first site is Dyngjusandur, located north of Vatnajökull, and the second site is the Lambahraun area, located south of Langjökull. At both sites, the sand origin is determined from field observations (wind directions from ventifacts), chemical and mineralogical analyses of rocks and sands. At Dyngjusandur, the sand is dominated by glass grains, a situation typical of sand plains in Iceland. Hyaloclastite ridges presently buried beneath Vatnajökull are the dominant source of the sand, and only large size plagioclase crystals (0.5 cm) in sands seem to be derived from the lava flows. Hyaloclastite ridges were crushed by glaciers and mechanically eroded sediments were washed out by melt‐water onto flood plains. The sand chemical composition is spatially homogeneous and similar to the average composition of neighboring sub‐aerial lava flows, reflecting efficient mixing of distinct sources below the glacier. The presence of sand north of Dyngjujökull can be taken as a way to explore the average chemical composition of non‐exposed volcanic material beneath the glacier. In the case of Lambahraun, prevailing winds indicate several potential sources of sand at the north of the sand sheet. Comparison of chemical and mineralogical analyses of sands and rock samples helped to refine the exact origin. In contrast with the first site, the sand is dominated by crystals and is chemically consistent with a mixture of material derived from the lava flows of Eldborgir and Skersli shield volcanoes. Analysis of the contact between the lava flows and the glacier reveals that basaltic sand grains formed as the result of recent advances of the glacier abrading the rocks. The direct interaction of glacial and fluvio‐glacial activity with basaltic plains appears to be necessary to produce a large amount of sands in a relatively short period of time (<4000 years). This site appears to be an excellent natural laboratory for further studies concerning the sand evolution and physical sorting processes in basaltic material, which have important implications for understanding aeolian processes on Mars. Copyright © 2011 John Wiley & Sons, Ltd.     

Energy approach to unsaturated cyclic strength of sand

The mechanical response to cyclic loading of saturated cohesionless soils is usually investigated by means of effective stress method considering pore water pressure changes that lead to reduced strength and stiffness. On the other hand, the behavior of partially saturated sands is different from the behavior of saturated sand deposits. The development of negative pore water pressures in particular makes it difficult to estimate the behavior of partially saturated sands. The response of partially saturated sands, however, can be examined in a physically understandable manner by investigating their energy characteristics independently of pore pressure behavior. To establish a general framework for understanding the behavior of partially saturated sand, a total of 52 resonant column and dynamic torsional shear tests were conducted under undrained conditions. The effects of factors such as the amplitude of shear strain, relative density, saturation ratio and confining pressure on the dynamic characteristics of the sand and on energy dissipation were studied. The use of the energy concept in the evaluation of partially saturated soils is shown to be a promising method for the evaluation of the cyclic behavior of partially saturated sands.     

Pore-scale imbibition experiments in dry and prewetted porous media

Two-phase imbibition behavior of immiscible fluids was studied in dry and prewetted porous media using a laser-induced fluorescence technique. Imbibition was first investigated in two-dimensional (2-D) systems under conditions comparable to those for a study of drainage [Ovdat H, Berkowitz B. Pore-scale study of drainage displacement under combined capillary and gravity effects in index-matched porous media. Water Resources Research 2006;42:W06411. doi:10.1029/2005WR004553] in the capillary-dominated regime. The effect of initial wetting saturation (IWS) was then explored in 2-D and 3-D porous media under the combined effect of gravity, capillary and viscous forces, within and outside the capillary-dominated regime. Parameters that describe maximum vertical advance, volumetric fraction, total surface area and specific surface area of the invading fluid were used to quantify the behavior. Comparison of 2-D drainage and imbibition patterns demonstrates significant qualitative differences under analogous viscosity ratio, buoyancy number, and capillary number values. However, quantitative analyses show strong pore-scale similarities between these patterns. Invasion structures in 3-D, prewetted (IWS ≈ 8% of the pore volume) porous media are ramified, with lateral branching and regions containing trapped residual fluid. These structures are qualitatively and quantitatively different from the compact, branchless structures that develop in dry (IWS = 0) porous media.