大型湖泊水库平面及垂向二维流场与水质数值模拟
|
摘要
由于实际工程中湖泊、水库、河流、河口区域边界几何形状的不规则，以及水动力学方程的非线性性质，使得理论分析解难以作出，只能求助于数值解。数值模拟具有成本低、速度快、资料完备、具有改变各种参数、模拟真实条件及理想条件的能力。计算水力学中所求解的是非线性偏微分(积分)方程组，其数值方法的现有数学理论尚不够充分，严格的稳定性分析、误差估计和收敛性证明等理论工作的发展还跟不上数值模拟方法的进展。所以，需依靠对非线性数学方程的数值实验及数值解与试验或典型算例结果的比较和物理特性分析，验证数值计算结果，进一步改进计算方法。大多数生产实践问题要求进行二维水流水质数值模拟。在海岸、河口、湖泊、大型水库等广阔水域地区，水平尺度远大于垂向尺度，水力参数(如流速等)在垂直方向变化要小于水平方向的变化，其流场可用沿水深的平均流动量来表示，可采用平面二维水流流场与水质数值模拟技术。而在另外一些水域，如窄深潮汐通道、窄深型河谷水库如三峡、窄深河口地区，有关参量(如流速、温度、浓度等)的垂向变化要比水平横向的变化为大，应采用垂向二维水流水质数值模拟技术。
西安理工大学博士学位论文
本文主要工作:对二维非恒定流，既处理了不规则域，又研究了各种算法;不仅
做了平面二维，而且做了垂向二维;不仅计算了水流流场，也计算了水质变化;包括
了算例算法的研究和应用，也结合了Gls技术;平面二维模型采用了边界拟合坐标变
换，适应不规则的区域形状;垂向二维模型结合口坐标变换，可精确跟踪自由表面与
非规则的地形边界。对水平尺度远大于垂向尺度，水力参数(如流速等)在垂直方向
变化要小于水平方向变化的广阔水域地区，其流场可用沿水深的平均流动量来表示。
考虑自由水面与水深的影响，建立坐标变换、交替方向、交错网格、线性递推、隐格
式的平面二维水流与水质模型。对宽深比较小，有关参量(如流速、温度、浓度)
的垂向变化比水平横向的变化为大的窄深型河谷水库等区域，建立垂向二维水流与
水温数学模型。垂向做叮变换将可变域变为规则区域，使计算精确拟合自由面与底
部边界。采用算子分裂法、交错网格法、有限差分法进行计算。
所建立的平面二维模型和垂向二维模型均在实际工程中得到应用，其中:
在大庆市洪水模拟中采用平面二维非恒定流模型，ro天洪水可在1小时内计
算完毕，实现洪水演进的快速预报，为指挥防洪提供科学依据;采用贴体坐标技术，
可完全适应不规则的区域、河流边界;同时将GIS、遥感技术和洪水模拟技术有机
地结合起来，实现实时动态地显示、率定模型参数，显示洪水演进过程、实时地查
询洪水演进过程的参数等。
在博斯腾湖大水体研究中采用平面二维水流与水质模型，可完全适应不规则的
湖岸边界，考虑表面风应力，对吞吐量、风力、泵站位置、水位四种因素的多种组
合情况下的博斯腾湖风生环流、矿化度分布及出流水质进行了模拟计算，获得了一
些有意义的研究成果。
糯扎渡水库属峡谷型高坝大库，对糯扎渡水电站，采用垂向二维水流数学模型，
模拟了水电站机组启动过程的非恒定水流流动，模拟结果是合理的。这对于中国西
南地区类似的峡谷型高坝大库的水电站与三峡水电站垂向二维非恒定流模拟具有
参考价值。
In practical projects, the boundary shape irregularity of the lake, reservoir, river and estuary area and the non-linear character of hydrodynamic equation make it difficult to get the theoretical solution of the equation, so the numerical solution has to be the only a choice. The numerical simulation is of some advantages such as low cost, quick speed, complete data, and being capable to change parameters and simulate both the real and the ideal conditions. The non-linear partial differential (integral) equation often needs to be solved in the computational hydraulics, however, the mathematical theory of its numerical method has yet not been sufficient at present, for example, the development of some theories such as the strict stability analysis, error estimate and convergence proof has been unable to keep pace with the development of the numerical simulation method. Therefore, it is necessary to test the computation results and further modify the computation method based on the comparison between the nume
rical solution of the linear or non-linear equations and the trial value or the computation results of typical cases. The two-dimension numerical simulation of water flow and water quality is required in most practical projects, in which two types are
divided according to the different water area. In the broad water area such as the coast, estuary, lake and large reservoir, the horizontal scale is far larger than the vertical scale, the variation of the hydraulic parameters (for instance, velocity) is smaller in the vertical direction than in the horizontal direction, and the flow field can be expressed as the average flowing value along the water depth, so the two-dimension plane numerical simulation technique of the water flow and water quality may be used in these areas. But, in some other water areas such as the narrow-deep tide channel, narrow-deep valley reservoir and narrow-deep estuary area, the vertical varies of parameters such as the velocity, temperature and concentration is far greater than that in the horizontal direction, so the two-dimension vertical numerical simulation technique should be adopted in these areas.
The main contents in this paper are as follows: simulated irregular boundary and researched varies of algorithms; studied the two-dimension plane and two-dimension vertical computing; simulated water flow filed and water quality; researched algorithms and it's examples and combined with GIS tech. In two-dimension plane model the boundary fitted coordinate transformation were adopted to fit in with irregular area shape. In two-dimension vertical model a coordinate transformation were used to following the tracks of free surface and irregular topography boundary. In this paper, in order to fit in with the irregular boundary shape, and further to realize the quick forecast of the water flow and water quality, the two-dimension plane water flow and water quality model were built, in which the techniques of the coordinate transformation, alternate direction, staggered grid, linear recursion and implicit schemes were included. As for the reservoir of which the width to depth ratio is smaller, and the vertical variation of the parameters related such as the velocity, temperature and concentration is greater than their variation in the horizontal direction, the vertical two- dimension water flow model were set up, in which the water depth may change with the time and space, and when the coordinate transformation was made along the vertical direction, the variable area could be changed
into the regular one. In this model, the operator splitting technique, staggered grid technique and finite difference method were used.
The two-dimension plane model and two-dimension vertical model developed have been used in practical projects, in which:
In the flood simulation of Daqing City, the two-dimension plane unsteady flow model has been built to fit in with the irregular river boundary and realize the rapid forecast of the water flow. In this model, the geographic information sys
西安理工大学博士学位论文
本文主要工作:对二维非恒定流，既处理了不规则域，又研究了各种算法;不仅
做了平面二维，而且做了垂向二维;不仅计算了水流流场，也计算了水质变化;包括
了算例算法的研究和应用，也结合了Gls技术;平面二维模型采用了边界拟合坐标变
换，适应不规则的区域形状;垂向二维模型结合口坐标变换，可精确跟踪自由表面与
非规则的地形边界。对水平尺度远大于垂向尺度，水力参数(如流速等)在垂直方向
变化要小于水平方向变化的广阔水域地区，其流场可用沿水深的平均流动量来表示。
考虑自由水面与水深的影响，建立坐标变换、交替方向、交错网格、线性递推、隐格
式的平面二维水流与水质模型。对宽深比较小，有关参量(如流速、温度、浓度)
的垂向变化比水平横向的变化为大的窄深型河谷水库等区域，建立垂向二维水流与
水温数学模型。垂向做叮变换将可变域变为规则区域，使计算精确拟合自由面与底
部边界。采用算子分裂法、交错网格法、有限差分法进行计算。
所建立的平面二维模型和垂向二维模型均在实际工程中得到应用，其中:
在大庆市洪水模拟中采用平面二维非恒定流模型，ro天洪水可在1小时内计
算完毕，实现洪水演进的快速预报，为指挥防洪提供科学依据;采用贴体坐标技术，
可完全适应不规则的区域、河流边界;同时将GIS、遥感技术和洪水模拟技术有机
地结合起来，实现实时动态地显示、率定模型参数，显示洪水演进过程、实时地查
询洪水演进过程的参数等。
在博斯腾湖大水体研究中采用平面二维水流与水质模型，可完全适应不规则的
湖岸边界，考虑表面风应力，对吞吐量、风力、泵站位置、水位四种因素的多种组
合情况下的博斯腾湖风生环流、矿化度分布及出流水质进行了模拟计算，获得了一
些有意义的研究成果。
糯扎渡水库属峡谷型高坝大库，对糯扎渡水电站，采用垂向二维水流数学模型，
模拟了水电站机组启动过程的非恒定水流流动，模拟结果是合理的。这对于中国西
南地区类似的峡谷型高坝大库的水电站与三峡水电站垂向二维非恒定流模拟具有
参考价值。
In practical projects, the boundary shape irregularity of the lake, reservoir, river and estuary area and the non-linear character of hydrodynamic equation make it difficult to get the theoretical solution of the equation, so the numerical solution has to be the only a choice. The numerical simulation is of some advantages such as low cost, quick speed, complete data, and being capable to change parameters and simulate both the real and the ideal conditions. The non-linear partial differential (integral) equation often needs to be solved in the computational hydraulics, however, the mathematical theory of its numerical method has yet not been sufficient at present, for example, the development of some theories such as the strict stability analysis, error estimate and convergence proof has been unable to keep pace with the development of the numerical simulation method. Therefore, it is necessary to test the computation results and further modify the computation method based on the comparison between the nume
rical solution of the linear or non-linear equations and the trial value or the computation results of typical cases. The two-dimension numerical simulation of water flow and water quality is required in most practical projects, in which two types are
divided according to the different water area. In the broad water area such as the coast, estuary, lake and large reservoir, the horizontal scale is far larger than the vertical scale, the variation of the hydraulic parameters (for instance, velocity) is smaller in the vertical direction than in the horizontal direction, and the flow field can be expressed as the average flowing value along the water depth, so the two-dimension plane numerical simulation technique of the water flow and water quality may be used in these areas. But, in some other water areas such as the narrow-deep tide channel, narrow-deep valley reservoir and narrow-deep estuary area, the vertical varies of parameters such as the velocity, temperature and concentration is far greater than that in the horizontal direction, so the two-dimension vertical numerical simulation technique should be adopted in these areas.
The main contents in this paper are as follows: simulated irregular boundary and researched varies of algorithms; studied the two-dimension plane and two-dimension vertical computing; simulated water flow filed and water quality; researched algorithms and it's examples and combined with GIS tech. In two-dimension plane model the boundary fitted coordinate transformation were adopted to fit in with irregular area shape. In two-dimension vertical model a coordinate transformation were used to following the tracks of free surface and irregular topography boundary. In this paper, in order to fit in with the irregular boundary shape, and further to realize the quick forecast of the water flow and water quality, the two-dimension plane water flow and water quality model were built, in which the techniques of the coordinate transformation, alternate direction, staggered grid, linear recursion and implicit schemes were included. As for the reservoir of which the width to depth ratio is smaller, and the vertical variation of the parameters related such as the velocity, temperature and concentration is greater than their variation in the horizontal direction, the vertical two- dimension water flow model were set up, in which the water depth may change with the time and space, and when the coordinate transformation was made along the vertical direction, the variable area could be changed
into the regular one. In this model, the operator splitting technique, staggered grid technique and finite difference method were used.
The two-dimension plane model and two-dimension vertical model developed have been used in practical projects, in which:
In the flood simulation of Daqing City, the two-dimension plane unsteady flow model has been built to fit in with the irregular river boundary and realize the rapid forecast of the water flow. In this model, the geographic information sys
引文
[1] 周雪漪编著．计算水力学．清华大学出版社，1995
[2] 冯民权．淹没式温排水试验及三维紊动扩散数值模拟．武汉水利电力学院研究生论文，1988
[3] 郑邦民，赵昕．计算水动力学[M]．武汉：武汉大学出版社，2002年8月
[4] 朱家鲲．计算流体力学．科学出版社，1985
[5] 吴江航，韩庆书．计算流体力学的理论、方法及应用．科学出版社，1988
[6] 曹志芳．洪泛区水沙调度及洪灾损失快速评估模型研究．武汉大学博士学位论文，2001
[7] 黄东．明渠非恒定流数值模拟应用研究．武汉水利电力大学博士论文，1998
[8] 郑国栋．河渠非恒定流数值模拟研究与应用．武汉水利电力大学硕士学位论文
[9] 谭维炎著．计算浅水动力学—有限体积法的应用．清华大学出版社，1998
[10] 陶文铨．数值传热学．西安：西安交通大学出版社，1988
[11] 陶文铨．计算传热学的最近进展．西安：西安交通大学出版社
[12] Patanker,S.V.著，张政译．传热与流体流动的数值计算．科学出版社，1992
[13] C. J. Chen. The finite analytic method in flows and heat transfer. The university of Iowa, USA, 1987
[14] 陈景仁著．流体力学及传热学．北京：国防工业出版社，1987
[15] 陈景仁(美)．湍流模型及有限分析法．上海交通大学出版社，1988
[16] 李炜．对流扩散方程与边界层层流流动的混合有限分析．武汉水利电力大学学报，1994
[17] 曹祖德，王运洪著．水动力泥沙数值模拟．天津大学出版社，1994
[18] 黄胜，卢启苗．河口动力学．北京：水利电力出版社，1995，6
[19] 郑邦民，槐文信，齐鄂荣编著．洪水水力学．湖北科学技术出版社，2000年5月
[20] 沈鹤芳．洞庭湖君山蓄洪区洪水演进规律．西安理工大学硕士学位论文
[21] 李健生主编．中国江河防洪丛书．总论卷．中国水利水电出版社，1999
[22] 中华人民共和国水利部．中国’98大洪水．中国水利水电出版社，1999
[23] Liggett, J.A., and Woolhiser, D.A.. Difference solutions of the shallow-water equation. J.Engrg.mech. Div., ASCE, 1967, 93(2):39-71
[24] Chow, V.T., and Ben—Zvi, A.. Hydrodynamic modelling of two-dimensional water flow. J.Hydr.Div., ASCE, 1973, 99(11):2023-2040
[25] Xanthopoulos, T., and koutitas, C.. Numerical simulation of a two-dimensional flood wave propagation due to dam failure. J.Hydr.Res., Delft, the Netherland, 1976,14(4): 321-331
[26] Katopodes, N.D., and strelkoff, R.. Computing two-dimensional dam break flood waves. J.Hydr.Drv., ASCE, 1978,104(9): 1269—1288
[27] Garcia, R., and kahawita, R. A.. Numerical solution of the St Venant equations with the Mccormack finite difference scheme. Int.J.Numer. Methods in Fluids, 1986(6): 507-527
[28] Fennema, R.J., and Chaudhry, M.H.. Explicit methods for 2D transient free—surface flows. J.Hydr. Engrg., ASCE, 1990,116(8): 1013-1034
[29] Ponce, V.M., Li, R.M., and Simons, D.B.. Applicability of kinematic and diffusion models. J.Hydr. Engrg., ASCE, 1978,104(3): 353-360
[30] Akan, A.O., and Yen, B.C.. Diffusion—wave flood routing in channel networks. J.Hydr.Div., ASCE, 1981,107(6): 719-731
[31] 井上和也．水路非定常流数值计算法．水工学八应用研究，1986
[32] 刘树坤等．永定河泛区洪水灾害的数值模拟研究．非工程防洪措施资料选编．水利部海河水利委员会科技情报处，1988
[33] 刘树坤，袁小勇．二维溃坝洪水数值模拟的研究．水利水电科学研究院科学研究论文集第29集
[34] 刘树坤等．小清河分洪区洪水演进的数值模拟．水科学进展，1991；3
[35] 程玉慧，陶建华．天津市永定新河防洪非工程措施的研究
[36] 周孝德等．滞洪区二维洪水演进及洪灾风险分析．西安理工大学学报，1996
[37] 章汉毅等．洪潮遭遇情况下的水动力学计算．武汉水利电力大学学报，2000，33(5)
[38] Akanbi A. A etal. Model for flood propagation on initially dry land. Journal of Hydraulic Engineering, 1988,114(7)
[39] 李大鸣等．河道洪水演进的二维水流数学模型．天津大学学报，1998，31(4)
[40] T.Tucciarelli and D.Termini. Finite-element modelling of flood plain flow. J.Hydraulic Engineering, 2000, 126 (6)
[41] 谭维炎，胡四一．二维浅水流动的一种普适的高性能格式．水科学进展，1991，2(3)：154-161
[42] 谭维炎，胡四一．二维浅水明流的一种二阶高性能格式．水科学进展，1992，3(2)：89-95
[43] Zhao DH, Shen HW, Tabios GQ, Lai Js, Tan Wy. Finite-volume two-dimensional unsteady-flow model for river basins. Journal of Hydraulic Engineering, ASCE 1994, 120(7): 863-883
[44] 胡四一，谭维炎．无结构网格上二维浅水流动的数值模拟．水科学进展，1995，6(1)
[45] 谭维炎，胡四一等．长江中游洞庭湖防洪系统水流模拟．水科学进展，1996，7(4)
[46] 程晓陶等．分蓄洪区洪水演进数值模型．自然灾害学报，1996,5(1)
[47] 王船海，李光炽．流域洪水模拟．水利学报，1996，3
[48] 徐高洪，雒文生．有限分析法在浅水湖泊水流水质计算中的应用．水文，1996(6)
[49] 陈小红．湖泊水库垂向二维水温分布预测．武汉水利电力学院学报，1992，25(4)
[50] 赵昕，郑邦民．浅水热对流传输问题的瞬态解．武汉水利电力学院，1984年2月
[51] 郑邦民，赵昕．湖泊型温度流的数值分析．水利学报，1986年第3期
[52] 郑邦民．温度流的有限元数值模型．武汉水利电力学院学报，1981(1)
[53] 李炜主编．环境水力学进展[M]．武汉：武汉水利电力大学出版社，1999：318-330
[54] 马生伟，蔡启铭．浅水湖泊风生流的迎风有限元数值模型研究
[J]．水科学进展，2000，11(1)：70-74
[55] 叶守泽．水库水环境模拟预测与评价，中国水利水电出版社，1994
[56] 陈小红，刘美南，林艳珊(1997)．水库垂向二维水质分布研究．水利学报，1997(4)，9-16
[57] 雒文生，周志军．水库垂直二维湍流与水温水质耦合模型[J]．水电能源科学，1997，15(3)：1-7
[58] 江春波，张庆海，高忠信．河道立面二维非恒定水温及污染物分布预报模型．水利学报，2000(9)，20-24
[59] 谈广鸣，夏军强，郑邦民．挟沙水流紊动扩散系数的研究．水利学报，1998(12)
[60] 博斯腾湖可持续水资源管理应用研究[R]．武汉水利电力大学，新疆巴音郭愣自治州水利局，1999：44-53
[61] 李桂芬，倪浩清．水库风壅水面及风成波的试验研究．水利水电科学研究院科学研究论文集第7集．水利出版社，1982
[62] 清华大学水力学教研组编．水力学(下册)[M]．北京：高等教育出版社，1981：273-286
[63] 陈凯麒，李平衡，密小斌．温排水对湖泊、水库富营养化影响的数值模拟[J]．水利学报．1999(1)：22-26
[64] K.Mahmood and V.Yevjevich. Unsteady flow in open channels. Water Resources Publications, 1975
[65] 余常昭．水环境中污染物扩散输移原理与水质模型．中国环境科学出版社，1989
[66] 傅国伟．河流水质数学模型及其模拟计算，中国环境科学出版社，1987
[67] 金忠青．N-S方程的数值解和紊流模型，南京：河海大学出版社，1989
[68] P.J.罗奇著．计算流体力学．北京：科学出版社，1983
[69] W.罗迪著．环境问题的紊流模型．北京：水利电力出版社，1985
[70] 郑邦民．三维紊流对流扩散的数值分析．中国学者论文集锦．中国力学学会编．北京：北京大学出版社，1991
[71] 陈小红．分层型水库水温水质模型预测研究．武汉水利电力学院博士学位论文，1991年6月
[72] 李玉梁．水环境中污染物扩散输移原理与水质模型．北京：中国环境科学出版社，1989
[73] Harleman D R F. Hydrothermal Analysis of Lakes and Reservoirs. J.Hyd. Div., ASCE,Vol.108, No.Hy3, 1982
[74] [美]H.B费希尔著．内陆及近海水域中的混合．清华大学水力学教研组译．北京：水利电力出版社，1987
[75] 金士博N著．水环境数学模型．北京：中国建筑工业出版社，1987
[76] Thompson J F. Boundary-fitted coordinate system for numerical solution of partial differential equations[J]. J Comput Phys, 1982,47:1-108
[77] 姚孟电厂白山水库冷却池温排水数值模拟．武汉水利电力学院，1984年5月
[78] Roache,R.J.著．钟锡昌，刘学宗译．计算流体动力学．科学出版社，1983
[79] 汪德灌，杨艳艳．边界拟和坐标法的应用．河海大学学报，1989年第8期
[80] 韩国其，汪德灌，许协庆．风生环流的准三维数值模拟．河海大学学报，1989年第3期
[81] 万庆等．洪水灾害系统分析与评估．科学出版社，1999
[82] 孙桂华，王善序，王金銮，刘金平编译．洪水风险分析制图实用指南．水利电力出版社，1992
[83] Arvind kumar etal．Flood flow simulation in low land areas by two—dimensional model．第33回水理讲演会论文集，1989
[84] 徐正凡等编．明渠非恒定流基本方程及其数值解．武汉水利电力学院，1984
[85] 渭河下游三门峡库区洪水灾害风险分析及灾损评估．水利部西北水利科学研究所，1995
[86] Aura, R.A., and Wang, J.D.. Two—dimensional flood routing on steep slopes. J.Hydr.Engrg.,ASCE, 1984,110(8): 1121-1135
[87] Hromadka Ⅱ,T.V., and Lai, C.. Solving the two-dimensional diffusion flow equations. Proc. of the spec. conf. Sponsored by the Hydr. Div. of the ASCE, America Society of Civil Engineers,New York,N.Y., 1985: 555-561.
[88] A.M.Wasantha Lal. Weighted implicit finite-volume model for overland flow .J.of Hydraulic Engineering,vol. 124,No.9,1998
[89] Cappalaere, B.. Accurate diffusive wave routing. J.Hydr. Engrg., ASCE, 1997, 123(3): 174-181
[90] 张金凤．渭河下游洪水模拟模型及洪灾损失随机评估模型的研究．西安理工大学硕士学位论文
[91] 傅湘，纪昌明．洪灾损失评估指标的研究．水科学进展，2000(4)
[92] 唐允吉等．三门峡渭河库区洪灾损失研究．西安理工大学，1999
[93] 马宗晋．自然灾害评估指标体系的研究．自然灾害学报，1993，2(3)：1-7
[94] 姜彤，许明柱等．洪灾易损性概念模式．中国减灾，1997，7(2)
[95] 施国庆．洪灾损失率及其确定方法探讨．水利经济，1990(2)
[96] 李翔等．我国灾害经济统计评估系统及其指标体系的研究．自然灾害学报，1993，(1)：5-15
[97] 刘树坤等．黄河滩区及分滞洪区风险分析和减灾对策．黄河水
利出版社，1997
[98] 陈秀万．洪水灾害损失评估系统——遥感与GIS技术应用研究．中国水利水电出版社，1999
[99] 陈静生．中国水环境重金属研究．中国环境科学出版社，1992
[100] 周孝德．河流中重金属污染物迁移转化数学模型．陕西科技出版社，1996
[101] 武汉水利电力学院水力学教研室编．水力计算手册．北京：水利电力出社，1983：612-635
[102] 谢任之编著．溃坝水力学．山东科学技术出版社，1993：20-100
[103] 孙国浩，吴媛娇．溃坝洪水计算概况．水文，1987(1)：54-59
[104] 林秉南．突泄坝过程线简化分析．清华大学学报，1980(1)：17-31
[105] 吴葱葱，郭洪宇．两种土坝溃坝数学模型的研究与比较．武汉水电大学学报，2000，33(4)：17-20
[106] 朱伯芳．库水温度估算．水利学报，1985(2)
[107] 陈永灿，张宝旭，李玉梁．密云水库垂向水温模型研究．水利学报，1998(9)
[108] 范乐年，柳新之．湖泊、水库和冷却池水温预报通用模型．水利水电科学研究论文集——第17集(冷却水)．北京：水电出版社，1984
[109] 刘少文，方子云．水库水温模型初步研究．人民长江，1987(6)
[110] 李怀恩，沈晋．一维垂向水库水温数学模型研究与黑河水库水温预测．陕西机械学院学报，1990(4)
[111] 殷学鹏，杨传智．垂向一维水质模型及在龙滩水库水质预测中的应用．水利水电环境，1991(1)
[112] 沈晋，沈冰，李怀恩，李长兴，周孝德著．环境水文学．合肥：安徽科学技术出版社，1992
[113] 李怀恩．水库水温和水质预测研究述评．陕西机械学院学报，1987(4)
[114] 蒋红．水库水温计算方法探讨．水力发电学报，1999(2)
[115] 蒋红．水库、湖泊一维水温数学模型研究及瀑布沟水库水温预报计算．水利水电环境，1991(1)
[116] 钱小蓉，廖红，顾恒岳等．水库水温预测模型研究．重庆大学学报，1997(3)
[117] 陈新民，夏佳，罗国煜．黄河下游悬河决口灾害的风险分析与评价．水利学报，2000(10)
[118] 白玉川，万艳春，黄本胜等．河网非恒定流数值模拟的研究进展．水利学报，2000(12)
[119] 苑希民，刘树坤，陈浩．基于人工神经网络的多泥沙洪水预报．水科学进展，1999,10(4)
[120] 马千程，闾国年，施毅．GIS支持下计算格网自动生成技术．水科学进展，1999，10(1)
[121] 谭维炎．浅水动力学的回顾和当代前沿问题，水科学进展，1999,10(3)
[122] 魏文礼，王德意，刘玉玲．正交曲线网格生成技术的一点改进．武汉水利电力大学学报，2000,33(2)
[123] 傅湘，纪昌明．洪水保险的若干问题探讨．武汉水利电力大学学报，2000,33(1)
[124] 仇劲卫，李娜，程晓陶等．天津市城区暴雨沥涝仿真模拟系统．水利学报，2000(11)
[125] 郝瑞霞，周力行，陈惠泉．天然河道湍浮力流动的数值模拟Ⅰ：用隔离法实现三维复杂边界的处理．水利学报，1999(5)
[126] 郝瑞霞，陈惠泉，周力行．天然河道湍浮力流动的数值模拟Ⅱ：数学模型及其在冷却水工程中的应用实例．水利学报，1999(6)
[127] Prandtl,L.et.al.著，郭永怀、陆士嘉译．流体力学概论．科学出版社，1981
[128] 郑邦民，黄克中，魏良琰．流体力学．水利电力出版社，1989
[129] 郑邦民．流体力学中一般曲线坐标下的张量．武汉水利电力学院，1988年8月
[130] 郑邦民，赵明登．三维一般曲线坐标变换、网格生成及流场、温度场计算．水利学报，1992(7)
[131] 何少苓，王连祥，熊涛．平面二维不恒定流数值计算的若干经验．水利学报，1986(9)
[132] 刘智，林秉南，何少苓．三角形网格在二维不恒定流计算中的应用．水利学报，1987(9)
[133] 汪德灌，陈新建，杨艳艳．长江口南支流场的平面二维有限差分模拟．河海大学学报，1987(5)
[134] 汪德灌，杨艳艳．人工粘性、动边界及其在河口二维水流计算中的应用．水动力学研究与进展，1988(3)
[135] 汪德灌．潮汐河口河床演变准恒定二维模式．水利学报，1989(8)
[136] 周建军．处理不规则边界的边界对称点法．水利学报，1988(7)
[137] 周建军，林秉南，王连祥．平面二维泥沙数学模型研究及应用．水利学报，1993(11)
[138] 忻孝康．Poisson方程网格生成方法的一点改进．水动力学研究与进展，1992，7(增刊)
[139] 耿兆铨．非恒定水流数值解空间失稳问题分析．水利学报，1985(7)
[140] 耿兆铨．二维非恒定流有限元计算模式．华东水利学院学报，1984(4)
[141] 史金松．计算二维非恒定流的守恒型有限元法．水利学报，
1987(10)
[142] 马铁犹．计算流体动力学．北京航空学院出版社，1986
[143] 程心一．计算流体动力学．科学出版社，1984
[144] 郑邦民．溢流体形的数值模拟．中国科学A辑，1985(3)
[145] 郑邦民，齐鄂荣．天然河道非恒定流条件下阻力系数变化规律的研究．水动力学研究与进展，1990年9月
[146] 韩国其，汪德灌，许协庆．天然水流三维数值模拟的进展．河海大学科技进展．1989(3)
[147] 韩国其，汪德灌．宽浅型明渠非恒定流的数值模拟．河海大学学报．1990(3)
[148] 杨国录．河流数学模型．海洋出版社，1993
[149] 曹志先．水沙两相立面二维数学模型及数值方法研究．武汉水利电力学院博士论文，1991
[150] 吉庆丰．蒙特卡罗方法及其在水动力学中的应用研究．武汉水利电力大学博士学位论文，1999
[151] Sam S.Y.Wang, Advances in computational experimentational of river flows and sedimentation processes, Mississippi, 1994
[152] John D. Anderson, JR. Computational Fluid Dynamics. McGraw—Hill, 2002
[2] 冯民权．淹没式温排水试验及三维紊动扩散数值模拟．武汉水利电力学院研究生论文，1988
[3] 郑邦民，赵昕．计算水动力学[M]．武汉：武汉大学出版社，2002年8月
[4] 朱家鲲．计算流体力学．科学出版社，1985
[5] 吴江航，韩庆书．计算流体力学的理论、方法及应用．科学出版社，1988
[6] 曹志芳．洪泛区水沙调度及洪灾损失快速评估模型研究．武汉大学博士学位论文，2001
[7] 黄东．明渠非恒定流数值模拟应用研究．武汉水利电力大学博士论文，1998
[8] 郑国栋．河渠非恒定流数值模拟研究与应用．武汉水利电力大学硕士学位论文
[9] 谭维炎著．计算浅水动力学—有限体积法的应用．清华大学出版社，1998
[10] 陶文铨．数值传热学．西安：西安交通大学出版社，1988
[11] 陶文铨．计算传热学的最近进展．西安：西安交通大学出版社
[12] Patanker,S.V.著，张政译．传热与流体流动的数值计算．科学出版社，1992
[13] C. J. Chen. The finite analytic method in flows and heat transfer. The university of Iowa, USA, 1987
[14] 陈景仁著．流体力学及传热学．北京：国防工业出版社，1987
[15] 陈景仁(美)．湍流模型及有限分析法．上海交通大学出版社，1988
[16] 李炜．对流扩散方程与边界层层流流动的混合有限分析．武汉水利电力大学学报，1994
[17] 曹祖德，王运洪著．水动力泥沙数值模拟．天津大学出版社，1994
[18] 黄胜，卢启苗．河口动力学．北京：水利电力出版社，1995，6
[19] 郑邦民，槐文信，齐鄂荣编著．洪水水力学．湖北科学技术出版社，2000年5月
[20] 沈鹤芳．洞庭湖君山蓄洪区洪水演进规律．西安理工大学硕士学位论文
[21] 李健生主编．中国江河防洪丛书．总论卷．中国水利水电出版社，1999
[22] 中华人民共和国水利部．中国’98大洪水．中国水利水电出版社，1999
[23] Liggett, J.A., and Woolhiser, D.A.. Difference solutions of the shallow-water equation. J.Engrg.mech. Div., ASCE, 1967, 93(2):39-71
[24] Chow, V.T., and Ben—Zvi, A.. Hydrodynamic modelling of two-dimensional water flow. J.Hydr.Div., ASCE, 1973, 99(11):2023-2040
[25] Xanthopoulos, T., and koutitas, C.. Numerical simulation of a two-dimensional flood wave propagation due to dam failure. J.Hydr.Res., Delft, the Netherland, 1976,14(4): 321-331
[26] Katopodes, N.D., and strelkoff, R.. Computing two-dimensional dam break flood waves. J.Hydr.Drv., ASCE, 1978,104(9): 1269—1288
[27] Garcia, R., and kahawita, R. A.. Numerical solution of the St Venant equations with the Mccormack finite difference scheme. Int.J.Numer. Methods in Fluids, 1986(6): 507-527
[28] Fennema, R.J., and Chaudhry, M.H.. Explicit methods for 2D transient free—surface flows. J.Hydr. Engrg., ASCE, 1990,116(8): 1013-1034
[29] Ponce, V.M., Li, R.M., and Simons, D.B.. Applicability of kinematic and diffusion models. J.Hydr. Engrg., ASCE, 1978,104(3): 353-360
[30] Akan, A.O., and Yen, B.C.. Diffusion—wave flood routing in channel networks. J.Hydr.Div., ASCE, 1981,107(6): 719-731
[31] 井上和也．水路非定常流数值计算法．水工学八应用研究，1986
[32] 刘树坤等．永定河泛区洪水灾害的数值模拟研究．非工程防洪措施资料选编．水利部海河水利委员会科技情报处，1988
[33] 刘树坤，袁小勇．二维溃坝洪水数值模拟的研究．水利水电科学研究院科学研究论文集第29集
[34] 刘树坤等．小清河分洪区洪水演进的数值模拟．水科学进展，1991；3
[35] 程玉慧，陶建华．天津市永定新河防洪非工程措施的研究
[36] 周孝德等．滞洪区二维洪水演进及洪灾风险分析．西安理工大学学报，1996
[37] 章汉毅等．洪潮遭遇情况下的水动力学计算．武汉水利电力大学学报，2000，33(5)
[38] Akanbi A. A etal. Model for flood propagation on initially dry land. Journal of Hydraulic Engineering, 1988,114(7)
[39] 李大鸣等．河道洪水演进的二维水流数学模型．天津大学学报，1998，31(4)
[40] T.Tucciarelli and D.Termini. Finite-element modelling of flood plain flow. J.Hydraulic Engineering, 2000, 126 (6)
[41] 谭维炎，胡四一．二维浅水流动的一种普适的高性能格式．水科学进展，1991，2(3)：154-161
[42] 谭维炎，胡四一．二维浅水明流的一种二阶高性能格式．水科学进展，1992，3(2)：89-95
[43] Zhao DH, Shen HW, Tabios GQ, Lai Js, Tan Wy. Finite-volume two-dimensional unsteady-flow model for river basins. Journal of Hydraulic Engineering, ASCE 1994, 120(7): 863-883
[44] 胡四一，谭维炎．无结构网格上二维浅水流动的数值模拟．水科学进展，1995，6(1)
[45] 谭维炎，胡四一等．长江中游洞庭湖防洪系统水流模拟．水科学进展，1996，7(4)
[46] 程晓陶等．分蓄洪区洪水演进数值模型．自然灾害学报，1996,5(1)
[47] 王船海，李光炽．流域洪水模拟．水利学报，1996，3
[48] 徐高洪，雒文生．有限分析法在浅水湖泊水流水质计算中的应用．水文，1996(6)
[49] 陈小红．湖泊水库垂向二维水温分布预测．武汉水利电力学院学报，1992，25(4)
[50] 赵昕，郑邦民．浅水热对流传输问题的瞬态解．武汉水利电力学院，1984年2月
[51] 郑邦民，赵昕．湖泊型温度流的数值分析．水利学报，1986年第3期
[52] 郑邦民．温度流的有限元数值模型．武汉水利电力学院学报，1981(1)
[53] 李炜主编．环境水力学进展[M]．武汉：武汉水利电力大学出版社，1999：318-330
[54] 马生伟，蔡启铭．浅水湖泊风生流的迎风有限元数值模型研究
[J]．水科学进展，2000，11(1)：70-74
[55] 叶守泽．水库水环境模拟预测与评价，中国水利水电出版社，1994
[56] 陈小红，刘美南，林艳珊(1997)．水库垂向二维水质分布研究．水利学报，1997(4)，9-16
[57] 雒文生，周志军．水库垂直二维湍流与水温水质耦合模型[J]．水电能源科学，1997，15(3)：1-7
[58] 江春波，张庆海，高忠信．河道立面二维非恒定水温及污染物分布预报模型．水利学报，2000(9)，20-24
[59] 谈广鸣，夏军强，郑邦民．挟沙水流紊动扩散系数的研究．水利学报，1998(12)
[60] 博斯腾湖可持续水资源管理应用研究[R]．武汉水利电力大学，新疆巴音郭愣自治州水利局，1999：44-53
[61] 李桂芬，倪浩清．水库风壅水面及风成波的试验研究．水利水电科学研究院科学研究论文集第7集．水利出版社，1982
[62] 清华大学水力学教研组编．水力学(下册)[M]．北京：高等教育出版社，1981：273-286
[63] 陈凯麒，李平衡，密小斌．温排水对湖泊、水库富营养化影响的数值模拟[J]．水利学报．1999(1)：22-26
[64] K.Mahmood and V.Yevjevich. Unsteady flow in open channels. Water Resources Publications, 1975
[65] 余常昭．水环境中污染物扩散输移原理与水质模型．中国环境科学出版社，1989
[66] 傅国伟．河流水质数学模型及其模拟计算，中国环境科学出版社，1987
[67] 金忠青．N-S方程的数值解和紊流模型，南京：河海大学出版社，1989
[68] P.J.罗奇著．计算流体力学．北京：科学出版社，1983
[69] W.罗迪著．环境问题的紊流模型．北京：水利电力出版社，1985
[70] 郑邦民．三维紊流对流扩散的数值分析．中国学者论文集锦．中国力学学会编．北京：北京大学出版社，1991
[71] 陈小红．分层型水库水温水质模型预测研究．武汉水利电力学院博士学位论文，1991年6月
[72] 李玉梁．水环境中污染物扩散输移原理与水质模型．北京：中国环境科学出版社，1989
[73] Harleman D R F. Hydrothermal Analysis of Lakes and Reservoirs. J.Hyd. Div., ASCE,Vol.108, No.Hy3, 1982
[74] [美]H.B费希尔著．内陆及近海水域中的混合．清华大学水力学教研组译．北京：水利电力出版社，1987
[75] 金士博N著．水环境数学模型．北京：中国建筑工业出版社，1987
[76] Thompson J F. Boundary-fitted coordinate system for numerical solution of partial differential equations[J]. J Comput Phys, 1982,47:1-108
[77] 姚孟电厂白山水库冷却池温排水数值模拟．武汉水利电力学院，1984年5月
[78] Roache,R.J.著．钟锡昌，刘学宗译．计算流体动力学．科学出版社，1983
[79] 汪德灌，杨艳艳．边界拟和坐标法的应用．河海大学学报，1989年第8期
[80] 韩国其，汪德灌，许协庆．风生环流的准三维数值模拟．河海大学学报，1989年第3期
[81] 万庆等．洪水灾害系统分析与评估．科学出版社，1999
[82] 孙桂华，王善序，王金銮，刘金平编译．洪水风险分析制图实用指南．水利电力出版社，1992
[83] Arvind kumar etal．Flood flow simulation in low land areas by two—dimensional model．第33回水理讲演会论文集，1989
[84] 徐正凡等编．明渠非恒定流基本方程及其数值解．武汉水利电力学院，1984
[85] 渭河下游三门峡库区洪水灾害风险分析及灾损评估．水利部西北水利科学研究所，1995
[86] Aura, R.A., and Wang, J.D.. Two—dimensional flood routing on steep slopes. J.Hydr.Engrg.,ASCE, 1984,110(8): 1121-1135
[87] Hromadka Ⅱ,T.V., and Lai, C.. Solving the two-dimensional diffusion flow equations. Proc. of the spec. conf. Sponsored by the Hydr. Div. of the ASCE, America Society of Civil Engineers,New York,N.Y., 1985: 555-561.
[88] A.M.Wasantha Lal. Weighted implicit finite-volume model for overland flow .J.of Hydraulic Engineering,vol. 124,No.9,1998
[89] Cappalaere, B.. Accurate diffusive wave routing. J.Hydr. Engrg., ASCE, 1997, 123(3): 174-181
[90] 张金凤．渭河下游洪水模拟模型及洪灾损失随机评估模型的研究．西安理工大学硕士学位论文
[91] 傅湘，纪昌明．洪灾损失评估指标的研究．水科学进展，2000(4)
[92] 唐允吉等．三门峡渭河库区洪灾损失研究．西安理工大学，1999
[93] 马宗晋．自然灾害评估指标体系的研究．自然灾害学报，1993，2(3)：1-7
[94] 姜彤，许明柱等．洪灾易损性概念模式．中国减灾，1997，7(2)
[95] 施国庆．洪灾损失率及其确定方法探讨．水利经济，1990(2)
[96] 李翔等．我国灾害经济统计评估系统及其指标体系的研究．自然灾害学报，1993，(1)：5-15
[97] 刘树坤等．黄河滩区及分滞洪区风险分析和减灾对策．黄河水
利出版社，1997
[98] 陈秀万．洪水灾害损失评估系统——遥感与GIS技术应用研究．中国水利水电出版社，1999
[99] 陈静生．中国水环境重金属研究．中国环境科学出版社，1992
[100] 周孝德．河流中重金属污染物迁移转化数学模型．陕西科技出版社，1996
[101] 武汉水利电力学院水力学教研室编．水力计算手册．北京：水利电力出社，1983：612-635
[102] 谢任之编著．溃坝水力学．山东科学技术出版社，1993：20-100
[103] 孙国浩，吴媛娇．溃坝洪水计算概况．水文，1987(1)：54-59
[104] 林秉南．突泄坝过程线简化分析．清华大学学报，1980(1)：17-31
[105] 吴葱葱，郭洪宇．两种土坝溃坝数学模型的研究与比较．武汉水电大学学报，2000，33(4)：17-20
[106] 朱伯芳．库水温度估算．水利学报，1985(2)
[107] 陈永灿，张宝旭，李玉梁．密云水库垂向水温模型研究．水利学报，1998(9)
[108] 范乐年，柳新之．湖泊、水库和冷却池水温预报通用模型．水利水电科学研究论文集——第17集(冷却水)．北京：水电出版社，1984
[109] 刘少文，方子云．水库水温模型初步研究．人民长江，1987(6)
[110] 李怀恩，沈晋．一维垂向水库水温数学模型研究与黑河水库水温预测．陕西机械学院学报，1990(4)
[111] 殷学鹏，杨传智．垂向一维水质模型及在龙滩水库水质预测中的应用．水利水电环境，1991(1)
[112] 沈晋，沈冰，李怀恩，李长兴，周孝德著．环境水文学．合肥：安徽科学技术出版社，1992
[113] 李怀恩．水库水温和水质预测研究述评．陕西机械学院学报，1987(4)
[114] 蒋红．水库水温计算方法探讨．水力发电学报，1999(2)
[115] 蒋红．水库、湖泊一维水温数学模型研究及瀑布沟水库水温预报计算．水利水电环境，1991(1)
[116] 钱小蓉，廖红，顾恒岳等．水库水温预测模型研究．重庆大学学报，1997(3)
[117] 陈新民，夏佳，罗国煜．黄河下游悬河决口灾害的风险分析与评价．水利学报，2000(10)
[118] 白玉川，万艳春，黄本胜等．河网非恒定流数值模拟的研究进展．水利学报，2000(12)
[119] 苑希民，刘树坤，陈浩．基于人工神经网络的多泥沙洪水预报．水科学进展，1999,10(4)
[120] 马千程，闾国年，施毅．GIS支持下计算格网自动生成技术．水科学进展，1999，10(1)
[121] 谭维炎．浅水动力学的回顾和当代前沿问题，水科学进展，1999,10(3)
[122] 魏文礼，王德意，刘玉玲．正交曲线网格生成技术的一点改进．武汉水利电力大学学报，2000,33(2)
[123] 傅湘，纪昌明．洪水保险的若干问题探讨．武汉水利电力大学学报，2000,33(1)
[124] 仇劲卫，李娜，程晓陶等．天津市城区暴雨沥涝仿真模拟系统．水利学报，2000(11)
[125] 郝瑞霞，周力行，陈惠泉．天然河道湍浮力流动的数值模拟Ⅰ：用隔离法实现三维复杂边界的处理．水利学报，1999(5)
[126] 郝瑞霞，陈惠泉，周力行．天然河道湍浮力流动的数值模拟Ⅱ：数学模型及其在冷却水工程中的应用实例．水利学报，1999(6)
[127] Prandtl,L.et.al.著，郭永怀、陆士嘉译．流体力学概论．科学出版社，1981
[128] 郑邦民，黄克中，魏良琰．流体力学．水利电力出版社，1989
[129] 郑邦民．流体力学中一般曲线坐标下的张量．武汉水利电力学院，1988年8月
[130] 郑邦民，赵明登．三维一般曲线坐标变换、网格生成及流场、温度场计算．水利学报，1992(7)
[131] 何少苓，王连祥，熊涛．平面二维不恒定流数值计算的若干经验．水利学报，1986(9)
[132] 刘智，林秉南，何少苓．三角形网格在二维不恒定流计算中的应用．水利学报，1987(9)
[133] 汪德灌，陈新建，杨艳艳．长江口南支流场的平面二维有限差分模拟．河海大学学报，1987(5)
[134] 汪德灌，杨艳艳．人工粘性、动边界及其在河口二维水流计算中的应用．水动力学研究与进展，1988(3)
[135] 汪德灌．潮汐河口河床演变准恒定二维模式．水利学报，1989(8)
[136] 周建军．处理不规则边界的边界对称点法．水利学报，1988(7)
[137] 周建军，林秉南，王连祥．平面二维泥沙数学模型研究及应用．水利学报，1993(11)
[138] 忻孝康．Poisson方程网格生成方法的一点改进．水动力学研究与进展，1992，7(增刊)
[139] 耿兆铨．非恒定水流数值解空间失稳问题分析．水利学报，1985(7)
[140] 耿兆铨．二维非恒定流有限元计算模式．华东水利学院学报，1984(4)
[141] 史金松．计算二维非恒定流的守恒型有限元法．水利学报，
1987(10)
[142] 马铁犹．计算流体动力学．北京航空学院出版社，1986
[143] 程心一．计算流体动力学．科学出版社，1984
[144] 郑邦民．溢流体形的数值模拟．中国科学A辑，1985(3)
[145] 郑邦民，齐鄂荣．天然河道非恒定流条件下阻力系数变化规律的研究．水动力学研究与进展，1990年9月
[146] 韩国其，汪德灌，许协庆．天然水流三维数值模拟的进展．河海大学科技进展．1989(3)
[147] 韩国其，汪德灌．宽浅型明渠非恒定流的数值模拟．河海大学学报．1990(3)
[148] 杨国录．河流数学模型．海洋出版社，1993
[149] 曹志先．水沙两相立面二维数学模型及数值方法研究．武汉水利电力学院博士论文，1991
[150] 吉庆丰．蒙特卡罗方法及其在水动力学中的应用研究．武汉水利电力大学博士学位论文，1999
[151] Sam S.Y.Wang, Advances in computational experimentational of river flows and sedimentation processes, Mississippi, 1994
[152] John D. Anderson, JR. Computational Fluid Dynamics. McGraw—Hill, 2002