乌玛高速公路中卫段风沙环境及输移规律

　　摘要：作为宁夏自治区南北向交通运输通道，乌玛高速公路对于带动区域经济发展有着重要意义。为保证公路正常运行，该文采用野外定位监测和室内分析方法，揭示了公路沿线风沙活动规律。结果表明：(1)区域0～5m/s等级的风速频次占全年风速的90.45%，起沙风向主要集中在西北、东和偏南方向，合成输沙方向为偏东南方向;研究区输沙势小于200VU，属于低风能环境。(2)沿乌玛高速公路自东北向西南方向，所有样品均以细沙和中沙为主，粒径范围集中在1.0～2.76Φ。除5号采样点20cm，50cm和100cm深度沉积物的百分含量分布曲线呈双峰态，其余样品全部为单峰态，表明研究区域表层沉积物经过了充分分选。随着深度的增加，不同粒径范围百分含量的变化规律并非十分明显，即区域风能环境比较稳定。(3)所有样品分选性较好或很好，偏度为正偏或极正偏，粒径偏粗。(4)区域内风沙流主要集中在30cm以下，占总输沙量的89.46%，地表0～20cm高度内的输沙量占总输沙量的63.81%。基于防护体系走向与区域合成输沙方向垂直时可取得较好防护效益，建议乌玛高速防护体系走向设计为东北—西南，且固沙措施的高度不低于20cm。

　　关键词：乌玛高速公路;输沙势;粒度特征;风沙流

　　乌玛高速公路是内蒙古乌海至青海玛沁高速公路的简称，设计时速100km/h，采用双向四车道标准建设，途径内蒙、宁夏、甘肃和青海四省份，是“一带一路”战略重要的运输通道，其建成与通车对于缓解宁夏自治区运输压力有着重要意义。但穿越的腾格里沙漠由于流动沙丘广布，其对乌玛高速公路的建设、维护及安全运营带来了极大挑战。

　　高速公路论文范例： 城市高速铁路明挖隧道装配式结构设计方案研究

　　虽我国有包括塔克拉玛干石油公路、京新高速公路等穿越沙漠或边缘公路的既有建设和维护经验，但前者或公路设计等级略低，或风沙环境略弱，加之区域环境的异质性，故以往既得经验无法直接应用在乌玛高速公路防护体系的设计之中[1-2]。自包兰铁路在腾格里沙漠东南缘修建之日起，研究者对区域的风沙环境、沙丘类型及移动规律等进行了详细分析，为包兰铁路沙坡头段的建设和正常运行提供了理论依据[3-8]。

　　后期随着包兰铁路防护体系的建设与运营，铁路防护措施的防护效益、防护机理以及周边生态环境效益逐渐得到明晰[9-16]，成为我国防沙治沙工程的标杆。包兰铁路沙坡头段防护体系的有效性为乌玛高速公路的建设和维护提供可借鉴的现有经验，充分验证了乌玛高速公路的建设可行性。但受区域风况环境、地形等因素影响，乌玛高速公路所经区域的风能环境、沉积物及沉积环境均有别于前者，甚至存在明显差异。因此亟需对乌玛高速公路沿线风能环境、沉积物特征等风沙输移规律进行研究，继而为乌玛高速公路的建设及沿线风沙防护体系的设计提供参考和指导，具有重要的科学和实践意义。

　　1研究区概况和数据来源

　　腾格里沙漠是我国第四大沙漠，常年受蒙古高压控制，干旱少雨，区域内流动沙丘广布，沙丘类型以格状沙丘为主[1，3，17]，线性沙垄和金字塔沙丘偶有分布。乌玛高速公路腾格里沙漠段共计21km，自东北向西南方向延伸。本文采样点共设置5个，两两相距4km。

　　其中，1号采样点东侧和东北侧为太阳能光伏电厂，2号、3号、4号采样点为流动沙丘，5号采样点靠近包兰铁路，属于包兰铁路防护体系内部，周边沙丘已经固定，植被覆盖度超过30%，植被以人工灌木花棒(Hedysarumscoparium)等为主。区域风速观测采用定位观测，仪器型号为HOBOU30，风速风向统一架设高度为2.0m。风沙流观测采用中科院沙漠与沙漠化重点实验室自制的口琴式集沙仪，仪器高60cm，共30层，单个进样口径2cm×2cm。沙样采集采用人工剖挖法，采集深度分别为5，10，20，50，100cm，使用环刀取样。

　　2区域风速分布

　　最大平均风速值出现在8月，7月次之，最小风速值出现在6月，仅为1.46m/s。总体来看，5月、6月的风速较小，到了7月风速开始突变，增加2.5倍左右，持续到9月结束，1—4月、10—12月风速变化较为平稳，平均为2.34m/s。与平均风速变化趋势相比，月平均起沙风(>5m/s)的分布规律略有差别，1—6月平均起沙风速变化幅度较小，7月、8月起沙风速达到最大值，之后开始减小。

　　从风速分布情况看，0～5m/s等级风速占全年风速的90.45%，仅有不足10%左右分布在其他风速等级，并且随着风速的增加，高风速频次的比重逐渐减小。进一 步分析可知，0～5m/s的风速频率最大月份为6月，接近20%，而本月平均风速为全年最小，可见6月的风速特征是频次高、风速低。5～6，6～7m/s等级的风速主3输沙势变化输沙势是衡量区域风沙活动强度及风沙地貌演变的重要指标，也是目前应用较为广泛的方法[1，3，7]。2019年输沙势为37.93VU，合成输沙势9.35VU，合成输沙方向107.46°，方向变率0.25。2020年输沙势较2019年有所减小，为19.17VU，合成输沙势7.02VU，合成输沙方向117.91°，方向变率。

　　4沉积物粒度特征

　　4.1沉积物空间分布

　　1号采样点5cm深度处沉积物的粒径主要分布在1.26～1.76Φ，百分含量为65.03%。随着深度的增加，上述范围内的颗粒百分含量有增加趋势，分别为70.68%，84.69%，74.03%，78.46%，且分布范围扩大至1.5～2.76Φ，但平均粒径呈减小趋势。进一步分析不同深度沉积物百分含量分布曲线可知5cm和50cm沉积物的百分含量分布曲线峰形较窄，其余深度沉积物分布曲线峰形较宽，但均为单0.37。

　　进一步分析可知2019年和2020年输沙势主要集中在西北方向、东方向和偏南方向，方向相对稳定，但强度略有变化。其中，2019年上述三个方向的输沙势分别为5.69VU，4.56VU和2.11VU，2020年则为4.16VU，2.42VU和1.32VU。可见研究区各个风向的输沙势虽有变化，但风能环境基本稳定，属于低风能环境，且风向较为稳定，主要集中西北、东和偏南方向。

　　峰态分布。从累计百分比曲线可以看出，5cm和50cm深度处沉积物的累计百分比增长速率最快，斜率最大，表明该深度的土壤粒径分布最为集中。与1号采样点沉积物百分含量分布曲线相比，2号采样点沉积物分布曲线也为单峰态分布，但峰态较1号窄，粒径分布范围也存在差异。其中，5cm深度处沉积物粒径范围集中在2.0～2.51Φ，百分含量为74.71%，其余深度分布范围集中在1.76～2.76Φ，分布范围较1号点宽，百分含量集中在70.77%～85.33%。与1号采样点累计百分比曲线相比，2号采样点所有深度沉积物的累计百分比曲线相似性较1号高。

　　5风沙流结构

　　风沙流结构是指沙量随高度的变化特征，是沙粒运动的基本形式和过程。风沙流结构可以反映沙粒的跃移高度和运动轨迹，是揭示区域风沙活动的重要参数，更是设置风沙工程的重要参考标志之一[3]。风沙流结构与地形、沙粒粒径组成、湿度密切相关，但气流的波动和速度对其影响更大。依据野外多梯度口琴式集沙仪的观测数据可知：研究区内的风沙流主要集中在30cm以下，占总输沙量的89.46%。地表0—20cm高度内的输沙量占总输沙量的63.81%，进一步拟合输沙量与高度的关系可知两者符合指数函数。

　　6讨论与结论

　　表层沉积物的粒径分布曲线及特征参数对于研究区域的风况环境和沉积环境有着重要指示意义[18]。研究结果表明：沿乌玛高速公路自东北向西南方向，沉积物粒径范围集中在1.0～2.76Φ，以细沙和中沙为主，不含粉沙和砾石。除5号采样点20cm、50cm和100cm深度沉积物的百分含量分布曲线呈双峰态，其余样品全部表现为单峰态，说明区域内沉积物经过了充分分选和运移，局地沉积环境由早期的两组主导风向转变为近期的单组风向，主要原因是包兰铁路防护体系的影响。随着深度的增加，1号采样点的粗沙百分含量减小，细沙和极细沙含量随深度增加呈先增加后减少趋势。

　　2号采样点中沙百分含量呈现上升趋势，细沙百分含量呈下降趋势。3号采样点和4号采样点沉积物以中沙和细沙为主，随着深度的增加，4号采样点中沙百分含量逐渐减小，细沙百分含量增加。5号采样点沉积物以细沙为主，各粒级百分含量随深度无明显变化规律。区域0～5m/s等级的风速频次占全年风速的90.45%，仅有不足10右的频次分布在其他等级，起沙风向集中在西北、东和偏南方向，合成输沙方向为偏东南方向，属于低风能环境。区域内沙粒的运移高度集中在30cm以下，且超过50%的沙粒在20cm高度以下运移，即本区域内沙粒的跳跃高度不超过20cm[3]，建议区域固沙措施的设置高度不低于20cm。

　　综上可知，沿乌玛高速公路腾格里沙漠段，沉积物粒径分布范围、百分含量及累计曲线存在一定的差异性，区域风况环境发生过变化。但从固定气象观测设备的分析结论可知上述差异性较小，对乌玛高速公路的影响范围和程度也弱，故在防护体系的设置中，设计者也应恰当考虑，建议收集当地长时间序列的气象资料，研究其风能环境变化趋势，并在公路沿线增加固定气象观测点，获取更多更详细的风况资料，继而为防护体系的设置和高速公路的安全运营保驾护航。

　　参考文献：

　　[1]黄磊，张余，何明珠，等.腾格里沙漠东南缘风沙环境及其退交通干线生态防护体系的影响[J].水土保持通报，2018，38(2)：135-139.

　　[2]左合君，董智，魏江生，等.沙漠地区高速公路防沙体系效益分析[J].水土保持学报，2005，12(6)：222-225.

　　[3]张克存.风沙活动特征及其工程防沙主要措施研究[D].北京：中国科学院大学，2006.

　　[4]罗雅曦，刘任涛，张静，等.腾格里沙漠草方格固沙林土壤颗粒组成、分形维数及其对土壤性质的影响[J].应用生态学报，2019，30(2)：525-535.

　　[5]张正偲，董治宝，管梦鸾.腾格里沙漠东南缘反向沙丘形态演化过程[J].中国沙漠，2018，38(4)：709-715.

　　作者：王金国1，安志山2，3，张克存2，屈建军2，何明珠2