用户名: 密码: 验证码:
基于GIS的河流动力学模型及风险图制作应用研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着计算机技术的迅速发展及人们对流体运动基本理论的进一步认识,河流动力学数学模型在科学研究、工程设计、环境和生态影响评估等相关领域得到广泛的应用。
     本文在传统的河流动力学数学模型的基础上,构造适应复杂边界条件的方程离散格式和数值方法,考虑了非均匀沙沉速、床沙级配调整及横向冲刷坍塌变形等关键因素,针对科研、生产及学科发展的现实和未来需求,耦合了基于GIS的前后处理系统。
     对于FVM方法从空气动力学移植到浅水流动的模拟后,由于采用“平底”的近似计算在水深梯度大的地方所产生的虚假流动,本文采用三点高斯积分法来计算各边的平均法向通量,避免了由水压力与底坡项不平衡所导致的动量不守恒问题,不但避免了由于大的水深梯度所产生的虚拟流动现象,也真实反应了状态向量在控制体内的分布。
     对于生成计算域地形边界条件所需的重要数据源-地形大断面,本文采用道格拉斯-普克数据抽稀算法提取断面特征节点,根据断面节点的干湿属性识别主槽和嫩滩分界点,根据生产堤和控导工程特征识别出嫩滩和老滩分界点,使得精确快速生成大断面标志性数据成为可能。
     河流动力学数学模型的可视化是以大量的区域空间信息为基础的,GIS的发展为数学模型可视化提供了强大的技术支撑。本文利用ArcGIS软件及GIS的二次开发接口,实现了模型与GIS的高效集成,设计和开发出基于GIS的前后处理系统,在地形高程的生成精度、信息查询速度和可视化程度方面有了明显提高。
     将经过验证的模型系统应用于东平湖蓄滞洪区风险图制作,对整个湖区进行了水动力学数值模拟计算。针对50年一遇、100年一遇、1000年一遇不同频率洪水,根据蓄滞洪区相应的调度预案,仿真模拟了洪水的演进路线、到达时间及淹没水深,预测了东平湖区各处在遭受超标准洪水时的危险程度。
With the development of information technology and fundamental fluid theory, hydrodynamic mathematical model has become even more important in different areas including scientific research, engineering design, environment and ecology impaction assessment.
     Based on traditional hydrodynamic model, the paper constructed the methodology scheme and method to suit for the complicated boundary, concerned about the non-uniform setting velocity, bedload grain size adjustment and lateral deformation. Aims to meet the future requirement of research and practical cases, the established model is combined with pre-processing and post-processing system based on GIS.
     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.
     As for the cross-section measured data treatment which guarantee to the correct boundary condition, Douglas-Peucker method is used for deriving the characteristic points to classify the main channel, floodplain and temporary dike effectively.
     By using ArcGIS and secondary developing with ArcGIS, the model’s pre-processes and post processes which including digital elevation generation system and flow field visualization system ,not only improve the efficiency but also improve the precision of digital elevation generation, information searching speed and visualization technology obviously.
     After the model verification, sets of scenario studies are carried out in the floods risk map of Dongpinghu Lake including different flood probability of occurring once per 50 years, 100 years and 1000 years. The model is applied to justify the flood risk map in terms of the flood proceeding route, flood arrival time, flooding area and flow velocity distribution analysis.
引文
[1]朱庆平等,关于黄河数学模型系统建设的思考.人民黄河.2004,6:3-5.
    [2]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,MS 38677,2001.
    [3]秦崇仁,改进的系列模型延伸法及其在航道淤积研究中的应用.海洋学报,1992,14(6): 113-119.
    [4]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 :367-378.
    [5]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.:36-39.
    [6]Molls T,Chaudhry M H. Depth averaged open channel flow model,Journal of Hydraulic Engineering,ASCE,1995,121(6):453-456.
    [7]Price R K,Ahmad K,Holz P.Hydroinformatics concepts.Hydroinformatics tools, Marsake,J.et al(eds.),Kluwer Academic Publishers,Netherlands,1998.47~76.
    [8]Yang C.T.Sediment Transport-Theory and Practice, The McGraw-Hill Companies,Inc.1996
    [9]秦崇仁,肖波,高学平.波浪水流共同作用下人工岛周围局部冲刷的研究.海洋学报,1994,16(3): 5.
    [10]胡春宏,周文浩.中国水科院十年来泥沙研究综述.泥沙研究.1999,(6).
    [11]李义天,谢鉴衡,王明甫.武汉水利电力大学泥沙研究进展.泥沙研究.1999,(6).
    [12]汤立群,唐洪武,陈国祥.河海大学泥沙研究进展与瞻望.泥沙研究.1999,(6).
    [13]朱自强,等著.应用计算流体力学.北京航空航天大学出版社.2001.1.
    [14]韩其为,何明民.泥沙数学模型中冲淤计算的几个问题,水利学报.1988,(5)
    [15]韩其为,何明民.水库淤积与河床演变的数学模型.泥沙研究,1987.(3)
    [16]韩其为,何明民,陈永灿.恢复饱和系数研究及挟沙能力的初步分析,中国水利水电科学研究院.1995.
    [17]韩其为,何明民.恢复饱和系数初步研究,泥沙研究.1997,(3):12-15.
    [18]韩其为.水库不平衡输沙的初步研究,水库泥沙报告汇编,黄河泥沙研究协调小组编印.1973.
    [19]韦直林.天然河道平面二维水沙数学模型研究.武汉水利电力大学学报.1992,(6).
    [20]张世奇.二维动边界潮汐输沙及河床变形计算.泥沙研究.1988,(4).
    [21]刘树坤,宋玉山,程晓陶,等.黄河滩区及分滞洪区风险分析和减灾对策.黄河水利出版社,1999.
    [22]张红武,姚文艺,张俊华,钱意颖.黄科院 90 年代泥沙研究进展.泥沙研究.1999,(6).
    [23]余欣,杨明等.基于 MPI 的黄河数学模型并行计算技术研究及初步实现.人民黄河.2005(3) .
    [24]朱庆平,芮孝芳等.黄河下游基于 GIS 的二维水沙数学模型水沙构件设计.人民黄河.2005(3).
    [25]白玉川,王尚毅.流速和潮位变化作用下波浪在近岸区的传播.《海洋学报》.1996,18,(3) .
    [26]张俊华,张红武,李远发等.水库泥沙模型异重流运动相似条件的研究.应用基础与工程科学学报.1997,(3) .
    [27]姚文艺,彭乃志,杨道富.水土保持治理对河道系统影响研究的若干基础理论问题.泥沙研究.1999,(4):20-25.
    [28]Yang Ming,Zhang Junhua,Zhang Hongwu.Numerical Modeling of Dongzhuang Reservoir on Jing River.Journal of Hydrodynamics, Ser.B Vol.13 No.2 2001.
    [29]Ikeda,Sawai K.Bend theory of river meanders,Part Ⅰ :linear development.Journal of Fluid Mechanics,112:363-377.
    [30]Hasegawa K.Universal bank erosion coefficient for meandering rivers.Journal of Hydraulic Engineering,ASCE,1989,115(6):744-765.
    [31]夏军强.河岸冲刷机理研究及数值模拟.清华大学水利水电系博士论文.2002.4.
    [32]Engelund F.Flow and bed topography in channel bends.Journal of Hydraulic Div.ASCE,100,1974,(11):1631-1648.
    [33]夏军强,王光谦.考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟.水利学报.2002, (6).
    [34]张红武,张清.黄河水流挟沙能力的计算公式.人民黄河.1992,(11).
    [35]王士强.黄河泥沙冲淤数学模型研究.水科学进展.1996(9).
    [36]韦直林.河道水流泥沙问题的一种有限元解法..武汉水利电力学院学报 .1990 年 6 月.
    [37]王新宏.复式断面明渠平面准二维流场的数值模拟.中国水利学会水力学专业委员会2002 年年会.2002,(4) .
    [38]曹如轩.高含沙水流挟沙力的初步研究.水利水电技术,1979,(5): 55-61.
    [39]钱意颖,曲少军.黄河泥沙冲淤数学模型 .郑州:黄河水利出版社.1997.
    [40]韩其为,何明民.论非均匀悬移质二维不平衡输沙方程及其边界条件.水利学报.1997(1)
    [41]李义天,等著.河道平面二维水沙数学模型.中国水利水电出版社.2001.
    [42]余明辉,杨国录,等.非均匀沙水流挟沙力公式的初步研究.泥沙研究. P25~28,2001.6.
    [43]张红武,吕昕.弯道水力学.水利电力出版社,1993.
    [44]杨国录著.河流数学模型.武汉水利电力学院出版社.1992.
    [45]钱意颖,等.编著.黄河泥沙冲淤数学模型.黄河水利出版社.1998.8.
    [46]韩其为,何明民.底层泥沙交换和状态概率及推悬比研究.水利学报.1999,(10).
    [47]Odgaard.Meandering flow model development.Journal of Hydraulic Engineering,ASCE,1986,112(12):1117-1136..
    [48]Shimizu Y,Yamaguchi H,Itakura T.Three dimensional computation of flow and bed deformation.Journal of Hydraulic Engineering,ASCE,1990,116(9):1090-1108.
    [49]Велцканов.М.А..Дцнамика Русловых Потоков, Гос.Изд.Техн-тоор.Лит.1955.
    [50]Ming-Hseng Tseng.Explicit Finite Volume Nin-Oscillatory Schemes For 2d Transient. Free-Surface Flows. International Journal for Numerical Method Fluids. VOL.30. p831~843.1999.
    [51]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.
    [52]Mosselman ,E. Disscussion of “river width adjustment I: processes and mechanisms”, Journal of Hydraulic Engineering,Vol.126 No.2,2000 .
    [53]Sangmook Shin. Reynolds-Averaged Navier-Stokes Computation of Tip Clearance Flow In a Compressor Cascade Using an Unstructured Grid.Blacksburg, Virginia.2001.
    [54]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.
    [55]Engelund, F.and E, Hansen. Amonograph on Sediment Transport in Alluvial Streams. Teknisk Forlag, Copenhagen, 1972, p.62..
    [56]Osman A M and Thorne C R. Riverbank stability analysis I: Theory. ASCE, Journal of Hydraulic Engineering.1998.114 (2).
    [57]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.
    [58]谭维炎著.计算浅水动力学-有限体积法的应用.清华大学出版社.1998.10
    [59]Dartzi Pan,Jen-Chien Cheng . A Second-Order Upwind Finite-Volume Method For The Euler Solution On Unstructured Triangular Meshes. International Journal for Numerical Method in Fluids,Vol.16,1079~1098,1993.
    [60]De-yu Zhong,Hong-Wu Zhang. A Two Dimensional Numerical Model for Simulation ofFlood Propagation in Initially Dry Bed. Flood Defence’2002 Wu et al (eds)? 2002 Science Press,New York Ltd.,ISBN 1-880132-54-0.
    [61]Tan,W.Y.Shallow water hydrodynamics.Elservier, Amsterdam, the Netherlands.1992.
    [62]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.
    [63]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) .
    [64]E. Erwin Holly, Numerical Calculation of Flow Field, World Scientific, 1984.
    [65]Biker E. W. Littoral Drift as Function of Waves and Current[C].Proc. 11th Conf. Coastal Engineering,1968.
    [66]I.Demirdzic .A Collocated Finite Volume Method For Predicting Flows At All Speeds. International Journal for Numerical Method Fluids.VOL.16.1029~1050.1993.
    [67]Osman, A. M., and Thorne, C. R. Riverbank Stability Analysis Ⅰ: Theory & Ⅱ Application.ASCE,Journal of Hydraulic Engineering,1988,114(2):134-172.
    [68]Stanley Osher and Fred Solomon.Upwing Difference Schemes for Hyperbolic Systems Of Conservation Laws. Mathematics of Computation.Vol38.NO.158.339~374.1982.
    [69]Francisco Alcrudo .Flux Difference Splitting For 1d Open Channel Flow Equations. International Journal for Numerical Method Fluids.VOL.14.1009~1018.1992.
    [70]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.
    [71]Thorne C R and Osman A M. Riverbank stability analysis II: Application. ASCE, Journal of Hydraulic Engineering.1998.114 (2) .
    [72]Zhao, D. H., Shen, H. W., Lai, J. S., and Tabios, G. Q. (1996). “Approximate Riemann solvers in FVM for 2D hydraulic shock wave modeling.” J. Hydr. Engrg., ASCE, 122, 692-702.
    [73]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..
    [74]Pan H,Damodaran M.Parallel computation of viscous incompressible flows using Godunov-projection method on overlapping grids.Int J Numer Meth Fl,2002,39:441-463.
    [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]Slagter W J,Florie C J L,Venis A C J.Advances in 3D forging process modeling.Journal ofComputational Physics,1991,92:82~105.
    [77]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.
    [78]潘承洞,潘承彪.代数数论[M].济南:山东大学出版社,2001.
    [79]Pan H,Damodaran M.Computation of unsteady viscous incompressible flows in...FVM algorithm for three-dimensional viscoelastic flows.Computational Mechanics,2003,30:265-280
    [80]陈新一.一种多点迭代方法.甘肃教育学院学报(自然科学版),2001,15(1):13-16.
    [81]李国英.治理黄河思辨与践行.黄河水利出版社.2003.8.
    [82]韩其为,何明民.水库淤积与河床演变的数学模型,泥沙研究.1987 年 3 期.
    [83]谢鉴衡,魏良琰.河流泥沙数学模型的回顾与展望.泥沙研究.1987
    [84]张瑞瑾,编,河流动力学.中国工业出版社,1961.
    [85]涂启华,安催花,李世滢,余欣等,著.《水利动能手册(泥沙分册)》,中国水利水电出版社,2004.
    [86]曹如轩.高含沙水流挟沙能力初步研究.水利水电技术.1979,(5).
    [87]谢鉴衡,河流模拟.水利电力出版社,1990.
    [88]刘月兰, 韩少发, 吴知.黄河下游河道冲淤计算方法.泥沙研究, 1987,(3) .
    [89]谢鉴衡,魏良琰.河流泥沙数学模型回顾与展望.泥沙研究.1987(3).
    [90]刘善建,麦乔威.黄河下游治理开发中遇到的泥沙问题.泥沙研究.1958(4).
    [91]Van Rijin L C.Sediment Transport. Part I: Bedload Transport.ASCE.Journal of Hydraulic Engineering, 1984.110(10) .
    [92]舒安平.水流挟沙能力公式的检验与评述.人民黄河.1993,(1).
    [93]韦直林等.黄河下游变动河床洪水演进与洪水位预报水动力学数学模型研究.武汉水利电力学院河流数值模拟研究室.1990,3.
    [94]王光谦,张仁,惠遇甲,王桂仙,王兴奎.清华大学泥沙研究回顾.泥沙研究.1999,(6).
    [95]万洪涛,周虎成,吴应湘,万庆.黄河下游花园口~夹河滩河段二维洪水模拟.水科学进展. Vol.13,No.2,P215~222,2002.3.
    [96]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).
    [97]Van-Rijin L C and Wijinberg K M. One-dimensional modeling of individual waves and wave-induced longshore currents in the surf zone. Coastal Engineering, 1996.28:121-145.
    [98]Van-Rijin L C. Sedimentation of dredged channels by currents and waves. Journal of Waterway, Port, Coastal and Ocean Engineering, 1986, 112(5):541-559.
    [99]朱庆平.数字黄河工程规划项目概述.人民黄河.2003.8.
    [100]安新代,王道席,王祥辉.数字黄河工程规划项目概述.人民黄河.2003.8.
    
    [101]李琦 , 杨超伟 , 陈爱军 .WebGIS 中的地理关系数据模型研究. 中国图像图形学报,2000,5A(2):119~123.
    [102]李青元,刘晓东,营代勇.WebGIS 矢量空间数据压缩方法探讨.中国图像图形学报,2001,6A(12):1 125~1129.
    [103]唐泽圣,徐志强.二维点集三角剖分的动态生成与修改.计算机辅助设计与图形学报,1990,3:1-8.
    [104]王净,江刚武.无拓扑矢量数据快速压缩算法的研究与实现.测绘学报.2003,32(2):173~177.
    [105]马正飞 , 殷翔著. 数学计算方法与软件的工程应用 .p111 ~ p115. 化学工业出版社.2002.12.
    [106]D.J.Wheeler Commentary:Lin king Environmental Models with GIS for G lobal Change Research[J].PE&RS.1993,59(1 0):1497~1501.
    [107]王尚毅,顾元蒳,郭传镇.河口工程泥沙数学模型.海洋出版社.1990.120~192.
    [108]杨明,秦崇仁,刘振纹等.基于 Pelnard-Considere 方法的岸线预报模型.第六届全国泥沙基本理论会议论文集,2005.
    [109]李大鸣.多重尺度水沙数学模型的理论及应用.天津大学博士论文. 1996.7.
    [110]员瑛,秦崇仁.潮汐河口半潮平均流速沿深度的分布及悬移质输沙量的计算.海洋科学.1980(04),P10-12.
    [111]江春波,邢秀英,张庆海.用分步有限元法求解三维不可压缩流动. 清华大学学报,2000,Vol.40(8):110-113.
    [112]李国英.黄河中下游水沙的时空调度理论与实践.水利学报,2004.8
    [113]李国英.黄河第三次调水调沙试验的设计、实施、效果与认识.黄河报.2004(36) .
    [114]水利部.洪水风险图编制导则.2005.3.
    [115]美国国家科学基金委员会.防洪减灾总报告及研究规划[M].谭徐明等编译.科学出版社.1997.
    [116]刘树坤,杜一等.全民防洪减灾手册[M] .辽宁人民出版社.1993 年.
    [117]徐乾清.对未来防洪减灾形势和对策的一些思考.水科学进展,1999(3):235~241.
    [118]陈仕亮主编.风险管理[M] .1994 年 8 月第 1 版.成都,西南财经大学出版社.1994 年13~47.
    [119]赵世来等.东平湖水库运用与移民开发[M],郑州:黄河水利出版社,2000.187.
    [120]张锡令,等.南水北调东线可行性研究阶段东平湖水库工程地质勘察报告[R].天津:水利电力部天津勘测设计院.1988.
    [121]胡一三.黄河防洪.黄河水利出版社.1996 年.
    [122]东平湖志.山东大学出版社.1993 年.
    [123]山东省防汛抗旱应急预案.山东省政府.2005.12.
    [124]东平湖防汛资料手册(内部文件).黄河防汛办公室.2002.6.
    [125]张旭.东平湖治理和防洪运用刍议.山东黄河河务局东平湖管理局科研报告.2005.4.

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

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

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