JP7105809B2 - 周波数領域における地震速度反転のための振幅に依存しない勾配の計算 - Google Patents

周波数領域における地震速度反転のための振幅に依存しない勾配の計算 Download PDF

Info

Publication number
JP7105809B2
JP7105809B2 JP2019564802A JP2019564802A JP7105809B2 JP 7105809 B2 JP7105809 B2 JP 7105809B2 JP 2019564802 A JP2019564802 A JP 2019564802A JP 2019564802 A JP2019564802 A JP 2019564802A JP 7105809 B2 JP7105809 B2 JP 7105809B2
Authority
JP
Japan
Prior art keywords
seismic
velocity
points
objective function
computer
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.)
Active
Application number
JP2019564802A
Other languages
English (en)
Other versions
JP2020521140A (ja
JP2020521140A5 (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 JP2020521140A publication Critical patent/JP2020521140A/ja
Publication of JP2020521140A5 publication Critical patent/JP2020521140A5/ja
Application granted granted Critical
Publication of JP7105809B2 publication Critical patent/JP7105809B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • 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/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • 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/362—Effecting static or dynamic corrections; Stacking
    • 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/52—Move-out correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/614—Synthetically generated data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

本願は、２０１７年５月２２日に出願された米国出願第６２／５０９，３００号、及び２０１８年４月１８日に出願された米国出願第１５／９５６，２６１号の優先権を主張し、その内容全体は参照により本明細書に組み込まれる。
本開示は、地震データ処理に関する。
速度反転（ｖｅｌｏｃｉｔｙ ｉｎｖｅｒｓｉｏｎ、速度インバージョン）は、典型的には、ターゲット領域（ｔａｒｇｅｔ ｄｏｍａｉｎ、目標ドメイン）の地下構造のイメージング（画像化）に使用することが可能な、ターゲット領域の信頼性のある地震速度モデルを構築するために使用される。広く使用されている速度反転法としては、記録された地震データ中の進行時間情報を利用する光線断層撮影法（ｒａｙ ｔｏｍｏｇｒａｐｈｙ）である。しかしながら、この光線断層撮影法は通常、手動で到着イベントを選択する必要があり、これは時間がかかり、実際には、しばしば３次元モデル構築を困難にする。さらに、光線ベースの反転は複雑な領域、例えば、複雑な地下構造の領域を処理することができない。波動方程式に基づく反転方法も使用される。しかしながら、完全な波形の反転には、速度反転を極小にするサイクルスキップの問題がある。相互相関に基づくミスフィット関数（ｍｉｓｆｉｔ ｆｕｎｃｔｉｏｎ、不適合関数）を用いた反転は、勾配（ｇｒａｄｉｅｎｔ、傾斜）に大きなアーチファクトを有し、真の速度モデルに収束する速度を遅くする。したがって、既存の反転方法は、時間がかかるか、又は信頼性がない。
本開示は、周波数領域における地震速度反転のための振幅に依存しない勾配の計算に係る、コンピュータにより実施される方法、コンピュータプログラム製品及びコンピュータシステムを含む、方法及びシステムについて説明する。
一実施では、一の領域に関連付けられた地震データが受信され、この領域は複数のポイントによって表される１又は複数の地下層を含み、各ポイントは地震速度に関連付けられる。複数のポイントにおける地震速度は、複数の勾配値に基づいて地震速度を反復的に更新することによって特定され、各勾配値はポイントに対応し、ポイントの位置における目的関数の勾配を評価することによって特定される。特定された地震速度に基づいて１つ以上の地下層の地震画像が生成され表示される。
前述の実施は、コンピュータにより実施される方法と、コンピュータにより実施される方法を実行するためのコンピュータ読取可能な命令を記憶する非一時的なコンピュータ読取可能媒体と、コンピュータにより実施される方法／非一時的なコンピュータ読取可能媒体に記憶された命令を実行するように構成されたハードウェアプロセッサと相互動作可能に結合されたコンピュータメモリを含むコンピュータにより実施されるシステムとを用いて実施することができる。
本開示で説明される主題は、周波数領域において地震データを効率的かつ効果的に処理することによって、ターゲット領域の地震速度を確実に特定することができる。説示されたアプローチは、地震データの振幅変化の影響とは無関係である目的関数の勾配に基づいて、速度を反復的に更新する。説示されたアプローチは初期の速度推定の品質とは無関係に、真の速度に迅速に収束することができる。特定された地震速度は、ターゲット領域の地震画像を生成するために使用することができる。生成された地震画像は、炭化水素坑井の掘削パラメータを特定する等、効果的なオイル及びガス探査に使用することができる。他の利点は、当業者には明らかであろう。
本願に記載された主題における１又は複数の実施の詳細は、添付図面及び引き続く明細書に説明される。本主題の他の特徴、態様及び利点は、明細書、図面及び特許請求の範囲から明らかになるであろう。
図１は、いくつかの実施に係る、周波数領域における地震速度反転のための振幅に依存することなく勾配を計算するための例示的な方法を示すフローチャートである。
図２は、いくつかの実施に係る、速度モデルを反復的に更新するための例示的な方法を示すフローチャートである。
図３は、いくつかの実施に係る、真の速度モデルを示した図である。
図４Ａ～図４Ｂは、いくつかの実施に係る、計算された勾配を示す図である。
図５Ａ～図５Ｃは、いくつかの実施に係る、ガウスモデル上で説明されたアプローチを使用したものを示す図である。
図６は、いくつかの実施に係る、本開示で説明されるような、説明されるアルゴリズム、方法、機能、プロセス、フロー及び手順に関連するコンピューティング機能を提供するために使用される例示的なコンピュータシステムを示すブロック図である。
種々の図面における類似の参照番号及び記号は、同様の要素を指す。
以下の詳細な説明は周波数領域における地震速度反転のための振幅に依存しない勾配計算について説明し、当業者が１又は複数の特定の実施の文脈において開示された主題を作成し使用することを可能とするために提示される。開示された実施の様々な修正、変更、及び置換を行うことができることは当業者に容易に明らかとなるであろう。開示された実施において定義された一般的な原理は、本開示の範囲から逸脱することなく、他の実施及び応用（アプリケーション）に適用されてよい。したがって、本開示は、説明された又は図示された実施に限定されることを意図するものではなく、開示された原理及び特徴と一致する最も広い範囲が与えられるべきである。
複数の震源位置において地震波をターゲット領域に送出し、複数の受振器位置で反射波を記録することによって、１又は複数の地下層を含むターゲット領域について、地震データを収集することができる。地震データは、ターゲット領域の地下構造を分析するための速度モデル、例えば３次元速度モデルを構築するために使用することができる。速度モデルは、ターゲット領域内の異なる位置における地震波速度を含むことができる。
高レベルでは、説示されたアプローチがタイムシフト（又は進行時間（ｔｒａｖｅｌｔｉｍｅ））情報を使用して、周波数領域で記録された地震データから、ターゲット領域における地震速度を特定する。説示されたアプローチは、目的関数（又はミスフィット関数）の勾配に基づいて、速度モデルを反復的に更新する。初期速度推定値の品質とは無関係に、真の速度への迅速な収束（例えば、真の速度から所定の閾値内の最終推定速度は、閾値の反復回数内で見つけることができる）のため、目的関数は、記録された地震データと、反復中の推定速度を使用して生成された合成地震データとの間の相互相関に基づく。目的関数は、進行時間情報を容易に抽出できるように、周波数領域で表現される。目的関数の勾配を導出する際には、振幅情報と進行時間情報とを分離し、振幅情報（すなわち、振幅又はパワースペクトルの速度反転への影響）を無視する。結果としての勾配は、記録された地震データに埋め込まれた運動学的情報（ｋｉｎｅｍａｔｉｃａｌ ｉｎｆｏｒｍａｔｉｏｎ）に純粋に依存する。説示されたアプローチは、良好なグローバル収束の自動速度反転を可能にする。
図１は、いくつかの実施に係る、周波数領域における地震速度反転のための振幅に依存することなく勾配を計算するための例示的な方法１００を示すフローチャートである。提示を明確にするために、以下の説明は、本開示における他の図の文脈（コンテキスト）で方法１００を一般的に説明する。例えば、方法１００は、図６に示すコンピュータシステムによって実行することができる。しかしながら、方法１００は、必要に応じて、システム、環境、ソフトウェア及びハードウェア、又はシステム、環境、ソフトウェア、及びハードウェアの組み合わせによって実行され得ることが理解されるであろう。いくつかの実施では、方法１００の様々なステップを、並列して、組み合わせて、ループで、又は任意の順序で実行することができる。
方法１００は、一の領域に関連する地震データが受信されるブロック１０２より開始する。この領域は、複数のポイント（地点）又は位置（ロケーション）によって表される１又は複数の地下層を含むことができる。各ポイントは、地震波がポイントを伝播するときの地震波速度を表す地震速度に関連付けることができる。例えば、領域は、長さ２キロメートル（ｋｍ）、幅２ｋｍ、及び深さ２ｋｍの３次元（３Ｄ）領域とすることができる。その領域を１０メートル（ｍ）×１０ｍ×１０ｍの立方体に分割すると、その領域は２０１３立方体コーナーポイントで表すことができる。言い換えれば、この領域の地震速度は、これら２０１３のポイント又は位置での速度で表すことができる。他の方法を使用して、領域を表すポイント又は位置のセットを見つけることもできる。
地震データは地表面の上又は下に配置された多数の受振器（例えば、ジオホン（ｇｅｏｐｈｏｎｅ）又はハイドロホン（ｈｙｄｒｏｐｈｏｎｅ））によって受信することができる。震源は地震波をターゲット領域に送ることができ、受振器はターゲット領域によって反射された波を記録することができる。震源は、例えばトラックによって牽引され、異なる位置で地震波を発生させることができる。例えば、震源は受振器が反射波を記録するために第１の位置で第１のショットを発射することができ、震源は、第２のショットを発射するために第２の位置に移動することができる。１つのショットに対応する各受振器での記録データは、トレース（ｔｒａｃｅ）と呼ばれる。例えば、震源が１０万の異なる場所で発射し、１０００の受振器がある場合、結果として得られる地震データセットは１０８のトレースを得ることができる。
ブロック１０４では、複数の勾配値に基づいて地震速度を反復的に更新することによって、ターゲット領域内の複数のポイントにおける地震速度が特定される。各勾配値はポイントに対応し、ポイントの位置における目的関数の勾配を評価することによって特定される。図２は、ブロック１０４の詳細な動作を示す。
図２は、いくつかの実施に係る、速度モデルを反復的に更新するための例示的な方法２００を示すフローチャートである。提示を明確にするために、以下の説明は一般に、本開示における他の図の文脈で方法２００を説明する。例えば、方法２００は、図６に示すコンピュータシステムによって実行することができる。しかしながら、方法２００は、必要に応じて、システム、環境、ソフトウェア及びハードウェア、又はシステム、環境、ソフトウェア及びハードウェアの組み合わせによって実行され得ることが理解されるであろう。いくつかの実施では、方法２００の様々なステップは、並列で、組み合わせて、ループで、又は任意の順序で実行することができる。
方法２００は、反復カウンタｋがゼロとして初期化されるブロック２０２より開始する。
ブロック２０４において、ターゲット領域の初期速度モデルｖ０が生成される。例えば、ターゲット領域を２０１３のポイントで表すことができれば、ポイントごとに初期速度が推定され、初期速度モデルは２０１３の初期速度が含まれる。いくつかの実施では、各ポイントにおける初期速度推定値をノーマルムーブアウト（ｎｏｒｍａｌ ｍｏｖｅｏｕｔ、ＮＭＯ）速度又は他の速度推定値とすることができる。
ブロック２０６において、ターゲット領域の合成地震データ（又は合成データトレース）が、ｖｋ、すなわちｋ回目の反復からの推定速度モデルを使用して計算される。例えば、初期合成地震データは、初期速度モデルｖ０に基づいて生成される。いくつかの実施では、合成トレースが、速度モデルｖｋを使用して音響波動場伝播を数値的にシミュレートすることにより計算することができる。
ブロック２０８において、式（１）の目的関数の値が計算され、目的関数は、記録されたデータトレースと合成トレースとの間の相互相関に基づく。
Figure 0007105809000001
 ここで、上線付きのＣ（以下、単に「Ｃ」と示す）は、Ｃ（ｘ，ω；ｘｓ）であり、以下のように表される。
Figure 0007105809000002
 また、
Figure 0007105809000003
 である。ここで、τはタイムシフトであり、ωは角周波数であり、ｄ（ｘ，ω；ｘｓ）及びｄ（ｘ，ｔ；ｘｓ）は、それぞれ周波数領域及び時間領域における震源位置ｘｓ及び受振器位置ｘにおける記録されたデータトレースであり、それぞれ、周波数領域及び時間領域における震源位置ｘｓ及び受振器位置ｘにおける合成データトレースであり、＊は共役オペレーション（ｃｏｎｊｕｇａｔｅ ｏｐｅｒａｔｉｏｎ、共役動作）を表す。１０８の記録されたトレースの初期の例では、ｋ回目の反復での速度モデルｖｋに基づいて、１０８個の合成トレースを生成することができる。１０８個の記録されたトレースと１０８個の合成トレースに基づき、式（１）の目的関数の値Ｊを計算することができる。
ブロック２１０において、目的関数の値Ｊが所定の閾値より小さいかどうかの特定が行われる。Ｊが閾値未満である場合には、方法２００は、ブロック２１８に進み、速度モデルｖｋが最終速度モデルとして出力される。Ｊが閾値以上である場合には、方法２００は、速度モデルをさらに更新するためにブロック２１２に進む。
ブロック２１２において、式（１）の目的関数の勾配が、ターゲット領域を表す指定されたポイント又は位置について計算される。言い換えれば、ターゲット領域を表す各ポイントについて勾配を計算することができる。目的関数の勾配を導出するために、ｄ（ｘ，ω；ｘｓ）及びｄ（ｘ，ｔ；ｘｓ）は、それぞれ、下記式Ｃ及びＤのように表すことができる。
Figure 0007105809000004
Figure 0007105809000005
 ここで、Ａｄ及びＡｐは、それぞれ、記録されたトレース及び合成トレースの振幅を表し、φｄ及びφｐは、それぞれ、記録されたトレース及び合成トレースの位相を表す。したがって、式（２）のＣ（ｘ，ω；ｘｓ）は、以下のように表すことができる。
Figure 0007105809000006
 ここで、
Figure 0007105809000007
 また、
Figure 0007105809000008
 である。
周波数領域のターゲットポイントｘ’での式（１）の目的関数の勾配は、次のように表すことができる。
Figure 0007105809000009
 ここで、ｘ’はターゲットポイント（デカルト座標など）の位置を表し、（（∂Ｃ）／∂ω）＊は残差を表し、（∂Ｃ）／∂ｖは勾配の計算において残差が逆伝播される波の経路を与える。式（４）を波の経路（∂Ｃ）／∂ｖに代入すると、次の式が得られる。
Figure 0007105809000010
式（６）は振幅及び位相に関連する成分を混合するため、式（５）の勾配は、後述するように、信頼できる勾配を提供することができない。説示されたアプローチでは、振幅変化は無視され（すなわち、式（６）の１／ＡＣ・∂ＡＣ／∂ｖを無視し）、周波数に依存しないタイムシフト、すなわち、Δφ（ｘ，ω；ｘｓ）＝ｉωΔτ（ｘ；ｘｓ）が想定される。その結果、Ｒｙｔｏｖ近似を用いることにより、式（６）は以下のようになる。
Figure 0007105809000011
式（７）を式（５）に代入すると、勾配は次のように表すことができる。
Figure 0007105809000012
 ここで、Δτは記録データと合成データの間のタイムシフトであり、以下を満たす。
Figure 0007105809000013
陰関数導関数（ｉｍｐｌｉｃｉｔ ｆｕｎｃｔｉｏｎ ｄｅｒｉｖａｔｉｖｅ）の規則に従って、式（８）の∂Δτ／∂ｖは、次のように表すことができる。
Figure 0007105809000014
 ここで、Ｅは下記式に示すものであって、
Figure 0007105809000015
 ターゲットポイントの位置ｘ’とは独立であり、ｆ（ω）は震源によって生成される震源ウェーブレット（例えば、リッカーウェーブレット（Ｒｉｃｋｅｒ ｗａｖｅｌｅｔ））を表し、Ｇ（ｙ，ω；ｘ）は、位置ｘにおけるインパルス震源による位置での観測された波動場を示すグリーン（Ｇｒｅｅｎ）関数である。上述したアプローチでは、振幅変化が無視されるため（つまり、Ａｄ（ｘ，ω；ｘｓ）＝Ａｐ（ｘ，ω；ｘｓ））、またｄ（ｘ，ω；ｘｓ）はｐ（ｘ，ω；ｘｓ）のタイムシフトバージョンである、すなわち、ｐ（ｘ，ω；ｘｓ）＝ｅｉωΔτ・ｄ（ｘ，ω；ｘｓ）であるため、Δτを特定する必要がないことに留意すべきである。したがって、式（９）は以下のようになる。
Figure 0007105809000016
 ここで、Ｅは下記に示すものである。
Figure 0007105809000017
 いくつかの実施では、勾配は、有限差分法又は他の方法を使用して、式（８）及び（１０）によって計算することができる。
新しい残差は次のように表すことができることに留意すべきである。
Figure 0007105809000018
 式（１１）のように、加算によりΔτは自動的に尊重され、速度誤差情報を明確に示す。Δτ又は式（１１）の符号は、推定速度が真の速度より高いか低いかを示す。式（８）の勾配は、時間領域において振幅及び位相情報を明示的に分離できない式（１２）及び（１３）における以下の従来の公式からは取得できないことに留意すべきである。
Figure 0007105809000019
Figure 0007105809000020
要約すると、各ポイントの位置ｘ’について、勾配ｇ（ｘ’）は、式（８）及び（１０）を使用して計算できる。例えば、ターゲット領域が２０１３ポイントある場合、２０１３個の勾配が計算される。
ブロック２１４において、式（１４）に示すように、最急降下法（ｓｔｅｅｐｅｓｔ ｄｅｓｃｅｎｔ ｍｅｔｈｏｄ）を使用して各ポイントの速度が更新される。
Figure 0007105809000021
 ここで、ｖｋ（ｘ’）は、ポイントｘ’のｋ回目の反復での推定速度であり、γｋ（ｘ’）＝－ｇ（ｘ’）は、式（１）の目的関数の最急降下方向であり、αｋは、ｋ回目の反復のステップ長であり、直線探索（ラインサーチ）法又はそれらの方法によって取得できる。
ブロック２１６において、反復カウンタｋが１加算され、方法２００はブロック２０６に戻り、最後の反復で得られた推定速度を使用して合成データトレースを計算し、速度を反復的に更新する。
ブロック２１８において最終速度が見つかると、ターゲット領域を表す指定されたポイントに最終速度を含めることによって、３Ｄ速度モデルを生成することができる。場合によっては、３Ｄ速度モデルを使用して、ターゲット領域の地震画像を生成することができる。地震画像は地下の地質学的特徴を示すことができ、炭化水素坑井の潜在的な位置を識別すること、及び炭化水素坑井の掘削パラメータを特定すること等、効果的なオイル及びガス探査のために使用することができる。地震画像及び掘削パラメータは、ユーザインタフェース上に表示することができる。
図３は、いくつかの実施に係る、真の速度モデル３００を示す。説示されたアプローチを用いて、数値実験を行った。数値実験では、１つの震源３０２と１つの受振器３０４の単純なモデルが使用される。実験における２Ｄのターゲット領域は、長さ２ｋｍ、幅２ｋｍである。１０ｍの間隔を使用して、ターゲット領域を２０１×２０１のグリッドに分割し、２０１×２０１個のポイントによって表すことができる。図３の水平軸（横軸）及び垂直軸（縦軸）は、それぞれｘ方向及びｙ方向のグリッドインデックスを表す。震源３０２は（０．１ｋｍ、１．０ｋｍ）の位置にあり、受振器３０４は（１．９ｋｍ、１．０ｋｍ）の位置にある。バー３０６は、グレースケールカラーと速度値の連続範囲との間のマッピングを示す。速度モデル３００は、速度が２０００メートル／秒である等速モデルである。
図４Ａ～図４Ｂは、いくつかの実施に係る、計算された勾配４００ａ及び４００ｂを示す。勾配４００ａ及び４００ｂは、図３に示すターゲット領域、震源、及び受振器について式（８）を使用して計算される。勾配４００ａ及び４００ｂは、それぞれ、初期推定速度が真の速度よりも低い場合及び高い場合の、ターゲット領域における２０１×２０１ポイントの勾配を示す。図４Ａ及び図４Ｂにおいて、白色領域４０２は負の勾配を示し、黒色領域４０４は正の勾配を示す。図４Ａの初期推定速度は真の速度よりも低いが、図４Ｂの初期推定速度は真の速度よりも高いので、２つの反対の勾配符号は反対の速度誤差を正しく示している。
図５Ａ～図５Ｃは、いくつかの実施に係る、ガウスモデル（Ｇａｕｓｓｉａｎ ｍｏｄｅｌ）上での説明されたアプローチの使用を示す。２次元のターゲット領域は、長さ３ｋｍ、深さ３ｋｍである。数値試験では、ターゲット領域が１０ｍの間隔を使用して３０１×３０１のグリッドに分割される。図５Ａ～図５Ｃの水平軸及び垂直軸は、それぞれの方向のグリッドインデックスを表す。バー５０６は、グレースケールカラーと速度値の連続範囲との間のマッピングを示す。この試験では、３１ショットが頂部に均等に分配される。各ショットに、３０１個の受振器が底部に配置される。ショット間隔及び受振器間隔は、それぞれ１００ｍ及び１０ｍである。リッカーウェーブレットは１０Ｈｚの支配的な周波を有する震源波形として使用される。図５Ａは、等速２ｋｍ／ｓを含む初期速度モデル５００ａを示す。図５Ｃは、均一なバックグラウンド速度２ｋｍ／ｓと、平均速度１．７ｋｍ／ｓの低速度異常５０２を含む真の速度モデル５００ｃを示す。図５Ｂは、異常５０４が中央で概ね回復している、説示されたアプローチによって得られた最終速度モデルを示す。
図６は、いくつかの実施に係る、本開示で説明されるような、説明されるアルゴリズム、方法、関数、プロセス、フロー、及び手順に関連するコンピューティング機能を提供するために使用される例示的なコンピュータシステム６００のブロック図である。図示されたコンピュータ６０２は、サーバ、デスクトップコンピュータ、ラップトップ／ノートブックコンピュータ、無線データポート、スマートフォン、パーソナルデータアシスタント（ＰＤＡ）、タブレットコンピューティングデバイス、又はコンピューティングデバイスの物理的又は仮想インスタンス（又は両方）を含むこれらのデバイス内の１又は複数のプロセッサ等の任意のコンピューティングデバイスを包含することが意図される。さらに、コンピュータ６０２は（ユーザ情報を受け入れることができるキーパッド、キーボード、又はタッチスクリーン等の）入力デバイスと、（例えば、グラフィカルユーザインタフェース（ＧＵＩ）上でデジタルデータ、ビジュアル、又はオーディオ情報（又は情報の組合せ）を伝達する）コンピュータ６０２の動作に関連する情報を伝達する出力デバイスと、を含むコンピュータを備えることができる。
コンピュータ６０２は、クライアント、ネットワークコンポーネント、サーバ、データベースとしての役割、又はこれらの役割の組合せで、本開示で説明される主題を実行するための役割を果たし得る。図示のコンピュータ６０２は、ネットワーク６３０と通信可能に結合されている。いくつかの実施では、コンピュータ６０２の１又は複数のコンポーネントが、クラウドコンピューティングベース、ローカル、グローバル、又はこれらの環境の組合せを含む環境内で動作するように構成され得る。
高レベルでは、コンピュータ６０２は、説示した主題に関連するデータ及び情報を受信し、送信し、処理し、記憶し、又は管理するように動作可能な電子コンピューティングデバイスである。いくつかの実施によれば、コンピュータ６０２はまた、アプリケーションサーバ、電子メールサーバ、ウェブサーバ、キャッシングサーバ、ストリーミングデータサーバ、又はサーバの組合せを含むか、又はそれらと通信可能に結合され得る。
コンピュータ６０２は、（例えば、別のコンピュータ６０２上で実行される）クライアントアプリケーションからネットワーク６３０を介してリクエストを受信し、適切なソフトウェアアプリケーションを使用して受信したリクエストを処理することによって、受信したリクエストに応答することができる。さらに、リクエストはまた、内部ユーザ（例えば、コマンドコンソールから）、外部又はサードパーティ、他の自動化アプリケーション、ならびに任意の他の適切なエンティティ、個人、システム、又はコンピュータからコンピュータ６０２に送信され得る。
コンピュータ６０２の各コンポーネントは、システムバス６０３を使用して通信することができる。いくつかの実施では、コンピュータ６０２のコンポーネントのいずれか又は全ては、ハードウェア又はソフトウェアの両方（又はハードウェア及びソフトウェアの組合せ）を含み、アプリケーションプログラミングインタフェース（ＡＰＩ）６１２又はサービスレイヤ６１３（又はＡＰＩ６１２及びサービスレイヤ６１３の組合せ）を使用して、システムバス６０３を介して、互いに又はインタフェース６０４（又は両方の組合せ）とインタフェース接続することができる。ＡＰＩ６１２は、ルーチン、データ構造、及びオブジェクトクラスの仕様を含むことができる。ＡＰＩ６１２は、コンピュータ言語に依存しないものであっても、依存するものであってもよく、完全なインタフェース、単一の機能、又はＡＰＩのセットを指す場合もある。サービスレイヤ６１３は、コンピュータ６０２又はコンピュータ６０２に通信可能に結合された他の構成要素（図示されているか否かにかかわらず）にソフトウェアサービスを提供する。コンピュータ６０２の機能は、このサービスレイヤを使用するすべてのサービスコンシューマにとってアクセス可能であってもよい。サービスレイヤ６１３によって提供されるようなソフトウェアサービスは、定義されたインタフェースを介して再使用可能な定義された機能を提供する。例えば、インタフェースは、ＪＡＶＡ（登録商標）、Ｃ＋＋、又は拡張マークアップ形式（ＸＭＬ）フォーマット又はフォーマットの組み合わせでデータを提供するコンピューティング言語の組み合わせで記述されたソフトウェアであってもよい。コンピュータ６０２の統合コンポーネントとして示されているが、代替の実施は、コンピュータ６０２の他のコンポーネント、又はコンピュータ６０２に通信可能に結合されている他のコンポーネント（図示されているかどうかにかかわらず）に関して、スタンドアロンコンポーネントとしてＡＰＩ６１２又はサービスレイヤ６１３を示すことができる。さらに、ＡＰＩ６１２又はサービスレイヤ６１３の任意の又は全ての部分は、本開示の範囲から逸脱することなく、別のソフトウェアモジュール、企業アプリケーション、又はハードウェアモジュールの子モジュール又はサブモジュールとして実施され得る。
コンピュータ６０２は、インタフェース６０４を含む。図６では単一のインタフェース６０４として示されているが、コンピュータ６０２の特定のニーズ、要望、又は特定の実施に従って、２つ以上のインタフェース６０４を使用することができる。インタフェース６０４は、分散環境においてネットワーク６３０（図示されているか否かにかかわらず）に接続されている他のシステムと通信するために、コンピュータ６０２によって使用される。一般に、インタフェース６０４はソフトウェア又はハードウェア（又はソフトウェアとハードウェアの組合せ）でエンコードされたロジックを含み、ネットワーク６３０と通信するように動作可能である。より具体的には、インタフェース６０４がネットワーク６３０又はインタフェースのハードウェアが図示のコンピュータ６０２の内部及び外部で物理信号を通信するよう動作可能であるように、通信に関連する１又は複数の通信プロトコルをサポートするソフトウェアを備えることができる。
コンピュータ６０２は、プロセッサ６０５を含む。図６では単一のプロセッサ６０５として示されているが、コンピュータ６０２の特定のニーズ、要望、又は特定の実施に従って、２つ以上のプロセッサを使用することができる。一般に、プロセッサ６０５は命令を実行し、データを操作して、コンピュータ６０２の動作、及び本開示で説明する任意のアルゴリズム、方法、機能、プロセス、フロー、及び手順を実行する。
コンピュータ６０２はまた、ネットワーク６３０（図示されているか否かにかかわらず）に接続することができるコンピュータ６０２又は他の構成要素（又は両方の組合せ）のためのデータを保持することができるデータベース６０６を含む。例えば、データベース６０６は、本開示と一致するデータを格納するインメモリ又は従来のデータベースであり得る。いくつかの実施では、データベース６０６がコンピュータ６０２及び説明した機能の特定のニーズ、要望、又は特定の実施に従って、２つ以上の異なるデータベースタイプ（例えば、ハイブリッドインメモリ及び従来のデータベース）の組合せとすることができる。図６では単一のデータベース６０６として示されているが、コンピュータ６０２及び説明された機能の特定のニーズ、要望、又は特定の実施に従って、（同じ又はタイプの組み合わせの）２つ以上のデータベースを使用することができる。データベース６０６はコンピュータ６０２の不可欠なコンポーネントとして示されているが、代替の実施では、データベース６０６がコンピュータ６０２の外部にあってもよい。例えば、データベース６０６は、地震データを保持することができる。
コンピュータ６０２はまた、コンピュータ６０２又はネットワーク６３０（図示されているか否かにかかわらず）に接続することができる他のコンポーネント（又は両方の組合せ）のためのデータを保持することができるメモリ６０７を含む。例えば、メモリ６０７は、本開示と一致するデータを記憶するランダムアクセスメモリ（ＲＡＭ）、リードオンリメモリ（ＲＯＭ）、光学的なもの、磁気的なもの等であり得る。いくつかの実施では、メモリ６０７がコンピュータ６０２及び説明した機能の特定のニーズ、要望、又は特定の実施に従って、２つ以上の異なるタイプのメモリの組合せ（たとえば、ＲＡＭと磁気記憶装置の組合せ）とすることができる。図６では単一のメモリ６０７として示されているが、コンピュータ６０２及び説明した機能の特定のニーズ、要望、又は特定の実施に応じて、２つ以上のメモリ６０７（同じ又はタイプの組合せ）を使用することができる。メモリ６０７はコンピュータ６０２の不可欠なコンポーネントとして示されているが、代替の実施ではメモリ６０７がコンピュータ６０２の外部にあってもよい。
アプリケーション６０８は、特に本開示で説明する機能に関して、コンピュータ６０２の特定のニーズ、要望、又は特定の実施に従って機能を提供するアルゴリズムソフトウェアエンジンである。例えば、アプリケーション６０８は、１又は複数のコンポーネント、モジュール、又はアプリケーションとして機能することができる。さらに、単一のアプリケーション６０８として示されているが、アプリケーション６０８はコンピュータ６０２上の複数のアプリケーション６０８として実施されてもよい。さらに、コンピュータ６０２に不可欠なものとして図示されているが、代替の実施では、アプリケーション６０８はコンピュータ６０２の外部にあってもよい。
コンピュータ６０２を含むコンピュータ・システムに関連するか、又はその外部にある任意の数のコンピュータ６０２が存在し、各コンピュータ６０２はネットワーク６３０を介して通信する。さらに、用語「クライアント」、「ユーザ」、及び他の適切な用語は、本開示の範囲から逸脱することなく、必要に応じて交換可能に使用され得る。さらに、本開示は、複数のユーザが１つのコンピュータ６０２を使用できること、又は１人のユーザが複数のコンピュータ６０２を使用できることを想定している。
主題の説明された実施は、単独で又は組み合わせて、１又は複数の特徴を含むことができる。
例えば、第１の実施は、一の領域に関連する地震データを受信することを含む方法である。この領域は複数のポイントによって表される１又は複数の地下層を含み、各ポイントは、地震速度に関連付けられる。複数のポイントにおける地震速度は、複数の勾配値に基づいて地震速度を反復的に更新することによって特定され、各勾配値はポイントに対応し、ポイントの位置における目的関数の勾配を評価することによって特定される。特定された地震速度に基づいて、１つ以上の地下層の地震画像が表示される。
上述の及び他の説明された実施はそれぞれ、オプションで、以下に示す機能の１つ又は複数を含むことができる。
以下に示す機能のいずれかと組み合わせ可能な第１の特徴として、本方法は、ＮＭＯ速度に基づいて各ポイントにおける初期の地震速度を特定することをさらに含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第２の特徴として、地震速度を反復的に更新することは、最終の反復により特定された複数のポイントにおける地震速度に基づいて、反復で使用される合成地震データを特定することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第３の特徴として、目的関数は、下記式Ａに示すものであり、
Figure 0007105809000022
 ここで、下記式Ｂが成り立ち、
Figure 0007105809000023
 ここで、ｄ（ｘ，ω；ｘｓ）は震源位置ｘｓ及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘｓ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第４の特徴として、地震速度を反復的に更新することは、目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第５の特徴として、地震速度は、最急降下法を使用して反復的に更新される。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第６の特徴として、目的関数の勾配は、ポイントの位置に関する地震データの振幅の変化とは無関係である。
第２の実施は、コンピュータメモリと、コンピュータメモリと相互運用可能に結合された１又は複数のハードウェアプロセッサとを含むシステムである。１又は複数のハードウェアプロセッサは、一の領域に関連する地震データを受信する動作を実行するように構成されている。領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる。複数の勾配値に基づいて地震速度を反復的に更新することにより複数のポイントにおける地震速度が特定される。複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定される。特定された地震速度に基づいて、１又は複数の地下層の地震画像が表示される。
上述の及び他の説明された実施はそれぞれ、オプションで、以下に示す機能の１つ又は複数を含むことができる。
以下に示す機能のいずれかと組み合わせ可能な第１の特徴として、動作は、ＮＭＯ速度に基づいて各ポイントの初期の地震速度を特定することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第２の特徴として、地震速度を反復的に更新することは、最終の反復により特定された複数のポイントにおける地震速度に基づいて、反復で使用される合成地震データを特定することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第３の特徴として、目的関数は、下記式Ａに示すものであり、
Figure 0007105809000024
 ここで、下記式Ｂが成り立ち、
Figure 0007105809000025
 ここで、ｄ（ｘ，ω；ｘｓ）は震源位置ｘｓ及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘｓ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第４の特徴として、地震速度を反復的に更新することは、目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第５の特徴として、地震速度は、最急降下法を使用して反復的に更新される。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第６の特徴として、目的関数の勾配は、ポイントの位置に関する地震データの振幅の変化とは無関係である。
第３の実施は、一の領域に関連する地震データを受信する動作を含む動作を実行するためにコンピュータシステムによって実行可能な１又は複数の命令を格納する非一時的なコンピュータ読取可能媒体である。この領域は、複数のポイントで表される１又は複数の地下層を備え、各ポイントは、地震速度に関連付けられる。複数の勾配値に基づいて地震速度を反復的に更新することにより複数のポイントにおける地震速度が特定され、前記複数の勾配値は、一のポイントに対応し、且つポイントの位置で目的関数の勾配を評価することにより特定される。特定された地震速度に基づいて、１又は複数の地下層の地震画像が表示される。
上述の及び他の説明された実施はそれぞれ、オプションで、以下に示す機能の１つ又は複数を含むことができる。
以下に示す機能のいずれかと組み合わせ可能な第１の特徴として、上記動作は、ＮＭＯ速度に基づいて各ポイントの初期の地震速度を特定することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第２の特徴として、地震速度を反復的に更新することは、最終の反復により特定された複数のポイントにおける地震速度に基づいて、反復で使用される合成地震データを特定することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第３の特徴として、目的関数は、下記式Ａに示すものであり、
Figure 0007105809000026
 ここで、下記式Ｂが成り立ち、
Figure 0007105809000027
 ここで、ｄ（ｘ，ω；ｘｓ）は震源位置ｘｓ及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘｓ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第４の特徴として、地震速度を反復的に更新することは、目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に含む。
前述の又は以下に示す機能のいずれかと組み合わせ可能な第５の特徴として、目的関数の勾配は、ポイントの位置に関する地震データの振幅の変化とは無関係である。
本明細書で記述された主題及び機能的な操作の実施は、デジタル電子回路、有形に具体化されたコンピュータのソフトウェア又はファームウェア、コンピュータハードウェアに、又はそれらの一又は複数の組み合わせで実施でき、これらデジタル電子回路、ソフトウェア、ファームウェア、及びコンピュータハードウェアは、本明細書に開示される構造及びそれらの構造的な同等物を含む。本明細書で記述された主題に係る実施は、１又は複数のコンピュータプログラム、つまり、コンピュータプログラム命令の一又は複数のモジュールとして実施されることができ、コンピュータプログラム命令は、データ処理装置により実施されるシステムによる実行のために、有形の非一時的なコンピュータ読取可能なコンピュータ記憶媒体上にエンコードされ、又はコンピュータ若しくはコンピュータにより実施されるシステムの動作を制御する。代替的に又は追加的に、プログラム命令は、人工的に生成された伝播信号、例えば、マシン生成の電気的、光学的、又は電磁気的な信号にエンコードされ、この信号は、データ処理装置による実行のために受振器装置への送信用の情報をエンコードするために生成される。コンピュータ記憶媒体は、機械により読取可能な記憶装置、機械により読取可能な記憶基板、ランダム若しくはシリアルアクセスのメモリデバイス、又はコンピュータ記憶媒体の組み合わせであり得る。
「リアル－タイム」、「リアルタイム」、「リアル（ファースト）タイム（ＲＦＴ）」、「ほぼリアルタイム（ＮＲＴ）」、「準リアルタイム」、又は（当業者が理解する）類似の用語は、個々が、実質的に同時に動作及び応答が生じることを知覚するように、該動作及び応答が時間的に近接していることを意味する。例えば、データにアクセスするための個々の動作に続くデータの表示（又は表示の開始）に対する応答の時間差は、１ｍｓ未満、１ｓｅｃ未満、又は５ｓｅｃ未満であり得る。要求されたデータが、即時に表示（又は表示開始）される必要はないが、記述されたコンピューティングシステムの処理制限と、データを例えば収集し、正確に評価し、分析し、処理し、格納し、又は送信するために必要な時間とを考慮して、意図的な遅延なしに表示（又は表示が開始）される。
用語「データ処理装置」、「コンピュータ」、又は「電子コンピュータデバイス」（又は当業者によって理解される同等の用語）は、データ処理ハードウェアを指し、あらゆる種類の装置、デバイス、及びマシンを包含し、これらの種類は、データを処理するためのものであり、例示として、プログラマブルプロセッサ、コンピュータ、又は複数のプロセッサ若しくはコンピュータを含む。コンピュータは、例えば、中央処理装置（ＣＰＵ）、ＦＰＧＡ（フィールドプログラマブルゲートアレイ）、又はＡＳＩＣ（特定用途向け集積回路）を包含する専用ロジック回路であるか、又はこれらを更に含むことができる。いくつかの実施では、データ処理装置、又は専用ロジック回路（又はデータ処理装置と専用ロジック回路との組み合わせ）は、ハードウェア系又はソフトウェア系（又はソフトウェア系及びハードウェア系の両方の基づく組み合わせ）であることができる。装置は、オプションとして、コンピュータプログラムの実行環境を生成するコード、例えば、プロセッサファームウェア、プロトコルスタック、データベース管理システム、オペレーティングシステム、又は実行環境の組み合わせを構成するコードを含むことができる。本願は、例えば、ＬＩＮＵＸ、ＵＮＩＸ（登録商標）、ＷＩＮＤＯＷＳ（登録商標）、ＭＡＣＯＳ、ＡＮＤＲＯＩＤ（登録商標）、ＩＯＳ、又はオペレーティングシステムの組み合わせといった従来のオペレーティングシステムを伴う又は伴わないデータ処理装置を使用することを意図している。
コンピュータプログラムは、プログラム、ソフトウェア、ソフトウェアアプリケーション、モジュール、ソフトウェアモジュール、スクリプト、又はコードとして参照され又は記述されるものであって、コンパイル言語若しくはインタープリタ言語を含むプログラミング言語、又は宣言型若しくは手続き型の言語の形式で記載できる。コンピュータプログラムは、コンピューティング環境で使用するための任意の形式で展開でき、スタンドアロンプログラム、モジュール、コンポーネント、サブルーチン、又は他の適切なユニットとして含む。コンピュータプログラムは、ファイルシステム内のファイルに対応できるが、必ずしも対応する必要はない。プログラムは、他のプログラム又はデータを保持するファイルの一部、例えばマークアップ言語ドキュメントに格納された一又は複数のスクリプトに保存でき、該他のプログラムは、問題のプログラム専用の単一ファイル内に、或いは複数の連携したファイル、例えば一又は複数のモジュール、サブプログラム、又はコードの一部に格納できる。コンピュータプログラムは、一又は複数のコンピュータ上に展開でき、該コンピュータは、一のサイトに位置し、又は複数のサイトに分散されて配置され、これらは通信ネットワークによって相互接続される。様々な図面に例示されたプログラムの部分は、個々のモジュールとして示され、様々なオブジェクト、方法、又は他のプロセスを介して、種々の特徴及び機能を実施しているが、プログラムは、それらの代わりに、いくつかのサブモジュール、サードパーティのサービス、コンポーネント、ライブラリ、及びそれらを、必要に応じて含むことができる。逆に、必要に応じて、様々なコンポーネントの機能及び特徴は、単一のコンポーネントに組み合わせることができる。計算上の特定を為すために使用されたしきい値は、静的に、動的に、又は静的と動的との両方で特定できる。
本明細書で説示された方法、プロセス、又は論理フローは、一又は複数のプログラム可能なコンピュータによって実行でき、該コンピュータは、入力データを操作して出力データを生成することによって機能を実行する一又は複数のコンピュータプログラムを実行する。方法、プロセス、又は論理フローは、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実行でき、また装置も、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実施できる。
コンピュータプログラムの実行のための適切なコンピュータは、汎用又は特定用途のマイクロプロセッサに基づくことができる。一般的には、ＣＰＵは、読み取り専用メモリ（ＲＯＭ）又はランダムアクセスメモリ（ＲＡＭ）あるいはその両方から命令及びデータを受け取る。コンピュータの重要な要素は、命令を行い又は実行するためのＣＰＵ、及び命令及びデータを保存するための一又は複数のメモリデバイスである。一般に、コンピュータは、データを保存するための一又は複数の大容量記憶装置、例えば磁気、光磁気ディスク、又は光ディスクを含み、或いはこれらの大容量記憶装置に動作可能に結合され、データを受信し、転送し、又はこの両方を行う。しかし、コンピュータは、そのようなデバイスを持つ必要はない。さらには、コンピュータは、別のデバイス、例えば携帯電話、携帯情報端末（ＰＤＡ）、モバイルオーディオ若しくはビデオプレーヤー、ゲームコンソール、全地球測位システム（ＧＰＳ）受信機、又はポータブルストレージデバイス、例えばユニバーサルシリアルバス（ＵＳＢ）、フラッシュドライブ等に組み込むことができる。
コンピュータプログラム命令及びデータを格納するための適切な（必要に応じて一時的又は非一時的な）コンピュータ読取可能媒体は、例示として半導体メモリデバイス、磁気デバイス、及び光学メモリデバイスを含む、あらゆる形態のメモリ、媒体、及びメモリデバイスを含む。あらゆる形態のメモリデバイスは：半導体メモリデバイス（例えば、消去可能プログラマブルリードオンリーメモリ（ＥＰＲＯＭ）、電気的消去可能プログラマブル読み取り専用メモリ（ＥＥＰＲＯＭ）、及びフラッシュメモリデバイス）と磁気ディスク（例えば、内部／ハードディスク又はリムーバブルディスク）と、光磁気ディスクと、光学的メモリディスク（例えば、ＣＤ－ＲＯＭ、ＤＶＤ＋／－Ｒ、ＤＶＤ－ＲＡＭ、及びＤＶＤ－ＲＯＭディスク）と；を含む。メモリは、様々なオブジェクト又はデータを格納でき、オブジェクト又はデータは、キャッシュ、クラス、フレームワーク、アプリケーション、バックアップデータ、ジョブ、Ｗｅｂページ、Ｗｅｂページテンプレート、データベーステーブル、動的情報を格納するリポジトリ、及び他の適切な情報を包含しており、他の適切な情報は、任意のパラメータ、任意の変数、任意のアルゴリズム、任意の命令、任意のルール、任意の制約、又は任意の参照を含む。さらに、メモリは、ログ、ポリシー、セキュリティ又はアクセスのデータ、もしくはレポートファイルといった他の適切なデータを含むことができる。プロセッサ及びメモリは、特定用途ロジック回路によって補完され又は組み込むことができる。
ユーザとのインタラクションを提供するために、本明細書で記載された主題の実施は、ユーザに情報を表示するためのディスプレイデバイス、例えばＣＲＴ（陰極線管）、ＬＣＤ（液晶ディスプレイ）、ＬＥＤ（発光ダイオード）又はプラズマモニターと、ユーザがコンピュータに入力を提供するキーボード及びポインティングデバイス、例えばマウス、トラックボール、トラックパッドとを有するコンピュータ上で実施できる。入力が、タブレットコンピュータの感圧性表面、静電容量式若しくは電気式のセンシングを使用したマルチタッチスクリーンといった、タッチスクリーンを用いてコンピュータに提供することもできる。他の種類のデバイスが、ユーザとのインタラクションを提供するために使用できる。同様に、例えば、ユーザに提供されるフィードバックは、例えば、視覚フィードバック、聴覚フィードバック、触覚フィードバック、及びあらゆる形態の感覚のフィードバックであり得る。ユーザからの入力は、音響入力、音声入力、触覚入力を包含する、あらゆる形式で受信できる。さらに、コンピュータは、ユーザによって使用されるクライアントコンピューティングデバイスからのドキュメントを受信し及び該ドキュメントをデバイスに送信することによって（例えば、Ｗｅｂブラウザから受けたリクエストに応答してユーザのクライアントデバイス上においてＷｅｂブラウザにＷｅｂページを送ることにより）ユーザとインタラクションできる。
用語「グラフィカルユーザインタフェース」又は「ＧＵＩ」は、単数又は複数で使用して、一又は複数のグラフィカルユーザインタフェース及び特定のグラフィカルユーザインタフェースのディスプレイの各々を説明できる。したがって、ＧＵＩは、Ｗｅｂブラウザ、タッチスクリーン、又はコマンドラインインタフェース（ＣＬＩ）を含むが、これらに限定されることなく、任意のグラフィカルユーザインタフェースを提示することができ、コマンドラインインタフェースは、情報を処理すると共に情報の結果をユーザに効率的に提示する。一般的に、ＧＵＩは、いくつかのユーザインタフェース（ＵＩ）要素、例えばインタラクティブフィールド、プルダウンリスト、ボタンを含むことができ、これらのうちのいくつか又は全ては、Ｗｅｂブラウザに関連付けられる。これら及び他のＵＩ要素は、Ｗｅｂブラウザの機能に関連付けられまたＷｅｂブラウザの機能を表すことができる。
本明細書で記述された主題の実施は、バックエンドコンポーネントを（例えばデータサーバとして）含むコンピューティングシステム、ミドルウェアコンポーネントを（例えばアプリケーションサーバとして）含むコンピューティングシステム、フロントエンドコンポーネントを（例えばクライアントコンピュータとして）含むコンピューティングシステム、及び、一又は複数のこのようなバックエンド、ミドルウェア、又はフロントエンドコンポーネント（例えばクライアントコンピュータは、ユーザが本明細書においで記述される主題の実施とインタラクションできるグラフィカルユーザインタフェース又はＷｅｂブラウザを有する）の任意の組み合わせを含むコンピューティングシステムにおいて実施できる。システムのコンポーネントは、有線又は無線のデジタルデータ通信（又はデータ通信の組み合わせ）、例えば通信ネットワークの任意の形式又は媒体によって相互接続されることができる。通信ネットワークの例示は、ローカルエリアネットワーク（ＬＡＮ）、無線アクセスネットワーク（ＲＡＮ）、メトロポリタンエリアネットワーク（ＭＡＮ）、広域ネットワーク（ＷＡＮ）、マイクロ波アクセスの世界的な相互運用性（ＷＩＭＡＸ）、例えば８０２．１１ａ／ｂ／ｇ／ｎ又は８０２．２０（又は８０２．１１ｘと８０２．２０の組み合わせ、又は本開示に整合する他のプロトコル）を用いる無線ローカルエリア（ＷＬＡＮ）、インターネットの全て若しくは一部、別の通信システム、又は通信ネットワークとの組み合わせを含む。ネットワークは、ネットワークアドレス間、例えば、インターネットプロトコル（ＩＰ）パケット、フレームリレーフレーム、非同期転送モード（ＡＴＭ）セル、音声、又はビデオとデータ通信することができる。
コンピューティングシステムは、クライアント及びサーバを含むことができる。クライアント及びサーバは、一般的には、互いに離れており、また典型的には、通信ネットワークを介してインタラクションする。クライアント及びサーバの関係は、それぞれのコンピュータ上で動作すると共に互にクライアント・サーバ関係にあるコンピュータプログラムのおかげで生じる。
本明細書は多くの特定の実施の詳細を含む一方で、これらは、請求される得る範囲の制限として解釈されるべきではなく、むしろ、特定のコンセプトの特定の実施固有の特徴の説明として解釈されるべきである。個別の実施の観点で本明細書において説明された機能は、組み合わせて又は単一の実施で、実現できる。逆に、単一の実施の観点で既述された様々な特徴は、複数の実施で、個別に、又は任意のサブコンビネーションで実施されることもできる。さらには、既述の特徴は特定の組み合わせで動作するものとして説示され、最初はそのようなものとして請求されているが、一又は複数の特徴は、請求された組み合わせから場合によっては削除され、またサブコンビネーション又はサブコンビネーションの変形であり得る。
主題の特定の実施が説明された。記載された実施の置換、変更、及び他の実施は、当業者には明らかである以下の請求の範囲内である。操作は特許請求の範囲に又は特定の順序で図面に描かれている一方で、これは、望ましい結果を達成するために、示された特定の順序で又は順番でそのような操作が実行されること又は全ての図示された操作が実行される（いくつかの操作は随意的と見なされる）ことを要求するものとして理解されるべきではない。ある状況では、マルチタスク又は並列処理（又はマルチタスク及び並列処理の組み合わせ）を行うことが有利であり、適切と思われる場合に実行される。
さらに、既述の実施における様々なシステムモジュール及びコンポーネントの分離又は統合が、全ての実施においてそのような分離又は統合を必要とするものとして理解されるべきではない。また、記載されたプログラムコンポーネント及びシステムは、一般的には、単一のソフトウェア製品に統合でき、又は複数のソフトウェア製品にパッケージ化できることが理解されるべきである。
これに従って、既述の例示的な実施は、本願を定義又は制約しない。他の変更、置換、及び変更も、本開示の範囲及び精神から逸脱することなく可能である。
さらには、請求された任意の実施は、少なくとも、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータシステムに適用可能であると考えられ、非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータにより読取可能な命令を格納し、またコンピュータシステムは、ハードウェアプロセッサに相互に動作可能に結合されたコンピュータメモリを含み、ハードウェアプロセッサは、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法又は命令を実行するように構成される。
以下、本発明の実施の態様の例を列挙する。
［第１の局面］
一の領域に関連する地震データを受信するステップであって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、ステップと；
複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定するステップであって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定される、ステップと；
前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示するステップと；を備える、
方法。
［第２の局面］
ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定するステップを更に備える、
第１の局面に記載の方法。
［第３の局面］
前記地震速度を反復的に更新することは、最終の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定することを更に備える、
第１の局面に記載の方法。
［第４の局面］
前記目的関数は、下記式Ａに示すものであり、
Figure 0007105809000028
ここで、下記式Ｂが成り立ち、
Figure 0007105809000029
ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す、
第３の局面に記載の方法。
［第５の局面］
前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
第４の局面に記載の方法。
［第６の局面］
前記地震速度は、最急降下法を使用して反復的に更新される、
第１の局面に記載の方法。
［第７の局面］
前記目的関数の勾配は、前記ポイントの位置に関する前記地震データの振幅の変化には依存しない、
第１の局面に記載の方法。
［第８の局面］
コンピュータメモリと；
前記コンピュータメモリと相互運用可能に結合され、所定の動作を実行するように構成された１又は複数のプロセッサであって、前記動作は：
一の領域に関連する地震データを受信する処理であって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、処理と；
複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定する処理であって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定される、処理と；
前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示する処理と；を備える、前記１又は複数のプロセッサと；を備える、
システム。
［第９の局面］
前記動作は、ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定する処理を更に備える、
第８の局面に記載のシステム。
［第１０の局面］
前記地震速度を反復的に更新することは、最終の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定することを更に備える、
第８の局面に記載のシステム。
［第１１の局面］
前記目的関数は、下記式Ａに示すものであり、
Figure 0007105809000030
ここで、下記式Ｂが成り立ち、
Figure 0007105809000031
ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す、
第１０の局面に記載のシステム。
［第１２の局面］
前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
第１１の局面に記載のシステム。
［第１３の局面］
前記地震速度は、最急降下法を使用して反復的に更新される、
第８の局面に記載のシステム。
［第１４の局面］
前記目的関数の勾配は、前記ポイントの位置に関する前記地震データの振幅の変化には依存しない、
第８の局面に記載のシステム。
［第１５の局面］
所定の動作を実行するためにコンピュータシステムによって実行可能な１又は複数の命令を格納する非一時的なコンピュータ読取可能媒体であって、前記動作は：
一の領域に関連する地震データを受信する処理であって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、処理と；
複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定する処理であって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定される、処理と；
前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示する処理と；を備える、
非一時的なコンピュータ読取可能媒体。
［第１６の局面］
前記動作は、ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定する処理を更に備える、
第１５の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１７の局面］
前記地震速度を反復的に更新することは、最終の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定することを更に備える、
第１５の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１８の局面］
前記目的関数は、下記式Ａに示すものであり、
Figure 0007105809000032
ここで、下記式Ｂが成り立ち、
Figure 0007105809000033
ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表す、
第１７の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１９の局面］
前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
第１８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第２０の局面］
前記目的関数の勾配は、前記ポイントの位置に関する前記地震データの振幅の変化には依存しない、
第１５の局面に記載の非一時的なコンピュータ読取可能媒体。
６００ コンピュータシステム
６０２ コンピュータ
６０５ プロセッサ
６０７ メモリ

Claims (14)

  1. 一の領域に関連する地震データを受信するステップであって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、ステップと；
    複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定するステップと、
    直前の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定するステップであって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定され、
    前記目的関数は、下記式Ａに示すものであり、
    Figure 0007105809000034
    ここで、下記式Ｂが成り立ち、
    Figure 0007105809000035
    ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表し、前記目的関数は、前記地震速度と前記合成地震データとの間の相互相関に基づく、ステップと；
    前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示するステップと；を備える、
    方法。
  2. ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定するステップを更に備える、
    請求項１に記載の方法。
  3. 前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
    請求項１に記載の方法。
  4. 前記地震速度は、最急降下法を使用して反復的に更新される、
    請求項１に記載の方法。
  5. 前記目的関数の勾配は、前記ポイントの位置に関する前記地震データの振幅の変化には依存しない、
    請求項１に記載の方法。
  6. コンピュータメモリと；
    前記コンピュータメモリと相互運用可能に結合され、所定の動作を実行するように構成された１又は複数のプロセッサであって、前記動作は、
    一の領域に関連する地震データを受信する処理であって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、処理と；
    複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定する処理と、
    直前の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定する処理であって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定され、
    前記目的関数は、下記式Ａに示すものであり、
    Figure 0007105809000036
    ここで、下記式Ｂが成り立ち、
    Figure 0007105809000037
    ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表し、前記目的関数は、前記地震速度と前記合成地震データとの間の相互相関に基づく、処理と；
    前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示する処理と；を備える、前記１又は複数のプロセッサと；を備える、
    システム。
  7. 前記動作は、ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定する処理を更に備える、
    請求項６に記載のシステム。
  8. 前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
    請求項６に記載のシステム。
  9. 前記地震速度は、最急降下法を使用して反復的に更新される、
    請求項６に記載のシステム。
  10. 前記目的関数の勾配は、前記ポイントの位置に関する前記地震データの振幅の変化には依存しない、
    請求項６に記載のシステム。
  11. 所定の動作を実行するためにコンピュータシステムによって実行可能な１又は複数の命令を格納する非一時的なコンピュータ読取可能媒体であって、前記動作は、
    一の領域に関連する地震データを受信する処理であって、前記領域は、複数のポイントで表される１又は複数の地下層を備え、前記複数のポイントの各ポイントは、地震速度に関連付けられる、処理と；
    複数の勾配値に基づいて前記地震速度を反復的に更新することにより前記複数のポイントにおける前記地震速度を特定する処理と、
    直前の反復により特定された前記複数のポイントにおける前記地震速度に基づいて、前記反復で使用される合成地震データを特定する処理であって、前記複数の勾配値は、一のポイントに対応し、且つ当該ポイントの位置で目的関数の勾配を評価することにより特定され、
    前記目的関数は、下記式Ａに示すものであり、
    Figure 0007105809000038
    ここで、下記式Ｂが成り立ち、
    Figure 0007105809000039
    ここで、ｄ（ｘ，ω；ｘ ｓ ）は震源位置ｘ ｓ 及び受振器位置ｘにおいて受信した地震データを表し、ｐ（ｘ，ω；ｘ ｓ ）は合成地震データを表し、＊は共役オペレーションを表し、ωは角周波数を表し、前記目的関数は、前記地震速度と前記合成地震データとの間の相互相関に基づく、処理と；
    前記特定された地震速度に基づいて、前記１又は複数の地下層の地震画像を表示する処理と；を備える、
    非一時的なコンピュータ読取可能媒体。
  12. 前記動作は、ノーマルムーブアウト（ＮＭＯ）速度に基づいて前記各ポイントの初期の地震速度を特定する処理を更に備える、
    請求項１１に記載の非一時的なコンピュータ読取可能媒体。
  13. 前記地震速度を反復的に更新することは、前記目的関数の値が所定の閾値よりも小さい場合に、反復を停止することを更に備える、
    請求項１１に記載の非一時的なコンピュータ読取可能媒体。
  14. 前記地震速度は、最急降下法を使用して反復的に更新される、
    請求項１１に記載の非一時的なコンピュータ読取可能媒体。
JP2019564802A 2017-05-22 2018-05-22 周波数領域における地震速度反転のための振幅に依存しない勾配の計算 Active JP7105809B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201762509300P 2017-05-22 2017-05-22
US62/509,300 2017-05-22
US15/956,261 2018-04-18
US15/956,261 US11269097B2 (en) 2017-05-22 2018-04-18 Computing amplitude independent gradient for seismic velocity inversion in a frequency domain
PCT/US2018/033828 WO2018217706A1 (en) 2017-05-22 2018-05-22 Computing amplitude independent gradient for seismic velocity inversion in a frequency domain

Publications (3)

Publication Number Publication Date
JP2020521140A JP2020521140A (ja) 2020-07-16
JP2020521140A5 JP2020521140A5 (ja) 2021-06-17
JP7105809B2 true JP7105809B2 (ja) 2022-07-25

Family

ID=64271536

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019564802A Active JP7105809B2 (ja) 2017-05-22 2018-05-22 周波数領域における地震速度反転のための振幅に依存しない勾配の計算

Country Status (7)

Country Link
US (1) US11269097B2 (ja)
EP (1) EP3631522A1 (ja)
JP (1) JP7105809B2 (ja)
CN (1) CN110914718B (ja)
CA (1) CA3064293A1 (ja)
SA (1) SA519410566B1 (ja)
WO (1) WO2018217706A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11048001B2 (en) * 2018-03-30 2021-06-29 Cgg Services Sas Methods using travel-time full waveform inversion for imaging subsurface formations with salt bodies
CN111461988B (zh) * 2020-04-02 2022-09-23 清华大学 一种基于多任务学习的地震速度模型超分辨率技术
US12044124B2 (en) 2021-02-05 2024-07-23 Saudi Arabian Oil Company Method and system for real-time hole cleaning using a graphical user interface and user selections
US20220260746A1 (en) * 2021-02-12 2022-08-18 Saudi Arabian Oil Company Method and system for determining permeability and lost circulation
CN113011006B (zh) * 2021-02-25 2021-10-22 中国科学院声学研究所 一种基于互相关函数脉冲波形匹配的目标深度估计方法
US11768303B2 (en) * 2021-04-22 2023-09-26 Saudi Arabian Oil Company Automatic data enhancement for full waveform inversion in the midpoint-offset domain
US11971513B2 (en) 2021-05-21 2024-04-30 Saudi Arabian Oil Company System and method for forming a seismic velocity model and imaging a subterranean region
US11867857B2 (en) * 2021-07-13 2024-01-09 Saudi Arabian Oil Company Method and system for updating a seismic velocity model
CN113777651B (zh) * 2021-08-16 2024-05-07 华能西藏雅鲁藏布江水电开发投资有限公司 人工地震动合成方法、装置、电子设备及存储介质
US12037857B2 (en) 2021-11-30 2024-07-16 Saudi Arabian Oil Company Method and system for determining hole cleaning efficiency based on wellbore segment lengths

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010531453A (ja) 2007-06-26 2010-09-24 シン，チャンス 地下構造の映像化のためにラプラスドメインでの波形逆算を利用した速度分析方法
US20120314538A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2569246B2 (ja) * 1992-03-02 1997-01-08 戸田建設株式会社 切羽前方探査方式
US5570321A (en) 1994-03-03 1996-10-29 Atlantic Richfield Company Seismic velocity model optimization method using simulated annearling to determine prestack travel-times
US5583825A (en) * 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
MY129095A (en) 2001-02-13 2007-03-30 Exxonmobil Upstream Res Co Method for spectral balancing of near-and far-offset seismic data.
US7072767B2 (en) * 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
KR100966904B1 (ko) 2008-05-06 2010-06-30 재단법인서울대학교산학협력재단 라플라스-푸리에 영역의 파형 역산을 이용한 지하 구조의 영상화 장치, 방법 및 기록매체
CN101630018B (zh) 2008-07-16 2011-12-07 中国石油天然气股份有限公司 一种控制全声波方程反演过程的地震勘探数据处理方法
US20100142316A1 (en) 2008-12-07 2010-06-10 Henk Keers Using waveform inversion to determine properties of a subsurface medium
US10288751B2 (en) * 2009-10-08 2019-05-14 Westerngeco L.L.C. Migration velocity analysis using seismic data
US9244181B2 (en) 2009-10-19 2016-01-26 Westerngeco L.L.C. Full-waveform inversion in the traveltime domain
US20120316791A1 (en) * 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion by non-linear model update
PT2732312T (pt) 2011-07-12 2021-06-23 Eni Spa Análise de velocidade de migração de equação de onda usando distorção de imagem
US10871584B2 (en) * 2013-06-17 2020-12-22 Westerngeco L.L.C. Seismic data processing
EP3069171B1 (en) * 2013-11-12 2021-08-25 Westerngeco LLC Seismic data processing
US20150293246A1 (en) * 2014-04-09 2015-10-15 Thomas A. Dickens Frequency-domain augmented time-domain full wavefield inversion
US9470811B2 (en) * 2014-11-12 2016-10-18 Chevron U.S.A. Inc. Creating a high resolution velocity model using seismic tomography and impedance inversion
WO2016154404A1 (en) * 2015-03-26 2016-09-29 Schlumberger Technology Corporation Seismic waveform inversion
GB2538807B (en) * 2015-05-29 2019-05-15 Sub Salt Solutions Ltd Method for improved geophysical investigation
GB2538804A (en) 2015-05-29 2016-11-30 Sub Salt Solutions Ltd Improved method for inversion modelling
CN105005076B (zh) 2015-06-02 2017-05-03 中国海洋石油总公司 基于最小二乘梯度更新速度模型的地震波全波形反演方法
WO2018013257A1 (en) * 2016-07-13 2018-01-18 Exxonmobil Upstream Research Company Joint full wavefield inversion of p-wave velocity and attenuation using an efficient first order optimization

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010531453A (ja) 2007-06-26 2010-09-24 シン，チャンス 地下構造の映像化のためにラプラスドメインでの波形逆算を利用した速度分析方法
US20120314538A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion

Also Published As

Publication number Publication date
JP2020521140A (ja) 2020-07-16
US11269097B2 (en) 2022-03-08
CN110914718B (zh) 2022-07-15
EP3631522A1 (en) 2020-04-08
US20180335530A1 (en) 2018-11-22
SA519410566B1 (ar) 2022-09-21
CA3064293A1 (en) 2018-11-29
CN110914718A (zh) 2020-03-24
WO2018217706A1 (en) 2018-11-29

Similar Documents

Publication Publication Date Title
JP7105809B2 (ja) 周波数領域における地震速度反転のための振幅に依存しない勾配の計算
US11635539B2 (en) Imaging shallow heterogeneities based on near-surface scattered elastic waves
EP3073296B1 (en) Full waveform inversion method for seismic data processing using preserved amplitude reverse time migration
JP7145866B2 (ja) 指向性波動場イメージングを用いた地球物理学的画像生成
JP7019713B2 (ja) 共通焦点法を用いる連続的な地震貯留層モニタリング
JP2020523589A (ja) リバースタイムマイグレーションからの角度ドメイン共通画像ギャザーのモデリング
US11275190B2 (en) Generating diffraction images based on wave equations
JP2020521126A (ja) 地震データにおけるノイズの抑制
US11340368B2 (en) Generating a velocity model and a density model of a subsurface structure of a reservoir
CN113348385A (zh) 弹性波逆时迁移中的avo成像条件
JP2020519889A (ja) 全波形音波動場分離のための処理方法
CN110799856B (zh) 使用与地下方位角和反射角相关的全波形反演产生速度模型
US11016205B2 (en) Generating target-oriented acquisition-imprint-free prestack angle gathers using common focus point operators
JP7254778B2 (ja) リバースタイムマイグレーションにおける偽像除去
US10908307B2 (en) Evaluating processing imprint on seismic signals
US20190187310A1 (en) Subsalt Imaging Tool for Interpreters
US20220317325A1 (en) System and method for estimating one-way propagation operators
US11460595B2 (en) Unified continuous seismic reservoir monitoring
US11656378B2 (en) Seismic imaging by visco-acoustic reverse time migration

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210430

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210430

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220301

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220524

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20220524

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220712

R150 Certificate of patent or registration of utility model

Ref document number: 7105809

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150