林业科学  2018, Vol. 54 Issue (7): 73-83   PDF    
DOI: 10.11707/j.1001-7488.20180708
0

文章信息

孙兴林, 张宇清, 张举涛, 秦树高, 周金星
Sun Xinglin, Zhang Yuqing, Zhang Jutao, Qin Shugao, Zhou Jinxing
青藏铁路路基对风沙运动规律影响的数值模拟
Numerical Simulation on the Influence of Subgrade of Qinghai-Tibet Railway on Wind-Sand Movement
林业科学, 2018, 54(7): 73-83.
Scientia Silvae Sinicae, 2018, 54(7): 73-83.
DOI: 10.11707/j.1001-7488.20180708

文章历史

收稿日期:2017-05-15
修回日期:2017-10-14

作者相关文章

孙兴林
张宇清
张举涛
秦树高
周金星

青藏铁路路基对风沙运动规律影响的数值模拟
孙兴林1, 张宇清1,2, 张举涛1, 秦树高1,2, 周金星1,2     
1. 北京林业大学水土保持学院 北京 100083;
2. 水土保持国家林业局重点实验室 北京 100083
摘要:【目的】通过Fluent软件进行模拟计算,研究青藏铁路原有路基及新提出的一种输沙型路基的流场分布情况以及路基积沙特征,并讨论了轨枕式挡沙墙与PE阻沙网对路基风场及积沙特征的影响,探索路基对风沙运动规律的影响,为青藏铁路防沙治沙提供理论基础。【方法】采用数值模拟的方法,选择3种路基(普通路基、通风路基、输沙路基)以及2种沙障(轨枕式挡沙墙、PE阻沙网)的组合形式,通过模拟流场及积沙分布特征,评价沙障的防沙效果并分析路基受风蚀影响的程度。【结果】3种路基的风速流场具有相同的功能分区,且风速对流场的影响较小;在路基迎风侧,普通路基对风速的降低作用最强,极易引起沙粒在路基坡脚和边坡堆积,路基顶部风速明显增加,大风时可能会影响行车安全,路基背风侧存在回流现象,部分沙粒在重力作用下会堆积在背风侧边坡处;普通路基迎风侧边坡积沙情况严重,沙埋铁轨现象频发,布设有通风管的路基会受到积沙影响,导致管道堵塞,影响其冻土保护功能,输沙路基利用狭管效应,使风沙流加速通过,防止迎风侧边坡积沙和风沙上道,且一旦输沙管道发生积沙,可通过自然风力清除,若必须人工清理也便于操作,可节省大量人力物力;当路基上风向布设沙障时,沙障对风速有极强的减弱作用,路基和沙障之间存在一个减速区,该区域速度呈中间大两端小分布规律,风速无法达到起沙风速,沙粒沉积在沙障前后,路基边坡积沙明显减少。【结论】3种路基在无防护措施情况下,均有不同程度的积沙,输沙路基积沙量最小,且便于清理;轨枕式挡沙墙积沙量迎风侧大于背风侧,PE阻沙网积沙量背风侧大于迎风侧,并且沙障间距影响着阻沙效果,适当调整沙障间距,使此区间内的风速不足以起动沙粒,提高沙障阻沙效果;模拟结果表明,输沙路基相对于普通路基和通风路基更适用于风沙地区,但目前的结果还停留在理论阶段,需要在高原环境中进一步进行试验验证。
关键词:青藏铁路    风沙流    Fluent软件    数值模拟    路基    
Numerical Simulation on the Influence of Subgrade of Qinghai-Tibet Railway on Wind-Sand Movement
Sun Xinglin1, Zhang Yuqing1,2 , Zhang Jutao1, Qin Shugao1,2, Zhou Jinxing1,2    
1. School of Soil and Water Conservation, Beijing Forestry University Beijing 100083;
2. Key Laboratory of State Forestry Administration on Soil and Water Conservation, Beijing Forestry University Beijing 100083
Abstract: 【Objective】Qinghai Tibet railway is an important link in the development of the western region of China. As a result of the special natural environment, the frequent occurrence of wind sand disaster has become one of the most important factors affecting the safe operation of the railway. In this study, we used the Fluent software to simulate the different types of Qinghai-Tibet railway subgrade with sand barrier layout and analyzed its wind speed flow field and sediment characteristics. Our aim is to provide a theoretical basis for reasonable layout of sand prevention measures to reduce the accumulation of sand in the subgrade and is to ensure the safe operation of the railway.【Method】With the numerical simulation method, three kinds of subgrade(Ordinary subgrade, Ventilated subgrade, Sand-transmitting subgrade) and two sand barriers(Sleeper typed retaining wall, Sand blocking PE nets) combinations were selected to simulate the flow field and sand distribution, for evaluating the sand prevention effect of sand barrier and analyzing the affecting extent of the subgrade by wind erosion.【Result】Simulation result show that the wind speed field of three types of subgrade all had the same functional partition, and the wind speed had little effect on the flow field. On the windward side, the effect of ordinary subgrade on reducing wind speed was the strongest. Sands were easy to accumulate on the slope and the foot of the slope. The wind speed was significantly increased at the top of subgrade which may affect traffic safety. There was a reflow phenomenon at leeward, leading to the accumulation of sands on the leeward slope under gravity. The sand accumulation on the windward slope of common subgrade was serious and buried rail frequent. The ventilated subgrade would be affected by sand sediment, limiting the permafrost protection. The sand-transmitting subgrade prevented sand accumulation on the windward slope and rail. Once the occurrence of sand deposition on the transport pipeline, could be removed by natural wind. It was also easy to clean up and save a lot of manpower and resources. When the sand-break wall laying before the subgrade, the sand-break wall on the wind speed had a strong weakening effect. There was a deceleration zone between the subgrade and the sand-break wall, the wind speed could not reach the sand-driving wind speed. Sands deposited before and after sand-break wall, the amount of sediment accumulation on subgrade slope was greatly reduced.【Conclusion】Three kinds of subgrade in the absence of protective measures, there are different degrees of sand accumulation. The amount of sediment of sand-transmitting subgrade is minimal, and easy to clean up. The sand accumulation of Sleeper typed retaining wall on the wind side is greater than the leeward side. The sand accumulation of Sand blocking PE nets on the leeward side is greater than the wind. The spacing of sand barrier affects the effect of sand blocking. Appropriate adjusting sand barrier spacing to improve sand blocking effect. According to the simulation, the sand-transmitting subgrade is more suitable for the sandy area. However, the result are still in the theoretical stage and need to be further tested in Plateau environment.
Key words: Qinghai-Tibet railway    wind-blown sand flow    Fluent software    numerical simulation    subgrade    

青藏铁路东起西宁, 西至拉萨, 全长1 956 km, 是世界海拔最高的高原铁路, 途经河谷、退化草原、戈壁等多种地貌类型(谢胜波等, 2014), 穿越海拔4 000 m以上地段达960 km, 最高点海拔5 072 m。青藏高原地理位置特殊, 自然灾害频发, 其中风沙灾害是影响铁路安全运行的最重要因素之一(邱青云等, 2004)。青藏铁路多处经过风沙地段, 沙害面积约880 km2, 风沙段长度总计442.6 km, 占线路总长的22.5%(Zhang et al., 2010; 张克存等, 2011; 薛智德等, 2010), 其中轻度沙害路段282 km, 中度沙害路段105 km, 严重沙害路段53 km(于云江等, 2001; 张克存等, 2010)。目前, 在沙害严重路段已设立系统的防沙措施并取得良好成效, 但随着时间的推移, 部分防沙措施已逐渐失去作用, 沙害问题仍十分严峻(牛清河等, 2009)。

关于风沙流的研究手段主要有2种, 一种是传统的野外观测及风洞试验, 另一种是近年来兴起的数值模拟。与传统手段相比, 数值模拟不仅可以减少野外的工作量, 而且节约成本, 又可以快速得到结果, 因此越来越多地应用于风沙运动的研究。李驰等(2010)利用Fluent软件对风沙环境下铁路路基的风蚀破坏规律进行了数值模拟研究, 分析了路基高度和边坡比对流场的影响, 结果显示风速减弱区和恢复区主要受路基高度影响。石龙等(2014)基于Fluent软件, 对兰新铁路路堤周围风沙两相流运动特性进行了模拟研究, 揭示了迎风坡积沙量大于背风坡的形成机制。李晓军等(2016)对兰新铁路既有挡风墙周围风沙两相流运动特性进行了数值模拟, 阐述了挡风墙背风侧风速廓线变化规律, 提出了挡风墙建设高度的合理建议。然而, 对于青藏铁路路基流场及积沙的研究仍然十分薄弱, 路基周围流场及积沙分布规律仍然不很清楚。本文通过Fluent软件进行模拟计算, 研究青藏铁路原有路基及新提出的一种输沙型路基的风场分布情况以及路基积沙特征, 并讨论了轨枕式挡沙墙与PE (Polyethylene)阻沙网对路基风场及积沙特征的影响, 探索路基对风沙运动规律的影响, 为青藏铁路防沙治沙提供理论基础。

1 模拟方法

数值模拟方法包括前处理、计算、后处理3个步骤。前处理使用Gambit 2.4建立几何模型及网格划分, 计算过程采用Fluent 16.0, 后处理使用Tecplot 360。具体步骤见图 1

图 1 风沙运动数值模拟方法 Figure 1 Simulation methods of wind-sand movement numerical simulation
1.1 几何模型

风沙流是风成沙在多种作用力下的三维运动现象, 需要考虑多个因素在3个空间方向上的变化, 十分复杂, 由于风沙流在运动过程中主要受到水平和垂直方向上力的影响, 因此采用二维模型进行模拟。

本研究对3种不同的路基进行模拟。第1种是普通路基, 根据《铁路路基设计规范》(TB10001—2005), 按比例对普通长直型铁路路基进行建模, 具体规格见图 2a。第2种是通风路基, 根据刘志强等(2006)的研究, 通风路基规格如图 2b, 其中区域Ⅰ为道碴层, 区域Ⅱ为碎石砂砾。通风管在区域Ⅱ中, 管中心线距地面1 m, 管直径0.4 m。第3种是输沙型路基, 根据现有普通长直路基创新而来, 主要适用于主风向单一并与铁路线夹角较大或垂直的严重风沙路段。具体规格如图 2c所示, 其中:进口段的底面斜坡坡度为1/8, 集风段底面斜坡坡度为1/15, 输沙段底面为水平面, 出口段的底面斜坡坡度为1/10, 集沙池在沿路基延伸方向上的底面斜坡坡度依地势分段设置坡降, 其长度约为输沙洞长度的1/4~2/5。

图 2 3种路基截面示意 Figure 2 Subgrade sectional views a:普通路基 Ordinary subgrade; b:通风路基 Ventilated subgrade; c:输沙路基 Sand-transmitting subgrade.
1.2 边界条件设置

计算域左侧边界是风沙流入口, 定义为速度进口边界条件(VELOCITY_INLET), 右侧是风沙流出口边界, 定义为完全发展出流边界条件(OUT_FLOW), 计算区域顶部采用对称边界条件(SYMMETRY), 路基模型及地面定义为无滑移壁面(WALL), 介质类型为FLUID。

本模拟中入口速度通过udf函数自定义, 垂直风速廓线表示为:

$ V = {V_1}{\left({z/{z_1}} \right)^\alpha } 。$ (1)

式中:z代表任意高度; z1取值为2, 代表 2 m高度; V1代表z1高度处的风速; V代表z高度处的风速; α代表粗糙度指数, 取值0.2。

1.3 参数设置

本文对风沙流进行非稳态模拟, 采用欧拉多相流模型, 由于风沙流经过路基时会产生涡旋流, 假定湍流完全发展, 使用标准k-ε湍流模型, 经计算得湍流强度I=0.02。

空气密度ρg=1.225 kg·m-3, 沙粒密度ρs= 2 650 kg·m-3, 粒径分布主要在0.05~0.25 mm之间, 取0.1 mm, 沙粒相体积分数取0.02。空间离散化都使用一阶迎风格式, 时间步长取0.001 s, 流场求解采用SIMPLE算法。

1.4 基本方程

1) 质量守恒方程  对于不可压缩的流体流动, 密度为常数, 表达式如下:

$ \frac{\partial }{{\partial t}}\left({{\alpha _q}{\rho _q}} \right) + \nabla \cdot \left({{\alpha _q}{\rho _q}{{\vec v}_q}} \right) = 0 。$ (2)

式中:αq表示第q相体积分数, 风沙流是气固2相流, 因此气体相与沙粒相体积分数之和为1, 即αg+αs=1;ρq为第q相的密度; ${{{\vec v}_q}}$为第q相的矢量速度。

2) 动量守恒方程  气体相:

$ \begin{array}{l} \frac{\partial }{{\partial t}}\left({{\alpha _{\rm{g}}}{\rho _{\rm{g}}}{{\vec v}_{\rm{g}}}} \right) + \nabla \cdot \left({{\alpha _{\rm{g}}}{\rho _{\rm{g}}}{{\vec v}_{\rm{g}}}{{\vec v}_{\rm{g}}}} \right) = \\ \;\;\; - {\alpha _{\rm{g}}}\nabla p + \nabla \cdot {\overline{\overline \tau } _{\rm{g}}} + {\alpha _{\rm{g}}}{\rho _{\rm{g}}}\vec g + {f_{{\rm{sg}}}} ;\end{array} $ (3)

沙粒相:

$ \begin{array}{l} \frac{\partial }{{\partial t}}\left({{\alpha _{\rm{s}}}{\rho _{\rm{s}}}{{\vec v}_{\rm{s}}}} \right) + \nabla \cdot \left({{\alpha _{\rm{s}}}{\rho _{\rm{s}}}{{\vec v}_{\rm{s}}}{{\vec v}_{\rm{s}}}} \right) = \\ - {\alpha _{\rm{s}}}\nabla p + \nabla {p_{\rm{s}}} + \nabla \cdot {\overline{\overline \tau } _{\rm{s}}} + {\alpha _{\rm{s}}}{\rho _{\rm{s}}}\vec g + {f_{{\rm{sg}}}} 。\end{array} $ (4)

式中: ${\vec g}$为重力加速度; p为2相共享的相同压力; ps是沙粒相压力; ${\overline{\overline \tau } _{\rm{g}}}、{\overline{\overline \tau } _{\rm{s}}}$分别为气体相和沙粒相的应力-应变张量; fsg为气体相和固体相间的相互作用力。

本文使用标准k-ε湍流模型, 由于沙粒相体积分数只有0.02, 因此气体相的湍流运动在风沙流中占主导作用, 假设湍流完全发展, 且不考虑源相的影响。

3) 湍流动能k方程:

$ \rho \frac{{{\rm{D}}k}}{{{\rm{D}}t}} = \frac{\partial }{{\partial {x_i}}}\left[ {\left({{\mu _l} + \frac{{{\mu _t}}}{{{\sigma _k}}}} \right)\frac{{{{\rm{δ }}_\mathit{k}}}}{{{\rm{ δ }}{x_i}}}} \right] + {G_k} + {G_b} - \rho \varepsilon ;$ (5)

4) 湍流耗散率ε方程:

$ \begin{array}{l} \rho \frac{{{\rm{D}}\varepsilon }}{{{\rm{D}}\varepsilon }} = \frac{\partial }{{\partial {x_i}}}\left[ {\left({{\mu _l} + \frac{{{\mu _t}}}{{{\sigma _\varepsilon }}}} \right)\frac{{{{\rm{ δ }}_\varepsilon }}}{{{\rm{ δ }}{x_i}}}} \right] + \\ {C_{1\varepsilon }}\frac{\varepsilon }{k}\left({{G_k} + {G_{3\varepsilon }}{G_b}} \right) - {C_{2\varepsilon }}\rho \frac{{{\varepsilon ^2}}}{k} ;\end{array} $ (6)

式中:μ=μt+μl, μt= $\rho {C_u}\frac{{{k^2}}}{\varepsilon }$。其中:Gk是由层流速度梯度而产生的湍流动能; Gb是由浮力产生的湍流动能; μ是有效黏性系数; μl是层流黏性系数; μt是湍流黏性系数; C1εC2εG3εσk, σε为经验常数, Cu为湍流常数。

2 结果与分析 2.1 路基平面流场速度分布特征

根据观测资料, 青藏铁路部分路段持续时间较长的输沙风速为15 m·s-1左右, 因此在模拟时分别测试了高度2 m处风速为10、15、20 m·s-1时的流场情况, 沙粒相体积分数取0.02, 3种路基的流场速度分布见图 3

图 3 不同风速下3种路基流场分布特征 Figure 3 Distribution characteristics of three subgrade flows for different wind speeds a-c:普通路基风速10、15、20 m·s-1时的流场 The flow field of ordinary subgrade at 10, 15 and 20 m·s-1; d-f:通风路基风速10、15、20 m·s-1时的流场 The flow field of ventilated subgrade at 10, 15 and 20 m·s-1; g-i:输沙路基风速10、15、20 m·s-1时的流场 The flow field of sand-transmitting subgrade at 10, 15 and 20 m·s-1. x:水平位置 Horizontal position; y:高度 Hight.所有风向与图a相同, 下同。All wind directions are the same as Figure a, the same as below.

图 3可知, 不同风速的流场没有明显差异。3种路基的流场存在相似分区, 风沙流在平坦戈壁上流动, 途经路基阻碍时, 在迎风侧坡脚处速度降低形成减速区(A); 路基迎风坡在垂直方向上对气流有抬升作用, 速度逐渐增大, 形成抬升区(B); 路基上方由于气流聚集, 风速增大, 形成高速区(C); 在路基背风侧有明显的涡旋流, 风速降低, 部分区域风速变为负值, 该区域风向与来流风向相反, 形成回流区(D); 风沙流远离路基之后, 气流开始扩散, 速度缓慢上升, 最终恢复到初始状态, 为恢复区(E); 此外, 针对通风路基和输沙路基, 由于路基内部存在管道, 气流从开阔地带流入狭窄的管道内, 空气不能大量堆积, 因此加速流过管道, 形成加速区(F)。

2.2 路基边坡特征位置流场速度变化情况

从上述分析中可知, 风速对流场变化的影响较小, 本部分着重分析风沙流速度为15 m·s-1时, 路基边坡特征位置处的速度变化情况。在野外调查中, 笔者发现积沙主要位于路基迎风侧坡脚、铁轨内和背风侧坡脚, 因此分析这3个断面0~10 m的风速廓线(卢玉华等, 2012)。由图 4可见, 在迎风侧坡脚处, 3种路基的风速廓线变化趋势一致, 输沙路基的风速随高度增幅最大, 说明在迎风坡输沙路基对气流阻碍作用最弱, 沙粒最不易遇阻沉积。路基顶部选取铁轨中轴位置, 在0~2 m范围内, 气流穿过通风路基和输沙路基内部的管道, 风速激增, 管道中部风速高于近管壁处风速, 并且输沙路基的管道内部风速更大; 2~5 m处是路基的实体部分, 风速值为0。路基背风侧存在明显的涡旋流, 流场情况较为复杂, 风速都呈现上升、下降、再上升的趋势; 在0~5 m范围内, 普通路基的风速始终低于5 m·s-1, 因此最可能导致沙粒在坡脚沉积。

图 4 铁路路基特征断面风速廓线(风速15 m·s-1) Figure 4 Wind speed profile for the feature section of railway subgrades(wind speed 15 m·s-1) a:迎风坡脚 Windward slope; b:路基顶部 Subgrade top; c:背风坡脚 Leeward slope.

边坡侵蚀和铁轨磨蚀主要是由于近地层气流中沙粒对边坡和铁轨的冲击(邹维, 2012), 因此分析路基特征表面30 cm高度处的风速变化情况。从表 1可以看出, 普通路基4个特征位置处的风速, 都比入口速度(15 m·s-1)低, 说明在路基表层有一层低速区域, 速度低于起沙风速, 所以普通路基边坡积沙较为严重。通风路基迎风侧坡脚风速减小约64%, 路肩风速增加约19%, 气流在迎风坡上升时, 速度显著增大。输沙路基迎风侧坡脚风速减小约36%, 路肩风速增加约25%, 在3种路基中对气流的阻碍作用最小, 因此不易在路基边坡形成积沙。背风侧坡脚风速减小约20%, 由于输沙路基输沙段促使气流加速, 所以在背风侧坡脚风速仍较高, 经模拟证明该路基在单风向风沙较大的区域具有良好的适用性。

表 1 距路基床面30 cm高度特征位置风速变化 Tab.1 Variations of wind speed for the subgrade bed height of 30 cm
2.3 路基边坡积沙分布特征

青藏高原的风向与铁路的走向有一定夹角, 而且很多地段处于多风向地带, 不同的风向会影响路基的积沙情况。本文主要考虑风向与铁路垂直时(风沙流速度为15 m·s-1)的坡面积沙情况, 这对多风向的实际情况具有一定的借鉴意义, 并能反映出边坡积沙的普遍规律。

图 5是风沙流经过时路基周围积沙分布图(沙粒的体积分数表示该区域沙粒体积占总体积的比值, 其值越接近1表示该区域积沙量越大。不同颜色代表沙粒不同的体积分数, 红色最大, 蓝色最小, 红色表示沙粒已沉积在此处, 其他颜色代表沙粒仍在继续运动)。从图 5可以看出, 在风沙流经过路基时, 3种路基均有不同程度的积沙, 且路基迎风坡积沙量明显大于背风坡。普通路基边坡积沙最为严重, 迎风侧路堤和路肩都存在积沙现象, 路基近地层气流中沙粒含量较大, 导致风沙流持续吹蚀边坡及铁轨。通风路基的迎风坡脚和背风坡脚均有部分积沙, 通风管上部的边坡没有积沙, 但是通风管内积沙严重, 可能会造成通风管堵塞, 长此以往必定会使其保护冻土的作用减弱, 甚至失效。输沙路基边坡没有积沙, 当风沙流经过路基时, 一些沙粒沉积在输沙段内, 表层沙粒继续运动, 最终在积沙池内沉积, 其余沙粒随气流移动, 离开路基。整个过程中路基边坡以及顶部无沙粒沉积, 可见输沙路基在单风向风沙较严重的地区, 可以起到减轻风蚀的作用, 并且人工清理积沙池内的积沙也相对便利。

图 5 不同路基积沙分布(风速15 m·s-1) Figure 5 Distribution of sand deposition under different subgrades(wind speed 15 m·s-1) a.普通路基 Ordinary subgrade; b.通风路基 Ventilated subgrade; c.输沙路基 Sand-transmitting subgrade;图 6, 8同此 The same as in Fig. 6, 8.图中色阶表示沙粒相体积分数(m3·m-3), 下同。Color scale indicates the volume fraction of the sand(m3·m-3), and the same as below.
图 6 布设沙障时铁路路基流场分布特征(风速15 m·s-1) Figure 6 Distribution characteristics of flow field in railway subgrades with sand barriers(wind speed 15 m·s-1) b, c中的PE阻沙网、栅栏沙障(轨枕式挡沙墙)设置同a。Laying of sand blocking PE nets, sand fence (sleeper-type sand barrier)in b, c are the same as in a.
图 8 布设沙障时路基积沙分布(风速15 m·s-1) Figure 8 Distribution of sand deposition on subgrades with sand barriers(wind speed 15 m·s-1)
2.4 沙障对路基流场及积沙情况的影响

根据青藏铁路南山口路段的实际调查情况, 该路段在路基上风向布设了多种沙障组合, 一共5个分区(姜鑫贵等, 2016)。本部分主要分析轨枕式挡沙墙(以下简称栅栏沙障)和PE阻沙网对路基流场及积沙的影响。

本模拟中栅栏沙障高度2.1 m, 厚度0.3 m, 由废旧轨枕并排埋设而成。PE阻沙网高度为1.5 m, 整个沙障均匀透风, 孔隙度约为40%。从图 6可以看出, PE阻沙网和栅栏沙障布设在上风向时, 风沙流遇阻会有明显的减速过程, 2种沙障对迎风侧风场的影响范围均小于背风侧, 并且PE阻沙网对风场的影响范围略小于栅栏沙障, 这是因为PE阻沙网的高度小于栅栏沙障, 而且PE阻沙网具有一定孔隙度, 气流可以穿过阻沙网(唐玉龙, 2013)。风沙流离开沙障影响范围后, 速度开始缓慢增加, 当抵达路基迎风坡时, 又遇阻下降。与无沙障布设时相比, 风沙流在路基迎风坡的速度明显降低, 但由于沙障将大部分沙粒阻挡, 并且障后风速不足以起动沉积在背风侧的沙粒, 因此路基边坡的积沙量大大降低(杨印海等, 2010)。从模拟结果来看, 2种沙障之间存在一个低速区, 由于阻沙网对风沙流的阻碍作用, 阻沙网背风侧风沙流减速沉降, 随后在平坦开阔的戈壁地表流速有所回升, 但经过栅栏沙障时速度再次下降, 所以这个区域的风速呈现两端小、中间大的趋势, 这个规律在栅栏沙障与路基之间也是存在的, 只是由于两者间距较小, 所以并不明显。

图 7是2种沙障的积沙分布情况, 栅栏沙障迎风侧积沙量大于背风侧, PE阻沙网背风侧积沙量大于迎风侧, 辛国伟等(2016)以及郝晓杰(2012)也发现了这一结论。图 9为南山口沙障积沙情况, 可以看出沙障两侧实际积沙形态与本文模拟基本一致。从沙障的积沙情况来看, 栅栏沙障迎风侧积沙量较大, 并且还有部分沙粒继续以跃移和悬移的形式运动, 但由于风速很低, 最后沙粒会在迎风侧堆积。PE阻沙网迎风侧没有明显沙粒堆积, 大部分沙粒会通过孔隙直接穿越到阻沙网背风侧, 只有少部分沙粒会随气流攀升, 在PE阻沙网上部穿过, 最终沉积在背风侧。对比2种沙障的积沙分布特征可以看出, 栅栏沙障的拦沙效果明显比PE阻沙网好, PE阻沙网背风侧积沙量较少, 堆积高度很低, 有部分沙粒会在风力作用下继续运动。

图 7 2种沙障周围沙粒分布(风速15 m·s-1) Figure 7 Distribution of sand around two sand barriers(wind speed 15 m·s-1) a.轨枕式挡沙墙 Sleeper-type sand barrier; b.PE阻沙网 Sand blocking PE nets.
图 9 青藏铁路南山口段沙障积沙情况 Figure 9 Sand accumulation situation in Nanshankou section of Qinghai-Tibet Railway 栅栏沙障 Sand fence:a.迎风侧 Windward side; b.背风侧 Leeward side. PE阻沙网 Sand blocking PE: c.迎风侧 Windward side; d.背风侧 Leeward side.

图 8为布设沙障时路基积沙情况, 与图 5相比, 3种路基边坡的积沙量都相对减少。普通路基迎风坡的积沙情况有很大改观, 只有极少的沙粒堆积在坡脚, 风沙流沿着边坡抬升, 经过铁轨后在背风坡脚沉积, 其余沙粒仍处在运动状态, 随气流搬运远离路基。在通风路基上风向布设沙障后, 有部分沙粒在边坡移动, 但并未堆积, 且通风管内的积沙也明显减少, 贴近管道下部堆积的沙粒层变薄, 保温功能不会受到太大影响。输沙路基是这3种路基中积沙最少的, 整个路基结构几乎没有积沙, 积沙池内沙量较少, 积沙池后也有一小部分积沙, 人工清理相对便利。可以看出沙障的存在限制了沙粒的运动, 大量沙粒被拦截, 只有少部分沙粒被搬运至路基附近, 起到了很好的防沙效果。

3 讨论

本研究模拟结果显示, 随着风速的增加, 3种路基周围流场形态差异不大, 都具有明显的分区。张克存等(2010)对沱沱河路段路基流场结构进行了风洞模拟试验研究, 也发现气流在经过路基时形成明显分区。姚立强等(2013)对临策铁路戈壁段4 m高度路基的流场实测结果与本文模拟的流场分布相近。路基边坡的蚀积特征与气流速度有很好的对应关系, 在路基上方的高速区产生强烈风蚀, 沙粒主要堆积在路基两侧坡脚低速区(张军平等, 2011)。笔者通过模拟发现路基迎风侧坡脚的积沙量大于背风侧, 而不同种类的路基边坡积沙情况也有所不同。尤全刚等(2011)也指出路基背风坡坡度越小越易造成沙粒沉积, 并且随着风速增加沙粒沉积位置逐渐从背风坡中部向坡脚移动。风沙路段路基边坡的损害主要是风蚀破坏, 路基高度在一定程度上影响风蚀以及沙粒的堆积位置(苑宝华等, 2012; 刘伟伟等, 2006)。笔者提出的输沙路基主要起输导作用, 气流经过输沙管道时先迅速增大再缓慢减小, 管道的长宽比会影响风速增幅(洪新民等, 2017), 可根据实际情况适当调整以保证风速大于起沙风, 携沙气流始终从输沙管通过, 因此风沙流对路基边坡的风蚀及对铁轨的磨蚀作用有所减轻, 边坡积沙量明显减少。

高立式沙障的防风固沙作用主要体现在降低风速和促使沙粒从空气中分离并沉积(刘世海等, 2010)。阻沙栅栏是青藏铁路沿线使用最多的防护措施之一, 主要是因为其积沙量较多, 防沙效益明显(谢胜波等, 2014)。笔者的模拟结果表明, 栅栏沙障对风速的减弱作用更强, 影响范围也大于PE阻沙网, 而且栅栏沙障拦沙能力强于PE阻沙网, PE阻沙网输沙率较高, 障前几乎没有积沙, 障后沙粒的堆积高度较低, 更容易输送至路基处, 因此在布设防护措施时可适当增加栅栏沙障的数量, 提高防沙效益。李凯崇等(2015)对几种不同的挡沙墙进行了风洞试验研究, 其结果也表示轨枕式挡沙墙对风能的削弱能力远大于高立式PE网沙障, 并且PE网沙障输沙率大于其他高立式沙障。马学喜等(2016)发现不同结构的HDPE (High density PE)阻沙网都能阻拦大量流沙, 积沙主要分布于阻沙网的下风侧, 而上风侧积沙量相对较少。上述研究结果与本文模拟的结果近似, 说明模拟的结果基本可以反映路基及沙障的流场和积沙分布的实际情况。

本研究中, PE阻沙网后积沙范围在2H内, 并且阻沙网前几乎没有积沙, 这与郝晓杰等(2012)张克存等(2016)的研究略有出入, 特别是障后积沙范围差别较大。郝晓杰的研究指出PE阻沙网后积沙范围在10H左右, 张克存等人(2016)利用Leica C10三维激光扫描仪精确测量了北麓河段两侧典型防沙措施及防护体系积沙形态特征和风沙堆积量, 其结果表明尼龙网阻沙栅栏障后积沙范围可达5H左右但障前并无积沙。本研究模拟的风速为固定值, 但实际情况下风速会不断变化, 风速较小的时候沙粒会堆积在靠近阻沙网的位置, 大风会把靠近阻沙网的沙粒继续向后输送并沉积, 这会导致PE阻沙网后积沙范围越来越大, 最终达到阈值。同时, 由于研究所选的材料及疏透度不同, 也可能会造成研究结果的不一致。

尽管数值模拟结果显示笔者设计的输沙铁路路基及沙障的组合形式明显减少了路基周围积沙量, 然而目前并未在青藏铁路实际采用, 尚无法实际验证其疏导效果, 需要在今后的研究中继续探索其优缺点, 在实践中进一步检验和完善。

4 结论

3种路基在无防护措施的情况下, 均有不同程度的积沙, 普通路基积沙情况最严重, 迎风坡积沙量明显大于背风坡; 通风路基边坡积沙较少, 部分沙粒沉积在通风管, 破坏通风管的保温功能; 输沙路基边坡几乎没有积沙, 并且输沙管内的积沙可以通过风力作用输移到积沙池内集中处理, 节省大量人力物力, 在实际应用中可以更好地避免风沙危害带来的损失。

在布设有沙障的情况下, 3种路基边坡的积沙量都明显减少, 沙障在一定程度上阻挡了风沙流, 减少风沙对路基的危害。沙障与路基之间存在一个低速区, 该区域速度两端低, 中间高, 根据这个规律, 可以适当调整沙障布设的位置, 使二者之间的过渡区域风速无法达到起沙风速, 保证沙障背风侧积沙不被再次起动, 提高沙障的防沙效益, 减轻风沙对路基的侵蚀。

模拟结果表明, 输沙路基相对于普通路基和通风路基更适合布设在青藏铁路风沙路段, 其独特的设计可有效地增加来流风速, 防止风沙上道, 并且有专门的积沙池收集堆积的沙粒, 但沙粒不能完全沉积在积沙池内, 会有少量沙粒沿积沙池攀升, 堆积在积沙池后。输沙路基的设计为青藏铁路防沙治沙提供了一种新的思路, 从研究结果来看, 此设计行之有效, 可在实践中进行尝试, 为防治沙害做出贡献。

参考文献(References)
北京林业大学. 2013. 输沙铁路路基. 中国专利, 2013050800677430.
(Beijing Forestry University. 2013. Sand-transmitting railway subgrade. Chinese Patent: 2013050800677430. [in Chinese])
郝晓杰, 熊治文, 蒋富强, 等. 2012. 青藏铁路不同防沙栅栏的布设位置研究[J]. 铁道标准设计, 56(4): 19-23.
(Hao X J, Xiong Z W, Jiang F Q, et al. 2012. Research on the positions of various kinds of sand barriers in Qinghai-Tibet railway[J]. Railway Standard Design, 56(4): 19-23. [in Chinese])
洪新民, 郭文华, 熊安平. 2017. 山区峡谷风场分布特性及地形影响的数值模拟[J]. 长安大学学报:自然科学版, 37(5): 56-64.
(Hong X M, Guo W H, Xiong A P. 2017. Numerical simulation of distribution characteristic of wind fields and terrain's influence in mountain canyon[J]. Journal of Chang'an University:Natural Science Edition, 37(5): 56-64. [in Chinese])
姜鑫贵, 周金星, 刘玉国, 等. 2016. 青藏铁路南山口段沙害防治措施的效果[J]. 林业科学, 52(10): 55-63.
(Jiang X G, Zhou J X, Liu Y G, et al. 2016. Effectiveness of sand control measures in Nanshankou Area along Qinghai-Tibet railway[J]. Scientia Silvae Sinicae, 52(10): 55-63. DOI:10.11707/j.1001-7488.20161007 [in Chinese])
李驰, 黄浩, 孙兵兵. 2010. 风沙环境下沙漠路基风蚀破坏数值模拟研究[J]. 岩土力学, 31(2): 378-382.
(Li C, Huang H, Sun B B. 2010. Numerical simulation studies on wind erosion damage for desert subgrade in sandstorm environment[J]. Rock & Soil Mechanics, 31(2): 378-382. [in Chinese])
李凯崇, 薛春晓, 刘贺业, 等. 2015. 不同类型挡沙墙风沙防护机理的风洞实验研究[J]. 铁道工程学报, 32(1): 17-21.
(Li K C, Xue C X, Liu H Y, et al. 2015. Wind tunnel test on sand-preventing mechanism of different kinds of sand-barriers[J]. Journal of Railway Engineering Society, 32(1): 17-21. [in Chinese])
李晓军, 蒋富强. 2016. 风沙流对戈壁地区挡风墙响应规律的数值模拟分析[J]. 铁道标准设计, 60(3): 47-51.
(Li X J, Jiang F Q. 2016. Numerical simulation analysis of response law of wind-blown sand flow around wind-break wall in Gobi Area[J]. Railway Standard Design, 60(3): 47-51. [in Chinese])
刘世海, 冯玲正, 许兆义. 2010. 青藏铁路格拉段高立式沙障防风固沙效果研究[J]. 铁道学报, 32(1): 133-136.
(Liu S H, Feng L Z, Xu Z Y. 2010. Study on effect of wind erosion controlling in the Geerm-Lahsa section of the Qinghai-Tibet railway[J]. Journal of the China Railway Society, 32(1): 133-136. [in Chinese])
刘伟伟, 崔立臣. 2006. 风沙流速与沙漠路基合理填土高度分析[J]. 青岛理工大学学报, 27(3): 140-143.
(Liu W W, Cui L C. 2006. The analysis between sand flow rate and filling height of roadbed in desert regions[J]. Journal of Qingdao Technological University, 27(3): 140-143. [in Chinese])
刘志强, 赖远明, 张淑娟, 等. 2006. 青藏铁路通风路基三维温度场的数值分析[J]. 土木工程学报, 39(1): 70-74.
(Liu Z Q, Lai Y M, Zhang S J, et al. 2006. Numerical analysis for the ventilated embankment in Qing-Tibet railway[J]. China Civil Engineering Journal, 39(1): 70-74. [in Chinese])
卢玉华, 曾雪兰, 黄润兰. 2012. 障碍物对山坡地形风场的影响分析[J]. 可再生能源, 30(4): 33-36.
(Lu Y H, Zeng X L, Huang R L. 2012. Influence analysis of obstacle on wind field over hilly terrain[J]. Renewable Energy Resources, 30(4): 33-36. [in Chinese])
马学喜, 李生宇, 何志辉, 等. 2016. 高立式阻沙网的阻沙效果对比[J]. 中国沙漠, 36(4): 925-931.
(Ma X X, Li S Y, He Z H, et al. 2016. Comparison of the sand-blocking effects of HDPE windbreaks[J]. Journal of Desert Research, 36(4): 925-931. DOI:10.7522/j.issn.1000-694X.2015.00120 [in Chinese])
牛清河, 屈建军, 张克存, 等. 2009. 青藏铁路典型路段风沙灾害现状与机械防沙效益估算[J]. 中国沙漠, 29(4): 596-603.
(Niu Q H, Qu J J, Zhang K C, et al. 2009. Status of aeolian-sand disaster and estimation of mechanical sand-controlling benefit at typical sections of Qinghai-Tibet railway[J]. Journal of Desert Research, 29(4): 596-603. [in Chinese])
邱青云, 魏庆朝, 吴西良. 2004. 青藏铁路唐拉段地质灾害的评价分析[J]. 中国安全科学学报, 14(5): 3-5.
(Qiu Q Y, Wei Q C, Wu X L. 2004. Evaluation analysis on geological disasters at Tang-La section of Qing-Zang railway[J]. China Safety Science Journal, 14(5): 3-5. [in Chinese])
石龙, 蒋富强, 韩峰. 2014. 风沙两相流对铁路路堤响应规律的数值模拟研究[J]. 铁道学报, 36(5): 82-87.
(Shi L, Jiang F Q, Han F. 2014. Numerical simulation of response law of wind-blown sand flow around the railway embankment[J]. Journal of the China Railway Society, 36(5): 82-87. [in Chinese])
唐玉龙. 2013. 青藏铁路西格段戈壁风沙流防治体系研究[J]. 中国沙漠, 33(1): 72-76.
(Tang Y L. 2013. Wind-blown sand prevention system at Gobi Desert Area along Xining-Golmud section of Qinghai-Tibet railway[J]. Journal of Desert Research, 33(1): 72-76. DOI:10.7522/j.issn.1000-694X.2013.00010 [in Chinese])
谢胜波, 屈建军. 2014. 青藏铁路沿线地形、气候、水文特征及其对沙害的影响[J]. 干旱区资源与环境, 28(10): 157-163.
(Xie S B, Qu J J. 2014. The effects of landform and the climatic and hydrological characteristics on the sand hazards along Qinghai-Tibet railway[J]. Journal of Arid Land Resources & Environment, 28(10): 157-163. [in Chinese])
谢胜波, 屈建军, 刘冰, 等. 2014. 青藏铁路沙害及其防治研究进展[J]. 中国沙漠, 34(1): 42-48.
(Xie S B, Qu J J, Liu B, et al. 2014. Advances in research on the sand hazards and its controls along the Qinghai-Tibet railway[J]. Journal of Desert Research, 34(1): 42-48. DOI:10.7522/j.issn.1000-694X.2013.00087 [in Chinese])
辛国伟, 程建军, 景文宏, 等. 2016. 来流廓线对风沙流场和风沙堆积影响的数值模拟—以挡沙墙为例[J]. 干旱区研究, 33(3): 672-679.
(Xin G W, Cheng J J, Jing W H, et al. 2016. Numerical simulation of the influence of incoming flow profile on sand flow field and aeolian Sand deposition:a case study around sand retaining wall[J]. Arid Zone Research, 33(3): 672-679. [in Chinese])
薛智德, 刘世海, 许兆义, 等. 2010. 青藏铁路措那湖沿岸防风固沙工程效益[J]. 北京林业大学学报, 32(6): 61-65.
(Xue Z D, Liu S H, Xu Z Y, et al. 2010. Efficiency of wind erosion control measures along Cuona Labe on Qinghai-Tibet railway, western China[J]. Journal of Beijing Forestry University, 32(6): 61-65. [in Chinese])
杨印海, 蒋富强, 王锡来, 等. 2010. 青藏铁路错那湖段沙害防治措施研究[J]. 中国沙漠, 30(6): 1256-1262.
(Yang Y H, Jiang F Q, Wang X L, et al. 2010. Prevention and management of blown-sand hazards along Qinghai-Tibet railway in Cuonahu Lake Area[J]. Journal of Desert Research, 30(6): 1256-1262. [in Chinese])
尤全刚, 薛娴, 王涛, 等. 2011. 戈壁地区公路防沙措施防沙效应的风洞试验[J]. 中国沙漠, 31(3): 550-557.
(You Q G, Xue X, Wang T, et al. 2011. Wind tunnel experiment of the effect of sand prevention measures along Gobi highway[J]. Journal of Desert Research, 31(3): 550-557. [in Chinese])
于云江, 鲁春霞. 2001. 青藏线客城区段沙害成因与防治途径的研究[J]. 干旱区地理, 24(3): 236-242.
(Yu Y J, Lu C X. 2001. Studies on forming cause of sand disaster and sand-controlling approach at Ke-Cheng section in Qinghai-Tibet railway[J]. Arid Land Geography, 24(3): 236-242. [in Chinese])
于云江, 史培军, 鲁春霞. 2001. 青藏铁路客城区段风沙流特点及沙害防治措施的研究[J]. 自然灾害学报, 10(1): 30-36.
(Yu Y J, Shi P J, Lu C X. 2001. Research on characteristics of wind-sand fluid and sandy disaster reduction at Kecheng(K454-459) section along Qing-Zang railway[J]. Journal of Natural Disasters, 10(1): 30-36. [in Chinese])
苑宝华, 张立刚, 杨有海, 等. 2013. 戈壁地区高速铁路路基边坡防护措施研究[J]. 路基工程, 30(2): 65-68.
(Yuan B H, Zhang L G, Yang Y H, et al. 2013. Study on slope protection measures for subgrade of high-speed railway in Gobi region[J]. Subgrade Engineering, 30(2): 65-68. [in Chinese])
姚立强, 左合君, 李钢铁, 等. 2013. 临策铁路戈壁段路堤二维风速流场特征研究[J]. 铁道标准设计, 57(9): 4-9.
(Yao L Q, Zuo H J, Li G T, et al. 2013. Research on the characteristics of two-dimensional wind velocity field of embankment of Gobi section of Linhe-Ceke railway[J]. Railway Standard Design, 57(9): 4-9. [in Chinese])
张军平, 王引生, 蒋富强. 2011. 兰新铁路戈壁地区路基周围风沙流运动特征数值分析[J]. 中国铁道科学, 32(4): 14-18.
(Zhang J P, Wang Y S, Jiang F Q. 2011. Numerical analysis on the features of sand flow movement around the embankment of Lan-Xin railway in Gobi region[J]. China Railway Science, 32(4): 14-18. [in Chinese])
张克存, 安志山, 庞营军, 等. 2016. 青藏铁路北麓河路段风沙防护体系阻沙效益[J]. 中国沙漠, 36(5): 1216-1222.
(Zhang K C, An Z S, Pang Y J, et al. 2016. Benefit of wind-blown sand prevention system in the Beiluhe section of the Qinghai-Tibet railway[J]. Journal of Desert Research, 36(5): 1216-1222. DOI:10.7522/j.issn.1000-694X.2015.00148 [in Chinese])
张克存, 牛清河, 屈建军, 等. 2010. 青藏铁路沱沱河路段流场特征及沙害形成机理[J]. 干旱区研究, 27(2): 303-308.
(Zhang K C, Niu Q H, Qu J J, et al. 2010. Study on the characteristics of flow field and the mechanism of wind-blown sand disasters in the Tuotuohe region along the Qinghai-Tibet railway[J]. Arid Zone Research, 27(2): 303-308. [in Chinese])
张克存, 屈建军, 牛清河, 等. 2011. 青藏铁路沿线阻沙栅栏防护机理及其效应分析[J]. 中国沙漠, 31(1): 16-20.
(Zhang K C, Qu J J, Niu Q H, et al. 2011. Protective mechanism and efficiency of sand-blocking fences along Qinghai-Tibet railway[J]. Journal of Desert Research, 31(1): 16-20. [in Chinese])
朱红钧, 林元华, 谢龙汉. 2010. FLUENT流体分析及仿真实用教程[M]. 北京: 人民邮电出版社.
(Zhu H J, Lin Y H, Xie L H. 2010. FLUENT fluid analysis and simulation practical tutorial[M]. Beijing: Posts & Telecom Press. [in Chinese])
邹维. 2012. 风沙运动规律的分析和研究[J]. 中国水土保持, 32(5): 44-46.
(Zou W. 2012. Analysis and study of wind blow sand movement[J]. Soil and Water Conservation in China, 32(5): 44-46. [in Chinese])
Zhang K C, Qu J J, Liao K T, et al. 2010. Damage by wind-blown sand and its control along Qinghai-Tibet railway in China[J]. Aeolian Research, 1(3): 143-146.