JP3973113B2 - 楕円波の識別方法とその装置および地震探査システム - Google Patents
楕円波の識別方法とその装置および地震探査システム Download PDFInfo
- Publication number
- JP3973113B2 JP3973113B2 JP19767596A JP19767596A JP3973113B2 JP 3973113 B2 JP3973113 B2 JP 3973113B2 JP 19767596 A JP19767596 A JP 19767596A JP 19767596 A JP19767596 A JP 19767596A JP 3973113 B2 JP3973113 B2 JP 3973113B2
- Authority
- JP
- Japan
- Prior art keywords
- wave
- elliptic
- receiver
- received
- ellipticity
- 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
Links
- 238000000034 method Methods 0.000 title claims description 45
- 238000004458 analytical method Methods 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000005755 formation reaction Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 230000001902 propagating effect Effects 0.000 claims description 10
- 239000000463 material Substances 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 6
- 230000000644 propagated effect Effects 0.000 claims 1
- 230000033001 locomotion Effects 0.000 description 29
- 239000010410 layer Substances 0.000 description 16
- 230000009466 transformation Effects 0.000 description 16
- 238000001914 filtration Methods 0.000 description 15
- 238000006243 chemical reaction Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 5
- 239000002245 particle Substances 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000010287 polarization Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 229930091051 Arenine Natural products 0.000 description 1
- 241001441571 Hiodontidae Species 0.000 description 1
- 241000146377 Lestoros Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Filters And Equalizers (AREA)
Description
【発明の属する技術分野】
本発明は、多軸の受信機によって受信された電波の複数の軸に沿って測定された成分の組み合わされた処理によって、物質媒体で伝搬する電波の中から楕円波を識別することができるろ波方法に関する。
【0002】
【従来の技術】
本発明による方法は、種々の分野で用途を見いだすことができる。本発明は、地層に適用されたディスターバンス(変動)に応じて地層によって放射される波が3軸ジオホンすなわち"トリホン"のような地面に設置された受信機によって受信される地球物理学的な探査、および、地震探査に特に適用することができる。
【0003】
特に、地震探査において、表面波のような波は深く浸透しないことがよく分かっている。従って、それらは地球物理学的な探査にとって有効な情報は提供せず、ノイズと考慮される。それらのエネルギーは、有効な信号のエネルギーと比較して往々にして高い。従って表面波はろ波しなければならない。
【0004】
2つのタイプの表面波は弾性を有する一様な半無限媒体の表面に沿って伝搬する。2つのタイプの波とは、粒子の動きが伝搬平面に残るレイリ(Rayleigh)波と粒子の動きが伝搬平面に直角なラブ(Love)波とである。レイリ波としての粒子の動きは三角関数(レトログレード、反時計)方向をカバーする楕円に追随する(以下、この波を楕円波という。)。楕円の軸線は垂直方向であり、主軸と第2の軸との間の比は媒体の弾性パラメータに依存するが、1.4乃至1.9の範囲内になる。
【0005】
地層は、その表面の近くでウエザー層と称される低速度層(以下、この低速度層をウエザー層という。)を有し、この低速度層でフロア上の反射とこの低速度層の上部との間でいくつかの周波数について正干渉によって形成された(疑似レイリ波と称される)案内波(以下、この案内波を疑似レイリ波という。)が形成される。またこれらの案内波は、源と受信機を含む垂直方向の平面において伝搬する楕円波である。これらの波のいくつかについては粒子の動きは、レイリ波の特性として反時計回り(レトログレード)であり、他の波に関しては時計回り(プログレード)である。
【0006】
基本軸の比はウエザー層の厚さ、および、この層と下層との間の速度コントラストによって変化する。レイリ波とウエザー層の種々の案内波モードによって形成された全体は非常に分散的であり、明らかな速度は、周波数で変化し、広い範囲をカバーすることができる。その結果、表面波はサイズモグラム(地震動記録)の大きな部分をカバーし、それによって大きなエネルギーによる有効な情報を隠す。これらの波のすべては伝搬平面に含まれる。それらは楕円形であり、一般に直線的に極性化される直接、および、反射波から識別される。その基本軸は垂直方向、および、水平方向に接近している。
【0007】
本体波と同様に、表面波は過渡的である。疑似レイリ波モードの各々によるディスターバンスは、その特定の速度で受信点を通過する。作用、および、反作用の原理によって、ディスターバンスの時間積分はゼロであり、要するに、1つの点で記録されたディスターバンスはゼロ領域の小さい波である。
【0008】
これらを合計するためにレイリタイプの表面波は、認識基準を利用できるときに容易に配置することができる基本部分と、本体波とは異なり、認識基準を与える楕円部分と、受信点に対応するトレースでそれらの持続時間が1つまたはいくつかの振動に制限される過渡部分と、時間ウインドウで制限された周波数定数を有する分散部分とを有する。
さらに詳細には表面波については次の書類において示されている。
【0009】
1953年トルストイI.、ウスディンE.(Tolstoy I.,Usdin E.)による"層状に形成された弾性液体媒体の分散性特性:レイ理論"地球物理学第18巻844乃至870ページまたは1996年ムーニーH.、ボルトB.(Mooney H., Bolt B.)による"1つの表面層において第1の3つのレイリモードの分散特性"に記載されている。
【0010】
従来のろ波方法は、直列の相互接続されているジオホンの"ストリング"を地震受信機として使用し、その各信号を重ねることから成る。さらに有効には、表面波の最も長い波長に等しい距離または少なくとも2つのトレースの間の距離に等しい距離をカバーしなければならない。これは、数十メートルのケーブル、大きなスペースと重量、および、大きな位置決め、および、メインテナンス作業を必要とする。このろ波タイプは、表面波と反射信号との間の明らかな水平方向の速度差に基づいて位置決めの方向にのみ有効である。3Dにおいて使用可能な等方性の地面のフィルタの製造は実際に実施するには重すぎる。
【0011】
地面フィルタを使用した後にサイズモグラム(地震動記録)に残る表面波をろ波する他の公知の方法は、ある時間制限の後、すなわち、表面波の到着から、ジオホンによってピックアップされたものを無音とするか、ローカットフィルタを使用することによってある値以下のすべての周波数を削除することによることから成る。これらのフィルタは分離された周波数領域への記録の分解、および、ノイズのある周波数領域の削除に基づいている。
【0012】
複数の軸線の地震受信機を使用し、少なくとも2つの直交軸線に沿って信号を処理して表面波の識別を行うことができる公知の地震探査方法がある。
【0013】
例としての米国特許第US−A−4,458,341号は、レイリ波を検出するに適した処理装置と関連する2軸または3軸のピックアップによって地層からくる地震信号を受信する地震探査方法を示している。信号ピックアップの半径方向X、および、垂直方向V成分は、この処理装置のろ波、エンベロープ検出、および、比較回路に提供することができ、この処理装置は、その特別の特徴からレイリ波の受信ウインドウを決定して受信信号からの減算を制御することができる。
【0014】
米国特許第4,935,905号は、地層からくる波を受けるために2軸または3軸ジオホンを使用する地震探査方法を示す。ピックアップされる種々の成分は極座標で表される。レイリ波の成分は、平均値方法を適用することによってこれらの成分から得られ減算することができる。
【0015】
複数軸のジオホンによって表面波をろ波する方法は、1995年グラスゴーの57回EAEG会議の"拡張アブストラクト"で示されている。
【0016】
また特許第FR−A−2,648,567号は、電波分析技術を必要とする、例えば脳波の分析に適用される信号デジタル処理方法を示している。
【0017】
混合時間周波数ろ波:
周波数ろ波は、表面波(15Hzまでの最も低い周波数あるいはそれ以下の周波数)の周波数成分のフーリエドメインにおける削除、および、保存した周波数領域からの信号の再構成が続くトレースの高周波分解である。勿論これらの時間または周波数領域のノイズの削除は、これらの領域に属する有効な信号部分の削除を招く。削除すべきノイズから認識される領域と信号を引き出すことが望ましい場合残りの部分との間で重複を制限するために、2つの二次元の時間周波数分解を実行することは公知である。
【0018】
本発明による方法は、"ウェーブレット解析"と称される"時間−周波数"信号分析を使用し、その場合、信号は、幾何学数列によってサイズが変化するウェーブレットの配列によって信号が分析され、各ウェーブレットは、その長さの部分と等しいステップでトレースに沿って解釈される。この技術は、例えば、1984/1985年グーピロードP.とグロスマンA.とモーレットJ.(Goupillaud P., Grossmann A., Morlet J.)による地質探査第23巻の85乃至102ページの"地震信号分析におけるサイクルオクターブ、および、関連変換"に示されている。
【0019】
ウェーブレット変換の分析信号は、ゼロ平均値の信号である。例えば、リッカーインパルス
【0020】
【数1】
を使用することができる。
【0021】
ウェーブレット変換において、信号は一定の形状を有する。ウェーブレット変換による時間・周波数分解はディスクリート(サンプル)オペレーションである。与えられた寸法に関して、信号はせいぜい1つの期間(サイクル)に等しい一定のステップでシフトされる。信号のサイズは変化する。2(多くても、すなわち、1オクターブ)までの一定の幾何学的な比例によってスケール係数Eを変化することにより拡張または圧縮される。時間Tでのサンプル、および、ウェーブレット変換のスケール係数Eは、時間Tの中心のスケール係数Eのウェーブレットによるトレースの結果である。
【0022】
もしh(t)がウェーブレットのベースで、aがウェーブレットのスケール係数で、bによって変位されるとき、
【0023】
【数2】
ここで、aは、スケール係数またはパラメータであり、bは時間トランスレーションパラメータである。信号s(t)のウェーブレット変換は次のように表される。
【0024】
【数3】
逆ウェーブレット変換は、係数の範囲内について、対応するウェーブレットによるウェーブレット変換の各項の積の合計である。
【0025】
【数4】
【0026】
【課題を解決するための手段】
本発明による方法は、(例えば地層のような)物質媒体で伝搬する他の波から少なくとも1つの伝搬平面で伝搬する楕円波を識別することができる。それは、少なくとも2つの直交方向の波の成分に関する少なくとも1つの受信機(前記地層と結合されたジオホンまたは加速度計のような1つまたはいくつかの検出器を有する)による検出を含む。
【0027】
さらに、本方法は、受信された波に応じて受信機によって生じた信号のウェーブレット変換の決定、および、その各振幅の比に基づいた基準によってこれらのウェーブレット変換の振幅の選択とを有することを特徴とする。
【0028】
1つの実施例によれば、ろ波すべき楕円波の楕円率を有する楕円率の範囲を選択する段階と、選択された解析ウェーブレットの一連のスケール係数(f)と時間(T)とに対応したウェーブレット変換の一連の連続したサンプル(TOn(Z),TOn(X))を得るために、少なくとも2つの直行する方向で受信する信号のウェーブレット変換を行う段階と、
下記の比を決定する段階と、
r1=TOn(Z)/TO(n-1)(X),および,
r2=−TOn(Z)/TO(n+1)(X),
ここで(n−1,n,n+1)は連続整数、
前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する段階と、逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する段階とを有し、この再構成の楕円波フィールドは受信された信号から減算される。
【0029】
この識別方法は、表面波またはチューブ波と称される波をろ波するように地面の近傍または地中に1つまたはいくつかの、多軸受信機を設定することによって地層を伝搬する楕円波をろ波するために使用される。
【0030】
この識別方法は、音響または地震波が地層に出力され、地下の地層の地震探査方法に用いられ、地層の断層によって反射された波が少なくとも2つの異なる方向で、その波の成分を検出するために適した1つまたはいくつかの受信機によって受信され、ピックアップされた信号が記録され、その記録された信号は、受信された波に応じて各受信機に発生された信号のウェーブレット変換を決定し、かつ、各振幅の比に基づいた基準によってこのウェーブレット変換から生じる信号の振幅を選択することによって、楕円波を消去するために処理される。
【0031】
地震探査方法の1つの実施例によれば、ろ波すべき楕円波の楕円率を有する楕円率の範囲を選択する段階と、選択された解析ウェーブレットの一連のスケール係数(f)と時間(T)とに対応したウェーブレット変換の一連の連続したサンプル(TOn(Z),TOn(X))を得るために、少なくとも2つの直行する方向で受信する信号のウェーブレット変換を行う段階と、
下記の比を決定する段階と、
r1=TOn(Z)/TO(n-1)(X),および,
r2=−TOn(Z)/TO(n+1)(X),
ここで(n−1,n,n+1)は連続整数、
前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する段階と、逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する段階によって記録信号を処理し、再構成された楕円波フィールドを受信された信号から減算する。
【0032】
本発明による装置は、(例えば、地球物理学の応用における2軸または3軸のジオホンのような)少なくとも2つの直交する方向において波の成分を検出するのに適した、物質媒体(例えば、地層のような)と結合した少なくとも1つの受信機によって、物質媒体中を伝搬する音響または地震波から楕円波を識別することができる。
【0033】
それは、受信された波に応じて受信機に生じた信号のウェーブレット変換を決定し、および、その信号の各振幅の比に基づいた基準によってこれらのウェーブレット変換から生じる信号の振幅を選択するのに適した処理装置を有することを特徴とする。
【0034】
処理装置は、少なくとも2つの直行する方向で受信する信号のウェーブレット変換を行い、選択された解析ウェーブレットの一連のスケール係数(f)と時間(T)とに対応したウェーブレット変換の一連の連続したサンプル(TOn(Z),TOn(X))を得る装置と、比率
r1=TOn(Z)/TO(n-1)(X),および,
r2=−TOn(Z)/TO(n+1)(X),
ここで(n−1,n,n+1)は連続整数、
を決定する装置と、前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する装置と、逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する装置とを有し、再構成された楕円波フイールドは、その後、受信された信号から減算によって削除される。
【0035】
また、本発明は地層において波を発生する地震源と、地層に結合されると共に、少なくとも2つの直交する方向で受信された波の成分を検出するために適した少なくとも1つの受信機と、地震波の送出−受信サイクルの実施を制御し、各受信機によって受信された信号を記録する制御装置とを有する地層を地震で探査するシステムに関する。
【0036】
この装置は、それらを中央制御、および、記録研究室またはコンピュータセンタに伝達する前に信号を予め処理する局所集録ユニットに含まれる上述したような楕円波識別装置を有することを特徴とする。
【0037】
本発明による楕円波フィルタは、地面の各場所に1つの2成分受信機を使用することによって表面波のエネルギーろ波を行い、地震トレースを時間−スケールブロックに分解する。この方法は、選択された各受信点に3成分受信機によって1つ、または、いくつかのジオホンストリングスを備えた、通常使用されている重い地表受信装置を置換する方法を提供する。
【0038】
地表の3軸受信機によってつくられる信号に関して提案した本発明の方法で得られたろ波品質は、非常に長い長さにわたってジオホンストリングの記録を積み重ねることによって従来のように得ることができるものよりも良好である。それは楕円波の真の楕円性を考慮し、真の楕円性は、その偏波の性質に関して簡単な仮説が立てられることなく、実際には非常に変化する。楕円波のろ波は、高速の表面波の減衰を可能にする。また、1つの多成分受信機によって記録することは、距離をおいた記録の積み重ねによる高周波の減衰を回避する。
【0039】
【発明の実施の形態】
以下、本発明の一実施例を図面を参照して詳細に説明する。
【0040】
本発明による方法は、楕円波を検出し識別するために上述したこの従来の技術分野の信号波変換技術の開発に関する。
【0041】
その効果に対して所定の期間(すなわち、1サイクルの)の1/4でサンプルされた成分X、および、Zのウェーブレット変換TO(X)、および、TO(Z)を決定した後、楕円波の存在の可能性は、(基本である)楕円波の通過を通してZのn番目のサンプルのXの(n−1)番目のサンプルに対する変換比率
r=TOn(Z)/TO(n-1)(X)並びに
−TOn(Z)/TO(n+1)(X)
を計算することによって調べることができる。実際にこの比rはこの楕円波の楕円率a/bに等しいことをチェックすることができる。
【0042】
従って、比rの値がこの楕円率の範囲fに含まれる波のウェーブレット変換を選択することによって、その楕円率が所定の楕円率の範囲fに含まれる楕円波の存在を検出することが可能になる。
【0043】
表面波のない地震波の再構成は、比
r=TOn(Z)/TO(n-1)(X)または
r=−TOn(Z)/TO(n+1)(X)
が楕円率の範囲fに含まれる時間−スケールブロックを除去することによって行うことができる。比rが決定された楕円率の範囲にある時間−スケールブロックは表面波の影響を評価することができる。この影響は、最初のサイズモグラム(地震動記録)から減算することができる。
【0044】
次に、フィルタ構造について記述する。
【0045】
楕円波は次のようにろ波することができる。
【0046】
1)2つまたは3つの成分の記録の成分X、および、Zの波のウェーブレット変換を行う。
【0047】
2)ろ波すべき波の楕円率を画定する楕円率の範囲を選択する。
【0048】
3)各トラックに関して(すなわち、解析ウェーブレットのスケール係数の各値において)また各時間サンプル毎に(すなわち、この波の各移動による)次の比を計算する。
【0049】
r1=TOn(Z)/TO(n-1)(X)
r2=−TOn(Z)/TO(n+1)(X)
4)サンプルの確認による選択
a)r1が選択された楕円率の範囲に含まれるサンプルTOn(Z)、及び
、TO(n-1)(X)、および、
b)r2が選択された楕円率の範囲に含まれるサンプルTOn(Z)、及び
、TO(n+1)(X)を確認によって選択し、および、選択された楕円率
の範囲に含まれないサンプルをゼロにする。
【0050】
5)逆ウェーブレット変換によって楕円波フィールドを再構成する。
【0051】
オリジナルのサイズモグラム(地震動記録)から得られる楕円波フィールドを減算することは楕円波のないサイズモグラム(地震動記録)を得るためにはそれで十分である。
【0052】
本発明を実施する装置は、少なくとも1つの多成分音響または地震波受信機Rと組み合わされる(図1,2参照)。この受信機Rは、少なくとも2つのピックアップを有し、その軸線は、検出される楕円偏波平面に関連して選択される2つの直交する方向を向いている。例えば、地震探査において、通常、軸線の1つが地震源、および、受信機結合場所を含む平面に配置された3軸の受信機が使用される。この装置は、受信機Rによって生成された種々の信号をサンプリング、および、デジタル化するアナログデジタル変換器2を含む処理装置1を有する。垂直方向のピックアップの軸線を含む平面の可能な包囲性の変位と、使用する地震源の場所、および、受信機Rの場所を含む垂直方向の基準平面に関する水平方向のピックアップの軸線を含む平面の可能な方位角の変位を計算に入れるために、処理装置1は、この基準平面で座標X,Zを計算する計算エレメント3と受信された波の間の楕円波を識別し、上述したアルゴリズムによってそれらを削除するフィルタプロセッサ4とを有する。実際には、プログラムされたコンピュータのフィルタプロセッサ4は、識別アルゴリズム段階を実行するために使用されることが好ましい。
【0053】
各受信機Rによって生成された信号の集録、および、処理の機能のすべてを満足するために、本発明によれば、フランス特許出願に示されたような局所集録ユニットRTU(図2参照)を使用することが有利である。各地域的な局所集録ユニットRTUi(図2参照)は、少なくとも1つの地表地震受信機Ri,Rjに隣接してそれに接続して配置される。管理プロセッサに加えて、種々の地震信号の組み合わせ、データ圧縮、地震設備(地震受信機、および、電子装置等)の展開段階中にリアルタイムの品質の制御のような伝送の前に、これらの局所集録ユニットRTUは、リアルタイム延長機能を実行するためにプログラムされた(DSPタイプの)補完的なコンピュータプロセッサを有する。各局所集録ユニットRTUはケーブル、および、無線装置を通してそれらを中央制御、および、記録研究室に搬送する前に地震源Sによって地面に搬送される地震波に応じて層の構造から受信された地震データを獲得して予め処理する。
【0054】
本発明による方法を実施するために適した追加の処理プログラムは、この場合、このような環境において少なくとも1つの多軸受信機Rに接続されたこのような地表ユニットに含まれる。
【0055】
また処理装置は、間に配置されたハイパスフィルタ(図示せず)を有する。
【0056】
また中央の地上の研究室か、または、計算センターのいずれかにおいて、地上で獲得した後、各受信機の3つの成分の内少なくとも2つの成分を記録し、ろ波することができる。
【0057】
次に、実験上のバリデーションについて説明する。地震の記録は地震装置で行われ、その地震装置は9個の相互接続されたジオホン、または、3Cを備えた1つの3軸地震受信機を有する従来のジオホンストリングと、重量低下地震源と、2ミリ秒ステップ集録アセンブリとを有する。
【0058】
楕円波ろ波の適用例:
図3乃至図6は、本発明の方法を適用したウェーブレット変換の逆変換によるサイズモグラム(地震動記録)の再構成の結果を示す。この変換は時間の1/4の期間でサンプル化される。半オクターブのサンプリングのスケールは、24msのベースサイクルの1/4の十分な再構成、および、2.5オクターブの分析信号圧縮を行う。良好な再構成を行うために、サンプリング間隔は、最も狭い信号の期間の1/8に等しいことが望ましい。
【0059】
次に、楕円波の判定について述べる。
【0060】
表面波を構成する周波数バンドは、有効な信号よりさらに制限される。これらの波は、制限された数のオクターブで判定され、楕円波の判定結果はオリジナルのサイズモグラム(地震動記録)から減算される。
【0061】
周波数分析は異なる周波数バンドでフィルタを並置することによって実行される。楕円率分析は、分離された楕円率の範囲を有するサイズモグラム(地震動記録)をフィルタすることによって実行することができる。
【0062】
図4(A)〜(D)は分析は570m/sでの到達が主として(−0.8と2.2)の範囲内の楕円率の範囲でフィルタされたZ成分のサイズモグラム(地震動記録)であり、それぞれ、(−2.2,−0.8),(−0.8,0),(0,0.8),(0.8,2.2)の楕円率の範囲のフィルタによるものである。負の値は、これらの楕円波が粒子の時計方向の動き(天体運行方向の動き、プログレイドモーション)を生じることを意味する。これらはセザワ波として知られているM2波である。また、説明した楕円波の主軸は水平方向にある。これらの分離したフィルタは種々の記録された波の楕円特性を分析することを可能にしてその性質を特定することを可能にする。
【0063】
次に、地表フィルタとの比較について述べる。
【0064】
図5(A)乃至(C)において、それぞれ、ジオホンストリングで得られる記録と、1つの3−Cピックアップの成分Zと、楕円波フィルタからの信号が示されており、例えば、570m/sの疑似レイリ波ジオホンストリングによっては短か過ぎるので、ろ波されないことを示している。また、これは2200m/sの表面波の場合である。楕円波フィルタは、これらの領域でさらに有効である。さらに、反射フィールドの周波数の内容は、1つの3−Cピックアップから得られるろ波サイズモグラム(地震動記録)でさらに広く見える。ろ波サイズモグラム(地震動記録)5(C)において、低速の波が残り、その偏波は実際にはリニヤである。
【0065】
この波は一定の疑似期間を有し、ウェーブレット時間−周波数フィルタによって容易に除去される。
【0066】
図6(A)乃至(D)は高周波のイベントによって示される空気波で、使用される楕円波フィルタによって減衰される尺度を示している。(A)は未処理成分Zのサイズモグラム(地震動記録)、(B)は楕円波フィルタ経由後のサイズモグラム(地震動記録)、(C)はウェーブレット時間−周波数フィルタ経由後のサイズモグラム(地震動記録)、(D)はロ−パスフィルタ経由後のサイズモグラム(地震動記録)である。
【0067】
これは上述した例によって確認されるように、もしソース受信機方向、および、垂直方向の記録が利用できるのであれば、表面波の楕円波、レイリ波、および、疑似レイリ波は楕円波フィルタによってろ波することができる。
【0068】
地表の3軸受信機Rによって製造された信号で提出されたウェーブレットフィルタ方法の結果は、20メートルにわたって広がる9個のジオホンによって行われる記録を積み重ねることによる従来のものより良好である。ウェーブレットフィルタは、高速度表面波の減衰を行うことができ、1つのピックアップの記録は、スペースアウト記録によって高周波の減衰を防止することができる。
【0069】
本発明による方法は、地球物理学の応用の範囲内で、他の波に関して楕円特性の波を識別するために、例えば、直線偏波、または、他の楕円波を識別することができる。このような楕円波は、表面波またはチューブ波である。
【0070】
地球物理学的なデータの処理に適用される実施モードを説明した。本発明による楕円波の識別、および、ろ波方法は、もちろん電気、電子物理的な信号に、音響または地震波とは別の他の波のタイプに、及び他の分野の応用にも使用することができる。
【0071】
【発明の効果】
以上説明したように本発明は、少なくとも2つの直交方向で波の成分を少なくとも1つの受信機で検出し、それに応じて発生する信号のウェーブレット変換の決定と、その各振幅の比に基づいた基準によって振幅を選択することにより、地震波の表面波の楕円波をろ波検出でき、この楕円波のフィールドを受信信号から減ずることにより、容易、かつ、正確に地震波の主要波を分離できる効果がある。
【図面の簡単な説明】
【図1】 本発明の媒体内伝搬楕円波ろ波方法が適用された装置の主要部のブロック図である。
【図2】 本方法を実施するように設計された装置が受信した信号を予め処理する局所集録ユニットの配置の説明図である。
【図3】 図3は、(A)はオリジナルのサイズモグラム(地震動記録)、(B)は楕円波、(C)は(A)と(B)の2つのサイズモグラム(地震動記録)の間の差を計算することによって得られるサイズモグラム(地震動記録)を、示す図である。
【図4】 図4は分離された楕円率の範囲にろ波されたZ成分のサイズモグラム(地震動記録)を示す図で、(A)は(−2.2,−0.8),(B)は(−0.8,0),(C)は(0,0.8)、(D)は(0.8,2.2)の楕円率の範囲のろ波を経たものを示す。
【図5】 図5は、(A)はグラウンドフィルタを形成する9つのジオホンのストリングで得られたサイズモグラム(地震動記録)、(B)は1つの3軸受信機で得られるサイズモグラム(地震動記録)、(C)は本発明による楕円波ろ波の結果のサイズモグラム(地震動記録)を、示す図である。
【図6】 図6はZ成分のサイズモグラム(地震動記録)であって、(A)は処理前のサイズモグラム(地震動記録)、(B)は楕円波フィルタにろ波されたときのサイズモグラム(地震動記録)、(C)はウェーブレット時間−周波数フィルタを経た後のサイズモグラム(地震動記録)、(D)は空中波のろ波後のサイズモグラム(地震動記録)を、示す図である。
【符号の説明】
1 処理装置
2 アナログデジタル変換器(A/Dコンバータ)
3 計算エレメント
4 フィルタプロセッサ
R 受信機
RTU 局所集録ユニット
S 地震源
Claims (16)
- 1つの平面を形成する少なくとも2つの直交方向の波の成分に対応する信号を提供する少なくとも1つの多成分受信機(R)によって受信され、物質媒体中の少なくとも1つの伝搬平面で伝搬する、楕円波フィールドを識別する楕円波識別方法であって、
識別すべき楕円波の楕円率を画定する楕円率の範囲を選択する段階と、
解析ウェーブレットの異なるスケール係数と時間に対応したウェーブレット変換の一連のサンプルを形成するために、前記波の成分に対応する信号のウェーブレット変換を各々の波の成分に対して決定する段階と、
お互いに関して移動した前記一連のサンプルのサンプル振幅の比を形成することにより、選択された楕円率の範囲に依存する基準によって、各々の前記一連のサンプルの中からサンプルを選択する段階とを有する楕円波識別方法。 - 1つの平面を形成する少なくとも2つの直交方向(X、Z)の波の成分に対応する信号を提供する少なくとも1つの多成分受信機(R)によって受信され、物質媒体中の少なくとも1つの伝搬平面で伝搬する、楕円波フィールドを識別する楕円波識別方法であって、
識別すべき楕円波の楕円率を画定する楕円率の範囲を選択する段階と、
選択された解析ウェーブレットの一連のスケール係数(f)と時間(T)に対応したウェーブレット変換の一連の連続したサンプル(TO n (Z)、TO n (X))を得るために、少なくとも2つの直行方向(X、Z)で受信した信号のウェーブレット変換(TO(X)、TO(Z))を形成する段階と、
移動した前記サンプルの比
r1=TOn(Z)/TO(n-1)(X)、および、
r2=−TOn(Z)/TO(n+1)(X)、
ここでn−1,n,n+1は連続整数、
を決定する段階と、
前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する段階と、
逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する段階とを有する楕円波識別方法。 - 少なくとも1つの前記受信機は地層に結合されたジオホンまたは加速度計である複数の方向性検知器を含む請求項1に記載の楕円波識別方法。
- 少なくとも1つの前記受信機は地層に結合されたジオホンまたは加速度計である複数の方向性検知器を含む請求項2に記載の楕円波識別方法。
- 少なくとも1つの前記受信機は、楕円表面波をろ波するために地面の表面近傍に設定された少なくとも1つの受信機に結合されたジオホンまたは加速度計である複数の方向性検知器を含む請求項1に記載の楕円波識別方法。
- 少なくとも1つの前記受信機は、楕円表面波をろ波するために地面の表面近傍に設定された少なくとも1つの受信機に結合されたジオホンまたは加速度計である複数の方向性検知器を含む請求項2に記載の楕円波識別方法。
- 少なくとも1つの前記受信機は、少なくとも1つの井戸に沿って伝搬する楕円波をろ波するために、地層に堀込んだ少なくとも1つの井戸内に設定された少なくとも1つの受信機に結合されたジオホンまたは加速度計である複数の方向性検知器を含む請求項2に記載の楕円波識別方法。
- 前記再構成の楕円波フィールドを少なくとも1つの前記受信機により受信された信号から減算する請求項1に記載の楕円波識別方法。
- 前記再構成の楕円波フィールドを少なくとも1つの前記受信機により受信された信号から減算する請求項2に記載の楕円波識別方法。
- 地下の地層に対して音響または地震波を出力する段階と、
少なくとも2つの異なる方向で受信した波の成分を検出する少なくとも1つの多成分受信機(R)によって地下の地層の断層によって反射された楕円波を受信する段階と、
受信された波に応じて少なくとも1つの受信機によって受信した信号を記録する段階と、
少なくとも1つの伝搬平面で伝搬する、受信した楕円波を識別する段階とを有し、
前記識別する段階は、
識別すべき楕円波の楕円率を画定する楕円率の範囲を選択する段階と、
解析ウェーブレットの異なるスケール係数と時間に対応したウェーブレット変換の一連のサンプルを形成するために、前記波の成分に対応する信号のウェーブレット変換を各々の波の成分に対して決定する段階と、
お互いに関して移動した前記一連のサンプルのサンプル振幅の比を形成することにより、選択された楕円率の範囲に依存する基準によって、各々の前記一連のサンプルの中からサンプルを選択する段階とを有する地下の地層の地震探査方法。 - 地下の地層に対して音響または地震波を出力する段階と、
少なくとも2つの異なる直交方向(X、Z)で受信した波の成分を検出する少なくとも1つの多成分受信機(R)によって地下の地層の断層によって反射された楕円波を受信する段階と、
受信された波に応じて少なくとも2つの直行方向で少なくとも1つの多成分受信機(R)によって受信した信号を記録する段階と、
2つの異なる直交方向により形成された少なくとも1つの伝搬平面で伝搬する、受信した楕円波を識別する段階とを有し、
前記識別する段階は、
識別すべき楕円波の楕円率を画定する楕円率の範囲を選択する段階と、
解析ウェーブレットの異なるスケール係数と時間に対応したウェーブレット変換の一連の連続したサンプル(TO n (Z)、TO n (X))を形成するために、前記波の成分に対応する信号のウェーブレット変換(TO(X)、TO(Z))を各々の波の成分に対して決定する段階と、
移動した前記サンプルの比
r1=TOn(Z)/TO(n-1)(X)、および、
r2=−TOn(Z)/TO(n+1)(X)、
ここでn−1,n,n+1は連続整数、
を決定する段階と、
前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する段階と、
逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する段階とを有する地下の地層の地震探査方法。 - 前記再構成の楕円波フィールドを受信された楕円波から減算する段階を更に有する請求項11に記載の地震探査方法。
- 少なくとも2つの直交方向の波の成分を検出する物質媒体に結合された少なくとも1つの受信機によって受信され、物質媒体中で伝搬する音響または地震波の中から少なくとも1つの伝搬平面で伝搬する楕円波を識別する楕円波識別装置であって、
解析ウェーブレットの異なるスケール係数と時間に対応したウェーブレット変換の一連のサンプルを形成するために、前記波の成分に対応する信号のウェーブレット変換を各々の波の成分に対して決定し、
お互いに関して移動した各々の異なる前記一連のサンプルのサンプル振幅の比を形成することにより、選択された楕円率の範囲に依存する基準によって、各々の前記一連のサンプルの中からサンプルを選択する、処理アセンブリを有する楕円波識別装置。 - 前記処理アセンブリは、
解析ウェーブレットの異なるスケール係数と時間に対応したウェーブレット変換の一連の連続したサンプル(TO n (Z)、TO n (X))を形成するために、前記波の成分に対応する信号のウェーブレット変換(TO(X)、TO(Z))を各々の波の成分に対して決定する手段と、
移動した前記サンプルの比
r1=TOn(Z)/TO(n-1)(X)、および、
r2=−TOn(Z)/TO(n+1)(X)、
ここでn−1,n,n+1は連続整数、
を決定する手段と、
前記比の値r1,r2が選択された楕円率の範囲にあるサンプルを選択する手段と、
逆ウェーブレット変換を行うことによって楕円波フィールドを再構成する手段とを有する請求項13に記載の楕円波識別装置。 - 前記再構成の楕円波フィールドを受信された楕円波から減算する手段を更に有する請求項13に記載の楕円波識別装置。
- 前記再構成の楕円波フィールドを受信された楕円波から減算する手段を更に有する請求項14に記載の楕円波識別装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR95/09197 | 1995-07-26 | ||
FR9509197A FR2737308B1 (fr) | 1995-07-26 | 1995-07-26 | Methode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu |
Publications (3)
Publication Number | Publication Date |
---|---|
JPH0968578A JPH0968578A (ja) | 1997-03-11 |
JP3973113B2 true JP3973113B2 (ja) | 2007-09-12 |
JP3973113B6 JP3973113B6 (ja) | 2008-07-02 |
Family
ID=
Also Published As
Publication number | Publication date |
---|---|
NO963116D0 (no) | 1996-07-25 |
US5781502A (en) | 1998-07-14 |
EP0756182B1 (fr) | 2004-07-14 |
NO963116L (no) | 1997-01-27 |
FR2737308A1 (fr) | 1997-01-31 |
EP0756182A1 (fr) | 1997-01-29 |
JPH0968578A (ja) | 1997-03-11 |
RU2165093C2 (ru) | 2001-04-10 |
FR2737308B1 (fr) | 1997-09-19 |
DE69632892T2 (de) | 2004-10-28 |
CA2182111C (fr) | 2007-04-03 |
DE69632892D1 (de) | 2004-08-19 |
CA2182111A1 (fr) | 1997-01-27 |
NO318319B1 (no) | 2005-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7660200B2 (en) | Method for coherence-filtering of acoustic array signal | |
RU2523734C2 (ru) | Система и способ сбора сейсмических данных | |
EP2375268B1 (en) | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers | |
CA2658307C (en) | Method for deghosting marine seismic streamer data with irregular receiver positions | |
US4794573A (en) | Process for separating upgoing and downgoing events on vertical seismic profiles | |
US9188688B2 (en) | Flexural wave attenuation | |
US5781502A (en) | Method and device for filtering elliptical waves propagating in a medium | |
Keilis-Borok | Recording, identification, and measurement of surface wave parameters | |
Stone | Wavelet estimation | |
US4601022A (en) | Seismic exploration using non-impulsive vibratory sources activated by stationary, Gaussian codes, and processing the results in distortion-free final records particularly useful in urban areas | |
RU96115391A (ru) | Способ и устройство для селекции эллиптических волн, распространяющихся в среде | |
Marzetta et al. | A hydrophone vertical seismic profiling experiment | |
Satoh et al. | Nonlinear behaviour of scoria soil sediments evaluated from borehole records in eastern Shizuoka prefecture, Japan | |
JP3973113B6 (ja) | 楕円波の識別方法とその装置および地震探査システム | |
CN115826039A (zh) | 一种时间切片分类模型训练方法、系统及应用方法、系统 | |
CA2610997C (en) | Method for coherence-filtering of acoustic array signal | |
Iio | Seismic source spectrum of microearthquakes | |
US3736556A (en) | Method and apparatus for increasing seismic signal-to-noise ratio | |
Filson | Array seismology | |
Xu et al. | Ore-body imaging by crosswell seismic waveform inversion: A case study from Kambalda, Western Australia | |
Mammo et al. | Imaging of the weathered zone and estimation of Q in sediments | |
Pant et al. | A multicomponent offset VSP scale model investigation | |
Khusnani et al. | Relocation Study of Flores Sea Hypocenter (Mw= 7.3) Based on Single Station Estimation Using ObsPy | |
Tura et al. | Crosswell seismic reflection/diffraction tomography: A reservoir characterization application | |
Rogers et al. | A comparison of ground response in the Los Angeles region from nuclear explosions and the 1971 San Fernando earthquake |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051130 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060228 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060705 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060920 |
|
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: 20070509 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070608 |
|
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: 20100622 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |