(57)【要約】連続的な地震貯留層モニタリングの方法は、一定の期間内に、貯留槽に関連付けられた複数の地震データセットを受信するステップを含み、複数の地震データセットは期間内の異なる時間に受信した地震データに対応する。貯留槽は複数の反射体を含み、各反射体は反射係数を有する。複数の地震データセットのそれぞれに対して、複数の反射体の反射係数は、複数の反射体と複数の取得表面位置とに関連付けられた共通焦点（ＣＦＰ）演算子を繰り返し更新するステップにより特定される。異なる地震データセットに対応する反射係数は、期間内の反射係数の変化を特定するために比較される。反射係数の変化が表示される。【特許請求の範囲】
【請求項１】
一定の期間内に、それぞれが反射係数を有する複数の反射体を含む貯留槽、に関連付けられた複数の地震データセットであって前記期間内の異なる時間に受信した地震データに対応する複数の地震データセットを受信するステップと；
前記複数の地震データセットのそれぞれに対して、前記複数の反射体と複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップにより、前記複数の反射体の反射係数を特定するステップと；
前記期間内の前記反射係数の変化を特定するために、異なる地震データセットに対応する前記反射係数を比較するステップと；
前記反射係数の前記変化を表示するステップと；を備える、
方法。
【請求項２】
各共通焦点演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す、
請求項１に記載の方法。
【請求項３】
前記複数の反射体と前記複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新された共通焦点演算子を
【数６】
2020511656.tif 000007
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間シフト取得に基づいて特定された残余伝播時間である、
請求項１に記載の方法。
【請求項４】
前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期共通焦点演算子を特定するステップをさらに備える、
請求項３に記載の方法。
【請求項５】
前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える、
請求項１に記載の方法。
【請求項６】
前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える、
請求項１に記載の方法。
【請求項７】
前記貯留層が炭化水素液を備える、
請求項１に記載の方法。
【請求項８】
コンピュータメモリと；
前記コンピュータメモリと相互運用可能に接続された１つ以上のハードウェアプロセッサとを備え；
前記ハードウェアプロセッサは、
一定の期間内に、それぞれが反射係数を有する複数の反射体を含む貯留槽、に関連付けられた複数の地震データセットであって前記期間内の異なる時間に受信した地震データに対応する複数の地震データセットを受信するステップと；
前記複数の地震データセットのそれぞれに対して、前記複数の反射体と複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップにより、前記複数の反射体の反射係数を特定するステップと；
前記期間内の前記反射係数の変化を特定するために、異なる地震データセットに対応する前記反射係数を比較するステップと；
前記反射係数の前記変化を表示するステップと；を備える操作を実行する、
システム。
【請求項９】
各共通焦点演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す、
請求項８に記載のシステム。
【請求項１０】
前記複数の反射体と前記複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新された共通焦点演算子を
【数７】
2020511656.tif 000008
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間シフト取得に基づいて特定された残余伝播時間である、
請求項８に記載のシステム。
【請求項１１】
前記操作は、前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期共通焦点演算子を特定するステップをさらに備える、
請求項１０に記載のシステム。
【請求項１２】
前記操作は、前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える、
請求項８に記載のシステム。
【請求項１３】
前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える、
請求項８に記載のシステム。
【請求項１４】
前記貯留層が炭化水素液を備える、
請求項８に記載のシステム。
【請求項１５】
コンピュータシステムによって実行可能な１つ以上の指令が格納された非一時的なコンピュータ読み取り可能な記録媒体であって；
一定の期間内に、それぞれが反射係数を有する複数の反射体を含む貯留槽、に関連付けられた複数の地震データセットであって前記期間内の異なる時間に受信した地震データに対応する複数の地震データセットを受信するステップと；
前記複数の地震データセットのそれぞれに対して、前記複数の反射体と複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップにより、前記複数の反射体の反射係数を特定するステップと；
前記期間内の前記反射係数の変化を特定するために、異なる地震データセットに対応する前記反射係数を比較するステップと；
前記反射係数の前記変化を表示するステップと；を備える操作を実行する、
非一時的なコンピュータ読み取り可能な記録媒体。
【請求項１６】
各共通焦点演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す、
請求項１５に記載の非一時的なコンピュータ読み取り可能な記録媒体。
【請求項１７】
前記複数の反射体と前記複数の取得表面位置とに関連付けられた共通焦点演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新された共通焦点演算子を
【数８】
2020511656.tif 000009
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間シフト取得に基づいて特定された残余伝播時間である、
請求項１５に記載の非一時的なコンピュータ読み取り可能な記録媒体。
【請求項１８】
前記操作は、前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期共通焦点演算子を特定するステップをさらに備える、
請求項１７に記載の非一時的なコンピュータ読み取り可能な記録媒体。
【請求項１９】
前記操作は、前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える、
請求項１５に記載の非一時的なコンピュータ読み取り可能な記録媒体。
【請求項２０】
前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える、
請求項１５に記載の非一時的なコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、２０１７年３月１６日に出願された米国特許出願第６２／４７２，４０７号に基づく優先権を主張し、当該米国特許出願のすべての記載内容を援用する。
【０００２】
本開示は、地震データ処理に関する。
【背景技術】
【０００３】
一定期間にわたる連続的な貯留層モニタリングは、取得された大量のデータセットを処理するために高速かつ正確な画像を生成することができる地震深度画像化法に関与する。連続的な貯留層モニタリングのための地震画像化の重要なステップは、直接的に観測できない地球の固体内部の波動場の推定である。表面に置かれた受信機の開境界に沿って記録された地震データに基づく標準的な推定は、画像化の前に高解像度の地震波速度モデルを利用できない限り、エネルギーが複雑な地下内をどのように伝播するかを説明するには一般的に不十分であり、さもなければ、地下の複数の散乱波（多重）を正確に予測することができない。場合によっては、十分な精度の速度モデルを作成するために、表土全体に対して厳密な深度指向の速度推定法が用いられる。
【発明の概要】
【０００４】
本開示は、共通の焦点（ＣＦＰ）方法を使用する連続的な地震貯留層モニタリングのためのコンピュータ実装方法、コンピュータプログラム製品、及びコンピュータシステムを含む方法及びシステムを説明する。
【０００５】
ある実施において、複数の地震データセットは、一定の期間内に貯留槽に関連付けられて受信され、前記複数の地震データセットは前記期間内の異なる時間に受信した地震データに対応する。前記貯留槽は複数の反射体を含み、各反射体は反射係数を有する。前記複数の地震データセットのそれぞれに対して、前記複数の反射体の反射係数は、前記複数の反射体と複数の取得表面位置とに関連付けられた共通焦点（ＣＦＰ）演算子を繰り返し更新するステップにより特定される。異なる地震データセットに対応する前記反射係数は、前記期間内の前記反射係数の変化を特定するために比較される。前記反射係数の前記変化は表示される。
【０００６】
前述の実施は、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータにより実施されるシステムを用いて実施可能である。非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータ読取可能な命令を格納する。コンピュータにより実施されるシステムは、非一時的なコンピュータ読取可能媒体に格納された命令／コンピュータにより実施される方法を実行するように構成されたハードウェアプロセッサと相互に利用可能に結合されたコンピュータメモリを含む。
【０００７】
本開示で述べる主題は、大量の地震データセットを効率的かつ効果的に処理して、地震画像を迅速かつ正確に生成することにより、四次元の連続的な地震貯留層モニタリングを可能にする。ここで述べる手法は、詳細な速度モデルを導き出さず、かつ、取得した三次元データセット全体を処理しない、堅牢な反復演算子更新手順に基づく手法である。一定期間内の生成された地震画像と、特定された反射係数の変化とを用いて、油井の掘削パラメータを特定する等、効果的な石油及びガス探査が可能になる。他の利点は、当業者にとっては明白となろう。
【０００８】
本明細書の主題の１つ以上の実施の詳細を、添付する図面及び説明に記載する。本主題の他の特徴、側面及び利点は、かかる説明及び図面並びに請求項から明らかになる。
【図面の簡単な説明】
【０００９】
【図１】ある実施による共通焦点（ＣＦＰ）法を用いて連続的な地震貯留層モニタリングを行うための例示の方法のフローチャートである。
【００１０】
【図２】ある実施による１つの反射体に対するＣＦＰ集束演算子を繰り返し更新するための例示の方法のフローチャートである。
【００１１】
【図３Ａ】ある実施による１つの目標反射体に対する初期集束演算子を用いた差分時間シフト（ＤＴＳ）取得のセットの図である。
【００１２】
【図３Ｂ】ある実施による１つの目標反射体に対する更新された集束演算子を用いたＤＴＳ取得のセットの図である。
【００１３】
【図４Ａ】ある実施による１つの目標反射体に対する初期集束演算子を示す図である。
【００１４】
【図４Ｂ】ある実施による１つの目標反射体に対する１回の反復後に更新された集束演算子を示す図である。
【００１５】
【図４Ｃ】ある実施による１つの目標反射体に対する初期集束演算子と更新された集束演算子との間の差を示す図である。
【００１６】
【図５Ａ】ある実施による目標地平線における反射係数の正規化された二乗平均平方根（ＮＲＭＳ）振幅を示す図である。
【００１７】
【図５Ｂ】ある実施による調査２と調査１７との間のＮＲＭＳ振幅差を示す図である。
【００１８】
【図６】ある実施による本開示に記載されているような記述されたアルゴリズム、方法、機能、プロセス、フロー、及び手順に関連付けられた計算機能を提供するために使用される例示のコンピュータシステムのブロック図である。
【発明を実施するための形態】
【００１９】
種々の図面における類似の参照番号及び記号は、同様の要素を指す。
【００２０】
以下の詳細な説明は、共通焦点（ＣＦＰ）法を用いた連続的な地震貯留層モニタリングについて述べ、当業者が１つ以上の特定の実施において、開示の主題を作り上げて、これを用いることができるように提示する。開示の実施に対し、様々な改変、変更、及び置き換えを行うことができ、それらは当業者には容易に明白になろう。開示の実施において定義される一般原理は、本開示の範囲から逸脱することなく、他の実施及び応用に適用されてもよい。したがって、本開示は、記述又は図示の実施に限定されるものではなく、開示の原理及び特徴に一致する最も広い範囲が付与されるべきである。
【００２１】
貯留層内の変化をモニタするために、一定期間にわたって収集された地震データを、地震画像の生成に用いることができる。例えば、貯留層の地震反射データを、一定期間にわたって月に１回取得できる。一定期間にわたる貯留層の変化、例えば炭化水素液の組成や位置の変化、を特定するための複数の三次元（３Ｄ）地震画像を生成するために、取得した複数の地震データセットを用いることができる。大まかに言えば、ここで述べる手法は、連続する３Ｄ地震データセットにＣＦＰ法を用いて、詳細な速度モデルを導き出さずに、ＣＦＰ取得と差分時間シフト（又は、差分時間偏移、差分時間ずれ）（ＤＴＳ）取得とに基づいてＣＦＰ演算子を繰り返し更新する。
【００２２】
図１は、いくつかの実施による、ＣＦＰ法を用いて連続的な地震貯留層モニタリングを行うための例示の方法１００のフローチャートである。説明を明確にするために、以下の説明は、概して、本開示の中での他図との関連において方法１００を説明する。例えば、図６で説明するコンピュータシステムにより、方法１００を実行できる。しかし、方法１００は、例えば、任意の適切なシステム、（システム）環境、ソフトウェア、及びハードウェアによって、又は、必要に応じて、システムと、（システム）環境と、ソフトウェアと、ハードウェアとの組み合わせによって実行してもよいことが分かるであろう。実施によっては、方法１００の様々な手順を、並行して、組み合わせて、ループで、又は任意の順序で、実行できる。
【００２３】
方法１００は、貯留層を識別するブロック１０２から開始する。例えば、貯留層は地球内部の１つ以上の地下層を含むことができ、貯留層の地理的境界を、使用者からの入力に基づいて識別できる。貯留層は複数の反射体を含むことができ、各反射体は反射係数を有することができる。例えば、貯留層は、長さ５キロメートル（ｋｍ）、幅５ｋｍ、深さ５ｋｍの３Ｄ領域であり得る。各反射体を長さ５メートル（ｍ）、幅５ｍ、深さ５ｍの立方体としてモデル化する場合、貯留層を１０９個の反射体で表すことができる。反射体の位置は、立方体の中心又は立方体内の他のポイントで表すことができる。実施によっては、貯留層は立方体に分割でき、立方体の角は貯留層内で反射体を表すことができる。例えば、領域が５ｍ×５ｍ×５ｍの立方体に分割される場合、貯留層を１００１３個の反射体で表せる。他の方法を用いて、貯留層を表すために反射体のセットを特定できる。
【００２４】
ブロック１０４において、一定期間内に、貯留層に関連付けられた複数の地震データセットが受信され、その複数の地震データセットは、その期間内の異なる時間に受信された地震データに対応する。例えば、貯留層の地震データを１８か月にわたり毎月１回取得でき、合計１８セットの地震データを取得できる。地震データ取得の期間中、多数の受信機（例えば、ジオフォン（受振器）やハイドロフォン（水中聴音器））を地表の上又は下に配置できる。地震源は地震波を地球内部に送り、受信機は地球内部の各地下層で反射した波を記録できる。地震源は、例えばトラックで牽引され、様々な位置で地震波を発生させることができる。例えば、地震源は、受信機が反射波を記録するために第１の位置で最初のショットを発射し、その地震源は第２の位置へ移動して第２のショットを発射する。単発ショットに対応する各受信機で記録されたデータは、トレース（出力記録）と呼ばれる。例えば、地震源が１００，０００箇所の異なる位置でショットを発射し、１，０００個の受信機がある場合、結果の地震データセットは１０８個のトレースを有することができる。
【００２５】
ブロック１０６において、各地震データセットに対する貯留層内の反射体の反射係数は、反射体及び複数の取得表面位置に関連付けられたＣＦＰ集束演算子を繰り返し更新することによって特定できる。取得表面位置は、地震源の位置と地震の受信機の位置とを含むことができる。１００，０００箇所のショット位置と１，０００個の受信機の前述の貯留層の実施例では、取得表面位置の合計は１０１，０００箇所である。
【００２６】
集束演算子の繰り返しの更新は、式（１）及び式（２）に示すように、ＣＦＰ法に基づく一方向グリーン関数を用いて、データ駆動を行うことができる。
【数１】
2020511656.tif 000002
ここで、ＣＦＰ演算子のＦ（ｉ）（ｚ０，ｚｍ）は、ｉ番目の反復における取得表面位置ｚ０と目標反射体ｚｍとの間の伝播時間又は移動時間を表し、ΔＦ（ｉ）（ｚ０，ｚｍ）は、取得表面位置ｚ０と目標反射体ｚｍとの間の移動時間をΔＴ（ｉ）／２ごとに更新し、ωは角周波数である。式（１）及び式（２）は周波数領域で表され、演算子Ｆ（ｉ）（ｚ０，ｚｍ）は、取得表面位置ｚ０と目標反射体ｚｍとの間の移動時間に相当する時間において無限の値を持ち、その他はゼロである、ディラックデルタ関数のフーリエ変換であることに注意されたい。繰り返し更新する演算は、式（３）及び式（４）により時間領域で表わすこともできる。
【数２】
2020511656.tif 000003
ここで、＊はコンボリューション（畳み込み）、δはディラックデルタ関数、ｔは初期グリーン関数の移動時間である。場合によっては、後述するように、式（２）又は式（４）のΔＴ（ｉ）は、差分時間シフト（ＤＴＳ）取得における時間ゼロから取り出された時間偏差である。
【００２７】
実施によっては、初期ＣＦＰ演算子Ｆ（０）（ｚ０，ｚｍ）は、ｚ０とｚｍとの間の移動時間の大まかな推定値であり得る。典型的な実施では、目標地平線において画成された通常の反射走（ＮＭＯ）速度を用いて、初期演算子を導き出すことができる。後述するように、ＣＦＰ取得において、目標反射体以外からの反射イベントがさらに存在する可能性があることに注意されたい。等移動時間の原理は、演算子の運動学が正しい場合、各ＣＦＰ取得における目標反射の移動時間が、それぞれの集束演算子の時間反転と同じになることを確約する。等移動時間の原理に従って、ＤＴＳ取得は、後述するように、集束演算子をトレースごとにそれぞれのＣＦＰ取得に時間関連付けることで生成できる。ＤＴＳ取得における自動取り出しに基づいて、ｉ番目の繰り返しにおける集束演算子は、目標反射体の各地下グリッドポイントに対して、等移動時間の原理が満たされるまで（すなわち、ＤＴＳ取得の目標イベントが平坦になるまで）式（１）又は式（３）を用いて繰り返し更新される。式（２）又は式（４）のΔＴ（ｉ）は、ＤＴＳ取得における時間ゼロから取り出された時間偏差である。目標指向のスタックされた画像は、最後の集束演算子からのＤＴＳ取得をスタックすることによって導き出すことができる。
【００２８】
図２は、いくつかの実施による、１つの反射体に対するＣＦＰ集束演算子を繰り返し更新するための例示の方法２００のフローチャートである。言い換えれば、方法２００は各反射体に対して実行され、１０９個の反射体でできている先の貯留層の例では、方法２００を１０９回実行することができる。説明を明確にするために、以下の説明は、概して、本開示中の他図において方法２００を説明する。例えば、図６で説明するコンピュータシステムにより方法２００を実行できる。しかし、方法２００は、システム、（システム）環境、ソフトウェア、ハードウェアによって、又は、必要に応じて、システムと、（システム）環境と、ソフトウェアと、ハードウェアの組み合わせによって実行してもよいことが分かるであろう。実施によっては、方法２００の様々なステップを、並行して、組み合わせて、ループで、又は任意の順序で実行できる。
【００２９】
方法２００は、反復カウンタｉがゼロとして初期化され、目標反射体ｚｍが選択されたブロック２０２から開始する。
【００３０】
ブロック２０４において、取得表面位置ｚ０と目標反射体ｚｍとの間の移動時間を記述する演算子Ｆ（０）（ｚ０，ｚｍ）は、各取得表面位置に対して初期化される。１０１，０００箇所の取得表面位置の先の貯留層の実施例では、初期演算子Ｆ（０）（ｚ０，ｚｍ）は１０１，０００個のｚ０とｚｍとのペア（組）に対して特定される。実施によっては、初期演算子を導き出すために、目標地平線で画成された通常の反射走（ＮＭＯ）速度を用いることができる。例えば、目標反射体の位置と取得表面位置とが判っているため、初期移動時間は、ＮＭＯ速度と、目標反射体と取得表面位置との間の距離と、に基づいて導き出せる。実施によっては、Ｆ（０）（ｚ０，ｚｍ）は、周波数ドメインにおいてｅｉωＴＮＭＯに等しくできる、又は時間ドメインにおいてδ（ｔ＋ＴＮＭＯ）に等しくできる、ここで、ＴＮＭＯはＮＭＯ速度に基づいて推定される初期移動時間である。
【００３１】
ブロック２０６において、ＣＦＰ取得はＦ（ｉ）（ｚ０，ｚｍ）を用いて生成される。最初の反復において、ＣＦＰ取得はＦ（０）（ｚ０，ｚｍ）を用いて生成される。実施によっては、地震データは取得表面位置に基づいて分類される。例えば、地震データは受信機の位置に基づいて分類される。前述の貯留層の例では、１０８個の地震データのトレースは、以下の順序で並べ替えることができる。即ち、第１の受信機で受信した１００，０００個のトレース、第２の受信機で受信した１００，０００個のトレース他。最初の受信機からの１００，０００個のそれぞれのトレースは、対応する演算子Ｆ（ｉ）（ｚ０，ｚｍ）と関連付けられ、ここで、ｚ０は特定のトレースに対応する地震源の位置である。演算子Ｆ（ｉ）（ｚ０，ｚｍ）と関連付けることにより、各トレースは演算子の移動時間に基づいてタイムシフト（又は修正）される。最初の受信機で演算子Ｆ（ｉ）（ｚ０，ｚｍ）を１００，０００個のトレースに対応させた後、結果のトレースは１つのＣＦＰトレースになるまで加算される。同様に、２番目の受信機で、演算子Ｆ（ｉ）（ｚ０，ｚｍ）を１００，０００個のトレースに対応させた後、２番目のＣＦＰトレースになるまで加算される。１，０００個の受信機で演算子をトレースに適用すると、１，０００個のＣＦＰトレースが生成される。実施によっては、例えば１，０００個のＣＦＰトレースが１つのＣＦＰ取得を形成するように、複数のＣＦＰトレースが１つのＣＦＰ取得を形成できる。実施によっては、１つのＣＦＰトレースは１つのＣＦＰ取得である。
【００３２】
ブロック２０８において、演算子Ｆ（ｉ）（ｚ０，ｚｍ）を用いてＣＦＰ取得を関連付けることによりＤＴＳ取得を生成することができる。各ＣＦＰ取得は、対応する演算子Ｆ（ｉ）（ｚ０，ｚｍ）と関連付けることができ、ここで、ｚ０は特定のＣＦＰ取得に対応する受信機の位置である。前述の貯留層の例では、１，０００個のＣＦＰ取得に基づいて、１，０００個のＤＴＳ取得を生成できる。
【００３３】
ブロック２１０において、ＤＴＳ取得に基づいて、双方向の残余移動時間ΔＴ（ｉ）を特定できる。例えば、ＤＴＳ取得から反射イベントを識別でき、時間ゼロからの反射イベントの時間偏差に基づいて双方向残余移動時間ΔＴ（ｉ）を特定できる。前記貯留層の例に対しては、図３Ａ（後述する）は、１から１０００までのインデックスを有する１，０００個の受信機に対応する１，０００個のＤＴＳ取得を示す。赤線３０２は反射イベントを示す。各ＤＴＳ取得に対して、赤線３０２と時間ゼロとの間の時間差に等しくなるように１つのΔＴ（ｉ）を特定でき、さらに、図３Ａから１，０００個の受信機に対して１，０００個のΔＴ（ｉ）を特定できる。実施によっては、すべての取得表面位置に対する残余移動時間の完全なセットを生成するために、ＤＴＳ取得から得られる双方向の残余移動時間を補間できる。例えば、１０１，０００箇所の取得表面位置に対する１０１，０００個のΔＴ（ｉ）を生成するために、図３から得られる１，０００個のΔＴ（ｉ）を補間できる。ドロネー三角形分割法のような補間方法を用いることができる。場合によっては、地震データを震源位置に基づいて分類でき、さらに、１００，０００箇所の震源位置に対応する１００，０００個のΔＴ（ｉ）を生成するために分類されたデータ上で、ステップ２０６からステップ２１０を実行できる。
【００３４】
ブロック２１２において、ＤＴＳ取得からの（すべての）ΔＴ（ｉ）がゼロであるかゼロに近い（例えば、ゼロから所定の閾値内）かどうかを特定する。言い換えれば、ブロック２０２は、ＤＴＳ取得が時間ゼロで平坦な反射イベントを含むかどうかを特定する。時間ゼロでの平坦なイベントは、最終の集束演算子が見つかり、反復を停止できることを示す。ここで、方法２００はブロック２１８へ進む。ＤＴＳ取得が時間ゼロで平坦なイベントを含まない場合、方法２００はブロック２１４へ進み、演算子を繰り返し更新する。
【００３５】
ブロック２１４において、各取得表面位置に対する式（３）及び式（４）（又は、式（１）及び（２））に基づいて、演算子Ｆ（ｉ＋１）（ｚ０，ｚｍ）を特定できる。例えば、１，０００個の受信機に対応する１，０００個のΔＴ（ｉ）は、１，０００個の受信機位置に対するＣＦＰ演算子を更新するために用いられ、１００，０００個の震源位置に対応する１００，０００個のΔＴ（ｉ）は、１００，０００個の震源位置に対するＣＦＰ演算子を更新するために用いられる。
【００３６】
ブロック２１６において、反復カウンタｉは１だけ増加し、方法２００はブロック２０６へ戻り、受信機によって記録された元の地震データに、更新された演算子を適用して、ＣＦＰ及びＤＴＳ取得を生成する。
【００３７】
ブロック２１８では、最終の集束演算子を見つけた後（ＤＴＳ取得が平坦な反射イベントを示すとき）、最終の演算子からのＤＴＳ取得を追加することにより、ＤＴＳスタックを生成できる。前述の貯留層の例では、１，０００個のＤＴＳ取得を追加して１つのＤＴＳスタックを生成する。
【００３８】
ブロック２２０において、目標反射体ｚｍの反射係数を、時間ゼロでのＤＴＳスタックの値であると特定できる。
【００３９】
方法２００を各反射体ｚｍに適用した後、地震画像を、地震データセットから生成でき、その地震画像は、貯留層内のすべての反射体における反射係数を含む。実施によっては、生成された地震画像を３Ｄ画像とすることもできる。
【００４０】
図１に戻り、ブロック１０８において、異なる地震データセットに対応する反射体の反射係数（又は地震画像）を比較して、期間内の反射係数の変化を特定できる。例えば、月ごとに取得される地震データにより地震画像を生成できる。貯留層内の炭化水素液の組成や位置の変化のような、貯留層内の変化をモニタするために、１８か月間における１８個の地震画像を比較できる。
【００４１】
図３Ａ乃至図５Ｂは、説明された手法の利用を示し、ＣＯ２が貯留層に注入されるプロジェクトにおいて取得された時間経過のフィールド（領域）地震データセットから、効率的で精確な３Ｄの目標指向の地震画像を生成するものである。このプロジェクトでは、密集したショット分布（ｘ方向とｙ方向との両方向で１０ｍ間隔）及び地面下７０ｍの深さに約１，０００個の埋設された受信機を用いて、毎月フル３Ｄ地震探査を実施する。受信機は、地表面近くの複雑さの影響を軽減し、調査間の再現性を高められるように、埋設される。集束演算子が導出されるＣＦＰ及びＤＴＳ取得を生成するために相互関係性を採用する。ＣＯ２注入の前に取得したベースライン（基準値）調査は、ＣＦＰ及びＤＴＳ取得を用いて集束演算子を生成するために用いられる。この手法を時間経過データセットへ適用することにより、ＣＯ２注入履歴と一致する、貯留層レベルにおける変化が明らかになる。
【００４２】
図３Ａは、いくつかの実施による、１つの目標反射体に対する初期集束演算子を用いたＤＴＳ取得のセット３００ａを示す。図３Ａにおいて、横軸はＤＴＳ取得インデックスを表し、縦軸は秒単位の時間を表す。前述のように、ＤＴＳ取得のセット３００ａは、ＮＭＯ速度から導出される、初期集束演算子Ｆ（０）（ｚ０，ｚｍ）（図４Ａに示され、後述する）を用いて得られた約１，０００個のＤＴＳ取得を含む。ＤＴＳ取得のセット３００ａは、赤線３０２で示される自動選択された反射イベントを含む。このプロジェクトの最大オフセットは３ｋｍであるため、初期集束演算子からのＤＴＳ取得は、時間ゼロにおいて、ほぼフラットなイベントを含むように見える一方、残余時間の変動は、反射イベント３０２から計算でき、集束演算子の更新に用いることができる。 
【００４３】
図３Ｂは、いくつかの実施による、１つの目標反射体に対する更新された集束演算子を用いたＤＴＳ取得のセット３００ｂを示す。図３Ｂにおいて、横軸はＤＴＳ取得インデックスを表し、縦軸は秒単位の時間を表す。ＤＴＳ取得のセット３００ｂは、１回の反復に対して初期集束演算子を更新した後、更新された演算子Ｆ（１）（ｚ０，ｚｍ）（図４Ｂに示され、後述する）を用いた約１，０００個のＤＴＳ取得を含む。ＤＴＳ取得のセット３００ｂは、時間ゼロでの平坦なイベントである反射イベント３０４を含む。
【００４４】
図４Ａは、いくつかの実施による、１つの目標反射体に対する初期集束演算子４００ａを示す。図４Ａは、目標反射体とすべての取得表面位置との間の初期集束演算子を示す。水平軸と垂直軸は、それぞれ取得表面のｘ方向とｙ方向とのグリッド（格子）インデックスを表す。緑色の点４０２は目標反射体の位置を示す。カラーバー４０８は、カラースペクトルの色と秒単位の移動時間の連続範囲との間のマッピング（写像）を描く。例えば、符号４０６で示す赤色の領域は、約１．４秒の移動時間を持つ集束演算子を示し、それに対し、符号４０４で示す紺青色の領域は、約０．６秒の移動時間を持つ集束演算子を示す。同様に、図４Ｂは、いくつかの実施による、１つの目標反射体に対する１回の反復後に更新された集束演算子４００ｂを示す。図４Ｃは、いくつかの実施による、１つの目標反射体に対する初期集束演算子４００ａと更新された集束演算子４００ｂとの間の差４００ｃを示す。カラーバー４１０は、カラースペクトルの色と秒単位の移動時間の連続範囲との間のマッピングを描く。例えば、符号４１４で示す赤色の領域は、約０の移動時間差を示し、それに対し、符号４１６で示す紺青色の領域は、約１５０ミリ秒の移動時間差を示す。青色の点４１２は、埋設された受信機の位置を示す。図４Ｃにおいて青色の点で示す埋設された受信機はまばらなので、すべての取得表面位置における更新された残余移動時間を生成するために、補間及び外挿の手段としてドロネー三角形分割法が用いられる。
【００４５】
図５Ａは、いくつかの実施による、目標地平線における反射係数の正規化された二乗平均平方根（ＮＲＭＳ）振幅５００ａを示す。図５Ａにおいて、水平軸及び垂直軸は、それぞれ、目標地平線のｘ方向及びｙ方向のグリッドインデックスを表す。図５Ａは、２番目の月に取得された調査２からの反射係数のＮＲＭＳ振幅を示す。黒色の四角５０２、５０４、５０６、及び５０８は、４箇所のＣＯ２注入位置を示す。カラーバー５１０は、カラースペクトルの色と、ＮＲＭＳ振幅値の連続範囲との間のマッピングを描く。例えば、符号５１４で示す赤色の領域は、１に近い値のＮＲＭＳ振幅を示し、それに対し、符号５１２で示す紺青色の領域は、ゼロに近い値のＮＲＭＳ振幅を示す。異なる調査間の表土速度は同じであると想定されるため、調査２（調査２及びそれに続く調査は、表面の部分的な統合により再現性が得られたため、ベースライン調査として用いられる）から導き出された同じ集束演算子は地震画像を生成するための異なる調査に適用される。
【００４６】
図５Ｂは、いくつかの実施による、調査２と調査１７との間のＮＲＭＳ振幅差５００ｂを示す。調査１７（不図示）は調査２の１５か月後に行われる。カラーバー５１６は、カラースペクトルの色とＮＲＭＳ振幅差の連続範囲との間のマッピングを描く。例えば、符号５１８で示す橙色の領域は、値が１０のＮＲＭＳ振幅差を示し、それに対して符号５２０で示す水色の領域は、値が約−７のＮＲＭＳ振幅差を示す。赤色の円５２２の異常は、ＣＯ２注入に関連すると考えられるが、貯留層シミュレータを用いて生成された予測と一致し、ＣＯ２プルームとの境界をより明確に画成する。
【００４７】
図６は、例示のコンピュータシステム６００のブロック図であり、ある実装に従い、本開示に記載されているように、記述されたアルゴリズム、方法、機能、プロセス、フロー、及び手順に関連付けられた計算機能を提供するために使用される。図示されたコンピュータ６０２は、サーバ、デスクトップコンピュータ、ラップトップ／ノートブックコンピュータ、無線データポート、スマートフォン、パーソナルデータアシスタント（ＰＤＡ）、タブレットコンピューティングデバイス、又はこれらの内部にある一又は複数のプロセッサといった任意のコンピューティングデバイスを包含することを意図しており、コンピューティングデバイスは、コンピューティングデバイスの物理インスタンス又は仮想インスタンス（又はこれらの両方）を含む。加えて、コンピュータ６０２は、入力デバイス（例えば、キーパッド、キーボード、又はタッチスクリーン）と、コンピュータ６０２の操作に関連付けられた情報を伝達する（例えば、グラフィカルユーザインタフェース（ＧＵＩ）上のデジタルデータ、ビジュアルデータ、オーディオ情報（情報の組み合わせ）を伝達する）出力デバイスと、を含むコンピュータを備えてもよい。
【００４８】
コンピュータ６０２は、本開示で説明される主題を実行するために、クライアント、ネットワークコンポーネント、サーバ、データベース、又は役割の組み合わせとして役割を果たし得る。図示のコンピュータ６０２は、ネットワーク６３０に通信可能に結合されている。いくつかの実施では、コンピュータ６０２の一又は複数のコンポーネントは、クラウドコンピューティングベースの環境、ローカル環境、グローバル環境、又は組み合わせの環境を含む環境内において動作するように構成できる。
【００４９】
高レベルでは、コンピュータ６０２は、記載されている主題に関連付けられたデータ及び情報を受信し、送信し、処理し、格納し、又は管理するように動作可能な電子式計算デバイスである。いくつかの実施によれば、コンピュータ６０２は、アプリケーションサーバ、電子メールサーバ、ウェブサーバ、キャッシングサーバ、ストリーミングデータサーバ、又はサーバの組み合わせを含むか通信可能に結合されていてもよい。
【００５０】
コンピュータ６０２は、ネットワーク６３０を介してクライアントアプリケーション（例えば、別のコンピュータ６０２上で動作する）からリクエストを受けることができ、適切なソフトウェアアプリケーションを用いて受信リクエストを処理することによって受信リクエストに応答することができる。加えて、リクエストは、コンピュータ６０２へ、内部ユーザ（例えば、コマンドコンソール）、外部ユーザ又はサードパーティ、又は自動化されたアプリケーションから、同様に、あらゆる他の適切なエンティティ、個人、システム、又はコンピュータから送ってもよい。
【００５１】
コンピュータ６０２の各コンポーネントは、システムバス６０３を用いて通信することができる。いくつかの実施では、コンピュータ６０２のうちのいずれか又はすべてのコンポーネント、ハードウェア又はソフトウェア（又はハードウェア及びソフトウェアの組み合わせ）は、アプリケーションプログラミングインタフェース（ＡＰＩ）６１２、又はサービスレイヤー６１３（又はＡＰＩ６１２及びサービスレイヤー６１３の組み合わせ）を用いてシステムバス６０３上で互いに又はインタフェース６０４（又は両者の組み合わせ）とインタフェース接続してもよい。ＡＰＩ６１２は、ルーチン、データ構造、及びオブジェクトクラスの仕様を含んでもよい。ＡＰＩ６１２は、コンピュータ言語に依存しない又は依存してもよく、完全なインタフェース、又は単一の機能を示し、或いは一組のＡＰＩさえも示す。サービス層６１３は、コンピュータ６０２又はコンピュータ６０２に通信可能に結合された他のコンポーネント（図示されているか否かにかかわらず）にソフトウェアサービスを提供する。コンピュータ６０２の機能は、このサービス層を用いてすべてのサービス消費者にアクセスしてもよい。サービス層６１３によって提供されるサービスといったソフトウェアサービスは、定義されたインタフェースを介して再利用可能な定義された機能を提供する。例えば、インタフェースは、Ｊａｖａ（登録商標）、Ｃ＋＋、又はコンピュータ言語の組み合わせで記述されたソフトウェアであって、これらの言語は、拡張マークアップ言語（ＸＭＬ）フォーマット、又はフォーマットの組み合わせで、データを提供する。コンピュータ６０２の統合コンポーネントとして示されているが、代替の実施は、コンピュータ６０２の他のコンポーネント又はコンピュータ６０２に通信可能に結合された他のコンポーネント（図示されているか否かにかかわらず）に関連して、スタンドアロンコンポーネントとしてＡＰＩ６０２又はサービス層６１３を例示してもよい。さらには、ＡＰＩ６１２又はサービス層６１３のいずれか又はすべての部分は、この開示の範囲から逸脱することなく、別のソフトウェアモジュール、企業アプリケーション、又はハードウェアモジュールの子モジュール又はサブモジュールとして実現してもよい。
【００５２】
コンピュータ６０２は、インタフェース６０４を含む。図６では単一のインタフェース６０４として示されているが、コンピュータ６０２の特定のニーズ、要求、又は特定の実施に従って２以上のインタフェース６０４を使用してもよい。インタフェース６０４は、分散環境においてネットワーク６３０に通信可能に結合された別のシステム（図示されているか否かにかかわらず）と通信するために、コンピュータ６０２によって使用される。一般に、インターフェ−ス６０４は、ソフトウェア又はハードウェア（又はソフトウェア及びハードウェアの組み合わせ）でエンコードされたロジックを含み、ネットワーク６３０と通信するように動作可能である。より具体的には、インタフェース６０４は、ネットワーク６３０又はインタフェースのハードウェアが図示のコンピュータ６０２の内側及び外側において物理信号を送受信するように動作可能であるように、通信に関連付けられた一又は複数の通信プロトコルをサポートするソフトウェアを備えていてもよい。
【００５３】
コンピュータ６０２は、プロセッサ６０５を含む。図６では単一のプロセッサ６０５として示されているが、コンピュータ６０２の特定のニーズ、要望、又は特定の実施に従って、２以上のプロセッサ６０５を使用してもよい。一般的には、プロセッサ６０５は命令を実行し、データを操作して、コンピュータ６０２の動作、並びに本開示において説明された任意のアルゴリズム、方法、機能、プロセス、フロー、及び手順を実行する。
【００５４】
コンピュータ６０２は、コンピュータ６０２又はネットワーク６３０に結合することができる別のコンポーネント（又は両者の組み合わせ）（図示されているか否かにかかわらず）のためのデータを保持できるデータベース６０６をも含む。例えば、データベース６０６は、本開示に整合するデータを格納するインメモリ又は従来のデータベースであり得る。いくつかの実施では、データベース６０６は、コンピュータ６０２及び記述された機能の特定のニーズ、要望、又は特定の実施に従って、二又はそれより多くの異なるデータベースタイプの組み合わせ（例えば、インメモリと従来型とのハイブリッドデータベース）であり得る。図６では単一のデータベース６０６として示されているが、（同じ又は複数のタイプの組み合わせの）二又はそれより多くのデータベースは、コンピュータ６０２及び記述された機能の特定のニーズ、要望、又は特定の実施に従って使用することができる。データベース６０６は、コンピュータ６０２の不可欠の構成要素として示されているが、代替の実施では、データベース６０６は、コンピュータ６０２の外部にあってもよい。例えば、データベース６０６は、連続的な地震貯留層の監視のための地震データセットを保持することができる。
【００５５】
また、コンピュータ６０２は、メモリ６０７を含み、このメモリ６０７は、コンピュータ６０２又はネットワーク６３０に結合することができる別のコンポーネント（又は両者の組み合わせ）（図示されているか否かにかかわらず）のためのデータを保持できる。例えば、メモリ６０７は、ランダムアクセスメモリ（ＲＡＭ）、リードオンリーメモリ（ＲＯＭ）、光学、磁気などの、本開示に整合するデータを格納できる。いくつかの実施では、メモリ６０７は、コンピュータ６０２及び記述された機能の特定のニーズ、要望、又は特定の実施に係る２以上の異なるタイプのメモリの組み合わせ（例えば、ＲＡＭ及び磁気ストレージの組み合わせ）であり得る。図６では単一のメモリ６０７として示されているが、（同じ又は複数のタイプの組み合わせの）二又はより多くのメモリ６０７は、コンピュータ６０２及び記述された機能の特定のニーズ、要望、又は特定の実施に従って使用することができる。メモリ６０７は、コンピュータ６０２の不可欠なコンポーネントとして示されているが、代替の実施では、メモリ６０７はコンピュータ６０２の外部にあってもよい。
【００５６】
アプリケーション６０８は、特に本開示で記述された機能に関して、コンピュータ６０２の特定のニーズ、要望、又は特定の実施に係る機能を、提供するアルゴリズムソフトウェアエンジンである。例えば、アプリケーション６０８は、一又は複数のコンポーネント、モジュール、又はアプリケーションとして機能できる。さらに、単一のアプリケーション６０８として示されているが、アプリケーション６０８は、コンピュータ６０２上の複数のアプリケーション６０８として実施してもよい。加えて、代替の実施では、コンピュータ６０２と一体として示されているが、アプリケーション６０８は、コンピュータ６０２の外部にあってもよい。
【００５７】
任意の数のコンピュータ６０２が、このコンピュータ６０２を含むコンピュータシステムに関連付けられ、又は該コンピュータシステムの外部にあり、各コンピュータ６０２はネットワーク６３０を介して通信する。さらに、用語「クライアント」、「ユーザ」、及び他の適切な用語は、本開示の範囲から逸脱することなく、必要に応じて、交換可能に使用されてもよい。さらに、本開示では、多くのユーザが１台のコンピュータ６０２を使用できることや、１人のユーザが複数のコンピュータ６０２を使用できることを想定している。
【００５８】
主題の記載された実施は、単独で又は組み合わせて、一又は複数の特徴を含むことができる。
【００５９】
例えば、第１の実施では、方法は、一定の期間内に、貯留槽に関連付けられた複数の地震データセットを受信するステップを含み、前記複数の地震データセットは前記期間内の異なる時間に受信した地震データに対応する。前記貯留槽は複数の反射体を含み、各反射体は反射係数を有する。前記複数の地震データセットのそれぞれに対して、前記複数の反射体の反射係数は、前記複数の反射体と複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップにより特定される。異なる地震データセットに対応する前記反射係数は、前記期間内の前記反射係数の変化を特定するために比較される。前記反射係数の前記変化が表示される。
【００６０】
前述及び他の説明された実施は、各々、オプションとして、以下の特徴のうちの一又は複数を含むことができる。
【００６１】
１つ目の特徴は、以下の特徴のいずれかと組み合わせ可能であって、各ＣＦＰ演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す。
【００６２】
２つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の反射体と前記複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新されたＣＦＰ演算子を
【数３】
2020511656.tif 000004
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間移動（ＤＴＳ）取得に基づいて特定された残余伝播時間である。
【００６３】
３つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期ＣＦＰ演算子を特定するステップをさらに備える。
【００６４】
４つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える。
【００６５】
５つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える。
【００６６】
６つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記貯留層が炭化水素液を備える。
【００６７】
第２の実施では、システムは、コンピュータメモリと、前記コンピュータメモリと相互運用可能に接続された１つ以上のハードウェアプロセッサとを備える。前記１つ以上のハードウェアプロセッサは、一定の期間内に、貯留槽に関連付けられた複数の地震データセットを受信するステップを含む操作を実行するように構成され、前記複数の地震データセットは、前記期間内の異なる時間に受信した地震データに対応する。前記貯留槽は複数の反射体を含み、各反射体は反射係数を有する。前記複数の地震データセットのそれぞれに対して、前記複数の反射体の反射係数は、前記複数の反射体と複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップにより特定される。異なる地震データセットに対応する前記反射係数は、前記期間内の前記反射係数の変化を特定するために比較される。前記反射係数の前記変化が表示される。
【００６８】
前述及び他の説明された実施は、各々、オプションとして、以下の特徴のうちの一又は複数を含むことができる。
【００６９】
１つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、各ＣＦＰ演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す。
【００７０】
２つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の反射体と前記複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新されたＣＦＰ演算子を
【数４】
2020511656.tif 000005
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間移動（ＤＴＳ）取得に基づいて特定された残余伝播時間である。
【００７１】
３つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記操作は、前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期ＣＦＰ演算子を特定するステップをさらに備える。
【００７２】
４つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記操作は、前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える。
【００７３】
５つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える。
【００７４】
６つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記貯留層が炭化水素液を備える。
【００７５】
第３の実施は、一定の期間内に、貯留槽に関連付けられた複数の地震データセットを受信するステップを含む操作を実行するコンピュータシステムによって実行可能な１つ以上の指令が格納された非一時的なコンピュータ読み取り可能な記録媒体であって、前記複数の地震データセットは、前記期間内の異なる時間に受信した地震データに対応する。前記貯留槽は複数の反射体を含み、各反射体は反射係数を有する。前記複数の地震データセットのそれぞれに対して、前記複数の反射体の反射係数は、前記複数の反射体と複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップにより特定される。異なる地震データセットに対応する前記反射係数は、前記期間内の前記反射係数の変化を特定するために比較される。前記反射係数の前記変化が表示される。
【００７６】
前述及び他の説明された実施は、各々、オプションとして、以下の特徴のうちの一又は複数を含むことができる。
【００７７】
１つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、各ＣＦＰ演算子は、特定の反射体と特定の取得表面位置との間の伝播時間を示す。
【００７８】
２つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の反射体と前記複数の取得表面位置とに関連付けられたＣＦＰ演算子を繰り返し更新するステップは、特定の反射体ｚｍ及び特定の取得表面位置ｚ０に対して、（ｉ＋１）番目に更新されたＣＦＰ演算子を
【数５】
2020511656.tif 000006
を用いて計算するステップを備え、
ΔＴ（ｉ）は、ｉ番目の反復における差分時間移動（ＤＴＳ）取得に基づいて特定された残余伝播時間である。
【００７９】
３つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記操作は、前記特定の反射体ｚｍ及び前記特定の取得表面位置ｚ０に対して、前記特定の反射体ｚｍと前記特定の取得表面位置ｚ０との間の伝播時間の初期推定に基づいて初期ＣＦＰ演算子を特定するステップをさらに備える。
【００８０】
４つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記操作は、前記複数の地震データセットのそれぞれに対して、前記複数の反射体の前記特定された反射係数に基づいて、前記貯留層の三次元地震画像を生成するステップをさらに備える。
【００８１】
５つ目の特徴は、既述又は以下の特徴のいずれかと組み合わせ可能であって、前記複数の取得表面位置は、地震源の位置及び地震受信機の位置を備える。
【００８２】
本明細書で記述された主題及び機能的な操作の実施は、デジタル電子回路、有形に具体化されたコンピュータのソフトウェア又はファームウェア、コンピュータハードウェアに、又はそれらの一又は複数の組み合わせで実施でき、これらデジタル電子回路、ソフトウェア、ファームウェア、及びコンピュータハードウェアは、本明細書に開示される構造及びそれらの構造的な同等物を含む。本明細書に記載された主題における実施は、一又は複数のコンピュータプログラム、つまり、コンピュータプログラム命令の一又は複数のモジュールとして実施されることができ、コンピュータプログラム命令は、データ処理装置による実行のために、有形の非一時的なコンピュータ読取可能なコンピュータ貯蔵媒体上にエンコードされ、又はコンピュータ若しくはコンピュータにより実施されるシステムの動作を制御する。代替的に又は追加的に、プログラム命令は、人工的に生成された伝播信号、例えば、マシン生成の電気的、光学的、又は電磁気的な信号にエンコードされ、この信号は、データ処理装置による実行のために適切な受振器装置への送信用の情報をエンコードするために生成される。コンピュータ記憶媒体は、機械により読取可能な記憶装置、機械により読取可能な記憶基板、ランダム若しくはシリアルアクセスのメモリデバイス、又はコンピュータ記憶媒体の組み合わせであり得る。
【００８３】
「リアル−タイム」、「リアルタイム」、「リアル（ファースト）タイム（ＲＦＴ）」、「ほぼリアルタイム（ＮＲＴ）」、「準リアルタイム」、又は（当業者が理解する）類似の用語は、個々が、実質的に同時に動作及び応答が生じることを知覚するように、該動作及び応答が時間的に近接していることを意味する。例えば、データにアクセスするための個々の動作に続くデータの表示（又は表示の開始）に対する応答の時間差は、１ｍｓ（ミリ秒）未満、１秒未満、又は５秒未満であってもよい。要求されたデータが、即時に表示（又は表示開始）される必要はないが、記述されたコンピューティングシステムの処理制限と、データを例えば収集し、正確に評価し、分析し、処理し、格納し、又は送信するために必要な時間とを考慮して、意図的な遅延なしに表示（又は表示開始）される。
【００８４】
用語「データ処理装置」、「コンピュータ」、又は「電子コンピュータデバイス」（又は当業者によって理解される同等のもの）は、データ処理ハードウェアを指し、あらゆる種類の装置、デバイス、及びマシンを包含し、これらの種類は、データを処理するためのものであり、例示として、プログラマブルプロセッサ、コンピュータ、又は複数のプロセッサ若しくはコンピュータを含む。装置は、例えば、中央処理装置（ＣＰＵ）、ＦＰＧＡ（フィールドプログラマブルゲートアレイ）、又はＡＳＩＣ（特定用途向け集積回路）を包含する専用ロジック回路であるか、又はこれらを更に含むことができる。いくつかの実施では、データ処理装置、又は専用ロジック回路（又はデータ処理装置と専用ロジック回路との組み合わせ）は、ハードウェア系又はソフトウェア系（又はソフトウェア系及びハードウェア系の両方の基づく組み合わせ）であってもよい。装置は、オプションとして、コンピュータプログラムの実行環境を生成するコード、例えば、プロセッサファームウェア、プロトコルスタック、データベース管理システム、オペレーティングシステム、又は実行環境の組み合わせを構成するコードを含むことができる。本開示は、従来のオペレーティングシステムを持つ又は持たないデータ処理装置の使用を想定しており、オペレーティングシステムは、例えば、ＬＩＮＵＸ、ＵＮＩＸ（登録商標）、ＷＩＮＤＯＷＳ（登録商標）、ＭＡＣＯＳ、ＡＮＤＲＯＩＤ（登録商標）、ＩＯＳ、又は複数のオペレーティングシステムの組み合わせである。
【００８５】
コンピュータプログラムは、プログラム、ソフトウェア、ソフトウェアアプリケーション、ユニット、モジュール、ソフトウェアモジュール、スクリプト、又はコードとして参照され又は記述されてもよいものであるが、コンパイル言語若しくはインタープリタ言語を含むプログラミング言語、又は宣言型若しくは手続き型の言語の形式で記載できる。コンピュータプログラムは、スタンドアロンプログラムとして、又はモジュール、コンポーネント、サブルーチン、又はコンピューティング環境で使用するために適した他のユニットとして含む、任意の形式で展開することができる。コンピュータプログラムは、ファイルシステム内のファイルに対応してもよいが、必ずしも対応する必要はない。プログラムは、他のプログラム又はデータを保持するファイルの一部、例えばマークアップ言語ドキュメントに格納された一又は複数のスクリプトに保存でき、該他のプログラムは、問題のプログラム専用の単一ファイル内に、或いは複数の連携したファイル、例えば一又は複数のモジュール、サブプログラム、又はコードの一部に格納できる。コンピュータプログラムは、一又は複数のコンピュータ上に展開でき、該コンピュータは、一のサイトに位置し、又は複数のサイトに分散されて配置され、これらは通信ネットワークによって相互接続される。様々な図面に例示されたプログラムの部分は、様々なオブジェクト、方法、又は他のプロセスを通して特徴及び機能を実施する個々のモジュールとして示されているが、プログラムは、それらの代わりに、いくつかのサブユニット、サブモジュール、サードパーティのサービス、コンポーネント、ライブラリなどを、必要に応じて含んでもよい。逆に、必要に応じて、様々なコンポーネントの機能及び特徴は、単一のコンポーネントに組み合わせることができる。計算上の特定を為すために使用されたしきい値は、静的に、動的に、又は静的と動的との両方で特定できる。
【００８６】
本明細書に記載された方法、プロセス、又は論理フローは、一又は複数のプログラム可能なコンピュータによって実行でき、該コンピュータは、入力データを操作して出力データを生成することによって機能を実行する一又は複数のコンピュータプログラムを実行する。方法、プロセス、又は論理フローは、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実行でき、また装置も、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実施できる。
【００８７】
コンピュータプログラムの実行に適したコンピュータは、汎用又は特定用途のマイクロプロセッサに基づくことができる。一般的には、ＣＰＵは、リードオンリーメモリ（ＲＯＭ）、又はランダムアクセスメモリ（ＲＡＭ）、又はその両方から命令及びデータを受け取るであろう。コンピュータの重要な要素は、命令を行い又は実行するためのＣＰＵ、及び命令及びデータを保存するための一又は複数のメモリデバイスである。一般に、コンピュータは、データを保存するための一又は複数の大容量記憶装置、例えば磁気、光磁気ディスク、又は光ディスクを含み、或いはこれらの大容量記憶装置に動作可能に結合され、データを受信し、転送し、又はこの両方を行う。しかし、コンピュータは、そのようなデバイスを持つ必要はない。さらには、コンピュータは、ほんの数例を挙げると、別のデバイス、例えば携帯電話、携帯情報端末（ＰＤＡ）、モバイルオーディオ若しくはビデオプレーヤー、ゲームコンソール、全地球測位システム（ＧＰＳ）受信機、又はポータブルストレージデバイス、例えばユニバーサルシリアルバス（ＵＳＢ）フラッシュドライブ、に組み込まれることができる。
【００８８】
コンピュータプログラム命令及びデータを格納するのに適したコンピュータ読取可能媒体（適宜、一時的又は非一時的）は、あらゆる形態の不揮発性のメモリ、媒体、及びメモリデバイスを含み、これらは、手段として、例えば、半導体メモリデバイス（例えば、消去可能プログラマブルリードオンリーメモリ（ＥＰＲＯＭ）、電気的消去可能プログラマブル読み取り専用メモリ（ＥＥＰＲＯＭ）、及びフラッシュメモリデバイス）、磁気デバイス（例えば、内部ハードディスク又はリムーバブルディスク）、光磁気ディスク、光学メモリデバイス（例えば、ＣＤ−ＲＯＭ、ＤＶＤ＋／−Ｒ、ＤＶＤ−ＲＡＭ、ＤＶＤ−ＲＯＭディスク）を含む。メモリは、様々なオブジェクト又はデータを格納してもよく、オブジェクト又はデータは、キャッシュ、クラス、フレームワーク、アプリケーション、バックアップデータ、ジョブ、Ｗｅｂページ、Ｗｅｂページテンプレート、データベーステーブル、動的情報を格納するリポジトリ、及び任意の他の適切な情報を包含しており、他の適切な情報は、任意のパラメータ、任意の変数、任意のアルゴリズム、任意の命令、任意のルール、任意の制約、又は任意の参照を含む。さらに、メモリは、ログ、ポリシー、セキュリティ又はアクセスのデータ、レポートファイルといった任意の他の適切なデータを含んでいてもよい。プロセッサ及びメモリは、特定用途ロジック回路によって補完され又は組み込むことができる。
【００８９】
ユーザとのインタラクションを提供するために、本明細書で記載された主題の実施は、ユーザに情報を表示するためのディスプレイデバイス、例えばＣＲＴ（陰極線管）、ＬＣＤ（液晶ディスプレイ）、ＬＥＤ（発光ダイオード）又はプラズマモニターと、ユーザがコンピュータに入力を提供するキーボード及びポインティングデバイス、例えばマウス、トラックボール、トラックパッドとを有するコンピュータ上で実施できる。入力は、タブレットコンピュータの感圧性表面、又は静電容量式若しくは電気式のセンシングを使用したマルチタッチスクリーンといった、タッチスクリーンを用いてコンピュータに提供してもよい。他の種類のデバイスもまた、ユーザとのインタラクション用に提供するために使用できる。例えば、ユーザに提供されるフィードバックは、視覚フィードバック、聴覚フィードバック、又は触覚フィードバックといった、あらゆる形態の感覚のフィードバックであり得る。ユーザからの入力は、音響入力、音声入力、触覚入力を包含する、あらゆる形式で受信できる。さらに、コンピュータは、ユーザによって使用されるデバイスからのドキュメントを受信し及び該ドキュメントをデバイスに送信することによって、例えば、Ｗｅｂブラウザから受けたリクエストに応答してユーザのクライアントデバイス上においてＷｅｂブラウザにＷｅｂページを送ることにより、ユーザとインタラクションできる。
【００９０】
用語「グラフィカルユーザインタフェース」又は「ＧＵＩ」は、単数又は複数で使用して、一又は複数のグラフィカルユーザインタフェース及び特定のグラフィカルユーザインタフェースのディスプレイの各々を説明できる。したがって、ＧＵＩは、Ｗｅｂブラウザ、タッチスクリーン、又はコマンドラインインタフェース（ＣＬＩ）を含むが、これらに限定されることなく、任意のグラフィカルユーザインタフェースを提示することができ、コマンドラインインタフェースは、情報を処理すると共に情報の結果をユーザに効率的に提示する。一般的に、ＧＵＩは、複数のユーザインタフェース（ＵＩ）要素、例えばインタラクティブフィールド、プルダウンリスト、ボタンを含むことができ、これらのうちのいくつか又はすべては、Ｗｅｂブラウザに関連付けられる。これら及び他のＵＩ要素は、Ｗｅｂブラウザの機能に関連付けられまたＷｅｂブラウザの機能を表すことができる。
【００９１】
本明細書で記述された主題の実施は、バックエンドコンポーネント（例えばデータサーバとして）を含むコンピューティングシステム、ミドルウェアコンポーネント（例えばアプリケーションサーバとして）を含むコンピューティングシステム、フロントエンドコンポーネント（例えばユーザが本明細書においで記述される主題の実施とインタラクションできるグラフィカルユーザインタフェース又はＷｅｂブラウザを有するクライアントコンピュータとして）を含むコンピューティングシステム、又は、このようなバックエンド、ミドルウェア、及びフロントエンドコンポーネントの一又は複数の任意の組み合わせにおいて実施できる。システムのコンポーネントは、有線又は無線のデジタルデータ通信（又はデータ通信の組み合わせ）、例えば通信ネットワークの任意の形式又は媒体によって相互接続されることができる。通信ネットワークの例示は、ローカルエリアネットワーク（ＬＡＮ）、無線アクセスネットワーク（ＲＡＮ）、メトロポリタンエリアネットワーク（ＭＡＮ）、広域ネットワーク（ＷＡＮ）、マイクロ波アクセスの世界的な相互運用性（ＷＩＭＡＸ）、例えば８０２．１１ａ／ｂ／ｇ／ｎ又は８０２．２０（又は８０２．１１ｘと８０２．２０の組み合わせ、又は本開示に整合する他のプロトコル）を用いる無線ローカルエリア（ＷＬＡＮ）、インターネットのすべて若しくは一部、又は通信ネットワークの組み合わせを含む。ネットワークは、ネットワークアドレス間のデータ、例えば、インターネットプロトコル（ＩＰ）パケット、フレームリレーフレーム、非同期転送モード（ＡＴＭ）セル、音声、又はビデオと通信できる。
【００９２】
コンピューティングシステムは、クライアント及びサーバを含むことができる。クライアント及びサーバは、一般的には、互いに離れており、また典型的には、通信ネットワークを介してインタラクションする。クライアント及びサーバの関係は、それぞれのコンピュータ上で動作すると共に互にクライアント・サーバ関係にあるコンピュータプログラムのおかげで生じる。
【００９３】
本明細書は多くの特定の実施の詳細を含む一方で、これらは、請求される得る範囲の制限として解釈されるべきではなく、むしろ、特定の概念の特定の実施固有の特徴の説明として解釈されるべきである。個別の実施の観点で本明細書において説明された機能は、組み合わせて又は単一の実施で、実現できる。逆に、単一の実施の観点で既述された様々な特徴は、複数の実施で、個別に、又は任意の適切なサブコンビネーションで実施されることもできる。さらには、既述の特徴は特定の組み合わせで動作するものとして説示され、最初はそのようなものとして請求されているが、一又は複数の特徴は、請求された組み合わせから場合によっては削除され、またサブコンビネーション又はサブコンビネーションの変形であり得る。
【００９４】
主題の特定の実施が説明された。記載された実施の置換、変更、及び他の実施は、当業者には明らかである以下の請求の範囲内である。操作は特許請求の範囲に又は特定の順序で図面に描かれている一方で、これは、望ましい結果を達成するために、示された特定の順序で又は順番でそのような操作が実行されること又はすべての図示された操作が実行される（いくつかの操作は随意的と見なされる）ことを要求するものとして理解されるべきではない。ある状況では、マルチタスク又は並列処理（又はマルチタスク及び並列処理の組み合わせ）を行うことが有利であり、適切と思われる場合に実行される。
【００９５】
さらに、既述の実施における様々なシステムモジュール及びコンポーネントの分離又は統合が、すべての実施においてそのような分離又は統合を必要とするものとして理解されるべきではなく、また、記載されたプログラムコンポーネント及びシステムは、一般的には、単一のソフトウェア製品に統合でき、又は複数のソフトウェア製品にパッケージ化できることが理解されるべきである。
【００９６】
これに従って、既述の例示的な実施は、本開示を定義し又は制約しない。他の変更、置換、及び変更も、本開示の範囲及び精神から逸脱することなく可能である。
【００９７】
さらには、請求された任意の実施は、少なくとも、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータシステムに適用可能であると考えられ、非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータにより読取可能な命令を格納し、またコンピュータシステムは、ハードウェアプロセッサに相互に動作可能に結合されたコンピュータメモリを含み、ハードウェアプロセッサは、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法又は命令を実行するように構成される。
【符号の説明】
【００９８】
１００  方法
２００  方法
３００ａ  初期集束演算子を用いたＤＴＳ取得のセット
３００ｂ  更新された集束演算子を用いたＤＴＳ取得のセット
４００ａ  初期集束演算子
４００ｂ  更新された集束演算子
４００ｃ  初期集束演算子と更新された集束演算子との間の差
５００ａ  反射係数の正規化されたＮＲＭＳ振幅
５００ｂ  調査２と調査１７との間のＮＲＭＳ振幅差
６００  コンピュータシステム
【図１】
2020511656.tif 000010
【図２】
2020511656.tif 000011
【図３Ａ】
2020511656.tif 000012
【図３Ｂ】
2020511656.tif 000013
【図４Ａ】
2020511656.tif 000014
【図４Ｂ】
2020511656.tif 000015
【図４Ｃ】
2020511656.tif 000016
【図５Ａ】
2020511656.tif 000017
【図５Ｂ】
2020511656.tif 000018
【図６】
2020511656.tif 000019
【国際調査報告】
2020511656.tif 000020
2020511656.tif 000021