内容文本文稿2009北京9004521.pdf

全国第五届研究生数学建模竞赛 全国第五届研究生数学建模竞赛 题 目 唐家山堰塞湖泄洪问题的研究 摘 要 本文研究的是唐家山 的堰塞湖问题，结合数字高程 地图和 所提供的数据分别建立模型研究了唐家山堰塞湖形成之后 湖水高程、蓄水量、溃坝情形、溃坝灾害等一系列问题。 针对第一问，首先对数字高程地图进行等高图像分析求解了堰塞湖不 同高程水位对应的湖区面积，建立了蓄水量体积与堰塞湖水位高程的离散 化模型，然后建立了神经网络模型和多元线性回归模型研究了北川天气预 报 3 雨量与堰塞湖入库流量的关系，继而求解得到不同降雨量下 堰塞湖水位高程。针对第二问，首先研究泄洪过程和溃坝过程的区别与联 系，从而分别建立正交多项式逼近和 模型得到溃坝时的溃口流量随时 间变化的关系，继而分析求解得到溃坝时溃口宽度、深度、水位高程和水 流速度随时间变化的关系。针对第三问，综合数字高程地图和行政区域地 图，在数字地图中查找地势相对较低区域,进而得到洪水下泄区域及被淹没 区域。确定洪水的流速之后，利用数字地图计算了洪水到达各被淹没区域 的时间，淹没范围，以及淹没之后的安全区域，并据此制订了初步的撤离 方案。针对第四问，建立安全最优化模型定量评价了决策部门的堰塞体导 流渠掘进方案并定性分析了决策部门制定的撤离方案。 关键词：图像分析 神经网络 正交拟合 溃坝 参赛队号 9004521 参赛 （由组委会填写） 1 问题重述 1.建立唐家山堰塞湖以水位高程为自变量的蓄水量的数学模型。并以该地区天气预 报的降雨情况的50%，80%，100%，150%为实际降雨量建立模型预计自5 月25 日起至6 月12 日堰塞湖水位 上升的高度(不计及泄洪)。 2.利用这次唐家山堰塞湖泄洪时科技人员记录下了大量宝贵的数据，在合理的 下建立堰塞湖蓄水漫顶后在水流作用下发生溃坝的数学模型，包含缺口宽度、深度、水 流速度、水量、水位高程，时间等变量。 3.根据数字地图给出坝体发生溃塌，造成堰塞湖内1/3 的蓄水突然下泻时的洪水水 流速度及淹没区域（包括洪水到达各地的时间），并在此基础上考虑洪水淹没区域中人 口密集区域的人员撤离方案。 4.请根据建立的数学模型分析当时所采取对策的正确性和改进的可能性。讨论为应 对 后次生山地灾害 (不限堰塞湖) ，科技工作中应该设法解决的关键问题，并提出 有关建议。 2 问题分析 1.研究不同高程水位对应的堰塞湖蓄水量，即研究不规则体的体积问题，可以考虑 引入面积对高度的微元进行积分这一思想进行求解。 中的实际数值可以作为求 解时的重要 和校正点。在研究降雨量和堰塞湖水位的关系之时，认为将这一问题转 化为首先研究降雨量和库容变化量之间的关系更为直观，考虑到两者关系的复杂性，并 不一定需要建立明确的函数关系式。 2.现实中唐家山并没有发生溃坝，无法获得溃坝的直接数据。这里也十分缺乏利用 成型模型求解溃坝问题的大量参数。但是唐家山的泄洪的规模也是非常大的，可以借助 研究唐家山的泄洪过程分析这一泄洪和溃坝的区别与 ，继而通过拟合、 等 得出唐家山的溃坝过程。 3.结合数字地图可以求得不同位置的海拔高度，根据水顺势而流的原理就可以推算 洪水的流经区域。结合泻出的洪水量就可以结合问题1 的思想计算求解淹没范围。 3 重要符号说明 V −−堰塞湖蓄水量； S −−堰塞湖实际剖面面积； SK −−堰塞湖图像面积； h −−水位高程； w −−溃口宽度值； g −−溃口深度值； q −−水流量大小； f −−溃口中水深； d −−掘进导流槽 高度； K −−面积系数即堰塞湖实际面积与图像面积的比例系数