Air Channel Formation, Size, Spacing, and Tortuosity During Air Sparging

Characterizing mass transfer during in situ air sparging requires knowledge of the size, shape, and interfacial area of air channels. These characteristics were determined by analysis of digital images of air channels passing through submerged glass beads having particle size in the sand range. Pore-scale channeling occurred in all cases. The analysis showed that the air channels were narrower, more tortuous, more closely spaced, and moved nearly vertically through the coarser media. In the finer media, air channels had larger diameter, were spaced further apart, and passed nearly horizontally through the media. The mean diameter of the channels varied between 2.8 and 8.1 mm, and the mean spacing varied between 8.3 and 19.4 mm. Estimates of the area of the air-water interface per unit volume of soil (a₀), computed using data from the digital images and an assumed arrangement of channels, ranged from 0.02 to 0.2 mm²/mm³. Larger a₀ were obtained for coarser media and uniformly graded media. These estimates of a₀ compare well with published values for common packed-column materials and for unsaturated soils.     

Laboratory Study of Air Sparging: Air Flow Visualization

Laboratory flow visualization experiments, using glass beads as the porous medium, were conducted to study air sparging, an innovative technology for subsurface contaminant remediation. The purpose of these experiments was to observe how air flows through saturated porous media and to obtain a basic understanding of air plume formation and medium heterogeneity effects. The experiments indicate that air flow occurring in discrete, stable channels is the most probable flow behavior in medium to fine grained water saturated porous media and that medium heterogeneity plays an important role in the development of air channels. Several simulated scales of heterogeneities, from pore to field, have been studied. The results suggest that air channel formation is sensitive to the various scales of heterogeneities. Site-specific hydrogeologic settings have to be carefully reviewed before air sparging is applied to remediate sites contaminated by volatile organic compounds.     

Influence of water saturation on horizontal and vertical motion at a porous soil interface induced by incident SV wave

Recent field observations have indicated that water saturation of soils may strongly affect the vertical ground motion. A study is therefore carried out to investigate the effect of saturation on horizontal and vertical motion at an interface of porous soils with potential contributions directed to site evaluation based on field observations of both the horizontal and vertical motion. The problem described in this paper corresponds to an SV wave incident at the interface between the overlying soil and the underlying rock formation. The soils are modeled as partially water-saturated porous material with a small amount of air inclusions, while the rock are approximately regarded as ordinary one-phase solid. Theoretical formulation is developed for the computation of amplitudes of horizontal and vertical interface motion, which are expressed as functions of the degree of saturation, the angle of incidence as well as the frequency. Numerical results are given for a typical sand to illustrate the influence of saturation on the interface motion in two directions and their ratios. The present study demonstrates that the effect of water saturation may be substantial on both the horizontal and vertical motion as well as on their ratios, implying the importance of such effects in the interpretation of field observations.     

Radar determination of the spatial structure of hydraulic conductivity

Spatial variability of hydraulic conductivity exerts a predominant control on the flow of fluid through porous media. Heterogeneities influence advective pathways, hydrodynamic dispersion, and density-dependent dispersion; they are, therefore, a key concern for studies of ground water resource development, contaminant transport, and reservoir engineering. Ground-penetrating radar contributes to the remote, geophysical characterization of the macroscale variability of natural porous media. On a controlled excavation of a glacial-fluvial sand and gravel deposit in the Fanshawe Delta area (Ontario, Canada), the hydraulic conductivity field of a 45 x 3 m vertical exposure was characterized using constant-head permeameter measurements performed on undisturbed horizontal sediment cores. Ground-penetrating radar data were collected along the excavation face in the form of both reflection and common midpoint surveys. Comparison of geostatistical analyses of the permeameter measurements and the radar data suggests thatthe horizontal correlation structure of radar stack velocity can be used to directly infer the horizontal correlation structure of hydraulic conductivity. The averaging nature of the common midpoint survey is manifest in the vertical correlation structure of stack velocity, making it less useful. Radar reflection data do not exhibit a spatial structure similar to that of hydraulic conductivity possibly because reflections are a result of material property contrasts rather than the material properties themselves.     

Rapid injection of particles and gas into non-fluidized granular material, and some volcanological implications

In diatremes and other volcanic vents, steep bodies of volcaniclastic material having differing properties (particle size distribution, proportion of lithic fragments, etc.) from those of the surrounding vent-filling volcaniclastic material are often found. It has been proposed that cylindrical or cone-shaped bodies result from the passage of “debris jets” generated after phreatomagmatic explosions or other discrete subterranean bursts. To learn more about such phenomena, we model experimentally the injection of gas-particulate dispersions through other particles. Analogue materials (glass beads or sand) and a finite amount of compressed air are used in the laboratory. The gas is made available by rapidly opening a valve—therefore the injection of gas and coloured particles into a granular host is a brief (<1 s), discrete event, comparable to what occurs in nature following subterranean explosions. The injection assumes a bubble shape while expanding and propagating upwards. In reaction, the upper part of the clastic host moves upward and outward above the ‘bubble’, forming a ‘dome’. The doming effect is much more pronounced for shallow injection depths (thin hosts), with dome angles reaching more than 45°. Significant surface doming is also observed for some full-scale subterranean blasts (e.g. buried nuclear explosions), so it is not an artefact of our setup. What happens next in the experiments depends on the depth of injection and the nature of the host material. With shallow injection into a permeable host (glass beads), the compressed air in the “bubble’ is able to diffuse rapidly through the roof. Meanwhile the coloured beads sediment into the transient cavity, which is also closing laterally because of inward-directed granular flow of the host. Depending on the initial gas pressure in the reservoir, the two-phase flow can “erupt” or not; non-erupting injections produce cylindrical bodies of coloured beads whereas erupting runs produce flaring upward or conical deposits. Changing the particle size of the host glass beads does not have a large effect under the size range investigated (100–200 to 300–400 μm). Doubling the host thickness (injection depth) requires a doubling of the initial gas pressure to produce similar phenomena. Such injections—whether erupting or wholly subterranean—provide a compelling explanation for the origin and characteristics of multiple cross-cutting bodies that have been documented for diatreme and other vent deposits.     

Evidence for global circuit current flow through water droplet layers

Foggy air and clear air have appreciably different electrical conductivities. The conductivity gradient at horizontal droplet boundaries causes droplet charging, as a result of vertical current flow in the global atmospheric electrical circuit. The charging is poorly known, as both the current flow through atmospheric water droplet layers and the air conductivity are poorly characterised experimentally. Surface measurements during three days of continuous fog using new instrument techniques show that a shallow (of order 100 m deep) fog layer still permits the vertical conduction current to pass. Further, the conductivity in the fog is estimated to be approximately 20% lower than in clear air. Assuming a fog transition thickness of one metre, this implies a vertical conductivity gradient of order 10 fS m^?2 at the boundary. The actual vertical conductivity gradient at a cloud boundary would probably be greater, due to the presence of larger droplets in clouds compared to fog, and cleaner, more conductive clear air aloft.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

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.     

In Situ Vertical Circulation Column: Containment System for Small-Scale DNAPL Field Experiments

The in situ vertical circulation column (ISVCC) is a cylindrical containment system consisting of an instrumented steel cylinder used for experimental ground water studies in sandy aquifers. Vertical flow is imposed inside the ISVCC. Although vertical wells are an option, the ISVCC installed in the Borden Aquifer is instrumented with horizontal wells and monitoring ports to avoid creating vertical preferential flow paths. The cylinder was driven downward into the aquifer using a small backhoe equipped with a vibrating plate. The ISVCC penetrates the 2.3-m-thic sand aquifer and is keyed 20 cm into the underlying clay aquitard. The cylinder was installed inside a 2 m X 2 m steel sheet pile enclosure so that the enclosed segment of aquifer could be conveniently dewatered and then excavated to allow installation of the horizontal wells. The dispersivity of the column was comparable to literature values for long sand-packed laboratory columns.
Pure phase DNAPL (tetrachloroethene and 1,1,1-trichloroethane) was slowly pumped into two ports in the center of the column. Following this DNAPL injection, an aqueous solution of vitamin B₁₂ and reduced titanium was circulated through the column to promote degradation of the solvents. Processes observed in the ISVCC included DNAPL distribution, dissolution, and degradation, and geochemical evolution of the aquifer.
The ISVCC provides a convenient means for testing in situ technologies in the experimental stage or for selection of proven technologies to find the most effective at a specific site. It is inexpensive, easy to install, and maximizes control over flow distribution in a heterogeneous aquifer. Its application will be restricted where low hydraulic conductivity beds are present in the aquifer.     

A stochastic model for air injection into saturated porous media

Air injection into porous media is investigated by laboratory experiments and numerical modelling. Typical applications of air injection into a granular bed are aerated bio-filters and air sparging of aquifers. The first stage of the dynamic process consists of air injection into a fixed or a quasi-fixed water-saturated granular bed. Later stages could include stages of movable beds as well, but are not further investigated here. A series of laboratory experiments were conducted in a two-dimensional box of the size 60 cm × 38 cm × 0.55 cm consisting of glass walls and using glass beads of diameter 0.4–0.6 mm as granular material. The development of the air flow pattern was optically observed and registered using a digital video camera. The resulting transient air flow pattern can be characterized as channelled flow in a fixed porous medium with dynamic tree-like evolution behaviour. Attempts are undertaken to model the air injection process. Multiphase pore-scale modelling is currently disregarded since it is restricted to very small scales. Invasion percolation models taking into account gravity effects are usually restricted to slow processes. On the other hand a continuum-type two-phase flow modelling approach is not able to simulate the observed air flow pattern. Instead a stochastic continuum-type approach is discussed here, which incorporates pore-scale features on a subscale, relevant for the immiscible processes involved. Consequently, the physical process can be modelled in a stochastic manner only, where the single experiment represents one of many possible realizations. However, the present procedure retains realistic water and air saturation patterns and therefore produces similar finger lengths and widths as observed in the experiments. Monte Carlo type modelling leads to ensemble mean water saturation and the related variance.     

STC加筋砂水平循环剪切与竖向激振特性试验研究

随着中国汽车工业的发展,大量废旧轮胎带来的"黑色污染"问题日益显著。提出一种采用废旧轮胎柱(Scrap tire columns,STC)的加筋土结构,并初步探究其作为基础减隔振材料的可行性。通过室内水平循环剪切试验和竖向激振试验研究STC加筋砂(STCRS)的水平循环剪切和竖向激振特性。结果表明,最大剪应变为1%时STCRS的等效阻尼比未加筋前增加约10%,等效动剪切模量减小20%~25%,水平向减振性能得以提高;STCRS的竖向加速度衰减呈现出速度快、幅值大的特点,竖向减振效果较未加筋砂显著提高。STC加筋砂作为基础减振材料是可行的,且为原形废旧轮胎的资源化利用提供新思路。     

Laboratory and Numerical Investigations of Air Sparging Using MTBE as a Tracer

Air sparging experiments were conducted in a laboratory column to investigate air flow and mass transfer behavior in different types of sand at different air injection rates. Methyl tertiary butyl ether (MTBE) was applied as a tracer, and by measuring the volatilization and the mean air content during the experiments, the air flow pattern and its influence on mass transfer were assessed. The experimental results showed large differences among the sand types. In fine sand, the mean air content was high and the volatilization of MTBE was rapid with total recovery after a few hours. In coarse sand, the mean air content was low and the volatilization of MTBE was limited. The results indicate two different air flow distributions. In fine-grained materials, a uniform air distribution can be expected compared to coarse-grained materials where isolated air channels will limit the mass transfer. Afterwards, the experiments were simulated using the numerical multiphase flow code T2VOC, and the results compared to those obtained in the laboratory. The experiments with fine sand were simulated well, while for coarser sand types the volatilization was highly overestimated. The differences between model and laboratory results were mainly attributed to the nonuniformity of the air saturation and the neglection of kinetics in the mass transfer formulation.     

Particle tracking experiments in match-index-refraction porous media

A low-cost, noninvasive, three-dimensional (3D), particle tracking velocimetry system was designed and built to investigate particle movement in match-index-refraction porous media. Both a uniform load of the glass beads of the same diameter and a binary load of the glass beads of two diameters were used. The purpose of the experiments is to study the effect of the two loads on the trajectories, velocity distribution, and spreading of small physical particles. A total of 35 particles were released and tracked in the uniform load and 46 in the binary load. The 3D trajectory of each particle was recorded with two video camcorders and analyzed. It is found that the particle's velocity, trajectory, and spreading are very sensitive to its initial location and that the smaller pore size or heterogeneity in the binary load increases the particles' velocity and enhances their spreading as compared with the uniform load. The experiments also verified the previous finding that the distribution of the particle velocities are lognormal in the longitudinal direction and Gaussian in two transverse directions and that the particle spreading is much larger along the longitudinal direction than along the traverse directions.     

核电厂进水口直立式翼墙结构地基三维地震响应分析

以某核电站海域工程进水口直立翼墙为背景,运用FLAC~(3D)有限差分程序对直立式翼墙结构在地震作用时的动力响应进行模拟,结合PL-Finn液化后大变形本构模型,研究翼墙结构在动力荷载作用下的动力响应规律。从结构的位移时程、结构变形、超孔压比、液化区域等方面定量评价翼墙护岸结构的安全性。分析结果表明:砂土液化后发生流动使结构出现规律性残余变形,且随地震强度增加而变大。由地震惯性力和砂土液化共同引起的水平和竖向变形,在SL1作用下翼墙结构顶部水平残余变形0.05 m,竖向残余变形为0.07 m。在SL2作用下翼墙结构顶部水平残余变形0.26 m,竖向残余变形为0.16 m;与基底输入地震动相比,在翼墙结构顶部水平和竖直加速度放大4～5倍,且越靠近翼墙顶部处加速度呈现出明显放大效应。     

On the Use of Crushed Recycled Glass Instead of Silica Sand in Dual‐Media Filters

Crushed recycled glass was evaluated as an alternative to silica sand in dual‐media filters. Pilot scale inline filtration experiments were carried out using raw waters from three different water sources with turbidities between 6.0 and 14.0 NTU. Two physically identical filter columns were operated in parallel in the experiments. One filter consisted of 62.5 cm silica sand and 41.5 cm anthracite coal, whereas the other filter contained 62.5 cm crushed recycled glass plus 41.5 cm anthracite coal. The total bed depth was 104 cm for both filters. The properties of the media were as follows: Glass effective size = 0.77 mm, uniformity coefficient = 1.41. Sand effective size = 0.79 mm, uniformity coefficient = 1.33. Coal effective size = 1.45 mm, uniformity coefficient = 1.39. Experiments were repeated five times as follows: (i) Without the use of a coagulant, (ii–iii) with 5 and 10 mg/L of alum, and (iv–v) with 5 and 10 mg/L of ferric chloride. The filtration rate used was 11.5 m/h. Turbidity, particle counts, and head losses were measured and compared as functions of time. The following were observed: (i) Effluent turbidities and particle counts of the two filters were very close, i.e., essentially the same effluent quality was obtained when crushed glass was used instead of silica sand. (ii) In the majority of the tests, the filter with crushed glass generated both a smaller clean‐bed head loss and smaller clogging head losses than those of the filter containing sand. It is concluded that crushed glass may be a good alternative to silica sand in dual‐media filtration.     

Capillary Crack Imbibition: A Theoretical and Experimental Study Using a Hele-Shaw Cell

—?We study the filling of horizontal cracks with constant aperture driven by capillary forces. The physical model of the crack consists of a narrow gap between two flat glass plates (Hele-Shaw cell). The liquid enters the gap through a hole in the bottom plate. The flow is driven purely by the force acting on the contact lines between solid, liquid, and gas. We developed a theoretical model for this type of flow on the basis of Darcy's law; it allows for the consideration of different surface conditions.¶We run the experiment for two surface conditions: Surfaces boiled in hydrogen peroxide to remove initial contamination, and surfaces contaminated with 2-propanol after boiling in hydrogen peroxide. The flow rate depends on the gap aperture and on the interaction of the liquid with the air and the solid surfaces: The smaller the aperture, the lower the flow rate due to viscous resistance of the liquid. The flow rate is also reduced when the glass surfaces are contaminated with 2-propanol. The contact line force per unit length is approximately 60% higher on clean glass surfaces than it is on glass surfaces with the 2-propanol contamination. These experimental results are in agreement with our theoretical model and are confirmed by independent measurements of the liquid-solid interaction in capillary rise experiments under static conditions with the same Hele-Shaw cell.¶Another aspect of this study is the distribution of the liquid for the different surface conditions. The overall shape is a circular disk, as assumed in the theoretical model. However, a pronounced contact line roughness develops in case of the surfaces contaminated with 2-propanol, and air bubbles are trapped behind the contact line. A further analysis of the flow regime using the capillary number and the ratio of the viscosities of the involved fluids (water and air) reveals that the experiments take place in the transition zone between stable displacement and capillary fingering, i.e., neither viscous nor capillary fingers develop under the conditions of the experiment. The contact line roughness and the trapped air bubbles in the contaminated cell reflect local inhomogeneities of the surface wettability.     

Experimental investigation of the velocity of a sand cloud blowing over a sandy surface

The velocity of a wind‐blown sand cloud is important for studying its kinetic energy, related erosion, and control measures. PDA (particle dynamics analyser) measurement technology is used in a wind tunnel to study the probability distribution of particle velocity, variations with height of the mean velocity and particle turbulence in a sand cloud blowing over a sandy surface. The results suggest that the probability distribution of the particle velocity in a blowing sand cloud is stochastic. The probability distribution of the downwind velocity complies with a Gaussian function, while that of the vertical velocity is greatly complicated by grain impact with the bed and particle–particle collisions in the air. The probability distribution of the vertical velocity of ?ne particles (0·1–0·3 mm sands) can be expressed as a Lorentzian function while that of coarse particles (0·3–0·6 mm sands) cannot be expressed by a simple distribution function. The mean downwind velocity is generally one or two orders greater than the mean vertical velocity, but the particle turbulence in the vertical direction is at least two orders greater than that in the downwind direction. In general, the mean downwind velocity increases with height and free‐stream wind velocity, but decreases with grain size. The variation with height of the mean downwind velocity can be expressed by a power function. The particle turbulence of a blowing sand cloud in the downwind direction decreases with height. The variations with height of the mean velocity and particle turbulence in the vertical direction are very complex. It can be concluded that the velocity of a sand cloud blowing over a sandy surface is mainly in?uenced by wind velocity, grain impact with the bed and particle–particle collisions in the air. Wind velocity is the primary factor in?uencing the downwind velocity of a blowing sand cloud, while the grain impact with the bed and particle–particle collisions in the air are the primary factors responsible for the vertical velocity. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling offshore sand wave evolution

We present a two-dimensional vertical (2DV) flow and morphological numerical model describing the behaviour of offshore sand waves. The model contains the 2DV shallow water equations, with a free water surface and a general bed load formula. The water movement is coupled to the sediment transport equation by a seabed evolution equation. Using this model, we investigate the evolution of sand waves in a marine environment. As a result, we find sand wave saturation for heights of 10–30% of the average water depth on a timescale of decades. The stabilization mechanism, causing sand waves to saturate, is found to be based on the balance between the shear stress at the seabed and the principle that sediment is transported more easily downhill than uphill. The migration rate of the sand waves decreases slightly during their evolution. For a unidirectional steady flow the sand waves become asymmetrical in the horizontal direction and for a unidirectional block current asymmetrical in the vertical. A sensitivity analysis showed the slope effect of the sediment transport plays an important role herein. Furthermore, the magnitude of the resistance at the seabed and the eddy viscosity influence both the timescale and height of sand waves. The order of magnitudes found of the time and spatial scales coincide with observations made in the southern bight of the North Sea, Japan and Spain.     

垂直线电流源的三维电阻率成像

在油井注水和深井注浆的电法监测中,地表布设平面的电位测量网,而将钻井的钢套管作为线电流源来供电,据此可确定流体的运动方向和分布特征.数据处理中需解决垂直线电流源的三维电阻率反演.本文叙述了相应的反演理论和偏导数计算方法.在某油田的试验中,以生产井为电极圈的中心,每圈均匀布设18个电极,共6圈,在地表得到了108个电位数据,获得500～2 000m之间5个不同深度层的真电阻率分布图像,并对深度层为1 500～1 600m的图像进行了初步解释.结果表明,该方法可以对井下某深度层中的流体运动和分布可进行有效的监测.     

Particular solutions to the problem of horizontal flow of water and air through porous media near a wetting front

Particular solutions to the problem of horizontal flow of water and air through homogeneous porous media are derived and regularity properties of the solutions are presented. It is found that a singularity occurs in the solutions at the wetting fronts. Effects of air flow on water flow are discussed.