JPH0356115B2 - - Google Patents
Info
- Publication number
- JPH0356115B2 JPH0356115B2 JP59118913A JP11891384A JPH0356115B2 JP H0356115 B2 JPH0356115 B2 JP H0356115B2 JP 59118913 A JP59118913 A JP 59118913A JP 11891384 A JP11891384 A JP 11891384A JP H0356115 B2 JPH0356115 B2 JP H0356115B2
- Authority
- JP
- Japan
- Prior art keywords
- standard deviation
- sludge
- meter
- water
- sewage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 52
- 239000010802 sludge Substances 0.000 claims description 35
- 239000010865 sewage Substances 0.000 claims description 28
- 238000005273 aeration Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000005856 abnormality Effects 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 8
- 238000012806 monitoring device Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 239000005416 organic matter Substances 0.000 description 18
- 239000007787 solid Substances 0.000 description 10
- 230000001364 causal effect Effects 0.000 description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 6
- 229910052760 oxygen Inorganic materials 0.000 description 6
- 239000001301 oxygen Substances 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004062 sedimentation Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 230000001932 seasonal effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000001139 pH measurement Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- -1 as shown in Figure 1 Substances 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
- Y02W10/00—Technologies for wastewater treatment
- Y02W10/10—Biological treatment of water, waste water, or sewage
Landscapes
- Activated Sludge Processes (AREA)
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は都市下水や産業廃水のように有機物を
含む汚水を浄化する下水処理システムの監視装置
に関するものである。
含む汚水を浄化する下水処理システムの監視装置
に関するものである。
一般に都市下水の下水処理システムは、第１図
に示すように有機物質を含む下水を管路１より曝
気槽２に導入し、曝気槽２には、沈殿池３から引
抜いた返送汚泥が道路４を介して返送されるとと
もに、散気管５より曝気槽内にブロワー６から圧
送された空気が供給され、管路１からの流入下水
と管路４からの返送汚泥は曝気槽内で攪拌混合さ
れ、下水中の有機物質は活性汚泥により沈降性の
良い汚泥に変換して沈殿池３に導かれる。
に示すように有機物質を含む下水を管路１より曝
気槽２に導入し、曝気槽２には、沈殿池３から引
抜いた返送汚泥が道路４を介して返送されるとと
もに、散気管５より曝気槽内にブロワー６から圧
送された空気が供給され、管路１からの流入下水
と管路４からの返送汚泥は曝気槽内で攪拌混合さ
れ、下水中の有機物質は活性汚泥により沈降性の
良い汚泥に変換して沈殿池３に導かれる。
沈殿池３においては処理水と汚泥に沈降分離さ
れ、清澄水は管路７から放流され、沈殿池に沈殿
した汚泥の大部分は返送汚泥ポンプ８により引抜
かれて管路４を介して曝気槽内に返送され残りの
汚泥は引抜泥汚ポンプ９により排出される。
れ、清澄水は管路７から放流され、沈殿池に沈殿
した汚泥の大部分は返送汚泥ポンプ８により引抜
かれて管路４を介して曝気槽内に返送され残りの
汚泥は引抜泥汚ポンプ９により排出される。
以上述べた下水処理システムにおいて、各種水
質データを計測、収集し、表示することは公知で
ある。
質データを計測、収集し、表示することは公知で
ある。
また、前記水質データをシステムの状態監視に
利用するため、上限値、下限値を計測値と共に表
示したり、あらかじめ設定された値を越えたり、
下回つたりした場合に警報を発するようにした装
置も公用されている。
利用するため、上限値、下限値を計測値と共に表
示したり、あらかじめ設定された値を越えたり、
下回つたりした場合に警報を発するようにした装
置も公用されている。
ところが、このような従来の監視装置では、処
理水の水質の悪化といつたような事故に対して
は、オペレータが水質項目間の因果関係を経験等
にもとづいて考えながら各水質項目を点検して行
く必要があり、未熟練オペレータでは正常に復旧
させるのが非常に困難であつた。
理水の水質の悪化といつたような事故に対して
は、オペレータが水質項目間の因果関係を経験等
にもとづいて考えながら各水質項目を点検して行
く必要があり、未熟練オペレータでは正常に復旧
させるのが非常に困難であつた。
そこで、水質項目間の因果関係をあらかじめ設
定しておき、ある水質項目を表示する際、その因
果関係に従つてその水質項目に影響する項目を索
引して表示させることにすることを提案したもの
として特開昭56−111087がある。
定しておき、ある水質項目を表示する際、その因
果関係に従つてその水質項目に影響する項目を索
引して表示させることにすることを提案したもの
として特開昭56−111087がある。
しかしながら該提案のように、個々の下水処理
システムごとの特徴や、プロセスの季節的な変動
を考えずに、その因果関係を画一的に設定してお
くことは、危険であり、しかも因果関係自体を把
握することもきわめて困難であると言わざるを得
ない。
システムごとの特徴や、プロセスの季節的な変動
を考えずに、その因果関係を画一的に設定してお
くことは、危険であり、しかも因果関係自体を把
握することもきわめて困難であると言わざるを得
ない。
本発明は、このような事情に鑑みてなされたも
ので、個々の下水処理システムの特徴や季節変動
が処理プロセスに与える影響を考慮して、下水処
理システムの異常を検知し、その原因を表示する
下水処理システムの監視装置を提供することを目
的とする。
ので、個々の下水処理システムの特徴や季節変動
が処理プロセスに与える影響を考慮して、下水処
理システムの異常を検知し、その原因を表示する
下水処理システムの監視装置を提供することを目
的とする。
上記目的を達成するため、本発明では、他の水
質項目（以後、システム変数と称する）が、ある
システム変数に影響を与える度合（以後、寄与率
と称する）を、自己回帰モデルから算出し、その
寄与率の大きい順に因果関係が強いとみなして、
該寄与率の大きなシステム変数の順に、そのシス
テム変数の現在値が異常がないかどうかを検知し
て表示するような下水処理システムの表示装置を
構成したことを特長としている。
質項目（以後、システム変数と称する）が、ある
システム変数に影響を与える度合（以後、寄与率
と称する）を、自己回帰モデルから算出し、その
寄与率の大きい順に因果関係が強いとみなして、
該寄与率の大きなシステム変数の順に、そのシス
テム変数の現在値が異常がないかどうかを検知し
て表示するような下水処理システムの表示装置を
構成したことを特長としている。
具体的には、流入下水量、流入水のPH、流入水
有機物濃度、溶存酸素濃度、曝気槽出口の水温、
曝気槽内混合液浮遊物質濃度（以下、MLSSとい
う）、放流水の懸濁物質、放流水のPH、放流水有
機物濃度、返送汚泥量、返送汚泥濃度、余剰汚泥
引抜量、曝気風量を測定する測定装置と、これら
測定装置の検出信号からプロセスの状態を表わす
統計的自己回帰モデルの係数を推定演算する係数
演算装置と、推定した自己回帰モデルから清澄水
の有機物濃度、懸濁物質濃度およびMLSSに影響
を及ぼすシステム変数を選別する選択装置と、前
記測定装置の出力信号から標準偏差を算出する標
準偏差演算装置と、前記選択装置および前記標準
偏差演算装置の出力信号を記憶する記憶装置と、
標準偏差を用いて計測値の検定を行い計測値が異
常であれば前記選択装置によつて選択されたシス
テム変数について点検を実施する異常検出装置
と、その点検結果を表示する表示装置とを備え
て、異常事態の発生と、その推定原因の調査結果
を知らしめようとするものである。
有機物濃度、溶存酸素濃度、曝気槽出口の水温、
曝気槽内混合液浮遊物質濃度（以下、MLSSとい
う）、放流水の懸濁物質、放流水のPH、放流水有
機物濃度、返送汚泥量、返送汚泥濃度、余剰汚泥
引抜量、曝気風量を測定する測定装置と、これら
測定装置の検出信号からプロセスの状態を表わす
統計的自己回帰モデルの係数を推定演算する係数
演算装置と、推定した自己回帰モデルから清澄水
の有機物濃度、懸濁物質濃度およびMLSSに影響
を及ぼすシステム変数を選別する選択装置と、前
記測定装置の出力信号から標準偏差を算出する標
準偏差演算装置と、前記選択装置および前記標準
偏差演算装置の出力信号を記憶する記憶装置と、
標準偏差を用いて計測値の検定を行い計測値が異
常であれば前記選択装置によつて選択されたシス
テム変数について点検を実施する異常検出装置
と、その点検結果を表示する表示装置とを備え
て、異常事態の発生と、その推定原因の調査結果
を知らしめようとするものである。
理解を容易にするため、自己回帰モデルについ
て説明する。
て説明する。
下水処理装置において、現在のプロセスの状態
は過去のプロセス状態の線形結合によつてその大
部分を表現できる。いま時刻ｎにおけるプロセス
の状態をｋ次元の全変数ベクトルＸ（ｎ）で表わ
すと、その自己回帰表現は次のようになる。
は過去のプロセス状態の線形結合によつてその大
部分を表現できる。いま時刻ｎにおけるプロセス
の状態をｋ次元の全変数ベクトルＸ（ｎ）で表わ
すと、その自己回帰表現は次のようになる。
Ｘ（ｎ）＝M
〓m=1
Ａ（ｍ）・Ｘ（ｎ−ｍ）＋Ｕ（ｎ） ……(1)
但しＸ（ｎ−ｍ）は時刻ｎよりｍ時点前の全変
数ベクトル、Ｕ（ｎ）は白色雑音ベクトル、Ａ
（ｍ）は自己回帰モデルの回帰係数、Ｍは自己回
帰モデルの最適次数である。
数ベクトル、Ｕ（ｎ）は白色雑音ベクトル、Ａ
（ｍ）は自己回帰モデルの回帰係数、Ｍは自己回
帰モデルの最適次数である。
(1)式の回帰係数Ａ（ｍ）の要素ij（ｍ）は次の連
立一次方程式の解として求められる。M 〓m=1 K 〓j=1 aij(m)・Rjh(l-m)＝Rih(l) ……(2) （ｈ＝１、２…、ｋ：１＝１、２、……Ｍ） ただし、Rih（ｌ）はＸのｉ、ｈ行の要素の相
互相関関数である。また白色雑音ベクトルＵ（ｎ）
の要素をｉ（ｎ）とすると、その分散σ2〓iはつぎの
ようになる。
立一次方程式の解として求められる。M 〓m=1 K 〓j=1 aij(m)・Rjh(l-m)＝Rih(l) ……(2) （ｈ＝１、２…、ｋ：１＝１、２、……Ｍ） ただし、Rih（ｌ）はＸのｉ、ｈ行の要素の相
互相関関数である。また白色雑音ベクトルＵ（ｎ）
の要素をｉ（ｎ）とすると、その分散σ2〓iはつぎの
ようになる。
σ2〓i＝Rii（ｏ）−M
〓m=1 K
〓j=1
aij（ｍ）・Rij（ｍ） ……(3)
なおモデルの最適次数Ｍは予測誤差を示す(4)式
のMFPE（Ｍ）を最小にする値である。
のMFPE（Ｍ）を最小にする値である。
MFPE（Ｍ）＝｛１＋Ｍ・ｋ＋１／Ｎ｝k・
｛１−Ｍ・ｋ＋１／Ｎ｝-k・‖dM‖ ……(4)
ただし、Ｎはデータ数、‖dM‖はＵ（ｎ）の分
数共分散行列推定値である。またMFPEは多次
元最終予測誤差Multiple Final Prediction
Errorの頭文字である。
数共分散行列推定値である。またMFPEは多次
元最終予測誤差Multiple Final Prediction
Errorの頭文字である。
このようにして自己回帰係数、モデル次数およ
び白色雑音の分散（固有ノイズ）が求められ、自
己回帰モデルが作成される。
び白色雑音の分散（固有ノイズ）が求められ、自
己回帰モデルが作成される。
つぎに、下水処理装置の被制御変数すなわち放
流水有機物濃度、放流水懸濁物質濃度、MLSSを
一定に保つために、この自己回帰モデルを構成す
る多数のシステム変数の中で、被制御変数に寄与
するシステム変数を選別する必要がある。
流水有機物濃度、放流水懸濁物質濃度、MLSSを
一定に保つために、この自己回帰モデルを構成す
る多数のシステム変数の中で、被制御変数に寄与
するシステム変数を選別する必要がある。
いまＫ個のシステム変数が伝達要素によつて結
ばれた閉レープにおいて、aij(f)を変数Xi(f)と変
数Xj(f)とを結ぶ伝達要素のxj(f)からxi(f)への周
波数応答関数とし、Ui(f)をxi(f)の内部雑音の周
波数領域における表現とすると、 xi(f)＝K 〓j=1 aij(f)・xj(f)＋Ui(f) ……(5) 但しｊ≠１ がえられる。
ばれた閉レープにおいて、aij(f)を変数Xi(f)と変
数Xj(f)とを結ぶ伝達要素のxj(f)からxi(f)への周
波数応答関数とし、Ui(f)をxi(f)の内部雑音の周
波数領域における表現とすると、 xi(f)＝K 〓j=1 aij(f)・xj(f)＋Ui(f) ……(5) 但しｊ≠１ がえられる。
ここで
〔bij(f)〕＝b11(f)
b21(f)
〓
〓
bk1(f) b12(f)
b22(f)
〓
〓
bk2(f) …
…
…
…
b1k(f)
b2k(f)
〓
〓
bkk(f)＝１−a11(f)
−a21(f)
〓
〓
−ak1(f)−a12(f)
−a22(f)
〓
〓
−ak2(f) …
…
…
…
−a1k(f)
−a2k(f)
〓
〓
１−akk(f) ……(6)
とおけば
ｉ(f)＝K
〓j=1
bij(f)・Uj(f) ……(7)
ここでbij(f)はｊ番目の変数の固有ノイズUj(f)
が、フイードバツクループを通してｉ番目システ
ム変数ｉ(f)へ及ぼす影響を示す。
が、フイードバツクループを通してｉ番目システ
ム変数ｉ(f)へ及ぼす影響を示す。
(7)式のパワースペクトル領域における表現式は
Pii(f)＝K
〓j=1
｜bij(f)｜2・Ｐ（Uj）(f) ……(8)
ただしＰ（Uj）(f)は周波数ｆにおける固有ノイ
ズUjのパワースペクトル密度である。さらに周
波数ｆにおけるxi(f)のパワースペクトル密度のう
ちUj(f)に寄与する部分をqij(f)とすると、 qij(f)＝｜ij(f)｜2P（Uj）(f) ……(9) となる。
ズUjのパワースペクトル密度である。さらに周
波数ｆにおけるxi(f)のパワースペクトル密度のう
ちUj(f)に寄与する部分をqij(f)とすると、 qij(f)＝｜ij(f)｜2P（Uj）(f) ……(9) となる。
本出願では、このqij(f)を寄与率と称する。
つまり、例えばあるシステム変数Ａに対するシ
ステム変数Ｂ、Ｃ、Ｄの寄与率が、それぞれ35
％、40％、25％と求まると、システム変数Ｃ、
Ｂ、Ｄの順で、システム変数Ａに影響を与えてい
ることになる。
ステム変数Ｂ、Ｃ、Ｄの寄与率が、それぞれ35
％、40％、25％と求まると、システム変数Ｃ、
Ｂ、Ｄの順で、システム変数Ａに影響を与えてい
ることになる。
以上が、自己回帰モデルの説明である。
また計測値の変動域は、標準偏差を求めること
によつて容易に推定できる。
によつて容易に推定できる。
すなわち標準偏差Ｓは(10)式で得られ、
ただし＝１／Ｎ・N
〓l=1
xi(l)
計測値は（11）式で得られるｉの範囲になると
考えて良い。
考えて良い。
−3S≦xi≦＋3S ……（11）
したがつて、刻々得られる計測値が、（11）式
の変動域に含まれているか否かによつて、正常、
異常の判定を下すことが可能となる。
の変動域に含まれているか否かによつて、正常、
異常の判定を下すことが可能となる。
以上が発明の原理である。
第２図は、上記した原理に基づいて成された本
発明の一実施例を示す構成図である。
発明の一実施例を示す構成図である。
第２図において１から９までは第１図の同一符
号と同じものを意味し、１０は流入下水量測定
計、１１は流入水のPH測定計、１２は流入水懸濁
物質濃度計、１３は流入水有機物濃度計、１４は
溶存酸素濃度計測計、１５はMLSS濃度計、１６
は水温計、１７は曝気槽流出水有機物濃度計、１
８は曝気槽流出水MLSS計、１９は放流水懸濁物
質濃度計、２０は放流水有機物濃度計、２１は放
流水のPH測定計、２２は返送汚泥量計測計、２３
は返送汚泥濃度計、２４は余剰汚泥引抜量測定
計、２５は余剰汚泥濃度計、２６は曝気風量測定
計、２７は係数演算装置、２８は選択装置、２９
は各計測値の標準偏差を算出する標準偏差演算装
置、３０は前記選択装置によつて選択された変
数、および前記標準偏差演算装置の出力を記憶し
ておく記憶装置、３１は異常検出装置、３２は表
示装置である。
号と同じものを意味し、１０は流入下水量測定
計、１１は流入水のPH測定計、１２は流入水懸濁
物質濃度計、１３は流入水有機物濃度計、１４は
溶存酸素濃度計測計、１５はMLSS濃度計、１６
は水温計、１７は曝気槽流出水有機物濃度計、１
８は曝気槽流出水MLSS計、１９は放流水懸濁物
質濃度計、２０は放流水有機物濃度計、２１は放
流水のPH測定計、２２は返送汚泥量計測計、２３
は返送汚泥濃度計、２４は余剰汚泥引抜量測定
計、２５は余剰汚泥濃度計、２６は曝気風量測定
計、２７は係数演算装置、２８は選択装置、２９
は各計測値の標準偏差を算出する標準偏差演算装
置、３０は前記選択装置によつて選択された変
数、および前記標準偏差演算装置の出力を記憶し
ておく記憶装置、３１は異常検出装置、３２は表
示装置である。
本発明の最適な実施に当つては、まず当該処理
場の統計的モデルおよび計測値の標準偏差を得る
ために、適切な期間（約１ケ月程度）、計測値を
採取し、統計的モデル、標準偏差を決定する。そ
してこれ以降は、刻々、放流水の水質に関わる重
要な水質項目について、通常の変動域内にあるか
否かを監視し、それを逸脱した時には、記憶装置
に記憶されている因果関係に従つて、その異常の
発生原因の調査を行う。
場の統計的モデルおよび計測値の標準偏差を得る
ために、適切な期間（約１ケ月程度）、計測値を
採取し、統計的モデル、標準偏差を決定する。そ
してこれ以降は、刻々、放流水の水質に関わる重
要な水質項目について、通常の変動域内にあるか
否かを監視し、それを逸脱した時には、記憶装置
に記憶されている因果関係に従つて、その異常の
発生原因の調査を行う。
管路１より曝気槽２に流入する下水の水量、
PH、流入水懸濁物質濃度、および流入水有機物濃
度は流入下水量測定計１０、PH測定計１１、懸濁
物質濃度計１２、有機物濃度計１３によりそれぞ
れ検出され、その検出信号は、演算装置２７に与
えられる。
PH、流入水懸濁物質濃度、および流入水有機物濃
度は流入下水量測定計１０、PH測定計１１、懸濁
物質濃度計１２、有機物濃度計１３によりそれぞ
れ検出され、その検出信号は、演算装置２７に与
えられる。
また、曝気槽２内に流入した下水と沈殿槽３か
ら返送された返送汚泥との混合水の溶存酸素濃
度、MLSSおよび水温は溶存酸素濃度測定計１
４、MLSS計１５、水温計１６によりそれぞれ検
出され、その検出信号は係数演算装置２７に与え
られる。また曝気槽２内で混合水が活性汚泥処理
された処理水の有機物濃度およびMLSSは、有機
物濃度計１７およびMLSS計１８によりそれぞれ
検出され、その検出信号は係数演算装置２７に与
えられる。また処理水が沈殿槽３で沈降分離され
清澄された放流水は管路７より装置外に放流さ
れ、その放流水の懸濁物質濃度、有機物濃度およ
びPHが懸濁物質濃度計１９、有機物濃度計２０お
よびPH測定計２１でそれぞれ検出され、係数演算
装置２７に与えられる。
ら返送された返送汚泥との混合水の溶存酸素濃
度、MLSSおよび水温は溶存酸素濃度測定計１
４、MLSS計１５、水温計１６によりそれぞれ検
出され、その検出信号は係数演算装置２７に与え
られる。また曝気槽２内で混合水が活性汚泥処理
された処理水の有機物濃度およびMLSSは、有機
物濃度計１７およびMLSS計１８によりそれぞれ
検出され、その検出信号は係数演算装置２７に与
えられる。また処理水が沈殿槽３で沈降分離され
清澄された放流水は管路７より装置外に放流さ
れ、その放流水の懸濁物質濃度、有機物濃度およ
びPHが懸濁物質濃度計１９、有機物濃度計２０お
よびPH測定計２１でそれぞれ検出され、係数演算
装置２７に与えられる。
さらに、返送汚泥量、返送汚泥濃度、余剰汚泥
量、余剰汚泥濃度、および曝気風量が返送汚泥量
測定計２２、返送汚泥濃度計２３、余剰汚泥引抜
量測定計２４、余剰汚泥濃度計２５および風量計
２６によつて検出され、係数演算装置２７に与え
られる。
量、余剰汚泥濃度、および曝気風量が返送汚泥量
測定計２２、返送汚泥濃度計２３、余剰汚泥引抜
量測定計２４、余剰汚泥濃度計２５および風量計
２６によつて検出され、係数演算装置２７に与え
られる。
このようにして、係数演算装置２７の自己回帰
モデル演算装置２８では、入力したそれぞれの時
系列の検出信号から、前記(2)、(3)、(4)式に基づい
たモデルの最適次数Ｍ、モデルの回帰係数Ａ（ｍ）
及び白色雑音の分散を演算すると共に、(1)式の自
己回帰モデルを演算し、選択装置２９では自己回
帰モデル演算装置２８からの出力により(9)式、(10)
式に基づいて、被制御変数すなわち放流水の有機
物濃度、懸濁物質濃度、および曝気槽内のMLSS
濃度に対するシステム変数の影響度合を演算し、
その因果関係を記憶装置３０へ記憶させておく。
標準偏差演算装置２９は、前記係数演算装置２７
に入力した測定値信号より(10)式から標準偏差を算
出し、記憶装置３０へ出力し、記憶させておく。
モデル演算装置２８では、入力したそれぞれの時
系列の検出信号から、前記(2)、(3)、(4)式に基づい
たモデルの最適次数Ｍ、モデルの回帰係数Ａ（ｍ）
及び白色雑音の分散を演算すると共に、(1)式の自
己回帰モデルを演算し、選択装置２９では自己回
帰モデル演算装置２８からの出力により(9)式、(10)
式に基づいて、被制御変数すなわち放流水の有機
物濃度、懸濁物質濃度、および曝気槽内のMLSS
濃度に対するシステム変数の影響度合を演算し、
その因果関係を記憶装置３０へ記憶させておく。
標準偏差演算装置２９は、前記係数演算装置２７
に入力した測定値信号より(10)式から標準偏差を算
出し、記憶装置３０へ出力し、記憶させておく。
このようにして、該処理場の自己回帰モデル及
び標準偏差が得られた後は、順次、計測される放
流水質や処理場の運転において重要な水質項目
（放流水有機物濃度、放流水懸濁物質濃度、放流
水PHおよびMLSS等）について監視を実施でき
る。
び標準偏差が得られた後は、順次、計測される放
流水質や処理場の運転において重要な水質項目
（放流水有機物濃度、放流水懸濁物質濃度、放流
水PHおよびMLSS等）について監視を実施でき
る。
監視のために、異常検出装置３１は記憶装置３
０に記憶させていた水質項目の標準偏差を（11）
式に適用し、該水質項目の計測値と比較すること
により、その正常、異常の判定を行う。「異常」
の判定が出た場合、記憶装置３０に記憶されてい
る因果関係に基づき、異常の発生原因となつた水
質項目の調査を実施し、その結果を表示装置３２
へ出し、表示を行う。
０に記憶させていた水質項目の標準偏差を（11）
式に適用し、該水質項目の計測値と比較すること
により、その正常、異常の判定を行う。「異常」
の判定が出た場合、記憶装置３０に記憶されてい
る因果関係に基づき、異常の発生原因となつた水
質項目の調査を実施し、その結果を表示装置３２
へ出し、表示を行う。
第３図に、ある稼働中の下水処理場で実測した
データを説明する。
データを説明する。
第３図は、放流水有機物濃度に対する寄与率を
周波数帯ごとに表わしたもので、その平均値が右
端の帯グラフに示されている。
周波数帯ごとに表わしたもので、その平均値が右
端の帯グラフに示されている。
この例では、放流水有機物濃度に変動に対し、
影響の大きな順は、それ自身の29.6％は別とし
て、流入水有機物濃度16.9％、溶存酸素濃度16.8
％、流入水量16.4％の順でつづいている。
影響の大きな順は、それ自身の29.6％は別とし
て、流入水有機物濃度16.9％、溶存酸素濃度16.8
％、流入水量16.4％の順でつづいている。
表示装置３２は、これらの順で表示されるよう
にしてもよいし、パーセントをともなつて一挙に
表示してもよい。
にしてもよいし、パーセントをともなつて一挙に
表示してもよい。
ここでは前記自己回帰モデルを作成して、その
モデルに従つて、一定時間、状態監視を行うよう
にしたが、一定周期たとえば１ケ月ごとにデータ
を採取して自己回帰モデルおよび標準偏差を修正
するようにいてもよい。
モデルに従つて、一定時間、状態監視を行うよう
にしたが、一定周期たとえば１ケ月ごとにデータ
を採取して自己回帰モデルおよび標準偏差を修正
するようにいてもよい。
本発明は以上述べたように、当該処理システム
の特徴、季節的な変動をも考慮した上で、下水処
理システムの異常発生を、すばやく検知し、その
原因となつた水質項目の調査をも実施、表示する
ことができるので、下水処理システムの自動化、
省力化に大なる効果がある。
の特徴、季節的な変動をも考慮した上で、下水処
理システムの異常発生を、すばやく検知し、その
原因となつた水質項目の調査をも実施、表示する
ことができるので、下水処理システムの自動化、
省力化に大なる効果がある。
第１図は、公知の下水処理装置のブロツク図、
第２図は本発明の実施例のブロツク図、第３図は
実測データの一例を示す図である。 １……管路、２……曝気槽、３……沈殿池、４
……管路、５……散気管、６……ブロワー、７…
…管路、８……返送汚泥ポンプ、９……引抜汚泥
ポンプ、１０……流入下水量測定計、１１……流
入水のPH測定計、１２……流入水懸濁物質濃度
計、１３……流入水有機物濃度計、１４……溶存
酸素量測定計、１５……MLSS濃度計、１６，１
７……水温計、１８……MLSS濃度計、１９……
放流水懸濁物質濃度計、２０……放流水有機物濃
度計、２１……放流水のPH測定計、２２……返送
汚泥量測定計、２３……返送汚泥濃度計、２４…
…余剰汚泥引抜量測定計、２５……余剰汚泥濃度
計、２６……曝気風量計、２７……係数演算装
置、２８……選択装置、２９……標準偏差演算装
置、３０……記憶装置、３１……異常検出装置、
３２……表示装置。
第２図は本発明の実施例のブロツク図、第３図は
実測データの一例を示す図である。 １……管路、２……曝気槽、３……沈殿池、４
……管路、５……散気管、６……ブロワー、７…
…管路、８……返送汚泥ポンプ、９……引抜汚泥
ポンプ、１０……流入下水量測定計、１１……流
入水のPH測定計、１２……流入水懸濁物質濃度
計、１３……流入水有機物濃度計、１４……溶存
酸素量測定計、１５……MLSS濃度計、１６，１
７……水温計、１８……MLSS濃度計、１９……
放流水懸濁物質濃度計、２０……放流水有機物濃
度計、２１……放流水のPH測定計、２２……返送
汚泥量測定計、２３……返送汚泥濃度計、２４…
…余剰汚泥引抜量測定計、２５……余剰汚泥濃度
計、２６……曝気風量計、２７……係数演算装
置、２８……選択装置、２９……標準偏差演算装
置、３０……記憶装置、３１……異常検出装置、
３２……表示装置。
Claims (1)
- 【特許請求の範囲】 １ 活性汚泥処理法による下水処理システムの処
理状況を監視する装置において、 (a) 流入下水、曝気槽、沈殿池、放流水、返送汚
泥および引抜汚泥の各プロセス状態の時系列計
測値を入力し、それぞれのプロセスの状態を表
わす自己回帰モデルの係数を算出する係数演算
装置と、 (b) 前記係数演算装置の出力を入力し、あるプロ
セス計測値に影響を及ぼす他のプロセス計測値
を選別する選択装置と、 (c) 前記計測値から標準偏差を算出する標準偏差
演算装置と (d) 前記選択装置及び前記標準偏差記憶の出力を
記憶する記憶装置と (e) 前記プロセス計測値が、前記記憶装置に記憶
されている前記標準偏差の範囲を外れたことを
検出する異常検出装置と、 (f) 前記異常検出装置が前記標準偏差の範囲を外
れたことを検出したときに、前記記憶装置に記
憶されていた前記選択装置の出力を表示する表
示器 とを備えたことを特徴とする下水処理システムの
監視装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59118913A JPS60261593A (ja) | 1984-06-08 | 1984-06-08 | 下水処理システムの監視装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59118913A JPS60261593A (ja) | 1984-06-08 | 1984-06-08 | 下水処理システムの監視装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60261593A JPS60261593A (ja) | 1985-12-24 |
| JPH0356115B2 true JPH0356115B2 (ja) | 1991-08-27 |
Family
ID=14748269
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59118913A Granted JPS60261593A (ja) | 1984-06-08 | 1984-06-08 | 下水処理システムの監視装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60261593A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015534031A (ja) * | 2012-09-21 | 2015-11-26 | ローズマウント インコーポレイテッド | 通風圧およびプロセス変数を用いた火炎不安定性のモニタリングの方法、システム及び装置 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017225918A (ja) * | 2016-06-21 | 2017-12-28 | 株式会社日立製作所 | 水処理システム |
-
1984
- 1984-06-08 JP JP59118913A patent/JPS60261593A/ja active Granted
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015534031A (ja) * | 2012-09-21 | 2015-11-26 | ローズマウント インコーポレイテッド | 通風圧およびプロセス変数を用いた火炎不安定性のモニタリングの方法、システム及び装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60261593A (ja) | 1985-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6862190B2 (ja) | プロセス診断装置、プロセス診断方法及びプロセス診断システム | |
| CN103534658B (zh) | 工序监视诊断装置 | |
| DE69033328T2 (de) | Unterstützungsverfahren und -vorrichtung für den Betrieb einer Anlage | |
| US6808630B2 (en) | System and method for AI controlling waste-water treatment by neural network and back-propagation algorithm | |
| JP4468269B2 (ja) | プロセス監視装置及びその方法 | |
| JP6661426B2 (ja) | プロセス診断装置、プロセス診断方法及びコンピュータプログラム | |
| CN109205781A (zh) | 控制系统及控制方法 | |
| WO2023163054A1 (ja) | 情報処理方法、情報処理装置、プログラム及び排水処理システム | |
| CN121559862A (zh) | 一种基于物联网的农污处理控制方法 | |
| CN118878053A (zh) | 基于神经网络模型的工艺药剂预测投加系统和方法 | |
| CN119143209A (zh) | 基于物联网的矿井疏干水处理设备智能运维管理系统 | |
| Klanderman et al. | Fault isolation for a complex decentralized waste water treatment facility | |
| JPH0356115B2 (ja) | ||
| CN114355846A (zh) | 一种基于sbr仿真模型的造纸污水处理过程故障诊断方法 | |
| JPH0356116B2 (ja) | ||
| JP4088075B2 (ja) | 水処理プロセスハイブリッド水質計測装置及びこれを有する水処理システム | |
| JPS6182896A (ja) | 下水処理制御装置 | |
| JPH07219623A (ja) | 計測器の異常検知装置 | |
| JP3866083B2 (ja) | 水処理プラント制御システム | |
| EP4453679B1 (en) | A method of monitoring processes in an urban water cycle | |
| Haimi | Data-derived soft sensors in biological wastewater treatment-With application of multivariate statistical methods | |
| JPH1133580A (ja) | プロセスの運転支援装置 | |
| JP2000010627A (ja) | 計測器の異常値検出方法およびその装置 | |
| JPH07275882A (ja) | 下水処理制御装置 | |
| Kim et al. | Nonparametric method for changepoint detection using Empirical likelihood |