(57)【特許請求の範囲】
【請求項１】
波力エネルギー変換器（ＷＥＣ）であって、
波に応答して互いに対して移動できる第１および第２の物体と、
前記第１および第２の物体の相対運動を有効エネルギーに変換するために、前記第１の物体と前記第２の物体との間に接続される動力取出（ＰＴＯ）デバイスと、
前記ＰＴＯと、前記第１および第２の物体との間の力を選択的に制御する手段を含むプログラマブル制御コンピュータと、
前記ＷＥＣ上に到来する波を検知し、前記ＷＥＣに対する前記到来する波の作用を予測する信号を生成するために、前記到来する波に相当する信号を前記プログラマブル制御コンピュータに供給するための、前記ＷＥＣの外部のセンサと、
前記第１および第２の物体の実際の移動を示す信号を生成し、前記実際の移動を示す信号を、処理のために前記プログラマブル制御コンピュータに供給するための、前記第１および第２の物体の移動に応答する手段とを備え、
前記プログラマブル制御コンピュータは、ＷＥＣによって前記波から抽出される平均エネルギーを増加させるように、前記ＰＴＯと、前記第１および第２の物体との間の前記力を調整するための前記ＰＴＯに印加される制御信号を生成するために、前記第１および第２の物体の前記実際の移動に相当する信号および前記予測信号に応答する手段を含む波力エネルギー変換器（ＷＥＣ）。
【請求項２】
前記第１の物体は、前記波と同相で移動するフロートであり、前記第２の物体は、前記波と異なる位相で移動し、前記第１および第２の物体の実際の移動を示す信号を生成するための、前記第１および第２の物体の移動に応答する前記手段は、前記第１および第２の物体の実際の相対速度を検知し、対応する信号を生成する手段を含む請求項１に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項３】
前記第２の物体は、本質的に固定して保持される請求項２に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項４】
前記センサは、前記ＷＥＣに瞬間的に衝撃を与えることを意図された到来する波を検知するように位置決めされる請求項１に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項５】
前記センサは、前記ＷＥＣの上流に配置される請求項４に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項６】
前記第１および第２の物体の前記実際の相対速度を検知し、対応する信号を生成する前記手段は、前記第１および第２の物体の位置の変化を、時間の関数として検知する手段、および、前記第１および第２の物体の相対速度を計算する手段を含む請求項２に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項７】
前記実際の相対速度を検知する前記手段は、前記第１および第２の物体の実際の相対速度を示す信号を生成することができる、前記第１および第２の物体に結合されたデバイスを含む請求項６に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項８】
前記実際の相対速度を検知する前記手段は、前記第１および第２の物体に取り付けられた(attach)機械的手段、電気機械的手段、光学手段、および電子手段のうちの少なくとも１つを含む請求項７に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項９】
前記ＰＴＯに印加される制御信号は、
ＦＰＴＯ＝［ＦＡ＋ＦＣ］＝［（ＺＰＴＯ−β）（ＶＰ）＋（β）（ＶＡＣＴ）］＝［ＦＰＴＯ（Ｐ）−β（ＶＰ−ＶＡＣＴ）］
として定義され、
式中、
ＦＰＴＯは、前記ＰＴＯと、前記第１および第２の物体との間に加えられる力を表し、
ＦＰＴＯ（Ｐ）は、前記予想される波情報に基づく、前記ＰＴＯと、前記第１および第２の物体との間の力を表し、
ＺＰＴＯは、選択された条件についての前記ＰＴＯのインピーダンスであり、
ＶＰは、前記ＷＥＣの予想速度であり、
ＶＡＣＴは、ＷＥＣ／ＰＴＯの実際の速度であり、
βは、支配的な波条件について項（ＺＰＴＯ−β）を最小にするように選択される請求項１に記載の波力エネルギー変換器（ＷＥＣ）。
【請求項１０】
波に応答して互いに対して移動できる第１および第２の物体ならびに前記物体間に結合された動力取出（ＰＴＯ）デバイスを有する波力エネルギー変換器（ＷＥＣ）の動作を制御する方法であって、
（ａ）到来する波を検知し、前記到来する波に関する予測情報を生成するステップと、
（ｂ）前記予測情報に対する前記ＷＥＣの応答に関する情報を予めプログラムされたコンピューティング装置に前記予測情報を供給するステップであって、それにより、信号ＦＰＴＯ（Ｐ）＝（ＶＰ）（ＺＰＴＯ）を含む、前記予測情報に対する前記ＷＥＣの前記応答を示す信号を生成する、供給するステップと、
（ｃ）選択された波条件について、項（ＺＰＴＯ−β）を最小にするβの特定の値を選択するステップと、
（ｄ）（ＺＰＴＯ−β）（ＶＰ）に等しい力ＦＡを生成するステップと、
（ｅ）前記ＷＥＣの前記第１の物体と前記第２の物体との間の前記実際の相対速度ＶＡＣＴを検知し、前記情報を前記コンピューティング装置に供給するステップと、
（ｆ）（β）（ＶＡＣＴ）に等しい力ＦＣを生成するステップと、
（ｇ）前記ＰＴＯと、前記ＷＥＣの前記第１および第２の物体との間に加えられた力である信号ＦＰＴＯを生成するステップと含み、
ＦＰＴＯ＝［ＦＡ＋ＦＣ］＝［（ＶＰ）（ＺＰＴＯ−β）＋（ＶＡＣＴ）（β）］＝［ＦＰＴＯ（Ｐ）−β（ＶＰ−ＶＡＣＴ）］であり、ＦＰＴＯ（Ｐ）は、前記予測波情報に基づいて、前記ＰＴＯおよび前記ＷＥＣに加えられる力であり、
ＶＰは、前記予測波に応答する前記ＷＥＣの予測速度であり、
ＺＰＴＯは、選択された条件についての前記ＰＴＯのインピーダンスである方法。
【請求項１１】
前記第１の物体は、前記波と同相で移動するフロートであり、前記第２の物体は、前記波と異なる位相で移動し、前記第１の物体と前記第２の物体との間に接続される動力取出（ＰＴＯ）デバイスは、前記第１および第２の物体の相対運動を有効エネルギーに変換する請求項１０に記載の方法。
【請求項１２】
前記第２の物体は、本質的に固定して保持される請求項１１に記載の方法。
【請求項１３】
到来する波を検知するセンサが、前記ＷＥＣの上流に配置される請求項１０に記載の方法。
【請求項１４】
前記ＷＥＣの前記第１の物体と第２の物体との間の前記実際の相対速度ＶＡＣＴは、前記２つの物体間に接続されたデバイスによって検知される請求項１０に記載の方法。
【請求項１５】
波力エネルギー変換器（ＷＥＣ）システムであって、
波に応答して互いに対して移動できる第１および第２の物体と、
前記第１および第２の物体の相対運動を有効エネルギーに変換するために、前記第１の物体と前記第２の物体との間に接続される動力取出（ＰＴＯ）デバイスと、
前記第１および第２の物体の前記相対位置および前記相対運動のうちの少なくとも１つを制御するための力を、前記ＰＴＯが前記ＷＥＣに加えることを選択的に可能にするプログラマブル制御コンピュータと、
前記ＷＥＣ上に到来する波を検知するための、前記ＷＥＣの外部のセンサ、および、前記ＷＥＣに対する前記到来する波の作用を予測する信号を生成するために、前記到来する波に関する情報を前記プログラマブル制御コンピュータに供給するための、前記検知された波に応答する手段と、
前記第１および第２の物体の実際の相対速度を示す信号を生成し、前記信号を、処理のために前記プログラマブル制御コンピュータに供給するための、前記第１および第２の物体の移動に応答する手段とを備え、
前記制御コンピュータは、前記ＷＥＣによって前記波から抽出される平均エネルギーを増加させる傾向がある前記ＰＴＯに対する制御信号を生成するために、前記第１および第２の物体の前記実際の相対速度に相当する信号および前記予測信号に応答するようにプログラムされる波力エネルギー変換器（ＷＥＣ）システム。
【請求項１６】
波力エネルギー変換器（ＷＥＣ）システムであって、
波に応答して互いに対して移動できる第１および第２の物体と、
前記第１および第２の物体の相対運動を有効エネルギーに変換するために、前記第１の物体と前記第２の物体との間に接続され、同様に、前記ＷＥＣに対して制御用力を加えることが可能な動力取出（ＰＴＯ）デバイスと、
前記ＰＴＯに結合され、前記第１および第２の物体の前記相対位置および前記相対運動のうちの少なくとも１つを制御するための力を、前記ＷＥＣが前記ＷＥＣに加えることを選択的に可能にする手段を含むプログラマブル制御コンピュータと、
到来する波および前記ＷＥＣに対する前記到来する波の予測衝撃に応答し、それにより、前記予測衝撃に相当する信号を前記プログラマブル制御コンピュータに供給する、第１の手段、および、前記ＷＥＣの前記実際の応答を反映する条件に応答し、それにより、前記予測信号値と、前記ＷＥＣが受ける前記実際の値との差を補償するための実際の条件に相当する信号を前記プログラマブル制御コンピュータに供給し、また、前記ＰＴＯに印加される対応する制御信号を生成する、第２の手段とを備える波力エネルギー変換器（ＷＥＣ）システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、その内容が参照により本明細書に組込まれる高度波力エネルギー変換器の制御についての２００８年５月３０日に出願された仮出願シリアルナンバー第６１／１３０，５２９号からの優先権を主張する。
【０００２】
本発明は、エネルギーを生成するときの波力エネルギー変換器（ＷＥＣ）の効率を上げるために、ＷＥＣの応答を制御する装置および方法に関する。本発明はまた、波力エネルギー変換デバイス用の制御システムに関し、より詳細には、波エネルギーをより効率的に捕捉する制御アルゴリズムに関する。
【背景技術】
【０００３】
対象となる波力エネルギー変換器（ＷＥＣ）は、一般に、水塊内に設置され、かつ、波の力にさらされると、互いに対して移動できる第１および第２の物体を含む。通常、フロートまたはシェルと呼ばれる一方の物体は、波と同相で移動するように設計され、支柱またはスパーと呼ばれる他の物体は、波と異なる位相で移動する（２物体型システムをもたらす）か、または、相対的固定して保持される（１物体型システムと呼ばれるものをもたらす）ように設計される。パワー変換メカニズム（または、動力取出、ＰＴＯ、デバイス）は、２つの物体間に接続されて、２つの物体の相対運動をある形態の有効エネルギー（たとえば、電気）に変換する。通常、対象となるＷＥＣは、動力取出（ＰＴＯ）デバイスの動作を調整するための、ＰＴＯの一部であるか、または、ＰＴＯから分離したコントローラおよびデータ・プロセッサを含む。
【０００４】
本発明を実施するときに使用するためのＰＴＯは、第１および第２の物体の相対運動を有効エネルギー（たとえば、電力）に変換することが可能な、ＷＥＣの第１の物体と第２の物体との間に接続された任意の適したデバイスまたはコンポーネント（たとえば、モータ／発電機、リニア発電機、ラックおよびピニオン、あるいは任意の他のメカニズムまたはリンク機構）を含んでもよい。動力取出（ＰＴＯ）の性能および機能は、効率的な方式で波パワーを利用することにとって重要である。一般に、ＷＥＣのＰＴＯの動的挙動は、環境因子と制御因子の組合せによって支配される。環境因子は、ＷＥＣの物体に作用する波誘導力および風誘導力を含む。制御因子は、ＷＥＣの第１の物体および第２の物体に、また、両者の間に加えられる力（複数可）を調整するために使用されてもよいＷＥＣのコントローラおよびデータ・プロセッサによって支配される。
【０００５】
ＷＥＣシステムは、ＷＥＣおよびＰＴＯの応答を予測するために使用されうる数学的表現を生成するためにモデル化されうる。異なるＷＥＣおよびＰＴＯシステムは、異なる方法で応答する可能性がある。特定のＷＥＣおよびそのＰＴＯが種々の刺激に応答する方法は予測されうる。そのため、２つの物体（たとえば、フロートおよびスパー）の相対運動に対するＰＴＯの応答は、定義されるかまたはモデル化されうる。同様に、ＰＴＯに供給されるエネルギーに対するＷＥＣの応答（移動）もまた定義されるかまたはモデル化されうる。モデルが良ければ良いほど、ＷＥＣおよびＰＴＯについての数学的に予測される応答が正確である。数学的予測（複数可）は、システムの効率を上げるために、ＰＴＯの同調および／またはＰＴＯのポジショニングおよび／またはＰＴＯへの／からの電力の注入／抽出を制御するために使用されてもよい。これは、ＷＥＣのコンポーネント（たとえば、フロート、スパー、ＰＴＯ）に、また、それらの間に結合されたコントローラおよびデータ・プロセッサの組込みを仮定する。波の移動に応答して互いに対して移動する２つの物体間に接続され、２つの物体の相対運動を有効エネルギーに変換するＰＴＯの挙動は、効率的な方式で波パワーを利用するときに重要である。
【０００６】
通常、電力は、ＰＴＯから抽出される。しかし、時として、２つの物体（フロートおよびスパー）の相対運動から得られる平均電力を最適化するために、ＰＴＯに電力を供給し、２つの物体の相対運動をＰＴＯに制御させるかまたは調整させることが望ましい。先に述べたように、２つの物体は、移動するように設計されてもよく、または、他の物体が固定して保持された状態で、一方の物体だけが移動するように設計されてもよい。
【０００７】
ＷＥＣによって生成される電力を最適化するために、ＷＥＣシステムは、到来する波の選択された特性を検知し、到来する波についての「予測（predictive）」情報をＷＥＣに提供するために、ＷＥＣからある距離のところで「上流に（upstream）」配置された波センサを含んでもよい。移動する浮遊物体の特性のために、最適なパワー捕捉には、数秒後に到達する波に対するＷＥＣおよびＰＴＯの動きが必要となるため、予測情報が望まれる。そのため、対象となるシステムは、ＷＥＣに衝突することが予測される波の性質を確定するために、到来する波の振幅（高さ）および周波数などの選択された特性を検知する波検知（一般に、ＷＥＣの外部）および予測サブシステムを含む。波予測サブシステムは、その信号を制御コンピュータに供給し、制御コンピュータは、その後、ＷＥＣの平均電力出力を最適化しようとしてＰＴＯを制御するまたは調整する（ＷＥＣから電力を供給すること、および／または、抽出することを含んでもよい）。
【０００８】
しかし、予測される波ΨＰが、ＷＥＣに衝突する実際の波ΨＡと異なる可能性がある点で、問題が存在する。差は、いくつかの異なる因子による可能性がある。１つの因子は、波予測サブシステムの波検出サブシステムが、サンプリングの不正確さを受け、予測誤差をもたらすことである。同様に、ＷＥＣに入る実際の波を、予測される波と異なるようにさせる、多くの異なる流れおよび大気条件が存在する。予測される波ΨＰと実際の波ΨＡとの差は、「波誤差（wave error）」として定義されてもよい。予測される波ΨＰと異なる実際の波ΨＡが、ＷＥＣに衝突する場合、また、衝突するとき、ＷＥＣコンポーネント（たとえば、ＰＴＯおよび浮遊物体）は、衝突する波に対する応答を最大にするための、その最適位置に位置決めされていない可能性がある（また、実際には、最適位置からかなり遠くに位置決めされている可能性がある）。その結果として、ＷＥＣの応答が効率的でなくなる。
【０００９】
知られている従来技術の方法では、予測される波に対応する情報および予測される波がＷＥＣに衝突するという仮定に対して完全な信頼が置かれる。これは、ＷＥＣに実際に衝突する波が、予測される波と異なるときに著しい誤差をもたらす。したがって、従来技術のスキームは、満足できるものではない。さらに、十分に正確な波予測を提供しうる経済的に実行可能な、知られているシステムが存在しない。そのため、知られている予測システムに関する問題が、波誤差が存在するときに（たとえば、ＷＥＣに実際に衝突する波が、予測される波と異なるときに）起こる。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】仮出願シリアルナンバー第６１／１３０，５２９号
【発明の概要】
【課題を解決するための手段】
【００１１】
出願人の発明は、一部には、波力エネルギー変換器(wave energy converter)（ＷＥＣ）に衝突する「実際の（actual）」波（ＷＥＣによって見られる実際の力を決定する）が、期待される、または、「予測される（predicted）」波と著しく異なる可能性があるため、到来する波に関連する予測情報にだけ頼ることは、最適な結果をもたらさないという認識に存在する。そのため、ＷＥＣコンポーネントを位置決めし移動させるために、予測情報にだけ頼ることは、ＷＥＣコンポーネントが、その最適位置（複数可）に位置決めされず、最適運動を受けない可能性があるため、「波誤差(wave error)」を補償せず、最適以下の結果をもたらす。
【００１２】
出願人の発明は、ＷＥＣを制御するために、「予測（predictive）」信号および「実際の」信号の組合せを使用する方法および装置を含む。信号の組合せは、一部には、波センサから得られる予測情報への、また、一部には、ＷＥＣのコンポーネントの実際の応答（たとえば、相対速度）への信頼を含む。これらの信号の適切な選択および組合せは、誤差の小さい改良型制御システムおよびより大きな電力効率をもたらすＷＥＣシステムをもたらす。
【００１３】
本発明を具現化するＷＥＣシステムは、到来する波の振幅（および／または、周波数などの他の選択された特性）を検知するための、一般にＷＥＣの外部にあるセンサ、および、データ処理能力を含む制御コンピュータに、対応する「予測」信号を供給する波検知および応答性サブシステムを含む。ＷＥＣの選択された実際の条件（たとえば、速度）を表す信号はまた、制御コンピュータに供給され、制御コンピュータは、ＰＴＯによって捕捉され、かつ、負荷に利用可能にされる平均波パワーが最大になるように、ＷＥＣのコンポーネントに対して適切な信号（力）を生成するために、予測される情報と実際の情報を処理するようにプログラムされる。
【００１４】
制御コンピュータのプログラミングは、以下で説明されるように、ＷＥＣの数学的モデルおよび種々の条件に対するその応答に基づいてもよい。本発明によれば、ＷＥＣの数学的モデリングに基づいて、（電力がＷＥＣによって生成されているときに）ＷＥＣによってＰＴＯに加えられる力である力Ｆ１またはＦＰＴＯ、および／または、（電力がＷＥＣによって消費されているときに）ＰＴＯによってＷＥＣに加えられる力は、
ＦＰＴＯ＝ＦＰＴＯ（Ｐ）−β（ＶＰ−ＶＡＣＴ）        式１
として確定される。
式中、
（ａ）ＦＰＴＯ（Ｐ）は、予測される波に基づいて、ＷＥＣによってＰＴＯに加えられる力、および／または、ＰＴＯによってＷＥＣに加えられる力であり、
（ｂ）ＶＰは、予測される波に応答するＷＥＣの予測される速度であり、
（ｃ）ＶＡＣＴは、ＷＥＣを構成する複数の物体（または、１つの物体）の実際の（瞬時の）相対運動を反映する選択された１つの点（または、複数の点）を検知することによって確定されうるＷＥＣの「実際の」または瞬時の速度であり、
（ｄ）βは、
ｉ−ＰＴＯの選択された複素機械インピーダンスであるＺＰＴＯ、および、
ｉｉ−βが、ωに関して一定であるか、または、ωに関して因果関数であるという制約を受けて、ＺＰＴＯの選択された条件（たとえば、支配的な波周波数についてのＺＰＴＯの値）について、項（ＺＰＴＯ−β）（ＶＰ）を最小にすることになるβの値
を計算することによって確定される。
【００１５】
例によれば、本発明を具現化する方法は、
（ｉ）ＷＥＣ幾何形状および支配的な波条件に固有であるΓβ（ω）と呼ばれる数学的関数を得ること、
（ｉｉ）波予測サブシステムから波予測信号を得ること、
（ｉｉｉ）実際のＰＴＯ／ＷＥＣの速度および／または位置の瞬時測定値を得ること、
（ｉｖ）所望のＰＴＯ力（因果的力、ＦＣ）のある部分を、実際のＰＴＯの速度に定数（β）を乗じた値に比例するものとして計算することであって、βの値は、ＺＰＴＯの選択された条件（たとえば、支配的な波周波数についてのＺＰＴＯの値）について、項（ＺＰＴＯ−β）を最小にするように選択され、ＺＰＴＯは、図９（ＷＥＣの特定の例についての複素ＰＴＯインピーダンスの実数部および虚数部を示す）に示す特性を有してもよい、計算すること、
（ｖ）所望のＰＴＯ力（因果的力、ＦＡ）の第２の部分を、上記ステップ（ｉ）からのΓβ（ω）関数およびステップ（ｉｉ）からの波予測に応じた、うまく定義された数学的演算の結果であるとして計算すること、
（ｖｉ）２つのＰＴＯ力を加算することであって、それにより、所望の全体の力（Ｆ１＝ＦＰＴＯ＝ＦＡ＋ＦＣ）を得る、加算すること、および、
（ｖｉｉ）その力をＰＴＯに加えること
を含む。
【００１６】
添付図面では、同じ参照文字は同じコンポーネントを示す。
【図面の簡単な説明】
【００１７】
【図１】本発明を実施するのに適するＷＥＣの略ブロック図である。
【図２】本発明による「外側ループ（outer loop）」制御システムのフローチャートである。
【図３】本発明を実施するときに使用するための「波予測（wave prediction）」システムのフローチャートである。
【図４】本発明によるＦｐｔｏを計算する「内側ループ（inner−loop）」制御システムのフローチャートである。
【図５】本発明の一態様による情報の処理および制御を示すフロー図である。
【図６】ＷＥＣの実際の速度を得るための例証的な図である。
【図７】ＷＥＣの実際の速度を得るための例証的な図である。
【図８】本発明および従来技術による、時間領域の最適な力応答伝達関数の波形状表現の図である。
【図９】対象となるＷＥＣ幾何形状についての、最適ＰＴＯインピーダンス（ＺＰＴＯ）の実数部および虚数部の関数形態を示す図である。
【図１０】βの所与の値についての、関数Γ０（ω）およびΓβ（ω）の周波数領域のグラフである。
【図１１】ＷＥＣシステムに本発明の概念を適用することの利益を示すことを意図された例証的な波形図である。
【図１２】ＷＥＣシステムに本発明の概念を適用することの利益を示すことを意図された例証的な波形図である。
【図１３】ＷＥＣシステムに本発明の概念を適用することの利益を示すことを意図された例証的な波形図である。
【図１４】ＷＥＣシステムに本発明の概念を適用することの利益を示すことを意図された例証的な波形図である。
【図１５】波誤差に関して、従来技術の方法に比べて、本発明の方法を使用して改善された特性を示す、数値モデルの結果を示すグラフである。
【発明を実施するための形態】
【００１８】
図１は、本発明を実施するときに使用するための波力エネルギー変換器（ＷＥＣ）システムの主要なコンポーネントを示す。ＷＥＣは、第１の物体またはフロート１００および第２の物体またはスパー２００を含む。第１の物体（たとえば、フロート）は、概して波と同相で移動するように設計され、第２の物体（たとえば、スパー）は、概して波と異なる位相で移動するか、または、機械的グラウンドに相対的に固定されて保持されるように設計される。動力取出（ＰＴＯ）デバイス３００は、第１の物体と第２の物体との間に結合され、第１および第２の物体の相対運動を有効エネルギーに変換し、第１および第２の物体から動力を受取るか、または、第１および第２の物体へ動力を伝達することが可能である。ＰＴＯは、この機能を実施しうる任意の適したデバイスであってよい。２つの物体は、ＰＴＯと共に、波力エネルギー変換器すなわちＷＥＣを構成する。
【００１９】
到来する波６００に応答する、上流の波検出システム５００からの信号は、波予測サブシステム７００に供給され、波予測サブシステム７００の信号は、コンピュータ・コントローラ４００と共に機能するデータ・プロセッサ４０２に供給される。コンピュータ／コントローラ４００は、プログラミング能力および拡張型データ格納能力４０４を含んでもよく、それらは全て、本発明を実施するために、ＷＥＣに関するデータを格納するため、また、信号を生成するための種々の計算を補助するためにプロセッサ４０２にとってアクセス可能である。たとえば、ＷＥＣシステムに実際に衝突する波の高さに関するデータが格納される。サブシステム７００は、コンピュータ４００と独立であってよく、または、コンピュータ４００の一部であってよい。実際の波センサは、入ってくる波の検知を可能にする任意の好都合な地点に配置されてもよい。波検出システム５００および波予測サブシステム７００は、到来する波（複数可）を解析して、対応する期待されるまたは予測される速度（ＶＰ）および予測情報に基づいてＷＥＣ／ＰＴＯが受けるべきである力［ＦＰＴＯ（Ｐ）］を確定し、対応する信号をプロセッサ４０２に供給するように機能する。
【００２０】
ＷＥＣの「実際の」速度（すなわち、第１および第２の物体の相対運動）を示す信号もまた、コントローラ４００に適用される。ＷＥＣ（またはＰＴＯ）の実際のまたは瞬時の速度は、フロートおよびスパーの相対速度を指す。すなわち、フロートおよびスパーは共に、互いに対して移動してもよく、または、一方が固定され、他方が移動してもよい。速度は、直線または回転運動を含む、ＷＥＣを構成する物体の任意の相対運動を指してもよい。図１は、第１および第２の物体（１００、２００）に結合され、ＰＴＯおよびコンピュータ／プロセッサ４００、４０２に結合された出力を有する位置または速度センサ８００を示す。センサ８００は、物体１００、２００の位置を検知するために使用されてもよく、その信号は、第１および第２の物体の実際の相対速度（ＶＡＣＴ）を計算するために（デバイス８００またはＰＴＯまたはプロセッサによって）使用される。あるいは、実際の相対速度は、ＰＴＯ内に含まれるコンポーネントから直接得られてもよい。ＶＡＣＴは、ＷＥＣまたはＰＴＯのＶＡＣＴを指すために交換可能に使用される。
【００２１】
そのため、コンピュータ・コントローラ４００およびプロセッサ４０２は、データベースおよびデータ処理能力を含み、ＰＴＯ３００に制御信号を選択的に供給するために、到来する（「予測される」）データおよび「実際の」データを処理するようにプログラムされる。本発明によれば、ＷＥＣ内に存在する（「予測される」と対照的な）「実際の」または「瞬時の（instantaneous）」条件（たとえば、速度）は、（たとえば、内部または外部検知メカニズムによって）検知され、「予測される」波（複数可）に関する情報を処理するように同様に機能するコンピュータ・コントローラ４００に給送されて、ＷＥＣ／ＰＴＯに対して、これらの信号に対応する力が生成される。
【００２２】
先に説明したように、「予測される」信号は、いくつかの因子のために誤っている可能性がある。そのため、予測される情報は、パワー捕捉を最適化するための所望の条件にＷＥＣコンポーネントをセットアップするために、非常に望ましくかつ必要であるが、予測される情報に関連する誤差がほぼ常に存在することになること、および、これらの誤差の影響が最小にされるべきであることを、出願人は認識した。出願人の発明は、一部には、ＷＥＣ／ＰＴＯの実際の条件を検知すること、および、実際の情報を検知された情報と組合せることが、より最適な動作を生成し、より大量の平均パワー捕捉およびより大きな効率をもたらすために使用されうることに存在する。
【００２３】
本発明は、幅広いクラスの波力エネルギー変換器と共に使用するために適用可能であり、幅広いクラスの波力エネルギー変換器は、波に応答する、かつ／または、海底に対して強固な接続を有してもよい、かつ／または、ヒーブで、または、別の自由度（サージ、スウェイ、ピッチ、ロール、ヨー）で応答してもよく、または、これらの自由度の任意の組合せで応答してもよく；かつ／または、振動性水柱または水はしけなどのデバイスを含んでもよい、１つまたは複数の物体を有するＷＥＣを含む。本発明はまた、直接リンク機構あるいは空気圧または水圧リンク機構を含んでもよい幅広いクラスの動力取出デバイスと共に使用するために適用可能である。
【００２４】
説明を容易にするために、本発明は、図１に示すように、（ａ）概して波と同相で移動するように設計された第１の物体と、（ｂ）概して波と異なる位相で移動するように設計された第２の物体と、ここで、これらの２つの物体のうちの少なくとも一方は一般に流体力学的波力に応答性がある、（ｃ）２つの物体間に結合されたＰＴＯであって、（ｉ）一般に２つの物体の相対運動を有効エネルギーに変換し、（ｉｉ）ＷＥＣの効率の増加を達成するために、２つの物体が互いに対して移動するように選択的にさせるデバイスと見なされてもよい、ＰＴＯとを有するＷＥＣを参照して述べられる。
【００２５】
本発明をよりよく理解するために、数値モデリングに適用しうるＷＥＣシステムの簡略化した数学的表現が、以下の式２および式３に示される。図１を参照して、これらの式は、スパー２００が固定して保持され、フロート１００が概して波と同相で移動する傾向があると仮定して得られる。
【００２６】
【数１】
0005474056.tif 000002
式中、
ＸＷＥＣは、スパーに対するフロートの位置を示し、
Ｆｗａｖｅは、フロートに作用する励起力を指し、
ＦＰＴＯは、ＰＴＯによってフロートおよびスパーに加えられる力を指し、
Ｂは、ＰＴＯ力と、フロートとスパーの相対運動との間の抵抗性関係を指定するある比例定数である。
【００２７】
式２は、単一物体波力エネルギー変換器についてのダイナミック運動方程式である。ＷＥＣに対する２つの力、波の力（ＦＷＡＶＥ）および動力取出デバイスの力（ＦＰＴＯ）が存在する。ダイナミック運動方程式（式２）の解は、ＦＰＴＯ（式３）の数学的表現の知識を必要とする。
【００２８】
式２および式３を周波数領域で書き直すことが望ましい。
ＶＰＴＯ（ＺＩ＋ＺＰＴＯ）＝ＦＷＡＶＥ                      式４
【００２９】
式２および式３から式４を導出するために、ダイナミック変数ＸＷＥＣは、ＰＴＯの速度（ＶＰＴＯ）で置換され、慣性項
【数２】
0005474056.tif 000003
は、より一般的なＶＰＴＯＺＩで置換され、ここで、ＺＩは固有の機械インピーダンスであり、式３は、以下のように書き直される。
ＦＰＴＯ＝ＺＰＴＯＶＰＴＯ                                式５
【００３０】
式４内の全ての量は、複素数でかつ周波数依存性がある可能性がある。（複素数が物理量を表すのに使用されるとき、複素数の実数部が、物理量を表すと考えられる。）
【００３１】
一般に、ＺＰＴＯは、波とＰＴＯとの間のパワー伝達を最大にするように選択されることに留意されたい。さらに、最も一般的なクラスのＺＰＴＯ値の場合、式５におけるＦｐｔｏの計算は、波の予測知識を必要とする。ＰＴＯインピーダンスＺＰＴＯが定義されると、式４は、知られている数学的方法によって解かれてれもよい。速度ＶＰＴＯは、以下のように、ＦＷＡＶＥ、ＺＩ、およびＺＰＴＯの関数として表されてもよい。
【数３】
0005474056.tif 000004
式６で表されるＶＰＴＯについての解を用いて、ＦＰＴＯについて解くことが望ましい。
【数４】
0005474056.tif 000005
予測情報に基づいて、ＰＴＯ力が、
ＦＰＴＯ（Ｐ）＝ＺＰＴＯＶＰ                            式８
として表されてもよいことを、出願人は認識した。
式中、
（ｉ）ＦＰＴＯ（Ｐ）は、予測される波に基づく、ＷＥＣによってＰＴＯに加えられる力、および／または、ＰＴＯによってＷＥＣに加えられる力であり、
（ｉｉ）ＶＰは、ＷＥＣの予測される速度であり、
（ｉｉｉ）ＰＴＯのＺＰＴＯは、パワー伝達を最適化するように選択される。
ＶＰは、それ自体、ＺＰＴＯの関数であるため、式８は、数学的に複雑であることが留意されるべきである。しかし、ＦＰＴＯのこの表現は、本発明の性質を明確にするのに有用である。
ＦＰＴＯ（Ｐ）についての式８は、
ＦＰＴＯ（Ｐ）＝ＺＰＴＯＶＰ＝ＺＰＴＯＶＰ−βＶＰ＋βＶＰ＝（ＺＰＴＯ−β）ＶＰ＋βＶＰ      式９
として書き直されてもよいことを、出願人はさらに認識した。上記式９において、βが一定の制約を満たす（たとえば、βが因果作用素である）場合、（ＺＰＴＯ−β）ＶＰ項は、非因果的力（ＦＡ）を表し、βＶＰは、因果的力（ＦＣ）を表す。ここで、（ａ）本明細書で使用される非因果的力（ＦＡ）は、確定される（計算される）ために、ある程度の予測情報（たとえば、将来のある時間における波の高さ）を必要とする力として定義されてもよく、（ｂ）本明細書で使用される因果的力（ＦＣ）は、予測情報を全く使用することなく計算される可能性がある力として定義されてもよい。さらに、βは因果作用素（たとえば、定数）であり、因果作用素は、定義上、ＰＴＯの速度（ＶＰ）などの量によって乗算されるために、将来の（予測される）知識を全く必要としないことを、出願人は認識した。予測知識が全く必要でないため、ＶＰをＶａｃｔｕａｌで置換することが可能である。したがって、因果的力ＦＣ＝βＶＰは、βＶＡＣＴＵＡＬで置換されてもよく、誤差が減少した状態で、必要とされる力ＦＰＴＯは、
ＦＰＴＯ（Ｎ）＝ＦＩ＝（ＺＰＴＯ−β）ＶＰ＋βＶＡＣＴＵＡＬ        式１０
として表されてもよい。
式中、ＦＰＴＯ（Ｎ）（またはＦＰＴＯ）は、本発明による、ＷＥＣによってＰＴＯに加えられる力、および／または、ＰＴＯによってＷＥＣに加えられる力である。
【００３２】
式１０では、因果的力ＦＣは、βＶＡＣＴを計算することによって確定され、非因果的力（ＦＡ）は、（ＺＰＴＯ−β）ＶＰとして確定される。ここで、（ａ）ＶＡＣＴは、ＷＥＣの実際の速度であり、（ｂ）（ＺＰＴＯ−β）が最小になるため、非因果的力（ＦＡ）が、波誤差に対して比較的感度が低くなるようにさせるβの値が選択されうる。
【００３３】
ＦＰＴＯ（Ｎ）が、
ＦＰＴＯ（Ｎ）＝（ＶＰ）（ＺＰＴＯ）−β（ＶＰ−ＶＡＣＴ）      式１１
として書き直されてもよいことに留意されたい。
【００３４】
そのため、式１１は、波誤差が存在しない場合、ＶＰ＝ＶＡＣＴおよび項β（ＶＰ−ＶＡＣＴ）は、ゼロになり、予測される値の補正についての必要性が存在しない。波誤差が存在する場合、本発明を実施することは、波誤差を補償し、パワー捕捉の損失を低減する。
【００３５】
ＷＥＣを動作させるためのこれらの制御信号の実施態様は、波誤差が存在するときに、著しく大きな量の平均パワーの生成をもたらす。単一周波数の場合、生成される平均パワーＰＡＶＧは、
【数５】
0005474056.tif 000006
として表されてもよい。式中、
【数６】
0005474056.tif 000007
は、ＦＡの複素共役を示し、Ｒｅ［］は、括弧内の量の実数部を示す。
【００３６】
本発明の実施態様は、図５を参照して説明されてもよい。図５のブロック７０１〜７１３は、プロセッサ４０２およびプログラマブル・コントローラ４００の一部であってよいことに留意されたい。検知された波高さ情報は、ＷＥＣのモデリングに関する情報を含み、波高さの変化に対するＷＥＣの反応を計算することが可能なデータベースおよびプロセッサ７０１に適用される。プロセッサ７０１は、その後、（ａ）ＰＴＯ力信号ＦＰＴＯ（Ｐ）、および、（ｂ）ＷＥＣの期待される／予測される速度を示す信号ＶＰを生成しうる。ＰＴＯのインピーダンス（ＺＰＴＯ）は、ＷＥＣの種々の条件について確定され、波およびＺＰＴＯに相当するデータは、データベース・セクション７０３（実際には、７０１の一部であってよい）または任意の他の適したサイトに格納される。予測される速度信号（ＶＰ）は、その後、信号（ＺＰＴＯ）ＶＰを生成するために、ソフトウェア／ハードウェア７０３によって処理されてもよい。信号（ＺＰＴＯ）ＶＰは、その後、積（ＺＰＴＯ−β）（ＶＰ）を生成するために、ソフトウェア／ハードウェア７０５によって修正されうる。βの適切な値は、項（ＺＰＴＯ−β）（ＶＰ）を最小にするように選択されてもよい。たとえば、波の支配的な周波数においてＺＰＴＯの値に近く、信号（ＺＰＴＯ−β）（ＶＰ）を最小にするように機能するβの値が選択されてもよい。βの値またはβの任意の他の値を最小にすることは、所定の条件に応じて、使用するために選択されてもよいことに留意されたい。本発明の好ましい実施形態では、βは、選択された間隔にわたって定数であり、一方、ＺＰＴＯは、波周波数の関数として変わる。そのため、βの最良値は、波の支配的な周波数においてＺＰＴＯを消去する値である。βの選択された値について、積（ＶＰ）（β）は、処理ステップ（ボックス７０７を参照）によって得られうる。積（ＶＰ）（β）は、ＦＰＴＯの予測される因果的力の部分を表すと考えられうる。次のステップは、ＦＰＴＯ（Ｐ）から積（ＶＰ）（β）を引いて、ＦＰＴＯ（Ｐ）の非因果的成分ＦＡ、ＦＡ＝［ＦＰＴＯ（Ｐ）−（ＶＰ）（β）］を、ボックス７０９によって示す処理ステップからの出力として得ることである（ボックス７０９を参照）。ＰＴＯ力の実際の因果的成分ＦＣは、処理ステップ７１１で示すように、選択されたβでＷＥＣ速度の実際の値（ＶＡＣＴ）を乗じて、［ＶＡＣＴ（β）］の積を生成することによって得られる。実際の因果的力［ＦＣ＝（ＶＡＣＴ）（β）］は、その後、ボックス７１３に示すように、予め決定された非因果的力［ＦＰＴＯ（Ｐ）−（ＶＰ）（β）］に加算されて、修正されたＰＴＯ力を生成してもよく、波予測誤差が存在するときに、より大きな平均パワーの生成をもたらすことになる。図に示すように、ボックス７１３で示す処理ステップの出力は、
｛Ｆ１＝［ＦＡ＋ＦＣ］＝［ＦＰＴＯ（Ｐ）−（ＶＰ）（β）］＋（ＶＡＣＴ）（β）＝［ＦＰＴＯ（Ｐ）−β（ＶＰ−ＶＡＣＴ）］｝
をもたらす。
【００３７】
本発明に従ってＷＥＣシステムを動作させることは、波誤差が存在するとき（通常の動作状態である）、システムによって生成される平均パワーの著しい増加をもたらす。
【００３８】
例によれば、本発明を実施するために使用される所望の値を計算するための種々のステップは、図２、３、および４に示される。
【００３９】
図２は、検知された到来する波が、その特性およびその状態の任意の変化を確定するために処理されること（ボックス２１を参照）、および、このサンプリング・プロセスが、選択された間隔で起こること（ボックス２９を参照）を示す簡略化したフローチャートである。波の場に存在する周波数範囲が評価される（ボックス２３を参照）。波の位相（示さず）が解析されてもよい。プロセスにおける別のステップ（ボックス２５を参照）は、波の場の周波数範囲にわたって、式２６（以下を参照）において周波数領域で関数Γβ（ω）を最小にするβを選択することである。さらなるステップ（ボックス２７を参照）は、式２７に従って時間領域で関数Γβ（ｔ）を確定することを含む（以下を参照）。別のステップ（ボックス２９を参照）は、到来する波が、選択された間隔で（たとえば、５分ごとに）サンプリングされることを示す。サンプリング時間間隔は、波の変化レートに応じて、増減されてもよい。こうして、図２は、本発明を実施するために使用されてもよい「外側ループ（outer−loop）」制御システムと呼ばれてもよいものを示す。図２は、本発明が、支配的な波特性に依存する方法で最もよく実施されることを強調する。すなわち、支配的な波条件に依存する因果的および非因果的制御を確定するための手段／プロセスが含まれる。
【００４０】
図３は、検知された到来する波の振幅が、その特性およびその状態の任意の変化を確定するために処理されることを示すフローチャートである。１つのステップ（ボックス３１を参照）では、到来する波の高さが、ＷＥＣから上流の地点で測定される。波伝播のモデルに基づいて（ボックス３３を参照）、検知された到来する波に相当する、ＷＥＣにおける波の高さについての予測が行われうる（ボックス３５を参照）。情報は、その後、プロセッサで処理される（図４のボックス４５を参照）。
【００４１】
図４は、本発明を実施する基本処理ステップを示すフローチャートである。ステップは、予測される波の選択された特性の収集および確定、ならびに、因果的力（ＦＣ）および非因果的力（ＦＡ）の計算、ＦＰＴＯの計算、およびＰＴＯへのＦＰＴＯの適用を含む。図４に詳述されるように、ＦＡは、図３のボックス３５に示すステップからの出力を取得し、図４のボックス４５および４６に示すようにその情報を処理することによって得られてもよい。これは、ＦＡを生成するために、以下の式２８に示す畳込み積分から非因果的力を計算することを含む。同時に、図４は、ＰＴＯの実際の速度（ＶＡＣＴ）が得られ（ボックス４７を参照）、因果的項ＦＣを計算するためにシステムに適用されることを示す。ＦＡおよびＦＣについての情報は、加算されて（ボックス４９を参照）、本発明による、ＰＴＯに適用される力Ｆ１に相当する信号が生成される。
【００４２】
したがって、本発明を具現化する方法は、
（ｉ）ＷＥＣ幾何形状および支配的な波条件に固有の数学的関数（Γβと呼ぶ）を得ること、
（ｉｉ）波予測サブシステムから波予測信号を得ること、
（ｉｉｉ）ＰＴＯ速度の瞬時測定値を得ること、
（ｉｖ）所望のＰＴＯ力のある部分を、ＰＴＯの速度に定数（β）を乗じた値に比例するものとして計算することであって、βの値は、ＷＥＣの特定の例についての複素ＰＴＯインピーダンスの実数部および虚数部を示す図９を参照して、以下に述べるように選択される、計算すること、
（ｖ）所望のＰＴＯ力の第２の部分を、ステップ（ｉ）からのΓβ関数およびステップ（ｉｉ）からの波予測に応じた、うまく定義された数学的演算の結果であるとして計算すること、
（ｖｉ）２つのＰＴＯ力を加算することであって、それにより、所望の全体の力を得る、加算すること、および、
（ｖｉｉ）その力をＰＴＯに加えること
を含んでもよい。
【００４３】
簡単な波力エネルギー変換システムとより複雑な波力エネルギー変換システムの両方に適用可能な簡略化された数学的モデルがここで説明される。ＷＥＣのこの数学的モデルは、以下の部分、すなわち、（ａ）固有の機械インピーダンスＺＥを有する波応答性（相互作用）物体、（ｂ）その機械インピーダンスがＺＰＴＯとして選択されるＰＴＯデバイス、および、式
【数７】
0005474056.tif 000008
によって支配される流体力学的波励起力ＦＷＡＶＥに対する応答を表すと解釈されてもよい。
式中、ＶＰＴＯは波応答性物体の速度であり、ＺＥは固有の機械インピーダンスであり、ＦＷＡＶＥは、波応答性物体に作用する流体力学的波励起力である。流体力学的波励起力ＦＷＡＶＥはまた、
ＦＷＡＶＥ≡ηχＥ                                          式１３
として表されてもよい。
式中、
ηは、ＷＥＣが存在しないときに測定されることになるＷＥＣの場所における水の瞬時表面高さであり、
χＥは、波励起力応答である。波励起力応答は、水深と物体の幾何形状の関数である。式１２において、ＶＰＴＯは、ＦＷＡＶＥが流体力学的トルクを表すように回転自由度を表してもよいため、角速度を表すことになることに留意されたい。
【００４４】
ＰＴＯにおける力ＦＰＴＯは、
ＦＰＴＯ＝ＺＰＴＯＶＰＴＯ                                    式１４
に等しい。
ＰＴＯによって吸収される瞬時パワーは、
【数８】
0005474056.tif 000009
として表されうる。
式中、
【数９】
0005474056.tif 000010
は、ＦＰＴＯの複素共役である。
【００４５】
図１に示すタイプの波吸収体の場合、式１２、式１３、式１４、および式１５は、スパー２００が固定して保持される場合、または、スパー２００が移動することを許容される場合について導出されてもよい。フロート１００が波と同相で移動し、スパー２００が固定して保持される場合、機械インピーダンスＺＥは、
【数１０】
0005474056.tif 000011
として表されうる。
式中、ωは角周波数であり、Ｍは、変位した質量および流体力学的に付加された質量を含む物体１００の質量であり、Ｂは、物体１００の流体力学的放射減衰であり、Ｃは、（ｇ）（Ｓ）（ρ）として定義される流体静力学的回復力であり、ここで、ρは水の密度であり、ｇは重力加速度であり、Ｓは物体１００の水面の面積である。量Ｍ、Ｂ、およびＣは、物体の形状および水深の関数であり、適切なモデリングによって確定されてもよい。
【００４６】
スパー２００が固定して保持されない場合、式１２、式１３、式１４、および式１５は、依然として有効である。機械システムの固有インピーダンスは、もはや式１６によって指定されず、知られている技法を使用して確定されてもよい。
【００４７】
既に説明したように、ＺＰＴＯが最適パワー変換効率を提供するために、ＰＴＯインピーダンスＺＰＴＯは、
【数１１】
0005474056.tif 000012
であるように、固有インピーダンスＺＥの複素共役に等しくなるように選択されてもよい。ＺＰＴＯについてこの表現が与えられると、最適なＰＴＯ速度は、
【数１２】
0005474056.tif 000013
として表されてもよい。式中、Ｒｅ［］は、括弧内の量の実数部を示す。ＦＥ＝ηχＥであることに留意し、
【数１３】
0005474056.tif 000014
を定義する。以下の表現が、最適なＰＴＯ速度について得られてもよい。
【数１４】
0005474056.tif 000015
式２０は、
【数１５】
0005474056.tif 000016
の解が、波の高さの予備知識を必要とするため、因果的でない表現を定義することに留意されたい。数学的には、これは、この式が、畳込みとして表されなければならないからであり、関数Λが、時間の負値について非ゼロ値を有するため、ＶＰＴＯの計算は、時間の正値について（すなわち、将来において）、波の高さηの知識を必要とする。
【００４８】
波の高さが、
【数１６】
0005474056.tif 000017
が計算されるように、十分に遠い将来にわたってわかっている場合、ＰＴＯは、ＶＰＴＯが、強制的に
【数１７】
0005474056.tif 000018
になるように制御されてもよい。式２０で表されるこのタイプの制御は、ＰＴＯの速度ＶＥを制御するためにフィードバック・ループが使用されるため、フィードバック則と呼ばれてもよい。
【数１８】
0005474056.tif 000019
を生成するために、Γ０＝ΛＺＰＴＯを定義することによって、「フィードフォワード（feedforward）」則を構築することも可能である。これは、最適なＰＴＯ力について以下の表現を生成するために使用されうる。
【数１９】
0005474056.tif 000020
式２２は、式２０のように非因果的なＷＥＣについての制御則を定義する。この方法は、制御アルゴリズム計算の出力が、閉ループ制御が無い状態でＰＴＯに加えられる力であるため、フィードフォワード則と呼ばれてもよい。この方法はまた、式２２における
【数２０】
0005474056.tif 000021
の計算が、将来の時間における波の高さの知識を必要とするため、非因果的である。式２０および式２２は共に、周波数領域の乗算として書かれる。これらの式を時間領域で表現するために、基本的な作用素（式２０の場合Λ、また、式２２の場合Γ０）が、フーリエ変換によって周波数領域から時間領域へ変換されなければならず、また、その後、乗算が、畳込み積分として実行されなければならない。周波数領域におけるΓ０の明示的な表現は、以下で式２３に示される。時間領域においてΓ０についての表現をもたらすフーリエ変換演算は、以下で式２４に明示的に示される。時間領域におけるΓ０についての明示的な表現および波の高さの時系列データセットη（ｔ）によって、畳込み積分を実行し、そして、ＰＴＯに加えられる力（ＦＰＴＯ）を計算することが可能である。この畳込み積分は、式２５に明示的に示される。式２５における積分の制限が、全ての時間にわたって拡張されるべきであることが留意されるべきである。しかし、図８に示すように、Γ０関数は、比較的小さい値に減少し、積分は、したがって、時間ｔｓ−ＴΓ＋およびｔｓ＋ＴΓ−において打切られてもよい。
【数２１】
0005474056.tif 000022
Γ０（ｔ）＝∫Γ０（ω）ｅ−ｉωｔｄω                  式２４
【数２２】
0005474056.tif 000023
本発明の重要な態様は、従来技術の最適フィードフォワード力
【数２３】
0005474056.tif 000024
（式２２で定義される）が、共に加算されて最適フィードフォワード力になり、かつ、所望の特性を個々に有する２つの力（ＦＡおよびＦＣ）の和として表現されてもよい新しい最適フィードフォワード力Ｆ１によって置換されてもよいという出願人の認識である。βＶＡＣＴとして表現され、以下の式２９で定義されうる力ＦＣは、因果的である（したがって、波予測誤差に対する感度が低い）という所望の特性を有する。ＦＰＴＯ（Ｐ）−βＶＰとして表現され、以下の式２８で定義されうる非因果的力ＦＡは、（予測される波情報だけから計算される従来技術による最適ＰＴＯ力である）Ｆ０に比べて波予測誤差に対する感度が低いという所望の特性を有する。理想化された（ゼロ波誤差の）条件下で、Ｆ０とＦ１は同値である。波予測誤差が存在するという期待される条件下で、Ｆ１は、Ｆ１がＦ０に比べて波予測誤差に対する感度が低いという所望の特性を有する。結果として、Γβと呼ばれ、かつ、以下の定義
【数２４】
0005474056.tif 000025
を有する新しい関数が、Γ０関数のための置換として導入される。
【００４９】
実際には、量βは、因果的関係の一部であり、Γβ（ω）について結果として得られる表現は、対象となる周波数範囲内でできる限り小さくさせられる。図９は、特定のＷＥＣシステムについての機械インピーダンスＺＥの例を示す。最適ＰＴＯインピーダンス
【００５０】
【数２５】
0005474056.tif 000026
は、ＺＥの複素共役である。この知識によって、βは、対象となる周波数範囲内でΓβ（ω）をできる限り小さくするように選択されてもよい。βの適切な選択は、対象となる周波数範囲内でΓβがΓ０より小さいことを保証することになると結論付けられてもよい。対象となる波の周波数に応じてβの異なる値が望ましい可能性があると結論付けられてもよい。Γβ（ω）が周波数領域で定義されると、時間領域でΓβ（ｔ）を表現するためにフーリエ変換を実行することが望ましい。このフーリエ変換演算は、（式２４と類似する）式２７で定義される。ＦＡの計算をもたらす畳込み積分は、式２８に明示的に示される。
Γβ（ｔ）＝∫Γβ（ω）ｅ−ｉωｔｄω                  式２７
【００５１】
【数２６】
0005474056.tif 000027
ＦＣは、
ＦＣ＝βＶＡＣＴ（ｔｓ）                                式２９
として計算されてもよい。ＦＡとＦＣの和は、本発明に従って、ＦＰＴＯまたはＦ１をもたらす。
Ｆ１＝ＦＡ＋ＦＣ                                      式３０
本発明による、フィードフォワード制御則の規範的概要は、以下に与えられる。
１）波条件について１回実施するタスク
ａ．ＷＥＣの固有機械インピーダンスＺＥを計算する
ｂ．波励起力応答χＥを計算する
ｃ．ＰＴＯの最適機械インピーダンス
【００５２】
【数２７】
0005474056.tif 000028
を計算する
ｄ．βが一定であるか、または、βが因果的作用素であるという制約を受けて、対象となる周波数範囲にわたって（たとえば、波の周波数範囲にわたって）Γβ（ω）（式２６による）を最小にすることになるβの値を選択する
ｅ．式２６に従ってΓβ（ω）を計算する
ｆ．式２７に従って時間領域でΓβ（ｔ）を計算する
ｇ．Γβ（ｔ）の最小サポート、最大サポートＴΓ−、ＴΓ＋を確定する（図８、９、１０を参照）
２）全ての時間ステップｔｓにおいて実施するタスク
ａ．将来にわたる波の高さＴΓ−秒の予測を確定する
［ここで、ｔｓ−ＴΓ＋からｔｓ＋ＴΓ−までの時間の関数として波の高さηについての表現が存在することに留意されたい］
ｂ．畳込み積分（式２８）から非因果的力の項を計算する
ｃ．瞬時の実際のＰＴＯ速度ＶＡＣＴを確定する
ｄ．因果的力の項ＦＣ＝βＶＡＣＴ（ｔｓ）を計算する
ｅ．Ｆ１＝ＦＡ＋ＦＣを計算する
ｆ．Ｆ１は、本発明による最適なＰＴＯ力である。ＰＴＯにＦ１を加えるようにＰＴＯを制御する
図６は、速度検知手段を有するＷＥＣの略図である。波応答性物体１００および２００は、波に応答して移動し、動力取出システム（ＰＴＯ）３００の運動をもたらす。ＰＴＯは、水圧ピストン３０１、水圧チュービング３０２、および水圧モータ３０３からなる非整流型直線−回転変換メカニズムからなる。回転モータ３０３は、機械−電気変換器（発電機）３０５に関する軸３０４上で同じ場所に配置される。負荷回路３０６は、発電機３０５の両端に接続された抵抗器（負荷）３０８、および、負荷の両端に接続された電圧検知デバイス３０９を含む。１つの物体１００（たとえば、フロートおよびスパーの一方）と他の物体２００（たとえば、フロートおよびスパーの他方）との間の相対速度と、発電機３０５によって生成され、電圧センサ３０９によって検知される電圧との間に直線関係が存在する。測定電圧と、フロートとスパーとの間の相対速度との間の関係は、ＰＴＯのコンポーネントの特性の関数として確定されてもよい。ケーブル３１１は、電圧信号を制御コンピュータ４００に送信し、したがって、制御コンピュータがＰＴＯの実際の速度を確定することを可能にする。そのため、図６は、フロートとスパーとの間の相対速度を測定する手段を示す。
【００５３】
図７は、位置検知手段を有するＷＥＣの略図を示す。波応答性物体１００および２００は、波に応答して移動する。コンポーネント５５０および５５１は、位置検知メカニズムを備える。知られている位置検知メカニズムの場合、５５１は、物体１００と物体２００との間の張力下のワイヤを表す。コンポーネント５５０は、物体１００と物体２００との間の相対位置が変化するにつれて、ワイヤ５５１を取込みうる（または、繰り出しうる）張力をかけたスプールを含む。コンポーネント５５０は、繰り出されたワイヤ５５１の長さを検出する手段を含み、この情報をケーブル３１１を介して制御コンピュータ４００に送信し、したがって、制御コンピュータが、物体１００と物体２００との間の実際の距離を確定することを可能にする。制御コンピュータ４００は、この信号を識別して、ＰＴＯの実際の速度を確定してもよい。図６および７は、任意の適した実際の速度の検知手段が、本発明を実施するために使用されてもよいことを示すことを意味される。
【００５４】
図８は、力応答伝達関数Γ０（ｔ）およびΓβ（ｔ）の略図を示す。これらの関数Γ０（ｔ）およびΓβ（ｔ）は、ｔ＝０軸の近くで大きさが比較的大きくなる傾向があり、また、ｔが＋無限大または−無限大に近づくにつれて、大きさがゼロに向かう傾向がある。その結果、ｔの正値および負値が存在し、それらの値を超えると、関数Γ０（ｔ）およびΓβ（ｔ）の大きさが十分にゼロに近くなるため、実際には、ｔのこれらの値を超えて畳込み積分（式２８）を実行することは必要でない。これらの制限は、ｔの負値についてＴΓ−として、ｔの正値についてＴΓ＋として示される。
【００５５】
図９は、波力エネルギー変換器についての最適ＰＴＯインピーダンス（ＺＰＴＯ）の例を示す。実数部および虚数部は、上部プロットおよび底部プロットにそれぞれ示される。虚数部は、所与のＷＥＣ設計について０．１１Ｈｚの近傍でゼロを通過することが見てわかる。
【００５６】
図１０は、従来技術の関数Γ０（ω）、および、ベータの実数でかつ一定値についての本発明による対応するΓβ（ω）の周波数領域のグラフを示す。関数Γβ（ω）は、
【数２８】
0005474056.tif 000029
として式２６で定義される。関数Γβ（ω）は、
（ａ）そのＷＥＣについての
【数２９】
0005474056.tif 000030
のグラフが図９に示すようになるような特定の波力エネルギー変換器の特性、および、
（ｂ）約５００ｋＮｓ／ｍのβの値
に基づく。β＝５００ｋＮｓ／ｍの場合、
【数３０】
0005474056.tif 000031
の値が、約０．２Ｈｚでゼロであることが、
【数３１】
0005474056.tif 000032
のグラフ（図９）から見てわかる。同様に、関数Γβ（ω）もまた、０．２Ｈｚで最小になる。そのため、本発明の実施に関連して選択されたβのこの値は、約０．２Ｈｚの周波数の波について、波誤差に対するＷＥＣシステムの感度を最小にする。他の支配的な波周波数においてＷＥＣシステムの感度を最小にするために、βの他の値が、本発明に従って選択されてもよい。
【００５７】
図１１、１２、１３、および１４は、本発明を示し、その説明を補助するために使用される波形図である。
【００５８】
図１１は、従来技術を使用した、最適制御の中間結果を示す。この図は、ゼロ波誤差の条件に当てはまり、ゼロ波誤差の条件は、期待されるまたは予測される波がＷＥＣに衝突する波と同一である条件として定義される。
ａ−図１１の波形Ａは、ＷＥＣの場所に衝突する波の高さ（Ｅ）を示す。
ｂ−図１１の波形Ｂは、波形Ａに応答するＷＥＣ（ＷＥＣの１つまたは複数の移動可能部分）の速度を示す。
ｃ−図１１の波形Ｃは、ＷＥＣのＰＴＯのＰＴＯ力ＦＰＴＯを示す。実際には、波の高さおよびＦＰＴＯは、移動可能物体（ＷＥＣ）の速度を定義する。
ｄ−図１１の波形Ｄは、ＷＥＣによる生成電力（Ｐ）を示す。
Ｐは、速度×（ＦＰＴＯ）に等しく、速度は（ＶＡＣＴ）である。
【００５９】
図１２は、本発明を使用した、最適制御の中間結果を示す。この図もまた、ゼロ波誤差の条件に当てはまる。図１２の波形ＡおよびＢは、ＷＥＣの場所に衝突する波の高さ（Ｅ）およびＷＥＣの対応する速度をそれぞれ示す。図１２の波形Ｃは、ＷＥＣのＰＴＯのＰＴＯ力ＦＰＴＯを示し、ＦＰＴＯは、因果的力ＦＣおよび非因果的力ＦＡによって表現されてもよい。
ここで、
ＦＰＴＯ＝（ＶＰ）（ＺＰＴＯ−β）＋（ＶＰ）（β）      式Ｂ
ＦＡ＝（ＶＰ）（ＺＰＴＯ−β）                        式Ｂ１、および、
ＦＣ＝（ＶＰ）（β）                                式Ｂ２
（サブプロット３−波形Ｃの）因果的力および非因果的力が、加算されてサブプロット３における最適力になることに留意されるべきである。図１２の波形Ｄは、ＷＥＣによる生成電力（Ｐ）を示す。Ｐは、（ＦＣ＋ＦＡ）にＶＰを乗じた値に等しく、図１１のシステムの電力に等しい。因果的電力寄与および非因果的電力寄与が、加算されてサブプロット４の波形Ｄにおける最適電力プロファイルになることに留意されるべきである。
【００６０】
図１３は、波誤差の存在下で従来技術を使用した制御を示す。そのため、この図は、ＷＥＣに衝突する実際の波が、波形Ａに示すように、予測される波と異なる条件を示す。この条件の場合、ＦＰＴＯおよび実際の速度が、図１２の条件の場合と異なるであろうことが立証される。そのため、予測される速度ＶＰは、実際の速度ＶＡＣＴと異なることになる。しかし、従来技術によれば、ＦＰＴＯ＝（ＶＰ）（ＺＰＴＯ）を維持し続けること、すなわち、因果的力と非因果的力を識別しないこと、また、ＦＰＴＯ＝（ＶＰ）（ＺＰＴＯ）を維持すること、すなわち、因果的力と非因果的力を識別せず、ＦＣ＝（ＶＰ）（β）およびＦＡ＝（ＶＰ）（ＺＰＴＯ−β）を維持することは、波誤差が存在する（すなわち、実際の波が予測される波と異なる）とき、ＷＥＣによってずっと小さな電力の生成をもたらす。たとえば、図１１および１２の波形Ｄで示す５６３ｋＷに対して、図１３の波形Ｄで示す４１２ｋＷの電力が生成される。力の因果的部分は、予測速度に比例することが留意されるべきである。
【００６１】
図１４は、波誤差の存在下で本発明を使用した制御を示す。そのため、この図は、ＷＥＣに衝突する実際の波が、波形Ａに示すように、予測される波と異なる条件を示す。しかし、本発明によれば、図１４に示す条件の場合、ＦＰＴＯは、因果的力（ＦＣ）および非因果的力（ＦＡ）を有するように処理され、ここで、ＦＣは、ＶＡＣＴにβを乗じた値に等しいようにセットされ、ＦＡ＝（ＶＰ）（ＺＰＴＯ−β）である。この処理配置構成によるＷＥＣによって生成される電力は、従来技術の未補正配置構成の場合に比べて著しく大きい。これは、図１３の波形Ｄに示す平均電力に対して図１４の波形Ｄに示す生成平均電力（５１１ｋＷ）を比較することによって立証されてもよい。力の因果的部分は、予測される速度ではなく、実際の／観測されたＷＥＣ速度に基づくことが留意されるべきである。
【００６２】
図１５は、特定の波条件下の特定のＷＥＣ構成についての、本発明の方法と「知られている」方法を比較する、変換電力についての結果の要約を示す。グラフのｘ軸は、波誤差の程度を表し、ｙ軸は、７秒の波の存在下で所与の波力エネルギー変換システムによって生成される電力を表す。ｙ軸は、決して予測波情報に依存しない制御則である「比例制御（proportional control）」によって生成される電力に対して正規化される。そのため、「１」未満の値は、制御則について失敗を示す。制御則は、「比例制御」の制御より大きな電力を生成しない場合、利益を全く提供しない。グラフの３つの曲線は、３つの異なる制御スキームのそれぞれにおいて生成される電力を示す。波誤差が全く存在しない場合（たとえば、グラフの左側で、ｘ軸がゼロである場合）、３つの制御スキームが全て一致していることに留意されたい。「Λ」曲線で示す制御スキームは、「知られている」制御スキームを表す。Γ１０００およびΓ２００曲線は、考えられる制御スキームのうちの２つについての電力の結果を示す。Γ１０００制御則が、ＷＥＣ位置の１０％の誤差が存在するまで、「比例制御」に対して利益を提供していることに留意されたい。Λ制御則は、５％未満の波誤差についてだけ、「比例制御」に対して利益を提供する。したがって、「Λ」制御則についての結果と比べて、Γ１０００の場合、１００％を超える改善が存在する。
【００６３】
【表１】
0005474056.tif 000033
【図１】
0005474056.tif 000034
【図２】
0005474056.tif 000035
【図３】
0005474056.tif 000036
【図４】
0005474056.tif 000037
【図５】
0005474056.tif 000038
【図６】
0005474056.tif 000039
【図７】
0005474056.tif 000040
【図８】
0005474056.tif 000041
【図９】
0005474056.tif 000042
【図１０】
0005474056.tif 000043
【図１１】
0005474056.tif 000044
【図１２】
0005474056.tif 000045
【図１３】
0005474056.tif 000046
【図１４】
0005474056.tif 000047
【図１５】
0005474056.tif 000048