WO2010055642A1

WO2010055642A1 - 水質シミュレーション方法及び装置

Info

Publication number: WO2010055642A1
Application number: PCT/JP2009/005993
Authority: WO
Inventors: 加藤文隆; 加藤敏朗; 味埜俊; 佐藤弘泰; 庄司仁
Original assignee: Nippon Steel Corp; University of Tokyo NUC
Current assignee: Nippon Steel Corp; University of Tokyo NUC
Priority date: 2008-11-14
Filing date: 2009-11-10
Publication date: 2010-05-20
Anticipated expiration: 2011-05-14
Also published as: CN102209688B; JPWO2010055642A1; BRPI0921825A2; JP5435658B2; CN102209688A; KR101233623B1; KR20110074577A

Abstract

　本発明は、生物反応槽に流入する廃水中の化合物成分の各成分濃度を分析する分析工程と、前記各成分濃度の分析値を各成分のＣＯＤ濃度に換算するＣＯＤ換算工程と、化学量論パラメーター及び反応速度式パラメーターを設定するパラメーター設定工程と、前記生物反応槽の溶存酸素濃度を測定する溶存酸素濃度測定工程と、前記各成分のＣＯＤ濃度、パラメーター及び溶存酸素濃度を用いて、前記生物学的反応槽内で前記生物学的好気処理した後の処理水における各成分のＣＯＤ濃度を算出する計算工程と、当該算出された生物学的処理水の各成分のＣＯＤ濃度を、各成分濃度に換算する成分濃度換算工程とを、含む水質シミュレーション方法及び装置を提供する。

Description

水質シミュレーション方法及び装置

　本発明は、産業廃水や工場廃水を含む廃水の処理に関し、より詳しくは、生物分解性のある化合物成分を含む廃水の生物学的好気処理プロセスにおける水質シミュレーション方法及び装置に関する。　本願は、２００８年１１月１４日に、日本国に出願された特願２００８－２９２５１９号に基づき優先権を主張し、その内容をここに援用する。

　近年の閉鎖性水域における水質規制は強化傾向にあり、下水道、産業廃水、工場廃水における対策が急務である。水質規制値を遵守するためには、どの成分をどこまで処理すれば基準を満たせるかを予測することが必要である。これまで下水道分野では、活性汚泥法等の生物学的処理プロセスによる下水の処理方法が取られてきた。その運転管理は主に管理者の経験に基づいて行われていることが多く、安定した処理水質を得ることは困難であった。たとえば、流入水質が変動した場合に、変更すべき操作条件、操作量は下水処理場ごとに異なっていた。

　そこで、管理者の経験に依存しない水質予測及び運転支援ツールとして、ＩＷＡ（国際水協会）の提唱する活性汚泥モデル（Ａｃｔｉｖａｔｅｄ　Ｓｌｕｄｇｅ　Ｍｏｄｅｌ；ＡＳＭ）が提案されている。活性汚泥モデルは、大きく次の手順から構成されている：
　（１）流入廃水のＣＯＤ濃度を、溶解性不活性有機物、易分解性有機物、浮遊不活性有機物、遅分解性有機物等の性質によって分画し、それぞれＣＯＤ濃度ベースの変数と設定する；
　（２）従属栄養生物の増殖や自己分解等のプロセスごとに、変数間の化学量論及びプロセスの反応速度式を設定する；
　（３）化学量論係数及び反応速度定数のパラメーターを、酸素消費速度試験又は水質の実測データからのキャリブレーションにより決定する；
　（４）シミュレーションを実行し、生物処理槽、処理水のＣＯＤ濃度等が算出される。

　この活性汚泥モデルは管理ツールとして提案されており、活性汚泥モデルを用いた下水処理管理システムが提案されている（以下、特許文献１参照）。
　一方で、下水と異なる成分をもつ産業廃水や工場廃水等についても、生物学的処理プロセスによる産業廃水や工場廃水の処理が行われている。しかしながら、これらの廃水に対して、水質シミュレーション方法に活性汚泥モデルが適用された事例はない。
　特に、製鉄所コークス工場から発生する廃水は安水と呼ばれ、フェノール、チオ硫酸、チオシアン酸を主成分とするが、これらは下水には含まれない成分であり、活性汚泥モデルの適用性について検討した事例は見られない。

特開２００３－３０００９３公報

味埜俊、活性汚泥モデル、日本、株式会社環境新聞社、２００５年1月３１日Ｊ．Ｓ．Ｃｅｃｈ、Ｊ．Ｃｈｕｄｏｂａ　ａｎｄ　Ｐ．Ｇｒａｕ、Ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　Ｋｉｎｅｔｉｃ　Ｃｏｎｓｔａｎｔｓ　ｏｆ　Ａｃｔｉｖａｔｅｄ　Ｓｌｕｄｇｅ　Ｍｉｃｒｏｏｒｇａｎｉｓｍｓ、Ｗａｔｅｒ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ、Ｖｏｌ．１７、ｐｐ．２５９－２７２、１９８４

　前記したように、生物学的好気処理における水質シミュレーション方法として活性汚泥モデルが一般的であるが、産業廃水や工場廃水の生物学的好気処理に適用された事例はない。この原因として、以下の２点が挙げられる。
　（１）活性汚泥モデルでは、下水のような複数成分が混合した廃水を生物学的好気処理する際の流入廃水と処理水のＣＯＤ濃度のみ予測するが、生物分解性を有する化合物成分の濃度は予測できないことが挙げられる。例えば、コークス製造過程で発生する安水では、安水中に含まれるフェノールが廃水基準項目の一つであり、フェノール濃度をシミュレーションにより予測する必要があるが、活性汚泥モデルでは特に安水のような複数成分が混合した廃水を生物学的好気処理した後の処理水中フェノール濃度を予測することが困難であるため、産業廃水や工場廃水に活性汚泥モデルを適用する意義は小さかった。
　（２）産業廃水や工場廃水には、溶解性遅分解性の有機物（界面活性剤、例えば、直鎖アルキルベンゼンスルホン酸等）、または、無機物（例えば、チオ硫酸、チオシアン酸）が含まれる。しかしながら、活性汚泥モデルで対象としている都市下水中には分解速度の非常に遅い溶解性遅分解性成分はそれほど含まれておらず、モデルの酸素収支に与える影響を考慮していなかった。すなわち、そもそも溶解性遅分解性成分の概念が活性汚泥モデルには含まれていなかったため、産業廃水や工場廃水に活性汚泥モデルを適用することは困難であった。

　そこで本発明は、生物分解性を有する化合物成分を含む産業廃水や工場廃水などの廃水を、生物反応槽内で生物学的好気処理するプロセスにおいて、新たな活性汚泥モデルを構築し、廃水中の当該化合物の濃度が生物学的好気処理後に、どのように変化するかを推定することができる水質シミュレーション方法及び装置を提供することを目的とする。

　本発明は、廃水中の成分ごとに分画し、成分ごとに成分を分解する細菌の種類および変数および反応プロセスを設定することで、活性汚泥モデルを利用した生物学的好気処理プロセスのシミュレーションが可能となり、処理水中の成分濃度を求め、さらに、溶解性遅分解性成分を含む廃水であっても、廃水中の成分ごとに分画するため、溶解性遅分解性成分という概念を用いること無しに、活性汚泥モデルを利用した生物学的好気処理プロセスのシミュレーションが可能となる、生物分解性のある化合物成分を含む産業廃水や工場廃水に適用可能な水質シミュレーション方法を提供することを目的とする。

　本発明は、具体的には、以下の［１］～［６］である：
　［１］生物分解性を有する化合物成分を含む廃水を生物反応槽内で生物学的好気処理するプロセスにおける水質シミュレーション方法であって、
　前記生物反応槽に流入する前記廃水中の前記化合物成分の各成分濃度を分析する分析工程と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、各成分濃度の分析値を各成分のＣＯＤ濃度に換算するＣＯＤ換算工程と、
　化学量論パラメーターである増殖収率、反応速度式パラメーターである飽和定数、最大比増殖速度、並びに前記化合物成分を分解する細菌の種類および濃度を設定するパラメーター設定工程と、
　前記生物反応槽の溶存酸素濃度を測定する溶存酸素濃度測定工程と、
　前記各成分のＣＯＤ濃度、前記増殖収率、前記飽和定数、前記最大比増殖速度、並びに前記化合物成分を分解する細菌の種類及び濃度、並びに溶存酸素濃度を用いて、下記演算式（１）：

｛式中、Ｃ_iは、各成分濃度であり、ｉは、各成分を表す通し番号であり、Ｐ_ijは、化学量論パラメーターであり、ｊは、各プロセスを表す通し番号であり、そしてρ_jは、反応速度式（反応速度式パラメーターを含む速度式）である。｝により、前記生物学的反応槽内で前記生物学的好気処理した後の処理水における各成分のＣＯＤ濃度を算出する計算工程と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、当該算出された生物学的処理水の各成分のＣＯＤ濃度を、各成分濃度に換算する成分濃度換算工程とを、含む水質シミュレーション方法。

　［２］前記パラメーター設定工程は、
（ア）前記流入する廃水及び処理水において、事前に時系列的に別途採取したそれぞれの生物分解性を有する化合物成分濃度と溶存酸素濃度を用いて、キャリブレーションによって決定する方法、
（イ）前記各成分それぞれに対して、各成分単独で前記生物学的好気処理を行って連続的に溶存酸素濃度を測定し、当該測定値から酸素消費速度を算出し、当該算出された酸素消費速度のデータから決定する方法、のいずれかの方法を用いて、前記化学量論パラメーター及び前記反応速度パラメーターを設定する、前記［１］に記載の水質シミュレーション方法。

　［３］上記廃水が、コークス製造工程で発生する安水であり、上記化合物成分のうち、有機物成分はフェノール、並びに、無機物成分はチオ硫酸及びチオシアン酸であり、チオ硫酸、チオシアン酸のシミュレーションについては、上記ＣＯＤ換算工程の代わりにチオ硫酸及びチオシアン酸濃度を硫黄換算する硫黄換算工程と、
　上記生物学的反応槽で上記生物学的好気処理した後の処理水における各成分の硫黄濃度を算出する硫黄計算工程と、
　当該算出された生物学的処理水の各成分の硫黄濃度を、チオ硫酸濃度及びチオシアン酸濃度に再換算する硫黄再換算工程と
　を有する、前記［１］又は［２］に記載の水質シミュレーション方法。

　［４］生物分解性を有する化合物成分を含む廃水を生物反応槽内で生物学的好気処理するプロセスで用いる水質シミュレーション装置であって、
　前記生物反応槽に流入する前記廃水中の前記化合物成分の各成分濃度を分析する分析手段と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、各成分濃度の分析値を各成分のＣＯＤ濃度に換算するＣＯＤ換算手段と、
　化学量論パラメーターである増殖収率、反応速度式パラメーターである飽和定数、最大比増殖速度、並びに前記化合物成分を分解する細菌の種類および濃度を設定するパラメーター設定手段と、
　前記生物反応槽の溶存酸素濃度を測定する溶存酸素濃度測定手段と、
　前記各成分のＣＯＤ濃度、前記増殖収率、前記飽和定数、前記最大比増殖速度、並びに前記化合物成分を分解する細菌の濃度及び溶存酸素濃度を用いて、下記演算式（１）：

｛式中、Ｃ_iは、各成分濃度であり、ｉは、各成分を表す通し番号であり、Ｐ_ijは、化学量論パラメーターであり、ｊは、各プロセスを表す通し番号であり、そしてρ_jは、反応速度式（反応速度式パラメーターを含む速度式）である。｝により、前記生物学的反応槽で前記生物学的好気処理した後の処理水における各成分のＣＯＤ濃度を算出する計算手段と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、当該算出された生物学的処理水の各成分のＣＯＤ濃度を、各成分濃度に換算する成分濃度換算手段とを、含む水質シミュレーション装置。

　［５］前記パラメーター設定手段は、
（ア）前記流入する廃水及び処理水において、事前に時系列的に採取したそれぞれの生物分解性を有する化合物成分濃度と溶存酸素濃度を用いて、キャリブレーションによって決定する手段、
（イ）前記各成分それぞれに対して、各成分単独で前記生物学的好気処理を行って連続的に溶存酸素濃度を測定し、当該測定値から酸素消費速度を算出し、当該算出された酸素消費速度のデータから決定する手段、のいずれかの手段を用いて、前記化学量論パラメーター及び前記反応速度式パラメーターを設定する、前記［４］に記載の水質シミュレーション装置。

　［６］上記廃水が、コークス製造工程で発生する安水であり、上記化合物成分のうち、有機物成分はフェノール、並びに、無機物成分はチオ硫酸及びチオシアン酸であり、チオ硫酸、チオシアン酸のシミュレーションについては、上記ＣＯＤ換算手段の代わりにチオ硫酸及びチオシアン酸濃度を硫黄換算する硫黄換算手段と、
　上記生物学的反応槽で上記生物学的好気処理した後の処理水における各成分の硫黄濃度を算出する硫黄計算手段と、
　当該算出された生物学的処理水の各成分の硫黄濃度を、チオ硫酸濃度及びチオシアン酸濃度に再換算する硫黄再換算手段と
　を有する、前記［４］又は［５］に記載の水質シミュレーション装置。

　本発明の奏する効果は以下のとおりである。
　生物分解性を有する化合物成分を含む産業廃水や工場廃水に対して、流入廃水中の各成分濃度をＣＯＤ濃度に変換し、活性汚泥モデルを適用することで、処理水の各成分濃度が予測できる。また、活性汚泥モデルでは対象外であった溶解性遅分解性成分についても、その成分を分解する細菌の種類を設定することで、処理水の成分濃度を計算することが可能である。また、各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係を求める際には、ＣＯＤ_Mnの分析結果を用いることができるので、ＣＯＤ_Crのように危険かつ有害性のある薬品を用いた分析を行わなくてもよい。

本発明に係る、廃水の生物学的処理プロセスの水質シミュレーション方法を示す構成図である。フェノール濃度とＣＯＤ_Mnの相関関係を示す図である。チオ硫酸濃度とＣＯＤ_Mnの相関関係を示す図である。チオシアン酸濃度とＣＯＤ_Mnの相関関係を示す図である。酸素消費速度試験装置の構成図である。本発明に係る実施例１におけるシミュレーション結果を示す図である。本発明に係る実施例２におけるシミュレーション結果を示す図である。本発明に係る実施例３におけるシミュレーション結果を示す図である。

廃水に含まれる成分は水に溶解しない浮遊成分と、溶解性成分に分けられる。これは、例えばろ紙を用いて分離を行うことで、ろ紙に残ったものを浮遊成分、ろ液を溶解性成分として得られる。本発明では、生物分解と関係が深い溶解性成分の生物分解を対象としたシミュレーション方法である。浮遊成分はほとんど生物分解されず、最終的には沈降分離によって除去されるため、本発明では対象外とする。
溶解性成分は、生物分解性成分及び生物分解性がほとんど無い難分解性成分に分けられる。難分解性成分は、もともと生物分解が困難であり、濃度としても小さい。さらに、生物分解性成分は、分析による定性、定量が困難な未知成分が含まれるが、廃水中の割合として少量なため無視しても構わない。例えば安水では、フェノール、チオ硫酸、チオシアン酸以外の成分がこれに相当する。以上のことから、廃水における難分解性成分、及び、未知成分は無視しても、本発明の生物分解性を有する化合物成分の予測は可能である。
　図１は、本発明による廃水の水質シミュレーション方法のフローを例示した図である。また、図中の各「工程」は各「手段」に替えることができ、廃水の水質シミュレーション装置を例示した図でもある。

　ここでは廃水として、製鉄所におけるコークス工場から発生する安水を例として示すが、本発明の適用廃水は安水に限定されるものではなく、食品工業、化学工業、薬品工業、塗装工業、繊維工業、染色工業等の産業廃水や工場廃水に対しても、生物分解性を有する化合物成分を含有している限り適用可能である。例えば、安水中には、生物分解性を有する有機物成分としてフェノール、無機物成分としてチオ硫酸及びチオシアン酸を含んでいる。

　また、図１の処理水とは、連続処理の場合は生物学的好気処理後の流出水、バッチ処理の場合は所定時間経過後の槽内水のことである。
　また、生物学的好気処理とは、基質と細菌が存在し、ばっ気装置を有する反応槽において、ばっ気を行うことで溶存酸素を供給しながら、混合・撹拌を行うことで基質と細菌を接触・反応させることで、基質を分解処理する方法である。

　図１に示すように、本発明による水質シミュレーション方法１は、生物学的好気処理プロセスに流入する産業廃水や工場廃水の成分濃度を分析する分析工程２と、各成分濃度とＣＯＤ濃度との相関関係３を用いて各成分をＣＯＤ濃度９に換算するＣＯＤ換算工程４と、化学量論パラメーター（増殖収率）及び反応速度式パラメーター（飽和定数及び最大比増殖速度）、並びに化合物成分を分解する細菌の種類及び濃度のパラメーター５を設定するパラメーター設定工程６と、生物反応槽の溶存酸素濃度７を測定する溶存酸素濃度測定工程８と、パラメーター５、溶存酸素濃度７及びＣＯＤ濃度９を用いて、前述の演算式（１）により、生物学的処理水の各成分のＣＯＤ濃度１１を算出する計算工程１０と、各成分濃度とＣＯＤ濃度との相関関係３を用いて各成分のＣＯＤ濃度１１を各成分濃度に再換算することで、処理水の各成分濃度１３を算出する成分濃度換算工程１２とを備えている。

　分析工程２は、成分濃度を分析する方法について指定しないが、毎回同じ分析方法を用いることが望ましい。たとえば、フェノール濃度についてはＪＩＳ　Ｋ０１０２吸光光度法、チオ硫酸、チオシアン酸濃度についてはイオンクロマト分析、ＭＬＳＳ濃度については前記した非特許文献１等により分析できる。

相関関係３は、各成分濃度とＣＯＤ濃度から得られるものである。例として、図２Ａ～図２Ｃにフェノール濃度、チオ硫酸、チオシアン酸とＣＯＤ_Mn濃度との実測値に基づく相関関係を示す。その結果、フェノール、チオ硫酸、チオシアン酸１ｍｇあたりの酸素消費量は、それぞれ２．１ｍｇ、０．４８ｍｇ、１．０ｍｇであった。また、相関関係を求める際には、分析工程でＣＯＤ濃度を測定した方法と同じＣＯＤ_Ｃｒ又はＣＯＤ_Ｍｎのいずれの分析方法を用いることが望ましく、分析方法は一貫して同じ方法を用いる。

　このように、既知成分濃度と、ＣＯＤ濃度との相関関係を求める際には、ＣＯＤ_Ｃｒを用いなくとも、ＣＯＤ_Ｍｎの分析結果を用いることができるので、ＣＯＤ_Ｃｒのように危険かつ有害性のある薬品を用いた分析を行わなくても良いという利点がある。
　また、この他にも、各成分濃度とＣＯＤ濃度の相関関係を簡易に得る方法として、成分の理論的酸素消費量を化学反応式から算出する方法もある。たとえば、フェノール、チオ硫酸、チオシアン酸１ｍｇあたりの理論的酸素消費量は、以下の化学反応式（１）～（３）：

を用いて、それぞれ、２．３８ｍｇ、０．５７ｍｇ、１．１ｍｇと算出される。

　ＣＯＤ換算工程４は、分析工程２で得られた各成分濃度を、相関関係３を用いることで、ＣＯＤ濃度に換算するものである。ここで、ＣＯＤ濃度に換算するのは、活性汚泥モデルにおいて、最も優れた尺度として採用されているためである。すなわち、前記した非特許文献１にあるように、有機物や生物あるいは利用された酸素といった一連の電子等量の関係をＣＯＤだけで表現できるとともに、酸素収支をとることができるといった利点があるためである。ＣＯＤ濃度に換算せずに各成分濃度で計算することも可能ではあるが、上記の利点を活用できないため、現実的ではない。

　チオ硫酸、チオシアン酸のような硫黄化合物は、ここで硫黄換算してシミュレーションを行い、後述する成分濃度換算工程１２で各成分濃度に再変換することも可能である。これにより、チオ硫酸、チオシアン酸を硫黄成分として分けて考えることで、ＣＯＤとしてカウントされない硫酸イオン成分等が算出可能となる。これにより、硫黄の物質収支を追跡することが可能となり、硫黄成分の挙動が把握できる。但し、硫黄換算することで硫黄の物質収支を取る場合であっても、酸素の物質収支は取れているので、シミュレーションに影響することはない。また、硫黄化合物を分解する細菌は、一般的には硫黄酸化細菌と考えられており、それらは有機物を分解する従属栄養細菌と異なり、互いに作用を及ぼしあうことがないと考えられるため、硫黄化合物をその他の成分と分けて考えることにより、精度が高まる。その場合、チオ硫酸、チオシアン酸１ｍｇの硫黄換算値は、前記化学反応式（２）～（３）より、それぞれ０．５７ｍｇ、０．５５ｍｇと算出される。

　また、成分濃度からＣＯＤへの換算は、理論的酸素消費量又はＣＯＤ_Cr又はＣＯＤ_Mnのいずれか一つを一貫して用いることで、後述する成分濃度換算工程１２と同じとする。
　パラメーター５としては、各成分の化学量論パラメーターである増殖収率及び反応速度式パラメーターである飽和定数、最大比増殖速度並びに前記化合物成分を分解する細菌の種類及び濃度を用いる。
　また、生物学的処理において生物自身の自己分解が与える影響を加味するため、反応速度パラメーターに好気条件下における生物の比内生呼吸速度を設定しても良いが、通常は影響が少ないので、設定しなくてもよい。

　パラメーター設定工程６は、大きく分けて（Ａ）実験を行って、ある条件での実験データを解析する方法と、（Ｂ）文献値を拾う方法が挙げられ、前者によりパラメーター設定を行う方が望ましい。
　上記（Ａ）の方法について以下に述べる。この方法は、（ア）前記流入する廃水及び処理水において、事前に時系列的に別途採取したそれぞれの生物分解性を有する化合物成分濃度と溶存酸素濃度を用いて、キャリブレーションによって決定する方法、（イ）前記各成分それぞれに対して、各成分単独で前記生物学的好気処理を行って連続的に溶存酸素濃度を測定し、当該測定値から酸素消費速度を算出し、当該算出された酸素消費速度のデータから決定する方法のうち（ア）又は（イ）のいずれかの方法を用いて、事前にパラメーターを設定する。

　上記（ア）の時系列実測値データについて、具体的なデータ取得方法は指定しないが、濃度変化がわかるような時系列データを採取する。例えば、１週間の処理水水質をシミュレーションするのであれば、その間の３点以上、たとえば、２日、４日、６日経過時のデータを取得することが望ましい。　また、時系列実測値データは安定したものが望ましく、たとえば上記の場合、流入する安水の含有成分濃度、組成などが１週間の間で大きく変動しているデータは望ましくない。変動の大きさの基準は、たとえば、シミュレーション予測の精度を±５％以内に設定する場合、シミュレーション生物処理時間（滞留時間）の時間間隔で流入する安水の含有成分濃度の変動幅が５％以内であるものを用いることができる。

　上記（ア）のキャリブレーションによって決定する方法は、特に指定しないが、例えば、市販のシミュレーションソフトＡＱＵＡＳＩＭによって、時系列実測値、計算値及び時系列実測値の標準偏差から、下記演算式（２）：

｛式中、χ²（ｐ）は、対象とするモデルパラメーターｐのχ²の値であり、ｙ_meas,mは、ｍ番目の時系列データの実測値（成分濃度または酸素消費速度）であり、ｙ_m（ｐ）は、モデルパラメーターｐの値を仮定した時の、ｍ番目の時系列データの計算値（成分濃度または酸素消費速度）であり、σ_meas,mは、ｍ番目の時系列データの実測値の標準偏差又は実測値全体の標準偏差であり、そしてｎは、データポイントの数である。｝によりχ²の値を求め、χ²の値が最小となったときのパラメーターの値を用いる方法がある。パラメーターの値を求める方法は特に指定しないが、シンプレックス法、モンテカルロ法、遺伝的アルゴリズム等の数値解析手法によるものが望ましい。

　上記（イ）の酸素消費速度試験の酸素消費速度データからパラメーターを決定する方法は、図３の試験装置を用いる。図３で、酸素消費速度試験装置２１に、生物学的処理プロセスで用いる微生物汚泥２２、アンモニア性窒素が硝酸性窒素に硝化される際の酸素消費を抑制するため硝化阻害剤２３を添加し、撹拌装置２４で酸素消費速度試験装置２１の液を混合する。栄養塩の不足による微生物の活性の低下を防ぐため、栄養塩２５を添加してもよい。次に、酸素消費速度試験装置２１に対象とする成分２６を添加し、溶存酸素濃度計２７によって溶存酸素濃度の経時変化を測定する。溶存酸素濃度の記録はデータ記録装置２８で行ってもよい。測定中は溶存酸素濃度を一定以上に制御するため、制御値を下回った時に空気供給装置２９で空気を供給してもよい。ｐＨを一定に制御するため、ｐＨ計３０でｐＨを測定し、酸・アルカリ供給装置３１で酸又はアルカリを供給し、制御してもよい。酸素消費速度試験装置２１の温度を制御するため、ヒーター３２を備えた恒温水槽３３を用いてもよい。

　対象とする成分（以下、成分Ａと表記する）を添加して得られた溶存酸素濃度及び酸素消費速度の経時変化からパラメーターを決定する方法は特に指定しないが、たとえば増殖収率については、成分ＡのＣＯＤ換算濃度と、酸素消費量から、下記演算式（３）：

｛式中、Ｙ_Aは、成分Ａを分解する細菌の増殖収率であり、Ｓ_A,CODは、成分ＡのＣＯＤ換算濃度であり、そしてＯ_2measは、酸素消費量の実測値である。｝によって求めることができる。

　増殖収率をより厳密に求める方法には、内生呼吸による酸素消費速度を事前に求め、上記酸素消費量の実測値から内生呼吸による酸素消費速度を差し引くことで、増殖に使われた酸素消費速度を求めることができるが、本発明では特にその方法を指定しない。

　最大比増殖速度、飽和定数の値を決定する方法については特に指定しないが、例えば、前記した非特許文献２にある方法で求めることができる。すなわち、酸素消費速度試験で用いる微生物汚泥を1時間、ばっ気した後に、微生物汚泥を様々な希釈率で廃水と混合し、酸素消費速度を測定する方法である。測定中は溶存酸素濃度を一定に保つように、ばっ気による制御を行う。次に、測定された酸素消費速度を汚泥濃度で除すことで比酸素消費速度を求め、比酸素消費速度から内生呼吸による酸素消費速度を引いたもの（比基質酸化速度）を求める。比内生呼吸速度とは、生物自身の自己分解時の内生呼吸による酸素消費速度であり、その求め方は、酸素消費速度試験において、成分Ａを添加せずに行った時の酸素消費速度から求められる。この時、比増殖速度と比基質酸化速度の間には、下記演算式（４）：

｛式中、μ_Aは、成分Ａを分解する細菌の比増殖速度であり、Ｙ_Aは、成分Ａを分解する細菌の増殖収率であり、そしてｒ_OXは、比基質酸化速度である。｝に示す関係が成り立つ。

　一方で、比増殖速度μ_Aと最大比増殖速度μ_A,maxの間には、下記演算式（５）：

｛式中、μ_A,maxは、成分Ａを分解する細菌の最大比増殖速度であり、Ｋ_Aは、成分ＡのＣＯＤ濃度に対する飽和定数であり、Ｓ_A,CODは、成分ＡのＣＯＤ換算濃度であり、Ｓ_O2は、溶存酸素濃度であり、そしてＫ_O2は、溶存酸素濃度に対する飽和定数である。｝に示すミカエリス－メンテン式が成り立つ。

　ここで、酸素消費速度試験では溶存酸素濃度を高濃度に維持するため、Ｓ_O2／（Ｋ_O2＋Ｓ_O2）の値はほぼ１と見なせる。したがって、比増殖速度は成分濃度のみの関数とみなせる。以上より、演算式（４）で得られた比増殖速度μ_Aの値と成分濃度とのプロットが得られ、そのプロットに演算式（５）の式をフィッティングすることで、微生物の比最大増殖速度及び成分Ａに対する飽和定数が求まる。

　化合物成分を分解する細菌の種類および濃度を決定する方法については、特に指定しないが、種類については、基本的に一成分に対し一種類の細菌がその分解に対応することが望ましい。さらに、実験データや文献の知見から、複数成分を1種類の細菌が分解する可能性がある場合、実際の現象に即した計算予測や計算予測精度の向上のため、複数成分に1種類の細菌が対応するよう設定してもよい。但し、一成分を複数種類の細菌が分解することもあり得るが、それらは同じような分解作用を示すことが多く、また、細菌の種類を増やすことで計算量が増加するため、一成分に対し一種類の細菌を対応させることが望ましい。細菌の濃度については、化合物成分を用いた酸素消費速度が、下記演算式（６）：

｛式中、ＯＵＲは、酸素消費速度であり、そしてＸ_Aは、成分Ａを分解する細菌の濃度であり、そしてＹ_Aは、成分Ａを分解する細菌の増殖収率であり、Ｋ_Ａは、成分Ａ濃度のＣＯＤ換算値に対する飽和定数であり、Ｓ_Aは、成分ＡのＣＯＤ換算濃度であり、Ｋ_O2は、溶存酸素濃度に対する飽和定数であり、Ｓ_O2は、溶存酸素濃度である。｝で表されることが多く、この式を用いる方法がある。

　すなわち、演算式（６）で、最大比増殖速度、成分Ａに対する飽和定数は、前記の方法で決定され、酸素消費速度試験では溶存酸素濃度を高濃度に維持するため、Ｓ_O2／（Ｋ_O2＋Ｓ_O2）の値はほぼ１と見なせる。このことから、酸素消費速度は成分Ａを分解する細菌の濃度のみの関数として表される。したがって、演算式（６）に、酸素消費速度試験で得られるＯＵＲ実測値を代入することで、Ｘ_Ａを決定することができる。

　上記（Ｂ）の文献値を拾う方法については、例えば、前記した非特許文献１に典型的なパラメーター値が記述されており、好気条件下での従属栄養生物の収率は０．６３とされている。この値をパラメーター値として設定してもよいが、対象とする廃水によって値が異なる可能性があるため、なるべく上記（Ａ）の方法により設定するほうが望ましい。

　溶存酸素濃度７は、生物学的処理プロセスに溶存酸素濃度計を設置し、生物反応槽の溶存酸素濃度を測定した値を用いる方法、又は、パラメーター設定工程６で溶存酸素濃度をパラメーターとして扱い、キャリブレーションにより決定した値を用いる方法のいずれでもよい。

　計算工程１０は、ＣＯＤ濃度９、溶存酸素濃度７及びパラメーター５を用いて、処理水の各成分のＣＯＤ濃度１１を求める。計算に用いるモデルは、ＩＷＡ活性汚泥モデルをベースとした前述の演算式（１）により、処理水の各成分のＣＯＤ濃度１１を算出する。処理水の各成分のＣＯＤ濃度１１は成分濃度換算工程１２で、相関関係３を用いることで、処理水の各成分の成分濃度値１３に再変換することができる。

さらに、計算工程１０では、ｐＨ、アルカリ度、アンモニア濃度等の影響を含めたシミュレーションを行ってもよい。その計算方法の指定はしないが、例えば、Ｃ_iに水素イオン濃度、溶存二酸化炭素濃度、アンモニア濃度等を追加し、Ｐ_ijに水素イオン濃度、溶存二酸化炭素濃度、アンモニア濃度等の化学量論パラメーターを追加し、ρ_jに反応速度式を追加して、計算を行ってもよい。また、生物学的処理において生物自身の自己分解が与える影響を加味することで予測精度を上げるため、反応速度パラメーターに好気条件下における生物の比内生呼吸速度を設定してもよいが、通常は影響が少ないので、設定しなくてもよい。設定する場合は、下記演算式（７）：

｛式中、Ｘ_Ａは成分Ａを分解する細菌の濃度であり、そしてＰ_Ａは成分Ａを分解するプロセスにおける化学量論パラメーターであり、そしてρ_Ａは成分Ａを分解するプロセスにおける反応速度論パラメーターであり、そしてｂはＸ_Ａが自己分解するときの比内生呼吸速度定数である。｝によって求めることができる。

実施例１：フェノール、チオ硫酸、チオシアン酸を対象としたバッチ試験のシミュレーション（１成分を１種類の細菌が分解する場合）
　以下、本発明の実施の形態を図に基づいて説明する。図１は、生物学的処理プロセスの水質シミュレーション方法の構成と計算の流れを示したブロック図である。
　以下、生物分解性のある化合物成分としてフェノール、チオ硫酸、チオシアン酸を対象としたバッチ試験のシミュレーション方法について説明する。バッチ試験は、１Ｌの反応容器に、フェノール濃度１００ｍｇ／Ｌ、チオ硫酸濃度１００ｍｇ／Ｌ、チオシアン酸濃度１０ｍｇ／Ｌとなるよう成分を添加し、ＭＬＳＳ濃度は５０００ｍｇ／Ｌで行った。試験中は、溶存酸素濃度計により溶存酸素濃度を測定し、ｐＨ計によりｐＨを測定した。また、溶存酸素濃度を３．２５ｍｇ／Ｌ、ｐＨ７．５に制御した。

　分析工程２は、生物学的処理プロセスに流入する廃水中の成分濃度を分析する工程であるが、ここでは上記のバッチ試験装置中の初期成分濃度であり、各成分の添加量から算出するとともに、フェノール濃度についてはＪＩＳ　Ｋ０１０２吸光光度法、チオ硫酸、チオシアン酸濃度についてはイオンクロマト分析、ＭＬＳＳ濃度については前記した非特許文献１の遠心分離法により確認を行った。ＣＯＤ換算工程４は、分析工程２で得られたフェノール、チオ硫酸、チオシアン酸の濃度をＣＯＤ濃度に換算するＣＯＤ換算工程である。ＣＯＤ換算には、図２Ａ～図２Ｃに示すフェノール濃度、チオ硫酸濃度及びチオシアン酸濃度とＣＯＤ_Mn濃度との相関関係３により換算を行い、フェノール、チオ硫酸、チオシアン酸の濃度を、処理水の各成分のＣＯＤ濃度１１に換算した。

　パラメーター設定工程６は、計算工程１０で用いるパラメーター設定を行った。
　パラメーター５は、増殖収率、最大比増殖速度、飽和定数、及び後述する計算工程１０で用いるフェノール分解細菌濃度、チオ硫酸分解細菌濃度及びチオシアン酸分解細菌濃度を設定した。設定方法は、酸素消費速度試験の酸素消費速度データから決定する方法、及び、キャリブレーションによる方法を検討した結果、いずれも妥当な値が得られたが、今回のケースでは酸素消費速度試験から得られたパラメーター値に多少のばらつきが見られたため、キャリブレーションによる方法を用いた。キャリブレーションには市販のシミュレーションソフトＡＱＵＡＳＩＭを用いた。具体的には、実測値、計算値及び実測値の標準偏差から、演算式（２）に示されるχ²の値を求め、χ²の値が最小となったときのパラメーターの値を算出した。実測値は、事前に行った成分ごとのバッチ試験から、成分濃度の経時変化を得た。計算値は、後述する計算工程１０の方法により求めた。

　溶存酸素濃度７は、バッチ試験の制御値３．２５ｍｇ／Ｌを用いた。
　計算工程１０は、廃水の各成分のＣＯＤ濃度９、溶存酸素濃度７及びパラメーター５を用いて、バッチ試験のフェノール、チオ硫酸、チオシアン酸濃度の経時変化のシミュレーションを行い、処理水の各成分のＣＯＤ濃度１１を求めた。計算方法は、前記演算式（１）を用いた。以下の表１に、計算に用いる各成分濃度Ｃ_i及び化学量論パラメーター、反応速度パラメーターを示す。

　以下の表２は、プロセスρ_jにおける化学量論の関係を表しており、例えば、プロセスρ₁のフェノール分解細菌の増殖では、フェノール分解細菌（Ｘ_Phe）に対してＰ₁₁＝－１／Ｙ_Pheの割合でフェノール濃度が変化し、Ｐ₄₁＝－（１－Ｙ_Phe）／Ｙ_Pheの割合でＳ_O2が変化することを表す。

　以下の表３は、計算に用いる反応速度式を示す。

　例えば、プロセスρ₁のフェノール分解細菌の増殖では、フェノール濃度の増加速度は、下記演算式（７）：

｛式中、Ｓ_Pheは、フェノール濃度であり、ｔは、時間であり、Ｓ_O2は、溶存酸素濃度であり、Ｋ_O2は、溶存酸素濃度に対する飽和定数であり、Ｋ_Pheは、フェノール濃度に対する飽和定数であり、Ｘ_Pheは、フェノール分解細菌濃度であり、そしてμ_Pheは、フェノール分解細菌の最大比増殖速度であり、Ｙ_Pheは、フェノール分解細菌の増殖収率である。｝で表される。

　成分濃度換算工程１２は、ＣＯＤ換算工程４と同様に、図２Ａ～図２Ｃに示すフェノール、チオ硫酸、チオシアン酸濃度とＣＯＤ濃度との相関関係３により、処理水の各成分ＣＯＤ濃度１１を各成分濃度に再換算した。

　上記のシミュレーションの結果を図４に示す。図４は、バッチ試験の各成分濃度の経時変化を示す。以上述べたように、本実施例によれば、生物分解性のある化合物成分を含む廃水の生物学的好気処理プロセスにおいて、活性汚泥モデルを適用することで、処理水の各成分濃度を予測できる。さらに、元々フェノール、チオ硫酸、チオシアン酸の各成分毎に計算を行っているため、単成分であっても２つの成分の組み合わせであっても、処理水の各成分濃度を予測できる。

実施例２：フェノール、チオ硫酸、チオシアン酸を対象としたバッチ試験のシミュレーション（２成分を1種類の細菌が分解する場合）
　以下、実施例１においてチオ硫酸、チオシアン酸の２成分を1種類の細菌（硫黄成分分解細菌）が分解すると設定した時のシミュレーション方法について説明する。ここで、チオ硫酸、チオシアン酸分解細菌を硫黄成分分解細菌としたのは、実際の処理では、硫黄成分分解細菌は硫黄酸化細菌が相当すると推測されたためである。

　シミュレーションは図１と同じ流れで行い、生物分解性のある化合物成分、バッチ試験方法、分析工程２、パラメーター設定工程６、溶存酸素濃度７及び成分濃度換算工程１２は実施例１と同じ方法で行った。但し、ＣＯＤ換算工程４、成分濃度換算工程１２は硫黄換算工程に変更し、表１の成分濃度、化学量論パラメーター、反応速度パラメーターを、以下の表４に示すように設定し、また、計算工程１０ではチオ硫酸分解細菌濃度及びチオシアン酸分解細菌濃度を、硫黄成分分解細菌濃度にまとめ、以下の表４に示す：

　また、以下の表５に示すように、プロセスρ₂、ρ₃における化学量論の関係では、チオ硫酸及びチオシアン酸の分解による細菌増殖が硫黄成分分解細菌に発生することとした。

　また、以下の表６に示すように、プロセスρ₂、ρ₃における反応速度式の細菌濃度の項を硫黄成分分解細菌に置き換えた。

　上記のシミュレーションの結果を図５に示す。図５は、バッチ試験の各成分濃度の経時変化を示す。図４と比べて、チオ硫酸について分解速度が増加した。これは、実施例１ではチオ硫酸分解細菌（Ｘ_S2O3）とチオシアン酸分解細菌（Ｘ_SCN）を分けて各々の分解細菌濃度を算出したが、本実施例ではチオ硫酸とチオシアン酸の２成分を硫黄成分分解細菌１種類で分解すると設定したため、チオ硫酸分解細菌（Ｘ_S2O3）とチオシアン酸分解細菌（Ｘ_SCN）の増殖が硫黄成分分解細菌濃度（Ｘ_S）に集約されたことが主な原因である。また、チオシアン酸については、パラメーター設定行程６でキャリブレーションにより硫黄成分分解細菌濃度（Ｘ_S）を設定した際に、実施例１のＸ_ＳＣＮ初期値より若干低い結果となり、分解速度がわずかに低下したが大きな違いは見られなかった。
ここで、実施例１と本実施例のいずれが良いかは、実験値と計算値との比較や、実際の現象に近い反応が計算できているか等によって判断し、選択することが望ましい。本実施例においては、計算値の精度向上が可能となったと考えられた。

　以上述べたように、本実施例によれば、生物分解性のある化合物成分を分解する細菌の種類および濃度を決定する方法については、基本的に一成分に対し一種類の細菌がその分解に対応することが望ましいが、実際の現象に即した計算予測や計算予測精度の向上のためには、複数成分に1種類の細菌が対応するよう設定してもよい。また、チオ硫酸、チオシアン酸濃度を硫黄濃度に換算することで、硫黄の物質収支を追跡することが可能となり、硫黄成分の挙動が把握できるようになった。また、比内生呼吸速度を考慮した計算を行ったが、結果は無視しうる範囲であり、明確な差異は見られなかった。

実施例３：クレゾールを対象とした酸素消費速度試験のシミュレーション
　以下、生物分解性のある有機物成分としてｏ－クレゾールを対象とした酸素消費速度試験のシミュレーション方法について説明する。
　酸素消費速度試験は、１Ｌの反応容器に、ｏ－クレゾール濃度１５ｍｇ／Ｌとなるよう成分を添加し、ＭＬＳＳ濃度は３３０ｍｇ／Ｌで行った。試験中は、溶存酸素濃度計により溶存酸素濃度を測定し、ｐＨ計によりｐＨを測定した。また、溶存酸素濃度を３．２５ｍｇ／Ｌ、ｐＨ７．５に制御した。

　分析工程２は、生物学的処理プロセスに流入する廃水中の成分濃度を分析する工程である。ここでは上記の酸素消費速度試験装置中の初期成分濃度であるが、ｏ－クレゾール濃度については、ＪＩＳ　Ｋ０１０２吸光光度法により確認を行った。ＣＯＤ換算工程４では、分析工程２で得られたｏ－クレゾールの濃度をＣＯＤ濃度に換算するＣＯＤ換算工程である。ＣＯＤ換算には、下記化学反応式（４）：

からｏ－クレゾールの理論的酸素消費量を算出した。その結果、ｏ－クレゾール１ｍｇあたりの理論酸素消費量は２．５２ｍｇであった。

　この相関関係３により換算を行い、ｏ－クレゾール濃度をＣＯＤ濃度９に換算した。パラメーター設定工程６は、計算工程１０で用いるパラメーター設定を行った。
　パラメーター５は、増殖収率、最大比増殖速度、飽和定数、及び後述する計算工程１０で用いるクレゾール分解細菌濃度を設定した。設定方法は、増殖収率をｏ－クレゾールのＣＯＤ換算濃度と、酸素消費速度試験から得られた酸素消費量から演算式（３）によって求めた。具体的には、ｏ－クレゾールを１３ｍｇ／Ｌとなるよう添加したので、ＣＯＤ濃度は相関関係３から３２．７ｍｇ／Ｌとなる。

　この時の酸素消費量を求める際には、酸素消費速度試験での溶存酸素濃度の減少を常に見る必要がある。しかしながら、試験中には溶存酸素濃度が３．２５ｍｇ／Ｌにするために、ばっ気を行うため、その時の正確な酸素消費速度を求めることはできない。そこで、ここでは、ばっ気時の酸素消費速度を、ばっ気前後における無ばっ気時の酸素消費速度の値から線形近似を行って求めた。たとえば、時間ｔ_１における酸素消費速度を１０ｍｇ／Ｌ／ｈｒ、時間ｔ₂における酸素消費速度を７．５ｍｇ／Ｌ／ｈｒとし、時間ｔ₁から時間ｔ₂の間、ばっ気がされていた時、その間の時間Ｔにおける酸素消費速度を、下記演算式（８）：

｛式中、ＯＵＲ_Tは、時間Ｔにおける酸素消費速度であり、ｔ₁は、ばっ気前の時間であり、ｔ₂は、ばっ気後の時間であり、そしてＴは、ばっ気期間において、酸素消費速度を求めたい時間である。｝を用いて求めた。但し、これに限らず、その他の近似方法を使ってもよい。

　このようにして求めた酸素消費速度を時間で積分することで、酸素消費量は１３．９ｍｇ／Ｌであったことから、増殖収率は（３２．７－１３．９）／３２．７＝０．５７と算出された。また、分解細菌はクレゾール分解細菌が対応すると考えた。最大比増殖速度、飽和定数、及びクレゾール分解細菌濃度は、キャリブレーションによる方法で求め、実測値、計算値及び実測値の標準偏差から、演算式（２）に示されるχ²の値を求め、χ²の値が最小となったときのパラメーターの値を算出した。酸素消費速度及び成分濃度の実測値は、本実施例での酸素消費速度試験中に実測して求めた。計算値は、後述する計算工程１０の方法により求めた。

　溶存酸素濃度７は、バッチ試験の制御値３．２５ｍｇ／Ｌを用いた。
　計算工程１０は、ＣＯＤ濃度９、溶存酸素濃度７及びパラメーター５を用いて、ｏ－クレゾール濃度の経時変化についてシミュレーションを行い、処理水のｏ－クレゾールのＣＯＤ濃度１１を求めた。計算方法は、演算式（１）を用いた。以下の表７は、計算に用いる各成分濃度及びパラメーターを示す：

　以下の表８は、計算に用いる化学量論を示す：

　そして以下の表９は、計算に用いた反応速度式を示す：

　成分濃度換算工程１２は、ＣＯＤ換算工程４と同様に、ｏ－クレゾール濃度とＣＯＤ濃度との相関関係３により、処理水のｏ－クレゾールＣＯＤ濃度１１をｏ－クレゾール濃度１３に再換算した。
　上記のシミュレーションの結果を図６に示す。図６は、ｏ－クレゾール濃度の経時変化を示す。以上述べたように、本実施例によれば、フェノール、チオ硫酸、チオシアン酸以外の成分にも本発明が適用可能であり、相関関係３に理論的酸素消費量を用いても成分濃度を求めることができる。

　本発明によれば、活性汚泥モデルでは対象外であった溶解性遅分解性成分についても、その成分を分解する細菌の種類を設定することで、処理水の成分濃度を計算することが可能であり、産業上有用である。

　１　　水質シミュレーション方法の概要図
　２　　分析工程
　３　　各成分とＣＯＤ濃度との相関関係
　４　　ＣＯＤ換算工程
　５　　パラメーター
　６　　パラメーター設定工程
　７　　溶存酸素濃度
　８　　溶存酸素濃度測定工程
　９　　ＣＯＤ濃度
　１０　　計算工程
　１１　　処理水の各成分ＣＯＤ濃度
　１２　　成分濃度換算工程
　１３　　処理水の各成分濃度
　２１　　酸素消費速度試験装置
　２２　　微生物汚泥
　２３　　硝化阻害剤
　２４　　撹拌装置
　２５　　栄養塩
　２６　　対象成分
　２７　　溶存酸素濃度計
　２８　　データ記録装置
　２９　　空気供給装置
　３０　　ｐＨ計
　３１　　酸・アルカリ供給装置
　３２　　ヒーター　３３　　恒温水槽

Claims

　生物分解性を有する化合物成分を含む廃水を生物反応槽内で生物学的好気処理するプロセスにおける水質シミュレーション方法であって、
　前記生物反応槽に流入する前記廃水中の前記化合物成分の各成分濃度を分析する分析工程と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、各成分濃度の分析値を各成分のＣＯＤ濃度に換算するＣＯＤ換算工程と、
　化学量論パラメーターである増殖収率、反応速度式パラメーターである飽和定数、最大比増殖速度、並びに前記化合物成分を分解する細菌の種類および濃度を設定するパラメーター設定工程と、
　前記生物反応槽の溶存酸素濃度を測定する溶存酸素濃度測定工程と、
　前記各成分のＣＯＤ濃度、前記増殖収率、前記飽和定数、前記最大比増殖速度、並びに前記化合物成分を分解する細菌の種類及び濃度、並びに溶存酸素濃度を用いて、下記演算式（１）：

｛式中、Ｃ_iは、各成分濃度であり、ｉは、各成分を表す通し番号であり、Ｐ_ijは、化学量論パラメーターであり、ｊは、各プロセスを表す通し番号であり、そしてρ_jは、反応速度式（反応速度式パラメーターを含む速度式）である。｝により、前記生物学的反応槽内で前記生物学的好気処理した後の処理水における各成分のＣＯＤ濃度を算出する計算工程と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、当該算出された生物学的処理水の各成分のＣＯＤ濃度を、各成分濃度に換算する成分濃度換算工程とを、含む水質シミュレーション方法。
　前記パラメーター設定工程は、
（ア）前記流入する廃水及び処理水において、事前に時系列的に別途採取したそれぞれの生物分解性を有する化合物成分濃度と溶存酸素濃度を用いて、キャリブレーションによって決定する方法、
（イ）前記各成分それぞれに対して、各成分単独で前記生物学的好気処理を行って連続的に溶存酸素濃度を測定し、当該測定値から酸素消費速度を算出し、当該算出された酸素消費速度のデータから決定する方法、のいずれかの方法を用いて、前記化学量論パラメーター及び前記反応速度式パラメーターを設定する、請求項１に記載の水質シミュレーション方法。
　上記廃水が、コークス製造工程で発生する安水であり、上記化合物成分のうち、有機物成分はフェノール、並びに、無機物成分はチオ硫酸及びチオシアン酸であり、
　チオ硫酸、チオシアン酸のシミュレーションについては、上記ＣＯＤ換算工程の代わりにチオ硫酸及びチオシアン酸濃度を硫黄換算する硫黄換算工程と、
　上記生物学的反応槽で上記生物学的好気処理した後の処理水における各成分の硫黄濃度を算出する硫黄計算工程と、
　当該算出された生物学的処理水の各成分の硫黄濃度を、チオ硫酸濃度及びチオシアン酸濃度に再換算する硫黄再換算工程と
　を有する、請求項１又は２に記載の水質シミュレーション方法。
　生物分解性を有する化合物成分を含む廃水を生物反応槽内で生物学的好気処理するプロセスで用いる水質シミュレーション装置であって、
　前記生物反応槽に流入する前記廃水中の前記化合物成分の各成分濃度を分析する分析手段と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、各成分濃度の分析値を各成分のＣＯＤ濃度に換算するＣＯＤ換算手段と、
　化学量論パラメーターである増殖収率、反応速度式パラメーターである飽和定数、最大比増殖速度、並びに前記化合物成分を分解する細菌の種類および濃度を設定するパラメーター設定手段と、
　前記生物反応槽の溶存酸素濃度を測定する溶存酸素濃度測定手段と、
　前記各成分のＣＯＤ濃度、前記増殖収率、前記飽和定数、前記最大比増殖速度、並びに前記化合物成分を分解する細菌の種類及び濃度、並びに溶存酸素濃度を用いて、下記演算式（１）：

｛式中、Ｃ_iは、各成分濃度であり、ｉは、各成分を表す通し番号であり、Ｐ_ijは、化学量論パラメーターであり、ｊは、各プロセスを表す通し番号であり、そしてρ_jは、反応速度式（反応速度式パラメーターを含む速度式）である。｝により、前記生物学的反応槽で前記生物学的好気処理した後の処理水における各成分のＣＯＤ濃度を算出する計算手段と、
　前記各成分濃度と、ＣＯＤ_Cr、ＣＯＤ_Mn及び理論的酸素消費量から選ばれる一つのＣＯＤ濃度との相関関係をもとに、当該算出された生物学的処理水の各成分のＣＯＤ濃度を、各成分濃度に換算する成分濃度換算手段と
を、含む水質シミュレーション装置。
　前記パラメーター設定手段は、
（ア）前記流入する廃水及び処理水において、事前に時系列的に別途採取したそれぞれの生物分解性を有する化合物成分濃度と溶存酸素濃度を用いて、キャリブレーションによって決定する手段、
（イ）前記各成分それぞれに対して、各成分単独で前記生物学的好気処理を行って連続的に溶存酸素濃度を測定し、当該測定値から溶存酸素濃度を算出し、当該算出された酸素消費速度のデータから決定する手段、のいずれかの手段を用いて、前記化学量論パラメーター及び前記反応速度式パラメーターを設定する、請求項４に記載の水質シミュレーション装置。
　上記廃水が、コークス製造工程で発生する安水であり、上記化合物成分のうち、有機物成分はフェノール、並びに、無機物成分はチオ硫酸及びチオシアン酸であり、
　チオ硫酸、チオシアン酸のシミュレーションについては、上記ＣＯＤ換算手段の代わりにチオ硫酸及びチオシアン酸濃度を硫黄換算する硫黄換算手段と、
　上記生物学的反応槽で上記生物学的好気処理した後の処理水における各成分の硫黄濃度を算出する硫黄計算手段と、
　当該算出された生物学的処理水の各成分の硫黄濃度を、チオ硫酸濃度及びチオシアン酸濃度に再換算する硫黄再換算手段と
　を有する、請求項４又は５に記載の水質シミュレーション装置。