JP3804129B2 - Densitometer - Google Patents

Densitometer Download PDF

Info

Publication number
JP3804129B2
JP3804129B2 JP30019896A JP30019896A JP3804129B2 JP 3804129 B2 JP3804129 B2 JP 3804129B2 JP 30019896 A JP30019896 A JP 30019896A JP 30019896 A JP30019896 A JP 30019896A JP 3804129 B2 JP3804129 B2 JP 3804129B2
Authority
JP
Japan
Prior art keywords
signal
microwave
frequency
phase difference
oscillator
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 - Fee Related
Application number
JP30019896A
Other languages
Japanese (ja)
Other versions
JPH10142168A (en
Inventor
照彦 矢上
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Tokyo Keiki Inc
Original Assignee
Toshiba Corp
Tokyo Keiki Inc
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 Toshiba Corp, Tokyo Keiki Inc filed Critical Toshiba Corp
Priority to JP30019896A priority Critical patent/JP3804129B2/en
Publication of JPH10142168A publication Critical patent/JPH10142168A/en
Application granted granted Critical
Publication of JP3804129B2 publication Critical patent/JP3804129B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
  • Measuring Phase Differences (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、マイクロ波を用いて被測定流体の濃度を測定する濃度計に関し、特にその精度向上をはかることのできる濃度計に関する。
【0002】
【従来の技術】
従来、この種の濃度計としては、図9に示されたものがある。この濃度計は、配管を流れる被測定流体にマイクロ波を入射し、被測定物質を伝搬してきた受信波の位相遅れを測定し、さらには、基準位相遅れ(同一配管に例えば水を流したときの位相遅れ)との差から被測定物質の濃度を測定するものである。
【0003】
具体的な装置構成としては、検出用の配管10に、互いに対向するようにして、送信アンテナ12と受信アンテナ14が取り付けられている。高周波発振器60から出力される例えば1.6GHzのマイクロ波信号が送信アンテナ12に送られて、送信アンテナ12から発射されたマイクロ波が配管10内の被測定流体の中を伝搬し、受信アンテナ14で受信される。受信アンテナ14からの受信信号はアンプ61を通った後ミキサ62に送られ、同時に、高周波発振器60からのマイクロ波信号は、別のミキサ64に送られて、両ミキサ62、64において、1.6GHz−60MHzの周波数を持つマイクロ波信号を発振する高周波発振器66からの出力と混合され、それぞれバンドパスフィルタ68、70を通ってそのビート信号が抽出され、アンプ72、74、リミッタ76、78等を介して直交同期検波器80に送られる。直交同期検波器80で2つの入力の位相差Δθの正弦成分sinΔθと余弦成分cosΔθが出力され、それぞれA/D変換器82、84で、A/D変換されて位相差演算回路86にて、
【0004】
【数1】
Δθ=tan-1(sinθ/cosθ)
が演算されて、位相差、即ち、被測定流体を伝搬したマイクロ波の位相遅れが求められる。
【0005】
【発明が解決しようとする課題】
かかる従来の濃度計では、特に、直交同期検波器80において電圧オフセット誤差が発生するために、その誤差を抑えるために、リミッタ76、78を用いて、電圧を一定に抑える必要があり、回路が複雑になり、コストも高くなるという問題がある。
【0006】
また、2つの高周波発振器60,66を用いているので、これらの出力が個別にゆらぐとこれにより、誤差が生じるという問題がある。
さらに、上記誤差以外にも、アナログ回路を多数用いているので、各アナログ回路の持つ誤差要因が重なり、結果として得られる位相差には、電圧オフセット誤差、位相リニアリティ誤差、位相オフセット誤差、周波数オフセット誤差等が含まれているため、これらの誤差を校正するために、予め、電力、位相、周波数をパラメータとして、これらを変化させたときの誤差を求めて、補正カーブを作成しておき、測定時に、この補正カーブを用いて出力された結果の補正を行っている。
【0007】
しかしながら、かかる補正カーブを作成するには時間と手間がかかる上に、アナログ回路毎にその特性が異なるために、各濃度計ごとに補正カーブを作成しなければならないという問題がある。また、補正カーブを用意しても、アナログ回路の経年変化や温度変化によって誤差特性が変化してしまうという問題もある。温度変化に関しては、装置をヒータで一定温度に加熱した容器の中に収納することで対応することもできるが、ヒータ等を必要とするために、コストがかかるという問題もある。
【0008】
本発明は、かかる問題点に鑑みなされたもので、請求項1ないし請求項3の記載の発明は、補正データを予め用意することなく、精度良く、濃度を測定することができる濃度計を提供することをその目的とする。
【0009】
【課題を解決するための手段】
上記目的を達成するために本発明のうち請求項1記載の発明は、被測定流体を伝搬したマイクロ波と被測定流体を伝搬しない信号との位相差を検出することにより、被測定流体の濃度を測定する濃度計において、基準周波数の信号を発振する基準発振器と、該基準周波数からフェーズロックされた周波数f1のマイクロ波信号を出力する第1マイクロ波発振器と、該基準周波数からフェーズロックされた前記マイクロ波信号と異なる周波数f2のマイクロ波信号を出力する第2マイクロ波発振器と、前記第1マイクロ波発振器に接続され、マイクロ波を被測定流体に放射する送信器と、被測定流体を伝搬したマイクロ波を受波する受信器と、前記受信器からのマイクロ波信号を前記第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換する周波数変換器と、前記周波数変換器からの低周波信号を矩形波に変換する波形整形回路と、被測定流体を伝搬しない信号と、前記波形整形回路からの矩形波信号とが入力され、両信号の位相差を、前記基準発振器からの基準周波数の信号を用いて計数する位相差演算回路と、を備えることを特徴とする。
【0010】
また、請求項2記載の発明は、請求項1記載の濃度計において、前記第1マイクロ波発振器からのマイクロ波信号を、前記第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換する第2周波数変換器と、前記第2周波数変換器からの低周波信号を矩形波信号に変換する第2波形整形回路と、を備え、該第2波形整形回路からの矩形波信号を前記被測定流体を伝搬しない信号として、前記位相差演算回路に入力することを特徴とする。
【0011】
また、請求項3記載の発明は、請求項1記載の濃度計において、基準発振器からの基準周波数の信号を|f1−f2|に分周した信号を、前記被測定流体を伝搬しない信号として、前記位相差演算回路に入力することを特徴とする。
【0012】
【発明の実施の形態】
以下、図面を用いて本発明の実施の形態について説明する。図1は、本発明の第1の実施の形態を示すブロック図であり、図において、被測定流体が流れる検出用の配管10に、互いに対向するようにして、送信アンテナ12と受信アンテナ14が取り付けられている。尚、送信アンテナ12と受信アンテナ14が取り付けられる配管10の部分は絶縁体となっている。
【0013】
水晶発振器等の基準発振器18からは分周器20を介してPLLシンセサイザ22、PLLシンセサイザ24が接続されており、シンセサイザ22は後述のスイッチ26を介して、送信アンテナ12に接続されている。基準発振器18から出力される、例えば60MHzの信号は、分周器20で1/6に分周された後、シンセサイザ22でフェイズロックされた1.6GHzのマイクロ波信号にされて、送信アンテナ12に送られ、送信アンテナ12からマイクロ波が発射される。また、シンセサイザ24では、シンセサイザ22と非常に近い、1.6GHz−333kHzのフェーズロックされた信号が出力される。
【0014】
受信アンテナ14からの受信信号は、スイッチ30を介してミキサ(周波数変換器)32に入力されて、そこで、シンセサイザ24からの信号と混合され、バンドパスフィルタ36を通って333kHzの低周波信号に変換される。また同時に、シンセサイザ22からの信号は、直接、ミキサ(周波数変換器)34に入力されて、シンセサイザ24からの信号と混合され、バンドパスフィルタ38を通って333kHzの低周波信号に変換される。これら低周波信号に変換された信号は、それぞれアンプ40、42及びコンパレータ44、46、位相差演算回路48に接続される。コンパレータ44、46は、アンプ40、42から送られてくる333kHzの正弦波信号のゼロクロスを検出するコンパレータであり、これにより正弦波を矩形波に整形して、位相差演算回路48に送るものである。
【0015】
位相差演算回路48は、ASICで構成されるデジタルシグナルプロセッサであり、基準発振器18からの信号をクロック信号とし、コンパレータ44及びコンパレータ46から入力されてくるデジタル波形の時間差を計数することによって、シンセサイザ22から被測定流体を透過してきた信号と、シンセサイザ22から被測定流体を透過せずに直接送られてきた信号の位相差を検出するものである。検出された位相差は、中央演算回路50へさらに送られる。
【0016】
本実施の形態では、シンセサイザ22と送信アンテナ12間、受信アンテナ14とミキサ32間に、それぞれ、スイッチ26、スイッチ30が介挿されており、スイッチ26の一つの接点S1は送信アンテナ12に接続され、スイッチ30の一つの接点S2は受信アンテナ14に接続され、これらS1−S2の間にある送信アンテナ12、受信アンテナ14を含むシンセサイザ22から位相差演算回路48までのラインは、測定ラインSLとなっている。また、スイッチ26の他方の接点R1と、スイッチ30の他方の接点R2をつなぐ、シンセサイザ22から位相差演算回路48までの間は、被測定流体を通過しない、参照ラインRLを形成しており、この参照ラインRLには、被測定流体を伝搬してきた場合と同程度の減衰を与える減衰器28が設けられている。スイッチ26、スイッチ30の切換は、中央演算回路50からの制御信号に基づいて高速に行われ、測定ラインSLまたは参照ラインRLのいずれかにシンセサイザ22からのマイクロ波信号が選択して、送出される。
【0017】
以上のように構成された濃度計において、まず、スイッチ26、30によりシンセサイザ22からのマイクロ波信号が測定ラインSLに伝搬するように切り換えられている状態では、1.6GHzのマイクロ波信号が送信アンテナ12に送られて、送信アンテナから発射されたマイクロ波が配管10内の被測定流体の中を伝搬し、受信アンテナ14で受信されて、ミキサ32、バンドパスフィルタ36、アンプ40、コンパレータ44を経て位相差演算回路48に測定信号Sとして入力される。位相差演算回路48には、同時に、シンセサイザ22から被測定流体を通らずに、ミキサ34、バンドパスフィルタ38、アンプ42、コンパレータ46を経た基準信号RFも入力されて、基準信号RFからの測定信号Sの位相遅れ(位相差)が求められる。即ち、位相差は、
【0018】
【数2】
Δθ1=θs−θ0+θe1
で表される。ここで、Δθ1は位相差演算回路48で求められた位相差であり、θsは測定信号の位相、θ0は基準信号の位相、θe1は、測定系のドリフト、オフセット誤差、温度変化による誤差、アナログ回路のばらつき等による誤差成分である。
【0019】
次に、スイッチ26、30を切り換えて、シンセサイザ22からのマイクロ波信号が参照ラインRLを伝搬するように切り換えると、マイクロ波信号が減衰器28を通過し、ミキサ32、バンドパスフィルタ36、アンプ40、コンパレータ44を経て位相差演算回路48に参照信号Rとして入力される。同時に基準信号RFも位相差演算回路48に入力されるのは同じである。このとき位相差演算回路48で求められる位相遅れ(位相差)は、
【0020】
【数3】
Δθ2=θr−θ0+θe2
で表される。ここで、Δθ2は位相差演算回路48で求められた位相遅れであり、θrは参照ラインRLを通過してきた参照信号Rの位相、θ0は基準信号の位相、θe2は、測定系のドリフト、オフセット誤差、温度変化による誤差、アナログ回路のばらつきによる誤差成分である。
【0021】
ここで、測定ラインSLを通った測定信号Sと参照ラインRLを通った参照信号Rは、時分割によって交互に得られるものの、測定ラインSLと参照ラインRLとが同じアナログ回路を共有しており、測定系のドリフト、オフセット誤差、温度による誤差、アナログ回路のばらつきによる誤差に関して対等であると考えられる。即ち、
【0022】
【数4】
θe1 =θe2
が成り立つと考えられる。従って、Δθ1とΔθ2との差を求めると、
【0023】
【数5】

Figure 0003804129
となる。θrは、参照ラインRLの電気長から予め既知であるため、θsに相当する値を得ることができる。さらに、同一配管10に水を流してそのときのθsに相当する値を求めることにより、両者の差から、被測定流体の濃度が測定される。
【0024】
図2は位相差演算回路48に入力される信号の関係を表しており、図2(a)は基準信号RF、測定信号S及びクロック信号CLとの関係を表し、図2(b)は基準信号RF、参照信号R及びクロック信号CLとの関係を表す。基準信号RFと測定信号Sとの位相差は、2つの信号の立ち上がり間に存在するクロック数(taとする)を計数すると共に、1つの信号の立ち上がりと立ち上がりとの間に存在するクロック数(tcとする)を計数することで、Δθ1が求められる。即ち、
【0025】
【数6】
Δθ1=2π・ta/tc
である。また、同様に、基準信号RFと参照信号Rとの位相差は、2つの信号の立ち上がり間に存在するクロック数(tdとする)を計数すると共に、1つの信号の立ち上がりと立ち上がりとの間に存在するクロック数(tcとする)を計数することで、Δθ2が求められる。即ち、
【0026】
【数7】
Δθ2=2π・td/tc
である。
【0027】
図2(c)、(d)には、さらに、スイッチ26、30の漏れがあるときの図2(a)、(b)対応図を示す。スイッチ26、30で測定ラインSLと参照ラインRLを切り換えているが、スイッチ26、30のアイソレーションが不十分であると、他方のラインからの漏れが干渉して、ベクトル合成される。図3(a)に信号をベクトル表示した説明図を示すと、今測定ラインSLに接続されて、測定信号S1が得られているとすると、参照ラインRLからの漏れ信号Eが微小ながら加算されて、位相がθe3だけ大きくなり誤差となる。測定信号S1から位相が180度ずれたところS2では、参照ラインRLからの漏れ信号Eにより、位相がθe3だけ小さくなる。この位相誤差θe3の位相依存性を示すと図3(b)のようになり、180度位相が異なるところでは、位相誤差θe3は絶対値が同じで符号が反転することがわかる。従って、このことを利用して、図2(c)において、基準信号RFの立ち上がりと測定信号Sの立ち上がりの間に存在するクロック数(ta)と、基準信号RFの立ち上がりと測定信号Sの立ち下がりの間に存在するクロック数(tb)とをそれぞれ求め、その平均をとる。
【0028】
【数8】
Δθ3=2π・(ta+tb)/2tc
同様に、参照信号Rについても行う。
【0029】
【数9】
Δθ4=2π・(td+te)/2tc
Δθ3とΔθ4の差をとることでθs−θrを求めることができる。言い換えれば、測定信号S及び参照信号Rのぞれぞれの位相90度の時点を基準にして、位相差θs−θrを求めることにより、スイッチ26、30の漏れによる誤差の影響を受けないようにすることができる。
【0030】
以上のように、シンセサイザ22とシンセサイザ24の出力する周波数を非常に近いものとし(333kHz)、その差をデジタルシグナルプロセッサで処理できる程度の周波数となるように設定することにより、位相差演算回路48をデジタルシグナルプロセッサとすることができるので、従来のアナログ回路を多数用いていたものに比べて、アナログ回路の特性に起因する誤差が少なくなり、精度を向上させることができる。
【0031】
また、シンセサイザ22、24からのマイクロ波信号は基準発振器18からの信号にフェーズロックされ、且つ位相差演算回路48は、基準発振器18からの信号を用いて位相差を検出することから、例えば、基準発振器18の信号がゆらいでも、全てが同一の影響を受けるため、結果として位相差演算回路48から出力される位相差は、基準発振器18のゆらぎの影響を受けずにすむ。
【0032】
図6及び図7は、それぞれθsとθrの時間変化を表したもので、0.18度から0.19度の範囲での変動が観測されるのに対して、図8に示したように、その差をとることにより、0.06度の範囲内の変動に抑えることができた。
図4は、本発明の他の実施の形態を表すブロック図である。この実施の形態では、基準信号RFを基準発振器18から分周器52で分周したものとする点で、第1の実施の形態と異なっている。
【0033】
即ち、基準発振器18からの60MHzの信号は、分周器52で333kHzの信号に分周されて、直接、位相差演算回路48に入力されて、測定信号S、または参照信号Rとの位相差Δθ1、Δθ2(またはΔθ3、Δθ4)が求められる。この構成においても、上記実施の形態と同様に位相差を求めることができると共に、ミキサ32、バンドパスフィルタ36、アンプ40のアナログ回路1系統だけでよいため、回路構成が単純になり、製造コストが安価になるという効果を有している。
【0034】
図5は、本発明のさらに他の実施の形態を表すブロック図である。この実施の形態では、参照ラインRLに設けられた減衰器を異なる減衰率を持つ複数の減衰器28−1、28−2、28−3から構成し、適宜、スイッチ54、55でいずれかの減衰器を選択するように構成すると共に、アンプ40からの出力をピーク検出回路56、A/D変換器58を介して中央演算回路50に送出している。
【0035】
本実施の形態では、測定ラインSLが選択されている状態で、ピーク検出回路56でアンプ40から出力される波形のピークを検出し、このピークをA/D変換器58でA/D変換する。中央演算回路50では、ピークの大きさによって、減衰器28−1、28−2、28−3のいずれを選択するかを決定し、スイッチ54、55を切り換える。こうして、参照ラインRLが選択されたときに、選択された減衰器の減衰を受け、測定信号Sと同じ程度の電力を持つ参照信号Rが得られるようにしている。
【0036】
電力によってもオフセット誤差が変化するので、このように、測定信号Sの電力に応じて減衰量が変化するように構成することにより、測定信号Sと参照信号Rがそれぞれ包含する誤差θe1とθe2とをできる限り同じ値にすることができ、さらに精度良いものとすることができる。
【0037】
【発明の効果】
以上説明したように、請求項1記載の発明によれば、受信器からのマイクロ波信号を第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換し、さらに、この低周波信号を矩形波に変換した後、位相差演算回路で、被測定流体を伝搬しない信号と、前記波形整形回路からの矩形波信号とが入力され、両信号の位相差を、前記基準発振器からの基準周波数の信号を用いて計数することにより、従来の濃度計に比べて、アナログ回路を少なくすることができるため、アナログ回路の特性に起因する誤差が少なくなり、精度を向上させることができる。
【0038】
また、第1マイクロ波発振器と第2マイクロ波発振器からのマイクロ波信号は基準発振器からの信号にフェーズロックされ、且つ位相差演算回路は、基準発振器からの信号を用いて位相差を検出することから、例えば、基準発振器の信号がゆらいでも、全てが同一の影響を受けるため、結果として位相差演算回路から出力される位相差は、基準発振器のゆらぎの影響を受けずにすむ、という従来にない優れた効果が得られる。
【0039】
また、請求項2記載の発明によれば、被測定流体を伝搬しない信号として、前記第1マイクロ波発振器からのマイクロ波信号を、前記第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換し、さらに、この低周波信号を矩形波信号に変換したものとし、位相差演算回路で、前記被測定流体を伝搬してきた波形整形回路からの矩形波信号との間での位相差を計数することにより、従来の濃度計に比べて、アナログ回路を少なくすることができるため、アナログ回路の特性に起因する誤差が少なくなり、精度を向上させることができる。
【0040】
また、請求項3記載の発明によれば、被測定流体を伝搬しない信号として、基準発振器からの基準周波数の信号を|f1−f2|に分周した信号とすることにより、被測定流体を伝搬した信号を処理するアナログ回路1系統だけでよいため、回路構成が単純になり、製造コストが安価になるという効果を有している。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態を示すブロック図である。
【図2】図1の位相差演算回路で処理される信号波形を表しており、(a)は基準信号RF、測定信号S及びクロック信号CLとの関係を表し、図2(b)は基準信号RF、参照信号R及びクロック信号CLとの関係を表し、図2(c)(d)はスイッチからの漏れがあるときの図2(a)(b)対応図である。
【図3】(a)は、スイッチからの漏れがあるときの測定信号Sと位相誤差θe3とをベクトル表示した説明図、(b)は、位相誤差θe3の位相依存性を示す説明図である。
【図4】本発明の第2の実施の形態を示すブロック図である。
【図5】本発明の第3の実施の形態を示すブロック図である。
【図6】θsの実測値の時間変化を表すグラフである。
【図7】θrの実測値の時間変化を表すグラフである。
【図8】θs−θrの実測値の時間変化を表すグラフである。
【図9】従来の濃度計のブロック図を示す。
【符号の説明】
12 送信アンテナ(送信器)
14 受信アンテナ(受信器)
18 基準発振器
22 シンセサイザ(第1マイクロ波発振器)
24 シンセサイ(第2マイクロ波発振器)
32 ミキサ(周波数変換器)
34 ミキサ(第2周波数変換器)
44 コンパレータ(波形整形回路)
46 コンパレータ(第2波形整形回路)
48 位相差演算回路[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a densitometer that measures the concentration of a fluid to be measured using a microwave, and more particularly to a densitometer that can improve accuracy.
[0002]
[Prior art]
Conventionally, this type of densitometer is shown in FIG. This densitometer measures the phase lag of the received wave that has propagated through the substance to be measured by injecting the microwave into the fluid to be measured that flows through the pipe. Furthermore, the reference phase lag (for example, when water flows through the same pipe) In this case, the concentration of the substance to be measured is measured from the difference from the phase delay.
[0003]
As a specific device configuration, a transmission antenna 12 and a reception antenna 14 are attached to the detection pipe 10 so as to face each other. For example, a microwave signal of 1.6 GHz output from the high-frequency oscillator 60 is sent to the transmission antenna 12, and the microwave emitted from the transmission antenna 12 propagates through the fluid to be measured in the pipe 10 to receive the reception antenna 14. Received at. The reception signal from the reception antenna 14 passes through the amplifier 61 and then is sent to the mixer 62. At the same time, the microwave signal from the high-frequency oscillator 60 is sent to another mixer 64, and in both mixers 62 and 64, 1.. It is mixed with the output from a high-frequency oscillator 66 that oscillates a microwave signal having a frequency of 6 GHz-60 MHz, and the beat signal is extracted through band-pass filters 68, 70, respectively, and amplifiers 72, 74, limiters 76, 78, etc. To the quadrature synchronous detector 80. A quadrature synchronous detector 80 outputs a sine component sinΔθ and a cosine component cosΔθ of two input phase differences Δθ, which are A / D converted by A / D converters 82 and 84, respectively, and a phase difference calculation circuit 86
[0004]
[Expression 1]
Δθ = tan -1 (sinθ / cosθ)
Is calculated to obtain the phase difference, that is, the phase delay of the microwave propagated through the fluid to be measured.
[0005]
[Problems to be solved by the invention]
In such a conventional densitometer, in particular, since a voltage offset error occurs in the quadrature synchronous detector 80, in order to suppress the error, it is necessary to keep the voltage constant by using the limiters 76 and 78. There is a problem that it becomes complicated and costs increase.
[0006]
In addition, since the two high-frequency oscillators 60 and 66 are used, there is a problem that an error occurs when these outputs fluctuate individually.
In addition to the above errors, since many analog circuits are used, the error factors of each analog circuit overlap, and the resulting phase difference includes voltage offset error, phase linearity error, phase offset error, frequency offset Since errors etc. are included, in order to calibrate these errors, the power, phase, and frequency are used as parameters to obtain the error when these are changed, and a correction curve is created and measured. Sometimes the correction result is corrected using this correction curve.
[0007]
However, it takes time and labor to create such a correction curve, and the characteristics are different for each analog circuit, so that there is a problem that a correction curve must be created for each densitometer. In addition, even if a correction curve is prepared, there is a problem that the error characteristic changes due to aging or temperature change of the analog circuit. The temperature change can be dealt with by storing the apparatus in a container heated to a constant temperature with a heater, but there is also a problem that it costs because it requires a heater or the like.
[0008]
The present invention has been made in view of such problems, and the inventions according to claims 1 to 3 provide a densitometer capable of measuring the concentration with high accuracy without preparing correction data in advance. The purpose is to do.
[0009]
[Means for Solving the Problems]
In order to achieve the above object, the invention according to claim 1 of the present invention is characterized in that the concentration of the fluid to be measured is detected by detecting the phase difference between the microwave that has propagated through the fluid to be measured and the signal that does not propagate through the fluid to be measured. A reference oscillator that oscillates a signal at a reference frequency, a first microwave oscillator that outputs a microwave signal having a frequency f 1 phase-locked from the reference frequency, and a phase-lock from the reference frequency A second microwave oscillator that outputs a microwave signal having a frequency f 2 different from the microwave signal; a transmitter that is connected to the first microwave oscillator and radiates the microwave to the fluid to be measured; and a fluid to be measured A receiver that receives the microwave propagated through the microwave, and a microwave signal from the receiver is mixed with a second microwave signal from the second microwave oscillator; a frequency converter for converting to a low frequency signal of f 1 −f 2 |, a waveform shaping circuit for converting the low frequency signal from the frequency converter into a rectangular wave, a signal that does not propagate through the fluid to be measured, and the waveform shaping A rectangular wave signal from the circuit, and a phase difference calculation circuit that counts a phase difference between the two signals using a reference frequency signal from the reference oscillator.
[0010]
The invention according to claim 2 is the densitometer according to claim 1, wherein the microwave signal from the first microwave oscillator is mixed with the second microwave signal from the second microwave oscillator, and | a second frequency converter for converting into a low frequency signal of f 1 −f 2 |, and a second waveform shaping circuit for converting the low frequency signal from the second frequency converter into a rectangular wave signal, A rectangular wave signal from a two-waveform shaping circuit is input to the phase difference calculation circuit as a signal that does not propagate through the fluid to be measured.
[0011]
According to a third aspect of the present invention, there is provided the densitometer according to the first aspect , wherein a signal obtained by dividing a reference frequency signal from a reference oscillator into | f 1 −f 2 | is a signal that does not propagate through the fluid to be measured. As input to the phase difference calculation circuit.
[0012]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing a first embodiment of the present invention. In FIG. 1, a transmitting antenna 12 and a receiving antenna 14 are arranged so as to face each other in a detection pipe 10 through which a fluid to be measured flows. It is attached. The portion of the pipe 10 to which the transmission antenna 12 and the reception antenna 14 are attached is an insulator.
[0013]
A PLL synthesizer 22 and a PLL synthesizer 24 are connected from a reference oscillator 18 such as a crystal oscillator via a frequency divider 20, and the synthesizer 22 is connected to the transmission antenna 12 via a switch 26 described later. A signal of 60 MHz, for example, output from the reference oscillator 18 is frequency-divided into 1/6 by the frequency divider 20 and then converted to a 1.6 GHz microwave signal phase-locked by the synthesizer 22. And a microwave is emitted from the transmission antenna 12. The synthesizer 24 outputs a 1.6 GHz-333 kHz phase-locked signal that is very close to the synthesizer 22.
[0014]
A received signal from the receiving antenna 14 is input to a mixer (frequency converter) 32 through a switch 30, where it is mixed with a signal from the synthesizer 24, passes through a bandpass filter 36, and becomes a low frequency signal of 333 kHz. Converted. At the same time, the signal from the synthesizer 22 is directly input to the mixer (frequency converter) 34, mixed with the signal from the synthesizer 24, and converted to a low frequency signal of 333 kHz through the band pass filter 38. The signals converted into these low frequency signals are connected to amplifiers 40 and 42, comparators 44 and 46, and a phase difference calculation circuit 48, respectively. The comparators 44 and 46 are comparators for detecting a zero crossing of the 333 kHz sine wave signal sent from the amplifiers 40 and 42, thereby shaping the sine wave into a rectangular wave and sending it to the phase difference calculation circuit 48. is there.
[0015]
The phase difference calculation circuit 48 is a digital signal processor composed of an ASIC, and uses the signal from the reference oscillator 18 as a clock signal and counts the time difference between the digital waveforms input from the comparator 44 and the comparator 46, thereby synthesizer. The phase difference between the signal that has passed through the fluid to be measured from 22 and the signal that has been sent directly from the synthesizer 22 without passing through the fluid to be measured is detected. The detected phase difference is further sent to the central processing circuit 50.
[0016]
In this embodiment, between the transmission and the synthesizer 22 antenna 12, between the receiving antenna 14 and the mixer 32, respectively, switch 26, and switch 30 is interposed, one contact S 1 of the switch 26 to the transmitting antenna 12 One contact S 2 of the switch 30 is connected to the receiving antenna 14, and the line from the synthesizer 22 including the transmitting antenna 12 and the receiving antenna 14 between these S 1 -S 2 to the phase difference calculation circuit 48 is The measurement line SL. Further, the other contact R 1 of the switch 26, connecting the other contact R 2 of the switch 30, between the synthesizer 22 to the phase difference calculation circuit 48 does not pass through the fluid to be measured, to form a reference line RL The reference line RL is provided with an attenuator 28 that gives the same degree of attenuation as when the fluid to be measured has propagated. Switching of the switch 26 and the switch 30 is performed at high speed based on the control signal from the central processing circuit 50, and the microwave signal from the synthesizer 22 is selected and transmitted to either the measurement line SL or the reference line RL. The
[0017]
In the densitometer configured as described above, first, in a state in which the microwave signal from the synthesizer 22 is switched by the switches 26 and 30 so as to propagate to the measurement line SL, a 1.6 GHz microwave signal is transmitted. The microwave transmitted to the antenna 12 and emitted from the transmission antenna propagates through the fluid to be measured in the pipe 10 and is received by the reception antenna 14. The mixer 32, the band pass filter 36, the amplifier 40, and the comparator 44 are received. Then, the signal is input to the phase difference calculation circuit 48 as the measurement signal S. At the same time, the reference signal RF that has passed through the mixer 34, the band-pass filter 38, the amplifier 42, and the comparator 46 without passing through the fluid to be measured from the synthesizer 22 is also input to the phase difference calculation circuit 48, and measurement from the reference signal RF is performed. The phase lag (phase difference) of the signal S is obtained. That is, the phase difference is
[0018]
[Expression 2]
Δθ 1 = θs−θ 0 + θe 1
It is represented by Here, Δθ 1 is the phase difference obtained by the phase difference calculation circuit 48, θs is the phase of the measurement signal, θ 0 is the phase of the reference signal, θe 1 is due to drift of the measurement system, offset error, and temperature change. It is an error component due to errors, variations in analog circuits, and the like.
[0019]
Next, when the switches 26 and 30 are switched so that the microwave signal from the synthesizer 22 propagates through the reference line RL, the microwave signal passes through the attenuator 28, and the mixer 32, the bandpass filter 36, and the amplifier 40 and the comparator 44, and then input to the phase difference calculation circuit 48 as a reference signal R. At the same time, the reference signal RF is also input to the phase difference calculation circuit 48. At this time, the phase lag (phase difference) obtained by the phase difference calculation circuit 48 is:
[0020]
[Equation 3]
Δθ 2 = θr−θ 0 + θe 2
It is represented by Here, Δθ 2 is the phase lag obtained by the phase difference calculation circuit 48, θr is the phase of the reference signal R that has passed through the reference line RL, θ 0 is the phase of the reference signal, and θe 2 is the phase of the measurement system. It is an error component due to drift, offset error, temperature change error, and analog circuit variation.
[0021]
Here, the measurement signal S passing through the measurement line SL and the reference signal R passing through the reference line RL are alternately obtained by time division, but the measurement line SL and the reference line RL share the same analog circuit. It is considered that measurement system drift, offset error, temperature error, and analog circuit variation error are equivalent. That is,
[0022]
[Expression 4]
θe 1 = θe 2
Is considered to hold. Therefore, when the difference between Δθ 1 and Δθ 2 is obtained,
[0023]
[Equation 5]
Figure 0003804129
It becomes. Since θr is known in advance from the electrical length of the reference line RL, a value corresponding to θs can be obtained. Further, the concentration of the fluid to be measured is measured from the difference between them by flowing water through the same pipe 10 and obtaining a value corresponding to θs at that time.
[0024]
FIG. 2 shows the relationship of signals input to the phase difference calculation circuit 48, FIG. 2 (a) shows the relationship between the reference signal RF, the measurement signal S, and the clock signal CL, and FIG. 2 (b) shows the reference. The relationship among the signal RF, the reference signal R, and the clock signal CL is represented. The phase difference between the reference signal RF and the measurement signal S counts the number of clocks (ta) that exist between the rising edges of the two signals, and the number of clocks that exist between the rising edges of one signal (referred to as ta) [Delta] [theta] 1 is obtained by counting (tc). That is,
[0025]
[Formula 6]
Δθ 1 = 2π · ta / tc
It is. Similarly, the phase difference between the reference signal RF and the reference signal R counts the number of clocks (td) existing between the rising edges of two signals, and between the rising edges of one signal. By counting the number of existing clocks (denoted by tc), Δθ 2 is obtained. That is,
[0026]
[Expression 7]
Δθ 2 = 2π · td / tc
It is.
[0027]
FIGS. 2C and 2D are diagrams corresponding to FIGS. 2A and 2B when the switches 26 and 30 are leaking. Although the measurement lines SL and the reference lines RL are switched by the switches 26 and 30, if the isolation of the switches 26 and 30 is insufficient, leakage from the other line interferes and vector synthesis is performed. When an explanatory diagram vector representation signal in FIG. 3 (a), summed connected to now the measurement line SL, When the measurement signals S 1 is obtained, the leakage signal E from the reference line RL is with small As a result, the phase is increased by θe 3 and an error occurs. At S 2 where the phase is shifted from the measurement signal S 1 by 180 degrees, the phase is reduced by θe 3 due to the leakage signal E from the reference line RL. The phase dependence of the phase error θe 3 is shown in FIG. 3B, and it can be seen that the phase error θe 3 has the same absolute value and the sign is inverted where the phase is different by 180 degrees. Therefore, by utilizing this fact, in FIG. 2C, the number of clocks (ta) existing between the rising edge of the reference signal RF and the rising edge of the measurement signal S, the rising edge of the reference signal RF, and the rising edge of the measurement signal S are obtained. The number of clocks (tb) existing during the fall is obtained and averaged.
[0028]
[Equation 8]
Δθ 3 = 2π · (ta + tb) / 2tc
Similarly, the process is performed for the reference signal R.
[0029]
[Equation 9]
Δθ 4 = 2π · (td + te) / 2tc
By taking the difference between Δθ 3 and Δθ 4 , θs−θr can be obtained. In other words, by obtaining the phase difference θs−θr based on the time point of 90 degrees of each of the measurement signal S and the reference signal R, it is not affected by errors due to leakage of the switches 26 and 30. Can be.
[0030]
As described above, the frequency output from the synthesizer 22 and the synthesizer 24 is very close (333 kHz), and the difference is set so as to be a frequency that can be processed by the digital signal processor. Can be used as a digital signal processor, so that the error due to the characteristics of the analog circuit is reduced and the accuracy can be improved as compared with the case where many conventional analog circuits are used.
[0031]
Further, since the microwave signals from the synthesizers 22 and 24 are phase-locked to the signal from the reference oscillator 18 and the phase difference calculation circuit 48 detects the phase difference using the signal from the reference oscillator 18, for example, Even if the signal of the reference oscillator 18 fluctuates, all are affected by the same effect, and as a result, the phase difference output from the phase difference calculation circuit 48 is not affected by the fluctuation of the reference oscillator 18.
[0032]
FIGS. 6 and 7 show the time changes of θs and θr, respectively, and fluctuations in the range of 0.18 degrees to 0.19 degrees are observed, as shown in FIG. By taking the difference, it was possible to suppress the fluctuation within the range of 0.06 degrees.
FIG. 4 is a block diagram showing another embodiment of the present invention. This embodiment is different from the first embodiment in that the reference signal RF is divided from the reference oscillator 18 by the frequency divider 52.
[0033]
That is, the 60 MHz signal from the reference oscillator 18 is frequency-divided into a 333 kHz signal by the frequency divider 52 and directly input to the phase difference calculation circuit 48, and the phase difference from the measurement signal S or the reference signal R. Δθ 1 and Δθ 2 (or Δθ 3 and Δθ 4 ) are obtained. Also in this configuration, the phase difference can be obtained in the same manner as in the above embodiment, and only one analog circuit system of the mixer 32, the band pass filter 36, and the amplifier 40 is required, so that the circuit configuration is simplified and the manufacturing cost is reduced. Has the effect of becoming inexpensive.
[0034]
FIG. 5 is a block diagram showing still another embodiment of the present invention. In this embodiment, the attenuator provided in the reference line RL is composed of a plurality of attenuators 28-1, 28-2, 28-3 having different attenuation rates, and any one of the switches 54, 55 is used as appropriate. The attenuator is selected, and the output from the amplifier 40 is sent to the central processing circuit 50 via the peak detection circuit 56 and the A / D converter 58.
[0035]
In the present embodiment, the peak of the waveform output from the amplifier 40 is detected by the peak detection circuit 56 with the measurement line SL selected, and this peak is A / D converted by the A / D converter 58. . The central processing circuit 50 determines which one of the attenuators 28-1, 28-2, 28-3 is selected according to the size of the peak, and switches the switches 54, 55. Thus, when the reference line RL is selected, the reference signal R having the same power as that of the measurement signal S is obtained by being attenuated by the selected attenuator.
[0036]
Since the offset error also changes depending on the power, the errors θe 1 and θe included in the measurement signal S and the reference signal R are configured by changing the attenuation amount according to the power of the measurement signal S in this way. 2 can be set to the same value as much as possible, and the accuracy can be further improved.
[0037]
【The invention's effect】
As described above, according to the first aspect of the present invention, the microwave signal from the receiver is mixed with the second microwave signal from the second microwave oscillator, and the low frequency of | f 1 −f 2 | After converting the low-frequency signal into a rectangular wave, a signal that does not propagate through the fluid to be measured and a rectangular wave signal from the waveform shaping circuit are input to the phase difference calculation circuit. By using the reference frequency signal from the reference oscillator to count the phase difference of the analog circuit, the analog circuit can be reduced compared to the conventional densitometer, so there is less error due to the characteristics of the analog circuit. Thus, the accuracy can be improved.
[0038]
The microwave signals from the first microwave oscillator and the second microwave oscillator are phase-locked to the signal from the reference oscillator, and the phase difference calculation circuit detects the phase difference using the signal from the reference oscillator. Therefore, for example, even if the signal of the reference oscillator fluctuates, all are affected by the same effect, and as a result, the phase difference output from the phase difference calculation circuit is not affected by the fluctuation of the reference oscillator. There is no excellent effect.
[0039]
According to a second aspect of the present invention, a microwave signal from the first microwave oscillator is mixed with a second microwave signal from the second microwave oscillator as a signal that does not propagate through the fluid to be measured. , | F 1 −f 2 | is converted into a low-frequency signal, and this low-frequency signal is converted into a rectangular wave signal. The phase difference calculation circuit uses the waveform shaping circuit that has propagated the fluid to be measured. By counting the phase difference from the square wave signal, the analog circuit can be reduced compared to conventional densitometers, so errors due to the characteristics of the analog circuit are reduced and accuracy is improved. Can be made.
[0040]
According to the third aspect of the present invention, the signal to be measured is obtained by dividing the signal of the reference frequency from the reference oscillator into | f 1 −f 2 | as the signal that does not propagate through the fluid to be measured. Since only one analog circuit for processing the signal propagated through the circuit is required, the circuit configuration is simplified and the manufacturing cost is reduced.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a first embodiment of the present invention.
2 shows signal waveforms processed by the phase difference calculation circuit of FIG. 1, wherein (a) shows the relationship between the reference signal RF, the measurement signal S, and the clock signal CL, and FIG. 2 (b) shows the reference. The relationship between the signal RF, the reference signal R, and the clock signal CL is shown, and FIGS. 2C and 2D are diagrams corresponding to FIGS. 2A and 2B when there is leakage from the switch.
FIG. 3A is an explanatory diagram in which a measurement signal S and a phase error θe 3 when there is leakage from the switch are displayed as vectors, and FIG. 3B is an explanatory diagram showing the phase dependence of the phase error θe 3 . It is.
FIG. 4 is a block diagram showing a second embodiment of the present invention.
FIG. 5 is a block diagram showing a third embodiment of the present invention.
FIG. 6 is a graph showing a change over time of an actually measured value of θs.
FIG. 7 is a graph showing a change over time of an actually measured value of θr.
FIG. 8 is a graph showing a time change of an actual measurement value of θs−θr.
FIG. 9 shows a block diagram of a conventional densitometer.
[Explanation of symbols]
12 Transmitting antenna (transmitter)
14 Receiving antenna (receiver)
18 Reference oscillator 22 Synthesizer (first microwave oscillator)
24 synthesizer The (second microwave oscillator)
32 Mixer (frequency converter)
34 Mixer (second frequency converter)
44 Comparator (waveform shaping circuit)
46 Comparator (second waveform shaping circuit)
48 Phase difference calculation circuit

Claims (3)

被測定流体を伝搬したマイクロ波と被測定流体を伝搬しない信号との位相差を検出することにより、被測定流体の濃度を測定する濃度計において、
基準周波数の信号を発振する基準発振器と、
該基準周波数からフェーズロックされた周波数f1のマイクロ波信号を出力する第1マイクロ波発振器と、
該基準周波数からフェーズロックされた前記マイクロ波信号と異なる周波数f2のマイクロ波信号を出力する第2マイクロ波発振器と、
前記第1マイクロ波発振器に接続され、マイクロ波を被測定流体に放射する送信器と、
被測定流体を伝搬したマイクロ波を受波する受信器と、
前記受信器からのマイクロ波信号を前記第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換する周波数変換器と、
前記周波数変換器からの低周波信号を矩形波に変換する波形整形回路と、
被測定流体を伝搬しない信号と、前記波形整形回路からの矩形波信号とが入力され、両信号の位相差を、前記基準発振器からの基準周波数の信号を用いて計数する位相差演算回路と、
を備えることを特徴とする濃度計。
In a densitometer that measures the concentration of the fluid under measurement by detecting the phase difference between the microwave that has propagated through the fluid under measurement and the signal that does not propagate through the fluid under measurement,
A reference oscillator that oscillates a signal of a reference frequency;
A first microwave oscillator for outputting a microwave signal having a frequency f 1 phase-locked from the reference frequency;
A second microwave oscillator for outputting a microwave signal having a frequency f 2 different from the microwave signal phase-locked from the reference frequency;
A transmitter connected to the first microwave oscillator for emitting microwaves to a fluid to be measured;
A receiver for receiving the microwave propagated through the fluid to be measured;
A frequency converter for mixing the microwave signal from the receiver with the second microwave signal from the second microwave oscillator and converting it to a low frequency signal of | f 1 −f 2 |;
A waveform shaping circuit for converting a low-frequency signal from the frequency converter into a rectangular wave;
A signal that does not propagate through the fluid to be measured and a rectangular wave signal from the waveform shaping circuit are input, and a phase difference calculation circuit that counts the phase difference between the two signals using a reference frequency signal from the reference oscillator;
A densitometer characterized by comprising:
請求項1記載の濃度計において、
前記第1マイクロ波発振器からのマイクロ波信号を、前記第2マイクロ波発振器からの第2マイクロ波信号と混合し、|f1−f2|の低周波信号に変換する第2周波数変換器と、
前記第2周波数変換器からの低周波信号を矩形波信号に変換する第2波形整形回路と、
を備え、該第2波形整形回路からの矩形波信号を前記被測定流体を伝搬しない信号として、前記位相差演算回路に入力することを特徴とする濃度計。
The densitometer according to claim 1, wherein
A second frequency converter for mixing the microwave signal from the first microwave oscillator with the second microwave signal from the second microwave oscillator and converting it to a low frequency signal of | f 1 −f 2 |; ,
A second waveform shaping circuit for converting a low frequency signal from the second frequency converter into a rectangular wave signal;
And a rectangular wave signal from the second waveform shaping circuit is input to the phase difference calculation circuit as a signal that does not propagate through the fluid to be measured.
請求項1記載の濃度計において、
基準発振器からの基準周波数の信号を|f1−f2|に分周した信号を、前記被測定流体を伝搬しない信号として、前記位相差演算回路に入力することを特徴とする濃度計。
The densitometer according to claim 1, wherein
A densitometer, wherein a signal obtained by dividing a reference frequency signal from a reference oscillator into | f 1 −f 2 | is input to the phase difference calculation circuit as a signal that does not propagate through the fluid to be measured.
JP30019896A 1996-11-12 1996-11-12 Densitometer Expired - Fee Related JP3804129B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP30019896A JP3804129B2 (en) 1996-11-12 1996-11-12 Densitometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP30019896A JP3804129B2 (en) 1996-11-12 1996-11-12 Densitometer

Publications (2)

Publication Number Publication Date
JPH10142168A JPH10142168A (en) 1998-05-29
JP3804129B2 true JP3804129B2 (en) 2006-08-02

Family

ID=17881921

Family Applications (1)

Application Number Title Priority Date Filing Date
JP30019896A Expired - Fee Related JP3804129B2 (en) 1996-11-12 1996-11-12 Densitometer

Country Status (1)

Country Link
JP (1) JP3804129B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102474924B (en) 2009-09-29 2013-08-14 松下电器产业株式会社 High-frequency heating device and high-frequency heating method

Also Published As

Publication number Publication date
JPH10142168A (en) 1998-05-29

Similar Documents

Publication Publication Date Title
CN107271938B (en) A kind of frequency mixer vectorial property measurement method
US5210539A (en) Linear frequency sweep synthesizer
US9893716B2 (en) Resonant circuit calibration
JPH0382248A (en) 90× phase shifter and data receiver
JPH09119977A (en) Filling-level measuring method of liquid
US4447782A (en) Apparatus for automatic measurement of equivalent circuit parameters of piezoelectric resonators
EP0913031B1 (en) Method and system for tuning resonance modules
JP3804129B2 (en) Densitometer
JP5109492B2 (en) Radar equipment
US4812849A (en) Nonquadrature correction circuit
JP3804128B2 (en) Densitometer
US6449307B1 (en) Apparatus for measuring electrical characteristics of circuit
US2480128A (en) Frequency measuring system
US6246727B1 (en) Method and system for tuning resonance modules
JP3629252B2 (en) Microwave concentration measuring device
JP2623931B2 (en) Phase-synchronized transceiver
JP2013205101A (en) Radio terminal distance measuring system, and distance measuring system
US2820898A (en) Distance measuring equipment utilizing frequency modulation
TW202142891A (en) Six-port self-injection-locked radar
JPH02280083A (en) Measuring instrument for water vapor content and rainfall quantity
CN110476057B (en) Flow measurement system
JP3474509B2 (en) Microwave concentration measurement device
SU681388A1 (en) Digital phase meter
JP2004177370A (en) Microwave density measuring apparatus
CN214412705U (en) 5MHz, 10MHz outer reference self-adaptation frequency reference device

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060117

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060317

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060411

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060501

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20090519

Year of fee payment: 3

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

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

Free format text: PAYMENT UNTIL: 20090519

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

Free format text: PAYMENT UNTIL: 20090519

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20100519

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20110519

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20110519

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20120519

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20120519

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20130519

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20130519

Year of fee payment: 7

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

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

Free format text: PAYMENT UNTIL: 20130519

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

Free format text: PAYMENT UNTIL: 20130519

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20140519

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees