王光谦－流域泥沙过程.ppt

流域泥沙过程模拟 流域泥沙过程的基本规律与建模方式 河网编码与计算机制 坡面产沙模型 河网汇流与输沙模型 重力侵蚀模型 一维河道水动力学模型 流域泥沙过程模拟的应用 1.流域泥沙过程的基本规律与建模方式 从峁顶至峁坡、沟坡和各级沟道，各部位存在极限含沙量，且渐次增大 由于重力侵蚀的滞后和高含沙量对挟沙力的增强作用，沙峰相对于洪峰有滞后和拖长现象，且随着流域面积增大和沟道级别的增加，现象更加显著 由于重力侵蚀的随机性和沟道内的冲淤变化，含沙量与流量的关系极为散乱 单场洪水中沟道的冲淤变化使含沙量低时颗粒组成细，含沙量高时颗粒组成粗，且这一差异随沟道等级的增大而愈明显 2.河网编码与计算机制 将超一维的河网抽象为二叉树 直接定位能力: 第n层的起始编码为2n-1 节点k的子节点编号为2k和2k+1 与河网的对应关系： 节点k—河段 节点2k—上游干流河段 节点2k+1—上游支流 分级分区河网的处理： 按子流域级别由高至低、支流由上游至下游的顺序排列各河网树 计算任务搜索在整个森林中进行 每株树中允许生成多个计算任务（子树） 将待算子树的河段编码按上下游顺序展平到数组中，作为计算任务包发送到空闲的计算节点 计算节点运行专业模型，依次计算任务包中的各河段 3.坡面产沙模型 单位面积上泥沙颗粒所占的面积比为： 单层泥沙颗粒个数为: 单层泥沙总质量为: 单位时间流出观察面的单层泥沙质量为: 定义e为单位时间单位面积上侵蚀产沙的质量，则： 4.河网汇流与输沙模型 为正确描述沟道中复杂的水沙现象，沟道系统的水沙运动模拟应满足以下要求： 考虑洪水涨落不同阶段沟道内泥沙的冲淤变化； 考虑沟道水流运动对重力侵蚀发生的影响作用，并将沟道作为滑坡和崩塌带来的干物质进入水流参与水沙运动的入口； 考虑沟道内泥沙级配的调整过程。 沟道不具备实测断面资料 只有通过由数字地形模型提取河网的过程获得的河段长度、河底平均坡降等基本参数 因此，水流运动模拟采用扩散波方法的形式，并略作改进以获得水深、流速等水力参数 水深形式的圣维南方程 贡日证明了马斯京根法可在一定条件下二阶逼近扩散波方程 为控制负反应，并同时满足收敛性要求，时间步长确定的情况下，取空间步长 V形断面假定C=4U/3 水力坡降取为 应用于V形断面顺直河道，与水动力学方法和马斯京根法对比，作为验证 河道参数 L=100km, S0=0.00025,m=4,n=0.035 为区分侵蚀与产沙的不同概念，反映大的河网范围内不同位置的冲淤分布，以及各河段的冲淤变化过程，沟道输沙模型采用不平衡输沙模式。 从水槽实验与部分渠道实测数据的 关系出发，得到 5.重力侵蚀模型 黄土区气候条件影响： 暴雨使浅沟和切沟迅速发展，为重力侵蚀提供地形条件 水流通过裂缝进入垂直节理或者渗入中下层土体，较高的含水量不仅破坏土壤中盐类的胶结作用，而且造成正的孔隙水压力，对沟坡稳定造成不利影响 坡脚沟道中的洪水在大暴雨时涨水很快，水位及流速都很高，这将直接导致对坡脚的掏刷 黄土区地形条件影响： 黄土湿陷性影响： 坡脚土体起动切应力 坡顶裂缝中的水压力 沟坡抗滑力的确定 非饱和土强度分析 附加粘聚力 的处理： 6.一维河道水动力学模型 为了能精确计算有断面资料部分河道的水沙运动，干流河段不再用精确扩散波模型，而采用一维水动力学法计算。 关键问题的处理 采用水沙方程的非耦合解法 离散采用Preissmann四点偏心隐格式 挟沙力计算采用张红武公式 不平衡输沙采用时空错位法求解 “二级悬河”这一特殊地形采用堰流加水库模式处理 含沙量与冲淤量在横断面方向根据水力特性分布 7.流域泥沙过程模拟的应用 河床变形方程 非平衡悬沙输移方程 运动方程式 水流连续方程 岔巴沟 多沙区 实验验证结果分析 模型参数比较稳定，且能够正确反映坡面产沙率的数值和趋势。 对于不同组实验，模型参数中的k值有一定变化，原因主要是不同实验的土壤特性差别较大，而参数k是与土壤特性有关的参数，因此也有一定变化。 三组数据都存在某些输沙率突然增大，远离其他点群趋势的情况，据原实验资料，出现这些较远偏离点时，在实验现象上一般都伴有细沟的形成、沟壁的坍塌滑塌等重力侵蚀现象。因此，要想模拟这些偏离点的产沙规律，还必须专门分析细沟侵蚀的产生以及重力侵蚀，得到相应的模型。 沟间侵蚀模型的验证 精确扩散波水流模型 马斯京根-贡日法 忽略动量方程的惯性项，两式联立得扩散波方程 马斯京根-贡日法 这一条件为 算法无条件稳定 主要要素的处理 水力要素的计算采用四点平均 主要要素的处理 有学者曾经指出，如果在扩散波方程系数的求解中，忽略 项的