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