基于MATLAB图像处理算法的风沙两相流PIV图像的分割方法

为了深入理解沙粒起动的力学机制,需要从风洞实验中同步地获得风沙边界层中沙粒群和气流的瞬时速度信息。基于MATLAB平台的风沙两相流PIV图像分割的方法,利用沙粒群和气流示踪粒子特征粒径的显著差异,对风沙两相流PIV图像中的气固两相流进行图像分割,可分别获得沙粒群和气流示踪粒子的图像,并进一步对分离后的两相粒子按粒径进行分级分离,得到各级粒径粒子PIV图像,为进一步恢复气固两相流的流场奠定了基础。     

几种常用风沙颗粒测速方法对比

沙粒运动速度是风沙颗粒运动学和风沙流（一种特殊的气固两相流）研究中的一个重要参数,受到风沙物理学研究者们的长期关注。风沙颗粒运动速度研究遇到的关键问题是测量方面的困难,目前已经尝试应用了多种测量技术。对几种常用测试手段——光电管测速计、高速摄影、粒子动态分析仪（PDA）和粒子图像测速仪（PIV）的测量结果进行了对比分析。光电管测速计、高速摄影法、以及粒子动态分析仪在研究颗粒速度的分布方面比较有用,而粒子图像测速仪则是一种全场测速技术,是测量特定区域内平均速度场时空变化的有用工具。上述几种测试方法所得的结果既有相似之处,但也存在较大差别。根据这些测速手段所得到的沙粒速度的概率分布特征是相似的,但测量结果所反映的沙粒平均速度与风速、粒径以及其他可能影响因子间的关系存在明显差异。目前的测量技术均不能直接测量贴近层或沙床表面的沙粒速度,而这些数据在研究地表颗粒的起动机制时是相当重要的。虽然用不同的测量技术所得出的沙粒平均速度随高度的变化关系均为幂函数,但这些幂函数的估算结果差别很大。基于对比分析,我们认为,各种测速技术均有其独特的测量原理,各有其优劣,但对其测量结果的可靠性尚无定论。所以,为了充分运用不同的测速技术,需要做大量的工作对其进行验证标定,研发一种可靠的技术来测量风沙流中的沙粒速度仍然是目前风沙物理研究的重要任务。     

风沙流起动阶段沙粒输运特征

为了探究风沙流起动过程中沙粒输运特征,利用PTV测量技术在风洞中对风沙流起动过程进行了测量,分析了沙粒空间分布、沙粒平均水平速度、输沙率、沙粒数密度和输沙通量随时间的变化规律。结果表明：风沙流起动时间大约为1.5 s。起动过程中,输沙率随时间迅速增加,气流中沙粒总数目随时间的变化可表示为指数函数,沙粒数密度和输沙通量随高度的变化均可近似表示为负指数衰减函数。在t=1.0 s时刻的沙粒平均水平速度大于相同高度处以后时刻的沙粒平均水平速度,同一高度处t=1.5 s以后的沙粒数密度大于t=0.5 s、1.0 s时刻的沙粒数密度,同一高度处t=1.5 s以后的输沙通量大于t=1.0 s时刻的输沙通量。沙粒数密度随高度的衰减率一般随时间的增加而减小,并在t=1.5 s后逐渐接近稳定值。     

风沙物理学中两个焦点问题研究现状与未来研究思路刍议

风沙物理学经过数十年的发展与完善,可以说基本上自成体系。但由于学科的起步较晚,加之现象的复杂性,目前仍有很多重要的理论问题没能彻底解决,至少还未达成共识,成为学术界长期以来争论的焦点,从而影响了学科的继续发展。对“沙粒起动机制”和“风成沙纹的形成机理”两个焦点问题在进行研究现状剖析的基础上,浅析未来研究思路。提出沙粒起动主要驱动力是驻点升力的观点,以空气动力学理论为基础的沙粒起动机制研究将是未来该领域的发展趋势;关于沙纹的形成,认为沙粒粒径的不均匀性以及所导致的非平衡运动是前提和基础,风沙流运动过程中类似左右沙丘运动的蚀积转换规律是使沙纹能够维持其稳定形态的必要条件,而沙纹有规律、整齐划一的分布与排列可能与背风面气流小尺度的涡流有关。     

粗糙床面风廓线的转折特征及其物理意义

为了进一步认识粗糙床面空气动力学粗糙度的物理意义,理解空气动力学粗糙度对动量传递和跃移起动的影响机制,从风洞试验测定的粗糙床面风廓线转折特征入手,分析了粗糙床面的空气动力学性质并提出了内边界层动量传递及近壁区沙粒起动的可能机制。结果表明,细高粗糙元（方向比率在4~20之间）和孔隙粗糙元（孔隙度在0.15~0.75）风廓线呈现4个转折段,对应的湍流垂直分层为近壁区-尾涡层（Z≤0.15H~0.5H）、内边界层-尾涡层过渡层（0.15H≤Z≤0.75H）、内边界层（0.3H~0.75H≤Z≤1.2H~6H）和内边界层与外边界层过渡层等。粗矮粗糙元（粗糙元的方向比率在0.4~1.25之间）风廓线存在几个转折段,对应气流垂直分层为近壁区-尾涡层（Z≤1H~1.5H）和内边界层（1H~1.5H≤Z≤7H~35H）等。细高粗糙元和孔隙粗糙元覆盖21组床面（侧影盖度在0.007~0.50,粗糙元的高度在10~100 mm）的内边界层内空气动力学粗糙度在0.07~30.74 mm之间,比内边界层以上或以下过渡层的空气动力学粗糙度高几个数量级到数倍;内边界层摩阻风速的0.50~1.66 m\5s-1之间,是内边界层以下过渡层的1.5~10倍、内边界层以上过渡层的摩阻风速的1.1~2.8倍。内边界层的空气动力学粗糙度和摩阻风速分别代表了粗糙床面对气流阻力特征和湍流切应力。内边界层以下过渡层湍流切应力与粗糙元之间光滑地表所受切应力关系不大,而近壁层切应力与光滑地表所受切应力直接相关。近壁层猝发过程上抛运动和内边界层湍流猝发过程下扫运动耦合关系是内边界层的一部分动量传递到近壁区并导致沙粒起动的可能机制。     

PIV技术及其在风沙边界层研究中的应用

为了考察粒子图像测速度技术在风沙环境风洞中的测量精度及在风沙边界层研究中的应用潜力,通过筛选适当的示踪颗粒,借助PIV测量系统重新测量了风沙环境风洞中的风廓线,并获得了风沙边界层内跃移沙粒的速度和浓度分布规律。实验结果表明：PIV测得的风速廓线与标准风速廓线仪所测结果相当吻合(R2≥0.99);沙粒跃移的平均水平速度和相对浓度（灰度）沿距离沙面高度分别呈幂函数和负指数分布,沙粒速度随高度和自由风速的增加而增大,相对浓度随高度的增加而快速衰减,风速越小衰减越快,风速越大衰减越慢,这一结果与前人的结论一致。PIV系统为将来能够进行更加精确的风沙运动微观机理研究提供了技术保证。     

风沙床面颗粒起动临界受力平衡模型的对比分析

风沙颗粒起动的受力平衡模型是风沙起动机制研究的重要内容。根据已有的研究成果,论述风沙颗粒起动的临界条件,介绍几种代表性模型。认为,影响自然状况下颗粒起动过程的因素是复杂的,尽管对主要影响因素的作用还未获得深入的认识,但是已提出的风沙颗粒起动的受力平衡模型对于认识风沙起动过程和机理具有重要的指导和启发作用,然而这些模型往往包含了较多的简化、假设以及有待于通过实验测定的系数,所以,不能够在实践中直接应用。而且,这些模型都是在特定的条件下得出来的,其应用范围具有局限性。有鉴于此,关于风沙颗粒起动的机制研究仍有较多的创新空间。     

稳态风沙流中瞬态输沙特征

风沙流中沙粒运动在来流风速不变时也会表现出非稳态特征。在风洞内利用粒子图像测速系统（PIV）测量了风沙运动的时间序列,并基于PIV测量技术提出风沙流中沙粒平均直径、数密度、平均水平速度和输沙通量等参数在某一时刻的计算方法,其中输沙通量的计算考虑沙粒大小垂向分布的影响。结果表明：来流风速不变时,沙粒平均直径、数密度、平均水平速度和输沙通量随时间具有明显的波动性;沙粒平均直径和平均水平速度的标准偏差一般随高度增加而增加,沙粒数密度和输沙通量标准偏差随高度增大而减小;这些参数的相对标准偏差均随高度增加而增大。     

粒子图像测速(PIV)技术在风沙环境风洞中的应用

粒子图像测速技术是一种全新的非接触、瞬时、全场流速测量方法,广泛应用于流体力学实验中。笔者将其引入风沙运动研究是一个尝试,旨在为风沙物理和防沙工程研究提供新的测量手段。通过分析一系列的风洞实验结果,我们认为,在风沙流中颗粒速度分布、浓度分布、风沙地貌形成演化、防沙工程原理及应用方面,粒子图像测速技术具有一定优势,可以克服传统测量手段的某些不足。     

海岸沙丘沙运动特征若干问题的研究——以闽江口南岸为例

在海岸带,沙粒粒径和水分含量是决定沙粒起动风速的最主要因素。闽江口南岸海滩沙的平均粒径为2.0Φ,水分含量多在0.3%~1.0%之间。运用Johnson的经验公式,可计算出海滩沙的起动风速为6~7m/s。这一结果与实际观测值十分吻合。由于海岸环境因素的影响,使得同等条件下海岸带沙粒的起动风速明显大于内陆沙漠沙,而风沙流的输沙率明显小于内陆沙漠沙。海岸带是风速剧变带,在引用陆地气象台站的风速记录分析海岸风况和研究风沙运动规律时,必须对风速进行订正。     

The dynamic effects of moisture on the entrainment and transport of sand by wind

There is an increasing awareness of the influence of surface moisture on aeolian entrainment and transport of sediment. Existing wind tunnel studies have shown the impact of a limited range of moisture contents on entrainment thresholds but similar investigations are lacking in the field. The research reported here investigated the influence of changes in surface moisture content on sand entrainment and transport on a meso-tidal beach in Anglesey, North Wales.High frequency (1 Hz) wind velocities measured with hot-wire anemometers were combined with grain impact data from a Sensit monitor and mass flux measurements from a standard sand trap. Surface and near-surface moisture contents were assessed gravimetrically from surface sand scrapes and also directly by using a ThetaProbe. Critical threshold values for entrainment were specified using a modified form of the time fraction equivalence method (Stout, J.E., Zobeck, T.M., 1996a. Establishing the threshold condition for soil movement in wind-eroding fields. Proceedings of the International Conference on Air Pollution from Agricultural Operations. MWPS C-3, Kansas City, 7–9 February 1996, pp. 65–71).Results indicate a time-dependent change in dominant control of the sand transport system from moisture to wind speed, dependent upon the moisture content of the surface sediment. This interchange between controlling parameters on both entrainment and transport was very sensitive to prevailing moisture conditions and took place over a period of minutes to hours. Under conditions experienced in the experiments presented here, the critical moisture threshold for sediment entrainment was determined to be between 4% and 6%, higher than the 1–4% specified in previous wind tunnel experiments. Furthermore, a moisture content of nearly 2% (where moisture was adhered to transported sediment) appeared to have little or no impact on the rate of sand flux.     

塔克拉玛干沙漠腹地湍流能量耗散率和结构函数参数特征

利用塔中气象站垂直探测系统的湍流观测资料,计算并分析了塔克拉玛干沙漠腹地不同稳定层结条件的湍涡特征长度尺度、能量耗散率及湍流结构函数参数特征。结果表明：湍涡特征长度尺度在弱不稳定或近中性条件时最大,随不稳定程度的增强有明显减小的趋势,随着稳定程度的增加有先迅速减小后又缓慢增加的趋势,且热量特征长度尺度比动量特征长度尺度总体要大一个数量级;无因次湍流热量耗散率的拟合函数形式与其他试验一致,但经验系数略有差异;利用Kaimal公式和湍流谱方法间接计算的稳定条件湍流能量耗散率及湍流结构函数参数具有较好的一致性。     

Numerical simulation of aeolian dusty sand transport in a marginal desert region at the early entrainment stage

Aeolian dusty sand transport in the marginal region of a desert is described numerically from first suspension motion at the early entrainment stage to the unsteady state within a moderate range (1000 m long × 500 m high). A two-dimensional model is built for wind-blown dust flow, and the calculations are carried out using Fluent software. The simulation results describe an integrated picture of aeolian dusty sand transport including uplift, suspension, diffusion, deposition and its space–time concentration. According to the features of a sand–dust storm near the surface, a volume concentration expression of ejection grains is developed as a boundary condition in the simulation approach. The model is verified by comparing results with both experimental data from a wind tunnel and an analytic solution. Uniform dust sizes and R–R distributions are used in the simulation. Gas–solid two-phase flow patterns with these grains are obtained in the downwind space, including the turbulence intensity, gas phase stream functions and solid volume concentration distributions. The influence of wind velocities and grain sizes is analyzed. From the simulation results, spatial distributions of dust volume concentration in the early entrainment stage are described clearly. Different from coarse sands, there is a clear band of uniformly saturated dust concentration in the region directly above the surface.     

湍流度随高度的变化及其对城市宏观地形的依赖

分别构建广州主建成区垂直比例尺为1﹕2 000、1﹕1 000和1﹕500的3个建筑物模型,利用大型边界层风洞,在西北和东南两风向下,基于中性流模拟分析了复杂城市地形下湍流度随高度的变化及其对宏观地形的依赖。结果表明：风廓线指数α与不同高度的湍流度之间的关系密切,利用现有模型,根据4类粗糙度边界层和不同垂直比例尺,可确定相应的湍流度随高度变化模型的主要系数,预测精度高。城市地形下最大湍流度面发育在0~0.2 h之间狭窄的范围内。用湍流度形态指数β来表征湍流度随高度的变化,无论城市屋脊还是平坦地形,随着风程区的延伸,廓线的指数α升高,湍流度形态指数β降低。表明同一高度湍流度值具有由迎风区、丘顶区向背风区增高,沿风程逐渐增大的规律,对地形部位和风程的依赖性强,与来流翻越简单地形时的特征一致。     

风沙流中风速纵向脉动的实验研究

在4种风速条件下对风沙流中风速纵向脉动特性进行了风洞实验研究。研究发现，风速纵向脉动很大．不容忽视；脉动绝对大小随进口风速提高有增强趋势，随着高度的增加则变化不大；垂线脉动强度平均值随主流速度基本成正向线性关系；以脉动强度与当地风速之比表示的相对脉动强度随高度增加而降低，并且不同风速条件下相对脉动强度在床面附近相差较大，远离床面则趋向一致。另外，采用幂函数形式可较好地描述风沙流中平均风速的垂线分布。     

METEOROLOGICAL CONDITIONS DURING SLUSH-FLOW RELEASE AND THEIR GEOMORPHOLOGICAL IMPACT IN NORTHWESTERN ICELAND: A CASE STUDY FROM THE BÍLDUDALUR VALLEY

ABSTRACT. This paper examines triggering factors and geomorphic significance of slush flows in the Bíldudalur valley, northwestern Iceland. The area is prone to release slush flows from two confined gullies, and at least ten flows have been reported since the beginning of the twentieth century. Despite their short path (600 m) and their moderate magnitude (from 6000 to 8000 m³), slush flows in the Bíldudalur valley represent a serious threat for the local community that is situated within the runout and deposition zones. With the help of meteorological data, the release of known slush flows is examined, highlighting the role of heavy rainfall and rapid snow-melt during winter cyclonic activity. The geomorphological impact of slush flows is assessed through the characteristics of the landforms produced during the 1997 and 1998 slush-flow events. It appears that the most obvious characteristics of slush flows in the Bíldudalur valley are the entrainment and deposition of debris, spatially differentiated. Chaotic sedimentation occurs chiefly in the middle part of a clearly concave cone, even if the flows continue beyond the cone.     

塔克拉玛干沙漠湍流尺度

利用塔克拉玛干沙漠腹地塔中地区近地层32 m高度涡动相关系统2009年1、4、7、10月的湍流观测数据,对流动沙漠地区不同稳定度下的湍流特征进行了分析。结果表明:(1)流动沙漠区在不同季节、不同稳定度条件下的近地层湍流欧拉时间自相关系数R_u(t)、R_v(t)、R_w(t)衰减存在明显的时间差异。春季R_w(t)的衰减最快,夏季3个方向的欧拉时间自相关系数衰减时间相近;(2)欧拉时间自相关系数主要表现出强衰减、逐渐衰减和缓慢衰减3种类型,沙尘天气的欧拉时间自相关系数衰减总体上要强于晴天;(3)湍流的时间积分尺度L_t(z)最小,L_t(x)和L_t(y)大小相近,在沙尘暴发生时,水平风速占主导,而垂直气流很难发展;(4)湍流的空间尺度在不同季节和不同稳定度条件下明显不同:L(x)平均值6.02 m,L(y)平均值6.90 m,L(z)平均值0.73 m;湍流水平空间尺度大多数情况下呈扁平的椭球体,且水平方向呈明显的椭圆形状。     

A New Approach for the Health Assessment of River Systems Based on Interconnected Water System Networks

Interconnected river system networks is a national water conservancy strategy in China and focus of research. Here we discuss the classification system, material and energy exchange between rivers and lakes, various dynamic flows and ecological functions of river-lake interconnected relationships. We then propose a novel method for the health assessment of river systems based on interconnected water system networks. In a healthy river system there is “material and energy exchange” and it is the first and foremost relationship of material and energy exchange between rivers and lakes. There are unobstructed various “flows” between rivers and lakes including material flows (water, dissolved substances, sediments, organisms and contaminants), energy flows (water levels, flow and flow velocity), information flows (information generated with water flows, organisms and human activities) and value flows (shipping, power generation, drinking and irrigation). Under the influences of nature and human activity, various flows are connected by river-lake interconnection to carry material and energy exchange between rivers and lakes to achieve river-lake interactions. The material and energy exchange between rivers and lakes become one of the approaches and the direct driving forces of changes in river-lake interconnected relationships. The benignant changes in river-lake interconnected relationship tend to be in relatively steady state and in ideal dynamic balance.     

大气水汽同位素组成的短期变异特征

Stable isotopes of atmospheric water vapor reveal rich information on water movement and phase changes in the atmosphere. Here we presented two nearly continuous time-series of δD and δ18O of atmospheric water vapor (δv) measured at hourly intervals in surface air in Beijing and above a winter wheat canopy in Shijiazhuang using in-situ meas-urement technique. During the precipitation events, the δv values in both Beijing and Shiji-azhuang were in the state of equilibrium with precipitation water, revealing the influence of precipitation processes. However, the δv departures from the equilibrium state were positively correlated with local relative humidity. Note that the δv tended to enrich in Beijing, but deplete in Shijiazhuang during the precipitation events, which mainly resulted from the influence of transpiration processes that enriched the δv in Shijiazhuang. On seasonal time-scale, the δv values were log-linear functions of water vapor mixing ratios in both Beijing and Shijiazhuang. The water vapor mixing ratio was an excellent predictor of the δv by the Rayleigh distillation mechanisms, indicating that air mass advection could also play an important role in deter-mining the δv. On a diurnal time-scale, the δv reached the minimum in the early afternoon hours in Beijing which was closely related to the atmospheric processes of boundary layer entrainment. During the peak of growing season of winter wheat, however, the δv reached the minimum in the early morning, and increased gradually through the daytime, and reached the maximum in the late afternoon, which was responsible by the interaction between boundary layer entrainment and the local atmospheric processes, such as transpiration and dew for-mation. This study has the implications for the important role of vegetation in determining the surface δv and highlights the need to conduct δv measurement on short-term (e.g. diurnal) time scales.