(57)【要約】【課題】培養槽における酸素移動速度を計算するためのシミュレータを提供する。【解決手段】本発明のシミュレータは、溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、を備える。【選択図】図１
2019129726.tif 000002
【特許請求の範囲】
【請求項１】
培養槽における酸素移動速度を計算するためのシミュレータであって、
溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、
飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、
通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、
反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、
酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、
を備えたことを特徴とするシミュレータ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＯＴＲ（Ｏｘｙｇｅｎ  Ｔｒａｎｓｆｅｒ  Ｒａｔｅ；酸素移動速度）をシミュレートするシステム、装置、プログラムに関するものである。
【背景技術】
【０００２】
製造現場の大型生産機などの培養槽内で、微生物や動植物細胞を培養するためには、培養槽内における各種パラメータを決定する必要がある。研究所などの小型の培養槽（ラボスケール）で得られた結果を基に、製造現場の大型生産機にも適用できるように、各種パラメータ（大型の培養槽の装置寸法や運転条件など）を決定することをスケールアップという。
【０００３】
培養槽が完全に同じ性能になるようにスケールアップすることは難しい。これは、例え幾何学的に相似になるよう培養槽を設計しても、液混合特性や物質移動特性などの特性を異なるスケールの装置同士で完全に一致させられないからである。
【０００４】
したがって、スケールアップにあたっては、何が重要な因子であるかを見抜き、少なくともその因子については小型培養槽と大型生産機で出来るだけ一致するように各種条件を設定することが重要である。
【０００５】
好気的な培養では、多くの場合、培養槽内の気相から液相への酸素の移動速度が細胞量や目的物質の生産性と非常に高く相関することが知られる。したがって、特に、好気的な培養では培養槽の酸素移動速度が重要なスケールアップ因子となる。
【０００６】
＜酸素移動速度ＯＴＲの計算式＞
培養槽における気相から液相への酸素移動速度は以下の式で示すことができる。
【０００７】
【数１】
2019129726.tif 000003
【０００８】
ＯＴＲ  ：酸素移動速度      ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ｋＬａ  ：酸素移動容量係数  ［ｈ−１］
ＤＯ＊  ：飽和溶存酸素濃度  ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯ    ：溶存酸素濃度      ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量      ［ｇ／ｍｏｌ］
【０００９】
ここで、ＯＴＲは、Ｏｘｙｇｅｎ  Ｔｒａｎｓｆｅｒ  Ｒａｔｅ（酸素移動速度）を意味しており、物理的な酸素溶解速度を示す。
好気的な培養で小型培養槽の結果を大型生産機で再現するためには、両者の酸素移動速度ＯＴＲを一致させることが極めて重要である。
【００１０】
＜スケールアップ＞
スケールアップの手段として、小型実験機と大型生産機との両方において、酸素移動容量係数ｋＬａを一致させることが、従来技術では行われてきた。しかしながら、数式（１）からもわかるように、酸素移動速度ＯＴＲを決める因子は、酸素移動容量係数ｋＬａの他にも、溶存酸素濃度ＤＯと、飽和溶存酸素濃度ＤＯ＊がある。これらの因子が全て一致しないと、スケールアップしたときに酸素移動速度ＯＴＲが一致しない。
【００１１】
スケールアップでよく行われるのは、更に、培養中に溶存酸素濃度をオンラインで測定し、それがラボスケールと一致するように撹拌速度等を制御する方法である。これを、溶存酸素濃度制御（ＤＯ制御）という。
【００１２】
スケールアップ時に、酸素移動容量係数ｋＬａを合わせたうえでＤＯ制御を行うと、数式（１）のうち酸素移動容量係数ｋＬａ、溶存酸素濃度ＤＯが一致するので、飽和溶存酸素濃度ＤＯ＊がおおむね同じであれば酸素移動速度ＯＴＲも一致することになる。
【００１３】
しかし、上記の方法では、溶存酸素濃度ＤＯを実測しながら調整する必要があり、所与の性能の培養槽において所望の溶存酸素濃度ＤＯが得られるかを事前に予測していない。つまり、実際にスケールアップを行ってみたら、所望の溶存酸素濃度ＤＯが維持できず、結局、所望の酸素移動速度ＯＴＲを実現できないという事態がありうる。
【００１４】
なお、撹拌速度等を変更してＤＯ制御を行った場合、酸素移動容量係数ｋＬａも変化しているはずであり、結果としてＤＯがラボスケールと一致したとしても酸素移動速度ＯＴＲが一致する保証はない。酸素移動速度ＯＴＲの観点からは、ＤＯ制御は必ずしもスケールアップにおいて有効な手段とは限らないことに注意が必要である。
【００１５】
酸素移動速度ＯＴＲは測定が可能であり（非特許文献１）、大型生産機と小型実験機の双方で各種の操作条件とＯＴＲの関係を予め測定しておくことで、小型実験機と酸素移動速度ＯＴＲが一致するような大型生産機の操作条件を探索することは技術的には可能である。
【００１６】
しかしながら、実生産で現に使用中の大型生産機では、多数の実験によって酸素移動速度ＯＴＲを測定することはコストや装置使用スケジュールの面から困難であり、理論計算によって酸素移動速度ＯＴＲを精度よく予測する方法が望まれる。また、新設の大型生産機を設計する段階でも、やはり、酸素移動速度ＯＴＲを事前に予測することが不可欠である。
【００１７】
以上の要望より、従来においても培養槽の酸素移動速度ＯＴＲの予測計算は行われてきた。しかしながら、酸素移動速度ＯＴＲには各種パラメータが相互に関係していて複雑な計算が要求されるために、様々な仮定をおいての簡易計算が行われてきたのが実情である。
【００１８】
酸素移動速度ＯＴＲの計算のためには、数式（１）の右辺に現れる各パラメータに数値を与える必要がある。例えば、非特許文献２では、溶存酸素濃度ＤＯ＝０と仮定している。また、非特許文献３では、溶存酸素濃度ＤＯ＝ＤＯｃｒｉｔ（細胞の生産活動のために許容できる最低値）と仮定している。さらに、非特許文献２、非特許文献３のいずれも、飽和溶存酸素濃度ＤＯ＊は通気ガスである大気の酸素濃度（２１％）と培養液にかかる圧力より計算している。
【００１９】
しかしながら、溶存酸素濃度ＤＯは、実際には必ずしも０（ゼロ）ではなく、一般に臨界溶存酸素濃度ＤＯｃｒｉｔよりも大きい値を維持するように培養操作が実施されることが多い。さらに、培養液として接触している気相の酸素濃度（排ガスの酸素濃度）は、細胞の呼吸活動によって酸素が消費される結果、通気ガスの酸素濃度とは一致しないので、気相の酸素濃度が通気ガスと同じ（通気ガスが大気の場合、２１％）であることを前提として、飽和溶存酸素濃度ＤＯ＊を計算することは適切ではない。従来技術では、以上のような適切ではない仮定の結果として、酸素移動速度ＯＴＲが過大に計算される。このため、実際にスケールアップを行ってみると大型生産機に適合する条件が得られない場合も少なくない。
【先行技術文献】
【非特許文献】
【００２０】
【非特許文献１】田口久治／永井史郎（１９８５）“微生物培養工学”（共立出版）Ｐ．１６７
【非特許文献２】川瀬義矩（１９９３）“生物反応工学の基礎”（化学工業社）ｐ．２２４−２２５
【非特許文献３】吉田敏臣（１９９８）“培養工学”（コロナ社）ｐ．５８−ｐ．６２
【非特許文献４】小林猛、本田裕之（２００２）”生物化学工学”（東京化学同人）ｐ．６２
【非特許文献５】石崎文彬（訳）（１９８８）“発酵工学の基礎”（学会出版センター）ｐ．１７６
【非特許文献６】日本化学会（編）（２００４）化学便覧改訂５版  基礎編ＩＩ（丸善）ｐ．ＩＩ−１４４
【非特許文献７】化学工学会（編）（１９９９）化学工学便覧第６版（丸善）ｐ．１９
【発明の概要】
【発明が解決しようとする課題】
【００２１】
本発明では、特に、培養液中の溶存酸素濃度ＤＯ、飽和溶存酸素濃度ＤＯ＊を事前に予測しつつ、ラボスケールからスケールアップ時に実現可能な酸素移動速度ＯＴＲを推定するシミュレータを提供することを目的のひとつとする。
【課題を解決するための手段】
【００２２】
気相に関する物質収支を解き、特に、排ガス酸素濃度を算出し、それに応じた飽和溶存酸素濃度ＤＯ＊を算出する。そして、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲ（Ｏｘｙｇｅｎ  Ｕｐｔａｋｅ  Ｒａｔｅ）との間で定常状態（酸素移動速度と酸素摂取速度が等しい）が成立すると見做せることを仮定したシミュレータを提供する。
【００２３】
具体的には、培養槽における酸素移動速度を計算するためのシミュレータは、溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、を備えればよい。
【発明の効果】
【００２４】
本発明のシミュレータによれば、溶存酸素濃度ＤＯ、飽和溶存酸素濃度ＤＯ＊を推定しつつ、与えられた酸素移動容量係数ｋＬａ下での酸素移動速度ＯＴＲを予測可能となる。これにより、スケールアップ時にも使用可能な酸素移動速度ＯＴＲが算出できるようになり、スケールアップの際の培養槽の条件をどのように設定すべきかを事前に検討することが可能となる。よって、スケールアップにおける試行錯誤を大幅に削減できる。
本発明の他の目的、特徴及び利点は添付図面に関する以下の本発明の実施例の記載から明らかになるであろう。
【図面の簡単な説明】
【００２５】
【図１】培養槽に酸素を供給する概念図を示す。
【図２】本発明の一実施例によるＯＴＲシミュレータの概念図を示す。
【図３】本発明の一実施例によるＯＴＲシミュレータの機能ブロック図を示す。
【図４】本発明の一実施例によるＯＴＲシミュレータのフローチャートを示す。
【図５】本発明の別の実施例によるＯＴＲシミュレータの機能ブロック図を示す。
【図６】酸素移動速度ＯＴＲをシミュレートするにあたって、各パラメータとの相関関係を示す。
【図７】本発明の一実施例によるＯＴＲシミュレータの活用例のイメージ（スケールアップ）を示す。
【図８】培養槽のサイズの影響について、本発明の一実施例によるＯＴＲシミュレータでの計算結果を示す。
【図９】同じ酸素移動容量係数ｋＬａであっても、通気量が異なる場合の影響について、本発明の一実施例によるＯＴＲシミュレータでの計算結果を示す。
【発明を実施するための形態】
【００２６】
＜酸素移動速度ＯＴＲの計算式＞
培養槽における気相から液相への酸素移動速度は、数式（１）（再掲）で示したように、酸素移動容量係数ｋＬａ、飽和溶存酸素濃度ＤＯ*、溶存酸素濃度ＤＯによって定まる。
【数２】
2019129726.tif 000004
【００２７】
ｋＬａは、酸素移動容量係数［ｈ−１］である。酸素移動速度ＯＴＲの大きさを決める定数であり、気相から液相への酸素の移動しやすさを表しており、装置形状、通気条件、撹拌条件、培養液物性によって決まると言われている。
【００２８】
＜飽和溶存酸素濃度ＤＯ＊＞
ＤＯ＊は、飽和溶存酸素濃度［ｇ−Ｏ２／Ｌ］であり、単位体積あたりの培養液内に溶存する酸素の最大限の濃度を示しており、その値はＨｅｎｒｙの法則に従って定まる。Ｈｅｎｒｙの法則とは、温度が一定のとき、一定量の液体に溶解する気体の質量はその気体の圧力（分圧）に比例するという法則である。すなわち、気相の酸素分圧が変化すると、飽和溶存酸素濃度ＤＯ＊も変化することになる。酸素の分圧は、気相の酸素濃度に依存する。
【００２９】
ここで、培養槽内の気相の酸素濃度が、単純に培養槽に導入される通気ガスの酸素濃度に等しければ、飽和溶存酸素濃度ＤＯ＊の計算は、それほど複雑にはならない。このため、従来は、飽和溶存酸素濃度ＤＯ＊は通気ガスの酸素濃度のみから定まるとして計算が行われていた。現実には、培養槽内の気相の酸素濃度は通気ガスの酸素濃度から変化するため従来法は必ずしも妥当ではない。
【００３０】
培養槽内の気相の酸素分圧は、通気ガスに含まれる酸素の量と気相から培養液に移動する酸素の量、気相に培養液から排出されたガス（一般に、排出ガスの大半は二酸化炭素ＣＯ２であるが、水蒸気も加わる）の量、および気相の全圧（培養槽が加圧されている場合にはその圧力、更に液深が大きい場合は静水圧が加わる）によって定まる。
【００３１】
要約すると、酸素移動速度ＯＴＲが決定されると通気量との関係により気相酸素濃度が決定される。
【００３２】
気相酸素濃度が決まると、飽和溶存酸素濃度ＤＯ＊が決まり、飽和溶存酸素濃度ＤＯ*との関係で酸素移動速度ＯＴＲが決定される。すなわち、酸素移動速度ＯＴＲと気相酸素濃度と飽和溶存酸素濃度ＤＯ＊とは依存関係を持ったループを形成している。これが飽和溶存酸素濃度ＤＯ＊の計算を非常に複雑にしている。
【００３３】
言い換えると、酸素移動速度ＯＴＲの計算のため、飽和溶存酸素濃度ＤＯ＊を事前に予測する上で最も問題なのは、飽和溶存酸素濃度ＤＯ＊を定める気相中の酸素分圧が酸素移動速度ＯＴＲの関数であり、かつ、酸素移動速度ＯＴＲが飽和溶存酸素濃度ＤＯ＊の関数であるという点である。
【００３４】
発明者は、この問題を解決して酸素移動速度ＯＴＲ、飽和溶存酸素濃度ＤＯ＊の推定値を得るには、酸素移動速度ＯＴＲと飽和溶存酸素濃度ＤＯ＊の関係を表現する非線形方程式につき、収束計算（繰り返し計算）を実施すれば良いことを見出した。
【００３５】
＜溶存酸素濃度ＤＯ＞
溶存酸素濃度ＤＯ（Ｄｉｓｓｏｌｖｅｄ  Ｏｘｙｇｅｎ）は、培養液内に溶存する酸素の濃度を示しており、単位は、［ｇ−Ｏ２／Ｌ］である。
【００３６】
酸素ガスを溶媒に溶かして溶存酸素にするために、通気撹拌を伴った培養槽（通気撹拌型培養槽）を設けてもよい。通気撹拌型培養槽とは、例えば、培養槽の底部に設けた管から酸素を含む気体（例えば空気）を供給して、培養槽の中央に下垂した撹拌軸の先に付けた羽根などの翼で培養液を撹拌するもので、酸素が微小な気泡になることにより、培養槽内に溶存することができる。
【００３７】
＜溶存酸素濃度ＤＯと酸素摂取速度ＯＵＲ＞
溶存酸素濃度ＤＯと酸素摂取速度ＯＵＲの関係について説明する。
培養液中の細胞は、培養液内に溶解している酸素を取り込むが、その活性は溶存酸素濃度ＤＯの関数で表現できる場合が多いことが知られる。単位細胞量当たりの酸素摂取速度（酸素比摂取速度）は、例えば以下のような酵素反応のミカエリス・メンテン式と同じ形式の数式（Ｍｏｎｏｄ式；モノー式）で表現できる。
【００３８】
【数３】
2019129726.tif 000005
【００３９】
ここで、
ＱＯ２      ：酸素比摂取速度      ［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
ＱＯ２,ｍａｘ：最大比酸素摂取速度  ［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
ＤＯ      ：溶存酸素濃度        ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
Ｋｍ，Ｏ２  ：飽和定数          ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
【００４０】
Ｋｍ，Ｏ２は飽和定数であり、最大酸素比摂取速度の１／２の速度を与える時の溶存酸素濃度に相当する。
【００４１】
酸素比摂取速度ＱＯ２を表現する式は、溶存酸素濃度ＤＯの関数であればどのような形でもよい。例えば、溶存酸素濃度ＤＯが臨界溶存酸素濃度ＤＯｃｒｉｔより低い場合、酸素比摂取ＱＯ２は溶存酸素濃度ＤＯと比例して減少する形でもよい（非特許文献４：小林猛、本田裕之（２００２）”生物化学工学”（東京化学同人）ｐ．６２参照）。
【００４２】
細胞の酸素摂取速度ＯＵＲは、酸素比摂取速度ＱＯ２を用いて、以下の式で示される。
【００４３】
【数４】
2019129726.tif 000006
【００４４】
ここで、
ＯＵＲ  ：酸素摂取速度      ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＱＯ２：酸素比摂取速度    ［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
Ｘ    ：細胞濃度          ［ｇ−ｄｒｙ  ｃｅｌｌ／Ｌ−ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量      ［ｇ／ｍｏｌ］
【００４５】
＜酸素摂取速度ＯＵＲと気相酸素濃度＞
酸素摂取速度ＯＵＲが定まると、更に、呼吸商（Ｒｅｓｐｉｒａｔｏｒｙ  Ｑｕｏｔｉｅｎｔ；ＲＱ）より、ＣＯ２生成速度（Ｃａｒｂｏｎ  ｄｉｏｘｉｄｅ  Ｅｖｏｌｕｔｉｏｎ  Ｒａｔｅ；ＣＥＲ）が定まる。呼吸商ＲＱとは、ある時間において生体内で栄養素が分解されてエネルギーに変換するまでの酸素消費量に対する二酸化炭素排出量の体積比のことである。呼吸率、呼吸係数とも呼ばれる。
【００４６】
【数５】
2019129726.tif 000007
ここで、
ＣＥＲ  ：ＣＯ２生成速度  ［ｍｍｏｌ−ＣＯ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＲＱ    ：呼吸商           ［−］
ＯＵＲ  ：酸素摂取速度     ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
【００４７】
酸素摂取速度ＯＵＲとＣＯ２生成速度ＣＥＲが定まると、培養槽に導入する通気ガスに含まれる酸素と二酸化炭素の量との物質収支より、培養槽内気相に含まれる酸素濃度が定まる。
【００４８】
以上をまとめると、溶存酸素濃度ＤＯが決定されると数式（２）、数式（３）より酸素摂取速度ＯＵＲが決定される。更に、酸素摂取速度ＯＵＲが決定されると、数式（４）より二酸化炭素生成速度ＣＥＲが決定される。また、酸素摂取速度ＯＵＲとＣＯ２生成速度ＣＥＲは、培養槽内気相酸素濃度にも影響を与える。培養槽内気相酸素濃度が決まると飽和溶存酸素濃度ＤＯ＊が決まり、数式（１）より飽和溶存酸素濃度ＤＯ*との関係で酸素移動速度ＯＴＲが決定される。
【００４９】
すなわち、溶存酸素濃度ＤＯは、数式（１）より酸素移動速度ＯＴＲの決定に直接関与しているが、同時に、酸素摂取速度ＯＵＲを介して間接的にも酸素移動速度ＯＴＲの決定に関与している。
【００５０】
よって、酸素移動速度ＯＴＲを予測してスケールアップを成功させるためには、飽和溶存酸素濃度ＤＯ*に加え、溶存酸素濃度ＤＯを事前に予測する手段が望まれる。
【００５１】
この問題を解決するには、溶存酸素濃度ＤＯによって定まる複数の因子に関して方程式を立て、方程式を解くことが可能になるような何らかの条件を与えた上で、解としての溶存酸素濃度ＤＯの値を得ればよい。
【００５２】
具体的には、本実施例の培養槽における酸素移動速度を計算するためのシミュレータは、溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、を備えればよい。
【００５３】
培養槽内の液相における酸素の物質収支は以下の式で表される。
【００５４】
【数６】
2019129726.tif 000008
【００５５】
ここで、
ｔ：時間［ｈ］
ＤＯ    ：溶存酸素濃度    ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＯＴＲ  ：酸素移動速度    ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＵＲ  ：酸素摂取速度    ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＭＷＯ２：酸素の分子量    ［ｇ／ｍｏｌ］
【００５６】
一般の培養では、数式（５）の左辺、すなわち溶存酸素濃度ＤＯの経時変化はそれほど急激ではないため、ｄ［ＤＯ］／ｄｔ≒０とおくことが出来る。このとき、酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲの定常状態が成り立っている。酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲの定常状態仮定は、実際に、排気ガス分析法で培養槽のｋＬａを測定する際に用いられており（非特許文献５：石崎文彬（訳）（１９８８）“発酵工学の基礎”（学会出版センター）ｐ．１７６を参照）、現実的な設定である。
【００５７】
酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲという条件を与えることで、培養槽周りの物質収支に関する連立方程式の自由度が０（ゼロ）になり、溶存酸素濃度ＤＯについて解が得られるようになる。
【００５８】
要約すると、本発明は、以下の２点を計算手順として含む。
（１）培養槽内気相についての物質収支に関する方程式を解く。特に、培養槽内気相の酸素濃度を算出し、それに応じた飽和溶存酸素濃度ＤＯ＊を算出する。
（２）溶存酸素濃度ＤＯに関する方程式をたてて、酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲとなるような定常状態仮定を置くことで、当該溶存酸素濃度ＤＯを算出する。
【００５９】
本シミュレータで使用している理論について以下に説明する。
まず、ガス物性について説明する。
気体の溶解度は、以下の式で示される。
【００６０】
【数７】
2019129726.tif 000009
ここで、溶媒が水（Ｈ２Ｏ）で、気体が酸素（Ｏ２）であるとき、定数ＡＳ＝−６６．７３５３８、ＢＳ＝８７．４７５４７、ＣＳ＝２４．４５２６４である（非特許文献６：日本化学会（編）（２００４）化学便覧改訂５版  基礎編ＩＩ（丸善）ｐ．ＩＩ−１４４）。
また、蒸気圧は、以下のＡｎｔｏｉｎｅ式（アントワン式）により示される。
【００６１】
【数８】
2019129726.tif 000010
ここで、蒸気圧Ｐ、温度Ｔ［Ｋ］、アントワン定数ＡＡ、ＢＡ、ＣＡである。
水の蒸気圧を求める場合、ＡＡ＝２３．１９６４、ＢＡ＝３８１６．４４、ＣＡ＝−４６．１３である（非特許文献７：化学工学会（編）（１９９９）化学工学便覧第６版（丸善）ｐ．１９）。
【００６２】
数式（６）、数式（７）の定数に上述した（［００５９］、［００６０］で述べた）数値を代入し、さらに気体の溶解度をモル分率から［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］に換算すると、以下の数式（８）のようになる。なお、酸素の溶解度はきわめて低いため、酸素が溶解した溶液の密度は溶媒の密度に等しいものとして換算を行っている。また、数式（８）では純水に関する数値を代入しているが、培養液に関する数値を代入してもよい。
【数９】
2019129726.tif 000011
【００６３】
ここで、
ｘ２，Ｏ２  ：気相中の酸素分圧＝大気圧での飽和溶存酸素モル分率［−］
Ｔ          ：絶対温度    ［Ｋ］
ＤＯ＊ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＭＷＯ２    ：酸素の分子量［ｇ／ｍｏｌ］
ＭＷＨ２Ｏ  ：水の分子量  ［ｇ／ｍｏｌ］
Ｄ          ：溶媒の密度  ［ｇ−ｂｒｏｔｈ／Ｌ−ｂｒｏｔｈ］
ｐｇ，Ｈ２Ｏ：水の蒸気圧  ［Ｐａ］
【００６４】
＜気相部の物質収支式＞
培養槽に入る直前の通気ガスについて説明する。
培養槽に入る直前の通気ガスが水蒸気を含まない乾燥気体とすると、通気ガスの酸素流量Ｆｇ，Ｏ２，ｆｅｅｄ、ＣＯ２流量Ｆｇ，ＣＯ２，ｆｅｅｄ、窒素流量Ｆｇ，Ｎ２，ｆｅｅｄ、は、通気ガス流量Ｆｇ，ｆｅｅｄとそれぞれのガス成分濃度Ｃｇとから、以下の計算式により算出される。
【００６５】
【数１０】
2019129726.tif 000012
【００６６】
ここで、
Ｆｇ，ｆｅｅｄ      ：通気ガス流量              ［Ｌ−ｇａｓ／ｍｉｎ］
Ｆｇ，Ｏ２，ｆｅｅｄ  ：通気ガス中の酸素流量      ［Ｌ−Ｏ２ ｇａｓ／ｍｉｎ］
Ｆｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２流量    ［Ｌ−ＣＯ２ ｇａｓ／ｍｉｎ］
Ｆｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素流量      ［Ｌ−Ｎ２ ｇａｓ／ｍｉｎ］
Ｃｇ，Ｏ２，ｆｅｅｄ：通気ガス中の酸素濃度      ［％］
Ｃｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２濃度    ［％］
Ｃｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素濃度      ［％］
【００６７】
培養槽に入った直後、通気ガスに飽和水蒸気が加わった入口ガスについて説明する。入口ガスの各成分の流量、濃度、分圧は数式（１０）〜数式（１２）で示される。入口ガスの酸素分圧より、培養液底部の飽和溶存酸素濃度は数式（１３）で示される。
【００６８】
【数１１】
2019129726.tif 000013
【００６９】
【数１２】
2019129726.tif 000014
【００７０】
【数１３】
2019129726.tif 000015
【００７１】
【数１４】
2019129726.tif 000016
ここで、
Ｆｇ, ｉｎ            ：入口ガス流量          [Ｌ−ｇａｓ／ｍｉｎ]
Ｆｇ，ｆｅｅｄ        ：通気ガス流量          [Ｌ−ｇａｓ／ｍｉｎ］
Ｆｇ, Ｏ２, ｉｎ      ：入口ガス中の酸素流量    [Ｌ−Ｏ２ ｇａｓ／ｍｉｎ]  
Ｆｇ，Ｏ２，ｆｅｅｄ  ：通気ガス中の酸素流量    [Ｌ−Ｏ２ ｇａｓ／ｍｉｎ］
Ｆｇ, ＣＯ２, ｉｎ  ：入口ガス中のＣＯ２流量    [Ｌ−ＣＯ２ ｇａｓ／ｍｉｎ]  
Ｆｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２流量    [Ｌ−ＣＯ２ｇａｓ／ｍｉｎ］
Ｆｇ, Ｎ２, ｉｎ      ：入口ガス中の窒素流量    [Ｌ−Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ，Ｎ２，ｆｅｅｄ  ：通気ガス中の窒素流量    [Ｌ−Ｎ２ｇａｓ／ｍｉｎ］
Ｆｇ, Ｈ２Ｏ, ｉｎ   ：入口ガス中の水蒸気流量    [Ｌ−Ｈ２Ｏ ｇａｓ／ｍｉｎ]
Ｃｇ, Ｏ２, ｉｎ    ：入口ガス中の酸素濃度    [％]
Ｃｇ, ＣＯ２, ｉｎ   ：入口ガス中のＣＯ２濃度    [％]
Ｃｇ, Ｎ２, ｉｎ    ：入口ガス中の窒素濃度     [％]
Ｃｇ, Ｈ２Ｏ, ｉｎ   ：入口ガス中の水蒸気濃度  [％]
Ｐｂｏｔｔｏｍ      ：培養液底部の全圧        [Ｐａ]
ｐｇ, Ｏ２, ｉｎ    ：入口ガス中の酸素分圧    [Ｐａ]
ｐｇ, ＣＯ２, ｉｎ   ：入口ガス中のＣＯ２分圧   [Ｐａ]
ｐｇ, Ｎ２, ｉｎ  ：入口ガス中の窒素分圧    [Ｐａ]
ｐｇ, Ｈ２Ｏ, ｉｎ  ：入口ガス中の水蒸気分圧    [Ｐａ]
ｐｇ，Ｈ２Ｏ  ：水の蒸気圧                [Ｐａ］
ＤＯ＊ａｔｍ  ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯ*ｂｏｔｔｏｍ  ：培養液底部の飽和溶存酸素濃度  [ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ]
Ｐａｔｍ        ：大気圧                        [Ｐａ]
【００７２】
ガスに関する物質移動について説明する。
酸素移動速度ＯＴＲは、例えば容量が５０Ｌ以下の小型培養槽では数式（１）で表される（再掲）。大型培養槽では、培養液全体の酸素移動速度ＯＴＲは数式（１）を修正して、数式（１４）で示される。また、大型培養槽の培養液の最表面と最底部では、数式（１）について、培養液の最表面と最底部それぞれの飽和溶存酸素濃度ＤＯ*と溶存酸素濃度ＤＯを適用することで、数式（１５）で示される。
【００７３】
【数１５】
2019129726.tif 000017
【数１６】
2019129726.tif 000018
【００７４】
【数１７】
2019129726.tif 000019
【００７５】
ここで、
ｋＬａ  ：酸素移動容量係数                  ［ｈ−１］
ＯＴＲ  ：培養液全体の酸素移動速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＴＲｂｏｔｔｏｍ：培養液底部の酸素移動速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＴＲｔｏｐ  ：培養液表面の酸素移動速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＤＯ*ｂｏｔｔｏｍ：培養液底部の飽和溶存酸素濃度［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯｂｏｔｔｏｍ：培養液底部の溶存酸素濃度      ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯ＊ｔｏｐ   ：培養液表面の飽和溶存酸素濃度  ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯｔｏｐ    ：培養液表面の溶存酸素濃度      ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＭＷＯ２    ：酸素の分子量                  ［ｇ／ｍｏｌ］
【００７６】
排ガスについて説明する。排ガスとは、入口ガスから消費された酸素Ｏ２を引き、発生した二酸化酸素ＣＯ２を足したものである。なお、排ガスは水蒸気で飽和しているものとする。
【００７７】
【数１８】
2019129726.tif 000020
【００７８】
【数１９】
2019129726.tif 000021
【数２０】
2019129726.tif 000022
【００７９】
【数２１】
2019129726.tif 000023
【００８０】
【数２２】
2019129726.tif 000024
【００８１】
ここで、
Ｆｇ, ｏｕｔ          ：排ガス流量                [Ｌ−ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｏｕｔ    ：排ガス中の酸素流量        [Ｌ−Ｏ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｉｎ      ：入口ガス中の酸素流量      [Ｌ−Ｏ２ ｇａｓ／ｍｉｎ]
ＯＴＲ          ：培養液全体の酸素移動速度［ｍｍｏｌ−Ｏ2／（Ｌ−ｂｒｏｔｈ・ｈ）］
Ｖｌ                  ：培養液量                  [Ｌ−ｂｒｏｔｈ]
Ｆｇ, ＣＯ２, ｏｕｔ  ：排ガス中のＣＯ２流量    [Ｌ−ＣＯ２ ｇａｓ／ｍｉｎ]
Ｆｇ, ＣＯ２, ｉｎ    ：入口ガス中のＣＯ２流量    [Ｌ−ＣＯ２ ｇａｓ／ｍｉｎ]
ＣＥＲ          ：培養液全体のＣＯ２生成速度［ｍｍｏｌ−ＣＯ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
Ｆｇ, Ｎ２, ｏｕｔ    ：排ガス中の窒素流量      [Ｌ−Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｎ２, ｉｎ      ：入口ガス中の窒素流量    [Ｌ−Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｈ２Ｏ, ｏｕｔ  ：排ガス中の水蒸気流量        [Ｌ−Ｈ２Ｏ  ｇａｓ／ｍｉｎ]
Ｃｇ, Ｏ２, ｏｕｔ    ：排ガス中の酸素濃度    [％]
Ｃｇ, ＣＯ２, ｏｕｔ  ：排ガス中のＣＯ２濃度      [％]
Ｃｇ, Ｎ２, ｏｕｔ    ：排ガス中の窒素濃度    [％]
Ｃｇ, Ｈ２Ｏ, ｏｕｔ  ：排ガス中の水蒸気濃度   [％]
Ｐｔｏｐ              ：培養液表面の全圧      [Ｐａ]
ｐｇ, Ｏ２, ｏｕｔ    ：排ガス中の酸素分圧    [Ｐａ]
ｐｇ, ＣＯ２, ｏｕｔ  ：排ガス中のＣＯ２分圧    [Ｐａ]
ｐｇ, Ｎ２, ｏｕｔ    ：排ガス中の窒素分圧    [Ｐａ]
ｐｇ, Ｈ２Ｏ, ｏｕｔ  ：排ガス中の水蒸気分圧   [Ｐａ]
ｐｇ，Ｈ２Ｏ        ：水の蒸気圧            [Ｐａ］
ＤＯ＊ｔｏｐ          ：培養液表面の飽和溶存酸素濃度  ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＤＯ＊ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
Ｐａｔｍ              ：大気圧                [Ｐａ]
【００８２】
細胞のガス代謝モデルについて説明する。細胞の酸素消費が酸素濃度に関するＭｏｎｏｄ（モノー）式に従うと、微生物による酸素比摂取速度ＱＯ２は数式（２）（再掲）で示される。大型培養槽では培養液の最表面と最底部で溶存酸素濃度ＤＯが異なることを考慮すると、以下の式で示される。
【００８３】
【数２３】
2019129726.tif 000025
【００８４】
【数２４】
2019129726.tif 000026
ここで、
ＱＯ２,ｍａｘ：最大比酸素摂取速度  ［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
Ｋｍ，Ｏ２    ：飽和定数            ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＱＯ２,ｂｏｔｔｏｍ：培養液底部の酸素比摂取速度［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
ＤＯｂｏｔｔｏｍ  ：培養液底部の溶存酸素濃度  ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
ＱＯ２，ｔｏｍ  ：培養液表面の酸素比摂取速度［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
ＤＯｔｏｐ      ：培養液表面の溶存酸素濃度  ［ｇ−Ｏ２／Ｌ−ｂｒｏｔｈ］
【００８５】
微生物の酸素摂取速度ＯＵＲは、酸素比摂取速度ＱＯ２を用いて、数式（３）（再掲）で示されるが、数式（２１）で示したように、大型培養槽では培養液の最表面と最底部でＱＯ２が異なることを考慮すると、数式（２２）で示される。
【００８６】
【数２５】
2019129726.tif 000027
【００８７】
【数２６】
2019129726.tif 000028
【００８８】
ここで、
ＯＵＲｂｏｔｔｏｍ：培養液底部の酸素摂取速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＱＯ２,ｂｏｔｔｏｍ：培養液底部の酸素比摂取速度［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
ＯＵＲｔｏｐ：培養液表面の酸素摂取速度    ［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＱＯ２, ｔｏｐ：培養液表面の酸素比摂取速度［ｇ−Ｏ２／（ｇ−ｄｒｙ  ｃｅｌｌ・ｈ）］
Ｘ            ：細胞濃度                ［ｇ−ｄｒｙ  ｃｅｌｌ／Ｌ−ｂｒｏｔｈ］
ＭＷＯ２      ：酸素の分子量            ［ｇ／ｍｏｌ］
【００８９】
また、ＣＯ２生成速度ＣＥＲ（Ｃａｒｂｏｎ  ｄｉｏｘｉｄｅ  Ｅｖｏｌｕｔｉｏｎ  Ｒａｔｅ）は、培養液の最表面と最底部それぞれのＯＵＲではなく、培養槽全体での酸素摂取速度ＯＵＲを用いれば、数式（４）で示される（再掲）。
【数２７】
2019129726.tif 000029
【００９０】
次に、定常状態仮定について説明する。本実施例における定常状態仮定とは、一定の過渡状態（後述）が発生した後に、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲとが等しくなるという定常状態を仮定することである。
【００９１】
【数２８】
2019129726.tif 000030
また、本実施例のシミュレータにおいては、大型培養槽の場合、培養液底部と培養液表面とでは酸素移動速度ＯＴＲ、酸素摂取速度ＯＵＲの値に差異が生じることに着目して、以下の計算式を用いる。
【００９２】
【数２９】
2019129726.tif 000031
ここで、
ＯＴＲｂｏｔｔｏｍ：培養液底部の酸素移動速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＵＲｂｏｔｔｏｍ：培養液底部の酸素摂取速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＴＲｔｏｐ    ：培養液表面の酸素移動速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
ＯＵＲｔｏｐ    ：培養液表面の酸素摂取速度［ｍｍｏｌ−Ｏ２／（Ｌ−ｂｒｏｔｈ・ｈ）］
【００９３】
定常状態を仮定できる理由を、概念的に説明する。
【００９４】
まず、供給過剰の例について説明する。供給過剰とは、物理的な酸素溶解速度である酸素移動速度ＯＴＲが、細胞による酸素消費速度である酸素摂取速度ＯＵＲよりも多くなった過渡的な場合である（酸素移動速度ＯＴＲ＞酸素摂取速度ＯＵＲ）。この場合、培養液内においては、溶存酸素濃度ＤＯが増加していく。溶存酸素濃度ＤＯが増加すると、飽和溶存酸素濃度ＤＯ＊と溶存酸素濃度ＤＯとの差が小さくなる。当該差が小さくなると、数式（１）（再掲）により、酸素移動速度ＯＴＲが小さくなる。すなわち、当初多かった酸素移動速度ＯＴＲが減少して、酸素摂取速度ＯＵＲに等しくなるように作用する。
【００９５】
【数３０】
2019129726.tif 000032
【００９６】
次に、消費過剰の例について説明する。消費過剰とは、酸素移動速度ＯＴＲが、細胞による酸素摂取速度ＯＵＲよりも少ない場合である（酸素移動速度ＯＴＲ＜酸素摂取速度ＯＵＲ）。この場合、培養液内においては、溶存酸素濃度ＤＯが減少していく。溶存酸素濃度ＤＯが減少すると、数式（１）により、飽和溶存酸素濃度ＤＯ＊と溶存酸素濃度ＤＯとの差が大きくなる。当該差が大きくなると、数式（１）により、酸素移動速度ＯＴＲが大きくなる。すなわち、当初少なかった酸素移動速度ＯＴＲが増加して、酸素摂取速度ＯＵＲに等しくなるように作用する。
よって、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲの値に差異があっても、系は定常状態へ向かう。
【００９７】
本実施例のシミュレータでは、上述した一連の式（数式１〜数式２４）を連立非線形方程式として収束計算（繰り返し計算）で解く。求解には各種の数値計算ソフトウェアが利用できる。例えば、ＥＱＵＡＴＲＡＮ（商標）−Ｇ（イコートラン；オメガシミュレーション社製）を用いてもよい。なお、発明者の知る限りこの様な考え方を培養槽シミュレーションに適用した例はない。
【００９８】
本実施例のシミュレータによると、排ガスの酸素濃度に応じた飽和溶存酸素濃度ＤＯ＊を算出していることで、従来よりも正確に酸素移動速度ＯＴＲを推算できる。
【００９９】
特に、酸素移動容量係数ｋＬａが同じであっても通気量が異なる場合に酸素移動速度ＯＴＲが変化する状況を予測可能である。
【０１００】
また、従来、酸素移動容量係数ｋＬａが不足してスケールアップ不能と判断されてきた条件でも、酸素移動速度ＯＴＲが十分なのであればスケールアップできる可能性について検討できる。
【０１０１】
そのほか、細胞の呼吸特性の違いによる酸素移動速度ＯＴＲの変化を計算できる。
【実施例１】
【０１０２】
図３に、本実施例によるＯＴＲシミュレータ１０００の機能ブロック図を示す。なお、本実施例の機能ブロックは、ハードウェアでもソフトウェアでもハードウェア及びソフトウェアの組み合わせでも実現可能である。
【０１０３】
ＯＴＲシミュレータ１０００は、通気量設定手段１０１０、Ｏ２濃度計算手段１０２０と、
ＯＴＲ計算手段１０３０と、データ調整手段１０４０と、ＯＵＲ計算手段１０５０と、ＤＯ＊計算手段１０６０と、ＤＯ計算手段１０７０と、操作条件ＤＢ１１１０と、入力変数ＤＢ１１２０と、出力変数ＤＢ１１３０と、制御手段１２００とから構成される。
【０１０４】
通気量設定手段１０１０は、通気量を設定する手段である。
【０１０５】
Ｏ２濃度計算手段１０２０は、他の手段やデータベース（ＤＢ：ＤａｔａＢａｓｅ）を用いて、酸素濃度を計算する手段である。
【０１０６】
ＯＴＲ計算手段１０３０は、他の手段やＤＢを用いて、酸素移動速度ＯＴＲを計算する手段である。
【０１０７】
データ調整手段１０４０は、他の手段やＤＢで使用されているデータ（値）を調整する手段であり、本実施例においては、ＯＴＲ計算手段で計算された酸素移動速度ＯＴＲの値と、ＯＵＲ計算手段で計算された酸素摂取速度ＯＵＲの値とが定常状態仮定になるようにする手段である。一実施例によると酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲの計算式が成立するように構成されてもよい。
【０１０８】
ＯＵＲ計算手段１０５０は、他の手段やＤＢを用いて、酸素摂取速度ＯＵＲを計算する手段である。
【０１０９】
ＤＯ＊計算手段１０６０は、他の手段やＤＢを用いて、飽和溶存酸素濃度ＤＯ＊を計算する手段である。
【０１１０】
ＤＯ計算手段１０７０は、他の手段やＤＢを用いて、溶存酸素濃度ＤＯを計算する手段である。
【０１１１】
装置操作条件ＤＢ１１１０は、収束計算をするにあたり必要となる計算式（方程式）に適用するための各パラメータを記憶するＤＢである。
【０１１２】
必要なパラメータは、例えば、以下のようなものである。
（１） 培養槽に関する諸元
・培養液量、培養槽深さ
・通気量
・温度、培養液表面圧力
・通気ガスの酸素濃度
・酸素移動容量係数（ｋＬａ）
（２）細胞に関する諸元
・細胞濃度
・細胞の酸素消費・ＣＯ２生成パラメータ
−  酸素に対する飽和定数
−  最大比酸素摂取速度
−  呼吸商
【０１１３】
入力変数ＤＢ１１２０は、収束計算をするにあたり必要となる入力条件を記憶するデータベース（ＤＢ）である。
【０１１４】
出力変数ＤＢ１１３０は、収束計算によって算出された値を格納するデータベース（ＤＢ）である。
【０１１５】
制御手段１２００は、ＯＴＲシミュレータ内部の各手段や各ＤＢを制御したり、ＯＴＲシミュレータ外部と信号を送受信したりする。
また、ＯＴＲシミュレータ１０００は、ユーザ端末１３００に接続されている。ユーザは、ユーザ端末１３００を介して、ＯＴＲシミュレータを作動させるために必要な変数を入力する。入力された変数は、入力変数ＤＢ１１２０に記憶される。また、ユーザ端末１３００は、ＯＴＲシミュレータを作動させた結果として得られた出力を表示する。この出力は、ＯＴＲシミュレータ１０００の出力変数ＤＢ１１３０に記憶された値である。
【０１１６】
酸素移動速度ＯＴＲの計算方法の一例を図４のフローチャートを用いて説明する。
【０１１７】
Ｓ３０１０では、排ガス酸素濃度を計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【０１１８】
【数３１】
2019129726.tif 000033
Ｓ３０２０では、飽和溶存酸素濃度ＤＯ＊を計算する。特に、上述した連立方程式において、以下の計算式について計算する。
【０１１９】
【数３２】
2019129726.tif 000034
【０１２０】
Ｓ３０３０では、溶存酸素濃度ＤＯを計算する。上述した連立方程式において、例えば、以下の式を満たす溶存酸素濃度ＤＯを計算する。
【０１２１】
【数３３】
2019129726.tif 000035
【０１２２】
Ｓ３０４０では、酸素摂取速度ＯＵＲを計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【０１２３】
【数３４】
2019129726.tif 000036
【０１２４】
Ｓ３０５０では、定常状態を仮定して酸素移動速度ＯＴＲを計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【０１２５】
【数３５】
2019129726.tif 000037
【０１２６】
Ｓ３０６０では、Ｓ３０１０−Ｓ３０６０で計算された各値が収束しているか否かを判定する。もし、収束していれば、Ｓ３０８０へ進み、収束していなければ、Ｓ３０１０へ戻り、計算を続ける。
【０１２７】
Ｓ３０８０では、収束した値を表示する。特に、酸素移動速度ＯＴＲの計算結果を表示する。
【０１２８】
なお、上述した連立方程式を計算することによって、様々な数値が算出されるが、本実施例において、一番重要な値は、酸素移動速度ＯＴＲの値である。
【０１２９】
なお、上述したフローチャートでは、Ｓ３０１０−Ｓ３０８０が逐次的に計算されるように説明したが、各変数は相互に依存しているために、上述した処理を並列に実行するように構成されている。
【０１３０】
本実施例のシミュレータは、以下のような点を考慮している。
・気相に関する物質収支から、気相部（培養排気）の酸素分圧を求める。
・細胞による酸素消費およびＣＯ２生成による培養排気の酸素分圧低下を考慮したうえで、ＯＴＲを計算する。
・（主に大型培養槽で行われる）加圧の影響を考慮する。
・液深による静水圧を考慮する。
・水蒸気分圧を考慮する。
・培養液の最表面と最深部の溶存酸素濃度分布を考慮する。
【０１３１】
また、本実施例のシミュレータは、以下のような点を前提とする。
・定常状態シミュレータであり、系の時間発展は計算しない。
・ガスホールドアップによる培養液高さの変化は考慮していない。
・細胞濃度は一定としている。
・細胞濃度の培養液内での分布は考慮しない（培養液の最表面と最底部の細胞濃度を同一と仮定）。
・ｋＬａの培養液内での分布は考慮しない（培養液の最表面と最底部のｋＬａを同一と仮定）。
・培養槽のいたるところで定常状態（酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲ）が成立していることを仮定している。
・通気は培養槽の底部から行うものとする。
・通気は乾燥気体（水蒸気を含まない）とする。
・気相部は培養槽に入ると直ちに水蒸気が飽和するものとする。
・細胞によるＣＯ２消費は考慮しない。
【実施例２】
【０１３２】
図５は、本発明の第二実施例によるＯＴＲシミュレータの機能ブロック図を示す。第１実施例（図３）に示した実施例と比較すると、ユーザ端末１３００が、撹拌装置、センサ等に置き換わっている。第１実施例では、ユーザ端末１３００で所定の値を入力していたが、第２実施例では、センサが検知した値を使って、シミュレートすることができる。別の実施例として、センサが検知した値とユーザが入力する値を任意に選択できるように構成されてもよい。
【０１３３】
本実施例のシミュレータの活用例を示す。
培養槽のサイズの影響について、本実施例のシミュレータでの計算結果を図８に示す。ここで、培養槽は、ミニジャー（小型槽）と５０ｍ３タンク（大型槽）との２つについてシミュレーションをした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］（小型槽）ｏｒ３０［ｍ３］（大型槽）
・培養液高さ＝０．１２（小型槽）ｏｒ４．５４［ｍ］（大型槽）
・温度＝３０［℃］
・圧力＝大気圧（小型槽）ｏｒ０．０３ＭＰａＧ（大型槽）
・通気気体＝大気（酸素濃度２０．９％）
・通気量＝１［ｖｖｍ］
同じ酸素移動容量係数ｋＬａであっても、５０ｍ３タンクの方が酸素移動速度ＯＴＲが大きくなることが、シミュレーション結果からわかる。
【０１３４】
本実施例のシミュレータの活用例を示す。
同じ酸素移動容量係数ｋＬａであっても、通気量が異なる場合の影響について、本実施例のシミュレータでの計算結果を図９に示す。ここで、培養槽はミニジャー（小型槽）とし、通気量は０．１、０．５、１．０、２．０ｖｖｍとした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］
・培養液高さ＝０．１２［ｍ］
・温度＝３０［℃］
・圧力＝大気圧
・通気気体＝大気（酸素濃度２０．９％）
同じ酸素移動容量係数ｋＬａであっても、通気量が低下するのに伴って酸素移動速度ＯＴＲが低下することが、シミュレーション結果からわかる。
【０１３５】
図７は、本実施例のＯＴＲシミュレータの活用例のイメージ（スケールアップ）を示す。
（１）ラボスケールの実験で、ＯＴＲ、剪断力と生産性の関係についてデータを収集。
（２）（１）のデータより等高線マップを作成。
（３）ＯＴＲシミュレータで、実機の操作可能領域を推算。
（４）（２）と（３）を重ね、生産性を最大化する操作条件と、そのとき期待される生産性を予測。
ここで、図６の出力変数の中で、ＯＴＲｔｏｔａｌが、培養槽の性能評価に用いられ、それ以外の変数は、参考値である。なお、ＯＴＲｔｏｔａｌは、培養槽全体のＯＴＲを示しており、単に、ＯＴＲと表記してもよい。
なお、図６の装置操作条件の入力変数である培養槽底部圧力Ｐｂｏｔｔｏｍは、培養槽直径を入力して簡易計算するように構成してもよい。すなわち、培養槽を円筒形として培養液量から液深を求め、液深・培養槽密度・培養液上部圧力から培養液底部の静水圧を求めるように構成してもよい。また、パラメータの入力変数（培養槽の装置パラメータ）である酸素移動容量係数ｋＬａは、攪拌速度、通気流量、装置形状、培養液物性等の関数であるが、ユーザが数値を入力するように構成してもよい。
【０１３６】
以上のように本発明の実施態様について説明したが、上述の説明に基づいて当業者にとって種々の代替例、修正又は変形が可能であり、本発明はその趣旨を逸脱しない範囲で前述の種々の代替例、修正又は変形を包含するものである。
【図１】
2019129726.tif 000038
【図２】
2019129726.tif 000039
【図３】
2019129726.tif 000040
【図４】
2019129726.tif 000041
【図５】
2019129726.tif 000042
【図６】
2019129726.tif 000043
【図７】
2019129726.tif 000044
【図８】
2019129726.tif 000045
【図９】
2019129726.tif 000046