JP6219239B2 - 水処理プラント - Google Patents

Description

本発明は，主として下水処理場の処理水の水質を制御する水処理制御システムを備えた水処理プラントに関する。

環境問題やコスト削減への対応が必須となった昨今，下水処理場においても，公共水域へ放流する処理水の水質向上、さらなる省エネ化、ICTを活用した維持管理性の向上が求められている。

下水処理場では，活性汚泥と呼ばれる微生物縣濁液により，下水中の有機物や窒素などを除去する。ブロワにより活性汚泥に空気を吹き込む反応槽を好気槽と呼ぶが，好気槽では有機物は微生物による同化・異化反応により摂取・消費され，除去される。流入下水中の窒素の多くはアンモニア性窒素の形で含まれるが，これは酸素の存在下で硝化菌により硝酸性窒素に酸化される。この硝酸性窒素の一部は返送汚泥中に残留し，上流側に返送される。その際，窒素ガスへと還元する脱窒反応が生じて，窒素成分が除去される。一方，硝化不足により放流水中にアンモニア性窒素が残留すると，放流水域の水生生物への影響や溶存酸素（ＤＯ）の消費が懸念されることから，放流先の環境維持には硝化反応の適切な管理が要求される。そのためには電力を多く消費するブロワによる風量供給を適正に制御する必要がある。風量供給量が十分でない場合は，硝化不足による環境への悪影響を引きおこす。あるいは，風量供給量が過剰な場合，硝化完了後も風量が無駄に供給されるため消費電力が増大する。

下水処理の制御には，好気槽の下流側の末端に設置したＤＯ計のＤＯを制御指標として用いるＤＯ制御がある。好気槽の下流側の末端ＤＯを一定に保つようにブロワ風量を制御することで，微生物の活性を維持し，有機物除去や硝化反応を制御する（例えば，［非特許文献１］）。

近年，活性汚泥中のアンモニア性窒素濃度を計測するアンモニア計の精度の向上や，個別の生物反応槽に適した小容量のブロワの制御性が向上したことから，生物反応槽へのブロワ風量の制御に，アンモニア計を用いた制御方式が検討されている（［非特許文献２］、［非特許文献３］）。

［特許文献１］の方法では，生物反応槽へ流入する下水の流量，各好気槽への風量，およびそれぞれの好気槽に設置したアンモニア計の測定値から，あらかじめ設定した各好気槽間の硝化量に必要な各好気槽への風量を演算し，過不足ない風量で硝化を進行させる。

特開２０１２−１７０８８３号公報

「下水道施設計画・設計指針と解説」２００９年版，発行所日本下水道協会 遠藤和広: アンモニア計とDO計を用いた送風量制御システムの開発, 第47回下水道研究発表会講演集, pp.918-920 (2010) 奥大典： 効率的散気風量制御による電力量の削減−第２報 実設備における効果検証, 第50回下水道研究発表会講演集, pp.799-802 (2013)

［非特許文献１］の方法において，ＤＯは微生物の反応活性に関係するパラメータであるが，硝化反応で考慮すべきアンモニア性窒素そのものではない。そのため，流入流量や流入水質の変動により，風量不足による処理水質の悪化あるいは過剰風量となる課題があった。

次にアンモニア計を用いた制御の場合の課題について述べる。なお、一般に、イオン電極式のアンモニア計はアンモニア濃度が低下すると消耗品である電極の寿命が短くなるため、好気槽の最後段ではなく、少し上流側の中段に設置される。

［非特許文献２］の方法は、好気槽の中段に設置したアンモニア計に基づいて好気槽の後段のＤＯ値の目標値を設定するカスケード型のフィードバック（ＦＢ）制御である。処理途中あるいは処理された水質に基づいた制御のため、流入水質が変動した場合、その影響がセンサ位置に到達するまで変動は考慮されず、風量不足による処理水質の悪化あるいは過剰風量となる課題があった。詳細について述べると、例えば、日間変動が大きい場合や急な降雨で一時的にアンモニアが希釈された下水が流入した場合である。中段のアンモニア計で計測した濃度が大きいと、全体風量は増加する。その際、希釈されている上流側では過剰処理となり、下流側に到達する前に処理水の目標値に達することがある。その結果、その後流下した下流側では最小風量下でも過剰処理となり、過剰曝気となる可能性が考えられる。一度中段のアンモニア計位置での目標値から大きく外れると、その後も風量が振動し、処理が安定しない可能性もある。もう一つの課題は維持管理性に関する課題である。制御に用いる中段地点のアンモニア濃度とDOの目標値の関係式は、オペレータ側で処理特性に合わせて試行錯誤的に、且つ、活性汚泥性状の季節変動に合わせて継続的に調整する必要があり、維持管理は必ずしも容易ではない。

［非特許文献３］の方法は、流入側と好気槽中段にアンモニア計を設置し、好気槽中段のアンモニア計によるＦＢ制御に、流入側のアンモニア計によるフィードフォワード（ＦＦ）制御を加えることで、水質負荷変動への追従を高めている。しかし、流入側のアンモニア計で計測したアンモニア濃度は現時刻のみの値である、既に好気槽へ流入した流体における流入時のアンモニア濃度（過去の流入アンモニア濃度）とは異なる。好気槽全体の風量を演算するために用いる値としては十分ではなく、必ずしも水質負荷変動への追従性を確保できず、風量不足による処理水質の悪化あるいは過剰風量となる課題があった。

［非特許文献１］，［非特許文献２］，そして、［非特許文献３］の方法は、制御のためのパラメータが手動で与えられるものである。活性汚泥の特性は時間とともに変化していくため、パラメータは試行錯誤的に調整する必要があり、維持管理の労力が増大する課題があった。

［特許文献１］の方法は，好気各槽にアンモニア計を設置することで、各反応槽間で処理されたアンモニア濃度と風量の関係を計測し、その処理性能に基づいて必要風量を演算する。上流側のアンモニア濃度の計測値から必要風量を演算するＦＦ制御の連結であり、制御の安定性を確保するためには、［特許文献１］のように全ての反応槽にアンモニア計を設置する必要がある。通常の処理場での設置が推奨される好気槽末端のＤＯ計に加えて，全ての反応槽に高価なアンモニア計を設置することは，コスト的に不利であり、その維持管理が増大するといった課題があった。

上記課題を解決するため、本発明は、流入する被処理水を酸化処理する好気槽と、前記好気槽に空気を送るブロワと、好気槽内を上流側から下流側へと流れる流下流量を推定する流下流量推定部と、前記被処理水の水質を推定する被処理水水質推定部と、前記好気槽内に設置された好気槽水質推定部と、前記ブロワの風量を計測するブロワ風量計測部と、前記ブロワの風量を演算するブロワ風量演算部と、を備えた水処理プラントにおいて、前記被処理水水質推定部及び前記好気槽水質推定部で推定する水質が前記ブロワから酸素を吹き込むことで変動する水質であり、前記好気槽水質推定部より上流側の好気槽が上流側好気槽、前記好気槽水質推定部より下流側の好気槽が下流側好気槽であって、前記ブロワ風量演算部は、前記被処理水水質推定部から前記好気槽水質推定部までの好気槽へ一定の制御周期毎に流入する前記被処理水を仮想流体塊ｉ（ｉ＝１〜Ｎ）とし、仮想流体塊ｉを目標水質に制御するための必要風量である必要累積風量を、前記必要累積風量と目標処理水質である目標処理アンモニア濃度の減少量との関係を規定する処理特性関数に基づき求めると共に、時刻ｔにおける仮想流体塊ｉへの累積風量を、少なくとも前記上流側好気槽への曝気風量及び風量配分密度に基づき求め、求めた前記必要累積風量と前記累積風量との差分を各仮想流体塊ｉの残りの滞留時間で除することで仮想流体塊ｉに必要な風量を求め、仮想流体塊ｉのＮ個の総和を前記上流側好気槽に必要な上流側風量とし、前記好気槽推定部により推定された水質及び目標水質である目標処理アンモニア濃度に基づき、前記下流側好気槽に必要な下流側風量を求め、求めた前記上流側風量と前記下流側風量を合算して求めた全風量をブロワの風量とすることを特徴とすることを特徴とするものである。

更に、本発明は水処理プラントにおいて、前記被処理水水質推定部及び前記好気槽水質推定部で推定する水質がアンモニア性窒素濃度であることを特徴とするものである。

更に、本発明は水処理プラントにおいて、前記推定する前記好気槽水質推定部に加えて、溶存酸素濃度を推定する第二好気槽水質推定部を備えたことを特徴とするものである。

更に、本発明は水処理プラントにおいて、前記ブロワ風量演算部は、前記被処理水水質推定部と、前記好気槽水質推定部で推定した水質の値と、前記上流側風量と前記下流側風量を合算して求めた全風量とに基づき、前記処理特性関数を更新することを特徴とするものである。

更に、本発明は水処理プラントにおいて、前記処理特性関数を時系列で表示する処理特性表示部を備えたことを特徴とするものである。

本発明によれば，維持管理性を向上し、下水処理の水質を適切に制御しつつ消費エネルギーを抑制できる。

実施例１の水処理プラントの構成図 流下流量と好気槽アンモニア濃度の目標値の関係 仮想流体塊による上流側風量演算 処理特性関数 処理特性関数の更新の様子 抽出した処理特性の時間変化

本発明の各実施例を図面により説明する。

図１は本発明の実施例１の構成図である。

本実施例は循環式硝化脱窒法の下水処理プラントに下水処理制御システムを適用した例である。上流側から最初沈殿池１、無酸素槽２、好気槽３、最終沈殿池４が連通し、好気槽３はブロワ５と連通する。最初沈殿池１では、流入する下水１００を重力沈降により上澄液である初沈越流水１０１と沈降物である初沈汚泥に分離する。最終沈殿池４では、流入する活性汚泥１０２を上澄液である処理水１０３と沈降物である返送汚泥１０４に分離する。返送汚泥１０４は初沈越流水１０１と混合し、活性汚泥１０２として無酸素槽２に流入する。好気槽３の末端からは活性汚泥１０２の一部が無酸素槽２へと循環液１０５として循環する。ブロワ５からは好気槽３に空気１０６が送られる。

無酸素槽２は流下流量推定部である流下流量計１０と流入水質推定部である流入アンモニア計１１が設置され、好気槽３へ流入する流入水の流量とアンモニア濃度を計測する。ここでのアンモニア濃度は溶存酸素と同様にブロワから酸素を吹き込むことで変動する水質である。好気槽３には好気槽水質推定部である好気槽アンモニア計１２が設置され、好気槽３のアンモニア濃度を計測する。ブロワ５と好気槽３を連通する配管にはブロワ風量計測部である風量計１３が設置され、好気槽３へ送られる空気の風量を計測する。

流下流量計１０、流入アンモニア計１１、好気槽アンモニア計１２、風量計１３での計測値は、ブロワ風量演算部２０に送信される。ブロワ風量演算部２０の演算結果はブロワ風量制御部２１に送信され、ブロワ５の風量はブロワ風量演算部２０で演算した風量に制御される。

実施例１での風量演算方法について説明する。実施例１では、好気槽３のうち、好気槽アンモニア計１２よりも上流側の好気槽を上流側好気槽、下流側の好気槽を下流側好気槽とする。

実施例１では、流入アンモニア濃度等の計測値および水質-必要風量関係である処理特性関数を用いてフィードフォワード演算で上流側好気槽に必要な上流側風量を演算し、好気槽アンモニア計等の計測値に基づいたフィードバック演算で下流側好気槽に必要な下流側風量を演算する。以下では下流側風量演算、上流側風量演算の順に説明する。

上流側、下流側の好気槽に対して個別の風量を適用できない場合、時刻ｔの全風量Ｑ_Ｂ（ｔ）［ｍ^３／ｈ］に対してＤ_ｕｐ［−］を好気槽全体に対する上流側好気槽への風量の配分比として、上流側好気槽への風量Ｑ_Ｂｕｐ（ｔ）、下流側好気槽への風量Ｑ_Ｂｄｗｎ（ｔ）は式（１）で表される。
下流側風量演算は、好気槽アンモニア濃度の目標値ＮＨ_{ｍｄ_ｔｇｔ}（ｔ）［ｍｇ−Ｎ／Ｌ］と好気槽アンモニア濃度の計測値ＮＨ_ｍｄ（ｔ）［ｍｇ−Ｎ／Ｌ］に対してＰＩ制御によるフィードバック演算を実施する。ＮＨ_{ｍｄ_ｔｇｔ}（ｔ）は、例えば図２のように流下流量Ｑ［ｍ^３／ｈ］に対するグラフより算出する。ここでＮＨ_{ｏｕｔ_ｔｇｔ}［ｍｇ−Ｎ／Ｌ］は処理水アンモニア濃度の目標値で、管理者の所望の水質で、本実施例では１．０ｍｇ−Ｎ／Ｌとした。流下流量Ｑが大きくなると滞留時間は短くなるため、同じ風量で処理できるアンモニア量は少なくなる。したがって、流下流量Ｑが大きい場合は好気槽アンモニア濃度の目標値ＮＨ_{ｍｄ_ｔｇｔ}を低減させることで、処理水アンモニア濃度の目標値ＮＨ_{ｏｕｔ_ｔｇｔ}を達成できるようにする。

図２より好気槽アンモニア濃度の目標値NH_{md_tgt}(t)は変動する流下流量Qによる関数として変動幅を持った値として計算される。この目標値に基づいて演算される時刻t+Δtの下流側風量Q_Bdwn(t+Δt) [m³/h]を式(2)に示す。Δt [h]は制御周期である。
C_par(Z)は離散時間の並列PIコントローラの伝達関数で、比例項のパラメータをP、積分項のパラメータをI、サンプリング時間をT_s [min](=60Δt)として式(3)で表される。
ここでの流下流量は生物反応槽を流下する流量である。本実施例では循環式硝化脱窒法を対象としているため、流下流量は流入流量、循環流量、返送流量の合計値となる。

上流側風量演算では、一定の制御周期Δｔ毎に計算系（ここでは上流側のアンモニア計１１から中間点のアンモニア計１２までの生物反応槽）に１次元的に流入する仮想流体塊を考え、これをラグランジュ的に追跡する。概念図を図３に示す。灰色で示した仮想流体塊は時刻ｔ_０で計算系に流入し、時刻ｔ_Ｎで計算系後端に到達している。流下流量Ｑ_ｉｎ（ｔ）によって仮想流体塊の大きさは異なる。時刻ｔに上流側からｉ番目の仮想流体塊ｉの位置をＸ_ｖｃ，ｉ（ｔ）［ｍ］とすると、Ｘ_ｖｃ，１（ｔ）およびＸ_ｖｃ，ｉ（ｔ＋Δｔ）は式（４）（５）で表される。
ここでＳ［ｍ^２］は生物反応槽流下方向の断面積である。それぞれの仮想流体塊は流入時に計測・演算される固有の値を保持しつつ流下していく。流入時のアンモニア濃度の計測値をＮＨ_ｉｎ（ｔ）［ｍｇ−Ｎ／Ｌ］とすると、仮想流体塊ｉに対応する流入アンモニア濃度ＮＨ_ｉｎ，ｉ（ｔ）［ｍｇ−Ｎ／Ｌ］は、式（６）で表される。
仮想流体塊ｉを目標水質に制御するための必要風量である必要累積風量Ｖ_{Ｂ_ｔｇｔ，ｉ}［ｍ^３］は図４の処理特性関数と式（７）から算出される。
処理特性関数は上流側好気槽末端に到達した仮想流体塊の情報から構築する。ΔＮＨ_ｔｇｔ（ｔ）［ｍｇ−Ｎ／Ｌ］は目標処理アンモニア濃度減少量であり、上流側のアンモニア計１１の計測値（ＮＨ _ｉｎ （ｔ））と中間点のアンモニア計１２の位置でのアンモニア濃度の目標値（ＮＨ _{ｍｄ_ｔｇｔ} （t））との差分である。

仮想流体塊ｉが計算系へ流入した時刻をｔ_０，ｉ、上流側好気槽(中間点のアンモニア計１２までの好気槽)への曝気風量をＱ_Ｂｕｐ（ｔ）［ｍ^３／ｈ］、位置Ｘ_ｉ（ｔ）での風量配分密度をＤ（Ｘ_ｉ（ｔ））［−］とすると、時刻ｔにおける仮想流体塊ｉへの累積風量Ｖ_Ｂ，ｉ（ｔ）［ｍ^３］は式（８）で表される。
ここで風量配分密度Ｄ（Ｘ_ｉ（ｔ））は上流側各好気槽へ曝気風量の配分率を表す関数で、上流側好気槽全体での配分比の平均を１とする。Ｖ_ｒｔ，ｉ（ｔ）［−］は上流側好気槽の全体積Ｖ_ａｌｌ［ｍ^３］に対する仮想流体塊ｉの体積比で、式（９）で表される。
仮想流体塊ｉは図４で求める必要累積風量Ｖ_{Ｂ_ｔｇｔ，ｉ}を持つため、累積風量Ｖ_Ｂ，ｉ（ｔ）との差分を残りの滞留時間で除した値が流体塊ｉに必要な風量となる。したがって、計算系がＮ個の仮想流体塊で満たされているとすると、時刻ｔ＋Δｔの上流側風量の演算値Ｑ_Ｂｕｐ（ｔ＋Δｔ）［ｍ^３／ｈ］は式（１０）となる。
Ｌ_ａｌｌ［ｍ］は上流側好気槽の全長である。中括弧内が流体塊ｉに必要な曝気風量であるが、実施設の運用にしたがって上限値、下限値を設けることで、分母が小さくなることによる過大値や負値を避けることができる。

ここでｉ＝１は現時刻ｔの流入アンモニア濃度に基づき求めた値であり、ｉ≧２は過去の流入アンモニア濃度に基づき求めた値である。すなわち時刻ｔ＋Δｔの上流側風量の演算値Ｑ_Ｂｕｐ（ｔ＋Δｔ）は、前記流入水水質推定部で推定した水質である流入アンモニア濃度の現時刻の値に基づき演算した必要風量と過去の値に基づき演算した必要風量を用いて演算したブロワ風量となる。

演算で求めた時刻t+Δtの上流側風量と下流側風量を合算して求めた全風量Q_B(t+Δt)を用いることで、水質の目標値への安定した追随で省エネを図れる風量を求めることができる。本実施例では、上流側風量と下流側風量をそのまま合算したが、例えば下流側のフィードバック要素を重視するなど、要求に応じてその重み付けを変化させても良い。また、複数のブロワや、上流側/下流側好気槽へ連通する弁の制御で、上流側/下流側好気槽へ個別に送る風量を制御しても良い。その場合、式(1)による全風量を上流側/下流側風量に分配する必要はなく、上流側風量と下流側風量を本実施例の方法により、それぞれ演算すればよい。

本実施例で、処理特性関数として見える化された処理特性は、活性汚泥性状の季節変動にともない変動するが、維持管理性の向上にはこの処理特性関数の更新方法が重要となる。本実施例の実運用により、上流側好気槽の末端に到達した仮想流体塊Ｎには、実測値に基づいた流入アンモニア濃度、好気槽アンモニア濃度、累積風量が蓄積する。この情報に基づいて処理特性関数の現在値を統計的に更新できる。さらにその時間変化により、異常時に対する気づきが早まる。上述の制御方式では、実績値に基づいた制御パラメータの自動更新により制御精度を担保すると同時に維持管理性を向上することが可能となる。

図５に更新の様子を示す。更新前の関数を用いて制御を実施していく毎に、流入アンモニア濃度、好気槽アンモニア濃度、累積風量から実績値の情報が蓄積していく。任意のタイミング（一定間隔でも良いし、オペレータの判断でも良い）で、実績値群に対して例えば回帰分析による近似曲線を引くことで実績値に基づいた処理特性関数に更新できる。ここでの近似曲線は、実測値を確からしく近似する曲線であればｙ＝ａｘ＋ｂやｙ＝ａｘ^−１＋ｂなどのｘの正負の累乗で表される単項式でも良いし、これらを組み合わせた多項式でもよい。また、指数や対数、三角関数などの数学関数との組み合わせでも良い。また必ずしも関数ではなく不連続で段階的な対応表でも良い。その場合は、あらかじめ用意したデータベースあるいはオペレータの判断に基づいた対応関係で処理特性関数を更新すれば良い。

処理特性の時間変化については、例えば図６の処理特性関数から求めた指標の時間変化を操作画面に表示すればよい。図６は指標として、アンモニア濃度を２５ｋｇ処理する際に必要となる風量を選び、時間変化をプロットした。一般的に低温期に硝化反応速度は小さくなるが、ここでは温度の影響は取り除いてある。これより、Ａ処理場では１０月に急激に風量が増加していることから、硝化菌量が減るなど何らかの異変があったことが推定される。処理特性を可視化して連続してモニタリングすることにより、この異変を即座に察知でき、図６のように、短期間の上昇とすることができる。またＢ処理場は通年を通じてＡ処理場より必要風量が大きい。Ａ処理場はＢ処理場より効率よく運転されていることを表しており、Ａ処理場の運用を参考に運用改善を検討するべきであることが示唆される。

以上より、本実施例の下水処理プラントによる制御システムを適用することで、処理水質の安定と省エネおよび維持管理性の向上を図ることが出来る。

本実施例では、好気槽に設置した流下流量計１０を用いたが、制御対象となる領域の生物反応槽を流下する流速が分かれば良いので、流下流量計の代替として、例えば、流入下水流量計と返送汚泥流量計と循環流量計で計測した流量の合計値としても良い。

本実施例では、適用例として循環式硝化脱窒法を取り上げたが、標準法やＡＯ法、Ａ２Ｏ法、準高度処理法やこれらにステップフィードを適用した方式など、好気槽で硝化を実施する方式であれば全ての方式に適用できる。その場合、例えば標準法やＡＯ法では、流下流量は流入下水流量と返送汚泥流量の合計となるため、流下流量計１０の代替として、流入下水流量計と返送汚泥流量計の合計値としても良い。

本実施例では上流側好気槽でのフィードフォワード制御に対して、流入アンモニア計１１によるアンモニア濃度を用いたが、下水１００や初沈越流水１０１部分に設置したアンモニア計やあるいはUV計など相関関係などからアンモニア濃度を推定する他の代替手段でも良い。あるいは、流量変動、日間変動、季節変動からのアンモニア濃度の推定値でも良い。これらの場合、好気槽に流入するアンモニア濃度を求める場合、下水１００や初沈越流水１０１のアンモニア濃度に対して、下水流量、返送汚泥流量、循環流量および、好気槽アンモニア計１２によるアンモニア濃度から推定される返送汚泥、循環液に含まれるアンモニア濃度を用いて、流入アンモニア濃度を推定すれば良い。

本実施例では下流側風量の演算としてPI制御を用いたが、PID制御や他のフィードバック制御手法でも良い。

本実施例では、好気槽アンモニア計によるフィードバック制御を実施したが、実運用では施設の経験的な運転方針から、好気槽の下流側のDOに下限値を設け、好気槽下流側に設置したDO計によるDO値が下限値を下回る場合、下限値を目標とするDO制御としてもよい。

本実施例では、好気槽水質推定部としてアンモニア計を用いたが、DO計でも良い。その場合、式(2)は目標DO値に基づくPI制御等のフィードバック制御としても良い。ここでの目標DO値は一定の設定値でも良く、あるいは時間や季節や流量などに応じて変動させても良い。時間の場合は、負荷の大きい昼間はDO設定値を大きく、小さい夜間はDO設定値を小さくすることが考えられる。季節の場合は、水温が低く反応速度が小さい冬季はDO設定値を大きく、大きい夏季はDO設定値を小さくすることが考えられる。流量の場合は、流量が大きく滞留時間が短い期間はDO設定値を大きく、滞留時間が長い期間はDO設定値を小さくすることが考えられる。

本実施例では、好気槽水質推定部としてアンモニア計を用いたが、その下流側に第二好気槽水質測定部としてDO計を設置しても良い。その場合、式(2)は目標DO値に基づくPI制御等のフィードバック制御としても良い。ここでの目標DO値は好気槽アンモニア濃度の計測値NH_md(t)と対応させるカスケード制御としても良い。例えば、好気槽アンモニア濃度の計測値NH_md(t)が大きい場合は、処理水アンモニア濃度の目標値NH_{out_tgt}までの差分が大きいため、DO設定値を大きく、好気槽アンモニア濃度の計測値NH_md(t)が小さい場合は、DO設定値を小さくすることが考えられる。

本実施例では、最初沈殿池１を設置したが、最初沈殿池１を除いた構成でも良い。また、最終沈殿池４の代替として、活性汚泥と処理水の分離に膜を用いる膜分離活性汚泥法でも良い。その場合、例えば好気槽３に膜を浸漬させても良い。
本実施例では、被処理水として下水、流入水水質で推定する水質としてアンモニア性窒素を例としてあげたが、これらに限定されるものではなく、ブロワ風量を制御する制御システムを備えた水処理プラントに関わる発明として、流入水質推定部で推定する水質はブロワから酸素を吹き込むことで酸化などにより変動する水質であれば良い。

１．最初沈殿池
２．無酸素槽
３．好気槽
４．最終沈殿池
５．ブロワ
１００．下水
１０１．初沈越流水
１０２．活性汚泥
１０３．処理水
１０４．返送汚泥
１０５．循環液
１０６．空気
１０．流入流量計
１１．流入アンモニア計
１２．好気槽アンモニア計
１３．風量計
２０．ブロワ風量演算部
２１．ブロワ風量制御部

Claims (6)

1. 流入する被処理水を酸化処理する好気槽と、
   前記好気槽に空気を送るブロワと、
   前記好気槽内を上流側から下流側へと流れる流下流量を推定する流下流量推定部と、
   前記被処理水の水質を推定する被処理水水質推定部と、
   前記好気槽内に設置された好気槽水質推定部と、
   前記ブロワの風量を計測するブロワ風量計測部と、
   前記ブロワの風量を演算するブロワ風量演算部と、を備えた水処理プラントにおいて、
   前記被処理水水質推定部及び前記好気槽水質推定部で推定する水質が前記ブロワから酸素を吹き込むことで変動する水質であり、前記好気槽水質推定部より上流側の好気槽が上流側好気槽、前記好気槽水質推定部より下流側の好気槽が下流側好気槽であって、
   前記ブロワ風量演算部は、
   前記被処理水水質推定部から前記好気槽水質推定部までの好気槽へ一定の制御周期毎に流入する前記被処理水を仮想流体塊ｉ（ｉ＝１〜Ｎ）とし、仮想流体塊ｉを目標水質に制御するための必要風量である必要累積風量を、前記必要累積風量と目標処理水質である目標処理アンモニア濃度の減少量との関係を規定する処理特性関数に基づき求めると共に、時刻ｔにおける仮想流体塊ｉへの累積風量を、少なくとも前記上流側好気槽への曝気風量及び風量配分密度に基づき求め、求めた前記必要累積風量と前記累積風量との差分を各仮想流体塊ｉの残りの滞留時間で除することで仮想流体塊ｉに必要な風量を求め、仮想流体塊ｉのＮ個の総和を前記上流側好気槽に必要な上流側風量とし、
   前記好気槽水質推定部により推定された水質及び目標水質である目標処理アンモニア濃度に基づき、前記下流側好気槽に必要な下流側風量を求め、
   求めた前記上流側風量と前記下流側風量を合算して求めた全風量をブロワの風量とすることを特徴とする水処理プラント。
2. 請求項１に記載の水処理プラントにおいて、
   前記被処理水水質推定部及び前記好気槽水質推定部で推定する水質がアンモニア性窒素濃度であることを特徴とする水処理プラント。
3. 請求項１又は請求項２に記載の水処理プラントにおいて、
   前記風量配分密度は、前記上流側好気槽への曝気風量の配分率を表す関数であることを特徴とする水処理プラント。
4. 請求項２又は請求項３に記載の水処理プラントにおいて、
   前記推定する前記好気槽水質推定部に加えて、溶存酸素濃度を推定する第二好気槽水質推定部を備えたことを特徴とする水処理プラント。
5. 請求項１から請求項４のうち、いずれか１項に記載の水処理プラントにおいて、
   前記ブロワ風量演算部は、
   前記被処理水水質推定部で推定した水質の値と、前記好気槽水質推定部で推定した水質の値と、前記上流側風量と前記下流側風量を合算して求めた全風量とに基づき、前記処理特性関数を更新することを特徴とする水処理プラント。
6. 請求項５に記載の水処理プラントにおいて、
   前記処理特性関数を時系列で表示する処理特性表示部を備えたことを特徴とする水処理プラント。