泥石流区桥墩损毁机制与控制研究

为了保证安全,泥石流多发区的道路常采用桥梁的方式通过。位于泥石流区的桥梁的桥墩在泥石流暴发时受到泥石流的冲击,常会造成不同程度的损毁。泥石流冲击力是泥石流区桥墩受到的主要荷载,目前在计算泥石流冲击力时一般将其视为均质体,但是泥石流为固液两相流体,因此得到的计算结果误差较大。为了提高精确度,本文对泥石流冲击力计算公式进行了修正,先分别运用动力学中的弹性碰撞理论和牛顿第二定律对泥石流固相和液相的冲击力进行计算,然后将其进行叠加,得到总的泥石流冲击力。泥石流区桥墩破坏的模式主要有冲击破坏、倾覆破坏和滑动破坏3种形式。为了减小泥石流对桥墩的破坏能,可在桥墩受到泥石流冲击的部位设置消能层,通过算例发现,设置消能层后冲击力减小了26.94%。     

Extended Roof snap-off for a continuous nonwetting fluid and an example case for supercritical CO2

A dominant mechanism for residual trapping of a nonwetting fluid in porous media during imbibition is snap-off or the disconnection of a continuous stream of the nonwetting fluid when it passes through pore constrictions and when a criterion based on capillary pressure imbalance is met. While quasi-static criteria for Roof snap-off have been defined for pores based on the imbalance between capillary pressure across the front/tail meniscus and local capillary pressure at the pore throat, and expressed in terms of pore body to pore throat ratio for simplification, we extended the previous quasi-static snap-off criterion by considering the local capillary pressure imbalance between the pore body and the pore throat for both circular and noncircular pores when the wetting film exists. We then used the criterion to analyze results from computational fluid dynamics (CFD) simulations of multi-phase flow with supercritical CO₂ as the nonwetting fluid and water as the wetting fluid. The extended criterion successfully described most situations we modeled. Furthermore, we compared fluid interface shape for a noncircular 3D pore predicted by the minimum surface energy (MSE) theory against 3D CFD simulations. While the fluid interface shape at the pore throat for 3D simulation was consistent with the shape predicted by MSE theory, the shape could not be successfully predicted by the MSE theory at the upstream and downstream pore body. Moreover, film flow existed for the noncircular pore at the downstream pore body.     

Dissolution of residual non-aqueous phase liquids in porous media: pore-scale mechanisms and mass transfer rates

Experiments designed to elucidate the pore-scale mechanisms of the dissolution of a residual non-aqueous phase liquid (NAPL), trapped in the form of ganglia within a porous medium, are discussed. These experiments were conducted using transparent glass micromodels with controlled pore geometry, so that the evolution of the size and shape of individual NAPL ganglia and, hence, the pore-scale mass transfer rates and mass transfer coefficients could be determined by image analysis. The micromodel design permitted reasonably accurate control of the pore water velocity, so that the mass transfer coefficients could be correlated in terms of a local (pore-scale) Peclet number. A simple mathematical model, incorporating convection and diffusion in a slit geometry was developed and used successfully to predict the observed mass transfer rates. For the case of non-wetting NAPL ganglia, water flow through the corners in the pore walls was seen to control the rate of NAPL dissolution, as recently postulated by Dillard and Blunt [Water Resour. Res. 36 (2000) 439–454]. Break-up of doublet non-wetting phase ganglia into singlet ganglia by snap-off in pore throats was also observed, confirming the interplay between capillarity and mass transfer. Additionally, the effect of wettability on dissolution mass transfer was demonstrated. Under conditions of preferential NAPL wettability, mass transfer from NAPL films covering the solid surfaces was seen to control the dissolution process. Supply of NAPL from the trapped ganglia to these films by capillary flow along pore corners was observed to result in a sequence of pore drainage events that increase the interfacial area for mass transfer. These observations provide new experimental evidence for the role of capillarity, wettability and corner flow on NAPL ganglia dissolution.     

Numerical simulation of 2D buoyant jets in ice-covered and temperature-stratified water

A two-dimensional (2D) unsteady simulation model is applied to the problem of a submerged warm water discharge into a stratified lake or reservoir with an ice cover. Numerical simulations and analyses are conducted to gain insight into large-scale convective recirculation and flow processes in a cold waterbody induced by a buoyant jet. Jet behaviors under various discharge temperatures are captured by directly modeling flow and thermal fields. Flow structures and processes are described by the simulated spatial and temporal distributions of velocity and temperature in various regions: deflection, recirculation, attachment, and impingement. Some peculiar hydrothermal and dynamic features, e.g. reversal of buoyancy due to the dilution of a warm jet by entraining cold ambient water, are identified and examined. Simulation results show that buoyancy is the most important factor controlling jet behavior and mixing processes. The inflow boundary is treated as a liquid wall from which the jet is offset. Similarity and difference in effects of boundaries perpendicular and parallel to flow, and of buoyancy on jet attachment and impingement, are discussed. Symmetric flow configuration is used to de-emphasize the Coanda effect caused by offset.     

Seasonal transition of hydrological processes in a slow‐moving landslide in a snowy region

A comprehensive understanding of seasonal hydrological dynamics is required to describe the influence of pore‐water pressure on the stability of landslides in snowy regions. This study reports on the results of continuous meteorological and hydrological observations over 2 years on a landslide body comprising Neogene sedimentary rocks in northern Japan, where a thick (3–5 m) seasonal snowpack covers the land surface. Monitoring of the volumetric water content in shallow unsaturated zones (<0.8 m depth) and pore‐water pressure in saturated bedrock at depths of 2.0 and 5.2 m revealed clear seasonality in hydrological responses to rainfall and meltwater supply. During snow‐free periods, both the shallow soil moisture and deep pore‐water pressure responded rapidly to intense rainwater infiltration. In contrast, during snowmelt, the deep pore pressure fluctuated in accordance with the daily cycle of meltwater input, without notable changes in shallow moisture conditions. During occasional foehn events that cause intense snow melting in midwinter, meltwater flows preferentially through the layered snowpack, converging to produce a localized water supply at the ground surface. This episodically triggers a significant rise in pore‐water pressure. The seasonal differences in hydrological responses were characterized by a set of newly proposed indices for the magnitude and quickness of increases in the pressure head near the sliding surface. Under snow‐covered conditions, the magnitude of the pressure increase tends to be suppressed, probably owing to a reduction in infiltration caused by a seasonal decrease in the permeability of surface soils, and effective pore‐water drainage through the highly conductive colluvial layer. Deep groundwater flow within bedrock remained in a steady upwelling state, enhanced by increasing moisture in shallow soils under snow cover, reflecting the convergence of subsurface water from surrounding hillslopes.     

Laboratory experiments in a baroclinic annulus with heating and cooling on the horizontal boundaries

Abstract

Experiments have been performed in a cylindrical annulus with horizontal temperature gradients imposed upon the horizontal boundaries and in which the vertical depth was smaller than the width of the annulus. Qualitative observations were made by the use of small, suspended, reflective flakes in the liquid (water).

Four basic regimes of flow were observed: (1) axisymmetric flow, (2) deep cellular convection, (3) boundary layer convective rolls, and (4) baroclinic waves. In some cases there was a mix of baroclinic and convective instabilities present. As a “mean” interior Richardson number was decreased from a value greater than unity to one less than zero, axisymmetric baroclinic instability of the Solberg type was never observed. Rather, the transition was from non-axisymmetric baroclinic waves, to a mix of baroclinic and convective instability, to irregular cellular convection.     

A practical model for mobile, residual, and entrapped NAPL in water-wet porous media

Flow of nonvolatile nonaqueous phase liquid (NAPL) and aqueous phases that account for mobile, entrapped, and residual NAPL in variably saturated water-wet porous media is modeled and compared against results from detailed laboratory experiments. Residual saturation formation in the vadose zone is a process that is often ignored in multifluid flow simulators, which might cause an overestimation of the volume of NAPL that reaches the ground water. Mobile NAPL is defined as being continuous in the pore space and flows under a pressure gradient or gravitational body force. Entrapped NAPL is defined as being occluded by the aqueous phase, occurring as immobile ganglia surrounded by aqueous phase in the pore space and formed when NAPL is replaced by the aqueous phase. Residual NAPL is defined as immobile, nonwater entrapped NAPL that does not drain from the pore spaces and is conceptualized as being either continuous or discontinuous. Free NAPL comprises mobile and residual NAPL. The numerical model is formulated on mass conservation equations for oil and water, transported via NAPL and aqueous phases through variably saturated porous media. To account for phase transitions, a primary variable switching scheme is implemented for the oil-mass conservation equation over three phase conditions: (1) aqueous or aqueous-gas with dissolved oil, (2) aqueous or aqueous-gas with entrapped NAPL, and (3) aqueous or aqueous gas with free NAPL. Two laboratory-scale column experiments are modeled to verify the numerical model. Comparisons between the numerical simulations and experiments demonstrate the necessity to include the residual NAPL formation process in multifluid flow simulators.     

水库蓄放水对库底岩石介质弹性波速影响的数值模拟

为了解水库水体加、卸载及渗透过程中,库底岩石介质有效应力变化对弹性波速的影响,以及弹性波速变化与孔隙压力、有效应力变化之间的定量关系,以四川省紫坪铺水库为例,在前期已建立库区地质构造与水文地质结构模型,并推导出岩石介质变形与流体渗流耦合数学模型的基础上,进一步建立了有效应力对岩石弹性波速影响的数学模型,利用有限元方法定量计算了库底岩石介质的孔隙压力、附加有效应力和弹性波速的变化。结果表明:在水库水体加、卸载及渗透过程中,地下岩石介质弹性波速增量值的变化并不明显,最大变幅仅为±0.013km/s,变化范围主要集中在库底5km以上的局部区域;弹性波速增量值与附加有效应力变化同步,与水库水位变化形态相似;不同观测点弹性波速变化量的大小与埋藏位置、深度及附加有效应力的变化幅度有关,同一观测点,弹性纵波波速变化幅度大于弹性横波波速。文中数值模拟的结果与一些学者通过环境噪声方法获得的紫坪铺水库附近区域相对波速变化的反演结果在变化形态上大体一致,但变化幅度略有差异。     

A review of visualization techniques of biocolloid transport processes at the pore scale under saturated and unsaturated conditions

Field and column studies of biocolloid transport in porous media have yielded a large body of information, used to design treatment systems, protect water supplies and assess the risk of pathogen contamination. However, the inherent “black-box” approach of these larger scales has resulted in generalizations that sometimes prove inaccurate. Over the past 10–15 years, pore scale visualization techniques have improved substantially, allowing the study of biocolloid transport in saturated and unsaturated porous media at a level that provides a very clear understanding of the processes that govern biocolloid movement. For example, it is now understood that the reduction in pathways for biocolloids as a function of their size leads to earlier breakthrough. Interception of biocolloids by the porous media used to be considered independent of fluid flow velocity, but recent work indicates that there is a relationship between them. The existence of almost stagnant pore water regions within a porous medium can lead to storage of biocolloids, but this process is strongly colloid-size dependent, since larger biocolloids are focused along the central streamlines in the flowing fluid. Interfaces, such as the air–water interface, the soil–water interface and the soil–water–air interface, play a major role in attachment and detachment, with significant implications for risk assessment and system design. Important research questions related to the pore-scale factors that control attachment and detachment are key to furthering our understanding of the transport of biocolloids in porous media.     

Influence of bottom frictional effects in sill regions upon lee wave generation and implications for internal mixing

A cross-sectional nonhydrostatic model using idealized sill topography is used to examine the influence of bottom friction upon unsteady lee wave generation and flow in the region of a sill. The implications of changes in shear and lee wave intensity in terms of local mixing are also considered. Motion is induced by a barotropic tidal flow which produces a hydraulic transition, associated with which are convective overturning cells, wave breaking, and unsteady lee waves that give rise to mixing on the lee side of the sill. Calculations show that, as bottom friction is increased, current profiles on the shallow sill crest develop a highly sheared bottom boundary layer. This enhanced current shear changes the downwelling of isotherms downstream of the sill with an associated increase in the hydraulic transition, wave breaking, and convective mixing in the upper part of the water column. Both short and longer time calculations with wide and narrow sills for a number of sill depths and buoyancy frequencies confirm that increasing bottom friction modifies the flow and unsteady lee wave distribution on the downstream side of a sill. Associated with this increase in bottom friction coefficient, there is increased mixing in the upper part of the water column with an associated decrease in the vertical temperature gradient. However, this increase in mixing and decrease in temperature gradient in the upper part of the water column is very different from the conventional change in near-bed temperature gradient produced by increased bottom mixing that occurs in shallow sea regions as the bottom drag coefficient is increased.     

湖-气热传输模型及参数敏感性研究

在一维涡扩散模型的基础上,发展了考虑冬天结冰和夏天蒸发等水体相变问题及对流混合过程的湖-气热传输模型.模型采用焓代替温度作为预报变量,这也更加方便处理水的相变问题.模型的计算结果与以色列Kinneret湖的观测资料进行对比,验证了模型的合理性,说明了加入对流混合过程是符合湖-气热传输的实际物理机制的.对风速和湖面摩擦速度作了敏感性实验,证明他们对湖面温度和蒸发潜热有着重要的影响.     

Geophysical evidence on the structure of the Taupo Volcanic Zone and its hydrothermal circulation

The Taupo Volcanic Zone (TVZ) of New Zealand is characterised by extensive volcanism and by high rates of magma production. Associated with this volcanism are numerous high-temperature (> 250 °C) geothermal systems through which the natural heat output of 4200 ± 500 MW is channelled. Outside the geothermal fields the heat flow is negligible. The average heat flux from the central 6000 km² of the TVZ, which contains most of the geothermal fields, is 700 mW/m³. This heat flux appears to be more concentrated along the eastern margin of the TVZ.Schlumberger resistivity measurements (AB/2 of 500 m and 1000 m) have identified 17 distinct geothermal fields with natural heat outputs greater than 20 MW. An additional six, low-heat-output geothermal fields also occur, and may represent formerly more active systems now in decline. Two extinct fields have also been identified. The average spacing between fields is 10–15 km. The distribution of geothermal fields does not appear to be directly associated with individual volcanic features except for the geothermal system that occurs within Lake Taupo and which occupies the vent of the 1800 yr.B.P. Taupo eruption. The positions of the geothermal fields do not appear to have varied for at least the last 200,000 years. These data are consistent with a model of large-scale convection occurring throughout the TVZ, in which the geothermal fields represent the upper portion of the rising, high-temperature, convective plumes. The majority of the recharge to the convection system is provided by the downward movement of cold meteoric water between the fields which suppresses the heat flow in these regions.Gravity measurements indicate that to a depth of about 2.5 km the upper layers of the TVZ consist of low-density pyroclastic infill. A seismic refraction interface with velocity change from 3.2 km/s to 5.5 km/s occurs at a similar depth. The cross-sectional area of the convection plumes (identified electrically) appears to increase at depths of 1–2 km, consistent with a decrease in permeability at the depth at which the velocity and density increase.The seismicity is dominated by swarm activity which accounts for about half of all earthquakes and is highly variable in both space and time. The small number of seismic events (and swarms) that have well determined depths show a cut off of seismicity at depths of 7–9 km. The depth of the transition from brittle to ductile behaviour of the rocks is identified with the transition from a regime where heat is transported by (hydrothermal) convection and pore pressures are near-hydrostatic to a regime where heat transport is dominantly conductive and pore pressures are lithostatic. Within the convective region, temperatures are moderated by the circulation of water so that the depth of the transition from convective to conductive heat transfer can be linked to the bottom of the seismogenic zone. Rocks must become ductile within about 1 km of the bottom of the overlying convective zone.Seismic refraction studies suggest that the crust beneath the TVZ is highly thinned with a seismic velocity of about 7.5 km/ s, typical of the upper mantle, occurring at depth of 15 km. Seismological studies indicate the upper mantle is highly attenuating beneath the TVZ. Conductive heat transfer between the bottom of the convective system, at about 8 km, and the base of the material with crustal velocities, at 15 km, is not able to provide all the heat that is discharged at the surface. Repeated intrusion from the mantle may provide the additional heat transport required.     

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.     

The effects of rock fragments on hydrologic and hydraulic responses along a slope

A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd.     

A semi-analytical model for computation of capillary entry pressures and fluid configurations in uniformly-wet pore spaces from 2D rock images

A novel semi-analytical model for computation of capillary entry pressures and associated fluid configurations in arbitrary, potentially non-convex, 2D pore space geometries at uniform wettability is developed. The model computes all possible centre positions of circular arcs, and physically sound criteria are implemented to determine the set of these arcs that correspond to geometrically allowed interfaces. Interfaces and pore boundary segments are connected to form closed boundaries of identified geometrical regions. These regions are classified as either oil regions, located in the wider parts of the pore space, or as water regions located in pore space constrictions. All possible region combinations are identified and evaluated for each radius value in an iterative procedure to determine the favourable entry radius and corresponding configuration based on minimisation of free energy. The model has been validated by comparison with known analytical solutions in idealised pore geometries. In cases where different analytical solutions are geometrically possible, the model generates several oil and water regions, and the valid solution is determined by the region combination that corresponds to the most favourable entry pressure, consistent with the analytical solution. Entry pressure radii and configurations are computed in strongly non-convex pore spaces extracted from an image of Bentheimer sandstone, which demonstrates that the model captures successfully well-known characteristics of capillary behaviour at different wetting conditions. The computations also demonstrate the importance of selecting the fluid configuration of minimum change in free energy. In some cases, a merged region formed by a combination of oil and water regions corresponds to the favourable entry configuration of oil, whereas in other cases, an individual oil region may correspond to the favourable oil entry configuration. It is also demonstrated that oil entry configurations may constitute merged regions for weakly water-wet conditions and individual oil regions for strongly water-wet conditions in the same pore space. The computations show that the ratio of pore area to perimeter multiplied by the cosine of the contact angle under-predicts the entry pressure radii in Bentheimer sandstone pore spaces. An alternative formula is proposed for prediction of entry radii for nonzero contact angles based on the entry radii obtained for zero contact angles.     

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: One-dimensional soil thaw with conduction and advection

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.     

Heat flow in the Basin and Range province and thermal effects of tectonic extension

In regions of tectonic extension, vertical convective transport of heat in the lithosphere is inevitable. The resulting departure of lithosphere temperature and thickness from conduction-model estimates depends upon the mechanical mode of extension and upon how rapidly extension is (and has been) taking place. Present knowledge of these processes is insufficient to provide adequate constraints on thermal models. The high and variable regional heat flow and the intense local heat discharge at volcanic centers in the Basin and Range province of the United States could be accounted for by regional and local variations in extensional strain rate without invoking anomalous conductive heat flow from the asthenosphere. Anomalous surface heat flow typical of the province could be generated by distributed extension at average rates of about 1/2 to 1%/m.y., similar to rates estimated from structural evidence. To account for higher heat flow in subregions like the Battle mountain High, these rates would be increased by a factor of about 3, and locally at active bimodal volcanic centers, by an order of magnitude more.     

Freeze Core Sampling to Validate Time‐Lapse Resistivity Monitoring of the Hyporheic Zone

A freeze core sampler was used to characterize hyporheic zone storage during a stream tracer test. The pore water from the frozen core showed tracer lingered in the hyporheic zone after the tracer had returned to background concentration in collocated well samples. These results confirmed evidence of lingering subsurface tracer seen in time‐lapse electrical resistivity tomographs. The pore water exhibited brine exclusion (ion concentrations in ice lower than source water) in a sediment matrix, despite the fast freezing time. Although freeze core sampling provided qualitative evidence of lingering tracer, it proved difficult to quantify tracer concentration because the amount of brine exclusion during freezing could not be accurately determined. Nonetheless, the additional evidence for lingering tracer supports using time‐lapse resistivity to detect regions of low fluid mobility within the hyporheic zone that can act as chemically reactive zones of importance in stream health.     

PREFERENTIAL WATER FLOW IN A FROZEN SOIL — A TWO-DOMAIN MODEL APPROACH

Earlier modelling studies have shown the difficulty of accurately simulating snowmelt infiltration into frozen soil using the hydraulic model approach. Comparison of model outputs and field measurements have inferred the occurrence of rapid flow even during periods when the soil is still partly frozen. A one-dimensional, physically based soil water and heat model (SOIL) has been complemented with a new two-domain approach option to simulate preferential flow through frozen layers. The ice is assumed to be first formed at the largest water filled pore upon freezing. Infiltrating water may be conducted rapidly through previously air-filled pores which are not occupied by ice. A minor fraction of water is slowly transferred within the liquid water domain, which is absorbed by the solid particles. A model validation with field measurements at a location in the middle-east of Sweden indicated that the two-domain approach was suitable for improving the prediction of drainage during snowmelting. In particular, the correlation between simulated and observed onset of drainage in spring was improved. The validation also showed that the effect of the high flow domain was highly sensitive to the degree of saturation in the topsoil during freezing, as well as to the hydraulic properties at the lower frost boundary regulating the upward water flow to the frozen soil and ice formation.