JP3143500B2 - 火炎分析器及び火炎特性決定方法 - Google Patents

火炎分析器及び火炎特性決定方法

Info

Publication number
JP3143500B2
JP3143500B2 JP03239682A JP23968291A JP3143500B2 JP 3143500 B2 JP3143500 B2 JP 3143500B2 JP 03239682 A JP03239682 A JP 03239682A JP 23968291 A JP23968291 A JP 23968291A JP 3143500 B2 JP3143500 B2 JP 3143500B2
Authority
JP
Japan
Prior art keywords
value
flame
linearity
determining
threshold
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.)
Expired - Lifetime
Application number
JP03239682A
Other languages
English (en)
Other versions
JPH04244922A (ja
Inventor
シー.コーミアー ケネス
Original Assignee
ファイァーライ、インコーポレイテッド
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 ファイァーライ、インコーポレイテッド filed Critical ファイァーライ、インコーポレイテッド
Publication of JPH04244922A publication Critical patent/JPH04244922A/ja
Application granted granted Critical
Publication of JP3143500B2 publication Critical patent/JP3143500B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N5/00Systems for controlling combustion
    • F23N5/02Systems for controlling combustion using devices responsive to thermal changes or to thermal expansion of a medium
    • F23N5/08Systems for controlling combustion using devices responsive to thermal changes or to thermal expansion of a medium using light-sensitive elements
    • F23N5/082Systems for controlling combustion using devices responsive to thermal changes or to thermal expansion of a medium using light-sensitive elements using electronic means
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2223/00Signal processing; Details thereof
    • F23N2223/08Microprocessor; Microcomputer
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2223/00Signal processing; Details thereof
    • F23N2223/10Correlation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2229/00Flame sensors
    • F23N2229/08Flame sensors detecting flame flicker
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2229/00Flame sensors
    • F23N2229/22Flame sensors the sensor's sensitivity being variable

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Control Of Combustion (AREA)

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明はボイラー、火炉もしくは
類似の燃焼装置と共に使用する火炎センサに関し、より
詳細には多バーナシステムにおける火炎の存在及び特性
の表示を与えるようなセンサに関する。
【0002】
【従来の技術】大型ボイラ及び火炉は複数の火炎を発生
する数台のバーナを使用している。バーナの電気制御シ
ステムには火炎の存在を検出し火炎特性に関する情報を
提供する機構が含まれていることが多い。このような情
報はコントロールシステムにおいてバーナの安全運転を
調整するのに使用される。このようなコントロールシス
テムには単もしくは多バーナ装置内のバーナ火炎の存否
を検出する火炎スキャナが組み込まれている。バーナが
オンとされてバーナのスロートから燃料が噴射される
と、スキャナが火炎を監視して火炎の状態、強さ及び種
類を表わす信号を発生する。従って、火炎スキャナは被
精査バーナからの火炎と隣接バーナの火炎と他の背景状
態とを識別できることが必要である。
【0003】さまざまなスキャナが火炎に向けて火炎か
らの光の強さに比例する振幅を有する電気信号を発生す
る光センサを利用している。センサ信号の振幅は、帯域
濾波もしくは高域濾波された後、バーナ火炎のオンオフ
状態間を識別するのに使用される。しかしながら信号の
大きさはダンパ位置、火炎とセンサの近接度、燃料の種
類、及び燃料のBTU含有量等のいくつかの変数に依存
する。同様に、多バーナシステム内の他の火炎によりセ
ンサ信号内に広範に変動する背景信号成分が生成され
る。振幅依存火炎センサの主要な問題点は火炎オフ及び
火炎オン状態においてセンサ信号の大きさが変動するこ
とである。その結果、火炎オン及び火炎オフ間のセンサ
信号振幅の差が小さ過ぎて火炎状態間の識別を行うため
の信頼度の高い閾値を設定できなくなることがある。
【0004】スキャナが異なる火炎状態間を適切に識別
できなくなると、コントロールシステムは全バーナを誤
遮断したりオペレータがバーナを始動できなくなること
がある。さらに、被感知近接バーナが点火される表示と
して背景信号成分が誤って解釈されることにより、誤っ
た決定がなされることもある。このような状況では、近
接バーナ火炎は消炎することもあるが、火炎スキャナは
バーナ火炎がオンであることを示す信号をコントロール
システムへ発生する。この誤表示により燃料弁は開放さ
れたままとされバーナ室内に爆発ガスが蓄積されること
になる。従って、コントロールシステムには被感知バー
ナ火炎により生じる信号を多バーナシステム内の他の火
炎による信号から識別する機構を設けなければならな
い。
【0005】
【発明が解決しようする課題】火炎分析器が燃焼装置か
らの放射を検出して、例えば火炎の存在等の火炎特性を
感知する。センサが検出された放射を示す電気信号を発
生する。信号は、実施例においてフーリエ変換により、
信号内に存在する周波数成分スペクトルの大きさを表わ
す複数の振幅値へ変換される。これらの周波数成分は火
炎力の時間変化により生成される。火炎特性は複数の周
波数成分振幅値分布の形状により決定される。
【0006】好ましくは、特性は複数の振幅値に対する
対数値を引き出すことにより決定される。周波数成分振
幅の対数値分布の線型性の度合いが計算される。線型度
は積分線型誤差、勾配差及び線型度回帰相関等の一つ以
上のパラメータにより定義される。
【0007】火炎の存在を決定するために、計算された
各パラメータに対する一連の値がそのパラメータに対す
る別々の閾値と比較される。閾値よりも上及び下の値の
量が表にされる。実施例において、それぞれの閾値より
も低いパラメータ値の量が平均化されて第1の平均が求
められる。同様に、それぞれの閾値よりも上のパラメー
タ値の量が平均化されて第2の平均が求められる。
【0008】第1の平均が第1の基準値を越える場合に
は火炎が存在しないと決定され、第2の平均値が第2の
閾値を越える場合には火炎が存在すると決定される。
【0009】
【課題を解決するための手段】本発明の一般的目的はセ
ンサの火炎に対する近接度、煙道ダンパ位置、及び燃料
の種類とそのBTU含有量により火炎感知が影響を受け
ることのない、火炎特性決定装置及び方法を提供するこ
とである。
【0010】本発明のもう一つの目的は火炎のフリッカ
周波数スペクトル成分の時間による形状変化に基づいた
火炎検出技術を提供することである。
【0011】さらにもう一つの目的は火炎センサ信号の
周波数スペクトルの分析、特に周波数成分振幅の対数値
分布の線型性分析を行うことである。
【0012】
【実施例】図1に本発明に従ったバーナ火炎分析器10
の電子回路の実施例を示す。硫化鉛検出器等のセンサ1
2が被検出火炎から放出される火放射を受けるように配
置されている。センサ12は所望の火炎からの放射を受
けるように示されているが、代表的には多バーナ燃焼装
置内の他の火炎からの放射も受けるようにされている。
センサ12の一方の導線は負バイアス電圧源(−V)に
接続され、他方の導線は定利得プリアンプ14の反転入
力に接続されている。プリアンプ14の非反転入力は抵
抗器16を介して接地され、プリアンプの出力は帰還抵
抗器17を介してその反転入力に接続されている。
【0013】プリアンプ14の出力はアンプ20及び付
属部品からなる回路の自動利得制御部に接続されてい
る。特に、プリアンプ14の出力は抵抗器18及びコン
デンサ19の直列接続を介して第1のアンプ20の反転
入力に接続されている。コンデンサ19はセンサ20に
印加される直流バイアス電圧及びプリアンプ14のオフ
セット電圧が第1のアンプに印加されないように切り離
す。アンプ20の非反転入力は抵抗器21を介して接地
されている。第1のアンプ20の出力は固定抵抗器22
及びホトレジスタ24を介して非反転入力に接続されて
いる。ホトレジスタ24は発光ダイオード25からの光
を受光し、ホトレジスタ24の抵抗値は発光ダイオード
25を流れる電流に逆比例する。
【0014】抵抗器26は第1のアンプ20の出力を第
2のアンプ28の非反転入力へ接続し、その非反転入力
は抵抗器29により接地されている。第2のアンプ2
8、ダイオード30及び31、及び帰還抵抗器32によ
り第1のアンプ20の出力信号が全波整流される。全波
整流された信号は第3のアンプ34により形成されるロ
ーパスフィルタへ抵抗器33を介して接続され、整流さ
れた信号はその反転入力へ加えられる。第3のアンプ3
4の非反転入力は抵抗器35を介して接地される。第3
のアンプ34の出力は抵抗器36とコンデンサ38の並
列接続を介してその反転入力に接続されている。ローパ
ス濾波により第1のアンプ20からの交番信号の振幅に
比例する直流信号が生成される。
【0015】この直流信号は抵抗器39及び41を介し
て差動アンプ40の非反転入力へ加えられる。差動アン
プ40は第3のアンプ34の出力を、抵抗器42を介し
て差動アンプ40の反転入力に加えられる基準電圧V
REF により定まる設定点と比較する。差動アンプ40の
出力は帰還抵抗器43を介して非反転入力へ接続され
る。差動アンプの出力は基準電圧VREF とそれ自体が第
1のアンプ20の信号出力に比例する第3のアンプ34
からの直流電圧との差に比例するエラー電圧を与える。
【0016】エラー電圧は発光ダイオード(LED)2
5のアノードに接続された抵抗器44を介してエラー電
流信号へ変換される。LED25を駆動するこのエラー
電流信号により第1のアンプ20の負帰還利得制御が行
われる。その結果、プリアンプ14からの信号の交流振
幅が低下するため、第1のアンプ20の利得が増大す
る。制御利得の変化率はローパスフィルタの抵抗器36
及びコンデンサ38により形成されるRC回路網の時定
数により定まる。この時定数は火炎分析に使用される最
低周波数よりも5倍緩く選定される。回路設計により5
0の利得制御を行って異なるタイプの火炎及び点火条件
により生成される広範な火炎に対して良好な信号/ノイ
ズ比が維持される。
【0017】第3のアンプ34の反転入力は抵抗器46
を介して第1のアンプ20の出力ノード48に接続され
ている。ノード48は自動利得制御アンプ回路の出力を
形成し、マイクロコンピュータ50のアナログ入力に接
続されている。マイクロコンピュータはマイクロプロセ
ッサの他に、アナログ入力が接続されたアナログ/デジ
タルコンバータを含む集積回路である。マイクロコンピ
ュータ50は1組のデータ・アドレス及びコントロール
バス51,52,53がそれぞれ接続された並列入出力
ポートも含んでいる。前記したように、火炎分析を行う
のにマイクロコンピュータは読取専用メモリ(ROM)
54内に記憶されたプログラムを実行する。マイクロコ
ンピュータ50のアナログ入力に受信される信号で表わ
されるセンサ12からのデータはランダムアクセスメモ
リ(RAM)56内に記憶される。さらに、ランダムア
クセスメモリ56はマイクロプロセッサ50が行う分析
の中間及び最終結果の記憶場所も提供する。処理結果
は、入出力インターフェイス回路58を介して、燃焼装
置のバーナ制御回路等の外部装置へ与えられる。ROM
54,RAM56及び入出力インターフェイス回路56
はバス51〜53へ接続されている。
【0018】図1に示す火炎分析器10の動作中に、セ
ンサ12はバーナ組立体からの放射を電気信号へ変換す
る。センサ12の出力は時間と共に変動する直流信号で
あり、それは火炎力に比例する。信号の時間変動部分は
コンデンサ19を介して直流成分から引き離され、第1
のアンプ20の出力信号は火炎力の変化の差に等しくな
るようにされる。この出力信号は一般的に“火炎フリッ
カ”と呼ばれる、火炎の形状もしくは強さの時間変化に
より生成される。火炎フリッカは火炎の存在及び安定度
等のいくつかの特性を決定するのに使用できる。
【0019】ノード48におけるセンサ信号の時間変動
部分はマイクロコンピュータ50のアナログ入力に加え
られ、アナログ信号の大きさを表わす10桁の10進数
へデジタル化される。マイクロコンピュータは規則的間
隔で割込みが行われ、アナログ/デジタルコンバータの
出力をサンプルしてデジタルサンプルをRAM56内の
リング型バッファに格納するためのソフトウェアルーチ
ンを実行する。例えば、マイクロコンピュータ50は毎
秒300の火炎信号サンプルを取得するように割込ま
れ、リングバッファは定期的に取り込まれるデータサン
プルが格納される600の記憶場所を有している。RA
M56内のもう一つの記憶場所には最近のデジタル数が
記憶されているリング型バッファの記憶場所のポインタ
が記憶されている。このポインタは後のデータ処理ステ
ップにおいて被処理データをリングのどこへ入れるかを
示す表示として使用される。
【0020】リング型バッファに256のデータサンプ
ルが格納されると、マイクロコンピュータ50は連続的
に背景分析タスクを実行開始する。図2Aを参照とし
て、分析の第1ステップ60においてRAM56に記憶
されたデータは時間領域から周波数領域へ変換される。
こうすることにより、従来の高速フーリエ変換ソフトウ
ェアルーチンを使用してデータサンプルに256点変換
を行って、0〜149Hzの火炎データのベクトル表示
周波数スペクトルを表わす、128の複素数が得られ
る。これらの複素数はRAM56内に一時記憶される。
変換が完了すると、プログラムの実行はステップ62へ
進みそこでマイクロコンピュータ50は各複素周波数ベ
クトルの大きさを計算する。これは複素数の実数部の二
乗プラス虚数部の二乗の平方根をとって行うことができ
る。各ベクトルの大きさはセンサ信号の周波数成分の振
幅を表わす。周波数領域への変換を行うのにフーリエ分
析が使用されるが、他の技術を使用することもできる。
ステップ64においてこれらの各振幅値に対する底eの
対数が計算されRAM56内のアレイに記憶される。次
に、ステップ66において対数値をデジタルローパス濾
波することによりそのデータが平滑化される。
【0021】図3は活性火炎により生じるセンサ信号に
対する対数周波数成分振幅値の分布を表わすグラフであ
る。0〜149Hzの周波数成分は火炎形状の時間変
化、すなわち火炎フリッカにより生じる。このグラフは
振幅が高い周波数に対して数学的に予測可能な、レシオ
メトリック関係で低減することを示している。対数振幅
値分布全体にわたって、一般的に線型関係が存在する。
【0022】被感知火炎が消炎してセンサ12が、多バ
ーナシステムの他の火炎等の、背景火炎源からの放射を
検出する場合には、周波数成分対数振幅値の分布は図4
のグラフに示すものと同様なものとなる。この場合、振
幅は周波数と共に低減するが、振幅値の対数と周波数と
の関係は線型とはならなくなる。従って、火炎が存在す
る場合と消炎する場合とでは周波数スペクトルデータの
数学的関係が違ってくる。
【0023】火炎の存在、種類及び状態は火炎信号の周
波数スペクトルに基づいてマイクロコンピュータ50に
より決定される。フーリエ変換された振幅値の対数が導
出されてRAM56に記憶されると、対数振幅の周波数
による変化(勾配)の連続性及び均一性が(例えば、0
〜100Hzの)所定帯域幅にわたってテストされる。
連続度及び均一度は量子化されて火炎安定度に比例する
数が得られる。スペクトルの線型性はいくつかの要因に
より変動する火炎信号振幅及び火炎力に対応する火炎信
号利得に無関係であることが判っている。簡単に言え
ば、火炎サイズはスペクトル周波数成分の形状分布に影
響を及ぼさない。従って、センサ12からの火炎信号の
振幅ではなくスペクトル線型性に依存する本システムに
より、ダンパ位置、燃料圧、噴霧圧、燃料充填率、空燃
比、BTU含有量、燃料タイプ、及び他の変数が分析に
及ぼす影響が著しく低減される。バーナ火炎がオンとさ
れて比較的安定しておれば、周波数成分振幅の勾配分布
は実質的に線型且つ均一になる。
【0024】図3のグラフの縦座標は対数であるため、
点火火炎のスペクトルは周波数領域内のセンサ信号を次
式で表わせることを示唆している。
【0025】
【数2】S(f)=Aekf (1)
【0026】直線に対して、この式のeを底とする対数
は次のようになり、
【0027】
【数3】 Log〔S(f)〕=kf+Log〔A〕 (2)
【0028】ここに、S(f)は周波数の関数としての
信号振幅、Aは0Hz(DC)における電圧振幅、eは
1の逆自然対数、kはデータの勾配(減衰定数)、fは
周波数Hzである。
【0029】消炎した火炎に対するセンサ信号周波数成
分分布を表わす図4は、周波数対対数振幅値の非線型す
なわち断片線型グラフを示している。このスペクトルデ
ータの粗適合により、センサ信号は次の形式を有するこ
とが示唆され、
【0030】
【数4】 S(f)=Aeaf+Bebf (3)
【数5】 Log〔S(f)〕=af+Log〔A〕+bf+Log〔B〕 (4) ここに、“a”は低周波データの勾配、“b”は高周波
データの勾配である。火炎オフ状態では、係数“a”は
低周波黒体すなわち対流放射から引き出され、係数
“b”は高周波ホワイトノイズすなわち他の隣接バーナ
火炎から引き出される。
【0031】図2Aの火炎分析プログラムフロー図の残
りの部分により、対数振幅値対周波数分布の線型度が決
定される。ステップ68で開始して、マイクロコンピュ
ータ50は最小二乗法及び周波数スペクトルデータを次
のような形式の3次多項式へ適合させる行列式を使用す
るルーチンを実行する。
【0032】
【数6】 s=a+bf+cf2 +df3 (5)
【0033】バーナ火炎がオンで安定しておれば、周波
数分布は線型に漸減して係数c及びdはゼロに近ずき、
直線式sL =a+bf.となる。しかしながら、火炎が
不安定となったり消えたりする場合、係数c及びdは一
層重要となる。これらの係数により次の3つのパラメー
タを決定することができ、積分線型誤差(E)、勾配差
(md )及び線型度回帰相関(R)、これらによって勾
配の連続度及び均一度が定まる。
【0034】第1のパラメータ、積分線型誤差、により
周波数成分振幅分布を一つの勾配係数により満足に記述
できるかどうかが決定される。このパラメータを決定す
るには、ステップ70において3次式(5)の導関数を
計算して比較的低い第1の周波数f1 (例えば、20H
z)における勾配が求められる。ステップ72において
この勾配を使用して第1の周波数におけるスペクトルグ
ラフの接線の公式が引き出され、次にこの公式を投影し
てそのY座標値が決定される。これにより、次式が得ら
れ、
【0035】
【数7】 st =at +bt 1 (6)
【0036】ここに、at は座標値、bt は接線の勾配
である。これを図5に示し、ここで実線は周波数成分振
幅値を表わし、破線は20Hzのf1 における周波数成
分振幅分布への接線を表わす。
【0037】全スペクトルにわたって周波数成分振幅値
の分布を一つの勾配係数により満足に描くことができる
(すなわち、スペクトルが線型である)場合には、全て
の周波数成分振幅値に対して(6)式が適合できなけれ
ばならない。積分線型誤差(E)、すなわち接線式の適
合度は周波数成分振幅値曲線と接線間の0〜100Hz
までの面積を求めることにより決定される。この面積は
次式で与えられる。
【0038】
【数8】
【0039】これを3次多項式で表わすと次のようにな
る。
【0040】
【数9】
【0041】積分線型誤差(E)はステップ74におい
て計算され、RAM56に格納される。
【0042】図2の分析プログラムの実行は次にステッ
プ76へ進み、そこで第2のパラメータ、勾配差、が周
波数スペクトル勾配連続性のもう一つの表示として計算
される。こうすることにより、第1の周波数点(20H
z)における勾配計算に使用したのと同じ方法で、第2
の周波数(例えば、80Hz)におけるスペクトルデー
タの勾配が計算される。2つのスペクトル勾配間の差が
計算されてRAM56に記憶される。
【0043】勾配の連続性及び均一性を示す最後のパラ
メータは周波数成分対数振幅分布に対する線型度回帰相
関である。2次元データアレイの所与数nの対データ点
(f 1 ,s1 ),(f2 ,s2 ),…(fn ,sn )に
対して、線型度回帰相関Rは次式で表わされ、
【0044】
【数10】 ここに、Si は周波数成分fi の対数振幅である。線型
度回帰相関が1であれば、線型度及び相関は完全であ
る。しかしながら、線型度回帰相関値が1よりも小さけ
れば、非線型で相関が悪いことを示す。ステップ78で
計算される相関値はRAM56に記憶される。
【0045】図2Bで開始される火炎分析背景プログラ
ムの残部は3つの火炎スペクトル形状パラメータ値を解
釈して火炎の存在及び安定性に関する決定に達する。背
景分析プログラムの連続ループにより1組の新しいスペ
クトル形状パラメータが毎秒数回計算される。各パラメ
ータに対する一連の値がRAM56内の別々のリング型
バッファ内に保持される。各バッファのサイズは、ユー
ザが設定する、失炎応答時間(FFRT)を指定した構
成パラメータにより決定される。FFRTは火炎がオン
であるかオフであるかを分析器10が決定するための最
大時間を定める。各火炎形状パラメータバッファ内のサ
ンプル記憶場所数はFFRTに単位時間当りの高速フー
リエ変速数を乗じた値に等しい。
【0046】各火炎スペクトル形状パラメータバッファ
内に記憶されたデータは平均化されて標準偏差が計算さ
れる。各算術平均及び標準偏差によりパラメータのガウ
ス分布曲線が定められる。一般的に“Lower−To
il Test”と呼ばれる統計法を適用して、そのパ
ラメータの臨界閾値よりも上及び下のデータの%が決定
される。火炎オン決定に達するのにデータサンプルの6
0%を閾値以上とし、且つ火炎オフの決定を行うのにデ
ータサンプルの60%を閾値以下とすることにより形状
パラメータ閾値に関するヒステリシスが与えられる。
【0047】図2Bの分析プログラムプログラムのフロ
ー図を参照として、ステップ80においてRAM56内
の対応するリングバッファ内に記憶されている積分線型
誤差に対する算術平均及び標準偏差がマイクロコンピュ
ータ50により計算される。次に、これらの統計値を使
用して、積分線型誤差に対する所定の閾値よりも上及び
下のバッファ値の%が決定される。ステップ84〜86
及び88〜90において同様な統計処理を行って、勾配
差に対するこのような%及びこれら別々の閾値に対する
線型度回帰相関値が引き出される。こうして求められる
%はRAM56に記憶される。
【0048】特定燃焼装置に対する分析器10の設定中
にパラメータ閾値を経験的に決定した。その時に、火炎
を点火させ一連の火炎スペクトル分析を実施した。この
一連の分析に対して、積分線型誤差、勾配差、線型度回
帰相関パラメータの最大値を求めた。次に、火炎を消炎
して、さらに一連の火炎スペクトル分析を実施した。次
に、3つの線型度パラメータの最大及び最小値を求め
た。各パラメータに対する最大及び最小値の中点がパラ
メータ閾値となる。
【0049】ステップ92における火炎分析プログラム
実行に戻って、3つのパラメータ値の閾値よりも上の%
を平均化し、3つのパラメータ値の閾値よりも下の%を
平均化した。“上”の%の平均をステップ94において
平均化して60%よりも大きいかどうかを決定し、大き
い場合にはプログラム実行をステップ98へ分岐してフ
ラグを設定し火炎がオンであることを表示する。次に、
プログラム実行をステップ60へ戻し、新しく取得した
データを使用して再度分析する。平均化された“上”の
%がステップ94において60%よりも大きくない場合
には、プログラムはステップ95へ進む。ここで“下”
の%の平均を調べて60%よりも大きいかどうかを決定
し、大きい場合にはプログラム実行はステップ96へ分
岐して火炎オンフラグをリセットしステップ60へ戻る
前に消炎したことを表示する。ステップ94及び95に
おいて行われるテストがいずれも真でない場合には、火
炎オンフラグ状態を変えることなくプログラムは直接ス
テップ60へ戻り、前に決定した状態がそのまま残され
る。
【0050】もう一つの背景ソフトウェアルーチンが火
炎オンフラグを定期的に調べ、フラグ状態を示す信号を
I/Oインターフェイス回路58を介して適切な外部装
置へ送る。
【図面の簡単な説明】
【図1】本発明を組み入れた火炎分析器の電子回路図。
【図2】火炎分析ソフトウェアのフロー図。
【図3】図1の火炎分析器により生成される安定バーナ
火炎に対する信号の周波数成分スペクトルを示す図。
【図4】被感知火炎が消炎する場合に背景放射により生
成される信号の、図3と類似の、スペクトルを示す図。
【図5】火炎信号分析における1ステップのグラフ表示
の図面。
【符号の説明】
10 バーナ火炎分析器 12 センサ 14 定利得プリアンプ 16 抵抗器 17 抵抗器 18 抵抗器 19 コンデンサ 20 アンプ 21 抵抗器 22 固定抵抗器 24 ホトレジスタ 25 発光ダイオード 26 抵抗器 28 アンプ 29 抵抗器 30 ダイオード 31 ダイオード 32 帰還抵抗器 33 抵抗器 34 アンプ 35 抵抗器 36 抵抗器 38 コンデンサ 39 抵抗器 40 差動アンプ 41 抵抗器 42 抵抗器 43 帰還抵抗器 46 抵抗器 48 出力ノード 50 マイクロコンピュータ 51 コントロールバス 52 コントロールバス 53 コントロールバス 54 ROM 56 RAM 58 インターフェイス回路
フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01J 1/42 F23M 11/04 102

Claims (21)

    (57)【特許請求の範囲】
  1. 【請求項1】 火炎により生じる放射を検出して放射を
    示す電気信号を発生するセンサと、電気信号を火炎力の
    時間変化により生じる電気信号の複数の周波数成分から
    なるスペクトルへ変換する手段を具備し、各周波数成分
    は振幅を有し、さらに、スペクトル全体にわたって周波
    数成分振幅分布の線型度を決定する手段と、線型度に応
    答して火炎の存否を決定する手段、を具備する火炎分析
    器。
  2. 【請求項2】 請求項記載の火炎分析器において、前
    記電気信号変換手段は電気信号に対してフーリエ変換を
    行う手段を具備する火炎分析器。
  3. 【請求項3】 請求項記載の火炎分析器において、前
    記線型度決定手段は少くともスペクトルの2つの周波数
    において周波数成分振幅分布の勾配差を決定する手段を
    具備する火炎分析器。
  4. 【請求項4】 請求項記載の火炎分析器において、前
    記線型度決定手段は次式で与えられる周波数成分振幅分
    布に対する線型度回帰相関(R)を決定し、 【数1】 ここに、iは周波数成分番号、Siは周波数成分f1
    の振幅である火炎分析器。
  5. 【請求項5】 請求項記載の火炎分析器において、前
    記線型度決定手段は、1組の周波数成分に対して各々が
    周波数成分振幅分布の所与の点における接線の公式によ
    り定まるデータ値を引き出す手段と、前記変換手段によ
    り生じる振幅と1組の周波数成分の各要素に対するデー
    タ値との差を積分して線型度を示す第1の値を生成する
    手段、を具備する火炎分析器。
  6. 【請求項6】 請求項記載の火炎分析器において、前
    記線型度決定手段はさらに、スペクトルの線型性を示す
    第2の値となる周波数成分振幅分布の2つの位置におけ
    る勾配差を決定する手段と、スペクトルの線型性を示す
    第3の値となる周波数成分振幅分布に対する線型度回帰
    相関を決定する手段、を具備する火炎分析器。
  7. 【請求項7】 請求項記載の火炎分析器において、前
    記火炎特性決定手段は、複数の第1の値を第1の閾値と
    比較して第1の閾値よりもそれぞれ上及び下である第1
    の値の量を決定する第1の手段と、複数の第2の値を第
    2の閾値と比較して第2の閾値よりもそれぞれ上及び下
    である量を決定する第2の手段と、複数の第3値を第3
    の閾値と比較して第3の閾値よりもそれぞれ上及び下で
    ある量を決定する第3の手段と、それぞれの閾値よりも
    下の第1、第2及び第3の値の量を平均化して第1の平
    均を生成する第1の手段と、それぞれの閾値よりも上の
    第1、第2及び第3の値の量を平均化して第2の平均を
    生成する第2の手段と、第1の平均が第1の基準値を越
    える時に火炎の消炎表示を発生し第2の平均が第2の基
    準値を越える時に火炎の存在表示を発生する手段、を具
    備する火炎分析器。
  8. 【請求項8】 火炎の存在を検出する装置において、該
    装置は、火炎により生じる放射を検出して放射を示す電
    気信号を発生するセンサと、電気信号を増幅する自動利
    得制御増幅器と、前記増幅器からの電気信号を複数の信
    号サンプルへデジタル化する手段と、複数の信号サンプ
    ルを記憶する手段と、時間領域からの信号サンプルを周
    波数領域へ変換して複数の周波数成分振幅値を生成する
    手段と、電気信号の各周波数成分振幅値の対数値を引き
    出す手段と、対数値分布の線型度を決定する手段と、線
    型度を評価して火炎が存在しているかどうかを決定する
    手段、を具備する火炎の存在検出装置。
  9. 【請求項9】 請求項記載の火炎分析器において、前
    記線型度決定手段は対数値分布に沿った2つの位置にお
    ける勾配差を決定する、火炎分析器。
  10. 【請求項10】 請求項記載の火炎分析器において、
    前記線型度決定手段は対数値分布に対する線型回帰相関
    を決定する、火炎分析器。
  11. 【請求項11】 請求項記載の火炎分析器において、
    前記線型度決定手段は、対数値分布に沿った所与の点に
    おける接線を定義する手段と、対数値と定義された線上
    の点との一連の差を積分して対数値の線型性を示す第1
    の値を生成する手段、を具備する火炎分析器。
  12. 【請求項12】 請求項11記載の火炎分析器におい
    て、前記線型度決定手段は、対数値分布上の2点間の勾
    配差を決定して線型度を示す第2の値を生成する手段
    と、対数値分布に対する線型回帰相関を決定して線型度
    を示す第3の値を生成する手段、を具備する火炎分析
    器。
  13. 【請求項13】 請求項12記載の火炎分析器におい
    て、線型度評価手段は、複数の第1の値を第1の閾値と
    比較して第1の閾値よりも下の第1の値の量及び第1の
    閾値よりも上の第1の値の量を決定する第1の手段と、
    複数の第2の値を第2の閾値と比較して第2の閾値より
    も下の第2の値の量及び第2の閾値よりも上の第2の値
    の量を決定する第2の手段と、複数の第3の値を第3の
    閾値と比較して第3の閾値よりも下の第3の値の量及び
    第3の閾値よりも上の第3の値の量を決定する第3の手
    段と、それぞれの閾値よりも下の第1、第2及び第3の
    値の量を平均化して第1の平均を生成する第1の手段
    と、それぞれの閾値よりも上の第1、第2及び第3の値
    の量を平均化して第2の平均を生成する第2の手段と、
    第1の平均が第1の基準値を越える場合に火炎の消炎表
    示を発生し第2の平均が第2の基準値を越える場合に火
    炎の存在表示を発生する手段、を具備する火炎分析器。
  14. 【請求項14】 火炎特性決定方法において、該方法
    は、火炎により生じる周波数における放射を検出して、
    放射を表示する電気信号を発生し、電気信号を時間領域
    から周波数領域へ変換して火炎形状の時間変化により生
    じる複数の周波数成分に対する振幅値を生成し、振幅値
    分布の線型度を決定し、線型度を使用して火炎が存在す
    るか消えているかを決定する、ことからなる火炎特性決
    定方法。
  15. 【請求項15】 請求項14記載の方法において、前記
    電気信号変換ステップは電気信号のフーリエ変換を実施
    することからなる火炎特性決定方法。
  16. 【請求項16】 請求項14記載の方法において、前記
    線型度決定ステップは各振幅値に対する対数値を引き出
    し、対数値分布の線型度を決定することからなる、火炎
    特性決定方法。
  17. 【請求項17】 請求項14記載の方法において、前記
    線型度決定ステップは振幅値分布上の2点における勾配
    差を引き出すことからなる、火炎特性決定方法。
  18. 【請求項18】 請求項14記載の方法において、前記
    線型度決定ステップは振幅値分布に対する線型回帰相関
    を決定する手段を具備する、火炎特性決定方法。
  19. 【請求項19】 請求項14記載の方法において、前記
    線型度決定ステップは、振幅値分布の所与の点における
    接線の公式から周波数成分のデータ値を引き出し、振幅
    値と所与の周波数帯域内の各周波数成分に対するデータ
    値との差を積分して線型度を示す第1の値を生成する、
    ことからなる火炎特性決定方法。
  20. 【請求項20】 請求項19記載の方法において、前記
    線型度決定ステップは、振幅値分布上の2点間の勾配差
    を決定して線型度を表わす第2の値を生成し、振幅値分
    布に対する線型回帰相関を決定して線型度を示す第3の
    値を生成する、ことからなる火炎回帰決定方法。
  21. 【請求項21】 請求項20記載の方法において、前記
    線型度決定ステップは、複数の第1の値を第1の閾値と
    比較して第1の閾値よりも下の第1の値の量及び第1の
    値よりも上の第1の値の量を決定し、複数の第2の値を
    第2の閾値と比較して第2の閾値よりも下の第2の値の
    量及び第2の閾値よりも上の第2の値の量を決定し、複
    数の第3の値を第3の閾値と比較して第3の閾値よりも
    下の第3の値の量及び第3の閾値よりも上の第3の値の
    量を決定し、それぞれの閾値よりも下の第1、第2及び
    第3の値の量を平均化して第1の平均を生成し、それぞ
    れの閾値よりも上の第1、第2及び第3の閾値の量を平
    均化して第2の平均を生成し、前記線型度を使用して火
    炎特性を決定するステップにより第1の平均が第1の基
    準値を越える場合に火炎の消炎表示が発生され、第2の
    平均が第2の基準を越える場合に火炎の存在表示が発生
    される、火炎特性決定方法。
JP03239682A 1990-09-19 1991-09-19 火炎分析器及び火炎特性決定方法 Expired - Lifetime JP3143500B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US584690 1990-09-19
US07/584,690 US5077550A (en) 1990-09-19 1990-09-19 Burner flame sensing system and method

Publications (2)

Publication Number Publication Date
JPH04244922A JPH04244922A (ja) 1992-09-01
JP3143500B2 true JP3143500B2 (ja) 2001-03-07

Family

ID=24338421

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03239682A Expired - Lifetime JP3143500B2 (ja) 1990-09-19 1991-09-19 火炎分析器及び火炎特性決定方法

Country Status (7)

Country Link
US (1) US5077550A (ja)
EP (1) EP0476601A1 (ja)
JP (1) JP3143500B2 (ja)
AR (1) AR245821A1 (ja)
BR (1) BR9103943A (ja)
CA (1) CA2050805C (ja)
MX (1) MX9101134A (ja)

Families Citing this family (58)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9019457D0 (en) * 1990-09-06 1990-10-24 Dresser Holmes Limited Flame monitoring apparatus and method
US5365223A (en) * 1991-10-28 1994-11-15 Honeywell Inc. Fail-safe condition sensing circuit
US5236328A (en) * 1992-09-21 1993-08-17 Honeywell Inc. Optical flame detector performance tester
US5434560A (en) * 1993-05-11 1995-07-18 Detector Electronics Corporation System for detecting random events
US5549469A (en) * 1994-02-28 1996-08-27 Eclipse Combustion, Inc. Multiple burner control system
US5798946A (en) * 1995-12-27 1998-08-25 Forney Corporation Signal processing system for combustion diagnostics
US5864714A (en) * 1996-04-16 1999-01-26 Comsys Communication & Signal Processing Ltd. Communication system which dynamically switches sizes of sample buffer between first size for quick response time and second size for robustness to interrupt latency
US6045353A (en) 1996-05-29 2000-04-04 American Air Liquide, Inc. Method and apparatus for optical flame control of combustion burners
US5829962A (en) * 1996-05-29 1998-11-03 L'air Liquide, Societe Anonyme Pour L'etude Et, L'exploitation Des Procedes Georges Method and apparatus for optical flame control of combustion burners
US6389330B1 (en) 1997-12-18 2002-05-14 Reuter-Stokes, Inc. Combustion diagnostics method and system
AU768582B2 (en) * 1998-06-02 2003-12-18 Hochiki Kabushiki Kaisha Flame detection device and flame detection method
US6341519B1 (en) 1998-11-06 2002-01-29 Reuter-Stokes, Inc. Gas-sensing probe for use in a combustor
US6277268B1 (en) 1998-11-06 2001-08-21 Reuter-Stokes, Inc. System and method for monitoring gaseous combustibles in fossil combustors
GB2344883B (en) * 1998-12-16 2003-10-29 Graviner Ltd Kidde Flame monitoring methods and apparatus
US6127932A (en) * 1998-12-23 2000-10-03 Carrier Corporation Optical flame sensor having opaque hollow tube
US6024561A (en) * 1999-01-20 2000-02-15 Autoflame Engineering Limited Monitoring for the presence of a flame in a burner
US6478573B1 (en) * 1999-11-23 2002-11-12 Honeywell International Inc. Electronic detecting of flame loss by sensing power output from thermopile
US6261086B1 (en) 2000-05-05 2001-07-17 Forney Corporation Flame detector based on real-time high-order statistics
US6702571B2 (en) 2001-09-05 2004-03-09 Gas Technology Institute Flex-flame burner and self-optimizing combustion system
US7128818B2 (en) * 2002-01-09 2006-10-31 General Electric Company Method and apparatus for monitoring gases in a combustion system
US7255285B2 (en) * 2003-10-31 2007-08-14 Honeywell International Inc. Blocked flue detection methods and systems
CN1894543A (zh) * 2003-12-11 2007-01-10 Abb公司 用于改进的火焰扫描器鉴别的信号处理技术
WO2005111556A2 (en) * 2004-05-07 2005-11-24 Walter Kidde Portable Equipment, Inc. Flame detector with uv sensor
US8085521B2 (en) 2007-07-03 2011-12-27 Honeywell International Inc. Flame rod drive signal generator and system
US7764182B2 (en) * 2005-05-12 2010-07-27 Honeywell International Inc. Flame sensing system
US7800508B2 (en) * 2005-05-12 2010-09-21 Honeywell International Inc. Dynamic DC biasing and leakage compensation
US8066508B2 (en) * 2005-05-12 2011-11-29 Honeywell International Inc. Adaptive spark ignition and flame sensing signal generation system
US7768410B2 (en) * 2005-05-12 2010-08-03 Honeywell International Inc. Leakage detection and compensation system
US8310801B2 (en) * 2005-05-12 2012-11-13 Honeywell International, Inc. Flame sensing voltage dependent on application
US8300381B2 (en) 2007-07-03 2012-10-30 Honeywell International Inc. Low cost high speed spark voltage and flame drive signal generator
US8875557B2 (en) * 2006-02-15 2014-11-04 Honeywell International Inc. Circuit diagnostics from flame sensing AC component
US7806682B2 (en) * 2006-02-20 2010-10-05 Honeywell International Inc. Low contamination rate flame detection arrangement
US7728736B2 (en) * 2007-04-27 2010-06-01 Honeywell International Inc. Combustion instability detection
US8926317B2 (en) * 2008-12-15 2015-01-06 Exxonmobil Research And Engineering Company System and method for controlling fired heater operations
US9366433B2 (en) * 2010-09-16 2016-06-14 Emerson Electric Co. Control for monitoring flame integrity in a heating appliance
US8457835B2 (en) * 2011-04-08 2013-06-04 General Electric Company System and method for use in evaluating an operation of a combustion machine
US9494320B2 (en) 2013-01-11 2016-11-15 Honeywell International Inc. Method and system for starting an intermittent flame-powered pilot combustion system
US10208954B2 (en) 2013-01-11 2019-02-19 Ademco Inc. Method and system for controlling an ignition sequence for an intermittent flame-powered pilot combustion system
US20140202549A1 (en) 2013-01-23 2014-07-24 Honeywell International Inc. Multi-tank water heater systems
US10670302B2 (en) 2014-03-25 2020-06-02 Ademco Inc. Pilot light control for an appliance
US20150277463A1 (en) 2014-03-25 2015-10-01 Honeywell International Inc. System for communication, optimization and demand control for an appliance
US10288286B2 (en) 2014-09-30 2019-05-14 Honeywell International Inc. Modular flame amplifier system with remote sensing
US10042375B2 (en) 2014-09-30 2018-08-07 Honeywell International Inc. Universal opto-coupled voltage system
US10402358B2 (en) 2014-09-30 2019-09-03 Honeywell International Inc. Module auto addressing in platform bus
US10678204B2 (en) 2014-09-30 2020-06-09 Honeywell International Inc. Universal analog cell for connecting the inputs and outputs of devices
US9799201B2 (en) 2015-03-05 2017-10-24 Honeywell International Inc. Water heater leak detection system
US9920930B2 (en) 2015-04-17 2018-03-20 Honeywell International Inc. Thermopile assembly with heat sink
US10132510B2 (en) 2015-12-09 2018-11-20 Honeywell International Inc. System and approach for water heater comfort and efficiency improvement
US10119726B2 (en) 2016-10-06 2018-11-06 Honeywell International Inc. Water heater status monitoring system
US10473329B2 (en) 2017-12-22 2019-11-12 Honeywell International Inc. Flame sense circuit with variable bias
US11236930B2 (en) 2018-05-01 2022-02-01 Ademco Inc. Method and system for controlling an intermittent pilot water heater system
US10935237B2 (en) 2018-12-28 2021-03-02 Honeywell International Inc. Leakage detection in a flame sense circuit
US10969143B2 (en) 2019-06-06 2021-04-06 Ademco Inc. Method for detecting a non-closing water heater main gas valve
US11656000B2 (en) 2019-08-14 2023-05-23 Ademco Inc. Burner control system
US11739982B2 (en) 2019-08-14 2023-08-29 Ademco Inc. Control system for an intermittent pilot water heater
JP7690290B2 (ja) * 2021-01-13 2025-06-10 関西電力株式会社 輝度判定システム
CN115035677B (zh) * 2022-06-17 2023-01-20 新创碳谷集团有限公司 一种氧化炉炉内火灾火焰的探测系统及灭火的安保系统
TWI869088B (zh) * 2023-11-30 2025-01-01 財團法人金屬工業研究發展中心 燃燒器火焰量測系統及其方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866420A (en) 1988-04-26 1989-09-12 Systron Donner Corp. Method of detecting a fire of open uncontrolled flames

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2748846A (en) * 1952-01-25 1956-06-05 Honeywell Regulator Co Combustion safeguard apparatus
US2762033A (en) * 1955-09-23 1956-09-04 Electronics Corp America Fire detectors
DE1960218A1 (de) * 1969-12-01 1971-06-03 Rainer Portscht Temperaturstrahlungsdetektor zur automatischen Brandentdeckung oder Flammenueberwachung
US3689773A (en) * 1971-02-01 1972-09-05 Bailey Miters & Controls Ltd Flame monitor system and method using multiple radiation sensors
CH537066A (de) * 1971-04-08 1973-05-15 Cerberus Ag Flammen-Detektor
IT1011594B (it) * 1973-03-20 1977-02-10 Electricite De France Dispositivo per il controllo della fiamma di bruciatori
CH558577A (de) * 1973-09-25 1975-01-31 Cerberus Ag Verfahren zur flammen-detektion und vorrichtung zur durchfuehrung dieses verfahrens.
JPS5435426A (en) * 1977-08-24 1979-03-15 Showa Yuka Kk Apparatus for monitoring flame from flare stack
US4410266A (en) * 1980-08-25 1983-10-18 Bsc Industries Corp. Method and apparatus for combustion control and improved optical pyrometer related thereto
US4370557A (en) * 1980-08-27 1983-01-25 Honeywell Inc. Dual detector flame sensor
US4769775A (en) * 1981-05-21 1988-09-06 Santa Barbara Research Center Microprocessor-controlled fire sensor
JPS59101830A (ja) * 1982-12-01 1984-06-12 Canon Inc 転写装置
US4553031A (en) * 1983-09-06 1985-11-12 Firetek Corporation Optical fire or explosion detection system and method
US4691196A (en) * 1984-03-23 1987-09-01 Santa Barbara Research Center Dual spectrum frequency responding fire sensor
US4665390A (en) * 1985-08-22 1987-05-12 Hughes Aircraft Company Fire sensor statistical discriminator
US4709155A (en) * 1984-11-22 1987-11-24 Babcock-Hitachi Kabushiki Kaisha Flame detector for use with a burner
US4701624A (en) * 1985-10-31 1987-10-20 Santa Barbara Research Center Fire sensor system utilizing optical fibers for remote sensing

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866420A (en) 1988-04-26 1989-09-12 Systron Donner Corp. Method of detecting a fire of open uncontrolled flames

Also Published As

Publication number Publication date
JPH04244922A (ja) 1992-09-01
BR9103943A (pt) 1992-05-26
EP0476601A1 (en) 1992-03-25
CA2050805C (en) 1998-02-10
AR245821A1 (es) 1994-02-28
US5077550A (en) 1991-12-31
MX9101134A (es) 1992-05-04
CA2050805A1 (en) 1992-03-20

Similar Documents

Publication Publication Date Title
JP3143500B2 (ja) 火炎分析器及び火炎特性決定方法
US20080266120A1 (en) Combustion instability detection
AU675755B2 (en) Flame detector self diagnostic system
WO1997024560A1 (en) Signal processing system for combustion diagnostics
KR860008509A (ko) 패턴 인지 자동 조정 제어기
US5755819A (en) Photodiode array for analysis of multi-burner gas combustors
EP2431663B1 (en) Control for monitoring flame integrity in a heating appliance
US6700495B2 (en) Flame monitor for an oil- and gas-operated burner
JP3357132B2 (ja) 被検物質の濃度測定のための分光測光法およびこの方法を実施するための分光光度計
CN1099660C (zh) 探测火焰的方法及实施该方法的火焰报警器
JP3938276B2 (ja) 炎感知器および炎検知方法
US7076403B2 (en) Apparatus and method for dynamic smoothing
JP2002525544A (ja) 燃焼室の煤負荷を求める方法および装置
JPH02302520A (ja) バーナの燃焼制御装置
US20070207423A1 (en) Preemptive tripping of a fossil fuel fired burner using the quality function
JP2540990B2 (ja) バ―ナの燃焼制御装置
CN120823680B (zh) 一种基于多波段红外信号融合的火焰识别方法及火焰探测器
JPH07117238B2 (ja) バーナの燃焼制御方法
JPH0472126B2 (ja)
KR20240106796A (ko) 보일러의 배기 폐쇄 검출 장치 및 방법
JP2622382B2 (ja) 火炎検出器
US20140170574A1 (en) Flame instability detector
JPS6082933A (ja) 火炎検出装置
JPH01142419A (ja) 火炎検知方法
JP3814047B2 (ja) 燃焼状況検知方法

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071222

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081222

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091222

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101222

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101222

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111222

Year of fee payment: 11

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111222

Year of fee payment: 11