洪水排险论文.docVIP

洪水排险 摘要： 洪水险情是危害人民财产安全的重要不利因素，本文针对洪水排险范方案的设计问题，在满足洪水量的存储的条件下以损失最小为方案目标，通过图论Dijkstra算法 首先根据题意整理数据并导入matlab中求解1-20各区域的面积s及p=1.2时的各区域体积，根据k/s求解出单位面积的损失量并排序。 对于无海拔因素下的排险问题，使得所有区域全部受损失的最小洪水量即是使得1-15以p=1.2为水高，16-20以1.0为水高时的洪水量，此时满足所有区域全部受损失且洪水量最小。针对泄洪最优方案的决策，通过图论Dijkstra算法 对于有海拔因素的排险方案，要充分运用由于海拔因素所导致的储水量大幅增加的条件，通过局部破堤的方式比较破堤后综合单位体积损失量与未破堤时1-20的单位体积损失，确定最优方案。 最后的加固方案可假设全部最大加固是区域体积的变化所导致的单位损失量排序的变化，之后通过matlab 最优路径求法得到泄洪最优路径方案并比较加固方案中有效加固费用及损失费加和和未加固方案中损失费。选择费用最小方案。 关键词：matlab编程 图论Dijkstra算法 1 问题重述 近段时间以来，由于我国南方地区遭受持续的强暴雨袭击，许多江河流域出现百年不遇的险情。为减少损失，将采取破堤的方法排除险情，而泄洪必然影响下游沿岸区域。 位于广西 该区域内又分成20个矩形小区。下图是该区域的示意图： 该区域与河流之间有大堤隔离，各小区域之间有高度为1.2m的小堤相互隔离，假设决堤口可选在大堤或小堤的任何地方，决堤口数目不受限制。但一经决口，就不能再补合。从河流经大堤决口流入小区域的洪水量按决口数成比例分配。如在小区域之间小堤开一决口，则假设该两小区域之间的这段小堤不复存在。若水位高于小堤，则将自动向邻近最低的一个小区泄洪。若这样的小区有几块时，则平均泄洪。 各小区域被完全淹没时土地、房屋和财产等损失总数k（单位：百万元）如下表： 小区域 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 损失 1.4 7.0 5.8 3.3 2.0 9.4 0.9 6.0 7.2 1.6 3.0 4.1 4.1 5.3 4.4 1.0 1.1 2.1 1.3 4.2 当洪水淹没一个小区域且水位高于该小区域高度p（单位m）时，该小区域的损失为该小区域的k和p的函数。 损失= 1、假定各小区域的海拔高度相同 2、假定各小区域的海拔高度不同，各小区域的海拔高度h（单位m）如下表： 小区域 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 海拔 3.6 4.0 4.7 4.4 3.8 3.3 3.2 2.5 5.0 4.4 3.0 3.5 2.4 3.8 3.8 3.7 2.8 2.5 3.0 4.8 注：此时各小区域间的小堤指相对高度，例如小区域1、2间事实上有海拔5.2的小堤。 你要做的事情是： 一、在海拔高度相同的假设下，求： 1、整个区域全部受损失的最小洪水量； 2、当洪水量为20*106，40*106（m3）时，分别制定泄洪方案，使总损失最小（在一种方案中，决堤同时进行）。须计算出该方案的损失数。 二、在海拔高度不同的假设下，求： 1、整个区域全部受损失的最小洪水量； 2、当洪水量为40*106，80*106（m3）时，分别制定泄洪方案，使总损失最小。 三、在海拔高度不同的假设下，可以对小堤增高加固，增高加固的费用（不计厚度）为：1（百元/平方米），最高可将原有小堤提高0.5（m），制定泄洪增高加固方案，使总损失最小。 2.问题的分析 确定洪水排险方案的依据是损失量最小，因此在初始整理数据时要首先对各区域单位体积损失进行排序，同时要在区域均可破堤的连通路径假设下通过Dijkstra算法 而在有海拔高度因素的方案决策中，突破点应为由于海拔因素所导致的局部区域储水量显著增加，进而单位体积损失显著减少，此时单位体积损失排序也会随时发生巨大变化，当重新排序过后，再次通过Dijkstra算法 在加固方案选择的问题中，加固导致的实质改变因素是区域的体积，进而是单位体积损失的改变，因此依然通过排序新的单位体积损失量，并推导最优路径，线性规划甄选，模拟演绎检验来确定最优方案。 3.模型的假设 1． 假设题目中所给出的数据真实可信。 2． 假设每个小区域的地面是水平的，没有凸凹不平现象。 3． 假设区域全部受损失指该区域被淹没1米以及1米以上。 4． 假设在各个区域海拔高度相同时，打通相邻区域间的小堤后，相邻区域内洪水水位相同。 5． 假设在未人为打通堤坝的情况下，大堤和小堤均不会自行决口或溃堤。 6． 假设大堤各个决口处流入洪水量相同