用户名: 密码: 验证码:
基于GIS的二维水沙数学模型及其在黄河下游应用的研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
针对黄河下游河道水沙运动的特点,以“IT”技术为平台,研制了基于GIS的黄河下游二维水沙数学模型。
     通过对现有国内外水沙模型的比较,选取了适合多沙河流水沙特点的方程离散格式及数值方法:采用Osher及LSS格式计算对流项水沙数值通量,借助高斯积分确定扩散项中控制体形心处的梯度矢量,通过预测、校正二步格式求解方程组,以保证水流计算格式的解具有良好的物理基础,符合解的特征构造,并具有守恒性、逆风性、单调保持性和对间断的高分辨率等优良特性。对于由水压力与底坡项不平衡所导致的动量不守恒问题,采用三点高斯积分法来计算各边的平均法向通量。本文对非均匀沙沉速、床沙级配调整、横向冲刷坍塌变形等关键技术问题进行了深入分析,研制的二维水沙模型不但适用于低含沙水流的计算,而且适用于多沙河流的模拟计算。通过对ArcGIS软件的二次开发,设计和开发出了基于GIS的二维数学模型的前后处理系统。基于GIS的前后处理系统,不仅极大的提升了前期工作效率,而且在地形高程的生成精度、信息查询速度和可视化程度方面有了明显提高。
     为检验所建模型的合理性和精度,利用国内外试验资料对其进行了测试,测试案例包括恒定流、非恒定流、溃坝、压缩水流、弯道、复式河道及复杂边界等不同情况。此外,利用黄河下游不同时期的实测资料和专题测验资料,对泥沙沉速模块、动床阻力模块、分组挟沙力模块等进行了测试,在此基础上,又将所建立的模型应用于黄河下游实测调水调沙洪水的模拟计算,结果表明,计算精度满足防洪调度的实际要求。
A 2D water and sediment numerical model in the lower Yellow River was developed aimed at flow and sediment transport characteristics and based on GIS platform.By comparing the development of numerical modeling technologies in the world, equation discrete schemes and methodologies suited for hyper-concentration flow are selected .The Osher and LSS schemes are used to calculate the convection term of water and sediment. The Gauss theory is adopted to induce the control volume central gradient of the diffusion term. Solutions of the equations with forecasting and adjusting two-steps method may ensure results have real physical means and fit well with characteristic construction of solution, which has flux conservation, upwind, monotonic keeping and high resolving power for discontinuity . As for the unbalance momentum induced from water pressure non-uniform with river bed slope, the Gauss theory is used to calculate averaged flux of each side, which not only avoids the error flow phenomena for steep water depth gradient, but also shows the distribution of status vectors correctly. The paper studies the key problems including non-uniform sediment settling velocity, suspend sediment transport capacity, bed sediment composition adjustment, movable bed roughness variety, turbulence viscous coefficient, sediment viscous coefficient and transverse bank scour simulating. By using ArcGIS and secondary developing with GIS, the model's pre-processes and post processes include grid generation system, digital elevation generation system and flow field visualization system.A typical testing procedure is used to test the flow module and sediment module, including processes of steady and unsteady flow, dam break flow, impressed flow, compounded river channel and complex boundary conditions. Data of conventional measurement in the Lower Yellow River and special observations for this study are used to calibrate the sediment settling velocity, movable bed roughness coefficient and gradations sediment transport capacity etc. Sets of scenario studies are carried out. Several typical flood processes are verified and results prove that the model can be satisfied for the flood management in the Lower Yellow River.The innovation achievements of the paper may be summarized as follows: Based on numerical method and sediment transport law with developed IT technology, realized the combination of IT technology with water science; Numerical methodologies of hyper-concentration flow and research on key problems of hyper-concentration flow fundamental theory improve the modeling precision for hyper-concentration flow. Based on GIS platform, not only improve the efficiency but also improve the precision of digital elevation generation, information searching speed and visualization technology obviously.
引文
[1] 李国英,治理黄河思辨与践行,郑州,黄河水利出版社,2003.8
    [2] 芮孝芳,产汇流理论,北京:水利电力出版社,1995
    [3] 朱庆平等,关于黄河数学模型系统建设的思考,黄河报,2004(6).
    [4] Ruth Flower, Cost and Quality of Higher Education--Reaction to the National Commission on the Cost of Higher Education. The Education Digest 1998, October, Vol. 64, No.2.
    [5] Price R K, Ahmad K, Holz P. Hydroinformatics concepts. Hydroinformatics tools, Marsake, J. et al(eds.), Kluwer Academic Publishers, Netherlands, 1998. 47~76
    [6] Thompson G. Hydroinformatics: a cornerstone to information technology in business. Hydroinformatics'98, Babovic & Larsen(eds.), Balkema, Rotterdam, 1998. 9~14
    [7] De Hoop S., Van Der Meij L., Van Hekken M., and Vijlbrief T., Integrated 3D Modeling within a GIS[A]. In: Proceedings of the 6th International GIS Workshop "Advanced Geographic Data Modeling", AGDM'94[C], 1994, Delft
    [8] 王光谦,张仁,惠遇甲,王桂仙,王兴奎.清华大学泥沙回顾.泥沙研究,1999,(6)
    [9] 胡春宏,周文浩.中国水科院十年来泥沙研究综述.泥沙研究,1999,(6)
    [10] 李义天,谢鉴衡,王明甫.武汉水利电力大学泥沙研究进展.泥沙研究,1999,(6)
    [11] 汤立群,唐洪武,陈国祥.河海大学泥沙研究进展与瞻望.泥沙研究,1999,(6)
    [12] 朱自强,等著.应用计算流体力学.北京航空航天大学出版社.2001.1
    [13] 刘儒勋,舒其望.著.计算流体力学的若干新方法.科学出版社.2003.4.
    [14] 韩其为,何明民.泥沙数学模型中冲淤计算的几个问题,水利学报,1988,(5)
    [15] 韩其为,何明民.水库淤积与河床演变的数学模型.泥沙研究,1987,(3)
    [16] 韩其为,何明民,陈永灿.恢复饱和系数研究及挟沙能力的初步分析,中国水利水电科学研究院,1995
    [17] 韩其为,何明民.恢复饱和系数初步研究,泥沙研究,1997,(3)
    [18] 韩其为.水库不平衡输沙的初步研究,水库泥沙报告汇编,黄河泥沙研究协调小组编印,1973
    [19]韦直林.天然河道平面二维水沙数学模型研究.武汉水利电力大学学报,1992,(6)
    [20] 张世奇。二维动边界潮汐输沙及河床变形计算[J]。泥沙研究,1988,(4).
    [21] 刘树坤,宋玉山,程晓陶,等.黄河滩区及分滞洪区风险分析和减灾对策,黄河水利出版社,1999.
    [22] 张红武,姚文艺,张俊华,钱意颖.黄科院90年代泥沙研究进展.泥沙研究,1999,(6)
    [23] 李炜。不可压缩常粘性流体运动的数值模拟及混合有限分析解法,武汉水利电力大学研究生教材,1994.
    [24] 夏军强.河岸冲刷机理研究及数值模拟.清华大学水利水电系博士论文,2002.4
    Ikeda S, Parker G, Sawai K. Bend theory of river meanders, Part Ⅰ:linear development[J].Journal of Fluid Mechanics, 112: 363-377.
    [25] Hasegawa K. Universal bank erosion coefficient for meandering rivers[J], Journal of HydraulicEngineering, ASCE, 1989, 115(6): 744-765.
    [26] Engelund EFlow and bed topography in channel bends[J]. Journal of Hydraulic Div. ASCE, 100, 1974, (11): 1631-1648.
    [27] Odgaard A J. Meander flow model Ⅰ:development[J]. Journal of Hydraulic Engineering, ASCE, 1986, 112(12):1117-1136.
    [28] Shimizu Y, Yamaguchi H, Itakura T. Three dimensional computation of flow and bed deformation [J]. Journal of Hydraulic Engineering, ASCE, 1990, 116(9): 1090-1108.
    [29] 夏军强,王光谦.考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟.水利学报,2002,(6)
    [30] 1955
    [31] 韩其为,何明民.底层泥沙交换和状态概率及推悬比研究,水利学报,1999,(10)
    [32] 韩其为,何明民.论非均匀悬移质二维不平衡输沙方程及其边界条件,水利学报,1997(1)
    [33] D. H. Zhao, H. W. Shen, et al. Finite-Volume Two-Dimension Unsteady-Flow Model For River Basins. Journal of the Hydraulic Engineering, Vol. 120, No. 7, P863~883, 1994.
    [34] 谭维炎著.计算浅水动力学—有限体积法的应用.清华大学出版社.1998.10
    [35] Dartzi Pan, Jen-Chien Cheng. A Second-Order Upwind Finite-Volume Method For The Euler Solution On Unstructured Triangulaur Meshes. Intemational Joumal for Numerical Method in Ruids, Vol. 16, 1079~1098,1993.
    [36] De-yu Zhong, Hong-Wu Zhang, A Two Dimensional Numerical Model for Simulation of Flood Propagation in Initially Dry Bed. Rood Defence' 2002 Wu et al(eds)(?)2002 Science Press, New York Ltd., ISBN 1-880132-54-0.
    [37] Ming-Hseng Tseng. Explicit Finite Volume Nin-Oscillatory Schemes For 2d Transient Free-Surface Flows. Intemational Journal for Numerical Method Fluids. VOL. 30. p831~843. 1999
    [38] 万洪涛,周虎成,吴应湘,万庆.黄河下游花园口~夹河滩河段二维洪水模拟.水科学进展.Vol.13,No.2,P215~222,2002.3.
    [39] Nagata, N. Hosoda, T. Muramoto Y.. Numerical analysis of river channel processes with bank erosion Journal of Hydraulic Engineering, Vol. 126 No. 4, 2000 2413~2431. 1997
    [40] Mosselman, E. Disscussion of "river width adjustment I: processes and mechanisms", Journal of Hydraulic Engineering, Vol. 126 No. 2, 2000
    [41] Sangmook Shin. Reynolds-Averaged Navier-Stokes Computation of Tip Clearance Flow In a Compressor Cascade Using an Unstructured Grid. Blacksburg, Virginia. 2001.
    [42] K. Willcox and J. Peraire. Aeroelastic Computations in the Time Domain Using Unstructured Meshes. International Journal for Numerical Method Fluids. VOL. 40. p2413~2431. 1997
    [43] Engelund, Eand E, Hansen. Amonograph on Sediment Transport in Alluvial Streams. Teknisk Forlag, Copenhagen, 1972, p. 62.
    [44] Osman A M and Thorne C R. Riverbank stability analysis Ⅰ: Theory. ASCE, Journal of Hydraulic Engineering. 1998. 114(2)
    [45] Osman, A. M., and Thorne, C. R. Riverbank Stability Analysis Ⅰ: Theory & Ⅱ Application. ASCE, Joumal of Hydraulic Engineering, 1988, 114(2): 134-172.
    [46] Stanley Osher and Fred Solomon. Upwing Difference Schemes for Hyperbolic Systems Of Conservation Laws. Mathematics of Computation. Vol38. NO. 158. 339~374. 1982
    [47] Francisco Alcrudo. Flux Difference Splitting For 1d Open Channel Flow Equations. International Journal for Numerical Method Fluids. VOL. 14.1009~1018. 1992
    [48] J. S. Mathur. The Simulation Of In viscid Compressible Flows Using An Upwind Kinetic Method on Unstructured Grids. International Journal for Numerical Method Fluids.VOL. 15. p59~82. 1992
    [49] I. Demirdzic. A Collocated Finite Volume Method For Predicting Flows At All Speeds. International Journal for Numerical Method Fluids. VOL. 16. 1029~1050. 1993
    [50] Tan, W. Y, Shallow water hydrodynamics. Elservier, Amsterdam, the Netherlands. 1992.
    [51] K. Anastasiou, C. T. Chan. Solution of the 2D Shallow Water Equations Using the Finite Volume Methods on Unstructured Triangular Meshes. International Journal for Numerical Method in Fluids, Vol. 24, 1225~1245, 1997.
    [52] Lin G-F.; Lai J. -S.; Guo W. -D. Finite-volume component-wise TVD schemes for 2D shallow water equations, Advances in Water Resources, August 2003, vol. 26, no. 8, pp. 861-873(13)
    [53] E. Erwin Holly, Numerical Calculation of Flow Field, World Scientific, 1984.
    [54] Biker E. W. Littoral Drift as Function of Waves and Current[C]. Proc. 11th Conf. Coastal Engineering, 1968.
    [55] 钱意颖,等.编著.黄河泥沙冲淤数学模型,黄河水利出版社,1998.8
    [56] 黄金池,万兆惠.黄河下游河床平面变形模拟研究.水利学报,1999,(2)
    [57] 王尚毅,启元棪.黄河“揭底冲刷”问题的初步研究,泥沙研究,1982,(2)
    [58] 万兆惠,宋天成.“揭河底”冲刷现象的分析(齐璞,赵文林,杨美卿编),黄河高含沙水流运动规律及应用前景,科学出版社,1998
    [59] 李义天,等.著.河道平面二维水沙数学模型.中国水利水电出版社,2001
    [60] 涂启华,安催花,李世滢,余欣,等,著.《水利动能手册(泥沙分册)》,中国水利水电出版社(待出版),2004
    [61] 芮孝芳,由流路长度分布律和坡度分布律确定地貌单位线,水科学进展,P548~550,2005.5
    [62] 芮孝芳,石朋.基于地貌扩散和水动力扩散的流域瞬时单位线研究,水科学进展,439-4442002,13(4):
    [63] Wu W M, Rodi and W and Wenka T. 3D Numerical Modeling of Flow and Sediment and Sediment Transport in Open Channels. ASCE. Journal of Hydraulic Engineering, 2000, 126(1)
    [64] Bagnold, R. A.. The Noture of Sediment and of Bed Load Transpat in Water, Proc. Rayal Society Ser. A. Val. 332, 1973
    [65] 张瑞瑾,编,河流动力学,中国工业出版社,1961
    [66] 余明辉,杨国录,等.非均匀沙水流挟沙力公式的初步研究.泥沙研究.P25~28,2001.6
    [67] 曹如轩.高含沙水流挟沙能力初步研究,水利水电技术,1979,(5)。
    [68] Yang, C. T. Sediment Transport-Theory and Practice, The McGraw-Hill Companies, Inc. 1996
    [69] 舒安平.水流挟沙能力公式的检验与评述,人民黄河,1993,(1)
    [70] 吴保生,龙毓骞.黄河输沙能力的对比分析,全国泥沙基本理论研究学术讨论会论文集,中国水利学会泥沙专业委员会,1992
    [71] 王士强.黄河泥沙冲淤数学模型研究.水科学进展,1996 (9)
    [72] 张红武,张清.黄河水流挟沙能力的计算公式,人民黄河,1992,(11)
    [73] 杨国录.著.河流数学模型.武汉水利电力学院出版社.1992
    [74] VALERIO CALEFFI, ALESSANDRO VALIANI, ANDREA ZANNI. Finite volume method for simulating extreme flood events in natural channels. Journal of Hydraulic Research Vol. 41, No. 2(2003), p167~177.[75] Wu W M, Rodi and W and Wenka T. 3D Numerical Modeling of Flow and Sediment and Sediment Transport in Open Channels. ASCE. Journal of Hydraulic Engineering, 2000, 126(1)
    [76] Thorne C R and Osman A M. Riverbank stability analysis Ⅱ: Application. ASCE, Journal of Hydraulic Engineering. 1998. 114(2)
    [77] 张书农,唐存本.水槽试验的侧壁影响问题[J].高等学校自然科学学报.土木、建筑、水利版,Vol.1,No.4,1965.10.
    [78] Van Rijin L C. Sediment Transport. Part Ⅰ: Bedload Transport. ASCE. Journal of Hydraulic Engineering, 1984. 110(10)
    [79] Klingeman, P. C., Kehe, S. M., and Owusu, Y. A. (1984): Streambank Erosion and Channel Scour Manipulation Using Pockfill Dikes and Gabions, WRRI Report for Project No. 373909, Oregon State Univ., Corvallis, Oregon.
    [80] 闵卫东,唐泽圣.二维任意域内点集的Delaunay三角剖分的研究[J].计算机学报,1995,18(5):357-364.
    [81] 武晓波,王世新,肖春生.Delaunay三角网的生成算法研究.测绘学报,1999,28(1):28-35.
    [82] 普雷帕拉塔FP,沙莫斯MI.计算几何导论[M].北京:科学出版社,1990
    [83] Sibson R. Locally Equiangular T triangulates ions. Computer Journal [J]. 1978, 21 (3):243-245.
    [84] Shamos M I. and Hoey D. Closest-point Problems[C]. In: Proceedings of the 16th Annual Symposium on the Foundations of Computer Science, 1975. 151-162
    [85] Lewis B A. and Robinson J S. Triangulation of Planar Regions with Applications [J]. The Computer Journal, 1978, 21(4): 324-332.
    [86] Lee D T. and Schachter B J. Two Algorithms for Constructing a Delaunay Triangulation [J]. Int. J. of Computer and Information Sciences, 1980, 9(3).
    [87] Tsai V J D. Delaunay Triangulations in T IN Creation: an Overview and a Lineartime Algorithm[J]. Int. J. of GIS, 1993, 7(6): 501-524.
    [88] Law son C L. Software for C'Surface Interpolation [J]. Mathematical Software Ⅲ. J. Rice, Ed. New York: Academic Press, 1977.
    [89] Bowyer A. Computing Dirichlet Tessellations [J]. Computer Journal, 1981(24): 162-166.
    [90] Watson D F. Computing the n-dimension Delaunay Tesselation with Application to Voronoi Polytopes [J]. Computer Journal, 1981 (24): 167-172.
    [91] Sloan S W. A fast Algorithm for Constructing Delaunay Triangulations in the Plane [J]. Advanced Engineering Software, 1987(9): 34-55.
    [92] Macedonia G and Pareschi M T. An Algorithm for the Triangulation of Arbitrarily Distributed Points: Applications to Volume Estimate and Terrain Fitting, Computers & Geosciences, 1991(17):859-874.
    [93] De Floriani L and Puppo E. An on-line Algorithm for Constrained Delaunay Triangulation[J]. CV GIP: Graphical Models and Image Processing, 1992, 54(3): 290-300.
    [94] Tsai V J D. Delaunay Triangulations in T IN Creation: an Overview and a Lineartime Algorithm [J]. Int. J. of GIS, 1993, 7(6): 501-524.
    [95] Green P J. and Sibson R. Computing Dirich let Tesselations in the Plane [J]. The Computer Journal, 1978, 21(2): 168-173.
    [96] Brassel K E. and Reif D. Procedure to Generate Thiessen Polygons [J]. Geophysical Analysis, 1979, (11): 289-303.
    [97] McCaullagh M T. and Ross C G Delaunay Triangulation of a Random Data Set for Iirarithmic Mapping [J]. The Cartographic Journal, 1980(17): 93-99.
    [98] Maus A. Delaunay Triangulation and the Convex Hull of n Points in Expected Linear Time [J]. BIT, 1984(24): 151-163.
    [99] 胡恩球等.有限元网格生成方法发展综述[J].计算机辅助设计与图形学学报,1997,9(4):378-383.
    [100] 胡恩球.电磁器件有限元仿真及超导单极电机电磁场数值分析[学位论文].武汉:武汉华 中理工大学,1995.
    [101] 管镭,孟宪琦,魏生民.Delaunay三角网格化的算法及实现[J].西北工业大学学报,1996,14(1):138-142.
    [102] 唐泽圣,徐志强.二维点集三角剖分的动态生成与修改[J].计算机辅助设计与图形学报,1990,3:1-8.
    [103] 周晓云,朱心雄.散乱数据点三角剖分方法综述[J].工程图学学报,1993,1:48-54.
    [104] 芮孝芳,朱庆平.分布式流域水文模型研究中的几个问题,水利水电科技进展,2002(22).
    [105] 马正飞,殷翔.著.数学计算方法与软件的工程应用.p111~p115.化学工业出版社.2002.12.
    [106] Wu W M, Rodi W, Wenka T, 3D Numerical modeling of flow and sediment and sediment transport in open channels[J], Journal of Hydraulic Engineering, ASCE, 2000, 126(1): 4-15.
    [107] Molls T, Chaudhry M H. Depth averaged open channel flow model[J], Journal of Hydraulic Engineering, ASCE, 1995, 121(6): 453-456.
    [108] Lien H C, et al. Bend flow simulation using 2-d depth-averaged model[J], Journal of Hydraulic Engineering, ASCE, 1999, 125(10): 1097-1106.
    [109] Teision Ch, Latteux B. A Depth-Integrated Bidimensional Model of Suspended Sediment Transport[C]。Proc. of the 3rd Int. Symp. on River sedimentation, The Univ. of Mississippi, 1986.
    [110] Christopher G. Koutitas. Mathematical Models in Coastal Engineering[M]。Pentech Press, 1988.
    [111] Sui Liang HUANG and Onyx W. H. Wai, A review of heavy metal pollution in the Pearl River estuary, Journal of Hydrodynamics, Ser.B, Vol. 16, No4. pp. 367-378.
    [112] Sui Liang Huang (2003), Adsorption of cadmium ions onto the Yellow River sediments, Water Quality Research Journal of Canada, Vol. 38, No. 2. pp413~432.
    [113] NCCHE, Verification of NCCHE 3d numerical model with Prescribed Forcing Solution Technicque, 1998. 10
    [114] 谢鉴衡,魏良琰.河流泥沙数学模型的回顾与展望.泥沙研究,1987
    [115] Cole P, Miles G V. Two Dimensional of Mud Transport[J]。ASCE, 1983, 109(1).
    [116] Hurther, D.and Lemmin, U.. Discussion for "Equilibrium Near-Bed Concentration of suspended sediment", J. of Hydra. Eng. ASCE Vol 127 No 5, p430-423
    [117] Yafei Jia, Sam S. Y. Wang. Two-dimensional Hydrodynamic and sediment Transport Model For Unsteady Open Channel Flows Over Loose Bed, NCCHE, the University of Mississippi University, MS38677, 2001
    [118] 朱庆平,余欣.关于黄河数学模型系统建设的思考,人民黄河,2005(3).
    [119] 朱庆平,芮孝芳等.黄河下游基于GIS的二维水沙数学模型水沙构件设计,人民黄河,2005(3).
    [120] 朱庆平,芮孝芳等.黄河下游基于GIS的二维水沙数学模型测试与率定,人民黄河,2005(3).

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700