JP2021186761A - 学習モデル生成装置、推測装置および散気量制御装置 - Google Patents

Abstract

【課題】将来の分離膜の状態を推測して安定した膜ろ過運転を行う。【解決手段】学習モデル生成装置（１）は、膜ろ過運転中に計測される、膜ろ過圧および散気量を含む運転データから導出される入力データを取得する入力データ取得部（２１）と、取得された入力データを入力とした機械学習により、分離膜の状態を推測するための学習モデル（３１）を生成する学習部（１３）と、を備える。【選択図】図３

Description

本発明は、被処理水中に浸漬して配置された分離膜の膜面散気を行いながら、分離膜を透過した処理水を得る膜ろ過処理において適用される学習モデル生成装置等に関する。

特許文献１には、散気量を制御する制御時刻において、或る過去の時点からの膜間差圧の変化量、変化率または上昇速度に基づき、事前に設定した閾値や有機物濃度から選定した目標上昇速度と比較し、散気量を決定する技術が開示されている。

特許第６３４２１０１号

膜分離装置による膜ろ過処理においては、膜表面のファウリングの進行などに伴い、膜間差圧が急上昇する現象（TMP jump）が発生することがある。しかしながら、特許文献１に開示された技術では、制御時刻において上述のファウリングの進行が起こっている場合、換言すれば、分離膜に異常がある場合、将来、膜間差圧が急上昇する現象が発生し、膜分離装置の運転に支障をきたすおそれがある。

本発明の一態様は、将来の分離膜の状態を推測して安定した膜ろ過運転を行うための学習モデル生成装置などを実現することを目的とする。

上記の課題を解決するために、本発明の一態様に係る学習モデル生成装置は、被処理水中に浸漬して配置された分離膜と、前記分離膜の膜面散気を行う散気装置とを備え、前記散気装置により前記膜面散気を行いながら前記分離膜を透過した処理水を得る膜分離装置にて行われる膜ろ過運転中に計測される、膜ろ過圧および散気量を含む運転データから導出される入力データを取得する入力データ取得部と、前記取得された入力データを入力とした機械学習により、前記分離膜の状態を推測するための学習モデルを生成する学習部と、を備える。

前記の構成によれば、膜ろ過圧および散気量を含む運転データに基づき分離膜の状態を推測する学習モデルを生成することができる。分離膜の状態が推測できれば、それに基づき散気量を制御することができるため、分離膜が正常な状態（例えば、膜間差圧の急上昇が起こる可能性を低減した状態）を維持しながら安定した膜ろ過運転を行うことができる。

また、本発明の一態様に係る学習モデル生成装置では、前記膜ろ過運転は、間欠運転であり、前記膜ろ過圧から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における、膜ろ過圧の最大値、膜ろ過圧の最小値、膜ろ過圧の標準偏差値、膜ろ過圧の平均値、および、膜間差圧、並びに、前記単位期間以前の所定期間における、膜間差圧の変動速度、膜間差圧の変動量、および、膜間差圧の変動率、の少なくともいずれかを含み、前記散気量から導出される前記入力データは、前記単位期間における前記散気量の平均値、および、前記所定期間における前記散気量の積算値、の少なくともいずれかを含んでもよい。

前記の構成によれば、膜ろ過圧および散気量から導出される種々の特徴量を有する入力データの少なくともいずれかを用いて学習モデルを生成するので、当該学習モデルの推測精度を向上させることができる。

また、本発明の一態様に係る学習モデル生成装置では、前記運転データは、さらに、前記膜ろ過運転中に計測された膜ろ過流量を含んでもよい。

前記の構成によれば、学習モデルは、さらに、膜ろ過流量から導出される入力データを用いて生成されるので、当該学習モデルの推測精度を向上させることができる。

また、本発明の一態様に係る学習モデル生成装置では、前記膜ろ過運転は、間欠運転であり、前記膜ろ過流量から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における膜ろ過流量の平均値、および、前記単位期間以前の所定期間における前記膜ろ過流量の積算値、の少なくともいずれかを含んでもよい。

前記の構成によれば、学習モデルは、さらに、膜ろ過流量から導出される、単位期間における膜ろ過流量の平均値と、所定期間における膜ろ過流量の積算値とのいずれかを用いて生成されるので、当該学習モデルの推測精度を向上させることができる。

また、本発明の一態様に係る学習モデル生成装置では、前記入力データと、当該入力データに対する前記分離膜の状態を示すラベルとが対応付けられた教師データを生成する教師データ生成部をさらに備え、前記学習部は、前記生成された教師データを用いた教師あり学習により前記学習モデルを生成してもよい。

前記の構成によれば、教師あり学習により学習モデルが生成されるので、精度の高い推測が可能となる。

また、本発明の一態様に係る学習モデル生成装置では、前記ラベルは、前記分離膜の状態が正常となる場合の前記入力データに対応付けられる正常ラベル、および、前記分離膜の状態が異常となる場合の前記入力データに対応付けられる異常ラベルを含んでもよい。

前記の構成によれば、分離膜の状態の推測において、分離膜が正常となるかまたは異常となるかを精度高く推測することができる。そのため、推測結果に基づいて、分離膜が正常となるように散気量を適切に（必要最小限に近い量に）制御することができる。

また、本発明の一態様に係る学習モデル生成装置では、前記ラベルは、前記分離膜の状態が正常と異常との中間状態となる場合の前記入力データに対応付けられる中間ラベルをさらに含んでもよい。

前記の構成によれば、分離膜が正常となるか異常となるかに加え、正常と異常との中間状態となるとの推測が可能となるので、分離膜の状態をさらに精度高く推測することができる。

また、本発明の一態様に係る学習モデル生成装置では、前記学習部は、教師なし学習により、学習結果としてクラスタを含む前記学習モデルを生成してもよい。

前記の構成によれば、十分な教師データを用意することができない状況においても、簡易に学習モデルを生成することができる。

また、本発明の一態様に係る学習モデル生成装置では、前記学習部は、前記学習モデルを用いて前記分離膜の状態を推測したときに前記学習モデルに入力されたデータを取得し、該取得したデータを入力とした機械学習により前記学習モデルを更新してもよい。

前記の構成によれば、推測工程で入力した入力データを用いた機械学習により学習モデルを更新するので、推測を行うたびに学習モデルの推測精度を向上させることができる。

上記の課題を解決するために、本発明の一態様に係る推測装置は、被処理水中に浸漬して配置された分離膜と、前記分離膜の膜面散気を行う散気装置とを備え、前記散気装置により前記膜面散気を行いながら前記分離膜を透過した処理水を得る膜分離装置にて行われる膜ろ過運転中に計測される、膜ろ過圧および散気量を含む運転データに基づき機械学習により生成された、前記分離膜の状態を推測するための学習モデルにアクセスするアクセス部と、前記膜ろ過運転中に計測された前記運転データから導出される入力データを取得する入力データ取得部と、前記アクセスした学習モデルを用いて、前記取得された入力データから、前記分離膜の状態を推測する推測部と、を備える。

前記の構成によれば、膜ろ過圧および散気量を含む運転データに基づき分離膜の状態を推測することができる。分離膜の状態が推測できれば、それに基づき散気量を制御することができるため、分離膜が正常な状態（例えば、膜間差圧の急上昇が起こる可能性を低減した状態）を維持しながら安定した膜ろ過運転を行うことができる。

また、本発明の一態様に係る推測装置は、前記膜ろ過運転は、間欠運転であり、前記入力データ取得部は、運転期間と当該運転期間に続く休止期間とから成る単位期間毎に前記運転データから導出された前記入力データを取得し、前記推測部は、前記単位期間毎に、前記取得された入力データから前記状態を推測してもよい。

前記の構成によれば、膜ろ過運転中の単位期間毎に分離膜の状態が推測されるので、分離膜の状態の急な変化を迅速に確認することができる。

また、本発明の一態様に係る推測装置は、前記膜ろ過運転は、間欠運転であり、前記膜ろ過圧から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における、膜ろ過圧の最大値、膜ろ過圧の最小値、膜ろ過圧の標準偏差値、膜ろ過圧の平均値、および、膜間差圧、並びに、前記単位期間以前の所定期間における、膜間差圧の変動速度、膜間差圧の変動量、および、膜間差圧の変動率、の少なくともいずれかを含み、前記散気量から導出される前記入力データは、前記単位期間における前記散気量の平均値、および、前記所定期間における前記散気量の積算値、の少なくともいずれかを含んでもよい。

前記の構成によれば、膜ろ過圧および散気量から導出される種々の特徴量を有する入力データの少なくともいずれかを用いて生成された学習モデルを用いて推測を行うので、推測精度を向上させることができる。

また、本発明の一態様に係る推測装置は、前記運転データは、さらに、前記膜ろ過運転中に計測された膜ろ過流量を含んでもよい。

前記の構成によれば、さらに膜ろ過流量から導出される入力データを用いて生成された学習モデルを用いて推測を行うので、推測精度を向上させることができる。

また、本発明の一態様に係る推測装置は、前記膜ろ過運転は、間欠運転であり、前記膜ろ過流量から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における膜ろ過流量の平均値、および、前記単位期間以前の所定期間における前記膜ろ過流量の積算値、の少なくともいずれかを含んでもよい。

前記の構成によれば、さらに、膜ろ過流量から導出される、単位期間における膜ろ過流量の平均値と、所定期間における膜ろ過流量の積算値とのいずれかを用いて生成され他学習モデルを用いて推測を行うので、推測精度を向上させることができる。

また、本発明の一態様に係る推測装置は、前記学習モデルは、前記分離膜の状態を示すラベルを含む教師データを用いた教師あり学習により生成されてもよい。

前記の構成によれば、教師あり学習により生成された学習モデルを用いるので、精度の高い推測が可能となる。

また、本発明の一態様に係る推測装置は、前記ラベルは、前記分離膜の状態が正常となることを示すラベル、および、前記分離膜の状態が異常となることを示すラベルを含んでもよい。

前記の構成によれば、分離膜の状態の推測において、分離膜が正常となるかまたは異常となるかを精度高く推測することができる。そのため、推測結果に基づいて、分離膜が正常となるように散気量を適切に（必要最小限に近い量に）制御することができる。

また、本発明の一態様に係る推測装置は、前記ラベルは、前記分離膜の状態が正常と異常との中間状態となることを示すラベルをさらに含んでもよい。

前記の構成によれば、分離膜が正常となるか異常となるかに加え、正常と異常との中間状態となるとの推測が可能となるので、分離膜の状態をさらに精度高く推測することができる。

また、本発明の一態様に係る推測装置は、前記学習モデルは、教師なし学習により生成される、学習結果としてクラスタを含む学習モデルであってもよい。

前記の構成によれば、十分な教師データを用意することができない状況においても、簡易に分離膜の状態の推測を行うことができる。

また、本発明の一態様に係る推測装置は、前記学習モデルの出力値は、外れ値か否かであり、前記推測部は、前記出力値が外れ値であるとき、前記分離膜の状態は異常になると推測し、前記出力値が外れ値でないとき、前記分離膜の状態は正常になると推測してもよい。

前記の構成によれば、外れ値検出が可能な学習モデルを用いて簡易に分離膜の状態を推測することができる。

また、本発明の一態様に係る散気量制御装置は、前記推測装置が前記推測した前記状態に応じて前記散気量のレベルを決定し、該決定したレベルで散気するように前記散気装置を制御してもよい。

前記の構成によれば、推測した分離膜の状態に応じて、自動的に散気量のレベルを決定するので、膜ろ過運転を行う作業者が散気量のレベルを調整する必要が無い。結果として、作業者の作業負担を低減させることができる。

また、本発明の一態様に係る散気量制御装置は、前記推測装置が推測した前記状態と、自装置が直近に決定したレベルとに応じて、前記散気量の今回のレベルを決定してもよい。

前記の構成によれば、直近の散気量レベルをも考慮して、散気量のレベルを決定するので、直近の散気量レベルを考慮しない構成と比べて、より適切な散気量のレベルを決定することができる。

本発明の各態様に係る学習モデル生成装置、推測装置および散気量制御装置は、コンピュータによって実現してもよく、この場合には、コンピュータを前記学習モデル生成装置、推測装置および散気量制御装置が備える各部（ソフトウェア要素）として動作させることにより前記学習モデル生成装置、推測装置および散気量制御装置をコンピュータにて実現させる学習モデル生成装置、推測装置および散気量制御装置の制御プログラム、およびそれを記録したコンピュータ読み取り可能な記録媒体も、本発明の範疇に入る。

本発明の一態様によれば、将来の分離膜の状態を推測して安定した膜ろ過運転を行うことができる。

本発明の実施形態１に係る分離膜の状態推測システムの概要を示す図である。 膜ろ過運転における膜ろ過圧の経時変化を示す図である。 学習モデル生成装置および推測装置の要部構成の一例を示すブロック図である。 入力データの一具体例を示す図である。 分離膜の状態を判定するための所定の条件について説明する図である。 教師データの一具体例を示す図である。 学習モデル生成処理の流れの一例を示すフローチャートである。 学習モデル更新処理の流れの一例を示すフローチャートである。 推測処理および散気量制御処理の流れの一例を示すフローチャートである。 教師データの一具体例を示す図である。 本発明の実施形態２に係る学習モデル生成装置および推測装置の要部構成の一例を示すブロック図である。 学習モデル生成処理の流れの一例を示すフローチャートである。 推測処理および散気量制御処理の流れの一例を示すフローチャートである。 本発明の参考形態に係る処理概要を示す模式図である。 参考形態に係る回帰モデル生成装置および推測装置の要部構成の一例を示すブロック図である。 推測処理および散気量制御処理の流れの一例を示すフローチャートである。

〔実施形態１〕
＜分離膜の状態推測システムの概要＞
図１は、本実施形態に係る分離膜の状態推測システム１００の概要を示す図である。状態推測システム１００は、機械学習により生成される学習モデルを用いて膜ろ過運転に用いられる分離膜９３の状態を推測したうえで、推測結果に応じて分離膜９３へ供給する散気量を制御するシステムである。分離膜９３の「状態」とは、どの程度汚染されている状態にあるかを意味する。正常な状態とは、汚染度が低く、膜間差圧が急上昇する現象が当面は発生し難い状態である。

状態推測システム１００は、学習モデル生成装置１、推測装置２、記憶装置３、運転データ取得装置４、入力データ算出装置５、散気量制御装置８、および膜分離装置９０を含み、また、記憶装置７を含んでもよい。

なお、学習モデル生成装置１、推測装置２、記憶装置３、運転データ取得装置４、入力データ算出装置５、および記憶装置７の実装方法および所在は限定されるものではないが、好ましい典型例としては、運転データ取得装置４および散気量制御装置８はＰＬＣ（Programmable Logic Controller）として実装され、推測装置２、入力データ算出装置５、および記憶装置７はエッジコンピューティングであり、学習モデル生成装置１および記憶装置３はクラウドコンピューティングである。

（膜分離装置９０）
膜分離装置９０は、被処理水中に対し分離膜を用いたろ過を行い、分離膜を透過した処理水を得る膜ろ過運転を行う装置である。当該処理水は、ろ過により不純物質が除去された被処理水と表現することもできる。

膜分離装置９０は、膜分離槽９１、分離膜９３、散気管９４、散気装置９５、ろ過水配管９６およびろ過ポンプ９７を含む。膜分離槽９１は、被処理水９２を貯留する。分離膜９３は、被処理水９２中に浸漬して配置され、被処理水９２をろ過する。ろ過水配管９６は、分離膜９３を介して膜分離槽９１に接続され、被処理水９２が分離膜９３によりろ過された処理水を流通する。ろ過ポンプ９７は、ろ過水配管９６を介して分離膜９３に接続され、処理水を流出させる。散気装置９５は、分離膜９３に付着した不純物質を剥離するための空気を供給する。換言すれば、散気装置９５は分離膜９３の膜面散気を行う。散気管９４は、分離膜９３の真下に配置され、散気装置９５から供給された空気により、分離膜９３の下方から上方に向かって流れる気泡を供給する。

膜分離槽９１は、膜分離槽９１へ流入する被処理水９２を受け入れて貯留できればよく、コンクリート、ステンレスまたは樹脂などの水漏れしない材質で形成されていればよい。また、膜分離槽９１の構造は、水漏れしない構造であればよい。

分離膜９３は、中空糸膜や平膜などの固体と液体とを分離可能な膜であればよい。分離膜９３としては、例えば、逆浸透（ＲＯ：Reverse Osmosis）膜、ナノろ過（ＮＦ：Nano Filtration）膜、限外ろ過（ＵＦ：Ultra Filtration）膜、精密ろ過（ＭＦ：Micro Filtration）膜などが挙げられるが、これに限定されない。

散気管９４は気泡を供給できる能力があればよく、その材質としては、ガラス、ステンレス、焼結金属または樹脂などを使用することができる。散気装置９５は、ブロワなどの空気を圧送できる装置であればよい。

（運転データ取得装置４）
運転データ取得装置４は、各種センサ等を用いて、膜ろ過運転中に計測される運転データを取得する。本実施形態に係る運転データは、少なくとも、膜ろ過圧および散気量を含み、さらに好ましくは、膜ろ過流量を含む。膜ろ過圧は、例えば、分離膜９３とろ過ポンプ９７との間のろ過水配管９６に配置された圧力計から取得される。散気量は、散気装置９５が供給する空気量であり散気装置９５から直接取得される。膜ろ過流量は、例えば、ろ過水配管９６に配置された流量計から取得される。運転データ取得装置４は、取得した運転データを入力データ算出装置５へ送信する。

（入力データ算出装置５）
入力データ算出装置５は、受信した運転データから、学習モデル生成装置１および推測装置２へ入力するための入力データを導出する。入力データは、運転データの特徴量を表すデータであり、運転データそのものである場合と、運転データに演算を施して得られる場合とがある。そして、学習モデルを生成するフェーズでは、入力データ算出装置５は、算出した入力データを学習モデル生成装置１へ直接送信するか、または、入力データを格納するための記憶装置７へ送信する。一方、分離膜９３の状態を推測するフェーズでは、入力データ算出装置５は、算出した入力データを推測装置２へ送信する。なお、入力データの詳細については後述する。

（学習モデル生成装置１）
学習モデル生成装置１は、受信した入力データを入力とした機械学習により、分離膜９３の状態を推測するための学習モデルを生成し、記憶装置３に記憶する。学習モデル生成の詳細については後述する。

（記憶装置３）
記憶装置３は、学習モデル生成装置１が生成した学習モデルを記憶する。なお、記憶装置３は、学習モデル以外のプログラムやデータを記憶してもよい。

（推測装置２）
推測装置２は、記憶装置３に記憶された学習モデルにアクセスし、当該学習モデルを用いて、入力データ算出装置５から受信した入力データから、分離膜９３の状態を推測する。分離膜９３の状態推測の詳細については後述する。

（散気量制御装置８）
散気量制御装置８は、推測装置２による推測結果に応じて散気装置９５の散気量のレベル（以下、単に「散気量レベル」と表記する）を決定し、決定した散気量レベルで散気するように散気装置９５を制御する。一例として、散気量制御装置８は、推測装置２による推測結果が「正常」である場合は散気量レベルを現在値より下げ、一方、推測装置２による推測結果が「異常」である場合は散気量レベルを現在値より上げる。

散気量レベルの上げ幅および下げ幅は、固定値であってもよいし、可変値であってもよい。後者の場合、散気量制御装置８は、一例として、直近の散気量レベルの上げ下げ状況に応じて今回の可変値を決定してもよい。

例えば、推測結果として「正常」を取得した場合、直近に散気量レベルを上げた回数に応じて下げ幅を小さくしたり、直近に散気量レベルを下げた回数に応じて下げ幅を大きくする。逆に、例えば、推測結果として「異常」を取得した場合、直近に散気量レベルを上げた回数に応じて上げ幅を大きくしたり、直近に散気量レベルを下げた回数に応じて下げ幅を小さくする。なお、以降、散気量制御装置８は、この例が適用されたものとして説明する。

＜膜ろ過運転のサイクル＞
図２は、膜分離装置９０が実行する膜ろ過運転中に計測される膜ろ過圧の経時変化を示す図である。図２を用いて、膜ろ過運転のサイクル（以下では、「単位期間」と表記することがある）について説明する。膜ろ過運転のサイクルは、膜ろ過運転を実行する運転期間（例えば約５分）と、当該運転期間に続く、膜ろ過運転を実行しない休止期間（例えば約１分）とから成る。膜ろ過運転は、このサイクルが繰り返されて成る間欠運転である。

状態推測システム１００では、一例として、運転期間が終了して休止期間に入る毎に、運転データ取得装置４が当該運転期間中に取得した運転データを用いて、入力データ算出装置５が入力データを導出することが好ましい。そして、これに続けて、推測装置２が分離膜９３の状態を推測し、その推測結果に応じて散気装置９５を制御する。これにより、状態推測システム１００では、膜ろ過運転の１サイクル毎に、分離膜９３の状態を推測したうえで、散気装置９５の散気量を適切に制御することができる。

＜学習モデル生成装置１の要部構成＞
図３は、学習モデル生成装置１および推測装置２の要部構成の一例を示すブロック図である。

学習モデル生成装置１は制御部１０を備えている。制御部１０は、学習モデル生成装置１の各部を統括して制御するものであり、一例として、プロセッサおよびメモリにより実現される。この例において、プロセッサは、ストレージ（不図示）にアクセスし、ストレージに格納されているプログラム（不図示）をメモリにロードし、当該プログラムに含まれる一連の命令を実行する。これにより、制御部１０の各部が構成される。

当該各部として、制御部１０は、入力データ取得部１１、教師データ生成部１２および学習部１３を含んでいる。

（入力データ取得部１１）
入力データ取得部１１は、入力データ算出装置５から直接に入力データを取得するか、または、入力データ算出装置５が算出した入力データを格納している記憶装置７から入力データを取得する。そして、入力データ取得部１１は、取得した入力データを教師データ生成部１２へ出力する。

（入力データの具体例）
図４は、運転データから導出される入力データの一具体例を示す図である。入力データ算出装置５は、運転データである膜ろ過圧から、一例として、膜ろ過圧の最大値、膜ろ過圧の最小値、膜ろ過圧の標準偏差値、膜ろ過圧の平均値、膜間差圧、膜間差圧の変動速度、膜間差圧の変動量、および膜間差圧の変動率を算出する。

膜ろ過圧の最大値（以下、「最大膜ろ過圧」と表記する）は、膜ろ過運転の或るサイクル（以下、「注目サイクル」と表記する）における膜ろ過圧の最大値である。膜ろ過圧の最小値（以下、「最小膜ろ過圧」と表記する）は、注目サイクルにおける膜ろ過圧の最小値である。膜ろ過圧の標準偏差値は、注目サイクルにおける膜ろ過圧の標準偏差値である。膜ろ過圧の平均値（以下、「平均膜ろ過圧」と表記する）は、注目サイクルにおける膜ろ過圧の平均値である。

膜間差圧（ＴＭＰ：Trans Membrane Pressure）は、分離膜９３における、被処理水９２側にかかる圧力と、処理水側にかかる圧力との差分である。膜間差圧の変動速度（以下、単に「変動速度」と表記する）は、注目サイクルにおける所定時点からの所定期間（以下、「Ｐ」と表記する）における膜間差圧の傾き（ΔＴＭＰ／ΔＴ）として算出される。なお、Ｐは数時間から数日の間で適宜選択される。一例として、変動速度は、Ｐにおける膜間差圧の経時変化における回帰モデル（直線回帰）の傾きとして算出されてもよい。このとき、変動速度は、負の値をとらないものとしてもよい。膜間差圧の変動量（以下、単に「変動量」と表記する）は、Ｐにおける変動量を指す。一例として、変動量は、所定時点におけるＴＭＰの値と、Ｐが経過した時点におけるＴＭＰの値との差分として算出される。膜間差圧の変動率（以下、単に「変動率」と表記する）は、Ｐにおける変動率を指す。一例として、変動率は、変動速度を膜間差圧で除算することにより算出される（ΔＴＭＰ／（ＴＭＰ×ΔＴ））。

また、入力データ算出装置５は、運転データである散気量から、一例として、散気量の平均値および散気量の積算値を算出する。散気量の平均値（以下、「平均散気量」と表記する）は、注目サイクルにおける散気量の平均値である。散気量の積算値（以下、積算散気量と表記する）は、Ｐにおける散気量の積算値であり、一例として、Ｐにおける平均散気量の積分値として算出される。

また、入力データ算出装置５は、運転データである膜ろ過流量から、一例として、膜ろ過流量の平均値および膜ろ過流量の積算値を算出する。膜ろ過流量の平均値（以下、「平均膜ろ過流量」と表記する）は、注目サイクルにおける膜ろ過流量の平均値である。膜ろ過流量の積算値（以下、「積算膜ろ過流量」と表記する）は、Ｐにおける膜ろ過流量の積算値であり、一例として、Ｐにおける平均膜ろ過流量の積分値として算出される。

なお、図示していないが、入力データには、入力データの導出元である運転データが取得された時刻を示す時刻情報が対応付けられているものとする。

（教師データ生成部１２）
教師データ生成部１２は、入力データと、当該入力データに対する分離膜９３の状態を示すラベルとが対応付けられた教師データを生成する。そして、教師データ生成部１２は、生成した教師データを学習部１３へ出力する。

ラベルは、例えば、将来の分離膜９３の状態が正常となる入力データに対応付ける「正常ラベル」、および、将来の分離膜９３の状態が異常となる入力データに対応付ける「異常ラベル」を含む。ラベルの種類は２種類に限られるものではなく、３種類以上であってもよい。

入力データとラベルとの対応付けは、熟練作業者等の判断による手作業で行ってもよいし、自動で行ってもよい。手作業で行う場合、教師データ生成部１２は、対応付けを行うための入出力インタフェースを提供すればよい。以下では、自動で対応付ける方法の一例について説明する。

教師データ生成部１２は、まず、注目する入力データの導出元である運転データが取得された時刻（以下、「現時刻」と表記する）における運転状況を特定する。特定条件の一例として、例えば、次の２つの条件が好ましい。（第１条件）変動速度が第１所定値未満であること。（第２条件）膜間差圧が第２所定値未満であること。これらの条件は、分離膜９３の状態が正常であることの典型条件である。これらの条件を充足すれば、分離膜９３の状態は正常とみなすことができる。なお、第１所定値および第２所定値はそれぞれ、例えば、０．０８ｋＰａ／ｈ、および、１０ｋＰａであるが、これらに限定されるものではない。

次に、教師データ生成部１２は、現時刻から第１所定時間前の時刻を示す時刻情報が対応付けられている入力データ（以下、「第１入力データ」と表記する）を特定するとともに、現時刻から第２所定時間前の時刻を示す時刻情報が対応付けられている入力データ（以下、「第２入力データ」と表記する）を特定する。第１所定時間および第２所定時間はそれぞれ、例えば、３時間および２４時間であるが、これらに限定されるものではない。

そして、教師データ生成部１２は、現時刻における第１条件の充足性を、第１入力データに対応付ける。つまり、第１入力データには、例えば、「３時間後の変動速度が第１所定値未満であるか否か」が対応付けられる。

また、教師データ生成部１２は、現時刻における第２条件の充足性を、第２入力データに対応付ける。つまり、第２入力データには、例えば、「２４時間後の膜間差圧が第２所定値未満であるか否か」が対応付けられる。

そして、教師データ生成部１２は、注目する入力データを入れ替えて以上の処理を実行することを、処理対象となる全ての入力データについて繰り返す。その結果、入力データの各々に対して、将来の運転状況として、（１）第１所定時間後の第１条件の充足性、および、（２）第２所定時間後の第２条件の充足性、が特定される。

そして、教師データ生成部１２は、上記（１）（２）の両方を充足する入力データについては、将来的に分離膜９３の状態を正常にするものと判断し、「正常ラベル」を対応付ける。一方、教師データ生成部１２は、上記（１）（２）の一方でも充足しない入力データについては、将来的に分離膜９３の状態を正常にしないものと判断し、「異常ラベル」を対応付ける。以上により、教師データ生成部１２は、入力データとラベルとの対応付けを自動で行う。

なお、特定条件は上記第１条件および上記第２条件に限られるものではない。また、特定条件の数は２つに限られるものではなく、１つでもよいし、３つ以上でもよい。特定条件が多いほど、複数の尺度に基づく正確なラベル付けが可能となるため、推測装置２の推測精度が向上することが期待される。

図５は、上述した特定条件について説明する図である。グラフ８１は、各入力データに対応付けられた第２条件の充足性、すなわち、２４時間後の膜間差圧が１０ｋＰａ未満であるか否かを可視化したものである。グラフ８２は、各入力データに対応付けられた第１条件の充足性、すなわち、３時間後の変動速度が０．０８ｋＰａ／ｈ未満であるか否かを可視化したものである。なお、図５に示す各グラフでは、条件を満たす場合を「正常」（図５の灰色点）、条件を満たさない場合を「異常」（図５の黒色点）と表記している。

第１条件および第２条件と異なる特定条件の一例として「膜間差圧の経時変化において、傾きが急激に増大していない」との条件を用いてもよい。グラフ８３は、当該条件に基づき、各入力データの膜間差圧において傾きが急激に増大しているか否かを可視化したものである。

図６は、教師データの一具体例を示す図である。教師データ生成部１２は、第１条件および第２条件の両方を充足する入力データに正常ラベルを対応付ける。つまり、３時間後の変動速度が０．０８ｋＰａ／ｈ未満、かつ、２４時間後の膜間差圧が１０ｋＰａ未満の入力データに正常ラベルを対応付ける。一方、少なくとも一方の条件を充足しない入力データに異常ラベルを対応付ける。図６では、各レコードの右端にラベルが対応付けられている。

（学習部１３）
図３に戻り、学習部１３について説明する。学習部１３は、教師データ生成部１２が生成した教師データを入力とした機械学習により、分離膜９３の状態を推測するための学習モデル３１を生成する。学習部１３は、生成した学習モデル３１を記憶装置３に格納する。学習部１３は、教師データ生成部１２から取得した複数の教師データを用い、ニューラルネットワーク（ＮＮ：Neural Network）などの既知のアルゴリズムにより学習モデル３１を生成する。これにより、入力データを入力すると、将来の分離膜９３の状態が正常となる確率（または、異常となる確率）を出力する学習モデル３１が生成される。

＜推測装置２の要部構成＞
推測装置２は制御部２０を備えている。制御部２０は、推測装置２の各部を統括して制御するものであり、一例として、プロセッサおよびメモリにより実現される。この例において、プロセッサは、ストレージ（不図示）にアクセスし、ストレージに格納されているプログラム（不図示）をメモリにロードし、当該プログラムに含まれる一連の命令を実行する。これにより、制御部２０の各部が構成される。

当該各部として、制御部２０は、入力データ取得部２１、およびアクセス部２２を含む。

（入力データ取得部２１）
入力データ取得部２１は、入力データ算出装置５から入力データを取得し、アクセス部２２へ出力する。当該入力データは、直近の膜ろ過運転において計測された運転データから導出されたものであることが好ましい。

（アクセス部２２）
アクセス部２２は、記憶装置３に記憶された学習モデル３１にアクセスする。アクセス部２２は、推測部２３を含む。

推測部２３は、アクセス部２２がアクセスした学習モデル３１を用いて、入力データ取得部２１から取得された入力データから、分離膜９３の状態を推測する。具体的には、推測部２３は、入力データを学習モデル３１へ入力した結果として学習モデル３１から出力される、分離膜９３の状態が正常となる確率（または、異常となる確率）を取得する。当該確率は、現在の運転条件で膜ろ過運転を継続した場合における、将来の分離膜９３の状態が正常となる確率である。

そして、推測部２３は、学習モデル３１から得られた確率値に基づき、将来の分離膜９３の状態が、正常となるか、異常となるかを推測する。具体的には、推測部２３は、確率値が閾値以上である場合、将来の分離膜９３の状態が正常となると推測する。一方、推測部２３は、確率値が閾値未満である場合、将来の分離膜９３の状態が異常となると推測する。閾値は例えば５０（％）であるが、この例に限定されない。なお、推測すべき分離膜９３の状態は、「正常」「異常」の２状態に限られるものではなく、３状態以上であってもよい。３状態であれば、例えば、「正常」「正常と異常との中間」「異常」である。

そして、推測部２３は、この推測結果を散気量制御装置８へ出力する。なお、推測部２３は、学習モデル３１から得られた確率値そのものを推測結果として散気量制御装置８へ出力してもよい。

なお、制御部２０（入力データ取得部２１または推測部２３）は、推測に用いた入力データを用いて学習モデル生成装置１に学習モデル３１を更新させる（つまり、再学習させる）ために、当該入力データを記憶装置７に格納してもよい。

＜学習モデル生成処理の流れ＞
図７は、学習モデル生成装置１が実行する学習モデル生成処理の流れの一例を示すフローチャートである。なお、図７に示す学習モデル生成処理は、新たに学習モデル３１を生成する処理であり、学習モデル３１を更新する処理ではない。この例では、上述した、（第１条件）変動速度が第１所定値未満であること、（第２条件）膜間差圧が第２所定値未満であること、の両方を充足するとき、分離膜９３の状態が正常であるものとする。また、入力データとラベルとの対応付けを自動で行うものとする。

まず、入力データ取得部１１は、入力データ算出装置５または記憶装置７から、複数の入力データを取得する（ステップＳ１、以下、「ステップ」の記載を省略）。入力データ取得部１１は、取得した複数の入力データを教師データ生成部１２へ出力する。

教師データ生成部１２は、各入力データについて、上述した方法により、将来の運転状況を特定する（Ｓ２）。具体的には、教師データ生成部１２は、入力データの各々について、（１）第１所定時間後の第１条件の充足性、および、（２）第２所定時間後の第２条件の充足性を特定する。

次に、教師データ生成部１２は、入力データの各々に、将来の分離膜９３の状態を示すラベルを対応付けて、教師データを生成する（Ｓ３）。具体的には、上記（１）（２）の両方を充足する入力データについては「正常ラベル」を対応付け、一方、上記（１）（２）の一方でも充足しない入力データについては「異常ラベル」を対応付ける。そして、教師データ生成部１２は、生成した教師データを学習部１３へ出力する。

次に、学習部１３は、教師データ生成部１２が生成した教師データから、分離膜９３の状態が正常となる確率（または、異常となる確率）を出力する学習モデル３１を生成する（Ｓ４）。最後に、学習部１３は、生成した学習モデル３１を記憶装置３に格納する（Ｓ５）。以上で、学習モデル生成処理は終了する。

＜学習モデル更新処理の流れ＞
図８は、学習モデル生成装置１が実行する学習モデル更新処理の流れの一例を示すフローチャートである。図８に示す学習モデル更新処理は、図７に示す学習モデル生成処理により学習モデル３１が生成された後、膜ろ過運転の実行により新たに取得した入力データに基づき、学習モデル３１を更新する処理である。なお、この例においても、上述した、（第１条件）変動速度が第１所定値未満であること、（第２条件）膜間差圧が第２所定値未満であること、の両方を充足するとき、分離膜９３の状態が正常であるものとする。また、入力データとラベルとの対応付けを自動で行うものとする。

入力データ取得部１１は、入力データ算出装置５または記憶装置７から、推測装置２が推測に用いた入力データを取得するまで待機している（Ｓ１１）。入力データを取得すると（Ｓ１１でＹＥＳ）、入力データ取得部１１は、当該入力データを教師データ生成部１２へ出力する。

教師データ生成部１２は、取得した入力データの導出元である運転データが取得された時刻における運転状況を特定する（Ｓ１２）。すなわち、教師データ生成部１２は、当該時刻における第１条件および第２条件の充足性を特定する。

続いて、教師データ生成部１２は、当該時刻から第１所定時間前の時刻を示す時刻情報が対応付けられた入力データに第１条件の充足性を対応付け、また、当該時刻から第２所定時間前の時刻を示す時刻情報が対応付けられた入力データに第２条件の充足性を対応付ける。なお、これらの入力データは過去に取得されたものであり、記憶装置７に格納されている。

教師データ生成部１２は、（１）第１所定時間後の第１条件の充足性、および、（２）第２所定時間後の第２条件の充足性の両方が対応付けられた入力データについて、上記（１）（２）の充足性に応じてラベル付けを行い、教師データとして生成する（Ｓ１３）。そして、教師データ生成部１２は、生成した教師データを学習部１３へ出力する。

学習部１３は、教師データ生成部１２が今回生成した教師データに基づき再学習を行い、学習モデル３１を更新する（Ｓ１４）。そして、学習モデル更新処理はＳ１１へ戻る。

このように学習モデルを更新することにより、推測装置２による推測処理を、最新の運転状態に即したものとすることができる。

＜推測処理および散気量制御処理の流れ＞
図９は、推測装置２が実行する推測処理および散気量制御装置８が実行する散気量制御処理の流れの一例を示すフローチャートである。

入力データ取得部２１は、入力データ算出装置５から入力データを取得するまで待機している（Ｓ２１）。入力データを取得すると（Ｓ２１でＹＥＳ）、入力データ取得部２１は、当該入力データをアクセス部２２へ出力する。

アクセス部２２は、入力データを取得すると、記憶装置３に記憶されている学習モデル３１へアクセスする（Ｓ２２）。アクセス部２２に含まれる推測部２３は、入力データを学習モデル３１へ入力し（Ｓ２３）、学習モデル３１から、例えば、分離膜９３の状態が正常となる確率を取得する（Ｓ２４）。

そして、推測部２３は、一例として、取得した確率値が閾値以上であるか否かを判定する（Ｓ２５）。閾値以上であると判定した場合（Ｓ２５でＹＥＳ）、推測部２３は、分離膜９３の状態が正常となると推測する（Ｓ２６）。一方、取得した確率の値が閾値未満であると判定した場合（Ｓ２５でＮＯ）、推測部２３は、分離膜９３の状態が異常となると推測する（Ｓ２７）。推測部２３は、推測結果を散気量制御装置８へ出力する。

散気量制御装置８は、取得した推測結果と、直近の散気量レベルとに基づき、次の散気量レベルを決定する（Ｓ２８）。そして、散気量制御装置８は、決定した散気量レベルで散気するように散気装置９５を制御する（Ｓ２９）。そして、処理はＳ２１へ戻る。

＜効果＞
以上のように、本実施形態に係る学習モデル生成装置１は、膜ろ過運転中に計測される運転データから導出される入力データを取得する入力データ取得部１１を備える。また、学習モデル生成装置１は、当該入力データを入力とした機械学習により、将来の分離膜９３の状態を推測するための学習モデル３１を生成する学習部１３を備える。

また、本実施形態に係る推測装置２は、膜ろ過運転中に計測された運転データから導出される入力データを取得する入力データ取得部２１を備える。また、推測装置２は、学習モデル３１にアクセスするアクセス部２２を備える。また、推測装置２は、学習モデル３１を用いて、入力データから、将来の分離膜９３の状態を推測する推測部２３を備える。

これにより、学習モデル生成装置１は、分離膜９３の状態を推測する学習モデル３１を生成し、また、推測装置２は、学習モデル３１を用いた分離膜９３の状態を推測することができる。分離膜９３の状態が推測できれば、それに基づき散気量を制御することができるため、分離膜９３が正常な状態（例えば、膜間差圧の急上昇が起こる可能性を低減した状態）を維持しながら安定した膜ろ過運転を行うことができる。

また、学習モデル生成装置１は、入力データと、将来の分離膜９３の状態を示すラベルとが対応付けられた教師データを生成する教師データ生成部１２を備える。学習部１３は、教師データを用いた教師あり学習により学習モデル３１を生成する。これにより、教師あり学習による学習モデル３１が生成されるので、精度の高い推測が可能となる。

また、ラベルは、正常ラベルおよび異常ラベルを含む。これにより、分離膜９３の状態の推測において、分離膜９３が正常となるかまたは異常となるかを精度高く推測することができる。そのため、推測結果に基づいて、分離膜９３が正常となるように散気量を適切に制御することができる。

また、学習部１３は、推測装置２が使用した入力データを入力とした機械学習により学習モデルを更新する。これにより、分離膜９３の状態の推測に用いた入力データにより、学習モデルを更新するので、推測を行うたびに学習モデル３１の推測精度を向上させることができる。

また、入力データ取得部２１は、運転期間と休止期間とから成る膜ろ過運転のサイクル毎に入力データを取得する。また、推測部２３は、サイクル毎に、将来の分離膜９３の状態を推測する。これにより、膜ろ過運転のサイクル毎に、将来の分離膜９３の状態が推測されるので、分離膜９３の状態の急な変化を迅速に確認することができる。

また、散気量制御装置８は、推測装置２が推測した分離膜９３の状態に応じて散気量レベルを決定し、決定した散気量レベルで散気するように散気装置９５を制御する。これにより、推測した分離膜９３の状態に応じて、自動的に散気量レベルを決定するので、膜ろ過運転を行う作業者が散気量レベルを調整する必要が無い。結果として、作業者の作業負担を低減させることができる。

散気量制御装置８は、さらに、自装置が直近に決定した散気量レベルにも応じて、今回の散気量レベルを決定する。これにより、直近の散気量レベルを考慮しない構成と比べて、より適切な散気量レベルを決定することができる。

膜ろ過圧から導出される入力データは、あるサイクルにおける、最大膜ろ過圧、最小膜ろ過圧、膜ろ過圧の標準偏差値、平均膜ろ過圧、および膜間差圧、並びに、当該サイクル以前のＰにおける変動速度、変動量および変動率の少なくともいずれかを含む。また、散気量から導出される入力データは、当該サイクルにおける平均散気量および積算散気量の少なくともいずれかを含む。また、運転データはさらに膜ろ過流量を含み、膜ろ過流量から導出される入力データは、当該サイクルにおける平均膜ろ過流量および積算膜ろ過流量の少なくともいずれかを含む。これらにより、学習モデル３１を用いた推測精度を向上させることができる。

＜実施形態１の変形例＞
上述したとおり、ラベルは、正常ラベルおよび異常ラベルの２種類に限られず、「正常と異常との中間状態」を示すラベル（以下、「中間ラベル」と表記する）を含んでいてもよい。「中間状態」とは、「正常」と言い切れる状態ではなく、かつ、「異常」と言い切れる状態でもない状態のことを意味する。

図１０は、教師データの一具体例を示す図である。図１０に示す各データは図６に示す各データと同一である。図１０に示す例では、教師データ生成部１２は、３時間後の変動速度が０．０４ｋＰａ／ｈ未満、かつ、２４時間後の膜間差圧が１０ｋＰａ未満の入力データに正常ラベルを対応付ける。また、教師データ生成部１２は、３時間後の変動速度が０．０４ｋＰａ／ｈ以上で０．０８ｋＰａ未満、かつ、２４時間後の膜間差圧が１０ｋＰａ未満の入力データに中間ラベルを対応付ける。また、教師データ生成部１２は、３時間後の変動速度が０．０８ｋＰａ／ｈ以上、かつ、２４時間後の膜間差圧が１０ｋＰａ未満の入力データに異常ラベルを対応付ける。また、教師データ生成部１２は、３時間後の変動速度がどのような値であったとしても、２４時間後の膜間差圧が１０ｋＰａ以上の入力データには異常ラベルを対応付ける。

また、推測部２３は、分離膜９３の状態として、「正常」、「異常」、「中間」のいずれかを推測する構成であってもよい。例えば、推測部２３は、確率値が第１閾値以上である場合、分離膜９３の状態が正常となると推測する。また、推測部２３は、確率値が第１閾値未満かつ第２閾値以上である場合、分離膜９３の状態は中間であると推測する。また、推測部２３は、確率値が第２閾値未満である場合、分離膜９３の状態は異常となると推測する。第１閾値は例えば７０％、第２閾値は例えば４０％であるが、これに限定されない。また、散気量制御装置８は、推測結果が「中間」である場合、現在の散気量レベルを維持してもよい。

また、推測部２３は、学習モデル３１が出力した、分離膜９３が正常となる確率値そのものを、散気量制御装置８へ出力する構成であってもよい。散気量制御装置８は、当該確率値に応じて、散気量レベルを上げるかまたは下げるか、および、散気量レベルの上げ幅または下げ幅を決定してもよい。例えば、散気量制御装置８は、確率が５０％以上である場合、散気量レベルを下げ、確率値が大きくなるほど散気量レベルの下げ幅を大きくする構成であってもよい。一方、確率が５０％未満である場合、散気量レベルを上げ、確率値が小さくなるほど散気量レベルの上げ幅を大きくする構成であってもよい。

〔実施形態２〕
本発明の他の実施形態について、以下に説明する。なお、説明の便宜上、実施形態１にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を繰り返さない。

図１１は、本実施形態に係る学習モデル生成装置１Ａおよび推測装置２Ａの要部構成の一例を示すブロック図である。

実施形態１では、学習モデル生成装置１は、いわゆる教師あり学習を行い、学習モデル３１を生成した。一方、本実施形態に係る学習モデル生成装置１Ａは、教師なし学習を行い、学習モデル３１Ａを生成する。

＜学習モデル生成装置１Ａの要部構成＞
以降、学習モデル生成装置１が備える部材と同名の部材については、相違点のみを説明する。学習モデル生成装置１Ａは制御部１０Ａを備えている。制御部１０Ａは、入力データ取得部１１および学習部１３Ａを含む。

学習部１３Ａは、入力データ取得部１１から取得した入力データを用いた教師なし学習により、学習モデル３１Ａを生成する。具体的には、学習部１３Ａは、入力データの分布を分けるための境界を定めることにより、学習結果として、分離膜９３の状態が正常となる入力データが分類されるクラスタ（以下、「正常クラスタ」と表記する）と、分離膜９３の状態が異常となる入力データが分類されるクラスタ（以下、「異常クラスタ」と表記する）とを作成する。そして、作成したクラスタを学習モデル３１Ａとして、記憶装置３Ａに格納する。

教師なし学習では、ラベル付けされていない入力データを大量に学習させることにより、入力データがどのような分布をしているか学習することができる。クラスタ分析のアルゴリズムとしては、例えば、Ｋ−ｍｅａｎｓ等などの既知のアルゴリズムを用いればよい。

＜推測装置２Ａの要部構成＞
以降、推測装置２が備える部材と同名の部材については、相違点のみを説明する。推測装置２Ａは制御部２０Ａを備えている。制御部２０Ａは、入力データ取得部２１、およびアクセス部２２Ａを含む。アクセス部２２Ａは、記憶装置３Ａに記憶された学習モデル３１Ａにアクセスする。アクセス部２２Ａは推測部２３Ａを含む。

推測部２３Ａは、アクセス部２２Ａがアクセスした学習モデル３１Ａに、入力データ取得部２１が取得した入力データを入力する。そして、推測部２３Ａは、学習モデル３１Ａから出力される出力値として、当該入力データが正常クラスタに属することを示す「正常」、または、当該入力データが異常クラスタに属することを示す「異常」のいずれかを取得する。推測部２３Ａは、取得した出力値を散気量制御装置８へ出力する。

散気量制御装置８は、取得した出力値が「正常」である場合、散気量レベルを下げる。一方、散気量制御装置８は、取得した出力値が「異常」である場合、散気量レベルを上げる。

＜学習モデル生成処理の流れ＞
図１２は、学習モデル生成装置１Ａが実行する、学習モデル生成処理の流れの一例を示すフローチャートである。なお、図１２に示す学習モデル生成処理は、新たに学習モデル３１Ａを生成する処理である。

入力データ取得部１１は、入力データ算出装置５または記憶装置７から、複数の入力データを取得する（Ｓ３１）。入力データ取得部１１は、取得した複数入力データを学習部１３Ａへ出力する。

学習部１３Ａは、取得した入力データの分布に境界を定めることにより、正常クラスタおよび異常クラスタを作成する（Ｓ３２）。そして、学習部１３Ａは、これらのクラスタを含む学習モデル３１Ａを生成し（Ｓ３３）、記憶装置３Ａへ格納する（Ｓ３４）。なお、学習部１３は、取得した教師データも記憶装置３Ａへ格納する。以上で学習モデル生成処理は終了する。

なお、図示してはいないが、学習モデル生成装置１Ａが実行する学習モデル更新処理について以下に説明する。

入力データ取得部１１は、入力データ算出装置５から入力データを取得すると、学習部１３Ａへ出力する。学習部１３Ａは、取得した入力データと、記憶装置３Ａに記憶されている入力データとから成る分布を分けるための境界を改めて定めることにより、正常クラスタおよび異常クラスタを更新する。そして、学習部１３Ａは、当該更新されたクラスタを含む学習モデル３１Ａを生成し、記憶装置３Ａに記憶されている学習モデル３１Ａを、当該生成した学習モデル３１Ａで上書きする。

＜推測処理および散気量制御処理の流れ＞
図１３は、推測装置２Ａが実行する推測処理および散気量制御装置８が実行する散気量制御処理の流れの一例を示すフローチャートである。

入力データ取得部２１は、入力データ算出装置５から入力データを取得するまで待機している（Ｓ４１）。入力データを取得すると（Ｓ４１でＹＥＳ）、入力データ取得部２１は、当該入力データをアクセス部２２Ａへ出力する。

アクセス部２２Ａは、入力データを取得すると、記憶装置３Ａに記憶されている学習モデル３１Ａへアクセスする（Ｓ４２）。推測部２３Ａは、アクセス部２２Ａがアクセスした学習モデル３１Ａに、入力データ取得部２１が取得した入力データを入力し（Ｓ４３）、学習モデル３１Ａから出力された出力値である「正常」または「異常」を取得する（Ｓ４４）。推測部２３Ａは、取得した出力値を散気量制御装置８へ出力する。

散気量制御装置８は、取得した出力値と、直近の散気量レベルとに基づき、次の散気量レベルを決定する（Ｓ４５）。そして、散気量制御装置８は、決定した散気量レベルで散気するように散気装置９５を制御する（Ｓ４６）。そして、処理はＳ４１へ戻る。

＜効果＞
本実施形態に係る学習モデル生成装置１Ａによれば、学習部１３Ａは、教師なし学習により、学習結果として、正常クラスタおよび異常クラスタを含む学習モデル３１Ａを生成する。これにより、十分な教師データを用意することができない状況においても、簡易に学習モデルを生成することができる。そのため、学習モデル生成装置１Ａによる学習および推測装置２Ａによる推測は、学習モデル生成装置１による学習および推測装置２による推測の代替として、または、その前段として実行することができる。

＜外れ値検出の例＞
学習モデル３１Ａは、学習結果として正常クラスタのみを作成し、いわゆる外れ値検出を行うものであってもよい。この例の場合、学習モデル３１Ａは、入力された入力データと、正常クラスタ内の代表点（例えば重心）との距離を特定する。そして、学習モデル３１Ａは、当該距離が所定値以上である場合は入力データが外れ値であることを出力し、所定値未満である場合は入力データが正常値である（外れ値でない）ことを出力する。そして、外れ値が出力されたとき、推測部２３Ａは分離膜９３の状態が異常になると推測し、正常値が出力されたとき、推測部２３Ａは分離膜９３の状態が正常になると推測すればよい。なお、外れ値検出のアルゴリズムとしては、局所外れ値因子法（ＬＯＦ）、１クラスサポートベクターマシン（ＯＣ−ＳＶＭ）等を用いればよい。

〔参考形態〕
本発明の参考形態について、以下に説明する。なお、説明の便宜上、実施形態１、２にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を繰り返さない。

本参考形態は、膜ろ過運転にかかる運転コストを最適化することを目的として、回帰分析による長期的な推測をパラメータを変えながら繰り返すシミュレーションを行うものである。運転コストとは、散気装置９５による散気にかかるエネルギーのコスト（以下、「エネルギーコスト」と表記する）と、汚染された分離膜９３の薬洗にかかるコスト（以下、「薬洗コスト」と表記する）とを含む。薬洗コストとは、薬洗に用いる薬品の仕入れ費と、薬洗を行う作業員の人件費とを含む。

〔処理１〕図１４を参照しながら、まず、本参考形態の回帰分析による長期的な推測について概説する。本参考形態の回帰分析は、入力データ算出装置５または記憶装置７から得られる入力データを説明変数とし、入力データに対応付けられた時刻から所定のｎ時間後（ｎは正の整数）の分離膜９３の膜間差圧に関するデータ（以下、「膜間差圧関連データ」と表記する）を目的変数とするものである。膜間差圧関連データとは、一例として、実施形態１にて説明した、膜間差圧そのもの、膜間差圧の変動速度、膜間差圧の変動量、および膜間差圧の変動率、の少なくともいずれかである。この回帰分析により、散気量の現在値を維持したときのｎ時間後の膜間差圧を推測する処理を行う。そして、推測されたｎ時間後の膜間差圧関連データを用いて入力データのうちの膜間差圧関連データを更新したデータ（以下、「更新データ」と表記する。）を生成し、更新データについて再び回帰分析を行う。この処理をＮ回繰り返す（Ｎは２以上の整数）。つまり、「回帰分析によりｎ時間後の膜間差圧関連データを推測するとともに、入力データに含まれる膜間差圧関連データをｎ時間後の膜間差圧関連データに変更することにより入力データを更新する」という処理を、Ｎ回実行する。

具体的には、繰り返しの１回目として、入力データを説明変数とした回帰分析により、入力データに対応付けられた時刻からｎ時間後の膜間差圧関連データを推測したうえで、当該膜間差圧関連データを用いて入力データのうちの膜間差圧関連データを更新した更新データＵ_{（１，１）}を生成する。

繰り返しのＸ回目として（Ｘは２以上Ｎ未満の整数）、更新データＵ_{（１，Ｘ−１）}を説明変数とした回帰分析により、入力データに対応付けられた時刻からＸ×ｎ時間後の膜間差圧関連データを推測したうえで、当該膜間差圧関連データを用いて更新データＵ_{（１，Ｘ−１）}のうちの膜間差圧関連データを更新した更新データＵ_{（１，Ｘ）}を生成する。

繰り返しのＮ回目として、更新データＵ_{（１，Ｎ−１）}を説明変数とした回帰分析により、入力データに対応付けられた時刻からＮ×ｎ時間後の膜間差圧関連データを推測する。当該膜間差圧関連データを用いて更新データＵ_{（１，Ｎ−１）}のうちの膜間差圧関連データを更新した更新データＵ_{（１，Ｎ）}を生成してもよい。

以上により、入力データに対応付けられた時刻からｎ時間後、２ｎ時間後、・・・Ｎ×ｎ時間後の、合計Ｎ個の膜間差圧を推測する。以上により、散気量の現在値を維持したときのＮ×ｎ時間後までの膜間差圧の経時変化を推測する。

〔処理２〕次に、本参考形態では、入力データのうち散気量に関するデータ（データの一部、以下、「散気量関連データ」と表記する）を変更したデータ（以下、「シミュレーションデータ」と表記する）を用いて、上記処理１を行う。散気量関連データとは、一例として、実施形態１、２にて説明した、散気量の平均値および散気量の積算値の少なくともいずれかである。この処理を、散気量関連データを変えながらＭ回行う（Ｍは２以上の整数）。つまり、互いに異なるＭ個のデータ（入力データ、および、Ｍ−１個のシミュレーションデータ）に対して、上記処理１を行う。つまり、処理２では、上記処理１を、入力データに含まれるデータの一部を異ならせてＭ回実行することにより、膜間差圧のＮ×ｎ時間後までの経時変化についてＭ個の推測結果を得る。

Ｍ回中の１回目の処理は、入力データに対する上記処理１であり、具体的には上述したとおりである。

Ｍ回中のＹ回目の処理として（Ｙは２以上Ｍ以下の整数）、入力データのうち散気量関連データを変更したシミュレーションデータＳ_Ｙに対して上記処理１を行う。具体的には、上記処理１における繰り返しの１回目として、シミュレーションデータＳ_Ｙを説明変数とした回帰分析により、シミュレーションデータＳ_Ｙに対応付けられた時刻からｎ時間後の膜間差圧関連データを推測したうえで、当該膜間差圧関連データを用いてシミュレーションデータＳ_Ｙのうちの膜間差圧関連データを更新した更新データＵ_{（Ｙ，１）}を生成する。上記処理１における繰り返しのＸ回目として、更新データＵ_{（Ｙ，Ｘ−１）}を説明変数とした回帰分析により、シミュレーションデータＳ_Ｙに対応付けられた時刻からＸ×ｎ時間後の膜間差圧関連データを推測したうえで、当該膜間差圧関連データを用いて更新データＵ_{（Ｙ，Ｘ−１）}のうちの膜間差圧関連データを更新した更新データＵ_{（Ｙ，Ｘ）}を生成する。上記処理１における繰り返しのＮ回目として、更新データＵ_{（Ｙ，Ｎ−１）}を説明変数とした回帰分析により、シミュレーションデータＳ_Ｙに対応付けられた時刻からＮ×ｎ時間後の膜間差圧関連データを推測する。

以上の処理１および処理２を行うことにより、Ｍ個の散気量の各々について、散気量の現在値を維持したときのＮ×ｎ時間後までの膜間差圧の経時変化を推測することができる。このＭ個の推測結果から、適切な運転コストとするための、薬洗のタイミングおよび散気量を決定する。

本参考形態は、ファウリングが進行した膜ろ過運転後期において膜間差圧の上昇を抑え続けるだけでは必ずしもトータルの運転コストを抑制できていない状況が発生することに鑑み、散気量の制御による分離膜９３の長期延命化を回避し、トータルの運転コストの最適化を図るものである。

図１５は、本参考形態に係る回帰モデル生成装置６および推測装置２Ｂの要部構成の一例を示すブロック図である。

＜回帰モデル生成装置６の要部構成＞
回帰モデル生成装置６は制御部６０を備えている。制御部６０は、入力データ取得部６１、対応付け部６２および回帰モデル生成部６３を含む。

入力データ取得部６１は、入力データ算出装置５または記憶装置７から入力データを取得し、取得した入力データを対応付け部６２へ出力する。

対応付け部６２は、入力データの各々に、当該入力データに対応付けられた時刻からｎ時間後の膜間差圧関連データを対応付ける。ｎの値は例えば１２または２４であるが、この例に限定されない。

対応付け部６２は、ｎ時間後の膜間差圧関連データが対応付けられた入力データを、回帰モデル生成部６３へ出力する。なお、ｎ時間後の膜間差圧関連データがまだ存在しないことにより、膜間差圧関連データが対応付けられなかった入力データは、膜間差圧関連データが取得可能となるまで対応付け部６２に保持されればよい。

回帰モデル生成部６３は、入力データを説明変数とし、ｎ時間後の膜間差圧関連データを目的変数とする回帰モデル３２を生成し、記憶装置３Ｂへ格納する。

＜推測装置２Ｂの要部構成＞
以降、推測装置２が備える部材と同名の部材については、相違点のみを説明する。推測装置２Ｂは制御部２０Ｂを含む。制御部２０Ｂは、入力データ取得部２１、アクセス部２２Ｂ、およびコスト算出部２４を含む。

アクセス部２２Ｂは、記憶装置３Ｂに記憶された回帰モデル３２にアクセスする。アクセス部２２Ｂは、推測部２３Ｂを含む。推測部２３Ｂは、まず、上述した処理１を行う。具体的には、推測部２３Ｂは、アクセス部２２Ｂがアクセスした回帰モデル３２に、入力データ取得部２１が取得した入力データを入力することにより、回帰モデル３２からｎ時間後の膜間差圧関連データを取得する。続いて、推測部２３Ｂは、取得されたｎ時間後の膜間差圧関連データを用いて入力データを更新した更新データを回帰モデル３２に入力することにより、回帰モデル３２から２ｎ時間後の膜間差圧関連データを取得する。この処理をＮ回繰り返し、Ｎ×ｎ時間後までの膜間差圧関連データを取得する。

そして、推測部２３Ｂは、上述した処理２を行う。具体的には、推測部２３Ｂは、入力データとシミュレーションデータについて、上述した処理１による回帰モデルを用いた推測を行い、各々について、Ｎ×ｎ時間後までの膜間差圧関連データを取得する。これにより、推測部２３Ｂは、Ｎ×ｎ時間後までの膜間差圧の経時変化について、Ｍ個の推測結果を得ることができる。

そして、推測部２３Ｂは、（Ａ）入力データ、（Ｂ）入力データに基づき推測されたＮ×ｎ時間後までの膜間差圧の経時変化、（Ｃ）シミュレーションデータＳ_２〜Ｓ_Ｍ、（Ｄ）シミュレーションデータＳ_２〜Ｓ_Ｍの各々に基づき推測されたＮ×ｎ時間後までの膜間差圧の経時変化、をコスト算出部２４に出力する。なお、上記（Ｃ）（Ｄ）は、それぞれＭ−１個存在する。

コスト算出部２４は、推測部２３Ｂが推測した経時変化に基づき、膜ろ過運転にかかる予想コストを算出する。具体的には、推測部２３Ｂから取得した入力データおよびシミュレーションデータの散気量関連データから、エネルギーコストを算出する。また、コスト算出部２４は、推測部２３Ｂから取得した膜間差圧の経時変化から、薬洗時期を特定する。薬洗時期は、例えば、膜間差圧の急上昇（TMP jump）が発生した時点としてもよい。そして、コスト算出部２４は、特定した薬洗時期の数に基づき、薬洗コストを算出する。

また、コスト算出部２４は、算出したエネルギーコストおよび薬洗コストの組み合わせに基づき、適切な組み合わせを特定する。コスト算出部２４は、例えば、あらかじめ定められたエネルギーコストの条件と、薬洗コストの条件とを参照し、これら条件に最も近い組み合わせを特定する。コスト算出部２４は、特定した組み合わせに対応する散気量関連データを、散気量制御装置８へ出力する。

本参考形態の散気量制御装置８は、コスト算出部２４から取得した散気量関連データに基づき散気量レベルを決定し、決定した散気量レベルで散気するように散気装置９５を制御する。

＜推測処理および散気量制御処理の流れ＞
図１６は、推測装置２Ｂが実行する推測処理および散気量制御装置８が実行する散気量制御処理の流れの一例を示すフローチャートである。

入力データ取得部２１は、入力データ算出装置５から入力データを取得するまで待機している（Ｓ５１）。入力データを取得すると（Ｓ５１でＹＥＳ）、入力データ取得部２１は、当該入力データをアクセス部２２Ｂへ出力する。

アクセス部２２Ｂは、入力データを取得すると、記憶装置３Ｂに記憶されている回帰モデル３２へアクセスする。続いて、推測部２３Ｂは、回帰分析処理を実行する（Ｓ５２）。具体的には、推測部２３Ｂは、アクセス部２２Ｂがアクセスした回帰モデル３２に、入力データ取得部２１が取得した入力データを入力することにより、回帰モデル３２からｎ時間後の膜間差圧関連データを取得する。

続いて、推測部２３Ｂは、回帰分析処理の実行回数である回帰分析回数がＮ回に到達したか否かを判定する（Ｓ５３）。Ｎ回に到達していない場合（Ｓ５３でＮＯ）、推測部２３Ｂは、データ更新処理を実行する（Ｓ５４）。具体的には、推測部２３Ｂは、取得したｎ時間後の膜間差圧関連データを用いて、入力データのうちの膜間差圧関連データを更新した更新データを生成する。そして、推測部２３Ｂは、更新データを用いてＳ５２の処理を再度実行する。なお、これ以降のＳ５２の処理の実行対象は、直近のＳ５４の処理で生成した更新データである。つまり、推測部２３Ｂは、生成した更新データを回帰モデル３２へ入力する。推測部２３Ｂは、Ｓ５３の処理において回帰分析回数がＮ回に到達したと判定するまで、Ｓ５４の処理と、これに続くＳ５２の処理との実行を繰り返す。回帰分析回数がＮ回に到達すると、推測部２３Ｂは、入力データについて、散気量の現在値を維持したときのＮ×ｎ時間後までの膜間差圧関連データを取得することとなる。

回帰分析回数がＮ回に到達した場合（Ｓ５３でＹＥＳ）、推測部２３Ｂは、シミュレーションデータの生成回数がＭ−１回に到達したか否かを判定する（Ｓ５５）。なお、Ｓ５５に初めて達したときは当該生成回数は０回であるため、推測部２３Ｂは、当該生成回数がＭ−１回に到達していないと判定する。

シミュレーションデータの生成回数がＭ−１回に到達していない場合（Ｓ５５でＮＯ）、推測部２３Ｂは、シミュレーションデータ生成処理を実行する（Ｓ５６）。具体的には、推測部２３Ｂは、入力データのうちの散気量関連データを変更し、シミュレーションデータを生成する。そして、推測部２３Ｂは、生成したシミュレーションデータに対してＳ５２からＳ５４までの処理を実行する。これにより、推測部２３Ｂは、生成したシミュレーションデータについて、散気量の現在値を維持したときのＮ×ｎ時間後までの膜間差圧関連データを取得する。また、推測部２３Ｂは、Ｓ５５の処理においてシミュレーションデータの生成回数がＭ−１回に到達したと判定するまで、Ｓ５６の処理と、これに続くＳ５２からＳ５４までの処理との実行を繰り返す。当該生成回数がＭ−１回に到達すると、推測部２３Ｂは、Ｎ×ｎ時間後までの膜間差圧関連データを、Ｍ個取得することとなる。

シミュレーションデータの生成回数がＭ−１回に到達した場合（Ｓ５５でＹＥＳ）、推測部２３Ｂは、Ｎ×ｎ時間後までの膜間差圧関連データの各々から、膜間差圧の経時変化を生成する。そして、推測部２３Ｂは、（Ａ）入力データ、（Ｂ）入力データに基づき推測されたＮ×ｎ時間後までの膜間差圧の経時変化、（Ｃ）シミュレーションデータＳ_２〜Ｓ_Ｍ、（Ｄ）シミュレーションデータＳ_２〜Ｓ_Ｍの各々に基づき推測されたＮ×ｎ時間後までの膜間差圧の経時変化、をコスト算出部２４に出力する。

続いて、コスト算出部２４は、予想コスト特定処理を実行する（Ｓ５７）。具体的には、コスト算出部２４は、推測部２３Ｂから取得した入力データおよびシミュレーションデータの散気量関連データから、エネルギーコストを算出する。また、コスト算出部２４は、推測部２３Ｂから取得した膜間差圧の経時変化から、薬洗時期を特定する。そして、コスト算出部２４は、特定した薬洗時期の数に基づき、薬洗コストを算出する。また、コスト算出部２４は、算出したエネルギーコストおよび薬洗コストの組み合わせに基づき、適切な組み合わせを特定する。コスト算出部２４は、例えば、あらかじめ定められたエネルギーコストの条件と、薬洗コストの条件とを参照し、これら条件に最も近い組み合わせを特定する。

続いて、コスト算出部２４は、散気量関連データ出力処理を実行する（Ｓ５８）。具体的には、コスト算出部２４は、予想コスト特定処理にて特定された最適な組み合わせに対応する入力データまたはシミュレーションデータを特定する。そして、コスト算出部２４は、特定した入力データまたはシミュレーションデータの散気量関連データを、散気量制御装置８へ出力する。

続いて、散気量制御装置８は、散気量制御処理を実行する（Ｓ５９）。具体的には、散気量制御装置８は、コスト算出部２４から取得した散気量関連データに基づき散気量レベルを決定し、決定した散気量レベルで散気するように散気装置９５を制御する。

＜効果＞
膜間差圧の急上昇が起こる可能性を低減した状態を維持しながら膜ろ過運転を行うことは、コスト削減や省エネルギーにつながらない場合がある。この問題に対し、本参考形態に係る推測装置２Ｂによれば、散気にかかるエネルギーコストと、薬洗コストとが適切なものとなるときの散気装置９５の動作条件を特定することができる。例えば、適切なタイミングで薬洗を行うことで、散気にかかるエネルギーコストを抑えることができ、結果として運転コストを抑えることが可能になる。つまり、推測装置２Ｂによれば、膜ろ過運転にかかるトータルコストを最適化することができる。

＜変形例＞
上述では、入力データを変更してシミュレーションデータＳ_２〜Ｓ_Ｍを生成するにあたり、変更対象となるデータが散気量関連データであるものとして説明したが、これに限定されるものではなく、入力データのうち制御可能なものであればよい。例えば、膜ろ過運転における膜ろ過流量が制御可能である場合、平均膜ろ過流量および積算膜ろ過流量を変更してもよい。

コスト算出部２４は、運転コストの選択を、膜ろ過運転の作業者による操作に基づき行ってもよい。換言すれば、適切な運転コストの指標を、作業者が指定できるようになっていてもよい。

図１６に示すフローチャートにおいて、繰り返し処理のネスト構造は、この例に限定されない。つまり、図１６では、回帰分析処理を内側のループで実行し、シミュレーションデータ生成処理を外側のループで実行しているが、これに代えて、シミュレーションデータ生成処理を内側のループで実行し、回帰分析処理を外側のループで実行してもよい。

〔実施形態および参考形態に共通の変形例〕
入力データは、膜ろ過圧から算出される、最大膜ろ過圧、最小膜ろ過圧、膜ろ過圧の標準偏差値、平均膜ろ過圧、膜間差圧、膜間差圧の変動速度、膜間差圧の変動量、および膜間差圧の変動率のうち、少なくともいずれかを含んでいればよい。また、入力データは、散気量から算出される平均散気量および積算散気量のうち、少なくともいずれかを含んでいればよい。

入力データは、図４に示すデータ以外のデータを含んでいてもよい。例えば、入力データは、被処理水９２のデータを含んでいてもよい。当該データとしては、平均水温、粘度、有機物濃度（排水有機物濃度、全有機炭素濃度、ＵＶ２６０など）、ｐＨ、汚泥濃度（ＭＬＳＳ）、浮遊懸濁物（ＳＳ）、溶存酸素濃度（ＤＯ）、酸化還元電位（ＯＲＰ）、アンモニウムイオンの濃度、および硝酸イオンの濃度、などが挙げられるが、これらに限定されない。また、入力データは、膜ろ過圧を膜ろ過流量で除算した膜ろ過抵抗値およびその変動値を含んでいてもよい。

推測装置２、２Ａ、２Ｂは、散気量制御装置８の機能を兼ね備えていてもよい。この場合、散気量制御装置８は不要である。

状態推測システム１００は、推測装置２、２Ａ、２Ｂの推測結果を運転管理者へ通知する機能を有していてもよい。これにより、膜ろ過運転の安定運転を支援することができる。

学習モデル３１、３１Ａおよび回帰モデル３２は、推測装置２、２Ａ、２Ｂの記憶部（不図示）に格納されていてもよい。

学習モデル生成装置１、１Ａおよび回帰モデル生成装置６は、学習モデルまたは回帰モデルの更新処理において、入力データを推測装置２、２Ａ、２Ｂから取得してもよい。

学習モデル生成装置１および推測装置２、学習モデル生成装置１Ａおよび推測装置２Ａ、並びに、回帰モデル生成装置６および推測装置２Ｂは、それぞれ、別体または別システムとして構成され通信可能に接続されてもよいし、一体または１つのシステムとして構成されてもよい。

〔ソフトウェアによる実現例〕
学習モデル生成装置１、１Ａ、推測装置２、２Ａ、２Ｂおよび回帰モデル生成装置６の制御ブロック（制御部１０、１０Ａ、２０、２０Ａ、２０Ｂ、６０）は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ソフトウェアによって実現してもよい。

後者の場合、学習モデル生成装置１、１Ａ、推測装置２、２Ａ、２Ｂおよび回帰モデル生成装置６は、各機能を実現するソフトウェアであるプログラムの命令を実行するコンピュータを備えている。このコンピュータは、例えば１つ以上のプロセッサを備えていると共に、上記プログラムを記憶したコンピュータ読み取り可能な記録媒体を備えている。そして、上記コンピュータにおいて、上記プロセッサが上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記プロセッサとしては、例えばＣＰＵ（Central Processing Unit）を用いることができる。上記記録媒体としては、「一時的でない有形の媒体」、例えば、ＲＯＭ（Read Only Memory）等の他、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムを展開するＲＡＭ（Random Access Memory）などをさらに備えていてもよい。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明の一態様は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。

本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。

１、１Ａ 学習モデル生成装置
２、２Ａ 推測装置
８ 散気量制御装置
１１ 入力データ取得部
１２ 教師データ生成部
１３、１３Ａ 学習部
２１ 入力データ取得部
２２、２２Ａ アクセス部
２３、２３Ａ 推測部
３１、３１Ａ 学習モデル
９０ 膜分離装置
９２ 被処理水
９３ 分離膜
９５ 散気装置

Claims (21)

1. 被処理水中に浸漬して配置された分離膜と、前記分離膜の膜面散気を行う散気装置とを備え、前記散気装置により前記膜面散気を行いながら前記分離膜を透過した処理水を得る膜分離装置にて行われる膜ろ過運転中に計測される、膜ろ過圧および散気量を含む運転データから導出される入力データを取得する入力データ取得部と、
   前記取得された入力データを入力とした機械学習により、前記分離膜の状態を推測するための学習モデルを生成する学習部と、を備える学習モデル生成装置。
2. 前記膜ろ過運転は、間欠運転であり、
   前記膜ろ過圧から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における、膜ろ過圧の最大値、膜ろ過圧の最小値、膜ろ過圧の標準偏差値、膜ろ過圧の平均値、および、膜間差圧、並びに、前記単位期間以前の所定期間における、膜間差圧の変動速度、膜間差圧の変動量、および、膜間差圧の変動率、の少なくともいずれかを含み、
   前記散気量から導出される前記入力データは、前記単位期間における前記散気量の平均値、および、前記所定期間における前記散気量の積算値、の少なくともいずれかを含む、請求項１に記載の学習モデル生成装置。
3. 前記運転データは、さらに、前記膜ろ過運転中に計測された膜ろ過流量を含む、請求項１または２に記載の学習モデル生成装置。
4. 前記膜ろ過運転は、間欠運転であり、
   前記膜ろ過流量から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における膜ろ過流量の平均値、および、前記単位期間以前の所定期間における前記膜ろ過流量の積算値、の少なくともいずれかを含む、請求項３に記載の学習モデル生成装置。
5. 前記入力データと、当該入力データに対する前記分離膜の状態を示すラベルとが対応付けられた教師データを生成する教師データ生成部をさらに備え、
   前記学習部は、前記生成された教師データを用いた教師あり学習により前記学習モデルを生成する、請求項１から４のいずれか１項に記載の学習モデル生成装置。
6. 前記ラベルは、前記分離膜の状態が正常となる場合の前記入力データに対応付けられる正常ラベル、および、前記分離膜の状態が異常となる場合の前記入力データに対応付けられる異常ラベルを含む、請求項５に記載の学習モデル生成装置。
7. 前記ラベルは、前記分離膜の状態が正常と異常との中間状態となる場合の前記入力データに対応付けられる中間ラベルをさらに含む、請求項６に記載の学習モデル生成装置。
8. 前記学習部は、教師なし学習により、学習結果としてクラスタを含む前記学習モデルを生成する、請求項１から４のいずれか１項に記載の学習モデル生成装置。
9. 前記学習部は、前記学習モデルを用いて前記分離膜の状態を推測したときに前記学習モデルに入力されたデータを取得し、該取得したデータを入力とした機械学習により前記学習モデルを更新する、請求項１から８のいずれか１項に記載の学習モデル生成装置。
10. 被処理水中に浸漬して配置された分離膜と、前記分離膜の膜面散気を行う散気装置とを備え、前記散気装置により前記膜面散気を行いながら前記分離膜を透過した処理水を得る膜分離装置にて行われる膜ろ過運転中に計測される、膜ろ過圧および散気量を含む運転データに基づき機械学習により生成された、前記分離膜の状態を推測するための学習モデルにアクセスするアクセス部と、
    前記膜ろ過運転中に計測された前記運転データから導出される入力データを取得する入力データ取得部と、
    前記アクセスした学習モデルを用いて、前記取得された入力データから、前記分離膜の状態を推測する推測部と、を備える推測装置。
11. 前記膜ろ過運転は、間欠運転であり、
    前記入力データ取得部は、運転期間と当該運転期間に続く休止期間とから成る単位期間毎に前記運転データから導出された前記入力データを取得し、
    前記推測部は、前記単位期間毎に、前記取得された入力データから前記状態を推測する、請求項１０に記載の推測装置。
12. 前記膜ろ過運転は、間欠運転であり、
    前記膜ろ過圧から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における、膜ろ過圧の最大値、膜ろ過圧の最小値、膜ろ過圧の標準偏差値、膜ろ過圧の平均値、および、膜間差圧、並びに、前記単位期間以前の所定期間における、膜間差圧の変動速度、膜間差圧の変動量、および、膜間差圧の変動率、の少なくともいずれかを含み、
    前記散気量から導出される前記入力データは、前記単位期間における前記散気量の平均値、および、前記所定期間における前記散気量の積算値、の少なくともいずれかを含む、請求項１０または１１に記載の推測装置。
13. 前記運転データは、さらに、前記膜ろ過運転中に計測された膜ろ過流量を含む、請求項１０から１２のいずれか１項に記載の推測装置。
14. 前記膜ろ過運転は、間欠運転であり、
    前記膜ろ過流量から導出される前記入力データは、運転期間と当該運転期間に続く休止期間とから成る単位期間における膜ろ過流量の平均値、および、前記単位期間以前の所定期間における前記膜ろ過流量の積算値、の少なくともいずれかを含む、請求項１３に記載の推測装置。
15. 前記学習モデルは、前記分離膜の状態を示すラベルを含む教師データを用いた教師あり学習により生成される、請求項１０から１４のいずれか１項に記載の推測装置。
16. 前記ラベルは、前記分離膜の状態が正常となることを示すラベル、および、前記分離膜の状態が異常となることを示すラベルを含む、請求項１５に記載の推測装置。
17. 前記ラベルは、前記分離膜の状態が正常と異常との中間状態となることを示すラベルをさらに含む、請求項１６に記載の推測装置。
18. 前記学習モデルは、教師なし学習により生成される、学習結果としてクラスタを含む学習モデルである、請求項１０から１４のいずれか１項に記載の推測装置。
19. 前記学習モデルの出力値は、外れ値か否かであり、
    前記推測部は、
    前記出力値が外れ値であるとき、前記分離膜の状態は異常になると推測し、
    前記出力値が外れ値でないとき、前記分離膜の状態は正常になると推測する、請求項１８に記載の推測装置。
20. 請求項１０から１９のいずれか１項に記載の推測装置が推測した前記状態に応じて前記散気量のレベルを決定し、該決定したレベルで散気するように前記散気装置を制御する、散気量制御装置。
21. 前記推測装置が推測した前記状態と、自装置が直近に決定したレベルとに応じて、前記散気量の今回のレベルを決定する、請求項２０に記載の散気量制御装置。

Images