(57)【特許請求の範囲】
【請求項１】
定在波変動を表現する表面波の基礎方程式を境界要素法を用いて空間離散化することにより導出した定在波変動のモデルに、鋳型内溶鋼の体積変動に対応する変数を追加することにより得られた、鋳型内幅方向の湯面高さ分布の変動モデルを用いた逐次計算アルゴリズムによって、連続鋳造鋳型内の任意位置に配置された複数の湯面レベル計による観測値を用いて、鋳型内湯面変動の体積変動成分と定在波変動成分とをリアルタイムに求めることにより、鋳型内幅方向の湯面高さ分布を算出することを特徴とする、連続鋳造鋳型内の湯面変動の状態推定方法。
【請求項２】
前記逐次計算アルゴリズムとして、カルマンフィルタを用いることを特徴とする、請求項１に記載の連続鋳造鋳型内の湯面変動の状態推定方法。
【請求項３】
連続鋳造鋳型内の任意位置に配置された複数の湯面レベル計と、
定在波変動を表現する表面波の基礎方程式を境界要素法を用いて空間離散化することにより導出した定在波変動のモデルに、鋳型内溶鋼の体積変動に対応する変数を追加することにより得られた、鋳型内幅方向の湯面高さ分布の変動モデルを用いた逐次計算アルゴリズムによって、前記複数の湯面レベル計による観測値を用いて、鋳型内湯面変動の体積変動成分と定在波変動成分とをリアルタイムに求めることにより、鋳型内幅方向の湯面高さ分布を算出する、湯面形状推定部と
を含むことを特徴とする、連続鋳造鋳型内の湯面変動の状態推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、連続鋳造鋳型内の湯面変動の状態を推定する方法に関する。ここで、「湯面変動の状態」には、湯面全体の上下動や、定在波が生じている湯面状態が含まれる。
【背景技術】
【０００２】
連続鋳造機には、湯面レベル制御を目的として、通常、一つの渦流センサ（以下において、「湯面レベル計」とも称する。）が設置されている。この渦流センサの観測値には、溶鋼収支の不均衡による湯面全体の上下動（鋳型内溶鋼の体積変動。以下において、「体積変動」と称することがある。）に加えて、湯面の波立ちによる定在波や、観測雑音等の外乱が含まれている。体積変動成分、定在波変動成分、および、外乱成分が重畳している渦流センサの観測値に基づいて注湯制御手段を操作すると、系への外乱となり湯面変動を増幅するおそれがある。それゆえ、高精度に湯面レベル制御を行うためには、湯面レベル制御の対象となる鋳型内溶鋼の体積変動だけを精度良く抽出する必要がある。ところが、通常、一つの渦流センサの観測値から体積変動成分だけを抽出することは困難であるため、鋳型内溶鋼の体積変動を抑制することは困難である。さらに、高品質な鋳片製造のためには、鋳型内幅方向における湯面高さ分布（以下において、「湯面形状」と称することがある。）の変動を抑制することが必要である。しかしながら、一つの渦流センサの観測値から制御対象となる湯面形状の変動を検出することは困難であるため、その抑制のための適正な操業条件を把握することは困難である。
【０００３】
連続鋳造機における湯面レベル制御に関する技術として、例えば特許文献１には、連続鋳造機の鋳型内の所定位置における湯面レベル変動の定在波変動分の定在波周波数を、鋳型幅から算出し、その周波数で定在波をｓｉｎ関数とｃｏｓ関数で記述して、その記述された定在波の係数を測定した湯面レベル変動からオンライン推定することにより、定在波変動分を求める鋳型内の湯面定在波変動検出方法が開示されている。また、特許文献１には、上記鋳型内の湯面定在波変動検出方法を用いて定在波変動分を算出し、算出された定在波変動分を湯面レベル変動から差し引いて体積変動分を求めて、その体積変動分を用いて湯面レベルを制御する鋳型内の湯面レベル制御方法が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−２４１１５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に記載の技術では、定在波変動を、三角関数を用いて表現し、その振幅と位相を未知パラメータとしてリアルタイムに推定している。しかしながら、この方法では、湯面レベル計の観測点における定在波の振幅と位相が求まるにとどまり、湯面形状を直接算出していないため、鋳型内の湯面形状をリアルタイムに推定することはできなかった。
【０００６】
そこで本発明は、鋳型内の湯面形状をリアルタイムに推定することが可能な、連続鋳造鋳型内の湯面変動の状態推定方法を提供することを課題とする。
【課題を解決するための手段】
【０００７】
本発明者らは、鋭意検討の結果、鋳型内湯面変動を体積変動と定在波変動との重ね合わせとして表現し、定在波変動を表面波の基礎方程式により表現する定在波変動のモデルを用いた逐次計算アルゴリズムにより、複数の湯面レベル計による鋳型内の湯面レベル観測結果を用いて、鋳型内湯面変動として体積変動と定在波変動とをリアルタイムに推定可能であることを見出した。定在波変動のモデルとしては、例えば、表面波の基礎方程式を、境界要素法を用いて空間離散化することにより構築した、常微分方程式による状態空間モデルを用いることができる（以下において、このモデルを「表面波モデル」と称する。）。湯面変動には、定在波変動に加えて体積変動が重畳して観測されるので、本発明では、体積変動に対応する変数を表面波モデルに追加する（本発明において、この拡大したモデルを「湯面形状変動モデル」と称する。）。また、逐次計算アルゴリズムとしては、例えば、カルマンフィルタを用いる。従来の単一の湯面レベル計の観測値からは、鋳型内の任意位置での湯面高さを検出することは困難であったが、上記形態にすることにより、複数の湯面レベル計による鋳型内の湯面レベル観測結果を用いて、定在波変動と体積変動とをリアルタイムに検出することが可能になるので、鋳型内の任意位置での湯面高さを算出できるようになる。
本発明は、このような知見に基づいて完成させた。以下、本発明について説明する。
【０００８】
本発明は、定在波変動を表面波の基礎方程式により表現する定在波変動のモデルを用いた逐次計算アルゴリズムによって、連続鋳造鋳型内の任意位置に配置された複数の湯面レベル計による観測値を用いて、鋳型内湯面変動の体積変動成分と定在波変動成分とをリアルタイムに求めることにより、上記任意位置における湯面高さを算出することを特徴とする、連続鋳造鋳型内の湯面変動の状態推定方法である。
【０００９】
体積変動成分と定在波変動成分とをリアルタイムに求めることにより、鋳型内幅方向における湯面高さ分布をリアルタイムに推定することが可能になる。その結果、鋳型内の湯面形状をリアルタイムに推定することが可能になる。
【００１０】
また、上記本発明において、定在波変動のモデルとして、表面波の基礎方程式を、境界要素法を用いて空間離散化することにより導出したモデルを用いることができる。このような形態であっても、鋳型内の湯面形状をリアルタイムに推定することが可能である。
【００１１】
また、上記本発明において、逐次計算アルゴリズムとして、カルマンフィルタを用いることができる。このような形態であっても、鋳型内の湯面形状をリアルタイムに推定することが可能である。
【発明の効果】
【００１２】
本発明によれば、鋳型内の湯面形状をリアルタイムに推定することが可能な、連続鋳造鋳型内の湯面変動の状態推定方法を提供することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明にかかる連続鋳造鋳型内の湯面変動の状態推定方法を実施可能なシステムの構成例を説明する図である。
【図２】本発明にかかる連続鋳造鋳型内の湯面変動の状態推定方法を説明するフローチャートである。
【図３】境界要素法による空間離散化を模式的に説明する図である。
【図４】実施例における湯面レベル計の取り付け位置を説明する図である。
【図５】鋳型幅方向中央付近に設置した湯面レベル計の実測値および本発明による推定値を説明する図である。
【図６】鋳型幅方向西側に設置した湯面レベル計の実測値および本発明による推定値を説明する図である。
【図７】鋳型幅方向東側に設置した湯面レベル計の実測値および本発明による推定値を説明する図である。
【図８】本発明により推定した体積変動成分を説明する図である。
【図９】本発明により推定した定在波変動成分を説明する図である。
【図１０】６．２秒時点における湯面形状の実測値および本発明による湯面形状の推定値を説明する図である。
【図１１】６．９秒時点における湯面形状の実測値および本発明による湯面形状の推定値を説明する図である。
【発明を実施するための形態】
【００１４】
本発明の実施の形態について、以下に説明する。なお、以下に説明する形態は本発明の例示であり、本発明は以下に説明する形態に限定されない。
【００１５】
はじめに、本発明で構築する問題を説明する。鋳型は、幅方向の長さに対して厚み方向の長さが十分に短いので、本発明では厚み方向の湯面変動は考慮しない。また、本発明において、自由表面上では圧力一定とし、湯面変動の振幅は幅方向の長さに対して微小であると仮定する。
連続鋳造時に、鋳型内に浸漬ノズルで供給される溶湯が、浸漬ノズル吐出孔に付着物がなく鋳型の幅方向の中心を境にして左右均等に供給されている場合、上記仮定は十分妥当な仮定である。また、鋳型の幅方向の長さが厚み方向の長さの５倍以上である場合、厚み方向の振幅が幅方向の振幅に対して微小であるという仮定は、十分妥当な仮定である。
また、本発明では、連続鋳造鋳型を底面が十分深い容器であるとみなし、鋳型内の湯面変動を矩形容器内表面波の基礎方程式により表現する。実際には鋳型は底面を持たないが、溶鋼プール深さは湯面変動の振幅に対して十分大きいので、底の深い容器として近似してよい。底面の深さの選び方については後述する。
【００１６】
１．表面波の基礎方程式
粘性のない非圧縮流体の渦なし流れを仮定すると、矩形容器内の表面波の基礎方程式は、式（１）〜（４）で表現される（神部勉  編著、流体力学、裳華房、ｐ．１２２〜１２６）。
【００１７】
【数１】
0006372217.tif 000002
【００１８】
【数２】
0006372217.tif 000003
【００１９】
【数３】
0006372217.tif 000004
【００２０】
【数４】
0006372217.tif 000005
ただし、Ωは鋳型内における液体の存在領域であり、後述する図３に示したように、鋳型の幅方向（ｘ方向）が−Ｗ／２〜Ｗ／２、且つ、鋳型の深さ方向（ｚ方向）が０〜ｈの２次元平面とする。また、Γｆは液体領域と空気との境界線を表す自由表面であり、φは速度ポテンシャル（ｍ２／ｓ）であり、ηは自由表面の高さ（ｍ）であり、ｇは重力加速度（ｍ／ｓ２）であり、ｈは底面の深さ（ｍ）であり、Ｗは水平方向の鋳型幅（ｍ）である。
【００２１】
２．境界要素法による基礎方程式の空間離散化
ここでは、上記の基礎方程式（１）〜（４）を、境界要素法を用いて空間離散化する場合について説明する。境界要素法とは、元の境界値問題を境界積分方程式に変換した後に領域の境界を有限個の要素に離散化し、境界上の変数のみを未知変数とする連立方程式を解く問題に帰着させる手法である。境界要素法の利点は境界上の変数のみが未知変数になることにより、求めるべき未知変数の数が少なくなることである。以下に示す方法では、自由表面上の速度ポテンシャルφと自由表面の高さ（湯面高さ）ηのみが未知変数となる。ここで、自由表面上の高さの具体的な値が湯面高さηである。
なお、境界要素法による自由表面および境界上での離散化点数は任意でよいが、ここでは自由表面上では等間隔にＭ点に分割し、他の３つの壁面もそれぞれ等間隔にＭ点ずつ分割するものとする。分割の様子を図３に示す。また、節点（境界要素上の代表点）に関して様々な選び方があるが、要素中心を節点として選べばよい。さらに未知変数φおよびηの値は、各領域（境界要素）上で一定値であると仮定する一定要素近似を用いる。すなわち、自由表面上のＭ個の選点をｘｊ（ｊ＝１、２、…、Ｍ）とすると、各領域上で下記式（５）および（６）が成り立つ。図３に示すように、境界領域を階段関数状に近似するものである。
【００２２】
【数５】
0006372217.tif 000006
【００２３】
【数６】
0006372217.tif 000007
【００２４】
表面波の基礎方程式（１）〜（４）を、境界要素法を用いて空間離散化することにより、自由表面上の点ｘｊ（ｊ＝１、２、…、Ｍ）における、速度ポテンシャルφｊ（ｊ＝１、２、…、Ｍ）および湯面高さηｊ（ｊ＝１、２、…、Ｍ）を状態変数とする２Ｍ次元の状態空間モデル（下記式（７））が導出される。境界要素法による空間離散化の手順については、例えば、「登坂宣好、中山司  著、境界要素法の基礎、日科技連、ｐ．１７０〜１８４、１９８７年）に詳しく述べられている。
【００２５】
【数７】
0006372217.tif 000008
ここで、Ａは境界要素法により空間離散化することで得られる２Ｍ次の正方行列である。なお、自由表面以外の境界上の変数の情報は行列Ａに織り込まれている。行列Ａは底面の深さｈに応じて変化し、ｈを大きくするにつれ行列Ａはある一定値に収束する。ｈの値は大きいほど好ましく、要求するシミュレーション精度に応じて適当なｈの値を選定すればよい。具体的には、鋳型幅Ｗに対し、鋳型深さｈは２倍以上とすることが望ましい。例えば、後述する本発明の実施例では、鋳型幅Ｗ＝１．４２５（ｍ）に対し、鋳型深さｈ＝３．０（ｍ）とした。
【００２６】
３．湯面形状変動モデル
鋳型内の湯面変動では、定在波変動に加えて体積変動も重畳する。そのため、上記式（７）の表面波モデルに体積変動を表現する状態変数ｄを追加した下記式（８）の拡大したモデルを考える。この拡大したモデルを、本発明では湯面形状変動モデルと称する。このモデルは２Ｍ＋１次の状態変数を持つ。
【００２７】
【数８】
0006372217.tif 000009
ただし、Ａａｕｇは下記式（９）で表される。
【００２８】
【数９】
0006372217.tif 000010
【００２９】
湯面レベル計の数をＬ個とすると、レベル計観測値ｙ’＝（ｙ１、ｙ２、…、ｙＬ）’は、定在波変動に体積変動ｄが重畳して観測され、下記式（１０）のように表わされる。
【００３０】
【数１０】
0006372217.tif 000011
ただし、η’＝（η１、η２、…、ηＭ）’である。Ｃは観測により定まるＬ行Ｍ列の行列であり、１’は全ての要素が１であるＬ次のベクトルである。
【００３１】
４．カルマンフィルタのアルゴリズム
ここでは、逐次計算アルゴリズムとしてカルマンフィルタを使用する場合について説明する。「樋口知之  編著、データ同化入門、朝倉書店、ｐ．４７〜５７、２０１１年」に基づいてカルマンフィルタを説明する。
カルマンフィルタは、以下の時変な線形・ガウス状態空間モデルを対象にする。
【００３２】
【数１１】
0006372217.tif 000012
【００３３】
【数１２】
0006372217.tif 000013
ただし、Ｒｎはｎ次元の実数ベクトルの集合であり、Ｒｍはｍ次元の実数ベクトルの集合である。また、Ｆｔは時変のｎ×ｎ行列、Ｇｔは時変のｎ×ｎ行列であり、Ｈｔは時変のｍ×ｎ行列である。また、ｖ’ｔ〜Ｎ（０、Ｑｔ）はＲｎの要素であり、ｗ’ｔ〜Ｎ（０、Ｒｔ）はＲｍの要素である。Ｎ（０、Ｑｔ）は平均が０であり且つ分散共分散行列がＱｔである多次元正規分布を表わし、Ｎ（０、Ｒｔ）は平均が０であり且つ分散共分散行列がＲｔである多次元正規分布を表わす。
カルマンフィルタのアルゴリズムは、状態ベクトルと分散共分散行列の初期条件ｘ’０｜０、Ｖ０｜０を与えた後に、予測とフィルタリングの手順を逐次的に繰り返す。
【００３４】
４．１．予測
時刻ｔ−１における推定値ｘ’ｔ−１｜ｔ−１および上記式（１１）で表わされるモデルを用いて、状態の予測値ｘ’ｔ｜ｔ−１および分散共分散行列の予測値Ｖｔ｜ｔ−１を以下のように計算する。
【００３５】
【数１３】
0006372217.tif 000014
【００３６】
【数１４】
0006372217.tif 000015
ただし、ｘ’ｔ｜ｔ−１は、時刻ｔ−１における、時刻ｔでの状態量の予測値を表わす。
【００３７】
４．２．フィルタリング
分散共分散行列の推定値Ｖｔ｜ｔとカルマンゲインＫｔを以下のように計算する。
【００３８】
【数１５】
0006372217.tif 000016
【００３９】
【数１６】
0006372217.tif 000017
計算したカルマンゲインＫｔ、および、観測データｙ’ｔを用いて、予測値ｘ’ｔ｜ｔ−１を以下の式（１７）のように修正して推定値とする。
【００４０】
【数１７】
0006372217.tif 000018
ただし、ｘ’ｔ｜ｔは時刻ｔにおける状態量の修正値を表わす。
このカルマンフィルタを上記式（８）および式（１０）から構成される湯面形状変動モデルに適用して、湯面形状をリアルタイムに推定する。
【００４１】
図面を参照しつつ、本発明の実施形態についての説明を続ける。図１は本発明を実施可能なシステムの構成例を説明する図である。図１に示した連続鋳造鋳型内の湯面変動状態推定システム１０では、連続鋳造鋳型１における湯面の状態が湯面レベル計２によって計測される。湯面レベル計２によって計測された操業データ３は、湯面形状推定部４へと送られるとともに、操業データベース５および制御量演算部６へと送られる。湯面形状推定部４では、操業データ３を用いて、上記湯面形状変動モデルにより湯面高さを算出する。湯面形状推定部４で算出された湯面高さの算出結果は、出力部７に出力される。出力部７は、操業データベース５に記憶された情報も出力できるように構成されている。出力部７に出力された湯面高さの算出結果は、制御量演算部６へと送られる。制御量演算部６では、湯面高さの算出結果と、操業データ３と、操業データベース５に記憶されているデータとを用いて、湯面変動を低減するための制御量が算出される。制御量演算部６で算出された制御量は、注湯制御手段８へと送られ、湯面変動を低減する制御が行われる。
【００４２】
図２は、本発明の処理の流れを説明するフローチャートである。図２に示したように、本発明により連続鋳造鋳型内の湯面変動の状態を推定する際には、湯面レベル計によって計測されたデータを含む操業データを取得するデータ取得工程（Ｓ２１）が行われる。その後、取得した操業データを用いて、上記湯面形状変動モデルにより湯面高さを算出する湯面高さ算出工程（Ｓ２２）が行われる。Ｓ２２では、例えば、逐次計算アルゴリズムとしてカルマンフィルタが用いられ、上記「４．１．予測」および「４．２．フィルタリング」が逐次的に繰り返される。そして、Ｓ２２で算出された湯面高さの結果は、出力工程（Ｓ２３）で出力される。
このようにして連続鋳造鋳型内の湯面変動の状態を推定する本発明によれば、鋳型内の湯面形状をリアルタイムに知ることができる。それゆえ、これを用いて湯面制御を行うことにより、高品質な鋳片を製造することが可能になる。また、本発明によれば、湯面レベル制御の制御対象である体積変動成分と定在波変動成分とを高精度に分離することが可能なので、湯面レベル制御性能を高めるフィードバック制御への応用が可能となる。このような効果を奏する本発明によれば、鋳片品質を大幅に向上させることが可能になるほか、歩留まりを減少させることが可能になるので、製造コストを大幅に低減することが可能になる。
【実施例】
【００４３】
実施例を参照しつつ、本発明についてさらに説明を続ける。
【００４４】
１．実施例の設定条件
（１）領域の設定
先にも述べたように、境界要素法による自由表面および境界上での離散化点数は任意で良いが、ここでは自由表面上では等間隔にＭ点分割し、他の３つの壁面もそれぞれ等間隔にＭ点ずつ分割するものとした。この分割の様子を図３に示す。本実施例では、要素中心を節点として選び、一定要素近似を用いた。また、自由表面上の空間離散化の点数Ｍは１４点とした。鋳型幅Ｗは１．４２５（ｍ）であった。上述のように、境界要素法により離散化して得られる式（７）の行列Ａは、底面の深さｈに応じて変化し、ｈを大きくするにつれ行列Ａはある一定値に収束する。本実施例では、鋳型深さｈを３．０（ｍ）とした。
【００４５】
（２）観測モデルの設定
本実施例では、鋳型厚み０．２７（ｍ）、鋳型幅Ｗ１．４２５（ｍ）である鋳型に３基の湯面レベル計を設置し、湯面レベル計直下の湯面高さを観測した。湯面レベル計の設置個所は、鋳型幅方向両端、および、鋳型幅方向中央から０．２６ｍ東側の位置とした。図４に、湯面レベル計の取り付け位置を示す。
それぞれの湯面レベル計により、湯面レベル計の設置位置に最も近い湯面上の節点を観測した。具体的には、１４点の節点のうち、左から１番目の節点（ｘ＝−０．６６（ｍ））、左から１０番目の節点（ｘ＝０．２５（ｍ））、および、左から１４番目の節点（ｘ＝０．６６（ｍ））の３点の湯面高さを観測した。この場合、上記式（１０）における行列Ｃは、３行１４列の行列であり、１行１列成分、２行１０列成分、および、３行１４列成分に１の値を持ち、その他の成分の値は０である。すなわち、観測値ｙ’（ｔ）＝（ｙ１（ｔ）、ｙ２（ｔ）、ｙ３（ｔ））’と湯面高さη’（ｔ）＝（η１（ｔ）、η２（ｔ）、…、η１４（ｔ））’は、下記式（１８）で表わされる。
【００４６】
【数１８】
0006372217.tif 000019
式（１８）は、以下の式（１９）のように表わすことができる。
【００４７】
【数１９】
0006372217.tif 000020
ただし、観測行列Ｈは、下記式（２０）で表される３行２９列の行列である。
【００４８】
【数２０】
0006372217.tif 000021
以下において、状態変数Ｘ（ｔ）≡［φ１（ｔ）、φ２（ｔ）、…、φ１４（ｔ）、η１（ｔ）、η２（ｔ）、…η１４（ｔ）、ｄ（ｔ）］’と表記する。
【００４９】
（３）時間離散化の設定
上記式（８）の湯面形状変動モデルを計算機でシミュレーションするために、適当な時間間隔で時間離散化を行う。この離散化の方法としては、たとえば双一次変換を用いればよい。双一次変換を用いることで、上記式（８）の離散時間版として下記式（２１）が得られる。ｙ’（ｔ）は、、式（２１）のＸ（ｔ）を用いて下記式（２２）で表わすことができる。
【００５０】
【数２１】
0006372217.tif 000022
【００５１】
【数２２】
0006372217.tif 000023
ここに、ＡＤおよびＨＤは、双一次変換による時間離散化により変換された行列である。双一次変換においては、下記式（２３）および式（２４）により変換される。
【００５２】
【数２３】
0006372217.tif 000024
【００５３】
【数２４】
0006372217.tif 000025
ここに、Ｉ２９は２９時の単位行列を表わし、Δｔは時間刻み幅を表わす。また、Ａａｕｇは上記式（９）で定義される行列である。
【００５４】
（４）カルマンフィルタの設定
カルマンフィルタは確率的なシステムを対象とするので、上記式（２１）および式（２２）に正規雑音ｖ’ｔ、ｗ’ｔが加わると仮定し、下記式（２５）および式（２６）の確率システムとしてモデル化する。ここで、ｖ’ｔは湯面変動に加わる外乱成分であり、吐出流の乱れにより生じる湯面変動の乱れに相当する。また、ｗ’ｔは湯面レベル計の測定誤差を表わす外乱成分である。ｖ’ｔおよびｗ’ｔは、いずれも湯面形状変動モデルで考慮していない外乱成分を表現する変数である。カルマンフィルタの設定では、これら外乱成分は正規分布に従う白色雑音であると仮定する。
【００５５】
【数２５】
0006372217.tif 000026
【００５６】
【数２６】
0006372217.tif 000027
ここで、Ｇは２９次の単位行列である。
【００５７】
システムノイズの分散共分散行列Ｑ、および、観測ノイズの分散共分散行列Ｒは、それぞれ、２９次および３次の正定値対称行列であり、カルマンフィルタにおける調整パラメータである。その値は適切に定める必要があり、試行錯誤的に選定すればよい。本実施例において、システムノイズの分散共分散行列Ｑは対角行列とし、その対角要素Ｑｉｉを、０．００５２（ｍ２）（ｉ＝１、２、…、２８）、および、０．０１２（ｍ２）（ｉ＝２９）とした。また、本実施例において、観測ノイズの分散共分散行列Ｒは対角行列とし、その対角要素Ｒｌｌを、０．００５２（ｍ２）（ｌ＝１、２、３）とした。
【００５８】
本実施例では、カルマンフィルタの予測ステップ（下記式（２７）〜（２８））とフィルタリングステップ（下記式（２９）〜（３１））の操作を交互に行うことにより、逐次的に湯面形状変動モデル内の状態変数を推定した。
【００５９】
＜予測ステップ＞
時刻ｔ−１における推定値Ｘｔ−１｜ｔ−１、および、湯面形状変動モデルを用いて、時刻ｔ−１における、次の時刻ｔの状態予測値Ｘｔ｜ｔ−１、および、分散共分散行列の予測値Ｖｔ｜ｔ−１を、以下のように計算する。
【００６０】
【数２７】
0006372217.tif 000028
【００６１】
【数２８】
0006372217.tif 000029
【００６２】
＜フィルタリングステップ＞
分散共分散行列の推定値Ｖｔ｜ｔ、および、カルマンゲインＫｔを、以下のように計算する。
【００６３】
【数２９】
0006372217.tif 000030
【００６４】
【数３０】
0006372217.tif 000031
予測ステップで計算した推定値ＨＸｔ｜ｔ−１と観測値ｙ’ｔとの誤差をフィードバックすることにより、上記式（２７）の予測値Ｘｔ｜ｔ−１を修正し、推定値Ｘｔ｜ｔを計算する。修正ゲインは、上記式（３０）で求めたカルマンゲインＫｔを用いる。
【００６５】
【数３１】
0006372217.tif 000032
【００６６】
［結果］
本実施例では、実機データを用いて、上記式（２７）〜（３１）で構成したカルマンフィルタのアルゴリズムを適用し、リアルタイムに湯面形状を検出した。３点の湯面レベル計観測値の実績値、および、本発明により算出した推定値の１０秒間の時系列データを、図５〜７に示す。また、本発明により算出した、体積変動成分および定在波変動成分の時系列データを、図８および図９に示す。さらに、リアルタイムに推定した湯面形状を、図１０および図１１に示す。
【００６７】
単一の湯面レベル計の観測値を用いる従来技術では、湯面形状を推定することは困難であったが、図５〜図１１に示したように、本発明によれば、鋳型内の湯面形状をリアルタイムに推定できることが確認された。また、図１０および図１１に示したように、本発明によれば、実績値を高精度に推定できた。
【００６８】
このように、本発明により湯面形状をリアルタイムに推定することができるので、本発明を用いることにより、高品質な鋳片製造のための適正な操業条件を把握することが可能になる。また、本発明によれば、湯面レベル制御の制御対象である体積変動成分と定在波変動成分とを高精度に分離することが可能なので、湯面レベル制御性能を高めるフィードバック制御への応用が可能となる。このような効果を奏する本発明によれば、鋳片品質を大幅に向上させることが可能になるほか、歩留まりを減少させることが可能になるので、製造コストを大幅に低減することが可能になる。
【符号の説明】
【００６９】
１…連続鋳造鋳型
２…湯面レベル計
３…操業データ
４…湯面形状推定部
５…操業データベース
６…制御量演算部
７…出力部
８…注湯制御手段
１０…連続鋳造鋳型内の湯面変動状態推定システム
【図１】
0006372217.tif 000033
【図２】
0006372217.tif 000034
【図３】
0006372217.tif 000035
【図４】
0006372217.tif 000036
【図５】
0006372217.tif 000037
【図６】
0006372217.tif 000038
【図７】
0006372217.tif 000039
【図８】
0006372217.tif 000040
【図９】
0006372217.tif 000041
【図１０】
0006372217.tif 000042
【図１１】
0006372217.tif 000043