"リアルタイムで、グローバルナビゲーション衛星システム-GNSSに基づいてセンチメートル単位の正確性で地震動又は構造物の振動を計測する方法であって、受信手段及び1Hz以上のサンプリング周波数におけるGNSS位相観測のための手段を備える受信機r、無線通信によってブロードキャストされる補正データを受信する手段、前記観測を処理する手段、以下の情報を貯蔵する記憶手段による方法であり、前記情報とは:-前記受信機の位置の複数の参照座標referencecoordinates、-リアルタイムでの無線通信によってブロードキャストされる補正データであって、少なくとも複数の位置推算表ephemerides、複数のクロック補正値、及び電離圏モデルを含む、ブロードキャストされる補正データ、-複数の処理結果であり、;当該方法は以下のステップを有する:-リアルタイムでの、少なくとも4つのGNSS源からの一対のGNSS位相観測値の受信と決定、及び、無線通信によりブロードキャストされた補正データの受信をすること、-1Hz以上のサンプリング周波数の連続的な時点t,t1において受信される一対の前記位相観測値のための時間内位相差phasedifferenceintimeの計算をすることであり、位相観測値の各対は前記複数のGNSS源の各々から来る、-バリオメータ位相式による前記時間内位相差の各々の表現をすることであって、該表現は、少なくとも4つのバリオメータ位相式の系を定義するためであり、各式は、位相観測値の各対のためのものであって、4つの未知数を含み、該未知数は、-前記連続的な時点t,t1の間に生じる三次元の変位の3つのデカルト成分Cartesiancomponents、-前記連続的な時点t,t1の間に生じる前記受信機のクロックエラーの変化、であり、-各バリオメータ位相式の重みづけ要素の計算をすること、-最小二乗法推定により、前記各4つの未知数に関し、少なくとも4つのバリオメータ位相式の系を解くこと。少なくとも前記少なくとも4つのGNSS源が、-複数の衛星の一つ以上の群constellationsの中の1つ以上の衛星、-1つ以上の擬似衛星、に属することを特徴とする、請求項1に記載の方法。少なくとも一つの信号源が擬似衛星である時に、前記擬似衛星の位置と、関連する一つのクロック補正データを得るステップを更に有することを特徴とする、請求項2に記載の方法。前記バリオメータ式が、の表現を有する一般的な位相観測式であり、ここで、sは、少なくとも4つのGNSS源の一つに関連し、rは受信機に関連し、λΔΦSrは1Hz以上のサンプリング周波数の連続する時点t,t1において受信される複数の位相観測値の間の位相差phasedifferenceであり、eSr●ΔξrcΔδtrは4つの未知数を備え、それらのうちの3つΔξrは連続する時点t,t1の間に生じる三次元の変位に関し、それらのうちの1つΔδtrは連続する時点t,t1の間に生じるクロックエラーの前記変化に関し、[ΔρSr]OR-cΔδtSは、無線通信を経て受信される前記補正データによって計算される既知の項であり、ΔεSrはノイズ要素である、請求項1乃至3のいずれか1項に記載の測定方法。前記バリオメータ位相式が、の種類であって、ここで、ΔTSr-ΔISrは、連続的な時点t,t1の間に生じ、前記補正データによって計算される大気屈折atmosphericrefractionの影響の変化を定義し、[ΔρSr]EtOIΔρSrは、連続的な時点t,t1の間に生じ、前記補正データによって計算される、固体地球潮汐solidEarthtideの効果の変化、海洋潮汐の効果の変化、及び相対論的効果を定義することを特徴とする、請求項1乃至4のいずれか1項に記載の測定方法。各バリオメータ位相式の前記重みづけが、以下の式から計算され、ここでZとは、前記受信機rの天頂角と前記少なくとも4つの衛星sのうちの1つの衛星との間の角度であり;バリオメータ位相式が、1つの擬似衛星から受信される単一の信号に基づいて計算されるならば、前記重みづけは1に等しいと仮定されることを特徴とする、請求項1乃至5のいずれか1項に記載の測定方法。連続的な時点t,t1の前記対を複数備える近接した時間間隔での受信機rの変位を計算するために、前記近接した時間間隔において、三次元の変位Δξrが足し合わされることを特徴とする、請求項1乃至6のいずれか1項に記載の測定方法。ゼロではない平均値を有する定誤差systematicerrorを排除するステップを更に有する、請求項1乃至7のいずれか1項に記載の測定方法。前記定誤差は、数分の近接した時間間隔において検知されることを特徴とする、請求項7に記載の方法。地震動又は構造物の振動をリアルタイムで計測する装置であって、請求項1に従う方法を実行する手段を備えることを特徴とする装置。前記GNSS受信機が二重周波数dualfrequencyのものであり、双方の周波数において受信される複数の信号に基づく前記方法を実行することが可能であることを特徴とする、請求項8に記載の装置。請求項1乃至9に記載の複数のステップを実行することが可能なプログラムコーディングのための手段を有するコンピュータプログラムであって、該プログラムはコンピュータ上で実行されることを特徴とするコンピュータプログラム。記録されたプログラムを備えるコンピュータ読み取り可能手段であって、前記コンピュータ読み取り可能手段は、前記プログラムがコンピュータ上で実行される際、請求項1乃至9のステップを実行することが可能なプログラムコーディングのための手段を有することを特徴とする、コンピュータ読み取り可能手段。本発明の目的は、とりわけ地震動や構造物の振動において、センチメートル単位での大きさを有する精度を伴う、運動の測定方法をもたらすことであり、該方法は、リアルタイムでも帰納的にも動作することが可能であり、少なくとも一つのGNSSの配置及び/又は擬似衛星からの少なくとも一つの周波数による、信号に基づく位相観測、及び、無線通信によるブロードキャスト補正データを用いた効果計算の補正を実行することが可能な単一の受信機によるものである。本発明が関連するのは、地震動又は構造物の振動を、グローバルナビゲーション衛星システム-GNSS及び/又は擬似衛星に基づいて測定するためのシステムであって、請求項1に従うものである。少なくとも一つのGNSSの配置及び/又は擬似衛星からの少なくとも一つの周波数に基づく位相観測を実行することが可能な受信機は、以降は、GNSS受信機と短く言及される。本発明の計算方法を、先行技術と差別化する基本的なコンセプトは、バリオメータ位相方程式が、連続的な観測値の複数の対のために計算され、当該観測値は、1Hzより大きいかそれに等しい周波数におけるものであり、全く同一の源によって受け取られるものであるが、この動作は、4つの未知数における4つの方程式の系を得るために、同時に少なくとも4つの源のために、繰り返される。4つの未知数の内1つの未知数は、クローズドシステムを得るために無線通信による補正データ放送によって決定される。実際、他の3つの未知数は、当該系を解き、前記連続的な観測値の対によって定義される時間間隔の間に生じる変位を決定することにより、決定される。有利なことに、このことが可能とするのは、無線通信による、リアルタイムで利用できる不正確な補正データブロードキャストを用い、単一の受信機を用いてセンチメートル単位での変位のリアルタイムでの測定値を得ることである。本発明の更なる目的は、グローバルナビゲーション衛星システム-GNSS及び/又は擬似衛星に基づき、構造物の地震動又は振動を測定する装置をもたらすことであり、それによって前記問題が解決され得る。本発明はまた、請求項9に従うグローバルナビゲーション衛星システム-GNSSに基づく構造物の地震動又は振動を測定するシステムにも関連する。有利なことに、本発明は、リアルタイムでもオフラインモードにおいても、ただ一つのGNSS受信機を用いて、センチメートル単位の正確性を伴って、構造物の地震による地震動変位及び振動を決定することを可能とする。とりわけリアルタイムにおける、これらの量の決定は、津波などの壊滅的な出来事のトリガーとなることができ、従って、全住民に警告し警報するための基礎的な情報を表す、地震による変位を適時に検知することを可能とする。更に本発明が可能とするのは、大きな地震の震央の近くに位置する地震計に共通して存在する飽和の問題を避けながら、地震モーメントとマグニチュードを決定することである。有利なことに、本願発明は、先行技術と比較して、リアルタイムでもオフラインでも、ただ一つの受信機と補正データブロードキャストを用いて、地震による変位を決定することを可能とする。更に、本願発明は、リアルタイムでもオフラインでも、橋、タワー、高層ビル等の大構造物の振動のモードに関連する変位の度合いを、センチメートルの大きさの単位の正確さで決定することにも応用される。複数の請求項は、本発明の好適実施例について述べており、本記述に不可欠な部分を構成する。本発明の更なる特徴と長所は、本発明を実施する方法の、好適だが排他的ではない実施例の詳細な記述からより明確になるが、以下に付随する図面を用いた非限定例として説明される。図面中の同一の参照番号及び文字は、同一の要素または構成要素を特定する。本発明の好適実施例の詳細な説明ただ一つのGNSS受信機の使用を構想する方法は、そのファームウェアに有利に実装することができ、それによって、リアルタイムに、しかし帰納的にもオフラインで、センチメートルの大きさの単位での正確性を伴うが、周波数の数と追跡される群constellationsの数に依存した変位の測定をもたらし、更なる処理の必要性は伴わない。本発明の方法は、GNSS受信機によって好ましく実装され、該GNSS受信機は、1Hz又はそれより高いサンプリング周波数を用いた受信と位相観測のための手段、前記位相観測の処理手段、及び記憶手段を備え、該記憶手段は、-受信機の位置の参照座標、-無線通信によってブロードキャストされる補正データであって、位置推算表、クロック補正、及び、好ましくは同一のGNSS衛星からリアルタイムで、同一の受信手段によって、位相観測値の受信と同時に受信されるブロードキャスト電離圏モデルを含む、ブロードキャストされる補正データ、-処理結果、を記憶する。当該方法は、以下の基本的なステップを有する:-視野にある少なくとも4つの衛星からのGNSS位相観測値の受信、-以下のバリオメータ位相式と呼ばれる種類の方程式による時間内位相差aphasedifferenceintimeの表現であって、該方程式は、連続的な時点t,t1において受信される、全く同一の源から各源、すなわち前記4つの衛星、及び/又は擬似衛星に対する、少なくとも一対の信号の方程式であり、ここで、[λΔΦSr]は、時間内位相差、又は、汎用的な周波数generalfrequency又は複数の周波数の組に基づく二つの連続的な期間における測定において得られたGNSS位相観測値の差分であり、eSr●ΔξrcΔδtrは、4つの未知のパラメータ、すなわち、3方向における変位Δξrと、クロックエラーの変化cΔδtrを含む項であり、[ΔρSr]OR-cΔδtSは、追跡されるGNSS衛星の位置推算表及び同期誤差のために利用できる情報この情報は、リアルタイムでナビゲーションメッセージにおいて利用可能であるに基づいて計算される周知の項であり、ΔεSrは、電気的ノイズである;-各バリオメータ式に対する数式w=cos2Zを用いた重みづけの計算であって、ここで、Zは衛星の天頂角から受信機に対するものである;-Δξrについての、4つの未知数における少なくとも4つの方程式の系を解くことであって、その3つの未知数は、受信機の3方向の変位を表し、4つめの未知数は、クロックエラーの変化に関連する;このようにして、前記連続的な時点の対によって定義され、1秒以下の時間間隔において生じるいかなる変位も決定され、連続的な時点t、t1の他の対のために計算される3方向の変位が合計されるが、その目的は、前記複数の時点の複数の対を備える割り当てられた時間間隔における、受信機の変位の推移courseを再構成することである。補正データは、多くの別の源、例えば前記GNSS衛星、及び/又は擬似衛星から伝送されることが可能であるが、それはリアルタイムで受信することを目的としており、同GNSSのアンテナと同時に位相観測値の決定のために必要な信号が伴う。その結果の正確性は、バリオメータ式に更なる追加の項を挿入し、生じる影響/擾乱のモデリングによって高めることができる。:ΔTSr-ΔISrは、大気伝搬の影響を考慮することができ、該影響は複数の周波数の適切な組み合わせによってモデル化され又は消去される。[ΔρSr]EtOIΔρSrは、地球と海洋の潮汐の効果と相対論的効果を考慮することができる。ΔmSrは、伝送される信号のマルチパスである。それらは、以下の数式に至る。この数式2は、位相観測式の一般的な数式から直接得ることができる。更なるステップは、何らかの定誤差の影響を除去するために実行されるが、該定誤差は、前記の連続的な複数の時点の対の間に生じる変位Δξrの計算において達成することの出来る正確性に比較すれば無視できるように見えるが、三次元の変位の合計においては、大きな影響を伴って蓄積される。我々は、位相観測式の周知の一般的な数式、例えば、Hoffman-Wellenhof他2008の数式について有利に考察してもよいが、当該数式は長さの単位で書かれ:ここで、ΦSrは、衛星sに対する受信機rの位相観測値である;λは位相の波長である;ρSrは衛星sと受信機rとの間の幾何学的距離である;cは光の速さである;δtrとδtSは、受信機rと衛星sのクロックエラーである;TSrとISrは、衛星sから受信機rへのパスに沿った対流圏及び電離圏の遅延である;NSrは、初期位相の曖昧度である;ρSrはその他の効果相対論的効果、位相中心の変化、フェーズワインドアップphasewindupの合計である;mSrとεSrは、マルチパス効果及びその誤差である。本発明によれば、単一の差分が、数式aによって述べられる位相観測の、二つの連続する時点t,t1の間の時間において計算される。高周波数での位相観測値が用いられる、すなわち、1Hz以上の周波数における位相観測値が用いられると仮定するならば、前記位相差分式の第二の表現が得られる:受信機の位置は、ECEFEarthCentredEarthFixedレファレンスシステムに固定されることが好ましい;そして、当該差分式bの右辺の第1項ΔρSrt,t1は、衛星と受信機の間の距離の変化にのみ依存し、地球の潮汐及び海洋荷重のとても小さな影響を別にすれば、衛星の軌道運動と地球の自転の双方によって決定される。前記第1項は、以下のように同等とされる。受信機が、二つの連続する時点t,t1の間でECEFレファレンスシステムに比較した変位Δξrt,t1を受けると、前記第1項ΔρSrt,t1は、衛星sと受信機rの間の見通し線に沿って投影された変位Δξrt,t1の影響も含み、当該変位は、前記二つの連続する時点の間で、一定値を維持すると仮定される。従って、前記第1項は、以下のように同等とされる。ここで、eSrは、衛星sと受信機rの間のベルソルversorであって、記号●が示すのは、ベルソルeSrとΔξrt,t1の間の内積である。右辺の第1項の第2の表現を置換すると、差分式bにおいて、差し当たり時間間隔に対する言及を省略して、バリオメータ式が得られる。:今や、項[ΔρSr]ORとΔδtSには特別な注意が払われるであろう。:現在の状態において、本発明によれば、一時的には1秒も離れていない連続的な時点の複数の対が考慮されるので、項[ΔρSr]ORとΔδtSとを計算するために、GNSS受信機からリアルタイムで得られるブロードキャスト位置推算表broadcastephemeridesとクロックエラーを用いることが可能であり、その際は、1ミリメートル未満の誤差を得る。とりわけ、これらのデータは、各時点における衛星の位置を計算するために必要なケプラー軌道パラメータと、衛星のクロックの同期誤差のドリフトの放物線モデルの係数を表すが、それは、これらのデータが正確なタイプの前記積に対して最小のドリフトを示して以降のものである。従って、このバリオメータ式は、以下の好ましい形態を前提とするが、それは2に一致する。ここで、λΔΦSrは観測値の差分であり、eSr●ΔξrcΔδtrは4つの未知のパラメータ、すなわち3つの未知数を定義する三次元ECEFシステムにおける変位Δξrと、クロックエラーの変化を含む項であり、[ΔρSr]OR-cΔδtSΔTSr-ΔISr[ΔρSr]EtOIΔρSrは、伝送された位置推算表と適切なモデルに基づいて計算される既知の項であり、ΔmSrはマルチパスであり、ΔεSrはノイズである。好ましい形態eにおけるバリオメータ式は、連続的な時点の各対のために、受信機の変位を決定するための最小二乗法による推定において用いられる機能的モデルを表している。低い仰角の衛星から来る観測値はノイズが多いことが知られており、このような訳で、推定方法のための確率論的モデルが想定するのは、観測値wλΔΦSrを、衛星の天頂角から受信機に対する角Zのコサイン二乗に等しくする重みづけの利用である。ここで、Zとは、衛星の天頂角から受信機に対するものである。最小二乗法による推定は、二つの一般的な連続する時点のために書かれ得る一連のバリオメータ式の系に基づく。;それらの数は、前記二つの連続する時点の間の、複数の衛星、及び又は擬似衛星の利用可能性に依存する。位相観測値の数を増やすことを目的とするなら、受信機が、複数の周波数に基づく複数の群constellationsからの信号を受信し、それによって、同時に多くの位相観測値を受信することが可能となることが望ましい。更に、複数の周波数に基づく複数の位相観測値の利用が可能とするのは、異なる複数の周波数に関連する複数のバリオメータ式の適切な線形結合によって電離圏効果ΔISrを排除することである。受信機の変位を検知するために、特定の時間間隔の間の前記好ましいバリオメータ式eによって推定される複数の変位が合計される。周知の項[ΔρSr]OR-cΔδtSΔTSr-ΔISr[ΔρSr]EtOIΔρSrの計算に伴い得る誤差は、連続する時点の単一の対にとっては小さいが、複数の連続する時間の対のために時間と共に合計されるなら、大きなものとなりがちである。;そして、累積する変位ΣΔξrの時系列が示すのは、低周波数の要素トレンドである。好ましくは、本発明の方法は、変位Δξrの時系列から、前記トレンドを直接排除するステップを備え、当該ステップは、運動学的に記載されるイベント地震性のイベント、構造の振動等の時間に先立つ適切な間隔を考慮し、低水準の多項式補間、例えば二次への補間を用いてなされる。;好ましくは、この多項式補間は、ロバスト推定、例えば、いわゆるルスーRousseewのリーストトリムドスクエアLeastTrimmedSquares種の推定を用いて実施されるべきである。リアルタイムでの、しかしオフラインでも用いられる、変位の決定方法の好適実施例は、少なくとも一つのGNSSの配置constellationからの少なくとも一つの周波数に基づく位相観測値を得ることの出来るGNSS受信機に基づく。前記GNSS受信機は、その位置が国際標準のWGS84システムにおいて数メートルの正確性を伴って知られなければならないがいかなるGNSS受信機によっても簡易に得ることができる、少なくとも以下のデータを得る。:-GNSS配置からの1秒以下の抽出間隔による一つの周波数に基づく位相観測値、-少なくとも、観測されるGNSS配置のブロードキャスト位置推算表とクロック補正値を含むナビゲーションメッセージ。従って、本発明による好適な計算モデルは、確率論的モデルを定義する、前記バリオメータ式と前記観測値の重みづけ要素に基づく:複数の変位が、得られる複数の係数の変化の最小二乗法による推定によって決定されるが、それは、前記確率論的モデルによる前記バリオメータ式の利用によるものである。明らかに当該方法のパフォーマンスは、周波数の数と、GNSS受信機が追跡することのできる群constellationsの数に依存する。というのも、これらの特徴量が、書かれることの出来るバリオメータ式の数と可能な組み合わせを決定するからである。これらの3つの入力データ、すなわち、位相観測値、ナビゲーションメッセージ、WGS84システムにおける位置は、動作している間利用可能であるが、本発明に従う方法は、それらのデータから始まって、以下のステップを有し、好ましくはGNSS受信機のファームウェアに実装される。:-ステップ1、二つの部分的な間隔から成る、全ての変位の計算のグローバルな時間間隔の定義:イベントが現れるのに先立つ間隔Δtaであって、運動学的に記載され地震性のイベント、構造の振動等、少なくとも1分の継続時間を有する、前記イベントが現れる間隔Δteであって;例えば、地震にとっては、地震性の現象が現れる間隔であり、一般的には1分から3分である;-ステップ2、前記時間間隔の間、連続的な時点の一つの組における少なくとも4つの衛星に関連する位相観測値の間の差分の計算であって、複数の時点の各対は、1秒以下の間隔を定義する;-ステップ3、連続する時点t,t1の各対に関連するバリオメータ式の右辺の[ΔρSr]OR-cΔδtSΔTSr-ΔISr[ΔρSr]EtOIΔρSrの項の計算であって、各衛星における、当該衛星の各位置推算表に基づく二つの同一の連続的な時点における各周波数のための、各配置constellationに関連する複数の位相観測値の変化量と、受信機の位置と、衛星のクロック補正、及び/又は、例えば当該計算に影響を与える適切な要素による計算;-ステップ4、バリオメータ式の右辺の変位の未知数の係数[eSr]を計算するための、連続する時点t,t1の各対のための当該手段の受信機-衛星の方向のベルソルの計算;-ステップ5、各衛星の連続する時点t,t1の各対のための各バリオメータ式λΔΦStt,t1のための重みづけ要素の計算であって、の式を用い、前記二つの連続する時点における当該手段の受信機-衛星の方向のベルソルに基づくステップ3;-ステップ6、三つの変位の未知数Δξr、及びクロックエラーΔδtrとクロックの正確性の変化に関連する未知数の、連続する時点t,t1の各対のための最小二乗法による推定;-ステップ7、ステップ1で定義される計算間隔が完了していないことの認証;仮にまだ完了していないのであれば、それはステップ2から、すなわち、連続する時点の一つの対における少なくとも4つの衛星に関連する位相観測値の計算から再開する、でなければ、-ステップ8、ポイント1で定義される全計算間隔における最小二乗法によって推定される複数の変位Δξrの合計;-ステップ9、仮にステップ1で定義される計算時間間隔Δtaに関連する累積変位ΣΔξrの時系列が、低周波数要素トレンドを示すのであれば、-ステップ10、低水準の多項式補間を用いて、ステップ1における全グローバルな計算間隔における前記トレンドの排除、でなければ、仮にステップ9において低周波の構成要素が特定されないのであれば、-ステップ11、総変位の計算であり、従って、事前に決められた計算間隔の終わりに利用できる計算。ステップ2の、連続する時点の対における少なくとも4つの源衛星、及び/又は擬似衛星に関する位相観測値の間の差分の前記計算は、ステップ1で定義される前記総計算間隔の連続する時点の各対のために、各配置constellationのために、各衛星のために、各観察される周波数のために行われるが、それは、バリオメータ式の左辺[λΔΦSrt,t1]を計算することを目的とする。;しかし、仮に、前記連続する時点の間の間隔において、一つの衛星と受信機との間で、いわゆるサイクルスリップcycleslip、すなわち、計測の中断を引き起こすイベントが起こるならば、その衛星、その周波数、及び連続する時点のその対に関連するバリオメータ式は考慮されない、すなわち、破棄される。仮に、ある時点において、少なくとも4つの衛星が同時に利用できないのであれば、当該方法は、ステップ2から再スタートする。様々な好適実施例において説明される要素及び特徴は組み合わせることが可能であるが、それは、本特許出願の保護範囲内にある。上記に示すバリオメータによるアプローチに従った、リアルタイムでの変位の決定方法は、以下を用いても適用することが可能である。:現在運転中の全てのGNSS群から、及び、認識段階にあるGNSS群から、及び/又は擬似衛星からの複数の観測値。現在運転中の全てのGNSS群によって、及び、認識段階にあるGNSS群によって利用可能とされる一つ以上の周波数及び、適用可能なそれらの組み合わせ。様々な種類の複数のGNSS受信機。ブロードキャストされる補正データの中で直接利用可能な、ブロードキャスト電離圏モデルに関連する情報、又は、ネットワークにおいて利用可能な電離圏モデルの情報であって、当該アルゴリズムの正確性を高めるための情報。複数の擬似衛星によって発生する複数の信号であって、現在運転中の一つ以上のGNSS群から来る信号、及び、認識段階にある一つ以上のGNSS群から来る信号に加えて、又はそれらに置き換えて用いられる複数の信号。更に、当該バリオメータアプローチによる複数の変位の決定方法は、ブロードキャスト補正データ及び正確なデータ衛星の正確な位置推算表及びクロックエラー等と共に、オフラインモードすなわちリアルタイムではなくで適用されることも可能である。しかし、リアルタイムでセンチメートル単位の測定を実行するためには、ブロードキャスト補正データを実行することが望ましい。ここに述べる方法に従う複数の変位の決定方法の機能性及び有効性は、第一に、ブロードキャスト補正データと、重要な地震性のイベントによって影響される、常置局から得られるGPS位相観測値を用いて、オフラインモードで得られる実験結果によって、正当性が立証された。同結果は、同一の又はそれ以上の測定の正確性を得ることの出来る先行技術の方法によって得られる実験結果とも比較され、センチメートル単位の大きさで極めてよく一致した。最終的には、当該方法は、リアルタイムでも成功裏に正当性が立証された。"