(57)【要約】      （修正有）【課題】容器壁内部の温度計算を行わなくとも、容器内に高温物質が存在する状態で、容器壁の厚みを直接的にかつ精度よく推定できる容器壁の厚み推定方法を提供する。【解決手段】容器の外壁面の温度計測点にて計測された温度ｈ(ｔi)が入力されるとともに、温度ｈ(ｔi)を基に熱流束を熱伝導率で除した物理量ｇ(ｔi)を算出する（Ｓ１０３）。そして、残存厚みの仮定値ｌ〜を設定し、変数ｖ(ｌ〜，ｔi)を求める（Ｓ１０４）。次に、式（１０）のＭＡ×ＶＢ＝Ｖｂを解くことにより、Ｂ0(ｌ〜)、Ｂ1(ｌ〜)、・・・、ＢN(ｌ〜)を求め（Ｓ１０５）、式（１１）のｐ(ｌ〜，ｔ)を算出する（Ｓ１０６）。設定値ε以下となるｐ(ｌ〜，ｔ)が得られるまで演算処理を繰り返し（Ｓ１０７）、設定値ε以下となったときの仮定値ｌ〜を容器壁の残存厚みｌとして決定する（Ｓ１０８）。【選択図】図３
2008063593.tif 000002
【特許請求の範囲】
【請求項１】
内壁面と外壁面とに温度差を有する容器における容器壁の厚み推定方法であって、
前記外壁面における温度ｈ(ｔ)及び熱流束を熱伝導率で除した物理量ｇ(ｔ)を取得する手順と、
容器壁の厚みの仮定値ｌ〜を設定する手順と、
内壁面（ｘ＝０）から外壁面方向をｘ軸とし、式（１０１）における容器壁の温度ｕ(ｘ，ｔ)の代替として変数ｖ(ｘ，ｔ)を定義して導入した式（１０２）から、変数ｖ(ｌ〜，ｔi)を算出する手順と、
【数１】
2008063593.tif 000017
変数ｗ(ｘ，ｔ)を定義して導入した式（１０３）から、変数ｗ(ｘ，ｔ)を式（１０４）により表わし、式（１０４）に基づき式（１０５）のＭＡ×ＶＢ＝Ｖｂを解くことによりＶＢを決定する手順と、
【数２】
2008063593.tif 000018
式（１０６）で表わされるｐ(ｌ〜，ｔ)が設定値以下となる仮定値ｌ〜を容器壁の厚みとして決定する手順と、
【数３】
2008063593.tif 000019
を有することを特徴とする容器壁の厚み推定方法。
【請求項２】
内壁面と外壁面とに温度差を有する容器における容器壁の厚み推定装置であって、
前記外壁面における温度ｈ(ｔ)及び熱流束を熱伝導率で除した物理量ｇ(ｔ)を取得する手段と、
容器壁の厚みの仮定値ｌ〜を設定する手段と、
内壁面（ｘ＝０）から外壁面方向をｘ軸とし、式（１０１）における容器壁の温度ｕ(ｘ，ｔ)の代替として変数ｖ(ｘ，ｔ)を定義して導入した式（１０２）から、変数ｖ(ｌ〜，ｔi)を算出する手段と、
【数４】
2008063593.tif 000020
変数ｗ(ｘ，ｔ)を定義して導入した式（１０３）から、変数ｗ(ｘ，ｔ)を式（１０４）により表わし、式（１０４）に基づき式（１０５）のＭＡ×ＶＢ＝Ｖｂを解くことによりＶＢを決定する手段と、
【数５】
2008063593.tif 000021
式（１０６）で表わされるｐ(ｌ〜，ｔ)が設定値以下となる仮定値ｌ〜を容器壁の厚みとして決定する手段と、
【数６】
2008063593.tif 000022
を備えたことを特徴とする容器壁の厚み推定装置。
【請求項３】
内壁面と外壁面とに温度差を有する容器における容器壁の厚み推定演算をコンピュータに実行させるコンピュータプログラムであって、
前記外壁面における温度ｈ(ｔ)及び熱流束を熱伝導率で除した物理量ｇ(ｔ)を取得する処理と、
容器壁の厚みの仮定値ｌ〜を設定する処理と、
内壁面（ｘ＝０）から外壁面方向をｘ軸とし、式（１０１）における容器壁の温度ｕ(ｘ，ｔ)の代替として変数ｖ(ｘ，ｔ)を定義して導入した式（１０２）から、変数ｖ(ｌ〜，ｔi)を算出する処理と、
【数７】
2008063593.tif 000023
変数ｗ(ｘ，ｔ)を定義して導入した式（１０３）から、変数ｗ(ｘ，ｔ)を式（１０４）により表わし、式（１０４）に基づき式（１０５）のＭＡ×ＶＢ＝Ｖｂを解くことによりＶＢを決定する処理と、
【数８】
2008063593.tif 000024
式（１０６）で表わされるｐ(ｌ〜，ｔ)が設定値以下となる仮定値ｌ〜を容器壁の厚みとして決定する処理と、
【数９】
2008063593.tif 000025
をコンピュータに実行させることを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内壁面と外壁面とに温度差を有する容器における容器壁の厚み推定方法、装置、コンピュータプログラムに関する。
【背景技術】
【０００２】
高炉、転炉、脱ガス炉、燃焼による鋼材加熱炉、石炭ガス化反応炉等の高温のガス反応又は液体反応を伴う反応容器や、混銑車、溶銑鍋、溶鋼鍋等の溶鉄を運搬する容器を使用しての操業を管理する場合、これら高温物質を取り扱う容器の壁の状況（例えば、内壁面の損耗状態）を観測し、管理する必要がある。
【０００３】
従来、容器壁の損耗状態は、容器内に溶鋼等の高温物質が存在しないときに、内壁面の状態を目視で観察することで管理されてきた。
【０００４】
しかしながら、容器内に高温物質が存在しないときでも、高温物質の排出直後には耐火物表面は５００℃以上の高温に熱せられている。上記のような目視による推定では、損耗状態を定量的な数値として捉えることは極めて困難であり、定性的な管理とならざるを得ない。また、容器内に高温物質が存在しないことを条件とした管理を余儀無くされるため、稼動中の高温物質の流出という事態を管理することができなかった。
【０００５】
ここで、容器壁の損耗状態は、容器壁の厚みによって判断することができる。例えば、損耗が均一な形状で１次元形状に近似できる場合、容器壁が熱的に定常状態にあれば、容器壁内部の温度分布は直線状になる。したがって、容器壁の厚みＬは、容器の外壁面で計測した熱流束Ｑ、容器壁の厚み方向の熱伝導率ｋx、容器の内壁面温度Ｔin、及び容器の外壁面温度Ｔoutを用いて次式（５１）より推定できる。
【０００６】
【数１】
2008063593.tif 000003
【０００７】
しかしながら、実際の容器壁の温度は、稼動（容器内に高温物質が存在する）・非稼動（容器内に高温物質が存在しない）の時間サイクルによって異なった値を示すため、容器の外壁面で計測した熱流束Ｑも非定常的に変化する。これに加え、容器壁内部の温度分布は曲線形状で非定常的に変化するため、上式で容器壁の厚みを推定すると、大きな誤差を引き起こすことになる。
【０００８】
一方、容器壁内部の熱伝導現象を非定常熱伝導逆問題と考えて、容器壁に設置した温度計測手段によって計測された温度データを基に、非定常熱伝導逆問題により容器壁内部の温度を計算し、容器壁の温度が溶鉄の凝固温度に一致する位置を検索することにより容器壁の厚みを推定する方法が提案されている（特許文献１を参照）。
【０００９】
【特許文献１】特開２００１−２３４２１７号公報
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、特許文献１に開示されている逆問題解析は、容器壁内部の温度を計算した結果から間接的に容器壁の厚みを推定するものである。そのため、特に容器内への高温物質の装入直後は温度計算結果が安定せず、容器壁の厚み推定精度の低下が懸念される。
【００１１】
本発明は上記のような点に鑑みてなされたものであり、従来法のように容器壁内部の温度計算を行わなくとも、容器内に高温物質が存在する状態で、容器壁の厚みを直接的にかつ精度よく推定できるようにすることを目的とする。
【課題を解決するための手段】
【００１２】
本発明による容器壁の厚み推定手法は、内壁面と外壁面とに温度差を有する容器における容器壁の厚みを推定するものであって、
前記外壁面における温度ｈ(ｔ)及び熱流束を熱伝導率で除した物理量ｇ(ｔ)を取得し、
容器壁の厚みの仮定値ｌ〜を設定し、
内壁面（ｘ＝０）から外壁面方向をｘ軸とし、式（１０１）における容器壁の温度ｕ(ｘ，ｔ)の代替として変数ｖ(ｘ，ｔ)を定義して導入した式（１０２）から、変数ｖ(ｌ〜，ｔi)を算出し、
変数ｗ(ｘ，ｔ)を定義して導入した式（１０３）から、変数ｗ(ｘ，ｔ)を式（１０４）により表わし、式（１０４）に基づき式（１０５）のＭＡ×ＶＢ＝Ｖｂを解くことによりＶＢを決定し、
式（１０６）で表わされるｐ(ｌ〜，ｔ)が設定値以下となる仮定値ｌ〜を容器壁の厚みとして決定するようにした。
【００１３】
【数２】
2008063593.tif 000004
【００１４】
【数３】
2008063593.tif 000005
【００１５】
【数４】
2008063593.tif 000006
【００１６】
ここで、内壁面と外壁面とに温度差を有する容器としては、高炉、転炉、脱ガス炉、燃焼による鋼材加熱炉、石炭ガス化反応炉等の高温のガス反応又は液体反応を伴う反応容器や、混銑車、溶銑鍋、溶鋼鍋等の溶鉄を運搬する容器等がある。
【発明の効果】
【００１７】
本発明によれば、従来法のように容器壁内部の温度計算を行わなくとも、容器内に高温物質が存在する状態で、容器壁の厚みを直接的にかつ精度よく推定することができる。
【発明を実施するための最良の形態】
【００１８】
以下、添付図面を参照して、本発明の好適な実施形態について説明する。
まず、本発明による容器壁の厚み推定手法の基本的な考え方について説明する。図１は、容器壁の一部を表わす図であり、ｘ＝０が容器の内壁面の位置である。同図において、容器壁の残存厚みｌ、容器内に存在する高温物質の温度ｆ(ｔ)＝ＵM、容器壁の温度ｕ(ｘ，ｔ)（外壁面の温度計測点にて計測された温度ｈ(ｔ)）、温度ｈ(ｔ)を基に算出した熱流束（又は外壁面の温度計測点にて計測された熱流束）を熱伝導率で除した物理量ｇ(ｔ)である。
【００１９】
（定式化）
式（１）は、非定常熱伝導方程式を表わす。なお。ｕtは∂ｕ／∂ｔを、ｕxxは∂2ｕ／∂ｘ2を表わす。式（１）において、αは熱拡散係数、ｕ(ｘ，０)＝ｕ0(ｘ)は容器壁の温度の初期値である。この場合、容器壁の温度の初期値ｕ(ｘ，０)＝ｕ0(ｘ)は未知である。
【００２０】
また、高温物質の温度ＵM、外気温度ｕa、熱拡散係数α、放射伝熱のステファンボルツマン係数σ、容器壁の熱伝導率λは正の定数である。
【００２１】
【数５】
2008063593.tif 000007
【００２２】
ここで、式（３）を導入してフーリエ展開することにより、式（１）の容器壁の温度ｕ(ｘ，ｔ)は式（４）のように求められる。
【００２３】
【数６】
2008063593.tif 000008
【００２４】
【数７】
2008063593.tif 000009
【００２５】
ｘ＝ｌとすると、式（５）が得られる。
【００２６】
【数８】
2008063593.tif 000010
【００２７】
ところが、コンピュータによる演算処理を実行する場合、式（５）の右辺において、特に第２項、第３項の計算は打ち切り誤差を引き起こしやすいという問題がある。
【００２８】
そこで、本発明においては、容器壁の温度ｕ(ｘ，ｔ)の代替として変数ｖ(ｘ，ｔ)を定義し、式（６）を導入する。式（６）において、変数の初期値ｖ(ｘ，０)＝ｘｇ(０)＋ｆ(０)は、内壁面の熱流束を熱伝導率で除した物理量の初期値ｇ(０)、及び、内壁面の温度の初期値（即ち、高温物質の温度）ｆ(０)をいずれも既知とできるので、既知とすることができる。したがって、変数ｖ(ｘ，ｔ)については、例えば後退差分法により直接計算することができる。
【００２９】
【数９】
2008063593.tif 000011
【００３０】
さらに、容器壁の温度ｕ(ｘ，ｔ)と変数ｖ(ｘ，ｔ)との差をｗ(ｘ，ｔ)と定義すると、式（７）のようになる。
【００３１】
【数１０】
2008063593.tif 000012
【００３２】
ここで、上述したのと同様にフーリエ展開することにより、式（７）の変数ｗ(ｘ，ｔ)は式（８）のように求められる。式（８）においては、式（４）と比較して明らかなように、第２項、第３項のない簡単な式とすることができる。
【００３３】
【数１１】
2008063593.tif 000013
【００３４】
また、ｗ(ｌ，ｔ)＝ｕ(ｌ，ｔ)−ｖ(ｌ，ｔ)である。そして、ｕ(ｌ，ｔ)は既知のｈ(ｔ)であり、また、ｖ(ｌ，ｔ)は式（６）から後退差分法により直接計算することができる。したがって、ｗ(ｌ，ｔ)が既知であるとして、式（８）からＢn(ｌ)の近似値を得ることができる。
【００３５】
（残存厚みｌを求めるための逆問題）
逆問題においては、容器壁の残存厚みｌは未知であり、したがって変数ｖ(ｌ，ｔ)は未知であるが、ｕ(ｌ，ｔ)は計測値ｈ(ｔ)として与えられる。式（８）から式（９）が得られる。
【００３６】
【数１２】
2008063593.tif 000014
【００３７】
以下述べるように、最適化計算により、容器壁の残存厚みｌの近似値を得ることができる。即ち、観察時間を（Ｔst，Ｔend）と設定し、Ｔst＜Ｔ1＜Ｔ2＜Ｔendとする。そして、Ｔ1＝ｔ1＜ｔ2＜・・・＜ｔM＝Ｔ2と均一格子にする。残存厚みの仮定値ｌ〜＞０は既知条件とする。なお、本明細書において、ｌ〜の表記は、ｌの上に〜が付されているものとする。
【００３８】
【数１３】
2008063593.tif 000015
【００３９】
ここで、外壁面の温度計測点にて計測された温度ｈ(ｔ)は既知で、変数ｖ(ｌ〜，ｔi)は仮定値ｌ〜を与えることにより式（６）から後退差分法により求められる。式（１０）のように、ＭＡはＭ×（Ｎ＋１）行列、ＶＢは（Ｎ＋１）×１ベクトル、ＶｂはＭ×１ベクトルであって、ＭＡ×ＶＢ＝Ｖｂを解くことにより、Ｂ0(ｌ〜)、Ｂ1(ｌ〜)、・・・、ＢN(ｌ〜)が求められる。
【００４０】
そして、実測によるｗ(ｌ，ｔ)と計算によるｗ(ｌ〜，ｔ)との差分を表わす式（１１）を定義する。ｔ∈（Ｔ2，Ｔend）としてＢi(ｌ〜)を式（１１）に代入すると、ｐ(ｌ〜，ｔ)が得られるので、ｐ(ｌ〜，ｔ)が０に近づくように残存厚みの仮定値ｌ〜を選択することにより、その仮定値ｌ〜を容器壁の残存厚みｌの近似値として求めることができる。
【００４１】
【数１４】
2008063593.tif 000016
【００４２】
図２は、本実施形態に係る容器壁の厚み推定装置の概略構成を示す図である。また、図３は、本実施形態における容器壁の厚み推定処理を説明するためのフローチャートである。
【００４３】
図２において、１０１は入力部であり、容器の外壁面の温度計測点にて計測された温度ｈ(ｔi)が入力される（ステップＳ１０１）。
【００４４】
１０２は熱流束算出部であり、入力部１０１に入力される温度ｈ(ｔi)を基に熱流束を算出し、熱伝導率で除した物理量ｇ(ｔi)を算出する（ステップＳ１０２）。なお、外壁面の温度計測点にて計測された熱流束が入力部１０１に入力されるようにしてもよく、その場合、熱流束算出部１０２は不要である。
【００４５】
１０３は演算部であり、上述した容器壁の厚み推定手法により容器壁の厚みを演算する。即ち、まず残存厚みの仮定値ｌ〜を設定する（ステップＳ１０３）。次に、仮定値ｌ〜を与えることにより、式（６）から変数ｖ(ｌ〜，ｔi)を求める（ステップＳ１０４）。そして、式（１０）のＭＡ×ＶＢ＝Ｖｂを解くことにより、Ｂ0(ｌ〜)、Ｂ1(ｌ〜)、・・・、ＢN(ｌ〜)を求め（ステップＳ１０５）、式（１１）のｐ(ｌ〜，ｔ)を算出する（ステップＳ１０６）。設定値ε以下となるｐ(ｌ〜，ｔ)が得られるまでステップＳ１０３〜Ｓ１０６を繰り返し（ステップＳ１０７）、設定値ε以下となったときの仮定値ｌ〜を容器壁の残存厚みｌとして決定する（ステップＳ１０８）。
【００４６】
１０４は出力部であり、演算部１０３により演算、推定された容器壁の残存厚みを、例えば不図示のディスプレイに表示等する。
【００４７】
（実施例）
本発明の手法（本法）による容器壁の厚みの推定結果について説明する。予め厚みが分かっている３種の溶鋼鍋（容器壁の厚み１００ｍｍ、１５０ｍｍ、２００ｍｍ）について、本法と、従来法（特許文献１に開示された手法）とにより容器壁の厚みを推定し、その結果を比較した。
【００４８】
本実施例において、溶鋼鍋（耐火物）の熱拡散係数αは０．００８６５ｍ2/Ｈｒ（熱伝導率：８．４９Ｗ/ｍ/Ｋ、比熱：１２１４．２Ｊ/ｋｇ/Ｋ、密度２９１０ｋｇ/ｍ3）であった。
【００４９】
また、図４〜６に、各溶鋼鍋の外壁面の温度及び熱流束のデータを示す。本実施例では、受鋼から１２００秒経過後、外壁面の温度計測点にて５秒おきに温度を５回計測し、各温度を基に熱流束を算出した。
【００５０】
図７に示すように、従来法では、特に容器壁が厚い溶鋼鍋において、残存厚みの推定値に誤差が生じやすくなっている。既述したように、従来法では、容器壁内部の温度を計算した結果から間接的に容器壁の厚みを推定するため、この誤差が生じたものと考えられる。
【００５１】
それに対して、本法では、いずれの溶鋼鍋においても、実際の厚みと推定厚みとが略一致しており、良好な推定値が得られた。
【００５２】
以上述べた本発明の目的は、上述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、システム或いは装置に供給し、そのシステム或いは装置のコンピュータ（又はＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。
【００５３】
この場合、記憶媒体から読み出されたプログラムコード自体が上述した実施形態の機能を実現することになり、プログラムコード自体及びそのプログラムコードを記憶した記憶媒体は本発明を構成することになる。プログラムコードを供給するための記憶媒体としては、例えば、フレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【図面の簡単な説明】
【００５４】
【図１】容器壁の一部を表わす図である。
【図２】本実施形態に係る容器壁の厚み推定装置の概略構成を示す図である。
【図３】本実施形態における容器壁の厚み推定処理を説明するためのフローチャートである。
【図４】容器壁の厚みが１００ｍｍの場合における外壁面の温度及び熱流束のデータを示す特性図である。
【図５】容器壁の厚みが１５０ｍｍの場合における外壁面の温度及び熱流束のデータを示す特性図である。
【図６】容器壁の厚みが２００ｍｍの場合における外壁面の温度及び熱流束のデータを示す特性図である。
【図７】本法と従来法とにより容器壁の厚みを推定した結果を比較する特性図である。
【符号の説明】
【００５５】
１０１  入力部
１０２  熱流束算出部
１０３  演算部
１０４  出力部
【図１】
2008063593.tif 000026
【図２】
2008063593.tif 000027
【図３】
2008063593.tif 000028
【図４】
2008063593.tif 000029
【図５】
2008063593.tif 000030
【図６】
2008063593.tif 000031
【図７】
2008063593.tif 000032