JP6520253B2

JP6520253B2 - 配水管網計算装置及び配水管網計算方法

Info

Publication number: JP6520253B2
Application number: JP2015053007A
Authority: JP
Inventors: 貴雅藤澤; 孝則林; 達斎外山
Original assignee: Meidensha Corp
Current assignee: Meidensha Corp
Priority date: 2015-03-17
Filing date: 2015-03-17
Publication date: 2019-05-29
Anticipated expiration: 2035-03-17
Also published as: JP2016173704A

Description

本発明は、上水道の配水管網の水理計算を行う分野の技術であり、特に、管網データの配水管の接続点（以下、節点）が多くなった場合の水理計算を行う技術に関する。

配水管網計算装置は自治体などが管理する上水道設備において上水道の配水管網の水理解析を行う装置である。この配水管網の水理解析を行うことにより、配水状態等をシミュレーションでき、設備機器の運転状態を変えたときの影響や、設備を追加したときの影響などを事前に検討することができる。そして、これにより、上水道設備の運転支援やエンジニアリングを行うことができる。

水理解析を行う方法としては特許文献１に開示された計算方法が挙げられる。この方法は配水管網の配水管内の流量を表す流量公式と当該配水管網の節点における流量の収支式（以下、節点方程式）とを適用したものである。特に、流量公式には配水管網計算法において最も一般的なヘーゼン・ウィリアムズ式が適用されている。これらの式により流量と節点エネルギー位を導くのが配水管網計算であって管路数及び節点数の非線形連立方程式を解くことが必要となっている。

また、特許文献１の計算方法では、流量公式をニュートン法によって線形近似を行い、この流量公式を節点方程式に代入してこの方程式の解として節点エネルギー位を算出する。

前記流量公式は下記の式となる。
Ｒ_ij＝０．２７８５３Ｃ×Ｄ^2.63×Ｌ^-0.54
ｑ_i＝Ｒ_ij（ｈ_i−ｈ_j）^0.54
Ｑ_ij＝ｑ_i＋０．５４Ｒ_ij｜ｈ_i−ｈ_j｜^-0.46×（（Ｈ_i−ｈ_j）−（Ｈ_j−ｈ_j））…（０）
Ｑ_ij：節点ｉ，ｊ間の管路流量
Ｄ：管径
Ｌ：管路長
Ｈ_i：節点ｉのエネルギー位
Ｈ_j：節点ｊのエネルギー位
ｈ_i，ｈ_j：Ｈ_i，Ｈ_jの仮定値
式（０）を節点方程式に代入し、節点エネルギー位を解とする連立方程式を解き、節点の収支式における誤差が所定の精度以下であれば水理計算を終了する。

特開平２−２６４３０８号公報

金田康正著，「並列数理処理シリーズ９並列数値処理−高速化と性能向上のために−」，コロナ社，２０１０年４月３０日，初版，ｐ．９７

少ない節点数であれば、特許文献１の計算手法のように、節点エネルギー位の連立方程式を解くことは可能である。しかし、節点が多くなってきた場合に計算が遅くなる。また、節点エネルギー位の連立方程式を解く場合、行列方程式を用いて解くことになるが、ほとんどの行列要素はゼロであり、計算的には無駄な処理を行うことが多い。

そこで、前記行列方程式において、非ゼロ要素のみをとり扱うことにより、処理を高速で行うための方法として、非特許文献１のBi-CGSTAB法（安定化双共役勾配法）が知られている。この方法は繰り返しの演算によって行列方程式の解を算出する。

また、式（１）において、ｈ_i，ｈ_jの初期値として「ゼロ」等の値を用いる。しかし、ｈ_i，ｈ_jの初期値として「ゼロ」を与えると、式（１）の｜ｈ_i−ｈ_j｜^-0.46の計算結果は無限大となり、計算を行うことができない。したがって、特許文献１の計算方法において、｜ｈ_i−ｈ_j｜の計算結果が「ゼロ」になる場合、微小な値（例えば１０^-5等）を用いて計算を行う必要がある。

Bi-CGSTAB法を用いて行列方程式を解く場合、｜ｈ_i−ｈ_j｜の計算結果が「ゼロ」になるときに微小な値（例えば１０^-5等）を用いると、節点の収支式の誤差が所定の許容値以下にはならない。これは、Bi‐CGSTAB法で繰り返し演算を行い、算出したエネルギー位のｈ_i，ｈ_jの計算結果がほぼ等しい場合、１０^-5の値で端数処理されるため、１０^-5以上の精度が出なくなるからである。

ここで、｜ｈ_i−ｈ_j｜の計算結果が「ゼロ」になるときに微小な値を１０^-20等とした場合、最終的な節点の収支式の誤差は所定の許容値以下となる。

しかしながら、初回計算時に｜ｈ_i−ｈ_j｜の計算結果が「ゼロ」になるときに使用する値として１０^-20を使用すると、桁落ちが発生し、Bi-CGSTAB法の計算が行えなくなる。

つまり、｜ｈ_i−ｈ_j｜の計算結果が「ゼロ」になるときに使用する微小な値によって節点方程式を解くことができないという問題が発生する。

本発明は、上記の事情に鑑みなされたもので、配水管網の節点エネルギー位を精度よく計算することを課題とする。

そこで、本発明の配水管網計算装置は、配水管網の水理解析を行う配水管網計算装置であって、配水管網の各節点における流量収支の絶対値の合計値が許容値を越えている場合に当該配水管網において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき当該各節点の節点エネルギー位の補正値を解とする方程式を作成する式作成手段と、前記隣り合う節点の節点エネルギー位の現在値の差がゼロである場合に前記方程式の変数である当該節点エネルギー位の差の累乗値を前記各節点の節点エネルギー位の補正回数に依る当該累乗値よりも極小な値に変更して当該方程式を補正する式補正手段と、前記補正された方程式の解として前記各節点の節点エネルギー位の補正値を算出する補正値演算手段と、前記算出された前記各節点の節点エネルギー位の補正値を当該各節点の節点エネルギー位の現在値に各々加算して当該現在値を補正するエネルギー位補正手段を備える。

また、本発明の配水管網計算方法は、配水管網の水理解析を行うコンピュータによる配水管網計算方法であって、配水管網の各節点における流量収支の絶対値の合計値が許容値を越えている場合に当該配水管網において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき当該各節点の節点エネルギー位の補正値を解とする方程式を作成する式作成ステップと、前記隣り合う節点の節点エネルギー位の現在値の差がゼロである場合に前記方程式の変数である当該節点エネルギー位の差の累乗値を前記各節点の節点エネルギー位の補正回数に依る当該累乗値よりも極小な値に変更して当該方程式を補正する式補正ステップと、前記補正された方程式の解として前記各節点の節点エネルギー位の補正値を算出する補正値演算ステップと、前記算出された前記各節点の節点エネルギー位の補正値を当該各節点の節点エネルギー位の現在値に各々加算して当該現在値を補正するエネルギー位補正ステップを有する。

尚、本発明は上記装置の各手段としてコンピュータを機能させるプログラムまたは上記方法の各ステップをコンピュータに実行させるプログラムの態様とすることもできる。

本発明によれば配水管網の節点エネルギー位を精度よく計算できる。

本発明の実施形態における配水管網計算装置のブロック構成図。同実施形態における節点、管路の接続関係の表示状態を説明した説明図。同実施形態における水理解析部のブロック構成図。同実施形態における水理解析の手順を説明したフローチャート。一つの配水池、複数の通常節点を有する管網モデルの一例。

以下に図面を参照しながら本発明の実施形態について説明する。

［本実施形態の概要］
図１に示された本実施形態の配水管網計算装置１は、図５に例示した配水管網１０の節点１２〜１７における流量収支の絶対値の合計値が許容値を越えている場合に配水管網１０において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき節点１２〜１７の節点エネルギー位の補正値を解とする方程式を作成する。次に、前記隣り合う節点の節点エネルギー位の現在値の差がゼロである場合に前記方程式の変数である当該節点エネルギー位の差の累乗値を節点１２〜１７の節点エネルギー位の補正回数に依る当該累乗値よりも極小な値に変更して当該方程式を補正する。次に、この補正した方程式の解として配水管網１０の節点１２〜１７における節点エネルギー位の補正値を算出する。そして、この節点１２〜１７の節点エネルギー位の補正値を節点１２〜１７の節点エネルギー位の現在値に各々加算して当該現在値を補正する。

［本実施形態の装置構成］
配水管網計算装置１は管網データ作成部２、水理解析部３、表示処理部４、出力表示部５を備える。これらの機能部２〜４はコンピュータのハードウェア資源（ＣＰＵ，ＲＡＭ，ＨＨＤ等）とソフトウエア資源（プログラム）との協働により実装される。

管網データ作成部２は出力表示部５（例えば、ディスプレイ）に出力表示された節点、管路の接続関係を設定する（図４のステップＳ１）。出力表示部５に表示された各節点の入力フォームには、「配水量」「需要量」等のデータの設定（入力）が行われる。一方、各管路の入力フォームには、「管路長」「管径」等のデータの設定（入力）が行われる。

水理解析部３は、例えば図５に示した配水管網１０の節点における節点エネルギー位の計算とその補正を行う。水理解析部３は、具体的な構成要素として図３に例示したように流量演算部３１、流量収支演算部３２、式作成部３３、式補正部３４、補正値演算部３５及びエネルギー位補正部３６を備える。

流量演算部３１は、配水管網１０の各管路における水の流量を計算する。

流量収支演算部３２は、配水管網１０の各節点の流量収支を計算してこの各流量収支の絶対値の合計値を算出する。

式作成部３３は、前記合計値が許容値を越えている場合に配水管網１０において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき配水管網１０の各節点における節点エネルギー位の補正値を解とする方程式を作成する。

式補正部３４は、前記隣り合う節点の節点エネルギー位の現在値の差がゼロである場合に前記方程式の変数である当該節点エネルギー位の差の累乗値を前記各節点の節点エネルギー位の補正回数に依る当該累乗値よりも極小な値に変更して当該方程式を補正する。

例えば、前記現在値の差がゼロであって前記補正回数が所定回数未満である場合、式補正部３４は前記累乗値をこの累乗値よりも指数関数的に極小な値に変更して前記方程式を補正する。

一方、前記現在値の差がゼロであって前記補正回数が所定回数以上である場合、式補正部３４は前記累乗値を前記指数関数的に極小な値よりもさらに指数関数的に極小な値に変更して前記方程式を補正する。

補正値演算部３５は、前記補正された方程式の解として配水管網１０における各節点の節点エネルギー位の補正値を算出する。

エネルギー位補正部３６は、前記各節点の節点エネルギー位の補正値を当該各節点の節点エネルギー位の現在値に各々加算して当該現在値を補正する。

表示処理部４は、前記水理解析の結果を出力表示部５に表示させる。例えば、前記水理解析結果を時系列データや各節点、管路毎のデータとしてグラフや数値等で表示させる。その他、例えば図２に示したように、出力表示部５に表示される節点５１〜５５や管路Ｒ１〜Ｒ５が前記水理解析の結果に基づき色分けされる。

［本実施形態の水理解析手順］
図４を参照しながら本実施形態の水理解析手順について説明する。

ここでは、図５に例示された１つの圧力既知の配水池（節点１１）と６つの通常節点（節点１２〜１７）を有する配水管網１０の水理解析の事例について説明する。

Ｓ１：配水管網１０の節点１１〜１７における節点エネルギー位の仮定値が管網データ作成部２によって設定される。初期値としてはゼロを設定してもよい。出力表示部５では表示処理部４によって前記仮定値の内容が表示される。具体的な設定例を以下に説明する。

節点１１〜１７間の管路１０２〜６０７の流量Ｑは流量公式（ヘーゼン・ウィリアムズ式）に基づき以下の式（１）で示される。Ｑ_ijは節点ｉ，節点ｊ間の管路における流量を意味する。

一方、節点方程式は以下の式（２）で示される。

式（１）においてＲ_ij＝０．２７８５３ＣＤ^2.63Ｌ^-0.54（Ｒ_ijは定数）とすると、管路１０２〜６０７の個々の流量Ｑ₁₂〜Ｑ₆₇は以下の式（３）〜（９）で示される。

図５に示された節点１２〜１７の個々の流量収支ｓ₂〜ｓ₇は式（２）に基づき以下の式（１０）〜（１５）で示される。

ここで、具体的な数値として、例えば以下の値を仮設定する。
Ｈ₁＝１００（圧力既知の配水池の節点エネルギー位）
Ｈ₂＝Ｈ₃＝Ｈ₄＝Ｈ₅＝Ｈ₆＝Ｈ₇＝０（各節点の節点エネルギー位の仮定値を０と設定）
Ｒ_ij＝５０（実際は管路長等により管路毎に異なるが簡略化のため全て同じ値を設定）
ｓ₃＝ｓ₄＝ｓ₅＝ｓ₆＝３０
Ｓ２：流量演算部３１は管路１０２〜６０７の水の流量を計算する。

ステップＳ１にて仮設定された値から式（３）〜（９）を計算すると図５に示された節点１１〜１７間の各流量Ｑ₁₂〜Ｑ₆₇は以下の式（１６）〜（２２）によって算出される。

Ｓ３：流量収支演算部３２は節点１２〜１６の流量収支を計算してこの各流量収支の絶対値の合計値を計算する。

具体的には式（１０）〜（１５）の計算により各節点の流量収支を計算する。

次いで上記の式（２３）〜（２８）で示される流量収支の絶対値の合計値を計算すると、以下の式（２９）のように算出される。

Ｓ４：ステップＳ３にて算出された流量収支の絶対値の合計値が所定の許容値ε以下であるかの判断を行い、当該合計値が許容値ε以下である場合には水理解析処理を終了（Ｅｎｄ）させる。一方、許容値εを超えている場合はステップＳ５を実行する。尚、許容値εは任意に設定される。

本事例では許容値ε＝０．０１に設定されている。ステップＳ３（式（２９）の計算）で算出された合計値は、７２１．１３となり、許容値ε＝０．０１を越えているので、次ぎのステップＳ５を実行する。

Ｓ５：式作成部３３は、前記合計値が許容値εを越えている場合に配水管網１０の各節点エネルギー位の補正値を解とする連立方程式（節点方程式）を作成する。

式（３）〜（９）は非線形式であるので、この連立方程式を解くには式を線形化する必要がある。線形化の方法としては以下の式（３０）による接線近似法（ニュートン法）が適用される。

ここで、Ｈ_i、Ｈ_jの仮定値をｈ_i、ｈ_jと設定し、ｑ_ij＝Ｒ_ij（ｈ_i−ｈ_j）と定義すると、式（３）は以下の式（３１）となる。

式（４）〜（９）についても同様の線形化を行うと各々以下の式（３３）〜（３８）となる。

ここで、Ｕ_ij＝０．５４Ｒ_ij｜ｈ_i−ｈ_j｜^-0.46と定義すると、式（３２）〜（３８）は各々以下の式（３９）〜（４５）となる。

そして、この式（３９）〜（４５）をそれぞれ流量収支の式（１０）〜（１５）に代入して整理すると以下の式（４６）〜（５１）となる。

Ｓ６：式補正部３４は、隣り合う節点の節点エネルギー位の現在値の差｜ｈ_i−ｈ_j｜が「０」でない場合、ステップＳ７を実行する。一方、前記差が「０」である場合、以下のステップＳ６０１〜Ｓ６０３で決定された極小な値をステップＳ５にて作成した節点方程式に代入してこの節点方程式を補正する。

ステップＳ６について具体的に説明すると、先ず、ステップＳ５にて得られた式（４６）〜（５１）において具体的な数値を代入して計算する。

式（５５）〜（６０）は、｜０−０｜^-0.46の値が無限大となるので、節点方程式（式（４６）〜式（５１））の解を算出できない。そこで、｜０｜^-0.46は所定の極小な値に置換される。

すなわち、繰り返し回数がＮ回未満である場合（Ｓ６０１）、前記節点方程式の変数「｜ｈ_i−ｈ_j｜^-0.46」がこの変数よりも指数関数的に極小な値「１０^-a」に置換される（Ｓ６０２）。一方、前記繰り返し回数がＮ回以上である場合（Ｓ６０１）、前記変数が「１０^-a」よりもさらに指数関数的に極小な値「１０^-b」に置換される（Ｓ６０３）。

「Ｎ」「ａ」「ｂ」の値は「０」以外の任意の自然数である（但し、ａ＜ｂ）。「ａ」の値は、例えば、初回のBi-CGSTAB法によって算出された計算値である。「ｂ」の値は、例えば、前記流量収支の絶対値の合計値が許容値ε以下となる値である。本事例ではＮ＝５、ａ＝５、ｂ＝２０に設定されている。

Ｓ７：補正値演算部３５は、ステップＳ６にて補正された連立方程式の解として、配水管網１０の各節点エネルギー位の補正値を算出する。具体的にはBi-CGSTAB法を用いて節点方程式を解き、節点エネルギー位の補正値を計算する。

先ず、式（５４）〜（６０）の計算を行う。但し、式（５５）〜（６０）はステップＳ６にて補正された式である。式（５４）〜（６０）の計算値はそれぞれ以下の式（６１）〜（６７）の値となる。

次いで、この式（６１）〜（６７）の値を式（４６）〜（５１）に各々代入すると以下の式（６８）〜（７３）の連立方程式が得られる。

この式（６８）〜（７３）の解として算出された以下の式（７４）〜（７９）の値が節点１２〜１７の節点エネルギー位Ｈ₂〜Ｈ₇の補正値ΔＨ₂〜ΔＨ₇となる。

Ｓ８：エネルギー位補正部３６は、ステップＳ７にて算出された節点１２〜１７の節点エネルギー位Ｈ₂〜Ｈ₇の補正値ΔＨ₂〜ΔＨ₇を節点１２〜１７の節点エネルギー位Ｈ₂〜Ｈ₇の現在値に各々加算して当該現在値を補正する。

本事例では補正前の節点１２〜１７の節点エネルギー位Ｈ₂〜Ｈ₇の現在値に式（７４）〜（７９）で示される補正値ΔＨ₂〜ΔＨ₇の値が各々加算される。補正された節点Ｈ₂〜Ｈ₇の節点エネルギー位は以下の式（８０）〜（８５）の値となる。

以上のステップＳ２〜Ｓ８からなる繰り返し回数１回分の処理（１Ｃ）が所定の許容値ε＝０．０１となるまで繰り返し実行される。

［本実施形態の効果］
以上のように配水管網計算装置１によれば、配水管網１０の各節点における流量収支の絶対値の合計値が許容値を超えて且つ隣り合う節点における節点エネルギー位の現在値の差がゼロである場合に当該各節点の節点エネルギー位の補正値を解とする方程式の変数である当該節点エネルギー位の差の累乗値が節点エネルギー位の補正回数に依る当該累乗値よりも極小な値に置換される。そして、この置換によって補正された方程式の解として算出された節点エネルギー位の補正値がこの節点エネルギー位の現在値に加算されることにより当該現在値が補正される。これにより、配水管網１０の節点が多くなった場合の水理計算にあたり、桁落ちを起こさせることなく配水管網１０の節点エネルギー位を算出できる。したがって、配水管網１０の節点エネルギー位を精度よく計算できる。

特に、行列の非ゼロ要素のみを扱うBi-CGSTAB法を用いた水理解析により計算速度を向上させることができる。また、コンピュータの計算資源を無駄なく使用することができる。さらに、従来の配水管網計算法よりも多くの節点数の水理解析を行うことができる。

また、前記補正の繰り返し回数が所定回数未満である場合に、前記累乗値よりも指数関数的に極小な値が前記流量収支の節点方程式に供されて算出された節点エネルギー位の補正値が前記節点エネルギー位の値に加算されて当該節点エネルギー位が補正される。一方、前記繰り返し回数が所定回数以上である場合に、前記指数関数的に極小な値よりもさらに指数関数的に極小な値が前記節点方程式に供されて算出された節点エネルギー位の補正値が前記節点エネルギー位の値に加算されて当該節点エネルギー位が補正される。このように補正の繰り返し回数に依り節点方程式の変数が任意に設定されることにより、前記繰り返し数に応じて節点エネルギー位の補正を任意に行える。そして、前記各節点における流量収支の絶対値の合計値が前記許容値以下となるまで前記補正が繰り返し実行されることにより、前記節点エネルギー位の精度がさらに向上する。

［本発明の他の態様］
尚、本発明はプログラムの態様とすることもできる。例えば、コンピュータを配水管計算装置１の機能部３３〜３６として機能させるプログラムや上記のステップＳ２〜Ｓ８の一部若しくは全てをコンピュータに実行させるプログラムが挙げられる。これらのプログラムはコンピュータが読み取り可能な周知の記録媒体（例えば、ＣＤ、ＤＶＤ等）またはネットワーク（例えば、インターネット、メール等）を介して提供できる。

１…配水管網計算装置
２…管網データ作成部
３…水理解析部、３３…式作成部（式作成手段）、３４…式補正部（式補正手段）、３５…補正値演算部（補正値演算手段）、３６…エネルギー位補正部（エネルギー位補正手段）
４…表示処理部

Claims

配水管網の水理解析を行う配水管網計算装置であって、
配水管網の各節点における流量収支の絶対値の合計値が許容値を越えている場合に当該配水管網において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき当該各節点の節点エネルギー位の補正値を解とする方程式を作成する式作成手段と、
前記隣り合う節点の節点エネルギー位の現在値の差がゼロであって当該各節点の節点エネルギー位の補正回数が所定回数未満である場合に前記方程式の変数である当該節点エネルギー位の差の累乗値（但し、冪指数は負）を１０^-a（但し、ａは０以外の自然数）の値に変更して当該方程式を補正する一方で当該現在値の差がゼロであって当該補正回数が所定回数以上である場合に当該累乗値を１０ ^-b （但し、ｂは０以外の自然数、ａ＜ｂ）の値に変更して当該方程式を補正する式補正手段と、
前記補正された方程式の解として前記各節点の節点エネルギー位の補正値を算出する補正値演算手段と、
前記算出された前記各節点の節点エネルギー位の補正値を当該各節点の節点エネルギー位の現在値に各々加算して当該現在値を補正するエネルギー位補正手段と
を備えたことを特徴とする配水管網計算装置。
配水管網の水理解析を行うコンピュータによる配水管網計算方法であって、
配水管網の各節点における流量収支の絶対値の合計値が許容値を越えている場合に当該配水管網において隣り合う節点の節点エネルギー位の差からこの隣り合う節点間の流量を算出する公式に基づき当該各節点の節点エネルギー位の補正値を解とする方程式を作成する式作成ステップと、
前記隣り合う節点の節点エネルギー位の現在値の差がゼロであって当該各節点の節点エネルギー位の補正回数が所定回数未満である場合に前記方程式の変数である当該節点エネルギー位の差の累乗値（但し、冪指数は負）を１０^-a（但し、ａは０以外の自然数）の値に変更して当該方程式を補正する一方で当該現在値の差がゼロであって当該補正回数が所定回数以上である場合に当該累乗値を１０ ^-b （但し、ｂは０以外の自然数、ａ＜ｂ）の値に変更して当該方程式を補正する式補正ステップと、
前記補正された方程式の解として前記各節点の節点エネルギー位の補正値を算出する補正値演算ステップと、
前記算出された前記各節点の節点エネルギー位の補正値を当該各節点の節点エネルギー位の現在値に各々加算して当該現在値を補正するエネルギー位補正ステップと
を有することを特徴とする配水管網計算方法。