Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7622673B2 - Quality factor analysis device and steelworks operation method - Google Patents
[go: Go Back, main page]

JP7622673B2 - Quality factor analysis device and steelworks operation method - Google Patents

Quality factor analysis device and steelworks operation method Download PDF

Info

Publication number
JP7622673B2
JP7622673B2 JP2022034779A JP2022034779A JP7622673B2 JP 7622673 B2 JP7622673 B2 JP 7622673B2 JP 2022034779 A JP2022034779 A JP 2022034779A JP 2022034779 A JP2022034779 A JP 2022034779A JP 7622673 B2 JP7622673 B2 JP 7622673B2
Authority
JP
Japan
Prior art keywords
factors
quality
factor
causal
probability distribution
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
JP2022034779A
Other languages
Japanese (ja)
Other versions
JP2023130229A (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.)
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 JP2022034779A priority Critical patent/JP7622673B2/en
Publication of JP2023130229A publication Critical patent/JP2023130229A/en
Application granted granted Critical
Publication of JP7622673B2 publication Critical patent/JP7622673B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Description

本開示は、品質要因解析装置及び製鉄所の操業方法に関する。 This disclosure relates to a quality factor analysis device and a method for operating a steelworks.

従来、操業条件が品質に影響する製造プロセスにおいて、操業条件などから製品の品質を予測する様々な手法が提案されている。例えば特許文献1は、操業データから製品品質を予測する装置であって、操業変数の局所領域ごとに操業と品質との関係を表す関数式を有することで、高精度に品質を予測する装置を開示する。 Conventionally, various methods have been proposed for predicting product quality from operational conditions in manufacturing processes where the operational conditions affect the quality. For example, Patent Document 1 discloses a device that predicts product quality from operational data, and predicts quality with high accuracy by having a function formula that expresses the relationship between operation and quality for each local region of the operational variables.

また、例えば特許文献2は、品質異常の原因推定支援方法であって、品質評価指標と製造工程との相関関係と、決定木学習によって選定された温度・圧力・水量の数値情報又は文字情報がオペレータに表示されることによって、品質異常の原因を、オペレータの経験度合いに左右されることなく推定する方法を開示する。 For example, Patent Document 2 discloses a method for supporting estimation of the cause of quality abnormalities, in which the correlation between quality evaluation indexes and the manufacturing process, and numerical or textual information on temperature, pressure, and water volume selected by decision tree learning are displayed to the operator, thereby estimating the cause of the quality abnormality regardless of the operator's level of experience.

また、例えば特許文献3は、複数製造工程に対する異常工程推定方法であって、工程ごとに正常状態における特徴を機械学習によって相関モデルを生成し、正常状態からの乖離を異常度合と判定することによって、異常工程を精度よく推定することができる方法を開示する。 For example, Patent Document 3 discloses a method for estimating abnormal processes for multiple manufacturing processes, which generates a correlation model for the characteristics in a normal state for each process through machine learning, and determines the degree of abnormality as a deviation from the normal state, thereby enabling accurate estimation of abnormal processes.

特開2012-137813号公報JP 2012-137813 A 特開2014-179060号公報JP 2014-179060 A 特開2019-49940号公報JP 2019-49940 A

ここで、上記の従来手法は、品質異常の原因を特定するために、製品品質データと製造条件データから品質と製造条件の相関に基づく関係に着目することで原因推定を行っている。しかしながら、相関関係があれば必ず因果関係があるということにはならない。 Here, the conventional method described above estimates the cause of quality abnormalities by focusing on the relationship based on the correlation between quality and manufacturing conditions from product quality data and manufacturing condition data. However, the existence of a correlation does not necessarily mean that there is a causal relationship.

例えば製造条件と品質の間に交絡がある場合、製造条件及び品質に関する変数間に多重共線性が生じて擬相関が起こり得るため、従来手法では因果関係が無いにも関わらず原因候補と推定することがあった。ここで、「製造条件と品質の間に交絡がある」とは、製造条件の変数と品質の変数の両方に影響を与える外部変数が存在することを意味する。このような場合に、従来手法によって品質異常の要因と推定された製造条件を変更しても、品質不良及び異常状態が改善しない又は悪化してしまうことがある。 For example, when there is confounding between manufacturing conditions and quality, multicollinearity can occur between variables related to manufacturing conditions and quality, resulting in spurious correlation, and conventional methods may presume that a causal relationship is a possible cause even when there is no causal relationship. Here, "confounding between manufacturing conditions and quality" means that there is an external variable that affects both the manufacturing condition variables and the quality variables. In such cases, even if manufacturing conditions that are presumed to be the cause of quality abnormalities by conventional methods are changed, the quality defect or abnormal state may not improve or may even worsen.

以上の問題を解決すべくなされた本開示の目的は、高精度に品質要因を推定可能な品質要因解析装置及び製鉄所の操業方法を提供することにある。ここで、品質要因は、品質異常及び品質向上などの品質の変化の原因となる因子である。 The purpose of this disclosure, which has been made to solve the above problems, is to provide a quality factor analysis device and a steelworks operation method that can estimate quality factors with high accuracy. Here, quality factors are factors that cause quality changes such as quality abnormalities and quality improvements.

本開示の一実施形態に係る品質要因解析装置は、
製品の品質特性を表す品質因子、前記品質特性を得るための製造プロセスの操作量を表す操作因子、及び、前記操作量と前記品質特性との関係に影響し得る中間的な因子を表す中間因子を含む因子を用いて前記品質特性の要因を解析する品質要因解析装置であって、
過去に製造された前記製品における前記品質因子、前記操作因子及び前記中間因子に関する製造データを蓄積するデータベースと、
前記製造データから前記品質因子、前記操作因子及び前記中間因子の値の範囲を表す初期の確率分布を算出する初期化部と、
前記因子の一部に対して計測又は設定によって確定された条件を読み取る条件読取部と、
前記確定された条件、前記初期の確率分布及び前記因子の間の関係に基づいて条件付確率を算出し、前記条件付確率に基づいて前記因子の確率分布を更新する因果推論部と、
前記因子の少なくとも1つについて、前記初期の確率分布及び更新された確率分布を表示させる表示制御部と、を備える。
A quality factor analysis device according to an embodiment of the present disclosure includes:
1. A quality factor analysis device that analyzes factors of quality characteristics using factors including quality factors that represent quality characteristics of a product, operation factors that represent operation amounts of a manufacturing process for obtaining the quality characteristics, and intermediate factors that represent intermediate factors that may affect a relationship between the operation amounts and the quality characteristics, comprising:
a database that accumulates manufacturing data relating to the quality factors, the operation factors, and the intermediate factors of the products that have been manufactured in the past;
an initialization unit that calculates an initial probability distribution representing a range of values of the quality factor, the operation factor, and the intermediate factor from the manufacturing data;
a condition reading unit that reads conditions determined by measurement or setting for some of the factors;
a causal inference unit that calculates a conditional probability based on the determined condition, the initial probability distribution, and a relationship between the factors, and updates the probability distribution of the factors based on the conditional probability;
and a display control unit that displays the initial probability distribution and the updated probability distribution for at least one of the factors.

本開示の一実施形態に係る製鉄所の操業方法は、
前記製品を鋼材とする、上記の品質要因解析装置が更新した前記操作因子の確率分布に基づいて、前記鋼材の製造条件を決定することを含む。
A method for operating a steel mill according to an embodiment of the present disclosure includes:
The method includes determining manufacturing conditions for the steel product based on the probability distribution of the operational factors updated by the quality factor analysis device.

本開示によれば、高精度に品質要因を推定可能な品質要因解析装置及び製鉄所の操業方法を提供することができる。 This disclosure provides a quality factor analysis device and a steelworks operation method that can estimate quality factors with high accuracy.

図1は、本開示の一実施形態に係る品質要因解析装置の概略を示す構成図である。FIG. 1 is a diagram illustrating an outline of a quality factor analysis device according to an embodiment of the present disclosure. 図2は、演算処理部の処理を示すフローチャートである。FIG. 2 is a flowchart showing the processing of the arithmetic processing unit. 図3は、初期化の処理を示すフローチャートである。FIG. 3 is a flowchart showing the initialization process. 図4は、データベースに格納されている製造データの例である。FIG. 4 is an example of the manufacturing data stored in the database. 図5は、因果接続表を例示する図である。FIG. 5 is a diagram illustrating a causal connection table. 図6は、因果ダイアグラムを例示する図である。FIG. 6 is a diagram illustrating a causal diagram. 図7は、因子の確率分布を表す図である。FIG. 7 is a diagram showing the probability distribution of factors. 図8は、条件付確率テーブルを例示する図である。FIG. 8 is a diagram illustrating a conditional probability table. 図9は、実施例における実験データのグラフである。FIG. 9 is a graph of experimental data in the examples. 図10は、実施例における従来手法による結果を示すグラフである。FIG. 10 is a graph showing the results of the conventional method in the examples. 図11は、実施例における因果接続表である。FIG. 11 is a causal connection table in the embodiment. 図12は、実施例における因果ダイアグラムと確率分布を示す図である。FIG. 12 is a diagram showing a causal diagram and a probability distribution in an embodiment. 図13は、実施例における条件付確率テーブルである。FIG. 13 is a conditional probability table in the embodiment. 図14は、実施例における品質要因解析装置の解析結果を示す図である。FIG. 14 is a diagram showing the analysis results of the quality factor analysis device in the embodiment. 図15は、実施例における因果ダイアグラムを示す図である。FIG. 15 is a diagram showing a causal diagram in the embodiment. 図16は、実施例における不良率100%と設定した場合を示す図である。FIG. 16 is a diagram showing a case where the defect rate is set to 100% in the embodiment. 図17は、実施例におけるMn(マンガン濃度)の確率分布を示す図である。FIG. 17 is a diagram showing the probability distribution of Mn (manganese concentration) in the examples. 図18は、実施例における不良率ゼロと設定した場合を示す図である。FIG. 18 is a diagram showing a case where the defect rate is set to zero in the embodiment. 図19は、実施例における不良率ゼロの場合のMnの確率分布を示す図である。FIG. 19 is a diagram showing the probability distribution of Mn when the defective rate is zero in the example. 図20は、製造コストの分布を示すグラフである。FIG. 20 is a graph showing the distribution of manufacturing costs.

以下、図面を参照して本開示の一実施形態に係る品質要因解析装置及び製鉄所の操業方法が説明される。 Below, a quality factor analysis device and a method for operating a steelworks according to one embodiment of the present disclosure will be described with reference to the drawings.

本実施形態に係る品質要因解析装置は、製品の品質特性を表す品質因子、品質特性を得るための製造プロセスの操作量を表す操作因子、及び、操作量と品質特性との関係に影響し得る中間的な因子を表す中間因子を含む因子を用いて品質特性の要因を解析する。品質要因解析装置が対象とする製品は限定されないが、本実施形態において、鉄鋼プロセスによって製造される鋼材を対象として例示する。品質要因解析装置が実行する処理は、製鉄所の操業方法の一部として行われる。 The quality factor analysis device according to this embodiment analyzes the factors of quality characteristics using factors including quality factors that represent the quality characteristics of a product, operation factors that represent the operation amount of the manufacturing process to obtain the quality characteristics, and intermediate factors that represent intermediate factors that can affect the relationship between the operation amount and the quality characteristics. The products that the quality factor analysis device targets are not limited, but in this embodiment, steel materials manufactured by a steelmaking process are used as an example. The processing executed by the quality factor analysis device is carried out as part of the operating method of the steelworks.

(構成)
図1は、本実施形態に係る品質要因解析装置10の構成例を示す図である。ここで、図1には、品質要因解析装置10の他に、品質要因解析装置10の使用者40と、解析結果などを使用者40に表示する表示装置20と、使用者40から入力を受け付ける入力装置30が示されている。
(composition)
Fig. 1 is a diagram showing an example of the configuration of a quality factor analysis device 10 according to this embodiment. In addition to the quality factor analysis device 10, Fig. 1 shows a user 40 of the quality factor analysis device 10, a display device 20 that displays analysis results and the like to the user 40, and an input device 30 that accepts input from the user 40.

品質要因解析装置10は、品質要因解析を実行するための各種演算処理を行う演算処理部101と、品質要因解析プログラム102を記憶する記憶部(ROM103)と、各種演算処理の結果及び演算値などを一時的に記憶する一時記憶部(RAM104)と、バス105と、過去に製造された製品における因子に関する製造データ(製造実績データ)を蓄積するデータベース110と、を備える。 The quality factor analysis device 10 includes a calculation processing unit 101 that performs various calculation processes to execute quality factor analysis, a memory unit (ROM 103) that stores a quality factor analysis program 102, a temporary memory unit (RAM 104) that temporarily stores the results of various calculation processes and calculation values, a bus 105, and a database 110 that accumulates manufacturing data (manufacturing performance data) related to factors in products manufactured in the past.

演算処理部101は、ROM103に記憶された品質要因解析プログラム102を読み込んで実行し、品質要因解析に必要な各種演算処理を実行する。演算処理部101は、1つ以上のプロセッサである。プロセッサは、例えば汎用のプロセッサ又は特定の処理に特化した専用プロセッサであるが、これらに限られず任意のプロセッサとすることができる。例えば品質要因解析装置10がコンピュータであって、演算処理部101はCPUで実現されてよい。 The arithmetic processing unit 101 reads and executes the quality factor analysis program 102 stored in the ROM 103, and executes various arithmetic processing required for quality factor analysis. The arithmetic processing unit 101 is one or more processors. The processor is, for example, a general-purpose processor or a dedicated processor specialized for a specific process, but is not limited to these and can be any processor. For example, the quality factor analysis device 10 may be a computer, and the arithmetic processing unit 101 may be realized by a CPU.

このような演算処理部101は、品質要因解析プログラム102の実行によって機能する機能ブロックとして、初期化部106、条件読取部107、因果推論部108と、表示制御部109と、を備える。 Such a calculation processing unit 101 includes an initialization unit 106, a condition reading unit 107, a causal inference unit 108, and a display control unit 109 as functional blocks that function when the quality factor analysis program 102 is executed.

(処理)
図2は、演算処理部101の処理を示すフローチャートである。演算処理部101は、概要として、初期化(ステップS201)、条件読取(ステップS202)、因果推論(ステップS203)、結果表示(ステップS204)及び終了判定(ステップS205)の処理を実行する。初期化(ステップS201)の処理は初期化部106で実行される。条件読取(ステップS202)の処理は条件読取部107で実行される。因果推論(ステップS203)の処理は因果推論部108で実行される。結果表示(ステップS204)及び終了判定(ステップS205)の処理は表示制御部109で実行される。演算処理部101は、入力装置30を介して使用者40からの開始指示を受け取ると、品質要因解析の処理を開始する。
(process)
2 is a flowchart showing the processing of the arithmetic processing unit 101. In summary, the arithmetic processing unit 101 executes the processes of initialization (step S201), condition reading (step S202), causal inference (step S203), result display (step S204), and end determination (step S205). The initialization process (step S201) is executed by the initialization unit 106. The condition reading process (step S202) is executed by the condition reading unit 107. The causal inference process (step S203) is executed by the causal inference unit 108. The result display (step S204) and end determination (step S205) are executed by the display control unit 109. When the arithmetic processing unit 101 receives a start instruction from the user 40 via the input device 30, it starts the quality factor analysis process.

(初期化)
初期化部106は、図3に記載のフローチャートに基づき処理(ステップS201-1からステップS201-4)を実行し、因子の値の範囲を表す初期の確率分布を算出する。
(Initialization)
The initialization unit 106 executes the process (steps S201-1 to S201-4) based on the flowchart shown in FIG. 3, and calculates an initial probability distribution that represents the range of the factor values.

初期化部106は製造データの読み込みを実行する(ステップS201-1)。初期化部106は、データベース110に格納されている過去に製造された多数の鋼材における加熱条件、圧延条件、冷却条件、素材の化学成分、走査型電子顕微鏡等により撮像し読み取った鋼材のミクロ組織の特徴量、強度などの鋼材の材質要求値、引張試験等による材質特性値などを因子として抽出する。ここで、本実施形態において素材はスラブである。また、過去に製造された多数の鋼材のデータである製造データは、図4に示すように製品ごとの因子をまとめた製造データとして管理されている。行は過去に製造された製品を表す。また、列は因子を表す。 The initialization unit 106 reads the manufacturing data (step S201-1). The initialization unit 106 extracts factors such as the heating conditions, rolling conditions, cooling conditions, chemical components of the material, features of the microstructure of the steel material imaged and read by a scanning electron microscope or the like, material requirements of the steel material such as strength, and material characteristic values obtained by tensile tests or the like for many steel materials manufactured in the past that are stored in the database 110. Here, in this embodiment, the material is a slab. Furthermore, the manufacturing data, which is data on many steel materials manufactured in the past, is managed as manufacturing data that summarizes factors for each product, as shown in FIG. 4. The rows represent products manufactured in the past. Furthermore, the columns represent factors.

上記のように、因子は、製品の品質特性を表す品質因子、品質特性を得るための製造プロセスの操作量を表す操作因子、及び、操作量と品質特性との関係に影響し得る中間的な因子を表す中間因子を含む。例えば、強度などの鋼材の材質要求値は、製品である鋼材の品質特性を表すものであって、品質因子に対応する。また、例えば加熱条件、圧延条件、冷却条件、素材の化学成分などは、鋼材の製造条件であって、操作因子又は中間因子(品質因子以外の因子)に対応する。本実施形態のような製鉄所の操業方法の一部として行われる品質要因解析において、品質要因解析装置10は、品質因子として、鋼材の強度、延性、靱性、表面疵、介在物及び内質の少なくとも1つを表す因子を用いてよい。ここで、内質は例えば中心偏析、ポロシティなどである。また、品質要因解析装置10は、操作因子として、鋼材の製造条件のうち加熱条件、圧延条件、冷却条件及び素材の化学成分の少なくとも1つを用いてよい。ただし、以下において、説明及び図示の便宜上、因子及び因子の値について数字(1、2、3…)及び文字(A、B、C)を用いて抽象的に表現することがある。 As described above, the factors include quality factors that represent the quality characteristics of the product, operation factors that represent the amount of operation of the manufacturing process to obtain the quality characteristics, and intermediate factors that represent intermediate factors that may affect the relationship between the amount of operation and the quality characteristics. For example, the material requirement value of the steel material, such as strength, represents the quality characteristics of the steel material, which is the product, and corresponds to the quality factor. In addition, for example, the heating conditions, rolling conditions, cooling conditions, and chemical components of the material are manufacturing conditions of the steel material, and correspond to the operation factors or intermediate factors (factors other than the quality factor). In the quality factor analysis performed as part of the operating method of the steelworks such as this embodiment, the quality factor analysis device 10 may use factors that represent at least one of the strength, ductility, toughness, surface defects, inclusions, and internal quality of the steel material as the quality factor. Here, the internal quality is, for example, center segregation, porosity, etc. In addition, the quality factor analysis device 10 may use at least one of the heating conditions, rolling conditions, cooling conditions, and chemical components of the material among the manufacturing conditions of the steel material as the operation factor. However, in the following, for ease of explanation and illustration, factors and factor values may be abstractly represented using numbers (1, 2, 3...) and letters (A, B, C).

初期化部106は、因子間の因果関係を表現した因果接続表を読み込む(ステップS201-2)。因果接続表は、例えば使用者40によって予め用意されて、データベース110に格納されている。図5は、因果接続表を例示する図である。因果接続表は、行が原因の因子、列が結果の因子を表し、行列が交差するセルに1が示されていれば因果関係があることを意味する。図5の例において、因子Aから因子B、因子Aから因子C、因子Bから因子Cへの因果関係がある。初期化部106は、読み込んだ因果接続表に基づき因果関係を図示した因果ダイアグラムを生成して、データベース110に記憶させてよい。図6は、因果ダイアグラムを例示する図である。図6の例において、因果ダイアグラムは、因果関係のある因子が矢印で結ばれており、矢印の元に原因となる因子が示され、矢印の先に結果となる因子が示されている。因果ダイアグラムは、表示制御部109によってデータベース110から読みだされ、表示装置20に表示されてよい。 The initialization unit 106 reads in a causal connection table expressing the causal relationships between factors (step S201-2). The causal connection table is prepared in advance by the user 40, for example, and stored in the database 110. FIG. 5 is a diagram illustrating a causal connection table. In the causal connection table, the rows represent causal factors and the columns represent result factors, and if a 1 is displayed in a cell where the rows and columns intersect, it means that there is a causal relationship. In the example of FIG. 5, there is a causal relationship from factor A to factor B, from factor A to factor C, and from factor B to factor C. The initialization unit 106 may generate a causal diagram illustrating the causal relationships based on the read causal connection table and store it in the database 110. FIG. 6 is a diagram illustrating a causal diagram. In the example of FIG. 6, the causal diagram shows causal factors connected by arrows, with the causal factor at the base of the arrow and the resultant factor at the tip of the arrow. The causal diagram may be read from the database 110 by the display control unit 109 and displayed on the display device 20.

初期化部106は、各因子のヒストグラムの計算をする(ステップS201-3)。初期化部106は、各因子の確率分布(発生確率)も算出する。図7は、因子の確率分布を表す図である。例えば図4における因子A、因子B及び因子Cの度数及び確率分布は、図7のような表形式でそれぞれ表現される。ここで、1列目が各因子の値の範囲、2列目が度数、3列目が確率分布を表す。 The initialization unit 106 calculates a histogram for each factor (step S201-3). The initialization unit 106 also calculates the probability distribution (occurrence probability) of each factor. Figure 7 is a diagram showing the probability distribution of factors. For example, the frequencies and probability distributions of factors A, B, and C in Figure 4 are each expressed in a table format as shown in Figure 7. Here, the first column shows the range of values for each factor, the second column shows the frequency, and the third column shows the probability distribution.

初期化部106は、条件付確率を取得する(ステップS201-4)。図8は、条件付確率テーブルを例示する図である。初期化部106は、図4の製造データ、図5又は図6で表現される因果関係及び図7で表現される各因子の確率分布に基づいて条件付確率テーブルを生成することによって、ステップS201-4の処理を実行できる。ここで、条件付確率は、ある因子について条件が満たされる(ある事象が起こる)場合における、別の因子の確率分布を表現したものである。図8の条件付確率テーブルは、因子Aがある事象が起こる因子、因子Bが別の因子であって、因子Aの値が特定の範囲である場合における、因子Bが値a、b又はcをとり得る確率を示す。図5に示したように、因子Aから因子Bへの因果関係が存在している。例えば因子Aが1.2以上1.4未満のとき、因子Bの値がaとなる条件付確率が0.5、bとなる条件付確率が0.5である。また、初期化部106は、データベース110から条件付確率テーブルを読み込むことによって、ステップS201-4の処理を実行してよい。つまり、条件付確率テーブルは初期化部106によって生成されるのでなく、予め準備されてデータベース110に格納されていてよい。このとき、条件付確率テーブルの値は、物理学等の科学理論、実験で得られた経験則、日々の操業から得た人の経験則等をもとに設定されてよい。 The initialization unit 106 acquires the conditional probability (step S201-4). FIG. 8 is a diagram illustrating a conditional probability table. The initialization unit 106 can execute the process of step S201-4 by generating a conditional probability table based on the manufacturing data in FIG. 4, the causal relationships represented in FIG. 5 or FIG. 6, and the probability distribution of each factor represented in FIG. 7. Here, the conditional probability represents the probability distribution of a factor when a condition for a factor is satisfied (a certain event occurs). The conditional probability table in FIG. 8 shows the probability that factor B can take the value a, b, or c when factor A is a factor for which a certain event occurs, factor B is another factor, and the value of factor A is within a specific range. As shown in FIG. 5, there is a causal relationship from factor A to factor B. For example, when factor A is 1.2 or more and less than 1.4, the conditional probability that factor B's value is a is 0.5, and the conditional probability that factor B's value is b is 0.5. Furthermore, the initialization unit 106 may execute the process of step S201-4 by reading the conditional probability table from the database 110. In other words, the conditional probability table may not be generated by the initialization unit 106, but may be prepared in advance and stored in the database 110. In this case, the values in the conditional probability table may be set based on scientific theories such as physics, rules of thumb obtained from experiments, rules of thumb obtained from people's daily operations, etc.

(条件読取)
条件読取部107は、因子の一部に対して計測又は設定によって確定された条件を読み取る。本実施形態において、条件読取部107は、条件読取の処理(図2のステップS202)として、入力装置30から入力された因子の観測情報を読み込む。観測情報(確定情報)は、例えば、因子Aの値が1.0であるというような確定された因子の値についての情報である。条件読取部107は、1つ以上の因子についての観測情報を読み込む。条件読取部107が、全ての因子についての観測情報を読み込む必要はない。ここで、複数の因子についての観測情報が読み込まれた場合に、各因子についてステップS202からステップS204の処理が実行される。
(Condition reading)
The condition reading unit 107 reads conditions that have been determined by measurement or setting for some of the factors. In this embodiment, the condition reading unit 107 reads observation information of the factors input from the input device 30 as the condition reading process (step S202 in FIG. 2). The observation information (determined information) is information about the determined value of a factor, for example, that the value of factor A is 1.0. The condition reading unit 107 reads observation information for one or more factors. It is not necessary for the condition reading unit 107 to read observation information for all factors. Here, when observation information for multiple factors is read, the processes of steps S202 to S204 are executed for each factor.

(因果推論)
因果推論部108は、条件読取部107が読み取った確定された条件、初期化部106によって算出された初期の確率分布及び因子の間の関係に基づいて条件付確率を算出し、条件付確率に基づいて因子の確率分布を更新する。本実施形態において、因果推論部108は、図2のステップS203に対応する処理として、ステップS202で読み込んだ観測情報に基づいて、以下に説明する式に従うベイズ推論によって、S201で算出した確率分布を更新する。
(Causal inference)
The causal inference unit 108 calculates a conditional probability based on the determined condition read by the condition reading unit 107, the initial probability distribution calculated by the initialization unit 106, and the relationship between the factors, and updates the probability distribution of the factors based on the conditional probability. In this embodiment, the causal inference unit 108 updates the probability distribution calculated in S201 by Bayesian inference according to the formula described below based on the observation information read in step S202, as a process corresponding to step S203 in Fig. 2.

因果推論部108は、ある因子(ここでは因子Uと呼ぶ)に対して、設定によって確定された条件(観測情報)が与えられた場合に、下記の式(1)に従うベイズ推論によって、因子Uの原因となる因子(ここでは因子Sと呼ぶ)の確率分布を更新する。更新された因子Sのヒストグラムは、因子Uが観測情報の通りである場合における、因子Uの原因の因子Sの値の推論値を意味する。 When a condition (observation information) determined by a setting is given for a certain factor (here called factor U), the causal inference unit 108 updates the probability distribution of a factor (here called factor S) that causes factor U by Bayesian inference according to the following formula (1). The histogram of the updated factor S represents the inferred value of the value of factor S that causes factor U when factor U is as shown in the observation information.

Figure 0007622673000001
Figure 0007622673000001

ここで、P(S|U=u)は、因子Uの値がuと確定したときの因子Sの事後確率である。P(U=u|S)は因子Sと因子Uの尤度関数である。P(S)は因子Sの事前確率である。また、P(U=u)は因子Uの観測情報である。 Here, P(S|U=u) is the posterior probability of factor S when the value of factor U is determined to be u. P(U=u|S) is the likelihood function of factors S and U. P(S) is the prior probability of factor S. Also, P(U=u) is the observed information of factor U.

また、因果推論部108は、ある因子(ここでは因子Uと呼ぶ)に対して、設定によって確定された条件(観測情報)が与えられた場合に、下記の式(2)に従うベイズ推論によって、因子Uの結果となる因子(ここでは因子Tと呼ぶ)の確率分布を更新する。更新された因子Tのヒストグラムは、因子Uが観測情報の通りである場合における、因子Uの結果の因子Tの値の推論値を意味する。 Furthermore, when a condition (observation information) determined by the settings is given for a certain factor (here called factor U), the causal inference unit 108 updates the probability distribution of a factor (here called factor T) that is the result of factor U by Bayesian inference according to the following formula (2). The histogram of the updated factor T represents the inferred value of the value of factor T that is the result of factor U when factor U is as per the observation information.

Figure 0007622673000002
Figure 0007622673000002

ここで、P(T|U=u)は、因子Uの値がuと確定したときの因子Tの事後確率である。P(T|U)は因子Tと因子Uの尤度関数である。また、P(U=u)は因子Uの観測情報である。 Here, P(T|U=u) is the posterior probability of factor T when the value of factor U is determined to be u. P(T|U) is the likelihood function of factors T and U. Also, P(U=u) is the observed information of factor U.

ここで、因果推論部108は、データベース110から図5又は図6で表現される因果関係を取得することによって、因子Uの原因の因子S及び因子Uの結果の因子Tを特定することができる。因果推論部108は、観測情報が与えられる因子Uに対する原因の因子S及び結果の因子Tの全てについて、上記の式(1)及び式(2)に従って、確率分布を更新する。ここで、因果推論部108は、因子間の因果関係が多段に設定されるような場合に、再帰的に(全ての原因と結果の因子をたどりながら)確率分布を更新する。 The causal inference unit 108 can identify the causal factor S of factor U and the result factor T of factor U by acquiring the causal relationships shown in FIG. 5 or FIG. 6 from the database 110. The causal inference unit 108 updates the probability distribution for all of the causal factors S and result factors T for factor U for which observation information is given, according to the above formulas (1) and (2). Here, the causal inference unit 108 updates the probability distribution recursively (by tracing all cause and result factors) when the causal relationships between factors are set at multiple stages.

(結果表示)
表示制御部109は、因子の少なくとも1つについて、初期の確率分布及び更新された確率分布を表示させる。本実施形態において、表示制御部109は、図2のステップS204に対応する処理として、ステップS203で更新された各因子の確率分布を、因果ダイアグラムとともに、表示装置20に表示させる。つまり、本実施形態において、表示制御部109は、因果ダイアグラムの表示によって因子の間の関係を視覚的に示し、使用者40が容易に因子の間の関係を把握できるようにする。
(Results display)
The display control unit 109 displays the initial probability distribution and the updated probability distribution for at least one of the factors. In this embodiment, the display control unit 109 displays the probability distribution of each factor updated in step S203, together with a causal diagram, on the display device 20 as a process corresponding to step S204 in Fig. 2. That is, in this embodiment, the display control unit 109 visually shows the relationship between the factors by displaying the causal diagram, allowing the user 40 to easily grasp the relationship between the factors.

(終了判定)
表示制御部109は、図2のステップS205に対応する分岐処理として、入力装置30から使用者40の終了指示があった場合に一連の処理を終了させる。そうでない場合に、再び図2のステップS202の処理が行われる。例えば使用者40の終了指示がなく、複数の因子についての観測情報が読み込まれた場合に、再度のステップS202の処理において別の因子についての観測情報が読み込まれる。
(End judgment)
The display control unit 109 ends a series of processes when an end instruction is received from the user 40 via the input device 30, as a branching process corresponding to step S205 in Fig. 2. If not, the process of step S202 in Fig. 2 is performed again. For example, when there is no end instruction from the user 40 and observation information on a plurality of factors is read, observation information on another factor is read in the process of step S202 again.

(操業方法)
品質要因解析装置10は、鋼材の機械的特性(例えば強度など)並びに内部及び表面の性質(例えば疵及び欠陥の有無など)の少なくとも1つを表す品質因子と、鋼材の製造条件である操作因子と、中間因子との間の関係を解析できる。品質要因解析装置10は、解析結果(更新した操作因子の確率分布)を上位システムに出力してよい。上位システムは、例えば鋼材の製造を管理するプロセスコンピュータを備えて構成され、製鉄所の操業方法を実行する。例えば品質因子に対して所望の結果を観測情報として与えた場合に、品質要因解析装置10が原因となる操作因子を特定して、確率分布を更新することができる。上位システムは、所望の品質の鋼材を製造するために、品質要因解析装置10が更新した操作因子の確率分布に基づいて、鋼材の製造条件を決定してよい。
(Operation method)
The quality factor analysis device 10 can analyze the relationship between quality factors representing at least one of the mechanical properties (e.g., strength, etc.) and internal and surface properties (e.g., the presence or absence of flaws and defects, etc.) of steel, operation factors that are manufacturing conditions of steel, and intermediate factors. The quality factor analysis device 10 may output the analysis result (probability distribution of the updated operation factors) to a higher-level system. The higher-level system is configured, for example, with a process computer that manages the manufacturing of steel, and executes the operation method of the steelworks. For example, when a desired result is given for the quality factor as observation information, the quality factor analysis device 10 can identify the causative operation factor and update the probability distribution. In order to manufacture steel of a desired quality, the higher-level system may determine the manufacturing conditions of the steel based on the probability distribution of the operation factors updated by the quality factor analysis device 10.

以下、本実施形態に係る品質要因解析装置10の効果を説明するために、厚板鋼材の製造に適用した実施例が説明される。 Below, an example of application to the production of thick steel plates will be described to explain the effects of the quality factor analysis device 10 according to this embodiment.

(実施例1)
鋼材商品の強度は、鉄鋼ミクロ組織の特徴である炭窒化物の析出形態、固相の種類等によって決まる。鉄鋼ミクロ組織の特徴は、製造の冷却過程おける鋼材温度変化量[℃]に影響される。ここで、様々に条件を振った実験データを製造データとして用いて、冷却時の鋼材温度変化量・析出物の大きさ・強度の関係性を解析する実験が行われた。
Example 1
The strength of steel products is determined by the precipitation form of carbonitrides, the type of solid phase, and other characteristics of the steel microstructure. The characteristics of the steel microstructure are affected by the amount of steel temperature change [℃] during the cooling process in manufacturing. Here, experiments were conducted to analyze the relationship between the amount of steel temperature change during cooling, the size of precipitates, and strength, using experimental data with various conditions as manufacturing data.

様々に条件を振った500サンプルに対して連続する識別番号であるサンプル連番が与えられた。図9は、横軸をサンプル連番として、縦軸をその識別番号を有するサンプルの因子の値とする実験データのグラフである。「A」は温度変化量である。「B」は析出物半径の逆数である。「X」は鋼材の強度である。これらの数値は-1から1の範囲に含まれるように正規化が行われている。強度について基準値(-0.12)以上を示すサンプルが品質合格、すなわち「良」と判定される。図9の例では、サンプルの総数が500であるところ、基準値を下回るサンプルが149個あり、29.8%の不良率であった。 500 samples under various conditions were given consecutive sample serial numbers, which are identification numbers. Figure 9 is a graph of experimental data with the sample serial numbers on the horizontal axis and the factor values of the samples with those identification numbers on the vertical axis. "A" is the amount of temperature change. "B" is the inverse of the precipitate radius. "X" is the strength of the steel. These values are normalized so that they fall within the range of -1 to 1. Samples that exhibit strength equal to or greater than the reference value (-0.12) are deemed to pass the quality test, i.e., "good." In the example of Figure 9, out of a total of 500 samples, 149 were below the reference value, resulting in a defect rate of 29.8%.

従来の相関に基づく方法では、まず、「A」及び「B」が強度である「X」に与える影響を把握するために、「X」を目的変数、「A」及び「B」を説明変数とする線形回帰モデルをつくる。図7のデータに基づいて、X=0.78×A+1.03×Bという回帰式が得られる。従来の方法では、「X」の最小値が-0.67であるため、「X」を1.0だけ上げて合格率(「良」である比率)を高めるように説明変数が調整される。このとき、1.0を「A」の回帰係数0.78で除した値が1.28であるため、「A」の狙い値を1.28だけ上げる。図10は、調整後の(「A」の狙い値を1.28だけあげた)結果である。図10における「X」の平均値は0.04であって、図9における調整前の「X」の平均値である0.05からほとんど変化しない。また、図10において、基準値を下回るサンプルが依然として108個あり、21.6%の不良率であった。 In the conventional correlation-based method, first, in order to grasp the influence of "A" and "B" on "X", which is the strength, a linear regression model is created with "X" as the objective variable and "A" and "B" as explanatory variables. Based on the data in FIG. 7, the regression equation X = 0.78 x A + 1.03 x B is obtained. In the conventional method, since the minimum value of "X" is -0.67, the explanatory variables are adjusted to increase the pass rate (the rate of "good"). In this case, since 1.0 divided by the regression coefficient of "A", 0.78, is 1.28, the target value of "A" is increased by 1.28. FIG. 10 shows the result after adjustment (the target value of "A" is increased by 1.28). The average value of "X" in FIG. 10 is 0.04, which is almost unchanged from the average value of "X" before adjustment, 0.05, in FIG. 9. Also, in Figure 10, there were still 108 samples that fell below the standard value, resulting in a defect rate of 21.6%.

従来の相関に基づく方法で意図通りに不良率を改善できない理由は、因子Bが因子Aから影響を受けるという現象を考慮できないためと考えられる。つまり、因子Aの変更が因子Bに影響しており、因子Aの変更による因子X(強度)の改善が、影響された因子Bによる因子Xへの改悪で打ち消され、結果として「X」が改善されていないと考えられる。 The reason why conventional correlation-based methods cannot improve the defect rate as intended is thought to be because they fail to take into account the phenomenon in which factor B is affected by factor A. In other words, a change in factor A affects factor B, and the improvement in factor X (strength) resulting from the change in factor A is cancelled out by the deterioration of factor X due to the affected factor B, resulting in no improvement in "X."

同じ例において品質要因解析装置10による要因解析が行われた。初期化部106は、図11に示される因果接続表と図9のデータを読み込み、因子A、因子B、因子X、因子Yに対する因果ダイアグラムと各因子の確率分布、条件付確率テーブルを算出した。図12は、本実施例における因果ダイアグラムと確率分布を示す図である。図13は本実施例における条件付確率テーブルである。ここで、因子Yは良・不良の判定を表す因子である。 In the same example, a factor analysis was performed by the quality factor analysis device 10. The initialization unit 106 read the causal connection table shown in FIG. 11 and the data in FIG. 9, and calculated a causal diagram for factors A, B, X, and Y, a probability distribution for each factor, and a conditional probability table. FIG. 12 shows the causal diagram and probability distribution in this embodiment. FIG. 13 shows the conditional probability table in this embodiment. Here, factor Y is a factor that represents a pass/fail judgment.

次に、条件読取部107は「因子Yにおける不良率が0%」という観測情報を取得する。因果推論部108は、ベイズ推論に基づいて「因子Yにおける不良率が0%」を観測情報としたときの各因子の確率分布を更新する。表示制御部109は、図14のように、更新の前後の確率分布を表示させる。図14によれば、不良率を0%にするためには、因子Aの平均値を0.15あげるだけでなく、因子Aの分布の幅を狭めて、かつ、因子Bの平均値を0.2下げる必要があることが視覚的に示されている。ここで、図14の「OK」が「良」に対応し、「NG」が「不良」に対応する。 Next, the condition reading unit 107 acquires the observation information that "the defect rate in factor Y is 0%." The causal inference unit 108 updates the probability distribution of each factor when "the defect rate in factor Y is 0%" is the observation information based on Bayesian inference. The display control unit 109 displays the probability distribution before and after the update, as shown in FIG. 14. FIG. 14 visually shows that in order to bring the defect rate to 0%, it is necessary not only to increase the average value of factor A by 0.15, but also to narrow the distribution width of factor A and to lower the average value of factor B by 0.2. Here, "OK" in FIG. 14 corresponds to "good" and "NG" corresponds to "bad."

以上のように、品質要因解析装置10は、従来技術に比べて高精度に品質要因を推定できる。 As described above, the quality factor analysis device 10 can estimate quality factors with higher accuracy than conventional techniques.

(実施例2)
本実施例では、厚板鋼材の強度・靱性の不良リスクを例として、製造プロセスの不良原因の解析と改善の方法が説明される。本実施例における因子は、図15の因果ダイアグラムで示されるように、圧延条件の因子として、鋼材の仕上幅[mm]、仕上厚[mm]、制御圧延中の圧延材の厚み[mm]を含み得る。また、本実施例における因子は、温度条件の因子として、加熱炉抽出温度、制御圧延温度、仕上温度、冷却温度を含み得る。また、本実施例における因子は、スラブの化学成分濃度の因子として、炭素、シリコン、マンガン、リン、サルファ-、銅、ニッケル、クロム、モリブテン、ニオブ、バナジウム、チタン、ボロン、アルミ、窒素、水素、酸素、カルシウムを含み得る。また、本実施例における因子は、製品の金属学因子として、炭化物析出量、鉄組織への元素固溶量、マルテンサイト分率を含み得る。また、本実施例における因子は、材質の因子として、引張強度、靱性を意味するシャルビー試験により算出される遷移温度、試験片ばらつき及び試験条件ばらつきを除いた引張強度と遷移温度を含み得る。また、本実施例における因子は、材質要求値の因子として、引張強度の上限値、引張強度の下限値、遷移温度の下限値を含み得る。また、本実施例における因子は、不良リスク因子として、引張強度の要求割れリスク、遷移温度の要求割れリスクを含み得る。ここで、図15において、黒塗りの四角は各因子の確率分布を簡略化して示したものである。
Example 2
In this embodiment, the risk of defects in the strength and toughness of a thick steel plate is taken as an example to explain a method of analyzing and improving the cause of defects in a manufacturing process. As shown in the causal diagram of FIG. 15, the factors in this embodiment may include the finishing width [mm] of the steel material, the finishing thickness [mm], and the thickness [mm] of the rolled material during controlled rolling as factors of rolling conditions. In addition, the factors in this embodiment may include the heating furnace extraction temperature, the controlled rolling temperature, the finishing temperature, and the cooling temperature as factors of temperature conditions. In addition, the factors in this embodiment may include carbon, silicon, manganese, phosphorus, sulfur, copper, nickel, chromium, molybdenum, niobium, vanadium, titanium, boron, aluminum, nitrogen, hydrogen, oxygen, and calcium as factors of chemical component concentrations of the slab. In addition, the factors in this embodiment may include the amount of carbide precipitation, the amount of element solid solution in the iron structure, and the martensite fraction as metallurgical factors of the product. In addition, the factors in this embodiment may include, as material factors, transition temperature calculated by a Charvy test, which means tensile strength and toughness, and tensile strength and transition temperature excluding test piece variation and test condition variation. In addition, the factors in this embodiment may include, as material requirement value factors, an upper limit value of tensile strength, a lower limit value of tensile strength, and a lower limit value of transition temperature. In addition, the factors in this embodiment may include, as defect risk factors, a risk of tensile strength failing the required value and a risk of transition temperature failing the required value. Here, in FIG. 15, the black squares show a simplified probability distribution of each factor.

本実施例では、(1)不良リスク診断、(2)不良要因の推定、(3)不良率ゼロの操業アクションの決定という3つ段階で、因果関係に基づいて品質要因を推定する。これらの3つの段階のそれぞれで、上記の実施形態における条件読取(ステップS202)、因果推論(ステップS203)及び結果表示(ステップS204)の処理が実行される。 In this embodiment, quality factors are estimated based on causal relationships in three stages: (1) defect risk diagnosis, (2) defect cause estimation, and (3) decision on operational actions to achieve zero defect rates. In each of these three stages, the processes of condition reading (step S202), causal inference (step S203), and result display (step S204) in the above embodiment are executed.

(不良リスク診断)
品質不良リスク診断は、圧延条件、温度条件、スラブの化学成分濃度の全ての条件が与えられたとき、その条件下で製造するとどれくらいの確率で不良が発生しそうか、というリスクを算出することである。上記の実施形態における初期化の処理の後に、条件読取の処理として、圧延条件、温度条件、スラブの化学成分濃度の各因子の条件が入力される。因子の値は1つの確定値ではなく、操業バラツキを考慮した分布として与えられる。上記の実施形態における因果推論の処理によって、金属学因子、材質の因子、材質要求値の因子、不良リスク因子が推論され、各因子の確率分布が更新される。図15は、このような推論によって、不合格の確率(NG)及び合格の確率(OK)が表示された様子を示す。
(Defect risk diagnosis)
The quality defect risk diagnosis is to calculate the risk of the occurrence of defects when all the conditions of rolling conditions, temperature conditions, and chemical component concentrations of the slab are given. After the initialization process in the above embodiment, the conditions of each factor of rolling conditions, temperature conditions, and chemical component concentrations of the slab are input as a condition reading process. The value of the factor is not a single definite value, but is given as a distribution considering the operational variation. By the causal inference process in the above embodiment, the metallurgical factor, material factor, material requirement value factor, and defect risk factor are inferred, and the probability distribution of each factor is updated. FIG. 15 shows the state in which the probability of rejection (NG) and the probability of acceptance (OK) are displayed by such inference.

図15は、現操業における圧延条件、温度条件、スラブの化学成分濃度のばらつき(確率分布)で、最終的な不良率が26.626%となり得ること(品質不良リスク診断の結果)を示している。 Figure 15 shows that the final defect rate could be 26.626% (result of quality defect risk diagnosis) due to the variations (probability distribution) in the rolling conditions, temperature conditions, and chemical component concentrations of the slabs in the current operation.

(不良要因の推定)
不良要因の推定では、不良要因を特定するために、条件読取の処理として、不良率が100%であると設定する。図16はこのような設定を行った場合を示す図である。上記の実施形態における因果推論の処理によって、入力された条件に基づき各因子の確率分布が更新される。結果表示の処理によって、各因子について、更新の前後の確率分布が表示される。
(Estimation of the cause of defects)
In estimating the cause of defects, in order to identify the cause of defects, the defect rate is set to 100% as a condition reading process. Fig. 16 is a diagram showing a case where such a setting is made. The process of causal inference in the above embodiment updates the probability distribution of each factor based on the inputted condition. The process of displaying the results displays the probability distribution of each factor before and after the update.

図17は、一例として、スラブの化学成分濃度の因子であるMn(マンガン濃度)の確率分布を示したものである。横軸は、-5.5~+5.5に正規化されたMn濃度である。縦軸は、各Mn濃度となる確率を示している。更新前を表すグラフのMn量の確率分布は0.5~1.5にピークがある。一方、更新後を表すグラフのMn量の確率分布は、-2.5~-1.5にピークがある。よって、不良の原因の一つはMn量が3単位ほど少ないためと推測できる。ここで、1単位は、横軸の正規化された濃度の1つの区間に対応する。 Figure 17 shows, as an example, the probability distribution of Mn (manganese concentration), which is a factor in the chemical component concentration of a slab. The horizontal axis is the Mn concentration normalized to -5.5 to +5.5. The vertical axis shows the probability of each Mn concentration. The probability distribution of the Mn amount in the graph representing the state before renewal has a peak at 0.5 to 1.5. On the other hand, the probability distribution of the Mn amount in the graph representing the state after renewal has a peak at -2.5 to -1.5. Therefore, it can be inferred that one of the causes of the defects is that the Mn amount is about 3 units too low. Here, 1 unit corresponds to one interval of the normalized concentration on the horizontal axis.

(不良率ゼロの操業アクション)
不良率ゼロの操業アクションでは、合格率が100%であるように、換言すると不良率がゼロであると設定する。図18はこのような設定を行った場合を示す図である。上記の実施形態における因果推論の処理によって、入力された条件に基づき各因子の確率分布が更新される。結果表示の処理によって、各因子について、更新の前後の確率分布が表示される。
(Zero defect rate operational action)
In the operation action of zero defect rate, the pass rate is set to 100%, in other words, the defect rate is set to zero. Fig. 18 is a diagram showing the case where such a setting is made. The process of causal inference in the above embodiment updates the probability distribution of each factor based on the input conditions. The process of displaying the results displays the probability distribution of each factor before and after the update.

図19は、不良率ゼロの場合のMnの確率分布を示す図である。横軸は、-5.5~+5.5に正規化されたMn濃度である。縦軸は、各Mn濃度となる確率を示している。更新前を表すグラフのMn量の確率分布は、平均-0.25、標準偏差2.15である。一方、更新後を表すグラフのMn量の確率分布は、平均0.63、標準偏差1.69である。このことから、不良率を減らすためには、Mn量の狙いを現行より0.89(=0.63-(-0.25))増やして、さらに0.46(=2.15-1.69)だけレンジを狭めるように制御する必要があることが分かる。 Figure 19 shows the probability distribution of Mn when the defect rate is zero. The horizontal axis is the Mn concentration normalized to -5.5 to +5.5. The vertical axis shows the probability of each Mn concentration. The probability distribution of the Mn amount in the graph representing the state before the update has an average of -0.25 and a standard deviation of 2.15. On the other hand, the probability distribution of the Mn amount in the graph representing the state after the update has an average of 0.63 and a standard deviation of 1.69. From this, it can be seen that in order to reduce the defect rate, it is necessary to increase the target Mn amount by 0.89 (= 0.63 - (-0.25)) from the current level, and further narrow the range by 0.46 (= 2.15 - 1.69).

ここで、本実施例において、スラブの化学成分濃度に基づいて製造コスト分布が更新されて、他の因子とともに表示される。製造コスト分布は、例えば図20のように表示される。横軸が正規化された費用で0~10の値で表現され、縦軸は確率を表す。更新前のグラフ(左軸の数値を使用)と更新後のグラフ(右軸の数値を使用)を比較すると、更新前と比べて更新後の値が大きい。そのため、Mnの狭レンジ化によって、コストが増大し得ることも推定することができる。 In this embodiment, the production cost distribution is updated based on the chemical component concentrations of the slab and displayed together with other factors. The production cost distribution is displayed, for example, as shown in Figure 20. The horizontal axis represents normalized costs expressed as values from 0 to 10, and the vertical axis represents probability. Comparing the graph before the update (using the values on the left axis) with the graph after the update (using the values on the right axis), the values after the update are larger than those before the update. Therefore, it can be estimated that narrowing the range of Mn may result in an increase in costs.

以上のように、品質要因解析装置10は高精度に品質要因を推定できる。また、品質要因解析装置10は、不良要因の推定も、不良率ゼロの操業アクションも、同様の操作によって容易に実行することが可能であり、さらにコストについての変動も解析することが可能である。 As described above, the quality factor analysis device 10 can estimate quality factors with high accuracy. Furthermore, the quality factor analysis device 10 can easily perform both estimation of defect factors and operational actions for zero defect rate with the same operations, and can also analyze fluctuations in costs.

本開示に係る実施形態について、諸図面及び実施例に基づき説明してきたが、当業者であれば本開示に基づき種々の変形又は修正を行うことが容易であることに注意されたい。従って、これらの変形又は修正は本開示の範囲に含まれることに留意されたい。例えば、各構成部又は各ステップなどに含まれる機能などは論理的に矛盾しないように再配置可能であり、複数の構成部又はステップなどを1つに組み合わせたり、或いは分割したりすることが可能である。本開示に係る実施形態は装置が備えるプロセッサにより実行されるプログラム又はプログラムを記録した記憶媒体としても実現し得るものである。本開示の範囲にはこれらも包含されるものと理解されたい。 Although the embodiments of the present disclosure have been described based on the drawings and examples, it should be noted that those skilled in the art would easily be able to make various modifications or amendments based on the present disclosure. Therefore, it should be noted that these modifications or amendments are included in the scope of the present disclosure. For example, the functions included in each component or step can be rearranged so as not to cause logical inconsistencies, and multiple components or steps can be combined into one or divided. The embodiments of the present disclosure can also be realized as a program executed by a processor included in the device or a storage medium having a program recorded thereon. It should be understood that these are also included in the scope of the present disclosure.

10 品質要因解析装置
20 表示装置
30 入力装置
40 使用者
101 演算処理部
102 品質要因解析プログラム
103 ROM
104 RAM
105 バス
106 初期化部
107 条件読取部
108 因果推論部
109 表示制御部
110 データベース
REFERENCE SIGNS LIST 10 Quality factor analysis device 20 Display device 30 Input device 40 User 101 Arithmetic processing unit 102 Quality factor analysis program 103 ROM
104 RAM
105 Bus 106 Initialization unit 107 Condition reading unit 108 Causal inference unit 109 Display control unit 110 Database

Claims (6)

製品の品質特性を表す品質因子、前記品質特性を得るための製造プロセスの操作量を表す操作因子、及び、前記操作量と前記品質特性との因果関係に影響し得る中間的な因子を表す中間因子を含む因子を用いて前記品質特性の要因を解析する品質要因解析装置であって、
過去に製造された前記製品における前記品質因子、前記操作因子及び前記中間因子に関する製造データと前記因子の間の因果関係を表現した因果接続表を蓄積するデータベースと、
前記因果接続表を読み込み、前記製造データから前記品質因子、前記操作因子及び前記中間因子の値の範囲を表す初期の確率分布を算出する初期化部と、
前記因子の一部に対して計測又は設定によって確定された条件を読み取る条件読取部と、
前記確定された条件、前記初期の確率分布及び前記因子の間の因果関係に基づいて条件付確率を算出し、前記条件付確率に基づいて前記因子の確率分布を更新する因果推論部と、
前記因子の少なくとも1つについて、前記初期の確率分布及び更新された確率分布を表示させる表示制御部と、を備える、品質要因解析装置。
1. A quality factor analysis device that analyzes factors of a quality characteristic using factors including quality factors that represent quality characteristics of a product, operation factors that represent operation amounts of a manufacturing process for obtaining the quality characteristics, and intermediate factors that represent intermediate factors that may affect a causal relationship between the operation amounts and the quality characteristics, comprising:
a database for storing manufacturing data relating to the quality factors, the operation factors, and the intermediate factors in the products manufactured in the past, and a causal connection table expressing causal relationships between the factors ;
an initialization unit that reads the causal connection table and calculates an initial probability distribution representing a range of values of the quality factors, the operation factors, and the intermediate factors from the manufacturing data;
a condition reading unit that reads conditions determined by measurement or setting for some of the factors;
a causal inference unit that calculates a conditional probability based on the determined condition, the initial probability distribution, and a causal relationship between the factors, and updates the probability distribution of the factors based on the conditional probability;
a display control unit that displays the initial probability distribution and an updated probability distribution for at least one of the factors.
前記表示制御部は前記因子の間の因果関係を視覚的に表示させる、請求項1に記載の品質要因解析装置。 The quality factor analysis device according to claim 1 , wherein the display control unit visually displays the causal relationships between the factors. 前記製品は鋼材であって、
前記因子の間の因果関係は、前記鋼材の機械的特性並びに内部及び表面の性質の少なくとも1つを表す前記品質因子と、前記鋼材の製造条件である前記操作因子と、前記中間因子との間の因果関係である、請求項1又は2に記載の品質要因解析装置。
The product is a steel material,
3. The quality factor analysis device according to claim 1 or 2, wherein the causal relationship between the factors is a causal relationship between the quality factor representing at least one of the mechanical properties and internal and surface properties of the steel material, the operational factor being a manufacturing condition of the steel material, and the intermediate factor.
前記品質因子は、前記鋼材の強度、延性、靱性、表面疵、介在物及び内質の少なくとも1つを表す、請求項3に記載の品質要因解析装置。 The quality factor analysis device according to claim 3, wherein the quality factor represents at least one of the strength, ductility, toughness, surface defects, inclusions, and internal quality of the steel material. 前記操作因子は、前記鋼材の製造条件のうち加熱条件、圧延条件、冷却条件及び素材の化学成分の少なくとも1つである、請求項3又は4に記載の品質要因解析装置。 The quality factor analysis device according to claim 3 or 4, wherein the operational factor is at least one of the heating conditions, rolling conditions, cooling conditions, and chemical components of the steel material among the manufacturing conditions of the steel material. 請求項3から5のいずれか一項に記載の品質要因解析装置が更新した前記操作因子の確率分布に基づいて、前記鋼材の製造条件を決定することを含む、製鉄所の操業方法。 A method for operating a steelworks, comprising determining the manufacturing conditions for the steel material based on the probability distribution of the operation factors updated by the quality factor analysis device described in any one of claims 3 to 5.
JP2022034779A 2022-03-07 2022-03-07 Quality factor analysis device and steelworks operation method Active JP7622673B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022034779A JP7622673B2 (en) 2022-03-07 2022-03-07 Quality factor analysis device and steelworks operation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022034779A JP7622673B2 (en) 2022-03-07 2022-03-07 Quality factor analysis device and steelworks operation method

Publications (2)

Publication Number Publication Date
JP2023130229A JP2023130229A (en) 2023-09-20
JP7622673B2 true JP7622673B2 (en) 2025-01-28

Family

ID=88025003

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022034779A Active JP7622673B2 (en) 2022-03-07 2022-03-07 Quality factor analysis device and steelworks operation method

Country Status (1)

Country Link
JP (1) JP7622673B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2026069849A1 (en) * 2024-09-25 2026-04-02 Jfeスチール株式会社 Information processing device, information processing method, and process execution method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009116842A (en) 2007-03-09 2009-05-28 Omron Corp Factor estimation support device, control method therefor, and factor estimation support program
JP2018116545A (en) 2017-01-19 2018-07-26 オムロン株式会社 Prediction model creation device, production facility monitoring system, and production facility monitoring method
JP2019057077A (en) 2017-09-20 2019-04-11 新日鐵住金株式会社 Prediction device, method, and program of product characteristics, and control system of manufacturing process
JP2022014878A (en) 2020-07-07 2022-01-20 Jfeスチール株式会社 Manufacturing specification determination support device, manufacturing specification determination support method, computer program, and computer readable recording medium

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009116842A (en) 2007-03-09 2009-05-28 Omron Corp Factor estimation support device, control method therefor, and factor estimation support program
JP2018116545A (en) 2017-01-19 2018-07-26 オムロン株式会社 Prediction model creation device, production facility monitoring system, and production facility monitoring method
JP2019057077A (en) 2017-09-20 2019-04-11 新日鐵住金株式会社 Prediction device, method, and program of product characteristics, and control system of manufacturing process
JP2022014878A (en) 2020-07-07 2022-01-20 Jfeスチール株式会社 Manufacturing specification determination support device, manufacturing specification determination support method, computer program, and computer readable recording medium

Also Published As

Publication number Publication date
JP2023130229A (en) 2023-09-20

Similar Documents

Publication Publication Date Title
KR102867502B1 (en) Quality abnormality analysis method, method of manufacturing metal material, and quality abnormality analysis device
JP7622673B2 (en) Quality factor analysis device and steelworks operation method
JP5195331B2 (en) Quality prediction apparatus, prediction method, program, and computer-readable recording medium in manufacturing process
KR20210110661A (en) Analysis system and analysis method
JP4365600B2 (en) Steel product quality design apparatus and steel product manufacturing method
JP2013084057A (en) Management method for product quality and management device for product quality
CN113834828A (en) Product quality analysis support system
JP4400253B2 (en) Quality influence factor analysis method, quality prediction method, quality control method, quality influence factor analysis device, quality prediction device, quality control device, quality influence factor analysis system, quality prediction system, quality control system, and computer program
JP4623020B2 (en) Optimal quality design support device
Wu et al. A quantitative causal analysis and optimization framework for inclusions of steel products
JP2009064054A (en) Product quality control method and control device
JP2013103227A (en) Method for controlling plate width of thick plate
Wang et al. A fast monitoring method for multiple operating batch processes with incomplete modeling data types
CN119446332A (en) A method and related equipment for online prediction and diagnosis of process data of steel
JP7540235B2 (en) METHOD FOR MANUFACTURING LEARNING MODEL USED IN QUALITY ASSESSMENT, LEARNING SYSTEM, QUALITY ASSESSMENT SYSTEM, AND QUALITY ASSESSMENT METHOD
JP2010207900A (en) Method of predicting deformation resistance in hot rolling
JP2002358120A (en) Batch plant operation management device
WO2022239157A1 (en) Learning control device for rolling process
JP7474220B2 (en) Deformation resistance prediction system, deformation resistance prediction method, and deformation resistance prediction device
JP7131964B2 (en) Estimation device, estimation system, estimation method and program
JP2013140127A (en) Method for predicting hardness in vicinity of welded part and method for maintaining vicinity of welded part
JP2009245141A (en) Output value prediction method, device thereof, and program for the method
CN117831659A (en) Method and device for online detection of quality of wide and thick plates, electronic equipment and storage medium
CN116240368B (en) Tempering embrittlement degree prediction method, tempering embrittlement degree prediction system, electronic equipment and storage medium
TW202609480A (en) Methods for generating quality prediction models for metallic materials, methods for predicting the quality of metallic materials, methods for estimating quality influencing factors of metallic materials, manufacturing methods for metallic materials, apparatus for generating quality prediction models for metallic materials, apparatus for predicting the quality of metallic materials, and apparatus for estimating quality influencing factors of metallic materials.

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231024

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240820

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241010

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: 20241217

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241230

R150 Certificate of patent or registration of utility model

Ref document number: 7622673

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150