河道泥沙粒度非线性特征及其应用.docVIP

河道泥沙粒度非线性特征及其应用 王协康 方 铎 曹叔尤  (四川大学高速水力学国家重点实验室 成都 610065) 摘 要 在分析河道泥沙级配特性的基础上，利用分形理论定义其粒度特征。针对床沙粗化实验中依据推移质输沙量确定粗化稳定的不足，提出了以粗化层及推移质粒度分维值判别床面粗化稳定的方法；并建立了粗化床面糙度分维与床面糙率的关系式。通过实验资料计算验证，得到较满意的结果。 关键词 泥沙级配 分形理论 粗化稳定 床面糙度. 河道泥沙的级配曲线通常画在半对数坐标纸上，以横坐标(对数坐标)表示粒径，纵坐标(自然坐标)表示小于某粒径沙重占总沙重的百分数。天然河道泥沙级配与很多因素有关，它除了与流域特性、河道上游来水来沙情况以及河段内河床河岸性质有关外，还与沙样类别(悬移质、推移质或河床质)以及河床所处的冲淤状态(冲刷过程、淤积过程或相对平衡过程)有关。此外，泥沙取样的随机性、颗分仪器精度及操作人为性等因素都在一定程度上对泥沙级配分析有影响。实测资料表明，泥沙级配曲线在半对数纸上可表现为正态分布、正偏态分布、负偏态分布、双峰型分布等多种形态［１］。 河道泥沙级配是进行河道演变分析、模型设计、冲淤计算等问题研究中必不可少的基本参变量，具有重要的应用价值。众所周知，表征泥沙级配特性的主要参数，一个是反映沙样粗细的中值粒径；另一个是反映沙样组成的非均匀性参数。如人们按照不同的要求采用不同的值，如、、、等。此外，不少学者对悬移质泥沙级配曲线提出了不同的曲线拟合方式。如谢葆玲采用差分判别法提出了型级配曲线的函数关系式［２］； 熊治平利用双曲正切函数进行描述泥沙级配曲线［３］。但对于级配范围宽的粗颗粒床沙和推移质泥沙，研究其级配曲线拟合的较少。为此利用分形理论研究河道泥沙粒度分维及其应用。 1 分形理论及泥沙级配分维粒度 分形(Fractal)一词是由法国数学家曼德布罗特于1967年首创，用来描述那些表面看来杂乱无章、变幻莫测而实质上潜在有某种内在规律性的几何图形或形状。1960年，曼德布罗特发现棉花价格随时间的变化隐含着这样一种规律，即日变化规律与月变化规律非常相似；此后不久又发现计算机通讯用的电话线的噪音也有类似的规律。至1967年，他将对自然界的这种现象的思考集中反映到“英国的海岸线有多长?”一文中［４］，指出在没有参照物的情况下，从空中拍摄的100km长的海岸线与放大的10km的海岸线的照片看上去非常相似。此外自然界中起伏的山峦、变幻的浮云也都具有这种在不同尺度上具有相同程度的不规则性和复杂性即自相似性。曼德布罗特将这种具有自相似性且有特征长度的图形称为分形。描述分形的一个重要参数是分形维数，简称分维,其定义可表示为 (1) 式中，为图形的测度，为对应所得的测量数，为该图形的分维数。 河道泥沙级配是在水沙及边界条件下形成的混合物的一种量度，其粒度组成是一种没有特征长度的结构图形，其分形是建立在一种统计分布上，也就是说具有统计意义上的自相似性。设泥沙混合物颗粒的粒径为，粒径小于的颗粒数为，如果泥沙混合物粒度组成具有分形结构，对比分维定义式(1)，则有 (2) 设为粒径小于颗粒累积重量，为总重量，若，则   (3) 由(2)式得，而粒度成分增加是与重量(质量)增加相关，即  (4) 由此可知，而幂律分布也等价于分形分布。 实际上，为粒径小于的颗粒累积重量百分数，也是表示泥沙级配曲线的量，这样就为我们提供了一种泥沙级配粒度组成分维计算的方法。如果和在双对数坐标图中具有很好的线性关系，则表明粒度成分具有分形结构，拟合线段斜率值为，利用可求粒度组成的分维值。 表1 粗化层级配粒度分维计算 序 号 回归方程 值 相关系数 说 明 1－1 0.70 0.90 2.30 近似粗化稳定 1－2 0.75 0.91 2.25 2－1 0.74 0.96 2.26 2－2 1.00 0.93 2.00 粗化稳定 2－3 1.01 0.96 1.99 3－1 1.07 0.96 1.93 3－2 1.06 0.97 1.94 粗化稳定 4－1 1.09 0.95 1.91 4－2 1.08 0.97 1.92 粗化稳定 5－1 1.27 0.