JP2017178651A - セメントクリンカーの製造条件の予測方法 - Google Patents

Abstract

【課題】短時間でかつ高い精度でセメントクリンカーの製造条件を予測することができる方法を提供する。【解決手段】ニューラルネットワーク、セメントクリンカーの製造運転における監視データの実測値が、学習データから得られる特定の範囲内となる学習データの群、及び、モニターデータを用いて、σＬ＜σＭとなる大きい学習回数でニューラルネットワークの学習を行った後に、学習回数を減らしながらニューラルネットワークの学習をσＬ≧σＭとなるまで繰り返し、かつ、解析度判定値が第一の設定値を満たさない場合であって、解析度判定値が予め定めた第二の設定値を満たす場合、監視データの実測値が特定の限定条件によって形成される数値範囲内に属する場合であれば、ニューラルネットワークの入力層に当該の監視データの実測値を入力し、出力層から評価データの推測値を出力するセメントクリンカーの製造条件の予測方法。【選択図】図１

Description

本発明は、ニューラルネットワークを用いたセメントクリンカーの製造条件の予測方法に関する。

セメント原料及び焼成用燃料には硫黄化合物が含有されており、セメントクリンカーの焼成運転中のロータリーキルン内では、これら硫黄化合物の一部が揮発している。この揮発成分は、ロータリーキルンからの排気ガスの一部として、サスペンションプレヒータに導かれるが、この間、特に、キルン窯尻とサスペンションプレヒータの最下段サイクロンの間で、揮発成分とセメント原料が接触して凝縮し、設備の内壁面上にコーティング（被膜）を形成する。コーティングが形成されると、通風障害やサイクロンの詰り等を引き起こし、操業の中断に至る場合もある。

近年、資源循環型社会の実現に向け、セメントクリンカー製造におけるリサイクル資源の大量活用と、低品位の原料及び燃料の使用量の増大等が進められている。このような中、セメントクリンカー製造において、セメント原料及び焼成用燃料に含まれる硫黄化合物の量も増大していく傾向にある。このため、サスペンションプレヒータにおいて、コーティングが形成されることに起因するトラブルの発生が助長されてきている。

サスペンションプレヒータにおいて、コーティングが形成されることに起因するトラブルの発生防止のため、セメント工場では数時間に１回の頻度で、サスペンションプレヒータ内の原料（以後、「プレヒータ原料」と略称する。）を採取し、当該プレヒータ原料中のＳＯ_３量を分析して、必要に応じてセメントクリンカー製造プロセスにおいて用いられる焼成用燃料中のＳＯ_３量を調整したり、ロータリーキルンやサスペンションプレヒータ（仮焼炉）等の運転条件を調整するなどの作業が行われている。

しかしながら、プレヒータ原料の採取及び分析は煩雑であると共に、１回の作業に数時間を要するため、ロータリーキルンやサスペンションプレヒータ等の運転状況をきめ細かく把握してセメントクリンカー製造プロセスを安定させる方法については、改善の余地があった。

プレヒータ原料のＳＯ_３量（結果）を所定の範囲内にするために、セメントクリンカー製造プロセスにおいて用いられる焼成用燃料や、ロータリーキルンやサスペンションプレヒータの運転条件（原因）等をどのようにすべきか、という逆問題に対するアプローチ方法の一つとして、ニューラルネットワークを用いた予測方法がある。

セメントの製造方法に関するニューラルネットワークを用いた予測方法として、例えば、特許文献１には、学習データとモニターデータを用いて、σ_Ｌ＜σ_Ｍとなるような、十分に大きい学習回数でニューラルネットワークの学習を行った後に、学習回数を減らしながらニューラルネットワークの学習をσ_Ｌ≧σ_Ｍとなるまで繰り返し、学習後の解析度判定値が予め定めた設定値未満である場合に、ニューラルネットワークの入力層に、セメント製造における監視データの実測値を入力して、学習後のニューラルネットワークの出力層から、セメントの品質または製造条件の評価に関連する評価データの推測値を出力する、セメントの品質または製造条件の予測方法が記載されている。

特許第５３２３２９０号公報

特許文献１に記載された方法の予測の対象は、セメントの品質に直接的に関係するセメントの製造条件であり、セメントクリンカー製造プロセスの安定化を支援するための、セメントクリンカーの製造条件（特に、プレヒータ原料中のＳＯ_３量）を予測の対象としたものではない。
本発明の目的は、高い精度で、プレヒータ原料中のＳＯ_３量を予測することができる方法を提供することである。

本発明者らは、上記課題を解決するために鋭意検討した結果、ニューラルネットワークを用いたセメントクリンカーの製造条件の予測方法であって、学習データとモニターデータを用いて、σ_Ｌ＜σ_Ｍ（この式の意味は後で説明する。）となるような、十分に大きい学習回数でニューラルネットワークの学習を行った後に、学習回数を減らしながらニューラルネットワークの学習をσ_Ｌ≧σ_Ｍとなるまで繰り返し、次いで、解析度判定値を算出して、得られた値が、予め定めた設定値未満である場合には、以下のことを行えば、上記目的を達成しうることを見出し、本発明を完成した。
（ａ）セメントクリンカーの製造運転における監視データの実測値が、学習データの監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値である場合には、得られたニューラルネットワークの入力層にセメントクリンカーの製造運転における監視データの実測値を入力して、出力層からプレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力する。
（ｂ）セメントクリンカーの製造運転における監視データの実測値が、上記範囲外の数値であれば、特定の学習データを新たに追加した学習データ群を用いて、ニューラルネットワークの学習を行った後に、得られたニューラルネットワークの入力層に、セメントクリンカーの製造運転における監視データの実測値を入力して、出力層からプレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力する。

すなわち、本発明は、以下の［１］〜［５］を提供するものである。
［１］ 入力層及び出力層を有するニューラルネットワークを用いたセメントクリンカーの製造条件の予測方法であって、上記入力層は、セメントクリンカーの製造運転における監視データの実測値を入力するためのものであり、上記出力層は、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力するためのものであり、上記監視データが、セメントクリンカーの原料に関するデータ、セメントクリンカーの焼成条件に関するデータ、及びセメントクリンカーに関するデータの中から選ばれる一種以上であり、
（Ａ）ニューラルネットワークの学習に使用する、監視データの実測値と評価データの実測値の組み合わせである学習データの複数からなる学習データ群の初期設定を行う工程と、
（Ｂ）学習回数の初期設定を行う工程と、
（Ｃ）設定された学習データ群を用いて、ニューラルネットワークの学習を、前工程で設定された学習回数行う工程と、
（Ｄ）学習データの監視データの実測値を、直近の工程（Ｃ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）、及び、ニューラルネットワークの学習結果の信頼性を確認するための監視データの実測値と評価データの実測値の組み合わせであるモニターデータの中の監視データの実測値を、直近の工程（Ｃ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの中の評価データの実測値との平均２乗誤差（σ_Ｍ）を算出し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍである場合、工程（Ｅ）を実施し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ＜σ_Ｍである場合、工程（Ｆ）を実施する工程と、
（Ｅ）直近の工程（Ｂ）で設定された学習回数および再設定された直近のニューラルネットワークの学習回数のいずれの学習回数よりも大きい学習回数を、工程（Ｂ）において新たな学習回数として再設定し、再度工程（Ｃ）〜（Ｄ）を実施する工程と、
（Ｆ）直近のニューラルネットワークの学習で実施された学習回数を減らした学習回数を、新たな学習回数として再設定する工程と、
（Ｇ）設定された学習データ群を用いて、ニューラルネットワークの学習を直近の工程（Ｆ）で設定された学習回数行う工程と、
（Ｈ）学習データの監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）、及び、モニターデータの中の監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの中の評価データの実測値との平均２乗誤差（σ_Ｍ）を算出し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍである場合、工程（Ｊ）を実施し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ＜σ_Ｍである場合、工程（Ｉ）を実施する工程と、
（Ｉ）直近の工程（Ｇ）におけるニューラルネットワークの学習回数が予め設定された回数を超えている場合、再度工程（Ｆ）〜（Ｈ）を実施し、直近の工程（Ｇ）におけるニューラルネットワークの学習回数が予め定めた数値以下の場合、工程（Ｍ）を実施する工程と、
（Ｊ）下記式（１）を用いて解析度判定値を算出し、該解析度判定値が予め定めた第一の設定値未満である場合、ニューラルネットワークの学習を終了し工程（Ｋ）を実施する工程と、上記解析度判定値が予め定めた第一の設定値以上である場合、工程（Ｍ）を実施する工程と、
（Ｋ）セメントクリンカーの製造運転における監視データの実測値が、設定された学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値である場合には、工程（Ｒ）を実施し、範囲外の数値である場合には、工程（Ｌ）を実施する工程と、
（Ｌ）監視データの実測値が、設定された学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値である学習データであって、上記学習データ群に含まれない１個以上の学習データを、上記学習データ群に追加して新たな学習データ群に設定した後、後、工程（Ｃ）以降を実施する工程と
（Ｍ）工程（Ｂ）を実施した回数の大きさについての判定を行い、該回数が予め設定した回数以下である場合は、工程（Ｎ）を実施し、該回数が予め設定した回数を超える場合には工程（Ｏ）を実施する工程と、
（Ｎ）学習条件の初期化を行って、再度工程（Ｂ）〜（Ｊ）を実施する工程と、
（Ｏ）工程（Ｊ）において算出した全ての解析度判定値のうち、最も小さい解析度判定値が予め定めた第二の設定値未満である場合、最も小さい解析度判定値を得ることができた工程（Ｊ）におけるニューラルネットワークを、学習済みのニューラルネットワークとした後、工程（Ｐ）を実施し、最も小さい解析度判定値が予め定めた第二の設定値以上である場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
（Ｐ）工程（Ｊ）において算出した全ての解析度判定値のうち、最も小さい解析度判定値を得ることができた工程（Ｊ）において、学習データとして使用した監視データの実測値と評価データの実測値の組み合わせについて無相関検定を行い、５％の有意水準で有意であると判断された監視データの種類が２種以上である場合、５％の有意水準で有意であると判断された監視データの全種類を座標軸とする座標空間に学習データとして使用した監視データの実測値をプロットし、座標空間において、プロットされた監視データ同士を結ぶことで形成される監視データの全てを包含する領域であって、該領域が最大となるように監視データ同士を結ぶことで形成される領域を、予測可能監視データ領域として設定した後、工程（Ｑ）を実施し、５％の有意水準で有意であると判断された監視データが０または１種類である場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
（Ｑ）セメントクリンカー製造運転における監視データの実測値が、工程（Ｐ）で設定した予測可能監視データ領域に含まれる場合、セメントクリンカーの製造条件の予測を高い精度で行うことができると判断し、工程（Ｋ）を実施し、セメントクリンカー製造運転における監視データの実測値が、上記予測可能監視データ領域に含まれない場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
（Ｒ）学習済みのニューラルネットワークの入力層に、セメントクリンカー製造運転における監視データの実測値を入力して、上記ニューラルネットワークの出力層から、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力して、セメントクリンカーの製造条件を予測する工程と、
を含むことを特徴とするセメントクリンカーの製造条件の予測方法。

（上記式（１）中、学習データの平均２乗誤差（σ_Ｌ）とは、学習データの監視データの実測値を学習後のニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）である。評価データの推測値の平均値とは、学習データの監視データの実測値を学習後のニューラルネットワークの入力層に入力して得られた評価データの推測値の平均値である。）

［２］ 上記工程（Ｌ）で追加する学習データの数が、設定された学習データ群における学習データの数の８％以上である前記［１］に記載のセメントクリンカーの製造条件の予測方法。
［３］ 上記解析度判定値の予め定めた第一の設定値が１０％以下であり、上記解析度判定値の予め定めた第二の設定値が上記第一の設定値よりも大きくかつ３０％以下である、前記［１］または［２］に記載のセメントクリンカーの製造条件の予測方法。
［４］ 上記ニューラルネットワークが、上記入力層と上記出力層の間に中間層を有する階層型のニューラルネットワークである、前記［１］〜［３］のいずれかに記載のセメントクリンカーの製造条件の予測方法。
［５］ 上記監視データの値を人為的に変動させて得られた上記評価データの推測値に基づいて、プレヒータ原料中のＳＯ_３量を最適化する、前記［１］〜［４］のいずれかに記載のセメントクリンカーの製造条件の予測方法。

本発明のセメントクリンカーの製造条件の予測方法を用いれば、セメントクリンカーの製造運転で得られる様々なデータに基づいて、プレヒータ原料中のＳＯ_３量を短時間で、且つ高い精度で予測することができる。
また、ニューラルネットワークの入力層に入力する、セメントクリンカーの製造運転における監視データの実測値が、ニューラルネットワークの学習に用いた監視データの実測値のデータ範囲の上限または下限付近である場合においても、短時間でかつ高い精度でセメントクリンカーの製造条件を予測することができる。
また、得られた推測値を基に、リアルタイムでセメントクリンカーの製造運転を最適化することが可能であり、セメントクリンカーの製造プロセスの安定化の向上を図ることができる。
さらに、ニューラルネットワークの学習を継続することによって、高い予測の精度を維持することができる。

本発明の予測方法の一例を示すフロー図である。 実施例１において設定した予測可能監視データ領域を示す図である。

以下、本発明について詳細に説明する。
本発明の予測方法は、セメントクリンカーの製造運転における監視データの実測値を入力するための入力層と、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力するための出力層を有するニューラルネットワークを用いて、セメントクリンカーの製造条件を予測する方法である。
本発明のニューラルネットワークは、入力層と出力層の間に中間層を有する階層型のニューラルネットワークであってもよい。

上記監視データは、セメントクリンカーの原料に関するデータ、セメントクリンカーの焼成条件に関するデータ、及びセメントクリンカーに関するデータの中から選ばれる一種以上のデータである。
上記評価データは、プレヒータ原料中のＳＯ_３量に関するデータである。

監視データの一つである「セメントクリンカーの原料に関するデータ」は、セメントクリンカーの調合原料の化学組成、水硬率、ふるい試験残分量、ブレーン比表面積（粉末度）、強熱減量や、キルンへの投入時から所定の時間前の時点（例えば、５時間前の１つの時点や、３時間前、４時間前、５時間前、及び６時間前の４つの時点のような複数の時点）のセメントクリンカーの原料（搬送中に向流する空気流によって微粒分等が抜き取られたセメントクリンカーの調合原料。以後、「セメントクリンカーの窯入原料」と称す。）の化学組成、水硬率、供給量や、廃棄物のような特殊な原料からなるセメントクリンカーの副原料の供給量や、調合原料のブレンディングサイロの貯留量（残量）や、調合原料のストレージサイロの貯留量（残量）や、原料ミルと調合原料のブレンディングサイロの間に位置するサイクロンの電流値（サイクロンの回転数を表し、サイクロンを通過する原料の速度と相関関係があるもの）や、セメントクリンカーの窯入原料と副原料を混合してなる原料の化学組成、水硬率、ブレーン比表面積、ふるい試験残分量、脱炭酸率、水分量、等が挙げられる。これらのデータは、１種を単独でまたは２種以上を組み合わせて用いられる。
ここで、セメントクリンカーの原料（調合原料または窯入原料）の化学組成とは、セメントクリンカーの原料中のＳｉＯ_２、Ａｌ_２Ｏ_３、Ｆｅ_２Ｏ_３、ＣａＯ、ＭｇＯ、ＳＯ_３、Ｎａ_２Ｏ、Ｋ_２Ｏ、Ｎａ_２Ｏｅｑ（全アルカリ）、ＴｉＯ_２、Ｐ_２Ｏ_５、ＭｎＯ、Ｃｌ、Ｃｒ、Ｚｎ、Ｐｂ、Ｃｕ、Ｎｉ、Ｖ、Ａｓ、Ｚｒ、Ｍｏ、Ｓｒ、Ｂａ、Ｆ等の含有率である。

監視データの一つである「セメントクリンカーの焼成条件に関するデータ」は、セメントクリンカーの原料の挿入量、キルン回転数、キルン平均トルク、落口温度、焼成帯温度、セメントクリンカー温度、ボトムサイクロンの出口温度、クリンカークーラー温度、プレヒーターのガスの流量（プレヒーターの温度と相関関係があるもの）、仮焼炉出口ガスのＯ_２濃度やＣＯ濃度等の仮焼炉出口ガス成分濃度、キルンＥＰ出口ガスのＳＯ_２濃度やＮＯ_Ｘ濃度等のキルンＥＰ出口ガス成分濃度、キルン窯尻ガスのＯ_２濃度やＮＯ_Ｘ濃度等のキルン窯尻ガス成分濃度、塩素バイパス抽気ガスのＳＯ_Ｘ濃度等の塩素バイパス抽気ガス成分濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度等のプレヒータトップサイクロン出口ガス成分濃度等が挙げられる。これらのデータは、１種を単独でまたは２種以上を組み合わせて用いられる。

監視データの一つである「セメントクリンカーに関するデータ」は、セメントクリンカーの鉱物組成、各鉱物の結晶学的性質（格子定数や結晶子径など）、２種以上の鉱物組成の比、化学組成、容重、ＳＯ_３含有量、Ｃｌ含有量等が挙げられる。これらのデータは、１種を単独でまたは２種以上を組み合わせて用いられる。

ここで、セメントクリンカーの鉱物組成とは、３ＣａＯ・ＳｉＯ_２（Ｃ_３Ｓ）、２ＣａＯ・ＳｉＯ_２（Ｃ_２Ｓ）、３ＣａＯ・Ａｌ_２Ｏ_３（Ｃ_３Ａ）、４ＣａＯ・Ａｌ_２Ｏ_３・Ｆｅ_２Ｏ_３（Ｃ_４ＡＦ）、ｆ．ＣａＯ、ｆ．ＭｇＯ等の含有率である。また「２種以上の鉱物組成の比」としては、例えば、Ｃ_３Ｓ／Ｃ_２Ｓの比が挙げられる。なお、セメントクリンカーの鉱物組成は、例えばＸＲＤ−リートベルト法によって得ることができる。
セメントクリンカーの化学組成とは、セメントクリンカー中のＳｉＯ_２、Ａｌ_２Ｏ_３、Ｆｅ_２Ｏ_３、ＣａＯ、ＭｇＯ、ＳＯ_３、Ｎａ_２Ｏ、Ｋ_２Ｏ、Ｎａ_２Ｏｅｑ（全アルカリ）、ＴｉＯ_２、Ｐ_２Ｏ_５、ＭｎＯ、Ｃｌ、Ｃｒ、Ｚｎ、Ｐｂ、Ｃｕ、Ｎｉ、Ｖ、Ａｓ、Ｚｒ、Ｍｏ、Ｓｒ、Ｂａ、Ｆ等の含有率である。
監視データは、セメントクリンカーの原料に関するデータ、セメントクリンカーの焼成条件に関するデータ、及びセメントクリンカーに関するデータの中から選ばれるいずれか一種のデータのみを用いてもよいが、これら３種のデータのうちの２種以上（複数）のデータを用いることが、評価データの予測の精度を高める観点から、好ましい。

本発明のセメントクリンカーの製造条件の予測方法において、対象となるセメントクリンカーは、特に限定されず、例えば、普通ポルトランドセメント、早強ポルトランドセメント、中庸熱ポルトランドセメント、低熱ポルトランドセメント等の各種ポルトランドセメント用や、エコセメント用のセメントクリンカーが挙げられる。
セメントクリンカーの製造工程は、原料工程、焼成工程の２工程に大別される。原料工程は、石灰石、粘土、珪石、酸化鉄原料などのセメントクリンカー原料を適当な割合で調合して、原料ミルで微粉砕し、セメントクリンカーの調合原料を得る工程である。焼成工程は、セメントクリンカーの調合原料を、サスペンションプレヒーターを経由してロータリーキルンに供給し、充分に焼成した後、冷却して、セメントクリンカーを得る工程である。

本発明では、セメントクリンカー製造運転における監視データと、プレヒータ原料中のＳＯ_３量に関する評価データの関係を、ニューラルネットワークによって予め学習し、その学習結果を用いて、上記監視データのみに基づいて、上記評価データを予測する。
以下、本発明の予測方法について、図１を参照しながら詳しく説明する。
［工程（Ａ）］
工程（Ａ）において、ニューラルネットワークの学習に使用する、監視データの実測値と評価データの実測値の組み合わせである学習データの複数からなる学習データ群の初期設定を行う。
学習データは、学習用のサンプルを用意し、該サンプルの監視データの実測値、及び評価データの実測値を測定することで得ることができる。学習データは、作業性の観点から、予め十分な数を用意することが好ましいが、後述する工程（Ｌ）において、学習データ群に新たな学習データを追加する場合等、必要に応じて、新たな学習用のサンプルを用意し、該サンプルから新たな学習データを得てもよい。

学習データ群は、複数の学習用のサンプルから得られた全ての学習データからなるものでもよく、全ての学習データから選択した複数（ただし、全てのデータより１つ以上のデータを削除したもの）の学習データからなるものであってもよい。
全ての学習データから複数の学習データを選択して、学習データ群とする場合、学習データの選択は、全ての学習データの中から、監視データの実測値が、本発明の予測方法の対象となるセメントクリンカーの製造運転における監視データの実測値と大きく異ならないものを選択することが好ましい。
具体的には、全ての学習データの中から、セメントクリンカー原料の構成が、本発明の予測方法の対象となるセメントクリンカーの製造運転におけるセメントクリンカー原料の構成と大きく異ならない学習データや、製造設備の変更や更新、さらには、製造工程の変更等が生じていない学習データを選択する。セメントクリンカー原料の構成が大きく異なる学習データや、製造設備の変更等が生じている学習データを選択した場合、該学習データがノイズとなって、予測精度の向上が困難となる場合がある。

学習データ群を構成する学習データの数は、特に限定されるものではないが、予測精度を高くする観点から、好ましくは５以上、より好ましくは５０以上、さらに好ましくは１００以上、特に好ましくは１５０以上である。上記数の上限は、特に限定されないが、例えば、１，０００である。
工程（Ａ）終了後、工程（Ｂ）を実施する。

［工程（Ｂ）］
工程（Ｂ）において、学習回数の初期設定を実施する。設定される学習回数は、特に限定されるものではないが、好ましくは、ニューラルネットワークの過学習（オーバーラーニング）が発生する程度に、十分に大きな回数である。具体的には、通常５千〜１００万回、好ましくは１万〜１０万回である。
工程（Ｂ）では、ニューラルネットワークの過学習が発生する学習回数、具体的にはσ_Ｌ＜σ_Ｍ（詳しくは後述する）となるような学習回数を設定することが好ましいが、後の工程において、学習回数の増減が行われるため、工程（Ｂ）において最初に設定される学習回数は、ニューラルネットワークの学習に通常行われる学習回数を用いても問題ない。
工程（Ｂ）終了後、工程（Ｃ）を実施する。

［工程（Ｃ）］
工程（Ｃ）では、設定された学習データ群を用いて、ニューラルネットワークの学習を、前工程で設定された学習回数行う。
ここで、「設定された学習データ群」とは、工程（Ａ）または工程（Ｌ）において設定された学習データ群であって、直近に設定されたものをいう。
また、「前工程で設定された学習回数」とは、工程（Ｂ）において設定される学習回数、または、工程（Ｅ）において再設定された新たな学習回数であって、直近の工程（工程（Ｂ）または工程（Ｅ））で設定された学習回数である。
具体的には、設定された学習データ群から、学習データの監視データの実測値をニューラルネットワークの入力層に入力して、出力層から出力された評価データの推測値と、該評価データの推測値に対応する学習データの評価データの実測値を比較評価してニューラルネットワークの修正することを、設定された学習回数行うことで、ニューラルネットワークの学習が行われる。
なお、学習回数を変更して、ニューラルネットワークの再学習を行う際には、前回の学習の結果得られたニューラルネットワークは初期化され、再度学習が行われる。
工程（Ｃ）終了後、工程（Ｄ）を実施する。

［工程（Ｄ）］
工程（Ｄ）では、σ_Ｌとσ_Ｍが算出される。σ_Ｌとσ_Ｍの大小関係から、学習がニューラルネットワークの過学習が発生する程度に十分に大きな回数行われたか否かを判断することができる。
具体的には、学習データの監視データの実測値を、直近の工程（Ｃ）において学習が行われたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）を算出する。次いで、モニターデータの監視データの実測値を、直近の工程（Ｃ）において学習が行われたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの評価データの実測値との平均２乗誤差（σ_Ｍ）を算出する。その後、算出されたσ_Ｌとσ_Ｍの数値を比較することで、ニューラルネットワークの学習が十分に大きな回数で行われたか判断することができる。
ここで、モニターデータとは、学習データを得るために用いられたサンプルとは別のサンプルから得られた、監視データの実測値及び評価データの実測値の組み合わせであり、ニューラルネットワークの信頼性を確認するためのデータである。
モニターデータ（監視データの実測値及び評価データの実測値の組み合わせ）のサンプルの数は、作業性の観点から、学習データのサンプル数の好ましくは５〜５０％、より好ましくは１０〜３０％である。

工程（Ｄ）で算出されたσ_Ｌとσ_Ｍの関係が、σ_Ｌ≧σ_Ｍである場合（図１の過学習判定における「Ｎｏ」）、直近に行った工程（Ｃ）の学習回数は、十分に大きな回数ではないと判断することができる。この場合、工程（Ｅ）を実施する。工程（Ｄ）で算出されたσ_Ｌとσ_Ｍの関係が、σ_Ｌ＜σ_Ｍである場合（図１の過学習判定における「Ｙｅｓ」）、直近に行った工程（Ｃ）の学習回数は、十分に大きな回数であったと判断することができる。この場合、工程（Ｆ）を実施する。

［工程（Ｅ）］
工程（Ｅ）では、直近の工程（Ｂ）で設定された学習回数および再設定された直近のニューラルネットワークの学習回数のいずれの学習回数よりも大きい学習回数を新たな学習回数として再設定する（例えば、直近の工程（Ｃ）で実施された学習回数に２．０を乗じた数を新たな学習回数として設定する。）。新たな学習回数を再設定した後、再度工程（Ｃ）〜（Ｄ）を実施する。

［工程（Ｆ）］
工程（Ｆ）では、直近のニューラルネットワークの学習で実施された学習回数を減らした学習回数を、新たな学習回数として再設定する（例えば、直近のニューラルネットワークの学習で実施された学習回数に０．９５を乗じた数を新たな学習回数として設定する。）。
なお、直近のニューラルネットワークの学習とは、より近い過去に実施された学習を指す。具体的には、工程（Ｃ）もしくは後述の工程（Ｇ）のうち、より近い過去に実施された学習を指す。工程（Ｆ）終了後、工程（Ｇ）を実施する。

［工程（Ｇ）］
工程（Ｇ）では、設定された学習データ群（工程（Ａ）または工程（Ｌ）において設定された学習データ群であって、直近に設定されたもの）を用いて、ニューラルネットワークの学習を直近の工程（Ｆ）で設定された学習回数行う。工程（Ｇ）で実施する内容は、ニューラルネットワークの学習を工程（Ｆ）において新たに設定された学習回数行う以外は、工程（Ｃ）と同じである。
工程（Ｇ）終了後、工程（Ｈ）を実施する。

［工程（Ｈ）］
工程（Ｈ）では、直近の工程（Ｇ）の学習において得られたニューラルネットワークを用いて終了判定を行う。具体的には学習データの監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）、及び、モニターデータの中の監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの中の評価データの実測値との平均２乗誤差（σ_Ｍ）を算出し、算出されたσ_Ｌとσ_Ｍの関係が、σ_Ｌ≧σ_Ｍである場合（図１の終了判定における「Ｙｅｓ」）、直近に行った工程（Ｇ）の学習回数は、もはや十分に大きな回数ではないと判断することができる。この場合、後述する工程（Ｊ）を実施する。算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ＜σ_Ｍである場合（図１の終了判定における「Ｎｏ」）、直近に行った工程（Ｇ）の学習回数は、いまだ十分に大きな回数であったと判断することができる。この場合、後述する工程（Ｉ）を実施する。

［工程（Ｉ）］
工程（Ｉ）では、直近の工程（Ｇ）におけるニューラルネットワークの学習回数が予め定めた数値を超えていたかどうかの判定を行う。工程（Ｉ）は、工程（Ｆ）〜（Ｈ）を無限に繰り返すことを回避するために行われる。工程（Ｉ）において直近に行った工程（Ｇ）におけるニューラルネットワークの学習回数が予め定めた数値を超えていた場合（図１における「Ｙｅｓ」）は、再度工程（Ｆ）〜（Ｈ）を実施する。工程（Ｉ）において直近に行った工程（Ｇ）の学習回数が予め定めた数値以下場合（図１における「Ｎｏ」）は、後述の工程（Ｍ）を実施する。
なお、上記予め定めた数値とは、特に限定されず、例えば、工程（Ｆ）で設定された学習回数の１００分の１の数値以下、もしくは、１以下または０以下等が挙げられる。

［工程（Ｊ）］
工程（Ｊ）では解析度判定値が予め定めた第一の設定値未満であるか否かによって、解析度の判定を行うことができる。解析度判定値は下記式（１）を用いて算出される。

上記式（１）中、学習データの平均２乗誤差（σ_Ｌ）とは、直近の工程（Ｈ）で算出された平均２乗誤差（σ_Ｌ）と同じである。評価データの推測値の平均値とは、学習データの監視データの実測値を、直近の工程（Ｇ）にて得られたニューラルネットワークの入力層に入力して得られた評価データの推測値の平均値である。
解析度の判定を行うことで学習を行ったニューラルネットワークを用いて、セメントクリンカーの製造条件（プレヒータ原料中のＳＯ_３量）の予測を高い精度で行うことができるか否かを判断することができる。解析度判定値が予め定めた第一の設定値未満（図１の第一の解析度判定における「Ｙｅｓ」）であれば、解析は十分であると判断され、ニューラルネットワークの学習は終了し、工程（Ｋ）を実施する。
解析度判定値が予め定めた第一の設定値以上（図１の第一の解析度判定における「Ｎｏ」）であれば、学習データを用いて学習を行ったニューラルネットワークをそのまま用いて、セメントの品質等の予測を高い精度で行うことはできないと判断され、工程（Ｍ）を実施する。

予め定めた第一の設定値は、特に限定されないが、より高い精度で予測を行う観点から、好ましくは１０％以下、より好ましくは８％以下、特に好ましくは７％以下の値である。
なお、工程（Ｂ）〜（Ｊ）は、工程（Ｊ）において解析度判定値が予め定めた第一の設定値未満となるか、あるいは、工程（Ｍ）において該回数が予め設定した回数を超えるまで繰り返される。工程（Ｊ）を実施するたびに得られる、解析度判定値及び学習済みのニューラルネットワークは、工程（Ｏ）において使用するため、データとして保存する。
また、予測精度の向上の観点から、工程（Ｌ）を実施した場合、工程（Ｌ）を実施する前に、工程（Ｊ）を実施するたびに得られた、解析度判定値及び学習済みのニューラルネットワークのデータは、工程（Ｏ）において使用せずに、破棄する。

［工程（Ｋ）］
工程（Ｋ）では、セメントクリンカーの製造運転における監視データの実測値（工程（Ｒ）において、学習が終了したニューラルネットワークの入力層に入力するもの）が、設定された学習データ群（工程（Ａ）または工程（Ｌ）において設定された学習データ群であって、直近に設定されたもの）における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値である場合（図１の学習データの判定における「Ｙｅｓ」）には、工程（Ｒ）を実施し、範囲外の数値である場合（図１の学習データの判定における「Ｎｏ」）には、工程（Ｌ）を実施する。
セメントクリンカーの製造運転における監視データの実測値が複数の種類ある場合、そのうち１種類でも範囲外の数値である場合には、工程（Ｌ）を実施する。
本工程において、セメントクリンカーの製造運転における監視データの実測値が、上記範囲内の数値であるか否かを判断することによって、学習データ群を構成する学習データの数を増加して、ニューラルネットワークの学習を再度行う必要があるか否かを判断することができる。これにより、セメントクリンカーの製造条件の評価に関連する評価データの推測値の予測精度をより向上させることができる。

なお、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値とは、該数値が、監視データの実測値の平均値−平均２乗誤差（σ_Ｇ）以上であり、かつ、監視データの実測値の平均値＋平均２乗誤差（σ_Ｇ）以下であることを意味する。また、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値とは、該数値が、監視データの実測値の平均値＋平均２乗誤差（σ_Ｇ）を超える、または、監視データの実測値の平均値−平均２乗誤差（σ_Ｇ）未満であることを意味する。

［工程（Ｌ）］
工程（Ｌ）では、監視データの実測値が、設定された学習データ群（工程（Ａ）または工程（Ｌ）において設定された学習データ群であって、直近に設定されたもの）における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値である学習データであって、上記学習データ群に含まれない１個以上の学習データを、上記学習データ群に追加して新たな学習データ群に設定した後、工程（Ｃ）以降を実施する（工程（Ｃ）〜（Ｔ）から適宜選択される工程を実施する）。
学習データの監視データの実測値が複数の種類ある場合、複数の種類の監視データの実測値の各々が、学習データの監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値である学習データを使用する。
上記範囲外の数値である、監視データの実測値が複数の種類ある場合、該複数の種類の監視データの実測値の各々について、学習データの監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値である学習データであって、上記学習データ群に含まれない１個以上の学習データを、上記学習データ群に追加して新たな学習データ群に設定する。
なお、追加する学習データの監視データの実測値の種類のうち、学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値であった監視データの種類の実測値については、特に考慮しなくてもよい。
上記学習データを追加することで、ニューラルネットワークの入力層に入力する、セメントクリンカーの製造運転における監視データの実測値が、ニューラルネットワークの学習に用いた学習データの監視データの実測値のデータ範囲の上限または下限付近である場合であっても、予測精度を向上することができる。

追加する学習データの数は、設定された学習データ群における学習データの数の、好ましくは８％以上、より好ましくは１０％以上、特に好ましくは１５％以上の数である。該数が８％以上であれば予測精度をより向上することができる。
なお、工程（Ｌ）は、工程（Ｋ）において、セメント製造における監視データの実測値が、設定された学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値となるまで、繰り返し実施される。

［工程（Ｍ）］
工程（Ｍ）では、学習回数を設定する工程（Ｂ）を実施した回数が予め設定した数値以下であるかどうかの判定を実施する。判定を実施することによって、工程（Ｂ）から工程（Ｊ）を無限に繰り返すことを回避することができる。
工程（Ｍ）において、工程（Ｂ）を実施した回数が予め設定した回数以下（図１の回数判定における「Ｙｅｓ」）である場合、学習条件の初期化を行って、再度工程（Ｂ）〜（Ｊ）を行い、該回数が予め設定した回数を超える場合（図１の回数判定における「Ｎｏ」）、工程（Ｏ）を実施する。
予め設定した回数は、特に限定されないが、通常、５回以上である。予め設定した回数の上限は、工程（Ｂ）から工程（Ｊ）を多大に繰り返すことを防ぐ観点から、好ましくは１００回以下である。

［工程（Ｎ）］
工程（Ｎ）では、学習条件の初期化を行って、再度工程（Ｂ）〜（Ｊ）を実施する。
学習条件の初期化の方法としては、例えば、ニューラルネットワークを構成するユニットの閾値やユニットを結合している重みをランダムで変更した上で、学習データを再入力する方法、学習データを得るためのサンプルの数を増やす、使用する監視データの種類を変更する、又は不適切な学習データを除外する等を行った上で、新たな学習データを入力する方法等が挙げられる。

［工程（Ｏ）］
工程（Ｏ）では、工程（Ｊ）において算出した全ての解析度判定値（ただし、工程（Ｊ）において、工程（Ｌ）を実施する前に、工程（Ｊ）を実施するたびに得られる、解析度判定値及び学習済みのニューラルネットワークのデータを破棄した場合、工程（Ｌ）を実施した後に、工程（Ｊ）において算出した全ての解析度判定値）のうち、最も小さい解析度判定値が予め定めた第二の設定値未満であるか否かによって、次の予測の実施の可否の判定を行うことができる。
工程（Ｏ）の判定を追加することで、工程（Ｊ）において、セメントクリンカーの製造条件（プレヒータ原料中のＳＯ_３量）の予測を高い精度で行うことはできないと判断された学習済みのニューラルネットワークであっても、次工程（Ｐ）〜（Ｑ）を実施することによって、セメントクリンカーの製造条件の予測を高い精度で行うことができるか否かを判断することができる。最も小さい解析度判定値が予め定めた第二の設定値未満（図１の第二の解析度判定における「Ｙｅｓ」）である場合、最も小さい解析度判定値を得ることができた工程（Ｊ）におけるニューラルネットワークを、学習済みのニューラルネットワークとして得た後、工程（Ｐ）を実施する。

最も小さい解析度判定値が予め定めた第二の設定値以上（図１の第二の解析度判定における「Ｎｏ」）であれば、学習データを用いて学習を行ったニューラルネットワークを用いて、セメントクリンカーの製造条件の予測を高い精度で行うことはできないと判断して予測を終了する。
予め定めた第二の設定値は、上記第一の設定値よりも大きいものである。また、上限は、より高い精度で予測を行う観点から、好ましくは３０％以下、より好ましくは２５％である。

［工程（Ｐ）］
工程（Ｐ）では、次工程（Ｑ）で用いられる予測可能監視データ領域を設定する。
最初に、工程（Ｊ）において算出した全ての解析度判定値（ただし、工程（Ｊ）において、工程（Ｌ）を実施する前に、工程（Ｊ）を実施するたびに得られる、解析度判定値及び学習済みのニューラルネットワークのデータを破棄した場合、工程（Ｌ）を実施した後に、工程（Ｊ）において算出した全ての解析度判定値）のうち、最も小さい解析度判定値を得ることができた工程（Ｊ）において、学習データとして使用した監視データの実測値と評価データの実測値の組み合わせについて無相関検定を実施する。無相関検定において５％の有意水準で有意であると判断された監視データの種類が２種以上である場合（図１の無相関検定における「Ｙｅｓ」）、５％の有意水準で有意であると判断された監視データの全種類を座標軸とする座標空間を作成する。
例えば、５％の有意水準で有意であると判断された監視データが、塩素バイパスの抽気ガスのＳＯ_ｘ濃度とキルンＥＰ出口ガスのＳＯ_２濃度の二種類である場合、塩素バイパスの抽気ガスのＳＯ_ｘ濃度をｘ軸とし、キルンＥＰ出口ガスのＳＯ_２濃度をｙ軸とする座標空間を作成する。

次いで、学習データとして使用した監視データの実測値のうち、５％の有意水準で有意であると判断された種類の監視データの実測値を全て、座標空間にプロットし、座標空間においてプロットされた監視データ同士を結ぶことで予測可能監視データ領域を設定する。該予測可能監視データ領域は、プロットされた監視データの全てを包含する領域であって、該領域が最大となるように監視データ同士を結ぶことで形成される領域である。予測可能監視データ領域を設定した後、工程（Ｑ）を実施する。
５％の有意水準で有意であると判断された監視データが０または１種類である場合（図１の無相関検定における「Ｎｏ」）、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する。

［工程（Ｑ）］
工程（Ｑ）では、セメントクリンカー製造における監視データの実測値と工程（Ｐ）で設定された座標空間を用いて、セメントクリンカー製造運転における監視データの実測値と工程（Ｏ）の学習済みのニューラルネットワークによって、セメントクリンカーの製造条件の予測を高い精度で行うことができるか否かを判定することができる。
セメントクリンカーの製造条件の予測に使用される、セメントクリンカー製造運転における監視データの実測値が、工程（Ｑ）で設定した予測可能監視データ領域に含まれる場合（図１の座標判定における「Ｙｅｓ」）、セメントクリンカーの製造条件の予測を高い精度で行うことができると判断し、工程（Ｋ）を実施する。セメントクリンカー製造運転における監視データの実測値が、上記予測可能監視データ領域に含まれない場合（図１の座標判定における「Ｎｏ」）、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する。
なお、セメントクリンカー製造運転における監視データが、工程（Ｐ）で設定された座標空間の座標軸として用いられていない種類の監視データの実測値（無相関検定において５％の有意水準で有意であると判断されなかった監視データの種類）を含む場合、当該の座標軸として用いられていない種類の監視データからは、監視データの実測値に何ら制限は与えない。

［工程（Ｒ）］
工程（Ｒ）では、学習済みのニューラルネットワークの入力層に、セメントクリンカー製造運転における監視データの実測値を入力して、学習済みのニューラルネットワークの出力層から、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力することで、セメントクリンカーの製造条件を予測することができる。

本発明の予測方法によれば、ニューラルネットワークの入力層に入力する、セメントクリンカー製造運転における監視データの実測値から、高い精度でセメントクリンカーの製造条件（プレヒータ原料中のＳＯ_３量）を予測することができる。

ニューラルネットワークは、予測の精度を高い状態に維持するために、評価データの推測値と、該推測値に対応する実測値の乖離の大きさを定期的に点検し、その点検結果に基づいて、ニューラルネットワークを更新することが好ましい。更新の周期は、好ましくは１か月に一回、より好ましくは１週間に一回、特に好ましくは１日に一回である。

本発明のセメントクリンカーの製造条件の予測方法によれば、ニューラルネットワークを用いることによって、監視データを入力するだけで、評価データ（プレヒータ原料中のＳＯ_３量）の推測値を、１時間以内に得ることができる。
また、得られた評価データの推測値に基づいて、セメントクリンカー製造の途中における、サスペンションプレヒータでのコーティングトラブルの発生の可能性を早期に察知し、原料工程及び焼成工程等における諸条件の最適化を行うことにより、セメントクリンカー製造プロセスを安定させることができる。
具体的には、プレヒータ原料中のＳＯ_３量の推測値に増加等の異常が認められた場合、原料の調合、ロータリーキルン及び仮焼炉（サスペンションプレヒータ）の焼成条件の調整等を行うことで、プレヒータ原料中のＳＯ_３量を目的のものにすることができる。
また、評価データの推測値に基いて、セメントクリンカーの製造工程の管理目標値を修正することも可能である。例えば、プレヒータ原料中のＳＯ_３量が目標値に達しないと予測される場合、燃料炭種のＳＯ_３量を分析して、最適な燃料炭種のＳＯ_３量の管理目標値を確認することで、セメントクリンカー製造プロセスを安定させることができる。

さらに、セメントクリンカーの製造工程を制御するコンピュータと、本発明のセメントクリンカーの製造条件の予測方法を実施するために用いるコンピュータを接続することによって、評価データに基づいて監視データを人為的に変動させるための制御システムを自動化することもできる。
本発明において、ニューラルネットワークによる演算を行うためのソフトウェアとしては、例えば、ＯＬＳＯＦＴ社製の「Ｎｅｕｒａｌ Ｎｅｔｗｏｒｋ Ｌｉｂｒａｒｙ」（商品名）等が挙げられる。

以下、実施例により本発明を説明する。
［実施例１］
学習用のサンプルとして、異なる４７点のサンプリング時間における、プレヒータボトムサイクロンから採取したプレヒータ原料のＳＯ_３量を、蛍光Ｘ線分析法を用いて測定して、学習データ（評価データの実測値）とした。
また、上記４７点のサンプリング時間における、ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度を測定して、学習データ（監視データの実測値）とした。
また、モニター用のサンプルとして、前記４７点のサンプリング時間とは異なる５点のサンプリング時間における、プレヒータ原料のＳＯ_３量を測定して、モニターデータ（評価データの実測値）とした。
さらに、上記５点のサンプリング時間における、ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度を測定して、モニターデータ（監視データの実測値）とした。

上記学習データを用いてニューラルネットワークの学習を行った。ニューラルネットワークとしては、入力層、中間層及び出力層を有する階層型のニューラルネットワークを用いた。ニューラルネットワークの初期学習回数を１０，０００回に設定した。得られたニューラルネットワークを用いて、σ_Ｌとσ_Ｍを算出したところ、σ_Ｌとσ_Ｍの関係はσ_Ｌ＜σ_Ｍであった。その後、ニューラルネットワークを初期化し、上記学習データとモニターデータを用いて、ニューラルネットワークの学習を、前記学習回数に０．９５を乗じた数の学習回数（端数切捨て）行うことを、学習後のニューラルネットワークを用いて算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍとなるまで繰り返した。
σ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍとなった後、解析度判定値を算出したが、予め定めた第一の設定値である６％未満とはならず、最も小さい解析度判定値は８．０％であった。最も小さい解析度判定値を算出した際のニューラルネットワークを学習済みのニューラルネットワークとした。
なお、第二の解析度判定値の設定値は２０％とした。

学習データの評価データの実測値であるプレヒータ原料のＳＯ_３量と、学習データの監視データの実測値であるボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度の各々について、無相関検定を行った。学習データの評価データの実測値と、各々の学習データの監視データの実測値との相関関係を表１に示す。
学習データの監視データの実測値のうち、５％の有意水準（相関係数が０．２７８を超えるもの、または、−０．２７８未満のもの）で有意であった学習データは、キルンＥＰ出口ガスのＳＯ_２濃度と塩素バイパス抽気ガスのＳＯ_ｘ濃度の２つであった。

上記無相関検定より有意であると判定されたキルンＥＰ出口ガスのＳＯ_２濃度および塩素バイパス抽気ガスのＳＯ_ｘ濃度について、キルンＥＰ出口ガスのＳＯ_２濃度をｘ軸とし、塩素バイパス抽気ガスのＳＯ_ｘ濃度をｙ軸とする座標空間に、学習データ（監視データの実測値）におけるキルンＥＰ出口ガスのＳＯ_２濃度と塩素バイパス抽気ガスのＳＯ_ｘ濃度のデータをプロットした。プロットされた学習データ同士を結ぶことで形成される学習データの全てを包含する領域であって、該領域の面積が最大となるよう形成される領域を予測可能監視データ領域として設定した。（図２ａ）

得られた学習済みのニューラルネットワークの学習に使用された、監視データの実測値（ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）について、それぞれ平均値と平均２乗誤差を算出した。
また、学習済みのニューラルネットワークを用いて、セメントクリンカーの製造条件を予測するためのサンプル（以下、「予測用サンプル」ともいう。）として、上記の学習用サンプル及びモニター用サンプルとは異なる１点のサンプリング時間における、ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度を測定して、予測用サンプルの監視データの実測値とした。
結果を表２に示す。なお、本明細書中、ｐｐｍは質量基準である。

予測用サンプルの監視データ（ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）の実測値のうち、仮焼炉出口ガスのＯ_２濃度の実測値が、学習データの仮焼炉出口ガスのＯ_２濃度の実測値の平均値±平均２乗誤差（σ）の数値範囲外の数値であった。具体的には、予測用サンプルの仮焼炉出口ガスのＯ_２濃度の実測値（２．６質量％）は、学習データの仮焼炉出口ガスのＯ_２濃度の実測値の平均値（３．４質量％）から、平均２乗誤差（σ）を減算した数値（２．９質量％）よりも小さいものである。
このため、上記学習済みのニューラルネットワークと、セメントクリンカーの製造条件における監視データの実測値を用いて、評価データの推測値（プレヒータ原料のＳＯ_３量）を予測した場合、予測精度が低くなるおそれがある。

そこで、仮焼炉出口ガスのＯ_２濃度の実測値が、上記学習データの実測値の平均値±平均２乗誤差（σ）の範囲外の数値であって、かつ、上記学習データ又はモニターデータとは異なる新たな学習データ５個（前回のニューラルネットワークの学習に使用した学習データ群の約２５％に相当する個数）を新たな学習データとして、上記４７個の学習データに追加して、合計５２個の学習データを、新たな学習データ群として使用して、上述したニューラルネットワークの学習と同様に、ニューラルネットワークの学習を再度行った。
なお、新たな学習用データの仮焼炉出口ガスのＯ_２濃度以外の監視データ（ボトムサイクロンの出口温度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）の実測値は、学習データ（新たな学習データを追加する前の学習データ）の仮焼炉出口ガスのＯ_２濃度以外の監視データ（ボトムサイクロンの出口温度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）の実測値の平均値±平均２乗誤差（σ）の範囲外の数値と範囲内の数値が混在するものであった。
再度のニューラルネットワークの学習に使用した学習データの仮焼炉出口ガスのＯ_２濃度（監視データ）の実測値の平均値は３．４９質量％であり、平均２乗誤差（σ）は０．５６質量％であり、予測サンプルの仮焼炉出口ガスのＯ_２濃度の実測値は、依然として、学習データの仮焼炉出口ガスのＯ_２濃度の実測値の平均値±平均２乗誤差（σ）の数値範囲外の数値であった。

そこで、仮焼炉出口ガスのＯ_２濃度（監視データ）の実測値が、上記学習データの仮焼炉出口ガスのＯ_２濃度の実測値の平均値±σの範囲外の数値であって、かつ、前回のニューラルネットワークの学習に使用した学習データ及びモニターデータとは異なる新たな学習データ１７個（追加した学習データの合計個数（２２個）は、最初のニューラルネットワークの学習に使用した学習データの個数の４７％に相当する個数）を、新たな学習データとして、上記５２個の学習データに追加した合計６９個の学習データを、新たな学習データ群として使用して、上述したニューラルネットワークの学習と同様に、ニューラルネットワークの学習を再度行った。
再度のニューラルネットワークの学習に使用した学習データの監視データの実測値（ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）について、それぞれ平均値と平均２乗誤差を算出した。結果を表３に示す。なお、表３において、参考として予測用サンプルの監視データの実測値を示す。
再度のニューラルネットワークの学習に使用した学習データの仮焼炉出口ガスのＯ_２濃度（監視データ）の実測値の平均値は３．３９質量％であり、平均２乗誤差（σ）は０．８８質量％であり、予測用サンプルの仮焼炉出口ガスのＯ_２濃度（監視データ）の実測値は、学習データの仮焼炉出口ガスのＯ_２濃度の実測値の平均値±平均２乗誤差（σ）の範囲内の数値であった。

ニューラルネットワークの学習に使用した学習データ（６９個）の評価データの実測値であるプレヒータ原料のＳＯ_３量と、監視データの実測値であるボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度の各々について、無相関検定を行った。
プレヒータ原料のＳＯ_３量と、ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度の各々の相関係数を表４に示す。
学習データの監視データの実測値のうち、５％の有意水準（相関係数が０．２３５を超えるもの、または、−０．２３５未満のもの）で有意であった学習データは、前回の無相関検定同様にキルンＥＰ出口ガスのＳＯ_２濃度と塩素バイパス抽気ガスのＳＯ_ｘ濃度であった。

上記無相関検定より有意であると判定されたキルンＥＰ出口ガスのＳＯ_２濃度および塩素バイパス抽気ガスのＳＯ_ｘ濃度について、キルンＥＰ出口ガスのＳＯ_２濃度をｘ軸とし、塩素バイパス抽気ガスのＳＯ_ｘ濃度をｙ軸とする座標空間に学習データ（監視データの実測値）におけるキルンＥＰ出口ガスのＳＯ_２濃度と塩素バイパス抽気ガスのＳＯ_ｘ濃度のデータをプロットした。
プロットした学習データ同士を結ぶことで形成される学習データの全てを包含する領域であって、該領域の面積が最大となるよう形成される領域を予測可能監視データ領域として再設定した。（図２ｂ）
上記座標空間に、予測用サンプルの監視データ（キルンＥＰ出口ガスのＳＯ_２濃度、塩素バイパス抽気ガスのＳＯ_ｘ濃度）の実測値をプロットしたところ、予測可能監視データ領域範囲内であった。

得られた学習済みのニューラルネットワークの入力層に、予測用サンプルの監視データ（ボトムサイクロンの出口温度、仮焼炉出口ガスのＯ_２濃度、仮焼炉出口ガスのＣＯ濃度、キルンＥＰ出口ガスのＳＯ_２濃度、キルンＥＰ出口ガスのＮＯ_ｘ濃度、プレヒータトップサイクロン出口ガスのＣＯ濃度、及び塩素バイパス抽気ガスのＳＯ_ｘ濃度）の実測値を入力し、プレヒータ原料のＳＯ_３量の推測値を得た。該推測値は４．２４±０．６９（偏差は３σを示す。）質量％であった。
また、予測用サンプルのプレヒータ原料のＳＯ_３量の実測値は４．４４質量％であった。

［比較例１］
単回帰分析により、プレヒータ原料のＳＯ_３量と塩素バイパス抽気ガスのＳＯ_ｘ濃度の関係より、予測用サンプルのプレヒータ原料のＳＯ_３量を出力したところ、推測値は６．０２±２．７９（偏差は３σを示す。）質量％であった。

実施例１で得られた評価データの推測値（４．２４±０．６９質量％）と、比較例１で得られた評価データの推測値（６．０２±２．７９質量％）を比較すると、実施例１で得られた推測値は、比較例１で得られた推測値よりも、実測値（４．４４質量％）に近く、かつ、偏差が小さいことから、より信頼性が高いことがわかる。

Claims (5)

1. 入力層及び出力層を有するニューラルネットワークを用いたセメントクリンカーの製造条件の予測方法であって、上記入力層は、セメントクリンカーの製造運転における監視データの実測値を入力するためのものであり、上記出力層は、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力するためのものであり、
   上記監視データが、セメントクリンカーの原料に関するデータ、セメントクリンカーの焼成条件に関するデータ、及びセメントクリンカーに関するデータの中から選ばれる一種以上であり、
   （Ａ）ニューラルネットワークの学習に使用する、監視データの実測値と評価データの実測値の組み合わせである学習データの複数からなる学習データ群の初期設定を行う工程と、
   （Ｂ）学習回数の初期設定を行う工程と、
   （Ｃ）設定された学習データ群を用いて、ニューラルネットワークの学習を、前工程で設定された学習回数行う工程と、
   （Ｄ）学習データの監視データの実測値を、直近の工程（Ｃ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）、及び、ニューラルネットワークの学習結果の信頼性を確認するための監視データの実測値と評価データの実測値の組み合わせであるモニターデータの中の監視データの実測値を、直近の工程（Ｃ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの中の評価データの実測値との平均２乗誤差（σ_Ｍ）を算出し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍである場合、工程（Ｅ）を実施し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ＜σ_Ｍである場合、工程（Ｆ）を実施する工程と、
   （Ｅ）直近の工程（Ｂ）で設定された学習回数および再設定された直近のニューラルネットワークの学習回数のいずれの学習回数よりも大きい学習回数を、工程（Ｂ）において新たな学習回数として再設定し、再度工程（Ｃ）〜（Ｄ）を実施する工程と、
   （Ｆ）直近のニューラルネットワークの学習で実施された学習回数を減らした学習回数を、新たな学習回数として再設定する工程と、
   （Ｇ）設定された学習データ群を用いて、ニューラルネットワークの学習を直近の工程（Ｆ）で設定された学習回数行う工程と、
   （Ｈ）学習データの監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）、及び、モニターデータの中の監視データの実測値を、直近の工程（Ｇ）の学習で得られたニューラルネットワークの入力層に入力して得られた評価データの推測値とモニターデータの中の評価データの実測値との平均２乗誤差（σ_Ｍ）を算出し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ≧σ_Ｍである場合、工程（Ｊ）を実施し、算出されたσ_Ｌとσ_Ｍの関係がσ_Ｌ＜σ_Ｍである場合、工程（Ｉ）を実施する工程と、
   （Ｉ）直近の工程（Ｇ）におけるニューラルネットワークの学習回数が予め設定された回数を超えている場合、再度工程（Ｆ）〜（Ｈ）を実施し、直近の工程（Ｇ）におけるニューラルネットワークの学習回数が予め定めた数値以下の場合、工程（Ｍ）を実施する工程と、
   （Ｊ）下記式（１）を用いて解析度判定値を算出し、該解析度判定値が予め定めた第一の設定値未満である場合、ニューラルネットワークの学習を終了し工程（Ｋ）を実施する工程と、上記解析度判定値が予め定めた第一の設定値以上である場合、工程（Ｍ）を実施する工程と、
   （Ｋ）セメントクリンカーの製造運転における監視データの実測値が、設定された学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲内の数値である場合には、工程（Ｒ）を実施し、範囲外の数値である場合には、工程（Ｌ）を実施する工程と、
   （Ｌ）監視データの実測値が、設定された学習データ群における、監視データの実測値の平均値±平均２乗誤差（σ_Ｇ）の範囲外の数値である学習データであって、上記学習データ群に含まれない１個以上の学習データを、上記学習データ群に追加して新たな学習データ群に設定した後、工程（Ｃ）以降を実施する工程 と
   （Ｍ）工程（Ｂ）を実施した回数の大きさについての判定を行い、該回数が予め設定した回数以下である場合は、工程（Ｎ）を実施し、該回数が予め設定した回数を超える場合には工程（Ｏ）を実施する工程と、
   （Ｎ）学習条件の初期化を行って、再度工程（Ｂ）〜（Ｊ）を実施する工程と、
   （Ｏ）工程（Ｊ）において算出した全ての解析度判定値のうち、最も小さい解析度判定値が予め定めた第二の設定値未満である場合、最も小さい解析度判定値を得ることができた工程（Ｊ）におけるニューラルネットワークを、学習済みのニューラルネットワークとした後、工程（Ｐ）を実施し、最も小さい解析度判定値が予め定めた第二の設定値以上である場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
   （Ｐ）工程（Ｊ）において算出した全ての解析度判定値のうち、最も小さい解析度判定値を得ることができた工程（Ｊ）において、学習データとして使用した監視データの実測値と評価データの実測値の組み合わせについて無相関検定を行い、５％の有意水準で有意であると判断された監視データの種類が２種以上である場合、５％の有意水準で有意であると判断された監視データの全種類を座標軸とする座標空間に学習データとして使用した監視データの実測値をプロットし、座標空間において、プロットされた監視データ同士を結ぶことで形成される監視データの全てを包含する領域であって、該領域が最大となるように監視データ同士を結ぶことで形成される領域を、予測可能監視データ領域として設定した後、工程（Ｑ）を実施し、５％の有意水準で有意であると判断された監視データが０または１種類である場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
   （Ｑ）セメントクリンカー製造運転における監視データの実測値が、工程（Ｐ）で設定した予測可能監視データ領域に含まれる場合、セメントクリンカーの製造条件の予測を高い精度で行うことができると判断し、工程（Ｋ）を実施し、セメントクリンカー製造運転における監視データの実測値が、上記予測可能監視データ領域に含まれない場合、セメントクリンカーの製造条件を予測することはできないと判断して予測を終了する工程と、
   （Ｒ）学習済みのニューラルネットワークの入力層に、セメントクリンカー製造運転における監視データの実測値を入力して、上記ニューラルネットワークの出力層から、プレヒータ原料中のＳＯ_３量に関する評価データの推測値を出力して、セメントクリンカーの製造条件を予測する工程と、
   を含むことを特徴とするセメントクリンカーの製造条件の予測方法。
   

   

   （上記式（１）中、学習データの平均２乗誤差（σ_Ｌ）とは、学習データの監視データの実測値を学習後のニューラルネットワークの入力層に入力して得られた評価データの推測値と学習データの評価データの実測値との平均２乗誤差（σ_Ｌ）である。評価データの推測値の平均値とは、学習データの監視データの実測値を学習後のニューラルネットワークの入力層に入力して得られた評価データの推測値の平均値である。）
2. 上記工程（Ｌ）で追加する学習データの数が、設定された学習データ群における学習データの数の８％以上である請求項１に記載のセメントクリンカーの製造条件の予測方法。
3. 上記解析度判定値の予め定めた第一の設定値が１０％以下であり、上記解析度判定値の予め定めた第二の設定値が上記第一の設定値よりも大きくかつ３０％以下である、請求項１または２に記載のセメントクリンカーの製造条件の予測方法。
4. 上記ニューラルネットワークが、上記入力層と上記出力層の間に中間層を有する階層型のニューラルネットワークである、請求項１〜３のいずれか１項に記載のセメントクリンカーの製造条件の予測方法。
5. 上記監視データの値を人為的に変動させて得られた上記評価データの推測値に基づいて、プレヒータ原料中のＳＯ_３量を最適化する、請求項１〜４のいずれか１項に記載のセメントクリンカーの製造条件の予測方法。
   
   

Images