通用河网二维水流模拟模式的研究.pdfVIP

通用河网二维水流模拟模式研究 王船海，向小华 （河海大学 水文水资源及水利工程科学国家重点实验室，南京210098） 摘 要：通过对流域河网二维水流模拟问题的分析，提出了一个通用化的河网二维模拟概化 模式，该模式将河网二维概化为：“树状”河网计算单元、“环状”河网计算单元、“十字型” 河网计算单元等河网二维计算单元的组合，通过这些河网计算单元的求解，实现全流域型河 网一二维耦合模式的求解模型，并籍由此建立了通用化的河网二维水流数学模型软件系统。 通过 “树状”河网二维计算单元求解模型的推导与实现，验证了概化模式的可行性，结果表 明该思路可以应用于大型流域河网一、二维水流模拟，具有很好的实用价值。 关键词：流域河网二维水流模型，“树状”河网计算单元、“环状”河网计算单元、“十字型” 河网计算单元、河网二维计算单元，面向对象 由于数学模型理论及计算技术的飞速发展，对于流域型大范围一维河网水流的数值模拟 [1],[2] 已经很成熟 ，并且广泛应用于各行各业的实际问题中。河道二维水流模拟计算的研究， 在近二十年来也已日趋成熟，主要有ADI[3]法、SIMPLE[4]及“矩阵追赶法”[5][6]等算法，这 些算法成功的解决了一些实际问题,但这些问题一般局限在局部小范围区域，且以单一河道 模拟计算为多。随着实际应用需求的不断增加，迫切要求研究出流域性范围内的河网二维水 流模拟以及实现全流域河网一、二维耦合算法。 1． 问题的提出 河网二维水流模拟计算除了单一河道二维模拟所存在难度[6][7] 外，还由于河网二维是由 若干个单一河道组成，河道分汊口处的模拟精度对于整个区域的精度起着至关重要的作用， 尤其是当关注的重点在河口区域处时，因此如何解决河道分汊口的模拟问题是河网二维水流 模拟的关键所在。目前常用的一些河道二维水流模拟模型算法一般主要是用在单一河道或单 一分汊口的河网二维模拟，对于多分汊口的河网二维模拟研究不多；对于多分汊口的河网二 维模拟主要有如下几方面的难点： （1）多分汊河口区域的模拟问题 目前常用的河道二维算法中，ADI法由于 “轴化现象” 的限制[8]，相应采用的时间步长也很短，需要很多的计算机时，效率不高，且精度也难以满 足要求，因此ADI法一般的仅仅适用于宽广水域。SIMPLE算法实际上是一种迭代方法，对 由多河口组成的河网二维计算问题，由于河道分汊口处的流态复杂，因而其迭代不易收敛， 对于求解全流域型的河网问题迭代这种算是不可行的。非分步的全耦合的矩阵追赶法[5][9] ， 在河道的分汊口的处理上有全隐连接、半隐连接等方法，基本解决了单分汊口的河网二维问 题，但对于多分汊河口目前还无能为力。 （2 ）河网一、二维水流模拟的耦合问题 河道二维水流计算的上下边界一般的情况均是 上下游水位已知的边界条件，由于计算范围的问题，往往是通过一维模型的计算提供。对于 - 1 - 流域型河网二维水流计算问题，由于边界众多，更需要一维河网模型提供计算边界，这就提 出了一个问题，能否实现一二维模型的耦合，希望通过一二维的耦合可以解决河网二维水流 计算的边界条件问题，同时还可以解决边界条件的类型问题（流量边界条件问题），目前的 二维水流模拟算法是无法解决上述问题。 本文的目的通过对流域型河网二维水流模拟问题的分析，提出通用化的河网二维水流模 拟算法，解决上述难点问题。 2 ． 流域河网二维水流计算思路 流域河网一维水流模拟计算中，是以一条河道作为最小计算单元, 以水位节点方程作为流 域控制方程，其求解的关键是如何构建流域节点水位方程，而在构建该方程时的关键是获得 了河道 （最小计算单元）的首末断面流量与河道首末节点水位的线性关系[1] 。采用类似的思 想，将流域需求解的河网二维区域概化为由若干个河网二维计算单元的组合，对于河网二维 计算单元