(57)【特許請求の範囲】
【請求項１】
処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を予測すべく、前記細菌による前記処理対象物質の最大反応速度の値をパラメータに用いるシミュレーション方法であって、
生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で前記最大反応速度の値をパラメータに用い、しかも、前記関数は、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴とするシミュレーション方法。
【請求項２】
前記処理対象物質が窒素成分であり、前記生物処理工程が硝化工程又は脱窒工程のいずれかである請求項１記載のシミュレーション方法。
【請求項３】
前記窒素成分がアンモニア性窒素であり、前記硝化工程におけるアンモニア酸化細菌の酸化反応の最大反応速度をＶAOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）、硝化工程においてアンモニア酸化細菌１個あたりに単位時間に負荷される前記アンモニア性窒素の量をＬAOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）としたときに、下記式（４）
【数１】
0004278701.tif 000048
〔ただし、ｆAOBは、１．０×１０2≦ＬAOB≦３．５×１０4において、｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）−２．５×１０3｝≦ＶAOB≦｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）＋２．５×１０3｝となるように定められる関数である。〕
となる関数ｆAOBでＶAOBとＬAOBとを関係付けてパラメータに用いる請求項２記載のシミュレーション方法。
【請求項４】
前記窒素成分が亜硝酸性窒素であり、前記硝化工程における亜硝酸酸化細菌の酸化反応の最大反応速度をＶNOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）、硝化工程において亜硝酸酸化細菌１個あたりに単位時間に負荷される前記亜硝酸性窒素の量をＬNOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）としたときに、下記式（５）
【数２】
0004278701.tif 000049
〔ただし、ｆNOBは、１．０×１０3≦ＬNOB≦１．２×１０6において、｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）−１．０×１０5｝≦ＶNOB≦｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）＋１．０×１０5｝となるように定められる関数である。〕
となる関数ｆNOBでＶNOBとＬNOBとを関係付けてパラメータに用いる請求項２記載のシミュレーション方法。
【請求項５】
前記窒素成分が硝酸性窒素であり、前記脱窒工程における硝酸還元細菌の還元反応の最大反応速度をＶNARB（ｆｇ・ｃｏｐｙ-1・ｈ-1）、脱窒工程において硝酸還元細菌１個あたりに単位時間に負荷される前記硝酸性窒素の量をＬNARB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）としたときに、下記式（６）
【数３】
0004278701.tif 000050
〔ただし、ｆNARBは、５．０≦ＬNARB≦２．５×１０3において、｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）−１．７×１０2｝≦ＶNARB≦｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）＋７０｝となるように定められる関数である。〕
となる関数ｆNARBでＶNARBとＬNARBとを関係付けてパラメータに用いる請求項２記載のシミュレーション方法。
【請求項６】
前記窒素成分が亜硝酸性窒素であり、前記脱窒工程における亜硝酸還元細菌の還元反応の最大反応速度をＶNIRB（ｆｇ・ｃｏｐｙ-1・ｈ-1）、脱窒工程において亜硝酸還元細菌１個あたりに単位時間に負荷される前記亜硝酸性窒素の量をＬNIRB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）としたときに、下記式（７）
【数４】
0004278701.tif 000051
〔ただし、ｆNIRBは、５．０≦ＬNIRB≦３．５×１０3において、｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）−２．５×１０2｝≦ＶNIRB≦｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）＋２．５×１０2｝となるように定められる関数である。〕
となる関数ｆNIRBでＶNIRBとＬNIRBとを関係付けてパラメータに用いる請求項２記載のシミュレーション方法。
【請求項７】
前記窒素成分とともに被処理水に塩化物イオンが含有されており、該被処理水を細菌によって生物学的に硝化処理又は脱窒処理した後の処理水の水質を、前記被処理水に含有されている窒素成分の濃度に基づいて計算された最大反応速度の値をパラメータに用い、しかも、前記被処理水中の塩化物イオン濃度を変数とする関数の計算結果に基づいて前記最大反応速度の値を計算することにより、前記塩化物イオン濃度との間に関数関係を有する状態で、前記最大反応速度の値をパラメータに用いて予測する請求項３乃至６のいずれか１項に記載のシミュレーション方法。
【請求項８】
処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を予測すべく、前記細菌による前記処理対象物質の最大反応速度の値がパラメータに用いられてシミュレーションが実施されるシミュレーション装置であって、
生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で前記最大反応速度の値がパラメータに用いられており、しかも、前記関数は、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴とするシミュレーション装置。
【請求項９】
処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を、前記細菌による前記処理対象物質の最大反応速度の値をパラメータに用いたシミュレーションによって予測しつつ前記生物処理工程を実施する生物処理方法であって、
生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で最大反応速度の値を前記パラメータに用いており、しかも、前記関数が、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴とする生物処理方法。
【請求項１０】
処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を、前記細菌による前記処理対象物質の最大反応速度の値がパラメータに用いられたシミュレーションによって予測しつつ前記生物処理工程が実施される生物処理装置であって、
生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で最大反応速度の値が前記パラメータに用いられており、しかも、前記関数が、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴とする生物処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シミュレーション方法、シミュレーション装置、生物処理方法、ならびに、生物処理装置に関する。
【背景技術】
【０００２】
現在、水処理、例えば、下水処理、工場排水処理などにおいて、処理の効率化、処理能力の高度化、水処理に用いられるエネルギーの省力化、水処理のコストの低コスト化などが進められている。
例えば、種々の条件下でのプロセスの挙動を経験に基づいて予測し、水処理施設の運転条件の設定を行なう手法は、種々の条件下でのプロセスの挙動を定量的に予測することが困難であるため、水処理が非効率的であり、エネルギーの浪費、コストの増加をもたらすという欠点があり改良が求められている。
このようなことから、経験に基づく種々の条件下でのプロセスの挙動の予測に代えて、細菌群の増殖や死滅などの反応を計算するシミュレーションが導入され、より定量的な予測を行なうことが試みられている（例えば、特許文献１〜４参照）。
【０００３】
この内、特許文献２には、下水処理プロセスを構成する単位装置を部品としてモデル化する、下水処理プロセスシミュレータシステムが開示されている。
このシステムは、一般に「ＡＳＭ２」などと呼ばれる、ＩＡＷＱ（現「ＩＷＡ」）活性汚泥モデルＮｏ．２を有し、該ＩＡＷＱ（現「ＩＷＡ」）活性汚泥モデルＮｏ．２によりモデル成分入力値に基づいて部品毎にモデル成分出力値を求めるシミュレータと、下水処理プロセスへ流入する流入水の水質をオンラインで計測する計測手段とを有している。
さらに、このシステムは、計測値とモデル成分入力値を含む相関式を用いた変換手段を有し、計測手段からの計測値を変換手段の相関式を用いてモデル成分入力値に変換する演算手段とを備えている。
そして、前記特許文献２に記載の発明は、流入水の水質と処理水質に基づき、反応速度論定数をキャリブレーションするものである。
【０００４】
しかし、生物処理における処理水の水質は、種々の要因に影響を受け、たとえ処理水質に基づくキャリブレーションを実施して、一旦は、実際の処理水質に近いシミュレーション結果が得られたとしても、時間の経過とともに、実際の処理水質とシミュレーション結果との乖離が大きくなる。
したがって、頻繁にキャリブレーションをやり直さなければ処理水質の予測精度を低下させてしまうこととなる。
すなわち、従来のシミュレーション方法やシミュレーション装置においては、予測精度の低下を抑制しつつキャリブレーションの手間を削減することが困難であるという問題を有している。
また、このようなシミュレーションによって処理水質の予測を実施しつつ生物処理工程を実施する生物処理方法や生物処理装置においては、予測精度の低下を防止することが困難なことから、実際の処理水質を確認する作業などの手間を削減することが困難である。
【０００５】
【特許文献１】特開平８−３２３３９３号公報
【特許文献２】特開２０００−１０７７９６号公報
【特許文献３】特公平７−１０６３５７号公報
【特許文献４】特開平９−４７７８５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
本発明は、高い予測精度を維持しつつキャリブレーションの手間を削減させ得るシミュレーション方法ならびにシミュレーション装置の提供を課題としている。
また、本発明は、要する手間を削減し得る生物処理方法や生物処理装置の提供を課題としている。
【課題を解決するための手段】
【０００７】
本発明者らは、前記課題を解決すべく鋭意検討を行った結果、従来は、定数として扱われていた細菌による処理対象物質の最大反応速度の値が、生物処理工程において単位数量の細菌に単位時間あたりに負荷される前記処理対象物質の量や単位数量の細菌に単位時間あたりに処理された前記処理対象物質の量によって変化することを見出した。
また、これらと関数関係を有する状態で最大反応速度の値をパラメータに用いることで従来よりも予測精度に優れたシミュレーション結果が得られることを見出し本発明の完成に到ったのである。
【０００８】
すなわち、本発明に係るシミュレーション方法は、処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を予測すべく、前記細菌による前記処理対象物質の最大反応速度の値をパラメータに用いるシミュレーション方法であって、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で前記最大反応速度の値をパラメータに用い、しかも、前記関数は、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴としている。
【００１０】
また、本発明にかかるシミュレーション装置は、処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を予測すべく、前記細菌による前記処理対象物質の最大反応速度の値がパラメータに用いられてシミュレーションが実施されるシミュレーション装置であって、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で前記最大反応速度の値がパラメータに用いられており、しかも、前記関数は、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴としている。
【００１２】
また、本発明にかかる生物処理方法は、処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を、前記細菌による前記処理対象物質の最大反応速度の値をパラメータに用いたシミュレーションによって予測しつつ前記生物処理工程を実施する生物処理方法であって、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で最大反応速度の値を前記パラメータに用いており、しかも、前記関数が、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴としている。
【００１４】
さらに、本発明にかかる生物処理装置は、処理対象物質を含有する被処理水を、前記処理対象物質を分解する細菌によって生物学的に処理する生物処理工程後の処理水の水質を、前記細菌による前記処理対象物質の最大反応速度の値がパラメータに用いられたシミュレーションによって予測しつつ前記生物処理工程が実施される生物処理装置であって、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量との間に関数関係を有する状態で最大反応速度の値が前記パラメータに用いられており、しかも、前記関数が、前記最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数であることを特徴としている。
【００１６】
なお、「最大反応速度の値をＶとし、生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量をＬとしたときに、Ｌの値の増大とともにＶの値を増大させる関数」とは、横軸にＬの値、縦軸にＶの値をとってこの関数をグラフ化したときに、全体が右上がりの状態になっていることを意図するものであり、局所的にＬの値の増大によってＶの値が低下する場合をも含むことを意図するものである。
【００１７】
また、この「生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量」とは、「生物処理工程に外部から導入される処理対象物質の量」ならびに「生物処理工程において他の細菌などによって産生される処理対象物質の量」などを合計した「処理対象物質の量」を意図している。
【００１８】
また、この「生物処理工程において細菌１個あたりに単位時間に処理された処理対象物質の量」とは、「生物処理工程に外部から導入される処理対象物質の量の内で、分解された処理対象物質の量」ならびに「生物処理工程において他の細菌などによって産生される処理対象物質の量の内で、分解された処理対象物の量」などを合計した「分解された処理対象物質の量」を意図している。
この「分解された処理対象物の処理量」は、「生物処理工程に外部から導入される処理対象物質の量」と「生物処理工程において他の細菌などによって産生される処理対象物質の量」との合計から「処理液に残留する処理対象物質の量」を減じることで求めることができる。
【００１９】
なお、本明細書中においては、特別の記載がない限りにおいては、「最大反応速度」との用語は、一個の細菌が単位時間あたりに分解する処理対象物質の最大質量を意図しており、例えば、この「最大反応速度」は、生物処理工程において単位時間あたりの処理対象物質の濃度変化の量：ΔＳ（ｆｇ／ｍ3／ｈ）と、生物処理工程において用いられている単位体積あたりの細菌数：ｎ（ｃｏｐｉｅｓ／ｍ3）とを求めて下記式（４）を計算して求めることができる。
【００２０】
【数１３】
0004278701.tif 000002
【００２１】
また、単位時間あたりに分解される処理対象物質の濃度変化量については、処理対象物質を含む溶液と処理対象物質を分解する細菌を含む汚泥とを混合した試料を作製して、該試料中の処理対象物質の濃度の経時変化を測定するなどして求めることができる。
【００２２】
また、「生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量」は、例えば、生物処理工程に単位時間あたりに単位体積あたりに外部から導入される処理対象物質の量と生物処理工程において他の細菌などによって単位時間あたりに単位体積あたりに産生される処理対象物質の量との合計量：Ｑ（ｆｇ／ｍ3／ｄａｙ）と、生物処理工程において用いられている単位体積あたりの細菌数：ｎ（ｃｏｐｉｅｓ／ｍ3）とを求めて下記式（５）を計算して求めることができる。
【００２３】
【数１４】
0004278701.tif 000003
【００２４】
また、「生物処理工程において細菌１個あたりに単位時間に処理された処理対象物質の量」も、上記式（５）における“Ｑ（ｆｇ／ｍ3／ｄａｙ）”を、生物処理工程に外部から導入される処理対象物質の量と生物処理工程において他の細菌などによって産生される処理対象物質の量との合計から処理液に残留する処理対象物質の量を減じたものに置き換えて上記式（５）における“Ｌ（ｆｇ／ｃｏｐｙ／ｄａｙ）”を計算することで求めることができる。
【００２５】
そして、「最大反応速度と生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量との間の関係を表す関数ｆ（ｘ）がｇ（ｘ）とｈ（ｘ）との間を通る関数である」とは、ｘ＞０において「ｈ（ｘ）＞ｆ（ｘ）＞ｇ（ｘ）の関係を満たしている」ことを意図している。
【００２６】
なお、ｈ（ｘ）＞ｆ（ｘ）＞ｇ（ｘ）の関係については、生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量（Ｌ：ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が実質上とり得る値の範囲において満足していればよく、例えば、ｘが無限大となるような場合においてもｈ（ｘ）＞ｆ（ｘ）＞ｇ（ｘ）の関係を満たしていることを必要とするものではない。
この生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量（Ｌ）が実質上とり得る値の範囲とは、処理対象物質が主としてアンモニアや硝酸等の窒素化合物で、その他に有機化合物等があまり含有されていない窒素含有排水などを処理するためのアンモニア酸化細菌によるアンモニア酸化の工程においては、通常、１００〜４０００（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
また、例えば、亜硝酸酸化細菌による亜硝酸酸化の工程においては、通常、１０００〜６００００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
また、硝酸還元細菌による硝酸還元の工程においては、通常、５〜７０（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
さらに、亜硝酸還元細菌による亜硝酸還元の工程においては、通常、５〜１２０（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
【００２７】
また、下水等、アンモニア等の窒素化合物以外に有機化合物も含有されている場合をも含めると、生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量（Ｌ）が実質上とり得る値の範囲は、アンモニア酸化細菌によるアンモニア酸化の工程においては、通常、１００〜３５０００（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
また、亜硝酸酸化細菌による亜硝酸酸化の工程においては、通常、１０００〜１２００００００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
また、硝酸還元細菌による硝酸還元の工程においては、通常、５〜２５００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
さらに、亜硝酸還元細菌による亜硝酸還元の工程においては、通常、５〜３５００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である。
【００２８】
なお、細菌１個あたりに単位時間に処理されている処理対象物質の量が、生物処理工程において細菌１個あたりに単位時間に負荷されている処理対象物質の量を超えることは通常の生物処理工程においては発生しない。
そして、硝化、脱窒といった生物処理工程において細菌１個あたりに単位時間に負荷される処理対象物質（窒素成分）の量が、アンモニア酸化細菌によるアンモニア酸化の工程において、１００〜３５０００（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）である場合、亜硝酸酸化細菌による亜硝酸酸化の工程において、１０００〜１２００００００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である場合、硝酸還元細菌による硝酸還元の工程において、５〜２５００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である場合、及び、亜硝酸還元細菌による亜硝酸還元の工程において、５〜３５００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）である場合には、通常、細菌１個あたりに単位時間に処理された処理対象物質の量に基づいたシミュレーション結果との間に大きな差が発生される可能性は低い。
【発明の効果】
【００２９】
本発明によれば、生物処理工程において単位数量の細菌に単位時間あたりに負荷される前記処理対象物質の量などが変動した場合などにおいても細菌による処理対象物質の最大反応速度をより正確に定義できシミュレーションの予測精度を向上させ得る。
したがって、シミュレーション方法ならびにシミュレーション装置における予測精度の低下を抑制しつつキャリブレーションの手間を削減させ得る。
また、実際の処理水質の確認頻度を低減させることができ、生物処理方法や生物処理装置の運転に要する手間を削減させ得る。
【発明を実施するための最良の形態】
【００３０】
本実施形態の生物学的水処理のシミュレーション方法について、窒素成分を処理対象物質として含有する排水を被処理水として、硝化・脱窒などの生物処理工程を実施する場合を例に説明する。
【００３１】
図１は、本実施形態のシミュレーション方法により各生物処理工程後の処理水の水質を予測しつつ排水の生物学的処理が実施される生物処理装置を示す概略ブロック図である。
図中の符号１は、排水を一連の処理工程に導入するための第一連絡管であり、符号２は、該第一連絡管１により排水（被処理水）が導入される硝化槽を表している。
【００３２】
図中の符号３は、硝化槽２から排出される処理液が流通される第二連絡管であり、符号４は、この第二連絡管３を通じて硝化槽２の処理液が被処理水として流入される脱窒槽を表している。
【００３３】
図中の符号５は、脱窒槽４から排出される処理液が流通される第三連絡管であり、符号６は、この第三連絡管５を通じて脱窒槽４の処理液が被処理水として流入される再曝気槽を表している。
【００３４】
図中の符号７は、再曝気槽６から排出される処理液が流通される第四連絡管であり、符号８は、この第四連絡管７を通じて再曝気槽６の処理液が被処理水として流入される沈殿槽を表している。
【００３５】
さらに、図中の符号９は、沈殿槽８から沈殿分離された上澄液が分離液として系外に排出される第五連絡管を表している。
また、ここでは、詳述しないがこの図１に示す生物処理装置には、前記沈殿槽８で沈殿分離された汚泥の一部を系外に余剰汚泥として排出するための連絡管（以下「汚泥引抜き配管」ともいう）ならびに、前記沈殿槽８で沈殿分離された汚泥の一部を硝化槽２に返送するための連絡管（以下「返送汚泥配管」ともいう）とが備えられている。
【００３６】
前記硝化槽２に第一連絡管１を通じて導入される被処理水（排水）には、通常、アンモニア性窒素が含有されている。
そして、前記硝化槽２には、アンモニア酸化細菌、亜硝酸酸化細菌などを含む活性汚泥が収容されており、前記被処理水の導入によって硝化槽２に導入される処理対象物質であるアンモニア性窒素は、前記活性汚泥と前記被処理水とにより硝化槽２に形成されてなる混合相中において細菌によって分解（酸化）される。
【００３７】
前記アンモニア酸化細菌としては、例えば、ニトロソモナス（Ｎｉｔｒｏｓｏｍｏｎａｓ）、ニトロソコッカス（Ｎｉｔｒｏｓｏｃｏｃｃｕｓ）などが挙げられ、前記亜硝酸酸化細菌としては、例えば、ニトロバクター（Ｎｉｔｒｏｂａｃｔｅｒ）、ニトロスピラ（Ｎｉｔｒｏｓｐｉｒａ）などが挙げられる。
【００３８】
また、酸化能力を発揮する細菌として、被処理水中の有機物質の分解において酸化能力を発揮する細菌が活性汚泥中にさらに含有されていてもよく、このような細菌としては、例えば、バチルス（Ｂａｃｉｌｌｕｓ）属細菌、ズーグレア（Ｚｏｏｇｌｏｅａ）属細菌、マイクロコッカス（Ｍｉｃｒｏｃｏｃｃｕｓ）属細菌などが挙げられる。
【００３９】
この硝化槽２には、流入される被処理水と活性汚泥とによってこの硝化槽２に形成される混合相について各種測定値が実施される。
この硝化槽２においては、この測定値をパラメータとして脱窒槽４に向けて流下される処理水の水質を予測するシミュレーション装置が設けられている。
このシミュレーション装置には、汚泥に含まれている細菌による硝化反応の最大反応速度の値がパラメータに用いられており、しかも、この最大反応速度の値は、この硝化槽２においてアンモニア酸化細菌に負荷されるアンモニア性窒素の量などと関数関係を有する状態でパラメータに用いられている。
このシミュレーション装置によるシミュレーション方法については、後段において詳述する。
【００４０】
前記脱窒槽４には、硝酸還元細菌、亜硝酸還元細菌、亜酸化窒素還元細菌等を含む活性汚泥が収容されている。
そして、この脱窒槽４には、前記硝化槽２の処理水（脱窒槽４における被処理水）の導入によって処理対象物質である硝酸性窒素、亜硝酸性窒素、および、亜酸化窒素などが導入される。
この導入された処理対象物質は、前記活性汚泥と前記被処理水とにより脱窒槽４に形成される混合相中において前記細菌によって分解（還元）される。
そして、この細菌によって亜硝酸性窒素や硝酸性窒素などを窒素ガスに還元し、該窒素ガスを大気中に放散することで被処理水から窒素成分が除去される。
【００４１】
前記脱窒槽４の活性汚泥には、処理対象物質に対して還元能力を発揮する細菌として、例えば、脱窒活性に関与する細菌を含有させることができ、例えば、アルカリジェネス（Ａｌｃａｌｉｇｅｎｅｓ）属細菌、アゾアルカス（Ａｚｏａｒｃｕｓ）属細菌、パラッコッカス（Ｐａｒａｃｏｃｃｕｓ）属細菌、シュードモナス（Ｐｓｅｕｄｏｍｏｎａｓ）属細菌などが挙げられる。
【００４２】
この脱窒槽４においては、硝化槽２から流入される被処理水と活性汚泥とによって形成される混合相について各種測定が実施される。
そして、この測定値をパラメータとして、再曝気槽６に向けて流下される処理水の水質を予測するシミュレーション装置が設けられている。
このシミュレーション装置には、汚泥に含まれている細菌による脱窒反応の最大反応速度の値がパラメータに用いられており、しかも、この最大反応速度の値は、この脱窒槽４において細菌に負荷される処理対象物質の量と関数関係を有する状態で用いられている。
【００４３】
このシミュレーション装置によるシミュレーション方法については、後段において詳述する。
【００４４】
前記再曝気槽６には、脱窒槽４の処理水（脱窒槽４から導入される被処理水）中に残留する有機物を好気条件下において分解除去するための曝気手段（図示せず）が備えられている。
【００４５】
前記沈殿槽８は、前記再曝気槽６から導入される被処理水に含有される活性汚泥などの固形成分を液体成分と沈殿分離するために必要な平均滞留時間を確保すべく、被処理水の流入量に対する十分な容積を有するよう形成されている。
【００４６】
次いで、シミュレーション装置ならびにシミュレーション方法について説明する。
【００４７】
このシミュレーション装置に用いられる水質のシミュレーションモデルとしては、最大反応速度の値と、生物処理工程において細菌１個あたりに単位時間（例えば、一日）に負荷される処理対象物質の量との関数関係を、ＩＷＡ（世界水協会）により作成された、ＡＳＭ１、ＡＳＭ２、ＡＳＭ２ｄ、ＡＳＭ３などに組み込んで改良したものを用いることができる。
なかでも、ＡＳＭ３は、モデルの拡張が容易であり、ＡＳＭ３に基づいたシミュレーションモデルを用いることが好ましい。
【００４８】
例えば、このＡＳＭ３に基づいたモデルに対して、最大反応速度の値と生物処理工程において細菌１個あたりに一日に負荷される前記処理対象物質の量との間に所定の関係が保たれるように関数をシミュレーション装置に組み込んでシミュレーションを実施する。
【００４９】
ここで、この最大反応速度の値と生物処理工程において細菌１個あたりに一日に負荷される前記処理対象物質の量との間に設ける関数の定義について説明する。
この関数については、最大反応速度の値と生物処理工程において細菌１個あたりに一日に負荷される前記処理対象物質の量との相関関係を求めるべく複数のデータを採取し、該データ群による回帰分析を実施するなどして定めることができる。
【００５０】
この最大反応速度：Ｖ（ｆｇ／ｃｏｐｙ／ｈ）については、生物処理工程において単位時間あたりの処理対象物質の濃度変化の量：ΔＳ（ｆｇ／ｍ3／ｈ）と、生物処理工程において用いられている単位体積あたりの細菌数：ｎ（ｃｏｐｉｅｓ／ｍ3）とを求めて下記式（４）
【数１５】
0004278701.tif 000004
を計算して求めることができる。
【００５１】
なお、細菌数（ｃｏｐｙ）の単位については、換算係数を用いて、生物量（ｇ−ＣＯＤCr）の単位に置き換えた状態でシミュレーションに採用することも可能であり、例えば、細菌数：１ｃｏｐｙ＝３．９６５×１０-10ｍｇ−ＣＯＤCr＝３．９６５×１０-13ｇ−ＣＯＤCrなる換算係数で細菌数（ｃｏｐｙ）を生物量（ｇ−ＣＯＤCr）の単位に換算して採用することができる。
【００５２】
この生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量：Ｌ（ｆｇ／ｃｏｐｙ／ｄａｙ）については、生物処理工程において単位時間あたりに単位体積あたりに負荷される処理対象物質の量：Ｑ（ｆｇ／ｍ3／ｄａｙ）と、生物処理工程において用いられている単位体積あたりの細菌数：ｎ（ｃｏｐｉｅｓ／ｍ3）とを求めて下記式（５）
【数１６】
0004278701.tif 000005
を計算して求めることができる。
【００５３】
なお、生物処理工程において単位時間あたりに単位体積あたりに負荷される処理対象物質の量（Ｑの値）については、１日間の実測を必要とするものではなく、例えば、半日の負荷の値を実測してその値を倍にして１日あたりの値に換算することもでき、逆に、２日以上の実測を実施して得られた実測値を測定日数で除して１日あたりの値に換算することもできる。
【００５４】
そして、生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量：Ｌ（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）の値が変化した場合に最大反応速度の値：Ｖ（ｆｇ・ｃｏｐｙ-1・ｈ-1）の値がどのように変化するかを、例えば、数点〜数十点のデータを採取して、該データをもとに分析を実施する。
なお、データ採取の点数は、少なくとも４点以上とするのが好ましい。
【００５５】
通常、これらのデータは、生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量をｘ軸（横軸）とし、最大反応速度の値をｙ軸（縦軸）とする平面にプロットすると右肩上がりの点群を形成する。
なお、最大反応速度の値については、後段において説明する測定方法によって定め得ることができる。
【００５６】
すなわち、これらの点群は、下記式（２）
【数１７】
0004278701.tif 000006
におけるｙ1の値をｘの値の増大とともに増大させる関数ｇ（ｘ）と
下記式（３）
【数１８】
0004278701.tif 000007
におけるｙ2の値をｘの値の増大とともに増大させ、しかも、ｙ2＞ｙ1となる関数ｈ（ｘ）との間に、その殆どを存在させることができる。
そして、このｇ（ｘ）とｈ（ｘ）との間を通る関数ｆ（ｘ）を適宜設定して、生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量（Ｌ値）との間に下記式（１）
【数１９】
0004278701.tif 000008
となる関数関係を有する状態で、この最大反応速度の値：Ｖ（ｆｇ・ｃｏｐｙ-1・ｈ-1）をパラメータとしてシミュレーション装置に用いることにより処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
したがって、この関数ｆ（ｘ）については、通常、生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量の値（Ｌ）の増大にともなって、最大反応速度の値（Ｖ）を増大させる関数となる。
そして、一般的な生物処理においてシミュレーション結果を実際の処理水質により近似させ得る点においてこの関数ｆ（ｘ）としては、下記一般式（６）
【数２０】
0004278701.tif 000009
（ただし、ａ、ｂ、ｃは定数）などで定義される関数から選定することが好ましい。
【００５７】
次いで、硝化、脱窒の各生物処理工程におけるシミュレーション装置ならびにシミュレーション方法について、より詳しく説明する。
【００５８】
（硝化工程）
前記硝化槽２においては、主として、アンモニア性窒素がアンモニア酸化細菌によって亜硝酸性窒素に分解されるアンモニア酸化と、該アンモニア酸化によって形成された亜硝酸性窒素の一部又は全部が亜硝酸酸化細菌によって硝酸性窒素に酸化される亜硝酸酸化との二つの工程が実施されることから、それぞれについてパラメータを設定する。
【００５９】
（アンモニア酸化について）
硝化槽２における、アンモニア性窒素の酸化反応は、例えば、従来のモデルでは、下記式（７）のようにして与えられてきた。
【００６０】
【数２１】
0004278701.tif 000010
ただし、
μNH4：最大比増殖速度（１／ｄａｙ）
ＹNH4：増殖収率（ｇ−ＣＯＤCr／ｇ ＮＨ4−Ｎ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNH4：溶解性のアンモニア濃度（ｇ ＮＨ4−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNH4：アンモニア飽和係数（ｇ ＮＨ4−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸNH4：硝化槽内のアンモニア酸化細菌濃度（ｇ−ＣＯＤCr／ｍ3）
【００６１】
そして、このμNH4、ＹNH4については定数であることから、単位アンモニア酸化細菌量が一日あたりに分解する処理対象物質の最大質量（最大反応速度）である（μNH4／ＹNH4）の値は、従来のモデルにおいては定数として与えられてきた。
【００６２】
一方で本実施形態にかかるシミュレーション装置には、上記式（７）のモデルに代えて、下記式（８）を組み込んで用いる。
【００６３】
【数２２】
0004278701.tif 000011
ただし、
ＬAOB：アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量（ｆｇ  ＮＨ4−Ｎ／ｃｏｐｙ／ｄａｙ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNH4：溶解性のアンモニア濃度（ｇ ＮＨ4−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNH4：アンモニア飽和係数（ｇ ＮＨ4−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸAOB：硝化槽内のアンモニア酸化細菌濃度（ｃｏｐｉｅｓ／ｍ3）
なお、細菌数（ｃｏｐｙ）に代えて、生物量（ｇ−ＣＯＤCr）をシミュレーションに採用することが可能である点については先に述べたとおりである。
【００６４】
そして、単位アンモニア酸化細菌量が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度（式（７）における（μNH4／ＹNH4））をアンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量を変数とする関数（ｆAOB（ＬAOB））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【００６５】
なお、この式（８）に用いる関数を、変数をｘ（ただし、ｘ＞０）としたときに下記式（９）
【数２３】
0004278701.tif 000012
におけるｙAOB1の値をｘの値の増大とともに増大させる関数ｇAOB（ｘ）と下記式（１０）
【数２４】
0004278701.tif 000013
におけるｙAOB2の値をｘの値の増大とともに増大させ、しかも、ｙAOB2＞ｙAOB1となる関数ｈAOB（ｘ）との間を通る関数とすることによってシミュレーションの予測精度をいっそう向上させ得る。
通常、アンモニアや硝酸等の窒素化合物以外に有機化合物等があまり含有されていない被処理水（処理対象物が主として窒素成分である排水）の生物学的処理に用いられているアンモニア酸化細菌によるアンモニア酸化の反応において、このアンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量：ＬAOB（ｆｇ  ＮＨ4−Ｎ／ｃｏｐｙ／ｄａｙ）が１００≦ＬAOB≦４０００の範囲では、前記式（９）における関数ｇAOB（ｘ）を下記式（１１）
【数２５】
0004278701.tif 000014
（ただし、ａ1＝７．０×１０2、ｂ1＝５．５×１０3、ｃ1＝−７．０×１０）
とし、前記式（１０）における関数ｈAOB（ｘ）を下記式（１２）
【数２６】
0004278701.tif 000015
（ただし、ａ2＝７．０×１０2、ｂ2＝５．５×１０3、ｃ2＝８．０×１０）で与えられる関数とすることでシミュレーションの予測精度をいっそう向上させ得る。
【００６６】
すなわち、例えば、このアンモニア酸化細菌による生物処理工程における最大反応速度を、以下の式（１３）
【数２７】
0004278701.tif 000016
ただし、
ＶAOB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬAOB：アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量（ｆｇ  ＮＨ4−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
により定義付けてシミュレーション装置のパラメータとして用いることができる。
【００６７】
また、下水なども含めて、一般に窒素成分が処理対象物質として含有されている被処理水の生物学的処理に用いられているアンモニア酸化細菌全般におけるアンモニア酸化反応の予測を行う場合においては、前記式（１１）におけるａ1、ｂ1、ｃ1の値を、それぞれ、ａ1＝４．０×１０3、ｂ1＝１．０×１０4、ｃ1＝−２．５×１０3とし、前記式（１２）におけるａ2、ｂ2、ｃ2の値を、それぞれ、ａ2＝４．０×１０3、ｂ2＝１．０×１０4、ｃ2＝２．５×１０3とすることで従来のシミュレーションに比べて精度の高い予測値を得ることができる。
【００６８】
したがって、ＶAOBとＬAOBとを関係付ける関数については、アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量：ＬAOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が１．０×１０2≦ＬAOB≦３．５×１０4において、最大反応速度：ＶAOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）−２．５×１０3｝≦ＶAOB≦｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）＋２．５×１０3｝となるように定められることが好ましい。
この好ましい態様によれば、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得る。
【００６９】
また、窒素成分とともに被処理水に塩化物イオンが含有されている場合には、この塩化物イオン濃度によっても最大反応速度が影響される。
しかも、アンモニア酸化細菌は、塩化物イオン濃度の増大とともにアンモニア酸化反応の最大反応速度を増大させる傾向にあることから被処理水中のアンモニア性窒素のみならず塩化物イオンの量との間にも所定の関数関係を有する状態で最大反応速度の値をパラメータに用いることでシミュレーションの精度をより向上させることができる。
【００７０】
したがって、上記式（８）のモデルに代えて、下記式（１４）を組み込んで用いることでシミュレーションの精度をより向上させることができる。
【００７１】
【数２８】
0004278701.tif 000017
ただし、
ＬAOB：アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量（ｆｇ  ＮＨ4−Ｎ／ｃｏｐｙ／ｄａｙ）
ＤCL：硝化工程に流入される被処理水中の塩化物イオン濃度（ｍｇ／ｌ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNH4：溶解性のアンモニア濃度（ｇ ＮＨ4−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNH4：アンモニア飽和係数（ｇ ＮＨ4−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸAOB：硝化槽内のアンモニア酸化細菌濃度（ｃｏｐｉｅｓ／ｍ3）
なお、細菌数（ｃｏｐｙ）に代えて、生物量（ｇ−ＣＯＤCr）をシミュレーションに採用することが可能である点については先に述べたとおりである。
【００７２】
そして、一個のアンモニア酸化細菌が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度：ＶAOB）をアンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量を変数とする関数（ｆAOB（ＬAOB））と塩化物イオン濃度を変数とする関数（ｋAOB（ＤCL））との積（ＶAOB＝ｆAOB（ＬAOB）・ｋAOB（ＤCL））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【００７３】
この塩化物イオン濃度を変数とする関数（ｋAOB（ＤCL））は、硝化工程に流入される被処理水中の塩化物イオン濃度：ＤCL（ｍｇ／ｌ）がＤCL≦４．０×１０3の範囲において、（１＋４．４×１０-5・ＤCL）≦ｋAOB（ＤCL）≦（１＋１．６４×１０-4・ＤCL）となるように定められることが好ましい。
【００７４】
（亜硝酸酸化について）
硝化槽２における、亜硝酸性窒素の酸化反応は、例えば、従来のモデルでは、下記式（１５）のようにして与えられてきた。
【数２９】
0004278701.tif 000018
ただし、
μNO2：最大比増殖速度（１／ｄａｙ）
ＹNO2：増殖収率（ｇ−ＣＯＤCr／ｇ ＮＯ2−Ｎ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNO2：溶解性の亜硝酸濃度（ｇ ＮＯ2−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNO2：亜硝酸飽和係数（ｇ ＮＯ2−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸNO2：硝化槽内の亜硝酸酸化細菌濃度（ｇ−ＣＯＤCr／ｍ3）
【００７５】
そして、このμNO2、ＹNO2については定数であることから、単位亜硝酸酸化細菌量が一日あたりに分解する処理対象物質の最大質量（最大反応速度）である（μNO2／ＹNO2）の値は、従来のモデルにおいては定数として与えられてきた。
【００７６】
一方で本実施形態にかかるシミュレーション装置には、上記式（１５）のモデルに代えて、下記式（１６）を組み込んで用いる。
【００７７】
【数３０】
0004278701.tif 000019
ただし、
ＬNOB：亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｇ ＮＯ2−Ｎ／ｃｏｐｙ／ｄａｙ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNO2：溶解性の亜硝酸濃度（ｇ ＮＯ2−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNO2：亜硝酸飽和係数（ｇ ＮＯ2−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸNOB：硝化槽内の亜硝酸酸化細菌濃度（ｃｏｐｉｅｓ／ｍ3）
なお、細菌数（ｃｏｐｙ）に代えて、生物量（ｇ−ＣＯＤCr）をシミュレーションに採用することが可能である点については先に述べたとおりである。
【００７８】
そして、単位亜硝酸酸化細菌量が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度（式（１５）における（μNO2／ＹNO2））を、亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量を変数とする関数（ｆNOB（ＬNOB））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【００７９】
なお、この式（１６）に用いる関数を、変数をｘ（ただし、ｘ＞０）としたときに下記式（１７）
【数３１】
0004278701.tif 000020
におけるｙNOB1の値をｘの値の増大とともに増大させる関数ｇNOB（ｘ）と下記式（１８）
【数３２】
0004278701.tif 000021
におけるｙNOB2の値をｘの値の増大とともに増大させ、しかも、ｙNOB2＞ｙNOB1となる関数ｈNOB（ｘ）との間を通る関数とすることによってシミュレーションの予測精度をいっそう向上させ得る。
通常、アンモニアや硝酸等の窒素化合物以外に有機化合物等があまり含有されていない被処理水（処理対象物が主として窒素成分である排水）がアンモニア酸化細菌により生物学的処理された後の処理水の亜硝酸酸化に用いられている亜硝酸酸化細菌によって、亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNOB（ｆｇ  ＮＯ2−Ｎ／ｃｏｐｙ／ｄａｙ）が１０００≦ＬNOB≦６００００となる範囲で亜硝酸性窒素の酸化を実施させる場合には、前記式（１７）における関数ｇNOB（ｘ）を下記式（１９）
【数３３】
0004278701.tif 000022
（ただし、ａ3＝２．５×１０4、ｂ3＝９．５×１０4、ｃ3＝−１．７×１０3）
とし、前記式（１８）における関数ｈNOB（ｘ）を下記式（２０）
【数３４】
0004278701.tif 000023
（ただし、ａ4＝２．５×１０4、ｂ4＝９．５×１０4、ｃ4＝２．２×１０3）
として、このｇNOB（ｘ）とｈNOB（ｘ）との間を通って、そして、亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量を変数とする関数によってこの亜硝酸酸化細菌による亜硝酸性窒素の酸化反応の最大反応速度を定義することにより、シミュレーションの予測精度をいっそう向上させ得る。
【００８０】
すなわち、例えば、この亜硝酸酸化細菌による生物処理工程における最大反応速度を、以下の式（２１）
【数３５】
0004278701.tif 000024
ただし、
ＶNOB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNOB：亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｆｇ ＮＯ2−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
により定義付けてシミュレーションのパラメータとして用いることができる。
【００８１】
また、下水なども含めて、一般に窒素成分が処理対象物質として含有されている被処理水の生物学的処理に用いられている亜硝酸酸化細菌全般における亜硝酸酸化反応の予測を行う場合においては、前記式（１９）におけるａ3、ｂ3、ｃ3の値を、それぞれ、ａ3＝２．５×１０5、ｂ3＝４．０×１０5、ｃ3＝−１．０×１０5とし、前記式（２０）におけるａ4、ｂ4、ｃ4の値を、それぞれ、ａ4＝２．５×１０5、ｂ4＝４．０×１０5、ｃ4＝１．０×１０5とすることで従来のシミュレーションに比べて精度の高い予測値を得ることができる。
【００８２】
したがって、ＶNOBとＬNOBとを関係付ける関数については、亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が１．０×１０3≦ＬNOB≦１．２×１０6となる範囲において、最大反応速度：ＶNOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）−１．０×１０5｝≦ＶNOB≦｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）＋１．０×１０5｝となるように定められることが好ましい。
この好ましい態様によれば、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得る。
【００８３】
また、窒素成分とともに被処理水に塩化物イオンが含有されている場合には、この塩化物イオン濃度によっても最大反応速度が影響される。
しかも、亜硝酸酸化細菌は、塩化物イオン濃度の増大とともに亜硝酸酸化反応の最大反応速度を低下させる傾向にあることから被処理水中の亜硝酸性窒素のみならず塩化物イオンの量との間にも所定の関数関係を有する状態で最大反応速度の値をパラメータに用いることでシミュレーションの精度をより向上させることができる。
【００８４】
したがって、上記式（１６）のモデルに代えて、下記式（２２）を組み込んで用いることでシミュレーションの精度をより向上させることができる。
【００８５】
【数３６】
0004278701.tif 000025
ただし、
ＬNOB：亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｇ ＮＯ2−Ｎ／ｃｏｐｙ／ｄａｙ）
ＤCL：硝化工程に流入される被処理水中の塩化物イオン濃度（ｍｇ／ｌ）
ＳO2：溶存酸素濃度（ｇ Ｏ2／ｍ3）
ＳNO2：溶解性の亜硝酸濃度（ｇ ＮＯ2−Ｎ／ｍ3）
ＳALK：アルカリ度（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＫO2：溶存酸素飽和係数（ｇ Ｏ2／ｍ3）
ＫNO2：亜硝酸飽和係数（ｇ ＮＯ2−Ｎ／ｍ3）
ＫALK：アルカリ度飽和係数（ｍｏｌｅ  ＨＣＯ3／ｍ3）
ＸNOB：硝化槽内の亜硝酸酸化細菌濃度（ｃｏｐｉｅｓ／ｍ3）
なお、細菌数（ｃｏｐｙ）に代えて、生物量（ｇ−ＣＯＤCr）をシミュレーションに採用することが可能である点については先に述べたとおりである。
【００８６】
そして、一個の亜硝酸酸化細菌が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度：ＶNOB）を亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量を変数とする関数（ｆNOB（ＬNOB））と塩化物イオン濃度を変数とする関数（ｋNOB（ＤCL））との積（ＶNOB＝ｆNOB（ＬNOB）・ｋNOB（ＤCL））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【００８７】
この塩化物イオン濃度を変数とする関数（ｋNOB（ＤCL））は、硝化工程に流入される被処理水中の塩化物イオン濃度：ＤCL（ｍｇ／ｌ）がＤCL≦４．０×１０3の範囲において、（１−８．７×１０-5・ＤCL）≦ｋNOB（ＤCL）≦（１−４．０×１０-5・ＤCL）となるように定められることが好ましい。
【００８８】
（脱窒工程について）
前記脱窒槽４においては、主として、硝酸性窒素が硝酸還元細菌によって亜硝酸性窒素に還元される亜硝酸還元と、前記硝酸還元細菌によって形成された亜硝酸性窒素が亜硝酸還元細菌によって窒素に還元される亜硝酸還元との二つの工程が実施されることから、それぞれについてパラメータを設定する。
【００８９】
（硝酸還元について）
前記脱窒槽４における、硝酸還元細菌による硝酸性窒素の還元反応については、変数をｘ（ただし、ｘ＞０）としたときに下記式（２３）
【数３７】
0004278701.tif 000026
におけるｙNARB1の値をｘの値の増大とともに増大させる関数ｇNARB（ｘ）と下記式（３）
【数３８】
0004278701.tif 000027
におけるｙNARB2の値をｘの値の増大とともに増大させ、しかも、ｙNARB2＞ｙNARB1となる関数ｈNARB（ｘ）との間を通る関数で硝酸還元細菌による硝酸性窒素の還元反応の最大反応速度を定義することにより、シミュレーションの予測精度をいっそう向上させ得る。
通常、アンモニアや硝酸等の窒素化合物以外に有機化合物等があまり含有されていない被処理水（処理対象物が主として窒素成分である排水）がアンモニア酸化細菌及び亜硝酸酸化細菌により生物学的処理（硝化工程）された後の処理水の硝酸還元に用いられている硝酸還元細菌によって、硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素量：ＬNARB（ｆｇ  ＮＯ3−Ｎ／ｃｏｐｙ／ｄａｙ）が、５．０≦ＬNARB≦７０の範囲で硝酸性窒素の還元を実施させる場合には、前記式（２３）における関数ｇNARB（ｘ）を下記式（２５）
【数３９】
0004278701.tif 000028
（ただし、ａ5＝１．０×１０2、ｂ5＝８．５×１０2  、ｃ5＝−２．０）
とし、前記式（２４）における関数ｈNARB（ｘ）を下記式（２６）
【数４０】
0004278701.tif 000029
（ただし、ａ6＝１．０×１０2、ｂ6＝８．５×１０2、ｃ6＝１．８）
として、このｇNARB（ｘ）とｈNARB（ｘ）との間を通って、そして、硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素を変数とする関数によってこの硝酸還元細菌による硝酸性窒素の還元反応の最大反応速度を定義することにより、シミュレーションの予測精度をいっそう向上させ得る。
【００９０】
すなわち、例えば、この硝酸還元細菌による生物処理工程における最大反応速度を、以下の式（２７）
【数４１】
0004278701.tif 000030
ただし、
ＶNARB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNARB：硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素量（ｆｇ ＮＯ3−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
により定義付けてシミュレーションのパラメータとして用いることができる。
【００９１】
また、下水なども含めて、一般に窒素成分が処理対象物質として含有されている被処理水の生物学的処理に用いられている硝酸還元細菌全般における硝酸還元反応の予測を行う場合においては、前記式（２５）におけるａ5、ｂ5、ｃ5の値を、それぞれ、ａ5＝２．２×１０2、ｂ5＝７．０×１０2、ｃ5＝−１．７×１０2とし、前記式（２６）におけるａ6、ｂ6、ｃ6の値を、それぞれ、ａ6＝２．２×１０2、ｂ6＝７．０×１０2、ｃ6＝７０とすることで従来のシミュレーションに比べて精度の高い予測値を得ることができる。
【００９２】
したがって、ＶNARBとＬNARBとを関係付ける関数については、硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNARB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が５．０≦ＬNARB≦２５００となる範囲において、最大反応速度：ＶNARB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）−１．７×１０2｝≦ＶNARB≦｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）＋７０｝となるように定められることが好ましい。
この好ましい態様によれば、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得る。
【００９３】
また、窒素成分とともに被処理水に塩化物イオンが含有されている場合には、この塩化物イオン濃度によっても最大反応速度が影響される。
しかも、硝酸還元細菌は、塩化物イオン濃度の増大とともに硝酸還元反応の最大反応速度を低下させる傾向にあることから被処理水中の硝酸性窒素のみならず塩化物イオンの量との間にも所定の関数関係を有する状態で最大反応速度の値をパラメータに用いることでシミュレーションの精度をより向上させることができる。
【００９４】
したがって、一個の硝酸還元細菌が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度：ＶNARB）を、硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量を変数とする関数（ｆNARB（ＬNARB））と塩化物イオン濃度を変数とする関数（ｋNARB（ＤCL））との積（ＶNARB＝ｆNARB（ＬNARB）・ｋNARB（ＤCL））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【００９５】
この塩化物イオン濃度を変数とする関数（ｋNARB（ＤCL））は、硝化工程に流入される被処理水中の塩化物イオン濃度：ＤCL（ｍｇ／ｌ）がＤCL≦４．０×１０3の範囲において、（１−７．９×１０-5・ＤCL）≦ｋNARB（ＤCL）≦（１−１．０×１０-5・ＤCL）となるように定められることが好ましい。
【００９６】
（亜硝酸還元について）
前記脱窒槽４における、亜硝酸還元細菌による亜硝酸性窒素の還元反応については、変数をｘ（ただし、ｘ＞０）としたときに下記式（２８）
【数４２】
0004278701.tif 000031
におけるｙNIRB1の値をｘの値の増大とともに増大させる関数ｇNIRB（ｘ）と下記式（２９）
【数４３】
0004278701.tif 000032
におけるｙNIRB2の値をｘの値の増大とともに増大させ、しかも、ｙNIRB2＞ｙNIRB1となる関数ｈNIRB（ｘ）との間を通る関数で亜硝酸還元細菌による亜硝酸性窒素などの還元反応の最大反応速度を定義することにより、シミュレーションの予測精度をいっそう向上させ得る。
通常、アンモニアや硝酸等の窒素化合物以外に有機化合物等があまり含有されていない被処理水（処理対象物が主として窒素成分である排水）がアンモニア酸化細菌及び亜硝酸酸化細菌により生物学的処理（硝化工程）された後の処理水の還元に用いられている亜硝酸還元細菌によって、亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNIRB（ｆｇ  ＮＯ2−Ｎ／ｃｏｐｙ／ｄａｙ）が、５．０≦ＬNIRB≦１２０の範囲で亜硝酸性窒素の還元を実施させる場合には、前記式（２８）における関数ｇNIRB（ｘ）を下記式（３０）
【数４４】
0004278701.tif 000033
（ただし、ａ7＝６．０×１０1、ｂ7＝３．５×１０2  、ｃ7＝−７．０）
とし、前記式（２９）における関数ｈNIRB（ｘ）を下記式（３１）
【数４５】
0004278701.tif 000034
（ただし、ａ8＝６．０×１０1、ｂ8＝３．５×１０2、ｃ8＝３．０）
として、このｇNIRB（ｘ）とｈNIRB（ｘ）との間を通って、そして、亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素を変数とする関数によってこの亜硝酸還元細菌による亜硝酸性窒素などの還元反応の最大反応速度を定義することにより、シミュレーションの予測精度をいっそう向上させ得る。
【００９７】
すなわち、例えば、この亜硝酸還元細菌による生物処理工程における最大反応速度を、以下の式（３２）
【数４６】
0004278701.tif 000035
ただし、
ＶNIRB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNIRB：亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｆｇ ＮＯ2−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
により定義付けてシミュレーションのパラメータとして用いることができる。
【００９８】
また、下水なども含めて、一般に窒素成分が処理対象物質として含有されている被処理水の生物学的処理に用いられている硝酸還元細菌全般における硝酸還元反応の予測を行う場合においては、前記式（３０）におけるａ7、ｂ7、ｃ7の値を、それぞれ、ａ7＝７．０×１０2、ｂ7＝２．５×１０3、ｃ7＝−２．５×１０2とし、前記式（３１）におけるａ8、ｂ8、ｃ8の値を、それぞれ、ａ8＝７．０×１０2、ｂ8＝２．５×１０3、ｃ8＝２．５×１０2とすることで従来のシミュレーションに比べて精度の高い予測値を得ることができる。
【００９９】
したがって、ＶNIRBとＬNIRBとを関係付ける関数については、亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNIRB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が５．０≦ＬNIRB≦３．５×１０3となる範囲において、最大反応速度：ＶNIRB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）−２．５×１０2｝≦ＶNIRB≦｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）＋２．５×１０2｝となるように定められることが好ましい。
この好ましい態様によれば、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得る。
【０１００】
また、窒素成分とともに被処理水に塩化物イオンが含有されている場合には、この塩化物イオン濃度によっても最大反応速度が影響される。
しかも、亜硝酸還元細菌は、塩化物イオン濃度の増大とともに亜硝酸還元反応の最大反応速度を低下させる傾向にあることから被処理水中の亜硝酸性窒素のみならず塩化物イオンの量との間にも所定の関数関係を有する状態で最大反応速度の値をパラメータに用いることでシミュレーションの精度をより向上させることができる。
【０１０１】
したがって、一個の亜硝酸還元細菌が一日あたりに分解する処理対象物質の最大質量（すなわち、最大反応速度：ＶNIRB）を、亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量を変数とする関数（ｆNIRB（ＬNIRB））と塩化物イオン濃度を変数とする関数（ｋNIRB（ＤCL））との積（ＶNIRB＝ｆNIRB（ＬNIRB）・ｋNIRB（ＤCL））とすることにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させることができ、キャリブレーションの手間を削減させ得る。
【０１０２】
この塩化物イオン濃度を変数とする関数（ｋNIRB（ＤCL））は、硝化工程に流入される被処理水中の塩化物イオン濃度：ＤCL（ｍｇ／ｌ）がＤCL≦４．０×１０3の範囲において、（１−６．０×１０-5・ＤCL）≦ｋNIRB（ＤCL）≦（１−５．０×１０-6・ＤCL）となるように定められることが好ましい。
【０１０３】
なお、上記に説明したように、主としてアンモニアや硝酸等の窒素化合物が処理対象物質であり、他に有機性化合物を多く含有していない被処理水の生物処理の工程における細菌の反応については、例えば、硝化槽内において処理対象物質がアンモニア性窒素で、細菌がアンモニア酸化細菌である場合においては、前記式（１１）に示す関数ｇAOB（ｘ）と、前記式（１２）に示す関数ｈAOB（ｘ）の間を通る関数、例えば、前記式（１３）に示す関数によって最大反応速度を定義することで通常、精度の高いシミュレーションを実施させることができる。
また、処理対象物質が亜硝酸性窒素で、細菌が亜硝酸酸化細菌である場合においては、前記式（１９）に示す関数ｇNOB（ｘ）と、前記式（２０）に示す関数ｈNOB（ｘ）の間を通る関数、例えば、前記式（２１）に示す関数によって最大反応速度を定義してシミュレーション装置に組み込むことで、通常、精度の高いシミュレーションを実施させることができる。
【０１０４】
また、脱窒槽内において処理対象物質が硝酸性窒素で、細菌が硝酸還元細菌である場合においては、前記式（２５）に示す関数ｇNARB（ｘ）と、前記式（２６）に示す関数ｈNARB（ｘ）の間を通る関数、例えば、前記式（２７）に示す関数によって最大反応速度を定義することで、通常、精度の高いシミュレーションを実施させることができる。
また、脱窒槽内において処理対象物質が亜硝酸性窒素で、細菌が亜硝酸還元細菌である場合においては、前記式（３０）に示す関数ｇNIRB（ｘ）と、前記式（３１）に示す関数ｈNIRB（ｘ）の間を通る関数、例えば、前記式（３２）に示す関数によって最大反応速度を定義することで、通常、精度の高いシミュレーションを実施させることができる。
【０１０５】
したがって、生物処理工程に硝化工程を有する場合は、前記式（１３）に示す関数および前記式（２１）に示す関数によってアンモニア酸化ならびに亜硝酸酸化の最大反応速度を定義してシミュレーション装置に組み込み、生物処理工程にさらに脱窒工程を有する場合は、前記式（１３）に示す関数、前記式（２１）に示す関数に加え前記式（２７）に示す関数と前記式（３２）に示す関数によってアンモニア酸化、亜硝酸酸化、硝酸還元及び亜硝酸還元の最大反応速度を定義してシミュレーション装置に組み込むことで、より精度の高いシミュレーションを実施させることができる。
なお、要すれば、下記に示す方法により、最大反応速度と細菌１個あたりに一日に負荷される処理対象物質の量との関係を求めてシミュレーション装置に組み込むようにしてもよい。
また、下記方法により細菌１個あたりに一日に負荷される処理対象物質の量を求めるとともに処理後に残留する処理対象物質の量を求めて、細菌１個あたりに一日に処理された処理対象物質の量を計算し、最大反応速度と細菌１個あたりに一日に処理された処理対象物質の量との関係を示す関数をシミュレーション装置に組み込むようにしてもよい。
【０１０６】
（最大反応速度）
（アンモニア酸化速度）
硝化槽の活性汚泥（硝化汚泥）中に含まれる細菌によるアンモニア酸化速度の測定（アンモニア酸化速度試験）は、以下のように行なう。
【０１０７】
５００ｍｌ容三角フラスコに、希釈水［１ｌあたりの組成：炭酸水素ナトリウム  ２４０ｍｇ、ＢＯＤ−Ａ液〔ＪＩＳ  Ｋ  ０１０２の２１の項に従う、緩衝液（ｐＨ７．２）〕１ｍｌ、ＢＯＤ−Ｂ液（ＪＩＳ  Ｋ  ０１０２の２１の項に従う、硫酸マグネシウム溶液）  １ｍｌ、ＢＯＤ−Ｃ液（ＪＩＳ  Ｋ  ０１０２の２１の項に従う、塩化カルシウム溶液）  １ｍｌ、ＢＯＤ−Ｄ液〔ＪＩＳ  Ｋ  ０１０２の２１の項に従う、塩化鉄（ＩＩＩ）溶液〕１ｍｌ、残部  水］  ３９０ｍｌを入れ、１０００ｍｇ−Ｎ／ｌの塩化アンモニウム水溶液  １０ｍｌを添加して、混合物Ａを調製する。ついで、５００ｍｌ容三角フラスコ内の混合物Ａを、恒温槽中、３０℃に維持し、撹拌しながら、１０分間以上曝気し、溶液Ａを得る。ここで、溶液ＡのｐＨを測定する。
【０１０８】
硝化汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、硝化汚泥試料を得る。
【０１０９】
その後、得られた硝化汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記溶液Ａと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、塩化アンモニウムに由来するアンモニア性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１１０】
また、前記硝化汚泥試料と前記溶液Ａとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、ＪＩＳ  Ｋ  ０１０２の１４の項に従い、汚泥濃度を測定する。
【０１１１】
塩化アンモニウムに由来するアンモニア性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１１２】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４２．５の項に従う）を行なうことにより、アンモニア性窒素の量を測定する。分析結果より単位時間あたりのアンモニア性窒素の量の変化を算出し、これをアンモニア酸化の最大反応速度とすることができる。
【０１１３】
（亜硝酸酸化速度）
一方、硝化汚泥中に含まれる細菌による亜硝酸酸化速度の測定（亜硝酸酸化速度試験）は、以下のように行なう。
【０１１４】
５００ｍｌ容三角フラスコに、前記希釈水  ３９０ｍｌを入れ、１０００ｍｇ−Ｎ／ｌの亜硝酸ナトリウム水溶液  １０ｍｌを添加して、混合物Ｂを調製する。ついで、５００ｍｌ容三角フラスコ内の混合物Ｂを、恒温槽中、３０℃に維持し、撹拌しながら、１０分間以上曝気し、溶液Ｂを得る。ここで、溶液ＢのｐＨを測定する。
【０１１５】
硝化汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、硝化汚泥試料を得る。
【０１１６】
その後、得られた硝化汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記溶液Ｂと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、亜硝酸ナトリウムに由来する亜硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１１７】
また、前記硝化汚泥試料と前記溶液Ｂとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１１８】
亜硝酸ナトリウムに由来する亜硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１１９】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．１．２の項に従う）を行なうことにより、亜硝酸性窒素の量を測定する。分析結果より単位時間あたりの亜硝酸性窒素の量の変化を算出し、これを亜硝酸酸化の最大反応速度として求めることができる。
【０１２０】
（硝酸還元速度）
脱窒槽の活性汚泥（脱窒汚泥）中に含まれる細菌による硝酸還元速度の測定（硝酸還元速度試験）は、以下のように測定する。
５００ｍｌ容三角フラスコに、前記希釈水３８０ｍｌを入れ、恒温槽中、３０℃に維持し、撹拌しながら、１０００ｍｇ−Ｎ／ｌの硝酸ナトリウム水溶液  １０ｍｌを添加し、ついで、メタノール水溶液〔５０００ｍｇメタノール／ｌ〕  １０ｍｌを添加して、混合物Ｃを調製する。その後、前記混合物Ｃに、散気球により窒素ガスを吹き込み、１０分間脱気する。ここで、前記混合物ＣのｐＨを測定する。前記５００ｍｌ容三角フラスコに、シリコン栓をし、撹拌しながら窒素ガスを吹き込み、該５００ｍｌ容三角フラスコ内の気相中の空気を除く。その後、１ｌ容テドラーバックに窒素ガスを吹き込み、前記５００ｍｌ容三角フラスコのシリコン栓に設置する。
【０１２１】
脱窒汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、脱窒汚泥試料を得る。
【０１２２】
その後、得られた脱窒汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記混合物Ｃと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、硝酸ナトリウムに由来する硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１２３】
また、前記脱窒汚泥試料と前記混合物Ｃとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１２４】
硝酸ナトリウムに由来する硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１２５】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．２．５の項に従う）を行なうことにより、硝酸性窒素の量を測定する。
分析結果より単位時間あたりの硝酸性窒素の量の変化を算出し、これを硝酸還元反応の最大反応速度として求めることができる。
【０１２６】
（亜硝酸還元速度）
一方、脱窒槽の活性汚泥（脱窒汚泥）中に含まれる細菌による亜硝酸還元速度の測定（亜硝酸還元速度試験）は、以下のように測定する。
【０１２７】
５００ｍｌ容三角フラスコに、前記希釈水３８０ｍｌを入れ、恒温槽中、３０℃に維持し、撹拌しながら、１０００ｍｇ−Ｎ／ｌの亜硝酸ナトリウム水溶液  １０ｍｌを添加し、ついで、メタノール水溶液〔５０００ｍｇメタノール／ｌ〕  １０ｍｌを添加して、混合物Ｄを調製する。その後、前記混合物Ｄに、散気球により窒素ガスを吹き込み、１０分間脱気する。ここで、前記混合物ＤのｐＨを測定する。前記５００ｍｌ容三角フラスコに、シリコン栓をし、撹拌しながら窒素ガスを吹き込み、該５００ｍｌ容三角フラスコ内の気相中の空気を除く。その後、１ｌ容テドラーバックに窒素ガスを吹き込み、前記５００ｍｌ容三角フラスコのシリコン栓に設置する。
【０１２８】
脱窒汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、脱窒汚泥試料を得る。
【０１２９】
その後、得られた脱窒汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記混合物Ｄと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、亜硝酸ナトリウムに由来する亜硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１３０】
また、前記脱窒汚泥試料と前記混合物Ｄとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１３１】
亜硝酸ナトリウムに由来する亜硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、活性汚泥濃度を測定する。
【０１３２】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．１．２の項に従う）を行なうことにより、亜硝酸性窒素の量を測定する。分析結果より単位時間あたりの亜硝酸性窒素の量の変化を算出し、これを亜硝酸還元反応の最大反応速度として求めることができる。
【０１３３】
（細菌１個あたりに一日に負荷される処理対象物質の量）
生物処理工程において細菌１個あたりに一日に負荷される前記処理対象物質の量は、生物処理工程に単位時間あたりに単位体積あたりに外部から導入される処理対象物質の量と、生物処理工程において単位時間あたりに単位体積あたりに産生される処理対象物質の量との合計量と、生物処理工程において用いられている単位体積あたりの細菌数とを求めることで計算により求め得る。
この内、生物処理工程に単位時間あたりに単位体積あたりに外部から導入される処理対象物質の量については、例えば、被処理水の流入量と該被処理水中の処理対象物質濃度とを測定して、一日あたりに硝化槽や脱窒槽などへ流入される処理対象物質の量を計算により求めて、硝化槽や脱窒槽などの容積で除して求めることができる。また、生物処理工程において産生される処理対象物質の量は、通常、この処理対象物質の元になる物質が外部から生物処理工程に導入される量と同等とみなすことができ、アンモニア酸化細菌が収容されている硝化槽における亜硝酸性窒素の産生量は、硝化槽に外部から導入されるアンモニア性窒素の量と同等とみなすことができる。さらに、単位体積あたりの細菌数については、硝化槽や脱窒槽などの汚泥を分析して細菌数を求めて硝化槽や脱窒槽などの容積で除して求めることができる。
【０１３４】
この内、被処理水中の処理対象物質濃度については、例えば、処理対象物質がアンモニア性窒素である場合には、ＪＩＳ  Ｋ  ０１０２にしたがって、インドフェノール青吸光光度法、中和滴定法、イオン電極法、イオンクロマトグラフ法などにより測定することができる。
また、亜硝酸性窒素については、ＪＩＳ  Ｋ  ０１０２にしたがって、吸光光度法、イオンクロマトグラフ法などにより測定することができる。
また、硝酸性窒素については、ＪＩＳ  Ｋ  ０１０２にしたがって、還元蒸留−吸光光度法、銅・カドミウムカラム還元−吸光光度法、イオンクロマトグラフ法などにより測定することができる。
さらに、全窒素の測定方法としては、例えば、総和法、ケルダール窒素法、還元蒸留ケルダール法、紫外線吸光光度法などが挙げられる（「下水試験方法」、上巻、１９９７年版、財団法人 日本下水道協会参照）。
【０１３５】
また、活性汚泥（硝化汚泥または脱窒汚泥）中に含まれる各種細菌の細菌数は、例えば、活性汚泥（硝化汚泥または脱窒汚泥）に含まれる細菌に由来するＤＮＡを抽出して、リアルタイムＰＣＲに代表される各種の核酸検出方法により定量することができる。
【０１３６】
活性汚泥（硝化汚泥または脱窒汚泥）に含まれる細菌に由来するＤＮＡの抽出は、土壌からのＤＮＡの抽出に用いられる手法、例えば、活性汚泥中の細菌を物理的手段（ビーズなど）により破砕し、ＤＮＡを抽出することなどにより行なわれうる。ＤＮＡの単離には、特に限定されないが、例えば、商品名：ＦａｓｔＤＮＡ  ＳＰＩＮ  Ｋｉｔ  ｆｏｒ  Ｓｏｉｌ〔キュービオジェン（Ｑｂｉｏｇｅｎｅ）社製〕、商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）などが用いられうる。具体的には、例えば、商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）を用いた場合、以下のように活性汚泥からＤＮＡを単離することができる。
【０１３７】
各反応槽から活性汚泥（硝化汚泥または脱窒汚泥）を採取し、２ｍｌ容のマイクロ遠心チューブに入れる。また、活性汚泥の固形物濃度（ＭＬＳＳ濃度）が、１５００〜２０００ｍｇ／ｌである場合、２ｍｌの活性汚泥、２０００〜３０００ｍｇ／ｌである場合、１．５ｍｌの活性汚泥、３０００〜５０００ｍｇ／ｌである場合、１ｍｌの活性汚泥、５０００〜７０００ｍｇ／ｌである場合、０．７ｍｌの活性汚泥、７０００〜１００００ｍｇ／ｌである場合、０．５ｍｌの活性汚泥を採取する。
【０１３８】
その後、前記マイクロ遠心チューブに入れた活性汚泥を、２０６３０×ｇ（１５００ｒｐｍ）、２分間、４℃の遠心分離および２０６３０×ｇ（１５０００ｒｐｍ）、３０秒間、４℃の遠心分離に供する。ついで、得られた活性汚泥を、予め６５℃に加温した４５０μｌのＬｙｓｉｓ  Ｓｏｌｕｔｉｏｎ  ＢＢ〔商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）に添付〕に懸濁させる。
【０１３９】
その後、得られた懸濁物を、Ｂｅａｄｓ  Ｔｕｂｅ〔商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）に添付〕に移す。また、元のチューブを、予め６５℃に加温した４５０μｌのＬｙｓｉｓ  Ｓｏｌｕｔｉｏｎ  ＢＢで洗浄し、洗浄後に得られた懸濁物を前記Ｂｅａｄｓ  Ｔｕｂｅに移す。前記Ｂｅａｄｓ  Ｔｕｂｅ中の懸濁物に、５０μｌのＬｙｓｉｓ  Ｓｏｌｕｔｉｏｎ  ２０Ｓ〔商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）に添付〕を添加し、混和させる。
【０１４０】
その後、得られた混合物を、６５℃で１５分間維持し、ついで、ビーズ式破砕機〔商品名：Ｂｅａｄｓ  Ｂｅａｔｅｒ（株式会社ニッポンジーン製）〕に供して、３０００ｒｐｍで、９０秒間、Ｂｅａｄｓ  Ｂｅａｔｉｎｇを行なう。その後、得られた産物を、６５℃で４０分間、穏やかに混合しながら維持し、ついで、１２０００×ｇ、１分間、２０℃で遠心分離し、上清  約６６０μｌを２ｍｌのチューブに回収する。前記上清に、４４０μｌのＰｕｒｉｆｉｃａｔｉｏｎ  Ｓｏｌｕｔｉｏｎ〔商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）に添付〕を添加し、混和させる。その後、６００μｌのクロロホルムを添加し、穏やかに撹拌し、ついで、１２０００×ｇ、１５分間、２０℃の遠心分離にて、水層９００μｌを２ｍｌチューブに回収する。得られた産物に、等量（９００μｌ）のＰｒｅｃｉｐｉｔａｔｉｏｎ  Ｓｏｌｕｔｉｏｎ〔商品名：ＩＳＯＩＬ  ｆｏｒ  Ｂｅａｄｓ  Ｂｅａｔｉｎｇ（株式会社ニッポンジーン製）に添付〕を添加し、混合する。得られた産物を遠心分離し、得られた沈殿物を、Ｗａｓｈ  Ｓｏｌｕｔｉｏｎで洗浄する。その後、得られた沈殿物に１ｍｌの７０容量％  エタノールと２μｌの商品名：Ｅｔｈａｃｈｉｎｍａｔｅ（株式会社ニッポンジーン製）とを添加して、エタノール沈殿を行ない、ＤＮＡの沈殿物を得る。得られたＤＮＡの沈殿物に、２００μｌのＴＥ緩衝液（組成：１０ｍＭ  Ｔｒｉｓ−ＨＣｌ（ｐＨ８．０）、０．１ｍＭ  ＥＤＴＡ）を添加し、該ＤＮＡを溶解させ、ＰＣＲ用ＤＮＡ試料を得ることができる。
【０１４１】
そして、リアルタイムＰＣＲにより、活性汚泥中に含まれる各種細菌の細菌数を定量することができる。
この場合、活性汚泥から抽出した核酸試料と、活性汚泥中に含まれる細菌に適したプライマー対およびプローブとを用いて、定量対象となる細菌の核酸の増幅に適したＰＣＲ条件（温度、時間、サイクル）で反応を行なうことにより、活性汚泥中に含まれる各種細菌および細菌数が定量される。前記プライマー対としては、例えば、アンモニア酸化細菌数定量用として、ＣＴＯ  １８９ｆＡ／Ｂ、ＣＴＯ  １８９ｆＣ、ＲＴ１ｒなど、亜硝酸酸化細菌の一種であるＮｉｔｒｏｓｐｉｒａ数定量用として、ＮＳＲ１１１３ｆ、ＮＳＲ１２６４ｒなどが挙げられる。また、前記プローブとしては、アンモニア酸化細菌定量用として、ＴＭＰ１、亜硝酸酸化細菌の一種であるＮｉｔｒｏｓｐｉｒａ数定量用として、ＮＳＲ１１４３Ｔａｑなどが挙げられる。
【０１４２】
具体的には、例えば、アンモニア酸化細菌、亜硝酸酸化細菌および真正細菌それぞれの細菌数は、アンモニア酸化細菌について、ＣＴＯ  １８９ｆＡ／Ｂ〔ＧＧＡＧＲＡＡＡＧＣＡＧＧＧＧＡＴＣＧ（配列番号：１）〕およびＣＴＯ  １８９ｆＣ〔ＧＧＡＧＧＡＡＡＧＴＡＧＧＧＧＡＴＣＧ（配列番号：２）〕をフォワードプライマー〔例えば、ＣＴＯ  １８９ｆＡ／Ｂ：ＣＴＯ  １８９ｆＣ＝２：１の混合物〕とし、ＲＴ１ｒ〔ＣＧＴＣＣＴＣＴＣＡＧＡＣＣＡＲＣＴＡＣＴＧ（配列番号：３）〕をリバースプライマーとするプライマー対と、ＴＭＰ１〔ＣＡＡＣＴＡＧＣＴＡＡＴＣＡＧＲＣＡＴＣＲＧＣＣＧＣＴＣ（配列番号：４）〕プローブとを有するプライマー／プローブセット；亜硝酸酸化細菌について、ＮＳＲ１１１３ｆ〔ＣＣＴＧＣＴＴＴＣＡＧＴＴＧＣＴＡＣＣＧ（配列番号：５）〕とＮＳＲ１２６４ｒ〔ＧＴＴＴＧＣＡＧＣＧＣＴＴＴＧＴＡＣＣＧ（配列番号：６）〕とのプライマー対と、ＮＳＲ１１４３Ｔａｑ〔ＡＧＣＡＣＴＣＴＧＡＡＡＧＧＡＣＴＧＣＣＣＡＧＧ（配列番号：７）〕のプローブとを有するプライマー／プローブセットを用いたＴａｑＭａｎ法により定量化されうる。
【０１４３】
なお、それぞれのプローブとしては、５’末端をＦＡＭ（６−ｃａｒｂｏｘｙｆｌｕｏｒｅｓｃｅｉｎ）、３’末端をＴＡＭＲＡ（６−ｃａｒｂｏｘｙｔｅｔｒａｍｅｔｈｙｌｒｈｏｄａｍｉｎｅ）で標識したものが挙げられる。
【０１４４】
また、硝酸還元細菌数の定量用プライマー対として、ｎａｒＨ５０Ｆ、ｎａｒＨｒ３Ｂなどが挙げられる。
具体的には、硝酸還元細菌定量用として、ｎａｒＨ５０Ｆ 〔ＡＡＲＴＧＹＡＴＣＧＧＹＴＧＣＣＡ（配列番号：８）〕をフォワードプライマーとし、ｎａｒＨｒ３Ｂ〔ＴＣＣＣＡＲＫＣＣＴＴＧＧＧＲＴＡＧ（配列番号：９）〕をリバースプライマーとするプライマー対としたＳＹＢＲ Ｇｒｅｅｎ法などにより定量化されうる。
【０１４５】
また、亜硝酸還元細菌数の定量用プライマー対として、ｎｉｒＫ８７６、ｎｉｒＫ１０４０およびｃｄ３ａＦ、Ｒ３ｃｄなどが挙げられる。
具体的には、亜硝酸還元細菌定量用として、ｎｉｒＫ８７６ 〔ＡＴＹＧＧＣＧＧＶＡＹＧＧＣＧＡ（配列番号：１０）〕をフォワードプライマーとし、ｎｉｒＫ１０４０ 〔ＧＣＣＴＣＧＡＴＣＡＧＲＴＴＲＴＧＧＴＴ（配列番号：１１）〕をリバースプライマーとするプライマー対、ｃｄ３ａＦ 〔ＧＴＳＡＡＣＧＴＳＡＡＧＧＡＲＡＣＳＧＧ（配列番号：１２）〕をフォワードプライマーとし、Ｒ３ｃｄ〔ＧＡＳＴＴＣＧＧＲＴＧＳＧＴＣＴＴＧＡ（配列番号：１３）〕をリバースプライマーとするプライマー対としたＳＹＢＲ Ｇｒｅｅｎ法により定量化されうる。
【０１４６】
前記リアルタイムＰＣＲには、濃度が１０7、１０6、１０5、１０4、１０3、または１０2コピー／５μｌ／１反応であるサンプルを用いて作成された検量線が用いられる。具体的には、前記検量線は、リアルタイムＰＣＲの対象となる遺伝子に対応する核酸をＰＣＲにより増幅し、その後、プラスミドベクターにクローニングした組み換えプラスミドを用いて作成されたものである。例えば、亜硝酸酸化細菌の場合、ＮＳＲ１１１３ｆとＮＳＲ１２６４ｒとのプライマー対を用いて、亜硝酸酸化細菌に特異的な１６Ｓ  ｒＲＮＡ遺伝子を増幅させ、得られた産物をプラスミドベクターにクローニングして得られた組み換えプラスミドを検量線作成用のスタンダードとして用いた。前記リアルタイムＰＣＲは、濃度（遺伝子数）が既知のスタンダードサンプルと、試料から精製したＤＮＡ（１ｎｇ、若しくは、１０ｎｇ）を鋳型として用いて行なわれる。通常、前記リアルタイムＰＣＲは、１６Ｓ  ｒＲＮＡ遺伝子や水質の浄化に係わる酵素をコードする遺伝子、例えば、アンモニア酸化細菌の１６Ｓ  ｒＲＮＡ遺伝子〔Ａｐｐｌ．Ｅｎｖｉｒｏｎ．Ｍｉｃｒｏｂｉｏｌ．、２００１年発行、第６７巻、第９７２頁〜第９７６頁〕、亜硝酸酸化細菌の一種であるニトロスピラの１６Ｓ  ｒＲＮＡ遺伝子〔Ｅｎｖｉｒｏｎ．Ｓｃｉ．Ｔｅｃｈｎｏｌ．、２００３年発行、第３７巻、第３４３頁〜第３５１頁〕、硝酸還元細菌のｎａｒＨ遺伝子〔Ｓｙｓｔｅｍ．Ａｐｐｌ．Ｍｉｃｒｏｂｉｏｌ．、２０００年発行、第２３巻、第４７頁〜第５７頁〕、亜硝酸還元細菌のｎｉｒＳ（シトクロムｃｄ１タイプの亜硝酸還元酵素）遺伝子〔ＦＥＭＳ Ｍｉｃｒｏｂｉｏｌｏｇｙ Ｅｃｏｌｏｇｙ ２００４年発行、第４９巻、第４０１頁〜第４１７頁〕、ｎｉｒＫ（銅含有タイプの亜硝酸還元酵素）遺伝子〔Ｊ．Ｍｉｃｒｏｂｉｏｌ．Ｍｅｔｈｏｄｓ．、２００４年発行、第５９巻、第３２７頁〜第３３５頁〕などをターゲットとして行なわれる。
【０１４７】
スタンダードサンプルそれぞれの濃度のＣｔ値（閾値とＰＣＲの増幅曲線が交わる点）を算出し、Ｃｔ値と濃度との関係から検量線を作成する。一方、試料ＤＮＡについてもＣｔ値を求め、スタンダードサンプルから作成した検量線に当てはめることにより、リアルタイムＰＣＲに使用した試料ＤＮＡ  １ｎｇ若しくは、１０ｎｇあたりの前記遺伝子のコピー数を求める。最終的に、例えば、活性汚泥１ｍｌあたりの細菌数は、下記式（３３）
【数４７】
0004278701.tif 000036
から導かれる。
【０１４８】
なお、リアルタイムＰＣＲのターゲットとする遺伝子は、細菌の種類により１個の細菌細胞に存在する数が異なるので、上記式（３３）から求めた細菌数を遺伝子数で割ることにより、正確な細菌数を算出することができる。
【０１４９】
リアルタイムＰＣＲ法により得られた細菌数から、下記式（３４）
【数４８】
0004278701.tif 000037
に基づき、細菌（アンモニア酸化細菌、亜硝酸酸化細菌、硝酸還元細菌、亜硝酸還元細菌）の濃度（ｍｇ−ＣＯＤCr／ｌ）を求める。
なお、前記１細菌あたりの乾燥重量には、測定値あるいは文献値〔文献値の例：０．２８ｐｇ／細胞、Ａｐｐｌ．Ｅｎｖｉｒｏｎ．Ｍｉｃｒｏｂｉｏｌ（２００２），６８，２４５−２５３）などの値が用いられる。
また、前記換算係数は、下記式（３５）
Ｃ5Ｈ7Ｏ2Ｎ＋５Ｏ2  →  ５ＣＯ2＋２Ｈ2Ｏ＋ＮＨ3  ・・・（３５）
から、細菌の酸素消費量（５×３２）を分子量（１１３）で割ることにより計算できる。この場合、前記の換算係数は、１．４１６である。
これを細胞１つあたりに換算すると、１細胞 （１細菌）＝０．２８ｐｇ×１．４１６ｍｇ−ＣＯＤCr／ｍｇ−ＳＳ＝３．９６５×１０-10ｍｇ−ＣＯＤCr＝３．９６５×１０-13ｇ−ＣＯＤCrとなる。
【０１５０】
通常、主としてアンモニアや硝酸等の窒素化合物が処理対象物質であり、他に有機性化合物を多く含有していない被処理水の生物処理の工程における細菌の反応については、前記式（１３）、（２１）、（２７）、（３２）に示す関数によって最大反応速度を定義し、パラメータとして用いることにより精度の高いシミュレーションを実施させることができるが、上記に説明した方法によって、より実際の生物処理工程に合致した関数を求めてシミュレーション装置に組み込むことでより精度の高いシミュレーションを実施させることができ、キャリブレーションの手間をよりいっそう削減させることができる。
【０１５１】
なお、このシミュレーション装置に組み込まれる設定項目としては、上記に説明した以外に、例えば、被処理水の溶存酸素（ｍｇ−Ｏ2／ｌ）、溶解性不活性有機物質（ｍｇ−ＣＯＤCr／ｌ）、易分解性有機物質（ｍｇ−ＣＯＤCr／ｌ）、アルカリ度（ｍｏｌｅＨＣＯ3／ｌ）、浮遊不活性有機物質（ｍｇ−ＣＯＤCr／ｌ）、遅分解性有機物質（ｍｇ−ＣＯＤCr／ｌ）、従属栄養細菌（ｍｇ−ＣＯＤCr／ｌ）、従属栄養細菌の細胞内貯蔵有機物質（ｍｇ−ＣＯＤCr／ｌ）、浮遊物質（ｍｇ−ＳＳ／ｌ）などが挙げられる。
【０１５２】
また、このシミュレーション装置は、適宜設定を変更して、他の生物学的処理におけるシミュレーションを実施させることも可能である。
【０１５３】
このようなシミュレーション装置でのシミュレーションを実施しつつ生物処理を実施することにより、予想外に処理対象物質を多く含んだ処理液を次段の処理に流下させてしまって正常な状態への復旧に手間取ったり、あるいは、過剰に時間を掛けて処理を実施してしまったりするなど生物処理に係る手間が必要以上に増大することを抑制させることができ、生物処理方法を効率の良いものとし得る。
したがって、生物処理装置の運転に要する手間も削減させることができ、運転コストの低減なども図ることができる。
【０１５４】
なお、このシミュレーション方法やシミュレーション装置、あるいは、これらを用いた生物処理方法ならびに生物処理装置は、硝化、脱窒にかかるもののみならず、処理対象物質が有機ハロゲン化合物であれば、デハロコッコイデス（Ｄｅｈａｌｏｃｏｃｃｏｉｄｅｓ）属細菌などの脱ハロゲン能を持つ細菌を用いた生物処理工程や、処理対象物質が硫酸などの硫黄化合物であれば、硫酸還元菌、例えば、デスルフォトマキュラム（Ｄｅｓｕｌｆｏｔｏｍａｃｕｌｕｍ）属細菌、デスルフォバクター（Ｄｅｓｕｌｆｏｂａｃｔｅｒ）属細菌、デスルフォバクテリウム（Ｄｅｓｕｌｆｏｂａｃｔｅｒｉｕｍ）属細菌、デスルフォヴィクリオ（Ｄｅｓｕｌｆｏｖｉｃｒｉｏ）属細菌などを用いた生物処理工程においても本実施形態おいて説明した硝化、脱窒における場合と同様の効果を期待できるものである。
【実施例】
【０１５５】
（実験例１）
主としてアンモニア及び硝酸が処理対象物質であり、他に有機性化合物を多く含有していない被処理水に対する硝化・脱窒の生物処理工程に用いられている汚泥を対象として最大反応速度の値と、生物処理工程において細菌１個あたりに一日に負荷される処理対象物質の量との相関関係調査を実施した。
【０１５６】
（アンモニア酸化細菌の反応速度）
アンモニア酸化細菌の反応速度は、以下のようにして測定を実施する。
５００ｍｌ容三角フラスコに、希釈水［１ｌあたりの組成：炭酸水素ナトリウム  ２４０ｍｇ、ＢＯＤ−Ａ液〔ＪＩＳ  Ｋ  ０１０２の２１の項に従う、緩衝液（ｐＨ７．２）〕１ｍｌ、ＢＯＤ−Ｂ液（ＪＩＳ  Ｋ  ０１０２の２１の項に従う、硫酸マグネシウム溶液）  １ｍｌ、ＢＯＤ−Ｃ液（ＪＩＳ  Ｋ  ０１０２の２１の項に従う、塩化カルシウム溶液）  １ｍｌ、ＢＯＤ−Ｄ液〔ＪＩＳ  Ｋ  ０１０２の２１の項に従う、塩化鉄（ＩＩＩ）溶液〕１ｍｌ、残部  水］  ３９０ｍｌを入れ、１０００ｍｇ−Ｎ／ｌの塩化アンモニウム水溶液  １０ｍｌを添加して、混合物Ａを調製する。ついで、５００ｍｌ容三角フラスコ内の混合物Ａを、恒温槽中、３０℃に維持し、撹拌しながら、１０分間以上曝気し、溶液Ａを得る。ここで、溶液ＡのｐＨを測定する。
【０１５７】
硝化汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、硝化汚泥試料を得る。
【０１５８】
その後、得られた硝化汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記溶液Ａと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、塩化アンモニウムに由来するアンモニア性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１５９】
また、前記硝化汚泥試料と前記溶液Ａとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、ＪＩＳ  Ｋ  ０１０２の１４の項に従い、汚泥濃度を測定する。
【０１６０】
塩化アンモニウムに由来するアンモニア性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１６１】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４２．５の項に従う）を行なうことにより、アンモニア性窒素の量を測定する。分析結果より単位時間あたりのアンモニア性窒素の量の変化を算出し、この変化量をアンモニア酸化の最大反応速度とする。
【０１６２】
上記測定方法により、汚泥濃度の異なる試料を約５０点測定し、最大反応速度の値と、生物処理工程においてアンモニア酸化細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を調査した。
得られた結果を、最大反応速度の値を縦軸に、アンモニア酸化細菌１個あたりに一日に負荷される処理対象物質の量を横軸としたグラフにプロットした。
結果を図２に示す。
【０１６３】
この図２からも、最大反応速度の値が定数ではないこと、ならびに、アンモニア酸化細菌１個あたりに一日に負荷される処理対象物質の量の増大にともなって最大反応速度の値が増大していることがわかる。
さらに、最大反応速度をＶAOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）とし、細菌１個あたりに一日に負荷されるアンモニア性窒素量をＬAOB（ｆｇ ＮＨ4−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）とした際に、アンモニア酸化細菌１個あたりに一日に負荷される処理対象物質の量：ＬAOBが１００〜４０００（ｆｇ ＮＨ4−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）の範囲においては、ＶAOBとＬAOBとが下記式（３６）
【数４９】
0004278701.tif 000038
に良く一致する関係にあることがわかる。
【０１６４】
（亜硝酸酸化細菌の反応速度）
亜硝酸酸化細菌の反応速度は、以下のようにして測定を実施する。
５００ｍｌ容三角フラスコに、前記希釈水  ３９０ｍｌを入れ、１０００ｍｇ−Ｎ／ｌの亜硝酸ナトリウム水溶液  １０ｍｌを添加して、混合物Ｂを調製する。ついで、５００ｍｌ容三角フラスコ内の混合物Ｂを、恒温槽中、３０℃に維持し、撹拌しながら、１０分間以上曝気し、溶液Ｂを得る。ここで、溶液ＢのｐＨを測定する。
【０１６５】
硝化汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、硝化汚泥試料を得る。
【０１６６】
その後、得られた硝化汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記溶液Ｂと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、亜硝酸ナトリウムに由来する亜硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１６７】
また、前記硝化汚泥試料と前記溶液Ｂとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１６８】
亜硝酸ナトリウムに由来する亜硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１６９】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．１．２の項に従う）を行なうことにより、亜硝酸性窒素の量を測定する。分析結果より単位時間あたりの亜硝酸性窒素の量の変化を算出し、これを亜硝酸酸化の最大反応速度として求める。
【０１７０】
上記測定法により、汚泥濃度の異なる試料を約５０点測定し、最大反応速度の値と、生物処理工程において亜硝酸酸化細菌１個あたりに一日に負荷される処理対象物質（アンモニア性窒素および亜硝酸性窒素）の量との相関関係を調査した。
得られた結果を、最大反応速度の値を縦軸に、亜硝酸酸化細菌１個あたりに一日に負荷される処理対象物質の量を横軸としたグラフにプロットした。
結果を図３に示す。
【０１７１】
この図３からも、最大反応速度の値が定数ではないこと、ならびに、亜硝酸酸化細菌１個あたりに一日に負荷される処理対象物質の量の増大にともなって最大反応速度の値が増大していることがわかる。
【０１７２】
さらに、最大反応速度をＶNOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）とし、亜硝酸酸化細菌１個あたりに一日に負荷されるアンモニア性窒素と亜硝酸性窒素との合計量をＬNOB（ｆｇ−Ｎ（ＮＨ4−Ｎ＋ＮＯ2−Ｎ）・ｃｏｐｙ-1・ｄａｙ-1）とした際に、亜硝酸酸化細菌１個あたりに一日に負荷される処理対象物質の量：ＬNOBが１０００〜６００００（fｇ・ｃｏｐｙ-1・ｄａｙ-1）の範囲においては、ＶNOBとＬNOBとが下記式（３７）
【数５０】
0004278701.tif 000039
に良く一致する関係にあることがわかる。
ここでＬNOBをアンモニア性窒素と亜硝酸性窒素との合計量としたのは、硝化槽内でアンモニア性窒素はアンモニア酸化細菌により亜硝酸性窒素に転換され、亜硝酸性窒素としては、通常、この合計量が負荷されるためである。
【０１７３】
（硝酸還元細菌の反応速度）
脱窒汚泥中に含まれる細菌による硝酸還元速度の測定（硝酸還元速度試験）は、以下のようにして実施する。
５００ｍｌ容三角フラスコに、前記希釈水３８０ｍｌを入れ、恒温槽中、３０℃に維持し、撹拌しながら、１０００ｍｇ−Ｎ／ｌの硝酸ナトリウム水溶液  １０ｍｌを添加し、ついで、メタノール水溶液〔５０００ｍｇメタノール／ｌ〕  １０ｍｌを添加して、混合物Ｃを調製する。その後、前記混合物Ｃに、散気球により窒素ガスを吹き込み、１０分間脱気する。ここで、前記混合物ＣのｐＨを測定する。前記５００ｍｌ容三角フラスコに、シリコン栓をし、撹拌しながら窒素ガスを吹き込み、該５００ｍｌ容三角フラスコ内の気相中の空気を除く。その後、１ｌ容テドラーバックに窒素ガスを吹き込み、前記５００ｍｌ容三角フラスコのシリコン栓に設置する。
【０１７４】
脱窒汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、脱窒汚泥試料を得る。
【０１７５】
その後、得られた脱窒汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記混合物Ｃと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、硝酸ナトリウムに由来する硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１７６】
また、前記脱窒汚泥試料と前記混合物Ｃとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１７７】
硝酸ナトリウムに由来する硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、汚泥濃度を測定する。
【０１７８】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．２．５の項に従う）を行なうことにより、硝酸性窒素の量を測定する。
分析結果より単位時間あたりの硝酸性窒素の量の変化を算出し、これを硝酸還元反応の最大反応速度として求めることができる。
【０１７９】
上記測定方法により、汚泥濃度の異なる試料を約５０点測定し、最大反応速度の値と、生物処理工程において硝酸還元細菌１個あたりに一日に負荷される処理対象物質量との相関関係を調査した。
得られた結果を、最大反応速度の値を縦軸に、硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素の量を横軸としたグラフにプロットした。
結果を図４に示す。
【０１８０】
この図４からも、最大反応速度の値が定数ではないこと、ならびに、硝酸還元細菌１個あたりに一日に負荷される処理対象物質の量の増大にともなって最大反応速度の値が増大していることがわかる。
【０１８１】
さらに、最大反応速度をＶNARB（ｆｇ・ｃｏｐｙ-1・ｈ-1）とし、硝酸還元細菌１個あたりに一日に負荷される全窒素量をＬNARB（ｆｇＮ・ｃｏｐｙ-1・ｄａｙ-1）とした際に、硝酸還元細菌１個あたりに一日に負荷される処理対象物質の量：ＬNARBが５〜７０（fｇ・ｃｏｐｙ-1・ｄａｙ-1）の範囲においては、ＶNARBとＬNARBとが下記式（３８）
【数５１】
0004278701.tif 000040
に良く一致する関係にあることがわかる。
【０１８２】
ここで、実施例では流入する被処理水中の全窒素量はアンモニア性窒素量と亜硝酸性窒素量と硝酸性窒素量の合計量にほぼ等しいことと、硝化槽においてほぼ全てのアンモニア性窒素および亜硝酸性窒素が硝酸性窒素に転換されていたことから、硝酸還元細菌１個あたりに負荷される窒素量は被処理水中の全窒素量（Ｔ−Ｎ）に等しいとして計算を行った。
仮に、全窒素量がアンモニア性窒素量と亜硝酸性窒素量と硝酸性窒素量の合計量に等しくない場合、または、硝化槽内でアンモニア性窒素および亜硝酸性窒素が硝酸性窒素まで完全に硝化されない場合は、脱窒槽に流入する硝酸性窒素量をＬNARBとすることが望ましい。
【０１８３】
（亜硝酸還元細菌の反応速度）
脱窒汚泥中に含まれる細菌による亜硝酸還元速度の測定（亜硝酸還元速度試験）は、以下のように測定する。
【０１８４】
５００ｍｌ容三角フラスコに、前記希釈水３８０ｍｌを入れ、恒温槽中、３０℃に維持し、撹拌しながら、１０００ｍｇ−Ｎ／ｌの亜硝酸ナトリウム水溶液  １０ｍｌを添加し、ついで、メタノール水溶液〔５０００ｍｇメタノール／ｌ〕  １０ｍｌを添加して、混合物Ｄを調製する。その後、前記混合物Ｄに、散気球により窒素ガスを吹き込み、１０分間脱気する。ここで、前記混合物ＤのｐＨを測定する。前記５００ｍｌ容三角フラスコに、シリコン栓をし、撹拌しながら窒素ガスを吹き込み、該５００ｍｌ容三角フラスコ内の気相中の空気を除く。その後、１ｌ容テドラーバックに窒素ガスを吹き込み、前記５００ｍｌ容三角フラスコのシリコン栓に設置する。
【０１８５】
脱窒汚泥を５００ｍｌ採取し、遠心分離器によって固液分離する。得られた産物から、上澄み液を除去し、前記希釈水  ５０ｍｌに撹拌し、分散させる。得られた産物の全容量を、前記希釈水を用いて、１００ｍｌとなるように調整し、脱窒汚泥試料を得る。
【０１８６】
その後、得られた脱窒汚泥試料  １００ｍｌを、前記５００ｍｌ容フラスコに入れ、前記混合物Ｄと混合する。混合と同時に、５ｍｌ容シリンジで、前記混合により得られた試料  ５ｍｌをサンプリングし、フィルター（アドヴァンテック社製、商品名：ガラス繊維ろ紙ＧＦ−７５、孔径：０．３μｍ）で濾過する。なお、亜硝酸ナトリウムに由来する亜硝酸性窒素が消費されるまで（約２時間）、一定時間毎に、サンプリングを行なう。
【０１８７】
また、前記脱窒汚泥試料と前記混合物Ｄとの混合直後に、前記サンプリングとは別に、３０ｍｌの試料を採取し、汚泥濃度を測定する。
【０１８８】
亜硝酸ナトリウムに由来する亜硝酸性窒素が消費された後、残りの混合物のｐＨを測定すると共に、活性汚泥濃度を測定する。
【０１８９】
各サンプリング時点での試料を、イオンクロマトグラフ分析法（ＪＩＳ  Ｋ  ０１０２の４３．１．２の項に従う）を行なうことにより、亜硝酸性窒素の量を測定する。分析結果より単位時間あたりの亜硝酸性窒素の量の変化を算出し、これを亜硝酸還元反応の最大反応速度として求める。
【０１９０】
上記測定方法により、汚泥濃度の異なる試料を約５０点測定し、最大反応速度の値と、生物処理工程において亜硝酸還元細菌１個あたりに一日に負荷される処理対象物質量との相関関係を調査した。
得られた結果を、最大反応速度の値を縦軸に、亜硝酸還元細菌１個あたりに一日に負荷される処理対象物質の量を横軸としたグラフにプロットした。
結果を図５に示す。
【０１９１】
この図５からも、最大反応速度の値が定数ではないこと、ならびに、亜硝酸還元細菌１個あたりに一日に負荷される処理対象物質の量の増大にともなって最大反応速度の値が増大していることがわかる。
【０１９２】
さらに、最大反応速度をＶNIRB（ｆｇ・ｃｏｐｙ-1・ｈ-1）とし、亜硝酸還元細菌１個あたりに一日に負荷される全窒素量をＬNIRB（ｆｇＮ・ｃｏｐｙ-1・ｄａｙ-1）とした際に、亜硝酸還元細菌１個あたりに一日に負荷される処理対象物質の量：ＬNIRBが５〜１２０（fｇ・ｃｏｐｙ-1・ｄａｙ-1）の範囲においては、ＶNIRBとＬNIRBとが下記式（３９）
【数５２】
0004278701.tif 000041
に良く一致する関係にあることがわかる。
【０１９３】
ここで、実施例では流入する被処理水中の全窒素量はアンモニア性窒素量と亜硝酸性窒素量と硝酸性窒素量の合計量にほぼ等しいことと、硝化槽においてほぼ全てのアンモニア性窒素が硝酸性窒素または亜硝酸性窒素に転換されていたことから、亜硝酸還元細菌１個あたりに負荷される窒素量は被処理水中の全窒素量（Ｔ−Ｎ）に等しいとしてＬNIRBの算出を行った。
仮に、全窒素量がアンモニア性窒素量と亜硝酸性窒素量と硝酸性窒素量の合計量に等しくない場合、または、硝化槽内でアンモニアが完全に硝化されない場合は、脱窒槽に流入する硝酸性窒素と亜硝酸性窒素の合計量をＬNIRBとすることが望ましい。
【０１９４】
（実験例２）
前記実験例１よりも、高負荷な状態で、しかも、実設備を模擬させた状態で最大反応速度と細菌１個あたりに一日に負荷される処理対象物質の量との関係を調査した。
【０１９５】
（設備）
用いた実験設備の概要を図６に示す。
実験設備は、有効容積３５リットルの反応槽と、該反応槽で生物処理させる被処理水（原水）を調整する原水タンク、該原水タンクから被処理水を反応槽に供給するための原水ポンプ、反応槽から処理水を排出するための処理水ポンプ、反応槽中において散気を行うためのエアポンプ（ブロワ）、脱窒工程時に反応槽にメタノールを供給するためのメタノールポンプ、反応槽中の槽内水のｐＨ調整を行うｐＨ調整ポンプ等によって構成されている。
また、この実験設備は、各機器の運転が制御盤によりタイマーや設定値等で制御され、運転記録がデータロガーに記録されるように構成されている。
【０１９６】
（生物処理工程の設定）
前記反応槽の運転状態を下記（１）〜（４）のように切り替えることにより硝化工程、脱窒工程、酸化工程、沈殿工程のサイクル運転を行った。
尚、反応槽への原水の投入量は、１４リットル／１サイクルとし、５６リットル／日とした。
（１）硝化工程（２時間）
・活性汚泥含有液が所定容量入っている反応槽に、原水の投入を開始すると同時に、ブロワによる散気を行い、硝化工程を開始する。
・ブロワは反応槽内に設置した溶存酸素濃度計（ＤＯ計）の指示値と連動させ、溶存酸素濃度（ＤＯ）が３．０±０．１ｍｇ／ｌの範囲内で運転できるよう自動でＯＮ−ＯＦＦさせる。
・硝化によるｐＨの低下を避けるため、反応槽内に設置したｐＨ計と連動させ、ｐＨが７．７以上を維持できるよう、自動で２％苛性ソーダをｐＨ調整ポンプで添加する。
・反応槽の温度センサで水温を測定し、ヒーターで３０℃前後を維持できるよう加温する（この反応槽内の温度維持については、脱窒工程、酸化工程も同じ）。
（２）脱窒工程（２時間）
・ブロワによる散気を停止し、反応槽内を撹拌機により混合する。
・それと同時に５％に調整したメタノールをポンプにより注入する。
（３）酸化工程（１時間）
・酸化工程の開始と同時にブロワによる散気を再開し、残留するメタノールを処理する。
・硝化工程と同様にＤＯを３．０±０．１ｍｇ／ｌでコントロールする。
（４）沈殿工程（１時間）
・酸化工程終了時に、撹拌機、ブロワ、ヒーター、苛性ソーダの添加を停止して反応槽内を静置し、固液分離を行う。
・４５分後に処理水排出ポンプを稼働させて、指定水位まで処理水を排出する。
この（１）〜（４）の工程（６時間／１サイクル）を１日４サイクル実施、設定の負荷の窒素含有排水を処理した。
【０１９７】
（被処理水の調整）
反応槽に供給する被処理水（原水）は以下の試薬を用いて、窒素やリンの組成が火力発電所排水処理設備を模擬した濃度になるよう調整を行い、窒素負荷を変動させつつ反応槽に供給し処理後の処理水に残留する窒素の量を測定した。
（使用試薬）
ＮＨ4−Ｎ：  ３０  ｍｇ／ｌ（硫酸アンモニウム）
ＮＯ3−Ｎ：  ２０  ｍｇ／ｌ（硝酸ナトリウム）
ＰＯ4−Ｐ：  ２    ｍｇ／ｌ（リン酸二水素ナトリウム）
ＮａＨＣＯ3：３６０ｍｇ／ｌ（炭酸水素ナトリウム）
【０１９８】
なお、窒素負荷の設定を増加する場合は、原水と同様の基質の比率で作成した濃縮基質を作成し、原水投入と並行して濃縮基質を反応槽内に投入した。
濃縮基質の投入量が原水の投入量の１／５０以下になるよう基質濃度を調整し、負荷上昇時にも反応槽内の水力学的滞留時間がほとんど変化しないよう配慮した。
この反応槽への負荷の与え方をどのように変化させたかを図７に示す。
反応槽における容積負荷は、０．０８ｇ−Ｎ／ｌ／ｄを基準状態とし、一時的にこの基準負荷に対して２倍の負荷、ならびに２．５倍の負荷となるように変動させた。
【０１９９】
この負荷変動によって、沈殿工程時に排出される処理水中の窒素成分がどのように変化するかを調査した。
結果を、図８に示す。
この図８にも示されているように、処理水に残留する各態窒素（ＮＨ4−Ｎ、ＮＯ2−Ｎ、ＮＯ3−Ｎ）の濃度は運転期間を通して１ｍｇ／ｌ以下となっており、反応槽に流入させた窒素のほぼ全量が脱窒処理されていることがわかった。
【０２００】
この期間における反応槽内の活性汚泥中の菌数の変化を調査した結果を、図９、図１０に示す。
また、槽内水のＭＬＳＳ濃度の変動状況について調査した結果を図１１に示す。
これらの図からもわかるように、それぞれの細菌は運転期間中で約１オーダーの範囲で大きく変動していた。
ＭＬＳＳの変動幅が２０００〜５０００ｍｇ／ｌの範囲であったのに比べると、その中に含まれる細菌が大きく変動していたことがわかる。
【０２０１】
この菌数で窒素負荷を除して、細菌１個あたりに一日に負荷される窒素の量：Ｌ（fｇ・ｃｏｐｙ-1・ｄａｙ-1）と、最大反応速度を：Ｖ（ｆｇ・ｃｏｐｙ-1・ｈ-1）とを求めた。
結果を、図１２〜図１５に示す。図１２は、アンモニア酸化細菌について、図１３は、亜硝酸酸化細菌について、図１４は硝酸還元細菌について、そして、図１５は、亜硝酸還元細菌について細菌１個あたりに一日に負荷される窒素の量：Ｌ（fｇ・ｃｏｐｙ-1・ｄａｙ-1）と、最大反応速度を：Ｖ（ｆｇ・ｃｏｐｙ-1・ｈ-1）とを求めた結果を示すものである。
このように、実験例１と同様に、高い負荷においても、細菌１個あたりに一日に負荷される処理対象物質の量の増大にともなって最大反応速度の値が増大していることがわかる。
また、その増大する傾向が下記式（４０）
【数５３】
0004278701.tif 000042
（ただし、ａ、ｂ、及びｃは、それぞれ独立した定数である。）
となる形式の関数に良く一致することがわかる。
したがって、この式（４０）に示すような関数を設けてシミュレーションを行うことでその予測精度を向上させ得ることがわかる。
【０２０２】
また、先の実験例１の結果（図２〜図５）と、この図１２〜図１５に示す結果と、さらに、有機化合物が含まれた窒素含有排水を処理している排水処理槽から採取した汚泥について同様の実験（実験例３）を行った結果を、図１６〜図１９に示す。
【０２０３】
この図１６からも、ＶAOBとＬAOBとを関係付ける関数については、アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量：ＬAOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が１．０×１０2≦ＬAOB≦３．５×１０4において、最大反応速度：ＶAOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）−２．５×１０3｝≦ＶAOB≦｛４．０×１０3・ＬAOB／（１．０×１０4＋ＬAOB）＋２．５×１０3｝となるように定められることでシミュレーションの予測精度を向上させ得ることがわかる。
また、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得ることがわかる。
【０２０４】
また、図１７からは、ＶNOBとＬNOBとを関係付ける関数について、亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNOB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が１．０×１０3≦ＬNOB≦１．２×１０6となる範囲において、最大反応速度：ＶNOB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）−１．０×１０5｝≦ＶNOB≦｛２．５×１０5・ＬNOB／（４．０×１０5＋ＬNOB）＋１．０×１０5｝となるように定められることでシミュレーションの予測精度を向上させ得ることがわかる。
また、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得ることがわかる。
【０２０５】
また、図１８からは、ＶNARBとＬNARBとを関係付ける関数について、硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素量：ＬNARB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が５．０≦ＬNARB≦２５００となる範囲において、最大反応速度：ＶNARB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）−１．７×１０2｝≦ＶNARB≦｛２．２×１０2・ＬNARB／（７．０×１０2＋ＬNARB）＋７０｝となるように定められることでシミュレーションの予測精度を向上させ得ることがわかる。
また、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得ることがわかる。
【０２０６】
さらに、図１９からは、ＶNIRBとＬNIRBとを関係付ける関数については、亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量：ＬNIRB（ｆｇ・ｃｏｐｙ-1・ｄａｙ-1）が５．０≦ＬNIRB≦３．５×１０3となる範囲において、最大反応速度：ＶNIRB（ｆｇ・ｃｏｐｙ-1・ｈ-1）が｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）−２．５×１０2｝≦ＶNIRB≦｛７．０×１０2・ＬNIRB／（２．５×１０3＋ＬNIRB）＋２．５×１０2｝となるように定められることでシミュレーションの予測精度を向上させ得ることがわかる。
また、下水なども含めて窒素成分が処理対象物質として含有されている被処理水の生物学的処理全般に当該関数を適応し得ることがわかる。
【０２０７】
すなわち、生物処理工程に硝化工程を有する場合は、下記式（４１）
【数５４】
0004278701.tif 000043
ただし、
ＶAOB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬAOB：アンモニア酸化細菌１個あたりに一日に負荷されるアンモニア性窒素量（ｆｇ  ＮＨ4−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
に示す関数によってアンモニア酸化の最大反応速度を定義するとともに、下記式（４２）
【数５５】
0004278701.tif 000044
ただし、
ＶNOB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNOB：亜硝酸酸化細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｆｇ ＮＯ2−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
に示す関数によって亜硝酸酸化の最大反応速度を定義してシミュレーション装置に組み込み、生物処理工程にさらに脱窒工程を有する場合は、前記式（４１）に示す関数、前記式（４２）に示す関数に加え、下記式（４３）
【数５６】
0004278701.tif 000045
ただし、
ＶNARB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNARB：硝酸還元細菌１個あたりに一日に負荷される硝酸性窒素量（ｆｇ ＮＯ3−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
に示す関数によって硝酸還元の最大反応速度を定義し、下記式（４４）
【数５７】
0004278701.tif 000046
ただし、
ＶNIRB：最大反応速度（ｆｇ・ｃｏｐｙ-1・ｈ-1）
ＬNIRB：亜硝酸還元細菌１個あたりに一日に負荷される亜硝酸性窒素量（ｆｇ ＮＯ2−Ｎ・ｃｏｐｙ-1・ｄａｙ-1）
に示す関数によって亜硝酸還元の最大反応速度を定義してシミュレーション装置に組み込むことで、より精度の高いシミュレーションを実施させることができる。
【０２０８】
これら結果をもとに、従来のＩＷＡのＡＳＭ３（活性汚泥モデルＮｏ．３）計算プログラムの硝化工程及び脱窒工程に関する最大反応速度に代えて、上記式（４１）乃至（４４）に示す、細菌１個あたりに一日に負荷される窒素の量によって変化する最大反応速度の関係式を組み込んで、本発明の活性汚泥モデル（新規モデル）の実証を行った。
【０２０９】
この新規モデルに基づくＡＳＭ計算用のＰｅｔｅｒｓｏｎ’ｓマトリクスを表１（表は、一部のみ）に示す。
【０２１０】
【表１】
0004278701.tif 000047
【０２１１】
このマトリクスは、図２０の比較表に示すとおり最大反応速度が一定の従来モデルと相違している。
これら新規モデル及び従来モデル（ＡＳＭ３）を用いて、窒素負荷を基準負荷から２．５倍に増加させた、直後と、負荷上昇後の処理量の履歴の影響を十分受けていると考えられる４８時間後の処理水の水質についてシミュレートし、実測データとの比較を行った。
結果を、図２１（負荷増加直後）、図２２（負荷増加後４８時間）に示す。
図中のプロットは、実験の分析値、実線は細菌１個あたりに一日に負荷される窒素の量によって変化する最大反応速度の関係式を組み込んだシミュレータ（新規モデル）により計算した結果、破線は従来のＡＳＭ３にて計算した結果をそれぞれ示す。
実験の水質分析値は、負荷上昇直後は硝化工程でアンモニアが消費されるのに約１２０分、脱窒工程で硝酸が処理されるのにも１２０分要していたが、負荷上昇後４８時間が経過すると、アンモニアの消費には約８０分、硝酸の消費は約６０分で終了しており、処理速度の向上が確認された。
亜硝酸の生成は硝化脱窒を含めほとんど認められなかった。
新規モデルのシミュレータによる計算結果は、負荷上昇直後および負荷上昇後４８時間の計算結果を良好に再現できており、また、水質予測精度は、目標値である分析値±２０％を達成した。
一方、従来のＡＳＭ３の計算結果は、硝化工程は４８時間の間の硝化菌の増殖により処理速度の向上を再現できているが、脱窒工程は従属栄養微生物の増殖のみでは処理速度の向上が再現できておらず、活性の変化を再現できていないものと考えられた。
以上のように、本発明によればシミュレーションの予測精度を向上させ得ることがわかる。
【０２１２】
（実験例４）
（細菌１個あたりに一日に処理された処理対象物質量適用事例）
実験例１乃至３において採取されたデータを基にして、“細菌１個あたりに一日に処理された処理対象物質（窒素）の量”を計算して求め、図１６〜図１９と同様にグラフ化した結果を、図２３〜図２６に示す。
なお、実験例１〜３のように、“細菌１個あたりに一日に負荷される処理対象物質（窒素）の量”をファクターとして利用すれば、処理水質の予測計算を行う場合に、生物処理の結果を反映させる必要がないため計算速度を早めることが出来る。
一方で、この実験例４のように“細菌１個あたりに一日に処理された処理対象物質（窒素）の量”をファクターとして利用する場合には、生物処理後の処理水の水質を測定し、その測定に基づいて“細菌１個あたりに一日に処理された処理対象物質（窒素）の量”を算出して、結果を最大反応速度に反映させるべくフィードバックさせる必要があることから結果の予測に手間を要することとなる。
しかし、負荷と処理量とにおいて差が生じやすい（処理水に処理対象物質を残留させやすい）場合などにおいては、“細菌１個あたりに一日に処理された処理対象物質の量”をファクターとして利用する方が、精度の高いシミュレーション結果を得られやすいという利点がある。
なお、通常の生物処理工程においては、図８や、図１６〜図１９と図２３〜図２６との対比に示されているように、細菌１個あたりの一日の負荷と処理量とには大きな差がないため、計算速度を向上させる点において“細菌１個あたりに一日に負荷される処理対象物質の量”をファクターとして利用することが好適である。
【０２１３】
（実験例５）
（塩化物イオンの影響調査）
実験例１にて作製した溶液Ａ、Ｂ、Ｃ、Ｄにそれぞれ、２０００ｍｇ／ｌ、４０００ｍｇ／ｌの濃度となるように塩化物イオン濃度を加えた以外は、実験例１の最大反応速度の測定と同様に最大反応速度を測定した。
そして、塩化物イオンの影響がない状態での最大反応速度の値をそれぞれ、ＶAOB0、ＶNOB0、ＶNARB0、ＶNIRB0とし、これらを１００とした場合の最大反応速度の値（ＶAOB、ＶNOB、ＶNARB0、ＶNIRB）を塩化物イオン２０００ｍｇ／ｌ及び４０００ｍｇ／ｌ加えた場合について求めた。
結果をそれぞれ、図２７〜図３０に示す。
【０２１４】
この図２７からも、アンモニア酸化細菌によるアンモニア酸化反応においては、（１＋４．４×１０-5・ＤCL）≦ｋAOB（ＤCL）≦（１＋１．６４×１０-4・ＤCL）となる関数：ｋAOB（ＤCL）（ただし、「ＤCL」は塩化物イオン濃度（ｍｇ／ｌ）を表す）を設定して、最大反応速度：ＶAOBと、細菌１個あたりに単位時間に負荷されるアンモニア性窒素の量：ＬAOBとのあいだに、ＶAOB＝ｆAOB（ＬAOB）・ｋAOB（ＤCL）の関数を設定することにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させうることがわかる。
【０２１５】
また、図２８からは、亜硝酸酸化細菌による亜硝酸酸化反応においては、（１−８．７×１０-5・ＤCL）≦ｋNOB（ＤCL）≦（１−４．０×１０-5・ＤCL）となる関数：ｋNOB（ＤCL）（ただし、「ＤCL」は塩化物イオン濃度（ｍｇ／ｌ）を表す）を設定して、最大反応速度：ＶNOBと、細菌１個あたりに単位時間に負荷される亜硝酸性窒素の量：ＬNOBとのあいだに、ＶNOB＝ｆNOB（ＬNOB）・ｋNOB（ＤCL）の関数を設定することにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させうることがわかる。
【０２１６】
また、図２９からは、硝酸還元細菌による硝酸還元反応においては、（１−７．９×１０-5・ＤCL）≦ｋNARB（ＤCL）≦（１−１．０×１０-5・ＤCL）となる関数：ｋNARB（ＤCL）（ただし、「ＤCL」は塩化物イオン濃度（ｍｇ／ｌ）を表す）を設定して、最大反応速度：ＶNARBと、細菌１個あたりに単位時間に負荷される硝酸性窒素の量：ＬNARBとのあいだに、ＶNARB＝ｆNARB（ＬNARB）・ｋNARB（ＤCL）の関数を設定することにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させうることがわかる。
【０２１７】
さらに、図３０からは、亜硝酸還元細菌による亜硝酸還元反応においては、（１−６．０×１０-5・ＤCL）≦ｋNIRB（ＤCL）≦（１−５．０×１０-6・ＤCL）となる関数：ｋNIRB（ＤCL）（ただし、「ＤCL」は塩化物イオン濃度（ｍｇ／ｌ）を表す）を設定して、最大反応速度：ＶNIRBと、細菌１個あたりに単位時間に負荷される亜硝酸性窒素の量：ＬNIRBとのあいだに、ＶNIRB＝ｆNIRB（ＬNIRB）・ｋNIRB（ＤCL）の関数を設定することにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させうることがわかる。
【０２１８】
なお、この図２７から図３０に示された実験例５の結果については、実験例１のように“細菌１個あたりに一日に負荷される処理対象物質の量”をファクターとして利用する場合のみならず、実験例４のように“細菌１個あたりに一日に処理された処理対象物質の量”をファクターとして利用する場合にも同様に反映させうるものである。
すなわち、アンモニア酸化細菌１個あたりに単位時間に処理されるアンモニア性窒素の量を“ＬAOB”、亜硝酸酸化細菌１個あたりに単位時間に処理される亜硝酸性窒素の量を“ＬNOB”、硝酸還元細菌１個あたりに単位時間に処理される硝酸性窒素の量を“ＬNARB”、及び亜硝酸還元細菌１個あたりに単位時間に処理される亜硝酸性窒素の量を“ＬNIRB”としたときにおいても、（１＋４．４×１０-5・ＤCL）≦ｋAOB（ＤCL）≦（１＋１．６４×１０-4・ＤCL）となる関数：ｋAOB（ＤCL）、（１−８．７×１０-5・ＤCL）≦ｋNOB（ＤCL）≦（１−４．０×１０-5・ＤCL）となる関数：ｋNOB（ＤCL）、（１−７．９×１０-5・ＤCL）≦ｋNARB（ＤCL）≦（１−１．０×１０-5・ＤCL）となる関数：ｋNARB（ＤCL）、及び（１−６．０×１０-5・ＤCL）≦ｋNIRB（ＤCL）≦（１−５．０×１０-6・ＤCL）となる関数：ｋNIRB（ＤCL）（ただし、「ＤCL」は塩化物イオン濃度（ｍｇ／ｌ）を表す）を設定して、最大反応速度：ＶAOB、ＶNOB、ＶNARB、及びＶNIRBとの間に、ＶAOB＝ｆAOB（ＬAOB）・ｋAOB（ＤCL）、ＶNOB＝ｆNOB（ＬNOB）・ｋNOB（ＤCL）、ＶNARB＝ｆNARB（ＬNARB）・ｋNARB（ＤCL）、及びＶNIRB＝ｆNIRB（ＬNIRB）・ｋNIRB（ＤCL）の関数を設定することにより、処理水の水質の予測精度を従来のシミュレーションに比べて向上させうる。
【０２１９】
以上のように、本発明によれば予測精度の低下を抑制しつつキャリブレーションの手間を削減させ得るシミュレーション方法などを実施させうることがわかる。
【図面の簡単な説明】
【０２２０】
【図１】生物処理装置を示す概略ブロック図。
【図２】実験例１におけるアンモニア酸化細菌のアンモニア酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図３】実験例１における亜硝酸酸化細菌の亜硝酸酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図４】実験例１における硝酸還元細菌の硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図５】実験例１における亜硝酸還元細菌の亜硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図６】生物処理装置をモデル化した実験設備を示す概略図。
【図７】窒素の負荷変動を示すグラフ。
【図８】生物処理工程（実験後）の処理水質を示すグラフ。
【図９】実験期間における反応槽中の菌数変化を示すグラフ。
【図１０】実験期間における反応槽中の菌数変化を示すグラフ。
【図１１】実験期間における反応槽中のＭＬＳＳ変化を示すグラフ。
【図１２】実験例２におけるアンモニア酸化細菌のアンモニア酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１３】実験例２における亜硝酸酸化細菌の亜硝酸酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１４】実験例２における硝酸還元細菌の硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１５】実験例２における亜硝酸還元細菌の亜硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１６】実験例1〜３を網羅したアンモニア酸化細菌のアンモニア酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１７】実験例1〜３を網羅した亜硝酸酸化細菌の亜硝酸酸化反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１８】実験例1〜３を網羅した硝酸還元細菌の硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図１９】実験例1〜３を網羅した亜硝酸還元細菌の亜硝酸還元反応における最大反応速度と、細菌１個あたりに一日に負荷される処理対象物質の量との相関関係を示すグラフ。
【図２０】本発明のシミュレーションと従来のシミュレーションとにおける相違を示す比較表。
【図２１】本発明のシミュレーションと従来シミュレーションとによるシミュレーション結果比較グラフ。
【図２２】本発明のシミュレーションと従来シミュレーションとによるシミュレーション結果比較グラフ。
【図２３】アンモニア酸化細菌のアンモニア酸化反応における最大反応速度と、細菌１個あたりに一日に処理された処理対象物質の量との相関関係（実験例４の結果）を示すグラフ。
【図２４】亜硝酸酸化細菌の亜硝酸酸化反応における最大反応速度と、細菌１個あたりに一日に処理された処理対象物質の量との相関関係（実験例４の結果）を示すグラフ。
【図２５】硝酸還元細菌の硝酸還元反応における最大反応速度と、細菌１個あたりに一日に処理された処理対象物質の量との相関関係（実験例４の結果）を示すグラフ。
【図２６】亜硝酸還元細菌の亜硝酸還元反応における最大反応速度と、細菌１個あたりに一日に処理された処理対象物質の量との相関関係（実験例４の結果）を示すグラフ。
【図２７】塩化物イオン濃度によるアンモニア酸化細菌の最大反応速度変化（実験例５の結果）を示すグラフ。
【図２８】塩化物イオン濃度による亜硝酸酸化細菌の最大反応速度変化（実験例５の結果）を示すグラフ。
【図２９】塩化物イオン濃度による硝酸還元細菌の最大反応速度変化（実験例５の結果）を示すグラフ。
【図３０】塩化物イオン濃度による亜硝酸還元細菌の最大反応速度変化（実験例５の結果）を示すグラフ。
【符号の説明】
【０２２１】
２：硝化槽、４：脱窒槽、６：再曝気槽、８：沈殿槽
【要約】
【課題】  本発明は、予測精度の低下を抑制しつつキャリブレーションの手間を削減させ得るシミュレーション方法ならびにシミュレーション装置の提供、ならびに、要する手間を削減し得る生物処理方法や生物処理装置の提供を課題としている。
【解決手段】  細菌による処理対象物質の分解反応における最大反応速度の値と、生物処理工程において細菌１個あたりに単位時間に負荷される前記処理対象物質の量または細菌１個あたりに単位時間に処理された前記処理対象物質の量との間に所定の関数関係を有する状態でパラメータに用いることで課題を解決する。
【選択図】図１【図１】
0004278701.tif 000052
【図２】
0004278701.tif 000053
【図３】
0004278701.tif 000054
【図４】
0004278701.tif 000055
【図５】
0004278701.tif 000056
【図６】
0004278701.tif 000057
【図７】
0004278701.tif 000058
【図８】
0004278701.tif 000059
【図９】
0004278701.tif 000060
【図１０】
0004278701.tif 000061
【図１１】
0004278701.tif 000062
【図１２】
0004278701.tif 000063
【図１３】
0004278701.tif 000064
【図１４】
0004278701.tif 000065
【図１５】
0004278701.tif 000066
【図１６】
0004278701.tif 000067
【図１７】
0004278701.tif 000068
【図１８】
0004278701.tif 000069
【図１９】
0004278701.tif 000070
【図２０】
0004278701.tif 000071
【図２１】
0004278701.tif 000072
【図２２】
0004278701.tif 000073
【図２３】
0004278701.tif 000074
【図２４】
0004278701.tif 000075
【図２５】
0004278701.tif 000076
【図２６】
0004278701.tif 000077
【図２７】
0004278701.tif 000078
【図２８】
0004278701.tif 000079
【図２９】
0004278701.tif 000080
【図３０】
0004278701.tif 000081
【配列表】
0004278701000001.app