10队_最优洪水排险方案问题.docxVIP

最优洪水排险方案问题[十队]参赛队员姓名学 号学院及专业年级联系电话李融学院∕金融工程201015828619680廖菡计学院∕金融统计201013550140554李重融学院∕金融理财201013550140595摘要作为一种自然灾害，洪水会对人们的生命、财产安全等造成严重损害。本文在规定条件下计算给出一系列测度数据，讨论了合理的泄洪方案。首先，由于假设地区内各小区域海拔相同，且小堤之间不开口，为求得使其全部受损的最小洪水量，先从单个区域入手分析，将某区域被“灌满”视作整个泄洪过程中的一“均衡”状态，状态的改变对应着每一区域储水量的变化。将20个区域视作位置不同的点，抽象出一赋权图，构造带权邻接矩阵，引入图论知识。邻接矩阵中的点表示某一均衡状态下从区域溢水灌满区域所需剩余水量。水量变化对应矩阵的改变，运用MATLAB编写程序，计算出满足条件的最小洪水量，同时搜索定位不同情况下所对应的洪水路径。针对第二问中泄洪方案的制定，由于同样有第一问中的假设条件，且洪水量固定，而各区域完全被淹没时的损失不同，因而应使单位体积损失小的区域尽可能承载下所有洪水。计算分析每一区域的体积损失量大小，并根据上述思想对逐个区域分析讨论，得到在给定洪水量时的最优泄洪方案。第三问由于假定各小区域的海拔高度不同，则为了求得使区域全部受损的最小洪水量，只需关注海拔最高的区域9，区域9决定了其他区域的水位高低与损失大小，此时问题容易求得。第四问与第二问类似，其思想方法是相同的，但因此时各小区域海拔不同，所以有必要进一步讨论不同区域储水量与损失的关系，也即因为每个区域的实际有效围墙高度是不同的，所以需要计算出新的有效体积损失率。通过明确区域间相对高度，容易计算出有效围墙高度，从而易求得有效体积损失率，进而采用第二问的方法进行查找推算，通过MATLAB编程求得最终精确结果与合理泄洪方案。针对问题四中各小区域小堤固定不变的问题，问题五考虑了小堤可以增高加固的情况。本文在问题五中所给模型的基础上，逐步考虑各个区域与四周区域的小堤能否加固，在满足泄洪要求并考虑总损失和费用最小的前提下，得到了洪水量分别为40*、80*时最佳加固方案和泄洪方案，得到的总损失分别为7.1901百万元和23.50617百万元。关键词：邻接矩阵体积损失率泄洪方案一、问题重述近段时间以来，我国南方地区遭受持续的强暴雨袭击，许多江河流域出现百年不遇的险情。为减少损失，将采取破堤的方法排除险情，而泄洪必然影响下游沿岸区域。位于广西省的某区域三面环山，山很高，一面是河流，用堤坝隔离，使该区域成为封闭区域。当上游发生泄洪时，洪水极易破堤淹没区域，造成人员和财产的损失。为减少总的损失，人们采取破堤排除险情。该区域内又分成20个矩形小区，题目给出该区域的示意图。该区域与河流之间有大堤隔离，各小区域之间有高度为1.2m的小堤相互隔离，假设决堤口可选在大堤或小堤的任何地方，决堤口数目不受限制。但一经决口，就不能再补合。从河流经大堤决口流入小区域的洪水量按决口数成比例分配。如在小区域之间小堤开一决口，则假设该两小区域之间的这段小堤不复存在。若水位高于小堤，则将自动向邻近最低的一个小区泄洪。若这样的小区有几块时，则平均泄洪。各小区域被完全淹没时土地、房屋和财产等损失总数k（单位：百万元）如题目中表格所示。当洪水淹没一个小区域且水位高于该小区域高度p（单位m）时，该小区域的损失为该小区域的k和p的函数（题目已给出）。各小区域海拔高度（单位m）假设为相同与不同两种情况，在后一种情况中，题目给出各小区域海拔高度资料。此时各小区域间的小堤指相对高度。建立数学模型分析研究以下问题：一、在海拔高度相同的假设下，求：1、整个区域全部受损失的最小洪水量；2、当洪水量为20*106，40*106（m3）时，分别制定泄洪方案，使总损失最小（在一种方案中，决堤同时进行）。须计算出该方案的损失数。二、在海拔高度不同的假设下，求：1、整个区域全部受损失的最小洪水量；2、当洪水量为40*106，80*106（m3）时，分别制定泄洪方案，使总损失最小。三、在海拔高度不同的假设下，可以对小堤增高加固，增高加固的费用（不计厚度）为：1（百元/平方米），最高可将原有小堤提高0.5（m），制定泄洪、增高加固方案，使总损失最小。二、问题分析第一问要求在各小区域海拔高度相同的情况下，求解出使得整个区域全部受损失的最小洪水量。假设小区域之间没有人工开口破堤的情况，则洪水将渐次填充每一区域，直至所有区域全部受损。将某区域或片区（几个区域均灌满后可视作一个片区）被灌满时视作一均衡状态，每一状态下各区域水量不同，因而将小区域抽象为不同的点，得到以变化储水量为权的赋权图，构造邻接矩阵，用图论知识加以描述求