Underestimation of roughness in rough rock joints

Numerous studies have been made to improve Barton's shear strength model for the quantification of rock joints. However, in these previous studies, the roughness and shear strength of the rock joint have been underestimated especially for relatively high undulated profiles (joint roughness coefficient (JRC) >14). The main factors of roughness underestimation in rough rock joints are investigated for the proper quantification of rock joint roughness. The aliasing effect and the roughness characteristics are analyzed by using artificial joint profiles and natural rock joint profiles. A 3D camera scanner is adopted to verify the main source of underestimation when using conventional measurement methods. Shear strength tests are carried out by using two types of shear apparatus to study the roughness mobilization characteristics, which may also affect the roughness underestimation. The results of joint roughness assessment, such as aliasing and undulation of waviness, show that the roughness can be underestimated in relatively rough joint profiles (JRC>14). At least two components of roughness parameters are needed to properly represent the joint roughness, for example, the amplitude and the inclination angle of joint asperity. Roughness mobilization is affected by both the normal stress and the asperity scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Study of the relationship between non-dimensional roughness length and wave age, effected by wave directionality

Relationship between the non-dimensional roughness length and inverse of wave age has been discussed without consideration of wave directions, though wind wave field consists of various directional component waves. In this study we observe wave heights by an array of four wave gauges at the Hiratsuka Tower of (Independent Administrative Institution) National Research Institute for Earth Science and Disaster Prevention (NIED), Japan, and discuss the effect of wave directionality. As a result, the data sets were classified into two different groups according to the directional wave spectrum distribution. In case 1 only swell and wind waves exist and in case 2 there exist wave components from several directions. It is shown that the case of multiple-directional component waves (case 2) may affect the non-dimensional roughness length and friction velocity.     

有关粉尘释放模型的应力分配模式存在问题的讨论

准确模拟粉尘释放通量对评价粉尘气溶胶的气候效应、沙漠化和沙尘暴防治具有重要意义。敏感性试验、室内和野外风洞实验表明被广泛使用的DPM（Dust Production Model）粉尘释放模型的2个应力分配模式——Marticorena模式和 Alfaro模式存在明显问题：由于目前空气动力学粗糙度模式还不能考虑粗糙元间隙率和粗糙元分布状况空气动力学粗糙度的影响,以空气动力学粗糙度为主要参数的Marticorena模式和 Alfaro模式还不能考虑粗糙元间隙率和粗糙元分布状况对应力分配的影响,而相关实验表明这些因素对应力分配存在显著影响;Marticorena模式和 Alfaro模式预测的应力分配系数和起动摩阻风速存在显著差异,目前实验结果还不能判断二者正确与否;Marticorena模式和 Alfaro模式预测的能够有效控制风蚀的应空气动力学粗糙度与实际观测相矛盾。针对上述问题提出了通过风洞实验改进当前应力分配模式的主要途径。     

一种描述结构面剖面线粗糙度的新方法

采用法向矢量单位圆描述结构面剖面线粗糙度,从各微分段的角度关系阐述粗糙度,进而将二维问题转化为一维问题处理并提出"角度粗糙度"的概念.考虑到各微分段的实际长度对粗糙度的贡献,采用加权均值与加权方差定量描述角度粗糙度;角度粗糙度越大表明该剖面越粗糙.对规则剖面线与不规则剖面线采用"角度粗糙度"进行描述,所得结果跟已有的剖...     

Atmospheric surface layer characteristics with MONTBLEX data

Data obtained from the 30m high MONTBLEX tower installed at the Indian Institute of Technology, Kharagpur are described. Data on wind speed and direction, temperature and humidity were recorded during the first week of July 1989 at six levels on the tower. They indicate some of the atmospheric surface layer characteristics. Using two levels of tower data involving wind speed and temperature a profile method was used for computing the surface fluxes of heat and momentum. A method for calculating the surface roughness length using the data was also used to obtain the vertical windspeed profile. The computed wind profile was compared with observations, while the surface roughness length was compared with values obtained by a least square fit. The computed surface fluxes were compared with theoretical values obtained by an energy budget method which uses only wind data at one level. The agreement was reasonably good.     

Roughness length and drag coefficient at two MONTBLEX-90 tower stations

Using MONTBLEX-90 mean velocity data, roughness lengths and drag coefficients are estimated at Jodhpur and Kharagpur. At Jodhpur, since the surface is not uniform the roughness length is estimated separately in three different subsectors within the range of prevailing wind directions and averages to 1.23 cm in the sector between 200° and 230° which is relatively flat with no obstacles on the ground. At Kharagpur, where the terrain is more nearly homogeneous, the average value (for all prevailing wind directions) is 1.94 cm. The drag coefficient C_D at Jodhpur shows variation both with the roughness subsector and with wind speed, the average over all directions increasing rapidly as themean wind speed Ū₁₀ at 10m height drops according to the power lawC _D = 0.05 Ū ₁₀ ^t-1.09 in trie range 0.5 < Ū₁₀ < 7 m s⁻¹. At Kharagpur, the drag coefficient is smaller than at Jodhpur by nearly 50% for the same range of wind speeds (> 3 ms⁻¹).     

月表虹湾Hurst指数分析

月球遥感是在登月技术还不完善的情况下,利用遥感数据,通过各种分析来研究月球演化、内部结果和表面特征的方法。粗糙度是地表隆升、沉降、侵蚀和火山喷射物覆盖在行星表面留下的记录,因此对行星表面粗糙度的定量分析可以精确反映其地质演化过程。而Hurst指数则是粗糙度的精确反映。文中利用LOLA制作的DEM,以月球正面虹湾为例,选取水平方向的12条剖面和垂直方向的15条剖面对其各粗糙度参数进行计算,并结合月表成熟度对Hurst指数与底层年龄的相关性进行了分析,得出如下结论:(1)在1km的剖面长度上,虹湾地区均方根高程的平均值在3m左右;在0.2～3km研究尺度范围内,虹湾地区的坡度不超过2°,说明该地表起伏度小。(2)虹湾地区的Hurst指数跨度比较大,尤其在垂直方向(南北)上,在0.4到0.9之间,表明其构造样式丰富。(3)在0.2～3km研究尺度范围内,Hurst指数自北往南增大、自西往东增大,表示在虹湾地区粗糙度自西北向东南方向增大。(4)虹湾地区地质年龄为自西北向东南方向变小,与Hurst指数增大方向一致。(5)在月表年轻撞击坑地区月表成熟度与Hurst指数大致呈负相关关系。     

汶川地震地表破裂面形貌特征

准确描述破裂面形貌对于我们理解地震断层作用是非常重要的,破裂面的形貌特征包含许多关于地震和断层机制的有用信息。在Mw7.9 2008汶川地震中断层活动产生了两个新鲜的破裂面,八角庙破裂面和沙坝破裂面。我们使用3D便携式激光扫描仪(Tri mble GX)对两个破裂表面进行测量,在野外微观尺度上研究了破裂面形貌特征。通过能谱密度和均方值两个方法分析破裂面形貌,新鲜的破裂面表现为自相仿性,能谱密度和均方值均与剖面长度存在幂律关系。在能谱密度与空间频率的对数图中,能谱密度曲线存在明显的拐点,该拐点所对应的波长称为"特征波长",表明单一分形不能准确描述破裂面形貌。八角庙破裂面在平行滑动方向上的特征波长为7 mm,在垂直方向上特征波长略大一些(区域Ⅰ为10 mm,区域Ⅱ为9 mm);沙坝破裂面在平行滑动方向上的特征波长为8 mm,但垂直方向上特征波长略小(6 mm)。均方值曲线的最小二乘拟合直线的斜率为Hurst指数,该指数依赖于剖面线方向并描述破裂面形貌的各向异性,H指数的最小值和最大值分别与平行擦痕和垂直擦痕方向对应,这与野外断层面擦痕测量结果一致。沙坝破裂面的H指数极坐标图中存在次级H指数峰值(对应剖面线方向为85°和160°),这揭示破裂面上存在一组隐匿擦痕。该组隐匿擦痕为汶川地震之前断层活动中形成的,但这还不足以推测上一次断层活动的时间和规模。另外,通过比较新鲜节理面和破裂面表明H指数是否大于0.8反映了断层类型。在整个空间频率域上,能谱密度曲线斜率(-α)和均方值曲线斜率(H)的线性拟合关系为α=1.22+1.72H,并不严格满足两者之间理论关系式,α=1+2H。这个差异是由于测量信号噪音、破裂面的多分形性和分析方法的差异造成的。     

Wind tunnel experiments on the turbulent transmission over the near surface layer of different surfaces

To understand the generation and development mechanism of soil erosion, data of the wind profile from wind tunnel experiments were used, and the characteristics of the turbulent transmission with free stream velocity and the properties of the surfaces in arid and semi-arid regions were studied by the atmospheric aerodynamic method. The result showed that the turbulent momentum fluxes are proportional to the velocity and the roughness of the surfaces. In addition, they are closely related to the density of the roughness elements; the turbulent exchange coefficient over the near surface layer increases with height.     

流动型态对曼宁糙率系数的影响研究

根据明渠水流流动型态的概念,通过对实验资料的分析,得到均匀流时糙率系数随水深增加而减小,形成M1型水面线的非均匀流时糙率系数随水深、水力坡度的增大而增大的规律。提出的糙率系数的二步计算法,可以解决因非均匀流水深沿程变化,难以建立糙率系数与水深关系的困难,为在水力计算中修正糙率系数提供了依据。     

沙粒当量粗糙度的分形表达

沙粒当量粗糙度k_s的准确确定是计算沙粒阻力的关键。针对沙粒轮廓剖面具有自相似精细结构的特点,提出采用分形维数对其复杂性及不规则性进行定量描述;通过粘沙平整床面、初始静平整床面和清水冲刷粗化床面的沙粒阻力试验研究k_s的分形表达方法,建立了静平整床面时k_s的分形表达式,并对公式进行了验证、比较和讨论;结果表明,公式引入的中值粒径d₅₀、分形维数D和弗劳德数Fr 3个参数,可定量地综合反映床沙颗粒大小、级配组成、随机排列以及水流条件对k_s的影响。     

基于功率谱密度的评估岩石节理粗糙度新方法

为进一步提高结构面粗糙特征描述的精确性和可靠性,更全面地体现出细微形貌的影响,将轮廓线凹凸起伏分解为幅值、倾角和曲率3个特征,分别代表轮廓线的高度变化、角度变化和弯曲程度变化。以Barton建议的10条标准轮廓线为例,计算其幅度、角度和弯曲度的特征曲线,获得了三者的分布规律。利用可表征微细观形貌的功率谱密度函数研究了3种特征曲线的频谱特征,然后从功率谱密度对数曲线中提取了可以描述形貌特征的分形参数D_j（斜率相关）和粗糙度参数A_j（截距相关）,在此基础上研究了D_j和A_j与节理粗糙度系数（joint roughness coefficient,简称JRC）的关系。同时,将该指标推广到了三维形貌的表征形式,并计算了大理岩剪切面的三维粗糙度。结果表明：新指标评估精确度高,可反映结构面各向异性特征,能有效表征结构面粗糙特征,对研究结构面剪切力学性质有重要意义。     

大型浅水湖泊水动力模型不确定性和敏感性分析

选取国内外常用的水动力学模型(EFDC)和典型的浅水湖泊(太湖),采用拉丁超立方取样(LHS),研究湖泊水动力模块中常用的5个重要参数(风拖曳系数、床面粗糙高度、涡粘性系数、紊流扩散系数以及风遮挡系数)对湖体水位和流速的影响。结果表明:针对大型浅水湖泊,湖泊岸线形状和湖底地形、湖泊周围地形、湖泊水面风场对模拟结果产生决定性影响。尤其是在湖湾区和周边地形比较复杂的地区,风场参数对水动力模拟结果不确定性的贡献率最大。在垂向上,表层流速受到参数不确定性的影响最大,底层次之,中层最小。床面粗糙高度对水动力模拟结果不确定性贡献率较风场参数要小,水体涡粘滞系数和扩散系数影响则更小。故在选择大型浅水湖泊水动力模型参数时,要充分考虑湖泊岸线和周围地形,着重率定风场参数以及床面粗糙高度。     

近底层潮流速分布模式比较

考虑到潮流运动的非恒定性和非线性特征,基于潮流雷诺应力分布式,应用Prandtl混合长度和von Karman自相似理论,并结合潮流运动方程,推导了近底层潮流速垂向分布的双对数模式.然后,将该分布模式应用于英国大陆南部西Solent水道和江苏岸外西洋海域流速的实测资料,并与对数线性分布的计算成果相比较,获得了比较满意的结果.最后,对几种分布模式进行了讨论和比较,结果表明:在近岸水域采用对数或对数线性流速分布模式计算摩阻流速和粗糙高度往往会带来较大的偏差,双对数模式能使计算的结果更为合理.     

基于植被分布的河道糙率分区及率定方法

天然河流的河道综合糙率呈现出空间上的差异性和随水位（或流量）变化的动态性,但目前缺乏相关参数化方法来定量描述河道糙率的动态变化规律。尝试通过参数化方法开展受河道植被影响显著河流的糙率反演研究,用以提升模型精度。基于植被分布将河道断面划分为若干糙率不同的子区,通过率定河道断面各分区的糙率,从而反演糙率—水位曲线。在此基础上通过分析河道植被覆盖情况与河道断面特点对糙率曲线变化的影响,推求了基于分区糙率的河道断面综合糙率计算公式,从而定量描述分区糙率与综合糙率的关系。以漓江干流为例,采用该方法率定漓江干流（大溶江至阳朔段）水动力模型。结果表明:漓江干流综合糙率随水位在0.022~0.180间变化;在1.5 m的临界水深下,断面可划分为底床植被区（n=0.210）与非植被区（n=0.006）,能较好地反演糙率—水位曲线并获得理想的水位模拟效果。漓江底床植被繁茂是糙率随水位变化的根本原因,断面边滩的坡度变化是糙率与水位曲线梯度变化的主要驱动因素,两者的共同作用使得糙率随着水位呈现两段式的非线性变化。     

地表粗糙度参数化研究综述

粗糙度反映了地表的起伏程度,是许多陆面过程的关键影响因子。然而,人们对地表系统认识的不足,造成现有的各种地表粗糙度参数化方案均存在一定的问题。从地表测量技术、粗糙度相关参数和遥感研究3个方面对地表粗糙度参数化研究现状进行了综述。针式廓线法和激光廓线法是当前主流的地表测量技术,而三雏激光扫描和摄影测量技术已展示出了较大的潜力。基于统计方法和基于分形理论的粗糙度相关参数具有截然不同的物理意义,但地表复杂的多尺度特性使得很难用一类参数进行描述。光学遥感与微波遥感均具有广阔的前景,其中前者需要注意与经典粗糙度参数化方案相结合,后者则需拓展在下一代遥感平台上的技术与方法。还针对不同尺度粗糙度参数的比较与转换,地表粗糙度的空间异质性与各向异性,以及三维表面粗糙度参数化等当前地表粗糙度参数化研究中的一些关键问题进行了讨论。     

岩石结构面粗糙度系数尺寸效应的推拉试验研究

粗糙度系数是结构面抗剪强度的主要影响因素,然而由于结构面表面形态的复杂性,粗糙度系数尺寸效应研究并未获得较大进展。总结了结构面粗糙度系数的3种获取手段：标准剖面对比法、理论公式法、试验反分析法。在此基础上分析了3种方法在研究粗糙度系数尺寸效应方面存在的问题和困难。为了研究结构面粗糙度系数与试样尺寸的相关度,对中砂、硅粉、水泥、非引气型萘系减水剂等原材料的配比进行了研究,获得了与天然钙质板岩物理力学特性相类似的岩石模型材料,然后采用研发的结构面制作模具及其制备工艺制作了8组共176对具有不同尺寸和表面起伏粗糙程度的结构面,并利用改进的高精度岩石结构面推拉仪对结构面粗糙度系数进行了推拉试验研究和数据统计分析,结果表明：模型结构面粗糙度系数的统计均值随试样尺寸的增加而降低,但特定结构面粗糙度系数的尺寸效应规律需要根据结构面的具体表面形貌进行测试;Barton理论公式计算的结构面粗糙度系数尺寸效应变化规律与推拉试验测试规律总体上一致,但试验值与理论值有差异,且结构面试样尺寸越小,二者的差异就越大;具有特定表面形貌的模型结构面粗糙度系数也有差异,工程大尺寸岩体结构面粗糙度系数需要根据表面形貌和分布特征进行综合判定。     

微凸体对抛光岩石表面静摩擦系数影响试验研究

对典型岩石摩擦滑动试验装置进行了改进,以8种硅酸盐岩作为研究对象,在低正应力条件下对岩石摩擦滑动过程中的静摩擦系数进行了试验研究。对滑动表面形貌进行了表征,从统计学角度分析了粗糙度对静摩擦系数的影响,并基于滑动表面微凸体（凹凸体、微观粗糙度）对这一影响作进一步分析。研究表明：滑动面为抛光面时,千枚岩、石英岩、岩屑砂岩和石英砂岩的静摩擦系数为0.38～0.47,砾岩、含砾粗粒石英砂岩、中粒岩屑砂岩和中粒石英砂岩的静摩擦系数为0.83～1.07;在粗糙度中,轮廓最大谷深Rm和轮廓最大峰高Rp两个参数导致摩擦滑动中产生不同静摩擦系数,且随着这两个参数的增加,静摩擦系数呈指数规律增加;千枚岩、石英岩、岩屑砂岩和石英砂岩滑动表面的微凸体数量较少,砾岩、含砾粗粒石英砂岩、中粒岩屑砂岩和中粒石英砂岩滑动表面的微凸体数量较多,随着微凸体数量的增加岩石摩擦滑动过程中的静摩擦系数增加。     

Development of deflation lag surfaces

A series of wind tunnel tests were carried out to investigate the development of deflation lags in relation to the non-erodible roughness element concentration. Glass spheres (18 mm in diameter) were placed along the complete length of the wind tunnel working section in regular staggered arrays using three different spacings (d=18, 30 and 60 mm) and completely covered with a 0.27-mm erodible sand. A pre-selected free stream velocity above threshold (8m s^?1) was established above the surface and the sediment transport measured at 2-s intervals using a wedge-shaped trap in which an electronic balance is incorporated. Throughout each test, the emerging lag surface was periodically photographed from above at two locations upwind of the trap. The photographs were electronically scanned and analysed to calculate the lag element coverage and location, as well as mean height and frontal area for each time period. Test results indicate that lag development has a profound effect on both the sediment flux and wind profile characteristics. Initially, there is an increase in sediment flux above that for a rippled sand bed because of increased erosion around and reduced kinetic energy loss in highly elastic collisions with the emerging roughness elements. With further emergence, a dynamic threshold is reached whereupon the sediment flux decreases rapidly, tending towards zero. At this point, the supply of grains to the air stream through fluid drag follows a reduction in aerodynamic roughness and, therefore, surface shearing stress. At least as important is the lesser potential for grain ejection at impact because of reduced momentum imparted from the air stream during saltation. Although recent shear stress partitioning models indicate when particle movement may commence on varying surfaces, our experimental results demonstrate that this partitioning has a further direct bearing upon the saltation flux ratio.     

随机介质模型的模拟与混合型随机介质

讨论了随机介质模型的基本概念及指数型和高斯型椭圆自相关函数所描述的随机介质模型的特点, 并提出了混合型随机介质模型的概念, 该随机介质模型能更加灵活、准确地描述实际介质.通过选择在水平方向和垂直方向上的自相关长度a、b以及粗糙度r, 可以产生出各种不同形式的混合型随机介质模型.模拟结果显示, 混合型随机介质模型能更加灵活地描述实际介质, 具有适应性强, 使用方便、灵活, 能有效地模拟油气藏细节的优点.