JPWO2010058584A1 - 配合及び配船計画作成システム、方法及びプログラム - Google Patents

Abstract

このシステムは、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；を作成する配合及び配船計画作成システムであって：前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成部と；作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成部と；作成された前記配船計画を格納するデータベース部と；を備える。

Description

本発明は、複数銘柄の配合原材料（以下、単に原材料とも称する）を入荷して混合する配合計画及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を作成するのに好適な配合及び配船計画作成システム、方法及びプログラムに関する。
本願は、２００８年１１月２１日に、日本に出願された特願２００８−２９８５９６号に基づき優先権を主張し、その内容をここに援用する。

鉄鋼を始めとする多くの産業においては、契約した購入目標量に沿い購入した種々の性状を有する多種類の原材料を混合して、この混合後の原材料を使用して鉄鋼等の製品、半製品を製造する。この際、鉄鋼等の製品、半製品の性状を一定範囲内に納めつつ生産量を満足し、かつ、各原材料が製造現場で在庫切れを発生させないように、輸送計画（船の場合は配船計画と呼ぶ）及び配合計画を立てることが求められている。そして、配船計画及び配合計画を作成するに際しては、費用を重要な指標として計画の良し悪しが判断される。このため、輸送費用（船を雇う費用であるフレートと、船が港で契約期間以上に停泊した場合に支払う滞船料）、及び原材料の購入費用のミニマム化（最小化）が求められる。

しかし、原材料を荷揚げ（入荷）する船舶の配船計画が定まらないと、揚港に荷揚げされる原材料の銘柄と量が決定されない。この結果、原材料を荷揚げされた製鉄所で、どの銘柄をどれだけの量使えるかを決定できない、つまり配合計画を作成できない。一方、船舶の配船計画を作成するためには、原材料を荷揚げする揚港である製鉄所で、どの銘柄をどれだけの量使うかが決まっている必要がある。このように、配合計画と配船計画とはお互いに深く相関し合っているため、そのいずれかを独立して計画することが非常に難しい。

関連する技術として、特許文献１に開示の原料輸送配船計画用推論装置では、原材料の使用計画及び原料船の年間稼動計画を既知データとして読み込んだ後に、在庫切れを起こしそうな銘柄を優先的に船と紐付ける。この工程を繰り返すことで配船計画が作成される。

また、特許文献２に開示の物流計画作成装置では、シミュレート部で制約条件に基づいて各輸送手段の運行をシミュレートし、原料在庫量推移算出部で上記シミュレート結果に基づいて原料銘柄毎の在庫量の推移を算出する。この計算の結果が評価値算出部で評価される。

日本国特開平８−２７２４０２号公報 日本国特開平１１−３１０３１３号公報

しかしながら、特許文献１及び特許文献２はいずれも、原材料の使用計画及び原料船の年間稼動計画を既知データとして読み込んで配船計画を作成する。つまり、配合計画と配船計画とを相互に連係させて作成することは考慮されていない。

また、現実の輸送に際しては、異なる種別の傭船契約に基づいて運用される船舶、例えば連続航海船、不定期船、スポット船が使用されるが、特許文献１、２ではそういった船舶の種類までも含めた輸送費用のミニマム化は、考慮されていない。

更に、連続航海船及び不定期船といった契約では、対象船舶を必ず（最優先で）配船をする必要がある。一方、スポット船は、引取目標量、在庫状況に応じて適切に雇う船の大きさ、船数を決定する必要がある。しかし、特許文献１、２ではそういった傭船契約種別を考慮した船の船団構成までを考慮した送費用のミニマム化は、行われていない。つまり、例えば、１つの連続航海船と最大積載量７５０００トン（ｔｏｎ、ｔ）の２つのスポット船で配船するのが良いのか、或いは１つの連続航海船と最大積載量５００００トンの３つのスポット船で配船するのが良いのか、などの選択肢を考慮した配船計画の方法が、実際の業務では必要となる。雇う船の最大積載量、航路によって船の輸送費用は大きく変わってくるが、これらを考慮した配船は、上記文献では行われていない。

更に、配合計画で輸送費用を考慮せずに計画を作成した場合には、必要以上に輸送費用が高い原材料を使用する配合計画が作成される畏れがある。このような計画方法の場合、配合計画決定後にどのように輸送を工夫しても、輸送費用を最適化することは困難である。例えば、性状がほぼ同一の原材料Ｘ，Ｙがあり、揚港（製鉄所）Ａ，Ｂでは原材料Ｘ，Ｙどちらでの使用も可能な場合を考慮する。揚港Ａに原材料Ｘを輸送する費用が２０＄／トン、原材料Ｙを輸送する費用４０＄／トン、揚港Ｂに原材料Ｘを輸送する費用が４０＄／トン、原材料Ｙを輸送する費用２０＄／トンであるとする。この場合、本来揚港Ａで原材料Ｘ、揚港Ｂで原材料Ｙを使用する計画を立てる方が輸送費用の観点からより有利な計画である。しかし、輸送費用が考慮されていない計画方法を用いた場合、揚港Ａで原材料Ｙ、揚港Ｂで原材料Ｘを使用する計画が作成される畏れがある。
加えて、実操業においては、使用予定銘柄の在庫状況が厳しい場合には、性状の近い銘柄（性状が規定した範囲に含まれる銘柄；一定の化学性質を共通して備える銘柄；互いに置き換えても使用可能な銘柄）を代替として使用することで、在庫切れの抑止が行われている。また、この代替使用を積極的に行うことで、フレートが高い船でしか輸送できない銘柄に変わり、性状の近い銘柄（性状が規定した範囲に含まれる銘柄）でフレートのより安い船で手配できる銘柄を輸送することで、輸送費用の削減を行っている。

このように、費用を考慮する場合には、配合原材料の購入費用はもちろん、その輸送費用も考慮して、計画を作成することが求められる。しかし、特許文献１及び特許文献２はいずれも、原材料の使用計画の段階において、上記事情までは考慮していない。

本発明は以上のような状況に鑑みてなされたものであり、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を、相互に連係させて一括して作成できるようにすることを一つの目的とする。また、船舶の契約種類までも含めて輸送費用のミニマム化が考慮できるようにすることを別の一つの目的とする。
更に、性状の近い銘柄（性状が規定した範囲に含まれる銘柄）を考慮することで、銘柄毎に個別に考慮するのに比べて更なる在庫切れの抑制と輸送費用のミニマム化を可能にすることを目的とする。

（１） 上記目的を達成するための、本発明の一態様にかかるシステムは、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；を作成する配合及び配船計画作成システムであって：前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成部と；作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成部と；作成された前記配船計画を格納するデータベース部と；を備える配合及び配船計画作成システムである。
（２） 上記（１）の配合及び配船計画作成システムでは、前記配合計画の作成と、前記配船計画の作成と、において、前記引取目標量および前記配合原材料の在庫状況に関して共通の制約条件が用いられてもよい。
（３） 上記（１）の配合及び配船計画作成システムにおいて、前記第１の配合計画作成部は：前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込むデータ取込み部と；前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する数式モデル設定部と；設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートするシミュレータと；前記シミュレータによるシミュレーション結果である配合計画を出力する出力部と；を備えてもよい。
（４） 上記（３）の配合及び配船計画作成システムにおいて、前記最適化計算部が、前記配合原材料の入荷量と前記引取目標量との関係に関して予め構築された目的関数に更に基づいて最適化計算を行ってもよい。
（５） 上記（３）または（４）の配合及び配船計画作成システムにおいて、前記データ取込み部により取り込まれる前記輸送費用には；船舶別・積港別・揚港別のフレートの情報と；前記銘柄別・揚港別フレートの情報と；が含まれてもよい。
（６） 上記（１）の配合及び配船計画作成システムにおいて、前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、前記配船計画作成部は：前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；前記船舶財源リストに含まれる前記船舶の運航制約、前記揚地での前記配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；を備えてもよい。
（７） 上記（１）の配合及び配船計画作成システムにおいて、前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、前記配船計画作成部は：前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われたグループ化配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；前記船舶財源リストに含まれる前記船舶の運航制約、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；を備えてもよい。
（８） 上記（１）（３）（７）のいずれかの配合及び配船計画作成システムは、 作成された前記配船計画に基づいて、第２の配合計画を作成する第２の配合計画作成部を更に備えてもよく、前記データベース部は前記第２の配合計画を更に格納してもよい。
（９） 上記（８）の配合及び配船計画作成システムにおいて、前記第２の配合計画作成部は：作成された前記配船計画に基づく前記配合原材料の入荷予定、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込む第２のデータ取込み部と；前記入荷予定を用い、前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する第２の数式モデル設定部と；設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う第２の最適化計算部と；前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートする第２のシミュレータと；前記シミュレータによるシミュレーション結果である配合計画を出力する第２の出力部と；を備えてもよい。
（１０） 上記（９）の配合及び配船計画作成システムでは、前記第１の配合計画の前記需給バランス制約と、前記第２の配合計画の前記需給バランス制約と、において、需要に関して同一の制約条件を設定し；前記第１の配合計画の前記需給バランス制約において、供給に関して、引取目標量を制約条件として設定し；前記第２の配合計画の前記需給バランス制約において、供給に関して、船舶による原材料の入荷量を制約条件として設定してもよい。
（１１） 上記（９）の配合及び配船計画作成システムでは、前記第１の配合計画の前記目的関数において、銘柄別・揚港別見做しフレートを用い；前記第２の配合計画の前記目的関数において、船舶別・積港別・揚港別フレートを用いてもよい。
（１２） 上記（１）または（３）の配合及び配船計画作成システムでは、作成された前記配合計画において、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われてもよい。
（１３） 上記目的を達成するための、本発明の別の一態様にかかる方法は、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；を作成する配合及び配船計画作成方法であって：前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成工程と；作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成工程と；作成された前記配船計画をデータベースに格納するデータベース格納工程と；を備える。
（１４） 上記目的を達成するための、本発明の別の一態様にかかるプログラムは、上記（１）の配合及び配船計画作成システムの各部分としてコンピュータを機能させるためのプログラムである。

本発明の各態様によれば、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を、相互に連係させ、一括して作成することができる。特に、傭船契約の種別の異なる船舶を考慮して、輸送費用に関して構築された目的関数を用意し、この目的関数に基づいて最適化計算を行うことができる。このことにより、船舶の傭船契約の種類（連続航海船、不定期船、スポット船）、船団の構成を決める船舶を雇うか、雇わないか（船財源リストの作成）、までも含めた、輸送費用のミニマム化のための計画立案が可能になる。更には、配合計画段階においても、輸送費用をミニマム化することを考慮した計画の作成が可能となる。
加えて、性状の近い銘柄（性状が規定した範囲に含まれる銘柄）を考慮することで、銘柄毎に個別に考慮するのに比べて更なる在庫切れの抑制と輸送費用のミニマム化のための最適化が可能になる。

本発明の各実施形態に係る配合及び配船計画作成システムの概略構成を示す図である。 第１の、または第２の配合計画作成装置を含むシステム構成例を示す図である。 第１の、または第２の配合計画作成装置の基本的な構成を示すブロック図である。 第１の、または第２の配合計画作成装置の構成を示すブロック図である。 第１の、または第２の配合計画作成装置による配合計画作成処理を説明するためのフローチャートである。 配合計画作成の概要を説明するための図である。 利用するフレートを設定するテーブルの例を示す図である。 入荷量が引取目標量から一定幅以上離れないという制約を説明するための図である。 入荷量が引取目標量から一定幅以上離れないという制約を説明するための図である。 配合計画作成の手順を説明するための図である。 非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に代えて線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）を導入したときの処理を示すフローチャートである。 配合計画を旬毎に作成した例を示す図である。 配船計画作成装置による配船計画作成処理を説明するためのフローチャートである。 取り込みデータのうちの船舶リストを説明するための図である。 取り込みデータのうちの船舶運航状況を説明するための図である。 取り込みデータのうちのフレートリストを説明するための図である。 取り込みデータのうちの滞船料のリストを説明するための図である。 船舶の選択処理を説明するためのフローチャートである。 抽出した連続航海船について計画作成期間における積地と揚地の組み合わせのパターンを作成している様子を示す図である。 スポット船で補うべき原材料輸送量を説明するための図である。 スポット船の航路リストを説明するための図である。 時刻と在庫量との関係を示す図である。 引取量が引取目標量から一定幅以上離れないという制約を説明するための図である。 引取量が引取目標量から一定幅以上離れないという制約を説明するための図である。 マクロ最適化とミクロ最適化との関係を模式的に示した図である。 積地における負荷の平準化を目的とする目的関数について説明するための図である。 揚地における負荷の平準化を目的とする目的関数について説明するための図である。 配船計画の例を示す図である。 本発明の配合計画装置、配船計画装置として機能しうるコンピュータ装置のハードウェア構成例を示す図である。

以下、添付図面を参照して、本発明の好適な実施形態について説明する。本実施形態では、複数の製鉄所に、世界中に点在する山元（積地）から鉱石や石炭等の配合原材料を船舶による輸送で入荷し、それらを混合する例を説明する。すなわち、この例では、複数銘柄の配合原材料が、複数の供給元から船舶にて複数の供給先に輸送され、複数の供給先に入荷する。そして、この配合原材料が、各供給先において配合されて使用される。この際に、この複数銘柄の配合原材料の配合計画、及びこの複数銘柄の配合原材料は、積地である複数の供給元から揚地である複数の供給先に船舶によって輸送される。この船舶による輸送の際の配船計画を作成するのに好適な配合及び配船計画作成について、以下説明する。

（システム構成）
図１は、本発明の実施形態１から４に係る配合及び配船計画作成システムの概略構成を示す図である。なお、図１には例えば第２の配合計画作成装置３００が示されているが、実施形態１においては、これは使用されない。以下、複数の実施形態に概ね共通する構成、工程について、その概要を説明する。その後に、各実施形態に含まれる要素、工程について、詳細に記載する。
１００は第１の配合計画作成装置であり、原材料の引取目標量に基づいて、複数銘柄の原材料を入荷して混合する配合計画を作成する。ここでは、製鉄所毎に原材料を配合する配合計画を作成する。これにより、各製鉄所において原材料を日々どれだけ使用していくかを示す、使用予定量が計画される。この第１の配合計画作成装置１００が、本発明でいう第１の配合計画作成部として機能するものである。

２００は配船計画作成装置であり、第１の配合計画作成装置１００により作成された配合計画に基づいて、複数銘柄の原材料（鉱石や石炭等）を複数の積地（世界中に点在する山元）から複数の揚地（製鉄所）に輸送する配船計画を作成する。本実施形態では、製鉄所毎の輸送費用平準化ではなく、全製鉄所合計での輸送費用をミニマム化する配船計画を作成することを目的としている。更には、輸送費用に加えて、原材料の購入費用を含めた費用をミニマム化する配船計画を作成することを目的としている。この配船計画作成装置２００が、本発明でいう配船計画作成部として機能するものである。

３００は第２の配合計画作成装置であり、配船計画作成装置２００により作成された配船計画に基づいて、複数銘柄の原材料を入荷して混合する配合計画を作成する。この第２の配合計画作成装置３００が、本発明でいう第２の配合計画作成部として機能するものである。

４００はデータベースであり、各装置１００〜３００が計画を作成する上で使用するデータや、各装置１００〜３００が作成した計画をコンピュータ５００が参照可能、更には修正可能な方式で格納する。これにより、各装置１００〜３００が計画を作成する上で使用するデータや、各装置１００〜３００が作成した計画の一括管理が可能となり、最新情報の共有化を図ることができる。

５００はコンピュータであり、データベース４００に格納されたデータを参照、更新したり、データベース４００にデータを格納したりする。なお、ここではコンピュータ５００を一つしか図示していないが、本実施形態では、複数台のコンピュータがＬＡＮやインターネットを介して接続される。コンピュータ５００には、例えばプロセスコンピュータ等と称される上位コンピュータや、各所（製鉄所、本社、船会社、山元等）に設置されたデータベース４００にアクセス可能なコンピュータ端末等がある。

なお、図１に示すシステム構成は一例に過ぎず、それに限定されるものではない。例えば各装置１００〜３００がそれぞれ一つの装置として構成されるかのように図示したが、装置１００〜３００がそれぞれ複数の機器からなるコンピュータシステムで構成されるようにしてもかまわない。また、例えば第１の配合計画作成装置１００及び第２の配合計画作成装置３００が同様のアルゴリズムで配合計画を作成するような場合には、一つのコンピュータシステムが第１の配合計画作成装置１００及び第２の配合計画作成装置３００として機能するように構成してもよい。更には、本発明でいう第１の配合計画作成部、配船計画作成部、データベース部、第２の配合計画作成部が一つの装置として実現される場合も本発明の範疇にあるものとする。

（第１の実施形態）
（第１の配合計画作成装置１００）
第１の配合計画作成装置１００では、原材料の引取目標量に基づいて、複数銘柄の原材料を入荷して混合する配合計画を作成する。引取目標量とは、山元（積地）別、銘柄別の引取目標量（引取予定量）である。各山元とは、銘柄毎に、例えば年間といった期間中で、どれだけの量を引き取るかについて契約している。この期間中の引取量を、相当する期間に含まれる旬の数で割れば、旬毎の引取目標量が得られる。なお、本明細書において、旬は月を３つに分割した期間の単位を指す。

図２は、第１の配合計画作成装置１００を含むシステム構成例を示す図である。図２に示すように、第１の配合計画作成装置１００は、配合計画を作成するに際して、操業者によるデータ設定を受け付けるか、或いはデータベース４００からデータを取り込む。この際、取り込まれるデータには、例えば以下のものが含まれる：配合計画を立案する上で必要となる計画作成期間、引取目標量、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、船舶を利用する際の輸送費用情報を含む制約条件、及び前提条件。原材料の性状に関するデータは、例えば、石炭化度、各成分含有比率など、物理的、化学的な性状を含む。

第１の配合計画作成装置１００は、多種類（複数銘柄）の原材料を入荷して混合する混合計画を、シミュレーションを実行して作成する。第１の配合計画作成装置１００は、原材料の需給バランス制約、混合後の性状制約を満たすように、配合計画として各銘柄の使用量（配合割合）、入荷量を求める。詳細は後述するが、第１の配合計画作成装置１００は、ＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法に則って構築された、原材料の需給バランス制約を表す数式モデル（「需給バランスモデル」とも称する）、及び、混合後の性状制約を表す数式モデル（「性状モデル」とも称する）を設定することにより配合計画の最適化を図る。
数式モデルの設定とは、数式モデルの展開とも呼ばれる一連の工程であり、以下に説明するような工程を含み、本実施形態の装置の各部、あるいは方法によって行われる。本実施形態において、数式モデルは、船舶数や港数などの各条件の変化に対応できるように、抽象的な形式で予め構築・定式化されている。この数式モデルに対して、各配列の添え字の最大数（例えば船舶数を表す）や、式中の係数及び定数の値などを、計画立案の条件に沿って具体的に定める。

表示部３０３では、第１の配合計画作成装置１００で求められた各銘柄の使用量（配合割合）、入荷量、在庫推移グラフ、各種帳票を表示する。

操業者評価部３０４では、求められた配合計画を様々な観点（例えば、在庫推移、性状等）から操業者が評価し、満足のいく結果でなければ、必要に応じて配合割合等を修正する。その際に、必要に応じて目的関数の重みや評価の指標を変えたり、数式モデルを設定する対象期間・計画確定期間を変えたりする。また、全部の或いは指定した処理のみ使用量の固定をする等、操業者の意志を反映させた数式モデルの設定が可能である。そして、第１の配合計画作成装置１００で再度配合計画を作成し直す。

図３は、第１の配合計画作成装置の基本的な構成を示すブロック図である。図３に示すように、第１の配合計画作成装置１００は、シミュレータ（在庫推移シミュレータ３１１、性状シミュレータ３１２を含む）、本発明でいう数式モデル設定部として機能するモデル設定部（需給バランスモデル設定部３１３、性状モデル設定部３１４を含む）、本発明でいう最適化計算部として機能する計画部３１５を含んで構成され、更に入出力部を併せ持つ。

在庫推移シミュレータ３１１は、各原材料の需給状態（在庫推移）を計算するシミュレータである。性状シミュレータ３１２は、原材料を混合した後の性状を計算するシミュレータである。在庫推移シミュレータ３１１、性状シミュレータ３１２が互いに連動することで、原材料の在庫推移、混合後の性状を計算する。

本実施形態においては、例えば以下の情報を含む、入力データ３１６に基づいて数式モデルの設定処理を行う：配合計画を立案する上で必要となる計画作成期間、引取目標量、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、及び船舶を利用する際の輸送費用情報等。配合計画の立案開始日時から予め設定された最適化期間分を対象として、予め設定した時間精度に基づいて、ＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理整数計画法等に則り、需給バランスモデル設定部３１３にて需給バランス制約（在庫制約）を表す数式モデルが設定され、性状モデル設定部３１４により性状制約を表す数式モデルが設定される。

需給バランスモデル設定部３１３、性状モデル設定部３１４により設定された数式モデルを用いて、在庫を切らさないようにするとともに要求される性状を満足し、かつ、費用（原材料の購入費用及び輸送費用）をミニマム化して配合計画を作成するように、計画部３１５が最適化計算を行う。この最適化計算の結果に基づいて、在庫推移シミュレータ３１１、性状シミュレータ３１２に対する計算指示を作成する。この計算指示を受けて、在庫推移シミュレータ３１１が、在庫推移をシミュレートし、性状シミュレータ３１２が計画に従って製造される製品・半製品の性状をシミュレートする。例えば、鉄鋼業においては、石炭を混合して焼き固めたコークス（製品）、鉄鉱石を還元して得られる銑鉄を精錬した溶鋼を凝固させたスラブ（半製品）等の性状をシミュレートする。

かかる第１の配合計画作成装置１００によれば、従来のように予め決められたルールに基づいて計算指示が行われるのではなく、計画部３１５により行われた最適化計算の結果に基づいた計算指示を在庫推移シミュレータ３１１、性状シミュレータ３１２に出力する。このため、そのときの事象に応じた最適な計算指示を確実に行うことが可能となる。

例えば、図７に示すように予め設定された計画確定期間分について、在庫推移シミュレータ３１１、性状シミュレータ３１２によるシミュレーションが、行われる。このシミュレーションが終了すると、立案開始日が更新され、更新前の計画確定期間の最終状態、つまり更新後の立案開始日での在庫推移、性状の情報に基づいて、新たな最適化期間分の需給バランスモデル設定部３１３により在庫制約を表す数式モデルが設定され、性状モデル設定部３１４により性状制約を表す数式モデルが設定される。設定後の数式モデルは、計画部３１５に与えられる。このように、計画部３１５に、在庫推移、性状の情報に基づいて設定された数式モデルが与えられると、これを用いて、最適化計算が実行される。

以上のように、各シミュレータ（在庫推移シミュレータ３１１、性状シミュレータ３１２）と、各モデル設定部（需給バランスモデル設定部３１３、性状モデル設定部３１４）と、計画部３１５とを連動させた詳細なシミュレーションを実行することで、最適な配合計画を作成することができる。すなわち、本実施形態において行われるシミュレーションは、従来のような所定のルールに基づくシミュレーションではなく、最適化計算の結果に基づいて行われる。このため、１回のシミュレーションを実行するだけでも、理論的な最適解を確実に得ることが可能となる。この構成により、従来のようにシミュレーション結果を評価してシミュレーションを何回も繰り返して実行する必要がなく、シミュレーション結果３１７を迅速に、かつ、高精度に作成することができる。したがって、配合計画を作成する対象が大規模であっても実用時間内に作成することが十分に可能である。

また、従来法では、計画作成期間が長くなると、計算対象となる期間が長くなり、これにつれて、問題規模が急激に大きくなるため、求解が不可能になる問題があった。しかし、本手法では全計画作成期間を、これより短い最適化期間に分割することで、問題規模を小さくすることができる。このため、全体としての計画作成期間が長くなっても、問題を解くことが可能となる。上述のようにして得られたシミュレーション結果３１７を配合計画として出力する。

また、需給バランスモデル設定部３１３、性状モデル設定部３１４により設定するモデルの規模が非常に大きい場合や、制約条件が非常に多くて複雑な場合でも、在庫推移シミュレータ３１１、性状シミュレータ３１２に記載された需給バランス制約、性状制約のうち、配合計画作成に影響が大きく、重要な部分のみを抽出し、この抽出された部分のみを、需給バランスモデル設定部３１３、性状モデル設定部３１４に取り込んでもよい。これによって、需給バランスモデル設定部３１３、性状モデル設定部３１４の数式モデルの規模を適切な範囲にして、実用的な時間内で最適化計算を行うことができる。在庫推移シミュレータ３１１、性状シミュレータ３１２は、考慮すべき需給バランス制約、性状制約を全て記載することができるので、１回のシミュレーションを実行して作成された配合計画は現実に実行可能となることが保証される。

上述したように、本実施形態においては、シミュレータ（在庫推移シミュレータ３１１、性状シミュレータ３１２）とモデル設定部（需給バランスモデル設定部３１３、性状モデル設定部３１４）と計画部３１５とを連動させて配合計画を作成する。このため、以下の効果が得られる：（１）シミュレーションを繰り返して実行せずに配合計画を作成することができる；（２）配合計画作成に影響が大きい重要な部分の制約のみを計画部３１５に取り込むことで、計算時間を短縮することができる；（３）大規模な問題を解くことが可能になる。

以下、第１の配合計画作成装置１００の構成及びこの装置を用いて実行する配合計画作成方法の各ステップをより詳細に説明する。図４は、図３を用いて説明した第１の配合計画作成装置１００の基本的な構成に対する、配合計画作成装置の詳細な構成を示す図である。また、図５は、この装置を用いて実行する配合計画作成方法の各ステップを示すフローチャートである。

配合計画作成の概要として、例えば図６に示すように、以下のような計算、調整の工程が含まれる：複数ある製鉄所（揚港）ａ〜ｃでの原材料（銘柄）の需給バランスを取ること；（各銘柄Ａ〜Ｎの在庫を切らさない等）要求される性状を満足させること；かつ、費用（原材料の購入費用及び輸送費用）をミニマム化すること；前記の条件を満たすような配合計画として、製鉄所ａ〜ｃ毎の各銘柄Ａ〜Ｎの使用量（配合割合）、入荷量を決定すること。ここで、製鉄所毎の使用量の合計量である予定使用量は、入力データとして与えられる。このため、配合割合（％）＝使用量／予定使用量×１００となる。このため、使用量、配合割合の一方が決定されれば、他方が決定されることとなる。

（１）入力データの取り込み（図４の入力データ取込み部３５１、図５のステップＳ３０１）
本処理に必要な情報（引取目標量、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、船舶を利用する際の輸送費用情報等）をオンラインにて読み込み、必要に応じて操業者が修正を加える。

ここで、入力データ取込み部３５１により取り込まれる引取目標量は、山元（積地）別、銘柄別の引取目標量（引取予定量）を表す情報である。各山元とは原材料の銘柄毎に、例えば年間等の期間において、どれだけの量の原材料を引き取るか（引取量）について契約している。この期間全体の引取量を月数で割れば、月毎の引取目標量が得られる。この引取目標量に近づけるように、入荷することが求められる。しかし、年間で数万トン程度の引取量の上下へのぶれは山元との交渉により、許容範囲内となる。また、契約によっては、所定の銘柄については所定の期間は引取しないといった契約も考えられる。そういった例外的な条件を、計算に含めてもよい。例えば、原材料Ａのある月での引取目標量が５万トンの場合で、引取目標量からの上下へのぶれが年間で６万トン（月当たり５千トン）の場合を考える。この場合、当該月の入荷予定について、上限（入荷予定量上限）５万トン＋５千トン、下限（入荷予定量下限）５万トン−５千トンの間の範囲が、許容範囲となる。

原材料の在庫状況は、計画作成期間の初日における製鉄所別、銘柄別の在庫量（トン数）を表す情報である。

原材料の性状は、原材料毎の成分等の性状を表す情報である。例えば、原材料である鉄鉱石の性状としては、Fe₂O₃、Fe₃O₄、SiO₂、Al₂O₃等の原材料についての、性状情報が含まれる。

原材料の購入費用情報は、山元（積地）別、銘柄別の原材料の単価（＄／トン（ｔｏｎ、ｔ））を表す情報である。

船舶を利用する際の輸送費用情報は、図１５、１６に示すように、船舶リストにリストアップされている船舶を利用する場合のフレート、及び、船舶リストにリストアップされている船舶を利用する場合の積揚港別の滞船料、を表す情報である。また、船舶を利用する際の輸送費用情報には、銘柄別・揚港（揚地）別見做しフレートを表す情報も含まれる。輸送費用は、本来上述した船舶別・積港別・揚港別フレートにより、一意的に定まるものである。しかし、第１の配合計画作成装置１００の作業段階においては、入荷に関して、原材料を積載する船舶が決定されていない。これら積載する船舶が未決定の原材料に対して、原材料の輸送費用を見積もるために、銘柄別・揚港別のフレートが必要となる。ここで、銘柄別・揚港別フレートは、原材料を積載する船舶の選択等によって本来はフレートが異なるため、一意に決定することができない。そこで、銘柄別・揚港別フレートに替え、銘柄別・揚港別に、概算である見做しフレートの情報を取得する。銘柄別・揚港別見做しフレートとしては、例えば経験等に基づいて設定した銘柄別・揚港別フレート、或いは過去実績から統計的手法、例えば銘柄別・揚港別に過去実績のフレートを収集し、その平均値を銘柄別・揚港別フレートとして見做した銘柄別・揚港別見做しフレートが予めリストアップされている。

以上説明した入力データ取込み部３５１及びステップＳ３０１が、本発明でいう第１の配合計画作成部のデータ取込み部及びそれによる処理の例である。

（２）配合計画作成期間の設定（図４の計画作成期間設定部３５２、図５のステップＳ３０２）
配合計画を作成する期間を設定する。この作成期間は立案者の必要に応じて任意の期間を設定可能とする。ここでは、一例として１０日間分を立案する。

（３）配合計画作成時間精度の設定（図４の時間精度設定部３５３、図５のステップＳ３０３）
配合計画を作成する時間精度並びにシミュレーション精度を設定する。この時間精度並びにシミュレーション精度は、立案者の必要に応じて個別に任意の精度を設定可能とする。例えば立案の細かな精度を必要とする計画作成期間の前半では精度を細かくし、粗い計画で十分な計画作成期間の後半では精度を粗くすることで、十分な精度と短い計算時間での、効率的な計画作成が可能になる。

（４）最適化期間の設定（図４の最適化期間設定部３５４、図５のステップＳ３０４）
配合計画を作成する最適化期間を設定する。この最適化期間は立案者の必要に応じて個別に任意の対象期間を設定可能とする。ここでは、一例として、計画作成期間を通して、最適化期間を３日間とする。

（５）計画確定期間の設定（図４の計画確定期間設定部３５５、図５のステップＳ３０５）
配合計画を確定する計画確定期間を設定する。この計画確定期間は、立案者の必要に応じて個別に任意の期間を設定可能とする。例えば、立案の細かな精度を必要とする計画作成期間の前半では計画確定期間を短くし、粗い計画で十分な計画作成期間の後半では計画確定期間を長くする。このことで、十分な精度を持ちながら短い計算時間で、効率的な計画作成が可能になる。ここでは、一例として、計画確定期間を１日に設定する。この場合は、数式モデルに対する解に基づいてシミュレートした結果得られる配合計画に対しては、計画作成期間を通して最初の１日分を確定する。

（６）配合計画の需給バランス制約を数式モデルに設定（図３の基本的構成図の需給バランスモデル設定部３１３、図４の需給バランスモデル設定部３５６、図５のステップＳ３０６）
入力データ取込み部３５１により取り込まれたデータの全部又は一部に基づいて、設定された最適化期間分を、設定された時間精度で、需給バランス制約を数式モデルに対して設定する。

各銘柄の使用量を表す変数を下記の（式１）に示すように定義する。また、銘柄の在庫量を表す変数を下記の（式２）に示すように定義する。また、各銘柄の入荷量を表す変数を下記の（式３）に示すように定義する。なお、式中の「所」は、本実施形態では揚地と対応する製鉄所を示す。

需給情報を基に設定した数式モデル、つまり需給バランス制約モデルを以下に示す。ここでは、各銘柄の在庫量が、一定の安全在庫量と呼ばれる値以上であることが要求される。この場合の制約は、下記の（式４）と表される。

また、各銘柄の在庫量は、前日の在庫量、前日の入荷量、前日の使用量より決定される。この場合の関係を表す制約式は、下記の（式５）と表される。つまり、当日の在庫量は、前日の在庫量と当日に入荷（荷揚）する量を足した値から、当日の使用量を引いた値となる。

また、各銘柄の使用量のある日の合計は、当該日の全銘柄合計に対して予定された使用量と一致する必要がある。この場合の関係を表す制約式は、下記の（式６）と表される。

また、各種原材料の購買に対する要因等に鑑み、操業者は、目標とする配合割合を設定し、この与えた目標とする配合割合に近い配合割合を実現する配合計画が作成されることを求める。つまり配合割合が操業者の想定と大きくかけ離れると、想定した購買量を満たせなくなったり、購買量を越えたり、また操業設備に無理な操業を及ぼすことが想定される。このため、目標として与えた配合割合に近い配合割合が出力されることが必要となる。上記機能を実現するための制約を以下に示す。つまり、銘柄の使用量から使用目標量（目標とする配合割合）（定数）を引いた値を、使用目標量からの溢れ量の変数として定義する。ここで、使用量と使用目標量が近い量を取る程良い計画であるため、この溢れ量は少ない程良い。上記理由のため、後述するように、この溢れ量が、目的関数の項目として追加され、最適化によってミニマム化される。同様に銘柄の使用目標量から使用量を引いた値を、使用目標量からの不足量の変数として定義する。ここで、使用量と使用目標量は近い量を取る程良い計画であるため、この不足量は少ない程良い。上記理由のため、後述するようにこの不足は、目的関数の項目として追加され、最適化によってミニマム化される。この場合、各銘柄の使用量、使用目標量、溢れ量、不足量との関係を表す制約式は下記の（式７）と表される。つまり、溢れが生じる場合は使用量から溢れ量を引き、不足が生じる場合は不足量を足すと、使用目標量と一致する。

更に、前日の配合割合とその翌日の配合割合が大きく乖離すると、操業に困難を来たす。つまり、別原材料を使用するための段取り時間の増加や、設備の故障の原因となる。このため、前日の配合割合とその翌日の配合割合が大きく乖離することがない配合計画が求められる。上記機能を実現するため、銘柄の当該日の使用量と前日の使用量との差について、上限量を表す変数を下記の（式８）に示すように定義する。

この変数を用いて上記を実現するための制約を下記の（式９）に示す。つまり、銘柄の当該日の使用量から当該日前日の使用を引いた値は、当該日の使用量と当該日前日の使用量との差以下とする。ここで、当該日の使用量と当該日前日の使用量とが近い量を取る程、好適な計画である。上記理由のため、後述するように、この溢れ量は、目的関数の項目として追加され、最適化によってミニマム化される。同様に、銘柄の当該日前日の使用量から当該日の使用を引いた値に関しても、制約式として定式化する。

また、各銘柄の入荷量は、入荷予定量として与えられた量の範囲内に入っていることが要求される。この場合の関係を表す制約式は、下記の（式１０）、（式１１）と表される。つまり、当該月に入荷する入荷量の合計は当該月の入荷予定量上限以下、入荷予定量下限以上になる必要がある。

また、第１の配合計画作成装置１００では、引取目標量に基づいて、配合計画を作成する。後述するように、最適化（ステップＳ１０３〜Ｓ１０６）及びシミュレーション（Ｓ１０７）を含む一連の工程は、複数ループ反復して実行できる。初回ループでは、ステップＳ３０１で取り込んだデータに基づいて、後述する数式モデルの設定を行い、次回ループ以降では、シミュレータでのシミュレーション結果を反映させて、数式モデルの設定を行う。引取目標量制約に関しての数式モデル設定では、例えば以下の条件を考慮する：最適化する引取量（入荷量）が引取目標量から一定幅以上離れないこと；及び、引取の可否（前述したように所定の銘柄については所定の期間は引取しないといった事情もありうる）。ここで、引取量が引取目標量から一定幅以上離れないという制約を定式化する方法として、例えば図８Ａに示すように、単に旬毎（或いは月毎）の引取目標量に対して、上下の許容範囲を設定し、入荷量がその許容範囲に含まれることを制約条件とする方法、または構成が考えられる。しかし、その場合、例えば入荷量が許容範囲の下限を満たしているが旬毎の引取目標量を下回る状況が続いたような場合、年間の累積入荷量は、年間の引取目標量を下回ることもありうる。そこで、図８Ｂに示すように、旬毎（或いは月毎）に立案開始日から該当する旬までの引取目標量累積及び引取量累積を算出し、引取目標量累積と引取量累積との差を小さくする（最小とする、上下限値を越えないようにする等）という制約を設定するのが好適である。上記制約式を定式化するために、旬毎の引取目量標累積からの溢れ量、不足量の変数を定義する。

また、各旬の引取量累積の変数を定義する。

まず、各銘柄の引取量累積を表す制約式は下記の（式１２）と表される。つまり、引取量累積は、立案開始日から当該旬までの期間に入荷する量の合計となる。

各銘柄の引取目標量累積と溢れ量、不足量との関係を表す制約式は下記の（式１３）と表される。つまり、引取累積量から溢れ量を引くか、或いは不足量を足すと引取目標累積量と一致する。ここで、引取累積量と引取目標累積量は近い量を取る程好適な計画であるため、この溢れ量、及び不足量は少ない程良い。上記理由のため、後述するように、この溢れ量、及び不足量は、目的関数の項目として追加され、最適化によってミニマム化される。

なお、上述した需給バランス制約は一例であり、他の制約に替えたり、他の制約を加えたりしてもよい。

（７）配合計画の性状制約を数式モデルに設定（図３の基本構成図の性状モデル設定部３１４、図４の線形化部３５７ａを含む性状モデル設定部３５７、図５のステップＳ３０７、Ｓ３０７ａ）
入力データ取込み部３５１により取り込まれたデータの全部又は一部に基づき、設定された最適化期間及び時間精度を用いて、性状制約を数式モデルに設定する。鉄鉱石の配合計画を作成する場合に用いられる原材料の性状としては、例えば、以下のものが挙げられる：鉄分、SiO₂、Al₂O₃、SiO₂、等。石炭の配合計画を作成する場合の性状としては、例えば、以下のものが挙げられる：ＣＳＲ（熱間反応後強度）、ＤＩ（コークス強度）、ＶＭ（揮発分）、膨張圧等。これらの性状が、要求される性状制約を満たす必要がある。混合後の性状モデルの一例を（式１４）に示した。なお、（式１４）では下限値Ｓを有する例を示すが、上限値を有する場合や、上限値及び下限値の両方を有する場合もありうる。
ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓ・・・（式１４）
ｘ_A〜ｘ_N：原材料（銘柄）Ａ〜Ｎの配合割合
Ｓ：下限値（定数）

ここで、多くの性状について、性状モデルに含まれる数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）は、下式（式１５）に示すように、配合割合に対して線形となる。
ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）
＝Ｗ_A×Ｘ_A＋Ｗ_B×Ｘ_B＋・・・＋Ｗ_N×Ｘ_N・・・（式１５）
Ｗ_A〜Ｗ_N：銘柄毎の銘柄iに含まれる当該成分の性状
例えば、SiO₂に関して、銘柄Ａの配合割合が４０％、銘柄ＡのSiO₂成分が１%、銘柄Ｂの配合割合が６０％、銘柄ＢのSiO₂成分が２%の条件で混合した場合、混合後のSiO₂成分の性状は、１×０．４＋２×０．６＝１．６％となる。

ところが、性状によっては、その性状を表す数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）が非線形となることがある。この場合、次に述べるように、線形化部３５７ａで、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に代えて線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）を導入して数式モデルを定式化する。

線形化部３５７ａでの処理について説明する。ある性状を表す数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）が非線形である場合、それに代えて線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）を導入する。この線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）は、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）の下限をなすもの、すなわち（式１６）の関係が成立するものを考える。なお、（式１６）は常に成立する必要はなく、必要な範囲で成立していればよい。
ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）・・・（式１６）

例えば線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）として、（式１７)に示す加重平均を考える。加重平均は、単一銘柄を１００％使用した場合の性状を非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）から求め、配合割合を乗算して、使用銘柄分足し合わせた値である。
加重平均製鉄所＝Σ[配合割合（＝使用量（製鉄所、銘柄）／使用量合計（製鉄所））×単一銘柄１００％時性状銘柄]・・・（式１７）

説明を簡単にするため、銘柄Ａの配合割合が９０％、銘柄Ｃの配合割合が１０％の例を考える。この場合、線形の数式ｆ´（９０、０、１０、・・・、０）となる加重平均は、下式で表される。
ｆ´（９０、０、１０、・・・、０）
＝０．９×ｆ（１００、０、・・・、０）＋０．１×ｆ（０、０、１００、・・・０）

過去の実績等からこの加重平均が（式１６）を満たすことがわかっていれば、この加重平均を線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）として利用することができる。すなわち、加重平均≧Ｓを制約とすれば、（式１４）が成立するものとみなすことで、定式化できる可能性が得られる。

線形化部３５７ａでは、（式１４）´に示すように、線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する下限値として、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する下限値Ｓよりも小さな仮下限値Ｓ´＝Ｓ−ｓ（ｓ：オフセット値）を設定して、数式モデルを設定する。
ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓ´・・・（式１４）´

以上は、混合後の性状制約が下限値を有する場合を例に説明した。なお、上述した性状制約は一例であり、他の制約に替えたり（混合後の性状制約が上限値を有する場合を含む）、他の制約を加えたりしてもよい。

以上説明した需給バランスモデル設定部３５６（図３の需給バランスモデル設定部３１３）及びステップＳ３０６、並びに、性状モデル設定部３５７（図３の性状モデル設定部３１４）及びステップＳ３０７、Ｓ３０７ａが、本発明でいう第１の配合計画作成部の数式モデル設定部及びそれによる処理例である。

（８）固定化抽出処理（図４の固定化抽出処理部３５８、図５のステップＳ３０８）
図９に示すように、配船計画の項目である積港、積銘柄、積量、揚港、揚銘柄、揚量のうち固定化されているもの、すなわち変更できないものを抽出する。予め与えられる条件によって、各傭船に対して「積港」及び「積銘柄」が固定化されている場合は、船舶別・積港別・揚港別フレート（図１５を参照）を用いる。つまり、これらの傭船は原材料を積載する傭船が決定されているため、船舶別・積港別・揚港別フレートを用いることで、最適化によって原材料を入荷する製鉄所が決定された時点で、正確な輸送費用計算が可能となる。

また、いずれの項目も固定化されていない場合や、「積港」だけが固定化されている場合は、銘柄別・揚港別見做しフレートを用いる。つまり、積港と積銘柄が決定されていない場合には、当該船舶に関する積港と積銘柄自体を変更可能にすることで、後述の最適化を用いて、当該積港と積銘柄より輸送費用の安い積港と積銘柄の有無を検討することが可能となる。この場合は、銘柄別・揚港別見做しフレートを用いることで当該船舶に関して、当該積港と積銘柄より輸送費用の安い積港と積銘柄に、当該船舶の積港と積銘柄を変更させた計画を後述する最適化により評価、検討させる。これにより輸送費用のより安い計画を作成することを可能としている。なお、同一の傭船に関しては、固定化が最もされていないレコードの状態をこの傭船の固定化状況と考える。この固定化抽出処理は、図５に示したタイミングである必要はなく、例えば配合計画作成を開始するときに行われるようにしてもよい。

ここで、本実施形態は、配船計画が存在しない場合、つまり立案開始日以降の配船計画が全く立てられておらず、最初から配合、配船計画を作成して行く場合を例にしてあるため、上記固定化の対象となる船舶情報は存在しない。しかし、本発明はこれを変形した実施形態にも適用できる。例えば、立案開始日から３ヶ月分の配合・配船計画を立案し、１月計画した後に、２ヶ月分は前回立案を参考にし、１ヶ月を全く新規に作成することで、新たな３ヶ月分の配合・配船計画を立案する、等のローリングを行いながら、配合・配船計画が立案されてもよい。このような運用においては、配合計画を立てる際には、計画作成期間の立案開始日に近い日付の一部に対して、配船計画が存在し、またその傭船の一部が固定化されており、計画作成期間の上記に対する残りの期間に対しては、配船計画が全く存在しない。この変形した実施形態において、配合計画を実施する際には、固定化対象となる船舶が存在することとなる。このような実施形態においては、上記固定化抽出処理を用いることで、既存の条件を尊重しつつ、配船及び配合計画の漸次向上を図ることができる。

（９）配合計画数式モデルを目的関数に基づいて最適化（図３の計画部３１５、図４の配合計画求解部３５９、図５のステップＳ３０９）
上記のように設定された線形及び整数制約式でなる需給バランスモデル、性状モデルを併せて配合計画数式モデルとし、予め構築した目的関数に基づきＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法により最適化問題として問題を解くことにより、最適な使用量、入荷量を計算する。

ここでは、目的関数に関して線形式を用いた場合の例を示す。本実施形態では、費用（原材料の購入費用及び輸送費用）のミニマム化を目的としており、目的関数Ｊの一例を（式１７）に示す。目的関数を用いて求解するに際して、購入費用情報及びステップＳ３０８において設定された輸送費用情報を用いる。

ここで、本実施形態は、第１の配合計画では、配船計画が存在しない場合、つまり立案開始日以降の配船計画が全く立てられておらず、最初から配合、配船計画を作成して行く場合を例にしてあるため、上記固定化の対象となる船舶情報は存在しない。このため、（式１７）において、船舶別・積港別・揚港別フレートに関する項（右辺第２項）は存在しないこととなる。

また、第１の配合計画作成装置１００では、原材料の入荷量と引取目標量との関係、すなわち最適化する原材料の入荷量が引取目標量から一定幅以上離れないことを目的としている。
例えば、旬単位の入荷量と旬毎の引取目標量との差のミニマム化を目的とした目的関数を構築する。或いは、引取量累積（入荷量累積）と引取目標量累積の差のミニマム化を目的とした目的関数を構築する。具体的には、各銘柄の入荷量を旬単位（或いは月単位）に集計し、それまでの累積を考える（入荷量累積）。また、各銘柄の引取目標量を旬単位（或いは月単位）に集計し、それまでの累積を目標値として設定する（引取目標量累積）。そして、入荷量累積と引取目標量累積の差のミニマム化を目的とするよう目的関数を構築する。つまり、引取目標量累積及び引取量累積を考え、引取目標量累積と引取量累積との差を小さくすることも目的としている。本実施例では、累積を考え、旬毎の引取目標累積量からの溢れ量、不足量の合計量をミニマム化する項目を目的関数に追加する。

なお、（式１７）、（式１８）は目的関数の一例であり、他の目的関数に替えたり、他の目的関数を加えたりしてもよい。

例えば、与えた目標とする配合割合に近い配合割合に配合計画を近づける必要があり、更に前日の配合割合とその翌日の配合割合が大きく乖離することがない配合計画を作成する必要がある場合は、使用目標量からの溢れ量、不足量、及び当該日の使用量と当該日前日の使用量との差をミニマム化する項目を目的関数に追加する。

以上の設定した式（数式モデル）を混合整数計画法にて解くことにより、需給バランスモデル、性状モデルを併せた配合計画数式モデルに対する最適解が得られる。

以上説明した配合計画求解部３５９（計画部３１５）及びステップＳ３０９が、本発明でいう第１の配合計画作成部の最適化計算部及びそれによる処理の例である。

（１０）最適化計算による求解結果の判定（図４の求解結果判定部３６０、図５のステップＳ３１０、Ｓ３１１）
（式１４）´を用いた最適化計算による求解結果が、非線形の数式を含む数式モデルｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓを満たすか否かを判定する。その結果、非線形の数式を含む数式モデルｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓを満たせば、この求解結果を、後述する性状シミュレータ３６２に対する計算指示としてシミュレーションを実行させる。上記求解結果が非線形の数式を含む数式モデルｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓを満たさなければ、線形の数式を含む数式モデルｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓ´を調整する（図５のステップＳ３１１）。具体的には、仮下限値Ｓ´を微増させる。

図１０は、ステップＳ３０７〜Ｓ３１０の処理、すなわち非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に代えて線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）を導入したときの処理を示すフローチャートである。ステップＳ４０１において、需給バランスモデル、性状モデル（非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に代えて線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）を導入して構築したもの）、目的関数Jに基づいて最適化計算を実行する。

この場合に、（式１４）´に示したように、線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する下限値として、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する下限値Ｓよりも小さな仮下限値Ｓ´＝Ｓ−ｓ（ｓ：オフセット値）を設定する。

次にステップＳ４０２において、線形の数式を含む数式モデルｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓ´を用いた最適化計算による求解結果が、非線形の数式を含む数式モデルｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）≧Ｓを満たすか否かを判定する。すなわち、ステップＳ４０１の最適化計算による求解結果（各銘柄Ａ〜Ｎの使用量（配合割合））を（式１４）に代入し、（式１４）が成立するか否かを判定する。

ステップＳ４０２の結果、（式１４）が成立すれば、本処理を終了する（図１０のステップＳ４１２に移行する）。それに対して、（式１４）が成立しなければ、ステップＳ４０３に進んで、仮下限値Ｓ´を予め設定された増減幅で微増させて、再度ステップＳ４０１の処理を実行する。すなわち、（式１４）が成立するまで、仮下限値Ｓ´を微増させて、最適化計算による求解を繰り返す収束計算を実行する。

なお、本実施形態では、混合後の性状制約が下限値を有する場合を例にして説明したが、上限値を有する場合も同様である。この場合、線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）は、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）の上限をなすものを考える。また、ステップＳ４０１では、線形の数式ｆ´（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する上限値として、非線形の数式ｆ（ｘ_A、ｘ_B、ｘ_C、・・・、ｘ_N）に対する上限値よりも大きな仮上限値を設定する。

（１１）求解した解に基づいて在庫推移をシミュレーション（図３の在庫推移シミュレータ３１１、図４の在庫推移シミュレータ３６１、図５のステップＳ３１２）
上記配合計画数式モデルに対する解、及び、入力データ取込み部３５１により取り込まれたデータの全部又は一部に基づいて、配合の全部或いは一部を対象として、設定した計画確定期間分について、設定した計画作成精度で、シミュレーションを実行する。このシミュレーションでは、配合計画数式モデルには組込むことができなかった制約条件、例えば一定の規則に基づかない条件など、定式化が難しいもの、及び、操業のルール等も組み込んでシミュレートする。このことで、配合計画数式モデルに対する求解結果として出された解を実操業で問題なく使用可能な配合計画に変更する。これにより、実操業で求められる時間精度と、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

また、数式モデルでは取扱うことが難しい制約の一例として、配合割合が変わった場合の設備の段取りに掛かる段取時間等をシミュレーションに取り込み、正確にシミュレートすることで、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

（１２）求解した解に基づいて性状をシミュレーション（図３の性状シミュレータ３１２、図４の性状シミュレータ３６２、図５のステップＳ３１３）
上記配合計画数式モデルに対する解、在庫推移シミュレータ３６１によりシミュレーションされた在庫推移、及び、入力データ取込み部３５１により取り込まれたデータの全部或いは一部に基づいて、配合の全部或いは一部を対象として、設定した計画確定期間分について、設定した計画作成精度で、性状のシミュレートをおこなう。シミュレーションの結果として、原材料の混合後の性状結果が得られる。このシミュレーションでは、配合計画数式モデルには組み込むことができなかった制約条件、操業のルール等も組み込んでシミュレートすることで、配合計画数式モデルに対する求解結果として出された解を実操業で問題なく使用可能な配合計画に変更する。これにより、実操業で求められる時間精度と、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

以上説明した在庫推移シミュレータ３６１（在庫推移シミュレータ３１１）及びステップＳ３１２、並びに、性状シミュレータ３６２（性状シミュレータ３１２）及びステップＳ３１３が、本発明でいう第１の配合計画作成部のシミュレータ及びそれによる処理の例である。

（１３）配合計画の確定（図４の確定部３６３、図５のステップＳ３１４）
上記在庫推移シミュレーション、性状シミュレーションにより導き出された配合計画のうちで設定した計画確定期間分を確定する。図７に示すように、本実施形態では計画確定期間を１日と設定しているので、作成した配合計画の最初の１日分を確定する。作成した配合計画のうちで上記計画確定期間に入らなかった部分については、その計画は確定せずに破棄する。

（１４）計画作成期間分、或いは計画確定期間分の計画が確定したか判定（図４の判定部３６４、図５のステップＳ３１５）
このステップの実行時点までに確定した計画確定期間が予め設定した計画作成期間の全体を含んでいるかを判断する。本実施形態では、計画作成期間が１０日間であるので、第１０ループで計画を確定した時点で計画確定期間分の計画が確定する。このため第１０ループで計画を確定終了した時点で１０日分の配合計画を作成して、処理を終了する。

（１５）立案開始日の更新（図４の更新部３６５、図５のステップＳ３１６）
このステップの実行時点で確定した計画確定期間が予め設定した計画作成期間の全体を含んでいない場合、上記配合計画のうちで確定した配合計画期間直後の日時を新たな立案開始日として設定する。本実施形態では、図７に示すように、第１ループでは当初１日目０時であった立案開始日を２日目０時に、第２ループでは当初２日目０時であった立案開始日を３日目０時に更新する。

（１６）配合計画の出力（図４の出力部３６６、図５のステップＳ３１７）
以上のようにして作成した配合計画は、出力部３６６により、表示部３０３に画面表示されたり、データベース４００を含む外部機器にデータ送信されたりする。

以上説明した出力部３６６及びステップＳ３１７が、本発明でいう第１の配合計画作成部の出力部及びそれによる処理の例である。

以上のように、現在の在庫推移状態に応じて、需給バランス制約、性状制約について、まず所定の最適化期間分について、計画作成時間精度で数式モデルを設定し、設定した配合計画数式モデルを目的関数に基づいて求解し、求解した解に基づいて、在庫推移、混合後の性状をシミュレートし、シミュレーション結果から求められた配合計画のうちで、設定した計画確定期間分を確定し、計画確定期間直後の日時を新たな立案開始日時とすることにより、新たな計画対象期間分の配合計画を確定する一連の処理を順次、予め定めた回数だけ、繰り返して実行する。このことで、所望する計画作成期間分の配合計画を作成することができる。これにより、任意の時間精度を必要とする配合計画を高速かつ詳細に最適化することができ、しかも得られた結果をそのままで実操業に適用できる。

なお、上記実施形態では、図１１に示すように、配合計画は一定の期間（例えば旬）毎に作成される。また、複数の性状α、βについて性状モデルが非線形となることがある。
図１１において、Ａは性状制約を満たしている（（式１４）が成立している）ことを、Ｂは性状制約を満たしていないことを意味する。すなわち、図１１の例では、性状αについて複数旬（４月上旬及び下旬）で性状違反が発生しており、同様に性状βについて複数旬（４月上旬及び下旬）で性状違反が発生している。

この場合に、各旬及び各性状について別個に図１０で説明した収束計算を行うと、以下の問題が生ずる。具体的にいえば、４月上旬で性状αについて収束計算を行い、続いて性状βについて収束計算を行い、また、４月下旬で性状αについて収束計算を行い、続いて性状βについて収束計算を行うのでは、計算処理に時間がかかってしまう。

そこで、対象の旬及び性状についてまとめて図１０で説明した収束計算を行うようにする。例えば４月上旬及び下旬で性状α、βについてまとめて収束計算を行う（図１０のステップＳ４０３で性状α、βの仮下限値の微増（或いは仮上限値の微減）を同時に行う）ことにより、高速化を図ることができる。

また、上記実施形態では、図１０のステップＳ４０３で仮下限値Ｓ´を微増（或いは仮上限値を微減）させた後、再度ステップＳ４０１の処理を実行すると説明した。この場合に、収束計算が変化しても設定に変化のない数式モデル、具体的には上述した需給バランスモデルや元々線形の性状モデルは保持しておく。そして、仮下限値を微増（或いは仮上限値を微減）させて再度ステップＳ４０１の処理を実行する場合に、収束計算に従って変化のある数式モデル、具体的には仮下限値を微増させた（或いは仮上限値を微減させた）数式モデルのみ変更するような仕組とすることにより、高速化を図ることができる。

また、配合計画（例えば使用量（配合割合））として、年次計画、期計画、月次計画といった長期間の計画を立案することが多い。このように長期の配合計画を予め作成し、その配合計画を基準の配合計画とし、第１の配合計画作成装置１００により作成されたより短期の配合計画が、基準となる配合計画から一定幅以上離れないようにすることも重要となる。

そこで、（式１７）に示したような費用（原材料の購入費用及び輸送費用）に関して構築された目的関数Ｊに加え、予め作成された基準となる配合計画と一定幅以上離れないようにする条件に関して構築された目的関数J´を用いるようにしてもよい。目的関数Ｊ´の一例を（式２０）に示す。
Ｊ´＝Σ（｜基準配合割合（銘柄）−配合割合（銘柄、日）｜）→ミニマム化・・・（式２０）
基準配合割合：基準となる配合計画における配合比

上記例では、月次計画において、期計画を基準となる配合計画として、日々の配合計画を作成する場合の一例を示した。この場合は、配合割合（銘柄、日）と基準配合割との差の銘柄毎、日毎に合計したものをミニマム化する。他の例として、期計画を立案する場合、年次計画を基準となる配合計画として計画を作成しても良い。この場合、月次計画では配合割合（銘柄、月）を決定するとした場合は、配合割合（銘柄、月）と基準配合割との差を銘柄毎、月毎に合計した値をミニマム化する。

なお、基準となる配合計画は、例えば過去の実績に基づいて作成され、その作成手法はどのようなものであってもよい。もちろん、本発明を適用した配合計画作成手法により長期間の計画を予め作成しておき、それを基準となる配合計画としてもよい。

（配船計画作成装置２００）
配船計画作成装置２００は、データベース４００から、例えば以下のデータを取り込む：原材料の使用予定量（第１の配合計画作成装置１００により作成された配合計画による使用量）、引取目標量、傭船契約の種別の異なる船舶がリストアップされた船舶リスト、船舶リストにリストアップされている船舶の運航状況、原材料の在庫状況、原材料の単価を表す購入費用情報、船舶リストにリストアップされている船舶を利用する場合の輸送費用情報。配船計画作成装置２００は、取り込んだデータに基づいて、例えば３ヶ月（９旬）分の配船計画を作成する。ここで、旬は月を３つに分割した期間の単位を指す。配船計画として、具体的には、連続航海船、不定期船、スポット船に対する積揚地（積揚港）、積揚銘柄、積揚量、寄港順、着岸バース、入出港タイミング、及び雇うべきスポット船の船数と船型（船の最大積載量を基に定義される船舶の大きさ）等を決定する。

ここで、使用予定量は、第１の配合計画作成装置１００により作成された配合計画による使用量である。使用予定量は、入荷量が設定された引取目標量にできるだけ近い量となり、尚且つ輸送費用を安くする条件を考慮して決定されている。このため、第１の配合計画作成装置１００により作成された配合計画に基づく使用量を、配船計画作成装置２００の入力データとして用いることで、配船計画を立てる際にも、輸送費用を考慮せず計画された使用量を用いるより、より安い輸送費用での配船計画の立案が可能となる。

例えば、性状がほぼ同一の原材料Ｘ，Ｙがあり、揚港（製鉄所）Ａ，Ｂでは原材料Ｘ，Ｙどちらでの使用も可能な場合を考える。揚港Ａに原材料Ｘを輸送する費用が２０＄／トン、原材料Ｙを輸送する費用４０＄／トンであり、揚港Ｂに原材料Ｘを輸送する費用が４０＄／トン、原材料Ｙを輸送する費用２０＄／トンである例を考慮する。この場合、配合計画を作成する際に輸送費用が考慮されていないと、揚港Ａで原材料Ｙ、揚港Ｂで原材料Ｘを使用する計画を立ててしまう畏れがある。しかし、輸送費用の観点からは、揚港Ａで原材料Ｘ、揚港Ｂで原材料Ｙを使用する計画がより好ましい。本実施形態では、上記のような状況を回避することが可能となる。

図１に示すように、配船計画作成装置２００において、２０１は船舶の運航、積地、揚地での設備、ヤード等を模擬したシミュレータである。シミュレータ２０１は、後述するマクロ最適化部２０２、ミクロ最適化部２０３により決定された入力情報を受け、これに基づいて詳細なシミュレーションを実行する。この入力情報は、連続航海船、不定期船、スポット船に対する積揚地（積揚港）、積揚銘柄、積揚量、寄港順、着岸バース、入出港タイミング、及び雇うべきスポット船の船数と船型（船の最大積載量を基に定義される船舶の大きさ）の情報を含む。本シミュレータは、在庫推移シミュレータ及び船舶運航状況推移シミュレータにより構成される。

在庫推移シミュレータは、各製鉄所における原材料の在庫推移を計算する。この在庫推移シミュレータでは、各製鉄所の原材料の使用予定量、船舶の原材料の銘柄毎の荷揚げ時刻を考慮し、詳細に原材料の銘柄毎の在庫推移を計算する。例えば、船舶に複数銘柄が積載され、１銘柄を荷揚げした後、２銘柄目を荷揚げする時点で、ヤード能力が溢れていた場合、ヤード上の原材料の在庫量が減り、ヤード能力に余裕が出来るまで時間を空けて荷揚げをする必要が生じる場合がある。このような事情等を考慮して、荷揚げ時刻に対応させて在庫の推移が正確にシミュレートされる。

船舶運航状況推移シミュレータは、積揚港の沖着日時（ＥＴＡ：Estimated Time Of Arrival）、積揚港着岸の日時（ＥＴＢ：Estimated Time Of Berthing）、積揚港出港の日時（ＥＴＤ：Estimated Time Of Departure）を含む船舶の運航状況の推移を計算する。この船舶運航状況推移シミュレータでは、荷積能力、荷揚能力の他に、他岸壁に船舶が存在するかどうか（存在する場合には着岸できない）等、他船舶との干渉等も考慮し、詳細に船舶の運航状況をシミュレートする。例えば、船舶の荷揚げに使用するアンローダの基数は、荷揚げする銘柄が積載されているハッチの位置、同一揚港の別バースで荷役している船舶があるか、ないか等を考慮する。この荷揚げに使用するアンローダの基数により荷揚能力が影響される。一例として、１基のアンローダで荷揚げする場合は、１５００ｔ／ｈで１００％能力で荷揚げを行える。また、２基のアンローダの場合は、１５００ｔ／ｈ×２基で７０％能力で荷揚げを行える。上記船舶運航状況推移シミュレータは、これらアンローダ基数等の条件による荷揚能力の変化を、シミュレーションに取り込み、正確にシミュレートする。このことで、実操業に求められる細かな制約まで考慮した、具体的な生産・物流計画の立案を可能とする。

２０２はマクロ最適化部であり、製鉄所の配合計画（上記の原材料の使用予定量）に支障をきたさないこと、及び、積み出し可能量を守ることを前提に、輸送費用のうちのフレートの合計金額を最も安価にすることを一つの目的として、以下の各条件を決定するように最適化を行う：連続航海船、不定期船、スポット船に対する積揚地（積揚港）、積揚銘柄、積揚量、寄港順、及び雇うべきスポット船の船数と船型（船の最大積載量を基に定義される船舶の大きさ）。マクロ最適化部２０２は、本発明でいう船舶財源リスト作成部として機能する船舶財源リスト作成部２０２ａ、本発明でいう数式モデル設定部として機能する数式モデル設定部２０２ｂ、本発明でいう最適化計算部として機能する最適化計算部２０２ｃを備え、例えば９旬分を旬精度に演算する。

２０３はミクロ最適化部であり、マクロ最適化部２０２により最適化された計画において滞船料の合計金額を最も安価にする着岸バース、入出港タイミングを決定するように最適化を行って、シミュレータ２０１に対する指示を算出する。ミクロ最適化部２０３は、本発明でいう数式モデル設定部として機能する数式モデル設定部２０３ａ、本発明でいう最適化計算部として機能する最適化計算部２０３ｂを備え、例えば１旬分を分精度に演算する。

２０４は本発明でいうデータ取込み部として機能するデータ取込み部であり、データベース４００から原材料の使用予定量、引取目標量、傭船契約の種別の異なる船舶がリストアップされた船舶リスト、船舶リストにリストアップされている船舶の運航状況、原材料の在庫状況、原材料の購入費用、船舶リストにリストアップされている船舶を利用する場合の輸送費用、積地での船舶停泊状況、荷積能力状況、設備修理・休止予定、揚地での船舶停泊状況、荷揚能力状況、設備修理・休止予定等のデータを取り込む。

２０５は本発明でいう出力部として機能する出力部であり、シミュレータ２０１によるシミュレーション結果として作成された配船計画、具体的には、連続航海船、不定期船、スポット船に対する積揚地（積揚港）、積揚銘柄、積揚量、寄港順、着岸バース、入出港タイミング、及び雇うべきスポット船の船数と船型（船の最大積載量を基に定義される船舶の大きさ）を画面表示したり、データベース４００を含む外部機器にデータ送信したりする。

以下、本実施形態に係る配船計画作成装置２００による配船計画作成処理の詳細を説明する。図１２は、配船計画作成装置２００を用いた配船計画作成方法における各処理のステップを説明するためのフローチャートである。本実施形態では、ユーザが設定した立案開始日から３ヶ月（９旬）を計画作成期間として配船計画を作成する。

（１）データの取込み（ステップＳ１０１）
配船計画作成装置２００のデータ取込み部２０４は、データベース４００から原材料の使用予定量、引取目標量、契約の種別の異なる船舶がリストアップされた船舶リスト、船舶リストにリストアップされている船舶の運航状況、原材料の在庫状況、原材料の購入費用、船舶リストにリストアップされている船舶を利用する場合の輸送費用、積地での船舶停泊状況、荷積能力状況、設備修理・休止予定、揚地での船舶停泊状況、荷揚能力状況、設備修理・休止予定等のデータを取り込む。

ここで、原材料の使用予定量は、第１の配合計画作成装置１００で作成された配合計画から算出される、計画作成期間における製鉄所（揚地）別、原材料の銘柄別の使用予定量を表す情報である。原材料はその銘柄毎に、品質・性状等に違いがあるため、それぞれ銘柄毎に使用予定量を決め、配合される。

引取目標量は、山元（積地）別、銘柄別の引取目標量（引取予定量）を表す情報である。各山元とは銘柄毎に例えば年間どれだけの量を引き取るかについて契約しており、それを旬数で割れば旬毎の引取目標量が得られる。この引取目標量に近づけるように、配船することが求められる。ただし、配船計画との関係で、引取目標量からの数万トン程度の上下へのぶれは山元との交渉により、許容範囲内となる。また、契約によっては、所定の銘柄については所定の期間は引取しないといった契約も考えられる。このような契約に関する具体的な情報をデータに含めてもよい。

船舶リストは、図１３に示すように、契約の種別の異なる船舶、ここでは具体的に連続航海船、不定期船、スポット船をリストアップした情報である。連続航海船の傭船契約は、契約期間において連続航海する契約である。このため、必ず（最優先で）配船することが求められる。不定期船の傭船契約は、契約期間において契約した航海数又は契約した航海期間のみ航海する契約である。このため、契約した航海数又は契約した航海期間内で必ず（最優先で）配船することが求められる。スポット船は、本実施形態の使用、または実行時点では未契約である。連続航海船及び不定期船を配船しても、目標とする引取量を満たせない、或いは揚地の在庫を充足できない場合に、これらスポット船にスポット的な航海を依頼することができる。連続航海船については、取込みデータ中に、傭船コード（船一隻一隻を特定するコード）、契約区分、契約期間（開始日及び終了日）、最大積載量、船名が記載される。不定期船については、傭船コード、契約区分、契約期間（開始日及び終了日）、契約の内容（契約した航海数又は契約した航海期間）、最大積載量、船名が記載される。これら連続航海船及び不定期船は船舶を個別にリストアップしているが、スポット船については、船舶の航行できる地域名と、船舶の船型（船の最大積載量を基に定義される船舶の大きさ）でリストアップし、傭船コード（地域名と大きさが記述される）、契約区分、最大積載量が記載される。

なお、船型とは、船の最大積載量を基に定義される船舶の大きさを表す。スポット船の船型を表すＰｍａｘはパナマ運河を通過できる船舶（一般にこの船型はパナマックスと呼ばれる）、Ｃａｐｅはケープ岬を通過できる船舶（一般にこの船型はケープサイズと呼ばれる）、ＶＬ（Very Large）はこれらより大きい大型船であることを意味する。ここで、通常パナマックスとは、長さ９００フィート以内、幅１０６フィート以内の船で、最大積載可能量が６万〜８万トンクラスの船を指す。また通常ケープサイズとは、最大積載か能力が１５万〜１７万トンクラスの船を指す。スポット船については、配船計画を立てる段階で、この航行できる地域と船舶の大きさを基に、必要な船数、船型（船の最大積載量を基に定義される船舶の大きさ）を決定する。配船計画がある程度確定される段階になって、実際の船会社と交渉して、上記船型にマッチする船を契約する手続きが取られる。このため、配船計画を立てる段階では、まず未契約の状態（船会社と交渉する前の段階）で、必要な船数、船型（船の最大積載量を基に定義される船舶の大きさ）を決定することが求められる。

船舶運航状況は、図１４に示すように、船舶リストにリストアップされている各船舶の運航状況の実績及び確定している予定を表す情報である。積荷から揚荷までを一つの航海として取り扱い、航海Ｎｏが付される。各航海について、積−積−揚、積−揚−揚のように、積港及び揚港は１港の場合も、複数港の場合もある。船舶リストにリストアップされている各船舶について、航海Ｎｏ．、積揚ドック区分、積揚連番、積揚港コード、バースコード、積揚港沖着の日時（ＥＴＡ）、積揚港着岸の日時（ＥＴＢ）、積揚港出港の日時（ＥＴＤ）、航海時間が記載される。

例えば連続航海船Ａの航海Ｎｏ．３とは、２００８年３月７日２０時に積港（Ｘ１港）沖に着き、２００８年３月１２日２０時に積港（Ｘ１港）のコード「１」で表されるバースに着岸し、２００８年３月１４日２０時に積港（Ｘ１港）を出港した後、４６９２０分航海して、２００８年４月１６日１０時に揚港（Ｂ港）沖に着き、２００８年４月１６日１３時に揚港（Ｂ港）のコード「１１」で表されるバースに着岸し、２００８年４月１８日１４時に揚港（Ｂ港）を出港する航海である。また、この船は契約区分が連続航海であるため、連続航海船Ａは、積−揚−積−揚の順に連続的に航海、寄港している。つまり連続航海船Ａは、航海Ｎｏ．２の最後の揚港（Ｄ港）を２００８年２月２２日９時に出港した後、２０８２０時間航海して、２００８年３月７日２０時に航海Ｎｏ．３の最初の積港（Ｘ１港）の沖に着いている。

立案開始日（配船計画を立案する対象期間の始めの日）が、立案を実行する日に対して、将来である場合は、第１の配合計画作成装置１００で作成された配合計画から、原材料の在庫状況が算出される、立案開始日における製鉄所（揚地）別、銘柄別の在庫状況を表す情報である。また、立案開始日が立案を実行する日に対して、過去の場合は、各製鉄所がデータベース４００にインプットした原材料の銘柄別の実績在庫状況を表す情報である。

原材料の購入費用情報は、山元（積地）別、銘柄別の原材料の単価（＄／トン（ｔｏｎ、ｔ））を表す情報である。

輸送費用情報は、船舶リストにリストアップされている船舶を利用する場合のフレート、及び、船舶リストにリストアップされている船舶を利用する場合の積揚港別の滞船料を表す情報である。

図１５には、フレートのリストの例を示す。同図に示すように、船舶リストにリストアップされている各船舶について、傭船コード、積港、１揚港、２揚港、３揚港、フレート（＄／トン）が記載されている。例えば連続航海船Ａは、積港Ｘ１から揚港Ａまで航海した場合のフレートが１６．００であり、積港Ｘ１から揚港Ａを経て揚港Ｂまで航海した場合のフレートが１６．２４である。なお、フレートのリストからもわかるように、一般的には、連続航海船を利用した方が不定期船やスポット船を利用するよりもフレートが安い。

図１６には、滞船料のリストの例を示す。同図に示すように、船舶リストにリストアップされている各船舶について、傭船コード、揚ラン（ｔ／Ｄａｙ）、デスデマレート（＄／日）が記載されている。揚ラン（Ｄｉｓｃｈａｒｇｉｎｇ Ｒａｔｅ）とは、契約上の基準となる荷揚能力であり、１日で荷役ができる量を表す。その荷揚能力で荷を揚げると仮定した場合と比較して、実際の揚げ時間が早くなった場合は、デスデマレートで設定された金額を船会社から受け取ることができる。逆に遅くなればデスデマレートに設定された金額を船会社に支払うこととなる。デスデマレート（Ｄｅｓｐａｔｃｈ／Ｄｅｍｕｒａｇｅ Ｒａｔｅ）とは、積・揚港において、早出（Ｄｅｓｐａｔｃｈ）、滞船（Ｄｅｍｕｒａｇｅ）が発生した場合に行われる、課金または払い戻しについて、契約上の料金レートを併せて指す用語である。本明細書の各式中では、デマレート（Ｄｅｍｕｒａｇｅ Ｒａｔｅ）の記載も用いる。例えば、連続航海船Ａが１００００トンの荷揚げをした際に、ＥＴＡから１１時間後にＥＴＤした場合を考える。揚ランは２００００（ｔ／Ｄａｙ）であるため、荷揚には１２時間かかる見込みとなる。この場合に１１時間で荷揚げすれば、デスデマレート１６２５０（＄／日）で規定された金額の１時間分＝１６２５０／２４＄を船会社より受け取る。逆に、ＥＴＡから１３時間後にＥＴＤした場合は、デスデマレート１６２５０（＄／日）で規定された金額の１時間分＝１６２５０／２４＄を船会社に支払う。

（２）船舶財源リストの作成（ステップＳ１０２）
マクロ最適化部２０２の船舶財源リスト作成部２０２ａは、ステップＳ１０１で取り込んだ船舶リスト（図１３を参照）から、配船計画の以降の処理の対象となる、あるいは対象となる可能性のある船舶を選択し、船舶財源リストを作成する。

図１７は、船舶の選択処理を説明するためのフローチャートである。船舶財源リスト作成部２０２ａは、まず、船舶リスト（図１３を参照）及び船舶運航状況（図１４を参照）に基づいて、計画作成期間において運航予定の未定部分がある連続航海船を抽出する（ステップＳ２０１）。例えば立案開始日を２００８年３月１日として３ヶ月分の配船計画を作成するとしたならば、図１４に示すように、連続航海船Ａは２００８年４月１８日以降が未定となっているので、連続航海船Ａは抽出される。

そして、抽出した連続航海船Ａについて計画作成期間における積地と揚地の組み合わせのパターンを全て作成する（ステップＳ２０２）。このとき、積地と揚地との距離等に基づいて特定の条件を設け、この条件を満たすパターンを全て作成しても良い。この場合、例えば明らかに不適な運行距離を持つパターン等を予め排除でき、シミュレーションの効率を上げられる。図１８は、既に確定している航海Ｎｏ．３（図中「Ａ−３」）に続けて、積港Ｘ２から揚港Ａ（航海Ｎｏ．４）、積港Ｘ１から揚港Ｂ（航海Ｎｏ．５）のパターンを作成している様子を示す図である。パターンの作成に際して、各時刻は、標準的な航海時間（港間距離及びこの船舶Ａの標準ノット）や標準的な積揚時間を使用して求めるようにしている。例えば、航海Ｎｏ．４における揚港Ａの沖着時刻は、［航海Ｎｏ．４における積港Ｘ２の沖着時刻］＋［標準積時間］＋［港Ｘ２と港Ａとの距離］／［船舶Ａの標準ノットで求めることができる。もちろん連続航海船Ａについて計画作成期間における積地と揚地の組み合わせは複数あるので、それら全て（あるいは上記の特定条件に合致するパターンを全て）のパターンを作成する。他の連続航海船についても同様の作業を行う。

次に、船舶リスト（図１３を参照）及び船舶運航状況（図１４を参照）に基づいて、計画作成期間において利用可能で未定部分がある不定期船を抽出する（ステップＳ２０３）。例えば、図１３に示すように、不定期船５の配船予定年月は計画作成期間から外れているので、不定期船５は抽出されない。そして、連続航海船の場合と同様に、抽出した各不定期船について計画作成期間における積地と揚地の組み合わせのパターンを全て（あるいは特定条件に合致するパターンを全て）作成する（ステップＳ２０４）。

次に、船舶リスト（図１３を参照）に基づいて、スポット船の候補を抽出する（ステップＳ２０５）。具体的には、まず計画作成期間における総引取目標量を計算する。また、ステップＳ２０１、Ｓ２０２で抽出した連続航海船及び不定期船の最大積載量の合計を計算する。これにより、スポット船で補うべき輸送量を、総引取目標量から、計画作成期間にはいる連続航海船及び不定期船の最大積載量の合計を減算することで、算出することができる（図１９を参照）。このスポット船で補うべき輸送量に基づいて、各スポット船の最大積載量を参照し、何隻のスポット船が必要となるかを計算し、各スポット船の最少船数を求める。例えばスポット船で補うべき輸送量が２５００００トンである場合、豪州-ＰｍａｘＳｐｏｔであれば２５００００÷７５０００＝３．３３で４隻必要となり、４隻の豪州-ＰｍａｘＳｐｏｔを契約するスポット船の候補とする。同様に、２隻の豪州-ＣａｐｅＳｐｏｔ、１隻の豪州-ＶＬＳｐｏｔ、４隻のカナダ-ＰｍａｘＳｐｏｔ、２隻のカナダ-ＣａｐｅＳｐｏｔ、１隻のカナダ-ＶＬＳｐｏｔ、１隻の豪州-ＰｍａｘＳｐｏｔ及び１隻の豪州-ＣａｐｅＳｐｏｔ、・・・のようにスポット船の最少船数が求まる。ここで求めたスポット船の最少船数は、当該傭船コードのスポット船のみで引取を補った場合に必要となる最少のスポット船数となる。後述するように、この最少船数より多くのスポット船が必要になる場合がある。

次に、船舶リスト（図１３を参照）及び船舶運航状況（図１４を参照）に基づいて、スポット船の候補を抽出する。ここでは、船舶運航状況で確定された予定がある場合には、当該船舶をスポット船の候補として抽出し、更に船舶リストの契約区分が未契約の傭船コードのそれぞれに対して、予め設定した日にち毎に、計画作成期間分のスポット船の候補を作成する。図２０に、各傭船コードに対するスポット船の候補を作成する間隔を１０日としたスポット船の航路リストの例を示す。

ここで、上記等間隔で作成した船数と、上記計算した最少船数を比較して、上記等間隔で作成した船数の方が少ない場合には、全てのスポット船の候補を雇ったとしても、引取目標量を満足する引取量を実現することが難しい場合がある。このため、上記計算した最少船数より船数が多くなるように、スポット船の候補を作成する間隔を狭めて、スポット船の候補を作成する。そして、連続航海船の場合と同様に、作成した各スポット船の候補について、計画作成期間における積地と揚地の組み合わせのパターンを全て（あるいは特定条件に合致する組み合わせのパターンを全て）作成する（ステップＳ２０６）。ここで、後述するマクロ最適化において、上記スポット船の各候補について雇う、雇わないが判断され、必要となる船型（船の最大積載量を基に定義される船舶の大きさ）、船数分のスポット船が決定される。例えば、豪州-ＰｍａｘＳｐｏｔ-航海Ｎｏ．３が、候補として作成された後、マクロ最適化にておいて、雇わないと計画されることもある。

（３）マクロ数式モデルの設定（ステップＳ１０３）
マクロ最適化部２０２の数式モデル設定部２０２ｂは、ステップＳ１０２で作成した船舶の運航制約、揚地での原材料の需給バランス制約、引取目標量制約を表すよう構築された数式モデルを設定する。設定を受ける数式モデルは、例えばＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法に則ったモデルとして構築（定式化）されている。ここでは、例としてＭＩＰの定式化に基づいた数式モデルを示す。ここで数式モデルの設定とは、船舶数や港数などの変化に対応できるように抽象的な形式で構築されている基礎数式モデルに対して、各配列の添え字の最大数（例えば船舶数を表す）や、式中の係数の値などを、実際の計画に沿って具体的に定めることを言う。

まず、当該船が、当該積港を選択する、選択しないかを示す変数を定義する。この変数は選択する場合を示す１、選択しない場合を示す０のいずれかの値を取る整数変数とする。後述する最適化によって得られるこの変数の値に基づいて、当該積港を選択するか、選択しないかが判断される。

例えば、図１４に示す連続航海船Ａが候補となる船として挙げられ、この船の航海Ｎｏ．４（図中「Ａ−４」）において、当該船舶の寄航可能な積港がＸ１、Ｘ２の２つある場合には、各積港に対応するように、以下の２つの整数変数を定義する。ここで、これらの整数変数の第３の添え字であるＥＴＡは、ステップＳ１０２で計算された沖着時刻である。

もし、最適化の結果としてＸ１に寄航することが選択された場合は、変数は以下の値を取ることとなる。

また、当該船が、当該積港、当該揚港、当該寄航順（揚港の何番目に寄ったかを表す数字、例えば積港Ｘ１、揚港Ａ、揚港Ｂの順に寄港した場合、揚港Ｂは寄航順２とする）を選択するか、選択しないか、を示す整数変数を定義する。つまり、当該積港に寄った後、当該揚港に、当該寄航順で寄港することを選択する場合、この変数は１の値を取る。一方、このような積港、揚港、寄航順の組み合わせを選択しない場合、この変数は０の値を取る。ここで扱う例では、最大２揚港まで寄航できる例を提示するが、寄航できる揚港数、寄航できる積港数は、それ以上の値を取っても構わない。

また、当該船が、当該積港で、当該銘柄を荷積する量を示す変数を定義する。

また、当該船が、当該積港、当該揚港、当該寄航順で、当該銘柄を荷揚げする量を示す変数を定義する。

更に、当該日、当該銘柄の、当該揚港での在庫量を示す変数を定義する。

「各船舶の積載量が最大積載量を超えないこと」、「積載量は全部荷揚げすること」、等の条件を示す制約式は、基礎となる数式モデルとして予め構築しておく。後述するように、最適化（ステップＳ１０３〜Ｓ１０６）及びシミュレーション（Ｓ１０７）を含む一連の工程は、複数ループ反復して実行できる。初回ループの最適化では、ステップＳ１０１で取り込んだデータに基づいて船舶の運航制約を数式モデルに設定する。第２次ループ（ステップＳ１０３〜Ｓ１０７）以降の最適化ではシミュレータ２０１が前回のループで行ったシミュレーション結果を反映させて、数式モデルを設定する。ここでは、各船舶の積載量が最大積載量を超えないこと、積載量は全部荷揚げすること、等を数式モデルとして設定する。

各船舶の積載量が最大積載量を超えないという制約は、下記の制約式（式２１）と表される。

積載量は全部荷揚げするという制約は、下記の制約式（式２２）と表される。

また、揚地での原材料の需給バランス制約としては、図２１に示すような、「各銘柄の在庫量が常に安全在庫量以上確保されている」という制約条件が、数式モデルとして構築されている。初回ループの最適化では、ステップＳ１０１で取り込んだデータに基づいて、更に次ループ（ステップＳ１０３〜Ｓ１０７）以降はシミュレータ２０１でのシミュレーション結果を反映させて、この数式モデルが設定される。

まず、各銘柄の在庫量の推移を表す制約式は、下記の（式２３）と表される。つまり、当日の在庫量から前日の在庫量と当日に荷揚する量を引いた値は、当日の使用予定量となる。

各銘柄の在庫量が常に安全在庫以上確保されているという制約は、下記の制約式（式２４）と表される。

また、揚港で荷揚げされた原材料は、ヤードに積上げられるが、この荷揚げされる原材料の在庫量は、ヤード能力の上限以下になっていないと着岸できない。但し、ヤードに積上げられた原材料は、日にちが経過すれば、つまり船舶を待たせれば荷揚げすることができる。しかし、荷揚げ量がヤード能力を一定量以上超えていると、滞船時間が無視できないほど大きくなる。そこで、例えば「荷揚げ量について、ヤード能力の１％程度の超過まで許容する」という制約を定式化する。この制約は、下記の制約式（式２５）と表される。

また、引取目標量制約について、初回ループでは、ステップＳ１０１で取り込んだデータに基づいて、更に次ループ（ステップＳ１０３〜Ｓ１０７）以降はシミュレータ２０１でのシミュレーション結果を反映させて、設定が行われる。引取目標量制約に関して、例えば、最適化する引取量（荷積量）が引取目標量から一定幅以上離れないこと、引取の可否（前述したように所定の銘柄については所定の期間は引取しないといった事情もありうる）、等が数式モデルに構築されている。ここで、引取量が引取目標量から一定幅以上離れないという制約を考える場合に、例えば図２２Ａに示すように、単に旬毎（或いは月毎）に引取目標量に対して上下限値を設定し、荷積量がその上下限値を超えないことを制約とすることが考えられる。しかしながら、その場合、例えば荷積量が下限値を満たしているが引取目標量を下回る状況が続いたような場合、年間で蓄積すると、引取割れが発生してしまうこともありうる。そこで、図２２Ｂに示すように、旬毎（或いは月毎）にそれまで引取目標量累積及び引取量累積を考え、引取目標量累積と引取量累積との差を小さくする（最小とする、上下限値を越えないようにする等）制約を設定するのが好適である。上記制約式を定式化するために、旬毎の引取目標累積量からの溢れ量、不足量を示す変数を定義する。

また、各旬の引取量累積を示す変数を定義する。

まず、各銘柄の引取量累積を表す制約式は下記の（式２６）と表される。つまり、引取量累積は、立案開始日から当該旬までの期間にＥＴＡが入っている船舶（航海）の荷揚量の合計となる。

各銘柄の引取目標量累積と溢れ量、不足量との関係を表す制約式は下記の（式２７）と表される。つまり、引取累積量から、溢れが生じている場合は溢れ量を引き、不足が生じている場合は不足量を足すと、引取目標累積量と一致する。ここで、引取累積量と引取目標累積量は近い量を取る程良い計画であるといえる。つまり、この溢れ量、及び不足量は少ない程良い。上記理由のため、後述するようにこの溢れ量、及び不足量は、目的関数の項目として追加され、最適化によってミニマム化される。

ここで、目的関数としてフレートの合計金額のミニマム化を定式化するために、寄港順を示す整数変数を導入する。この寄港順変数は、特定の船舶が、積港として特定の積港、第一揚港として当該揚港１、第二揚港として特定の揚港２の組み合わせを選択する場合を示す１、この組み合わせの寄港順を選択しない場合を示す０、のいずれかの値を取る。

この論理関係を混合整数計画法の定式として記述する方法が、一般的に良く知られており、以下のように定式化することができる。

（４）マクロ数式モデル及び目的関数に基づいて最適化（ステップＳ１０４）
マクロ最適化部２０２の最適化計算部２０２ｃは、ステップＳ１０３で設定した数式モデルを用いて、輸送費用に関して構築された目的関数（評価関数）に基づいて最適化計算を行う。最適化計算に際しては、例えばＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法により最適化問題として問題を解く。

ここでの最適化計算では、輸送費用のうちフレートの合計金額のミニマム化を目的とした目的関数を用い、下記の変数の値が決定される。これにより、フレートの合計金額を最も安価にする船型（船の最大積載量を基に定義される船舶の大きさ）、船数、積揚地（積揚港）、積揚銘柄、積揚量が選定される。

ここで、船舶に掛かるフレートは、積港から第一港目に寄航する揚港までの基準フレートと、上記に更に他の港に寄航した際に掛かる多港揚追加フレートの合計、ここでは、第一揚港から第二揚港に余分に寄航した際に発生する多港揚追加フレートの合計と、積載した量との積となる。

例えば、図１５より連続航海船Ａ−航海Ｎｏ．１がＸ１港から第一揚港Ａに７５０００トンの荷を運んだ際には、基準フレート１６．００となり、フレート（雇船費用）は１６．００×７５０００＝１,２００,０００となる。また、第二揚港としてＡに寄った後、Ｂによると多港揚追加フレートは（１６．２４−１６．００）＝０．２４となり、この際の雇船費用は（１６．００＋０．２４）×７５０００＝１,２１８,０００となる。

以上より、マクロ最適化で用いる目的関数（以下マクロ目的関数と呼ぶ）を式で表すと、下記の式（式３１）を得る。

ここで、マクロ最適化では、引取目標量累積及び引取量累積を考え、引取目標量累積と引取量累積との差を小さくすることも目的としている。このため、旬毎の引取目標累積量からの溢れ量、不足量の合計量をミニマム化する項目を目的関数に追加する。このため、目的関数を表す（式３１）を下記の式（式３２）に変更する。
マクロ最適化では、全体として傭船に関する問題を最適化する。

なお、フレートに関して目的関数を構築することを説明したが、フレートの合計金額及び原材料の購入費用の合計金額のミニマム化を目的とした目的関数としてもよい。既述したように原材料の引取目標量は契約により定められており、原材料の購入費用に大幅な変動はないが、その中でも原材料の購入費用の合計金額のミニマム化が可能になる。

上記の項目（３）、（４）で説明した如く、ミニマム化すべき式が目的関数、満足すべき各式が制約式として定式化されている。上記の制約式は線形等式、或いは不等式で表現されており、上記の目的関数は１次式で表されている。また、変数の中に整数となるべき変数が存在するモデルとして数式モデル、目的関数が構築されている。このように定式化された問題は、混合整数計画問題として一般に良く知られており、本問題は（解析的に）最適化することが可能である。

マクロ最適化に際しては、時間精度を旬精度として演算する。最適化期間は、最初のループ（ステップＳ１０３〜Ｓ１０７）では９旬、次ループ（ステップＳ１０３〜Ｓ１０７）では８旬、・・・、最後のループ（ステップＳ１０３〜Ｓ１０７）では１旬とする。時間精度は、旬精度として演算が行われる。そして、最適化期間（９旬〜１旬）のうちの最初の１旬を計画確定期間とし、その計画確定期間での演算結果をミクロ最適化部２０３に出力する。

（５）ミクロ数式モデルの設定（ステップＳ１０５）
ミクロ最適化部２０３の数式モデル設定部２０３ａは、マクロ最適化部２０２で求めた計画確定期間の配船計画における船舶の運航制約のうち、滞船制約、及び、揚地での原材料の需給バランス制約を表す数式モデルを設定する。ここで用いられる数式モデルは、例えばＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法に則って構築されている。ここでは、例としてＭＩＰの定式化に基づいた数式モデルを示す。

マクロ最適化により、寄航する揚港が決定されている。ここで、揚港には船が着岸するために複数のバース（岸壁）が存在するため、当該揚港のいずれのバースに着岸するかを選択する変数を定義する。この変数は当該バースを選択する場合を示す１、選択しない場合を示す０のいずれかの値を取る整数変数とする。

また、当該船が当該バースに着岸するために沖待を開始する時刻（ＥＴＡ）を示す変数を定義する。時刻はＭＩＰで定式化するための変数として直接定義できないため、立案開始日からの経過分として定義する。つまり、立案開始日が１月１日０時０分の場合で、ＥＴＡが１月１日１時１０分の場合は、７０という値を取るとして定義する。また、この変数は整数変数ではなく、連続値を取る変数として定義する。

同様にＥＴＢを示す変数を定義する。

同様にＥＴＤを示す変数を定義する。

更に、当該分、当該銘柄の、当該揚港での在庫量を示す変数を定義する。

船舶の滞船制約についても、初回ループではステップＳ１０１で取り込んだデータに基づいて、更に次ループ（ステップＳ１０３〜Ｓ１０７）以降はシミュレータ２０１でのシミュレーション結果を反映させて設定する。船舶の滞船制約に関しては、揚港での船舶運行条件（ＥＴＢ＞ＥＴＡ、ＥＴＤ＞ＥＴＢ＋荷揚時間等）、バースの条件（許容されるＬＯＡ（全長）、ＤＲＡＦＴ（全深）、ＢＥＡＭ（全幅）、積揚能力、ヤード能力等）、等の条件が数式モデルとして構築されている。ここでは、これらの数式モデルが設定される。

マクロ最適化で揚港が決定された船舶は、当該揚港の何れかのバースに着岸する必要がある。このため、この制約は、下記の制約式（式３３）と表される。つまり当該船に対して、着岸可能なバースの内で、必ず一つのバースが選択（変数の値が１）される必要がある。

ＥＴＢはＥＴＡ以降となる必要がある。この制約は、下記の制約式（式３４）と表される。

ＥＴＤはＥＴＢ＋荷揚時間以降となる必要がある。ここで、マクロ最適化により当該バースでの荷揚量は決定されているため、当該バースでの荷揚時間は、当該バースでの標準的な荷揚能力を用いると、荷揚時間＝荷揚量／荷揚能力となる。以上より、上記制約は、下記の制約式（式３５）と表される。

また、上記マクロ最適化と類似の、揚地での原材料の需給バランス制約が、ミクロ最適化（ステップＳ１０５〜Ｓ１０７）でも用いられる。この需給バランス制約条件も、初回ループのミクロ最適化ではステップＳ１０１で取り込んだデータに基づいて、次ループ以降のミクロ最適化ではシミュレータ２０１でのシミュレーション結果を反映させて、設定される。ミクロ最適化でも、図２１に示すような、「各銘柄の在庫量が常に安全在庫量以上確保されている」という制約条件を表すように構築された数式モデルが設定される。ただし、ミクロ最適化では、用いられる時間精度がマクロ最適化と異なる。このため、ここで用いられる制約条件は、「当該時刻の在庫量から１分前の在庫量と当該時刻に荷揚する量を引いた値が、当該時刻１分間の使用予定量である」となる。

また、揚港で荷揚げされた原材料は、ヤードに積上げられるが、この荷揚げされた原材料の在庫量は、ヤード能力の上限以下になっていないと着岸できない。つまりＥＴＢ時点での在庫量はヤード能力上限以下となる必要がある。この制約は、下記の制約式（式３７）と表される。

（６）ミクロ数式モデル及び目的関数に基づいて最適化（ステップＳ１０６）
ミクロ最適化部２０３の最適化計算部２０３ｂは、ステップＳ１０５で設定した数式モデルを用いて、輸送費用に関して構築された目的関数（評価関数）に基づいて最適化計算を行う。最適化計算に際しては、例えばＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法により最適化問題として問題を解く。

ここでの最適化計算では、滞船料の合計金額のミニマム化を目的とした目的関数を用い、船舶が着岸する／しないを表すδ（船舶、バース）、ＥＴＡ時刻を表すＥＴＡ（船舶、バース）、ＥＴＢ時刻を表すＥＴＢ（船舶、バース）、ＥＴＤ時刻を表すＥＴＤ（船舶、バース）といった変数を決定する。これにより、輸送費用を最も安価にするバース、入出港タイミングが選定される。

ここで、船舶に掛かる滞船料は、ＥＴＤ−ＥＴＡと契約上の基準停泊時間とを比較し、基準停泊時間より、停泊が長い、つまり、ＥＴＤ−ＥＴＡ＞基準停泊時間の場合には、デスデマレートとして契約された費用を支払い、逆の場合には、デスデマレートとして契約された費用を受取ることとなる。基準停泊時間は、契約上設定された荷揚能力である揚ランを用いて揚量／揚ランで計算される。例えば、連続航海船Ａ−航海Ｎｏ．１が揚港で１００００トンの荷揚げを行い、ＥＴＡからＥＴＤまで１１時間掛かった場合は、基準停泊時間＝１００００／２００００＝０．５日、１２時時間より１時間早いため、デスデマレートで設定された１６２５０／２４の金額を船会社より受取ることとなる。以上より、ミクロ最適化で用いる目的関数（以下ミクロ目的関数と呼ぶ）を式で表すと、下記の式（式３８）を得る。

上記式では定数部分が含まれているが、ミニマム化において定数部分は影響を与えないため、下記の式（式３９）が目的関数となる。

上記の項目（５）、（６）で説明した如く、ミニマム化すべき式が目的関数、満足すべき各式が制約式として定式化されている。上記の制約式は、線形等式、或いは不等式で表現され、上記の目的関数は１次式で表されている。また、変数の中に整数となるべき変数が存在するモデルとして数式モデル、目的関数が構築されている。このように定式化された問題は、混合整数計画問題として一般に良く知られており、本問題は（解析的に）最適化することが可能である。

ミクロ最適化に際しては、最適化期間を１０日（１旬）とし、時間精度を分精度として演算する。

（７）シミュレーション（ステップＳ１０７）
シミュレータ２０１は、ミクロ最適化部２０３で求めた数式モデルに対する解に基づいてシミュレーションを実行して、計画確定期間（１旬）の配船計画を確定する。シミュレーションの時間精度は分精度とする。このシミュレーションでは、マクロ数式モデル、ミクロ数式モデルには組み込むことができなかった制約等も組み込むことで、実際に求められる細かな制約までも考慮した配船計画を作成することが可能となる。

例えば、マクロ・ミクロ最適化で取扱うことが難しい制約の一例として１隻の船舶の荷揚げに使用するアンローダの基数がある。この基数は、荷揚げする銘柄が積載されているハッチの位置、同一揚港の別バースで荷役している船舶があるか、ないか等により変わって来る。この荷揚げに使用するアンローダの基数により荷揚能力は変わって来る。例えば、「１基で荷揚げする場合は、１５００ｔ／ｈで１００％能力で揚げられる」、または、「２基の場合は、１５００ｔ／ｈ×２基で７０％能力で揚げられる」等の状況が例示できる。マクロ・ミクロ最適化ではこれらのアンローダ基数まで考慮されていないため、最適化で計算された時間を、シミュレータによりアンローダ基数まで考慮して、最適化の時間のずれ等をシミュレーションに取込み、正確にシミュレートすることで、実操業に求められる細かな制約まで考慮した生産・物流計画の立案が可能となる。

シミュレータ２０１では、ミクロ最適化部２０３で船舶の入出港タイミングの入れ替え等による時間調整があった場合、それを波及的に反映させて時刻修正する。特に連続航海船では、ある港で時間調整があった場合その後の航海にも波及的に影響するので、シミュレータ２０１で時刻修正を行い、その後のマクロ最適化部２０２での処理に反映させるようにしている。

（８）立案開始日の更新（ステップＳ１０９）
ステップＳ１０８において計画作成期間（３ヶ月（９旬））分の計画が確定したかどうかを判定する。まだ確定していない場合、計画が確定した旬の次旬の初日、例えばＮ旬の計画が確定したならばＮ＋１旬の初日を立案更新日として更新し（ステップＳ１０９）、ステップＳ１０３に戻る。ステップＳ１０３から始まる次ループでは、計画が確定した旬（Ｎ旬）における在庫推移や船舶の運航状況を更新して、次旬（Ｎ＋１旬）の計画を確定させる。これを繰り返すことにより、計画作成期間（３ヶ月）分の計画が確定することになる（図２３を参照）。

（９）配船計画の出力（ステップＳ１１０）
以上のようにして作成した配船計画は、出力部２０５により、不図示のモニタに画面表示されたり、データベース４００を含む外部機器にデータ送信されたりする。

以上述べたように、マクロ最適化部２０２及びミクロ最適化部２０３では、まず入力データもしくは前ループから得られる、初期値（初期条件）に基づいて数式モデルを設定し、最適化計算を行い、シミュレータ２０１に対する指示を算出する。シミュレータ２０１は、計画確定期間（１旬）についてシミュレーションを終了すると、計画確定期間の最終状態での原材料の在庫推移、船舶の運航状況の推移の情報をマクロ最適化部２０２及びミクロ最適化部２０３に与える。マクロ最適化部２０２及びミクロ最適化部２０３は、その与えられた情報に基づいて数式モデルを設定し、最適化計算を行い、シミュレータに対する指示を算出する。このようにシミュレータ２０１と最適化部２０２、２０３を連動させることにより、計画作成期間（３ヶ月（９旬））の配船計画を作成することができる。

本実施形態に係る配船計画作成装置（方法）によれば、マクロ最適化部２０２及びミクロ最適化部２０３により行われた最適化計算の結果に基づいた計算指示をシミュレータ２０１（在庫推移シミュレータ、船舶運航状況推移シミュレータ）に出力する。このように、最適化計算の結果に基づいてシミュレーションが行われるものであるので、理論的な最適解を確実に得ることが可能となる。これにより、従来のようにシミュレーション結果を評価してシミュレーションを何回も繰り返して実行する必要がなく、シミュレーション結果を迅速かつ高精度に作成することができる。

また、シミュレータ２０１の規模が非常に大きい場合や制約条件が非常に多くて複雑な場合でも、シミュレータ２０１に記載された制約のうち、配船計画の作成に影響が大きい重要な部分のみを数式モデルに取り込むようにすることで、数式モデル設定部２０２ｂ、１０３ａの規模を適切な範囲にして、実用的な時間内で最適化計算を行うようにすることができる。シミュレータ２０１には、考慮すべき制約を全て記載することができるので、１回のシミュレーションを実行して作成された配船計画は現実に実行可能となることが保証される。

また、配船計画を作成する場合、ブラジル等の遠方より輸送される銘柄は、２ヶ月或いは3ヶ月に一度しか入荷されないといったこともあるため、長期間を考慮して配船計画を立てる必要がある。一方で、中国等頻繁に輸送される銘柄では数日で搬送される銘柄も存在する。更にバースの管理は、滞船料が発生することもあり、分単位で行われるため、分精度の計画が要求される。これらの要求に対して、マクロ最適化部２０２で船型（船の最大積載量を基に定義される船舶の大きさ）、船数、積揚地（積揚港）、積揚銘柄、積揚量、寄港順を選定し、ミクロ最適化部２０３で使用バース、入出港タイミングを選定するように演算の分担を行った。このため、負荷を抑えるとともに、高精度で求解可能となる。すなわち、マクロ最適化とミクロ最適化を連動させ、繰り返し実行することで、長期間（例えば３ヶ月）で特定する必要がある使用可能な船、積揚地、銘柄、量を長期間で考慮すると同時に、細かな時間精度が要求される使用バース、入出港タイミングは、短期間の（細かい）時間精度で最適化することを可能とした。

また、配船計画作成の際に、船型（船の最大積載量を基に定義される船舶の大きさ）、船数、積揚地（積揚港）、積揚銘柄、積揚量をユーザが個別に固定できるようにしてもよい。例えば所定の船舶を使用する、所定の積港を利用する等が予め決まっているような事情もあるからである。特に、引取目標量として設定した量に基づいて配船計画を立案した後で、山元との交渉が進み、引取量が確定されるが、この際には、引取量、積地（積港）、積銘柄、積量（荷積量）は、契約の都合上変更することが許されない。しかし、揚地に関しては、在庫状況を判断して、揚地、揚銘柄、揚量を変更する余地が残される場合が多い。このため、積地（積港）、積銘柄、積量を一括して固定化できるような操作を可能にすれば、ユーザにとって利便性が高くなる。

なお、ここまで述べた配船計画作成手法では、マクロ最適化部２０２の最適化計算部２０２ｃで輸送費用に関して構築された目的関数（評価関数）に基づいて最適化計算を行う例を説明したが、他の目的関数を加えてもよい。

例えば積地における負荷を平準化するために、同一の積地に入出港する船舶が同時期に集中したり、逆に船舶が入出港しない期間が続いたりすることを避ける、すなわち同一の積地では計画作成期間中にできるだけ均等に配船することが求められる。

そこで、図２４に示すように、積地毎に全銘柄の引取量を旬単位（或いは月単位）に集計し、それまでの累積を考える（引取量累積）。また、積地毎に全銘柄の引取目標量を旬単位（或いは月単位）に集計し、それまでの累積を目標値として設定する（引取目標量累積）。そして、引取量累積と引取目標量累積の差のミニマム化を目的とした目的関数を構築する。

これにより、各旬（或いは各月）で積地の引取量が均等に近づき、換言すれば、均等配船が可能になる。

同様に、揚地における負荷を平準化するために、同一の揚地に入出港する船舶が同時期に集中したり、逆に船舶が入出港しない期間が続いたりすることを避ける、すなわち同一の揚地では計画作成期間中にできるだけ均等に配船することが求められる。

そこで、図２５に示すように、揚地毎に全銘柄の荷揚量を旬単位（或いは月単位）に集計し、それまでの累積を考える（荷揚量累積）。また、揚地毎に標準荷揚能力量を旬単位（或いは月単位）に集計し、それまでの累積を目標値として設定する（揚地標準荷揚能力量累積）。そして、その差を残荷揚量と定義し、この残荷揚量のミニマム化を目的とした目的関数を構築する。

これにより、各旬（或いは各月）で揚地の荷揚量が均等に近づき、換言すれば、均等配船が可能になる。即ち、マクロ最適化においても、滞船を抑制することが可能となる。

以上により、複数銘柄の原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を相互に連係させて一括で作成することが可能となる。また、これら計画においては、配合計画・配船計画を通して全体として輸送費用をミニマム化することが可能となる。

（第２の実施形態）
第２の実施形態として、第２の配合計画作成装置３００を導入する例を説明する。第２の配合計画作成装置３００は、配船計画作成装置２００により作成された配船計画に基づいて、複数銘柄の配合原材料を入荷して混合する配合計画を作成する。なお、第２の配合計画作成装置３００の構成や基本的な処理動作は第１の配合計画作成装置１００と同様であり、ここでも図２〜１１を参照して説明する。

図２は、第２の配合計画作成装置３００を含むシステム構成例を示す図である。図２に示すように、第２の配合計画作成装置３００は、配合計画を作成するに際して、配合計画を立案する上で必要となる以下のような制約条件、前提条件のデータを操業者が設定するか、或いはデータベース４００から取り込む：計画作成期間、原材料の入荷予定（配船計画作成装置２００により作成された配船計画による入荷量）、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、船舶を利用する際の輸送費用情報を含む。

第２の配合計画作成装置３００は、多種類（複数銘柄）の原材料を入荷して混合する混合計画を、シミュレーションを実行して作成する。この混合計画は、原材料の需給バランス制約、混合後の性状制約を満たすような、各銘柄の使用量（配合割合）を含む。第２の配合計画作成装置３００では、後述するように、ＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法に則って構築された、原材料の需給バランス制約を表す数式モデル（「需給バランスモデル」とも称する）、及び、混合後の性状制約を表す数式モデル（「性状モデル」とも称する）を、各条件に沿って設定し、求解することにより配合計画の最適化を図る。

ここで、第１の配合計画作成装置１００では、入荷量は、入荷量の累積量が引取目標累積量と、大幅に外れないように、配合計画を作成した。しかし、第２の配合計画作成装置３００において、入荷量は、配船計画作成装置２００により船舶毎の荷揚量として予め決定されているため、既知の量となる。このため、第２の配合計画作成装置３００では、入荷量は求める量ではなく、固定されている量となる。また、上記理由により引取目標量を考慮する必要はない。

また、第１の実施形態においては、第１の配合計画作成装置１００により作成された使用計画を、配船計画作成装置２００の入力データとして、配船計画を作成している。この場合、第１の配合計画作成装置１００において、引取目標累積量に基づいて計画を作成している。このとき、例えば引取目標累積量が原材料Ｘに対して、５万トンと設定されており、使用計画、入荷計画共に引取目標累積量と同じ５万トンが累積量として、計画された場合を考える。上記使用計画を配船計画作成装置の入力データとして使用した場合、船舶リストに登録された船の最大積載量が４．５万トンの船１隻しか存在しない場合に、入荷量は４．５万トンとなり、下記の如く、在庫切れを発生する可能性がある。
在庫量＝４．５万トン − ５万トン ＝ −０．５万トン

上記の例は極端な例ではあるが、第１の配合計画作成装置１００では、引取目標量を基に配合計画を作成し、作成された使用計画を基に配船計画を作成しているため、引取目標量と船舶の最大積載量がマッチしていない場合、上記の例のような現象は起こりえる。

そこで、より良い配合計画を作成するために、引取目標量ではなく、配船計画作成装置１００により決定された配船計画を入力データとして、第２の配合計画作成装置３００により、配合計画を修正する。これにより、上記のような微妙な引取目標量と船舶の最大積載量がマッチしていないために発生する可能性のある在庫切れを防ぐことが可能となる。
本実施形態では、船舶の運賃は考慮するが、バースでの着岸時刻など詳細な運航状況まで考慮することなく、引取目標量を元に第１の配合計画作成装置により計画された使用量を、配船計画側の在庫推移シミュレータ、及び船舶運航状況推移シミュレータによるシミュレーションの結果に現れる、詳細な配船事情まで考慮して、配合計画を再度更新できる。この作用により、配合計画の精度の向上に大きな効果が得られる。

表示部３０３では、第２の配合計画作成装置３００で求められた各銘柄の使用量（配合割合）、入荷量、在庫推移グラフ、各種帳票を表示する。

操業者評価部３０４では、求められた配合計画を様々な観点（例えば、在庫推移、性状等）から操業者が評価し、満足のいく結果でなければ、必要に応じて配合割合等を修正する。その際に、必要に応じて目的関数の重みや評価の指標を変えたり、数式モデルを設定する対象期間・計画確定期間を変えたりする。また、全部の或いは指定した処理のみ使用量の固定をする等、操業者の意志を反映させた数式モデルの設定が可能である。そして、第２の配合計画作成装置３００で再度配合計画を作成し直す。

図３は、第２の配合計画作成装置の基本的な構成を示すブロック図である。図３に示すように、第２の配合計画作成装置３００は、シミュレータ（在庫推移シミュレータ３１１、性状シミュレータ３１２を含む）、本発明でいう数式モデル設定部として機能するモデル設定部（需給バランスモデル設定部３１３、性状モデル設定部３１４）、本発明でいう最適化計算部として機能する計画部３１５を含んで構成され、更に入出力部を併せ持つ。

在庫推移シミュレータ３１１は、各原材料の需給状態（在庫推移）を計算するシミュレータである。性状シミュレータ３１２は、原材料を混合した後の性状を計算するシミュレータである。在庫推移シミュレータ３１１、性状シミュレータ３１２が互いに連動することで、原材料の在庫推移、混合後の性状を計算する。

本実施形態においては、例えば以下の情報を含む、入力データ３１６に基づいて数式モデルの設定処理を行う：計画作成期間、原材料の入荷予定、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、及び船舶を利用する際の輸送費用情報。配合計画の立案開始日時から予め設定された最適化期間分を対象として、予め設定した時間精度に基づいて、ＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理整数計画法等に則り、需給バランスモデル設定部３１３にて需給バランス制約（在庫制約）を表す数式モデルが設定され、性状モデル設定部３１４により性状制約を表す数式モデルが設定される。

需給バランスモデル設定部３１３、性状モデル設定部３１４により設定された数式モデルを用いて、在庫を切らさないようにするとともに要求される性状を満足し、かつ、費用（原材料の購入費用及び輸送費用）をミニマム化して配合計画を作成するように、計画部３１５により最適化計算を行う。この最適化計算の結果に基づいて、在庫推移シミュレータ３１１、性状シミュレータ３１２に対する計算指示を算出する。この計算指示を受けて、在庫推移シミュレータ３１１が在庫推移をシミュレートし、性状シミュレータ３１２が計画に従って製造される製品・半製品の性状をシミュレートする。例えば、鉄鋼業においては、石炭を混合して焼き固めたコークス（製品）、鉄鉱石を還元して得られる銑鉄を精錬した溶鋼を凝固させたスラブ（半製品）等の性状をシミュレートする。

かかる第２の配合計画作成装置３００によれば、従来のように予め決められたルールに基づいて計算指示が行われるのではなく、計画部３１５により行われた最適化計算の結果に基づいた計算指示を在庫推移シミュレータ３１１、性状シミュレータ３１２に出力する。このため、そのときの事象に応じた最適な計算指示を確実に行うことが可能となる。

また、例えば、図７に示すように予め設定された計画確定期間分について、在庫推移シミュレータ３１１、性状シミュレータ３１２によるシミュレーションが、行われる。またこのシミュレーションが終了すると、立案開始日が更新され、更新前の計画確定期間の最終状態、つまり更新後の立案開始日での在庫推移、性状の情報に基づいて、新たな最適化期間分の需給バランスモデル設定部３１３により在庫制約を表す数式モデルが設定され、性状モデル設定部３１４により性状制約を表す数式モデルが設定され、計画部３１５に与えられる。この在庫推移、性状の情報が与えられると、計画部３１５は最適化計算を実行する。

第２の配合計画作成装置３００においても、第１の実施形態で述べたように、シミュレータ（在庫推移シミュレータ３１１、性状シミュレータ３１２）とモデル設定部（需給バランスモデル設定部３１３、性状モデル設定部３１４）と計画部３１５とを連動させて配合計画を作成するようにした。このため以下の作用効果が得られる。（１）シミュレーションを繰り返して実行せずに配合計画を作成することができる。（２）配合計画作成に影響が大きい重要な部分のみを計画部３１５に取り込むようにすることで計算時間を短縮することができるとともに、（３）大規模問題を解くことが可能になる。

以下、第２の配合計画作成装置３００の構成及びこの装置を用いて実行する配合計画作成方法の各ステップをより詳細に説明する。図４は、図３を用いて説明した第２の配合計画作成装置３００の基本的な構成に対する、配合計画作成装置の詳細な構成を示す図である。また、図５は、この装置を用いて実行する配合計画作成方法の各ステップを示すフローチャートである。

配合計画作成の概要として、例えば図６に示すように、、以下のような計算、調整の工程が含まれる：複数ある製鉄所（揚港）ａ〜ｃでの原材料（銘柄）の需給バランスを取ること；（各銘柄Ａ〜Ｎの在庫を切らさない等）要求される性状を満足させること；かつ、費用（原材料の購入費用及び輸送費用）をミニマム化すること；前記の条件を満たすような配合計画として製鉄所ａ〜ｃ毎の各銘柄Ａ〜Ｎの使用量（配合割合）、入荷量を決定すること。ここで、製鉄所毎に使用量の合計量である予定使用量は、入力データとして与えられる。このため、配合割合（％）＝使用量／予定使用量×１００となる。このため、使用量、配合割合の一方が決定されれば、他方が決定されることとなる。

（１）入力データの取込み（図４の入力データ取込み部３５１、図５のステップＳ３０１）
本処理に必要な情報（原材料の入荷予定、原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報、船舶を利用する際の輸送費用情報等）をオンラインにて読み込み、必要に応じて操業者が修正を加える。

ここで、入力データ取込み部３５１により取り込まれる原材料の入荷予定には、配船計画（船舶１隻毎の積港、積港着の日時、積銘柄、積量、揚港、揚港着の日時、揚銘柄、揚量を含む項目について計画）といった、入荷量を表す情報が含まれる。例えば、図２６に示す配船計画では、図１３に示したような船舶リストにリストアップされている各船舶の運航予定が組まれている。配船計画では、船舶リストにリストアップされた船舶について、積港、積港着の日時、銘柄、積量、揚港、揚港着の日時、銘柄、揚量を含む項目について計画が立案されている。

例えば、図２６に示す配船計画では、連続航海船Ａの航海Ｎｏ．１は、２００７年１１月１９日２１時に積港（Ｘ１港）沖に着き、２００７年１２月１３日２１時に積港（Ｘ１港）のコード「１」で表されるバースに着岸し、２００７年１２月１４日２１時に積港（Ｘ１港）を出港する。この際に原材料の銘柄Ａを４００００トン、銘柄Ｂを３５０００トン積載する。その後、１６９８０分航海して、２００７年１２月２６日１６時に揚港（Ａ港）沖に着き、２００７年１２月２７日１時に揚港（Ａ港）のコード「４」で表されるバースに着岸し、２００７年１２月２８日１６時に揚港（Ａ港）を出港する。この際に原材料の銘柄Ａを２５０００トン、銘柄Ｂを１５０００トン荷揚げする。その後、３０６０分航海して、２００７年１２月３０日１９時に揚港（Ｂ港）沖に着き、２００７年１２月３０日１９時に揚港（Ｂ港）のコード「１３」で表されるバースに着岸し、２００８年１月１日２３時に揚港（Ｂ港）を出港する。この際に原材料の銘柄Ａを１５０００トン、銘柄Ｂを２００００トン荷揚げする。上記例示した連続航海船Ａの航海Ｎｏ．１においては、揚港（Ａ港）に２００７年１２月２７日、原材料の銘柄Ａを２５０００トン、銘柄Ｂを１５０００トン、揚港（Ｂ港）に２００７年１２月３０日、原材料の銘柄Ａを１５０００トン、銘柄Ｂを２００００トン、入荷したとする、配合計画作成装置の入力データとなる。

原材料の在庫状況、原材料の性状、原材料の単価を表す購入費用情報については、第１の実施形態で説明した第１の配合計画作成装置１００の場合と同様である。

また、船舶を利用する際の輸送費用情報については、配船計画作成装置２００により、将来に亘り計画された配船計画が、第２の配合計画作成装置３００に入力データとして渡される。このため、銘柄別・揚港別見做しフレートを使用する必要はなく、船舶別・積港別・揚港別フレートを使用することができる。

以上説明した入力データ取込み部３５１及びステップＳ３０１が、本発明でいう第２の配合計画作成部のデータ取込み部及びそれによる処理の例である。

（２）配合計画作成期間の設定（図４の計画作成期間設定部３５２、図５のステップＳ３０２）
配合計画を作成する期間を設定する。この作成期間は立案者の必要に応じて任意の期間を設定可能とする。ここでは、一例として１０日間分を立案する。

（３）配合計画作成時間精度の設定（図４の時間精度設定部３５３、図５のステップＳ３０３）
配合計画を作成する時間精度並びにシミュレーション精度を設定する。この時間精度並びにシミュレーション精度は、立案者の必要に応じて個別に任意の精度を設定可能とする。例えば立案の細かな精度を必要とする計画作成期間の前半では精度を細かくし、粗い計画で十分な計画作成期間の後半では精度を粗くすることで、十分な精度と短い計算時間での効率的な計画作成が可能になる。

（４）最適化期間の設定（図４の最適化期間設定部３５４、図５のステップＳ３０４）
配合計画を作成する最適化期間を設定する。この最適化期間は立案者の必要に応じて個別に任意の対象期間を設定可能とする。ここでは、一例として計画作成期間を通して、最適化期間を３日間とする。

（５）計画確定期間の設定（図４の計画確定期間設定部３５５、図５のステップＳ３０５）
配合計画を確定する計画確定期間を設定する。この計画確定期間は、立案者の必要に応じて個別に任意の期間を設定可能とする。例えば、立案の細かな精度を必要とする計画作成期間の前半では計画確定期間を短くし、粗い計画で十分な計画作成期間の後半では計画確定期間を長くする。このことで、十分な精度を持ちながら短い計算時間で、効率的な計画作成が可能になる。ここでは、一例として、計画確定期間を１日に設定する。この場合は、数式モデルに対する解に基づいてシミュレートした結果得られる配合計画に対しては、計画作成期間を通して最初の１日分を確定する。

（６）配合計画の需給バランス制約を数式モデルに設定（図３の需給バランスモデル設定部３１３、図４の需給バランスモデル設定部３５６、図５のステップＳ３０６）
入力データ取込み部３５１により取込まれたデータの全部又は一部に基づいて、設定された最適化期間分を設定された時間精度で需給バランス制約に基づいて数式モデルを設定する。

第１の配合計画作成装置１００の場合と同様に、各銘柄の在庫量が一定の安全在庫量と呼ばれる値以上であることが要求される。この場合の制約は、上記の（式４）と表される。

また、各銘柄の在庫量は、前日の在庫量、前日の入荷量、前日の使用量より決定される。この場合の関係式を表す制約は、上記の（式５）と表される。つまり、当日の在庫量は、前日の在庫量と当日に入荷（荷揚）する量を足した値から、当日の使用量を引いた値となる。

また、各銘柄の使用量のある日の合計は、当該日の全銘柄合計に対して予定された使用量と一致する必要がある。この場合の関係を表す制約式は、上記の（式６）と表される。

また、各種原材料の購買に対する要因等に鑑み、操業者は、目標とする配合割合を設定し、この与えた目標とする配合割合に近い配合割合を実現する配合計画が作成されることを求める。つまり配合割合が操業者の想定と大きくかけ離れると、想定した購買量を満たせなくなったり、購買量を越えたり、また操業設備に無理な操業を及ぼすことが想定される。このため、目標として与えた配合割合に近い配合割合が出力されることが必要となる。上記機能を実現するための制約を以下に示す。つまり、銘柄の使用量から使用目標量（目標とする配合割合）（定数）を引いた値を、使用目標量からの溢れ量の変数として定義する。ここで、使用量と使用目標量が近い量を取る程、良い計画であるため、この溢れ量は少ない程良い。上記理由のため、後述するように、この溢れ量が、目的関数の項目として追加され、最適化によってミニマム化される。同様に銘柄の使用目標量から使用量を引いた値を、使用目標量からの不足量の変数として定義する。ここで、使用量と使用目標量は近い量を取る程良い計画であるため、この不足量は少ない程良い。上記理由のため、後述するようにこの不足は、目的関数の項目として追加され、最適化によってミニマム化される。この場合、各銘柄の使用量、使用目標量、溢れ量、不足量との関係を表す制約式は上記の（式７）と表される。つまり、溢れが生じる場合は使用量から溢れ量を引き、不足が生じる場合は不足量を足すと、使用目標量と一致する。

更に、前日の配合割合とその翌日の配合割合が大きく乖離すると、操業に困難を来たす。つまり、別原材料を使用するための段取り時間の増加や、設備の故障の原因となる。このため、前日の配合割合とその翌日の配合割合が大きく乖離することがない配合計画が、作成されることを求める。これを実現するための制約を上記の（式９）に示す。

なお、上述した需給バランス制約は一例であり、他の制約に替えたり、他の制約を加えたりしてもよい。

（７）配合計画の性状制約を用いて数式モデルを設定（図３の性状モデル設定部３１４、図４の線形化部３５７ａを含む性状モデル設定部３５７、図５のステップＳ３０７、Ｓ３０７ａ）
入力データ取込み部３５１により取込まれたデータの全部又は一部に基づいて、設定された最適化期間及び時間精度を用いて、性状制約を用いて数式モデルを設定する。鉄鉱石の配合計画を作成する場合、性状としてはに用いられる原材料の性状としては、例えば、以下のものが挙げられる：鉄分、SiO₂、Al₂O₃、SiO₂、等。石炭の配合計画を作成する場合の性状としては、、以下のものが挙げられる：ＣＳＲ（熱間反応後強度）、ＤＩ（コークス強度）、ＶＭ（揮発分）、膨張圧等。これらの性状が、要求される性状制約を満たす必要がある。混合後の性状モデルの一例を上記の（式１４）に示した。なお、（式１４）では下限値Ｓを有する例を示すが、上限値を有する場合、上限値及び下限値の両方を有する場合もありうる。

ここで、多くの性状について、性状モデルに含まれる数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）は、上記の（式１５）に示すように、配合割合に対して線形となる。例えば、SiO₂に関して、銘柄Ａの配合割合が４０％、銘柄ＡのSiO₂成分が１%、銘柄Ｂの配合割合が６０％、SiO₂成分が２%の条件で混合した場合、混合後のSiO₂成分の性状は、１×０．４＋２×０．６＝１．６％となる。

ところが、性状によっては、その性状を表す数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）が非線形となることがある。この場合、次に述べるように、線形化部３５７ａで、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に代えて線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）を導入して数式モデルを定式化する。

線形化部３５７ａでの処理について説明する。ある性状を表す数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）が非線形である場合、それに代えて線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）を導入する。この線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）は、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）の下限をなすもの、すなわち上記の（式１６）の関係が成立するものを考える。なお、（式１６）は常に成立する必要はなく、必要な範囲で成立していればよい。

例えば線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）として、上記の（式１７)に示す加重平均を考える。加重平均は、単一銘柄を１００％使用した場合の性状を非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）から求め、配合割合を乗算して、使用銘柄分足し合わせた値である。

説明を簡単にするため、銘柄Ａの配合割合が９０％、銘柄Ｃの配合割合が１０％の例を考える。この場合、線形の数式ｆ´（９０、０、１０、・・・、０）となる加重平均は、下式で表される。
ｆ´（９０、０、１０、・・・、０）
＝０．９×ｆ（１００、０、・・・、０）＋０．１×ｆ（０、０、１００、・・・０）

過去の実績等から、この加重平均が（式１６）を満たすことがわかっていれば、この加重平均を線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）として利用することができる。すなわち、加重平均≧Ｓを制約とすれば、（式１４）が成立するものとみなすことで、定式化できる可能性が得られる。

線形化部３５７ａでは、上記の（式１４）´に示すように、線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）に対する下限値として、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に対する下限値Ｓよりも小さな仮下限値Ｓ´＝Ｓ−ｓ（ｓ：オフセット値）を設定して数式モデルを設定する。

以上は、混合後の性状制約が下限値を有する場合を例に説明した。なお、上述した性状制約は一例であり、他の制約に替えたり（混合後の性状制約が上限値を有する場合を含む）、他の制約を加えたりしてもよい。

以上説明した需給バランスモデル設定部３５６（需給バランスモデル設定部３１３）及びステップＳ３０６、並びに、性状モデル設定部３５７（図３の性状モデル設定部３１４）及びステップＳ３０７、Ｓ３０７ａが、本発明でいう第２の配合計画作成部の数式モデル設定部及びそれによる処理例である。

（８）固定化抽出処理（図４の固定化抽出処理部３５８、図５のステップＳ３０８）
図９に示すように、配船計画の項目である積港、積銘柄、積量、揚港、揚銘柄、揚量のうち固定化されているもの、すなわち変更できないものを抽出する。予め与えられる条件によって、各傭船に対して「積港」及び「積銘柄」が固定化されている場合は、船舶別・積港別・揚港別フレート（図１５を参照）を用いる。つまり、これらの傭船は原材料を積載する傭船が決定されているため、船舶別・積港別・揚港別フレートを用いることで、最適化によって原材料を入荷する製鉄所が決定された時点で、正確な輸送費用計算が可能となる。

また、いずれの項目も固定化されていない場合や、「積港」だけが固定化されている場合は、銘柄別・揚港別見做しフレートを用いる。つまり、積港と積銘柄が決定されていない場合には、当該船舶に関する積港と積銘柄自体を変更可能にすることで、後述の最適化を用いて、当該積港と積銘柄より輸送費用の安い積港と積銘柄の有無を検討することが可能となるこの場合は、銘柄別・揚港別見做しフレートを用いることで、当該船舶に関して、当該積港と積銘柄より輸送費用の安い積港と積銘柄に、当該船舶の積港と積銘柄を変更させることを後述する最適化により計画させる。これにより輸送費用のより安い計画を作成することを可能としている。なお、同一の傭船に関しては、固定化が最もされていないレコードの状態をこの傭船の固定化状況と考える。この固定化抽出処理は、図５に示したタイミングである必要はなく、例えば配合計画作成を開始するときに行われるようにしてもよい。

（９）配合計画数式モデルを目的関数に基づいて最適化（図３の計画部３１５、図４の配合計画求解部３５９、図５のステップＳ３０９）
上記構築された線形及び整数制約式でなる需給バランスモデル、性状モデルを併せて配合計画数式モデルとし、予め構築した目的関数に基づきＬＰ（線形計画法）、ＭＩＰ（混合整数計画法）、ＱＰ（２次計画法）等の数理計画法により最適化問題として問題を解くことにより、最適な使用量、入荷量を計算する。

ここでは、目的関数に関して線形式を用いた場合の例を示す。本実施形態では、費用（原材料の購入費用及び輸送費用）のミニマム化を目的としており、目的関数Ｊの一例を（式４０）に示す。目的関数を用いて求解するに際して、購入費用情報及びステップＳ３０８において設定された輸送費用情報を用いる。

なお、（式４０）は目的関数の一例であり、他の目的関数に替えたり、他の目的関数を加えたりしてもよい。例えば、与えられた目標とする配合割合に近い配合割合に配合計画を近づける必要があり、更に前日の配合割合とその翌日の配合割合が大きく乖離することがない配合計画を作成する必要がある場合を考える。この場合は、使用目標量からの溢れ量、不足量、及び当該日の使用量と当該日前日の使用量との差をミニマム化する項目を目的関数に追加する。

以上の定式化した式（数式モデル）を展開、設定して、混合整数計画法にて解くことにより、需給バランスモデル、性状モデルを併せた配合計画数式モデルに対する最適解が得られる。

以上説明した配合計画求解部３５９（計画部３１５）及びステップＳ３０９が、本発明でいう第２の配合計画作成部の最適化計算部及びそれによる処理の例である。

（１０）最適化計算による求解結果の判定（図４の求解結果判定部３６０、図５のステップＳ３１０、Ｓ３１１）
（式１４）´を用いた最適化計算による求解結果が、非線形の数式を含む数式モデルｆ（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓを満たすか否かを判定する。その結果、非線形の数式を含む数式モデルｆ（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓを満たせば、この求解結果に従い、後述する性状シミュレータ３６２に対する計算指示を作成して、シミュレーションを実行させる。非線形の数式を含む数式モデルｆ（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓを満たさなければ、線形の数式を含む数式モデルｆ´（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓ´を調整する（図５のステップＳ３１１）。具体的には、仮下限値Ｓ´を微増させる。

図１０は、ステップＳ３０７〜Ｓ３１０の処理、すなわち非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に代えて線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）を導入したときの処理を示すフローチャートである。ステップＳ４０１において、需給バランスモデル、性状モデル（非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に代えて線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）を導入して構築したもの）、目的関数Jに基づいて最適化計算を実行する。

この場合に、（式１４）´に示したように、線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）に対する下限値として、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に対する下限値Ｓよりも小さな仮下限値Ｓ´＝Ｓ−ｓ（ｓ：オフセット値）を設定する。

次にステップＳ４０２において、線形の数式を含む数式モデルｆ´（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓ´を用いた最適化計算による求解結果が、非線形の数式を含む数式モデルｆ（ｘA、ｘB、ｘC、・・・、ｘN）≧Ｓを満たすか否かを判定する。すなわち、ステップＳ４０１の最適化計算による求解結果（各銘柄Ａ〜Ｎの使用量（配合割合））を（式１４）に代入し、（式１４）が成立するか否かを判定する。

ステップＳ４０２の結果、（式１４）が成立すれば、本処理を終了する（図１０のステップＳ４１２に移行する）。それに対して、（式１４）が成立しなければ、ステップＳ４０３に進んで、仮下限値Ｓ´を予め設定された増減幅で微増させて、再度ステップＳ４０１の処理を実行する。すなわち、（式１４）が成立するまで、仮下限値Ｓ´を微増させて、最適化計算による求解を繰り返す収束計算を実行する。

なお、本実施形態では、混合後の性状制約が下限値を有する場合を例にして説明したが、上限値を有する場合も同様である。この場合、線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）は、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）の上限をなすものを考える。また、ステップＳ４０１では、線形の数式ｆ´（ｘA、ｘB、ｘC、・・・、ｘN）に対する上限値として、非線形の数式ｆ（ｘA、ｘB、ｘC、・・・、ｘN）に対する上限値よりも大きな仮上限値を設定する。

（１１）求解した解に基づいて在庫推移をシミュレーション（図３の在庫推移シミュレータ３１１、図４の在庫推移シミュレータ３６１、図５のステップＳ３１２）
上記配合計画数式モデルに対する解、及び、入力データ取込み部３５１により取込まれたデータの全部又は一部に基づいて、配合の全部或いは一部を対象として、設定した計画確定期間分について、設定した計画作成精度でシミュレーションを実行する。このシミュレーションでは、配合計画数式モデルには組込むことができなかった制約条件、例えば一定の規則に基づかない条件など、定式化が難しいもの、及び、操業のルール等も組み込んでシミュレートする。このことで、配合計画数式モデルに対する求解結果として出された解を実操業で問題なく使用可能な配合計画に変更する。これにより、実操業で求められる時間精度と、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

また、数式モデルでは取扱うことが難しい制約の一例として、配合割合が変わった場合の設備の段取りに掛かる段取時間等をシミュレーションに取込み、正確にシミュレートすることで、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

（１２）求解した解に基づいて性状をシミュレーション（図３の性状シミュレータ３１２、図４の性状シミュレータ３６２、図５のステップＳ３１３）
上記配合計画数式モデルに対する解、在庫推移シミュレータ３６１によりシミュレーションされた在庫推移、及び、入力データ取込み部３５１により取込まれたデータの全部或いは一部に基づいて、配合の全部或いは一部を対象として、設定した計画確定期間分について、設定した計画作成精度で、性状のシミュレートをおこなう。シミュレーションの結果として、原材料の混合後の性状結果が得られる。このシミュレーションでは、配合計画数式モデルには組み込むことができなかった制約条件、操業のルール等も組み込んでシミュレートすることで、配合計画数式モデルに対する求解結果として出された解を実操業で問題なく使用可能な配合計画に変更する。これにより、実操業で求められる時間精度と、実操業に求められる細かな制約まで考慮した配合計画の立案が可能となる。

以上説明した在庫推移シミュレータ３６１（在庫推移シミュレータ３１１）及びステップＳ３１２、並びに、性状シミュレータ３６２（性状シミュレータ３１２）及びステップＳ３１３が、本発明でいう第２の配合計画作成部のシミュレータ及びそれによる処理の例である。

（１３）配合計画の確定（図４の確定部３６３、図５のステップＳ３１４）
上記在庫推移シミュレーション、性状シミュレーションにより導き出された配合計画のうちで設定した計画確定期間分を確定する。図７に示すように、本実施形態では計画確定期間を１日と設定しているので、作成した配合計画の最初の１日分を確定する。作成した配合計画のうちで上記計画確定期間に入らなかった部分については、その計画は確定せずに破棄する。

（１４）計画作成期間分、或いは計画確定期間分の計画が確定したか判定（図４の判定部３６４、図５のステップＳ３１５）
このステップの実行時点までに確定した計画確定期間が予め設定した計画作成期間の全体を含んでいるかを判断する。本実施形態では、計画作成期間が１０日間であるので、第１０ループで計画を確定した時点で計画確定期間分の計画が確定する。このため第１０ループで計画を確定終了した時点で１０日分の配合計画を作成して、処理を終了する。

（１５）立案開始日の更新（図４の更新部３６５、図５のステップＳ３１６）
このステップの実行時点で確定した計画確定期間が予め設定した計画作成期間の全体を含んでいない場合、上記配合計画のうちで確定した配合計画期間直後の日時を新たな立案開始日として設定する。本実施形態では、図７に示すように、第１ループでは当初１日目０時であった立案開始日を２日目０時に、第２ループでは当初２日目０時であった立案開始日を３日目０時に更新する。

（１６）配合計画の出力（図４の出力部３６６、図５のステップＳ３１７）
以上のようにして作成した配合計画は、出力部３６６により、表示部３０３に画面表示されたり、データベース４００を含む外部機器にデータ送信されたりする。

以上説明した出力部３６６及びステップＳ３１７が、本発明でいう第２の配合計画作成部の出力部及びそれによる処理の例である。

以上のように、現在の在庫推移状態に応じて、需給バランス制約、性状制約について、まず所定の最適化期間分について、計画作成時間精度で数式モデルを設定し、設定した配合計画数式モデルを目的関数に基づいて求解し、求解した解に基づいて、在庫推移、混合後の性状をシミュレートし、シミュレーション結果から求められた配合計画のうちで、設定した計画確定期間分を確定し、計画確定期間直後の日時を新たな立案開始日時とすることにより、新たな計画対象期間分の配合計画を確定する一連の処理を順次、予め定めた回数だけ、繰り返して実行する。このことで、所望する計画作成期間分の配合計画を作成することができる。これにより、任意の時間精度を必要とする配合計画を高速かつ詳細に最適化することができ、しかも得られた結果をそのままで実操業に適用できる。

なお、上記実施形態では、図１１に示すように、配合計画は一定の期間（例えば旬）毎に作成される。また、複数の性状α、βについて性状モデルが非線形となることがある。図１１において、Ａは性状制約を満たしている（（式１４）が成立している）ことを、Ｂは性状制約を満たしていないことを意味する。すなわち、図１１の例では、性状αについて複数旬（４月上旬及び下旬）で性状違反が発生しており、同様に性状βについて複数旬（４月上旬及び下旬）で性状違反が発生している。

この場合に、各旬及び各性状について別個に図１０で説明した収束計算を行うと、以下の問題が生ずる。具体的にいえば、４月上旬で性状αについて収束計算を行い、続いて性状βについて収束計算を行い、また、４月下旬で性状αについて収束計算を行い、続いて性状βについて収束計算を行うのでは、計算処理に時間がかかってしまう。

そこで、第１の実施形態でも説明したように、対象の旬及び性状についてまとめて図１０で説明した収束計算を行うようにする。例えば４月上旬及び下旬で性状α、βについてまとめて収束計算を行う（図１０のステップＳ４０３で性状α、βの仮下限値の微増（或いは仮上限値の微減）を同時に行う）ことにより、高速化を図ることができる。

また、上記実施形態では、図１０のステップＳ４０３で仮下限値Ｓ´を微増（或いは仮上限値を微減）させた後、再度ステップＳ４０１の処理を実行すると説明した。この場合に、収束計算が変化しても設定に変化のない数式モデル、具体的には上述した需給バランスモデルや元々線形の性状モデルは保持しておく。そして、仮下限値を微増（或いは仮上限値を微減）させて再度ステップＳ４０１の処理を実行する場合に、収束計算に従って変化のある数式モデル、具体的には仮下限値を微増させた（或いは仮上限値を微減させた）数式モデルのみ変更するような仕組とすることにより、高速化を図ることができる。

また、配合計画（例えば使用量（配合割合））として、年次計画、期計画、月次計画といった長期間の計画を立案することが多い。このように長期の配合計画を予め作成し、その配合計画を基準の配合計画とし、第１の配合計画作成装置１００により作成されたより短期の配合計画が、基準となる配合計画から一定幅以上離れないようにすることも重要となる。

そこで、（式１７）に示したような費用（原材料の購入費用及び輸送費用）に関して構築された目的関数Jに加え、予め作成された基準となる配合計画と一定幅以上離れないようにすることに関して構築された目的関数J´を用いるようにしてもよい。目的関数Ｊ´の一例を（式２０）に示す。

上記例では、月次計画において、期計画を基準となる配合計画として、日々の配合計画を作成する場合の一例を示した。この場合は、配合割合（銘柄、日）と基準配合割との差の銘柄毎、日毎に合計したものをミニマム化する。他の例として、期計画を立案する場合、年次計画を基準となる配合計画として計画を作成しても良い。この場合、月次計画では配合割合（銘柄、月）を決定するとした場合は、配合割合（銘柄、月）と基準配合割との差を銘柄毎、月毎に合計した値をミニマム化する。

なお、基準となる配合計画は、例えば過去の実績に基づいて作成され、その作成手法はどのようなものであってもよい。もちろん、本発明を適用した配合計画作成手法により長期間の計画を予め作成しておき、それを基準となる配合計画としてもよい。

以上により、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を相互に連係させて一括で作成することが可能となる。また、これら計画においては、配合計画・配船計画を通して全体として輸送費用をミニマム化することが可能となる。更に、第１の配合計画作成装置１００により使用計画を作成し、この使用計画を入力データとして配船計画を作成する場合より、より確実に在庫切れを防ぐことができる配合計画・配船計画を作成できる。

（第３の実施形態）
上記第２の実施形態では、配船計画作成装置２００において、各銘柄の在庫に関する制約（式２４）として、個別の銘柄毎に安全在庫量を切らないことがある。しかし、実際の操業においては、性状が近い銘柄（規定した範囲に含まれる化学性質を共通して備える複数の銘柄：互いに置き換えても使用可能な銘柄）は、互いに置き換えて使用することで、融通を利かせている。例えば、原材料Ａ，Ｂの在庫量がそれぞれ５万トン、０トンの場合、使用計画において原材料Ａが０トン、Ｂが２万トンの場合、Ｂ銘柄のみ考慮すれば在庫切れが発生する。しかし、原材料Ａ，Ｂの性状が、互いに置き換えても使用可能な関係である場合は、原材料Ｂを２万トン使用する代わりにＡを２万トン使用する運用が行われている。こうすることで、在庫切れを防止できる。
例えば石炭の場合、「石炭化度が７０％以下；７０％より大きく９０％以下；９０％より大きい」等の、原材料の物理的、または化学的な性状に関する条件を用いて、銘柄のグループ化を行う。 また、このように銘柄をグループ化し、一つのものとして取り扱うことにより、変数を少なくして計算量を減らすことができる。

本実施形態においては、上記運用を実現する配合計画・配船計画の作成を実現することを目的とする。このため、配船計画作成装置２００において、各銘柄の在庫に関する制約を個別の銘柄毎に安全在庫量を切らないものとして取り扱うのではなく、相互に代替可能な性状を持つ銘柄は、グループ化して取り扱う。つまり、グループ化された銘柄においては、グループとしての銘柄の在庫量が、グループとしての安全在庫量を切らないものとして、在庫制約を作成する。

つまり、配船計画作成装置２００において、（式２４）の代わりに、各グループ化された銘柄の在庫量が常にグループとしての安全在庫以上確保されていることを表す制約式を用いる。この制約式は、下記の（式４２）と表される。

このように銘柄をグループ化して取扱うことで、フレートが高い船でしか輸送できない銘柄に変わり、同一銘柄でフレートのより安い船で手配できる銘柄を輸送することが可能となり、輸送費用を抑制できる。例えば、性状がほぼ同一の原材料Ｘ，Ｙがあり、揚港（製鉄所）Ａでは原材料Ｘ，Ｙどちらでの使用も可能な場合を考える。揚港Ａに原材料Ｘを輸送する費用が２０＄／トン、原材料Ｙを輸送する費用４０＄／トン、揚港Ｂに原材料Ｘを輸送する費用が４０＄／トン、原材料Ｙを輸送する費用２０＄／トン、初期在庫量は全て０、他製鉄所との引取の関係等何らかの理由により、揚港Ａで原材料Ｘの使用が０トン、Ｙの使用が５万トン、揚港Ｂで原材料Ｘの使用が５万トン、Ｙの使用が０トンの配合計画が第１の配合計画作成装置１００により作成された場合を考える。この場合、配船計画において個別に銘柄の在庫を維持する場合には、揚港Ａに原材料Ｙを５万トン、揚港Ｂに原材料Ｘを５万トン入荷する配船が、在庫制約より導かれる。この場合は、４０＄／トン×５万トン×２の輸送費用が必要となる。しかし、代替可能な性状を持つ銘柄の群をグループとして考える場合には、揚港Ａに原材料Ｘを５万トン、揚港Ｂに原材料Ｙを５万トン入荷する計画を配船計画で立てた場合も、揚港Ａ，Ｂ共にＸとＹの合計の在庫量が０を下回らないこととなり、在庫制約を満たす。この場合の輸送費用は、２０＄／トン×５万トン×２となり、グループを考慮しない場合に比べて、半分の費用となる。

以上により、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を相互に連係させて一括で作成することが可能となる。また、これら計画においては、配合計画・配船計画を通して全体として輸送費用をミニマム化することが可能となる。更に、第１の配合計画作成装置１００により使用計画を作成し、この使用計画に関する入力データとして、代替可能な性状を持つ銘柄をグループとして考えることにより、在庫を切らさない配船計画を作成する場合と比較して、より輸送費用が安い配合計画・配船計画を作成することが可能となる。

（第４の実施形態）
第３の実施形態において、第１の配合計画作成装置１００により使用計画を作成し、この使用計画を入力データとして、代替可能な性状を持つ銘柄をグループとして考え、在庫を切らさない配船計画を作成する実施形態を説明した。

しかし、実際の操業を行う上では、最終的な配合計画・配船計画としては、個別銘柄毎の在庫が確保されていることが必要となる場合が多い。つまり、上記第３の実施形態において、グループ間での入替を行った場合には、第１の配合計画作成装置１００により作成された使用計画と、配船計画作成装置により作成された配船計画では、代替可能な銘柄グループとしての在庫切れは発生しない。しかし、個別の銘柄としては、銘柄間の使用と入荷にずれが生じるため、個別銘柄で在庫を見た場合に在庫が切れているように見え、操業者が扱い難いという問題が発生する。この問題を解決するための実施形態として、配船計画作成装置２００により作成された配船計画を入力データとして、第２の配合計画作成装置３００により、配合計画を作成することで、入荷に沿った配合計画を作成することが可能となる。この場合、代替可能な性状を持つ銘柄間での入替が行われているため、第２の配合計画作成装置３００においても、性状制約に無理なく、配合計画を作成することが可能となる。更に、第１の配合計画作成装置１００では、船舶が決定していない状態で、引取目標量を基に、使用計画を作成しているため、使用計画としては精度が良くない。しかし、第２の配合計画作成装置３００では、配船計画作成装置により船舶が決定した状態、つまり入荷が決定した状態で使用計画を作成するため、より精度が良い使用計画が作成できる。

以上の構成または方法により、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を相互に連係させて一括で作成することが可能となる。また、これら計画においては、配合計画・配船計画を通して全体として輸送費用をミニマム化することが可能となる。また、第１の配合計画作成装置１００により使用計画を作成し、この使用計画を入力データとして配船計画を作成する場合と比較して、より輸送費用が安い配合計画・配船計画を作成することが可能となる。さらに、第２の配合計画作成装置３００により使用計画を作成することで、より精度の良い使用計画を作成することが可能となる。

図２７には、本発明の配合計画装置１００、３００或いは配船計画作成装置２００として機能しうるコンピュータ装置２５００のハードウェア構成例を示す。装置全体を制御する中央処理装置であるＣＰＵ２５０１、各種入力条件や結果等を表示する表示部２５０２、結果等を保存するハードディスク等の記憶部２５０３、制御プログラム、各種アプリケーションプログラム、データ等を記憶するＲＯＭ（リードオンリーメモリ）２５０４、ＣＰＵ２５０１が処理を行うときに用いる作業領域であるＲＡＭ（ランダムアクセスメモリ）２５０５、及びキーボード、マウス等の入力部２５０６等により構成される。

また、上述した実施形態の機能を実現するべく各種のデバイスを動作させるように、この各種デバイスと接続された装置或いはシステム内のコンピュータに対し、上記実施形態の機能を実現するためのソフトウェアのプログラムコードを供給し、そのシステム或いは装置のコンピュータ（ＣＰＵ或いはＭＰＵ）に格納されたプログラムに従って上記各種デバイスを動作させることによって実施したものも、本発明の範疇に含まれる。この場合、上記ソフトウェアのプログラムコード自体が上述した実施形態の機能を実現することになり、そのプログラムコード自体、及びそのプログラムコードをコンピュータに供給するための手段、例えば、かかるプログラムコードを格納した記録媒体は本発明を構成する。プログラムコードを記憶する記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。

本発明は、以下の態様を更に含む。
（１） 複数の供給元から船舶にて複数の供給先に輸送されて入荷した複数銘柄の配合原材料を、各供給先において配合して使用する際に、当該複数銘柄の配合原材料の配合計画、及び当該複数銘柄の配合原材料を、積地である複数の供給元から揚地である複数の供給先に輸送する船舶の配船計画を作成する配合及び配船計画作成システムであって、配合原材料の、供給元毎及び銘柄毎に予め設定された引取目標量に基づいて、複数銘柄の配合原材料を入荷して混合する配合計画を作成する第１の配合計画作成手段と、前記第１の配合計画作成手段により作成された配合計画に基づいて、複数銘柄の配合原材料を複数の積地から複数の揚地に輸送する配船計画を作成する配船計画作成手段と、前記配船計画作成手段により作成された配船計画を格納するデータベース手段とを備えたことを特徴とする配合及び配船計画作成システム。
（２） 前記第１の配合計画作成手段は、配合原材料の引取目標量、配合原材料の在庫状況、配合原材料の性状、配合原材料の購入費用情報、船舶を利用する際の輸送費用情報を含むデータを取り込むデータ取込み手段と、配合原材料の需給状態及び混合後の性状を計算するシミュレータと、配合原材料の需給バランス制約を表す数式モデル、及び、混合後の性状制約を表す数式モデルを構築する数式モデル構築手段と、前記数式モデル構築手段により構築された数式モデルを用いて、配合原材料の購入費用及び輸送費用に関して構築された目的関数に基づいて最適化計算を行い、前記シミュレータに対する指示を算出する最適化計算手段と、前記シミュレータによるシミュレーション結果である配合計画を出力する出力手段とを備えたことを特徴とする上記（１）に記載の配合及び配船計画作成システム。
（３） 前記第１の配合計画作成手段の最適化計算手段では、更に、配合原材料の入荷量と引取目標量との関係に関して構築された目的関数に基づいて最適化計算を行うことを特徴とする上記（２）に記載の配合及び配船計画作成システム。
（４） 前記第１の配合計画作成手段のデータ取込み手段により取り込む輸送費用情報には、船舶別・積港別・揚港別フレートの情報と、銘柄別・揚港別フレートの情報とが含まれることを特徴とする上記（２）又は（３）に記載の配合及び配船計画作成システム。
（５） 前記第１の配合計画作成手段は、作成された配合計画による配合原材料の使用予定量を算出するものであり、前記配船計画作成手段は、前記第１の配合計画作成手段により作成された配合計画による配合原材料の使用予定量、引取目標量、契約の種別の異なる船舶がリストアップされた船舶リスト、前記船舶リストにリストアップされている船舶の運航状況、配合原材料の在庫状況、配合原材料の購入費用情報、前記船舶リストにリストアップされている船舶を利用する場合の輸送費用情報を含むデータを取り込むデータ取込み手段と、配合原材料の在庫推移を計算する在庫推移シミュレータ、及び、船舶運航状況の推移を計算する船舶運航状況推移シミュレータにより構成されるシミュレータと、前記船舶運航状況に基づいて前記船舶リストから船舶を選択し、必要な船舶財源を作成する船舶財源作成手段と、少なくとも前記船舶財源作成手段により作成された船舶の運航制約、揚地での配合原材料の需給バランス制約、及び引取目標量制約を表す数式モデルを構築する数式モデル構築手段と、前記数式モデル構築手段により構築された数式モデルを用いて、少なくとも輸送費用に関して構築された目的関数に基づいて最適化計算を行い、前記シミュレータに対する指示を算出する最適化計算手段と、前記シミュレータによるシミュレーション結果である配船計画を出力する出力手段とを備えたことを特徴とする上記（１）乃至（４）のいずれか１つに記載の配合及び配船計画作成システム。
（６） 前記配船計画作成手段は、前記第１の配合計画作成手段により作成された配合計画による配合原材料の使用予定量、引取目標量、契約の種別の異なる船舶がリストアップされた船舶リスト、前記船舶リストにリストアップされている船舶の運航状況、前記性状が近い銘柄をグループ化して取り扱った配合原材料の在庫状況、配合原材料の購入費用情報、前記船舶リストにリストアップされている船舶を利用する場合の輸送費用情報を含むデータを取り込むデータ取込み手段と、前記性状が近い銘柄をグループ化して取り扱う配合原材料の在庫推移を計算する在庫推移シミュレータ、及び、船舶運航状況の推移を計算する船舶運航状況推移シミュレータにより構成されるシミュレータと、前記船舶運航状況に基づいて前記船舶リストから船舶を選択し、必要な船舶財源を作成する船舶財源作成手段と、少なくとも前記船舶財源作成手段により作成された船舶の運航制約、揚地での前記性状が近い銘柄をグループ化して取り扱う配合原材料の需給バランス制約、及び引取目標量制約を表す数式モデルを構築する数式モデル構築手段と、前記数式モデル構築手段により構築された数式モデルを用いて、少なくとも輸送費用に関して構築された目的関数に基づいて最適化計算を行い、前記シミュレータに対する指示を算出する最適化計算手段と、前記シミュレータによるシミュレーション結果である配船計画を出力する出力手段とを備えたことを特徴とする上記（１）乃至（４）のいずれか１つに記載の配合及び配船計画作成システム。
（７） 前記配船計画作成手段により作成された配船計画に基づいて、複数銘柄の配合原材料を入荷して混合する配合計画を作成する第２の配合計画作成手段と、前記第２の配合計画作成手段により作成された配合計画を格納するデータベース手段とを備えたことを特徴とする上記（１）乃至（６）のいずれか１つに記載の配合及び配船計画作成システム。
（８） 前記第２の配合計画作成手段は、前記配船計画作成手段により作成された配船計画による配合原材料の入荷予定、配合原材料の在庫状況、配合原材料の性状、配合原材料の購入費用情報、船舶を利用する際の輸送費用情報を含むデータを取り込むデータ取込み手段と、配合原材料の需給状態及び混合後の性状を計算するシミュレータと、配合原材料の需給バランス制約を表す数式モデル、及び、混合後の性状制約を表す数式モデルを構築する数式モデル構築手段と、前記数式モデル構築手段により構築された数式モデルを用いて、配合原材料の購入費用及び輸送費用に関して構築された目的関数に基づいて最適化計算を行い、前記シミュレータに対する指示を算出する最適化計算手段と、前記シミュレータによるシミュレーション結果である配合計画を出力する出力手段とを備えたことを特徴とする上記（７）に記載の配合及び配船計画作成システム。
（９） 前記第１の配合計画作成手段により作成された配合計画を、性状が近い銘柄をグループ化して取り扱うことを特徴とする上記（１）乃至（８）のいずれか１つに記載の配合及び配船計画作成システム。
（１０） 複数の供給元から船舶にて複数の供給先に輸送されて入荷した複数銘柄の配合原材料を、各供給先において配合して使用する際に、当該複数銘柄の配合原材料の配合計画、及び当該複数銘柄の配合原材料を、積地である複数の供給元から揚地である複数の供給先に輸送する船舶の配船計画を作成する配合及び配船計画作成方法であって、第１の配合計画作成手段が、配合原材料の、供給元毎及び銘柄毎に予め設定された引取目標量に基づいて、複数銘柄の配合原材料を入荷して混合する配合計画を作成するステップと、配船計画作成手段が、前記第１の配合計画作成手段により作成された配合計画に基づいて、複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を作成するステップと、前記配船計画作成手段により作成された配船計画をデータベース手段に格納するステップとを有することを特徴とする配合及び配船計画作成方法。
（１１） 上記（１）乃至（９）のいずれか１つに記載の配合及び配船計画作成システムの各手段としてコンピュータを機能させるためのプログラム。

本発明によれば、複数銘柄の配合原材料を入荷して混合する配合計画、及び複数銘柄の原材料を複数の積地から複数の揚地に輸送する配船計画を、相互に連係させ、一括して作成することができる。特に、傭船契約の種別の異なる船舶を考慮して、輸送費用に関して構築した目的関数を用意し、この目的関数に基づいて最適化計算を行うことができる。このことにより、船舶の傭船契約の種類（連続航海船、不定期船、スポット船）、船団の構成を決める船舶を雇うか、雇わないかまでも含めた、輸送費用のミニマム化のための計画立案が可能になる。更には、配合計画段階においても、輸送費用をミニマム化することを考慮した計画の作成が可能となる。

１００：第１の配合計画作成装置
２００：配船計画作成装置
２０１：シミュレータ
２０２：マクロ最適化部
２０２ａ：船舶財源リスト作成部
２０２ｂ：数式モデル設定部
２０２ｃ：最適化計算部
２０３：ミクロ最適化部
２０３ａ：数式モデル設定部
２０３ｂ：最適化計算部
２０４：データ取込み部
２０５：出力部
３００：第２の配合計画作成装置
３０３：表示部
３０４：操業者評価部
３１２：性状シミュレータ
３１３：需給バランスモデル設定部
３１４：性状モデル設定部
３１５：計画部
３５１：入力データ取込み部
３５２：計画作成期間設定部
３５３：時間精度設定部
３５４：最適化期間設定部
３５５：計画確定期間設定部
３５６：需給バランスモデル設定部
３５７：性状モデル設定部
３５７ａ：線形化部
３５８：固定化抽出処理部
３５９：配合計画求解部
３６０：確認部
３６１：在庫推移シミュレータ
３６２：性状シミュレータ
３６３：確定部
３６４：判定部
３６５：更新部
３６６：出力部
４００：データベース
５００：コンピュータ

（１） 上記目的を達成するための、本発明の一態様にかかるシステムは、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
を作成する配合及び配船計画作成システムであって：
前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成部と；
作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成部と；
作成された前記配船計画を格納するデータベース部と；
を備え、
前記第１の配合計画作成部は：
前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込むデータ取込み部と；
前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する数式モデル設定部と；
設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートするシミュレータと；
前記シミュレータによるシミュレーション結果である配合計画を出力する出力部と；
を備え、
前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
前記配船計画作成部は：
前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；
前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；
前記船舶財源リストに含まれる前記船舶の運航制約、前記揚地での前記配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；
設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；
前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；
を備えることを特徴とする配合及び配船計画作成システムである。
（２） 上記目的を達成するための、本発明の別の一態様にかかる方法は、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
を作成する配合及び配船計画作成システムであって：
前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成部と；
作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成部と；
作成された前記配船計画を格納するデータベース部と；
を備え、
前記第１の配合計画作成部は：
前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込むデータ取込み部と；
前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する数式モデル設定部と；
設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートするシミュレータと；
前記シミュレータによるシミュレーション結果である配合計画を出力する出力部と；
を備え、
前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
前記配船計画作成部は：
前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われたグループ化配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；
前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；
前記船舶財源リストに含まれる前記船舶の運航制約、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；
設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；
前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；
を備えることを特徴とする配合及び配船計画作成システムである。
（３） 上記（１）の配合及び配船計画作成システムでは、前記配合計画の作成と、前記配船計画の作成と、において、前記引取目標量および前記配合原材料の在庫状況に関して共通の制約条件が用いられてもよい。
（４） 上記（１）または（２）の配合及び配船計画作成システムにおいて、前記最適化計算部が、前記配合原材料の入荷量と前記引取目標量との関係に関して予め構築された目的関数に更に基づいて最適化計算を行ってもよい。
（５） 上記（１）または（４）の配合及び配船計画作成システムにおいて、前記データ取込み部により取り込まれる前記輸送費用には；船舶別・積港別・揚港別のフレートの情報と；前記銘柄別・揚港別フレートの情報と；が含まれてもよい。
（６） 上記（２）の配合及び配船計画作成システムは、 作成された前記配船計画に基づいて、第２の配合計画を作成する第２の配合計画作成部を更に備えてもよく、前記データベース部は前記第２の配合計画を更に格納してもよい。
（７） 上記（６）の配合及び配船計画作成システムにおいて、前記第２の配合計画作成部は：作成された前記配船計画に基づく前記配合原材料の入荷予定、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込む第２のデータ取込み部と；前記入荷予定を用い、前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する第２の数式モデル設定部と；設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う第２の最適化計算部と；前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートする第２のシミュレータと；前記シミュレータによるシミュレーション結果である配合計画を出力する第２の出力部と；を備えてもよい。
（８） 上記（７）の配合及び配船計画作成システムでは、前記第１の配合計画の前記需給バランス制約と、前記第２の配合計画の前記需給バランス制約と、において、需要に関して同一の制約条件を設定し；前記第１の配合計画の前記需給バランス制約において、供給に関して、引取目標量を制約条件として設定し；前記第２の配合計画の前記需給バランス制約において、供給に関して、船舶による原材料の入荷量を制約条件として設定してもよい。
（９） 上記（７）の配合及び配船計画作成システムでは、前記第１の配合計画の前記目的関数において、銘柄別・揚港別見做しフレートを用い；前記第２の配合計画の前記目的関数において、船舶別・積港別・揚港別フレートを用いてもよい。
（１０） 上記（１）または（２）の配合及び配船計画作成システムでは、作成された前記配合計画において、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われてもよい。
（１１） 上記目的を達成するための、本発明の別の一態様にかかる方法は、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
を作成する配合及び配船計画作成方法であって：
前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画作成部が第１の配合計画を作成する第１の配合計画作成工程と；
作成された前記第１の配合計画に基づいて、配船計画作成部が前記配船計画を作成する配船計画作成工程と；
作成された前記配船計画をデータベースに格納するデータベース格納工程と；
を備え、
前記第１の配合計画作成工程は：
前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータをデータ取込み部が取り込むデータ取込み工程と；
前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ数式モデル設定部が設定する数式モデル設定工程と；
設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算部が最適化計算を行う最適化計算工程と；
前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレータがシミュレートするシミュレーション工程と；
前記シミュレーション工程によるシミュレーション結果である配合計画を出力部が出力する出力工程と；
を備え、
前記第１の配合計画作成工程は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
前記配船計画作成工程は：
前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータをデータ取込み部が取り込むデータ取込み工程と；
前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを船舶財源リスト作成部が作成する船舶財源リスト作成工程と；
前記船舶財源リストに含まれる前記船舶の運航制約、前記揚地での前記配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを数式モデル設定部が設定する数式モデル設定工程と；
設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算部が最適化計算を行う最適化計算工程と；
前記在庫状況の推移を在庫推移シミュレータがシミュレートする在庫推移シミュレーション工程と、前記運航状況の推移を船舶運航状況推移シミュレータがシミュレートする船舶運航状況推移シミュレーション工程と、を含み、前記最適化計算の結果に基づいて動作する、シミュレーション工程と；
前記シミュレーション工程によるシミュレーション結果である前記配船計画を出力部が出力する出力工程と；
を備える。
（１２） 上記目的を達成するための、本発明の別の一態様にかかる方法は、複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
を作成する配合及び配船計画作成方法であって：
前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画作成部が第１の配合計画を作成する第１の配合計画作成工程と；
作成された前記第１の配合計画に基づいて、配船計画作成部が前記配船計画を作成する配船計画作成工程と；
作成された前記配船計画をデータベースに格納するデータベース格納工程と；
を備え、
前記第１の配合計画作成工程は：
前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータをデータ取込み部が取り込むデータ取込み工程と；
前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ数式モデル設定部が設定する数式モデル設定工程と；
設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算部が最適化計算を行う最適化計算工程と；
前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレータがシミュレートするシミュレーション工程と；
前記シミュレーション工程によるシミュレーション結果である配合計画を出力部が出力する出力工程と；
を備え、
前記第１の配合計画作成工程は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
前記配船計画作成工程は：
前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われたグループ化配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータをデータ取込み部が取り込むデータ取込み工程と；
前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを船舶財源リスト作成部が作成する船舶財源リスト作成工程と；
前記船舶財源リストに含まれる前記船舶の運航制約、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを数式モデル設定部が設定する数式モデル設定工程と；
設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を最適化計算部が行う最適化計算工程と；
性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の前記在庫状況の推移を在庫推移シミュレータがシミュレートする在庫推移シミュレーション工程と、前記運航状況の推移を船舶運航状況推移シミュレータがシミュレートする船舶運航状況推移シミュレーション工程と、を含み、前記最適化計算の結果に基づいて動作する、シミュレーション工程と；
前記シミュレーション工程によるシミュレーション結果である前記配船計画を出力部が出力する出力工程と；
を備える。
（１３） 上記目的を達成するための、本発明の別の一態様にかかるプログラムは、上記（１）または（２）の配合及び配船計画作成システムの各部分としてコンピュータを機能させるためのプログラムである。

Claims (14)

1. 複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
   入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
   を作成する配合及び配船計画作成システムであって：
   前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成部と；
   作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成部と；
   作成された前記配船計画を格納するデータベース部と；
   を備えることを特徴とする配合及び配船計画作成システム。
2. 前記配合計画の作成と、前記配船計画の作成と、において、前記引取目標量および前記配合原材料の在庫状況に関して共通の制約条件が用いられることを特徴とする請求項１に記載の配合及び配船計画作成システム。
3. 前記第１の配合計画作成部は：
   前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込むデータ取込み部と；
   前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する数式モデル設定部と；
   設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
   前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートするシミュレータと；
   前記シミュレータによるシミュレーション結果である配合計画を出力する出力部と；
   を備えることを特徴とする請求項１に記載の配合及び配船計画作成システム。
4. 前記最適化計算部が、前記配合原材料の入荷量と前記引取目標量との関係に関して予め構築された目的関数に更に基づいて最適化計算を行うことを特徴とする請求項３に記載の配合及び配船計画作成システム。
5. 前記データ取込み部により取り込まれる前記輸送費用には；
   船舶別・積港別・揚港別のフレートの情報と；
   前記銘柄別・揚港別フレートの情報と；
   が含まれることを特徴とする請求項３又は４に記載の配合及び配船計画作成システム。
6. 前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
   前記配船計画作成部は：
   前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、前記配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；
   前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；
   前記船舶財源リストに含まれる前記船舶の運航制約、前記揚地での前記配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；
   設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
   前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；
   前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；
   を備えることを特徴とする請求項１に記載の配合及び配船計画作成システム。
7. 前記第１の配合計画作成部は、更に、作成された前記配合計画に従った前記配合原材料の使用予定量を算出し、
   前記配船計画作成部は：
   前記配合原材料の前記使用予定量、前記配合原材料の前記引取目標量、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われたグループ化配合原材料の在庫状況、前記配合原材料の購入費用、複数の種別の傭船契約に基づいて運用される複数の船舶がリストアップされた船舶リスト、各々の前記船舶の運航状況、及び、輸送費用、を含むデータを取り込むデータ取込み部と；
   前記運航状況に基づいて前記船舶リストから配船対象の候補となる前記船舶を選択し、船舶財源リストを作成する船舶財源リスト作成部と；
   前記船舶財源リストに含まれる前記船舶の運航制約、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の需給バランス制約、及び、引取目標量制約を表す数式モデルを設定する数式モデル設定部と；
   設定された前記数式モデルを用いて、前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う最適化計算部と；
   性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われた前記グループ化配合原材料の前記在庫状況の推移をシミュレートする在庫推移シミュレータと、前記運航状況の推移をシミュレートする船舶運航状況推移シミュレータと、を含み、前記最適化計算の結果に基づいて動作する、シミュレータと；
   前記シミュレータによるシミュレーション結果である前記配船計画を出力する出力部と；
   を備えることを特徴とする請求項１に記載の配合及び配船計画作成システム。
8. 作成された前記配船計画に基づいて、第２の配合計画を作成する第２の配合計画作成部を更に備え、
   前記データベース部は前記第２の配合計画を更に格納することを特徴とする請求項１，３，７のいずれか一項に記載の配合及び配船計画作成システム。
9. 前記第２の配合計画作成部は：
   作成された前記配船計画に基づく前記配合原材料の入荷予定、前記配合原材料の在庫状況、前記配合原材料の性状、前記配合原材料の購入費用、及び、前記配合原材料の輸送費用を含むデータを取り込む第２のデータ取込み部と；
   前記入荷予定を用い、前記配合原材料の需給バランス制約、及び、前記配合原材料の混合後の性状制約を表す数式モデルをそれぞれ設定する第２の数式モデル設定部と；
   設定された前記数式モデルを用いて、前記購入費用及び前記輸送費用に関して予め構築された目的関数に基づいて最適化計算を行う第２の最適化計算部と；
   前記最適化計算の結果に基づいて動作し、前記配合原材料の需給状態の推移、及び、前記配合原材料の混合後の前記性状の推移をシミュレートする第２のシミュレータと；
   前記シミュレータによるシミュレーション結果である配合計画を出力する第２の出力部と；
   を備えることを特徴とする請求項８に記載の配合及び配船計画作成システム。
10. 前記第１の配合計画の前記需給バランス制約と、前記第２の配合計画の前記需給バランス制約と、において、需要に関して同一の制約条件を設定し；
    前記第１の配合計画の前記需給バランス制約において、供給に関して、引取目標量を制約条件として設定し；
    前記第２の配合計画の前記需給バランス制約において、供給に関して、船舶による原材料の入荷量を制約条件として設定する；
    ことを特徴とする請求項９に記載の配合及び配船計画作成システム。
11. 前記第１の配合計画の前記目的関数において、銘柄別・揚港別見做しフレートを用い；
    前記第２の配合計画の前記目的関数において、船舶別・積港別・揚港別フレートを用いる、ことを特徴とする請求項９に記載の配合及び配船計画作成システム。
12. 作成された前記配合計画において、性状が規定した範囲に含まれる複数の前記銘柄の前記配合原材料がグループ化されて取り扱われることを特徴とする請求項１又は３に記載の配合及び配船計画作成システム。
13. 複数の銘柄の配合原材料を、複数の積地から複数の揚地に輸送する配船計画と；
    入荷した前記配合原材料を、それぞれの前記揚地で混合する配合計画と；
    を作成する配合及び配船計画作成方法であって：
    前記積地毎及び前記銘柄毎に予め設定された前記配合原材料の引取目標量に基づいて、第１の配合計画を作成する第１の配合計画作成工程と；
    作成された前記第１の配合計画に基づいて、前記配船計画を作成する配船計画作成工程と；
    作成された前記配船計画をデータベースに格納するデータベース格納工程と；
    を備えることを特徴とする配合及び配船計画作成方法。
14. 請求項１に記載の配合及び配船計画作成システムの各部分としてコンピュータを機能させるためのプログラム。