(57)【要約】【課題】総括酸素移動容量係数の測定時間が不要で、季節の変動にも短時間で関係式が得られる予測精度の高い運転支援装置を得る。【解決手段】本発明の運転支援装置は、下水処理場の水質データを蓄積するデータ蓄積部（１）と、施設条件設定部（２）と、運転条件を設定する運転条件設定部（３）と、風量を設定する風量設定部（６）と、水質を予測する予測処理部（７）と、最適な運転をするための操作変数を選択する最適運転選択部（８）と、最適運転選択部からの出力値を実プラントの運転設備に送る操作変数の設定値を設定する操作変数設定部（９）とを備え、データ蓄積部（１）と予測処理部（７）との間に、関連する変数間の時間差を調整する時間差調整部（４）と、溶存酸素濃度に影響を与える変数から主成分モデルを作成し、そのモデルから溶存酸素濃度を演算する主成分処理部（５）とを備えた構成である。【選択図】　　図１
2004066120.tif 000002
【特許請求の範囲】
【請求項１】
最初沈殿池、反応槽、最終沈殿池などの処理工程をもつ下水処理場の水質データを蓄積するデータ蓄積部（１）と、土木構造などの処理場仕様を設定する施設条件設定部（２）と、運転条件を設定する運転条件設定部（３）と、前記反応槽に注入する空気の風量を設定する風量設定部（６）と、前記データ蓄積部（１）のデータを基に水質を予測する予測処理部（７）と、最適な運転をするための操作変数を選択する最適運転選択部（８）と、前記最適運転選択部からの出力値を実プラントの運転設備に送る操作変数の設定値を設定する操作変数設定部（９）とを備えた運転支援装置において、
前記データ蓄積部（１）と予測処理部（７）との間に関連する変数間の時間差を調整する時間差調整部（４）と、溶存酸素濃度に影響を与える変数から主成分モデルを作成し、そのモデルから溶存酸素濃度を演算する主成分処理部（５）とを備えたことを特徴とする運転支援装置。。
【請求項２】
前記時間差調整部（４）は、溶存酸素濃度とそれに影響を与える変数間の時間差を演算する時間差演算部（４１）と、前記時間差演算部で演算した時間差から変数間の時間差を修正する時間差修正部（４２）とを備えたことを特徴とする請求項１記載の運転支援装置。
【請求項３】
前記時間差演算部（４１）は、溶存酸素濃度とそれに影響を与える変数の相関関係を演算する相互相関演算部（４１１）と、変数に一定の入力を与えた時の応答を演算するステップ応答演算部（４１２）の少なくとも１つを備えたことを特徴とする請求項２記載の運転支援装置。
【請求項４】
前記主成分処理部（５）は、溶存酸素濃度に影響を与える変数の主成分モデルを作成する主成分モデル作成部（５１）と、変数の主成分得点を演算する主成分得点演算部（５２）と、前記主成分得点演算部で作成した主成分得点から溶存酸素濃度へ変換する式を作成するデータ変換式作成部（５３）と、前記主成分モデル作成部で作成されたモデルと前記データ変換式作成部で作成された変換式と前記風量設定部で設定した値と前記データ蓄積部で蓄積された値とから溶存酸素濃度を演算する主成分溶存酸素演算部（５４）とを備えたことを特徴とする請求項１から３のいずれか１項に記載の運転支援装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は都市下水や産業廃水のように、汚水を浄化する活性汚泥法、嫌気・好気汚泥処理法による運転支援装置に関する。
【０００２】
【従来の技術】
都市下水処理、工場廃水などの処理プロセスは、微生物の動作状態、気候などから運転方法を決定する必要がある。図２は従来の運転支援装置とその周辺との接続状況を示すブロック図である。図において、１はデータ蓄積部、２は施設条件設定部、３は運転条件設定部、６は風量設定部、７は予測処理部、８は最適運転選択部、９は操作変数設定部、１０は総括酸素移動容量係数測定部、１１は溶存酸素演算部である。
先ず、下水処理場の処理概要を説明する。
下水処理場では、流入した下水は、最初沈殿池１００、反応槽１０１、最終沈殿池１０２の順に流れて処理される。汚濁物質を含む下水は、最初沈澱池１００に導入され、汚濁物質の中の沈降しやすいものを沈降分離して上澄水を反応槽１０１に流出する。反応槽１０１には最終沈澱池１０２の汚泥の一部が返送汚泥ポンプによって返送されており、反応槽１０１はその返送汚泥と最初沈殿池１００の上澄水を処理する。反応槽１０１では、ブロワー（図示せず）から圧送された空気が曝気槽内の散気管によって放出されており、汚濁物質は活性汚泥により吸着、分解されて最終沈澱池１０２に導かれる。最終沈澱池１０２では活性汚泥を沈降分離し、沈降汚泥は余剰汚泥ポンプ（図示せず）により汚泥処理系（図示せず）に排出され、清澄水は処理水として滅菌槽（図示せず）を経て放流される。
次にこの運転支援システムについて説明する。
データ蓄積部１は、下水処理場の水質予測に必要な情報を収集・蓄積し、総括酸素移動容量係数設定部１０へデータを出力する。ここで収集されるデータはリン酸態リン濃度、硝酸態窒素濃度、アンモニア態窒素濃度、Ｔｏｔａｌ−ＣＯＤ、溶解性ＣＯＤなど水質情報と流量、風量、溶存酸素濃度などの運転情報である。
施設条件設定部２は、水質の予測に必要な反応槽の体積、管路などの処理場土木構造を設定し、予測処理部７へ出力する。
運転条件設定部３は、プラントを動作させるために必要な運転指標である返送率、ＳＲＴなどをオペレータの判断で設定する。
風量設定部６は、現在の風量を継続した時の水質状態を把握するために現在の風量を入力、または種々の運転方法を検討するため、最適水質になると思われる風量を入力する。
総括酸素移動容量係数設定部１０は、溶存酸素濃度と酸素利用速度の測定値により総括酸素移動容量係数を設定する。好気タンクの所定位置で、時間を変えて数回混合液を採取し、溶存酸素濃度と酸素利用速度を測定する。好気タンク内において、酸素利用速度と総括酸素移動容量係数には、次の（１）式が成り立つ。
【０００３】
【数１】
2004066120.tif 000003
【０００４】
Ｃｓ：好気タンク内混合液の飽和溶存酸素濃度（ｍｇ／ｌ）
Ｃ：好気タンク内混合液の溶存酸素濃度（ｍｇ／ｌ）
ｒｒ：好気タンク内混合液の酸素利用速度（ｍｇ／ｌ・ｈ）
ＫＬａ：総括酸素移動容量係数（１／ｈ）
本測定では好気タンク内は定常状態のため、ｄｃ／ｄｔ＝０である。つまり、グラフの縦軸に溶存酸素濃度、横軸に酸素利用速度をプロットすると直線が得られ、その傾きの逆数が総括酸素移動容量係数となる。
溶存酸素演算部１１は、総括酸素移動容量係数と風量から溶存酸素濃度を演算する。演算は（２）式にて行う。
【０００５】
【数２】
2004066120.tif 000004
【０００６】
Ｑ：必要送風量［Ｎｍ３／ｍｉｎ］
ηＤ：酸素吸収効率０．１［−］
ρ：空気の密度１．２９［ｋｇ／Ｎｍ３］
α：清水に対するＫＬａ係数０．８５［−］
β：清水に対する酸素飽和濃度係数０．９５［−］
ＫＬａ：総括酸素移動容量係数［１／ｈ］
Ｖ：対象槽の反応槽容積［ｍ３］
Ｃｓ：酸素飽和濃度２０℃の時［ｍｇ／ｌ］
Ｔ：水温［℃］
予測処理部７は、データ蓄積部１に蓄積された分析データのうち、下水処理場の入口に流入する水質データと、施設条件設定部２で設定された施設の値と、運転条件設定部３で設定された条件と、溶存酸素演算部１１で演算した溶存酸素濃度を入力し、反応槽内、処理水の窒素・リンなどの水質予測を行う。その予測計算に用いる式は１０数個あり、その１つを例示すると、例えば発酵生成物濃度の計算はつぎの（３）式を用いて行われる。
【０００７】
【数３】
2004066120.tif 000005
【０００８】
記号の意味はつぎのとおりである。
ＳＡ　（ｉ）：ｉ時点の対象タンクの発酵生成物濃度（ｇＣＯＤ／ｍ３）
ＳＡ（ｉ）Ｒ：ｉ時点の対象タンクの化学反応による変化量を考慮に入れた発酵生成物
濃度（ｇＮ／ｍ３）
ＳＡ（ｉ）ｉｎ：ｉ時点の対象タンクに流入する発酵生成物濃度（ｇＣＯＤ／ｍ３）
ＳＡ（ｉ）ｏｕｔ：ｉ時点の対象タンクから流出する発酵生成物濃度（ｇＣＯＤ／ｍ３）
Ｖ：対象タンクの体積（ｍ３）
Ｑｉｎ：対象タンクへ流入する量（ｍ３／ｈ）
Ｑｏｕｔ：対象タンクから流出する量（ｍ３／ｈ）
Ｔ：対象タンクから流入・流出する時間（ｈ）
ｉ：シミュレーション予測間隔（ｈ）
但し各変数は次の（４）式〜（１１）式のように計算される。
【０００９】
【数４】
2004066120.tif 000006
【００１０】
記号の意味はつぎのとおりである。
ＳＡ（ｉ）：ｉ時点の対象タンクで増殖・減少した発酵生成物（ｇＣＯＤ／ｍ３）
ρ５：ＳＡによる好気的増殖（ｇＣＯＤ／ｄ）
ρ７：ＳＡによる無酸素的増殖と脱窒（ｇＣＯＤ／ｄ）
ρ８：発酵（ｇＣＯＤ／ｄ）
ρ１０：ポリヒドロキシアルカノエートの蓄積（ｇＣＯＤ／ｄ）
ρ１５：ポリヒドロキシアルカノエートの分解（ｇＣＯＤ／ｄ）
データ蓄積部１から入力される分析データは以下の通りである。
ＳＯ２：溶存酸素濃度（ｇＯ２／ｍ３）
ＳＮＨ４：溶解性のアンモニア濃度（ｇＮ／ｍ３）
ＳＮＯ３：溶解性の硝酸性窒素濃度（ｇＮ／ｍ３）
ＳＰＯ４：無機溶解性リン酸性リン濃度（ｇＰ／ｍ３）
ＳＡＬＫ：アルカリ度（ｍｏｌｅ　ＨＣＯ３／ｍ３）
ＳＦ：易分解性有機物濃度（ｇＣＯＤ／ｍ３）
ＳＡ：発酵生成物濃度（ｇＣＯＤ／ｍ３）
ＸＡＵＴ：硝化菌濃度（ｇＣＯＤ／ｍ３）
ＸＨ：従属栄養微生物（ｇＣＯＤ／ｍ３）
ＸＰＨＡ：ポリヒドロキシアルカノエート（ｇＣＯＤ／ｍ３）
ＸＰＰ：ポリリン酸（ｇＰ／ｍ３）
ＸＰＡＯ：リン酸蓄積細菌（ｇＣＯＤ／ｍ３）
モデルパラメータはつぎのとおりである。
μＨ：基質の最大増殖速度（１／ｄ）
ηＮＯ３：硝酸性窒素濃度の無酸素状態の加水分解による減少係数（−）
ＫＡ：発酵生成物濃度飽和係数（ｇＯ２／ｍ３）
ＫＯ２：酸素飽和係数（ｇＯ２／ｍ３）
ＫＦ：易分解性基質増殖飽和係数（ｇＣＯＤ／ｍ３）
ＫＮＯ３：硝酸性窒素飽和係数（ｇＮ／ｍ３）
ただし、このシミュレーションでは亜硝酸性窒素と硝酸性窒素を同じ物質として扱う。
ＫＮＨ４：アンモニア飽和係数（ｇＮ／ｍ３）
ＫＰ：リン酸性リン飽和係数（ｇＰ／ｍ３）
ＫＡＬＫ：アルカリ度飽和係数（ｍｏｌｅ　ＨＣＯ３／ｍ３）
Ｋｆｅ：ＳＦの発酵飽和濃度（ｇＣＯＤ／ｍ３）
ｑｆｅ：発酵の最大速度（１／ｄ）
ＫＰＰ：ポリリン酸の飽和濃度（ｇＰＰ／ｇＰＡＯ）
ｑＰＨＡ：ポリヒドロキシアルカノエート蓄積速度（ｇＣＯＤ／（ｇＰＡＯ・ｄ）
ｂＰＨＡ：ポリヒドロキシアルカノエート分解速度（１／ｄ）
予測処理部７では、上記のように硝酸性窒素だけでなく、リン酸性リン濃度、アルカリ度、アンモニア性窒素濃度なども同時に計算される。その予測結果を最適運転選択部８に出力する。
最適運転選択部８では、予測処理部７の結果がプラント運転に適しているかの判断をコスト、水質などから検討を行う。その結果がオペレータの満足したものであれば、求められた操作変数を操作変数設定部９に出力する。予測処理部７の結果が水質の基準値を超えたり、コストがかかるなど不満足なものであれば、施設条件設定部２、運転条件設定部３、風量設定部６の値を再設定し、再度予測処理部７にて水質予測を行う。
操作変数設定部９では、最適運転選択部８で求められた最適な操作変数（送風量、ポンプ台数など）を実プラントへ設定する。
【００１１】
【発明が解決しようとする課題】
しかしながら、従来の運転支援装置では、反応槽に送風する空気の風量から溶存酸素濃度の算出に必要な総括酸素移動容量係数を測定することに多くの時間を要していた。そのため、一度測定した総括酸素移動容量係数を長期間使用することが多く、季節変動などへの対応が困難な状況となり、水質の予測精度に問題があった。
そこで、本発明は総括酸素移動容量係数の測定に時間を要することがなく、季節の変動にも短時間で関係式が得られる予測精度の高い運転支援装置を提供することを目的とする。
【００１２】
【課題を解決するための手段】
上記問題を解決するため、本発明は、次のように構成している。
請求項１記載の運転支援装置は、最初沈殿池、反応槽、最終沈殿池などの処理工程をもつ下水処理場の水質データを蓄積するデータ蓄積部（１）と、土木構造などの処理場仕様を設定する施設条件設定部（２）と、運転条件を設定する運転条件設定部（３）と、前記反応槽に注入する空気の風量を設定する風量設定部（６）と、前記データ蓄積部（１）のデータを基に水質を予測する予測処理部（７）と、最適な運転をするための操作変数を選択する最適運転選択部（８）と、前記最適運転選択部からの出力値を実プラントの運転設備に送る操作変数の設定値を設定する操作変数設定部（９）とを備えた運転支援装置において、前記データ蓄積部（１）と予測処理部（７）との間に関連する変数間の時間差を調整する時間差調整部（４）と、溶存酸素濃度に影響を与える変数から主成分モデルを作成し、そのモデルから溶存酸素濃度を演算する主成分処理部（５）とを備えたものである。
本構成によれば、風量から溶存酸素濃度を自動的に演算することができるようになり、季節変動に対応した主成分モデルおよび変換式が作成できるため、予測精度を向上させることができる。
請求項２記載の運転支援装置は、前記時間差調整部（４）は、溶存酸素濃度とそれに影響を与える変数間の時間差を演算する時間差演算部（４１）と、前記時間差演算部で演算した時間差から変数間の時間差を修正する時間差修正部（４２）とを備えたものである。
本構成によれば、溶存酸素とそれ以外の変数間の時間遅れを自動的に調整することができ、溶存酸素濃度算出時の精度が向上する。
請求項３記載の運転支援装置は、前記時間差演算部（４１）は、溶存酸素濃度とそれに影響を与える変数の相関関係を演算する相互相関演算部（４１１）と、変数に一定の入力を与えた時の応答を演算するステップ応答演算部（４１２）の少なくとも１つを備えたものである。
本構成によれば、溶存酸素濃度を自動的に風量、活性汚泥浮遊物質、生物化学的酸素要求量などから演算することができるため、季節変動への対応も可能となり、溶存酸素濃度を計算する精度が向上し、溶存酸素濃度の計算時間の短縮も可能となる。
請求項４記載の運転支援装置は、前記主成分処理部（５）は、溶存酸素濃度に影響を与える変数の主成分モデルを作成する主成分モデル作成部（５１）と、変数の主成分得点を演算する主成分得点演算部（５２）と、前記主成分得点演算部で作成した主成分得点から溶存酸素濃度へ変換する式を作成するデータ変換式作成部（５３）と、前記主成分モデル作成部で作成されたモデルと前記データ変換式作成部で作成された変換式と前記風量設定部で設定した値と前記データ蓄積部で蓄積された値とから溶存酸素濃度を演算する主成分溶存酸素演算部（５４）とを備えたものである。
本構成によれば、溶存酸素濃度とそれに影響を与える変数間の時間差を求めることができ、変数間の関係式を構築する時に時間調整ができるため、風量、活性汚泥浮遊物質、生物化学的酸素要求量などから溶存酸素濃度を計算する精度が向上する。
【００１３】
【発明の実施の形態】
以下、本発明の実施の形態を図に基づいて説明する。
図１は、本発明の運転支援装置の構成とその周辺との接続状況を示すブロック図である。図において、４は時間差演算部４１と時間差修正部４２とからなる時間差調整部、５は主成分処理部である。
時間差演算部４１は相互相関演算部４１１とステップ応答演算部４１２とからなり、主成分処理部５は主成分モデル作成部５１、主成分得点変換部５２、データ変換式作成部５３、主成分溶存酸素演算部５４からなる。なお、同じ部品には同じ符号を付しており、従来技術の項で説明したものは説明を省略する。
つぎに動作について述べる。
▲１▼　データ蓄積部１は、下水プロセスの反応槽から処理水までの水質予測に必要な情報を収集・蓄積し、時間差演算部４へデータを出力する。収集されるデータは従来と同じである。
▲２▼　施設条件設定部２および運転条件設定部３は、従来と同じ機能である。
▲３▼　時間差調整部４は、データ処理を行う前に変数間の時間差を調整する。つまり、風量が増加して溶存酸素濃度が増加するまでは、数十分の時間がかかるため、データ処理においてデータを同一に取り扱うには時間差を算出し、時間差の影響を除去することが必要となる。
時間差演算部４１は、相互相関演算部４１とステップ応答演算部４１２からなる。相互相関演算部４１１は、溶存酸素濃度と生物化学的酸素要求量、活性汚泥浮遊物質、送風量の相互相関を算出する演算部である。つまり、風量が増加して溶存酸素濃度が増加するまで時間がかかるため、データ処理においてデータを同一に取り扱うには時間差の算出が必要となる。時間差は相互相関関数を使用して求める。長さＮの多変量時系列［ｙ１（ｊ），・・・ｙＮ（ｊ）］（ｊ＝１，・・・，ｌ）が観測された時μ（ｉ）、相互共分散関数Ｃｋ（ｉ，ｊ）および相互相関関数Ｒｋ（ｉ，ｊ）の推定値は（１２）式〜（１４）式を用いて求められる。
先ず、各時系列データの平均値を算出する。
【００１４】
【数５】
2004066120.tif 000007
【００１５】
つぎに、相互共分散関数を（１３）式にて算出する。
【００１６】
【数６】
2004066120.tif 000008
【００１７】
相互相関関数は（１４）式にて算出する。
【００１８】
【数７】
2004066120.tif 000009
【００１９】
相互共分散関数、相互相関関数と呼ばれる。
溶存酸素濃度と送風量、活性汚泥浮遊物質、生物化学的酸素要求量などの相互相関関数は上記の式で表現することができる。溶存酸素濃度と送風量の相互相関関数の値が最も高い点が、変数間に相関があるとみなせるため、溶存酸素濃度と送風量の時間差の調整を行う場合の時間となる。
ステップ応答演算部４１２は、風量、活性汚泥浮遊物質、生物化学的酸素要求量などと溶存酸素濃度のステップ応答を演算する。ステップ応答は風量をステップ状に変化させ、その時の溶存酸素濃度の変化を見る。風量を変化させて溶存酸素濃度が変化し、しばらくすると一定になる。一定になった時の時間をステップ応答時間とする。
時間差修正部４２は、相互相関演算部４１１またはステップ応答演算部４１２から算出した時間差から、対象データ間の時間差の影響を削除するようにデータ修正を行う。つまり、溶存酸素濃度と送風量のデータを時間差から（１５）式のように調整する。
【００２０】
【数８】
2004066120.tif 000010
【００２１】
ＤＯ（ｉ）：ｉ時刻の溶存酸素濃度
ＡＱ（ｉ）：ｉ時刻の風量
ｔ：溶存酸素濃度と送風量の相互相関係数が最も高い時間差
▲４▼　主成分処理部５は、風量、活性汚泥浮遊物質、生物化学的酸素要求量などからそれらの代表値である主成分得点を演算し、主成分得点から溶存酸素濃度への変換式を作成して溶存酸素濃度を計算する。主成分処理部５は、主成分モデル作成部５１と主成分得点作成部５２とデータ変換式作成部５３と主成分溶存酸素演算部５４からなる。主成分モデル作成部５１は、溶存酸素濃度へ影響を与える風量、活性汚泥浮遊物質、生物化学的酸素要求量などを対象データとして主成分分析を行う。主成分分析は、ｐ種類の変量をＸ１、Ｘ２、・・・Ｘｐ、任意の係数ａ１、ａ２、・・・ａｐを用いて次のような１次式による総合的指標（合成変量）をつくる。
【００２２】
【数９】
2004066120.tif 000011
【００２３】
ここで、ａ＝（ａ１、・・・・、ａｐ）’、Ｘ＝（Ｘ１、・・・・、Ｘｐ）’である。
このように結合されたｚは、ｐ個の変量Ｘをよく代表していなければならない。そのための基準として主成分ｚの分散の最大化を行う。
以下にこの基準（主成分の分散の最大化）を用いた時の主成分分析の定式化について説明する。
ｐ変量から（１６）式により合成するときｚの分散は、（１５）式となる。
【００２４】
【数１０】
2004066120.tif 000012
【００２５】
ここにｓｊｋはｘｊｋとｘｋとの共分散、Ｓはｓｊｋを（ｊ，ｋ）要素にもつ分散共分散行列である。（１６）式のｚは、ｐ次元空間の中で原点ＯからあるＯＺ方向にｚ軸をとることを意味するが、そのときｚ座標のスケールをｘ１、ｘ２・・・ｘｐ軸を同じにとることにすれば、係数ａ１、ａ２・・・ａｐはそれぞれ直線ＯＺの方向余弦になり、（１８）式を満たす。
【００２６】
【数１１】
2004066120.tif 000013
【００２７】
したがって、問題は（１８）式の制約条件のもとで（１７）式を最大化することになる。このような制約条件つき最大化問題は、ラグランジェ乗数λを用いてつぎの（１９）式を最大化する問題に変形される。
【００２８】
【数１２】
2004066120.tif 000014
【００２９】
（１９）式をａの各要素で偏微分してゼロとおけば、（２０）式のような行列Ｓの固有値問題を得る。
【００３０】
【数１３】
2004066120.tif 000015
【００３１】
Ｓは対称行列であり、任意のａに対してａ’Ｓａ＝Ｖ（ａ’ｘ）≧０より非負値であるから、ｐ個の実数で非負の固有値λ１≧λ２≧・・・・・≧λｐ≧０をもつ。
各固有値ベクトルをａ１、ａ２、ａ３、・・・、ａｐとすれば（２０）式より
【００３２】
【数１４】
2004066120.tif 000016
【００３３】
となり、λｊはｚｊ＝ａｊ’ｘの分散に等しくなる。
したがって分散を最大にする主成分は、最大固有値λ１に対応する固有ベクトルａ１＝（ａ１１、ａ２１、・・・、ａｐ１）’の要素の係数として（２２）式となる。
【００３４】
【数１５】
2004066120.tif 000017
【００３５】
（２２）式を第１主成分と呼ぶ。第１主成分ｚ１の分散はλ１で、ｐ個の変量をよく代表する主成分になっている。
主成分得点演算部５２は、前記主成分モデル作成部で求められた係数ａ１、ａ２・・・・ａｐと現在時刻ｎの風量、活性汚泥浮遊物質、生物化学的酸素要求量などから主成分得点ｚ（ｎ）を算出する。
【００３６】
【数１６】
2004066120.tif 000018
【００３７】
データ変換式作成部５３は、前記主成分得点演算部にて求めた主成分得点と溶存酸素濃度の関係式を作成する。つまり、主成分得点と溶存酸素濃度の回帰式を作成することである。回帰式は１次回帰、２次回帰、指数回帰など処理場に適した方法を選択する。一例として、一次回帰について説明する。主成分得点と溶存酸素濃度の一次回帰は（２４）式のように表すことができる。
【００３８】
【数１７】
2004066120.tif 000019
【００３９】
Ｚ：主成分得点
ＤＯ：溶存酸素濃度
一次回帰は、ＤＯ軸に平行な直線に沿っての距離ｄｉの平方の和（２５）式を最小にするようなａ，ｂの値を求めることになる。
【００４０】
【数１８】
2004066120.tif 000020
【００４１】
主成分溶存酸素演算部５４は、前記主成分モデル作成部５１で作成したモデルに対し、データ蓄積部に蓄積した活性汚泥浮遊物質、生物化学的酸素要求量などと、風量設定部で設定した風量を設定することで、それらの主成分得点ｚを（２６）式のように求める。
【００４２】
【数１９】
2004066120.tif 000021
【００４３】
Ｚ：主成分得点
ＡＱ：風量
そして、前記データ変換式作成部５２で演算した回帰モデルを適用して、主成分得点から溶存酸素濃度を（２４）式のように演算する。
▲５▼　　風量設定部６は、従来と同じ機能である。
▲６▼　予測処理部７は、データ蓄積部１に蓄積されたデータと、施設条件設定部２で設定された土木構造などのプラント仕様データと、運転条件設定部３で設定されたプラントの操作変数と、主成分処理部５で計算された溶存酸素濃度を入力し、水質モデルをベースとした処理場内の水質（窒素、リンなど）の予測を行い、その予測結果を最適運転選択部８に出力する。詳細は従来技術と同様のため省略する。
▲７▼　最適運転選択部８および操作変数設定部９は、従来と同じ機能である。
【００４４】
【発明の効果】
以上述べたように、本発明によれば、風量、活性汚泥浮遊物質、生物化学的酸素要求量などから溶存酸素濃度を計算する計算式に主成分モデルと変換式を使用することによって、溶存酸素濃度と風量、活性汚泥浮遊物質、生物化学的酸素要求量などの関係式を導き出すことができ、総括酸素移動容量係数を求める必要がなくなり、溶存酸素濃度の計算時間が短縮される。さらに、蓄積データを使用するため季節変動にも対応した溶存酸素濃度と風量、活性汚泥浮遊物質、生物化学的酸素要求量などとの関係式が構築可能となることにより、予測精度が向上するという効果がある。
【図面の簡単な説明】
【図１】本発明の運転支援装置を示すブロック図である。
【図２】従来の運転支援装置を示すブロック図である。
【符号の説明】
１　データ蓄積部
２　施設条件設定部
３　運転条件設定部
４　時間差調整部
４１　時間差演算部
４１１　　相互相関演算部
４１２　　ステップ応答演算部
４２　時間差修正部
５　主成分処理部
５１　主成分モデル作成部
５２　主成分得点変換部
５３　データ変換式作成部
５４　主成分溶存酸素演算部
６　風量設定部
７　予測処理部
８　最適運転選択部
９　操作変数設定部
１０　総括酸素移動容量係数
１１　溶存酸素演算部
【図１】
2004066120.tif 000022
【図２】
2004066120.tif 000023