"複数のポンプの運転を制御するポンプ制御システムであって、配水流量と吐出圧とに応じた運転すべきポンプの台数を定めた運転台数判定テーブルと、前記ポンプの台数と前記配水流量と前記吐出圧とに応じた各ポンプへの前記配水流量の配分を定めた流量配分テーブルと、を記憶する記憶部と、前記運転台数判定テーブルを参照し、前記配水流量と前記吐出圧とに基づいて、運転すべきポンプの台数を決定する運転台数判定部と、前記運転すべきポンプの台数と、前記流量配分テーブルにより定められた前記配水流量と前記吐出圧とに基づいて、前記各ポンプへの配水流量の配分を計算する配水流量計算部と、前記運転台数判定テーブルに記憶された前記運転すべきポンプの台数と、前記配水流量計算部が求めた前記各ポンプへの配水流量の配分とに基づいて、前記各ポンプの目標回転数を設定する目標回転数計算部と、を備えたことを特徴とするポンプ制御システム。請求項1に記載のポンプ制御システムであって、前記運転台数判定テーブルにより定められる前記ポンプの台数、および前記流量配分テーブルにより定められる前記配水流量は、前記各ポンプの特性を示す特性曲線と、あらかじめ定められた運転順序によって運転される各ポンプの全体の消費エネルギーが最小となるような前記ポンプの運転台数、および前記各ポンプに配分される配水流量を前記配水流量および前記吐出圧ごとに求めたものである、ことを特徴とするポンプ制御システム。請求項1または2に記載のポンプ制御システムであって、前記ポンプの配水流量と前記吐出圧と前記各ポンプの回転数と前記各ポンプの消費電力とを計測し、計測されたこれらの情報に基づいて、前記各ポンプの特性を推定するポンプ特性推定部を備える、ことを特徴とするポンプ制御システム。請求項3に記載のポンプ制御システムであって、前記ポンプ特性推定部は、計測した前記ポンプの配水流量と前記吐出圧と前記各ポンプの回転数と前記各ポンプの消費電力とに基づいて定められる前記各ポンプの特性を示す特性曲線および性能曲線を定めるためのパラメータを本来の設定値とは異なる設定値に変更して前記各ポンプを運転したときの配水流量と吐出圧と各ポンプの回転数と消費電力とを計測し、計測したこれらの情報に基づいて、前記各ポンプの特性を推定する、ことを特徴とするポンプ制御システム。請求項2のポンプ制御システムにおいて、ポンプの特性曲線として、請求項3記載の推定方法で推定されるポンプの特性を利用することを特徴とするポンプ制御システム請求項2のポンプ制御システムにおいて、ポンプの特性曲線として、請求項4記載の推定方法で推定されるポンプの特性を利用することを特徴とするポンプ制御システム請求項4~6のいずれか1項に記載のポンプ制御システムであって、前記ポンプ特性推定部が推定する各ポンプの特性とは、ポンプの劣化特性を示すものである、ことを特徴とするポンプ制御システム。複数のポンプの運転を制御するポンプ制御システムであって、前記複数のポンプの配水流量を計測するための流量センサと、前記複数のポンプの吐出圧を計測するための圧力センサと、前記配水流量と前記吐出圧との組み合わせに対して、最小消費電力を与えるポンプの運転台数と流量配分とを決定し、前記各ポンプの特性曲線および性能曲線に基づいて、前記流量センサが計測した配水流量と前記圧力センサが計測した吐出圧とに応じた運転すべきポンプの台数を定めた運転台数判定テーブルと、前記ポンプの台数と前記配水流量と前記吐出圧とに応じた各ポンプへの前記配水流量の配分を定めた流量配分テーブルと、を作成する最適化部と、前記運転台数判定テーブルと前記流量配分テーブルとを記憶する記憶部と、前記運転台数判定テーブルを参照し、前記配水流量と前記吐出圧とに基づいて、運転すべきポンプの台数を決定する運転台数判定部と、前記運転すべきポンプの台数と、前記流量配分テーブルにより定められた前記配水流量と前記吐出圧とに基づいて、前記各ポンプへの配水流量の配分を計算する配水流量計算部と、前記運転すべきポンプの台数と、前記配水流量計算部が求めた前記各ポンプへの配水流量の配分とに基づいて、前記各ポンプの目標回転数を設定する目標回転数計算部と、を備えたことを特徴とするポンプ制御システム。請求項8のポンプ制御システムにおいて、前記各ポンプに対して定められた運転順序のパターンのそれぞれについて、前記配水流量と前記吐出圧の1日における実績時系列値に対して、任意の時刻におけるポンプの最適の運転台数と流量配分とを計算してポンプの消費電力を算出し、算出した前記消費電力をすべての時刻について加算することにより前記ポンプ制御システム全体の1日の消費電力を算出し、算出した前記1日の消費電力が最小となる運転順序を、最適なポンプの運転順序として採用する消費エネルギー計算部、を備えたことを特徴とするポンプ制御システム。請求項9のポンプ制御システムにおいて、前記各ポンプの回転数を計測する回転数センサと、前記流量センサが計測した前記配水流量と前記圧力センサが計測した前記吐出圧との組み合わせと前記各ポンプの回転数と前記消費電力とを所定の期間蓄積したデータに基づいて、前記各ポンプの性能曲線および効率曲線を推定するポンプ特性推定部と、を備えたことを特徴とするポンプ制御システム。上述した特許文献1は、ポンプの効率特性曲線を考慮することで省エネ運転を実現するポンプ運転台数の増減の方法が記載されているが、特性の異なる複数のポンプを活用する場合、運転台数だけでなく、各ポンプが受け持つ配水流量も同時に決定し運転しなければならない。従来技術はこの点の配慮がなく、更なる省エネ運転の余地があるという問題がある。また、ポンプ特性の経時劣化について配慮がなされていないため、時間ともに効率点での運転ができなくなり必要以上に電力を消費するという問題がある。また、上述した特許文献2では、ポンプのトータル軸駆動力を考慮することで最適な台数切り替えのタイミング流量を決定しているが、技術の適用が吐出圧一定制御、あるいは、推定末端圧一定制御に限定されており、より高度な制御には適用できないという問題がある。また、ポンプのトータル消費エネルギーと必ずしも一致しないトータル軸駆動力に基づいて台数切り替え流量解列流量を決定しているため厳密な消費エネルギー最小化は実現できないという問題がある。また、台数切り替えの運転状態は、流量のみならず圧力にも依存するが流量のみで台数切り替えタイミング解列タイミングを決定しており、省エネを実現する厳密な切り替えが行えないという問題がある。また、ポンプ特性の経時劣化について配慮がなされていないため、時間ともに効率点での運転ができなくなり必要以上に電力を消費するという問題がある。本発明は、こうした事情を考慮してなされたもので、省エネルギーを実現する運転台数だけでなく各ポンプが受け持つ最適な配水流量も同時に求めることで消費エネルギーの最小化を実現するポンプ制御システムを提供することにある。また、オンラインでポンプの特性曲線を推定し、それをポンプの効率運転に利用することで、時間によらずポンプの最適効率点運転を維持可能なポンプ制御システムを提供することにある。本発明第1の発明のポンプ制御システムは、複数のポンプの運転を制御するポンプ制御システムであって、配水流量と吐出圧とに応じた運転すべきポンプの台数を定めた運転台数判定テーブルと、前記ポンプの台数と前記配水流量と前記吐出圧とに応じた各ポンプへの前記配水流量の配分を定めた流量配分テーブルと、を記憶する記憶部と、前記運転台数判定テーブルを参照し、前記配水流量と前記吐出圧とに基づいて、運転すべきポンプの台数を決定する運転台数判定部と、前記運転すべきポンプの台数と、前記流量配分テーブルにより定められた前記配水流量と前記吐出圧とに基づいて、前記各ポンプへの配水流量の配分を計算する配水流量計算部と、前記運転台数判定テーブルに記憶された前記運転すべきポンプの台数と、前記配水流量計算部が求めた前記各ポンプへの配水流量の配分とに基づいて、前記各ポンプの目標回転数を設定する目標回転数計算部と、を備えたことを特徴とする。本発明第2の発明のポンプ制御システムは、第1の発明において、前記運転台数判定テーブルにより定められる前記ポンプの台数、および前記流量配分テーブルにより定められる前記配水流量は、前記各ポンプの特性を示す特性曲線と、あらかじめ定められた運転順序によって運転される各ポンプの全体の消費エネルギーが最小となるような前記ポンプの運転台数、および前記各ポンプに配分される配水流量を前記配水流量および前記吐出圧ごとに求めたものである、ことを特徴とする。本発明第3の発明のポンプ制御システムは、前記ポンプの配水流量と前記吐出圧と前記各ポンプの回転数と前記各ポンプの消費電力とを計測し、計測されたこれらの情報に基づいて、前記各ポンプの特性を推定するポンプ特性推定部を備える、ことを特徴とする。本発明第4の発明のポンプ制御システムは、前記ポンプ特性推定部は、計測した前記ポンプの配水流量と前記吐出圧と前記各ポンプの回転数と前記各ポンプの消費電力とに基づいて定められる前記各ポンプの特性を示す特性曲線および性能曲線を定めるためのパラメータを本来の設定値とは異なる設定値に変更して前記各ポンプを運転したときの配水流量と吐出圧と各ポンプの回転数と消費電力とを計測し、計測したこれらの情報に基づいて、前記各ポンプの特性を推定する、ことを特徴とする。本発明第5の発明のポンプ制御システムは、第2の発明においてポンプの特性曲線として、第3の発明で推定されるポンプの特性を利用することを特徴とする。本発明第6の発明のポンプ制御システムは、第2の発明においてポンプの特性曲線として、第4の発明で推定されるポンプの特性を利用することを特徴とする。本発明では、ポンプ運転を規定する配水流量と吐出圧の2つの変数をベースにしてポンプのトータル消費エネルギー最小となるポンプ運転台数とポンプの流量配分を決定するのでポンプ制御システムの省エネルギー化が実現できる。また、ポンプの特性を、計測データから推定し、最新の推定特性に基づいてポンプ運転台数、ポンプ流量配分を決定するのでポンプ特性の経時変化によらず最適効率点での運転を維持できる。また、本来のパラメータを異なるパラメータでポンプを運転することで、オンラインにおいて、より幅広い運転領域での運転が実現できる。そこで計測したデータを活用することでポンプ特性を高精度に推定できるので、最適効率点運転の信頼性を向上でき、電力消費の増加を極限まで抑制できる。本発明の実施形態を、図面を用いて説明する。図1から図7を参照し、実施例1を説明する。図1は、実施例1のポンプ制御システムの構成図である。制御システムは、配水管1、配水池11、吐出圧を計測する圧力センサ2、配水流量を計測する流量センサ3、可変速ポンプ4、5、6、ポンプ回転数を計測する回転数センサ7、8、9、目標吐出圧、計測吐出圧、計測配水流量を入力として目標吐出圧を実現するような各ポンプの目標回転数を計算するポンプ制御システム100、計測回転数が目標回転数に一致するよう各ポンプへの入力信号を調整するPID制御装置12、13、14を有して構成されている。図1では、ポンプを並列に3台配置した例を示しているが、これ以外の台数を設置することもできる。配水流量増加に対してポンプの運転台数は増加する。ここでは、流量が小さい少ないときはポンプ4のみを運転し、流量が増加するにしたがって、ポンプ5、ポンプ6の順にポンプを追加し、運転台数を増加させるようになっている。流量センサ3により計測される配水流量は、水需要に応じて時々刻々変化する。ポンプ制御システム100は、この流量変動に対して、圧力センサ2による吐出圧計測値が目標吐出圧に一致するよう各ポンプの目標回転数回転数0の場合はそのポンプは停止とみなすを算出する。目標吐出圧を実現する回転数の組み合わせは無数に存在するが、ポンプ制御システムは、最小消費エネルギーを与える回転数組み合わせを算出する。これを実現するため、ポンプ制御システムは、運転台数判定手段101、運転台数判定テーブル102、流量配分計算手段103、流量配分テーブル104、目標回転数計算手段105を備えている。なお、ポンプ制御システムが有する上述した各手段は、実際にはそれらの各機能を実現するためのソフトウェアプログラムにより実行される。これらのソフトウェアプログラムは、例えば、上述した各部を含むモジュール構成となっており、実際のハードウェアとしては、CPU等の制御部がHDD等の記録装置からこれらのソフトウェアプログラム、あるいは上述した各テーブルを読み出して実行することにより、上記各部が主記憶装置上にロードされ、運転台数判定手段101、流量配分計算手段103、目標回転数計算手段105の各部が主記憶装置上に生成されるようになっている。運転台数判定手段101では、予め構築された運転台数判定テーブル102を参照し、運転すべきポンプの台数を決定する。運転台数判定テーブル102の一例を図2に示す。横軸、縦軸はそれぞれ、吐出圧H、配水流量Qである。流量小流量が少ない場合がポンプ1台運転領域ポンプ4のみ運転、流量中ぐらい流量が中程度の場合が2台運転領域ポンプ4、5の2台運転、流量大流量が多い場合が3台運転領域になる。各領域の境界は、実践と点線の2つが存在する。実線は、配水流量が増加してゆき運転台数が増えるときに活用する台数切り替えラインであり、点線は配水流量が減少して運転台数が減るとき活用する台数切り替えラインである。このように流量の増減で切り替えラインを変えるのは、運転が境界上に留まったときポンプ運転のオンオフが頻繁に起こるのを防ぐためのものである。運転台数判定手段101では、最新の配水流量Qと吐出圧Hの計測値を運転台数判定テーブル上にプロットしたとき、そのプロットがどの台数領域にあるかにより運転台数を決定する。もし、実線切り替え線と点線の切り替え線の間の狭い領域にプロットが存在した場合は、1時刻前の運転台数を参照し、その運転台数を最新時刻の運転台数とする。このテーブルは、ポンプのトータル消費エネルギーが最小となる運転台数を与えるものであるが、テーブル作成方法については後述する。流量配分計算手段103では、ポンプを複数台運転する場合の各ポンプが受け持つ配水流量を、流量配分テーブル104を参照して算出する。図3に流量配分テーブル104の一例を示す。上段aがポンプ2台運転時に1台のポンプポンプ4の流量配分D1を求めるテーブル31であり、下段bがポンプ3台運転時に2台のポンプポンプ4、5の流量配分D1、D2を求めるテーブル32、33である。ポンプ運転台数は、運転台数判定手段101で算出された情報を活用する。ポンプ2台運転の場合は、配水流量計測値と吐出圧計測値をベースにテーブル31を参照し、ポンプ4への流量配分D1を求める。配水流量計測値をQとしたとき、ポンプ4、ポンプ5の受け持つ流量Q1、Q2は次式で与えられる。Q1=QD1Q2=QH=a2N2-b2Q2H=a3N2-b3Q2ここで、Hはポンプ揚程吐出圧、Nはポンプ回転数、Qはポンプ配水流量、ai,biは正の定数である。流量配分計算手段103で求められたポンプ4,5,6の流量をQ1、Q2、Q3、ポンプ4,5,6の回転数をN1、N2、N3、吐出圧計測値をH0とし、上記性能曲線に代入すると次式が得られる。H0=a1N12-b1Q12H0=a2N22-b2Q22H0=a3N32-b3Q32これらの式を回転数N1、N2、N3について解いて次式を得る。N1=H0b1Q12a112N2=H0b1Q22a212N3=H0b1Q32a312もし近似した性能曲線にモデル化誤差がなければ、この計算値を各ポンプの目標回転数として与えれば、目標の吐出圧を達成できる。すなわち、計測吐出圧が目標圧に等しくなる。しかしながら、モデル化誤差は少なからず存在するので、目標の吐出圧は達成されず、計測吐出圧と目標吐出圧の間にわずかな差が生じる。モデル化誤差に起因する目標値と計測値のずれを補償するためフィードバック制御を併用する。これは、目標吐出圧と計測吐出圧の差に基づき、例えばPID制御則により回転数補正量ΔNを計算するものである。上の式で計算されるN1、N2、N3にΔNを加算し、ポンプ4,5,6の目標回転数N10、N20、N30を次のように算出する。N10=N1ΔNN20=N2ΔNN30=N3ΔNこの回転数計算値を、ポンプの速度制御を行うPID制御装置に送り、ポンプの回転数制御を行う。なお、2台運転の場合は式78を、1台運転の場合は式7を活用して目標回転数を計算すれば良い。次に、図2と図3のテーブルデータを作成するための方法を、図4、図5に基づいて説明する。図4は、テーブルデータを作成するための処理手順を示すものである。様々な配水流量Qと吐出圧Hの組み合わせに対して、次に述べる最適化問題を解いて、その組み合わせに対する最小エネルギーを与える運転台数と流量配分を決定する。以下に、ポンプを1台運転する場合、2台同時運転する場合、3台同時に運転する場合のポンプの消費電力を計算する方法を述べる[1]1台運転の場合次式に基づいて、ポンプ1台運転時ポンプ4運転時の消費電力Pを決定する。P=kQH10ηp=-c1QN2d1QNH=a1N2-b1Q2ここに、Pはポンプ4の消費電力kW、Qは配水流量m3min、Hはポンプ吐出圧m、ηpはポンプ効率、ηmはモータ効率で定数、ηiはインバータ効率で定数、Nはポンプ回転数rpm、k,a1,b1,c1,d1は正の定数である。式はポンプ4の効率曲線を近似したもの、式はポンプ4の性能曲線を近似したものである。Q,Hが与えられると、12式から回転数Nが求められる。それを1110式に代入することで、ポンプの消費電力Pが求められる。[2]2台同時運転の場合次式に基づいて、ポンプ2台運転時ポンプ4,5運転時の消費電力Pが計算できる。P=P1P2P1=kQ1H14ηp1=-c1Q1N12d1Q1N1H=a1N12-b1Q12P2=kQ2H17ηp2=-c2Q2N22d2Q2N2H=a2N22-b2Q22Q=Q1Q2ここに、Pはポンプのトータル消費電力kW、P1はポンプ4の消費電力kW、P2はポンプ5の消費電力kW、Qは配水流量m3min、Q1はポンプ4の配水流量m3min、Q2はポンプ5の配水流量m3min、Hはポンプ吐出圧m、ηp1はポンプ4のポンプ効率、ηp2はポンプ5のポンプ効率、ηmはモータ効率で定数、ηiはインバータ効率で定数、N1はポンプ4の回転数rpm、N2はポンプ5の回転数rpm、k,ai,bi,ci,diは正の定数である。配水流量Qを満足するQ1,Q2の組み合わせは無数に存在し、それに対応してトータル消費電力Pが決まる。ここでは、Pを最小にするQ1,Q2を最適化計算で求めて、これに対応する消費電力P2台運転時の最小の消費電力を2台運転時のポンプの消費電力とする。[3]3台同時運転の場合次式に基づいて、ポンプ3台運転時ポンプ4,5、6運転時の消費電力Pが計算できる。P=P1P2P3P1=kQ1H22ηp1=-c1Q1N12d1Q1N1H=a1N12-b1Q12P2=kQ2H25ηp2=-c2Q2N22d2Q2N2H=a2N22-b2Q22P3=kQ3H28ηp3=-c3Q3N32d2Q3N3H=a3N32-b3Q32Q=Q1Q2Q3ここに、Pはポンプのトータル消費電力kW、P1はポンプ4の消費電力kW、P2はポンプ5の消費電力kW、P3はポンプ6の消費電力kW、Qは配水流量m3min、Q1はポンプ4の配水流量m3min、Q2はポンプ5の配水流量m3min、Q3はポンプ6の配水流量m3min、Hはポンプ吐出圧m、ηp1はポンプ4のポンプ効率、ηp2はポンプ5のポンプ効率、ηp3はポンプ6のポンプ効率、ηmはモータ効率で定数、ηiはインバータ効率で定数、N1はポンプ4の回転数rpm、N2はポンプ5の回転数rpm、N3はポンプ6の回転数rpm、k,ai,bi,ci,diは正の定数である。配水流量Qを満足するQ1,Q2,Q3の組み合わせは無数に存在し、それに対応してトータル消費電力Pが決まる。ここでは、Pを最小にするQ1,Q2,Q3を最適化計算で求めて、これに対応する消費電力P3台運転時の最小の消費電力を3台運転時のポンプの消費電力とする。以上の方法を活用すれば、配水流量Qと吐出圧Hに対応する1台運転時のポンプの消費電力、2台運転時のポンプのトータル消費電力最小消費電力、3台運転時のポンプのトータル消費電力最小消費電力を求めることができる。この3つの消費電力で最小値を与えるものが最も効率の高いポンプの運転方法運転台数と流量配分になる。図4の最適解DB作成部42では、以上の方法で、様々なQ,Hに対して最も効率的な運転方法運転台数と流量配分をオフラインで求めておき、図5に示すテーブルに記憶しておく。このテーブルは、配水流量Qと吐出圧Hの様々な値図では、Qは5.0から30.0まで0.1刻みの値、Hは25.0から40.0まで0.1刻みの値に対して、最適消費電力最小となるポンプの運転台数と流量配分を記憶したものである。このテーブルデータを活用して、図2の運転台数判定テーブル、図3の流量配分テーブルを作成することができる。なお、N台からN1台あるいはN1台からN台N=1,2の最適切り替えラインは1本になるが、その前後に2本のラインを設定し、上方にある線をN台からN1台への切り替えライン、下方にある線をN1台からN台への切り替えラインとするようにして図2の台数判定テーブルを作成する。以上述べた方法は、ポンプ運転が、流量増加に対してポンプ4→ポンプ4、5の2台同時運転→ポンプ4、5、6の3台同時運転と台数変化する場合の最適運転台数、流量配分を求めたものである。後述するポンプの特性変化によっては、この運転順序ポンプ4→ポンプ4、5→ポンプ4、5、6が最も効率的とは言えなくなる場合がある。すなわち、ポンプ6が単体として最も効率が良くなった場合、流量が小さい領域でポンプ6を運転するほうが一日の電力消費を少なくできる場合がある。様々な運転順序ポンプ3台を利用する場合は、6通りの運転順序があるを想定して運転台数判定テーブル、流量配分テーブルを求めておき、一日の運転を行ったときの最小消費電力を与える運転順序を決定しポンプ制御に活用する必要がある。運転順序は以下の6通りである[1]ポンプ4→ポンプ4,5→ポンプ4、5、6[2]ポンプ4→ポンプ4,6→ポンプ4、5、6[3]ポンプ5→ポンプ4,5→ポンプ4、5、6[4]ポンプ5→ポンプ5、6→ポンプ4、5、6[5]ポンプ6→ポンプ4,6→ポンプ4、5、6[6]ポンプ6→ポンプ5、6→ポンプ4、5、6図6は、あるポンプの運転順序に対して1日の消費エネルギーを算出するための方法の説明図である。まず、[1]の運転順序に対して、運転台数判定テーブル、流量配分計算テーブルを作成し、運転台数、流量配分計算手段61に記憶しておく。一日の配水流量と吐出圧の実績時系列Q、Htは時刻に対して、時刻tにおける最適の運転台数X、流量配分Dを計算し、それを利用して消費エネルギー計算手段は、ポンプのトータル消費電力Ptを計算し、さらに、すべての時刻でPを加算することで1日の消費エネルギーを算出する。これを[2]から[6]の運転順序に対して行い、最小の消費エネルギーを与える運転順序を最適な運転順序として採用する。本来の順序[1]と異なる順序が得られたなら、その順序で図1に示すポンプ運転制御を行うようにする。なお、ポンプ消費エネルギー計算は、式10から31を活用して行うことができる。図7は、第1実施例のポンプ制御システムにおいて、目標回転数を決定する制御プログラムのフローを示すもので、所定の周期例えば1分周期で実行されるようになっている。ステップ701では、計測配水流量、計測吐出圧、図2のポンプ運転台数判定テーブルを活用してポンプ運転台数を決定する。ステップ702では、ポンプの運転台数が2台以上の場合のみ、計測配水流量、計測吐出圧、図3のポンプ流量配分計算テーブルを活用してポンプ流量配分を決定し、さらに、それに基づいて各ポンプの受け持つ配水流量を計算する。最後にステップ703では、前述の方法、すなわち式789を活用して各ポンプの目標回転数を計算し処理を終了する。この計算値は、ポンプの回転数を制御するPID制御装置に送られ、ポンプの回転数が制御される。以上で実施例1の説明を終了する。次に、図8から図11に基づいて、様々な計測データからポンプ特性を推定して運転台数判定テーブルと流量配分テーブルを作成する第2実施例について説明する。図8は、ポンプの性能曲線と効率曲線を示すものである。実線はポンプが新しい状態の特性、点線はポンプが経時劣化した場合の特性の一例を示している。図でNは、規格化された回転数実回転数を定格回転数で割ったものを表している。図で示すような特性劣化が発生した場合、最適効率点での運転が維持できなくなり、余分に電力消費するという問題が発生する。最適効率点での運転を維持するためには、劣化した特性を推定してポンプ運転制御に反映する、すなわち、劣化特性に基づいて運転台数判定テーブル、流量配分テーブルを作成して、ポンプ運転制御を行う必要がある。テーブル作成の方法を以下で説明する。図9は、図1のポンプ制御システムでポンプの運転を行った場合の、様々な時刻に同時に計測された配水流量計測値Qと吐出圧計値Hの組み合わせQ,Hを、ポンプの特性平面上にプロットしたものである。これらの計測データを活用して、ポンプの性能曲線と効率曲線を推定できる。図10は、ポンプ特性データを推定して流量配分テーブル、運転台数判定テーブルを作成する第2実施例の全体構成図である。データベースDB106には、例えば、10分おきに、圧力センサにより計測される吐出圧H、流量センサにより計測される配水流量Q、回転数センサにより計測される各ポンプの回転数N、各ポンプの消費電力の計測値が入力、記憶されるようになっている。記憶データの一例を、図11に示す。計測時点におけるポンプの運転台数分だけ回転数と消費電力が記憶されるようになっている。なお、DBには過去1週間分のデータが格納されるようになっている。ポンプ特性推定手段107では、以下の方法でDBに格納されたデータから各ポンプの性能曲線と効率曲線を推定する。ここでは、運転順序がポンプ4→ポンプ4,5→ポンプ4、5、6であるとする。ポンプ1台運転時ポンプ4運転時のデータをDBから抽出し、それを配水流量Qi、吐出圧Hi、回転数Ni、消費電力Pii=1,2...nとする。ポンプ4の性能曲線、特性曲線が次式で近似できるとする。H=a1N2-b1Q232ηp=-c1QN2d1QNまた、消費電力に関して次式が成立する。P=kQHここに、Pはポンプ消費電力kW、Qは配水流量m3min、Hはポンプ吐出圧m、ηpはポンプ効率、ηmはモータ効率で定数、ηiはインバータ効率で定数、Nはポンプ回転数rpm、k,a1,b1,c1,d1は正の定数である。式からηpを消去して次式を得る。kQH=-c1QN2d1QNポンプ特性を推定することは、式3233内のパラメータa1,b1,c1,d1を推定することである。パラメータa1,b1,c1,d1は、n個の計測データQi、Hi、Ni、Pii=1,2...n活用し、これらのデータが式3235に最もフィッティングするよう最小二乗法を活用して行うことができる。得られたパラメータと式をもって、ポンプ4の最新の特性とする。次に、ポンプ2台同時運転時のデータをDBから抽出し、それを配水流量Qi、吐出圧Hi、ポンプ4の回転数N1i、ポンプ5の回転数N2i、ポンプ5の消費電力P2ii=1,2...mとする。これらのデータを活用してポンプ5の特性を推定する。ポンプ5の性能曲線、効率曲線が次式で近似できるとする。H=a2N22-b2Q2236ηp2=-c2Q2N22d2Q2N2また、消費電力に関して次の式が成立する。P2=kQ2Hここに、P2はポンプ5の消費電力kW、Q2はポンプ5の配水流量m3min、Hはポンプ吐出圧m、ηp2はポンプ5の効率、ηmはモータ効率で定数、ηiはインバータ効率で定数、N2はポンプ5の回転数rpm、k,a2,b2,c2,d2は正の定数である。式からηp2を消去して次式を得る。kQ2H=-c2Q2N22d2Q2N2計測データである吐出圧Hi、ポンプ4の回転数N1iを式32に代入して次式を得る。Hi=a1N1i2-b1Q2この式からa1,b1は推定済みであるので、ポンプ4の流量Qが次式で求められる。Q=12このQと計測されたトータル配水流量Qiから、ポンプ5の配水流量Q2iを次式で計算できる。Q2i=Qi-QこのQ2iと計測データである吐出圧Hi、ポンプ4の回転数N1i、ポンプ5の回転数N2i、ポンプ5の消費電力P2ii=1,2...mのmセットのデータが、式に最もフィットするよう最小二乗法を活用してパラメータa2,b2,c2,d2を推定する。得られたパラメータと式をもって、ポンプ5の最新の特性とする。同様の考え方を用いて、3台同時運転時の計測データを活用してポンプ6の性能曲線と効率特性を推定できる。ここでは、ポンプ4,5の性能特性曲線が既知として、それらの計測データを活用してポンプ4,5の配水流量を計算する。計測されるトータル配水流量からポンプ4,5の配水流量を減算してポンプ6の配水流量が計算できる。これを利用して前述と同様の方法で、ポンプ6の特性を推定できる。図12のテーブル作成手段108では、得られた最新のポンプ特性を活用して図4の方法で、運転台数判定テーブルと流量配分テーブルを作成する。作成されたテーブルは、ポンプ制御システム100にセットされポンプ制御が行われる。以上、最新のポンプ特性データを活用してポンプ制御を実施することで、ポンプの最適効率点での運転を維持でき、特性劣化に伴う電力消費の増加をある程度抑制できる。上記処理は、オンラインでもオフラインでも行うことができる。次に図12に基づいて第3実施例を説明する。図12はより広い運転領域でのデータ計測を可能にしたポンプ制御システムの全体構成図である。図9に示すようにポンプの最適運転では、計測が行われる運転領域が限られるため、それらのデータを活用してもポンプの性能曲線、効率曲線を正しく推定するのは困難な場合があると考えられる。この問題に対応し、より広い運転領域でのデータ計測を可能にするため目標回転数計算手段の処理内容を変更している。図12の目標回転数計測手段1051では、式から9に代わり次式で目標回転数を計算する。N10=N1×補正係数N20=N2×補正駅数N30=N3×補正係数補正係数は、1より大きな定数で、例えば、1.1から1.3程度の値に設定する。本来より大きな回転数でポンプを運転するものである。Nが大きくなるので、図9のbでQNのより小さい領域での運転や、図9のaでHのより大きな領域での運転が可能になり、計測できる運転領域を広げることができる。例えば、1日だけこのような補正処理を行い、データ計測を行う。このデータに従来の計測データを加えた、より広い領域のデータを利用してポンプ特性推定を行えばポンプ特性の推定精度を向上できる。なお、本補正で目標より少し高い吐出圧が実現されるのでエネルギー消費は一時的に増大することになる。次に図13と図14に基づいて第4実施例を説明する。図13はより広い運転領域でのデータ計測を可能にしたポンプ制御システムの全体構成図であり、第3実施例とは別の実施形態を与えるものである。具体的には、従来と異なる修正運転台数判定テーブル1021を用いる。図14はこのテーブルの一例を示すものである。1台運転から2台運転への切り替えラインを本来のものより上方にシフトした修正切り替えラインとすることで、1台運転ポンプの運転領域を拡大するものである。運転領域を拡大することで特性推定精度を向上できる。次に図15と図16に基づいて第5実施例を説明する。図15はより広い運転領域でのデータ計測を可能にしたポンプ制御システムの全体構成図であり、第3、4実施例とは別の実施形態を与えるものである。図15では、流量配分計算を従来と異なる方法を用いて行っている。流量配分補正計算手段1031では、図16に示すように、図3で与えられる流量配分テーブルの検索値に補正係数を乗じて、最終的な流量配分を求めるものである。例えば、図16内の補正係数k、k1を1より小さくすれば、ポンプ4の流量配分D1を本来の流量配分より小さくでき、補正係数を1より大きくすれば、ポンプ4の配水配分D1を本来の流量配分より大きくできる。補正係数k2を1より小さくすれば、ポンプ5の流量配分D2を本来の流量配分より小さくでき、補正係数を1より大きくすれば、ポンプ5の配水配分D2を本来の流量配分より大きくできる。これらの係数によりポンプ6の流量配分をもちろん操作することができる。このように補正係数を操作することでポンプ4,5,6の配水流量のレンジを広げることができる。得られた計測データを活用してポンプ推定を行う場合、その推定精度を向上できる。上述した各実施例においては、配水流量と吐出圧の2つの変数をベースにしてポンプのトータル消費エネルギー最小となるポンプ運転台数とポンプの流量配分を決定したり、ポンプの特性を、計測データから推定し、最新の推定特性に基づいてポンプ運転台数、ポンプ流量配分を決定したり、本来のパラメータを異なるパラメータでポンプを運転することで、オンラインにおいて、より幅広い運転領域での運転等を行ったが、消費エネルギー計算手段によって得られた消費電力データの推移を検証することによりポンプの劣化レベルを推定し、ポンプのメンテナンス時期が到来したことをディスプレイ装置不図示等の出力媒体に出力させ、事業者にアナウンスさせることとしてもよい。例えば、ある一定の期間の消費エネルギーの増分が所定の閾値を超えた場合、あるいは消費エネルギーに基づいて算出される電力料金の増分が所定の閾値を超えた場合に、メンテナンス時期が到来した旨をアナウンスさせることとしてもよい。この場合、事業者は適切なタイミングでポンプを交換すべきことを容易に知ることができる。"