(57)【要約】  本発明は、開放プレート型反応器内の化学反応を最適化する方法に関するものであり、前記反応に対する反応器の動的モデルを確立すること、積分ソフトウェアツールを前記動的モデルに適用して、その方程式を解くこと、達成されるべき目標、考慮されるべき制約条件、ならびに変数反応器および反応データの評価から一組の反応器パラメータを決定して、最適化すること、最適化されたパラメータセットに適合する反応器を構築すること、その動作を検証するためにそうして構築された反応器内で反応が実行されるときに物理的パラメータを測定することからなる。
2008529757.tif 000002
【特許請求の範囲】
【請求項１】
反応器内の化学反応を最適化する方法であって、
・熱交換ユーティリティ流体を循環させるため２つの熱交換側室（２７、２８）の間に形成された反応室（２４）を有する少なくとも１つのブロックを板と板の間に定める積層板（２０、２１、２２、２３）と、前記反応室（２４）に１つまたは複数の試薬の連続流を供給するための手段と、前記２つの側室（２７、２８）にユーティリティ流体の連続流を供給するための手段とにより構成される開放型反応器を使用することと、
・前記反応室（２４）および前記反応室（２４）とユーティリティ流体を循環させるための前記側室（２７、２８）との間の熱交換のモデルに基づいて与えられた反応に対する前記反応器の動的モデル（３６）であって、質量平衡およびエネルギー平衡方程式、および制約方程式を含む前記動的モデル（３６）を確立することと、
・積分ソフトウェアツール（３７）を前記動的モデル（３６）に適用して、上記方程式を解くことと、
・達成すべき目標条件（４７）、従うべき制約条件（４８）、および前記反応器（４９）および／または前記反応（５３）に関する変数データの評価に基づいて前記反応器の一組の寸法および／または運転パラメータを決定し、最適化することと、
・前記一組の最適化されたパラメータに応じて前記型のプレート型反応器を組み立てることと、
・その運転を検証するために、前記組み立てられた反応器内で前記反応を行わせるときに温度および圧力などの物理的パラメータの測定を実行することからなることを特徴とする方法。
【請求項２】
前記反応器を組み立てた後、
・前記組み立てられた反応器に基づいて前記反応器の新しい動的モデルを確立することと、
・前記積分ソフトウェアツール（３７）をこの新しい動的モデルに適用して、上記方程式を解くことと、
・再び前記反応器の前記一組の寸法および／または運転パラメータを最適化することと、
・必要ならば、前記新しい最適化された一組のパラメータに応じて、前記組み立てられた反応器を修正することとを含むことを特徴とする請求項１に記載の方法。
【請求項３】
前記反応室（２４）の前記モデルは、前記反応室をそれぞれ完全に攪拌された液状媒体を収めた一連の連続的個別セルに分けることを含むことを特徴とする請求項１または請求項２に記載の方法。
【請求項４】
それぞれの個別セル内の前記流体媒体は、均質であり、その状態およびその変化は、質量平衡、エネルギー平衡、圧力平衡、および容積制約方程式により定められることを特徴とする請求項３に記載の方法。
【請求項５】
前記反応室（２４）のそれぞれの個別セルに収められた流体媒体は、連続相および分散相を含む二相媒体であることを特徴とする請求項３に記載の方法。
【請求項６】
前記二相媒体の状態および変化は、前記二相媒体の密度、比熱、熱伝導率、および粘度を推定する方程式、および質量平衡、エネルギー平衡、圧力平衡、および容積制約方程式からそれぞれの個別セルにおいて定義されることを特徴とする請求項５に記載の方法。
【請求項７】
前記熱交換モデルは、前記側室（２７、２８）のそれぞれを、完全に攪拌された流体媒体をそれぞれ含む一連の個別セルに分けることと、それぞれの側室（２７、２８）内の個別セルの個数は前記反応室（２４）内の個別セルの個数に等しいことと、前記反応室（２４）の前記対応する個別セル内の前記反応流体の流れ方向に関してそれぞれの個別セル内の前記ユーティリティ流体の流れ方向を定義することとを含むことを特徴とする請求項３から６のいずれか一項に記載の方法。
【請求項８】
前記熱交換モデルは、さらに、前記反応流体、前記熱交換側室（２７、２８）を定めるプレート（２０、２１、２２、２３）、および前記側室（２７、２８）内に循環する前記ユーティリティ流体に関係する熱平衡方程式、ならびに前記側室（２７、２８）のさまざまな上記の個別セル内の前記ユーティリティ流体に対する質量平衡方程式も含むことを特徴とする請求項７に記載の方法。
【請求項９】
前記反応器の前記動的モデルは、前記反応器の入口および出口ゾーン、それぞれの反応室（２４）の入口および出口ゾーン、ならびに前記反応器のさまざまなブロック間の遷移ゾーンのモデルも含むことを特徴とする請求項１から８のいずれかに記載の方法。
【請求項１０】
完全に攪拌された内容物を有する個別セルの個数を定義するために、前記反応器のさまざまな部分における前記流体の滞留時間分布の実験研究に基づいて、前記反応器の流体力学的挙動をモデル化することを特徴とする請求項３から９のいずれか一項に記載の方法。
【請求項１１】
前記反応室（２４）のそれぞれの個別セル内の反応速度、前記反応の成分の生成速度、および熱発生率に対する方程式に基づく前記反応のモデルをさらに含むことを特徴とする請求項３から１０のいずれか一項に記載の方法。
【請求項１２】
前記二相媒体の物理的特性に基づき、及び前記分散相の液滴のサイズに基づき、前記反応媒体の前記２つの相の間の二相反応媒体中の物質移動をモデル化することを含むことを特徴とする請求項１から１１のいずれかに記載の方法。
【請求項１３】
前記側室（２７、２８）を定める前記プレートと前記反応室（２４）および前記側室（２７、２８）内で流れる流体との間の熱伝達係数の定義に基づく、及び前記反応流体および前記ユーティリティ流体の境膜係数の推定値に基づく、熱伝達のモデルをさらに含むことを特徴とする請求項１から１２のいずれかに記載の方法。
【請求項１４】
実験測定結果に基づく前記反応器内の反応流体圧力降下のモデルをさらに含むことを特徴とする請求項１から１３のいずれかに記載の方法。
【請求項１５】
前記反応器の構成要素ならびに前記反応およびユーティリティ流体の物理的特性の推定値をさらに含むことを特徴とする請求項１から１４のいずれかに記載の方法。
【請求項１６】
一連の個別セルの個数、前記反応室（２４）への試薬供給点の個数、前記反応媒体に関する前記ユーティリティ流体の前記流れ方向、前記ユーティリティ流体の性質および分布、前記反応室（２４）の容積、および／またはユーティリティ流体を循環させる前記側室（２７、２８）の容積などのパラメータを調節することにより、与えられた目標条件（４７）および制約条件（４８）について前記反応器の寸法決定を最適化することであることを特徴とする請求項３から１５のいずれか一項に記載の方法。
【請求項１７】
前記反応流体および前記ユーティリティ流体の温度、圧力、組成、および／または流量を含む動作パラメータを調節することにより、与えられた目標条件（４７）および制約条件（４８）について前記反応器の動作を最適化することであることを特徴とする請求項１から１６のいずれかに記載の方法。
【請求項１８】
複数の流体入口ゾーンを備える反応器または反応室（２４）について、前記方法は、これらのさまざまなゾーン内に注入される前記流体の前記流量、前記温度、前記圧力、および／または前記組成を調節することであることを特徴とする請求項１７に記載の方法。
【請求項１９】
前記二相反応媒体は、擬似均質媒体であると考えられ、前記２つの相の前記物理的特性は、前記二相媒体の前記物理的特性を前記相のそれぞれの特性およびさらにはその割合に関連付ける関連付け関係により決定されることを特徴とする請求項５から１８のいずれか一項に記載の方法。
【請求項２０】
前記分散相の液滴のサイズ、存在する前記相の熱力学的特性、および／または前記流れ条件の関数として、前記２つの相の間の平衡状態に対する前記相混合規則に適用することであることを特徴とする請求項１９に記載の方法。
【請求項２１】
前記反応室（２４）の前記個別セル内の前記分散相の前記液滴の前記サイズの変化に追随することであることを特徴とする請求項１９または請求項２０に記載の方法。
【請求項２２】
さらに、安全性条件および／または環境制約条件を満たす動作パラメータに対する限界値を定義することでもあることを特徴とする請求項１６から２１のいずれか一項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、１つの反応室と熱交換ユーティリティ流体を循環させるための２つの側室を有する、少なくとも１つのブロックを板と板の間に定める積層板により構成されるプレート型開放反応器内の化学反応を最適化する方法に関するものである。
【背景技術】
【０００２】
現在のところ、運転を最適化し、一般に化学反応の収量を改善する試みがなされているが、これらの反応は、通常、所定の量の試薬が入れられ、混ぜ合わされ、その後、反応を促進するように最適化された環境条件（圧力、温度、ｐＨなど）にさらされる容器を備えるバッチ型反応器、つまり、不連続に動作する反応器内で行われる。
【０００３】
このような反応器の大きな欠点は、反応により発生する熱を除去するのが困難である点であり、これが、反応を制御し、最適化する際の制約となっている。
【発明の開示】
【発明が解決しようとする課題】
【０００４】
本発明の具体的な目的は、開放型の反応器を最適化することにより化学反応を制御し、最適化する、しかも不連続反応器またはバッチ型反応器の上述の欠点を持たない、単純で、効果的な解決手段を提示することである。
【０００５】
この目的を達成するために、本発明では、反応器内の化学反応を最適化する方法を提示するが、この方法は、
【０００６】
・熱交換ユーティリティ流体を循環させるため２つの熱交換側室の間に形成された反応室を有する少なくとも１つのブロックを板と板の間に定める積層板と、反応室に１つまたは複数の試薬の連続流を供給するための手段と、２つの側室にユーティリティ流体の連続流を供給するための手段とにより構成される開放型反応器を使用することと、
・反応室および前記反応室とユーティリティ流体を循環させるための側室との間の熱交換のモデルに基づいて与えられた反応に対する反応器の動的モデルであって、質量平衡およびエネルギー平衡方程式、および制約方程式を含む動的モデルを確立することと、
・積分ソフトウェアツールを動的モデルに適用して、上記方程式を解くことと、
・達成すべき目標条件、従うべき制約条件、および反応器および／または反応に関する変数データの評価に基づいて反応器の一組の寸法および／または運転パラメータを決定し、最適化することと、
・その一組の最適化されたパラメータに応じて前記型のプレート型反応器を組み立てることと、
・その運転を検証するために、組み立てられた反応器内で前記反応を行わせるときに温度および圧力などの物理的パラメータの測定を実行することからなることを特徴とする。
【０００７】
２つのユーティリティ流体循環室の間に少なくとも１つの反応室を定める積層板を備える連続または開放型反応器を使用すると、反応室のいずれかの側で（複数の）循環するユーティリティ流体の温度および流量を適宜選択し、さらに反応室内の反応流体の流れ方向に関してこれらの流体の流れ方向（並流、向流、または交差流）を選択することにより、バッチ型の不連続または密閉型反応器内で観察される熱放散の問題を解決することが可能である。
【０００８】
また、反応器の寸法または最適な構成を最適化し、化学反応に対し制御パラメータを用意し、反応器内の温度分布を正確に識別することにより、反応器の運転条件または動作条件を最適化することが可能である。
【０００９】
本発明と不連続反応器内の反応を最適化する従来技術の方法との本質的な違いとなって現れる、反応器内の温度の変動の仕方を正確に決定することは、問題の反応について反応器の動的モデリングを作成することにより、また反応器の一組の寸法および／または運転パラメータを達成すべき目標条件および従うべき制約条件に応じて最適化することにより可能である。
【００１０】
そこで、前記型のプレート型反応器は、あらかじめ最適化されているその一組の寸法および／または運転パラメータに応じて組み立てられ、この反応器内で化学反応が実行され、化学反応最適化を検証するため物理的パラメータの実験測定が行われる。
【００１１】
本発明の方法は、さらに、反応器を組み立てた後、
・組み立てられた反応器に基づいて反応器の新しい動的モデルを確立することと、
・積分ソフトウェアツールをこの新しい動的モデルに適用して、上記方程式を解くことと、
・再び反応器の前記一組の寸法および／または運転パラメータを最適化することと、
・必要ならば、その新しい最適化された一組のパラメータに応じて、前記組み立てられた反応器を修正することとであってよい。
【００１２】
反応器の動的モデルは、反応器および化学反応の特性を示す、質量およびエネルギー平衡方程式および制約方程式を含む。
【００１３】
本発明の他の特徴によれば、反応室のモデルは、反応室をそれぞれ完全に攪拌された液状媒体を収めた一連の連続的個別セルにさらに分けることを含む。ユーティリティ流体循環室は、反応室内の個別セルの個数に等しい数の個別セルにさらに分けられることで同様にモデル化される。
【００１４】
熱交換のモデリングでも、反応室の前記個別セルのいずれかの側に置かれている側室に対応する個別セル内の（複数の）ユーティリティ流体の流れ方向（並流、向流、交差流）に関して反応室のそれぞれの個別セル内の流体の流れ方向を考慮する。このモデルは、さらに、反応流体、熱交換側室を定めるプレート、および側室内に循環するユーティリティ流体に関係する熱平衡方程式、ならびに反応媒体およびさまざまな室の個別セル内のユーティリティ流体に対する質量平衡方程式も含む。
【００１５】
それぞれの反応室は、反応室の個別セルに関連付けられ、また例えば、試薬供給手段と反応室供給点との間の接続要素、および反応室出口点と反応の生成物を回収するための手段との間の接続要素である、入口ゾーンと出口ゾーンに接続される。２つの反応室が直列につなぎ合わされている場合、第１の反応室の出口要素は第２の反応室の入口要素と一致する。
【００１６】
都合のよいことに、反応器の動的モデルは、反応器の入口および出口ゾーン、それぞれの反応室の入口および出口ゾーン、ならびに反応器のさまざまなブロック間の遷移ゾーンのモデルを含む。
【００１７】
この動的モデルでは、さらに、反応室のそれぞれの個別セルに対し供給点を設けることが可能である。
【００１８】
反応器の動的モデルを構成する連立方程式に適用される積分ソフトウェアツールは、これらの方程式を解き、研究対象の化学反応に対する、「反応状態」と呼ばれる、反応器の挙動をシミュレートするために使用される。
【００１９】
与えられた目標条件および制約条件について、反応器の寸法パラメータは、反応器の可変特性を調節することにより決定され、最適化されるが、前記可変特性は、反応室内の供給点の個数、反応媒体に関する（複数の）ユーティリティ流体の（複数の）流れ方向、（複数の）ユーティリティ流体の性質および分布、反応室の総容積、および／またはユーティリティ流体循環に対する側室の容積を含み、与えられた目標条件をおよび制約条件について、両方とも反応器の運転パラメータは、反応媒体および（複数の）ユーティリティ流体の温度、圧力、組成、および／または流量を含む可変反応データまたは動作パラメータを調節することにより決定され、最適化される。
【００２０】
この方法は、例えば、１つまたは複数の試薬をさまざまな個別セルに供給する速度を調節して反応の収量を高め、また安全および／または環境制約条件を満たすように動作パラメータに対する限界値を定義することである。
【００２１】
均質な媒体中の化学反応では、反応器の動的モデルは、質量平衡、エネルギー平衡、圧力平衡、および容積制約条件に対する方程式からなる反応室の動的モデルを含む。これらの方程式は、それぞれの個別セル内の流体媒体の状態および変化量を定義するために使用される。
【００２２】
これを行うために、完全に攪拌された内容物を有する反応室の個別セルにおいて、特性量の大きさ（温度、圧力、流量、組成など）、物理的特性（密度、粘度など）、および物理化学的現象（反応、物質移動など）は、それぞれのセル内で一様であり、ユーティリティ流体の物理的特性は、与えられた温度について一様で、一定であると考えられる。
【００２３】
二相媒体（液体−液体または液体−気体）における化学反応では、反応媒体は、連続相と分散相の混合であり、２つの混和液体の混合である反応媒体では、反応媒体の状態および反応媒体中の変量は、反応媒体の密度、比熱、熱伝導率、および粘度を推定する方程式から、また質量平衡、エネルギー平衡、圧力平衡、および容積制約条件に対する方程式から、それぞれの個別セルにおいて定義される。
【００２４】
完全に攪拌されている内容物を有する反応室の個別セル内において、特性量の大きさ（温度、圧力、流量、組成）、物理的特性（密度、粘度など）、および物理化学的現象（反応、物質移動など）は、それぞれのセル内で一様であり、相の温度および圧力は、同じであり、反応器媒体は、擬似均質媒体であるとみなされ、２つの相の物理的特性は、媒体の物理的特性を相のそれぞれの特性、さらにその割合に関連付けることを可能にする関連付けの関係により決定されると考えられ、その際に相の分布はそれぞれの個別セル内で一様であると仮定される。
【００２５】
また、二相平衡規則を、分散相の液滴のサイズ、存在する相の熱力学的特性、および／または流れ条件に応じて相混合に適用することも可能である。
【００２６】
ザウター径により、分散相の液滴のサイズを決定することが可能である。分散相中の液滴のサイズの変化量は、さらに、反応室のさまざまな個別セル内で追跡することもできる。
【００２７】
反応器の動的モデルに追加し、反応器の動的モデルを精密化できる他のモデルは、
・反応器のさまざまな部分における流体の滞留時間分布の実験研究に基づいて、反応器の流体力学的挙動をモデル化し、完全に攪拌された内容物を有する個別セルの個数を定義できるようにすること、
・（複数の）反応室のそれぞれの個別セル内の反応速度、反応成分の生成速度、および発熱速度の方程式に基づく反応モデル、
・二相媒体の物理的特性に基づく、また分散相の液滴のサイズに基づく、反応媒体の２つの相の間の二相媒体中の反応の物質移動モデル、
・側室を定めるプレートと反応室および側室内で循環する流体との間の熱伝達係数の定義に基づく、また反応流体およびユーティリティ流体に対する推定境膜係数に基づく、熱伝達モデル、および
・実験測定結果に基づく、反応器内の反応流体圧力降下のモデルなどである。
【００２８】
この方法は、さらに、反応器の構成要素ならびに反応およびユーティリティ流体の物理的特性を推定することも含みうる。
【００２９】
本発明の方法は、与えられた動作条件に対するプレート型反応器の挙動を予測する利点を有し、またこれにより、与えられた反応または一組の反応について反応器の最適な構成または寸法決定および最適な動作条件または運転を定義することも可能である。この方法では、さらに、時間の経過による前記挙動の変化を記述することが可能であり、反応器の制御の最適化、安全性の研究、および非生産過渡期の管理に関してさまざまな方法によりこのような変化を表すことが可能である。
【００３０】
本発明の他の利点および特徴は、非制限的な実施例を使い、付属の図面を参照しつつなされている以下の説明を読むと理解されるであろう。
【発明を実施するための最良の形態】
【００３１】
最初に図１を参照すると、これは、互いの間に平行な室を定義するプレート２０、２１、２２、２３からなる積層板を備える開放または連続反応器の図である。プレート２１および２２は、「遷移」プレートと呼ばれ、これらは、間に、１つまたは複数の他の試薬に対する試薬供給点２５ａおよび１つまたは複数の供給点２５ｂ、および反応の（複数の）生成物に対する出口点２６を有する反応室２４を定める。プレート２０および２３は、断熱プレートであり、プレート２１および２２と連携し、熱交換ユーティリティ流体を循環させるための左側室２７および右側室２８を定め、それぞれの側室２７、２８はユーティリティ流体用の入口２９および出口３０を備える。
【００３２】
左側室２７内で循環するユーティリティ流体は、右側室２８内で循環する流体と同じでも、異なっていてもよく、通常使用されるユーティリティ流体は、水またはエチレングリコールである。
【００３３】
プレート２０、２１、２２、２３からなる積層板は、反応ブロックを形成し、プレート２０および２３は、それを外部環境から、または直列に接続されている１つまたは複数の反応ブロックから絶縁するために使用される。
【００３４】
室２４、２７、および２８は、一般に、前記室内に流体流の通り道を定める手段を取り付けられる。例えば、反応室２４は、図１および２に示されているように、反応媒体が反応室２９内部でジグザグ経路を辿るように制約するバッフルを形成する挿入物を備える。
【００３５】
室２４、２７、および２８は、適当な手段（図に示されていない）により連続的供給を受ける。
【００３６】
最適な動作条件の下で、プレート型反応器の総括伝熱係数は、２０００ワット／平方メートル／ケルビン（Ｗ．ｍ−２．Ｋ−１）から５０００  Ｗ．ｍ−２．Ｋ−１までの範囲とすることができる。反応室２４内の温度は、例えば、反応室２４のいずれかの側で循環する（複数の）ユーティリティ流体の温度および／または流量を適宜選択することにより、さらに反応室内の反応媒体の組成、温度、および／または流量を適宜選択することにより、および／またはそれぞれの側室内のユーティリティ流体の流れ方向を反応室内の反応媒体に関して適宜選択することにより制御される。図１に示されているように、ユーティリティ流体は、反応媒体が反応室２４内で流れるのと同じ方向に２つの側室２７および２８内で流れる場合に並流として流れ、反対方向に流れる場合に向流として流れ、一方の側室内では反応媒体と同じ方向に、他方の側室内では反対方向に、流れる場合に交差流として流れる。この実施例では、反応媒体がジグザグ経路を辿って水平方向に流れており、ユーティリティ流体が垂直方向に流れているため、反応媒体とユーティリティ流体は交差するが、並列構成など他の流れの構成も採用することが可能である。
【００３７】
本発明の方法の本質的な特徴は、与えられた反応に対する反応器の動的モデルを確立することにある。
【００３８】
このモデルは、滞留時間分布の実験調査を使って、反応室をそれぞれ完全に攪拌されている複数の個別セルにさらに分け、それにより、流れを分析し、デッドボリュームおよび好ましい通り道を見つけ出し、反応室２４内の混合物を特徴付けられるようにすることを含む。
【００３９】
滞留時間分布の実験分析により、平均滞留時間
【００４０】
【数１】
2008529757.tif 000003
【００４１】
および分散（σ２）に関して特徴付けることが可能である。これらの特徴に基づき、混合の程度および実際の反応器の流体力学的挙動を表すために必要な個別セルの個数を推定することが可能である。例えば、この数値は、式
【００４２】
【数２】
2008529757.tif 000004
【００４３】
を使用して計算される。
【００４４】
実際、この数値は、特に反応室２４の供給速度、および反応媒体の粘度に依存する。３つの反応ブロックを含むプレート型反応器に対する個別セルの個数は、例えば、２５リットル／時（Ｌ．ｈ−１）の速度で供給される水のみからなる反応媒体の場合には７０であり、５０Ｌ．ｈ−１で供給される水の場合には１００であり、８０Ｌ．ｈ−１で供給される水の場合には１３０であり、５０Ｌ．ｈ−１の速度で供給されるエチレングリコールからなる反応媒体の場合には２５である。
【００４５】
図１および２の実施例では、反応室は、１２個の個別セルを有する。これら１２個の個別セルは、反応室２４を通り供給点２５ａ、２５ｂから前記反応室２４の出口点２６まで進む反応媒体の水平方向ジグザグの流れの順序に従って１から１２までの番号が振られる。この流れは、反応室２４内に存在する上述の挿入物により強いられる。
【００４６】
１２個のセルは、それぞれの行が３個のセルからなる４つの行を構成する。試薬は、セル１内の供給点２５ａから反応室２４内に導入される。反応媒体は、水平方向に、セル１からセル２に進み、次いで、セル３に進み、セル３を出て、その下に垂直に配置されているセル４内に入り、次いで水平方向に進んで、セル５および６内に入る。反応媒体は、セル６を出て、その下に垂直に配置されているセル７内に入り、次いで、水平方向に進み、セル８および９内に入る。最後に、反応媒体は、セル９を出て、その下に垂直に配置されているセル１０内に入り、次いで、水平方向に進み、セル１１および１２内に入り、点２６を経由して反応室２４を出る。１つまたは複数の他の試薬に対する他の供給点２５ｂも、第１のセルまたは第１のセル以外の個別セル内に用意することができる。
【００４７】
一変更形態では、反応媒体は、側室２７および２８内を流れているユーティリティ流体に関して平行に反応室２４内で流れることができる。
【００４８】
同様に、側室２７および２８は、ユーティリティ流体の流れについて複数の個別セルにさらに分けられ、それぞれの個別反応器セルは、２つの個別ユーティリティ流体流セルの間に配置される。
【００４９】
図１では、左および右側室２７および２８は、個別のユーティリティ流体流セルのそれぞれの系列Ａ、Ｂ、Ｃ、Ｄ、Ｅ、Ｆ、Ｇ、Ｈ、Ｉ、Ｊ、Ｋ、Ｌ、およびＡ’、Ｂ’、Ｃ’、Ｄ’、Ｅ’、Ｆ’、Ｇ’、Ｈ’、Ｉ’、Ｊ’、Ｋ’、Ｌ’を含む。セルを表す文字のアルファベット順は、個別反応セルの個数の昇順と一致し、左側室２７では下方向で、右側室２８では上方向のユーティリティ流体の流れとは一致しない。
【００５０】
反応室２４を複数の個別セルにさらに分けることにより、それぞれの個別反応セル内を流れる反応媒体が反応室２４を定める２つの遷移プレート２１および２２を通して作用する、２つの対応する側部個別セル内を流れるユーティリティ流体の熱的な影響を決定することが可能である。例えば、熱は、プレート２１を通して個別ユーティリティ流体フローセルＥと個別反応セル５との間で、またプレート２２を通して個別ユーティリティ流体フローセルＥ’と個別反応セル５との間で、交換されうる。
【００５１】
図１に示されている実施形態では、セル５の熱的特性は、反応媒体の流体流れ方向でセル５の直前に置かれているセル、つまりセル４の熱的特性により決まる。同様に、セルＥおよびＥ’の熱的特性は、ユーティリティ流体流れ方向でそれぞれＥおよびＥ’の前に配置されるセルＢおよびＨ’によりそれぞれ決定される。
【００５２】
個別反応セルは、連続的に攪拌される反応器と考えられ、したがって、反応器および化学反応の特性および平衡状態に対する方程式を、個別反応セルの個数の関数となっている多数の方程式を使用する数学的モデルで一般化することが可能であり、この数値は、式
ｎ＝［セルの個数］×（セル１つ当たりの方程式の個数）
で与えられる。
【００５３】
数学的モデルの中の方程式の個数は、また、ユーティリティ流体フローチャンバ内で使用されるゾーンの個数の関数でもある。「ユーティリティゾーン」という用語は、専用のユーティリティ流体供給点および出口点を有するユーティリティ流体フローチャンバの一部を指し示すために使用される。したがって、３つのユーティリティゾーンを有するユーティリティ流体フローチャンバは、ユーティリティ流体用の３つの供給点および３つの出口点を備える。そこで、数学的モデル内の方程式の個数は、式
ｎ＝［セルの個数］×（セル１つ当たりの方程式の個数）＋（ゾーン１つ当たりの方程式の個数）×［ユーティリティゾーンの個数］
で与えられる。
【００５４】
考察対象のモデルに応じて、つまり、均質媒体または二相媒体について、方程式の個数は、同じでなく、モデルの方程式の総数は以下の式で定義される。
・均質モデル：
ｎ＝［（ｎｃｅｌｌ／ｂｌｏｃｋ＋３）×ｎｂｌｏｃｋ＋１］×（１６＋ｎｃｏｎｓ）＋３×ｎｚｏｎｅ／ｂｌｏｃｋ×ｎｂｌｏｃｋ
・二相モデル：
ｎ＝［（ｎｃｅｌｌ／ｂｌｏｃｋ＋３）×ｎｂｌｏｃｋ＋１］×（２１＋２×ｎｃｏｎｓ）＋３×ｎｚｏｎｅ／ｂｌｏｃｋ×ｎｂｌｏｃｋ
ただし、ｎｃｏｎｓは、化学反応の成分の個数であり、ｎｂｌｏｃｋは、反応ブロックの個数であり、ｎｚｏｎｅ／ｂｌｏｃｋは、反応ブロック１つ当たりのユーティリティゾーンの個数であり、ｎｃｅｌｌ／ｂｌｏｃｋは、反応ブロック１つ当たりの個別セルの個数である。
均質または単相媒体での化学反応については、特性方程式は、例えば、以下のとおりである。
・反応室のそれぞれの個別セルについては、総質量平衡、構成要素質量平衡、エネルギー平衡、圧力平衡、容積モデル、モルエンタルピーモデル、容積制約条件、モル容積モデル、熱的環境の熱平衡、遷移プレートの熱平衡、ユーティリティ流体の熱平衡、断熱プレートの熱平衡、およびユーティリティ流体の質量平衡、および
・ユーティリティ流体流れ内のそれぞれのユーティリティゾーンについては、ユーティリティ流体温度制約条件およびユーティリティ流体流量制約条件。
【００５５】
二相または不均質媒体での化学反応については、特性方程式は、例えば、以下のとおりである。
・反応室のそれぞれの個別セルについては、連続相および分散相の総質量平衡、連続相および分散相の構成要素質量平衡、エネルギー平衡、圧力平衡、連続相および分散相の容積モデル、連続相および分散相のモルエンタルピーモデル、連続相および分散相のモル容積モデル、連続相および分散相の容積制約条件、熱的環境の熱平衡、遷移プレートの熱平衡、ユーティリティ流体の熱平衡、断熱プレートの熱平衡、およびユーティリティ流体の質量平衡、および
・ユーティリティ流体流れ内のそれぞれのユーティリティゾーンについては、ユーティリティ流体温度制約条件およびユーティリティ流体流量制約条件。
【００５６】
本発明の動的反応器モデルは、反応室の動的モデルおよび熱交換の動的モデルに基づいて作成され、それぞれの動的モデルは矩形により図３に線図として表されている個別反応セル５に関係する方程式を含む。
【００５７】
矢印３２および３３は、それぞれ、反応室２４の個別セル５内の反応媒体および側室２７および２８の個別側部セルＥおよびＥ’内のユーティリティ流体の流れ方向を表すが、この実施例の２つの側室２７および２８内のユーティリティ流体流れ方向は、図１に示されているのとは反対に、反応室２４内の反応媒体の流れに関して並流となっている。
【００５８】
インデックスｐおよびｕは、それぞれ、反応媒体およびユーティリティ流体を定義するために使用される。Ｆｐは、個別セル内の反応媒体のモル流量を表し、モル／秒（ｍｏｌ．ｓ−１）で示され、ＴｐおよびＴｕは、それぞれ、°Ｋで表される反応媒体、および個別セル内のユーティリティ流体の温度であり、Ｈｐは、ジュール／モル（Ｊ．ｍｏｌ−１）で表された個別セル内の反応媒体のモルエンタルピーであり、Ｖｐは、立方メートル（ｍ３）で表された個別反応室の容積であり、ｕｐは、個別セル内の反応媒体中のモルの数であり、ｘｐは、個別セル内の試薬または成分（ｉ）のモル分率であり、上付き文字ｋ−１およびｋは個別反応セルの番号（図１および２に示されている反応室２４の実施形態では１から１２まで）または個別ユーティリティ流体フローセルの文字（図１の左および右側室２７および２８についてＡからＬおよびＡ’からＬ’まで）を示す。
【００５９】
反応室２４の個別セルの内容物が完全に攪拌されていると仮定されている場合、アレニウスの式を適用できる均質単相液体媒体中で化学反応が生じ、反応室の動的モデルは、例えば、
・個別セル内の反応媒体に対する総質量平衡（ｍｏｌ．ｓ−１）：
【００６０】
【数３】
2008529757.tif 000005
【００６１】
・個別セル内の反応媒体の成分（ｉ）に対する構成要素質量平衡（ｍｏｌ．ｓ−１）：
【００６２】
【数４】
2008529757.tif 000006
【００６３】
・個別セル内の反応媒体に対する圧力平衡方程式（パスカル（Ｐａ））：
【００６４】
【数５】
2008529757.tif 000007
【００６５】
・個別セル内の反応媒体に対する容積モデル方程式（ｍ３）：
【００６６】
【数６】
2008529757.tif 000008
【００６７】
・個別セル内の反応媒体のモルエンタルピー方程式（Ｊ．ｍｏｌ−１）：
【００６８】
【数７】
2008529757.tif 000009
【００６９】
・個別セル内の反応媒体に対する容積制約条件またはモル流量方程式（ｍ３またはｍｏｌ．ｓ−１）：
【００７０】
【数８】
2008529757.tif 000010
【００７１】
を含み、この段階では、個別反応セルはすべて、注入されており、そこで、セルの注入が済むと、
【００７２】
【数９】
2008529757.tif 000011
【００７３】
となるが、
ただし、Ｆｆは、単位ｍｏｌ．ｓ−１で表される供給速度であり、Δｎは、モル／立方メートル／秒（ｍｏｌ．ｍ−３．ｓ−１）を単位とする生成速度を表し、αは、化学量論係数であり、ｒは、ｍｏｌ．ｍ−３．ｓ−１を単位とする反応速度であり、ｋ°は、前指数因子であり、Ｅａは、Ｊ．ｍｏｌ−１を単位とする活性化エネルギーであり、Ｒは、理想気体定数であり、Ｃは、ｍｏｌ．ｍ−３を単位とする試薬濃度であり、Ｖは、ｍ３を単位とする容積であり、ｍＶおよびｍＨは、それぞれ、モジュールのモデルおよびエンタルピーのモデルであり（ｍ３およびＪ．ｍｏｌ−１）、Ｖｍｌは、立方メートル／モル（ｍ３．ｍｏｌ−１）を単位とするモル容積であり、上付き文字「ｃｅｌｌ」は、物理セルを表し、上付き文字βは、反応次数であり、下付き文字ｊは、反応の番号であり、下付き文字ｉは、考察されている成分であり、ｔは、ｓ（秒）を単位とする時間であり、Ｔは、°Ｋを単位とする温度であり、Ｐは、圧力であり、ΔＰは、式
【００７４】
【数１０】
2008529757.tif 000012
【００７５】
によって定義されるパスカル（Ｐａ）を単位とする圧力降下であるが、
ただし、ファニングまたは摩擦係数ｆは、層流（レイノルズ数Ｒｅ＜２２００）である個別反応セルを通る反応媒体の流れに対しｆ＝５．０４６４  Ｒｅ−０．５３２８であるか、または乱流（Ｒｅ  ＞  ２２００）に対しｆ＝２．１７３４７  Ｒｅ−０．４２３１６であり、
Ｌは、メートル（ｍ）を単位とする反応室２４を通る反応媒体が辿る経路の長さであり、ｐは、キログラム／立方メートル（ｋｇ．ｍ−３）を単位とする反応媒体の密度であり、ｖは、メートル／秒（ｍ．ｓ−１）を単位とする反応媒体の速度であり、ｄｈは、メートル（ｍ）を単位とする個別反応セルの相当水力直径である。
【００７６】
プレート型反応器が、一連の連続攪拌反応器で構成され、ユーティリティ流体の物理的特性が与えられた温度に関し一定で、均質であると仮定された場合、ユーティリティ流体の動的モデルは、例えば、式として、
・ユーティリティ流体流れに対する連続方程式（立方メートル／秒（ｍ３．ｓ−１））
【００７７】
【数１１】
2008529757.tif 000013
【００７８】
・ユーティリティ流体の温度に対する連続方程式（ケルビン（Ｋ））
【００７９】
【数１２】
2008529757.tif 000014
【００８０】
を含むが、
ただし、上付き文字「ｉｎ」および「ｏｕｔ」は、ユーティリティ流体フローチャンバの入口および出口をそれぞれ表し、「ｎｃｒ」は、ユーティリティ流体フローチャンバ内のセルの行毎のセルの個数を表し（図１では、ｎｃｒは、側室２７および２８に関して３に等しい）、「ｎｃ」は、ユーティリティ流体フローチャンバ内の最後のセルを表す（図１では、側室２７において、ｎｃは、Ｌを表し、ｎｃ−１は、Ｋを表し、ｎｃ−２は、Ｊを表し、側室２８において、ｎｃは、Ａ’を表し、ｎｃ−１は、Ｂ’を表し、ｎｃ−２は、Ｃ’を表す）。
【００８１】
動的モデルは、さらに、試薬供給手段とセル１内の供給点２５ａとの間にある接続要素３４（図２）内の（複数の）試薬の挙動に対するモデル、およびセル１２からの出口点２６と試薬の生成物を回収するための手段または第２の反応室（図に示されておらず、第２の反応ブロックを構成する）内のセルへの供給点または入口との間の接続要素３５内の反応媒体の挙動に対するモデルも含む。プレート型反応器の動的モデルを作成する場合に接続要素３４および３５を考慮することが重要であるが、それは、これらの要素の全容積が、反応室２４の容積と比較して無視できず、化学反応の収量に対し一定の影響を持ちうるからである。さらに、このような接続要素内では、ユーティリティ流体との熱交換は制限されている。
【００８２】
遷移プレート２１、２２、ユーティリティ流体フローチャンバ２７、２８、断熱プレート２０、２３、および反応媒体の熱的環境（反応室２４内の挿入物の存在）などの、複数のプレート型反応器要素が熱交換に関わっている。
【００８３】
個別反応セル内の反応器媒体の熱交換モデルは、２つのプレート２１および２２により定められた矩形、および図１のプレート２０、２１、および２２、２３によりそれで定められる２つのユーティリティ流体流れ側部セルＥ、Ｅ’により図４において線図として表される個別反応セル５に適用される特性方程式を含む。矢印３２および３３は、それぞれ、反応室２４の個別セル５内の反応媒体および側室２７、２８の個別側部セルＥ、Ｅ’内のユーティリティ流体の流れ方向を表すが、２つの側室２７、２８内のユーティリティ流体の流れ方向は、図１に示されているのとは反対に、反応室２４内の反応媒体の流れに関して並流となっている。
【００８４】
個別側部セル内のユーティリティ流体の温度は、ユーティリティ流体の流れ方向の先行する個別側部セルの温度、およびまた反応ブロック内のユーティリティゾーンの存在に左右される。そのため、反応媒体に関してユーティリティ流体のさまざまな流れ方向について考えられる熱的構成すべてを表すことが可能である。
【００８５】
下付き文字ｕＬおよびｕＲは、左側室２７内および右側室２８内のユーティリティ流体をそれぞれ表すために使用され、他のシンボルは、上述のものと同じ意味を有する。
【００８６】
遷移および断熱プレート、ならびにユーティリティ流体の物理的特性は、与えられた温度において個別セルに対し一定であり、それぞれのプレートの温度は、壁の真ん中で計算されると仮定される場合、反応媒体熱交換のモデルは、式として、例えば、
・個別反応セル内の反応媒体に対するエネルギー平衡方程式（ジュール／秒（Ｊ．ｓ−１））：
【００８７】
【数１３】
2008529757.tif 000015
【００８８】
・左側室２７の側部セル内のユーティリティ流体に対する熱平衡方程式（Ｊ．ｓ−１、同じ式が右側室２８内の側部セルに使用される）：
【００８９】
【数１４】
2008529757.tif 000016
【００９０】
・個別反応セルと左側室２７の側部セルとの間に配置されているプレート２１に対する熱平衡方程式（Ｊ．ｓ−１、同じ式が個別反応セルと右側室２８の側部セルとの間に配置されているプレート２２に使用される）：
【００９１】
【数１５】
2008529757.tif 000017
【００９２】
を含むが、
ただし、ＦｐおよびＦｕは、それぞれ、ｍｏｌ．ｓ−１を単位とする反応媒体のモル流量、およびｍ３．ｓ−１を単位とするユーティリティ流体の容積流量を表し、ｅは、ｍ（メートル）を単位とする反応室２４を定めるプレート２１、２２の厚さであり、λは、ジュール／秒／メートル／ケルビン（Ｊ．ｓ−１．ｍ−１．Ｋ−１）を単位とする熱伝導率であり、ｈは、ジュール／秒／平方メートル／ケルビン（Ｊ．ｓ−１．ｍ−２．Ｋ−１）を単位とする熱伝達係数であり、Ｔは、Ｋを単位とする温度であり、Ｃｐは、ジュール／キログラム／ケルビン（Ｊ．ｋｇ−１．Ｋ−１）を単位とする媒体の比熱であり、ρは、ｋｇ．ｍ−３を単位とする媒体の密度であり、Δｑは、Ｊ．ｍ−３．ｓ−１を単位とする反応により生じる熱であり、ΔＨｒは、Ｊ．ｍｏｌ−１を単位とする反応のモル熱であり、Ａは、ｍ２を単位とする熱交換の面積であり、下付き文字ｐｕＬおよびｐｕＲは、それぞれ、反応室２４と左側室２７との間のプレート２１、および反応室２４と右側室２８との間のプレート２２を表す。
【００９３】
熱交換に対する上述のモデルは、プレート型反応器の熱慣性を含むように拡張される。熱慣性モデルは、構造、およびプレート型反応器のさまざまな要素の組成を考慮に入れるために使用される。例えば、これは、さまざまな流体の熱的環境を考慮に入れるという点で先行する方程式と異なる方程式である、
・個別反応セル内の反応媒体の環境に対するエネルギー平衡方程式（Ｊ．ｓ−１）：
【００９４】
【数１６】
2008529757.tif 000018
【００９５】
・左側室２７の側部セル内のユーティリティ流体環境に対する熱平衡方程式（Ｊ．ｓ−１、同じ式が右側室２８内の側部セルに使用される）：
【００９６】
【数１７】
2008529757.tif 000019
【００９７】
・個別反応セル内の反応媒体の環境に対する熱平衡方程式（Ｊ．ｓ−１）：
【００９８】
【数１８】
2008529757.tif 000020
【００９９】
・左側室２７の断熱プレート２０の環境に対する熱平衡方程式（Ｊ．ｓ−１、同じ式が右側室２８内の断熱プレート２３に使用される）：
【０１００】
【数１９】
2008529757.tif 000021
【０１０１】
を含むが、
ただし、熱伝達係数（単位Ｊ．ｓ−１．ｍ−２．Ｋ−１）は、式
【０１０２】
【数２０】
2008529757.tif 000022
【０１０３】
で表され、
下付き文字ｔｅｐおよびｔｅｕは、それぞれ、反応媒体の熱的環境およびユーティリティ媒体の熱的環境を表し、Ｈｆは、Ｊ．ｍｏｌ−１を単位とする供給流量のモルエンタルピーである。
【０１０４】
反応媒体とユーティリティ流体との間の熱交換は、それぞれの遷移プレート２１、２２について別々に考慮され、動的モデルにおけるこれらのプレート（例えば、ステンレス製）による熱慣性が含まれる。
【０１０５】
反応媒体の熱的環境の温度に依存する、個別反応セル内の反応媒体の環境に対する熱平衡方程式では、反応室２４の構造、および前記反応室内に存在する挿入物（例えば、ポリエーテルエーテルケトンでできているもの）による熱慣性を考慮に入れることが可能である。与えられた個別セルについて、前記熱的環境の熱伝達面積および質量は、挿入物の幾何学的特徴に基づき、また反応室２４の構造に基づき、計算される。
【０１０６】
ユーティリティ流体の熱的環境の温度に依存する、左および右側室２７および２８の個別セル内のユーティリティ流体の環境に対する熱平衡方程式では、動的モデルにおける断熱プレート２０および２３（例えば、ステンレス製）による熱慣性を考慮に入れることができる。熱的環境の特性（熱伝達面積、質量）は、これらの断熱プレートの幾何学的特徴から計算される。
【０１０７】
二相媒体において実行される化学反応では、反応器の動的モデルは、相のそれぞれについて確立された、上述の特性方程式および平衡方程式（エネルギー、熱、質量など）のすべてを含み、連続相は液体であり、分散相は液体または気体である。反応室２４内の反応媒体、側室２７および２８内のユーティリティ流体、（複数の）反応室２４の入口および出口要素３４および３５内の反応媒体の挙動、圧力降下、個別セル１から１２内の反応媒体の熱交換、および反応器の熱慣性の動的モデルは、与えられた二相媒体における反応に対するプレート型の開放反応器の動的モデルを確立するように適合される。
【０１０８】
このために、相のそれぞれにおける温度および圧力は、同じであると考えられ、反応媒体は、対応する二相媒体の物理的特性が二相媒体の物理的特性をその相のそれぞれの特性に、さらにそれぞれの割合に関連付ける関連付けの関係により決定される擬似均質媒体であると仮定され、相は、それぞれの個別セル内に均質分布すると仮定される。関連付けの関係の例は、
・二相媒体の推定密度：
【０１０９】
【数２１】
2008529757.tif 000023
【０１１０】
・二相媒体の推定比熱：
【０１１１】
【数２２】
2008529757.tif 000024
【０１１２】
・二相媒体の推定熱伝導率：
【０１１３】
【数２３】
2008529757.tif 000025
【０１１４】
・二相媒体の推定粘度：
【０１１５】
【数２４】
2008529757.tif 000026
【０１１６】
で与えられるが、
ただし、ｎｉは、モルを単位とする相のモル数を表し、Ｍｉは、キログラム／モル（ｋｇ．ｍｏｌ−１）を単位とする相のモル質量であり、ｍは、キログラム（ｋｇ）を単位とする二相媒体の全質量であり、ｉは、連続または分散相を指定する下付き文字であり、λは、Ｊ．ｍ−１．ｓ−１．Ｋ−１を単位とする媒体の熱伝導率であり、μは、パスカル秒（Ｐａ．ｓ）を単位とする媒体の粘度である。
【０１１７】
二相媒体中の反応に対する反応器の動的モデルの精度は、反応媒体が安定していること、したがって２つの相の間で平衡状態にあることに基づく。数学的モデルで使用される平衡規則は、連続相内の分散相の液滴のサイズの関数となっている。液滴のサイズおよび分布は、液滴母集団を評価するためのソフトウェアツールを使って決定される。
【０１１８】
母集団を評価するソフトウェアツールは、時間の経過とともに液滴サイズの初期分布がどのように変化するかを予測するために使用され、したがって、プレート型反応器の反応室２４のさまざまな個別反応セルにおけるザウター径を予測するために使用される。この変化量は、液滴が出現、消失する速度に依存し、この速度は、凝集および破裂の現象に直接関連する。凝集および破裂の現象は、二相媒体の物理的および輸送特性に依存する相関関係からプログラムに導入される。
【０１１９】
物質移動モデルは、分散相がザウター径に等しい平均直径を有する球状液滴の形態で存在するという仮定に依存する。２つの相の間の物質移動は、Ｗｈｉｔｍａｎの二重境膜説に基づいてモデル化される。
【０１２０】
本発明の最適化法の動的モデルでは、個別反応セル毎に物質移動特性が計算される。物質移動には相関関係があるため、物質移動特性のすべて（物質移動係数、界面面積、溶解度）は、液体−液体二相媒体の物理的特性およびザウター径に関連する。したがって、与えられた液体−液体二相媒体について、プレート型反応器の挙動を正確に表現するには、反応室のそれぞれの個別セルに対するザウター径の正しい推定値があるだけでよい。
【０１２１】
図５ａは、本発明の最適化法における本質的工程を示す図であり、プレート型反応器の動的モデル３６は、上述のように反応器に対する複数の特性方程式および平衡方程式を含む。
【０１２２】
積分ソフトウェアツール３７は、与えられた反応に対する反応器の動的モデルに基づきプレート型反応器の挙動をシミュレートするために使用される。このツールを使用することで、反応器３８の流体力学的記述、反応媒体、（複数の）ユーティリティ流体、および反応器の構成要素（挿入物、プレート、．．．）の材質の物理的特性３９、化学反応の反応熱および反応速度のモデル４０、さまざまな流体の間の熱伝達４１、および二相媒体中の反応に対する液体−液体物質移動４２または液体−気体物質移動４２を含むデータベースと連結することにより動的モデルの代数方程式と微分方程式の両方の連立方程式を解くことができる。
【０１２３】
動的モデルを組み込むために、このツールではさまざまなモデルに関わる多数のパラメータまたは入力を必要とする。これらのパラメータは、反応の個数、試薬の個数、試薬の性質、その流量、その温度、その圧力、その供給持続時間、反応の特性、化学量論係数、反応次数、速度定数に対する前指数因子、活性化エネルギー、ならびにユーティリティ流体の流量および温度などの反応に対する動作条件４３、例えば反応ブロック毎の個別セルの個数、行毎の個別セルの個数、反応ブロックの個数、反応ブロック毎のユーティリティゾーンの個数、プレート型反応器の寸法、外部媒体の温度、熱交換面積、プレートの枚数、その厚さ、その密度、その熱容量、その熱伝導率、挿入物の寸法、およびユーティリティ流体の性質を含む反応器の基本構成の記述４４、および絶対誤差の許容範囲または積分に使用可能な時間などの他の積分パラメータを含む。
【０１２４】
当業者であれば、この種類のツールおよびデータベースを知っている。データベースは、例えば、Ｂｉｐｐｈｙ（登録商標）である。積分ソフトウェアツールは、例えば、ＤＩＳＣｏであり、これは、５０００個を超える特性方程式および平衡方程式を含む動的モデルを高速に、高い精度で解くことができる。
【０１２５】
ＤＩＳＣｏ（Ｓｏｆｔｗａｒｅ  Ｃｏｍｐｏｎｅｎｔｓ社のＤｏ  Ｉｎｔｅｇｒａｔｅ）は、予測子−修正子法を使用するＧｅａｒの後退微分法に基づいて代数および微分連立方程式を積分するためのツールである。このツールは、方程式または変数を区別することなく、連立方程式を全体として処理すること、追加の方程式、したがって追加の変数を基本連立方程式に組み込み、これにより数学的モデルの自由度を高めることが可能であること、イベントを検出する自動化手順、イベント（手順の開始、力学）を探し、管理する場合に特に有用な矛盾のない初期条件を計算すること、および計算時間を大幅に短縮する疎行列を使って連立方程式を処理することが可能であることなどの多くの利点を有する。
【０１２６】
積分ソフトウェアツール３７は、与えられた化学反応に対する「反応器状態」４５と呼ばれる化学反応器の挙動をシミュレートするために使用される。例えば、このシミュレーションにより、化学反応の収量、反応時の温度変化などを予測することが可能になる。
【０１２７】
図５ｂでは、運用者は、まず、反応器４６の最適なサイズまたは構成を定義しようとする。このために、運用者は、生産性（反応収量）を改善すること、環境にやさしいこと（廃棄物が少量であること）、または低コストであることなどの達成すべき（複数の）目標条件４７を決定または評価し、従うべき、また例えば、安全性、例えば超えてはならない反応媒体の最高温度を指定すること、または所定の最大または最小量の試薬による生産性または収益性、または有毒な、または有害な、またはリサイクルもしくは生物分解が困難な廃棄物に対する限度が規定されている環境に関係しうる制約条件４８を設定し、さらに反応器の構造内の可変特性４９を設定する。
【０１２８】
これらの可変特性４９は、例えば、反応ブロックの個数、ブロック１つ当たりのユーティリティゾーンの個数、供給点の個数および位置、反応室の寸法、反応室内の反応媒体の流れに関するユーティリティ流体の並流、向流、または交差流の流れ方向、および反応室内の挿入物の個数、性質、および幾何学的形状を含む。
【０１２９】
一変更形態では、達成されるべき目標条件は、複数の定量的基準を含む数学的関数とすることも可能である。例えば、生産ラインに対する目標条件を規定する関数は、生産性基準および環境基準に関して式
ｆ＝ｑｐ×ｋｐ−ｑｘ×ｋｘ
により定義することができ、
ただし、ｑｐは、反応から予期される生成物の量であり、ｋｐは、生成物の単価であり、ｑｘは、価値のない反応生成物の量であり、ｋｘは、この生成物を処理するか、または除去するコストである。
【０１３０】
運用者は、さらに、これらの目標条件および制約条件に対する許容される変量、つまり公認されたマージンを定義することもできる。
【０１３１】
本発明の方法では、上で定義された反応器状態４５、上で指定された動作条件４３からの反応器の寸法または構成、および前記目標条件４７、制約条件４８、および変数４９を最適化する（４６）ことが可能である。最適化は、複数の連続するシミュレーション５０を実行することにより行われ、このため、プレート型反応器の構造の可変特性すべてに従うことにより、制約条件、および場合によっては目標条件および制約条件に対する認められたマージン条件を満たしつつ、可能な限り多数の目標条件を満たすか、または目標関数を得ることができる。
【０１３２】
反応の大きさは、与えられた反応について反応器の寸法を決定したり、または構成したりするのに最適なパラメータを定義する（４６）ことを可能にするＳＱＰ（登録商標）（Ｓｕｃｃｅｓｓｉｖｅ  Ｑｕａｄｒａｔｉｃ  Ｐｒｏｇｒａｍｍｉｎｇ）などの第２のソフトウェアツール５１を使って最適化される。
【０１３３】
図５ｃでは、運用者は、次いで、与えられた反応について反応器の動作を最適化するため最適な動作条件５２を決定しようとする。達成すべき（複数の）目標条件４７および従うべき制約条件４８は、反応器の寸法を最適化している間と変わりない。例えばさまざまな試薬の流量、温度、圧力、および組成、反応室２４内に試薬が供給される順序、ならびにユーティリティ流体の流量および温度を含む試薬変数５３である変数のみが変化する。
【０１３４】
本発明の方法では、反応器のすでに定義されている最適な寸法決定またはサイズ、上述の動作条件、および上述の目標条件、制約条件、および変数に基づいて、与えられた反応に対する反応器の最適な動作条件を決定することにより反応器の動作を最適化することが可能である。最適化は、複数の連続するシミュレーション５０を実行することにより行われ、そのため、反応に対する変数データのすべてに従うことにより、制約条件、および場合によってはさらに目標条件および制約条件に対する認められたマージン条件に従いながら、上述の目標条件または目標関数のできるだけ多くの条件を達成することが可能である。
【０１３５】
反応器の動作は、さらに、上述のソフトウェアツール５１を使って最適化される。最適な動作条件を定めることで、例えば、化学反応の反応器制御、安全性、および感度に関して比較対照研究を行うことが可能である。
【０１３６】
与えられたサイズおよび与えられた反応に対するプレート型反応器の挙動を最適化する実施例については、図６から１２に例示されている。
【０１３７】
図６および７は、チオ硫酸ナトリウムを酸化する反応を追跡するグラフである。横軸１４９は、リットルを単位として反応室の容積を表す。左の縦座標軸１５０は、温度を℃で示しており、右の縦座標軸１５１は、最大１までのスケールで反応の収量を示している。
【０１３８】
曲線１５２は、酸化反応の収量がどのように変化するかを示しており、曲線１５３、１５４、および１５５は、それぞれ、反応媒体、ユーティリティ流体、および反応室を定めるプレートの温度の変化を示している。
【０１３９】
チオ硫酸ナトリウムは、以下の反応を使用して過酸化水素により酸化され、反応生成物はトリチオン酸ナトリウム、硫酸ナトリウム、および水である。
２Ｎａ２Ｓ２Ｏ３＋４Ｈ２Ｏ２→Ｎａ２Ｓ３Ｏ６＋Ｎａ２ＳＯ４＋４Ｈ２Ｏ
【０１４０】
反応は、強い発熱反応であり、３つの反応ブロックを有するプレート型反応器内の均質な液体媒体中で生じる。滞留時間分布の上述の実験研究において、プレート型反応器を一連の９１個の個別セルとして扱うことが可能であり、反応器はそれぞれ２７個の個別セルを含む３つの反応ブロック、および２つの反応ブロックの間に１０個の個別セルを含む遷移要素を備える。この反応に対する反応器の動的モデルにおける方程式の個数は、１９２０である。
【０１４１】
図６では、反応媒体は、低濃度の試薬を有し、そこにチオ硫酸ナトリウムは０．６２モル／リットル（ｍｏｌ．Ｌ−１）の濃度で４０リットル／時（Ｌ．ｈ−１）の速度により注入され、過酸化水素は１．２８ｍｏｌ．Ｌ−１の濃度で１０Ｌ．ｈ−１の速度により注入されており、ユーティリティ流体は、それぞれのブロックにおいて、３．４ｍ３．ｈ−１の流量で供給される温度１４℃の水であった。
【０１４２】
運用者は、単一の目標条件、つまり、最適化反応の収量を改善すること、ユーティリティ流体の流量である反応変数データ、および４５℃以下に保たれなければならない反応媒体に対する温度制約条件を設定する。
【０１４３】
図７では、反応媒体は、より高い濃度の試薬を有し、チオ硫酸ナトリウムは０．７５ｍｏｌ．Ｌ−１の濃度で４０Ｌ．ｈ−１の速度により注入されており、過酸化水素は１．５９ｍｏｌ．Ｌ−１の濃度で１０Ｌ．ｈ−１の速度により注入されており、ユーティリティ流体は第１の反応ブロックでは１．８ｍ３．ｈ−１の流量で、また２つの他の反応ブロックでは０．９ｍ３．ｈ−１の流量で供給される温度１４℃の水である。
【０１４４】
本発明の方法では、ユーティリティ流体流量などの反応変数データに従うことによりチオ硫酸ナトリウム酸化反応についてプレート型反応器を最適化することが可能である。
【０１４５】
次いで、チオ硫酸ナトリウムを酸化する反応が、ユーティリティ流体の最適化された流量で動作するように組み立てられたプレート型反応器内で実行され、温度および圧力などの物理的パラメータの測定が、反応器の動作を検証するためにこの反応の実行時に行われる。
【０１４６】
図６および７では、菱形点１５６は、反応に対する収量の実験測定結果であり、正方形点１５７および三角形点１５８は、それぞれ、反応媒体およびユーティリティ流体の温度の実験測定結果である。
【０１４７】
これらの実験結果は、対応する曲線に近く、そこで、プレート型反応器においてこの化学反応のシミュレーションを検証し、それにより、動的モデルの精度を評価することができる。
【０１４８】
シミュレーションは、さらに、センサの数が小さいことや、その配置の仕方のせいで実験測定結果に表れない図７の温度ピーク１５８（９０℃を超える）を正確に特定するために使用され、最高測定反応は６０℃未満である。
【０１４９】
図８から１０は、無水酢酸を加水分解する反応を追跡するグラフである。横軸１５９は、個別反応セルの個数を表す。この実施形態では、プレート型反応器は、３０個の個別セルを有し、これらのセルは３つの反応ブロックに分けられ、一方の反応室から次のブロックへの通路は縦の破線１６０により線図として表されている。縦座標軸１６１は、温度を℃で表している。曲線１６２、１６３、および１６４は、それぞれ、反応媒体、ユーティリティ流体、および反応室を定める遷移プレートの温度の変化を示している。
【０１５０】
無水酢酸の加水分解により、以下の反応を介して酢酸が生じる。
ＣＨ３−ＣＯ−Ｏ−ＣＯ−ＣＨ３＋Ｈ２Ｏ→２ＣＨ３−ＣＯ−ＯＨ
【０１５１】
運用者は、最大生産性の目標条件、８０℃未満であるべき反応媒体の温度制約条件、および反応器の寸法決定を最適化する反応供給点の個数および配置である反応器構造可変特性を設定した。
【０１５２】
図８では、温度３０℃、流量４０Ｌ．ｈ−１で無水酢酸が反応室のセル１に完全に注入された。水も、温度３０℃、流量１０Ｌ．ｈ−１でセル１内に注入された。ユーティリティ流体は、反応が強い発熱反応であるため温度１５℃、流量１０，８００Ｌ．ｈ−１であった。最適化変数に選択された値では、温度制約条件を満たすことができないことがわかる。
【０１５３】
図９では、温度３０℃、流量４０Ｌ．ｈ−１で無水酢酸が反応室のセル１に完全に注入された。水は、温度３０℃、流量５Ｌ．ｈ−１でセル１（第１のブロックの第１のセル）に、温度３０℃、流量５Ｌ．ｈ−１でセル１１（第２の反応ブロックの第１のセル）に注入された。ユーティリティ流体は、温度１５℃、流量１０，８００Ｌ．ｈ−１であった。
【０１５４】
そこで、水が、反応室に２回、２つの異なる供給点のところで導入された。これにより、反応媒体の最高温度を少し下げることが可能であり、図８の１３８℃から図９の１０４℃になった。反応器内の試薬供給点の個数および配置は、しかしながら、温度制約条件に従っていないため、正しくなかった。
【０１５５】
図１０では、温度３０℃、流量４０Ｌ．ｈ−１で無水酢酸が反応室のセル１に完全に注入された。水は、温度３０℃、流量３．３３Ｌ．ｈ−１でセル１内に、温度３０℃、流量３．３３Ｌ．ｈ−１でセル１１内に、温度３０℃、流量３．３３Ｌ．ｈ−１でセル２１（第１のブロック第１のセル）内に注入された。ユーティリティ流体は、温度１５℃、流量１０，８００Ｌ．ｈ−１であった。
【０１５６】
この場合、水は、反応室に３回、３つの異なる供給点のところで導入された。これにより、反応媒体の最高温度を８０℃の制約条件よりも低くすることが可能である。そこで最適化方法により、生産性目標条件を満たし、温度制約条件に従う適切な寸法決定を識別することが可能であった。
【０１５７】
図１１および１２は、ネオペンチルグリコールを生成する反応を追跡するグラフである。横軸１６５は、リットルを単位として反応室の容積を表す。この実施では、反応室は、それぞれ約０．０３５リットル（Ｌ）の容積を有する１０個の個別セルからなる４つのブロックにさらに分けられる４０個の個別セルを備えていた。左側の縦座標軸１６６は、温度を℃で表し、右側の縦座標軸１６７は、反応の収量をパーセンテージで表している。
【０１５８】
曲線１６８は、反応の収量がどのように変化するかを示しているが、曲線１６９、１７０、および１７１は、それぞれ、反応媒体、ユーティリティ流体、および反応室を定めるプレートの温度の変化を示している。
【０１５９】
ネオペンチルグリコールを生成する化学式は以下の通りである。
２（ホルムアルデヒド）＋（２−メチルプロパノール）＋（水酸化ナトリウム）→（ネオペンチルグリコール）＋（ギ酸ナトリウム）
【０１６０】
図１１および１２では、縦の破線１７２および１７３は、それぞれ、水酸化ナトリウム供給点および２−メチルプロパノール供給点を表す。そこで、水酸化ナトリウムは、反応室の容積が、反応室内の第２のセル（個別セル１つ当たり０．０３５Ｌ）に対応する、０．０５Ｌに近くなったときに１回だけ導入され、２−メチルプロパノールは、第３のセル、第１５のセル、および第２８のセルと３回導入された。
【０１６１】
図１１では、反応が発熱反応であれば、ユーティリティ流体は温度５０℃、流量１．５ｍ３．ｈ−１であった。全流量９．１４Ｌ．ｈ−１に対して、試薬は、温度１８℃で、２−メチルプロパノールは、濃度２ｍｏｌ．ｋｇ−１で注入された（３回の注入は、同一流量３．０５Ｌ．ｈ−１であった）。達成された転化率は、８９．４％であり、反応中の反応媒体の最高温度は、６６．６℃であった。
【０１６２】
運用者は、ネオペンチルグリコールを生成する反応からの収量を高める目標条件、それぞれの供給点について２−メチルプロパノールの供給速度である反応の活動を最適化する変数データ、および６５℃を超えるべきでない反応媒体の最高温度である従うべき制約条件を設定した。
【０１６３】
結果は、図１２に示されており、そこでは、ユーティリティ流体は、温度５０℃、流量１．５ｍ３．ｈ−１となっており、試薬は、温度１８℃で供給され、２−メチルプロパノールは、全流量９．１４Ｌ．ｈ−１に対し濃度２ｍｏｌ．ｋｇ−１で供給される。２−メチルプロパノールの３回の注入は、異なる流量で行われる、つまり、流量はセル３への第１の供給では５．９１Ｌ．ｈ−１、セル１１への第２の供給では２．７３Ｌ．ｈ−１、セル２１への第３の供給では０．５０Ｌ．ｈ−１である。達成された転化率は、９５．６％であり、反応中の反応媒体の最高温度は、６５℃に到達し、したがって課されている制約条件を満たしている。
【０１６４】
そのため、この最適化方法では、この場合複数の地点で試薬を反応室に供給する速度である反応変数データに従って温度制約条件を満たしつつ、反応に対する目標収量を達成することが可能である。
【図面の簡単な説明】
【０１６５】
【図１】本発明の方法で使用される開放プレート型反応器の分解斜視線図である。
【図２】複数の個別セルにさらに分けられた図１の反応器の反応室の線図である。
【図３】図２の反応室の個別セルの線図である。
【図４】２つの個別ユーティリティ流体循環セルに関連付けられている個別反応セルの線図である。
【図５ａ】本発明の最適化方法の主要工程を示す図である。
【図５ｂ】本発明の最適化方法の主要工程を示す図である。
【図５ｃ】本発明の最適化方法の主要工程を示す図である。
【図６】試薬濃度および流量、ならびにユーティリティ流体流量が、チオ硫酸ナトリウム酸化反応の温度および収量に及ぼす影響を示すグラフである。
【図７】試薬濃度および流量、ならびにユーティリティ流体流量が、チオ硫酸ナトリウム酸化反応の温度および収量に及ぼす影響を示すグラフである。
【図８】試薬供給点の個数および配置が、無水酢酸加水分解反応の温度に及ぼす影響を示すグラフである。
【図９】試薬供給点の個数および配置が、無水酢酸加水分解反応の温度に及ぼす影響を示すグラフである。
【図１０】試薬供給点の個数および配置が、無水酢酸加水分解反応の温度に及ぼす影響を示すグラフである。
【図１１】複数の供給点における試薬供給速度が、ネオペンチルグリコールを生成する反応の温度および収量に及ぼす影響を示すグラフである。
【図１２】複数の供給点における試薬供給速度が、ネオペンチルグリコールを生成する反応の温度および収量に及ぼす影響を示すグラフである。
【符号の説明】
【０１６６】
Ｂ、Ｅ、Ｅ’、Ｈ’  個別ユーティリティ流体フローセル
１、２、３、４、５、６、７、８、９、１０、１１、１２  セル
２０、２１、２２、２３  プレート
２４  反応室
２５ａ  試薬供給点
２５ｂ  供給点
２６  出口点
２７  左側室
２８  右側室
２９  入口
３０  出口
３７  積分ソフトウェアツール
３８  反応器
３９  物理的特性
４０  モデル
４１  熱伝達
４２  液体−液体物質移動または液体−気体物質移動
４３  動作条件
４４  基本構成の記述
４５  「反応器状態」
４６  制約条件
４７  目標条件
４９  可変特性
５０  シミュレーション
５１  第２のソフトウェアツール
５２  動作条件
１４９  横軸
１５０  左の縦座標軸
１５１  右の縦座標軸
１５２、１５３、１５４、および１５５  曲線
１５６  菱形点
１５７  正方形点
１５８  三角形点
１５９  横軸
１６０  破線
１６２、１６３、および１６４  曲線
１６５  横軸
１６６  左側の縦座標軸
１６７  右側の縦座標軸
１６８、１６９、１７０、および１７１  曲線
【図１】
2008529757.tif 000027
【図２】
2008529757.tif 000028
【図３】
2008529757.tif 000029
【図４】
2008529757.tif 000030
【図５ａ】
2008529757.tif 000031
【図５ｂ】
2008529757.tif 000032
【図５ｃ】
2008529757.tif 000033
【図６】
2008529757.tif 000034
【図７】
2008529757.tif 000035
【図８】
2008529757.tif 000036
【図９】
2008529757.tif 000037
【図１０】
2008529757.tif 000038
【図１１】
2008529757.tif 000039
【図１２】
2008529757.tif 000040