渭河下游泥沙数学模型的改进
|
摘要
近年来，由于水沙条件发生变化，以及河槽边界条件的改变，使得渭河出现了许多新的防洪问题，其特点是：“流量小、水位高、灾害大”。频繁发生的异常洪水位现象给渭河下游防洪增加了困难。针对这种现状，有必要系统地研究渭河下游的洪水行进规律，认清渭河现有河床中发生的洪水特性，这对于渭河下游的河道治理、设计洪水位的合理确定及对渭河下游的洪水预报研究都具有重要的意义。因此，为促进渭河下游洪水的预测预报和风险管理，确保防洪安全，迫切需要利用数学模型来研究渭河下游洪水演进的各种现象和问题，为洪水的监测、防洪和抢险指挥提供参考。
目前，在黄河中下游上得到应用的一、二维水沙模型有多个，基本上都是由描述水沙运动的四个基本方程构成，即水流连续方程、水流运动方程、泥沙连续方程与河床变形方程，不同之处在于各模型在求解基本方程时，对某些参数的具体取值和计算方法有所不同，如反映不平衡输沙规律中的挟沙力问题、动床阻力、悬移质泥沙与床沙的交换机理等。然而，与黄河中下游相比，渭河下游发生高含沙洪水的频率高，洪水历时长，含沙量大，高含沙水流的静态剪切应力作用明显，造成河床断面形态变形剧烈，尤其在横向变形上显得更为突出。而现有的渭河下游数学模型，主要是沿用黄河中下游的一维数学模型，模型结构反映更多的是黄河中下游一般高含沙挟沙水流的输沙特性，并不能合理地反映出渭河下游的输沙特性。因此，有必要针对渭河下游泥沙问题，在挟沙规律和横向冲淤变形等方面去改进现有的渭河下游一维恒定流泥沙数学模型。本研究论文建立了常规一维数学模型，并对其作了部分改进，以期达到适应渭河特点的目的。工作主要表现在以下几方面。
1、在泥沙计算中，考虑到行洪过程中泥沙在冲刷和淤积情况下，悬移质泥沙与床沙交换机理的不同，首次在冲淤计算中采用物理含义更明确的双值挟沙力公式，并对双值挟沙力公式中的系数计算问题进行了详细的分析及开发出与该公式应用相配套的应用方法。
2、在河床演变部分，对渭河下游特有的“贴边淤积”现象进行了机理探讨，分析研究表明，造成渭河下游河槽急剧淤积萎缩变形的重要原因之一为高含沙洪水的静态剪切应力，它是浑水的流型从牛顿体转变成凝胶体后的产物。
3、为了使数学模型在模拟过程中能反映出高含沙洪水静态剪切应力造成的“贴边淤积”，作者根据实测资料建立经验公式对“贴边淤积”进行了概化处理。
4、最后利用1994年高含沙小洪水过程对模型进行验证，模拟结果表明模型具有一定的模拟能力。相比较而言，经本论文改进后的模型在反映物理成因概念方面较同类型模型清晰、准确，也更接近渭河下游实际情况。
In recent years, many flood control problems resulting from small river flux with high river stage and lengthen in the time of flood wave travel have happened because of the condition change of river water , sediment and variation of channel boundary in the downstream reach of Weihe River, which brings the case of small discharge but a heavy flood damage. The frequently occurring unordinary flood stage increases the difficulty of flood control practice in the river. For the river regulation, the rational determination of design flood table and building a study base of flood forecast for the lower reaches of Weihe River, it is quite important to study systematically the flood travel and to recognize the flood character appearing under the existing river bed conditions of Weihe River. In order to improve the forecasting risk management and make a safe flood control, there is a active demand to explore the many phenomena and probing problems existing in hyper-concentration sediment flood travel by mathematical model in the lower reaches of Weihe River, which could supply a reference for monitoring flood and directing actions on flood control.
At the present, many 1-D and 2-D mathematical models being used in the middle and lower reaches of Yellow River build with the flow equation, the continuous equation, the sediment continuous equation and the channel bed deformation equation. The differences among them come from the special parameters used and the method of model operation, such as sediment-carrying capacity of flow, movable-bed resistance, and the way of exchange between the suspended load and the bed load and so on, which reflect non-equilibrium transportation of sediment in stream channel. In contrast with Yellow River, the lower reaches of Wei River happen more frequently the hyper-concentration flood. The flood has the character of long wave travel time, hyper-concentration sediments and more manifest static shearing stress which results in abrupt change of the river bed cross section, especially in the lateral deformation. Now, the mathematical model for the lower reaches of Wei River is mainly the one of 1-D mathematical model used in middle and lower reaches of Yellow River. So the model more reflects the character of river sediment delivery of general muddy flow in Yellow River but that of hyper-concentration flow in lower Wei River. In order to explain the characteristic of Wei River, it is necessary to improve the model according to the character of sediment delivery and the lateral deformation in the Wei River. This thesis builds up a common 1-D mechanical model with modifications as follows so as to hopefully reach the goal.
1.In the aspect of sediment computation, considering the existing of different exchange mechanics between suspended load and bed load under the condition of water scouring and silting process in flood travel, this thesis adopt firstly the marginal scouring and marginal silting sediment concentration capacity formula which possess more definite physical meaning and developed a methods for its application.
2.In the part of evolution of river bed, the specific phenomena of integral depositing at the channel boundary often occurred at downstream river reach of Weihe River has been probed into from a physical base, and the results of this analysis and study indicated that one of the most important causes that cause the river channel of the downstream channel of Weihe River deform rapidly is the static yield stress of hyper-concentration flood. And this static share stress is a product of muddy water stream pattern changing from Newtonian into gel.
3. For fitting the characteristic of hyper-concentration flood with mathematic model, in the process of simulation, the integral depositing was been made a generalized disposal through building empirical formula according to observed data.
4. Lastly, the model verification has been carried out using the hydrograph of observed flood of year 1994, and the results of the simulation indicated that the model has preferably simulation capability. Comparison with the same kind models, the modified model in this thesis is more distinct, exact and practicable at the aspect of reflecting the physical base conception.
目前，在黄河中下游上得到应用的一、二维水沙模型有多个，基本上都是由描述水沙运动的四个基本方程构成，即水流连续方程、水流运动方程、泥沙连续方程与河床变形方程，不同之处在于各模型在求解基本方程时，对某些参数的具体取值和计算方法有所不同，如反映不平衡输沙规律中的挟沙力问题、动床阻力、悬移质泥沙与床沙的交换机理等。然而，与黄河中下游相比，渭河下游发生高含沙洪水的频率高，洪水历时长，含沙量大，高含沙水流的静态剪切应力作用明显，造成河床断面形态变形剧烈，尤其在横向变形上显得更为突出。而现有的渭河下游数学模型，主要是沿用黄河中下游的一维数学模型，模型结构反映更多的是黄河中下游一般高含沙挟沙水流的输沙特性，并不能合理地反映出渭河下游的输沙特性。因此，有必要针对渭河下游泥沙问题，在挟沙规律和横向冲淤变形等方面去改进现有的渭河下游一维恒定流泥沙数学模型。本研究论文建立了常规一维数学模型，并对其作了部分改进，以期达到适应渭河特点的目的。工作主要表现在以下几方面。
1、在泥沙计算中，考虑到行洪过程中泥沙在冲刷和淤积情况下，悬移质泥沙与床沙交换机理的不同，首次在冲淤计算中采用物理含义更明确的双值挟沙力公式，并对双值挟沙力公式中的系数计算问题进行了详细的分析及开发出与该公式应用相配套的应用方法。
2、在河床演变部分，对渭河下游特有的“贴边淤积”现象进行了机理探讨，分析研究表明，造成渭河下游河槽急剧淤积萎缩变形的重要原因之一为高含沙洪水的静态剪切应力，它是浑水的流型从牛顿体转变成凝胶体后的产物。
3、为了使数学模型在模拟过程中能反映出高含沙洪水静态剪切应力造成的“贴边淤积”，作者根据实测资料建立经验公式对“贴边淤积”进行了概化处理。
4、最后利用1994年高含沙小洪水过程对模型进行验证，模拟结果表明模型具有一定的模拟能力。相比较而言，经本论文改进后的模型在反映物理成因概念方面较同类型模型清晰、准确，也更接近渭河下游实际情况。
In recent years, many flood control problems resulting from small river flux with high river stage and lengthen in the time of flood wave travel have happened because of the condition change of river water , sediment and variation of channel boundary in the downstream reach of Weihe River, which brings the case of small discharge but a heavy flood damage. The frequently occurring unordinary flood stage increases the difficulty of flood control practice in the river. For the river regulation, the rational determination of design flood table and building a study base of flood forecast for the lower reaches of Weihe River, it is quite important to study systematically the flood travel and to recognize the flood character appearing under the existing river bed conditions of Weihe River. In order to improve the forecasting risk management and make a safe flood control, there is a active demand to explore the many phenomena and probing problems existing in hyper-concentration sediment flood travel by mathematical model in the lower reaches of Weihe River, which could supply a reference for monitoring flood and directing actions on flood control.
At the present, many 1-D and 2-D mathematical models being used in the middle and lower reaches of Yellow River build with the flow equation, the continuous equation, the sediment continuous equation and the channel bed deformation equation. The differences among them come from the special parameters used and the method of model operation, such as sediment-carrying capacity of flow, movable-bed resistance, and the way of exchange between the suspended load and the bed load and so on, which reflect non-equilibrium transportation of sediment in stream channel. In contrast with Yellow River, the lower reaches of Wei River happen more frequently the hyper-concentration flood. The flood has the character of long wave travel time, hyper-concentration sediments and more manifest static shearing stress which results in abrupt change of the river bed cross section, especially in the lateral deformation. Now, the mathematical model for the lower reaches of Wei River is mainly the one of 1-D mathematical model used in middle and lower reaches of Yellow River. So the model more reflects the character of river sediment delivery of general muddy flow in Yellow River but that of hyper-concentration flow in lower Wei River. In order to explain the characteristic of Wei River, it is necessary to improve the model according to the character of sediment delivery and the lateral deformation in the Wei River. This thesis builds up a common 1-D mechanical model with modifications as follows so as to hopefully reach the goal.
1.In the aspect of sediment computation, considering the existing of different exchange mechanics between suspended load and bed load under the condition of water scouring and silting process in flood travel, this thesis adopt firstly the marginal scouring and marginal silting sediment concentration capacity formula which possess more definite physical meaning and developed a methods for its application.
2.In the part of evolution of river bed, the specific phenomena of integral depositing at the channel boundary often occurred at downstream river reach of Weihe River has been probed into from a physical base, and the results of this analysis and study indicated that one of the most important causes that cause the river channel of the downstream channel of Weihe River deform rapidly is the static yield stress of hyper-concentration flood. And this static share stress is a product of muddy water stream pattern changing from Newtonian into gel.
3. For fitting the characteristic of hyper-concentration flood with mathematic model, in the process of simulation, the integral depositing was been made a generalized disposal through building empirical formula according to observed data.
4. Lastly, the model verification has been carried out using the hydrograph of observed flood of year 1994, and the results of the simulation indicated that the model has preferably simulation capability. Comparison with the same kind models, the modified model in this thesis is more distinct, exact and practicable at the aspect of reflecting the physical base conception.
引文
[1] 李勇，董雪娜，张晓华等．黄河水沙特性变化研究[M] ．郑州：黄河水利出版社，2004．205～245
[2] 张仁．潼关高程升高及其解决方法．三门峡水利枢纽运用四十周年论文集[M] ．2000
[3] 周文浩，陈建国，李慧梅．潼关高程及遏制渭河下游淤积的对策．三门峡水利枢纽运用四十周年论文集[M] ．2000
[4] 胡建华，宋红霞，曹廷立．渭河下游当前防洪形势．西北水资源与水工程[J] ，2000年第4期
[5] 陈建国，胡春宏，戴清．渭河下游近期河道萎缩特点及治理对策．泥沙研究[J] ，2002年第6期
[6] 冯普林，石长伟，张广林．渭河“2003”洪水灾害及其减灾措施的分析．中国水利水电科学研究院学报[J] ，2004年第1期
[7] 曹如轩．三门峡水库淤积上延机理的研究．泥沙研究[J] ，2001年第2期
[8] 秦毅，郑学萍等．渭河下游洪水变化特点及成因．干旱区资源与环境[J] ，2006年第3期
[9] 谢鉴衡，魏良琰．河流泥沙数学模型的回顾与展望．泥沙研究[J] ，1987年第3期
[10] Feldman, A. D. HEC Models For Water Resources System Simulation: Thee ry and Experience. The Hydraulic Engineering Ceter, Davis, California, 1981
[11] Chang, H. H. Fluvial processes in river engineering. Wiley Inter—science. New York, N. Y, 429. 1988
[12] Chang, H. H. , Harrison, L. L. , Lee, W. , Tu, S. , Numerical Modeling For Sediment—Pass—Through Reservoirs, J. Hydr. Engrg. , ASCE, 122(7), 1996
[13] Yang, C. T. , Molinas, A. , and Song, C. S"GSTARS—Generalized Stream Model for alluvial river Simulation. "Twelve selected computer stream sedimentation Models developed in the U. S. , S. Fan ed. , Energy Regulatory Commission, Washington, D. C. 1988
[14] Borah, D. K. , Bordoloi, P. K. "Stream bank erosion and bed evolution model. "Sediment Transport Modeling, S. Wang ed. , ASCE, New York. 1989
[15] Borah, D. K. , Dashputre, M. S. Field evaluation of the sediment transport model STREAM with a bank erosion component. Proc. Hydr. Engrg. '94, G. V Cotroneo and R. R. Rumer, eds. , ASCE, New York, N. Y, 1994
[16] Osman, A. M. , Thorne, C. R, Riverbank stability analysis. Ⅰ: Theory. J Hydr. Engrg. , ASCE, 114(2), 1988
[17] Thome, C. R. , Osman, A. M. , Riverbank stability analysis. Ⅱ: Applications. J. of Hydr. Engrg. , 114(2), 1988
[18] Fan, S. , ed. Twelve selected computer stream sedimentation models developed in the United States. Fed. Energy Regulatory Commission, Washington, D. C. 1988
[19] 李义天，尚全民．一维不恒定流泥沙数学模型研究．泥沙研究[J] ，1998年第1期
[20] 杨国录．河流数学模型[M] ．海洋出版社，1993
[21] 王兆印，宋振琪．欧美泥沙运动研究述评．第二届全国泥沙基本理论研究学术讨论会论文集[C] ，北京，1995
[22] ASCE Task Committee on Hydraulics, Bank Mechanics and Modeling of River Width Adjustment. River width adjustment. Ⅱ: Modeling. J. Hydr. Engrg. , ASCE, 124(9). 1998
[23] 许协庆，朱鹏程．河床变形问题的特征线解．水利学报[J] ，1963第5期
[24] 窦国仁．潮汐水流中悬沙运动及冲淤计算．水利学报[J] ，1963第6期
[25] 韩其为．悬移质不平衡输沙的研究．河流泥沙国际学术讨论会论文集[C] ，北京，1980
[26] 林秉南．关于一维动床数学模型的讨论．三峡工程泥沙问题研究成果汇编(160—180m蓄水位方案)[C] ，水利电力部科学技术司，1988
[27] 夏军强，王光谦，吴保生．游荡型河流演变及其数值模拟[M] ．北京：中国水利水电出版社，2005
[28] 王光谦，张红武，夏军强．游荡型河流演变及模拟[M] ．北京：中国科学出版社，2005
[29] 梁志勇，伊学良．冲积河流河床横向变形的初步数学模拟．泥沙研究[J] ，1991年第4期
[30] 周建军，林秉南，王连祥．平面二维泥沙数学模型的研究及应用．水利学报[J] ，1993年第11期
[31] 梁国亭，张仁．黄河小北干流一维泥沙冲淤数学模型．人民黄河[J] ，1996年第9期
[32] 王新宏，程文，曹如轩，钱善琪．黄河中游一维悬移质泥沙数学模型．西安理工大学学报[J] ，1996年第3期
[33] 程文，王新宏等．黄河龙潼段非恒定流泥沙数学模型．西安理工大学学报[J] ，1998年第2期
[34] 王新宏．冲积河道纵向冲淤和横向变形数值模拟研究及应用：博士学位论文[D] ．西安：西安理工大学，2000
[35] 黄金池，万兆惠．黄河下游河床平面变形模拟研究．水利学报[J] ，1999年第2期
[36] 夏军强．河岸冲刷机理研究及数值模拟：博士学位论文[D] ．北京：清华大学，2002
[37] 于清来，窦国仁．高含沙河流泥沙数学模型研究．水利水运科学研究[J] ，1999年第2期
[38] 张红艺，杨明，张俊华，邓洪亮．高含沙水库泥沙运动数学模型的研究及应用．水利学报[J] ，2001年第11期
[39] 赵振武．渭河下游洪水演进数值模拟及防洪减灾对策研究：硕士学位论文[D] ．武汉：武汉大学，2003
[40] 张强．渭河下游非恒定流泥沙数学模型研究：硕士学位论文[D] ．西安：西安理工大学，2003
[41] 杨方社．潼关高程对渭河下游冲淤影响的数值模拟研究：硕士学位论文[D] ．西安：西安理工大学，2004
[42] 黄修山．渭河下游河道淤积萎缩对洪水演进规律的影响研究：硕士学位论文[D] ．西安：西安理工大学，2005
[43] 钱宁，张仁，周志德．河床演变学[M] ．北京：科学出版社，1987
[44] 钱意颖等，高含沙均质水流基本特性的实验研究．见：齐璞等主编，黄河高含沙水流运动规律及应用前景[M] ．北京：科学出版社，1993
[45] 钱宁．高含沙水流运动[M] ．北京：清华大学出版社，1989
[46] 钱宁，万兆惠．泥沙运动力学[M] ．北京：科学出版社，1983
[47] 杨美卿，钱宁．紊动对细泥沙浆液絮凝结构的影响．水利学报[J] ，1986年第8期
[48] 费祥俊，朱程清．高含沙水流运动中的宾汉切应力．泥沙研究[J] ，1991年第4期
[49] 陈立，王明甫．水流强度对高含沙水流流变参数影响的试验研究．泥沙研究[J] ，1993年第2期
[50] 曹如轩等．高含沙异重流阻力规律．第二次河流泥沙国际学术讨论会论文集[M] ．北京：水利电力出版社，1983
[51] 郑学萍．静态高含沙浑水剪切力及对渭河洪水影响的初步研究：硕士学位论文[D] ．西安：西安理工大学，2006
[52] 胡春宏，吉祖稳．复式断面边界剪切应力分布规律研究．泥沙研究[J] ，1999年第6期
[53] 张瑞瑾．河流泥沙运力学[M] ．北京：中国水利水电出版社，1989
[54] 曹如轩．高含沙水流挟沙力初步研究．水利水电技术[J] ，1979年第5期
[55] 舒安平．水流挟沙力公式的验证与评述．水利学报[J] ，1993年第1期
[56] 吴保生等．黄河水流挟沙能力公式的若干修正．人民黄河[J] ，1993年第7期
[57] 张红武．黄河水流挟沙力的计算公式．人民黄河[J] ，1992年第1期
[58] 陈雪峰，陈立，李义天．高、中、低浓度挟沙水流挟沙力公式的对比分析．武汉水利电力大学学报[J] ，1999年第5期
[59] 沙玉清．泥沙运动学引论[M] ．西安：陕西科学技术出版社，1996
[60] 邓贤艺，曹如轩，钱善琪．水流挟沙力双值关系研究．水利水电技术[J] ，2000年第9期
[61] 曹如轩，程文，钱善琪．高含沙洪水揭河底冲刷初探．人民黄河[J] ，1997年第2期
[62] 费祥俊．黄河中下游含沙水流粘度的计算模型．泥沙研究[J] ，1991年第2期
[63] 吴保生，马吉星．床沙代表粒径与输移泥沙中值粒径的关系．泥沙研究[J] ，2002年第2期
[64] 钟德钰，王士强，王光谦．河流冲泻质挟沙力研究．泥沙研究[J] ，1998年第3期
[65] 钟德钰，王士强，王光谦．水流冲泻质挟沙机理探讨．泥沙研究[J] ，1998年第3期
[66] 钱宁，万兆惠．高含沙水流运动研究述评．水利学报[J] ，1985年第5期
[67] 万兆惠．高含沙水流运动研究续评．全国泥沙基本理论研究学术讨论会论文集[C] ．1992年第2卷
[68] 陈界仁．高含沙水流运动立面二维数学模型．河海大学学报[J] ，1994年第4期
[69] 于清来，窦国仁．高含沙水流泥沙数学模型研究．水利水运科学研究[J] ，1999年第5期
[70] 黄金池，万兆惠．高含沙水流二维数值模拟．水动力学研究与进展[J] ，2001年第3期
[71] 尹学良．清水冲刷河床粗化研究．水利学报[J] ，1963年第1期
[72] 王士强．沙波运动与推移质测验．泥沙研究[J] ，1988年第4期
[73] 王士强．沙波变化与床沙交换调整．泥沙研究[J] ，1992年第4期
[74] 刘金梅，王士强，王光谦．冲积河流长距离冲刷不平衡输沙过程初步研究．沙波变化与床沙交换调整．水利学报[J] ，2002年第2期
[75] 刘金梅，王光谦，王士强．沙质河道冲刷不平衡输沙机理及规律研究．沙波变化与床沙交换调整．水利学报[J] ，2003年第5期
[2] 张仁．潼关高程升高及其解决方法．三门峡水利枢纽运用四十周年论文集[M] ．2000
[3] 周文浩，陈建国，李慧梅．潼关高程及遏制渭河下游淤积的对策．三门峡水利枢纽运用四十周年论文集[M] ．2000
[4] 胡建华，宋红霞，曹廷立．渭河下游当前防洪形势．西北水资源与水工程[J] ，2000年第4期
[5] 陈建国，胡春宏，戴清．渭河下游近期河道萎缩特点及治理对策．泥沙研究[J] ，2002年第6期
[6] 冯普林，石长伟，张广林．渭河“2003”洪水灾害及其减灾措施的分析．中国水利水电科学研究院学报[J] ，2004年第1期
[7] 曹如轩．三门峡水库淤积上延机理的研究．泥沙研究[J] ，2001年第2期
[8] 秦毅，郑学萍等．渭河下游洪水变化特点及成因．干旱区资源与环境[J] ，2006年第3期
[9] 谢鉴衡，魏良琰．河流泥沙数学模型的回顾与展望．泥沙研究[J] ，1987年第3期
[10] Feldman, A. D. HEC Models For Water Resources System Simulation: Thee ry and Experience. The Hydraulic Engineering Ceter, Davis, California, 1981
[11] Chang, H. H. Fluvial processes in river engineering. Wiley Inter—science. New York, N. Y, 429. 1988
[12] Chang, H. H. , Harrison, L. L. , Lee, W. , Tu, S. , Numerical Modeling For Sediment—Pass—Through Reservoirs, J. Hydr. Engrg. , ASCE, 122(7), 1996
[13] Yang, C. T. , Molinas, A. , and Song, C. S"GSTARS—Generalized Stream Model for alluvial river Simulation. "Twelve selected computer stream sedimentation Models developed in the U. S. , S. Fan ed. , Energy Regulatory Commission, Washington, D. C. 1988
[14] Borah, D. K. , Bordoloi, P. K. "Stream bank erosion and bed evolution model. "Sediment Transport Modeling, S. Wang ed. , ASCE, New York. 1989
[15] Borah, D. K. , Dashputre, M. S. Field evaluation of the sediment transport model STREAM with a bank erosion component. Proc. Hydr. Engrg. '94, G. V Cotroneo and R. R. Rumer, eds. , ASCE, New York, N. Y, 1994
[16] Osman, A. M. , Thorne, C. R, Riverbank stability analysis. Ⅰ: Theory. J Hydr. Engrg. , ASCE, 114(2), 1988
[17] Thome, C. R. , Osman, A. M. , Riverbank stability analysis. Ⅱ: Applications. J. of Hydr. Engrg. , 114(2), 1988
[18] Fan, S. , ed. Twelve selected computer stream sedimentation models developed in the United States. Fed. Energy Regulatory Commission, Washington, D. C. 1988
[19] 李义天，尚全民．一维不恒定流泥沙数学模型研究．泥沙研究[J] ，1998年第1期
[20] 杨国录．河流数学模型[M] ．海洋出版社，1993
[21] 王兆印，宋振琪．欧美泥沙运动研究述评．第二届全国泥沙基本理论研究学术讨论会论文集[C] ，北京，1995
[22] ASCE Task Committee on Hydraulics, Bank Mechanics and Modeling of River Width Adjustment. River width adjustment. Ⅱ: Modeling. J. Hydr. Engrg. , ASCE, 124(9). 1998
[23] 许协庆，朱鹏程．河床变形问题的特征线解．水利学报[J] ，1963第5期
[24] 窦国仁．潮汐水流中悬沙运动及冲淤计算．水利学报[J] ，1963第6期
[25] 韩其为．悬移质不平衡输沙的研究．河流泥沙国际学术讨论会论文集[C] ，北京，1980
[26] 林秉南．关于一维动床数学模型的讨论．三峡工程泥沙问题研究成果汇编(160—180m蓄水位方案)[C] ，水利电力部科学技术司，1988
[27] 夏军强，王光谦，吴保生．游荡型河流演变及其数值模拟[M] ．北京：中国水利水电出版社，2005
[28] 王光谦，张红武，夏军强．游荡型河流演变及模拟[M] ．北京：中国科学出版社，2005
[29] 梁志勇，伊学良．冲积河流河床横向变形的初步数学模拟．泥沙研究[J] ，1991年第4期
[30] 周建军，林秉南，王连祥．平面二维泥沙数学模型的研究及应用．水利学报[J] ，1993年第11期
[31] 梁国亭，张仁．黄河小北干流一维泥沙冲淤数学模型．人民黄河[J] ，1996年第9期
[32] 王新宏，程文，曹如轩，钱善琪．黄河中游一维悬移质泥沙数学模型．西安理工大学学报[J] ，1996年第3期
[33] 程文，王新宏等．黄河龙潼段非恒定流泥沙数学模型．西安理工大学学报[J] ，1998年第2期
[34] 王新宏．冲积河道纵向冲淤和横向变形数值模拟研究及应用：博士学位论文[D] ．西安：西安理工大学，2000
[35] 黄金池，万兆惠．黄河下游河床平面变形模拟研究．水利学报[J] ，1999年第2期
[36] 夏军强．河岸冲刷机理研究及数值模拟：博士学位论文[D] ．北京：清华大学，2002
[37] 于清来，窦国仁．高含沙河流泥沙数学模型研究．水利水运科学研究[J] ，1999年第2期
[38] 张红艺，杨明，张俊华，邓洪亮．高含沙水库泥沙运动数学模型的研究及应用．水利学报[J] ，2001年第11期
[39] 赵振武．渭河下游洪水演进数值模拟及防洪减灾对策研究：硕士学位论文[D] ．武汉：武汉大学，2003
[40] 张强．渭河下游非恒定流泥沙数学模型研究：硕士学位论文[D] ．西安：西安理工大学，2003
[41] 杨方社．潼关高程对渭河下游冲淤影响的数值模拟研究：硕士学位论文[D] ．西安：西安理工大学，2004
[42] 黄修山．渭河下游河道淤积萎缩对洪水演进规律的影响研究：硕士学位论文[D] ．西安：西安理工大学，2005
[43] 钱宁，张仁，周志德．河床演变学[M] ．北京：科学出版社，1987
[44] 钱意颖等，高含沙均质水流基本特性的实验研究．见：齐璞等主编，黄河高含沙水流运动规律及应用前景[M] ．北京：科学出版社，1993
[45] 钱宁．高含沙水流运动[M] ．北京：清华大学出版社，1989
[46] 钱宁，万兆惠．泥沙运动力学[M] ．北京：科学出版社，1983
[47] 杨美卿，钱宁．紊动对细泥沙浆液絮凝结构的影响．水利学报[J] ，1986年第8期
[48] 费祥俊，朱程清．高含沙水流运动中的宾汉切应力．泥沙研究[J] ，1991年第4期
[49] 陈立，王明甫．水流强度对高含沙水流流变参数影响的试验研究．泥沙研究[J] ，1993年第2期
[50] 曹如轩等．高含沙异重流阻力规律．第二次河流泥沙国际学术讨论会论文集[M] ．北京：水利电力出版社，1983
[51] 郑学萍．静态高含沙浑水剪切力及对渭河洪水影响的初步研究：硕士学位论文[D] ．西安：西安理工大学，2006
[52] 胡春宏，吉祖稳．复式断面边界剪切应力分布规律研究．泥沙研究[J] ，1999年第6期
[53] 张瑞瑾．河流泥沙运力学[M] ．北京：中国水利水电出版社，1989
[54] 曹如轩．高含沙水流挟沙力初步研究．水利水电技术[J] ，1979年第5期
[55] 舒安平．水流挟沙力公式的验证与评述．水利学报[J] ，1993年第1期
[56] 吴保生等．黄河水流挟沙能力公式的若干修正．人民黄河[J] ，1993年第7期
[57] 张红武．黄河水流挟沙力的计算公式．人民黄河[J] ，1992年第1期
[58] 陈雪峰，陈立，李义天．高、中、低浓度挟沙水流挟沙力公式的对比分析．武汉水利电力大学学报[J] ，1999年第5期
[59] 沙玉清．泥沙运动学引论[M] ．西安：陕西科学技术出版社，1996
[60] 邓贤艺，曹如轩，钱善琪．水流挟沙力双值关系研究．水利水电技术[J] ，2000年第9期
[61] 曹如轩，程文，钱善琪．高含沙洪水揭河底冲刷初探．人民黄河[J] ，1997年第2期
[62] 费祥俊．黄河中下游含沙水流粘度的计算模型．泥沙研究[J] ，1991年第2期
[63] 吴保生，马吉星．床沙代表粒径与输移泥沙中值粒径的关系．泥沙研究[J] ，2002年第2期
[64] 钟德钰，王士强，王光谦．河流冲泻质挟沙力研究．泥沙研究[J] ，1998年第3期
[65] 钟德钰，王士强，王光谦．水流冲泻质挟沙机理探讨．泥沙研究[J] ，1998年第3期
[66] 钱宁，万兆惠．高含沙水流运动研究述评．水利学报[J] ，1985年第5期
[67] 万兆惠．高含沙水流运动研究续评．全国泥沙基本理论研究学术讨论会论文集[C] ．1992年第2卷
[68] 陈界仁．高含沙水流运动立面二维数学模型．河海大学学报[J] ，1994年第4期
[69] 于清来，窦国仁．高含沙水流泥沙数学模型研究．水利水运科学研究[J] ，1999年第5期
[70] 黄金池，万兆惠．高含沙水流二维数值模拟．水动力学研究与进展[J] ，2001年第3期
[71] 尹学良．清水冲刷河床粗化研究．水利学报[J] ，1963年第1期
[72] 王士强．沙波运动与推移质测验．泥沙研究[J] ，1988年第4期
[73] 王士强．沙波变化与床沙交换调整．泥沙研究[J] ，1992年第4期
[74] 刘金梅，王士强，王光谦．冲积河流长距离冲刷不平衡输沙过程初步研究．沙波变化与床沙交换调整．水利学报[J] ，2002年第2期
[75] 刘金梅，王光谦，王士强．沙质河道冲刷不平衡输沙机理及规律研究．沙波变化与床沙交换调整．水利学报[J] ，2003年第5期