JP2020519889A - 全波形音波動場分離のための処理方法 - Google Patents

全波形音波動場分離のための処理方法 Download PDF

Info

Publication number
JP2020519889A
JP2020519889A JP2019562287A JP2019562287A JP2020519889A JP 2020519889 A JP2020519889 A JP 2020519889A JP 2019562287 A JP2019562287 A JP 2019562287A JP 2019562287 A JP2019562287 A JP 2019562287A JP 2020519889 A JP2020519889 A JP 2020519889A
Authority
JP
Japan
Prior art keywords
fws
wave
data
aliasing
velocity
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
Application number
JP2019562287A
Other languages
English (en)
Other versions
JP6943979B2 (ja
Inventor
スン,シュエカイ
ビー. アヤディウノ,クリス
ビー. アヤディウノ,クリス
プランチャート,カルロス
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Saudi Arabian Oil Co
Original Assignee
Saudi Arabian Oil Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Saudi Arabian Oil Co filed Critical Saudi Arabian Oil Co
Publication of JP2020519889A publication Critical patent/JP2020519889A/ja
Application granted granted Critical
Publication of JP6943979B2 publication Critical patent/JP6943979B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32Transforming one recording into another or one representation into another
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/40Transforming data representation
    • G01V2210/46Radon transform
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/57Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Image Processing (AREA)
  • Complex Calculations (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Controls And Circuits For Display Device (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

全波形音波(FWS)波動場分離の方法は、FWSデータを受信するステップと;前記受信したFWSデータにアンチエイリアシング線形ラドン変換を実行するステップと;スラントウィンドウを用いて、波成分に対応するラドン変換済みFWSデータを抽出するステップと;前記抽出されたラドン変換済みFWSデータに逆ラドン変換を実行することによって、前記波成分の信号を特定するステップと;を含む。

Description

本願は、2017年5月11日に出願された米国特許出願第15/592,979号の優先権を主張し、その内容全体は参照により本明細書に組み込まれる。
本願は、全波形音波(フルウェーブフォーム音波、FWS)データ処理に関する。
全波形音波(FWS)デバイスは、FWSデータを収集するために、坑井において使用することができ、該FWSデータは、地層の岩相特性及び流体特性についての情報を提供するために処理できる。地震及び垂直地震プロファイル(VSP)といった他の地球物理学的測定とは異なり、FWSデータは、典型的には、短いサンプリング時間間隔を用いて測定されるが、オフセット(offset)間隔が大きいために、FWSデータのエイリアシング効果の一因になる。既存のFWS処理方法は、エイリアシング効果を大目にみて、他の分野において開発された処理技術をFWSデータに単に適用しており、この結果、FWS波動場の分離結果における不確実性及び不正確さをもたらす。
本開示は、全波形音波(FWS)波動場分離のための処理方法に関連するコンピュータシステム、コンピュータプログラム製品、及びコンピュータにより実施される方法を含む、方法及びシステムについて説明する。
いくつかの実施では、FWSデータが受信される。アンチエイリアシング線形ラドン変換が、受信したFWSデータに対して実行される。波成分に対応するラドン変換済みFWSデータが、スラントウィンドウを用いて抽出される。該波成分の信号は、抽出されたラドン変換済みFWSデータに対して逆ラドン変換を実行することによって特定される。
前述の実施は、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータにより実施されるシステムを用いて実施可能である。非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータ読取可能な命令を格納する。コンピュータにより実施されるシステムは、非一時的なコンピュータ読取可能媒体に格納された命令/コンピュータにより実施される方法を実行するように構成されたハードウェアプロセッサと相互に利用可能に結合されたコンピュータメモリを含む。
本願において説明された主題は、いくつかの実施においては、アンチエイリアシング線形ラドン変換(AALRT)及び適応的傾斜抽出(ASE)を使用することによって、エイリアシングアーチファクトのないFWS波成分を効率的に分離できる。AALRTは、FWSデータの効率的な局所的な補間を可能にするために、指定された速度値によって特定されたアンチエイリアシングオフセット間隔を使用する。望まれない波成分の一部又は全てを含むことなく対象波の抽出を可能にするために、ASEは、FWSオフセット範囲に基づきスラントウィンドウを適応的に調整できる。手動の抽出方法と比較して、ASEは時間を節約でき、人間の先入観を低減できる。他の利点は、当業者には明らかであろう。
本願に記載された主題における一又は複数の実施の詳細は、添付図面及び引き続く明細書に説明される。本主題の他の特徴、態様及び利点は、明細書、図面及び特許請求の範囲から明らかになる。
図1は、本願のいくつかの実施に係る、全波形音波(FWS)波動場分離のための例示的な方法のフローチャートである。
図2A〜図2Cは、本願のいくつかの実施に係る、圧縮波のためのFWS波動場分離を示すものである。
図3A〜図3Cは、本願のいくつかの実施に係る、合成FWSギャザーのFWS波動場分離を示すものである。 図4A〜図4Cは、本願のいくつかの実施に係る、合成FWSギャザーのFWS波動場分離を示すものである。
図5A〜図5Dは、本願のいくつかの実施に係る、合成されていない実際のFWSギャザーのためのFWS波動場分離を示すものである。
図6は、本願のいくつかの実施に係る、本願において記述されたアルゴリズム、記述された方法、記述された機能、記述されたプロセス、記述されたフロー、及び記述された手順に関連付けられた計算機能を提供するために使用された例示的なコンピュータシステムを示すブロック図である。
種々の図面における類似の参照番号及び記号は、同様の要素を指す。
以下の詳細な説明は、全波形音波(FWS)波面分離の処理方法について説明し、1又は複数の特定の実施の観点において開示された主題を、当業者が作成し及び使用できるように提示される。開示された実施の様々な修正、変更、及び並べ替えを行うことができることは当業者には容易に明らかであり、また定義された一般原理は、本開示の範囲から逸脱することなく、他の実施及び応用(アプリケーション)に適用されるようにしてもよい。これ故に、本願は、説示または例示された実施に限定されることを意図するものではなく、開示された原理及び特徴と整合する最も広い範囲が与えられるべきである。
FWSデバイスは、FWSデータを収集するために、坑井において使用でき、該FWSデータは、地層の岩相特性及び流体特性についての情報を提供するために処理することができる。FWSデータは、線形又は準線形の屈折波、例えば、圧縮波、せん断波、ストンリー波の情報を含むことができ、該屈折波は、様々な波速度によって特徴付けられる。線形ラドン(Radon)変換(又はスラントスタック(傾斜重合))が、FWS波動場分離のために様々な波成分を抽出するために使用できる。ラドン変換は地震及びVSPの処理に既に使用されているが、FWSデータ処理に対する使用は制限されてきた。地震を記録することと比較すると、空間方向(つまり、オフセット方向)におけるFWSサンプルは、時間(time)方向のサンプルよりかなり粗(sparse)である。粗な空間サンプルに起因して、地震処理のために開発されたラドン変換の方法論を直接的にFWSデータに対して適用することを妨げるエイリアシング(aliasing)効果が、存在する。
線形ラドン変換は、t=px+τの関係に基づき、データをoffset(x)−time(t)ドメイン(領域)からslowness(p)−intercept(τ)ラドンパラメータドメインに変換する。ここで、オフセットxは、震源(ソース)送信器と受振器との距離を表すことができ、τは、ゼロオフセット(つまりx=0)における移動時間を示し、slowness(スローネス)のpは、p=1/vによって、速度vに関連付けることができる。FWSデータの場合、密な時間サンプル及び粗な空間サンプルに起因して、通常、エイリアシング(疑信号)は、slowness(p)方向に生じる。ナイキストの法則に基づき、ラドンパラメータドメインにおいてスローネスサンプリングを制御するアンチエイリアシング(anti−aliasing)条件は、次のように表現できる。
Figure 2020519889
ここで、fmaxはFWSデータの最大有効周波数であり、△xは、オフセット間隔であり、|x|max及び|x|minは、絶対値における最大及び最小のオフセット距離である。いくつかのFWSデバイスでは、FWS測定におけるサンプリング時間間隔(つまり、FWSデータを記録するときのサンプリング時間間隔)及びオフセット間隔は、それぞれ、2e−5秒(s)及び0.1524メートル(m)である。ナイキストの法則は、最大有効周波数fmax=1/(2×2e−5)=25000Hzであることを示す。式(1)によると、スローネスp<1/(25000×0.1524)=0.0003s/mであり、速度v>1/0.0003=3810m/sである。換言すれば、上記のFWS測定の設定の場合におけるエイリアシングアーチファクトの無い速度範囲は、3810m/sよりも大きく、ほとんどの地層のせん断速度及びストンリー波の速度を超えている。いくつかの他の場合には、FWSデータは、さらにいっそう短いサンプリング時間間隔及びより大きなオフセット間隔、例えば1e−5s及び0.3048m、を使用可能なFWSデバイスによって入手でき、波速度は、エイリアシング(疑信号)を回避するために15240m/sより大きくすべきである。換言すると、FWS波の多くの成分は、エイリアシングアーチファクトの下では区別できないことがあり得る。
高レベルにおいては、説明されたアプローチは、FWS波動場の特徴を利用することによってFWS波動場を分離するための処理フローを説明する。説明したアプローチは、波動場を分離するために、アンチエイリアシング線形ラドン変換(Anti−Aliasing Linear Radon Transform、AALRT)及び適応的傾斜抽出(Adaptive Slant Extraction、ASE)を使用する。AALRTは、アンチエイリアシングオフセット間隔に基づきFWSデータを補間する。AALRTは、実際の地層速度範囲に基づきFWSエイリアシング効果に柔軟に取り組むことができ、この結果として、波面分離のための高品質なラドン画像を発生できる。ASEは、他の波成分情報を含めずに実際のFWS構成に適応させることができるスラントウィンドウ(slant window、傾斜されたウィンドウ)を用いて各波成分を抽出できる。
いくつかの実施では、FWSデータが受信される。受信したFWSデータは、共通ショットギャザー(gather、波形出力の集合)であり得る。アンチエイリアシング線形ラドン変換を、受信したFWSデータに対して実行できる。波成分に対応するラドン変換済みFWSデータは、スラントウィンドウを用いて抽出できる。波成分の信号は、抽出されたラドン変換済みFWSデータに対して逆ラドン変換を実行することによって特定できる。波成分は、圧縮波、せん断波、又はストンリー波のいずれかであってよい。いくつかの実施においては、アンチエイリアシング線形ラドン変換を実行することは、アンチエイリアシングオフセット間隔に基づき一組のアンチエイリアシングオフセット値及び一組のスローネス値を特定すること、受信したFWSデータ及び一組のアンチエイリアシングオフセット値に基づき補間されたFWSデータを生成すること、及び一組のアンチエイリアシングオフセット値及び一組のスローネス値に基づき、補間されたFWSデータに対してラドン変換を実行すること、を含むことができる。アンチエイリアシングオフセット間隔は、△x=0.5×△t×Vminであり、ここで、△tは、受信したFWSデータのサンプリング時間間隔であり、Vminは、波成分の速度の下限である。いくつかの実施においては、速度解析は、波成分の中心速度、速度範囲、中心時間、及び時間範囲を特定するために実行できる。スラントウィンドウは、波成分の中心速度及び中心時間に基づき特定された点を中心にして置くことができ、またスラントウィンドウは、波成分の速度範囲及び時間範囲に基づき特定されたサイズを持つことができる。いくつかの実施においては、抽出されたラドン変換済みFWSデータに対して逆ラドン変換を実行することは、一組の実際のオフセット値に基づき逆ラドン変換を実行することを含むことができる。
図1は、いくつかの実施に係る、FWS波動場分離のための例示的な方法100のフローチャートである。説明を明確にするために、以下の説明は、一般に、本願における他の図面の文脈にて方法100を説明する。例えば、方法100は、図6に記載されたコンピュータシステム、若しくは任意の適切なシステム、環境、ソフトウェア、及びハードウェア、又はシステム、環境、ソフトウェア、ハードウェアの適切な組み合わせによって実行することができる。いくつかの実施においては、方法100の様々なステップは、並行して、組み合わせて、ループで、又は任意の順序で実行することができる。
高レベルにおいては、方法100は主に、(1)情報チェック及び前処理、(2)アンチエイリアシング線形ラドン変換、(3)適応的傾斜抽出、及び(4)逆ラドン変換、の各ステップを含む。ブロック102においては、方法100は、対象領域のFWSデータ(又はトレース)を受信し及び前処理することを含む。例えば、対象領域は、1又は複数の地下の層を含むことができる。震源及び受振器を含むFWSデバイスは、対象領域内における掘削孔内へ下ろすことができる。FWSの震源は、音波又はショット(shot)を掘削孔及び隣接の岩層に送ることができ、受振器は、屈折波又は反射波を測定及び記録することができる。受振器は、サンプリング時間間隔に合わせて受信波形をサンプリングし、サンプリングされたデータを記録することができる。場合によっては、震源は、受振器が反射波又は屈折波を記録するために、異なる場所で音波を発生することができる。例えば、FWSデバイスは、掘削孔内の第1の深さまで下ろすことができ、FWS震源は、受振器が屈折波又は反射波を記録するための第1の波形又はショットを送出することができる。FWSデバイスは、次いで、掘削孔内の第2の深さに移されて、第2の波形又はショットを送出することができる。単一の震源波伝送又はショットに対応する一の受振器で記録されたデータは、トレースと称する。いくつかの実施においては、ブロック102において受信したFWSデータは、多数の受振器において記録された単一のショットからのトレースを含む共通ショットギャザーである。
受信したFWSデータは、ブロック104及び106の処理ステップにおいて使用できる情報を生成するために、前処理を行うことができる。前処理は、インタラクティブな速度解析のステップ、線形ムーブアウト補正のステップ、及び他の前処理のステップを含むことができる。前処理はまた、受信したFWSデータの時間範囲、オフセット範囲、及び周波数範囲を特定することができる。例えば、周波数範囲は、受信したFWSトレースのフーリエ変換を実行すると共にフーリエドメインにおいて最大信号周波数を識別することによって特定される。インタラクティブな速度解析は、波成分毎に時間及び速度の情報を特定できる。速度解析は、波成分毎の複数のタイプ情報、例えば中央時間及び中央速度、並びに時間と速度との両方の上限及び下限を含む速度ファイルを生成できる。中心速度及び中心時間は、沈んだ(dipped)波動イベントを平坦化する線形ムーブアウト補正だけでなく、ラドンドメインにおける波成分の中心位置(τ、p)を特定する際にも使用できる。速度及び時間のための上限及び下限(つまり、速度範囲及び時間範囲)、加えて(τ、p)は、ブロック106においてASEのスラントウィンドウを特定するために使用できる。
ブロック104においては、AALRTが、受信したFWSデータに対して実行される。AALRTにおいては、従来採用されていたトレース補間の代わりに、ラドンパラメータドメインにおいて局所化(ローカライズ)された補間のスキームが、利用可能なデータ情報を圧縮するために使用される。式(1)におけるアンチエイリアシング条件に基づき、アンチエイリアシングオフセット間隔は。△x=0.5*△t*Vminとして定義でき、ここで、△tは、FWSデータを記録するために受信波をサンプリングする際のサンプリング時間間隔であり、また、Vminは、波動速度の下限である。例えば、ブロック102においては、速度の下限が、圧縮波、せん断波、及びストンリー波の各々について特定することができ、Vminは、3つの下限値のうちの最小値とすることができる。ここで、△xは、Vminより大きい速度を持つ波成分のための理論的なオフセット間隔であり、該波成分は、有効に測定され且つ有効に処理されると共に効果的にエイリアシングアーチファクトから免れていることに留意すべきである。アンチエイリアシングオフセットの広がりxは、最初に測定されたオフセットから最後に測定されたオフセットまでの範囲にあり、またΔxの間隔を空けて生成できる。例えば、実際のオフセットの広がり(つまり、FWSデータ獲得における震源と複数の受振器との間の実距離)が、オフセット間隔0.15mを持つベクトルx=[0.15、0.3、0.45、0.6]mとして表現される。アンチエイリアシングオフセット間隔△x=0.05mを用いると、アンチエイリアシングオフセットの広がりは、ベクトルx=[0.15、0.2、0.25、0.3、0.35、0.4、0.45、0.5、0.55、0.6]mとして表現できる。
受信したFWSデータを補間するために、FWS共通ショットギャザーを形成し、周波数ドメイン(F−Xドメインとして知られている)に変換することができる。周波数ドメインにおけるFWSギャザーは、下式に基づき、現実のオフセットの広がりxにおける勾配g及び切片iといった線形表現を特定するために使用できる。
Figure 2020519889
ここで、Diag(x)は、その対角要素が実際のオフセットの広がりxである対角行列であり、Diag(1)は、Diag(x)と同じサイズを有すると共にその対角要素が1である対角行列であり、g及びiは各々ベクトルxと同数の要素を持つ列ベクトルであり、またd(f)は、周波数fを中心としたフィルタリング処理されたFWSデータである。例えば、d(f)は、時間方向に沿って、トレース毎に共通ショットギャザーに対して離散フーリエ変換を実行することによって発生されることができ、またd(f)は、全トレースの部分集合又は全トレースからであって特定周波数fでのフーリエ変換済みデータを含むことができる。式(2)は、劣決定系逆問題(under−determined inversion)であり得、多くの正則化手法が、gとiの滑らかな推定値を生成するために使用できる。いくつかの実施では、式(2)は、離散フーリエ変換からの周波数ポイント毎に解くことができる。
g及びiが特定された後に、アンチエイリアシングオフセットの広がりx上における補間されたデータd(f)は、局所的な補間を実行することによって生成できる。いくつかの実施では、局所的な補間は、xをxに置き換えると共にアンチエイリアシングオフセットの広がりxのためにg及びiを補間することによって、式(2)を順方向に適用することになる。例えば、補間されたデータd(f)は、以下のように特定できる。
Figure 2020519889
ここで、Diag(x)は、その対角要素がアンチエイリアシングオフセットの広がりxである対角行列であり、Diag(1)は、Diag(x)と同サイズを有すると共にその対角要素が1である対角行列であり、g及びiは各々、ベクトルxと同数の要素を有する列ベクトルである。
ベクトルg及びiは、近接の原理を用いて、g及びiに基づき補間できる。換言すれば、アンチエイリアシングオフセットにおける勾配及び切片は、該アンチエイリアシングオフセットに最も近い実際のオフセットにおける勾配及び切片によって特定できる。例えば、実際のオフセットの広がりx=[0.15、0.3、0.45]mの場合、g及びiは、それぞれ、g=[g、g、g及びi=[i、i、iとして表すことができる。アンチエイリアシングオフセットの広がりが、x=[0.15、0.2、0.25、0.3、0.35、0.4、0.45]mの場合、アンチエイリアシングオフセット0.2mが実際のオフセット0.15mに最も近いので、実際のオフセット0.15mにおける勾配及び切片であるg及びiは、アンチエイリアシングオフセット0.2mにおける勾配及び切片として使用できる。つまり、g=[g、g、g、g、g、g、g]及びi=[i、i、i、i、i、i、i]となる。
アンチエイリアシングオフセットの広がりxに対応するラドンパラメータドメインにおけるslownessサンプリング値pは、式(1)のアンチエイリアシング条件にxを代入することによって特定できる。換言すると、slownessサンプリングベクトルpにおけるi番目の要素は、以下の通り特定できる。
Figure 2020519889
ここで、xa(i)は、アンチエイリアシングオフセットベクトルxのi番目の要素であり、fmaxは、ブロック102において特定されたFWSデータの最大信号周波数である。
アンチエイリアシングオフセットの広がりxにおける補間済みFWSデータのラドン変換は、次の反転形式として定式化できる。
Figure 2020519889
a(n)は、アンチエイリアシングオフセットベクトルxのn番目の要素であり
、pa(j)は、slownessサンプリングベクトルpのj番目の要素であり、Lj、nは、行列Lの(j、n)番目の要素である。式(3)への正則化された解は、周波数fに関するラドン変換済みデータを以下のように生成できる。
Figure 2020519889
ここで、Lは、Lのエルミート転置を示し、μは、ノイズの影響を抑制する事前白色化(pre−whitening)係数である。いくつかの実施では、μは、LLの最大絶対値のパーセンテージ(例えば、0.1%〜1%の間のパーセンテージ)として特定できる。いくつかの実施では、式(4)は、周波数毎のラドン変換済みデータを生成するために使用され、またラドン画像は、逆フーリエ変換を用いて、全周波数の部分集合又は全てのラドン変換済みデータを、時間ドメインへ変換することによって生成される。換言すれば、式(4)は、補間済みFWSデータのラドン画像を生成するために使用することができる。式(3)は、他のラドン変換タイプと同様な形式及び解を有するものの、式(3)は、次の2つの側面に関して他のものと異なる:(1)xは、所望の速度範囲を十分にモデル化できるアンチエイリアシングオフセットの広がりとして計算されまた採用されること;(2)xにわたってデータ情報を補償するために、局所化された補間は、dを効率的に取得でき、また周波数及びオフセット上の局所化された詳細は、該プロセスにおいて大幅に保護されること。
ブロック106では、スラントウィンドウは、ブロック104において生成されたラドン画像から波成分を抽出するために使用することができる。クロスハッチされた(cross−hatched)エネルギクラスターが、有限の空間オフセットに起因してラドン画像に典型的に観察され、これ故に、通常用いられる長方形ウィンドウは、クロスハッチされた対象信号をラドン画像から抽出するには非効率的である。FWS波動場分離の場合には、長方形の抽出は、対象波成分上において(図3Bに示されるように)不完全な又は混在した抽出を引き起こすことがあり、したがって分離品質を低下させる。クロスハッチされた対象信号が、所定のFWS構成に従って効率的にキャプチャするために、ASEは、FWS波動場分離のためにスラントウィンドウを使用できる。
ラドンパラメータドメインにおいて(τ、p)を中心とする対象波を抽出するために、ウィンドウ長の半分nτが特定できる。例えば、nτは、スラントウィンドウがラドンドメインにおいて対象波を囲むことができるように特定できる。いくつかの実施では、ブロック102において実行された速度及び時間の解析が、スラントウィンドウを特定するために使用できる。τ方向におけるスラントウィンドウの上下の境界は、次の式で特定できる。
Figure 2020519889
ここで、absは、絶対値を取る操作を示し、min及びmaxは、それぞれ、最小化操作と最大化操作を示し、min(abs(x))*p及びmax(abs(x))*pにおける演算子*は、要素毎の乗算を示し、また(τ、p)は、ブロック102において特定される。式(5)の結果としての上限及び下限は、ベクトルpと同サイズのベクトルであり、また上限及び下限の境界ベクトルにおける要素は、対応するslownessサンプリング値pにおける境界値を表す。式(5)では、アンチエイリアシングオフセット変数xによれば、クロスハッチされた対象信号の適応的な記述が可能になり、また最小化操作及び最大化操作は、抽出のためにウィンドウエリアを拡大できる。いくつかの実施では、対象波の速度範囲又は時間範囲が利用可能な場合、抽出ウィンドウをさらに絞り込むことができる。例えば、ブロック102において特定された対象波の速度範囲に基づき、式(5)における完全なベクトルp及びxを使用することに代えて、特定された速度範囲に対応するp及びxの部分集合が、スラントウィンドウを導き出すために使用できる。いくつかの実施では、半ウィンドウ長nτが、ブロック102において特定された対象波の時間範囲に基づき特定できる。
ブロック108では、対象波成分の信号は、ラドンパラメータドメイン内の抽出された成分に逆ラドン変換を実行することによって回復する(recover)ことができる。例えば、逆ラドン変換は、式(3)を用いて、実際のオフセットの広がりxに対して実行できる。いくつかの実施では、逆ラドン変換の結果は、複素数であり得、該複素数の実数部を用いると共に虚数部を削除することで、波成分を回復することができる。
図2A〜図2Cは、いくつかの実施に係る、圧縮波のためのFWS波動場分離を示す。図2Aは、0.3048mの間隔を有する8つの受振器上に記録されると共に5000m/秒の波速度を持つ圧縮波をシミュレートする合成FWSギャザー200aを例示する。データ記録のためのサンプリング時間間隔は1e−5秒であり、合成FWSギャザーの主要な周波数は10kHzである。図2Aの横軸及び縦軸は、それぞれ、オフセット指標及び時間を表す。図2Bは、合成FWSギャザー200aに対して従来のラドン変換を直接に適用した後のラドン画像200bを例示する。図2Bは、従来のラドン変換が散乱エネルギー写像(map)になる可能性があることを表しており、波成分が区別できないことを示す。図2Cは、合成FWSギャザー200aにAALRTを適用した後のラドン画像200cを例示する。ドット202は、圧縮波成分を表し、この圧縮波成分は2e−3s/mのslowness、つまり、波速度5000m/sの逆数を有する。また、図2Cは、アンチエイリアシングオフセットxを採用することと局所的な補間を実行することとを含むAALRTを適用した後、ラドン画像200cは、ラドン画像200bと比較すると、slowness(p)−intercept(τ)ドメインにおいてよりフォーカスされるようになることを示す。例えば、従来のラドン変換は、−2e−5s/m〜2e−5s/mの範囲のslownessサンプリング値のラドン画像200bを生成する一方で、AALRTは、−1e−3s/m〜1e−3s/mの範囲のアンチエイリアシングslownessサンプリング値のラドン画像200cを生成する。換言すれば、1000m/sを超える速度の波成分は、AALRTを用いてラドンパラメータドメインにおいて、有効にモデル化されまた特徴付けられ得る。
図3A〜図4Cは、いくつかの実施に係る、合成FWSギャザーのFWS波動場分離を例示する。合成ギャザーは、0.3048mの間隔の8つの受振器に記録されると共に異なる主要な周波数を持つ圧縮波、せん断波、ストンリー波をシミュレートする。図3A〜図3Cは、ラドンドメインにおける波動場分離を示し、図4A〜図4Cは、図3A〜図3Cに対応するoffset(x)−time(t)ドメインにおける波動場を示す。図4A〜図4Cは、それぞれ、図3A〜図3Cの逆ラドン変換を実行することによって生成することができる。図3Aは、ラドンドメインにおける圧縮波302及びせん断波304を含むトータル波動場300aを示し、該ラドンドメインでは、圧縮波の線形ムーブアウト補正が既に実行されており、圧縮波302が分離されるべき対象波である。図4Aは、オフセット時間ドメインにおける合成FWSギャザーのトータル波動場400aを示しており、トータル波動場400aは、圧縮波402及びせん断波404を含む。図3Bは、小さな長方形ウィンドウには対象波成分302の完全な特徴を取り込む能力がないことを示しているが、大きなウィンドウ306は、せん断波304といった他の波成分に容易に切り込むことができ、図4Bに示されるように、混ざりものがある(inpure)分離の結果をもたらし、対象の圧縮波402に加えてせん断波404が存在する。図3Cは、ASEのスラントウィンドウ308が、せん断波304に切り込むことなく対象成分302を抽出できることを例示する。図4Cは、ASEが、せん断波404ではなく対象の圧縮波402を含む混じりけのない(クリーンな)分離の結果を生成することを例示する。
図5A〜図5Dは、いくつかの実施に係る、実際のFWSギャザーのためのFWS波動場分離を示す。図5Aは、実際のFWSショットギャザー500aを示しており、縦軸及び横軸は、FWSギャザーの時間サンプル指標及びオフセット指標を表す。FWSギャザー500aは、浅い圧縮波成分508を除いて、せん断波及びストンリー波がエイリアシング効果に起因して区別できないことを示す。またノイズは、特に深いストンリー波の波動場分離のためには好ましくない要因である。図5B〜図5Dは、FWSギャザー500aへ説明されたアプローチを使用することによる分離結果を示す。図5B〜図5Dは、それぞれ、分離された圧縮波502、せん断波504、及びストンリー波506を示し、該説明されたアプローチは、エイリアシング及びノイズといった実際上の脅威にもかかわらず高品質の分離結果を生成できることを示す。
図6は、いくつかの実施に係る、本願で説明されるような記述されたアルゴリズム、方法、機能、プロセス、フロー、及び手順に関連付けられた計算機能を提供するために使用された例示的なコンピュータシステム600のブロック図である。図示されたコンピュータ602は、サーバ、デスクトップコンピュータ、ラップトップ/ノートブックコンピュータ、無線データポート、スマートフォン、パーソナルデータアシスタント(PDA)、タブレットコンピューティングデバイス、これらのデバイスのうちの1つ若しくは複数のプロセッサ、又は任意の他の適切な処理デバイスといった任意のコンピューティングデバイスを包含することを意図しており、該コンピューティングデバイスの物理インスタンス又は仮想インスタンス(又はその両方)を含む。加えて、コンピュータ602は、キーパッド、キーボード、タッチスクリーン、又はユーザ情報を受け入れ可能な他のデバイスといった入力デバイスと、コンピュータ602の動作に関連付けられた情報を伝達する出力デバイスとを含むコンピュータを包含していてもよく、デジタルデータ、視覚情報、または音声情報(若しくは情報の組み合わせ)、又はグラフィカルユーザインタフェース(GUI)を含む。
コンピュータ602は、クライアント、ネットワークコンポーネント、サーバ、データベースあるいは他の持続性のもの、又は本願に説明された主題を実行するためのコンピュータシステムの任意の他のコンポーネント(又は役割の組み合わせ)としての役割を果たすことができる。例示のコンピュータ602は、ネットワーク630と通信可能に結合される。いくつかの実施では、コンピュータ602の1又は複数のコンポーネントは、クラウドコンピューティングベースの、ローカルな、グローバルな、又は他の環境(又は環境の組み合わせ)を含む環境内で動作するように構成できる。
高レベルでは、コンピュータ602は、説明された主題に関連付けられたデータ及び情報を受け、送出し、処理し、保存し、又は管理するように動作可能な電子的計算デバイスである。いくつかの実施によれば、コンピュータ602は、アプリケーションサーバ、電子メールサーバ、ウェブサーバ、キャッシングサーバ、ストリーミングデータサーバ、又は他のサーバ(又はサーバの組み合わせ)を含んでいても、或いは通信可能に結合されていてもよい。
コンピュータ602は、(例えば、別のコンピュータ602上で動作する)クライアントアプリケーションからネットワーク630を介してリクエストを受けることができ、受信したリクエストを適切なソフトウェアアプリケーションを用いて処理することによって受信したリクエストに応答することができる。加えて、リクエストは、コンピュータ602へ、(例えば、コマンドコンソールから又は他の適切なアクセス方法によって)内部ユーザ、外部ユーザ又はサードパーティ、他の自動化されたアプリケーション、更には、その他の適切なエンティティ、個人、システム、又はコンピュータから送信されるようにしてもよい。
コンピュータ602の各コンポーネントは、システムバス603を用いて通信できる。いくつかの実施では、コンピュータ602のコンポーネントのいずれか又は全ては、ハードウェア又はソフトウェアの両方(又はハードウェア及びソフトウェアの組み合わせ)が、アプリケーションプログラミングインタフェース(API)612又はサービスレイヤ613(又は、API612及びサービスレイヤ613の組み合わせ)を用いて、システムバス603を介して、互いに又はインタフェース604(又は両方の組み合わせ)とインタフェース接続することができる。API612は、ルーチン、データ構造、及びオブジェクトクラスの仕様を含むことができる。API612は、コンピュータ言語に依存しないか又は依存するかのいずれかであり、また完全なインタフェース、単一の関数、又は一組のAPIを指していてもよい。サービスレイヤ613は、コンピュータ602へソフトウェアサービスを、又はコンピュータ602に通信可能に結合された他のコンポーネントを(図示されているかどうかに関わらず)提供する。コンピュータ602の機能は、このサービスレイヤを用いて全てのサービスコンシューマにアクセス可能であってもよい。サービスレイヤ613によって提供されたサービスといったソフトウェアサービスは、定義済みのインタフェースを介して再利用可能な定義済み機能を提供する。例えば、インタフェースは、JAVA(登録商標)、C++、又は、拡張マークアップ言語(XML)フォーマット若しくは他の適切なフォーマットでデータを提供する他の適切な言語で書かれたソフトウェアであり得る。コンピュータ602の統合コンポーネントとして図示される一方で、代替の実施は、API612又はサービスレイヤ613を、コンピュータ602の他のコンポーネント又はコンピュータ602に通信可能に結合された他のコンポーネントに(図示されているかどうかに関わらず)関連してスタンドアロンコンポーネントとして例示することができる。さらには、API612又はサービスレイヤ613のいずれか又は全ての部分は、本願の範囲から逸脱することなく、別のソフトウェアモジュール、エンタープライズアプリケーション、又はハードウェアモジュールの子モジュール又はサブモジュールとして実施されるようにしてもよい。
コンピュータ602は、インタフェース604を含む。図6には単一のインタフェース604として図示されているが、2つ又はそれより多くのインタフェース604を、コンピュータ602の特定のニーズ、要望、又は特定の実施に従って使用することができる。インタフェース604は、分散環境においてネットワーク630に(図示されているかどうかに関わらず)接続される他のシステムと通信するために、コンピュータ602によって使用される。一般にインタフェース604は、ソフトウェア又はハードウェア(又はソフトウェア及びハードウェアの組み合わせ)でエンコードされたロジックを含み、またネットワーク630と通信するように動作可能である。より具体的には、インタフェース604は、ネットワーク630又はインタフェースのハードウェアが図示のコンピュータ602の内側及び外側において物理信号を通信するように動作可能であるように、通信に関連付けられた1又は複数の通信プロトコルをサポートするソフトウェアを含むことができる。
コンピュータ602は、プロセッサ605を含む。図6においては単一のプロセッサ605として示されるが、2又はそれより多くのプロセッサを、コンピュータ602の特定のニーズ、要望、又は特定の実施に従って使用することができる。一般に、プロセッサ605は、命令を実行し、また、本願に記載されるように、コンピュータ602の動作、及び任意のアルゴリズム、方法、機能、プロセス、フロー、若しくは手順を実行するためにデータを処理する。
また、コンピュータ602は、コンピュータ602又はネットワーク630に(図示されているかどうかに関わらず)接続可能な他のコンポーネント(又はこの両方の組み合わせ)のためのデータを保持することができるデータベース606を含む。例えば、データベース606は、インメモリの、従来型の、又は他のタイプのデータベースであることができ、該データベースは、本願に整合するデータを格納する。いくつかの実施では、データベース606は、説明された機能及びコンピュータ602の特定のニーズ、要望、又は特定の実施に従う、2以上の異なるデータベースタイプの組み合わせ(例えば、インメモリと従来型とのハイブリッド型データベース)であり得る。図6において単一のデータベース606として示されるが、(同じタイプ又はタイプの組み合わせの)2以上のデータベースは、説明された機能及びコンピュータ602の特定のニーズ、要望、又は特定の実施に従って、使用できる。データベース606は、コンピュータ602の不可欠な構成要素として示されるが、代替の実施では、データベース606は、コンピュータ602の外部にあってもよい。例えば、データベース606は、FWSデータを保持することができる。
また、コンピュータ602は、ネットワーク630に(図示されているかどうかに関わらず)接続できるコンピュータ602又は他のコンポーネント(又は両方の組み合わせ)のためのデータを保持可能なメモリ607を含む。例えば、メモリ607は、ランダムアクセスメモリ(RAM)、読み取り専用メモリ(ROM)、光学的なもの、磁気的なもの等であることができ、本願に整合するデータを格納する。いくつかの実施では、メモリ607は、説明された機能及びコンピュータ602の特定のニーズ、要望又は特定の実施に従う、2又はより多くの異なるタイプのメモリの組み合わせ(例えば、RAMと磁気記憶装置との組み合わせ)であることができる。図6には単一のメモリ607として示されるが、(同じタイプ又はタイプの組み合わせの)2以上のメモリ607が、説明された機能及びコンピュータ602の特定のニーズ、要望、又は特定の実施に従って使用できる。メモリ607は、コンピュータ602の不可欠な構成要素として示されるが、代替の実施では、メモリ607はコンピュータ602の外部にあってもよい。
アプリケーション608は、コンピュータ602の特定のニーズ、要望、又は特定の実施に従う機能を、特に本願において説明された機能に関して提供するアルゴリズムソフトウェアエンジンである。例えば、アプリケーション608は、1又は複数のコンポーネント、モジュール、又はアプリケーションとして働くことができる。さらに、単一のアプリケーション608として示されているが、アプリケーション608は、コンピュータ602上において複数のアプリケーション608として実施されていてもよい。加えて、コンピュータ602と一体として図示されているが、代替の実施では、アプリケーション608は、コンピュータ602の外部にあってもよい。
コンピュータ602を含むコンピュータシステムに関連付けられ又は外部にあるコンピュータ602が任意の数あってよく、各コンピュータ602は、ネットワーク630を介して通信する。さらに、語句「クライアント」、「ユーザ」、及び他の適切な用語は、本願の範囲から逸脱することなく、適宜交換可能に使用されてもよい。さらには、本願は、多くのユーザが1台のコンピュータ602を使用してよいこと、又は1人のユーザが複数のコンピュータ602を使用してよいことを想定している。
本明細書で記述された主題及び機能的な操作の実施は、デジタル電子回路、有形に具体化されたコンピュータのソフトウェア又はファームウェア、コンピュータハードウェアに、又はそれらの一又は複数の組み合わせで実施でき、これらデジタル電子回路、ソフトウェア、ファームウェア、及びコンピュータハードウェアは、本明細書に開示される構造及びそれらの構造的な同等物を含む。本明細書で記述された主題に係る実施は、1又は複数のコンピュータプログラム、つまり、コンピュータプログラム命令の一又は複数のモジュールとして実施されることができ、コンピュータプログラム命令は、データ処理装置により実施されるシステムによる実行のために、有形の非一時的なコンピュータ読取可能なコンピュータ記憶媒体上にエンコードされ、又はコンピュータ若しくはコンピュータにより実施されるシステムの動作を制御する。代替的に又は追加的に、プログラム命令は、人工的に生成された伝播信号、例えば、マシン生成の電気的、光学的、又は電磁気的な信号にエンコードされ、この信号は、データ処理装置による実行のために受振器装置への送信用の情報をエンコードするために生成される。コンピュータ記憶媒体は、機械により読取可能な記憶装置、機械により読取可能な記憶基板、ランダム若しくはシリアルアクセスのメモリデバイス、又はコンピュータ記憶媒体の組み合わせであり得る。
「リアル−タイム」、「リアルタイム」、「リアル(ファースト)タイム(RFT)」、「ほぼリアルタイム(NRT)」、「準リアルタイム」、又は(当業者が理解する)類似の用語は、個々が、実質的に同時に動作及び応答が生じることを知覚するように、該動作及び応答が時間的に近接していることを意味する。例えば、データにアクセスするための個々の動作に続くデータの表示(又は表示の開始)に対する応答の時間差は、1ms未満、1sec未満、又は5sec未満であり得る。要求されたデータが、即時に表示(又は表示開始)される必要はないが、記述されたコンピューティングシステムの処理制限と、データを例えば収集し、正確に評価し、分析し、処理し、格納し、又は送信するために必要な時間とを考慮して、意図的な遅延なしに表示(又は表示が開始)される。
用語「データ処理装置」、「コンピュータ」、又は「電子コンピュータデバイス」(又は当業者によって理解される同等の用語)は、データ処理ハードウェアを指し、あらゆる種類の装置、デバイス、及びマシンを包含し、これらの種類は、データを処理するためのものであり、例示として、プログラマブルプロセッサ、コンピュータ、又は複数のプロセッサ若しくはコンピュータを含む。コンピュータは、例えば、中央処理装置(CPU)、FPGA(フィールドプログラマブルゲートアレイ)、又はASIC(特定用途向け集積回路)を包含する専用ロジック回路であるか、又はこれらを更に含むことができる。いくつかの実施では、データ処理装置、又は専用ロジック回路(又はデータ処理装置と専用ロジック回路との組み合わせ)は、ハードウェア系又はソフトウェア系(又はソフトウェア系及びハードウェア系の両方の基づく組み合わせ)であることができる。装置は、オプションとして、コンピュータプログラムの実行環境を生成するコード、例えば、プロセッサファームウェア、プロトコルスタック、データベース管理システム、オペレーティングシステム、又は実行環境の組み合わせを構成するコードを含むことができる。本願は、例えば、LINUX、UNIX(登録商標)、WINDOWS(登録商標)、MACOS、ANDROID(登録商標)、IOS、又は他の任意の適切なオペレーティングシステムといった従来のオペレーティングシステムを伴う又は伴わないデータ処理装置を使用することを意図している。
コンピュータプログラムは、プログラム、ソフトウェア、ソフトウェアアプリケーション、モジュール、ソフトウェアモジュール、スクリプト、又はコードとして参照され又は記述されるものであって、コンパイル言語若しくはインタープリタ言語を含むプログラミング言語、又は宣言型若しくは手続き型の言語の形式で記載でき、またコンピュータプログラムは、コンピューティング環境で使用するための任意の形式で展開でき、スタンドアロンプログラム、モジュール、コンポーネント、サブルーチン、又は他の適切なユニットとして含む。コンピュータプログラムは、ファイルシステム内のファイルに対応できるが、必ずしも対応する必要はない。プログラムは、他のプログラム又はデータを保持するファイルの一部、例えばマークアップ言語ドキュメントに格納された一又は複数のスクリプトに保存でき、該他のプログラムは、問題のプログラム専用の単一ファイル内に、或いは複数の連携したファイル、例えば一又は複数のモジュール、サブプログラム、又はコードの一部に格納できる。コンピュータプログラムは、一又は複数のコンピュータ上に展開でき、該コンピュータは、一のサイトに位置し、又は複数のサイトに分散されて配置され、これらは通信ネットワークによって相互接続される。様々な図面に例示されたプログラムの部分は、個々のモジュールとして示され、様々なオブジェクト、方法、又は他のプロセスを介して、種々の特徴及び機能を実施しているが、プログラムは、それらの代わりに、いくつかのサブモジュール、サードパーティのサービス、コンポーネント、ライブラリ、及びそれらを、必要に応じて含むことができる。逆に、必要に応じて、様々なコンポーネントの機能及び特徴は、単一のコンポーネントに組み合わせることができる。計算上の特定を為すために使用されたしきい値は、静的に、動的に、又は静的と動的との両方で特定できる。
本明細書で説示された方法、プロセス、又は論理フローは、一又は複数のプログラム可能なコンピュータによって実行でき、該コンピュータは、入力データを操作して出力データを生成することによって機能を実行する一又は複数のコンピュータプログラムを実行する。方法、プロセス、又は論理フローは、特定用途のロジック回路、例えばCPU、FPGA、又はASICとして実行でき、また装置も、特定用途のロジック回路、例えばCPU、FPGA、又はASICとして実施できる。
コンピュータプログラムの実行のための適切なコンピュータは、汎用又は特定用途のマイクロプロセッサ、これらの両方、又は別の種類CPUに基づくことができる。一般的には、CPUは、読み取り専用メモリ(ROM)又はランダムアクセスメモリ(RAM)あるいはその両方から命令及びデータを受け取る。コンピュータの重要な要素は、命令を行い又は実行するためのCPU、及び命令及びデータを保存するための一又は複数のメモリデバイスである。一般に、コンピュータは、データを保存するための一又は複数の大容量記憶装置、例えば磁気、光磁気ディスク、又は光ディスクを含み、或いはこれらの大容量記憶装置に動作可能に結合され、データを受信し、転送し、又はこの両方を行う。しかし、コンピュータは、そのようなデバイスを持つ必要はない。さらには、コンピュータは、別のデバイス、例えば携帯電話、携帯情報端末(PDA)、モバイルオーディオ若しくはビデオプレーヤー、ゲームコンソール、全地球測位システム(GPS)受信機、又はポータブルストレージデバイス、例えばユニバーサルシリアルバス(USB)、フラッシュドライブ等に組み込むことができる。
コンピュータプログラム命令及びデータを格納するための適切な(必要に応じて一時的又は非一時的な)コンピュータ読取可能媒体は、例示として半導体メモリデバイス、磁気デバイス、及び光学メモリデバイスを含む、あらゆる形態のメモリ、媒体、及びメモリデバイスを含む。あらゆる形態のメモリデバイスは:半導体メモリデバイス、例えば、消去可能プログラマブルリードオンリーメモリ(EPROM)、電気的消去可能プログラマブル読み取り専用メモリ(EEPROM)、及びフラッシュメモリデバイスと;磁気ディスク、例えば、内部/ハードディスク又はリムーバブルディスクと;光磁気ディスクと;CD−ROM、DVD+/−R、DVD−RAM、及びDVD−ROMディスクと;を含む。メモリは、様々なオブジェクト又はデータを格納でき、オブジェクト又はデータは、キャッシュ、クラス、フレームワーク、アプリケーション、バックアップデータ、ジョブ、Webページ、Webページテンプレート、データベーステーブル、動的情報を格納するリポジトリ、及び他の適切な情報を包含しており、他の適切な情報は、任意のパラメータ、任意の変数、任意のアルゴリズム、任意の命令、任意のルール、任意の制約、又は任意の参照を含む。さらに、メモリは、ログ、ポリシー、セキュリティ又はアクセスのデータ、レポートファイル、同様の他のものといった他の適切なデータを含むことができる。プロセッサ及びメモリは、特定用途ロジック回路によって補完され又は組み込むことができる。
ユーザとのインタラクションを提供するために、本明細書で記載された主題の実施は、ユーザに情報を表示するためのディスプレイデバイス、例えばCRT(陰極線管)、LCD(液晶ディスプレイ)、LED(発光ダイオード)又はプラズマモニターと、ユーザがコンピュータに入力を提供するキーボード及びポインティングデバイス、例えばマウス、トラックボール、トラックパッドとを有するコンピュータ上で実施できる。入力が、タブレットコンピュータの感圧性表面、静電容量式若しくは電気式のセンシングを使用したマルチタッチスクリーン、又は別のタイプのタッチスクリーンといった、タッチスクリーンを用いてコンピュータに提供することもできる。他の種類のデバイスが、ユーザとのインタラクションを提供するために使用できる。同様に、例えば、ユーザに提供されるフィードバックは、例えば、視覚フィードバック、聴覚フィードバック、触覚フィードバック、及びあらゆる形態の感覚のフィードバックであり得る。ユーザからの入力は、音響入力、音声入力、触覚入力を包含する、あらゆる形式で受信できる。さらに、コンピュータは、ユーザによって使用されるクライアントコンピューティングデバイスからのドキュメントを受信し及び該ドキュメントをデバイスに送信することによって(例えば、Webブラウザから受けたリクエストに応答してユーザのクライアントデバイス上においてWebブラウザにWebページを送ることにより)ユーザとインタラクションできる。
用語「グラフィカルユーザインタフェース」又は「GUI」は、単数又は複数で使用して、一又は複数のグラフィカルユーザインタフェース及び特定のグラフィカルユーザインタフェースのディスプレイの各々を説明できる。したがって、GUIは、Webブラウザ、タッチスクリーン、又はコマンドラインインタフェース(CLI)を含むが、これらに限定されることなく、任意のグラフィカルユーザインタフェースを提示することができ、コマンドラインインタフェースは、情報を処理すると共に情報の結果をユーザに効率的に提示する。一般的に、GUIは、いくつかユーザインタフェース(UI)要素、例えばインタラクティブフィールド、プルダウンリスト、ボタンを含むことができ、これらのうちのいくつか又は全ては、Webブラウザに関連付けられる。これら及び他のUI要素は、Webブラウザの機能に関連付けられまたWebブラウザの機能を表すことができる。
本明細書で記述された主題の実施は、バックエンドコンポーネントを例えばデータサーバとして含むコンピューティングシステム、ミドルウェアコンポーネントを、例えばアプリケーションサーバとして含むコンピューティングシステム、フロントエンドコンポーネントを、例えばクライアントコンピュータとして含むコンピューティングシステム、及び、一又は複数のこのようなバックエンド、ミドルウェア、又はフロントエンドコンポーネントの任意の組み合わせを含むコンピューティングシステムにおいて実施でき、クライアントコンピュータは、ユーザが本明細書においで記述される主題の実施とインタラクションできるグラフィカルユーザインタフェース又はWebブラウザを有する。システムのコンポーネントは、有線又は無線のデジタルデータ通信(又はデータ通信の組み合わせ)、例えば通信ネットワークの任意の形式又は媒体によって相互接続されることができる。通信ネットワークの例示は、ローカルエリアネットワーク(LAN)、無線アクセスネットワーク(RAN)、メトロポリタンエリアネットワーク(MAN)、広域ネットワーク(WAN)、マイクロ波アクセスの世界的な相互運用性(WIMAX)、例えば802.11a/b/g/n又は802.20(又は802.11xと802.20の組み合わせ、又は本開示に整合する他のプロトコル)を用いる無線ローカルエリア(WLAN)、インターネットの全て若しくは一部、別の通信システム、又は1又は複数の位置におけるシステム(又は通信ネットワークとの組み合わせ)を含む。ネットワークは、例えば、インターネットプロトコル(IP)パケット、フレームリレーフレーム、非同期転送モード(ATM)セル、音声、ビデオ、データ、又はネットワークアドレス間の他の情報(又は通信タイプの組み合わせ)と通信してよい。
コンピューティングシステムは、クライアント及びサーバを含むことができる。クライアント及びサーバは、一般的には、互いに離れており、また典型的には、通信ネットワークを介してインタラクションする。クライアント及びサーバの関係は、それぞれのコンピュータ上で動作すると共に互にクライアント・サーバ関係にあるコンピュータプログラムのおかげで生じる。
本明細書は多くの特定の実施の詳細を含む一方で、これらは、任意の発明の範囲又は請求される得る範囲の制限として解釈されるべきではなく、むしろ、特定の発明の特定の実施固有の特徴の説明として解釈されるべきである。個別の実施の観点で本明細書において説明された機能は、組み合わせて又は単一の実施で、実現できる。逆に、単一の実施の観点で既述された様々な特徴は、複数の実施で、個別に、又は任意のサブコンビネーションで実施されることもできる。さらには、既述の特徴は特定の組み合わせで動作するものとして説示され、最初はそのようなものとして請求されているが、一又は複数の特徴は、請求された組み合わせから場合によっては削除され、またサブコンビネーション又はサブコンビネーションの変形であり得る。
主題の特定の実施が説明された。記載された実施の置換、変更、及び他の実施は、当業者には明らかである以下の請求の範囲内である。操作は特許請求の範囲に又は特定の順序で図面に描かれている一方で、これは、望ましい結果を達成するために、示された特定の順序で又は順番でそのような操作が実行されること又は全ての図示された操作が実行される(いくつかの操作は随意的と見なされる)ことを要求するものとして理解されるべきではない。ある状況では、マルチタスク又は並列処理(又はマルチタスク及び並列処理の組み合わせ)を行うことが有利であり、適切と思われる場合に実行される。
さらに、既述の実施における様々なシステムモジュール及びコンポーネントの分離又は統合が、全ての実施においてそのような分離又は統合を必要とするものとして理解されるべきではない。また、記載されたプログラムコンポーネント及びシステムは、一般的には、単一のソフトウェア製品に統合でき、又は複数のソフトウェア製品にパッケージ化できることが理解されるべきである。
これに従って、既述の例示的な実施は、本願を定義し又は制約しない。他の変更、置換、及び変更も、本開示の範囲及び精神から逸脱することなく可能である。
さらには、請求された任意の実施は、少なくとも、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータシステムに適用可能であると考えられ、非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータにより読取可能な命令を格納し、またコンピュータシステムは、ハードウェアプロセッサに相互に動作可能に結合されたコンピュータメモリを含み、ハードウェアプロセッサは、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法又は命令を実行するように構成される。
600 コンピュータシステム
602 コンピュータ
605 プロセッサ
607 メモリ

Claims (20)

  1. 全波形音波(FWS)データを受信するステップと;
    前記受信したFWSデータにアンチエイリアシング線形ラドン変換を実行するステップと;
    スラントウィンドウを用いて、波成分に対応するラドン変換済みFWSデータを抽出するステップと;
    前記抽出されたラドン変換済みFWSデータに逆ラドン変換を実行することによって、前記波成分の信号を特定するステップと;を備える、
    方法。
  2. 前記波成分は、圧縮波、せん断波、又はストンリー波のうちの1つである、
    請求項1に記載の方法。
  3. 前記アンチエイリアシング線形ラドン変換を実行するステップは:
    アンチエイリアシングオフセット間隔に基づき、一組のアンチエイリアシングオフセット値及び一組のslowness値を特定するステップと;
    前記受信したFWSデータ及び前記一組のアンチエイリアシングオフセット値に基づき、補間されたFWSデータを生成するステップと;
    前記一組のアンチエイリアシングオフセット値及び前記一組のslowness値に基づき、前記補間されたFWSデータにラドン変換を実行するステップと;を備える、
    請求項1に記載の方法。
  4. 前記アンチエイリアシングオフセット間隔は△x=0.5*△t*Vminであり、
    ここで、△tは、前記受信したFWSデータのサンプリング時間間隔であり、Vminは、前記波成分の速度の下限である、
    請求項3に記載の方法。
  5. 前記波成分の中心速度、速度範囲、中心時間、及び時間範囲を特定するために、速度解析を実行するステップを更に備える、
    請求項1に記載の方法。
  6. 前記スラントウィンドウは、前記波成分の前記中心速度及び前記中心時間に基づき特定された点を中心としており、前記スラントウィンドウは、前記波成分の前記速度範囲及び前記時間範囲に基づき特定されたサイズを有する、
    請求項5に記載の方法。
  7. 前記抽出されたラドン変換FWSデータに逆ラドン変換を実行するステップは、一組の実際のオフセット値に基づき前記逆ラドン変換を実行するステップを備える、
    請求項1に記載の方法。
  8. 前記受信したFWSデータは、共通ショットギャザーである、
    請求項1に記載の方法。
  9. コンピュータメモリと;
    前記コンピュータメモリに相互利用可能に結合され、所定の動作を実行するように構成された1又は複数のハードウェアプロセッサと;を備え、
    前記動作は、
    全波形音波(FWS)データを受信する処理と;
    前記受信したFWSデータにアンチエイリアシング線形ラドン変換を実行する処理と;
    斜されたウィンドウを用いて、波成分に対応するラドン変換済みFWSデータを抽出する処理と;
    前記抽出されたラドン変換済みFWSデータに対して逆ラドン変換を実行することによって、前記波成分の信号を特定する処理と;を備える、
    システム。
  10. 前記波成分は、圧縮波、せん断波、又はストンリー波のうちの1つである、
    請求項9に記載のシステム。
  11. 前記アンチエイリアシング線形ラドン変換を実行することは:
    アンチエイリアシングオフセット間隔に基づき、一組のアンチエイリアシングオフセット値及び一組のslowness値を特定する処理と;
    前記受信したFWSデータ及び前記一組のアンチエイリアシングオフセット値に基づき、補間されたFWSデータを生成する処理と;
    前記一組のアンチエイリアシングオフセット値及び前記一組のslowness値に基づき、前記補間されたFWSデータにラドン変換を実行する処理と;を備える、
    請求項9に記載のシステム。
  12. 前記アンチエイリアシングオフセット間隔は△x=0.5*△t*Vminであり、ここで、△tは、前記受信したFWSデータのサンプリング時間間隔であり、Vminは、前記波成分の速度の下限である、
    請求項11に記載のシステム。
  13. 前記動作は、前記波成分の中心速度、速度範囲、中心時間、及び時間範囲を特定するために速度解析を実行することを更に備える、
    請求項9に記載のシステム。
  14. 前記スラントウィンドウは、前記波成分の前記中心速度及び前記中心時間に基づき特定された点を中心としており、
    前記スラントウィンドウは、前記波成分の前記速度範囲及び前記時間範囲に基づき特定されるサイズを有する、
    請求項13に記載のシステム。
  15. 非一時的なコンピュータ読取可能媒体であって、コンピュータシステムによって実行可能な1又は複数の命令を所定の動作を実行するために格納し、前記動作は:
    全波形音波(FWS)データを受信する処理と;
    前記受信したFWSデータにアンチエイリアシング線形ラドン変換を実行する処理と;
    スラントウィンドウを用いて、波成分に対応するラドン変換済みFWSデータを抽出する処理と;
    前記抽出されたラドン変換済みFWSデータに逆ラドン変換を実行することによって、前記波成分の信号を特定する処理と;を備える、
    非一時的なコンピュータ読取可能媒体。
  16. 前記波成分は、圧縮波、せん断波、又はストンリー波のうちの1つである、
    請求項15に記載の非一時的なコンピュータ読取可能媒体。
  17. 前記アンチエイリアシング線形ラドン変換を実行することは:
    アンチエイリアシングオフセット間隔に基づき、一組のアンチエイリアシングオフセット値及び一組のslowness値を特定する処理と;
    前記受信したFWSデータ及び前記一組のアンチエイリアシングオフセット値に基づき、補間されたFWSデータを生成する処理と;
    前記一組のアンチエイリアシングオフセット値及び前記slowness値に基づき、前記補間されたFWSデータにラドン変換を実行する処理と;を備える、
    請求項15に記載の非一時的なコンピュータ読取可能媒体。
  18. 前記アンチエイリアシングオフセット間隔は△x=0.5*△t*Vminであり、ここで、△tは前記受信したFWSデータのサンプリング時間間隔であり、Vminは、前記波成分の速度の下限である、
    請求項17に記載の非一時的なコンピュータ読取可能媒体。
  19. 前記動作は、前記波成分の中心速度、速度範囲、中心時間、及び時間範囲を特定するために、速度解析を実行する処理をさらに備える、
    請求項15に記載の非一時的なコンピュータ読取可能媒体。
  20. 前記スラントウィンドウは、前記波成分の前記中心速度及び前記中心時間に基づき特定された点を中心としており、
    前記スラントウィンドウは、前記波成分の前記速度範囲及び前記時間範囲に基づき特定されたサイズを有する、
    請求項19に記載の非一時的なコンピュータ読取可能媒体。
JP2019562287A 2017-05-11 2018-05-07 全波形音波動場分離のための処理方法 Active JP6943979B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US15/592,979 US10564304B2 (en) 2017-05-11 2017-05-11 Processing methodology for full-waveform sonic wavefield separation
US15/592,979 2017-05-11
PCT/US2018/031399 WO2018208676A1 (en) 2017-05-11 2018-05-07 Processing methodology for full-waveform sonic wavefield separation

Publications (2)

Publication Number Publication Date
JP2020519889A true JP2020519889A (ja) 2020-07-02
JP6943979B2 JP6943979B2 (ja) 2021-10-06

Family

ID=62218341

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019562287A Active JP6943979B2 (ja) 2017-05-11 2018-05-07 全波形音波動場分離のための処理方法

Country Status (7)

Country Link
US (2) US10564304B2 (ja)
EP (1) EP3622327B1 (ja)
JP (1) JP6943979B2 (ja)
CN (1) CN110832355B (ja)
CA (1) CA3063129C (ja)
SA (1) SA519410515B1 (ja)
WO (1) WO2018208676A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10564304B2 (en) * 2017-05-11 2020-02-18 Saudi Arabian Oil Company Processing methodology for full-waveform sonic wavefield separation
US11353612B2 (en) 2019-03-11 2022-06-07 Saudi Arabian Oil Company Nonstationary maximum likelihood method to estimate dispersion spectra for full waveform sonic logging
CN111856565B (zh) * 2019-04-26 2022-04-22 中国石油化工股份有限公司 一种利用自适应分析时窗提取地震属性的方法、存储介质
US11500116B2 (en) 2019-05-15 2022-11-15 Saudi Arabian Oil Company Identifying characteristics of a subterranean region using vector-based wavefield separation of seismic data from the subterranean region
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum
CN115951408A (zh) * 2023-01-03 2023-04-11 上海勘测设计研究院有限公司 地下地层的声波探测的慢度提取方法、装置、介质及设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09211142A (ja) * 1996-01-31 1997-08-15 Schlumberger Overseas Sa 音波検層方法及びシステム
US20020138202A1 (en) * 2001-01-23 2002-09-26 Martinez Ruben D. Weighted slant stack for attenuating seismic noise
US20030220744A1 (en) * 2002-05-21 2003-11-27 Schonewille Michel Albert Method for processing seismic data to attenuate multiples
US20040223412A1 (en) * 2003-04-25 2004-11-11 Schlumberger Technology Corporation, Incorporated In The State Of Texas Borehole sonic data processing method
US20050174885A1 (en) * 2005-03-04 2005-08-11 Robinson John M. Radon transformations for removal of noise from seismic data
JP2006122084A (ja) * 2004-10-26 2006-05-18 Hitachi Medical Corp Ct装置、画像処理装置、画像処理方法、プログラム、記録媒体。
US20090067286A1 (en) * 2007-09-12 2009-03-12 Schlumberger Technology Corporation Dispersion extraction for acoustic data using time frequency analysis
US20120004849A1 (en) * 2010-03-22 2012-01-05 Schlumberger Technology Corporation Efficient windowed radon transform

Family Cites Families (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4648039A (en) 1983-12-30 1987-03-03 Schlumberger Technology Corporation Compressional/shear wave separation in vertical seismic profiling
US4703460A (en) 1984-12-28 1987-10-27 Schlumberger Technology Corporation Method and apparatus for acoustic dipole direct shear wave well logging
US5970023A (en) * 1999-02-19 1999-10-19 Baker Hughes Incorporated Reducing aliasing artifacts in seismic data processing using sharpened (TAU-P-Q) transforms
EP1254383B1 (en) 2000-01-21 2005-08-24 Schlumberger Holdings Limited System and method for seismic wavefield separation
US6449560B1 (en) 2000-04-19 2002-09-10 Schlumberger Technology Corporation Sonic well logging with multiwave processing utilizing a reduced propagator matrix
US6625541B1 (en) 2000-06-12 2003-09-23 Schlumberger Technology Corporation Methods for downhole waveform tracking and sonic labeling
US6614716B2 (en) 2000-12-19 2003-09-02 Schlumberger Technology Corporation Sonic well logging for characterizing earth formations
US6636810B1 (en) * 2002-05-24 2003-10-21 Westerngeco, L.L.C. High-resolution Radon transform for processing seismic data
GB2409304B (en) 2003-12-19 2007-11-14 Westerngeco Ltd Processing geophysical data
US7492664B2 (en) 2005-10-31 2009-02-17 Baker Hughes Incorporated Method for processing acoustic reflections in array data to image near-borehole geological structure
US7953556B2 (en) * 2007-03-09 2011-05-31 Fairfield Industries Incorporated Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique
US7526385B2 (en) 2007-06-22 2009-04-28 Schlumberger Technology Corporation Method, system and apparatus for determining rock strength using sonic logging
US8521433B2 (en) 2008-04-09 2013-08-27 Schlumberger Technology Corporation Methods and systems for processing acoustic waveform data
US7916576B2 (en) 2008-07-16 2011-03-29 Westerngeco L.L.C. Optimizing a seismic survey for source separation
US8559270B2 (en) * 2008-08-15 2013-10-15 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US8395966B2 (en) 2009-04-24 2013-03-12 Westerngeco L.L.C. Separating seismic signals produced by interfering seismic sources
US8811114B2 (en) 2010-08-23 2014-08-19 Baker Hughes Incorporated Imaging of formation structure ahead of the drill-bit
US8880348B2 (en) 2010-09-21 2014-11-04 Baker Hughes Incorporated Radon migration of acoustic data
US9091788B2 (en) * 2011-11-14 2015-07-28 Cggveritas Services Sa Device and method for de-blending simultaneous shooting data with apex shifted radon transform
CN103149592A (zh) * 2013-03-07 2013-06-12 天津城市建设学院 一种变偏移距vsp波场分离方法
US10577915B2 (en) 2014-01-16 2020-03-03 Schlumberger Technology Corporation Sonic logging for assessing well integrity
US9772419B2 (en) * 2014-03-10 2017-09-26 Read As Decomposing full-waveform sonic data into propagating waves for characterizing a wellbore and its immediate surroundings
US10392934B2 (en) 2014-10-03 2019-08-27 Schlumberger Technology Corporation Method and apparatus for processing waveforms
CN104570095A (zh) * 2014-12-24 2015-04-29 李姗 一种基于Radon变换消除斜缆虚反射的方法
CN104932010B (zh) * 2015-06-09 2018-08-07 中国海洋石油集团有限公司 一种基于近道镶边稀疏Radon变换的绕射波分离方法
US10564304B2 (en) * 2017-05-11 2020-02-18 Saudi Arabian Oil Company Processing methodology for full-waveform sonic wavefield separation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09211142A (ja) * 1996-01-31 1997-08-15 Schlumberger Overseas Sa 音波検層方法及びシステム
US20020138202A1 (en) * 2001-01-23 2002-09-26 Martinez Ruben D. Weighted slant stack for attenuating seismic noise
US20030220744A1 (en) * 2002-05-21 2003-11-27 Schonewille Michel Albert Method for processing seismic data to attenuate multiples
US20040223412A1 (en) * 2003-04-25 2004-11-11 Schlumberger Technology Corporation, Incorporated In The State Of Texas Borehole sonic data processing method
JP2006122084A (ja) * 2004-10-26 2006-05-18 Hitachi Medical Corp Ct装置、画像処理装置、画像処理方法、プログラム、記録媒体。
US20050174885A1 (en) * 2005-03-04 2005-08-11 Robinson John M. Radon transformations for removal of noise from seismic data
US20090067286A1 (en) * 2007-09-12 2009-03-12 Schlumberger Technology Corporation Dispersion extraction for acoustic data using time frequency analysis
US20120004849A1 (en) * 2010-03-22 2012-01-05 Schlumberger Technology Corporation Efficient windowed radon transform

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MAELAND, EINAR: "Sampling, aliasing, and inverting the linear Radon transform", GEOPHYSICS, vol. 69, no. 3, JPN6021002905, 2004, US, pages 859 - 861, ISSN: 0004435179 *
TRAD, DANIEL, ET AL.: "Latest views of the sparse Radon transform", GEOPHYSICS, vol. 68, no. 1, JPN6021002906, 2003, US, pages 386 - 399, XP055218664, ISSN: 0004435180, DOI: 10.1190/1.1543224 *

Also Published As

Publication number Publication date
JP6943979B2 (ja) 2021-10-06
WO2018208676A1 (en) 2018-11-15
US20180329097A1 (en) 2018-11-15
CN110832355B (zh) 2021-10-22
EP3622327A1 (en) 2020-03-18
CA3063129C (en) 2022-07-26
US10585200B2 (en) 2020-03-10
US10564304B2 (en) 2020-02-18
CN110832355A (zh) 2020-02-21
SA519410515B1 (ar) 2022-12-13
EP3622327B1 (en) 2023-03-29
US20190011585A1 (en) 2019-01-10
CA3063129A1 (en) 2018-11-15

Similar Documents

Publication Publication Date Title
JP6943979B2 (ja) 全波形音波動場分離のための処理方法
JP7105809B2 (ja) 周波数領域における地震速度反転のための振幅に依存しない勾配の計算
JP2020511645A (ja) 大規模地震データセットの適応性があり堅牢なノイズの除去のための自動化されたシステム及び方法
US20180136353A1 (en) Use of wavelet cross-correlation for virtual source denoising
US10288755B2 (en) Seismic processing workflow for broadband single-sensor single-source land seismic data
US20210181363A1 (en) Mapping near-surface heterogeneities in a subterranean formation
JP2020521126A (ja) 地震データにおけるノイズの抑制
CN113348385A (zh) 弹性波逆时迁移中的avo成像条件
EP3969943A1 (en) Identifying characteristics of a subterranean region using vector-based wavefield separation of seismic data from the subterranean region
CN110799856B (zh) 使用与地下方位角和反射角相关的全波形反演产生速度模型
EP3586170B1 (en) Surface-scattered noise reduction
JP7358385B2 (ja) 地震画像の向上
JP7254778B2 (ja) リバースタイムマイグレーションにおける偽像除去
US20190146111A1 (en) Applying orthogonalization filtering to wavefield separation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200127

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210119

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210202

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20210430

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210609

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: 20210810

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210909

R150 Certificate of patent or registration of utility model

Ref document number: 6943979

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150