洪水淹没分析报告.docx
- 文档编号:29309396
- 上传时间:2023-07-22
- 格式:DOCX
- 页数:18
- 大小:687.69KB
洪水淹没分析报告.docx
《洪水淹没分析报告.docx》由会员分享,可在线阅读,更多相关《洪水淹没分析报告.docx(18页珍藏版)》请在冰豆网上搜索。
洪水淹没分析报告
存档编号
华北水利水电学院
NorthChinaUniversityofWaterResourcesandElectricPower
课程设计
题目:
洪水淹没分析
学院资源与环境学院
专业地理信息系统
姓名黄福想
学号200900931
指导教师李小根
完成时间2012年6月1日
目录
摘要I
关键词:
ArcEngine;三角网(TIN);ArcGIS;I
FloodAnalysisReportII
AbstractII
第1章前言-1-
1.1研究课题背景-1-
1.2研究的内容-1-
第2章系统架构-2-
2.1实现的基本原理及技术路线-2-
2.1.1.实现的基本原理-2-
2.1.2.技术路线-2-
2.2系统的实现-3-
2.2.1.应用地理信息系统的优越性-3-
2.2.2.洪水淹没面积的计算-4-
2.2.3.系统模块-4-
2.2.4.系统功能实现关键技术-7-
第3章结束语-7-
3.1系统的不足-8-
3.2核心代码:
-8-
洪水淹没分析报告
摘要
水淹没范围的确定及可视化是防洪和洪灾损失评估的核心环节,三维GIS能够在空间延展度上直观的反映地形地貌的要素特征。
基于arcengine开发的三维GIS系统通过洪水水位线和水淹范围能够形象的展示洪水淹没情况,并在此基础上分析水淹影响,统计相关的数据,计算水淹面积,为三维防洪提供一定的决策支持。
关键词:
ArcEngine;三角网(TIN);ArcGIS;
FloodAnalysisReport
Abstract
Submergedinwatertodeterminethescopeandvisualizationisacorepartoffloodcontrolandflooddamageassessment,3DGIStospatiallyextendedintuitivereflectthecharacteristicsoftheelementsofthetopography.Rangedisplaytovisualizefloodsubmergedbyfloodwaterlevelandfloodingbasedonarcenginedevelopedthree-dimensionalGISsystem,andonthisbasistoanalyzetheimpactofflooding,thestatisticaldatarelatedtocalculationofthefloodedareatoprovideacertaindegreeofdecision-makingthree-dimensionalfloodcontrolsupported....
Keywords:
ArcEngine;TIN;ArcGIS;
第1章前言
1.1研究课题背景
我国是一个自然灾害十分频繁的国家,洪涝灾害一直严重威胁着人民生命财产安全和社会的稳定与发展。
近几十年来,自然资源的开发利用不断扩大,城乡经济建设飞速发展,洪水出现的频率及其造成的损失也不断的增加。
因此,快速、准确、科学地模拟、预测和显示洪水淹没范围,以便发挥防洪工程效益,并以非工程措施来减轻洪水危害,对防洪减灾具有重要意义。
别是对于城市和蓄滞洪区,如果能够预先获知洪水的淹没范围和水深的分布情况,对于挽救人民的生命财产和减少洪涝灾害损失都具有十分重要的价值。
获取洪涝灾害范围和该范围内相对准确的水深分布对洪涝灾害评估和经济损失评估也具有重要的意义。
1.2研究的内容
洪水淹没是一个很复杂的过程,受多种因素的影响,其中洪水特性和受淹区的地形地貌是影响洪水淹没的主要因素。
对于一个特定防洪区域而言,洪水淹没可能有两种形式,一种是漫堤式淹没,即堤防并没有溃决,而是由于河流中洪水水位过高,超过堤防的高程,洪水漫过堤顶进入淹没区;另一种是决堤式淹没,即堤防溃决,洪水从堤防决口处流入淹没区。
无论是漫堤式淹没还是决堤式淹没,洪水的淹没都是一个动态的变化的过程。
针对目前防洪减灾的应用需求,对于洪水淹没分析的要求可以概化为两种情况,一是在某一洪水水位条件下,它最终会造成多大的淹没范围和怎样的水深分布,这种情况比较适合于堤防漫顶式的淹没情况。
另外一种情况是在给定某一洪量条件下,它会造成多大的淹没范围和怎样的水深分布,这种情况比较适合于溃口式淹没。
对于第一种情况,需要有维持给定水位的洪水源,这在实际洪水过程中是不可能发生的,处理的办法是可以根据洪水水位的变化过程,取一个合适的洪水水位值作为淹没水位进行分析。
对于第二种情况,当溃口洪水发生时,溃口大小是在变化的,导致分流比也在变化。
另外一般都会采取防洪抢险措施,溃口大小与分流比在抢险过程中也在变化,洪水淹没并不能自然地发生和完成,往往有人为防洪抢险因素的作用,如溃口的堵绝,蓄滞洪区的启用等。
这种情况下要直接测量溃口处进入淹没区的流量是不大可能的,因为堤防溃决的位置不确定,决口的大小也在变化,测流设施要现场架设是非常困难也是非常危险的。
所以实际应用时,考虑使用河道流量的分流比来计算进入淹没区的洪量。
第2章系统架构
2.1实现的基本原理及技术路线
2.1.1.实现的基本原理
通过对数字高程模型数据的获取和分析,建立数字高程模型(DEM)。
在建立的DEM模型的基础上,并在数据库的支持下,进行洪水淹没分析。
洪水淹没分析可以根据给定的洪水位进行计算,得到洪水淹没分析结果。
2.1.2.技术路线
利用GIS技术进行相关处理,应用自身功能,求的低于一定高程的洪水淹没范围,从而计算洪水淹没面积。
由于3D点多采用的沿等高线采点,因此在等高线密集的地方采点较多,等高线稀疏的地方采点少,这样能很好的表示地形的高低起伏,又能有效地减少数据冗余。
针对一个特定地区的洪水淹没分析,为了减少数据量和便于分析,一般根据洪水风险,预先圈定一个最大的可能淹没范围,这样处理对于防洪减灾来说是合理的。
系统的框架图如下图1所示:
图1系统框架图
2.2系统的实现
2.2.1.应用地理信息系统的优越性
1)地理信息系统是管理空间信息的最佳选择
就传统意义上的数据库而言,无论是FoxPRO,还是Oracle,Sybase,它们管理的信息类型主要是文字和数字,对图形的管理功能则十分薄弱。
众所周知,地理信息系统的数据库包括图形和属性两部分,图形数据库主要负责点(point)、线(line,arc)、面(surface)、体(body)等图形实体的管理,而各类图形实体的属性管理则由FoxPRO,Oracle等数据库来承担。
事实上,人类对地球三维空间中任意实体的表达,应该是图形和属性的并集,这样,才能对物体进行较为精确的描述。
所以,只有地理信息系统才能对空间目标进行有效的描述。
2)地理信息系统不仅是一个管理系统,而且还是一个分析系统
对地理信息系统而言,不仅图形实体具有属性信息、拓扑关系,而且还有大量的分析模型。
所以,人们不仅能够得到图形信息,而且还可以图文并茂地对地学现象进行分析研究,为决策提供依据。
例如,利用GIS系统的图形和数字高程模型(DEM),可以圈定环境的污染范围和程度,得到地形的坡向图,计算矿体的储量。
2.2.2.洪水淹没面积的计算
1)给定洪水水位下的淹没分析。
选定洪水水源入口,设定洪水水位、选出洪水水位以下的三角单元,从洪水入口单元开始进行三角格网连通性分析,能够连通的所有单元即组成淹没范围,得到连通的三角单元,对连通的每个单元计算水深W,即得到洪水淹没水深分布,图3是基于三角单元格网模型在给定洪水水位条件下的淹没分析。
2)给定洪量(Q)条件下的淹没分析。
进行灾前预评估分析时可以根据可能发生的情况给定一个洪量,或者取洪水频率对应的流量的百分数。
评估分析时洪量Q值可以根据流量过程曲线和溃口的分流比计算得到,有条件的地方,可以实测,不能实测的可以根据上下游水文站的流量差,并考虑一定区间来水的补给误差计算得到。
2.2.3.系统模块
系统的开始界面
图2系统的开始界面
TIN是地表或空间现象起伏变化的数字表示,包含了地球表面许多有用的地理信息,其中就包括水文信息。
通过TIN来生成地表水淹没分析统计模型具有最形象化的特点如下图所示:
图3淹没分析
系统统计表单模块如下图4所示:
图4系统表单
2.2.4.系统功能实现关键技术
该系统中三维场景都是通过在TIN的表面进行内插实现的。
对于场景中的水位线,它是通过数据文件中的记录的最高水位线,通过接口查询功能在contournngrid22获取相应的数据,从这些查询到的特征数据中来获取几何,为内插提供几何数据。
场景中现实的洪水淹没范围是由已有的二维的数据在表面内插所得。
图1系统淹没分析
第3章结束语
基于DEM的洪水淹没分析可以解决上述两种洪水最终淹没范围和水深分布的问题,但由于DEM数据量大,对于较大范围的洪水淹没分析,在目前的计算机硬件技术水平上还不能较快地计算出结果,这对于防洪减灾决策实施等方面,这种计算速度是不能忍受的。
格网模型的思想很早就已经提出,并且在各个领域得到广泛的应用,如有限元计算的离散单元模型,目前所能见到的较先的洪水模拟演进模型(如陆吉康教授的水动力洪水演进模型)也是一种格网化的模型。
基于空间展布式社会经济数据库的洪涝灾害损失评估模型也是基于格网模型的思想(见李纪人等人的相关文章)。
由于格网本身对模型概化的优越性,同时考虑到与洪水演进和洪涝灾害损失评估模型更好地结合,所以采用基于格网的洪水淹没分析模型是比较好的选择。
3.1系统的不足
没能做成web版发布到网上,有些功能不够完善
3.2核心代码:
publicpartialclassForm1:
Form
{
publicForm1()
{
InitializeComponent();
}
privateI3DPropertiesp3D=null;//I3DProperties控?
制?
访¤?
问¨º三¨y维?
特¬?
性?
的Ì?
接¨®口¨²
IFeaturepFeature=null;
IFeatureLayerpLayer=null;
intindex;
ITinptin;
DataRowdr;
DataTabledt=newDataTable();
privatevoidForm1_Load(objectsender,EventArgse)
{
stringpath=Application.StartupPath;
//MessageBox.Show(path);
path=path+"\\图ª?
.sxd";
this.axSceneControl1.LoadSxFile(path);
textBox1.Text="5";
textBox2.Text="5";
this.MouseWheel+=newSystem.Windows.Forms.MouseEventHandler(axSceneControl_OnMouseWheel);
}
//实º¦Ì现?
滚?
轮?
缩?
放¤?
privatevoidaxSceneControl_OnMouseWheel(objectsender,MouseEventArgse)
{
try
{
System.Drawing.PointpSceLoc=axSceneControl1.PointToScreen(this.axSceneControl1.Location);
System.Drawing.PointPt=this.PointToScreen(e.Location);
if(Pt.X
{
return;
}
doublescale=0.2;
if(e.Delta<0)scale=-0.2;
ICamerapCamera=axSceneControl1.Camera;
IPointpPtObs=pCamera.Observer;
IPointpPtTar=pCamera.Target;
pPtObs.X+=(pPtObs.X-pPtTar.X)*scale;
pPtObs.Y+=(pPtObs.Y-pPtTar.Y)*scale;
pPtObs.Z+=(pPtObs.Z-pPtTar.Z)*scale;
pCamera.Observer=pPtObs;
axSceneControl1.SceneGraph.RefreshViewers();
}
catch
{
}
}
doublezz=1100-888.216895;
privatevoidtimer1_Tick(objectsender,EventArgse)
{
if((double)(pFeature.get_Value(index))>1800)//高?
程¨¬值¦Ì超?
过y1800则¨°山¦?
谷¨¨被À?
淹¨ª没?
{
pFeature.set_Value(index,1100);
toolStripStatusLabel1.Text="淹¨ª没?
完ª¨º成¨¦";
timer1.Enabled=false;
}
else
{
pFeature.set_Value(index,double.Parse(pFeature.get_Value(index).ToString())+Convert.ToDouble(textBox2.Text));//高?
程¨¬值¦Ì加¨®10
p3D.BaseOption=esriBaseOption.esriBaseExpression;
p3D.BaseExpressionString=pFeature.get_Value(index).ToString();
p3D.Apply3DProperties(pLayer);
axSceneControl1.SceneGraph.RefreshViewers();
toolStripStatusLabel1.Text="洪¨¦水?
淹¨ª没?
进?
行D中D......";
ISurfacepSurface=ptinasISurface;
zz=zz+Convert.ToDouble(textBox2.Text);
doublearea=pSurface.GetSurfaceArea(zz,esriPlaneReferenceType.esriPlaneReferenceBelow);
doublevolum=pSurface.GetVolume(zz,esriPlaneReferenceType.esriPlaneReferenceBelow);
dr=dt.NewRow();
dr["高?
程¨¬(m)"]=zz.ToString();
dr["面?
积y(m*m)"]=area.ToString();
dr["容¨Y积y(m*m*m)"]=volum.ToString();
dt.Rows.Add(dr);
//dt.Rows[0]=Z.ToString();
dataGridView1.DataSource=dt;
}
}
privatevoidcheckBox1_CheckedChanged(objectsender,EventArgse)
{
if(checkBox1.Checked)
{
timer2.Enabled=true;//设¦¨¨置?
计?
时º¡À器¡Â运?
行D
}
else
{
timer2.Enabled=false;
}
}
privatevoidtimer2_Tick(objectsender,EventArgse)
{
axSceneControl1.Camera.Rotate(Convert.ToDouble(textBox1.Text));//按ã¡ä照?
给?
定¡§的Ì?
角?
度¨¨旋y转Áa
axSceneControl1.Refresh();
}
privatevoidbutton1_Click(objectsender,EventArgse)
{
try
{
stringsFileName;
SaveFileDialogpSaveFile=newSaveFileDialog();
pSaveFile.Title="保À¡ê存ä?
图ª?
片?
";
pSaveFile.Filter="BMP图ª?
片?
(*.bmp)|*.bmp|JPEG图ª?
片?
(*.jpg)|*.jpg|TIF图ª?
片?
(*.tif)|*.tif";
pSaveFile.ShowDialog();
sFileName=pSaveFile.FileName;
if(pSaveFile.FilterIndex==1)
{
axSceneControl1.SceneViewer.GetScreenShot(esri3DOutputImageType.BMP,sFileName);
}
else
if(pSaveFile.FilterIndex==2)
{
axSceneControl1.SceneViewer.GetScreenShot(esri3DOutputImageType.JPEG,sFileName);
}
MessageBox.Show("成¨¦功|保À¡ê存ä?
图ª?
片?
至¨¢:
"+sFileName);
}
catch
{
MessageBox.Show("出?
现?
错䨪误¨®返¤¦Ì回?
");
}
}
privatevoidbutton2_Click(objectsender,EventArgse)
{
IScenepScene=this.axSceneControl1.Scene;//IScene提¬¨¢供?
访¤?
问¨º和¨ª控?
制?
场?
景¡ã的Ì?
接¨®口¨²
IEnumLayerpEnumLayer=pScene.get_Layers(null,true);
pLayer=(IFeatureLayer)pEnumLayer.Next();
toolStripStatusLabel1.Text="洪¨¦水?
淹¨ª没?
进?
行D中D......";
while(pLayer!
=null)
{
//MessageBox.Show(pLayer.Name);
if(pLayer.Name=="水?
面?
")
{
IFeatureClasspFeaClass=pLayer.FeatureClass;
pFeature=pFeaClass.GetFeature(0);
index=pFeature.Fields.FindField("高?
程¨¬");
timer1.Enabled=true;
break;
}
else
{
pLayer=(IFeatureLayer)pEnumLayer.Next();
}
}
ILayerExtensionspLayerExts=pLayerasILayerExtensions;
for(inti=0;i { I3DPropertiespDDD=pLayerExts.get_Extension(i)asI3DProperties; if(pDDD! =null) { p3D=pDDD; } } ITinLayerpTinLayer; pTinLayer=axSceneControl1.Scene.get_Layer (2)asITinLayer; ptin=pTinLayer.Dataset; dt.Columns.Add("高? 程¨¬(m)",typeof(string)); dt.Columns.Add("面? 积y(m*m)",typeof(string)); dt.Columns.Add("容¨Y积y(m*m*m)",typeof(string)); //IWorkspaceFactorypWSFact=newTinWorkspaceFactoryClass(); //IWorkspacepWS=pWSFact.OpenFromFile(Application.StartupPath,0); //ITinWorkspacepTinWS=pWSasITinWorkspace; //ITinpTin=pTinWS.OpenTin("dtm_tin"); } privatevoidcheckBox2_CheckedChanged(objectsender,EventArgse) { } publicvoidaxSceneControl1_OnMouseDown(objectsender,ISceneControlEvents_OnMouseDownEvente) { } //boolflag=false; privatevoidbutton1_Click_1(objectsender,EventArgse) { timer1.Enabled=false; //if(flag==false) //{ //timer1.Enabled=false; //button1.Text="继¨¬续? "; //flag=true; //} //else //{ //timer1.Enabled=true; //flag=false; //button1.Text="停ª¡ê止1"; //} } privatevoidbutton4_Click(objectsender,EventArgse) { timer1.Enabled=true; } privatevoidbutton3_Click(objectsender,EventArgse) { try { stringsFileName; SaveFileDialogpSaveFile=newSaveFileDialog(); pSaveFile.Title="保À¡ê存ä?
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 洪水 淹没 分析 报告