(57)【特許請求の範囲】
【請求項１】
水質データ及び運転状態記憶装置と、運転条件を設定する運転条件設定装置と、プラントの化学反応に関する定数を設定するパラメータ設定装置と、処理場の土木構造やプロセスフローを設定するプラント土木条件設定装置と、未来の流入水状態を設定する未来データ設定装置と、前記各設定装置の設定値に基づき、プロセス各部の水質予測演算を行なうシミュレーション演算装置とを備えた水質シミュレーション装置において、
前記未来データ設定装置は、過去の運転データに基づき変動する流入水状態を予測演算する流入水変動係数演算装置を有し、この演算結果と現時点のデータとに基づき未来の流入水量と水質データとを設定するものとし、
さらに、前記未来データ設定装置の設定値に基づき、流入水の有機成分としては、全TOC濃度と溶解性TOC濃度、または全CODcr濃度と溶解性CODcr濃度、窒素成分としては、全窒素濃度と溶解性窒素濃度と硝酸態窒素濃度とアンモニア態窒素濃度、りん成分としては、全りん濃度と溶解性りん濃度とりん酸態りん濃度、および浮遊物質濃度の測定結果と一致するようにシミュレーションで定義している有機物への分画比やその分画された各有機物に含まれる物質の含有比を演算する流入水質演算装置と、前記分画された各有機物に含まれる物質の含有比を前記シミュレーション演算装置に設定する流入水質設定装置とを備え、
前記流入水変動係数演算装置は、日変動数値を作成する日変動演算装置と、最初沈殿池の上流側に接続された管路に蓄積した物質の量を物質濃度として入力設定するための管路蓄積物質設定装置と、前記管路内に蓄積または流出される物質の濃度変化特性を表す所定のパラメータの値を入力設定するための管路パラメータ設定装置と、予想降雨量設定装置と、前記各設定装置の設定値に基づき、管路に蓄積した物質が未来の流量および水質に及ぼす影響を所定の演算式により演算する降雨影響演算装置とを備えたことを特徴する水質シミュレーション装置。
【請求項２】
請求項１に記載の水質シミュレーション装置において、前記水質データ及び運転状態記憶装置に対する水質データの入力値は、オンライン水質測定装置の測定値とすることを特徴とする水質シミュレーション装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、都市下水や産業廃水のように、汚水を浄化する活性汚泥法、嫌気・好気汚泥処理法などの汚水廃水処理プラントに利用される水質シミュレーション装置に関する。
【０００２】
【従来の技術】
近年、閉鎖性水域の富栄養化が問題となっている。これまで、ＣＯＤの総量規制により閉鎖性水域に流入する有機物量を減らし水質向上に向けた施策が実施されてきたが、これまでの取り組みだけでは、目標とした水質基準を十分達成できない。水域中に存在する窒素やりんは、プランクトンの栄養源となり赤潮被害の原因となることが知られており、そこで、水質保全のために水域へ流入する窒素とりんの総量規制が実施されようとしている。
【０００３】
これに伴い有機物のみならず窒素やりんが除去できる、より高度な水質処理が求められている。そして、窒素やりん除去が可能な下水処理プロセスの開発や下水処理プラント等の導入が進められている。前述のような下水処理プロセスにおいては、水中の有機物の状況や気候（特に雨水の流入）などから水質を予測して運転方法を決定する必要がある。
【０００４】
下水処理施設の運転は、流入水量や水質が経時的に変動し、かつ水理学的滞留時間が長いプラントであることから運転制御が難しく、これまで、運転員の経験と勘に基づいてポンプなどの操作量を決定している。一方で、近年International Water Associationから発表された活性汚泥のシミュレーション技術を使って下水処理場の運転状態を模擬するものがある。即ち、下水処理プラントの運転状況の把握や運転方針の検討を行うために、シミュレータが使用されることがある。
【０００５】
このシミュレータを使用するためには、以下の４項目に関して事前に決定する必要がある。まず1つ目は、下水処理プラントの水槽の数と水路の設定やポンプの流量や溶存酸素濃度や温度などの運転条件である。２つ目は、測定した流入水質（全TOC濃度、溶解性TOC濃度、全CODcr濃度、溶解性CODcr濃度、全窒素濃度、溶解性窒素濃度、全りん濃度、溶解性りん濃度、りん酸態りん濃度、硝酸態窒素濃度、アンモニア態窒素濃度）をシミュレータが物質として扱う公知の状態変数（SI、SA、SF、SPO4、SNH4、SNO3、SN2、XFeP、XS、XAUT、XPHA、XPP、XPAO、XH、XI）に分画する必要がある。３つ目は、活性汚泥が行う生物反応の反応速度や収率を示すパラメータを設定する必要がある。４つ目は、各水槽内の物質の初期値を設定する必要がある（Henze, M, et al.：Activated Sludge Model No.2、IAWQ Scientific and Technical Report No.3(1995)参照）。
【０００６】
特に、前記２つ目の項目は、シミュレーションを開始する段階で、入力値はその後のシミュレーションの精度に大きく影響する。
【０００７】
下水道の流入水質データを提供する先行技術としては、日変動データを算出するための係数と日代表値の積から日変動データを計算する手法や、雨天日データを作成するために雨天日の係数と水質分析データから雨天日の水質データを算出する手法が公表されている。
【０００８】
例えば、下記特許文献１は、種々の水質シミュレーションを高精度に行なうことを目的として、図６（特許文献１の図１の部番を一部変更した図）に示すような装置を開示している。図６は、水質シミュレーション装置とその周辺との接続状況を示すブロック図であり、本件発明と共通する構成部分も多いので、本図について、特許文献１の記載に基づいて、以下にその概要を述べる。
【０００９】
図６において、１は最初沈殿池、２は反応槽、３は最終沈殿池である。先ず、下水処理プラントの概要について述べる。下水処理プラントでは、流入した下水は、最初沈殿池１、反応槽２、最終沈殿池３の順に流れて処理される。汚濁物質を含む下水は、最初沈殿池１に導入され、汚濁物質の中の沈降しやすいものを沈降分離して上澄水を反応槽２に流出する。反応槽２には最終沈殿池３の汚泥の一部が返送汚泥ポンプによって返送されており、反応槽２はその返送汚泥と最初沈殿池１の上澄水を処理する。反応槽２では、ブロワー（図示せず）から圧送された空気が曝気槽内の散気管によって放出されており、汚濁物質は活性汚泥により吸着、分解されて最終沈殿池３に導かれる。最終沈殿池３では活性汚泥を沈降分離し、沈降汚泥は余剰汚泥ポンプ（図示せず）により汚泥処理系（図示せず）に排出され、清澄水は処理水として滅菌槽（図示せず）を経て放流される。
【００１０】
次に、図６における水質シミュレーション装置について述べる。図６のシミュレーション装置は、特許文献１の記載によれば、「下水処理場の水質データを蓄積するデータ蓄積装置１ａと、土木構造などの処理場仕様を設定する処理場仕様設定装置２ａと、運転条件を設定する運転条件設定装置３ａと、水質モデルパラメータを設定するモデルパラメータ設定装置４ａと、日変動データを算出するための係数を設定する係数設定装置７ａと、係数設定装置の係数と日代表値との積から日変動データを計算する日変動データ計算装置８ａと、各装置からのデータを受け取り処理場の水質を予測するシミュレーション計算装置９ａとを備えており、データ蓄積装置１ａと係数設定装置７ａとの間に日変動の係数を自動的に演算する係数演算装置５ａと、データ蓄積装置１ａからのデータに基づき雨天日の水質データを作成しシミュレーション計算装置に送る雨天日データ作成装置６ａとを設けた構成」とし、
さらに、「前記係数演算装置５ａは、日変動のパターン係数を計算する基本係数演算装置５１と補間係数を計算する補間係数演算装置５２とを有し、また、前記雨天日データ作成装置６ａは、雨天日の雨量によってパターン係数を計算する雨天日係数作成装置６１および雨天日データを計算する雨天日データ計算装置６２とを有する構成とした」ことを特徴としている。
【００１１】
上記構成により、「日変動のデータおよび雨天日のデータを数分単位で作成することができ、種々の状況に対応したシミュレーションを極めて高精度に行える効果がある」旨、特許文献１に記載されている。（詳細は、特許文献１参照）。
【００１２】
また、特許文献１とは異なる先行技術として、下記特許文献２は、「流入水質の計測方法としてオンラインＵＶ計を用い、流入量が既定値を超えた場合、雨天モードと判断し、流入水質濃度の換算方法を切り替える手法」を開示している（詳細は、特許文献２参照）。
【００１３】
【特許文献１】
特開２００２−１３７０号公報（第１〜２頁、第５頁、図１）
【特許文献２】
特開２０００−１０７７９６号公報（第５頁、図７）
【００１４】
【発明が解決しようとする課題】
上述した特許文献１および２の方法によれば、運転員の経験や勘によることなく、晴天日と雨天日の水質データの切替が自動的に設定できるようになるため運転員の負担軽減を図ることができるものの、下記のような問題点がある。
【００１５】
現実の降雨時流入水質変化は、降雨初期には下水管や地表に蓄えられた汚濁物質が一気に流入するため流入水質濃度が上昇し、降雨が続くと雨による希釈効果で流入水質濃度が低下する。前記特許文献１における実施形態の説明によれば、雨天時の水質は、晴天時の水質に雨天日と晴天日の水質データの差を加えているが、雨の降り始めからの経過時間に関して考慮されていないので降雨が継続した場合の流入水質変化が現実と乖離する部分が生ずる問題がある。
【００１６】
前記特許文献２において、降雨モードと晴天モードを切り替える理由は、雨天の場合、流入量の増加によって下水管や地表に蓄えられた汚濁物質が一気に流入するためＵＶ計データとＢＯＤ濃度の相関が変化し、下水処理シミュレータへの流入有機物濃度に誤差が生じるためである。この方法でも同様に、降雨が継続した場合の雨の希釈効果による流入水質濃度低下を再現できない問題がある。
【００１７】
この発明は、上記従来の問題点に鑑みてなされたもので、この発明の課題は、諸設定を自動的に行い使用者の負担を軽減すると共に、降雨時の流入水質変化の現実とシミュレーション演算結果との誤差を低減させ、精度の高いシミュレーションを可能とした水質シミュレーション装置を提供することにある。
【００１８】
【課題を解決するための手段】
前述の課題を解決するために、この発明は、水質データ及び運転状態記憶装置と、運転条件を設定する運転条件設定装置と、プラントの化学反応に関する定数を設定するパラメータ設定装置と、処理場の土木構造やプロセスフローを設定するプラント土木条件設定装置と、未来の流入水状態を設定する未来データ設定装置と、前記各設定装置の設定値に基づき、プロセス各部の水質予測演算を行なうシミュレーション演算装置とを備えた水質シミュレーション装置において、
前記未来データ設定装置は、過去の運転データに基づき変動する流入水状態を予測演算する流入水変動係数演算装置を有し、この演算結果と現時点のデータとに基づき未来の流入水量と水質データとを設定するものとし、
さらに、前記未来データ設定装置の設定値に基づき、流入水の有機成分としては、全TOC濃度と溶解性TOC濃度、または全CODcr濃度と溶解性CODcr、窒素成分としては、全窒素濃度と、溶解性窒素濃度と硝酸態窒素濃度とアンモニア態窒素濃度、りん成分としては、全りん濃度と溶解性りん濃度とりん酸態りん濃度、および浮遊物質濃度の測定結果と一致するようにシミュレーションで定義している有機物への分画比やその分画された各有機物に含まれる物質の含有比を演算する流入水質演算装置と、前記分画された各有機物に含まれる物質の含有比を前記シミュレーション演算装置に設定する流入水質設定装置とを備え、
前記流入水変動係数演算装置は、日変動数値を作成する日変動演算装置と、最初沈殿池の上流側に接続された管路に蓄積した物質の量を物質濃度として入力設定するための管路蓄積物質設定装置と、前記管路内に蓄積または流出される物質の濃度変化特性を表す所定のパラメータの値を入力設定するための管路パラメータ設定装置と、予想降雨量設定装置と、前記各設定装置の設定値に基づき、管路に蓄積した物質が未来の流量および水質に及ぼす影響を所定の演算式により演算する降雨影響演算装置とを備えたものとする（請求項１の発明）。
【００１９】
上記未来データ設定装置および流入水質設定装置を備えた構成によれば、降雨が継続した場合の流入水質変化の現実との誤差が低減され、精度の高いシミュレーションを可能となる。なお、上記において、管路パラメータ設定装置における所定のパラメータとは、後述するQstayinmax、Qstayoutmax、V1、V2、Kstayin、Kstayoutである。また、上記降雨影響演算装置における所定の演算式とは、後述する［数１］である。詳細は、実施例とともに詳述する。
【００２０】
前記請求項１の発明の実施態様としては、下記請求項２の発明が好ましい。即ち、請求項１に記載の水質シミュレーション装置において、前記水質データ及び運転状態記憶装置に対する水質データの入力値は、オンライン水質測定装置の測定値とする（請求項２の発明）。
【００２１】
上記オンライン水質測定装置により、例えば３０分毎に計測すれば、一日１回もしくは数回の計測に比べて、降雨影響が精度よく反映された自動化水質シミュレーションが可能となる。
【００２２】
【発明の実施の形態】
図面に基づき、本発明の実施例について以下にのべる。
【００２３】
図１は、本発明の水質シミュレーション装置５の構成と、対象としている下水処理プラントとの接続状況を示した実施例の構成図である。図１において、図６に示した部材と同一機能を有する部材には、同一番号を付してその詳細説明を省略する。
【００２４】
図１において、４はオンライン水質測定装置であり、サンプリング地点a、b、c、dの水質をオンラインで測定するものである。測定項目としては、全TOC濃度、溶解性TOC濃度、全CODcr濃度、溶解性CODcr濃度、全窒素濃度、溶解性窒素濃度、全りん濃度、溶解性りん濃度、りん酸態りん濃度、硝酸態窒素濃度、アンモニア態窒素濃度、浮遊物質濃度などである。
【００２５】
水質データおよび運転状態記憶装置６は、前記の水質、プラントに流入する水量、各種ポンプ流量、水温などのデータを蓄積する。未来データ設定装置７は、流入水変動係数演算装置８により演算された未来の流入水量データと水質データとを流入水質演算装置９およびシミュレーション演算装置１４に設定する。
【００２６】
前記流入水変動係数演算装置８は、過去の運転データや予想降雨量から未来の流量や水質を予測する演算装置である。流入水変動係数演算装置８において、８１は、日変動演算装置であり、過去の運転データから１週間毎の流入水量や水質を水質データおよび運転状態記憶装置６から抽出し、移動平均を計算する演算装置である。８２は、管路蓄積物質設定装置であり、最初沈殿池１の上流側に接続された管路に蓄積した物質の量を物質濃度として、装置の使用者が入力設定するための装置である。
【００２７】
８３は、管路パラメータ設定装置であり、前記管路内に蓄積または流出される物質の濃度変化特性を表す所定のパラメータの値を入力設定するための装置である。所定のパラメータとは、後述するQstayinmax、Qstayoutmax、V1、V2、Kstayin、Kstayoutであり、装置の使用者が入力する。この管路パラメータの設定により、管路内に蓄積または流出される物質の濃度が変化する。８４は、予想降雨量設定装置であり、シミュレーション演算装置１４により計算させる期間の降雨量を天気予報等の情報から設定する装置である。８５は、降雨影響演算装置であり、未来の時系列データを算出する装置である。即ち、降雨影響演算装置８５は、前記日変動演算装置８１から、後述するQin,CINiを設定し、前記管路蓄積物質設定装置８２からC1i,C2i濃度を設定し、前記管路パラメータ設定装置８３からQstayinmax,Qstayoutmax,V1,V2,Kstayin,Kstayoutを設定し、予想降雨量設定装置８４からQrainを設定し、後述する［数１］に従って演算するものである。なお、iは、水質項目の違いを表している。具体的な水質項目は、全TOC濃度、溶解性TOC濃度、全CODcr濃度、溶解性CODcr濃度、全窒素濃度、溶解性窒素濃度、全りん濃度、溶解性りん濃度、りん酸態りん濃度、硝酸態窒素濃度、アンモニア態窒素濃度、浮遊物質濃度などである。
【００２８】
未来の時系列データは、後述する連立微分方程式から計算できる。この連立微分方程式においては、最初沈殿池１の上流側に接続された管路を図４に示すように模擬している。即ち、最初沈殿池１には、下廃水と雨水とが、仮想管路１８を経由して最初沈殿池１に流入するが、仮想管路１８中の流れが淀んだ部分を仮想管路滞留域１９として設定して、模擬管路を構成している。この仮想管路滞留域１９内の水質は、常時は淀んでいるために濃度が高いが、降雨により、徐々に淀みがなくなり、仮想管路１８内の水と仮想管路滞留域１９内の水とが混ざり合って、濃度が低減し、仮想管路１８内の濃度に近づく。詳細は、後述するデータにより述べる。
【００２９】
前記未来の時系列データの計算は、下記の連立微分方程式から計算できる。
【００３０】
【数１】
0004026057.tif 000002
ただし、前記式において、各記号の意味は次の通りである。
C1i：仮想管路１８の水質iの濃度(g/m3)
C2i：仮想管路滞留域１９の水質iの濃度(g/m3)
CINi：下廃水の水質iの濃度(g/m3)
V1：仮想管路１８の体積(m3)
V2：仮想管路滞留域１９の体積(m3)
Qin：下廃水の流量(m3/h)
Qrain：雨水の流量(m3/h)
Qstayin：仮想管路１８から仮想管路滞留域１９への流量(m3/h)
Qstayout：仮想管路滞留域１９から仮想管路１８への流量(m3/h)
Qstayinmax：仮想管路１８から仮想管路滞留域１９への飽和流量(m3/h)
Qstayoutmax：仮想管路滞留域１９から仮想管路１８への飽和流量(m3/h)
Kstayin：Qstayinの飽和係数
Kstayout：Qstayoutの飽和係数
前記［数１］において、式（１）は図４における仮想管路１８の水質を、式（２）は仮想管路滞留域１９の水質変化を数式で表現したものである。
【００３１】
次に、図１において、９は、流入水質演算装置であり、水質データおよび運転状態記憶装置６に蓄積された過去の水質の時系列データと、未来データ設定装置７により設定された未来の時系列データとを基にして、図２に示した全CODcr濃度や全窒素濃度や全りん濃度やりん酸態りん濃度や硝酸態窒素濃度やアンモニア態窒素濃度などの流入水質もしくは、図３に示した全TOC濃度や全窒素濃度や全りん濃度やりん酸態りん濃度や硝酸態窒素濃度やアンモニア態窒素濃度などの流入水質から、シミュレーションで定義されている各有機物質への分画比と、各有機物質にふくまれる窒素やりんの含有比とを最小二乗法などの最適化手法を用いて全有機物濃度と全窒素濃度と全りん濃度を近似的に一致させる装置である。
【００３２】
なお、この手法および装置は、公知の技術であり、図２および図３に示す記号は、前述のように、シミュレータが物質として扱う公知の状態変数（SI、SA、SF、SPO4、SNH4、SNO3、SN2、XFeP、XS、XAUT、XPHA、XPP、XPAO、XH、XIなど）である。
【００３３】
次に、図１における流入水質設定装置１０は、流入水質演算装置９で演算した数値をシミュレーション演算装置１４に設定する。運転条件設定装置１１は、返送率、SRT、制御目標溶存酸素濃度など、水処理プロセスを運転する上で不可欠な運転条件を設定する。パラメータ設定装置１２は、シミュレーションで使用するモデルの反応速度定数や化学量論係数などのパラメータを設定する。シミュレーション演算装置１４は、未来データ設定装置７、流入水質設定装置１０、運転条件設定装置１１、パラメータ設定装置１２、プラント土木条件設定装置１３からのデータを入力値として水処理プラントの各部の水質を予測する。
【００３４】
次に、図５に、図１の実施例によりシミュレーションした演算例として、未来の時系列データのグラフの一例を示す。図５は、横軸に時間を示し、縦軸には、最初沈殿池に流入する水質濃度（C1濃度）、具体的にはCODcr濃度と、降雨量（Qrain）の経時変化を示す。図示（ａ）で示すC1濃度は、（ｂ）で示す降雨の初期には下水管や地表に蓄えられた汚濁物質が一気に流入するため流入水質濃度が上昇し、降雨が続くと雨による希釈効果で流入水質濃度が低下し、降雨の終焉に伴って降雨開始前の濃度に戻る経過を示している。
【００３５】
なお、図５においては図示を省略しているが、図４における仮想管路滞留域１９の水質濃度（C2濃度）は、降雨開始後、仮想管路滞留域１９の水が仮想管路１８の水と混合して、徐々にC1濃度に接近する。
【００３６】
【発明の効果】
上記のとおり、この発明によれば、水質データ及び運転状態記憶装置と、運転条件を設定する運転条件設定装置と、プラントの化学反応に関する定数を設定するパラメータ設定装置と、処理場の土木構造やプロセスフローを設定するプラント土木条件設定装置と、未来の流入水状態を設定する未来データ設定装置と、前記各設定装置の設定値に基づき、プロセス各部の水質予測演算を行なうシミュレーション演算装置とを備えた水質シミュレーション装置において、
前記未来データ設定装置は、過去の運転データに基づき変動する流入水状態を予測演算する流入水変動係数演算装置を有し、この演算結果と現時点のデータとに基づき未来の流入水量と水質データとを設定するものとし、
さらに、前記未来データ設定装置の設定値に基づき、流入水の有機成分としては、全TOC濃度と溶解性TOC濃度、または全CODcr濃度と溶解性CODcr濃度、窒素成分としては、全窒素濃度と溶解性窒素濃度と硝酸態窒素濃度とアンモニア態窒素濃度、りん成分としては、全りん濃度と溶解性りん濃度とりん酸態りん濃度、および浮遊物質濃度の測定結果と一致するようにシミュレーションで定義している有機物への分画比やその分画された各有機物に含まれる物質の含有比を演算する流入水質演算装置と、前記分画された各有機物に含まれる物質の含有比を前記シミュレーション演算装置に設定する流入水質設定装置とを備え、
前記流入水変動係数演算装置は、日変動数値を作成する日変動演算装置と、最初沈殿池の上流側に接続された管路に蓄積した物質の量を物質濃度として入力設定するための管路蓄積物質設定装置と、前記管路内に蓄積または流出される物質の濃度変化特性を表す所定のパラメータの値を入力設定するための管路パラメータ設定装置と、予想降雨量設定装置と、前記各設定装置の設定値に基づき、管路に蓄積した物質が未来の流量および水質に及ぼす影響を所定の演算式により演算する降雨影響演算装置とを備えたものとしたので、
降雨時の流入水質変化の現実とシミュレーション演算結果との誤差を低減させ、精度の高いシミュレーションを可能とした自動化水質シミュレーション装置を提供することができる。
【図面の簡単な説明】
【図１】  本発明の実施例に関わる水質シミュレーション装置の構成図
【図２】  各有機物質の分画比と含有する窒素やりんの比を求める模式図
【図３】  各有機物質の分画比と含有する窒素やりんの比を求める図２とは異なる模式図
【図４】  最初沈殿池の上流側に接続された管路の模擬管路構成図
【図５】  図１の実施例によりシミュレーションした演算例を示す図
【図６】  従来の水質シミュレーション装置の一例を示す構成図
【符号の説明】
１：最初沈殿池、２：生物反応槽、３：最終沈殿池、４：オンライン水質測定装置、５：水質シミュレーション装置、６：水質データおよび運転状態記憶装置、７：未来データ設定装置、８：流入水変動係数演算装置、９：流入水質演算装置、１０：流入水質設定装置、１１：運転条件設定装置、１２：パラメータ設定装置、１３：プラント土木条件設定装置、１４：シミュレーション演算装置、８１：日変動演算装置、８２：管路蓄積物質設定装置、８３：管路パラメータ設定装置、８４：予想降雨量設定装置、８５：降雨影響演算装置。
【図１】
0004026057.tif 000003
【図２】
0004026057.tif 000004
【図３】
0004026057.tif 000005
【図４】
0004026057.tif 000006
【図５】
0004026057.tif 000007
【図６】
0004026057.tif 000008