(57)【要約】【課題】制約条件を基に解を探索する順番を設定し、その後制約条件の終点条件内若しくは終点条件近傍で解を探索して解の中から最適解を決定する最適解探索方法を提供する。【解決手段】本発明に係る最適解探索方法は、制約条件を満たす入力条件に対して最適化計算を行って最適解を求める最適解探索方法において、制約条件は始点条件と終点条件を有し、入力ステップで入力条件を入力し、初期条件設定ステップで始点条件と終点条件を設定し、それらの条件の組み合わせを踏まえて探索パターンを作成し、その探索パターンの順番と最初の探索パターンと探索点を設定し、最適化計算ステップで解の探索及び保存をし、探索パターン設定ステップで新たな解とこれまでの最適解との差を算出し、その差を基に次回の探索パターン及び次回探索点を設定し、解探索条件を満たすと最適解決定ステップで保存された解の中から最適解を決定し、出力ステップで最適解を出力する。【選択図】図１
2014209273.tif 000002
【特許請求の範囲】
【請求項１】
制約条件を満たすような入力条件に対して、評価関数に基づいて最適化計算を行って最適解を求める最適解探索方法において、
前記制約条件は、解の探索において必須とされる始点条件と、解の探索において遵守した方が望ましいとされる終点条件とを有し、
前記入力条件を入力する入力ステップと、
前記始点条件を設定する始点条件設定工程と、前記終点条件を設定する終点条件設定工程と、これら設定された条件の組み合わせを踏まえて作成される探索パターンの順番を設定する探索パターン順番設定工程と、前記設定された探索パターンのうち、最適化計算において解を探索する際に、最初に使用する探索パターンを設定する初期探索パターン設定工程と、最適化計算において解を探索する際に、最初の探索する点を設定する初期探索点設定工程とで構成され、前記設定された条件を最適化計算を行うための初期の条件として設定する初期条件設定ステップと、
前記制約条件を満たす解を発見する解探索工程と、前記解探索工程で新たに発見された解を保存する解保存工程とを有し、解探索条件を満たすまで解の計算を行う最適化計算ステップと、
前記最適化計算ステップにおいて発見された解とこれまでに得られた最適解とを比較し、且つそれらの解の差を算出する解比較工程と、前記解比較工程で算出された解の差を基に、次回の解探索の際に使用する探索パターンを求める次回探索パターン選択工程と、前記次回探索パターン選択工程で求められた次回の探索パターンを設定する次回探索パターン設定工程と、前記次回探索パターン設定工程で設定された新たな探索パターンにおいて探索点を設定する新探索点設定工程とで構成され、前記設定した探索パターンを前記最適化計算ステップに戻す探索パターン設定ステップと、
前記解探索条件を満たした際に、前記解保存工程で保存された解の中から最適解を決定する最適解決定ステップと、
前記最適解決定ステップで決定された最適解を出力する出力ステップと、
を有することを特徴とする最適解探索方法。
【請求項２】
前記初期条件設定ステップは、前記始点条件と終点条件との間に、複数の中間点条件を設定する中間点条件設定工程を有することを特徴とする請求項１に記載の最適解探索方法。
【請求項３】
前記最適化計算ステップは、粒子群最適化方法に基づいた集団的降下法を用いるものとされ、
前記探索パターン設定ステップは、これまでの探索パターンから新たな探索パターンに変更して前記新探索点を設定する際に、これまでの最適解を解保存工程において保存し、この最適解とその周辺に解を探索する粒子を散布することを特徴とする請求項１に記載の最適解探索方法。
【請求項４】
前記探索パターン設定ステップは、前記制約条件の範囲を緩める緩和条件確認部が備えられており、
前記緩和条件確認部は、前記制約条件によって当該制約条件の範囲を緩めることが必要とされた場合、前記制約条件の範囲を緩める緩和条件を作成することを特徴とする請求項１又は３に記載の最適解探索方法。
【請求項５】
前記出力ステップは、前記最適化計算ステップで計算された解の結果、及び計算時の前記制約条件をオペレータに提示することを特徴とする請求項１〜４のいずれかに記載の最適解探索方法。
【請求項６】
請求項１〜５のいずれかに記載された最適解探索方法を、転炉の副原料の投入量の算出に用いることを特徴とする最適解探索方法。
【請求項７】
制約条件を満たすような入力条件に対して、評価関数に基づいて最適化計算を行って最適解を求める最適解探索装置において、
前記制約条件は、解の探索において必須とされる始点条件と、解の探索において遵守した方が望ましいとされる終点条件とが有するものであり、
前記入力条件を入力する入力部が備えられている入力手段と、
前記始点条件を設定する始点条件設定部と、前記終点条件を設定する終点条件設定部と、これら設定された条件の組み合わせを踏まえて作成される探索パターンの順番を設定する探索パターン順番設定部と、前記設定された探索パターンのうち、最適化計算において解を探索する際に、最初に使用する探索パターンを設定する初期探索パターン設定部と、最適化計算において解を探索する際に、最初の探索する点を設定する初期探索点設定部とが備えられ、前記設定された条件を最適化計算を行うための初期の条件として設定する初期条件設定手段と、
前記制約条件を満たす解を発見する解探索部と、前記解探索部で新たに発見された解を保存する解保存部とが備えられ、解探索条件を満たすまで解の計算を行う最適化計算手段と、
前記最適化計算手段において発見された解とこれまでに得られた最適解とを比較し、且つそれらの解の差を算出する解比較部と、前記解比較部で算出された解の差を基に、次回の解探索の際に使用する探索パターンを求める次回探索パターン選択部と、前記次回探索パターン選択部で求められた次回の探索パターンを設定する次回探索パターン設定部と、前記次回探索パターン設定部で設定された新たな探索パターンにおいて探索点を設定する新探索点設定部とが備えられ、前記設定した探索パターンを前記最適化計算手段に戻す探索パターン設定手段と、
前記解探索条件を満たした際に、前記解保存部で保存された解の中から最適解を決定する最適解決定部が備えられている最適解決定手段と、
前記最適解決定手段で決定された最適解を出力する出力部が備えられている出力手段と、
を有することを特徴とする最適解探索装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ある条件下において予測モデルを用いて、転炉に投入される副原料の投入量等の解を計算する際に、解を評価する評価関数から最も良い解を短時間に予測する最適解探索方法、及び最適解探索装置に関する。
【背景技術】
【０００２】
従来から、転炉で行われる溶銑の精錬では、溶鋼中の炭素（Ｃ）、りん（Ｐ）およびマンガン（Ｍｎ）等の各濃度に対し目標値が設定されており、酸素（気酸、固酸）の供給量や副原料の投入量を適切にすることで、この目標値を満たすように各成分の調整が行われる。副原料の投入は、溶鋼の品質を考慮しつつ、副原料のコストを考慮することによって行われるが、その投入量の決定は、長年の操業実績から得られた手法が用いられている。
【０００３】
その一方で、様々な数学的手法や予測モデルを適用して、副原料の投入量の決定を行う技術も開発されている。
例えば、特許文献１は、溶銑に投入する副原料の量の最適解を求める技術に関するものであり、制約条件を満たす解を、評価関数を使用して発見する評価手段と、制約条件として第１条件と、第２条件とを記憶しており、前記第１条件を制約条件として評価関数を使用して解を発見させ、発見した解に基づいて新たな条件を生成し、生成した条件を制約条件として解を発見させることを所定の終了条件を満たすまで繰り返し実行し、発見した解のうち、前記第２条件に最も近い解を最適解として決定する最適解を探索する方法を開示する。
【０００４】
また、転炉に関係する技術ではないものの、ファジー等式制約や有限幅等式制約と確率密度を用いてある制約条件が無い場合の最適解が制約条件の外にある場合でも制約条件内で最適な解を探索する最適解探索方法に関する技術が特許文献２に開示されている。
すなわち、特許文献２は、不等式制約によって探索空間内に生成される閉集合の境界線を等式制約として分離して、ファジー等式制約もしくは有限幅等式制約とし、そこにメンバーシップ関数によって関係づけられる確率密度もしくは一定の密度で探索点を配置し、別途等確率密度で探索点が配置される不等式部分と併せて、境界線内側近傍もしくは境界線近傍部分の探索点の配置密度を上げることによって最適点の発見確率を上げる手段を最適化探索の各個体部に持たせる最適化装置を開示する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１３−４１３２１号公報
【特許文献２】特開２０１０−２５０３９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
前述した如く、転炉で行われる溶銑の精錬では、様々な副原料が吹練条件に応じて投入されることとなるが、投入量の決定は非常に複雑で難しいものとなっている。
そこで、特許文献１のような予測モデルを用いて最適解を探索することで、副原料の投入量を決定することが考えられる。その場合、計算に用いる制約条件の選定が問題となる。例えば、制約溶銑の品質の調整が容易ではない工程や、データに大きなノイズ要素が存在する工程の最適化計算をするとき、溶銑の品質をできる限り維持するために、制約条件に始点条件（必ず実現しなくてはいけない条件）と終点条件（出来れば実現したい条件）を加えることがある。
【０００７】
しかしながら、特許文献１の技術を用いて、終点条件内に存在する解、若しくは終点条件に近い解を求める場合、始点条件から最適化計算し始めると、評価関数によっては始点条件の境界線近傍に解が集中してしまう虞がある。このような問題を解決する場合、終点条件から最適化計算を始め、終点条件に解が存在しない場合には制約条件を逐次緩和させて再計算する方法があるが、探索回数が増加してしまうため、最適解を探索する時間が長くなる。また、特許文献１は、過去の解情報をリセットせずに次の解探索を行っている。
この過去の解情報は最適化計算を行う上で重要な情報となり得る。
【０００８】
一方で、最適解を探索するにあたり、粒子群最適化（ＰＳＯ：Particle Swarm Optimization）のような集団的降下法を使用することが考えられる。この場合、制約条件が変化する状況下においては、リセットされない過去の解情報が逆に解の探索を妨げる虞がある。なぜならば、制約条件を限定したり、制約条件を緩和したりする際に、前段階の範囲を探索していた粒子は、限定された条件部分や追加された条件部分へ辿り着くまでに、無駄な解の探索（解の存在する可能性が低い範囲での探索）を行うためである。
【０００９】
また、特許文献１は、制約条件の範囲を限定するのみであるが、範囲を限定するだけでは、状況の表現が的確に表現できないことがある。例えば、制約条件が複数存在し、且つトレードオフ等の関係によって制約条件間に明確な優先順位が存在する場合には、範囲を限定するだけではなく、範囲を緩和する方が状況を的確に表現できる。
ゆえに、実操業下においては、短時間で、且つ効率的に最適解を探索することが求められているため、特許文献１の技術を改良する要望が挙げられている。
【００１０】
そこで、特許文献１の技術に、特許文献２の最適化技術を適用することが考えられる。特許文献２は、見つけにくい解をより、短時間で見つけることができる点（解）と、制約条件が無い場合の最適点（最適解）とが制約条件の範囲外にある場合を想定して最適化を行っており、確率密度によって探索点を絞り、境界線内側近傍もしくは境界線近傍部分の探索点の配置密度を上げることによって最適点の発見確率を上げている。
【００１１】
しかしながら、特許文献２の技術は、制約条件を１つとして考えているが、実操業では一つの変数に、可能な限り実現したい条件（終点条件）と、必ず実現しなくてはいけない条件（始点条件）とが存在することが多くある。そのため、特許文献２の技術では、制約条件外に最適点があることから評価関数も境界線近傍に行くことが多く、終点条件内に解を向かわせることは困難である。また、実操業下での解の探索を考えた場合、始点条件が操業条件の限界点であるとすると、ノイズの大きさによっては始点条件近傍に解が集中することを好まないことがある。ゆえに、特許文献２に記載された技術をそのまま転炉の副原料投入量決定の最適計算を適用することはできない。
【００１２】
そこで、本発明は上記問題点を鑑み、制約条件をもとに解を探索する順番を設定し、その後、制約条件の終点条件内、若しくは終点条件近傍で解を探索して、これら解の中から最適解を決定する最適解探索方法、及び最適解探索装置を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上述の目的を達成するため、本発明においては以下の技術的手段を講じた。
本発明に係る最適解探索方法は、制約条件を満たすような入力条件に対して、評価関数に基づいて最適化計算を行って最適解を求める最適解探索方法において、前記制約条件は、解の探索において必須とされる始点条件と、解の探索において遵守した方が望ましいとされる終点条件とを有し、前記入力条件を入力する入力ステップと、前記始点条件を設定する始点条件設定工程と、前記終点条件を設定する終点条件設定工程と、これら設定された条件の組み合わせを踏まえて作成される探索パターンの順番を設定する探索パターン順番設定工程と、前記設定された探索パターンのうち、最適化計算において解を探索する際に、最初に使用する探索パターンを設定する初期探索パターン設定工程と、最適化計算において解を探索する際に、最初の探索する点を設定する初期探索点設定工程とで構成され、前記設定された条件を最適化計算を行うための初期の条件として設定する初期条件設定ステップと、前記制約条件を満たす解を発見する解探索工程と、前記解探索工程で新たに発見された解を保存する解保存工程とを有し、解探索条件を満たすまで解の計算を行う最適化計算ステップと、前記最適化計算ステップにおいて発見された解とこれまでに得られた最適解とを比較し、且つそれらの解の差を算出する解比較工程と、前記解比較工程で算出された解の差を基に、次回の解探索の際に使用する探索パターンを求める次回探索パターン選択工程と、前記次回探索パターン選択工程で求められた次回の探索パターンを設定する次回探索パターン設定工程と、前記次回探索パターン設定工程で設定された新たな探索パターンにおいて探索点を設定する新探索点設定工程とで構成され、前記設定した探索パターンを前記最適化計算ステップに戻す探索パターン設定ステップと、前記解探索条件
を満たした際に、前記解保存工程で保存された解の中から最適解を決定する最適解決定ステップと、前記最適解決定ステップで決定された最適解を出力する出力ステップと、を有することを特徴とする。
【００１４】
好ましくは、前記初期条件設定ステップは、前記始点条件と終点条件との間に、複数の中間点条件を設定する中間点条件設定工程を有するとよい。
好ましくは、前記最適化計算ステップは、粒子群最適化方法に基づいた集団的降下法を用いるものとされ、前記探索パターン設定ステップは、これまでの探索パターンから新たな探索パターンに変更して前記新探索点を設定する際に、これまでの最適解を解保存工程において保存し、この最適解とその周辺に解を探索する粒子を散布するとよい。
【００１５】
好ましくは、前記探索パターン設定ステップは、前記制約条件の範囲を緩める緩和条件確認部が備えられており、前記緩和条件確認部は、前記制約条件によって当該制約条件の範囲を緩めることが必要とされた場合、前記制約条件の範囲を緩める緩和条件を作成するとよい。
好ましくは、前記出力ステップは、前記最適化計算ステップで計算された解の結果、及び計算時の前記制約条件をオペレータに提示するとよい。
【００１６】
本発明に係る最適解探索方法は、転炉の副原料の投入量の算出に用いるとよい。
本発明に係る最適解探索装置は、制約条件を満たすような入力条件に対して、評価関数に基づいて最適化計算を行って最適解を求める最適解探索装置において、前記制約条件は、解の探索において必須とされる始点条件と、解の探索において遵守した方が望ましいとされる終点条件とが有するものであり、前記入力条件を入力する入力部が備えられている入力手段と、前記始点条件を設定する始点条件設定部と、前記終点条件を設定する終点条件設定部と、これら設定された条件の組み合わせを踏まえて作成される探索パターンの順番を設定する探索パターン順番設定部と、前記設定された探索パターンのうち、最適化計算において解を探索する際に、最初に使用する探索パターンを設定する初期探索パターン設定部と、最適化計算において解を探索する際に、最初の探索する点を設定する初期探索点設定部とが備えられ、前記設定された条件を最適化計算を行うための初期の条件として設定する初期条件設定手段と、前記制約条件を満たす解を発見する解探索部と、前記解探索部で新たに発見された解を保存する解保存部とが備えられ、解探索条件を満たすまで解の計算を行う最適化計算手段と、前記最適化計算手段において発見された解とこれまでに得られた最適解とを比較し、且つそれらの解の差を算出する解比較部と、前記解比較部で算出された解の差を基に、次回の解探索の際に使用する探索パターンを求める次回探索パターン選択部と、前記次回探索パターン選択部で求められた次回の探索パターンを設定する次回探索パターン設定部と、前記次回探索パターン設定部で設定された新たな探索パターンにおいて探索点を設定する新探索点設定部とが備えられ、前記設定した探索パターンを前記最適化計算手段に戻す探索パターン設定手段と、前記解探索条件を満たした際に、前記解保存部で保存された解の中から最適解を決定する最適解決定部が備えられている最適解決定手段と、前記最適解決定手段で決定された最適解を出力する出力部が備えられている出力手段と、を有することを特徴とする。
【発明の効果】
【００１７】
本発明の最適解探索方法、及び最適解探索装置よれば、制約条件をもとに解を探索する順番を設定し、その後、制約条件の終点条件内、若しくは終点条件近傍で解を探索して、これら解の中から最適解を決定することになるため、最も良い解を短時間に算出することが可能となる。
【図面の簡単な説明】
【００１８】
【図１】本発明の最適解探索装置のシステム構成を示した図である。
【図２】粒子が解を探索する際に設定される始点条件及び終点条件を模式的に示した図である。
【図３】粒子が解を探索する探索パターンの種類を示した図である。
【図４】本発明の最適解探索方法を示したフローチャート図である。
【図５】最適解探索方法の変形例を示したフローチャート図である。
【図６】転炉の概略を示した図である。
【発明を実施するための形態】
【００１９】
以下、本発明に係る最適解探索方法、および最適解探索装置１の実施の形態を、図をもとに説明する。
本発明に係る最適解探索方法、および最適解探索装置１の詳細を述べる前に、本発明の技術が適用される転炉３０について、図をもとに説明する。
図６に示す如く、脱りん処理を行う転炉３０は、気体酸素を溶銑に吹き込む上吹きランス３１と炉底から酸素又は不活性ガスを溶銑に吹き込む羽口３３を備えた上底吹き型であって、上吹きランス３１からの気体酸素により酸素を供給し、羽口３３からの酸素又は不活性ガスにより溶湯を攪拌するものである。また、転炉３０は、供給装置３２を備えている。この供給装置３２は、副原料（生石灰、固体酸素源等）を供給するものであって、例えば、ホッパーやシュート等である。
【００２０】
図１は、前述した転炉３０において、投入する副原料の量を正確、且つ短時間で算出するための最適解を探索する装置のシステム構成を示したものである。なお、最適解を探索するにあたっては、本実施形態では、粒子群最適化法（Particle Swarm Optimization、ＰＳＯ法）等の集団的降下法を用いている。
最適解探索装置１は、転炉３０における入力条件を入力する入力手段２と、最適化計算を行うための初期の条件を設定する初期条件設定手段３と、解探索条件を満たすまで解の計算を行う最適化計算手段１０と、最適化計算手段１０での計算に用いる探索パターンを設定し、設定した探索パターンを前記最適化計算ステップに送る探索パターン設定手段１３と、解探索条件を満たした際に、後述の解保存部１２で保存された解の中から最適解を決定する最適解決定手段１８と、最適解決定手段１８で決定された最適解を出力する出力手段２０と、を備えている。
【００２１】
まず、入力手段２では、最適化計算を行うため、転炉３０における初期の条件（制約条件を満たすような入力条件）を入力部に入力する。本実施形態の入力条件は、転炉３０に投入する副原料の量を決定するために必要な溶銑の成分（溶銑燐比率、溶銑温度、溶銑配合率など）、目標値（溶鋼燐比率、溶鋼温度など）及び、それ以外の情報（炉令、休炉時間など）を入力する。
【００２２】
初期条件設定手段３では、最適化計算を行うための初期の条件、すなわち制約条件の始点条件（解の探索において遵守することが必須とされる条件）及び終点条件（解の探索において遵守した方が望ましいとされる条件）の設定と、これら設定された条件（始点条件及び終点条件）の組み合わせを踏まえて作成される探索パターンの順番の設定と、初期の探索パターンに選ばれた条件の設定、及び初期の探索点の設定が行われる。
【００２３】
また、初期条件設定手段３では、始点条件と終点条件の間に位置し、且つ可能であれば遵守した方が望ましいとされる中間点条件の設定も必要な場合に行われる。この中間点条件を設定することで、より多くの探索パターンが設定でき、制約条件間の関係を精度良く表した最適化計算が可能となる。
本実施形態の初期条件設定手段３では、探索パターンの探索順番の設定は、始点条件と終点条件に中間点条件を加えた３パターンの組み合わせを踏まえて作成されている。
【００２４】
図１に示すように、この初期条件設定手段３は、始点条件を設定する始点条件設定部４と、終点条件を設定する終点条件設定部６と、これら設定された条件の組み合わせを踏まえて作成される探索パターンの順番を設定する探索パターン順番設定部７と、設定された探索パターンのうち、最適化計算において解を探索する際に、最初に使用する探索パターンを設定する初期探索パターン設定部８と、最適化計算において解を探索する際に、最初の探索する点を設定する初期探索点設定部９と、を備えている。
【００２５】
また、本実施形態の初期条件設定手段３は、始点条件と終点条件との間に、複数の中間点条件を設定する中間点条件設定部５を有する。
転炉３０における初期の制約条件を設定するにあたっては、以下の方法を用いる。
一般的な最適化問題（Ｐ） は、不等式制約、等式制約、上下限制約を有しており、以下の式（１）のように定義できる。
【００２６】
【数１】
2014209273.tif 000003
【００２７】
ここで、ｘ＝（ｘ１，・・・，ｘｎ）はｎ次元決定変数ベクトル、ｆ（ｘ）は評価関数、ｇｊ（ｘ）≦α  はｑ個の不等式制約、ｈｊ（ｘ）＝βはｍ−ｑ個の等式制約であり、ｆ，ｇｊ，ｈｊは線形あるいは非線形の実数値関数である。ｌｉ，ｕｉはそれぞれｎ個の決定変数ｘｉの下限値、上限値である。
本実施形態の評価関数ｆ（ｘ）は、コストを算出する関数であり、ｘは、副原料の投入量を示すことになる。また、例えば、ｇ（ｘ）が溶鋼燐比率を求める関数であり、ｈ（ｘ）が溶鋼温度を求める関数となる。
【００２８】
また、本実施形態で設定する始点条件、終点条件、及び中間点条件については、不等式制約、等式制約、及び決定変数の上下限値である。
ここで、始点条件、終点条件、及び中間点条件と、この３つの条件を踏まえて組合せられる探索パターンについて、説明する。
図２は、始点条件と終点条件についての模式図である。始点条件は、副原料Ａ及び副原料Ｂの限定（規制）が緩和された条件であり、終点条件よりも粒子の探索範囲が広くなっている。一方で、終点条件は、副原料Ａ及び副原料Ｂの限定を厳しくして、粒子の探索範囲を狭めている。
【００２９】
なお、中間点条件は、図示していないが、始点条件と終点条件との間に存在し、副原料Ａ及び副原料Ｂの限定を終点条件より若干緩和して設定される。これら３つの条件は、それぞれの操業条件や操業における目標値によって決定される。制約条件を満たす解を少なくとも１つ以上発見したい場合には、始点条件を許容限界まで緩和すると望ましい。
図３に示すように、探索パターンは、上述した始点条件、中間点条件、及び終点条件を用いて、制約条件のパターン分けを行ったものである。
【００３０】
本実施形態では、転炉３０の操業においてスラグがスロッピングすることを抑えるための合計投入量制約（副原料Ｂ制約）と、転炉３０の炉内壁を保護すべく、スラグ中に存在する副原料Ａの成分比率を一定以上に守るための副原料Ａ制約とを用いている。
合計投入量制約（副原料Ｂ制約）については、投入できる限界量を始点条件に設定し、できるだけ守りたい合計投入量を終点条件に設定している。なお、副原料Ｂ制約の中間点条件の設定は行っていない。
【００３１】
副原料Ａ制約については、スラグ中の副原料Ａの成分比率が許容される限界点を始点条件に設定し、スラグ中の副原料Ａの成分比率が十分に満足する点を終点条件に設定している。なお、副原料Ａ制約の中間点条件は、始点条件と終点条件の間に設定している。
制約条件は他にも存在するが、条件の規制を緩和することが難しい制約条件については、あえて始点条件、及び終点条件を設定しなくてもよい。このような場合の制約条件は、終点条件のみが存在する、もしくは始点条件と終点条件が同じ条件と見なすことができる。
【００３２】
表１に示すように、本実施形態では、探索パターンを６つに設定している。
図３（ｆ）に示すように、パターン６は、２つの制約条件で始点条件となっており、パターン１〜パターン５の範囲を含むものとなっている。すなわち、パターン６は、粒子の探索範囲が６つの探索パターンのうち、最も緩和されたものである。図３（ｂ）〜図３（ｅ）に示すように、パターン番号が小さくなるにつれて、粒子の探索範囲は小さくなり、制約条件の規制が厳しくなる。
【００３３】
図３（ａ）に示すように、パターン１は、最も厳しい規制であり、２つの制約条件で終点条件となる。
【００３４】
【表１】
2014209273.tif 000004
【００３５】
これら６個のパターンに基づいて、探索パターンの順番（粒子が解を探索する範囲の順番）を設定する。本実施形態においては、表２に示す探索の順番で解の探索を行う。なお、探索パターンの順番は、表２に示す順番だけでなく、転炉３０の操業下において任意に設定することができる。
【００３６】
【表２】
2014209273.tif 000005
【００３７】
本実施形態の探索順は、最初に解を探索する「探索順１」に、探索領域が最も広いパターン６を配置し、最後に解を探索する「探索順６」にパターン１を配置する。「探索順１」と「探索順６」との間には、順次制約条件が縮小するように探索パターンを配置する。
このとき、探索パターンの順番の設定にあたっては、副原料Ａ制約と合計投入量制約（副原料Ｂ制約）とが干渉し合う関係にならないようにする。
【００３８】
なお、制約条件が干渉し合う関係のときに探索パターンの順番を設定する場合については、後述する変形例で説明する。
次に、最初に探索する範囲として設定されたパターン６内にそれぞれの粒子を割り振り、割り振られた地点を最適化計算を行うための粒子の初期探索地点と設定する。なお、ＰＳＯ法における粒子の探索範囲は、通常、各ＰＳＯ変数で設定される範囲である。
【００３９】
最適化計算手段１０では、制約条件を満たすまで、ＰＳＯ法など最適化手法で解の計算を繰り返し行われ、制約条件を満たす解を発見した際には解保存部１２で保存する。
図１に示すように、最適化計算手段１０は、制約条件を満たす解を発見する解探索部１１と、解探索部１１で新たに発見された解を保存する解保存部１２とを備えている。
解探索部１１では、転炉３０における制約条件を満たす解を発見されるまで探索する。そして、解探索部１１で発見された解は、解保存部１２に保存される。解保存部１２には、複数個、解が保存されていて、それら複数の解は最適解の候補とされる。
【００４０】
ここで、解探索部１１で解の探索を行い、解探索部１１で発見された解を解保存部１２に保存するまで過程を、詳細に述べる。
各変数の制約条件を始点条件として、ＰＳＯ法を用いて最適化計算を開始する。
ＰＳＯ法は集団的降下法の一つであり、集団的降下法での探索ルーチンは一般的に以下のように記述できる。
【００４１】
（１）初期化: 集団に属する解をランダムに発生させる。
（２）評価: 発生した全ての解を評価する。
（３）終了判定: 終了条件を満足すると、探索を終了する。
（４）各解に対して、
（ａ） 生成: 各解と集団の情報に基づき、新しい解を生成する。
【００４２】
（ｂ） 評価: 生成された新しい解を評価する。
（ｃ） 更新: 新しい解が古い解より良ければ、最適解を新しい解に置き換える。
（５）終了判定へ戻る。
なお、（２）の「評価」では、発生した全ての解が評価関数の値と制約条件とを満たしているか否かについて、評価している。（２）の「評価」で制約条件を満たす解が発見されなかった場合、（４）（ａ）の「生成」に移り、制約条件を満たす解が発見されるまで探索を繰り返す。また、（１）「初期化」については、上述した粒子の初期探索地点の設定に該当する。
【００４３】
探索パターン設定手段１３は、過去の解の探索結果をリセットする際に、解情報を全てリセットしてしまうと、これまでの最適解まで失ってしまい、解が最後まで発見できない可能性や効率的な解探索ができない問題を回避するものであり、過去の解の探索結果をリセットする際にも、最適解（最適点）の情報を記憶するようになっている。更に、その最適解の周辺に粒子を重点的に配置することで効率的な探索を可能にしている。
【００４４】
探索パターン設定手段１３では、これまでの探索パターンから新たな探索パターンに変更して新探索点を設定する際に、これまでの最適解を解保存部１２において保存し、この最適解とその周辺に解を探索する粒子を散布する。
図１に示すように、探索パターン設定手段１３は、最適化計算手段１０において発見された解とこれまでに得られた最適解とを比較し、且つそれらの解の差を算出する解比較部１４と、解比較部１４で算出された解の差を基に、次回の解探索の際に使用する探索パターンを求める次回探索パターン選択部１５と、次回探索パターン選択部１５で求められた次回の探索パターンを設定する次回探索パターン設定部１６と、次回探索パターン設定部１６で設定された新たな探索パターンにおいて探索点を設定する新探索点設定部１７と、を備えている。
【００４５】
ここで、探索パターン設定手段１３において、次回の探索パターンを設定し、設定された探索パターンの中で新たな探索点を設定する過程を詳細に説明する。
まず、最適化計算手段１０で発見した解が下位の探索順のパターン（例えば、表１のパターン１）まで到達しているかを判断する。このとき、その解が下位の探索順のパターンに到達したと判断された場合、そのパターンまで探索順を進める。なお、その解が下位の探索順のパターンに到達しておらず、且つ現状のパターンのままの場合、最適化計算を終えるための終了条件の判定を行う。
【００４６】
終了条件を満たしていないと判断された場合、現状のパターンで解の探索を継続する。なお、終了条件を満たしていると判断された場合、最適解決定手段１８へ移り、最適解の決定を行う。
なお、本実施形態の終了条件については、「指定回数の探索を行ったかどうか」と「一定回数以内に最適解が更新されているかどうか」を終了条件としている。
【００４７】
表３に示すように、解が下位の探索順のパターンまで到達していると判断された場合、到達したパターンの制約条件を変更する。具体的には、発見された最適解のパターンより緩い条件（始点条件に近い条件）のパターン、及び探索していたパターンは、探索済みとする。ゆえに、より緩い条件のパターンにおける解の探索が実施されないようになっている。
【００４８】
【表３】
2014209273.tif 000006
【００４９】
表４に示すように、探索中のパターンで解が発見されなかった場合、それより厳しい条件のパターンでは解が存在しないと判断され、この厳しい条件のパターン、及び探索中のパターンを探索済みとする。ゆえに、より厳しい条件のパターンにおける解の探索が実施されないようになっている。
【００５０】
【表４】
2014209273.tif 000007
【００５１】
到達したパターンの制約条件を変更した後、最適解以外の解情報を消去し、発見された暫定の最適解を１つの粒子の探索座標に再設定し、その最適解の周辺の探索点に指定された数の粒子を探索座標に設定する。また、残りの粒子の探索点を該当探索順番のパターン内でランダムに再設定する。
このようにすることで、範囲縮小によってパターンを変更した後に一度見つけた解が見つからないといった問題を無くし、効率的な最適解探索が行える。なお、暫定最適解の近傍に配置されなかった残りの粒子の探索点については、設定されたパターンの制約条件の範囲内で再度ランダムに再設定される。
【００５２】
最適解の近傍の上限値と下限値を算出するにあたっては、式（２）、及び式（３）を用いる。
最適解近傍の上限値 ＝ ｍｉｎ（探索座標上限 ，最適解探索座標＋（探索座標上限−
探索座標下限）×近傍範囲（％）／２)  ・・・(２)  
最適解近傍の下限値 ＝ ｍａｘ（探索座標下限 ，最適解探索座標−（探索座標上限−
探索座標下限）×近傍範囲（％）／２)  ・・・(３)
粒子の探索点が再設定された後、発見された解が最下位の探索順のパターンに到達しているか否かを判断する。探索終了条件まで最適化計算を行う。
【００５３】
発見された解が最下位の探索順のパターンに（例えば、表１の「探索順６」のパターン１）まで到達している場合、最終の探索順であるため、パターンは変更されない。この場合、探索を終了する探索終了条件まで最適化計算を行う。
最適解決定手段１８では、解探索条件を満たした際に、最適解決定部１９において、解保存部１２で保存された解の中から最も評価関数の良い解を副原料の投入量の最適解とし
、最適化計算を終了する。
【００５４】
出力手段２０では、最適解決定手段１８で決定された副原料の投入量の最適解及び当該最適解に関する制約条件を出力部から出力する。
副原料の投入量の最適解を出力するにあたっては、その最適解を転炉３０の設備の近傍に配備されたモニタなどの表示器に示し、転炉３０の操業に携わるオペレータに通達する。オペレータに最適解を通達することによって、転炉３０の操業の安定化や操業効率の向上を図ることができる。
【００５５】
以下、本発明の最適解探索装置１を用いて最適解を探索する方法について、図をもとに詳しく説明する。
図４は、最適解探索装置１を用いて最適解を探索する方法を示したフローチャートである。
図４に示すように、ステップＳ１−１では、転炉３０に投入する副原料の量を決定するために必要な溶銑の成分（溶銑燐比率、溶銑温度、溶銑配合率など）、目標値（溶鋼燐比率、溶鋼温度など）及び、それ以外の情報（炉令、休炉時間など）を入力する（入力ステップ）。
【００５６】
ステップＳ１−２では、始点条件設定工程、終点条件設定工程及び中間点条件設定工程で、最適化計算を行うための初期の条件、すなわち制約条件の始点条件、終点条件及び中間点条件の設定が行われ、これら設定された条件（始点条件、終点条件及び中間点条件）の組み合わせを踏まえて作成される探索パターンの順番の設定が探索パターン順番設定工程で行われる（初期条件設定ステップ）。
【００５７】
ステップＳ１−３では、初期探索パターン設定工程で、ステップＳ１−２で設定した探索順番の１番目にあたるパターンの制約条件を設定する（初期条件設定ステップ）。
ステップＳ１−４では、初期探索点設定工程において、最適化計算を行うために、ステップＳ１−３で設定されたパターン内にそれぞれの粒子を割り振り、割り振られた地点を粒子の初期探索地点と設定する（初期条件設定ステップ）。
【００５８】
ステップＳ１−５及びステップＳ１−６では、解探索工程でＰＳＯ法を用いて転炉３０における制約条件を満たす解を探索する（最適化計算ステップ）。
ステップＳ１−７では、では、解探索工程で発見された解を解保存部１２に保存する（解保存工程、最適化計算ステップ）。
ステップＳ１−８では、解比較工程で、探索によって発見された解が下位の探索順のパターンまで到達しているかを判断する（探索パターン設定ステップ）。解が下位の探索順のパターンに到達している場合、ステップＳ１−１１に移行する（ｙｅｓ）。また、解が下位の探索順のパターンに到達しておらず、現状のパターンのままの場合、ステップＳ１−９に移行する（ｎｏ）。
【００５９】
ステップＳ１−９では、最適化計算を終えるための探索終了条件を満たしているかどうかを判断する（探索パターン設定ステップ）。探索終了条件を満たしていない場合、ステップＳ１−５（解探索工程）に戻り、現状のパターンで解の探索を継続する（ｎｏ）。探索終了条件を満たしている場合、ステップＳ１−１０に移行する（ｙｅｓ）。
ステップＳ１−１０では、終了条件を満たした後、現状の探索パターンの中でもっとも評価関数の良かった解を最適解として決定する（最適解決定ステップ）。
【００６０】
ステップＳ１−１１では、ステップＳ１−８で解が下位の探索順パターンまで到達していることが分かった場合、ここで到達したパターンに制約条件を次回探索パターン選択工程で変更する（探索パターン設定ステップ）。
ステップＳ１−１２では、次回探索パターン設定工程で最適解以外の解情報を消去する（探索パターン設定ステップ）。
【００６１】
ステップＳ１−１３では、新探索点設定工程で発見された暫定の最適解を１つの粒子の探索座標に設定し、指定された粒子数だけその近傍にランダムに探索座標を設定する（探索パターン設定ステップ）。
ステップＳ１−１４では、残りの粒子の探索点を該当探索順番のパターン内で再設定する（探索パターン設定ステップ）。
【００６２】
ステップＳ１−１５では、解が最下位の探索順のパターンまで到達しているか否かを判断する（探索パターン設定ステップ）。解が最下位の探索順のパターンまで到達している場合、ステップＳ１−１６に移行する（ｙｅｓ）。解が最下位の探索順のパターンまで到達していない場合、ステップＳ１−５（解探索工程）に戻り、現状のパターンで解の探索を継続する（ｎｏ）。
【００６３】
ステップＳ１−１６では、終了条件を満たすまで、最適化計算を繰り返す（最適化計算ステップ）。
ステップＳ１−１７では、終点条件を満たした際に、最適解決定工程において、解保存部１２で保存された解の中から最も評価関数の良い解を副原料の投入量の最適解とする（最適解決定ステップ）。
【００６４】
ステップＳ１−１８では、最適解決定ステップで決定された副原料の投入量の最適解及び当該最適解に関する制約条件を転炉３０の設備の近傍に配備されたモニタなどの表示器に出力する（出力ステップ）。
以上述べた最適解探索方法、及び最適解探索装置１を用いることで、制約条件をもとに解を探索する順番を設定し、その後、制約条件の終点条件内、若しくは終点条件近傍で解を探索して、これら解の中から最適解を決定することになるため、最も良い解を短時間に算出することが可能となる。
【００６５】
[変形例]
ところで、上述した最適解探索装置１、及び最適解探索方法は、始点条件、及び終点条件を有する２つの制約条件が干渉し合う関係ではない場合、始点条件と終点条件の組合せ（制約条件）を基に作成される探索パターンの順番を設定し、そして最初の探索パターンで解を探索した後、その解を基に次回の探索パターンを再設定して再び解を探索し、これら解の中から最適解を決定する技術であることを述べた。
【００６６】
しかしながら、上述の最適解探索装置１、及び最適解探索方法は、２つ制約条件が干渉し合う場合についても採用することができる。
この場合、探索パターン設定手段１３に、制約条件の範囲を緩める緩和条件確認部を備えておくとよい。この緩和条件確認部は、制約条件によって当該制約条件の範囲を緩めることが必要とされた場合、前記制約条件の範囲を緩める緩和条件を作成する。
【００６７】
具体的には、制約条件間で干渉し合う関係の場合、それら制約条件の間に優先順位がつけられることがある。最適解を探索する際に、制約条件を緩和する緩和ループ（緩和条件確認部）を加えて、探索パターンの探索順を変更することで、制約条件間で干渉し合う関係の場合においても本発明の最適解探索装置１、及び最適解探索方法を採用することが可能となる。
【００６８】
例えば、制約条件Ａと制約条件Ｂが干渉する関係にあり、且つ制約条件Ａの重要度が高い場合（優先権がある場合）、探索パターンの順番は表６に示すように設定する。なお、変形例の探索パターンについては、表５に示すものを採用しているが、上述した始点条件と終点条件の組合せを基に作成される探索パターンであれば、いずれのものを採用することもできる。
【００６９】
【表５】
2014209273.tif 000008
【００７０】
【表６】
2014209273.tif 000009
【００７１】
以下、上述した探索パターン、及び探索パターンの順番で最適化計算を行う場合について、図５を基に説明する。
なお、図５に示されているステップＳ１−１〜ステップＳ１−５、ステップＳ１−７〜ステップＳ１−１８については、図４に示すステップステップＳ１−１〜ステップＳ１−５、ステップＳ１−７〜ステップＳ１−１８と同様であるので、説明は省略する。
【００７２】
つまり、変形例の特徴である緩和条件確認部に該当する、ステップＳ２−６、ステップＳ２−１９〜ステップＳ２−２２に着目して、説明する。
まず、図５に示すステップＳ１−１〜ステップＳ１−５については、上述したように進めてゆく。
ステップＳ２−６では、ステップＳ１−５で解を探索した結果、解が発見されているか否かを判断する（解探索工程）。解が発見されている場合、表６に示す探索順１，５，６の探索パターン（パターン３，パターン２，パターン１）の中から最適解を見つける（最適化計算ステップ）。このとき、探索順２〜４のパターン（パターン６，パターン５，パターン４）は、優先順の関係から自動的に「探索済」となる。そして、ステップＳ１−７に移行する（ｙｅｓ）。また、解が発見されていない場合、ステップＳ２−１９に移行する（ｎｏ）。
【００７３】
ステップＳ２−１９では、表６に示す探索順１で解が見つからない場合、探索順２〜４に進むための緩和条件を満たしているか否かを判断する（最適化計算ステップ）。例えば、一定回数以上の探索を行っても解が発見されない場合などが挙げられる。緩和条件を満たす場合、ステップＳ２−２０に移行する（ｙｅｓ）。緩和条件を満たさない場合、ステップＳ１−５（解探索工程）に戻って解探索を再度行う（ｎｏ）。
【００７４】
ステップＳ２−２０では、探索順２のパターン６に変更する（最適化計算ステップ）。
ステップＳ２−２１では、粒子を再度散布するために解情報を消去する（最適化計算ステップ）。
ステップＳ２−２２では、パターン６内に粒子を再散布し、その後、ステップＳ１−５に戻って解を探索する（最適化計算ステップ）。
【００７５】
このように、図５に示す制約条件を緩和する緩和ループ（ステップＳ２−６，ステップＳ２−１９〜ステップＳ２−２２）を加えて、探索パターンの探索順を変更することで、複数の段階に分けて解を探索することで、より良い最適解を決定することができる。
そして、図５に示すステップＳ１−７〜ステップＳ１−１８については、上述した実施形態のように進める。
【００７６】
上記した変形例に係る最適解探索方法、及び最適解探索装置１であっても、制約条件をもとに解を探索する順番を設定し、その後、制約条件の終点条件内、若しくは終点条件近傍で解を探索して、これら解の中から最適解を確実に決定することが可能となる。
ところで、本発明は上記各実施形態に限定されるものではなく、発明の本質を変更しない範囲で各部材の形状、構造、材質、組み合わせなどを適宜変更可能である。また、今回開示された実施形態において、明示的に開示されていない事項、例えば、運転条件や操業条件、各種パラメータ、構成物の寸法、重量、体積などは、当業者が通常実施する範囲を逸脱するものではなく、通常の当業者であれば、容易に想定することが可能な事項を採用
している。
【符号の説明】
【００７７】
１    最適解探索装置
２    入力手段（入力部）
３    初期条件設定手段
４    始点条件設定部
５    中間点条件設定部
６    終点条件設定部
７    探索パターン順番設定部
８    初期探索パターン設定部
９    初期探索点設定部
１０  最適化計算手段
１１  解探索部
１２  解保存部
１３  探索パターン設定手段
１４  解比較部
１５  次回探索パターン選択部
１６  次回探索パターン設定部
１７  新探索点設定部
１８  最適解決定手段
１９  最適解決定部
２０  出力手段（出力部）
３０  転炉
３１  上吹きランス
３２  供給装置
３３  羽口
【図１】
2014209273.tif 000010
【図３】
2014209273.tif 000011
【図４】
2014209273.tif 000012
【図５】
2014209273.tif 000013
【図６】
2014209273.tif 000014
【図２】
2014209273.tif 000015