塔克拉玛干沙漠南缘策勒绿洲干季输沙通量的垂向分布.docxVIP

塔克拉玛干沙漠南缘策勒绿洲干季输沙通量的垂向分布 沙流是风和气的两相流。它是沙子物理的核心内容，也是沙地形、沙漠化和防沙工程的基本理论之一。风沙流中输沙通量的垂向分布是沙粒微观运动的宏观表现,对于正确估计输沙率、验证已有的数值模型具有重要意义。目前国内外对输沙通量的垂向分布有定性和定量两种方式的描述。定性描述方面,学者们基于地表处于吹蚀和堆积不同过程而具有的特征,提出结构数及特征值以区别这两个过程,并在风沙环境研究中得到了广泛应用。定量描述方面,研究者大多采用经验函数拟合[8,9,10,11,12,13,14,15,16,17,18,19,20,21],从现有的结论来看,描述模型主要为简单的指数函数或幂函数[9,10,11,12,13,14]、修正的指数函数或幂函数,幂函数加指数函数的复合函数描述,或者利用指数函数和幂函数进行分段描述。定量描述的研究多集中在对模型预测精度的探讨,以使其能够更为准确地描述输沙通量随高度的分布,这在防沙治沙的实践中非常必要,因为研究者常需要通过积分高度-输沙率关系式的方法来估算整个截面的输沙率。另一方面,为了研究风沙流结构的影响机制,模型本身的参数应该反映丰富的风沙流信息,尽管学者们对此越来越关注,然而到目前为止,关于模型的预测精度与模型参数包涵的风沙流信息程度,尚缺乏系统的比较研究。此外,风洞实验的研究结果表明,对于均匀沙和非均匀沙,输沙通量的垂向分布遵循不同的规律;在海岸带,自然混合沙的观测结果表现出类似于风洞中非均匀沙的垂向分布规律,但这种现象很少被内陆干旱区的风沙工作者所关注。 塔克拉玛干沙漠南缘深居亚洲内陆,干燥度大,风沙频繁,沙粒粒径极细,地表物质易于流动,特别是频繁发生于沙漠———绿洲过渡带的风沙活动,对绿洲安全构成极大的威胁,但关于该区域风沙流结构的系统研究却罕见报道。本研究通过野外实地观测,研究该区域风沙流中输沙通量的垂直分布,系统比较4种常见的输沙通量垂向分布模型在拟合精度和反映风沙流信息程度两方面的差异,并初步探讨在塔克拉玛干沙漠南缘特定物质组成情况下输沙通量垂向分布的适用模型。 1 研究区域和方法 1.1 地表沙粒粒度分析 野外观测在塔克拉玛干沙漠南缘的策勒绿洲西侧外缘进行(图1A),具体位置为37°00′59″N和80°43′16″E,距离中国科学院策勒荒漠草地生态系统国家野外研究站约2km,地表为相对平坦的裸露沙地。根据2009年对策勒绿洲西部外围过渡带的沙粒粒度分析结果,沙物质组成基本以极细沙为主,其中沙漠化程度最为严重的样区粉沙、极细沙和细沙的含量依次为24.4%、59.1%和14.7%。该区域年均降水量35mm,主要集中在5—7月,年均蒸发量2 595mm,干燥度20.8,属极端干旱内陆荒漠气候。起沙风主要集中在春夏两季,以W风为主,WNW次之,合成输沙势为27.98VU,合成输沙方向为92.10°,风向变率指数为0.94,为单风向低风能环境。 1.2 输沙通量垂向分布的观测系统设计 于2010年4—6月在选定地点进行风速、风向与输沙量同步观测,观测期间地表沙粒呈自然干燥状态。风速观测采用长春气象仪器研究所研制的DETI可移动测风系统(图1B),风速风向传感器的设置高度均为2m,风速传感器的分辨率为0.1m·s-1,测量范围为0.8~30m·s-1,风向传感器的分辨率为5.625°,测量范围为0°~360°,数据采集间隔为1min。输沙通量观测采用20cm和40cm两个高度的阶梯式积沙仪,进沙口大小为2cm×2cm,这种积沙仪在风洞中测试积沙效率可以达到90%。较低风速下采用20cm高积沙仪,高风速下采用40cm高积沙仪,起沙时将积沙仪埋入沙地,使最底层的进沙口与地表齐平,进沙口正对气流方向,观测时将铁盖打开并记录观测时间,积沙时间为10~15min,同步风速风向值由测风系统自动记录和保存。每一次积沙过程利用一个积沙仪采集沙样,记为一组,共采集到28组,对采集到的沙样用0.01g电子天平称重。 选取4种最为常用的描述输沙通量垂向分布的模型(表1),利用野外观测数据,通过代数拟合程度和拟合参数反映的风沙流信息来评价这4种模型。在代数拟合程度方面,利用Akaike Information Criterion(AIC)标准进行评价,该标准同时考虑了拟合方程的逼近效果和参数个数两方面的因素被学者所认可,其计算方法为: 式中:n为积沙仪层数;qi0为各层实测值;qi为计算的理论值;p为方程中参数个数;AIC值小的方程为较理想的拟合方程。在拟合参数所能反映的风沙流信息方面,利用逐步线性回归,建立拟合参数与输沙率、平均风速和平均跃移高度等的关系,若相关程度越高,则其包涵的风沙流的信息越丰富。综合以上两个方面的评价结果,确定能够描述输沙通量垂向分布的最佳方程。方程拟合在sigm