"燃焼ガスにより蓄熱煉瓦を加熱して蓄熱する燃焼期間と、当該蓄熱煉瓦に冷風を通して当該蓄熱煉瓦との熱交換により熱風を生成して高炉に供給する送風期間とを含む期間を1サイクルとしてそれぞれが稼働する複数の熱風炉の操業の状態を減風減温状態にした後、当該複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための炉別サイクル別の投入熱量を推算する熱風炉操業推算装置であって、前記複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための、目標送風温度、目標送風流量、目標送風時間、及び炉別サイクル別投入熱量を少なくとも入力とし、前記複数の熱風炉におけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度、炉別サイクル別送風流量、及び炉別サイクル別送風時間を少なくとも出力とする計算を行う熱風炉シミュレーション手段と、前記熱風炉シミュレーション手段による計算の結果に基づいて、前記複数の熱風炉を前記減風減温の状態にすることによる不利益が小さいほど評価が高くなることを示す評価値を算出する評価値算出手段と、前記評価値算出手段により算出された評価値に基づいて、前記炉別サイクル別投入熱量を決定する炉別サイクル別投入熱量決定手段と、を有し、前記炉別サイクル別の投入熱量決定手段は、前記評価値算出手段により算出された評価値に基づいて、前記炉別サイクル別投入熱量の変更終了条件を満足したか否かを判断し、当該変更終了条件を満足したと判断するまで、前記炉別サイクル別投入熱量を繰り返し変更して、前記熱風炉シミュレーション手段による計算の入力とし、前記炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の前記炉別サイクル別投入熱量を出力することを特徴とする熱風炉操業推算装置。前記複数の熱風炉における炉別残熱量を算出する炉別残熱量算出手段を更に有し、前記炉別サイクル別投入熱量は、前記炉別残熱量を基準とした投入熱量であり、前記炉別残熱量算出手段は、前記炉別サイクル別投入熱量、前記目標送風温度、前記目標送風流量、前記目標送風時間のうち、何れか1つのみを変更して前記熱風炉シミュレーション手段における計算の入力とし、入力した値と、当該熱風炉シミュレーション手段により得られた値とが等しくなり、且つ、前記熱風炉シミュレーション手段により得られた炉別サイクル別珪石煉瓦下端温度が基準温度と等しくなるという条件を満たす熱風炉を判別する判別手段と、前記条件を満たしたときに前記熱風炉シミュレーション手段により得られた値を用いて、前記判別した熱風炉の残熱量を算出する算出手段と、を更に有し、前記判別した熱風炉については、前記条件を満たしたときに前記熱風炉シミュレーション手段により得られた値を固定値として、前記判別手段による判別と、前記算出手段による算出と、を、前記判別手段により全ての熱風炉が前記条件を満たすと判別されるまで繰り返し行うことを特徴とする請求項1に記載の熱風炉操業推算装置。前記炉別残熱量算出手段は、前記炉別残熱量を繰り返し算出し、前記炉別サイクル別投入熱量は、前記炉別残熱量算出手段により算出された最新の炉別残熱量を基準とした投入熱量であることを特徴とする請求項2に記載の熱風炉操業推算装置。前記熱風炉シミュレーション手段は、前記炉別サイクル別送風温度、前記炉別サイクル別送風流量、及び前記炉別サイクル別送風時間に加えて、炉別サイクル別珪石煉瓦下端温度、及び炉別サイクル別排ガス温度を少なくとも出力とする計算を行い、前記評価値は、前記複数の熱風炉を減風減温の状態にすることにより前記高炉における出銑量が減ったことに起因するコストの損失と、前記複数の熱風炉の状態を減風減温の状態から目標の状態に変更するために使用するガスのコストと、が低いほど評価が高くなる評価値であり、前記評価値算出手段は、前記炉別サイクル別珪石煉瓦下端温度及び炉別サイクル別排ガス温度に基づく制約式を満足する範囲で評価が最も高くなる前記評価値を算出することを特徴とする請求項1~3の何れか1項に記載の熱風炉操業推算装置。前記複数の熱風炉は、スタッガードパラレル方式で操業されることを特徴とする請求項1~4の何れか1項に記載の熱風炉操業推算装置。燃焼ガスにより蓄熱煉瓦を加熱して蓄熱する燃焼期間と、当該蓄熱煉瓦に冷風を通して当該蓄熱煉瓦との熱交換により熱風を生成して高炉に供給する送風期間とを含む期間を1サイクルとしてそれぞれが稼働する複数の熱風炉の操業の状態を減風減温状態にした後、当該複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための炉別サイクル別の投入熱量を推算する熱風炉操業推算方法であって、前記複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための、目標送風温度、目標送風流量、目標送風時間、及び炉別サイクル別投入熱量を少なくとも入力とし、前記複数の熱風炉におけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度、炉別サイクル別送風流量、炉別サイクル別送風時間、及び炉別サイクル別珪石煉瓦下端温度を少なくとも出力とする計算を行う熱風炉シミュレーション工程と、前記熱風炉シミュレーション工程による計算の結果に基づいて、前記複数の熱風炉を前記減風減温の状態にすることによる不利益が小さいほど評価が高くなることを示す評価値を算出する評価値算出工程と、前記評価値算出工程により算出された評価値に基づいて、前記炉別サイクル別投入熱量を決定する炉別サイクル別投入熱量決定工程と、を有し、前記炉別サイクル別の投入熱量決定工程は、前記評価値算出工程により算出された評価値に基づいて、前記炉別サイクル別投入熱量の変更終了条件を満足したか否かを判断し、当該変更終了条件を満足したと判断するまで、前記炉別サイクル別投入熱量を繰り返し変更して、前記熱風炉シミュレーション工程による計算の入力とし、前記炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の前記炉別サイクル別投入熱量を出力することを特徴とする熱風炉操業推算方法。前記複数の熱風炉における炉別残熱量を算出する炉別残熱量算出工程を更に有し、前記炉別サイクル別投入熱量は、前記炉別残熱量を基準とした投入熱量であり、前記炉別残熱量算出工程は、前記炉別サイクル別投入熱量、前記目標送風温度、前記目標送風流量、前記目標送風時間のうち、何れか1つのみを変更して前記熱風炉シミュレーション工程における計算の入力とし、入力した値と、当該熱風炉シミュレーション工程により得られた値とが等しくなり、且つ、前記熱風炉シミュレーション工程により得られた炉別サイクル別珪石煉瓦下端温度が基準温度と等しくなるという条件を満たす熱風炉を判別する判別工程と、前記条件を満たしたときに前記熱風炉シミュレーション工程により得られた値を用いて、前記判別した熱風炉の残熱量を算出する算出工程と、を更に有し、前記判別した熱風炉については、前記条件を満たしたときに前記熱風炉シミュレーション工程により得られた値を固定値として、前記判別工程による判別と、前記算出工程による算出と、を、前記判別工程により全ての熱風炉が前記条件を満たすと判別されるまで繰り返し行うことを特徴とする請求項6に記載の熱風炉操業推算方法。前記炉別残熱量算出工程は、前記炉別残熱量を繰り返し算出し、前記炉別サイクル別投入熱量は、前記炉別残熱量算出工程により算出された最新の炉別残熱量を基準とした投入熱量であることを特徴とする請求項7に記載の熱風炉操業推算方法。前記熱風炉シミュレーション工程は、前記炉別サイクル別送風温度、前記炉別サイクル別送風流量、及び前記炉別サイクル別送風時間に加えて、炉別サイクル別珪石煉瓦下端温度、及び炉別サイクル別排ガス温度を少なくとも出力とする計算を行い、前記評価値は、前記複数の熱風炉を減風減温の状態にすることにより前記高炉における出銑量が減ったことに起因するコストの損失と、前記複数の熱風炉の状態を減風減温の状態から目標の状態に変更するために使用するガスのコストと、が低いほど評価が高くなる評価値であり、前記評価値算出工程は、前記炉別サイクル別珪石煉瓦下端温度及び炉別サイクル別排ガス温度に基づく制約式を満足する範囲で評価が最も高くなる前記評価値を算出することを特徴とする請求項6~8の何れか1項に記載の熱風炉操業推算方法。前記複数の熱風炉は、スタッガードパラレル方式で操業されることを特徴とする請求項6~9の何れか1項に記載の熱風炉操業推算方法。請求項1~5の何れか1項に記載の熱風炉操業推算装置の各手段としてコンピュータを機能させることを特徴とするコンピュータプログラム。ところで、高炉では、高炉における還元反応が適切に行われなくなり、高炉側の操業が不安定になること等、様々な要因で出銑量を一時的に低下させる場合がある。このような場合には、熱風炉から供給する熱風を減風減温させる熱風の減風減温の少なくとも何れか一方を行うことが望まれる以下の説明では、熱風炉から供給する熱風を減風減温させている状態を、必要に応じて減風減温状態と称する。ここで、減風とは、熱風炉から高炉へ供給する熱風の送風流量を低減させることであり、減温とは、熱風炉から高炉へ供給する熱風の温度を低減させることである。低コストで銑鉄の生産をするために、原料燃料を低コストのものにすると共に、使用する燃料が出来るだけ少なくなるように高炉の操業を行うことが求められる状況下では、高炉からの出銑量を一時的に低下させる頻度が多くなる。そうすると、熱風炉を減風減温状態にするメリットが増大する。熱風炉への投入熱量をより小さくすることができるからである。このようにして出銑量を低下させ、熱風炉の状態を減風減温状態にした後、出銑量を低下させなければならない要因が取り除かれ、出銑量を増加させる場合、熱風炉から供給する熱風を増風増温させる熱風の増風増温の少なくとも何れか一方を行う必要がある以下の説明では、熱風炉から供給する熱風を増風増温させている状態を、必要に応じて増風増温状態と称する。効率よく銑鉄を製造するためには、増加させる出銑量に応じた流量温度の熱風をいち早く高炉に供給する必要がある。しかしながら、前述した従来の技術では、熱風炉の状態を、減風減温状態から、高炉の要求に応じた状態にするまでの時間サイクル数を正確に把握することができなかった。このため、従来は、出銑量を一時的に低下させる場合であっても、熱風炉への投入熱量を大きく低下させずに、出銑量の増加を待つようにしていた。よって、減風減温状態のときの熱風炉への投入熱量を十分に低減させることができず、熱風炉に無駄な熱量を投入することになり、省コスト化を図ることが困難になる。以上のように、従来の技術では、熱風炉の状態を、減風減温状態から、高炉の要求に応じた状態にするまでの時間を正確に把握することができないという問題点があった。本発明は、このような問題点に鑑みてなされたものであり、熱風炉の状態を、減風減温状態から、高炉の要求に応じた状態にするまでのサイクル数時間及び各サイクルにおける投入熱量を正確に把握することを目的とする。本発明の熱風炉操業推算装置は、燃焼ガスにより蓄熱煉瓦を加熱して蓄熱する燃焼期間と、当該蓄熱煉瓦に冷風を通して当該蓄熱煉瓦との熱交換により熱風を生成して高炉に供給する送風期間とを含む期間を1サイクルとしてそれぞれが稼働する複数の熱風炉の操業の状態を減風減温状態にした後、当該複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための炉別サイクル別の投入熱量を推算する熱風炉操業推算装置であって、前記複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための、目標送風温度、目標送風流量、目標送風時間、及び炉別サイクル別投入熱量を少なくとも入力とし、前記複数の熱風炉におけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度、炉別サイクル別送風流量、及び炉別サイクル別送風時間を少なくとも出力とする計算を行う熱風炉シミュレーション手段と、前記熱風炉シミュレーション手段による計算の結果に基づいて、前記複数の熱風炉を前記減風減温の状態にすることによる不利益が小さいほど評価が高くなることを示す評価値を算出する評価値算出手段と、前記評価値算出手段により算出された評価値に基づいて、前記炉別サイクル別投入熱量を決定する炉別サイクル別投入熱量決定手段と、を有し、前記炉別サイクル別の投入熱量決定手段は、前記評価値算出手段により算出された評価値に基づいて、前記炉別サイクル別投入熱量の変更終了条件を満足したか否かを判断し、当該変更終了条件を満足したと判断するまで、前記炉別サイクル別投入熱量を繰り返し変更して、前記熱風炉シミュレーション手段による計算の入力とし、前記炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の前記炉別サイクル別投入熱量を出力することを特徴とする。本発明の熱風炉操業推算方法は、燃焼ガスにより蓄熱煉瓦を加熱して蓄熱する燃焼期間と、当該蓄熱煉瓦に冷風を通して当該蓄熱煉瓦との熱交換により熱風を生成して高炉に供給する送風期間とを含む期間を1サイクルとしてそれぞれが稼働する複数の熱風炉の操業の状態を減風減温状態にした後、当該複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための炉別サイクル別の投入熱量を推算する熱風炉操業推算方法であって、前記複数の熱風炉の操業を前記高炉の要求に応じた状態に変更するための、目標送風温度、目標送風流量、目標送風時間、及び炉別サイクル別投入熱量を少なくとも入力とし、前記複数の熱風炉におけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度、炉別サイクル別送風流量、炉別サイクル別送風時間、及び炉別サイクル別珪石煉瓦下端温度を少なくとも出力とする計算を行う熱風炉シミュレーション工程と、前記熱風炉シミュレーション工程による計算の結果に基づいて、前記複数の熱風炉を前記減風減温の状態にすることによる不利益が小さいほど評価が高くなることを示す評価値を算出する評価値算出工程と、前記評価値算出工程により算出された評価値に基づいて、前記炉別サイクル別投入熱量を決定する炉別サイクル別投入熱量決定工程と、を有し、前記炉別サイクル別の投入熱量決定工程は、前記評価値算出工程により算出された評価値に基づいて、前記炉別サイクル別投入熱量の変更終了条件を満足したか否かを判断し、当該変更終了条件を満足したと判断するまで、前記炉別サイクル別投入熱量を繰り返し変更して、前記熱風炉シミュレーション工程による計算の入力とし、前記炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の前記炉別サイクル別投入熱量を出力することを特徴とする。本発明のコンピュータプログラムは、前記熱風炉操業推算装置の各手段としてコンピュータを機能させることを特徴とする。本発明によれば、複数の熱風炉におけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度、炉別サイクル別送風流量、及び炉別サイクル別送風時間を少なくとも出力とする計算を行う。次に、その計算の結果に基づいて、複数の熱風炉を前記減風減温の状態にすることによる不利益が小さいほど評価が高くなることを示す評価値を算出する。次に、その算出した評価値に基づいて、複数の熱風炉の操業を高炉の要求に応じた状態に変更するための炉別サイクル別投入熱量を決定する。このようにして炉別サイクル別投入熱量を決定するに際し、炉別サイクル別投入熱量の変更終了条件を満足したか否かを評価値に基づいて判断し、当該変更終了条件を満足したと判断するまで、炉別サイクル別投入熱量を繰り返し変更して、前述して計算を繰り返し行い、炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の炉別サイクル別投入熱量を出力する。よって、熱風炉の状態を、減風減温の状態から、高炉の要求に応じた状態にするまでのサイクル数時間及び各サイクルにおける投入熱量を正確に把握することができる。以下、図面を参照しながら、本発明の一実施形態を説明する。[熱風炉100の構成]図1は、熱風炉100の概略構成の一例を示す図である。尚、各図では、説明の都合上、必要な部分のみを、必要に応じて簡略化して示す。図1において、熱風炉100は、不図示の高炉に熱風を供給するための蓄熱式熱交換器であり、高炉への送風に熱を与える蓄熱室101と、蓄熱室101を加熱するための燃焼室102と、熱風の温度調節を行うための混冷室103と、を有している。燃焼室102では、ガス供給ダクト112から吹き込まれるBFGとCOGとLDGとの混合ガス燃焼ガス及び燃焼空気供給ダクト113から吹き込まれる燃焼空気とを燃焼バーナ108で燃焼させ、この燃焼ガスを蓄熱室101の内部に積層された蓄熱煉瓦の間を通過させて加熱して熱を蓄える。図1に示す例では、この蓄熱煉瓦として、下側から順に、粘土煉瓦109、ハイアルミナ煉瓦110、シリカを主成分とする珪石煉瓦111が積層されており、これらの粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111には、上下方向に延びる複数の通過口が形成されている。ガス供給ダクト112には、ガス遮断弁130、ガスバタフライ弁131、及び燃焼ガス流量計132が設けられており、ガスバタフライ弁131を開閉することにより、燃焼室102に流入する燃焼ガスの流入量を調節することができる。燃焼空気供給ダクト113は、燃焼空気ファンから送風された空気を熱風炉100に送風する。燃焼空気供給ダクト113には、空気流量計127、空気バタフライ弁128、及び空気遮断弁129が設けられている。燃焼空気供給ダクト113には、燃焼ガスの流量に応じて、燃焼に必要な量の空気が流入されるようにしている。蓄熱室101の下端部には、ダクト114が設けられており、このダクト114は、N2、CO2等を含む燃焼ガスを排出するためのガス排出ダクト119と、ダクト114を介して蓄熱室101に冷風を供給するための冷風導入ダクト116と、に分岐される。ガス排出ダクト119には、煙道弁126が設けられている。冷風導入ダクト116には、送風弁124、及び送風バタフライ弁125が設けられており、送風バタフライ弁125を開閉させることにより、熱風炉100に流入する冷風の流入量を調節することができる。また、混冷室103には、高炉用の熱風を排出するための熱風排出ダクト117が接続されている。この熱風排出ダクト117には、熱風弁121が設けられている。また、冷風導入ダクト116の送風バタフライ弁125より上流側には、混冷室103に繋がるダクト118が設けられている。このダクト118には、冷風弁122と、冷風バタフライ弁123とが設けられている。図2は、熱風炉100における燃焼期間と送風期間の動作の概要の一例を示す図である。尚、以下の説明では、BFGガスとCOGガスとLDGガスとを燃焼ガスと総称する。図2aに示すように、燃焼期間において、蓄熱室101に熱を蓄える場合には、送風弁124、冷風弁122、及び熱風弁121を完全に閉じて、ガス供給ダクト112及び燃焼空気供給ダクト113を介して燃焼室102内に燃焼ガス及び燃焼空気を流入させる。これらの燃焼ガス及び燃焼空気はバーナ108によって燃焼され、この燃焼ガスは、蓄熱室101の粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111に形成された開口部を通って、粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111を蓄熱する。粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111を通過した燃焼ガスは、ガス排出ダクト119を介して排ガスとして煙道に排出される。ここで、通常は、珪石煉瓦111の最下部での最低温度は変態点温度以下とならないように管理される。また、粘土煉瓦109の最下部での温度の下限値は排ガス温度が高くならないようにできるだけ低く一定値に管理される。本実施形態では、珪石煉瓦111の下端部の最低温度の管理値を基準温度と記す。蓄熱室101への蓄熱が完了すると、図2bに示すように、煙道弁126、空気遮断弁129、及びガス遮断弁130を完全に閉じて、冷風導入ダクト116を介して蓄熱室101に冷風を流入させる。蓄熱室101に流入した冷風は、粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111に形成された開口部を通過して900~1300°Cに加熱された後、高炉用の熱風として熱風排出ダクト117から排出される。図3は、熱風炉の制御システムの概略構成の一例を示す図である。図3において、実線は信号の流れを示し、破線は、冷風、熱風、燃焼ガス、燃焼空気の流れを示している。図3では、1基の高炉に対して4基の熱風炉100a~100dを付帯させた場合を例に挙げて示している。また、これら4基の熱風炉100a~100dは、スタッガードパラレル方式で操業するものとする。図4は、スタッガードパラレル方式における操業スケジュールの概略の一例を説明する図である。図4に示す例では、左から右へ送風から燃焼への切替と、燃焼と、燃焼から送風への切替と、送風とをこの順番で行い、これらの期間を合わせた期間で1サイクルを構成するようにしている図4に示す1サイクル=切り替え時間401燃焼時間402切り替え時間401送風時間403の部分を参照。そして、熱風を供給する順番で前後隣接する2基例えば熱風炉1と熱風炉2の送風時間の一部をラップさせるようにする。更に、図4に示す例では、簡単のために送風時間と燃焼時間とを同じ長さにしているので、熱風を供給する順番で隣接しない2基例えば熱風炉1と熱風炉3については、一方の熱風炉が送風期間であるときに他方の熱風炉が燃焼期間となり、一方の熱風炉が燃焼期間であるときに他方の熱風炉が送風期間となるようにする。また、本実施形態では、全てのサイクルにおいて、送風時間が同じ時間である場合を例に挙げて説明する。[制御装置300の構成]図3の説明に戻り、熱風炉の制御システムは、熱風炉100a~100dの操業を制御する制御装置300を有している。制御装置300は、熱風炉制御計算機301と、入出力装置302と、流量調節計304と、温度調節計305と、開度調節計313a~313dと、を有している。熱風炉制御計算機301は、予め設定される入力データをもとに熱風炉100a~100dの操業を行うための演算を行う。熱風炉制御計算機301に対するデータの入力は、インターフェース部である入出力装置302を介して行われる。入出力装置302に入力される情報としては、例えば、以下のようなものがある。まず、熱風炉のオペレータによる入力装置312の操作に基づく情報として、以下の情報がある。熱風炉100a~100dの操業を、減風減温状態から、高炉からの要求に応じた状態に変更する際の目標値目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemただし、目標送風時間Btimemは固定値。設備負荷炉別サイクル別熱効率ηin、炉別サイクル別投入熱量Qinin、炉別サイクル別蓄熱量Qbfin、炉別サイクル別蓄熱分担比率αinの初期値。操業制約設備制約後述する制約条件の上下限値熱風炉100を構成する設備粘土煉瓦109、ハイアルミナ煉瓦110、珪石煉瓦111等ガス燃焼ガス燃焼空気等の既知データ物理定数設備定数。この他にも、後述する計算において予めオペレータが設定する情報も入出力装置302に入力される。ここで、添字i添字に,が付されているものについては、,の後に表示されているiは、熱風炉100を表す。本実施形態では、熱風炉100a~100dは4基なので、iは1~4となる。また、添字nは、サイクルを表し、添字mは、目標値であることを表す。また、操業実績を示す情報としては、例えば、以下のようなものがある。熱風炉100a~100dから高炉に排出される熱風の温度送風温度を測定する熱風温度計310により測定された送風温度BT。珪石煉瓦111の下端の温度珪石煉瓦下端温度を炉毎に測定する珪石煉瓦温度計315a~315dにより測定された炉別珪石煉瓦下端温度Tsii。送風機306から送風された冷風の流量冷風流量を測定する送風流量計307で測定された冷風流量=送風流量。送風機306から送風された冷風の温度を測定する冷風温度計308により測定された冷風温度。また、図示を省略するが、各熱風炉100a~100dに付帯するその他センサから出力されるその他の操業実績を示す情報も入出力装置302に入力される。本実施形態では、熱風炉制御計算機301は、入力された情報等をもとに、投入熱量計算プログラムを実行して、各熱風炉100a~100dに投入する投入熱量を熱風炉100a~100d毎に個別に計算し、その結果を用いて、燃焼ガスの流量を熱風炉100a~100d毎に個別に計算する。そして、熱風炉制御計算機301は、各熱風炉100a~100dに対応して設けられている流量調節計304a~304dに対して、燃焼ガスの流量を指示する。流量調節計304a~304dは、熱風炉制御計算機301から指示された流量が指示通りになるように、ガスバタフライ弁131a~131dの開閉動作を調節する。また、熱風炉制御計算機301は、オペレータによる入力装置312の操作に基づいて得られた温度を、熱風炉100a~100dから排出される熱風の温度の目標値として、温度調節計305に設定する。温度調節計305は、送風バタフライ弁125a~125dの開閉動作を調節する開度調節計313a~313dに対して、目標開度を指示する。例えば、2基の熱風炉100で送風している期間であれば、温度調節計305は、送風温度が高ければ、先行して送風する熱風炉100の送風バタフライ弁125の開度を開けて、後行して送風する熱風炉100の送風バタフライ弁125の開度を閉める。逆に、送風温度が低ければ、温度調節計305は、先行して送風する熱風炉100の送風バタフライ弁125の開度を閉めて、後行して送風する熱風炉100の送風バタフライ弁125の開度を開ける。先行して送風する熱風炉100は、熱交換が進んでいるため、熱風温度が低く、後行して送風する熱風炉100は、十分に蓄熱しているため、熱風温度が高いので、これらの送風バタフライ弁125の開度風量の配分を調節することによって送風温度を制御することができる。また、1基の熱風炉100で送風している期間中に送風温度が高い場合は冷風バタフライ弁123を開けて冷風を熱風に混入させて送風温度を調節する。1基の熱風炉100で送風している期間では、送風バタフライ弁の開度は一定にしておく。送風機306は、送風流量調節計309に設定された風量の冷風を送風する。送風流量調節計309の設定は、送風工場で変更される。また、熱風炉制御計算機301は、以上のような計算結果に基づく画面の表示を、入出力装置302を介して表示装置303に行わせる。特に、本実施形態では、熱風炉制御計算機301は、熱風炉100a~100dを減風減温状態にした後、増風増温を行って、熱風炉100a~100dの操炉条件を目標値目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemに変更するために必要な炉別サイクル別投入熱量Qininを算出し、オペレータに報知するための処理を行うことを特徴としている。したがって、以下では、このような処理の説明に必要な部分のみの構成及び動作を説明し、本実施形態の説明に不要な部分の構成及び動作の詳細な説明を省略する。尚、以下の説明では、前述した目標値を必要に応じて要求操炉条件と称する。熱風炉制御計算機301は、例えば、CPU、ROM、RAM、HDDを備えたコンピュータを用いることにより実現することができる。前述したプログラムは、例えば、HDDに記憶され、CPUにより実行される。[熱風炉制御計算機301の機能構成]概要図5は、熱風炉制御計算機301の機能的な構成の一例を示す図である。図5に示すように、熱風炉制御計算機301は、炉別残熱量算出部501と、炉別サイクル別投入熱量決定部502と、熱風炉シミュレーション部503と、評価値算出部504と、炉別サイクル別投入熱量表示指示部505と、を有する。図6は、熱風炉制御計算機301が行う処理の一例を概念的に説明する図である。図6において、Qinsは、減風減温状態にする前の各熱風炉100a~100dへの投入熱量総量である添字sは、減風減温状態にする前であることを表す。Qminは、各熱風炉100a~100dが最低限保っていなければならない最低熱量総量を表す。最低熱量Qminは、熱風炉100に備わる煉瓦の特性や、操業条件等によって定まるものである。Qrmは、最低熱量Qminを基準0としたときの各熱風炉100a~100dの熱量残熱量である。尚、図6では、残熱量Qrmが一定であるように示しているが、残熱量Qrmは、時々刻々と変化するものである。Qin1~3図6の斜線部は、それぞれ、増風増温を開始してから、1サイクル目、2サイクル目、3サイクル目における各熱風炉100a~100dへの投入熱量から高炉への送風で失われた各熱風炉の熱量を引いたいわゆる各熱風炉蓄熱量増量分である。図6に示す例では、増風増温を開始してから3サイクルで要求操炉条件に従う状態になっている。尚、図6に示すサイクル1~3は、基準となる熱風炉100例えば熱風炉100aにおけるサイクルを表す。例えば、基準となる熱風炉100を熱風炉100aとした場合、熱風炉100aのサイクル1の期間において、各熱風炉100a~100dに投入される投入熱量の総量がQin1となる。本実施形態では、熱風炉制御計算機301は、計算を開始するタイミングにおける各熱風炉100a~100dの残熱量Qrmiを計算し、この残熱量Qrmiを基準0として、増風増温を開始してから要求操炉条件に従う状態になるまでの間の、炉別サイクル別投入熱量Qininの最適値を計算する。熱風炉シミュレーション部503熱風炉シミュレーション部503は、熱風炉シミュレータを実行する。熱風炉シミュレータは、熱物理学を活用して、燃焼送風によって、熱風炉100で行われる熱交換を忠実に再現する熱風炉モデルを用いて構築される。図7は、熱風炉モデルの一例を概念的に示す図である。ここでは、蓄熱室101における熱風炉モデルの一例を示す。本実施形態で説明する熱風炉モデルでは、蓄熱室101の形が円筒であるとする。そして、その円筒を軸に沿って切った断面二次元平面における熱収支を表すモデル式を構築する。図7に示すように、例えば、珪石煉瓦111は、蓄熱煉瓦珪石701と、炉壁煉瓦702とにより表され、この炉壁煉瓦702は、更に、鉄皮702a、耐火煉瓦A~C702b~702d、キャスタブル702e、鉄皮702fにより構成される。珪石煉瓦111のモデル式では、これらについての前述した二次元平面の伝熱が表現される。他の粘土煉瓦109及びハイアルミナ煉瓦110についても、珪石煉瓦111と同様に、前述した二次元平面の伝熱が表現される。本実施形態の熱風炉モデルでは、燃焼時において、燃焼ガスと蓄熱煉瓦珪石701との間と、蓄熱煉瓦珪石701と炉壁煉瓦702との間の伝熱は、対流熱伝達と輻射熱伝達によるものとする。また、送風時において、冷風と蓄熱煉瓦珪石701との間と、蓄熱煉瓦珪石701と炉壁702との間と、炉壁702と大気との間の伝熱は、対流熱伝達によるものとする。また、炉壁702内は、熱伝導による伝熱がなされるものとする。また、蓄熱煉瓦珪石701は、通過口が存在することによる寸法換算を行う。また、高さ方向の熱伝導は無視するものとする。これらの仮定は、他の粘土煉瓦109及びハイアルミナ煉瓦110についても同じである。まず、ガスの熱収支を表すモデル式は、以下の1式で表される。また、蓄熱煉瓦の熱収支を表すモデル式は、以下の2式で表される。また、炉壁煉瓦の熱収支を表すモデル式は、以下の3式で表される。1式~3式における記号の意味は、以下の通りである。ρ:密度[kg/m3]Cp:比熱[J/kgK]S:断面積[m2]前述した円筒をその軸に沿って切った断面積から通過口の断面積の総和を引いた煉瓦の正味の切り口面積v:流速[m/s]T:温度[K]h:熱伝達率[W/m2K]A:接触面積[m2]L:高さ[m]炉底を0とするε:放射率[-]z:高さ方向の位置[m]t:時間[s]G:ガスCH:蓄熱煉瓦W:炉壁煉瓦Wb:炉壁煉瓦の構成煉瓦珪石煉瓦111であれば、耐火煉瓦珪石702a、耐熱煉瓦A~C702b~702d、キャスタブル702e、鉄皮702fWb-1、Wb1:Wbで特定される煉瓦の隣の煉瓦耐熱煉瓦AがWbであれば、Wb-1は耐火煉瓦珪石702aであり、Wb1は、耐熱煉瓦Bδ:ステファンボルツマン定数熱風炉シミュレーション部503は、1式~3式のモデル式を使って、以下に示す機能を実現するものである。図8は、熱風炉シミュレーション部503の詳細な機能構成の一例を示す図である。熱風炉シミュレーション部503は、燃焼時シミュレーション部801と、送風時シミュレーション部802とに大別される。燃焼時シミュレーション部801は、燃焼期間における各熱風炉100a~100bの状態をシミュレーションするためのものである。燃焼時シミュレーション部801は、使用ガス容量算出部801aと、モデル式計算部801bとを有する。<使用ガス容量算出部801a>使用ガス容量算出部801aは、炉別サイクル別投入熱量Qininを入力する。炉別サイクル別投入熱量Qininは、全投入熱量Qinに、炉別サイクル別蓄熱分担率αinΣαin=1を乗じることにより算出することもできる。よって、全投入熱量Qinと、炉別サイクル別蓄熱分担率αinとを入力することも、炉別サイクル別投入熱量Qininを入力することと等価である。また、使用ガス容量算出部801aは、BFG、COG、及びLDGの各使用ガスの使用比率[-]を入力する。この使用ガスの使用比率は、燃焼期間中に変化するものとし、1つの燃焼期間を分割した複数の期間毎に個別に設定することができる。さらに、使用ガス容量算出部801aは、BFG、COG、LDG、及び空気の各使用ガスのガスカロリーガス組成単位熱量[J/Nm3]を入力する。使用ガス容量算出部801aは、以上の入力した情報に基づいて、BFG、COG、LDG、及び空気の各使用ガスの容量[Nm3]を算出してモデル式計算部801bに出力する。<モデル式計算部801b>モデル式計算部801bは、使用ガス容量算出部801aから、BFG、COG、LDG、及び空気の各ガスの容量を入力する。モデル式計算部801bは、ガス燃焼温度TG[°C]、密度ρG[kg/m3]、比熱CpG[J/kgK]を入力する。ここで、添字Gは、ガスを表す。そして、モデル式計算部801bは、入力した情報を用いて、前述したモデル式1式~3式を使用して、燃焼開始から燃焼終了まで1つの燃焼期間の、熱風炉100a~100dにおける熱収支の計算を行い、燃焼終了の時点での以下の値を出力する。炉別サイクル別蓄熱煉瓦温度TCHin炉別サイクル別炉壁煉瓦の構成煉瓦温度TWbin炉別サイクル別排ガス温度Texin炉別サイクル別燃焼効率ηin尚、ガスボリュームVGは、1式の左辺のSG×vGに反映されるものである。また、炉別サイクル別排ガス温度Texinは、計算対象の燃焼期間において蓄熱されなかった熱量と、ガスボリュームVGと、ガスの成分とを用いて算出されるものである。以上の使用ガス容量算出部801aとモデル式計算部801bの処理は、サイクル毎に繰り返し行われる。尚、サイクル数は、使用ガス容量算出部801aが入力した炉別サイクル別投入熱量Qininで識別される。送風時シミュレーション部802は、送風期間における各熱風炉100a~100bの状態をシミュレーションするためのものである。送風時シミュレーション部802は、各炉運転内容指示部802aと、モデル式計算部802bと、送風流量温度決定部802cとを有する。<各炉運転内容指示部802a>各炉運転内容指示部802aは、モデル式計算部801bから、炉別サイクル別蓄熱煉瓦温度TCHinと炉別サイクル別炉壁煉瓦の構成煉瓦温度TWbinとを入力する度に燃焼時シミュレーション部801において、燃焼期間における燃焼終了時の計算が終了する度に計算を開始する。ここでは、1秒毎の計算結果が得られるように繰り返し計算が行われるものとする。各炉運転内容指示部802aは、後述する送風流量温度決定部802cで決定された炉別時間別通過送風流量BHSiと、炉別時間別混冷流量Airiを、計算対象の時間の炉別時間別通過送風流量BHSitと、炉別時間別混冷流量Airitとする。そして、各炉運転内容指示部802aは、炉別時間別通過送風流量BHSitをモデル式計算部802bに出力すると共に、炉別時間別混冷流量Airitを送風流量温度決定部802cに出力する。ここで、通過送風流量とは、図1に示すダクト114、蓄熱室101、燃焼室102を通って混冷室103に進入した冷風の流量である。また、混冷流量とは、図1に示すダクト118を通って混冷室103に進入した冷風の流量である。これらの通過送風流量と混冷流量との加算値が送風流量BVitになる。また、前記において、添字tは時間を表す。尚、最初の計算では、後述する送風流量温度決定部802cで炉別時間別通過送風流量BHSiと、炉別時間別混冷流量Airiが決定されていないので、各炉運転内容指示部802aは、炉別時間別通過送風流量BHSitと、炉別時間別混冷流量Airitの初期値を採用する。<モデル式計算部802b>モデル式計算部802bは、モデル式計算部801bから、蓄熱煉瓦の温度TCHinと炉壁煉瓦の構成煉瓦の温度TWbinとを入力する。また、モデル式計算部802bは、密度ρG、及び比熱CPGを入力する。また、モデル式計算部802bは、各炉運転内容指示部802aから、炉別時間別通過送風流量BHSitを入力する。ここで、炉別時間別通過送風流量BHSitに時間tを掛けた値がガスボリュームVGとなる。そして、モデル式計算部802bは、入力した情報を用いて、前述したモデル式1式~3式を使用して、送風期間における計算対象の時間tでの熱収支の計算を繰り返し行い、炉別時間別ガス温度TGitを、送風流量温度決定部802cに出力する。尚、燃焼時と送風時とでは、モデル式1式及び2式熱伝達率hGCHの値が異なる。<送風流量温度決定部802c>送風流量温度決定部802cは、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを入力する。また、送風流量温度決定部802cは、各炉運転内容指示部802aから、炉別時間別混冷流量Airitを入力する。また、送風流量温度決定部802cは、モデル式計算部802bから、炉別時間別のガスの温度TGitを入力する。また、送風流量温度決定部802cは、冷風温度TCWを入力する。そして、送風流量温度決定部802cは、入力した情報を用いて、以下の4式~6式により、目標送風温度BTm及び目標送風流量BVmが得られるように、炉別時間別送風流量BVi、炉別時間別通過送風流量BHSi、炉別時間別混冷流量Airiを算出する。そして、算出された炉別時間別通過送風流量BHSiと、炉別時間別混冷流量Airiは、前述したように、各炉運転内容指示部802aに出力される。すなわち、炉別時間別通過送風流量BHSiと、炉別時間別混冷流量Airiは、次の計算対象の時間における炉別の通過送風流量混冷流量となる。送風流量温度決定部802cは、送風期間の終了時送風終了時に得られた炉別時間別のガスの温度TGit、炉別時間別通過送風流量BHSit、炉別時間別混冷流量Airitを、4式~6式の右辺に与えて、炉別サイクル別送風温度BTkin及び炉別サイクル別送風流量BVkinを算出する。尚、本実施形態では、炉別サイクル別送風時間BTimekinは、固定値である。また、添字kは結果を表す。また、送風流量温度決定部802cは、送風期間の終了時送風終了時に得られた炉別サイクル別蓄熱煉瓦温度TCHinの中から珪石煉瓦下端部の温度を抽出し、抽出した温度を炉別サイクル別珪石煉瓦下端温度Tsiinとする。そして、送風流量温度決定部802cは、これらの炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinを出力する。また、送風流量温度決定部802cは、燃焼効率ηinを算出して出力する。燃焼効率ηinは、燃焼終了時の炉別サイクル別投入熱量Qininと、送風終了時に得られる送風期間中の炉別サイクル別送風温度BTin、炉別サイクル別送風流量BVinによって算出される炉別サイクル別煉瓦畜熱量Qbfinから、以下の7式で算出されるものである。ここで、炉別サイクル別煉瓦畜熱量Qbfinは、以下の8式で与えられる。8式の記号の意味は、前述した通りである。炉別残熱量算出部501図5の説明に戻り、炉別残熱量算出部501は、計算を開始する時刻における炉別残熱量Qrmiを算出する。図9は、炉別残熱量算出部501の詳細な機能構成の一例を示す図である。図9において、炉別残熱量算出部501は、第1の煉瓦下端温度判定部901と、第2の煉瓦下端温度判定部902と、第1の送風流量変化分決定部903と、第1の残熱量判定部904と、該当炉残熱量算出部905と、残熱量記憶部906と、第2の残熱量判定部907と、第2の送風流量変化分決定部908と、第3の送風流量変化分決定部909と、第3の残熱量判定部910と、を有する。第1の煉瓦下端温度判定部901は、珪石煉瓦温度計315a~315dにより測定された炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えているか否かを判定する。第2の煉瓦下端温度判定部902は、第1の煉瓦下端温度判定部901により、炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えていると判定されると、珪石煉瓦温度計315a~315dにより測定された炉別珪石煉瓦下端温度Tsiiのうち、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがあるか否かを判定する。この判定の結果、第2の煉瓦下端温度判定部902は、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがあると判定すると、当該炉別珪石煉瓦下端温度Tsiiを有する熱風炉iの残熱量Qrmiは、最低熱量Qminとなる0となる。よって、残熱量記憶部906は、当該熱風炉iの残熱量Qrmiとして0を記憶する。第1の送風流量変化分決定部903は、第2の煉瓦下端温度判定部902により、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがないと判定されると、例えば二分法に基づき、目標送風流量BVmの増加分ΔBVを決定する。そして、第1の送風流量変化分決定部903は、目標送風流量BVmに増加分ΔBVを加算したものを新たな目標送風流量BVmとする。第1の送風流量変化分決定部903は、この新たな目標送風流量BVmを炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。この出力に応じて、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinが、熱風炉シミュレーション部503から第1の残熱量判定部904に出力される。第1の残熱量判定部904は、熱風炉シミュレーション部503から出力された炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinを入力する。そして、第1の残熱量判定部904は、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがあるか否かを判定する。この判定の結果、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがなければ、第1の送風流量変化分決定部903は、例えば二分法に基づき、目標送風流量BVmの増加分ΔBVを更に決定して新たな目標送風流量BVmを決定する。そして、第1の送風流量変化分決定部903は、新たな目標送風流量BVmを熱風炉シミュレーション部503に出力する。そうすると、熱風炉シミュレーション部503から、新たな炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinが第1の残熱量判定部904に入力される。すなわち、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iが見つかるまで、第1の送風流量変化分決定部903と、第1の残熱量判定部904との処理を繰り返し行う。該当炉残熱量算出部905は、第1の残熱量判定部904により、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがあると判定されると、当該熱風炉iに対して熱風炉シミュレーション部503により算出された炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekinを用いて、当該熱風炉iの残熱量Qrmiを算出する。残熱量記憶部906は、算出された当該熱風炉iの残熱量Qrmiを記憶する。第2の残熱量判定部907は、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしたか否かを判定する。全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしていると、全ての熱風炉100a~100dの残熱量Qrmiの計算が終了したことになる。これら熱風炉100a~100dの残熱量Qrmiは、残熱量記憶部906から炉別サイクル別投入熱量決定部502図5を参照に出力される。第2の送風流量変化分決定部908は、第2の残熱量判定部907により、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしていない場合に、次の処理を行う。まず、第2の送風流量変化分決定部908は、第1の残熱量判定部904により、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしていると判定された熱風炉iに対して熱風炉シミュレーション部503により算出された炉別サイクル別送風流量BVkinを固定値とするように、熱風炉シミュレーション部503に指示する。このとき、第2の送風流量変化分決定部908は、複数の熱風炉iについて、炉別サイクル別送風流量BVkinを固定値とした場合には、その総和ΣBVkinを計算する。そして、第2の送風流量変化分決定部908は、現在の目標送風流量BVmから、固定値とされた炉別サイクル別送風流量BVkinの総和ΣBVkinを減算した送風流量BV ́BV ́=BVm-ΣBVkinを基準に、目標送風流量BVmの増減分ΔBVを決定する。そして、第2の送風流量変化分決定部908は、目標送風流量BVmに増加分ΔBVを加算又は減算したものを新たな目標送風流量BVmとし、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。この出力に応じて、熱風炉シミュレーション部503から、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinが、第1の残熱量判定部904に出力される。そして、前述したように、第1の残熱量判定部904は、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがあるか否かを判定する。第3の送風流量変化分決定部909は、第1の煉瓦下端温度判定部901により、炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えていないと判定すると、例えば二分法に基づき、目標送風流量BVmの減少分ΔBVを決定する。そして、第3の送風流量変化分決定部909は、目標送風流量BVmに減少分ΔBVを減算したものを新たな目標送風流量BVmとし、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。この出力に応じて、熱風炉シミュレーション部503から、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinが、第3の残熱量判定部910に出力される。第3の残熱量判定部910は、熱風炉シミュレーション部503から出力された炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinを入力する。そして、第3の残熱量判定部910は、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たしたか否かを判定する。この判定の結果、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たしていななければ、第3の送風流量変化分決定部909は、例えば二分法に基づき、目標送風流量BVmの減少分ΔBVを決定して新たな目標送風流量BVmを決定する。そして、第3の送風流量変化分決定部909は、新たな目標送風流量BVmを熱風炉シミュレーション部503に出力する。そうすると、熱風炉シミュレーション部503から、新たな炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinが第3の送風流量変化分決定部909に入力される。すなわち、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たすまで、第1の送風流量変化分決定部903と、第1の残熱量判定部904との処理を繰り返し行う。そして、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たすと、第3の残熱量判定部910は、そのときに熱風炉シミュレーション部503から出力された炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiinを第1の残熱量判定部904に出力する。そして、前述したように、第1の残熱量判定部904は、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがあるか否かを判定する。以上の第1の煉瓦下端温度判定部901と、第2の煉瓦下端温度判定部902と、第1の送風流量変化分決定部903と、第1の残熱量判定部904と、該当炉残熱量算出部905と、第2の残熱量判定部907と、第2の送風流量変化分決定部908と、第3の送風流量変化分決定部909と、第3の残熱量判定部910の処理は、全ての熱風炉100a~100dの残熱量Qrmiの計算が終了するまで繰り返し行われる。評価値算出部504図5の説明に戻り、評価値算出部504は、熱風炉シミュレーション部503から、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texinを入力する。そして、評価値算出部504は、入力した情報を用いて、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度よりも大きいという制約条件制約式と、炉別サイクル別排ガスの温度Texinが基準温度よりも小さいという制約条件制約式を満たす範囲で、最も小さな評価値Yを算出する。評価値Yは、以下の9式で定められるものである。Y=fx1,x299式の右辺は、決定変数x1、x2とする評価関数を表す。また、これら決定変数x1、x2は、それぞれ以下の10式、11式で表される。ここで、10式において、cost1は、出銑量単位重さ当たりの収益である。また、IRONmは、目標操業時の出銑量である。また、IRONnは、nサイクル目の炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekinに基づいて操業した場合の出銑量である。また、11式において、cost2は、ガス単位量当たりのコストである。以上のように、決定変数x1は、減風減温状態にしたことによる高炉の出銑機会の損失コストの損失を評価するものである。一方、決定変数x2は、熱風炉100a~100dにおいて畜熱量を増やすために使用するガスのコストを評価するためのものである。炉別サイクル別投入熱量決定部502炉別サイクル別投入熱量決定部502は、減風減温の状態から、要求操炉条件に応じた状態高炉からの要求に応じた状態に変更するための炉別サイクル別投入熱量Qininを決定するものであり、図5に示す各ブロックの中枢となる部分である。尚、要求操炉条件は、熱風炉100a~100dを、減風減温状態にする前の元の状態に戻すという条件でも、そうでない条件でもよい。まず、炉別サイクル別投入熱量決定部502は、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを入力すると、減風減温状態から、高炉からの要求に応じた状態に変更するための炉別サイクル別投入熱量Qininを決定する処理を開始する。また、炉別サイクル別投入熱量決定部502は、現在の炉別残熱量Qrmiの計算を炉別残熱量算出部501に指示して、炉別残熱量算出部501から現在の炉別残熱量Qrmiを入力する。また、炉別サイクル別投入熱量決定部502は、設備負荷炉別サイクル別熱効率ηin、炉別サイクル別投入熱量Qinin、炉別サイクル別蓄熱量Qbfin、炉別サイクル別蓄熱分担比率αinの初期値を入力する。これらは、現在の操業に応じたものであり、以下に説明する繰り返し計算の最初の計算で使用されるものである。ここで、炉別サイクル別投入熱量Qininは、現在の炉別残熱量Qrmiを基準とする値であるとする。また、炉別サイクル別投入熱量決定部502は、評価値算出部504から評価値Yを入力して、その履歴を記憶しておく。炉別サイクル別投入熱量決定部502は、前回、熱風炉シミュレーション部503に出力した炉別サイクル別投入熱量Qininと、評価値Yの履歴とに基づいて、より値が低い評価値Yを得るために、炉別サイクル別投入熱量Qininのそれぞれをどれだけ変更すればよいのかを、統計的に判断する。この判断のロジックは、GA、山登り法、及び線形計画法等、公知の方法で実現することができる。炉別サイクル別投入熱量Qininの変更量は、このロジックに応じて定まるものである。ここで、前述したように、炉別サイクル別投入熱量Qininは、現在の炉別残熱量Qrmiを基準とする値である。すなわち、炉別サイクル別投入熱量Qininは、現在の炉別残熱量Qrmiに対して、各熱風炉100a~100dにどれだけの熱量の熱を投入すればよいのかをサイクル毎に定めたものである。本実施形態では、炉別サイクル別投入熱量Qininは、現在の炉別残熱量Qrmiの入力があってから、炉別サイクル別投入熱量Qininの変更を行う。そして、炉別サイクル別投入熱量決定部502は、変更後の炉別サイクル別投入熱量Qininを、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。これにより、熱風炉シミュレーション部503で、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texinが得られ、これらを用いて、評価値算出部504で、評価値Yが得られる。炉別サイクル別投入熱量決定部502は、炉別サイクル別投入熱量Qininの変更終了条件を満足したか否かを判定する。この判定は、評価値Yが一定値以上変化しなくなったか否かすなわち収束したか否か、又は、繰り返しの計算の回数評価値Yの個数が所定回数になったか否か等により行うことができる。炉別サイクル別投入熱量決定部502は、炉別サイクル別投入熱量Qininの変更終了条件を満足したときに得られている最新の炉別サイクル別投入熱量Qininを、炉別サイクル別投入熱量表示指示部505に出力する。炉別サイクル別投入熱量表示指示部505炉別サイクル別投入熱量表示指示部505は、炉別サイクル別投入熱量決定部502から出力された炉別サイクル別投入熱量Qininを表示装置303図3を参照に表示するための表示データを生成して、入出力装置302を介して表示装置303に出力する。これにより、オペレータは、減風減温状態から、要求操炉条件目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemに戻すのに、何サイクル必要であり、各サイクルにおいて各熱風炉100a~100dにどのくらいの熱量を投入すればよいのかを知ることができる。尚、オペレータが、このことを把握できる情報を表示すれば、どのような情報を表示するようにしてもよい。[フローチャート]次に、図10のフローチャートを参照しながら、熱風炉制御計算機301の処理の一例を説明する。まず、ステップS1001において、炉別サイクル別投入熱量決定部502は、要求操炉条件に応じた状態高炉からの要求に応じた状態に変更するための炉別サイクル別投入熱量Qininを決定する計算を開始するか否かを判定する。この判定は、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを入力したか否かによって行われる。この判定の結果、計算を開始しないと判定した場合には、熱風炉制御計算機301が行うその他の公知の処理を行って、図10のフローチャートによる処理を終了する。一方、計算を開始すると判定した場合には、ステップS1002に進む。ステップS1002に進むと、炉別残熱量算出部501による炉別残熱量算出処理が行われる。この炉別残熱量算出処理により、現在の炉別残熱量Qrmiが得られる。この炉別残熱量算出処理の詳細については後述する。次に、ステップS1003において、炉別サイクル別投入熱量決定部502は、炉別サイクル別投入熱量Qininの初期値を設定する。炉別サイクル別投入熱量Qininは、現在の炉別残熱量Qrmiを基準とした値である。そして、炉別サイクル別投入熱量決定部502は、目標送風温度BTm、目標送風流量BVm、及び目標送風時間BTimemと、炉別サイクル別投入熱量Qininと、を熱風炉シミュレーション部503に出力する。次に、ステップS1004において、熱風炉シミュレーション部503による熱風炉シミュレータ実行処理が行われる。この熱風炉シミュレータ実行処理により、ここでは、少なくとも、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガス温度Texinが得られる。この熱風炉シミュレータ実行処理の詳細については後述する。次に、ステップS1005において、評価値算出部504は、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度よりも大きいという制約条件制約式と、炉別サイクル別排ガスの温度Texinが基準温度よりも小さいという制約条件制約式を満たす範囲で、最も小さな評価値Y9式を参照を算出する。次に、ステップS1006において、炉別サイクル別投入熱量決定部502は、炉別サイクル別投入熱量Qininの変更終了条件を満足したか否かを判定する。この判定の結果、変更終了条件を満足した場合には、後述するステップS1009に進む。一方、変更終了条件を満足していない場合には、ステップS1007に進む。ステップS1007に進むと、炉別残熱量算出部501による炉別残熱量算出処理が行われる。この炉別残熱量算出処理により、現在の炉別残熱量Qrmiが得られる。この炉別残熱量算出処理の詳細については後述する。次に、ステップS1008において、サイクル別投入熱量決定部502は、前回、熱風炉シミュレーション部503に出力した炉別サイクル別投入熱量Qininと、評価値Yの履歴とに基づいて、炉別サイクル別投入熱量Qininのそれぞれをどれだけ変更すればよいのかを、統計的に判断する。そして、サイクル別投入熱量決定部502は、判断した炉別サイクル別投入熱量Qininの変更量に応じて新たな炉別サイクル別投入熱量Qininを決定する。この新たな炉別サイクル別投入熱量Qininは、ステップS1007で得られた現在の炉別残熱量Qrmiを基準とした値である。サイクル別投入熱量決定部502は、この新たな炉別サイクル別投入熱量Qininと、目標送風温度BTm、目標送風流量BVm、及び目標送風時間BTimemと、を熱風炉シミュレーション部503に出力する。そして、ステップS1004に戻り、変更終了条件を満足するまで、ステップS1004~S1008の処理を繰り返し行う。そして、前述したように、変更終了条件を満足すると、ステップS1009に進む。ステップS1009に進むと、炉別サイクル別投入熱量決定部502は、炉別サイクル別投入熱量Qininの変更終了条件を満足したときに得られている最新の炉別サイクル別投入熱量Qininを、炉別サイクル別投入熱量表示指示部505に出力する。そして、炉別サイクル別投入熱量表示指示部505は、この最新の炉別サイクル別投入熱量Qininを表示装置303図3を参照に表示するための表示データを生成して、入出力装置302を介して表示装置303に出力する。そして、図10のフローチャートによる処理を終了する。次に、図11のフローチャートを参照しながら、図10のステップS1004の熱風炉シミュレータ実行処理の詳細を説明する。尚、この図11の熱風炉シミュレータ実行処理は、後述する図12のフローチャートにおいても実行されるものである。まず、ステップS1101において、使用ガス容量算出部801aは、炉別サイクル別投入熱量Qininを取得する。次に、ステップS1102において、使用ガス容量算出部801aは、計算を行うサイクルnとして最初のサイクルを示す1を指定する。次に、ステップS1103において、使用ガス容量算出部801aは、BFG、COG、LDG、及び空気の各使用ガスの容量を算出する。次に、ステップS1104において、モデル式計算部801bは、モデル式1式~3式を使用して、サイクルnにおける燃焼開始から燃焼終了までの、熱風炉100a~100dにおける熱収支の計算を行う。この計算により、サイクルnの燃焼終了の時点での炉別サイクル別蓄熱煉瓦温度TCHin、炉別サイクル別炉壁煉瓦の構成煉瓦温度TWbin、炉別サイクル別排ガス温度Texin、炉別サイクル別燃焼効率ηinが得られる。次に、ステップS1105において、モデル式計算部801bは、サイクルnの燃焼終了の時点での炉別サイクル別蓄熱煉瓦温度TCHin、炉別サイクル別炉壁煉瓦の構成煉瓦温度TWbinを各炉運転内容指示部802aに出力する。また、モデル式計算部801bは、サイクルnの燃焼終了の時点での炉別サイクル別排ガス温度Texinを、これらの情報を要求した炉別残熱量算出部501又は炉別サイクル別投入熱量決定部502に出力する。次に、ステップS1106において、送風流量温度決定部802cは、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを取得する。次に、ステップS1107において、各炉運転内容指示部802aは、計算対象の時間tとしてt=1[秒]を設定する。次に、ステップS1108において、各炉運転内容指示部802aは、計算対象の時間t=1の炉別時間別通過送風流量BHSitと、炉別時間別混冷流量Airitの初期値を決定し、モデル式計算部802bに与える。次に、ステップS1109において、モデル式計算部802bは、モデル式1式~3式を使用して、送風期間における計算対象の時間tでの熱収支の計算を行う。この計算により、送風期間における計算対象の時間tでの炉別時間別ガス温度TGit、炉別サイクル別蓄熱煉瓦温度TCHinが得られる。次に、ステップS1110において、モデル式計算部802bは、送風期間が終了したか否かを判定する。尚、前述したように、本実施形態では、目標送風期間BTimemを固定値としているので、この送風期間も固定値である。この判定の結果、送風期間が終了していない場合には、ステップS1111に進む。ステップS1111に進むと、送風流量温度決定部802cは、計算対象の時間tに1[秒]を加算する。次に、ステップS1112において、送風流量温度決定部802cは、4式~6式により、目標送風温度BTm及び目標送風流量BVmが得られるように、炉別時間別送風流量BVi、炉別時間別通過送風流量BHSi、炉別時間別混冷流量Airiを算出決定する。そして、送風流量温度決定部802cは、各炉運転内容指示部802aを介して、計算対象の時間tの炉別時間別通過送風流量BHSitと、炉別時間別混冷流量Airitをモデル式計算部802bに与える。そして、送風期間が終了するまでの1秒毎の炉別時間別送風流量BVi、炉別時間別通過送風流量BHSi、炉別時間別混冷流量Airiが得られるまで、ステップS1109~S1112の処理を繰り返し行う。そして、前述したように、ステップS1110において、送風期間が終了したと判定すると、ステップS1113に進む。ステップS1113に進むと、送風流量温度決定部802cは、サイクルnが、ステップS1101で取得した炉別サイクル別投入熱量Qininで定まるサイクルNとなったか否かを判定する。この判定の結果、サイクルnがサイクルNとなっていない場合には、ステップS1114に進む。ステップS1114に進むと、使用ガス容量算出部801aは、計算を行うサイクルnに1を加算する。そして、サイクルNの送風期間が終了するまで、ステップS1103~S1114の処理を繰り返し行う。そして、サイクルnがサイクルNになると、ステップS1115に進む。ステップS1115に進むと、送風流量温度決定部802cは、送風終了時に得られた炉別時間別のガスの温度TGit、炉別時間別通過送風流量BHSit、炉別時間別混冷流量Airitを、4式~6式の右辺に与えて、炉別サイクル別送風温度BTkin及び炉別サイクル別送風流量BVkinを算出する。また、送風流量温度決定部802cは、送風終了時に得られた炉別サイクル別蓄熱煉瓦温度TCHinの中から珪石煉瓦下端部の温度を抽出し、抽出した温度を炉別サイクル別珪石煉瓦下端温度Tsiinとする。また、送風流量温度決定部802cは、7式、8式により、燃焼効率ηinを算出する。そして、送風流量温度決定部802cは、これらの炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、燃焼効率ηinを、これらの情報を要求した炉別残熱量算出部501又は炉別サイクル別投入熱量決定部502に出力する。これにより、図11のフローチャートによる処理が終了する。次に、図12のフローチャートを参照しながら、図10のステップS1002、S1007における炉別残熱量算出処理の一例を説明する。まず、ステップS1201において、第1の煉瓦下端温度判定部901は、珪石煉瓦温度計315a~315dにより測定された炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えているか否かを判定する。この判定の結果、炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えていない場合には、後述するステップS1210に進む。一方、炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えている場合には、ステップS1202に進む。ステップS1202に進むと、第2の煉瓦下端温度判定部902は、珪石煉瓦温度計315a~315dにより測定された炉別珪石煉瓦下端温度Tsiiの中に、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがあるか否かを判定する。この判定の結果、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがある場合には、当該炉別珪石煉瓦下端温度Tsiiを有する熱風炉iの残熱量Qrmiは0となる。よって、ステップS1203~S1206の処理を省略して後述するステップS1207に進む。一方、基準温度と等しい炉別珪石煉瓦下端温度Tsiiがない場合には、ステップS1203に進む。ステップS1203に進むと。第1の送風流量変化分決定部903は、例えば二分法に基づき、目標送風流量BVmの増加分ΔBVを決定する。そして、第1の送風流量変化分決定部903は、目標送風流量BVmに増加分ΔBVを加算したものを新たな目標送風流量BVmとし、これを、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。次に、ステップS1204において、熱風炉シミュレーション部503による熱風炉シミュレータ実行処理が行われる。この熱風炉シミュレータ実行処理により、ここでは、少なくとも、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガス温度Texinが得られる。尚、この熱風炉シミュレータ実行処理の詳細は、図11に示した通りである。次に、ステップS1205において、第1の残熱量判定部904は、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがあるか否かを判定する。この判定の結果、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしている熱風炉iがない場合には、当該熱風炉iが見つかるまでステップS1203~S1205の処理を繰り返し行う。そして、当該熱風炉iが見つかると、ステップS1206に進む。ステップS1206に進むと、該当炉残熱量算出部905は、ステップS1204で算出された最新の炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekinを用いて、当該熱風炉iの残熱量Qrmiを算出する。次に、ステップS1207において、残熱量記憶部906は、当該熱風炉iの残熱量Qrmiを記憶する。次に、ステップS1208において、第2の残熱量判定部907は、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たしたか否かを判定する。この判定の結果、当該条件を満たしている場合には、図12のフローチャートによる処理を終了する。一方、当該条件を満たしていない場合には、ステップS1209に進む。ステップS1209に進むと、第2の送風流量変化分決定部908は、ステップS1205において条件を満たしたと判定された熱風炉iに対してステップS1204で算出された最新の炉別サイクル別送風流量BVkinを固定値とするように熱風炉シミュレーション部503に指示する。そして、第2の送風流量変化分決定部908は、目標送風流量BVmの増減分ΔBVを決定する。そして、第2の送風流量変化分決定部908は、目標送風流量BVmに増加分ΔBVを加算又は減算したものを新たな目標送風流量BVmとし、これを、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。そして、ステップS1204に戻り、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しいという条件を満たすまで、ステップS1204~S1209の処理を繰り返し行う。前述したように、ステップS1201において、炉別珪石煉瓦下端温度Tsiiの全てが基準温度を超えていないと判定されると、ステップS1210に進む、ステップS1210に進むと、第3の送風流量変化分決定部909は、例えば二分法に基づき、目標送風流量BVmの減少分ΔBVを決定する。そして、第3の送風流量変化分決定部909は、目標送風流量BVmに減少分ΔBVを減算したものを新たな目標送風流量BVmとし、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風時間BTimemと共に、熱風炉シミュレーション部503に出力する。次に、ステップS1211において、熱風炉シミュレーション部503による熱風炉シミュレータ実行処理が行われる。この熱風炉シミュレータ実行処理により、ここでは、少なくとも、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガス温度Texinが得られる。尚、この熱風炉シミュレータ実行処理の詳細は、図11に示した通りである。次に、ステップS1212において、第3の残熱量判定部910は、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たしたか否かを判定する。この判定の結果、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たしていない場合には、ステップS1210に戻り、この条件を満たすまで、ステップS1210~S1213の処理を繰り返し行う。一方、全ての熱風炉100a~100dについて、炉別サイクル別送風流量BVkinが、目標送風流量BVmと等しく、且つ、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度以上であるという条件を満たしている場合には、前述したステップS1205に進む。[まとめ]以上のように本実施形態では、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを少なくとも入力として、熱風炉100a~100dにおけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texinを計算する熱風炉シミュレーション部503を構築する。炉別サイクル別投入熱量決定部502が、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを熱風炉シミュレーション部503に与えると、熱風炉シミュレーション部503は、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekin、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texinを算出する。評価値算出部504は、これらの情報を用いて、炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度よりも大きいという制約条件と、炉別サイクル別排ガスの温度Texinが基準温度よりも小さいという制約条件を満たす範囲で、減風減温により高炉の出銑量が減ったことに起因するコストの損失と、熱風炉100a~100dの状態を減風減温状態から目標の状態に変更するために使用するガスのコストと、が低いほど評価が高くなる評価関数の値評価値Yを求める。炉別サイクル別投入熱量決定部502は、評価値Yに基づいて、熱風炉シミュレーション部503に与えた炉別サイクル別投入熱量Qininが最適であると見なせるかどうかを判定し、最適であると見なせる炉別サイクル別投入熱量Qininを出力する。したがって、熱風炉100a~100dの状態を、減風減温の状態から、高炉の要求に応じた状態にするまでのサイクル数時間及び各サイクルにおける投入熱量を正確に把握することができる。また、本実施形態では、現在の炉別残熱量Qrmiを基準にした炉別サイクル別投入熱量Qininを決定するようにした。したがって、炉別サイクル別投入熱量Qininをより正確に決定することができる。ここで、現在の炉別残熱量Qrmiは、以下のようにして計算する。まず、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemのうち、目標送風流量BVmを変更して熱風炉シミュレーション部503を与える。次に、熱風炉シミュレーション部503から得られた炉別サイクル別送風流量BVkinがこの目標送風流量BVmと等しくなり、且つ、熱風炉シミュレーション部503から得られた炉別サイクル別珪石煉瓦下端温度Tsiinが基準温度と等しくなるという条件を満たす熱風炉iを見つけ出す。次に、当該条件を満たしたときの熱風炉iに対する炉別サイクル別送風温度BTkinと、炉別サイクル別送風流量BVkin、炉別サイクル別送風時間BTimekinを用いて、当該熱風炉iの残熱量Qrmiを算出する。そして、当該熱風炉iの送風流量を、当該条件を満たしたときの熱風炉iに対する炉別サイクル別送風流量BVkinに固定して、熱風炉シミュレーション部503による計算を、全ての熱風炉100a~100dの残熱量Qrmiが得られるまで繰り返し行う。[変形例1]以上のように、本実施形態では、現在の炉別残熱量Qrmiを計算する際に、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemのうち、目標送風流量BVmを変更するようにした。しかしながら、炉別サイクル別投入熱量Qinin、目標送風温度BTm、又は目標送風時間BTimemを変更するようにしてもよい。このようにした場合、まず、第1の送風流量変化分決定部903、第2の送風流量変化分決定部908、及び第3の送風流量変化分決定部909で決定するのは、目標送風温度BTm、目標送風時間BTimem、熱風炉シミュレーション部503に入力する炉別サイクル別投入熱量Qininの変化分となる。また、第1の残熱量判定部904、第2の残熱量判定部907、第3の残熱量判定部910では、炉別サイクル別送風流量BVkin炉別サイクル別送風時間BTimekin熱風炉シミュレーション部503から出力された炉別サイクル別投入熱量Qininと、目標送風温度BTm目標送風時間BTimem熱風炉シミュレーション部503に入力した炉別サイクル別投入熱量Qininとを比較することになる。また、炉別サイクル別投入熱量Qininを変更する場合、該当炉残熱量算出部905は、熱風炉シミュレーション部503で得られる燃焼効率ηinと、炉別サイクル別投入熱量Qininとを掛け合わせることにより、当該熱風炉iの残熱量Qrmiを算出する下記の12式も参照。また、炉別サイクル別投入熱量Qininではなく、炉別サイクル別蓄熱量Qbfinを変更するようにしてもよい。このようにした場合、炉別サイクル別蓄熱量Qbfinと炉別サイクル別投入熱量Qininと、燃焼効率ηinとが、以下の12式の関係にあることを利用すればよい。Qinin×ηin=Qbfin12すなわち、該当炉残熱量算出部905は、当該熱風炉iの炉別サイクル別蓄熱量Qbfinから、当該熱風炉iの残熱量Qrmiが直接的に求められる。[変形例2]炉別残熱量算出部501は、計算を開始する時刻における各熱風炉100a~100dの残熱量Qrmiを算出するようにして、残熱量Qrmiを繰り返し算出するようにした。しかしながら、必ずしもこのようにする必要はない。例えば、所定の時間隔で各熱風炉100a~100dの残熱量Qrmiを算出してもよい。[変形例3]評価関数評価値と制約条件制約式とは、前述したものに限定されない。すなわち、熱風炉100a~100dを減風減温の状態にすることによる不利益が小さくなるほど評価が高くなるものとしていれば、どのような評価関数評価値を用いてもよい。また、適切な操業ができなくなることを制約条件制約式としていれば、どのような制約条件制約式を用いてもよい。[変形例4]また、熱風炉シミュレーション部503が実行する熱風炉シミュレータは、炉別サイクル別投入熱量Qinin、目標送風温度BTm、目標送風流量BVm、目標送風時間BTimemを少なくとも入力として、熱風炉100a~100dにおけるガスの熱収支と煉瓦の熱収支の計算を行って、炉別サイクル別送風温度BTkin、炉別サイクル別送風流量BVkin、及び炉別サイクル別送風時間BTimekin好ましくは、これらに加えて、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texinを出力とするシミュレータであれば、前述したものに限定されない。また、炉別サイクル別珪石煉瓦下端温度Tsiin、及び炉別サイクル別排ガスの温度Texin以外の操業制約設備制約を採用してもよい。この場合、必ずしも熱風炉シミュレータにより操業制約設備制約を得る必要はない。[請求項との対応]<請求項1、6>熱風炉操業推算装置は、例えば、熱風炉制御計算機301を用いることにより実現される。熱風炉シミュレーション手段は、例えば、熱風炉シミュレーション部503を用いることにより実現される。熱風炉シミュレーション工程は、例えば、図11のフローチャートによる処理を行うことにより実現される。評価値算出手段は、例えば、評価値算出部504を用いることにより実現される。評価値算出工程は、例えば、図10のステップS1005の処理を行うことにより実現される。炉別サイクル別投入熱量決定手段は、例えば、炉別サイクル別投入熱量決定部502を用いることにより実現される。炉別サイクル別投入熱量決定工程は、例えば、図10のステップS1006、S1008の処理を行うことにより実現される。前記評価値算出手段工程により算出された評価値に基づいて、前記炉別サイクル別投入熱量の変更終了条件を満足したか否かを判断し、当該変更終了条件を満足したと判断するまで、前記炉別サイクル別投入熱量を繰り返し変更して、前記熱風炉シミュレーション手段工程による計算の入力とし、前記炉別サイクル別投入熱量の変更終了条件を満足したと判断したときの最新の前記炉別サイクル別投入熱量を出力することは、例えば、図10のS1004~S1009の処理により実現される。<請求項2、7>炉別残熱量算出手段は、例えば、炉別残熱量算出部501を用いることにより実現される。炉別残熱量算出工程は、例えば、図12のフローチャートによる処理を行うことにより実現される。判別工程は、例えば、図12のステップS1205の処理を行うことにより実現される。判定手段は、例えば、このような処理を行う炉別残熱量算出部501により実現される。算出工程は、例えば、図12のステップS1206の処理を行うことにより実現される。算出手段は、例えば、このような処理を行う炉別残熱量算出部501により実現される。前記判別した熱風炉については、前記条件を満たしたときに前記熱風炉シミュレーション手段により得られた値を固定値として、前記判別手段による判別と、前記算出手段による算出と、を、前記判別手段により全ての熱風炉が前記条件を満たすと判別されるまで繰り返し行うことは、例えば、ステップS1203~S1206の処理を行うことにより実現される。<請求項3、8>前記炉別残熱量を繰り返し算出することは、例えば、図10において、ステップS1007が繰り返し算出されることにより実現される。<請求項4、9>前記炉別サイクル別送風温度、前記炉別サイクル別送風流量、及び前記炉別サイクル別送風時間に加えて、炉別サイクル別珪石煉瓦下端温度、及び炉別サイクル別排ガス温度を少なくとも出力とする計算を行うことは、ステップS1104において、燃焼終了の時点での炉別サイクル別排ガス温度Texinが算出され、ステップS1111において、送風終了の時点での炉別時間別ガス温度TGitnが算出されることにより実現される。尚、以上説明した本発明の実施形態は、コンピュータがプログラムを実行することによって実現することができる。また、プログラムをコンピュータに供給するための手段、例えばかかるプログラムを記録したCD-ROM等のコンピュータ読み取り可能な記録媒体、又はかかるプログラムを伝送する伝送媒体も本発明の実施の形態として適用することができる。また、前記プログラムを記録したコンピュータ読み取り可能な記録媒体などのプログラムプロダクトも本発明の実施の形態として適用することができる。前記のプログラム、コンピュータ読み取り可能な記録媒体、伝送媒体及びプログラムプロダクトは、本発明の範疇に含まれる。また、以上説明した本発明の実施形態は、何れも本発明を実施するにあたっての具体化の例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。"