JP2023125249A - Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate - Google Patents

Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate Download PDF

Info

Publication number
JP2023125249A
JP2023125249A JP2022029241A JP2022029241A JP2023125249A JP 2023125249 A JP2023125249 A JP 2023125249A JP 2022029241 A JP2022029241 A JP 2022029241A JP 2022029241 A JP2022029241 A JP 2022029241A JP 2023125249 A JP2023125249 A JP 2023125249A
Authority
JP
Japan
Prior art keywords
temperature
slab
heating furnace
extraction
skid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2022029241A
Other languages
Japanese (ja)
Inventor
知義 小笠原
Tomoyoshi Ogasawara
行宏 松原
Yukihiro Matsubara
寛人 後藤
Hiroto Goto
光俊 ▲剣▼持
Mitsutoshi Kemmochi
豪巨 深江
Hideo Fukae
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JFE Steel Corp
Original Assignee
JFE Steel Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by JFE Steel Corp filed Critical JFE Steel Corp
Priority to JP2022029241A priority Critical patent/JP2023125249A/en
Publication of JP2023125249A publication Critical patent/JP2023125249A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Control Of Metal Rolling (AREA)
  • Heat Treatments In General, Especially Conveying And Cooling (AREA)
  • Control Of Heat Treatment Processes (AREA)

Abstract

To provide a correction device and a correction method, etc. for a slab temperature model capable of accurately predicting a temperature distribution including a lengthwise direction and a thickness direction in at least a width central part of a slab at extraction of a heating furnace.SOLUTION: A correction device 10 for a slab temperature model corrects, at a model correction part 14, a heat release amount Qskd at a heating furnace skid included in a thermal conduction equation composing a slab temperature model, to be used in an extraction temperature calculation part 12, so that a deviation e is below a threshold value ε, the deviation e between: a temperature difference Tcal_diff of a temperature Tcal_skid at a position in a lengthwise direction corresponding to a heating furnace skid of a slab S and a temperature Tcal at a position in a lengthwise direction corresponding to an interval between heating furnace skids, outputted from a crude temperature calculation part 13; and a temperature difference Tact_diff of a temperature Tact_skid at the position in the lengthwise direction corresponding to the heating furnace skid of the slab S and a temperature Tact at the position in the lengthwise direction corresponding to the interval between heating furnace skids, acquired at a crude temperature acquisition part 11.SELECTED DRAWING: Figure 2

Description

本発明は、スラブを加熱炉で加熱し、加熱炉から抽出されたスラブを粗圧延する際における、加熱炉の抽出時のスラブ温度の予測に用いられるスラブ温度モデルの修正装置、修正方法、加熱炉の炉温制御装置、炉温制御方法、及び鋼板の製造方法に関する。 The present invention provides a slab temperature model correction device, correction method, and heating method used for predicting the slab temperature at the time of extraction from the heating furnace when a slab is heated in a heating furnace and the slab extracted from the heating furnace is roughly rolled. The present invention relates to a furnace temperature control device, a furnace temperature control method, and a steel plate manufacturing method.

製鉄所において、製鋼工程において鋳造されたスラブは、一旦、スラブヤードに保管されて冷えた状態の冷片、もしくは、直送されるDHCR(Direct Hot Charge)の熱片として、熱延工程の加熱炉に装入され、所定の温度(約1200度)まで加熱される。そして、加熱炉から抽出されたスラブは、粗圧延工程にて粗圧延された後、仕上圧延工程で仕上圧延され、薄い鋼板となる。そして、鋼板は、ランナウト冷却工程を経て冷却された後、巻取工程にてコイルに巻き取られる。 At a steelworks, slabs cast during the steelmaking process are stored in a slab yard as cold pieces, or sent directly to the heating furnace of the hot rolling process as DHCR (Direct Hot Charge) hot pieces. and heated to a predetermined temperature (approximately 1200 degrees). Then, the slab extracted from the heating furnace is rough rolled in a rough rolling process, and then finished rolled in a finish rolling process to become a thin steel plate. After the steel plate is cooled through a runout cooling process, it is wound into a coil in a winding process.

ここで、ウォーキングビーム式加熱炉では、スラブは複数本(例えば8本)のスキッドビームで支えられる。スキッドビームは、固定式と移動式の2種類からなり、移動式スキッドビームが上昇、移動、及び下降を繰り返すことでスラブは抽出側へと搬送される。スラブは、幅方向、長さ方向、及び厚さ方向に所定の大きさを有する直方体形状をなし、スラブの幅方向がスラブの搬送方向となるように加熱炉内のスキッドビーム上に載置され、移動式スキッドビームが上昇、移動、及び下降を繰り返すことでスラブは抽出側へと搬送される。 Here, in a walking beam type heating furnace, the slab is supported by a plurality of (for example, eight) skid beams. There are two types of skid beams, a fixed type and a movable type, and the movable skid beam repeatedly ascends, moves, and descends to transport the slab to the extraction side. The slab has a rectangular parallelepiped shape with predetermined dimensions in the width direction, length direction, and thickness direction, and is placed on a skid beam in the heating furnace so that the width direction of the slab is the conveyance direction of the slab. , the slab is conveyed to the extraction side by the movable skid beam repeatedly rising, moving, and descending.

このとき、加熱中のスラブはスキッドビームを介して熱が奪われることにより、スキッドビームの接触箇所が低温度となり、スキッドマークと呼ばれるスラブの長さ方向の温度変動が生じる。スキッドマークは、スラブの幅方向にスラブを横断するように延びるものがスラブの長さ方向に沿って所定ピッチで形成される。なお、加熱炉の内部は高温環境であり、スラブの表面のスケールも発生するため、放射温度計でスラブの表面温度を精度良く測定することは困難であるため、加熱炉内で前述のスキッドマークを直接測定することはできず、一般的には粗圧延機の出側での温度計測によりスキッドマークの温度変動を確認することができる。 At this time, heat is removed from the slab being heated via the skid beam, resulting in a low temperature at the contact point of the skid beam, resulting in temperature fluctuations in the length direction of the slab called skid marks. The skid marks extend across the slab in the width direction of the slab and are formed at predetermined pitches along the length of the slab. The inside of the heating furnace is a high-temperature environment and scale on the surface of the slab is generated, so it is difficult to accurately measure the surface temperature of the slab with a radiation thermometer. It is not possible to directly measure the temperature of skid marks, and temperature fluctuations in skid marks can generally be confirmed by measuring the temperature at the exit side of the roughing mill.

加熱炉の操業では、熱延コイルの材質の観点から抽出温度を規定範囲に管理することが求められる。特に、昇温不足となりがちなスキッドビームの接触箇所(以下、スキッド部という)の温度を規定温度まで上昇させて抽出することが必須となる。ただし、抽出温度を過剰に上昇させると、燃料原単位(GJ/ton-steel)の悪化を招くため、適切に加熱操業を行う必要がある。 When operating a heating furnace, it is required to control the extraction temperature within a specified range from the viewpoint of the material of the hot-rolled coil. In particular, it is essential to raise the temperature of the contact area of the skid beam (hereinafter referred to as skid part), which tends to be insufficiently heated, to a specified temperature before extraction. However, if the extraction temperature is increased excessively, the fuel consumption rate (GJ/ton-steel) will deteriorate, so it is necessary to perform the heating operation appropriately.

従来、加熱炉の操業方法として、例えば、特許文献1及び特許文献2に示すものが知られている。
特許文献1に示す圧延加熱炉燃焼制御方法は、スキッド間についてはスラブの厚み方向一次元の熱伝導差分方程式を、スキッド部についてはスラブの長手方向の熱伝導を係数導入によって簡略化した二次元(スラブ厚み方向×スラブ長手方向)の熱伝導差分方程式を解くことで、炉内スラブのスキッド間とスキッド部の将来温度分布を予測し、スラブの焼き上げ品質を評価することにより設定炉温を算出するものである。
BACKGROUND ART Conventionally, methods for operating a heating furnace are known, for example, as shown in Patent Document 1 and Patent Document 2.
The rolling heating furnace combustion control method disclosed in Patent Document 1 uses a one-dimensional heat conduction difference equation in the thickness direction of the slab between the skids, and a two-dimensional heat conduction difference equation in the longitudinal direction of the slab that is simplified by introducing coefficients for the skid portion. By solving the heat conduction difference equation (slab thickness direction x slab longitudinal direction), the future temperature distribution between the skids and the skid part of the slab in the furnace is predicted, and the set furnace temperature is calculated by evaluating the baking quality of the slab. It is something to do.

特許文献1に示す圧延加熱炉燃焼制御方法によれば、少ない計算量で、炉内の全スラブに対してスキッド間とスキッド部の温度分布予測を行うことができる。 According to the rolling heating furnace combustion control method disclosed in Patent Document 1, the temperature distribution between skids and in the skid portion can be predicted for all slabs in the furnace with a small amount of calculation.

また、特許文献2に示す連続式加熱炉内のスラブ温度予測方法は、連続式加熱炉の「抽出口におけるスラブ温度測定値」と「炉温、在炉時間等から計算されるスラブ抽出温度」との誤差をスラブ抽出毎に測定し、これと別途スラブ温度計算により求めた影響係数に対し逐次最小2乗推定を行って炉内複数ゾーンの総括熱吸収率をオンラインで逐次推定すると共に、この推定値に基づいて各ゾーンの総括熱吸収率を更新することで炉内スラブの温度予測精度を向上させるものである。 In addition, the slab temperature prediction method in a continuous heating furnace shown in Patent Document 2 is based on the "slab temperature measurement value at the extraction port" of the continuous heating furnace and the "slab extraction temperature calculated from the furnace temperature, furnace time, etc." The error is measured for each slab extraction, and the overall heat absorption rate of multiple zones in the furnace is sequentially estimated online by sequentially performing least squares estimation on the influence coefficient obtained by separately calculating the slab temperature. By updating the overall heat absorption rate of each zone based on the estimated value, the temperature prediction accuracy of the furnace slab is improved.

特開平7-258752号公報Japanese Patent Application Publication No. 7-258752 特開平6-264153号公報Japanese Patent Application Publication No. 6-264153

しかしながら、特許文献1に示す圧延加熱炉燃焼制御方法及び特許文献2に示す連続式加熱炉内のスラブ温度予測方法にあっては、以下の問題点があった。
即ち、特許文献1に示す圧延加熱炉燃焼制御方法の場合、スキッド間についてはスラブの厚み方向一次元の熱伝導差分方程式を、スキッド部についてはスラブの長手方向の熱伝導を係数導入によって簡略化した二次元の熱伝導差分方程式を解くことで、炉内スラブのスキッド間とスキッド部の将来温度分布を予測しているが、その予測された炉内スラブのスキッド間とスキッド部の将来温度分布が真に正しいか否かの検証さなされていない。従って、実際の炉内スラブのスキッド間とスキッド部の温度分布が、予測された炉内スラブのスキッド間とスキッド部の温度分布と乖離している可能性があり、スラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度の予測値の精度を担保することができない。
However, the rolling heating furnace combustion control method shown in Patent Document 1 and the slab temperature prediction method in a continuous heating furnace shown in Patent Document 2 have the following problems.
That is, in the case of the rolling heating furnace combustion control method shown in Patent Document 1, the one-dimensional heat conduction difference equation in the thickness direction of the slab is simplified between the skids, and the heat conduction in the longitudinal direction of the slab is simplified by introducing coefficients for the skid part. By solving the two-dimensional heat conduction difference equation, the future temperature distribution between the skids and the skid part of the furnace slab is predicted. It has not been verified whether or not it is truly correct. Therefore, the actual temperature distribution between the skids and the skid part of the slab in the furnace may deviate from the predicted temperature distribution between the skids and the skid part of the slab in the furnace. It is not possible to ensure the accuracy of the predicted value of the temperature at least in the width center of the slab during extraction in the furnace, including the length and thickness directions.

また、特許文献2に示す連続式加熱炉内のスラブ温度予測方法の場合、連続式加熱炉の「抽出口におけるスラブ温度測定値」と「炉温、在炉時間等から計算されるスラブ抽出温度」との誤差をスラブ抽出毎に測定しているが、一般的に加熱炉の抽出直後だとスラブ表面のスケールの影響で温度測定自体が困難であり、連続式加熱炉の適切な「抽出口におけるスラブ温度測定値」を得ることができない。このため、スラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度の予測値の精度を担保することができない。 In addition, in the case of the slab temperature prediction method in a continuous heating furnace shown in Patent Document 2, the slab temperature measurement value at the extraction port of the continuous heating furnace and the slab extraction temperature calculated from the furnace temperature, furnace time, etc. ” error is measured every time a slab is extracted, but it is generally difficult to measure the temperature immediately after extraction from the heating furnace due to the influence of scale on the slab surface. "slab temperature measurements" cannot be obtained. For this reason, it is not possible to ensure the accuracy of the predicted value of the temperature including the length direction and thickness direction at least in the width center portion of the slab at the time of extraction in the heating furnace using the slab temperature model.

従って、本発明はこれら従来の問題点を解決するためになされたものであり、その目的は、スラブ温度モデルを修正することで、加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を精度よく予測することができる、スラブ温度モデルの修正装置、修正方法、加熱炉の炉温制御装置、炉温制御方法、及び鋼板の製造方法を提供することにある。 Therefore, the present invention has been made to solve these conventional problems, and its purpose is to modify the slab temperature model so that the temperature at least in the longitudinal direction of the slab during extraction in the heating furnace is To provide a slab temperature model correction device and correction method, a heating furnace furnace temperature control device, a furnace temperature control method, and a steel plate manufacturing method, which can accurately predict the temperature distribution including the thickness direction. be.

上記課題を解決するために、本発明の一態様に係るスラブ温度モデルの修正装置は、スラブを加熱炉で加熱し、該加熱炉から抽出された前記スラブを粗圧延機で粗圧延する際における、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布の予測に用いられるスラブ温度モデルの修正装置であって、前記粗圧延機の出側に設置された粗温度測定計によって測定した粗圧延後の前記スラブの裏面の幅中央部における長さ方向の温度分布を取得する粗温度取得部と、前記加熱炉に装入される前記スラブの装入温度を初期値として、前記スラブの加熱炉在炉中の操業条件を用いて、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、前記スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する抽出温度算出部と、該抽出温度算出部から出力された、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を初期値として、前記加熱炉から抽出された前記スラブを粗圧延する際の操業条件を用いて、粗圧延後の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、熱伝導方程式に基づき算出し、出力する粗出側温度算出部と、該粗出側温度算出部から出力された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗出側温度算出部から出力された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差と、前記粗温度取得部で取得された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗温度取得部で取得された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差、の偏差の大きさが閾値未満となるように、前記抽出温度算出部で用いられる、前記スラブ温度モデルを構成する熱伝導方程式に含まれる加熱炉スキッド部での抜熱量を修正するモデル修正部とを備えていることを要旨とする。 In order to solve the above problems, a slab temperature model correction device according to one aspect of the present invention heats a slab in a heating furnace, and when rough rolling the slab extracted from the heating furnace in a rough rolling mill. , a slab temperature model correction device used for predicting the temperature distribution including the length direction and thickness direction in at least the width center portion of the slab during extraction in the heating furnace, the device comprising: a rough temperature acquisition unit that acquires the temperature distribution in the longitudinal direction at the width center portion of the back surface of the slab after rough rolling measured by an installed rough temperature measuring meter; Using the input temperature as an initial value, and using the operating conditions while the slab is in the heating furnace, calculate the temperature distribution including the length direction and thickness direction at least at the width center of the slab at the time of extraction from the heating furnace. , an extraction temperature calculation unit that calculates and outputs the calculation based on a heat conduction equation constituting the slab temperature model; and a length of the slab at least in the width center portion at the time of extraction of the heating furnace, which is output from the extraction temperature calculation unit. Using the temperature distribution including the width direction and the thickness direction as an initial value, and using the operating conditions when rough rolling the slab extracted from the heating furnace, the length of the slab after rough rolling at least at the center width A rough-out side temperature calculation unit that calculates and outputs the temperature distribution including the width direction and the thickness direction based on a heat conduction equation, and the width center of the back side of the slab output from the rough-out side temperature calculation unit. The temperature at the longitudinal position corresponding to the heating furnace skid part in the heating furnace skid part, and the temperature at the longitudinal position corresponding to the heating furnace skid part at the width center part of the back surface of the slab, which is output from the raw side temperature calculation part. temperature, the temperature at a longitudinal position corresponding to the heating furnace skid part at the width center of the back surface of the slab, acquired by the rough temperature acquisition unit, and the temperature at the longitudinal position corresponding to the heating furnace skid part, acquired by the rough temperature acquisition unit. , used in the extraction temperature calculation unit so that the magnitude of the deviation of the temperature difference between the width center part of the back surface of the slab and the temperature at the longitudinal position corresponding between the heating furnace skids is less than a threshold value, The present invention further comprises a model correction section for correcting the amount of heat removed in the heating furnace skid section included in the heat conduction equation constituting the slab temperature model.

また、本発明の別の態様に係る加熱炉の炉温制御装置は、前述のスラブ温度モデルの修正装置によって修正されたスラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を予測するスラブ温度予測部と、該スラブ温度予測部で予測された加熱炉の抽出時のスラブの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、前記加熱炉内における炉温を制御する炉温制御部とを備えていることを要旨とする。 Further, a furnace temperature control device for a heating furnace according to another aspect of the present invention uses a slab temperature model corrected by the above-mentioned slab temperature model correction device to control the temperature at least in the width center portion of the slab during extraction of the heating furnace. A slab temperature prediction unit that predicts temperature distribution including the length direction and thickness direction, and a length corresponding to the heating furnace skid portion at the width center of the slab at the time of extraction of the heating furnace predicted by the slab temperature prediction unit The present invention further comprises a furnace temperature control section that controls the furnace temperature in the heating furnace so that the temperature in the thickness direction at the direction position exceeds the minimum extraction temperature.

また、本発明の別の態様に係るスラブ温度モデルの修正方法は、スラブを加熱炉で加熱し、該加熱炉から抽出された前記スラブを粗圧延機で粗圧延する際における、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布の予測に用いられるスラブ温度モデルの修正方法であって、前記粗圧延機の出側に設置された粗温度測定計によって測定した粗圧延後の前記スラブの裏面の幅中央部における長さ方向の温度分布を取得する粗温度取得ステップと、前記加熱炉に装入される前記スラブの装入温度を初期値として、前記スラブの加熱炉在炉中の操業条件を用いて、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、前記スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する抽出温度算出ステップと、該抽出温度算出ステップで出力された、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を初期値として、前記加熱炉から抽出された前記スラブを粗圧延する際の操業条件を用いて、粗圧延後の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、熱伝導方程式に基づき算出し、出力する粗出側温度算出ステップと、該粗出側温度算出ステップで出力された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗出側温度算出ステップで出力された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差と、前記粗温度取得ステップで取得された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗温度取得ステップで取得された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差、の偏差の大きさが閾値未満となるように、前記抽出温度算出ステップで用いられる、前記スラブ温度モデルを構成する熱伝導方程式に含まれる加熱炉スキッド部での抜熱量を修正するモデル修正ステップとを含むことを要旨とする。 Further, in a method for correcting a slab temperature model according to another aspect of the present invention, a slab is heated in a heating furnace, and the slab extracted from the heating furnace is roughly rolled in a rough rolling mill. A method for correcting a slab temperature model used for predicting temperature distribution including length direction and thickness direction in at least the center width of the slab during extraction, the method comprising: a rough temperature model installed at the exit side of the rough rolling mill; a rough temperature acquisition step of acquiring the temperature distribution in the length direction at the width center of the back surface of the slab after rough rolling measured by a measuring meter; and setting the charging temperature of the slab charged to the heating furnace to an initial value. Using the operating conditions of the slab while it is in the heating furnace, the temperature distribution including the length direction and thickness direction at least in the width center of the slab at the time of extraction in the heating furnace is determined by the slab temperature model. an extraction temperature calculation step that calculates and outputs the extraction temperature based on a heat conduction equation that constitutes the heating furnace; Using the operating conditions when rough rolling the slab extracted from the heating furnace, with the temperature distribution including the direction as an initial value, the length direction and thickness of the slab after rough rolling at least at the center width thereof. a rough output side temperature calculation step of calculating and outputting the temperature distribution including the direction based on a heat conduction equation; The temperature difference between the temperature at the longitudinal position corresponding to the part and the temperature at the longitudinal position corresponding to between the heating furnace skids at the width center part of the back surface of the slab, which is output in the raw side temperature calculation step. , the temperature at the longitudinal position corresponding to the heating furnace skid part in the width center of the back surface of the slab, obtained in the rough temperature obtaining step, and the back surface of the slab, obtained in the rough temperature obtaining step. The slab temperature model used in the extraction temperature calculation step is adjusted such that the temperature difference between the temperature at the longitudinal position corresponding to between the heating furnace skids at the width center of the slab is less than a threshold value. The gist of the present invention is to include a model modification step of modifying the amount of heat removed from the heating furnace skid part included in the heat conduction equation.

また、本発明の別の態様に係る加熱炉の炉温制御方法は、前述のスラブ温度モデルの修正方法によって修正されたスラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を予測するスラブ温度予測ステップと、該スラブ温度予測ステップで予測された加熱炉の抽出時のスラブの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、前記加熱炉内における炉温を制御する炉温制御ステップとを含むことを要旨とする。 Further, a furnace temperature control method for a heating furnace according to another aspect of the present invention uses a slab temperature model corrected by the above-described slab temperature model correction method to A slab temperature prediction step that predicts the temperature distribution including the length direction and the thickness direction, and a length corresponding to the heating furnace skid portion at the width center of the slab at the time of extraction of the heating furnace predicted in the slab temperature prediction step. The method further comprises a furnace temperature control step of controlling the furnace temperature in the heating furnace so that the temperature in the thickness direction at the direction position exceeds the minimum extraction temperature.

また、本発明の別の態様に係る鋼板の製造方法は、前述の加熱炉の炉温制御方法によって加熱炉内の炉温を制御する炉温制御工程を含むことを要旨とする。 Moreover, the gist of the method for manufacturing a steel plate according to another aspect of the present invention includes a furnace temperature control step of controlling the furnace temperature in the heating furnace by the above-mentioned furnace temperature control method for the heating furnace.

本発明に係るスラブ温度モデルの修正装置、修正方法、加熱炉の炉温制御装置、炉温制御方法、及び鋼板の製造方法によれば、スラブ温度モデルを修正することで、加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を精度よく予測することができる、スラブ温度モデルの修正装置、修正方法、加熱炉の炉温制御装置、炉温制御方法、及び鋼板の製造方法を提供できる。 According to the slab temperature model correction device, correction method, heating furnace furnace temperature control device, furnace temperature control method, and steel plate manufacturing method according to the present invention, by correcting the slab temperature model, A slab temperature model correction device, correction method, furnace temperature control device for a heating furnace, and furnace temperature control method that can accurately predict the temperature distribution including the length direction and thickness direction in at least the central width of a slab. , and a method for manufacturing a steel plate.

本発明の一実施形態に係るスラブ温度モデルの修正装置及び加熱炉の炉温制御装置が適用される熱間圧延ラインの概略構成図である。1 is a schematic configuration diagram of a hot rolling line to which a slab temperature model correction device and a heating furnace furnace temperature control device according to an embodiment of the present invention are applied. 図1におけるスラブ温度モデルの修正装置の内部構成を示すブロック図である。FIG. 2 is a block diagram showing the internal configuration of the slab temperature model correction device in FIG. 1. FIG. 図1における加熱炉の炉温制御装置の内部構成を示すブロック図である。FIG. 2 is a block diagram showing the internal configuration of the furnace temperature control device of the heating furnace in FIG. 1. FIG. 粗温度測定計によって測定した粗圧延後のスラブの裏面の幅中央部における長さ方向の温度分布の一例を示すグラフである。It is a graph which shows an example of the temperature distribution of the length direction in the width center part of the back surface of the slab after rough rolling measured by the rough temperature meter. 粗出側温度算出部から出力された、粗圧延後のスラブの裏面の幅中央部における長さ方向の温度分布の一例を示すグラフである。It is a graph which shows an example of the temperature distribution in the length direction in the width center part of the back surface of the slab after rough rolling, which is output from the rough rolling side temperature calculation unit. スラブ温度モデルの修正装置における処理の流れを説明するためのフローチャートである。It is a flowchart for explaining the flow of processing in the slab temperature model correction device. 図6に示すフローチャートにおけるステップS4(モデル修正ステップ)における処理の流れの詳細を示すフローチャートである。7 is a flowchart showing details of the process flow in step S4 (model modification step) in the flowchart shown in FIG. 6. FIG. 加熱炉の炉温制御装置における処理の流れを説明するためのフローチャートである。It is a flowchart for explaining the flow of processing in the furnace temperature control device of the heating furnace. スラブの形状を説明するための図である。It is a figure for explaining the shape of a slab. 実施例において、粗圧延後のスラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置近傍の温度分布を示すグラフである。2 is a graph showing the temperature distribution in the vicinity of the longitudinal position corresponding to the heating furnace skid portion at the width center portion of the back surface of the slab after rough rolling in an example. 実施例において、スラブ温度モデル(補正係数η)を修正する前と修正した後の、粗出側温度算出部から出力された、粗圧延後のスラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置近傍の温度分布を比較して示すグラフである。In the example, the heating furnace skid portion at the width center of the back side of the slab after rough rolling, which is output from the rough-rolling side temperature calculation unit before and after the slab temperature model (correction coefficient η) is corrected. It is a graph showing a comparison of temperature distributions in the vicinity of corresponding longitudinal positions.

以下、本発明の実施の形態を図面を参照して説明する。以下に示す実施形態は、本発明の技術的思想を具体化するための装置や方法を例示するものであって、本発明の技術的思想は、構成部品の材質、形状、構造、配置等を下記の実施形態に特定するものではない。また、図面は模式的なものである。そのため、厚みと平面寸法との関係、比率等は現実のものとは異なることに留意すべきであり、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれている。 Embodiments of the present invention will be described below with reference to the drawings. The embodiments shown below illustrate devices and methods for embodying the technical idea of the present invention. It is not limited to the embodiments described below. Furthermore, the drawings are schematic. Therefore, it should be noted that the relationships, ratios, etc. between thickness and planar dimensions are different from those in reality, and the drawings also include portions where the relationships and ratios of dimensions are different.

図1には、本発明の一実施形態に係るスラブ温度モデルの修正装置及び加熱炉の炉温制御装置が適用される熱間圧延ラインの概略構成が示されている。
製鉄所において、製鋼工程において鋳造されたスラブSは、一旦、スラブヤードに保管されて冷えた状態の冷片、もしくは、直送されるDHCR(Direct Hot Charge)の熱片として、図1に示すように、熱間圧延ライン1における、加熱炉2に装入され、所定の温度(約1200度)まで加熱される。そして、加熱炉2内を搬送されて抽出されたスラブSは、複数台(本実施形態ではn台)の粗圧延機4~4によって粗圧延された後、図示しない仕上圧延機によって仕上圧延され、薄い鋼板となる。そして、この鋼板は、冷却設備5にて冷却された後、巻取り装置6によってコイル状に巻き取られる。
FIG. 1 shows a schematic configuration of a hot rolling line to which a slab temperature model correction device and a heating furnace furnace temperature control device according to an embodiment of the present invention are applied.
At a steelworks, the slab S cast during the steelmaking process is temporarily stored in a slab yard as a cold piece, or directly delivered as a hot piece of DHCR (Direct Hot Charge), as shown in Figure 1. Then, it is charged into a heating furnace 2 in a hot rolling line 1 and heated to a predetermined temperature (approximately 1200 degrees). Then, the slab S transported through the heating furnace 2 and extracted is rough rolled by a plurality of (n units in this embodiment) rough rolling mills 4 1 to 4 n , and then finished by a finishing mill (not shown). It is rolled into a thin steel plate. After this steel plate is cooled in the cooling equipment 5, it is wound up into a coil shape by the winding device 6.

この加熱炉2は、スラブSの搬送方向の上流側から下流側に向けて複数(本実施形態ではn個)の燃焼制御帯2~2を備えている。そして、各燃焼制御帯2~2には、スラブSを加熱するための燃焼バーナー3が設置されている。各燃焼制御帯2~2の燃焼バーナー3は、燃料流量制御部23に接続され、それぞれの燃料流量制御部23は、加熱炉の炉温制御装置20に接続されている。 The heating furnace 2 includes a plurality of (n in this embodiment) combustion control zones 2 1 to 2 n extending from the upstream side to the downstream side in the conveyance direction of the slab S. A combustion burner 3 for heating the slab S is installed in each combustion control zone 2 1 to 2 n . The combustion burners 3 of each combustion control zone 2 1 to 2 n are connected to a fuel flow rate control section 23, and each fuel flow rate control section 23 is connected to a furnace temperature control device 20 of the heating furnace.

ここで、加熱炉2は、ウォーキングビーム式加熱炉であり、スラブSは複数本(例えば8本)のスキッドビームで支えられる。スキッドビームは、固定式と移動式の2種類からなり、移動式スキッドビームが上昇、移動、及び下降を繰り返すことでスラブSは抽出側へと搬送される。スラブSは、図9に示すように、矢印xで示す長さ方向、矢印yで示す幅方向、及び矢印zで示す厚さ方向のそれぞれにおいて、幅w、長さl、及び厚さtを有する直方体形状をなしている。加熱炉2内においては、スラブSは、その表面Saが上側、裏面Sbが下側となり、図1に示すように、スラブSの幅方向がスラブSの搬送方向となるように加熱炉2内のスキッドビーム上に載置され、移動式スキッドビームが上昇、移動、及び下降を繰り返すことでスラブSは抽出側へと搬送される。 Here, the heating furnace 2 is a walking beam type heating furnace, and the slab S is supported by a plurality of (for example, eight) skid beams. There are two types of skid beams, a fixed type and a movable type, and the slab S is conveyed to the extraction side by repeating rising, moving, and lowering of the movable skid beam. As shown in FIG. 9, the slab S has a width w, a length l, and a thickness t in the length direction indicated by an arrow x, the width direction indicated by an arrow y, and the thickness direction indicated by an arrow z. It has a rectangular parallelepiped shape. Inside the heating furnace 2, the slab S is placed so that its front surface Sa is on the upper side and its back surface Sb is on the lower side, and the width direction of the slab S is the conveying direction of the slab S, as shown in FIG. The slab S is placed on a skid beam, and the movable skid beam repeats raising, moving, and lowering, thereby conveying the slab S to the extraction side.

加熱炉2中をスラブSが搬送される際には、スラブSはスキッドビームを介して熱が奪われることにより、スキッドビームの接触箇所が低温度となり、スキッドマークと呼ばれるスラブSの長さ方向の温度変動が生じる。スキッドマークは、スラブSの幅方向にスラブを横断するように延びるものがスラブの長さ方向に沿って所定ピッチで形成される。
加熱炉2における操業では、熱延コイルの材質の観点からスラブSの抽出温度を規定範囲に管理することが求められる。特に、昇温不足となりがちなスキッドビームの接触箇所(以下、加熱炉スキッド部という)の温度を規定温度まで上昇させて抽出することが必須となる。
When the slab S is conveyed through the heating furnace 2, heat is removed from the slab S through the skid beam, so that the contact point of the skid beam becomes low temperature, and the lengthwise direction of the slab S, called a skid mark, is Temperature fluctuations occur. The skid marks extend across the slab in the width direction of the slab S and are formed at predetermined pitches along the length of the slab.
In operation in the heating furnace 2, it is required to control the extraction temperature of the slab S within a specified range from the viewpoint of the material of the hot-rolled coil. In particular, it is essential to raise the temperature of the contact point of the skid beam (hereinafter referred to as the heating furnace skid section), which tends to be insufficiently heated, to a specified temperature for extraction.

このため、本実施形態においては、加熱炉2の装入時のスラブSの温度を測定する装入温度測定計7により加熱炉2に装入されるスラブSの装入温度を測定し、加熱炉の炉温制御装置20において、この加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件(操業条件としては、在炉時間、炉温、スキッド接触温度などが挙げられる)を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向と厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式に基づき算出し、この算出された加熱炉抽出時のスラブSの幅中央部における長さ方向と厚さ方向の温度分布から得られる、加熱炉2の抽出時のスラブSの加熱スキッド部における厚さ方向の温度が最低抽出温度(例えば、1200℃)を上回るように、各燃料流量制御部23を制御して各燃焼バーナー3への燃料流量を制御し、加熱炉2内の炉温を制御するようにしている。 Therefore, in this embodiment, the charging temperature of the slab S charged into the heating furnace 2 is measured by the charging temperature meter 7 that measures the temperature of the slab S when charging the heating furnace 2, and the charging temperature of the slab S charged into the heating furnace 2 is measured. In the furnace temperature control device 20 of the furnace, the charging temperature of the slab S charged into the heating furnace 2 is set as an initial value, and the operating conditions (operating conditions include the furnace time, The temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 is calculated using the heat conduction equation that constitutes the slab temperature model. The thickness at the heating skid part of the slab S at the time of extraction in the heating furnace 2, obtained from the calculated temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2. Each fuel flow control unit 23 is controlled to control the fuel flow rate to each combustion burner 3 so that the temperature in the horizontal direction exceeds the minimum extraction temperature (for example, 1200° C.), and the furnace temperature in the heating furnace 2 is controlled. I try to do that.

ここで、スラブ温度モデルを構成する熱伝導方程式は、後述する抽出温度算出部12で用いられる(1)式~(7)式と同一である。
しかしながら、加熱炉の炉温制御装置20において、加熱炉2の抽出時のスラブSの幅中央部の長さ方向と厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式に基づき算出した予測値は、実測値と乖離しており、加熱炉2の抽出時のスラブSの幅中央部の長さ方向と厚さ方向の温度分布を精度よく予測することができないという問題がある。
Here, the heat conduction equations constituting the slab temperature model are the same as equations (1) to (7) used in the extraction temperature calculating section 12, which will be described later.
However, in the furnace temperature control device 20 of the heating furnace, the temperature distribution in the length direction and thickness direction of the central width of the slab S at the time of extraction in the heating furnace 2 is calculated based on the heat conduction equation that constitutes the slab temperature model. The predicted value deviates from the actual value, and there is a problem that the temperature distribution in the length direction and thickness direction of the width center portion of the slab S at the time of extraction in the heating furnace 2 cannot be accurately predicted.

従って、本実施形態にあっては、図1に示すように、加熱炉の炉温制御装置20にスラブ温度モデルの修正装置10を接続し、この修正装置10において、炉温制御装置20で用いられるスラブ温度モデルを修正することで、加熱炉の抽出時のスラブの幅中央部における長さ方向及び厚さ方向の温度分布を精度よく予測することを可能としている。 Therefore, in this embodiment, as shown in FIG. By modifying the slab temperature model, it is possible to accurately predict the temperature distribution in the longitudinal and thickness directions at the central width of the slab during extraction in the heating furnace.

このスラブ温度モデルの修正装置10は、図1に示すように、加熱炉の炉温制御装置20に接続されている。また、スラブ温度モデルの修正装置10は、加熱炉2の入側に設置されて加熱炉2に装入されるスラブSの表面Sa及び裏面Sbの装入温度を測定する装入温度測定計7と、粗圧延機4~4の出側に設置されて粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を測定する粗温度測定計8とに接続されている。装入温度測定計7及び粗温度測定計8は、放射温度計により構成されている。 As shown in FIG. 1, this slab temperature model correction device 10 is connected to a furnace temperature control device 20 of a heating furnace. The slab temperature model correction device 10 also includes a charging temperature measuring meter 7 installed at the entrance side of the heating furnace 2 and measuring the charging temperature of the front surface Sa and the back surface Sb of the slab S charged into the heating furnace 2. and a rough temperature measuring meter 8 which is installed on the exit side of the rough rolling mills 4 1 to 4 n and measures the temperature distribution in the length direction at the width center of the back surface Sb of the slab S after rough rolling. There is. The charging temperature meter 7 and the crude temperature meter 8 are constituted by radiation thermometers.

スラブ温度モデルの修正装置10の内部構成は、図2に示されており、修正装置10は、粗温度測定計8に接続された粗温度取得部11と、装入温度測定計7に接続された抽出温度算出部12と、粗出側温度算出部13と、モデル修正部14と、出力部15とを備えている。 The internal configuration of the slab temperature model correction device 10 is shown in FIG. It includes an extraction temperature calculation section 12, a crude extraction temperature calculation section 13, a model correction section 14, and an output section 15.

スラブ温度モデルの修正装置10は、演算処理機能を有するコンピュータシステムであり、ハードウェアに予め記憶された各種専用のコンピュータプログラムの命令により、粗温度取得部11(後に述べるステップS1)、抽出温度算出部12(ステップS2)、粗出側温度算出部13(ステップS3)、モデル修正部14(ステップS4)、及び出力部15(ステップS5)の各機能をソフトウェア上で実行する。 The slab temperature model correction device 10 is a computer system having an arithmetic processing function, and the rough temperature acquisition unit 11 (step S1 described later) performs extraction temperature calculation based on instructions from various dedicated computer programs stored in advance in the hardware. The functions of the section 12 (step S2), the raw temperature calculation section 13 (step S3), the model correction section 14 (step S4), and the output section 15 (step S5) are executed on software.

修正装置10の粗温度取得部11は、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を取得する。粗温度測定計8では、粗圧延機4~4の出側でスラブSの裏面Sbの温度を定期的に測定し、その結果が図示しない記録部に記録され、粗温度取得部11がその結果を取得する。 The rough temperature acquisition unit 11 of the correction device 10 calculates the length at the width center of the back surface Sb of the slab S after rough rolling, which is measured by the rough temperature measuring meter 8 installed on the exit side of the rough rolling mills 4 1 to 4 n . Obtain the temperature distribution in the direction. The rough temperature measuring meter 8 periodically measures the temperature of the back surface Sb of the slab S on the outlet side of the rough rolling mills 4 1 to 4 n , and the results are recorded in a recording section (not shown). Get the result.

粗温度測定計8によって測定した粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布の一例が図4に示されている。
また、抽出温度算出部12は、加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件(操業条件としては、在炉時間、炉温、スキッド接触温度などが挙げられる)を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する。
An example of the temperature distribution in the length direction at the width center of the back surface Sb of the slab S after rough rolling measured by the rough temperature measuring meter 8 is shown in FIG.
In addition, the extraction temperature calculation unit 12 sets the charging temperature of the slab S to be charged into the heating furnace 2 as an initial value, and calculates operating conditions for the slab S while it is in the heating furnace (operating conditions include time in the furnace, temperature, skid contact temperature, etc.), the temperature distribution in the length direction and thickness direction at the center of the width of the slab S at the time of extraction in the heating furnace 2 is expressed in the heat conduction equation that constitutes the slab temperature model. Calculate and output based on

加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度が装入温度測定計7により測定され、この測定された加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度が抽出温度算出部12に入力される。そして、抽出温度算出部12は、このスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の前述した操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する次の(1)式~(7)式で表される熱伝導方程式に基づき算出し、出力する。 The charging temperature of the front surface Sa and back surface Sb of the slab S when charging into the heating furnace 2 is measured by the charging temperature measuring meter 7, and the charging temperature of the front surface Sa and back surface of the slab S when charging into the heating furnace 2 is measured. The Sb charging temperature is input to the extraction temperature calculation section 12. Then, the extraction temperature calculation unit 12 uses the charging temperature of the slab S as an initial value and the width center of the slab S at the time of extraction in the heating furnace 2, using the above-mentioned operating conditions while the slab S is in the heating furnace. The temperature distribution in the length direction and thickness direction in the section is calculated and output based on the heat conduction equation expressed by the following equations (1) to (7) that constitute the slab temperature model.

Figure 2023125249000002
Figure 2023125249000002

ここで、ρ:比重[kg/]、C:比熱[kcal/(kg・K)]、λ:熱伝導率[kcal/(m・hr・K)]、θ:温度[℃]、x:長さ方向座標[m]、y:厚さ方向座標[m]である。 Here, ρ: specific gravity [kg/ 3 ], C: specific heat [kcal/(kg・K)], λ: thermal conductivity [kcal/(m・hr・K)], θ: temperature [℃], x : Length direction coordinate [m], y: Thickness direction coordinate [m].

(1)式の熱伝導方程式を用いての温度分布の算出では、時間方向と空間方向に差分化してスラブSの表面Sa及び裏面Sbでの境界条件を設定することで、スラブ温度を時間刻み毎に抽出時まで算出して、出力する。これにより、加熱炉スキッド部により生じるスラブSの長さ方向の温度変動が得られる。
この熱伝導方程式に設定される境界条件は、加熱炉2からのスラブSへの輻射入熱と加熱炉スキッド部での抜熱量とからなる。
表面(上面)Saからの輻射入熱は、次の(2)式及び(3)式により表される。
In calculating the temperature distribution using the heat conduction equation of equation (1), the slab temperature is calculated in time increments by setting the boundary conditions on the front surface Sa and the back surface Sb of the slab S by differentiating them in the time and space directions. Calculates up to the time of extraction each time and outputs it. Thereby, temperature fluctuations in the length direction of the slab S caused by the heating furnace skid portion can be obtained.
The boundary conditions set in this heat conduction equation include the radiant heat input from the heating furnace 2 to the slab S and the amount of heat removed at the heating furnace skid portion.
The radiant heat input from the surface (upper surface) Sa is expressed by the following equations (2) and (3).

Figure 2023125249000003
Figure 2023125249000003

また、裏面(下面)Sbからの輻射入熱は、次の(4)式及び(5)式で表される。 Further, the radiant heat input from the back surface (lower surface) Sb is expressed by the following equations (4) and (5).

Figure 2023125249000004
Figure 2023125249000004

更に、加熱炉スキッド部での抜熱量は、次の(6)式及び(7)式で表される。 Furthermore, the amount of heat removed at the heating furnace skid portion is expressed by the following equations (6) and (7).

Figure 2023125249000005
Figure 2023125249000005

ここで、境界条件における変数について説明すると、qui及びqLiのそれぞれは上面熱負荷[kcal/(m・hr]、下面熱負荷[kcal/(m・hr]、ΦCGUi及びΦCGLiのそれぞれは上面総括熱吸収率[-]、下面総括熱吸収率[-]、θgU及びΦgLのそれぞれは上面炉温[℃]、θsiは1計算ステップ前の上面又は下面のスラブ温度[℃]、Vはシャドーファクター(スキッド影による補正)、Qskdは加熱炉スキッド部での抜熱量(スキッド熱負荷)[kcal/(m・hr]、αskdはスキッド接触熱伝達係数[kcal/(m・hr・K]、θはスラブ温度[℃]、θは加熱炉スキッド部での水温[℃]である。 Here, to explain the variables in the boundary conditions, q ui and q Li are the upper surface heat load [kcal/(m 2・hr], the lower surface heat load [kcal/(m 2・hr]), Φ CGUi and Φ CGLi are the top surface overall heat absorption rate [-], the bottom surface overall heat absorption rate [-], θ gU and Φ gL are the top furnace temperature [℃], and θ si is the slab temperature of the top or bottom surface before one calculation step. [°C], V i is the shadow factor (correction due to skid shadow), Q skd is the amount of heat removed at the skid part of the heating furnace (skid heat load) [kcal/(m 2 · hr]), α skd is the skid contact heat transfer coefficient [kcal/(m 2 ·hr·K], θ is the slab temperature [°C], and θ w is the water temperature in the heating furnace skid part [°C].

熱伝導方程式を示す(1)式~(7)式の情報は、図示しない記録部に記録されており、抽出温度算出部12は、その記録部から熱伝導方程式を示す(1)式~(7)式の情報を取得する。 The information of equations (1) to (7) representing the heat conduction equation is recorded in a recording section (not shown), and the extraction temperature calculation section 12 reads the information of equations (1) to (7) representing the heat conduction equation from the recording section. 7) Obtain expression information.

また、粗出側温度算出部13は、抽出温度算出部12から出力された、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を初期値として、加熱炉2から抽出されたスラブSを粗圧延する際の操業条件(操業条件としては、圧延荷重、圧延速度、冷却などが挙げられる)を用いて、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、熱伝導方程式に基づき算出し、出力する。
熱伝導方程式は、次の(8)式により表される。
In addition, the crude extraction side temperature calculation unit 13 uses the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 as an initial value, which is output from the extraction temperature calculation unit 12. Using the operating conditions (operating conditions include rolling load, rolling speed, cooling, etc.) when rough rolling the slab S extracted from the heating furnace 2, the width center of the slab S after rough rolling is Temperature distribution in the length direction and thickness direction is calculated based on the heat conduction equation and output.
The heat conduction equation is expressed by the following equation (8).

Figure 2023125249000006
Figure 2023125249000006

ここで、ρ:比重[kg/]、C:比熱[kcal/(kg・K)]、λ:熱伝導率[kcal/(m・hr・K)]、θ:温度[℃]、x:長さ方向座標[m]、y:厚さ方向座標[m]である。
この熱伝導方程式に設定される境界条件は、空冷区間、水冷区間、圧延区間のそれぞれに応じて以下のように設定される。
空冷区間の熱負荷qair[kcal/(m・hr]は、次の(9)式により表される。
Here, ρ: specific gravity [kg/ 3 ], C: specific heat [kcal/(kg・K)], λ: thermal conductivity [kcal/(m・hr・K)], θ: temperature [℃], x : Length direction coordinate [m], y: Thickness direction coordinate [m].
The boundary conditions set in this heat conduction equation are set as follows according to each of the air-cooled section, water-cooled section, and rolling section.
The heat load q air [kcal/(m 2 ·hr]) in the air-cooled section is expressed by the following equation (9).

Figure 2023125249000007
Figure 2023125249000007

ここで、ε:放射率[-]、θ:スラブSの表面温度[℃]、θair:大気温度[℃]である。
また、水冷区間の熱負荷qwater[kcal/(m・hr]は、次の(10)式により表される。
Here, ε: emissivity [-], θ S : surface temperature of slab S [°C], θ air : atmospheric temperature [°C].
Further, the heat load q water [kcal/(m 2 ·hr]) in the water cooling section is expressed by the following equation (10).

Figure 2023125249000008
Figure 2023125249000008

ここで、α:水冷熱伝達係数[kcal/(m・hr・K)]、θ:スラブSの表面温度[℃]、θWater:水温[℃]である。
また、圧延区間に関しては、加工発熱、ロール接触抜熱、摩擦発熱の3つの熱バランスを考慮する。
加工発熱ΔQ[kcal/(m・hr]に関しては次の(11)式を用いる。
Here, α w : water cooling heat transfer coefficient [kcal/(m 2 ·hr·K)], θ S : surface temperature of slab S [°C], θ Water : water temperature [°C].
In addition, regarding the rolling section, three types of heat balance are considered: process heat generation, roll contact heat removal, and friction heat generation.
Regarding the processing heat generation ΔQ [kcal/(m 3 ·hr]), the following equation (11) is used.

Figure 2023125249000009
Figure 2023125249000009

ここで、P:圧延圧力[kgf/mm]、hin:入側板厚[mm]、hout:出側板厚[mm]、A:熱の仕事当量427[kgfm/kcal]、τ:圧延接触時間[hr]である。
また、ロール接触抜熱q[kcal/(m・hr)]は次の(12)式で与えられる。
Here, P m : rolling pressure [kgf/mm 2 ], h in : inlet side plate thickness [mm], h out : outlet side plate thickness [mm], A: heat work equivalent 427 [kgfm/kcal], τ: The rolling contact time is [hr].
Further, the roll contact heat removal q r [kcal/(m 2 ·hr)] is given by the following equation (12).

Figure 2023125249000010
Figure 2023125249000010

ここで、λ:スラブSの熱伝導率[kcal/(m・hr・K)]、αrS:スラブSの温度拡散係数[m/hr]、τ:スラブSの接触温度[℃]、T:平衡温度[℃]である。平衡温度Tは、次の(13)式により算出される。 Here, λ S : Thermal conductivity of slab S [kcal/(m・hr・K)], α rS : Temperature diffusion coefficient of slab S [m 2 /hr], τ : Contact temperature of slab S [°C] , T m : Equilibrium temperature [°C]. The equilibrium temperature T m is calculated by the following equation (13).

Figure 2023125249000011
Figure 2023125249000011

ここで、αrS=λ/(ρ・C)、αrR=λ/(ρ・C)であり、λ:スラブSの熱伝導率[kcal/(m・hr・K)]、ρ:スラブSの比重[kg/m]、C:スラブSの比熱[kcal/(kg・K)]、λ:ロールの熱伝導率[kcal/(m・hr・K)]、ρ:ロールの比重[kg/m]、C:ロールの比熱[kcal/(kg・K)]、TR0:ロール温度[℃]、αrR:ロールの温度拡散率[m/hr]である。
また、摩擦発熱qμ[kcal/(m・hr)]は、次の(14)式で与えられる。
Here, α rS = λ S /(ρ S・C S ), α rR = λ R /(ρ R・C R ), and λ S : Thermal conductivity of the slab S [kcal/(m・hr・K)], ρ S : Specific gravity of slab S [kg/m 3 ], C S : Specific heat of slab S [kcal/(kg・K)], λ R : Thermal conductivity of roll [kcal/(m・hr・K)], ρ R : Specific gravity of roll [kg/m 3 ], C R : Specific heat of roll [kcal/(kg・K)], T R0 : Roll temperature [°C], α rR : Temperature diffusion of roll The rate is [m 2 /hr].
Further, frictional heat generation q μ [kcal/(m 2 ·hr)] is given by the following equation (14).

Figure 2023125249000012
Figure 2023125249000012

ここで、μ:摩擦係数[-]、P:圧延圧力[kgf/mm]、VRM:ロールとスラブSの相対速度[m/hr]、A:熱の仕事当量427[kgfm/kcal]である。
圧延圧力Pは次の(15)式により表される。
Here, μ: coefficient of friction [-], P m : rolling pressure [kgf/mm 2 ], V RM : relative speed between roll and slab S [m/hr], A: work equivalent of heat 427 [kgfm/kcal ].
The rolling pressure P m is expressed by the following equation (15).

Figure 2023125249000013
Figure 2023125249000013

ここで、P:圧延荷重[tonf]、R’:偏平ロール径[mm]、hin:入側板厚[mm]、hout:出側板厚[mm]、B:スラブSの板幅[mm]である。
また、相対速度VRMは次の(16)式により表される。
Here, P: rolling load [tonf], R': flat roll diameter [mm], h in : inlet side plate thickness [mm], h out : outlet side plate thickness [mm], B: plate width of slab S [mm] ].
Further, the relative speed VRM is expressed by the following equation (16).

Figure 2023125249000014
Figure 2023125249000014

ここで、V:ロール速度[m/hr]、f:先進率[-]、b:後進率[-]であり、後進率bは次の(17)式により表される。

Figure 2023125249000015
Here, V R is the roll speed [m/hr], f is the forward rate [-], and b is the reverse rate [-], and the reverse rate b is expressed by the following equation (17).
Figure 2023125249000015

粗出側温度算出部13から出力された、粗圧延後のスラブSの裏面の幅中央部における長さ方向の温度分布の一例が図5に示されている。
また、モデル修正部14は、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skid、の偏差e=Tcal_diff-Tact_diffの大きさが閾値ε(1℃)未満となるように、抽出温度算出部12で用いられる、スラブ温度モデルを構成する(1)式~(7)式で表される熱伝導方程式に含まれる加熱炉スキッド部での抜熱量Qskdを修正する。
FIG. 5 shows an example of the temperature distribution in the length direction at the width center of the back surface of the rough-rolled slab S, which is output from the rough-rolling side temperature calculation unit 13.
In addition, the model correction unit 14 calculates the temperature T cal_skid at a longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, which is output from the rough side temperature calculation unit 13, and the rough side temperature T cal_skid . The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to between the heating furnace skids at the width center portion of the back surface Sb of the slab S, which is output from the temperature calculation unit 13, and the rough temperature Temperature Tact_skid at a longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, acquired by the acquisition unit 11, and the back surface Sb of the slab S, acquired by the rough temperature acquisition unit 11. The magnitude of the deviation e=T cal_diff - T act_diff of the temperature difference T act_diff = T act - T act_skid from the temperature T act at the longitudinal position corresponding to between the heating furnace skids at the center width of is the threshold value ε(1 ℃), the amount of heat removed from the heating furnace skid part included in the heat conduction equation expressed by equations (1) to (7) that constitute the slab temperature model used in the extraction temperature calculation unit 12. Correct Q skd .

ここで、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度は複数あり、また、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度も複数ある。ここで、モデル修正部14は、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の複数の温度のうち特定の加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における前述の特定の加熱炉スキッド部とこの特定の加熱炉スキッド部に隣接する加熱炉スキッド部との間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidを算出する。 Here, there are a plurality of temperatures at the longitudinal position corresponding to the heating furnace skid part in the width center portion of the back surface Sb of the slab S, which are output from the roughing side temperature calculating section 13, and the roughing side temperature calculating section There are also a plurality of temperatures outputted from 13 at longitudinal positions corresponding to between the heating furnace skids at the width center portion of the back surface Sb of the slab S. Here, the model correction unit 14 selects a specific temperature among the plurality of temperatures at the longitudinal position corresponding to the heating furnace skid part in the width center of the back surface Sb of the slab S, which is output from the raw side temperature calculation unit 13. The temperature T cal_skid at the longitudinal position corresponding to the heating furnace skid portion, the above-mentioned specific heating furnace skid portion at the width center portion of the back surface Sb of the slab S, and this identification, which is output from the raw side temperature calculation unit 13. The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to the heating furnace skid part adjacent to the heating furnace skid part is calculated.

また、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度は複数あり、また、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度も複数ある。ここで、モデル修正部14は、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の複数の温度のうち前述の特定の加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における前述の特定の加熱炉スキッド部とこの特定の加熱炉スキッド部に隣接する加熱炉スキッド部との間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skidを算出する。 In addition, there are multiple temperatures at the longitudinal position corresponding to the heating furnace skid part in the width center of the back surface Sb of the slab S, which are acquired by the rough temperature acquisition unit 11. , there are also a plurality of temperatures at longitudinal positions corresponding to between the heating furnace skids at the width center portion of the back surface Sb of the slab S. Here, the model correction unit 14 selects the above-mentioned specific temperature among the plurality of temperatures at the longitudinal position corresponding to the heating furnace skid part in the width center portion of the back surface Sb of the slab S, acquired by the rough temperature acquisition unit 11. The temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion, the above-mentioned specific heating furnace skid portion at the width center of the back surface Sb of the slab S, and this specific heating acquired by the rough temperature acquisition unit 11. A temperature difference T act_diff = T act - T act_skid between the temperature T act at a corresponding longitudinal position between the furnace skid part and the adjacent heating furnace skid part is calculated.

そして、モデル修正部14は、算出された温度差Tcal_diff=Tcal-Tcal_skidと温度差Tact_diff=Tact-Tact_skidとの偏差e=Tcal_diff-Tact_diffの大きさが閾値ε(1℃)未満となるように、抽出温度算出部12で用いられる、スラブ温度モデルを構成する(1)式~(7)式で表される熱伝導方程式に含まれる加熱炉スキッド部での抜熱量Qskdを修正する。
加熱炉スキッド部での抜熱量は、前述したように、(7)式で表されるが、(7)式に補正係数ηを導入した次の(18)式における適切なηを求めることが加熱炉スキッド部での抜熱量Qskdの適切な修正となる。
Then, the model correction unit 14 determines that the magnitude of the deviation e=T cal_diff -T act_diff between the calculated temperature difference T cal_diff = T cal -T cal_skid and the temperature difference T act_diff =T act -T act_skid is a threshold value ε(1 ℃), the amount of heat removed from the heating furnace skid part included in the heat conduction equation expressed by equations (1) to (7) that constitute the slab temperature model used in the extraction temperature calculation unit 12. Correct Q skd .
As mentioned above, the amount of heat removed in the heating furnace skid part is expressed by equation (7), but it is possible to find an appropriate η in the following equation (18) by introducing the correction coefficient η into equation (7). This is an appropriate correction for the amount of heat removed Qskd at the heating furnace skid section.

Figure 2023125249000016
Figure 2023125249000016

この(18)式における適切なηを求める方法については、図7に示すステップS4(モデル修正ステップ)における処理の流れの詳細を示すフローチャートに基づき後述する。
そして、修正装置10の出力部15は、モデル修正部14によって加熱炉スキッド部での抜熱量Qskdを修正したスラブ温度モデルを構成する熱伝導方程式((1)式~(7)式、(18)式)の情報を、加熱炉の炉温制御装置20のスラブ温度予測部21に対し出力する。
A method for determining an appropriate η in equation (18) will be described later based on a flowchart showing details of the process flow in step S4 (model modification step) shown in FIG.
Then, the output unit 15 of the correction device 10 outputs heat conduction equations (Equations ( 1 ) to (7), (Equations (1) to (7), The information in equation 18) is output to the slab temperature prediction unit 21 of the furnace temperature control device 20 of the heating furnace.

また、加熱炉の炉温制御装置20は、図3に示すように、スラブ温度予測部21と、炉温制御部22とを備えている。炉温制御装置20は、演算処理機能を有するコンピュータシステムであり、ハードウェアに予め記憶された各種専用のコンピュータプログラムの命令により、スラブ温度予測部21(後に述べるステップS11)、及び炉温制御部22(ステップS12)の各機能をソフトウェア上で実行する。 Further, the furnace temperature control device 20 for the heating furnace includes a slab temperature prediction section 21 and a furnace temperature control section 22, as shown in FIG. The furnace temperature control device 20 is a computer system having an arithmetic processing function, and operates a slab temperature prediction section 21 (step S11 described later) and a furnace temperature control section according to instructions from various dedicated computer programs stored in advance in hardware. 22 (step S12) are executed on software.

スラブ温度予測部21は、スラブ温度モデルの修正装置10によって修正されたスラブ温度モデルの情報を修正装置10の出力部15から取得し、その情報を用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測する。 The slab temperature prediction unit 21 acquires information on the slab temperature model corrected by the slab temperature model correction device 10 from the output unit 15 of the correction device 10, and uses that information to estimate the value of the slab S at the time of extraction in the heating furnace 2. Predict the temperature distribution in the length direction and thickness direction at the width center.

具体的に述べると、スラブ温度予測部21は、装入温度測定計7により測定された加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度を初期値として、スラブSの加熱炉在炉中の前述した操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、修正装置10の出力部15から取得した修正されたスラブ温度モデルを構成する熱伝導方程式((1)式~(7)式、(18)式)に基づき算出し、この算出された値を予測値とする。 Specifically, the slab temperature prediction unit 21 uses the charging temperatures of the front surface Sa and back surface Sb of the slab S at the time of charging the heating furnace 2 measured by the charging temperature measuring meter 7 as initial values, and calculates the charging temperature of the slab S. Using the above-mentioned operating conditions while the heating furnace is in use, the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 is obtained from the output unit 15 of the correction device 10. The calculated value is calculated based on the heat conduction equations (Equations (1) to (7), and Equation (18)) constituting the revised slab temperature model, and this calculated value is used as the predicted value.

また、炉温制御部22は、スラブ温度予測部21で予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、加熱炉2内における炉温を制御する。 The furnace temperature control unit 22 also controls the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S at the time of extraction of the heating furnace 2, which is predicted by the slab temperature prediction unit 21. The furnace temperature in the heating furnace 2 is controlled so that the temperature exceeds the minimum extraction temperature.

具体的に述べると、炉温制御部22は、スラブ温度予測部21で予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度(例えば、1200℃)を上回る目標温度となるように、各燃料流量制御部23を制御して各燃焼バーナー3への燃料流量を制御し、加熱炉2内の炉温を制御する。 Specifically, the furnace temperature control unit 22 controls the thickness at the longitudinal position corresponding to the heating furnace skid portion in the width center of the slab S at the time of extraction of the heating furnace 2, which is predicted by the slab temperature prediction unit 21. Each fuel flow control unit 23 is controlled to control the fuel flow rate to each combustion burner 3 so that the temperature in the horizontal direction reaches a target temperature exceeding the minimum extraction temperature (for example, 1200°C), and the fuel flow rate in the heating furnace 2 is controlled. Control the furnace temperature.

これにより、加熱炉2から抽出された時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回り、昇温不足となりがちなスキッドビームの接触箇所を規定温度まで上昇させて抽出することができることになる。 As a result, when the slab S is extracted from the heating furnace 2, the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S exceeds the minimum extraction temperature, which tends to result in insufficient temperature rise. This means that the contact point of the skid beam can be raised to a specified temperature for extraction.

このように本実施形態に係るスラブ温度モデルの修正装置10によれば、粗温度取得部11において、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を取得し、抽出温度算出部12において、加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式((1)式~(7)式)に基づき算出し、出力する。そして、粗出側温度算出部13において、抽出温度算出部12から出力された、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を初期値として、加熱炉2から抽出されたスラブSを粗圧延する際の操業条件を用いて、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、熱伝導方程式((8)式~(17)式)に基づき算出し、出力する。そして、モデル修正部14において、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skid、の偏差e=Tcal_diff-Tact_diffの大きさが閾値ε(1℃)未満となるように、抽出温度算出部12で用いられる、スラブ温度モデルを構成する(1)式~(7)式で表される熱伝導方程式に含まれる加熱炉スキッド部での抜熱量Qskdを修正する((18)式)。 As described above, according to the slab temperature model correction device 10 according to the present embodiment, the rough temperature acquisition unit 11 calculates the rough temperature measured by the rough temperature measuring meters 8 installed on the outlet sides of the roughing mills 4 1 to 4 n . The temperature distribution in the longitudinal direction at the width center of the back surface Sb of the slab S after rolling is acquired, and the extraction temperature calculation unit 12 calculates the temperature distribution of the slab S by using the charging temperature of the slab S charged into the heating furnace 2 as an initial value. Using the operating conditions of S during heating furnace 2, the temperature distribution in the length direction and thickness direction at the width center of slab S at the time of extraction in heating furnace 2 is calculated using the heat conduction equation ( Calculate based on equations (1) to (7) and output. Then, in the crude extraction side temperature calculation unit 13, the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2, which is output from the extraction temperature calculation unit 12, is set as an initial value. Using the operating conditions when rough rolling the slab S extracted from the heating furnace 2, the temperature distribution in the length direction and thickness direction at the width center of the slab S after rough rolling is calculated using the heat conduction equation ((8 ) to (17)) and output. Then, in the model correction unit 14, the temperature T cal_skid at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S, which is output from the rough side temperature calculation unit 13, and the rough side temperature T cal_skid The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to between the heating furnace skids at the width center portion of the back surface Sb of the slab S, which is output from the temperature calculation unit 13, and the rough temperature Temperature Tact_skid at a longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, acquired by the acquisition unit 11, and the back surface Sb of the slab S, acquired by the rough temperature acquisition unit 11. The magnitude of the deviation e=T cal_diff - T act_diff of the temperature difference T act_diff = T act - T act_skid from the temperature T act at the longitudinal position corresponding to between the heating furnace skids at the center width of is the threshold value ε(1 ℃), the amount of heat removed from the heating furnace skid part included in the heat conduction equation expressed by equations (1) to (7) that constitute the slab temperature model used in the extraction temperature calculation unit 12. Correct Q skd (formula (18)).

これにより、スラブ温度モデルを用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測するに際し、当該スラブ温度モデルを修正することで、加熱炉2の抽出時のスラブの幅中央部における長さ方向及び厚さ方向の温度分布を精度よく予測することができることになる。 As a result, when predicting the temperature distribution in the length direction and thickness direction in the width center of the slab S at the time of extraction in the heating furnace 2 using the slab temperature model, by correcting the slab temperature model, the heating furnace It is possible to accurately predict the temperature distribution in the length direction and thickness direction at the width center portion of the slab during extraction in step 2.

また、このスラブ温度モデルの修正に介し、モデル修正部14において、粗出側温度算出部13から出力された情報と、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗温度取得部11で取得された情報とに基づいてスラブ温度モデルを修正しているので、加熱炉2の出側直後に設置される放射温度計によって加熱炉2の抽出直後のスラブSの温度を測定する必要はない。加熱炉2の抽出直後だとスラブSの表面のスケールの影響で温度測定自体が困難であるが、本実施形態に係る修正装置10によれば、このおそれはない。 In addition, through the correction of this slab temperature model, the model correction unit 14 uses the information output from the rough delivery side temperature calculation unit 13 and the rough temperature measurement meters installed on the exit sides of the rough rolling mills 4 1 to 4 n . Since the slab temperature model is corrected based on the information acquired by the rough temperature acquisition unit 11 measured in step 8, the radiation thermometer installed immediately after the exit side of the heating furnace 2 measures the temperature immediately after the extraction of the heating furnace 2. There is no need to measure the temperature of the slab S. Immediately after extraction from the heating furnace 2, temperature measurement itself is difficult due to the influence of scale on the surface of the slab S, but according to the repair device 10 according to the present embodiment, there is no such fear.

そして、本実施形態に係る加熱炉の炉温制御装置20によれば、スラブ温度予測部21において、スラブ温度モデルの修正装置10によって修正されたスラブ温度モデルを用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測し、炉温制御部22において、スラブ温度予測部21で予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、加熱炉2内における炉温を制御する。 According to the furnace temperature control device 20 of the heating furnace according to the present embodiment, the slab temperature prediction unit 21 uses the slab temperature model corrected by the slab temperature model correction device 10 to The temperature distribution in the length direction and the thickness direction at the width center of the slab S is predicted, and the width center of the slab S at the time of extraction of the heating furnace 2 predicted by the slab temperature prediction unit 21 is calculated in the furnace temperature control unit 22. The furnace temperature in the heating furnace 2 is controlled so that the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid exceeds the minimum extraction temperature.

これにより、加熱炉2から抽出された時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回り、昇温不足となりがちなスキッドビームの接触箇所を規定温度まで上昇させて抽出することができ、加熱炉スキッド部の昇温不足による熱延コイルの材質不良が防止できるとともに、過加熱操業による加熱原単位を削減することができる。 As a result, when the slab S is extracted from the heating furnace 2, the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S exceeds the minimum extraction temperature, which tends to result in insufficient temperature rise. It is possible to raise the contact point of the skid beam to a specified temperature and extract it, preventing material defects in the hot-rolled coil due to insufficient temperature rise of the skid part of the heating furnace, and reducing the heating unit consumption due to overheating operation. can.

次に、本発明の一実施形態に係るスラブ温度モデルの修正方法を、図6に示されたスラブ温度モデルの修正装置における処理の流れを説明するためのフローチャート、及び図7に示された図6に示すフローチャートにおけるステップS4(モデル修正ステップ)における処理の流れを示すフローチャートを参照して説明する。 Next, a method for correcting a slab temperature model according to an embodiment of the present invention will be explained by a flowchart for explaining the process flow in the slab temperature model correction apparatus shown in FIG. 6, and a diagram shown in FIG. This will be explained with reference to a flowchart showing the flow of processing in step S4 (model modification step) in the flowchart shown in FIG.

先ず、ステップS1において、修正装置10の粗温度取得部11は、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を取得する(粗温度取得ステップ)。 First, in step S1, the rough temperature acquisition unit 11 of the correction device 10 measures the back surface Sb of the slab S after rough rolling, which is measured by the rough temperature measuring meter 8 installed on the outlet side of the rough rolling mills 4 1 to 4 n . Acquire the temperature distribution in the length direction at the center width (rough temperature acquisition step).

次いで、ステップS2において、修正装置10の抽出温度算出部12は、加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件(操業条件としては、在炉時間、炉温、スキッド接触時間などが挙げられる)を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する(抽出温度算出ステップ)。 Next, in step S2, the extraction temperature calculation unit 12 of the correction device 10 sets the charging temperature of the slab S charged into the heating furnace 2 as an initial value, and calculates the operating conditions (operating conditions) while the slab S is in the heating furnace. The temperature distribution in the length direction and the thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 is calculated using Calculate and output based on the heat conduction equation that constitutes the model (extraction temperature calculation step).

加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度が装入温度測定計7により測定され、この測定された加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度が抽出温度算出部12に入力される。そして、抽出温度算出部12は、このスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の前述した操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する前述の(1)式~(7)式で表される熱伝導方程式に基づき算出し、出力する。 The charging temperature of the front surface Sa and back surface Sb of the slab S when charging into the heating furnace 2 is measured by the charging temperature measuring meter 7, and the charging temperature of the front surface Sa and back surface of the slab S when charging into the heating furnace 2 is measured. The Sb charging temperature is input to the extraction temperature calculation section 12. Then, the extraction temperature calculation unit 12 uses the charging temperature of the slab S as an initial value and the width center of the slab S at the time of extraction in the heating furnace 2, using the above-mentioned operating conditions while the slab S is in the heating furnace. The temperature distribution in the length direction and thickness direction in the section is calculated based on the heat conduction equations expressed by the above-mentioned equations (1) to (7) that constitute the slab temperature model, and is output.

(1)式の熱伝導方程式を用いての温度分布の算出では、時間方向と空間方向に差分化してスラブSの表面Sa及び裏面Sbでの境界条件を設定することで、スラブ温度を時間刻み毎に抽出時まで算出して、出力する。これにより、加熱炉スキッド部により生じるスラブSの長さ方向の温度変動が得られる。
この熱伝導方程式に設定される境界条件は、加熱炉2からのスラブSへの輻射入熱と加熱炉スキッド部での抜熱量とからなる。
表面(上面)Saからの輻射入熱は、前述した通りの(2)式及び(3)式により表される。
また、裏面(下面)Sbからの輻射入熱は、前述した通りの(4)式及び(5)式で表される。
更に、加熱炉スキッド部での抜熱量は、前述した通りの(6)式及び(7)式で表される。
In calculating the temperature distribution using the heat conduction equation of equation (1), the slab temperature is calculated in time increments by setting the boundary conditions on the front surface Sa and the back surface Sb of the slab S by differentiating them in the time and space directions. Calculates up to the time of extraction each time and outputs it. Thereby, temperature fluctuations in the length direction of the slab S caused by the heating furnace skid portion can be obtained.
The boundary conditions set in this heat conduction equation include the radiant heat input from the heating furnace 2 to the slab S and the amount of heat removed at the heating furnace skid portion.
The radiant heat input from the surface (upper surface) Sa is expressed by equations (2) and (3) as described above.
In addition, the radiant heat input from the back surface (lower surface) Sb is expressed by equations (4) and (5) as described above.
Furthermore, the amount of heat removed in the heating furnace skid portion is expressed by the above-mentioned equations (6) and (7).

次いで、ステップS3において、修正装置10の粗出側温度算出部13は、ステップS2で出力された、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を初期値として、加熱炉2から抽出されたスラブSを粗圧延する際の操業条件(操業条件としては、圧延荷重、圧延速度、冷却などが挙げられる)を用いて、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、熱伝導方程式に基づき算出し、出力する(粗出側温度算出ステップ)。
熱伝導方程式は、前述した通りの(8)式により表される。
また(8)式の熱伝導方程式に設定される境界条件は、前述した通り、空冷区間、水冷区間、圧延区間のそれぞれに応じて設定される。
空冷区間の熱負荷qair[kcal/(m・hr]は、前述の(9)式により表される。
Next, in step S3, the raw temperature calculation unit 13 of the correction device 10 calculates the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2, which is output in step S2. As an initial value, using the operating conditions (operating conditions include rolling load, rolling speed, cooling, etc.) when rough rolling the slab S extracted from the heating furnace 2, the slab S after rough rolling is The temperature distribution in the length direction and the thickness direction at the center width of is calculated based on the heat conduction equation and output (raw side temperature calculation step).
The heat conduction equation is expressed by equation (8) as described above.
Further, the boundary conditions set in the heat conduction equation of equation (8) are set according to each of the air-cooled section, water-cooled section, and rolling section, as described above.
The heat load q air [kcal/(m 2 ·hr]) in the air-cooled section is expressed by the above-mentioned equation (9).

また、水冷区間の熱負荷qwater[kcal/(m・hr]は、前述の(10)式により表される。
また、圧延区間に関しては、加工発熱、ロール接触抜熱、摩擦発熱の3つの熱バランスを考慮する。
加工発熱ΔQ[kcal/(m・hr]に関しては前述の(11)式を用いる。
また、ロール接触抜熱q[kcal/(m・hr)]は前述の(12)式で与えられる。
衡温度[℃]である。平衡温度Tは、前述の(13)式により算出される。
また、摩擦発熱qμ[kcal/(m・hr)]は、前述の(14)式で与えられる。圧延圧力Pは前述の(15)式により表される。また、相対速度VRMは前述の(16)式により表される。また、後進率bは前述の(17)式により表される。
Further, the heat load q water [kcal/(m 2 ·hr]) in the water-cooled section is expressed by the above-mentioned equation (10).
In addition, regarding the rolling section, three types of heat balance are considered: process heat generation, roll contact heat removal, and friction heat generation.
Regarding the processing heat generation ΔQ [kcal/(m 3 ·hr]), the above-mentioned equation (11) is used.
Further, the roll contact heat removal q r [kcal/(m 2 ·hr)] is given by the above-mentioned equation (12).
Equilibrium temperature [°C]. The equilibrium temperature T m is calculated by the above-mentioned equation (13).
Furthermore, the frictional heat generation q μ [kcal/(m 2 ·hr)] is given by the above-mentioned equation (14). The rolling pressure P m is expressed by the above-mentioned equation (15). Further, the relative speed VRM is expressed by the above-mentioned equation (16). Further, the backward movement rate b is expressed by the above-mentioned equation (17).

次いで、ステップS4において、モデル修正部14は、ステップS3で出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、ステップS3で出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidと、ステップS1で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、ステップS1で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skid、の偏差e=Tcal_diff-Tact_diffの大きさ(絶対値の大きさ)が閾値ε(1℃)未満となるように、ステップS2で用いられる、スラブ温度モデルを構成する(1)式~(7)式で表される熱伝導方程式に含まれる加熱炉スキッド部での抜熱量Qskdを修正する(モデル修正ステップ)。 Next, in step S4, the model correction unit 14 calculates the temperature T cal_skid at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S, which was output in step S3, and the temperature T cal_skid output in step S3. The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to between the heating furnace skids at the width center of the back surface Sb of the slab S, and the temperature difference between the slab S obtained in step S1 The temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S corresponds to the temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, which is obtained in step S1. The magnitude (absolute value) of the temperature difference T act_diff = T act - T act_skid from the temperature T act at the longitudinal position, e = T cal_diff - T act_diff , is less than the threshold ε (1°C). As shown in FIG. step).

加熱炉スキッド部での抜熱量は、前述したように、(7)式で表されるが、(7)式に補正係数ηを導入した前述の(18)式における適切なηを求めることが加熱炉スキッド部での抜熱量Qskdの適切な修正となる。 As mentioned above, the amount of heat removed in the heating furnace skid section is expressed by equation (7), but it is possible to find an appropriate η in equation (18), which introduces the correction coefficient η into equation (7). This is an appropriate correction for the amount of heat removed Qskd at the heating furnace skid section.

ステップS4での処理を、図7を参照して詳細に説明する。
先ず、ステップS41において、モデル修正部14は、ステップS3で出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、ステップS3で出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidと、ステップS1で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、ステップS1で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skid、の偏差e=Tcal_diff-Tact_diffを算出する。
The process in step S4 will be explained in detail with reference to FIG.
First, in step S41, the model correction unit 14 calculates the temperature T cal_skid at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S, which was output in step S3, and The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to between the heating furnace skids at the width center of the back surface Sb of the slab S, and the temperature difference between the slab S obtained in step S1 The temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S corresponds to the temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, which is obtained in step S1. The deviation e =T cal_diff - T act_diff of the temperature difference T act_diff = T act - T act_skid from the temperature T act at the longitudinal position is calculated.

ここで、モデル修正部14は、ステップS3で出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の複数の温度のうち特定の加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、ステップS3で出力された、スラブSの裏面Sbの幅中央部における前述の特定の加熱炉スキッド部とこの特定の加熱炉スキッド部に隣接する加熱炉スキッド部との間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidを算出する。 Here, the model correction unit 14 selects a specific heating furnace skid portion from among the plurality of temperatures at a longitudinal position corresponding to the heating furnace skid portion in the width center portion of the back surface Sb of the slab S, output in step S3. The temperature T cal_skid at the corresponding longitudinal position, the above-mentioned specific heating furnace skid portion at the width center of the back surface Sb of the slab S, and the heating furnace adjacent to this specific heating furnace skid portion, output in step S3. The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to the skid part is calculated.

また、モデル修正部14は、ステップS1で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の複数の温度のうち前述の特定の加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、ステップS1で取得された、スラブSの裏面Sbの幅中央部における前述の特定の加熱炉スキッド部とこの特定の加熱炉スキッド部に隣接する加熱炉スキッド部との間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skidを算出する。
そして、モデル修正部14は、算出された温度差Tcal_diff=Tcal-Tcal_skidと温度差Tact_diff=Tact-Tact_skidとの偏差e=Tcal_diff-Tact_diff(絶対値の大きさ)を算出する。
In addition, the model correction unit 14 selects the above-mentioned specific heating furnace skid portion from among the plurality of temperatures at the longitudinal position corresponding to the heating furnace skid portion in the width center portion of the back surface Sb of the slab S, acquired in step S1. The temperature T act_skid at the longitudinal position corresponding to the above-mentioned specific heating furnace skid portion at the width center of the back surface Sb of the slab S and the heating adjacent to this specific heating furnace skid portion obtained in step S1. The temperature difference T act_diff = T act - T act_skid between the temperature T act at the corresponding longitudinal position between the furnace skid portion and the temperature T act is calculated.
Then, the model correction unit 14 calculates the deviation e=T cal_diff - T act_diff (magnitude of absolute value) between the calculated temperature difference T cal_diff = T cal - T cal_skid and the temperature difference T act_diff = T act - T act_skid . calculate.

次いで、ステップS42において、モデル修正部14は、ステップS41で算出した偏差e(絶対値の大きさ)が閾値ε(=1℃)未満か否かを判定する。
ステップS42における判定結果がNO(偏差eが閾値ε未満でない)の場合、ステップS43に移行し、ステップS42における判定結果がYES(偏差eが閾値ε未満である)の場合、ステップS49に移行する。
ステップS43では、モデル修正部14は、前述の(18)式における補正係数ηを微小変化後の補正係数η+δに設定する。
Next, in step S42, the model correction unit 14 determines whether the deviation e (magnitude of absolute value) calculated in step S41 is less than a threshold value ε (=1° C.).
If the determination result in step S42 is NO (the deviation e is not less than the threshold ε), the process moves to step S43, and if the determination result in step S42 is YES (the deviation e is less than the threshold ε), the process moves to step S49. .
In step S43, the model correction unit 14 sets the correction coefficient η in the above-mentioned equation (18) to the correction coefficient η+δ after the slight change.

次いで、ステップS44において、モデル修正部14は、補正係数ηを微小変化後の補正係数η+δとした(18)式の加熱炉スキッド部での抜熱量Qskdを用いて、ステップS2における、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出し直し、そして、これに基づき、ステップS3での粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布の算出をし直し、ステップS41における、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skid’と、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcal’との温度差Tcal_diff’=Tcal’-Tcal_skid’を算出する。 Next, in step S44, the model correction unit 14 uses the amount of heat removed from the heating furnace skid part Q skd in equation (18), where the correction coefficient η is the correction coefficient η+δ after the slight change, to improve the heating furnace in step S2. The temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in step S3 is recalculated, and based on this, the length at the width center of the slab S after rough rolling in step S3 is calculated. The temperature distribution in the direction and the thickness direction is recalculated, and the temperature T cal_skid ' at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S and the back surface of the slab S in step S41 is calculated. The temperature difference T cal_diff ′=T cal ′ −T cal_skid ′ between the temperature T cal ′ at the longitudinal position corresponding to between the heating furnace skids at the center width of Sb is calculated.

次いで、ステップS45において、モデル修正部14は、ステップS41で算出した温度差Tcal_diffと、ステップS44で算出した温度差Tcal_diff’と、ステップ43で設定した補正係数の微小変化量δとから、影響係数∂e/∂η=(Tcal_diff’-Tcal_diff)/δを算出する。 Next, in step S45, the model correction unit 14 calculates, from the temperature difference T cal_diff calculated in step S41, the temperature difference T cal_diff ' calculated in step S44, and the minute change amount δ of the correction coefficient set in step S43, An influence coefficient ∂e/∂η=(T cal_diff '−T cal_diff )/δ is calculated.

次いで、ステップS46において、モデル修正部14は、前述の(18)式における補正係数ηを更新後補正係数η-e/(∂e/∂η)に更新する。 Next, in step S46, the model correction unit 14 updates the correction coefficient η in the above-mentioned equation (18) to the updated correction coefficient η−e/(∂e/∂η).

次いで、ステップS47において、モデル修正部14は、補正係数ηを更新後補正係数η-e/(∂e/∂η)とした(18)式の加熱炉スキッド部での抜熱量Qskdを用いて、ステップS2における、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出し直し、ステップS3での粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布の算出をし直し、ステップS41における、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skid’’と、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcal’’との温度差Tcal_diff’’=Tcal’’-Tcal_skid’’を算出する。 Next, in step S47, the model correction unit 14 uses the amount of heat removed Q skd in the heating furnace skid part of equation (18) where the correction coefficient η is the updated correction coefficient η−e/(∂e/∂η). Then, the temperature distribution in the length direction and the thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 in step S2 is recalculated, and the temperature distribution in the width center of the slab S after rough rolling in step S3 is calculated again. The temperature distribution in the length direction and the thickness direction is recalculated, and the temperature T cal_skid '' at the longitudinal position corresponding to the heating furnace skid part in the width center of the back surface Sb of the slab S in step S41 and the slab The temperature difference T cal_diff '' between the temperature T cal '' at the longitudinal position corresponding to between the heating furnace skids at the center width of the back surface Sb of S is calculated.

次いで、ステップS48において、モデル修正部14は、ステップS41で算出した偏差eを偏差e=Tcal_diff’’-Tact_diffに更新する。
そして、ステップ42に戻り、ステップS42において、モデル修正部14は、ステップS48で更新した偏差eが閾値ε(=1℃)未満か否かを判定する。
そして、ステップS42における判定結果がYESとなるまでステップS42からステップS48までを繰り返す。
ステップS42における判定結果がYESとなった場合、ステップS49に移行し、偏差eが閾値ε(=1℃)未満となったときの補正係数ηを補正係数とする。
Next, in step S48, the model correction unit 14 updates the deviation e calculated in step S41 to deviation e=T cal_diff ''-T act_diff .
Then, the process returns to step S42, and in step S42, the model correction unit 14 determines whether the deviation e updated in step S48 is less than the threshold value ε (=1° C.).
Then, steps S42 to S48 are repeated until the determination result in step S42 becomes YES.
If the determination result in step S42 is YES, the process moves to step S49, and the correction coefficient η when the deviation e becomes less than the threshold value ε (=1° C.) is set as the correction coefficient.

これにより、ステップS4での処理が終了し、前述の(18)式における適切な補正係数ηが求められ、加熱炉スキッド部での抜熱量Qskdの適切な修正となる。
そして、ステップS5において、修正装置10の出力部15は、モデル修正部14によって加熱炉スキッド部での抜熱量Qskdを修正したスラブ温度モデルを構成する熱伝導方程式((1)式~(7)式、(18)式)の情報を、加熱炉の炉温制御装置20のスラブ温度予測部21に対し出力する(出力ステップ)。
Thereby, the process in step S4 is completed, and an appropriate correction coefficient η in the above-mentioned equation (18) is determined, and the amount of heat removed Q skd in the heating furnace skid portion is appropriately corrected.
Then, in step S5, the output unit 15 of the correction device 10 outputs heat conduction equations (Equations (1) to (7) ) and (18) are output to the slab temperature prediction unit 21 of the furnace temperature control device 20 of the heating furnace (output step).

これにより、スラブ温度モデルの修正装置10における処理が終了する。 This completes the processing in the slab temperature model correction device 10.

このように本実施形態に係るスラブ温度モデルの修正方法によれば、粗温度取得ステップ(ステップS1)において、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を取得し、抽出温度算出ステップ(ステップS2)において、加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、スラブ温度モデルを構成する熱伝導方程式((1)式~(7)式)に基づき算出し、出力する。そして、粗出側温度算出ステップ(ステップS3)において、抽出温度算出ステップで出力された、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を初期値として、加熱炉2から抽出されたスラブSを粗圧延する際の操業条件を用いて、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、熱伝導方程式((8)式~(17)式)に基づき算出し、出力する。そして、モデル修正ステップ(ステップS4)において、粗出側温度算出ステップで出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、粗出側温度算出部13から出力された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、粗温度取得部11で取得された、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skid、の偏差e=Tcal_diff-Tact_diffの大きさが閾値ε(1℃)未満となるように、抽出温度算出ステップで用いられる、スラブ温度モデルを構成する(1)式~(7)式で表される熱伝導方程式に含まれる加熱炉スキッド部での抜熱量Qskdを修正する((18)式)。 As described above, according to the slab temperature model correction method according to the present embodiment, in the rough temperature acquisition step (step S1), the rough temperature measurement meter 8 installed on the outlet side of the roughing mills 4 1 to 4 n The temperature distribution in the longitudinal direction at the width center of the back surface Sb of the rough-rolled slab S is obtained, and in the extraction temperature calculation step (step S2), the charging temperature of the slab S to be charged into the heating furnace 2 is calculated. As an initial value, a slab temperature model is constructed using the operating conditions of the slab S while it is in the heating furnace, and the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2. Calculate based on the heat conduction equation (Equations (1) to (7)) and output. Then, in the crude extraction side temperature calculation step (step S3), the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2, which was output in the extraction temperature calculation step, is set to an initial value. Using the operating conditions when rough rolling the slab S extracted from the heating furnace 2, the temperature distribution in the length direction and thickness direction at the width center of the slab S after rough rolling is calculated using the heat conduction equation ( Calculate based on equations (8) to (17) and output. Then, in the model correction step (step S4), the temperature T cal_skid at the longitudinal position corresponding to the heating furnace skid part in the width center of the back surface Sb of the slab S, which was output in the rough output side temperature calculation step, and the rough The temperature difference T cal_diff = T cal - T cal_skid between the temperature T cal at the longitudinal position corresponding to between the heating furnace skids at the width center portion of the back surface Sb of the slab S, which is output from the outlet temperature calculation unit 13; The temperature T act_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, acquired by the rough temperature acquisition unit 11, and the temperature T act_skid of the slab S acquired by the rough temperature acquisition unit 11. The size of the deviation e=T cal_diff - T act_diff of the temperature difference T act_diff = T act - T act_skid between the temperature T act at the longitudinal position corresponding to between the heating furnace skids at the center width of the back surface Sb is the threshold value ε. (1°C) or less at the heating furnace skid, which is included in the heat conduction equation expressed by equations (1) to (7) that constitute the slab temperature model used in the extraction temperature calculation step. Correct the amount of heat Q skd (formula (18)).

これにより、スラブ温度モデルを用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測するに際し、当該スラブ温度モデルを修正することで、加熱炉2の抽出時のスラブの幅中央部における長さ方向及び厚さ方向の温度分布を精度よく予測することができることになる。 As a result, when predicting the temperature distribution in the length direction and thickness direction in the width center of the slab S at the time of extraction in the heating furnace 2 using the slab temperature model, by correcting the slab temperature model, the heating furnace It is possible to accurately predict the temperature distribution in the length direction and thickness direction at the width center portion of the slab during extraction in step 2.

また、このスラブ温度モデルの修正に介し、モデル修正ステップにおいて、粗出側温度算出ステップから出力された情報と、粗圧延機4~4の出側に設置された粗温度測定計8によって測定した粗温度取得ステップで取得された情報とに基づいてスラブ温度モデルを修正しているので、加熱炉2の出側直後に設置される放射温度計によって加熱炉2の抽出直後のスラブSの温度を測定する必要はない。加熱炉2の抽出直後だとスラブSの表面のスケールの影響で温度測定自体が困難であるが、本実施形態に係る修正方法によれば、このおそれはない。 In addition, through the correction of this slab temperature model, in the model correction step, the information output from the rough delivery side temperature calculation step and the rough temperature measurement meter 8 installed on the exit side of the roughing mills 4 1 to 4 n are used. Since the slab temperature model is corrected based on the information acquired in the measured crude temperature acquisition step, the radiation thermometer installed immediately after the exit side of the heating furnace 2 can be used to determine the temperature of the slab S immediately after extraction from the heating furnace 2. There is no need to measure temperature. Immediately after extraction from the heating furnace 2, temperature measurement itself is difficult due to the influence of scale on the surface of the slab S, but according to the correction method according to the present embodiment, there is no such fear.

次に、本発明の一実施形態に係る加熱炉の炉温制御方法を、図8に示された加熱炉の炉温制御装置における処理の流れを説明するためのフローチャートを参照して説明する。 Next, a method for controlling the furnace temperature of a heating furnace according to an embodiment of the present invention will be described with reference to a flowchart for explaining the flow of processing in the furnace temperature control device for a heating furnace shown in FIG.

先ず、ステップS11において、炉温制御装置20のスラブ温度予測部21は、スラブ温度モデルの修正方法によって修正されたスラブ温度モデルの情報を修正装置10の出力部15から取得し、その情報を用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測する(スラブ温度予測ステップ)。 First, in step S11, the slab temperature prediction unit 21 of the furnace temperature control device 20 acquires information on the slab temperature model corrected by the slab temperature model correction method from the output unit 15 of the correction device 10, and uses the information. The temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 is predicted (slab temperature prediction step).

具体的に述べると、スラブ温度予測部21は、装入温度測定計7により測定された加熱炉2の装入時のスラブSの表面Sa及び裏面Sbの装入温度を初期値として、スラブSの加熱炉在炉中の前述した操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、修正装置10の出力部15から取得した修正されたスラブ温度モデルを構成する熱伝導方程式((1)式~(7)式、(18)式)に基づき算出し、この算出された値を予測値とする。 Specifically, the slab temperature prediction unit 21 uses the charging temperatures of the front surface Sa and back surface Sb of the slab S at the time of charging the heating furnace 2 measured by the charging temperature measuring meter 7 as initial values, and calculates the charging temperature of the slab S. Using the above-mentioned operating conditions while the heating furnace is in use, the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 is obtained from the output unit 15 of the correction device 10. The calculated value is calculated based on the heat conduction equations (Equations (1) to (7), and Equation (18)) constituting the revised slab temperature model, and this calculated value is used as the predicted value.

次いで、ステップS12において、炉温制御装置20の炉温制御部22は、ステップS11で予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、加熱炉2内における炉温を制御する(炉温制御ステップ)。 Next, in step S12, the furnace temperature control unit 22 of the furnace temperature control device 20 determines the longitudinal position corresponding to the heating furnace skid portion at the width center of the slab S at the time of extraction in the heating furnace 2 predicted in step S11. The furnace temperature in the heating furnace 2 is controlled so that the temperature in the thickness direction at is higher than the minimum extraction temperature (furnace temperature control step).

具体的に述べると、炉温制御部22は、ステップS11で予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度(例えば、1200℃)を上回る目標温度となるように、各燃料流量制御部23を制御して各燃焼バーナー3への燃料流量を制御し、加熱炉2内の炉温を制御する。 Specifically, the furnace temperature control unit 22 controls the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S at the time of extraction of the heating furnace 2 predicted in step S11. Each fuel flow control unit 23 is controlled to control the fuel flow rate to each combustion burner 3 so that the temperature reaches a target temperature exceeding the minimum extraction temperature (for example, 1200° C.), and the furnace temperature in the heating furnace 2 is controlled. Control.

これにより、加熱炉2から抽出された時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回り、昇温不足となりがちなスキッドビームの接触箇所を規定温度まで上昇させて抽出することができることになる。 As a result, when the slab S is extracted from the heating furnace 2, the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S exceeds the minimum extraction temperature, which tends to result in insufficient temperature rise. This means that the contact point of the skid beam can be raised to a specified temperature for extraction.

このように、本実施形態に係る加熱炉の炉温制御方法によれば、スラブ温度予測ステップ(ステップS11)において、スラブ温度モデルの修正方法によって修正されたスラブ温度モデルを用いて加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測し、炉温制御ステップ(ステップS12)において、スラブ温度予測ステップで予測された加熱炉2の抽出時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、加熱炉2内における炉温を制御する。 As described above, according to the furnace temperature control method for a heating furnace according to the present embodiment, in the slab temperature prediction step (step S11), the slab temperature model corrected by the slab temperature model correction method is used to control the temperature of the heating furnace 2. The temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction is predicted, and in the furnace temperature control step (step S12), the temperature distribution of the slab S at the time of extraction of the heating furnace 2 predicted in the slab temperature prediction step is calculated. The furnace temperature in the heating furnace 2 is controlled so that the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid portion in the width center of S exceeds the minimum extraction temperature.

これにより、加熱炉2から抽出された時のスラブSの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回り、昇温不足となりがちなスキッドビームの接触箇所を規定温度まで上昇させて抽出することができ、加熱炉スキッド部の昇温不足による熱延コイルの材質不良が防止できるとともに、過加熱操業による加熱原単位を削減することができる。 As a result, when the slab S is extracted from the heating furnace 2, the temperature in the thickness direction at the longitudinal position corresponding to the heating furnace skid part in the width center of the slab S exceeds the minimum extraction temperature, which tends to result in insufficient temperature rise. It is possible to raise the contact point of the skid beam to a specified temperature and extract it, preventing material defects in the hot-rolled coil due to insufficient temperature rise of the skid part of the heating furnace, and reducing the heating unit consumption due to overheating operation. can.

また、本発明の一実施形態に係る鋼板の製造方法は、前述の加熱炉の炉温制御方法によって加熱炉2内の炉温を制御する炉温制御工程を含み、この炉温制御工程によって炉温が制御された加熱炉2によって加熱され、加熱炉2から抽出されたスラブSは、複数台(本実施形態ではn台)の粗圧延機4~4による粗圧延工程で粗圧延された後、仕上圧延工程で図示しない仕上圧延機によって仕上圧延され、薄い鋼板となる。そして、この鋼板は、冷却工程で冷却設備5にて冷却された後、巻取工程で巻取り装置6によってコイル状に巻き取られて製造される。 Further, the method for manufacturing a steel plate according to an embodiment of the present invention includes a furnace temperature control step of controlling the furnace temperature in the heating furnace 2 by the above-mentioned method for controlling the furnace temperature of the heating furnace. The slab S that is heated by the temperature-controlled heating furnace 2 and extracted from the heating furnace 2 is rough rolled in a rough rolling process by a plurality of (n units in this embodiment) rough rolling mills 4 1 to 4 n . After that, in a finish rolling process, the steel plate is finish rolled by a finish rolling mill (not shown) to become a thin steel plate. This steel plate is manufactured by being cooled in a cooling facility 5 in a cooling process, and then wound into a coil shape by a winding device 6 in a winding process.

以上、本発明の実施形態について説明してきたが、本発明はこれに限定されずに種々の変更、改良を行うことができる。 Although the embodiments of the present invention have been described above, the present invention is not limited thereto, and various changes and improvements can be made.

例えば、抽出温度算出部12(ステップS2)において、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出しているが、スラブSの少なくとも幅中央部における長さ方向と厚さ方向を含む温度分布を算出すればよく、スラブSの幅中央部における長さ方向及び厚さ方向に限らず、例えば、スラブSの幅方向、長さ方向及び厚さ方向の温度分布を算出するようにしてもよい。 For example, the extraction temperature calculation unit 12 (step S2) calculates the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2. It is sufficient to calculate the temperature distribution including the length direction and the thickness direction in the width center of the slab S. The temperature distribution in the horizontal direction may also be calculated.

また、粗出側温度算出部13(ステップS3)において、粗圧延後のスラブSの幅中央部における長さ方向と厚さ方向の温度分布を算出しているが、粗圧延後のスラブSの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を算出すればよく、スラブSの幅中央部における長さ方向及び厚さ方向に限らず、例えば、スラブSの幅方向、長さ方向及び厚さ方向の温度分布を算出するようにしてもよい。 In addition, in the rough-rolling side temperature calculation unit 13 (step S3), the temperature distribution in the length direction and thickness direction at the width center part of the slab S after rough rolling is calculated. It is sufficient to calculate the temperature distribution including the length direction and the thickness direction at least in the width center part, and is not limited to the length direction and thickness direction in the width center part of the slab S. The temperature distribution in the width direction and the thickness direction may be calculated.

また、スラブ温度予測部21(ステップS11)において、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を予測しているが、加熱炉2の抽出時のスラブSの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を予測すればよく、スラブSの幅中央部における長さ方向及び厚さ方向に限らず、例えば、スラブSの幅方向、長さ方向及び厚さ方向の温度分布を予測するようにしてもよい。算出するようにしてもよい。 In addition, the slab temperature prediction unit 21 (step S11) predicts the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2. What is necessary is to predict the temperature distribution including the length direction and the thickness direction at least at the width center of the slab S, and is not limited to the length direction and thickness direction at the width center of the slab S. The temperature distribution in the direction, length direction, and thickness direction may be predicted. It may be calculated.

本発明の効果を検証すべく、図1に示す熱間圧延ライン1において、粗圧延機4~4の出側に設置された粗温度測定計8によって粗圧延後のスラブSの裏面Sbの幅中央部における長さ方向の温度分布を測定した。 In order to verify the effects of the present invention, in the hot rolling line 1 shown in FIG . The temperature distribution in the length direction at the center of the width was measured.

その結果、図10に示すように、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tact_skidと、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tactとの温度差Tact_diff=Tact-Tact_skidが26℃であった。 As a result, as shown in FIG. 10, the temperature Tact_skid at the longitudinal position corresponding to the heating furnace skid portion at the width center of the back surface Sb of the slab S, and the heating furnace skid at the width center portion of the back surface Sb of the slab S. The temperature difference T act_diff = T act - T act_skid between the temperature T act at the corresponding longitudinal position in between was 26°C.

一方、加熱炉スキッド部での抜熱量Qskdの補正係数ηを0.20(修正前の補正係数)とした熱伝導方程式に基いて、修正装置10の抽出温度算出部12が、加熱炉2に装入されるスラブSの装入温度を初期値として、スラブSの加熱炉在炉中の操業条件を用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出し、出力した。そして、修正装置10の粗出側温度算出部13が、その出力された、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を初期値として、加熱炉2から抽出されたスラブSを粗圧延する際の操業条件を用いて、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を、熱伝導方程式に基づき算出し、出力した。 On the other hand, based on the heat conduction equation in which the correction coefficient η of the amount of heat removed Q skd at the heating furnace skid portion is 0.20 (correction coefficient before correction), the extraction temperature calculation unit 12 of the correction device 10 calculates the Using the charging temperature of the slab S charged in the heating furnace 2 as an initial value, and using the operating conditions of the slab S while it is in the heating furnace, the length direction and thickness at the width center of the slab S at the time of extraction in the heating furnace 2 are determined. The temperature distribution in the horizontal direction was calculated and output. Then, the rough-out side temperature calculation unit 13 of the correction device 10 uses the outputted temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction in the heating furnace 2 as an initial value, and performs heating. Using the operating conditions when rough rolling the slab S extracted from the furnace 2, the temperature distribution in the length direction and thickness direction at the width center of the rough rolled slab S is calculated based on the heat conduction equation. , outputted.

その結果、図11の破線で示すように、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skidと、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcalとの温度差Tcal_diff=Tcal-Tcal_skidが19℃であった。 As a result, as shown by the broken line in FIG. The temperature difference between the temperature T cal at the corresponding longitudinal position between the furnace skids, T cal_diff = T cal - T cal_skid , was 19°C.

このため、修正装置10のモデル修正部14は、温度差Tcal_diff=Tcal-Tcal_skid=19℃と、温度差Tact_diff=Tact-Tact_skid=26℃との偏差e=Tcal_diff-Tact_diff=-7℃を算出した。 Therefore, the model correction unit 14 of the correction device 10 calculates the deviation e=T cal_diff -T between the temperature difference T cal_diff = T cal - T cal_skid = 19°C and the temperature difference T act_diff = T act - T act_skid = 26° C . act_diff =-7°C was calculated.

修正装置10のモデル修正部14は、この算出した偏差e=-7℃(絶対値の大きさ)が閾値ε(=1℃)未満か否かを判定し、その結果、判定結果がNOであったため、補正係数η=0.20を微小変化後の補正係数0.20+0.05=0.25に設定した。 The model correction unit 14 of the correction device 10 judges whether the calculated deviation e=-7°C (absolute value magnitude) is less than the threshold value ε (=1°C), and as a result, the judgment result is NO. Therefore, the correction coefficient η=0.20 was set to the correction coefficient 0.20+0.05=0.25 after the slight change.

そして、モデル修正部14は、η=0.20を微小変化後の補正係数0.20+0.05=0.25とした加熱炉スキッド部での抜熱量Qskdを用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出し直し、そして、これに基づき、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布の算出をし直し、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skid’と、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcal’との温度差Tcal_diff’=Tcal’-Tcal_skid’を算出した。この結果、温度差Tcal_diff’は22℃となった。 Then, the model correction unit 14 extracts the heating furnace 2 using the amount of heat removed Q skd at the heating furnace skid portion with η = 0.20 as the correction coefficient 0.20 + 0.05 = 0.25 after the slight change. The temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of rough rolling is recalculated, and based on this, the temperature distribution in the length direction and thickness direction at the width center of the slab S after rough rolling is calculated. The distribution is recalculated, and the temperature T cal_skid ' at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S and the temperature T cal_skid ' at the length direction position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S is calculated. The temperature difference T cal_diff ′=T cal ′−T cal_skid ′ with respect to the temperature T cal ′ at the corresponding longitudinal position was calculated. As a result, the temperature difference T cal_diff ' was 22°C.

そして、モデル修正部14は、温度差Tcal_diff’=22℃と、温度差Tcal_diff=19℃と、補正係数の微小変化量δ=0.05とから、影響係数∂e/∂η=(Tcal_diff’-Tcal_diff)/δ=(22-19/0.05=60を算出した。 Then, the model correction unit 14 calculates the influence coefficient ∂e/∂η = ( T cal_diff ′−T cal_diff )/δ=(22−19/0.05=60) was calculated.

次いで、モデル修正部14は、加熱炉スキッド部での抜熱量Qskdにおける補正係数ηを更新後補正係数η-e/(∂e/∂η)=0.20-(-7)/60=0.31に更新した。 Next, the model correction unit 14 updates the correction coefficient η for the amount of heat removed Q skd in the heating furnace skid part by updating the correction coefficient η-e/(∂e/∂η)=0.20-(-7)/60= Updated to 0.31.

次いで、モデル修正部14は、補正係数ηを更新後補正係数η-e/(∂e/∂η)=0.31とした加熱炉スキッド部での抜熱量Qskdを用いて、加熱炉2の抽出時のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布を算出し直し、粗圧延後のスラブSの幅中央部における長さ方向及び厚さ方向の温度分布の算出をし直し、スラブSの裏面Sbの幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度Tcal_skid’’と、スラブSの裏面Sbの幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度Tcal’’との温度差Tcal_diff’’=Tcal’’-Tcal_skid’’を算出した。この結果、図11の実線で示すように、温度差Tcal_diff’’=26℃となった。 Next, the model correction unit 14 adjusts the heating furnace 2 by using the amount of heat removed Q skd at the heating furnace skid portion with the correction coefficient η updated as the correction coefficient η−e/(∂e/∂η)=0.31. Recalculate the temperature distribution in the length direction and thickness direction at the width center of the slab S at the time of extraction, and calculate the temperature distribution in the length direction and thickness direction at the width center of the slab S after rough rolling. The temperature T cal_skid '' at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface Sb of the slab S, and the length corresponding to the length between the heating furnace skids at the width center of the back surface Sb of the slab S. The temperature difference T cal_diff '' from the temperature T cal '' at the horizontal position was calculated: T cal '' - T cal_skid ''. As a result, as shown by the solid line in FIG. 11, the temperature difference T cal_diff '' was 26°C.

そして、モデル修正部14は、前述の偏差e=-7℃を偏差e=Tcal_diff’’-Tact_diff=26℃-26℃=0℃に更新した。
そして、モデル修正部14は、更新した偏差e=0℃が閾値ε(=1℃)未満か否かを判定し、判定結果がYESとなったので、判定結果がYESとなったときの更新後補正係数η-e/(∂e/∂η)=0.20-(-7)/60=0.31を補正係数ηとした。
これにより、モデル修正部14での処理が終了し、加熱炉スキッド部での抜熱量Qskd((18)式)における適切な補正係数η=0.31が求められ、加熱炉スキッド部での抜熱量Qskdの適切な修正となった。
Then, the model correction unit 14 updated the aforementioned deviation e=-7°C to the deviation e=T cal_diff ''-T act_diff =26°C-26°C=0°C.
Then, the model correction unit 14 determines whether the updated deviation e=0°C is less than the threshold value ε (=1°C), and since the determination result is YES, update when the determination result is YES. The post-correction coefficient η-e/(∂e/∂η)=0.20-(-7)/60=0.31 was set as the correction coefficient η.
As a result, the processing in the model correction unit 14 is completed, and an appropriate correction coefficient η = 0.31 for the amount of heat removed from the heating furnace skid Q skd (Equation (18)) is obtained, and the This resulted in an appropriate correction for the amount of heat removed Qskd .

従って、本実施例によれば、加熱炉スキッド部での抜熱量Qskdを粗温度測定値と適合するように調整できることが確認された。
また、本発明例に係る鋼板の製造方法によって製造された熱延コイルの長さ方向の材質は、本発明法を適用しない場合に比較して均一になり、本発明の効果が確認された。
Therefore, according to this example, it was confirmed that the amount of heat removed Q skd in the heating furnace skid part could be adjusted to match the crude temperature measurement value.
Further, the material quality in the longitudinal direction of the hot rolled coil manufactured by the method for manufacturing a steel plate according to the example of the present invention was more uniform than that in the case where the method of the present invention was not applied, confirming the effect of the present invention.

1 熱間圧延ライン
2 加熱炉
~2 燃焼制御帯
3 燃焼バーナー
~4 粗圧延機
5 冷却設備
6 巻取り装置
7 装入温度測定計
8 粗温度測定計
10 スラブ温度モデルの修正装置
11 粗温度取得部
12 抽出温度算出部
13 粗出側温度算出部
14 モデル修正部
15 出力部
20 加熱炉の炉温制御装置
21 スラブ温度予測部
22 炉温制御部
23 燃料流量制御部
S スラブ
Sa 表面
Sb 裏面
1 Hot rolling line 2 Heating furnace 2 1 ~ 2 n combustion control zone 3 Combustion burner 4 1 ~ 4 n rough rolling mill 5 Cooling equipment 6 Winding device 7 Charging temperature meter 8 Rough temperature meter 10 Slab temperature model Correction device 11 Crude temperature acquisition unit 12 Extraction temperature calculation unit 13 Crude side temperature calculation unit 14 Model correction unit 15 Output unit 20 Furnace temperature control device for heating furnace 21 Slab temperature prediction unit 22 Furnace temperature control unit 23 Fuel flow rate control unit S Slab Sa front side Sb back side

Claims (5)

スラブを加熱炉で加熱し、該加熱炉から抽出された前記スラブを粗圧延機で粗圧延する際における、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布の予測に用いられるスラブ温度モデルの修正装置であって、
前記粗圧延機の出側に設置された粗温度測定計によって測定した粗圧延後の前記スラブの裏面の幅中央部における長さ方向の温度分布を取得する粗温度取得部と、
前記加熱炉に装入される前記スラブの装入温度を初期値として、前記スラブの加熱炉在炉中の操業条件を用いて、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、前記スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する抽出温度算出部と、
該抽出温度算出部から出力された、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を初期値として、前記加熱炉から抽出された前記スラブを粗圧延する際の操業条件を用いて、粗圧延後の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、熱伝導方程式に基づき算出し、出力する粗出側温度算出部と、
該粗出側温度算出部から出力された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗出側温度算出部から出力された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差と、前記粗温度取得部で取得された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗温度取得部で取得された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差、の偏差の大きさが閾値未満となるように、前記抽出温度算出部で用いられる、前記スラブ温度モデルを構成する熱伝導方程式に含まれる加熱炉スキッド部での抜熱量を修正するモデル修正部とを備えていることを特徴とするスラブ温度モデルの修正装置。
When a slab is heated in a heating furnace and the slab extracted from the heating furnace is roughly rolled in a rough rolling mill, the length direction and the thickness direction of the slab at least in the width center portion at the time of extraction from the heating furnace. A slab temperature model correction device used for predicting temperature distribution including
a rough temperature acquisition unit that acquires the temperature distribution in the length direction at the width center portion of the back surface of the slab after rough rolling, measured by a rough temperature measuring meter installed on the exit side of the rough rolling mill;
Using the charging temperature of the slab charged into the heating furnace as an initial value, and using the operating conditions of the slab while it is in the heating furnace, the length of the slab at least in the center width at the time of extraction from the heating furnace is determined. an extraction temperature calculation unit that calculates and outputs a temperature distribution including a width direction and a thickness direction based on a heat conduction equation forming the slab temperature model;
The temperature distribution extracted from the heating furnace, with the temperature distribution including the length direction and thickness direction in at least the width center portion of the slab at the time of extraction from the heating furnace, which is output from the extraction temperature calculation unit, as an initial value. Using the operating conditions when rough rolling the slab, calculate the temperature distribution including the length direction and the thickness direction at least at the center width of the slab after rough rolling based on the heat conduction equation, and output the rough rolling temperature distribution. An outlet temperature calculation section;
The temperature at the longitudinal position corresponding to the heating furnace skid portion at the width center portion of the back surface of the slab, which is output from the rough output side temperature calculation unit, and the temperature of the slab, which is output from the rough output side temperature calculation unit. The temperature difference between the temperature at the longitudinal position corresponding to between the heating furnace skids at the width center of the back surface of the slab and the temperature at the heating furnace skid portion at the width center of the back surface of the slab, which is obtained by the rough temperature acquisition section. the temperature difference between the temperature at the corresponding longitudinal position and the temperature at the longitudinal position corresponding between the heating furnace skids at the width center portion of the back surface of the slab, which is acquired by the rough temperature acquisition unit; and a model correction unit that corrects the amount of heat removed in the heating furnace skid part included in the heat conduction equation constituting the slab temperature model used in the extraction temperature calculation unit so that the size is less than a threshold value. A slab temperature model correction device characterized by:
請求項1に記載のスラブ温度モデルの修正装置によって修正されたスラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を予測するスラブ温度予測部と、該スラブ温度予測部で予測された加熱炉の抽出時のスラブの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、前記加熱炉内における炉温を制御する炉温制御部とを備えていることを特徴とする加熱炉の炉温制御装置。 A slab whose temperature distribution including the length direction and thickness direction at least in the width center portion of the slab at the time of extraction in a heating furnace is predicted using the slab temperature model corrected by the slab temperature model correction device according to claim 1. The temperature prediction unit and the temperature in the thickness direction at a longitudinal position corresponding to the heating furnace skid portion in the width center of the slab at the time of extraction of the heating furnace predicted by the slab temperature prediction unit exceed the minimum extraction temperature. A furnace temperature control device for a heating furnace, comprising a furnace temperature control section that controls the furnace temperature in the heating furnace. スラブを加熱炉で加熱し、該加熱炉から抽出された前記スラブを粗圧延機で粗圧延する際における、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布の予測に用いられるスラブ温度モデルの修正方法であって、
前記粗圧延機の出側に設置された粗温度測定計によって測定した粗圧延後の前記スラブの裏面の幅中央部における長さ方向の温度分布を取得する粗温度取得ステップと、
前記加熱炉に装入される前記スラブの装入温度を初期値として、前記スラブの加熱炉在炉中の操業条件を用いて、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、前記スラブ温度モデルを構成する熱伝導方程式に基づき算出し、出力する抽出温度算出ステップと、
該抽出温度算出ステップで出力された、前記加熱炉の抽出時の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を初期値として、前記加熱炉から抽出された前記スラブを粗圧延する際の操業条件を用いて、粗圧延後の前記スラブの少なくとも幅中央部における長さ方向と厚さ方向とを含む温度分布を、熱伝導方程式に基づき算出し、出力する粗出側温度算出ステップと、
該粗出側温度算出ステップで出力された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗出側温度算出ステップで出力された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差と、前記粗温度取得ステップで取得された、前記スラブの裏面の幅中央部における加熱炉スキッド部に対応する長さ方向位置の温度と、前記粗温度取得ステップで取得された、前記スラブの裏面の幅中央部における加熱炉スキッド間に対応する長さ方向位置の温度との温度差、の偏差の大きさが閾値未満となるように、前記抽出温度算出ステップで用いられる、前記スラブ温度モデルを構成する熱伝導方程式に含まれる加熱炉スキッド部での抜熱量を修正するモデル修正ステップとを含むことを特徴とするスラブ温度モデルの修正方法。
When a slab is heated in a heating furnace and the slab extracted from the heating furnace is roughly rolled in a rough rolling mill, the length direction and the thickness direction of the slab at least in the width center portion at the time of extraction from the heating furnace. A method for modifying a slab temperature model used for predicting temperature distribution including
a rough temperature acquisition step of acquiring a temperature distribution in the length direction at the width center portion of the back surface of the slab after rough rolling, measured by a rough temperature measuring meter installed on the exit side of the rough rolling mill;
Using the charging temperature of the slab charged into the heating furnace as an initial value, and using the operating conditions of the slab while it is in the heating furnace, the length of the slab at least in the center width at the time of extraction from the heating furnace is determined. an extraction temperature calculation step of calculating and outputting a temperature distribution including a width direction and a thickness direction based on a heat conduction equation constituting the slab temperature model;
The temperature distribution extracted from the heating furnace, with the temperature distribution including the length direction and thickness direction in at least the width center portion of the slab at the time of extraction from the heating furnace, which is output in the extraction temperature calculation step, as an initial value. Using the operating conditions when rough rolling the slab, calculate the temperature distribution including the length direction and the thickness direction at least at the center width of the slab after rough rolling based on the heat conduction equation, and output the rough rolling temperature distribution. Outlet temperature calculation step;
The temperature at the longitudinal position corresponding to the heating furnace skid part at the width center of the back surface of the slab, which is output in the rough output side temperature calculation step, and the temperature of the slab, which is output in the rough output side temperature calculation step. The temperature difference between the temperature at the longitudinal position corresponding to between the heating furnace skids at the width center of the back surface of the slab and the temperature at the heating furnace skid portion at the width center of the back surface of the slab obtained in the rough temperature acquisition step. The temperature difference between the temperature at the corresponding longitudinal position and the temperature at the longitudinal position corresponding between the heating furnace skids at the center width of the back surface of the slab, which is obtained in the rough temperature acquisition step. and a model modification step of modifying the amount of heat removed from the heating furnace skid part included in the heat conduction equation constituting the slab temperature model used in the extraction temperature calculation step so that the magnitude is less than a threshold value. A method for modifying a slab temperature model featuring .
請求項3に記載のスラブ温度モデルの修正方法によって修正されたスラブ温度モデルを用いて加熱炉の抽出時のスラブの少なくとも幅中央部における長さ方向及び厚さ方向を含む温度分布を予測するスラブ温度予測ステップと、該スラブ温度予測ステップで予測された加熱炉の抽出時のスラブの幅中央部における加熱炉スキッド部に対応する長さ方向位置での厚さ方向の温度が最低抽出温度を上回るように、前記加熱炉内における炉温を制御する炉温制御ステップとを含むことを特徴とする加熱炉の炉温制御方法。 A slab whose temperature distribution including the length direction and thickness direction at least in the width center portion of the slab at the time of extraction in a heating furnace is predicted using the slab temperature model modified by the slab temperature model modification method according to claim 3. The temperature prediction step and the temperature in the thickness direction at a longitudinal position corresponding to the heating furnace skid portion in the width center of the slab during extraction of the heating furnace predicted in the slab temperature prediction step exceed the minimum extraction temperature. A furnace temperature control method for a heating furnace, comprising a furnace temperature control step of controlling the furnace temperature in the heating furnace. 請求項4に記載の加熱炉の炉温制御方法によって加熱炉内の炉温を制御する炉温制御工程を含むことを特徴とする鋼板の製造方法。 A method for manufacturing a steel sheet, comprising a furnace temperature control step of controlling a furnace temperature in a heating furnace by the furnace temperature control method for a heating furnace according to claim 4.
JP2022029241A 2022-02-28 2022-02-28 Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate Pending JP2023125249A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022029241A JP2023125249A (en) 2022-02-28 2022-02-28 Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022029241A JP2023125249A (en) 2022-02-28 2022-02-28 Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate

Publications (1)

Publication Number Publication Date
JP2023125249A true JP2023125249A (en) 2023-09-07

Family

ID=87887731

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022029241A Pending JP2023125249A (en) 2022-02-28 2022-02-28 Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate

Country Status (1)

Country Link
JP (1) JP2023125249A (en)

Similar Documents

Publication Publication Date Title
JPS6111289B2 (en)
JP4598586B2 (en) Cooling control method, apparatus, and computer program
JP2012040593A (en) Device for controlling finishing temperature in hot rolling
JP6035817B2 (en) Automatic combustion control method and apparatus for continuous heating furnace
JP6021450B2 (en) Heating furnace operation support system
CN113423517B (en) Cooling control method and cooling control device for thick steel plate, and manufacturing method for thick steel plate
JP4998655B2 (en) Combustion control method for continuous heating furnace
JP2023125249A (en) Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate
JP2023125248A (en) Correction device and correction method for slab temperature model, furnace temperature control device and furnace temperature control method for heating furnace and manufacturing method of steel plate
JP5493993B2 (en) Thick steel plate cooling control device, cooling control method, and manufacturing method
JP4598580B2 (en) Cooling control method, apparatus, and computer program
JP3710572B2 (en) Heating furnace control device
JP4349177B2 (en) Steel extraction temperature prediction method for continuous heating furnace
JPH0663039B2 (en) Temperature control device for heating furnace
JP2010001548A (en) Method for deciding arrangement interval of steel in heating furnace
JP4258341B2 (en) Manufacturing method of high-strength steel sheet with excellent material uniformity in the longitudinal direction of the steel sheet
JP2010247234A (en) Method, device and computer program for controlling cooling
JP2010167503A (en) Method, device and computer program for controlling cooling
JP2587528B2 (en) Estimating method of charging temperature of heating slab to heating furnace
JPH0759722B2 (en) Induction heating control method during subsequent induction heating of a slab previously gas-heated
JP3504118B2 (en) Slab baking control device for tunnel furnace
JPH06306453A (en) Method for eliminating skid mark in continuous heating furnace
JP6237492B2 (en) Heating furnace heating control device, combustion control method, rolled material manufacturing method
JP2011195932A (en) Method for heating billet
JP3007107B2 (en) Material heating curve determination method for heating furnace

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230926