洪涝模拟洪涝模拟

6IInChinawaterloggingischaracterizedbyhighfrequency,longtime,heavylosses,wideinfluentialareaandcomplexinteractionwithflooding,etc.Manystudieshavebeengenerallylimitedtolocalexperiencesandtativeysisforalongtime.Inthispaper,appliedflownumericalsimulation,combinedwithdrainagetheoryandhydrography,awaterloggingandfloodingnumericalmodelisestablished,whichcancalculatethedepth,areaandratioofwaterlogging.Theinteractionbetweenwaterloggingandfloodisyzedtatively.Themodelresultscanprovidetechnologicalsupportforthewaterloggingforecasting，assessmentandfloodandwaterloggingcontrolplan.Thisthesisincludesixchapters,andmaincontentofeachisasItsumsupthecharacteristicofwaterloggingandthecurrentresearch,andpointouttheshortagesoftheexitingwaterloggingnumericalmodelsinchapterone.Inchapter2,flownumericaltechnologyissummarizedly；basedonthefinitevolumemethodandunstructuredmesh,thetwo-dimensionalwaterloggingnumericalisestablishedsecondly；Forthelast,linkingaone-dimensionalflowmodeltothetwo-dimensionalmodel,acompositivewaterloggingandfloodingmodelisestablished.Lowlanddrainageissimulated;thedepth,areaandratioofwaterloggingaregained,thenegativeinfluenceoffloodonthewaterloggingisstudied,inchapter3.Inchapter4,theseriesofrainfallintheSishuibasininHunanprovinceiscalculatedbyhydrologicalmethod;theaveragedrainingmethodisappliedtocalculatingthedrainingInchapter5,applyingtwo-dimensionalmodellinkedone-dimensionalnumericalmodel,thedrainageandfloodingoftheXiangriver,theZIriver,theYuanriver,theLiriverandtheDongtinglakearesimulated,thedischargeoftheoutletindrainagebasin,thewaterlevelofSishuiriver,andtheflowfieldofDongtingLakearegained,andtheinfluenceofdrainageonfloodisyzed.chapter6sumsupthewholetheris,yzssomeproblemsofthemodelsandgivessomeconstructivesuggestionsforfuturestudy.:numericalmodelingofdrainage;two-dimensionalflowmodel;ratioofwaterlogging;SishuiriverandDongtinglake;two-dimensionallinkedwithone-dimensional;theinteractionbetweenwaterloggingandflooding.1章绪论涝灾研究综述前全球各种造成的损失，洪涝占40%，热带气旋占20%，干旱占15%，占15%，其余占10%[1]。世界上洪涝最为频繁而又严重的国家之一，洪涝灾害对社会经济发展影响巨大。半个多世纪以来，中国平均每年因洪涝受灾耕地流域圩区，珠江流域的珠角洲等地区[2]。洪涝多年平均损失在800～550～750亿元。可见与洪灾相比，涝灾造成的损失往往更为严重。农田水利学及学等多个学科。长期以来各地针对当地的涝灾特点开展了大量研1.1表 我国各涝渍易损地区涝灾防治研[2]、鄱阳湖湖步的成果，各数学模型概况见表1.2。1.2模型等4个基础分析模型，GIS1998洪水部分与平原湖区相比，城市经济较发达，对涝灾的观测与较齐全，在泵、闸等，其中排水是非常重要的组成部分；在圩垸排涝中，主要涉及农业的堤垸区，垸内地势低洼，一般低于大堤3～5m，部分地段低于外湖（河）水面达5～8m。外江外湖水易通过潜、渗流入垸田，抬高垸内水位，酿成涝渍[10]。的水文效应[32]。以长江干流为例，城陵矶河段，每增加流量6000m3/s，高水位时抬1m3/s1.3m。1991年长江洪水不大，沿岸内涝严重，大量电排提水入江，在高洪时期，造成很不利的影响，江汉平原以排水流量4700m3/s计算，可抬高关水位0.6m左右[33]。另外，根据1990年以前的1931年、1935年、1954年等典型年洪水进行计算分析，由于排涝水量的长江流域易涝区分布面广，全流域内易涝耕地约7000万亩，其下游平原区占93％。区内地面高程45～3m，普遍低于当地洪水位数米至十余米，涝灾既广泛又频原湖区涝渍的影响，从上个世纪八十年始已有较多研究成果见诸于期[36~38]。但到目前为止，就这一问题还时有不同看法。大多研究者认为，三峡工平原湖区提供可靠电源，有效保证电排设施的正常使用[37~42]。另有一些研究者则认排水量减少，加重排涝负担[43~45]。鉴于目前在这一问题上的，就三峡工程防洪与排涝相互影响并在一定程度、一定意义上相互。协调好防洪与排涝的关系，关系到人民生命安全和社会经济发展。利用资料分析和数学模型模的分析与研究，以为下一步的涝灾治理提供理论依据。本研究方法及内水利学上的平均排除法对199871日~912日湖南省四水流域降雨后的排涝199871日~1998912日各流域净雨过程、各流域汇流出口流量199871199871日~1998912日各流2章洪涝数学模型水流数值模拟20世纪七、八十年代伴随着电子计算机发展起来的，一维数学模型比较，须采用数学模型来研究解决上述问题。二维数学模型[47]于河口、海岸水流泥沙运动的研究，HansonＷ.是早期进行20502080、90年代计算机软硬件水平的飞速发 特征[51]：在计算机普遍应用之前河流模拟的数值计算主要是利用特征（MethodofCharacteristics，MOC）理论采用图解法等进行手工计算。因特征求有限差分方法（FiniteDifferenceMethod，FDM）[52]：是计算机数值模拟最早采代替连续的求解域。FDM把控制方程中的微差商代替进行离散，从而建立代数FDM一般只适用于矩形网格，用于天然河道计算时一般需采用坐标变换。矩形网格计算格式复杂、计算量及量较大，大型系数矩阵较难求解等缺点[57-63]。有限体积法（FiniteVolumeMethod，FVM）[64]：又称有限控制容积积分法，该整个计算区域，都严格满足物理守恒定律，不存在守恒量的误差[65-66]。由于有限体洞庭湖洪涝模型图 洪、涝物理模型2.1。洪是河道中峰高量大涝关系见表2.1。表 洪与涝关系研究洪线涝面Delaunay三角形化方法生成较优非结构网格；采取有限体积法对平面二维水流QA Q(Q2)gAZ

(2- (2-（m3/s；x（m；A（m2；t/s；（m3/s2；z为水位（m；n为糙率系数；R为断面水力半径涝区一般边界复杂，本文主要对Delaunay三角形化方法进行研究，目的是为数2080-90年代得到了迅速的发展，目前己成为生成非结构化网格的一种主要方法[75-78]。根据Delaunay三角形化方法，本文按照以下步骤生R(k)最大的在序列的最顶端。（4）R(k)max(R(k))>为设定Delaunay三角形化方法连接新点与有关的网格点，生成一组新的DelaunayDelaunayR(k)并按大小排序，R(k)最大的在（5）max(R(k))<，结束，否则返回到（4在网格生成过程中，对网格生成质量和速度有较大影响的几个问题是：（1）判断（2项”来控制网格的疏密；（3）通过Laplace光顺化方法移动网格点，使它位于包ZMN gn2Mu2h4 gn2Mu2h4txyfHvghxtx2y2 2 2N

gn2Nu2v

h fHuh

t

x y2

其中：H为水深（m）uvxy方向的流速（m/s，M=uh，N=vh；Z为水位（m；nManning糙率系数；t为紊动粘性系数；q为单位面积上水流的源汇强度（ms。图2.2。图 节点布置3 (φφutφ)dA 3

j1 Vpj为控制体中界面（边）uhvh。由式（4）可知，为了由此获得关于节点P0处变量的代数方程，关键是如何将界面上的对流项与扩散项的积分结果表示成P0及与P0相邻节点上之值，这也是不同作者非结构化网非稳态项离散φd uh dA j

P0

P0P0处变量。对流项的离散j上的对流项为CjCjφudA(uhu)jAjFjuuj

K

{0.5*[(gz)P

ddjd

ddj]gzPjzP0}dd )K0.5*[(V V)])( auP au

a0a由式（2-9）可以看出界面流速公式由三部分组成：线性化插值部分uPoj、在P0Pj单元的平均水位梯度作用下所引起的速度变化及假定水位是呈线性变化时，P0和Pj间平均水位梯度作用下的速度变化。FuCDSF(uCDSuUCD j (2-扩散项的离散DjtuhdA（tuh）j nnP0、Pj连线方向的分量jh )Ah ) uPjuP0dh )Ah )SNddDSNdd

t j j j

j

2.3水位梯度项的离散根据Gauss定理，可得水位梯度项ghzi

3ghz j 33GghzidAghz j动量方程的建立Na0uP0ajuPj

ta (dA)min(Fd d

u ub F(uCDSuUCD)h( j)AghzASSNSNj

t

c

a

hV3P3P a j

PP

欠松弛处理16，速度的欠松弛因子直接组织到代数方程组的求解过程中，最终生a0uP

NNaubj

1

.au0

压力修正方程的离散SIMPLEddjz'z' u'KjghPj P0 u*u度P0

u j应能使（j

uapz apz'b0 j z ajk

g

hs

bpP d

z'uP0u*gh P z*

P zP P0P 非结构化同位网格SIMPLE算法的计算步骤①给定初始速度场和水位（压力）场；②由给定的速度场和水位场计算离散动量方程组的系数和③求解动量方程组，得出速度场u*和v*④计算水位校正方程中的系数和源项；'⑥按照（3.32）修正水位和速度，获得本层次上满足连续方程的速度和水位场；⑦反复迭代，直式是最好的选择[92-93]，见图2.4、图2.5。2.4二维模 一维模2.5小3199897871345亿元，其中300194亿元。1950~199546合下的涝灾情况，为涝灾和评估提供技术支持。大通湖垸排涝图 大通湖垸示意面积169.03万亩，合1127.43km2。模型z(hu)(hv) (3- uuuvugzgn2uu2v

2u (3-

t

h4

x

y2vuvvvgzgn2vu2v

2v (3-

t

h4

x

y2g为重力加速度;ν为水流紊动粘性系数；n为糙率系数[6]。定解条件包括边界条件和洞庭湖，内有大通湖（见图3.1。考虑到内湖湖水为渍堤所隔开，堤内外流动特性不动边界处理：在计算过程中，计算域分节点有时会被“淹没”，有时会“干hdrying时水域处理，不参与计算；水深大于某值hflooding时作为水域处理，参与计开边界，主要开边界为：（1）明山大电排（抽排流量为126m3/s，大东口电排（抽排90m3/s中；（2）五门闸（底板高程24.5m，自流排水入东洞庭湖。大通湖水系主要外排泵站概化情况见表3.1。3.1沥涝，导致农田低产。旱作物的排涝历时一般为2天，水稻的排涝历时一般为3天。6、7月份为洞庭湖地区涝灾多发期，这个时期属于农事活动繁忙时期，5月底至6月初抢收冬小麦，该时期发生涝灾将会影响小麦收割和其它旱地作物播种；7月棉模型大通湖垸排3210mm3日排出，同时，大通湖垸最3日降雨3308.7mm。本文分别对这两种雨型所形成的沥涝情况进行模拟计算。按不同外主要外排泵站按照设计流量抽排；（2）当外河水位较高时，扬程增大，参考大通湖垸4599个三角形网格，三角形最大40m2，最小允许角度为29。,3.2。3.20.02~0.05；动边界系数：hdrying=0.005m，hflooding=0.010m；计算时间步长：60s。紊动粘性系数：采用Smagorinsky公式确定：uuxyx v v22 ytCS据经验取0.28。3.2耐淹历时淹没水深成熟（6月底~7月初3开花、结铃（7月1拔节抽穗（7月1苗期（6月下旬~7月上旬5开花期（7月中下旬2全期（5月底~8月初23210mmT=1440△t1（10cm弓、千山红和南大膳的的局部地区积水甚至超过0.5m。T=2880△t23.4可以看出，3210mm20.1m的区域结铃，造成棉花减产；超过玉米的耐淹历时（1天，影响玉米拔节抽穗；超过大豆耐淹历时（2天，影响大豆开花；超过甘薯耐淹历时（2天，影响甘薯生长。积水0.5m的区域主要分布在八百弓、千山红、南大膳、明山头的局部地区，如图 3日210mm降雨及设计抽排流量情况下第1日末致涝分布3.43日210mm降雨及设计抽排流量情况下第2日末致涝分布图主要泵站按照设计流量抽排在计算时段T=4320△t即第3过0.5m的区域主要分布在八百弓、千山红、南大膳、明山头的局部地区，这些地方可知，在南大膳地区水稻致涝面积逐日增大，3日末已扩大到3210mm雨型下各天的致涝率（33可以水深，致涝率为24.16%。因此，内湖满蓄之后，大通湖垸设计排涝流量不足，已有排涝能力造成24.16%的致涝率，不能满足排涝要求。效率下降。将3日210mm降雨第3日述3种外洪水位遭遇情况下的沥涝情况进3.43.4权值分别为：24.16%、24.38%25.47%。当外洪水位抬高使大通湖垸外排泵站效率下降20%时，垸内致涝率加重0.22%；当明山外洪水位达驼峰水位时，垸内致涝率加重1.31%。可见外洪水位抬高对内涝造成明显的不利影响。表 3日210mm降雨在设计抽排流量下的致涝率逐日计算深表 3日210mm降雨在不同外洪条件下致涝率计算率值20%3308.7mm3210mm雨型。外排泵站按照设计抽排流量抽排时第3日末的水深分布图见图3.6。千山红农场、南大膳、金盆镇旱作物致涝面积、水稻致涝面积均有所增加。3日308.7mm3.53.5可见，在三种不同外洪水位工况下致涝率值分别为：30.77%、31.05%和31.88%。当外洪水位抬高使大通湖垸外排泵1.11%3210mm降雨情况排时，致涝率由24.16%增加到30.77%，上升了6.61%；当外水位较高外排泵站效率下降20%时，致涝率由24.38%增加到31.05%，上升了6.67%；当外水位高至明山驼峰水位明山被25.47%增加31.88%6.41%3308.7mm降雨与210mm降雨相比，致涝率平均增加约6.56%。3.5.3210mm3图 3日308.7mm降雨及设计抽排流量情况下第3日末致涝分布表 3日308.7mm降雨在不同外洪条件下致涝率计算值小6种计算工况下的沥涝情况。模型计算结果包括3日内任意时段的水深分布。3210mm3种不同外洪条件下3308.7mm3210mm两种雨型下3308.7mm3210mm降雨致涝率平均增加约模型计算结果及所采取的处理方法可为涝灾与评估提供技术支持。4章四水流域排涝计算20世纪七十年代以来，湖南省大搞电排建设，兴建了大量电排站，使得垸内涝涝区面雨量 进行面雨量计算。1998年7~9月四水流域面雨量见表4.2。单位：717273747576777879710711712713714715716717718719720721722723724725726727728729730731818283848586878889810811812813814815816817818819820821822823824825826827828829830831919293949596979899910911912净产水量分别进行扣损计算。见图4.3。 4.3R1等于时段净雨量P减去时段水面蒸发量E，然后可由蒸发皿的蒸发量折算确定。R1=P-S湖面-E湖 净产水深R2；R2=P-S水田-E水田 R3=P-I0- 式中I0为初损，且I0=Im-Pa，Im为最大初损量，Pa为前期影响雨量，均可通过降

水量计算见表4.3~4.6。排涝入江流量涝水位。蓄涝电排工作方式可考虑两种：（1）汛期采取随时降雨随时电排的方式，即水滞蓄于垸内，随后排出；（2）当圩垸或涝区蓄水（如内湖）到电排水位时，才开机m3/s、800m3/s、600m3/s、355m3/s。湘湘 湖 水 总 耕地量3(m1998-7-001998-7-0001998-7-001998-7-0001998-7-1998-7-001998-7-001998-7-01998-7-001998-7-001998-7-001998-7-001998-7-411998-7-01998-7-001998-7-001998-7-001998-7-001998-7-001998-7-001998-7-001998-7-001998-7-001998-7-01998-7-001998-7-001998-7-01998-7-001998-7-001998-7-001998-7-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-01998-8-001998-8-01998-8-001998-8-01998-8-00表 湘水流域排涝流量计算时 000000000000000000001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-001998-8-1998-8-01998-8-001998-8-001998-8-001998-9-001998-9-0001998-9-0001998-9-1998-9-1998-9-001998-9-001998-9-001998-9-001998-9-001998-9-001998-9-04.400000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004.5总总净加水量3流旱时 沅湖水及耕4.6总总净加水量3流旱时 澧湖水及耕天然径流与排涝后径流01998-6-

1998-7-

1998-7-

1998-7-

1998-8-

1998-8-

1998-8-

1998-9-

1998-9-

1998-9-图 湘江株洲天然流量、排涝后流量对比

01998-24

1998-7-

1998-14

1998-24

1998-8-

1998-13

1998-23

1998-9-

1998-12

1998-224.5实测实测沅水径流排涝后沅水径流01998-7- 1998-7- 1998-7- 1998-7- 1998-8- 1998-8- 1998-8- 9-94.6实测实测澧水径排涝后澧水径流01998-6-241998-7-41998-7-141998-7-241998-8-31998-8-131998-8-231998-9-21998-9-121998-9-4.7峰流量明显增大，其中，199872448897m3/s小5模型计算湘、资、沅、澧四水流域总面积262823km2。其中湘江流域株洲水文站以上71919km2，干流长670km，距河口约130km26704km2，距河口44km；沅水流域桃源水文站以上集水面积85223km2距河口68km。四水流域计算范围见图5.1及表5.1。5.15.1距河口长度5.2—— 5.25.3 图 洞庭湖一二维耦合模型平面模型199871日~年月日实测洪水资料、年实测地形，进行5.8水位水位

14

1998-24

13

1998-225.91998-6-141998-7-41998-7-241998-8-131998-9-21998-9-5.10 1998-6-141998-7-41998-7-241998-8-131998-9-21998-9-5.115.9~5.11可知，四水一维河道洪水模拟结果与实测结果拟合较好，说明入湖边界。计算结果见图5.12~5.15。0

1998-7-241998-8- 5.1201998-6-141998-7-41998-7-241998-8-131998-9-21998-9-5.13沅江沅江河口流量过01998-6-

1998-7-

1998-7-

1998-8-

1998-9-

1998-9-5.1401998-6-141998-7-41998-7-241998-8-131998-9-21998-9- 5.155.16。5.16参系数取值：糙率采取二维糙率，深槽取0.020~0.028，浅滩取0.025~0.033；动边uuxyx v v22 ytCS式中，u、vx、y方向垂线平均流速，Δ为网格间距，Cs为计算常数，根据经验取0.28。题。图5.17~5.32为计算中某些时刻的水位分布图和流场图。Surfaceelevation[m]Above32.165-32.163-32.161-32.16-32.159-32.157-32.155-32.154-32.153-32.151-32.15-32.148-32.146-32.145-Below32.145UndefinedValue0:00:001998-07-08TimeStep336of5.17t=7dSurfaceelevation[m]Above31.596-31.594-31.593-31.591-31.59-31.588-31.587-31.586-31.584-31.583-31.581-31.58-31.578-31.576-Below31.576UndefinedValue0:00:001998-07-15TimeStep672ofSurfaceelevation[m]Above31.848-31.844-31.84-31.836-31.832-31.828-31.824-31.82-31.816-31.812-31.808-31.804-31.8-31.796 Below31.796UndefinedValue0:00:001998-07-22TimeStep1008ofSurfaceelevation[m]Above33.206-33.204-33.202-33.2-33.198 33.196-33.194-33.192-33.19-33.188-33.186-33.184-33.182-33.18- UndefinedValue0:00:001998-07-29TimeStep1344ofSurfaceelevation[m] 33.237-33.235-33.232-33.23-33.228-33.225-33.222-33.22-33.217-33.215-33.213-33.21-33.208-33.205-Below33.205UndefinedValue0:00:001998-08-05TimeStep1680ofSurfaceelevation[m]Above33.297-33.296-33.294-33.292-33.291-33.29-33.288-33.286-33.285-33.284-33.282-33.28-33.279-33.278-BelowSurfaceelevation[m]Above33.702-33.7-33.698 33.696-33.694-33.692-33.69-33.688-33.686-33.684-33.682-33.68-33.678-33.676-Below33.676UndefinedValue0:00:001998-08-19TimeStep2352ofSurfaceelevation[m]Above33.324-33.322-33.321-33.319-33.318-33.317-33.315-33.313-33.312-33.311-33.309-33.307-33.306-33.305-BelowCurrentspeed80.060.040.020-0.020Below-0.02UndefinedValue0:00:001998-07-08TimeStep336of5.25T=7dCurrentspeed[m/s]Above0.3-0.275 0.25-0.225-0.2-0.175 0.15-0.125-0.1-0.075 0.05-0.025-0--0.025 Below-Currentspeed0.40.360.320.260.120.080.040-0.040-0.08---0.12--Below-0.12UndefinedValue0:00:001998-07-22TimeStep1008ofCurrentspeed0.360.320.260.120.080.040-0.040-0.08---0.12---0.16--Below-Currentspeed0.440.40.360.320.260.120.080.040-0.040-0.08--Below-0.08UndefinedValue0:00:001998-08-05TimeStep1680ofCurrentspeed[m/s]Above0.325-0.3-0.275 0.25-0.225-0.2-0.175 0.15-0.125-0.1-0.075 0.05-0.025-0- Undefined0:00:00199808-12TimeStep2016ofCurrentspeed0.40.360.320.260.120.080.040-0.040-0.08---0.12--Below-0.12UndefinedValue0:00:001998-08-19TimeStep2352ofCurrentspeed[m/s]Above0.3-0.275 0.25-0.225-0.2-0.175 0.15-0.125-0.1-0.075 0.05-0.025-0--0.025 Below-0.025UndefinedValue0:00:001998-08-26TimeStep2688of排涝后洪水199871日~199812日各计算时刻下排涝之后各流域净雨、各流域汇流出口流量过程、四水河道水基于相同的网格，对水深场作差可以得出任一计算时刻任一计算网格内的水位5.33~.36。Surfaceelevation0.00480.00440.0040.00360.00320.00280.00240.0020.00160.00120.00080.00040-0.00040Undefined0:00:001998-07-08TimeStep336of5.337d四水流域从年月日起持续排涝至年8日，湖面水位平均可抬高Surfaceelevation[m]Above0.0078-0.0072-0.0066-0.006-0.0054-0.0048-0.0042-0.0036-0.003-0.0024-0.0018-0.0012-0.0006-0- Undefined0:00:001998-07-15TimeStep672of5.3414d14d 年月日起持续排涝至 年月15日，湖面水位平均可抬高Surfaceelevation0.0880.080.0720.0640.0560.0480.040.0320.0240.0160.0080-0.0080-0.016--Below-0.016UndefinedValue0:00:001998-07-22TimeStep1008of5.3514d14d排涝至1998-7-22，湖面水位平均可抬高0.042m。Surfaceelevation0.0120.0110.010.0090.0080.0070.0060.0050.0040.0030.0020.0010-0.0010Below-0.001UndefinedValue0:00:001998-07-29TimeStep1344of5.3614d14d5.31998-7-51998-8-131d表 排涝后对洞庭湖湖面平均水位抬高统计5.3199875日~1998724日水位差值统计可知，排涝期间洞庭湖湖面水位抬高平均抬高0.0380m。四水洞庭湖洪涝数学模洪水演算，模拟了四水流域、洞庭湖区间的洪涝形成过程，模型计算结果包括199871日~1998912日各流域净雨过程、各流域汇流出口流量过程、四水河道水位、流速、洞庭湖流场、水深等，见图5.37和图5.38。Surfaceelevation[m]Above34.3234.24-34.16-34.08-34-33.92 33.84-33.76-33.68-33.6-33.52-33.44-33.36-33.28-33.2-Below33.2UndefinedValue0:00:001998-07-27TimeStep144ofCurrentspeed2.221.810.80-0.20-0.4--Below-0.4UndefinedValue0:00:001998-07-27TimeStep144of图 3d210mm降雨排涝后洞庭湖流场澧四水流域和洞庭湖区综合规划、防洪与调度等服务。洞庭湖排涝标准为3d210mm降雨，利用本模型，模拟湘、资、沅、澧四水流域同时遭遇3d210mm降雨6章结论与展望主要本从涝灾研究现状出发，提出了我国涝灾研究中的不足；然后以河流模拟理了模拟，分析了四水排涝对洞庭湖洪水位的影响。的研究成果包括以下几部分：于平原湖区一般边界复杂，本采用了Delaunay三角形化方法剖分网格，为数值24.16%25.47%，定量研究了外洪水位变化对6章结论与展展 ，刘立成等．洪涝[M]．：气象朱建强.易涝易渍农田排水应用基础研究[M].：科学雷声隆,丘传忻,.排涝工程[M].：科学技术国家防汛抗旱总指挥部,水利部水文水资源.中国水旱[M].：中朱建强，欧光华.四湖流域农田排水有关问题研究[J].刘章勇江汉平原涝渍地域农业生态环境评价与综合治理对策[J朱建强,李靖.江汉平原水综合防治研究[J].长江大学学报，2005，25（3）：8-毛德华,夏军.洞庭湖区洪涝的形成机制分析[J].大学学报，2004，51（2）：199-刘辉,蒋孝平.洞庭湖区洪涝成因与治理对策[J].湖南水利王克林,章春华,易爱军.洞庭湖洪涝形成机理与生态减灾和流域管理对策[J].应用生态刘晓清,毛德华,曾桐秀,洞庭湖区洪涝的预报模式[J].自然学报李景保,杨燕,许树辉.洞庭湖区1991～2000年大洪涝特点与成因分析[J].湖南师范、文，洞庭湖洪涝的地学背景及治理策略[J]，湖南水利陈红根,林等.洪涝对鄱阳湖区血吸虫病的影响及其控制策略研究[J].中国血吸,刘影,.鄱阳湖区洪涝的生态环境因素与生态减灾对策[J].江西师范大学刘群红,叶滢.对鄱阳湖区涝渍问题的[J].江西师范大学学报高介福，朱宏富.鄱阳湖地区的综合开发与治理[M].江西人民.].杨世伦.太湖流域洪涝的激发机制和减灾策略探讨[J].学.1991年洪涝及今后治理措施[J].，崔广臣，杨国顺.三江平原洪涝及治理措施[J].水利水电科技进展,坤，，赵正文.三江平原典型垦区的治涝措施[J].中国农村水利水电陈丙咸等.GIS的流域洪涝数字模拟和灾情损失评估的研究[J].李娜,仇劲卫,,陆吉康.市城区暴雨沥涝仿真模拟系统的研究[J].自然学等.城市暴雨内涝数学模型的研究与应用[J赵思健,,熊利亚.利用空间分析建立简化的城市内涝模型[J].自然学报杨克俭,李清泉,.遥感制图技术在洪涝监估、治理中的应用[J].地图童庆禧1998年洪水监测中的作用[J].费永法.平原坡水区除涝治理工程对流域洪水影响的分析[J].水利部长江水利.长江流域水旱[M].：水利水电闵骞，鄱阳湖高水还湖对水位影响的计算与分析[J]方子云,,袁弘任.论三峡工程对水库上游、江汉平原洪、涝、渍及长江河口生态环境的影响[J].环境科学与技术,1986,34,14~17.镇英明.三峡工程对四湖地区水土资源开发的影响.人民长江,1987,10,22-许正甫.三峡水利枢纽对四湖地区排涝的作用[J].人民长江1986,4,秦文凯,府,,为.三峡建坝前后洞庭湖的淤积[J].学报(自然科学版),1998,38(1):84~87.周国华,唐承丽,朱翔,周海燕.[J].水土保持学报,2002,朱翔.洞庭湖水灾机制和减灾分析[J].自然学报.1999,8(1):张英.江汉平原和洞庭湖区洪涝治理的探讨[J].人民长江范大学学报（自然科学版）.1999,20(4):376~379.,三峡工程与江汉平原农业持续发展[J].生态学杂志1995,16(2廖小红.三峡工程对洞庭湖区滨湖城市排涝的影响[J].杨国录，吴卫民.一维河流数值模拟算法的概述.泥沙研究，1995：34-朱家鲲.计算流体力学.：科学杨国录.河流数学模型，海洋..陆金浦，关治.偏微分方程数值解法.：瓘.计算水力学理论与应用.：河海大学Yanenko.N.N.NumericalMathodsinFluidDynamics.Moscow:MIRPublisher,Leedertes.J.J.AWaterQualitySimulationmodelforWell–MixedEstuariesandCoastalSeas.PrinciplesofComputation.TherandCorporationRM-6230-RC,1970(1).Peaceman.D.W.andRachfordjr.H.H.TheNumercalSolutionofParabolicandEllipticDifferentialEquationgs.J.Soc.Ind.Appl.Math,1955.3.计算水力学.：谭维炎.计算浅水动力学—有限体积法的应用.库克.有限元分析的概念和应用.何穷，程耿东译.：科学.有限元基础.：李开泰等.有限元方法及其应用.西安：西安交通大学姜尚礼，厐之垣 有限元方法及其理论基础，：人民教育章本照.流体力学的有限元方法.：机械工业谭维炎，计算浅水动力学—有限体积法的应用，，黄素逸.计算流体力学基础.：,[美]S.V.帕坦卡著.传热与流体流动的数值计算.张政译.：科学程心一.计算流体动力学.：科学，谭维炎一维不恒定流计算的三种高性能差分格式，谭维炎无结构网格上二维浅水流动的数值模拟..计算水动力学（讲义）.水利电力学院林秉南.明渠不恒定流研究的现状和发展.水利水电科技进展.：水利苑明顺.计算流体力学[M].苏明德，黄素逸.计算流体力学基础[M]. Abbot,M.B.andIonescu,F.,1967.Onthenumericalcomputationofnearlyhorizontalflows.J.Hyd.Res.,5,97-117.19911（282~87，卢朝阳柯如 用Delaunay三角形化实现的举行边界表面描述算法.计算机学报，唐泽圣，徐志强 二维点集三角形剖分的动态生成与修改 计算机辅助设计与图形学学报1990，2（3：1~8，郭烈锦.非结构网格快速生成技术[J].2000，34（1：18~21，张俨彬，陶文铨.一种改进的Delaulay三角形剖分方法 石油大学学报25（2：10-SHIMADAK，GOSSARDDC.Bubblemesh:automatedtriangularmeshingofnon-manifoldgeometrybyspherepacking[C]//ProceedingsoftheThirdSymposiumonSolidModelingandApplications，NewYork,USA:ACMPress,1995:409-419.].，河流模拟[M].:水利电力，程文辉.河道二维非恒定流计算.河海大学学报，19876月39-，登，.河道平面二维水沙数学模型.中国水利水电DarwishMS.Anew-highresolutionschemebasedonthenormalizedvariableNumerHeatTransfer,PartB,1993,24:353-ChoiSK,NamHY,ChoM.Evaluationofahigher-orderboundedconvectionscheme:three-dimensionalnumericalexperiment[J].NumerHeatTransfer,PartB,1995,28:23-28JasakH,WellerHG,GosmanAD.HighresolutionNVDdifferencingschemeforarbitrarilyunstructuredmeshes[J].IntJourNumernetho