JPH06117989A - 粒子分布ピッチ測定方法および装置 - Google Patents
粒子分布ピッチ測定方法および装置Info
- Publication number
- JPH06117989A JPH06117989A JP4265127A JP26512792A JPH06117989A JP H06117989 A JPH06117989 A JP H06117989A JP 4265127 A JP4265127 A JP 4265127A JP 26512792 A JP26512792 A JP 26512792A JP H06117989 A JPH06117989 A JP H06117989A
- Authority
- JP
- Japan
- Prior art keywords
- particle
- light
- optical fiber
- calculated
- data
- 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.)
- Pending
Links
Landscapes
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
(57)【要約】
【目的】 高炉内装入物の移動層内の粒子分布ピッチを
簡易かつ迅速に測定する。 【構成】 投光用ファイバおよび受光用ファイバからな
り、検出視野Dが装入物粒子直径より小さい光ケ−ブル
2個38,39を、粒子直径程度の距離ΔL離して装備
したプロ−ブ3を、該ΔLの方向を垂直方向として高炉
に差し込み、上流側光ケ−ブル38の受光信号に対する
下流側光ケ−ブル39の受光信号の位相差から装入物降
下速度Vmを算出し、一方上流側光ケ−ブル38のみの
受光信号の周期Tphmを算出し、降下速度Vmと周期Tp
hmより装入物粒子分布ピッチPthmを算出する。
簡易かつ迅速に測定する。 【構成】 投光用ファイバおよび受光用ファイバからな
り、検出視野Dが装入物粒子直径より小さい光ケ−ブル
2個38,39を、粒子直径程度の距離ΔL離して装備
したプロ−ブ3を、該ΔLの方向を垂直方向として高炉
に差し込み、上流側光ケ−ブル38の受光信号に対する
下流側光ケ−ブル39の受光信号の位相差から装入物降
下速度Vmを算出し、一方上流側光ケ−ブル38のみの
受光信号の周期Tphmを算出し、降下速度Vmと周期Tp
hmより装入物粒子分布ピッチPthmを算出する。
Description
【0001】
【産業上の利用分野】本発明は、粒子群中の粒子分布の
測定に関し、例えば高炉内に炉頂部から装入されたコ−
クス,鉱石などにより形成される移動層内の粒子分布ピ
ッチを測定する方法および装置に関する。
測定に関し、例えば高炉内に炉頂部から装入されたコ−
クス,鉱石などにより形成される移動層内の粒子分布ピ
ッチを測定する方法および装置に関する。
【0002】
【従来技術】通常、移動層内の粒子径又は粒子分布を求
めるには、ファイバスコープなどにより移動層内を直接
観察し、その画像を画像処理装置によって処理して求め
る方法が知られている。ファイバスコープは多数の光フ
ァイバの素線を整然と精密に配置して束ねたものであ
る。
めるには、ファイバスコープなどにより移動層内を直接
観察し、その画像を画像処理装置によって処理して求め
る方法が知られている。ファイバスコープは多数の光フ
ァイバの素線を整然と精密に配置して束ねたものであ
る。
【0003】
【発明が解決しようとする課題】しかし、従来技術によ
れば以下に示すような問題を生じる。
れば以下に示すような問題を生じる。
【0004】(1) 計測系がファイバスコープおよび2次
元画像処理装置によって構成されるため全体的にコスト
高となる。
元画像処理装置によって構成されるため全体的にコスト
高となる。
【0005】(2) 2次元画像処理の演算に時間がかか
り、リアルタイムでの粒子径又は粒子分布の変化の検出
に使用することは困難である。
り、リアルタイムでの粒子径又は粒子分布の変化の検出
に使用することは困難である。
【0006】(3) 粒子の輪郭を求める画像処理ソフトは
複雑であり、特に粒子が重なり合った場合、重なり部分
を除外することは困難であり精度上も問題がある。ま
た、高炉内装入物の場合、装入物(コ-クス,鉱石)の形状
がランダムであり表面凹凸も大きいので、粒子輪郭を正
確に摘出しにくい。
複雑であり、特に粒子が重なり合った場合、重なり部分
を除外することは困難であり精度上も問題がある。ま
た、高炉内装入物の場合、装入物(コ-クス,鉱石)の形状
がランダムであり表面凹凸も大きいので、粒子輪郭を正
確に摘出しにくい。
【0007】(4) ファイバスコ-プの視野が、粒子径又
は粒子分布ピッチの数倍である必要があり、検出部端部
が大型になる。
は粒子分布ピッチの数倍である必要があり、検出部端部
が大型になる。
【0008】そこで本発明は、上記問題を解決して粒子
分布ピッチを簡易かつ迅速に測定する装置を提供するこ
とを目的とする。
分布ピッチを簡易かつ迅速に測定する装置を提供するこ
とを目的とする。
【0009】
【課題を解決するための手段】本発明では、粒子群への
光の投射と光ファイバ(38)を介した粒径代表値(数10mm)
よりも小さい径(D=10mm)の視野での粒子群の反射光の
検出を、粒子群と光ファイバ(38)が相対的に移動する状
態で、連続もしくは所定周期(Ts)で行なって反射光強度
を表わす受光信号Iaを得、受光信号Iaのレベル変動
周期T(Tphm)を算出し、該周期T(Tphm)と前記相対移動
の速度Vより、粒子群中の粒子分布ピッチT・V(Tphm・
Vm)を算出する。
光の投射と光ファイバ(38)を介した粒径代表値(数10mm)
よりも小さい径(D=10mm)の視野での粒子群の反射光の
検出を、粒子群と光ファイバ(38)が相対的に移動する状
態で、連続もしくは所定周期(Ts)で行なって反射光強度
を表わす受光信号Iaを得、受光信号Iaのレベル変動
周期T(Tphm)を算出し、該周期T(Tphm)と前記相対移動
の速度Vより、粒子群中の粒子分布ピッチT・V(Tphm・
Vm)を算出する。
【0010】なお、カッコ内の記号は、図面に示し後述
する実施例装置の対応要素又は対応事項を示す。
する実施例装置の対応要素又は対応事項を示す。
【0011】本発明の好ましい実施例では、次のように
して前記相対移動の速度V(Vm)を測定する。すなわち、
前記光の反射点から、前記視野径(D)より大きく粒径代
表値の数十倍以内の所定距離ΔL(20mm)離れた位置への
光の投射と光ファイバ(39)を介した(数10mm)よりも小さ
い径(D)の視野での粒子群の反射光の検出を連続もしく
は所定周期(Ts)で行なって反射光強度を表わす第2受光
信号Ibを得、受光信号Iaに対する第2受光信号Ib
の位相差Tdを算出し、位相差Tdと距離ΔLより速度
V=ΔL/Tdを算出する。
して前記相対移動の速度V(Vm)を測定する。すなわち、
前記光の反射点から、前記視野径(D)より大きく粒径代
表値の数十倍以内の所定距離ΔL(20mm)離れた位置への
光の投射と光ファイバ(39)を介した(数10mm)よりも小さ
い径(D)の視野での粒子群の反射光の検出を連続もしく
は所定周期(Ts)で行なって反射光強度を表わす第2受光
信号Ibを得、受光信号Iaに対する第2受光信号Ib
の位相差Tdを算出し、位相差Tdと距離ΔLより速度
V=ΔL/Tdを算出する。
【0012】
【作用】光ファイバ(38)を介して、粒径代表値(数10mm)
よりも小さい径(D=10mm)の視野で、粒子群の反射光を
検出し反射光強度を表わす受光信号Iaを得るので、受
光信号Iaは、粒子群と光ファイバ(38)の相対移動速度
Vに反比例し粒子分布ピッチに比例する周期T(Tphm)の
レベル変化を示す。相対移動速度Vが概知であると、又
は測定により相対移動速度Vを得ると、粒子分布ピッチ
T・V(Tphm・Vm)が得られる。
よりも小さい径(D=10mm)の視野で、粒子群の反射光を
検出し反射光強度を表わす受光信号Iaを得るので、受
光信号Iaは、粒子群と光ファイバ(38)の相対移動速度
Vに反比例し粒子分布ピッチに比例する周期T(Tphm)の
レベル変化を示す。相対移動速度Vが概知であると、又
は測定により相対移動速度Vを得ると、粒子分布ピッチ
T・V(Tphm・Vm)が得られる。
【0013】よって、2次元画像処理によって粒子輪郭
を切出すような複雑な処理は必要なく高速演算処理が可
能である。この粒子分布ピッチT・V(Tphm・Vm)に補正
を施すことにより粒子径(Dpm)を算出することができ
る。
を切出すような複雑な処理は必要なく高速演算処理が可
能である。この粒子分布ピッチT・V(Tphm・Vm)に補正
を施すことにより粒子径(Dpm)を算出することができ
る。
【0014】前述の好ましい実施例によれば、第1組の
光ファイバ(39)の受光信号Iaは、粒子分布ピッチの算
出用および相対速度Vの算出用に兼用される。すなわ
ち、受光信号Iaを得る位置から、視野径(D)より大き
く粒径代表値の数十倍以内の所定距離ΔL(20mm)離れた
位置の反射光を、光ファイバ(39)を介して粒径代表値
(数10mm)よりも小さい径(D)の視野で検出し反射光強度
を表わす第2受光信号Ibを得るので、第2受光信号I
bも、相対移動に連動して粒子分布ピッチに相当するレ
ベル変化を示し、しかも、受光信号Iaに対して、前記
所定距離ΔL(20mm)および相対速度(Vm)に対応した位相
差Td(Tm)すなわち遅れ(Tm)を有するものとなる。よっ
て、位相差Td(Tm)と前記所定距離ΔLより、粒子の相
対移動速度V=ΔL/Td(Vm=ΔL/Tm)が求まる。
所定距離ΔL(20mm)が、視野径(D)より大きく粒径代表
値の数十倍以内であるので、群単位(層単位)ではな
く、この所定距離ΔL単位すなわち粒径代表値の数十倍
以内の単位での相対移動速度測定となり、比較的に短周
期の速度変動を把握しうる。ΔLを粒径代表値の数倍以
内に設定することにより、粒子群中の粒子ピッチもしく
は該ピッチの数倍の群単位で、すなわちきわめて高分解
能で、降下速度を測定しうる。
光ファイバ(39)の受光信号Iaは、粒子分布ピッチの算
出用および相対速度Vの算出用に兼用される。すなわ
ち、受光信号Iaを得る位置から、視野径(D)より大き
く粒径代表値の数十倍以内の所定距離ΔL(20mm)離れた
位置の反射光を、光ファイバ(39)を介して粒径代表値
(数10mm)よりも小さい径(D)の視野で検出し反射光強度
を表わす第2受光信号Ibを得るので、第2受光信号I
bも、相対移動に連動して粒子分布ピッチに相当するレ
ベル変化を示し、しかも、受光信号Iaに対して、前記
所定距離ΔL(20mm)および相対速度(Vm)に対応した位相
差Td(Tm)すなわち遅れ(Tm)を有するものとなる。よっ
て、位相差Td(Tm)と前記所定距離ΔLより、粒子の相
対移動速度V=ΔL/Td(Vm=ΔL/Tm)が求まる。
所定距離ΔL(20mm)が、視野径(D)より大きく粒径代表
値の数十倍以内であるので、群単位(層単位)ではな
く、この所定距離ΔL単位すなわち粒径代表値の数十倍
以内の単位での相対移動速度測定となり、比較的に短周
期の速度変動を把握しうる。ΔLを粒径代表値の数倍以
内に設定することにより、粒子群中の粒子ピッチもしく
は該ピッチの数倍の群単位で、すなわちきわめて高分解
能で、降下速度を測定しうる。
【0015】本発明の他の目的および特徴は図面を参照
した以下の実施例の説明より明らかになろう。
した以下の実施例の説明より明らかになろう。
【0016】
【実施例】図1に、本発明装置の一実施例を示す。この
実施例は、高炉内装入物(コ−クス,鉱石等)の各移動
層(コ−クス層,鉱石層)内の粒子分布ピッチを測定し
この粒子分布ピッチより粒子径を算出するものである。
この測定装置は、装置全体の制御を行うマイクロコンピ
ュータ(CPU)1,演算処理用のマイクロコンピュー
タ(CPU)2,高炉内に直接挿入されるプローブ3,
このプロ−ブ3を介して高炉装入物に投光するレーザ光
源(He−Neレーザ)4,プロ−ブ3を介して反射光
を受光し電気信号に変換して出力する第1受光素子(上
流側)5および第2受光素子(下流側)6,レーザ光源
4を発光付勢するレーザドライバ7,受光素子5および
6の光電変換信号にノイズ除去,増幅および校正の処理
を施こす増幅器8および9,インターフェイス10,操
作表示ボード11,RAM12,およびROM13など
から構成されている。
実施例は、高炉内装入物(コ−クス,鉱石等)の各移動
層(コ−クス層,鉱石層)内の粒子分布ピッチを測定し
この粒子分布ピッチより粒子径を算出するものである。
この測定装置は、装置全体の制御を行うマイクロコンピ
ュータ(CPU)1,演算処理用のマイクロコンピュー
タ(CPU)2,高炉内に直接挿入されるプローブ3,
このプロ−ブ3を介して高炉装入物に投光するレーザ光
源(He−Neレーザ)4,プロ−ブ3を介して反射光
を受光し電気信号に変換して出力する第1受光素子(上
流側)5および第2受光素子(下流側)6,レーザ光源
4を発光付勢するレーザドライバ7,受光素子5および
6の光電変換信号にノイズ除去,増幅および校正の処理
を施こす増幅器8および9,インターフェイス10,操
作表示ボード11,RAM12,およびROM13など
から構成されている。
【0017】図2に、プローブ3の拡大正面(先端面)
を、図3にプロ−ブ3の拡大側断面を、また図4に図3
中に示す光ファイバ38,39の拡大横断面を示す。
を、図3にプロ−ブ3の拡大側断面を、また図4に図3
中に示す光ファイバ38,39の拡大横断面を示す。
【0018】図2に示すように、プロ−ブ3の先端面に
は2つの窓開口があり、その内側に、図3に示すよう
に、保護ガラス35がある。外筒31は2重管構造の中
空筒体であって、内壁板と外壁板の間に筒状の仕切板3
2があり、この仕切板32で区分された内側空間33に
プロ−ブ後端(図示せず)から冷却ガスが供給される。
冷却ガスは内側空間33を通ってプロ−ブ先端に達しそ
こで折り返して外側空間34に入り外側空間34を通っ
てプロ−ブ後端に向かう。この冷却ガス流により、プロ
−ブ3の内部は、炉内挿入時に300℃以下の温度に保
持される。この冷却により、レンズ36,37および光
ファイバケ−ブル38,39の光学特性が、実質上一定
に維持される。外筒31の内壁板の内側の空間に第1光
ファイバケ−ブル38と第2光ファイバケ−ブル39が
挿入されており、それらの先端がレンズ36,37に対
向している。
は2つの窓開口があり、その内側に、図3に示すよう
に、保護ガラス35がある。外筒31は2重管構造の中
空筒体であって、内壁板と外壁板の間に筒状の仕切板3
2があり、この仕切板32で区分された内側空間33に
プロ−ブ後端(図示せず)から冷却ガスが供給される。
冷却ガスは内側空間33を通ってプロ−ブ先端に達しそ
こで折り返して外側空間34に入り外側空間34を通っ
てプロ−ブ後端に向かう。この冷却ガス流により、プロ
−ブ3の内部は、炉内挿入時に300℃以下の温度に保
持される。この冷却により、レンズ36,37および光
ファイバケ−ブル38,39の光学特性が、実質上一定
に維持される。外筒31の内壁板の内側の空間に第1光
ファイバケ−ブル38と第2光ファイバケ−ブル39が
挿入されており、それらの先端がレンズ36,37に対
向している。
【0019】図4に示すように、光ファイバケ−ブル3
8,39の光ファイバは、光投射用のもの38e,39
eと受光用のもの38r,39rに区分されており、図
1に示すように、両ケ−ブルの光投射用の光ファイバ3
8e,39eはレ−ザ光源4に結合されており、光源4
のレ−ザ光をレンズ36,37に案内する。第1光ファ
イバケ−ブル38の受光用の光ファイバ38rは第1受
光素子5に結合され、第2光ファイバケ−ブル39の受
光用の光ファイバ39rは第2受光素子6に結合されて
いる。
8,39の光ファイバは、光投射用のもの38e,39
eと受光用のもの38r,39rに区分されており、図
1に示すように、両ケ−ブルの光投射用の光ファイバ3
8e,39eはレ−ザ光源4に結合されており、光源4
のレ−ザ光をレンズ36,37に案内する。第1光ファ
イバケ−ブル38の受光用の光ファイバ38rは第1受
光素子5に結合され、第2光ファイバケ−ブル39の受
光用の光ファイバ39rは第2受光素子6に結合されて
いる。
【0020】プロ−ブ3は、その先端(保護ガラス35
側)が、高炉の炉壁を通して高炉内の装入物に差し込ま
れ、しかも、レンズ36がレンズ37より上方でそれら
の中心を結ぶ線が垂直になるように、姿勢が設定され
る。この状態でレ−ザ光源4が発したレ−ザ光は、光投
射用の光ファイバ38e,39eを介して高炉内の装入
物に投射され、装入物による反射光が受光用の光ファイ
バ38r,39rを介して受光素子5,6で検出され
る。
側)が、高炉の炉壁を通して高炉内の装入物に差し込ま
れ、しかも、レンズ36がレンズ37より上方でそれら
の中心を結ぶ線が垂直になるように、姿勢が設定され
る。この状態でレ−ザ光源4が発したレ−ザ光は、光投
射用の光ファイバ38e,39eを介して高炉内の装入
物に投射され、装入物による反射光が受光用の光ファイ
バ38r,39rを介して受光素子5,6で検出され
る。
【0021】この実施例では、高炉内の装入物の粒子直
径が数十mmであり代表値が30mm前後であるので、プロ
ーブ3の先端における光ファイバケ−ブル38,39の
受光視野直径Dは、粒子分布ピッチの分解能がある受光
信号を得るために、粒子直径代表値より小さい10mmに
している。上下の受光の干渉を防止するために、受光視
野(36,37)間の最短距離ΔLminは視野直径D以
上にするのが好ましくこの実施例ではΔLmin=10mm
に設定している。また、光ファイバケ−ブル38,39
の受光視野(36,37)の垂直方向距離ΔLは、装入
物粒子分布ピッチ相当の装入物区分で降下速度を測定す
るために、20mmにしている。このように小型(Dが小
&ΔLが短))に、しかも一個のプローブ3で上下2点
の光反射測定をする構成としたので、コークス,鉱石各
層(数百mm)内の粒子個々の動きを検知することができ
る。
径が数十mmであり代表値が30mm前後であるので、プロ
ーブ3の先端における光ファイバケ−ブル38,39の
受光視野直径Dは、粒子分布ピッチの分解能がある受光
信号を得るために、粒子直径代表値より小さい10mmに
している。上下の受光の干渉を防止するために、受光視
野(36,37)間の最短距離ΔLminは視野直径D以
上にするのが好ましくこの実施例ではΔLmin=10mm
に設定している。また、光ファイバケ−ブル38,39
の受光視野(36,37)の垂直方向距離ΔLは、装入
物粒子分布ピッチ相当の装入物区分で降下速度を測定す
るために、20mmにしている。このように小型(Dが小
&ΔLが短))に、しかも一個のプローブ3で上下2点
の光反射測定をする構成としたので、コークス,鉱石各
層(数百mm)内の粒子個々の動きを検知することができ
る。
【0022】なお、図4に模式的に示すように光ファイ
バケ−ブル38,39内では、投光用光ファイバ38
e,39eおよび受光用光ファイバ38r,39rがラ
ンダムに分布している。一本のファイバの直径は125μm
である。この集合構成により、小視野(D)内の反射光
を効率良く収集しうる。なお、これら投,受光用ファイ
バの配置はこれに限らず、例えば投光用光ファイバを受
光用光ファイバ群の横や周囲などに配置するようにして
もよい。投光用光ファイバグル−プと受光用光ファイバ
グル−プとを塊で分離した場合、装入物の光反射方向が
受光用光ファイバグル−プの受光端から外れると受光が
無くなるなど受光効率が低下するので、それを防止する
工夫、例えば投光,受光ファイバ各々にレンズを設け、
その焦点を一致させるなどの工夫を施した方が好まし
い。
バケ−ブル38,39内では、投光用光ファイバ38
e,39eおよび受光用光ファイバ38r,39rがラ
ンダムに分布している。一本のファイバの直径は125μm
である。この集合構成により、小視野(D)内の反射光
を効率良く収集しうる。なお、これら投,受光用ファイ
バの配置はこれに限らず、例えば投光用光ファイバを受
光用光ファイバ群の横や周囲などに配置するようにして
もよい。投光用光ファイバグル−プと受光用光ファイバ
グル−プとを塊で分離した場合、装入物の光反射方向が
受光用光ファイバグル−プの受光端から外れると受光が
無くなるなど受光効率が低下するので、それを防止する
工夫、例えば投光,受光ファイバ各々にレンズを設け、
その焦点を一致させるなどの工夫を施した方が好まし
い。
【0023】図5に、図1に示すCPU1の処理内容を
示す。CPU1は、操作・表示ボ−ド11より、粒子分
布ピッチ測定が指定されると図5に示す処理に進む。こ
の処理の最初には、粒子分布ピッチ測定に関連するフラ
グやレジスタ等をイニシャライズする(ステップ1:以
下カッコ内ではステップという語を省略する)。次に、
操作・表示ボ−ド11の指示に応答して増幅器8,9の
増幅率を設定しオフセット調整を行う(2)。なお、オ
ペレ−タは、プロ−ブ3の先端面全面に基準反射板を当
てて、調整指示を入力する。これに応答してCPU1が
増幅器8の出力信号レベルと増幅器9の出力信号レベル
の差が許容範囲内になるように、増幅器8,9の増幅率
等を調整する(2)。
示す。CPU1は、操作・表示ボ−ド11より、粒子分
布ピッチ測定が指定されると図5に示す処理に進む。こ
の処理の最初には、粒子分布ピッチ測定に関連するフラ
グやレジスタ等をイニシャライズする(ステップ1:以
下カッコ内ではステップという語を省略する)。次に、
操作・表示ボ−ド11の指示に応答して増幅器8,9の
増幅率を設定しオフセット調整を行う(2)。なお、オ
ペレ−タは、プロ−ブ3の先端面全面に基準反射板を当
てて、調整指示を入力する。これに応答してCPU1が
増幅器8の出力信号レベルと増幅器9の出力信号レベル
の差が許容範囲内になるように、増幅器8,9の増幅率
等を調整する(2)。
【0024】オペレ−タは、調整を終えるとプロ−ブ3
を高炉炉壁の穴に通してプロ−ブ3の先端を炉内装入物
中に所要深さまで差し込み、プロ−ブ3の姿勢を、レン
ズ36が37の上方でしかもそれらの中心を結ぶ線が垂
直となる、予め定められた姿勢に設定する。そして操作
・表示ボ−ド11でCPU1にスタ−トを指示する。 −1024×3回分のデ−タ蓄積− このスタ−ト指示に応答してCPU1が図5のステップ
3以下の処理を開始する。すなわち、サンプリング回数
をカウントするためのレジスタmに数値1を示すデ−タ
をセットし(3)、タイマTsをスタートする(4)。
タイマ時限Tsは本実施例では 400/1024 sec であり、
これがサンプリング周期となる。CPU1は次に、増幅
器8,9の出力信号をデジタル変換して読込む。すなわ
ち、上流側の光ファイバ38rを介して第1受光素子
(上流側)5が受光した光強度(受光信号)Iaと、下
流側の光ファイバ39rを介して第2受光素子(下流
側)6が受光した光強度(受光信号)Ibを読み込む
(5)。次に、受光強度IaおよびIbを記憶するための
バッファメモリMa,Mb(RAM12内に割り当ててい
る)のデータを1アドレス分シフトさせる(6,7)。
メモリMa,Mbはそれぞれ1024×3個のアドレスを有
し、データシフトではアドレス2のデータをアドレス1
に書込み、アドレス3のデータをアドレス2に書込み、
・・・,アドレス1024×3のデータをアドレス1024×3-1
に書込む。そしてこのシフトで開いたアドレス1024×3
に最新の受光強度Ia,Ibをメモリする(8,9)。
を高炉炉壁の穴に通してプロ−ブ3の先端を炉内装入物
中に所要深さまで差し込み、プロ−ブ3の姿勢を、レン
ズ36が37の上方でしかもそれらの中心を結ぶ線が垂
直となる、予め定められた姿勢に設定する。そして操作
・表示ボ−ド11でCPU1にスタ−トを指示する。 −1024×3回分のデ−タ蓄積− このスタ−ト指示に応答してCPU1が図5のステップ
3以下の処理を開始する。すなわち、サンプリング回数
をカウントするためのレジスタmに数値1を示すデ−タ
をセットし(3)、タイマTsをスタートする(4)。
タイマ時限Tsは本実施例では 400/1024 sec であり、
これがサンプリング周期となる。CPU1は次に、増幅
器8,9の出力信号をデジタル変換して読込む。すなわ
ち、上流側の光ファイバ38rを介して第1受光素子
(上流側)5が受光した光強度(受光信号)Iaと、下
流側の光ファイバ39rを介して第2受光素子(下流
側)6が受光した光強度(受光信号)Ibを読み込む
(5)。次に、受光強度IaおよびIbを記憶するための
バッファメモリMa,Mb(RAM12内に割り当ててい
る)のデータを1アドレス分シフトさせる(6,7)。
メモリMa,Mbはそれぞれ1024×3個のアドレスを有
し、データシフトではアドレス2のデータをアドレス1
に書込み、アドレス3のデータをアドレス2に書込み、
・・・,アドレス1024×3のデータをアドレス1024×3-1
に書込む。そしてこのシフトで開いたアドレス1024×3
に最新の受光強度Ia,Ibをメモリする(8,9)。
【0025】タイマTsがオーバ(Tsが経過)すると
(10)、CPU1はフラグレジスタMFに、「1」
(1024×3回のデ−タIa,Ibの読込みとメモリM
a,Mbへの書込みが終了した)がセットされているか
否かをチェックし(11)、セットされていないと、レ
ジスタmが1024×3以上であるかチェックし(12)、1
024×3未満であるとレジスタmを1インクレメントし
(13)、レジスタmが1024×3以上となるまでステッ
プ4〜13の処理を繰り返す。これにより1024×3個の
受光強度Ia,IbがメモリMa,Mbに蓄積される。
(10)、CPU1はフラグレジスタMFに、「1」
(1024×3回のデ−タIa,Ibの読込みとメモリM
a,Mbへの書込みが終了した)がセットされているか
否かをチェックし(11)、セットされていないと、レ
ジスタmが1024×3以上であるかチェックし(12)、1
024×3未満であるとレジスタmを1インクレメントし
(13)、レジスタmが1024×3以上となるまでステッ
プ4〜13の処理を繰り返す。これにより1024×3個の
受光強度Ia,IbがメモリMa,Mbに蓄積される。
【0026】−デ−タの更新− ステップ12でレジスタmが1024×3以上となるとフラ
グMFに「1」をセットし、レジスタnに0をセットす
る(14,15)。そして、演算処理用のCPU2に粒
子速度の演算を指示し(16)、ステップ4に戻る。そ
してステップ4〜10でデータを更新する。すなわち最
新のデ−タを読込み(5)、デ−タシフト(6,7)に
よりメモリMa,Mbの最も古いデ−タ(アドレス1の
デ−タ)を捨てて、最新のデ−タをメモリMa,Mbに
書込む(8,9)。このデ−タの更新をTs周期で繰返
す。したがって、メモリMa,Mbには、最近から1024
×3×Ts前までの、1024×3回分のサンプリングデータ
Ia,Ibが常に格納されていることになる。
グMFに「1」をセットし、レジスタnに0をセットす
る(14,15)。そして、演算処理用のCPU2に粒
子速度の演算を指示し(16)、ステップ4に戻る。そ
してステップ4〜10でデータを更新する。すなわち最
新のデ−タを読込み(5)、デ−タシフト(6,7)に
よりメモリMa,Mbの最も古いデ−タ(アドレス1の
デ−タ)を捨てて、最新のデ−タをメモリMa,Mbに
書込む(8,9)。このデ−タの更新をTs周期で繰返
す。したがって、メモリMa,Mbには、最近から1024
×3×Ts前までの、1024×3回分のサンプリングデータ
Ia,Ibが常に格納されていることになる。
【0027】−CPU2への演算指示− タイマTsがオ−バする毎にCPU1は、ステップ11
でフラグMFは「1」であるので(ステップ14でセッ
ト)、レジスタnを1インクレメントし(17)、nが
64以上であるか否かをチェックする(18)。レジス
タnが64となるまではステップ4〜11−17の処理
を繰り返し、64となるとレジスタnを0に再セットし
CPU2に粒子速度の演算を指示する(16)。すなわ
ちステップ17,18は粒子速度の演算タイミングを判
定する処理でありレジスタnが64になる毎に、つまり
64回のデ−タ読取を行なう毎に(64×Ts周期
で)、CPU1はCPU2に「演算」を指示する。
でフラグMFは「1」であるので(ステップ14でセッ
ト)、レジスタnを1インクレメントし(17)、nが
64以上であるか否かをチェックする(18)。レジス
タnが64となるまではステップ4〜11−17の処理
を繰り返し、64となるとレジスタnを0に再セットし
CPU2に粒子速度の演算を指示する(16)。すなわ
ちステップ17,18は粒子速度の演算タイミングを判
定する処理でありレジスタnが64になる毎に、つまり
64回のデ−タ読取を行なう毎に(64×Ts周期
で)、CPU1はCPU2に「演算」を指示する。
【0028】−CPU2の処理内容− 図6に、図1に示すCPU2の処理内容を示す。CPU
1から速度演算の指示があると(21)、CPU2はフ
ラグNFが「1」(1回以上の演算を実行済)であるか
否かをチェックする(22)。フラグNFが1でなけれ
ば、1をセットし(23)、レジスタNに0をセットす
る(24)。そして、遅延時間Tdを算出する(2
7)。遅延時間Tdは、受光素子5の受光信号(上流側
反射光強度)に対する受光素子6の受光信号(下流側反
射光強度)の位相差である。これらの受光信号は、高炉
内の、プロ−ブ3先端に対向する装入物の降下により、
例えば図11の(a)および(b)に示すように変動す
る。レンズ36(受光素子5)の前方にあった装入物は
ΔL降下するとレンズ37(受光素子6)の前方に到達
するので、装入物の降下速度をVと仮定すると、受光素
子5の受光信号(上流側)より、ΔL/Vの時間(T
d)の後に、受光素子5の受光信号と同様なレベル変化
が受光素子6の受光信号に現われる。すなわち受光素子
5の受光信号に対する受光素子6の受光信号の位相差T
dが遅れ時間であり、Td=ΔL/Vであり、これよ
り、V=ΔL/Tdが求まる。したがって遅延時間Td
を算出することにより装入物の降下速度Vを求めること
ができる。「遅延時間の算出」(27)では、この遅延
時間Tdを算出する。
1から速度演算の指示があると(21)、CPU2はフ
ラグNFが「1」(1回以上の演算を実行済)であるか
否かをチェックする(22)。フラグNFが1でなけれ
ば、1をセットし(23)、レジスタNに0をセットす
る(24)。そして、遅延時間Tdを算出する(2
7)。遅延時間Tdは、受光素子5の受光信号(上流側
反射光強度)に対する受光素子6の受光信号(下流側反
射光強度)の位相差である。これらの受光信号は、高炉
内の、プロ−ブ3先端に対向する装入物の降下により、
例えば図11の(a)および(b)に示すように変動す
る。レンズ36(受光素子5)の前方にあった装入物は
ΔL降下するとレンズ37(受光素子6)の前方に到達
するので、装入物の降下速度をVと仮定すると、受光素
子5の受光信号(上流側)より、ΔL/Vの時間(T
d)の後に、受光素子5の受光信号と同様なレベル変化
が受光素子6の受光信号に現われる。すなわち受光素子
5の受光信号に対する受光素子6の受光信号の位相差T
dが遅れ時間であり、Td=ΔL/Vであり、これよ
り、V=ΔL/Tdが求まる。したがって遅延時間Td
を算出することにより装入物の降下速度Vを求めること
ができる。「遅延時間の算出」(27)では、この遅延
時間Tdを算出する。
【0029】図7および図11を参照して、図6に示す
「遅延時間Tdを算出」(27)の内容を説明する。先
に説明したように、メモリMaとMbには、CPU1が
計測処理を開始しから1024×3×Ts以降は、最近から10
24×3×Ts前までの、1024×3回分のサンプリングデ−
タIa,Ibが常に格納され、Ts周期で新たなサンプ
リングと、メモリMa,Mbの最も古いデ−タを捨てて
新たに得たデ−タを書込むデ−タ更新が繰返えされる。
したがってこれらのメモリMa,Mbには例えば図11
の(c)に示すようにデ−タが存在する。「遅延時間T
dを算出」(27)では、CPU2は、メモリMaのア
ドレス1〜1024×3の、1024×3個のデ−タの内の前1/3
すなわちアドレス1〜1024の1024個のデ−タ(デ−タ
群)を、遅延時間Ts検出のための比較対象デ−タと
し、このデ−タ群に相当するデ−タ群がメモリMbのど
のアドレスにあるかを探索する。
「遅延時間Tdを算出」(27)の内容を説明する。先
に説明したように、メモリMaとMbには、CPU1が
計測処理を開始しから1024×3×Ts以降は、最近から10
24×3×Ts前までの、1024×3回分のサンプリングデ−
タIa,Ibが常に格納され、Ts周期で新たなサンプ
リングと、メモリMa,Mbの最も古いデ−タを捨てて
新たに得たデ−タを書込むデ−タ更新が繰返えされる。
したがってこれらのメモリMa,Mbには例えば図11
の(c)に示すようにデ−タが存在する。「遅延時間T
dを算出」(27)では、CPU2は、メモリMaのア
ドレス1〜1024×3の、1024×3個のデ−タの内の前1/3
すなわちアドレス1〜1024の1024個のデ−タ(デ−タ
群)を、遅延時間Ts検出のための比較対象デ−タと
し、このデ−タ群に相当するデ−タ群がメモリMbのど
のアドレスにあるかを探索する。
【0030】すなわちまず、遅延時間Tdが0であると
仮定して(271)、メモリMaのアドレス1〜102
4のデ−タと、メモリMbのアドレス1〜1024のデ
−タの、等しいアドレスのデ−タ同志の差の2乗Dtを
算出し(273〜275)、これらの値の総和IDtを
算出して(273〜278)、得た総和IDtを、メモ
リMd(RAM12の一領域に割当てられている)のア
ドレス0に書込む(279)。
仮定して(271)、メモリMaのアドレス1〜102
4のデ−タと、メモリMbのアドレス1〜1024のデ
−タの、等しいアドレスのデ−タ同志の差の2乗Dtを
算出し(273〜275)、これらの値の総和IDtを
算出して(273〜278)、得た総和IDtを、メモ
リMd(RAM12の一領域に割当てられている)のア
ドレス0に書込む(279)。
【0031】次にCPU2は、遅延時間Tdが1(時間
表示では1×Ts)であると仮定して(280,28
1)、メモリMaのアドレス1〜1024のデ−タと、
メモリMbのアドレス(1+1)〜(1024+1)の
デ−タの、MaのアドレスAmaとMbのアドレスAm
b=Ama+1のデ−タ同志の差の2乗Dtを算出し(2
73〜275)、これらの値の総和IDtを算出して
(273〜278)、得た総和IDtを、メモリMdの
アドレス1に書込む(279)。
表示では1×Ts)であると仮定して(280,28
1)、メモリMaのアドレス1〜1024のデ−タと、
メモリMbのアドレス(1+1)〜(1024+1)の
デ−タの、MaのアドレスAmaとMbのアドレスAm
b=Ama+1のデ−タ同志の差の2乗Dtを算出し(2
73〜275)、これらの値の総和IDtを算出して
(273〜278)、得た総和IDtを、メモリMdの
アドレス1に書込む(279)。
【0032】以下同様にして、仮定する遅延時間を1単
位づつ大きくして、仮定する遅延時間をt(時間表示で
はt×Ts)とすると、メモリMaのアドレス1〜10
24のデ−タと、メモリMbのアドレス(1+t)〜
(1024+t)のデ−タの、MaのアドレスAmaとM
bのアドレスAmb=Ama+tのデ−タ同志の差の2乗D
tを算出し(273〜275)、これらの値の総和ID
tを算出して(273〜278)、得た総和IDtを、
メモリMdのアドレスtに書込む(279)。これを、
t=1024×2まで行なう(280)。すなわち、メ
モリMaのアドレス1〜1024のデ−タ(図11のc
のメモリMaの左1/3)と、メモリMbのアドレス1
024×2+1〜1024×3のデ−タ(図11のcの
メモリMbの左1/3)の対比演算まで行なう。
位づつ大きくして、仮定する遅延時間をt(時間表示で
はt×Ts)とすると、メモリMaのアドレス1〜10
24のデ−タと、メモリMbのアドレス(1+t)〜
(1024+t)のデ−タの、MaのアドレスAmaとM
bのアドレスAmb=Ama+tのデ−タ同志の差の2乗D
tを算出し(273〜275)、これらの値の総和ID
tを算出して(273〜278)、得た総和IDtを、
メモリMdのアドレスtに書込む(279)。これを、
t=1024×2まで行なう(280)。すなわち、メ
モリMaのアドレス1〜1024のデ−タ(図11のc
のメモリMaの左1/3)と、メモリMbのアドレス1
024×2+1〜1024×3のデ−タ(図11のcの
メモリMbの左1/3)の対比演算まで行なう。
【0033】以上により、メモリMdには、アドレス0
〜1024に、メモリMaのアドレス1〜1024の1024個
のデ−タ群と、メモリMbのアドレス1〜1024の1024個
の第0デ−タ群,アドレス2〜1025の第1デ−タ群,・
・・,アドレス1+t〜1024+tの第tデ−タ群,・・
・,アドレス1024×2+1〜1024×3の第1024デ−タ群、
全1025個のデ−タ群、それぞれとの、不一致度を表わす
デ−タ(IDt)が格納されたことになる。
〜1024に、メモリMaのアドレス1〜1024の1024個
のデ−タ群と、メモリMbのアドレス1〜1024の1024個
の第0デ−タ群,アドレス2〜1025の第1デ−タ群,・
・・,アドレス1+t〜1024+tの第tデ−タ群,・・
・,アドレス1024×2+1〜1024×3の第1024デ−タ群、
全1025個のデ−タ群、それぞれとの、不一致度を表わす
デ−タ(IDt)が格納されたことになる。
【0034】CPU2はここで、これら不一致度を表わ
すデ−タ群の中の最低値を探索して、それが格納されて
いるアドレスtmを特定する(282)。このアドレス
tmが、求める遅延時間を表わす。このアドレスtm
を、遅延時間Td=tm×Tsに変換する(283)。
すデ−タ群の中の最低値を探索して、それが格納されて
いるアドレスtmを特定する(282)。このアドレス
tmが、求める遅延時間を表わす。このアドレスtm
を、遅延時間Td=tm×Tsに変換する(283)。
【0035】図6を再度参照する。遅延時間Tdを算出
(27)するとCPU2は、バッファメモリMc(RA
M12の一領域に割り当てられている)のデータをシフ
トする(28a)。メモリMcは128個のアドレスを
有し、データシフトではアドレス2のデータをアドレス
1に書込み、アドレス3のデータをアドレス2に書込
み、・・・、アドレス128のデータをアドレス127
に書込む。そしてこれにより空いたアドレス128に、
今回算出した遅延時間Tdをメモリする(28b)。
(27)するとCPU2は、バッファメモリMc(RA
M12の一領域に割り当てられている)のデータをシフ
トする(28a)。メモリMcは128個のアドレスを
有し、データシフトではアドレス2のデータをアドレス
1に書込み、アドレス3のデータをアドレス2に書込
み、・・・、アドレス128のデータをアドレス127
に書込む。そしてこれにより空いたアドレス128に、
今回算出した遅延時間Tdをメモリする(28b)。
【0036】次に、受光信号に基づいて信号波形のピー
クピッチを分析する「ピークピッチ分析」を行う(2
9)。受光信号は受光面の直前が粒子であるか粒子間の
空隙であるかにより強度が変化するので、その周期は粒
子速度が一定であれば粒子径に依存(粒子分布ピッチに
対応)する。すなわち、図12に示すように、粒子径が
小さい場合は頻繁に信号強度の強弱が繰り返されるが、
粒子径が大きい場合はその周期は長くなる。この「ピー
クピッチ分析」(29)では粒子分布ピッチと周期の関
係に着目して受光信号を分析するものであり、その内容
を図8および図9に示す。
クピッチを分析する「ピークピッチ分析」を行う(2
9)。受光信号は受光面の直前が粒子であるか粒子間の
空隙であるかにより強度が変化するので、その周期は粒
子速度が一定であれば粒子径に依存(粒子分布ピッチに
対応)する。すなわち、図12に示すように、粒子径が
小さい場合は頻繁に信号強度の強弱が繰り返されるが、
粒子径が大きい場合はその周期は長くなる。この「ピー
クピッチ分析」(29)では粒子分布ピッチと周期の関
係に着目して受光信号を分析するものであり、その内容
を図8および図9に示す。
【0037】図8,図9および図13を参照して、図6
に示す「ピークピッチ分析」(29)の内容を説明す
る。まず図8を参照すると、メモリMaのデータの上ピ
ーク値をアドレスに対応付けてメモリMfにメモリする
(291)。なお以下に示す処理は、上流側の第1受光
素子5が発生する受光信号に基づいて行うが、下側の受
光素子6が発生する受光信号に基づいて行ってもよい。
次に、メモリMfの、データがないアドレス,すなわち
デ−タが0値のアドレスに対しては、該アドレスより小
さい側と大きい側のピ−ク値の値を用いて補間法により
デ−タを算出して書き込む(292)。次に、メモリM
fのデ−タを平滑化してMfに更新メモリした後(29
3)、データを「0」又は「1」に2値化してレジスタ
Mgに書き込む(294)。このデ−タの「1」の連な
り(領域)はある層(たとえば鉱石層)を、「0」の連
なり(領域)は他の層(コ−クス層)を意味する。次
に、各領域Ri毎のピークピッチTiを算出する(29
5)。
に示す「ピークピッチ分析」(29)の内容を説明す
る。まず図8を参照すると、メモリMaのデータの上ピ
ーク値をアドレスに対応付けてメモリMfにメモリする
(291)。なお以下に示す処理は、上流側の第1受光
素子5が発生する受光信号に基づいて行うが、下側の受
光素子6が発生する受光信号に基づいて行ってもよい。
次に、メモリMfの、データがないアドレス,すなわち
デ−タが0値のアドレスに対しては、該アドレスより小
さい側と大きい側のピ−ク値の値を用いて補間法により
デ−タを算出して書き込む(292)。次に、メモリM
fのデ−タを平滑化してMfに更新メモリした後(29
3)、データを「0」又は「1」に2値化してレジスタ
Mgに書き込む(294)。このデ−タの「1」の連な
り(領域)はある層(たとえば鉱石層)を、「0」の連
なり(領域)は他の層(コ−クス層)を意味する。次
に、各領域Ri毎のピークピッチTiを算出する(29
5)。
【0038】図9に「ピ−クピッチTiの算出」(29
5)の内容を示す。ここではまずメモリMaのデータを
微分し直流成分をカットしてメモリMiに格納する(2
951)。次にメモリMiのデータを平滑化しこれをメ
モリMiに更新メモリした後(2952)、メモリMi
上のデ−タの隣り合うピーク間の距離T1,T2,T3,
・・,をそれぞれ算出してそれらの平均値Tiを算出す
る(2953)。この平均値Tiが降下する粒子の周期
を示す。そして、平均値Tiに、後述するステップ34
で求めレジスタに保持している粒子速度Vmを掛けて粒
子分布ピッチPthiを算出する(2954)。ここで、
図13に示すように検出基準となる粒子A,Bの間には
粒子C,Dが割り込むため算出された粒子分布ピッチP
thiはそのまま粒子の大きさ(直径)を示すものではな
い。そこで、粒子分布ピッチPthiに補正を加えて粒子
直径Dpを求める(2955)。粒子分布ピッチと粒径
の関係を図13に示す。CPU2は次に、処理対象の領
域Riに対応付けてメモリMhに粒子分布ピッチPthiと
粒子直径Dpiを記憶し(2956)、領域Ri対応のレ
ジスタMgの内容が「1」か否かチェックし(295
7)すなわち処理対象領域Riが「0」領域(コ−クス
層)か「1」領域(鉱石層)かをチェックして、「1」
領域であればメモリMe1のデータを1アドレス分シフト
する(2958a)。メモリMe1は128個のアドレス
を有し、データシフトではアドレス2のデータをアドレ
ス1に書込み、,・・,アドレス128のデータをアド
レス127に書込む。そしてこれにより空いたアドレス
128に最新の平均値Ti(2953で算出した値)を
メモリする(2958b)。また、レジスタMgの内容
が「0」であったときには、メモリMe0のデータを1ア
ドレス分シフトして(2959a)そしてこれにより空
いたアドレス128に最新の平均値Tiをメモリする
(2959b)。以上の処理を、図8に示す領域R2〜
R8(i=2〜8)に対して行なう。
5)の内容を示す。ここではまずメモリMaのデータを
微分し直流成分をカットしてメモリMiに格納する(2
951)。次にメモリMiのデータを平滑化しこれをメ
モリMiに更新メモリした後(2952)、メモリMi
上のデ−タの隣り合うピーク間の距離T1,T2,T3,
・・,をそれぞれ算出してそれらの平均値Tiを算出す
る(2953)。この平均値Tiが降下する粒子の周期
を示す。そして、平均値Tiに、後述するステップ34
で求めレジスタに保持している粒子速度Vmを掛けて粒
子分布ピッチPthiを算出する(2954)。ここで、
図13に示すように検出基準となる粒子A,Bの間には
粒子C,Dが割り込むため算出された粒子分布ピッチP
thiはそのまま粒子の大きさ(直径)を示すものではな
い。そこで、粒子分布ピッチPthiに補正を加えて粒子
直径Dpを求める(2955)。粒子分布ピッチと粒径
の関係を図13に示す。CPU2は次に、処理対象の領
域Riに対応付けてメモリMhに粒子分布ピッチPthiと
粒子直径Dpiを記憶し(2956)、領域Ri対応のレ
ジスタMgの内容が「1」か否かチェックし(295
7)すなわち処理対象領域Riが「0」領域(コ−クス
層)か「1」領域(鉱石層)かをチェックして、「1」
領域であればメモリMe1のデータを1アドレス分シフト
する(2958a)。メモリMe1は128個のアドレス
を有し、データシフトではアドレス2のデータをアドレ
ス1に書込み、,・・,アドレス128のデータをアド
レス127に書込む。そしてこれにより空いたアドレス
128に最新の平均値Ti(2953で算出した値)を
メモリする(2958b)。また、レジスタMgの内容
が「0」であったときには、メモリMe0のデータを1ア
ドレス分シフトして(2959a)そしてこれにより空
いたアドレス128に最新の平均値Tiをメモリする
(2959b)。以上の処理を、図8に示す領域R2〜
R8(i=2〜8)に対して行なう。
【0039】図8を再度参照する。各領域Riのピーク
ピッチTiの算出処理(295)が終了すると、両端の
領域R1,R9を除く各領域(R2〜R8)の厚みTH
Riを算出する(296)。例えば領域R8の厚みTH
R8は、領域R7の始端の「0」の前の「1」のアドレ
ス(領域R8の終端アドレス)と、領域R9の終端の
「0」のアドレスとの差に、Tsならびに後述するステ
ップ34で求めレジスタに保持している粒子速度Vmを
乗算して求める。次にCPU2は、メモリMfのデー
タ,層厚THRi,粒子分布ピッチPthiおよび粒子直
径Dpiを出力する(297)。
ピッチTiの算出処理(295)が終了すると、両端の
領域R1,R9を除く各領域(R2〜R8)の厚みTH
Riを算出する(296)。例えば領域R8の厚みTH
R8は、領域R7の始端の「0」の前の「1」のアドレ
ス(領域R8の終端アドレス)と、領域R9の終端の
「0」のアドレスとの差に、Tsならびに後述するステ
ップ34で求めレジスタに保持している粒子速度Vmを
乗算して求める。次にCPU2は、メモリMfのデー
タ,層厚THRi,粒子分布ピッチPthiおよび粒子直
径Dpiを出力する(297)。
【0040】図6を再度参照する。フラグレジスタAF
の内容が「1」(128回以上の遅延時間算出とメモリ
Mcへの書込みおよび128回以上のピークピッチの算
出とメモリMe0,Me1への書込みを実行済)であるかを
チェックし(30)、「1」でないと、CPU1より再
度の演算指示が到来するのを待ち、これが到来するとま
たステップ22〜30を同様に実行する。このようにし
てCPU1より演算の指示がある毎に、ステップ22〜
30の遅延時間算出とメモリMcへの書込みおよびピー
クピッチの算出とメモリMe0,Me1への書込みを行な
う。これを一回行なうごとにレジスタNの内容を1イン
クレメントし(32)、レジスタNの内容Nが128に
なると、すなわち128回の遅延時間算出とメモリMc
への書込みおよび128回のピークピッチの算出とメモ
リMe0,Me1への書込みが終了するとき、すなわち、6
4×Ts×128の時間の間の、64×Ts周期で算出
した、128個の遅延時間デ−タのメモリMcへの格納
および128個のピークピッチデ−タのメモリMe0,M
e1への格納が終了するとき、フラグレジスタAFに
「1」をセットする(26)。
の内容が「1」(128回以上の遅延時間算出とメモリ
Mcへの書込みおよび128回以上のピークピッチの算
出とメモリMe0,Me1への書込みを実行済)であるかを
チェックし(30)、「1」でないと、CPU1より再
度の演算指示が到来するのを待ち、これが到来するとま
たステップ22〜30を同様に実行する。このようにし
てCPU1より演算の指示がある毎に、ステップ22〜
30の遅延時間算出とメモリMcへの書込みおよびピー
クピッチの算出とメモリMe0,Me1への書込みを行な
う。これを一回行なうごとにレジスタNの内容を1イン
クレメントし(32)、レジスタNの内容Nが128に
なると、すなわち128回の遅延時間算出とメモリMc
への書込みおよび128回のピークピッチの算出とメモ
リMe0,Me1への書込みが終了するとき、すなわち、6
4×Ts×128の時間の間の、64×Ts周期で算出
した、128個の遅延時間デ−タのメモリMcへの格納
および128個のピークピッチデ−タのメモリMe0,M
e1への格納が終了するとき、フラグレジスタAFに
「1」をセットする(26)。
【0041】そしてフラグレジスタAFに「1」をセッ
トした後は、遅延時間を算出し(27)これをメモリM
cに書込み(28a,28b)、ピークピッチ分析の処
理が終了する(29)と、メモリMcの全デ−タ(12
8個)の平均値Tmを算出し(33)、この平均値Tm
とΔLより、降下速度Vm=ΔL/Tmを算出する(3
4)。そして、粒子径を出力する処理を行う(35)。
トした後は、遅延時間を算出し(27)これをメモリM
cに書込み(28a,28b)、ピークピッチ分析の処
理が終了する(29)と、メモリMcの全デ−タ(12
8個)の平均値Tmを算出し(33)、この平均値Tm
とΔLより、降下速度Vm=ΔL/Tmを算出する(3
4)。そして、粒子径を出力する処理を行う(35)。
【0042】図10を参照して、「粒子径出力」(3
5)の内容を説明する。メモリMe1,メモリMe0に蓄積
された128個分のピークピッチTiの平均値Tphm1,
Tphm0を算出し(351)、これに対応する粒子分布ピ
ッチの平均値Pthm1(=Tphm1×Vm),Pthm0(=Tph
m0×Vm)とこれに補正を加えた粒子直径の平均値Dpm
1,Dpm0を求めて、これらPthm1,Pthm0,Dpm1,Dp
m0と粒子の降下速度Vmを出力する(352〜35
4)。操作・表示ボ−ド11の表示手段に降下速度,粒
子分布ピッチおよび粒子径が更新表示される。
5)の内容を説明する。メモリMe1,メモリMe0に蓄積
された128個分のピークピッチTiの平均値Tphm1,
Tphm0を算出し(351)、これに対応する粒子分布ピ
ッチの平均値Pthm1(=Tphm1×Vm),Pthm0(=Tph
m0×Vm)とこれに補正を加えた粒子直径の平均値Dpm
1,Dpm0を求めて、これらPthm1,Pthm0,Dpm1,Dp
m0と粒子の降下速度Vmを出力する(352〜35
4)。操作・表示ボ−ド11の表示手段に降下速度,粒
子分布ピッチおよび粒子径が更新表示される。
【0043】粒子径出力の処理が終了するとステップ2
1に戻り、ステップ21〜35の処理を繰り返す。CP
U1が、64×Ts周期で演算指示をCPU2に与える
ので、64×Ts周期で、遅延時間の算出(27,28
a,28b)と降下速度Vm,粒子分布ピッチPthm1,
Pthm0および粒子直径Dpm1,Dpm0の算出(29,33
〜35)、および算出した降下速度,粒子分布ピッチお
よび粒子径の更新表示、が繰返えされる。
1に戻り、ステップ21〜35の処理を繰り返す。CP
U1が、64×Ts周期で演算指示をCPU2に与える
ので、64×Ts周期で、遅延時間の算出(27,28
a,28b)と降下速度Vm,粒子分布ピッチPthm1,
Pthm0および粒子直径Dpm1,Dpm0の算出(29,33
〜35)、および算出した降下速度,粒子分布ピッチお
よび粒子径の更新表示、が繰返えされる。
【0044】上述のプローブ3を図14に示すように、
高炉シャフト部の所要高さ位置で炉壁14に通し、水平
方向の挿入長を調整することにより、高炉の周辺部およ
び中心部での装入物の粒子径を測定しうる。このような
測定により、例えば図15に示すような、高炉内装入物
の粒子径分布を得ることができる。
高炉シャフト部の所要高さ位置で炉壁14に通し、水平
方向の挿入長を調整することにより、高炉の周辺部およ
び中心部での装入物の粒子径を測定しうる。このような
測定により、例えば図15に示すような、高炉内装入物
の粒子径分布を得ることができる。
【0045】なお、本実施例では受光素子5の受光信号
と受光素子6の受光信号との相関判定から遅延時間Td
を算出し、算出した遅延時間Tdと所定間隔ΔLから粒
子の降下速度Vmを算出したが、これに限らず他の速度
センサを使用して求めてもよい。この場合には、本実施
例で示したように上流側および下流側に光ファイバ3
8,39および受光素子5,6をそれぞれ設ける必要は
なく、粒子径出力用に光ファイバとその受光素子の一組
を設けることで足りる。また、粒子の降下速度が一定と
仮定できる対象においては、周期と粒子径は比例するの
で周期から粒子径の相対変化を検出することができる。
と受光素子6の受光信号との相関判定から遅延時間Td
を算出し、算出した遅延時間Tdと所定間隔ΔLから粒
子の降下速度Vmを算出したが、これに限らず他の速度
センサを使用して求めてもよい。この場合には、本実施
例で示したように上流側および下流側に光ファイバ3
8,39および受光素子5,6をそれぞれ設ける必要は
なく、粒子径出力用に光ファイバとその受光素子の一組
を設けることで足りる。また、粒子の降下速度が一定と
仮定できる対象においては、周期と粒子径は比例するの
で周期から粒子径の相対変化を検出することができる。
【0046】また、本実施例では粒子径を決定する周期
を時間領域で平均化して算出したが、該周期を周波数領
域で平均化、例えば周波数演算処理機構で周波数解析を
行い平均周波数(強度が最大となる周波数)を求めても
よい。平均周波数を求めるにあたっては条件により最低
周波数や最高周波数近傍にピークがきて装入粒子径や降
下速度から推算して、例えば図16の(a)に示すよう
に、妥当な周波数にピークを持たない場合も生じる。し
かし、これらはソフト処理で対処できる。例えば低周波
数にピークが検出された場合は、ハイパスフィルタをか
けたり微分して低周波成分を除去することにより、図1
6の(b)に示すような、装入粒子径と良好な対応がと
れる平均周波数が求まる。
を時間領域で平均化して算出したが、該周期を周波数領
域で平均化、例えば周波数演算処理機構で周波数解析を
行い平均周波数(強度が最大となる周波数)を求めても
よい。平均周波数を求めるにあたっては条件により最低
周波数や最高周波数近傍にピークがきて装入粒子径や降
下速度から推算して、例えば図16の(a)に示すよう
に、妥当な周波数にピークを持たない場合も生じる。し
かし、これらはソフト処理で対処できる。例えば低周波
数にピークが検出された場合は、ハイパスフィルタをか
けたり微分して低周波成分を除去することにより、図1
6の(b)に示すような、装入粒子径と良好な対応がと
れる平均周波数が求まる。
【0047】
【発明の効果】以上の通り本発明によれば、検出部は投
受光の伝送のみでよいので多数の光ファイバを整然と配
置したファイバスコープと比べ検出部の小型化,低価格
化が図れる。また、光ファイバ(38)の受光面から遠
く離れた粒子に関しては投受光量が減衰するため粒子の
重なりなどの影響は少ない。
受光の伝送のみでよいので多数の光ファイバを整然と配
置したファイバスコープと比べ検出部の小型化,低価格
化が図れる。また、光ファイバ(38)の受光面から遠
く離れた粒子に関しては投受光量が減衰するため粒子の
重なりなどの影響は少ない。
【0048】また画像処理のような複雑な処理は必要な
く高速演算処理が可能であり、移動層の運転上に重要な
指標となる層内の粒子径およびその分布の変化をリアル
タイムに検出できる。
く高速演算処理が可能であり、移動層の運転上に重要な
指標となる層内の粒子径およびその分布の変化をリアル
タイムに検出できる。
【0049】本発明の好ましい実施例によれば、受光信
号Iaは、装入物粒子分布ピッチの算出用および装入物
の降下速度の算出用に兼用されるので実用的である。こ
こで得られる降下速度は、粒子の流れ速度であり、高炉
内各所で細かに降下速度を測定しうる。従来の層単位の
検出にくらべて、速度変化に対する計測の遅れが実質上
なくなり、また、コ−クスと鉱石の混合層であっても、
正確に速度が検出され、信頼性が高い測定値が得られ
る。
号Iaは、装入物粒子分布ピッチの算出用および装入物
の降下速度の算出用に兼用されるので実用的である。こ
こで得られる降下速度は、粒子の流れ速度であり、高炉
内各所で細かに降下速度を測定しうる。従来の層単位の
検出にくらべて、速度変化に対する計測の遅れが実質上
なくなり、また、コ−クスと鉱石の混合層であっても、
正確に速度が検出され、信頼性が高い測定値が得られ
る。
【図1】 本発明装置の一実施例の構成概要を示すブロ
ック図である。
ック図である。
【図2】 図1に示すプローブ3の先端面を示す拡大正
面図である。
面図である。
【図3】 図1に示すプローブ3の拡大縦断面図であ
り、図2のIII−III線断面を示す。
り、図2のIII−III線断面を示す。
【図4】 図3に示す光ファイバケ−ブル38,39の
拡大横断面図である。
拡大横断面図である。
【図5】 図1に示すCPU1の処理動作を示すフロー
チャートである。
チャートである。
【図6】 図1に示すCPU2の処理動作を示すフロー
チャートである。
チャートである。
【図7】 図6に示す「遅延時間Tdを算出」(27)
の処理内容を示すフローチャートである。
の処理内容を示すフローチャートである。
【図8】 図6に示す「ピークピッチ分析」(29)の
処理内容を示すフローチャートである。
処理内容を示すフローチャートである。
【図9】 図8に示す「各領域RiのピークピッチTi
の算出」(295)の処理内容を示すフローチャートで
ある。
の算出」(295)の処理内容を示すフローチャートで
ある。
【図10】 図6に示す「粒子径出力」(35)の処理
内容を示すフローチャートである。
内容を示すフローチャートである。
【図11】 図4に示す光ファイバケ−ブル38,39
の反射光受光強度を示すグラフであり、(a)は光ファ
イバケ−ブル38のものを、(b)は光ファイバケ−ブ
ル39のものを示す。(c)は、光ファイバケ−ブル3
8,39が受光し図1に示す受光素子5,6が電気信号
に変換した2組の受光強度デ−タの、それらを格納し
た、図1のRAM12に割当てられたメモリMa,Mb
上でのデ−タ分布を示すブロック図である。
の反射光受光強度を示すグラフであり、(a)は光ファ
イバケ−ブル38のものを、(b)は光ファイバケ−ブ
ル39のものを示す。(c)は、光ファイバケ−ブル3
8,39が受光し図1に示す受光素子5,6が電気信号
に変換した2組の受光強度デ−タの、それらを格納し
た、図1のRAM12に割当てられたメモリMa,Mb
上でのデ−タ分布を示すブロック図である。
【図12】 粒子径を変化させた場合の出力強度を示す
グラフである。
グラフである。
【図13】 降下する粒子を示す断面図である。
【図14】 図1に示すプローブ3が挿入された高炉1
4を示す断面図である。
4を示す断面図である。
【図15】 図9に示すプローブ3先端位置と、測定し
た装入物粒子径の関係を示すグラフである。
た装入物粒子径の関係を示すグラフである。
【図16】 周波数と強度との関係を示すグラフであ
り、(a)はピークをもたない場合,(b)はピークを
もつ場合を示す。
り、(a)はピークをもたない場合,(b)はピークを
もつ場合を示す。
1:CPU 2:演算用CPU 3:プローブ 31:外筒 32:仕切板 33:内側空間 34:外側空間 35:保護ガラス 36:レンズ 37:レンズ 38:光ファイバケ−ブル 38e:光投射用
光ファイバ 38r:受光用光ファイバ 39:光ファイバ
ケ−ブル 39e:光投射用光ファイバ 39r:受光用光
ファイバ 4:レーザ光源 5:受光素子 6:受光素子 14:高炉
光ファイバ 38r:受光用光ファイバ 39:光ファイバ
ケ−ブル 39e:光投射用光ファイバ 39r:受光用光
ファイバ 4:レーザ光源 5:受光素子 6:受光素子 14:高炉
Claims (4)
- 【請求項1】粒子群への光の投射と光ファイバを介した
粒径代表値よりも小さい径の視野での粒子群の反射光の
検出を、粒子群と光ファイバが相対的に移動する状態
で、連続もしくは所定周期で行なって反射光強度を表わ
す受光信号Iaを得、 受光信号Iaのレベル変動周期Tを算出し、 該周期Tと前記相対移動の速度Vより、粒子群中の粒子
分布ピッチT・Vを算出する、粒子分布ピッチ測定方
法。 - 【請求項2】以下により前記相対移動の速度Vを測定す
る、請求項1記載の粒子分布ピッチ測定方法:前記光の
反射点から、前記視野径より大きく粒径代表値の数十倍
以内の所定距離ΔL離れた位置への光の投射と光ファイ
バを介した粒径代表値よりも小さい径の視野での粒子群
の反射光の検出を連続もしくは所定周期で行なって反射
光強度を表わす第2受光信号Ibを得、 前記受光信号Iaに対する第2受光信号Ibの位相差T
dを算出し、該位相差Tdと前記所定距離ΔLより速度
V=ΔL/Tdを算出する。 - 【請求項3】粒子群に光を投射し粒子群の反射光を粒径
代表値よりも小さい径の視野で受けるための光ファイ
バ,および、該光ファイバを支持する支持部材、を含
み、粒子群に差し込まれる測定プロ−ブ;光ファイバが
受けた反射光をその強度を表わす電気信号に変換する変
換手段;電気信号のレベル変動周期Tを検出する周期検
出手段;および、 粒子群と測定プロ−ブの間の相対移動速度と前記周期T
より粒子分布ピッチP=T・Vを算出する演算手段;を
備える、粒子分布ピッチ測定装置。 - 【請求項4】前記測定プローブは更に、前記光ファイバ
と同様に粒子群に光を投射し粒子群の反射光を装入物粒
径代表値よりも小さい径の視野で受けるための第2組の
光ファイバを有し請求項3に記載の光ファイバと第2組
の光ファイバを、粒子群に対するそれらの測定視野の距
離を所定距離ΔLとして支持し;装置は更に、 第2組の光ファイバが受けた反射光をその強度を表わす
第2電気信号に変換する第2変換手段;請求項3に記載
の電気信号に対する第2電気信号の位相差Tdを検出す
る位相差検出手段;および、 前記所定距離ΔLと位相差Tdより粒子と測定プロ−ブ
の間の相対移動速度V=ΔL/Tdを算出する速度演算
手段;を備える請求項3記載の粒子分布ピッチ測定装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4265127A JPH06117989A (ja) | 1992-10-02 | 1992-10-02 | 粒子分布ピッチ測定方法および装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4265127A JPH06117989A (ja) | 1992-10-02 | 1992-10-02 | 粒子分布ピッチ測定方法および装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH06117989A true JPH06117989A (ja) | 1994-04-28 |
Family
ID=17412998
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4265127A Pending JPH06117989A (ja) | 1992-10-02 | 1992-10-02 | 粒子分布ピッチ測定方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH06117989A (ja) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5926038A (ja) * | 1982-08-03 | 1984-02-10 | Nippon Steel Corp | 高炉内装入物検出方法および装置 |
| JPS5938313A (ja) * | 1982-08-25 | 1984-03-02 | Nippon Steel Corp | 高炉内融着帯根部検出方法 |
| JPS6237748B2 (ja) * | 1980-07-07 | 1987-08-13 | Tokyo Shibaura Electric Co |
-
1992
- 1992-10-02 JP JP4265127A patent/JPH06117989A/ja active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6237748B2 (ja) * | 1980-07-07 | 1987-08-13 | Tokyo Shibaura Electric Co | |
| JPS5926038A (ja) * | 1982-08-03 | 1984-02-10 | Nippon Steel Corp | 高炉内装入物検出方法および装置 |
| JPS5938313A (ja) * | 1982-08-25 | 1984-03-02 | Nippon Steel Corp | 高炉内融着帯根部検出方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU608807B2 (en) | Apparatus for monitoring bloodstream | |
| JPS6161616B2 (ja) | ||
| JPH0428005Y2 (ja) | ||
| US4993835A (en) | Apparatus for detecting three-dimensional configuration of object employing optical cutting method | |
| JPH0450611A (ja) | 距離測定装置 | |
| CN101310189B (zh) | 用于测量物体相对于表面的运动速度的光学设备 | |
| CA1168437A (en) | Process and device for the contact free measurement of a dimension | |
| EP0565320A2 (en) | Method and an apparatus for sensing wheel parameters in a wheel balancing machine | |
| KR19990081876A (ko) | 거리 및 위치 결정용 장치 | |
| KR0125442B1 (ko) | 표면 거칠기 광탐지 방법 및 장치 | |
| JPS6239379B2 (ja) | ||
| JP2517378B2 (ja) | 自由落下物体の大きさおよび/または形を決定するための方法および配置 | |
| JP2729014B2 (ja) | 高炉装入物の降下速度測定方法および装置 | |
| JP2501237B2 (ja) | 鋼管端部の外径・肉厚測定装置 | |
| JPH0621775B2 (ja) | 表面検査装置 | |
| JPH0635921B2 (ja) | スペックル測長計 | |
| JP3097414B2 (ja) | 形状認識装置 | |
| RU2052771C1 (ru) | Устройство для контроля изделий | |
| JPH085466A (ja) | 表面温度分布測定方法及び該方法に用いられる装置 | |
| JPH0621776B2 (ja) | 鮮映性測定方法 | |
| JPH0735515A (ja) | 対象物の直径測定装置 | |
| JPH0781960B2 (ja) | 連続状物体の欠点検出装置 | |
| JPH038487B2 (ja) | ||
| SU945648A1 (ru) | Способ измерени геометрических размеров прозрачных труб | |
| JPS5838808A (ja) | 棒状物体の形体測定方法 |