化学産業における最適化技術の適用
藤田薫，江本 源一，竹下 聡彦，佐中 俊哉

石油化学およびその川下に位置する樹脂・樹脂加エなどを含めた化学産業において，プラントオペレーション，生産計画・スケジューリング，物流計画などに線形・非線形計画法や混合整数計画法などの数理計画法や制約プログラミング手法などの最適化技術がサプライチェーンの効率的運用のためのコア技術として通用されている．本稿では，三菱化学㈱における事例をもとに，化学産業における最適化技術の適用について紹介する．

キーワード：サプライチェーンマネジメント，APS，生産計画，スケジューリング，物流，最適化技術，制約プログラミング
Illlll…………ll………l…lll川Illllll‖‖＝‖‖‖＝‖‖＝‖‖‖‖‖＝‖‖‖‖‖‖＝‖‖‖‖‖‖‖……lllll……lll刷Ill川l‖＝州Il川……lll………ll………lll………ll……ll川Ill……llllllllll………lll………lllll川…………lllll…l州Illllllllllll川……l…

′一【、＼発が重要な課題であるとともに，日々変化する需要に1．はじめに対して柔軟に対応できる効率的なサプライチェーンマ

化学産業は，川上においては原油を精製して得られネジメントの仕組みが不可欠となってきている．

るナフサを主原料としたエチレン，プロピレン，ポリ

このように化学産業全体を伸撤すると，川上およびエチレンなどの石油化学製品を製造する大規模な石油川下を含め，また，そのサプライヤやカスタマまでも化学コンプレックス，川下においては機能性樹月旨，樹含めたグローバルな視点でのサプライチェーン全体の脂加工品などを製造する多数の多品種小量生産プラン効率的運用がますます重要な課題となってきている．
本稿では，総合化学メーカである三菱化学㈱におけトから構成されている（図1）．

石油化学コンプレックスは，分解，反応，圧縮，蒸る適用事例をもとに，化学産業における最適化技術適留などの単位操作をもとに多様な製品を製造する大型用の実際について概説する．紙面の都合上，それぞれの多数のプラントで構成されており，最大効率で生産の詳細は参考文献を参照されたい．

活動を行うための適切な原料調達計画，生産計画，お2．化学産業における最適化技術の適用事

よび物流計画の最適化が要求されている．

例

末端消費者の近くに位置する川下においては，昨今の消費者ニーズの多様化に応えたスピーディな商品開2．1全体概要
オペレーションズ・リサーチにおける代表的な最適′′「、＼

石油化学コンプレックス化技術は，化学産業における様々な分野において活用

樹脂・樹脂加エほかされている．すなわち，プラントオペレーションのリ

石油精製 石油化学

ジスティクスにおける計画の最適化，他ネ土とのアライ．鴎．

畠車㌔

アルタイム最適化，生産計画・スケジューリング，ロアンスや設備統廃合などの経営意思決定支援，サプライチェーン最適化などのビジネス分野に活用されている．また，反応メカニズムの推定や反応条件の最適化，エチレン
プロピレン
ポリエチレンほか

包装資材

機能性樹脂
フイルム・シートほか

図1化学産業における流通

ゲノム創薬などの研究開発分野にも通用されているが，本分野については別紙に譲ることとする．
図2に三菱化学㈱におけるビジネス分野への最適化
技術通用の変遷を示す．プラントオペレーションの領域においては，石油化学プラントを対象に，90年代
ふじた かおる，えもと げんいち，たけした としひこ，
さなか としや

前半には，LP（Linear Programming）による多変三菱化学㈱技術・生産センター技術部

数制御システムと連携した反応器や蒸留塔などのユニ〒712−8054倉敷市潮通3−10

ットプロセスの最適化がなされるようになったが，90

2003年8月号© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

（3）5498085

90009505

図2 最適化技術適用の変遷

年代後半には，NLP（Non−1inear Programming）

図3 化学産業における最適化システムの概念図

ベースのプラントワイドのリアルタイム最適化技術
RTO（Rea卜time Optimization）が適用され，プラテム）において行われる．比較的大型の連続プラント

ント全系にわたりオンラインで最適化が実現されるよにおいては，ユニットプロセスの安定化のための多変うになった．石油化学コンプレックスにおける生産計数制御，およびプラントワイドに利益を最大化するリ画立案支援システムとしては，80年代からLPによアルタイム最適化制御が適用されている．バッチプラる最適化が行われていたが，主要なプラント群ごとのントなどの品種切替のある多品種製造プラントでは，部分最適化であった．90年代には，全プラントを網生産スケジューリングシステムによるスケジュール立羅し，かつ高精度で高機能の大規模なMILP（Mixed案支援が行われている．本社と連携した生産計画シス

Integer Linear Programming）による最適化が実現

テムが構築され，事業所全体の生産バランス最適化がされた．2000年には，このようなシステムをベース実施されている．また，物流の領域においては，日々としたマルチサイトの生産計画最適化も実現された．の出荷計画だけでなく，戦略，計画レベルの物流拠点非線形プロセスに対しては，MINLP（Mixed Inte−配置，輸送計画の最適化を支援するシステムが構築さgerNonLlinear Programming）の技術も利用されてれている．このように，本社，事業所，各プラントのいる．ロジスティクスの領域では，MIP（Mixed各階層において，戦略，計画，オペレーションのそれIntegerProgramming）を用いた生産・物流拠点配置

ぞれの意思決定を支援する垂直統合がなされつつ，かの最適化，日々の業務支援である配車計画などに活用

つサプライヤからカスタマまでのサプライチェーンをされている．化学産業における多品種製造プロセスの意識した水平統合がなされている．

生産スケジュールの立案は，その複雑さから，ほとん

次節より，プラントオペレーション，生産計画，多どが熟練者による手作業に依存していたが，近年のス

品種プラントの生産スケジューリング，物流計画の最ピード経営の要請により，化学産業に適した生産計適化に関する適用事例について述べる．

画・スケジューリングシステムを開発し，順次適用を開始した．いわゆる先進的スケジューリング技術
APS（AdvancedPlanningandScheduling）［1，2］が2．2 プラントオペレーションのリアルタイム最適化

プラントオペレーションの最適化の事例として，石適用され始めた．さらに，企業間のコラボレーション，油化学コンプレックスにおける電気・蒸気を供給するバーチャル統合などを意識したサプライチェーンの最発電プラントのリアルタイム最適化について紹介する．
図4に電気・蒸気を供給する発電プラントのプロセ適化なども実現されようとしている．
上述のような最適化技術を通用した化学産業におけスフローの一例を示す．本発電プラントは，複数のポるビジネスプロセス最適化システムの概念図を図3にイラ（1B，2B，3B，・・・，nB），複数のタービン（1示す［3］．

T，2T，3T，・・・，mB）からなり，蒸気圧カレベル

各プラントレベルにおいては，プロセスの状態量が秒単位で収集され自動制御がDCS（分散型制御シス
550（4）

は，ボイラ主蒸気（SPl）のほか，各プラントに複数の圧力レベルの蒸気（SP2，SP3，SP4，SP5）© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

オペレーションズ・リサーチ

SP2

図4 発電プラントのプロセスフロー

図5 石油化学コンプレックスのモテルを供給している．
このような発電プラントにおける最適化システムで／【、＼、は，制約条件（2）式を満たす範囲で目的関数（1）式を最2．3 生産計画の最適化小化する問題に定式化し，非線形計画法の代表的手法

石油化学コンプレックスにおいては，図5に示すよである逐次2次計画法（SQP：SuccessiveQuadratic

うに，多数のプラントが複雑に組み合わされ，留分やProgramming）［4］を用いて最通解を求めている．

エネルギの授受が行われているため，プロダクトミッ

Minimize f（z）

（1）

Subjectto h（z）＝0，ZL≦z≦zU（2）

三菱化学㈱水島事業所においては，1，200品目以上

目的関数′（z）は，主としてボイラの燃料コストかにおよぶ原料，燃料，ユーティリティ，および製品のらなる発電プラント全体の運転コストを表現している．

バランスを最適化し事業所全体として利益を最大化す

発電プラントのポイラ，タービンなどの特性は，物質るための生産計画最適化システムを構築し，1998年収支，熱収支などの厳密な非線形の物理モデルで表さより運用している［10］．
本システムは，年間および月間の生産計画策定に利れ，（2）式の等式制約の形で表現される．等式制約式の数は約10，000である．ここで，Zはプロセス変数，用されているのに加え，最適な生産バランスの計画値

zL，㍑Uはそれぞれプロセス変数の下限値，上限値であを各プラントに毎日オンラインで提供している．また，る．

原料やユーティリティの供給制約が発生した場合など

例えば，タービンの発電量モデルは（3）式にて表現さの突発生産対応，プラント新増設検討などケーススタれる．タービン効率は（4）式にて定義される．

r＼クスの観点からの適切な生産管理が不可欠である．

ディとしての機能も重視し開発した．

E＝ヴ（H。−Hs。。t）F（3）

生産管理上，非線形特性を無視できない場合がある

り＝（Hn−H云。t）／（H。−Hsent）（4）

ため，生産量に応じて非線形に変化する生産特性を区

ここで，E：発電量，符：タービン効率，ダ：蒸気分的に線形化してモデル化した．また，プラントやタ流量，茸n：タービン入口エンタルピ，筏ut：タービービンの起動／停止の判断も最適化の対象とするため，

ン出口エンタルビ，Hs。。t：等エントロピ変化した場

混合整数型線形計画問題（MILP）として定式化した．

合のタービン出口エンタルビ．

大規模最適化問題を効率的に解くために，モデルの定本システムにより，蒸気，電気の需要変化に対応し式化や解法に工夫を加えた結果，モデルの規模は1年て，複数のポイラ，タービンの負荷配分の最適化が約間（365日間）のマルチピリオドの生産計画問題で，10分周期で行われ，15秒周期で制御演算される下位式の数で約780，000，整数（0−1）変数の数で約2，800の多変数制御システムと協調し，オンラインかつリアとなりCPU450MHzのパソコンで約1時間と実用ルタイムで最適化制御を行うことで燃料使用量の削減的な時間で最適解を得ることができる．米国GAMSなど大幅な省エネを実現した［5，6］．

Development社の最適化ソフトGAMSによりモデル発電プラントにおける事例を紹介したが，その他の主要化学プラントにおいても，このようなリアルタイムの最適化制御が行われている［7〜9］．

化し最適化ソルバとしてCPLEXを採用した［11］．
加えて，各プラントの年間の定期修理計画，言い換えると稼動計画を立案する機能も開発した．石油化学
プラントは，法令等により定期的に点検・修理を行う

2003年8月号© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

（5）551

必要がある．各プラントの定期修理期間，すなわち停止期間は点検や工事に必要な期間であるが，多数のプ本システムは，石油化学部門の生産管理の効率化，高度化に不可欠な情報システムとなっている．

ラントの定期修理が集中しすぎると工事要員の制約上，2．4 多品種プラントの生産スケジューリング実現不可能となる場合や工事コストが増大する場合が機能性樹脂や樹脂加工品などの複数品種を製造する発生する．逆に，定期修理が離散しすぎると石油化学プラントの生産管理においては，顧客からのオーダやコンプレックス内の生産バランスが取れなくなる．し販売予測を基に生産・物流コストを考慮し，生産拠点たがって，生産バランスの制約と要員の制約を同時にを決定する生産計画問題とともに，各製品を日々どの考慮した上で，各プラントの定期修理開始日を意志決設備で，どのようなシーケンスで生産すべきかの生産定変数として，生産コストと工事要員コストからなるスケジューリングの良否が当該事業の生産性を左右す総コストを最小化する最適化問題として定式化した．

る重要な要素となっている．

ここでも，プラントの稼動／停止の変数，定期修理開

そこで，このような多品種プラントのニーズに対応始日を表現する整数（0−1）変数を導入することによすべく，化学産業に見られるバッチ反応条件，品種切り，混合整数型線形計画問題（MILP）として定式化替条件，要員制約など複雑な制約条件を考慮し，かつした．

最適化を指向する生産スケジューリングシステムを開

4ヶ月間の定期修理計画の最適化問題は，式の数で約400，000，整数（0−1）変数の数で約2，000の数理発した［12］．
生産スケジューリングの問題は，一般的には膨大なモデルとなった．計算時間は，CPU450MHzのパソ組合せ最適化問題となるが，今回のシステムは様々なコンで，約90分で解くことができる．図6に最適な制約条件を満足する実行可能なスケジュールを即座に定期修理計画の結果の一例を示す．このような大規模立案すること，またユーザが希望する評価基準に基づな石油化学コンプレックスにおける定期修理計画の事

き最適化されたスケジュールを立案することを重視し，例は，これまで世界的にも報告例がない．

スケジューリングのエンジンとして制約プログラミン

図7は本生産計画最適化システムの構成図である．

グの手法を用い開発した．制約プログラミングは，各本システムにおいては，原料，製品バランスなどの入変数間の制約伝播によりあらかじめ探索すべき組合せ

出力関係は，パソコンのWindowsベースのGUIでの数を制約条件などの情報をもとに効率的に絞り込む

定義することができ，各プラントの生産特性を表現す手法であるが，本手法を利用することにより短時間でる生産量と原料，ユーティリティなどの関係式および実行可能解および近似最適解を求めることが可能とな

そのパラメータは，ネットワークを介して事業所内のった．制約プログラミングのためのC＋＋ソフトウエア

各製造課のパソコンから設定することができる．また，コンポーネントとして仏国ILOG社のILOG Solver

最適化された計画の計算結果を参照することもできる．

およびILOGSchedulerを利用した［13］．
上述の生産スケジューリングシステムを，樹脂製造プラントにおいて，原料を調合する調合工程と，製品を製造する製品製造工程の2工程を対象として適用し蔭哺﹂∧小ト

図6 最適定期修理計画552（6）図7 生産計画最適化システム© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

オペレーションズ・リサーチ

′■、、＼

た事例について以下に述べる．調合工程は，ブレンダ

2．5 物流計画の最適化に配合すべき原料を投入して混合する工程であり，製化学製品の出荷計画として，サイロからタンクロー品製造工程は6系列の設備からなり，調合された原料リへのバルク品の樟込み計画，サイロから取り出されを処理して製品を連続運転で製造する工程であるが品る製品の包装計画，および配車計画などがあり，これ種切替が必要である．調合工程は昼間勤務の作業要員らの業務を支援するスケジューリングシステムを開発で対応しており，製品製造工程はシフト勤務の作業員し1999年より運用している［14］．

で対応し24時間連続運転となっている．また品種切例えば，タンクローリの配車計画においては，品質替作業には作業要員が必要であるが，要員数の制約か管理の観点から前荷，後荷の製品品種の組合せでの洗ら同時に複数系列の切替作業はできない場合がある．

浄の度合いや禁止の条件，車両の運行状況など全てのまた，切替時間は切替前後の品種の組合せに依存する制約を考慮し，日々受注されるオーダに対して車両をため，製品の納期や在庫の制約を満たした上で，より効率的に迅速に割り付けている．

良いスケジュールを立案することが要求される．製造

これらの問題は，制約を満たし，かつコスト最小とする製品の種類が約400，原料の種類が1製品に対しなるように整数計画問題として定式化した．本システて平均約15，オーダ数が約100，スケジューリング対

ムは，CPU700MHzのパソコンで，約1分で最適解

象期間を1ヶ月間とし1時間刻みでスケジューリングを導出することが可能であり，日々の出荷業務に活用計算を実施した場合，CPUIGHzのパソコンで，10されている．

一方，戦略，計画レベルの問題として，図9に示す秒以内で計算できる．
図8は上段に製品毎に製造する設備をあらかじめ決ような物流に係わる輸送費（串云送費，配送費），保管定してスケジューリングしていた従来の方法と，下段費，荷役費の総コストが最小となるように，製造プラに今回開発したスケジューリングシステムを利用し，ントの拠点，ストックポイントなどの物流拠点の配置

製品毎の代替設備を考慮した上でスケジュール期間のを最適化する問題，いわゆる物流ネットワークの設計最小化を評佃基準として最適化した結果を比較した図問題にも取り組んだ．

である．この場合，スケジュール期間が約15％短縮

このような物流合理化検討を効率的に行うために物されることからスループット向上に有用であることが示された．

なお，飛び込みオーダなどの突発に即座に対応できるよう，使いやすいエーザインターフェイスを開発し，計画立案者の手垂加こよるスケジュールの微調整なども容易に行えるようにしたこと，また，品種の追加や切′′、＼替時間の変更などに伴うデータベースのメンテナンスを容易にしたことも実用上重要なポイントである．

図9 物流ネットワークの設計

ル期廟

条件設定

物流ネットワーク最適化結果

弘ざ鶉轍添

図8 スケジューリング結果
2003年8月号

図10 物流ネットワーク最適化システム© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

（7）553

流ネットワーク最適化システムを開発した［15］．
本システムにおいては，上述の総コストが最小とな

http：／／www，i−emu．COm／を参照されたい．

るように数理計画法を用いて物流拠点を最適化してい参考文献

る．拠点配置の判断を伴うため，整数（0−1）変数を［1］票田充，村松健児：生産スケジューリング，朝倉書店，導入し整数計画法（MIP）を採用した．

また，図10に示すようなエーザインタフエースを介して，条件設定および最適化計算の結果を確認することができる．物流拠点や顧客の住所から拠点間の距離を自動的に算出したり地理情報を可視化したりするための地理情報システム，各種輸送モードに対応する輸送コスト計算を行うための各種料金テーブルを組込んでいるため，効率的かつ高度な物流合理化検討が行えるようになった．

（2002）．

［2］西岡靖之：APS，日本プラントメンテナンス協会，（2001）．

［3］藤田薫：石油化学プラントにおける最適化技術の適用，2001年度日本オペレーションズ・リサーチ学会秋季研究発表会アブストラクト集，（2001），pp．30−31．
［4］茨木俊秀，福島雅夫：最適化の手法，共立出版，（1996）．

［5］Emoto，G．，A．Tsuda，T．Takeshita，M．T．Moni
S．Nakagawa and K．Fujita：IntegratedlargeLSCale
multivariablecontrolandreal．timeoptimizationofa

3．おわりにpowerplant，InProceedingsofIEEEConf．OnControI

化学産業には連続プロセス，バッチプロセス，加工Applications，（1998）．

プロセスなど多様なプロセスが存在するが，それぞれ［6］藤軋江本，竹下，津田，小野寺，中川：石油化学工場ににオペレーションから計画，学研各のレベルにわたり，おける発電プラントのエネルギー最適化制御，計装，43−

様々な意思決定問題が存在する．本稿では，これらの5，（2000），pp．34−37．

問題に数理計画法や制約プログラミングなどの最適化技術がサプライチェーンマネジメントのコア技術として活用されていることを概説した．これらの事例は最適化技術だけで成し得たものではなく，情報システム［7］小河，江本：化学プラントにおけるプラントワイド最適化制御の実現，計測と制御，35−10，（1996），pp．788−792．
［8］高津春雄編著：プロセス制御，コロナ社（2003）．

［9］伊藤利昭編著：化学産業における制御，コロナ社，（2002）．

部門や利用部門，ソフトウエアベンダ，大学・研究機［10］竹下，朝倉，M．Turkay，藤田：大規模石油化学工場に関等との連携による幅広い技術の融合が重要であるこおける生産計画最適化システム，2000年度日本オペレーとは言うまでもない．

ションズ・

化学産業には多数の企業がサプライチェーンを構成している．今後は，企業内にとどまらず，企業間のコラボレーションによる，より高度なサプライチェーンの最適化などの新しいビジネスモデルの構築を目指し挑戦していきたい．本稿が，オペレーションズ・リサーチに携わる研究者，技術者の一肋となれば幸いであ

リサーチ学会秋季研究発表会アブストラクト集，（2000），pp．216−217．

［11］BrookeA．，D．KendrickandA．Meeraus：GAMSauser sguide，（1998）．
［12］竹下，藤田，佐中，松川：化学プロセスにおける生産ス
ケジュー
リングシステムの開発，スケジューリング・シンポジウム2002講演論文集，（2002），pp．56−61．

［13］ILOGSA，ILOGOptimizationSuiteWhitePaper，る．

（2001）．

［14］佐中，西森，黒田，R．Raman，藤田：製品サイロ出荷ス
謝辞 本稿中の生産スケジューリングシステムの開発に協力いただいた松川公司氏（三菱化学エンジニアリング㈱），および物流ネットワーク最適化システムのケジュール最適化，2000年度日本オペレーションズ・リ
サーチ学会秋季研究発表会アブストラグト集，（2000），
pp．98−99．

開発に協力いただいた小塩隆之氏（三菱化学物流㈱）［15］佐中，小塩：物流ネットワーク最適化システムの開に感謝します．なお，生産スケジューリングシステム

発，第11匝ほ十測自動制御学会中国支部学術講演会論文

は登録商標Dr．Plannerとして1）リースされているの集，（2002），pp．152−153．

で，詳細は関連のホームページ

554（8）© 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

オペレーションズ・リサーチ