(57)【要約】【課題】目標処理水条件を満たす適切な生物反応槽の土木構造，曝気条件，計装，運転条件等を支援することが可能とすることにある。【解決手段】反応槽内の汚泥濃度を設定する手段と、汚泥濃度情報中の不活性固形物濃度を設定する手段と、汚泥濃度情報と前記不活性固形物濃度から微生物濃度を演算する手段と、汚泥濃度情報中の自家栄養菌濃度を演算する手段と、汚泥濃度情報中のリン蓄積菌濃度を演算する手段と、汚泥濃度情報中の他家栄養菌濃度を演算する手段と、処理プロセスの運転条件設定手段と、最終沈殿池の汚泥沈降指標を入力するＳＶＩ設定手段を有するように構成した。【選択図】  図１
2005034717.tif 000002
【特許請求の範囲】
【請求項１】
反応槽内の汚泥濃度情報を入力する手段と、前記汚泥濃度情報に基づいて微生物濃度情報及び不活性固形物濃度情報を演算する手段を有することを特徴とする情報演算処理装置。
【請求項２】
活性汚泥法により処理するプロセスをシミュレーションするシステムであって、流入水の水量と水質を設定する流入条件設定手段と、処理プロセスの寸法や生物反応槽の分割を設定する土木構造設定手段と、反応槽内の汚泥濃度を設定する手段と、汚泥濃度情報中の不活性固形物濃度を設定する手段と、汚泥濃度情報と前記不活性固形物濃度から微生物濃度を演算する手段と、汚泥濃度情報中の自家栄養菌濃度を演算する手段と、汚泥濃度情報中のリン蓄積菌濃度を演算する手段と、汚泥濃度情報中の他家栄養菌濃度を演算する手段と、処理プロセスの運転条件設定手段と、最終沈殿池の汚泥沈降指標を入力するＳＶＩ設定手段と、生物反応槽，最終沈殿池の処理水及び返送汚泥の水質と汚泥濃度を演算するモデル演算手段と、該モデル演算手段にて演算された水質データと汚泥濃度を格納するデータベースと、該データベースから水質データと汚泥濃度を取り出すデータ編集手段と、データの表示手段を具備したことを特徴とする水質情報演算処理装置。
【請求項３】
請求項１乃至２において、汚泥濃度情報と不活性固形物濃度情報との相対的な差に基づいて微生物濃度情報を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項４】
請求項１乃至２において、好気的固形物滞留時間を用いて自家栄養菌濃度を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項５】
請求項１乃至２において、汚泥濃度情報中のリン濃度と、不活性固形物中のリン濃度と、微生物中のリン濃度を用いてリン蓄積菌濃度を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項６】
請求項１乃至２において、微生物濃度と自家栄養菌濃度とリン蓄積菌濃度との差に基づいて他家栄養菌濃度を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項７】
請求項１乃至２において、生物反応を演算する生物モデル演算手段と、流体を演算する輸送モデル演算手段，曝気装置の溶存酸素供給能力を演算する溶存酸素モデル演算手段を有することを特徴とする水質情報演算処理装置。
【請求項８】
請求項７において、前記生物モデル演算手段は、汚泥濃度情報の時間変化を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項９】
請求項１乃至８において、前記輸送モデル演算手段は、前記ＳＶＩ設定手段にて設定されたＳＶＩを用いて最終沈殿池の汚泥沈降速度を表現する手段と、返送汚泥濃度と余剰汚泥濃度を演算する手段を有することを特徴とする水質情報演算処理装置。
【請求項１０】
請求項１，８乃至９において、前記生物モデル演算手段にて演算された汚泥濃度情報と、前記輸送モデル演算手段にて演算された返送汚泥濃度及び余剰汚泥濃度を用いて、単位時間当たりの反応槽汚泥質量，返送汚泥質量，余剰汚泥質量を演算する手段と、固形物滞留時間を演算する手段を有することを特徴とする水質情報演算処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、水質情報演算処理装置に関し、特に、活性汚泥プロセスの設計や運転を支援するために、水質のシミュレーションを行うのに好適な水質情報演算処理装置に関する。
【０００２】
【従来の技術】
下水や工場排水などの汚水は活性汚泥と呼ばれる微生物群によって浄化されており、この処理を活性汚泥プロセスと称している。現在稼動している下水処理場は活性汚泥プロセスが主流で、有機物の除去を主目的とした標準活性汚泥法を採用している。標準活性汚泥法は、流入下水中の有機物を生物反応槽で活性汚泥に摂取あるいは酸化分解させた後、後段に設けられた最終沈殿池で活性汚泥を沈降させて上澄み液を放流する。
【０００３】
近年、閉鎖性水域での窒素・リンの総量規制を受け、下水処理場ではこれらの栄養塩類を除去する高度処理方式の導入が進められている。下水中からの窒素やリンを除去する高度処理方式は、大別して物理化学的方法と生物学的方法がある。生物学的方法は既存施設である標準活性汚泥法を改造して構築可能であることから、下水処理場に導入されつつある。標準活性汚泥法は生物反応槽全体に空気を供給しているので、生物反応槽は常に溶存酸素が存在する好気状態になっている。これに対して、生物学的な窒素・リン除去機構を利用した下水高度処理とは、生物反応槽に溶存酸素の存在しない嫌気状態を作り出し、好気状態と組合せてリン除去，窒素除去を実現する方法である。
【０００４】
生物学的窒素除去は、好気状態での硝化工程と無酸素状態での脱窒工程から構成する。一般的に、硝化反応を行う細菌を硝化菌あるいは自家栄養菌（以下、自家栄養菌と称す）と、脱窒反応を行う細菌を脱窒菌あるいは他家栄養菌（以下、他家栄養菌と称す）と称している。生物学的リン除去は、嫌気状態でのリン放出工程と好気状態でのリン過剰摂取工程を利用して下水中のリンを除去する処理方法である。このようにリン過剰摂取を行う能力を有する細菌を一般的にリン蓄積菌と呼ばれている。
【０００５】
このような下水高度処理は、有機物，窒素・リン除去に関連する他家栄養菌，自家栄養菌，リン蓄積菌、の生息環境を適切に維持することによって性能が発揮される。また、活性汚泥中の各種微生物の生息環境は流入条件や運転条件に大きく左右され、標準活性汚泥法に比べると、処理施設の設計や維持管理は複雑である。
【０００６】
しかし、これまで有機物，窒素・リンの複雑な反応過程を算出し提示できる方法がなく、下水高度処理の設計や運転は経験と勘に依存していた。そのため、経験していない流入下水水質と量，生物反応槽構成，運転条件に対しては予測できない事態が発生し、その都度対策を講じるという問題があった。
【０００７】
一方、生物反応をモデル化し、数値シミュレーションによって活性汚泥プロセスの特性を評価する方法が提案されている。生物反応のモデルの例として、海外では国際水協会（ＩＷＡ）が発表した活性汚泥モデルＡＳＭ  Ｎｏ．１（１９８６），Ｎｏ．２（１９９５），Ｎｏ．２ｄ（１９９８），Ｎｏ．３（１９９９）が提案されている。現在、窒素，リン除去モデルとして活性汚泥モデルＮｏ．２（以下ＡＳＭ−Ｎｏ．２と称す）とＮｏ．２ｄ（以下ＡＳＭ−Ｎｏ．２ｄと称す）が最も利用されている。しかし、これらの活性汚泥モデルは、生物反応を数学的に記述された構造であり、シミュレーションソフトウェアを提供するわけではない。また、特開２０００−１６７５８５号公報，特開２００２−１３７０号公報、のように、生物モデルを用いた下水処理プロセスシミュレータによって水質を計算する方法が提案されている。
【０００８】
【特許文献１】
特開２０００−１６７５８５号公報
【特許文献２】
特開２００２−１３７０号公報
【０００９】
【発明が解決しようとする課題】
従来技術のなかで国際水協会（ＩＷＡ）が発表した活性汚泥モデルＡＳＭ−Ｎｏ．２，ＡＳＭ−Ｎｏ．２ｄでは、硝化，脱窒，リン放出，リン摂取，有機物分解，菌体増殖・自己分解などの反応プロセスの反応速度は、それぞれ、他家栄養菌，自家栄養菌，リン蓄積菌の濃度によって計算されている。そのために、前記の反応プロセスの反応速度を正確に計算するためには、活性汚泥中の他家栄養菌，自家栄養菌，リン蓄積菌の濃度を把握しなければならない。国際水協会の活性汚泥モデルでは、重クロム酸カリウムによるＣＯＤ（化学的酸素要求量；以下ＣＯＤｃｒと略す）や呼吸速度や反応速度とを用いて、他家栄養菌，自家栄養菌，リン蓄積菌の濃度を表している。単位はｇＣＯＤ／ｍ３  である。しかし、呼吸速度や反応速度を計測するには特殊な機器と技術を要するだけではなく、他の反応を抑制するための様々な前処理が必要であり、処理場では容易に実施できない問題点がある。なお、前記の前処理とは、例えば、他家栄養菌の呼吸速度を測定する際は、硝化反応の抑制処理である。
【００１０】
一方、下水処理場において、日常の汚泥管理は、反応槽内の汚泥濃度（ＭＬＳＳ）（単位：ｍｇ／Ｌ）によって行われている。その理由は、手分析または自動計測器によって容易にＭＬＳＳを得ることができるからである。ＭＬＳＳの手分析測定は、反応槽の混合液を濾紙で濾過し、濾紙に残留した固形物を乾燥して、その固形物の重量を計測するという簡単な方法であるために、日常的に行われている。また、ＭＬＳＳの自動計測器の性能の向上に伴い、ＭＬＳＳ計による連続測定を行う処理場も多い。下水処理場において、例えば、ＳＲＴ（固形物滞留時間；単位：日）やＡ−ＳＲＴ（好気的固形物滞留時間；単位：日），ＳＶＩ（汚泥容量指標；単位：ｍＬ／ｇＳＳ）などの管理指標は全て、ＭＬＳＳを基準にしたものである。なお、ＳＲＴとは活性汚泥が余剰汚泥として引かれるまでの平均滞留時間をいい、ＳＶＩとは、ＳＶ３０（汚泥容量；単位：％）に基づき、汚泥１ｇが占める容量ｍＬ数で表したものをいう。前記ＳＶ３０［％］とは、生物反応槽内の混合液を１リットルのメスシリンダーに入れ、３０分静置したのち沈降した汚泥量をパーセントで表したものである。最終沈殿池での汚泥沈降状況はＳＶ３０、またはＳＶＩによって管理されている。ＳＶＩが大きいほど汚泥の沈降性が悪く、沈降しない汚泥は処理水と一緒に系外へ流出される。ＳＶＩは重要な指標であり、その定義を式（１）に示す。
【００１１】
2005034717.tif 000003
しかし、国際水協会の活性汚泥モデルのように、他家栄養菌，自家栄養菌，リン蓄積菌の濃度をｇＣＯＤ／ｍ３  で表すと、ＭＬＳＳ［ｍｇ／Ｌ］を計算できないため管理指標であるＳＲＴ，Ａ−ＳＲＴ，ＳＶＩを計算提示することができない。
【００１２】
また、低水温時、自家栄養菌の系内保持時間を高めるにはＡ−ＳＲＴを大きくする必要があり、そのためには、ＭＬＳＳを高める必要がある。一方、ＭＬＳＳを高めると空気量が増大するため、エネルギー節減の観点からは低く保つことが望ましい。このような相反する条件を適切に管理するにはＡ−ＳＲＴを維持できる最低ＭＬＳＳを求める必要がある。さらに、ＭＬＳＳは処理プロセスの維持管理に用いられるだけではなく、反応タンク容量の設計因子としても使用されている。このように、活性汚泥プロセスにとっては、ＭＬＳＳはなくてはならない必須の指標である。このように、活性汚泥をＭＬＳＳで表現することは非常に重要であるにもかかわらず、国際水協会の活性汚泥モデルでは表現できないという課題があった。
【００１３】
また、これらのことについて、従来の下水処理シミュレータでは、国際水協会の活性汚泥モデルと同様に微生物濃度をｇＣＯＤ／ｍ３  で扱うために、ＭＬＳＳを表現することはできなく、上記の課題を解決することができない。
【００１４】
そこで、本発明の目的は上記従来技術の問題点を克服し、ＭＬＳＳを基準に汚泥中の微生物濃度と不活性固形物濃度を定義し、生物モデルに反映し、目標処理水条件を満たす適切な生物反応槽の土木構造，曝気条件，計装，運転条件を支援することにある。
【００１５】
【課題を解決するための手段】
上記も目的を達成するために、本発明は、反応槽内の汚泥濃度情報を入力する手段と、汚泥濃度情報に基づいて微生物濃度情報及び不活性固形物濃度情報を演算する手段を有するように構成した。
【００１６】
或いは、活性汚泥プロセスの設計や運転を支援するために、水質シミュレーションを行う活性汚泥プロセスのシミュレーションシステムであって、流入水の水量と水質を設定する流入条件設定手段と、処理プロセスの寸法や生物反応槽の分割を設定する土木構造設定手段と、ＭＬＳＳを設定する手段と、該ＭＬＳＳ中の不活性固形物濃度を設定する手段と、ＭＬＳＳと前記不活性固形物濃度から微生物濃度を演算する手段と、ＭＬＳＳ中の自家栄養菌濃度を演算する手段と、ＭＬＳＳ中のリン蓄積菌濃度を演算する手段と、ＭＬＳＳ中の他家栄養菌濃度を演算する手段と、処理プロセスの運転条件設定手段と、最終沈殿池の汚泥沈降指標を入力するＳＶＩ設定手段と、生物反応槽，最終沈殿池の処理水及び返送汚泥の水質と汚泥濃度を演算するモデル演算手段と、該モデル演算手段にて演算された水質データと汚泥濃度を格納するデータベースと、該データベースから水質データと汚泥濃度を取り出すデータ編集手段と、データの表示手段を具備したことを特徴とする。
【００１７】
また、前記微生物濃度演算手段は、ＭＬＳＳと不活性固形物濃度との差に基づいて微生物濃度を演算する手段を有することを特徴とする。
【００１８】
また、前記自家栄養菌濃度演算手段は、Ａ−ＳＲＴを用いて自家栄養菌濃度を演算する手段を有することを特徴とする。
【００１９】
また、前記リン蓄積菌濃度演算手段は、ＭＬＳＳ中のリン濃度と、不活性固形物中のリン濃度と、微生物中のリン濃度を用いてリン蓄積菌濃度を演算する手段を有することを特徴とする。
【００２０】
また、前記他家栄養菌濃度演算手段は、微生物濃度と自家栄養菌濃度とリン蓄積菌濃度との差に基づいて他家栄養菌濃度を演算する手段を有することを特徴とする。
【００２１】
また、前記モデル演算手段は、生物反応を演算する生物モデル演算手段と、流体を演算する輸送モデル演算手段，曝気装置の溶存酸素供給能力を演算する溶存酸素モデル演算手段を有することを特徴とする。
【００２２】
また、前記生物モデル演算手段は、ＭＬＳＳの時間変化を演算する手段を有することを特徴とする。
【００２３】
また、前記輸送モデル演算手段は、前記ＳＶＩ設定手段にて設定されたＳＶＩを用いて最終沈殿池の汚泥沈降速度を表現する手段と、返送汚泥濃度と余剰汚泥濃度を演算する手段を有することを特徴とする。
【００２４】
また、前記データ編集手段は、前記生物モデル演算手段にて演算されたＭＬＳＳと、前記輸送モデル演算手段にて演算された返送汚泥濃度及び余剰汚泥濃度を用いて、単位時間当たりの反応槽汚泥質量，返送汚泥質量，余剰汚泥質量を演算する手段と、ＳＲＴを演算する手段を有することを特徴とする。
【００２５】
上記構成の効果に言及すると、混合液の浮遊物質は生物反応を行う微生物と生物反応に関与しない不活性固形物からなる。そのために、反応槽の活性汚泥濃度を表すＭＬＳＳを用いて他家栄養菌，自家栄養菌，リン蓄積菌濃度を定義することによって、処理水質を計算できるだけではなく、反応後のＭＬＳＳ，返送汚泥濃度と余剰汚泥濃度も計算可能となる。また、返送汚泥濃度と余剰汚泥濃度を用いて、単位時間当たりの返送汚泥質量と余剰汚泥質量を計算することもできる。さらに、ＭＬＳＳ，反応槽容積と単位時間当たり余剰汚泥質量を用いて、ＳＲＴやＡ−ＳＲＴを算出することができる。
【００２６】
また、流入水量，流入水質，季節，運転条件変動時、処理水質を満足できるＭＬＳＳを提示し、プロセスの維持管理を支援することができる。
【００２７】
このように、従来技術のように放流水の水質を計算しただけではなく、本発明のようにＭＬＳＳを計算することによって、土木，機械，計装，汚泥状況の整合を提示でき、汚濁物質の除去性能を向上させることができ、最適設計と最適運転を容易に行える。
【００２８】
【発明の実施の形態】
図１は本発明を活性汚泥プロセスのシミュレーションに適用した一実施例である。
【００２９】
図１において、嫌気槽１ａ，無酸素槽１ｂ，好気槽１ｃの３槽の生物反応槽１と、好気槽１ｃから無酸素槽１ｂに循環する嫌気−無酸素−好気法を説明する。流入汚水は沈砂池（図示せず）にて土砂，ゴミなど大きな固形物を沈降除去したのち、最初沈殿池（図示せず）に流入する。最初沈殿池では固形物を沈降除去し、有機物，アンモニア性窒素，リンなどを含む上澄み液は生物反応槽１への流入水３として送られる。生物反応槽１には最初沈殿池からの流入水３と返送汚泥管５からの返送汚泥（活性汚泥）が流入し、撹拌混合が行われる。一方、生物反応槽１には、送気管１２，送気装置１３を介してブロワ１１から空気が送気される。さらに、生物反応槽１では、循環ポンプ８によって循環汚泥管９を介して好気槽１ｃから無酸素槽１ｂへ汚泥を循環する。ここで、嫌気槽１ａには溶存酸素（ＤＯ）と硝酸性窒素（ＮＯ３  ）がともに存在しない状態であり、主にりん（ＰＯ４  −Ｐ）の放出反応が進行する。
【００３０】
嫌気槽１ａにおいて活性汚泥は体内に貯蔵していたりんを加水分解して下水中に放出する。又、活性汚泥はりん放出と同時に有機物を吸着し菌体内に貯蔵する。この生物反応によって、嫌気槽１ａではリンが増加し、有機物が減少する。
【００３１】
無酸素槽１ｂには嫌気槽１ａから流出した下水と返送汚泥の混合水、及び好気槽１ｃからの循環水が流入し、機械撹拌設備（図示せず）で撹拌混合される。無酸素槽１ｂでは好気槽１ｃから硝酸性窒素（ＮＯ３  ）が循環され、かつ溶存酸素の存在しない環境になり、硝酸性窒素は還元されて窒素ガス（Ｎ２  ）として大気中に放出される。これを脱窒反応と呼んでいる。
【００３２】
好気槽１ｃの底部には送気装置（散気装置）１３が設置されており、ブロワ１１から送気管（空気管）１２を介して送気された空気が送気装置（散気装置）１３によって散気し、好気槽１ｃ内の下水と活性汚泥からなる混合液を撹拌するとともに酸素を供給する。送気管（空気管）１２に設けられる開閉弁１４は好気槽への風量を調節する。ブロワ１１は、事前に設定した風量を送風する制御方式や、好気槽１ｃに設けられた溶存酸素濃度（ＤＯ）計１６の値が所定値に維持するように送風する制御方式によって運転される。好気槽１ｃ内の混合液の汚濁物質は、酸素供給により活発化した活性汚泥の働きにより処理される。例えば、活性汚泥は有機物を吸着し、供給された空気中の酸素を吸収して有機物を酸化分解して炭酸ガスと水にする。また、りんは好気条件において、活性汚泥中の菌体内に貯蔵され、嫌気槽１ａで放出された以上に摂取されて流入下水中の濃度よりも減少する。また、アンモニア性窒素は硝酸性窒素に酸化される。これを硝化反応と呼んでいる。なお、これら有機物，リン，アンモニア性窒素などの汚濁物質の一部は活性汚泥の増殖にも利用される。ＭＬＳＳ計１７は、生物反応槽内の混合液濃度を計測し、汚泥管理に使われる。
【００３３】
このような生物反応後の処理水は最終沈殿池２に導かれる。最終沈殿池２では、活性汚泥を重力沈降させ、上澄み液を塩素消毒した後、放流管１０によって放流する。最終沈殿池２の沈降汚泥の一部は返送ポンプ４によって返送汚泥管５を介して生物反応槽１へ送られ、残りの汚泥は余剰ポンプ６によって余剰汚泥管７を介して系外へ排出される。返送ポンプ４は、返送汚泥量の制御や、返送汚泥量と流入下水量の比率の制御によって運転される。循環ポンプ８は、循環液量の制御や、循環液量と流入下水量の比率の制御によって運転される。余剰ポンプ６は、余剰汚泥量の制御や、余剰泥量と流入下水量の比率の制御によって運転される。
【００３４】
以上に説明した嫌気−無酸素−好気法を対象としたシミュレータ２０の構成について、引き続き図１にて説明する。
【００３５】
シミュレータ２０はデータ設定装置３０，微生物濃度演算装置４０（ここで、微生物濃度情報とは、微生物濃度を含む、微生物濃度と実質的に同値な情報を示す。）と、モデル演算装置５０，データ編集手段６０，係数設定装置６５，入出力装置７０，データベース９０、から構成される。ここでは、プラントの設計に適用した例について説明する。データ設定装置３０はシミュレーションに必要なデータをキーボード７１、またはマウス７２を用いて入力し、モニタ７３に表示される。流入条件設定手段３１は流入汚水量，流入水質の濃度および水温の設定を行う。ここで水質とは例えば、有機物（易分解性と難分解性），アンモニア性窒素，全窒素，リン，浮遊物濃度，アルカリ度，溶存酸素，硝酸性窒素などである。データはプラント入力手段８０によって入力されたプラントの２４時間実測データでもよいし、日平均値と２４時間変動パターンにて作成した２４時間データでもよいし、２４時間を通して一定値としてもよい。土木構造設定手段３２は生物反応槽及び最終沈殿池の有効幅，有効長さおよび有効水深の寸法データ設定，生物反応槽の槽分割設定および配管設定を行う。前記の配管設定は例えば本実施例の循環汚泥管９，送気管（空気管）１２，返送汚泥管５の設置状況を設定する。運転条件設定手段３３はブロワ１１から生物反応槽１への送風量，好気槽１ｃから無酸素槽１ｂへの循環汚泥量，沈殿池２から嫌気槽１ａへの返送汚泥量，余剰汚泥量などの運転条件を設定する。
【００３６】
ＭＬＳＳ初期値設定手段３４は手分析またはＭＬＳＳ計によって測定された生物反応槽ＭＬＳＳを設定する（ここで、汚泥濃度情報とは、汚泥濃度を含む、汚泥濃度と実質的に同値な情報を示す。）。
【００３７】
不活性固形物濃度設定手段３５は、ＭＬＳＳ中で生物反応に関与しない不活性固形物濃度の設定を行う（ここで、不活性固形物濃度情報とは、不活性固形物濃度を含む、不活性固形物濃度と実質的に同値な情報を示す。）。通常の活性汚泥プロセスでは、この不活性固形物は主に流入水中の無機物と微生物自己分解時の生成物であるが、凝集剤などの化学薬品を添加する活性汚泥プロセスの場合、化学反応によって生成する塩類も不活性固形物となる。微生物の自己分解時は約１０％の不活性固形物が生成すると言われているので、これを参考に不活性固形物濃度を定義してもよいし、実データの解析および実験から得た結果を設定してもよい。不活性固形物濃度設定手段３５の一実施例として、式（２）に、アルミニウムを含む凝集剤（例えば、硫酸アルミニウム，ＰＡＣ）を添加したプロセスにおける不活性固形物濃度を計算する実施例を示す。
【００３８】
2005034717.tif 000004
ただし、ｆＡｌ：ＳＳ当たりのＡｌ含有量、ｆｓｓ＿Ａｌ：Ａｌ当たりの固形物発生量
ＳＳ（汚泥固形物量）当たりのＡｌ含有量は、汚泥分析値を参照してもよいし、凝集剤の注入率から得た検量線で設定してもよい。また、Ａｌ当たりの固形物発生量ｆｓｓ＿Ａｌは、凝集剤の種類によって異なり、文献値を参照してもよいし、実験で得られた実測値でもよい。また、一般的に下水処理場で用いられている凝集剤は、アルミニウム（Ａｌ）や鉄（Ｆｅ）塩などがある。式（２）の実施例には、アルミニウム（Ａｌ）を含む凝集剤を用いた例を示したが、硫酸第一鉄溶液や塩化第二鉄溶液のような鉄を含む凝集剤も用いた場合も同様に設定できる。
【００３９】
ＳＶＩ設定手段３６は、手分析よって測定されたＳＶＩを設定する。
【００４０】
このようなデータ設定装置３０から設定されたシミュレーション条件はデータベース９０に格納される。また、モニタ７３に設定内容をグラフィックなどにより表示する。
【００４１】
微生物濃度演算装置４０はデータベース９０のシミュレーション条件に基づき、微生物濃度演算手段４１と、自家栄養菌濃度演算手段４２と、リン蓄積菌濃度演算手段４３と、他家栄養菌濃度演算手段４４を用いて、生物反応に関与する微生物濃度の初期値を演算し、その結果をデータベース９０に格納する。
【００４２】
微生物濃度演算手段４１は、ＭＬＳＳ初期値設定手段３４で設定したＭＬＳＳと、不活性固形物濃度設定手段３５で設定した不活性固形物濃度を用いて、ＭＬＳＳ中の生物反応を行う微生物濃度設定を行う。微生物濃度は式（３）によって得られる。
【００４３】
2005034717.tif 000005
自家栄養菌濃度演算手段４２は、例えば式（３）で得られた微生物濃度を用いて自家栄養菌濃度を式（４）で計算する。
【００４４】
2005034717.tif 000006
ただし、φ：微生物当たりの最大自家栄養菌量［ｇＳＳ／ｇＳＳ］
θＸＡ：硝化に必要なＡ−ＳＲＴ［日］
また、自家栄養菌の硝化に必要なＡ−ＳＲＴであるθＸＡは、下水道施設計画・設計指針に基づいて式（５）で計算してもよいし、実験的に決定してもよい。
【００４５】
θＸＡ＝δ×２０．６×ｅｘｐ（−０．０６２７×Ｔ）                  …（５）
ただし、Ｔ：水温［℃］
δ：流入水Ｔ−Ｎの変動に対する補正係数
次に、図２により自家栄養菌濃度演算手段４２の演算フローを説明する。
【００４６】
図２は、自家栄養菌濃度演算手段４２の演算フローの一実施例である。ステップＳ１では、流入条件設定手段３１によって設定された水温と式（５）を用いて、硝化に必要なＡ−ＳＲＴを計算する。ステップＳ２では、ＭＬＳＳ初期値設定手段３４で得られたＭＬＳＳと、土木構造設定手段３２で得られた反応槽容積と、運転条件設定手段３３で得られた余剰汚泥量に基づきＡ−ＳＲＴの実績値を演算する。ステップＳ３では、式（３）で得られた微生物濃度と、ステップＳ１で得られた硝化に必要なＡ−ＳＲＴと、ステップＳ２で得られたＡ−ＳＲＴ実績値を用いて、自家栄養菌濃度を式（４）で計算する。
【００４７】
ここで、図１に戻り、リン蓄積菌濃度演算手段４３を説明する。
【００４８】
リン蓄積菌濃度演算手段４３は、ＭＬＳＳに含まれているリン濃度と、例えば式（２）で得られた不活性固形物濃度と、例えば式（３）で得られた微生物濃度を用いてリン蓄積菌濃度設定を算出する。一例を式（６）に示す。
【００４９】
2005034717.tif 000007
ただし、ｆｐｐ：リン蓄積菌当たりのポリリン酸含有量［ｇＰ／ｇＳＳ］
また、式（６）のなか、Ｓｐ，Ｂｐ，Ｉｐはそれぞれ、式（７），（８），
（９）で定義できる。
【００５０】
2005034717.tif 000008
ただし、ｆｓｐ：ＭＬＳＳ当たりのＰ含有量
ｆｘｐ：微生物当たりのＰ含有量
ｆｉｐ：不活性固形物当たりのＰ含有量
ｆｓｐ，ｆｘｐ，ｆｉｐ，ｆｐｐは文献値を参照してもよいし、実験で得られた経験値を用いてもよい。もちろん、これらの値は活性汚泥プロセスによって異なる。
【００５１】
次に、図３２によりリン蓄積菌濃度演算手段４３の演算フローを説明する。
【００５２】
図３はリン蓄積菌濃度演算手段４３の演算フローの一実施例を示す図である。ステップＳ１１では、ＭＬＳＳ初期値設定手段３４で得られたＭＬＳＳを用いて、ＭＬＳＳに含まれているリン濃度を式（７）で計算する。ステップＳ１２では、式（３）で得られた微生物濃度を用い、微生物に含まれているリン濃度を式（８）で計算する。ステップＳ１３では、式（２）で得られた不活性固形物濃度を用い、不活性固形物に含まれているリン濃度を式（９）で計算する。ステップＳ１４では、ステップＳ１１で得られたＭＬＳＳに含まれているリン濃度と、ステップＳ１２で得られた微生物に含まれているリン濃度と、ステップＳ１３で得られた不活性固形物に含まれているリン濃度を用いて、リン蓄積菌濃度を式（６）で計算する。
【００５３】
ここで、再び図１に戻り、他家栄養菌濃度演算手段４４を説明する。
【００５４】
図１の他家栄養菌濃度演算手段４４は、例えば、式（３）で得られた微生物濃度と、式（４）で得られた自家栄養菌濃度と、式（６）で得られたリン蓄積菌濃度を用いて他家栄養菌濃度を算出する。一例を式（１０）に示す。
【００５５】
2005034717.tif 000009
このように、微生物濃度演算装置４０から算出されたシミュレーション条件はデータベース９０に格納される。また、モニタ７３に算出内容をグラフィックなどにより表示する。
【００５６】
モデル演算装置５０はデータベース９０のシミュレーション条件に基づき、生物モデル演算手段５１と、輸送モデル演算手段５２と、溶存酸素モデル演算手段５３を用いて、生物反応槽，最終沈殿池，返送汚泥、及び余剰汚泥の水質，汚泥濃度及び流量を計算し、その結果をデータベース９０に格納する。
【００５７】
生物モデル演算手段５１は、生物反応モデルによって変化する水質、及びＭＬＳＳを計算する。これらの生物反応モデルには国際水協会が発表したＡＳＭ−Ｎｏ．２，ＡＳＭ−Ｎｏ．２ｄなどの公知のモデルを適用してもよいし、化学反応式から作成したモデルや実験的に求めたモデルを適用してもよい。ここで、ＡＳＭ−Ｎｏ．２の他家栄養菌増殖の反応速度を式（１１）に示す。
【００５８】
【数１】
2005034717.tif 000010
【００５９】
ただし、ＳＯ２：溶存酸素［ｇＯ２／ｍ３］，ＳＦ  ：容易に分解する有機物［ｇＣＯＤ／ｍ３］，ＳＡ：発酵産物（酢酸）［ｇＣＯＤ／ｍ３］，ＳＮＨ４：アンモニア態窒素［ｇＮ／ｍ３］，ＳＮＯ３：硝酸態窒素［ｇＮ／ｍ３］，ＳＰＯ４：溶解性リン［ｇＰ／ｍ３］，ＳＡＬＫ：アルカリ度［ｇＣＯＤ／ｍ３］，ＸＨ：他家栄養菌［ｇＣＯＤ／ｍ３］，μＨ：最大増殖速度［１／ｄ］，ＫＯ２：ＳＯ２の飽和係数［ｇＯ２／ｍ３］，ＫＦ：ＳＦ  の飽和係数［ｇＣＯＤ／ｍ３］，ＫＮＨ４：ＳＮＨ４の飽和係数［ｇＮ／ｍ３］，ＫＰＯ４：ＳＰＯ４の飽和係数［ｇＰ／ｍ３］，ＫＡＬＫ：ＳＡＬＫの飽和係数［ｇＣＯＤ／ｍ３］
発明が解決しようとする課題のなかにも述べたように、国際水協会の活性汚泥モデルでは、微生物濃度をｇＣＯＤ／ｍ３  で扱うために、例えば、式（１０）で得られた他家栄養菌濃度［ｍｇ／Ｌ］は、そのままでは式（１１）のＸｈには適用できない。
【００６０】
次に、図４を用いて本発明の生物モデル演算手段５１の一実施例をを説明する。図４に本発明の微生物濃度［ｍｇ／Ｌ］をＡＳＭ−Ｎｏ．２に適用する実施例を示す。ステップＳ２１では、式（４）で得られた自家栄養菌濃度［ｍｇ／Ｌ］と、式（６）で得られたリン蓄積菌濃度［ｍｇ／Ｌ］と、式（１０）で得られた他家栄養菌濃度［ｍｇ／Ｌ］の単位をｍｇ／ＬからｇＣＯＤ／ｍ３に変換する。ｍｇ／Ｌの単位をｇＣＯＤ／ｍ３に変換する一例を式（１２）に示す。
【００６１】
ＳＭ［ｇＣＯＤ／ｍ３］＝Ｓ［ｍｇ／Ｌ］×θ［ｇＣＯＤ／ｇ］    …（１２）
ただし、ＳＭ，Ｓ：微生物濃度、θ：単位変換値
なお、単位変換値θは微生物によって異なり、文献値を参照してもよいし、各自の実験データを用いてもよい。
【００６２】
ステップ２２では、ステップ２１で得られた他家栄養菌濃度［ｇＣＯＤ／ｍ３］，自家栄養菌濃度［ｇＣＯＤ／ｍ３］，リン蓄積菌濃度［ｇＣＯＤ／ｍ３］をＡＳＭ−Ｎｏ．２に適用し、水質及び微生物濃度，慢分解性固形物濃度（単位：ｇＣＯＤ／ｍ３  ）の変化を計算する。前記の慢分解性固形物は、微生物の自己分解時に生成する固形物質である。
【００６３】
ステップ２３では、ステップ２２で得られた、反応後の他家栄養菌濃度［ｇＣＯＤ／ｍ３］と、反応後の自家栄養菌濃度と、反応後のリン蓄積菌濃度と、反応後の慢分解性固形物濃度の単位をｇＣＯＤ／ｍ３からｍｇ／Ｌに変換する。ｇＣＯＤ／ｍ３の単位をｍｇ／Ｌに変換する一例を式（１３）に示す。
【００６４】
Ｓ［ｍｇ／Ｌ］＝ＳＭ［ｇＣＯＤ／ｍ３］／θ［ｇＣＯＤ／ｇ］    …（１３）
ステップ２４では、式（１３）で得られた、反応後の他家栄養菌濃度［ｍｇ／Ｌ］と、反応後の自家栄養菌濃度［ｍｇ／Ｌ］と、反応後のリン蓄積菌濃度［ｍｇ／Ｌ］と、反応後の慢分解性固形物濃度［ｍｇ／Ｌ］を用いて、反応後のＭＬＳＳを式（１４）で計算する。
【００６５】
2005034717.tif 000011
このように、生物モデル演算手段５１によってＭＬＳＳの濃度変化は計算できる。また、生物モデル演算手段５１は生物反応槽に適用するだけではなく、最終沈殿池にも同様に適用することで、生物反応による最終沈殿池の汚泥濃度変化も計算できる。
【００６６】
ここで、再び、図１を用いて、輸送モデル演算手段５２を説明する。
【００６７】
図１の輸送モデル演算手段５２は、流入汚水量，返送汚泥量，余剰汚泥量、及び循環汚泥量の流体に基づいて、変化する水質、及び汚泥濃度を計算する。また、最終沈殿池では、微生物濃度演算装置（ＳＶＩ設定手段）４０で得られたＳＶＩを用いて、汚泥沈降速度を表現し、沈降汚泥の汚泥濃度の変化を計算する。
【００６８】
溶存酸素モデル演算手段５３は送風量から好気槽１ｃに供給される溶存酸素を、例えば総括酸素移動容量係数の概念を利用して計算する。
【００６９】
これらモデル演算装置５０の計算結果はデータベース９０に格納される。
【００７０】
係数設定装置６５は、例えば、微生物濃度演算装置４０に必要な係数データと、生物モデル演算手段５１に必要な係数データをキーボード７１、またはマウス７２を用いて入力し、モニタ７３に表示される。前記の微生物濃度演算装置４０に必要な係数データは、例えば、式（２）のｆｓｓ＿Ａｌ，式（６）のｆｐｐなどであり、前記の生物モデル演算手段５１に必要な係数データは、例えば、式（１１）のμＨ，ＫＯ２，ＫＮＨ４などである。
【００７１】
データ編集手段６０は、データベース９０に格納されたデータを棒グラフ，トレンドグラフ，計算結果一覧，除去率，物質収支などの形式で編集しモニタ７３に出力する。
【００７２】
次に、図５に本発明のシミュレーションの動作を示す一実施例のフローチャートである。ステップＳ３１では、流入条件設定手段３１によって流入水量，流入水質（有機物，アンモニア性窒素，リン，ＳＳ，アルカリ度など）の濃度および水温を設定する。ステップＳ３２では、土木構造設定手段３２によって生物反応槽及び最終沈殿池の寸法と、生物反応槽の槽分割と、配管設定を行う。ステップＳ３３では、運転条件設定手段３３によって、ブロワ１１の制御条件（送風量，ＤＯ目標），返送ポンプ４の制御条件（返送汚泥量，返送率の目標値，タイマー引抜き），余剰ポンプ６の制御条件（余剰汚泥量，余剰率の目標値，タイマー引抜き），循環ポンプ８の制御条件（循環汚泥量，循環率の目標値）を設定する。ステップＳ３４ではＭＬＳＳ初期値設定手段３４によってＭＬＳＳの初期値を設定する。ステップＳ３５では式（２）を用いて不活性固形物濃度を設定する。ステップＳ３６では、ＳＶＩ設定手段３６によってＳＶＩを設定する。ステップＳ３７では式（３）を用いて微生物濃度の初期値を設定する。ステップＳ３８では、式（４）を用いて自家栄養菌濃度初期値を演算する。ステップＳ３９では、式（６）を用いてリン蓄積菌濃度初期値を演算する。ステップＳ４０では、式（１０）を用いて他家栄養菌濃度初期値を設定する。ステップＳ４１では、ステップ３７，３８，３９で設定した自家栄養菌濃度，リン蓄積菌濃度，他家栄養菌濃度の単位を式（１２）によってｍｇ／ＬからｇＣＯＤ／ｍ３に変換する。
【００７３】
以上の条件設定に基づいて、ステップＳ４２ではシミュレーションを実行し、有機物，窒素，リン，微生物，慢分解性固形物の濃度を計算する。ステップＳ４３では、ステップＳ４２の計算結果のうち、自家栄養菌濃度，リン蓄積菌濃度，他家栄養菌濃度，慢分解性固形物濃度の単位を式（１３）によってｇＣＯＤ／ｍ３  からｍｇ／Ｌに変換する。ステップＳ４４では、式（１４）を用いてＭＬＳＳ，返送汚泥濃度，余剰汚泥濃度の変化を計算する。
【００７４】
図６に、図５の実施例のステップＳ３４，Ｓ３５，Ｓ３６の設定画面例を示す。
【００７５】
図７には、係数設定装置６５の実施例を示す。図７には、式（６）のリン蓄積菌濃度の演算に必要な係数データを設定する画面例を示す。図５の実施例ステップ３９では、実施例図７の係数データを用いて、式（６）によってリン蓄積菌濃度を計算する。
【００７６】
図８に、微生物濃度演算装置（ＳＶＩ設定手段）４０で入力したＳＶＩを用いて最終沈殿池の汚泥沈降速度を表す実施例を示す。ＳＶＩが小さいほど、汚泥の沈降速度が大きく、逆に、ＳＶＩが大きいほど、汚泥の沈降速度が小さくなる。このように、ＳＶＩを用いて汚泥の沈降速度を表現することで、最終沈殿池の沈降汚泥の汚泥濃度を表現することができ、返送汚泥濃度と余剰汚泥濃度を計算できる。
【００７７】
図９には、図５の実施例ステップ４４で計算された返送汚泥濃度の２４ｈトレンドグラフの実施例を示す。このように、常時に変動している汚泥濃度（微生物と不活性固形物を含む）を計算することができる。
【００７８】
図１０にデータ編集手段６０の一実施例を示す。図１０には、図５の実施例で得られ反応槽流出ＭＬＳＳ，返送汚泥濃度と余剰汚泥濃度を用いて、単位時間当たりの汚泥質量を計算した例を示す。単位時間当たりの汚泥質量は式（１５）によって得られる。
【００７９】
汚泥質量［ｋｇ／ｈ］＝汚泥濃度［ｍｇ／Ｌ］×流量［ｍ３／ｈ］    …（１５）
なお、式（１５）の実施例には１ｈ当たりの汚泥質量を示したが、１日当たりの汚泥質量も同様に計算できる。図１０の実施例に示すように、生物反応槽と最終沈殿池において、各々の汚泥収支を取ることで、ＳＲＴやＡ−ＳＲＴを計算することができ、ＳＲＴやＡ−ＳＲＴによる維持管理が可能となる。なお、ＳＲＴは式（１６）によって得られる。
【００８０】
2005034717.tif 000012
図１１に、本発明による維持管理の実施例を示す。図１１には、余剰汚泥量を減らしてシミュレーションをした結果の実施例を示す。余剰汚泥量を減らすと、ＭＬＳＳは経過時間とともに増加し、放流水Ｔ−ＮはＭＬＳＳの増加に伴い、減少する。図１１には、ｔ時間後、ＭＬＳＳ２３００ｍｇ／Ｌに達するとき、放流水Ｔ−Ｎは１０ｍｇ／Ｌを満足できる例を示す。本発明による維持管理は、例えば、図１１の結果に基づいて、反応槽内のＭＬＳＳを２３００ｍｇ／Ｌに維持するように余剰汚泥量を操作する。このように、流入水量，流入水質，運転条件や水温の変動時、図１１のように処理水質を満足できるＭＬＳＳを提示し、処理プロセスの維持管理に反映できる。
【００８１】
以上のシミュレーションを、プロセスの土木構造や計装設備，汚泥仕様，運転条件を変えて試行し、プロセス内の水質と汚泥の変動を把握できる。この結果、プロセス全体が最適な除去性能を発揮するために、どこの土木構造や運転条件を改善すべきか容易に把握できる。
【００８２】
このように、処理水質を計算できるだけではなく、反応後のＭＬＳＳ，返送汚泥濃度と余剰汚泥濃度も計算可能となる。また、返送汚泥濃度と余剰汚泥濃度を用いて、単位時間当たりの返送汚泥質量と余剰汚泥質量を計算することもできる。さらに、ＭＬＳＳ，反応槽容積と単位時間当たり余剰汚泥質量を用いて、ＳＲＴやＡ−ＳＲＴを算出することができる。また、流入水量，流入水質，季節，運転条件変動時、処理水質を満足できるＭＬＳＳを提示し、プロセスの維持管理を支援することができる。
【００８３】
したがって、土木，機械，計装，汚泥状況の整合を提示することで、汚泥濃度条件と運転の最適化を図れるので、流入条件や曝気や返送や余剰など運転条件と組み合わせてシミュレーションを実行し、有機物，窒素，リンの除去性能を容易に高めることができる。
【００８４】
【発明の効果】
以上説明したとおり、本発明によれば、目標処理水条件を満たす適切な生物反応槽の土木構造，曝気条件，計装，運転条件等を支援することが可能となる。
【図面の簡単な説明】
【図１】本発明の活性汚泥プロセスのシミュレーションシステムの一実施例を示す構成図。
【図２】本発明の自家栄養菌濃度を演算する手段を示す実施例。
【図３】本発明のリン蓄積菌濃度を演算する手段を示す実施例。
【図４】本発明の生物モデルを演算する手段を示す実施例。
【図５】本発明の一実施例によるシミュレーション手順を示すフロー図。
【図６】本発明のＭＬＳＳを設定する画面例。
【図７】本発明の係数設定装置を示す画面例。
【図８】本発明の汚泥沈降速度を示す実施例。
【図９】本発明による汚泥濃度の計算結果を示す実施例。
【図１０】本発明による汚泥質量の計算結果を示す実施例。
【図１１】本発明による維持管理を示す実施例。
【符号の説明】
１…生物反応槽、２…最終沈殿池、３…流入水、４…返送ポンプ、５…返送汚泥管、６…余剰ポンプ、７…余剰汚泥管、８…循環ポンプ、９…循環汚泥管、１０…放流管、１１…ブロワ、１２…送気管、１３…送気装置、１４…開閉弁、１６…溶存酸素濃度（ＤＯ）計、１７…ＭＬＳＳ計、２０…シミュレータ、３０…データ設定装置、３１…流入条件設定手段、３２…土木構造設定手段、３３…運転条件設定手段、３４…ＭＬＳＳ初期値設定手段、３５…不活性固形物濃度設定手段、３６…ＳＶＩ設定手段、４０…微生物濃度演算装置、４１…微生物濃度演算手段、４２…自家栄養菌濃度演算手段、４３…リン蓄積菌濃度演算手段、４４…他家栄養菌濃度演算手段、５０…モデル演算装置、５１…生物モデル演算手段、５２…輸送モデル演算手段、５３…溶存酸素モデル演算手段、６０…データ編集手段、６５…係数設定装置、７０…入出力装置、７１…キーボード、７２…マウス、７３…モニタ、８０…プラント入力手段、９０…データベース。
【図１】
2005034717.tif 000013
【図２】
2005034717.tif 000014
【図３】
2005034717.tif 000015
【図４】
2005034717.tif 000016
【図５】
2005034717.tif 000017
【図６】
2005034717.tif 000018
【図７】
2005034717.tif 000019
【図８】
2005034717.tif 000020
【図９】
2005034717.tif 000021
【図１０】
2005034717.tif 000022
【図１１】
2005034717.tif 000023