(57)【要約】【課題】ポンプシステムにおいて、流量計を用いずに精度よく流量制御を行えるようにする。【解決手段】揚程（H）に基づいてポンプ（10）の流量（Q）を推定する第１流量推定モデル（M1）、及びモータ（20）の軸動力（L）を示す軸動力指標に基づいて流量（Q）を推定する第２流量推定モデル（M2）のうちから、管路抵抗曲線の変動に応じて、推定精度が高い方の流量推定モデルを選択し、流量（Q）の推定を行う制御部（40）を設ける。【選択図】図５
2015121153.tif 000002
【特許請求の範囲】
【請求項１】
モータ（20）で駆動されるポンプ（10）を備えたポンプシステムにおいて、
揚程（H）に基づいて上記ポンプ（10）の流量（Q）を推定する第１流量推定モデル（M1）、及び上記モータ（20）の軸動力（L）を示す軸動力指標に基づいて上記流量（Q）を推定する第２流量推定モデル（M2）のうちから、管路抵抗曲線の変動に応じて、推定精度が高い方の流量推定モデルを選択し、上記流量（Q）の推定を行う制御部（40）を備えたことを特徴とするポンプシステム。
【請求項２】
請求項１において、
上記揚程（H）と上記軸動力指標とを同一の流量（Q）において比較できるように、上記第１流量推定モデル（M1）及び上記第２流量推定モデル（M2）の少なくとも一方が正規化されており、
上記制御部（40）は、所定の回転速度（Np）における、上記揚程（H）の変化率と上記軸動力指標の変化率とを比較し、変化率が大きい方の推定モデル（M1,M2）を選択することを特徴とするポンプシステム。
【請求項３】
請求項１又は請求項２において、
上記制御部（40）は、選択した推定モデル（M1,M2）に基づいて、上記モータ（20）の回転速度を制御することを特徴とするポンプシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、モータで駆動されるポンプを備えたポンプシステムに関するものである。
【背景技術】
【０００２】
ビルなどの空気調和装置では、冷水や温水によって熱を搬送するためにポンプシステムが設けられることが多い。このようなポンプシステムでは、冷水や温水の流量制御が重要である。冷水や温水の流量を知るためには、例えば流路に電磁流量計が設けられる（例えば特許文献１を参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００３−２８４８８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、電磁流量計は一般的には高価であり、コストの観点から、常に採用できるとは限らない。
【０００５】
本発明は上記の問題に着目してなされたものであり、ポンプシステムにおいて、流量計を用いずに精度よく流量制御を行えるようにすることを目的としている。
【課題を解決するための手段】
【０００６】
上記の課題を解決するため、第１の発明は、
モータ（20）で駆動されるポンプ（10）を備えたポンプシステムにおいて、
揚程（H）に基づいて上記ポンプ（10）の流量（Q）を推定する第１流量推定モデル（M1）、及び上記モータ（20）の軸動力（L）を示す軸動力指標に基づいて上記流量（Q）を推定する第２流量推定モデル（M2）のうちから、管路抵抗曲線の変動に応じて、推定精度が高い方の流量推定モデルを選択し、上記流量（Q）の推定を行う制御部（40）を備えたことを特徴とする。
【０００７】
この構成では、管路抵抗の変化に応じて、推定精度がより高い流量推定モデルが選択される。
【０００８】
また、第２の発明は、
第1の発明において、
上記揚程（H）と上記軸動力指標とを同一の流量（Q）において比較できるように、上記第１流量推定モデル（M1）及び上記第２流量推定モデル（M2）の少なくとも一方が正規化されており、
上記制御部（40）は、所定の回転速度（Np）における、上記揚程（H）の変化率と上記軸動力指標の変化率とを比較し、変化率が大きい方の推定モデル（M1,M2）を選択することを特徴とする。
【０００９】
この構成では、揚程（H）の変化率と上記軸動力（L）の変化率とを比較し、変化率が大きい方の流量推定モデルを、推定精度がより高い流量推定モデルとして選択する。
【００１０】
また、第３の発明は、
第１又は第２の発明において、
上記制御部（40）は、選択した推定モデル（M1,M2）に基づいて、上記モータ（20）の回転速度を制御することを特徴とする。
【発明の効果】
【００１１】
第１の発明によれば、ポンプシステムにおいて、流量計を用いずに精度よく流量推定を行うことが可能になる。その結果、ポンプシステムにおいて、精度よく流量制御を行えるようになる。
【００１２】
また、第２の発明によれば、変化率が大きなモデルで推定を行うので、流量変化を高感度に検出できる。すなわち、より精度よく流量推定を行うことが可能になる。
【００１３】
また、第３の発明によれば、より精度よく流量の制御できる。
【図面の簡単な説明】
【００１４】
【図１】図１は、本発明の実施形態１におけるポンプシステムの概略構成を示す側面図である。
【図２】図２は、ポンプの特性曲線、及び管路抵抗曲線の一例を示す。
【図３】図３は、平均的な第１流量推定モデル及び第２流量推定モデルを示す曲線と、それぞれのモデルの補正例を示す。
【図４】図４は、実施形態１における流量推定を説明するフローチャートである。
【図５】図５は、本実施形態１における流量推定モデル選択の動作を説明するフローチャートである。
【図６】図６は、Ｌ−Ｑモデル及びＩ−Ｑモデルを示す曲線を２種類の回転速度について示す。
【発明を実施するための形態】
【００１５】
以下、本発明の実施形態について図面を参照しながら説明する。なお、以下の実施形態は、本質的に好ましい例示であって、本発明、その適用物、あるいはその用途の範囲を制限することを意図するものではない。
【００１６】
《発明の実施形態１》
〈概略構成〉
図１は、本発明の実施形態１におけるポンプシステム（1）の概略構成を示す側面図である。ポンプシステム（1）は、例えばビルなどの空調設備（図示は省略）に組み込まれて、冷水や温水を循環させる。空調設備は、室温の制御などを行うために、ポンプシステム（1）の制御や、冷水や温水を循環させる管路の途中に設けられたバルブの開閉制御などを行う。そのためため、空調設備は、マイクロコンピュータ等によって構成された制御装置（上位制御装置と呼ぶ）を備えている。当該上位制御装置は、ポンプシステム（1）の制御装置（後述）と、例えば、バルブの開閉状態の情報や、ポンプの流量を指示する命令等をやりとりできるように構成されている。
【００１７】
〈ポンプシステムの構成〉
ポンプシステム（1）は、図１に示すように、ポンプ（10）、モータ（20）、電力変換装置（30）、及び制御装置（40）を備えている。
【００１８】
−ポンプ−
本実施形態のポンプ（10）には、遠心ポンプを採用している。ポンプ（10）の吸入口（11）及び吐出口（12）には、配管との接続部が設けられ、空調設備の配管（図示は省略）に接続される。
【００１９】
−モータの構成−
モータ（20）は、ポンプ（10）のインペラ（図示を省略）を駆動するようになっている。モータ（20）の形式には限定はないが、例えば、誘導モータや、永久磁石埋込型モータ（いわゆるＩＰＭモータ）など、種々のモータを採用できる。なお、この例では、モータ（20）の駆動軸は、上記インペラに直結されており、ポンプ（10）の回転速度とモータ（20）の回転速度とは等しいものとする。
【００２０】
−電力変換装置−
電力変換装置（30）は、商用電源から交流電力が供給され、該交流電力を所定の交流電力に電力変換して、モータ（20）に供給する。電力変換装置（30）は、コンバータ回路、及びインバータ回路（何れも図示を省略）を備えている。コンバータ回路は、商用電源からの交流電力を直流電力に変換する。また、インバータ回路は、複数のスイッチング素子を備え、該スイッチング素子のオンオフを切り替えて、上記電力変換を行う。このオンオフの制御は、制御装置（40）が行う。電力変換装置（30）は、制御装置（40）とともに、ケーシングなどに収容されて、例えばモータ（20）の側面などに取り付けられる。勿論、電力変換装置（30）と制御装置（40）とは、別々に設置してもよいし、これらの設置場所も、モータ（20）の側面には限定されない。
【００２１】
−制御装置の構成−
制御装置（40）は、マイクロコンピュータ（41）、及び当該マイクロコンピュータ（41）を動作させるためのプログラムを含んでいる。
【００２２】
制御装置（40）では、マイクロコンピュータ（41）が、上記インバータ回路の出力を制御して、モータ（20）の運転を制御する。モータ（20）の制御には、例えばｄ−ｑ軸ベクトル制御が用いられる。制御装置（40）は、ｄ−ｑ軸ベクトル制御を行うために、モータ（20）の回転速度（Np）の情報を持っている。なお、回転速度（Np）の検出は、速度センサを設ける等、種々の方法を採用できる。
【００２３】
また、制御装置（40）は、ポンプ（10）の流量を推定し、その推定結果に基づいて、モータ（20）の回転速度（N）を制御するようになっている。制御装置（40）は、本発明の制御部の一例である。以下では、制御装置（40）における流量推定の方法について説明する。
【００２４】
−流量推定−
図２は、ポンプ（10）の特性曲線、及び管路抵抗曲線（以下、単に抵抗曲線、或いはシステムカーブとも呼ぶ）の一例を示す。同図において、Ｎ０〜Ｎ４は、ポンプ（10）の特性曲線であり、回転速度毎のポンプ（10）の特性を表している。また、Ｒ１、及びＲ２が抵抗曲線であり、ポンプ（10）が接続された管路の抵抗を表している。例えば、空調設備の管路の途中に設けられたバルブが開状態の場合の抵抗曲線がＲ１であるとすると、該バルブが閉じられることによって管路の抵抗は大きくなる。そうすると、この空調設備における抵抗曲線は、Ｒ１から、例えばＲ２へと変化する。
【００２５】
ポンプ（10）は、特性曲線と抵抗曲線の交点で動作する。例えば、回転速度がＮ０（この例では定格回転数とする）であり、管路の抵抗曲線がＲ１である場合（上記バルブが開状態の場合）は、ポンプ（10）はＡのポイント（運転点と呼ぶ）で動作する。ここで、運転点Ａの流量をＱＡとする。そして、ポンプシステム（1）で本来必要な流量がＱＢであるとすると、運転点Ａでの運転では流量が過剰である。そのため、ポンプシステム（1）では、インバータ回路を介してモータ（20）を制御してポンプ（10）の回転速度をＮ２に変化させる必要がある（図２の運転点Ｂを参照）。
【００２６】
一般的にポンプの特性として、ＱＢ／ＱＡ＝Ｎ２／Ｎ０が成り立つので、現在の運転点Ａの流量ＱＡが分かれば、流量ＱＡに基づいて、ポンプ（10）の回転速度がＮ０×ＱＢ／ＱＡとなるようにモータ（20）を制御すればよいことになる。すなわち、現在の流量（QA）を精度よく推定できれば、高価な流量計を用いずに、ポンプシステム（1）の流量を精度よく制御できるのである。
【００２７】
−流量推定モデル−
本実施形態では、ポンプ（10）の流量を推定する２つの推定モデルを用意し、これらの推定モデルのうちから、推定の精度が良い方を、抵抗曲線の変動に応じて選択し、流量の推定を行う。
【００２８】
本実施形態では、一方の流量推定モデル（第１流量推定モデル（M1））は、揚程（H）と流量（Q）との関係（Ｈ−Ｑ特性）に基づく流量推定モデル（以下ではＨ−Ｑモデルとも呼ぶ）である。また、もう一方の流量推定モデル（第２流量推定モデル（M2））は、軸動力（L）と流量（Q）との関係（Ｌ−Ｑ特性）に基づく流量推定モデル（以下では、Ｌ−Ｑモデルとも呼ぶ）である。本実施形態では、この流量推定は、マイクロコンピュータ（41）が主に行う。
【００２９】
第１流量推定モデル（M1）は、揚程（H）を流量（Q）の２次関数で表した流量推定モデルであり、揚程（H）は、次の式で表せる。
【００３０】
Ｈ＝ａ０Ｑ２＋ｂ０Ｑ＋ｃ０＋Δｃ０  ・・・式（１）
ここで、ａ０、ｂ０、及びｃ０は、ポンプ（10）を作動させた際の流量（Q）の測定値等に基づいて算出できる値であり、ポンプ性能を測定できる計測器が整った環境で測定した複数台のポンプ（10）の平均値を代表特性として用いて決定している。これらの値は、ポンプ（10）の設計時（開発時）に決定する定数である。一方、Δｃ０はポンプ（10）の個体ごとに決定する値であり、ポンプ（10）の製造時（或いは出荷の際）に決定する。第１流量推定モデル（M1）の各係数は、ポンプ（10）の回転速度に応じて定めることができる。本実施形態では、ポンプ（10）の定格回転速度（N0）時における係数に基づいて、式（１）を作成し、これをマイクロコンピュータ（41）のプログラム用のデータとして制御装置（40）に組み込んである。
【００３１】
より具体的には、ポンプシステム（1）の設計時（開発時）に複数台の、同モデルのポンプ（10）について実験（定格運転）を行って、ａ０、ｂ０、及びｃ０の値をそれぞれのポンプ（10）について求め、求めたそれらの値の平均値を代表特性として、上記プログラム用データとして組み込む。
【００３２】
Δｃ０は、工場でポンプシステム（1）を製造した際（或いは出荷の際）に、個体ごとにΔｃ０の測定を行って、求めた値を上記プログラム用データとして組み込む。すなわち、開発時に作成した平均的な第１流量推定モデル（M1）を個体ごとに補正するのである。これにより、ポンプシステム（1）では、制御装置（40）に組み込まれる第１流量推定モデル（M1）は個体ごとに異なる。
【００３３】
また、第２流量推定モデル（M2）は、ポンプ（10）の軸動力（L）を流量（Q）の２次関数で表した流量推定モデルである。軸動力（L）は、本発明の軸動力指標の一例であり、次の式で表せる。
【００３４】
Ｌ＝ｄ０Ｑ２＋ｅ０Ｑ＋ｆ０＋Δｆ０  ・・・式（２）
ここで、d０、e０、及びf０は、ポンプ（10）を作動させた際の流量（Q）の測定値等に基づいて算出する値であり、複数台のポンプ（10）の平均値から決定している。これらの値は、ポンプ（10）の設計時（開発時）に決定する定数である。本実施形態では、ポンプ（10）の定格回転速度（N0）時における係数に基づいて、式（２）を作成し、これをマイクロコンピュータ（41）のプログラム用のデータとして制御装置（40）に組み込んである。
【００３５】
より具体的には、ポンプシステム（1）の設計時（開発時）に複数台の、同モデルのポンプ（10）について、実験を行って、ｄ０、ｅ０、及びｆ０の値をそれぞれのポンプ（10）について求め、求めたそれらの値の平均値を代表特性として、上記プログラム用データとして組み込む。
【００３６】
Δｆ０は、工場でポンプシステム（1）を製造した際（或いは出荷の際）に、個体ごとにΔｆ０の測定を行って、求めた値を上記プログラム用データとして組み込む。すなわち、開発時に作成した平均的な第２流量推定モデル（M2）を個体ごとに補正するのである。つまり、ポンプシステム（1）では、制御装置（40）に組み込まれる第２流量推定モデル（M2）は、個体ごとに異なっている。図３は、平均的な第１流量推定モデル（M1）及び第２流量推定モデル（M2）を示す曲線と、それぞれのモデルの補正例を示す。
【００３７】
そして、本実施形態では、制御装置（40）は、後に詳述するように、２つの推定モデル（M1,M2）の中から、推定精度の良いと考えられる方を選択して流量の推定を行う。推定モデル（M1,M2）の選択は、管路の途中にあるバルブが操作される等して抵抗曲線が変動すると、その度にやり直すようになっている。
【００３８】
−流量推定モデルの選択−
本実施形態では、Ｈ−Ｑモデル、及びＬ−Ｑモデルを比較できるように、何れか一方を正規化し、正規化したモデルで比較して、流量（Q）の変化に対する変化率が大きい方のモデルが、流量推定モデルとして精度が高いものと考える。これは、変化率が大きなモデルの方が、流量変化を高感度に検出できるので、より精度よく流量推定を行えるからである。以下では、正規化の手順を説明する。
【００３９】
揚程（H）は、締切運転（すなわちＱ＝０）において最大楊程を示すので、Ｑ＝０のときが楊揚程（H）の最大値（Hmax）である。Hmaxは、次式（３）で表せる。
【００４０】
Ｈｍａｘ＝ｃ０＋Δｃ０  ・・・（３）
一方、軸動力（L）は、流量（Q）の増加に対してピークを持つ特性を示す（図３を参照）。そのため、次の微分値（式（４））が０となる流量（Q）において、軸動力（L）は最大値となる。
【００４１】
ｄＬ／ｄＱ＝２ｄ０Ｑ＋ｅ０  ・・・（４）
つまり、Ｑ＝−ｅ０／２ｄ０の時に、Ｌが最大となる。よって、Ｌの最大値（Lmax）は、次の式（５）で表せる。
【００４２】
Ｌｍａｘ＝−ｅ０２／４ｄ０＋ｆ０＋Δｆ０  ・・・（５）
流量（Q）の推定に使用するモデルを選択するにあたり、Ｈ−Ｑモデル、及びＬ−Ｑモデルのそれぞれの最大値を基準にして、次のように正規化を行う。
【００４３】
まず、本実施形態では、係数Ｔを次の式（６）のように定義する。
【００４４】
Ｔ＝Ｌｍａｘ／Ｈｍａｘ  ・・・（６）
そして、Ｈ−Ｑ特性（第１流量推定モデル）のモデル式（１）の両辺をＴ倍した次の式（７）を作る。すなわち、モデル式（１）を正規化するのである。
【００４５】
ＨＴ＝Ｔ（ａ０Ｑ２＋ｂ０Ｑ＋ｃ０＋Δｃ０）
＝Ａ０Ｑ２＋Ｂ０Ｑ＋Ｃ０＋ΔＣ０  ・・・（７）
ただし、Ａ０＝Ｔ・ａ０，Ｂ０＝Ｔ・ｂ０，Ｃ０＝Ｔ・ｃ０，ΔＣ０＝Ｔ・Δｃ０である。
【００４６】
式（７）をＱで微分すると次の式（８）が得られる。
【００４７】
ｄＨＴ／ｄＱ＝２Ａ０Ｑ＋Ｂ０  ・・・（８）
これによりＬ−Ｑ特性の微分値であるｄＬ／ｄＱと、最大値で正規化したＨ−Ｑ特性の微分値であるｄＨＴ／ｄＱの値を比較して大きい方のモデルで流量推定を行えばよいことになる。より具体的には、２ｄ０Ｑ＋ｅ０＞２Ａ０Ｑ＋Ｂ０が成立する場合はＬ−Ｑモデル、すなわち第２流量推定モデル（M2）を用いる。
【００４８】
また、２ｄ０Ｑ＋ｅ０≦２Ａ０Ｑ＋Ｂ０が成立する場合はＨ−Ｑモデル、すなわち第１流量推定モデル（M1）を用いる。
【００４９】
実際には、選択する推定モデル（M1,M2）は、式（４）、及び式（８）から定まる次の式（９）における、Ｑについての解を境界として決まる。
【００５０】
２ｄ０Ｑ＋ｅ０＝２Ａ０Ｑ＋Ｂ０  ・・・（９）
そこで、マイクロコンピュータ（41）における計算では、式（９）の解Ｑを求め、Ｑの値に応じて、どの流量推定モデルを使用するか決定すればよい。
【００５１】
なお、本実施形態は、初期状態の流量推定モデルは、本実施形態では、Ｈ−Ｑモデル（第１流量推定モデル（M1））とする。これは、図３から分かるように、Ｈの値に対してＱの値が一意に決まるからである。もっとも、式（４）の値が０となるＱを上限としてならば、Ｈの値に対してＱの値が一意に決まるので、初期状態において、Ｌ−Ｑモデルを使用することも可能である。
【００５２】
〈ポンプシステムの動作〉
図４は、本実施形態における流量推定を説明するフローチャートである。このフローにおける処理は、マイクロコンピュータ（41）が行う。
【００５３】
ポンプシステム（1）の運転が開始されると、制御装置（40）は、ステップ（S01）において、流量推定モデルの選択と、選択した流量推定モデルに基づいて定格時の流量（Q0）の算出を行う。図５は、本実施形態における流量推定モデル選択の動作を説明するフローチャートである。
【００５４】
図５に示すように、ステップ（S11）において、制御装置（40）は、現在のポンプ（10）の回転速度（Np）を検出する。制御装置（40）では、マイクロコンピュータ（41）が上記インバータ回路を制御しているので、マイクロコンピュータ（41）は、モータ（20）の回転速度（Np）を把握している。ステップ（S12）では、制御装置（40）は、現在の回転速度（Np）と定格回転速度（N0）との比ｋ＝Ｎｐ／Ｎ０を求める。
【００５５】
ステップ（S13）では、回転速度（Np）時における流量（Q）を求めるために、回転速度（Np）の時の揚程であるＨ（Ｎｐ）と、回転速度（Np）の時の流量であるＱ（Ｎｐ）とによって、定格運転時のＨ−Ｑモデルのモデル式（１）を書き換える。具体的には、揚程（H）がポンプ（10）の回転速度（N）の２乗に反比例すること、及び流量（Q）がポンプ（10）の回転速度（N）に比例することを利用し、定格時の揚程であるＨ（Ｎ０）と、定格時の流量であるＱ（Ｎ０）とを、回転速度（Np）の時の揚程であるＨ（Ｎｐ）、及び回転速度（Np）の時の流量であるＱ（Ｎｐ）でそれぞれ表して、式（１）に代入する。
【００５６】
詳しくは、Ｈ（Ｎ０）とＨ（Ｎｐ）の関係、及びＱ（Ｎ０）とＱ（Ｎｐ）の関係は、速度比（k）を用いると、次の式のように表せ、これらの関係式を用いる。
【００５７】
Ｈ（Ｎ０）＝Ｈ（Ｎｐ）／ｋ２  ・・・（１０）
Ｑ（Ｎ０）＝Ｑ（Ｎｐ）／ｋ  ・・・（１１）
式（１０）、式（１１）、及びＨ−Ｑモデルの式（１）によって、現在の流量（Qp）を算出することができる。
【００５８】
そして、ステップ（S14）では、ステップ（S13）で求めたＱ（Ｎｐ）と速度比（k）とを用いて定格時の流量であるＱ（Ｎ０）を算出する。
【００５９】
ステップ（S15）では、ステップ（S14）で求めた定格時のＱ（Ｎ０）を用いて、ｄＬ／ｄＱとｄＨＴ／ｄＱとを比較する（式（４）及び式（８）を参照）。ここで定格時のＱ（Ｎ０）を用いるのは、式（４）及び式（８）の各係数（ｄ０，ｅ０等）が、定格運転において定められた値だからである。
【００６０】
ステップ（S15）での比較の結果、ｄＬ／ｄＱ≦ｄＨＴ／ｄＱの場合には、ステップ（S16）に進む。ステップ（S16）では、定格時の流量（Q0）として、ステップ（S14）においてＨ−Ｑモデルによって算出した値を保持（例えばメモリーに格納）する。すなわち、ステップ（S16）に進むことで、流量推定モデルとして第１流量推定モデル（M1）が選択される。
【００６１】
一方、ｄＬ／ｄＱ＞ｄＨＴ／ｄＱの場合にはステップ（S17）に進む。ステップ（S17）に進んだ場合には、Ｌ-Ｑモデルを用いて回転速度（Np）の時のＱ（Ｎｐ）のを推定する。これには、軸動力（L）が回転速度（N）の３乗に比例することを利用する。定格時の軸動力であるＬ（Ｎ０）は、次の式で表せる。
【００６２】
Ｌ（Ｎ０）＝Ｌ（Ｎｐ）／ｋ３  ・・・（１２）
Ｑ（Ｎｐ）は、式（１１）、式（１２）、及びＬ−Ｑモデル（第２流量推定モデル（M2））の式（２）を利用して推定できる。次に、制御装置（40）は、推定したＱ（Ｎｐ）と速度比（k）の値から定格時の流量であるＱ（Ｎ０）を算出する。ステップ（S17）では、制御装置（40）は、Ｌ−Ｑモデルによって算出した流量（Q0）を保持（例えばメモリーに格納）する。すなわち、ステップ（S17）に進むことで、流量推定モデルとして第２流量推定モデル（M2）選択される。
【００６３】
そして、図４に戻って、ステップ（S02）では、制御装置（40）は、空調設備の上位制御装置と通信を行って、管路上のバルブの開閉状況などの抵抗曲線が変化する変化があったかどうかを確認する。例えば、抵抗曲線の変化があった場合には、ステップ（S03）に進む。ステップ（S03）では、制御装置（40）は、図５に示したフローを再度実行して、流量推定モデルの選択と、選択した流量推定モデルに基づく定格時の流量（Q0）の算出とを行ってステップ（S04）に進む。また、抵抗曲線の変化が無かった場合には、直ちにステップ（S04）に進む。
【００６４】
ステップ（S04）では、流量推定後、抵抗曲線が変化していない場合に流量（Q）は、回転速度（N）に比例するので、その比例係数Ｃｐを求める。具体的には、Ｃｐは、次式で表せる。
【００６５】
Ｃｐ＝Ｑ（Ｎ０）／Ｎ０＝Ｑ（Ｎｐ）／Ｎｐ
そして、上位制御装置から必要流量（Qreq）が与えられると、制御装置（40）は、Ｎｒｅｑ＝Ｑｒｅｑ／Ｃｐを計算し、上記インバータ回路を介してモータ（20）を回転速度（Nreq）に制御することによって、必要流量（Qreq）に制御する。
【００６６】
〈本実施形態の効果〉
以上のように、本実施形態によれば、揚程（H）に基づいてポンプ（10）の流量（Q）を推定する第１流量推定モデル（M1）、及びモータ（20）の軸動力（L）に基づいて流量（Q）を推定する第２流量推定モデル（M2）のうちから、管路抵抗曲線の変動に応じて一方を選択するようにした。したがって、ポンプシステムにおいて、流量計を用いずに精度よく流量制御を行うことが可能になる。
【００６７】
また、Ｈ−Ｑモデル、及びＬ−Ｑモデルは、ポンプシステム（1）の個体ごとに補正されているので、この点においても、高い精度での流量推定が可能になる。
【００６８】
《発明の実施形態２》
本発明の実施形態２では、第２流量推定モデル（M2）の他の例を説明する。具体的に、本実施形態では、モータ（20）の電流値（以下、モータ電流（I）という）に基づいて流量（Q）を推定するモデルを第２流量推定モデル（M2）として採用する。
【００６９】
〈モータ電流（I）と流量（Q）との関係を示すモデル〉
軸動力（L）は、モータ電流（I）、モータ（20）の端子電圧の実効値（以下、モータ電圧（V）という）、力率、及びモータ効率を用いて、次の式（１３）で表せる。
【００７０】
軸動力＝√３×モータ電流×モータ電圧×力率×モータ効率  ・・・（式１３）
同期モータ（例えば、いわゆるＩＰＭモータ）でポンプ（10）を駆動する場合は、モータ効率は、ほぼ一定となり、モータ（20）の力率も高くほぼ一定の値をとる。この結果、軸動力（L）は、モータ電流（I）とモータ電圧（V）の積にほぼ比例する。
【００７１】
また、同じ回転速度（N）のときに、ポンプ(10)を駆動する場合のモータ電圧（V）は、あまり大きく変化しないことから、モータ電流（I）は、軸動力（L）とほぼ比例関係になる。すなわち、モータ電流（I）は、軸動力（L）を示す軸動力指標の一例である。
【００７２】
以上より、モータ電流（I）と流量（Q）の関係（Ｉ−Ｑ特性）に基づく流量推定モデル（以下、Ｉ−Ｑモデルと呼ぶ）は、次の式（１４）で表すことができる。
【００７３】
Ｉ＝ｇ０Ｑ２＋ｈ０Ｑ＋ｉ０  ・・・（式１４）
なお、本実施形態でも、ポンプ（10）の個体差を考慮して、次の式（１５）に示すように、補正用のパラメータとしてΔｉ０を導入する。
【００７４】
Ｉ＝ｇ０Ｑ２＋ｈ０Ｑ＋ｉ０＋Δｉ０  ・・・（式１５）
このΔｉ０は、ポンプ（10）の個体ごとに決定する値であり、ポンプ（10）の製造時（或いは出荷の際）に決定する。なお、以下では、式（１５）もＩ−Ｑモデルと呼ぶ。
【００７５】
図６は、Ｌ−Ｑモデル及びＩ−Ｑモデルを示す曲線を２種類の回転速度（N）について示す。
【００７６】
以上のように、Ｉ−Ｑモデルは、同期モータ（例えばＩＰＭモータ）で駆動する場合、Ｌ−Ｑモデルの代替として使用できるので、第２流量推定モデル（M2）として使用することができる。すなわち、揚程（H）に基づいてポンプ（10）の流量（Q）を推定する第１流量推定モデル（M1）、及びモータ電流（I）に基づいて流量（Q）を推定する第２流量推定モデル（M2）のうちから、管路抵抗曲線の変動に応じて、推定精度が高くなるように一方の流量推定モデルを選択し、流量（Q）の推定を行うのである。２つのモデルの選択は、実施形態１と同様に、例えばＨ−Ｑモデルを正規化することで、Ｉ−Ｑモデルと比較できるようにし、上記変化率を比較することで実現できる。このように推定モデル（M1,M2）を選択することで、本実施形態においても、ポンプシステムにおいて、流量計を用いずに精度よく流量制御を行うことが可能になる。
【００７７】
《その他の実施形態》
なお、Ｈ−Ｑモデル（式（１）参照）を正規化する代わりに、Ｌ−Ｑモデル（式（２）参照）を正規化して、上記変化率を比較するように制御装置（40）を構成してもよい。
【００７８】
また、モータ（20）やポンプ（10）の形式は例示である。
【００７９】
また、電力変換装置（30）は、直流電源から供給された電力を所定の交流電力に変換するように構成してもよい。
【産業上の利用可能性】
【００８０】
本発明は、モータで駆動されるポンプを備えたポンプシステムとして有用である。
【符号の説明】
【００８１】
１        ポンプシステム
１０      ポンプ
２０      モータ
４０      制御部
【図１】
2015121153.tif 000003
【図２】
2015121153.tif 000004
【図３】
2015121153.tif 000005
【図４】
2015121153.tif 000006
【図５】
2015121153.tif 000007
【図６】
2015121153.tif 000008