数学建模课程培训(二)解析.pptx

数学建模中的物理 1 第五届全国研究生数学建模竞赛 A题：汶川地震中唐家山堰塞湖泄洪问题 堰塞湖，是由于地震而造成的山体滑坡，堵截河谷或河床后贮水而形成的湖泊。如果遇到强余震、暴雨，可能会发生溃坝，对下游百姓的生命财产造成威胁。同时，由于堰塞湖水位不断上升，也会对上游造成淹没的危险 针对唐家山堰塞湖泄洪问题，加强对震后次生灾害规律的研究，为国家抗震救灾提供更有力的科学支撑 2 问题的提出 1.建立唐家山堰塞湖以水位高程为自变量的蓄水量的数学模型。并以该地区天气预报的降雨情况的50%，80%，100%，150%为实际降雨量建立模型预计自5月25日起至6月12日堰塞湖水位每日上升的高度(不计及泄洪)。 2.在合理的假设下建立堰塞湖蓄水漫顶后在水流作用下发生溃坝的数学模型，建议包含缺口宽度、深度、水流速度、水量、水位高程，时间等变量。 3 问题的提出 3. 根据数字地图给出坝体发生溃塌，造成堰塞湖内1/3的蓄水突然下泻时的洪水水流速度及淹没区域（包括洪水到达各地的时间），并在此基础上考虑洪水淹没区域中人口密集区域的人员撤离方案。 4.请根据你们所建立的数学模型分析当时所采取对策的正确性和改进的可能性。讨论为应对地震后次生山地灾害 (不限堰塞湖) ，科技工作中应该设法解决的关键问题，并提出有关建议。 4 堰塞湖基本情况 唐家山堰塞坝位于北川县城上游6km的通口河上，下距苦竹坝约1km 集水面积3550km2 堰塞湖总容积约 3.15亿立方米 堰塞坝坝体长803.4m，宽611.8m，面积30.72万平方米，体积2037万立方米 坝址区通口河为不对称的V型河谷，右岸较陡，坡度60°左右；左岸坡较缓，坡度约为30° 坡体右侧、中部和左侧分布3条沟槽，右侧沟槽最高点高程752.2m，宽20～40m，贯通上下游。堰塞坝最高点高程793.9m（黄海高程，下同）；堰塞坝底高程667.4m。 5 问题一分析： 建立唐家山堰塞湖以水位高程为自变量的蓄水量的数学模型。 6 堰塞湖示意图 W：堰塞湖内蓄水量，即总库容，单位：亿立方米 Hb0：堰塞湖底部高程，常数(667.4米) H(t)：坝前水位高程，单位：米 h(t)：堰塞湖内水深，单位：米 b(t)：泄流槽的宽度，单位：米 ？ 7 数据采样 坝前水位 714.01 713.41 713.42 713.44 713.49 713.5 713.54 713.57 蓄水量 0.857 0.834 0.835 0.835 0.837 0.838 0.839 0.840 坝前水位 713.64 713.68 713.79 714.13 714.51 717.51 718.06 719.48 蓄水量 0.843 0.844 0.849 0.860 0.876 1.000 1.020 1.120 坝前水位 720.25 725.80 727.94 732.51 737.56 740.51 742.80 742.96 蓄水量 1.139 1.400 1.450 1.700 2.071 2.328 2.469 2.486 坝前水位 742.80 742.58 742.11 741.66 741.12 740.55 739.59 739.37 蓄水量 2.475 2.475 2.420 2.389 2.351 2.307 2.237 2.222 坝前水位 738.93 738.85 738.71 738.5 737.78 737.53 736.78 736.63 蓄水量 2.190 2.185 2.175 2.159 2.116 2.094 2.055 2.008 来自附件中《中国水利一线记者》和《水信息网唐家山堰塞湖抢险全记录》 蓄水量单位：亿立方米 坝前水位单位：米 8 数据采样分析 根据报道中搜集的所有数据，可以基本确定总蓄水量W 与水深hL 存在指数关系 9 模型拟合 模型检验 因为最高点高程为最高点高程752.2m，代入上式可得堰塞湖最大库容量为： 基本资料中给出的唐家山堰塞湖总容量3.15亿立方米，模型与实际的绝对误差为0.17亿立方米，相对误差为5.4% 10 问题1.2 以该地区天气预报的降雨情况的50%，80%，100%，150%为实际降雨量建立模型预计自5月25日起至6月12日堰塞湖水位每日上升的高度。 分析问题 研究降雨量对唐家山堰塞湖水位的影响主要体现在唐家山堰塞湖蓄水量的改变 唐家山堰塞湖蓄水量主要来源 每日正常入库的水流量 降雨后雨水入库的水流量 （假设为常数C） 11 降水量分析 前24小时直接落在堰塞湖面的降雨 周围山体的降雨经过一定的延时后汇入堰塞湖 上游地区的积雨经过一定的延时后汇入堰塞湖 直接降雨 量 设为Ra(t)，假设堰塞湖面积为S1，前24小时的降雨量为J1(t