JPH0854473A - 総和された共センサ地震記録における過渡ノイズを適応的に抑制する方法 - Google Patents
総和された共センサ地震記録における過渡ノイズを適応的に抑制する方法Info
- Publication number
- JPH0854473A JPH0854473A JP15466895A JP15466895A JPH0854473A JP H0854473 A JPH0854473 A JP H0854473A JP 15466895 A JP15466895 A JP 15466895A JP 15466895 A JP15466895 A JP 15466895A JP H0854473 A JPH0854473 A JP H0854473A
- Authority
- JP
- Japan
- Prior art keywords
- common
- geophone
- pressure
- signal
- aggregate
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 25
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 230000001052 transient effect Effects 0.000 claims abstract description 7
- 230000002159 abnormal effect Effects 0.000 claims abstract description 5
- 230000000712 assembly Effects 0.000 claims description 3
- 238000000429 assembly Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims 2
- 238000001514 detection method Methods 0.000 claims 1
- 239000000523 sample Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000001427 coherent effect Effects 0.000 description 3
- 241000283153 Cetacea Species 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 241000238557 Decapoda Species 0.000 description 1
- 241001417494 Sciaenidae Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- KWIUHFFTVRNATP-UHFFFAOYSA-N glycine betaine Chemical compound C[N+](C)(C)CC([O-])=O KWIUHFFTVRNATP-UHFFFAOYSA-N 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 230000000007 visual effect 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/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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【目的】 好ましくない異常な過渡現象に対処するノイ
ズ除去技術を提供する。 【構成】 共に設置された異なる種類のセンサからの地
震データは共通速度受信器集合体と共通圧力受信器集合
体とにフォーマットされる。速度サイン振幅と圧力サイ
ン振幅との間の第１の比は固定された解析ウインド内で
計測される。第２の比は、制限された範囲の重み付けさ
れた領域に関する圧力サインと速度サインとの間で計測
される。第１および第２の比は結合して演算子が形成さ
れる。重み付け領域ウインド内部の圧力サインおよび速
度サインは演算子の存在下で結合されて過渡現象のない
時間スケールのデータが規定される。
ズ除去技術を提供する。 【構成】 共に設置された異なる種類のセンサからの地
震データは共通速度受信器集合体と共通圧力受信器集合
体とにフォーマットされる。速度サイン振幅と圧力サイ
ン振幅との間の第１の比は固定された解析ウインド内で
計測される。第２の比は、制限された範囲の重み付けさ
れた領域に関する圧力サインと速度サインとの間で計測
される。第１および第２の比は結合して演算子が形成さ
れる。重み付け領域ウインド内部の圧力サインおよび速
度サインは演算子の存在下で結合されて過渡現象のない
時間スケールのデータが規定される。
Description
【０００１】
【産業上の利用分野】海老、ニベ科の魚、鯨などによっ
て発生される環境ノイズのような過渡ノイズにより汚染
された浅水の海底地震データを処理するための改良され
た方法に関するものである。
て発生される環境ノイズのような過渡ノイズにより汚染
された浅水の海底地震データを処理するための改良され
た方法に関するものである。
【０００２】
【従来の技術】海底地震探査の過程の間、数百もの地震
センサの長い列が適切なケーブルにより水中を牽引され
る。６〜８秒毎のような、周期的に、地下の地層を音響
探査すべく、音響源が波場を放出するためにトリガさ
れ、これにより、センサにより検出される反射音響波場
が発生される。反射波場を表す、検出された機械的な信
号は、センサにより電気的な信号に変換される。電気的
な信号はケーブル内の導体を介して牽引している船に伝
送され、ここで信号は地下の地層の地勢を表示するため
の将来の解析のために記録される。音響源の各トリガリ
ングからの地震信号は、デジタル化され、そして、１つ
の時間系列がケーブルに接続された各センサまたはセン
サグループに対応する、複数の離散的な時間系列を含む
組の形で記録される。音響源の各トリガリングから生じ
た離散的な時間系列の組は、単一の多トレース視覚地震
記録として表示される。これらのトレースは可変的な振
幅または可変的な密度または可変的に色付けされたトレ
ースの形をとる。
センサの長い列が適切なケーブルにより水中を牽引され
る。６〜８秒毎のような、周期的に、地下の地層を音響
探査すべく、音響源が波場を放出するためにトリガさ
れ、これにより、センサにより検出される反射音響波場
が発生される。反射波場を表す、検出された機械的な信
号は、センサにより電気的な信号に変換される。電気的
な信号はケーブル内の導体を介して牽引している船に伝
送され、ここで信号は地下の地層の地勢を表示するため
の将来の解析のために記録される。音響源の各トリガリ
ングからの地震信号は、デジタル化され、そして、１つ
の時間系列がケーブルに接続された各センサまたはセン
サグループに対応する、複数の離散的な時間系列を含む
組の形で記録される。音響源の各トリガリングから生じ
た離散的な時間系列の組は、単一の多トレース視覚地震
記録として表示される。これらのトレースは可変的な振
幅または可変的な密度または可変的に色付けされたトレ
ースの形をとる。
【０００３】２０〜２５０フィート程度の深さにおける
浅水では、所謂ベイケーブルが使用される。ある技術的
な理由により、一般的なベイケーブルは異なる種類の共
に設置されたセンサ対を含んでいる。通常は圧力センサ
が、ジンバル装着された粒子速度センサと対になってお
り、これらセンサの双方は習慣的に海底上に直接置かれ
る。浅水操作は公知であり、冗長な説明や付随した例示
が必要のない程、地球物理探査分野の当業者には日常的
なことである。
浅水では、所謂ベイケーブルが使用される。ある技術的
な理由により、一般的なベイケーブルは異なる種類の共
に設置されたセンサ対を含んでいる。通常は圧力センサ
が、ジンバル装着された粒子速度センサと対になってお
り、これらセンサの双方は習慣的に海底上に直接置かれ
る。浅水操作は公知であり、冗長な説明や付随した例示
が必要のない程、地球物理探査分野の当業者には日常的
なことである。
【０００４】地震記録上に現れる好ましくない波場はノ
イズであるとみなされる。多くの異なる音源が、一連の
浅水海底探査の間に受信される地震信号に影響を与える
ノイズを生成する。鯨の会話は地球物理学者には厄介な
ノイズであるが、環境学者の耳には音楽となる。
イズであるとみなされる。多くの異なる音源が、一連の
浅水海底探査の間に受信される地震信号に影響を与える
ノイズを生成する。鯨の会話は地球物理学者には厄介な
ノイズであるが、環境学者の耳には音楽となる。
【０００５】海底と水面との間でバウンドする音響信号
による残響は所望の反射波場を著しくゆがめてしまう。
米国特許出願第０８／１０１，９７９号であって、現在
はＷ．Ｈ．Ｄｒａｇｏｓｅｔに与えられ且つ本発明の譲
受人に譲渡された米国特許第５，３６５，４９２号であ
るもののような、そのような汚染を取り除くための種々
の方法が知られている。この特許には、互いに異なる種
類のセンサからの出力をスケーリングするための方法が
開示されている。さらに、総和の前にそれぞれの信号の
ノイズ成分を減じる、適応的なノイズ相殺処理が、教示
されている。
による残響は所望の反射波場を著しくゆがめてしまう。
米国特許出願第０８／１０１，９７９号であって、現在
はＷ．Ｈ．Ｄｒａｇｏｓｅｔに与えられ且つ本発明の譲
受人に譲渡された米国特許第５，３６５，４９２号であ
るもののような、そのような汚染を取り除くための種々
の方法が知られている。この特許には、互いに異なる種
類のセンサからの出力をスケーリングするための方法が
開示されている。さらに、総和の前にそれぞれの信号の
ノイズ成分を減じる、適応的なノイズ相殺処理が、教示
されている。
【０００６】競合する地震探査船により発生された地震
ショットはコヒーレントなノイズを生成する。Ｒ．Ａ．
Ｍａｒｓｃｈａｌｌ等に対して１９９０年６月２６日に
与えられた米国特許第４，９３７，７９４号には、共通
受信器集合体に対して共通ショット集合体をフォーマッ
トすることにより、地震記録上のコヒーレントなノイズ
を抑制する方法が説明されている。共通受信器集合体か
らの地震トレース対は、ディファランシャル・ノーマル
・ムーヴアウトに関して補正され、ＲＭＳ信号パワーに
逆比例した重み付けがされ、そして圧縮された共通受信
器集合体を発生させるために組み合わされる。圧縮され
た共通受信器集合体は、後続する処理用の共通中間点集
合体として再度フォーマットされる。
ショットはコヒーレントなノイズを生成する。Ｒ．Ａ．
Ｍａｒｓｃｈａｌｌ等に対して１９９０年６月２６日に
与えられた米国特許第４，９３７，７９４号には、共通
受信器集合体に対して共通ショット集合体をフォーマッ
トすることにより、地震記録上のコヒーレントなノイズ
を抑制する方法が説明されている。共通受信器集合体か
らの地震トレース対は、ディファランシャル・ノーマル
・ムーヴアウトに関して補正され、ＲＭＳ信号パワーに
逆比例した重み付けがされ、そして圧縮された共通受信
器集合体を発生させるために組み合わされる。圧縮され
た共通受信器集合体は、後続する処理用の共通中間点集
合体として再度フォーマットされる。
【０００７】近表面のばらつきによるノイズは、１９９
４年３月８日にＲｏｎａｌｄ Ｅ．Ｃｈａｍｂｅｒｓに
与えられ且つ本発明の譲受人に譲渡された米国特許第
５，２９３，３５２号の教示にしたがって取り除かれま
たは最小にされる。ここには、未処理の共通ショット集
合体の組が共通受信器集合体に対して再分類されること
が教示されている。共通受信器集合体からの波場エンベ
ロープは、移動させられたデータの組を得るために、近
表面速度の半分を使用して移動させられる。移動された
データの組は共通ショット集合体に戻されて再分類さ
れ、コヒーレントなノイズが減じられたデータの組を得
るために、元の未処理の共通ショット集合体から引かれ
る。
４年３月８日にＲｏｎａｌｄ Ｅ．Ｃｈａｍｂｅｒｓに
与えられ且つ本発明の譲受人に譲渡された米国特許第
５，２９３，３５２号の教示にしたがって取り除かれま
たは最小にされる。ここには、未処理の共通ショット集
合体の組が共通受信器集合体に対して再分類されること
が教示されている。共通受信器集合体からの波場エンベ
ロープは、移動させられたデータの組を得るために、近
表面速度の半分を使用して移動させられる。移動された
データの組は共通ショット集合体に戻されて再分類さ
れ、コヒーレントなノイズが減じられたデータの組を得
るために、元の未処理の共通ショット集合体から引かれ
る。
【０００８】当初は地上の地震振動データで使用される
ことを意図したものであるが、浅い海底の探査に適応す
ることができる方法では、有効な地震信号のレベルの基
準モデルが、地震データ記録サイクルの間の多数の時間
ウインドのそれぞれについて確立される。モデルは、平
均値の組を形成する多数の掃引から、各時間ウインドに
ついての絶対的な大きさの平均値を得ることで選択され
る。各組の中央値が選択され、且つ、各時間ウインドに
ついての基準モデルを提供するために、適切な係数によ
って調整される。次いで、記録が作られる。記録からの
各データサンプルは、サンプルを含む時間ウインドに対
応する基準モデルと比較される。データサンプルの振幅
値が基準モデルのレベルを超えた場合、サンプルは総和
の前に圧縮される。本発明の譲受人に譲渡された、Ｒａ
ｌｐｈ Ａ．Ｌａｎｄｒｕｍへの米国特許第４，３４
４，１５８号を参照のこと。
ことを意図したものであるが、浅い海底の探査に適応す
ることができる方法では、有効な地震信号のレベルの基
準モデルが、地震データ記録サイクルの間の多数の時間
ウインドのそれぞれについて確立される。モデルは、平
均値の組を形成する多数の掃引から、各時間ウインドに
ついての絶対的な大きさの平均値を得ることで選択され
る。各組の中央値が選択され、且つ、各時間ウインドに
ついての基準モデルを提供するために、適切な係数によ
って調整される。次いで、記録が作られる。記録からの
各データサンプルは、サンプルを含む時間ウインドに対
応する基準モデルと比較される。データサンプルの振幅
値が基準モデルのレベルを超えた場合、サンプルは総和
の前に圧縮される。本発明の譲受人に譲渡された、Ｒａ
ｌｐｈ Ａ．Ｌａｎｄｒｕｍへの米国特許第４，３４
４，１５８号を参照のこと。
【０００９】特に厄介なタイプのノイズ、つまり公知の
方法では取り扱うことが容易でないものは、生物学的な
ものである。図１はこのようなノイズの例を示したもの
で、これはＧｒａｎｄ Ｂｅｔａｉｌすなわちニベ科の
魚あるいは同様な魚によるものと思われる。ノイズバー
ストは、１０のような反転した矢印状の回折パターンで
ある。所望の波場は、１２により示した広く緩やかにカ
ーブした双曲線の波形である。
方法では取り扱うことが容易でないものは、生物学的な
ものである。図１はこのようなノイズの例を示したもの
で、これはＧｒａｎｄ Ｂｅｔａｉｌすなわちニベ科の
魚あるいは同様な魚によるものと思われる。ノイズバー
ストは、１０のような反転した矢印状の回折パターンで
ある。所望の波場は、１２により示した広く緩やかにカ
ーブした双曲線の波形である。
【００１０】
【発明が解決しようとする課題】本発明の目的は、好ま
しくない異常な過渡現象に対処するノイズ除去技術を提
供することにある。
しくない異常な過渡現象に対処するノイズ除去技術を提
供することにある。
【００１１】
【課題を解決するための手段】上記目的を達成するた
め、本発明によれば、複数の速度応答形および圧力応答
形の共に配置された共センサにより同時に検出された地
震信号の時間スケールの記録から好ましくない異常な過
渡現象を適応的に抑制する方法が提供される。反射波場
が水中環境内で発生され、速度および圧力の時間スケー
ルのサインを得るために、波場はそれぞれの共センサに
より検出される。固定された解析ウインドの内部では、
速度サインを圧力サインの振幅に対してスケーリングす
るための、スケール係数ｋが、反復的に規定される。圧
力応答形センサからのサインは、第１の共通圧力受信器
集合体を形成するためにグループ化される。検出された
速度応答形のサインは、第２の共通の速度受信器集合体
を形成するためにグループ化される。比Ｒ1 が、第２お
よび第１の集合体の間で規定される。それぞれの共通受
信器集合体は複数の重み付け領域に仕切られる。第１の
重み付け領域内において、比Ｒ2 が、当該領域内の圧力
サインと速度サインとの間で規定される。比Ｒ1 および
Ｒ2 は等化演算子を形成するために結合される。過渡現
象のない時間スケールのデータは、等化演算子の存在下
で共通圧力受信器重み付け領域内のデータを共通速度受
信器重み付け領域内のスケーリングされたデータと結合
することにより規定される。
め、本発明によれば、複数の速度応答形および圧力応答
形の共に配置された共センサにより同時に検出された地
震信号の時間スケールの記録から好ましくない異常な過
渡現象を適応的に抑制する方法が提供される。反射波場
が水中環境内で発生され、速度および圧力の時間スケー
ルのサインを得るために、波場はそれぞれの共センサに
より検出される。固定された解析ウインドの内部では、
速度サインを圧力サインの振幅に対してスケーリングす
るための、スケール係数ｋが、反復的に規定される。圧
力応答形センサからのサインは、第１の共通圧力受信器
集合体を形成するためにグループ化される。検出された
速度応答形のサインは、第２の共通の速度受信器集合体
を形成するためにグループ化される。比Ｒ1 が、第２お
よび第１の集合体の間で規定される。それぞれの共通受
信器集合体は複数の重み付け領域に仕切られる。第１の
重み付け領域内において、比Ｒ2 が、当該領域内の圧力
サインと速度サインとの間で規定される。比Ｒ1 および
Ｒ2 は等化演算子を形成するために結合される。過渡現
象のない時間スケールのデータは、等化演算子の存在下
で共通圧力受信器重み付け領域内のデータを共通速度受
信器重み付け領域内のスケーリングされたデータと結合
することにより規定される。
【００１２】
【実施例】本発明は、地震データにおけるノイズ、それ
に限定される訳ではないが、特に生物学的なノイズを減
衰させるための方法の教示を意図したものである。本方
法はコンピュータ処理により効率的に実施されるもので
ある。
に限定される訳ではないが、特に生物学的なノイズを減
衰させるための方法の教示を意図したものである。本方
法はコンピュータ処理により効率的に実施されるもので
ある。
【００１３】図２は共通圧力受信器（ハイドロホン）集
合体Ｈの表示であり、図３は共通速度受信器（ジオホ
ン）集合体Ｇの表示である。これらの図面において、垂
直軸は秒で表した２方向反射時間であり、水平軸はトレ
ース間が２５メートルである横方向のセンサ離間距離で
ある。本明細書を通して、用語「速度受信器」または
「速度応答形受信器」は「ジオホン」と互換的に使用さ
れ、また用語「圧力受信器」または「圧力応答形受信
器」は「ハイドロオン」と互換的に使用される。本開示
において、用語「サイン」は、時間の関数として大きさ
で表現される波形の特徴を意味し、それを区別するもの
である。。
合体Ｈの表示であり、図３は共通速度受信器（ジオホ
ン）集合体Ｇの表示である。これらの図面において、垂
直軸は秒で表した２方向反射時間であり、水平軸はトレ
ース間が２５メートルである横方向のセンサ離間距離で
ある。本明細書を通して、用語「速度受信器」または
「速度応答形受信器」は「ジオホン」と互換的に使用さ
れ、また用語「圧力受信器」または「圧力応答形受信
器」は「ハイドロオン」と互換的に使用される。本開示
において、用語「サイン」は、時間の関数として大きさ
で表現される波形の特徴を意味し、それを区別するもの
である。。
【００１４】前述した米国特許出願第０８／１０１，９
４９号により教示されている反復総和／自己相関技術を
使用して、スケール係数ｋが決定される。スケール係数
はハイドロホン信号振幅に対する平均的なジオホン信号
振幅をスケーリングするためのものである。他のノイズ
除去技術がない場合には、集合体の規格化された総和は
Σ（τ）＝１／２（Ｈ＋Ｋ＊Ｇ）である。
４９号により教示されている反復総和／自己相関技術を
使用して、スケール係数ｋが決定される。スケール係数
はハイドロホン信号振幅に対する平均的なジオホン信号
振幅をスケーリングするためのものである。他のノイズ
除去技術がない場合には、集合体の規格化された総和は
Σ（τ）＝１／２（Ｈ＋Ｋ＊Ｇ）である。
【００１５】スケール係数ｋを決定するための方法は、
次のようにまとめることができる。ハイドロホン圧力サ
インは、ジオホン速度サインの最も良い見積りを生成す
るためにフィルタ処理される。見積られたジオホン速度
サインは、ほぼ純粋なノイズサインを生成するために、
埋め込まれたノイズが加わった実際のジオホン速度サイ
ンから引かれる。次いで、ノイズサインは、純粋なジオ
ホン速度サインを得るために、実際のジオホン速度サイ
ンから引かれる。スケール係数は、反復的に純粋なジオ
ホン速度サインに適用され、ハイドロホン圧力サインに
加えられ且つ自己相関化される。スケール係数の測度
は、自己相関トレースのメインローブおよび第１サイド
ローブの考察から適切に決定される。反復は、適正測度
が最大になるまで進められる。
次のようにまとめることができる。ハイドロホン圧力サ
インは、ジオホン速度サインの最も良い見積りを生成す
るためにフィルタ処理される。見積られたジオホン速度
サインは、ほぼ純粋なノイズサインを生成するために、
埋め込まれたノイズが加わった実際のジオホン速度サイ
ンから引かれる。次いで、ノイズサインは、純粋なジオ
ホン速度サインを得るために、実際のジオホン速度サイ
ンから引かれる。スケール係数は、反復的に純粋なジオ
ホン速度サインに適用され、ハイドロホン圧力サインに
加えられ且つ自己相関化される。スケール係数の測度
は、自己相関トレースのメインローブおよび第１サイド
ローブの考察から適切に決定される。反復は、適正測度
が最大になるまで進められる。
【００１６】図２のハイドロホントレースの固定された
解析ウインド１４内では、検出された圧力サインが、第
１の共通ハイドロオン集合体を形成するために、一緒に
グループ化される。図３の速度受信器トレースの固定さ
れた解析ウインド１６内では、ジオホンサインが、第２
の共通ジオホン集合体を形成するために、一緒にグルー
プ化される。固定された解析ウインド内のハイドロホン
信号に対するジオホン信号の測定された平均トレース振
幅が、比Ｒ1 を規定するために使用される。
解析ウインド１４内では、検出された圧力サインが、第
１の共通ハイドロオン集合体を形成するために、一緒に
グループ化される。図３の速度受信器トレースの固定さ
れた解析ウインド１６内では、ジオホンサインが、第２
の共通ジオホン集合体を形成するために、一緒にグルー
プ化される。固定された解析ウインド内のハイドロホン
信号に対するジオホン信号の測定された平均トレース振
幅が、比Ｒ1 を規定するために使用される。
【００１７】固定された解析ウインドは、両方の集合体
に共通であるが全体として集合体の信号レベルを表す記
録の部分から好適に選択された横方向および縦方向領域
にわたって設定される。解析ウインドの寸法は任意であ
る。図２および図３の場合には、解析ウインドは約２．
０秒から約３．０秒まで垂直に延びており、また各集合
体について横方向に約５０トレースを取り囲んでいる。
に共通であるが全体として集合体の信号レベルを表す記
録の部分から好適に選択された横方向および縦方向領域
にわたって設定される。解析ウインドの寸法は任意であ
る。図２および図３の場合には、解析ウインドは約２．
０秒から約３．０秒まで垂直に延びており、また各集合
体について横方向に約５０トレースを取り囲んでいる。
【００１８】次のステップでは、２つの集合体は、図２
および図３の１８および２０のような重み付け領域ウイ
ンドに仕切られ、これらは両集合体に共通である。重み
付け領域ウインドは寸法が制限されており、１〜５トレ
ースの幅で且つ０．２５秒、０．５０秒あるいは１．５
秒のような任意の長さである。寸法は、所望のノイズ除
去分解能に依存する。各重み付け領域ウインドの内部に
おいて、各集合体についての重み付け領域ウインド内の
圧力信号Ｈi および速度信号Ｇi の平均振幅から、ハイ
ドロホン対ジオホン振幅比Ｒ2iが計算される（ｉは以下
に説明する移動走査インデックスである）。
および図３の１８および２０のような重み付け領域ウイ
ンドに仕切られ、これらは両集合体に共通である。重み
付け領域ウインドは寸法が制限されており、１〜５トレ
ースの幅で且つ０．２５秒、０．５０秒あるいは１．５
秒のような任意の長さである。寸法は、所望のノイズ除
去分解能に依存する。各重み付け領域ウインドの内部に
おいて、各集合体についての重み付け領域ウインド内の
圧力信号Ｈi および速度信号Ｇi の平均振幅から、ハイ
ドロホン対ジオホン振幅比Ｒ2iが計算される（ｉは以下
に説明する移動走査インデックスである）。
【００１９】動作においては、重み付け領域ウインド
は、各集合体が完全に走査されるまで、空間的（横方
向）および時間的（垂直方向）にそれぞれの集合体を通
って進んでいく。よって、Ｒ1 は単一の値であるが、各
走査ｉに対して計算された新しい値である、Ｒ2iは値の
アレイとなる。ここで、ｉ＝１、…、ｎであり、ｎは走
査の任意の番号である。よって、第１の走査について
は、ハイドロオン対ジオホン比はＲ2,1 となり、またそ
の走査についての等化演算子はｍ1 となる。
は、各集合体が完全に走査されるまで、空間的（横方
向）および時間的（垂直方向）にそれぞれの集合体を通
って進んでいく。よって、Ｒ1 は単一の値であるが、各
走査ｉに対して計算された新しい値である、Ｒ2iは値の
アレイとなる。ここで、ｉ＝１、…、ｎであり、ｎは走
査の任意の番号である。よって、第１の走査について
は、ハイドロオン対ジオホン比はＲ2,1 となり、またそ
の走査についての等化演算子はｍ1 となる。
【００２０】各走査についての等化演算子ｍi は ｍi ＝Ｒ1 ＊Ｒ2i （１） から計算される。
【００２１】ノイズのない時間スケールのデータｆ
（τ）は、等化演算子の存在下で、共通ハイドロホン受
信器重み付け領域内のサインを、共通速度受信器重み付
け領域内のスケーリングされたサインと次のように結合
することで得られる。 ｆ（τi）＝（Ｈi ＋ｋ＊Ｇi ＊ｍi ）÷（１＋ｍi ） （２）
（τ）は、等化演算子の存在下で、共通ハイドロホン受
信器重み付け領域内のサインを、共通速度受信器重み付
け領域内のスケーリングされたサインと次のように結合
することで得られる。 ｆ（τi）＝（Ｈi ＋ｋ＊Ｇi ＊ｍi ）÷（１＋ｍi ） （２）
【００２２】図４は、図２および図３からのデータの総
和に対する（２）の適用の結果得られたものである。図
２におけるノイズバースト２２、および図３における２
４と２６のようなノイズの多いトレースが、図４におい
ては十分に減じられている。
和に対する（２）の適用の結果得られたものである。図
２におけるノイズバースト２２、および図３における２
４と２６のようなノイズの多いトレースが、図４におい
ては十分に減じられている。
【００２３】用途によっては、等化演算子は、本方法を
信号増幅における迅速な変化に応答させるべく、ｍi 2
のように累乗されてもよい。
信号増幅における迅速な変化に応答させるべく、ｍi 2
のように累乗されてもよい。
【００２４】動作の最も良いモードは、次のように仮定
された数値例により示され得る。第１の例では、前述し
た米国特許出願番号第０８／１０１，９４９号の方法に
より決定されたスケール係数ｋが１．５であると仮定す
る。固定された解析ウインドに対しては、Ｒ1 ＝２．０
であった。選択された重み付け領域ウインド（ｉ＝ｘ）
においては、Ｒ2xは０．５であった。これらの量を上記
の（２）式に代入すると、ｆ（τx）は０．５（Ｈ＋
１．５Ｇx ）となり、これは規格化された、ノイズのな
い総和Σ（τ）と同一である。よって、走査ｉ＝ｘに関
しては、ハオドロホン集合体またはジオホン集合体のい
ずれもｘ番目の走査領域内にノイズバーストは発見され
なかった。
された数値例により示され得る。第１の例では、前述し
た米国特許出願番号第０８／１０１，９４９号の方法に
より決定されたスケール係数ｋが１．５であると仮定す
る。固定された解析ウインドに対しては、Ｒ1 ＝２．０
であった。選択された重み付け領域ウインド（ｉ＝ｘ）
においては、Ｒ2xは０．５であった。これらの量を上記
の（２）式に代入すると、ｆ（τx）は０．５（Ｈ＋
１．５Ｇx ）となり、これは規格化された、ノイズのな
い総和Σ（τ）と同一である。よって、走査ｉ＝ｘに関
しては、ハオドロホン集合体またはジオホン集合体のい
ずれもｘ番目の走査領域内にノイズバーストは発見され
なかった。
【００２５】第２の例としての走査ｉ＝ｙに対しては、
上記のＲ1 ＝２．０およびＲ2y＝５．０とした。つま
り、ハイドロホン集合体は、走査ｘに対するものの１０
倍の振幅を有するノイズ過渡現象と遭遇した。この場合
には、次のようになる。 ｆ（τy ）＝（Ｈy ／１１）＋１．５＊（１０／１１）
Ｇy ｙ番目の走査に関するハイドロホン重み付け領域におけ
るノイズ過渡現象は、ジオホントレースが寄与する総和
に対してその約０．１倍だけしか寄与しない。
上記のＲ1 ＝２．０およびＲ2y＝５．０とした。つま
り、ハイドロホン集合体は、走査ｘに対するものの１０
倍の振幅を有するノイズ過渡現象と遭遇した。この場合
には、次のようになる。 ｆ（τy ）＝（Ｈy ／１１）＋１．５＊（１０／１１）
Ｇy ｙ番目の走査に関するハイドロホン重み付け領域におけ
るノイズ過渡現象は、ジオホントレースが寄与する総和
に対してその約０．１倍だけしか寄与しない。
【００２６】最後の例では、走査ｉ＝ｚに対して、Ｒ1
は２．０であるが、Ｒ2Z＝０．０５である。この重み付
け領域におけるジオホントレースで受信されたノイズ過
渡現象の振幅は通常の１０倍大きい。走査ｉ＝ｚに対し
ては次のようになる。 ｆ（τz ）＝（Ｈz ／１．１）＋１．５＊（０．１／
１．１）Ｇz よって、ｚ番目のジオホン重み付け領域におけるジオホ
ントレースに関する過渡現象の効果は最小限にされる。
は２．０であるが、Ｒ2Z＝０．０５である。この重み付
け領域におけるジオホントレースで受信されたノイズ過
渡現象の振幅は通常の１０倍大きい。走査ｉ＝ｚに対し
ては次のようになる。 ｆ（τz ）＝（Ｈz ／１．１）＋１．５＊（０．１／
１．１）Ｇz よって、ｚ番目のジオホン重み付け領域におけるジオホ
ントレースに関する過渡現象の効果は最小限にされる。
【００２７】上記の技術は、地震データを処理するため
に使用される通常の汎用コンピュータにより実施するこ
とができ、コンピュータは十分なメモリおよびサブルー
チンの形でＣＰＵ内にある演算プロセッサ、あるいはア
レイプロセッサのような周辺機器を含むものである。
に使用される通常の汎用コンピュータにより実施するこ
とができ、コンピュータは十分なメモリおよびサブルー
チンの形でＣＰＵ内にある演算プロセッサ、あるいはア
レイプロセッサのような周辺機器を含むものである。
【図１】生物学的なノイズの一例を示す図である。
【図２】未処理の共通ハイドロホン−受信器集合体の地
震記録を示す図である。
震記録を示す図である。
【図３】未処理の共通速度−受信器集合体の地震記録を
示す図である。
示す図である。
【図４】本発明の教示にしたがった図２と図３の記録の
処理後の地震記録を示す図である。
処理後の地震記録を示す図である。
１４ 解析ウインド １６ 解析ウインド １８ 重み付け領域ウインド ２０ 重み付け領域ウインド
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ロナルド イー．チェインバース アメリカ合衆国，77082 テキサス，ヒュ ーストン，アシュフォード パイン 12635
Claims (6)
- 【請求項１】 複数の共に設置された速度応答形および
圧力応答形の共センサにより同時に検出された地震信号
の時間スケールの記録から好ましくない異常な過渡現象
を適応的に抑制する方法であって、 ａ）ハイドロホンおよびジオホンの時間スケールのサイ
ンを供給すべく、前記圧力応答形および速度応答形の共
センサによる検出のために水中環境において反射地震波
場を発生させる工程と、 ｂ）振幅においてハイドロホンサインと一致させるべ
く、ジオホンサインをスケーリングするためのスケール
係数ｋを規定する工程と、 ｃ）共通ハイドロホン集合体を形成すべく前記圧力応答
形のセンサからの検出したサインをグループ化し、且
つ、共通ジオホン集合体を形成すべく前記速度応答形の
センサからの検出したサインをグループ化する工程と、 ｄ）両集合体に共通の固定された解析ウインドを規定す
る工程と、 ｅ）前記固定された解析ウインド内の共通ジオホン集合
体と共通ハイドロホン集合体との間の比Ｒ1 を規定する
工程と、 ｆ）共通ハイドロホン集合体および共通ジオホン集合体
を複数の共通重み付け領域に仕切る工程と、 ｇ）第１の重み付け領域内において、ハイドロホンサイ
ンとジオホンサインとの間の比Ｒ2,1 を規定する工程
と、 ｈ）Ｒ1 とＲ2,1 との結合から等化演算子ｍ1 を形成す
る工程と、 ｉ）過渡現象のない時間スケールのデータｆ（τ）を得
るべく、前記等化演算子の存在下で共通ハイドロホン重
み付け領域内のハイドロホンサインを共通ジオホン重み
付け領域内のスケーリングされたジオホンサインと結合
する工程と、を具備する方法。 - 【請求項２】 重み付け領域の残りについて工程ｇ）〜
工程ｉ）を繰り返す工程を更に具備する請求項１記載の
方法。 - 【請求項３】 前記等化演算子がｍi ＝Ｒ1 ＊Ｒ2iによ
り計算される請求項２記載の方法。 - 【請求項４】 前記等化演算子がｍi n ＝（Ｒ1 ＊
Ｒ2i）n により計算され、ｎが１以上の整数である請求
項２記載の方法。 - 【請求項５】 工程ｉ）が式ｆ（τi）＝（Ｈi ＋ｋ＊
Ｇi ＊ｍi ）／（１＋ｍi ）の式により実施され、Ｈi
は共通ハイドロホン集合体データおよびｋ＊Ｇi は重み
付け領域ウインド内のスケーリングされた共通ジオホン
集合体データを表す請求項４記載の方法。 - 【請求項６】 地震信号内の異常な過渡現象を減衰させ
るノイズ除去方法であって、 圧力応答形の地震信号と速度応答形の地震信号とを受信
し、且つ、受信したこれら地震信号を共通圧力信号集合
体Ｈと共通速度信号集合体Ｇとにそれぞれフォーマット
する工程と、 前記集合体の両方に共通な選択された空間寸法および時
間寸法を有する領域で固定された解析ウインドを設定
し、且つ、前記ウインド内において、速度信号集合体Ｇ
と圧力信号集合体Ｈとの間の平均振幅比Ｒ1 を決定する
工程と、 それぞれの集合体を、制限された横方向範囲および縦方
向範囲を有する共通重み付け領域に仕切る方法と、 各重み付け領域内で圧力信号集合体Ｈi および速度信号
集合体Ｇi の平均振幅の間の比Ｒ2iを決定する工程であ
って、ｉは走査インデックスである、ものと、 Ｒ1 とＲ2iとを結合して等化演算子ｍi を形成する工程
と、 等化演算子ｍi の存在下で圧力信号集合体Ｈi をスケー
リングされた速度信号集合体ｋ＊Ｇi と総和することに
より、各重み付け領域ウインド内で過渡現象のない時間
スケールのデータｆ（τi ）を規定する工程と、を具備
する方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US08/263056 | 1994-06-21 | ||
| US08/263,056 US5442591A (en) | 1994-06-21 | 1994-06-21 | Method for adaptively suppressing noise transients in summed co-sensor seismic recordings |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0854473A true JPH0854473A (ja) | 1996-02-27 |
| JP2857351B2 JP2857351B2 (ja) | 1999-02-17 |
Family
ID=23000194
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP15466895A Expired - Fee Related JP2857351B2 (ja) | 1994-06-21 | 1995-06-21 | 総和された共センサ地震記録における過渡ノイズを適応的に抑制する方法 |
Country Status (10)
| Country | Link |
|---|---|
| US (1) | US5442591A (ja) |
| EP (1) | EP0689064B1 (ja) |
| JP (1) | JP2857351B2 (ja) |
| AU (1) | AU681509B2 (ja) |
| CA (1) | CA2152207C (ja) |
| CO (1) | CO4340574A1 (ja) |
| DE (1) | DE69513090D1 (ja) |
| NO (1) | NO308495B1 (ja) |
| NZ (1) | NZ272370A (ja) |
| OA (1) | OA10220A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018185206A (ja) * | 2017-04-25 | 2018-11-22 | 国立大学法人東京海洋大学 | 海底地質探査システム、海底地質探査方法および海底地質探査プログラム |
| JP2021526635A (ja) * | 2018-05-30 | 2021-10-07 | サウジ アラビアン オイル カンパニー | 海底受振器の地震データのノイズ対策 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| US5737433A (en) * | 1996-01-16 | 1998-04-07 | Gardner; William A. | Sound environment control apparatus |
| WO2002097472A2 (en) | 2001-05-25 | 2002-12-05 | Exxonmobil Upstream Research Company | Multiple suppression for ocean bottom seismic data |
| US6894948B2 (en) * | 2003-01-29 | 2005-05-17 | Pgs Exploration (Uk) Limited | Method for processing dual sensor seismic data to attenuate noise |
| US7953556B2 (en) * | 2007-03-09 | 2011-05-31 | Fairfield Industries Incorporated | Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique |
| US20120253681A1 (en) * | 2011-03-31 | 2012-10-04 | Chevron U.S.A. Inc. | System and method for processing seismic data |
| CN104730580B (zh) * | 2013-12-23 | 2017-06-23 | 中国石油化工股份有限公司 | 地震资料异常振幅压制方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4344158A (en) * | 1979-04-27 | 1982-08-10 | Western Geophysical Co. Of America | Noise-suppression method |
| US4937794A (en) * | 1985-05-06 | 1990-06-26 | Western Atlas International, Inc. | Seismic noise suppression method |
| US5293352A (en) * | 1993-01-07 | 1994-03-08 | Western Atlas International, Inc. | Method for removing noise due to near surface scatterers |
| US5365492A (en) * | 1993-08-04 | 1994-11-15 | Western Atlas International, Inc. | Method for reverberation suppression |
| US5408441A (en) * | 1993-08-20 | 1995-04-18 | Halliburton Company | Method for seismic exploration in artic regions |
-
1994
- 1994-06-21 US US08/263,056 patent/US5442591A/en not_active Expired - Lifetime
-
1995
- 1995-06-15 NZ NZ272370A patent/NZ272370A/en unknown
- 1995-06-19 EP EP95304251A patent/EP0689064B1/en not_active Expired - Lifetime
- 1995-06-19 DE DE69513090T patent/DE69513090D1/de not_active Expired - Lifetime
- 1995-06-20 OA OA60677A patent/OA10220A/en unknown
- 1995-06-20 CA CA002152207A patent/CA2152207C/en not_active Expired - Fee Related
- 1995-06-20 CO CO95026986A patent/CO4340574A1/es unknown
- 1995-06-20 AU AU21785/95A patent/AU681509B2/en not_active Ceased
- 1995-06-20 NO NO952453A patent/NO308495B1/no not_active IP Right Cessation
- 1995-06-21 JP JP15466895A patent/JP2857351B2/ja not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2018185206A (ja) * | 2017-04-25 | 2018-11-22 | 国立大学法人東京海洋大学 | 海底地質探査システム、海底地質探査方法および海底地質探査プログラム |
| JP2021526635A (ja) * | 2018-05-30 | 2021-10-07 | サウジ アラビアン オイル カンパニー | 海底受振器の地震データのノイズ対策 |
Also Published As
| Publication number | Publication date |
|---|---|
| NO308495B1 (no) | 2000-09-18 |
| CA2152207C (en) | 2001-02-27 |
| AU681509B2 (en) | 1997-08-28 |
| AU2178595A (en) | 1996-01-04 |
| JP2857351B2 (ja) | 1999-02-17 |
| NO952453L (no) | 1995-12-22 |
| DE69513090D1 (de) | 1999-12-09 |
| EP0689064A3 (en) | 1998-04-22 |
| CO4340574A1 (es) | 1996-07-30 |
| US5442591A (en) | 1995-08-15 |
| NO952453D0 (no) | 1995-06-20 |
| CA2152207A1 (en) | 1995-12-22 |
| EP0689064B1 (en) | 1999-11-03 |
| OA10220A (en) | 1997-09-19 |
| NZ272370A (en) | 1997-04-24 |
| EP0689064A2 (en) | 1995-12-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5365492A (en) | Method for reverberation suppression | |
| US8775091B2 (en) | Marine seismic surveying employing interpolated multi-component streamer pressure data | |
| US5995905A (en) | Source signature determination and multiple reflection reduction | |
| CN102016643B (zh) | 用于衰减双传感器地震拖缆中的低频噪声的方法 | |
| EP0681193B1 (en) | Method of attenuating coherent noise in marine seismic data | |
| US5524100A (en) | Method for deriving water bottom reflectivity in dual sensor seismic surveys | |
| US5621700A (en) | Method for attenuation of reverberations using a pressure-velocity bottom cable | |
| US6894948B2 (en) | Method for processing dual sensor seismic data to attenuate noise | |
| EP2108980B1 (en) | Method for deghosting marine seismic streamer data with irregular receiver positions | |
| US5621699A (en) | Apparatus and method of calibrating vertical particle velocity detector and pressure detector in a sea-floor cable with in-situ passive monitoring | |
| CN101014881A (zh) | 用于压力传感器和质点运动传感器记录的地震数据中水底多次波的衰减的系统 | |
| JP2857351B2 (ja) | 総和された共センサ地震記録における過渡ノイズを適応的に抑制する方法 | |
| US6246637B1 (en) | Method and system for combining three component seismic data | |
| Frazer et al. | Seabed sediment attenuation profiles from a movable sub-bottom acoustic vertical array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19981028 |
|
| LAPS | Cancellation because of no payment of annual fees |