【手続補正書】
【提出日】令和3年1月15日(2021.1.15)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
コンピュータにより実施される方法であって、
受信したショットギャザーを共通受振器ギャザーに分類するステップと；
合成データから、地震波動場の直接到着波のターゲット三次元（３Ｄ）振幅スペクトルを計算するステップと；
受振器毎にフィールドデータの地震波動場の直接到着波の３Ｄ振幅スペクトルを計算するステップと；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを計算するステップと；
フィルタリングされた下向き地震波動場を生成するために、前記計算された整合フィルタを、前記共通受振器ギャザーから分離された下向き地震波動場に適用するステップと；  平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行するステップと；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向き地震波動場との間の相互相関のグラフィカルユーザインタフェースを表示するステップであって、前記相互相関は地表付近の構造の低速度撮像による地震画像に関連する、ステップと；を備える、
コンピュータにより実施される方法。
【請求項２】
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択するステップと；
表面震源から埋設された受振器への前記一定のＰ波速度モデルを更新するステップと；
前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換するために、３Ｄフーリエ変換を実行するステップと；を更に備える、
請求項１に記載の方法。
【請求項３】
前記フィールドデータの前記３Ｄ振幅スペクトルと前記合成データの前記ターゲット３Ｄ振幅スペクトルとの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定するステップと；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを再計算するステップと；を更に備える、
請求項１に記載の方法。
【請求項４】
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離するステップを更に備える、
請求項１に記載の方法。
【請求項５】
各下向き地震波動場の時間ゲート処理を行うステップと；
各時間ゲート処理された下向き地震波動場に空間依存の開口フィルタリングを実行するステップと；を更に備える、
請求項４に記載の方法。
【請求項６】
前記上向き地震波動場に対してグラウンドロールの除去を実行するステップを更に備える、
請求項１に記載の方法。
【請求項７】
仮想震源ショットギャザーのスタックを行うステップと；
前記仮想震源ショットギャザーの出力を開始するステップと；を更に備える、
請求項１に記載の方法。
【請求項８】
コンピュータシステムによって実行可能な一又は複数の命令を格納する非一時的なコンピュータ読取可能媒体であって、
受信したショットギャザーを共通受振器ギャザーに分類する処理と；
合成データから、地震波動場の直接到着波のターゲット３次元（３Ｄ）振幅スペクトルを計算する処理と；
受振器毎に、フィールドデータにおける地震波動場の直接到着波の３Ｄ振幅スペクトルを計算する処理と；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを計算する処理と；
フィルタリングされた下向き地震波動場を生成するために、前記計算された整合フィルタを、前記共通受振器ギャザーから分離された下向き地震波動場に適用する処理と；
平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行する処理と；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向きの地震波動場との間の相互相関のグラフィカルユーザインタフェースを表示する処理であって、前記相互相関は地表付近の構造の低速度撮像による地震画像に関連する、処理と；をコンピュータに実行させるための前記命令を格納する、
非一時的なコンピュータ読取可能媒体。
【請求項９】
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択する処理と；
表面震源から埋設された受振器へ前記一定のＰ波速度モデルを更新する処理と；
３Ｄフーリエ変換を実行して、前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換する処理と；をコンピュータに実行させるための前記命令を更に格納する、  請求項８に記載の非一時的なコンピュータ読取可能媒体。
【請求項１０】
前記フィールドデータの前記３Ｄ振幅スペクトル及び前記合成データの前記ターゲット３Ｄ振幅スペクトルの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定する処理と；
前記フィールドデータの前記３Ｄ振幅スペクトルから前記ターゲット応答への整合フィルタを再計算する処理と；をコンピュータに実行させるための前記命令を更に格納する、
請求項８に記載の非一時的なコンピュータ読取可能媒体。
【請求項１１】
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離する処理をコンピュータに実行させるための前記命令を更に格納する、
請求項８に記載の非一時的なコンピュータ読取可能媒体。
【請求項１２】
各下向き地震波動場の時間ゲート処理を行う処理と；
各時間ゲート処理された下向き地震波動場に対して空間依存の開口フィルタリングを実行する処理と；をコンピュータに実行させるための前記命令を更に格納する、
請求項１１に記載の非一時的なコンピュータ読取可能媒体。
【請求項１３】
前記上向き地震波動場に対してグラウンドロールの除去を実行する処理をコンピュータに実行させるための前記命令を更に格納する、
請求項８に記載の非一時的なコンピュータ読取可能媒体。
【請求項１４】
仮想震源ショットギャザーをスタックする処理と；
前記仮想震源ショットギャザーの出力を開始する処理と；をコンピュータに実行させるための前記命令を更に格納する、
請求項８に記載の非一時的なコンピュータ読取可能媒体。
【請求項１５】
コンピュータにより実施されるシステムであって、
コンピュータメモリと；
前記コンピュータメモリと相互に操作可能に接続されたハードウェアプロセッサと；を備え、
前記ハードウェアプロセッサは、
受信したショットギャザーを共通受振器ギャザーに分類し；
合成データから地震波動場の直接到着波のターゲット３次元（３Ｄ）振幅スペクトルを計算し；
受振器毎にフィールドデータにおける地震波動場の直接到着波の３Ｄ振幅スペクトルを計算し；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを計算し；
フィルタリングされた下向き地震波動場を生成するために、前記共通受振器ギャザーから分離された下向き地震波動場に前記計算された整合フィルタを適用し；
平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行し；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向きの地震波動場との間の相互相関であって、地表付近の構造の低速度撮像による地震画像に関連する前記相互相関のグラフィカルユーザインタフェースを表示する；ように構成された、
コンピュータにより実施されるシステム。
【請求項１６】
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択し；
前記一定のＰ波速度モデルを表面震源から埋設された受振器へ更新し；
３Ｄフーリエ変換を実行して、前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換する；ように更に構成された、
請求項１５に記載のシステム。
【請求項１７】
前記フィールドデータの前記３Ｄ振幅スペクトルと合成データの前記ターゲット３Ｄ振幅スペクトルとの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定し；
前記フィールドデータの前記３Ｄ振幅スペクトルから前記ターゲット応答への整合フィルタを再計算する；ように更に構成された、
請求項１５に記載のシステム。
【請求項１８】
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離する；ように更に構成された、
請求項１５に記載のシステム。
【請求項１９】
各下向き地震波動場の時間ゲート処理を行い；
各時間ゲート処理された下向き地震波動場に対して空間依存の開口フィルタリングを実行する；ように更に構成された、
請求項１８に記載のシステム。
【請求項２０】
前記上向き地震波動場に対してグラウンドロールの除去を実行し；
仮想震源ショットギャザーをスタックし；
前記仮想震源ショットギャザーの出力を開始する；ように更に構成された、
請求項１５に記載のシステム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０１４９
【補正方法】変更
【補正の内容】
【０１４９】
さらには、請求された任意の実施は、少なくとも、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータシステムに適用可能であると考えられ、非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータにより読取可能な命令を格納し、またコンピュータシステムは、ハードウェアプロセッサに相互に動作可能に結合されたコンピュータメモリを含み、ハードウェアプロセッサは、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法又は命令を実行するように構成される。
以下、本発明の実施の態様の例を列挙する。
［第１の局面］
コンピュータにより実施される方法であって：
受信したショットギャザーを共通受振器ギャザーに分類するステップと；
合成データから、地震波動場の直接到着波のターゲット三次元（３Ｄ）振幅スペクトルを計算するステップと；
受振器毎にフィールドデータの地震波動場の直接到着波の３Ｄ振幅スペクトルを計算するステップと；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを計算するステップと；
フィルタリングされた下向き地震波動場を生成するために、前記計算された整合フィルタを、前記共通受振器ギャザーから分離された下向き地震波動場に適用するステップと；
平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行するステップと；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向き地震波動場との間の相互相関を計算するステップと；を備える、
コンピュータにより実施される方法。
［第２の局面］
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択するステップと；
表面震源から埋設された受振器への前記一定のＰ波速度モデルを更新するステップと；
前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換するために、３Ｄフーリエ変換を実行するステップと；を更に備える、
第１の局面に記載の方法。
［第３の局面］
前記フィールドデータの前記３Ｄ振幅スペクトルと前記合成データの前記ターゲット３Ｄ振幅スペクトルとの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定するステップと；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを再計算するステップと；を更に備える、
第１の局面に記載の方法。
［第４の局面］
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離するステップを更に備える、
第１の局面に記載の方法。
［第５の局面］
各下向き地震波動場の時間ゲート処理を行うステップと；
各時間ゲート処理された下向き地震波動場に空間依存の開口フィルタリングを実行するステップと；を更に備える、
第４の局面に記載の方法。
［第６の局面］
前記上向き地震波動場に対してグラウンドロールの除去を実行するステップを更に備える、
第１の局面に記載の方法。
［第７の局面］
仮想震源ショットギャザーのスタックを行うステップと；
前記仮想震源ショットギャザーの出力を開始するステップと；を更に備える、
第１の局面に記載の方法。
［第８の局面］
コンピュータシステムによって実行可能な一又は複数の命令を格納する非一時的なコンピュータ読取可能媒体であって：
受信したショットギャザーを共通受振器ギャザーに分類する処理と；
合成データから、地震波動場の直接到着波のターゲット３次元（３Ｄ）振幅スペクトルを計算する処理と；
受振器毎に、フィールドデータにおける地震波動場の直接到着波の３Ｄ振幅スペクトルを計算する処理と；
前記フィールドデータの前記３Ｄ振幅スペクトルから前記ターゲット応答への整合フィルタを計算する処理と；
フィルタリングされた下向き地震波動場を生成するために、前記計算された整合フィルタを、前記共通受振器ギャザーから分離された下向き地震波動場に適用する処理と；
平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行する処理と；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向きの地震波動場との間の相互相関を計算する処理と；をコンピュータに実行させるための前記命令を格納する、
非一時的なコンピュータ読取可能媒体。
［第９の局面］
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択する処理と；
表面震源から埋設された受振器へ前記一定のＰ波速度モデルを更新する処理と；
３Ｄフーリエ変換を実行して、前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換する処理と；をコンピュータに実行させるための前記命令を更に格納する、
第８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１０の局面］
前記フィールドデータの前記３Ｄ振幅スペクトル及び前記合成データの前記ターゲット３Ｄ振幅スペクトルの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定する処理と；
前記フィールドデータの前記３Ｄ振幅スペクトルから前記ターゲット応答への整合フィルタを再計算する処理と；をコンピュータに実行させるための前記命令を更に格納する、
第８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１１の局面］
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離する処理をコンピュータに実行させるための前記命令を更に格納する、
第８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１２の局面］
各下向き地震波動場の時間ゲート処理を行う処理と；
各時間ゲート処理された下向き地震波動場に対して空間依存の開口フィルタリングを実行する処理と；をコンピュータに実行させるための前記命令を更に格納する、
第１１の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１３の局面］
前記上向き地震波動場に対してグラウンドロールの除去を実行する処理をコンピュータに実行させるための前記命令を更に格納する、
第８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１４の局面］
仮想震源ショットギャザーをスタックする処理と；
前記仮想震源ショットギャザーの出力を開始する処理と；をコンピュータに実行させるための前記命令を更に格納する、
第８の局面に記載の非一時的なコンピュータ読取可能媒体。
［第１５の局面］
コンピュータにより実施されるシステムであって：
コンピュータメモリと；
前記コンピュータメモリと相互に操作可能に接続されたハードウェアプロセッサと；を備え、
前記ハードウェアプロセッサは：
受信したショットギャザーを共通受振器ギャザーに分類し；
合成データから地震波動場の直接到着波のターゲット３次元（３Ｄ）振幅スペクトルを計算し；
受振器毎にフィールドデータにおける地震波動場の直接到着波の３Ｄ振幅スペクトルを計算し；
前記フィールドデータの前記３Ｄ振幅スペクトルからターゲット応答への整合フィルタを計算し；
フィルタリングされた下向き地震波動場を生成するために、前記共通受振器ギャザーから分離された下向き地震波動場に前記計算された整合フィルタを適用し；
平滑化された下向き地震波動場を生成するために、前記フィルタリングされた下向き地震波動場の時間依存平滑化を実行し；
前記共通受振器ギャザーから分離された上向き地震波動場と前記平滑化された下向きの地震波動場との間の相互相関を計算する；ように構成された、
コンピュータにより実施されるシステム。
［第１６の局面］
初期の一定のＰ波速度モデルとして、オフセットに近い初期の地震波動場の直接到着波を選択し；
前記一定のＰ波速度モデルを表面震源から埋設された受振器へ更新し；
３Ｄフーリエ変換を実行して、前記更新された一定のＰ波速度モデルを周波数−波数ドメインに変換する；ように更に構成された、
第１５の局面に記載のシステム。
［第１７の局面］
前記フィールドデータの前記３Ｄ振幅スペクトルと合成データの前記ターゲット３Ｄ振幅スペクトルとの計算差と、許容値との比較が前記許容値より小さい又は等しいことを特定し；
前記フィールドデータの前記３Ｄ振幅スペクトルから前記ターゲット応答への整合フィルタを再計算する；ように更に構成された、
第１５の局面に記載のシステム。
［第１８の局面］
前記共通受振器ギャザーを下向き地震波動場及び上向き地震波動場に分離する；ように更に構成された、
第１５の局面に記載のシステム。
［第１９の局面］
各下向き地震波動場の時間ゲート処理を行い；
各時間ゲート処理された下向き地震波動場に対して空間依存の開口フィルタリングを実行する；ように更に構成された、
第１８の局面に記載のシステム。
［第２０の局面］
前記上向き地震波動場に対してグラウンドロールの除去を実行し；
仮想震源ショットギャザーをスタックし；
前記仮想震源ショットギャザーの出力を開始する；ように更に構成された、
第１５の局面に記載のシステム。