(57)【特許請求の範囲】
【請求項１】
通気攪拌培養槽内の溶質ガス濃度分布を最適に制御する方法において、前記槽内の１つ以上の点で溶存酸素濃度を計測し、前記槽のガス移動容量係数（ｋｌａ）の分布を下記の式によって求め、前記分布を用いて溶存酸素濃度分布を演算し、前記計測された溶存酸素濃度を用いて前記溶存酸素濃度分布を補正し、前記補正された溶存酸素濃度分布のうち培養条件から外れている領域があるかどうかを判断し、培養条件から外れている領域がある場合は溶存酸素濃度分布があらかじめ定められた値になるまで前記培養槽の攪拌回転数、通気量、通気ガス中の酸素分圧あるいは内圧を調整することを特徴とする通気攪拌培養槽の運転制御方法。
ｋｌａ=ｆ（Ｓｃ，ε，ν，α，Ｒｂ）
ただし、ｋｌａ：ガス移動容量係数（１／ｈ）、ｋｌ：液境膜移動係数（ｍ／ｓ）、a ：単位体積当たりの気液接触面積（１／ｍ）、Ｓｃ：シュミット数 （＝ν／Ｄｇ）、ε：乱流エネルギー散逸速度（ｍ／ｓ３）、ν：動粘度（ｍ２／ｓ）、Ｄｇ：液相中のガスの拡散係数（ｍ２／ｓ）、α：ガスホールドアップ、Ｒｂ：気泡半径（ｍ）である。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は溶質ガスを液中に溶解させまたは液中で発生したガスを液の外に追い出す必要のある通気攪拌槽の運転制御方法に関する。
【０００２】
【従来の技術】
攪拌槽内での化学反応や培養槽内での細胞の呼吸などにより、酸素の消費や二酸化炭素の生成が起きる槽では、最適な条件で化学反応や培養を行うために、槽内の全ての領域で酸素濃度または二酸化炭素濃度が最適な範囲に収まるように、消費した酸素の補給もしくは生成した二酸化炭素の除去が必要である。
【０００３】
従来、通気撹拌槽内の溶存酸素や溶存二酸化炭素の濃度分布を知るためにはCytotechnology,22,p87-94,1996,"Homogenisation and oxygen transfer rates in large agitated and sparged animal cell bioreactors:Some implications for growth and production"．Alvin W.Nienow etc. で論じられているように、実際の通気撹拌槽を用いて濃度分布を予想するしかなく、全く未知の通気撹拌槽内の濃度分布を定量的に予測することは不可能であった。
【０００４】
一方、攪拌槽内部の流れに関する数値解析では、kla（ガス移動容量係数）の分布まで考慮に入れた物質濃度分布の計算はされていなかった。
【０００５】
また、溶存酸素濃度を制御する方法として、特許第1552563号がある。これは溶存酸素濃度を測定するセンサーを培養槽に設置して、その指示値に合わせて培養槽の運転条件を制御する方法である。
【０００６】
【発明が解決しようとする課題】
槽内の全ての領域において、溶質ガスの濃度を最適にすることを目的とした槽の詳細構造及び運転条件の検討を行うには、槽の構造や運転条件の変化が反映されるような方法で濃度を求める必要がある。従来のような実際の槽を用いて濃度分布を測定する方法では、構造を変化させた効果を知るためには実際に槽を製作するしか方法はなかった。また、非接触で定量的に濃度を測定することはできないので、槽内分布の正確な測定は困難という問題もある。
【０００７】
また上記特許のような場合は、測定点における値しか評価できないため、槽内分布が適切かどうかは判断できない、あるいは不明であった。すなわち攪拌槽の溶質ガス濃度を制御するのに槽内の特定の点における測定値だけを直接用いる方法では、測定している点では濃度が最適であっても、他の点では溶質ガス濃度の過不足が生じている可能性がある。したがって、溶質ガス濃度を槽内の全ての領域で最適な範囲に収めるためには槽内全体の溶質ガス濃度分布を何らかの方法で求める必要がある。
【０００８】
本発明の目的は、従来は定量的に予測することが困難であった溶質ガス濃度の分布を計算評価し、その評価結果によって運転条件を決定して制御する通気攪拌槽の運転制御方法を提供することにある。
【０００９】
【課題を解決するための手段】
本発明は、通気攪拌培養槽内の溶質ガス濃度分布を最適に制御する方法であって、前記槽内の１つ以上の点で溶存酸素濃度を計測し、前記槽のガス移動容量係数（ｋｌａ）の分布を下記の式によって求め、前記分布を用いて溶存酸素濃度分布を演算し、前記計測された溶存酸素濃度を用いて前記溶存酸素濃度分布を補正し、前記補正された溶存酸素濃度分布のうち培養条件から外れている領域があるかどうかを判断し、培養条件から外れている領域がある場合は溶存酸素濃度分布があらかじめ定められた値になるまで前記培養槽の攪拌回転数、通気量、通気ガス中の酸素分圧あるいは内圧を調整することを特徴とする。
ｋｌａ=ｆ（Ｓｃ，ε，ν，α，Ｒｂ）
ただし、ｋｌａ：ガス移動容量係数（１／ｈ）、ｋｌ：液境膜移動係数（ｍ／ｓ）、a ：単位体積当たりの気液接触面積（１／ｍ）、Ｓｃ：シュミット数 （＝ν／Ｄｇ）、ε：乱流エネルギー散逸速度（ｍ／ｓ３）、ν：動粘度（ｍ２／ｓ）、Ｄｇ：液相中のガスの拡散係数（ｍ２／ｓ）、α：ガスホールドアップ、Ｒｂ：気泡半径（ｍ）である。
【００１２】
∂Ｃ／∂ｔ＝−ｄｉｖ（Ｃ・（Ｕｘ、Ｕｙ、Ｕｚ））＋Ｄｏ２∇２ Ｃ＋ｋｌａ（Ｃ＊−Ｃ）＋Ｇ
ただし、
Ｃ：液中の溶存酸素濃度（ｍｏｌ／ｍ３）
（Ｕｘ、Ｕｙ、Ｕｚ）：流速ベクトル（ｍ／ｓ）
Ｄｏ２：液相中の酸素の拡散係数（ｍ２／ｓ）
ｋｌａ：ガス移動容量係数（１／ｓ）
ｃ＊：飽和酸素濃度（ｍｏｌ／ｍ３）（＝ｐ／Ｈ）
ｐ：水圧（ａｔｍ）
Ｈ：ヘンリー定数（ａｔｍ・ｍ３／ｍｏｌ）
Ｇ：液中で単位体積・単位時間当たりに細胞が消費する溶存酸素のモル数（ｍｏｌ／ｓ・ｍ３）
また、培養中の培養槽における複数の点での溶存酸素の測定データを用いて上記の計算方法により溶存酸素濃度分布の計算を行い、その結果を利用して培養槽の溶存酸素濃度が槽内の全ての領域で培養条件を満たしているように運転制御をおこなう。
【００１３】
【発明の実施の形態】
本発明の実施例について図１、２により説明する。図１は培養槽１０の断面を表わしている。３枚羽根の２段プロペラで撹拌し、槽底に設置されたリングスパージャー通気口４から通気を行っている。１は上段のプロペラ翼を、２は下段のプロペラ翼を示し、３は攪拌機１３（図２）によって攪拌するプロペラ１、２の回転軸である。また、運転条件は、液容量１０ｍ３、撹拌回転数２５ｒｐｍ、通気量０．０３ＶＶＭである。本実施例では槽の構造及び運転条件だけから溶存酸素濃度分布を求めているので、どのような通気攪拌槽でも構造と運転条件を決定すれば、溶存酸素濃度分布を計算する事が可能である。例えばガス移動容量係数（ｋｌａ）の分布は、次式によって計算することができる。
【００１４】
kla=kl・a
kl=0.301So2-1/2(εν)1/4
a=3α/Rb
ただし、
kl:液境膜移動係数(m/s)
So2:酸素のシュミット数 (=ν/Do2)
ε:乱流エネルギー散逸速度(m2/s3)
ν:動粘度(m2/s)
Do2:液相中の酸素の拡散係数(m2/s)
a:単位体積当たりの気液接触面積(1/m)
α:ガスホールドアップ
Rb:気泡半径(m)
上記で求めたｋｌａの分布を用いて、次の式を有限差分法により計算することで、溶存酸素濃度分布を求めることができる。
【００１５】
∂C/∂t = -div(C・(Ux,Uy,Uz))+DO2∇2C+kla(C*-C)+G
ただし、
C  :液中の溶存酸素濃度(mol/m3)
(Ux,Uy,Uz)  :流速ベクトル(m/s)
kla :ガス移動容量係数(1/s)
C*  :飽和酸素濃度(mol/m3)  (=p/H)
p   :水圧(atm)
H   :ヘンリー定数(atm・m3/mol)
G :液中で単位体積・単位時間当たりに細胞が消費する溶存酸素のモル数(mol/s・m3)
次に酸素濃度分布の計算を利用して制御を行う培養槽の運転制御の実施例について述べる。通気培養槽に取り付けた１つ以上の溶存酸素計１２で溶存酸素を測定し信号を電子計算機１６に送る。電子計算機１６では、まず通気攪拌槽の構造及び運転条件を境界条件及び方程式の係数として流れの数値解析をした結果からkla分布を求め、そのkla分布を用いて溶存酸素濃度分布を計算する。その後、溶存酸素計１２から信号として送られた溶存酸素濃度で補正を行う。溶存酸素濃度の値が培養条件からはずれた領域が存在した場合、電子計算機は攪拌回転数、通気量、通気ガス中の酸素分圧、及び内圧のいずれかもしくは全部を変更して溶存酸素濃度分布を計算し直す。溶存酸素の最大値及び最小値が培養条件に適合する値が得られるまでパラメータの変更及び計算を繰り返した後、攪拌機１３、酸素分離器１５、圧力調整器１４を操作して計算で求めた攪拌回転数、通気量、酸素分圧及び内圧に実際の培養槽の運転条件を変更する。この操作を続けることで溶存酸素濃度分布が最適な運転条件で培養を行うことができる。
【００１６】
計算機１６におけるこれらの処理フローの例を図3に示す。ステップ３０では溶存酸素計１２からの信号を取り込む。ステップ３２ではガス移動容量係数ｋｌａを計算し、そしてステップ３０で取り込んだ溶存酸素により補正する。この演算された溶存酸素濃度分布の値が培養条件から外れた領域があるかどうかをステップ３４で判定し、外れている領域がある場合は、ステップ３６で攪拌回転数、通気量、通気ガス中の酸素分圧、および内圧のいずれか、または全部を変更して溶存酸素濃度分布を計算し直す。パラメータ変更計算を繰り返したあとで、実際の培養槽の運転条件をステップ３８で変更する。ステップ４０では溶存酸素が所定の範囲内にあるかどうかを判断し、満足していれば処理を終了する。満足していない場合はステップ３２以降の処理を繰り返す。このように予め定められた範囲に入るように制御することによって、培養成生物の品質の向上にもつながる。
【００１７】
また、溶存濃度が所定の範囲内にあるようにするガスは、酸素に限定されるものではなく、通気撹拌培養槽内で消費され、または生成される結果として、気泡から液中に溶け込むあるいは液中から気泡内に追い出される溶質ガスであればどのようなものでも適用できる。具体的には図３と同様な計算を、高濃度になった場合は菌の活動を阻害する二酸化炭素(菌の活動に伴い菌から液中に放出される)について行うことで、溶存二酸化炭素濃度分布を求めること及び溶存二酸化炭素濃度を無害な濃度に保つよう培養槽を制御することも可能である。
【００１８】
物質の保存則を表す微分方程式を離散化して溶質ガス濃度分布を数値解析する方法であってもよい。例えば、(溶質ガス濃度の時間変化)=(流れによる溶質ガスの流入流出量)+(拡散による溶質ガスの流入流出量)+(気液間のガス移動量)+(溶質ガスが微生物や細胞の活動又は化学反応により消費又は生成される量)、関係を利用すればよい。
【００１９】
【発明の効果】
本発明によれば、従来不可能であった通気撹拌槽内の溶質ガス濃度分布を構造及び運転条件だけから計算し、最適な値に制御することができる。したがって、槽構造の変更や新規設計を行った場合に、実際の槽を製作し溶質ガス濃度分布を測定することなく計算のみで予測することが出来る。
【００２０】
また、本発明によれば通気撹拌槽内の全ての点における溶質ガス濃度の最大値及び最小値が培養条件に収まるように、運転条件を制御することが可能であるため、培養生産物の品質維持に効果がある。
【図面の簡単な説明】
【図１】本発明の一実施例の数値解析による２段プロペラ翼を使用した通気撹拌槽の断面と溶存酸素濃度分布を示す図である。
【図２】本発明による通気撹拌培養槽の運転制御の一実施例を示すブロック構成図である。
【図３】溶存酸素濃度分布計算を制御用電子計算機で行なわせる場合の処理フロー図の一例を示す図である。
【符号の説明】
1…上プロペラ翼、2…下プロペラ翼、３…軸、４…リングスパージャー通気口、1０…培養槽、１2…溶存酸素計(一つ以上)、１３…攪拌機、１４…圧力調整器、１５…酸素分離器、１６…制御用電子計算機
【図１】
0004103288.tif 000002
【図２】
0004103288.tif 000003
【図３】
0004103288.tif 000004