JP6662109B2 - Product state prediction device, product state control device, product state prediction method and program - Google Patents

Product state prediction device, product state control device, product state prediction method and program Download PDF

Info

Publication number
JP6662109B2
JP6662109B2 JP2016041358A JP2016041358A JP6662109B2 JP 6662109 B2 JP6662109 B2 JP 6662109B2 JP 2016041358 A JP2016041358 A JP 2016041358A JP 2016041358 A JP2016041358 A JP 2016041358A JP 6662109 B2 JP6662109 B2 JP 6662109B2
Authority
JP
Japan
Prior art keywords
product
state
predicted
calculating
predicted value
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.)
Active
Application number
JP2016041358A
Other languages
Japanese (ja)
Other versions
JP2017157094A (en
Inventor
小林 俊介
俊介 小林
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.)
Nippon Steel Corp
Original Assignee
Nippon 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 Nippon Steel Corp filed Critical Nippon Steel Corp
Priority to JP2016041358A priority Critical patent/JP6662109B2/en
Publication of JP2017157094A publication Critical patent/JP2017157094A/en
Application granted granted Critical
Publication of JP6662109B2 publication Critical patent/JP6662109B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Control Of Heat Treatment Processes (AREA)
  • Heat Treatment Of Strip Materials And Filament Materials (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)
  • Control Of Metal Rolling (AREA)

Description

本発明は、例えば鋼板の冷却工程において、温度予測モデルにより鋼板の冷却停止温度の予測値を算出するのに利用して好適な製品の状態予測装置、製品の状態制御装置、製品の状態予測方法及びプログラムに関する。   The present invention relates to a product state prediction device, a product state control device, and a product state prediction method suitable for use in calculating a predicted value of a cooling stop temperature of a steel plate by using a temperature prediction model in a steel plate cooling step, for example. And programs.

鋼板の熱間圧延工程の下流の冷却工程では、所定の温度まで鋼板を冷却することで製品に必要な機械特性を得る。
鋼板の冷却制御は、製品の機械特性に直結し、また、歩留まりに影響を及ぼすため、その精度向上が求められる。冷却制御では、水冷による鋼板の温度変化を伝熱モデル計算により推定し、所望の冷却停止温度となるように、冷却水量や鋼板の搬送速度等を決定する。伝熱モデルは熱流体力学や伝熱工学の知見に則り作成されるが、鋼板の表面性状や設備の経時変化等、モデル化が難しいファクターも多く、伝熱モデルによる制御だけでは冷却停止温度を目標値に精度良く一致させることは困難である。
In the cooling step downstream of the hot rolling step of the steel sheet, the steel sheet is cooled to a predetermined temperature to obtain the mechanical properties required for the product.
Since the cooling control of the steel sheet is directly connected to the mechanical properties of the product and affects the yield, it is required to improve the accuracy. In the cooling control, a temperature change of the steel sheet due to water cooling is estimated by a heat transfer model calculation, and a cooling water amount, a transport speed of the steel sheet, and the like are determined so as to obtain a desired cooling stop temperature. Although the heat transfer model is created based on the knowledge of thermo-fluid dynamics and heat transfer engineering, there are many factors that are difficult to model, such as the surface properties of steel sheets and changes over time in equipment. It is difficult to accurately match the target value.

そこで、冷却停止温度のモデル計算値と実績値との差(モデル化誤差)を学習モデルによって予測し、冷却停止温度のモデル計算値を補償する方法が採用されている。
特許文献1には、冷却工程に供する当該厚鋼板について、冷却工程における厚鋼板の温度変化挙動を予測するための厚鋼板の温度予測モデルを用いて、冷却工程における当該厚鋼板の冷却停止温度の予測値を算出する予測値算出工程、スラブ毎に過去の実績データを蓄積したデータベースから、当該厚鋼板と製造条件が類似する厚鋼板の過去の実績データを抽出する抽出工程、前記抽出工程において抽出した前記過去の実績データに基づいて、線形回帰モデル式をたて、当該厚鋼板の冷却停止温度の予測値の誤差を推定する推定工程、前記予測値算出工程において算出した当該厚鋼板の冷却停止温度の予測値と、前記推定工程において推定した当該厚鋼板の冷却停止温度の予測値の誤差とから、前記冷却停止温度の修正値を算出する修正値算出工程、を有することが開示されている。
Therefore, a method is employed in which a difference (modeling error) between the model calculated value of the cooling stop temperature and the actual value is predicted by a learning model, and the model calculated value of the cooling stop temperature is compensated.
Patent Literature 1 discloses that, for a steel plate to be subjected to a cooling process, a temperature prediction model of the steel plate for predicting a temperature change behavior of the steel plate in the cooling process is used to calculate a cooling stop temperature of the steel plate in the cooling process. A predictive value calculating step of calculating a predictive value, an extracting step of extracting past actual data of a thick steel sheet having similar manufacturing conditions to the thick steel sheet from a database storing past actual data for each slab, and extracting in the extracting step. An estimation step of estimating an error of a predicted value of a cooling stop temperature of the thick steel plate based on the past actual data, and estimating an error of a predicted value of a cooling stop temperature of the thick steel plate. A correction value calculation for calculating a correction value of the cooling stop temperature from a predicted value of the temperature and an error of the predicted value of the cooling stop temperature of the thick steel plate estimated in the estimation step. Step, it is disclosed that has a.

特許第5682484号公報Japanese Patent No. 5682484

International Journal of Pharmaceutics, Vol.421, No.2, pp.269-274,2011International Journal of Pharmaceutics, Vol.421, No.2, pp.269-274,2011 Applied Intelligence, Vol.17, Issue 1, pp.49-60,2002Applied Intelligence, Vol.17, Issue 1, pp.49-60,2002 鉄と鋼, 15, 2509-2514, 1981Iron and steel, 15, 2509-2514, 1981

上述したようにデータベースから実績データを抽出するときに、製造条件の類似の判定は、例えば式(1)の重み付きユークリッド距離関数に基づいて行われる。ここで、xi:=(xi,1,・・・,xi,mTは製造条件、qjは製造条件係数、iは鋼板(製造No)を示す添字(i=0:予測対象材)である。ユークリッド距離diが小さい順に実績データを所定の数だけ抽出する。 As described above, when the performance data is extracted from the database, the similarity of the manufacturing conditions is determined based on, for example, the weighted Euclidean distance function of Expression (1). Here, x i : = (x i, 1 ,..., X i, m ) T is a manufacturing condition, q j is a manufacturing condition coefficient, and i is a suffix (i = 0: prediction) indicating a steel plate (manufacturing number). Target material). Euclidean distance d i actual data in ascending order of a predetermined number extracted.

Figure 0006662109
Figure 0006662109

このようにして抽出される実績データを近傍教師データとして局所回帰モデルを生成する。例えば式(2)、式(3)で表わされる局所重み付き回帰モデルを生成する。ここで、βは回帰パラメータ、yはモデル化誤差実績値、y^(^はyの上に付されているものとする)はモデル化誤差予測値、nは近傍教師データ数である。また、wiは重み係数である。 A local regression model is generated using the performance data extracted in this way as neighborhood teacher data. For example, a local weighted regression model represented by Expressions (2) and (3) is generated. Here, β is a regression parameter, y is a modeled error actual value, y ^ (^ is attached to y) is a modeled error predicted value, and n is the number of neighboring teacher data. W i is a weight coefficient.

Figure 0006662109
Figure 0006662109

具体的な手段として、局所重み付き線形重回帰(式(4)、式(5))、局所重み付きPLS(Partial Least Squares)回帰(式(6)〜式(8))が提案されている。   As specific means, local weighted linear multiple regression (Equations (4) and (5)) and local weighted PLS (Partial Least Squares) regression (Equations (6) to (8)) have been proposed. .

Figure 0006662109
Figure 0006662109

Figure 0006662109
Figure 0006662109

局所重み付き線形重回帰は、モデル化誤差予測値y^を製造条件xの線形結合で表わし(式(4))、重み付き誤差2乗和を最小化するような係数aregを求める方法である(式(5))。この方法では、製造条件同士に相関がある場合、回帰係数が不安定になる多重共線性の問題がある。
そこで、局所重み付きPLS回帰では、多重共線性が発生しないよう、変数同士の共分散が最大化される潜在変数zを生成する(式(8))。Ctrは行列XTWyyTWXの固有ベクトルからなる行列である。X:=(x1,・・・,xnT、W:=diag(w1,・・・,wn)、y:=(y1,・・・,ynTである。その後、潜在変数zに対して線形重回帰を行うことにより(式(6)、式(7))、多重共線性の問題を回避している。
The local weighted linear multiple regression is a method in which a modeled error prediction value y ^ is represented by a linear combination of manufacturing conditions x (Equation (4)), and a coefficient a reg that minimizes the weighted error sum of squares is obtained. (Equation (5)). This method has a problem of multicollinearity in which regression coefficients become unstable when there is a correlation between manufacturing conditions.
Therefore, in the local weighted PLS regression, a latent variable z that maximizes the covariance between variables is generated so that multicollinearity does not occur (Equation (8)). C tr is a matrix composed of eigenvectors of the matrix X T Wy T WX. X: = (x 1, ··· , x n) T, W: = diag (w 1, ···, w n), y: = (y 1, ···, y n) is T. Thereafter, the problem of multicollinearity is avoided by performing linear multiple regression on the latent variable z (Equations (6) and (7)).

以上の局所重み付き回帰の手順において、重み係数wiの設定では、一般的に、式(9)の形が用いられる。αはパラメータであり、正の実数である。 In the above-described procedure of local weighted regression, the setting of the weighting coefficient w i generally uses the form of Expression (9). α is a parameter and is a positive real number.

Figure 0006662109
Figure 0006662109

特許文献1の他にも、局所重み付き回帰を利用する非特許文献1や非特許文献2でも同様の形が用いられている。
これは、ユークリッド距離diが0〜∞をとる実数であることに対し、重み係数wiは1〜0となり、要求点(ここでは予測対象材)に近いほど重視し、かつ大きくなりすぎないことが局所重み付き回帰手法と親和性が高いためと考えられる。
In addition to Patent Literature 1, similar forms are used in Non-Patent Literatures 1 and 2 that use local weighted regression.
This is because the Euclidean distance d i is a real number taking 0 to ∞, the weight coefficient w i is 1-0, and the request point (here the prediction target material) focused closer to, and not too large This is thought to be due to the high affinity with the local weighted regression method.

しかしながら、予測対象材が特殊な鋼板であったり、新規の鋼種であったりすると、類似の製造条件を持つ実績データが近傍にないことがある。図5は、近傍教師データの分布と重み係数との関係を示す図である。類似の製造条件を持つ実績データが近傍にない場合、図5に示すようにすべてのiについてwi≒0となり、回帰計算が発散して、予測不可能となる。
一方、回帰計算が発散しないようにする対処法として、式(9)のパラメータαを小さくして、図6に示すように、wiが1→0となる勾配を緩やかにすることが考えられる。しかしながら、逆に近傍に実績データが多くある予測対象材の場合、すべてのiについてwi≒1となり、局所重み付き回帰の意義が失われ、予測精度の低下が懸念される。
However, when the material to be predicted is a special steel plate or a new steel type, there may be no actual data having similar manufacturing conditions in the vicinity. FIG. 5 is a diagram illustrating the relationship between the distribution of the neighborhood teacher data and the weighting factor. If there is no actual data having similar manufacturing conditions in the vicinity, as shown in FIG. 5, w i ≒ 0 for all i, the regression calculation diverges, and prediction becomes impossible.
On the other hand, as a remedy the regression calculation to prevent divergence, by decreasing the parameter α of the formula (9), as shown in FIG. 6, it is conceivable to gentle gradients w i is 1 → 0 . However, conversely, in the case of a prediction target material having a lot of actual data in the vicinity, w i ≒ 1 for all i, the significance of the locally weighted regression is lost, and there is a concern that the prediction accuracy may be reduced.

鋼板の冷却工程では様々な条件の鋼板を冷却するため、近傍の実績データがどのような分布となっているかを事前に知ることは困難である。そのため、一律のパラメータαで重み係数を管理することには無理がある。
非特許文献1では、α(σd)=α´/σdとしてdiのばらつきσdの関数としているが、diの平均μdと比較して非常にσdが小さい場合や大きい場合は依然として上記の問題が発生する懸念がある。
In the steel plate cooling process, since the steel plates under various conditions are cooled, it is difficult to know in advance what distribution the result data in the vicinity has. Therefore, it is impossible to manage the weight coefficient with the uniform parameter α.
In Non-Patent Document 1, α (σ d) = α'/ σ but d is a function of the variation sigma d of d i as much sigma d or when large or small compared to the average mu d of d i If There is still concern that the above problems may occur.

本発明は上記のような点に鑑みてなされたものであり、回帰計算が発散せず、局所重み付き回帰の意義を保ったまま高い精度で、鋼板の冷却停止温度等の製品の状態のモデル化誤差を予測できるようにすることを目的とする。   The present invention has been made in view of the above points, the regression calculation does not diverge, high accuracy while maintaining the significance of local weighted regression, a model of the state of the product such as the cooling stop temperature of the steel sheet It is an object of the present invention to be able to predict a conversion error.

上記の課題を解決するための本発明の要旨は、以下のとおりである。
[1] 製品の製造工程において、状態予測モデルにより製品の状態の予測値を算出する製品の状態予測装置であって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出する抽出手段と、
前記抽出手段で抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算する計算手段と、
前記計算手段で求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償する補償手段とを備え、
前記計算手段は、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出手段で抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、式(101)により設定することを特徴とする製品の状態予測装置。
The gist of the present invention for solving the above problems is as follows.
[1] A product state prediction device that calculates a predicted value of a product state using a state prediction model in a product manufacturing process,
When calculating a predicted value of the state of the product to be predicted, extracting means for extracting neighborhood teacher data based on a distance function from a database storing actual data,
Using the neighborhood teacher data extracted by the extraction means, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the local weighted regression model Calculating means for calculating the error;
Compensation means for compensating for the predicted value of the state of the product to be predicted calculated by the state prediction model, based on the error of the predicted value of the state of the product to be predicted determined by the calculation means,
Said computing means, the weight coefficient w i of the local weighted regression model (i subscripts indicating the product), the distance d i by the distance function near the teacher data extracted by the extraction means, the distance d i by using the average value mu d and standard deviation sigma d, the state prediction device products and setting the equation (101).

Figure 0006662109
Figure 0006662109

[2] 前記データベースは、実績データとして、製品の状態実績値と前記状態予測モデルを用いて算出した製品の状態の予測値との差を製造条件と紐付けて蓄積し、
前記抽出手段は、前記予測対象製品の製造条件と類似する製造条件を持つ実績データを、前記距離関数に基づいて近傍教師データとして抽出することを特徴とする[1]に記載の製品の状態予測装置。
[3] 前記製品は鋼板、前記製造工程は冷却工程、前記製品の状態は鋼板の冷却停止温度、前記状態予測モデルは温度予測モデルであることを特徴とする[1]又は[2]に記載の製品の状態予測装置。
[4] [1]又は[2]に記載の製品の状態予測装置と、
前記補償手段で補償した前記予測対象製品の状態の予測値が、予め前記予測対象製品毎に定められた目標値と一致するように、前記製造工程に用いられる製造設備の操作量を制御する制御手段とを備えたことを特徴とする製品の状態制御装置。
[5] 前記製品は鋼板、前記製造工程は冷却工程、前記製品の状態は鋼板の冷却停止温度、前記状態予測モデルは温度予測モデル、前記操作量は冷却水量及び鋼板の搬送速度のうち少なくともいずれかであることを特徴とする[4]に記載の製品の状態制御装置。
[6] 製品の製造工程において、状態予測モデルを用いて製品の状態の予測値を算出する製品の状態予測方法であって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出するステップと、
前記抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算するステップと、
前記求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償するステップとを有し、
前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求めるときに、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、前記式(101)により設定することを特徴とする製品の状態予測方法。
[7] 製品の製造工程において、状態予測モデルを用いて製品の状態の予測値を算出するためのプログラムであって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出する抽出手段と、
前記抽出手段で抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算する計算手段と、
前記計算手段で求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償する補償手段としてコンピュータを機能させ、
前記計算手段は、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出手段で抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、前記式(101)により設定することを特徴とするプログラム。
[2] The database accumulates, as performance data, a difference between a product status performance value and a product status prediction value calculated using the status prediction model in association with manufacturing conditions,
The product state prediction according to [1], wherein the extraction unit extracts actual data having manufacturing conditions similar to those of the product to be predicted as neighborhood teacher data based on the distance function. apparatus.
[3] The product according to [1] or [2], wherein the product is a steel plate, the manufacturing process is a cooling process, the state of the product is a cooling stop temperature of the steel plate, and the state prediction model is a temperature prediction model. Product state prediction device.
[4] A product state prediction device according to [1] or [2],
Control for controlling the operation amount of the manufacturing equipment used in the manufacturing process such that the predicted value of the state of the prediction target product compensated by the compensation unit matches a target value predetermined for each of the prediction target products. And a means for controlling the state of the product.
[5] The product is a steel plate, the manufacturing process is a cooling process, the state of the product is a cooling stop temperature of the steel plate, the state prediction model is a temperature prediction model, and the operation amount is at least one of a cooling water amount and a conveyance speed of the steel plate. The product state control device according to [4], wherein:
[6] A product state prediction method for calculating a predicted value of a product state using a state prediction model in a product manufacturing process,
When calculating a predicted value of the state of the product to be predicted, extracting neighborhood teacher data based on a distance function from a database in which actual data is accumulated;
Using the extracted neighborhood teacher data, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the error is calculated by the local weighted regression model. Steps to
Compensating for the predicted value of the state of the product to be predicted calculated by the state prediction model, based on the error of the predicted value of the state of the product to be predicted calculated,
When calculating the error of the predicted value of the state of the product to be predicted calculated by the state prediction model, the weighting coefficient w i (i is a subscript indicating the product) of the locally weighted regression model is calculated by using the extracted neighborhood teacher data. wherein the distance d i by the distance function, the distance mean value of d i with the mu d and standard deviation sigma d, the state estimation method of the product and setting by the formula (101) of.
[7] A program for calculating a predicted value of a product state using a state prediction model in a product manufacturing process,
When calculating a predicted value of the state of the product to be predicted, extracting means for extracting neighborhood teacher data based on a distance function from a database storing actual data,
Using the neighboring teacher data extracted by the extracting means, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the local weighted regression model Calculating means for calculating the error;
Based on the error of the predicted value of the state of the product to be predicted determined by the calculating means, the computer functions as compensation means for compensating the predicted value of the state of the product to be predicted calculated by the state prediction model,
Said computing means, the weight coefficient w i of the local weighted regression model (i subscripts indicating the product), the distance d i by the distance function near the teacher data extracted by the extraction means, the distance d i by using the average value mu d and standard deviation sigma d, programs and setting by the formula (101).

本発明によれば、近傍教師データの分布に合わせて重み係数が減少するように設定されるので、回帰計算が発散せず、局所重み付き回帰の意義を保ったまま高い精度で、鋼板の冷却停止温度等の製品の状態のモデル化誤差を予測することができる。   According to the present invention, since the weight coefficient is set so as to decrease in accordance with the distribution of the neighborhood teacher data, the regression calculation does not diverge, and the cooling of the steel sheet is performed with high accuracy while maintaining the significance of the locally weighted regression. It is possible to predict a modeling error of a product state such as a stop temperature.

実施形態の冷却制御システムの構成例を示す図である。It is a figure showing the example of composition of the cooling control system of an embodiment. 局所回帰モデル生成部の機能構成を示すブロック図である。FIG. 4 is a block diagram illustrating a functional configuration of a local regression model generation unit. 実施形態における近傍教師データの分布と重み係数との関係を示す図である。It is a figure showing distribution of neighborhood teacher data and a relation with a weight coefficient in an embodiment. 数値実験のための機能構成を示す図である。FIG. 3 is a diagram showing a functional configuration for a numerical experiment. 従来の近傍教師データの分布と重み係数との関係を示す図である。It is a figure which shows the relationship between the distribution of conventional neighborhood teacher data, and a weight coefficient. 従来の近傍教師データの分布と重み係数との関係を示す図である。It is a figure which shows the relationship between the distribution of conventional neighborhood teacher data, and a weight coefficient.

以下、添付図面を参照して、本発明の好適な実施形態について説明する。本実施形態では、本発明を冷却工程における鋼板の冷却停止温度の予測及び制御に適用した例を述べる。
図1は、冷却設備を含む冷却制御システムの構成例を示す図である。
冷却設備では、搬送テーブル3により鋼板1が搬送され、冷却ノズル2から水を噴射して鋼板1を冷却する。
冷却ノズル2は、冷却プロセスコンピュータ5から出力される指示値に基づいて、冷却水量を制御する。冷却水量の制御方式としては、流量弁開度調節、ノズルON/OFF制御等、一般的な冷却水量の制御方式が適用可能である。
また、搬送テーブル3は、冷却プロセスコンピュータ5から出力される指示値に基づいて、鋼板の搬送速度を制御する。鋼板の搬送速度の制御方式としては、モータドライバによる制御等、一般的な搬送速度制御方式が適用可能である。
冷却設備の出側において、温度計4により鋼板の冷却停止温度が測定され、冷却プロセスコンピュータ5に入力される。温度測定方式としては、放射温度計やサーモグラフィ等、一般的な温度測定方式が適用可能である。
Hereinafter, preferred embodiments of the present invention will be described with reference to the accompanying drawings. In this embodiment, an example in which the present invention is applied to prediction and control of a cooling stop temperature of a steel sheet in a cooling process will be described.
FIG. 1 is a diagram illustrating a configuration example of a cooling control system including a cooling facility.
In the cooling equipment, the steel plate 1 is transported by the transport table 3, and water is injected from the cooling nozzle 2 to cool the steel plate 1.
The cooling nozzle 2 controls the amount of cooling water based on an instruction value output from the cooling process computer 5. As a control method of the cooling water amount, a general cooling water amount control method such as flow rate valve opening adjustment, nozzle ON / OFF control and the like can be applied.
Further, the transfer table 3 controls the transfer speed of the steel sheet based on the instruction value output from the cooling process computer 5. As a control method of the transfer speed of the steel sheet, a general transfer speed control method such as control by a motor driver can be applied.
At the outlet side of the cooling equipment, the temperature at which the steel plate is stopped to cool is measured by the thermometer 4 and input to the cooling process computer 5. As the temperature measuring method, a general temperature measuring method such as a radiation thermometer or a thermography can be applied.

冷却プロセスコンピュータ5は、冷却停止温度計算部7、減算器8、データベース9、局所回帰モデル生成部10、及び冷却設定計算部11を備える。冷却プロセスコンピュータ5には、上位の圧延プロセスコンピュータ6から予測対象材の製造条件が入力される。製造条件としては、例えば鋼板サイズ、鋼種、冷却停止温度の目標値等があるが、外気温、圧延ロール粗度等、冷却工程に影響する情報として一般的である要素は含み得る。本実施形態では、冷却プロセスコンピュータ5が本発明でいう鋼板の冷却停止温度予測装置、鋼板の冷却制御装置として機能する。   The cooling process computer 5 includes a cooling stop temperature calculator 7, a subtractor 8, a database 9, a local regression model generator 10, and a cooling setting calculator 11. The cooling process computer 5 receives the production conditions of the material to be predicted from the upper rolling process computer 6. The manufacturing conditions include, for example, a steel sheet size, a steel type, a target value of a cooling stop temperature, and the like, but may include elements that are general as information affecting the cooling step, such as an outside air temperature and a rolling roll roughness. In the present embodiment, the cooling process computer 5 functions as the steel plate cooling stop temperature prediction device and the steel plate cooling control device according to the present invention.

冷却停止温度計算部7は、鋼板の製造条件と、冷却設備における水量実績値及び鋼板の搬送速度実績値に基づいて、伝熱モデルからなる温度予測モデルにより鋼板の冷却停止温度の計算値(「冷却停止温度計算値」と呼ぶ)を算出する。
減算器8は、温度計4で測定される冷却停止温度実績値と冷却停止温度計算部7から出力される冷却停止温度計算値との差(「モデル化誤差実績値」と呼ぶ)をとる。
The cooling stop temperature calculation unit 7 calculates a cooling stop temperature of the steel sheet by a temperature prediction model including a heat transfer model based on the manufacturing conditions of the steel sheet, the actual value of the water volume in the cooling facility, and the actual value of the transport speed of the steel sheet (“ "Calculated cooling stop temperature").
The subtracter 8 calculates the difference between the actual cooling stop temperature measured by the thermometer 4 and the calculated cooling stop temperature output from the cooling stop temperature calculator 7 (referred to as a “modeling error actual value”).

データベース9は、実績データとして、減算器8から出力されるモデル化誤差実績値を製造条件と紐付けて蓄積する。また、データベース9は、予測対象材の製造条件が入力されると、その製造条件と類似する製造条件を持つ実績データを近傍教師データとして抽出し、局所回帰モデル生成部10に出力する機能を有する。本実施形態では、データベース9が本発明でいう抽出手段として機能する。
データベースから実績データを抽出するときに、製造条件の類似の判定は、例えば式(1)の重み付きユークリッド距離関数に基づいて行われる。なお、本実施形態ではユークリッド距離を例とするが、マハラノビス距離等、多変数系の距離の定義として公知であるものは適用可能である。
The database 9 accumulates the modeled error actual value output from the subtractor 8 in association with the manufacturing condition as actual data. Further, the database 9 has a function of, when the production conditions of the material to be predicted are input, extracting actual data having production conditions similar to the production conditions as neighborhood teacher data, and outputting the data to the local regression model generation unit 10. . In the present embodiment, the database 9 functions as an extraction unit according to the present invention.
When the performance data is extracted from the database, the similarity of the manufacturing conditions is determined based on, for example, the weighted Euclidean distance function of Expression (1). In the present embodiment, the Euclidean distance is taken as an example, but a known variable such as a Mahalanobis distance can be applied to the definition of a multivariable distance.

局所回帰モデル生成部10は、データベース9から抽出される近傍教師データを用いて、温度予測モデルにより算出する予測対象材の冷却停止温度の予測値の誤差(モデル化誤差予測値と呼ぶ)を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルによりモデル化誤差予測値を計算する。本実施形態では、局所回帰モデル生成部10が本発明でいう計算手段として機能する。   The local regression model generation unit 10 obtains an error (referred to as a modeling error prediction value) of the predicted value of the cooling stop temperature of the prediction target material calculated by the temperature prediction model using the neighborhood teacher data extracted from the database 9. A local weighted regression model is generated, and a modeled error prediction value is calculated using the local weighted regression model. In the present embodiment, the local regression model generation unit 10 functions as a calculation unit in the present invention.

冷却設定計算部11は、局所回帰モデル生成部10から出力されるモデル化誤差予測値に基づいて、温度予測モデルにより算出する予測対象材の冷却停止温度の予測値を補償する。そして、冷却設定計算部11は、補償した予測対象材の冷却停止温度の予測値が、予め予測対象材毎に定められた目標値と一致するように、冷却水量及び鋼板の搬送速度のうち少なくともいずれかの指示値を補正した上で、冷却水量及び鋼板の搬送速度の指示値を冷却設備に出力する。本実施形態では、冷却設定計算部11が本発明でいう補償手段、制御手段として機能する。なお、本実施形態では冷却水量と鋼板の搬送速度を操作量とする例としたが、それ以外にも冷却ノズル2の高さ等、冷却設備の操作量として一般的である要素は含み得る。   The cooling setting calculation unit 11 compensates for the predicted value of the cooling stop temperature of the prediction target material calculated by the temperature prediction model, based on the modeling error predicted value output from the local regression model generation unit 10. Then, the cooling setting calculation unit 11 performs at least one of the cooling water amount and the transport speed of the steel sheet such that the compensated predicted value of the cooling stop temperature of the prediction target material matches a target value predetermined for each prediction target material. After correcting any one of the indicated values, the instruction value of the cooling water amount and the conveying speed of the steel sheet is output to the cooling equipment. In the present embodiment, the cooling setting calculation unit 11 functions as a compensation unit and a control unit in the present invention. In the present embodiment, an example in which the amount of cooling water and the transport speed of the steel plate are set as the operation amounts is described. However, other elements that are common as the operation amounts of the cooling equipment, such as the height of the cooling nozzle 2, may be included.

図2は、局所回帰モデル生成部10の機能構成を示すブロック図である。
101は入力部であり、予測対象材の製造条件と、データベース9から抽出される近傍教師データとを入力する。
FIG. 2 is a block diagram illustrating a functional configuration of the local regression model generation unit 10.
Reference numeral 101 denotes an input unit for inputting production conditions of a material to be predicted and neighboring teacher data extracted from the database 9.

102は重み係数設定部であり、式(10)により、局所重み付き回帰モデルの重み係数wiを設定する。ここで、μdは近傍教師データのユークリッド距離diの平均値、σdは近傍教師データのユークリッド距離diの標準偏差である。 Reference numeral 102 denotes a weight coefficient setting unit that sets a weight coefficient w i of a locally weighted regression model by Expression (10). Here, mu d average of the Euclidean distance d i in the vicinity of the teacher data, the sigma d is the standard deviation of the Euclidean distance d i in the vicinity of the teacher data.

Figure 0006662109
Figure 0006662109

103は局所重み付き回帰計算部であり、近傍教師データと重み係数wiに基づいて、局所重み付き線形重回帰(式(4)、式(5))や、局所重み付きPLS(Partial Least Squares)回帰(式(6)〜式(8))等により局所重み付き回帰モデルを生成する。 Reference numeral 103 denotes a local weighted regression calculation unit based on the neighborhood teacher data and the weighting coefficient w i , a local weighted linear multiple regression (equations (4) and (5)), and a locally weighted PLS (Partial Least Squares). ) Generate a locally weighted regression model by regression (equations (6) to (8)) and the like.

104はモデル化誤差計算部であり、局所重み付き回帰計算部103で生成した局所重み付き回帰モデルに予測対象材の製造条件を入力し、予測対象材のモデル化誤差予測値を算出する。   Reference numeral 104 denotes a modeling error calculation unit that inputs the manufacturing conditions of the material to be predicted into the local weighted regression model generated by the local weighted regression calculation unit 103 and calculates a modeled error prediction value of the material to be predicted.

105は出力部であり、モデル化誤差計算部104で算出した予測対象材のモデル化誤差予測値を冷却設定計算部11に出力する。   Reference numeral 105 denotes an output unit, which outputs a modeling error prediction value of the prediction target material calculated by the modeling error calculation unit 104 to the cooling setting calculation unit 11.

ここで、式(10)は、近傍教師データのユークリッド距離di、ユークリッド距離diの平均値μd及び標準偏差σdを変数とする誤差関数erf()を含み、正規分布の累積分布関数を平均値μdに対して左右反転した関数となる。これにより、図3に示すように、近傍教師データがどのような分布となっていても、重み係数wiがその分布に合わせた減少関数となる。このように近傍教師データの分布に合わせて重み係数wiが減少するように設定されるので、回帰計算が発散せず、局所重み付き回帰の意義を保ったまま高い精度で冷却停止温度のモデル化誤差を予測することができる。 Here, Equation (10) includes an error function erf () for the Euclidean distance d i in the vicinity of the teacher data, the average value mu d and standard deviation sigma d of the Euclidean distance d i and variables, the cumulative distribution function of the normal distribution a function that mirror-reversed with respect to the average value mu d. Thus, as shown in FIG. 3, regardless of the distribution of the neighboring teacher data, the weighting coefficient w i becomes a decreasing function adapted to the distribution. As described above, since the weighting coefficient w i is set to decrease in accordance with the distribution of the neighborhood teacher data, the regression calculation does not diverge, and the cooling stop temperature model is maintained with high accuracy while maintaining the significance of the locally weighted regression. The prediction error can be predicted.

本発明は、冷却工程における鋼板の冷却停止温度の予測及び制御に適用できるだけでなく、汎用的に適用可能であって、ある製造設備にて原材料から製品を製造する製造工程において、状態予測モデルにより製品の状態の予測値を算出し、その予測値に基づき製造設備を制御する場合に広く効果を有する。
すなわち、実施形態における鋼板は製品の一例、冷却工程及び設備は製造工程及び設備の一例、鋼板の冷却停止温度は製品の状態の一例、温度予測モデルは状態予測モデルの一例である。そして、本発明を一般的な製造工程に適用する場合、状態制御システムは例えば図1と同様の構成とすればよく、実施形態における冷却停止温度計算部7を状態計算部とし、冷却設定計算部11を設定計算部とし、冷却水量や鋼板の搬送速度の指示値を製造設備の操作量の指示値と置き換えればよい。
製造条件と操作量の実績値から、状態計算部にて製品の状態を状態予測モデルにより予測し、この予測値である状態計算値と状態実績値との誤差であるモデル化誤差実績値をデータベース9に蓄積する。そして、予測対象製品の状態の予測値を算出する際、実績データを蓄積したデータベース9から、距離関数に基づいて近傍教師データを抽出し、その近傍教師データを用いて、局所モデル生成部10にて前記モデル化誤差を求める局所重み付き回帰モデルを生成し、この局所重み付き回帰モデルに製造条件を与えて計算した予測対象製品のモデル化誤差の予測値に基づいて、前記状態予測モデルにより算出した予測対象製品の状態の予測値を補償する。このように補償することで、状態の予測値を実績値に近づけることができる。また、この補償された状態の予測値が予測対象製品毎に定められた目標値と一致するように、設定計算部にて製造設備の操作量の指示値を計算して製造設備を制御することですることで、製品の状態の実績値をより目標値に近づけることができる。
The present invention can be applied not only to the prediction and control of the cooling stop temperature of the steel sheet in the cooling process, but also to a general purpose, and in a manufacturing process of manufacturing a product from raw materials at a certain manufacturing facility, a state prediction model is used. The present invention is widely effective when calculating a predicted value of a product state and controlling a manufacturing facility based on the predicted value.
That is, the steel sheet in the embodiment is an example of a product, the cooling process and equipment are an example of a manufacturing process and equipment, the cooling stop temperature of the steel sheet is an example of a product state, and the temperature prediction model is an example of a state prediction model. When the present invention is applied to a general manufacturing process, the state control system may have, for example, the same configuration as that of FIG. 1. The cooling stop temperature calculator 7 in the embodiment is a state calculator, and the cooling setting calculator is 11 is a setting calculation unit, and the instruction values of the cooling water amount and the transport speed of the steel plate may be replaced with the instruction values of the operation amount of the manufacturing equipment.
The state calculation unit predicts the state of the product from the actual values of the manufacturing conditions and the manipulated variables using a state prediction model, and a modeled error actual value, which is an error between the state calculation value and the actual state value, is a database. 9 is accumulated. Then, when calculating the predicted value of the state of the product to be predicted, the neighborhood teacher data is extracted from the database 9 storing the actual data based on the distance function, and the local model generation unit 10 is used by using the neighborhood teacher data. A local weighted regression model for calculating the modeling error, and calculating the local weighted regression model based on the predicted value of the modeling error of the product to be predicted calculated by giving manufacturing conditions to the local weighted regression model. The predicted value of the predicted state of the product to be predicted is compensated. By compensating in this way, the predicted value of the state can be made closer to the actual value. In addition, the setting calculation unit calculates the operation amount instruction value of the manufacturing equipment and controls the manufacturing equipment so that the predicted value of the compensated state matches the target value determined for each product to be predicted. By doing so, the actual value of the state of the product can be made closer to the target value.

例えば本発明を熱間圧延における板幅の予測及び制御に適用する場合、第1の実施形態における冷却設備はエッジャ、サイジングプレス等の鋼板幅加工設備となる。また、製造条件は鋼板サイズ、鋼種、装置入側鋼板温度等、板幅に影響を与えると考えられる要素となる。また、冷却ノズル2、搬送テーブル3はエッジャならば油圧圧下装置、サイジングプレスならば機械クランク等の板幅制御のためのアクチュエータとなる。また、温度計4はラインセンサ等の板幅測定装置となる。また、冷却プロセスコンピュータ5は圧延プロセスコンピュータ、上位の圧延コンピュータ6は加熱プロセスコンピュータもしくはビジネスコンピュータとなる。また、冷却停止温度計算部7は例えば非特許文献3に記載のモデルに基づき、装置出側板幅を計算する板幅計算部となる。また、冷却設定計算部11は例えば非特許文献3に記載のモデルに基づき、装置出側板幅が所望の値となるようにエッジャやサイジングプレスの開度、速度を設定する幅圧下設定計算部となる。   For example, when the present invention is applied to prediction and control of a sheet width in hot rolling, the cooling equipment in the first embodiment is a steel sheet width processing equipment such as an edger and a sizing press. Further, the manufacturing conditions are factors considered to affect the sheet width, such as the steel sheet size, the steel type, and the temperature of the steel sheet on the apparatus entrance side. In addition, the cooling nozzle 2 and the transfer table 3 are actuators for controlling a plate width such as a hydraulic pressure reduction device for an edger and a mechanical crank for a sizing press. Further, the thermometer 4 is a plate width measuring device such as a line sensor. The cooling process computer 5 is a rolling process computer, and the upper rolling computer 6 is a heating process computer or a business computer. Further, the cooling stop temperature calculating unit 7 is a plate width calculating unit that calculates a device outlet side plate width based on a model described in Non-Patent Document 3, for example. Further, the cooling setting calculation unit 11 is based on a model described in Non-Patent Document 3, for example, and a width reduction setting calculation unit that sets the opening degree and speed of the edger or sizing press so that the apparatus outlet plate width becomes a desired value. Become.

本発明を適用した手法による効果を数値実験により検証した。
図4に、数値実験のための機能構成を示す。パーソナルコンピュータでデータベース9及び局所回帰モデル生成部10を再現して、数値実験を行う。工場で操業実績を蓄積している工場データベースからパーソナルコンピュータに製造条件を入力し、データベース9及び局所回帰モデル生成部10によりモデル化誤差予測値を出力する。
工場データベースにはモデル化誤差実績値も保存されているので、データベース9及び局所回帰モデル生成部10によるモデル化誤差予測値と比較し、その精度を評価した。図1と違い、モデル化誤差予測値は実操業に反映されないが、モデル化誤差を精度良く予測することで、冷却停止温度制御の精度も向上すると考えられる。
The effect of the method to which the present invention was applied was verified by numerical experiments.
FIG. 4 shows a functional configuration for a numerical experiment. The database 9 and the local regression model generation unit 10 are reproduced by a personal computer to perform a numerical experiment. Manufacturing conditions are input to a personal computer from a factory database storing operation results at the factory, and a modeling error prediction value is output by the database 9 and the local regression model generation unit 10.
Since the factory database also stores the modeled error actual values, the factory database was compared with the modeled error predicted values by the database 9 and the local regression model generator 10 to evaluate the accuracy. Unlike FIG. 1, the modeling error prediction value is not reflected in the actual operation, but it is considered that by accurately predicting the modeling error, the accuracy of the cooling stop temperature control is improved.

特許文献1、非特許文献1、及び非特許文献2に記載の内容に基づいて、局所回帰モデル生成部10では局所重み付きPLS回帰を用いるものとし、重み係数を1とした場合(比較法1)、式(9)で求めた重み係数とした場合(比較法2)、式(10)で求めた重み係数とした場合(発明法)のそれぞれについてモデル化誤差の予測精度を評価した。
数値実験は6カ月分10115本の操業実績について行い、精度評価には誤差平均とばらつきの両方を評価できる指標である2乗誤差平均RMSEを用いた。表1に、比較法1、2及び発明法の精度評価の結果を示す。
Based on the contents described in Patent Literature 1, Non-Patent Literature 1, and Non-Patent Literature 2, the local regression model generation unit 10 uses the locally weighted PLS regression, and the weight coefficient is set to 1 (Comparative method 1 ), The weighting factor determined by equation (9) (Comparative method 2), and the weighting factor determined by equation (10) (invention method), the prediction accuracy of the modeling error was evaluated.
Numerical experiments were performed on the operation results of 10115 lines for 6 months, and the accuracy evaluation was performed using the square error average RMSE, which is an index capable of evaluating both the error average and the variation. Table 1 shows the results of the accuracy evaluations of the comparative methods 1 and 2 and the inventive method.

Figure 0006662109
Figure 0006662109

表1によれば、比較法1では、重み係数を全て1としたため、予測対象材に近い近傍教師データを重視しながらの回帰ができず、モデル化誤差の予測精度は良くない。比較法2では、回帰計算が発散して、モデル化誤差予測値を計算できない鋼板があった。これらに対し、発明法では、回帰計算が発散することなく、良好な精度でモデル化誤差を予測できている。   According to Table 1, in the comparative method 1, since all the weighting factors are 1, regression cannot be performed while emphasizing the neighborhood teacher data close to the material to be predicted, and the prediction accuracy of the modeling error is not good. In Comparative Method 2, there was a steel sheet whose regression calculation diverged and could not calculate the modeling error prediction value. On the other hand, in the invention method, the modeling error can be predicted with good accuracy without diverging the regression calculation.

本発明を適用した鋼板の冷却停止温度予測装置、鋼板の冷却制御装置は、例えばCPU、ROM、RAM等を備えたコンピュータ装置により実現される。また、図1では、冷却停止温度計算部7、局所回帰モデル生成部10、冷却設定計算部11が1つのプロセスコンピュータ(冷却プロセスコンピュータ5)で動作する例を説明したが、それぞれ別個のコンピュータ装置で動作し、ネットワークを介して入出力を行うような形態でもよい。また、製造条件は圧延プロセスコンピュータ6からの入力となっているが、上位のビジネスコンピュータからの入力となる等、熱間圧延における加速冷却制御の形態として一般的である構成に対して本発明は適用可能である。
また、本発明は、本発明の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータがプログラムを読み出して実行することによっても実現可能である。
The steel plate cooling stop temperature prediction device and the steel plate cooling control device to which the present invention is applied are realized by, for example, a computer device having a CPU, a ROM, a RAM, and the like. FIG. 1 illustrates an example in which the cooling stop temperature calculation unit 7, the local regression model generation unit 10, and the cooling setting calculation unit 11 operate on one process computer (cooling process computer 5). , And input / output via a network. In addition, the present invention is applied to a configuration that is generally used as a form of accelerated cooling control in hot rolling, for example, although the manufacturing conditions are input from the rolling process computer 6 and are input from an upper business computer. Applicable.
The present invention can also be realized by supplying software (program) for realizing the functions of the present invention to a system or apparatus via a network or various storage media, and a computer of the system or apparatus reads and executes the program. It is feasible.

1:鋼板、2:冷却ノズル、3:搬送テーブル、4:温度計、5:冷却プロセスコンピュータ、7:冷却停止温度計算部、8:減算器、9:データベース、10:局所回帰モデル生成部、11:冷却設定計算部、101:入力部、102:重み係数設定部、103:局所重み付き回帰計算部、104:モデル化誤差計算部、105:出力部   1: steel plate, 2: cooling nozzle, 3: transport table, 4: thermometer, 5: cooling process computer, 7: cooling stop temperature calculator, 8: subtractor, 9: database, 10: local regression model generator, 11: cooling setting calculation unit, 101: input unit, 102: weight coefficient setting unit, 103: regression calculation unit with local weight, 104: modeling error calculation unit, 105: output unit

Claims (7)

製品の製造工程において、状態予測モデルにより製品の状態の予測値を算出する製品の状態予測装置であって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出する抽出手段と、
前記抽出手段で抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算する計算手段と、
前記計算手段で求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償する補償手段とを備え、
前記計算手段は、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出手段で抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、式(101)
Figure 0006662109
により設定することを特徴とする製品の状態予測装置。
In a product manufacturing process, a product state prediction device that calculates a predicted value of the state of the product using a state prediction model,
When calculating a predicted value of the state of the product to be predicted, extracting means for extracting neighborhood teacher data based on a distance function from a database storing actual data,
Using the neighborhood teacher data extracted by the extraction means, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the local weighted regression model Calculating means for calculating the error;
Compensation means for compensating for the predicted value of the state of the product to be predicted calculated by the state prediction model, based on the error of the predicted value of the state of the product to be predicted determined by the calculation means,
Said computing means, the weight coefficient w i of the local weighted regression model (i subscripts of an article), the distance d i by the distance function near the teacher data extracted by the extraction means, the distance d i Using the average value μ d and the standard deviation σ d , equation (101)
Figure 0006662109
An apparatus for predicting the state of a product, wherein the apparatus is set by:
前記データベースは、実績データとして、製品の状態実績値と前記状態予測モデルを用いて算出した製品の状態の予測値との差を製造条件と紐付けて蓄積し、
前記抽出手段は、前記予測対象製品の製造条件と類似する製造条件を持つ実績データを、前記距離関数に基づいて近傍教師データとして抽出することを特徴とする請求項1に記載の製品の状態予測装置。
The database, as actual data, accumulates the difference between the actual state value of the product and the predicted value of the state of the product calculated using the state prediction model in association with the manufacturing conditions,
2. The product state prediction according to claim 1, wherein the extraction unit extracts actual data having manufacturing conditions similar to those of the product to be predicted as neighborhood teacher data based on the distance function. 3. apparatus.
前記製品は鋼板、前記製造工程は冷却工程、前記製品の状態は鋼板の冷却停止温度、前記状態予測モデルは温度予測モデルであることを特徴とする請求項1又は2に記載の製品の状態予測装置。   3. The product state prediction according to claim 1, wherein the product is a steel sheet, the manufacturing process is a cooling step, the product state is a cooling stop temperature of the steel sheet, and the state prediction model is a temperature prediction model. 4. apparatus. 請求項1又は2に記載の製品の状態予測装置と、
前記補償手段で補償した前記予測対象製品の状態の予測値が、予め前記予測対象製品毎に定められた目標値と一致するように、前記製造工程に用いられる製造設備の操作量を制御する制御手段とを備えたことを特徴とする製品の状態制御装置。
A product state prediction device according to claim 1 or 2,
Control for controlling the operation amount of the manufacturing equipment used in the manufacturing process such that the predicted value of the state of the prediction target product compensated by the compensation unit matches a target value predetermined for each of the prediction target products. And a means for controlling the state of the product.
前記製品は鋼板、前記製造工程は冷却工程、前記製品の状態は鋼板の冷却停止温度、前記状態予測モデルは温度予測モデル、前記操作量は冷却水量及び鋼板の搬送速度のうち少なくともいずれかであることを特徴とする請求項4に記載の製品の状態制御装置。   The product is a steel sheet, the manufacturing step is a cooling step, the state of the product is a cooling stop temperature of the steel sheet, the state prediction model is a temperature prediction model, and the operation amount is at least one of a cooling water amount and a conveyance speed of the steel sheet. The product state control device according to claim 4, wherein: 製品の製造工程において、状態予測モデルを用いて製品の状態の予測値を算出する製品の状態予測方法であって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出するステップと、
前記抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算するステップと、
前記求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償するステップとを有し、
前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求めるときに、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、式(101)
Figure 0006662109
により設定することを特徴とする製品の状態予測方法。
In a product manufacturing process, a product state prediction method for calculating a predicted value of the state of the product using a state prediction model,
When calculating a predicted value of the state of the product to be predicted, extracting neighborhood teacher data based on a distance function from a database in which actual data is accumulated;
Using the extracted neighborhood teacher data, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the error is calculated by the local weighted regression model. Steps to
Compensating for the predicted value of the state of the product to be predicted calculated by the state prediction model, based on the error of the predicted value of the state of the product to be predicted calculated,
When calculating the error of the predicted value of the state of the product to be predicted calculated by the state prediction model, the weighting coefficient w i (i is a subscript indicating the product) of the locally weighted regression model is calculated by using the extracted neighborhood teacher data. the distance and the distance d i by a function, by using the average value mu d and standard deviation sigma d of the distance d i in the formula (101)
Figure 0006662109
A method for predicting the state of a product, characterized in that:
製品の製造工程において、状態予測モデルを用いて製品の状態の予測値を算出するためのプログラムであって、
予測対象製品の状態の予測値を算出する際に、実績データを蓄積したデータベースから、距離関数に基づいて近傍教師データを抽出する抽出手段と、
前記抽出手段で抽出した近傍教師データを用いて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値の誤差を求める局所重み付き回帰モデルを生成し、当該局所重み付き回帰モデルにより当該誤差を計算する計算手段と、
前記計算手段で求めた前記予測対象製品の状態の予測値の誤差に基づいて、前記状態予測モデルにより算出する前記予測対象製品の状態の予測値を補償する補償手段としてコンピュータを機能させ、
前記計算手段は、前記局所重み付き回帰モデルの重み係数wi(iは製品を示す添字)を、前記抽出手段で抽出した近傍教師データの前記距離関数による距離diと、前記距離diの平均値μd及び標準偏差σdとを用いて、式(101)
Figure 0006662109
により設定することを特徴とするプログラム。
In a product manufacturing process, a program for calculating a predicted value of the state of the product using a state prediction model,
When calculating a predicted value of the state of the product to be predicted, extracting means for extracting neighborhood teacher data based on a distance function from a database storing actual data,
Using the neighborhood teacher data extracted by the extraction means, a local weighted regression model for calculating an error of a predicted value of the state of the product to be predicted calculated by the state prediction model is generated, and the local weighted regression model Calculating means for calculating the error;
Based on the error of the predicted value of the state of the product to be predicted determined by the calculating means, a computer functions as a compensation means for compensating the predicted value of the state of the product to be predicted calculated by the state prediction model,
Said computing means, the weight coefficient w i of the local weighted regression model (i subscripts indicating the product), the distance d i by the distance function near the teacher data extracted by the extraction means, the distance d i Using the average value μ d and the standard deviation σ d , equation (101)
Figure 0006662109
A program characterized by being set by:
JP2016041358A 2016-03-03 2016-03-03 Product state prediction device, product state control device, product state prediction method and program Active JP6662109B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016041358A JP6662109B2 (en) 2016-03-03 2016-03-03 Product state prediction device, product state control device, product state prediction method and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016041358A JP6662109B2 (en) 2016-03-03 2016-03-03 Product state prediction device, product state control device, product state prediction method and program

Publications (2)

Publication Number Publication Date
JP2017157094A JP2017157094A (en) 2017-09-07
JP6662109B2 true JP6662109B2 (en) 2020-03-11

Family

ID=59809821

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016041358A Active JP6662109B2 (en) 2016-03-03 2016-03-03 Product state prediction device, product state control device, product state prediction method and program

Country Status (1)

Country Link
JP (1) JP6662109B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022014876A (en) * 2020-07-07 2022-01-20 Jfeスチール株式会社 Manufacturing specification determination support device and manufacturing specification determination support method for steel plate, computer program, and computer readable recording medium, and method for manufacturing steel plate

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108070710B (en) * 2017-08-29 2019-03-29 东北大学 A kind of steel plate temperature control process for quenching based on roller quenching machine
JP7255330B2 (en) * 2019-04-10 2023-04-11 富士電機株式会社 Information processing device, information processing method, program
JP2022003158A (en) * 2020-06-23 2022-01-11 三菱重工業株式会社 Control method of heat treatment furnace, heat treatment system, and program
JP7472867B2 (en) 2021-07-08 2024-04-23 Jfeスチール株式会社 Method for determining heat treatment conditions for steel, method for manufacturing steel, information processing device, and program for determining heat treatment conditions for steel

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4488964B2 (en) * 2005-06-27 2010-06-23 新日本製鐵株式会社 Process operating state control method and computer program
JP5682484B2 (en) * 2010-09-16 2015-03-11 新日鐵住金株式会社 Thick steel plate cooling control method, cooling control device, and thick steel plate manufacturing method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022014876A (en) * 2020-07-07 2022-01-20 Jfeスチール株式会社 Manufacturing specification determination support device and manufacturing specification determination support method for steel plate, computer program, and computer readable recording medium, and method for manufacturing steel plate
JP7287416B2 (en) 2020-07-07 2023-06-06 Jfeスチール株式会社 Thick steel plate manufacturing specification determination support device, manufacturing specification search method, computer program, computer-readable recording medium, and thick steel plate manufacturing method

Also Published As

Publication number Publication date
JP2017157094A (en) 2017-09-07

Similar Documents

Publication Publication Date Title
JP6662109B2 (en) Product state prediction device, product state control device, product state prediction method and program
TWI554340B (en) Rolling simulation device
JP6881137B2 (en) Product status predictors and methods, manufacturing process control systems, and programs
JP6922675B2 (en) Manufacturing process state predictors, methods and programs, and manufacturing process control systems
JP6834209B2 (en) Product status prediction device, product status control device, product status prediction method and program
Zheng et al. Distributed model predictive control for plant-wide hot-rolled strip laminar cooling process
JP5012660B2 (en) Product quality prediction and control method
JP6439780B2 (en) Magnetic property prediction device and magnetic property control device for electrical steel sheet
JP6146553B1 (en) Steel plate temperature control device and temperature control method
JP6915428B2 (en) Manufacturing process control systems, methods and programs
KR102122143B1 (en) Steel plate temperature control device and temperature control method
Akram et al. Quality monitoring and process adjustment by integrating SPC and APC: a review
Zheng et al. An approach to model building for accelerated cooling process using instance-based learning
JP5962290B2 (en) Steel heat transfer coefficient prediction apparatus and cooling control method
Pian et al. A hybrid soft sensor for measuring hot-rolled strip temperature in the laminar cooling process
JP5195331B2 (en) Quality prediction apparatus, prediction method, program, and computer-readable recording medium in manufacturing process
Prinz et al. Online parameter estimation for adaptive feedforward control of the strip thickness in a hot strip rolling mill
Gao et al. Using transient measurements in iterative steady-state optimizing control
JP7102962B2 (en) Control set value determination device, control set value determination method, and program
JPH04232214A (en) Method and device for controlling cooling of steel plate
JP2010235972A (en) Manufacturing controller and manufacturing method for high tension steel sheet
JP5007630B2 (en) Product quality control method and control device
Song et al. A Digital Twin Model for Automatic Width Control of Hot Rolling Mill
JP2019057077A (en) Prediction device, method, and program of product characteristics, and control system of manufacturing process
Idzik et al. Decision Support for the Optimization of Continuous Processesusing Digital Shadows

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160404

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20160404

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181105

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190904

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190910

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191007

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200114

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200127

R151 Written notification of patent or utility model registration

Ref document number: 6662109

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151