JP2659461B2 - 反響抑制のための改良された方法 - Google Patents
反響抑制のための改良された方法Info
- Publication number
- JP2659461B2 JP2659461B2 JP7506413A JP50641395A JP2659461B2 JP 2659461 B2 JP2659461 B2 JP 2659461B2 JP 7506413 A JP7506413 A JP 7506413A JP 50641395 A JP50641395 A JP 50641395A JP 2659461 B2 JP2659461 B2 JP 2659461B2
- Authority
- JP
- Japan
- Prior art keywords
- signature
- pressure
- sine
- noise
- sign
- 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 - Lifetime
Links
- 238000000034 method Methods 0.000 title claims description 34
- 230000001629 suppression Effects 0.000 title description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 21
- 239000002245 particle Substances 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 9
- 238000005311 autocorrelation function Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000005314 correlation function Methods 0.000 claims description 2
- 238000002592 echocardiography Methods 0.000 claims description 2
- 238000013016 damping Methods 0.000 claims 1
- 238000007670 refining Methods 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 13
- 230000003044 adaptive effect Effects 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 238000010304 firing Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 101150112492 SUM-1 gene Proteins 0.000 description 2
- 101150096255 SUMO1 gene Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002596 correlated effect 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
- 230000007423 decrease Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/42—Devices characterised by the use of electric or magnetic means
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Measuring Volume Flow (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【発明の詳細な説明】 この発明は、対になった地震センサにより受信された
データを処理するための改良された方法を教示するもの
であり、この対になったセンサは種類が異なる。さらに
詳細には、この方法は比較的浅い海域での海洋地震探査
時に時々見られる水柱反響干渉の抑制に関連する。
データを処理するための改良された方法を教示するもの
であり、この対になったセンサは種類が異なる。さらに
詳細には、この方法は比較的浅い海域での海洋地震探査
時に時々見られる水柱反響干渉の抑制に関連する。
海洋地震調査の過程で、音響センサの長いひもが水中
を牽引される。定期的に、音響ソースは、反射地震波の
場を生成すべく、水面下の地層に音波をあて、その反射
地震波の場は、センサによって検出されて電気信号に変
換され、この電気信号は、ケーブルを通して船に送信さ
れる。対象とされた水面下の地層の地形図を表示するた
めに、信号は記録されて後に処理される。
を牽引される。定期的に、音響ソースは、反射地震波の
場を生成すべく、水面下の地層に音波をあて、その反射
地震波の場は、センサによって検出されて電気信号に変
換され、この電気信号は、ケーブルを通して船に送信さ
れる。対象とされた水面下の地層の地形図を表示するた
めに、信号は記録されて後に処理される。
海洋環境では、反射波の場は、それ自身を流体力学的
圧力の変動としてまたは粒子速度の変動として表現す
る。センサは異なる種類であってもよい。すなわち、そ
れらは圧力変動を測定する(水中聴音器)ように、また
は粒子速度を測定する(地中聴検器)ように設計されて
よい。場における克服すべき問題点によっては、一方ま
たは他方または両方の型式のセンサが、同じひも上に設
置されてもよい。
圧力の変動としてまたは粒子速度の変動として表現す
る。センサは異なる種類であってもよい。すなわち、そ
れらは圧力変動を測定する(水中聴音器)ように、また
は粒子速度を測定する(地中聴検器)ように設計されて
よい。場における克服すべき問題点によっては、一方ま
たは他方または両方の型式のセンサが、同じひも上に設
置されてもよい。
25〜200フィートのオーダの比較的浅い海では、いわ
ゆるベイ・ケーブルが使用される。深い海で行われるよ
うに船の後ろで連続的に牽引される代わりに、ここでは
センサひもは海底に直に据えつけられる。慣習的に、ジ
ンバル装着された地中聴検器が、海床の粒子速度を測定
するために選択されるセンサである。数百個または数千
個ものセンサが配備されることもあり、センサの電気的
出力は適切なデータ記録計へ多重化される。
ゆるベイ・ケーブルが使用される。深い海で行われるよ
うに船の後ろで連続的に牽引される代わりに、ここでは
センサひもは海底に直に据えつけられる。慣習的に、ジ
ンバル装着された地中聴検器が、海床の粒子速度を測定
するために選択されるセンサである。数百個または数千
個ものセンサが配備されることもあり、センサの電気的
出力は適切なデータ記録計へ多重化される。
図1は、海水本体14の下の海底16に置いた海底ケーブ
ル12を取り扱う船10を示す。周知のように海面18は、反
射性の空気/水界面を形成する。海底もまた通常は音響
インピーダンスに依存する良好な反射物である。速度セ
ンサ(地中聴検器)20および圧力センサ22(水中聴音
器)は海底に一緒に配置されて、ケーブル12内の別々の
電気チャネルに接続され、このチャネルを通してその信
号が船10上の記録保管・処理装置(図示せず)に送信さ
れる。単純化するためにセンサを2個だけ示す。定期的
に、ソース13は音波の場15を発生し、それは大地17へと
伝播し、そこで水面下の地層で反射され、24のような反
射波の場として戻ってくる。ソース13は、船10によって
または別の発砲用ボート(図示せず)によって発火され
る。
ル12を取り扱う船10を示す。周知のように海面18は、反
射性の空気/水界面を形成する。海底もまた通常は音響
インピーダンスに依存する良好な反射物である。速度セ
ンサ(地中聴検器)20および圧力センサ22(水中聴音
器)は海底に一緒に配置されて、ケーブル12内の別々の
電気チャネルに接続され、このチャネルを通してその信
号が船10上の記録保管・処理装置(図示せず)に送信さ
れる。単純化するためにセンサを2個だけ示す。定期的
に、ソース13は音波の場15を発生し、それは大地17へと
伝播し、そこで水面下の地層で反射され、24のような反
射波の場として戻ってくる。ソース13は、船10によって
または別の発砲用ボート(図示せず)によって発火され
る。
上昇する圧縮性反射波の場の到来24は海底の地中聴検
器20に打ち当たり、業界基準により、図2の26に示すよ
うに正方向の電気的インパルスを生成する。反射事象は
上昇を継続し、空気/水界面18に打ち当たり、そこで下
方へと再反射され、180゜の位相反転の後、ゴースト反
射として上面から地中聴検器に打ち当たる。通常、地中
聴検器の上面に印加される圧縮性パルスは負方向の電気
的インパルスを生成する。しかし、水面での位相反転の
ために、第一のゴースト反射は、図2で任意の振幅スケ
ールについてブロットした第二の正方向の電気的インパ
ルス30を生成する。海底面と水面のなめらかさの音響イ
ンピーダンスに依存して、理髪店内の相対する鏡の中で
多重反射が見られるのと全く同様に、波の場は海底面と
水面の間を何回も行ったり来たりして反射(反響)する
ことがある。第二のゴースト反射即ちマルチプル32は負
のパルスである。以降のマルチプルは、交互に正極性と
負極性を繰り返す。干渉性の多重反射は、記録されたデ
ータを台無しにする。反響効果は、潮位の低い海岸の砂
地や流砂のような湿地において観測されてきた。
器20に打ち当たり、業界基準により、図2の26に示すよ
うに正方向の電気的インパルスを生成する。反射事象は
上昇を継続し、空気/水界面18に打ち当たり、そこで下
方へと再反射され、180゜の位相反転の後、ゴースト反
射として上面から地中聴検器に打ち当たる。通常、地中
聴検器の上面に印加される圧縮性パルスは負方向の電気
的インパルスを生成する。しかし、水面での位相反転の
ために、第一のゴースト反射は、図2で任意の振幅スケ
ールについてブロットした第二の正方向の電気的インパ
ルス30を生成する。海底面と水面のなめらかさの音響イ
ンピーダンスに依存して、理髪店内の相対する鏡の中で
多重反射が見られるのと全く同様に、波の場は海底面と
水面の間を何回も行ったり来たりして反射(反響)する
ことがある。第二のゴースト反射即ちマルチプル32は負
のパルスである。以降のマルチプルは、交互に正極性と
負極性を繰り返す。干渉性の多重反射は、記録されたデ
ータを台無しにする。反響効果は、潮位の低い海岸の砂
地や流砂のような湿地において観測されてきた。
水中聴音器22は上昇する圧縮性反射波の場の到来を圧
縮性のパルスとして観測する。再び業界基準により、水
中聴音器は圧縮性パルスを、図2と同じ任意の振幅スケ
ールでプロットされた図3の36のような正方向の電気的
インパルスに変換する。水面からの第一のゴースト反射
は、疎密パルスとして下方へと進行し、水中聴音器22
は、これを負方向の電気的インパルス40として観測す
る。海底面と水面との間の第二の反射による第二のゴー
スト反射42は正である。以降のマルチプルは、交互に負
極性と正極性を示す。
縮性のパルスとして観測する。再び業界基準により、水
中聴音器は圧縮性パルスを、図2と同じ任意の振幅スケ
ールでプロットされた図3の36のような正方向の電気的
インパルスに変換する。水面からの第一のゴースト反射
は、疎密パルスとして下方へと進行し、水中聴音器22
は、これを負方向の電気的インパルス40として観測す
る。海底面と水面との間の第二の反射による第二のゴー
スト反射42は正である。以降のマルチプルは、交互に負
極性と正極性を示す。
パルス間の時間間隔τは、見かけの水深に水中減速度
(1/v)を乗算し2倍した積である。図1の28や38のよ
うな垂直入射に対しては、見かけの深さは真の水深であ
る。そうでなく、波の場44とパス45に関しては、見かけ
の深さは水深と水中減速度の2倍の積を入射角度αのコ
サイン(余弦)で除算したものである。
(1/v)を乗算し2倍した積である。図1の28や38のよ
うな垂直入射に対しては、見かけの深さは真の水深であ
る。そうでなく、波の場44とパス45に関しては、見かけ
の深さは水深と水中減速度の2倍の積を入射角度αのコ
サイン(余弦)で除算したものである。
このように、地中聴検器が第一の反射波の到来と第一
のゴーストの到来を同じ極性の電気的インパルスとして
観測するのに対して、水中聴音器は第一の反射波の到来
と第一のゴーストの到来を逆極性の電気的インパルスと
して観測する。したがって、一緒に配置した別々の種類
の2個のセンサを使用してそれらのサインを結合するこ
とにより、建設的に第一の到来エネルギーを強化するこ
とができ、かつ以降の反響性のマルチプルを破壊的に相
殺することができる。第二の利点は、別々のセンサ型式
により観測されるようなランダム・ノイズが必ずしも同
一性のあるものまたは相関性のあるものではないという
点に現れる。
のゴーストの到来を同じ極性の電気的インパルスとして
観測するのに対して、水中聴音器は第一の反射波の到来
と第一のゴーストの到来を逆極性の電気的インパルスと
して観測する。したがって、一緒に配置した別々の種類
の2個のセンサを使用してそれらのサインを結合するこ
とにより、建設的に第一の到来エネルギーを強化するこ
とができ、かつ以降の反響性のマルチプルを破壊的に相
殺することができる。第二の利点は、別々のセンサ型式
により観測されるようなランダム・ノイズが必ずしも同
一性のあるものまたは相関性のあるものではないという
点に現れる。
水中聴音器の信号と地中聴検器の信号との併合を成功
させるには、一方のセンサ型式からの信号が他方の型式
からの信号に適合するようにスケーリングされる必要が
ある。それらの伝達関数は、位相、振幅、周波数、減衰
において整合しなければならない。
させるには、一方のセンサ型式からの信号が他方の型式
からの信号に適合するようにスケーリングされる必要が
ある。それらの伝達関数は、位相、振幅、周波数、減衰
において整合しなければならない。
最も簡単明瞭なスケーリング方法は、2個のセンサか
らの記録された地震データのセットを分離し、機器の利
得効果を除去し、両方のデータ・セットに対して同一の
振幅復元関数を使用して、信号の振幅をバランスさせる
ことである。それぞれの第一到来の間の振幅比率がスケ
ール・ファクタである。しかし、ノイズや他の擾乱の存
在がその方法を過度に単純化させている。
らの記録された地震データのセットを分離し、機器の利
得効果を除去し、両方のデータ・セットに対して同一の
振幅復元関数を使用して、信号の振幅をバランスさせる
ことである。それぞれの第一到来の間の振幅比率がスケ
ール・ファクタである。しかし、ノイズや他の擾乱の存
在がその方法を過度に単純化させている。
1981年2月24日に付与され、この発明の譲受人に譲渡
された米国特許第4,253,164号の中でE.M.Hall.Jr.は、
加速度計または水中聴音器の伝達関数を地中聴検器の伝
達関数に整合される電気的ネットワークを教示してい
る。
された米国特許第4,253,164号の中でE.M.Hall.Jr.は、
加速度計または水中聴音器の伝達関数を地中聴検器の伝
達関数に整合される電気的ネットワークを教示してい
る。
1992年11月10日にF.J.Barrに付与された米国特許第5,
163,028号は、水中聴音器を地中聴検器に整合させるた
めの伝達関数を決定論的に計算し、且つセンサの海底へ
の結合上の不完全さを補償するアルゴリズムを教示して
いる。この方法は、選択された一緒に配置されたセンサ
の直上に位置する発砲から発生される、特殊な校正波の
場の使用を必要とする。校正発砲は、通常、地震探査を
始めるすぐ直前に発火される。
163,028号は、水中聴音器を地中聴検器に整合させるた
めの伝達関数を決定論的に計算し、且つセンサの海底へ
の結合上の不完全さを補償するアルゴリズムを教示して
いる。この方法は、選択された一緒に配置されたセンサ
の直上に位置する発砲から発生される、特殊な校正波の
場の使用を必要とする。校正発砲は、通常、地震探査を
始めるすぐ直前に発火される。
1990年12月18日にやはりF.J.Barrに付与された米国特
許第4,979,150号に、水柱反響を減ずるために地中聴検
器と共に配置された水中聴音器で使用するスケール・フ
ァクタを導き出すためのコンピュータ・プログラムが教
示されている。センサは海底より上の水中に一点に配置
してもよく、あるいは海底に配置することもできる。水
または海底物質の音響インピーダンスに由来するスケー
ル・ファクタは、決定論的にまたは統計的に計算するこ
とができる。前者の場合、前述の参考文献と同様にセン
サの直上で発砲される校正発砲を使用する。統計的な方
法には、決定論的方法に必要とされる発砲とセンサの位
置への特別な要求条件はない。統計的な方法は、それぞ
れのセンサにより記録された波列の自己相関関数および
相互相関関数の種々の組合せの間の比を繰り返し計算す
ることを含む。光路の角度に基づく、波の場の指向性に
関する補正が、導入されてもよい。
許第4,979,150号に、水柱反響を減ずるために地中聴検
器と共に配置された水中聴音器で使用するスケール・フ
ァクタを導き出すためのコンピュータ・プログラムが教
示されている。センサは海底より上の水中に一点に配置
してもよく、あるいは海底に配置することもできる。水
または海底物質の音響インピーダンスに由来するスケー
ル・ファクタは、決定論的にまたは統計的に計算するこ
とができる。前者の場合、前述の参考文献と同様にセン
サの直上で発砲される校正発砲を使用する。統計的な方
法には、決定論的方法に必要とされる発砲とセンサの位
置への特別な要求条件はない。統計的な方法は、それぞ
れのセンサにより記録された波列の自己相関関数および
相互相関関数の種々の組合せの間の比を繰り返し計算す
ることを含む。光路の角度に基づく、波の場の指向性に
関する補正が、導入されてもよい。
後に導入される概念である、適応性ノイズ・フィルタ
リングは、Bernard Widrow他によりProceedig of the I
EEE,v.63.n.12,1975年12月,1692〜1716頁で発表され
た、“Adaptive Noise Canceling:Principles and Appl
ications"の中で説明されている。
リングは、Bernard Widrow他によりProceedig of the I
EEE,v.63.n.12,1975年12月,1692〜1716頁で発表され
た、“Adaptive Noise Canceling:Principles and Appl
ications"の中で説明されている。
スケール・ファクタを計算する上記の既知の方法に対
する反対理由には、第一にランダム・ノイズの問題が含
まれる。特に、速度センサ(地中聴検器)は水中聴音器
に比較して非常にノイズが多い。統計的手法に関して
は、地中聴検器ノイズが相関処理の結果を必然的に汚染
する。第二に、相互相関比または自己相関比からスケー
ル・ファクタを決定する繰り返しプロセスにおいて、最
も適したスケール・ファクタを特定する収束係数が何ら
示唆されない。さらに、Barrにより教示されたような2
個のセンサからのデータの単純合計は、合計結果のSN比
を水中聴音器のそれに比較して低減させ得る。
する反対理由には、第一にランダム・ノイズの問題が含
まれる。特に、速度センサ(地中聴検器)は水中聴音器
に比較して非常にノイズが多い。統計的手法に関して
は、地中聴検器ノイズが相関処理の結果を必然的に汚染
する。第二に、相互相関比または自己相関比からスケー
ル・ファクタを決定する繰り返しプロセスにおいて、最
も適したスケール・ファクタを特定する収束係数が何ら
示唆されない。さらに、Barrにより教示されたような2
個のセンサからのデータの単純合計は、合計結果のSN比
を水中聴音器のそれに比較して低減させ得る。
「水中速度」または「速度サイン」に対して繰り返し
参照がなされる。別の方法で限定しない限り、「速度」
という用語は、そこを地震波の場が通過することによっ
て引き起こされる、媒体(大地または水)の粒子速度を
意味する簡略語である。速度という用語が他の意味で使
用される場合は、そのように限定されるだろう。圧力サ
イン、ノイズ・サイン、速度サインは、時間の関数とし
てそれぞれのパラメータの振幅レベルの量的な変動を参
照する。集合的な用語「ノイズ・サイン」は、あらゆる
望ましくない信号の全てを含むのに対して、「圧力」サ
インと「速度」サインは、求められている、探査に役立
つ地震信号である。
参照がなされる。別の方法で限定しない限り、「速度」
という用語は、そこを地震波の場が通過することによっ
て引き起こされる、媒体(大地または水)の粒子速度を
意味する簡略語である。速度という用語が他の意味で使
用される場合は、そのように限定されるだろう。圧力サ
イン、ノイズ・サイン、速度サインは、時間の関数とし
てそれぞれのパラメータの振幅レベルの量的な変動を参
照する。集合的な用語「ノイズ・サイン」は、あらゆる
望ましくない信号の全てを含むのに対して、「圧力」サ
インと「速度」サインは、求められている、探査に役立
つ地震信号である。
ノイズによってゆがめられない、水中聴音器と地中聴
検器のインパルス応答の間のスケール・ファクタを統計
的に決定する改良された方法への要求がある。さらに、
最適なスケール・ファクタを特定し、合計されたデータ
上の地中聴検器ノイズの影響を減じるために、収束係数
を明確に定義する方法に対する要求がある。
検器のインパルス応答の間のスケール・ファクタを統計
的に決定する改良された方法への要求がある。さらに、
最適なスケール・ファクタを特定し、合計されたデータ
上の地中聴検器ノイズの影響を減じるために、収束係数
を明確に定義する方法に対する要求がある。
この発明は、水柱反響を抑制するために地震データを
処理する改良された方法を提供する。この方法による
と、反射波の場は水域内で引き起こされる。反射波の場
は、埋め込まれたノイズを含む、圧力サインと速度サイ
ンとによって特徴づけられる。圧力サインと速度サイン
は、選択された併置された場所で同時に検出される。地
中聴検器のノイズ・サインは、速度サインと圧力サイン
とから適応的に推定される。ノイズ・サインは、きれい
な精製された速度サインを形成するために、逆極性で速
度サインへ加算される。選択されたスケール・ファクタ
が精製された速度サインに適用され、そして、圧力サイ
ンはスケーリングされ且つ精製された速度サインと合計
される。合計されたサインは、自己相関化される。
処理する改良された方法を提供する。この方法による
と、反射波の場は水域内で引き起こされる。反射波の場
は、埋め込まれたノイズを含む、圧力サインと速度サイ
ンとによって特徴づけられる。圧力サインと速度サイン
は、選択された併置された場所で同時に検出される。地
中聴検器のノイズ・サインは、速度サインと圧力サイン
とから適応的に推定される。ノイズ・サインは、きれい
な精製された速度サインを形成するために、逆極性で速
度サインへ加算される。選択されたスケール・ファクタ
が精製された速度サインに適用され、そして、圧力サイ
ンはスケーリングされ且つ精製された速度サインと合計
される。合計されたサインは、自己相関化される。
バリマックス関数が、合計サインのために計算され
る。上述のステップは何回も繰り返され、バリマックス
関数が最も密接に1に近づくまで各繰り返しにおいてス
ケール・ファクタを増加させる。
る。上述のステップは何回も繰り返され、バリマックス
関数が最も密接に1に近づくまで各繰り返しにおいてス
ケール・ファクタを増加させる。
この発明のある側面では、ノイズ抽出ステップは、埋
め込まれたノイズから信号を適応的に相殺するステップ
を含む。
め込まれたノイズから信号を適応的に相殺するステップ
を含む。
操作の機構および方法の両者に関してこの発明の特徴
であると信じられている新規な特徴は、その目的や利点
と共に、以下の詳細な説明と図面とからよりよく理解さ
れよう。ここで、この発明は説明だけの目的で例を通し
て説明されるが、この発明の限界を規定することを意図
したものではない。
であると信じられている新規な特徴は、その目的や利点
と共に、以下の詳細な説明と図面とからよりよく理解さ
れよう。ここで、この発明は説明だけの目的で例を通し
て説明されるが、この発明の限界を規定することを意図
したものではない。
図1は、浅い海底に据えつけられた異なる種類の2個
の地震センサの配置を示す。
の地震センサの配置を示す。
図2は、一次反射とゴースト反射への速度センサの応
答を示す。
答を示す。
図3は、一次反射とゴースト反射への圧力センサの応
答を示す。
答を示す。
図4は、適応性ノイズ回収論理と、水中聴音器から地
中聴検器への適応性フィルタとを図示する。
中聴検器への適応性フィルタとを図示する。
図5は、速度センサ信号を圧力センサ信号にスケーリ
ングする方法を定義するための操作シーケンスを示す。
そして、 図6は現在好ましい操作モードを図示している流れ図
である。
ングする方法を定義するための操作シーケンスを示す。
そして、 図6は現在好ましい操作モードを図示している流れ図
である。
一見したところ、以前に指摘したように、一緒に配置
された水中聴音器の振幅等化出力をもつ地中聴検器出力
の合計が多重反射を相殺するように思われる。問題の解
決はいくぶん複雑である。地中聴検器は、水中聴音器と
比較して本質的にノイズが多い。そのノイズは、何とか
して除去されなければならない。図2と図3で示したサ
インはディラック関数として理想化されているのに対し
て、実際には、信号は、センサの伝達関数を表している
過渡的波形を繰り込んだ、理想化された波形から成る。
された水中聴音器の振幅等化出力をもつ地中聴検器出力
の合計が多重反射を相殺するように思われる。問題の解
決はいくぶん複雑である。地中聴検器は、水中聴音器と
比較して本質的にノイズが多い。そのノイズは、何とか
して除去されなければならない。図2と図3で示したサ
インはディラック関数として理想化されているのに対し
て、実際には、信号は、センサの伝達関数を表している
過渡的波形を繰り込んだ、理想化された波形から成る。
引き続き考察を進める中で、データは、反射進行時間
の関数としてインデックスを付された、離散的なサンプ
ルの形であることが理解されるべきである。
の関数としてインデックスを付された、離散的なサンプ
ルの形であることが理解されるべきである。
地中聴検器の速度サインの最良の推定を適応的に行う
ために、水中聴音器の圧力サインをフィルタリングする
ことから始める。その後、最良に推定された地中聴検器
サインが、埋め込まれたノイズ・サインをも含む対応す
る実際の地中聴検器サインから引き算される。そのステ
ップは地中聴検器サインを相殺するという傾向があり、
それによってほぼ純粋のノイズ・サインを回収すること
を可能にする。その間に、少なくとも概念的に、オリジ
ナルの、実際のノイズの多い地中聴検器信号は、一時的
に、保管される。ノイズ・サインが回収された後、埋め
込まれたノイズを加えたオリジナルの実際の地中聴検器
サインは保管庫から取り出され、対応するノイズ・サイ
ンがそれから引き算され、きれいな精製された地中聴検
器信号が提供される。水中聴音器信号と合計されるきれ
いな地中聴検器信号に、スケール・ファクタを繰り返し
適用する。各繰り返し後に、合計が自己相関化される。
スケール・ファクタの大きさは、自己相関トレースのメ
イン・ローブと第一サイド・ローブの検討により適切に
決定される。適合性の大きさが最大にされるまで、繰り
返しは続く。
ために、水中聴音器の圧力サインをフィルタリングする
ことから始める。その後、最良に推定された地中聴検器
サインが、埋め込まれたノイズ・サインをも含む対応す
る実際の地中聴検器サインから引き算される。そのステ
ップは地中聴検器サインを相殺するという傾向があり、
それによってほぼ純粋のノイズ・サインを回収すること
を可能にする。その間に、少なくとも概念的に、オリジ
ナルの、実際のノイズの多い地中聴検器信号は、一時的
に、保管される。ノイズ・サインが回収された後、埋め
込まれたノイズを加えたオリジナルの実際の地中聴検器
サインは保管庫から取り出され、対応するノイズ・サイ
ンがそれから引き算され、きれいな精製された地中聴検
器信号が提供される。水中聴音器信号と合計されるきれ
いな地中聴検器信号に、スケール・ファクタを繰り返し
適用する。各繰り返し後に、合計が自己相関化される。
スケール・ファクタの大きさは、自己相関トレースのメ
イン・ローブと第一サイド・ローブの検討により適切に
決定される。適合性の大きさが最大にされるまで、繰り
返しは続く。
図4は、各データ・サンプルにサンプル時間インデッ
クスjの関数として作用する、適応性ノイズ回収ループ
の機構を概略的に示す。水中聴音器−地中聴検器変換フ
ィルタ47はブロック46と48を含み、これらは、各サンプ
ルに対して、周知のウィーナー・フィルタのような適応
性フィルタを構成する。フィルタは、振幅、位相、周波
数、減衰に関して速度サインの推定をもたらす。水中聴
音器の信号はフィルタ47のブロック46の入力であり、推
定された見かけの水深WDjとフィルタ振幅応答Ajとを繰
り込んで、推定された地中聴検器信号SUM1jを生成す
る。SUM1jは、合計連結点52でノイズの多い地中聴検器
入力トレースと減算的に結合させられる。出力は、ほぼ
純粋な地中聴検器ノイズεjであり、フィルタ47のブロ
ック50に入力される。フィードバック利得μを適用す
る。フィードバック利得は0から約0.01の狭い範囲内に
ある。利得が0である場合、フィルタは機能を停止す
る。利得があまりに大きい場合、不安定になる。したが
って、フィードバック利得は望ましい適応度を提供する
ように調整される。新しい水中聴音器−地中聴検器変換
フィルタを作成するために、見かけの水深が増加され
る。繰込みプロセスは、新しいフィルタを使用してブロ
ック48で繰り返される。オリジナルのフィルタとアップ
デートされたフィルタとの間の差異ΔWDが地中聴検器信
号の推定を改善する場合には、試行した水深増加は永久
的なエントリになる。
クスjの関数として作用する、適応性ノイズ回収ループ
の機構を概略的に示す。水中聴音器−地中聴検器変換フ
ィルタ47はブロック46と48を含み、これらは、各サンプ
ルに対して、周知のウィーナー・フィルタのような適応
性フィルタを構成する。フィルタは、振幅、位相、周波
数、減衰に関して速度サインの推定をもたらす。水中聴
音器の信号はフィルタ47のブロック46の入力であり、推
定された見かけの水深WDjとフィルタ振幅応答Ajとを繰
り込んで、推定された地中聴検器信号SUM1jを生成す
る。SUM1jは、合計連結点52でノイズの多い地中聴検器
入力トレースと減算的に結合させられる。出力は、ほぼ
純粋な地中聴検器ノイズεjであり、フィルタ47のブロ
ック50に入力される。フィードバック利得μを適用す
る。フィードバック利得は0から約0.01の狭い範囲内に
ある。利得が0である場合、フィルタは機能を停止す
る。利得があまりに大きい場合、不安定になる。したが
って、フィードバック利得は望ましい適応度を提供する
ように調整される。新しい水中聴音器−地中聴検器変換
フィルタを作成するために、見かけの水深が増加され
る。繰込みプロセスは、新しいフィルタを使用してブロ
ック48で繰り返される。オリジナルのフィルタとアップ
デートされたフィルタとの間の差異ΔWDが地中聴検器信
号の推定を改善する場合には、試行した水深増加は永久
的なエントリになる。
さて、図5を参照する。図5で文字Gを囲んでいる円
は、ノイズの多い地中聴検器入力信号(速度サイン)を
表す。文字Hを囲んでいる円は、対応する水中聴音器入
力信号(圧力サイン)を表す。図4の結合点48からの、
信号により汚染されていない、ほぼ純粋なノイズ・サイ
ンεjは、合計結合点54で、ノイズの多い地中聴検器入
力信号“G信号+ノイズ”に、逆の符号を付されて供給
され、もって、同図で“精製されたG信号”と呼ばれる
きれいな地中聴検器信号出力が提供される。
は、ノイズの多い地中聴検器入力信号(速度サイン)を
表す。文字Hを囲んでいる円は、対応する水中聴音器入
力信号(圧力サイン)を表す。図4の結合点48からの、
信号により汚染されていない、ほぼ純粋なノイズ・サイ
ンεjは、合計結合点54で、ノイズの多い地中聴検器入
力信号“G信号+ノイズ”に、逆の符号を付されて供給
され、もって、同図で“精製されたG信号”と呼ばれる
きれいな地中聴検器信号出力が提供される。
56で、スケール・ファクタが地中聴検器信号に適用さ
れ、合計結合点58で水中聴音器信号に加算され(適用シ
ーケンスは逆でも良い)、最終出力が生成される。スケ
ール・ファクタは、今説明されるようにブロック60で得
られる。
れ、合計結合点58で水中聴音器信号に加算され(適用シ
ーケンスは逆でも良い)、最終出力が生成される。スケ
ール・ファクタは、今説明されるようにブロック60で得
られる。
共通レシーバで分類され且つ振幅をバランスさせられ
た、対応する水中聴音器と地中聴検器とのデータセット
の所定の集合は、合理的にきれいな高水準信号によって
特徴づけられるトレースのウインドゥを選択する。選択
されたトレースによって表された、精製された地中聴検
器信号に対するスケール・ファクタを推定し、推定され
たスケール・ファクタを適用し、対応する水中聴音器ト
レースと合計する。合計を自己相関化する。スケール・
ファクタの望ましい範囲がカバーされるまで、各繰り返
しにおいてスケール・ファクタを増加して、収束係数を
計算し、繰り返す。集合全体に適用されるべきスケール
・ファクタとして、収束係数が最大になるような、すな
わち、1に最も密接に接近するような数値に対するスケ
ール・ファクタを選択する。
た、対応する水中聴音器と地中聴検器とのデータセット
の所定の集合は、合理的にきれいな高水準信号によって
特徴づけられるトレースのウインドゥを選択する。選択
されたトレースによって表された、精製された地中聴検
器信号に対するスケール・ファクタを推定し、推定され
たスケール・ファクタを適用し、対応する水中聴音器ト
レースと合計する。合計を自己相関化する。スケール・
ファクタの望ましい範囲がカバーされるまで、各繰り返
しにおいてスケール・ファクタを増加して、収束係数を
計算し、繰り返す。集合全体に適用されるべきスケール
・ファクタとして、収束係数が最大になるような、すな
わち、1に最も密接に接近するような数値に対するスケ
ール・ファクタを選択する。
収束係数は、バリマックス関数 から最も好都合に計算される。式中、xi,i=1,2,3,……
nは、自己相関関数の中央ローブと主要なサイド・ロー
ブとが期待される領域を含む、自己相関関数に沿った選
択された時間ゲート内での、離散的なデータ・サンプル
のシーケンスである。xが離散的なスパイクである場
合、図2または図3に示すように、VM=1である。波の
場の過渡現象がスパイク状でなくなるにつれて、VMの値
は減少する。
nは、自己相関関数の中央ローブと主要なサイド・ロー
ブとが期待される領域を含む、自己相関関数に沿った選
択された時間ゲート内での、離散的なデータ・サンプル
のシーケンスである。xが離散的なスパイクである場
合、図2または図3に示すように、VM=1である。波の
場の過渡現象がスパイク状でなくなるにつれて、VMの値
は減少する。
操作の現在好ましいモードは、図6の流れ図に要約さ
れている。水中聴音器と地中聴検器のデータ・セット
は、共通のレシーバ集合にそれぞれ分類される。データ
は、通常の態様での幾何学的な広がりに対し、逆の利得
で調整され、振幅をバランスさせられ、補正される。現
在部分的に処理されたデータはノイズ・バーストに対し
て処理され、f−kフィルタリングされ、サインは繰り
出される。前述したのはすべて、もちろん、業界基準の
手順である。
れている。水中聴音器と地中聴検器のデータ・セット
は、共通のレシーバ集合にそれぞれ分類される。データ
は、通常の態様での幾何学的な広がりに対し、逆の利得
で調整され、振幅をバランスさせられ、補正される。現
在部分的に処理されたデータはノイズ・バーストに対し
て処理され、f−kフィルタリングされ、サインは繰り
出される。前述したのはすべて、もちろん、業界基準の
手順である。
ブロック54で、図4と図5に関して説明したように、
適応性ノイズ相殺処理がデータ・セットに適用される。
ブロック61で、図5との関係で説明したように、相関係
数VMの最大値によって示されている最良推定スケール・
ファクタを使用して、データ・セットが互いに他方に対
してスケーリングされる。
適応性ノイズ相殺処理がデータ・セットに適用される。
ブロック61で、図5との関係で説明したように、相関係
数VMの最大値によって示されている最良推定スケール・
ファクタを使用して、データ・セットが互いに他方に対
してスケーリングされる。
この方法はベイ・ケーブル操作に関連して記述されて
きたが、どのような種類の海洋地震探査や、反響が湿気
の多い浜砂で遭遇されるような、いわゆる移行地帯にも
適用され得るということが理解されるべきである。地中
聴検器信号を水中聴音器のそれにスケーリングする操作
が、繰り返し言及されてきた。逆の処理がなされ得る、
すなわち、水中聴音器信号が地中聴検器信号に整合する
ようにスケーリングされてもよいということが理解され
るべきである。説明と例示の目的で特定の側面が説明さ
れたが、それはこの発明の範囲と精神とを限定するもの
ではない。
きたが、どのような種類の海洋地震探査や、反響が湿気
の多い浜砂で遭遇されるような、いわゆる移行地帯にも
適用され得るということが理解されるべきである。地中
聴検器信号を水中聴音器のそれにスケーリングする操作
が、繰り返し言及されてきた。逆の処理がなされ得る、
すなわち、水中聴音器信号が地中聴検器信号に整合する
ようにスケーリングされてもよいということが理解され
るべきである。説明と例示の目的で特定の側面が説明さ
れたが、それはこの発明の範囲と精神とを限定するもの
ではない。
Claims (15)
- 【請求項1】複数の反響を抑制するために地震データを
処理する改良された方法であって、 (a)水環境中に反射地震波の場を生成するステップで
あって、前記地震波の場は、圧力サインと、埋め込まれ
たノイズ・サインを含む速度サインとにより特徴付けら
れている、ものと、 (b)前記反射地震波の場の圧力サインを選択された位
置で検出するステップと、 (c)前記反射地震波の場の速度サインを、前記選択さ
れた位置に併置した場所で同時に検出するステップと、 (d)前記圧力サインをフィルタリングするステップ
と、 (e)前記フィルタリングされた圧力サインを前記速度
サインに減算的に結合することにより、前記ノイズ・サ
インを抽出するステップと、 (f)抽出されたノイズ・サインを前記速度サインに逆
の符号で加算することにより、精製された速度サインを
生成しするステップと、 (g)望ましいスケール・ファクタを前記精製された速
度サインに適用すると共に、前記圧力サインを前記精製
され且つスケーリングされた速度サインに合計するステ
ップと、 (h)合計されたサインに対する自己相関関数を生成す
るステップと、 を具備する方法。 - 【請求項2】請求項1で規定される方法であって、 自己相関関数用の収束係数を計算するステップと、 ステップ(g)からステップ(h)を繰り返すステップ
と、 前記収束係数が1に最も密接に接近するまで、各繰り返
しにおいて前記望ましいスケール・ファクタを増加する
ステップと、 を具備する方法。 - 【請求項3】請求項2で規定される方法であって、 抽出するステップが、前記埋め込まれたノイズ・サイン
から前記速度サインを適応的に相殺するステップを含む
方法。 - 【請求項4】請求項3で規定される方法であって、 ステップ(d)が、減算的に結合するステップに先立
ち、前記圧力サインを水深の関数としてフィルタリング
するステップを含む方法。 - 【請求項5】請求項4で規定される方法であって、 フィルタリングのステップが、望ましいフィードバック
利得ファクタを前記ノイズ・サインに導入するステップ
を含む方法。 - 【請求項6】請求項5で規定される方法であって、 前記フィードバック利得ファクタを調節して望ましい適
応度を提供するステップを具備する方法。 - 【請求項7】請求項2で規定される方法であって、 前記収束係数が、バリマックス関数によって定義される
方法。 - 【請求項8】請求項7で規定される方法であって、 前記バリマックス関数VMが、下式であって、 式中、x1=1,2,3,……nは、そのメイン・ローブおよび
第一のサイド・ローブを含む前記自己相関関数に沿った
選択された時間ゲート内から取られる離散的なデータ・
サンプルのシーケンスであるものから計算される方法。 - 【請求項9】地震データ・セットから複数の反響を抑制
する改良された方法であって、 (a)反射地震波の場を湿潤な環境の中で生成するステ
ップと (b)前記反射波の場の圧力サインを、前記環境の選択
された位置で測定するステップと、 (c)前記環境の中で実質的に同一の選択された位置
で、前記反射波の場の粒子速度サインを同時に測定する
ステップであって、前記粒子速度サインは不必要なノイ
ズ・サインを含んでいる、ものと、 (d)前記圧力サインをフィルタリングするステップ
と、 (e)前記フィルタリングされた圧力サインを前記粒子
速度サインに減算的に適用することにより、前記ノイズ
・サインを分離するステップと、 (f)前記分離されたノイズ・サインをそれから引き算
することによって前記粒子速度サインを精製するステッ
プと、 (g)精製された粒子速度サインを前記圧力サインに関
してスケーリングするステップと、 (h)前記圧力サインを前記粒子速度サインと合計する
ステップと、 (i)合計されたサインについて自己相関関数を生成す
るステップと、 (j)前記自己相関関数についての収束係数を計算する
ステップと、 を具備する方法。 - 【請求項10】請求項9で規定される方法であって、 ステップ(g)が、精製された粒子速度サインに関して
前記圧力サインをスケーリングするステップを含む方
法。 - 【請求項11】請求項10で規定される方法であって、 前記収束係数をバリマックス関数から計算するステップ
を具備する方法。 - 【請求項12】請求項9で規定される方法であって、 ステップ(g)〜ステップ(j)を繰り返すステップ
と、 前記相関係数が最も密接に1に接近するまで、各繰り返
しについて前記スケール・ファクタを増加するステップ
と、 を具備する方法。 - 【請求項13】請求項10で規定される方法であって、 ステップ(g)〜ステップ(j)を繰り返すステップ
と、 前記相関係数が最も密接に1に接近するまで、各繰り返
しの後に前記スケール・ファクタを増加するステップ
と、 を具備する方法。 - 【請求項14】請求項9で規定される方法であって、 ステップ(d)が、振幅、位相、周波数および減衰に関
して前記粒子速度サインを推定するために、前記圧力サ
インを変換するステップを含む方法。 - 【請求項15】請求項14で規定される方法であって、 ステップ(d)が、見かけの水深の関数として前記圧力
サインをフィルタリングするステップを含む方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/101,949 | 1993-08-04 | ||
US08/101,949 US5365492A (en) | 1993-08-04 | 1993-08-04 | Method for reverberation suppression |
US101,949 | 1993-08-04 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08501395A JPH08501395A (ja) | 1996-02-13 |
JP2659461B2 true JP2659461B2 (ja) | 1997-09-30 |
Family
ID=22287331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP7506413A Expired - Lifetime JP2659461B2 (ja) | 1993-08-04 | 1994-07-13 | 反響抑制のための改良された方法 |
Country Status (11)
Country | Link |
---|---|
US (1) | US5365492A (ja) |
EP (1) | EP0666992B1 (ja) |
JP (1) | JP2659461B2 (ja) |
AU (1) | AU671036B2 (ja) |
BR (1) | BR9405543A (ja) |
CA (1) | CA2145281C (ja) |
CO (1) | CO4290391A1 (ja) |
DE (1) | DE69420250D1 (ja) |
NO (1) | NO311779B1 (ja) |
RU (1) | RU2116657C1 (ja) |
WO (1) | WO1995004945A1 (ja) |
Families Citing this family (56)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9321125D0 (en) * | 1993-10-13 | 1993-12-01 | Geco As | Method of processing reflection data |
US5448531A (en) * | 1994-05-05 | 1995-09-05 | Western Atlas International | Method for attenuating coherent noise in marine seismic data |
US5442591A (en) * | 1994-06-21 | 1995-08-15 | Western Atlas International | Method for adaptively suppressing noise transients in summed co-sensor seismic recordings |
US5521881A (en) * | 1994-09-02 | 1996-05-28 | Exxon Production Research Company | Method of processing seismic data having multiple reflection noise |
US5774416A (en) * | 1995-04-07 | 1998-06-30 | Pgs, Tensor, Inc. | Method and device for attenuating water column reverberations using co-located hydrophones and geophones in ocean bottom seismic processing |
US5621699A (en) * | 1995-07-07 | 1997-04-15 | Pgs Ocean Bottom Seismic, Inc. | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring |
US5572483A (en) * | 1995-07-25 | 1996-11-05 | Western Atlas International, Inc. | Method of reducing noise in seismic signals by adaptive filtering of a noise reference |
US5724306A (en) * | 1995-12-29 | 1998-03-03 | Western Atlas International, Inc. | Method for correcting dual sensor data for imperfect geophone coupling using production seismic data |
US5737433A (en) * | 1996-01-16 | 1998-04-07 | Gardner; William A. | Sound environment control apparatus |
FR2743896B1 (fr) * | 1996-01-23 | 1998-04-10 | Geophysique Cie Gle | Procede de traitement de calibration d'une paire de capteurs hydrophone/geophone et procede de prospection sismique mettant en oeuvre ce traitement |
FR2743897B1 (fr) * | 1996-01-23 | 1998-04-10 | Geophysique Cie Gle | Procede de prospection sismique marine au moyen d'un couple de capteurs hydrophone et geophone |
US5754492A (en) * | 1996-02-12 | 1998-05-19 | Pgs Tensor, Inc. | Method of reverberation removal from seismic data and removal of dual sensor coupling errors |
US5696734A (en) * | 1996-04-30 | 1997-12-09 | Atlantic Richfield Company | Method and system for eliminating ghost reflections from ocean bottom cable seismic survey signals |
US5621700A (en) * | 1996-05-20 | 1997-04-15 | Schlumberger Technology Corporation, Geco-Prakla Div. | Method for attenuation of reverberations using a pressure-velocity bottom cable |
US5793702A (en) * | 1996-05-21 | 1998-08-11 | Western Atlas International, Inc. | Method for measuring the water bottom reflectivity |
GB9612470D0 (en) * | 1996-06-14 | 1996-08-14 | Geco As | Multiple attenuation method |
US5774417A (en) * | 1996-10-25 | 1998-06-30 | Atlantic Richfield Company | Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing |
US5818795A (en) * | 1996-10-30 | 1998-10-06 | Pgs Tensor, Inc. | Method of reduction of noise from seismic data traces |
GB9800741D0 (en) * | 1998-01-15 | 1998-03-11 | Geco As | Multiple attenuation of multi-component sea-bottom data |
GB2337591B (en) | 1998-05-20 | 2000-07-12 | Geco As | Adaptive seismic noise and interference attenuation method |
US6151275A (en) * | 1998-09-11 | 2000-11-21 | Pgs Tensor, Inc. | Method of dual wavefield reinforcement |
US6246637B1 (en) | 1998-09-28 | 2001-06-12 | Pgs Tensor, Inc. | Method and system for combining three component seismic data |
US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
US6725164B1 (en) * | 1999-03-17 | 2004-04-20 | Input/Output, Inc. | Hydrophone assembly |
GB9906456D0 (en) | 1999-03-22 | 1999-05-12 | Geco Prakla Uk Ltd | Method and system for reducing effects of sea surface ghost contamination in seismic data |
US20020118602A1 (en) | 2001-02-27 | 2002-08-29 | Sen Mrinal K. | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
CA2446987C (en) | 2001-05-25 | 2009-12-15 | Exxonmobil Upstream Research Company | Trapped water bottom multiple and peg-leg multiple suppression for ocean bottom seismic data |
US6738715B2 (en) | 2001-09-14 | 2004-05-18 | Exxonmobil Upstream Research Company | Method for attenuating noise in seismic data |
GB2389183B (en) * | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
US6884935B2 (en) | 2002-03-05 | 2005-04-26 | Robert H. Whidden | Method of transmitting electrical power |
US6903277B2 (en) | 2002-03-05 | 2005-06-07 | Robert H Whidden | Conduit for use in the transmission of electrical power |
RU2308743C2 (ru) * | 2002-03-20 | 2007-10-20 | Инпут/Аутпут, Инк. | Способ фильтрации сейсмических сигналов и устройство для его осуществления |
US7239577B2 (en) | 2002-08-30 | 2007-07-03 | Pgs Americas, Inc. | Apparatus and methods for multicomponent marine geophysical data gathering |
GB2415258B8 (en) * | 2002-08-30 | 2007-06-07 | Pgs Americas Inc | Method for multicomponent marine geophysical data gathering |
GB2395305B (en) * | 2002-11-15 | 2006-03-22 | Westerngeco Seismic Holdings | Processing seismic data |
US6894948B2 (en) * | 2003-01-29 | 2005-05-17 | Pgs Exploration (Uk) Limited | Method for processing dual sensor seismic data to attenuate noise |
US7561493B2 (en) * | 2003-05-30 | 2009-07-14 | Fairfield Industries, Inc. | Method and apparatus for land based seismic data acquisition |
US7310287B2 (en) | 2003-05-30 | 2007-12-18 | Fairfield Industries Incorporated | Method and apparatus for seismic data acquisition |
GB2405473B (en) * | 2003-08-23 | 2005-10-05 | Westerngeco Ltd | Multiple attenuation method |
GB2410551B (en) * | 2004-01-30 | 2006-06-14 | Westerngeco Ltd | Marine seismic acquisition system |
US8534959B2 (en) | 2005-01-17 | 2013-09-17 | Fairfield Industries Incorporated | Method and apparatus for deployment of ocean bottom seismometers |
US7319636B2 (en) * | 2005-03-14 | 2008-01-15 | Westerngeco, L.L.C. | Calibration of pressure gradient recordings |
US8127706B2 (en) * | 2005-05-02 | 2012-03-06 | Fairfield Industries Incorporated | Deck configuration for ocean bottom seismometer launch platforms |
US7466625B2 (en) * | 2006-06-23 | 2008-12-16 | Westerngeco L.L.C. | Noise estimation in a vector sensing streamer |
US8593907B2 (en) * | 2007-03-08 | 2013-11-26 | Westerngeco L.L.C. | Technique and system to cancel noise from measurements obtained from a multi-component streamer |
US20080270035A1 (en) * | 2007-04-26 | 2008-10-30 | Ahmet Kemal Ozdemir | System and Technique to Remove Perturbation Noise from Seismic Sensor Data |
US7957906B2 (en) * | 2008-05-07 | 2011-06-07 | Pgs Geophysical As | Method for attenuating low frequency noise in a dual-sensor seismic streamer |
US8611191B2 (en) * | 2008-05-22 | 2013-12-17 | Fairfield Industries, Inc. | Land based unit for seismic data acquisition |
US9025413B2 (en) | 2009-12-07 | 2015-05-05 | Pgs Geophysical As | Method for full-bandwidth source deghosting of marine seismic streamer data |
US8902699B2 (en) | 2010-03-30 | 2014-12-02 | Pgs Geophysical As | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers |
US9182512B2 (en) | 2011-10-26 | 2015-11-10 | Westerngeco L.L.C. | Processing multi-component seismic data |
US9448315B2 (en) * | 2011-12-27 | 2016-09-20 | Cgg Services Sa | Device and method for denoising ocean bottom data |
CN103487794B (zh) * | 2012-06-13 | 2016-01-06 | 中国科学院声学研究所 | 一种基于小波包变换的水底混响抑制方法 |
CA2910781A1 (en) * | 2013-04-29 | 2014-11-06 | Schlumberger Canada Limited | Deghosting with adaptive operators |
CN104597500B (zh) * | 2013-10-31 | 2017-05-10 | 中国石油天然气集团公司 | 一种水陆检波器地震数据匹配方法 |
US11249214B2 (en) * | 2018-05-30 | 2022-02-15 | Saudi Arabian Oil Company | Noise suppression of seafloor geophone seismic data |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4253164A (en) * | 1978-10-30 | 1981-02-24 | Western Geophysical Co. Of America | Multi-purpose seismic transducer |
US4910716A (en) * | 1989-01-31 | 1990-03-20 | Amoco Corporation | Suppression of coherent noise in seismic data |
EG19158A (en) * | 1989-08-25 | 1996-02-29 | Halliburton Geophys Service | System for attenuation of water-column reverberation |
US5235554A (en) * | 1991-03-11 | 1993-08-10 | Halliburton Geophysical Services, Inc. | Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations |
US5257241A (en) * | 1991-05-08 | 1993-10-26 | Atlantic Richfield Company | Method and system for acquisition of 3-dimensional marine seismic data |
US5163028A (en) * | 1991-09-27 | 1992-11-10 | Halliburton Geophysical Services, Inc. | Method for correcting impulse response differences of hydrophones and geophones as well as geophone coupling to the water-bottom in dual-sensor, bottom-cable seismic operations |
US5251183A (en) * | 1992-07-08 | 1993-10-05 | Mcconnell Joseph R | Apparatus and method for marine seismic surveying utilizing adaptive signal processing |
-
1993
- 1993-08-04 US US08/101,949 patent/US5365492A/en not_active Expired - Lifetime
-
1994
- 1994-07-13 RU RU95111130A patent/RU2116657C1/ru active
- 1994-07-13 BR BR9405543-2A patent/BR9405543A/pt not_active Application Discontinuation
- 1994-07-13 AU AU73280/94A patent/AU671036B2/en not_active Ceased
- 1994-07-13 EP EP94923411A patent/EP0666992B1/en not_active Expired - Lifetime
- 1994-07-13 WO PCT/US1994/007739 patent/WO1995004945A1/en active IP Right Grant
- 1994-07-13 DE DE69420250T patent/DE69420250D1/de not_active Expired - Lifetime
- 1994-07-13 JP JP7506413A patent/JP2659461B2/ja not_active Expired - Lifetime
- 1994-07-13 CA CA002145281A patent/CA2145281C/en not_active Expired - Fee Related
- 1994-07-29 CO CO94033344A patent/CO4290391A1/es unknown
-
1995
- 1995-03-31 NO NO19951270A patent/NO311779B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
DE69420250D1 (de) | 1999-09-30 |
RU95111130A (ru) | 1997-06-20 |
CO4290391A1 (es) | 1996-04-17 |
WO1995004945A1 (en) | 1995-02-16 |
AU7328094A (en) | 1995-02-28 |
EP0666992A1 (en) | 1995-08-16 |
NO951270D0 (no) | 1995-03-31 |
CA2145281A1 (en) | 1995-02-16 |
CA2145281C (en) | 1998-08-11 |
US5365492A (en) | 1994-11-15 |
EP0666992A4 (en) | 1998-04-29 |
AU671036B2 (en) | 1996-08-08 |
RU2116657C1 (ru) | 1998-07-27 |
NO951270L (no) | 1995-05-18 |
JPH08501395A (ja) | 1996-02-13 |
NO311779B1 (no) | 2002-01-21 |
EP0666992B1 (en) | 1999-08-25 |
BR9405543A (pt) | 1999-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2659461B2 (ja) | 反響抑制のための改良された方法 | |
Holland et al. | High-resolution geoacoustic inversion in shallow water: A joint time-and frequency-domain technique | |
EP0680616B1 (en) | Method for deriving water bottom reflectivity in dual sensor seismic surveys | |
EP0414344B1 (en) | Marine seismic reflection geophysical surveying | |
US5621699A (en) | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring | |
EP0515188B1 (en) | Method in marine seismics of removing multiple data | |
MXPA04000857A (es) | Metodo para procesar datos sismicos de sensor doble para atenuar ruido. | |
CA2491340A1 (en) | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers | |
WO2001001170A1 (en) | Dual sensor signal processing method for on-bottom cable seismic wave detection | |
US5991238A (en) | Weighted backus filter method of combining dual sensor traces | |
US5963507A (en) | Method for improving the accuracy of ocean bottom reflectivity estimations using the inverse backus filter | |
EP0851243B1 (en) | Measuring sea floor reflectivity | |
GB2328017A (en) | Seismic source signature estimation | |
Rauch | On the role of bottom interface waves in ocean seismo-acoustics: a review | |
GB2151359A (en) | Method for determining source and receiver statics in marine seismic exploration | |
Nicolas et al. | Geoacoustical parameters estimation with impulsive and boat-noise sources | |
EP0689064B1 (en) | Methods of adaptively reducing noise transients in co-sensor seismic signals | |
Barr et al. | A dual-sensor bottom-cable 3-D survey in the Gulf of Mexico | |
Holland et al. | High resolution geoacoustic inversion in shallow water: A joint time and frequency domain technique | |
Trabant | Single Channel Seismic Reflection Systems | |
IE52273B1 (en) | Method for seismic exploration by vertical seismic profiling and installation for its implementation |