【手続補正書】
【提出日】平成15年7月25日(2003.7.25)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の名称】楕円波の識別方法とその装置および地震探査システム
【特許請求の範囲】
【請求項１】少なくとも２つの直交方向の波の成分について、少なくとも１つの受信機（Ｒ）によって検出される、物質媒体中の少なくとも１つの伝搬平面で伝搬する楕円波を識別する楕円波識別方法において、
受信された波に応じて受信機によって発生された信号波の変換を決定する段階と、
前記変換された信号波の各振幅の比に基づいた基準によってこれらの変換された信号波の振幅を選択する段階とを有することを特徴とする楕円波識別方法。
【請求項２】ろ波すべき楕円波の楕円率を有する選択ウインドウ（ｆ）を選択する段階と、
一連のスケール係数ｆと選択された分析波の時間（Ｔ）に対応した一連の連続したサンプル（ＴＯn（Ｚ）、ＴＯn（Ｘ））を得るために、少なくとも２つの直行する方向（Ｚ、Ｘ）で受信した信号波を変換する段階と、
次に、比
ｒ1＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）、および、
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ）、
ここでｎ−１，ｎ，ｎ＋１は連続整数、
を決定する段階と、
前記比の値ｒ1，ｒ2を選択ウインドウに含むサンプルを有効にする段階と、
逆の信号波変換を行うことによって楕円波フィールドを再構成する段階とを有することを特徴とする請求項１に記載の楕円波識別方法。
【請求項３】前記再構成の楕円波フィールドが受信された信号から減算される請求項１または２に記載の楕円波識別方法。
【請求項４】前記媒体が地質構造である請求項１乃至３のいずれか１項に記載の楕円波識別方法。
【請求項５】各受信機は前記構造の媒体に結合されたジオホンまたは加速度計からなるいくつかの方向性検知器を有する請求項４に記載の楕円波識別方法。
【請求項６】楕円表面波をろ波するために地面の表面近傍に少なくとも１つの受信機を設定する段階を有する請求項１乃至５のいずれか１項に記載の方法。
【請求項７】井戸に沿って伝搬する楕円波をろ波するために、この地質構造に堀込んだ少なくとも１つの井戸内に、少なくとも１つの受信機を設定する段階を有する請求項１乃至６のいずれか１項に記載の方法。
【請求項８】地下層に対して音響または地震波を出力する段階と、
少なくとも２つの異なる方向で受信した波の成分（Ｘ，Ｙ，Ｚ）を検出するのに適した１つ、または、いくつかの受信機によって地質構造の断層によって反射された波を受信して信号を記録する段階と、
受信された波に応じて受信機に発生した信号波の変換の決定と、その各振幅の比に基づいた基準によってこれらの変換された信号波の振幅の選択とを含む、楕円波を無くすように記録された信号を処理する段階とを有する地層構造の地震探査方法。
【請求項９】前記記録された信号を処理する段階が、ろ波すべき楕円波の楕円率を含む選択ウインドウ（ｆ）を選択する段階と、
一連のスケール係数（ｆ）と、選択された分析波の時間（Ｔ）に対応した一連の連続したサンプル（ＴＯn（Ｚ）、ＴＯn（Ｘ））を得るために、少なくとも２つの直交する方向（Ｚ、Ｘ）で受けた信号の信号波変換を行う段階と、
下記の比を決定する段階と、
ｒ1＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）、および、
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ）、
ここで（ｎ−１，ｎ，ｎ＋１）は連続整数、
前記比の値ｒ1，ｒ2を選択ウインドウに含むサンプルを有効にする段階と、
逆の信号波変換を行うことによって楕円波フィールドを再構成する段階とを有することを特徴とする請求項８に記載の地震探査方法。
【請求項１０】前記再構成された楕円波フィールドが前記受信された信号から減算されることを特徴とする請求項８または９に記載の地震探査方法。
【請求項１１】受信された波に応じて受信機（Ｒ）によって生じた信号波の変換の決定と、その各振幅の比に基づいた基準によってこれらの変換された信号波の振幅の選択とに適する処理装置（１）を有する、少なくとも２つの直交方向で波の成分を検出するのに適し、前記媒体に結合された少なくとも１つの受信機（Ｒ）によって受信された物質媒体中を伝搬する音響または地震波の間で少なくとも１つの伝搬平面で伝搬する楕円波を識別する楕円波識別装置。
【請求項１２】前記処理装置は、少なくとも２つの直交する方向で受けられる信号の信号波変換を行う手段と、
一連のスケール係数（ｆ）と、選択された分析波の時間（Ｔ）に対応した一連の連続したサンプル（ＴＯn（Ｚ），ＴＯn（Ｘ））を得る手段と、
下記の比を決定する手段と
ｒ1＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）、および、
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ）
ここにおいてｎ−１，ｎ，ｎ＋１は連続整数、
前記比の値ｒ1，ｒ2をろ波すべき楕円波の楕円率を有する選択ウインドウに含むサンプルを有効にする手段と、
逆の信号波変換を行うことによって楕円波フィールドを再構成する手段とを有する請求項１１に記載の楕円波識別装置。
【請求項１３】前記処理装置（１）が、受信された信号から再構成された前記楕円波フィールドを減算する手段を有する請求項１１または１２に記載の楕円波識別装置。
【請求項１４】前記処理装置（１）はプログラムされたデジタルコンピュータを有することを特徴とする請求項１１または１２に記載の楕円波識別装置。
【請求項１５】地質構造中に波を発生する地震源と、
少なくとも２つの直交方向に受信された波の成分を検出するために適した前記地質構造と結合した少なくとも１つの受信機（Ｒ）と、
地震波の送出−受信のサイクルの実行を制御する制御装置と、
各受信機によって受信された信号の記録とを有し、請求項１１または１２に記載の楕円波識別装置を有する地質構造の地震探査システム。
【請求項１６】少なくとも一つの受信機に近接して設けられた少なくとも１つの地域的な局所集録ユニット（ＲＴＵ）を有し、処理手段を有する楕円波識別装置が前記局所集録ユニットに含まれている請求項１５に記載の地震探査システム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、多軸の受信機によって受信された電波の複数の軸に沿って測定された成分の組み合わされた処理によって、物質媒体で伝搬する電波の中から楕円波を識別することができるろ波方法に関する。
【０００２】
【従来の技術】
本発明による方法は、種々の分野で用途を見いだすことができる。本発明は、地質構造に適用されたディスターバンス（変動）に応じて地質的な特徴によって放射される波が３軸ジオホンすなわち"トリホン"のような地面に設置された受信機によって受信される地球物理学的な探査、および、地震探査に特に適用することができる。
【０００３】
特に、地震探査において、表面波のような波は深く浸透しないことがよく分かっている。従って、それらは地球物理学的な探査にとって有効な情報は提供せず、ノイズと考慮される。それらのエネルギーは、有効な信号のエネルギーと比較して往々にして高い。従って表面波はろ波しなければならない。
【０００４】
２つのタイプの表面波は弾性を有する一様な半無限媒体の表面に沿って伝搬する。２つのタイプの波とは、粒子の動きが伝搬平面に残るレイリ（Ｒａｙｌｅｉｇｈ）波と粒子の動きが伝搬平面に直角なラブ（Ｌｏｖｅ）波とである。レイリ波としての粒子の動きは三角関数（レトログレード、反時計）方向をカバーする楕円に追随する。楕円の軸線は垂直方向であり、主軸と第２の軸との間の比は媒体の弾性パラメータに依存するが、１．４乃至１．９の範囲内になる。
【０００５】
地質学的な構成は、その表面の近くでウエザー層と称される低速度層を有し、この低速度層でフロア上の反射とこの低速度層の上部との間でいくつかの周波数について正干渉によって形成された（疑似レイリ波と称される）案内波が形成される。またこれらの案内波は、源と受信機を含む垂直方向の平面において伝搬する楕円波である。これらの波のいくつかについては粒子の動きは、レイリ波の特性としてレトログレードであり、他の波に関しては直接またはプログレード、すなわち時計回りである。
【０００６】
基本軸の比はウエザー層の厚さ、および、この層と下層との間の速度コントラストによって変化する。ウエザー層のレイリ波と種々の案内波モードによって形成された全体は非常に分散的であり、明らかな速度は、周波数で変化し、広い範囲をカバーすることができる。その結果、表面波はシスモグラムの大きな部分をカバーし、それによって大きなエネルギーによる有効な情報を隠す。これらの波のすべては伝搬平面に含まれる。それらは楕円形であり、一般に直線的に極性化される直接、および、反射波から識別される。その基本軸は垂直方向、および、水平方向に接近している。
【０００７】
本体波と同様に、表面波は過渡的である。疑似レイリモードの各々によるディスターバンスは、その特定の速度で受信点を通過する。作用、および、反作用の原理によって、ディスターバンスの時間積分はゼロであり、要するに、１つの点で記録されたディスターバンスはゼロ領域の小さい波である。
【０００８】
これらを合計するためにレイリタイプの表面波は、認識基準を利用できるときに容易に配置することができる基本部分と、本体波とは異なり、認識基準を与える楕円部分と、受信点に対応するトレースでそれらの持続時間が１つまたはいくつかの振動に制限される過渡部分と、時間ウインドウで制限された周波数定数を有する分散部分とを有する。
さらに詳細には表面波については次の書類において示されている。
【０００９】
１９５３年トルストイＩ．、ウスディンＥ．(Tolstoy I.,Usdin E.)による"層状に形成された弾性液体媒体の分散性特性：レイ理論"地球物理学第１８巻８４４乃至８７０ページまたは１９９６年ムーニーＨ．、ボルトＢ．(Mooney H., Bolt B.)による"１つの表面層において第１の３つのレイリモードの分散特性"に記載されている。
【００１０】
従来のろ波方法は、直列の相互接続されているジオホンの"ストリング"を地震受信機として使用し、その各信号を重ねることから成る。さらに有効には、表面波の最も長い波長に等しい距離または少なくとも２つのトレースの間の距離に等しい距離をカバーしなければならない。これは、数十メートルのケーブル、大きなスペースと重量、および、大きな位置決め、および、メインテナンス作業を必要とする。このろ波タイプは、表面波と反射信号との間の明らかな水平方向の速度差に基づいて位置決めの方向にのみ有効である。３Ｄにおいて使用可能な等方性の地面のフィルタの製造は実際に実施するには重すぎる。
【００１１】
地面フィルタを使用した後にシスモグラムに残る表面波をろ波する他の公知の方法は、ある時間制限の後、すなわち、表面波の到着から、ジオホンによってピックアップされたものを無音とするか、ローカットフィルタを使用することによってある値以下のすべての周波数を削除することによることから成る。これらのフィルタは分離されたドメインへの記録の分解、および、ノイズのあるドメインの削除に基づいている。
【００１２】
複数の軸線の地震受信機を使用し、少なくとも２つの直交軸線に沿って信号を処理して表面波の識別を行うことができる公知の地震探査方法がある。
【００１３】
例としての米国特許第ＵＳ−Ａ−４，４５８，３４１号は、レイリ波を検出するに適した処理装置と関連する２軸または３軸のピックアップによって層の構成からくる地震信号を受信する地震探査方法を示している。信号ピックアップの半径方向Ｘ、および、垂直方向Ｖ成分は、この処理装置のろ波、エンベロープ検出、および、比較回路に提供することができ、この処理装置は、その特別の特徴からレイリ波の受信ウインドウを決定して受信信号からの減算を制御することができる。
【００１４】
米国特許第４，９３５，９０５号は、層の構成からくる波を受けるために２軸または３軸ジオホンを使用する地震探査方法を示す。ピックアップされる種々の成分は極座標で表される。レイリ波の成分は、平均値方法を適用することによってこれらの成分から得られ減算することができる。
【００１５】
複数軸のジオホンによって表面波をろ波する方法は、１９９５年グラスゴーの５７回ＥＡＥＧ会議の"拡張アブストラクト"で示されている。
【００１６】
また特許第ＦＲ−Ａ−２，６４８，５６７号は、電波分析技術を必要とする、例えば脳波の分析に適用される信号デジタル処理方法を示している。
【００１７】
混合時間周波数ろ波：
周波数ろ波は、表面波（１５Ｈｚまでの最も低い周波数あるいはそれ以下の周波数）の周波数成分のフーリエドメインにおける削除、および、セイブドドメインからの信号の再構成が続くトレースの高周波分解である。勿論ノイズのこれらの時間または周波数ドメイン特徴の削除は、これらのドメインに属する有効な信号部分の削除を招く。削除すべきノイズから認識されるドメインと信号を引き出すことが望ましい場合残りの部分との間で重複を制限するために、２つの二次元の時間周波数分解を実行することは公知である。
【００１８】
本発明による方法は、"波分析"と称される"時間−周波数"信号分析を使用し、その場合、信号は、幾何学数列によってサイズが変化する波の配列によって信号が分析され、各波は、その長さの部分と等しいステップでトレースに沿って解釈される。この技術は、例えば、１９８４／１９８５年グーピロードＰ．とグロスマンＡ．とモーレットＪ．（Goupillaud P., Grossmann A., Morlet J.）による地質探査第２３巻の８５乃至１０２ページの"地震信号分析におけるサイクルオクターブ、および、関連変換"に示されている。
【００１９】
信号波変換の分析信号は、ゼロ平均値の信号である。例えば、リッカーインパルス
【００２０】
【数１】
1997068578.tif 000001
を使用することができる。
【００２１】
信号波変換において、信号は一定の形状を有する。信号波変換による時間・周波数分解はディスクリート（サンプル）オペレーションである。与えられた寸法に関して、信号はせいぜい１つの期間（サイクル）に等しい一定のステップでシフトされる。信号のサイズは変化する。２（多くても、すなわち、１オクターブ）までの一定の幾何学的な比例によってスケール係数Ｅを変化することにより拡張または圧縮される。時間Ｔでのサンプル、および、信号波変換のスケール係数Ｅは、時間Ｔの中心のスケール係数Ｅの波によるトレースの結果である。
【００２２】
もしｈ（ｔ）が波のベースで、ａが信号波のスケール係数で、ｂによって変位されるとき、
【００２３】
【数２】
1997068578.tif 000002
ここで、ａは、スケール係数またはパラメータであり、ｂは時間トランスレーションパラメータである。信号の波変換は次のように表される。
【００２４】
【数３】
1997068578.tif 000003
逆変換は、係数の範囲内について、対応する波による信号波変換の各項の積の合計である。
【００２５】
【数４】
1997068578.tif 000004
【００２６】
【課題を解決するための手段】
本発明による方法は、（例えば地質構造のような）物質媒体で伝搬する他の波から少なくとも１つの伝搬平面で伝搬する楕円波を識別することができる。それは、少なくとも２つの直交方向の波の成分に関する少なくとも１つの受信機（前記構造と結合されたジオホンまたは加速度計のような１つまたはいくつかの検出器を有する）による検出を含む。
【００２７】
さらに、本方法は、受信された波にじて受信機によって生じた信号波の変換の決定、および、その各振幅の比に基づいた基準によってこれらの変換された信号波の振幅の選択とを有することを特徴とする。
【００２８】
１つの実施例によれば、ろ波すべき楕円波の楕円率を有する選択ウインドウ（ｆ）を選択する段階と、一連のスケール係数（ｆ）と、選択された分析波の時間（Ｔ）とに対応した一連の連続したサンプル（ＴＯn（Ｚ），ＴＯn（Ｘ））を得るために、少なくとも２つの直行する方向で受信する信号の信号波変換を行う段階と、
下記の比を決定する段階と、
ｒ1＝ＴＯn，（Ｚ）／ＴＯ(n-1)（Ｘ），および，
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ），
ここで（ｎ−１，ｎ，ｎ＋１）は連続整数、
前記比の値ｒ1，ｒ2を選択ウインドウに含むサンプルを有効にする段階と、逆の信号波変換を行うことによって楕円波フィールドを再構成する段階とを有し、この再構成の楕円波フィールドは受信された信号から減算される。
【００２９】
この識別方法は、表面波またはチューブ波と称される波をろ波するように地面の近傍または地中に１つまたはいくつかの、多軸受信機を設定することによって地質構造中を伝搬する楕円波をろ波するために使用される。
【００３０】
この識別方法は、音響または地震波が地質構造に出力され、断層構造によって反射された波が少なくとも２つの異なる方向で、その波の成分を検出するために適した１つまたはいくつかの受信機によって受信され、ピックアップされた信号が記録され、その記録された信号は、受信された波に応じて各受信機に発生された信号の信号波変換を決定し、かつ、各振幅の比に基づいた基準によってこの信号波変換から生じる信号の振幅を選択することによって、楕円波を消去するために処理される。
【００３１】
地震探査方法の１つの実施例によれば、ろ波すべき楕円波の楕円率を有する選択ウインドウ（ｆ）を選択する段階と、少なくとも２つの直交する方向で受信される信号の信号波変換を行い一連のスケール係数（ｆ）と、選択された分析波の時間とに対応した一連の連続したサンプルＴＯ（Ｘ）およびＴＯ（Ｚ）を得る段階と、比
ｒ1＝ＴＯｎ（Ｚ）／ＴＯ(n-1)（Ｘ）、および、
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ），
ここにおいて（ｎ−１，ｎ，ｎ＋１）は連続整数、
を決定する段階と、前記比の値ｒ1，ｒ2を選択ウインドウに含むサンプルを有効にする段階と、逆の波変換を行うことによって楕円波フィールドを再構成する段階によって記録信号を処理し、再構成された楕円波フィールダを受信された信号から減算する。
【００３２】
本発明による装置は、（例えば、地質学的な用途における２軸または３軸のような）少なくとも２つの直交する方向において波の成分を検出するのに適した、媒体と結合した少なくとも１つの受信機によって、（例えば、地質構造のような）物質媒体中を伝搬する音響または地震波から楕円波を識別することができる。
【００３３】
それは、受信された波に応じて受信機に生じた信号の信号波変換を決定し、および、その信号波の各振幅の比に基づいた基準によってこれらの変換された信号波の振幅を選択するのに適した処理装置を有することを特徴とする。
【００３４】
処理装置は、少なくとも２つの直交する方向で受けられる信号の信号波変形を行い、一連のスケール係数（ｆ）と、選択された分析波の時間（Ｔ）とに対応した一連の連続したサンプル（ＴＯn（Ｚ）、ＴＯn（Ｘ））を得る装置と、比率
ｒ1＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）、および、
ｒ2＝−ＴＯn／ＴＯ(n+1)（Ｘ）
ここにおいて（ｎ−１，ｎ，ｎ＋１）は連続整数、
を決定する装置と、前記比の値ｒ1，ｒ2をろ波すべき楕円波の楕円率を有する選択ウインドウに含むサンプルを有効にする装置と、逆の信号波変換を行うことによって楕円波フィールドを再構成する装置とを有し、再構成された楕円波フイールドは、その後、受信された信号から減算によって削除される。
【００３５】
また、本発明は地質構造において波を発生する地震源と、地質構造に結合されると共に、少なくとも２つの直交する方向で受信された波の成分を検出するために適した少なくとも１つの受信機と、地震波の送出−受信サイクルの実施を制御し、各受信機によって受信された信号を記録する制御装置とを有する地質的な構造を地震で探査する装置に関する。
【００３６】
この装置は、それらを中央制御、および、記録研究室またはコンピュータセンタに伝達する前に信号を予め処理する局所集録ユニットに含まれる上述したような楕円波識別装置を有することを特徴とする。
【００３７】
本発明による楕円波フィルタは、地面の各場所に１つの２成分受信機を使用することによって表面波のエネルギーろ波を行い、地震トレースを時間−スケールブロックに分解する。この方法は、選択された各受信点に３成分受信機によって１つ、または、いくつかのジオホンストリングスを備えた、通常使用されている重い地表受信装置を置換する方法を提供する。
【００３８】
地表の３軸受信機によってつくられる信号に関して提案した本発明の方法で得られたろ波品質は、非常に長い長さにわたってジオホンストリングの記録を積み重ねることによって従来のように得ることができるものよりも良好である。それは楕円波の真の楕円性を考慮し、真の楕円性は、その偏波の性質に関して簡単な仮説が立てられることなく、実際には非常に変化する。楕円波のろ波は、高速の表面波の減衰を可能にする。また、１つの多成分受信機によって記録することは、距離をおいた記録の積み重ねによる高周波の減衰を回避する。
【００３９】
【発明の実施の形態】
以下、本発明の一実施例を図面を参照して詳細に説明する。
【００４０】
本発明による方法は、楕円波を検出し識別するために上述したこの従来の技術分野の信号波変換技術の開発に関する。
【００４１】
その効果に対して所定の期間（すなわち、１サイクルの）の１／４でサンプルされた成分Ｘ、および、Ｚの信号波変形ＴＯ（Ｘ）、および、ＴＯ（Ｚ）を決定した後、楕円波の存在の可能性は、（基本である）楕円波の通過を通してＺのｎ番目のサンプルのＸの（ｎ−１）番目のサンプルに対する変換比率
ｒ＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）並びに
−ＴＯn（Ｚ）／ＴＯ(N+1)（Ｘ）
を計算することによって調べることができる。実際にこの比ｒはこの楕円波の楕円率ａ／ｂに等しいことをチェックすることができる。
【００４２】
従って、比ｒの値がこのウインドウｆに含まれる波の信号波変換を選択することによって、その楕円率が所定のウインドウｆに含まれる楕円波の存在を検出することが可能になる。
【００４３】
表面波のない地震波の再構成は、比
ｒ＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）または
ｒ＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ）
がウインドウｆに含まれる時間−スケールブロックを除去することによって行うことができる。比が決定されたウインドウにある時間−スケールブロックは表面波の影響を評価することができる。この影響は、最初のシスモグラムから減算することができる。
【００４４】
次に、フィルタ構造について記述する。
【００４５】
楕円波は次のようにろ波することができる。
【００４６】
１）２つまたは３つの成分の記録の成分Ｘ、および、Ｚの波の信号波変換を行う。
【００４７】
２）ろ波すべき波の楕円率を画定するセレクションウインドウを選択する。
【００４８】
３）各トラックに関して（すなわち、分析波のスケール係数の各値において）また各時間サンプル毎に（すなわち、この波の各移動による）次の比を計算する。
【００４９】
ｒ1＝ＴＯn（Ｚ）／ＴＯ(n-1)（Ｘ）
ｒ2＝−ＴＯn（Ｚ）／ＴＯ(n+1)（Ｘ）
４）有効化によるサンプルの選択
ａ）ｒ1が選択ウインドウに含まれるサンプルＴＯ（Ｚ）n、および、
ＴＯ(n-1)（Ｘ）(n-1)、および、
ｂ）ｒ2が選択ウインドウに含まれるサンプルＴＯn（Ｚ）、および、
ＴＯ(n-1)（Ｘ）を確認によって、および、無効サンプルをゼロにすることによってサンプルを選択する。
【００５０】
５）逆の信号波変換によって楕円波フィールドを再構成する。
【００５１】
元のシスモグラムから得られる楕円波フィールドを減算することは楕円波のないシスモグラムを得るためにはそれで十分である。
【００５２】
本発明を実施する装置は、少なくとも１つの多成分音響または地震波受信機Ｒと組み合わされる（図１，２参照）。この受信機Ｒは、少なくとも２つのピックアップを有し、その軸線は、検出される楕円偏波平面に関連して選択される２つの直交する方向を向いている。例えば、地震探査において、通常、軸線の１つが地震源、および、受信機結合場所を含む平面に配置された３軸の受信機が使用される。この装置は、受信機Ｒによって生成された種々の信号をサンプリング、および、デジタル化するアナログデジタル変換器２を含む処理装置１を有する。垂直方向のピックアップの軸線を含む平面の可能な包囲性の変位と、使用する地震源の場所、および、受信機Ｒの場所を含む垂直方向の基準平面に関する水平方向のピックアップの軸線を含む平面の可能な方位角の変位を計算に入れるために、処理装置１は、この基準平面で座標Ｘ，Ｚを計算する計算エレメント３と受信された波の間の楕円波を識別し、上述したアルゴリズムによってそれらを削除するフィルタプロセッサ４とを有する。実際には、プログラムされたコンピュータのフィルタプロセッサ４は、識別アルゴリズム段を実行するために使用されることが好ましい。
【００５３】
各受信機Ｒによって生成された信号の集録、および、処理の機能のすべてを満足するために、本発明によれば、フランス特許出願に示されたような局所集録ユニットＲＴＵ（図２参照）を使用することが有利である。各地域的な局所集録ユニットＲＴＵｉ（図２参照）は、少なくとも１つの地表地震受信機Ｒｉ，Ｒｊに隣接してそれに接続して配置される。管理プロセッサに加えて、種々の地震信号の組み合わせ、データ圧縮、地震設備（地震受信機、および、電子装置等）の展開段階中にリアルタイムの品質の制御のような伝送の前に、これらの局所集録ユニットＲＴＵは、リアルタイム延長機能を実行するためにプログラムされた（ＤＳＰタイプの）補完的なコンピュータプロセッサを有する。各局所集録ユニットＲＴＵはケーブル、および、無線装置を通してそれらを中央制御、および、記録研究室に搬送する前に地震源Ｓによって地面に搬送される地震波に応じて層の構造から受信された地震データを獲得して予め処理する。
【００５４】
本発明による方法を実施するために適した追加の処理プログラムは、この場合、このような環境において少なくとも１つの多軸受信機Ｒに接続されたこのような地表ユニットに含まれる。
【００５５】
また処理装置は、間に配置されたハイパスフィルタ（図示せず）を有する。
【００５６】
また中央の地上の研究室か、または、計算センターのいずれかにおいて、地上で獲得した後、各受信機の３つの成分の内少なくとも２つの成分を記録し、ろ波することができる。
【００５７】
次に、実験上のバリデーションについて説明する。地震の記録は地震装置で行われ、その地震装置は９個の相互接続されたジオホン、または、３Ｃを備えた１つの３軸地震受信機を有する従来のジオホンストリングと、重量低下地震源と、２ミリ秒ステップ集録アッセンブリとを有する。
【００５８】
楕円波ろ波の適用例：
図３乃至図６は、本発明の方法を適用した信号波変換の逆変換よるシスモグラムの再構成の結果を示す。この変換は時間の１／４の期間でサンプル化される。  半オクターブのサンプリングのスケールは、２４ｍｓのベースサイクルの１／４の十分な再構成、および、２．５オクターブの分析信号圧縮を行う。良好な再構成を行うために、サンプリング間隔は、最も狭い信号の期間の１／８に等しいことが望ましい。
【００５９】
次に、楕円波の判定について述べる。
【００６０】
表面波を構成する周波数バンドは、有効な信号よりさらに制限される。これらの波は、制限された数のオクターブで判定され、楕円波の判定結果は元のシスモグラムから減算される。
【００６１】
周波数分析は異なる周波数バンドでフィルタを並置することによって実行される。楕円分析は、分離楕円ウインドウを有するシスモグラムをフィルタすることによって実行することができる。
【００６２】
図４（Ａ）〜（Ｄ）は分析は５７０ｍ／ｓでの到達が主として（−０．８と２．２）の範囲内の楕円率のウインドウでフィルタされたＺ成分のシスモグラムであり、それぞれ、（−２．２，−０．８），（−０．８，０），（０，０．８），（０．８，２．２）の楕円率のウインドウのフィルタによるものである。負の値は、これらの楕円波が粒子の時計方向すなわち天体運行方向であるプログレイドモーションを生じることを意味する。これらはセサワ波として知られているＭ２波である。  また、説明した楕円波の主軸は水平方向にある。これらの分離したフィルタは種々の記録された波の楕円特性を分析することを可能にしてその性質を特定することを可能にする。
【００６３】
次に、地表フィルタとの比較について述べる。。
【００６４】
図５（Ａ）乃至（Ｃ）において、それぞれ、ジオホンストリングで得られる記録と、１つの３−Ｃピックアップの成分Ｚと、楕円波フィルタからの信号が示されており、例えば、５７０ｍ／ｓの仮想レイリ波ジオホンストリングによっては短か過ぎるので、ろ波されないことを示している。また、これは２２００ｍ／ｓの表面波の場合である。楕円波フィルタは、これらのドメインでさらに有効である。さらに、反射フィールドの周波数の内容は、１つの３−Ｃピックアップから得られるろ波シスモグラムでさらに広く見える。ろ波シスモグラム５（Ｃ）において、低速の波が残り、その偏波は実際にはリニヤである。
【００６５】
この波は一定の仮想期間を有し、波の時間−周波数フィルタによって容易に除去される。
【００６６】
図６（Ａ）乃至（Ｄ）は高周波のイベントによって示される空気波で、使用される楕円波フィルタによって減衰される尺度を示している。（Ａ）は未処理成分Ｚのシスモグラム、（Ｂ）は楕円フィルタ経由後のシスモグラム、（Ｃ）は信号波の時間−周波数フィルタ経由後のシスモグラム、（Ｄ）はロ−パスフィルタ経由後のシスモグラムである。
【００６７】
これは上述した例によって確認されるように、もしソース受信機方向、および、垂直方向の記録が利用できるのであれば、楕円表面波、レイリ波、および、仮想レイリ波は楕円波フィルタによってろ波することができる。
【００６８】
地表の３軸受信機Ｒによって製造された信号で提出された小波フィルタ方法の結果は、２０メートルにわたって広がる９個のジオホンによって行われる記録を積み重ねることによる従来のものより良好である。楕円波フィルタは、高速度表面波の減衰を行うことができ、１つのピックアップの記録は、スペースアウト記録によって高周波の減衰を防止することができる。
【００６９】
本発明による方法は、地質学的な用途の範囲内で、他の波に関して楕円特性の波を識別するために、例えば、直線偏波、または、他の楕円波を識別することができる。このような楕円波は、表面波またはチューブ波である。
【００７０】
地質学的なデータの処理に適用される実施モードを説明した。本発明による楕円波の識別、および、ろ波方法は、もちろん電子、電気物理的な信号に音響または地震波とは別の他の波のタイプ他の分野の用途にも使用することができる。
【００７１】
【発明の効果】
以上説明したように本発明は、少なくとも２つの直交方向で波の成分を少なくとも１つの受信機で検出し、それに応じて発生する信号の信号波変換の決定と、その各振幅の比に基づいた基準によって振幅を選択することにより、地震波の楕円表面波をろ波検出でき、この楕円波のフィールドを受信信号から減ずることにより、容易、かつ、正確に地震波の主要波を分離できる効果がある。
【図面の簡単な説明】
【図１】本発明の媒体内伝搬楕円波ろ波方法が適用された装置の主要部のブロック図である。
【図２】本方法を実施するように設計された装置が受信した信号を予め処理する局所集録ユニットの配置の説明図である。
【図３】図３は従来の２つのシスモグラムの間の差を計算することによって得られるシスモグラムで、（Ａ）は元のシスモグラム、（Ｂ）は楕円波、（Ｃ）は他のシスモグラムを示す図である。
【図４】図４は分離された楕円率ウインドウにろ波されたＺ成分のシスモグラムを示す図で、（Ａ）は（−２．２，−０．８），（Ｂ）は（−０．８，０），（Ｃ）は（０，０．８）、（Ｄ）は（０．８，２．２）の楕円率のウインドウのろ波を経たものを示す。
【図５】図５は、本発明による楕円波ろ波の結果を示す図で、（Ａ）はグラウンドフィルタを形成する９つのジオホンのストリングで得られたシスモグラム、（Ｂ）は１つの３軸受信機で得られるシスモグラム、（Ｃ）は他の第３のシスモグラムを示す図である。
【図６】図６はＺ成分のシスモグラムであって、（Ａ）は処理前のシスモグラム、（Ｂ）は楕円率にろ波されたときのシスモグラム、（Ｃ）は小波の時間−周波数ろ波を経た後のシスモグラム、（Ｄ）は空中波のろ波後のシスモグラムを示す図である。
【符号の説明】
１  処理装置
２  アナログデジタル変換器（Ａ／Ｄコンバータ）
３  計算エレメント
４  フィルタプロセッサ
Ｒ  受信機
ＲＴＵ  局所集録ユニット
Ｓ  地震源