黄河粗泥沙界限确定方法讨论.docVIP

冲积河流纵向冲淤和横向变形数学模型研究及应用 王新宏 张根广 曹如轩 沈 晋 （西安理工大学水利水电学院 西安 710048） 摘 要 本冲积河流中，河槽的纵向形态和横向形态都处于不断的冲淤调整之中。河流变形的数值模拟需要兼顾河槽的纵向冲淤和横向变形两个方面。许多河道变形的数学模型仅能够模拟河道的纵向冲淤变形，无法模拟因河岸（或滩地）坍塌造成的河槽横向变形过程。本文以冲积河流为对象，通过对河槽冲淤调整过程中河岸（或滩地）坍塌下滑的土动力学机理、岸壁泥沙坍塌的数量及其物理特性和冲淤特性等的分析研究，提出了一套可以同时反映河槽纵向冲淤和横向调整变形的河道变形模拟方法。以此为基础建立了一个可以模拟多沙河流上水库及河道冲淤演变过程的准二维泥沙数学模型。利用三门峡水库1960～1990年的实测资料，对本文建立的准二维泥沙数学模型进行了测试和验证，计算结果与实测值相符合。利用所建模型对小浪底水库不同运用方案的水库泥沙冲淤过程进行了预测计算，并研究了三门峡水库运用方式对潼关高程的影响。 关键词 河床演变 横向变形 河岸侵蚀 水库泥沙 数学模型 1 前言 冲积河流的冲淤调整是通过两方面进行的。一方面是通过床沙质来量和水流挟沙能力的对比关系使河床产生纵向冲淤变形，另一方面是通过河岸抗冲力和水流冲刷力之间的对比关系使河岸发生冲淤，导致河流产生横向变形。 实际上，河床纵向下切和河槽横向展宽是冲积河道冲刷过程中的普遍现象。河道冲刷过程中，导致河槽两岸（简称河岸，在本文中特指冲积河道的主槽两岸）失去稳定的外部原因是水流的冲刷作用，主要表现在两个方面：一方面，水流冲刷床面使河岸（或边滩）的相对高度增加，可导致河岸稳定性降低，这是河床的垂向水力侵蚀作用；另一方面，水流淘刷河岸坡脚使河岸的坡角变陡，同样可导致河岸的稳定性降低，这是河岸的横向(或侧向)水力侵蚀作用。导致河岸失去稳定的内部原因是岸壁土体内部的剪切应力超过了土体的抗剪强度。 在河道冲刷过程中，随着河床的冲刷下切，河岸高度不断增加，岸坡不断变陡。由于河岸坡脚处土体含水量通常较大，抗剪强度较小，加之土体内孔隙水不断析出，从而使土体产生垂向变形，导致滩地出现平行于主槽的裂缝。随着裂缝的出现和发展，作用于岸壁土体底部的重力将不断增加，其结果又促进变形进一步发展。这样，当剪切应力大于土体剪切强度时，河岸坍塌现象便发生了。河岸发生坍塌的地方，往往是水流顶冲河岸的地方。河岸泥沙坍塌进入主槽中后，常形成拦水坎，水流变得更为湍急。在水流的冲刷作用下，拦水沙坎逐渐崩溃，大量的泥核、土块等泥沙团粒因来不及溶解而随水流一起被带走。经过一段流程后，随着泥沙团粒的不断粉碎，最后完全溶解。 由此可见，河槽的纵向冲淤和横向变形是冲积河道冲淤调整过程中不可分割的两个方面，两者既相互联系，又相互区别。冲积河流河床变形的模拟必须同时考虑这两个方面。 然而，目前大多数泥沙数学模型仅能够模拟和预测河道的纵向冲淤变形，不模拟河岸的冲淤过程，无法预测河槽横向宽度随时间发生的调整变化[1]。天然河道中除坡降等表征纵向形态的物理量之外，河道宽度等表征横向形态的物理量也随着时间的推移处于不断的调整变化之中，河道的横向变形与纵向变形之间具有相互影响和相互制约的关系。在河流冲淤变形的数值模拟中，不考虑河岸的冲淤过程，是一个致命的缺陷。 为了能够模拟河道的横向变形过程，国内外已有一些学者对河岸冲淤变形的物理过程和力学机理进行了探讨，并通过理论分析或实测资料分析，建立了一些可用于预测河宽调整变化的泥沙数学模型。总体来看，国外的模型其研究对象大多数仅限于边界条件极为简单的模型小河或渠道，由于水流冲刷床面和坡脚造成的岸壁侵蚀过程，而且不考虑漫滩水流对河岸侵蚀过程的影响。国内的摸型其研究对象大多为天然河流，如黄河。除个别模型考虑了河道横向摆动的内在力学机理外，其它模型(尤其是1D模型)都是采用一些经验方法将冲淤面积分配到横断面上[1]。 本文以冲积河流为研究对象，在已有的一维泥沙数学模型研究的基础上，通过对河槽冲淤调整过程中河岸（或滩地）坍塌下滑的土动力学机理、岸壁泥沙坍塌的数量及其物理特性和冲淤特性等的分析研究，提出了一套可以同时反映河槽纵向冲淤和横向调整的河道变形模拟方法。以此为基础建立一个可以描述多沙河流上水库及河道冲淤演变过程的准二维泥沙数学模型。 2 数学模型的建立 本文建立的准二维泥沙数学模型由水动力学模块、泥沙输移模块和河岸变形模块三大模块组成。以下介绍所采用的基本理论和计算方法。 2.1 水流流场的模拟[2] 一维水流模型求解方法简单，但只能给出各个断面的平均水力要素，无法满足河道横向冲淤变形计算的需要。平面二维水流模型虽然能给出水力要素沿纵向和沿横向的详细变化，但由于其求解方法相对较为复杂，不适用于长历时长河段的河