(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-06-10
(45)【発行日】2022-06-20
(54)【発明の名称】ＯＴＲシミュレータ
(51)【国際特許分類】
C12M 1/00 20060101AFI20220613BHJP
C12N 1/00 20060101ALI20220613BHJP
C12M 1/34 20060101ALI20220613BHJP
C12M 3/02 20060101ALN20220613BHJP
【ＦＩ】
C12M1/00 D
C12N1/00 D
C12M1/34 Z
C12M3/02
【請求項の数】
8
(21)【出願番号】P 2018012709
(22)【出願日】2018-01-29
(65)【公開番号】P2019129726
(43)【公開日】2019-08-08
【審査請求日】2021-01-06
(73)【特許権者】
【識別番号】000214272
【氏名又は名称】長瀬産業株式会社
(74)【代理人】
【識別番号】110000855
【氏名又は名称】特許業務法人浅村特許事務所
(72)【発明者】
【氏名】石井 伸佳
【審査官】藤澤 雅樹
(56)【参考文献】
【文献】特開２０１０－１２４７２２（ＪＰ，Ａ）
【文献】特開２００６－２９６４２３（ＪＰ，Ａ）
【文献】特開２０１７－２０５１４７（ＪＰ，Ａ）
【文献】特開平０２－２１１８６６（ＪＰ，Ａ）
【文献】特開昭５４－０９６８６８（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１０／０１２００８２（ＵＳ，Ａ１）
【文献】米国特許出願公開第２０１７／０２５４８０１（ＵＳ，Ａ１）
【文献】国際公開第２０１５／０２５７２４（ＷＯ，Ａ１）
【文献】中国特許出願公開第１０６７７０５０９（ＣＮ，Ａ）
【文献】Biotechnology Advances (2009) Vol.27，pp.153-176
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｃ１２Ｍ １／００－３／１０
Ｃ１２Ｎ １／００－７／０８
ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）
ＣＡｐｌｕｓ／ＭＥＤＬＩＮＥ／ＥＭＢＡＳＥ／ＢＩＯＳＩＳ（ＳＴＮ）
(57)【特許請求の範囲】
【請求項1】
培養槽における酸素移動速度ＯＴＲを計算するためのシミュレータであって、（１）培養液における気相から液相への酸素移動速度を、酸素移動容量係数と培養液の飽和溶存酸素濃度と培養液の溶存酸素濃度の関数として算出する式と、
（２）単位細胞量当たりの細胞の酸素摂取速度である比酸素摂取速度を、培養液の溶存酸素濃度の関数として算出する式と、
（３）培養液全体における細胞の酸素摂取速度を、比酸素摂取速度と細胞濃度の積で算出する式と、
（４）培養液全体における細胞の二酸化炭素生産速度を、酸素摂取速度と呼吸商の積で算出する式と、
（５）培養液の溶存酸素濃度の変化速度を、酸素移動速度と酸素摂取速度の差で算出する式と、
（６）培養液の気体溶解度を培養液の温度・圧力・組成の関数として算出する式と、
（７）蒸気圧を培養液の温度・圧力・組成の関数として算出する式と、
（８）培養液の酸素溶解度を培養液の温度・圧力・組成の関数として算出する式と、その式を用いた飽和溶存酸素濃度の単位を変換する式と、水蒸気について培養液の温度・圧力・組成の関数として算出する式と、
（９）培養液最深部より通気される酸素、二酸化炭素、不活性気体の各種成分それぞれの標準体積流量を、培養液に通気される気体全体の標準体積流量と、培養液に通気される気体内の各種成分それぞれの濃度との積より算出する式と、
（１０）培養液最深部より通気された気体に水蒸気が加わった入口ガス流量を算出する式と、
（１１）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度を表す式と、
（１２）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度と、培養液最深部の全圧より、培養液最深部における各種気相部成分の分圧を求める式と、
（１３）培養液最深部の酸素分圧より、培養液最深部の飽和溶存酸素濃度を求める式と、
（１４）培養液全体の酸素移動速度を、培養液最深部と培養液最浅部それぞれの飽和溶存酸素濃度と、培養液最深部と培養液最浅部それぞれの溶存酸素濃度との関数として算出する式と、
（１５―１）培養液最深部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最深部の飽和溶存酸素濃度と培養液最深部の溶存酸素濃度の関数として算出する式と、
（１５－２）培養液最浅部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最浅部の飽和溶存酸素濃度と培養液最浅部の溶存酸素濃度の関数として算出する式と、
（１６）培養液外に排出される気体全体の標準体積流量を、酸素、二酸化炭素、不活性気体、水蒸気の各種成分の標準体積流量の和として算出する式と、
（１７－１）培養液外に排出される酸素の標準体積流量を、培養液に通気された酸素の標準体積流量と培養液全体の酸素移動速度から算出する式と、
（１７－２）培養液外に排出される二酸化炭素の標準体積流量を、培養液に通気された二酸化炭素の標準体積流量と培養液全体の二酸化炭素生産速度から算出する式と、
（１７－３）培養液外に排出される不活性気体の標準体積流量が、培養液に通気された不活性気体の標準体積流量と等しいとする式と、
（１７－４）培養液外に排出される水蒸気の標準体積流量を、培養液最浅部における水蒸気分圧と全圧、培養液外に排出される気体全体の標準体積流量より算出する式と、
（１８）培養液外に排出される気体中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部濃度を、各種成分それぞれの標準体積流量と、培養液外に排出される気体全体の体積流量より算出する式と、
（１９）培養液最浅部における気相中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部分圧を、培養液最浅部における全圧と、各種成分それぞれの気相部濃度の積より算出する式と、
（２０）培養液最浅部の酸素分圧より、培養液最浅部の飽和溶存酸素濃度を求める式と、
（２１－１）培養液最深部における細胞の比酸素摂取速度を、培養液最深部における溶存酸素濃度の関数として算出する式と、
（２１－２）培養液最浅部における細胞の比酸素摂取速度を、培養液最浅部における溶存酸素濃度の関数として算出する式と、
（２２－１）培養液最深部における細胞の酸素摂取速度を、培養液最深部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２２－２）培養液最浅部における細胞の酸素摂取速度を、培養液最浅部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２３）培養液全体における酸素移動速度が、培養液全体における細胞の酸素摂取速度と等しいとする式と、
（２４－１）培養液最深部における細胞の酸素摂取速度が培養液最深部における酸素移動速度と等しいとする式と、
（２４－２）培養液最浅部における細胞の酸素摂取速度が培養液最浅部における酸素移動速度と等しいとする式と、
を連立非線形方程式として収束計算する機能を有する手段を備えるシミューレータ。
【請求項2】
少なくとも酸素移動速度が収束したことを判定する手段を更に備えたことを特徴とする、請求項１に記載のシミュレータ。
【請求項3】
（１）培養液における気相から液相への酸素移動速度（ＯＴＲ）を、酸素移動容量係数（ｋＬａ）と培養液の飽和溶存酸素濃度（ＤＯ＊）と培養液の溶存酸素濃度（ＤＯ）の関数として算出する式は、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（５）培養液の溶存酸素濃度（ＤＯ）の変化速度（ｄ［ＤＯ］／ｄｔ）を、酸素移動速度（ＯＴＲ）と酸素摂取速度（ＯＵＲ）の差で算出する式は、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（６）培養液の気体溶解度（ｌｎ［ｘ２］）を培養液の温度・圧力・組成の関数として算出する式は、
により計算され、
（７）水蒸気圧（Ｐ）を培養液の温度・圧力・組成の関数として算出する式は、定数ＡＡ＝－６６．７３５３８、定数ＢＡ＝８７．４７５４７、定数ＣＡ＝２４．４５２６４を用いて、
により計算され、
（８）培養液の酸素溶解度を培養液の温度・圧力・組成の関数として算出する式と、その式を用いた飽和溶存酸素濃度の単位を変換する式と、水蒸気について培養液の温度・圧力・組成の関数として算出する式は、培養液の密度（D）を用いて、
により計算され、
（１１）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度（Ｃｇ,Ｏ２,ｉｎ、Ｃｇ,ＣＯ２,ｉｎ、Ｃｇ,Ｎ２,ｉｎ、Ｃｇ,Ｈ２Ｏ,ｉｎ）を表す式は、入口ガス流量（Ｆｇ, ｉｎ）、入口ガス中の酸素流量（Ｆｇ, Ｏ２, ｉｎ）、入口ガス中のＣＯ２流量（Ｆｇ, ＣＯ２, ｉｎ）、入口ガス中の窒素流量（Ｆｇ, Ｎ２, ｉｎ）、入口ガス中の水蒸気流量（Ｆｇ, Ｈ２Ｏ, ｉｎ ）を用いて、
により計算され、
（１３）培養液最深部の酸素分圧（ＤＯ＊ ａｔｍ）より、培養液最深部の飽和溶存酸素濃度（ＤＯ* ｂｏｔｔｏｍ）を求める式は、酸素の気相部成分の分圧（ｐｇ,Ｏ２,ｉｎ）を用いて、
により計算され、
（１４）培養液全体の酸素移動速度（ＯＴＲ）を、培養液最深部と培養液最浅部それぞれの飽和溶存酸素濃度（ＤＯ* ｂｏｔｔｏｍ、ＤＯｂｏｔｔｏｍ）と、培養液最深部と培養液最浅部それぞれの溶存酸素濃度（ＤＯ＊ ｔｏｐ、ＤＯｔｏｐ）との関数として算出する式は、酸素移動容量係数（ｋＬａ）と、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（１５）培養液最深部における気相から液相への酸素移動速度（ＯＴＲｂｏｔｔｏｍ）を、酸素移動容量係数（ｋＬａ）と培養液最深部の飽和溶存酸素濃度（ＤＯ* ｂｏｔｔｏｍ）と培養液最深部の溶存酸素濃度（ＤＯｂｏｔｔｏｍ）の関数として算出する式と、
培養液最浅部における気相から液相への酸素移動速度（ＯＴＲｔｏｐ）を、酸素移動容量係数（ｋＬａ）と培養液最浅部の飽和溶存酸素濃度（ＤＯ＊ ｔｏｐ）と培養液最浅部の溶存酸素濃度（ＤＯｔｏｐ）の関数として算出する式は、酸素移動容量係数（ｋＬａ）と、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（１７）培養液外に排出される酸素の標準体積流量（Ｆｇ,Ｏ２,ｏｕｔ）を、培養液に通気された酸素の標準体積流量（Ｆｇ,Ｏ２,ｉｎ）と培養液全体の酸素移動速度（ＯＴＲ）から算出する式と、
培養液外に排出される二酸化炭素の標準体積流量（Ｆｇ,ＣＯ２,ｏｕｔ）を、培養液に通気された二酸化炭素の標準体積流量（Ｆｇ,ＣＯ２,ｉｎ）と培養液全体の二酸化炭素生産速度（ＣＥＲ）から算出する式と、
培養液外に排出される不活性気体の標準体積流量（Ｆｇ,Ｎ２,ｏｕｔ）が、培養液に通気された不活性気体の標準体積流量（Ｆｇ,Ｎ２,ｉｎ）と等しくする式と、
培養液外に排出される水蒸気の標準体積流量（Ｆｇ,Ｈ２Ｏ,ｏｕｔ）を、培養液最浅部における水蒸気分圧（ｐｇ,Ｈ２Ｏ,ｏｕｔ）と全圧（Ｐｔｏｐ）、培養液外に排出される気体全体の標準体積流量（Ｆｇ,ｏｕｔ）より算出する式は、培養液量（Ｖ１）を用いて、
により計算され、
（１８）培養液外に排出される気体中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部濃度（Ｃｇ,Ｏ２,ｏｕｔ、Ｃｇ,ＣＯ２,ｏｕｔ、Ｃｇ Ｎ２,ｏｕｔ、Ｃｇ,Ｈ２Ｏ,ｏｕｔ）を、各種成分それぞれの標準体積流量（Ｆｇ,Ｏ２,ｏｕｔ、Ｆｇ,ＣＯ２,ｏｕｔ、Ｆｇ,Ｎ２,ｏｕｔ、Ｆｇ,Ｈ２Ｏ,ｏｕｔ）と、培養液外に排出される気体全体の体積流量（Ｆｇ,ｏｕｔ）より算出する式は、
により計算され、
（２０）培養液最浅部の酸素分圧（ＤＯ＊ ａｔｍ）より、培養液最浅部の飽和溶存酸素濃度（ＤＯ＊ ｔｏｐ）を求める式は、排ガス中の酸素分圧（ｐｇ,Ｏ２,ｏｕｔ）と、大気圧（Ｐａｔｍ）を用いて、
により計算され、
（２１）培養液最深部における細胞の比酸素摂取速度（ＱＯ２,ｂｏｔｔｏｍ）を、培養液最深部における溶存酸素濃度（ＤＯｂｏｔｔｏｍ）の関数として算出する式と、
培養液最浅部における細胞の比酸素摂取速度（ＱＯ２,ｔｏｐ）を、培養液最浅部における溶存酸素濃度（ＤＯｔｏｐ）の関数として算出する式は、最大比酸素摂取速度（ＱＯ２,ｍａｘ）と、飽和定数（Ｋｍ，Ｏ２）を用いて、
により計算される、請求項１に記載のシミュレータ。
【請求項4】
（２）単位細胞量当たりの細胞の酸素摂取速度である比酸素摂取速度（ＱＯ２）を、培養液の溶存酸素濃度（ＤＯ）の関数として算出する式は、最大比酸素摂取速度（ＱＯ２,ｍａｘ）と飽和定数（Ｋｍ，Ｏ２）とを用いて、
により計算される、請求項３に記載のシミュレータ。
【請求項5】
（３）培養液全体における細胞の酸素摂取速度（ＯＵＲ）を、比酸素摂取速度（ＱＯ２）と細胞濃度の積（Ｘ）で算出する式は、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（４）培養液全体における細胞の二酸化炭素生産速度（ＣＥＲ）を、酸素摂取速度（ＯＵＲ）と呼吸商（ＲＱ）の積で算出する式は、
により計算され、
（９）培養液最深部より通気される酸素、二酸化炭素、不活性気体の各種成分それぞれの標準体積流量（Ｆｇ，Ｏ２，ｆｅｅｄ、Ｆｇ，ＣＯ２，ｆｅｅｄ、Ｆｇ，Ｎ２，ｆｅｅｄ）を、培養液に通気される気体全体の標準体積流量（Ｆｇ，ｆｅｅｄ）と、培養液に通気される気体内の各種成分それぞれの濃度（Ｃｇ，Ｏ２，ｆｅｅｄ、Ｃｇ，ＣＯ２，ｆｅｅｄ、Ｃｇ，Ｎ２，ｆｅｅｄ）との積より算出する式は、
により計算され、
（１０）培養液最深部より通気された気体（Ｆｇ，ｆｅｅｄ）に水蒸気（Ｆｇ,Ｈ２Ｏ,ｉｎ）が加わった入口ガス流量（Ｆｇ, ｉｎ）を算出する式は、
により計算され、
（１２）培養液最深部における各種気相部成分の分圧（ｐｇ,Ｏ２,ｉｎ、ｐｇ,ＣＯ２,ｉｎ、ｐｇ,Ｎ２,ｉｎ、ｐｇ,Ｈ２Ｏ,ｉｎ）は、培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度（Ｃｇ,Ｏ２,ｉｎ、Ｃｇ,ＣＯ２,ｉｎ、Ｃｇ,Ｎ２,ｉｎ、ＣｇＨ２Ｏ,ｉｎ）と、培養液最深部の全圧（Ｐｂｏｔｔｏｍ）の積より、
により計算され、
（１６）培養液外に排出される気体全体の標準体積流量（Ｆｇ,ｏｕｔ）を、酸素、二酸化炭素、不活性気体、水蒸気の各種成分の標準体積流量（Ｆｇ,Ｏ２,ｏｕｔ、Ｆｇ,ＣＯ２,ｏｕｔ、Ｆｇ,Ｎ２,ｏｕｔ、Ｆｇ,Ｈ２Ｏ,ｏｕｔ）の和として算出する式は、
により計算され、
（１９）培養液最浅部における気相中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部分圧（ｐｇ,Ｏ２,ｏｕｔ、ｐｇ,ＣＯ２,ｏｕｔ、ｐｇ,Ｎ２,ｏｕｔ、ｐｇ,Ｈ２Ｏ,ｏｕｔ）を、培養液最浅部における全圧（Ｐｔｏｐ）と、各種成分それぞれの気相部濃度（Ｃｇ,Ｏ２,ｏｕｔ、Ｃｇ,ＣＯ２,ｏｕｔ、Ｃｇ Ｎ２,ｏｕｔ、Ｃｇ,Ｈ２Ｏ,ｏｕｔ）の積より算出する式は、水の蒸気圧（ｐｇ，Ｈ２Ｏ）を用いて、
により計算され、
（２２）培養液最深部における細胞の酸素摂取速度（ＯＵＲｂｏｔｔｏｍ）を、培養液最深部における細胞の比酸素摂取速度（ＱＯ２,ｂｏｔｔｏｍ）と細胞濃度（Ｘ）の積より算出する式と、
培養液最浅部における細胞の酸素摂取速度（ＯＵＲｔｏｐ）を、培養液最浅部における細胞の比酸素摂取速度（ＱＯ２,ｔｏｐ）と細胞濃度（Ｘ）の積より算出する式は、酸素の分子量（ＭＷＯ２）を用いて、
により計算され、
（２３）培養液全体における酸素移動速度（ＯＵＲ）が、培養液全体における細胞の酸素摂取速度（ＯＵＲ）と等しくする式は、
により計算され、
（２４）培養液最深部における細胞の酸素摂取速度（ＯＵＲｂｏｔｔｏｍ）が培養液最深部における酸素移動速度（ＯＴＲｂｏｔｔｏｍ）と等しくする式と、
培養液最浅部における細胞の酸素摂取速度（ＯＵＲｔｏｐ）が培養液最浅部における酸素移動速度（ＯＴＲｔｏｐ）と等しくする式は、
により計算される、請求項４に記載のシミュレータ。
【請求項6】
前記シミュレータは、攪拌装置とセンサとを含む培養層と接続可能であるように構成されており、前記攪拌装置の翼先端速度と、酸素移動速度との関係に基づいて、等高線マップを生成することを特徴とする、請求項１から５のいずれか一項に記載のシミュレータ。
【請求項7】
培養槽における酸素移動速度を計算するためのプログラムであって、（１）培養液における気相から液相への酸素移動速度を、酸素移動容量係数と培養液の飽和溶存酸素濃度と培養液の溶存酸素濃度の関数として算出する式と、
（２）単位細胞量当たりの細胞の酸素摂取速度である比酸素摂取速度を、培養液の溶存酸素濃度の関数として算出する式と、
（３）培養液全体における細胞の酸素摂取速度を、比酸素摂取速度と細胞濃度の積で算出する式と、
（４）培養液全体における細胞の二酸化炭素生産速度を、酸素摂取速度と呼吸商の積で算出する式と、
（５）培養液の溶存酸素濃度の変化速度を、酸素移動速度と酸素摂取速度の差で算出する式と、
（６）培養液の気体溶解度を培養液の温度・圧力・組成の関数として算出する式と、
（７）蒸気圧を培養液の温度・圧力・組成の関数として算出する式と、
（８）培養液の酸素溶解度を培養液の温度・圧力・組成の関数として算出する式と、その式を用いた飽和溶存酸素濃度の単位を変換する式と、水蒸気について培養液の温度・圧力・組成の関数として算出する式と、
（９）培養液最深部より通気される酸素、二酸化炭素、不活性気体の各種成分それぞれの標準体積流量を、培養液に通気される気体全体の標準体積流量と、培養液に通気される気体内の各種成分それぞれの濃度との積より算出する式と、
（１０）培養液最深部より通気された気体に水蒸気が加わった入口ガス流量を算出する式と、
（１１）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度を表す式と、
（１２）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度と、培養液最深部の全圧より、培養液最深部における各種気相部成分の分圧を求める式と、
（１３）培養液最深部の酸素分圧より、培養液最深部の飽和溶存酸素濃度を求める式と、
（１４）培養液全体の酸素移動速度を、培養液最深部と培養液最浅部それぞれの飽和溶存酸素濃度と、培養液最深部と培養液最浅部それぞれの溶存酸素濃度との関数として算出する式と、
（１５－１）培養液最深部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最深部の飽和溶存酸素濃度と培養液最深部の溶存酸素濃度の関数として算出する式と、
（１５－２）培養液最浅部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最浅部の飽和溶存酸素濃度と培養液最浅部の溶存酸素濃度の関数として算出する式と、
（１６）培養液外に排出される気体全体の標準体積流量を、酸素、二酸化炭素、不活性気体、水蒸気の各種成分の標準体積流量の和として算出する式と、
（１７－１）培養液外に排出される酸素の標準体積流量を、培養液に通気された酸素の標準体積流量と培養液全体の酸素移動速度から算出する式と、
（１７－２）培養液外に排出される二酸化炭素の標準体積流量を、培養液に通気された二酸化炭素の標準体積流量と培養液全体の二酸化炭素生産速度から算出する式と、
（１７－３）培養液外に排出される不活性気体の標準体積流量が、培養液に通気された不活性気体の標準体積流量と等しいとする式と、
（１７－４）培養液外に排出される水蒸気の標準体積流量を、培養液最浅部における水蒸気分圧と全圧、培養液外に排出される気体全体の標準体積流量より算出する式と、
（１８）培養液外に排出される気体中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部濃度を、各種成分それぞれの標準体積流量と、培養液外に排出される気体全体の体積流量より算出する式と、
（１９）培養液最浅部における気相中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部分圧を、培養液最浅部における全圧と、各種成分それぞれの気相部濃度の積より算出する式と、
（２０）培養液最浅部の酸素分圧より、培養液最浅部の飽和溶存酸素濃度を求める式と、
（２１－１）培養液最深部における細胞の比酸素摂取速度を、培養液最深部における溶存酸素濃度の関数として算出する式と、
（２１－２）培養液最浅部における細胞の比酸素摂取速度を、培養液最浅部における溶存酸素濃度の関数として算出する式と、
（２２－１）培養液最深部における細胞の酸素摂取速度を、培養液最深部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２２－２）培養液最浅部における細胞の酸素摂取速度を、培養液最浅部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２３）培養液全体における酸素移動速度が、培養液全体における細胞の酸素摂取速度と等しいとする式と、
（２４－１）培養液最深部における細胞の酸素摂取速度が培養液最深部における酸素移動速度と等しいとする式と、
（２４－２）培養液最浅部における細胞の酸素摂取速度が培養液最浅部における酸素移動速度と等しいとする式と、
を連立非線形方程式を実行する手段としてコンピュータを機能させるためのプログラム。
【請求項8】
培養槽における酸素移動速度をコンピュータによって計算するための方法であって、（１）培養液における気相から液相への酸素移動速度を、酸素移動容量係数と培養液の飽和溶存酸素濃度と培養液の溶存酸素濃度の関数として算出する式と、
（２）単位細胞量当たりの細胞の酸素摂取速度である比酸素摂取速度を、培養液の溶存酸素濃度の関数として算出する式と、
（３）培養液全体における細胞の酸素摂取速度を、比酸素摂取速度と細胞濃度の積で算出する式と、
（４）培養液全体における細胞の二酸化炭素生産速度を、酸素摂取速度と呼吸商の積で算出する式と、
（５）培養液の溶存酸素濃度の変化速度を、酸素移動速度と酸素摂取速度の差で算出する式と、
（６）培養液の気体溶解度を培養液の温度・圧力・組成の関数として算出する式と、
（７）蒸気圧を培養液の温度・圧力・組成の関数として算出する式と、
（８）培養液の酸素溶解度を培養液の温度・圧力・組成の関数として算出する式と、その式を用いた飽和溶存酸素濃度の単位を変換する式と、水蒸気について培養液の温度・圧力・組成の関数として算出する式と、
（９）培養液最深部より通気される酸素、二酸化炭素、不活性気体の各種成分それぞれの標準体積流量を、培養液に通気される気体全体の標準体積流量と、培養液に通気される気体内の各種成分それぞれの濃度との積より算出する式と、
（１０）培養液最深部より通気された気体に水蒸気が加わった入口ガス流量を算出する式と、
（１１）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度を表す式と、
（１２）培養液最深部より通気された気体に水蒸気が加わった直後の酸素、二酸化炭素、不活性気体、水蒸気の各種成分の気相部濃度と、培養液最深部の全圧より、培養液最深部における各種気相部成分の分圧を求める式と、
（１３）培養液最深部の酸素分圧より、培養液最深部の飽和溶存酸素濃度を求める式と、
（１４）培養液全体の酸素移動速度を、培養液最深部と培養液最浅部それぞれの飽和溶存酸素濃度と、培養液最深部と培養液最浅部それぞれの溶存酸素濃度との関数として算出する式と、
（１５－１）培養液最深部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最深部の飽和溶存酸素濃度と培養液最深部の溶存酸素濃度の関数として算出する式と、
（１５－２）培養液最浅部における気相から液相への酸素移動速度を、酸素移動容量係数と培養液最浅部の飽和溶存酸素濃度と培養液最浅部の溶存酸素濃度の関数として算出する式と、
（１６）培養液外に排出される気体全体の標準体積流量を、酸素、二酸化炭素、不活性気体、水蒸気の各種成分の標準体積流量の和として算出する式と、
（１７－１）培養液外に排出される酸素の標準体積流量を、培養液に通気された酸素の標準体積流量と培養液全体の酸素移動速度から算出する式と、
（１７－２）培養液外に排出される二酸化炭素の標準体積流量を、培養液に通気された二酸化炭素の標準体積流量と培養液全体の二酸化炭素生産速度から算出する式と、
（１７－３）培養液外に排出される不活性気体の標準体積流量が、培養液に通気された不活性気体の標準体積流量と等しいとする式と、
（１７－４）培養液外に排出される水蒸気の標準体積流量を、培養液最浅部における水蒸気分圧と全圧、培養液外に排出される気体全体の標準体積流量より算出する式と、
（１８）培養液外に排出される気体中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部濃度を、各種成分それぞれの標準体積流量と、培養液外に排出される気体全体の体積流量より算出する式と、
（１９）培養液最浅部における気相中の酸素、二酸化炭素、不活性気体、水蒸気の各種成分それぞれの気相部分圧を、培養液最浅部における全圧と、各種成分それぞれの気相部濃度の積より算出する式と、
（２０）培養液最浅部の酸素分圧より、培養液最浅部の飽和溶存酸素濃度を求める式と、
（２１－１）培養液最深部における細胞の比酸素摂取速度を、培養液最深部における溶存酸素濃度の関数として算出する式と、
（２１－２）培養液最浅部における細胞の比酸素摂取速度を、培養液最浅部における溶存酸素濃度の関数として算出する式と、
（２２－１）培養液最深部における細胞の酸素摂取速度を、培養液最深部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２２－２）培養液最浅部における細胞の酸素摂取速度を、培養液最浅部における細胞の比酸素摂取速度と細胞濃度の積より算出する式と、
（２３）培養液全体における酸素移動速度が、培養液全体における細胞の酸素摂取速度と等しいとする式と、
（２４－１）培養液最深部における細胞の酸素摂取速度が培養液最深部における酸素移動速度と等しいとする式と、
（２４－２）培養液最浅部における細胞の酸素摂取速度が培養液最浅部における酸素移動速度と等しいとする式と、
を連立非線形方程式として収束計算するステップを含むことを特徴とする方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、ＯＴＲ（Ｏｘｙｇｅｎ Ｔｒａｎｓｆｅｒ Ｒａｔｅ；酸素移動速度）をシミュレートするシステム、装置、プログラムに関するものである。
【背景技術】
【0002】
製造現場の大型生産機などの培養槽内で、微生物や動植物細胞を培養するためには、培養槽内における各種パラメータを決定する必要がある。研究所などの小型の培養槽（ラボスケール）で得られた結果を基に、製造現場の大型生産機にも適用できるように、各種パラメータ（大型の培養槽の装置寸法や運転条件など）を決定することをスケールアップという。
【0003】
培養槽が完全に同じ性能になるようにスケールアップすることは難しい。これは、例え幾何学的に相似になるよう培養槽を設計しても、液混合特性や物質移動特性などの特性を異なるスケールの装置同士で完全に一致させられないからである。
【0004】
したがって、スケールアップにあたっては、何が重要な因子であるかを見抜き、少なくともその因子については小型培養槽と大型生産機で出来るだけ一致するように各種条件を設定することが重要である。
【0005】
好気的な培養では、多くの場合、培養槽内の気相から液相への酸素の移動速度が細胞量や目的物質の生産性と非常に高く相関することが知られる。したがって、特に、好気的な培養では培養槽の酸素移動速度が重要なスケールアップ因子となる。
【0006】
＜酸素移動速度ＯＴＲの計算式＞
培養槽における気相から液相への酸素移動速度は以下の式で示すことができる。
培養槽における気相から液相への酸素移動速度は以下の式で示すことができる。
【0007】
【数1】
【0008】
ＯＴＲ ：酸素移動速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ｋＬａ ：酸素移動容量係数 ［ｈ－１］
ＤＯ＊ ：飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
ｋＬａ ：酸素移動容量係数 ［ｈ－１］
ＤＯ＊ ：飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
【0009】
ここで、ＯＴＲは、Ｏｘｙｇｅｎ Ｔｒａｎｓｆｅｒ Ｒａｔｅ（酸素移動速度）を意味しており、物理的な酸素溶解速度を示す。
好気的な培養で小型培養槽の結果を大型生産機で再現するためには、両者の酸素移動速度ＯＴＲを一致させることが極めて重要である。
好気的な培養で小型培養槽の結果を大型生産機で再現するためには、両者の酸素移動速度ＯＴＲを一致させることが極めて重要である。
【0010】
＜スケールアップ＞
スケールアップの手段として、小型実験機と大型生産機との両方において、酸素移動容量係数ｋＬａを一致させることが、従来技術では行われてきた。しかしながら、数式（１）からもわかるように、酸素移動速度ＯＴＲを決める因子は、酸素移動容量係数ｋＬａの他にも、溶存酸素濃度ＤＯと、飽和溶存酸素濃度ＤＯ＊がある。これらの因子が全て一致しないと、スケールアップしたときに酸素移動速度ＯＴＲが一致しない。
スケールアップの手段として、小型実験機と大型生産機との両方において、酸素移動容量係数ｋＬａを一致させることが、従来技術では行われてきた。しかしながら、数式（１）からもわかるように、酸素移動速度ＯＴＲを決める因子は、酸素移動容量係数ｋＬａの他にも、溶存酸素濃度ＤＯと、飽和溶存酸素濃度ＤＯ＊がある。これらの因子が全て一致しないと、スケールアップしたときに酸素移動速度ＯＴＲが一致しない。
【0011】
スケールアップでよく行われるのは、更に、培養中に溶存酸素濃度をオンラインで測定し、それがラボスケールと一致するように撹拌速度等を制御する方法である。これを、溶存酸素濃度制御（ＤＯ制御）という。
【0012】
スケールアップ時に、酸素移動容量係数ｋＬａを合わせたうえでＤＯ制御を行うと、数式（１）のうち酸素移動容量係数ｋＬａ、溶存酸素濃度ＤＯが一致するので、飽和溶存酸素濃度ＤＯ＊がおおむね同じであれば酸素移動速度ＯＴＲも一致することになる。
【0013】
しかし、上記の方法では、溶存酸素濃度ＤＯを実測しながら調整する必要があり、所与の性能の培養槽において所望の溶存酸素濃度ＤＯが得られるかを事前に予測していない。つまり、実際にスケールアップを行ってみたら、所望の溶存酸素濃度ＤＯが維持できず、結局、所望の酸素移動速度ＯＴＲを実現できないという事態がありうる。
【0014】
なお、撹拌速度等を変更してＤＯ制御を行った場合、酸素移動容量係数ｋＬａも変化しているはずであり、結果としてＤＯがラボスケールと一致したとしても酸素移動速度ＯＴＲが一致する保証はない。酸素移動速度ＯＴＲの観点からは、ＤＯ制御は必ずしもスケールアップにおいて有効な手段とは限らないことに注意が必要である。
【0015】
酸素移動速度ＯＴＲは測定が可能であり（非特許文献１）、大型生産機と小型実験機の双方で各種の操作条件とＯＴＲの関係を予め測定しておくことで、小型実験機と酸素移動速度ＯＴＲが一致するような大型生産機の操作条件を探索することは技術的には可能である。
【0016】
しかしながら、実生産で現に使用中の大型生産機では、多数の実験によって酸素移動速度ＯＴＲを測定することはコストや装置使用スケジュールの面から困難であり、理論計算によって酸素移動速度ＯＴＲを精度よく予測する方法が望まれる。また、新設の大型生産機を設計する段階でも、やはり、酸素移動速度ＯＴＲを事前に予測することが不可欠である。
【0017】
以上の要望より、従来においても培養槽の酸素移動速度ＯＴＲの予測計算は行われてきた。しかしながら、酸素移動速度ＯＴＲには各種パラメータが相互に関係していて複雑な計算が要求されるために、様々な仮定をおいての簡易計算が行われてきたのが実情である。
【0018】
酸素移動速度ＯＴＲの計算のためには、数式（１）の右辺に現れる各パラメータに数値を与える必要がある。例えば、非特許文献２では、溶存酸素濃度ＤＯ＝０と仮定している。また、非特許文献３では、溶存酸素濃度ＤＯ＝ＤＯｃｒｉｔ（細胞の生産活動のために許容できる最低値）と仮定している。さらに、非特許文献２、非特許文献３のいずれも、飽和溶存酸素濃度ＤＯ＊は通気ガスである大気の酸素濃度（２１％）と培養液にかかる圧力より計算している。
【0019】
しかしながら、溶存酸素濃度ＤＯは、実際には必ずしも０（ゼロ）ではなく、一般に臨界溶存酸素濃度ＤＯｃｒｉｔよりも大きい値を維持するように培養操作が実施されることが多い。さらに、培養液として接触している気相の酸素濃度（排ガスの酸素濃度）は、細胞の呼吸活動によって酸素が消費される結果、通気ガスの酸素濃度とは一致しないので、気相の酸素濃度が通気ガスと同じ（通気ガスが大気の場合、２１％）であることを前提として、飽和溶存酸素濃度ＤＯ＊を計算することは適切ではない。従来技術では、以上のような適切ではない仮定の結果として、酸素移動速度ＯＴＲが過大に計算される。このため、実際にスケールアップを行ってみると大型生産機に適合する条件が得られない場合も少なくない。
【先行技術文献】
【非特許文献】
【0020】
【文献】田口久治／永井史郎（１９８５）“微生物培養工学”（共立出版）Ｐ．１６７
【文献】川瀬義矩（１９９３）“生物反応工学の基礎”（化学工業社）ｐ．２２４－２２５
【文献】吉田敏臣（１９９８）“培養工学”（コロナ社）ｐ．５８－ｐ．６２
【文献】小林猛、本田裕之（２００２）”生物化学工学”（東京化学同人）ｐ．６２
【文献】石崎文彬（訳）（１９８８）“発酵工学の基礎”（学会出版センター）ｐ．１７６
【文献】日本化学会（編）（２００４）化学便覧改訂５版 基礎編ＩＩ（丸善）ｐ．ＩＩ－１４４
【文献】化学工学会（編）（１９９９）化学工学便覧第６版（丸善）ｐ．１９
【発明の概要】
【発明が解決しようとする課題】
【0021】
本発明では、特に、培養液中の溶存酸素濃度ＤＯ、飽和溶存酸素濃度ＤＯ＊を事前に予測しつつ、ラボスケールからスケールアップ時に実現可能な酸素移動速度ＯＴＲを推定するシミュレータを提供することを目的のひとつとする。
【課題を解決するための手段】
【0022】
気相に関する物質収支を解き、特に、排ガス酸素濃度を算出し、それに応じた飽和溶存酸素濃度ＤＯ＊を算出する。そして、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲ（Ｏｘｙｇｅｎ Ｕｐｔａｋｅ Ｒａｔｅ）との間で定常状態（酸素移動速度と酸素摂取速度が等しい）が成立すると見做せることを仮定したシミュレータを提供する。
【0023】
具体的には、培養槽における酸素移動速度を計算するためのシミュレータは、溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、を備えればよい。
【発明の効果】
【0024】
本発明のシミュレータによれば、溶存酸素濃度ＤＯ、飽和溶存酸素濃度ＤＯ＊を推定しつつ、与えられた酸素移動容量係数ｋＬａ下での酸素移動速度ＯＴＲを予測可能となる。これにより、スケールアップ時にも使用可能な酸素移動速度ＯＴＲが算出できるようになり、スケールアップの際の培養槽の条件をどのように設定すべきかを事前に検討することが可能となる。よって、スケールアップにおける試行錯誤を大幅に削減できる。
本発明の他の目的、特徴及び利点は添付図面に関する以下の本発明の実施例の記載から明らかになるであろう。
本発明の他の目的、特徴及び利点は添付図面に関する以下の本発明の実施例の記載から明らかになるであろう。
【図面の簡単な説明】
【0025】
【図1】培養槽に酸素を供給する概念図を示す。
【図2】本発明の一実施例によるＯＴＲシミュレータの概念図を示す。
【図3】本発明の一実施例によるＯＴＲシミュレータの機能ブロック図を示す。
【図4】本発明の一実施例によるＯＴＲシミュレータのフローチャートを示す。
【図5】本発明の別の実施例によるＯＴＲシミュレータの機能ブロック図を示す。
【図6】酸素移動速度ＯＴＲをシミュレートするにあたって、各パラメータとの相関関係を示す。
【図7】本発明の一実施例によるＯＴＲシミュレータの活用例のイメージ（スケールアップ）を示す。
【図8】培養槽のサイズの影響について、本発明の一実施例によるＯＴＲシミュレータでの計算結果を示す。
【図9】同じ酸素移動容量係数ｋＬａであっても、通気量が異なる場合の影響について、本発明の一実施例によるＯＴＲシミュレータでの計算結果を示す。
【発明を実施するための形態】
【0026】
＜酸素移動速度ＯＴＲの計算式＞
培養槽における気相から液相への酸素移動速度は、数式（１）（再掲）で示したように、酸素移動容量係数ｋＬａ、飽和溶存酸素濃度ＤＯ*、溶存酸素濃度ＤＯによって定まる。
培養槽における気相から液相への酸素移動速度は、数式（１）（再掲）で示したように、酸素移動容量係数ｋＬａ、飽和溶存酸素濃度ＤＯ*、溶存酸素濃度ＤＯによって定まる。
【数2】
【0027】
ｋＬａは、酸素移動容量係数［ｈ－１］である。酸素移動速度ＯＴＲの大きさを決める定数であり、気相から液相への酸素の移動しやすさを表しており、装置形状、通気条件、撹拌条件、培養液物性によって決まると言われている。
【0028】
＜飽和溶存酸素濃度ＤＯ＊＞
ＤＯ＊は、飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ］であり、単位体積あたりの培養液内に溶存する酸素の最大限の濃度を示しており、その値はＨｅｎｒｙの法則に従って定まる。Ｈｅｎｒｙの法則とは、温度が一定のとき、一定量の液体に溶解する気体の質量はその気体の圧力（分圧）に比例するという法則である。すなわち、気相の酸素分圧が変化すると、飽和溶存酸素濃度ＤＯ＊も変化することになる。酸素の分圧は、気相の酸素濃度に依存する。
ＤＯ＊は、飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ］であり、単位体積あたりの培養液内に溶存する酸素の最大限の濃度を示しており、その値はＨｅｎｒｙの法則に従って定まる。Ｈｅｎｒｙの法則とは、温度が一定のとき、一定量の液体に溶解する気体の質量はその気体の圧力（分圧）に比例するという法則である。すなわち、気相の酸素分圧が変化すると、飽和溶存酸素濃度ＤＯ＊も変化することになる。酸素の分圧は、気相の酸素濃度に依存する。
【0029】
ここで、培養槽内の気相の酸素濃度が、単純に培養槽に導入される通気ガスの酸素濃度に等しければ、飽和溶存酸素濃度ＤＯ＊の計算は、それほど複雑にはならない。このため、従来は、飽和溶存酸素濃度ＤＯ＊は通気ガスの酸素濃度のみから定まるとして計算が行われていた。現実には、培養槽内の気相の酸素濃度は通気ガスの酸素濃度から変化するため従来法は必ずしも妥当ではない。
【0030】
培養槽内の気相の酸素分圧は、通気ガスに含まれる酸素の量と気相から培養液に移動する酸素の量、気相に培養液から排出されたガス（一般に、排出ガスの大半は二酸化炭素ＣＯ２であるが、水蒸気も加わる）の量、および気相の全圧（培養槽が加圧されている場合にはその圧力、更に液深が大きい場合は静水圧が加わる）によって定まる。
【0031】
要約すると、酸素移動速度ＯＴＲが決定されると通気量との関係により気相酸素濃度が決定される。
【0032】
気相酸素濃度が決まると、飽和溶存酸素濃度ＤＯ＊が決まり、飽和溶存酸素濃度ＤＯ*との関係で酸素移動速度ＯＴＲが決定される。すなわち、酸素移動速度ＯＴＲと気相酸素濃度と飽和溶存酸素濃度ＤＯ＊とは依存関係を持ったループを形成している。これが飽和溶存酸素濃度ＤＯ＊の計算を非常に複雑にしている。
【0033】
言い換えると、酸素移動速度ＯＴＲの計算のため、飽和溶存酸素濃度ＤＯ＊を事前に予測する上で最も問題なのは、飽和溶存酸素濃度ＤＯ＊を定める気相中の酸素分圧が酸素移動速度ＯＴＲの関数であり、かつ、酸素移動速度ＯＴＲが飽和溶存酸素濃度ＤＯ＊の関数であるという点である。
【0034】
発明者は、この問題を解決して酸素移動速度ＯＴＲ、飽和溶存酸素濃度ＤＯ＊の推定値を得るには、酸素移動速度ＯＴＲと飽和溶存酸素濃度ＤＯ＊の関係を表現する非線形方程式につき、収束計算（繰り返し計算）を実施すれば良いことを見出した。
【0035】
＜溶存酸素濃度ＤＯ＞
溶存酸素濃度ＤＯ（Ｄｉｓｓｏｌｖｅｄ Ｏｘｙｇｅｎ）は、培養液内に溶存する酸素の濃度を示しており、単位は、［ｇ－Ｏ２／Ｌ］である。
溶存酸素濃度ＤＯ（Ｄｉｓｓｏｌｖｅｄ Ｏｘｙｇｅｎ）は、培養液内に溶存する酸素の濃度を示しており、単位は、［ｇ－Ｏ２／Ｌ］である。
【0036】
酸素ガスを溶媒に溶かして溶存酸素にするために、通気撹拌を伴った培養槽（通気撹拌型培養槽）を設けてもよい。通気撹拌型培養槽とは、例えば、培養槽の底部に設けた管から酸素を含む気体（例えば空気）を供給して、培養槽の中央に下垂した撹拌軸の先に付けた羽根などの翼で培養液を撹拌するもので、酸素が微小な気泡になることにより、培養槽内に溶存することができる。
【0037】
＜溶存酸素濃度ＤＯと酸素摂取速度ＯＵＲ＞
溶存酸素濃度ＤＯと酸素摂取速度ＯＵＲの関係について説明する。
培養液中の細胞は、培養液内に溶解している酸素を取り込むが、その活性は溶存酸素濃度ＤＯの関数で表現できる場合が多いことが知られる。単位細胞量当たりの酸素摂取速度（酸素比摂取速度）は、例えば以下のような酵素反応のミカエリス・メンテン式と同じ形式の数式（Ｍｏｎｏｄ式；モノー式）で表現できる。
溶存酸素濃度ＤＯと酸素摂取速度ＯＵＲの関係について説明する。
培養液中の細胞は、培養液内に溶解している酸素を取り込むが、その活性は溶存酸素濃度ＤＯの関数で表現できる場合が多いことが知られる。単位細胞量当たりの酸素摂取速度（酸素比摂取速度）は、例えば以下のような酵素反応のミカエリス・メンテン式と同じ形式の数式（Ｍｏｎｏｄ式；モノー式）で表現できる。
【0038】
【数3】
【0039】
ここで、
ＱＯ２ ：酸素比摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＱＯ２,ｍａｘ：最大比酸素摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
Ｋｍ，Ｏ２ ：飽和定数 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＱＯ２ ：酸素比摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＱＯ２,ｍａｘ：最大比酸素摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
Ｋｍ，Ｏ２ ：飽和定数 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
【0040】
Ｋｍ，Ｏ２は飽和定数であり、最大酸素比摂取速度の１／２の速度を与える時の溶存酸素濃度に相当する。
【0041】
酸素比摂取速度ＱＯ２を表現する式は、溶存酸素濃度ＤＯの関数であればどのような形でもよい。例えば、溶存酸素濃度ＤＯが臨界溶存酸素濃度ＤＯｃｒｉｔより低い場合、酸素比摂取ＱＯ２は溶存酸素濃度ＤＯと比例して減少する形でもよい（非特許文献４：小林猛、本田裕之（２００２）”生物化学工学”（東京化学同人）ｐ．６２参照）。
【0042】
細胞の酸素摂取速度ＯＵＲは、酸素比摂取速度ＱＯ２を用いて、以下の式で示される。
【0043】
【数4】
【0044】
ここで、
ＯＵＲ ：酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２：酸素比摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
Ｘ ：細胞濃度 ［ｇ－ｄｒｙ ｃｅｌｌ／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
ＯＵＲ ：酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２：酸素比摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
Ｘ ：細胞濃度 ［ｇ－ｄｒｙ ｃｅｌｌ／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
【0045】
＜酸素摂取速度ＯＵＲと気相酸素濃度＞
酸素摂取速度ＯＵＲが定まると、更に、呼吸商（Ｒｅｓｐｉｒａｔｏｒｙ Ｑｕｏｔｉｅｎｔ；ＲＱ）より、ＣＯ２生成速度（Ｃａｒｂｏｎ ｄｉｏｘｉｄｅ Ｅｖｏｌｕｔｉｏｎ Ｒａｔｅ；ＣＥＲ）が定まる。呼吸商ＲＱとは、ある時間において生体内で栄養素が分解されてエネルギーに変換するまでの酸素消費量に対する二酸化炭素排出量の体積比のことである。呼吸率、呼吸係数とも呼ばれる。
酸素摂取速度ＯＵＲが定まると、更に、呼吸商（Ｒｅｓｐｉｒａｔｏｒｙ Ｑｕｏｔｉｅｎｔ；ＲＱ）より、ＣＯ２生成速度（Ｃａｒｂｏｎ ｄｉｏｘｉｄｅ Ｅｖｏｌｕｔｉｏｎ Ｒａｔｅ；ＣＥＲ）が定まる。呼吸商ＲＱとは、ある時間において生体内で栄養素が分解されてエネルギーに変換するまでの酸素消費量に対する二酸化炭素排出量の体積比のことである。呼吸率、呼吸係数とも呼ばれる。
【0046】
【数5】
ここで、
ＣＥＲ ：ＣＯ２生成速度 ［ｍｍｏｌ－ＣＯ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＲＱ ：呼吸商 ［－］
ＯＵＲ ：酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
【0047】
酸素摂取速度ＯＵＲとＣＯ２生成速度ＣＥＲが定まると、培養槽に導入する通気ガスに含まれる酸素と二酸化炭素の量との物質収支より、培養槽内気相に含まれる酸素濃度が定まる。
【0048】
以上をまとめると、溶存酸素濃度ＤＯが決定されると数式（２）、数式（３）より酸素摂取速度ＯＵＲが決定される。更に、酸素摂取速度ＯＵＲが決定されると、数式（４）より二酸化炭素生成速度ＣＥＲが決定される。また、酸素摂取速度ＯＵＲとＣＯ２生成速度ＣＥＲは、培養槽内気相酸素濃度にも影響を与える。培養槽内気相酸素濃度が決まると飽和溶存酸素濃度ＤＯ＊が決まり、数式（１）より飽和溶存酸素濃度ＤＯ*との関係で酸素移動速度ＯＴＲが決定される。
【0049】
すなわち、溶存酸素濃度ＤＯは、数式（１）より酸素移動速度ＯＴＲの決定に直接関与しているが、同時に、酸素摂取速度ＯＵＲを介して間接的にも酸素移動速度ＯＴＲの決定に関与している。
【0050】
よって、酸素移動速度ＯＴＲを予測してスケールアップを成功させるためには、飽和溶存酸素濃度ＤＯ*に加え、溶存酸素濃度ＤＯを事前に予測する手段が望まれる。
【0051】
この問題を解決するには、溶存酸素濃度ＤＯによって定まる複数の因子に関して方程式を立て、方程式を解くことが可能になるような何らかの条件を与えた上で、解としての溶存酸素濃度ＤＯの値を得ればよい。
【0052】
具体的には、本実施例の培養槽における酸素移動速度を計算するためのシミュレータは、溶液酸素濃度に基づいて、酸素摂取速度を計算する手段と、飽和溶存酸素濃度と溶存酸素濃度と酸素移動容量係数に基づいて、酸素移動速度を計算する手段と、通気量と酸素移動速度と酸素摂取速度と二酸化炭素生成速度と水蒸気圧に基づいて、反応槽内の気相の酸素濃度を計算する手段と、反応槽内の気相の酸素分圧に基づいて、飽和溶存酸素濃度を計算する手段と、酸素移動速度と酸素摂取速度との間で定常状態を確立するために、酸素移動速度を計算する手段と酸素摂取速度を計算する手段との間でデータを調整する調整手段と、を備えればよい。
【0053】
培養槽内の液相における酸素の物質収支は以下の式で表される。
【0054】
【数6】
【0055】
ここで、
ｔ：時間［ｈ］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＯＴＲ ：酸素移動速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＵＲ ：酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
ｔ：時間［ｈ］
ＤＯ ：溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＯＴＲ ：酸素移動速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＵＲ ：酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＭＷＯ２：酸素の分子量 ［ｇ／ｍｏｌ］
【0056】
一般の培養では、数式（５）の左辺、すなわち溶存酸素濃度ＤＯの経時変化はそれほど急激ではないため、ｄ［ＤＯ］／ｄｔ≒０とおくことが出来る。このとき、酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲの定常状態が成り立っている。酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲの定常状態仮定は、実際に、排気ガス分析法で培養槽のｋＬａを測定する際に用いられており（非特許文献５：石崎文彬（訳）（１９８８）“発酵工学の基礎”（学会出版センター）ｐ．１７６を参照）、現実的な設定である。
【0057】
酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲという条件を与えることで、培養槽周りの物質収支に関する連立方程式の自由度が０（ゼロ）になり、溶存酸素濃度ＤＯについて解が得られるようになる。
【0058】
要約すると、本発明は、以下の２点を計算手順として含む。
（１）培養槽内気相についての物質収支に関する方程式を解く。特に、培養槽内気相の酸素濃度を算出し、それに応じた飽和溶存酸素濃度ＤＯ＊を算出する。
（２）溶存酸素濃度ＤＯに関する方程式をたてて、酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲとなるような定常状態仮定を置くことで、当該溶存酸素濃度ＤＯを算出する。
（１）培養槽内気相についての物質収支に関する方程式を解く。特に、培養槽内気相の酸素濃度を算出し、それに応じた飽和溶存酸素濃度ＤＯ＊を算出する。
（２）溶存酸素濃度ＤＯに関する方程式をたてて、酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲとなるような定常状態仮定を置くことで、当該溶存酸素濃度ＤＯを算出する。
【0059】
本シミュレータで使用している理論について以下に説明する。
まず、ガス物性について説明する。
気体の溶解度は、以下の式で示される。
まず、ガス物性について説明する。
気体の溶解度は、以下の式で示される。
【0060】
【数7】
ここで、溶媒が水（Ｈ２Ｏ）で、気体が酸素（Ｏ２）であるとき、定数ＡＳ＝－６６．７３５３８、ＢＳ＝８７．４７５４７、ＣＳ＝２４．４５２６４である（非特許文献６：日本化学会（編）（２００４）化学便覧改訂５版 基礎編ＩＩ（丸善）ｐ．ＩＩ－１４４）。
また、蒸気圧は、以下のＡｎｔｏｉｎｅ式（アントワン式）により示される。
【0061】
【数8】
ここで、蒸気圧Ｐ、温度Ｔ［Ｋ］、アントワン定数ＡＡ、ＢＡ、ＣＡである。
水の蒸気圧を求める場合、ＡＡ＝２３．１９６４、ＢＡ＝３８１６．４４、ＣＡ＝－４６．１３である（非特許文献７：化学工学会（編）（１９９９）化学工学便覧第６版（丸善）ｐ．１９）。
【0062】
数式（６）、数式（７）の定数に上述した（［００５９］、［００６０］で述べた）数値を代入し、さらに気体の溶解度をモル分率から［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］に換算すると、以下の数式（８）のようになる。なお、酸素の溶解度はきわめて低いため、酸素が溶解した溶液の密度は溶媒の密度に等しいものとして換算を行っている。また、数式（８）では純水に関する数値を代入しているが、培養液に関する数値を代入してもよい。
【数9】
【0063】
ここで、
ｘ２，Ｏ２ ：気相中の酸素分圧＝大気圧での飽和溶存酸素モル分率［－］
Ｔ ：絶対温度 ［Ｋ］
ＤＯ＊ ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量［ｇ／ｍｏｌ］
ＭＷＨ２Ｏ ：水の分子量 ［ｇ／ｍｏｌ］
Ｄ ：溶媒の密度 ［ｇ－ｂｒｏｔｈ／Ｌ－ｂｒｏｔｈ］
ｐｇ，Ｈ２Ｏ：水の蒸気圧 ［Ｐａ］
ｘ２，Ｏ２ ：気相中の酸素分圧＝大気圧での飽和溶存酸素モル分率［－］
Ｔ ：絶対温度 ［Ｋ］
ＤＯ＊ ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量［ｇ／ｍｏｌ］
ＭＷＨ２Ｏ ：水の分子量 ［ｇ／ｍｏｌ］
Ｄ ：溶媒の密度 ［ｇ－ｂｒｏｔｈ／Ｌ－ｂｒｏｔｈ］
ｐｇ，Ｈ２Ｏ：水の蒸気圧 ［Ｐａ］
【0064】
＜気相部の物質収支式＞
培養槽に入る直前の通気ガスについて説明する。
培養槽に入る直前の通気ガスが水蒸気を含まない乾燥気体とすると、通気ガスの酸素流量Ｆｇ，Ｏ２，ｆｅｅｄ、ＣＯ２流量Ｆｇ，ＣＯ２，ｆｅｅｄ、窒素流量Ｆｇ，Ｎ２，ｆｅｅｄ、は、通気ガス流量Ｆｇ，ｆｅｅｄとそれぞれのガス成分濃度Ｃｇとから、以下の計算式により算出される。
培養槽に入る直前の通気ガスについて説明する。
培養槽に入る直前の通気ガスが水蒸気を含まない乾燥気体とすると、通気ガスの酸素流量Ｆｇ，Ｏ２，ｆｅｅｄ、ＣＯ２流量Ｆｇ，ＣＯ２，ｆｅｅｄ、窒素流量Ｆｇ，Ｎ２，ｆｅｅｄ、は、通気ガス流量Ｆｇ，ｆｅｅｄとそれぞれのガス成分濃度Ｃｇとから、以下の計算式により算出される。
【0065】
【数10】
【0066】
ここで、
Ｆｇ，ｆｅｅｄ ：通気ガス流量 ［Ｌ－ｇａｓ／ｍｉｎ］
Ｆｇ，Ｏ２，ｆｅｅｄ ：通気ガス中の酸素流量 ［Ｌ－Ｏ２ ｇａｓ／ｍｉｎ］
Ｆｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２流量 ［Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ］
Ｆｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素流量 ［Ｌ－Ｎ２ ｇａｓ／ｍｉｎ］
Ｃｇ，Ｏ２，ｆｅｅｄ：通気ガス中の酸素濃度 ［％］
Ｃｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２濃度 ［％］
Ｃｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素濃度 ［％］
Ｆｇ，ｆｅｅｄ ：通気ガス流量 ［Ｌ－ｇａｓ／ｍｉｎ］
Ｆｇ，Ｏ２，ｆｅｅｄ ：通気ガス中の酸素流量 ［Ｌ－Ｏ２ ｇａｓ／ｍｉｎ］
Ｆｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２流量 ［Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ］
Ｆｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素流量 ［Ｌ－Ｎ２ ｇａｓ／ｍｉｎ］
Ｃｇ，Ｏ２，ｆｅｅｄ：通気ガス中の酸素濃度 ［％］
Ｃｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２濃度 ［％］
Ｃｇ，Ｎ２，ｆｅｅｄ：通気ガス中の窒素濃度 ［％］
【0067】
培養槽に入った直後、通気ガスに飽和水蒸気が加わった入口ガスについて説明する。入口ガスの各成分の流量、濃度、分圧は数式（１０）～数式（１２）で示される。入口ガスの酸素分圧より、培養液底部の飽和溶存酸素濃度は数式（１３）で示される。
【0068】
【数11】
【0069】
【数12】
【0070】
【数13】
【0071】
【数14】
ここで、
Ｆｇ, ｉｎ ：入口ガス流量 [Ｌ－ｇａｓ／ｍｉｎ]
Ｆｇ，ｆｅｅｄ ：通気ガス流量 [Ｌ－ｇａｓ／ｍｉｎ］
Ｆｇ, Ｏ２, ｉｎ ：入口ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ]
Ｆｇ，Ｏ２，ｆｅｅｄ ：通気ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ］
Ｆｇ, ＣＯ２, ｉｎ ：入口ガス中のＣＯ２流量 [Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ]
Ｆｇ，ＣＯ２，ｆｅｅｄ：通気ガス中のＣＯ２流量 [Ｌ－ＣＯ２ｇａｓ／ｍｉｎ］
Ｆｇ, Ｎ２, ｉｎ ：入口ガス中の窒素流量 [Ｌ－Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ，Ｎ２，ｆｅｅｄ ：通気ガス中の窒素流量 [Ｌ－Ｎ２ｇａｓ／ｍｉｎ］
Ｆｇ, Ｈ２Ｏ, ｉｎ ：入口ガス中の水蒸気流量 [Ｌ－Ｈ２Ｏ ｇａｓ／ｍｉｎ]
Ｃｇ, Ｏ２, ｉｎ ：入口ガス中の酸素濃度 [％]
Ｃｇ, ＣＯ２, ｉｎ ：入口ガス中のＣＯ２濃度 [％]
Ｃｇ, Ｎ２, ｉｎ ：入口ガス中の窒素濃度 [％]
Ｃｇ, Ｈ２Ｏ, ｉｎ ：入口ガス中の水蒸気濃度 [％]
Ｐｂｏｔｔｏｍ ：培養液底部の全圧 [Ｐａ]
ｐｇ, Ｏ２, ｉｎ ：入口ガス中の酸素分圧 [Ｐａ]
ｐｇ, ＣＯ２, ｉｎ ：入口ガス中のＣＯ２分圧 [Ｐａ]
ｐｇ, Ｎ２, ｉｎ ：入口ガス中の窒素分圧 [Ｐａ]
ｐｇ, Ｈ２Ｏ, ｉｎ ：入口ガス中の水蒸気分圧 [Ｐａ]
ｐｇ，Ｈ２Ｏ ：水の蒸気圧 [Ｐａ］
ＤＯ＊ ａｔｍ ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ* ｂｏｔｔｏｍ ：培養液底部の飽和溶存酸素濃度 [ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ]
Ｐａｔｍ ：大気圧 [Ｐａ]
【0072】
ガスに関する物質移動について説明する。
酸素移動速度ＯＴＲは、例えば容量が５０Ｌ以下の小型培養槽では数式（１）で表される（再掲）。大型培養槽では、培養液全体の酸素移動速度ＯＴＲは数式（１）を修正して、数式（１４）で示される。また、大型培養槽の培養液の最表面と最底部では、数式（１）について、培養液の最表面と最底部それぞれの飽和溶存酸素濃度ＤＯ*と溶存酸素濃度ＤＯを適用することで、数式（１５）で示される。
酸素移動速度ＯＴＲは、例えば容量が５０Ｌ以下の小型培養槽では数式（１）で表される（再掲）。大型培養槽では、培養液全体の酸素移動速度ＯＴＲは数式（１）を修正して、数式（１４）で示される。また、大型培養槽の培養液の最表面と最底部では、数式（１）について、培養液の最表面と最底部それぞれの飽和溶存酸素濃度ＤＯ*と溶存酸素濃度ＤＯを適用することで、数式（１５）で示される。
【0073】
【数15】
【数16】
【0074】
【数17】
【0075】
ここで、
ｋＬａ ：酸素移動容量係数 ［ｈ－１］
ＯＴＲ ：培養液全体の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＴＲｂｏｔｔｏｍ：培養液底部の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＴＲｔｏｐ ：培養液表面の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＤＯ* ｂｏｔｔｏｍ：培養液底部の飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯｂｏｔｔｏｍ：培養液底部の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ＊ ｔｏｐ ：培養液表面の飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯｔｏｐ ：培養液表面の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量 ［ｇ／ｍｏｌ］
ｋＬａ ：酸素移動容量係数 ［ｈ－１］
ＯＴＲ ：培養液全体の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＴＲｂｏｔｔｏｍ：培養液底部の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＴＲｔｏｐ ：培養液表面の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＤＯ* ｂｏｔｔｏｍ：培養液底部の飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯｂｏｔｔｏｍ：培養液底部の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ＊ ｔｏｐ ：培養液表面の飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯｔｏｐ ：培養液表面の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量 ［ｇ／ｍｏｌ］
【0076】
排ガスについて説明する。排ガスとは、入口ガスから消費された酸素Ｏ２を引き、発生した二酸化酸素ＣＯ２を足したものである。なお、排ガスは水蒸気で飽和しているものとする。
【0077】
【数18】
【0078】
【数19】
【数20】
【0079】
【数21】
【0080】
【数22】
【0081】
ここで、
Ｆｇ, ｏｕｔ ：排ガス流量 [Ｌ－ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｉｎ ：入口ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ]
ＯＴＲ ：培養液全体の酸素移動速度［ｍｍｏｌ－Ｏ2／（Ｌ－ｂｒｏｔｈ・ｈ）］
Ｖｌ ：培養液量 [Ｌ－ｂｒｏｔｈ]
Ｆｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２流量 [Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ]
Ｆｇ, ＣＯ２, ｉｎ ：入口ガス中のＣＯ２流量 [Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ]
ＣＥＲ ：培養液全体のＣＯ２生成速度［ｍｍｏｌ－ＣＯ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
Ｆｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素流量 [Ｌ－Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｎ２, ｉｎ ：入口ガス中の窒素流量 [Ｌ－Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気流量 [Ｌ－Ｈ２Ｏ ｇａｓ／ｍｉｎ]
Ｃｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素濃度 [％]
Ｃｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２濃度 [％]
Ｃｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素濃度 [％]
Ｃｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気濃度 [％]
Ｐｔｏｐ ：培養液表面の全圧 [Ｐａ]
ｐｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素分圧 [Ｐａ]
ｐｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２分圧 [Ｐａ]
ｐｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素分圧 [Ｐａ]
ｐｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気分圧 [Ｐａ]
ｐｇ，Ｈ２Ｏ ：水の蒸気圧 [Ｐａ］
ＤＯ＊ ｔｏｐ ：培養液表面の飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ＊ ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
Ｐａｔｍ ：大気圧 [Ｐａ]
Ｆｇ, ｏｕｔ ：排ガス流量 [Ｌ－ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｏ２, ｉｎ ：入口ガス中の酸素流量 [Ｌ－Ｏ２ ｇａｓ／ｍｉｎ]
ＯＴＲ ：培養液全体の酸素移動速度［ｍｍｏｌ－Ｏ2／（Ｌ－ｂｒｏｔｈ・ｈ）］
Ｖｌ ：培養液量 [Ｌ－ｂｒｏｔｈ]
Ｆｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２流量 [Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ]
Ｆｇ, ＣＯ２, ｉｎ ：入口ガス中のＣＯ２流量 [Ｌ－ＣＯ２ ｇａｓ／ｍｉｎ]
ＣＥＲ ：培養液全体のＣＯ２生成速度［ｍｍｏｌ－ＣＯ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
Ｆｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素流量 [Ｌ－Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｎ２, ｉｎ ：入口ガス中の窒素流量 [Ｌ－Ｎ２ ｇａｓ／ｍｉｎ]
Ｆｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気流量 [Ｌ－Ｈ２Ｏ ｇａｓ／ｍｉｎ]
Ｃｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素濃度 [％]
Ｃｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２濃度 [％]
Ｃｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素濃度 [％]
Ｃｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気濃度 [％]
Ｐｔｏｐ ：培養液表面の全圧 [Ｐａ]
ｐｇ, Ｏ２, ｏｕｔ ：排ガス中の酸素分圧 [Ｐａ]
ｐｇ, ＣＯ２, ｏｕｔ ：排ガス中のＣＯ２分圧 [Ｐａ]
ｐｇ, Ｎ２, ｏｕｔ ：排ガス中の窒素分圧 [Ｐａ]
ｐｇ, Ｈ２Ｏ, ｏｕｔ ：排ガス中の水蒸気分圧 [Ｐａ]
ｐｇ，Ｈ２Ｏ ：水の蒸気圧 [Ｐａ］
ＤＯ＊ ｔｏｐ ：培養液表面の飽和溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＤＯ＊ ａｔｍ：気相中の酸素分圧＝大気圧での飽和溶存酸素濃度［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
Ｐａｔｍ ：大気圧 [Ｐａ]
【0082】
細胞のガス代謝モデルについて説明する。細胞の酸素消費が酸素濃度に関するＭｏｎｏｄ（モノー）式に従うと、微生物による酸素比摂取速度ＱＯ２は数式（２）（再掲）で示される。大型培養槽では培養液の最表面と最底部で溶存酸素濃度ＤＯが異なることを考慮すると、以下の式で示される。
【0083】
【数23】
【0084】
【数24】
ここで、
ＱＯ２,ｍａｘ：最大比酸素摂取速度 ［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
Ｋｍ，Ｏ２ ：飽和定数 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＱＯ２,ｂｏｔｔｏｍ：培養液底部の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＤＯｂｏｔｔｏｍ ：培養液底部の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
ＱＯ２，ｔｏｍ ：培養液表面の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＤＯｔｏｐ ：培養液表面の溶存酸素濃度 ［ｇ－Ｏ２／Ｌ－ｂｒｏｔｈ］
【0085】
微生物の酸素摂取速度ＯＵＲは、酸素比摂取速度ＱＯ２を用いて、数式（３）（再掲）で示されるが、数式（２１）で示したように、大型培養槽では培養液の最表面と最底部でＱＯ２が異なることを考慮すると、数式（２２）で示される。
【0086】
【数25】
【0087】
【数26】
【0088】
ここで、
ＯＵＲｂｏｔｔｏｍ：培養液底部の酸素摂取速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２,ｂｏｔｔｏｍ：培養液底部の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＯＵＲｔｏｐ：培養液表面の酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２, ｔｏｐ：培養液表面の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
Ｘ ：細胞濃度 ［ｇ－ｄｒｙ ｃｅｌｌ／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量 ［ｇ／ｍｏｌ］
ＯＵＲｂｏｔｔｏｍ：培養液底部の酸素摂取速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２,ｂｏｔｔｏｍ：培養液底部の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
ＯＵＲｔｏｐ：培養液表面の酸素摂取速度 ［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＱＯ２, ｔｏｐ：培養液表面の酸素比摂取速度［ｇ－Ｏ２／（ｇ－ｄｒｙ ｃｅｌｌ・ｈ）］
Ｘ ：細胞濃度 ［ｇ－ｄｒｙ ｃｅｌｌ／Ｌ－ｂｒｏｔｈ］
ＭＷＯ２ ：酸素の分子量 ［ｇ／ｍｏｌ］
【0089】
また、ＣＯ２生成速度ＣＥＲ（Ｃａｒｂｏｎ ｄｉｏｘｉｄｅ Ｅｖｏｌｕｔｉｏｎ Ｒａｔｅ）は、培養液の最表面と最底部それぞれのＯＵＲではなく、培養槽全体での酸素摂取速度ＯＵＲを用いれば、数式（４）で示される（再掲）。
【数27】
【0090】
次に、定常状態仮定について説明する。本実施例における定常状態仮定とは、一定の過渡状態（後述）が発生した後に、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲとが等しくなるという定常状態を仮定することである。
【0091】
【数28】
また、本実施例のシミュレータにおいては、大型培養槽の場合、培養液底部と培養液表面とでは酸素移動速度ＯＴＲ、酸素摂取速度ＯＵＲの値に差異が生じることに着目して、以下の計算式を用いる。
【0092】
【数29】
ここで、
ＯＴＲｂｏｔｔｏｍ：培養液底部の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＵＲｂｏｔｔｏｍ：培養液底部の酸素摂取速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＴＲｔｏｐ ：培養液表面の酸素移動速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
ＯＵＲｔｏｐ ：培養液表面の酸素摂取速度［ｍｍｏｌ－Ｏ２／（Ｌ－ｂｒｏｔｈ・ｈ）］
【0093】
定常状態を仮定できる理由を、概念的に説明する。
【0094】
まず、供給過剰の例について説明する。供給過剰とは、物理的な酸素溶解速度である酸素移動速度ＯＴＲが、細胞による酸素消費速度である酸素摂取速度ＯＵＲよりも多くなった過渡的な場合である（酸素移動速度ＯＴＲ＞酸素摂取速度ＯＵＲ）。この場合、培養液内においては、溶存酸素濃度ＤＯが増加していく。溶存酸素濃度ＤＯが増加すると、飽和溶存酸素濃度ＤＯ＊と溶存酸素濃度ＤＯとの差が小さくなる。当該差が小さくなると、数式（１）（再掲）により、酸素移動速度ＯＴＲが小さくなる。すなわち、当初多かった酸素移動速度ＯＴＲが減少して、酸素摂取速度ＯＵＲに等しくなるように作用する。
【0095】
【数30】
【0096】
次に、消費過剰の例について説明する。消費過剰とは、酸素移動速度ＯＴＲが、細胞による酸素摂取速度ＯＵＲよりも少ない場合である（酸素移動速度ＯＴＲ＜酸素摂取速度ＯＵＲ）。この場合、培養液内においては、溶存酸素濃度ＤＯが減少していく。溶存酸素濃度ＤＯが減少すると、数式（１）により、飽和溶存酸素濃度ＤＯ＊と溶存酸素濃度ＤＯとの差が大きくなる。当該差が大きくなると、数式（１）により、酸素移動速度ＯＴＲが大きくなる。すなわち、当初少なかった酸素移動速度ＯＴＲが増加して、酸素摂取速度ＯＵＲに等しくなるように作用する。
よって、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲの値に差異があっても、系は定常状態へ向かう。
よって、酸素移動速度ＯＴＲと酸素摂取速度ＯＵＲの値に差異があっても、系は定常状態へ向かう。
【0097】
本実施例のシミュレータでは、上述した一連の式（数式１～数式２４）を連立非線形方程式として収束計算（繰り返し計算）で解く。求解には各種の数値計算ソフトウェアが利用できる。例えば、ＥＱＵＡＴＲＡＮ（商標）－Ｇ（イコートラン；オメガシミュレーション社製）を用いてもよい。なお、発明者の知る限りこの様な考え方を培養槽シミュレーションに適用した例はない。
【0098】
本実施例のシミュレータによると、排ガスの酸素濃度に応じた飽和溶存酸素濃度ＤＯ＊を算出していることで、従来よりも正確に酸素移動速度ＯＴＲを推算できる。
【0099】
特に、酸素移動容量係数ｋＬａが同じであっても通気量が異なる場合に酸素移動速度ＯＴＲが変化する状況を予測可能である。
【0100】
また、従来、酸素移動容量係数ｋＬａが不足してスケールアップ不能と判断されてきた条件でも、酸素移動速度ＯＴＲが十分なのであればスケールアップできる可能性について検討できる。
【0101】
そのほか、細胞の呼吸特性の違いによる酸素移動速度ＯＴＲの変化を計算できる。
【実施例1】
【0102】
図３に、本実施例によるＯＴＲシミュレータ１０００の機能ブロック図を示す。なお、本実施例の機能ブロックは、ハードウェアでもソフトウェアでもハードウェア及びソフトウェアの組み合わせでも実現可能である。
【0103】
ＯＴＲシミュレータ１０００は、通気量設定手段１０１０、Ｏ２濃度計算手段１０２０と、
ＯＴＲ計算手段１０３０と、データ調整手段１０４０と、ＯＵＲ計算手段１０５０と、ＤＯ＊計算手段１０６０と、ＤＯ計算手段１０７０と、操作条件ＤＢ１１１０と、入力変数ＤＢ１１２０と、出力変数ＤＢ１１３０と、制御手段１２００とから構成される。
ＯＴＲ計算手段１０３０と、データ調整手段１０４０と、ＯＵＲ計算手段１０５０と、ＤＯ＊計算手段１０６０と、ＤＯ計算手段１０７０と、操作条件ＤＢ１１１０と、入力変数ＤＢ１１２０と、出力変数ＤＢ１１３０と、制御手段１２００とから構成される。
【0104】
通気量設定手段１０１０は、通気量を設定する手段である。
【0105】
Ｏ２濃度計算手段１０２０は、他の手段やデータベース（ＤＢ：ＤａｔａＢａｓｅ）を用いて、酸素濃度を計算する手段である。
【0106】
ＯＴＲ計算手段１０３０は、他の手段やＤＢを用いて、酸素移動速度ＯＴＲを計算する手段である。
【0107】
データ調整手段１０４０は、他の手段やＤＢで使用されているデータ（値）を調整する手段であり、本実施例においては、ＯＴＲ計算手段で計算された酸素移動速度ＯＴＲの値と、ＯＵＲ計算手段で計算された酸素摂取速度ＯＵＲの値とが定常状態仮定になるようにする手段である。一実施例によると酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲの計算式が成立するように構成されてもよい。
【0108】
ＯＵＲ計算手段１０５０は、他の手段やＤＢを用いて、酸素摂取速度ＯＵＲを計算する手段である。
【0109】
ＤＯ＊計算手段１０６０は、他の手段やＤＢを用いて、飽和溶存酸素濃度ＤＯ＊を計算する手段である。
【0110】
ＤＯ計算手段１０７０は、他の手段やＤＢを用いて、溶存酸素濃度ＤＯを計算する手段である。
【0111】
装置操作条件ＤＢ１１１０は、収束計算をするにあたり必要となる計算式（方程式）に適用するための各パラメータを記憶するＤＢである。
【0112】
必要なパラメータは、例えば、以下のようなものである。
（１） 培養槽に関する諸元
・培養液量、培養槽深さ
・通気量
・温度、培養液表面圧力
・通気ガスの酸素濃度
・酸素移動容量係数（ｋＬａ）
（２）細胞に関する諸元
・細胞濃度
・細胞の酸素消費・ＣＯ２生成パラメータ
－ 酸素に対する飽和定数
－ 最大比酸素摂取速度
－ 呼吸商
（１） 培養槽に関する諸元
・培養液量、培養槽深さ
・通気量
・温度、培養液表面圧力
・通気ガスの酸素濃度
・酸素移動容量係数（ｋＬａ）
（２）細胞に関する諸元
・細胞濃度
・細胞の酸素消費・ＣＯ２生成パラメータ
－ 酸素に対する飽和定数
－ 最大比酸素摂取速度
－ 呼吸商
【0113】
入力変数ＤＢ１１２０は、収束計算をするにあたり必要となる入力条件を記憶するデータベース（ＤＢ）である。
【0114】
出力変数ＤＢ１１３０は、収束計算によって算出された値を格納するデータベース（ＤＢ）である。
【0115】
制御手段１２００は、ＯＴＲシミュレータ内部の各手段や各ＤＢを制御したり、ＯＴＲシミュレータ外部と信号を送受信したりする。
また、ＯＴＲシミュレータ１０００は、ユーザ端末１３００に接続されている。ユーザは、ユーザ端末１３００を介して、ＯＴＲシミュレータを作動させるために必要な変数を入力する。入力された変数は、入力変数ＤＢ１１２０に記憶される。また、ユーザ端末１３００は、ＯＴＲシミュレータを作動させた結果として得られた出力を表示する。この出力は、ＯＴＲシミュレータ１０００の出力変数ＤＢ１１３０に記憶された値である。
また、ＯＴＲシミュレータ１０００は、ユーザ端末１３００に接続されている。ユーザは、ユーザ端末１３００を介して、ＯＴＲシミュレータを作動させるために必要な変数を入力する。入力された変数は、入力変数ＤＢ１１２０に記憶される。また、ユーザ端末１３００は、ＯＴＲシミュレータを作動させた結果として得られた出力を表示する。この出力は、ＯＴＲシミュレータ１０００の出力変数ＤＢ１１３０に記憶された値である。
【0116】
酸素移動速度ＯＴＲの計算方法の一例を図４のフローチャートを用いて説明する。
【0117】
Ｓ３０１０では、排ガス酸素濃度を計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【0118】
【数31】
Ｓ３０２０では、飽和溶存酸素濃度ＤＯ＊を計算する。特に、上述した連立方程式において、以下の計算式について計算する。
【0119】
【数32】
【0120】
Ｓ３０３０では、溶存酸素濃度ＤＯを計算する。上述した連立方程式において、例えば、以下の式を満たす溶存酸素濃度ＤＯを計算する。
【0121】
【数33】
【0122】
Ｓ３０４０では、酸素摂取速度ＯＵＲを計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【0123】
【数34】
【0124】
Ｓ３０５０では、定常状態を仮定して酸素移動速度ＯＴＲを計算する。特に、上述した連立方程式において、例えば、以下の計算式について計算する。
【0125】
【数35】
【0126】
Ｓ３０６０では、Ｓ３０１０－Ｓ３０６０で計算された各値が収束しているか否かを判定する。もし、収束していれば、Ｓ３０８０へ進み、収束していなければ、Ｓ３０１０へ戻り、計算を続ける。
【0127】
Ｓ３０８０では、収束した値を表示する。特に、酸素移動速度ＯＴＲの計算結果を表示する。
【0128】
なお、上述した連立方程式を計算することによって、様々な数値が算出されるが、本実施例において、一番重要な値は、酸素移動速度ＯＴＲの値である。
【0129】
なお、上述したフローチャートでは、Ｓ３０１０－Ｓ３０８０が逐次的に計算されるように説明したが、各変数は相互に依存しているために、上述した処理を並列に実行するように構成されている。
【0130】
本実施例のシミュレータは、以下のような点を考慮している。
・気相に関する物質収支から、気相部（培養排気）の酸素分圧を求める。
・細胞による酸素消費およびＣＯ２生成による培養排気の酸素分圧低下を考慮したうえで、ＯＴＲを計算する。
・（主に大型培養槽で行われる）加圧の影響を考慮する。
・液深による静水圧を考慮する。
・水蒸気分圧を考慮する。
・培養液の最表面と最深部の溶存酸素濃度分布を考慮する。
・気相に関する物質収支から、気相部（培養排気）の酸素分圧を求める。
・細胞による酸素消費およびＣＯ２生成による培養排気の酸素分圧低下を考慮したうえで、ＯＴＲを計算する。
・（主に大型培養槽で行われる）加圧の影響を考慮する。
・液深による静水圧を考慮する。
・水蒸気分圧を考慮する。
・培養液の最表面と最深部の溶存酸素濃度分布を考慮する。
【0131】
また、本実施例のシミュレータは、以下のような点を前提とする。
・定常状態シミュレータであり、系の時間発展は計算しない。
・ガスホールドアップによる培養液高さの変化は考慮していない。
・細胞濃度は一定としている。
・細胞濃度の培養液内での分布は考慮しない（培養液の最表面と最底部の細胞濃度を同一と仮定）。
・ｋＬａの培養液内での分布は考慮しない（培養液の最表面と最底部のｋＬａを同一と仮定）。
・培養槽のいたるところで定常状態（酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲ）が成立していることを仮定している。
・通気は培養槽の底部から行うものとする。
・通気は乾燥気体（水蒸気を含まない）とする。
・気相部は培養槽に入ると直ちに水蒸気が飽和するものとする。
・細胞によるＣＯ２消費は考慮しない。
・定常状態シミュレータであり、系の時間発展は計算しない。
・ガスホールドアップによる培養液高さの変化は考慮していない。
・細胞濃度は一定としている。
・細胞濃度の培養液内での分布は考慮しない（培養液の最表面と最底部の細胞濃度を同一と仮定）。
・ｋＬａの培養液内での分布は考慮しない（培養液の最表面と最底部のｋＬａを同一と仮定）。
・培養槽のいたるところで定常状態（酸素移動速度ＯＴＲ＝酸素摂取速度ＯＵＲ）が成立していることを仮定している。
・通気は培養槽の底部から行うものとする。
・通気は乾燥気体（水蒸気を含まない）とする。
・気相部は培養槽に入ると直ちに水蒸気が飽和するものとする。
・細胞によるＣＯ２消費は考慮しない。
【実施例2】
【0132】
図５は、本発明の第二実施例によるＯＴＲシミュレータの機能ブロック図を示す。第１実施例（図３）に示した実施例と比較すると、ユーザ端末１３００が、撹拌装置、センサ等に置き換わっている。第１実施例では、ユーザ端末１３００で所定の値を入力していたが、第２実施例では、センサが検知した値を使って、シミュレートすることができる。別の実施例として、センサが検知した値とユーザが入力する値を任意に選択できるように構成されてもよい。
【0133】
本実施例のシミュレータの活用例を示す。
培養槽のサイズの影響について、本実施例のシミュレータでの計算結果を図８に示す。ここで、培養槽は、ミニジャー（小型槽）と５０ｍ３タンク（大型槽）との２つについてシミュレーションをした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］（小型槽）ｏｒ３０［ｍ３］（大型槽）
・培養液高さ＝０．１２（小型槽）ｏｒ４．５４［ｍ］（大型槽）
・温度＝３０［℃］
・圧力＝大気圧（小型槽）ｏｒ０．０３ＭＰａＧ（大型槽）
・通気気体＝大気（酸素濃度２０．９％）
・通気量＝１［ｖｖｍ］
同じ酸素移動容量係数ｋＬａであっても、５０ｍ３タンクの方が酸素移動速度ＯＴＲが大きくなることが、シミュレーション結果からわかる。
培養槽のサイズの影響について、本実施例のシミュレータでの計算結果を図８に示す。ここで、培養槽は、ミニジャー（小型槽）と５０ｍ３タンク（大型槽）との２つについてシミュレーションをした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］（小型槽）ｏｒ３０［ｍ３］（大型槽）
・培養液高さ＝０．１２（小型槽）ｏｒ４．５４［ｍ］（大型槽）
・温度＝３０［℃］
・圧力＝大気圧（小型槽）ｏｒ０．０３ＭＰａＧ（大型槽）
・通気気体＝大気（酸素濃度２０．９％）
・通気量＝１［ｖｖｍ］
同じ酸素移動容量係数ｋＬａであっても、５０ｍ３タンクの方が酸素移動速度ＯＴＲが大きくなることが、シミュレーション結果からわかる。
【0134】
本実施例のシミュレータの活用例を示す。
同じ酸素移動容量係数ｋＬａであっても、通気量が異なる場合の影響について、本実施例のシミュレータでの計算結果を図９に示す。ここで、培養槽はミニジャー（小型槽）とし、通気量は０．１、０．５、１．０、２．０ｖｖｍとした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］
・培養液高さ＝０．１２［ｍ］
・温度＝３０［℃］
・圧力＝大気圧
・通気気体＝大気（酸素濃度２０．９％）
同じ酸素移動容量係数ｋＬａであっても、通気量が低下するのに伴って酸素移動速度ＯＴＲが低下することが、シミュレーション結果からわかる。
同じ酸素移動容量係数ｋＬａであっても、通気量が異なる場合の影響について、本実施例のシミュレータでの計算結果を図９に示す。ここで、培養槽はミニジャー（小型槽）とし、通気量は０．１、０．５、１．０、２．０ｖｖｍとした。
そして、培養槽パラメータは、以下のとおりである。
・培養液量＝１．５［Ｌ］
・培養液高さ＝０．１２［ｍ］
・温度＝３０［℃］
・圧力＝大気圧
・通気気体＝大気（酸素濃度２０．９％）
同じ酸素移動容量係数ｋＬａであっても、通気量が低下するのに伴って酸素移動速度ＯＴＲが低下することが、シミュレーション結果からわかる。
【0135】
図７は、本実施例のＯＴＲシミュレータの活用例のイメージ（スケールアップ）を示す。
（１）ラボスケールの実験で、ＯＴＲ、剪断力と生産性の関係についてデータを収集。
（２）（１）のデータより等高線マップを作成。
（３）ＯＴＲシミュレータで、実機の操作可能領域を推算。
（４）（２）と（３）を重ね、生産性を最大化する操作条件と、そのとき期待される生産性を予測。
ここで、図６の出力変数の中で、ＯＴＲｔｏｔａｌが、培養槽の性能評価に用いられ、それ以外の変数は、参考値である。なお、ＯＴＲｔｏｔａｌは、培養槽全体のＯＴＲを示しており、単に、ＯＴＲと表記してもよい。
なお、図６の装置操作条件の入力変数である培養槽底部圧力Ｐｂｏｔｔｏｍは、培養槽直径を入力して簡易計算するように構成してもよい。すなわち、培養槽を円筒形として培養液量から液深を求め、液深・培養槽密度・培養液上部圧力から培養液底部の静水圧を求めるように構成してもよい。また、パラメータの入力変数（培養槽の装置パラメータ）である酸素移動容量係数ｋＬａは、攪拌速度、通気流量、装置形状、培養液物性等の関数であるが、ユーザが数値を入力するように構成してもよい。
（１）ラボスケールの実験で、ＯＴＲ、剪断力と生産性の関係についてデータを収集。
（２）（１）のデータより等高線マップを作成。
（３）ＯＴＲシミュレータで、実機の操作可能領域を推算。
（４）（２）と（３）を重ね、生産性を最大化する操作条件と、そのとき期待される生産性を予測。
ここで、図６の出力変数の中で、ＯＴＲｔｏｔａｌが、培養槽の性能評価に用いられ、それ以外の変数は、参考値である。なお、ＯＴＲｔｏｔａｌは、培養槽全体のＯＴＲを示しており、単に、ＯＴＲと表記してもよい。
なお、図６の装置操作条件の入力変数である培養槽底部圧力Ｐｂｏｔｔｏｍは、培養槽直径を入力して簡易計算するように構成してもよい。すなわち、培養槽を円筒形として培養液量から液深を求め、液深・培養槽密度・培養液上部圧力から培養液底部の静水圧を求めるように構成してもよい。また、パラメータの入力変数（培養槽の装置パラメータ）である酸素移動容量係数ｋＬａは、攪拌速度、通気流量、装置形状、培養液物性等の関数であるが、ユーザが数値を入力するように構成してもよい。
【0136】
以上のように本発明の実施態様について説明したが、上述の説明に基づいて当業者にとって種々の代替例、修正又は変形が可能であり、本発明はその趣旨を逸脱しない範囲で前述の種々の代替例、修正又は変形を包含するものである。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
