JP3492347B2 - 水中音響伝播特性模擬のための装置、方法及びプログラム - Google Patents

水中音響伝播特性模擬のための装置、方法及びプログラム

Info

Publication number
JP3492347B2
JP3492347B2 JP2001363149A JP2001363149A JP3492347B2 JP 3492347 B2 JP3492347 B2 JP 3492347B2 JP 2001363149 A JP2001363149 A JP 2001363149A JP 2001363149 A JP2001363149 A JP 2001363149A JP 3492347 B2 JP3492347 B2 JP 3492347B2
Authority
JP
Japan
Prior art keywords
propagation
propagation loss
calculating
loss
parameter input
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
JP2001363149A
Other languages
English (en)
Other versions
JP2003161773A (ja
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.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries Ltd
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 Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2001363149A priority Critical patent/JP3492347B2/ja
Publication of JP2003161773A publication Critical patent/JP2003161773A/ja
Application granted granted Critical
Publication of JP3492347B2 publication Critical patent/JP3492347B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は水中音響伝播特性模
擬のための装置、方法及びプログラムに係り、特に少な
い計算量で模擬精度を向上することの可能な水中音響伝
播特性模擬のための装置、方法及びプログラムに関す
る。
【0002】
【従来の技術】水中を走行する水中走行体は、自らが発
信する音波の障害物による反射波を受信して障害物を認
識するアクティブソナー、又は障害物自体が発信する音
響信号を受信して障害物を認識するパッシブソナーを使
用して障害物の有無並びに位置及び大きさを認識してい
る。
【0003】従ってソナーの取扱いに習熟することが必
要であるが、水中走行体を実際に水中走行させて取扱い
習熟のための訓練、あるいはソナーの性能評価を行うこ
とは費用がかかるだけでなく、同一条件化で繰り返し訓
練又は評価を行うことが困難であるため、従来から水中
音響伝播特性模擬装置を使用して訓練又は評価が行われ
ている。
【0004】水中音響伝播特性は水中の音源から受信点
に至る音波の伝播経路及び音波の伝播時間並びに水中の
伝播損失によって決定される。そして従来の水中音響伝
播特性模擬装置においては、いわゆる音線法によって伝
播経路、伝播時間及び伝播損失によって求めていた。
【0005】
【発明が解決しようとする課題】しかしながら、音線法
によれば伝播経路及び伝播時間は比較的簡易な計算によ
り正確に求めることができるものの、音線が極めて疎あ
るいは密である領域において伝播損失を正確に求めるこ
とは困難である。
【0006】さらに、実際の伝播損失にゆらぎが存在す
るが、音線法ではこのゆらぎを模擬することは不可能で
あった。
【0007】音線が極めて疎あるいは密である領域にお
ける伝播損失を正確に求める方法としてはノーマルモー
ド法が周知であるが、ノーマルモード法により伝播経路
及び伝播時間を計算すると、計算量は膨大となる。
【0008】表1は音線法とノーマルモード法の比較結
果であって、音線法は伝播経路並びに伝播時間の計算に
は適しているものの、伝播損失を厳密に計算することは
困難である。これに対し、ノーマルモード法は伝播損失
を厳密に計算できるものの、伝播経路及び伝播時間の計
算には時間がかかることが判る。
【0009】
【表1】
【0010】本発明は、上記課題に鑑みなされたもので
あって、伝播経路及び伝播時間を音線法で算出し、伝播
損失をノーマルモード法で算出することにより少ない計
算量で模擬精度を向上することの可能な水中音響伝播特
性模擬のための装置、方法及びプログラムを提供するこ
とを目的とする。
【0011】
【課題を解決するための手段】第一の発明に係る水中音
響伝播特性模擬のための装置は、水質及び水中音源に関
するパラメータを入力するパラメータ入力手段と、パラ
メータ入力手段から入力されたパラメータに基づいて音
線法により伝播経路及び伝播時間を計算する伝播経路・
伝播時間計算手段と、パラメータ入力手段から入力され
たパラメータに基づいてノーマルモード法により前記伝
播経路・伝播時間計算手段で計算された伝播経路に対応
する伝播損失を計算する伝播損失計算手段と、伝播経路
・伝播時間計算手段で計算された伝播経路及び伝播時間
と伝播損失計算手段で計算された伝播損失を記憶する記
憶手段と、受信位置に関するパラメータを入力する受信
位置パラメータ入力手段と、該水中音源から発生された
音の、受信位置パラメータ入力手段で入力された受信位
置における受信信号を計算する受信信号計算手段を具備
する。
【0012】本発明にあっては、音源から受信機に伝播
する音波の伝播経路及び伝播時間は音線法により計算さ
れ、伝播損失はノーマルモード法により計算される。
【0013】第二の発明に係る水中音響伝播特性模擬の
ための装置は、伝播損失計算手段が、計算された伝播損
失の中心値と該中心値周りの変動成分に分解する伝播損
失分解手段と、伝播損失分解手段で分解された変動成分
の標準偏差を計算する標準偏差計算手段を具備する。
【0014】本発明にあっては、伝播損失が中心値とこ
の中心値周りの変動成分に分解され、それぞれを屈折に
起因する伝播損失及び干渉に起因する伝播損失とみな
す。
【0015】第三の発明に係る水中音響伝播特性模擬の
ための装置は、伝播損失計算手段が、該水中音源が発生
する音の周波数の1/10から1/100の周波数を使
用して伝播損失を計算する。
【0016】本発明にあっては、伝播損失は実際に水中
音源が発生する音の周波数の1/10から1/100の
周波数を使用して計算される。
【0017】
【発明の実施の形態】図1は本発明に係る水中音響伝播
特性模擬装置の構成図であって、ソナーの性能評価を行
う場合を示す。
【0018】即ち、水中音響伝播特性模擬装置10は一
種のディジタル計算機であって、バス100を中心とし
てCPU101、メモリ102、制御卓用インターフェ
イス103及びソナーインターフェイス104から構成
される。そして制御卓用インターフェイス103には制
御卓11が、ソナーインターフェイス104にはソナー
12が接続される。
【0019】制御卓11は試験条件を設定するためのも
のであり、音源の位置、深度、発生音圧及びソナー受信
機の座標等を設定する。メモリ102内に記憶されてい
る水中音響伝播特性模擬ルーチンはCPU101で実行
され、ソナーインターフェイス104を介してソナー1
2に対してソナー12の受信機が受信する音響信号を出
力する。
【0020】図2は水中音響伝播特性模擬ルーチンのフ
ローチャートであって、ステップ21で制御卓11から
模擬領域の水質に関するパラメータ(温度、塩分濃度
等)、並びに音源に関するパラメータ(音源深度、音源
音圧等)を設定する。そして、ステップ22で音線法に
より伝播経路及び伝播時間を算出する。
【0021】図3は模擬方法の説明図であって、右向き
に距離軸(r軸)、下向きに深度軸(d軸)をとる。な
お、原点は音源S直上の水面とする。
【0022】即ち、模擬水域を所定深度及び所定距離ご
とに格子状に分割し、音源Sと各格子点(ri,dj)を
結ぶ伝達経路について伝播時間τ(ri,dj)を算出す
る。
【0023】次にステップ23でノーマルモード法によ
り、音源Sから格子点(ri,dj)に至る伝達経路の伝
播損失の厳密解L(ri,dj)を算出する。
【0024】図4は音波周波数を変更した場合の伝播損
失の計算結果であって、音波周波数をf、2f及び4f
とした場合を示す。この図から判るように、音波の周波
数が高くなると伝播損失の深さ方向の変動周波数は高く
なるものの、伝播損失の深さ方向平均値の形状及び平均
値周りの変動幅は変化しない。これに対し、ノーマルモ
ード法による伝播損失の計算時間は、音波周波数のほぼ
二乗に比例して長くなる。
【0025】従って伝播損失の計算にあたっては、実際
に使用する音波の周波数の1/10〜1/100の周波
数を使用することが、伝播損失の計算精度を損なわずに
計算時間を短縮するうえで有利である。例えば伝播損失
計算に使用する音波周波数を実際の周波数の1/10と
すれば、計算時間を1/100に短縮することが可能と
なる。
【0026】ノーマルモード法によれば、音速勾配によ
る屈折に起因する伝播損失Lr だけでなく、直接波と水
面もしくは水底で反射した反射波の干渉に起因する伝播
損失Li も厳密に計算される。従って、格子点を特定す
れば屈折に起因する伝播損失Lr 及び干渉に起因する伝
播損失Li を重畳した伝播損失は、一義的に定まる。し
かし実際にはソナー受信機の受信信号は、ある値を中心
にゆらぎを有しており、一義的には定まるものではな
い。
【0027】このゆらぎを模擬するために本発明では、
図5の分解処理の説明図に示すように、伝播損失の厳密
解L(ri,dj)を統計処理して中心値と変動成分に分
解し、中心値を屈折に起因する伝播損失Lr(ri
j)と、変動成分を反射波の干渉に起因する伝播損失
i (ri,dj)とみなす。
【0028】即ち、ステップ24において、音源からの
距離rを固定し、深度を水面から水底まで変化させた場
合の伝播損失L(ri,dj)を中心値及び変動成分に分
解し、それぞれを屈折に起因する伝播損失Lr(ri,d
j)及び反射波の干渉に起因する伝播損失Li(ri
j)とし、さらに反射波の干渉に起因する伝播損失Li
(ri,dj)の標準偏差σ(ri,dj)を計算する。
【0029】ステップ25において、格子点(ri
j)毎に伝播時間τ(ri,dj)並びに屈折に起因す
る伝播損失Lr (ri,dj)及び反射波の干渉に起因す
る伝播損失Li (ri,dj)の標準偏差σ(ri,dj
を記憶してこのルーチンを終了する。
【0030】図6は水中音響伝播特性計算ルーチンの機
能線図であって、音線解析によって伝播経路及び伝播時
間が計算されて、そしてノーマルモード解析によって計
算された伝播損失の厳密解は統計処理されて屈折に起因
する伝播損失と干渉に起因する伝播損失の標準偏差に分
解されてメモリ102に設けられた模擬情報テーブルに
記憶される。そして、このテーブルからソナーの受信機
の位置に対応した座標に記憶されている伝播時間、屈折
に起因する伝播損失及び干渉に起因する伝播損失の標準
偏差を読み出すことにより、受信信号を決定することが
できる。
【0031】図7は水中音響伝播模擬ルーチンのフロー
チャートであって、ステップ71で制御卓11から受信
位置座標を入力する。
【0032】ステップ72ではメモリ102内の模擬情
報テーブルから受信位置座標に対応した伝播時間、屈折
に起因する伝播損失及び干渉に起因する伝播損失の標準
偏差を読み出し、ステップ73でソナー受信器の受信信
号の中心値及び標準偏差を計算する。そして、ステップ
74で受信信号を出力してこのルーチンを終了する。
【0033】
【発明の効果】第一の発明に係る水中音響伝播特性模擬
装置、方法及びプログラムによれば、伝播経路及び伝播
時間を音線法により計算し、伝播損失をノーマルモード
法により計算することにより、少ない計算量で高い精度
で伝播特性を模擬することが可能となる。
【0034】第二の発明に係る水中音響伝播特性模擬装
置、方法及びプログラムによれば、ノーマルモード法に
より計算された伝播損失は、統計処理により屈折に起因
する伝播損失である中心値と干渉に起因する伝播損失で
ある変動成分に分解することができるので、より現実的
な模擬を行うことが可能となる。
【0035】第三の発明に係る水中音響伝播特性模擬装
置、方法及びプログラムによれば、実際の音波周波数の
1/10〜1/100の周波数により伝播損失を計算す
ることができるので、伝播損失の計算時間を短縮するこ
とが可能となる。
【図面の簡単な説明】
【図1】本発明に係る水中音響伝播特性模擬装置の構成
図である。
【図2】水中音響伝播特性計算ルーチンのフローチャー
トである。
【図3】模擬方法の説明図である。
【図4】音波周波数を変更した場合の伝播損失の計算例
である。
【図5】音波周波数を変更した場合の伝播損失の計算結
果である。
【図6】分解処理の説明図である。
【図7】水中音響伝播模擬ルーチンのフローチャートで
ある。
【符号の説明】
10…水中音響伝播特性模擬装置 100…バス 101…CPU 102…メモリ 103、104…インターフェイス 11…制御卓 12…ソナー
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開2000−241533(JP,A) 特開 平6−308228(JP,A) 特開 平9−304527(JP,A) 特開 昭62−96878(JP,A) 特開 昭60−111928(JP,A) 特公 平1−29429(JP,B2) R.J.ユーリック著、土屋明訳、西 村実監修,水中音響の原理,日本,共立 出版,1988年12月 1日,P.121−128 (58)調査した分野(Int.Cl.7,DB名) G01S 3/80 - 3/86 G01S 5/18 - 5/30 G01S 7/52 - 7/64 G01S 15/00 - 15/96 G09B 9/56

Claims (9)

    (57)【特許請求の範囲】
  1. 【請求項1】 水質及び水中音源に関するパラメータを
    入力するパラメータ入力手段と、 前記パラメータ入力手段から入力されたパラメータに基
    づいて音線法により伝播経路及び伝播時間を計算する伝
    播経路・伝播時間計算手段と、 前記パラメータ入力手段から入力されたパラメータに基
    づいてノーマルモード法により前記伝播経路・伝播時間
    計算手段で計算された伝播経路に対応する伝播損失を計
    算する伝播損失計算手段と、 前記伝播経路・伝播時間計算手段で計算された伝播経路
    及び伝播時間と前記伝播損失計算手段で計算された伝播
    損失を記憶する記憶手段と、 受信位置に関するパラメータを入力する受信位置パラメ
    ータ入力手段と、 該水中音源から発生された音の、前記受信位置パラメー
    タ入力手段で入力された受信位置における受信信号を計
    算する受信信号計算手段を具備する水中音響伝播特性模
    擬のための装置。
  2. 【請求項2】 前記伝播損失計算手段が、 計算された伝播損失の中心値と該中心値周りの変動成分
    に分解する伝播損失分解手段と、 前記伝播損失分解手段で分解された変動成分の標準偏差
    を計算する標準偏差計算手段を具備する請求項1に記載
    の水中音響伝播特性模擬のための装置。
  3. 【請求項3】 前記伝播損失計算手段が、 該水中音源が発生する音の周波数の1/10から1/1
    00の周波数を使用して伝播損失を計算する請求項1又
    は2に記載の水中音響伝播特性模擬のための装置。
  4. 【請求項4】 水質及び水中音源に関するパラメータを
    入力するパラメータ入力段階と、 前記パラメータ入力段階で入力されたパラメータに基づ
    いて音線法により伝播経路及び伝播時間を計算する伝播
    経路・伝播時間計算段階と、 前記パラメータ入力段階で入力されたパラメータに基づ
    いてノーマルモード法により前記伝播経路・伝播時間計
    算段階で計算された伝播経路に対応する伝播損失を計算
    する伝播損失計算段階と、 前記伝播経路・伝播時間計算段階で計算された伝播経路
    及び伝播時間と前記伝播損失計算段階で計算された伝播
    損失を記憶する記憶段階と、 受信位置に関するパラメータを入力する受信位置パラメ
    ータ入力段階と、 該水中音源から発生された音の、前記受信位置パラメー
    タ入力段階で入力された受信位置における受信信号を計
    算する受信信号計算段階を具備する水中音響伝播特性模
    擬のための方法。
  5. 【請求項5】 前記伝播損失計算段階が、 計算された伝播損失の中心値と該中心値周りの変動成分
    に分解する伝播損失分解段階と、 前記伝播損失分解段階で分解された変動成分の標準偏差
    を計算する標準偏差計算段階を具備する請求項4に記載
    の水中音響伝播特性模擬のための方法。
  6. 【請求項6】 前記伝播損失計算段階が、 該水中音源が発生する音の周波数の1/10から1/1
    00の周波数を使用して伝播損失を計算する請求項4又
    は5に記載の水中音響伝播特性模擬のための方法。
  7. 【請求項7】 コンピュータを、 水質及び水中音源に関するパラメータを入力するパラメ
    ータ入力手段、 前記パラメータ入力手段から入力されたパラメータに基
    づいて音線法により伝播経路及び伝播時間を計算する伝
    播経路・伝播時間計算手段、 前記パラメータ入力手段から入力されたパラメータに基
    づいてノーマルモード法により前記伝播経路・伝播時間
    計算手段で計算された伝播経路に対応する伝播損失を計
    算する伝播損失計算手段、 前記伝播経路・伝播時間計算手段で計算された伝播経路
    及び伝播時間と前記伝播損失計算手段で計算された伝播
    損失を記憶する記憶手段、 受信位置に関するパラメータを入力する受信位置パラメ
    ータ入力手段、及び該水中音源から発生された音の、前
    記受信位置パラメータ入力手段で入力された受信位置に
    おける受信信号を計算する受信信号計算手段として機能
    させるための水中音響伝播特性模擬のためのプログラ
    ム。
  8. 【請求項8】 前記伝播損失計算手段を、 計算された伝播損失の中心値と該中心値周りの変動成分
    に分解する伝播損失分解手段、及び前記伝播損失分解手
    段で分解された変動成分の標準偏差を計算する標準偏差
    計算手段として機能させる請求項7に記載の水中音響伝
    播特性模擬のためのプログラム。
  9. 【請求項9】 前記伝播損失計算手段が、 該水中音源が発生する音の周波数の1/10から1/1
    00の周波数を使用して伝播損失を計算する請求項7又
    は8に記載の水中音響伝播特性模擬のためのプログラ
    ム。
JP2001363149A 2001-11-28 2001-11-28 水中音響伝播特性模擬のための装置、方法及びプログラム Expired - Fee Related JP3492347B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001363149A JP3492347B2 (ja) 2001-11-28 2001-11-28 水中音響伝播特性模擬のための装置、方法及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001363149A JP3492347B2 (ja) 2001-11-28 2001-11-28 水中音響伝播特性模擬のための装置、方法及びプログラム

Publications (2)

Publication Number Publication Date
JP2003161773A JP2003161773A (ja) 2003-06-06
JP3492347B2 true JP3492347B2 (ja) 2004-02-03

Family

ID=19173545

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001363149A Expired - Fee Related JP3492347B2 (ja) 2001-11-28 2001-11-28 水中音響伝播特性模擬のための装置、方法及びプログラム

Country Status (1)

Country Link
JP (1) JP3492347B2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110231778A (zh) * 2019-05-07 2019-09-13 中国科学院声学研究所 一种通用的uuv水下目标探测仿真方法及系统
CN110826166A (zh) * 2018-08-06 2020-02-21 中国科学院声学研究所 一种模型驱动的水下探测体系仿真系统

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0724366D0 (en) * 2007-12-14 2008-01-23 Univ York Environment modelling
JP5458756B2 (ja) * 2009-09-09 2014-04-02 日本電気株式会社 水中物体捜索支援装置、水中物体捜索支援方法、及びそのプログラム
CN102809744B (zh) * 2012-05-15 2015-02-25 哈尔滨工程大学 水下远场目标被动定位过程中的近场强干扰源抑制方法
KR101405247B1 (ko) 2014-02-17 2014-06-17 연세대학교 산학협력단 수중 음파 통신 채널의 경로 손실 산출 장치 및 방법
CN108107437B (zh) * 2017-09-28 2021-05-28 哈尔滨工程大学 一种利用简正波耦合干涉的海洋环境监测方法
CN112399302B (zh) * 2020-11-25 2022-11-18 维沃移动通信有限公司 一种穿戴式音频播放设备的音频播放方法、装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
R.J.ユーリック著、土屋明訳、西村実監修,水中音響の原理,日本,共立出版,1988年12月 1日,P.121−128

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110826166A (zh) * 2018-08-06 2020-02-21 中国科学院声学研究所 一种模型驱动的水下探测体系仿真系统
CN110826166B (zh) * 2018-08-06 2021-05-18 中国科学院声学研究所 一种模型驱动的水下探测体系仿真系统
CN110231778A (zh) * 2019-05-07 2019-09-13 中国科学院声学研究所 一种通用的uuv水下目标探测仿真方法及系统
CN110231778B (zh) * 2019-05-07 2020-06-16 中国科学院声学研究所 一种通用的uuv水下目标探测仿真方法及系统

Also Published As

Publication number Publication date
JP2003161773A (ja) 2003-06-06

Similar Documents

Publication Publication Date Title
JP3492347B2 (ja) 水中音響伝播特性模擬のための装置、方法及びプログラム
JP6737464B2 (ja) 試験システム、波形シミュレータ装置、試験方法およびプログラム
CN114460588B (zh) 基于主动式声学成像仪的高精度成像方法
CN117709207B (zh) 多波束测量船的测量布线设计方法、装置、设备和介质
US7002877B2 (en) Method and apparatus for active sonar performance prediction
JP2866930B2 (ja) 水中音響模擬装置
JP2005083932A (ja) 伝搬シミュレーション装置、伝搬シミュレーション方法、および伝搬シミュレーションプログラム
US20210364471A1 (en) Method for Creating an Evaluation Table for an Ultrasonic Inspection and Method for Ultrasonic Inspection
JP2007065760A (ja) シミュレーション装置及びシミュレーション方法
Marion et al. Toward a real-time simulation of ultrasound image sequences based on a 3-D set of moving scatterers
CN115616084A (zh) 一种大规模三维超声阵列数据快速仿真方法
Delabie et al. An acoustic simulation framework to support indoor positioning and data driven signal processing assessments
JP3482055B2 (ja) 高精度音線追跡装置および高精度音線追跡方法
US10712443B2 (en) Nonlinear intermodulation distance determination system
CN109188378B (zh) 一种超声多普勒回波信号仿真方法
Temsamani et al. Experimental validation for the diffraction effect in the ultrasonic field of piston transducers and its influence on absorption and dispersion measurements
Dogaru et al. An efficient sound propagation software simulator based on cellular automata
Rao et al. Simulation of ultrasound two‐dimensional array transducers using a frequency domain model
Schmitz et al. Merging software for sound reinforcement systems and for room acoustics
Sawada et al. Precise target strength pattern measurement in an indoor tank
CN113534165B (zh) 基于回波数据模拟器的多普勒计程仪算法性能验证装置
Nascimento et al. Enhancing the accuracy of parking assistant sensors with Bayesian filter
Furusawa et al. Exact simulation of fish school echoes and its applications
Pulvirenti et al. A perceptual evaluation of numerical errors in acoustic FEM simulation for sound quality applications
Tronco Gasparini et al. Development of a MATLAB environment software for simulation of ultrasonic field

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20031014

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

Free format text: PAYMENT UNTIL: 20081114

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20091114

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20101114

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20101114

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20111114

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20111114

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20121114

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20121114

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20131114

Year of fee payment: 10

LAPS Cancellation because of no payment of annual fees