JPS649077B2 - - Google Patents
Info
- Publication number
- JPS649077B2 JPS649077B2 JP55017087A JP1708780A JPS649077B2 JP S649077 B2 JPS649077 B2 JP S649077B2 JP 55017087 A JP55017087 A JP 55017087A JP 1708780 A JP1708780 A JP 1708780A JP S649077 B2 JPS649077 B2 JP S649077B2
- Authority
- JP
- Japan
- Prior art keywords
- model
- concentration
- treatment facility
- sewage treatment
- water flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired
Links
- 239000010865 sewage Substances 0.000 claims description 34
- 239000010802 sludge Substances 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 27
- 238000005273 aeration Methods 0.000 claims description 26
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 12
- 229910052760 oxygen Inorganic materials 0.000 claims description 12
- 239000001301 oxygen Substances 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 12
- 238000001125 extrusion Methods 0.000 claims description 10
- 239000000758 substrate Substances 0.000 claims description 9
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000010206 sensitivity analysis Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000813 microbial effect Effects 0.000 description 2
- 244000005700 microbiome Species 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
- Y02W10/00—Technologies for wastewater treatment
- Y02W10/10—Biological treatment of water, waste water, or sewage
Landscapes
- Activated Sludge Processes (AREA)
Description
【発明の詳細な説明】
本発明は、いわゆる活性汚泥法を採用している
下水処理場のごとく曝気槽を処理施設の主要装置
とする。下水処理施設の運転等のための電子計算
機シミユレーシヨンに用いられる下水処理施設の
シミユレーシヨン方法に関するものである。
下水処理場のごとく曝気槽を処理施設の主要装置
とする。下水処理施設の運転等のための電子計算
機シミユレーシヨンに用いられる下水処理施設の
シミユレーシヨン方法に関するものである。
活性汚泥法による下水処理施設は第１図に示す
ように構成されている。すなわち活性汚泥法にお
ける主要装置は曝気槽１と沈殿池２であり、曝気
槽１において流入下水中の基質を微生物の作用に
より汚泥に変換し、次に沈殿池２において沈降に
よる固液分離を行ない、その上澄み液を処理水３
として河川に放流するものである。なお、第１図
で４は返送汚泥管路、５は余剰汚泥管路、６は流
入下水管路、７は曝気用空気管路をそれぞれ示し
ている。前記流入下水の時間変化特性は日間変動
が大きいことと処理水の水質は所定の値以下に保
持されることが本処理施設の運転上の特徴であ
る。本処理施設の運転方法は運転技術者の経験に
負うところが多いが、一方では計算機シミユレー
シヨンを用いて運転状態の検討が行なわれる場合
も多い。後者すなわちシミユレーシヨンを行なう
場合、曝気槽と沈殿池の数学モデルを必要とする
が、主体は曝気槽であり、この曝気槽の数学モデ
ルは精度の良いものでなければならない。
ように構成されている。すなわち活性汚泥法にお
ける主要装置は曝気槽１と沈殿池２であり、曝気
槽１において流入下水中の基質を微生物の作用に
より汚泥に変換し、次に沈殿池２において沈降に
よる固液分離を行ない、その上澄み液を処理水３
として河川に放流するものである。なお、第１図
で４は返送汚泥管路、５は余剰汚泥管路、６は流
入下水管路、７は曝気用空気管路をそれぞれ示し
ている。前記流入下水の時間変化特性は日間変動
が大きいことと処理水の水質は所定の値以下に保
持されることが本処理施設の運転上の特徴であ
る。本処理施設の運転方法は運転技術者の経験に
負うところが多いが、一方では計算機シミユレー
シヨンを用いて運転状態の検討が行なわれる場合
も多い。後者すなわちシミユレーシヨンを行なう
場合、曝気槽と沈殿池の数学モデルを必要とする
が、主体は曝気槽であり、この曝気槽の数学モデ
ルは精度の良いものでなければならない。
曝気槽内の主な反応要因は微生物を含む活性汚
泥と基質の混合状態であり、通常用いられる指標
は生物化学的酸素要求量（基質濃度）BOD、浮
遊物濃度（汚泥濃度）SS、溶存酸素濃度DOの３
つである。混合状態のモデルは槽列モデルが一般
に広く用いられており、その段数を最小の１とし
たときに完全混合モデル、段数を無限大としたと
きに押出流モデルとなる。実際の混合状態はこれ
ら両極限の中間状態であると考えられており、多
段完全混合モデルとすることもある。汚泥増殖過
程を表わすものとしては種々の提案がなされてい
るが、ここでは微生物化学反応式としてMonod
の式を基本とすれば、完全混合モデルおよび押出
流モデルは次の各式で表わすことができる。
泥と基質の混合状態であり、通常用いられる指標
は生物化学的酸素要求量（基質濃度）BOD、浮
遊物濃度（汚泥濃度）SS、溶存酸素濃度DOの３
つである。混合状態のモデルは槽列モデルが一般
に広く用いられており、その段数を最小の１とし
たときに完全混合モデル、段数を無限大としたと
きに押出流モデルとなる。実際の混合状態はこれ
ら両極限の中間状態であると考えられており、多
段完全混合モデルとすることもある。汚泥増殖過
程を表わすものとしては種々の提案がなされてい
るが、ここでは微生物化学反応式としてMonod
の式を基本とすれば、完全混合モデルおよび押出
流モデルは次の各式で表わすことができる。
() 完全混合モデル
() 押出流モデル
ここに、Xc、XpとSc、Spはそれぞれ汚泥濃度
（SS）と基質濃度（BOD）で添字ＣとＰは各々
完全混合モデル（）と押出流モデル（）の場
合の値を表わす。XIとSIはそれぞれ入力の汚泥濃
度（SS）と基質濃度（BOD）、μ^は最大増殖率、
Ksは飽和定数、Ｙは収率係数、Kdは自己酸化率、
TRは曝気槽滞留時間、ｔは任意時刻である。
（SS）と基質濃度（BOD）で添字ＣとＰは各々
完全混合モデル（）と押出流モデル（）の場
合の値を表わす。XIとSIはそれぞれ入力の汚泥濃
度（SS）と基質濃度（BOD）、μ^は最大増殖率、
Ksは飽和定数、Ｙは収率係数、Kdは自己酸化率、
TRは曝気槽滞留時間、ｔは任意時刻である。
また、溶存酸素（DO）モデルとしては、次式
が一般的である。
が一般的である。
ここで、Ｃは溶存酸素濃度、Csは飽和溶存酸素
濃度、CIは溶存酸素濃度の入力値、KLaは総括酸
素移動係数、Ｆは曝気風量、rrは酸素消費速度で
あり、K1〜K4は係数である。ＳとＸは上記完全
混合モデル（）および押出流モデル（）にお
いて、それぞれScとXcおよびSpとXpとなる。
濃度、CIは溶存酸素濃度の入力値、KLaは総括酸
素移動係数、Ｆは曝気風量、rrは酸素消費速度で
あり、K1〜K4は係数である。ＳとＸは上記完全
混合モデル（）および押出流モデル（）にお
いて、それぞれScとXcおよびSpとXpとなる。
沈殿池のモデルは曝気槽から流入する流量と汚
泥濃度の関数として沈殿池引抜汚泥濃度が計算さ
れるもので、物理現象のモデル化または統計的モ
デル化によつて得られる。
泥濃度の関数として沈殿池引抜汚泥濃度が計算さ
れるもので、物理現象のモデル化または統計的モ
デル化によつて得られる。
以上のように下水処理場の流入下水特性の大き
な変動を考慮した非定常モデルすなわち(1)〜(3)式
で用いられている各係数（μ^、Ks、Ｙ、Kd、K1
〜K4）を時変としたモデルを用いるシミユレー
シヨンは従来行なわれておらず、概略演算による
シミユレーシヨンしか行なわれていなかつた。こ
のため、従来はシミユレーシヨンにより精度の高
い指標値を得ることができず、対象システムの
種々の運転状態に対応する処置方法の充分な検討
を行なうことができなかつた。
な変動を考慮した非定常モデルすなわち(1)〜(3)式
で用いられている各係数（μ^、Ks、Ｙ、Kd、K1
〜K4）を時変としたモデルを用いるシミユレー
シヨンは従来行なわれておらず、概略演算による
シミユレーシヨンしか行なわれていなかつた。こ
のため、従来はシミユレーシヨンにより精度の高
い指標値を得ることができず、対象システムの
種々の運転状態に対応する処置方法の充分な検討
を行なうことができなかつた。
本発明は、このような事情を背景としてなされ
たもので、曝気槽を処理施設の主要装置とする下
水処理施設のシミユレーシヨンを行なうに当り、
日間変動の大きな流入下水特性を考慮し合理的で
且つ高精度の指標値を得ることを可能とする下水
処理施設のシミユレーシヨン方法を提供すること
を目的としている。
たもので、曝気槽を処理施設の主要装置とする下
水処理施設のシミユレーシヨンを行なうに当り、
日間変動の大きな流入下水特性を考慮し合理的で
且つ高精度の指標値を得ることを可能とする下水
処理施設のシミユレーシヨン方法を提供すること
を目的としている。
すなわち、本発明の特徴とするところは、日間
変動の大きな流入下水特性を考慮し、曝気槽の数
学モデルを、流入下水流量の範囲を複数分割し、
これら分割された互いに異なる水量の複数モード
に対して個別の混合モデルを対応させ、これらを
時系列的に連結する非定常モデルとすることにあ
る。
変動の大きな流入下水特性を考慮し、曝気槽の数
学モデルを、流入下水流量の範囲を複数分割し、
これら分割された互いに異なる水量の複数モード
に対して個別の混合モデルを対応させ、これらを
時系列的に連結する非定常モデルとすることにあ
る。
以下、本発明の一実施例について詳細に説明す
る。
る。
流入下水特性を考慮した曝気槽の数学モデルを
構築しようとする場合、上記(1)〜(3)式における８
係数（μ^、Ks、Ｙ、Kd、K1〜K4）を時変とする
非定常モデルが考えられるが、これは現状では理
論的にも実験的にも実現不可能である。そこで、
流入下水流量を高水量と低水量との２領域に分
け、それぞれの混合状態が異なるものと考えて
各々個別の混合モデルを対応させ、これらを時系
列的に連結させることにした。
構築しようとする場合、上記(1)〜(3)式における８
係数（μ^、Ks、Ｙ、Kd、K1〜K4）を時変とする
非定常モデルが考えられるが、これは現状では理
論的にも実験的にも実現不可能である。そこで、
流入下水流量を高水量と低水量との２領域に分
け、それぞれの混合状態が異なるものと考えて
各々個別の混合モデルを対応させ、これらを時系
列的に連結させることにした。
基本的な混合モデルとしては、曝気槽の実際の
混合状態を上述の両極限状態モデルである単段完
全混合モデルと押出流モデルとを重み係数により
結合させたモデルによつて表わすことにする。
混合状態を上述の両極限状態モデルである単段完
全混合モデルと押出流モデルとを重み係数により
結合させたモデルによつて表わすことにする。
この場合、濃度Ｚは次式で表わされる。
Ｚ＝α・Zc＋（１−α）Zp ……(4)
ここで、濃度ＺはＸまたはＳである。また、α
は重み係数であり、任意時刻ｔにおける重み係数
α（ｔ）は次式で表わされる。
は重み係数であり、任意時刻ｔにおける重み係数
α（ｔ）は次式で表わされる。
α（ｔ）＝αH when Ｑ（ｔ）≧QM
αL when Ｑ（ｔ）＜QM ……(5)
ここで、Ｑ（ｔ）は流入下水流量、QMは流水下
水流量の使用実績データ群Ｑの中央値である。
水流量の使用実績データ群Ｑの中央値である。
以上により、(1)〜(3)式の８係数は一定とし、こ
れらの時間変化分を、重み係数を導入しこの重み
係数を流量との関係で２値とすることによつて処
理する簡易形非定常モデルを得たことになる。こ
のモデルの利点は、多段完全混合モデルを用いた
場合には流量が高水量のときと低水量のときとで
異なる段数となつてしまい、時系列的に連結する
ときに各段における初期値の計算が必要となるの
に対し、単段であるためにこの初期値計算が不要
であり、演算時間が短縮されることになるという
点である。
れらの時間変化分を、重み係数を導入しこの重み
係数を流量との関係で２値とすることによつて処
理する簡易形非定常モデルを得たことになる。こ
のモデルの利点は、多段完全混合モデルを用いた
場合には流量が高水量のときと低水量のときとで
異なる段数となつてしまい、時系列的に連結する
ときに各段における初期値の計算が必要となるの
に対し、単段であるためにこの初期値計算が不要
であり、演算時間が短縮されることになるという
点である。
次に、上述のモデルの構築のために、(1)〜(3)式
の８係数と重み係数の２値の合計10個の値を決定
しなければならない。
の８係数と重み係数の２値の合計10個の値を決定
しなければならない。
第２図は、これら係数の決定手順を示す流れ図
であり、以下この流れ図に沿つて説明する。
であり、以下この流れ図に沿つて説明する。
最初は重み係数の組み合わせである。この場
合、(1)〜(3)式をそのまま用いたのでは予め８係数
の値の決定することが必要となるが、曝気槽にお
ける実機データが示すところによると、汚泥濃度
の増殖項と自己酸化項の差は汚泥濃度に比して１
桁以上も小さな値となつていることに着目し、重
み係数の決定に際してはこれらの項を省略して汚
泥濃度のみのデータを用いて決定する。この場
合、完全混合モデルおよび押出流モデルとしては
次式を用いる。
合、(1)〜(3)式をそのまま用いたのでは予め８係数
の値の決定することが必要となるが、曝気槽にお
ける実機データが示すところによると、汚泥濃度
の増殖項と自己酸化項の差は汚泥濃度に比して１
桁以上も小さな値となつていることに着目し、重
み係数の決定に際してはこれらの項を省略して汚
泥濃度のみのデータを用いて決定する。この場
合、完全混合モデルおよび押出流モデルとしては
次式を用いる。
() 完全混合モデル
dXc／dt＝−（Xc−XI）／TR ……（1a）
() 押出流モデル
Xp＝XI（ｔ−TR） ……（2a）
重み係数は０〜１の値をとるものであるから、０
〜１の範囲をＮ等分すれば（Ｎ＋１）個の点が得
られる。２個の重み係数αH、αLの組合わせは（Ｎ
＋１）2個となり、この全数について、（1a）式と
（2a）式を用いた計算値と実測値とが最も良く適
合する場合を列挙法によつて選定する。その結果
得られた２係数の組合わせが最適係数の組合わせ
α* Hとα* Lである。
〜１の範囲をＮ等分すれば（Ｎ＋１）個の点が得
られる。２個の重み係数αH、αLの組合わせは（Ｎ
＋１）2個となり、この全数について、（1a）式と
（2a）式を用いた計算値と実測値とが最も良く適
合する場合を列挙法によつて選定する。その結果
得られた２係数の組合わせが最適係数の組合わせ
α* Hとα* Lである。
次に、上記(1)式と(2)式を用いて基質濃度
（BOD）と汚泥濃度（SS）に関する４係数を決
定するが、１つの係数についてＭ個の値を考える
ものとすれば、全数列挙法ではM4個の場合数と
なり、多すぎるので、逐次決定してゆくことにす
る。この場合、順序決定しなければならないの
で、４係数と個々の適合度を示す指標との感度解
析を行ない、感度の高い順に係数を決定する。こ
のようにすれば、4M個の場合数で４係数が決定
できる。ここで、適合度を示す指標I1は次式のご
とく定義する。
（BOD）と汚泥濃度（SS）に関する４係数を決
定するが、１つの係数についてＭ個の値を考える
ものとすれば、全数列挙法ではM4個の場合数と
なり、多すぎるので、逐次決定してゆくことにす
る。この場合、順序決定しなければならないの
で、４係数と個々の適合度を示す指標との感度解
析を行ない、感度の高い順に係数を決定する。こ
のようにすれば、4M個の場合数で４係数が決定
できる。ここで、適合度を示す指標I1は次式のご
とく定義する。
I1＝D
〓
〓
｛β（XA−XC／XA）2＋（１−β）・（SA−SC／SA）2
｝……(6) ここで、ＸとＳはそれぞれ汚泥濃度（SS）と
基質濃度（BOD）であり、添字ＡとＣはそれぞ
れ実測値と計算値を示している。βは重み係数、
Ｄはここで用いるデータの離散時間数である。
｝……(6) ここで、ＸとＳはそれぞれ汚泥濃度（SS）と
基質濃度（BOD）であり、添字ＡとＣはそれぞ
れ実測値と計算値を示している。βは重み係数、
Ｄはここで用いるデータの離散時間数である。
引き続いて、(3)式の４係数を決定する。これは
(1)式、(2)式の場合と同様な取り扱いによつて行な
うが、この場合、適合度を示す指標I2は次式を用
いる。
(1)式、(2)式の場合と同様な取り扱いによつて行な
うが、この場合、適合度を示す指標I2は次式を用
いる。
I2＝D
〓
（CA−CC／CA） ……(7)
ここで、Ｃは溶存酸素濃度であり、添字ＡとＣ
はそれぞれ実測値と計算値である。
はそれぞれ実測値と計算値である。
上述の係数決定のアルゴリズムを図示したもの
が第２図に示す流れ図である。
が第２図に示す流れ図である。
このようにして、10個の係数が決定されれば、
基質濃度（BOD）、汚泥濃度（SS）、溶存酸素濃
度（DO）の３指標を計算することが可能な曝気
槽非定常モデルが得られたことになる。
基質濃度（BOD）、汚泥濃度（SS）、溶存酸素濃
度（DO）の３指標を計算することが可能な曝気
槽非定常モデルが得られたことになる。
上述のようにして具体化された、下水処理場の
下水流入特性の大幅な変動を考慮した曝気槽の非
定常モデルを用い、第１図に示したシステムにお
ける流入下水６と返送汚泥４の時系列特性と曝気
用空気７量の時系列データを与えれば、曝気槽１
の出口における上記３指標値（BOD、SS、DO）
を逐次計算し、第３図ａ，ｂ，ｃにそれぞれ
MLSS（曝気槽内混合液浮遊物濃度（SS））、DO
（溶存酸素濃度）、BOD（生物化学的酸素要求量）
の時系列データの実測値（実線）と本実施例によ
る計算値（破線）の一例を示すように、高精度の
時系列データが得られる。なお、第３図ａ〜ｃで
は、MLSSの計算値は実測値とよく適合している
ことがよくわかるのに対し、DOとBODについて
は計算値と実測値の適合度がやや低いように見え
るが、これはDOとBODはMLSSに比して著しく
小さい値であるためにこのようになるのであつ
て、変化の傾向としてはよく一致していることが
わかる。この場合、沈殿池２のモデルを用いれ
ば、返送汚泥４の時系列特性が計算可能となる。
この対象システムにおける可変要因としては流入
下水特性、操作電動機や調整弁、さらに流入流量
分配比などがあるが、実機では実行不可能なこれ
ら可変要因の各種組合わせについて数値実験を試
みることができるので、各要因の各指標に対する
影響について定量的に把握することができる。
下水流入特性の大幅な変動を考慮した曝気槽の非
定常モデルを用い、第１図に示したシステムにお
ける流入下水６と返送汚泥４の時系列特性と曝気
用空気７量の時系列データを与えれば、曝気槽１
の出口における上記３指標値（BOD、SS、DO）
を逐次計算し、第３図ａ，ｂ，ｃにそれぞれ
MLSS（曝気槽内混合液浮遊物濃度（SS））、DO
（溶存酸素濃度）、BOD（生物化学的酸素要求量）
の時系列データの実測値（実線）と本実施例によ
る計算値（破線）の一例を示すように、高精度の
時系列データが得られる。なお、第３図ａ〜ｃで
は、MLSSの計算値は実測値とよく適合している
ことがよくわかるのに対し、DOとBODについて
は計算値と実測値の適合度がやや低いように見え
るが、これはDOとBODはMLSSに比して著しく
小さい値であるためにこのようになるのであつ
て、変化の傾向としてはよく一致していることが
わかる。この場合、沈殿池２のモデルを用いれ
ば、返送汚泥４の時系列特性が計算可能となる。
この対象システムにおける可変要因としては流入
下水特性、操作電動機や調整弁、さらに流入流量
分配比などがあるが、実機では実行不可能なこれ
ら可変要因の各種組合わせについて数値実験を試
みることができるので、各要因の各指標に対する
影響について定量的に把握することができる。
このように、上述の曝気槽非定常モデルを用い
てシミユレーシヨンを行なえば、精度の高い３指
標（BOD、SS、DO）の時系列データを得るこ
とが可能となるので、可変要因として流入下水特
性、操作電動機や調整弁、さらに流入流量配分比
などを種々の組合わせについて試行演算を行な
い、その結果を比較検討することによつて運転状
態として可能と思われるもの全てについて知るこ
とができる。これは設置されている処理施設の運
転において、指標的に問題となるような条件の場
合に対する処置方法も事前に検討することがで
き、長期間にわたつて安定な運転が可能となる。
また、実際に設置されている機器や装置について
も、数字上では容易に定格を変更することが可能
であるので、定格をこえる機器や装置の運転を試
行演算し、その結果から機器や装置の定格を改め
ること、すなわち機器や装置の新設計を行なうこ
とが可能となる。総じて、活性汚泥法の下水処理
施設の運転について、本シミユレーシヨンは実機
の運転に支障をきたすことなくその総合的把握を
するための極めて有用な一手段となる。
てシミユレーシヨンを行なえば、精度の高い３指
標（BOD、SS、DO）の時系列データを得るこ
とが可能となるので、可変要因として流入下水特
性、操作電動機や調整弁、さらに流入流量配分比
などを種々の組合わせについて試行演算を行な
い、その結果を比較検討することによつて運転状
態として可能と思われるもの全てについて知るこ
とができる。これは設置されている処理施設の運
転において、指標的に問題となるような条件の場
合に対する処置方法も事前に検討することがで
き、長期間にわたつて安定な運転が可能となる。
また、実際に設置されている機器や装置について
も、数字上では容易に定格を変更することが可能
であるので、定格をこえる機器や装置の運転を試
行演算し、その結果から機器や装置の定格を改め
ること、すなわち機器や装置の新設計を行なうこ
とが可能となる。総じて、活性汚泥法の下水処理
施設の運転について、本シミユレーシヨンは実機
の運転に支障をきたすことなくその総合的把握を
するための極めて有用な一手段となる。
なお、本発明は上述し且つ図面に示す実施例に
のみ限定されることなく、その要旨を変更しない
範囲内で種々変形して実施することができる。
のみ限定されることなく、その要旨を変更しない
範囲内で種々変形して実施することができる。
例えば、上記実施例では、(5)式における重み係
数α（ｔ）は２値としたが、これをｋ値（ｋ≧
３；ｋは正整数）すなわち流入下水流量の範囲を
３以上に分割し各々異なる重み係数を割り当てる
ようにすることも可能である。この場合には、０
〜１の範囲をＮ等分すれば、場合数は（Ｎ＋１）k
となり、ｋを大きくすればそれだけシミユレーシ
ヨン費用がかさむことになり実用上はこの費用に
制約が付されるので、あまり大きくすることは得
策でなくなる。
数α（ｔ）は２値としたが、これをｋ値（ｋ≧
３；ｋは正整数）すなわち流入下水流量の範囲を
３以上に分割し各々異なる重み係数を割り当てる
ようにすることも可能である。この場合には、０
〜１の範囲をＮ等分すれば、場合数は（Ｎ＋１）k
となり、ｋを大きくすればそれだけシミユレーシ
ヨン費用がかさむことになり実用上はこの費用に
制約が付されるので、あまり大きくすることは得
策でなくなる。
また、汚泥増殖過程の微生物化学反応式として
のMonodの式の変形は(1)、(2)式に示したもの以
外にも多種公表されているが、これら他の式を用
いることも可能である。この場合、決定すべき係
数の個数が変ることもあるが、感度解析の結果、
指標（(6)式）に対する感度の高い順に係数を決定
するという手順は全く同様である。
のMonodの式の変形は(1)、(2)式に示したもの以
外にも多種公表されているが、これら他の式を用
いることも可能である。この場合、決定すべき係
数の個数が変ることもあるが、感度解析の結果、
指標（(6)式）に対する感度の高い順に係数を決定
するという手順は全く同様である。
さらに適合度を示す指標は(6)、(7)式に示したも
の以外に、実測値と計算値との差の絶対値を用い
るなどしてもよい。
の以外に、実測値と計算値との差の絶対値を用い
るなどしてもよい。
以上詳述したように、本発明によれば、日間変
動の大きな流入下水特性を考慮し合理的で且つ高
精度の指標値を得ることを可能とする下水処理施
設のシミユレーシヨン方法を提供することができ
る。
動の大きな流入下水特性を考慮し合理的で且つ高
精度の指標値を得ることを可能とする下水処理施
設のシミユレーシヨン方法を提供することができ
る。
第１図は本発明方法の一実施例の適用される活
性汚泥法の下水処理施設の一例の概略構成を示す
概略構成図、第２図は同実施例における係数決定
アルゴリズムを示す流れ図、第３図ａ〜ｃは同実
施例を用いたシミユレーシヨン結果と実測値との
関係を示す図である。 １……曝気槽、２……沈殿池、３……処理水
（放流）、４……返送汚泥管路、５……余剰汚泥管
路、６……流入下水管路、７……曝気用空気管
路。
性汚泥法の下水処理施設の一例の概略構成を示す
概略構成図、第２図は同実施例における係数決定
アルゴリズムを示す流れ図、第３図ａ〜ｃは同実
施例を用いたシミユレーシヨン結果と実測値との
関係を示す図である。 １……曝気槽、２……沈殿池、３……処理水
（放流）、４……返送汚泥管路、５……余剰汚泥管
路、６……流入下水管路、７……曝気用空気管
路。
Claims (1)
- 【特許請求の範囲】 １ 下水処理施設における曝気槽の汚泥濃度、基
質濃度および溶存酸素濃度のうち少なくとも１つ
の濃度指標を求める方法において、 時間の経過によつて大きく変動する流入下水特
性について流入下水流量の範囲を高水量領域と低
水量領域とに分割すると共に予め定められた時間
ごとに流入下水流量が高水量領域か低水量領域か
を判断し、高水量領域の場合には単段完全混合モ
デルを用いて前記濃度指標を求め、前記低水量領
域の場合には押出流モデルを用いて前記濃度指標
を求めた後、これら濃度指標を重み係数により結
合してなる非定常モデルを用いて、シミユレーシ
ヨン作業を実行することを特徴とする下水処理施
設のシミユレーシヨン方法。 ２ 特許請求の範囲第１項記載の下水処理施設の
シミユレーシヨン方法において、重み係数を決定
するに当り、単段完全混合モデルと押出流モデル
について曝気槽出口側で計測した汚泥濃度のみの
データを用いて計算値と実測データとの適合性を
検討して列挙法により実測データに最も良く適合
するように決定することを特徴とする下水処理施
設のシミユレーシヨン方法。 ３ 特許請求の範囲第２項記載の下水処理施設の
シミユレーシヨン方法において、汚泥濃度と基質
濃度に関する単段完全混合モデルと押出流モデル
および溶存酸素モデルの少なくとも一方の係数を
決定するに当り、各係数の感度解析を行つて感度
の高い順に係数を決定するようにしたことを特徴
とする下水処理施設のシミユレーシヨン方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1708780A JPS56113395A (en) | 1980-02-14 | 1980-02-14 | Simulation of drainage disposer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1708780A JPS56113395A (en) | 1980-02-14 | 1980-02-14 | Simulation of drainage disposer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS56113395A JPS56113395A (en) | 1981-09-07 |
| JPS649077B2 true JPS649077B2 (ja) | 1989-02-16 |
Family
ID=11934196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1708780A Granted JPS56113395A (en) | 1980-02-14 | 1980-02-14 | Simulation of drainage disposer |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS56113395A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3030579U (ja) * | 1996-04-24 | 1996-11-01 | 株式会社長谷幸製作所 | ハンマー |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4866333B2 (ja) * | 2007-11-19 | 2012-02-01 | 住友化学株式会社 | 混合段数算出装置 |
| JP7122486B1 (ja) * | 2022-04-22 | 2022-08-19 | 株式会社 小川環境研究所 | 活性汚泥の処理水ｂｏｄ測定方法 |
-
1980
- 1980-02-14 JP JP1708780A patent/JPS56113395A/ja active Granted
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3030579U (ja) * | 1996-04-24 | 1996-11-01 | 株式会社長谷幸製作所 | ハンマー |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS56113395A (en) | 1981-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zaghloul et al. | Application of machine learning techniques to model a full-scale wastewater treatment plant with biological nutrient removal | |
| Sharafati et al. | The potential of new ensemble machine learning models for effluent quality parameters prediction and related uncertainty | |
| Olsson et al. | Wastewater treatment systems | |
| CN109205781B (zh) | 控制系统及控制方法 | |
| US6808630B2 (en) | System and method for AI controlling waste-water treatment by neural network and back-propagation algorithm | |
| Novotny et al. | Time series analysis models of activated sludge plants | |
| CN111762958A (zh) | 基于asm2d模型的污水处理厂深井曝气工艺优化方法及装置 | |
| Komulainen et al. | Estimation of effluent nutrients in municipal MBBR process | |
| JPS649077B2 (ja) | ||
| JP7494669B2 (ja) | 排水処理シミュレータのパラメータ値予測方法及び装置並びに下水処理プラントの制御方法及び装置 | |
| Herbay et al. | Water quality management with time varying river flow and discharger control | |
| Raut et al. | Determination of wastewater quality index of municipal wastewater treatment plant using fuzzy rule base | |
| Lee et al. | Dynamic modelling and simulation of activated sludge process using orthogonal collocation approach | |
| CN120449402A (zh) | 基于Modelica的机理-数据混合城市污水处理仿真模型的构建方法 | |
| JP4038718B2 (ja) | 下水処理プロセスのシミュレーション方法および下水処理プロセスシミュレータ | |
| Fiacco et al. | Sensitivity analysis of a nonlinear water pollution control model using an upper Hudson River data base | |
| Elshorbagy et al. | Development of an ASM1 dynamic simulation model for an activated sludge process in United Arab Emirates | |
| JP4016307B2 (ja) | 下水処理プロセスシミュレータ | |
| WO2024121918A1 (ja) | 水処理場運転支援システムおよび水処理場運転支援方法 | |
| Holenda | Development of modelling, control and optimization tools for the activated sludge process | |
| JPH07106357B2 (ja) | 水処理装置 | |
| Tyagi et al. | Neural model for the operational control of activated sludge processes | |
| Lacroix et al. | Computer simulation of activated sludge plant operation | |
| CN113867233B (zh) | 一种基于中试研究的颗粒污泥处理的控制方法及系统 | |
| Rehm et al. | Relationship of Soil Texture to Sulfur Oxidation 1 |