-
铺设沙障是目前最有效的固沙方法之一,可以大幅度降低地表风速,减缓沙漠流动,在保护沙漠工程、遏制沙漠侵蚀方面起着重要作用,其中草沙障在国内应用最为广泛。草沙障材料一般为稻草和麦草,芦苇杆也有少量应用,平面形状为正方形网格,边长主要在1~2 m之间,也被称为半隐蔽式沙障。沙障的防风效益很早就引起科研人员的关注,是评价沙障固沙能力的重要指标。沙障方格尺寸是影响沙障防风效能的重要参数,许多研究工作都集中于沙障尺寸方面的对比性研究,通过改变沙障网格的边长高度参数,实验各种尺寸沙障的防风效能[1-4]。沙障的材料也是影响沙障防风固沙能力的重要影响因素,目前已有大量的新材料被引入沙障工程中[5-8],并得到相应的防风效益研究,如尼龙网沙障[9-10]、PLA沙障、砾石[11]等。沙障的防风固沙原理与沙障周围的流场结构也是沙障的一个重要的研究方向[12-15]。目前,格状沙障铺设的网格形状都采用正方形,对于沙障的研究也都是基于正方形网格。从平面几何学角度来看,能够铺满地面的正多边形还包括三角形和六边形两种,非正多边形还包括密铺五边形,但对于这种多边形网格形状的沙障目前还鲜有人研究。有鉴于此,本文针对网格面积相同的三角形、四边形和六边形草沙障的防风效益进行观测和分析,为草沙障的工程施工优化设计提供理论支撑。
-
该试验区处于38°08′15″ E、102°49′57″ N,位于腾格里沙漠南缘,行政区划隶属于武威市。该地区属于典型的温带大陆性干旱气候,年平均风速3~4 m/s,年平均气温7.8 ℃,年平均降水量102.9 mm,年均蒸发量2 258.8 mm,日照时数3 181 h,终年盛行西南风,主要害风为西北风。主要植被为沙米(Agriophyllum squarrosum)、梭梭(Haloxylon ammodendron)、沙蒿(Artemisia desertorum)等。
-
2016年11月,在研究区选择经过平整后的沙地,植被覆盖率接近0,进行草沙障铺设,如图 1所示。为方便施工,沙障材料使用稻草编织的草帘。设计的沙障网格为正三角形、正方形和正六边形,其边长分别为1.52、1和0.62 m,网格尺寸面积为1 m2。选择质量较好孔隙度为10%以下的草帘作为沙障材料。为保证方格尺寸精确,制作正三角形和正六边形模具,用铁锹按模具位置在沙地表面开沟槽,将草帘立于沟槽中,两边埋沙,为保证沙障稳定,沙障埋入10 cm,草头高度为20 cm。
-
应用风杯式风速仪进行风速风向测量。风速观测高度为20、30、50和100 cm 4个高度。由于沙障高度为20 cm,沙障内部流场相对复杂,因此没有选择10 cm测量高度,观测点设置在每个方格的形心位置,沿西南方向排布。风速仪每两秒自动记录风速1次,选择风向为西南风向并且风速稳定的一段数据作为最终的分析数据。每种沙障的数据为200组,时间为6 min 40 s。草沙障监测点布置如图 2所示。
-
防风效能定义为降低的风速比率,按下式计算:
$$ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {E_h} = \left( {{V_h} - {U_h}} \right)/{V_h} \times 100\% {\kern 1pt} $$ (1) 式中:Eh为地表高度为h处的风速降低比率,%;Vh为无沙障地表高度为h处的平均风速,m/s;Uh为沙障内部高度为h处的平均风速,m/s。
-
下垫面高低不平的起伏使得近地表的流场结构变化很大,主要用地表粗糙度这一指标来衡量地表的起伏程度[16]。地表粗糙度定义为下垫面平均风速为0时的高度[17]。0风速一般出现在平均表面粗糙因素以上的某个高度,本实验选择在风速较为稳定的条件下测量地表粗糙度,粗糙度按下式进行计算:
$$ \lg {Z_0} = \left( {\lg {Z_2} - A\lg {Z_1}} \right)/\left( {1 - A} \right) $$ (2) 式中:Z0为地表粗糙度,cm;A=u2/u1,u1和u2分别为高度为Z1和Z2处的风速,m/s;Z1=100 cm,Z2=20 cm。
-
沙障可以有效的减弱地表风速,对地表风沙流有很强的削弱作用,当来流风沙进入沙障时,由于风速的下降及沙障地表产生的涡流结构,使沙粒可以很快的沉降下来;另一方面,沙障的存在可以有效减弱沙障内部地表风速,抑制沙障内部风沙流的形成,因此,沙障对来流风速的减弱作用是评判沙障固沙效果的重要指标。通过对沙障地表和裸露沙地风速间的比较,来评判沙障的防风效能的高低。
表 1为3种沙障的防风效益,由表 1可以看出,不同的沙障形态对地表风速的减弱程度各不相同,沙障的防风效能随测量高度的增加而减小。对于网格形状不同,单个网格面积为1 m2的沙障,三角形沙障和六边形沙障的防风效能都显著高于正方形沙障。高度为0.2 m时,正方形沙障的防风效能比三角形和六边形沙障分别低12%和8%,高度为0.3 m时,分别低11%和10%,高度为0.5 m时,分别低9%和4%。高度为1 m时,正方形沙障防风效能与三角形沙障相同,但显著好于六边形沙障。可以看出,四边形方格沙障的防风效能最差。网格面积相同的沙障,防风效能差异显著,这是因为网格形状不同的沙障周围流场结构也各不相同。目前对于地表粗糙元周围流场结构的研究都基于风洞环境和PIV系统,且地表障碍都采用小尺寸的缩尺度模型,高度10 cm以下的居多[18-19]。对于足尺寸沙障模型还无法精确的测出流场速度分布图,由于三维模型周围流场的复杂性,无法精确测量,目前没有人对大几何尺寸的三维模型周围流场结构进行实验研究。对于二维缩尺度沙障的周围流场结构,国内外学者基于风洞实验有较为深入的认识。研究表明,地表障碍对流场的影响越大,障后湍流强度越高,障碍对气流能量的削弱作用越大[20]。对于三角形网格沙障,相邻直立的壁面的间距沿着顺风向不断的变化,导致网格内部产生多个气流的加速区与减速区,流场结构相对于正方形网格沙障更为复杂,地形的复杂程度相对于正方形更高,因而可以比正方形网格沙障对来流风沙造成更大的扰动。对于六边形沙障,网格的交错排布也同样表现出了壁面间距不断收窄变宽的效应,加强了对风沙流的扰动。而对于正方形网格,壁面间距沿顺风向等距排布,地表地形的复杂程度相对简单,每个网格内部更容易形成相似的流场结构,虽然方便施工,但也减弱了对来流风的扰动程度。
表 1 多边形沙障防风效益对比
Table 1. Comparison in windproof efficiency among polygon sand barrier types
% 测量高度
Measuring height/m沙障类型Sand barrier type 三角形
Triangular正方形
Square六边形
Hexagonal0.2 67±16a 55±15b 63±13c 0.3 46±13a 35±16b 45±12a 0.5 25±9a 16±11b 30±8c 1 7±5a 7±3a 3±4b 注:不同字母表示不同类型沙障防风效能差异显著(P < 0.05)。Note: different letters mean significant difference in windproof efficiency among three types of sand barrier at P<0.05 level. -
风速廓线是指风速随高度变化。研究表明,平坦地面的风速廓线服从对数函数分布。当地面存在植被、沙障等障碍时,由于地面气流结构发生改变,风速廓线将受到地表湍流结构影响而发生变化。对距离沙障边界约10 m处的风速廓线进行监测,取裸露地表 1 m高度处4.5~5.5 m/s间的风速进行分析。图 3为3种沙障的风速廓线对比图。
如图 3所示,沙障的存在改变了地表风速廓线,有效的削弱了地表风速。无沙障存在时,地表风速随高度变化服从对数分布。当地表铺设沙障时,风速沿高度分布更接近指数函数分布。在0.2~0.5 m高度范围内,三角形沙障和六边形沙障表现出相近的变化特征,且各高度风速显著小于正方形沙障风速。高度为0.2~0.3 m时,3种沙障的风速增长较为迅速,增长速度大致相同。高度为0.5~1 m时,由于地表形态结构不同,风速增长速度具有显著的差异。风速增长速度为正方形网格沙障 < 三角形网格沙障 < 六边形网格沙障,但都大于无沙障时的风速增长速度。这是因为近地面受沙障影响较大,风速大幅度降低,随着高度的增加,气流被抬升而速度增加,导致垂直方向速度梯度大于无沙障情况。较快的风速增长梯度表明,地表沙障形成的流场结构具有更大的气流抬升效应,对上层流场影响相对较强,沙障对地表风沙流的能量消耗作用相对也更强。从图 3也可以看出,四边形的规则图形,形成的附面层气动结构对风速减弱作用小于三角形和六边形沙障,高度在1 m时,风速增长缓慢,接近裸露地表风速沿高度的增长速度。研究表明,裸露沙地的输沙行为主要集中于10 cm高度以下,风沙流进入高度为0.2 m高的沙障范围后,风沙流受到沙障的抬升作用而改变,输沙20~30 cm高度间的输沙率迅速增加[21],因此,距离地表 30 cm范围内的防风效能对降沙固沙有重要的意义,传统的正方形沙障降沙能力小于三角形和六边形沙障。
-
地表粗糙度是衡量沙障防风固沙能力的一个重要指标,研究表明裸露表的粗糙度随着风速的增加而增加,当风速大于9 m/s时,可以增长到10-2 m数量级,因此风速对粗糙度有较大影响[21]。本文选择裸露地表 1 m高度处,风速范围为1~6 m/s之间的数据进行分析,计算此风速范围内的地表粗糙度。3种沙障的地表粗糙度如图 4所示。
同一地形特征下粗糙度随风速的变化规律很早就引起研究人员的关注,粗糙度虽然用来表征地表的粗糙程度,理论上为一个表征特定地形的常量值,但野外观测表明,地表粗糙度随风速的增加而增加,原因为随着地表风速的提高,地表输沙量增加,导致地表粗糙度也相应提高。而铺设沙障的地表,输沙效应对粗糙度的影响相对较弱,对地表粗糙度的影响因素转变为草头在风力作用下的弯曲扰动和沙障孔隙在风力作用下的涡流结构变化。研究表明[21],铺设沙障后,低风速条件下,地表粗糙度随风速增长快速增加,当风速大于7 m/s时,地表粗糙度的增长速度大幅度减小。由于本文分析风速小于6 m/s,可以用线性模型进行回归分析。如图 4所示,对于同一类型沙障,地表粗糙度随风速的增大而增大。对3种沙障的地表粗糙度随风速变化进行线性拟合。结果表明,三角形沙障的地表粗糙度最大,六边形沙障地表粗糙度次之,正方形沙障地表粗糙度最小。3种沙障的地表粗糙度随风速变化的速度也不相同,三角形沙障的地表粗糙度随风速增长速度最快,六边形沙障次之,正方形沙障的粗糙度随风算变化最慢。裸露地表 1 m高处风速为4 m/s时,裸露地表粗糙度为0.166 cm,三角形沙障、正方形沙障和六边形沙障的地表粗糙度分别为11.5、9.3和10.4 cm,分别为裸露地表粗糙度的69.3、56和62.7倍。三角形沙障的地表粗糙度比正方形和六边形沙障的地表粗糙度分别高23.7%和10.5%。由于近地表湍流强度较高,0.2 m高度的风速变化很大,导致每组由1和0.2 m计算出的地表粗糙度有较大的离散程度。
-
传统的格状沙障平面图形都为正方形,原因是易于施工,铺设效率相对较高。而从几何角度上看,能够铺满平面的形状还有正六边形,正三角形,此外还有各种尺寸的密铺五边形。以往的研究表明,沙障的材料、网格尺寸及其形态对沙障的防风效能有显著影响,但都未涉及网格形状对沙障防风效能影响的研究。本文对相同材料及相同网格面积但网格形状不同的沙障进行野外试验,结果表明,沙障铺设平面几何形状对沙障的防风降沙能力有很大影响。
1) 对于单个网格面积都为1 m2的三角形,正方形和六边形草沙障,高度为0.2 m时,正方形沙障的防风效能比三角形和六边形沙障分别低12%和8%,高度为0.3 m时,分别低11%和10%,高度为0.5 m时,分别低9%和4%。高度为1 m时,正方形沙障防风效能与三角形沙障相同,显著高于六边形沙障。高度小于0.3 m时,三角形沙障和六边形沙障防风效能相近,差异不显著。近地表的防风效能,正方形沙障最差,三角形与六边形沙障防风效果相近,均好于四边形沙障。
2) 当地表铺设沙障后,沙障内部风速沿高度变化更接近指数函数函数。在0.2~0.5 m高度范围内,三角形沙障和六边形沙障表现出相近的变化特征,且各高度风速分布显著小于正方形沙障风速。高度为0.2~0.3 m时,三种沙障的风速增长较为迅速,增长速度大致相同。高度为0.5~1 m时,风速增长速度具有显著的差异。风速增长速度为正方形网格沙障 < 三角形网格沙障 < 六边形网格沙障,但都大于无沙障时的风速增长速度。
3) 风速小于6 m/s时,地表粗糙度随风速变化接近线性变化规律。裸露地表 1 m高处风速为4 m/s时,裸露地表粗糙度为0.166 cm三角形沙障、正方形沙障和六边形沙障的地表粗糙度分别为11.5、9.3和10.4 cm,分别为裸露地表粗糙度的69.3、56和62.7倍。三角形沙障的地表粗糙度比正方形沙障和六边形沙障的地表粗糙度分别高23.7%和10.5%。
Research on the windproof efficiency of polygonal straw sand barrier
-
摘要: 为了解正三角形、正方形和正六边形格状草沙障的防风固沙能力,通过测定沙障样地内不同高度风速,针对单个网格面积为1 m2的3种沙障,对沙障的防风效能、风速廓线特征和地表粗糙度进行研究。结果表明,正三角形和正六边形沙障的防风效能显著大于正方形沙障。高度为0.2 m时,正方形沙障的防风效能比三角形和六边形沙障分别低12%和8%,高度为0.3 m时,分别低11%和10%。0.2~0.3 m高度范围内,三角形沙障和六边形沙障的防风效能差异不显著。风速小于6 m/s时,地表粗糙度随风速变化接近线性变化规律。裸露地表 1 m高处,风速为4 m/s时,裸露地表粗糙度为0.166 cm,三角形沙障、正方形沙障和六边形沙障的地表粗糙度分别为11.5、9.3和10.4 cm,分别为裸露地表粗糙度的69.3、56和62.7倍。三角形沙障的地表粗糙度比正方形沙障和六边形沙障的地表粗糙度分别高23.7%和10.5%。铺设三角形沙障的地表粗糙度最大。Abstract: To understand the windproof benefit of triangular sand barrier, square sand barrier and hexagonal sand barrier, we studied respectively the windproof efficiency, wind speed profile characteristics and roughness of ground surface of straw sand barrier, whose single mesh grid was 1 m2 by measuring the wind speeds in different heights in three straw sand barrier types. Results showed that wind speed efficiency could be reduced by triangular sand barrier and hexagonal sand barrier significantly than square sand barrier. The windproof efficiency of square sand barrier was 12% and 8% lower than triangular sand barrier and hexagonal sand barrier respectively at 0.2 m height, and was 11% and 10% lower than triangular sand barrier and hexagonal sand barrier respectively at 0.3 m height. The difference of windproof efficiency in triangular and hexagonal sand barrier was not obvious at he height between 0.2-0.3 m. The variation of roughness with wind speed was linear when the wind speed was less than 6 m/s. When wind speed was 4 m/s at 1 m height above bare sand, the roughness of bare sand was 0.166 cm, and the roughness of triangular sand barrier, square sand barrier and hexagonal sand barrier was 11.5, 9.3 and 10.4 cm, respectively, and was 69.3, 56 and 62.7 times of bare sand. The roughness of triangular sand barrier was 23.7% and 10.5% higher than that of square sand barrier and hexagonal sand barrier, respectively. Consequently, the roughness of triangular straw sand barrier is the highest.
-
Key words:
- straw sand barrier /
- windproof efficiency /
- roughness /
- wind speed profile
-
表 1 多边形沙障防风效益对比
Table 1. Comparison in windproof efficiency among polygon sand barrier types
% 测量高度
Measuring height/m沙障类型Sand barrier type 三角形
Triangular正方形
Square六边形
Hexagonal0.2 67±16a 55±15b 63±13c 0.3 46±13a 35±16b 45±12a 0.5 25±9a 16±11b 30±8c 1 7±5a 7±3a 3±4b 注:不同字母表示不同类型沙障防风效能差异显著(P < 0.05)。Note: different letters mean significant difference in windproof efficiency among three types of sand barrier at P<0.05 level. -
[1] 庞营军, 屈建军, 谢胜波, 等.高立式格状沙障防风效益[J].水土保持通报, 2014, 34(5):11-14. http://d.old.wanfangdata.com.cn/Periodical/stbctb201405003 PANG Y J, QU J J, XIE S B, et al. Windproof efficiency of upright checkerboard sand-barriers[J]. Bulletin of Soil and Water Conservation, 2014, 34(5):11-14. http://d.old.wanfangdata.com.cn/Periodical/stbctb201405003 [2] 张登山, 吴汪洋, 田慧丽, 等.青海湖沙地麦草方格沙障的蚀积效应与规格选取[J].地理科学, 2014, 34(5):627-634. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dlkx201405016 ZHANG D S, WU W Y, TIAN H L, et al.Effects of erosion and deposition and dimensions selection of straw-checkerboard barriers in the desert of Qinghai Lake[J]. Scientia Geographica Sinica, 2014, 34(5):627-634. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dlkx201405016 [3] 屈建军, 凌裕泉, 俎瑞平, 等.半隐蔽格状沙障的综合防护效益观测研究[J].中国沙漠, 2005, 25(3):329-335. doi: 10.3321/j.issn:1000-694X.2005.03.005 QU J J, LING Y Q, ZU R P, et al.Study on comprehensive sand-protecting efficiency of semi-buried checkerboard sand-barriers[J]. Journal of Desert Research, 2005, 25(3):329-335. doi: 10.3321/j.issn:1000-694X.2005.03.005 [4] 黄强, 雷加强, 王雪芹.塔里木沙漠公路不同地貌部位的高立式沙障阻沙特征[J].干旱区地理, 2000, 23(3): 227-232. doi: 10.3321/j.issn:1000-6060.2000.03.006 HUANG Q, LEI J Q, WANG X Q. Sand-obstructing effect of the high sand barriers in the different morhpologic sections along the Traim Desert Highway[J]. Arid Land Geography, 2000, 23(3): 227-232. doi: 10.3321/j.issn:1000-6060.2000.03.006 [5] 马学喜, 李生宇, 王海峰, 等.固沙网沙障积沙凹曲面特征及其固沙效益分析[J].干旱区研究, 2016, 33(4):898-903. http://d.old.wanfangdata.com.cn/Periodical/ghqyj201604029 MA X X, LI S Y, WANG H F, et al.Concave surface features in sand-fixing net barriers and evaluation of sand-fixing benefits[J]. Arid Zone Research, 2016, 33(4):898-903. http://d.old.wanfangdata.com.cn/Periodical/ghqyj201604029 [6] 党晓宏, 高永, 虞毅, 等.新型生物可降解PLA沙障与传统草方格沙障防风效益[J].北京林业大学学报, 2015, 37(3):118-125. doi: 10.13332/j.1000-1522.20140245 DANG X H, GAO Y, YU Y, et al. Windproof efficiency with new biodegradable PLA sand barrier and traditional straw sand barrier[J]. Journal of Beijing Forestry University, 2015, 37(3):118-125. doi: 10.13332/j.1000-1522.20140245 [7] 蒙仲举, 任晓萌, 高永.半隐蔽式沙柳沙障的防风阻沙效益[J].水土保持通报, 2014, 34(3):178-206. http://d.old.wanfangdata.com.cn/Periodical/stbctb201403034 MENG Z J, REN X M, GAO Y. Effect of semi-buried Salix psammophila checkerboad on wind-preventing and sand-resisting[J]. Bulletin of Soil and Water Conservation, 2014, 34(3):178-206. http://d.old.wanfangdata.com.cn/Periodical/stbctb201403034 [8] 刘虎俊, 王继和, 李毅, 等.塑料网方格沙障对新月形沙丘迎风坡的风沙流影响[J].水土保持学报, 2011, 25(5):26-34. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201105006 LIU H J, WANG J H, LI Y, et al. Effects of plastic checkerboard sand-barrier on wind-sand flux of leeward of crescentic dune[J]. Journal of Soil and Water Conservation, 2011, 25(5):26-34. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201105006 [9] 张克存, 俎瑞平, 方海燕.不同孔隙度尼龙网对风沙流减弱作用的风洞模拟[J].水土保持学报, 2004, 18(4):5-11. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb200404002 ZHANG K C, ZU R P, FANG H Y. Simulation on abated effect of nylon net with different porosities on wind-blown sand in wind tunnel[J]. Journal of Soil and Water Conservation, 2004, 18(4):5-11. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb200404002 [10] 屈建军, 刘万贤, 雷家强, 等.尼龙网栅栏防沙效应的风洞模拟实验[J].中国沙漠, 2001, 21(3):62-66. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsm200103013 QU J J, LIU W X, LEI J Q, et al. Simulation experiments on sand-arresting effect of nylon net fence in wind tunnel[J]. Journal of Desert Research, 2001, 21(3):62-66. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsm200103013 [11] 张克存, 屈建军, 牛清河, 等.青藏铁路沿线砾石方格固沙机理风洞模拟研究[J].地球科学进展, 2010, 25(3):284-289. http://d.old.wanfangdata.com.cn/Periodical/dqkxjz201003006 ZHANG K C, QU J J, NIU Q H, et al. Simulative research on the mechanism of rocky checker board sand barriers along Qinghai-Tibet rail way in wind tunnel[J]. Advances in Earth Science, 2010, 25 (3):284-289. http://d.old.wanfangdata.com.cn/Periodical/dqkxjz201003006 [12] BO T L, MA P, ZHENG X J. Numerical study on the effect of semi-buried straw checkerboard sand barriers belt on the wind speed[J]. Aeolian Research, 2015, 16:101-107. doi: 10.1016/j.aeolia.2014.10.002 [13] 周娜, 张春来, 田金鹭, 等.半隐蔽式草方格沙障凹曲面形成的流场解析及沉积表征[J].地理研究, 2014, 33(11):2145-2156. http://d.old.wanfangdata.com.cn/Periodical/dlyj201411014 ZHOU N, ZHANG C L, TIAN J L, et al. Flow field controlling the concave surface of the semi-buried checkerboards and its characterization by grain sizes of sediments[J]. Geographical Research, 2014, 33(11):2145-2156. http://d.old.wanfangdata.com.cn/Periodical/dlyj201411014 [14] 李锦荣, 孙保平, 高永, 等.基于空气动力学的沙袋沙障气流场模拟[J].北京理工大学学报, 2010, 30(6):749-752. http://d.old.wanfangdata.com.cn/Periodical/bjlgdxxb201006027 LI J R, SUN B P, GAO Y, et al. Simulation of sand-bag barrier wind flow field based on aero dynamics[J]. Transactions of Beijing Institute of Technology, 2010, 30(6):749-752. http://d.old.wanfangdata.com.cn/Periodical/bjlgdxxb201006027 [15] CORNELIS W M, GABRIELS D. Optimal windbreak design for wind-erosion control[J]. Journal of Arid Environments, 2005, 61(2):315-332. doi: 10.1016/j.jaridenv.2004.10.005 [16] SALEH A, FRYREAR D W. Soil roughness for the revised wind erosion equation[J]. Journal of Soil and Water Conservation, 1999, 54(2):473-476. [17] 鞠英芹, 刘寿东, 马德粟, 等.农田与草地下垫面空气动力学粗糙度的研究[J].科技通报, 2016, 32(8):5-16. doi: 10.3969/j.issn.1001-7119.2016.08.002 JU Y Q, LIU S D, MA D S, et al. Study on aerodynamic roughness length on crop and grass underlying surfaces[J]. Bulletin of Science and Technology, 2016, 32(8):5-16. doi: 10.3969/j.issn.1001-7119.2016.08.002 [18] CHEN K, ZHU F R, NIU Z N, et al. Evaluation on shelter effect of porous windbreak fence through wind tunnel test[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2006, 42(5):636-640. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bjdxxb200605013 [19] DONG Z B, LUO W Y, QIAN G Q, et al. A wind tunnel simulation of the mean velocity fields behind upright porous fences[J]. Agricultural and Forest Meteorology, 2007, 146(1):82-93. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=508be126cf115e08e0c8508898508892 [20] DONG Z B, LUO W Y, QIAN G Q, et al. Evaluating the optimal porosity of fences for reducing wind erosion[J]. Sciences in Cold Arid Regions, 2011, 3(1):1-12. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hhqkx-e201101001 [21] ZHANG C L, LI Q, ZHOU N, et al. Field observations of wind profiles and sand fluxes above the windward slope of a sand dune before and after the establishment of semi-buried straw checkerboard barriers[J]. Aeolian Research, 2016, 20:59-70. doi: 10.1016/j.aeolia.2015.11.003 -