专利名称：基于环方程法和管网分区的并行水力计算子区域构造方法
技术领域：
本发明涉及一种基于环方程法和管网分区的并行水力计算子区域构造方法，属于市政工程技术领域。
背景技术：
目前，基于区域分解的并行计算方法所面临的一个主要问题，就是如何将分析对象有效地分解为计算负荷相当的子区域，以实现真正高效的并行计算。对于给水管网并行水力计算，现有的方法是运用数学计算方法和基于图论的分析，对用来描述管网构造的稀疏矩阵进行结构变换，来使其满足并行计算的需要。可见，这种分解方法并非从分析实际管网构造出发，因此，得到的子区域并不一定具有明确的物理意义。对于分区改造后的管网， 虽然可直接将其物理功能分区作为并行计算子区域，但同样存在局限性。原因在于，首先， 物理功能分区的划分还没有成熟的理论方法，多依靠经验而定，非常耗时耗力，所以已采用物理功能分区划分的管网数量非常有限；其次，物理功能分区的划分是从改善给水管网管理的目的出发，如便于分级用水量计量、优化调度、漏失控制等，未考虑优化水力计算的需求，所以物理功能分区的规模，即包含节点和管段数量通常不一致，即使规模类似，分区的数量也会成为限制因素，很难恰到好处的适合于并行计算。因为物理功能分区的规模，即包含节点和管段数量通常不一致，即使规模类似，分区的数量也会成为限制因素，很难恰到好处的适合于并行计算。分区数量过多时，多计算机的数量需求过高，增加了系统复杂度和成本；分区数量过少时，每个区域的复杂度过高，同样无法实现高效并行计算。
发明内容
针对上述问题，为了简化并行计算子区域的构造方法，本发明提出了一种基于环方程法和管网分区的并行水力计算子区域构造方法。本发明以已经按照拓扑结构演化进行分区的给水管网为分析对象，旨在对这些以满足管网运行管理为目的物理功能分区进行二次划分，构造能高效服务于基于环方程法的并行水力计算的子区域。此类子区域可定义为并行计算分区，只存在于水力模型中而非对管网进行真实的改造，其构造过程主要涉及确定并行计算机的理想运算负荷；依据运算负荷对现有物理功能分区进行合并与拆解；二次划分结束后所有子区域序号重新排列。这种方法有效地利用了给水管网的物理功能分区结构，为并行计算中的子区域划分策略提供了有力依据；同时使得管网分区后的单独模拟成为可能，且操作简单。由于并行水力计算和管网分区都是未来大型给水管网管理的关键技术，因此，本方法作为将两者高效结合起来的技术，具有很好的应用前景。本发明的技术方案如下一种基于环方程法和管网分区的并行水力计算子区域构造方法，其特征在于所述方法具体步骤如下(1)根据并行计算机的配置条件，确定理想的运算负荷；(2)由运算负荷估算每个并行计算子区域所能容纳的节点数量和管段数量总和N的容限值R，Nfflin彡R彡Nfflax ；(3)统计给水管网各个物理功能分区中节点和管段的数量Np ；(4)进行二次划分，将&低于容限值的区域进行合并，而大于容限值的区域则进一步分解，最终得到计算负荷彼此相当的并行计算子区域；(5)标识耦合区域将各个子区域间的连接管道作为耦合关联，耦合关联与各子区域边界处的管道若相互连接且能形成环路，则将其作为耦合区域，加以标识；(6)编号重排将所有子区域和耦合区域中的环路进行编号，编号时首先依次对各个子区域内部的所有环路进行顺序编号，最后再依次对各个耦合区域中的所有环路进行
编号；(7)构造系数矩阵以步骤(6)的编号作为系数矩阵新的行号和列号，进行相应的变形，即可生成符合区域分解并行计算要求的系数矩阵，使并行水力计算得以实现。与现有技术相比，本发明具有以下优点(1)并行计算子区域的划分具有明确的参照对象和准则；(2)每个单独运算的子区域都是实际管网中的一个物理功能分区或其一部分，或几个物理功能分区的合并，因此具有明确的物理意义，便于参照和分析；(3)子区域的序号重排简便易行，大幅简化稀疏矩阵的结构转换过程，且不产生注入元；(4)有效地将管网分区结果与并行计算结合起来。
图1是本发明基于环方程法和管网分区的并行水力计算的环路编号方法示意图。图中1-39为编号。
具体实施例方式如图1所示，本发明的一种基于环方程法和管网分区的并行水力计算子区域构造方法，具体步骤如下(1)根据并行计算机的配置条件，确定理想的运算负荷；(2)由运算负荷估算每个并行计算子区域所能容纳的节点数量和管段数量总和N
的容限值R，Nfflin彡R彡Nfflax ；(3)统计给水管网各个物理功能分区中节点和管段的数量Np ；(4)进行二次划分，将&低于容限值的区域进行合并，而大于容限值的区域则进一步分解，最终得到计算负荷彼此相当的并行计算子区域；(5)标识耦合区域将各个子区域间的连接管道作为耦合关联，耦合关联与各子区域边界处的管道若相互连接且能形成环路，则将其作为耦合区域，加以标识；(6)编号重排将所有子区域和耦合区域中的环路进行编号，编号时首先依次对各个子区域内部的所有环路进行顺序编号，最后再依次对各个耦合区域中的所有环路进行编号。编号实例示意请参见图1。(7)构造系数矩阵以步骤(6)的编号作为系数矩阵新的行号和列号，对原矩阵中的各元素进行重新排列，例如，原编号为1的环路所对应的元素排在矩阵的第一行第一列，而重新编号后若此环路编号变为7，则其对应元素的位置相应移至矩阵的第7行第7列，即可生成符合区域分解并行计算要求的系数矩阵，使并行水力计算得以实现。
权利要求
1. 一种基于环方程法和管网分区的并行水力计算子区域构造方法，其特征在于所述方法具体步骤如下(1)根据并行计算机的配置条件，确定理想的运算负荷；(2)由运算负荷估算每个并行计算子区域所能容纳的节点数量和管段数量总和N的容限值R，Nfflin彡R彡Nfflax ；(3)统计给水管网各个物理功能分区中节点和管段的数量Np；(4)进行二次划分，将Np低于容限值的区域进行合并，而大于容限值的区域则进一步分解，最终得到计算负荷彼此相当的并行计算子区域；(5)标识耦合区域将各个子区域间的连接管道作为耦合关联，耦合关联与各子区域边界处的管道若相互连接且能形成环路，则将其作为耦合区域，加以标识；(6)编号重排将所有子区域和耦合区域中的环路进行编号，编号时首先依次对各个子区域内部的所有环路进行顺序编号，最后再依次对各个耦合区域中的所有环路进行编号；(7)构造系数矩阵以步骤(6)的编号作为系数矩阵新的行号和列号，对原矩阵中的各元素进行重新排列，即可生成符合区域分解并行计算要求的系数矩阵，使并行水力计算得以实现。
全文摘要
本发明涉及一种基于环方程法和管网分区的并行水力计算子区域构造方法，属于市政工程技术领域。其构造过程主要涉及确定并行计算机的理想运算负荷；依据运算负荷对现有物理功能分区进行合并与拆解；二次划分结束后所有子区域中的环路序号重新排列。这种方法有效地利用了给水管网的物理功能分区结构，为并行计算中的子区域划分策略提供了有力依据；同时使得管网分区后的单独模拟成为可能，且操作简单。由于并行水力计算和管网分区都是未来大型给水管网管理的关键技术，因此，本方法作为将两者高效结合起来的技术，具有很好的应用前景。
文档编号G06Q50/00GK102567114SQ201010599439
公开日2012年7月11日 申请日期2010年12月21日 优先权日2010年12月21日
发明者付亚平, 刁克功, 刘姗姗, 刘子龙, 叶婉露, 周玉文, 常胜昆, 曾玉蛟, 杨小艳, 王昊, 王正吉, 翁窈瑶 申请人:北京工业大学