(57)【要約】【課題】製品品質を保ちつつバッチサイクル時間を短縮する。【解決手段】反応速度定数推定部２ａは、バッチ運転実績データに基づいて開始剤の反応速度定数と原料モノマーの反応速度定数とその温度依存係数を推定する。モデル調整部７は、バッチ運転実績データを用いて熱収支モデル１１ａのモデルパラメータを調整する。反応温度目標値軌道決定部８は、所望の終了時点における原料モノマー濃度の推定値が所望の濃度と一致するように反応温度目標値を決定する。冷媒入口温度時間パターン設定部３は、重合プロセスモデルと反応温度目標値とを用いて冷媒入口温度目標値を決定する。制御パラメータ調整部５は、重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラ（フィードフォワード制御実行部４およびカスケード制御実行部６）の制御パラメータを調整する。【選択図】  図９
2018156152.tif 000002
【特許請求の範囲】
【請求項１】
バッチ重合プロセスの過去のバッチ運転により得られた、反応器の内部の反応温度と反応の徐熱に用いる冷媒の入口温度と前記冷媒の出口温度と前記冷媒の循環流量との組からなる時系列データに基づいて、前記反応器に投入される開始剤および原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定する第１のステップと、
前記時系列データを用いて前記反応温度と前記冷媒の温度とに関する熱収支モデルのパラメータを調整する第２のステップと、
前記反応速度定数と前記温度依存係数とによって規定される重合反応モデルから計算される、所望の終了時点における前記原料の濃度の推定値が所望の濃度と一致するように、反応温度目標値の時間パターンを決定する第３のステップと、
前記重合反応モデルと前記熱収支モデルとからなる重合プロセスモデルを用いたシミュレーションにより、前記反応温度が前記反応温度目標値の時間パターンと一致するように前記冷媒の入口温度目標値の時間パターンを決定する第４のステップと、
前記重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラの制御パラメータを調整する第５のステップと、
バッチ重合プロセスにおいて前記冷媒の入口温度が前記入口温度目標値と一致し且つ前記反応温度が前記反応温度目標値と一致するように前記コントローラによって制御を行なう第６のステップとを含むことを特徴とする制御方法。
【請求項２】
請求項１記載の制御方法において、
前記第１のステップは、
前記時系列データから反応発熱量を推定する第７のステップと、
前記反応発熱量の推定値の時間積分値に基づいて、バッチ経過時間におけるバッチ重合プロセスの生成物質の濃度を推定する第８のステップと、
前記生成物質の濃度の推定値に基づいて、未反応で残留する前記原料の濃度を推定する第９のステップと、
前記時系列データと前記未反応の原料の濃度の推定値とに基づいて、前記開始剤および前記原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定する第１０のステップとを含むことを特徴とする制御方法。
【請求項３】
請求項１または２記載の制御方法において、
前記第３のステップは、前記反応温度の昇温速度の上限と前記反応温度目標値の上下限とを制約条件として、前記原料の濃度の推定値と前記所望の濃度との２乗誤差を最小とする最適化問題を解くことにより、前記反応温度目標値の時間パターンを決定するステップを含むことを特徴とする制御方法。
【請求項４】
バッチ重合プロセスの過去のバッチ運転により得られた、反応器の内部の反応温度と反応の徐熱に用いる冷媒の入口温度と前記冷媒の出口温度と前記冷媒の循環流量との組からなる時系列データに基づいて、前記反応器に投入される開始剤および原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定するように構成された反応速度定数推定部と、
前記時系列データを用いて前記反応温度と前記冷媒の温度とに関する熱収支モデルのパラメータを調整するように構成されたモデル調整部と、
前記反応速度定数と前記温度依存係数とによって規定される重合反応モデルから計算される、所望の終了時点における前記原料の濃度の推定値が所望の濃度と一致するよう、反応温度目標値の時間パターンを決定するように構成された反応温度目標値軌道決定部と、
前記重合反応モデルと前記熱収支モデルとからなる重合プロセスモデルを用いたシミュレーションにより、前記反応温度が前記反応温度目標値の時間パターンと一致するように前記冷媒の入口温度目標値の時間パターンを決定するように構成された冷媒入口温度時間パターン設定部と、
前記重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラの制御パラメータを調整するように構成された制御パラメータ調整部と、
バッチ重合プロセスにおいて前記冷媒の入口温度が前記入口温度目標値と一致し且つ前記反応温度が前記反応温度目標値と一致するように制御を行なうコントローラとを備えることを特徴とする制御装置。
【請求項５】
請求項４記載の制御装置において、
前記反応速度定数推定部は、
前記時系列データから反応発熱量を推定する反応発熱量推定部と、
前記反応発熱量の推定値の時間積分値に基づいて、バッチ経過時間におけるバッチ重合プロセスの生成物質の濃度を推定する生成物質濃度推定部と、
前記生成物質の濃度の推定値に基づいて、未反応で残留する前記原料の濃度を推定する原料濃度推定部と、
前記時系列データと前記未反応の原料の濃度の推定値とに基づいて、前記開始剤および前記原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定する反応速度定数・温度依存係数推定部とを備えることを特徴とする制御装置。
【請求項６】
請求項４または５記載の制御装置において、
前記反応温度目標値軌道決定部は、前記反応温度の昇温速度の上限と前記反応温度目標値の上下限とを制約条件として、前記原料の濃度の推定値と前記所望の濃度との２乗誤差を最小とする最適化問題を解くことにより、前記反応温度目標値の時間パターンを決定することを特徴とする制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、バッチ重合プロセスをモデル化した反応プロセスモデルを適応的に調整し、調整結果をバッチ重合プロセスのＰＩＤ制御に反映するＢ２Ｂ制御方式の制御方法および制御装置に関するものである。
【背景技術】
【０００２】
容積が数リットルの反応装置を使った試作実験などにより製造処方が確立された高機能性ポリマーを、スケールアップした数ｍ3の反応器で変種変量生産するバッチ重合プロセスを制御対象とする場合、新製品のため品質規格が厳しく、初回のバッチ運転から直ちに、品質を規定する反応温度を精密に制御することが求められる。このため、発明者は、プロセスモデルをベースにしたＢ２Ｂ（Batch to Batch）制御システムを提案した（特許文献１参照）。
しかしながら、従来のＢ２Ｂ制御では、加温開始からポリマーの抜き出しを完了するまでのバッチサイクル時間が長いという課題があった。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−６９０９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明は、上記課題を解決するためになされたもので、昇温過程のあるバッチ重合プロセスにおいて、製品品質を保ちつつバッチサイクル時間を短縮して生産性を向上することができる制御方法および制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の制御方法は、バッチ重合プロセスの過去のバッチ運転により得られた、反応器の内部の反応温度と反応の徐熱に用いる冷媒の入口温度と前記冷媒の出口温度と前記冷媒の循環流量との組からなる時系列データに基づいて、前記反応器に投入される開始剤および原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定する第１のステップと、前記時系列データを用いて前記反応温度と前記冷媒の温度とに関する熱収支モデルのパラメータを調整する第２のステップと、前記反応速度定数と前記温度依存係数とによって規定される重合反応モデルから計算される、所望の終了時点における前記原料の濃度の推定値が所望の濃度と一致するように、反応温度目標値の時間パターンを決定する第３のステップと、前記重合反応モデルと前記熱収支モデルとからなる重合プロセスモデルを用いたシミュレーションにより、前記反応温度が前記反応温度目標値の時間パターンと一致するように前記冷媒の入口温度目標値の時間パターンを決定する第４のステップと、前記重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラの制御パラメータを調整する第５のステップと、バッチ重合プロセスにおいて前記冷媒の入口温度が前記入口温度目標値と一致し且つ前記反応温度が前記反応温度目標値と一致するように前記コントローラによって制御を行なう第６のステップとを含むことを特徴とするものである。
【０００６】
また、本発明の制御方法の１構成例において、前記第１のステップは、前記時系列データから反応発熱量を推定する第７のステップと、前記反応発熱量の推定値の時間積分値に基づいて、バッチ経過時間におけるバッチ重合プロセスの生成物質の濃度を推定する第８のステップと、前記生成物質の濃度の推定値に基づいて、未反応で残留する前記原料の濃度を推定する第９のステップと、前記時系列データと前記未反応の原料の濃度の推定値とに基づいて、前記開始剤および前記原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定する第１０のステップとを含むことを特徴とするものである。
また、本発明の制御方法の１構成例において、前記第３のステップは、前記反応温度の昇温速度の上限と前記反応温度目標値の上下限とを制約条件として、前記原料の濃度の推定値と前記所望の濃度との２乗誤差を最小とする最適化問題を解くことにより、前記反応温度目標値の時間パターンを決定するステップを含むことを特徴とするものである。
【０００７】
また、本発明の制御装置は、バッチ重合プロセスの過去のバッチ運転により得られた、反応器の内部の反応温度と反応の徐熱に用いる冷媒の入口温度と前記冷媒の出口温度と前記冷媒の循環流量との組からなる時系列データに基づいて、前記反応器に投入される開始剤および原料のそれぞれの反応速度定数と前記原料の反応速度定数の温度依存係数とを推定するように構成された反応速度定数推定部と、前記時系列データを用いて前記反応温度と前記冷媒の温度とに関する熱収支モデルのパラメータを調整するように構成されたモデル調整部と、前記反応速度定数と前記温度依存係数とによって規定される重合反応モデルから計算される、所望の終了時点における前記原料の濃度の推定値が所望の濃度と一致するよう、反応温度目標値の時間パターンを決定するように構成された反応温度目標値軌道決定部と、前記重合反応モデルと前記熱収支モデルとからなる重合プロセスモデルを用いたシミュレーションにより、前記反応温度が前記反応温度目標値の時間パターンと一致するように前記冷媒の入口温度目標値の時間パターンを決定するように構成された冷媒入口温度時間パターン設定部と、前記重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラの制御パラメータを調整するように構成された制御パラメータ調整部と、バッチ重合プロセスにおいて前記冷媒の入口温度が前記入口温度目標値と一致し且つ前記反応温度が前記反応温度目標値と一致するように制御を行なうコントローラとを備えることを特徴とするものである。
【発明の効果】
【０００８】
本発明では、バッチ重合プロセスの過去のバッチ運転により得られた、反応温度と冷媒の入口温度と冷媒の出口温度と冷媒の循環流量との組からなる時系列データに基づいて、開始剤および原料のそれぞれの反応速度定数と原料の反応速度定数の温度依存係数とを推定し、前記時系列データを用いて反応温度と冷媒の温度とに関する熱収支モデルのパラメータを調整し、前記反応速度定数と前記温度依存係数とによって規定される重合反応モデルから計算される、所望の終了時点における原料の濃度の推定値が所望の濃度と一致するように、反応温度目標値の時間パターンを決定し、前記重合反応モデルと前記熱収支モデルとからなる重合プロセスモデルを用いたシミュレーションにより、冷媒の入口温度目標値の時間パターンを決定し、前記重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラの制御パラメータを調整し、バッチ重合プロセスにおいて冷媒の入口温度が入口温度目標値と一致し且つ反応温度が反応温度目標値と一致するようにコントローラによって制御を行なう。本発明によれば、昇温過程のあるバッチ重合プロセスにおいて、生成物質の品質を保ちつつバッチサイクル時間を短縮して生産性を向上することができる。
【図面の簡単な説明】
【０００９】
【図１】従来のバッチ重合プロセスの計装図である。
【図２】従来の制御装置の構成を示すブロック図である。
【図３】従来の制御装置の動作を説明するフローチャートである。
【図４】従来の制御系のブロック線図である。
【図５】従来の冷媒入口温度目標値の時間パターンの１例を示す図である。
【図６】従来の反応温度目標値の時間パターンの１例を示す図である。
【図７】従来のＩＩ2−ＰＤコントローラでプロセスをフィードバック制御する制御系のブロック線図である。
【図８】従来のモデル調整部の構成を示すブロック図である。
【図９】本発明の実施例に係る制御装置の構成を示すブロック図である。
【図１０】本発明の実施例に係る制御装置の動作を説明するフローチャートである。
【図１１】本発明の実施例に係る制御装置の反応速度定数推定部の動作を説明するフローチャートである。
【図１２】本発明の実施例におけるバッチ運転実績データの１例を示す図である。
【図１３】本発明の実施例のシミュレーション時に計算した伝熱係数を示す図である。
【図１４】本発明の実施例において反応発熱量とポリマー濃度と未反応の原料モノマーの濃度とを推定した結果を示す図である。
【図１５】本発明の実施例における原料モノマー濃度の推定値と真値とを示す図である。
【図１６】本発明の実施例における反応温度目標値の時間パターンの１例を示す図である。
【図１７】本発明の実施例における原料モノマー濃度と反応温度目標値の時間パターンとポリマー濃度とバッチサイクルの短縮時間の１例を示す図である。
【発明を実施するための形態】
【００１０】
図１は特許文献１に開示されたバッチ重合プロセスの計装図である。バッチ重合プロセスにおいては、完全混合槽型の反応器１００に、バルブ１０５〜１０７を介して原料モノマーと溶媒と助剤とを定量仕込み、反応温度まで加温する。反応器１００は、冷却のためのコイル１０１と、同じく冷却のためのジャケット１０２と、反応物質の反応温度を測定する温度センサ１０３と、反応物質を攪拌する攪拌機１０４とを備えている。
【００１１】
原料モノマーの加温後に、バルブ１０８を介して反応器１００に重合開始剤（触媒）を定量加えると、数分以内の誘導時間（induction period）だけ遅れて急激にラジカル重合反応が始まる。反応温度を一定に保つ主反応は約２時間で、その後さらに反応温度を昇温して反応を加速し、開始剤の添加から８時間程度で反応終点となる。生成されたポリマーは、反応器１００からバルブ１０９を介して抜き出される。加温開始から生成ポリマーの抜き出しを完了するまでのバッチ時間は概ね１２時間である。
【００１２】
重合反応による発熱量は、反応器１００に設けられたコイル１０１とジャケット１０２とを一定流量で循環する冷媒（冷水）により除熱される。冷媒は、入口側配管１１０からコイル１０１とジャケット１０２に供給され、出口側配管１１１に排出される。出口側配管１１１に排出された冷媒は、一部が冷却塔１１３によって冷却され、残りが入口側配管１１０に戻されるようになっている。入口側配管１１０には、冷媒入口温度を測定する温度センサ１１２が設けられている。また、入口側配管１１０には、バルブ１１４，１１５を介して蒸気と冷却塔１１３によって冷却された冷水とが供給される。
【００１３】
制御装置１１６は、冷媒入口温度が冷媒入口温度目標値の時間パターンと一致し、かつ反応温度が反応温度目標値の時間パターンと一致するように、冷媒入口温度制御と反応温度制御のカスケード制御を行なう。また、制御装置１１６は、バッチ経過時間に対応して冷媒入口温度の目標値を変更する、フィードフォワード制御（Feed-Forward Control）の機能を有している。
【００１４】
図２は特許文献１に開示された従来の制御装置１１６の構成を示すブロック図、図３は制御装置１１６の動作を説明するフローチャートである。制御装置１１６は、重合プロセスのダイナミクスモデルである重合プロセスモデルをあらかじめ記憶するモデル記憶部１と、バッチ重合プロセスの実験データから重合の反応速度定数を推定する反応速度定数推定部２と、反応温度を反応温度目標値の時間パターン｛ｔ，Ｔrset｝に一致させる冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝をシミュレーションにより定める冷媒入口温度時間パターン設定部３と、フィードフォワード制御を実行するフィードフォワード制御実行部４と、カスケード制御のための制御パラメータを調整する制御パラメータ調整部５と、カスケード制御を実行するカスケード制御実行部６と、重合プロセスモデルを適応的に調整するモデル調整部７とを備えている。
【００１５】
バッチ運転を開始する前に、モデル調整部７は、バッチ運転実績データを用いて重合プロセスモデルのモデルパラメータを調整し（図３ステップＳ１）、制御パラメータ調整部５は、重合プロセスモデルを線形近似した伝達関数モデルを用いて制御パラメータ（ＰＩＤパラメータ）を調整する（ステップＳ２）。上記のようなバッチ重合プロセスを制御対象とするバッチ運転が開始されると、フィードフォワード制御実行部４は、冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝に応じて冷媒入口温度目標値Ｔcisetを出力するフィードフォワード制御を実行し（ステップＳ３）、カスケード制御実行部６は、冷媒入口温度Ｔciが冷媒入口温度目標値と一致し且つ反応温度Ｔrが反応温度目標値と一致するようにカスケード制御（ＰＩＤ制御）を実行する（ステップＳ４）。以上のようなステップＳ３，Ｓ４の処理を、バッチ重合プロセスが終了するまで（ステップＳ５においてＹＥＳ）、一定間隔のバッチ経過時間ｔ毎に実行する。
【００１６】
重合プロセスモデルは、ラジカル重合の重合反応モデル１０と、熱収支モデル１１とからなる。重合反応モデル１０は、開始剤の重量濃度（残存量の仕込量に対する分率）ｘiの時間変化および原料モノマーの重量濃度ｘmの時間変化を表すものである。ここで、開始剤の反応速度定数ｋiおよび原料モノマーの反応速度定数ｋmは共に０より大である。そして、反応温度Ｔrと冷媒温度Ｔcとに関する熱収支モデル１１では、反応物質の熱容量をＣr、冷媒の熱容量をＣc、重合反応熱（−ΔＨr）＞０による発熱量をＱr、反応物質と冷媒との熱交換による除熱量をＱc1、循環冷媒による除熱量をＱc2とする。このとき、バッチ経過時間ｔにおける重合プロセスモデルは式（１）〜式（４）のようになる。
【００１７】
【数１】
2018156152.tif 000003
【００１８】
【数２】
2018156152.tif 000004
【００１９】
【数３】
2018156152.tif 000005
【００２０】
【数４】
2018156152.tif 000006
【００２１】
式（１）、式（２）が重合反応モデル１０を表し、式（３）、式（４）が熱収支モデル１１を表している。開始剤の反応速度定数は一定とし、原料モノマーの反応速度定数だけに温度依存性があるとする。この原料モノマーの反応速度定数をアレニウス式（Arrhenius rate expression）で表す。
【００２２】
【数５】
2018156152.tif 000007
【００２３】
ここで、Ｔr0は反応温度Ｔrの基準状態値、ｋm0は反応速度定数ｋmの基準状態値、ｂmは温度依存定数である。重合レートＲmは、原料モノマー濃度の変化速度の絶対値にモノマー仕込量Ｗmを乗じた式（６）で与えられる。
【００２４】
【数６】
2018156152.tif 000008
【００２５】
熱収支モデル１１の発熱量と除熱量は次のようになる。コイル１０１とジャケット１０２を循環する冷媒の入口温度Ｔciと冷媒の出口温度Ｔcoとの差が小さいので、冷媒温度ＴcはＴciとＴcoの算術平均（Ｔci＋Ｔco）／２とした。重合反応熱（−ΔＨr）＞０による発熱量Ｑr（ｔ）、反応物質と冷媒との熱交換による除熱量Ｑc1（ｔ）、循環冷媒による除熱量Ｑc2（ｔ）は、式（７）、式（８）、式（９）のようになる。
【００２６】
【数７】
2018156152.tif 000009
【００２７】
【数８】
2018156152.tif 000010
【００２８】
【数９】
2018156152.tif 000011
【００２９】
ここで、ｃcは冷媒の比熱、ｆcは冷媒の循環流量である。また、Ａtは全伝熱面積であり、コイル１０１の伝熱面積をＡc、ジャケット１０２の伝熱面積をＡjとすると、Ａt＝Ａc＋Ａjとなる。また、コイル１０１の総括伝熱係数をＵc、ジャケット１０２の総括伝熱係数をＵjとすると、総括伝熱係数Ｕc，Ｕjをコイル１０１とジャケット１０２のそれぞれの伝熱面積Ａc，Ａjで加重平均した総括伝熱係数Ｕ（ｔ）は、Ｕ（ｔ）＝（Ｕc（ｔ）Ａc＋Ｕj（ｔ）Ａj）／Ａtとなる。総括伝熱係数Ｕc，Ｕjは、モノマー濃度などから推算することができる。
【００３０】
反応速度定数推定部２は、原料モノマーの重量濃度（残存量の仕込量に対する分率）ｘmの実験データから、開始剤の反応速度定数ｋiおよび原料モノマーの反応速度定数ｋmを推定する。この反応速度定数ｋi，ｋmの推定方法は、文献「小河  守正，“バッチ重合プロセスのモデルベースＢ２Ｂ制御”，計測自動制御学会論文集，Ｖｏｌ．４６，Ｎｏ．３，ｐ．１３９−１４８，２０１０」に開示されているので、詳細な説明は省略する。反応速度定数推定部２の処理は、事前に１回だけ行っておけばよい。
【００３１】
冷媒入口温度時間パターン設定部３は、シミュレーション実行部３０と、時系列テーブル記憶部３１とを有する。シミュレーション実行部３０は、重合プロセスモデルを用いたシミュレーションにより、反応温度Ｔrが所望の反応温度目標値の時間パターン｛ｔ，Ｔrset｝と一致するように冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝を定める。時系列テーブル記憶部３１は、シミュレーション実行部３０が定めた冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝を記憶する。冷媒入口温度時間パターン設定部３の処理は、事前に１回だけ行っておけばよい。
【００３２】
制御パラメータ調整部５は、プロセス動特性モデル作成部５０と、ゲインスケジューリング部５１とを有する。特許文献１では、冷媒入口温度制御と反応温度制御のカスケード制御のアルゴリズムとしてＰＩＤを用いるが、ＰＩＤによる制御を設計するには、操作量（冷媒入口温度）と制御量（反応温度）との関係を表わす、線形のプロセス動特性モデルが必要になる。
【００３３】
プロセス動特性モデル作成部５０は、このような線形のプロセス動特性モデルを作成するものである。具体的には、プロセス動特性モデル作成部５０は、モデル記憶部１に記憶されている重合プロセスモデルを逐次線形化して状態方程式モデルを作成し、さらに状態方程式モデルを伝達関数モデルに変換する。ゲインスケジューリング部５１は、プロセス動特性モデル作成部５０が作成した伝達関数モデルを用いて、制御パラメータ（ＰＩＤパラメータ）を調整する。
【００３４】
以下、プロセス動特性モデル作成部５０の動作について説明する。バッチプロセスの状態変数は常に変化し、連続プロセスでは一般的な定常状態（すべての状態変数が一定の状態）が存在しない。そこで、バッチ経過時間ｔに対応して与えられる反応温度目標値軌道｛ｔ，Ｔrset｝に、反応温度Ｔrを完全に追従制御できると仮定する。そのときの重合レートは式（１）〜式（４）の重合反応モデルと式（５）、式（６）から定まり、さらに総括伝熱係数も決まる。
【００３５】
そして、反応温度目標値軌道｛ｔ，Ｔrset｝に沿った重合レートのもとで、式（１）〜式（４）における熱収支モデルの左辺＝０となる状態、すなわち反応温度Ｔrが目標値軌道｛ｔ，Ｔrset｝に常に一致するように冷媒入口温度Ｔciが調節されている状態を、仮想的な平衡状態と呼ぶことにする。この平衡状態のまわりで熱収支モデルだけを逐次線形化することにより、反応温度のプロセス動特性モデルを得る。
【００３６】
例えば、バッチ経過時間ｔの平衡状態における重合レートを＊Ｒmのように表わす。このとき、反応温度の微小変化ΔＴr（ｔ）による重合レート変化量ΔＲm（ｔ）は式（１０）で表わされる。
【００３７】
【数１０】
2018156152.tif 000012
【００３８】
ここで、＊Ｔrはバッチ経過時間ｔの平衡状態における反応温度、Δｋmは原料モノマーの反応速度定数ｋmの変化量である。バッチ経過時間ｔの平衡状態における重合反応熱（−ΔＨr）＞０による発熱量を＊Ｑrとすると、発熱量の変化量ΔＱr（ｔ）は、＊Ｑr＝＊Ｒm（−ΔＨr）より、式（１１）のようになる。
【００３９】
【数１１】
2018156152.tif 000013
【００４０】
逐次線形化したうえで、平衡状態からの変化量をΔＴr（ｔ）→Ｔr（ｔ）のように書き換えると、状態方程式モデルは式（１２）、式（１３）のように表わされる。
【００４１】
【数１２】
2018156152.tif 000014
【００４２】
【数１３】
2018156152.tif 000015
【００４３】
状態量ｘ（ｔ）、操作量ｕ（ｔ）、制御量ｙ（ｔ）は、それぞれ式（１４）、式（１５）、式（１６）のようになる。
【００４４】
【数１４】
2018156152.tif 000016
【００４５】
【数１５】
2018156152.tif 000017
【００４６】
【数１６】
2018156152.tif 000018
【００４７】
さらに、式（１２）、式（１３）の係数Ａ，ｂ，ｃは式（１７）のようになる。ここで、＊Ｕはバッチ経過時間ｔの平衡状態における総括伝熱係数である。
【００４８】
【数１７】
2018156152.tif 000019
【００４９】
次に、伝達関数モデルの作成方法について説明する。状態方程式モデルから、伝達関数モデルＰ（ｓ）は式（１８）のようになる。なお、ｓはラプラス演算子である。
【００５０】
【数１８】
2018156152.tif 000020
【００５１】
式（１８）の係数ｎ0，ｄ2，ｄ1は式（１９）のようになる。
【００５２】
【数１９】
2018156152.tif 000021
【００５３】
式（１８）の分母多項式の根の逆数である時定数Ｔp1，Ｔp2は式（２０）のようになる。
【００５４】
【数２０】
2018156152.tif 000022
【００５５】
ゲインＫpは式（２１）で表される。
【００５６】
【数２１】
2018156152.tif 000023
【００５７】
時定数Ｔp1，Ｔp2のうち動特性を支配する時定数をＴp1（｜Ｔp1｜＞Ｔp2＞０）とする。時定数Ｔp1が負のとき不安定プロセスとなり、不安定プロセスとなる条件は式（２２）のように表される。
【００５８】
【数２２】
2018156152.tif 000024
【００５９】
不安定プロセスとなるときは式（２１）で表わされるゲインＫpも負になる。重合反応の進行に伴い、不安定から安定なシステムに移行していくプロセスを、伝達関数モデルで見通しよく表現するために、ゲインと時定数は絶対値を取り常に正の値にする。すなわち、｜Ｋp｜→Ｋp、｜Ｔp1｜→Ｔp1とする。そして、プロセスむだ時間ＴLを付加した伝達関数モデルは式（２３）のようになる。
【００６０】
【数２３】
2018156152.tif 000025
【００６１】
制御システムを構成する温度センサ、操作部さらにカスケード制御２次ループの遅れなどが重なったものを、プロセスむだ時間ＴLとして付加することができる。式（２３）の分母第１項の±符号が、正のとき安定プロセス、負になると不安定プロセスである。ゲインＫpおよび時定数Ｔp1の計算値が負になる場合に、このゲインＫpおよび時定数Ｔp1を絶対値に置き換え、±符号を負にする。
不安定条件を物理的に解釈するために状態方程式モデルに戻ると、ａ11＜０のとき不安定プロセスになる。すなわち、式（２４）が成立する。
【００６２】
【数２４】
2018156152.tif 000026
【００６３】
バッチ経過時間ｔの平衡状態における反応物質と冷媒との熱交換による除熱量を＊Ｑc1、平衡状態における冷媒温度をＴcとすると、平衡状態＊Ｑr＝＊Ｑc1＝＊ＵＡt（＊Ｔr−＊Ｔc）から＊ＵＡt＝＊Ｑr／（＊Ｔr−＊Ｔc）なので、不安定条件は次の簡単な関係式で与えられる。
【００６４】
【数２５】
2018156152.tif 000027
【００６５】
このように、反応速度定数に温度催存性があり、反応温度Ｔrと冷媒温度Ｔcの温度差が式（２５）の条件を満たすとき、重合プロセスは不安定になる。この不安定さは、重合レートが高く発熱量が大きい反応初期において、ポリマーの生成により総括伝熱係数Ｕが急激に低下するため、冷媒温度Ｔcを下げ反応温度Ｔrとの温度差を大きくとり除熱量を確保している状態を意味する。
以上で、プロセス動特性モデル作成部５０による伝達関数モデルの作成処理が終了する。
【００６６】
次に、フィードフォワード制御実行部４とカスケード制御実行部６の動作について説明する。図４は特許文献１に開示された従来の制御系のブロック線図である。図４におけるＦはフィードフォワード制御実行部４が実現するフィードフォワードコントローラである。フィードフォワード制御実行部４は、冷媒入口温度時間パターン設定部３の時系列テーブル記憶部３１に予め記憶されている冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝を参照し、現在のバッチ経過時間ｔに対応する冷媒入口温度目標値Ｔcisetを出力する。図５に冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝の１例を示す。このように冷媒入口温度目標値の時間パターンは、バッチ経過時間ｔと冷媒入口温度目標値Ｔcisetとの組からなる時系列データである。
【００６７】
また、図４におけるＣ1はカスケード制御実行部６が実現する反応温度コントローラ（Ｉ−ＰＤコントローラまたはＩＩ2−ＰＤコントローラ）、Ｃ2は同じくカスケード制御実行部６が実現する冷媒入口温度コントローラ（ＰＩ−Ｄコントローラ）、Ｐ1は重合プロセス、Ｐ2は冷却プロセス、ｄ1，ｄ2は外乱である。
【００６８】
反応温度コントローラＣ1には、反応温度目標値Ｔrsetが与えられる。図６に反応温度目標値の時間パターン｛ｔ，Ｔrset｝の１例を示す。このように反応温度目標値の時間パターンは、バッチ経過時間ｔと反応温度目標値Ｔrsetとの組からなる時系列データである。反応温度コントローラＣ1は、反応温度Ｔrと反応温度目標値Ｔrsetとが一致するように操作量ｕ1を演算する。操作量ｕ1と冷媒入口温度目標値Ｔcisetとは加算され、この加算結果が冷媒入口温度目標値ｒ2として冷媒入口温度コントローラＣ2に与えられる。冷媒入口温度コントローラＣ2は、冷媒入口温度Ｔciと冷媒入口温度目標値ｒ2とが一致するように操作量ｕ2を演算する。この操作量ｕ2に応じて、図１に示したバルブ１１４，１１５の開度が決定される。
【００６９】
次に、カスケード制御実行部６が実現する反応温度コントローラＣ1について説明する。まず、反応温度コントローラＣ1としてＩ−ＰＤコントローラを用いる場合について説明する。カスケード制御２次ループの冷媒入口温度制御は、混合プロセスが制御対象なので、非常に速い応答特性を持つ。さらに、式（２３）の伝達関数モデルにおける、冷媒の熱容量による遅れに対応する時定数とむだ時間はそれぞれ０．２ｍｉｎ程度なので、時定数Ｔp2およびプロセスむだ時間ＴLは（Ｔp1／１０）に対して十分に小さい。そこで、時定数Ｔp2およびプロセスむだ時間ＴLを無視して、プロセス動特牲を式（２６）に示すように１次遅れ特性と見なす。
【００７０】
【数２６】
2018156152.tif 000028
【００７１】
Ｉ−ＰＤコントローラは、ＰＩＤパラメータを｛Ｋc，Ｔi，Ｔd｝、微分ゲインを１／γ、目標値をｒ（ｓ）、制御量をｙ（ｓ）、制御量偏差をｅ（ｓ）＝ｒ（ｓ）−ｙ（ｓ）、操作量をｕ（ｓ）とすると、式（２７）で表される。
【００７２】
【数２７】
2018156152.tif 000029
【００７３】
Ｋcは比例ゲイン、Ｔiは積分時間、Ｔdは微分時間である。また、目標値ｒ（ｓ）は図４の反応温度目標値Ｔrsetに相当し、制御量ｙ（ｓ）は図４の反応温度Ｔrに相当し、操作量ｕ（ｓ）は図４の操作量ｕ1に相当する。Ｉ−ＰＤコントローラの目標値応答特性Ｗc（ｓ）＝ｙ（ｓ）／ｒ（ｓ）は式（２８）のようになる。
【００７４】
【数２８】
2018156152.tif 000030
【００７５】
望ましい目標値応答特性Ｗd（ｓ）を２次臨界制動とし、応答の速さを規定するチューニングパラメータＴFを用いる。
【００７６】
【数２９】
2018156152.tif 000031
【００７７】
目標値応答特性Ｗc（ｓ），Ｗd（ｓ）の分母多項式の係数相等条件から、微分時間Ｔdがゼロとなる条件を付加して、比例ゲインＫcと積分時間Ｔiが求まる。そして、微分時間Ｔdは、経験的に積分時間Ｔiの１／８に設定する。このようにして得られるＰＩＤ設定則を表１に示す。ここで、ｑ≡ＴF／Ｔp1＜１である。
【００７８】
【表１】
2018156152.tif 000032
【００７９】
表１に示すＰＩＤパラメータ｛Ｋc，Ｔi，Ｔd｝および微分ゲイン１／γがカスケード制御実行部６に予め設定されており、カスケード制御実行部６は反応温度コントローラＣ1（Ｉ−ＰＤコントローラ）の動作を実現することができる。
【００８０】
次に、反応温度コントローラＣ1としてＩＩ2−ＰＤコントローラを用いる場合について説明する。Ｉ−ＰＤコントローラは良い制御性能を示すが、反応温度目標値Ｔrsetと反応温度Ｔrとの間に微小な定常偏差が生じることがある。この定常偏差は、バッチ経過時間を通して、ランプ状外乱が作用していることによる。プロセスの伝達関数が原点に極をもたない場合、フィードバック制御システムがランプ状外乱にオフセットフリーとなるには、内部モデル原理から、コントローラは自身の極にランプ状外乱の原点の２位極を含まなければならない。すなわち、２重積分動作が必要になる。
【００８１】
ＩＩ2−ＰＤコントローラでプロセスＰ（ｓ）をフィードバック制御する制御系のブロック線図を図７に示す。図７において、Ｋcは比例ゲイン、Ｔi1は積分時間、Ｔi2は２重積分時間、Ｔdは微分時間、１／γは微分ゲイン、ｄは外乱である。図７に示すＳ（ｓ）はカスケード制御実行部６が実現する目標値補償器であり、式（３０）のように表される。
【００８２】
【数３０】
2018156152.tif 000033
【００８３】
２重積分動作による目標値応答の行き過ぎを抑えるため、目標値ｒ（ｓ）を目標値補償器Ｓ（ｓ）に通し、目標値補償器Ｓ（ｓ）の出力ｒs（ｓ）をコントローラに与える。外乱ｄ（ｓ）は操作量ｕ（ｓ）に加法的に作用するプロセス入力等価外乱とする。ＩＩ2−ＰＤコントローラは、積分動作と２重積分動作が制御偏差ｅ（ｓ）＝ｒs（ｓ）−ｙ（ｓ）に作用し、比例動作と微分動作（微分ゲイン１／γ）は制御量ｙ（ｓ）だけに働く。これにより、１組のＰＩＩ2Ｄパラメータ｛Ｋc，Ｔi1，Ｔi2，Ｔd｝で、目標値追従性と外乱抑制性を同時に満たす２自由度制御システムを構成する。
【００８４】
プロセス伝達関数モデルとして式（２３）を用い、一般性のあるＰＩＩ2Ｄ設定則を与える。ランプ状外乱は、持続時間ｔrの間に操作量等価な外乱量ｄrだけ、定速度ａ≡ｄr／ｔrで変化する。
【００８５】
【数３１】
2018156152.tif 000034
【００８６】
この外乱に対する制御量ｙ（ｓ）の望ましい応答特性Ｗd（ｓ）を一般化臨界制御応答とする。制御システムの因果律を保つため、プロセス動特性モデルの相対次数ｊに対して一般化臨界制動の次数ｎ＝ｊ＋２とし、制御量応答をプロセスむだ時間だけ遅らせる。ここには、制御性能を調整する２つのパラメータがある。整定時間を決めるパラメータＴFと、制御量の減衰比を規定するパラメータαＫpである。
【００８７】
【数３２】
2018156152.tif 000035
【００８８】
ランプ状外乱に対する望ましい制御量応答を実現するＰＩＩ2Ｄ設定則を導く。ＩＩ2−ＰＤ制御に内在する目標値補償器Ｓ（ｓ）を図７のように目標値入力側に設けると、コントローラを次式のＰＩＩ2ＤコントローラＣ1（ｓ）としてＰＩＩ2Ｄ設定則を設計できる。
【００８９】
【数３３】
2018156152.tif 000036
【００９０】
外乱応答が望ましい制御量応答に一致する条件から、このＰＩＩ2Ｄコントローラは式（３４）のように定まる。
【００９１】
【数３４】
2018156152.tif 000037
【００９２】
式（３４）に、次数ｎ＝４とした式（３２）のＷd（ｓ）と式（２３）のＰ（ｓ）とを代入し、むだ時間をマクローリン（Maclaurin）展開近似すると、分数多項式表現は次のようになる。αは無次元で、ランプ状外乱速度ａが時間の逆数の次元ｈ−１を持つので、分数多項式の係数の次元はいずれも一致する。
【００９３】
【数３５】
2018156152.tif 000038
【００９４】
式（３３）のＰＩＩ2Ｄアルゴリズムの不完全微分項を無視したものと式（３５）を分数多項式で表し、係数相等条件から表２のＰＩＩ2Ｄ設定則を得る。
【００９５】
【表２】
2018156152.tif 000039
【００９６】
表２に示すＰＩＤパラメータ｛Ｋc，Ｔi1，Ｔi2，Ｔd｝および微分ゲイン１／γがカスケード制御実行部６に予め設定されており、カスケード制御実行部６は反応温度コントローラＣ1（ＩＩ2−ＰＤコントローラ）の動作を実現することができる。
冷媒入口温度コントローラＣ2については一般的なＰＩ−Ｄコントローラなので、説明は省略する。
【００９７】
次に、モデル調整部７の動作について説明する。Ｂ２Ｂ制御は、バッチ運転の実績データを用いて重合プロセスモデルを適応的に調整し、その結果をモデルベース制御に反映するものである。同品種の実績データを使うので、重合反応モデル１０を固定し、熱収支モデル１１を適応させる。熱収支モデル１１の主な変動要因は、反応物質粘度とコイル・ジャケット総括伝熱係数Ｕの推算誤差、および伝熱面の汚れによる総括伝熱係数Ｕの変化である。
【００９８】
図８はモデル調整部７の構成を示すブロック図である。モデル調整部７は、実績データ記憶部７０と、平滑化部７１と、２乗誤差算出部７２と、非線形最適化部７３と、モデルパラメータ調整部７４とを有する。特許文献１で用いるバッチ運転実績データは、反応器１００で生成しようとするポリマーと同種のポリマーを過去に生産したときの実績データ（前回の実績データ）であり、バッチ経過時間ｔ、冷媒入口温度Ｔci、反応温度Ｔrの時系列データ｛ｔ（ｋ），Ｔci（ｋ），Ｔr（ｋ）｝である。ここで、ｋはサンプリング周期τ（τは例えば１０ｓｅｃ）の離散時間ｔkの時点である。実績データ記憶部７０は、このバッチ運転実績データを予め記憶している。つまり、実績データ記憶部７０は、前回のバッチ運転における開始から終了までのｎサンプルの時系列データを記憶している。
【００９９】
入力手段となる平滑化部７１は、冷媒入口温度Ｔci（ｋ）を平滑化した冷媒入口温度Ｔciav（ｋ）を、重合反応モデル１０と熱収支モデル１１とからなる重合プロセスモデルの入力として与える。これにより、重合プロセスモデルからは、出力である反応温度Ｔrav（ｋ）が得られる。平滑化部７１は、ｎサンプルの冷媒入口温度Ｔciav（ｋ）を重合プロセスモデルに順次入力するので、重合プロセスモデルからは、ｎサンプルの反応温度Ｔrav（ｋ）が順次得られることになる。冷媒入口温度Ｔci（ｋ）を平滑化する理由は、冷媒入口温度Ｔci（ｋ）の変動が大きいからである。平滑化の手法としては、例えばローパスフィルタがある。
【０１００】
２乗誤差算出部７２は、ｎサンプルの反応温度Ｔrav（ｋ）とｎサンプルの反応温度の実績値Ｔr（ｋ）との２乗誤差Ｅを式（３６）のように計算する。
【０１０１】
【数３６】
2018156152.tif 000040
【０１０２】
非線形最適化部７３は、２乗誤差Ｅを最小とする非線形最適化問題を数値的に解き、重合プロセスモデルの伝熱特性適応パラメータλ∈｛λ1，λ2，λ3｝を得る。伝熱特性適応パラメータλ∈｛λ1，λ2，λ3｝は、次の意味を持つ。まず、総括伝熱係数Ｕは反応物質粘度μや攪拌機回転数Ｎなどにより変化する。伝熱特性適応パラメータλ1は、この総括伝熱係数Ｕを調整するためのものである。すなわち、総括伝熱係数Ｕは、伝熱特性適応パラメータλ1により式（３７）のように調整される。
【０１０３】
【数３７】
2018156152.tif 000041
【０１０４】
ｇ（）は関数である。具体的には、コイル１０１の総括伝熱係数Ｕc、ジャケット１０２の総括伝熱係数Ｕjは、それぞれ式（３８）、式（３９）のように調整される。
【０１０５】
【数３８】
2018156152.tif 000042
【０１０６】
【数３９】
2018156152.tif 000043
【０１０７】
反応物質粘度μは原料モノマー濃度ｘmにより変わる。伝熱特性適応パラメータλ2，λ3は、この反応物質粘度μを調整するためのものである。すなわち、反応物質粘度μは、伝熱特性適応パラメータλ2，λ3により式（４０）のように調整される。
【０１０８】
【数４０】
2018156152.tif 000044
【０１０９】
ｆ（）は関数である。具体的には、初期粘度μ0とモノマー濃度依存定数μ1と上記の次数ｎとパラメータλ2，λ3とを用いて、反応物質粘度μは式（４１）のように調整される。
【０１１０】
【数４１】
2018156152.tif 000045
【０１１１】
コイル１０１の総括伝熱係数Ｕcと総括伝熱係数Ｕcの境膜伝熱係数ｈcとは、式（４２）のような関係にある。
【０１１２】
【数４２】
2018156152.tif 000046
【０１１３】
ジャケット１０２の総括伝熱係数Ｕjと総括伝熱係数Ｕjの境膜伝熱係数ｈjとは、式（４３）のような関係にある。
【０１１４】
【数４３】
2018156152.tif 000047
【０１１５】
そして、境膜伝熱係数ｈcは、反応物質粘度μと攪拌機回転数Ｎを用いて式（４４）のように表される。
【０１１６】
【数４４】
2018156152.tif 000048
【０１１７】
μwは境膜面におけるモノマー粘度である。境膜伝熱係数ｈjは、式（４５）のように表される。
【０１１８】
【数４５】
2018156152.tif 000049
【０１１９】
なお、式（４２）〜式（４５）は実験的に求めた１例の式であって、これに限るものではない。
モデルパラメータ調整部７４は、非線形最適化部７３が求めた伝熱特性適応パラメータλ1を用いて、モデルパラメータである総括伝熱係数Ｕを式（３７）、式（３８）、式（３９）のように調整すると共に、非線形最適化部７３が求めた伝熱特性適応パラメータλ2，λ3を用いて、モデルパラメータである反応物質粘度μを式（４０）、式（４１）のように調整する。こうして、重合プロセスモデルが適応的に調整され、モデル調整部７の動作が終了する。
【０１２０】
次に、制御パラメータ調整部５のゲインスケジューリング部５１の動作について説明する。前述のとおり、制御パラメータ調整部５のプロセス動特性モデル作成部５０は、モデル記憶部１に記憶されている重合プロセスモデルを線形近似して伝達関数モデルを作成する。
【０１２１】
ゲインスケジューリング部５１は、プロセス動特性モデル作成部５０が作成した伝達関数モデルを用いて、制御パラメータ（ＰＩＤパラメータ）を調整する。反応温度コントローラＣ1としてＩ−ＰＤコントローラを用いる場合、ゲインスケジューリング部５１は、ＰＩＤパラメータ｛Ｋc，Ｔi，Ｔd｝を調整する。また、反応温度コントローラＣ1としてＩＩ2−ＰＤコントローラを用いる場合、ゲインスケジューリング部５１は、ＰＩＤパラメータ｛Ｋc，Ｔi1，Ｔi2，Ｔd｝を調整する。ゲインスケジューリングについては、例えば文献「D.E.Seborg et al.，“Process Dynamics and Control”，Second Edition，John Wiley & Sons，Inc.，pp.428-430，2003」に開示されているので、詳細な説明は省略する。以上の構成は特許文献１に開示されている。
【０１２２】
［実施例］
次に、本実施例の特徴的な構成について説明する。本実施例では、これまで述べてきた反応温度精密制御を活用して生産性を向上させる。すなわち、昇温反応過程のあるバッチ重合プロセスにおいて、製品品質を保ちつつバッチサイクル時間の短縮を図る。このために、制御装置に２つの機能を増強する。
【０１２３】
本実施例においてもバッチ重合プロセスの計装図は従来と同様であるので、図１の符号を用いて説明する。図９は本実施例の制御装置１１６ａの構成を示すブロック図であり、図２と同一の構成には同一の符号を付してある。
【０１２４】
本実施例の制御装置１１６ａは、重合プロセスのダイナミクスモデルである重合プロセスモデルをあらかじめ記憶するモデル記憶部１ａと、バッチ重合プロセスの実験データから重合の反応速度定数を推定する反応速度定数推定部２ａと、反応温度を反応温度目標値の時間パターン｛ｔ，Ｔrset｝に一致させる冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝をシミュレーションにより定める冷媒入口温度時間パターン設定部３と、フィードフォワード制御を実行するフィードフォワード制御実行部４と、カスケード制御のための制御パラメータを調整する制御パラメータ調整部５と、カスケード制御を実行するカスケード制御実行部６と、重合プロセスモデルを適応的に調整するモデル調整部７と、反応温度目標値の時間パターン｛ｔ，Ｔrset｝を決定する反応温度目標値軌道決定部８とを備えている。
【０１２５】
図１０は本実施例の制御装置１１６ａの動作を説明するフローチャートである。バッチ運転を開始する前に、反応速度定数推定部２ａは、バッチ運転実績データに基づいて開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定する（図１０ステップＳ１０）。
モデル調整部７は、バッチ運転実績データを用いて重合プロセスモデル（熱収支モデル１１ａ）のモデルパラメータを調整する（図１０ステップＳ１１）。
【０１２６】
反応温度目標値軌道決定部８は、所望の終了時点における原料モノマー濃度の推定値が所望のモノマー濃度と一致するように、反応温度目標値の時間パターン｛ｔ，Ｔrset｝を決定する（図１０ステップＳ１２）。
冷媒入口温度時間パターン設定部３は、重合プロセスモデルと反応温度目標値軌道決定部８が決定した反応温度目標値の時間パターン｛ｔ，Ｔrset｝とを用いて冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝を決定して記憶する（図１０ステップＳ１３）。
【０１２７】
制御パラメータ調整部５は、重合プロセスモデルを線形近似した伝達関数モデルを用いてコントローラ（フィードフォワード制御実行部４およびカスケード制御実行部６）の制御パラメータ（ＰＩＤパラメータ）を調整する（図１０ステップＳ１４）。
【０１２８】
バッチ重合プロセスを制御対象とするバッチ運転が開始されると、フィードフォワード制御実行部４は、冷媒入口温度目標値の時間パターン｛ｔ，Ｔciset｝に応じて冷媒入口温度目標値Ｔcisetを出力するフィードフォワード制御を実行し（図１０ステップＳ１５）、カスケード制御実行部６は、冷媒入口温度Ｔciが冷媒入口温度目標値と一致し且つ反応温度Ｔrが反応温度目標値と一致するようにカスケード制御（ＰＩＤ制御）を実行する（図１０ステップＳ１６）。以上のようなステップＳ１５，Ｓ１６の処理を、バッチ重合プロセスが終了するまで（図１０ステップＳ１７においてＹＥＳ）、一定間隔のバッチ経過時間ｔ毎に実行する。
【０１２９】
本実施例の増強機能の第一は、バッチ運転実績データに基づく反応速度定数の推定である。反応機構は明確でも実験データが不備なために、反応速度定数とその温度依存係数のパラメータを推定できないことがある。この場合に，バッチ運転実績データから反応発熱量を推算して、その数値積分値からポリマー（生成物質）の濃度を求める。このポリマー濃度から未反応で残留する原料モノマー濃度が決まるので、以下の方法でパラメータ推定できる。
【０１３０】
［反応速度定数推定の原理］
以下、重合反応モデルの反応速度定数をパラメータ推定する原理を数値例と共に示す。ポリマーの製造処方を開発する実験により、反応時間（バッチ経過時間）ｔ（ｎ）のときの原料モノマー濃度ｘm（ｎ）の実測値が与えられる。そのｎmax組のデータ｛ｔ（ｎ），ｘm（ｎ）｜ｎ＝１，２，・・・・，ｎmax｝から開始剤の反応速度定数ｋiおよび原料モノマーの反応速度定数ｋmを推定する。
【０１３１】
重合反応モデルは式（１）〜式（４）に示した。その反応速度定数（パラメータ）推定問題を定式化するために、重合反応モデルの状態変数とパラメータの表記法を次の式（４６）に示すように変える。重合反応モデルの状態変数の初期値は、ｘ1（０）＝１、ｘ2（０）＝１である。
【０１３２】
【数４６】
2018156152.tif 000050
【０１３３】
重合反応モデルは次の式（４７）のようになる。
【０１３４】
【数４７】
2018156152.tif 000051
【０１３５】
ここで、ｘ，θ，ｆ（ｘ，θ），ｇ（ｘ，θ），ｈは次のとおりである。
【０１３６】
【数４８】
2018156152.tif 000052
【０１３７】
パラメータ推定問題を一般化して定式化するために、状態変数と出力変数などを以下の式（４９）のように拡張する。ここでは、複数ｋmaxの出力変数の実験データであるチルダｙが与えられるものとしている。なお、以降では数式の頭に付した「〜」をチルダと呼ぶ。
【０１３８】
【数４９】
2018156152.tif 000053
【０１３９】
重合反応モデル解ｙ＝ｈとの残差ｅ＝チルダｙ−ｈを対角行列Ｗで重み付けした２次形式の目的関数Φを式（５０）のように定義する。
【０１４０】
【数５０】
2018156152.tif 000054
【０１４１】
この目的関数のｎ組目を式（５１）のように展開する。
【０１４２】
【数５１】
2018156152.tif 000055
【０１４３】
パラメータ推定問題は、後述するように非線形最適化問題として定式化できる。それをソルバーで数値的に解く際に、目的関数の勾配（gradient）∇Φを与え、効率的に解を得るようにする。ｎ組目の目的関数の勾配は、重み行列が対角行列なので、次の式（５２）のようになる。
【０１４４】
【数５２】
2018156152.tif 000056
【０１４５】
∂ｈ（ｎ）／∂θは直接求められないので、状態変数をパラメータで１次微分したモデル感度行列（model sensitibity matrix）Ｓ（ｎ）≡（∂ｘ（ｎ）／∂θ）Tを導入して次の式（５３）のように表す。
【０１４６】
【数５３】
2018156152.tif 000057
【０１４７】
このとき目的関数の勾配∇Φは次式のようになる。
【０１４８】
【数５４】
2018156152.tif 000058
【０１４９】
目的関数の勾配∇Φはｎmax組の実験データのそれぞれの勾配の総和なので次の式（５５）のようになる。
【０１５０】
【数５５】
2018156152.tif 000059
【０１５１】
次に、モデル感度行列は次の式（５６）のように定義され、ベクトル微分∂ｘ（ｎ）／∂θの転置であることに注意する。
【０１５２】
【数５６】
2018156152.tif 000060
【０１５３】
このモデル感度行列は解析的に求まらないので、式（５７）のように、時間に関する微分方程式で表現する。
【０１５４】
【数５７】
2018156152.tif 000061
【０１５５】
ｆi（ｘ，θ）に関してｘ，θは共に独立変数である。しかし、重合反応モデルの制約を受け、θが変化すると、重合反応モデルでは従属変数であるｘが変わる。このため、∂ｆi／∂θjはθjそのものの変化と、それによるｘの変化の両方を考慮しなければならない。そこで、次式のように連鎖律（chain rule）を適用する。
【０１５６】
【数５８】
2018156152.tif 000062
【０１５７】
式（５８）の右辺第１項がθjの変化によるｘの変化分によるもの、右辺第２項がθjそのものの変化分によるものである。ここで、偏微分の右添字はそれらの変数が一定であることを表している。たとえばθ≠jは、θj以外のθ1，θ2，・・・・，θj-1，θj+1，・・・・，θjmaxが一定であることを表している。
【０１５８】
そして、ｆx＝∂ｆ／∂ｘ，ｆθ＝∂ｆ／∂θとおくと、次式のように書くことができる。
【０１５９】
【数５９】
2018156152.tif 000063
【０１６０】
さらに次の式（６０）、式（６１）のように行列を定義する。
【０１６１】
【数６０】
2018156152.tif 000064
【０１６２】
【数６１】
2018156152.tif 000065
【０１６３】
これにより、式（６２）のようなモデル感度行列の微分方程式を得る。
【０１６４】
【数６２】
2018156152.tif 000066
【０１６５】
このモデル感度行列の微分方程式を重合反応モデルと共に制約に加えると、パラメータ推定問題は次の式（６３）のように、非線形最適化問題として定式化される。さらに、非線形最適化問題のソルバーには、式（５５）で表される目的関数の勾配∇Φを与える。
【０１６６】
【数６３】
2018156152.tif 000067
【０１６７】
本実施例では、実際の重合反応モデルの反応速度定数を、実験データに基づきパラメータ推定する。重合反応モデルは、すでに式（４７）に示したとおりである。２次形式の目的関数Φを次の式（６４）のように定義する。出力変数が一つなので、重み行列はスカラーでＷ＝１である。
【０１６８】
【数６４】
2018156152.tif 000068
【０１６９】
目的関数Φは、実験データであるチルダｙ（ｎ）＝チルダｘ2（ｎ）と推定パラメータに基づく重合反応モデル解ｙ（ｎ）＝ｈ（ｎ）＝ｘ2（ｎ）の残差ｅ（ｎ）＝チルダｙ（ｎ）−ｙ（ｎ）の２乗和である。ここで、残差ベクトルｅ＝［ｅ1  ｅ2  ｅnmax］Tである。
推定パラメータの単位変化に対する目的関数の変化を表す勾配は、式（５５）から次の式（６５）のようになる。
【０１７０】
【数６５】
2018156152.tif 000069
【０１７１】
ここで、モデル感度行列Ｓの微分方程式は次の式（６６）のようになる。
【０１７２】
【数６６】
2018156152.tif 000070
【０１７３】
ｆx，ｆθは式（６７）のようになる。
【０１７４】
【数６７】
2018156152.tif 000071
【０１７５】
また、∂ｈ／∂ｘは次のとおりである。
【０１７６】
【数６８】
2018156152.tif 000072
【０１７７】
モデル感度行列Ｓの微分方程式はＭＡＴＬＡＢ optimization toolboxの非線形最適化ソルバーfminconなどにより数値的に解ける。数値計算を効率化するため目的関数の勾配を与える。以上が、反応速度定数推定の原理である。次に、より実際的な推定処理について説明する。
【０１７８】
［バッチ運転実績データに基づく反応速度定数の推定］
バッチ運転実績データから、反応発熱量とポリマー濃度を推算する。その推算結果から原料モノマー濃度も定まるので、バッチ時間とモノマー濃度推定値の時系列データから、開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定する。
【０１７９】
［重合プロセスモデル］
従来と同様に、重合プロセスモデルは、ラジカル重合の重合反応モデル１０と、熱収支モデル１１ａとからなり、バッチ反応開始からの経過時間ｔを用いて式（６９）〜式（７３）で表せる。
【０１８０】
【数６９】
2018156152.tif 000073
【０１８１】
重合反応モデル１０は、開始剤の重量濃度（残存量の仕込量に対する分率）ｘiの時間変化および原料モノマーの重量濃度ｘmの時間変化を表すものである。開始剤の反応速度定数ｋiおよび原料モノマーの反応速度定数ｋmは共に０より大である。一方、本実施例の熱収支モデル１１ａは、式（３）、式（４）に、反応器１００のコイル１０１とジャケット１０２の伝熱壁の熱容量に関する式を追加したものである。
【０１８２】
ここで、Ｔrは反応温度、ＴMは伝熱壁温度、Ｔcは冷媒温度である。また、反応物質（反応液）の熱容量をＣr、伝熱壁の熱容量をＣM、冷媒の熱容量をＣc、重合反応熱（−ΔＨr）＞０による発熱量をＱr、反応物質と伝熱壁との熱交換による除熱量をＱc1、伝熱壁と循環冷媒との熱交換による除熱量をＱc2、循環冷媒による除熱量をＱc3としている。式（６９）、式（７０）が重合反応モデル１０を表し、式（７１）〜式（７３）が熱収支モデル１１ａを表している。
【０１８３】
開始剤の反応速度定数ｋiは一定とし、原料モノマーの反応速度定数ｋmだけに温度依存性があるとする。この原料モノマーの反応速度定数ｋmをアレニウス式で表すと、式（５）のようになる。Ｔr0は反応温度Ｔrの基準状態値、ｋm0は反応速度定数ｋmの基準状態値、ｂmは温度依存定数である。重合レートＲmは、原料モノマー濃度の変化速度の絶対値にモノマー仕込量Ｗmを乗じた式（６）で与えられる。
【０１８４】
熱収支モデル１１ａの発熱量と除熱量は次のようになる。コイル１０１とジャケット１０２を循環する冷媒の入口温度Ｔciと冷媒の出口温度Ｔcoとの差が数℃と小さいので、冷媒温度ＴcはＴciとＴcoの算術平均で近似する。重合反応熱（−ΔＨr）＞０による発熱量Ｑr（ｔ）、反応物質と伝熱壁との熱交換による除熱量Ｑc1（ｔ）、伝熱壁と循環冷媒との熱交換による除熱量Ｑc2（ｔ）、循環冷媒による除熱量Ｑc3（ｔ）は、それぞれ式（７４）、式（７５）、式（７６）、式（７７）のようになる。
【０１８５】
【数７０】
2018156152.tif 000074
【０１８６】
ここで、ｃcは冷媒の比熱、ｆcは冷媒の循環流量である。コイル１０１の反応物質側の伝熱面積Ａbrとジャケット１０２の反応物質側の伝熱面積Ａjrとを合わせた反応物質側の伝熱面積Ａr、コイル１０１の冷媒側の伝熱面積Ａbcとジャケット１０２の冷媒側の伝熱面積Ａjcとを合わせた冷媒側の伝熱面積Ａc、コイル１０１の反応物質側の伝熱係数Ｕbr（ｔ）とジャケット１０２の反応物質側の伝熱係数Ｕjr（ｔ）とを伝熱面積Ａbr，Ａjrで加重平均した反応物質側の伝熱係数Ｕr（ｔ）、コイル１０１の冷媒側の伝熱係数Ｕbcとジャケット１０２の冷媒側の伝熱係数Ｕjcとを伝熱面積Ａbc，Ａjcで加重平均した冷媒側の伝熱係数Ｕc（ｔ）は次の式（７８）のとおりである。原料モノマー濃度の影響を受けて、時変な伝熱係数Ｕbr（ｔ），Ｕjr（ｔ）は運転条件から推算することが可能である。
【０１８７】
【数７１】
2018156152.tif 000075
【０１８８】
［ポリマー濃度とモノマー濃度の推算］
本実施例で用いるバッチ運転実績データは、反応器１００で生成しようとするポリマーと同種のポリマーを過去に生産したときの実績データ（前回の実績データ）であり、バッチ経過時間ｔ、反応温度Ｔr、冷媒入口温度Ｔci、冷媒出口温度Ｔco、冷媒循環流量ｆcの時系列データ｛ｔ（ｋ），Ｔr（ｋ），Ｔci（ｋ），Ｔco（ｋ），ｆc（ｋ）｝である。ここで、ｋはサンプリング周期Δｔの離散時間ｔkの時点である。このバッチ運転実績データは、モデル調整部７の実績データ記憶部７０に記憶され、反応速度定数推定部２ａに与えられるようになっている。
【０１８９】
反応発熱量Ｑrは式（７９）のように２次ルンゲクッタ（Runge-kutta）近似して再帰的に求まる。
【０１９０】
【数７２】
2018156152.tif 000076
【０１９１】
バッチ反応の開始時間０から終了時間ｔbまでの反応発熱量Ｑrの時間積分が生成ポリマー全量に対応することから、バッチ経過時間ｔにおけるポリマー濃度であるハットｘpを次式から推定できる。なお、数式の頭に付した「∧」をハットと呼ぶ。ハットが付与された値は推定値であることを示している。
【０１９２】
【数７３】
2018156152.tif 000077
【０１９３】
そして、未反応の原料モノマーの濃度であるハットｘmは次の式（８１）のようになる。
【０１９４】
【数７４】
2018156152.tif 000078
【０１９５】
［反応速度定数と温度依存係数の推定法］
バッチ経過時間ｔと反応温度Ｔrと未反応の原料モノマー濃度の推定値ハットｘmの時系列｛ｔ（ｋ），Ｔr（ｋ），ハットｘm（ｋ）｝を用いて、開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定する。この推定は、任意のモデルパラメータによるモノマー濃度計算値であるハットｘm（ｋ）の誤差ｅ（ｋ）≡チルダｘm（ｋ）−ハットｘm（ｋ）の２次形式を最小化する最適化問題として、次の式（８２）のように定式化できる。
【０１９６】
【数７５】
2018156152.tif 000079
【０１９７】
ここで、状態変数ｘ（ｔ）≡［ｘi（ｔ）  ｘm（ｔ）］T、最適化変数θ≡［ｋi  ｋm0  ｂm］Tである。
【０１９８】
【数７６】
2018156152.tif 000080
【０１９９】
以上のようにして、反応速度定数推定部２ａは、開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定することができる。図１１は反応速度定数推定部２ａの処理の流れを説明するフローチャートである。図９に示すように、反応速度定数推定部２ａは、反応発熱量推定部２０と、ポリマー濃度推定部２１（生成物質濃度推定部）と、モノマー濃度推定部２２（原料濃度推定部）と、反応速度定数・温度依存係数推定部２３とから構成される。
【０２００】
反応速度定数推定部２ａの反応発熱量推定部２０は、バッチ運転実績データ｛ｔ（ｋ），Ｔr（ｋ），Ｔci（ｋ），Ｔco（ｋ），ｆc（ｋ）｝から反応発熱量Ｑr（ｋ）を推定する（図１１ステップＳ１００）。反応速度定数推定部２ａのポリマー濃度推定部２１は、反応発熱量Ｑr（ｋ）の推定値ハットＱr（ｋ）の時間積分値からバッチ経過時間ｔにおけるポリマー濃度ｘp（ｔ）を推定する（図１１ステップＳ１０１）。反応速度定数推定部２ａのモノマー濃度推定部２２は、ポリマー濃度ｘp（ｔ）の推定値ハットｘp（ｔ）から未反応で残留する原料モノマーの濃度ｘm（ｔ）を推定する（図１１ステップＳ１０２）。
【０２０１】
そして、反応速度定数推定部２ａの反応速度定数・温度依存係数推定部２３は、バッチ運転実績データ｛ｔ（ｋ），Ｔr（ｋ），Ｔci（ｋ），Ｔco（ｋ），ｆc（ｋ）｝に含まれる反応温度Ｔrの時系列Ｔr（ｋ）と未反応の原料モノマーの濃度の推定値ハットｘm（ｋ）と、原料モノマーの濃度ｘmの実験データとから、開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定する（図１１ステップＳ１０３）。以上で、反応速度定数推定部２ａの動作が終了する。従来と同様に、反応速度定数推定部２ａの処理は、事前に１回だけ行っておけばよい。
【０２０２】
［数値例］
バッチ運転実績データ｛ｔ（ｋ），Ｔr（ｋ），Ｔci（ｋ），Ｔco（ｋ），ｆc（ｋ）｝の１例を図１２（Ａ）、図１２（Ｂ）に示す。図１２（Ａ）、図１２（Ｂ）の横軸はバッチ経過時間ｔ［ｈ］、図１２（Ａ）の縦軸は温度［℃］、図１２（Ｂ）の縦軸は流量［ｋｇ／ｈ］である。このバッチ運転実績データは、以下に示す高機能性ポリマーの運転条件に基づき、式（６９）〜式（７３）の重合プロセスモデルを用いて、制御シミュレーションにより作成したものである。
【０２０３】
本実施例の数値例で用いる高機能性ポリマーの運転条件を表３に示す。
【０２０４】
【表３】
2018156152.tif 000081
【０２０５】
そして，重合反応モデルのパラメータを真値ｋi＝０．０５９２、ｋm0＝０．３８２、ｂm＝６９００（Ｔr0＝７０℃）として、表３の運転条件で制御シミュレーションを行なった。その結果得られた、反応温度、冷媒入口温度、冷媒出口温度、冷媒循環流量の時系列｛ｔ（ｋ），Ｔr（ｋ），Ｔci（ｋ），Ｔco（ｋ），ｆc（ｋ）｝をバッチ運転実績データとした。
【０２０６】
また、シミュレーション時に計算した伝熱係数を図１３（Ａ）、図１３（Ｂ）に示す。図１３（Ａ）はコイル１０１の反応物質側の伝熱係数Ｕbrとコイル１０１の冷媒側の伝熱係数Ｕbcを示し、図１３（Ｂ）はジャケット１０２の反応物質側の伝熱係数Ｕjrとジャケット１０２の冷媒側の伝熱係数Ｕjcを示している。図１３（Ａ）、図１３（Ｂ）の横軸はバッチ経過時間ｔ［ｈ］、縦軸は伝熱係数［ｋａｃｌ／ｍ2ｈＫ］である。
【０２０７】
図１２（Ａ）、図１２（Ｂ）の例では、バッチ経過時間ｔ＝２時間までは反応温度Ｔr＝７０℃の等温過程である。そして１時間でＴr＝８０℃まで昇温反応する。その後１時間はＴr＝８０℃を保持し、さらに１時間かけてＴr＝９０℃まで昇温反応し、等温での熟成反応に移る。このバッチ運転実績データを用いて、反応発熱量Ｑrとポリマー濃度ｘpと未反応の原料モノマーの濃度ｘmとを推定した結果を図１４（Ａ）、図１４（Ｂ）に示す。図１４（Ａ）、図１４（Ｂ）の横軸はバッチ経過時間ｔ［ｈ］、図１４（Ａ）の縦軸は熱量［ｋｃａｌ／ｈ］、図１４（Ｂ）の縦軸は濃度［ｗｔ．％］である。
【０２０８】
そして、反応温度Ｔrの時系列と未反応の原料モノマーの濃度の推定値ハットｘmの時系列とを用いることにより、開始剤の反応速度定数ｋi＝０．０６６１７、原料モノマーの反応速度定数ｋmの基準状態値ｋm0＝０．３８４７、反応速度定数ｋmの温度依存係数ｂm＝７３４１という推定結果を得た。なお、反応温度Ｔrの基準状態値Ｔr0＝７０℃である。
【０２０９】
これらの真値は、上記のようにｋi＝０．０５９２、ｋm0＝０．３８２、ｂm＝６９００である。推定値と真値によるモノマー濃度を比較し図１５（Ａ）、図１５（Ｂ）に示した。図１５（Ａ）は未反応の原料モノマーの濃度の真値ｘmと推定値ハットｘmとを示し、図１５（Ｂ）は真値ｘmと推定値ハットｘmとの誤差を示している。
【０２１０】
真値ｘmと推定値ハットｘmとの最大誤差は−０．４７ｗｔ．％で、許容誤差±２ｗｔ．％以内である。このように，重合反応の実験データが不備であっても、バッチ運転実績データから重合反応モデルのパラメータを推定することができる。重合反応モデルが不明確な対象の反応温度精密制御を設計する上で、反応速度定数推定部２ａの機能は有用である。
【０２１１】
［反応温度目標値軌道の決定］
次に、本実施例の増強機能の第二は、反応速度定数の温度依存性を利用して、バッチサイクル時間を短縮する反応温度目標値軌道の決定である。その目標値軌道で運転した場合の製品仕上がり品質を予測し、品質規格を満たすことを確認する。品質不良が予測される場合は、望ましいバッチサイクル時間を変更して目標値軌道の決定操作を繰り返す。そして、得られた目標値軌道に従って反応温度を精密制御する。
【０２１２】
多段階の昇温反応過程を持つバッチ重合プロセスにおいて、バッチ処方で規定された平均的な反応温度目標値の軌道｛ｔ，Ｔrset｝を変更することによってバッチサイクル時間を短縮する。すなわち、昇温速度を速めることによって反応を加速し、原料モノマーの濃度を速く終点濃度に到達させる。本実施例では、反応温度目標値の軌道｛ｔ，Ｔrset｝を、図１６（Ｂ）に示すような直線を多段につないだ折れ線とする。
【０２１３】
図１６（Ａ）は未反応の原料モノマーの濃度ｘmを示している。図１６（Ａ）、図１６（Ｂ）の横軸はバッチ経過時間ｔ、図１６（Ａ）の縦軸は濃度、図１６（Ｂ）の縦軸は温度である。図１６（Ｂ）では従来の平均的な反応温度目標値の軌道を点線１５０で示し、本実施例の反応温度目標値の軌道を実線１５１で示している。図１６（Ａ）の１６０は反応温度目標値が点線１５０の軌道を辿るときの原料モノマーの濃度ｘmを示し、１６１は反応温度目標値が実線１５１の軌道を辿るときの原料モノマーの濃度ｘmを示している。実線１５１の折れ点の座標を決めることにより、反応温度目標値の軌道を規定できる。
【０２１４】
［反応温度目標値軌道の決定法］
重合レートＲpとポリマー生産量ＰPは、原料モノマーの仕込量をＷmとしたとき、重合反応モデル１０により計算される原料モノマー濃度の推定値であるチルダｘmを用いて、次の式（８４）、式（８５）のように表すことができる。
【０２１５】
【数７７】
2018156152.tif 000082
【０２１６】
まず、バッチ終点近傍の時点ｔbにおける生産量を規定するモノマー濃度ｘmeを与える。そして、そのｔb時点の原料モノマー濃度の推定値チルダｘm（ｔb）が所望のモノマー濃度ｘmeにできるだけ一致するように、反応温度目標値の時間パターン｛ｔ（ｉ），Ｔrset（ｉ）｝を求める。上記のように反応速度定数推定部２ａにより開始剤の反応速度定数ｋiと原料モノマーの反応速度定数ｋmの基準状態値ｋm0と反応速度定数ｋmの温度依存係数ｂmとを推定することができれば、重合反応モデル１０によりｔb時点の未反応の原料モノマー濃度ｘm（ｔb）を推定することが可能である。
【０２１７】
反応温度目標値の時間パターン｛ｔ（ｉ），Ｔrset（ｉ）｝を求める方法は、次式の最適化問題として定式化できる。ここで、ｅ（ｔb）≡チルダｘm（ｔb）−ｘme、最適化変数φ≡［ｔ（ｉ）  Ｔrset（ｉ）］Tとし、制約条件として昇温速度の上限ｖrと目標値軌道座標の上限φuと目標値軌道座標の下限φlとを設定した。
【０２１８】
【数７８】
2018156152.tif 000083
【０２１９】
以上のようにして、反応温度目標値軌道決定部８は、反応温度目標値の時間パターン｛ｔ，Ｔrset｝を決定することができる（図１０ステップＳ１１）。
【０２２０】
［数値例］
表３の運転条件で説明した高機能性ポリマーについて、ポリマーの生産が完了するまでのバッチサイクル時間を短縮する反応温度目標値軌道を求めた。ｔb＝５．０ｈ、ｘme＝５．０ｗｔ．％のときの結果を図１７（Ａ）、図１７（Ｂ）、図１７（Ｃ）、図１７（Ｄ）に示す。図１７（Ａ）は原料モノマーの濃度ｘmを示し、図１７（Ｂ）は反応温度目標値Ｔrsetを示し、図１７（Ｃ）は生成されるポリマーの濃度ｘpを示し、図１７（Ｄ）はバッチサイクル時間の短縮時間を示している。図１７（Ａ）〜図１７（Ｃ）の横軸はバッチ経過時間ｔ［ｈ］、図１７（Ａ）、図１７（Ｃ）の縦軸は濃度［ｗｔ．％］、図１７（Ｂ）の縦軸は温度［℃］である。また、図１７（Ｄ）の横軸はポリマーの生産量［％］、図１７（Ｄ）の縦軸は時間［ｈ］である。
【０２２１】
図１７（Ｂ）では従来の平均的な反応温度目標値の軌道を点線１７０で示し、本実施例の反応温度目標値の軌道を実線１７１で示している。図１６（Ａ）の１８０は反応温度目標値が点線１７０の軌道を辿るときの原料モノマーの濃度ｘmを示し、１８１は反応温度目標値が実線１７１の軌道を辿るときの原料モノマーの濃度ｘmを示している。図１６（Ｃ）の１９０は反応温度目標値が点線１７０の軌道を辿るときのポリマーの濃度ｘpを示し、１９１は反応温度目標値が実線１７１の軌道を辿るときのポリマーの濃度ｘpを示している。
【０２２２】
図１７（Ｄ）によれば、従来のバッチサイクル時間８．０ｈに対して、ポリマーの生産量が９５％に達したときの短縮時間は０．５２ｈである。これは、図１７（Ｂ）の実線１７１で示した反応温度目標値の軌道を完全に実現できれば、０．５２／８．０×１００＝６．５％の生産性向上が期待できることを意味している。
【０２２３】
冷媒入口温度時間パターン設定部３とフィードフォワード制御実行部４と制御パラメータ調整部５とカスケード制御実行部６とモデル調整部７の動作は従来と同様である。ただし、本実施例では、モデル調整部７の調整対象が熱収支モデル１１でなく、熱収支モデル１１ａとなる。
【０２２４】
本実施例で説明した制御装置１１６ａは、ＣＰＵ（Central Processing Unit）、記憶装置及びインタフェースを備えたコンピュータと、これらのハードウェア資源を制御するプログラムによって実現することができる。本発明の制御方法を実現させるためのプログラムは、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、メモリカードなどの記録媒体に記録された状態で提供され、記憶装置に格納される。ＣＰＵは、記憶装置に格納されたプログラムに従って本実施例で説明した処理を実行する。
【産業上の利用可能性】
【０２２５】
本発明は、バッチ重合プロセス制御に適用することができる。
【符号の説明】
【０２２６】
１ａ…モデル記憶部、２ａ…反応速度定数推定部、３…冷媒入口温度時間パターン設定部、４…フィードフォワード制御実行部、５…制御パラメータ調整部、６…カスケード制御実行部、７…モデル調整部、８…反応温度目標値軌道決定部、１０…重合反応モデル、１１ａ…熱収支モデル、２０…反応発熱量推定部、２１…ポリマー濃度推定部、２２…モノマー濃度推定部、２３…反応速度定数・温度依存係数推定部、３０…シミュレーション実行部、３１…時系列テーブル記憶部、５０…プロセス動特性モデル作成部、５１…ゲインスケジューリング部、７０…実績データ記憶部、７１…平滑化部、７２…２乗誤差算出部、７３…非線形最適化部、７４…モデルパラメータ調整部、１１６ａ…制御装置。
【図１】
2018156152.tif 000084
【図２】
2018156152.tif 000085
【図３】
2018156152.tif 000086
【図４】
2018156152.tif 000087
【図５】
2018156152.tif 000088
【図６】
2018156152.tif 000089
【図７】
2018156152.tif 000090
【図８】
2018156152.tif 000091
【図９】
2018156152.tif 000092
【図１０】
2018156152.tif 000093
【図１１】
2018156152.tif 000094
【図１２】
2018156152.tif 000095
【図１３】
2018156152.tif 000096
【図１４】
2018156152.tif 000097
【図１５】
2018156152.tif 000098
【図１６】
2018156152.tif 000099
【図１７】
2018156152.tif 000100